弹性波叠前逆时深度偏移技术及其应用

研究了二维各向同性介质中两分量记录的基于弹性波标量波动方程的叠前逆时深度偏移方法.从弹性波波动方程出发,推导各向同性介质中弹性波逆时深度偏移技术中相关算法,并探讨其中有关问题的解决方案,实现二维多分量资料的叠前逆时深度偏移,最后进行模型试算和实际资料试处理.结果表明,叠前逆时深度偏移方法不受地层倾角限制,较适用于高陡构造地区或介质横向速度变化较大地区的多分量地震资料处理.     

各向异性介质中弹性波多分量联合逆时深刻偏移

当今石油工业界使用的大部分偏移算子都是在假设地球为均匀各向同性介质的前提下经过一定的近似得到的，其应用了一定的局限性。因为经过大量的研究已经证实地壳和上地幔是各向异性介质，不考虑介质各向异性的偏移算子必然带来一些在反射点归位方面的不可估计的错误。从弹性波动方程出发，推导出了适用于各向异性介质中进行弹必一波多分量联合偏移算子，并用数据实例进行证明。     

基于黏弹各向异性理论的地震波场数值模拟与分析

基于Carcione黏弹各向异性理论及具有任意倾斜对称轴的横向各向同性介质(TTI)的本构关系,利用二维三分量高阶交错网格有限差分和完全匹配层边界条件对黏弹TTI介质的地震波场信息进行了数值模拟。研究发现:1改变黏弹TTI介质的弹性参数,对地震波传播造成显著影响。其中,C_(13)主要控制qS1波中质点的偏振方向,进而影响了qS1波波形和反射同相轴的形态;而C_(66)则主要控制qS2波中质点的偏振方向,从而影响了qS2波波形和反射同相轴的形态;2利用层状黏弹TTI介质进行数值模拟,可以很好地观察到拟横波在遇到反射界面时发生的横波再分裂现象;同时由于介质黏弹性的存在,导致地震波出现显著的衰减。这将为进一步研究各向异性介质中复杂的波场信息提供帮助。     

用邦德矩阵法推导各向异性介质的弹性常数

弹性波各向异性介质是描述地球内部物质的比较符合实际的一种模型,而其弹性常数的个数是研究地震波传播特性的关键所在。横向各向同性介质（ＴＩ）及裂隙各向异性介质（ＥＤＡ）均高为具无穷次对称轴的介质,在此假设基础上,利用邦德矩阵法能够方便快捷地推导出具有五个相互独立常数的弹性刚度矩阵。     

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

基于Biot理论,给出了三维双相各向异性介质应力-速度弹性波方程交错网格任意偶阶精度有限差分解法,并对三维双相各向异性(TIH)介质中弹性波场进行了模拟.结果表明,弹性波在三维双相各向异性介质中传播时存在快纵波qP1、快横波qS1、慢横波qS2和慢纵波qP2,并清楚地观测到了横波分裂、横波分裂盲点、波面三分叉等特殊现象.快纵波在固相和流相中的相位相同;而慢纵波在固相和流相中的相位相反,即慢纵波在流相中振幅大,而在固相中的振幅较小.另外,三维双相各向异性介质中弹性波场的快纵波和快横波的耦合关系、波的类型、能量分布和相位等都是在三维空间中变化的.     

用有限单元法模拟各向异性介质中二维电性异常体的大地电磁响应

为了研究电导率各向异性介质中二维大电磁场的变化规律,首先从其满足的变分问题出发,利用有限单元法计算出大地电磁场,借助辅助场获得视电阻率,其次讨论TEz模式下各向异性介质和各向同性介质中视电阻率响应曲线的区别,接着探讨不同层面倾角、不同频率下各向异性异常体引起大地电磁场的变化,另外就高、低阻各向异性介质对大地电磁视电阻率影响的差异进行讨论,最后讨论平行层面和垂直层面的电阻率对各向异性介质平均电阻率的影响,从而获得各向异性大地电磁场响应特征。     

三维多分量资料处理关键技术进展

