同步挤压广义S变换在南海油气识别中的应用

受窗函数的影响,短时傅里叶变换、小波变换、S变换等传统时频分析方法的分辨率不高.为了提高信号时频变换分辨率,发展了同步挤压小波变换算法;该算法将信号的时频谱值向其中心频率位置挤压,改善了信号时频分析结果,使得信号时频分辨率得到了很大的提升.但由于传统算法时窗固定,处理缺乏灵活性;因此,本文对S变换窗函数进行扩展,提出了同步挤压广义S变换.通过自适应窗函数挤压信号在S域的时频谱值,提高了算法的灵活性和时频分析聚焦能力.运用同步挤压广义S变换对南海某工区实际地震数据进行分频处理,结果显示含油气层能量随着频率的增加而逐渐衰减.因此,使用同步挤压广义S变换对地震数据进行分析处理,不仅可以对储层含油气性进行准确的检测,还可以对地层进行精确的标定.     

联合广义S变换和压缩感知提高地震资料分辨率

提高地震资料的分辨率是获得高分辨率地震剖面的基础,常规的处理方法对于提高地震资料分辨率虽有一定的效果,但却十分有限,难以达到资料解释的要求.对此本文将广义S变换理论与压缩感知理论相结合,提出了一种新的叠后地震资料处理方法.实现步骤为:将叠后地震数据进行广义S变换;利用压缩感知理论构建L1/2范数最优化模型;半阈值迭代法求得最优稀疏解得到更精确的时频域信息;能量重新分配补偿时频域弱信号能量;广义S反变换得到处理后的地震信号.理论模型和实际资料处理表明,与常规广义S变换提高分辨率处理相比,本文方法处理结果具有更好的效果,地震资料分辨率有了明显提高.     

S变换在常规测井信息处理中的应用综述

为提高常规测井信息解决地质问题的能力,本文梳理了S变换的基本原理及性质,归纳总结了现今S变换用于地震资料处理的研究成果,分析展望了S变换用于测井信息处理的可能方面.研究表明,S变换是一种具有高时频分辨率的信号处理方法.在地震资料处理中,S变换可用于地震波的初至识别、地震信号滤波去噪、基于地震信息的沉积旋回划分、碳酸盐岩...     

局域波分解及其在地震信号时频分析中的应用

讨论非平稳、非线性局域波信号分解方法的物理意义和自适应性,数值实现三种局域波分解算法,研究局域波分解在地震信号时频分析中的应用方法和实践.通过Hilbert变换得到地震信号时频分布特征,与Fourier变换谱物理意义不同.采用一定的计算策略,改善瞬时频率计算精度,突出瞬时频率属性的物理意义;将局域波分解同Wigner分布结合,计算地震信号基于局域波分解的Wigner分布,抑制交叉项.模型和实际地震数据试算结果表明实现算法的正确性和有效性.研究表明：在地震信号时频分析中采用局域波分解求Hilbert谱具有自适应性,时域和频域的分辨率也较高;而基于局域波分解的地震信号Wigner分布,保留Wigner分布优良特性,抑制交叉项方法简单易实现;局域波分解方法适合处理地震信号,值得进一步研究和推广.     

二阶同步挤压S变换及其在地震谱分解中的应用

同步挤压S变换是一种处理非平稳信号的新时频变换方法,通过"挤压"信号的S变换结果,可以得到高分辨率的时频谱.但是,当信号的相位随时间呈非线性变化时,同步挤压S变换计算出的瞬时频率会出现误差,并造成时频谱分辨率降低.为了改善其对于该类信号的时频分析效果,我们利用时间和频率的二阶偏导数对瞬时频率计算式进行修正,提出二阶同步挤压S变换.合成信号处理结果表明,二阶同步挤压S变换的分辨率不但明显高于常用时频变换,在信号瞬时频率随时间呈二次或正弦变化的情况下,其时频挤压效果也好于同步挤压S变换.我们将二阶同步挤压S变换应用到天然气地震勘探资料的谱分解当中,结果表明,二阶同步挤压S变换可以很好地检测到与天然气相关的谱异常.因此,二阶同步挤压变换对于地震解释是一种很有潜力的方法.     

基于褶积原理和改进广义S变换的震动波衰减补偿试验研究

为了恢复震动波能量在传播过程中产生的衰减损耗,提出基于褶积原理求取品质因子Q的方法与改进广义S变换相结合的反Q滤波法。通过震动波衰减补偿模型试验,对试验数据进行改进广义S变换的时频特性分析,得出了信号的能量分布情况以及时间频率对应关系;采用基于褶积原理求取品质因子的方法,得到时变Q值;对试验数据进行反Q滤波处理,使震动波能量得到了补偿。结果表明本文提出的反Q滤波法提高了对震动波能量衰减补偿的效果,拓宽了地震资料的频带,提高了地震资料分辨率,有利于进行高分辨率地震勘探、深部信号增强和油气藏预测工作的开展。     

