华南地区时变重力场建模实验和异常分析

doi:10.12196/j.issn.1000-3274.2021.01.011

地震 ›› 2021, Vol. 41 ›› Issue (1): 141-152.doi: 10.12196/j.issn.1000-3274.2021.01.011

华南地区时变重力场建模实验和异常分析

杨锦玲^1,2, 陈石^1,3, 李红蕾^1,3, 张贝^1,3, 阮明明⁴

1.中国地震局地球物理研究所, 北京 100081;
2.福建省地震局, 福建福州 350003;
3.北京白家疃国家地球科学野外观测研究站, 北京 100095;
4.中国地震局地球物理勘探中心, 河南郑州 450002

收稿日期:2020-07-16 修回日期:2020-09-18 出版日期:2021-01-31 发布日期:2021-01-28
通讯作者: 陈石,研究员。E-mail:chenshi@cea-igp.ac.cn
作者简介:杨锦玲(1986-),女,福建泉州人,博士研究生,主要从事时变重力位场数据处理与反演解释研究。
基金资助:
科技部重点研发专项项目(2018YFC1504506; 2018YFC0603502), 国家自然科学基面上项目(41774090), 地震联合基金项目(U1939205)、中国地震局地球物理研究所基本科研业务费专项资助 (DQJB19B03)、中国地震局地震科技星火计划项目(XH18023Y)和中国地震局震情跟踪定向工作任务青年项目(2020010210)联合资助

Modeling Experiment and Anomaly Analysis of Time-varying Gravity Field in South China

YANG Jin-ling^1,2, CHEN Shi^1,3, LI Hong-lei^1,3, ZHANG Bei^1,3, RUAN Ming-ming⁴

1. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China;
2. Fujian Earthquake Agency, Fuzhou 350003, China;
3. Beijing Baijiatuan Earth Science National Observation and Research Station, Beijing 100095, China;
4. Geophysical Exploration Center, China Earthquake Administration, Zhengzhou 450002, China

Received:2020-07-16 Revised:2020-09-18 Online:2021-01-31 Published:2021-01-28

摘要/Abstract

摘要： 陆面时变重力测量是监测地壳内部密度变化和物质运移的重要手段。为确定华南时变重力观测网络的场源监测能力和重力场变化特征, 本文采用球面六面体单元构建重力场模型, 开展重力场建模实验, 对比不同建模方法与噪声条件下的局部重力场恢复效果, 并对2015—2017年来5期实测流动重力观测数据进行计算和分析。结果表明, 最小二乘配置方法的建模效果和抗噪声干扰能力较好。华南重力测网在地表等效可观测信号约±30×10^-8 m/s²的条件下, 场源空间分辨能力约为55 km; 实际重建的2015至2017年重力场累积变化基本趋势较一致, 重力负异常在不同时间尺度下表现为持续增强和积累, 且在局部出现正负波动的周期性变化。本文研究结果对评价陆地重力数据质量, 指导地震重力观测网络系统改造升级, 客观地开展重力场变化解释、分析强震孕育和物质变迁等科学问题具有重要的参考价值。

关键词: 时变重力场, 位场模型, 检测板实验, 最小二乘配置, 华南地区

Abstract: The terrestrial time-varying gravity measurement has been used for monitoring the density changes and tectonic movements associated with earthquakes. In order to determine the field source monitoring capability and gravity field characteristics of South China region, a gravity model was constructed based on spherical tesseroid units. To assess the quality of recovered gravity field by different modeling methods, gravity field source modeling experiments and noise tests were carried out. The variation of observed gravitational field was calculated and recovered from 2015 to 2017 using least square collocation method. Overall, it can be concluded that the least squares configuration method performs better gravity field recovering effect and anti-noise ability. The spatial resolution capability of the South China Gravity Survey Network is about 55 km, and the equivalent observable signal is about 30×10^-8 m/s². The accumulative temporal and spatial characteristics of the gravity field between 2015 and 2017 is actually consistent. Negative gravity anomalies continue to increase and accumulate at different time scales, and there are periodic changes in positive and negative fluctuations locally. The results have important value for evaluating the quality of terrestrial gravity data, guiding the transformation and upgrading of gravity observation network, while advancing our standing of the interpretation of gravity field changes, analyzing strong earthquakes and geophysical field changes.

Key words: Time-varying gravity field, Potential field model, Checkerboard tests, Least square collocation, South China region

中图分类号:

P315.7

杨锦玲, 陈石, 李红蕾, 张贝, 阮明明. 华南地区时变重力场建模实验和异常分析[J]. 地震, 2021, 41(1): 141-152.

YANG Jin-ling, CHEN Shi, LI Hong-lei, ZHANG Bei, RUAN Ming-ming. Modeling Experiment and Anomaly Analysis of Time-varying Gravity Field in South China[J]. EARTHQUAKE, 2021, 41(1): 141-152.

