松辽盆地沉积层结构的短周期地震背景噪声成像研究

使用位于松辽盆地内部的NECESSArray台阵连续两年背景噪声数据,通过波形互相关和多重滤波方法提取到2~14 s较短周期的Rayleigh波群速度和相速度频散曲线,基于快速行进（FMM）面波成像方法得到群速度和相速度成像结果,并采用最小二乘迭代线性方法反演获得了松辽盆地深至12 km的三维S波速度结构.本文成像结果显示：松辽盆地内部S波速度分布的横向不均匀性与该区域的构造单元呈现出良好的空间对应关系.从地表至下方的6 km深度,盆地北部比南部表现出更加强烈的低速异常,这一特征可能与盆地南北的沉积构造差异有关.中央坳陷区低速异常的边界与嫩江断裂走向相互平行,表明盆地基底断裂对盆地形成演化具有一定的控制作用.在垂直速度结构剖面中,2.9 km·s^-1的S波速度等值线与地震反射剖面显示的盆地基底深度大致对应.基于S波速度模型和盆地基底速度（2.9 km·s^-1）,我们获得精细的松辽盆地沉积层厚度模型,结果表明松辽盆地的沉积层厚度分布呈现出中间厚、四周薄的特征,中央坳陷区的沉积层厚度范围大约在3~6 km.     

基于深度学习的接收函数横波速度预测

远震接收函数中包含了大量台站下方速度间断面所产生的Ps转换波及其多次反射波的信息,被广泛应用于反演精细的地壳及上地幔速度结构.但当地壳结构复杂如存在沉积层或高速层时,转换波和多次波震相的到时和振幅发生变化,从而导致接收函数反演存在强非唯一性.深度学习作为一种高效的特征提取方法,能够建立接收函数与横波速度在空间上的映射关...     

面波频散与接收函数联合反演南北地震带北段壳幔速度结构

通过对南北地震带北段区域所布设的676个流动地震台站观测资料进行处理,联合反演面波频散与接收函数数据,获得了研究区内地壳厚度、沉积层厚度的分布情况以及地壳上地幔高分辨率S波速度结构成像结果.反演结果显示研究区地壳厚度从青藏高原东北缘向外总体逐渐变薄,秦岭造山带地壳厚度较同属青藏高原东北缘的北祁连块体明显减薄;鄂尔多斯盆地及河套盆地分布有非常厚的沉积层,阿拉善块体部分区域也有一定沉积层分布,沉积层与研究区内盆地位置较为一致;松潘—甘孜块体、北祁连造山带等青藏高原东北缘总体表现为S波低速异常;在中下地壳,松潘—甘孜块体下方的低速体比北祁连造山带下方的低速体S波速度值更小、分布深度更浅,更有可能对应于部分熔融的地壳;鄂尔多斯盆地在中下地壳以及上地幔内有着较大范围的高速异常一直延伸到120 km以下,而河套盆地地幔只在80 km以上部分有着高速异常的分布,此深度可能代表了河套盆地的岩石圈厚度,来自深部地幔的热物质上涌造成了该区域的岩石圈减薄;阿拉善块体在地壳和上地幔都表现出高低速共存的分布特征,暗示阿拉善块体西部岩石圈可能受青藏高原东北缘的挤压作用发生改造.     

接收函数和重力联合反演中国东北地区高精度莫霍面和波速比分布及地质意义

本文通过接收函数和重力数据联合反演获取了中国东北地区高精度的莫霍面和平均地壳波速比模型,新的联合反演方法综合考虑了接收函数在纵向上和重力数据在横向上的优势,在某种程度上可以降低单一接收函数数据易受盆地沉积层、莫霍面倾斜、介质不均匀分布等因素的干扰.在接收函数数据的预处理中,通过远震接收函数波形重建方法构建规则化的虚拟台站,提升整体的空间分辨率.在重力数据的滤波中,分别对莫霍面上方的浅部区域和岩石圈以下的深部背景场采用速度结构导向和向上延拓的滤波手段.新获得的莫霍面分布显示,海拉尔盆地、二连盆地、三江盆地、渤海湾盆地的莫霍面深度略浅于周围地区,而大兴安岭、小兴安岭、张广才岭的莫霍面略深.与前人结果相比,本研究得到的精细莫霍面模型与地形起伏和盆地边界有更好的对应关系,且在南北重力梯度线两侧的差异更明显.松辽盆地东部的莫霍面相对较浅,且存在高波速比异常,可能是由于松辽盆地区域的横向拉张引起的地幔热物质侵入地壳引起的.长白山火山地区下方地壳存在高波速比,意味着该处可能存在由于地幔热物质上涌引起的地壳部分熔融.     

利用远震接收函数研究黄海周边地区的地壳上地幔结构

