剪切波分裂中的快、慢波识别方法

本文改进了两种震相识别方法。用最大特征值法可识别快、慢剪切波到时,用波形识别算子法可确定快剪切波的到时及其序列。将两种方法结合起来,可比较和鉴别地震波的到时及其类型,从而可以避免将转换P波误认为快剪切波,达到最终确定剪切波到时及其序列的目的。以1985年禄劝Ms6.1地震为例,用上述方法进行了剪切波分裂的研究,结果与张性扩容各向异性假说是一致的。     

横向各向同性地层井孔伪瑞利波、弯曲波、螺旋波的低频极限速度

在多极源声测井中, 低频弯曲波或螺旋波被广泛应用于测量地层的横波速度. 前人的研究已证明在各向同性地层中井孔伪瑞利波、弯曲波及螺旋波的低频极限速度都等于地层横波速度. 大量的数值计算结果似乎表明此结论在横向各向同性（TI）地层情况下也能成立, 但缺乏理论证明. 本文在井孔平行于TI弹性地层对称轴的模型下, 考察了井内声源激发的波在流-固边界上的反射和透射情况, 阐述了非泄漏模式导波产生的必要条件并讨论了其速度上限值. 我们发现在各向异性参数满足一定条件的TI地层中,导波速度的低频极限值小于沿井孔方向上的横波速度. 通过对井孔导波的频散分析以及对时域全波列的数值模拟, 论证了在这类地层中进行多极源声测井时,不可能根据弯曲波或螺旋波的到达时间准确获取地层的横波速度值.     

新疆北天山西段P波、S波衰减特征研究

针对2008年1月至2011年12月发生在新疆北天山西段的地震事件,分析新源、温泉和巴仑台3个台站记录到的波形资料,采取尾波归一化方法,计算出P、S波在多个频率点f的衰减系数Q-1P、Q-1S,利用最小二乘法拟合得到Q-1=Q0fη,得到该区域P波衰减系数和S波衰减系数。结果表明,新源研究区关于P波、S波的衰减指数为-0.982 1、-0.978 1;温泉研究区关于P波、S波的衰减指数为-1.049 7、-0.773 8;巴仑台研究区关于P波、S波的衰减指数为-0.870 0、-0.956 8。与中国云南地区关于P波、S波衰减结果相比,新源研究区域与云南的禄劝研究区地下介质衰减特性比较一致。     

云南地区与频率有关的P波、S波衰减研究

分析了1999年7月-2003年12月发生在云南地区的由22个地震台站记录的5 668次地震的数字波形资料,从中选出保山、永胜、丽江、鹤庆、易门和禄劝6个地震台用扩展的尾波归一化方法计算P波、S波衰减.这6个台分为3个研究区,保山研究区内Qs-1=0.008 67f-0.86,Qp-1=0.011 55f-0.93;永胜、丽江、鹤庆研究区内Qs-1=0.018 24f-0.94,Qp-1=0.022 88f-0.92;禄劝、易门研究区内Qs-1=0.016 47f-0.91,Qp-1=0.028 26f-0.97.后两个研究区衰减结果比较接近,而保山研究区的P波、S波衰减明显偏低.与全球其它地区采用同样方法得到的P波、S波衰减结果对比表明:永胜、丽江、鹤庆研究区和禄劝、易门研究区的P波、S波衰减与日本关东地区很接近,而远高于韩国东南部的衰减;保山研究区的P波、S波衰减略高于韩国东南部.此外,还与前人得到的云南地区的尾波衰减结果进行了对比分析,发现本文S波衰减值与他们的结果相近.     

数字化地震波形小震P波、S波震源参数对比研究

目前,国内在数字化地震资料研究中,主要应用S波来研究震源参数.本文充分利用上海数字化地震台阵(网)资料,挑选了2001～2004年华东地区发生的9个震级较大的中等地震,反演了单台记录垂直向纵波(P波)及横波(S波)震源谱,根据布龙模型,计算每个地震单台拐角频率、地震矩、零频谱值、震源破裂半径、应力降等参数,求出每个地震的震源参数(P波和S波)的平均值及其均方差,对这9个地震的P波、S波震源参数计算结果进行了对比.结果表明,用P波段数据进行震源参数的研究是可行的.     

欧亚大陆均匀震级系统P波、S波和L波的最佳振幅曲线

地震的 P 波,S 波和 L 波的振幅曲线,在欧亚大陆均匀震级系统的基础上建立,最佳振幅曲线和地震实践中应用的标准振幅曲线作了比较。     

接收函数、面波与重力联合约束地壳厚度与波速比

提出了一种利用接收函数、面波频散和重力数据联合约束地壳厚度、v_P/v_S和平均P波速度的改进方法,并基于两种地壳模型对改进后的方法进行了验证。结果显示,改进后的方法不仅可以提高地壳厚度和波速比的估计精度,还能更准确地估计地壳平均P波速度。在此基础上,将该方法应用于华南地区两个固定台站的地壳结构分析,相关结果也证实了该方法在确定地壳结构中的可行性。     

用非零偏P-P波、P-SV波计算零偏反射纵波和反射横波