参考文献

[1] Bagnardi M, Poland M P, Carbone D, et al. Gravity changes and deformation at Kīlauea Volcano, Hawaii, associated with summit eruptive activity, 2009—2012[J]. Journal of Geophysical Research: Solid Earth, 2014, 119(9): 7288-7305.
[2] Carbone D, Poland M P, Diament M, et al. The added value of time-variable microgravimetry to the understanding of how volcanoes work[J]. Earth-Science Reviews, 2017, 169: 146-179.
[3] Hector B, Séguis L, Hinderer J, et al. Water storage changes as a marker for base flow generation processes in a tropical humid basement catchment (Benin): Insights from hybrid gravimetry[J]. Water Resources Research, 2015, 51(10): 8331-8361.
[4] Kennedy J, Ferré Ty P A, Güntner A, et al. Direct measurement of subsurface mass change using the variable baseline gravity gradient method[J]. Geophysical Research Letters, 2014, 41(8): 2827-2834.
[5] Kennedy J R, Ferré Ty P A. Accounting for time- and space-varying changes in the gravity field to improve the network adjustment of relative-gravity data[J]. Geophysical Journal International, 2016, 204(2): 892-906.
[6] 佘雅文, 付广裕, 王灼华, 等.重力与地形数据揭示的巴颜喀拉块体东缘垂向构造应力场[J]. 地球物理学报, 2017, 60(6): 2480-2492.
SHE Ya-wen, FU Guang-yu, WANG Zhuo-hua, et al. Vertical tectonic stress in eastern margin of Bayan Har block revealed by gravity and terrain data[J]. Chinese Journal of Geophysics, 2017, 60(6): 2480-2492.
[7] Chen S, Liu M, Xing L L, et al. Gravity increase before the 2015 M_W7.8 Nepal earthquake[J]. Geophysical Research Letters, 2016, 43(1): 111-117.
[8] Zhu Y Q, Zhan F B, Zhou J C, et al. Gravity Measurements and Their Variations before the 2008 Wenchuan Earthquake[J]. Bulletin of the Seismological Society of America, 2010, 100(5B): 2815-2824.
[9] Zhang Y, Chen S, Xing L L, et al. Gravity changes before and after the 2008 M_W7.9 Wenchuan earthquake at Pixian absolute gravity station in more than a decade[J]. Pure and Applied Geophysics, 2020, 177: 121-133.
[10] 申重阳, 李辉, 孙少安, 等. 重力场动态变化与汶川M_S8.0地震孕育过程[J]. 地球物理学报, 2009, 52(10): 2547-2557.
SHEN Chong-yang, LI Hui, SUN Shao-an, et al. Dynamic variations of gravity and the preparation process of the Wenchuan M_S8.0 earthquake[J]. Chinese Journal of Geophysics, 2009, 52(10): 2547-2557.
[11] Xuan S B, Jin S G, Chen Yong, et al. Insight into the preparation of the 2016 M_S6.4 Menyuan earthquake from terrestrial gravimetry-derived crustal density changes[J]. Scientific Reports, 2019, 9: 18227.
[12] 祝意青, 张勇, 张国庆, 等. 21世纪以来青藏高原大震前重力变化[J]. 科学通报, 2020, 65(7): 622-632.
ZHU Yi-qing, ZHANG Yong, ZHANG Guo-qing, et al. Gravity variations preceding the large earthquakes in the Qinghai-Tibet Plateau from 21st century[J]. Chinese Science Bulletin, 2020, 65(7): 622-632 (in Chinese).
[13] Moritz H. Least-squares collocation[J]. Reviews of geophysics, 1978, 16(3): 421-430.
[14] 阮明明, 陈石, 韩建成. 用最小二乘配置法构建局部重力场模型[J]. 地震学报, 2020, 42(1): 53-65.
RUAN Ming-ming, CHEN Shi, HAN Jian-Cheng. Regional gravity field model constructed by the least squares collocation[J]. Acta Seismologica Sinica, 2020, 42(1): 53-65 (in Chinese).
[15] 杜劲松, 陈超, 梁青, 等. 球冠体积分的重力异常正演方法及其与Tesseroid单元体泰勒级数展开方法的比较[J]. 测绘学报, 2012, 41(3): 339-346.
DU Jing-song, CHEN Chao, LIANG Qing, et al. Gravity anomaly calculation based on volume integral in spherical cap and comparison with the Tesseroid-Taylor series expansion approach[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(3): 339-346 (in Chinese).
[16] Uieda L, Barbosa V C F, Braitenberg C. Tesseroids: Forward-modeling gravitational fields in spherical coordinates[J]. Geophysics, 2016, 81(5): F41-F48.
[17] 胡敏章, 李辉, 刘子维, 等. 川滇地区2010—2013年重力变化及重力网的地震监测能力[J]. 大地测量与地球动力学, 2015, 35(4): 616-620.
HU Min-zhang, LI Hui, LIU Zi-wei, et al. The gravity change over Sichuan-Yunnan region in 2010-2013 and the earthquake monitoring ability of the gravimetric network[J]. Journal of Geodesy and Geodynamics, 2015, 35(4): 616-620 (in Chinese).
[18] 柯灏. 海洋无缝垂直基准构建理论和方法研究[D]. 武汉: 武汉大学, 2012.
KE Hao. Research on the theory and implementation method of marine seamless vertical reference surface[D]. Wuhan: Wuhan University, 2012 (in Chinese).
[19] 彭彬, 周艳莲, 高苹, 等. 气温插值中不同空间插值方法的适用性分析以江苏省为例[J]. 地球信息科学学报, 2011, 13(4): 539-548.
PENG Bin, ZHOU Yan-lian, GAO Ping, et al. Suitability Assessment of Different Interpolation Methods in the Gridding Process of Station Collected Air Temperature: A Case Study in Jiangsu Province, China[J]. Journal of Geo-Information Science, 2011, 13(4): 539-548 (in Chinese)
[20] Wu Y Q, Jiang Z S, Yang G H, et al. Comparison of GPS strain rate computing methods and their reliability[J]. Geophysical Journal International, 2011, 185(2): 703-717.
[21] 江在森, 刘经南. 应用最小二乘配置建立地壳运动速度场与应变场的方法[J]. 地球物理学报, 2010, 53(5): 1109-1117.
JIANG Zai-sen, LIU Jing-nan. The method in establishing strain field and velocity field of crustal movement using least squares collocation[J]. Chinese Journal of Geophysics, 2010, 53(5): 1109-1117 (in Chinese).
[22] 欧阳永忠. 海空重力测量数据处理关键技术研究[D]. 武汉: 武汉大学, 2013.
OUYANG Yong-zhong. On key technologies of data processing for air-sea gravity surveys[D]. Wuhan: Wuhan University, 2013 (in Chinese).
[23] Chen S, Zhuang J C, Li X Y, et al. Bayesian approach for network adjustment for gravity survey campaign: methodology and model test[J]. Journal of Geodesy, 2019, 93(5): 681-700.
[24] 高倩, 陈石. 基于欧拉反褶积方法计算川滇交界重力变化场源特征[J]. 地球物理学进展, 2015, 30(2): 503-509.
GAO Qian, CHEN Shi. Based on the euler deconvolution technique to determine the field source characteristics of gravity Variation in the boundary of Sichuan and Yunnan[J]. Progress in Geophysics, 2015, 30(2): 503-509 (in Chinese).

华南地区时变重力场建模实验和异常分析

Modeling Experiment and Anomaly Analysis of Time-varying Gravity Field in South China

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 10

编辑推荐

Metrics

本文评价

[1]	金阳, 金红林, 高原, 黄星. 基于GNSS数据的三种应变计算方法在青藏高原东南缘地区的应用[J]. 地震, 2021, 41(4): 57-67.
[2]	徐伟民, 陈石, 阮明明, 杨锦玲, 韩建成, 李红蕾, 卢红艳. 川滇地区陆地流动重力测网场源分辨能力评估[J]. 地震, 2021, 41(1): 191-204.
[3]	贾晓东, 翟丽娜, 陈石. 辽宁地区时变重力场源变化特征分析[J]. 地震, 2021, 41(1): 180-190.
[4]	张贝, 陈石, 李红蕾, 杨锦玲, 韩建成, 卢红艳. 时变重力场球面模型反演算法和模拟实验[J]. 地震, 2021, 41(1): 13-24.
[5]	魏文薪, 江在森, 刘晓霞, 邹镇宇. 川滇地区应变率场分布及其变化特征研究[J]. 地震, 2015, 35(4): 11-20.
[6]	熊小慧, 江在森, 武艳强, 刘晓霞, 魏文薪, 邹镇宇, 王启欣. 2008年汶川M_S8.0地震后GPS资料反映的龙门山断裂带南段地壳水平变形动态特征分析[J]. 地震, 2015, 35(3): 10-21.
[7]	章纯. 台湾地震对华南及华东地区应力场影响的数值模拟研究[J]. 地震, 2007, 27(3): 26-33.
[8]	顾瑾萍, 吕培苓, 彭美凤, 康春丽, 张凌空. 华南地区地震活动特征与台湾强震影响[J]. 地震, 2004, 24(1): 55-59.
[9]	刘特培, 秦乃岗, 郭贵安. 加卸载响应比理论在华南地区中强震的应用[J]. 地震, 2003, 23(1): 90-94.
[10]	江在森, 张希, 王双绪, 祝意青. 利用地形变观测量求解地壳水平应变场的方法^①[J]. 地震, 1999, 19(1): 41-48.