Compressibility of porous rocks: Part I. Measurements of Hungarian reservoir rock samples

Pore volume compressibility is one of the physical properties of a reservoir that must be specified in many reservoir-engineering calculations. In the presented research, the effect of compact pressure, temperature and porosity on compressibility was investigated. A total of twenty-two different cores were tested: five limestone, one friable sandstone, fourteen medium to hard sandstone, and two very dense sandstone. Core samples were placed in the test cell and subject to compacting pressure up to 10,000 psi. Runs were made at room temperature and at 52°C for limestone samples. Although there were some publications concerning measurement and study of the effect of pressure and temperature on pore volume compressibility of reservoir rocks, nothing has been published about compressibility of Hungarian reservoir rocks, except of the work of Tóth and Bauer (1988). The present study showed pore volume compressibility data for different Hungarian fields. The result of the study at high temperature (52°C) shows that pore compressibility increases with increasing temperature.     

高渗水合物储层的孔隙可压缩性反演

孔隙可压缩性与水合物储层物性相关.由于海域天然气水合物埋藏较浅,沉积物尚未成岩,海底水合物储层处于固结和完全未固结之间的状态,通过岩石物理推导证明,这两种状态下体积模量之差仅与孔隙可压缩性相关,因此孔隙可压缩性反应了岩石颗粒从悬浮态到正常压实的成岩进程,它可能与储层束缚水饱和度及渗透率密切相关.本文利用多口典型井数据分析了孔隙可压缩性与渗透率的关系,结果表明在海底浅层沉积物中渗透率越高,孔隙可压缩性越小,孔隙可压缩性对高渗储层的判识能力明显强于其他参数.然后,本文建立了两种状态下叠前地震反射特征的差异与岩石孔隙可压缩性的联系,提出双状态叠前反演方法,综合利用叠前地震数据以及测井资料反演得到了岩石孔隙可压缩性.实际应用效果表明,孔隙可压缩性较好的预测了高渗透率地层,气烟囱、粗粒的高含砂层等含气流体疏导通道渗透率较高,同时为水合物形成提供有利条件.     

基于数字岩心的含裂隙储层砂岩介电性质研究

了解储层岩石的介电特性在石油工业的各个方面都有重要的应用.小尺度裂隙是影响岩石介电性质的地质因素之一,获得裂隙对含裂隙岩石介电性质影响的定量关系具有重要的理论和实践意义.以含裂隙人造砂岩的三维微观数字结构为基础,通过基于三维有限差分算法计算的岩石介电性质与实验数据的对比验证数值计算方法的有效性.在此基础上,通过理论模型获得不同孔隙度基质的介电性质,并在不含裂隙人造砂岩的三维微观数字结构中人为添加以裂隙密度和纵横比为定量表征参数的裂隙,应用验证后的数值算法模拟随频率变化的含裂隙砂岩的介电性质,分析和研究不同孔隙度基质中定向排列裂隙对砂岩介电性质的影响.结果表明,当裂隙孔隙度随裂隙纵横比或裂隙密度发生改变时,含裂隙砂岩的介电性质与裂隙密度以及裂隙纵横比呈正相关关系,而当裂隙孔隙度保持不变时,含裂隙砂岩的介电性质随裂隙纵横比的减小而增大;裂隙参数的改变对不同基质孔隙度的含裂隙砂岩的介电性质的影响趋势较为一致,但随着基质孔隙度的减小,裂隙对砂岩介电性质的影响逐渐增大.裂隙参数和基质孔隙度对含裂隙砂岩介电性质影响的研究结果为基于介电特性的裂缝性油气储层的定量表征提供了依据,在油气勘探开发中具有重要的应用前景.     

Compressional- and shear-wave velocities as a function of confining stress in dry sandstones

