横向各向同性介质拟声波一阶速度—应力方程

基于弹性波的基本理论,首先推导了VTI介质拟声波一阶速度—应力方程,通过坐标变换,还得到了TTI介质拟声波一阶速度—应力方程,并从能量守恒的角度证明了上述方程比基于频散关系推导出的拟声波方程更加稳定。针对各向异性参数ε<δ的模型,对新推导的方程引入横波项,消除了残留的横波产生的数值模拟不稳定现象。通过对均匀二维模型和逆冲模型的正演模拟表明:新推导的方程不仅能很好地维持qP波的运动学特征,而且能够很好地适应任意旋转角度的TTI模型,并且在旋转角变化的TTI介质中数值稳定性要好于基于频散关系推导出的拟声波方程。     

利用不同角度域P波资料反演纵、横波速度

基于反射P波垂直入射假设的常规地震道反演方法不能得到可靠的速度或其他岩性信息，而利用叠前地震资料按部分角度道集叠加后，可以获得多个具有AVO特性的地震剖面，信噪比也得到显著改善，因此角度部分叠加道集可以作为弹性阻抗反演或纵、横波速度同步反演的基础道集记录。对每个角度范围的部分叠加剖面进行反演处理都会得到该角度范围内的阻抗或纵、横波速度值，但是不同角度范围内反演结果之间的关系以及如何求取零入射角的纵、横波速度值仍然是难以解决的问题。为此本文提出利用地震波传播速度在各向异性介质中随入射角变化的理论，由各个角度范围的纵、横波速度反演结果求取地震波零入射角的纵、横波速度的方法，实现了利用叠前P波资料同步反演地震纵、横波速度，进而可以求得纵、横波速度比和泊松比剖面，为解释人员进行岩性及含油气性解释提供可靠的依据。     

爆破源地震反射波三分量记录理论地震图研究

本文用各向异性反射率法计算爆破源反射地震的理论地震图。对比了各向同性介质的反射记录,讨论了水平叠层各向异性和裂隙各向异性对体波三分量记录特征的影响。结果表明,裂各向异性导致转换横波分裂,并使横向记录分量出现波动。     

Tau-p变换若干问题的讨论及其在反演声纳浮标资料中的应用

在时空域进行τ-p变换,由于接收空间范围有限,以及反射波时距曲线在远、近道被截断,必将引起假频和端点效应的干扰。采用相似系数滤波方法可以压制这种假频和端点效应的干扰。为了增加相似系数的平稳性在作相似滤波之前,需对相似系数进行低通滤波,以消除随机干扰。在τ-p变换中,参数选择尤为重要。参数τ可根据有效波出现的范围、不出现假频和计算速度等因素综合考虑。参数p(和△p)可根据最大和最小水平视速度、假频和计算速度等因素考虑。文中还讨论了有偏移距和倾斜界面模型在τ-p变换时出现的问题及τp变换在道内插和压制面波干扰中的应用。在τ-p域内应用τ-SUM法反演声纳浮标资料,可获得很好的效果。     

具有垂直和水平对称轴的横向各向同性模型的P波反射系数

对振幅随偏移距变化的分析来说,研究各向异性介质中的P波反射系数是重要的。直接对反射系数做数值估计时,数值解不能从解析角度反映各向异性对AVO特征的影响。 为解决这个问题,本文为计算具有垂直对称轴的横向各向同性介质中水平边界的P波反射系数提供一种改进了的近似方法,其解中有一个与以前文献中方程相同的描述反射系数低次项角度变化的AVO梯度项,但对大入射角来说,这个解更为确切。这样,精细的近似可以扩展到具有水平对称轴的横向各向同性(HTI)介质。这种情况通常是由一组垂直裂缝所引起的。对比HTI介质的两个垂直对称平面中入射P波的近似反射系数可以看出:AVO梯度的方位变化是横波分裂参数γ和各向异性参数的函数。其中,各向异性参数描述在具有对称轴的垂直平面中近似垂直传播的P波各向异性。     

基于Lebedev网格的TTI介质二维三分量正演模拟

