基于等效源的重力异常空间导数计算方法及其在川滇地区的应用

doi:10.12196/j.issn.1000-3274.2021.01.010

地震 ›› 2021, Vol. 41 ›› Issue (1): 129-140.doi: 10.12196/j.issn.1000-3274.2021.01.010

基于等效源的重力异常空间导数计算方法及其在川滇地区的应用

李红蕾^1,2,3, 韩建成^1,2, 陈石^1,2, 侍文^1,2, 张贝^1,2

1.中国地震局地球物理研究所, 北京 100081;
2.北京白家疃地球科学国家野外科学观测研究站, 北京 100095;
3.中国科学院精密测量科学与技术创新研究院大地测量与地球动力学国家重点实验室, 湖北武汉 430077

收稿日期:2020-10-23 修回日期:2020-11-13 出版日期:2021-01-31 发布日期:2021-01-28
通讯作者: 韩建成,助理研究员。E-mail:jchhan@cea-igp.ac.cn
作者简介:李红蕾(1987-),女,山东聊城人,助理研究员,主要从事重力数据处理、反演与解释研究。
基金资助:
科技部重点研发专项项目(2018YFC0603502, 2018YFC1504506, 2016YFC0601101); 国家自然科学基金青年项目(42004074); 国家自然科学基金面上项目(41774090); 国家自然科学基金地震联合基金项目(U1939205, U1839210); 大地测量与地球动力学重点实验室基金(SKLGED2019-1-1-E); 中国博士后科学基金(2018M641424); 中国地震科学实验场项目(2019CSES0105); 中国地震局地球物理研究所基本科研业务费项目(DQJB20Z04, DQJB19B03)

Calculative Method of Spatial Derivatives of Gravity Anomaly based on Equivalent Source and Its Application in Sichuan-Yunnan Region

LI Hong-lei^1,2,3, HAN Jian-cheng^1,2, CHEN Shi^1,2, SHI Wen^1,2, ZHANG Bei^1,2

1. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China;
2. Beijing Baijiatuan Earth Sciences National Observation and Research Station, Beijing 100095, China;
3. State Key Laboratory of Geodesy and Earth's Dynamics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430077, China

Received:2020-10-23 Revised:2020-11-13 Online:2021-01-31 Published:2021-01-28

摘要/Abstract

摘要： 在综合分析前人研究成果的基础上, 采用基于贝叶斯准则的等效源模型滤波增强方法来识别地质体边界。通过理论模型实验讨论了在相同的噪声水平情况下, 不同算法对于计算结果的影响, 模型实验结果表明, 与直接计算和传统高斯滤波增强技术相比, 新算法可以有效消除传统算法中普遍存在的对于噪声的过度敏感问题, 所得到的地质体边界形态失真更小。随后, 利用此方法对川滇地区的重力数据进行了处理, 并对该地区的重力水平总梯度(THD)、重力垂直梯度模(AVG)、重力解析信号模量(AAS)、重力水平梯度垂向导数(VGHD)应用效果进行了对比分析。研究结果提供了川滇地区的构造分区信息, 对该地区构造演化、边界定位及动力学过程分析具有重要的指导意义。

关键词: 等效源, 滤波增强, 贝叶斯优化, 川滇地区, 边界特征

Abstract: Based on the comprehensive analysis of previous research results, an equivalent source model filtering and enhancement method based on Bayesian criteria is used to identify the boundary of geological bodies. Through the theoretical model experiment,we discussed the noise level in the same circumstances, the influence of different algorithms for the calculation results. The model experimental results show that compared with the direct calculation and the traditional Gaussian filtering enhancement technique, the new algorithm can effectively eliminate the traditional algorithm in common to the problem of noise of overly sensitive geological boundary shape distortion smaller. Subsequently, this method was used to process the gravity data in Sichuan-Yunnan region, and the application effects of amplitude of vertical gradient (AVG), total horizontal gravity gradient (THD), amplitude of analytic signal (AAS), and vertical gradient of total horizontal derivative (VGHD) of this region were compared and analyzed. The results of this study provide the information of tectonic zoning in the Sichuan-Yunnan region, and have important guiding significance for tectonic evolution, seismic distribution and related dynamic analysis for this region.

Key words: The equivalent source, Filtering enhancement, Bayesian optimization, Sichuan and Yunnan region, Boundary characteristics

中图分类号:

P315.7

李红蕾, 韩建成, 陈石, 侍文, 张贝. 基于等效源的重力异常空间导数计算方法及其在川滇地区的应用[J]. 地震, 2021, 41(1): 129-140.

LI Hong-lei, HAN Jian-cheng, CHEN Shi, SHI Wen, ZHANG Bei. Calculative Method of Spatial Derivatives of Gravity Anomaly based on Equivalent Source and Its Application in Sichuan-Yunnan Region[J]. EARTHQUAKE, 2021, 41(1): 129-140.

参考文献

