利用接收函数计算地壳各向异性的可靠性分析及倾斜界面的影响*

使用接收函数研究壳幔速度间断面和速度结构已是常用的技术,但介质各向异性或倾斜的莫霍(Moho)界面,都会造成接收函数波形的复杂性。本文利用远震P波接收函数计算台站下方的地壳各向异性,通过信噪比测试和谐波分析两种方法来验证各向异性结果的可靠性。通过甘肃地震台网的两个台站记录,讨论各向异性和倾斜界面对接收函数的影响,结果显示,台站BYT(白银)下方具有各向异性,而台站WYT(渭源)下方由于可能存在倾斜界面,得到的各向异性结果则有待进一步分析。为了更好地认识倾斜界面对各向异性计算结果的影响,采用合成理论地震图,计算接收函数,然后利用合成接收函数进行各向异性分析。结果表明,计算得到的快波方向不会受到倾斜界面的影响,但是时间延迟会受到影响。     

利用远震接收函数研究安徽及周边地区地壳结构和各向异性

为了研究安徽及周边地区的地壳精细结构及其动力学变形机制,本文利用中国地震台网50个固定台站记录的远震波形数据,采用"水准量"反褶积法提取P波接收函数,采用H-κ叠加技术和莫霍面Ps横波分裂技术,获得了各台站下方的地壳厚度、纵横波速比和地壳各向异性信息.H-κ叠加结果显示:研究区域的地壳厚度和纵横波速比在空间分布上均存在明显的区域特征.华北地块地壳厚度较薄,平均厚度约29 km,且相对平坦,波速比较高,表明地壳中镁铁成分居多.大别造山带地壳较厚,平均为33.1 km,但波速比变化较大.造山带主体波速比约为1.75甚至更低,是长英质-中性岩石特征,与造山带中下地壳拆沉有关,而造山带北部的六安断裂带出现的高波速比与局部岩浆活动产生的基性火山岩有关.扬子地块地壳厚度介于前两者之间,平均为31.5 km,波速比局部特征明显.黄山区域波速比为1.79,反映了该区域中酸性花岗岩特征.地壳各向异性结果显示:整个研究区域的分裂时间为0.2～0.6 s.华北地块横波快波方向总体为NW-SE向,与SKS分裂方向和GPS方向一致,表明壳幔变形方向一致,再结合华北块体地壳厚度较薄和高波速比,推测华北地块在拉张应力下,上地幔物质上涌,侵蚀下地壳,导致地壳变薄和火山岩成分增多.大别造山带区域的快波方向表现为NW-SE向,这与SE向的中下地壳塑性流动或多期的拆离折返构造事件吻合.郯庐断裂带附近的快波方向为NE-SW甚至N-S向,反映了该断裂带在其东西两侧挤压作用下,近南北向拉张运动.     

直达Ps与接收函数Ps联合分析地壳各向异性

本文尝试直接利用地震图记录中的Ps转换波,而非传统的接收函数方法,检测地壳介质的各向异性.虽然远震记录中的P波对Ps转换波会有干扰,但P波与Ps转换波具有不同的偏振特性,最小能量法可以将Ps转换波的各向异性参数检测出来,对时间窗的选择也很强壮.本文利用IC台网（New China Digital Seismograph Network,中国新数字地震台网）中的KMI,ENH,LSA,BJT,XAN,WMQ六个台站接收的宽频远震记录进行了测试,并且对比了接收函数方法的结果.研究表明,利用最小能量法直接测量与利用旋转相关法测量接收函数中Ps转换波的结果基本一致.两种方法测量所得的快波偏振方向与地壳运动方向,最大水平主压应力方向,震源机制解压缩轴方向,最大主压应变率方向以及已有的各向异性研究结果基本一致.     

利用接收函数反演东北地区地壳结构

选取辽宁、吉林、黑龙江、内蒙东部共111个固定台站自2012年1月1日到2015年10月1日,震中距在2000～18000km范围内,震级在6.0以上的542个地震波形记录。采用最大熵反褶积法,经过仔细筛选共获得4683条接收函数,通过HK扫描法得到台站下方地壳厚度及波速比信息,得到了东北莫霍面分布。     

按方位叠加接收函数分析青藏高原东南缘的地壳各向异性

