珊溪水库地震与构造地震波谱时-频特征的对比研究

基于2002~2006年浙江地震台网记录的数字化波形资料,利用时-频分析方法,研究了珊溪水库诱发地震和构造地震波谱的时-频特征,发现两者的P波波谱存在差异:1)水库诱发地震P波主频低于构造地震;2)水库诱发地震P波频率成分比构造地震丰富,能量密度频谱随频率轴展布较宽,且出现多个能量较强的频率段;3)水库诱发地震P波能量强度的峰值在整个波列中出现时间较构造地震早且能量衰减快。     

基于频率拓展的表面多次波压制方法

为消除多次波预测中多余的子波效应,拓宽地震数据的频带,使预测多次波的走时、振幅和相位与实际的多次波相吻合,根据原始数据,衍生频率拓展道等多道地震数据;根据最小二乘原则,得到自适应滤波器和压制多次波后的有效波数据。结果表明:基于频率拓展的多次波压制方法,能够有效压制多次波,提高地震资料的信噪比;对水平层状模型测试结果频谱分析,自适应相减后得到的有效波同预测的多次波相比,频带拓宽,高频能量增强;通过对SMAART模型测试,基于频率拓展的多次波压制方法在噪声较强的情况下,也能得到很好的多次波压制效果。对于地质结构较复杂的情况,与抛物Radon变换滤波法相比,该方法多次波压制效果更好,并且不损害有效波能量。     

十字交叉排列面波压制方法及应用

考虑面波特征和古龙断陷区地震资料面波强的特点,提出面波压制的十字交叉排列法.该方法将地震资料中的炮点线和检波点线重排后,进行三维Fourier变换;在频率—波数域滤波后,进行Fourier反变换,实现面波压制处理.古龙断陷区地震资料面波压制结果表明,十字交叉排列面波压制方法具有很好的保幅性和实用性,能够为地震资料保幅处理提供较好的叠前预处理方法,可以在深层地震资料处理中推广应用.     

λ-f域抛物Radon变换多次波压制方法

空间截断效应降低传统Radon变换压制多次波的精度.提出一种λ-f域抛物Radon变换方法,与传统最小二乘抛物Radon变换相比,地震数据变换到λ-f域后表现为直线的映射,更便于设计滤波函数,从而有效消除传统多次波压制过程中空间截断效应的影响.另外,新变量λ的引入,使Radon正反算子成为与频率无关的矩阵,可大幅度提高计算效率.理论模型和实际海洋地震数据多次波压制试算结果表明,提出的λ-f域抛物Radon变换法可有效实现多次波压制,提高地震数据的信噪比,具有较强的实用性.     

2022-01-15汤加火山喷发的地震学震级测定北大核心CSCD

利用山东地震台网70个宽频带测震台站记录的2022-01-15南太平洋汤加海底火山喷发原始地震波形数据,从时间域和频率域分析火山喷发地震波的记录特征;通过对20 s优势周期面波进行研究,测算该次火山喷发的地震学震级与地震波辐射能量。研究表明,该次汤加火山喷发的面波震级为M_(S)5.674、能量震级为M_(e)5.704,与GFZ测定震级基本一致,比USGS测定震级小0.126。本次汤加火山爆炸指数(VEI)为5级,但本文测定的地震波辐射能量为8.155×10^(12)J,可见火山喷发以地震波辐射方式释放的能量较少,绝大部分能量以超过1000 km/h的冲击波形式穿过大气层迅速释放到空中,这种能量的快速释放会抑制深部应力的挤压与拉伸。     

基于波动理论和抛物Radon变换的表面多次波压制技术

基于波动方程预测相减的表面多次波压制方法,能较好压制地下复杂介质产生的多次波,但是计算成本较高。基于滤波方法压制多次波计算效率较高,对一次波和多次波存在较大时差的地震数据有明显的压制效果,但是对来自复杂介质的地震数据较难得到良好的压制效果。将波动方程法和滤波法结合,提出一种基于抛物Radon变换和波动方程的多次波压制方法。首先应用表面多次波形成机理预测多次波,利用对应的预测多次波的地震记录和初始地震记录,在抛物Radon域设计自适应滤波函数,通过准确滤波得到多次波在抛物Radon域的能量;然后通过反变换得到多次波,从原始地震记录中减去多次波,得到压制多次波后的结果。水平层状模型及SMAART模型测试结果表明,该方法结合两种方法的优点,能够准确、高效压制来自地下复杂介质产生的表面多次波。     

黄土塬斜坡带三维地震数据采集难点及对策

