基于接收函数与背景噪声联合反演的研究与应用

本文回顾了接收函数线性反演及接收函数与频散曲线联合反演的理论.为了检验两种方法在实测数据中的应用效果,我们选择云南区域台网位于四个不同次级块体的四个台站为例:首先,计算了四个台站的远震P波接收函数;然后,分别运用接收函数线性反演和接收函数与背景噪声数据联合反演两种方法分别得到台站下方的S波速度结构;另外,运用H-k扫描方法验证两种反演方法获取的莫霍面深度的可靠性.在反演过程中,当我们提供具有深一些莫霍面的初始模型时,联合反演获取的莫霍面深度与H-k扫描结果具有很好的一致性,同时联合反演也能有效地捕捉低速层,对反演的S波速度也有较好的约束.由此推知,加入背景噪声相速度数据后的联合反演对初始模型的依赖性低于单独用接收函数反演的结果.     

基于P波接收函数资料的华南大陆东部地壳结构研究

华南大陆是由扬子块体和华夏块体拼合而成,并经历了多期次的后期改造,研究块体间的深部结构差异对于理解华南大陆的构造演化历史具有重要意义.本文基于布设在华南大陆东部的流动宽频带地震台阵资料,计算了远震事件的P波径向接收函数,并利用Pms谐波分析方法获得了地壳各向异性分布.结果显示,研究区地壳各向异性具有明显的横向分段特征.JSF（江绍断裂带）附近与断裂带走向一致的地壳各向异性快波方向指示断裂带可能切穿整个地壳.在大别造山带、扬子块体和华夏块体南部,地壳各向异性快波方向与区域主压应力方向近垂直,指示该区域地壳变形可能受菲律宾海板块俯冲和太平洋板片俯冲后撤的远程效应的影响.此外,本文利用新的H-κ-c叠加方法获得了测线下方地壳厚度及平均波速比.结果显示,扬子块体的地壳厚度薄于两侧的华夏块体和大别造山带,速度比也较低（约1.7）,指示地壳以长英质和酸性岩石为主.Airy地壳均衡理论预测的地壳厚度指示29°N以南的地壳密度要低于以北的区域.以南的低密度及低波速比可能与菲律宾海板块北西向俯冲和太平洋板块西向俯冲后撤导致的地壳物质垂向分异拆沉有关.

     

用接收函数研究川滇地区国家地震台下地壳厚度及波速比

本文利用远震接收函数的方法,对川滇地区的昆明、腾冲、成都和攀枝花等4个国家地震台的台基下方不同方向的莫霍面深度及波速比进行了研究和分析。结果表明:昆明地震台台基下方的莫霍面深度基本在50km左右,波速比为1.62～1.69,地壳厚度和波速比不因方向不同而发生明显的变化;腾冲地震台台基下方的地壳厚度有着比较明显的方向性,东北方向厚为40.7km,东南方向为49.7km,两个方向的波速比相差也很大,差值达到0.2;成都地震台台基下方莫霍面的深度在40km左右,但是东北和西南方向要加深8km,两个方向波速比相差0.13;攀枝花地震台台基下方的地壳厚度比较稳定,厚度在60km左右,波速比变化也不明显。     

接收函数、面波与重力联合约束地壳厚度与波速比

提出了一种利用接收函数、面波频散和重力数据联合约束地壳厚度、v_P/v_S和平均P波速度的改进方法,并基于两种地壳模型对改进后的方法进行了验证。结果显示,改进后的方法不仅可以提高地壳厚度和波速比的估计精度,还能更准确地估计地壳平均P波速度。在此基础上,将该方法应用于华南地区两个固定台站的地壳结构分析,相关结果也证实了该方法在确定地壳结构中的可行性。     

接收函数揭示的兴蒙造山带西南部的地壳结构特征