传统交错网格有限差分法是研究地震波传播规律的一种较为常用的数值模拟方法,但是在交错网格中每个变量的不同分量都是交错定义的,对于没有定义的点需要进行变量插值,从而降低了模拟精度。为此,在前人的研究基础上,推导了TTI介质二维三分量的应力—速度方程,采用Lebedev网格对其进行了高精度的差分离散处理,避免了传统交错网格在处理各向异性介质时波场插值引起的误差,提高了模拟精度,并将多轴完全匹配层吸收边界（M-PML）引入了Lebedev网格。分别对单层TTI介质和含透镜体的复杂TTI介质模型进行了正演模拟,结果表明：①由于Lebedev网格对各向异性介质的弹性波方程做离散时无需进行波场插值,与传统交错网格有限差分方法相比,模拟精度更高;对TTI介质进行模拟时可以清晰地观察到纵波、快横波和慢横波,并且快、慢横波的偏振方向相反,在单炮记录中观测到的三种波的速度特征也符合波场传播规律。②引入多轴完全匹配层（M-PML）吸收边界条件后,在不影响模拟效果的情况下边界反射现象被有效地压制。     

各向异性介质中反射率法波场模拟及体波辐射图案

本文从一般各向异性的弹性动力学方程出发，采用反射率法来计算任意各向异性的黏弹性层状介质中的全波场合成地震记录。系统矩阵的特征解是在三维空间计算的，考虑了偏离径向平面的波传播。时间假频由复频率技术来压制，空间假频则通过使用Ｆｉｌｏｎ法估计两个波数的积分来避免。文中以裂隙介质和叠层－裂隙复合介质为例，给出了各向异性介质中合成地震记录计算的实用方法，并分析了横波分裂现象及振幅特征。     

各向异性介质弹性波高斯束偏移方法

弹性波高斯束偏移是一种兼顾成像精度与计算效率的多分量地震成像方法,但当前的研究多集中在各向同性介质声波、弹性波成像和各向异性介质声波成像,针对各向异性介质弹性波偏移的文献相对较少。为此,在前人的研究基础上,发展了各向异性介质弹性波射线追踪方程,以二维各向异性弹性波Kirchhoff-Helmholtz积分为基础,利用弹性动力学表征的格林张量推导了各向异性介质中弹性波波场正、反向延拓公式,提出了一种针对各向异性介质多分量地震数据的成像方法。通过在成像公式中引入权函数,很好地压制各向异性介质成像中不同波型引起的串扰。通过引入符号函数,解决了转换波成像过程中的极性反转问题。TTI介质洼陷模型和TTI介质多层构造模型试算结果表明：与各向同性介质弹性波高斯束偏移方法相比,各向异性介质弹性波高斯束偏移成像结果的信噪比较高、同相轴连续性较好、构造位置更准确;使用多分量地震记录,利用纵、横波波场信息成像,有效提高了成像分辨率;与传统的基于标量波场的成像方法相比,所提方法的适应性更好。     

复杂TTI介质稳定的纯qP波波场模拟方法

作为各向异性逆时偏移技术的基础,复杂各向异性介质情况下精确而稳定的波场模拟至关重要。本文给出一种针对复杂倾斜横向各向同性（tilted transversely isotropic,TTI）介质的稳定的纯qP波波场模拟方法。该方法基于Xu提出的伪微分算子分解思路,通过分析新算子的频散特性,引入旋转坐标系下的自共轭算子以保证稳定性,导出新的TTI介质一阶纯qP波控制方程;在Lebedev交错网格框架下推导了新方程的高阶有限差分形式,给出2次计算波场梯度的数值算法实施策略,以进一步保证波场模拟精度。均匀各向异性介质及复杂TTI介质模型的数值试算结果表明：新控制方程不受伪横波影响;相对于有限横波qP波波场模拟等算法,新算法可得到更稳定的各向异性纯qP波波场,即能更适应各向异性对称轴参数的空间变化,可应用于高精度纯qP波逆时偏移,改善对各向异性介质的成像质量。     

Tau-P变换及其初步应用

层状介质的共炮点或共中心点道集记录经τ-p变换后,直达波变成了一个能量团,反射波变成了椭圆,其交会处既是折射波的能量团,也是下一层在无穷处的反射。这样τ-p变换后就把反射波、折射波有机地联系起来,原来作为干扰处理的噪音(折射波等)就变为有用的信息。由于p是水平视慢度,因此纵,横波在τ-p域上有明显的差异。文中按照精确的时距曲线公式建立了具有多个反射层的共炮点理论记录,包括直达波、浅层折射波,各层的反射波及折射波。文中给出了一些消除线性干扰的实例,并就消除变换中的混叠现象作了尝试。     

