首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到16条相似文献,搜索用时 109 毫秒
1.
岩石的等效孔隙纵横比反演及其应用   总被引:4,自引:2,他引:2       下载免费PDF全文
通过融合Gassmann方程和由微分等效介质理论建立的干岩石骨架模型--DEM解析模型,本文提出根据纵波(和横波)速度反演岩石等效孔隙纵横比进行储层孔隙结构评价和横波速度预测的方法.首先,利用Gassmann方程和DEM解析模型建立岩石的纵、横波速度与密度、孔隙度、饱和度和矿物组分等各参数之间的关系;其次,将岩石孔隙等效为具有单一纵横比的理想椭球孔,应用非线性全局寻优算法来寻找最佳的等效孔隙纵横比使得理论预测与实际测量的弹性模量之间的误差最小;最后,将反演得到的等效孔隙纵横比代入到Gassmann方程和DEM解析模型中构建横波速度.实验室和井孔测量数据应用表明,反演得到的等效孔隙纵横比可准确反映储层的孔隙结构,对于裂缝型储层如花岗岩,其孔隙纵横比通常小于0.025,而对于孔隙型储层如砂岩,其孔隙纵横比通常大于0.08.只利用纵波与同时利用纵、横波反演得到的孔隙纵横比结果几乎完全一致,而且由纵波构建的横波与实测横波吻合良好,说明本文提出的等效孔隙纵横比反演及其横波速度预测方法是有效的.  相似文献   

2.
深层碳酸盐岩往往具有孔隙度低,孔隙形态复杂的特点.传统孔隙度反演方法仅考虑孔隙空间大小对地层岩石物性的影响,未考虑孔隙形态对岩石弹性特征的作用,在反演深层碳酸盐岩孔隙度时往往存在偏差.孔隙形态的影响可以利用含有孔隙纵横比的等效介质模型表示;通过引入DEM模型对声波孔隙度方程进行改进,定量体现孔隙形态的影响.将改进的声波反演与中子反演、密度反演相结合,得到井区受孔隙形态影响的地层孔隙度.首先利用中子、密度和纵波时差反演得到地层孔隙度初值;进而通过DEM公式获得孔隙纵横比初值;最后综合利用中子测井、密度测井和声波测井数据,实现受孔隙纵横比影响的地层孔隙度反演.应用该方法反演了西南某地MX19井深层白云岩储层的孔隙度,通过与岩心孔隙度测量数据进行对比,验证了反演结果的准确性.  相似文献   

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

4.
复杂孔隙储层往往同时发育孔缝洞等多种孔隙类型,这种孔隙结构的复杂性使得岩石的速度与孔隙度之间的相关性很差.经典的二维岩石物理模版只研究弹性参数与孔隙度和饱和度之间的定量关系,而不考虑孔隙结构的影响,用这样的模版来预测复杂孔隙储层的物性参数时带来很大偏差.本文首先证明多重孔隙岩石的干骨架弹性参数可以用一个等效孔隙纵横比的单重孔隙岩石物理模型来模拟;进而基于等效介质岩石物理理论和Gassmann方程,建立一个全新的三维岩石物理模版,用它来建立复杂孔隙岩石的弹性性质与孔隙扁度及孔隙度和饱和度之间的定量关系;在此基础上,预测复杂储层的孔隙扁度、孔隙度以及孔隙中所包含的流体饱和度.实际测井和地震反演数据试验表明,三维岩石物理模版可有效提高复杂孔隙储层参数的预测精度.  相似文献   

5.
石玉梅  谢韬 《地球物理学报》2017,60(12):4753-4765
地下固体介质弹性模量是固体地球定量表征的重要参数,是探测地球内部结构、地下流体分布、设计合理工程方案的关键之一.本文首先简要分析常用的流体饱和孔隙介质岩石物理模型,在此基础上通过数值计算,详细讨论了流体饱和多孔介质等效弹性模量对固体基质、孔隙度、孔隙结构和孔隙流体的敏感性,探讨了利用等效介质弹性模量反演固体基质弹性模量的可行性,建立了固体基质弹性模量的非线性方程系统,提出了该方程系统的数值迭代求解方法,并通过对三类岩石样本实验室测量数据的反演和结果分析,验证了方法的有效性.文中同时通过不同条件下对实验数据的反演,探讨了孔隙结构、孔隙流体和等效介质弹性模量对反演效果的影响,为方法的合理有效应用提供依据.  相似文献   