与纵波地震资料相比,由于多分量资料的矢量特征、各向异性和转换波传播路径非对称性,因此需要一套特殊的资料处理技术,解决多波地震资料处理中的难题,同时发挥矢量波场的优势,从而改善三分量地震资料处理的质量.利用多分量的矢量特征、偏振特征以及面波特征,压制随机和相干噪声,进行纵波和横波的波场合成与分离.针对转换波非对称传播路径问题,提出转换波共反射面元计算、建立近地表浅层横波速度模型的解决方案,从而解决转换波共反射面元和转换波静校正问题.各向异性对转换波的影响比纵波更加明显,因此如何从多波资料中反演各向异性参数,消除各向异性参数对成像质量的影响至关重要.在这个处理环节中,建立了4参数各向异性反演方法,实现转换波各向异性动校正和各向异性叠前时间偏移处理.同时在国内最早提出和实现了转换波弯曲射线各向异性叠前时间偏移技术,改善偏移能量的聚焦性,降低成像剖面低频化现象.     

横观各向同性介质中弹性波的交错格式差分解法

基于应力、速度混合变量弹性波方程及广义胡克定律，给出了一种横观各向同性介质中弹性波计算的高精度交错格式差分解法；本方法可适用于高泊松比材料，数值稳定性较好．均匀和非均匀横观各向同性介质中弹性波传播的数值模拟表明，本文方法计算耗时少、精度高，可有效地模拟复杂各向异性介质中弹性波的传播     

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

基于Biot理论，给出了三维双相各向异性介质应力一速度弹性波方程交错网格任意偶阶精度有限差分解法，并对三维双相各向异性（TIH）介质中弹性波场进行了模拟。结果表明，弹性波在三维双相各向异性介质中传播时存在快纵波qP1、快横波qS1、慢横波qS2和慢纵波qP2，并清楚地观测到。T横波分裂、横波分裂盲点、波面三分叉等特殊现象。快纵波在固相和流相中的相位相同；而慢纵波在固相和流相中的相位相反，即慢纵波在流相中振幅大，而在固相中的振幅较小。另外，三维双相各向异性介质中弹性波场的快纵波和快横波的耦合关系、波的类型、能量分布和相位等都是在三维空间中变化的。     

各向异笥二相介质中位错环的弹性场

将无奶大各向异性介质中无限小位错环弹性理论另以推广，结合运用Ｇｒｅｅｎ张量函数积分法提出计算各向异性二相介质中有限位弹性场的处理方案，并以γ／α－Ｆｅ二相系统中的有限位错不为例，对其在介所产生的总弹性场及计算，结果出这处位钽不构弹性场的复杂性，相界对位错环有吸引作用，所提出的理论方案可推广应用于任意复相系统中任意形状有限位错环弹性场及其所受像力的计算，所得结果对探讨位错环在实际材料和相变中的作用具     

横向各向同性介质中的地震波传播特性

在各向异性介质中 ,地震波的传播特性随传播方向而变化 ,常表现为地震波传播的相角和群角、相速度和群速度的不一致性。应用理论模型计算方法研究了具有垂直对称轴的横向各向同性介质中P波、SV波的传播规律和Thomsen参数对P波相速度的影响。计算表明 ,横向各向同性介质中的P波与SV波的相速度曲线不再是各向同性介质中的圆形。就相速度而言 ,在垂直方向附近 ,SV波的速度变化要比P波明显得多 ;Thomsen参数vS0 对P波相速度的影响非常小。研究结果对各向异性介质中的弹性参数反演、速度反演、NMO、DMO及时间域的偏移有着重要意义。     

粘弹性波动方程正演和偏移

理论和实际研究均表明大地介质属于非完全弹性介质,它是具有粘滞性的粘弹性体,地震波在大地中传播时,质点的震动能量转化为其他各种形式的能量,这种能量的转化使地震波在粘弹性介质中传播时,其高频成分很容易被吸收,振幅近似按指数规律衰减,波形和相位失真.作者从二维粘弹性波动方程出发,在频率波数域内导出了粘弹性波动方程的正演和偏移算法.根据设计的地质模型,分别对粘弹性介质和弹性介质中的点状绕射源进行正演模拟,模拟结果表明介质的粘滞性对地震波予以吸收和衰减;将粘滞性介质中点状绕射源的正演结果,分别用粘弹性波动方程偏移算法和弹性波动方程偏移算法进行偏移模拟,并予以对比分析,其结果进一步说明了对实际介质中地震波的吸收和衰减作用的处理效果,证明了所得算法的准确性和有效性.     

任意裂隙方位双相HTI介质的地震波方程及正演模拟

基于BISQ方程的高频极限,导出了适用于任意方位角的双相HTI介质波动方程,推导了该方程在交错网格空间中求解的高阶有限差分格式和对应的完全匹配层(PML)吸收边界条件,实现了该类介质的地震波场正演模拟。模拟结果表明,该方法能准确模拟地震波在双相HTI介质中的传播过程,得到高精度的波场快照和合成记录。分析得出:在双相各向异性介质中,弹性波的传播仍表现出各向异性特征;孔隙度、粘滞系数和固液耦合密度主要影响慢纵波的振幅和形态。     

