基于CT扫描建立三维非均质等效岩心模型的方法技术

本发明专利技术公开了基于CT扫描建立三维非均质等效岩心模型的方法，所述方法包括：孔隙定位模块，基于岩心CT扫描数据提取孔隙信息，根据岩心模型网格分辨率将孔隙的体素坐标转换为网格坐标；属性计算模块，利用网格粗化方法计算不同网格对应的孔隙等效半径，基于Kozeny

【技术实现步骤摘要】
基于CT扫描建立三维非均质等效岩心模型的方法


[0001]本专利技术属于油藏数值模拟领域，涉及基于CT扫描建立三维非均质等效岩心模型的方法。

技术介绍

[0002]岩心驱替实验是评价油水流动关系最重要的手段，在石油行业得到了极为广泛的应用。然而岩心驱替过程成本高、周期长且误差大，不利于进一步的分析与研究。数值模拟技术是近几十年随着计算机技术共同发展起来的一种手段，常被用来进行高温高压、微观尺度等物理模拟实验难以开展的研究。通过岩心数值模拟结果拟合物理模拟实验结果证明数值模拟方法的有效性，从而可以将数值模拟方法应用到其他复杂条件以及其他尺度的模拟中。CMG是成熟的商业软件，在建立准确模型的基础上，可以充分模拟实际油藏开发和岩心驱替实验，可进行结果预测和参数敏感性分析。然而，常规的数值模拟所用的等效均质模型无法充分表征非均质岩心内部的驱油动态特征，也无法观测岩心内部微观的、复杂的物理化学现象。基于核磁共振T2谱建立的随机单元等效岩心模型(专利号：202110715803.5)一定程度上考虑了岩心内部的非均质性，但由于其随机性，建立的模型与实际有差异。
[0003]岩心CT扫描作为一种无损的可视化手段，近年来在石油工程领域得到了广泛的应用。然而，其主要被用于建立数字岩心模型从而进行孔隙尺度的渗流模拟，在岩心尺度的应用较少。岩心CT扫描可获得无损的岩心三维图像，使得岩心内部微观孔隙结构可视化，通过岩心CT扫描可建立精确的岩心尺度模型。
[0004]因此，针对上述问题，本专利技术提出了基于CT扫描建立非均质岩心模型的方法，为数值模拟提供准确的模型基础，从而精确表征岩心内部的驱油动态特征。

技术实现思路