A laboratory study was carried out to investigate the influence of confining stress on compressional- and shear-wave velocities for a set of rock samples from gas-producing sandstone reservoirs in the Cooper Basin, South Australia. The suite of samples consists of 22 consolidated sublitharenites with helium porosity ranging from 2.6% to 16.6%. We used a pulse-echo technique to measure compressional- and shear-wave velocities on dry samples (cylindrical 4.6 × 2 cm) at room temperature and at elevated confining stress (≤ 60 MPa). Compressional- and shear-wave velocities in samples increase non-linearly with confining stress. A regression equation of the form V = A ? Be^?DP gives a good fit to the measured velocities with improved prediction of velocities at high confining stresses compared with equations suggested by other studies. The predicted microcrack-closure stresses of the samples show values ranging from 70 MPa to 95 MPa and insignificant correlation with porosity, permeability or clay content. There is a positive correlation between change in velocity with core porosity and permeability, but this association is weak and diminishes with increasing confining stress. Experimental results show that pore geometry, grain-contact type, and distribution and location of clay particles may be more significant than total porosity and clay content in describing the stress sensitivity of sandstones at in situ reservoir effective stress. The stress dependence of Cooper Basin sandstones is very large compared with data from other studies. The implication of our study for hydrocarbon exploration is that where the in situ reservoir effective stress is much less than the microcrack-closure stress of the reservoir rocks, the variation of reservoir effective stress could cause significant changes in velocity of the reservoir rocks. The velocity changes induced by effective stress in highly stress-sensitive rocks can be detected at sonic-log and probably surface-seismic frequencies.     

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

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

Reservoir parameter classification of a Miocene formation using a fractal approach to well logging,porosimetry and nuclear magnetic resonance

A methodology for rock classification is presented that considers lithology and reservoir parameters on the basis of a combined fractal analysis of well logs and mercury porosimetry results with nuclear magnetic resonance outcomes. A sandy‐shaly thinly‐bedded Miocene gas bearing formation in the Carpathian Foredeep is investigated. Fractal correlation dimensions D₂ calculated for standard logs are used to distinguish sandstone as the most homogeneous lithological group with the highest porosity. The fractal analysis also confirmed observations of gas accumulations in sandstone, shaly sandstone and sandy claystone as rocks of high porosity. The results of two laboratory methods are combined to improve reservoir properties assessment and evaluate movable media in pore space; this technique was based on the similarity of mercury porosimetry results plotted as cumulative intrusion volume versus pressure or pore diameter and also the curves of cumulative porosity and transverse relaxation time distributions using nuclear magnetic resonance. Close values of porosity from logs, recorded in situ and in laboratory measurements, provide the link between fractal analysis and porosimetry and nuclear magnetic resonance measurements.     

基于二维图像的数字岩心电导率计算方法研究

电导率是表征岩石电学性质的重要物理参数,在地质资源勘查和测井解释等领域发挥着巨大作用.快速、准确地确定岩石电导率具有重要的理论和实践意义.作为近年来发展的一种岩石物理数值模拟工具,数字岩心技术在定量计算电导率等物性参数方面应用广泛.三维微观结构的准确获取是数字岩心技术计算岩石电导率的关键,但传统获取岩石三维微观结构的方法较为复杂费时.为了方便快速地通过数字岩心技术计算岩石的电导率,本文研究了岩石二维与三维数字岩心的电导率联系.我们基于微米级X射线CT扫描得到的三个砂岩样品的微观结构信息建立了三维数字岩心,并通过有限元法计算的三维数字岩心电导率与实验数据的对比验证数值计算方法的有效性.随后我们数字地扩展了岩石的孔隙,产生了较大孔隙度的三维数字岩心样本,在此基础上,计算了三维数字岩心和相应二维数字岩心的电导率,并通过Archie公式分别拟合了电导率与孔隙度之间的关系,得到了相应的胶结系数.结果表明,三维数字岩心的胶结系数小于二维数字岩心的胶结系数,且二者的比值与岩石实测孔隙度呈线性负相关关系.以该联系为纽带,通过二维图像快速计算得到的电导率与孔隙度关系,确定了三维数字岩心的电导率与孔隙度关系,并进一步通过三维数字岩心的孔隙度计算其电导率.该方法计算得到的人工砂岩样品的电导率与其三维数字岩心电导率相关系数高于96%,验证了基于二维图像的数字岩心电导率计算方法的有效性.本文的研究结果为快速、准确地计算岩石电导率提供了新的思路,在油气勘探开发中有广阔的应用前景.     

多孔可变临界孔隙度模型及储层孔隙结构表征

