基于固体替换模型的有机质含量对页岩弹性性质的影响分析

页岩中的TOC(Total Organic Carbon,总有机碳)含量,对页岩的有效弹性模量以及与之相关的弹性波速度(P波和S波)有重要影响,建立弹性模量与TOC含量关系是页岩气甜点预测的重要手段之一.CS和SM两种固体置换理论主要针对孔隙度较大的砂岩,能否适用于孔隙度低、孔隙形态复杂和非均质性强的页岩目前尚未深入研究.鉴于目前已知的富有机质页岩的TOC赋存形态与裂缝以及孔隙形态类似,有关TOC含量对岩石弹性模量的影响可视为孔隙物质充填问题来研究.本文利用数字岩心技术,构造同一数字岩心不同TOC含量的样本群,基于CS和SM两种固体替换理论模型,通过有限元(FEM)数值模拟交叉验证,详细研究了两种固体替换方程对页岩的适用性和TOC含量对页岩弹性性质的影响.研究表明,由于实际岩心孔隙及TOC分布的非均质性,CS替换方程弹性模量预测值与FEM模拟结果存在差异,而SM替换方程预测值与FEM模拟结果基本一致,两种方程的预测差异揭示页岩非均质强度,利用SM替换方程中的参数α_1,α_2,β_1和β_2可详细分析实际岩心孔隙及TOC分布的非均质特征.     

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.     

龙马溪组页岩微观结构、地震岩石物理特征与建模

龙马溪组页岩是目前国内页岩气勘探的主要层位之一.由于岩石物理实验结果具有区域性,龙马溪组页岩的岩石特征与其地震弹性性质的响应规律需要开展相关的实验和理论研究工作予以明确.本研究基于系统的微观结构观察(扫描电镜和CT成像技术)和岩石物理实验来分析龙马溪组页岩样品地震弹性性质的变化规律,并依据微观结构特征建立相应的地震岩石物理表征模型.研究结果表明,石英含量对龙马溪组页岩的孔隙度以及有机碳(TOC)含量具有一定的控制作用,TOC和黄铁矿主要赋存于孔隙中;岩石骨架组成亦受控于石英或粘土含量,在石英含量大于40%(对应粘土含量小于30%)时,以石英、粘土共同作为岩石骨架,而粘土含量大于30%时,则以粘土作为岩石的骨架.因此,岩石骨架组成矿物、TOC含量、孔隙度共同制约龙马溪组页岩的地震弹性性质,富有机质储层岩石通常表现出低泊松比、低阻抗和低杨氏模量的特征,但由于支撑矿物的转换,某些富有机质页岩亦可表现为高阻抗特征.粘土矿物的定向排列仍然是造成页岩样品表现出各向异性的主要原因,各向异性参数与粘土含量具有指数关系.基于龙马溪组页岩的岩性特征及微观结构特征,可以利用自洽模型(SCA)、微分等效模量模型(DEM)和Backus平均模型的有效组合较为准确地建立龙马溪组页岩的地震岩石物理模型,实验结果和测井数据验证了模型的准确性.研究结果可为龙马溪组页岩气储层的测井解释和地震"甜点"预测提供依据.     

龙马溪组页岩数字岩心动态法弹性等效数值建模

数字岩心是计算岩石弹性性质的一类常用方法.龙马溪组页岩具有多矿物构成、复杂微结构和强非均质等特征,常规岩石物理的弹性等效解析建模局限性较大,目前流行的静态数值等效建模方法的精度有限.本文基于高分辨率的页岩数字岩心数据,采用多阈值分割方法将数字岩心分解为黏土、石英、孔隙、TOC、长石类、黄铁矿类等六种矿物类型;利用矿物组分等效模量法计算各类矿物的弹性模量;采用二元函数分水岭方法表征不同压力下的岩石孔隙变形和颗粒接触关系变化;通过取向分布函数(ODF)定量分析矿物颗粒展布造成的各向异性特征.最后基于Biot孔弹方程,采用不分裂卷积完全匹配层(CPML)旋转交错网格有限差分法模拟不同压力下弹性波在数字岩心中的传播.以未加压的数字岩心为参考模型,计算不同压力下弹性波走时的平均时间差,进而估算各压力点的数字岩心等效速度.与该岩心样品的超声实验测量速度比较,动态法数值计算结果略偏高,据此校正数值计算过程中表征岩石微结构及颗粒接触关系随压力变化的二元函数,有效改善动态法弹性等效数值建模精度.     

