泥质分布形式对泥质砂岩电性的影响规律研究

介绍了将混合理论用于结构泥质、将基于距离的升尺度(DBU)方法用于层状泥质电性研究的原理和方法,通过对实际岩心资料以及井资料的分析处理,证实了将混合理论用于结构泥质、将DBU方法用于层状泥质砂岩电性研究的合理性,最后利用上述方法进一步分析了泥质含量、泥质的分布形式等对泥质砂岩电阻率的影响及规律,研究结果表明:随着泥质含量的增加,泥质砂岩电阻率逐渐降低,且降低的幅度随含水饱和度的降低而增大;当测量电流方向与层状泥质垂直时,泥质对电性的影响较弱,但层状泥质和结构泥质的相对含量对结果的影响很大;当测量电流方向与层状泥质平行时,泥质对电性的影响较强,但层状泥质和结构泥质的相对含量对结果的影响不大.     

骨架导电的混合泥质砂岩通用孔隙结合电阻率模型研究

在水介质中顺序添加分散粘土颗粒、油珠、导电骨架颗粒、层状泥质,并对每一种成分进行连续积分,建立了一种适用骨架导电及含有分散粘土和层状泥质的泥质砂岩通用电阻率模型.通过对该模型的影响因素分析,发现泥质分布形式对模型计算的含水饱和度有很大影响;对应两个不同粘土颗粒电阻率或骨架颗粒电阻率的地层电导率之差,几乎与总含水饱和度无关,而对应两个不同层状泥质电阻率的地层电导率之差,随总含水饱和度增大而增大;骨架胶结指数变化对地层电导率与总含水饱和度关系曲线的影响最大,而粘土胶结指数变化对地层电导率与总含水饱和度关系曲线的影响最小;饱和度指数对地层电导率与总含水饱和度关系曲线的影响随总含水饱和度的增大而减小.通过一组骨架导电的人造岩样的试验,表明当地层水电阻率.小于颗粒电阻率时,该模型可以用于不含粘土的骨架导电的岩石.通过两组分散泥质砂岩岩样实验测量数据和一组层状泥质砂岩测井资料及实际测井资料的计算,表明本文给出的电阻率模型既适用于分散泥质砂岩地层解释又适用于层状泥质砂岩地层解释,同时,还适用于含有分散粘土和层状泥质的混合泥质砂岩地层解释.     

不同泥质分布形式泥质砂岩导电规律实验研究

本文利用人工制作的不同含量分散泥质和层状泥质砂岩岩心样品,测量不同矿化度和不同含油饱和度的岩心电阻率,从实验角度研究了不同泥质分布形式和含量的岩心导电规律,结果表明,泥质分布形式或含量不同,则泥质砂岩导电规律不同.基于层状泥质与分散泥质砂岩的并联导电实验规律,以及分散粘土和地层水混合物的导电规律可用HB电阻率方程描述,建立了考虑泥质分布形式影响的泥质砂岩电阻率模型.通过1组不同泥质分布形式泥质砂岩人造岩心实验数据的测试,表明该模型可以描述分散泥质砂岩、层状泥质砂岩和混合泥质砂岩地层的导电规律.分散泥质,层状泥质,人造岩样,实验测量,并联导电,HB方程,电阻率模型     

泥质砂岩液相渗透率计算新方法

液相渗透率描述了岩石的渗流特性,是评价储层与预测油气产能的重要参数.液相渗透率是指盐水溶液在岩石孔隙中流动且与岩石孔隙表面黏土矿物发生物理化学作用时所测得的渗透率;液相渗透率的实验测量条件更加接近实际地层泥质砂岩的条件,使得液相渗透率更能反映地层条件下泥质砂岩的渗流特性;然而,现有的液相渗透率评价模型较少,且模型未能揭示液相渗透率与溶液矿化度之间的关系.基于此,开展了液相渗透模型推导与计算方法研究;文中首先将岩石等效为毛管束模型,推导建立了液相渗透率与比表面、喉道曲折度、总孔隙度、黏土束缚水孔隙度等参数之间的关系;其次,根据岩石物理体积模型,推导建立了黏土束缚水孔隙度与阳离子交换容量、溶液矿化度等参数的关系;最终,将黏土束缚水孔隙度引入液相渗透率计算公式,建立了基于总孔隙度、阳离子交换容量、溶液矿化度、比表面、喉道曲折度等参数的液相渗透率理论计算模型.液相渗透率计算模型与两组实验数据均表明,液相渗透率随阳离子交换容量的增大而降低,随溶液矿化度的增大而增大.然而,液相渗透率理论计算模型的实际应用中喉道曲折度、比表面等参数求取困难,直接利用理论模型计算液相渗透率受到限制.在分析液相渗透率与孔隙渗透率模型的基础上,建立了液相渗透率与空气渗透率之间的转换模型,形成了利用转化模型计算液相渗透率的新方法.为进一步验证液相渗透率与空气渗透率转化模型的准确性,基于两组实验数据,利用转换模型计算了液相渗透率;液相渗透率计算结果与岩心测量液相渗透率实验结果对比显示,液相渗透率计算结果与实际岩心测量结果吻合较好,文中建立的液相渗透率与空气渗透率转化模型合理可靠.     