盖层各向异性对时移地震AVO的影响

油藏盖层各向异性在油藏开发前后对地震AVO曲线有相同程度的影响,而这种影响能否在2次地震AVO曲线相减后变得很小,可以被忽略？利用某油田数据建立了储层岩石物理模型,根据Gassmann方程获得了该油田储层的纵、横波速度与含油饱和度的关系曲线.设计了双层油藏模型,通过求解横向各向同性（VTI）介质界面的透射系数与反射系数的线性方程组,计算了油藏盖层存在各向异性时的P-P波和P-S转换波反射系数,分析了不同Thomsen参数对反射系数及其油藏开发前后反射系数差异的影响,并进行了定量分析.分析结果表明,小入射角情况下进行P-P波时移地震AVO分析,可以忽略盖层各向异性的影响;大入射角情况下进行P-P波时移地震AVO和多波时移地震AVO分析,盖层各向异性不能忽略.     

井间地震资料层位标定方法研究

制作高精度的深度域合成地震记录是实现井间地震资料层位标定的关键。介绍了采用深度域不对称子波算子求取合成地震记录的方法，即对已有的测井曲线，利用有效地震地质模型分析方法进行层位划分并建立深时关系；根据由波动方程导出的反射系数计算公式，求取频率域的反射系数；把井旁道记录转换到时间域，利用相关分析方法提取实际资料的地震子波，进而得到频率域的子波谱；利用频率域褶积模型计算合成记录，经傅里叶反变换、时深转换和重采样后得到深度域合成地震记录。实际应用表明，采用深度域方法制作的合成记录与井旁地震道吻合程度较高，可以用于井间地震资料的层位标定。     

拟合AVO属性反演

本文从Shuey的佐普里兹方程近似式出发，通过拟合随入射角变化的反射系数，获得较精确的零炮检距反射系数，提高了地震资料的信噪比；而且这种拟合方法不会从根本上改变原有地震资料CMP道集上反射振幅随炮检距变化的关系，进而反演出较准确的AVO属性剖面。通过理论模型和实际地震数据试算结果，表明该方法与以往采用的AVO反演的主要区别在于不仅能够反映炮检距上振幅能量随炮检距变化的关系，而且对入射角没有限制。     

Geochemical characteristics and diagenetic fluids of dolomite reservoirs in the Huanglong Formation, Eastern Sichuan Basin, China

Based on a comprehensive study of texture,diagenetic behavior and evolution of dolomite in the Huanglong Formation,trace （e.g.,Fe,Mn and Sr） and rare earth element （REE） geochemistry,andisotopic characteristics （e.g.,C,O and Sr）,four types of diagenetic fluids are identified in the Huanglong Formation dolomite reservoirs of the Eastern Sichuan Basin,China：1）：marine-derived pore waters in the marine diagenetic environment,2） sabkha compaction brine conserved in the early shallowburied diagenetic environment,3） strongly-oxidizing low-temperature meteoric water in the seepagesubsurface flow diagenetic environment,and 4） strongly reducing deeply seated mixed hot brine in the middle and deep burial diagenetic environment.The fluids developed hereditarily from one environment to another,which resulted in its respective characteristics.Fluid characteristics play an important role in the development of dolomite reservoirs：1） dolomitization by marine-derived pore water in the quasisyngenetic stage did not form an effective reservoir; 2） early diagenetic burial dolomitization by the sabkha compaction brine formed the basis for reservoir development; 3） meteoric water karstification in the paleo-epidiagenetic stage expanded both the distribution and the size of the reservoirs,and improved the reservoir quality; 4） deep-burial dissolution and tectonic fracturing in the reburial diagenetic stage further improved reservoir porosity and permeability.     

地震波AVO与地层岩性分析

介质的各向异性，孔隙率、流体饱和度及有效压力等是影响地震波振幅随炮检距变化（AVO）的重要因素，本文在现有的地震波反射系数表达式基础上，通过合理的假设，推导出弱横向各向同性（TI）介质中的反射系数近似表达式，在不降低精度的同时，表达式的形式大大简化（振幅是炮检距或入射角的显函数），物理意义明确，并且研究了双相介质（水－气，油－气）孔隙率，流体饱和度和上覆岩层压力变化对反射波AVO的影响，由于岩石孔隙流体饱和度对介质密度，纵波和横波速度及反射波振幅影响较大，本文尝试把饱和度变化对反射系数的影响表示成关于速度和密度值的显函数（饱和度修正项），从而使长期以来被人们定性分析的饱和度对速度，密度及反射系数的影响定量化，对于不同压力、不同孔隙流体情形，得到的饱和度修正项均适用。     