Numerical estimation of carbonate rock properties using multiscale images

Characterizing the pore space of rock samples using three‐dimensional (3D) X‐ray computed tomography images is a crucial step in digital rock physics. Indeed, the quality of the pore network extracted has a high impact on the prediction of rock properties such as porosity, permeability and elastic moduli. In carbonate rocks, it is usually very difficult to find a single image resolution which fully captures the sample pore network because of the heterogeneities existing at different scales. Hence, to overcome this limitation a multiscale analysis of the pore space may be needed. In this paper, we present a method to estimate porosity and elastic properties of clean carbonate (without clay content) samples from 3D X‐ray microtomography images at multiple resolutions. We perform a three‐phase segmentation to separate grains, pores and unresolved porous phase using 19 μm resolution images of each core plug. Then, we use images with higher resolution (between 0.3 and 2 μm) of microplugs extracted from the core plug samples. These subsets of images are assumed to be representative of the unresolved phase. We estimate the porosity and elastic properties of each sample by extrapolating the microplug properties to the whole unresolved phase. In addition, we compute the absolute permeability using the lattice Boltzmann method on the microplug images due to the low resolution of the core plug images. In order to validate the results of the numerical simulations, we compare our results with available laboratory measurements at the core plug scale. Porosity average simulations for the eight samples agree within 13%. Permeability numerical predictions provide realistic values in the range of experimental data but with a higher relative error. Finally, elastic moduli show the highest disagreements, with simulation error values exceeding 150% for three samples.     

Shales in the Qiongzhusi and Wufeng–Longmaxi Formations: a rock-physics model and analysis of the effective pore aspect ratio

The shales of the Qiongzhusi Formation and Wufeng–Longmaxi Formations at Sichuan Basin and surrounding areas are presently the most important stratigraphic horizons for shale gas exploration and development in China. However, the regional characteristics of the seismic elastic properties need to be better determined. The ultrasonic velocities of shale samples were measured under dry conditions and the relations between elastic properties and petrology were systemically analyzed. The results suggest that 1) the effective porosity is positively correlated with clay content but negatively correlated with brittle minerals, 2) the dry shale matrix consists of clays, quartz, feldspars, and carbonates, and 3) organic matter and pyrite are in the pore spaces, weakly coupled with the shale matrix. Thus, by assuming that all connected pores are only present in the clay minerals and using the Gassmann substitution method to calculate the elastic effect of organic matter and pyrite in the pores, a relatively simple rock-physics model was constructed by combining the self-consistent approximation (SCA), the differential effective medium (DEM), and Gassmann’s equation. In addition, the effective pore aspect ratio was adopted from the sample averages or estimated from the carbonate content. The proposed model was used to predict the P-wave velocities and generally matched the ultrasonic measurements very well.     

基于微结构-尺度双分解的页岩随机介质模拟和非均质特征分析

页岩储层矿物颗粒、孔/裂隙、干酪根等微观结构呈现明显的尺度化分布特征,常规的单结构单尺度随机介质模拟方法难以完整描述和重构微观尺度的页岩储层介质,本文提出了一种微结构-尺度双分解的随机介质模拟方法.基于龙马溪组页岩数字岩心,将岩心切片按照占比分解为脆性矿物、孔隙、干酪根及背景介质四种类型,对脆性矿物、孔隙和干酪根三种微结构进行尺度分解,通过优化随机介质模型参数,实现精确模拟不同尺度的微结构组分,再按占比进行微结构-尺度双合成.结果表明,微结构-尺度双分解随机介质模拟大幅度提高强非均质页岩储层介质的建模精度.     

数字岩心逆建模理论下的储层参数定量预测方法

数字岩心微观孔隙结构十分复杂,有限元模拟物性参数与弹性参数之间关系是非线性的,直接反演其物性参数准确度低、稳定性差.本文发展了一种数字岩石物理逆建模方法,实现了基于数字岩心的储层参数有效预测.从数字岩心基函数的构建出发,基于有限元方法,计算了一系列具有等间距物性参数值(孔隙度、泥质含量和含水饱和度)的数字岩心弹性参数(体积模量、剪切模量和密度),通过插值算法建立了数字岩心弹性参数三维数据集,从而实现了弹性模量的有限元数值解的快速构建;然后搜索弹性参数的单值等值面,通过等值面的空间交会得到交点,完成储层参数预测.测试结果表明:基于数字岩心逆建模理论的储层参数预测结果与实际模型一致,具有可行性,并且可以通过增加插值点数目提高预测的准确性;孔隙度和泥质含量预测结果稳定性很好,而含水饱和度对噪声的加入较为敏感.     

