Abstract:Downward continuation of potential field data plays an important role in gravity and magnetic data interpretation. Tikhonov regularization method is an effective method for the downward continuation of potential field. The key of this method is regularization parameter selection, which directly affect the precision and the computation time of the downward continuation. In order to effectively choose a regularization parameter, this paper proposes a new selection method based on the radial-average power spectrum of the potential field. The new method builds the relationship between the two cutoff wavenumbers which are respectively decided by the radial-average power spectrum and the Tikhonov regularization low-pass filter. Tests on theoretical gravity model and real aeromagnetic data shows that the Tikhonov downward continuation precision based on the regularization parameter selection with the proposed method is more accurate than the L-curve and C-norm methods. Moreover, the proposed method is easy to be realized.
曾小牛, 李夕海, 陈鼎新, 杨晓云. 基于径向谱的位场向下延拓正则参数选取方法[J]. 石油地球物理勘探, 2015, 50(4): 749-754.
Zeng Xiaoniu, Li Xihai, Chen Ding-xin, Yang Xiaoyun. New regularization parameter selection for potential field downward continuation based on the radial spectrum. OGP, 2015, 50(4): 749-754.
Liang Jinwen. The regularization method for downward continuation of potential field. Chinese Journal of Geophysics, 1989, 32(5): 600-608.
[2]
陈生昌, 肖鹏飞. 位场向下延拓的波数域广义逆算法. 地球物理学报, 2007, 50(6): 1816-1822] Chen Shengchang, Xiao Pengfei. Wavenumber domain generalize inverse algorithm for potential field downward continuation. Chinese Journal of Geophy-sics, 2007, 50(6): 1816-1822.
[3]
Pateka R, Karcol R, Kunirák D et al. Regcont: a Matlab based program for stable downward continuation of geophysical potential fields using Tikhonov regularization. computers & Geosciences, 2012, 49: 278-289.
[4]
Abedi M, Gholami A, Norouzi G H. A stable downward continuation of airborne magnetic data: A case study for mineral prospectivity mapping in Central Iran. computers & Geosciences, 2013, 52: 269-280.
[5]
Li Y G, Devriese S G R, Krahenbuhl R A et al. Enhancement of magnetic data by stable downward continuation for UXO application. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(6): 3605-3614.
[6]
Zeng Xiaoniu, Li Xihai, Su Juan et al. An adaptive iterative method for downward continuation of potential-field data from a horizontal plane. Geophysics, 2013, 78(4): J43-J52.
[7]
Zhou Jun, Shi Guiguo, Ge Zhilei. Study of iterative regularization methods for potential field downward continuation in geophysical navigation. Journal of Astronautics, 2011, 32(4): 787-794.
[8]
Zhang Henglei,Ravat D,Hu Xiangyun. An improved and stable downward continuation of potential field data: The truncated Taylor series iterative downward continuation method. Geophysics, 2013, 78(5): J75-J86.
[9]
Fedi M, Florio G. A stable downward continuation by using the ISVD method. Geophysical Journal International, 2002, 151(1): 146-156.
[10]
王彦国, 张凤旭, 王祝文等. 位场向下延拓的泰勒级数迭代法. 石油地球物理勘探, 2011, 46(4): 657-662] Wang Yanguo, Zhang Fengxu, Wang Zhuwen et al. Taylor series iteration for downward continuation of potential field. OGP, 2011, 46(4): 657-662.
[11]
徐世浙. 位场延拓的积分—迭代法. 地球物理学报, 2006, 49(4): 1176-1182] Xu Shizhe. The integral-iteration method for continuation of potential field. Chinese Journal of Geophysics, 2006,49(4): 1176-1182.
Liu Dongjia, Hong Tianqiu, Jia Zhihai et al. Wave number domain iteration method for downward of potential fields and its convergence. Chinese Journal of Geophysics, 2009, 52(6): 1599-1605.
Zeng Xiaoniu, Li Xihai, Niu Chao et al. The regularization-integral iteration method in wave number domain for downward continuation of potential fields. OGP, 2013, 48(4): 643-650.
Wang Bingzhu, Wang Shuoru. Spline function methods for upward continuation and downward continuation of 2D potential field. computing Techniques for Geophysical and Geochemical Exploration, 1998,20(2): 125-129.
[15]
Wang Bingzhu.2D and 3D potential-field upward continuation using splines. Geophysical Prospecting, 2006, 54(2): 199-209.
Yao Changli, Li Hongwei, Zheng Yuanman et al. Research on iteration method using in potential field transformations. Chinese Journal of Geophysics, 2012, 55(6): 2062-2078.
Ma Tao, Chen Longwei, Wu Meiping et al. The selection of regularization parameter in downward continuation of potential field based on L-curve method. Progress in Geophysics, 2013, 28(5): 2485-2494.
[18]
Pateka R, Richter F P, Karcol R et al. Regularized derivatives of potential fields and their role in semi-automated interpretation methods. Geophysical Prospecting, 2006, 57(4): 507-516.
[19]
Golub G H, Heath M, Wahba G. Generalized cross-validation as a method for choosing a good ridge parameter. Technometrics, 1979, 21: 215-223.
[20]
王彦飞. 反演问题的计算方法及其应用. 北京: 高等教育出版社, 2007.
[21]
Morozov V A. Methods for solving incorrectly posed problems. New York: Springer-Verlag, 1984.
[22]
Hansen P C, O'Leary D P. The use of the L-curve in the regularization of discrete ill-posed problems. SIAM J Sci comput, 1993, 14(6): 1487-1503.
[23]
Reginska T. A regularization parameter in discrete ill-posed problems. SIAM J Sci comput, 1996, 17(3): 740-749.
[24]
Spector A, Grant F S. Statistical models for interpreting aeromagnetic data. Geophysics, 1970, 35(2): 293-302.
[25]
Pawlowski R S. Preferential continuation for potential-field anomaly enhancement. Geophysics, 1995,60(2): 390-398.
[26]
Stefan M, Vijay D. Potential field power spectrum inversion for scaling geology. Journal of Geophysical Research, 1995, 100(B7): 12605-12616.
[27]
Ravat D, Pignatelli A, Nicolos I et al. A study of spectral methods of estimating the depth to the bottom of magnetic sources from near-surface magnetic anomaly data. Geophysical Journal International, 2007, 169(2): 421-434.
[28]
Guo Lianghui, Meng Xiaohong, Chen Zhaoxi et al. Preferential filtering for gravity anomaly separation. computers & Geosciences, 2013, 51: 247-254.