川西地区台阵环境噪声瑞利波相速度层析成像

2006年中国地震局地质研究所地震动力学国家重点实验室在川西地区(26°N~32°N,100°E~105°E)布设了由297台宽频带数字地震仪组成的流动观测台阵.利用该密集台阵29°N以北156个台站2007年1~12月份的地震环境噪声记录和互相关技术,我们得到了所有台站对的面波经验格林函数和瑞利波相速度频散曲线,并进一步反演得到了观测台阵下方2~35 s周期的瑞利波相速度分布图像．本文结果表明,观测台阵覆盖的川滇地块、松潘-甘孜地块和四川盆地的地壳速度结构存在显著差异,具体表现为：(1)短周期（2~8 s）相速度分布与地表构造特征相吻合,作为川滇地块、松潘-甘孜地块和四川盆地之间的边界断裂,龙门山断裂带和鲜水河断裂带对上述三个地块上地壳的速度结构具有明显的控制作用,四川盆地前陆低速特征表明相应区域存在较厚的（约10 km）沉积盖层;(2)中周期（12~18 s）相速度分布表明,川滇地块和松潘-甘孜地块中上地壳速度结构存在明显的不均匀横向变化,并形成了尺度不同且高、低速相间的分块结构,而四川盆地中地壳整体上已经表现出相对高速;(3)长周期（25~35 s）相速度分布表明,松潘-甘孜地块,特别是川滇地块中下地壳表现为广泛的明显低速异常,意味着它们的中下地壳相对软弱,而四川盆地的中下地壳呈现整体性的相对高速,意味着四川盆地具有相对坚硬的中下地壳,并且以汶川地震的震中为界,龙门山断裂带的地壳结构显示了北段为高速异常,南段为低速异常的分段特征．     

用Rayleigh面波方位各向异性研究中国大陆岩石圈形变特征

中国大陆地质构造历史非常复杂,岩石圈长期积累的形变较大,而利用地震面波传播的各向异性是研究岩石圈形变特征的强有力手段. 本文利用双台窄带通滤波-互相关方法与基于图像分析的相速度频散曲线提取技术,提取Rayleigh面波相速度频散资料,进而反演中国大陆及邻区20～120 s周期Rayleigh面波相速度方位各向异性空间分布图像. 检测板测试结果显示：中国大陆大部分区域的方位各向异性横向分辨率在5°左右. 各向异性研究结果表明：中国大陆地壳上地幔方位各向异性特征存在显著的空间差异,反映出形变特征的空间差异;104°E以东地区地壳上地幔各向异性弱于西部地区,表明其构造变形总体弱于西部地区. 青藏地块及其东缘地区地壳与上地幔顶部变形最为强烈. 但东部的局部地区如华南地块与珠江口地区、鄂尔多斯盆地西南缘以及秦岭-大别造山带,较强的各向异性显示这些区域在不同时期也经历了强变形. 青藏地块内中短周期快波方向自西向东顺时针旋转变化可能指示板块碰撞与挤压过程中软弱物质的流变方向. 青藏地块西部中下地壳和上地幔形变模式相似,可能处于壳幔耦合状态;而中东部及东缘地区地壳上地幔形变模式存在明显差异,壳幔似乎不具备垂直连贯的形变特征. 位于青藏地块北部的塔里木盆地、柴达木盆地以及祁连褶皱带同样经历了强变形. 包括四川盆地在内的上扬子地块快波方向的变化显示中地壳与下地壳上地幔形变模式不同,而形变特征一致的下地壳与上地幔应为强耦合. 大约以103°E为界,龙门山断裂带可分为南西段和北东段,南西段处于低速区,而北东段位于高速区,且方位各向异性强度明显大于南西段;2008年5月12日汶川M_S8.0级地震沿断裂带的单侧破裂模式除与北东段的高应力积累有关外,还可能与北东段地下介质物性存在密切关系,高速坚硬岩体的发育有利于应变能的积累与集中释放.     

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

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

基于图像分析的双台面波相速度频散曲线快速提取方法

