基于小波变换的地震动时-频特性研究

对美国NGA,采用一维连续小波变换得到每条记录的小波功率谱。研究了任意时间处小波功率谱最大值所对应的主频率。结果表明:地震动主频率随时间的增大逐渐减小;竖向地震动分量的主频比水平向分量随时间减小更快。分别采用线性函数模型、指数函数模型和指数三角函数模型,分析了场地、震级和震中距对主频变化曲线的影响,拟合了主频率在不同场地条件、震级、震中距等情况下随时间的变化曲线,水平向和竖向的频率时变曲线整体上都是随时间递减的,且竖向衰减得更快些,大多数情况下竖向记录的高频成分比水平向记录的相应成分要多。     

基于广义S变换的地震与爆炸识别

为了进一步增强区域台网针对天然地震与人工爆炸事件的识别能力,本文利用广义S变换方法,围绕河北省三河采石场爆破以及周边发生的天然地震波形记录展开研究,总结了地震与爆炸在时频谱中的差异性。研究表明,受震源机制影响,天然地震时频谱图的频带范围更宽广,能量团分布也较人工爆炸更为复杂。为了消除震中距的影响,本文将求解的时频谱保存为规格相等的灰度图像,通过滑动窗口计算图像的灰度一致性得到新的识别判据?谱图二阶矩。对三河地区已知地震和爆炸事件的测试表明,谱图二阶矩对单台波形记录的识别率可达91%,多台求取平均值之后的识别率超过98%,因此该判据具有良好的应用前景。      

一种实时场地校正方法及其在地震预警中的应用

地震预警系统对地震数据处理的实时性要求极高。其系统数据来源除布设在基岩的测震台站外,还有大量非基岩场地的强震动台站和地震烈度仪台站,其场地影响不容忽视。为了考虑震级估算和地震动场预测中的场地影响,需实时对各种场地条件下的地震波形进行校正。目前处理一般使用某个标量来表征场地放大效应。本文采用一种实时的场地校正方法,首先计算目标场地与参考场地的谱比,然后通过最小二乘法、双线性变换将谱比转化为因果递归的无限冲激响应(IIR)时域滤波器,之后可以应用该滤波器进行实时场地校正。该方法考虑了场地放大系数的频率依赖性,相比于标量校正提高了准确度。应用我国四川和日本部分强震动台站记录,验证并讨论了这种实时场地校正方法在地震预警中的应用效果。     

基于离散Gabor变换的磁暴识别

以地磁秒数据为研究对象,通过离散Gabor变换将时域的地磁数据转换至二维时频面,提取Gabor变换谱图的均值和方差作为特征值,使用支持向量机实现地磁正常数据与磁暴干扰数据的自动分类识别。对5个地磁台的200组地磁秒数据进行计算分析,结果表明该方法对测试样本数据的识别率可达94%。     

内蒙古中部地区形变主要干扰的时频响应特征分析

选取内蒙古中部地区定点形变观测中受降雨、大风、气压干扰以及地震波影响的典型数据,采用S变换,对该地区形变观测中4种主要干扰信号的时频响应特征进行分析与讨论。结果表明,降雨干扰信号主要集中在低频区域,其频率先增大后减小;大风干扰表现为高频干扰,时频域能量强度与观测曲线受风扰影响变化幅度成正比;气压干扰信号的优势频率分布在低频段内;地震波影响中VP垂直摆和JCZ地震计记录的地震波较为相似,但VP垂直摆记录的频段要小于JCZ地震计,VP垂直摆主要记录的是低频地震波。同时,4种干扰信号的时频域频率、能量强度随时间变化趋势与数据观测时域变化趋势相一致。     

基于压缩感知的小波阈值和CEEMD联合去噪方法

