基于数字岩心的碳酸盐岩孔隙结构对弹性性质的影响研究(下篇):储层孔隙结构因子表征与反演

碳酸盐岩复杂的孔隙结构如何影响其弹性性质一直是地球物理研究的难点问题,在此基础上如何半定量甚至是定量地对碳酸盐岩储层预测,特别是如何有效地获取孔隙结构参数相关的地震属性体一直是油气工业界追求的目标,本研究从数字岩心角度入手,联合测井以及地震数据尝试探究这一问题的解决方案.首先针对代表不同孔隙结构类型的有限数目的碳酸盐岩...     

基于弹性阻抗的孔隙结构与物性参数非线性同步反演

碳酸盐岩储层具有复杂的孔隙结构,其显著影响了岩石的弹性和地震响应特征,也影响物性参数预测精度.现有的岩石物理反演方法主要是针对碎屑岩储层,反演的对象主要是孔隙度和饱和度,采用的数据主要是弹性参数,使用的算法主要是线性近似法,缺少针对碳酸盐岩储层的孔隙结构与物性参数反演方法.针对这些问题,文章首次提出了一种基于弹性阻抗的适用于碳酸盐岩复杂孔隙储层的孔隙结构和物性参数贝叶斯非线性同步反演新方法,该方法把多孔介质岩石物理模型、Gassmann方程、AVO理论、贝叶斯理论和非线性反演算法结合起来,实现复杂孔隙储层的孔隙结构和物性参数的定量同步预测.模型正演揭示,孔隙结构参数即孔隙扁度对AVO反射系数和弹性阻抗的贡献度仅次于孔隙度,远远大于含水饱和度的影响.实际资料应用表明,文章提出的从叠前道集中直接反演物性参数和等效孔隙扁度新方法可以准确预测储层孔隙度和饱和度及其空间展布,并可评价有效储层的孔隙结构.     

地震尺度下碳酸盐岩储层的岩石物理建模方法(英文)

碳酸盐岩油藏的强非均质性以及孔隙结构的复杂性,使得作为连接油藏参数与地震参数重要桥梁的岩石物理模型,以及作为油藏预测和定量表征最有效工具的流体替换成为岩石物理建模的难点与重点。在碳酸盐岩储层复杂孔隙结构与地震尺度下碳酸盐岩储层非均质性分析基础上,研究采用岩石网格化方法,将地震尺度下非均质碳酸盐岩储层岩石划分为具有独立岩石参数的均质岩石子体,根据岩石孔隙成因与结构特征采用不同岩石物理模型分步计算岩石子块干岩石弹性模量,并根据不同孔隙连通性进行流体替换,计算饱和不同流体岩石弹性模量。基于计算的岩石子块弹性模量,采用Hashin-Shtrikman-Walpole弹性边界计算理论方法实现地震尺度下碳酸盐岩储层弹性参数计算。通过对含有不同类型孔隙组合碳酸盐岩储层模型的弹性模量进行计算与分析,明确不同孔隙对岩石弹性参数的影响特征,模拟分析结果与实际资料认识一致。     

微晶结构对碳酸盐岩地震弹性与储层物性特征变化规律的影响（英文）

碳酸盐岩储层是油气勘探的重要储层类型,不同地区的碳酸盐岩储层在沉积、成岩过程都存在着明显的差异,这种差异是引起其物性及地震弹性性质变化的的主要因素。为认识这种机制,在本文中我们对取自塔里木盆地不同地区碳酸盐岩储层不同层位的98块样品进行了岩石学特征、岩石微观结构特征、物性特征和地震弹性特征识别分析,并研究了样品物性对地震弹性性质的影响及其变化规律。研究结果表明,碳酸盐岩样品物性与地震弹性性质整体变化规律受微晶方解石结构特征的控制,而传统的岩石结构及孔隙结构分类不能完全反应上述特征的变化。依据形态特征可将样品微晶结构分为多孔微晶、紧密微晶体和致密微晶三种类型,随着上述三种微晶中晶体边界紧密缝合接触的程度增加,其微晶晶间孔隙度与孔喉半径逐渐减小,同时方解石晶体颗粒边界刚度特征与弹性均匀性的逐渐增强,致使样品渗透率与速度等地震弹性特征表现出随孔隙增大而逐渐减小的总体趋势。对于以致密微晶为主的低孔隙度样品(孔隙度φ5%),微晶孔隙对样品孔隙及渗透率贡献较小,同时微晶弹性性质接近岩石基质,致使岩石宏观物性特征与弹性性质受裂隙、溶蚀孔隙等的影响更为明显,样品物性与地震弹性性质的特征表征孔隙结构的控制作用。孔隙纵横比对孔隙类型也具有明显的指示作用,溶蚀孔隙的体积模量纵横比高于0.2,粒间孔体积模量纵横比介于0.1~0.2之间,多孔微晶与紧密微晶体积模量纵横比小于0.15,而致密微晶的体积模量纵横比接近于0.2。研究结果可为碳酸盐岩储层的分类以及相关储层的岩性及烃类地震检测提供依据。     

