简介:频率空间域地震波数值模拟具有独特的优势:可以同时模拟多源的波传播、每个频率之间独立并行地计算、计算频带选择灵活、不存在累计误差、容易模拟粘弹性介质中地震波传播。但是该方法的最大瓶颈是对于计算机内存的巨大需求。我们使用压缩存储系数矩阵的方法,极大地减少了计算机内存的需求量。同时为了减少短差分算子的数值频散,引用了频率空间域25点弹性波波动方程的差分格式,并使用了最小二乘意义下求出的优化差分系数。为了克服边界反射,采用了最佳匹配层吸收边界条件。数值模拟试验证明:用压缩存储系数矩阵及优化差分系数的频率空间域25点差分格式进行弹性波正演模拟,可以减少数值频散,提高计算精度。使用较大的网格间距,降低计算机内存需求,并保持较高的计算效率。该正演方法为后续弹性波偏移和弹性参数反演提供较好的基础。
简介:全波形反演是一种高精度的地震成像方法,可以对地下介质物性参数模型进行准确的重构。然而在实际应用中,尤其是在三维复杂介质反演中,计算成本太大是该方法的一个重要缺陷。将混叠震源技术引入到频率域全波形反演中可以大幅度地降低计算成本,提高反演效率。但是使用震源编码技术也带来了两个问题:一方面,参与编码的各个震源之间会产生“串扰噪声”,导致反演结果中出现假象;另一方面,基于震源编码的频率域全波形反演方法周围噪声较为敏感,使该方法对含噪数据反演质量较差。本文引入一种频率组编码方法来压制“串扰噪声”,并基于震源编码技术提出一种频率域自适应全波形反演方法,通过一个与频率相关的自适应选择机制,将常规频率域全波形反演方法和基于震源编码的全波形反演方法联合起来,在保证反演质量的同时也最大程度地提高了反演效率。
简介:频率域全波形反演充分利用全波场的振幅、相位以及频率信息,采用较少的频率便能反演得到精度很高的速度模型。本文以有限单元法为基础,对起伏地形条件下二维声波频率域全波形反演进行了研究。在正演算法中,针对截断边界问题,并考虑多频率联合反演中计算区域采用同一套剖分网格的需求,提出了一种适用于起伏地形的衰减边界条件算法。该算法的核心思想是在控制方程波数项中引入衰减因子,通过一定方式调节衰减因子使得声波在衰减层中充分衰减,达到压制截断边界影响的目的。根据指数衰减规律,文中推导出了一种新的衰减因子计算公式,并给出了不同频率条件下衰减层厚度计算公式;在反演算法中,采用共轭梯度法求解高斯牛顿反演迭代方程组,避免直接求解雅克比矩阵和Hessian矩阵带来的巨额计算量,并采用相同的反演模型,对比分析了不同初始模型和频率组合对全波形反演结果的影响。起伏地形模型数值模拟和全波形反演数值试验表明,本文提出的指数衰减边界条件算法和基于该算法的全波形反演算法具有很好的应用效果。
简介:在72届SEG会上,BP公司LindaHodgson等人提出了一种新的噪声压制方法——频率切片滤波(FSF)。在复杂的X-Y频率域中,FSF应用二维平滑滤波器直接去噪。它可以任意选择特定的频率处理范围进行有针对性的滤波,而数据的其余部分不受影响;任何适合于空间滤波的噪声都能够去除,特别是低频噪声的消除和剩余多次波能量的衰减。具体实施过程分为4个步骤:①应用一维快速傅立叶变换到目的层时窗;②检查X-Y频率域数据体,确定频率范围,有针对性地进行噪声衰减;③在每一个频率切片上进行平滑处理,数据的其余部分不改变;④进行傅立叶逆变换,得到滤波结果。以下选取北海2个油田的实例展示其良好的滤波效果。
简介:在地震频段内,中观尺度局域流是引起孔隙岩层速度频散与衰减的主要原因,研究中观尺度孔隙中局域流的过渡频率及其对应的衰减(逆品质因子)峰值频率之间的关系,有助于深入分析孔隙介质内部构造的细节。基于Biot孔隙弹性理论得出时间域内局域流流体通量,引入傅立叶变换得到频域上的流体通量。为避免不同单元体之间局域流的相互影响,建立双相孔隙介质单重单元体模型,并且在流体通量基础上选择含水层的过渡频率来等价局域流的过渡频率。此外,针对不同孔隙介质参数,分析了局域流整体过渡频率及相应的衰减峰值频率的变化特征。为了了解边界条件对计算结果造成的影响,建立了多重单元体模型,引出局部局域流过渡频率及相应的衰减峰值频率。数值模拟结果表明,在低饱和度状态下,两种频率的变化趋势相同,但变化幅度存在较大差异;在高饱和度状态下,两种频率的变化趋势与变化幅度都很接近。
简介:为研究环境温度对珠江黄埔大桥频率监测的影响,首先要对大桥模态频率进行连续识别。珠江黄埔大桥上架设的监测系统为强震动台阵,相较于其他健康监测系统测点较少,因此,应基于强震动台阵系统的特点,选取合适的方法对大桥频率进行识别。本文通过对比分析平均正则化功率谱法(ANPSD)、频域分解法(FDD)和协方差驱动的随机子空间法(Cov—SSI)的识别结果,择优应用于珠江黄埔大桥的频率自动识别中。采用珠江黄埔大桥强震动台阵记录的2013年4月至11月加速度响应数据进行频率识别,识别结果可用于观测和研究大桥频率在环境影响下的波动情况。
简介:分界面的频率属性分析是一项新的地震界面属性拾取和解释方法。与时间域反射波成像、地震振幅研究相比,地震信号频率域的研究较弱。本文系统分析了脉冲信号的各类频率意义,指出瞬时频率属性不能提供在到达时间位置上的准确的简谐成分频率信息;引用傅立叶积分的概念解释了时频分析方法原理,介绍了(俄罗斯)穆申教授的时频分析算法;给出了时频分析剖面、分频剖面、层位切片解释模型以及时频分析数据体的时间剖面显示方式(多频合成时间剖面、多频信息RGB三原色合成时间剖面);发展了三维地震数据体的层位地震数据时频分析技术。文后对三维数据体频率属性结果的应用给出了一些提示。