We study seismic data decomposition problem in this paper. Seismic data representation using Gaussian beams are proposed. This is an 10 quasi-norm constrained minimization problem. In solving the 10 quasi-norm minimization, scanning event in phase is proposed and a gradient descent method of rapid monotonic decrease is utilized. The sparse optimization method proposed in this paper possesses the ability to remove noise and is validated with numerical simulations.