To address the problems of irregular trace spacing and statics correction, simultaneous regularization and wave-equation statics (WE statics) are implemented by least-squares inversion. In general, inversion is found to be intractable in three dimensions, so series approximation is made to reduce significantly the number of required integrals. The resulting operator is suitable for direct inversion or for use with gradient methods. Real and synthetic data are used to determine the viability of the inversion. For synthetic data, even for severe velocity variation and topography, inversion converges to an acceptable solution, and aliasing is reduced significantly. Similarly, for real data, inversion is found to return an antialiased, regularized result with WE statics applied.