正则化方法在比值类位场边缘识别方法中的研究

doi:10.11720/wtyht.2019.1112

摘要
图/表
参考文献
相关文章
Metrics

全文: PDF(3283 KB) HTML
输出: BibTeX | EndNote (RIS)

摘要

位场边缘识别方法有很多种,其中归一化标准差法(NSTD)、倾斜角法(Ta)和Theta Map(cosθ)等方法属于比值类方法。比值类方法在计算过程中会出现分母接近于0或者等于0这种情况,致使计算结果不稳定,并产生错误的边缘识别结果。为此,对比值类边缘识别方法计算公式中的分母加一个大于零的正则化因子,不但解决了比值类方法的数值计算稳定性问题,而且提高了部分比值类边缘识别方法识别结果的精度。通过理论模型和实际资料检验了新方法的稳定性、精度以及有效性。正则化因子的引入同样可以改善以比值类方法为基础构建的二阶导数类边缘识别方法的识别效果,如倾斜角总水平导数(Ta-THDR)的识别效果。正则化这一思想不但可以解决比值类位场边缘识别方法的数值计算问题,而且可以解决比值类方法的数值计算问题。

	服务

	把本文推荐给朋友
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	何涛
	王万银
	黄金明
	张明华
	杨敏

关键词 ：边缘识别, 比值类方法, 正则化方法, 稳定性, 精度

Abstract：

The ratio methods are one sort of edge recognition methods by using potential field, which contains the Normalized Standard Deviation method (NSTD), Tilt Angle method (Ta) and Theta Map (cosθ). If the denominator of ratio method closes or even equals to zero in the process of calculation, the result obtaining from ratio methods is unstable and may even bear little resemblance to the true geology. In order to relief this problem, a regularization factor, which is greater than zero, is added in the denominator of the ratio methods’ formula, which not only enhances the numerical stability of the ratio methods but also improves the accuracy of some ratio edge recognition methods. The stability, accuracy and effectiveness of the new method is verified by testing synthetic models and calculating real data. Also, the introduction of the regularization factor also can improve the effect of recognizing edge by the second-derivative edge recognition methods, which is based on the ratio methods, such as the Total Horizontal Derivative of the Tilt Angle (Ta-THDR). The idea of regularization can not only solve the numerical calculation problem of the ratio methods of edge recognition for potential field, but also solve the numerical calculation problem of all ratio methods.

Key words： edge recognition ratio methods regularization method stability accuracy

收稿日期: 2018-12-03 出版日期: 2019-04-15

P631

基金资助:国家重点研发计划项目2017YFC0602200之课题“航空地球物理综合处理解释方法研究及软件开发”(2017YFC0602202);中国地质调查局发展研究中心“智能地质调查系统开发与推广”项目(121201004000150014)

作者简介: 何涛(1994-),男,在读硕士研究生,研究方向是重、磁方法理论及应用。Email: hetaochd@163.com

引用本文:

何涛, 王万银, 黄金明, 张明华, 杨敏. 正则化方法在比值类位场边缘识别方法中的研究[J]. 物探与化探, 2019, 43(2): 308-319.
Tao HE, Wan-Yin WANG, Jin-Ming HUANG, Ming-Hua ZHANG, Min YANG. The research of the regularization method in the ratio methods of edge recognition by potential field. Geophysical and Geochemical Exploration, 2019, 43(2): 308-319.

链接本文:

https://www.wutanyuhuatan.com/CN/10.11720/wtyht.2019.1112 或 https://www.wutanyuhuatan.com/CN/Y2019/V43/I2/308

Table1 使用模型参数说明

Fig.1 模型A重力异常

Fig.2 A模型比值类方法边缘识别结果对比

Fig.3 A模型倾斜角总水平导数边缘识别结果对比

Fig.4 模型B重力异常

Fig.5 B模型比值类方法边缘识别结果对比

Fig.6 B模型倾斜角总水平导数边缘识别结果对比

Fig.7 模型C重力异常

Fig.8 C模型比值类方法边缘识别结果对比

Fig.9 C模型倾斜角总水平导数边缘识别结果对比

Fig.10 模型D磁力异常图与化极磁异常a—原始磁力异常;b—化极磁力异常

Fig.11 D模型比值类方法边缘识别结果对比

Fig.12 D模型倾斜角总水平导数边缘识别结果对比