6.
孔隙纵横比是描述多孔岩石微观孔隙结构特征的重要参数,目前用于获取岩石完整孔隙纵横比分布的经典模型为David-Zimmerman(D-Z)孔隙结构模型,该模型假设岩石由固体矿物基质、一组纵横比相等的硬孔隙以及多组纵横比不等的微裂隙构成,并认为固体矿物基质和硬孔隙均不受压力影响,在此基础上,利用超声纵横波速度的压力依赖性反演岩石硬孔隙和各组微裂隙的孔隙纵横比及孔隙度.该方法的关键点在于以累积裂隙密度为桥梁,借助等效介质理论建立了岩石弹性模量和孔隙纵横比之间的内在联系.但在D-Z模型中,多重孔隙岩石累积裂隙密度的计算直接由单重孔隙裂隙密度公式实现,这种近似导致该模型在许多情况下难以获得良好的反演精度.为了完善经典D-Z模型,本文提出了一种基于虚拟降压的孔隙纵横比分布反演策略,通过多个假想降压过程实现累积裂隙密度的准确计算,并将基于DEM和MT的经典D-Z模型推广到KT和SCA中,结合四种等效介质理论建立了一套完整的反演流程.采用一系列砂岩和碳酸盐岩样品,测试了反演流程在实际岩芯孔隙纵横比提取中的应用效果,研究结果表明:与D-Z模型相比,本文方法的模拟结果与实际数据吻合更好,并同时适用于砂岩和碳酸盐岩;此外,通过分析四种等效介质理论的模拟结果发现,本文方法并不十分依赖于等效介质理论的选择,这些理论获得的孔隙结构参数随压力的变化趋势基本一致,数值上仅存在略微差异,且这种差异随着压力的增大逐渐消失.本文方法是经典D-Z孔隙结构模型的重要补充,对岩石孔隙结构表征、流体饱和岩石速度预测以及孔间喷射流效应的模拟具有十分重要的意义.  相似文献   

7.
碳酸盐岩组成复杂,储集层空间类型多,储层参数的规律性较差,AVO响应特征弱,尤其是孔隙形状对弹性参数的影响不可忽略.针对碳酸盐岩储层的特点,基于自洽模型和微分等效模型构建了多孔隙结构条件下多相介质的碳酸盐岩岩石物理模型.按照各相介质体积含量由小到大的顺序,先利用自洽模型计算孔隙度为50%时的双相等效介质,再通过微分等效模型将孔隙度调整至实际孔隙度,计算实际孔隙条件下的双相介质等效体,重复迭代该过程,最终构建出与实际岩石等效的介质体,计算获得等效岩石的弹性模量参数.通过实验对比了自洽模型、微分等效模型与自洽-微分等效模型的岩石弹性模量与孔隙度、孔隙纵横比之间的适用性规律.最后基于四川龙岗飞仙关组实际测井资料,结合自洽-微分等效模型进行了目标地层的横波预测.结果显示该模型的横波预测效果相对Xu-Payne模型预测,两者预测曲线比较接近,但在形态上更接近实际横波.SCA-DEM模型在孔隙度和孔隙纵横比上具有更好的适应性,为碳酸盐岩孔洞型储层的岩石物理建模提供了一种新参考.  相似文献   

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

9.
相对于常规砂岩,致密砂岩在岩石物理性质、力学性质等方面具有明显差异,并呈现出很强的非均质性.岩石物理模型能将储层参数与地震特性信息联系起来,因此可以作为致密砂岩储层参数与地震特性信息转换的桥梁.常规的岩石物理模型通常只考虑单一因素(例如非均匀性,单一孔隙,单一尺度等),建立的岩石物理模板并不适用于致密砂岩.本文针对高饱和气、微裂隙发育、非均质性强的致密砂岩储层,利用Voigt-Reuss-Hill模型计算混合矿物的弹性模量,采用微分等效介质(DEM)模型描述含裂隙、孔隙岩石的骨架弹性模量,基于Biot-Rayleigh波动方程构建了岩石物理弹性模板,给出了致密砂岩储层弹性参数与物性的关系.基于测井数据和实验数据对岩石物理弹性模板进行校正,并将校正后的岩石物理弹性模板结合叠前地震资料应用于川西地区储层孔隙度与裂隙含量预测.结果显示,反演裂隙含量、孔隙度与储层试气报告、测井孔隙度基本吻合,表明该模板能够较合理地应用于致密砂岩储层孔隙度及裂隙含量解释中.  相似文献   

