压力及流体对岩石孔隙结构与黏弹性的影响规律研究:理论模型及实验观测

地质成因和构造/热应力导致地壳岩石中的孔隙结构（裂隙和粒间孔）的变化.影响岩石黏弹性的因素包括压力、孔隙度、孔隙中包含的流体和孔隙几何形状等.相对于岩石中的硬孔隙,岩石黏弹性（衰减和频散）受软孔隙（裂隙）的影响更大.本文选取三块白云岩样本,进行了不同围压和流体条件下的超声波实验测量.利用CPEM（Cracks and Pores Effective Medium,裂隙和孔隙有效介质）模型获得了岩石高、低频极限的弹性模量,并通过Zener体（标准线性体）模型将CPEM模型拓展到全频带而得到CPEM-Zener模型,用该模型拟合岩石松弛和非松弛状态下的实验数据,本文得到平均裂隙纵横比和裂隙孔隙度以及纵波速度和品质因子随频率的变化关系.结果表明,饱水岩石的平均裂隙纵横比和裂隙孔隙度均高于饱油岩石,随着压差（围压和孔隙压力的差值）的增加,饱油岩石中的裂隙首先闭合.并且压差在70 MPa以内时,随着压差增大,岩石的平均裂隙纵横比和裂隙孔隙度在饱水和饱油时的差值增大,此时流体类型对于岩石裂隙的影响越来越显著,此外,对饱水岩石,平均裂隙纵横比随压差增加而增大,这可能是由于岩石中纵横比较小的裂隙会随压差增大而逐渐趋于闭合.在饱水和饱油岩石中,裂隙孔隙度和裂隙密度都随着压差增加而减小.通过对裂隙密度和压差的关系进行指数拟合,本文获得压差趋于0时的裂隙密度,且裂隙密度随孔隙度增大而增大,增大速率随压差增加而降低.针对饱水和饱油的白云岩样本,CPEM-Zener模型预测的纵波频散随压差增大而减小,此变化趋势和实验测得的逆品质因子随压差的变化关系基本一致,由此进一步验证了模型的实用性.本研究对岩石的孔隙结构和黏弹性分析以及声波测井、地震勘探的现场应用有指导意义.     

压力及流体对岩石孔隙结构与黏弹性的影响规律研究:理论模型及实验观测

地质成因和构造/热应力导致地壳岩石中的孔隙结构（裂隙和粒间孔）的变化.影响岩石黏弹性的因素包括压力、孔隙度、孔隙中包含的流体和孔隙几何形状等.相对于岩石中的硬孔隙,岩石黏弹性（衰减和频散）受软孔隙（裂隙）的影响更大.本文选取三块白云岩样本,进行了不同围压和流体条件下的超声波实验测量.利用CPEM（Cracks and Pores Effective Medium,裂隙和孔隙有效介质）模型获得了岩石高、低频极限的弹性模量,并通过Zener体（标准线性体）模型将CPEM模型拓展到全频带而得到CPEM-Zener模型,用该模型拟合岩石松弛和非松弛状态下的实验数据,本文得到平均裂隙纵横比和裂隙孔隙度以及纵波速度和品质因子随频率的变化关系.结果表明,饱水岩石的平均裂隙纵横比和裂隙孔隙度均高于饱油岩石,随着压差（围压和孔隙压力的差值）的增加,饱油岩石中的裂隙首先闭合.并且压差在70 MPa以内时,随着压差增大,岩石的平均裂隙纵横比和裂隙孔隙度在饱水和饱油时的差值增大,此时流体类型对于岩石裂隙的影响越来越显著,此外,对饱水岩石,平均裂隙纵横比随压差增加而增大,这可能是由于岩石中纵横比较小的裂隙会随压差增大而逐渐趋于闭合.在饱水和饱油岩石中,裂隙孔隙度和裂隙密度都随着压差增加而减小.通过对裂隙密度和压差的关系进行指数拟合,本文获得压差趋于0时的裂隙密度,且裂隙密度随孔隙度增大而增大,增大速率随压差增加而降低.针对饱水和饱油的白云岩样本,CPEM-Zener模型预测的纵波频散随压差增大而减小,此变化趋势和实验测得的逆品质因子随压差的变化关系基本一致,由此进一步验证了模型的实用性.本研究对岩石的孔隙结构和黏弹性分析以及声波测井、地震勘探的现场应用有指导意义.

     

