时变重力场球面模型反演算法和模拟实验

doi:10.12196/j.issn.1000-3274.2021.01.002

地震 ›› 2021, Vol. 41 ›› Issue (1): 13-24.doi: 10.12196/j.issn.1000-3274.2021.01.002

时变重力场球面模型反演算法和模拟实验

张贝^1,2, 陈石^1,2, 李红蕾^1,2, 杨锦玲^1,3, 韩建成^1,2, 卢红艳^1,2

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

收稿日期:2020-07-16 修回日期:2020-09-03 出版日期:2021-01-31 发布日期:2021-01-28
通讯作者: 陈石,研究员。E-mail:chenshi@cea-igp.ac.cn
作者简介:张贝(1985-),男,河南商丘人,助理研究员,主要从事计算地球动力学与数值模拟。
基金资助:
科技部重点研发专项项目(2018YFC1504506; 2018YFC0603502), 国家自然科学基金面上项目(41774090), 国家自然科学基金地震联合基金项目(U1939205), 中国地震局地球物理研究所基本科研业务费专项资助 (DQJB20X09), 中国地震科学试验场专项(2019CSES0105)联合资助

Time-varying Gravity Inversion using Spherical Model and Synthetic Test

ZHANG Bei^1,2, CHEN Shi^1,2, LI Hong-lei^1,2, YANG Jin-ling^1,3, HAN Jian-cheng^1,2, LU Hong-yan^1,2

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

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

摘要/Abstract

摘要： 时变重力场是研究地球内部介质物性变化的重要手段。本文提出了一种适用于地面流动重力测量获得的时变重力信号的场源反演方法, 该方法采用球坐标系下的六面体单元来模拟场源介质, 适合大尺度地震流动重力测量数据的等效源模型构建。通过引入重力时变信号的一阶光滑先验条件, 压制了时变重力信号中的短周期高频分量, 可用于提取与地震孕育相关的长周期信号。通过理论和模型实验证明了本文算法的可靠性和稳定性, 并使用南北地震带南段2014—2017年的流动重力实测数据进行了反演解释, 获得了地壳内部等效场源的视密度时变信号, 变化量级在正常地壳密度的±0.7‰之间, 其空间形态受川滇菱形块体边界控制。研究成果可用于时变重力场模型解释和深部场源特征提取, 可为地震重力前兆信号分析和相关研究提供完备的方法保障。

关键词: 时变重力场, 位场反演, 流动重力, 等效源, 密度变化

Abstract: Time-varying gravity field is an important means to study the physical properties of the earth. In this paper, a field source inversion method which is suitable for the time-varying gravity signals obtained from ground flow gravity measurements is proposed. The method uses hexahedral elements in spherical coordinate system to simulate the field source media and is suitable for constructing equivalent source models of large-scale seismic flow gravity measurements. In this method, the high-frequency components of time-varying gravity signals are suppressed by introducing a first-order smoothing priori, which can be used to extract long-period signals related to earthquake preparation. The theory and model experiments prove that the algorithm is reliable and stable. We also used the campaign gravity survey of the southern section of the north-south seismic zone from 2014 to 2017 for inversion. And the time-varying apparent density signal of the equivalent field source in the crust was obtained, which was within ±0.7‰ of the normal crust density, and its spatial morphology was controlled by the boundary of the Sichuan-Yunnan block. The results of this paper can be applied to the interpretation of time-varying gravity field model and the extraction of deep field source characteristics, which can provide a complete method for seismic gravity precursor signal analysis and related research.

Key words: Time-varying gravity field, Potential field inversion, Campaign gravity survey, Equivalent source, Density variation

中图分类号:

P315.7

张贝, 陈石, 李红蕾, 杨锦玲, 韩建成, 卢红艳. 时变重力场球面模型反演算法和模拟实验[J]. 地震, 2021, 41(1): 13-24.

ZHANG Bei, CHEN Shi, LI Hong-lei, YANG Jin-ling, HAN Jian-cheng, LU Hong-yan. Time-varying Gravity Inversion using Spherical Model and Synthetic Test[J]. EARTHQUAKE, 2021, 41(1): 13-24.

参考文献

