简介：在地震频段内,中观尺度局域流是引起孔隙岩层速度频散与衰减的主要原因,研究中观尺度孔隙中局域流的过渡频率及其对应的衰减（逆品质因子）峰值频率之间的关系,有助于深入分析孔隙介质内部构造的细节。基于Biot孔隙弹性理论得出时间域内局域流流体通量,引入傅立叶变换得到频域上的流体通量。为避免不同单元体之间局域流的相互影响,建立双相孔隙介质单重单元体模型,并且在流体通量基础上选择含水层的过渡频率来等价局域流的过渡频率。此外,针对不同孔隙介质参数,分析了局域流整体过渡频率及相应的衰减峰值频率的变化特征。为了了解边界条件对计算结果造成的影响,建立了多重单元体模型,引出局部局域流过渡频率及相应的衰减峰值频率。数值模拟结果表明,在低饱和度状态下,两种频率的变化趋势相同,但变化幅度存在较大差异;在高饱和度状态下,两种频率的变化趋势与变化幅度都很接近。

简介：全波形反演是一种高精度的地震成像方法，可以对地下介质物性参数模型进行准确的重构。然而在实际应用中，尤其是在三维复杂介质反演中，计算成本太大是该方法的一个重要缺陷。将混叠震源技术引入到频率域全波形反演中可以大幅度地降低计算成本，提高反演效率。但是使用震源编码技术也带来了两个问题：一方面，参与编码的各个震源之间会产生“串扰噪声”，导致反演结果中出现假象；另一方面，基于震源编码的频率域全波形反演方法周围噪声较为敏感，使该方法对含噪数据反演质量较差。本文引入一种频率组编码方法来压制“串扰噪声”，并基于震源编码技术提出一种频率域自适应全波形反演方法，通过一个与频率相关的自适应选择机制，将常规频率域全波形反演方法和基于震源编码的全波形反演方法联合起来，在保证反演质量的同时也最大程度地提高了反演效率。

简介：频率域全波形反演虽然克服了时间方向上的局部极小值问题,但是地下介质的复杂性使其在空间域仍然存在局部极小值缺陷。在优化梯度法基础上,本文采用预条件双共轭梯度稳定算法和多重网格方法计算反演中的波场传播和目标函数的梯度,在保证计算速度的同时,减小计算机内存的消耗。频率域波形反演和多重网格的多尺度性质有�

简介：频率域全波形反演充分利用全波场的振幅、相位以及频率信息，采用较少的频率便能反演得到精度很高的速度模型。本文以有限单元法为基础，对起伏地形条件下二维声波频率域全波形反演进行了研究。在正演算法中，针对截断边界问题，并考虑多频率联合反演中计算区域采用同一套剖分网格的需求，提出了一种适用于起伏地形的衰减边界条件算法。该算法的核心思想是在控制方程波数项中引入衰减因子，通过一定方式调节衰减因子使得声波在衰减层中充分衰减，达到压制截断边界影响的目的。根据指数衰减规律，文中推导出了一种新的衰减因子计算公式，并给出了不同频率条件下衰减层厚度计算公式；在反演算法中，采用共轭梯度法求解高斯牛顿反演迭代方程组，避免直接求解雅克比矩阵和Hessian矩阵带来的巨额计算量，并采用相同的反演模型，对比分析了不同初始模型和频率组合对全波形反演结果的影响。起伏地形模型数值模拟和全波形反演数值试验表明，本文提出的指数衰减边界条件算法和基于该算法的全波形反演算法具有很好的应用效果。

简介：基于曾新吾和MacBeth提出的时域内双源累积旋转方法,建立了频率相关介质中分析多分量VSP数据的横波分裂参数提取算法（DCTF）。该算法可以在频域中针对单个频率提取横波各向异性参数（快横波的极化方向以及快、慢横波间的时间差）,从而避免了目前常用方法中使用窄带通滤波可能带来的误差。通过对地震合成记录的数值分析,确定了该算法的可行性和正确性,并与目前常用方法的应用结果进行了比较。结果显示,频率相关横波分裂参数可以利用DCTF从地震四分量数据中直接提取。在地震频率范围内,含较大尺度裂缝时各向异性参数将表现出频率相关性,这意味着在地震频率范围会出现频散。随着频率的增加,各向异性有降低的趋势。