双台波形的互相关方法是面波相速度频散曲线提取的一种常用方法。传统的方法精度不高,效率较低。本文提出了一种基于图像分析的频散曲线快速提取方法,采用图像显示互相关振幅矩阵,且提出了快速追踪整条频散曲线的算法,并用Matlab编写了交互式的处理软件。基于这种方法编写的资料处理软件提高了双台相速度频散曲线提取的速度和精度,为大批量处理地震数据资料提供了坚实的软件基础。基于图像分析的方法可以更为清晰地展现出双台问的面波频散曲线的特征,对频散曲线的识别和分析也提供了更为有效的研究工具。     

基于长江安徽段气枪主动源实验庐江台阵的台站对双差地震速度成像

针对气枪主动源监测中相位加权叠加数据绝对到时拾取可能存在较大误差的问题,提出基于台站之间叠加信号波形的相似性,利用台站对波形的互相关得到更准确到时差的策略。在此基础上,将台站对双差地震速度成像方法直接用于台站对到时差数据,以确定地下介质的速度异常。本文利用长江安徽段气枪实验在庐江临时地震台阵接收到的数据进行了测试,发现基于波形互相关得到的到时差数据比绝对到时数据更准确。最后,采用台站对双差地震成像的方法,得到了庐江台阵覆盖区域的地壳三维速度异常模型,继而确定了该地区较大尺度的地下结构分布。     

基于伴随方法的线性台阵背景噪声面波和远震体波联合成像研究

伴随层析成像（Adjoint Tomography）通过求解全波方程来准确模拟地震波在复杂介质中的传播,并利用波形信息来反演地下结构,是新一代的高分辨率成像方法.其中3-D伴随层析成像需要庞大的计算资源,而2-D反演相对更具计算效率.面波和远震体波是研究地壳上地幔速度结构的重要方法,它们对S波速度及Moho面的敏感度不同,通过联合反演,可以得到更为准确的S波速度结构及Moho面.通过两种数据的高度互补性,本文提出基于伴随方法的线性台阵背景噪声面波和远震体波联合成像方法,同时约束台阵下方S波速度结构及Moho面形态.我们将该方法应用到符合华北克拉通岩石圈典型结构特征的理论模型上,测试结果表明联合反演方法优势明显,相比于面波伴随层析成像,能获得更高分辨率的S波速度结构,同时能精准约束Moho面形态.相比于体波伴随层析成像,联合反演能有效压制高频假象,降低波形反演过程中的非线性化程度.本研究有望提供一种更为高效精准的线性台阵成像方法,搭建联合伴随层析成像理论框架,提升岩石圈成像分辨率,并为后续其他类型波形数据的引入提供思路和方法.     

气枪震源资料反褶积方法及处理流程研究

不同工作条件下气枪震源产生的信号会存在细微差异,反褶积方法能在一定程度上消除由震源变化引起的记录信号变化。为了去除气枪震源子波信号,获取气枪源到台站之间的格林函数,通常需要选取一种恰当的方法对地震波形数据进行反褶积处理。频率域水准反褶积和时间域迭代反褶积是在接收函数等领域已被广泛使用的2种反褶积方法。本文以云南宾川主动源资料为例,对比了利用这2种方法处理气枪震源信号的效果,结果表明,在计算效率方面,频率域水准反褶积方法更具优势;在处理结果的信噪比方面,时间域迭代反褶积方法表现更好,P波初至也更清晰。此外,进一步讨论了在多炮资料的处理过程中反褶积和叠加等操作的顺序问题,最后提出了从气枪震源资料中提取气枪源到台站之间的格林函数的一般流程。     

2021年5月漾濞MS6.4地震序列发震机理的应力分析