碳酸盐岩、致密砂岩和页岩等储层具有孔隙类型多样、孔隙结构复杂和非均质性强等特征,属于典型的多重孔隙储层,孔隙结构表征是多重孔隙储层预测和流体识别的关键.现有的孔隙结构表征方法大多利用孔隙纵横比或者构建一种新参数来描述孔隙结构.岩石临界孔隙度模型是一种常用的岩石物理模型,具有一定的物理意义和地质含义.本文推导了岩石临界孔隙度与岩石孔隙结构(孔隙纵横比)之间的关系,进而利用极化(形状)因子建立临界孔隙度与弹性参数之间的关系,构建了能够包含多种孔隙类型的多孔可变临界孔隙度模型.利用多孔可变临界孔隙度模型由储层的弹性参数反演不同孔隙类型的体积含量.实验室测量数据和实际测井数据表明,多孔可变临界孔隙度模型能够适用于多重孔隙储层岩石物理建模和孔隙结构表征.     

利用热声发射技术测量岩石最高古温度的探索

根据岩石的热Kaiser效应,岩石能够记忆地质历史中经历过的最高温度.本次实验使用自行研制的岩石热声发射仪,利用人工加热的砂岩和灰岩样品证实了沉积岩存在热Kaiser效应,并探讨了根据热Kaiser效应测量沉积岩经历的最高温度的准确性.通过对塔里木盆地的系列埋深(2800~5300 m)的砂岩和泥岩样品及川西地区灰岩样...     

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

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

礁滩储层内部孔隙结构模型模拟与孔隙度预测

地下岩石是由岩石基质和孔隙流体组成的双相介质,其有效弹性参数受岩石基质、孔隙度、孔隙结构及孔隙流体的影响,因此为了得到孔隙度与岩石有效弹性参数之间的关系,必须消除其他因素对孔隙度的影响.本文首先引进等效体的概念和Eshelby椭球包体裂缝理论,然后在合理的假设前提下,运用Gassman流体替换方程,推导并建立了生物礁滩...     

POROSITY AND PORE STRUCTURE FROM ACOUSTIC WELL LOGGING DATA1

Wyllie's time-average equation and subsequent refinements have been used for over 20 years to estimate the porosity of reservoir rocks from compressional (P)-wave velocity (or its reciprocal, transit time) recorded on a sonic log. This model, while simple, needs to be more convincingly explained in theory and improved in practice, particularly by making use of shear (S)-wave velocity. One of the most important, although often ignored, factors affecting elastic velocities in a rock is pore structure, which is also a controlling factor for transport properties of a rock. Now that S-wave information can be obtained from the sonic log, it may be used with P-waves to provide a better understanding of pore structure. A new acoustic velocities-to-porosity transform based on an elastic velocity model developed by Kuster and Toksöz is proposed. Employing an approximation to an equivalent pore aspect ratio spectrum, pore structure for reservoir rocks is taken into account, in addition to total pore volume. Equidimensional pores are approximated by spheres and rounded spheroids, while grain boundary pores and flat pores are approximated by low aspect ratio cracks. An equivalent pore aspect ratio spectrum is characterized by a power function which is determined by compressional-and shear-wave velocities, as well as by matrix and inclusion properties. As a result of this more sophisticated elastic model of porous rocks and a stricter theory of elastic wave propagation, the new method leads to a more satisfactory interpretation and fuller use of seismic and sonic log data. Calculations using the new transform on data for sedimentary rocks, obtained from published literature and laboratory measurements, are presented and compared at atmospheric pressure with those estimated from the time-average equation. Results demonstrate that, to compensate for additional complexity, the new method provides more detailed information on pore volume and pore structure of reservoir rocks. Examples are presented using a realistic self-consistent averaging scheme to consider interactions between pores, and the possibility of extending the method to complex lithologies and shaly rocks is discussed.     

Application of alternative digital rock physics methods in a real case study: a challenge between clean and cemented samples

