In recent years, there has been a boom in research on lump solutions. In 1977, Manakov
et al [
16] first found the lump solution and interaction for the KP equation. In 1990, Glison
et al [
17] described the single lump solution and the N-lump solution, and confirmed that single lump solutions are only nonsingular for spectral parameters lying in certain regions of the complex plane. In 1996, Minzoni
et al [
18] used the group velocity argument to determine the propagation direction of the liner dispersive radiation generated as the lump evolves in the KdV equation, and further studied the evolution of the initial conditions of the lump-like. In 2000, Sipcic
et al [
19] studied the modified Zakharov Kuznetsov equation and confirmed that when two lumps interact, the initial energy exchange between them is followed by the emergence of a single collapsing lump and a radiation field behind it. In 2004, based on exact and numerical methods, Lu
et al [
20] analysed the interaction of two lump solitons described by the Kadomtsev–Petviashvili I equation. In 2009, Villarroel
et al [
21] derived a class of localized solutions of a (2+1)-dimensional nonlinear Schrödinger equation and studied their dynamical properties. Ma
et al [
22–
27] obtained a class of lump solutions of some nonlinear partial differential equations by the Hirota bilinear method. Wang
et al [
28] derived the lump solution when the period of complexiton solution went to infinite and investigated the dynamics of the lump solution of the Hirota bilinear equation in 2017. In 2018, Foroutan
et al [
29] studied the (3+1)-dimensional potential Yu–Toda–Sasa–Fukuyama equation by implementing the Hirota bilinear method, and acquired a type of the lump solution and five types of interaction solutions. In 2021, ${\rm{L}}\ddot{u}$
et al [
30] studied the one-lump-multi-stripe and one-lump-multi-soliton types interaction solutions to nonlinear partial differential equations.