有耗介质中脉冲响应的研究

本文详细论述了一维有耗介质中脉冲响应的正演理论及其模拟技术，深入研究了导电率、介电常数对脉冲响应的影响，并以平面脉冲电磁波在一维介质中传播的理论为基础编制了脉冲响应的仿真模拟软件，最后用实际模型验证了该模型技术的正确性和可靠性。     

粘弹VTI介质单程波正演模拟与照明分析

实际地下介质不仅表现各向异性的特征,同时还表现粘弹性的特征,对入射波和反射波的传播及能量分布造成一定的影响。为了更加准确地分析地震波在实际地下介质中的传播规律,利用改进的单程波动方程数值模拟方法,对粘弹VTI介质模型进行qP波正演模拟以及点震源和不同角度入射平面波震源的照明度分析,并对正演及照明结果进行了对比分析。通过两组模型试验,总结了地震波在粘弹VTI介质中与在弹性各向同性介质中传播的不同之处:一方面,介质存在粘弹性会对在其中传播的地震波能量造成吸收;另一方面,介质还存在各向异性,导致地震波在介质中传播速度随方向变化,对地震波的能量分布造成影响。     

用于有耗介质FDTD计算的改进的Mur吸收边界条件

从常规Mur吸收边界条件的推导出发,引入有耗介质的相应参数,使其可应用于截断边界电磁参数不连续和介质为有耗的情形。通过模拟公式计算,并与常规Mur吸收边界条件对比,说明了改进后的Mur吸收边界条件更具有优越性和普遍适用性。     

基于Kelvin黏声介质的最小二乘品质因子线性反演

品质因子Q值是描述介质黏性大小最重要的物理量,准确获取地下介质Q值分布,是黏性介质叠前属性反演、地震波正演模拟及黏性介质偏移成像的基础。为此,基于迭代反演思想,将地下介质的黏滞性引入最小二乘框架中,对品质因子Q值进行线性反演,形成了一种基于多次迭代策略的Q值反演方法;通过简单Q模型、二维盐丘Q模型试算,验证了方法的正确性及适应性。研究结果表明,利用最小二乘迭代反演方法不仅可以得到地下Q值扰动界面,而且可以得到Q绝对值分布。     

多孔介质中凝析油气体系相平衡规律的研究现状及其意义

凝析气藏相平衡过程发生于地下多孔介质中,其相平衡过程和常规PVT筒有很大的差别,因此对多孔介质中凝析油气相平衡的研究具有重要的意义。对国内外在多孔介质中,凝析油气相平衡的物理模拟技术和理论模型的研究现状进行分析认为,目前还没有能准确地测试凝析油气混合物,在真实岩心中的相态变化的物理模拟技术。在相平衡计算模型方面,考虑气固吸附、毛细凝聚、毛管力影响的相平衡模型已较成熟,凝析油与多孔介质吸附的相平衡计算模型还有待发展。目前国外尽管在多孔介质对凝析油气相平衡的影响上仍有争议,但由于该研究直接关系到凝析油气相平衡测试方法和理论、渗流理论及试井解释方法的发展,仍有必要对该课题进行深入而严密的研究。     

双相VTI介质地震波场特征及其影响因素

含流体性和各向异性是实际地下介质的两大主要表征,综合考虑这两大特征的双相各向异性介质模型更接近实际地下介质。研究这类介质模型中的地震波传播特征,有助于对实际地震波传播规律的认识,进而指导地震数据处理和解释。通过双相横向各向异性(VTI)介质波传播方程,采用伪谱法数值求解该方程进行地震波场模拟,并对波场特征及其影响因素进行详细分析。结果表明:利用伪谱法能够得到高精度的波场模拟结果,双相VTI介质中的弹性波包括快准纵波、准SV波和慢准纵波;快准纵波、准SV波的传播特征主要受固体骨架各向异性的影响,而慢准纵波的传播特征主要受弯曲度各向异性和流体黏滞系数的影响。     

双相黏弹VTI介质一阶速度—应力方程正演模拟及双程波照明研究

地层岩石具有各向异性、黏弹性及双相性等特性,因此建立更加精确的地下介质模型,研究地震波传播规律,对认识复杂介质中地震波的传播特征和实际资料的解释有着重要的实际意义.为此,基于黏弹性广义标准线性体(GSLS)模型,首先推导了双相黏弹VTI介质的一阶速度—应力方程并进行正演模拟,与单相弹性各向同性介质相比,双相黏弹VTI介...     

双相PTL介质地震波场模拟

从基于Biot理论建立的双相各向异性弹性波方程出发,推导出二维双相PTL介质弹性波方程。在此基础上,高阶实现了高阶有限差分数值模拟,并对双相PTL介质模型和实际测井数据模型进行模拟试算。基于模拟结果,分析了各类波的传播规律,讨论了弹性参数和耗散系数对地震波传播的影响。结果表明：在双相PTL介质中存在快纵波、慢纵波和SV波,并具有各向异性,两种纵波有明显的区别,快纵波的速度远大于慢纵波的速度;慢纵波具有很强的衰减性,能否观察到慢纵波是与含流体地层介质的耗散性质有关。     

利用多波地震资料初步探讨地下介质各向异性

在油气勘探和开发中,对地下介质各向异性的研究是个前沿课题.利用横波地震技术开发裂缝性油气田已成为当今的热门课题.文中利用我国海上首次采集的多波地震资料,对莺歌海盆地3个构造的地下岩石介质的各向异性进行了初步探讨,得出了很有意义的结论.     

利用多波地震资料初步探讨地下介质各向异性

在油气勘探和开发中，对地下介质各向异性的研究是个前沿课题，利用横波地震技术开发裂缝性油气田已成为当今的热门课题。中利用我国海上首次采集的多波地震资料，对莺歌海盆地3个构造的地下岩石介质的各向异性进行了初步探讨，得出了很有意义的结论。     

声波完全匹配层吸收边界条件的改进算法

传统位移形式的声波完全匹配层(PML)吸收边界条件以三分裂位移参量来构建,需要求解时间3阶导数,计算量大,计算时间长.为此,提出了一种改进的声波分裂PML吸收边界条件.阐述了传统和改进的PML吸收边界条件的构建原理,通过数值模拟讨论了改进和传统的PML吸收边界以及低阶Higdon吸收边界对边界反射的吸收效果.结果表明,低阶Higdon吸收边界条件吸收效果较差,在波场快照和模拟记录中存在较强的边界反射;改进的PML吸收边界条件与传统的PML吸收边界条件效果相当,均能有效地吸收衰减任意角度的边界反射波.改进算法用褶积近似运算代替时间3阶导数的求解,因此简化了计算过程,减少了计算量,是一种高效稳健的算法.     

交错网格有限差分正演模拟的联合吸收边界

三维声波方程交错网格有限差分正演模拟中的边界问题一直是热点问题。完全匹配层吸收边界（PML）具有较强且稳定的吸收效果,但必须具有一定的边界厚度才能吸收干净,这就增大了三维正演模拟的模型空间,即增加了运算量; Higdon边界能消除任意角度入射波的边界反射,也具有较强稳定性,但该高阶吸收边界离散化后过于复杂,而低阶时吸收效果不如PML边界。因此,基于对PML吸收层中的平面波传播规律的研究,重新推导PML最外层的Higdon吸收边界条件,得到含PML吸收系数的新的Higdon吸收边界条件。联合吸收边界不仅可使用较小厚度（相对于单纯PML边界）的PML层对分量进行衰减,而且在PML边界外层,能应用新推导的Higdon吸收边界条件对反射波进行匹配吸收。在相同吸收效果下,联合吸收边界大幅度降低了PML厚度,减小了运算量,得到精确的模拟结果。     

VTI介质准P波频率空间域组合边界条件研究

讨论了VTI介质中准P波频率空间域的边界条件。首先给出了VTI介质准P波波动方程，阐述了特征分析方法的基本原理，讨论了边界上的反射系数与入射角和度量准纵、横波各向异性强度因子的关系；然后利用特征分析法结合Kelvin-Christoffel方程，构造了准P波波动方程在不同边界和角点的频率域吸收边界条件，利用最佳匹配层法构造了衰减边界条件；最后利用数值模拟对构造的边界条件进行了验证。为了得到好的吸收效果，将吸收边界条件和衰减边界条件有机地结合起来，即先用最佳匹配层法衰减传播到边界的入射波能量，然后再用吸收边界条件吸收边界反射，最终使边界反射降低到可以忽略不计。数值模拟的共炮点记录说明了组合边界条件的良好效果。     

多轴卷积完全匹配层吸收边界条件

在通常情况下,完全匹配层（PML）吸收边界条件较其他吸收边界条件具有更优越的吸收性能,已广泛用于地震波正演模拟。但传统的PML吸收边界条件也存在一定缺陷,如不分裂卷积完全匹配层（C-PML）吸收边界在某些介质中不稳定,分裂的多轴完全匹配层（M-PML）吸收边界的吸收效果不好。为此,本文提出了多轴卷积完全匹配层（MC-PML）吸收边界条件,既可提高C-PML吸收边界的稳定性,又可改善M-PML吸收边界的吸收效率,且采用不分裂算法实现,可有效降低计算成本。正演模拟结果表明,MC-PML吸收边界条件具有更高的计算稳定性和较好的吸收效果。     

三维各向同性介质弹性波方程交错网格高阶有限差分法模拟

给出了三维各向同性介质中一阶应力一速度弹性波方程交错网格任意偶数阶精度有限差分格式，推导出了三维各向同性弹性介质完全匹配层吸收边界条件公式和相应的交错网格高阶有限差分格式。对三维French模型进行了弹性波模拟，结果表明，该方法模拟精度高，边界吸收效果好。在三维French模型的XOZ平面和yoz平面的弹性波场快照中可以见到断面反射波、侧面反射波和散射波等波场特征。多次侧面反射波和多次散射波说明不同三维构造所形成的波场是相互影响的。     

Lax-Wendroff集中质量有限元法地震波场模拟

为了高精度、高效地震波模拟的需要,本文在空间离散上运用集中质量三角网格有限元法求解弹性波方程,在时间离散上运用Lax-Wendroff方法获得时间四阶精度,提出了模拟弹性波传播的Lax-Wendroff集中质量有限元法(LWFEM).为了防止人工截断边界而引起的虚假反射,构造了二阶位移形式的PML吸收边界条件.在周期性网格中,构造LWFEM频散分析的一般特征值问题,得到了LWFEM的稳定性条件.在数值实验中通过与中心差分有限元法(CDFEM)、Runge-Kutta有限元法(RKFEM)以及Newmark谱元法(NSEM)等常见方法的对比,证实了LWFEM弹性波模拟的高效、高精度性及对复杂模型的适应性.     

基于完全匹配层构建新方法的2.5维声波近似方程数值模拟

传统的吸收边界条件只对有限入射角和有限频率范围内的波具有很好的吸收效果,而在有效吸收范围外的边界反射将对数值模拟结果造成污染。为此,在2.5维声波数值模拟中引入完全匹配层吸收边界条件,并提出了一种完全匹配层构建新方法。该方法将完全匹配层构建在研究区域的右边界和下边界的外部,同时改变相应的衰减因子,使之关于完全匹配层的中心呈对称分布。在采用有限差分法计算波场时,将完全匹配层外边界的网格点与相应的研究区域边界的网格点连接起来,使波在传入完全匹配层之后又传入到研究区域中。由于完全匹配层良好的吸收效果,波在完全匹配层中被完全吸收,基本不会再传入到研究区域。在构建相同数量的完全匹配层的情况下,完全匹配层构建新方法使波在传播过程中比以往的构建方法多通过一倍的完全匹配层介质,从而提高了对边界反射的吸收效果。通过均匀速度模型和Marmousi速度模型证实,完全匹配层构建新方法对边界反射的吸收效果较好。     

基于双相介质BISQ模型的地震波正演模拟

Biot流动机制和喷射流动机制是含流体孔隙介质中的两种重要的力学机制,同时包含这两种机制的BISQ模型能够更好地反映孔隙介质中弹性波的传播规律。本文从BISQ模型的基本方程出发,推导了完全匹配层(PML)吸收边界的交错网格高阶有限差分方法的表达式,对双相各向同性介质中地震波场进行了数值模拟;然后通过调整孔隙度、渗透率等储层参数,对比分析了储层参数对双相介质波场特征的影响。数值模拟结果表明,基于BISQ模型的正演模拟正确反映了双相各向同性介质的波场特征;储层参数取值变化对波场特征的影响也十分明显。     

完全匹配层人工边界条件中的衰减因子分析

为了得到地下地质体的地震波场响应特征,需要进行三维地震波场数值模拟,但其计算量很大。利用有限差分方法进行数值模拟时,减少计算量的途径有两种:一是差分格式设计;另一种是边界条件的选择。从完全匹配层(PML)吸收边界条件的原理出发,总结了目前PML人工边界条件中的常用衰减函数,说明了衰减函数设计的重要性。数值模拟结果表明,余弦型吸收衰减因子的波场模拟效果要优于其它吸收衰减因子。