针对构造演化历史甚为复杂的兴蒙造山带的地壳结构探测较为有限的问题,本研究利用2016-2018年布设在兴蒙造山带西南部的第3期NECsaids地震台阵和固定地震台及NECESSArray的流动台数据,采用时间域最大熵谱反褶积方法提取到研究区的24027条高质量P波接收函数,进而采用H-κ和共转换点（CCP）叠加方法以及分层剥离转换震相分析地壳各向异性的方法来辨识研究区的地壳结构特征.分析表明：研究区地壳各向异性整体呈近NW向和近EW向,与区域速度场和板块绝对运动方向一致,表明现今地壳结构主要受控于太平洋板块构造域;大兴安岭-太行山重力梯级带是明显的地壳结构差异过渡区,其西侧地壳厚度明显高于东侧且呈降低趋势,而平均波速比则呈现东西两侧盆地较高的分布特征;局部地区的高波速比和不同的各向异性特征则显示受到新生代火山和断裂活动的改造作用,二连盆地的厚地壳和复杂的各向异性特征则展现出古亚洲洋闭合及蒙古-鄂霍茨克和古太平洋等多期构造作用的可能影响.

     

利用接收函数方法确定青藏高原东北缘近地表S波速度

利用襄樊台2015—2018年记录的震中距30°—90°,震级5.5级以上（含5.5级）的高信噪比远震波形事件,通过时间域反褶积方法提取台站接收函数,并用H-k扫描叠加方法得到台站下方地壳厚度和波速比。结果表明,襄樊台下方地壳厚度为35 km,略高于大陆地壳平均厚度33 km,处于西部山区厚地壳与中东部薄地壳的过渡地带。泊松比偏低,存在因区域内地壳应力积累发生微小地震的可能性。     

用接收函数方法研究中国境内地壳结构

利用中国数字地震台网30个台站的高质量宽频带远震数据,采用H-k叠加搜索法对中国境内的地壳结构进行研究,获得了研究区内的地壳厚度和v_P/v_S分布特征.结果表明, 中国境内的v_P/v_S值介于1.6—1.9之间,地壳厚度变化剧烈,在29—81 km之间.100°—110°E之间存在一个地壳厚度陡变带, 将中国分为东西两个部分.东部地壳厚度相对均匀,为31—36 km, 西部地区地壳厚度相对较厚且变化较大,中部地区地壳厚度为34—49 km.总的看来,青藏高原地区地壳最厚,可达81 km;天山、准噶尔盆地和内蒙古地区地壳厚度次之;华南地区地壳最薄.另外,中国大陆地壳平均波速比为1.738(σ=0.253),比全球大陆平均波速比1.78(σ=0.269)低.较低的波速比可能暗示中国境内地壳低速层的存在或者铁镁质成分的缺失.      

基于多台站的接收函数和重力联合反演确定莫霍面起伏和地壳平均波速比

地壳厚度和波速比是研究地壳结构和组分的两个重要参数,可为区域构造研究提供重要约束。接收函数被广泛地用于确定地壳厚度和波速比,例如H-κ方法或H-κ-c方法,但是该方法只能确定台站下方的地壳厚度和速度比,当地震台站分布稀疏时,很难约束台站间的横向不均匀性。另一方面,重力数据也可用于莫霍面的起伏变化研究,它在横向上覆盖很好,有较高的分辨率,但在纵向上分辨率相对较低。为此,本研究提出了一种联合反演算法求解莫霍面深度和地壳波速比参数。联合反演算法综合考虑了接收函数在纵向上的较高分辨率和重力数据在横向上的较高分辨率,同时拟合区域内所有台站上的接收函数和区域重力数据。模型测试表明联合反演算法较单一的接收函数反演更精确,特别是对于地壳厚度的确定。     

用面波频散与远震接收函数联合反演对朝鲜半岛三维地壳成像

结合Cho等(2006b)的高频面波层析成像结果与朝鲜半岛80个点位上的远震P波接收函数,对朝鲜半岛南部的剪切波速度的三维变化进行了祥细的研究。这些接收函数是根据韩国气象局和韩国地球科学与矿产资源学院台网提供的高增益加速度、短周期和宽带数字数据集求出的。垂直剖面追踪了莫霍面深度的水平变化、近地表的低速度变化以及自地表速度到中地壳速度过渡带的厚度变化。根据获得的地壳结构提出了关于朝鲜半岛地壳演变的新观点。     

