利用面波和接收函数联合反演滇西地区壳幔速度结构

考虑到面波频散对介质S波速度、接收函数对界面深度的各自敏感性优势,综合利用面波和接收函数资料实现联合反演,求取滇西地区壳幔速度结构. 本文利用适配滤波频时分析技术处理覆盖滇西地区的长周期面波资料,获得105～1050s周期范围内的面波群速度频散,进而利用分格反演方法提取研究区内1°×1°网格纯路径频散;基于滇西地区宽频带三分量远震记录,经反褶积后得到台站下方的远震P波接收函数. 联立面波纯路径频散信息和接收函数资料建立系统方程,利用阻尼最小二乘法实现联合反演,从而获得滇西地区壳幔S波速度结构. 结果表明,滇西地区以红河断裂为界,东西两侧壳幔结构存在明显差异,断裂西侧约20km深度处存在一厚度为10km左右的低速层,而东侧并不明显;滇缅泰块体上的畹町、沧源一带属于上地幔低速区,而另一个地幔低速区则位于滇中块体上的康滇古隆起上,两处地幔低速区与大地高热流分布、强震活动具有较好的对应关系.     

汶川地震破裂带地壳速度结构研究

利用分布在汶川地震破裂带上三个固定台站两年的远震接收函数,结合该区域已有面波频散结果,首先确定了各台站下方不同方位地壳的厚度、波速比和纵波平均速度,进而构建初始模型,对各台站下方不同方位地壳速度结构进行了线性反演.结果表明:位于汶川地震破裂带两端单台下方地壳结构的横向非均匀性较强,各方位接收函数差异较大,而破裂带中部台...     

大三江盆地及邻区地壳结构研究

收集了大三江盆地及其邻区区域地震台网及多个流动台阵的连续波形及远震事件资料,采用背景噪声层析成像和接收函数叠加方法,分别获得了研究区三维S波速度结构、基底及莫霍面深度和泊松比.结果显示:浅层速度结构较好地反映了地表地形及地质特征,三江盆地呈明显的低速,虎林和勃利等小型盆地的S波速度也相对较低,而小兴安岭、张广才岭等则呈...     

山西断陷带地壳结构的接收函数研究

利用2006年8月到2008年3月北京大学在山西断陷带南部架设的两条东西向流动观测地震台阵记录的远震资料,提取各台站接收函数,然后进行倾斜叠加(Slant stack)和台阵偏移成像,获得了沿台阵横穿山西断陷带的地壳和上地幔的精细结构变化.研究结果显示,山西断陷带下面莫霍面存在明显不连续,莫霍面上隆约4~6 km,纵横波速比从两侧的1.75上升为山西断陷带内部2.0左右,且中、下地壳可能存在一个低速层.山西断陷带的构造模式沿相距140 km的两条剖面表现出明显差异:南端的临汾盆地为拉张作用下的纯剪切模式,向北转化为太原盆地的简单剪切模式.     

基于背景噪声成像方法的新疆呼图壁储气库地区近地表速度结构研究

近年来,背景噪声成像方法在恢复高频面波信号及获取近地表速度结构方面得到了广泛的应用,本文将该方法应用于准噶尔盆地南缘的呼图壁背斜地区的呼图壁储气库.采用储气库及其周边区域22个台站记录的连续背景噪声数据的垂直分量,通过噪声互相关方法获得了台站对之间的瑞利面波经验格林函数,并进一步提取了0.5~1.5 s的基阶瑞利面波群速度频散曲线.首先根据区域平均频散曲线得到了该地区地下数百米的平均一维横波速度结构,然后利用基于面波射线路径追踪的面波频散直接成像方法得到该地区深度为500 m以上的三维横波速度结构.反演结果显示该地区沉积层较厚,整体横波速度值较小（0.4~0.9 km·s^-1）.储气库在地表投影区域的横波速度值较小,这可能是由于抽注水、气引起的沉积岩石裂隙所导致.储气库东北和东南方向均有明显的相对高速区,推测是区域地下水位和地形起伏综合作用的结果.本研究获得的近地表三维速度结构为呼图壁储气库地区的上覆地层物性研究、区域微震精定位、场地效应的评估和去除浅层影响的深部介质成像等研究提供了重要基础.     

面波频散与体波接收函数联合反演研究回顾及展望

