位场各阶垂向导数换算的新正则化方法  

A new regularization method for calculating the vertical derivatives of the potential field

作　　者：曾小牛 李夕海[1] 贾维敏[1] 刘代志[1] 

机构地区：[1]第二炮兵工程大学,西安710025

出　　处：《地球物理学报》2015年第4期1400-1410,共11页Chinese Journal of Geophysics

基　　金：国家自然科学基金(41374154)资助

摘　　要：位场垂向导数大量应用于位场数据处理与解释中.当前广泛采用的位场各阶垂向导数换算方法为基于Laplace方程并结合波数域和空间域方法的具有递推特性的ISVD(integrated second vertical derivative)算法.本文在位场垂向导数换算的正则化方法和径向平均功率谱的基础上,提出一种位场各阶垂向导数换算的新正则化方法.新正则化方法仅需通过分析位场径向平均功率谱来确定一个截止波数,即可稳定换算位场各阶垂向导数.理论模型和实测数据实验结果表明:(1)新正则化方法物理意义明确、计算简单,且各阶垂向导数换算的稳定性和精度明显优于ISVD算法;(2)在用新正则化方法求得各阶垂向导数的基础上,利用泰勒级数法可以获得大深度、高精度的位场向下延拓结果.Vertical derivatives of the potential field are commonly used in processing and interpretation of potential field data. These derivatives can narrow the width of anomalies and locate source bodies more accurately. The higher the order of the derivatives used, the more pronounced the effect. But as derivative filters are a form of high-pass filters, noise in the data is enhanced similarly. Thus calculation of higher vertical derivatives of the potential field is unstable operations in the sense of mathematical physics definition and their incorrect evaluation can contribute to enlargement of instability of the whole method. Many approaches have been proposed to solve this problem. Recently, the ISVD （integrated second vertical derivative） algorithm which is based on the Laplace equation and combines the wavenumber domain method with the space domain method has become the main method for the calculation of the vertical derivatives of the potential field. One way to stabilize higher vertical derivative evaluation is the utilization of the Tikhonov regularization. The key of this method is the choice of the regularization parameter. Based on the regularization theory and the radially averaged power spectrum of the potential field, we propose a new regularization method for calculating the vertical derivatives of the potential field. First, we simply define a special wavenumber named the cutoff wavenumber to divide the potential field spectrum into a signal part and a noise part based on the radially averaged power spectrum of potential field data. Next, we use the conventional vertical derivative operator to process the signal and the Tikhonov regularization operator to suppress the noise. Then, the parameters of the improved operator are defined by the cutoff wavenumber which has an obvious physical sense. The new regularization method only needs determining an ad hoc cutoff wavenumber based on the analysis of the radially averaged power spectrum of the potential field. Moreover, the new regularized method can n

关 键 词：位场 垂向导数 ISVD算法 TIKHONOV正则化 径向平均功率谱 向下延拓 

分 类 号：P631[天文地球—地质矿产勘探]