超深碳酸盐岩油气储层岩石弹性性质高温高压超声实验研究

面向超深碳酸盐岩油气勘探开发中储层精细地震刻画问题,研发了高温(200℃)、高压(200 MPa)岩石弹性性质超声测量系统.收集塔里木盆地9块奥陶系超深碳酸岩储层岩石样品,通过高温高压条件下超声纵、横波速度测量发现,在围压保持69.0 MPa不变的情况下,当温度从25℃上升到200℃时,干岩石样品的纵、横波速度以及体积模量与剪切模量随着温度的升高而下降,纵波速度与体积模量最大降幅分别达到3.54%和7.54%,横波速度与剪切模量最大降幅分别达到4.57%和8.93%;在温度保持200℃不变情况下,当围压从69.0 MPa上升到137.9 MPa时,干岩石样品的纵、横波速度以及体积模量与剪切模量随着围压的升高而升高,纵波速度与体积模量最大升幅分别达到1.78%和3.39%,横波速度与剪切模量最大升幅分别达到1.95%和5.26%.在超深碳酸盐岩油气储层环境下,温度对碳酸盐岩岩石骨架弹性性质的影响可能超过围压对岩石骨架弹性性质的影响.岩石样品温度与围压循环实验表明,干岩石弹性性质与温度以及围压的变化关系不具可逆性,碳酸盐岩岩石样品弹性模量与温度之间总体上存在线性关系,其截距与斜率由碳酸盐岩...     

基于保幅波场分离的弹性波逆时偏移方法研究（英文）

欲实现基于弹性波方程的矢量波场逆时偏移纵、横波独立成像,必须在波场延拓过程中实现纵、横波场的分离,散度和旋度算子分离的纵、横波出现振幅与相位的畸变,导致输出成像结果的振幅失真。本文提出一种在弹性波场延拓过程中实现纵、横波保幅分离的方法,在传统的弹性波方程中加入纵波压力、纵波振动速度和横波振动速度方程,实现纵横波的矢量分解,再对分解后的矢量纵波和矢量横波做标量化合成得到保幅分离的纵、横波场,对保幅分离的纵、横波场应用成像条件,然后实现矢量波场逆时偏移的保幅纵横波成像。该方法可以保证分离后纵、横波的振幅与相位不变;同时,分解后的纵波压力和纵波振动速度可用于层间反射噪音压制和横波极性校正,提高多分量地震资料联合逆时偏移的纵、横波成像质量,从而实现保幅弹性波逆时偏移的目的,为叠前深度剖面应用于叠前反演工作奠定基础。     

基于贝叶斯框架的各向异性页岩储层岩石物理反演技术

页岩岩石物理建模旨在建立页岩矿物组分、微观结构、流体填充与岩石弹性参数的关系.对四川盆地龙马溪组页岩进行岩石物理建模研究,针对页岩黏土含量高、层间微裂缝发育等特点,利用Backus平均理论描述页岩黏土矿物弹性参数,利用Chapman理论计算与水平微裂缝有关的VTI各向异性,并利用Bond变换考虑地层倾角的影响.提出以黏土矿物纵、横波速度和孔隙纵横比为拟合参数进行岩石物理反演的方法,并引入贝叶斯框架减小反演的多解性.由已知的黏土矿物纵、横波速度和孔隙纵横比作为先验信息,并以测井纵、横波速度作为约束条件建立反演的目标函数,同时利用粒子群算法进行最优化搜索.计算结果表明,基于先验约束和粒子群算法的反演方法能够较准确地反演黏土矿物的弹性参数、孔隙形态参数以及裂缝密度等参数.计算得到的黏土纵、横波速度较高,并且在一定范围内变化,这可能与龙马溪组页岩的黏土矿物组分中具有较高弹性模量的伊利石含量较高有关,同时也与黏土定向排列等微观物性特征有关.反演得到的裂缝密度与纵波各向异性参数ε呈明显的正相关,而与横波各向异性参数γ相关性较小.另外,页岩各向异性参数与黏土垂向的纵横波速度有较强的相关性.     

基于累积法理论的弹性阻抗参数提取

弹性阻抗中蕴含着丰富的地下介质信息。从弹性阻抗中准确地提取地下介质信息用于储层预测是十分有意义的,传统的提取方法主要基于最小二乘理论的线性方法。本文引入累积法理论,将其用于弹性阻抗参数的提取,用层状模型和部分Marmousi2模型对该方法进行测试。从背景噪声分别为正态分布和均匀分布的弹性阻抗体中提取纵、横波速度和密度,并将其与最小二乘法提取的纵、横波速度和密度做定量对比分析。结果表明两种方法所得结果精度相当,但累积法具有计算量小、简洁直观的优势,因此,将其用于弹性参数的提取是有效可行的。      