面波频散和接收函数是研究地壳上地幔速度结构的重要方法.面波频散对绝对S波速度较为敏感,但对界面结构约束不强;而接收函数对间断面具有很好的分辨能力.二者优势互补,通过联合反演,可以得到更为准确的壳幔内部界面及绝对S波速度结构.近十几年来,随着噪声地震学的迅猛发展以及大量宽频带地震仪的架设,该联合反演方法得到了更进一步的应用.为了便于梳理该研究方向的脉络,本文分别回顾了面波频散与体波接收函数的方法原理及其在壳幔结构反演中的应用,并分析了两种方法的优势与不足.在此基础上,介绍了二者的联合反演方法,以及在地壳上地幔速度结构研究中的应用.展望了联合反演方法的进一步发展,以及加入多种地球物理观测数据对反演地壳上地幔结构的意义.     

基于矢量波数变换法(VWTM)的多道Rayleigh波分析方法

在近二十年来,多道面波分析法(MASW)由于其便捷、高效等特性在浅层地震勘探领域得到了广泛的应用.本文基于多道面波勘探的采集方式,提出了一种新的面波多道分析方法——矢量波数变换法(VWTM).该方法通过对震源的近似,基于水平层状模型得到台站与震源间近似格林函数,然后进行矢量波数变换得到含有高阶模态Rayleigh波(频率-相速度)频散能量图.本研究首先利用合成地震数据到频散能量图与理论频散曲线进行叠加分析该方法的有效性和正确性;然后与相移法进行对比分析,我们发现在频散能量图中VWTM法对基阶、高阶模态成像均具有更高的分辨率和成像质量;最后我们将其应用于实际多道瞬态面波探测中,通过与相移法进行对比分析,发现VWTM法是一种方便、实用、有效的Rayleigh波频散提取方法.VWTM法提取多模态的Rayleigh波频散特征具有巨大潜力,可为基阶、高阶面波频散联合反演提供丰富的高阶模态频散信息.     

体波接收函数与瑞利波椭圆率联合反演研究

近年来,结合体波接收函数或瑞利波椭圆率,面波频散已被广泛应用于岩石层结构的成像研究中。由于面波频散是其传播路径的总传播效应,因此联合反演方法需依赖于密集的地震台阵或较高的地震活动性,从而获得局部速度结构。但是,接收函数和瑞利波椭圆率均为具有局部敏感性的单台测量方法,并且能自然地相互合并进行联合反演。本文中,我们对利用瑞利波椭圆率和接收函数联合反演的可行性进行了探讨。我们通过正演模拟进行了敏感性检验,发现接收函数对速度跃变的界面较为敏感,但对长波结构显示较弱的敏感性,这几乎与瑞利波椭圆率互补。因此,利用两种单台测量方法进行联合反演为地震台站下方的速度结构提供了更严格的条件约束。本文提出了一种基于快速模拟退火方法的联合反演算法,对岩石层结构的瑞利波椭圆率和接收函数进行反演。该算法在印度克拉通和美国威利斯顿盆地的应用中,在减少反演的非唯一性方面证明了其有效性。但是,通过该联合反演可能无法获得有效的地壳平均速度,表明有必要结合面波频散(或其他类型的观察研究)、接收函数和瑞利波椭圆率来获得更准确的速度结构。     

接收函数方法及研究进展

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

利用瑞雷面波反演华北克拉通东北部边界的岩石圈结构

运用瑞雷面波相速度频散曲线分析和反演得到了华北克拉通东北部边界及其邻近区域岩石圈的精细S波速度结构.利用11个地震事件、60个台站的瑞雷面波波形资料,得到了周期从25-150 S的相速度频散曲线,并且通过线性反演方法得到了深度从40-300 km的S波速度结构.结果表明,该研究区域S波速度存在强烈的不均匀性.从东南部的...     

阿尼玛卿缝合带东段地壳结构的接收函数研究

沿四川红原至甘肃武威一线布设了20个宽频带地震台站， 在一年的观测时间里共接收到81次远震记录， 利用台站记录的远震P波波形数据和接收函数方法， 获得了沿测线的接收函数剖面、 每个台站下方的S波速度结构. 研究结果表明， 研究区地壳速度结构复杂， 整个地壳的平均S波速度偏低. 在四川阿坝弧形断裂——秦岭地轴北缘断裂之间， 地壳中10~40 km的深度范围内普遍存在低速层， 是阿尼玛卿古特提斯洋从闭合、 斜向碰撞到俯冲板块折返或逆冲岩片抬升等复杂地质过程所形成的构造特征. 沿剖面莫霍面深度约为50 km， 南边略深北边略浅.      

