VTI介质中Thomsen参数对相速度和群速度的影响

相速度和群速度对于研究弹性波在各向异性介质中的传播具有重要意义,同时对合理评价地震资料解释成果也具有一定的指导意义。深入研究Thomsen参数与群、相速度的关系有助于研究复杂地层中波的传播规律以及速度的变化趋势,从而使地震资料解释更加准确。本文首先简要推导了VTI介质中弹性波的群、相速度解析表达式,然后通过编程实现,直观显示了三维速度曲面的变化特征与Thomsen参数的联系,并通过二维及一维速度曲线变化加以验证。实验结果表明：VTI介质中群速度较相速度对Thomsen参数的变化更加敏感,其中随着 (qP波各向异性强度表征参数)与 (描述纵横波相速度的过渡性参数)差值的不断增大,各向异性程度明显增强,群速度剧烈变化,尤其是qSV波;随着 (横波分裂强度表征参数)值的不断增大,SH波群、相速度面的扁率也不断增大。     

基于HTI介质理论的纵波速度分析方法

随着地震勘探精度的不断提高,仅考虑VTI介质各向异性速度分析方法已不能满足宽方位地震资料速度分析和动校正的需要.针对宽方位纵波资料方位各向异性特点,讨论了HTI介质各向异性速度分析方法.根据方位各向异性理论,得出了速度和裂缝方位满足的椭圆方程.采用最小二乘拟合方法,求出了椭圆方程的参数解.根据椭圆方程求出每个方位的速度...     

垂直对称轴的横向各向同性介质中Thomsen参数对相速度和群速度的影响

相速度和群速度对于研究弹性波在各向异性介质中的传播具有重要意义;同时对合理评价地震资料解释成果也具有一定的指导意义。深入研究Thomsen参数与群、相速度的关系有助于研究复杂地层中波的传播规律以及速度的变化趋势,从而使地震资料解释更加准确。首先简要推导了各向同性(VTI)介质中弹性波的群、相速度解析表达式,然后通过编程实现,直观显示了三维速度曲面的变化特征与Thomsen参数的联系;并通过二维及一维速度曲线变化加以验证。实验结果表明:VTI介质中群速度较相速度对Thomsen参数的变化更加敏感;其中随着ε(q P波各向异性强度表征参数)与δ(描述纵横波相速度的过渡性参数)差值的不断增大,各向异性程度明显增强,群速度剧烈变化;尤其是q SV波;随着γ(横波分裂强度表征参数)值的不断增大,SH波群、相速度面的扁率也不断增大。     

裂隙介质中地震波传播特征近似

地震波传播的精确表达式能够精确表达地震波在介质中的传播特征,但由于表达式较复杂,为一些参数隐性表达式,这不利于参数反演,因此有必要研究其近似的线性表达式,以便进行参数反演,准确探测地层结构。在扰动法求解Christoffel方程,得到弱各向异性条件下qP、qS相速度、群速度近似表达式的基础上,依据介质模型坐标与观测坐标,各种介质弹性系数矩阵关系式推导出裂隙介质中地震波相速度、偏振向量及群速度的近似表达。数值计算结果表明,当各向异性系数在0.3以内时,对一般裂隙各向异性介质都适用,且具有较高的精度,在各向异性系数大于0.3时,纵波(P波)仍然具有较高的精度,但横波(S波)的误差增大到50%,主要是因为扰动解是在弱各向异性条件下推导的,因此在介质各向异性系数较小的情况下,为地震波速度分析及参数反演奠定了基础。     

声学各向异性介质中的S波

声学VTI介质是人为定义横波在对称轴方向的速度V_S0的一种模型。传统理论认为当V_S0=0时可以消除横波能量。论证声学VTI介质的存在,并讨论其特殊的性质,可以看到在声学VTI介质中SV波的相速度和群速度并不是恒为零。用高阶交错网格有限差分方法对不同的声学VTI介质模型进行全波场模拟,同时计算V_s0取不同值时的群速度和相速度x、z分量。通过对比分析得出波现象：SV波在垂直方向传播时其极化方向是平行于P波传播方向的;由于邻近方向能量的聚集,在对称轴附近SV波的振幅会有较大的值;P波和SV波的群速度的值在某些方向上可以很接近甚至相等。     

声学各向异性介质中的S波

声学VTI介质是人为定义横波在对称轴方向的速度V_S0的一种模型。传统理论认为当V_S0=0时可以消除横波能量。论证声学VTI介质的存在,并讨论其特殊的性质,可以看到在声学VTI介质中SV波的相速度和群速度并不是恒为零。用高阶交错网格有限差分方法对不同的声学VTI介质模型进行全波场模拟,同时计算V_s0取不同值时的群速度和相速度x、z分量。通过对比分析得出波现象：SV波在垂直方向传播时其极化方向是平行于P波传播方向的;由于邻近方向能量的聚集,在对称轴附近SV波的振幅会有较大的值;P波和SV波的群速度的值在某些方向上可以很接近甚至相等。     