简介：频率空间域地震波数值模拟具有独特的优势:可以同时模拟多源的波传播、每个频率之间独立并行地计算、计算频带选择灵活、不存在累计误差、容易模拟粘弹性介质中地震波传播。但是该方法的最大瓶颈是对于计算机内存的巨大需求。我们使用压缩存储系数矩阵的方法,极大地减少了计算机内存的需求量。同时为了减少短差分算子的数值频散,引用了频率空间域25点弹性波波动方程的差分格式,并使用了最小二乘意义下求出的优化差分系数。为了克服边界反射,采用了最佳匹配层吸收边界条件。数值模拟试验证明:用压缩存储系数矩阵及优化差分系数的频率空间域25点差分格式进行弹性波正演模拟,可以减少数值频散,提高计算精度。使用较大的网格间距,降低计算机内存需求,并保持较高的计算效率。该正演方法为后续弹性波偏移和弹性参数反演提供较好的基础。

简介：为了探明介质的非完全弹性对瑞利型槽波传播的影响以及利用瑞利型槽波品质因子QR反演煤层碳化程度、裂隙、煤层厚度等煤层属性,在本文中,我们将复速度引入瑞利型槽波频散方程,计算了煤层瑞利型槽波的品质因子QR曲线;采用控制变量法,分析了瑞利型槽波品质因子QR随煤层厚度、煤层纵横波品质因子以及围岩纵横波品质因子的变化规律。研究表明:瑞利型槽波品质因子QR曲线与群速度曲线变化趋势一致;埃里相（Airy-phase）频率附近瑞利型槽波的品质因子QR最小,且埃里相频率随煤层厚度增加而减小;瑞利型槽波品质因子QR随煤层横波品质因子QS2的增大而增大。利用瑞利型槽波品质因子QR曲线可以实现对瑞利型槽波的吸收衰减进行补偿、对煤层厚度进行预测以及岩性成像。更多还原

简介：本文以中观孔隙结构的White模型为基础，构建了部分饱和孔隙介质模型，利用Biot方程的建立思路和Johnson推导的体变模量，推导了部分饱和孔隙介质中的纵波方程，并以平面波为例，求取了方程的衰减系数，分析了地震勘探频带范围内地震波的衰减特性。结果表明：在部分饱和孔隙介质中，地震波在低频段也会发生明显的衰减和频散现象，频率越大，衰减越大；且第二纵波的衰减比第一纵波更为明显；这一结论弥补了Biot理论在描述地震勘探频带范围内波的衰减现象的不足。文中还研究了孔隙度、饱和度和模型内径尺寸对第纵波衰减特性的影响机理，主要表现在在地震勘探频带范围内，波的衰减随孔隙度的增大而增大，随含油气饱和度的增大而减小，当孔隙内径尺寸小于二分之一外径尺寸时，波的衰减随内径尺寸的增大而增大，当内径尺寸大于二分之一外径尺寸时，波的衰减随内径尺寸增大而减小。

简介：由于地震数据中包含的噪声在不同频率或者频带数据中的分布强度存在差异,使得全频带数据上进行的噪声衰减处理改变了地震反射波信号的动力学特征,干扰后期的地震资料解释、储层预测、油气检测等问题,提出边界和振幅特性保持自适应噪声衰减方法。首先应用小波包变换对全频带地震数据进行多频段划分,然后对分频段数据进行非线性各向异性倾角导向边界保持自适应滤波处理。在该方法中,由结构张量计算的扩散张量实现自适应地确定平滑滤波方向,加入的不连续结构置信度量和不连续性算子自适应地控制不连续结构特征的保持程度,引入的去相关滤波迭代停止准则自适应地确定滤波迭代次数。这些参数的引入具有减少处理人员的干预和人为的主观性,且执行简单的特点。对合成地震记录和实际地震记录处理结果表明,提议的方法能够自适应地衰减地震数据中噪声,同时既能保持地震反射波中有效的不连续性信息,也能有效地保持有效信号的频率分布规律。能够为后期的地震资料解释和分析提供高品质的基础数据。

简介：由于地震数据中包含的噪声在不同频率或者频带数据中的分布强度存在差异,使得全频带数据上进行的噪声衰减处理改变了地震反射波信号的动力学特征,干扰后期的地震资料解释、储层预测、油气检测等问题,提出边界和振幅特性保持自适应噪声衰减方法。首先应用小波包变换对全频带地震数据进行多频段划分,然后对分频段数据进行非线性各向异性倾角导向边界保持自适应滤波处理。在该方法中,由结构张量计算的扩散张量实现自适应地确定平滑滤波方向,加入的不连续结构置信度量和不连续性算子自适应地控制不连续结构特征的保持程度,引入的去相关滤波迭代停止准则自适应地确定滤波迭代次数。这些参数的引入具有减少处理人员的干预和人为的主观性,且执行简单的特点。对合成地震记录和实际地震记录处理结果表明,提议的方法能够自适应地衰减地震数据中噪声,同时既能保持地震反射波中有效的不连续性信息,也能有效地保持有效信号的频率分布规律。能够为后期的地震资料解释和分析提供高品质的基础数据。