中国东北地区高分辨率地壳结构:远震接收函数

利用分布在东北地区的国家地震局台网、NECESSArray台网、吉林大学在长白山及其周边地区布设的26个临时台站总计259个台站接收到的16,070条高质量的P波接收函数,采用H-k和CCP(Common Conversion Point,共转换点)叠加成像方法,获得该区高分辨率的地壳结构.观测结果显示,东北地区莫霍界面深度和地表高程总体呈镜像关系;西部大兴安岭—太行山重力梯级带附近存在莫霍界面深度陡变带;中部的松辽盆地地区受晚中生代的地壳拉伸作用影响,地壳厚度较薄,北部的小兴安岭地区和南部的华北北缘造山带可能同样受拉伸运动影响,具有较小的地壳厚度;松辽盆地莫霍界面深度由西向东逐渐减小,推测这与太平洋板块俯冲作用有关;东部地区莫霍界面结构比较复杂,依兰—伊通断裂与敦化—密山断裂之间出现复杂震相,可能与该区存在地幔物质的底侵作用有关;长白山火山地区地壳厚度较大,对应较高的波速比,推测在该区地壳内存在岩浆囊.     

利用接收函数方法研究西秦岭构造带及其邻区地壳结构

根据西秦岭构造带及其周边地区117个宽频带地震台站的高质量波形数据,利用远震P波接收函数的H-k叠加方法,求得地壳厚度和平均波速比.通过分析地壳厚度、波速比及其关系和接收函数CCP叠加剖面,研究了该区域的地壳结构特征.结果表明,研究区域内地壳结构差异大,呈过渡带特征.地壳厚度总体上呈北北西向分布,自西南向东北逐渐减小.羌塘块体地壳厚度为72km,渭河盆地附近为39km.西秦岭构造带的地壳厚度为42—56km,南北向莫霍界面平坦.研究区域P波与S波波速比平均为1.74,其中西秦岭构造带平均为1.72.较低的波速比主要分布在西秦岭构造带、祁连山块体、松潘—甘孜地块北部以及香山—天景山断裂区域,这可能是由于含长英质酸性岩组分的上地壳叠置增厚而导致的.该区域缺少超高波速比,表明这一区域发生岩浆底侵或上地壳熔融的可能性很小.综合分析表明,西秦岭构造带及邻区的地壳结构主要是由于青藏高原隆升并在向东北向扩张中受到周边块体的阻挡而引起的地壳构造变形所致.西秦岭构造带的莫霍界面变化和波速比分布与该构造带经历碰撞地壳增厚后的伸展走滑运动有关.     

接收函数方法的研究综述

从远震体波波形数据中提取的P波接收函数和S波接收函数已经成为研究台站下方地壳上地幔速度间断面的最有效的方法之一.P波接收函数已经被广泛应用于获取地壳内部S波速度结构、地壳厚度及物质成分组成、地幔过渡带的厚度变化以及岩石圈地幔的间断面等,而S波接收函数是P波接收函数的一个很好的补充,因为在Moho和地幔过渡带之间的深度范围内,Sp转换波比来自浅部间断面的多次波到达早,而在P波接收函数中,相同深度范围的间断面的Ps转换波往往被来自浅部间断面的多次波干扰或者淹没,因此S波接收函数是目前获取岩石圈地幔深度范围内速度间断面结构(如Moho和LAB)的比较有效的方法,比面波观测具有更高的分辨率.本文详细阐述了P波接收函数和S波接收函数的方法原理以及在地壳上地幔速度间断面的研究中所采用的研究思路.     

利用接收函数方法研究西秦岭构造带及其邻区地壳结构