黏声TTI介质中纯qP波逆时偏移成像方法

地下介质广泛存在黏滞性及各向异性特征,使得地震波在传播过程中发生相位畸变和振幅衰减。若在偏移成像时不对这些黏滞性及各向异性影响进行校正,则会引入成像噪声而降低成像分辨率。针对上述问题,首先将黏滞性引入到TTI介质波动方程,得到一种黏声TTI介质纯qP波波动方程。然后,基于新推导的黏声TTI介质波动方程,发展相应的逆时偏移成像方法。由于黏滞性补偿过程中能量会被放大,波场中指数增加的高频噪声易造成波场传播不稳定。鉴于此,将正则化算子引入到逆时能量补偿过程,以压制波场中的高频噪声。结果表明：新提出的黏声TTI介质逆时偏移成像算法可校正相位畸变,补偿能量衰减,实现对含黏滞性及各向异性介质的高精度成像。     

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

三维数值模拟是研究各向异性介质中复杂弹性波波场特征和传播规律的重要手段。根据Taylor展开式,推导出了交错网格一阶空间导数的任意偶数阶精度有限差分近似式和相应的差分系数计算式,给出了三维各向异性介质中弹性波一阶双曲型应力速度方程交错网格任意偶数阶精度差分格式和稳定性条件,并推导出了三维各向异性介质PML法吸收层边界条件公式和相应的交错网格差分格式。对方位各向异性介质模型和正交各向异性介质模型中弹性波的传播进行的三维数值模拟结果表明,弹性波在三维各向异性介质中传播时存在拟P波、拟SV波和拟SH波,并出现横波分裂、横波分裂盲区、波面三分叉等特殊现象,另外,弹性波场在空间是变化的,其拟P波、拟SV波和拟SH波的耦合关系比较复杂。     

基于Chebyshev配置点谱方法的多孔介质平板通道内的流体流动数值模拟

本文基于Chebyshev配置点谱方法,利用数值模拟研究了多孔介质平板通道内的流体流动问题。针对动量方程的离散,空间上采用Chebyshev配置点谱方法,时间上采用准隐式格式离散,结合改进的投影算法(IPS)将速度和压力的计算解耦为一系列椭圆方程(泊松方程或亥姆霍兹方程),转换椭圆方程为矩阵方程形式后利用二步求解法求解矩阵方程。通过MATLAB编程实现对多孔介质平板通道内的流体流动问题的数值模拟并验证了程序的准确性。在此基础上,讨论了达西数(Da),雷诺数(Re)以及孔隙率(ε)对多孔介质平板通道内流体的速度分布、边界层厚度及入口长度的影响。     

椭圆高斯光束在强非局域非线性介质中传输特性的变分研究

运用变分法研究傍轴椭圆高斯光束在强非局域非线性介质中的传输特性,得到了光束各参量的演化方程和各参量的演化规律的近似解、2个临界功率和2个由强非局域介质特性与初始束宽决定的特征参量ε1和ε2。在0ε1<1/3,0ε2<1/3的前提下,当初始功率等于某一临界功率时,傍轴椭圆高斯光束在传输时,可以保持某一方向的束宽不变,得到空间光孤子。     

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

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

垂直发育裂隙介质中PP波扰动法近似反射系数研究

为了实现裂隙介质中的地震AVO(amplitude variation with offset)反演,基于Christoffel方程和边界条件,推导了EDA介质中PP波、PS1波、PS2波精确反射系数表达式,并采用扰动法推导了极端弱各向异性介质中PP波的近似反射系数,通过弹性系数简化得到HTI、裂隙EDA介质中PP波的近似反射系数。模型计算结果表明:在弱各向异性条件下,HTI(horizontal transverse isotropy)介质中地震波的近似反射系数公式计算精度较高,相对误差在4%以内;当裂隙介质在弱各向异性(各向异性系数0.1)的情况下,当入射角在40°以内时,PP波近似反射系数与精确反射系数的绝对误差在10~(-4)以内,相对误差小于4%,但随着入射角的加大其计算误差有所增加。对各向异性较强烈(各向异性系数达到0.2)的裂隙介质,当入射角相对较小时(小于45°时),裂隙介质中PP波近似反射系数计算式对于强各向异性的介质仍然成立。通过对反射系数的近似研究,可以将裂隙介质中反射系数的非线性问题转为线性问题,进而利用这些特性进行参数反演,有利于提高反演速度。     

线电荷在无限长矩形腔内激发的电势

线电荷在充满电各向异性介质的无限长矩形腔中激发的电势,是各向异性介质中泊松方程的边值问题;主要应用本征函数法和分离变量法求解了边界条件为第一类齐次和第一类非齐次边界条件时,线电荷在无限长矩形腔内激发的电势;在令ε_(11)=ε_(22)=ε_(33)=ε的情况下,所得的结果可适用于电各向同性介质。     

