Derivative calculation is a basic method used in data processing and interpretation of potential fields,but the instability of the result restricts the application of higher-order derivative in practical data.Thus,this paper proposes an iterative method for higher derivative calculation of potential fields to resolve the problem of instability.The authors added a lowpass filter which contained the derivative order number and derivative operator to the higher derivative operator,and used the iterative method to perform successive approximation.On such a basis,the authors deduced the recursive formula of the iterative method of higher-order derivative.In addition,the authors proved the convergence of the iterative method and made a comparative analysis on filtering characteristics of derivative operator by using different methods. In model test and application,the comparative analysis of the results of the iterative method and conventional method for higher-order derivative indicates that the iterative method has advantages of stronger stability and higher calculation accuracy.
王彦国, 张瑾. 位场高阶导数的波数域迭代法[J]. 物探与化探, 2016, 40(1): 143-147.
WANG Yan-Guo, ZHANG Jin. The iterative method for higher-order derivative calculation of potential fields in wave number domain. Geophysical and Geochemical Exploration, 2016, 40(1): 143-147.
[1] Blakely R J.Potential theory in gravity and magnetic applications[M].Cambridge University Press,1995:324-327.[2] 骆燕,李晓禄,蔡文良,等.潮水地区航磁梯度初步分析[J].东华理工学院学报,2007,30(2):164-170.[3] 侯重初,施志群.实现重磁异常规格化总梯度的傅氏积分法[J].石油物探,1986,25(3):75-86.[4] 张凤旭,孟令顺,张凤琴,等.利用Hilbert变换计算重力归一化总梯度[J].地球物理学报,2005,48(3):704-709.[5] Thompson D T.EUlDPH-A new technique for making computer assisted depth estimates from magnetic data[J].Geophysics,1982,47:31-37.[6] Stavrev P,Reid A.Euler deconvolution of gravity anomalies from thick contact/fault structures with extended negative structural index[J].Geophysics,2010,75(6):151-158.[7] 黄临平,管志宁.利用磁异常总梯度模确定磁源边界位置[J].华东地质学院学报,1998,21(2):143-150.[8] 张季生,高锐,李秋生,等.欧拉反褶积与解析信号相结合的位场反演方法[J].地球物理学报,2011,54(6):1634-1641.[9] 骆遥,王明,罗峰,等.重磁场二维希尔伯特变换——直接解析信号解释方法[J].地球物理学报,2011,54(7):1912-1920.[10] 刘保华,张维冈,孟恩.重力异常垂向一阶导数的一种简便算法[J].青岛海洋大学学报,1995,25(2):233-238.[11] 汪炳柱.用样条函数法求重力异常二阶垂向导数和向上延拓计算[J].石油地球物理勘探,1996,31(3):415-422.[12] Lourenco J S,Morrison H F.Vector magnetic anomalies derived from measurements of single component of the field[J].Geophysics,1973,38:395-368.[13] Kevin L M,Juan H H.The complete gravity gradient tensor derived from the vertical component of gravity:A Fourier transform technique[J].Journal of Applied Geophysics,2001,46:159-l74.[14] 姚长利,管志宁,黄卫宁.位场转换的抽样分组法[J].石油地球物理勘探,1997,32(5):696-702.[15] 王万银,邱之云,杨永,等.位场边缘识别方法研究进展[J].地球物理学进展,2010,25(1):196-210.[16] Strakhov A V,Devitsyn V N.The reduction of observed values of a potential field to values at a constant level[J].Akad. Nauk USSR Izv. Fizika Zemli (Physics of the Solid Earth),1965,4:256-261.[17] 徐世浙.位场延拓的积分——迭代法[J].地球物理学报,2006,49(4):1176-1182.[18] 王彦国,张凤旭,王祝文,等.位场向下延拓的泰勒级数迭代法[J].石油地球物理勘探,2011,46(4):657-662.