根据西秦岭构造带及其周边地区117个宽频带地震台站的高质量波形数据, 利用远震P波接收函数的H-k叠加方法, 求得地壳厚度和平均波速比. 通过分析地壳厚度、 波速比及其关系和接收函数CCP叠加剖面, 研究了该区域的地壳结构特征. 结果表明, 研究区域内地壳结构差异大, 呈过渡带特征. 地壳厚度总体上呈北北西向分布, 自西南向东北逐渐减小. 羌塘块体地壳厚度为72 km, 渭河盆地附近为39 km. 西秦岭构造带的地壳厚度为42—56 km, 南北向莫霍界面平坦. 研究区域P波与S波波速比平均为1.74, 其中西秦岭构造带平均为1.72. 较低的波速比主要分布在西秦岭构造带、 祁连山块体、 松潘—甘孜地块北部以及香山—天景山断裂区域, 这可能是由于含长英质酸性岩组分的上地壳叠置增厚而导致的. 该区域缺少超高波速比, 表明这一区域发生岩浆底侵或上地壳熔融的可能性很小. 综合分析表明, 西秦岭构造带及邻区的地壳结构主要是由于青藏高原隆升并在向东北向扩张中受到周边块体的阻挡而引起的地壳构造变形所致. 西秦岭构造带的莫霍界面变化和波速比分布与该构造带经历碰撞地壳增厚后的伸展走滑运动有关.     

利用接收函数方法研究胶东半岛地区地壳结构

<正>1研究背景胶东半岛地处山东省东部,位于黄海、渤海之间,隶属华北块体并位于华北克拉通的东部边缘,华北地块东缘与苏鲁造山带山东段在此交会,在大地构造划分上属于华北克拉通与扬子克拉通的结合部位（晁洪太等,2001;周瑶琪等,2015）。华北克拉通破坏是中国东部地质演化的重要事件,胶东半岛位于华北克拉通遭遇破坏的东部边缘,     

接收函数方法及其新的进展

远震P波波形数据中包含了大量地震台站下方地壳和上地幔速度间断面所产生的PS转换波及其多次反射波的信息，由此提取的接收函数是了解地壳上地幔速度细结构的重要步骤手段之一．最近几年，接收函数和面波联合反演方法获得了较大的成功，两种分别对两种波形进行拟合，对反演的速度结构提供了有效的约束．地震勘探中的一些成熟技术被引进接收函数的数据处理，使其可用于地壳和上地幔主要速度界面的侧向变化研究．为增强接收函数的信噪比，将不同事件提取的接收函数进行分类，按方位角进行叠加，可以反应一些速度界面的横向变化．接收函数方法可用于PS转换震相的剪切波分裂研究．泊松比是推测地球内部物质构成的有效参数之一．接收函数方法分离出的转换波为获取泊松比提供了一条行之有效的便捷途径，单一的地震台远震记录中转换波可以估计点位下面地壳Vp/Vs值．另一方面，一些最新的反演技术被引入接收函数的反演．例如，格子收索方法力图解决接收函数反演中的不稳定性和非唯一性问题，格子搜索设计比较容易地结合先验约束并保证搜索是全空间的，避免了使用任何初始模型．该方法能保证在格子的间隔和参数的限度内获得全局最小，该方法被用于中东、南美及临近地区的深部结构．相邻算法是基于计算几何的概念而构建的一种非线性反演方法．该方法避免了此前一些方法(遗传算法、模拟退火)的一些缺陷，如大量样本的舍弃、过多参数的引进等，该算法对我国的五大连池和腾冲火山区的深部结构的反演，以及滇西地区的深部动力学研究获得了良好的效果．     

利用接收函数方法研究蒙古中南部地区地壳结构

2011年8月至2013年7月中国地震局地球物理研究所与蒙古科学院天文与地球物理研究中心在蒙古中南部区域布设了宽频带流动地震台阵,这为开展远东地区深部结构的精细探测提供了有利的数据基础.利用台阵记录的远震地震事件,采用P波接收函数的H-κ叠加分析和共转换点（CCP）叠加方法获得了台站下方的地壳厚度及平均波速比.结果显示研究区的地壳厚度介于39 km至45 km之间.整体上Moho面埋深从西北往东南方向逐渐变浅.在蒙古主线性构造两侧地壳厚度呈现区域性变化特征,东南部地区地壳厚度较薄,约为39 km,而西北部地区地壳较厚,达45 km,为此推测蒙古主构造线可能是地壳的一个陡变带.此外,研究地区地壳的平均波速比值（VP/VS）在1.70到1.79之间,均值为1.75,低于全球大陆的平均值1.78,这可能暗示着该区其地壳是缺少铁镁质的.研究还发现测线的西北与东南地区其地壳波速比值较高,推测是古生代铁镁质地壳的残留或是新生代岩浆底侵的反映.     

