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

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

一种新的弹性波数值模拟吸收边界条件

 在弹性波方程有限差分波场数值模拟中,由于计算模型的限制,不可避免地导致很强的人为边界反射的出现。为了使计算得到的边界波场值更接近真实的边界波场值,在波场数值模拟中人们利用边界条件来减少来自计算区域边界的人为反射能量。与声波吸收边界条件相比,弹性波边界条件要复杂得多,到目前为止弹性波边界条件算法研究仍然是该领域的一个重要课题。本文借鉴Clayton等处理吸收边界反射的思想,从弹性波波动方程出发,提出了一种新的在各向同性介质中进行有限差分弹性波场数值模拟的吸收边界条件。波场计算表明：该方法计算量小、应用简单、能有效地减少来自模型边界的反射波场,其计算效率与Clayton等处理吸收边界条件所花费的时间相当。     

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

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

各向同性介质弹性波交错网格有限差分正演模拟

在弹性波有限差分数值模拟方面，差分网格及边界条件是影响弹性波模拟成功与否的关键。从各向同性介质速度一应力方程出发，利用交错网格高阶有限差分和完全匹配层（PML）边界条件，提出了各向同性介质弹性波交错网格有限差分正演模拟方法。数值试验结果表明，该方法精度较高，数值频散较小，人为边界反射吸收较好，为后续的地震属性分析（如AVO／AVA）奠定基础。     

声波正演中一种新的边界条件——双重吸收边界条件

在声波波动方程有限差分波场数值模拟中，由于计算模型的限制，导致很强的人为边界反射的产生。为了使计算得到的边界波场值更接近真实的边界波场值，往往使用边界条件以减少来自计算区域边界的人为反射能量，使得声波在有限区域传播达到在无限介质中传播的效果。目前，在声波正演中，较常用的一种边界条件是基于标量波动方程旁轴近似理论的吸收边界条件，在此基础上诞生了许多新的边界条件算法。在前人成果的基础上，提出了一种新的吸收边界条件——双重吸收边界条件。该方法将标量波动方程旁轴近似与另一种边界方程通过加权的方式，使所计算的边界波场值比利用通常吸收边界方程计算的边界波场值更接近真实的波场。该方法计算量小，简单易实现，应用于声波正演能得到比利用常规吸收边界条件更好的吸收效果。     

弹性波数值模拟中的吸收边界条件

在地震波数值模拟过程中,需要采用吸收边界条件以吸收入为边界反射。本文通过对ＴＩ介质中二维及三维一阶弹性波方程的特性分析,分别得到了与原波动方程形式相同的吸收边界条件,并在二维弹性波数值模拟中进行了应用。     

优化系数的吸收边界条件方程

在波动方程有限差分的正演及偏移成像过程中，由于引入人为的边界条件，必然会导致产生边界反射，从而影响正演结果及偏移成像的效果。为此，人们提出了各种吸收边界条件方程，但其吸收效果均不理想。本文在最小二乘意义下，用共轭梯度导出了优化求解吸收边界条件方程中的系数，使得吸收边界条件方程尽可能与内部方程相逼近，从而使得边界反射波被全部或部分吸收掉。理论分析和偏移成像实践表明，优化系数的吸收边界条件方程对以各种     

声波数值模拟中改进的非分裂式PML边界条件

如今完全匹配层(PML)吸收边界条件被认为是最好的吸收边界条件,但常规PML吸收边界条件存在变量个数多、计算存储量大、编程复杂以及场源靠近截断边界时吸收性能有所下降等缺点。本文推导出一种适用于声波方程改进的非分裂式(NPML)吸收边界条件。该方法基于非分裂式PML吸收边界条件,采用复频移伸展函数,将PML介质中的波动方程分解为正常项和衰减项并单独进行迭代更新。分析和试验结果表明该方法应用于声波数值模拟时不但具有优良的吸收衰减性能,而且有计算方程简单、编程实现容易、PML辅助变量占内存小等优点。     

PML边界条件下二维粘弹性介质波场模拟

