声波测井时裸眼井中弹性波的传播

本文分析了简化模型下声波测井时裸眼井（内充流体一般是泥浆）中弹性波的传播,此波场可看成是轴向传播有不同波矢的分波的叠加。分析表明,对于井外只有单一均匀岩层且无限深的井,仅对如（18）式所示的特征频率成分才存在具有轴向传播最大相速度v_zm=c₂的所谓“首波”;而v_x > c₂及c₂ > v_z>b₂的分波却不存在,此处c₂和b₂分别是井外岩层的纵波和横波速度。这正解释了文献[3]中所报道的裸眼井的“过滤效应”。分析还表明,对于井外具有多层岩层的井,只在一定条件下才存在与上述特性相近的“首波”。本文结果对声速测井显然是重要的。     

横向各向同性地层斜井中正交偶极子激发声场的数值模拟

本文采用三维应力-速度有限差分(SV-FD)方法,数值模拟了横向各向同性(TI)地层对称主轴与井轴斜交情况下正交偶极子声源激发的井孔声场.主要解决了与倾斜角有关的三维空间弹性模量矩阵的推导,柱坐标系下应力-速度有限差分方程组的建立,井轴上场点奇异性与内边界处理等几个关键问题,提高了计算精度.在横向各向同性地层对称主轴与井轴平行的情况下,与实轴积分法所得结果进行了对比,验证了本文方法的正确性.计算了不同倾角情况下xx和yy两分量的弯曲波,并用频域加权相似法提取了弯曲波频散曲线,结果显示了横向各向同性介质中不同方向偏振的弯曲波传播的分裂现象,其频散曲线在低频段分裂,随频率增大而逐渐重合.弯曲波低频截止频率处速度与理论公式得到的横波速度基本符合.     

套管井应力集中下正交偶极声测井的声弹理论、数值分析与实例

油田受异常地应力作用而损坏的油井绝大部分都是套管井,针对这一实际情况,通过求解套管井应力集中情况下各层介质中应力场的分布,应用摄动积分方法,将裸眼井声弹性方程发展为适用于套管井井孔体系的情况.数值考察了井外无限远单轴水平压应力作用下,声弹方程中各层介质敏感系数和速度-应力系数随频率的变化,及二相互垂直偏振的偶极弯曲波相速度的频散曲线.结果表明,套管井声弹方程中的速度-应力系数主要由地层参数决定,内层介质的贡献较小,其中与三阶弹性模量相对应的项又起决定作用;井孔周围应力集中将引起平行和垂直单轴压应力方向偏振的二弯曲波发生速度分裂,同时随频率的增大两条频散曲线产生交叉,说明这种交叉仍然可用作判断套管井情况下地层各向异性是本征的还是由异常应力诱导的判据,用正交偶极声波测井技术在套管井中无损检测异常地应力仍然可行;但低速地层交叉点将明显向高频区移动,这又是与裸眼井情况重要的不同之处.本文还给出现场实测数据的处理实例,理论结果得到验证.     

VTI介质SH-SH波地震反演方法研究

近年来,随着横波可控震源技术的发展,国内外已经实现了纯横波地震勘探.相较于地震纵波,地震横波对横波阻抗、横波速度尤其是各向异性参数变化更为敏感,因此地震横波能够用来更好地估算上述地层参数.VTI(具有垂直对称轴的横向各向同性)是地层介质中广泛存在的一种各向异性形式,对振幅随偏移距变化(AVO)影响显著.本文提出了一种改进的VTI介质SH-SH波反射系数近似公式,新公式具有较高的精度且仅包含两项待求参数：横波阻抗和水平横波(SH波)速度.基于新方程建立了VTI介质SH-SH波反演方法,该方法相比VTI介质的PP波反演不确定性明显下降,同时降低了常规PP波各向异性反演对大角度数据的要求.为了获得独立的横波各向异性参数,进一步地提出了一种基于岩石物理关系的横波各向异性参数估算方法.合成地震数据测试和柴达木盆地九分量地震勘探实际地震数据应用结果表明,新方法能够准确地预测地层的横波阻抗、水平横波(SH波)速度、各向异性参数,为各向异性地层的岩性解释和油气储层预测提供了可靠的解决方案,从而深化了横波地震勘探的应用潜力.     