[1] Pilkington M, Keating P B. The utility of potential field enhamcements for remote predictive mapping[J]. Canadian Journal of Remote Sensing, 2009, 35(sup1): S1-S11.
[2] Fajklewicz Z J. Gravity vertical gradient measurements for the detection of small geologic and anthropogenic forms[J]. Geophysics, 1976, 41(5): 1016-1030.
[3] Grauch V J S, Cordell L. Limitations of determining density or magnetic boundaries from the horizontal gradient of gravity or pseudogravity data[J]. Geophysics, 1987, 52(1): 118-121.
[4] Wijns C, Perez C, Kowalczyk P. Theta map: Edge detection in magnetic data[J]. Geophysics, 2005, 70(4): L39-L43.
[5] Miller H G, Singh V. Potential field tilt—a new concept for location of potential field sources[J]. Journal of Applied Geophysics, 1994, 32(2-3): 213-217.
[6] Roest W R, Verhoef J, Pilkington M. Magnetic interpretation using the 3-D analytic signal[J]. Geophysics, 1992, 57(1): 116-125.
[7] Cooper G R J, Cowan D R. Edge enhancement of potential-field data using normalized statistics[J]. Geophysics, 2008, 73(3): H1-H4.
[8] Nabighian M N, Ander M E, Grauch V J S, et al. Historical development of the gravity method in exploration [J]. Geophysics, 2005, 70(6): 63ND-89ND.
[9] 胡聪. 利用重力资料进行边界和质心提取的方法研究及应用[D]. 长春: 吉林大学, 2012.
HU Cong. Study and application of boundary and mass center extraction with the data of gravity[D]. Changchun: Jilin University, 2012 (in Chinese).
[10] 马国庆. 位场(重&磁)及其梯度异常自动解释方法研究[D]. 长春: 吉林大学, 2013.
MA Guo-qing. The study on the automatic interpretation methods of potential field (gravity &magnetic) and its gradients[D]. Changchun: Jilin University, 2013 (in Chinese).
[11] 张凤旭, 张凤琴, 刘财, 等. 断裂构造精细解释技术三方向小子域滤波[J]. 地球物理学报, 2007, 50(5): 1543-1550.
ZHANG Feng-xu, ZHANG Feng-qin, LIU Cai, et al. A technique for elaborate explanation of faulted structures: Three-directional small subdomain filtering[J]. Chinese Journal of Geophysics, 2007, 50(5): 1543-1550 (in Chinese).
[12] 王彦国, 张瑾, 张凤旭. 基于迭代差分的视导数及其在边界识别中的应用[J]. 地球物理学进展, 2016, 31(2): 794-804.
WANG Yan-guo, ZHANG Jin, ZHANG Feng-xu. Apparent derivative based on iterative differential and its application to edge recognition of potential fields[J]. Progress in Geophysics, 2016, 31(2): 794-804 (in Chinese).
[13] Hansen R O, Miyazaki Y. Continuation of potential fields between arbitrary surfaces[J]. Geophysics, 1984, 49(6): 787-795.
[14] 王万银, 潘作枢. 偶层位曲面位场数据处理及转换方法[J]. 西安工程学院学报, 1996, 18(3): 69-76.
WANG Wan-yin, PAN Zuo-shu. Data processing and transformation methods of potential fields on an arbitrary surface with dipole layer potential technique[J]. Journal of Xi'an College of Geology, 1996, 18(3): 69-76 (in Chinese).
[15] Akaike H. Information theory and an extension of the maximum likelihood principle[M]∥Selected papers of Hirotugu Akaike. New York: Springer, 1998: 199-213.
[16] Nagy D, Papp G, Benedek J. The gravitational potential and its derivatives for the prism[J]. Journal of Geodesy, 2000, 74(7-8): 552-560.
[17] 陈石, 徐伟民, 蒋长胜. 中国大陆西部重力场变化与强震危险性关系[J]. 地震学报, 2015, 37(4): 575-587.
CHEN Shi, XU Wei-min, JIANG Chang-sheng. Relationship between gravity variation and seismic hazards in the Western China[J]. Acta Seismological Sinica, 2015, 37(4): 575-587 (in Chinese).
[18] Wang W, Pan Y, Qiu Z. A new edge recognition technology based on the normalized vertical derivative of the total horizontal derivative for potential field data[J]. Applied Geophysics, 2009, 6(3): 226-233.
[19] Phillips J D. Geosoft executable (GX's) developed by the US Geological Survey, version 2.0, with notes on GX development from Fortran code[R]. US Geological Survey, 2007.
[20] 闻学泽, 张培震, 杜方, 等. 2008年汶川8.0级地震发生的历史与现今地震活动背景[J]. 地球物理学报, 2009, 52(2): 444-454.
WEN Xue-ze, ZHANG Pei-zhen, DU Fang, et al. The background of historical and modern seismic activities of the occurrence of the 2008 M_S8.0 Wenchuan, Sichuan, earthquake[J]. Chinese Journal of Geophysics, 2009, 52(2): 444-454 (in Chinese).
[21] 孟小红, 石磊, 郭良辉, 等. 青藏高原东北缘重力异常多尺度横向构造分析[J]. 地球物理学报, 2012, 55(12): 3933-3941.
MENG Xiao-hong, SHI Lei, GUO Liang-hui, et al. Multi-scale analyses of transverse structures based on gravity anomalies in the northeastern margin of the Tibetan Plateau[J]. Chinese Journal of Geophysics, 2012, 55(12): 3933-3941 (in Chinese).
[22] 李媛媛, 杨宇山. 位场梯度的归一化标准差方法在地质体边界定位问题中的应用[J]. 地质科技情报, 2009, 28(5): 138-142.
LI Yuan-yuan, YANG Yu-shan. Derivative-based normalized standard deviation of potential field data in Geological contact mapping[J]. Geological Science and Technology Information, 2009, 28(5): 138-142 (in Chinese).
[23] 祝意青, 刘芳, 李铁明, 等. 川滇地区重力场动态变化及其强震危险含义[J]. 地球物理学报, 2015, 58(11): 4187-4196.
ZHU Yi-qing, LIU Fang, LI Tie-ming, et al. Dynamic variation of the gravity field in the Sichuan-Yunnan region and its implication for seismic risk[J]. Chinese Journal of Geophysics, 2015, 58(11): 4187-4196 (in Chinese).
[24] 楼海, 王椿镛. 川滇地区重力异常的小波分解与解释[J]. 地震学报, 2005, 27(5): 515-523.
LOU Hai, WANG Chun-yong. Wavelet analysis and interpretation of gravity data in Sichuan-Yunnan region, China[J]. Acta Seismologica Sinica, 2005, 27(5): 515-523 (in Chinese).

基于等效源的重力异常空间导数计算方法及其在川滇地区的应用

Calculative Method of Spatial Derivatives of Gravity Anomaly based on Equivalent Source and Its Application in Sichuan-Yunnan Region

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	张熔鑫, 邢会林, 舒涛, 刘骏标, 郭志伟, 王建超, 谭玉阳. 中国地震科学实验场断层系统的三维计算网格模型[J]. 地震, 2023, 43(3): 18-33.
[2]	张贝, 陈石, 李红蕾, 杨锦玲, 韩建成, 卢红艳. 时变重力场球面模型反演算法和模拟实验[J]. 地震, 2021, 41(1): 13-24.
[3]	王林海, 陈石, 张贝, 卢红艳. 2013年芦山M_S7.0地震前重力场源变化分析[J]. 地震, 2021, 41(1): 78-92.
[4]	徐伟民, 陈石, 阮明明, 杨锦玲, 韩建成, 李红蕾, 卢红艳. 川滇地区陆地流动重力测网场源分辨能力评估[J]. 地震, 2021, 41(1): 191-204.
[5]	郑秋月, 王青华, 刘东, 黄江培, 陈石. 基于时变重力数据的川滇地区场源模型反演和解释[J]. 地震, 2021, 41(1): 205-218.
[6]	陈学忠, 李艳娥. 川滇地区地震活动与地球自转速率变化之间的关系[J]. 地震, 2019, 39(1): 126-135.
[7]	巩浩波,李光科,廖欣,陈敏. 川滇地区井水位潮汐响应特性分析[J]. 地震, 2017, 37(1): 20-30.
[8]	王俊, 邵志刚, 孙小龙, 王熠熙, 向阳, 王喜龙. 川滇地区强震前流体异常特征与预测指标体系初探[J]. 地震, 2016, 36(4): 109-119.
[9]	郑秋月, 陈石, 李晓一, 徐伟民. 设定场源参数的区域重力微变化特征与模型实验[J]. 地震, 2016, 36(4): 186-195.
[10]	刘月, 吕晓健, 田勤俭. 基于“区域-时间-长度算法”分析川滇地区强震前地震活动性变化^*[J]. 地震, 2016, 36(2): 94-104.
[11]	魏文薪, 江在森, 刘晓霞, 邹镇宇. 川滇地区应变率场分布及其变化特征研究[J]. 地震, 2015, 35(4): 11-20.
[12]	董曼, 杨天青, 郑通彦, 陈通. 川滇地区走滑型地震烈度衰减特征研究[J]. 地震, 2015, 35(3): 147-154.
[13]	何康, 晏锐, 王春, 郑海刚, 李军辉, 方震. 基于Wigner-Ville分布的地磁前兆异常提取方法及其在川滇地区的应用[J]. 地震, 2014, 34(4): 88-99.
[14]	杨文, 刘杰. 2012年4月11日苏门答腊8.6级地震前后川滇地区波速变化研究[J]. 地震, 2013, 33(4): 248-256.
[15]	康春丽. 川滇地区长波辐射场变化与地震活动关系研究[J]. 地震, 2008, 28(3): 43-48.