Please wait a minute...
E-mail Alert Rss
 
物探与化探  2022, Vol. 46 Issue (1): 123-129    DOI: 10.11720/wtyht.2022.1403
  方法研究·信息处理·仪器研制 本期目录 | 过刊浏览 | 高级检索 |
重力异常曲化平的改进插值—迭代法
杨婧(), 郭良辉()
中国地质大学(北京) 地球物理与信息技术学院,北京 100083
A modified interpolation-iteration method for gravity anomaly continuation from undulating surface to plane
YANG Jing(), GUO Liang-Hui()
School of Geophysics and Information Technology, China University of Geosciences,Beijing 100083,China
全文: PDF(1920 KB)   HTML
输出: BibTeX | EndNote (RIS)      
摘要 

重力异常曲化平即将起伏观测面上的重力异常延拓到平面上,从而为频率域处理和反演提供平坦观测面的重力异常数据。本文在常规插值—迭代法基础上,给出重力异常曲化平的改进插值—迭代法,即在异常迭代修正过程中引入起伏观测面修正因子,加快曲化平迭代收敛,促进曲化平效果提升。理论模型试验表明本文方法适用于观测面起伏较大、延拓跨度大的复杂条件下的稳定、有效曲化平,效果优于常规插值—迭代法。川滇地区实际重力异常数据试验表明本文曲化平方法有效增强了异常信号和细节特征,为后续处理和解释提供可靠数据。

服务
把本文推荐给朋友
加入引用管理器
E-mail Alert
RSS
作者相关文章
杨婧
郭良辉
关键词 曲化平重力异常插值—迭代法川滇地区    
Abstract

Gravity anomaly continuation from undulating surface to plane can provide gravity data on a flat horizontal plane for frequency-domain data processing and inversion. Based on the theory of conventional interpolation-iteration methods, this study proposed a modified interpolation-iteration method by introducing a correction factor of the undulating observation surface in the iteration and correction process. The improved method accelerated the iterative convergence speed and promoted the continuation effects. The theoretical model-based tests show that this method can be used to achieve stable and effective large-span gravity anomaly continuation from greatly undulating surface to plane. The continuation results presented by this method are better than those obtained using conventional interpolation-iteration methods. The application of Bouguer gravity anomaly data of the Sichuan-Yunnan region demonstrates that the modified interpolation-iteration method effectively enhanced anomalous signals and details and can provide reliable data for subsequent processing and interpretation.

