基于贝叶斯理论的VTI介质多波叠前联合反演

以VTI介质Ruger近似反射系数公式为基础,研究了多波叠前多参数联合反演方法。该方法以VTI介质Ruger近似公式为AVO正演方程,联合应用转换波和纵波数据,利用最小二乘原理构建目标函数来反演密度比、速度比、各向异性参数差等5个参数。在反演过程中引入贝叶斯理论,假定先验信息服从高斯分布,待求参数服从改进的Cauchy分布,并对待求参数去除相关性。对某地区地层模型正演的多波数据和含噪声数据分别采用本文方法进行了反演,与单独利用纵波数据反演的结果相比,联合反演的结果精度更高、抗噪能力更强、稳定性更好。模型数据测试结果验证了本文方法的可行性和有效性。     

基于傅里叶级数展开的纵波方位各向异性裂缝预测

利用地震反射纵波的方位各向异性信息预测地下裂缝的发育和分布情况一直是地球物理学家研究的热点课题,常规方法是使用两项Ruger线性近似公式反演得到裂缝相对密度和走向,其计算结果存在多解性且只能计算入射角较小的情况。对任意对称面各向异性介质的纵波方位各向异性弱反射系数近似公式进行整理,结合Schoenberg线性滑移理论,得到纵波方位各向异性反射系数近似公式的傅里叶级数展开公式,给出了基于傅里叶级数展开的纵波方位地震数据裂缝预测方法。将该方法应用于辽河油田雷家地区沙河街组白云岩储层的裂缝预测,并与常规各向异性反演结果及实际成像测井数据进行对比,验证了方法的有效性。研究结果表明,该方法计算简单且不受入射角度的限制,预测的裂缝相对密度与常规反演方法预测的结果精度近似,与成像测井资料吻合,故可作为传统方法的验证和补充。     

纵波速度频散属性反演方法研究

利用与频率有关的地震属性进行储层特征描述一直是地球物理学家研究的热点,由于含烃类储层会造成地震波的衰减和速度频散,因此,反映速度频散的属性在一定程度上可以指示含烃类储层。基于弥散介质界面反射系数随频率变化的假设,通过研究纵波法向入射时反射振幅随频率的变化规律,推导了包含表征纵波速度频散程度属性的反射系数近似公式,构建了反演方程,通过小波域分频和贝叶斯三参数反演分别得到分频地震数据和密度相对变化量,最终反演得到频散属性数据体。模型实验和实际资料的应用结果表明,纵波速度频散属性对含气储层具有较好的指示效果。     

黏弹性介质纵波反射系数近似

本文从黏弹性介质精确Zoeppritz方程出发,以黏弹性介质为背景,根据介质的分解理论,利用弱黏弹性近似和相似介质近似推导了黏弹性介质反射系数近似公式,该式分为实部和虚部,在实部中既存在弹性项,也存在黏弹性项,更精确地反映了实际地层介质中波的传播规律。特别地,当品质因子为无穷大、衰减角为零时,该式与以完全弹性介质为背景得到的反射系数近似公式(Aki-Richards近似公式)完全吻合,证明了公式的合理性。最后利用四类AVO模型,采用不同的黏弹性参数对黏弹性介质反射系数近似公式的有效性进行验证,结果表明,与Aki-Richards近似公式相比,以黏弹性介质为背景推导出的纵波反射系数近似公式的精度更高。     

基于傅里叶级数分析的各向异性参数估计及裂缝预测