在P波入射时，地下同一点产生的P－P波和P－SV波是入射角、界面两侧纵波阻抗和横波阻抗的函数，并且两者是耦合的，利用两者的耦合关系，可在入射角确定的情况下提取与入射角无关的零偏纵波反射和横波反射。本文对此进行了理论推导，结论是：入射角为零时的纵波与横波反射系数是入射角非零时的P－P波和P－SV波反射系数的线性组合，进而可以用纵横波反射计算泊松比反向间面，把多波叠前反演转化到叠后进行，通过合成弹性参数反射剖面，直接反演弹性参数等，都很有意义。文中给出了有关的公式，并对具体实现方法和有关参数进行了讨论。     

前驱波还是气压波？

2008年5月11日山西地区5个地壳形变测项出现了“前驱波”信号,但其真实的物理本质究竟如何？至今仍存质疑。为此,本文借助小波变换方法对气压波和“前驱波”在不同时、频段的相关性进行了较系统的诊断;还进一步将“前驱波”与短时气压波动所产生的干扰实例进行了对比分析。结果表明:这些“前驱波”并非源自汶川M_S8.0地震,而主要是气压波动的干扰。      

澜沧-耿马地震的S波偏振及P波、S波振幅比变化特征

本文通过对1988年11月6日云南澜沧—耿马7.6、7.2级强震前后S波偏振角的计算分析,得出大震前后U_(SH)—U_(Sv)的三种变化图象,这种变化特征具有一定的前兆意义。与此同时,对大震前后P波、S波最大振幅比值分布和S波偏振角空间分布进行研究发现,北纬22°57′可能是澜沧、耿马两地区介质差异的分界面。     

基于纵横波解耦的三维弹性波逆时偏移

纵横波波场分离是弹性波偏移方法的必要条件,通过纵横波成像的差异可以获取更多地下介质的信息.目前所用的纵横波波场分离方法多采用Helmholtz分解,这样得到的波场不仅物理意义发生了变化,振幅和相位也会发生改变.本文采用纵横波解耦的弹性波方程,将其应用于三维介质,对比分析了纵横波解耦方法相对传统Helmholtz分解方法在相位、振幅上的优势.将该解耦的波场分离方法应用于弹性波逆时偏移,能得到相位、振幅和物理意义不受改变的偏移结果.但是该解耦方法分离得到的纵横波波场均为矢量场,将该波场分离方法用于弹性波逆时偏移,还需要解决矢量场如何得到标量成像结果的问题.本文引入了Poynting矢量,通过Poynting矢量对矢量波场进行标量化,这样就能得到保振幅、相位,且无极性反转的标量PP和PS成像结果.同时针对S波Poynting矢量求取不准确的问题,采用拟S波应力场和S波速度场得到了更加准确的S波Poynting矢量.理论计算证明了本文采用的3D波场解耦的矢量波场分离方法的正确性和引入Poynting矢量对矢量波场进行标量成像的有效性.     

逆子波域消除多次波方法研究

SRMA（与表面相关多次波的衰减）算法包含预测和相减两步.相减算法中,当多次波与反射波同相轴相交时,如何有效减去多次波、保留反射波,是面临的主要问题.通过分析非正交性对滤波器（逆子波）的影响,可以证明：逆子波因非正交性产生的误差呈近似的高斯分布.在此基础上,本文提出了在逆子波域（单道自适应相减滤波的滤波算子的集合）,利用其误差的概率分布特征,对逆子波进行估计,用逆子波的估计对逆子波进行校正来消除多次波的方法.其步骤为：首先用SRMA方法预测出表面多次波,并对每一单炮进行单道自适应相减,得到逆子波,形成逆子波域;其次,在逆子波域采用中值滤波,提取接近真实逆子波的逆子波估计;第三,在逆子波域用逆子波估计对畸变的逆子波进行校正;最后采用校正后的逆子波来衰减多次波.通过简单模型和SMARRT模型的测试,该方法不仅能够有效减去多次波,而且在相交的区域,能够保持反射波同相轴的连续性并恢复其正确的振幅.     

利用地方震、远震走时和面波数据联合反演日本俯冲带P波和S波层析成像

通过最新收集的大量高质量的地方震和远震事件的到时数据进行联合反演,我们确定了日本俯冲带约700km深度的P波和S波速度层析成像。我们还使用远震瑞利波的振幅和相速度,确定了日本及其附近海域下方20~150s周期基阶瑞利波的二维相速度图像。研究区精细三维S波层析成像可通过地方震和远震事件的S波到时,及瑞利波相速度数据进行联合反演得到。我们的反演结果揭示:一维原始速度模型中,俯冲太平洋板块和菲律宾海板块呈现明显的高速区。在板块上方的地幔楔和太平洋板块下方的地幔中存在显著的低速异常。俯冲板块和周围地幔之间速度有明显的差异,表明温度、水含量和/或部分熔融程度有显著的横向变化。地幔楔低速异常是由板块脱水作用和地幔楔拐角流造成。在日本东北太平洋板块下方显示片状的低速区,这可能反映了地幔深部热上涌以及地幔柱软流圈的俯冲作用。我们的结果表明不同的地震数据联合反演,对于得到地壳和地幔可靠的层析成像图像是非常有效和重要的。     