High‐resolution three‐dimensional images are used in digital rock physics to numerically compute rock physical properties such as permeability and elastic moduli. These images are not widely available, and their preparation is both expensive and time consuming. All of these issues highlight the importance of alternative digital rock physics methods that are based on two‐dimensional images and use different approaches to compute effective properties of three‐dimensional samples. In addition, the scale of study in both standard and alternative digital rock physics is very small, which applications of its results are questionable at wells or reservoir scale. The aim of this study is to use two‐dimensional images and alternative digital rock physics techniques for computing seismic wave velocity and permeability, which are compared with well and laboratory data. For this purpose, data from one well in a reservoir located in the southwestern part of Iran are used. First, two clean (carbonate) and two cemented (limy sandstone) samples were collected from well cores at different depths. Then, two‐dimensional images by scanning electron microscope and conventional microscope were captured. In the next step, two alternative digital rock physics methods, namely, empirical relations and conditional reconstruction, have been employed to compute P‐wave velocity and permeability of a three‐dimensional medium. Results showed that, in clean (mono‐mineral) samples, velocity values were reasonably close to well data. However, permeability values are underestimated compared with laboratory data because laboratory data were obtained at ambient pressure, whereas alternative digital rock physics results are more representative of reservoir pressure conditions. Nevertheless, permeability–porosity trends are valid for both samples. In the case of cemented samples, a two‐scale procedure, along with a method for two‐scale computation and grain‐cement segmentation, is presented and developed. Results showed that P‐wave velocity is overestimated probably due to random sampling in this method. However, velocity–porosity trends are in agreement with well data. Moreover, permeability results obtained for cemented samples were also similar to those obtained for the clean samples.     

基于沉积过程的数字岩石建模方法研究

本文提出模拟地层沉积及成岩过程的矿物沉积算法,建立数字岩石模型,并通过对比Micro-CT扫描图像和数值模型的局部孔隙度及平均渗流概率函数分布特征,评价建模的准确性.结果表明,由二维扫描提取的粒径信息作为输入参数,模拟矿物沉积过程建模得到的三维数字岩石模型,能够准确重构原始岩心的非均质性及渗流特性,成功应用于泥质砂岩、碳酸盐岩、页岩等存在多矿物或多尺度孔隙的数字岩石建模中.数字岩石物理是正在兴起的重要技术.数字岩石采用超高分辨率先进成像装备,采集和表征微纳尺度岩石结构,在岩石弹性、电性、核磁、渗流特性等数值计算中发挥重要作用.但是,由于三维直接成像在有限视域内难以表征足够的岩石非均质性,提取二维结构统计特征,利用统计或地质过程法重构具有代表性的三维岩石结构成为十分有价值的研究课题,而且,对业界大量存在的岩石薄片及电镜高清二维图像的深度开发应用也具有重要的现实意义.本文发展的新方法,复原沉积过程,较好地解决了孔隙尺度岩石物理定量研究中数值建模与理论计算的技术瓶颈.     

泥质含量及其分布形式对岩电关系影响的数字岩心仿真(英文)

由于泥质所造成的附加导电现象,泥质含量及其分布形式对电阻率增大系数I和含水饱和度Sw关系具有重要影响,由于岩石物理实验中岩心孔隙结构及其组分构成、分布的微观不可调性,因而泥质分布形式所造成的影响很难通过岩心实验来单独研究。基于数字岩心的格子气自动机方法是一种有效的微观数值模拟方法,本研究利用储层岩心薄片的骨架颗粒尺寸信息资料建立数字岩心模型,结合格子气自动机技术对数字岩心不同饱和流体情况下电的传输特性进行数值模拟研究,揭示了不同泥质含量和泥质分布形式对孔隙介质导电特性非阿尔奇现象产生的影响,建立饱和度指数和泥质含量之间的关系模型,其良好的吻合性表明该方法在岩石物理研究中是一种十分有效的研究方法,而新模型适于在非阿尔奇储层进行准确的饱和度评价。     

四川盆地侏罗系致密砂岩弹性特征及岩石物理建模

