To find exact solutions of nonlinear systems, there are many effective methods, for instance, the Lie group method,
[14-15] the inverse scattering transformation (IST),
[16] the Darboux transformation (DT),
[17] the bilinear method,
[18] and so on. To research lump solutions of NLEEs, there are a large number of works by using a quite simple and effective Hirota's bilinear method.
[19-22] Usually, lump solutions are defined and found for high dimensional equations, say, the $(2+1)$-dimensional Boussinesq equation,
[23] the $(2+1)$-dimensional Burgers equation
[24] and the $(3+1)$-dimensional Jimbo-Miwa equation.
[25]