各向异性介质吸收边界条件

针对地震波勘探数值模拟中遇到的人为边界反射问题,提出了适用各种各向异性介质的吸收边界条件,证明了其稳定性,在频率域推出了反射系数算法。最后将公式进行差分离散,并对各向同性介质,TI介质及带方位人射的EDA介质进行了数值计算。     

拟地震法电磁波场有限元法正演模型

在电法勘探正演模拟中,利用有限元法提出了电磁波在二维电阻率介质模型中的传播特征。我们将电磁波波动方程与声波波动方程类比,运用适当的近似,取等效方程中相同的主要项,利用适当的边界条件,则电磁波场问题可化为类似于地震波场问题,从而可用地震资料处理技术中的任何方法去完成。我们选用的是波动方程有限元方法。计算可在频域也可在时域进行去模拟二维电阻率分布的响应即电场视剖面。实例表明:正演结果所得到的二维电阻率成像剖面与实际的二维电阻率分布完全相吻合。结果表明利用有限元进行电阻率成像是切实可行的。     

无反射波动方程模拟堤防复杂构造的波传问题

为了更好地了解浅层复杂介质中的波场特征 ,以及更好解释野外堤防的探测资料 ,该文利用无反射波动方程和交错网格的有限差分法研究了堤防复杂构造中波传问题 ;并结合实际堤防中物性参数 ,模拟出自激自收地震剖面 ,有效地避免了层间的多次反射。讨论了无反射波动方程的特性 ,推导出多次反射满足的波动方程 ,从特征线分析的角度探讨了无反射波动方程采用交错网格的波传特性和数值计算的稳定性 ,以及利用套网格技术解决源点奇异性的有效性。     

弥散黏滞性波动方程的吸收边界算法

针对地震波数值模拟中因计算区域有界而产生的人工边界反射问题,提出了一种弥散黏滞性波动方程的吸收边界算法.在所研究的有界计算区域周围加入合适的吸收层,使得吸收层内的地震波随传播距离指数衰减而达到吸收边界反射的目的.利用傅里叶变换求得无界空间该波动方程的解；引入衰减函数构建有界空间相应的辅助方程,合理地选择衰减函数使得计算区域内辅助方程的解为原方程的解,而边界吸收层内的解呈指数衰减；运用有限差分法在均匀介质和均匀层状介质中求解弥散黏滞性波动方程,采用该方法对边界进行处理并与未加边界处理的波场快照和地震记录进行对比.实验结果表明,所提方法处理后的波场快照几乎无边界反射波存在,在空间位置为x=147 m、z=747 m的地震记录中,0.5s处边界反射波的振幅趋于0.所提方法也适用于声波方程和Stokes方程.     

不同边界条件下的二维声波方程数值模拟研究

地震波场数值模拟是研究波动现象的重要手段之一,它对于油气田的勘探和开发具有重要意义。数值模拟过程中,需要通过添加边界条件来尽可能消除由于截断所产生的边界反射。本文选取雷克子波作为震源项,分别建立均匀及层状地质模型,拟定合适的波场模拟参数,实现了不同边界条件下的二维声波方程数值模拟。利用数值模拟得到的波场快照和地震记录直观地对比分析不同边界条件对边界反射的消除效果,本文认为透明边界条件(以下简称TBC吸收边界条件)和Clayton-Engquist边界条件(以下简称CE吸收边界条件)都能够较好地消除边界反射。最后,本文提出了一种组合边界条件的方法。     

基于波动方程的Ansys模型黏弹性边界条件研究

 研究了应用Ansys 有限元软件进行地震波场数值模拟的边界条件问题,提出一种基于波动方程的Ansys 模型黏弹性边界条件构建方法。由波动方程推导并建立了黏弹性边界条件的理论基础,在Ansys 模型中选用带阻尼的Combin14 单元作为加载有限元,采用APDL 代码成功地将黏弹性边界条件施加于Ansys 模型中。对模型进行实例计算,在模型边界设置虚拟检波器,加载Ricker 子波震源,提取计算后的波场快照和时间记录进行对比,结果显示,所施加的黏弹性边界条件可以很好地吸收边界反射波,表明了该方法所构建的Ansys 模型黏弹性边界条件的有效性。