基于广义S变换的地震资料振幅谱补偿和相位谱校正方法研究（英文）

由于地下介质对地震波振幅的影响和地震波频散因素,地震波振幅和相位随时间、空间及频率的变化而发生改变,本文提出一种基于广义S变换的振幅谱补偿和相位谱校正新方法。该方法在S域中分为振幅谱补偿和相位谱校正两个步骤进行处理:振幅谱补偿是在地震记录可靠频带范围内恢复反射系数的振幅谱,其具体实现是在S域中利用谱模拟技术来拟合时变子波振幅谱,从而补偿由地层吸收所引起的振幅衰减;相位谱校正是消除子波剩余相位的影响,其具体实现是在S域中利用相位扫描来拾取随时间、空间和频率而变化的相位校正量,并由Parsimony准则来进行最佳相位判别。本文方法不需要直接求取Q值,能够适用于变Q值情况。理论模型和实际资料处理表明,该方法不仅能恢复地层反射系数的振幅谱,还可以有效消除子波剩余相位的影响,使子波接近或达到零相位,从而提高地震资料分辨率。     

带有低通滤波的广义S变换在探地雷达层位识别中的应用

同其它分析时变信号的方法相比,S变换具有独特的优点.但是S变换的时窗函数形式固定,限制了其在信号分析和处理中的应用.针对这一问题,本文对S变换进行了改进和推广,在S变换的时窗函数中引入调节参数,同时将低通滤波函数融入其中,提出了一种广义S变换.通过调节参数大小来调节广义S变换的时间分辨率和频率分辨率,从而使该变换具有很强的适用性.本文使用该变换对实测探地雷达数据进行了层位识别研究,取得了较好的效果.     

基于S变换的软阈值滤波在深地震反射数据处理中的应用

深地震反射原始单炮数据是非平稳的弱能量反射信号,信噪比较低.如何提高信噪比一直是深地震反射数据前处理中的一大难题.S变换是一种适用于分析非平稳信号的时频变换方法.同其他分析时变信号的方法相比,S变换的基本小波不必满足小波在时间域均值为零的容许性条件,它的时频分辨率与分析信号的频率有关,且其在时间域的积分可以得到傅里叶频谱,其反变换也简单.因此,S变换容易表示深地震反射信号复杂的时频特性.本文在S变换的基础上,利用软阈值滤波方法对深地震反射数据进行处理,实验结果表明,该方法有效地提高了信噪比,压制了有效频带范围内的混频干扰,突出了弱反射信号,使得波组信息更加丰富,有利于连续追踪有效反射波组和识别薄地层,特别是提高了深部Moho界面反射层位的分辨率,为深地震反射剖面后续处理和准确解释奠定了基础.     

2012年印尼8.6级地震应变地震波的Hilbert-Huang时频分析

本文通过端点效应压制的Hilbert-Huang变换, 对大同及沁源台布置的四分量钻孔应变仪记录的印尼8.6级地震激发的应变地震波形进行时频分析, 结果显示印尼8.6级地震的主震和8.2级余震的应变地震波序列各个震相具有不同的时频特征: ① 地震波到达之前的所谓“环境噪声”部分, 瞬时频率低, 瞬时振幅小; ② P波初至时, 高频成分突然增加, 振幅也随即增强; ③ S波到达时, 频率有所降低而振幅剧烈上升; ④ 面波到达时, 振幅进一步剧烈上升达到整个序列的极大值; ⑤ 尾波部分振幅逐渐降低, 但与噪声部分相比频率依然偏高, 振幅依然偏大。 本文也将应变地震波与地震仪记录的地震波进行对比, 虽然应变地震波与地震波波形和Fourier谱具有极高的相关系数, 但从Hilbert-Huang变换得到的边际谱上看, 应变地震波与地震波有显著的区别, 应变地震波比地震波记录的低频成分相对更多。 通过Hilbert谱, 有助于更好地了解非平稳信号的局部特征, 对于突变信号的地震波, Hilbert-Huang变换是一个较好的时频分析工具。     

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

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

地震分倒频处理技术

