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

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

利用混合法构建三维数字岩心(英文)

在弹性波有三维数字岩心描述了岩石的微观孔隙结构。X射线CT扫描是获取三维数字岩心最准确和直接的方法,但实验成本高。本文结合沉积过程模拟和模拟退火算法,提出了重建三维数字岩心的混合法,基于岩石二维图像重建三维数字岩心。利用岩石颗粒沉积算法构建初始数字岩心,作为模拟退火算法的初始状态。运用模拟退火算法调整岩石颗粒和孔隙的相对位置,使三维数字岩心与岩心二维图像具有相似的自相关函数,从而建立三维数字岩心。与传统模拟退火算法相比,该方法计算时间明显减小。运用局部孔隙度理论定量比较了重建数字岩心和岩心X射线微CT图像,两种数字岩心具有相似的均质性和孔隙连通性。利用有限元方法和格子玻尔兹曼方法分别模拟了重建三维数字岩心的地层因素和渗透率,数值模拟结果与实验结果相符。相比传统模拟退火算法,混合法重建数字岩心的传导特性更接近真实岩心的传导特性。     

CPU-GPU异构并行计算在数字岩心线弹性静力学有限元模拟中的应用

线弹性静力学有限元模拟是计算数字岩心弹性模量的有效方法之一,已被用于研究数字岩心的弹性模量与其微观结构、物质组成之间的关系.若要形成解决实际问题的能力必须计算足够多的数字岩心样本,而有限元模拟的计算量较大,因此并行计算对于该方法的成功应用非常重要.本文将数字岩心线弹性静力学有限元模拟算法分解为在CPU和GPU上执行的两个部分,由CPU负责协调控制,GPU负责大规模数值计算,实现CPU-GPU异构并行计算,获得计算效率提升.采用该并行算法计算孔隙数字岩心、裂缝数字岩心和裂缝—孔隙数字岩心的弹性模量,得到的弹性模量—孔隙度关系符合一般的岩石物理规律.CPU-GPU异构并行的线弹性静力学有限元模拟能够迅速计算大量数字岩心的弹性模量,提供相当于物理实验的"观测数据",对岩石物理学研究具有重要的意义.     

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

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

三维数字岩心技术岩石物理应用研究进展

三维数字岩心可在孔隙尺度上对岩石的微观结构进行精细表征,以数字岩心为载体的岩石微观结构分析与岩石物理属性模拟研究已成为岩石物理分析的重要发展方向之一.本文从数字岩心建模、数字岩心图像分析及岩石物理属性模拟三个方面对数字岩心技术进行了介绍.将数字岩心技术在岩石物理领域的应用分为两大类:基于三维数字图像的直接分析和基于三维数字岩心的岩石物理属性数值模拟.通过三维图像分析可获得孔隙结构、矿物成分、粒度分布等信息;通过岩石物理属性模拟可以研究储层岩石电性、弹性、渗流和核磁共振特征.多种手段相结合建立多尺度、多组分数字岩心,突破分辨率与样品尺度的矛盾限制是数字岩心技术未来的发展趋势.数字岩心技术已发展成为岩石物理实验的重要组成部分,如何将孔隙尺度得到的岩石物理参数升尺度到数字井筒、数字油藏中是数字岩心技术发展亟需解决的难题.     

三维数字岩心透射超声波模拟与速度精度分析（英文）

基于数字岩心的透射超声波模拟是数字岩石物理的重要组成部分,可用于研究波在孔隙岩石中的传播规律,计算岩石的等效速度等。在三维数字岩心中模拟超声波传播时,沿波传播方向在数字岩心两侧附加均匀的镶边层,设置震源和检波器。震源激发后检波器记录数字岩心的入射波和透射波,拾取两者波峰时间的差值,计算超声波在数字岩心中的传播速度(即数字岩心的速度)。为评估透射超声波模拟的精度,将一个数字岩心的孔隙分别饱和气、油和水,计算它们的速度。当震源频率较高时,速度随频率的降低而升高,这应当是波场散射的结果。当孔隙流体按气→油→水方向变化时,不同频率对应的速度具有相似的变化规律,并与线弹性静力学模拟输出的速度的变化规律相似,但数值有差异。考虑到线弹性静力学模拟已被广泛应用,若将其作为基准,则可认为透射超声波模拟具有较高的相对精度。     

基于计算岩石物理方法的页岩储层弹性参数提取

发展了应用数值计算方法获取页岩储层的速度、各向异性参数的计算岩石物理系列方法.该系列方法包括了大尺度精细地质模型数值建模、计算网格尺度的地球物理建模和地震波数值模拟提取岩石物理弹性参数.本文方法利用储层的统计数据而不是具体岩心的测量数据,可获得储层岩石物理弹性参数的变化规律.相比于基于岩心测试的岩石物理方法,本文方法可精细考虑实际储层的非均匀特征,可得到岩心测试难以求取的与尺寸效应高度相关的弹性参数,也避免了求取弹性参数变化规律时获取不同地质特征岩心的困难.本文发展了计算岩石物理方法,为计算岩石物理面临的大尺度地质建模和计算能力限制问题提供了有效的解决方案.文中以胜利罗家的页岩储层为例,求得了储层TOC含量从3%到21%变化情况下储层的P波、S波速度以及各向异性参数变化规律.     