利用毛管模型研究泥质砂岩电化学测井响应机理

自然电位和激发极化电位测井响应所涉及的离子导体激发极化电位的微观机理解释,主要依据双电层形变假说和浓差极化假说,缺少定量描述的数学模型和理论体系.本文利用孔隙介质的微观毛管模型,给出了毛管模型中双电层理论和阳离子交换量与Zeta电位的关系,推导出毛管中离子流量和电流强度表达式.由电荷守恒定律和物质守恒定律,推导出毛管中离子浓度分布的解析表达式,建立了描述含水泥质砂岩激发极化电位和自然电位的数学模型.从而系统地严格证明了含水泥质砂岩激发极化现象是在电流场和浓度梯度场的共同作用下,由孔隙中离子浓度浓差极化电位和双电层形变电位形成的.并且证明了描述泥质砂岩自然电位的数学方程和描述激发极化电位的数学方程及形成机理是一致的.计算结果表明：激发极化极化率随孔隙度和渗透率的增大而减小;极化率随溶液浓度的增加而减小,随阳离子交换量的增加而增加;证明了地层水浓度、阳离子交换量是影响自然电位大小的主要因素.     

油水饱和泥质砂岩流动电位的频散特性

油水饱和泥质砂岩中流动电位的研究对于揭示含油储层震电勘探和动电测井的机理有着重要的意义.本文首先从岩石孔隙的微观结构出发,构造了描述水润湿条件下油水饱和泥质砂岩储层的毛管模型.在模型中依据油水流动遵守的Navier-Stokes方程和电化学传质动力学理论,建立了描述油水饱和泥质砂岩流动电位的数学方程,并数学模拟了岩石储渗参数对流动电位频散特性的影响规律.研究结果表明:储层孔隙内流体受到的粘滞力与惯性力控制着水相和油相的流动,从而决定了流动电位的频散特性.随着孔隙度的增大,油水两相各自的有效渗透率均增大;而含水饱和度的升高使得水相有效渗透率增大,油相有效渗透率减小.在水润湿条件下,流动电位耦合系数随含水饱和度升高而增大,随束缚水饱和度的升高而减小.另外,流动电位相对耦合系数也随含水饱和度的升高而增大,但无频散现象.     

不同矿化度下泥质对岩石电性影响的逾渗网络研究

为了进一步认识泥质对岩石宏观电性的影响规律,本文利用非规整三维逾渗网络模型,通过数值模拟研究了不同矿化度下泥质对岩石电性影响的规律. 模拟结果表明：在中低矿化度下,泥质对岩石导电整体上呈现减阻作用,随着泥质含量的增加,电阻率降低的速度减慢;中等矿化度下泥质的减阻效果明显弱于低矿化度下的减阻效果;在高矿化度下泥质对岩石导电整体上呈现增阻作用. 在高矿化度、高含水饱和度下泥质对岩石电性的影响较小. 泥质起减阻、增阻作用的具体矿化度范围取决于储层的孔隙度、连通性以及地层温度等特性.     

砂岩储层AVO特征影响因素的不确定性研究

传统的地震AVO正演研究多采用参数固定的岩石物理模型,而实际地层属性参数在勘探范围内具有不确定性.本研究以目标地层岩芯样品的实验室测试数据为基础,通过样品孔隙度和干燥状态下纵、横波阻抗的高度线性关系对岩石物理模型进行了简化,并结合实验测量和测井解释建立了主要模型参数的概率密度函数.采用Monte-Carlo随机正演和Gassmann流体替换技术,对模型参数不确定性与地震AVO响应关系的研究表明,模型饱和不同流体时的AVO截距-斜率点束的分布形态和偏离饱水背景趋势线的程度受盖层泥岩速度的不确定性影响非常大,而储层砂岩孔隙度的不确定性也有明显影响.正确解释实际地震AVO异常的关键,是将岩石物理模型的参数合理地概率化,通过正演模拟来获得储层物性和流体信息的最大可能性.     