基于波场分离的弹性波逆时偏移

尽管叠前逆时偏移成像精度高,但仅针对单一纵波的成像也可能形成地下介质成像盲区,由于基于弹性波方程的逆时偏移成像可形成多波模式的成像数据,因此弹性波逆时偏移成像可提供更为丰富的地下构造信息.本文依据各向同性介质的一阶速度-应力方程组构建震源和检波点矢量波场,再利用Helmholtz分解提取纯纵波和纯横波波场,使用震源归一化的互相关成像条件获得纯波成像,避免了直接使用坐标分量成像而引起的纵横波串扰问题.针对转换波成像的极性反转问题,文中提出一种共炮域极性校正方法.为有效节约存储成本,也提出一种适用于弹性波逆时偏移的震源波场逆时重建方法,在震源波场正传过程中,仅保存PML边界内若干层的速度分量波场,进而逆时重建出所有分量的震源波场.本文分别对地堑模型和Marmousi2模型进行了弹性波逆时偏移成像测试,结果表明:所提出的共炮域极性校正方法正确有效,基于波场分离的弹性波逆时偏移成像的纯波数据能够对复杂地下构造准确成像.     

利用标准时频变换方法在强噪声环境下无偏拾取地震P波、S波到时

地震P波、S波到时是精确分析地震水平位置、深度与速度结构等的重要参数,如何准确拾取P波和S波到时是地震学的一项重要的基础工作.大数据量与强噪声环境给地震到时的自动拾取带来了很大挑战.在频率域中可将信号与噪声分离,但会造成震相的偏移.针对上述问题,本文在STA/LTA、AIC方法的基础上,引入了标准时频变换(Normal...     

利用地震转换波探测地壳、上地幔构造试验

当远震P波通过地壳、上地幔速度界面时,除了产生PP型折射波,还产生一系列PS型折射转换波.这类PS转换波与初至P波之间的到时差,可用来测定地壳和上地幔构造. 转换波测深法首次在我国鲁西夏津-高唐地区和京西昌平地区进行了试验应用.为提高转换波震相识别的可靠性,在资料分析中采用了偏振滤波的方法.工作结果获得了上述两地区的深部构造剖面图.转换波法与爆炸地震法所得结果之对比表明,两者的主要界面基本一致.     

基于单程波偏移算子的地表相关多次波成像

在常规地震资料处理中,多次反射波被视为噪声并从地震数据中去除,以免在之后的地震资料解释中造成误解.而事实上,多次波也是地震信号,是照明波场的一部分,能够对地下构造成像的精度做出贡献.本文分析了多次波在传统单程波叠前深度偏移中产生构造假象的机制和表现,为实现基于单程波偏移算子的多次波成像,修改了单程波叠前深度偏移的边界条件,即将输入的震源波场用包含多次波的记录来替代,输入的记录波场用预测出的表层相关多次波来替代,实现了基于单程波偏移算子的地表相关多次波成像,并从理论上给出了其成像依据.通过基于二范式最小能量差原则求取的匹配因子,将多次波成像结果与一次波成像结果进行匹配叠加,应用多次波成像来弥补一次波成像的不足.简单模型验证了基于单程波偏移算子的多次波成像方法的有效性,最后对Sigsbee2B模型进行了一次波与多次波联合成像试算,盐边界高陡构造成像质量得到了明显改善.     

地震短期、临震前兆应力波的观测与分析

本文介绍了对地应力波进行的初步观测及结果。提出了孕震源是一扰动源,也就是震前应力波产生的波源。提出了前兆应力波扩展与传播的最小阻力或最短波径定律以及地应力波传播的概略模型。     

矢量分离纵横波场的弹性波逆时偏移

弹性波逆时偏移是当前多分量地震资料相对准确的偏移算法,它能够形成多波模式的成像剖面,从而减少纵波勘探的多解性.本文首先依据各向同性介质中矢量分离纵横波场的速度-应力方程组,利用高阶交错网格有限差分数值方法求解弹性波方程,进而构建矢量的纵波和横波波场,不同于散度和旋度算子分离纵横波场的传统方法,文中提出的矢量分离纵横波场方法保持了原始波场的振幅和相位特征.文中也提出将震源归一化的内积成像条件应用于分离后的纯纵波和横波矢量场,由此得到的转换波成像避免了传统弹性波成像方法中出现的极性反转.水平层状和复杂构造模型测试表明,文中提出的基于矢量分离纵横波场的弹性波逆时偏移方法成像精度高,转换波成像PS和SP极性无反转,所形成的多种模式纯波剖面能够准确地对复杂地下构造成像.     

长三角地区地壳P、S波速度结构模型

利用收集到的瑞利面波频散数据和区域地震台网震相数据,分别采用面波层析成像和体波走时层析成像方法反演了长三角地区地下速度异常分布,获得了该区地壳P、S波速度结构.层析成像结果表明速度结构模型与地质特征基本一致,研究结果有助于深化对长三角地区沉积盖层和地壳结构以及强震活动深部环境的认识.