为刻画地震信号的局部层次特征,提高复杂地质体的地震分辨率,本文利用倒谱变换能够提高分辨率的特点,结合短时傅里叶变换,构造了倒时-频变换算法,并提出地震分倒频处理技术.文中对地震分倒频处理技术的公式推导、算法设计、分倒频处理技术提高地震纵向分辨率的机理进行了详细研究.设计出了适合三维地震资料处理的软件,并从理论上阐述了地震分倒频处理技术能够提高地震分辨率.利用倒时-频变换对多分量线性调谐信号模型做倒时-频分析,可以得到良好的时频分辨率.实际三维地震资料计算表明,该方法能够提高地震纵向分辨率,且计算效率高,适合大规模三维地震资料的处理.     

基于小波分频与径向道变换的联合压制面波方法（英文）

面波是地震记录上严重的干扰波。基于小波变换的面波衰减方法考虑了面波频率低的特点,但是当有效信号和面波干扰存在频率重叠时,其压制效果不理想。径向道变换考虑了面波与有效信号视速度差异,有利于去噪,但是在去除面波的同时也会损害有效信号的低频分量。本文根据面波和有效信号在视速度以及小波域能量的差别,将小波变换的局部分析能力与径向道变换的去噪优势相结合,提出了基于小波分频与径向道变换的联合面波压制方法。首先应用小波变换,将地震记录分解为不同频段,对出现面波的频段作径向道变换,然后再作低切滤波处理,最后利用径向道反变换后的记录与其他频段的记录进行小波重构得到去除面波的记录。两种方法的联合使用提高了小波分频去除面波的能力,同时也较好的保护了有效信号。通过模型数据实验分析和实际资料数据处理的结果表明,本文提出的方法具有较强的去噪能力和良好的保幅性能。     

S变换谱分解技术在深反射地震弱信号提取中的应用

在深反射地震资料处理中,当来自深部的有效弱信号和噪声干扰频带差异较小且难以区分时,传统滤波方法的应用会受到限制.谱分解方法是一种使用离散傅里叶变换,基于信号的频率-振幅谱等信息生成高分辨率地震图像的方法,通常用来识别介质物性横向分布特征,处理复杂介质内频谱变化和局部相位的不稳定性等问题,包括定位复杂断层和小尺度断裂等.S变换作为一种新的时频分析方法,具有自动调节分辨率的能力,近些年来被广泛应用到勘探地震、大地电磁等数据处理中,逐渐成为地球物理方法中噪声压制的有效方法之一.与常规石油反射地震资料相比,深反射主动源地震为了探测深部结构信息,常采用大药量激发方式、长排列观测系统等,导致深部有效信号基本湮灭在噪声干扰之中.针对深反射数据特点,本文结合谱分解和S变换技术,首先设计了简单的脉冲函数实验数据,证实S变换方法的有效性,同时说明谱分解方法的效果受所用时频分析方法影响较大,而其中决定分辨能力的变换窗函数的选取尤为重要.在此基础上,分别应用到深反射地震资料的单道和叠加剖面实际数据上,对比分析了传统变换谱分解和S变换谱分解的应用效果,单道资料对比结果表明:相比传统谱分解,S变换谱分解方法具有自动调节分辨率的能力,能够精确的标定深反射地震资料中弱信号不同时刻的频率分量;叠加剖面资料应用结果表明:由S变换谱分解得到的剖面结果与其他谱分解方法结果整体上具有较高的一致性,同时清晰地刻画出原叠加剖面上被噪声湮灭的低频细节特征,提高了剖面的分辨率及同相轴连续性;对比结果明显看出,Gabor变换谱分解方法得到的结果同相轴较为破碎,分析原因认为这是由Gabor变换的时频分解方法的定长窗函数所致,窗口大小不会随着信号频率的变化来调节长度,只能在处理的过程中根据一定的记录长度范围选取窗函数参数,而S变换谱分解方法在窗函数的选取时,通过时变信号的局部频率特征自动调节窗口长度,能够更好的刻画各个频段的细节特征,在深反射剖面成像应用中效果尤为明显.本文结果表明S变换谱分解技术在深地震叠加剖面上的应用有效地提高了来自深部弱反射信号的信噪比和分辨率,并刻画出了叠加剖面上所不具有的低频细节特征,在实际深反射地震资料处理中能有效保护低频弱信号获得更好的成像效果.本文为深地震反射资料中弱信号的保护处理找到一种有效的方法.     

S变换在面波去噪中的应用r

S变换是一种用于分析非平稳信号的时频变换方法, 可以很好地刻画地震信号的时频特性. 本文将S变换用于地震面波数据的噪声去除中, 首先介绍了S变换的理论基础, 然后设计了时频滤波和阈值滤波两种方法, 分别对天然地震面波数据和背景噪声数据进行去噪处理, 并与相位匹配滤波进行了比较． 结果表明, 面波数据经S变换去噪后, 群速度频散曲线的短周期部分得到改善, 能够连续追踪至6 s左右, 但长周期部分出现了缺失; S变换去噪的效果优于相位匹配滤波, 两者相结合会得到更加理想的结果．      