近年来围绕四川盆地侏罗系陆相致密砂岩已取得了勘探突破,其中川中—川西过渡带具备形成大气田的地质条件,但对该套致密砂岩弹性性质变化规律的研究还较少,致使利用地震方法进行"甜点"储层预测的精度不高.本文利用四川盆地侏罗系沙溪庙组32块样品开展了系统的声学测量,在此基础上,分析了样品弹性性质的变化规律.结合X射线衍射矿物组分分析、扫描电镜、铸体薄片和岩石薄片特征确定了不同成岩作用对岩石储集性能的影响.研究结果表明,研究区致密砂岩储层表现为孔隙型储层,受差异性成岩作用影响,黏土含量、钙质含量和硅质含量的差异以及它们分布特征之间的差异对岩石弹性性质造成了很大的影响.在研究区对岩石物性及弹性性质有明显影响的成岩作用包括早期的钙质胶结作用、压实作用和溶蚀作用,因此针对不同时期的成岩作用对岩石弹性及物性的影响,利用接触-胶结模型、微分等效模量模型和临界孔隙度校正的Hashin-Shtrikman上限模型建立了研究区致密砂岩的岩石物理模型.     

基于三相孔隙弹性理论的火山岩气层岩石物理响应特征分析

Unlike previous theories with velocity and/or elastic modulus averaging, we use a three-phase porous rock physics model developed by Santos for analyzing the seismic response of two immiscible fluids in saturated porous media. Considering reservoir reference pressure and coupling drag of two fluids in pores, the effects of frequency, porosity, and gas saturation on the phase velocities of the P-and S-waves are discussed in detail under field conditions. The effects of porosity and gas saturation on Vp/Vs are also provided. The data for our numerical experiments are from a sample of deep volcanic rock from Daqing. The numerical results show that the frequency dispersion effect can be ignored for deep volcanic rocks with low porosity and low permeability. It is concluded that for deep volcanic rocks the effect of gas content in pores on Vp/Vs is negligible but the effect of porosity is significant when there is a certain amount of water contained in the pores. The accurate estimate of lithology and porosity in this case is relatively more important.     

热处理对致密岩石物理性质的影响

致密气藏低孔低渗和超低含水饱和度等特征使其潜在水相圈闭损害严重,致密天然气产出表现为多尺度特征.选取泥页岩、致密砂岩和致密碳酸盐岩岩心,开展了100~600℃高温处理对岩心渗透率、孔隙度、重量、长度、直径和声速的影响实验.实验结果表明,碳酸盐岩、致密砂岩和泥页岩的热开裂阈值分别在300~400℃、300~500℃和500~600℃;高温处理后,岩心重量和密度降低,体积增加,泥页岩岩心孔隙度和渗透率提高幅度最显著,600℃处理后声波时差比常温时岩心声波时差提高了1.3倍.热处理消除了水相圈闭和粘土矿物膨胀损害,提高岩石孔隙度和渗透率,恢复或改善致密储层多尺度传质,有利于致密天然气资源开发,但同时高温使岩石破裂,扩展天然裂缝或产生新裂缝,导致工作液漏失,因此,热致裂给勘探开发致密天然气提出了机遇与挑战.     

Effect of supercritical CO2 on carbonates: Savonnières sample case study

CO₂ geosequestration is an efficient way to reduce greenhouse gas emissions into the atmosphere. Carbonate rock formations are one of the possible targets for CO₂ sequestration due to their relative abundance and ability to serve as a natural trapping reservoir. The injected supercritical CO₂ can change properties of the reservoir rocks such as porosity, permeability, tortuosity, and specific surface area due to dissolution and precipitation processes. This, in turn, affects the reservoir characteristics, i.e., their elastic properties, storage capacity, stability, etc. The tremendous progresses made recently in both microcomputed X‐ray tomography and high‐performance computing make numerical simulation of physical processes on actual rock microstructures feasible. However, carbonate rocks with their extremely complex microstructure and the presence of microporosity that is below the resolution of microcomputed X‐ray tomography scanners require novel, quite specific image processing and numerical simulation approaches. In the current work, we studied the effects of supercritical CO₂ injection on microstructure and elastic properties of a Savonnières limestone. We used microtomographic images of two Savonnières samples, i.e., one in its natural state and one after injection and residence of supercritical CO₂. A statistical analysis of the microtomographic images showed that the injection of supercritical CO₂ led to an increase in porosity and changes of the microstructure, i.e., increase of the average volume of individual pores and decrease in the total number of pores. The CO₂ injection/residence also led to an increase in the mean radii of pore throats, an increase in the length of pore network segments, and made the orientation distribution of mesopores more isotropic. Numerical simulations showed that elastic moduli for the sample subjected to supercritical CO₂ injection/residence are lower than those for the intact sample.