Key wordscontinuation from undulating surface to plane    gravity anomaly    interpolation-iteration method    Sichuan-Yunnan region
收稿日期: 2021-07-22      修回日期: 2021-10-18      出版日期: 2022-02-20
ZTFLH:  P631  
基金资助:国家重点研发计划项目“月球内部圈层结构与演化过程的研究”(2021YFA0715101);国家自然科学基金面上项目“全波形与重力联合反演方法及在南岭—武夷—东南沿海成矿带岩石圈三维精细结构成像的应用”(41974101);中央高校基本科研业务费专项资金;地质过程与矿产资源国家重点实验室科技部专项经费
通讯作者: 郭良辉
作者简介: 杨婧(1999-),女,硕士研究生,主要从事重震数据处理和反演研究工作。Email: 2010200015@cugb.edu.cn
引用本文:   
杨婧, 郭良辉. 重力异常曲化平的改进插值—迭代法[J]. 物探与化探, 2022, 46(1): 123-129.
YANG Jing, GUO Liang-Hui. A modified interpolation-iteration method for gravity anomaly continuation from undulating surface to plane. Geophysical and Geochemical Exploration, 2022, 46(1): 123-129.
链接本文:  
https://www.wutanyuhuatan.com/CN/10.11720/wtyht.2022.1403      或      https://www.wutanyuhuatan.com/CN/Y2022/V46/I1/123
Fig.1  插值—迭代法曲化平示意
Fig.2  本文改进的插值—迭代法算法流程
Fig.3  理论模型起伏观测面与重力异常
a—海拔0 m平面的理论重力异常;b—起伏观测面高程;c—起伏观测面的理论重力异常
Fig.4  不同曲化平方法的结果及与海拔0 m平面理论重力异常对比
a、b—常规插值-迭代法曲化平结果及与海拔0 m平面理论重力异常残差;c、d—本文改进的插值-迭代法曲化平结果及与海拔0 m平面理论重力异常残差
Fig.5  剖面A不同修正因子S曲化平的对比
Fig.6  川滇地区曲化平前后重力异常对比
a—川滇地区地形高程;b—去噪后的布格重力异常;c—本文改进插值—迭代法曲化平后布格重力异常;d—曲化平前后的异常残差
Fig.7  剖面B的对比
Fig.8  川滇地区曲化平前后垂向导数对比
a—布格重力异常的垂向导数;b—曲化平后布格重力异常的垂向导数
[1] 曾华霖. 重力场与重力勘探[M]. 北京: 地质出版社, 2005.
[1] Zeng H L. Gravity field and gravity exploration[M]. Beijing: Geological Publishing House, 2005.
[2] Blakely R J. Potential theory in gravity and magnetic applications [M]. Cambridge: Cambridge University Press, 1995:316-320.
[3] Guo L, Meng X, Chen Z, et al. Preferential filtering for gravity anomaly separation[J]. Computers & Geosciences, 2013,51:247-254.
doi: 10.1016/j.cageo.2012.09.012
[4] Oldenburg D W. The inversion and interpretation of gravity anomalies[J]. Geophysics, 1974,39(4):526-536.
doi: 10.1190/1.1440444
[5] Cui Y, Guo L. A wavenumber-domain iterative approach for rapid 3-D imaging of gravity anomalies and gradients[J]. IEEE Access, 2019,7(1):34179-34188.
doi: 10.1109/Access.6287639
[6] Pilkington M, Boulanger O. Potential field continuation between arbitrary surfaces-Comparing methods[J]. Geophysics, 2017,82(3):J9-J25.
doi: 10.1190/geo2016-0210.1
[7] Dampney C N G. The equivalent source technique[J]. Geophysics, 1969,34(1):39-53.
doi: 10.1190/1.1439996
[8] David A E. Equivalent sources used as an analytic base for processing total magnetic field profiles[J]. Geophysics, 1973,38(2):339-348.
doi: 10.1190/1.1440344
[9] Bhattacharyya B K, Chan K C. Reduction of magnetic and gravity data on an arbitrary surface acquired in a region of high topographic relief[J]. Geophysics, 1977,42(7):1411-1430.
doi: 10.1190/1.1440802
[10] Hansen R O, Miyazaki Y. Continuation of potential fields between arbitrary surfaces[J]. Geophysics, 1984,49(6):787-795.
doi: 10.1190/1.1441707
[11] Pilkington M, Urquhart W E S. Reduction of potential field data to a horizontal plane[J]. Geophysics, 1990,55(5):549-555.
doi: 10.1190/1.1442866
[12] Xia J, Sprowl D R, Adkins-Heljeson D. Correction of topographic distortions in potential-field data: A fast and accurate approach[J]. Geophysics, 1993,58(4):515-523.
doi: 10.1190/1.1443434
[13] Henderson R G, Cordell L. Reduction of unevenly spaced potential field data to a horizontal plane by means of finite harmonic series[J]. Geophysics, 1971,36(5):1046-1046.
[14] Guspi F. Frequency-domain reduction of potential field measurements to a horizontal plane[J]. Geoexploration, 1987,24(2):87-98.
doi: 10.1016/0016-7142(87)90083-4
[15] Pilkington M, Thurston J B. Draping corrections for aeromagnetic data: Line versus grid-based approaches[J]. Exploration Geophysics, 2001,32(2):95-101.
doi: 10.1071/EG01095
[16] 程振炎. 重磁场的有限元法曲化平[J]. 物探与化探, 1981,5(3):153-158.
[16] Cheng Z Y. The finite element method of gravity and magnetic field[J]. Geophysical and Geochemical Exploration, 1981,5(3):153-158.
[17] Xu S Z. The boundary element method in geophysics[M]. Tulsa: Society of Exploration Geophysicists, 2001:63-67.
[18] Xu S Z, Yang C H, Dai S K, et al. A new method for continuation of 3D potential fields to a horizontal plane[J]. Geophysics, 2003,68(6):1917-1921.
doi: 10.1190/1.1635045
[19] 徐世浙. 位场延拓的积分—迭代法[J]. 地球物理学报, 2006,49(4):1176-1182.
[19] Xu S Z. The interpolation-iteration method for continuation of potential fields[J]. Chinese Journal of Geophysics, 2006,49(4):1176-1182.
[20] 刘东甲, 洪天求, 廖旭涛, 等. 位场曲化平积分方程的迭代解[J]. 地球物理学报, 2012,55(10):3467-3476.
[20] Liu D J, Hong T Q, Liao X T, et al. Iterative solution of integral equation for potential field continuation from irregular surface to a horizontal plane[J]. Chinese Journal of Geophysics, 2012,55(10):3467-3476.
[21] 徐世浙, 余海龙. 位场曲化平的插值—迭代法[J]. 地球物理学报, 50(6):193-197.
[21] Xu S Z, Yu H L. The interpolation-iteration method for potential fields continuation from undulating surface to plane[J]. Chinese Journal of Geophysics, 2007,50(6):193-197.
[22] 郭良辉, 孟小红, 石磊, 等. 优化滤波方法及其在中国大陆布格重力异常数据处理中的应用[J]. 地球物理学报, 2012,55(12):4078-4088.
[22] Guo L H, Meng X H, Shi L, et al. Preferential filtering method and its application to Bouguer gravity anomaly of Chinese continent[J]. Chinese Journal of Geophysics, 2012,55(12):4078-4088.
[1] 杨荣祥, 王万银, 蔡梦轲, 王丁丁, 罗新刚. 基于卫星重力异常的渤海盆地秦南凹陷及邻区构造格局研究[J]. 物探与化探, 2023, 47(3): 584-596.
[2] 王润生, 武斌, 张海瑞, 于嘉宾, 董彦龙, 郭国强, 康一鸣. 山东省临沂凸起东北缘重力场特征及大地构造单元边界讨论[J]. 物探与化探, 2023, 47(2): 279-289.
[3] 孙争, 王俊, 丁鹏, 谭鑫. 一种重力异常向上延拓高度最优化确定方法[J]. 物探与化探, 2023, 47(1): 162-170.
[4] 石泽玉, 张志厚, 刘鹏飞, 范祥泰. 重力及其梯度异常正演的Moving-footprint大尺度模型分解方法[J]. 物探与化探, 2022, 46(3): 576-584.
[5] 王蓉, 熊杰, 刘倩, 薛瑞洁. 基于深度神经网络的重力异常反演[J]. 物探与化探, 2022, 46(2): 451-458.
[6] 冯旭亮, 魏泽坤. 基于界面反演增强的位场边缘识别方法[J]. 物探与化探, 2022, 46(1): 130-140.
[7] 邢锦程, 袁炳强, 张春灌, 冯旭亮, 段瑞锋, 薛健, 贾洪杨, 李想. 特立尼达盆地重力场特征及油气远景[J]. 物探与化探, 2021, 45(6): 1606-1616.
[8] 高维强, 史朝洋, 张利明, 冯旭亮. 山区地形改正密度逐次回归选取方法[J]. 物探与化探, 2021, 45(6): 1530-1538.
[9] 王正科, 卢琳, 刘花婷. 一种实测重力异常区域场的消除方法[J]. 物探与化探, 2021, 45(6): 1569-1577.
[10] 陈青, 孙帅, 丁成艺, 黄小宇, 陈浩, 申鹏, 罗港, 魏耀聪. iTilt-Euler法在重力数据处理及断裂解释中的应用[J]. 物探与化探, 2021, 45(6): 1578-1587.
[11] 王冠鑫, 罗锋, 周锡华, 闫方. 航空重力弱信号提取的卡尔曼滤波方法研究[J]. 物探与化探, 2021, 45(1): 76-83.
[12] 王学发, 王万银, 马杰, 杨敏, 梁建设, 邱春光, 王丁丁, 纪晓琳, 刘金兰. 基于重力资料的马达加斯加岛及邻区盆地构造单元分布特征[J]. 物探与化探, 2020, 44(4): 928-937.
[13] 许文强, 袁炳强, 刘必良, 姚长利. 多种重磁位场边缘识别方法及在南黄海北部断裂构造识别中的应用研究[J]. 物探与化探, 2020, 44(4): 962-974.
[14] 陈挺, 严迪, 杨剑, 廖国忠, 武斌, 冉中禹. 基于重力值标准差研究长宁6.0级地震活动性[J]. 物探与化探, 2020, 44(3): 480-488.
[15] 邱峰, 杜劲松, 陈超. 重力异常及其梯度张量DEXP定量解释方法的影响因素分析[J]. 物探与化探, 2020, 44(3): 540-549.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址:北京市学院路29号航遥中心 邮编:100083
电话:010-62060192;62060193 E-mail:whtbjb@sina.com