致密砂岩气藏裂隙-孔隙型弹性岩石物理模板研究:以川西坳陷A区为例

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

流体黏度对砂岩弹性模量频散与衰减影响规律的实验及理论验证

波诱导的孔隙流体流动是造成地震波速度频散与衰减的重要机制, 其中孔隙流体性质(尤其是黏度)是影响速度频散变化特征的关键因素之一.为了研究孔隙流体黏度对速度频散的影响, 本文在饱和不同黏度流体条件下对2块低孔隙度砂岩和3块中孔隙度砂岩样品进行了低频应力-应变岩石物理测试, 获得了岩石样品1~3 kHz的杨氏模量、泊松比和衰减曲线.实验结果发现: 两种砂岩样品都表现出随着有效压力增加, 杨氏模量增加, 频散程度减弱的特征; 同时随着流体黏度增大, 频散梯度增大, 特征频率向较低频率移动.为了更好地解释不同流体黏度条件下的实验数据, 本文采用现有经典岩石物理频散模型对实测速度频散进行了描述.研究结果表明, 经典的Gurevich喷射流模型、Biot模型、White模型、Gassmann方程和BISQ模型在数值上和趋势上均难以定量刻画实测低频数据, 其预测结果存在明显的特征频率且中间过渡频段较小, 而实测速度则随频率的增大呈缓慢递增趋势, 且具有较宽的中间过渡频段.对于中孔隙度砂岩样品, 考虑Biot宏观流作用、软孔与软孔以及软孔与硬孔喷射流效应的扩展Gurevich喷射流模型能够给出与实验数据较为吻合的预测结果; 对于低孔隙度砂岩, 扩展Gurevich喷射流模型的预测结果也能在趋势上与实验数据保持一致.实验数据和理论模型共同指出归一化参数E(f)/E(1 Hz)具有提高储层流体预测与流体识别精度的潜力和重要意义.

     

孔隙-裂隙性质和黏土含量对页岩油层系致密岩石弹性波耗散规律的影响研究

作为一种非常规油气资源,页岩油储量丰富、分布范围广,具有巨大的勘探开发潜力,是近年来油气产业关注的重点与热点.然而,页岩油储层具有岩石矿物组分多样、低孔低渗、孔隙结构复杂、非均质性强等特征,与常规油气资源存在明显差异.本研究选取鄂尔多斯盆地中生界延长组长7油层组的10块致密砂岩样本,基于X射线衍射分析得到各样本的矿物组分,开展不同围压和流体条件下的超声波实验观测,进而获得样本的纵、横波速度和纵波衰减逆品质因子.基于实验测量获得变压力条件下的孔隙度,结合线性外推的方法,估算各样本的裂隙孔隙度,代入EIAS（Equivalent Inclusion-Average Stress,等效嵌入体应力平均）模型,求取对应的裂隙纵横比和裂隙密度,分析页岩油储层孔隙-裂隙性质对纵波衰减的影响.结果显示相对于衰减,致密砂岩总孔隙度、裂隙纵横比、裂隙密度和衰减变化量（不同围压下的衰减观测值与最大围压下的衰减观测值的差）之间的相关性更加明显.基于薄片分析,结果显示致密样本存在孔内黏土包体、微裂隙包体和粒间孔的三重孔隙结构,因此本文引入三重孔隙结构模型,定量估算各样本的孔内黏土含量,进而分析孔内黏土含量及总黏土含量和纵波衰减之间的关系.结果显示孔内黏土含量是主导页岩油储层纵波衰减大小的主要因素之一,而非总黏土含量.本研究可为页岩油储层衰减特征分析、岩石物理模型构建及地震勘探方法研究提供理论支撑.

     

Water retention effects on elastic properties of Opalinus shale