利用体积分方程的弱化BCGS算法模拟各向异性地层中多分量低频海洋可控源电磁响应

基于层状单轴各向异性介质的并矢Green函数和体积分方程的弱化稳定型双共轭梯度(BCGS)算法模拟各向异性背景地层中多分量低频海洋可控源电磁响应,并通过算例分析发射频率、海底背景地层各向异性和异常体埋藏深度对响应的影响特征。结果表明：发射信号的不同频率、电场强度的不同分量对海底油气藏的水平分辨能力不同;当测线沿x方向布设时,E_zz分量的分辨能力最强,E_xx分量次之,E_yy分量的水平分辨能力最低;在海底背景地层水平电导率相同的情况下,背景地层垂向电导率的变化对E_xx分量和E_zz分量的影响规律不同,而E_yy分量基本不受地层垂向电导率变化的影响;异常体埋藏深度越浅,E_xx分量和E_zz分量对异常体的水平分辨能力越高,而无论异常体埋藏深浅E_yy分量均不能对其进行分辨。     

各向异性固体弹性介质中群速度的解析表达式

通过引入一个仅仅与介质弹性常数有关的矩阵M ,将Christoffel系数矩阵分解成M和一个与波传播方向有关的矢量 P两个部分的乘积 .由此推导出适用于计算各种各向异性固体弹性介质弹性波群速度的普适解析表达式 .群速度作为弹性波波矢的函数 .在特殊对称性方向上该解析表达式的结果与已有文献的结果一致 .利用Matlab软件数值求解Chritoffel方程 ,得到对应于任意方向上的相速度 .图 2 ,参 7.     

含直立裂缝粘弹性介质地震波场正演模拟

以Carcione各向异性粘弹性理论为基础,导出了适用于任意方位角的粘弹性方位各向异性(HTI)介质的波动方程,推导了该方程在交错网格空间中求解的高阶有限差分格式和对应的完全匹配层(PML)吸收边界条件,实现了该类介质的地震波场正演模拟。模拟结果表明,该方法能准确模拟地震波在含直立裂缝粘弹性介质中的传播过程,得到高精度的波场快照和合成记录。初步分析了介质的各向异性特征和衰减特征,研究了裂缝方位角和粘滞性参数对地震波场的影响机理。     

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

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

各向异性介质中P波反射系数

分析了横向各向同性和方位各向异性介质的成因及其本构关系。由此讨论弹性波在这两种各向异性介质中的传播特点,并提出了表征各向异性程度的广义各向异性参数,这些参数既可表征横向各向同性介质的各向异性,也可表征方位各向异性介质的各向异性。讨论了两种各向异性介质中存在水平界面时的反射系数,并将横向各向同性介质中的反射系数公式推广到方位各向异性主轴上。根据算例修正了近似式,着重讨论和分析了这两种各向异性介质中的地震波反射振幅与炮检距的关系（ａｍｐｌｉｔｕｄｅｖｅｒｓｕｓｏｆｓｅｔ,ＡＶＯ）以及对实际勘探的影响和指导意义。     

致密储层各向异性地震岩石物理建模及应用

针对致密砂岩储层低孔低渗及微裂缝发育的特点,提出一种适合致密孔隙裂缝型储层的岩石物理模型构建方法:首先,在岩石基质中利用有效介质理论添加含流体孔隙,然后利用结合Hudson理论和各向异性Gassmann理论推导得到的各向异性流体替代方程添加含流体裂缝。基于构建的岩石物理模型,定量分析裂缝密度、裂缝形状、充填流体等储层岩石的物理性质对饱含流体岩石的纵横波速度的影响。采用某研究区实际井资料对该方法进行试算,完成纵横波速度的估算和裂缝储层各向异性参数的提取,纵横波速度的估算值与其实测值的吻合程度较高,且各向异性参数能够较好地反映裂缝发育位置,证明了本文方法的有效性。     

扰动KdV方程解的先验估计

对于带微扰的KdV方程u_t+6uu_x+u_xx=εR(u)，(ε＞0)，在初值u₀(x)∈C∞(-∞,+∞)，当|x|→∞时指数衰减的条件下，分别构造出带两种不同扰动项的KdV方程的扰动孤立波解满足的能量关系式，并运用能量分析方法对扰动的孤立波解进行先验估计，得到如下结论：(1)R(u)=δ(εt)u, δ(s)∈C［0,+∞)，δ(0)=0，时，解在-∞＜x＜+∞，0≤εt≤T内一致有界；(2)R(u)=-Δ(εt)u_xxx，Δ(0)=0，Δ(s)∈C¹［0,+∞), 解在-∞＜x＜+∞，0≤εt≤T，0≤ε≤ε₁内一致有界。