新疆伽师强震群区三维地壳上地幔S波速度结构及其地震成因的探讨

１９９７年１月２１日至４月１６日（北京时间）新疆伽师地区连续发生了７次６级以上 强烈地震．为了深入研究该强震群形成的构造背景,我们在伽师强震群区及其邻域大约 １１ ０００ｋｍ^２的范围内布设了由 ３０台宽频带数字地震仪组成的流动地震台阵,台站的间距约５－ １０ｋｍ．利用远震体波接收函数的叠加偏移分析及其非线性反演技术,研究了台阵下方０－ １００ｋｍ深度范围内地壳、上地幔的三维Ｓ波速度结构．结果表明：１)Ｍｏｈｏ界面的深度在塔里 木盆地一侧为４０-５２ｋｍ,而在靠近天山一侧的褶皱变形区为６０-７６ｋｍ,南天山山前折皱变形 区的地壳厚度明显增大,地壳内部产生了明显变形;２)塔里木盆地北缘存在明显的不均匀性, 塔里木盆地与天山之间的接触变形关系显示了塔里木盆地向北北西方向的挤压作用,台阵下 方壳幔界面的地形与地表的地貌特征有较好的相关性,山前弧形褶皱带的形成与地下深部结 构及挤压作用有关;３)伽师强震群分布在壳幔界面梯度带的上方,该震群的成因可能与震源 区附近的隐伏断裂活动有关;４)由于震源处于褶皱变形区的地壳上部,相应的剪切模量较小, 这可能是伽师强震群应力降明显偏低的...     

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

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

龙门山断裂带中上地壳速度结构的短周期环境噪声成像

本文利用在龙门山断裂带周边布设的57个台站自2008年11月至2009年11月为期一年的垂直分量连续地震环境噪声数据,通过短周期地震环境噪声成像方法,获得了龙门山断裂带中北段地壳25km深度范围的S波精细速度结构.结果表明:(1)龙门山断裂带周边区域10km以上的速度结构与地表断裂的分布形态具有良好的一致性,速度结构控制了龙门山主要断层的深部延展特征;在15km及以下深度,S波速度结构呈现沿龙门山和沿岷山隆起走向的交叉构造格局,由此造成的速度结构差异可能影响了汶川地震的破裂过程;(2)速度结构随深度的分布特征为龙门山断裂带主要断层的深部延伸形态给出了良好的约束,结果进一步确认了龙门山断裂中段的高角度铲型断裂构造特征;(3)研究区的南端发现了龙门山断裂下方20km以下深度具有与松潘地块中地壳低速层相关的低速结构的迹象,这可能是汶川地震破裂带南段22km左右深度存在脆韧转换带的一个证据.研究结果显示出密集台阵和短周期环境噪声成像方法在地壳浅部精细结构和断层探测研究中具有巨大潜力.     

接收函数曲波变换去噪与偏移成像

增强接收函数偏移图像的垂向分辨率意味着提高参与叠加的接收函数的频率,但是采用高频接收函数通常伴随着对接收函数质量和参考速度模型的更高要求.通过叠加处理可去除部分接收函数中的随机噪声干扰,但同一台站的接收函数之间经常存在难以通过简单叠加消除的噪声信号.压制接收函数随机噪声的干扰可加强成像效果和提高图像分辨率,对推进叠加偏移成像质量的提高有重要的实际意义.本文利用在川西地区布设的31个流动台站所记录的远震波形数据,使用曲波变换去噪后信噪比增强的接收函数进行共转换点叠加(CCP),获得沿北纬31°线下方800km深度范围内速度间断面图像.研究结果表明:(1)对接收函数进行曲波变换去噪,可压制随机噪声,增强转换震相的追踪性,提高数据信噪比;(2)通过去噪处理,大幅提高接收函数用于偏移成像的主频率;(3)偏移结果确认了接收函数反演得到的松潘和川滇块体下方具有厚度约10～20km的过渡性Moho的认识;(4)上地幔过渡带的结果预示在龙门山断裂带以西的小范围内有可能存在下地壳或上地幔物质的拆沉.     

各向异性介质中S波接收函数成像方法

岩石圈各向异性主要由上地幔矿物晶格优势排列方向和上地壳内裂缝、裂隙的定向分布造成.在各向异性特征显著区域,利用SKS震相剪切波分裂获得的延迟时间高达1.5s以上.本文根据方位各向异性,利用广义反射透射系数矩阵方法正演S波接收函数,研究各向异性对不同反方位角接收函数转换震相走时的影响.我们发展了基于HTI模型各向异性走时校正的方法,成功的在单层和多层(快轴方向相同或不同)的各向异性介质中对齐不同反方位角接收函数的Moho面和LAB的转换震相走时.我们将该方法应用于在青藏高原东北缘的流动台站,试图在各向异性强度较大区域对实测数据转换震相走时的校正效果进行测试.结果表明:各向异性走时校正能够加强单台接收函数转换震相的可追踪性,能量增强的叠后转换震相在时深转换后更利于对界面深度的识别与判断;在考虑三维成像的情况下,我们的各向异性校正方法对提高成像结果的准确性有重要意义.     