利用宽频带流动台站(YSBSN)记录的远震波形数据和远震接收函数方法,反演了黄海东、西两侧地壳上地幔的S波速度结构.结果表明,莫霍面深度在30~38 km之间变化,位于中方一侧的JNN台下方地壳厚度最大,可以归因于华北板块和扬子板块的碰撞;韩方一侧的地壳厚度自北向南逐渐变厚,但仍然难以厘定朝鲜半岛南部潜在碰撞带的位置,这些问题的解决需要更大范围的流动台站观测.由于部分台站位于巨厚的沉积层和多孔的火山岩之上,与浅部构造的相关性使得接收函数表现出较大振幅的混响,从而影响了来自深部结构的转换震相.     

松辽盆地地壳精细结构研究

基于东北地区已有的宽频带流动台阵远震数据,利用波场延拓和分解的H-β网格搜索法,对松辽盆地的沉积层及地壳结构进行了深入分析。结果显示:松辽盆地的沉积层厚度为0.2—2.5 km,整体呈现中央坳陷区厚、边缘薄且西南地区最薄的分布特征;研究区地壳较薄,厚度介于24—34 km之间,其横向变化特征与沉积层厚度分布具有一定的对应性。依据沉积层和地壳的厚度计算了地壳伸展系数,其平均值接近于以往接收函数研究估测的岩石圈伸展因子。因此,本文推测松辽盆地在伸展构造过程中,其地壳和岩石圈的减薄以纯剪切模式为主。此外,松辽盆地具有较高的地壳平均波速比v_P/v_S,暗示盆地下方岩石圈地幔的减薄过程中可能存在岩浆的底侵作用。      

基于贝叶斯理论的接收函数与环境噪声联合反演

基于Bayes反演理论(Tarantola,1987,2005),在接收函数非线性复谱比反演方法基础上(刘启元等,1996),本文讨论了接收函数与地震环境噪声Rayleigh波相速度频散的联合反演.本文采用修正后的快速广义反射/透射系数方法(Pei et al., 2008,2009) 计算Rayleigh波相速度频散, 并引入地壳泊松比的全局性搜索.数值检验表明:(1)接收函数与环境噪声的联合反演能够有效地解决反演结果对初始模型依赖的问题,即使对地壳速度结构仅有非常粗略的初始估计(例如,垂向均匀模型),本文方法仍能给出模型参数的可靠估计;(2)由于环境噪声与接收函数在频带上的适配性明显优于地震面波,接收函数与环境噪声的非线性联合反演能更好地约束台站下方近地表的速度结构;对于周期范围为2~40s的环境噪声相速度频散,利用本文方法能够可靠推测台站下方0~80 km深度范围的S波速度结构, 其浅表速度结构的分辨率可达到1 km; (3)本文方法能够可靠地估计地壳泊松比,泊松比的全局性搜索有助于合理解释接收函数和环境噪声的面波频散数据.利用本文方法对川西台阵KWC05台站观测的接收函数与环境噪声的联合反演表明,该台站下方地壳厚度为44 km,上地壳具有明显的高速结构,24~42 km范围的中下地壳具有低速结构.该台站下方地壳的平均泊松比为0.262,壳内低速带的泊松比为0.27.     

利用接收函数研究渭河地堑及其周边地壳结构

利用接收函数方法对横跨秦岭造山带、渭河地堑及鄂尔多斯块体的15个地震观测台站下方的地壳结构进行研究分析,结果表明三种不同类型（造山带型、拉张盆地型和稳定克拉通型）的构造单元的地壳结构和物质组成存在明显的差异．秦岭北缘平均地壳厚度为37．8km,泊松比为0．247,相对偏低的泊松比表明地壳物质长英质组分增加．鄂尔多斯块体南缘平均地壳厚度为39．2km,泊松比为0．265,偏高的泊松比与鄂尔多斯下方古老的铁镁质结晶基底以及浅部沉积有关．通过接收函数正演计算表明低速的、厚度较大的松散沉积层对Mohorovicic不连续面（Moho）的震相具有较大影响,是渭河地堑内部台站的接收函数Moho转换震相不清楚的主要原因．S波速度结构反演结果表明渭河地堑上覆松散沉积层,其厚度约为4-8km,该沉积层使得位于渭河地堑内台站的接收函数Moho震相复杂．另外渭河地堑下方中下地壳位置存在一高速区域,该高速体可能与渭河断裂系统的活动有关．     

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

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

基于地震背景噪声的四川威远地区浅层速度结构成像研究