Shales play an important role in many engineering applications such as nuclear waste, CO₂ storage and oil or gas production. Shales are often utilized as an impermeable seal or an unconventional reservoir. For both situations, shales are often studied using seismic waves. Elastic properties of shales strongly depend on their hydration, which can lead to substantial structural changes. Thus, in order to explore shaly formations with seismic methods, it is necessary to understand the dependency of shale elastic properties on variations in hydration. In this work, we investigate structural changes in Opalinus shale at different hydration states using laboratory measurements and X-ray micro-computed tomography. We show that the shale swells with hydration and shrinks with drying with no visible damage. The pore space of the shale deforms, exhibiting a reduction in the total porosity with drying and an increase in the total porosity with hydration. We study the elastic properties of the shale at different hydration states using ultrasonic velocities measurements. The elastic moduli of the shale show substantial changes with variations in hydration, which cannot be explained with a single driving mechanism. We suggest that changes of the elastic moduli with variations in hydration are driven by multiple competing factors: (1) variations in total porosity, (2) substitution of pore-filling fluid, (3) change in stiffness of contacts between clay particles and (4) chemical hardening/softening of clay particles. We qualitatively and quantitatively analyse and discuss the influence of each of these factors on the elastic moduli. We conclude that depending on the microstructure and composition of a particular shale, some of the factors dominate over the others, resulting in different dependencies of the elastic moduli on hydration.     

致密砂岩气储层的岩石物理模型研究

根据鄂尔多斯盆地苏里格气田以往实测和新测的共17口井51块岩样超声波实验数据,得到304组不同孔隙度和不同含水饱和度下对应的纵横波速度、泊松比等弹性参数.重新优选计算体积模量和泊松比与含气饱和度的关系,表明苏里格气田上古生界二叠系石盒子组盒8致密砂岩储层的模型与Brie模型（e=2）相似度最高.由此建立的苏里格气田储层岩石物理模型,更好的表征了致密岩石储层物理参数随含气饱和度变化规律,为该区储层预测提供了理论依据.致密储层岩石物理模型研究成果应用于苏里格气田多波地震资料气水预测中,实际例子表明该模型适用于该区的储层和含气性预测,并取得了较好的效果.     

Forced oscillation measurements of seismic wave attenuation and stiffness moduli dispersion in glycerine-saturated Berea sandstone

Fluid pressure diffusion occurring on the microscopic scale is believed to be a significant source of intrinsic attenuation of mechanical waves propagating through fully saturated porous rocks. The so-called squirt flow arises from compressibility heterogeneities in the microstructure of the rocks. To study squirt flow experimentally at seismic frequencies the forced oscillation method is the most adequate, but such studies are still scarce. Here we present the results of forced hydrostatic and axial oscillation experiments on dry and glycerine-saturated Berea sandstone, from which we determine the dynamic stiffness moduli and attenuation at micro-seismic and seismic frequencies (0.004–30 Hz). We observe frequency-dependent attenuation and the associated moduli dispersion in response to the drained–undrained transition (∼0.1 Hz) and squirt flow (>3 Hz), which are in fairly good agreement with the results of the corresponding analytical solutions. The comparison with very similar experiments performed also on Berea sandstone in addition shows that squirt flow can potentially be a source of wave attenuation across a wide range of frequencies because of its sensitivity to small variations in the rock microstructure, especially in the aspect ratio of micro-cracks or grain contacts.     

基于储层砂岩微观孔隙结构特征的弹性波频散响应分析

储层砂岩微观孔隙结构特征不仅影响干燥岩石的弹性波传播速度,也决定了岩石介质中与流体流动相关的速度频散与衰减作用.依据储层砂岩微观结构特征及速度随有效压力变化的非线性特征,将其孔隙体系理想化为不同形状的硬孔隙(纵横比α0.01)与软孔隙(纵横比α0.01)的组合(双孔隙结构).基于孔弹性理论,给出软孔隙最小初始纵横比值(一定压力下所有未闭合软孔隙在零压力时的纵横比最小值)的解析表达式,并在此基础上利用岩石速度-压力实验观测结果给出求取介质中两类孔隙纵横比及其含量分布特征的方法.通过逐步迭代加入软孔隙的方法对基于特征纵横比的"喷射流"(squirt fluid)模型进行了扩展,以考虑复杂孔隙分布特征对岩石喷射流作用的影响及其可能引起的速度频散特征.相较于典型的喷射流作用速度频散模式,对于岩石中软孔隙纵横比及其对应含量在较宽的范围呈谱分布的一般情况,其速度频散曲线不存在明显的低频段和中间频段,速度随频率的增大呈递增趋势直至高频极限.这说明即使在地震频段,微观尺度下的喷射流作用仍起一定作用,同样会造成流体饱和岩石介质的地震速度与Gassmann方程预测结果有不可忽略的差异.本文是对现有喷射流模型的重要补充,也为利用实验数据建立不同频段间岩石弹性波传播速度的可能联系提供了理论依据.     