基于小波谱能量曲线EWT的面波压制算法

面波提取和压制是地震资料处理的关键问题之一,对地震资料处理的准确性具有重要影响.针对传统地震面波压制算法面波分离精度低、有效波损伤大的问题,提出一种基于小波谱能量曲线(Wavelet Spectrum Energy Curve,WSEC)经验小波变换(Empirical Wavelet Transform,EWT)的面波压制算法,该算法利用连续小波变换(Continuous Wavelet Transform,CWT)精细化的时频分析能力得到地震信号的WSEC,从频率和能量两个维度对信号进行自适应分解,提高面波提取精度.首先对地震记录进行CWT,根据小波谱计算各频率点的能量得到WSEC;基于WSEC极大值点的频率及ε邻域法确定频谱分割边界进行改进EWT,得到依能量和频率分解的各固有模态分量(Intrinsic Mode Function,IMF);根据各子频带的频率和能量确定面波IMF,最后对面波IMF进行带通滤波以保护其中的有效波.仿真实验表明,算法能根据地震信号的频率和能量自适应分解,实现面波的精准分离与提取;实际地震资料处理表明,算法能准确分离面波,能最大限度保护有效波,提高资料信噪比,是一种有效可靠的面波压制算法.     

基于S变换的高频全球导航卫星系统同震数据的震相识别研究

为研究高频全球导航卫星系统（GNSS）信号的P波和S波到时,本文利用2008年汶川M_S8.0地震震时部分站点记录到的1 Hz高频GNSS数据,采用广义S变换将同震信号进行二维时频域平面分解,以此对P波和S波到时进行识别。本文通过调整广义S变换中的调节因子λ_a和p得出,当λ_a＝1.05,p＝1.05时,变换所得图像中的P波和S波到时均较为明显。结果表明,S变换在高频GNSS震相识别中效果明显,可作为GNSS地震学应用中一项有效的技术手段。      

2013年海西5.0级地震波形时频特性研究

传统的傅里叶变换因缺乏信号局域性的信息,无法对地震波等非平稳信号进行全面描述,时频分析通常是分析非平稳信号较好的工具.本文利用大柴旦地震台记录到的2013年6月5日青海海西5.0级地震波形,采用赵-阿特拉斯-马克时频分布（ZAM）分别研究了该地震P波、S波段的时间-频率特征.研究结果很好地解释了地震波的非平稳性,较全面地展示了地震波的衰减特征.     

基于地震信号波形形态差异的面波噪声稀疏优化分离方法

实际地震信号通常可表示为具有波形特征差异的多种基本波形信号的线性组合,如叠前道集中的工频干扰噪声与有效波信号、面波噪声与体波信号等.选择单一数学变换方法,往往不易实现地震信号的稀疏表示.近年来发展的形态成分分析理论,通过联合多种数学变换,可实现对复杂信号的稀疏表示.本文根据单道地震记录中面波与体波信号波形结构特征的差异性,提出一种基于形态成分分析的面波噪声衰减方法.针对面波的低频、窄带以及频散特性选择一维平稳小波变换作为其稀疏表示字典,而针对体波波形的局部相关特性选择局部离散余弦变换作为其稀疏表示字典,建立基于双波形字典的形态成分分析模型,通过求解该稀疏优化问题获得最终的信噪分离结果.理论模型和实际地震资料处理证实该方法不仅能够衰减单炮地震记录中的强面波干扰噪声,同时能够更好地保护有效信号的波形特征与频谱带宽,为地震资料的后续处理和分析提供良好的数据基础.     

浅层地震反射资料的多阶振型面波反演

浅层地震反射波法和面波方法是两种相互独立发展的地震勘探方法,在各自的数据采集和处理中,对方都是作为干扰信号而存在. 本文利用浅层地震反射资料中被视为干扰的面波信号,通过成熟的多道面波勘探技术处理浅层地震反射资料,在频率-波数域中提取多阶振型面波的频散曲线,并基于该曲线反演浅地表S波速度结构. 这种方法充分开发利用了已有数据,无需单独的面波数据采集系统,同时为解释浅层地震反射资料提供了额外的信息约束. 结果表明:浅层地震反射资料中可提取出可靠的多阶振型面波频散曲线,并能给出稳定的反演结果,同时,面波反演的多解性可以通过高阶振型反演得以进一步约束;低速层的存在是观测频散曲线出现振型跳跃或呈“之”字形回折的必要条件而非充分条件.