横向各向同性地层中随钻声波测井模式波分析

针对横向各向同性地层随钻声波测井模型,通过模式分析的方法,考察了快速地层和慢速地层井孔内随钻单极子、偶极子和四极子声源激发的斯通利波、弯曲波和螺旋波的相速度频散和激发强度特征,计算了这些模式波对于地层弹性常数的灵敏度,并与电缆测井中的情况进行了比较.结果表明:随钻斯通利波在低频时对地层弹性常数中c₆₆的灵敏度较电缆测井中有了很大提高,可用于反演地层水平向横波速度;随钻偶极子最低阶弯曲波在低频时不能用于直接获取地层横波信息,但在慢速地层中频率较高(例如6 kHz)时却可以间接得到地层垂直向横波速度;随钻四极子螺旋波的特征与电缆测井中的类似,可用于获取地层垂直向横波速度.     

井孔压力与地层应力集中诱导的偶极弯曲波的分裂

利用在大形变基础上叠加小扰动的非线性波动理论——声弹性理论，数值考察了井外均匀各向同性固体介质同时受井孔静压力和井外单轴水平压力作用时，Stoneley波和两种偏振的偶极弯曲波的相速度频散曲线，并研究了计算中用到的灵敏系数随频率的变化. 数值结果表明，模式波相速度频散曲线的变化主要受灵敏系数中对应的3个三阶弹性常数的项控制；由井孔周围应力集中引起的沿单轴压力方向和垂直于该方向偏振的弯曲波发生分裂，同时随频率的增大，两条频散曲线产生交叉；井孔静压力的作用使两条频散曲线交叉点的位置发生变化，交叉点相速度值随井孔压力的增加而增大，但对应的频率值不变.     

水平单轴应力与横波反射系数

地球物理勘探中,地下应力的估测是非常重要的.本文利用应力诱导的弱各向异性的各向异性参数,借助已有的具有水平对称轴横向各向同性介质顾及高阶小量的横波反射系数近似公式,针对下介质为各向同性受单轴应力作用情况随应力的增加求取应力精度降低的问题,从横波反射系数近似公式的推导出发,联合SH波与SV波反射系数给出新的改进的单轴应力与横波反射系数的直观表达式;计算并对比了受单轴应力作用的各向同性介质的横波近似和严格的反射系数,并利用理论得到的横波反射系数对应力进行了估值.结果表明,利用改进后横波反射系数与应力的直观表达式,当应力增加时,应力估值精度有了明显的改善,这为利用横波反射数据直观估测应力提供理论依据.     

利用伪谱法模拟横向各向同性介质中的波

介质的弹性常数为三维四阶张量的分量,共有８１个,由于应力张量和应变张量的对称性及能量密度是应变的二次函数,一般各向异常性介质的独立弹性常数可减为２１个,如果介质具有较高的对称性,独立弹性常数的数目会更少。 对于地壳和上地幔,具有５个独立弹性常数的横向各向同性介质是一个非常好的近似,本研究中横向各向同性介质的对称轴方向可以是任意的（即对称轴可以不平等于铅直方向）,在此情况下,需要进行坐标变换,如果已知介质在某一坐标系（其坐标轴平行或垂直于介质的对称轴）中的弹性常数,我们能够容易地利用变换公式得到变换后新坐标系中的弹性常数。 本文提出了一种方案,利用伪谱法既能模拟横向各向同性介质中的平面波,也能模拟点源激发的波场。在勘探地球物理和地震学中,模拟横向各向同性介拮中传播的平面波及区域源产生的波是最重要的研究课题之一。然而在一般各向异性介质中,很难或不可能确定弹性波的相速度和偏振方向,但在横向各向同性介质中,则可以通过坐标变换来实现,这里我们所提出的方法可以用于横向各向同性介质中弹性波的模拟。     

VTI介质纯P波混合法正演模拟及稳定性分析

各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性.     

VTI介质P波非双曲时差分析

