The existence of multiple reflections brings difficulty to seismic data processing and interpretation in seismic reflection exploration. Parabolic Radon transform is widely used in multiple attenuation because it is easily implemented, highly robust and efficient. However, finite seismic acquisition aperture of seismic data causes energy diffusion in the Radon domain, which leads to multiple residuals. In this paper, we propose a sparse parabolic Radon transform with the nonconvex Lq1-Lq2(0<q1,q2<1) mixed regularization (SPRTLq1-Lq2) that constrains the sparsity of primary and multiple reflections to overcome the energy diffusion and improve the effect of multiple attenuation, respectively. This nonconvex mixed regularization problem is solved approximately by the alternating direction method of multipliers (ADMM) algorithm, and we give the convergence conditions of the ADMM algorithm. The proposed method is compared with least squares parabolic Radon transform (LSPRT) and sparse parabolic Radon transform based on L1 regularization (SPRTL1) for multiple attenuation in the synthetic data and field data. We demonstrate that it improves the sparsity and resolution of the Radon domain data, and better results are obtained.