[0005]本专利技术的目的在于：提供了基于CT扫描建立三维非均质等效岩心模型的方法，该方法可表征实际岩心内部孔隙结构，能够在模拟岩心水驱油过程中充分考虑岩心内部的非均质性，解决了岩心数值模拟模型与物理模拟实验对象之间的差异，从而表征岩心内部精确的驱油动态特征，有望推广到水驱后岩心尺度的提高采收率机理模拟。
[0006]本专利技术采用的技术方案如下：
[0007]基于CT扫描建立三维非均质等效岩心模型的方法，所述方法包括：
[0008]孔隙定位模块，基于岩心CT扫描数据提取孔隙信息，根据岩心模型网格分辨率将孔隙的体素坐标转换为网格坐标；
[0009]属性计算模块，利用网格粗化方法计算不同网格对应的孔隙等效半径，基于Kozeny
‑
carman方程计算孔隙所在网格对应的渗透率；
[0010]模型生成模块，通过设置无效网格切圆建立三维岩心模型，导入渗透率生成三维非均质等效岩心模型；
[0011]数值模拟模块，将三维非均质等效岩心模型导入CMG
‑
IMEX模拟器中，利用油水两
相模型模拟岩心驱替过程。
[0012]进一步地，基于岩心CT扫描数据获取孔隙信息，根据岩心模型网格分辨率将孔隙的体素坐标转换为网格坐标包括以下步骤：
[0013]通过AVIZO软件对岩心CT扫描数据进行滤波降噪与二值化处理，在此基础上通过拟合岩心孔隙度确定阈值从而进行阈值分割，显著区分孔隙与骨架，从而进行孔隙分析并导出孔隙的体素坐标以及等效半径；
[0014]根据岩心CT扫描图像在三维方向的体素最大值与拟建立岩心模型在三维方向的网格数量，计算二者之间的转换比例；
[0015]利用岩心CT扫描获取的岩心内部孔隙的体素坐标数据，按照转换比例将体素坐标转换为网格坐标进行孔隙定位。
[0016]进一步地，利用网格粗化方法计算不同网格对应的孔隙等效半径，基于Kozeny
‑
carman方程计算孔隙所在网格对应的渗透率包括以下步骤：
[0017]利用岩心CT扫描数据获取岩心内部孔隙体素的等效孔隙半径，利用网格粗化方法计算不同网格对应的等效孔隙半径；
[0018]根据孔隙所在的网格坐标，判断各个网格内的孔隙存在情况；
[0019]根据不同网格对应的孔隙存在情况，定义不同网格对应的等效孔隙半径；
[0020]若网格内部不存在孔隙，则定义网格的等效孔隙半径为0；
[0021]若网格内部存在多个孔隙，则通过调和平均方法对不同孔隙半径进行粗化，计算多个孔隙的平均等效孔隙半径作为网格所对应的孔隙等效半径；
[0022]根据岩心渗透率、孔隙度、平均等效孔隙半径，利用Kozeny
‑
carman方程计算对应的平均迂曲度；
[0023]根据岩心孔隙度、平均迂曲度以及不同网格对应的等效孔隙半径，利用Kozeny
‑
carman方程计算不同网格对应渗透率。
[0024]更进一步地，通过设置无效网格切圆方法建立三维岩心模型，导入渗透率生成三维非均质等效岩心模型的步骤包括：
[0025]根据两点间距离公式计算网格距离，根据网格间距判断网格是否在内切圆区域中：
[0026][0027]式中j、k分别表示网格的y、z坐标，nj、nk分别表示岩心模型在y、z方向网格数。
[0028]导出内切圆区域外的网格坐标，并将其定义为无效网格；
[0029]设置无效网格，将岩心模型方形端面切成圆形端面，根据实际岩心的直径和长度定义网格尺寸，建立三维岩心模型；
[0030]在三维岩心模型中输入网格对应的渗透率，生成三维非均质等效岩心模型。
[0031]综上所述，由于采用了上述技术方案，本专利技术的有益效果是：
[0032]基于CT扫描建立三维非均质等效岩心模型的方法，利用目标岩心的CT扫描数据，输出岩心内部孔隙坐标以及等效半径，结合网格粗化原理计算不同网格对应的孔隙等效半径，基于Kozeny
‑
carman方程计算不同网格对应的渗透率，采用三维非均质等效岩心模型精确表征实际岩心。与均质等效岩心模型相比，该方法能够充分考虑岩心内部的非均质性，解
决了岩心数值模拟模型与物理模拟实验对象之间的差异，从而精确表征岩心内部的驱油动态特征，解决了岩心水驱油动态特征可视化模拟难题，有望推广到水岩心尺度的提高采收率机理模拟。
附图说明
[0033]为了更清楚地说明本专利技术实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本专利技术的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图，其中：
[0034]图1是基于CT扫描建立三维非均质等效岩心模型的方法的流程图；
[0035]图2是本专利技术实施例一的岩心冠状面CT扫描成像图；
[0036]图3是本专利技术实施例一的岩心冠状面孔隙提取结果示意图；
[0037]图4是本专利技术实施例一的网格切圆方法示意图；
[0038]图5是本专利技术实施例一的三维非均质等效岩心模型示意图；
[0039]图6是本专利技术实施例一的三维非均本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于CT扫描建立三维非均质等效岩心模型的方法，其特征在于，所述方法包括：孔隙定位模块，基于岩心CT扫描数据提取孔隙信息，根据岩心模型网格分辨率将孔隙的体素坐标转换为网格坐标；属性计算模块，利用网格粗化方法计算不同网格对应的孔隙等效半径，基于Kozeny
‑
carman方程计算孔隙所在网格对应的渗透率；模型生成模块，通过设置无效网格切圆建立三维岩心模型，导入渗透率生成三维非均质等效岩心模型；数值模拟模块，将三维非均质等效岩心模型导入CMG
‑
IMEX模拟器中，利用油水两相模型模拟岩心驱替过程。2.根据权利要求1所述的基于CT扫描建立三维非均质等效岩心模型的方法，其特征在于，所述基于岩心CT扫描数据提取孔隙信息，按岩心模型网格分辨率将孔隙的体素坐标转换为网格坐标的步骤包括：根据岩心CT扫描数据在三维方向的体素最大值与拟建立三维非均质等效岩心模型在三维方向的网格数量，计算二者之间的转换比例；利用岩心CT扫描获取的岩心内部孔隙的体素坐标数据，根据体素坐标与网格坐标的转换比例，将体素坐标转换为网格坐标。3.根据权利要求1所述的基于CT扫描建立三维非均质等效岩心模型的方法，其特征在于，所述利用网格粗化方法计算不同网格对应的孔隙等效半径，基于Kozeny
‑
carman方程计算孔隙所在网格对应的渗透率的步骤包括：根据孔隙所在的网格坐标，判断各个网格内的孔隙存在情况；根据不同网格对应的孔隙存在情况，定义不同网格对应的等效孔隙半径；根据网格对应的等效孔隙半径，利用Kozeny
‑
carman方程计算各网格对应的渗透率。4.根据权利要求1所述的基于CT扫描建立三维非均质等效岩心模...

【专利技术属性】
技术研发人员：贾虎，周子力，代杨，张瑞，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：