基于散度和旋度纵横波分离方法的改进

纵、横波的分离是多波多分量地震资料处理中很重要的一步,其分离结果直接影响到后续数据处理的质量.各向同性介质中纵波为无旋场,横波为无散场,因此可以在频率-波数域利用散度和旋度算子对地震记录进行纵、横波分离,但是此处理过程必须知道地表处的纵、横波速度.本文给出了一种估算地表纵、横波速度的方法,可以在纵、横波速度值未知的情况下,将其估算出来.针对弹性波场进行散度和旋度运算时,纵、横波的相位和振幅比发生改变的问题,本文给出了相位和纵、横波振幅比的校正方法.     

基于纵横波模量两项反射系数近似方程的叠前地震反演

与叠后波阻抗反演相比,叠前AVO反演能获得较为丰富的信息,因此利用叠前反演提取地层弹性参数并进行流体识别是目前的研究热点.考虑到基于Aki-Richard近似式的三参数反演不稳定性较突出,文章拟从减少参数维数入手,提高反演结果的稳健性.本文首先从宗兆云等人推导的基于纵横波模量的Zeoppritz近似公式出发,合理地对公式中的纵波模量与密度进行幂指数拟合,推导出直接反演纵波模量(M)、横波模量(U)的两项式反射系数近似方程并进行反演.其次,建立纵横波模量与高灵敏度流体识别因子的联系,并求取高灵敏度流体识别因子.模型试算和实例应用均表明,新方法能较好地识别储层内的流体性质.     

气饱和超压页岩的岩石物理建模研究

页岩的岩石物理建模随着页岩气勘探的深入而逐渐受到重视.现有的建模方法对孔隙压力的研究并不全面.本文认为页岩气对页岩弹性的影响可分为两个方面.第一个方面称为静态效应,定义为气体压力对页岩骨架的作用.第二方面称为动态效应,定义为页岩气本身对页岩弹性的作用.通过将孔隙纵横比与有效应力的关系、气体弹性参数与孔隙压力的关系嵌入岩石物理建模流程,实现静态和动态效应的作用.孔隙纵横比与有效应力的关系通过实测纵、横波速度标定,用于预测不同压力状态下页岩的纵、横波速度.本文方法虽存在一定的局限性,但提供了一种分析页岩弹性参数与孔隙压力相互关系的方法.     

基于一阶速度-胀缩-旋转方程的纵横波地震偏移角度域共成像点道集提取方法（英文）

角度域共成像点道集(angle-domain common-image gathers,ADCIGs)是多波AVA反演和偏移速度分析的基础数据。纵、横波坡印廷矢量的准确求取是获得多分量地震资料ADCIGs的前提。由常规速度-应力弹性波方程得到的坡印廷矢量代表了纵、横两种波的混合波场的传播方向,而非纯纵波或纯横波的传播方向。一阶速度-胀缩-旋转弹性波方程能够求取纯纵波或纯横波的坡印廷矢量,获得单一类型波的传播方向信息。本文基于一阶速度-胀缩-旋转弹性波方程逆时偏移实现纵、横波ADCIGs的求取,具体思路为:利用胀缩振动速度矢量和胀缩标量分别求取炮点和接收点波场中纯纵波的坡印廷矢量,利用剪切振动速度矢量和旋转矢量求取纯横波的坡印廷矢量,取炮点波场中各网格点纵波能量最大时刻的纵波坡印廷矢量为入射纵波传播方向矢量,取接收点波场中各网格点纵(横)波能量最大时刻的纵(横)波坡印廷矢量为反射纵波(反射转换横波)传播方向矢量,据此求取入射角并最终得到纵、横波角道集。模型试算表明:本文方法可以准确求取界面处入射波和反射波的传播方向矢量,得到地震波入射角,实现纵、横波角道集的高精度提取。     

基于弹性波阻抗的拉梅参数反演与应用