在各向异性地壳中,来自Moho的P-to-S转换波（Pms）的到时不仅取决于入射角和地壳厚度,而且还随地震事件方位角而变化.地处青藏高原东南缘的川滇地区,地壳变形十分强烈.本文利用川滇地区的108个固定台站记录的远震三分量地震波形数据提取台站下方的P波接收函数,并把接收函数被校正到了同一参考震中距处（例如67°）.然后按后方位角10°为间隔将接收函数叠加成一道信号以增强信噪比,并从叠加信号里拾取不同后方位角对应的Pms相的观测到时.在快波极化方向和分裂时间构成的解的平面上,能使观测到时与理论到时之差最小的点即为所求的分裂参数的位置.合成地震图和实际观测数据的实验表明,这个方法不但稳定性较好,而且误差估计也较小.我们从108个台中获得了96个Pms相的分裂参数,结果表明,川滇地区地壳各向异性十分强烈,Pms相分裂时间在0.05 s±0.06 s到1.27 s±0.10 s之间,平均值为0.54 s±0.12 s.地壳各向异性的快波极化方向与地表GPS速度场的差异性表明,印支块体的上下地壳之间是解耦的,而川滇菱形块体北部、松藩—甘孜和四川盆地的上下地壳之间是耦合的.然而,川滇菱形块体南部,地壳变形主要受控于小江断裂和金沙江—红河断裂.

     

地壳分层各向异性介质接收函数及其粒子群反演

地壳不同深度介质的地震各向异性是研究地壳不同深度范围变形方式的重要依据.鉴于地壳介质的复杂性,如何从远震体波接收函数中提取不同深度的各向异性参数仍是一个有待深入研究的课题.在已有研究的基础上,本文利用广义反射-透射系数矩阵方法计算的合成地震图,研究了复杂地壳分层各向异性介质的接收函数随反方位角（back azimuth）变化及不同层位各向异性参数对接收函数波场的影响,为各向异性介质接收函数的解释提供了新的理论依据.通过引入粒子群优化理论,发展了分层各向异性介质接收函数全局反演算法.数值及观测数据的验证结果表明,在各向同性速度模型确定的前提下,我们的方法能够可靠地提取地壳分层各向异性参数;在反演中引入曲波变换去噪技术,对于正确解析不同层位的各向异性参数具有重要价值.     

利用接收函数研究邯郸及周边地区地壳结构

选取邯郸地震台网记录的高质量远震波形数据资料,利用H—Kappa叠加搜索方法得到晋冀鲁豫交界地区(35.0°—38.0°N,112.5°—116.5°E)地震台站下方接收函数、地壳厚度和泊松比分布特征。结果显示:晋冀鲁豫交界地区地壳厚度变化具有明显区域分布特征,以太行山脉为界,以西附近地区地壳厚度约38—41 km,以东的邯郸、邢台、河南地区地壳厚度约31—34 km,而山东地区地壳厚度略有增加,范围约32—35 km,总体呈现自西向东逐渐减薄的规律;泊松比均值约0.251,与以往研究结果相当。     

利用接收函数研究青藏高原东部地壳结构

青海玉树Ms7.1地震之后,根据地震断裂及余震活动的分布情况,于2010年5～11月在玉树地震断裂带周边地区布设由26套宽频带流动地震仪组成的临时地震台网,组成一个经过玉树地区的南北剖面和分布在玉树周边地区的密集台阵.其中南北剖面台站分布由北往南依次穿过昆仑地块、巴彦喀拉地块和羌塘地块,由15个流动地震观测台和固定台站MAD台组成;密集台阵由由分布在玉树周边的14个台站组成.     

接收函数方法估计Moho倾斜地区的地壳速度比

造山带地区的复合力系作用往往使Moho界面发生变形,局部表现为倾斜状态.为了得到这些区域精确地壳速度比结构,本文基于H-κ方法发展了H-κ-θ方法.该方法不仅考虑了倾斜Moho层的响应,同时利用径向和切向接收函数信息,增加了对扫描的约束.利用该方法对青藏高原东南缘地壳厚度和速度比结构进行研究,结果表明：研究区内地壳明显存在不均匀性,松潘—甘孜地体平均地壳厚度约为60 km,四川盆地西缘约为47 km,扬子地台约为43 km,三江块体和扬子地台东南缘已接近正常地壳厚度;松潘—甘孜地体与扬子地台相邻部位地壳平均地震波速度比（V_p/V_s）普遍偏高,且四川盆地西侧发现一绕盆地边缘的弧形高V_p/V_s异常区（>1.88）,我们推测该异常可能由青藏高原向东逃逸的地壳流体受到高强度的四川盆地阻挡,在其西侧堆积所致.     