接收函数方法及研究进展

远震P波波形数据中包含了大量在台站下方地壳上地幔速度间断面所产生的P-S转换波及其多次反射波的信息，是研究台站下方局部区域S速度分布理想的震相，由此产生的接收函数方法是反演台站下方S波速度结构的有效手段。接收函数方法可以通过波形反演拟合接收函数的径向分量，对观测台站下方地球介质的S波速度结构进行估计，也可以通过偏移叠加获得的接收函数道集（地震剖面图）追踪速度间断面。这种方法避免了对天然地震震源及其附近结构混响效应等复杂因素的影响，对S波速度的垂向分布敏感，垂向分辨率高。由于宽频带流动地震台阵的发展，用此方法还可获得研究区域速度结构的横向变化，横向分辨能力主要取决于台站的间距。本文回顾20年来接收函数研究的进展，探讨了方法研究的发展趋势，介绍了对地壳-上地幔结构的部分研究结果。     

接收函数方法研究进展

远震P波波形数据中包含了大量的地震台站下方地壳和上地慢速度间断面所产生的PS转换波及其多次反射波的信息，由此提取的接收函数是了解地壳上地慢速度精细结构的重要步骤之一。本介绍了目前在提高接收函数结果的稳定性和精度方面的研究进展，包括石油勘探中一些成熟的地震反射处理方法逐步渗透到接收函数的研究领域。这些处理方法以前所未有的分辨率展示了地壳上地幔结构的横向非均匀性。     

接收函数的H-k方法及其应用

接收函数的H-k方法是研究地球内部结构的重要方法,由于接收函数不依赖于震源参数,可以更为精确地确定观测台站下方的地壳结构;同时H-k方法简单明了,可以利用深度-波速比扫描同时获取台站下方的地壳厚度及地壳内部的波速比或泊松比。本文利用H-k方法,对桂东北地区的宽频带地震流动台网观测数据的地震数据进行应用,结果表明该区的地壳厚度变化大（28~34 km）,而且地壳内部的波速比（1.801~1.897）或泊松比（0.277~0.308）较高,与桂东北地区的基性、超基性岩大片出露的事实相吻合。      

东北地区至华北北缘地壳结构的区域差异: 地壳厚度与波速比的联合约束

本文通过对分布相对均匀的127个固定台站下方接收函数的H-κ叠加分析,并结合前人对97个线性密集流动台站的研究结果,获得了东北地区和华北克拉通北缘地壳厚度(H)与平均波速比(κ).结果表明研究区域地壳总体较薄,波速比变化复杂,地壳密度横向变化大,暗示着地壳在中—新生代经历了显著的不均匀破坏与改造.东北和华北北缘都存在明显的东西向差异.东北地区西侧兴蒙造山带地壳大致随着海拔增高逐渐增厚,H和κ分别主要在31~39 km和1.71~1.83之间变化,平均值分别为~35 km和~1.77;东侧吉黑褶皱带地壳厚度与海拔不成镜像关系,H和κ集中在28~37 km和1.72~1.89范围, 平均值分别为~33 km和~1.79.华北北缘西侧燕山带地壳由东往西逐渐增厚,H和κ主要在28~40 km和1.70~1.91范围内变化, 平均值分别为~34 km和~1.79.东侧辽东台隆地壳表现为中间厚四周薄,H和κ集中在29~35 km和1.71~1.83范围, 平均值分别为~32 km和~1.77.东北地区吉黑褶皱带相对薄的H和变化范围大的κ表明,该区域可能由于其自身的地壳结构复杂性和紧邻太平洋板块前缘从而在中新生代遭受到了与太平洋板块俯冲相关的更为强烈的地壳减薄与改造.华北北缘燕山带H和κ复杂的变化特征表明,该地区可能受到中亚造山带增生和太平洋板块俯冲的共同影响,从而发生了更为复杂的地壳改造变形.