«上一篇/Previous Article|本期目录/Table of Contents|下一篇/Next Article»

《地震研究》[ISSN:1000-0666/CN:53-1062/P]

卷:

期数:

2025年04期

页码:

550-560

栏目:

重力场及其变化观测与地球科学应用专辑

出版日期:

2025-05-01

文章信息/Info

Title:

Joint Inversion Method for Seismic Body Wave Traveltime and Gravity Data Based on Structural Constraint

作者:

谭智¹; 5; 姜金钟¹; 2; 石磊³; 4; 李永华³; 4

(1.云南省地震局,云南 昆明 650224; 2.中国地震局昆明地震预报研究所,云南 昆明 650224; 3.中国地震局地球物理研究所,北京 100081; 4.中国地震局震源物理重点实验室,北京 100081; 5.云南省地震局 下关地震监测中心站,云南 大理 671000)

Author(s):

TAN Zhi¹; 5; JIANG Jinzhong¹; 2; SHI Lei³; 4; LI Yonghua³; 4

(1.Yunnan Earthquake Agency,Kunming 650224,Yunnan,China)(2.Kunming Institute of Prediction,China Earthquake Administration,Kunming 650224,Yunnan,China)(3.Institute of Geophysics,China Earthquake Administration,Beijing 100081,China)(4.Key Laboratory of Earthquake Source Physics,China Earthquake Administration,Beijing 100081,China)(5.Xiaguan Earthquake Monitoring Center Station,Yunnan Earthquake Agency,Dali 671000,Yunnan,China)

关键词:

重震联合反演;  结构耦合;  交叉梯度;  顺序策略

Keywords:

joint inversion of gravity and seismic data;  structural coupling;  cross-gradient;  sequential strategy

分类号:

P315.726

DOI:

10.20015/j.cnki.ISSN1000-0666.2025.0056

摘要:

不同地球物理数据参与的联合反演,可有效减少反演的多解性,为地球内部结构研究提供多参数约束。重力与地震联合反演的开展常基于速度-密度经验关系,但由于地下结构复杂,部分地区岩石物理关系难以准确确定。针对这一问题,提出了一种基于交叉梯度约束的地震走时与重力联合反演方法。该方法将频率域重力反演和近震体波走时反演相结合,采用优化的顺序策略,实现了速度模型与密度模型的相似性耦合。理论模型和青藏高原东北缘实际观测数据试验都显示,联合反演有效提高了密度结构的垂向分辨率。重力数据的引入,更好地约束了地壳浅层(地下5 km以内)P波速度结构,明显改善了体波射线分布稀疏地区的横向分辨能力。

Abstract:

Joint inversion of different geophysical data can effectively reduce the non-uniqueness of results and provide more constraints for studying the Earth's interior structure.Joint inversion of gravity and seismic data is usually carried out based on velocity-density relationships.However,the actual subsurface structures are complex,making it difficult to accurately determine the velocity-density relationships in certain areas.To address this issue,we have developed a joint inversion method for seismic traveltime and gravity based on cross-gradient constraints.This method combines frequency-domain gravity inversion and body wave tomography,and uses an optimized sequential strategy to achieve the structural similarity coupling between the velocity model and the density model.Synthetic tests and real data tests of the Qinghai-Xizang Plateau suggest that this joint inversion can effectively enhance the vertical resolution capabilities of the density structure.Furthermore,incorporating gravity data better constrains the P-wave velocity structure in the shallow Crust(5 km underground)and significantly improves the horizontal resolution capabilities in areas with sparse body wave ray coverage.

参考文献/References:

董治平,张元生.2007.河西走廊中部地区三维速度结构研究[J].地球学报,28(3):270-276.
Dong Z P,Zhang Y S.2007.A study of the 3-D velocity structure in central Hexi corridor[J].Acta Geoscientica Sinica,28(3):270-276.(in Chinese)
甘露,吴庆举,黄清华,等.2022.大地电磁与接收函数结构约束联合反演[J].地球物理学报,65(11):4460-4470.
Gan L,Wu Q J,Huang Q H,et al.2022.Structure constrained joint inversion of magnetotelluric data and receiver function[J].Chinese Journal of Geophysics,65(11):4460-4470.(in Chinese)
李静,万晓云,徐铭. 2023.基于GRACE卫星数据和位错模型的2007年印尼明古鲁M_W8.4地震同震重力梯度信号对比分析[J].地震研究,46(4):521-528.
Li J,Wan X Y,Xu M.2023.Comparative analysis of the coseismic gravity gradient signals of the 2007 Bengkulu,Indonesia M_W8.4 earthquake based on GRACE satellite data and the dislocation model[J].Journal of Seismological Research,46(4):521-528.(in Chinese)
任纪舜,牛宝贵,王军,等.2013.1:500万国际亚洲地质图[J].地球学报,34(1):24-30.
Ren J S,Niu B G,Wang J,et al.2013.1:5 Million international geological map of Asia[J].Acta Geoscientica Sinica,34(1):24-30.(in Chinese)
王长在,吴建平,房立华,等.2011.2009年姚安地震序列定位及震源区三维P波速度结构研究[J].地震学报,33(2):123-133.
Wang C Z,Wu J P,Fang L H,et al.2011.Relocation of aftershocks of the 2009 Yao'an M_S6.0 earthquake and 3-D P-wave velocity structure around its source region[J].Acta Seismologica Sinica,33(2):123-133.(in Chinese)
王新胜,方剑,许厚泽.2013.青藏高原东北缘岩石圈三维密度结构[J].地球物理学报,56(11):3770-3778.
Wang X S,Fang J,Xu H Z.2013.3D density structure of lithosphere beneath northeastern margin of the Tibetan Plateau[J].Chinese Journal of Geophysics,56(11):3770-3778.(in Chinese)
夏思茹,石磊,李永华,等.2021.青藏高原东北缘地壳及上地幔顶部速度结构研究[J].地球物理学报,64(9):3194-3206.
Xia S R,Shi L,Li Y H,et al.2021.Velocity structures of the crust and uppermost mantle beneath the northeastern margin of Tibetan Plateau revealed by double-difference tomography[J].Chinese Journal of Geophysics,64(9):3194-3206.(in Chinese)
杨利荣,岳乐平,王洪亮,等.2016.祁连山及邻区第四纪地层区划与沉积序列[J].中国地质,43(3):1041-1054.
Yang L R,Yue L P,Wang H L,et al.2016.Quaternary stratigraphic realm and sedimentary sequence of the Qilian Mountain and adjacent areas[J].Geology in China,43(3):1041-1054.(in Chinese)
张济东,郑浩,王鸣谦,等.2022.基于交叉梯度函数的三维重震联合反演[J].工程地球物理学报,19(2):258-267.
Zhang J D,Zheng H,Wang M Q,et al.2022.Three-dimensional joint inversion of gravity and seismic data based on cross gradient function[J].Chinese Journal of Engineering Geophysics,19(2):258-267.(in Chinese)
张佩,宋晓东,熊奥林.2019.重震联合反演方法及其应用进展[J].地球物理学进展,34(5):1818-1825.
Zhang P,Song X D,Xiong A L.2019.Advancement of joint inversion of gravity and seismic data and its application[J].Progress in Geophysics,34(5):1818-1825.(in Chinese)
DZ/T0082—1993,区域重力调查规范[S].
DZ/T0082—1993,The standard for regional gravity survey[S].(in Chinese)
Aki K,Lee W H K.1976.Determination of three-dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes:1.A homogeneous initial model[J].Journal of Geophysical Research,81(23):4381-4399.
Astic T,Heagy L J,Oldenburg D W.2020.Petrophysically and geologically guided multi-physics inversion using a dynamic Gaussian mixture model[J].Geophysical Journal International,224(1):40-68.
Brocher T M.2005.Empirical relations between elastic wavespeeds and density in the Earth's crust[J].Bulletin of the Seismological Society of America,95(6):2081-2092.
Cui Y T,Guo L H.2019.A wavenumber-domain iterative approach for rapid 3-D imaging of gravity anomalies and gradients[J].IEEE Access,7:34179-34188.
Deng Q,Zhang P,Ran Y,et al.2003.Basic characteristics of active tectonics of China[J].Science in China(Series D),46(4):356-372.
Fang Y,Wang J,Meng X,et al.2022.An efficient cross-gradient joint inversion algorithm for gravity and magnetic data using a sequential strategy[J].IEEE Transactions on Geoscience and Remote Sensing,60:1-16.
Gallardo L A,Meju M A.2003.Characterization of heterogeneous near-surface materials by joint 2D inversion of DC resistivity and seismic data[J].Geophysical Research Letters,30(13):1658.
Gallardo L A,Meju M A.2004.Joint two-dimensional DC resistivity and seismic travel time inversion with cross-gradients constraints[J].Journal of Geophysical Research,109(B3):B03311.
Gao J,Zhang H.2018.An efficient sequential strategy for realizing cross-gradient joint inversion:method and its application to 2-D cross borehole seismic traveltime and DC resistivity tomography[J].Geophysical Journal International,213(2):1044-1055.
Haber E,Oldenburg D W.1997.Joint inversion:A structural approach[J].Inverse Problems,13(1):63-77.
Hansen P C,O'Leary D P.1993.The use of the L-curve in the regularization of discrete ill-posed problems[J].SIAM Journal on Scientific Computing,14(6):1487-1503.
Huang Y,Moorkamp M,Gao J,et al.2023.Seismogenic structure of the 2014 M6.5 Ludian earthquake from three-dimensional joint inversion of magnetotelluric data and seismic arrival times[J].Journal of Geophysical Research:Solid Earth,128(7):e2022JB026151.
Li G,Cai H,Li C.2019.Alternating joint inversion of controlled-source electromagnetic and seismic data using the joint total variation constraint[J].IEEE Transactions on Geoscience and Remote Sensing,57(8):5914-5922.
Li Y G,Oldenburg D W.1998.3-D inversion of gravity data[J].Geophysics,63(1):109-119.
Maceira M,Ammon C J.2009.Joint inversion of surface wave velocity and gravity observations and its application to central Asian basins shear velocity structure[J].Journal of Geophysical Research:Solid Earth,114(B2):B02314.
Molodtsov D M,Kashtan B,Rostov Y.2011.Joint inversion of seismic and magnetotelluric data with structural constraint based on dot product of image gradients[J].SEG Technical Program Expanded Abstracts,740-744.
Moorkamp M.2022.Deciphering the state of the lower crust and upper mantle with multi-physics inversion[J].Geophysical Research Letters,49(9):e2021GL096336.
Shin Y H,Xu H,Braitenberg C,et al.2007.Moho undulations beneath Tibet from GRACE-integrated gravity data[J].Geophysical Journal International,170(3):971-985.
Sun A,Zhao D,Gao Y,et al.2019.Crustal seismic imaging of Northeast Tibet using first and later phases of earthquakes and explosions[J].Geophysical journal international,217(1):405-421.
Sun A,Zhao D,Guo H.2021.Cracks and fluids in the Northeast Tibetan crust:New insight into seismotectonics[J].Physics of the Earth and Planetary Interiors,311:106634.
Tu X,Zhdanov M S.2021.Joint Gramian inversion of geophysical data with different resolution capabilities:Case study in Yellowstone[J].Geophysical Journal International,226(2):1058-1085.
Um E S,Commer M,Newman G A.2014.A strategy for coupled 3D imaging of large-scale seismic and electromagnetic data sets:Application to subsalt imaging[J].Geophysics,79(3):ID1-ID13.
Vatankhah S,Liu S,Renaut R A,et al.2022.An efficient alternating algorithm for the Lp-norm cross-gradient joint inversion of gravity and magnetic data using the 2-D fast Fourier transform[J].IEEE Transactions on Geoscience and Remote Sensing,60:1-16.
Wadati K.1928.Shallow and deep earthquakes[J].Geophysical Magazine,1:162-202.
Wang S,Liu B,Tian X,et al.2018.Crustal P-wave velocity structure in the northeastern margin of the Qinghai-Tibetan Plateau and insights into crustal deformation[J].Science China Earth Sciences,61(9):1221-1237.
Wang X,Liu J,Dai S,et al.2021.Fast numerical simulation of 2D gravity anomaly based on nonuniform fast Fourier transform in mixed space-wavenumber domain[J].Journal of Applied Geophysics,194:104465.

Wu P,Tan H,Lin C,et al.2020.Joint inversion of two-dimensional magnetotelluric and surface wave dispersion data with cross-gradient constraints[J].Geophysical Journal International,221(2):938-950.
Xiao Q,Zhang J,Wang J,et al.2012.Electrical resistivity structures between the Northern Qilian Mountains and Beishan Block,NW China,and tectonic implications[J].Physics of the Earth and Planetary Interiors,200-201:92-104.
Zhang H J,Thurber C H.2003.Double-difference tomography:The method and its application to the Hayward fault,California[J].Bulletin of the Seismological Society of America,93(5):1875-1889.

备注/Memo

备注/Memo:

收稿日期:2024-09-27.
基金项目:国家自然科学基金项目(42274118,42274080); 云南省地震局地震科技专项(2024ZX03).
第一作者简介:谭 智(1999-),助理工程师,主要从事重力与天然地震联合反演研究.E-mail:2893540431@qq.com.
通信作者简介:石 磊(1984-),研究员,博士,主要从事重力与天然地震联合反演研究.E-mail:shilei@cea-igp.ac.cn.

常用功能

更新日期/Last Update: 2025-05-01