水饱和裂纹对地壳岩样中地震波速及各向异性的影响

选择４种地壳岩石样品,经干燥或水饱和处理后在不同围压条件下测量了在其中传播的纵、横波的速度及其各向异性．在大气压条件下低孔隙度(＜１％岩样中,水饱和样品中的纵波速度明显地比干燥样品中的高,但横波速度的差别不大．因为在低孔隙度岩样中纵波速度对孔隙流体的反应比横波速度敏感,可以用泊松比的变化来反映随着围压的增加晶粒间流体对弹性波传播特性的影响．根据实验数据,按Ｏ’Ｃｏｎｎｅｌｌ模型分别计算了干燥和水饱和岩样中的裂纹密度,与通过实测体应变曲线得到的裂纹孔隙度十分吻合．利用横波的速度和偏振特性可以推断岩样中定向排列微裂纹的空间取向情况．研究表明,同时测量在岩样中传播的纵、横波的速度,通过Ｖｐ／Ｖｓ比值可以给出有关颗粒边界流体的证据,也可以估计岩样中的裂纹密度．     

水饱和裂纹对地壳岩样中地震波速及各向异性的影响

选择４种地壳岩石样品,经干燥或水饱和处理后在不同围压条件下测量了在其中传播的纵、横波的速度及其各向异性．在大气压条件下低孔隙度(＜１％岩样中,水饱和样品中的纵波速度明显地比干燥样品中的高,但横波速度的差别不大．因为在低孔隙度岩样中纵波速度对孔隙流体的反应比横波速度敏感,可以用泊松比的变化来反映随着围压的增加晶粒间流体对弹性波传播特性的影响．根据实验数据,按Ｏ’Ｃｏｎｎｅｌｌ模型分别计算了干燥和水饱和岩样中的裂纹密度,与通过实测体应变曲线得到的裂纹孔隙度十分吻合．利用横波的速度和偏振特性可以推断岩样中定向排列微裂纹的空间取向情况．研究表明,同时测量在岩样中传播的纵、横波的速度,通过Ｖｐ／Ｖｓ比值可以给出有关颗粒边界流体的证据,也可以估计岩样中的裂纹密度．     

The effect of stress on pore pressure in rocks and the mechanism of water table anomaly before earthquakes

Triaxial compressive experiments of porous rock samples were carried out under various confining pressures and initial pore pressures without drainage; axial strain and pore pressure were observed versus differential stress. The results of such experiments show that pore pressure increases with increase of differential stress at low differential stress; pore pressure decreases with increase of differential stress at medium and high differential stress. Pore pressure also increases with large amplitude decrease of differential stress at high differential stress. Based on such experiments, it is suggested that water table anomaly before an earthquake reflects the change of differential stress in crustal rocks. The anomalous behavior of water tables in the epicentral and peripheral areas before the great Tangshan earthquake of July 28, 1976 are explained by such suggestion. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 88–95, 1991. This study is supported by the Chinese Joint Seismological Science Foundation. Professor Yongtai Che gave us much help in applying fund support and supplying earthquake case histories.     

应力对岩石中孔隙流体压力的影响和地下水位震前异常变化机理.

在不同的围压和初始孔隙压力作用下,对孔隙岩石试件进行了不排液的三轴压缩试验,观测了轴向应变和孔隙压力随差应力的变化.试验结果表明:在低差应力下,孔隙压力随差应力的增加而增加.在中等和高差应力下孔隙压力随差应力的增加而降低.在高差应力下,孔隙压力也会随差应力的大幅度降低而升高.基于上述试验,假定地下水位的震前异常变化是地壳岩石中差应力变化的反映.根据这种假定解释了唐山大震前震中及外围地区地下水位异常变化的特征.      

平面P波在饱和半空间中洞室周围的散射(II):数值结果

本文通过数值计算研究了入射平面P波在饱和半空间中洞室周围散射问题,分析了入射波频率和角度、边界渗透条件、孔隙率、泊松比等参数对散射的影响。研究表明,平面P波入射情况下,饱和半空间和单相(干土)半空间中洞室附近地表位移幅值的差别很大,干土情况的水平位移幅值相对较大,饱和情况的竖向地表位移幅值相对较大;由于波在洞室附近的干涉,饱和情况与干土情况的地表位移出现相位漂移。随着孔隙率的增大,洞室附近水平地表位移幅值逐渐减小,竖向地表位移幅值则逐渐增大;当孔隙率较低时,边界渗透条件对地表位移幅值的影响很小,而当孔隙率较大时,边界渗透条件的影响则不可忽视,不透水情况下,水平和竖向地表位移幅值的峰值均相对较大;随着入射频率的增加,孔隙率的影响逐渐增大,而且不透水情况下孔隙率的影响相对较大。随着泊松比的增大,洞室附近水平地表位移幅值逐渐降低,竖向地表位移幅值则逐渐增大;泊松比较小时,边界渗透条件对位移幅值的影响较大,泊松比较大时,边界渗透条件对位移幅值的影响则较小;随着入射频率的增加,泊松比的影响逐渐增大。当孔隙率较小时,半空间地表和洞室表面孔隙水压幅值较小,但空间变化比较剧烈,随着孔隙率的增大,孔隙水压逐渐增大但空间变化逐渐平缓;随着入射频率的增加,孔隙水压幅值逐渐增大,且孔隙水压的空间变化逐渐变得复杂。     

碳酸盐岩跨频段岩石物理测量与理论建模——不同孔隙结构对碳酸盐岩频散与衰减的影响研究

碳酸盐岩孔隙结构类型复杂多样,当地震波经过含有不同孔隙结构的流体饱和岩石后往往会产生不同的波频散和衰减特征,这使得根据波的不同响应特征来推断碳酸盐岩的孔隙结构类型,甚至孔隙流体性质信息成为可能.本文针对白云岩、灰岩以及人工碳酸盐岩样品开展了跨频段（超声+低频）实验测量和理论建模,探索碳酸盐岩的孔隙结构类型和孔隙流体对模量频散和衰减的影响机制.首先根据铸体薄片、扫描电镜的图像对碳酸盐岩样品进行了孔隙结构类型分析,并将样品主要分为裂缝型、裂缝-孔隙型、孔洞型三类,然后测量了相应样品完全饱和流体后在不同围压下的模量频散与衰减.在完全饱和甘油并处于低围压时,裂缝型与孔洞型样品均出现一个衰减峰,分别位于1 Hz与100 Hz附近,而裂缝-孔隙型样品则具有两个衰减峰,一个在1 Hz附近,另一个在100 Hz附近.裂缝型样品（裂缝主导）的衰减峰相比孔洞型样品（中等刚度孔隙主导）对应的衰减峰在低围压下幅度更大,且对围压变化更敏感.在测量数据的基础上,建立了考虑纵横比分布的软孔隙和中等刚度孔隙的喷射流模型,认为该模型能一定程度上解释裂缝型、裂缝-孔隙型、孔洞型三种类型碳酸盐岩在测量频带的频散.以上研究加深了对不同孔隙类型主导的碳酸盐岩储层地震响应特征的认识,对储层预测工作的进一步精细化具有重要意义.     

碳酸盐岩跨频段岩石物理测量与理论建模——不同孔隙结构对碳酸盐岩频散与衰减的影响研究

碳酸盐岩孔隙结构类型复杂多样,当地震波经过含有不同孔隙结构的流体饱和岩石后往往会产生不同的波频散和衰减特征,这使得根据波的不同响应特征来推断碳酸盐岩的孔隙结构类型,甚至孔隙流体性质信息成为可能.本文针对白云岩、灰岩以及人工碳酸盐岩样品开展了跨频段（超声+低频）实验测量和理论建模,探索碳酸盐岩的孔隙结构类型和孔隙流体对模量频散和衰减的影响机制.首先根据铸体薄片、扫描电镜的图像对碳酸盐岩样品进行了孔隙结构类型分析,并将样品主要分为裂缝型、裂缝-孔隙型、孔洞型三类,然后测量了相应样品完全饱和流体后在不同围压下的模量频散与衰减.在完全饱和甘油并处于低围压时,裂缝型与孔洞型样品均出现一个衰减峰,分别位于1 Hz与100 Hz附近,而裂缝-孔隙型样品则具有两个衰减峰,一个在1 Hz附近,另一个在100 Hz附近.裂缝型样品（裂缝主导）的衰减峰相比孔洞型样品（中等刚度孔隙主导）对应的衰减峰在低围压下幅度更大,且对围压变化更敏感.在测量数据的基础上,建立了考虑纵横比分布的软孔隙和中等刚度孔隙的喷射流模型,认为该模型能一定程度上解释裂缝型、裂缝-孔隙型、孔洞型三种类型碳酸盐岩在测量频带的频散.以上研究加深了对不同孔隙类型主导的碳酸盐岩储层地震响应特征的认识,对储层预测工作的进一步精细化具有重要意义.

     

Velocity anisotropy and attenuation of shale in under- and overpressured conditions

The seismic velocity and attenuation of fully saturated shales were measured for the first time under overpressured conditions, using the ultrasonic reflection technique. Shale cores from naturally overpressured horizons in the North Sea were tested in the laboratory, at confining and pore pressures relevant to in situ conditions.
A single-frequency tone-burst pulse wave was used to determine the seismic wave velocities and quality factors of the shale samples, with errors less than 0.3% and 0.1 dB/cm, respectively, at a frequency of 0.75 MHz. Sample length changes with varying confining and pore pressure were measured and the pore pressure equilibration time was monitored for each sample.
The anisotropy of the seismic attributes ( V _p, V _s, Q _p and Q _s) was determined over a range of differential pressures from 5 to 60 MPa, with respect to the predominant foliation. The ultrasonic velocity data followed a transversely isotropic pattern depending on the direction of wave propagation with respect to the laminations. The Poisson's ratio was found to rise by 5% as the shale material progressed from a normally pressured to an overpressured state. The quality factor ( Q ) characteristics were interpreted in terms of pore geometry and connectivity as well as the directional permeability of the transversely isotropic shale material. The results were converted to bulk and shear loss modulus defects, and a positive bulk loss was observed for waves propagating perpendicular to the lamination plane even above differential pressures of 20 MPa. This indicates different levels of Biot-flow and squirt-flow attenuation mechanisms acting within the shale structure, depending on the wave propagation and vibration directions.     

Diffraction of plane SV waves by a shallow circular-arc canyon in a saturated poroelastic half-space

An analytical solution for the scattering and diffraction of incident plane SV waves by a shallow circular-arc canyon in a saturated poroelastic half-space is derived by the wave function expansion method. The solution is utilized to analyze the dependence of the computed surface motions on the incident frequencies, incident angles, porosity, boundary drainage and Poisson's ratio. It is shown that, depending on the incident angles, the surface displacement amplitudes around a canyon in a dry poroelastic half-space and saturated poroelastic half-space can be very different. The surface displacement amplitudes of an undrained saturated poroelastic half-space are close to those of a drained saturated poroelastic half-space. For low porosity, the surface displacement amplitudes of a saturated poroelastic half-space are almost identical to those of a dry poroelastic half-space, and drainage condition has little influence on the surface displacement amplitudes. But for high porosity, the effect of drainage condition becomes significant, and for the same porosity, the displacement amplitudes of an undrained saturated half-space will be larger than those of a drained saturated half-space. Poisson's ratio is also an important factor affecting the surface displacement amplitudes around the canyon, both in drained and undrained conditions, but leads to larger effects for an undrained saturated half-space than for a drained saturated half-space. Large pore pressures are found around the canyon and their amplitudes depend on the incident angles and frequencies. Below the surface, the amplitudes of pore pressures are less than they are at the surface, especially for high frequencies.     

Compaction and Failure in High Porosity Carbonates: Mechanical Data and Microstructural Observations

We investigated systematically the micromechanics of compaction in two carbonates of porosity above 30%, Majella grainstone and Saint Maximin limestone. The composition, grain size and pore surface area of these rocks were determined. Hydrostatic compression experiments were performed under dry and wet conditions beyond the onset of grain crushing. A significant water weakening effect was observed in both rocks. A set of conventional triaxial experiments was also performed on both rocks under dry conditions at confining pressures ranging from 3 to 31 MPa. Microstructural observations were carried out on the deformed samples. The mechanical behavior of these high porosity carbonates is dominated by shear-enhanced compaction associated in most cases with strain hardening. Stress-induced cracking and grain crushing are the dominant micromechanisms of deformation in both rocks. In Majella grainstone, compactive shear bands appeared at low confinement, in qualitative agreement with the deformation bands observed in the field. At higher confining pressures, compaction localization was inhibited and homogeneous cataclastic flow developed. In Saint Maximin limestone, compaction localization was observed at all confining pressures. An increasing number of compactive shear bands at various orientations appeared with increasing strain. These new data suggest that compaction localization is important in the mechanical compaction of high porosity carbonates.     

不同围压下岩石中泊松比的各向异性

由弹性波纵横波速比计算得到的泊松比,在利用地震波反演地下结构和物质组成时,可以提供比单纯地利用纵波或者横波波速更强的约束. 为了研究不同岩石中泊松比及其各向异性随围压的变化,选取19块典型的岩石样品,在不同围压下沿相互正交的3个方向同时测量纵、横波速度,并通过计算得到了泊松比及其各向异性. 结果表明,在大部分样品中泊松比各向异性与纵波速度各向异性有一定的相关性,但在有的样品中泊松比各向异性与纵波速度各向异性没有明显的关系,甚至表现出完全相反的变化规律. 由于泊松比也表现出很强的各向异性,在排除微裂隙的影响后有的仍可以达到30髎以上,在利用实际地震波传播数据通过泊松比反演地下物质结构时必须考虑其各向异性的影响.     

不同围压下岩石中泊松比的各向异性

由弹性波纵横波速比计算得到的泊松比,在利用地震波反演地下结构和物质组成时,可以提供比单纯地利用纵波或者横波波速更强的约束. 为了研究不同岩石中泊松比及其各向异性随围压的变化,选取19块典型的岩石样品,在不同围压下沿相互正交的3个方向同时测量纵、横波速度,并通过计算得到了泊松比及其各向异性. 结果表明,在大部分样品中泊松比各向异性与纵波速度各向异性有一定的相关性,但在有的样品中泊松比各向异性与纵波速度各向异性没有明显的关系,甚至表现出完全相反的变化规律. 由于泊松比也表现出很强的各向异性,在排除微裂隙的影响后有的仍可以达到30髎以上,在利用实际地震波传播数据通过泊松比反演地下物质结构时必须考虑其各向异性的影响.     

含油岩石弹性特征及其与油气的关系

本文用准噶尔盆地东部油区58口井共2257块油气储层岩芯标本,在常温常压和模拟地下温度压力下对V_P/V_S及衰减系数进行了精确测量,结合密度、孔隙度及岩石学研究获得个同岩性地层含油、水、气不同流体的V_P/V_S比、泊松比及不同深度的速度校正曲线;用回归分析获得V_P-V_S,密度-速度,密度-孔隙度的相关关系式;用数学统计获得稠油地区岩芯砂岩的速度特征值;用AVO技术进行模型计算和实例验证.为油气碳氢检测,测井速度标定,声阻抗及演,模型正演等提供了有价值的数据.     

含油岩石弹性特征及其与油气的关系

本文用准噶尔盆地东部油区58口井共2257块油气储层岩芯标本,在常温常压和模拟地下温度压力下对V_P/V_S及衰减系数进行了精确测量,结合密度、孔隙度及岩石学研究获得个同岩性地层含油、水、气不同流体的V_P/V_S比、泊松比及不同深度的速度校正曲线;用回归分析获得V_P-V_S,密度-速度,密度-孔隙度的相关关系式;用数学统计获得稠油地区岩芯砂岩的速度特征值;用AVO技术进行模型计算和实例验证.为油气碳氢检测,测井速度标定,声阻抗及演,模型正演等提供了有价值的数据.     

Acoustic and petrophysical relationships in low-shale sandstone reservoir rocks

This paper describes the measurements of the acoustic and petrophysical properties of two suites of low‐shale sandstone samples from North Sea hydrocarbon reservoirs, under simulated reservoir conditions. The acoustic velocities and quality factors of the samples, saturated with different pore fluids (brine, dead oil and kerosene), were measured at a frequency of about 0.8 MHz and over a range of pressures from 5 MPa to 40 MPa. The compressional‐wave velocity is strongly correlated with the shear‐wave velocity in this suite of rocks. The ratio V_P/V_S varies significantly with change of both pore‐fluid type and differential pressure, confirming the usefulness of this parameter for seismic monitoring of producing reservoirs. The results of quality factor measurements were compared with predictions from Biot‐flow and squirt‐flow loss mechanisms. The results suggested that the dominating loss in these samples is due to squirt‐flow of fluid between the pores of various geometries. The contribution of the Biot‐flow loss mechanism to the total loss is negligible. The compressional‐wave quality factor was shown to be inversely correlated with rock permeability, suggesting the possibility of using attenuation as a permeability indicator tool in low‐shale, high‐porosity sandstone reservoirs.     

Influence of pore pressure on velocity in low-porosity sandstone: Implications for time-lapse feasibility and pore-pressure study

As seismic data quality improves, time‐lapse seismic data is increasingly being called upon to interpret and predict changes during reservoir development and production. Since pressure change is a major component of reservoir change during production, a thorough understanding of the influence of pore pressure on seismic velocity is critical. Laboratory measurements show that differential pressure (overburden minus fluid pressure) does not adequately determine the actual reservoir conditions. Changes in fluid pressure are found to have an additional effect on the physical properties of rocks. The effective‐stress coefficient n is used to quantify the effect of pore pressure compared to confining pressure on rock properties. However, the current practice in time‐lapse feasibility studies, reservoir‐pressure inversion and pore‐pressure prediction is to assume that n= 1. Laboratory measurements, reported in both this and previous research show that n can be significantly less than unity for low‐porosity rocks and that it varies with porosity, rock texture and wave type. We report the results of ultrasonic experiments to estimate n for low‐porosity sandstones with and without microcracks. Our results show that, for P‐waves, n is as low as 0.4 at a differential pressure of 20 MPa (about 3000 psi) for a low‐porosity sandstone. Thus, in pore‐pressure inversion, an assumption of n= 1 would lead to a 150% underestimation of the pore pressure. Comparison of the effective‐stress coefficient for fractured and unfractured samples suggests that the presence of microfractures increases the sensitivity of P‐wave velocity to pore pressure, and therefore the effective‐stress coefficient. Our results show that the effective‐stress coefficient decreases with the differential pressure, with a higher differential pressure resulting in a lower effective‐stress coefficient. While the effective‐stress coefficient for P‐wave velocity can be significantly less than unity, it is close to one for S‐waves.     

Joint elastic‐electrical effective medium models of reservoir sandstones

Improvements in the joint inversion of seismic and marine controlled source electromagnetic data sets will require better constrained models of the joint elastic‐electrical properties of reservoir rocks. Various effective medium models were compared to a novel laboratory data set of elastic velocity and electrical resistivity (obtained on 67 reservoir sandstone samples saturated with 35 g/l brine at a differential pressure of 8 MPa) with mixed results. Hence, we developed a new three‐phase effective medium model for sandstones with pore‐filling clay minerals based on the combined self‐consistent approximation and differential effective medium model. We found that using a critical porosity of 0.5 and an aspect ratio of 1 for all three components, the proposed model gave accurate model predictions of the observed magnitudes of P‐wave velocity and electrical resistivity and of the divergent trends of clean and clay‐rich sandstones at higher porosities. Using only a few well‐constrained input parameters, the new model offers a practical way to predict in situ porosity and clay content in brine saturated sandstones from co‐located P‐wave velocity and electrical resistivity data sets.     

Velocities of elastic waves andv P/v S ratios in dry and water-saturated rock samples at high pressure

Elastic wave velocities and the ratiov _P/v _S were studied for dry and initially saturated samples of carbonate and crystalline rocks at pressures to 2 kbar. In initially saturated samplesv _P increases in crystalline rock, whereas in sedimentary rock it may either increase or decrease with increasing pressure. Under the same conditionsv _S remains approximately constant in crystalline rocks and decreases in sedimentary samples. The ratiov _P/v _S as a function of pressure increases in dry rocks and decreases in saturated ones. Saturated samples always have higherv _P/v _S ratios than dry ones. It appears that the geometry of the pore space influences the acoustic properties of saturated crystalline rocks.