碳酸盐岩储层各向异性岩石物理建模与孔隙结构分析（英文）

针对非均质性较强的碳酸盐岩储层,综合各向异性自相容近似理论与差分等效介质理论建立了一种各向异性岩石物理弹性参数计算模型。通过对碳酸盐岩孔隙结构实验结果分析,提出了碳酸盐岩弹性参数计算的孔隙结构模型。综合使用各向异性等效介质理论与碳酸盐岩孔隙结构模型,给出了碳酸盐岩非均质储层弹性参数计算的实现方法。通过对比灰岩样品测试数据与理论模型计算结果,表明本文提出的等效弹性理论模型能够描述碳酸岩盐储层速度与孔隙度变化关系,可作为碳酸盐岩储层地震数据解释的基本依据。     

致密碳酸盐岩跨频段岩石物理实验及频散分析

致密碳酸盐岩在成岩和后成岩过程中形成了复杂的孔隙结构特征,其速度等地震弹性参数不仅与孔隙度有关,而且还与孔隙结构特征密切相关.为了进一步研究致密碳酸盐岩内部流体相关的速度频散特征,针对致密碳酸盐岩进行实验室的频散测量与频散理论分析尤为重要.本研究选用了一块典型的致密型碳酸盐岩样品,在对样品进行了精细的包括CT扫描与镜下薄片的孔隙结构描述基础上,进行了实验室跨频段(从地震频段-超声频段)的频散测量与频散响应分析.比较实验室跨频段岩石物理频散测量可以获得如下认识:1)较之于典型的"喷射流"机制,改进的"喷射流"模型可以半定量地解释频散测量的结果,这大大提高了对致密碳酸盐岩频散响应的理解与认识;2)改进的"喷射流"模型还不能完全精确地匹配实验室频散测量结果,这说明除了微观尺度下的"喷射流"机制,还存在着其他控制频散与衰减的机制;3)本项工作对研究致密碳酸盐岩储层中不同频段地震波响应以及对储层预测与流体识别提供了理论依据.     

基于碳酸盐岩孔隙结构预测孔隙度方法研究

许多研究都已经证实在碳酸盐岩储层中孔隙结构对声波速度影响很大,因此在孔隙度反演时必须考虑孔隙结构的影响.本文通过对Gassmann方程的合理简化并引入Eshelby-Walsh干燥岩石椭球包体近似公式,推导出包含岩石孔隙结构参数、饱和岩石压缩系数与岩石基质压缩系数三个参数的计算孔隙度的新公式,该式岩石基质压缩系数是通过Russell流体因子和Gassmann-Boit-Geer tsma方程计算式建立目标函数反演获得的,然后通过测井交汇图技术把岩石孔隙结构参数与岩石基质压缩系数优选转换成纵横波速度与密度关系式,进而导出具体地区考虑碳酸盐岩孔隙结构孔隙度具体计算公式,测井与地震资料应用表明,基于碳酸盐岩孔隙结构的孔隙度预测方法的精度高于常规方法.     

碳酸盐岩孔隙结构参数构建与储层参数反演(英文)

碳酸盐岩储层孔隙结构相对碎屑岩更复杂,常用的岩石物理模型不能较好的描述其孔隙结构的变化规律,且岩石孔隙结构的差异较大程度上会影响岩石的弹性性质。本文首先利用岩石薄片分析了碳酸盐岩的微观孔隙结构。然后基于Gassmann方程和Eshelby-Walsh椭球包体裂缝理论,在合理的假设前提下给出了一种新的岩石物理建模方法,并且从中提取了一个参数来表征孔隙结构的变化规律。最后,基于全波列测井数据,我们利用该方法计算了单井的孔隙度,并与用常规方法预测的结果进行了比较,同时进行了地震储层参数反演。研究结果表明,孔隙结构对岩石的弹性性质的影响较大,且新的建模方法预测的孔隙度误差仅为0.74%。因此,该方法可有效的减小孔隙结构对计算各岩石弹性参数的影响并提高孔隙度的预测精度。     