[1] Tapley B D, Bettadpur S, Ries J C, et al. GRACE measurements of mass variability in the Earth system[J]. Science, 2004, 305(5683): 503-505.
[2] Seo K W, Wilson C R, Famiglietti J S, et al. Terrestrial water mass load changes from Gravity Recovery and Climate Experiment (GRACE)[J]. Water Resources Research, 2006, 42(5): 173-180.
[3] Syed T H, Famiglietti J S, Rodell M, et al. Analysis of terrestrial water storage changes from GRACE and GLDAS[J]. Water Resources Research, 2008, 44(2): 339-356.
[4] Han S C, Shum C K, Bevis M, et. al. Crustal dilatation observed by GRACE after the 2004 Sumatra-Andaman earthquake[J]. Science, 2006, 313(5787): 658-662.
[5] Steffen H, Wu P, Wang H S. Optimal locations for absolute gravity measurements and sensitivity of GRACE observations for constraining glacial isostatic adjustment on the northern hemisphere[J]. Geophysical Journal International, 2012, 190(3): 1483-1494.
[6] Ivins E R, James T S, Wahr J, et al. Antarctic contribution to sea level rise observed by GRACE with improved GIA correction[J]. Journal of Geophysical Research: Solid Earth, 2013, 118(6): 3126-3141.
[7] Crossley D, Hinderer J, Riccardi U. The measurement of surface gravity[J]. Reports on Progress in Physics, 2013, 76(4): 046101.
[8] Young W M, Lumley D. Feasibility analysis for time-lapse seafloor gravity monitoring of producing gas fields in the Northern Carnarvon Basin, offshore Australia[J]. Geophysics, 2015, 80(2): WA149-WA160.
[9] Vigouroux N, Williams-Jones G, Chadwick W, et al. 4D gravity changes associated with the 2005 eruption of Sierra Negra volcano, Galápagos[J]. Geophysics, 2008, 73(6): WA29-WA35.
[10] 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.
[11] Tanaka Y, Kato A, Sugano T, et al. Gravity changes observed between 2004 and 2009 near the Tokai slow-slip area and prospects for detecting fluid flow during future slow-slip events[J]. Earth Planets and Space, 2010, 62(12): 905-913.
[12] 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.
[13] Mouyen M, Chao B F, Hwang C, et al. Gravity monitoring of Tatun Volcanic Group activities and inference for underground fluid circulations[J]. Journal of Volcanology and Geothermal Research, 2016, 328: 45-58.
[14] Van Camp M, de Viron O, Watlet A, et al. Geophysics from terrestrial time-variable gravity measurements[J]. Reviews of Geophysics, 2017, 55: 938-992.
[15] 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.
[16] 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.
[17] 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.
[18] Heck B, Seitz K. A comparison of the tesseroid, prism and point-mass approaches for mass reductions in gravity field modelling[J]. Journal of Geodesy, 2007, 81(2): 121-136.
[19] Li Y G, Oldenburg D W. 3-D inversion of gravity data[J]. Geophysics, 1998, 63(1): 109-119.
[20] Schrama E J O, Wouters B, Lavallée D A. Signal and noise in Gravity Recovery and Climate Experiment (GRACE) observed surface mass variations[J]. Journal of Geophysical Research: Solid Earth, 2007, 112(B8): B08407.
[21] Akaike H. Likelihood and the Bayes procedure[A]. In: Bernardo JM, DeGroot MH, Lindley DV, Smith AFM (eds) Bayesian statistics[M]. University Press, Valencia, 1980, 143-166.

时变重力场球面模型反演算法和模拟实验

Time-varying Gravity Inversion using Spherical Model and Synthetic Test

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 8

编辑推荐

Metrics

本文评价

[1]	王林海, 陈石, 张贝, 卢红艳. 2013年芦山M_S7.0地震前重力场源变化分析[J]. 地震, 2021, 41(1): 78-92.
[2]	李红蕾, 韩建成, 陈石, 侍文, 张贝. 基于等效源的重力异常空间导数计算方法及其在川滇地区的应用[J]. 地震, 2021, 41(1): 129-140.
[3]	杨锦玲, 陈石, 李红蕾, 张贝, 阮明明. 华南地区时变重力场建模实验和异常分析[J]. 地震, 2021, 41(1): 141-152.
[4]	贾晓东, 翟丽娜, 陈石. 辽宁地区时变重力场源变化特征分析[J]. 地震, 2021, 41(1): 180-190.
[5]	刘芳, 祝意青, 梁伟锋, 杨姣, 赵云峰. 基于欧拉反褶积方法估算华北地区重力变化的等效源参数[J]. 地震, 2016, 36(4): 163-170.
[6]	徐伟民, 孙少波, 李晓一, 卢红艳, 郑秋月. 基于变差函数的插值方法计算华北地区重力场变化[J]. 地震, 2016, 36(4): 171-185.
[7]	郑秋月, 陈石, 李晓一, 徐伟民. 设定场源参数的区域重力微变化特征与模型实验[J]. 地震, 2016, 36(4): 186-195.
[8]	郑兵, 郝洪涛, 王双洪, 李菲菲, 马伶俐, 蒲晓霞, 蒋海涛, 苏琴. 四川地区1987—2012年流动重力资料的初步整理[J]. 地震, 2014, 34(1): 134-143.