程函方程的接收函数走时反演方法研究青藏高原东南缘莫霍面深度  

Moho depth of the SE Tibetan Plateau inverted by an Eikonal equation-based receiver function traveltime inversion method

作　　者：黄俊杰 刘少林 杨顶辉[2] 李孟洋 宋汉杰 汪文帅 HUANG JunJie;LIU ShaoLin;YANG DingHui;LI MengYang;SONG HanJie;WANG WenShuai(National Institute of Natural Hazards,Ministry of Emergency Management of China,Beijing 100085,China;Department of Mathematical Sciences,Tsinghua University,Beijing 100084,China;School of information Engineering,Ningxia University,Yinchuan 750021,China)

机构地区：[1]应急管理部国家自然灾害防治研究院,北京100085 [2]清华大学数学科学系,北京100084 [3]宁夏大学信息工程学院,银川750021

出　　处：《地球物理学报》2025年第5期1767-1781,共15页Chinese Journal of Geophysics

基　　金：应急管理部国家自然灾害防治研究院基本科研业务专项(ZDJ2024-26);应急管理部国家自然灾害防治研究院创新团队项目(J2223803);国家自然科学基金(42174111,U23A2029)联合资助.

摘　　要：Moho面深度模型可为认识构造演化提供关键信息.接收函数H-κ叠加方法是目前获取Moho面深度最常用的方法之一.然而,复杂的地壳结构会导致H-κ叠加方法产生较大偏差.为了获得高精度Moho面深度成像结果,本文发展了适用于复杂模型的Moho面深度成像方法,即基于程函方程的接收函数走时反演方法.该方法通过数值求解程函方程得到任意非均匀介质中的远震射线和不同震相走时,根据射线信息和不同震相走时可构建tPs-tPp扰动与Moho面深度h扰动之间的对应关系,通过优化方法反演深度h,即可获得任意复杂模型的Moho面深度模型.本文选取青藏高原东南缘Moho面为研究对象,通过和前人研究结果对比,证实了基于程函方程的接收函数走时反演方法成像结果的可靠性.Moho depth model can provide key information for the understanding of tectonic evolution.The receiver function H-κstacking method is currently one of the most commonly used techniques to obtain the Moho depth.However,the Moho depth obtained by the H-κstacking method usually deviate from the true Moho depth when the crust is highly heterogeneous.To obtain high-accuracy Moho depth in complex models,an eikonal equation-based receiver function traveltime inversion method is developed in this study.In this method,the eikonal equation is numerically solved to obtain teleseismic traveltime fields in arbitrary heterogeneous models.Based on the teleseismic traveltime fields,we accurately calculate the teleseismic rays and traveltimes of different phases.The perturbation of tPs-tPp can be related to the perturbation of the Moho depth h according to the information of the teleseismic rays and the traveltimes of teleseismic Pp and Ps phases.Subsequently,the Moho depth is derived by an optimization method.Using the developed eikonal equation-based receiver function traveltime inversion method,we image the Moho depth of the SE Tibetan Plateau.Comparing the results from this study with those from previous researches,we demonstrate the reliability of the imaging results from the developed method.

关 键 词：Moho面深度 程函方程 非均匀模型 接收函数 青藏高原东南缘 

分 类 号：P315[天文地球—地震学]