利用Ruger反射系数近似公式进行地震方位各向异性椭圆拟合,研究裂缝密度与走向,在实际应用时因为待拟合点中存在奇异点,且椭圆拟合存在90°旋转问题,因而计算的方位角无法准确指示椭圆的长短轴方向。为此,将HTI模型纵波反射系数近似公式展开成傅里叶级数形式,利用最小二乘拟合提取傅里叶级数系数,从而提高各向异性参数估计和裂缝预测的准确性。针对不同方位角范围的数据,利用二倍角公式将纵波反射系数傅里叶级数展开式进一步展开,避免拟合提取傅里叶级数系数时因为同一个cos~nφ项对应两个各向异性反射系数而导致拟合结果出错。将该方法应用于HTI介质模型和实际地震资料,结果表明:对于HTI介质模型,展开后提取的拟合傅里叶级数系数与理论公式计算结果一致;对于实际地震资料,展开后提取傅里叶级数系数与工区裂缝发育及井资料指示的结果吻合,效果优于椭圆拟合结果和未展开时的傅里叶级数系数拟合结果。     

基于佐普里兹方程的高精度叠前反演方法

本文利用佐普里兹方程及其近似公式计算了四类具典型AVO特征地层模型的反射系数曲线,计算结果表明:在小角度(小于30°)入射时,利用Aki-Richards、Shuey近似公式计算的反射系数精度较高,与佐普里兹方程计算的反射系数非常接近;当入射角大于30°时,利用近似公式计算的反射系数与佐普里兹方程计算的反射系数的差别较大。为此,本文提出了直接利用佐普里兹方程进行叠前反演,并通过引入广义线性反演理论和贝叶斯方法详细推导了叠前反演算法。理论模型试算结果表明,基于佐普里兹方程的叠前反演能够获得精度较高的反演结果。     

横向各向同性介质纵波方位各向异性物理模型研究

 Ruger公式表明：在各向异性介质中,当各向异性参数、纵波阻抗、横波速度、介质密度已知的情况下,反射系数是方位角和入射角的二元函数;若采用固定炮检距观测方式,只改变测线与裂缝方向之间的夹角,则反射系数可化简为方位角的一元函数,根据Ruger公式便可得到不同炮检距下的反射系数值。本文基于Ruger公式,运用物理模型模拟方法,在水中以50,100,150mm的炮检距观测了HTI介质的地震响应,了解了反射波振幅、反射旅行时随裂缝方位角变化的规律。研究表明：在不同炮检距情况下,裂缝层顶、底界面的反射纵波振幅均出现方位AVO响应,即当方位角为0°时,纵波反射振幅最大,当方位角为90°时,纵波反射振幅最小;裂缝对纵波旅行时也产生影响,即反射时间随方位角增大而增大。物理模拟实验结果与理论计算结果相近,从而为裂缝识别和预测提供了借鉴。     

一种新的各向异性弹性阻抗近似公式

弹性阻抗是声波阻抗的拓展，可以更好地反映储层的性质。忽略各向异性影响的各向同性弹性阻抗公式不能正确地描述各向异性介质中弹性阻抗随入射角和方位角的变化。为此，在方位各向异性（HTI）介质模型的基础上，利用Connolly推导弹性阻抗的思想，根据各向异性反射系数近似公式，推导出了适合非零炮检距的各向异性弹性阻抗近似公式（简称新近似公式）。通过数值模拟，对纵波弹性阻抗随入射角和方位角的变化特征进行了描述，分析了各向异性参数对不同方位纵波弹性阻抗的影响，以及新近似公式的精度。分析结果表明，在弱各向异性介质且横波分裂参数γ不为0时，新近似公式比Martins近似公式的精度高。     

几种P-SV转换波反射系数近似公式的比较

多波多分量勘探利用了纵波和转换波的AVO属性来了解复杂的构造、岩性、裂隙裂缝的情况,以及确定地下流体的性质等。影响AVO属性的因素有很多,其反演的多解性很强,将纵波和转换波的AVO属性综合起来应用,就会大大减少多解性,还可以得到更多的岩石物理参数,并能提高一些参数的精度。AVO近似方程可以对单一的纵波、转换波,或纵、横波联合进行AVO岩性参数拟合反演。为此,对P—SV转换波反射系数近似公式进行了讨论和总结,并用理论模型对近似公式的特点和精度进行了分析。     

