复杂地表的单程波动方程地震叠前正演

基于数学检波器和等时叠加原理,实现了复杂地表的单程波动方程地震叠前正演模拟。该方法采用虚拟的数学检波器接收地下反射地震信号,可灵活地将接收点布置在地表的任何地方,从而满足地表起伏的要求。此外,根据等时叠加原理,该方法采用单程波动方程进行波场延拓和成像,计算简单、快速。通过复杂正断层的数值模拟,得到了高信噪比的共炮地震记录;采用适用于起伏地形的深度偏移方法对该共炮记录进行了叠前深度偏移,实现了地震波的偏移归位。从而证明了本文提出的适用于起伏地表的单程波动方程地震叠前正演方法的正确性和准确性。     

三维波动方程正演及模型应用研究

为了真实准确地反映三维地质体的波场特征,在频率-波数域将二维波场延拓算子推广到三维空间,采用三维波动方程延拓方法实现了三维地质模型的快速叠后正演.该方法可以采用相位移加插值方法处理一定的横向变速情况,可以更加灵活方便地模拟地下复杂的三维地质体.首先进行了三维French模型数值模拟,得到了和实际物理模型实验结果相一致的正演记录,并对比分析了三维偏移剖面和二维偏移剖面的偏移效果;然后进行了三维缝洞地质模型的正演计算,得到了高信噪比的正演记录.模拟结果验证了三维正演和偏移方法的正确性,以及利用该方法进行缝洞地质体识别和研究的可行性.     

相位移加有限差分法波动方程正演模拟

本文提出用相位移加有限差分法来实现速度纵横向变化，复杂陡倾地质模型的正演模拟，为复杂波场的地震，地质解释提供了精确、实和的方法。     

火山岩波动方程正演模拟研究

松辽盆地火山岩埋藏较深,下伏地层成像精度极低,对构造解释造成了一定的影响.为了了解火山岩体对下伏地层的影响程度,采用波动方程正演模拟技术,研究了9种波阻抗组合和两类岩性组合情况下火山岩理论模型的地震响应特征,其波场快照能直观地展现地震波在火山岩体中的传播过程.对正演模拟结果进行了分析,认为火山岩体的底界面反射不清晰是火山岩体对下伏地层屏蔽作用所致;火山岩体厚度增大时,下伏连续地层的反射中断现象更明显,以致出现假断层和假构造现象.     

反偏移法三维波动方程正演

三维地震资料的处理与解释往往都需要三维正演模型给予验证，以得到准确的处理、解释结果。本文提出的三维正演方法是一种偏移算法的逆运算方法－有限差分反偏移法。这一算法具有算法简单、运算量小及深度方向分辨率高的特点，是一种快速、高精度的正演算法。反偏移法虽是偏移算法的逆运算法，却出现常规偏移算法中所没有的很多问题。为此，本文进行了三种不同反偏移方法试验，并采用了一种新的浮动坐标变换，使正演结果精度得到很大     

薄互层油藏模型井间地震弹性波方程正演模拟研究

 本文采用非均匀介质中弹性波方程交错网格高阶有限差分法进行数值模拟,并通过对正演过程中每一时刻的弹性波场求散度和旋度实现了弹性波场分离。其散度场以纵波场为主,旋度场以转换波和横波为主;采用模型多尺度网格化、参数化技术,建立了薄互层油藏2D井间地震弹性介质模型。数值模拟试验表明,弹性波方程数值模拟能够更加真实地模拟井间地震的波场。正演模型与偏移成像结果剖面吻合很好。     

基于黑油模型的时移地震正演模拟

时移地震的正演模拟是时移地震油藏监测可行性研究的主要手段，它将油藏数值模拟通过岩石物理参数属性转换和地震正演模拟联系起来，具有重要的理论意义和实用价值。时移地震正演模拟的整个过程为建立油藏初始模型、油藏数值模拟（建立油藏时间模型）、岩石物理计算（流体置换方程）、地震正演模拟，最后通过理论模型试算验证了方法的正确性和有效性。     

多次反射可控制的波动方程共炮记录正演模拟

本文首先把声波方程分解为两个分别适用于做正演和偏移的单向旁轴方程，然后用适合于做正演的单向方程业合成共炮记录。通过公共推导，提出合成共炮记录的时间一致性原理和两个过程－正过程和逆过程，即对于简单的水平各向同性层状介质，仅用正过程就能精确地合成共炮记录；而对于复杂介质，则需借鉴射线法两点问题和波动方程共炮记录叠前深度偏移原理，给定炮点和检波点，在时间一致性原理的基础上，应用逆过程方法，以单向声波方程     

波动方程法地震波正演数值模拟研究综述

波动方程法地震波正演数值模拟主要包括有限差分法、连续有限元法,谱方法以及新兴的间断有限元法。本文总结了这些方法近年来的新进展,探究了各方法在应对起伏地表等复杂构造和提高模拟效率等方面出现的新技术,并对该领域的发展方向和趋势进行了分析和展望。     

基于Qt的地震交互建模

地震交互建模技术在地震勘探中具有十分广泛的应用。本文以Qt为程序开发工具,详细论述了研发人工交互方式建立速度模型的过程,并给出了基于所建模型的有限差分地震正演模拟记录。     

基于波动方程数值模拟的地震构造精细解释技术

为了准确识别生物礁等复杂地质体的分布范围、内部储集层结构等细节信息,提出了一种新的基于波动方程数值模拟的地震构造精细解释技术。该技术采用频率-波数域单程波动方程进行地震数值模拟,通过对比模拟记录与偏移前、后的实际地震剖面的主要反射特征,逐次修改地震构造解释模型,从而得到最终的构造精细解释结果。该技术采用频率-波数域单程波场延拓算子进行计算,不仅计算效率高,而且波场模拟记录信噪比高,便于目标地质体有效反射波场特征的对比。该技术充分利用了偏移前、后地震剖面的反射特征,有助于降低地震构造解释的多解性,获得精细的地震构造解释模型。     

三维快速高精度地震波正演数值模拟方法及其应用

如何有效提高三维地震波正演数值模拟精度和计算效率一直是勘探地球物理学研究的重要问题。为了克服常规中心有限差分法较难快速提高差分精度的缺陷和一阶双曲型波动方程内存占用多、计算量大、引入变量较多的困难,采用高阶交错网格有限差分法直接求解三维地震波动方程,推导的高阶差分格式计算形式简单,可以推广于求解任意偶数阶时空导数,同时给出其稳定性条件。在人工边界处,对比了镶边法和常规旁轴近似法两种吸收边界条件。从三维似French模型的正演结果看出,采用的高阶交错网格差分算法在快速有效地提高数值模拟精度的同时,大大提高了计算效率,同时结合镶边法吸收边界条件还可有效压制边界反射,提高整个计算域内波场的信噪比。     

正演子波响应特征及逆时成像分析

采用一阶形式和二阶形式的地震波动方程,开展了二维地震子波正演响应理论研究,并以二维均匀介质速度模型为例,将不同相位的 Ricker子波作为输入子波,进行了数值模拟。结果表明,一阶形式的子波响应与输入子波波形差异较大,且存在较大的相位移,同时,其比值 谱与频率呈线性正比关系;而二阶形式的子波响应与输入子波波形差异较小,相位移较小,其比值谱几乎不随频率变化。因此二阶形式的波动 方程具有更好的保形特征,可以实现任意期望子波波形的正演数值模拟。在上述研究的基础上,讨论了逆时子波波场响应特征,并以一维层状 介质反射系数模型为例,进行了一阶形式和二阶形式波动方程的相关型和反射率型逆时偏移成像试验,结果表明,当震源波场与接收波场的子 波波形相一致时,低频逆时噪声能量最弱,成像效果最佳;反之,成像剖面中的反射界面与真实的地层界面位置会存在时差,同时,低频逆时 噪声能量较强。因此,保持子波波形的一致性是实现高精度相关型、反射率型逆时偏移成像的前提。     

基于块状有理Krylov方法的大地电磁三维模型降阶正演

三维大地电磁正演需要求解两个极化源在不同频率下的电磁场分布,计算量巨大。基于块状有理Krylov方法,文中实现了大地电磁三维模型降阶快速正演计算。所用算法的创新点包括：①将大地电磁的源项显式表示为平面电流源,将随频率变化的电场响应表示为一个传递函数与电流源常矢量的乘积,从而可通过构建有理Krylov子空间实现所有频率电场响应的快速求解,避免多次求解不同频率的大型稀疏线性方程组;②采用块状Krylov技术,将TE和TM极化源表示为块状源矢量,将求解两个极化源的正演响应简化为构建一个块状有理Krylov子空间。引入渐近收敛公式得到了块状有理Krylov方法的最优化单个重复极点,结合直接求解器,将大地电磁三维正演的计算量降为一次系数矩阵分解和几十次矩阵回代。该算法在保证正演精度的同时,极大地提高了正演速度。半空间模型和三维DTM1模型的正演数值结果表明,相比常规的逐个频率的正演求解方法,块状有理Krylov方法可显著提高正演速度。