方位界面及其波场数值模拟

研究各向异性介质是勘探地震学中的一个重要内容。方位界面是各向异性介质中的一种特殊现象。本文首先解释了方位界面的含义，然后以裂隙各向异性介质为例，讨论了方位界面存在的三种情况，并给出每种情况下的数值模拟动方程；最后，用有限方法作了波场数值模拟。模拟结果表明方位界面也是波阻抗界面。本文讨论问题的方法可以扩展到ＴＩＭ、ＰＴＬ等其它六方各向异性介质。     

各向异性介质中弹性波传播特性的伪谱法模拟研究

本文首先研究了各向异性介质中弹性波动方程伪谱法求解的稳定性条件，然后应用伪谱法模拟研究了各向异性介质中三种常见震源（垂直力源、压力源和剪切源）的辐射特征和三种各向异性模型的波传播快照特征；最后模拟研究了力源作用于裂隙诱导各向异性介质时ＶＳＰ记录的特征。所有这些工作表明用伪谱法研究各向异性介质中弹性波传播是有效的，与此同时，我们也获得了一些各向异性介质中波传播的有益认识。     

正交介质中地震波传播数值模拟与分析

建立了正交各向异性介质中的波动方程,应用高阶有限差分进行数值模拟,得到三分量地震波场快照,并在此基础上,分析了正交各向异性介质中地震波场特征,结果表明：在正交各向异性介质中,地震波场十分复杂,三个分量相互藕合,两个水平分量波的偏振方向与裂隙走向既不平行,也不垂直。     

双相各向异性介质地震波传播特征

实际裂隙性油气储层表现出双相各向异性特征。本文基于岩石物理理论求取裂隙介质各向异性等效弹性参数,结合各向异性BISQ理论,建立了裂隙诱导双相各向异性介质理论模型。从双相各向异性介质的弹性波方程出发,推导出平面波的Christoffel方程,求解该方程获得了三维双相完全各向异性介质中qP1、qP2、qSV和qSH波的波数方程,据此波数方程计算四类波的相速度和衰减因子。最后针对碎屑岩水平裂隙诱导双相各向异性介质,分析了泥质含量、裂隙密度和形状、孔隙度、含水饱和度、流体黏滞系数、渗透率等参数对弹性波沿不同方向传播时相速度和能量衰减特征的影响,可为储层的孔隙、裂隙结构及含油气性等分析提供理论依据。     

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

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

利用纵波地震资料确定方位各向异性的主方向

在方位各向异性介质中，方位各向异性的主方向就是似Ｐ波和似Ｓ波以纯Ｐ波和纯Ｓ波方式传播的方向。如果方位各向异性是由作用于一个方位各向同性背景中的定和垂直裂缝引起的，则这些主方向中的其中两个对应平行和垂直于裂缝的方向。横波穿过一个方位各向异性介质时，对与主方向有关的传播方向很敏感。因此，可以用多分量的一次Ｓ波或转换Ｓ波资料来获得主方向。除了采集费用高以外，地震工业界目前还没有完全开发出多分量资料处理和     

犹他州Bluebell—Altamont油田裂隙气藏P波和S波方位各向异性响应

本文利用Ｐ波和Ｓ波反射数据描述第三系ＵｐｐｅｒＧｒｅｅｎＲｉｖｅｒ组天然裂敢世故的方位各向异性并阐述两种数据体的优点和能力。目的是评估Ｐ波地震资料检测气充填裂隙的可行性，估算裂隙密度和方位，并与Ｓ波的结果作比较。垂直裂隙的Ｐ波响应必须根据不同炮－检方位（传播路径）与隙走向的角度 进行分析。两条多分量反射资料的测线埋置为近似平行和垂直ＵｐｐｅｒＧｒｅｅｎＲｉｖｅｒ组裂隙的主方向的该方向由 岩 分析和     

横向各向同性介质中弹性波模拟及波场特征研究

某些结晶或含有裂隙的介质表现出来与弹性波传播有关的各向异性，这一性质可以用弹性波动理论建立的质点运动方程和应力应变本构方法（即广义虎克定律）来描述，本文利用傅里叶数值方法（即伪谱法）求解该方程来模拟各向异性介质中弹性波的传播。     

弹性介质三参数等相位面法层析成像

本文给出了一种利用散射Ｐ波进行弹性介质密度与拉梅系数二维成像的等相位面方法。该方法首先从Ｐ波方程出发，在介质小扰动和远场条件下导出了Ｐ－Ｐ散射波的广义Ｒａｄｏｎ变换公式；而后依据Ｐ－Ｐ波形上每一点与其相应等相位面所穿过的介质参数扰动建立代数方程。     

含裂隙介质中地震波运动学问题的正演模拟

碳酸岩油气藏与裂隙分布有关。在含有定向分布的裂隙介质中,地震波的传播呈现各向异性的特征。本文通过对定向分布裂隙介质中的弹性参数的研究,发现裂隙介质具有横向各向同性介质的特性,而且裂隙内含有液体,对波速的各向异性影响更甚。文中介绍了各向异性介质中地震波的射线解,并利用运动学射线追踪方程计算了裂隙介质中波的射线轨迹及走时曲线,求出似 P波、似 SH 波和似 SV 波的波速。由于介质的各向异性,似 SH 波和似 SV 波的走时曲线表现不一致,这说明各向异性介质中 S 波存在分裂现象。此种 S 波分裂现象可以用于研究介质中的裂隙分布和指示裂隙油气藏。     

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

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