具有垂直对称轴的横向各向同性介质模型(VTI)是目前各向异性理论研究和多波多分量地震资料叠前成像处理中最常用的一种各向异性模型.VTI介质中反射 P波时距曲线一般不再是双曲线.基于不同的相速度近似公式会得到不同的时距关系式.文中对几种典型的非双曲时距曲线与射线追踪得到的准确时距曲线在不同各向异性强度下进行了对比，结果表明Muir等和Stovas等提出的非双曲时距公式由于过高地考虑了横波垂直速度的影响与精确的时距曲线有很大偏差；Tsvankin等提出的弱各向异性非双曲时距公式在ε-δ<0时误差增大；Alkhalifah等提出的非双曲时距公式在大炮检距任意各向异性强度下都具有较高的精度，适于在实际资料处理中应用.     

水平应力作用下套管井不同胶结状况导波的声弹效应

油田的套管完井常伴有异常地应力的作用,对于套管井,井孔-套管-水泥环-地层组合的系统使井周的应力集中更为复杂.为认识地应力作用下不同胶结状态的套管井导波声弹效应,本文基于非线性声学理论,实现了受应力作用不同胶结状态套管井井周应力和径向横波速度计算与分析;采用摄动求解推导了套管井导波的声弹性方程,模拟并分析由应力引起不同胶结状态套管井斯通利波、伪瑞利波、弯曲波相速度的变化.研究结果表明:不同胶结状态会改变井周的应力分布,当界面出现滑移时井周地层中的应力集中有所增强,从而影响导波的声弹效应.本文的研究结果为地应力作用的套管井的声场特性提供认识和理论依据.     

Effect of stress on reflection and refraction of plane wave at the interface between fluid and stressed rock

Reflection and refraction of plane wave through an interface between fluid and rock with elastic deformations on the basis of the acoustoelastic theory are considered. The effects of stress on the anisotropy and energy reflection and transmission coefficients are investigated. The incident wave plane can coincide with or deviate from planes of material symmetry. Elastic deformations are assumed to be locally homogeneous and to satisfy static boundary conditions. Numerical computations are carried out and comparisons are made with the results predicted in the presence and absence of initial stress. The changes in phase velocity, group velocity, and energy reflection and refraction coefficients due to the presence of stress are displayed graphically and discussed. The results show that the effect of stress depends on its magnitude, direction and form (uniaxial and biaxial).     

孔、裂隙并存地层中的声波测井理论及多极子声场特征

在油、气储层的勘探和开发中观察到的一个现象是储层岩石中普遍存在孔隙和裂隙.随着近年来孔、裂隙介质弹性波动理论的进展,我们可以将此理论应用于测井技术,以此来指导从声波测井中测量孔、裂隙地层的声学参数.本文计算了孔、裂隙地层里充流体井眼中的多极子声场,分析了声场随裂隙介质的两个主要参数(即裂隙密度和裂隙纵横比)的变化特征.井孔声场的数值计算表明裂隙密度可以大幅度地降低井中声波纵、横波的波速和振幅.随着裂隙密度的增加,在测井频段内也可以看到纵、横波速的频散现象(这种频散在孔隙地层中一般是观察不到的).本文还研究了多极子模式波 (即单极的Stoneley波、伪瑞利波以及偶极的弯曲波)随裂隙参数的变化特征.结果表明,这些模式波的振幅激发和速度频散都受裂隙密度的影响.裂隙密度越高影响越大.此外,裂隙还对模式波的传播造成较大的衰减.相对裂隙密度而言,裂隙纵横比是一个频率控制参数,它控制裂隙对声场影响的频率区间.本文的分析结果对裂缝、孔隙型地层的声波测井具有指导意义.     

含孔隙、裂隙地层随钻多极子声波测井理论

储层岩石中普遍存在孔隙与裂隙,对钻井中的测井声波产生重要影响.基于孔、裂隙介质弹性波理论,导出了随钻声波测井的井孔声场表达式.据此考察了地层裂隙密度与含气饱和度的变化时井孔内随钻多极子模式波（斯通利波、弯曲波和螺旋波）的速度、衰减与灵敏度以及地层纵、横波的响应特征.裂隙密度与含气饱和度对模式波的速度频散与衰减都有影响,且两参数的值越大,影响越大.具体来说,速度对裂隙密度更敏感,而衰减对含气饱和度更敏感.具有"艾里相"特征的随钻偶极和四极子波在地层含气时产生强烈衰减,可以作为判断地层含气的一个明显指示.理论模拟与实际测井数据分析结果符合较好.     