Fig.13 模型E重力异常

Fig.14 E模型比值类方法边缘识别结果对比

Fig.15 E模型倾斜角总水平导数边缘识别结果对比

Fig.16 实测重力异常

Fig.17 实测重力异常比值类方法边缘识别结果对比

Fig.18 实测重力异常倾斜角总水平导数边缘识别结果对比

[1]	Miller H G, Singh V . Potential field tilt-a new concept for location of potential field sources[J]. Journal of Applied Geophysics, 1994,32:213-217. doi: 10.1016/0926-9851(94)90022-1
[2]	Hood P , McClure D J . Gradient measures in ground magnetic prospecting[J]. Geophysics, 1965,30(3):403-410. doi: 10.1190/1.1439592
[3]	Bhattacharyya B K . Two-dimensional harmonic analysis as a tool for magnetic interpretation[J]. Geosciences, 1965,30(5):829-857. doi: 10.1190/1.1439658
[4]	Cordell L . Gravimetric expression of graben faulting in Santa Fe Country and the Espanola Basin, New Mexico[C]//New Mexico Geol. Soc. Guidebook, 1979,30th Field Conf: 59-64.
[5]	Wigns C, Perez C, Kowalczyk P . Theta map: Edge detection in magnetic data[J]. Geophysics, 2005,70(4):L39-L43. doi: 10.1190/1.1988184
[6]	Nabighian M N . The analytic signal of two-dimensional magnetic bodies with polygonal cross-section: its properties and use for automated anomaly interpretation[J]. Geophysics, 1972,37(3):507-517. doi: 10.1190/1.1440276
[7]	Coorper G R J, Cowan D R . Edge enhancement of potential-field data using normalized statistics[J]. Geophysics, 2008,73(3):H1-H4. doi: 10.1190/1.2837309
[8]	王万银 . 位场边缘识别技术研究[D]. 西安: 长安大学, 2009.
[8]	Wang W Y . The research on the edge recognition methods and techniques for potential field[D]. Xi’an:Chang’an University, 2009.
[9]	Coorper G R J, Cowan D R . Enhancing potential field data using filters based on the local phase[J]. Computers & Geosciences, 2006,32(10):1585-1591. doi: 10.1016/j.cageo.2006.02.016
[10]	马国庆, 黄大年, 于平 , 等. 改进的均衡滤波器在位场数据边界识别中的应用[J]. 地球物理学报, 2012,55(12):4288-4295. doi: 10.6038/j.issn.0001-5733.2012.12.040
[10]	Ma G Q, Huang D N, Yu P , et al. Application of improved balancing filters to edge identification of potential field data[J]. Chinese Journal of Geophysics, 2012,55(12):4288-4295.
[11]	王明, 郭志宏, 骆遥 . 扩展的倾斜角与倾斜角总水平导数[J]. 地球物理学进展, 2013,28(3), 1453-1463. doi: 10.6038/pg20130340
[11]	Wang M, Guo Z H, Luo Y , et al. Extended tilt angle and total horizontal derivatives[J]. Progress in Geophysic, 2013,28(3):1453-1463.
[12]	Coorper G R J . Reply to a discussion about the “Hyperbolic tilt angle method” by Zhou et al[J] Computers & Geosciences, 2013,52:496-497. doi: 10.1016/j.cageo.2012.11.007
[13]	马国庆, 杜晓娟, 李丽丽 . 位场边缘识别的新方法—增强型水平导数法[J]. 地球物理学进展, 2013,28(1):402-408. doi: 10.6038/pg20130145
[13]	Ma G Q, Du X J, Li L L . New edge detection method of potential field data-enhanced horizontal derivative method[J]. Progress in Geophysics, 2013,28(1):402-408.
[14]	Li L L, Huang D N, Han L G , et al. Optimised edge detection filters in the interpretation of potential field data[J]. Exploration Geophysics, 2014,45(3):171-176. doi: 10.1071/EG13059
[15]	顾廷杰, 吴燕冈, 袁园 , 等. 应用加强解析信号倾斜角进行位场数据的边界检测[J]. 地球物理学报, 2016,59(7), 2694-2702.
[15]	Gu T J, Wu Y G, Yuan Y , et al. Edge detection of potential field data using an enhanced analytic signal tilt angle[J]. Chinese Journal of Geophysics, 2016,59(7):2694-2702.
[16]	Yuan Y, Huang D N, Yu Q L . Edge detection of potential field data with improved structure tensor methods[J]. Journal of Applied Geophysical, 2014,172(2):461-472. doi: 10.1016/j.jappgeo.2014.06.013
[17]	袁园, 黄大年, 余青露 . 利用加强水平方向总水平导数对位场全张量数据进行边界识别[J]. 地球物理学报, 2015,58(7):2556-2565. doi: 10.6038/cjg20150730
[17]	Yuan Y, Huang D N, Yu Q L . Using enhanced directional total horizontal derivatives to detect the edges of potential field full tensor data: Chinese Journal of Geophysics, 2015,58(7):2556-2565.
[18]	Tikhonov A N, Glasko V B . Using of the regularization method in non-linear problems[J]. Zh. Vychislit. Matem. i matem. fiz., 1965,5(N):463-473. doi: 10.1016/0041-5553(65)90150-3
[19]	Pašteka R, Richter F P, Karcol R , et al. Regularized derivatives of potential and their role in semi-automated interpretation methods[J]. Geophysical Prospecting, 2009,57:507-516. doi: 10.1111/j.1365-2478.2008.00780.x
[20]	曾小牛, 李夕海, 贾维敏 , 等. 位场各阶垂向导数换算的新正则化方法[J]. 地球物理学报, 2015,58(4):1400-1410. doi: 10.6038/cjg20150426
[20]	Zeng X N, Li X H, Su J . An adaptive iteration method for downward continuation of potential-field data from a horizontal plane: Geophysics, 2015,58(4):1400-1410.
[21]	栾文贵 . 场位解析延拓的稳定化算法[J]. 地球物理学报, 1983,26(3):263-273. doi:
[21]	Luan W G . Stabilization algorithm for analytical continuation of potential field[J]. Chinese Journal of Geophysics, 1983,26(3):263-273.
[22]	梁锦文 . 位场向下延拓的正则化方法[J]. 地球物理学报, 1989,33(5):600-608.
[22]	Liang J W . Regularization method for downward continuation of potential field[J]. Chinese Journal of Geophysics, 1989,33(5):600-608.
[23]	Pašteka R, Karcol R, Kuŝnirák D , et al. Regcont: A Matlab based program for stable downward continuation of geophysical potential fields using Tikhonov regularization[J]. Computers & Geosciences, 2012(49):278-289.
[24]	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[J]. Computers & Geosciences, 2013, ( 52):269-280.
[25]	Li Y G, Sarah G R D, Richard A K , et al. Enhancemant of magnetic data by stable downward continuation for UXO application[J]. IEEE transactions on Geoscience and Remote Sensing, 2013,51(6):3605-3614. doi: 10.1109/TGRS.2012.2220146
[26]	赵亚博, 刘天佑 . 迭代Tikhonov正则化位场向下延拓方法及其在尕林铁矿的应用[J]. 物探与化探, 2015,39(4):743-748. doi: 10.11720/wtyht.2015.4.14
[26]	Zhao Y B, Liu T Y . The iterative regularization method for downward continuation of potential fields and its application to the Galinge iron deposit[J]. Geophysical and Geochemical Exploration, 2015,39(4):743-748.
[27]	Verduzco B, Fairhead J D, Green C M, Mackenzie C . The meter reader-New insights into magnetic derivatives for structural mapping[J]. The Leading Edge, 2004,23(2):116-119. doi: 10.1190/1.1651454
[28]	王万银, 邱之云, 杨永 , 等. 位场边缘识别方法研究进展[J]. 地球物理学进展, 2010,25(1):196-210. doi: 10.3969/j.issn.1004-2903.2010.01.027
[28]	Wang W Y, Qiu Z Y, Yang Y , et al. Some advances in the edge recognition of the potential field[J]. Progress in Geophysics, 2010,25(1):196-210.
[29]	Ma G Q, Huang D N, Cai Liu . Step-edge detection filters for the Interpretation of potential field data[J]. Pure and Applied Geophysics, 2016,173:795-803. doi: 10.1007/s00024-015-1053-6
[30]	Ma G. Q . Edge detection of potential field data using improved local phase filter[J]. Exploration Geophysics, 2013,44:36-41. doi: 10.1071/EG12022
[31]	周帅, 黄大年, 焦健 . 基于三维构造张量的位场边界识别滤波器[J]. 地球物理学报, 2016,59(10):3847-3858. doi: 10.6038/cjg20161028
[31]	Zhou S, Huang D N, Jiao J . Construction of potential Field Boundary recognition filter based on three-dimensional Construction of Zhang Liang[J]. Journal of Geophysics, 2016,59(10):3847-3858.
[32]	颜廷杰, 吴燕冈, 袁园 , 等. 应用加强解析信号倾斜角进行位场数据的边界检测[J]. 地球物理学报, 2016,59(7):2694-2702.
[32]	Yan T J, Wu Y G, Yuan Y , et al. Boundary Detection of potential Field data by strengthening the inclination Angle of Analytic signal[J]. Journal of Geophysics, 2016,59(7):2694-2702.
[33]	陈国强, 马国庆 . 位场数据解释的Theta-Depth法[J]. 地球物理学报, 2016,59(6):2225-2231. doi: 10.6038/cjg20160625
[33]	Chen G Q, Ma G Q . Theta-Depth method for potential Field data interpretation[J]. Journal of Geophysics, 2016,59(6):2225-2231.