中深层地质条件复杂,地震资料品质差,主要表现为:地震资料信噪比低、有效信号弱.如何在去噪的同时有效保留弱有效信号,获取高信噪比的地震数据成为地震数据处理的关键问题.传统小波阈值与互补集合经验模态分解(CEEMD)联合去噪方法相比单一方法可以获取更高品质的地震数据.基于压缩感知理论的去噪方法利用地震数据在变换域中的稀疏特性,通过设定稀疏基矩阵和测量矩阵,可以将地震数据去噪问题转化成求解最优化问题,通过最优解重构原始信号,实现对地震资料的去噪处理.该方法能够在有效衰减随机噪声的同时最大限度的保留有效信号.本文基于压缩感知理论开展小波阈值去噪方法研究,并在此基础上结合CEEMD方法对含噪较多的固有模态分量进行有针对性的随机噪声压制.通过对含噪数据开展不同方法的去噪结果对比可见,本文方法可以在保证高信噪比的基础上更为有效的保留弱有效信号,数值试算验证了该方法对弱有效信号地震数据去噪具有显著优势.     

基于压缩感知理论的地震数据降噪方法

针对地震数据在采集处理过程中存在的随机噪声,本文从压缩感知的角度,给出了一种地震数据降噪方法.其基本思路是:首先对含有随机噪声地震数据通过离散余弦变换进行稀疏表示,然后选取随机高斯矩阵为测量矩阵,并计算出传感矩阵,在地震数据重构阶段,采用正交匹配追踪算法对地震数据进行重构;通过实验方法对比,本文方法的降噪效果在峰值信噪比、信噪比、均方误差指标上均优于对比方法,证明了本文方法对地震数据中的随机非平稳噪声有较好的压制效果,提高了地震数据的信噪比.     

等效交错网格高阶有限差分法标量波波场模拟

密度是岩石物性参数的重要组成部分,对储层评价、岩性解释和油藏描述等具有重要作用,因此实际资料应用中含有密度信息的正演建模是必不可少的.作为波动理论应用最广泛的实际资料处理技术——逆时偏移和全波形反演,其正演建模一般采用经典的二阶标量波方程,且不考虑密度的空间变化.而一阶速度-应力控制方程通过交错网格有限差分法正演建模,其精度高于二阶方程但计算成本过高,很难在三维实际资料中应用.鉴于非均质正演在实际生产中所面临问题,本文借用交错网格的思想,充分考虑空间变密度对波传播的影响,提出基于等效交错网格的高阶有限差分建模算法并应用于非均质标量波正演.文中从数学上给出一阶方程和二阶方程在数值模拟中的等价性证明,并分析了震源、边界条件和稳定性.最后通过简单层状模型进行数值测试,对比不同建模方法以验证新方法的准确性;利用高阶精度方法对Sigsbee-2a模型做正演测试,验证本文提出的方法的稳定性.     

利用小波变换分析羽状流模型偏移剖面尺度特征

甲烷羽状流备受关注,因为它提供了海底气体泄漏的直接证据,并对天然气水合物的勘探识别起到了间接指示作用.为了深入研究冷泉活动区气泡羽状流的地震响应特征,根据已有研究成果——羽状流模型及其偏移剖面,通过小波变换分析单道地震信号及羽状流剖面的尺度特征.单道地震信号最佳尺度曲线表明不同气含量对应的最佳尺度也不同,而且不同深度的最佳尺度变化较大;羽状流最佳尺度剖面能够显示出羽状流区域的轮廓.羽状流气含量与最佳尺度呈较好的正相关,即随着气含量的增加,最佳尺度也随之增大.     

用振幅变化长短时均值比实现地震预警中P波自动拾取

地震P波到时快速、精确的自动拾取是实现地震预警的第一步。目前所有P波到时自动拾取方法都需要进一步考虑准确度问题,据此提出用"振幅变化"长短时均值比方法拾取P波到时。研究表明引进"振幅变化"作为特征函数比Allen的"振幅变化平方"更能突显地震波初至时刻的"突变"特征。对"振幅变化长短时均值比方法"和Allen的"振幅变化平方长短时均值比方法"进行抗噪音分析表明,在较低信噪比条件下,Allen方法不能准确地拾取P波到时,而新方法依然能较准确地拾取P波到时。在较高信噪比条件下,用两种方法拾取的160个M_S4.0地震P波到时,其误差大都在允许范围内,但新方法拾取的准确度略高。新方法可以看做是对长短时均值比方法的拓展。