摘 要:针对黄土塬斜坡带地震采集的技术难点,分析了斜坡效应及其主要影响因素,提出了地震波激发“能量窗”及“波场矢量叠加”的理论,采用“优选点位法”微测井技术,准确建立了地震地质模型,应用谱元法对激发－接收参数进行了系统的数值模拟。在此基础上,开展点－线结合的野外试验,量化分析试验的效果,确定了以超多井、小药量、线性组合、聚能层激发为特点的激发参数,以多检波器、矩形面积、中等基距组合、高叠加次数为特点的接收参数,获得了高分辨率的三维地震资料。      

Ｓｅｉｓｌｅｔ变换在海上拖缆多震源采集数据分离中的应用

多震源混合采集技术具有减小勘探周期、降低勘探成本、提高数据质量等方面的优势，近年来，海上拖缆多震源混合采集技术得到了发展和应用，由于海上拖缆采集每炮对应的排列位置随着炮点同步移动，导致多震源混合数据的分离往往是一个欠定的反演问题，Ｓｅｉｓｌｅｔ变换对地震数据具有更好的压缩特性，并且能够压制随机噪声，为此，分析了不同数据域海上拖缆多震源采集数据特征，将 Ｓｅｉｓｌｅｔ变换应用于海上拖缆多震源混合数据的分离，海上拖缆多震源合成和实际采集的混合波场分离结果表明，该方法能够取得更好的分离效果。      

闽粤地区同震形变波记录序列分析及小波能谱特征

基于闽粤地区10个地形变观测点40个测项的同震形变波记录,利用幅度比值和小波多尺度分析方法探讨其记录特点和传导机制,结果显示：1)同震形变波到时与震中距有关,波形记录含有d-P和d-S波,不同测项记录到的波形不一致,同测项的记录波形相似;2)幅度比值在一定程度上与震源方位有关,进一步反映其各向异性;3)通过互（自）相关计算和滤波处理,能识别出同震形变波的频带为002~0.04 Hz,峰值在0025 Hz附近;4)智利地震形变波序列小波能谱为64~254 s,能量集中在128~254 s;而尼泊尔地震形变波在不同测点区别较大,潮州台为16~32 s,厦门台为16~64 s。研究结果表明,形变记录中除面波外还包含体波成分,不同震源机制引起的同震形变波传导机制可能存在各向异性,这对于开展地震引发应变等研究有借鉴意义。     

地震勘探全波形反演的应用与发展分析

本文首先对20世纪80年代发展起来的全波形反演应用及其在勘探地球物理领域的发展进行了分析;其次,面对定量化、精细化的地震勘探要求,提出了将地震勘探全波形反演与其他数据处理环节或处理技术相结合的研究设想,并展望了全波形反演的发展趋势;最后,论述了全波形反演研究中地震波场数值模拟、反演初始速度模型获取、目标函数形式选择、寻优算法启用及各向异性介质中的应用等关键问题,并总结了通过Laplace域的全波形反演获取反演初始速度模型、结合射线追踪并充分发挥并行计算之于波动方程方法来模拟地震波场的巨大优势,及灵活选用反演目标函数形式和寻优算法更新速度模型参数来加快全波形反演方法的实用化进程。     

高精度广义KL变换波场分离与去噪

KL滤波是地震数据处理中的重要方法技术之一。针对在实际中如何有效地应用好KL变换的问题，结合有效波同相轴高精度拾取和排齐技术，采用奇异值分解（SVD）算法实现KL变换，使KL滤波适用于任意弯曲同相轴的分离与消除。理论模型计算和实际资料处理试验表明，使甬该方法能有效地分离相干噪声和随机噪声，提高地震记录的信噪比，简单实用．便于操作。     

碎屑岩储层模型的转换波正演模拟波场特征及其分辨率分析

基于东部某油田物性资料,建立典型碎屑岩储层模型,利用弹性波方程交错网格有限差分正演,对模型进行多波正演模拟。
通过分析正演模拟的炮集记录、叠前偏移剖面以及振幅属性,总结了多波正演的波场特点及其分辨率。分析结果表明,转换波反射频
率比纵波的低,能量随偏移距的增加而增加,零偏移距处能量为零;转换波在地层厚度为λ/4处发生调谐效应,调谐振幅属性变化规
律与纵波变化规律一致;综合考虑该区纵波、转换波波速以及不同深度的纵波、转换波中心频率变化,可知在该区目标层深度转换波
分辨率低于纵波分辨率。      

Morlet变换在超宽频带地震仪记录的汶川震前数据处理中的应用

介绍了Morlet小波变换的基本原理和特点。以汶川地震为例，选取2008-05-01~05-12的地震计垂直向连续波形数据，使用Morlet小波变换得到时频谱。分析时频图发现，在汶川地震发生前的数十h中，0.1~0.2 Hz频段的能量幅值出现明显的递增趋势，是非台风扰动引起的频率异常，与已有论文结论符合良好，证明Morlet小波变换是分析震前频率特征的有效方法。     