利用远震接收函数探测四川盆地及周边地区的地壳结构

青藏高原东南缘的龙门山断裂两侧具有陡峭的地形特征,在约50~100 km的水平距离内,地形高程从2000 m增加到4000 m,该区强烈的壳幔变形特征及地球动力学模式一直是研究的热点问题.本文从四川地区49个固定台站记录的远震资料提取了P波接收函数,获得了四川盆地及周边的地壳厚度和泊松比,并以此构建反演的初始模型.在线性反演的基础上,引入了分别拟合低频和高频接收函数的两步反演技术,用以反演台站下方的地壳S波速度结构.数字试验表明,该方法可以有效抑制接收函数反演的不唯一性,为了得到最优解,最后用Bootstrap重采样技术估计解的不确定性.结果表明,四川盆地的地壳厚度在40~46 km,松潘-甘孜块体北部的地壳厚度为46~52 km,而南部增厚到50~60 km.从四川盆地向西跨过龙门山断裂,地壳厚度增加了10~15 km.在四川盆地及周边地区,地壳泊松比在0.26~0.32之间,呈块体分布特征,高泊松比（0.28~0.32）主要沿龙门山断裂以及安宁河-小江断裂分布.地壳S波速度结构表明,来自青藏高原中部的中下地壳低速层可能受到了坚硬的四川盆地阻挡,改变原来的运动方向并沿龙门山断裂展布,由于低速层的囤积导致该区地形陡峭和下地壳增厚.

     

利用远震接收函数分析三峡库区重庆段地壳厚度变化

针对三峡库区重庆段地震监测台网和重庆区域台网2007年到2010年记录到的远震数据,利用接收函数方法,得到该区域的地壳厚度值及泊松比。结果显示27个地震台下方的地壳平均厚度为44.6km。最厚处为城口(CHK),达50.9km;最薄处为荣昌(ROC),仅38.9km,两者相差12km。研究区的地壳泊松比σ值约为0.27,三峡重庆库区东部奉节荆竹(JIZ)地壳泊松比高达0.302,武隆(WUL)地区的地壳泊松比较低,最低为0.228,两者相差0.074;研究区内泊松比的最大值偏离标准泊松体比值(0.25)20.8%以上,其最小值偏离标准泊松体比值-8.80%。布格重力异常反映了地下的密度与地壳厚度变化。文中得出的三峡库区重庆段地壳厚度变化与该区布格重力异常分析结果一致,重庆境内最小的负布格重力异常在地壳厚度最厚的城口—巫溪地区(厚49.7～50.9km),地壳最薄的荣昌—合川地区(厚38.9～41km)对应着渝西最大的负布格重力异常。     

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

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

利用接收函数方法研究四川地区地壳结构

采用接收函数反演和共转换点(CCP)偏移叠加成像方法, 利用四川数字地震台网宽频带的52个区域固定地震台站和布设的两条52个宽频带流动地震观测台站的远震地震波形数据资料, 对四川地区地壳结构进行研究。 结果表明, 四川地区的Moho面深度在青藏高原和四川盆地差异明显, 在川西高原地区地壳厚度为52~68 km, 在川滇地块地壳厚度为50~60 km, 在中地壳内存在不连续的低速层分布; 而在四川盆地地壳厚度为38~45 km, 地壳内没有低速层存在。 Moho面深度从川西高原的60多公里至四川盆地的约40 km, 在二者的交界处龙门山断裂带下面, 存在厚度约30 km左右宽的下降过渡带, 说明其下的Moho面可能受断层影响, 结构比较复杂; 在高原地区的上地壳界面和下地壳上界面比四川盆地的相应界面深; 高原地区在中地壳的上部有不连续的低速层分布, 在松潘—甘孜地块的上地壳下部存在向南东运动的脆性推覆体, 在羌塘—理塘地块的上地壳下部存在向南东和南运动的脆性物质流动。     

