共查询到16条相似文献，搜索用时 140 毫秒
1.
波场在含流体的孔隙介质中传播时会产生频散和衰减现象.波场的频散和衰减与孔隙介质的岩石物理属性有关,包括孔隙度、渗透率、流体属性等.现有的三维裂缝/软孔隙网络模型利用椭圆截面纵横比的变化模拟从扁裂缝、软孔隙到硬孔隙的多种情况,而未考虑同时包含孔隙和裂缝的全局性网络空间.为了更好地描述裂缝-孔隙空间,本文提出同时包含裂缝和孔隙的三维裂缝-孔隙网络模型,并给出渗透率的计算方法.通过体积平均法推导了三维裂缝/软孔隙网络模型和三维裂缝-孔隙网络模型的波动方程,利用平面波分析方法得到纵波频散/衰减曲线的表达式,同时应用数值模拟研究了总孔隙度、裂缝孔隙度、裂缝纵横比、裂缝数密度、孔隙流体黏度对纵波衰减和速度频散特征的影响.结果表明,在三维裂缝-孔隙网络模型下,总孔隙度、裂缝参数等对纵波频散衰减特征的影响与三维裂缝/软孔隙网络模型相似.具体表现为:纵波在高频段内出现频散和衰减现象.孔隙度的变化主要影响逆品质因子曲线峰值的大小;裂缝数密度主要控制速度显著变化的范围;裂缝纵横比对纵波速度和特征频率有显著影响. 相似文献
2.
碳酸盐岩、致密砂岩和页岩等储层具有孔隙类型多样、孔隙结构复杂和非均质性强等特征,属于典型的多重孔隙储层,孔隙结构表征是多重孔隙储层预测和流体识别的关键.现有的孔隙结构表征方法大多利用孔隙纵横比或者构建一种新参数来描述孔隙结构.岩石临界孔隙度模型是一种常用的岩石物理模型,具有一定的物理意义和地质含义.本文推导了岩石临界孔隙度与岩石孔隙结构(孔隙纵横比)之间的关系,进而利用极化(形状)因子建立临界孔隙度与弹性参数之间的关系,构建了能够包含多种孔隙类型的多孔可变临界孔隙度模型.利用多孔可变临界孔隙度模型由储层的弹性参数反演不同孔隙类型的体积含量.实验室测量数据和实际测井数据表明,多孔可变临界孔隙度模型能够适用于多重孔隙储层岩石物理建模和孔隙结构表征. 相似文献
3.
致密砂岩气藏具有裂缝发育和有效应力高的特征,研究不同有效压力下孔、裂隙介质地震波传播特征,有利于地震解释与地下储层的识别.但是前人的研究较少考虑岩石内部微观孔隙结构特征与孔隙、裂隙间流体流动的关系.本文首先通过选取四川盆地典型致密砂岩岩样,在不同有效压力下对岩石样本进行超声波实验测量.然后基于实验测得的纵、横波速度进行裂隙参数反演,得到不同有效压力下致密砂岩样本的裂隙孔隙度.再将裂隙孔隙度和样本岩石物理参数代入双重孔隙介质模型,模拟得到不同有效压力下饱水致密砂岩样本纵横波速度频散和衰减的变化规律.结果表明模型预测的速度频散曲线与纵波速度实验测量结果能够较好的吻合.最后统计分析了致密砂岩裂隙参数,得到了致密砂岩储层裂隙参数随有效压力及孔隙度变化特征.依据实际岩石物理参数建立模型,其裂隙参数三维拟合结果能够较好描述致密砂岩裂隙结构与孔隙度、应力的关联,可为实际地震勘探中预测储层裂缝性质提供基础依据. 相似文献
4.
储层砂岩微观孔隙结构特征不仅影响干燥岩石的弹性波传播速度,也决定了岩石介质中与流体流动相关的速度频散与衰减作用.依据储层砂岩微观结构特征及速度随有效压力变化的非线性特征,将其孔隙体系理想化为不同形状的硬孔隙(纵横比α0.01)与软孔隙(纵横比α0.01)的组合(双孔隙结构).基于孔弹性理论,给出软孔隙最小初始纵横比值(一定压力下所有未闭合软孔隙在零压力时的纵横比最小值)的解析表达式,并在此基础上利用岩石速度-压力实验观测结果给出求取介质中两类孔隙纵横比及其含量分布特征的方法.通过逐步迭代加入软孔隙的方法对基于特征纵横比的"喷射流"(squirt fluid)模型进行了扩展,以考虑复杂孔隙分布特征对岩石喷射流作用的影响及其可能引起的速度频散特征.相较于典型的喷射流作用速度频散模式,对于岩石中软孔隙纵横比及其对应含量在较宽的范围呈谱分布的一般情况,其速度频散曲线不存在明显的低频段和中间频段,速度随频率的增大呈递增趋势直至高频极限.这说明即使在地震频段,微观尺度下的喷射流作用仍起一定作用,同样会造成流体饱和岩石介质的地震速度与Gassmann方程预测结果有不可忽略的差异.本文是对现有喷射流模型的重要补充,也为利用实验数据建立不同频段间岩石弹性波传播速度的可能联系提供了理论依据. 相似文献
5.
地下多孔介质中的孔隙类型复杂多样,既有硬孔又有扁平的软孔.针对复杂孔隙介质,假设多孔介质中同时含有球型硬孔和两种不同产状的裂隙(硬币型、尖灭型裂隙),当孔隙介质承载载荷时,考虑两种不同类型的裂隙对于孔隙流体压力的影响,建立起Biot理论框架下饱和流体情况含混合裂隙、孔隙介质的弹性波动方程,并进一步求取了饱和流体情况下仅由裂隙引起流体流动时的含混合裂隙、孔隙介质的体积模量和剪切模量,随后,在此基础上讨论了含混合裂隙、孔隙介质在封闭条件下地震波衰减和频散的高低频极限表达式.最后计算了给定模型的地震波衰减和频散,发现地震波衰减曲线呈现"多峰"现象,速度曲线为"多频段"频散.针对该模型分析讨论了渗透率参数、裂隙纵横比参数以及流体黏滞性参数对于地震波衰减和频散的影响,表明三个参数均为频率控制参数. 相似文献
6.
通过融合Gassmann方程和由微分等效介质理论建立的干岩石骨架模型--DEM解析模型,本文提出根据纵波(和横波)速度反演岩石等效孔隙纵横比进行储层孔隙结构评价和横波速度预测的方法.首先,利用Gassmann方程和DEM解析模型建立岩石的纵、横波速度与密度、孔隙度、饱和度和矿物组分等各参数之间的关系;其次,将岩石孔隙等效为具有单一纵横比的理想椭球孔,应用非线性全局寻优算法来寻找最佳的等效孔隙纵横比使得理论预测与实际测量的弹性模量之间的误差最小;最后,将反演得到的等效孔隙纵横比代入到Gassmann方程和DEM解析模型中构建横波速度.实验室和井孔测量数据应用表明,反演得到的等效孔隙纵横比可准确反映储层的孔隙结构,对于裂缝型储层如花岗岩,其孔隙纵横比通常小于0.025,而对于孔隙型储层如砂岩,其孔隙纵横比通常大于0.08.只利用纵波与同时利用纵、横波反演得到的孔隙纵横比结果几乎完全一致,而且由纵波构建的横波与实测横波吻合良好,说明本文提出的等效孔隙纵横比反演及其横波速度预测方法是有效的. 相似文献
7.
碳酸盐岩组成复杂,储集层空间类型多,储层参数的规律性较差,AVO响应特征弱,尤其是孔隙形状对弹性参数的影响不可忽略.针对碳酸盐岩储层的特点,基于自洽模型和微分等效模型构建了多孔隙结构条件下多相介质的碳酸盐岩岩石物理模型.按照各相介质体积含量由小到大的顺序,先利用自洽模型计算孔隙度为50%时的双相等效介质,再通过微分等效模型将孔隙度调整至实际孔隙度,计算实际孔隙条件下的双相介质等效体,重复迭代该过程,最终构建出与实际岩石等效的介质体,计算获得等效岩石的弹性模量参数.通过实验对比了自洽模型、微分等效模型与自洽-微分等效模型的岩石弹性模量与孔隙度、孔隙纵横比之间的适用性规律.最后基于四川龙岗飞仙关组实际测井资料,结合自洽-微分等效模型进行了目标地层的横波预测.结果显示该模型的横波预测效果相对Xu-Payne模型预测,两者预测曲线比较接近,但在形态上更接近实际横波.SCA-DEM模型在孔隙度和孔隙纵横比上具有更好的适应性,为碳酸盐岩孔洞型储层的岩石物理建模提供了一种新参考. 相似文献
8.
孔隙纵横比是描述多孔岩石微观孔隙结构特征的重要参数,目前用于获取岩石完整孔隙纵横比分布的经典模型为David-Zimmerman(D-Z)孔隙结构模型,该模型假设岩石由固体矿物基质、一组纵横比相等的硬孔隙以及多组纵横比不等的微裂隙构成,并认为固体矿物基质和硬孔隙均不受压力影响,在此基础上,利用超声纵横波速度的压力依赖性反演岩石硬孔隙和各组微裂隙的孔隙纵横比及孔隙度.该方法的关键点在于以累积裂隙密度为桥梁,借助等效介质理论建立了岩石弹性模量和孔隙纵横比之间的内在联系.但在D-Z模型中,多重孔隙岩石累积裂隙密度的计算直接由单重孔隙裂隙密度公式实现,这种近似导致该模型在许多情况下难以获得良好的反演精度.为了完善经典D-Z模型,本文提出了一种基于虚拟降压的孔隙纵横比分布反演策略,通过多个假想降压过程实现累积裂隙密度的准确计算,并将基于DEM和MT的经典D-Z模型推广到KT和SCA中,结合四种等效介质理论建立了一套完整的反演流程.采用一系列砂岩和碳酸盐岩样品,测试了反演流程在实际岩芯孔隙纵横比提取中的应用效果,研究结果表明:与D-Z模型相比,本文方法的模拟结果与实际数据吻合更好,并同时适用于砂岩和碳酸盐岩;此外,通过分析四种等效介质理论的模拟结果发现,本文方法并不十分依赖于等效介质理论的选择,这些理论获得的孔隙结构参数随压力的变化趋势基本一致,数值上仅存在略微差异,且这种差异随着压力的增大逐渐消失.本文方法是经典D-Z孔隙结构模型的重要补充,对岩石孔隙结构表征、流体饱和岩石速度预测以及孔间喷射流效应的模拟具有十分重要的意义. 相似文献
9.
复杂孔隙储层往往同时发育孔缝洞等多种孔隙类型,这种孔隙结构的复杂性使得岩石的速度与孔隙度之间的相关性很差.经典的二维岩石物理模版只研究弹性参数与孔隙度和饱和度之间的定量关系,而不考虑孔隙结构的影响,用这样的模版来预测复杂孔隙储层的物性参数时带来很大偏差.本文首先证明多重孔隙岩石的干骨架弹性参数可以用一个等效孔隙纵横比的单重孔隙岩石物理模型来模拟;进而基于等效介质岩石物理理论和Gassmann方程,建立一个全新的三维岩石物理模版,用它来建立复杂孔隙岩石的弹性性质与孔隙扁度及孔隙度和饱和度之间的定量关系;在此基础上,预测复杂储层的孔隙扁度、孔隙度以及孔隙中所包含的流体饱和度.实际测井和地震反演数据试验表明,三维岩石物理模版可有效提高复杂孔隙储层参数的预测精度. 相似文献
10.
多形态拓展孔隙、裂隙并存理论描述了介质中孔隙与多形态裂隙体之间的相互作用,现已被用于岩石的孔隙结构反演．利用实验室加压条件下干燥和流体饱和岩石的纵、横波速测量数据,联合反演了潜山储层岩石的孔隙纵横比谱．反演结果表明：随着有效压力的增加,裂隙元素的纵横比减小,小纵横比裂隙的闭合降低了裂隙密度,导致低压段的弹性波速度显著增大．同时,在双对数刻度下,孔隙度与孔隙纵横比之间存在幂函数关系,孔隙结构指数随总孔隙度单调增加,此关系可用于判别潜山储层的有效性．本文的研究不仅为潜山储层岩石孔隙结构分析方法提供了实验验证,而且给出了一种潜山储层有效性的辅助判别方法． 相似文献
11.
R. Agersborg M. Jakobsen B. O. Ruud T. A. Johansen 《Studia Geophysica et Geodaetica》2007,51(1):89-118
An effective medium model for the stress-dependent seismic properties of fractured reservoirs is developed here on the basis
of a combination of a general theory of viscoelastic waves in rock-like composites with recently published formulae for deformation
of communicating and interacting cavities (interconnected pores/cracks and fractures at finite concentration) under drained
loading. The inclusion-based model operates with spheroidal cavities at two different length scales; namely, the microscopic
scale of the pores and (grain-boundary) cracks, and the mesoscopic scale of the fractures (controlling the flow of fluid).
The different cavity types can in principle have any orientation and aspect ratio, but the microscopic pores/cracks and mesoscopic
fractures were here assumed to be randomly and vertically oriented, respectively. By using three different aspect ratios for
the relatively round pores (representing the stiff part of the pore space) and a distribution of aspect ratios for the relatively
flat cracks (representing the compliant part of the pore space), we obtained a good fit between theoretical predictions and
ultrasonic laboratory measurements on an unfractured rock sample under dry conditions. By using a single aspect ratio for
the mesoscopic fractures, we arrived at a higher-order microstructural model of fractured porous media which represents a
generalization of the first-order model developed by Chapman et al. (2002,2003). The effect of cavity size was here modelled
under the assumption that the characteristic time for wave-induced (squirt) flow at the scale of a particular cavity (pore/crack
vs. fracture) is proportional with the relevant scale-size. In the modelling, we investigate the effect of a decreasing pore
pressure with constant confining pressure (fixed depth), and hence, increasing effective pressure. The analysis shows that
the attenuation-peak due to the mesoscopic fractures in the reservoir will move downward in frequency as the effective pressure
increases. In the range of seismic frequencies, our modelling indicates that the P-wave velocities may change by more than
20% perpendicular to the fractures and close to 10% parallel to the fractures. In comparison, the vertical S-wave velocities
change by only about 5% for both polarization directions (perpendicular and parallel to the fractures) when the effective
pressure increases from 0 to 15 MPa. This change is mainly due to the overall change in porosity with pressure. The weak pressure
dependence is a consequence of the fact that the S waves will only sense if the fractures are open or not, and since all the
fractures have the same aspect ratio, they will close at the same effective pressure (which is outside the analysed interval).
Approximate reflection coefficients were computed for a model consisting of the fractured reservoir embedded as a layer in
an isotropic shale and analysed with respect to variations in Amplitude Versus Offset and aZimuth (AVOZ) properties at seismic
frequencies for increasing effective pressure. For the P-P reflections at the top of the reservoir, it is found that there
is a significant dependence on effective pressure, but that the variations with azimuth and offset are small. The lack of
azimuthal dependence may be explained from the approximate reflection coefficient formula as a result of cancellation of terms
related to the S-wave velocity and the Thomson’s anisotropy parameter δ. For the P-S reflection, the azimuthal dependence
is larger, but the pressure dependence is weaker (due to a single aspect ratio for the fractures).
Finally, using the effective stiffness tensor for the fractured reservoir model with a visco-elastic finite-difference code,
synthetic seismograms and hodograms were computed. From the seismograms, attenuation changes in the P wave reflected at the
bottom of the reservoir can be observed as the effective pressure increases. S waves are not much affected by the fractures
with respect to attenuation, but azimuthal dependence is stronger than for P waves, and S-wave splitting in the bottom reservoir
P-S reflection is clearly seen both in the seismograms and hodograms. From the hodograms, some variation in the P-S reflection
with effective pressure can also be observed. 相似文献
12.
Seismic wave propagation in reservoir rocks is often strongly affected by fractures and micropores. Elastic properties of fractured reservoirs are studied using a fractured porous rock model, in which fractures are considered to be embedded in a homogeneous porous background. The paper presents an equivalent media model for fractured porous rocks. Fractures are described in a stress‐strain relationship in terms of fracture‐induced anisotropy. The equations of poroelasticity are used to describe the background porous matrix and the contents of the fractures are inserted into a matrix. Based on the fractured equivalent‐medium theory and Biot's equations of poroelasticity, two sets of porosity are considered in a constitutive equation. The porous matrix permeability and fracture permeability are analysed by using the continuum media seepage theory in equations of motion. We then design a fractured porous equivalent medium and derive the modified effective constants for low‐frequency elastic constants due to the presence of fractures. The expressions of elastic constants are concise and are directly related to the properties of the main porous matrix, the inserted fractures and the pore fluid. The phase velocity and attenuation of the fractured porous equivalent media are investigated based on this model. Numerical simulations are performed. We show that the fractures and pores strongly influence wave propagation, induce anisotropy and cause poroelastic behaviour in the wavefields. We observe that the presence of fractures gives rise to changes in phase velocity and attenuation, especially for the slow P‐wave in the direction parallel to the fracture plane. 相似文献
13.
Todd H. Skaggs 《Advances in water resources》2011,34(10):1335-1342
Critical path analysis (CPA) is a method for estimating macroscopic transport coefficients of heterogeneous materials that are highly disordered at the micro-scale. Developed originally to model conduction in semiconductors, numerous researchers have noted that CPA might also have relevance to flow and transport processes in porous media. However, the results of several numerical investigations of critical path analysis on pore network models raise questions about the applicability of CPA to porous media. Among other things, these studies found that (i) in well-connected 3D networks, CPA predictions were inaccurate and became worse when heterogeneity was increased; and (ii) CPA could not fully explain the transport properties of 2D networks. To better understand the applicability of CPA to porous media, we made numerical computations of permeability and electrical conductivity on 2D and 3D networks with differing pore-size distributions and geometries. A new CPA model for the relationship between the permeability and electrical conductivity was found to be in good agreement with numerical data, and to be a significant improvement over a classical CPA model. In sufficiently disordered 3D networks, the new CPA prediction was within ±20% of the true value, and was nearly optimal in terms of minimizing the squared prediction errors across differing network configurations. The agreement of CPA predictions with 2D network computations was similarly good, although 2D networks are in general not well-suited for evaluating CPA. Numerical transport coefficients derived for regular 3D networks of slit-shaped pores were found to be in better agreement with experimental data from rock samples than were coefficients derived for networks of cylindrical pores. 相似文献
14.
地震波传播激发的不同尺度的流固相对运动(宏观、中观和微观)是许多沉积岩地层中地震波频散和衰减的主要原因,然而野外观测和试验测量都难以对非均匀多孔介质孔隙压力弛豫物理过程进行精细刻画.通过数字岩石物理技术,本文建立了三个典型的数字岩心分别用于表征孔隙结构、岩石骨架和斑状饱和流体引起的非均质性,利用动态应力应变模拟技术计算数字岩心的位移和孔隙流体增量图像.通过分析和比较三个数字岩心的位移和孔隙压力增量图像,细致刻画了发生于非均匀含流体多孔介质内的宏观、中观和微观尺度的流固相对运动:1)宏观尺度的波致孔隙流体流动导致波长尺度上数字岩心不同区域的孔隙压力和位移差异;2)中观尺度的流体流动发生在软层与硬层之间、气层与液层之间;3)微观尺度的流体流动发生在孔隙内部或相邻孔隙之间.数值模拟试验也证明基于数字岩心的动态应力应变模拟技术可以从微观尺度上更好的理解波致孔隙流体流动发生的物理机理,从而为建立岩石骨架、孔隙流体、孔隙结构非均质性和弹性波频散-衰减特征的映射关系奠定基础. 相似文献
15.
Wave‐induced fluid flow plays an important role in affecting the seismic dispersion and attenuation of fractured porous rocks. While numerous theoretical models have been proposed for the seismic dispersion and attenuation in fractured porous rocks, most of them neglect the wave‐induced fluid flow resulting from the background anisotropy (e.g. the interlayer fluid flow between different layers) that can be normal in real reservoirs. Here, according to the theories of poroelasticity, we present an approach to study the frequency‐dependent seismic properties of more realistic and complicated rocks, i.e. horizontally and periodically layered porous rock with horizontal and randomly orienting fractures, respectively, distributed in one of the two periodical layers. The approach accounts for the dual effects of the wave‐induced fluid flow between the fractures and the background pores and between different layers (the interlayer fluid flow). Because C33 (i.e., the modulus of the normally incident P‐wave) is directly related to the P‐wave velocity widely measured in the seismic exploration, and its comprehensive dispersion and attenuation are found to be most significant, we study mainly the effects of fracture properties and the stiffness contrast between the different layers on the seismic dispersion and attenuation of C33. The results show that the increasing stiffness contrast enhances the interlayer fluid flow of the layered porous rocks with both horizontal and randomly orienting fractures and weakens the wave‐induced fluid flow between the fractures and the background pores, especially for the layered porous rock with horizontal fractures. The modelling results also demonstrate that for the considered rock construction, the increasing fracture density reduces the interlayer fluid flow while improves the dispersion and attenuation in the fracture‐relevant frequency band. Increasing fracture aspect ratio is found to reduce the dispersion and attenuation in the fracture‐relevant frequency band only, especially for the layered porous rock with horizontal fractures. 相似文献
16.
利用核磁共振技术对致密砂岩储层不同渗透率级别基质岩心和裂缝基质岩心不同驱替压力下CO2驱油特征进行了研究,简述核磁共振原理及实验方法。表明:致密砂岩储层特低、超低渗透基质岩心在初始CO2驱替压力下,岩心毛细孔隙和微毛细孔隙区间的油不同程度被采出,随着CO2驱替压力增大,特低、超低渗透基质岩心毛细孔隙区间油的采出程度不断增加且累积采出程度不同。裂缝致密砂岩储层岩心,裂缝和毛细孔隙区间的油在初始CO2驱替压力下,岩心裂缝中的油及毛细孔隙中的部分油被驱替出来,CO2驱替压力提高毛细孔隙、微毛细区间油的采出程度和累积采收程度较小。致密砂岩储层特低、超低渗透基质岩心和裂缝致密砂岩储层岩心,随着CO2驱替压力增大毛细孔隙区间的部分剩余油成正比增加进入到微毛细孔隙区间改变储层剩余油分布。核磁共振技术能够深入研究致密砂岩储层CO2不同驱替压力阶段,岩心裂缝、毛细孔隙区间、微毛细孔隙区间油的采出程度和剩余油分布情况,对于研究致密砂岩储层微观驱油机理具有较重要的价值。 相似文献