简介：尾波是持续时间比主要波型长得多的多重散射波，它包含了直达波之外的部分有用信号。在油气田的开采过程中，由储层物性参数的微小变化而导致的储层速度的微小变化对初至的影响很小，无法通过初至变化直接观测，但是因为尾波为多重散射波，故储层速度的微小变化会在尾波的传播过程中被放大，因此利用尾波可以观测到这种变化。本文通过实验与数值模拟研究了尾波干涉方法在储层微小时移差异监测方面的有效性。合成地震数据是基于部分MarmousiII模型采用有限差分波动方程正演方法计算获得，为了模拟时移地震中的储层波速变化，在模型中选取波速发生微小变化的目标区，计算速度扰动前后的合成地震数据，进而观测尾波变化。研究表明利用尾波干涉可检测储层物性的微小变化，为时移地震信号监测提供新的方法和手段，从而提高时移监测的准确度，也为开发生产工作提供可靠的指导。

简介：瑞利波勘探主要是建立在弹性介质分层半空间模型基础上的。当实际地层中包含孔隙介质层时，需要将孔隙介质简化为弹性介质来进行分析和处理。这种简化处理究竟对瑞利波勘探造成怎样的影响是本文所研究的问题。本文基于弹性介质与孔隙介质共同构成的分层半空间模型，首先推导了分层半空间中两种介质处于不同相对位置时的瑞利波频散方程，解决了不同阶数矩阵之间的变量传递问题然后，针对传递矩阵法在求解频散函数时可能出现的溢出问题，给出了一种可以有效提高计算范围的解决方案；同时，提出了一套新的数值算法用于复频散方程的快速求解。数值计算结果表明：当孔隙介质位于半空间表面时对低频瑞利波频散特性的影响最为显著，而在其它情况下对瑞利波频散特性的影响相对较小。

简介：这篇文章重点研究改进的Gabor小波（improvedGaborwavelet，IGW）变）并讨论了它在地震信号处理和解释中的应用ThispaperintroducesanimprovedGaborwaveletanditscompletetransform，andmainlyanalysestheirpropertiesanddiscussesapplicationsofthesepropertiesinseismicsignalprocessandinterpretation。改进的Gabor小波变换具有以下特性：1）IGWT把时域信号映射到时间一频率域，而传统Gabor小波变换把时间信号映射到时间一尺度域；2）IGWT可用于信号分频，通过固定变换的主频参数dominantfrequency，并变换能提取相应的子带信号，且其主频部分的信息与原信号相应频率部分的信息一致，通过调节变换的分辨率因子，变换能有效控制子带信号的带宽；3）用IGWT和IGWIT构建的滤波器有良好的时一频局部性，在指定时一频范围内能实现针对性滤波。文章用仿真实验和实际用例验证IGWT的这些特性，并在提高地震信号分辨率、地震信号分频和识别小断层等地震信号处理和解释等方面的应用中取得良好效果。

简介：阵列声波信号是典型的非线性、非平稳信号，Hilbert～Huang变换（HHT）是处理非平稳信号的一种比较新的时频分析方法。通过对信号进行经验模态分解（EMD）和对瞬时频率的求解，可以获得声波信号的时一频谱。其关键技术就是进行经验模态分解，任何非平稳的信号都可以分解为有限数目并且具有一定物理意义的固有模态函数。EMD方法可以理解为以声波信号极值特征尺度为度量的时频滤波过程。滤波器充分保留了声波信号本身的非线性和非平稳特征，在声波信号的滤波和去噪中具有很大的优势。文中介绍了HHT时频滤波的实现过程，并列举了一些声波测井波列实例，说明了该方法的有效性。

简介：阵列声波信号是典型的非线性、非平稳信号,Hilbert-Huang变换(HHT)是处理非平稳信号的一种比较新的时频分析方法。通过对信号进行经验模态分解(EMD)和对瞬时频率的求解,可以获得声波信号的时-频谱。其关键技术就是进行经验模态分解,任何非平稳的信号都可以分解为有限数目并且具有一定物理意义的固有模态函数。EMD方法可以理解为以声波信号极值特征尺度为度量的时频滤波过程。滤波器充分保留了声波信号本身的非线性和非平稳特征,在声波信号的滤波和去噪中具有很大的优势。文中介绍了HHT时频滤波的实现过程,并列举了一些声波测井波列实例,说明了该方法的有效性。