为了从应力的角度解释2021年5月漾濞M_S6.4地震序列的发震机理,本文开展了包括刻画断层结构,反演震源区应力场,研究应力场对断层的加载作用以及地震序列中较大地震对断层活动的库仑应力影响等4个方面的研究,具体的认识为,(1)相比于前人工作中利用区域台阵定位结果刻画的断层结构,结合布设在漾濞震源区的实时地震监测台阵捕捉的7905个余震位置和沿断层均匀分布的16个M3.5以上地震的震源机制,我们刻画了漾濞地震序列发震断层的精细结构：漾濞地震序列最主要的特征为NW向和NE向两组断层共同破裂,NW向断层为右旋走滑破裂,NE向断层为左旋走滑破裂;(2)利用震源区1°×1°范围内的90个历史震源机制反演了震源区的背景应力场,得到震源区为近南北向水平挤压,东西向水平拉张的走滑型应力状态,应力方位的不确定度很小,而主应力相对大小(R)的不确定度很大,为0.2～0.8;(3)分析了漾濞地震序列两组发震断层的稳定性,发现两组发震断层均为应力场作用下最易失稳的断层.此外,发生在两组断层上地震的震源机制分别与两个最易失稳断层的震源机制相似,说明了背景应力对断层行为的控制作用;(4)两组断...     

利用区域地震体波走时评价华北地区三维地壳速度结构模型

本研究利用国家地震台网131个地震台站2009—2016年记录的1 749次近震的初至P波和S波走时数据,与采用快速行进方法正演计算得到的华北地区4个三维地壳速度模型对应的走时数据进行对比,通过统计分析的方法,评价这4个速度模型与真实地下结构的近似程度。结果表明:4个速度模型在大范围内存在较高的一致性,在整个研究区内（111°E—119.5°E,37°N—42°N）,Shen等的模型（简称“S模型”）相对优于Fang等的模型（简称“F模型”）和Duan等的模型（简称“D模型”）,Laske等的Crust1.0模型（简称“C模型”）相对较差。我们认为该结果与上述几个模型所使用的数据及其分辨率有关。对于研究区域内的构造单元,D模型在燕山褶皱带西南部、太行山山前构造带西北部和沧县隆起区表现较好,F模型在太行山隆起区中部、沧县隆起北部、黄骅凹陷区和燕山褶皱带表现较好,S模型在西部地块、山西凹陷区、太行山山前构造带和冀中凹陷区表现较好,C模型无明显连片表现较好区域。      

合肥市地壳浅部三维速度结构及城市沉积环境初探

合肥市位于合肥盆地东南缘,东侧紧邻郯庐断裂带,多条大型隐伏断裂穿过市区.为进一步认识合肥城市下方隐伏断裂的空间展布、性质,以及城市复杂的沉积环境,本文利用布设在合肥市区的58套三分量短周期地震仪组成的台阵,获得了37天的三分量连续波形数据,通过基于射线追踪的面波走时直接成像方法反演得到了合肥市地壳浅部0.6~3.6 km的三维剪切波速度结构,速度结构图像展现了地壳浅部的横向不均匀性和纵向成层性,揭示出NNE、NWW和近EW三组不同走向的隐伏断裂在城市地下浅部的构造特征.取得以下认识：（1）合肥市南、北方向在浅地表（2 km以内）存在显著速度差异,速度分界线位置与已知的近EW向的蜀山断裂一致,断裂南侧呈现低速凹陷,北侧则为高速隆起.低速中心深度达2~3 km,速度异常与该断裂在合肥盆地东部演化过程中的构造反转沉积了不同地层有关;（2）合肥市区存在明显的高速异常带,其走向、位置与穿过该区域的郯庐断裂带西支主干断裂相符,其中五河—合肥断裂在市区北部以东呈现低速凹陷特征,低速区范围与肥东凹陷晚白垩纪以来的沉积构造边界一致,认为肥东凹陷的最大沉积厚度可达2 km以上;（3）合肥市中心跨郯庐断裂带西支主干断裂之间呈现明显的凹、隆相间的复杂构造,推测其是在多组断裂的共同拉伸作用下形成的小型沉积盆地,沉积中心位于郯庐断裂带内部,最大厚度可达3~4 km.由于其展布方向在不同深度与该区域断裂的走向具有明显的相关性,推测不同深度的沉积形态与郯庐断裂带在不同时期的构造演化过程有关.