利用地脉动信息约束沉积层区域台站下方速度结构反演

利用接收函数和地震面波频散联合反演台站下方速度结构,已成为一种常用技术.联合反演结果能同时匹配两个不同数据集合,使得解的非唯一性在一定程度上得到了有效抑制.然而对于现有地震台站分布,短周期面波由于受多种因素的影响,可以获得的有效频散资料较少,层析成像的横向分辨率较低,难以反映台站附近局部介质真实的频散特性,而且这些面波...     

接收函数的曲波变换去噪

压制横向非均匀地壳介质引起的散射波场对于基于水平分层介质模型的接收函数地壳结构成像及其地震各向异性研究至关重要.虽然通过数据叠加和低通滤波,在一定程度上能够压制散射波场,但也有可能导致不希望的波形畸变、信息丢失或数据分辨率降低.为了避免上述问题,本文将近年来快速发展的曲波变换理论用于远震接收函数的散射噪声压制.与勘探地震学不同,我们面临的主要问题在于地震台站和震源的空间缺失导致的接收函数空间不均匀采样.为此,我们将压缩感知理论与曲波变换去噪相结合,在对缺失数据进行波场重建的同时,实现散射噪声的压制.为了论证方法的可行性,本文进行了噪声压制和波场重建的理论检验.并将本文方法用于处理IRIS全球台网固定台站和川西台阵远震接收函数.结果表明:1)地壳介质横向不均匀引起的散射噪声可以得到有效压制,接收函数的信噪比得到提高,震相的可追踪性得到改善,从而利于进一步的接收函数反演和地震各向异性研究;2)缺失数据可以正确重建;3)本文的方法既可用于单台-多事件的数据集,也可用于单个事件-阵列观测的数据去噪,但单台-多事件数据集的结果优于阵列观测的情况.     

基于互相关相移的主动源地震面波频散成像方法

以多道面波分析技术(MASW)为代表的主动源面波方法已经在地表勘探领域获得越来越多的应用.准确的面波相速度频散成像是利用主动源面波法进行近地表横波速度反演的必要条件.本文提出了一种新的主动源面波频散成像方法—CCPS(Cross-Correlation and Phase-Shifting)法.首先,将单炮面波记录由时空域变换到频率空间域并消去振幅谱的影响,得到归一化频谱.然后,依次将每道信号与整个记录进行互相关的同时进行相速度扫描叠加.最后经过一次求和过程即得到最终的频散图像.通过理论数据、实际工程地震面波资料和石油勘探面波资料比较了普通相移法和CCPS法的频散成像效果.研究表明:新方法相对于普通相移法而言,明显提高了频散图像的质量,相速度聚焦性能和背景噪声压制能力有显著改善.     

利用双平面波干涉面波层析成像方法研究山西断陷盆地及鄂尔多斯地台三维速度结构

本文利用了北京大学在山西地堑的34个台站以及中国地震局台网中心在鄂尔多斯地区46个台站的远震波形数据,运用双平面波干涉的面波层析成像方法,提取瑞利面波相速度频散曲线,开展台阵覆盖区三维速度结构反演,据此分析了研究区地壳和上地幔的结构特征。结果表明,瑞利面波相速度分区特征显著,山西断陷盆地和渭河断陷盆地的相速度整体偏低。鄂尔多斯块体在莫霍面以下有明显高速异常,表明该地块为构造稳定的克拉通块体,鄂尔多斯块体的岩石圈下界面在120～140km的深度左右。与此相反,山西断陷盆地和渭河断陷盆地地下70～120km左右均有低速异常,显示这两个地区构造活动活跃,这也与该区域历史上多次发生强震是相符的。     

依舒断裂带北段莫霍面特征与泊松比

伊舒断裂带北段是黑龙江省地震跟踪监测重点区域之一,具有较强的地震活动性,为主要历史强震发生区域。利用该区域4个固定地震台记录到的远震地震事件资料,采用频率域反褶积方法,提取远震P波接收函数,再利用H—κ叠加方法计算台站下方莫霍面深度和泊松比。结果表明：受太平洋板块俯冲作用的影响,块体间相互挤压,双鸭山地震台下方莫霍面发生了错断,并且一侧莫霍面叠置在另外一侧上方,鹤岗地震台西侧莫霍面上方存在1个速度间断面;区域内泊松比值较高,为0.266—0.277,反映了上地幔铁镁物质上涌到下地壳中。     

