Research and Application of Three‐Dimensional Density Interface Inversion

Based on the Fourier transform, the Parker-Oldenburg algorithm in the frequency domain was extended to the three-dimensional case in which the density changes with depth. From this, a gravity interface inversion formula was derived in which the assumed density varies laterally and vertically. Iterative convergence is assured by fixing a particular depth as the datum plane below the surface to reduce the interface fluctuation. The results of an example of synthetic gravity data indicate that the proposed method has high precision and rapid convergence, with high practical value for the inversion of density interfaces. This method was also used to determine the Moho depth beneath northern China. The results were verified by seismic sounding data.