共查询到20条相似文献，搜索用时 31 毫秒
1.
V. Shulakova J. Sarout L. Pimienta M. Lebedev S. Mayo M.B. Clennell M. Pervukhina 《Geophysical Prospecting》2017,65(1):251-265
CO2 geosequestration is an efficient way to reduce greenhouse gas emissions into the atmosphere. Carbonate rock formations are one of the possible targets for CO2 sequestration due to their relative abundance and ability to serve as a natural trapping reservoir. The injected supercritical CO2 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 CO2 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 CO2. A statistical analysis of the microtomographic images showed that the injection of supercritical CO2 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 CO2 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 CO2 injection/residence are lower than those for the intact sample. 相似文献
2.
Application of alternative digital rock physics methods in a real case study: a challenge between clean and cemented samples 
下载免费PDF全文
下载免费PDF全文 Sadegh Karimpouli Sadegh Khoshlesan Erik H. Saenger Hamed Hooshmand Koochi 《Geophysical Prospecting》2018,66(4):767-783
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. 相似文献
3.
碳酸盐岩复杂的孔隙结构如何影响其弹性性质一直是地球物理研究的难点问题,在此基础上如何半定量甚至是定量地对碳酸盐岩储层预测,特别是如何有效地获取孔隙结构参数相关的地震属性体一直是油气工业界追求的目标,本研究从数字岩心角度入手,联合测井以及地震数据尝试探究这一问题的解决方案.首先针对代表不同孔隙结构类型的有限数目的碳酸盐岩样品获得其对应的高精度数字岩心数据体,为了获得更加可靠的具有地球物理含义的弹性性质随孔隙度变化的统计规律,我们通过子网格的技术,在有限数目的碳酸盐岩数字岩心数据体上获得了大量的数字岩心子网格样本,对于每个子网格样本可以分别获得其对应的数字岩心图像孔隙度、表征孔隙软硬程度的孔隙结构参数(γ)、以及基于有限元法模拟的弹性性质,由此基于数字岩心的研究思路,我们最终获得了基于孔隙结构因子表征与分类下的弹性性质与孔隙度的定量化解释量版.与此同时,在地震尺度上通过叠前地震资料获取的纵横波及密度属性体后,基于如上获得的定量化解释量版,我们最终获得了针对碳酸盐岩储层的新的属性体——孔隙结构参数(γ)属性体,这使得在地震尺度上预测碳酸盐岩储层的孔隙结构类型成为可能,也使利用地震数据在孔隙结构参数表征与分类下的碳酸盐岩储层反演精度的提高成为可能. 相似文献
4.
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. 相似文献
5.
Analytical approximations of bulk and shear moduli for dry rock based on the differential effective medium theory 总被引：2，自引：0，他引：2
Differential effective medium theory has been applied to determine the elastic properties of porous media. The ordinary differential equations for bulk and shear moduli are coupled and it is more difficult to obtain accurate analytical formulae about the moduli of dry porous rock. In this paper, in order to decouple these equations we first substitute an analytical approximation for the dry‐rock modulus ratio into the differential equation and derive analytical solutions of the bulk and shear moduli for dry rock with three specific pore shapes: spherical pores, needle‐shaped pores and penny‐shaped cracks. Then, the validity of the analytical approximations is tested by integrating the full differential effective medium equation numerically. The analytical formulae give good estimates of the numerical results over the whole porosity range for the cases of the three given pore shapes. These analytical formulae can be further simplified under the assumption of small porosity. The simplified formulae for spherical pores are the same as Mackenzie's equations. The analytical formulae are relatively easy to analyse the relationship between the elastic moduli and porosity or pore shapes and can be used to invert some rock parameters such as porosity or pore aspect ratio. The predictions of the analytical formulae for experimental data show that the formulae for penny‐shaped cracks are suitable to estimate the elastic properties of micro‐crack rock such as granite, they can be used to estimate the crack aspect ratio while the crack porosity is known and also to estimate the crack porosity evolution with pressure if the crack aspect ratio is given. 相似文献
6.
Ida L. Fabricius Christian Høier Peter Japsen Uffe Korsbech 《Geophysical Prospecting》2007,55(4):487-506
In impure chalk, the elastic moduli are not only controlled by porosity but also by contact‐cementation, resulting in relatively large moduli for a given porosity, and by admixtures of clay and fine silica, which results in relatively small moduli for a given porosity. Based on a concept of solids suspended in pore fluids as well as composing the rock frame, we model P‐wave and S‐wave moduli of dry and wet plug samples by an effective‐medium Hashin–Shtrikman model, using chemical, mineralogical and textural input. For a given porosity, the elastic moduli correspond to a part of the solid (the iso‐frame value) forming the frame of an Upper Hashin–Shtrikman bound, whereas the remaining solid is modelled as suspended in the pore fluid. The iso‐frame model is thus a measure of the pore‐stiffness or degree of cementation of the chalk. The textural and mineralogical data may be assessed from logging data on spectral gamma radiation, density, sonic velocity and water saturation in a hydrocarbon zone, whereas the iso‐frame value of a chalk may be assessed from the density and acoustic P‐wave logs alone. The iso‐frame concept may thus be directly used in conventional log‐analysis and is a way of incorporating sonic‐logging data. The Rigs‐1 and Rigs‐2 wells in the South Arne field penetrate the chalk at the same depth but differ in porosity and in water saturation although almost the entire chalk interval has irreducible water saturation. Our model, combined with petrographic data, indicates that the difference in porosity is caused by a higher degree of pore‐filling cementation in Rigs‐1. Petrographic data indicate that the difference in water saturation is caused by a higher content of smectite in the pores of Rigs‐1. In both wells, we find submicron‐size diagenetic quartz. 相似文献
7.
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. 相似文献
8.
碳酸盐岩复杂的孔隙结构如何影响其弹性性质一直是地球物理研究的难点问题,在此基础上如何半定量甚至是定量地对碳酸盐岩储层预测,特别是如何有效地获取孔隙结构参数相关的地震属性体一直是油气工业界追求的目标.本研究从数字岩心角度入手,联合测井以及地震数据尝试探究这一问题的解决方案,包括如下几个方面:(1)代表性碳酸盐岩储层样品获取;(2)CT扫描数字岩心数据体获取;(3)数字岩心数据的图像处理;(4)数字岩心数据的静态弹性模拟;(5)数字岩心子数据体的孔隙结构因子提取;(6)孔隙结构因子表征与分类下的弹性性质与孔隙度的定量化量版建立;(7)数字岩心-井-地震联合的孔隙度属性提取;(8)孔隙结构因子的地震属性体获取.本研究分为两篇系列文章上篇与下篇,上篇主要阐述如上提出的(1)—(4)方面,重点在于针对碳酸盐岩二值化图像处理的流程建立与验证,以及数字岩心静态弹性模拟的理论方面,这两方面是基于数字岩心获得精确的碳酸盐岩弹性性质模拟结果的关键所在;下篇主要阐述利用数字岩心数据获得孔隙结构因子的思路、理论与流程,以及为碳酸盐岩储层预测为目标而获得孔隙结构因子的地震属性体的实际应用方面.由于两篇文章共享数字岩心数据,同时所涉及的研究思路与流程形成一个有机整体,因此写成两篇系列文章而非两篇独立文章.本文为两篇系列文章的第一篇:上篇. 相似文献
9.
地震尺度下碳酸盐岩储层的岩石物理建模方法(英文) 总被引：3，自引：3，他引：0
碳酸盐岩油藏的强非均质性以及孔隙结构的复杂性,使得作为连接油藏参数与地震参数重要桥梁的岩石物理模型,以及作为油藏预测和定量表征最有效工具的流体替换成为岩石物理建模的难点与重点。在碳酸盐岩储层复杂孔隙结构与地震尺度下碳酸盐岩储层非均质性分析基础上,研究采用岩石网格化方法,将地震尺度下非均质碳酸盐岩储层岩石划分为具有独立岩石参数的均质岩石子体,根据岩石孔隙成因与结构特征采用不同岩石物理模型分步计算岩石子块干岩石弹性模量,并根据不同孔隙连通性进行流体替换,计算饱和不同流体岩石弹性模量。基于计算的岩石子块弹性模量,采用Hashin-Shtrikman-Walpole弹性边界计算理论方法实现地震尺度下碳酸盐岩储层弹性参数计算。通过对含有不同类型孔隙组合碳酸盐岩储层模型的弹性模量进行计算与分析,明确不同孔隙对岩石弹性参数的影响特征,模拟分析结果与实际资料认识一致。 相似文献
10.
The dependence of elastic moduli of shales on the mineralogy and microstructure of shales is important for the prediction of sweet spots and shale gas production. Based on 3D digital images of the microstructure of Longmaxi black shale samples using X-ray CT, we built detailed 3D digital images of cores with porosity properties and mineral contents. Next, we used finite-element (FE) methods to derive the elastic properties of the samples. The FE method can accurately model the shale mineralogy. Particular attention is paid to the derived elastic properties and their dependence on porosity and kerogen. The elastic moduli generally decrease with increasing porosity and kerogen, and there is a critical porosity (0.75) and kerogen content (ca. ≤3%) over which the elastic moduli decrease rapidly and slowly, respectively. The derived elastic moduli of gas- and oil-saturated digital cores differ little probably because of the low porosity (4.5%) of the Longmaxi black shale. Clearly, the numerical experiments demonstrated the feasibility of combining microstructure images of shale samples with elastic moduli calculations to predict shale properties. 相似文献
11.
Differential effective medium (DEM) theory is applied to determine the elastic properties of dry rock with spheroidal pores. These pores are assumed to be randomly oriented. The ordinary differential equations for bulk and shear moduli are coupled and it is more difficult to obtain accurate analytical formulae about the moduli of dry porous rock. In this paper, we derive analytical solutions of the bulk and shear moduli for dry rock from the differential equations by applying an analytical approximation for dry-rock modulus ratio, in order to decouple these equations. Then, the validity of this analytical approximation is tested by integrating the full DEM equation numerically. The analytical formulae give good estimates of the numerical results over the whole porosity range. These analytical formulae can be further simplified under the assumption of small porosity. The simplified formulae for spherical pores (i.e., the pore aspect ratio is equal to 1) are the same as Mackenzie's equations. The analytical formulae are relatively easy to analyze the relationship between the elastic moduli and porosity or pore shapes, and can be used to invert some rock parameters such as porosity or pore aspect ratio. The predictions of the analytical formula for the sandstone experimental data show that the analytical formulae can accurately predict the variations of elastic moduli with porosity for dry sandstones. The effective elastic moduli of these sandstones can be reasonably well characterized by spheroidal pores, whose pore aspect ratio was determined by minimizing the error between theoretical predictions and experimental measurements. For sandstones the pore aspect ratio increases as porosity increases if the porosity is less than 0.15, whereas the pore aspect ratio remains relatively stable (about 0.14) if the porosity is more than 0.15. 相似文献
12.
碳酸盐岩孔隙结构参数构建与储层参数反演(英文) 总被引：3，自引：0，他引：3
碳酸盐岩储层孔隙结构相对碎屑岩更复杂,常用的岩石物理模型不能较好的描述其孔隙结构的变化规律,且岩石孔隙结构的差异较大程度上会影响岩石的弹性性质。本文首先利用岩石薄片分析了碳酸盐岩的微观孔隙结构。然后基于Gassmann方程和Eshelby-Walsh椭球包体裂缝理论,在合理的假设前提下给出了一种新的岩石物理建模方法,并且从中提取了一个参数来表征孔隙结构的变化规律。最后,基于全波列测井数据,我们利用该方法计算了单井的孔隙度,并与用常规方法预测的结果进行了比较,同时进行了地震储层参数反演。研究结果表明,孔隙结构对岩石的弹性性质的影响较大,且新的建模方法预测的孔隙度误差仅为0.74%。因此,该方法可有效的减小孔隙结构对计算各岩石弹性参数的影响并提高孔隙度的预测精度。 相似文献
13.
An approach to determining the effective elastic moduli of rocks with double porosity is presented. The double‐porosity medium is considered to be a heterogeneous material composed of a homogeneous matrix with primary pores and inclusions that represent secondary pores. Fluid flows in the primary‐pore system and between primary and secondary pores are neglected because of the low permeability of the primary porosity. The prediction of the effective elastic moduli consists of two steps. Firstly, we calculate the effective elastic properties of the matrix with the primary small‐scale pores (matrix homogenization). The porous matrix is then treated as a homogeneous isotropic host in which the large‐scale secondary pores are embedded. To calculate the effective elastic moduli at each step, we use the differential effective medium (DEM) approach. The constituents of this composite medium – primary pores and secondary pores – are approximated by ellipsoidal or spheroidal inclusions with corresponding aspect ratios. We have applied this technique in order to compute the effective elastic properties for a model with randomly orientated inclusions (an isotropic medium) and aligned inclusions (a transversely isotropic medium). Using the special tensor basis, the solution of the one‐particle problem with transversely isotropic host was obtained in explicit form. The direct application of the DEM method for fluid‐saturated pores does not account for fluid displacement in pore systems, and corresponds to a model with isolated pores or the high‐frequency range of acoustic waves. For the interconnected secondary pores, we have calculated the elastic moduli for the dry inclusions and then applied Gassmann's tensor relationships. The simulation of the effective elastic characteristic demonstrated that the fluid flow between the connected secondary pores has a significant influence only in porous rocks containing cracks (flattened ellipsoids). For pore shapes that are close to spherical, the relative difference between the elastic velocities determined by the DEM method and by the DEM method with Gassmann's corrections does not exceed 2%. Examples of the calculation of elastic moduli for water‐saturated dolomite with both isolated and interconnected secondary pores are presented. The simulations were verified by comparison with published experimental data. 相似文献
14.
B. Gurevich 《Pure and Applied Geophysics》2000,157(5):811-827
15.
岩石孔隙结构是控制砂岩和碳酸盐岩的地震波速度和渗透率的重要参数之一。如果两种类似的岩石其孔隙度一定,而渗透率不同,那么它们的声波速度相差2km/s,而渗透率两者则可能相差近6个数量级,即从0.01mD到20mO。在本文中我们总结了由一个广义孔隙弹性理论简化的一个双参数弹性速度模型，以描述孔隙结构对弹性波的影响。由于矿物和储层的流体是给定的,我们用孔隙度和骨架的柔性(挠性)因子来确定速度模型,这个模型可以用于地震反演和储层表征,已改善孔隙度和储量的计算骨架柔性因子可用于岩石结构(PST)类型的定量分类,并可以利用叠前、叠后的两种地震资料将其与孔隙的连通性和渗透率联系起来。本项研究同时也有助于说明为什么振幅与偏移距(AVO)分析用于流体检测在某些情况下失败。这是由于孔隙结构对地震波的影响能够掩盖所有流体效应，特别是在碳酸岩中。 相似文献
16.
A fundamental understanding of flow in porous media at the pore-scale is necessary to be able to upscale average displacement processes from core to reservoir scale. The study of fluid flow in porous media at the pore-scale consists of two key procedures: Imaging - reconstruction of three-dimensional (3D) pore space images; and modelling such as with single and two-phase flow simulations with Lattice-Boltzmann (LB) or Pore-Network (PN) Modelling. Here we analyse pore-scale results to predict petrophysical properties such as porosity, single-phase permeability and multi-phase properties at different length scales. The fundamental issue is to understand the image resolution dependency of transport properties, in order to up-scale the flow physics from pore to core scale. In this work, we use a high resolution micro-computed tomography (micro-CT) scanner to image and reconstruct three dimensional pore-scale images of five sandstones (Bentheimer, Berea, Clashach, Doddington and Stainton) and five complex carbonates (Ketton, Estaillades, Middle Eastern sample 3, Middle Eastern sample 5 and Indiana Limestone 1) at four different voxel resolutions (4.4 µm, 6.2 µm, 8.3 µm and 10.2 µm), scanning the same physical field of view. Implementing three phase segmentation (macro-pore phase, intermediate phase and grain phase) on pore-scale images helps to understand the importance of connected macro-porosity in the fluid flow for the samples studied. We then compute the petrophysical properties for all the samples using PN and LB simulations in order to study the influence of voxel resolution on petrophysical properties. We then introduce a numerical coarsening scheme which is used to coarsen a high voxel resolution image (4.4 µm) to lower resolutions (6.2 µm, 8.3 µm and 10.2 µm) and study the impact of coarsening data on macroscopic and multi-phase properties. Numerical coarsening of high resolution data is found to be superior to using a lower resolution scan because it avoids the problem of partial volume effects and reduces the scaling effect by preserving the pore-space properties influencing the transport properties. This is evidently compared in this study by predicting several pore network properties such as number of pores and throats, average pore and throat radius and coordination number for both scan based analysis and numerical coarsened data. 相似文献
17.
Maxim Lebedev Moyra E.J. Wilson Vassily Mikhaltsevitch 《Geophysical Prospecting》2014,62(6):1253-1265
Petrophysical properties of carbonate reservoirs are less predictable than that of siliciclastic reservoirs. One of the main reasons for this is the physical and chemical interactions of carbonate rocks with pore fluids. Such interactions can significantly change the elastic properties of the rock matrix and grains, making the applicability of Gassmann's fluid substitution procedure debatable. This study is an attempt to understand the mechanisms of fluid‐rock interactions and the influence of these interactions on elastic parameters of carbonates. We performed precise indentation tests on Savonnières limestone at a microscale level under dry, distilled water, and n‐Decane saturated conditions. Our experiments display softening of the rock matrix after water saturation. We have found that mainly the ooid cortices, peloid nuclei and prismatic intergranular cement are affected by water flooding. We also observed a shear modulus reduction in Savonnières limestone in an experiment performed at ultrasonic frequencies. One of the most important results obtained in our experimental study is that the Gassmann fluid substitution theory might not always be applicable to predict the elastic moduli of fluid‐saturated limestones. 相似文献
18.
19.
碳酸盐岩储层孔隙类型多样,各种孔隙的尺寸变化范围可以跨越几个数量级,孔隙结构非常复杂,这种复杂孔隙结构和不均匀分布的多元孔隙空间使得储层电性呈现明显非阿尔奇特性.为了了解影响电阻率变化的控制因素,本次研究选取中三叠世雷口坡组的8块全直径碳酸盐岩岩样,开展了核磁共振、岩电实验、孔渗实验、压汞实验及薄片等实验,并利用数字图像分析法定量分析了孔隙结构特征.研究结果表明：①孔隙度是影响电阻率高低的重要因素,但并非唯一因素,除孔隙度以外,孔隙尺寸和数量、孔隙网络复杂程度远比吼道大小对电阻率的影响大;②在孔隙度一定的条件下,胶结指数m随储层中孤立大孔隙占比的增多而增大,当孔隙度增大到一定程度后,胶结指数m又随大孔隙占比的增多而减小,微裂缝起重要沟通作用;③在给定孔隙度时,以简单大孔隙为主的岩样表现为胶结指数m值较大,而以复杂孔隙网络、细小孔隙为主的岩样表现为胶结指数m值较小,具分散、孤立大孔隙的岩样,胶结指数m值最高;④依据孔隙几何参数与电阻率和胶结指数之间的关系,可以利用测井资料间接判别储层类型,从而提高储层有效性和含水饱和度评价精度. 相似文献
20.
本文对渤海湾盆地不同孔隙结构样品的孔、渗、核磁、岩电、压汞、X衍射及铸体薄片等配套岩石物理实验数据进行了综合分析,通过逐一考察同一套岩芯样品的地层因素与渗透率、压汞喉径均值、储层品质指数之间的实验关系,并分别与地层因素-孔隙度交绘图进行对比分析,发现储层渗透性与地层因素之间并非简单的单调函数关系,孔隙度相近但孔隙结构类型不同、渗透率差异明显的岩芯可以具有相近的地层因素,导电能力接近.在实验数据分析的基础上通过理论分析证明了这一实验关系的合理性,并指出孔隙度及导电能力相近的岩芯,其渗透率差异与喉径均值的平方比、孔隙曲折度及几何形态相关. 相似文献