地震S波速度是防震减灾中场地分类、强地面震动模拟建模等的重要参数.地震背景噪声成像方法可以重建地下浅层三维结构以及探测浅层速度结构变化,为页岩气开采提供参考,有效规避地震灾害风险.本文收集了2015年11月至2016年2月间四川威远地区50个流动台站记录到的垂直分量连续波形记录,利用波形互相关提取2~6 s的基阶Rayleigh波经验格林函数,采用一维地壳结构模型拟合相速度频散曲线,获得了该区域近地表5 km以内的三维S波速度模型.结果表明,威远地区深度2 km以内的S波速度为2.0~2.7 km·s^-1;2.0~5.0 km深度的S波速度横向分布不均匀,西北地区大于2.7 km·s^-1,东南地区在2.3~2.8 km·s^-1之间,分别与背斜构造和较厚的盆地沉积层相对应.3个剖面图均表明S波速度随着深度的增加而逐渐增大,并且在台站S29附近的页岩气田钻井底部与2016年1月7日四川威远M_L3.9地震震源位置较为接近,均处于松散的沉积盖层与较为坚硬的花岗岩基底的分界处,推测它们之间可能存在一定的联系.     

沉积层结构被动源探测方法及其在典型盆地的应用

准确的沉积层结构,对盆地地区油气资源调查与勘探、地震动场地效应评估和壳幔深部结构成像等研究均具有重要意义.近二十年来,随着地震观测技术的进步和大量流动台阵数据的积累,利用被动源资料探测沉积层结构的地震学方法研究取得了长足的发展.本文对这些主要方法进行了归纳总结,阐述了基本原理及相关进展.其中,基于远震体波数据的研究方法主要有：接收函数和转换函数、波场反延拓的H-β方法、以及P波质点运动分析等.基于近场体波资料约束沉积层结构通常采用的是高频波形拟合方法.此外,还简要介绍了谱比法、背景噪声面波成像和Rayleigh波Z/H幅度比,以及多种方法的联合反演.最后,对我国东部典型盆地地区,如松辽盆地和华北盆地下方沉积层结构已取得的研究进展进行了总结.     

沉积盆地地区地壳结构估计——预测反褶积方法消除接收函数多次波混响

接收函数方法被广泛地应用于地壳上地幔结构的研究中,H-κ叠加方法是其中最常用的方法之一.对于布设在基岩区台站计算的接收函数,H-κ叠加方法可以准确地估计台站下方地壳厚度和平均波速比,但是对于沉积盆地地区计算的接收函数,由于低速沉积层内会产生多次波混响,干扰甚至覆盖接收函数中莫霍面的转换波和多次波震相,从而影响H-κ叠加结果的准确性.为准确估计沉积盆地地区地壳结构,本文提出使用预测反褶积方法去除接收函数中低速沉积层内多次波混响,其中预测步长由接收函数归一化自相关函数获得,物理意义为沉积层内S波双程走时.合成接收函数和实测接收函数试验表明,本文提出的预测反褶积方法可以有效地去除沉积层多次波混响,并结合改进的H-κ叠加方法可以准确地估计沉积层下覆地壳厚度和平均波速比.相比于其他去除接收函数多次波混响的方法,本文提出的预测反褶积方法具有参数设定简单、运算量小、震相幅值较大等特点,适用于大批量数据处理.     

珠江三角洲芦苞—佛山—中山—三乡断面的地壳结构——密集台阵勘探新成果

为了研究珠江三角洲核心地区中深部地壳结构,在芦苞—佛山—中山—三乡地区开展了密集台阵勘探工作,使用接收函数方法处理反演,获得珠江三角洲核心地区40 km以浅Vs速度结构.结合研究区的地质、重力异常特征对反演成果分析获得初步成果:珠江三角洲地区地壳结构可分为上地壳、中地壳、下地壳和上地幔四个大层位,上地壳厚度在临海一侧变厚,在中山三乡镇一带,厚达14 km;莫霍面隆起、上地幔物质侵入下地壳和高密度岩浆侵入上地壳这个三个因素共同导致研究区重力异常高;获得了沿线中深部断裂构造空间展布特征信息,速度结构和断裂构造具有较好的相关性;沿测线地壳在15～20 km深的深度上不存在低速层;三水盆地基底深度在3 km左右,与该地区沉积层厚度一致;莫霍面深度北深南浅,深度在27.5～30 km之间,其波速特征明显,由3.6 km/s突变为4.1 km/s.     

接收函数方法的研究综述

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

天水-礼县地区地壳速度结构

1984—1985年,利用厂坝铅锌矿工业爆破,在天水—礼县地区布设测线进行了大范围的地震测深工作。对该地区的地壳速度结构的研究结果表明,该地区沉积层平均厚度为2.5km,速度为4.0km/s(P波);地壳平均厚度为43.68km,平均速度为6.20km/s;徽县—礼县地壳速度剖面可分为5层,其中在24—29km深处有一低速层,基底深度变化较大,在礼县地壳浅部发现一断层。对天水—礼县地区还进行了P波、S波联合反演,获得了该区P波与S波速度结构,其地壳范围内的平均波速比为1.73。     

大兴安岭造山带及两侧盆地的地壳速度结构