基于数字岩心的岩石电性微观数值模拟(英文)

本文利用x射线CT获取反映岩心微观结构的三维数字岩心,利用数学形态学中的开运算模拟了岩石的油驱水排驱过程中,不同含水饱和度下油和地层水在孔隙空间中的分布。利用有限元方法计算了岩石电阻率,进而得到岩石地层因素和电阻率指数,并考查了岩石润湿性对岩石电阻率指数的影响。数值模拟结果表明：基于数字岩心的水湿岩石地层因素和电阻率指数数值模拟结果与实验结果一致,拓展了岩石电阻率实验的能力;岩石润湿性对岩石电性有重要影响,在相同含水饱和度下,油湿岩石电阻率高于水湿岩石电阻率,油湿岩石饱和度指数远大于水湿岩石饱和度指数。     

过程法数字岩心建模与弹性参数模拟研究

在数字岩石物理中数字岩心是弹性参数模拟的基础,现有的数字岩心建模方法虽然较多,但一般不能满足弹性参数模拟对骨架所提出的要求.过程法以岩石粒径分布曲线为约束条件构建数字岩心,是一种比较灵活的数字岩心建模方法.本文根据前人研究结果在过程法的压实过程中用一面状像素层表示颗粒与颗粒间的接触边界,使相邻的骨架颗粒能够相互区分,提高了数字岩心骨架结构的精细程度.在弹性参数模拟中,当变化数字岩心中的颗粒、颗粒接触边界和胶结物的弹性模量时,数字岩心的弹性模量发生相应的变化.数字岩心弹性模量的变化量与组成物质的体积含量及其弹性模量赋值有关.这说明能够利用数字岩石物理研究各种地质和环境因素对岩石弹性性质的影响.进一步分析认为,目前的建模方法对数字岩心的定义还不够充分,利用数字岩石物理解决实际地质问题尚有难度.     

利用三维数字岩心计算龙马溪组页岩等效弹性参数（英文）

弹性参数在甜点区预测和页岩气的开发过程中扮演着重要的角色,因此研究等效弹性参数随页岩气储层属性的变化是一项很有意义的工作。研究中我们用X射线CT扫描技术获得了较为精确的页岩样品微观结构图像。从这些图像中,我们可以获得孔隙度和矿物的详细情况,据此,我们构建了三维数字岩心,并应用有限元法对弹性参数进行了数值模拟,其间深入考察了子样选取、网格划分、求解器类型以及边界条件等,该方法易于区别不同的矿物及其百分含量。本文重点研究孔隙度和干酪根含量对弹性参数的影响,计算结果表明,孔隙度和干酪根含量对弹性性质有较大的影响,当孔隙度和干酪根含量增加时,弹性模量降低,且当孔隙度小于0.75%左右、干酪根含量大于3%左右时弹性参数减小速率较缓。因为孔隙度仅仅为4.5%,孔隙中填充油或气对弹性参数的影响甚微。不同岩心样本具有不同的孔隙度和干酪根含量,传统岩石物理实验不仅昂贵而且费时,而数值模拟是基于数字岩心来计算弹性参数,更加经济、方便。本研究证实了将页岩样品的微观结构图像与弹性模量的计算相结合来预测页岩弹性参数的可行性。     

基于微CT技术的砂岩数字岩石物理实验

数字岩石物理技术可弥补传统岩石物理实验的诸多不足,为岩石物理学研究提供一个新平台.本文以常规砂岩为研究对象,利用微CT扫描结合先进的图像处理技术建立了具有真实孔隙结构特征的三维数字岩芯模型;应用Avizo软件内含的多种形态学算法进行数字岩芯孔隙结构量化及表征研究,统计获取了孔隙度、孔隙体积分布及孔径分布特征,建立了等价孔隙网络模型;将Avizo与多场耦合有限元软件Comsol完美对接,实现了孔隙尺度的渗流模拟并计算获得绝对渗透率,对于考虑固相充填孔隙的情况,模拟计算了岩石有效弹性参数,并与近似Gassmann方程良好验证.本文所提出的将Avizo与Comsol结合使用的方法丰富了现有的数字岩石物理研究手段,为其大规模发展提供了一条新途径.     

Seismic inversion for underground fractures detection based on effective anisotropy and fluid substitution