Measurements of the elastic and anelastic properties of sandstone flooded with supercritical CO2

The aim of this study was to investigate the effects of supercritical CO₂ (scCO₂) injection on the elastic and anelastic properties of sandstone at seismic and ultrasonic frequencies. We present the results of the low‐frequency and ultrasonic experiments conducted on water‐saturated sandstone (Donnybrook, Western Australia) flooded with scCO₂. The sandstone was cut in the direction perpendicular to a formation bedding plane and tested in a Hoek triaxial pressure cell. During the experiments with scCO_2, the low‐frequency and ultrasonic systems and the pump dispensing scCO₂ were held at a temperature of 42°C. The elastic parameters obtained for the sandstone with scCO₂ at seismic (0.1 Hz–100 Hz) and ultrasonic (～0.5 MHz) frequencies are very close to those for the dry rock. The extensional attenuation was also measured at seismic frequencies for the dry, water‐saturated, and scCO₂‐injected sandstones. The applicability of Gassmann's fluid substitution theory to obtained results was also tested during the experiments.     

The impact of different clay minerals on the anisotropy of clay matrix in shale

Although clay is composed of disconnected anisotropic clay platelets, many rock physics models treat the clay platelets in shale as interconnected. However, the clay matrix in shales can be modelled as anisotropic clay platelets embedded within a soft isotropic interplatelet region, allowing the influence of disconnected clay platelets on the elastic properties of the clay matrix to be analysed. In this model, properties of the interplatelet region are governed by its effective bulk and shear moduli, whereas the effective properties of the clay platelets are governed by their volume fraction, aspect ratio and elastic stiffness tensor. Together, these parameters implicitly account for variations in clay and fluid properties, as well as fluid saturation. Elastic stiffnesses of clay platelets are obtained from the literature, including both experimental measurements and first-principles calculations of the full anisotropic (monoclinic or triclinic) elastic stiffness tensors of layered silicates. These published elastic stiffness tensors are used to compile a database of equivalent transverse isotropic elastic stiffness tensors, and other physical properties, for eight common varieties of layered silicates. Clay matrix anisotropy is then investigated by examining the influence of these different elastic stiffnesses, and of varying model parameters, upon the effective transverse isotropic elastic stiffness tensor of the clay matrix. The relationship between the different clay minerals and their associated anisotropy parameters is studied, and their impact on the resulting anisotropy of the clay matrix is analysed.     

Effect of sub‐core scale heterogeneities on acoustic and electrical properties of a reservoir rock: a CO2 flooding experiment of brine saturated sandstone in a computed tomography scanner

The effect of sub‐core scale heterogeneity on fluid distribution pattern, and the electrical and acoustic properties of a typical reservoir rock was studied by performing drainage and imbibition flooding tests with CO₂ and brine in a laboratory. Moderately layered Rothbach sandstone was used as a test specimen. Two core samples were drilled; one perpendicular and the other parallel to the layering to allow injection of fluids along and normal to the bedding plane. During the test 3D images of fluid distribution and saturation levels were mapped by an industrial X‐ray CT‐scanner together with simultaneous measurement of electrical resistivity, ultrasonic velocities as well as amplitudes. The results showed how the layering and the flooding direction influenced the fluid distribution pattern and the saturation level of the fluids. For a given fluid saturation level, the measured changes in the acoustic and electrical parameters were affected by both the fluid distribution pattern and the layering orientation relative to the measurement direction. The P‐wave amplitude and the electrical resistivity were more sensitive to small changes in the fluid distribution patterns than the P‐wave velocity. The change in amplitude was the most affected by the orientation of the layering and the resulting fluid distribution patterns. In some instances the change due to the fluid distribution pattern was higher than the variation caused by the change in CO₂ saturation. As a result the Gassmann relation based on ‘uniform' or ‘patchy' saturation pattern was not suitable to predict the P‐wave velocity variation. Overall, the results demonstrate the importance of core‐imaging to improve our understanding of fluid distribution patterns and the associated effects on measured rock‐physics properties.