江西地区地震动非弹性衰减和场地响应特征研究

以2009年以来江西省地震台网记录到的86个M_L2.5以上地震事件的波形资料为研究对象,通过遗传算法获得介质的品质因子和台站场地响应,并在此基础上计算地震的震源谱参数。结果显示,江西地区Q值与频率f的关系式为Q（f）=323.1f^{0.505 9};大部分台站的场地响应在频率域表现为平坦型,整体上与台站基岩性质相符;近震级M_L与地震矩M₀在单对数坐标下呈线性关系,地震矩与拐角频率呈负相关特征,应力降与地震矩之间没有显著相关性;震源半径与应力降存在显著的双对数关系。     

基于S变换的瑞利面波频散分析

提取面波频散曲线是面波资料处理中最关键的一步.由于时频分析方法的局限,提出了利用S变换进行瑞利面波频散分析的方法,并给出了具体算法.该方法在时间-频率域中计算相邻两道面波记录同一频率的时间差,再利用道间距除以该时间差来得到该频率对应的相速度,这样避免了在道间距较大的情况下,传统的相位谱法可能造成相位差的缺陷.通过理论模型和实际资料对该方法进行了验证.结果表明,该方法能够提高瑞利面波频散曲线的提取精度,而且算法简单,具有一定的实用性.     

黑龙江地区上地幔各向异性研究

采用黑龙江数字地震台网2012-01～2014-07间23个宽频带地震台站记录的远震SKS波形资料,利用Splitlab 软件对黑龙江地区进行剪切波分裂的分析研究,使用最小能量法、旋转相关法和最小特征值法计算所有观测台站的SKS快波偏振方向和快、慢波的延迟时间,最终得到黑龙江省上地幔各向异性图像。结果显示,黑龙江地区上地幔地震各向异性比较明显,快波偏振方向与主张应力方向基本一致,与GPS得出的速度场方向相符,说明该地区的壳幔耦合可能存在垂直连贯性。     

盖州青石岭震群重复地震尾波Qc值变化特征

选取营口台和岫岩台记录的2012~2017年辽宁盖州青石岭M_L≥2.5地震波形资料,利用波形互相关方法拾取青石岭震群重复地震序列,并采用Sato单次散射模型计算单台记录的重复地震序列尾波Q_c值对频率的依赖关系及随时间的变化。结合盖州青石岭震群地区的地震活动特征可知,Q_c值变化所反映的区域应力变化与中强地震具有一定相关性;盖州地区在M_L≥4.0地震后Q_c值波动明显,存在局部调整并恢复的特点,局部应力积累不足以孕育更大的地震,而是以多次中等地震的形式释放局部积累的能量。     

2019-12-26湖北应城M4.9地震强震动特征分析

对2019-12-26湖北应城M4.9地震的强震加速度记录进行基线校正和滤波处理,分析地震动场的分布特征。结果表明,此次地震记录的峰值加速度（PGA）随距离衰减较快;场地效应、上盘效应和破裂方向性效应共同导致地震动产生巨大差异。对应城地震台的强震记录进行旋转,其顺断层方向记录与垂直断层向和垂直向记录差异明显,存在一段周期约为0.3 s的类简谐强震动,与以走滑为主的震源机制相吻合。     

基于小波特征和神经网络的天然地震与人工爆破自动识别

采用福建地区天然地震和人工爆破事件波形记录，通过一维离散小波变换（DWT）及4层小波包变换（WPT）对信号进行分解，提取出用于识别的4种波形小波特征：小波能量比特征、小波包能量比特征、小波包香农熵特征及小波包对数能量熵，此外还提取出P/S震相振幅比；采用BP神经网络对4种小波特征及分别加入P/S震相振幅比的组合特征进行识别效果检验，结果表明，单小波判据小波能量比特征识别效果好；双判据组合P/S震相振幅比和小波包对数能量熵的组合识别效果最好，可考虑作为实际天然地震与人工爆破在线自动识别系统的识别判据。     

使用支持向量机识别地震类型的影响因素分析

对山东地区2006~2017年3种地震事件--天然地震、爆破及塌陷的波形记录进行小波变换,对提取出的香农熵特征采用支持向量机LIBSVM方法进行分类识别,并设计一系列实验研究影响最终分类效果的因素。结果表明,5种影响因素--信号窗长度、小波分解方式、小波基类型、向量机算法类型、向量机核函数类型均对地震类型的分类识别结果产生一定的影响;识别率最高的3组处理方式均采用了2 000 s信号窗长度+db7小波基+υ-SVC算法的组合方式。所得的识别率较高的几种影响因素组合,在未来可应用于地震类型的实时识别,进一步提高地震类型的识别率和触发准确率。