Improving seismic interpretation: a high-contrast approximation to the reflection coefficient of a plane longitudinal wave

Linearized approximations of reflection and transmission coefficients set a foundation for amplitude versus offset（AVO） analysis and inversion in exploration geophysics.However,the weak properties contrast hypothesis of those linearized approximate equations leads to big errors when the two media across the interface vary dramatically.To extend the application of AVO analysis and inversion to high contrast between the properties of the two layers,we derive a novel nonlinearized high-contrast approximation of the PP-wave reflection coefficient,which establishes the direct relationship between PPwave reflection coefficient and P-wave velocities,S-wave velocities and densities across the interface.（A PP wave is a reflected compressional wave from an incident compressional wave（P-wave）.） This novel approximation is derived from the exact reflection coefficient equation with Taylor expansion for the incident angle.Model tests demonstrate that,compared with the reflection coefficients of the linearized approximations,the reflection coefficients of the novel nonlinearized approximate equation agree with those of the exact PP equation better for a high contrast interface with a moderate incident angle.Furthermore,we introduce a nonlinear direct inversion method utilizing the novel reflection coefficient equation as forward solver,to implement the direct inversion for the six parameters including P-wave velocities,S-wave velocities,and densities in the upper and lower layers across the interface.This nonlinear inversion algorithm is able to estimate the inverse of the nonlinear function in terms of model parameters directly rather than in a conventional optimization way.Three examples verified the feasibility and suitability of this novel approximation for a high contrast interface,and we still could estimate the six parameters across the interface reasonably when the parameters in both media across the interface vary about 50%.     

TTI介质一阶qP波方程交错网格伪谱法正演模拟

受TTI声介质强各向异性和对称轴极化倾角急剧变化的影响,传统有限差分法正演模拟易出现不稳定现象。本文基于TTI介质二阶耦合准纵(qP)波方程,推导了TTI介质一阶速度-应力方程,构建了交错网格中基于完全匹配层(PML)吸收边界条件的高阶有限差分格式;通过引入基于PML原理的一阶qP波方程伪谱法递推格式,实现了TTI介质qP波波场的稳定延拓。对二维TTI介质Wedge和BP模型的正演模拟试算结果表明:基于伪谱法递推格式实现的TTI介质正演模拟法,能较准确且稳定地模拟强各向异性及倾角变化剧烈区域的地震波场,验证了本文方法的有效性和较强适应性。     

各向异性含气砂岩模型正演及AVO响应特征分析

为认识和研究各向异性含气砂岩储层中弹性波的多波波场特征,以利于多波地震资料解释。在各向异性弹性波方程基础上,采用高阶交错网格有限差分技术模拟地震波在各向异性含气砂岩模型中的传播,得到该类储层的弹性波场。同时结合AVO分析技术分析了各向异性系数对多波波场特征的影响,比较了地震波在各向同性含气砂岩模型和各向异性含气砂岩模型中的波场响应异同。对比研究结果表明,地层的各向异性对波场特征和AVO响应有明显影响,并且这种影响在纵波和转换波上表现不一致;正演记录中的振幅特征与AVO异常响应是吻合的,说明将正演模拟与AVO分析相结合可为准确识别多波波场提供一种思路。     

三维TTI介质弹性波相、群速度的近似配方表征法

TTI介质弹性波相速度可以通过求解Christoffel方程得到,其精确式往往比较复杂,难以运用。为此,从Thomsen各向异性参数表征的TTI介质弹性波精确三维相速度表达式出发,利用近似配方法推导了TTI介质qP波和qSV波三维相速度的近似公式;将其代入Berryman的群速度表征中,得到了qP波和qSV波群速度的近似公式。理论分析和数值计算表明,基于近似配方法的相、群速度表达式在一定范围内与精确值吻合较好,能够适用于中强各向异性情况,且Thomsen参数ε与δ越接近,其精度越高;沿TTI介质对称轴方向和垂直对称轴方向的速度误差最小,而与对称轴呈45°方向误差相对较大;在相同各向异性强度下,qP波相、群速度误差一般要小于qSV波。