基于扫描电镜的页岩有机孔隙空间定量表征

对于富含有机质的页岩储层来说,孔隙度是评价油气储量的关键参数之一,有机孔隙的存在又对岩石孔隙空间具有重要贡献,因此,研究页岩储层有机孔隙空间是有必要的。本文结合氩离子抛光扫描电镜（SEM）资料和图像分割技术,以川南永川区块YY2井龙马溪组页岩储层为例,对纳米有机孔隙的形态特征和有机孔隙度进行分析研究,研究表明:YY2井龙马溪组有机孔隙较发育,SEM图像中有机孔隙直径差异很大,对龙马溪组6块岩样的SEM图像进行处理,得到的有机质面孔率在0.97％～8.90％之间,个别样品孔径较小,面孔率在1％以下。利用面孔率和有机质含量计算有机孔隙度,计算得到的有机孔隙度在0.12％～0.75％范围内,均值为0.39％,占据岩石孔隙度的10.66％。      

矿物组分对龙马溪组页岩动、静态弹性特征的影响

龙马溪组页岩是目前国内页岩气勘探的主要层位之一,其静态弹性性质是影响压裂效果的重要因素,而动态弹性性质则是页岩"甜点"地震预测的重要基础.但对龙马溪组页岩动、静态弹性特征相互关系及其影响因素缺少对比研究,致使运用动态弹性性质进行"甜点"预测时存在不确定性.在对龙马溪组页岩样品系统动、静态性质联合测量基础上,分析了页岩样品动、静态弹性性质的变化规律,并讨论了成岩作用与动、静态弹性性质变化规律的因果关系.研究结果表明,龙马溪组页岩上下两段成岩过程存在明显差异,致使上段页岩在结构上表现为以塑性黏土颗粒作为受力骨架,而下段页岩在结构上逐渐转变为脆性石英颗粒作为受力骨架.页岩岩石结构特征的变化控制了动、静态弹性特征的整体变化规律,表现为动、静态杨氏模量、峰值强度等随石英含量的增大表现出近似"V"型的变化形式,而与钙质含量变化呈现正相关关系,与黏土含量变化表现出负相关性.上段页岩宏观力学性质受微观石英+黏土颗粒集合体弹性性质控制,应力应变曲线表现出明显塑性段,动、静态杨氏模量比大于1.4,破裂易于发生在黏土与石英颗粒边界形成宏观单剪型破裂模式,样品脆性低,且脆性变化不受石英含量的影响.下段页岩宏观力学性质受微观石英颗粒集合体弹性性质控制,应力应变曲线表现为弹性变形,动、静态杨氏模量比小于1.4,破裂发生在石英颗粒内部形成宏观劈裂型破裂模式,样品脆性高,且脆性变化与石英含量呈明显正相关关系.研究结果可为龙马溪页气储层的测井解释和地震"甜点"预测提供依据.     

基于二维图像的数字岩心电导率计算方法研究

电导率是表征岩石电学性质的重要物理参数,在地质资源勘查和测井解释等领域发挥着巨大作用.快速、准确地确定岩石电导率具有重要的理论和实践意义.作为近年来发展的一种岩石物理数值模拟工具,数字岩心技术在定量计算电导率等物性参数方面应用广泛.三维微观结构的准确获取是数字岩心技术计算岩石电导率的关键,但传统获取岩石三维微观结构的方法较为复杂费时.为了方便快速地通过数字岩心技术计算岩石的电导率,本文研究了岩石二维与三维数字岩心的电导率联系.我们基于微米级X射线CT扫描得到的三个砂岩样品的微观结构信息建立了三维数字岩心,并通过有限元法计算的三维数字岩心电导率与实验数据的对比验证数值计算方法的有效性.随后我们数字地扩展了岩石的孔隙,产生了较大孔隙度的三维数字岩心样本,在此基础上,计算了三维数字岩心和相应二维数字岩心的电导率,并通过Archie公式分别拟合了电导率与孔隙度之间的关系,得到了相应的胶结系数.结果表明,三维数字岩心的胶结系数小于二维数字岩心的胶结系数,且二者的比值与岩石实测孔隙度呈线性负相关关系.以该联系为纽带,通过二维图像快速计算得到的电导率与孔隙度关系,确定了三维数字岩心的电导率与孔隙度关系,并进一步通过三维数字岩心的孔隙度计算其电导率.该方法计算得到的人工砂岩样品的电导率与其三维数字岩心电导率相关系数高于96%,验证了基于二维图像的数字岩心电导率计算方法的有效性.本文的研究结果为快速、准确地计算岩石电导率提供了新的思路,在油气勘探开发中有广阔的应用前景.     