岩石孔隙结构的地震信号特征

岩石孔隙结构是控制砂岩和碳酸盐岩的地震波速度和渗透率的重要参数之一。如果两种类似的岩石其孔隙度一定,而渗透率不同,那么它们的声波速度相差2km/s,而渗透率两者则可能相差近6个数量级,即从0.01mD到20mO。在本文中我们总结了由一个广义孔隙弹性理论简化的一个双参数弹性速度模型，以描述孔隙结构对弹性波的影响。由于矿物和储层的流体是给定的,我们用孔隙度和骨架的柔性(挠性)因子来确定速度模型,这个模型可以用于地震反演和储层表征,已改善孔隙度和储量的计算骨架柔性因子可用于岩石结构(PST)类型的定量分类,并可以利用叠前、叠后的两种地震资料将其与孔隙的连通性和渗透率联系起来。本项研究同时也有助于说明为什么振幅与偏移距(AVO)分析用于流体检测在某些情况下失败。这是由于孔隙结构对地震波的影响能够掩盖所有流体效应，特别是在碳酸岩中。     

基于岩石导电效率理论建立碳酸盐岩储层饱和度的计算方法(英文)

由于孔隙型碳酸盐岩储层的复杂孔隙结构和强非均质性,岩石导电效率与含水孔隙度之间理想的线性关系相对复杂或并不存在,故前人基于岩石导电效率理论建立含水饱和度的计算公式并不完全适用于碳酸盐岩储层。基于岩石导电效率理论,推导了岩石导电效率的计算公式,阐明了岩石导电效率与含水孔隙度之间线性关系相对复杂或并不存在的根本原因,发现了岩石导电效率与电阻率之间的幂函数关系,分析了所建立含水饱和度计算公式误差的主控因素,得出了岩石导电效率的计算精度是该方法是否能推广应用的关键。与Archie公式相比,基于岩石导电效率理论建立的含水饱和度计算公式能更准确地计算研究靶区孔隙型碳酸盐岩储层的含水饱和度。在伊拉克某油田和印尼某气田3口井碳酸盐岩储层中的应用,表明当计算岩石导电效率的相对误差不大于0.1时,所计算储层含水饱和度的绝对误差不大于0.1,基本满足孔隙型碳酸盐岩储层精细评价的需求。     

Nonlinear simultaneous inversion of pore structure and physical parameters based on elastic impedance

Carbonate reservoirs have complex pore structures, which not only significantly affect the elastic properties and seismic responses of the reservoirs but also affect the accuracy of the prediction of the physical parameters. The existing rockphysics inversion methods are mainly designed for clastic rocks, and the inversion objects are generally porosity and water saturation. The data used are primarily based on the elastic parameters, and the inversion methods are mainly linear approximations. To date, there has been a lack of a simultaneous pore structure and physical parameter inversion method for carbonate reservoirs. To solve these problems, a new Bayesian nonlinear simultaneous inversion method based on elastic impedance is proposed. This method integrates the differential effective medium model of multiple-porosity rocks, Gassmann equation,Amplitude Versus Offset(AVO) theory, Bayesian theory, and a nonlinear inversion algorithm to achieve the simultaneous quantitative prediction of the pore structure and physical parameters of complex porous reservoirs. The forward modeling indicates that the contribution of the pore structure, i.e., the pore aspect ratio, to the AVO response and elastic impedance is second only to that of porosity and is far greater than that of water saturation. The application to real data shows that the new inversion method for determining the pore structure and physical parameters directly from pre-stack data can accurately predict a reservoir's porosity and water saturation and can evaluate the pore structure of the effective reservoir.     

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

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

Quantitative geophysical pore‐type characterization and its geological implication in carbonate reservoirs

This paper discusses and addresses two questions in carbonate reservoir characterization: how to characterize pore‐type distribution quantitatively from well observations and seismic data based on geologic understanding of the reservoir and what geological implications stand behind the pore‐type distribution in carbonate reservoirs. To answer these questions, three geophysical pore types (reference pores, stiff pores and cracks) are defined to represent the average elastic effective properties of complex pore structures. The variability of elastic properties in carbonates can be quantified using a rock physics scheme associated with different volume fractions of geophysical pore types. We also explore the likely geological processes in carbonates based on the proposed rock physics template. The pore‐type inversion result from well log data fits well with the pore geometry revealed by a FMI log and core information. Furthermore, the S‐wave prediction based on the pore‐type inversion result also shows better agreement than the Greensberg‐Castagna relationship, suggesting the potential of this rock physics scheme to characterize the porosity heterogeneity in carbonate reservoirs. We also apply an inversion technique to quantitatively map the geophysical pore‐type distribution from a 2D seismic data set in a carbonate reservoir offshore Brazil. The spatial distributions of the geophysical pore type contain clues about the geological history that overprinted these rocks. Therefore, we analyse how the likely geological processes redistribute pore space of the reservoir rock from the initial depositional porosity and in turn how they impact the reservoir quality.     