利用远震接收函数研究辽宁地区的地壳厚度及泊松比

文中利用辽宁省数字地震遥测台网15个台站记录的远震P波波形资料,用频率域反褶积方法提取接收函数,由H-Kappa叠加方法得到了各台站下方的地壳厚度和泊松比。研究结果表明,辽宁地区的地壳泊松比在0.25～0.29之间,地壳厚度介于31～36km之间,东部褶隆带的地壳厚度从北向南由31km增至35km。松辽构造盆地地壳厚度变化不大,平均厚度为31km。辽西褶隆带与东部褶隆带的地壳厚度均比哈尔滨构造盆地厚约2～4km     

利用接收函数研究滇西北地区地壳厚度及泊松比特征

计算了滇西北地区的16个固定地震台和喜马拉雅台阵项目的103个流动台站(共119个地震台站)记录的238个5.8级以上的远震事件,从中挑选了5558个信噪比高、震相清晰的接收函数,采用人工读取震相到时的方法,获得了各个台站下方的地壳厚度和泊松比.结果显示:滇西北地区地壳厚度等值线呈ES向舌状突出,地壳厚度和泊松比横向变...     

利用接收函数研究六盘山地区地壳上地幔结构特征

利用六盘山地区宽频带流动地震台阵的远震体波记录,采用接收函数方法研究台阵下方的地壳上地幔结构,并采用接收函数振幅加权叠加方法对这一地区平均地壳厚度和泊松比进行估算.研究结果显示,青藏高原东北缘和鄂尔多斯的接触过渡带接收函数震相复杂,地壳变形强烈,青藏高原东北缘地壳平均厚度约为51.5km,六盘山下方地壳厚度在53.5km,鄂尔多斯西南缘地壳平均厚度约为50km,整个莫霍面呈下凹状.泊松比计算结果显示,六盘山东侧和西侧地壳泊松比值在正常范围内(0.25～0.26),六盘山下方的地壳泊松比值偏高(0.27～0.29),推测与地壳中存在部分熔融.泊松比值的横向变化,指示测线范围地壳物质组成和力学性质存在横向差异,反映在欧亚板块与印度板块碰撞作用影响下青藏高原下地壳物质存在向北东方向的流动.     

利用接收函数研究延边地震台下方地壳结构

收集整理2007年以来延边地震台记录的113个远震数字波形资料,采用远震接收函数反演延边地震台下方地壳结构,运用H-Kappa叠加方法,计算得到台站下方地壳厚度和泊松比.采用全球平均地壳模型作为初始模型,反演台站下方0-100 km的S波速结构.反演结果表明,延边地震台下方地壳厚度为30.8 km,波速比为1.84,泊松比较高,为0.29.在台站下方15-20 km及25-30 km处存在低速层.     

利用远震接收函数探测襄樊台下方地壳厚度及泊松比

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

利用接收函数方位变化研究太行山区地壳各向异性

探讨了利用接收函数方位变化研究各向异性复杂地区速度结构的可能性. 详细介绍了具有任意对称轴方向的横向各向同性系统中介质响应的计算过程,并将利用接收函数方位变化研究地壳各向异性的方法应用于华北太行山地区. 结果表明,研究台站下方地壳具有明显的各向异性特征,其中地表盖层表现为8％——15％的慢轴型各向异性;深部地壳各向异性略小,表现为3％——5％的快轴型各向异性特征. 不同台站获得的各向异性对称轴的方向及深度分布不具有明显的一致性. 区域构造环境的变化可能是导致地壳各向异性存在的原因.     

利用SV分量接收函数反演地壳横波速度结构

详细讨论了远震体波SV分量接收函数的特点及其在反演地壳S波速度结构中的优势．与径向接收函数类似,SV分量接收函数可通过对远震体波的SV分量直接反褶积P分量获得．研究分析表明：与径向接收函数相比,SV分量接收函数的振幅随震中距的变化更加稳定,波形简单且突出了对结构最敏感的PS转换波信息．理论数值实验显示：在反演地壳S波速度结构时,SV分量接收函数比径向接收函数具有更好的收敛性．作为实例,利用SV分量接收函数反演方法反演了海拉尔台下的S波速度结构．