Quantitative characterization of reservoir space in the Lower Silurian Longmaxi Shale,southern Sichuan,China

Based on the drilling data of the Upper Ordovician Wufeng Shale and the Lower Silurian Longmaxi Shale in southern Sichuan Basin,the construction of matrix pores and the development condition of fractures in a marine organic-rich shale are quantitatively evaluated through the establishment of the reservoir petrophysical models and porosity mathematical models.Our studies show that there are four major characteristics of the Longmaxi Shale confirmed by the quantitative characterization:(1)the pore volume of per unit mass is the highest in organic matter,followed in clay minerals,finally in brittle minerals;(2)the porosity of the effective shale reservoir is moderate and equal to that of the Barnett Shale,and the main parts of the shale reservoir spaces are interlayer pores of clay minerals and organic pores;(3)the porosity of the organic-rich shale is closely related to TOC and brittle mineral/clay mineral ratio,and mainly increases with TOC and clay mineral content;(4)fractures are developed in this black shale,and are mainly micro ones and medium-large ones.In the Longmaxi Shale,the fracture density increases from top to bottom,reflecting the characteristics with high brittle mineral content,high Young’s modulus,low Poisson's ratio and high brittleness at its bottom.     

Application of alternative digital rock physics methods in a real case study: a challenge between clean and cemented samples

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.     

基于数字岩心的碳酸盐岩孔隙结构对弹性性质的影响研究(上篇):图像处理与弹性模拟

碳酸盐岩复杂的孔隙结构如何影响其弹性性质一直是地球物理研究的难点问题,在此基础上如何半定量甚至是定量地对碳酸盐岩储层预测,特别是如何有效地获取孔隙结构参数相关的地震属性体一直是油气工业界追求的目标.本研究从数字岩心角度入手,联合测井以及地震数据尝试探究这一问题的解决方案,包括如下几个方面:(1)代表性碳酸盐岩储层样品获取;(2)CT扫描数字岩心数据体获取;(3)数字岩心数据的图像处理;(4)数字岩心数据的静态弹性模拟;(5)数字岩心子数据体的孔隙结构因子提取;(6)孔隙结构因子表征与分类下的弹性性质与孔隙度的定量化量版建立;(7)数字岩心-井-地震联合的孔隙度属性提取;(8)孔隙结构因子的地震属性体获取.本研究分为两篇系列文章上篇与下篇,上篇主要阐述如上提出的(1)—(4)方面,重点在于针对碳酸盐岩二值化图像处理的流程建立与验证,以及数字岩心静态弹性模拟的理论方面,这两方面是基于数字岩心获得精确的碳酸盐岩弹性性质模拟结果的关键所在;下篇主要阐述利用数字岩心数据获得孔隙结构因子的思路、理论与流程,以及为碳酸盐岩储层预测为目标而获得孔隙结构因子的地震属性体的实际应用方面.由于两篇文章共享数字岩心数据,同时所涉及的研究思路与流程形成一个有机整体,因此写成两篇系列文章而非两篇独立文章.本文为两篇系列文章的第一篇:上篇.     

龙马溪组页岩数字岩心LSM-RVM数值建模方法研究及TOC含量影响分析

页岩气储层中含有大量有机碳(TOC),其丰度与成熟度对页岩力学特性有重要影响.建立包含TOC的精细数值模型,将有助于探索页岩微结构与矿物组分含量对等效弹性模量的作用程度,是“甜点区”预测的重要理论基础.本文提出了一种离散数值建模方法,基于高精度成像技术,采用晶格弹簧-随机孔隙耦合模型(LSM-RVM)模拟包含多种矿物组分及不同成熟度干酪根的数字岩心,分析TOC成熟度及含量对弹性参数的影响.在该模型中,参数设置(数值阻尼与加载应变速率)至关重要,选取不当会对计算精度造成一定影响.研究结果表明,LSM-RVM能够生成符合TOC及多种矿物实际分布特征的数值模型,是一种精细数值建模方法.     

