An improved stable algorithm of the reduction-to-the-pole at low latitudes
Zhang Qi1,2, Zhang Yingtang1, Li Zhining1, Fan Hongbo1
1. Department of Vehicle and Electrical Enginee-ring, Shijiazhuang Campus, Army Engineering University, Shijiazhuang, Hebei 050003, China; 2. Unit 94019 PLA, Hetian, Xinjiang 848000, China
Abstract:As the routine reduction-to-the pole (RTP) operator at low latitudes is usually unstable in the frequency domain,we propose an improved stable RTP algorithm at low latitudes.First the high frequency noise in amplification zones is suppressed with the local-suppression and global-suppression RTP approaches based on routine RTP operator (RO) modifications.Then the introduced parameters are optimized by calculating the correlation coefficient between the RTP results and the norma-lized source strength.RTP results with the optimal parameters are obtained,and the accuracy of these RTP results is accordingly improved.Finally,to resolve the problem of high-frequency amplification effect in the optimal RTP operators,the Tikhonov regularization method is utilized to stabilize the operators.Experimental results on theoretical model indicate that the proposed algorithm eliminates the high-frequency noise in the original data and improves the signal-to-noise ratio of the RTP results,and the accuracy and stability of the low latitude RTP are greatly improved.
Baranov V.A new method for interpretation of aeromagnetic maps:pseudo-gravimetric anomalies.Geophysics,1957,22(2):359-382.
[2]
Silva J B C.Reduction to the pole as an inverse pro-blem and its application to low-latitude anomalies.Geo-physics,1986,51(2):369-382.
[3]
Leao J W D,Silva J B C.Discrete linear transformations of potential field data.Geophysics,1989,54(4):497-507.
[4]
Bhattacharyya B K.Two-dimensional harmonic analysis as a tool for magnetic interpretation.Geophysics,1965,30(5):829-857.
[5]
柴玉璞.从化极算法误差方程看各种波数域低纬度化极方法.石油地球物理勘探,2012,47(3):496-505.Chai Yupu.Discussion of wavenumber domain RTP methods at low latitudes.OGP,2012,47(3):496-505.
[6]
姚长利,管志宁,高德章等.低纬度磁异常化极方法——压制因子法.地球物理学报,2003,46(5):690-696.Yao Changli,Guan Zhining,Gao Dezhang et al.Reduction to the pole of magnetic anomalies at low latitude with suppression filter.Chinese Journal of Geophysics,2003,46(5):690-696.
[7]
Li X.Magnetic reduction-to-the-pole at low latitudes:Observations and considerations.The Leading Edge,2008,27(8):990-1002.
[8]
姚长利,黄卫宁,张聿文等.直接阻尼法低纬度磁异常化极技术.石油地球物理勘探,2004,39(5):600-606.Yao Changli,Huang Weining,Zhang Yuwen et al.Reduction to the pole of magnetic anomalies at low latitude with direct damping.OGP,2004,39(5):600-606.
[9]
林晓星,王平.一种改进的低纬度磁场化极方法——变频双向阻尼因子法.地球物理学报,2012,55(10):3477-3484.Lin Xiaoxing,Wang Ping.An improved method for reduction to the pole of magnetic field at low latitude -the method of frequency conversion bidirectional damping factor.Chinese Journal of Geophysics,2012,55(10):3477-3484.
[10]
Guo L,Shi L,Meng X.The antisymmetric factor me-thod for magnetic reduction to the pole at low latitudes.Journal of Applied Geophysics,2013,92:103-109.
[11]
Zhang H,Marangoni Y R,Hu X et al.NTRTP:A new reduction to the pole method at low latitudes via a nonlinear thresholding.Journal of Applied Geophy-sics,2014,111:220-227.
[12]
Mendona C A,Silva J B C.A stable truncated series approximation of the reduction-to-the-pole operator.Geophysics,1993,58(8):1084-1090.
[13]
张培琴,赵群友.低磁纬度区航磁异常变倾角磁方向转换方法.物探化探计算技术,1996,18(3):206-214.Zhang Peiqin,Zhao Qunyou.Methods of the magnetic direction transform of aeromagnetic anomalies with differential inclinations in low magnetic latitudes.Computing Techniques for Geophysical and Geoche-mical Exploration,1996,18(3):206-214.
[14]
张振宇,王刚,胡祥云等.张量CSAMT方法及对比实验.石油地球物理勘探,2017,52(4):869-874.Zhang Zhenyu,Wang Gang,Hu Xiangyun et al.Tensor CSAMT technical and experiments.OGP,2017,52(4):869-874.
[15]
Wilson H.Analysis of the magnetic gradient tensor.Defence Research Establishment Pacific:Canada Technical Memorandum,1985,85:13-47.
[16]
Zeng X,Liu D,Li X et al.An improved regularized downward continuation of potential field data.Journal of Applied Geophysics,2014,106:114-118.
[17]
Yin G,Zhang Y,Mi S et al.Calculation of the magne-tic gradient tensor from total magnetic anomaly field based on regularized method in frequency domain.Journal of Applied Geophysics,2016,134:44-54.
[18]
曹静杰.基于广义高斯分布和非Lem>p范数正则化的地震稀疏盲反褶积.石油地球物理勘探,2016,51(3):428-433.Cao Jingjie.Seismic sparse blind deconvolution based on generalized Gaussian distribution and non-convex Lem>p normal regularization.OGP,2016,51(3):428-433.
[19]
杨国权,丁鹏程,李振春等.应用线性程函方程和整形正则化的三维初至波旅行时层析.石油地球物理勘探,2017,52(2):264-272.Yang Guoquan,Ding Pengcheng,Li Zhenchun et al.3D first-arrival traveltime tomography with linearized eikonal and shaping regularization.OGP,2017,52(2):264-272.
[20]
Hansen P C,O'Leary D P.The use of the L-curve in the regularization of discrete ill-posed problems.SIAM Journal on Scientific Computing,1993,14(6):1487-1503.