SeisProc：基于ObsPy库的地震数据处理软件

在已有地震数据处理软件的基础上,为了进一步丰富和完善地震数据处理工具库,基于ObsPy模块库,设计并实现了一款简单易用、可拓展的地震数据处理软件SeisProc。该软件利用ObsPy提供的信号接口,实现了对多种格式的地震数据文件的预处理功能,应用PyQt5库设计GUI,实现数据预览、格式转换、波形绘制和预处理等功能。     

青藏高原东北缘上地幔多尺度层析成像

利用布设在青藏高原东北缘地区的甘肃宽频带地震台阵记录到的远震P波走时数据,采用小波域参数化和基于L1范数的稀疏约束反演算法的多尺度层析成像方法,得到了该地区400km深度范围内上地幔的P波速度结构.本文采用的多尺度层析成像方法可以自适应数据非均匀采样的情况,有效降低谱泄漏效应和反问题的多解性,明显提高解的分辨率和可靠性.层析成像结果表明青藏高原东北缘上地幔整体上显示为低速特征,扬子地块上地幔则显示为高速特征,两者之间上地幔存在清晰的块体边界带,该边界带位于东经104°—105°之间并且随深度的增加逐渐东移.该特征暗示了青藏高原上地幔物质向东扩张的机制,但在西秦岭上地幔顶部不存在物质运移的通道.青藏高原东北缘内部也具有明显的分区特征,松潘—甘孜地块上地幔P波速度整体呈低速特征,而柴达木地块的上地幔顶部具有相对高速特征,而在上地幔200km以下这两个地块间的差别逐渐减小.1654年天水地震和1879年武都地震都发生在扬子地块与青藏高原的碰撞交汇区,其震中下方上地幔显示为高低速转换结构.     

1976年唐山大地震的孕震环境:密集地震台阵观测得到的结果

1976-07-28唐山地区发生了震惊中外的7.8级大地震。为什么在华北古老克拉通内部的唐山地区能够发生如此的大地震一直是一个令人费解的问题。是否会在唐山地区再次出现同样的破坏性地震值得认真研究。利用流动地震观测台阵数据和接收函数反演方法,我们研究了唐山地区60 km深度范围内的三维地壳上地幔速度结构。结果表明:(1)由活动断裂切割的唐山断块与周围介质存在明显差异,围限唐山断块的断裂均为超壳的活动断裂;(2)唐山大震区中上地壳具有明显的非均匀壳内低速体;(3)该地区壳幔界面表现为明显的断块式隆升,与两侧相比,唐山菱形地块下方的上地幔顶部异常隆起的高度达到10 km左右,下伏的上地幔具有异常的非均匀结构;(4)唐山大震区可能有幔源物质较大规模的侵入,形成了中、上地壳内的低速体。由于较已往的研究结果有更高的空间分辨率,我们得到了一些以往尚未发现的有关唐山地区深部结构的异常特征;(5)首都圈地区内破坏性地震发生的地点绝非偶然,它们均与其相应的深部构造背景密切有关,这为强震发生地点的预测提供了可能。根据本文结果,我们认为,1976年唐山大地震的主因源于上地幔的垂向运动变形及壳幔之间物质及能量的交换,区域水平向应力场为次要作用。这有助于解释为什么能够在我国华北古老克拉通地区发生7级以上强震,在唐山地区再次发生7级以上大地震的可能性值得给予进一步的研究和关注。     

中国境内天山上地幔小尺度对流与造山作用

根据横跨天山的流动地震台阵记录的走时数据, 利用地震层析成像技术获得了库车-奎屯剖面660 km深度内上地幔P波速度结构. 由P波速度模型导出了相应的上地幔二维密度模 型. 利用有限单元方法联合网格-粒子方法, 研究了中国境内天山由于密度差异分布引起的上地幔小尺度对流. 给出了对流的基本形态特征, 并分析了上地幔对流与造山作用的相互影响关系. 结果表明: (1) 北天山—准噶尔盆地下方的上地幔中存在一个逆时针对流环, 对流环尺度为500 km左右; (2) 南天山—塔里木盆地一侧的上地幔中存在一个顺时针对流环, 但相对较弱; (3) 上地幔对流速度值在很大程度上受黏滞系数变化的影响, 但天山上地幔顶部对流速度值不应小于20 mm/a; (4) 塔里木盆地向北的水平推挤作用是天山新生代隆升的决定性因素, 同时, 上地幔小尺度对流对天山现今构造变形格局也起着重要作用; (5) 塔里木板块向北的运动对上地幔对流场的影响很小.     

浅析北斗应用系统在海南渔业信息建设中的应用

根据北斗系统的现有技术条件，充分利用成熟的信息处理手段，探讨北斗应用系统在海洋渔业发展建设中的关键技术及其实现。