10.
针对低孔渗储层渗透率主要受孔隙结构影响致使利用常规孔隙度和束缚水饱和度等参数预测渗透率精度低的难题,依据低孔渗岩石孔隙中流体渗流的特点,考虑岩石孔隙空间中对流体渗流贡献最大的那一部分孔隙,引入有效流动孔隙概念,以提高低孔渗岩石渗透率计算精度.考虑岩性、物性、电性变化,设计岩石物理实验,根据压汞实验数据计算岩样对应不同孔隙半径的进汞孔隙度与渗透率之间的相关系数,制作相关系数与孔隙半径交会图,将相关系数达到某一值(如0.8)对应的孔隙半径确定为有效流动孔隙半径下限,其对应的进汞孔隙度确定为有效流动孔隙度实验分析值.依据水流与电流流动相似性原理,从导电角度推导有效流动孔隙度的计算公式.对于含水低孔渗泥质岩石,将束缚水和粘土水看成不导电干骨架,采用能够描述孔喉比的等效岩石元素模型推导出只有可动水孔隙存在的低孔渗岩石地层因素的公式;根据有效流动孔隙的含义,只有可动水孔隙存在的低孔渗岩石的有效流动孔隙可以等效为具有相同岩石体积和地层因素且由骨架和弯曲毛细管组成的岩石结构中的毛管孔隙,根据弯曲毛管模型推导出地层因素与有效流动孔隙度关系式,结合两式可得到有效流动孔隙度的表达式.根据实验数据采用最优化方法确定公式中参数,从而给出有效流动孔隙度计算式.统计有效流动孔隙度实验值与岩心分析渗透率关系,回归得到低孔渗泥质岩石渗透率的计算式.对B区X1、X2井低孔渗储层进行处理,从实际应用效果方面证实该方法提高了低孔渗泥质储层渗透率计算精度.  相似文献   

11.
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.  相似文献   

12.
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.  相似文献   

13.
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.  相似文献   

14.
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.  相似文献   

15.
田志 《地球物理学报》2021,64(3):1119-1130
油藏岩石的孔隙连通性是反映流体渗流难易程度的重要参数,对渗透率、有效孔隙度等岩石物理参数的评价具有重要作用.连通的孔隙中,核磁共振(NMR)弛豫的交换会产生扩散耦合现象,可作为孔隙连通性的表征和探测方法.本文提出利用横向弛豫T 2-T 2脉冲序列测量岩石的扩散耦合现象.运用随机游走方法模拟多孔岩石的核磁共振响应特征,分析扩散耦合的影响因素,推导表征扩散耦合强度的弛豫交换速率计算公式.结果表明:孔隙间的扩散耦合强度与T 2-T 2脉冲序列的混合时间呈正相关性,基于双孔弛豫交换模型推导的弛豫交换速率计算公式能够准确表征双尺度孔隙系统的扩散耦合强度.在孔隙尺寸不满足快扩散条件时,会出现与扩散耦合无关的非对角峰信号.针对含多类型孔隙的碳酸盐岩模型,随混合时间的增加,扩散耦合强度变大,一维T 2谱的形态畸变程度加重,在T 2-T 2二维谱中,代表微裂缝、粒间小孔、溶蚀大孔的信号能量变化趋势不同,反映不同类型孔隙间的连通性存在差异.本文的分析与讨论丰富了核磁共振弛豫在岩石物理性质评价中的应用方向,对利用核磁共振评价复杂孔隙岩石的孔隙结构和连通性提供了新思路和新方法.  相似文献   

16.
Quantifying the effects of pore-filling materials on elastic properties of porous rocks is of considerable interest in geophysical practice. For rocks saturated with fluids, the Gassmann equation is proved effective in estimating the exact change in seismic velocity or rock moduli upon the changes in properties of pore infill. For solid substance or viscoelastic materials, however, the Gassmann theory is not applicable as the rigidity of the pore fill (either elastic or viscoelastic) prevents pressure communication in the pore space, which is a key assumption of the Gassmann equation. In this paper, we explored the elastic properties of a sandstone sample saturated with fluid and solid substance under different confining pressures. This sandstone sample is saturated with octadecane, which is a hydrocarbon with a melting point of 28°C, making it convenient to use in the lab in both solid and fluid forms. Ultrasonically measured velocities of the dry rock exhibit strong pressure dependency, which is largely reduced for the filling of solid octadecane. Predictions by the Gassmann theory for the elastic moduli of the sandstone saturated with liquid octadecane are consistent with ultrasonic measurements, but underestimate the elastic moduli of the sandstone saturated with solid octadecane. Our analysis shows that the difference between the elastic moduli of the dry and solid-octadecane-saturated sandstone is controlled by the squirt flow between stiff, compliant, and the so-called intermediate pores (with an aspect ratio larger than that of compliant pore but much less than that of stiff pores). Therefore, we developed a triple porosity model to quantify the combined squirt flow effects of compliant and intermediate pores saturated with solid or viscoelastic infill. Full saturation of remaining stiff pores with solid or viscoelastic materials is then considered by the lower embedded bound theory. The proposed model gave a reasonable fit to the ultrasonic measurements of the elastic moduli of the sandstone saturated with liquid or solid octadecane. Comparison of the predictions by the new model to other solid substitution schemes implied that accounting for the combined effects of compliant and intermediate pores is necessary to explain the solid squirt effects.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号