(Connolly的弹性阻抗(EI)公式是纵、横波速度、密度和入射角的函数,应用该公式的常规反演方法只能直接得到纵、横波速度和密度的信息,而且不同角度的弹性阻抗在数值上差别很大,不在同一尺度上,无法进行对比。本文针对常规反演方法的缺点,提出了基于Gray的Zoeppritz近似方程的新的弹性阻抗公式,并对这个公式进行了标准化以实现不同角度的弹性阻抗在数值上的统一,最后用标准化后的公式进行了反演,利用拉梅参数与弹性阻抗问的线性关系可从反演得到的弹性阻抗数据体中直接提取得到拉梅参数。应用实例表明用这种方法提取得到的弹性参数更加稳定、准确,而且能很好地反映储层信息。这种新的方法是对以Connolly公式为基础的传统方法的改进。     

岩石弹性波速度和饱和度、孔隙流体分布的关系

在实验室对6种砂岩进行了连续的进水和失水实验. 测 量了两个过程的纵波速度VP、横波速度VS进水和失水速率随饱和度SW的 变化；分析了低饱和度时流体对岩石弹性性质的影响. 实验表明，进水和失水过程显示不同 的纵、横波速度与饱和度关系，速度不仅与饱和度有关，还与不同饱和阶段的孔隙流体分布 有关，而且也是水和岩石骨架之间的物理及化学作用所致.     

基于Xu-White模型横波速度预测的误差分析

在砂泥岩剖面中,基于Gassmann方程和Xu-White岩石物理模型,利用常规测井资料及实验室岩石物理数据进行横波速度估算.模型中考虑了泥质砂岩中基质性质、泥质含量、孔隙度大小和孔隙形状以及孔隙饱含流体性质对岩石速度的影响.综合分析了模型中砂、泥岩和孔隙流体弹性参数以及孔隙纵横比等输入参数的误差对预测横波速度精度的影响.数值模型试验表明,在Xu-White模型中采用变化的孔隙纵横比估算出的横波速度远远比采用固定孔隙纵横比估算出的横波速度准确.     

云南省9个地震台台基岩石的力学性质

应用岩石力学实验方法测定了云南省9个地震台台基岩石的力学性质。其中包括密度、强度、静态扬氏模量、静态泊松比、纵波速度及阻抗、横波速度及阻抗、动态扬氏模量、动态剪切模量和动态泊松比。讨论了地震台台基岩石力学性质参数在测震学中的作用。     

三分量地震在薄互层油气藏勘探中的应用（英文）

本文是三分量地震在陆相薄互层地震勘探中应用成功的首例。大庆长垣喇嘛甸油田处于高含水开发后期,储层预测的主要难题是密井网条件下厚度2m以上砂体的边界识别问题,常规单一纵波地震解释存在多解性,引入横波信息可增加地震解释的可靠性。本文根据纵、横波对气藏的响应特征,利用纵、横波联合振幅属性、分频属性对气藏进行了分析,能够很好地反应气藏的边界,体现了横波受流体影响较小的优势。通过测井曲线岩石物理定量的分析了该地区主要是密度、λρ和μρ三种物性参数能够很好地反映储层岩性变化。在以地震数据分频为核心的高精度地震参考标准层解释基础上利用纵、横波联合反演得出的纵、横波速度比、λρ、μρ等属性与测井解释的储层具有很好的一致性,多波属性λρ沿层切片描述了开发小层的砂体平面展布,与井砂岩图具有较高的一致性,刻画了井间砂体边界变化的细节,并指示了剩余油挖潜的有利区域。     

中低孔渗储层岩石弹性参数的流体敏感性研究(英文)

本文从WXS凹陷中低孔渗储层岩石声波实验出发,以岩样的纵横波速度和密度为基础数据,求取出一系列的弹性参数,包括纵横波波速比、纵波波阻抗、横波波阻抗、泊松比、拉梅常数、剪切模量、体积模量、杨氏模量,等等。在前人的孔隙流体识别究基础上,综合相关理论和实验分析,构建了一个新的流体识别因子F。以饱和流体岩石弹性参数及其组合参数的相对变化量Ag/w和Ao/w为定量指标,评价各流体识别因子的流体识别效果,并采用交会图技术进行了验证。新流体敏感因子在传统较难分辨的孔隙流体"水"和"油"的区别上具有良好效果,有利于提高中低孔渗储层流体识别的成功率。     

纵横波弹性阻抗联立反演在GD地区的应用

在GD油田复杂油藏描述中应用叠前纵横波弹性阻抗反演,精确地进行了油藏岩性的划分。利用三个或三个以上部分叠加数据,进行纵横波弹性阻抗联立反演,既克服了因叠后地震反演结果单一而不能满足复杂储层描述的需求,又避免了由于叠前道集信噪比低造成反演结果不稳定的缺陷。本文论述了叠前弹性波阻抗反演的基本原理,结合GD地区实际资料,对反演过程中涉及的角道集子波提取、层位标定、横波速度预测、弹性参数提取与解释等关键步骤进行了详细研究,指出基于测井资料分析的多种弹性参数综合解释是提高叠前地震反演应用效果的关键。     

岩石物理弹性参数规律研究

根据辽东湾凹陷某区在地层条件和不同流体相态(气饱和、水饱和等)下岩石纵波速度、横波速度及密度等岩心测试数据,以及岩石矿物成分、孔隙度等常规岩心分析数据,统计分析了岩石弹性参数变化规律.采用有效流体模型、斑块饱和模型进行了纵、横波速度理论计算,并和实验测量结果比较,认为高孔、高渗岩石可以看作有效流体模型,低孔、低渗岩石更接近斑块饱和模型.这些规律和认识对于指导储层预测和油气检测及地震振幅综合解释有重要的意义.     

巷道前方煤体波速与瓦斯参数相关分析研究

为了获取掘进巷道前方煤体地震波速度和瓦斯参数间的对应关系,以阳泉五矿为实施矿井.采用矿井地震波速度超前探测方法对掘进煤巷进行探测,通过初至波拾取获得地震波原始记录中纵、横波速,并与巷道前方30m范围内平均瓦斯参数作相关分析.分析结果表明:一定条件下,煤体瓦斯含量越高,其地震波纵、横波速度越低.运用曲线拟合对纵、横波视速度与瓦斯绝对涌出量和瓦斯浓度进行拟合发现:纵、横波视速度与瓦斯绝对涌出量间存在ΔP=-4.6415Vp-2.7048Vs+16.1353的负线性关系;纵、横波视速度与瓦斯浓度间存在η=-0.9302Vp-0.5704Vs+3.3132的负线性关系,并对线性拟合结果进行验证表明:在误差允许范围内,计算值与实测值基本吻合.最后计算得到适用于太原组15#煤层及变质程度更高煤层的视速度瓦斯预警指标:Vp=2.00m/ms和Vs=0.73m/ms.     

部分叠加资料纵横波速度测井约束反演

本文提出一种角度部分叠加资料同时反演纵波和横波速度的方法,以角度部分叠加资料为基础,利用地震波振幅随入射角变化与弹性参数间的数学关系,基于非线性最优化理论,通过将最小平方问题转化为大型带状矩阵的求解问题,采用逐道外推技术依次求得该角度剖面的每个点纵波和横波速度值,利用地震波速度在各向异性介质中与入射角的关系,进一步得到地震波零入射角纵、横波速度,从而可以求得纵横波速度比及泊松比,为地震资料的岩性及含油气性解释提供了丰富的参数信息,为保证反演的稳定性,降低多解性,采用测井、叠后构造层位解释结果的共同约束,并做好叠前去噪以及振幅保持处理.     

基于应力偏量的横波纯应力逆时偏移成像方法

基于弹性波解耦延拓方程的波场分离方法不仅可以得到解耦的纵、横波质点振动速度场,还可以获得纵、横波应力场.针对横波纯应力场在利用单一分量成像时不具有明确物理意义的问题,本文将应力偏张量引入到横波应力场中,基于应力偏量第二不变量构建得到横波应力不变量并将其用于应力场的逆时偏移成像中,获得了可完整表征横波应力场的成像结果.模型试算表明,本文构建的横波应力不变量可以有效利用横波应力张量中的波场信息,并得到准确的弹性逆时偏移成像结果.     

基于纵横波分离FCT弹性波正演频散压制

有限差分方法被广泛应用于地震波数值模拟和传播.传统有限差分法采用Taylor级数展开实现空间偏导数的差分,但该方法会因为网格离散化而产生数值频散,降低地震波模拟的精度.优化差分系数正演方法能在一定程度上压制部分频散,然而纵、横波速度取值差异较大,在弹性波有限差分正演模拟中,在满足纵波最大速度确定的稳定性条件下,浅层低速横波波场往往会产生明显的频散现象.为了削弱弹性波场正演数值频散,提高数值模拟精度,本文首先采用优化差分网格系数降低数值频散,然后再采用通量校正传输(Flux-Correction Transport, FCT)法来进一步压制弹性波场有限差分数值频散.常规的FCT法是对弹性波场直接进行频散压制,但由于弹性波场中纵、横波速度差异明显,横波波场频散明显强于纵波,为了压制横波波场的数值频散,往往需要选取较大的频散压制参数,但这会使频散较弱的纵波产生假象.因此本文提出基于纵横波分离FCT弹性波正演频散压制方法,对分离之后的纵横波场分别选择合适的频散压制参数进行通量校正,可以有效压制数值频散,削弱纵波FCT产生的假象.通过理论分析和数值算例发现,本文方法能有效削弱弹性波场有限差分数值...