内蒙新巴尔虎左旗-黑龙江齐齐哈尔深地震测深剖面长630 km,跨越海拉尔盆地、大兴安岭造山带和松辽盆地.本文根据沿测线爆破地震的9炮记录截面图中,5个震相的到时资料,结合地震记录中的振幅信息,确定了沿剖面的二维纵波地壳速度结构,海拉尔盆地的地壳厚度为39.0~41.0 km,大兴安岭造山带西侧莫霍面深度为38.5~43.5 km.东侧的莫霍面深度为34.5~36.4 km.松辽盆地的莫霍面深度为32.4~36.2km.整个地壳形态东浅西深,松辽盆地最浅的莫霍面深度为32.4 km,大兴安岭西侧最深的莫霍面深度为43.5 km.最后讨论了本区的深部特征和盆山结构关系.     

呼包盆地周缘壳、幔结构研究

对呼包盆地周边七个台站的远震接收函数研究表明：研究区地壳厚度为43～46 km,地壳速度比和S波速度结构均无异常;阴山造山带岩石圈埋深为65～85 km,且上地幔S波平均速度偏低,是典型的异常地幔区.根据S波速度结构和地壳厚度随地形高程反比的变化关系,以及地表广泛出露的幔源玄武岩分布,推测该区地幔深部热物质上涌是阴山造山带隆升的原因之一.上升的软流层物质与地幔发生交代、侵蚀作用导致岩石圈减薄,S波速度降低.呼包盆地的形成可能与深部物质上涌造成的拉张效应有关.     

P波和S波接收函数的贝叶斯联合反演

S波接收函数对于研究岩石圈速度结构具有重要价值. 本文利用合成地震图技术研究了S波接收函数的动力学特征. 在接收函数非线性复谱比反演方法的基础上,发展了基于贝叶斯理论的P波和S波接收函数的非线性联合反演方法. 结果表明:(1)适用于S波接收函数反演的震中距范围约为55°~80°,S波接收函数反演要求所用远震事件的震级大于5级; (2)与陡变的岩石圈底部界面(LAB)相比,梯度带类型LAB上生成的S_LP转换波相对较弱,台站下方的沉积盖层有助于相对增强S_LP震相; (3)由于S波接收函数径向分量不符合δ脉冲,不依赖于等效震源假定的三分量接收函数多道最大或然性反褶积方法更适合S波接收函数的估计;(4)数值检验的结果表明,在初始模型速度参数偏离真实模型20%的情况下,本文的方法能够预测300 km深度范围内的P波和S波速度结构;(5)观测数据的反演结果表明,由于P波接收函数低频分量相对不足,本文的联合反演方法对于大于100 km深度上地幔的S波速度结构约束相对较弱.     

面波反演苏北—南黄海盆地地壳三维S波速度结构

区域地壳速度结构对于地震定位和地球动力学特征认识十分重要,一直是地震反演研究的主要内容。利用收集到的基于背景噪声面波成像方法的瑞利面波相速度数据,我们对苏北—南黄海地区地下三维S波速度结构开展了深度反演,进而探讨苏北—南黄海盆地地质结构与地震构造的关系。面波层析成像结果表明区域速度结构模型与地质特征基本一致,速度结构清晰刻画出了苏北—南黄海沉积盆地的主要沉积构造范围和基底起伏。反演结果还揭示了南黄海盆地中地壳内存在一显著的低速层,与历史强震集中分布区高度一致,表明浅源强震活动与壳内的低速层密切相关。     

基于面波走时的三维结构面波直接成像：方法综述与应用

利用面波频散数据进行各向同性和各向异性横波速度结构成像是研究全球和区域构造变形、获取地壳上地幔和近地表精细结构的一种非常有效的方法.传统的基于频散数据的面波成像通常采用两步法,即首先反演二维相速度或群速度分布图,然后再基于格点的纯路径频散反演格点下方的一维横波速度模型,之后再拼合成三维速度模型.在本文中我们较为详细地综述新发展的基于面波频散走时的三维面波直接成像方法体系,即面波一步法成像.该方法体系包括通过所有路径、不同频率的面波频散走时直接反演三维各向同性横波速度结构的方法（DSurfTomo）、直接反演三维方位各向异性横波速度结构的方法（DAzimSurfTomo）,以及直接反演三维径向各向异性横波速度结构的方法（DRadiSurfTomo）.新的面波直接成像方法计算不同周期的面波射线路径,从而更好地考虑复杂结构下面波射线路径弯曲对成像精度的影响.随后我们介绍面波直接成像方法的应用实例,包括地壳上地幔及地壳浅部的多尺度各向同性和各向异性横波速度结构的成像,这些成像研究为认识区域构造演化、孕震构造、断裂带浅部结构、成矿构造、城市地下结构等提供了重要的约束.最后我们讨论面波直接成像方法...