利用多种地震数据联合反演剪切波速度结构的可靠性检测

本文模拟使用青藏高原东南缘区域台网及国家台网的170个宽频台站基于背景噪声、天然地震面波、P波接收函数反演时的实际数据,对青藏高原东南缘假定的初始模型进行恢复,通过计算初始模型台站下方纯路径频散、提取各台站对间的瑞雷波频散曲线、计算理论接收函数以及反演剪切波速度结构来测试使用不同单项数据与联合使用多种数据反演对初始模型的恢复程度。结果表明,同时使用接收函数、基于噪声经验格林函数的群速度、相速度频散以及基于天然地震面波的相速度频散联合反演的剪切波速度结构,充分利用了几种数据的分辨率优势,清晰地分辨出中下地壳及上地幔顶部的低速层。此外,本文也分析了实际数据处理中出现的计算误差、随机噪声干扰对计算结果稳定性的影响。结果显示:对于面波频散,加入1%的误差后,联合反演的结果仍可很好地反映低速层的形态,但是当误差提升至5%后,对最终结果则产生了一定程度的影响;而在接收函数中加入4%的随机噪声时,虽然地幔低速层的上界面和下界面会略微受到随机噪声的影响,但是低速层的深度范围和速度值均得到了较好的恢复。     

欧亚大陆及西太平洋边缘海瑞利面波层析成像研究

利用亚洲、欧洲及西太平洋边缘海地区台网114个数字地震台,搜集了400W～1800E,400S～800N范围内近万个地震事件,并从中挑选出1982～2006年间震级在5.0～7.0之间、震源深度小于100㎞的1800余个事件.经过筛选,共得到周期在8～200s之间10450余条质量较高的瑞利面波频散曲线,这些面波所通过的大圆路径很好地覆盖了欧亚大陆及西太平洋.本文首先对原始记录进行预处理,测定面波频散曲线,包括截取Rayleigh面波,进行频率时域分析等,以得到实际频散数据.然后对研究区域进行网格划分,反演网格点附近的纯路径频散曲线.即采用网格化纯路径频散反演方法得到各块的纯频散数据.最后由得到的每块纯频散数据作为观测数据,并利用LSQR、SVD等广义线性反演方法,对每一块的纯频散曲线进行反演求取该单元的S波速度随深度的分布.     

地震背景噪声成像研究综述

地震背景噪声成像方法已成为21世纪地震学的伟大突破之一,其原理是,对2个台站记录的连续背景噪声信号进行互相关计算,得到台站间的格林函数,利用经验格林函数得到面波速度函数,获取面波频散特性,用地震层析成像方法对面波速度进行反演,得到地球内部的速度结构。其应用近年来日益广泛,涉及各向同性和各向异性的速度结构成像、大地震前后速度结构变化监测、体波联合成像、衰减结构成像、地震定位精度提高、噪声源分布和物理起源探究、强地面运动评估等。其优点是不需要等待天然地震或使用对环境构成威胁的人工爆破,所有台阵都可以当作源,拓宽了频带范围,可以获得较多短周期频散数据,提高反演分辨率。     

面波频散曲线提取方法对比分析

高频面波方法是以瑞雷波和勒夫波为研究对象,在水文、工程和环境研究中具有广泛的应用,其主要思想是利用面波的频散特性来反映浅层地质问题,提取高质量的频散曲线成了面波勘探中的一个重要步骤.本文通过对理论模型研究说明,利用τ-p变换与F-K变换两种方法均能有效地进行频散能量成像,其效果在高频范围差别不大,在低频范围F-K变换提取频散曲线更为精确;与瑞雷波记录相比,勒夫波记录信噪比较高,通过对同一地区采集的瑞雷波和勒夫波进行频散分析可知,瑞雷波记录频散能量成像效果差,频散能量出现间断、跳跃的现象,而勒夫波记录频散能量连续、集中,相速度分辨率较高,提取频散曲线简单清晰,可为反演地表横波精细速度结构提供更好的频散曲线数据.