Underground fractures play an important role in the storage and movement of hydrocarbon fluid. Fracture rock physics has been the useful bridge between fracture parameters and seismic response. In this paper, we aim to use seismic data to predict subsurface fractures based on rock physics. We begin with the construction of fracture rock physics model. Using the model, we may estimate P-wave velocity, S-wave velocity and fracture rock physics parameters. Then we derive a new approximate formula for the analysis of the relationship between fracture rock physics parameters and seismic response, and we also propose the method which uses seismic data to invert the elastic and rock physics parameters of fractured rock. We end with the method verification, which includes using well-logging data to confirm the reliability of fracture rock physics effective model and utilizing real seismic data to validate the applicability of the inversion method. Tests show that the fracture rock physics effective model may be used to estimate velocities and fracture rock physics parameters reliably, and the inversion method is resultful even when the seismic data is added with random noise. Real data test also indicates the inversion method can be applied into the estimation of the elastic and fracture weaknesses parameters in the target area.     

地震尺度下碳酸盐岩储层的岩石物理建模方法(英文)

碳酸盐岩油藏的强非均质性以及孔隙结构的复杂性,使得作为连接油藏参数与地震参数重要桥梁的岩石物理模型,以及作为油藏预测和定量表征最有效工具的流体替换成为岩石物理建模的难点与重点。在碳酸盐岩储层复杂孔隙结构与地震尺度下碳酸盐岩储层非均质性分析基础上,研究采用岩石网格化方法,将地震尺度下非均质碳酸盐岩储层岩石划分为具有独立岩石参数的均质岩石子体,根据岩石孔隙成因与结构特征采用不同岩石物理模型分步计算岩石子块干岩石弹性模量,并根据不同孔隙连通性进行流体替换,计算饱和不同流体岩石弹性模量。基于计算的岩石子块弹性模量,采用Hashin-Shtrikman-Walpole弹性边界计算理论方法实现地震尺度下碳酸盐岩储层弹性参数计算。通过对含有不同类型孔隙组合碳酸盐岩储层模型的弹性模量进行计算与分析,明确不同孔隙对岩石弹性参数的影响特征,模拟分析结果与实际资料认识一致。     

基于方位各向异性弹性阻抗的裂缝岩石物理参数反演方法研究

裂缝储层岩石物理参数的准确获得对地下裂缝预测具有重要意义,而叠前地震反演是获得裂缝岩石物理参数的有效手段.本文从裂缝岩石物理等效模型的构建出发,从测井数据上估测了裂缝岩石物理参数,通过推导含裂缝岩石物理参数的方位各向异性弹性阻抗公式,探讨了基于方位各向异性弹性阻抗的裂缝岩石物理参数地震反演方法.实际工区地震数据应用表明,基于方位各向异性弹性阻抗的裂缝岩石物理参数反演方法合理、可靠,可以降低裂缝岩石物理参数估测的不确定性,为地下裂缝预测提供有力的依据.     

岩石物理模版在储层定量解释中的应用

岩石物理模板采用测井解释的各类岩性矿物骨架点值,选取适合该地区的岩石物理模型,模拟在不同储层组合、不同孔隙及不同饱和度情况下,储层岩石物理参数变化引起的储层测井参数及地球物理响应特征的变化,定量地建立起储层参数同地球物理弹性参数间的解释关系模版.本文根据新场JS42气藏储层参数分析结果,尝试性地将岩石物理解释模板应用于储层定量解释中,对储层高产气区、含水区域进行定量解释,并预测了该气藏的气水界面,该预测结果与实钻井测试情况吻合,证实了该方法的科学性.     

地震岩石物理模型综述

由于石油勘探与开发技术的发展,最近十多年岩石物理学已经成为一门非常有实用价值的学科,它有机地把岩石的物理特性和地震特性联系起来.理论模型是进行岩石物理学研究的主要方法之一.本文归纳总结了基于岩石弹性性质的各种模型,主要包括层状模型、球形孔隙模型、包含体模型和接触模型,以及它们的假设条件、适用范围、局限性和近些年来国内外的应用现状,最后预测了地震岩石物理研究的发展趋势.     

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

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

泥质含量及其分布形式对岩电关系影响的数字岩心仿真(英文)

由于泥质所造成的附加导电现象,泥质含量及其分布形式对电阻率增大系数I和含水饱和度Sw关系具有重要影响,由于岩石物理实验中岩心孔隙结构及其组分构成、分布的微观不可调性,因而泥质分布形式所造成的影响很难通过岩心实验来单独研究。基于数字岩心的格子气自动机方法是一种有效的微观数值模拟方法,本研究利用储层岩心薄片的骨架颗粒尺寸信息资料建立数字岩心模型,结合格子气自动机技术对数字岩心不同饱和流体情况下电的传输特性进行数值模拟研究,揭示了不同泥质含量和泥质分布形式对孔隙介质导电特性非阿尔奇现象产生的影响,建立饱和度指数和泥质含量之间的关系模型,其良好的吻合性表明该方法在岩石物理研究中是一种十分有效的研究方法,而新模型适于在非阿尔奇储层进行准确的饱和度评价。