横向各向同性介质中的AVO

介质各向异性是影响振幅随炮检距变化的重要因素之一。本文将Ａｋｉ和Ｒｒｃｈａｒｄａ的方法进行推广，推导出两个弹性横向各向同性介质密接条件下平面波反射和透射系数，并给出弱弹性ＴＩ假设下的近似式。近似式的物理意义明确，而且精度明显于前人的结果。     

裂缝介质转换波方位各向异性研究

该文介绍了在弱各向异性条件下,基于Thomsen 裂缝理论---HTI 介质各向异性参数的转换波近似反射系数公式,通过AVO 正演,比较分析各向异性系数对于P-SV 波反射系数的影响及对P-P 波反射系数的影响。结果表明,应用联合P-P 波和P-SV 波的数据进行联合反演,能够更好地降低数据结果的不唯一性,从而得到更精确的裂缝参数。     

裂缝性储层AVO模型研究

预测储层中的有利裂缝发育带及裂缝发育的主方位在油气勘探开发中具有重要的意义。当储层中存在裂缝时，流体与致密岩石的巨大物理性质差异形成一定规模的物理界面，引起储层地震波反射特征的改变，成为地震剖面上反映裂缝存在的重要信息。AVO分析方法所表现的方向性，可从三维道集的地震资料中分析出不同方向上的地震属性参数变化，这些参数与地层的各向异性具有一定的对应关系。当储层存在一定方位的裂缝时，其弹性参数会表现出方向性差异，这种差异能被AVO方法有效地检测。通过对不同方位角裂缝储层AVO模型研究得出：当地震波传播方向与裂缝走向的夹角逐渐增大时，反射系数随入射角的增大而减小；含水平裂缝地层的曲线斜率逐渐减小，地层的AVO截距逐渐减小，斜率逐渐增大；垂直裂缝地层的曲线斜率逐渐增大，地层的AVO截距逐渐增大，斜率逐渐减小。     

裂缝介质中多方位AVO特征分析技术及应用举例

有利裂缝发育带及裂缝发育的主方位在储层预测中具有重要的意义。当储层中存在裂缝时，流体与致密岩石的巨大物理性质差异形成一定规模的物理界面，引起储层地震波反射特征的改变，成为地震剖面上反映裂缝存在的重要信息。地层中分布的裂缝，特别是在构造应力作用下形成的具有一定方位性分布的裂缝，其特征可以采用多方位AVO 方法进行检测。理论上讲，裂缝对地震属性参数的影响在垂直裂缝走向的方向上表现为最强，在水平裂缝走向的方向上表现为最弱，对AVO 属性参数的影响也是如此。因此，对于发育裂缝的储层采用多方位AVO方法进行矢量差异分析可以确定储层中裂缝发育的强度和密度。从理论公式入手讨论了不同入射角、方位角以及不同裂缝倾角情况下裂缝性储层的方位AVO响应，研究了裂缝性储层多方位AVO的变化特征与裂缝走向和倾向的关系，尤其是具有一定倾斜角度的裂缝对地震波方位AVO的响应特征，并用这一方法分析了新疆某油区的有利裂缝发育带，取得了较为理想的效果。     

裂缝介质中地震波方位AVO特征分析

预测储集层中的有利裂缝发育带及裂缝发育的主方位在油气勘探开发中具有重要的意义,裂缝的存在将导致介质属性的变化,不同的裂缝发育模式具有不同的表现特征,地震反射波振幅的方位变化就是一个重要的特征。对不同裂缝各向异性地质模型的地震波方位AVO特性进行了分析,建立地震波方位AVO的变化特征与裂缝走向和倾向的关系,尤其是具有一定倾斜角度的裂缝对地震波方位AVO的响应特征,为利用振幅变化信息预测裂缝提供了依据。将其应用于千米桥古潜山有利裂缝发育带的分析,取得了一些效果,促进了该方法在油气勘探中的应用。     

正交各向异性介质反射系数精确解

中国陆相沉积岩大多具有薄互层特征,当受到构造运动的影响时,会发育垂向或近似垂向的高角度裂缝,呈现出典型的正交各向异性,可看作是VTI和HTI各向异性的叠置。研究正交各向异性介质的AVO响应特征对裂缝型储层的精细刻画有重要的意义。针对VTI背景介质中发育的一组直立裂缝诱导的正交各向异性,采用Tsvankin提出的各向异性参数构建刚度系数矩阵,再根据Christoffel方程和边界条件推导了精确反射系数和透射系数的计算方法。理论模型试算表明：①若P波从低阻抗各向同性介质入射至高阻抗正交各向异性介质,背景介质VTI各向异性强度的增大会导致PP波反射系数增大、PS1波和PS2波的反射系数减小;而随着裂缝弱度的增强会导致PP波反射系数的减小、PS1波和PS2波的反射系数的增大。②若P波从高阻抗各向同性介质入射到低阻抗正交各向异性介质,PP波、PS1波和PS2波的反射系数变化规律与①相反。③PP波反射系数的方位各向异性远小于PS1波和PS2波。与Rüger反射系数近似公式计算结果的对比,验证了该反射系数计算方法的正确性。     

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.     

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%.