Numerical estimation of carbonate rock properties using multiscale images

Characterizing the pore space of rock samples using three‐dimensional (3D) X‐ray computed tomography images is a crucial step in digital rock physics. Indeed, the quality of the pore network extracted has a high impact on the prediction of rock properties such as porosity, permeability and elastic moduli. In carbonate rocks, it is usually very difficult to find a single image resolution which fully captures the sample pore network because of the heterogeneities existing at different scales. Hence, to overcome this limitation a multiscale analysis of the pore space may be needed. In this paper, we present a method to estimate porosity and elastic properties of clean carbonate (without clay content) samples from 3D X‐ray microtomography images at multiple resolutions. We perform a three‐phase segmentation to separate grains, pores and unresolved porous phase using 19 μm resolution images of each core plug. Then, we use images with higher resolution (between 0.3 and 2 μm) of microplugs extracted from the core plug samples. These subsets of images are assumed to be representative of the unresolved phase. We estimate the porosity and elastic properties of each sample by extrapolating the microplug properties to the whole unresolved phase. In addition, we compute the absolute permeability using the lattice Boltzmann method on the microplug images due to the low resolution of the core plug images. In order to validate the results of the numerical simulations, we compare our results with available laboratory measurements at the core plug scale. Porosity average simulations for the eight samples agree within 13%. Permeability numerical predictions provide realistic values in the range of experimental data but with a higher relative error. Finally, elastic moduli show the highest disagreements, with simulation error values exceeding 150% for three samples.     

基于固体替换模型的有机质含量对页岩弹性性质的影响分析

页岩中的TOC(Total Organic Carbon,总有机碳)含量,对页岩的有效弹性模量以及与之相关的弹性波速度(P波和S波)有重要影响,建立弹性模量与TOC含量关系是页岩气甜点预测的重要手段之一.CS和SM两种固体置换理论主要针对孔隙度较大的砂岩,能否适用于孔隙度低、孔隙形态复杂和非均质性强的页岩目前尚未深入研究.鉴于目前已知的富有机质页岩的TOC赋存形态与裂缝以及孔隙形态类似,有关TOC含量对岩石弹性模量的影响可视为孔隙物质充填问题来研究.本文利用数字岩心技术,构造同一数字岩心不同TOC含量的样本群,基于CS和SM两种固体替换理论模型,通过有限元(FEM)数值模拟交叉验证,详细研究了两种固体替换方程对页岩的适用性和TOC含量对页岩弹性性质的影响.研究表明,由于实际岩心孔隙及TOC分布的非均质性,CS替换方程弹性模量预测值与FEM模拟结果存在差异,而SM替换方程预测值与FEM模拟结果基本一致,两种方程的预测差异揭示页岩非均质强度,利用SM替换方程中的参数α_1,α_2,β_1和β_2可详细分析实际岩心孔隙及TOC分布的非均质特征.     

Rock physics modelling of porous rocks with multiple pore types: a multiple-porosity variable critical porosity model

A critical porosity model establishes the empirical relationship between a grain matrix and a dry rock by the concept of critical porosity. The model is simple and practical and widely used. But the critical porosity in the model is a fixed value that cannot relate to pore structure. The aim of this paper is to establish the theoretical relationship between critical porosity and pore structure by combining Kuster–Toksöz theory with the critical porosity model. Different from the traditional critical porosity model, critical porosity is not an empirical value but varied with pore shape and the ratio of bulk modulus versus shear modulus of the grain matrix. The substitution of the theoretical relationship into Kuster–Toksöz theory will generate the formulae for the bulk and shear moduli of multiple-porosity dry rocks, which is named the multiple-porosity variable critical porosity model. The new model has been used to predict elastic moduli for sandstone and carbonate rock. We compare the modelling results for P- and S-wave velocities and elastic moduli with the experimental data. The comparison shows that the new model can be used to describe the elastic properties for the rocks with multiple pore types.