裸眼井中弹性波传播的非对称模式的数值研究

本文从数值上研究了裸眼井中弹性波传播的非对称模式,给出了合成微地震图和导波（弯曲波）的频散曲线,发现在“硬”地层和“软”地层的井中,导波都是高度频散的,其最大相速度等于地层的横波速度,其截止频率低于对称模式的伪瑞利波的截止频率;在低频（2-3kHz）和长源距（3-4m）的条件下,由非对称的声源（如声偶极子）所产生的微地震图中,初至信号是以横波速度传播的,而以纵波速度传播的信号被抑制。本文的结果对研制横波速度测井仪是有意义的。     

横向各向同性介质和斜方介质各向异性参数的约束条件

Thomsen提出的横向各向同性(TI)介质各向异性参数(α0、β0、ε、δ、γ)是各向异性理论研究和实际资料处理中的常用参数,Thomsen参数的取值必须符合物理学定律和实际地学情况,随意的取值可能导致无意义乃至错误的结果.本文根据热力学定律和弹性常数的物理意义,结合大量的实测数据,提出常见TI介质的Thomsen参...     

套管井偶极弯曲波频散向高频偏移的特性

对套管井偶极弯曲模式波的频散特性进行了系统的数值考察、实例对比和分析.发现套管井弯曲波频散曲线随地层横波速度的降低,特别是地层横波速度小于2000 m·s^-1以下,会迅速移向高频区,偶极弯曲波基础模式主频散区（或截止频率）可出现在13 kHz以上,以致超出了现行低频偶极子声波测井仪的激发与接收频带,这是一过去没有被研究者注意到的现象,并进一步被现场实例所证实.研究表明控制套管井弯曲波频散曲线主频散区位置的主要是钢套管的厚度和地层横波速度.对地层横波速度大于井孔流体声速的快速地层,在钢套管壁厚一定（8 mm）的情况下,频散曲线主频散区可移至11 kHz以上,可能出现的最大可能频域位置是同一井孔内径,井外全钢时的频散曲线上等于、小于地层横波速度那一段,这对各种地层和套管参数都是适用的.对地层横波速度小于等于井孔流体声速（1500 m·s^-1）的慢速地层,弯曲波频散曲线随地层横波速度的降低移向高频区的特点更为明显,可能移至16 kHz以上;而套管厚度的影响,也比快速地层大的多,对地层横波速度小于1380 m·s^-1的慢速地层,无论用多高的频率激发,都不能在现行使用的各类套管井（壁厚6~12 mm）中用偶极声波测井仪测到弯曲模式波.     

含孔隙、裂隙致密介质中多极子声波的传播特征

非常规油气藏（如致密性地层及蕴藏油气的页岩地层）的重要特征是低孔、低渗,但裂隙或裂缝比较发育.为满足非常规勘探的需求,本文将孔、裂隙介质弹性波传播理论应用于多极子声波测井的井孔声场模拟,重点研究了致密介质中裂隙发育时多极子声波的传播机理以及衰减特征.井孔声场的数值计算结果表明裂隙的存在明显改变了弹性波和井孔模式波的频散、衰减和激发强度,尤其是井壁临界折射纵波的激发谱的峰值随着频率的增加逐渐降低,这与应用经典的Biot理论下的计算结果相反,且裂隙的存在也使得饱含水和饱含气时临界折射纵波激发强度的差异变大.井孔模式波的衰减与地层横波衰减和井壁流体交换有关,井壁开孔边界下致密地层裂隙发育还使得井孔斯通利波和艾里相附近的弯曲波对孔隙流体的敏感性增强,在井壁闭孔边界条件下引起井孔模式波衰减的主要因素是裂隙引起的地层横波衰减造成的,且在截止频率附近弯曲波的衰减与地层的横波衰减一致.数值计算结果为解释非常规油气地层的声学响应特征提供了参考.