CN110222368A - 一种利用二维切片计算岩心三维孔隙度和渗透率的方法 - Google Patents

一种利用二维切片计算岩心三维孔隙度和渗透率的方法 Download PDF

Info

Publication number
CN110222368A
CN110222368A CN201910364713.9A CN201910364713A CN110222368A CN 110222368 A CN110222368 A CN 110222368A CN 201910364713 A CN201910364713 A CN 201910364713A CN 110222368 A CN110222368 A CN 110222368A
Authority
CN
China
Prior art keywords
rock core
permeability
porosity
dimensional
dimension slicing
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.)
Granted
Application number
CN201910364713.9A
Other languages
English (en)
Other versions
CN110222368B (zh
Inventor
冯其红
赵蕴昌
王森
察鲁明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201910364713.9A priority Critical patent/CN110222368B/zh
Publication of CN110222368A publication Critical patent/CN110222368A/zh
Application granted granted Critical
Publication of CN110222368B publication Critical patent/CN110222368B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种利用二维切片计算岩心三维孔隙度和渗透率的方法，该方法包括：获取岩心二维切片的二值化图像；计算岩心每个二维切片二值化图像的孔隙度，并对得到的岩心二维切片二值化图像的孔隙度进行算术平均，得到岩心的平均二维孔隙度；计算岩心每个二维切片二值化图像的渗透率，并对得到的岩心二维切片二值化图像的渗透率进行算术平均，得到岩心的平均二维渗透率；由岩心的平均二维孔隙度和平均二维渗透率，利用孔隙度和渗透率的二维‑三维转换关系，计算得到岩心的三维孔隙度和三维渗透率。本发明将求取岩心渗透率时的三维数值模拟问题转化为多个二维数值模拟问题，减小了计算量，缩短了模拟时间，提高了模拟效率和收敛性。

Description

一种利用二维切片计算岩心三维孔隙度和渗透率的方法
技术领域
本发明属于油气田开发技术领域，涉及一种利用二维切片计算岩心三维孔隙度和渗透率的方法。
背景技术
孔隙度和渗透率表征了多孔介质储集空间的大小和输运能力的强弱，是油气田开发过程中重要的基础参数。准确获取岩心的孔隙度和渗透率对于油气资源量的估算和产能评价具有极其重要的作用。
除物模实验外，目前油气田开发领域常通过数值模拟方法计算岩心的三维孔隙度和渗透率。但发明人发现现有数值模拟方法直接计算岩心的三维孔隙度和渗透率时计算量过大，模拟时间过长，且易于出现不收敛的情形。
发明内容
本发明的目的在于提供一种利用二维切片计算岩心三维孔隙度和渗透率的方法。本发明能够利用岩心的二维切片图像快速准确的模拟得到岩心的二维孔隙度和二维渗透率，进而利用本发明提供的二维-三维转换关系快速计算岩心的三维孔隙度和三维渗透率。本发明解决了直接利用三维数值模拟求取岩心渗透率时模拟时间过长、收敛性差的问题。本发明的有益效果是简单、准确、高效的计算岩心的三维孔隙度和渗透率。
本发明一种利用二维切片计算岩心三维孔隙度和渗透率的方法，包括以下步骤：
步骤a：获取岩心二维切片的二值化图像；
步骤b：计算岩心每个二维切片二值化图像的孔隙度，并对得到的岩心二维切片二值化图像的孔隙度进行算术平均，得到岩心的平均二维孔隙度；
步骤c：计算岩心每个二维切片二值化图像的渗透率，并对得到的岩心二维切片二值化图像的渗透率进行算术平均，得到岩心的平均二维渗透率；
步骤d：由岩心的平均二维孔隙度和平均二维渗透率，利用孔隙度和渗透率的二维-三维转换关系，计算得到岩心的三维孔隙度和三维渗透率。
进一步，获取岩心二维切片的二值化图像，包括：
利用X射线CT对所述岩心进行扫描，得到所述岩心的三维多孔介质模型；
将所述岩心的三维多孔介质模型导入AutoCAD软件，利用AutoCAD软件中的生成二维截面功能以等间距方式对所述岩心的三维多孔介质模型进行切片，得到至少60张所述岩心的二维切片图像；
采用OTSU算法对所述岩心的二维切片图像进行二值化处理，将图像上孔隙所对应像素的灰度值设置为0，岩石骨架所对应像素的灰度值设置为255，得到所述岩心二维切片的二值化图像。
进一步，计算岩心每个二维切片二值化图像的孔隙度，并对得到的岩心二维切片二值化图像的孔隙度进行算术平均，得到岩心的平均二维孔隙度，包括：
针对所述岩心二维切片的二值化图像，计算每个二维切片图像的孔隙度；
二维切片图像的孔隙度φi的计算公式如下：
式中：φi为第i个二维切片图像的孔隙度；Ni,孔为第i个二维切片图像上孔隙所对应像素的个数；Ni为第i个二维切片图像所对应像素的总数；
对岩心每个二维切片图像的孔隙度进行算术平均，得到所述岩心的平均二维孔隙度
进一步，计算岩心每个二维切片二值化图像的渗透率，并对得到的岩心二维切片二值化图像的渗透率进行算术平均，得到岩心的平均二维渗透率，包括：
利用AutoCAD软件将所述岩心二维切片的二值化图像转换为DXF文件，得到岩心二维切片的孔隙结构模型；
将所述岩心二维切片的孔隙结构模型导入计算流体力学模拟软件中；
在计算流体力学模拟软件中将所述岩心二维切片的孔隙结构模型进行网格划分；
在计算流体力学模拟软件中将所述岩心二维切片的孔隙结构模型的一侧定义为入口端，另一侧定义为出口端；
在计算流体力学模拟软件中设置所述岩心二维切片的孔隙结构模型孔隙内流体为水；
在计算流体力学模拟软件中设置所述岩心二维切片的孔隙结构模型入口端压力为P1，出口端压力为0，在保持入口端和出口端压力不变的情况下利用相场方法、有限体积法或格子Boltzmann方法进行流动模拟，得到定压差条件下通过岩心二维切片的流体流量Q；
在本发明中，较优地，所述P1为0.1-10kPa；
根据下式计算岩心每个二维切片图像的渗透率ki：
式中：ki为第i个二维切片图像的渗透率，μm2；Q为通过第i个二维切片图像的流量，cm3/s；μ为流体的粘度，mPa·s；L为岩心长度，cm；A为第i个二维切片图像的横截面积，cm2；Δp为所述岩心二维切片的孔隙结构模型入口端和出口端压力之差，10-1MPa。
对岩心每个二维切片图像的渗透率进行算术平均，得到所述岩心的平均二维渗透率。
进一步，由岩心的平均二维孔隙度和平均二维渗透率，利用孔隙度和渗透率的二维-三维转换关系，计算得到岩心的三维孔隙度和三维渗透率，包括：
将所述岩心的平均二维孔隙度代入孔隙度的二维-三维转换关系，计算得到所述岩心的三维孔隙度φ；
孔隙度的二维-三维转换关系为：
式中：φ为岩心的三维孔隙度；为岩心的平均二维孔隙度；
将所述岩心的平均二维渗透率代入渗透率的二维-三维转换关系，计算得到所述岩心的三维渗透率k；
渗透率的二维-三维转换关系为：
式中：k为岩心的三维渗透率；为岩心的平均二维渗透率。
附图说明
图1为本发明所述的一种利用二维切片计算岩心三维孔隙度和渗透率的方法的流程示意图；
图2为本发明实施例中岩心的三维多孔介质模型示意图；
图3为本发明实施例中对岩心的三维多孔介质模型进行切片过程的示意图；
图4为本发明实施例中是岩心的二维切片图像的示意图；
图5为本发明实施例中岩心各个二维切片图像的孔隙度与平均二维孔隙度示意图；
图6为本发明实施例中模拟计算岩心二维切片图像渗透率的示意图；
图7为本发明实施例中岩心各个二维切片图像的渗透率与平均二维渗透率示意图。
具体实施方式
下面结合附图和具体实施例对本发明的实施方式进行详细说明。
本发明提供了一种利用二维切片计算岩心三维孔隙度和渗透率的方法，如图1所示，包括以下步骤：
步骤a：获取岩心二维切片的二值化图像；
步骤b：计算岩心每个二维切片二值化图像的孔隙度，并对得到的岩心二维切片二值化图像的孔隙度进行算术平均，得到岩心的平均二维孔隙度；
步骤c：计算岩心每个二维切片二值化图像的渗透率，并对得到的岩心二维切片二值化图像的渗透率进行算术平均，得到岩心的平均二维渗透率；
步骤d：由岩心的平均二维孔隙度和平均二维渗透率，利用孔隙度和渗透率的二维-三维转换关系，计算得到岩心的三维孔隙度和三维渗透率。
由附图1所示流程可以得知，在本发明实施例中，通过获取岩心二维切片的二值化图像，分别计算二维切片图像的孔隙度算术平均值和渗透率算术平均值，并结合孔隙度和渗透率的二维-三维转换关系，计算得到岩心的三维孔隙度和渗透率。可见，该发明将求取岩心渗透率时的三维数值模拟问题转化为多个二维数值模拟问题，减小了计算量，缩短了模拟时间，提高了模拟效率和收敛性。
在本发明步骤a中，所述获取岩心二维切片的二值化图像，包括：
利用X射线CT对所述岩心进行扫描，得到所述岩心的三维多孔介质模型；
将所述岩心的三维多孔介质模型导入AutoCAD软件，利用AutoCAD软件中的生成二维截面功能以等间距方式对所述岩心的三维多孔介质模型进行切片，得到至少60张所述岩心的二维切片图像；
采用OTSU算法对所述岩心的二维切片图像进行二值化处理，将图像上孔隙所对应像素的灰度值设置为0，岩石骨架所对应像素的灰度值设置为255，得到所述岩心二维切片的二值化图像。
在本发明步骤b中，所述计算岩心每个二维切片图像的孔隙度，并对得到的岩心二维切片图像的孔隙度进行算术平均，得到岩心的平均二维孔隙度，包括：
针对所述岩心二维切片的二值化图像，计算每个二维切片图像的孔隙度；
二维切片图像的孔隙度φi的计算公式如下：
式中：φi为第i个二维切片图像的孔隙度；Ni,孔为第i个二维切片图像上孔隙所对应像素的个数；Ni为第i个二维切片图像所对应像素的总数；
对岩心每个二维切片图像的孔隙度进行算术平均，得到所述岩心的平均二维孔隙度
在本发明步骤c中，计算岩心每个二维切片图像的渗透率，并对得到的岩心二维切片图像的渗透率进行算术平均，得到岩心的平均二维渗透率，包括：
利用AutoCAD软件将所述岩心二维切片的二值化图像转换为DXF文件，得到岩心二维切片的孔隙结构模型；
将所述岩心二维切片的孔隙结构模型导入计算流体力学模拟软件中；
在计算流体力学模拟软件中将所述岩心二维切片的孔隙结构模型进行网格划分；
在计算流体力学模拟软件中将所述岩心二维切片的孔隙结构模型的一侧定义为入口端，另一侧定义为出口端；
在计算流体力学模拟软件中设置所述岩心二维切片的孔隙结构模型孔隙内流体为水；
在计算流体力学模拟软件中设置所述岩心二维切片的孔隙结构模型入口端压力为P1，出口端压力为0，在保持入口端和出口端压力不变的情况下利用相场方法、有限体积法或格子Boltzmann方法进行流动模拟，得到定压差条件下通过岩心二维切片的流体流量Q；
在本发明中，较优地，所述P1为0.1-10kPa；
根据下式计算岩心每个二维切片图像的渗透率ki：
式中：ki为第i个二维切片图像的渗透率，μm2；Q为通过第i个二维切片图像的流量，cm3/s；μ为流体的粘度，mPa·s；L为岩心长度，cm；A为第i个二维切片图像的横截面积，cm2；Δp为所述岩心二维切片的孔隙结构模型入口端和出口端压力之差，10-1MPa。
对岩心每个二维切片图像的渗透率进行算术平均，得到所述岩心的平均二维渗透率。
在本发明步骤d中，由岩心的平均二维孔隙度和平均二维渗透率，利用孔隙度和渗透率的二维-三维转换关系，计算得到岩心的三维孔隙度和三维渗透率，包括：
将所述岩心的平均二维孔隙度代入孔隙度的二维-三维转换关系，计算得到所述岩心的三维孔隙度φ；
孔隙度的二维-三维转换关系为：
式中：φ为岩心的三维孔隙度；为岩心的平均二维孔隙度；
将所述岩心的平均二维渗透率代入渗透率的二维-三维转换关系，计算得到所述岩心的三维渗透率k；
渗透率的二维-三维转换关系为：
式中：k为岩心的三维渗透率；为岩心的平均二维渗透率。
孔隙度和渗透率的二维-三维转换关系是采用以下方法确定的：采集油田现场不同区块、不同层系、不同岩性的大量岩心，利用X射线CT扫描得到岩心的三维多孔介质模型，并对其进行切片处理计算岩心的平均二维孔隙度和平均二维渗透率。同时采用物理模拟实验方法测得岩心的三维孔隙度和三维渗透率，然后对结果进行统计分析，对岩心的三维孔隙度与平均二维孔隙度之间的关系曲线进行回归，对岩心的三维渗透率与平均二维渗透率之间的关系曲线进行回归，由此分别得到孔隙度和渗透率的二维-三维转换关系。
下面以一个具体实施例对本发明进行说明：
步骤a：获取岩心二维切片的二值化图像；
利用X射线CT对所述岩心进行扫描，得到所述岩心的三维多孔介质模型。图2为所述岩心的三维多孔介质模型，其中黑色球体为岩石颗粒，白色为孔隙，模型在三个方向的长度均为100μm。
将所述岩心的三维多孔介质模型导入AutoCAD软件，利用AutoCAD软件中的生成二维截面功能沿垂直方向等间距的对所述岩心的三维多孔介质模型进行切片(参见附图3)，得到80张所述岩心的二维切片图像；
采用OTSU算法对所述岩心的二维切片图像进行二值化处理，将图像上孔隙所对应像素的灰度值设置为0，岩石骨架所对应像素的灰度值设置为255，得到所述岩心二维切片的二值化图像。图4为所述岩心二维切片的二值化图像示意图，其中白色为岩石颗粒，灰色为孔隙。
步骤b：计算岩心每个二维切片图像的孔隙度，并对得到的岩心二维切片图像的孔隙度进行算术平均，得到岩心的平均二维孔隙度；
以图4中所述岩心二维切片的二值化图像为例，计算二维切片图像的孔隙度。该二维切片模型的像素总数为Ni＝100×100＝10000，其中孔隙所对应像素的个数Ni,孔＝2763，利用下式计算得到该二维切片图像的孔隙度为：
用该方法分别计算岩心每个二维切片图像的孔隙度，并对这80个二维切片图像的孔隙度进行算术平均，得到所述岩心的平均二维孔隙度(参见附图5)。
步骤c：计算岩心每个二维切片图像的渗透率，并对得到的岩心二维切片图像的渗透率进行算术平均，得到岩心的平均二维渗透率；
以图4中所述岩心二维切片的二值化图像为例，利用AutoCAD软件将所述岩心二维切片的二值化图像转换为DXF文件，得到岩心二维切片的孔隙结构模型；
将所述岩心二维切片的孔隙结构模型导入计算流体力学模拟软件中；
在计算流体力学模拟软件中将所述岩心二维切片的孔隙结构模型进行网格划分；
在计算流体力学模拟软件中将所述岩心二维切片的孔隙结构模型的左侧定义为入口端，右侧定义为出口端；
在计算流体力学模拟软件中设置所述岩心二维切片的孔隙结构模型孔隙内流体为水，密度1000kg/m3,粘度为1mPa·s；
在计算流体力学模拟软件中设置所述岩心二维切片的孔隙结构模型入口端压力为3kPa，出口端压力为0，在保持入口端和出口端压力不变的情况下利用相场方法进行流动模拟(参见附图6)，得到定压差条件下通过该二维切片的流体流量Q＝0.0243cm3/s；
利用下式计算岩心二维切片图像的渗透率ki：
用该方法分别计算岩心每个二维切片图像的渗透率，并对这80个二维切片图像的渗透率进行算术平均，得到所述岩心的平均二维渗透率如图7所示。
步骤d：对岩心的平均二维孔隙度和平均二维渗透率进行分析计算，得到岩心的三维孔隙度和三维渗透率。
所述岩心的平均二维孔隙度为0.31，根据孔隙度的二维-三维转换关系，计算得到所述岩心的三维孔隙度为：
式中：φ为岩心的三维孔隙度；为岩心的平均二维孔隙度。因此所述岩心的三维孔隙度为0.27。
所述岩心的平均二维渗透率为0.56μm2，根据渗透率的二维-三维转换关系，计算得到所述岩心的三维渗透率为：
式中：k为岩心的三维渗透率；为岩心的平均二维渗透率。因此所述岩心的三维渗透率为1.98μm2。
利用三维数值模拟方法计算得到所述岩心的渗透率为1.96μm2，模拟时间为24小时，然而在计算机设备相同的情况下，利用本发明所提供的方法计算所述岩心的渗透率时所需时间仅为1.0小时，计算效率提高了24倍。此外，由物理模拟实验方法测得该岩心的渗透率为1.981μm2，因此本发明的模拟精度更高。对于其它的实施例，三维数值模拟方法可能出现计算不收敛的情形，然而本发明所提供的方法均能得到收敛的结果，因此稳定性更好。
由此可见，本发明将求取岩心渗透率时的三维数值模拟问题转化为多个二维数值模拟问题，减小了计算量，缩短了模拟时间，提高了模拟效率和收敛性。
以上所述仅是对本发明的较佳实施方式而已，并非对本发明作任何形式上的限制，凡是依据本发明的技术实质对以上实施方式所做的任何简单修改，等同变化与修饰，均属于本发明技术方案的范围内。

Claims (6)

1.一种利用二维切片计算岩心三维孔隙度和渗透率的方法，其特征在于，包括：
步骤a：获取岩心二维切片的二值化图像；
步骤b：计算岩心每个二维切片二值化图像的孔隙度，并对得到的岩心二维切片二值化图像的孔隙度进行算术平均，得到岩心的平均二维孔隙度；
步骤c：计算岩心每个二维切片二值化图像的渗透率，并对得到的岩心二维切片二值化图像的渗透率进行算术平均，得到岩心的平均二维渗透率；
步骤d：由岩心的平均二维孔隙度和平均二维渗透率，利用孔隙度和渗透率的二维-三维转换关系，计算得到岩心的三维孔隙度和三维渗透率。
2.根据权利要求1所述一种利用二维切片计算岩心三维孔隙度和渗透率的方法，其特征在于，所述步骤a中获取岩心二维切片的二值化图像，包括：
利用X射线CT对所述岩心进行扫描，得到所述岩心的三维多孔介质模型；
将所述岩心的三维多孔介质模型导入AutoCAD软件，利用AutoCAD软件中的生成二维截面功能以等间距方式对所述岩心的三维多孔介质模型进行切片，得到至少60张所述岩心的二维切片图像；
采用OTSU算法对所述岩心的二维切片图像进行二值化处理，将图像上孔隙所对应像素的灰度值设置为0，岩石骨架所对应像素的灰度值设置为255，得到所述岩心二维切片的二值化图像。
3.根据权利要求1所述一种利用二维切片计算岩心三维孔隙度和渗透率的方法，其特征在于，所述步骤b中计算岩心每个二维切片二值化图像的孔隙度，并对得到的岩心二维切片二值化图像的孔隙度进行算术平均，得到岩心的平均二维孔隙度，包括：
针对所述岩心二维切片的二值化图像，计算每个二维切片图像的孔隙度；
二维切片图像的孔隙度φi的计算公式如下：
式中：φi为第i个二维切片图像的孔隙度；Ni,孔为第i个二维切片图像上孔隙所对应像素的个数；Ni为第i个二维切片图像所对应像素的总数；
对岩心每个二维切片图像的孔隙度进行算术平均，得到所述岩心的平均二维孔隙度
4.根据权利要求1所述一种利用二维切片计算岩心三维孔隙度和渗透率的方法，其特征在于，所述步骤c中计算岩心每个二维切片二值化图像的渗透率，并对得到的岩心二维切片二值化图像的渗透率进行算术平均，得到岩心的平均二维渗透率，包括：
利用AutoCAD软件将所述岩心二维切片的二值化图像转换为DXF文件，得到岩心二维切片的孔隙结构模型；
将所述岩心二维切片的孔隙结构模型导入计算流体力学模拟软件中；
在计算流体力学模拟软件中将所述岩心二维切片的孔隙结构模型进行网格划分；
在计算流体力学模拟软件中将所述岩心二维切片的孔隙结构模型的一侧定义为入口端，另一侧定义为出口端；
在计算流体力学模拟软件中设置所述岩心二维切片的孔隙结构模型孔隙内流体为水；
在计算流体力学模拟软件中设置所述岩心二维切片的孔隙结构模型入口端压力为P1，出口端压力为0，在保持入口端和出口端压力不变的情况下利用相场方法、有限体积法或格子Boltzmann方法进行流动模拟，得到定压差条件下通过岩心二维切片的流体流量Q；
根据下式计算岩心每个二维切片图像的渗透率ki：
式中：ki为第i个二维切片图像的渗透率，μm2；Q为通过第i个二维切片图像的流量，cm3/s；μ为流体的粘度，mPa·s；L为岩心长度，cm；A为第i个二维切片图像的横截面积，cm2；Δp为所述岩心二维切片的孔隙结构模型入口端和出口端压力之差，10-1MPa。
对岩心每个二维切片图像的渗透率进行算术平均，得到所述岩心的平均二维渗透率。
5.根据权利要求1所述一种利用二维切片计算岩心三维孔隙度和渗透率的方法，其特征在于，所述步骤d中由岩心的平均二维孔隙度和平均二维渗透率，利用孔隙度和渗透率的二维-三维转换关系，计算得到岩心的三维孔隙度和三维渗透率，包括：
将所述岩心的平均二维孔隙度代入孔隙度的二维-三维转换关系，计算得到所述岩心的三维孔隙度φ；
孔隙度的二维-三维转换关系为：
式中：φ为岩心的三维孔隙度；为岩心的平均二维孔隙度；
将所述岩心的平均二维渗透率代入渗透率的二维-三维转换关系，计算得到所述岩心的三维渗透率k；
渗透率的二维-三维转换关系为：
式中：k为岩心的三维渗透率；为岩心的平均二维渗透率。
6.根据权利要求1、4所述一种利用二维切片计算岩心三维孔隙度和渗透率的方法，其特征在于，所述岩心二维切片的孔隙结构模型入口端压力P1为0.1-10kPa。
CN201910364713.9A 2019-08-02 2019-08-02 一种利用二维切片计算岩心三维孔隙度和渗透率的方法 Active CN110222368B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910364713.9A CN110222368B (zh) 2019-08-02 2019-08-02 一种利用二维切片计算岩心三维孔隙度和渗透率的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910364713.9A CN110222368B (zh) 2019-08-02 2019-08-02 一种利用二维切片计算岩心三维孔隙度和渗透率的方法

Publications (2)

Publication Number Publication Date
CN110222368A true CN110222368A (zh) 2019-09-10
CN110222368B CN110222368B (zh) 2021-09-17

Family

ID=67820560

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910364713.9A Active CN110222368B (zh) 2019-08-02 2019-08-02 一种利用二维切片计算岩心三维孔隙度和渗透率的方法

Country Status (1)

Country Link
CN (1) CN110222368B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110599477A (zh) * 2019-09-16 2019-12-20 南昌工程学院 一种植草混凝土有效孔隙率的检测方法及系统
CN113418848A (zh) * 2021-05-31 2021-09-21 中国石油天然气股份有限公司 岩心的渗透参数的确定方法、装置、设备及存储介质
CN117494597A (zh) * 2023-10-31 2024-02-02 中国石油大学(北京) 一种三维数字岩心渗透率的确定方法、装置及设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100135536A1 (en) * 2008-11-29 2010-06-03 Jack Dvorkin Method for determining permeability of rock formation using computer tomograpic images thereof
CN104596911A (zh) * 2015-02-05 2015-05-06 山东大学 一种基于图像处理的透水混凝土孔隙特征的检测方法
CN104885124A (zh) * 2012-12-25 2015-09-02 因特维芯研究与发展有限公司 产生用于分析渗透性特征的多孔材料样本的三维特征模型的方法
KR20160069379A (ko) * 2014-12-08 2016-06-16 공주대학교 산학협력단 암석의 박편 영상을 이용한 암석의 투수율 및 공극률 산출 장치 및 그 방법
CN108133474A (zh) * 2017-12-22 2018-06-08 西安石油大学 基于岩心样品二维孔隙图像的渗透率预测方法
CN109916794A (zh) * 2017-12-13 2019-06-21 重庆润泽医药有限公司 基于工业ct的多孔材料孔隙率的无损检测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100135536A1 (en) * 2008-11-29 2010-06-03 Jack Dvorkin Method for determining permeability of rock formation using computer tomograpic images thereof
CN104885124A (zh) * 2012-12-25 2015-09-02 因特维芯研究与发展有限公司 产生用于分析渗透性特征的多孔材料样本的三维特征模型的方法
KR20160069379A (ko) * 2014-12-08 2016-06-16 공주대학교 산학협력단 암석의 박편 영상을 이용한 암석의 투수율 및 공극률 산출 장치 및 그 방법
CN104596911A (zh) * 2015-02-05 2015-05-06 山东大学 一种基于图像处理的透水混凝土孔隙特征的检测方法
CN109916794A (zh) * 2017-12-13 2019-06-21 重庆润泽医药有限公司 基于工业ct的多孔材料孔隙率的无损检测方法
CN108133474A (zh) * 2017-12-22 2018-06-08 西安石油大学 基于岩心样品二维孔隙图像的渗透率预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘馨泽: "基于VTK技术的三维岩心可视化研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110599477A (zh) * 2019-09-16 2019-12-20 南昌工程学院 一种植草混凝土有效孔隙率的检测方法及系统
CN113418848A (zh) * 2021-05-31 2021-09-21 中国石油天然气股份有限公司 岩心的渗透参数的确定方法、装置、设备及存储介质
CN117494597A (zh) * 2023-10-31 2024-02-02 中国石油大学(北京) 一种三维数字岩心渗透率的确定方法、装置及设备

Also Published As

Publication number Publication date
CN110222368B (zh) 2021-09-17

Similar Documents

Publication Publication Date Title
Pan et al. Pore-scale modeling of saturated permeabilities in random sphere packings
Al-Kharusi et al. Network extraction from sandstone and carbonate pore space images
RU2639727C2 (ru) Системы и способы цифрового анализа горных пород с определением сэо многофазного потока
Mascini et al. Event-based contact angle measurements inside porous media using time-resolved micro-computed tomography
CN110222368A (zh) 一种利用二维切片计算岩心三维孔隙度和渗透率的方法
Gerke et al. Improving watershed-based pore-network extraction method using maximum inscribed ball pore-body positioning
CN104619952A (zh) 具可靠多相渗透性决定数字岩石分析系统及方法
CN106837315B (zh) 裂缝性碳酸盐岩基质与裂缝耦合作用表征方法
Jiang et al. Pore geometry characterization by persistent homology theory
CN106886682B (zh) 用于单裂隙中溶质运移数值模拟的随机行走粒子追踪方法
CN104533370A (zh) 压裂水平井油藏、裂缝、井筒全耦合模拟方法
Van Marcke et al. An improved pore network model for the computation of the saturated permeability of porous rock
CN109900617B (zh) 一种基于声电成像测井图的裂缝性地层渗透率曲线计算方法
CN109241646B (zh) 基于椭圆堆叠和随机场的多因素二维土石混合体生成方法
CN112132795B (zh) 一种基于LiDAR点云的电塔受灾风险评估方法及系统
CN105095857A (zh) 基于关键点扰动技术的人脸数据增强方法
CN103971395B (zh) 一种纤维过滤介质结构的拟态化重建及其性能计算方法
Karimpouli et al. A hierarchical sampling for capturing permeability trend in rock physics
CN112861890B (zh) 一种储层评价模型构建方法及储层识别方法
Jin et al. Statistics of highly heterogeneous flow fields confined to three-dimensional random porous media
Hearn Perspectives in coral reef hydrodynamics
CN109632604B (zh) 一种孔隙尺度到岩心尺度聚合物驱相对渗透率粗化方法
Suzuki et al. Characterization of 3D printed fracture networks
CN110609046B (zh) 一种基于微ct水驱油图像的毛管力预测方法
CN112903555A (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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Wang Sen

Inventor after: Feng Qihong

Inventor after: Zhao Yunchang

Inventor after: Chaluming

Inventor before: Feng Qihong

Inventor before: Zhao Yunchang

Inventor before: Wang Sen

Inventor before: Chaluming

GR01 Patent grant
GR01 Patent grant