在波场模拟中较多使用的有限差分法和有限元法都或多或少地存在一些缺陷，如为了提高计算精度而使运算效率降低，在泊松比变化大的界面存在稳定性问题等。伪谱法是一种计算精度高且计算效率也较高的方法，它是对空间坐标通过快速傅里叶变换在时间域作差分运算，避免了求偏导数；不存在有限差分法和有限元法对高频成分限制的问题，可以实现全频带地震波模拟；对内存容量的需求也远远低于有限差分法和有限元法。为此，采用了伪谱法进行粘弹性介质波场模拟。首先给出了开尔芬粘弹性介质中的波动方程，推导出二维粘弹性介质最佳匹配层（PML）吸收边界条件及其相应的伪谱法计算公式；然后对均匀介质模型进行了模拟，结果表明，该方法的边界吸收效果很好；最后通过数值模拟分析了具有不同粘滞系数介质对地震波的吸收和衰减，结果表明，随着粘滞系数的增大，粘弹性反射波的主频向低频方向移动，高频吸收明显，有效频带变窄，振幅降低。此外，在伪谱法中引入匹配层边界条件后，只需对空间做一维傅里叶变换，大大提高了伪谱法的计算效率。     

基于CPML-RML组合边界条件粘弹TTI介质旋转交错网格有限差分正演模拟

实际地下介质中因充满流体、裂缝和裂隙等而同时表现出粘滞性特征和各向异性特征,传统的分裂式完全匹配层(SPML)吸收边界条件不能有效吸收低频和大角度入射波。为此,提出将非分裂式卷积完全匹配层(CPML)吸收边界条件与随机介质层(RML)边界条件相结合的策略来改善边界吸收效果,并将其应用于粘弹TTI介质有限差分正演模拟。为提高正演模拟的稳定性和精度,采用了旋转交错网格有限差分算法,并推导出了适用于粘弹TTI介质的CPML吸收边界条件公式及相应的旋转交错网格有限差分格式;采用随机介质建模理论并结合CPML吸收边界条件,建立了CPML-RML组合边界条件。利用二维均匀介质模型和二维复杂Hess模型对组合边界条件的吸收效果进行了测试,结果表明:相比SPML吸收边界条件与CPML吸收边界条件,CPML-RML组合边界条件在不增加计算量的情况下,具有更好的吸收效果,能有效减少人工边界反射、提高数值模拟精度;粘弹TTI介质中的地震波场表现出明显的振幅衰减和相位延迟,地震记录表现出浅部能量强、深部能量弱的特征,验证了CPML-RML组合边界条件对复杂介质的适应性及正确性。     

一阶速度—压力常分数阶黏滞声波方程及其数值模拟

与传统的整数阶黏滞波动方程相比,分数阶拉普拉斯算子黏滞方程能更准确地匹配目前广泛使用的常Q模型,而且分数阶黏滞波动方程中控制振幅衰减和相位变化的算子是显式分离的,这对于发展稳定的衰减补偿逆时偏移算法至关重要。首先基于时间域二阶位移形式的常分数阶拉普拉斯算子黏滞声波方程,推导了一阶速度-压力形式常分数阶拉普拉斯算子黏滞声波方程;为了模拟更加真实的振幅变化信息,在新的黏滞声波方程中考虑了密度空变的影响;为了避免由傅里叶变换的周期性而引入的虚假反射,提出了一种适用于分数阶黏滞声波方程的卷积型完全匹配层（CPML）吸收边界加载方法;最后采用交错网格伪谱法进行数值模拟。均匀介质中数值解与解析解的对比证实了该一阶速度-压力常分数阶黏滞声波方程能准确描述常Q模型,BP盐丘模型的地震波场模拟结果证实了其对复杂介质的适用性。     

多轴复频移近似完全匹配层弹性波模拟

在地震正演模拟中完全匹配层是现今吸收效果最好、应用最广的边界条件。由分裂法发展到各种非分裂方法,很大程度上提升了完全匹配层的计算效率并减少所需内存。近似完全匹配层是一种直接对波场进行变换的新型非分裂方法,相比卷积完全匹配层具有不改变方程形式且易于实现等优势,但无法吸收掠入波和瞬逝波,因此目前仍以卷积完全匹配层为主。本文在近似完全匹配层边界处理技术的基础上引入复频移技术,实现了对大角度入射波和极低频波的吸收;采用多轴技术提升近似完全匹配层的稳定性;引入新的衰减函数以整体提升近似完全匹配层吸收效果。最终使近似完全匹配层能够更稳定有效地吸收任意入射波。     

VTI介质频率——空间域准P波正演模拟

本文从VTI介质弹性波动方程出发，借助VTI介质弹性参数和Thomsen参数，结合Kelvin-Christorffel方程，推导了VTI介质中准P波动方程，并对VTI介质准P波进行了正演模拟。在正演模拟中，为了克服常规差分算子的数值频散，采用了25点优化差分算子；再依据最优化理论求取的优化系数建立了频率-空间域中准P波波动方程的差分格式；为了消除人为边界反射，根据特征分析方法并利用Kelvin-Christoffel方程，构建了VTI介质中准P波方程在不同边界和角点处的边界条件，再由准P波波动方程和边界条件，通过频率-空间域有限差分法，对准P波在均匀VTI介质，层状VTI介质和回陷模型中的传播过程进行了数值正演模拟。通过正演模拟，得到了单频波波场，时间切片和共炮点记录，为研究地震成像及反演等提供了依据。     

钻柱纵向振动仿真分析

国仙外在研究钻柱振动时都把钻柱下端当作铰支,假设激振源为周期性的,偏离实际。从建立符合实际的下端边界条件入手,着重研究了钻柱的纵向振动。用有限单元法,求出整体质量、整体刚度,整体阻尼矩阵采用的是比例阻尼矩阵,利用牛顿法得到钻柱的振动方程,另外通过钻头与岩石的相互作用模型得到了起初的钻柱下端边界条件,并利用Ｈｏｕｂｏｌｔ数字解法模型求解了振动微分方程。根据上述模型编制了用于分析钻柱纵向振动的计算机仿     

用声波全波列评价薄互储集层

针对孔隙地层薄层声波测井问题,通过声压－速度有限差分数值模拟,分析了不同性质、不同厚度的水平薄层和薄互层在声波全波列波形中引起的不同特征变化。全波列各组分波在任意分层界面处都会产生反射波,其中低频斯通利波反射最为突出,反射波的幅度主要取决于薄层与围岩地层之间的声阻抗差异。当薄层性质由软地层变为水层时,反射波幅度增大,纵波幅度衰减严重。通过波形滤波处理,可根据各组分波面上下的时差,幅度变化及斯通利波的反射系数变化,分辨薄互层的位置和大致厚度。     

基于优化时空域频散关系的声波方程有限差分最小二乘逆时偏移

相比常规逆时偏移,最小二乘逆时偏移（LSRTM）成像结果更趋真实,分辨率更高,其精确成像的关键之一是波动方程的精确、高效求解。为了提高LSRTM计算精度,本文通过优化时空域频散关系求取差分系数,进行有限差分正演模拟,频散分析和数值模拟结果表明采用该算法可提高数值模拟精度;利用混合吸收边界条件压制边界反射,可取得较好吸收效果;采用改进的共轭梯度法计算流程,能减少波动方程正演次数,提高计算效率。模型数据测试及分析表明,本文方法能有效提高LSRTM成像精度,精细刻画构造细节,提高迭代收敛速度。     

南马庄潜山内幕反射空校作图方法及效果

南马庄地区第三系地层埋藏深度较浅,且直接覆盖在下古生界灰岩(潜山)之上.由于潜山面上覆介质的速度是深度的线性函数,而潜山面以下(潜山内幕)介质的速度近于常数,所以对于潜山内幕的空校就不能沿用常规的空校方法.     

多道联合约束的匹配追踪强反射轴压制方法

针对强反射轴中及其附近有效信号难以识别和区分的问题,利用匹配追踪算法的稀疏分解表示能力,提出一种多道联合约束的匹配追踪(TC-MMP)强反射轴压制方法.对于多道匹配追踪算法分解信号过程中迭代终止条件难以确定的问题,引入残差比阈值,与固定迭代次数或固定残差门限相比,在一定程度上提高了信号分解的效率,由此形成了迭代条件约束...     

非均匀S偏振电磁波在导电媒质界面的反射系数曲线

基于非均匀电磁波在导电媒质中传播时其相移常数方向和振幅衰减常数方向的不一致性，利用电磁波在导电媒质界面的边界条件给出了非均匀S偏振电磁波的反射系数，并计算了反射系数随入射角的变化曲线。曲线显示当电磁波由电导率大的介质向电导率小的介质传播时，将发生类全反射。在相移常数临界角和衰减常数临界角附近存在2个峰值，且反射系数的数值小于1，这一结果与全反射光的结果明显不同。