[1]	冯旭亮, 魏泽坤. 基于界面反演增强的位场边缘识别方法[J]. 物探与化探, 2022, 46(1): 130-140.
[2]	王猛, 刘媛媛, 王大勇, 董根旺, 田亮, 黄金辉, 林曼曼. 无人机航磁测量在荒漠戈壁地区的应用效果分析[J]. 物探与化探, 2022, 46(1): 206-213.
[3]	奚小环, 侯青叶, 杨忠芳, 叶家瑜, 余涛, 夏学齐, 成杭新, 周国华, 姚岚. 基于大数据的中国土壤背景值与基准值及其变化特征研究——写在《中国土壤地球化学参数》出版之际[J]. 物探与化探, 2021, 45(5): 1095-1108.
[4]	崔志强, 胥值礼, 李飞. 塔西南高精度航磁油气地质构造调查[J]. 物探与化探, 2021, 45(4): 846-858.
[5]	游军, 张小明, 罗乾周, 史朝洋, 杨运军, 陈剑祥, 袁攀. 略阳白雀寺杂岩体以西隐伏深大断裂的推断依据及其意义[J]. 物探与化探, 2021, 45(3): 645-652.
[6]	单希鹏, 谢汝宽, 梁盛军, 余学中. 直升机TEM测量影响因素分析[J]. 物探与化探, 2021, 45(1): 178-185.
[7]	李忠平, 戴广凯, 张茂辉. DEM网格间距对重力远区地改精度的影响及效果[J]. 物探与化探, 2020, 44(6): 1399-1407.
[8]	屈利军, 王庆, 李波, 姚伟. 综合物探方法在湖南香花岭矿田三合圩矿区深部成矿规律研究中的应用[J]. 物探与化探, 2020, 44(6): 1313-1321.
[9]	何继善. 广域电磁法理论及应用研究的新进展[J]. 物探与化探, 2020, 44(5): 985-990.
[10]	张菲菲, 孙建伟, 韩波, 杜润林, 王万银. SAG-2M型与KSS31M型海洋重力仪比测结果分析[J]. 物探与化探, 2020, 44(4): 870-877.
[11]	陈晓晶, 虎新军, 李宁生, 安百州, 白亚东. 银川平原基于地球物理资料三维建模的深部地质构造研究[J]. 物探与化探, 2020, 44(2): 245-253.
[12]	赵敏, 盛勇, 戚良刚. 高精度重磁测量在覆盖区找矿中的应用——以无为县蔚山铁铜矿预查为例[J]. 物探与化探, 2019, 43(6): 1211-1216.
[13]	岳航羽, 张保卫, 王凯, 李金丽, 王广科, 陈德元. 一种适用于沿海滩涂区的浅层高精度地震探测技术[J]. 物探与化探, 2019, 43(6): 1225-1235.
[14]	毕炳坤, 常云真, 施强, 申随水. 综合物探在崤山东部浅覆盖区勘查银多金属矿床中的应用[J]. 物探与化探, 2019, 43(5): 976-985.
[15]	何祎, 雷晓东, 关伟, 李晨, 李巧灵. 北京副中心地区燕郊断裂空间展布特征[J]. 物探与化探, 2019, 43(3): 461-467.

Viewed

Full text

Abstract

Cited

Shared

Discussed