位场向下延拓的波数域广义逆算法

Wavenumber domain generalized inverse algorithm for potential field downward continuation

位场向下延拓是位场数据处理和反演中的重要运算,但是它的不稳定性影响了它在许多处理和反演方法技术中的应用.本文通过把位场向下延拓视为向上延拓的反问题,得到向下延拓的褶积型线性积分方程,再利用Fourier变换矩阵的正交对称特性,并结合矩阵的奇异值分解和广义逆原理,提出了一种稳定的不需要进行求逆运算的位场向下延拓广义逆方法——波数域广义逆算法,解决了位场大深度向下延拓的不稳定性问题.把这种方法用于三维理论模型数据和实际磁场数据的向下延拓获得了理想的结果.

Downward continuation is an important operation in potential field data processing and inversion, but its instability problem influences its application in many processing and inversion techniques. In this paper, taking the downward continuation of potential field as an inverse problem of upward continuation, we obtain a convolution type linear integral equation for downward continuation. Making use of the orthogonal symmetry characteristic of Fourier transform matrix, and combining the principles of singular value decomposition of matrix and generalized inverse, we proposed a stable generalized inverse method for downward continuation of potential field, called wavenumber domain generalized inverse algorithm, which doesn＇t need the computation for inverse matrix. It resolves the instability of potential field downward continuation of large depth. The applications of the method to the downward continuation of 3-D theoretical model data and real magnetic field data give ideal results.