利用地震,测井资料联合反演储层物性参数

地震、测井资料联合反演技术，是利用地震，测井数据的互补特性，从井间地震数据中反演地层的物性参数。共分三步进行：（１）查清井间地层构造形态，应用Ｂｏｒｎ近似反射线性反演技术，从叠前或叠后地震记录中得到反映层速度界面的反射系数剖面，效果相当于有或叠后的偏移剖面，并用作非线性反演的井间几何形态约束信息；（２）求取井间地震波参数，应用最优拟合牛顿法学线性反演技术，分别得到地层反射系数，波阻抗和介质速度参数     

海底多波多分量AVO蕴含弹性参数信息的定量分析

海底多波地震勘探利用海底电缆采集海底各岩层的弹性波的多波多分量信息,在利用随偏移距变化进行弹性参数反演中,这些分量对反演所作的贡献是不同的,体现在各反射波的反射系统与弹性参数的信息量及其组合关在联程度是不同的。     

裂缝性储层AVO模型研究

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

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

方位弹性阻抗傅里叶级数裂缝预测方法

目前裂缝预测方法主要分为两类:第一类方法基于方位地震数据的振幅或振幅导出属性(频率、衰减等)预测裂缝,仅能得到界面两侧地层的综合响应,无法确定上覆、下伏地层的裂缝信息;第二类方法基于方位AVO叠前地震反演获得不同方位的弹性参数,可直接求取目的层内的裂缝信息,但待反演参数维度高,且各向异性参数值较弹性参数值小得多,反演稳定性有待提升。为此,提出了基于方位弹性阻抗傅里叶级数的裂缝预测方法。基于Rüger的方位AVO反射系数近似方程,首先推导了方位弹性阻抗方程,将裂缝介质的界面信息转化为地层内部的弹性信息,方位弹性阻抗差异体现了地层的各向异性;其次,利用傅里叶级数展开方位弹性阻抗方程求取二阶傅里叶系数,通过方位弹性阻抗与余弦函数的互相关求取裂缝法向。模型测试表明,基于方位弹性阻抗傅里叶级数展开的裂缝预测方法能较好地预测裂缝发育情况以及裂缝法向,精度较高且抗噪性较强。实际数据预测结果表明,裂缝预测结果与成像测井结果一致,具有较高精度。     

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

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

小尺度体反射系数近似公式研究

AVO技术经过几十年的发展已成为石油勘探中的一种重要技术,针对非均质性油气藏的储层预测,AV0方法具有其他方法不可替代的作用。应用波动方程差分数值解方法．研究了横向变速小尺度体的AVO特征,讨论了横向小尺度体的宽度和速度变化对反射系数的影响,然后应用最小二乘法原理拟合出了反射系数变化量与小尺度体宽度和速度的数学表达式,并结合测试数据进行了误差分析。结果表明：拟合公式与测试数据之间的误差小于0．5％,较好地反映了小尺度体宽度及岩石弹性参数对反射系数的影响。该公式可为利用AVO异常来预测小尺度体的宽度和岩石弹性参数奠定基础,同时也为下一步横向小尺度变速情况下反射系数精确表达式的验证提供保障。     

用切变模量近似P-SV波反射系数及在多波AVO反演中的意义

本文从精确的P－SV反射系数公式出发,推导出了作为射线参数函数的近似公式。该式说明,在非法线入射时,P－SV波反射系数主要取决于两个弹性参数：①反射界面两侧介质体密度之差与平均体密度之比△ρ／ρ;②反射界面两侧切变模量差与平均体密度之比△μ／ρ。文中还给出了数值实例,以说明该式的计算精度,同时分析了其适用范围,并与Aki和Richards（1980）的近似式的数值结果作了比较。本文采用阻尼最小平方反演理论计算出的反演参数的结果表明,在AVO资料线性反演中,估算△μ／ρ比估算泊松比差△σ精度高。