基于Curvelet变换与POCS方法的三维数字岩心重建

随着页岩气勘探与开发的深入,研究页岩裂隙的三维空间展布成为页岩岩石物理研究的必要步骤之一.但由于仪器的限制,页岩切片在深度上具有不连续性,以及数字岩心纵向上成像最小间隔与横向分辨率的不一致成为影响裂隙表征和数字岩石物理模拟精度提高的重要因素.为了更好的研究裂隙在三维的空间展布,本文将curvelet稀疏变换与凸集投影(POCS)迭代算法有效结合,实现三维数字岩心重建.首先对X射线扫描砂岩得到的三维数据体进行隔片抽稀,利用本文方法实现三维数据体重建,重建结果与完整数据体具有很好的一致性,且优于现有方法(spgl1),验证了新方法的有效性与先进性.其次对聚焦离子束扫描电镜(FIB-SEM)得到的纳米级页岩二维切片在深度上进行了加密重建,获得纵向上成像最小间隔与横向分辨率基本一致的三维数字岩心,由于仪器限制引起的页岩切片深度上的不连续性得到减弱,裂隙展布更加清晰.砂岩CT图像以及页岩FIB-SEM成像数据的重建结果验证了本文方法的有效性与先进性.     

Multiscale study for stochastic characterization of shale samples

Characterization of shale reservoirs, which are typically of low permeability, is very difficult because of the presence of multiscale structures. While three-dimensional (3D) imaging can be an ultimate solution for revealing important complexities of such reservoirs, acquiring such images is costly and time consuming. On the other hand, high-quality 2D images, which are widely available, also reveal useful information about shales’ pore connectivity and size. Most of the current modeling methods that are based on 2D images use limited and insufficient extracted information. One remedy to the shortcoming is direct use of qualitative images, a concept that we introduce in this paper. We demonstrate that higher-order statistics (as opposed to the traditional two-point statistics, such as variograms) are necessary for developing an accurate model of shales, and describe an efficient method for using 2D images that is capable of utilizing qualitative and physical information within an image and generating stochastic realizations of shales. We then further refine the model by describing and utilizing several techniques, including an iterative framework, for removing some possible artifacts and better pattern reproduction. Next, we introduce a new histogram-matching algorithm that accounts for concealed nanostructures in shale samples. We also present two new multiresolution and multiscale approaches for dealing with distinct pore structures that are common in shale reservoirs. In the multiresolution method, the original high-quality image is upscaled in a pyramid-like manner in order to achieve more accurate global and long-range structures. The multiscale approach integrates two images, each containing diverse pore networks – the nano- and microscale pores – using a high-resolution image representing small-scale pores and, at the same time, reconstructing large pores using a low-quality image. Eventually, the results are integrated to generate a 3D model. The methods are tested on two shale samples for which full 3D samples are available. The quantitative accuracy of the models is demonstrated by computing their morphological and flow properties and comparing them with those of the actual 3D images. The success of the method hinges upon the use of very different low- and high-resolution images.     

基于时频电磁法的富有机质页岩层系勘探研究

我国地质条件复杂、多样,能否充分发挥出非地震勘探技术成本低、效率高的技术优势,检验电法勘探技术在页岩气勘探中的有效性,是当前在页岩气勘探起步阶段亟待解决的重要问题之一.本文介绍了在研究和总结含气页岩密度、极化率、电阻率等岩石物理特征基础上,在四川盆地南部筠连地区开展的物性调查、时频电磁法勘探试验工作.勘探研究结果表明,本地区分布的富有机质页岩层系-志留系龙马溪组(S1l)具备开展电法勘探工作的物性条件,时频电磁法具有勘探富有机质页岩层系的能力.     

Modelling the poroelasticity of rocks and ice

For many geophysical and astrophysical applications the relationship between physical bulk properties and porosity is needed. As part of the preparation for the 'Rosetta' comet rendezvous mission, a simple model for the porosity dependence of the elastic properties of granular media, i.e. the elastic moduli and the propagation velocities of elastic waves, has been developed based on textural properties and the contact stiffness of the constituent particles. It is shown that the derived relationships fit very well with sandstone data. The model is also consistent with data for snow and ice and is in agreement with the transformation mechanisms from snow to ice. A short review shows the relevance to another physical bulk property, the thermal conductivity.