CN116484766B - 基于ct扫描建立三维非均质等效岩心模型的方法 - Google Patents
基于ct扫描建立三维非均质等效岩心模型的方法 Download PDFInfo
- Publication number
- CN116484766B CN116484766B CN202310540845.9A CN202310540845A CN116484766B CN 116484766 B CN116484766 B CN 116484766B CN 202310540845 A CN202310540845 A CN 202310540845A CN 116484766 B CN116484766 B CN 116484766B
- Authority
- CN
- China
- Prior art keywords
- grid
- core
- equivalent
- pores
- calculating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000002591 computed tomography Methods 0.000 title claims abstract description 33
- 239000011148 porous material Substances 0.000 claims abstract description 102
- 230000035699 permeability Effects 0.000 claims abstract description 28
- 238000006073 displacement reaction Methods 0.000 claims abstract description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 12
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000004088 simulation Methods 0.000 abstract description 27
- 230000007246 mechanism Effects 0.000 abstract description 4
- 238000011084 recovery Methods 0.000 abstract description 4
- 230000006872 improvement Effects 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 abstract description 2
- 230000000007 visual effect Effects 0.000 abstract description 2
- 239000003921 oil Substances 0.000 description 9
- 230000009471 action Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/04—Context-preserving transformations, e.g. by using an importance map
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/08—Indexing scheme for image data processing or generation, in general involving all processing steps from image acquisition to 3D model generation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Fluid Mechanics (AREA)
- Computing Systems (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了基于CT扫描建立三维非均质等效岩心模型的方法，所述方法包括：孔隙定位模块，基于岩心CT扫描数据提取孔隙信息，根据岩心模型网格分辨率将孔隙的体素坐标转换为网格坐标；属性计算模块，利用网格粗化方法计算不同网格对应的孔隙等效半径，基于Kozeny‑carman方程计算孔隙所在网格对应的渗透率；模型生成模块，通过设置无效网格切圆建立三维岩心模型，导入网格渗透率生成三维非均质等效岩心模型；数值模拟模块，将三维非均质等效岩心模型导入CMG‑IMEX模拟器中，利用油水两相模型模拟岩心驱替过程。所述方法可表征岩心内部的驱油动态特征，解决了岩心水驱油动态特征可视化模拟难题，有望推广到岩心尺度的提高采收率机理模拟。
Description
技术领域
本发明属于油藏数值模拟领域，涉及基于CT扫描建立三维非均质等效岩心模型的方法。
背景技术
岩心驱替实验是评价油水流动关系最重要的手段，在石油行业得到了极为广泛的应用。然而岩心驱替过程成本高、周期长且误差大，不利于进一步的分析与研究。数值模拟技术是近几十年随着计算机技术共同发展起来的一种手段，常被用来进行高温高压、微观尺度等物理模拟实验难以开展的研究。通过岩心数值模拟结果拟合物理模拟实验结果证明数值模拟方法的有效性，从而可以将数值模拟方法应用到其他复杂条件以及其他尺度的模拟中。CMG是成熟的商业软件，在建立准确模型的基础上，可以充分模拟实际油藏开发和岩心驱替实验，可进行结果预测和参数敏感性分析。然而，常规的数值模拟所用的等效均质模型无法充分表征非均质岩心内部的驱油动态特征，也无法观测岩心内部微观的、复杂的物理化学现象。基于核磁共振T2谱建立的随机单元等效岩心模型(专利号：202110715803.5)一定程度上考虑了岩心内部的非均质性，但由于其随机性，建立的模型与实际有差异。
岩心CT扫描作为一种无损的可视化手段，近年来在石油工程领域得到了广泛的应用。然而，其主要被用于建立数字岩心模型从而进行孔隙尺度的渗流模拟，在岩心尺度的应用较少。岩心CT扫描可获得无损的岩心三维图像，使得岩心内部微观孔隙结构可视化，通过岩心CT扫描可建立精确的岩心尺度模型。
因此，针对上述问题，本发明提出了基于CT扫描建立非均质岩心模型的方法，为数值模拟提供准确的模型基础，从而精确表征岩心内部的驱油动态特征。
发明内容
本发明的目的在于：提供了基于CT扫描建立三维非均质等效岩心模型的方法，该方法可表征实际岩心内部孔隙结构，能够在模拟岩心水驱油过程中充分考虑岩心内部的非均质性，解决了岩心数值模拟模型与物理模拟实验对象之间的差异，从而表征岩心内部精确的驱油动态特征，有望推广到水驱后岩心尺度的提高采收率机理模拟。
本发明采用的技术方案如下：
基于CT扫描建立三维非均质等效岩心模型的方法，所述方法包括：
孔隙定位模块，基于岩心CT扫描数据提取孔隙信息，根据岩心模型网格分辨率将孔隙的体素坐标转换为网格坐标；
属性计算模块，利用网格粗化方法计算不同网格对应的孔隙等效半径，基于Kozeny-carman方程计算孔隙所在网格对应的渗透率；
模型生成模块，通过设置无效网格切圆建立三维岩心模型，导入渗透率生成三维非均质等效岩心模型；
数值模拟模块，将三维非均质等效岩心模型导入CMG-IMEX模拟器中，利用油水两相模型模拟岩心驱替过程。
进一步地，基于岩心CT扫描数据获取孔隙信息，根据岩心模型网格分辨率将孔隙的体素坐标转换为网格坐标包括以下步骤：
通过AVIZO软件对岩心CT扫描数据进行滤波降噪与二值化处理，在此基础上通过拟合岩心孔隙度确定阈值从而进行阈值分割，显著区分孔隙与骨架，从而进行孔隙分析并导出孔隙的体素坐标以及等效半径；
根据岩心CT扫描图像在三维方向的体素最大值与拟建立岩心模型在三维方向的网格数量，计算二者之间的转换比例；
利用岩心CT扫描获取的岩心内部孔隙的体素坐标数据，按照转换比例将体素坐标转换为网格坐标进行孔隙定位。
进一步地，利用网格粗化方法计算不同网格对应的孔隙等效半径，基于Kozeny-carman方程计算孔隙所在网格对应的渗透率包括以下步骤：
利用岩心CT扫描数据获取岩心内部孔隙体素的等效孔隙半径，利用网格粗化方法计算不同网格对应的等效孔隙半径；
根据孔隙所在的网格坐标，判断各个网格内的孔隙存在情况；
根据不同网格对应的孔隙存在情况，定义不同网格对应的等效孔隙半径；
若网格内部不存在孔隙，则定义网格的等效孔隙半径为0；
若网格内部存在多个孔隙，则通过调和平均方法对不同孔隙半径进行粗化，计算多个孔隙的平均等效孔隙半径作为网格所对应的孔隙等效半径；
根据岩心渗透率、孔隙度、平均等效孔隙半径，利用Kozeny-carman方程计算对应的平均迂曲度；
根据岩心孔隙度、平均迂曲度以及不同网格对应的等效孔隙半径，利用Kozeny-carman方程计算不同网格对应渗透率。
更进一步地，通过设置无效网格切圆方法建立三维岩心模型，导入渗透率生成三维非均质等效岩心模型的步骤包括：
根据两点间距离公式计算网格距离，根据网格间距判断网格是否在内切圆区域中：
式中j、k分别表示网格的y、z坐标，nj、nk分别表示岩心模型在y、z方向网格数。
导出内切圆区域外的网格坐标，并将其定义为无效网格；
设置无效网格，将岩心模型方形端面切成圆形端面，根据实际岩心的直径和长度定义网格尺寸，建立三维岩心模型；
在三维岩心模型中输入网格对应的渗透率，生成三维非均质等效岩心模型。
综上所述，由于采用了上述技术方案，本发明的有益效果是：
基于CT扫描建立三维非均质等效岩心模型的方法，利用目标岩心的CT扫描数据，输出岩心内部孔隙坐标以及等效半径，结合网格粗化原理计算不同网格对应的孔隙等效半径，基于Kozeny-carman方程计算不同网格对应的渗透率，采用三维非均质等效岩心模型精确表征实际岩心。与均质等效岩心模型相比，该方法能够充分考虑岩心内部的非均质性，解决了岩心数值模拟模型与物理模拟实验对象之间的差异，从而精确表征岩心内部的驱油动态特征，解决了岩心水驱油动态特征可视化模拟难题，有望推广到水岩心尺度的提高采收率机理模拟。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图，其中：
图1是基于CT扫描建立三维非均质等效岩心模型的方法的流程图；
图2是本发明实施例一的岩心冠状面CT扫描成像图；
图3是本发明实施例一的岩心冠状面孔隙提取结果示意图；
图4是本发明实施例一的网格切圆方法示意图；
图5是本发明实施例一的三维非均质等效岩心模型示意图；
图6是本发明实施例一的三维非均质等效岩心模型在不同驱替时间的含油饱和度结果图；
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明，即所描述的实施例只是本发明的一部分实施例，而不是全部的实施例。通常在此处描述和附图中示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
需要说明的是，术语“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
基于CT扫描建立三维非均质等效岩心模型的方法，包括以下步骤：
步骤1：基于岩心CT扫描数据提取孔隙信息，根据岩心模型网格分辨率将孔隙的体素坐标转换为网格坐标；
步骤2：利用网格粗化方法计算不同网格对应的孔隙等效半径，基于Kozeny-carman方程计算孔隙所在网格对应的渗透率；
步骤3：通过设置无效网格切圆方法建立三维岩心模型，导入渗透率生成三维非均质等效岩心模型；
步骤4：将三维非均质等效岩心模型导入CMG-IMEX模拟器中，利用油水两相模型模拟岩心驱替过程。
与均质等效岩心模型相比，该方法可表征实际岩心内部孔隙结构，能够在模拟岩心水驱油过程中充分考虑岩心内部的非均质性，解决了岩心数值模拟模型与物理模拟实验对象之间的差异，有望推广到水驱后岩心尺度的提高采收率机理模拟。
下面结合实施例对本发明的特征和性能作进一步的详细描述。
实施例一
本发明的较佳实施例，以致密砂岩水驱油为例，提供了基于CT扫描建立三维非均质等效岩心模型的方法，如图1所示，包括以下步骤：
步骤1：基于岩心CT扫描数据提取孔隙信息，根据岩心模型网格分辨率将孔隙的体素坐标转换为网格坐标进行孔隙定位；
步骤1.1：通过AVIZO软件对岩心CT扫描数据进行处理，获得岩心冠状面CT成像图如图2所示，通过CT图像提取岩心冠状面孔隙结果如图3所示，导出孔隙的体素坐标以及等效半径；
步骤1.2：根据岩心CT扫描数据在三维方向的体素最大值与拟建立岩心模型在三维方向的网格数量(本实例网格数量为30×14×14)，计算二者之间的转换比例，如表1；
表1
步骤1.3：利用导出的孔隙体素坐标数据，按照转换比例将体素坐标转换为网格坐标。
步骤2：利用网格粗化方法计算不同网格对应的孔隙等效半径，基于Kozeny-carman方程计算孔隙所在网格对应的渗透率；
步骤2.1：利用导出的孔隙体素的等效孔隙半径，利用网格粗化方法计算相应网格的等效孔隙半径：
根据孔隙所在的网格坐标，判断各个网格内的孔隙存在情况；
根据不同网格对应的孔隙存在情况，定义不同网格对应的等效孔隙半径；
若网格内部不存在孔隙，则定义网格的等效孔隙半径为0；
若网格内部存在多个孔隙，则通过调和平均方法对不同孔隙半径进行粗化，计算多个孔隙的平均等效孔隙半径作为网格所对应的孔隙等效半径；
步骤2.2：根据岩心渗透率、孔隙度、平均等效孔隙半径，利用Kozeny-carman方程计算对应的平均迂曲度；
平均等效孔隙半径通过调和平均方法计算：
式中表示岩心平均孔隙半径；n表示孔隙数量；ri表示第i个孔隙的孔隙半径。
结合岩心渗透率、孔隙度，利用Kozeny-carman方程计算对应的平均迂曲度，如表2；
式中表示平均迂曲度；/>表示岩心平均渗透率，D；/>表示平均等效孔隙半径，μm；/>表示岩心平均孔隙度。
表2
根据岩心孔隙度、平均迂曲度以及不同网格对应的等效孔隙半径，利用Kozeny-carman方程计算不同网格对应渗透率；
式中K表示网格对应渗透率，D；r表示网格对应的等效孔隙半径，μm。
步骤3：通过设置无效网格切圆方法建立三维岩心模型，导入渗透率生成三维非均质等效岩心模型；
步骤3.1：根据两点间距离公式计算网格距离，利用MATLAB根据网格间距判断网格是否在内切圆区域中：
式中j、k分别表示网格的y、z坐标，nj、nk分别表示岩心模型在y、z方向网格数。
导出内切圆区域外的网格坐标，并将其定义为无效网格；
步骤3.2：设置无效网格，将岩心模型方形端面切成圆形端面，如图4，根据实际岩心的直径和长度定义网格尺寸，建立三维岩心模型。实际岩心尺寸与网格尺寸如表3所示；
表3
步骤3.3：在三维岩心模型中输入网格对应的渗透率，生成三维非均质等效岩心模型，如图5所示；
步骤4：将三维非均质等效岩心模型导入CMG-IMEX模拟器中，利用油水两相模型模拟岩心驱替过程。
本实施例中，三维非均质等效岩心模型在不同驱替时间的含油饱和度结果如图6所示，由图6可知，在水驱油过程中，由于中孔隙和大孔隙中的阻力相对较小，注入水首先驱替中孔隙和大孔隙区域的原油，这也被称为微观指进现象，说明通过该三维非均质等效岩心模型可表征精确的岩心内部驱油动态特征。
需要说明的是，由于说明书附图不得着色和涂改，所以本发明中部分区别明显的地方比较难以显示，若有必要，可提供彩色图片。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明的保护范围，任何熟悉本领域的技术人员在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (2)
1.基于CT扫描建立三维非均质等效岩心模型的方法，其特征在于，所述方法包括：
基于岩心CT扫描数据提取孔隙信息，根据岩心模型网格分辨率将孔隙的体素坐标转换为网格坐标；所述基于岩心CT扫描数据提取孔隙信息，按岩心模型网格分辨率将孔隙的体素坐标转换为网格坐标的步骤包括：根据岩心CT扫描数据在三维方向的体素最大值与拟建立三维非均质等效岩心模型在三维方向的网格数量，计算二者之间的转换比例；利用岩心CT扫描获取的岩心内部孔隙的体素坐标数据，根据体素坐标与网格坐标的转换比例，将体素坐标转换为网格坐标；
利用网格粗化方法计算不同网格对应的孔隙等效半径，基于Kozeny-carman方程计算孔隙所在网格对应的渗透率；所述利用网格粗化方法计算不同网格对应的孔隙等效半径，基于Kozeny-carman方程计算孔隙所在网格对应的渗透率的步骤包括：根据孔隙所在的网格坐标，判断各个网格内的孔隙存在情况；根据不同网格对应的孔隙存在情况，定义不同网格对应的等效孔隙半径，若网格内部不存在孔隙，则定义网格的等效孔隙半径为0，若网格内部存在多个孔隙，通过调和平均方法对不同孔隙半径进行粗化，计算多个孔隙的平均等效孔隙半径作为网格对应的孔隙等效半径；
根据网格对应的等效孔隙半径，利用Kozeny-carman方程计算各网格对应的渗透率；通过设置无效网格切圆建立三维岩心模型，导入渗透率生成三维非均质等效岩心模型；所述通过设置无效网格切圆方法建立三维岩心模型，导入渗透率生成三维非均质等效岩心模型的步骤包括：根据两点间距离公式计算网格距离，根据网格间距判断网格是否在内切圆区域中；导出内切圆区域外的网格坐标，将其定义为无效网格；在三维岩心模型中输入网格对应的渗透率，生成三维非均质等效岩心模型；所述根据网格间距判断网格是否在内切圆区域中的公式为：
式中j、k分别表示网格的y、z坐标，nj、nk分别表示岩心模型在y、z方向网格数；
将三维非均质等效岩心模型导入CMG-IMEX模拟器中，利用油水两相模型模拟岩心驱替过程。
2.根据权利要求1所述的基于CT扫描建立三维非均质等效岩心模型的方法，其特征在于，所述根据网格对应的等效孔隙半径，利用Kozeny-carman方程计算各网格对应的渗透率的步骤包括：
根据岩心渗透率、孔隙度、平均等效孔隙半径，利用Kozeny-carman方程计算对应的平均迂曲度；
根据岩心孔隙度、平均迂曲度以及不同网格对应的等效孔隙半径，利用Kozeny-carman方程计算不同网格对应渗透率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310540845.9A CN116484766B (zh) | 2023-05-15 | 2023-05-15 | 基于ct扫描建立三维非均质等效岩心模型的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310540845.9A CN116484766B (zh) | 2023-05-15 | 2023-05-15 | 基于ct扫描建立三维非均质等效岩心模型的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116484766A CN116484766A (zh) | 2023-07-25 |
| CN116484766B true CN116484766B (zh) | 2024-02-06 |
Family
ID=87225111
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310540845.9A Active CN116484766B (zh) | 2023-05-15 | 2023-05-15 | 基于ct扫描建立三维非均质等效岩心模型的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116484766B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111199582A (zh) * | 2019-12-31 | 2020-05-26 | 中国石油大学(华东) | 一种连续孔隙组分表征的三维数字井筒构建方法 |
| CN111475982A (zh) * | 2020-04-27 | 2020-07-31 | 中国石油大学(华东) | 一种岩石内部磁场梯度加权几何平均值的三因素估算方法 |
| CN111739149A (zh) * | 2020-06-15 | 2020-10-02 | 中国石油大学(华东) | 一种岩石ct扫描图像的油水分布连续性修复方法 |
| CN113051841A (zh) * | 2021-02-10 | 2021-06-29 | 中国石油大学(华东) | 一种强非均质页岩油储层表观渗透率确定方法和装置 |
| CN115688583A (zh) * | 2022-11-02 | 2023-02-03 | 中国石油大学(北京) | 用于岩心的渗流模拟方法、处理器及机器可读存储介质 |
| CN115711842A (zh) * | 2022-10-13 | 2023-02-24 | 陕西理工大学 | 考虑水的瞬态渗流的纳米多孔介质页岩渗透率计算方法 |
| CN116012545A (zh) * | 2023-03-22 | 2023-04-25 | 中国石油大学(华东) | 多尺度数字岩心建模方法、系统、存储介质及应用 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11320357B2 (en) * | 2019-12-23 | 2022-05-03 | Chevron U.S.A. Inc. | System and method for estimation of rock properties from core images |
-
2023
- 2023-05-15 CN CN202310540845.9A patent/CN116484766B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111199582A (zh) * | 2019-12-31 | 2020-05-26 | 中国石油大学(华东) | 一种连续孔隙组分表征的三维数字井筒构建方法 |
| CN111475982A (zh) * | 2020-04-27 | 2020-07-31 | 中国石油大学(华东) | 一种岩石内部磁场梯度加权几何平均值的三因素估算方法 |
| CN111739149A (zh) * | 2020-06-15 | 2020-10-02 | 中国石油大学(华东) | 一种岩石ct扫描图像的油水分布连续性修复方法 |
| CN113051841A (zh) * | 2021-02-10 | 2021-06-29 | 中国石油大学(华东) | 一种强非均质页岩油储层表观渗透率确定方法和装置 |
| CN115711842A (zh) * | 2022-10-13 | 2023-02-24 | 陕西理工大学 | 考虑水的瞬态渗流的纳米多孔介质页岩渗透率计算方法 |
| CN115688583A (zh) * | 2022-11-02 | 2023-02-03 | 中国石油大学(北京) | 用于岩心的渗流模拟方法、处理器及机器可读存储介质 |
| CN116012545A (zh) * | 2023-03-22 | 2023-04-25 | 中国石油大学(华东) | 多尺度数字岩心建模方法、系统、存储介质及应用 |
Non-Patent Citations (6)
| Title |
|---|
| "New technique of True Effective Mobility (TEM-Function) in dynamic rock typing: Reduction of uncertainties in relative permeability data for reservoir simulation";Mirzaei-Paiaman, A 等;《Journal of Petroleum Science and Engineering》;第210-227页 * |
| "基于微观毛细管模型的流动模拟与实验研究";刘鹏飞;《中国优秀硕士学位论文全文数据库（工程科技Ⅰ辑）》;第1-75页 * |
| "基于核磁共振测井的储层渗透率计算方法综述";陈瑶 等;《能源与环保》;第44卷(第01期);第125-131页 * |
| "基于格子Boltzmann方法的多孔介质微观渗吸机理研究";刘洋;《中国博士学位论文全文数据库 (工程科技Ⅰ辑)》;第1-119页 * |
| "电磁场-渗流场耦合作用下离子液体多孔介质流动模型";贾虎 等;《力学学报》;第53卷(第2期);第2214-2224页 * |
| "缝洞型储层渗透率和饱和度计算及约束反演方法研究";魏帅帅;《中国博士学位论文全文数据库 (工程科技Ⅰ辑)》;第1-103页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116484766A (zh) | 2023-07-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107817199B (zh) | 一种致密砂岩多尺度孔隙模型的构建方法和应用 | |
| CN101556703B (zh) | 基于连续切片图像的网络模型建立方法 | |
| CN111428321B (zh) | 一种基于简化数字岩心的砾岩储层孔隙网络模型建模方法 | |
| CN109242985B (zh) | 一种从三维图像确定孔隙结构关键参数的方法 | |
| JP5588713B2 (ja) | シミュレーション装置及びそのコンピュータプログラム | |
| WO2014104909A1 (ru) | Способ получения характеристической трехмерной модели образца пористого материала для исследования свойств проницаемости | |
| CN106023153A (zh) | 一种测量水体中气泡的方法 | |
| CN110222368B (zh) | 一种利用二维切片计算岩心三维孔隙度和渗透率的方法 | |
| CN112861890B (zh) | 一种储层评价模型构建方法及储层识别方法 | |
| CN114387332B (zh) | 一种管道测厚方法及装置 | |
| CN116484766B (zh) | 基于ct扫描建立三维非均质等效岩心模型的方法 | |
| Li et al. | Description of gas hydrate using digital core technology | |
| CN113029899B (zh) | 一种基于显微图像处理的砂岩渗透率计算方法 | |
| CN112903555B (zh) | 考虑孔隙各向异性的多孔介质渗透率计算方法及装置 | |
| CN116342815B (zh) | 一种页岩孔隙空间的多尺度孔隙网络模型构建方法 | |
| CN114705606A (zh) | 一种基于网络化分析的岩石内部关键渗流节点的封堵方法 | |
| CN116413183A (zh) | 一种孔隙连通性评价方法和装置 | |
| Bemis et al. | Case study: a methodology for plume visualization with application to real-time acquisition and navigation | |
| CN110441204B (zh) | 一种基于数字岩心模拟的致密储层压裂液伤害数字化评价方法 | |
| MacLeod et al. | Ford's Gamma‐Gamma Village Simulation Revisited: Highlighting the Need for a New Middle‐Range Theory of Archaeological Types 1 | |
| CN105317433B (zh) | 基于level-set函数的颗粒堆微观孔道提取方法 | |
| CN111612045A (zh) | 一种获取目标检测数据集的通用方法 | |
| CN117036635B (zh) | 一种基于图像纹理分类的页岩多尺度数字岩心构建方法 | |
| Banowski et al. | Experiments on vertical gas-liquid pipe flows using ultrafast X-ray tomography | |
| CN114092489B (zh) | 多孔介质渗流通道提取、模型训练方法、装置及设备 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |