CN106202695B - 一种采用数字岩心模拟计算岩心渗透率的方法 - Google Patents

一种采用数字岩心模拟计算岩心渗透率的方法 Download PDF

Info

Publication number
CN106202695B
CN106202695B CN201610531428.8A CN201610531428A CN106202695B CN 106202695 B CN106202695 B CN 106202695B CN 201610531428 A CN201610531428 A CN 201610531428A CN 106202695 B CN106202695 B CN 106202695B
Authority
CN
China
Prior art keywords
core
permeability
grid
simulated domain
simulation
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
Application number
CN201610531428.8A
Other languages
English (en)
Other versions
CN106202695A (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.)
Cleanergy Aike (Shenzhen) Energy Technology Co. Ltd
Original Assignee
Cleanergy Aike (shenzhen) Energy Technology Co Ltd
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 Cleanergy Aike (shenzhen) Energy Technology Co Ltd filed Critical Cleanergy Aike (shenzhen) Energy Technology Co Ltd
Priority to CN201610531428.8A priority Critical patent/CN106202695B/zh
Publication of CN106202695A publication Critical patent/CN106202695A/zh
Application granted granted Critical
Publication of CN106202695B publication Critical patent/CN106202695B/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/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design

Landscapes

  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种采用数字岩心模拟计算岩心渗透率的方法，其包括如下步骤：a.采用波尔兹曼网格法模拟岩心多孔介质的物理模型；b.对该岩心多孔介质的物理模型进行物理场初始化及边界处理，该边界处理采用镜像边界作为波尔兹曼网格法的边界条件，以获得模拟区域及与该模拟区域具有相同构型的若干镜像区域；c.在该模拟区域内使用质量力驱使流体流动；d.根据处理后的该岩心多孔介质的物理模型，计算该模拟区域内的密度分布函数；及e.根据该质量力和该密度分布函数模拟计算该岩心渗透率。上述方法不仅提高了岩心渗透率计算的准确度，还缩短了该岩心渗透率的计算时间。

Description

一种采用数字岩心模拟计算岩心渗透率的方法
技术领域
本发明涉及一种采用数字岩心模拟计算岩心渗透率的方法。
背景技术
岩石的孔隙度、渗透率和含油饱和度在油气储量评估中至关重要，其数据准确性会直接影响到油气资源的储量评估，甚至油气公司的股价与金融借贷能力等。
目前，岩心特征参数的传统测试方法是通过实验室直接测量获得，虽然测试的结果准确度高，但是存在测试费用高、耗时长等问题。尤其是对于低渗透油藏，准确测量岩心渗透率可能需要几天甚至数月。
数字岩心模拟是当前石油天然气勘探开发领域的最新前沿技术之一。目前采用数字岩心模拟计算岩石渗透率主要有三种方法：孔隙网格法、有限元法和波尔兹曼网格法(Lattice Boltzmann Method,LBM)。
采用孔隙网络法进行岩石渗透率模拟(Evaluation Of The Reliability OfPrediction Of Petrophysical Data Through Imagery And Pore Network Modelling,PETROPHYSICS,50(4),322-334,2009)，其渗透率测量的准确性很大程上受制于岩石孔隙结构的简化。
有限元法对流体在岩石三维多孔介质几何模型进行数值模拟，即对流体在孔道中微观流动的不可压缩方程N-S方程进行求解，得到流体流动的速度场和压力分布场，从而可计算得到岩心渗透率(一种低渗透储层原油边界层的模拟方法，申请号：201410665585.9)，虽然其模拟计算速度快，但其模拟方法复杂。
采用波尔兹曼网格法进行岩石渗透率计算，通常设置流入和流出的密度分布函数差来产生压降以驱使流体流动，并使用循环边界条件(Pressure boundary condition ofthe lattice Boltzmann method for fully developed periodic flows.J Zhang，DYKwok,Physical Review E,2006,73(2))，然模拟体积小于岩石实际被扫描体积，计算结果对于被测岩石所在油气储层渗透率的代表性有限。综上，现有的岩石渗透率模拟计算方法存在计算准确性不高、计算耗时长、模拟方法复杂以及计算结果代表性有限等问题。
发明内容
鉴于以上内容，有必要提供一种准确、快速的采用数字岩心模拟计算岩心渗透率的方法。
一种采用数字岩心模拟计算岩心渗透率的方法，其包括如下步骤：
a.采用波尔兹曼网格法模拟岩心多孔介质的物理模型；
b.对该岩心多孔介质的物理模型进行物理场初始化及边界处理，该边界处理采用镜像边界作为波尔兹曼网格法的边界条件，以获得模拟区域及与该模拟区域具有相同构型的若干镜像区域；
c.在该模拟区域内使用质量力驱使流体流动；
d.根据处理后的该岩心多孔介质的物理模型，计算该模拟区域内的密度分布函数；及
e.根据该质量力和该密度分布函数模拟计算该岩心渗透率。
在模拟岩心多孔介质的物理模型之前还包括扫描岩心样品，并获取岩心多孔介质的横截面的二维图像，以及将该岩心多孔介质的横截面的二维图像转换为波尔兹曼网格参数，并根据该波尔兹曼网格参数建立该岩心多孔介质的物理模型。
该岩心渗透率的计算公式为：其中，k表示为渗透率，F表示为质量力，fi表示为密度分布函数，ei表示为速度函数，τ为弛豫时间，Δx表示为扫描分辨率。
该物理场初始化包括对该密度函数、该质量力和该弛豫时间进行初始化。
该质量力的范围为0.001-0.01N。
该模拟区域内设有n个网格，每一网格具有不同的密度分布函数，将每一网格的密度分布函数进行计算并叠加，以得到该模拟区域内的流体流速。
该模拟区域内的第i网格的密度分布函数是通过相邻网格的密度分布函数计算得到，并表示为fi＝f(fi-1,fi+1)，其中，n≧i≧1。
在流动方向上，使用该模拟区域内的第n-1网格(或第2网格)的密度分布函数计算该镜像区域内的未知的第n+1网格(或第0网格)的密度分布函数，垂直该流体流动的方向上，使用已知的该模拟区域内的第n-1网格(或第2网格)的密度分布函数代替该镜像区域内的未知的第n+1网格(或第0网格)的密度分布函数，其中，该模拟区域内的第n-1(或第2网格)网格与该镜像区域内的第n+1(或第0网格)网格呈镜像分布。
该采用数字岩心模拟计算岩心渗透率的方法，还包括将该模拟区域内的每一网格的密度分布函数进行流速计算并叠加，以得到该模拟区域内的流速，并用于该岩心渗透率的计算。
与该模拟区域边界相邻的镜像区域和该模拟区域呈镜像分布。
相较现有技术，上述采用数字岩心模拟计算岩心渗透率的方法，通过使用镜像边界作为波尔兹曼网格法的边界条件，以获得与模拟区域具有相同构型的镜像区域，并根据该模拟区域内的密度分布函数和质量力模拟计算，以求得岩心渗透率。本发明采用镜像边界为边界条件，其能够保持流体流动始终畅通。此外，采用质量力驱动流体流动，能够使流体流动快速达到平衡。因此，本发明的采用数字岩心模拟计算岩心渗透率的方法不仅能够提高岩心渗透率计算的准确率，还能够缩短模拟计算的时间。
附图说明
图1是本发明的一较佳实施例的采用模拟岩心模拟计算岩心渗透率的流程图。
图2是现有的波尔兹曼网格法采用循环边界条件的示意图。
图3是本发明的一较佳实施例的波尔兹曼网格法采用镜像边界条件的示意图。
图4是图2中的网格采用波尔兹曼网格法设置压降驱使流动的示意图。
图5是图3中的网格采用波尔兹曼网格法设置质量力驱使流动的示意图。
主要元件符号说明
步骤 S10，S20，S30，S40，S50，S60，S70
如下具体实施方式将结合上述附图进一步说明本发明。
具体实施方式
请参阅图1，本发明的一较佳实施例的采用模拟岩心模拟计算岩心渗透率的方法，其包括以下步骤：
S10、扫描岩心样品，并获得岩心多孔介质的横截面的二维图像。
在本实施例中，该扫描采用微计算机断层扫描技术(micro computedtomography，Micro-CT)进行扫描。
可以理解的，该Micro-CT扫描的分辨率的范围为0.5微米至3微米，该扫描后的岩心样品选取有效物理体积的范围为1毫米x 1毫米x 1毫米至1毫米x 1毫米x 3毫米。
在本实施例中，该Micro-CT扫描的分辨率为3微米，该扫描后的岩心样品选取有效物理体积为1毫米x 1毫米x 3毫米。
本领域技术人员能够理解的，该岩心样品的Micro-CT扫描图像是由该岩心样品图像和黑色背景组成的灰度图像，由于该岩心样品图像的灰度值与背景灰度值有明显差异。因此，对Micro-CT扫描图像进行二值化处理，可以得到相应的二值图像，并且可以进一步读取该岩心样品的横截面的二维图像。
S20、三维重构该岩心样品的横截面的二维图像，并读取表征该岩心多孔介质的参数。
可以理解的，三维重构该岩心样品的横截面的二维图像是为了获得该岩心多孔介质的三维结构的数字文件。
进一步的，图像识别并分辨该岩心样品的三维结构的数字文件，以读取表征该岩心多孔介质的参数。
表征该岩心多孔介质的参数包括基质和孔隙度。
本领域技术人员能够理解的，该基质例如是，但不局限于中孔中渗岩心。该基质还可选用低孔低渗岩心、中孔低渗岩心、高孔高渗岩心等。
在本实施例中，该基质为中孔中渗岩心，该孔隙度为11.8％。
可以理解的，三维重构该岩心多孔介质的横截面的二维图像是为了将所述该岩心多孔介质的横截面的二维图像转换为波尔兹曼网格参数，并根据该波尔兹曼网格参数建立该岩心多孔介质的物理模型。该表征岩心多孔介质的参数也即表示该波尔兹曼网格参数。
S30、将表征该岩心多孔介质的参数设置为初始条件，并根据该初始条件采用波尔兹曼网格法模拟该岩心多孔介质的物理模型。
S40、对该岩心多孔介质的物理模型进行物理场初始化及边界处理，该边界处理采用镜像边界作为波尔兹曼网格法的边界条件，以获取模拟区域及与该模拟区域具有相同构型的若干镜像区域。
该物理场初始化包括对该密度函数、该质量力和该弛豫时间进行初始化。
在本实施例中，该模拟区域是真实模拟的，该镜像区域是虚拟的。
可以理解的，相邻的镜像区域呈镜像分布，也即相邻两镜像区域呈对称分布。
如图3所示，在本实施例中，该若干镜像区域包括与该模拟区域呈镜像分布的若干虚拟的第一镜像区域及与该第一镜像区域呈镜像分布的虚拟的第二镜像区域。
可以理解的，与该模拟区域边界相邻的镜像区域与该模拟区域呈镜像分布。
该模拟区域内设有n个网格，每一网格具有不同的密度分布函数。
可以理解的，该模拟区域内的第i网格的密度分布函数是通过相邻网格的密度分布函数计算得到，可简单表示为fi＝f(fi-1,fi+1)，其中，n≧i≧1。
可以理解的，该相邻网格包括直接与该网格边界相邻的网格及与该网格斜对角对接的网格。该网格是三维立体的(图中未示)。
将该模拟区域内的每一网格密度分布函数进行计算并叠加，以得到该模拟区域内的流体流速，并用于计算该岩心的渗透率。
可以理解的，将该模拟区域内的密度分布函数进行流速计算，以得到该模拟区域内的流速，并用于计算该岩心的渗透率。
每一镜像区域内设有与该模拟区域相对应的网格，每一网格具有不同的密度分布函数。
如图2所示，现有的波尔兹曼网格法的边界条件采用循环边界条件，其流体流动不连续。
进一步的，如图3所示，本发明采用镜像边界条件作为波尔兹曼网格法的边界条件。因此，可以始终保持流体流动的畅通性，从而提高了渗透率计算的准确度。
可以理解的，该流体密度分布函数不限于本发明图3所示的呈斜对角分布，该方向可以是水平分布、竖直分布或与水平边界存在倾斜角度分布(例如，该倾斜角度为45度，也即斜对角分布)。
S50、在该模拟区域内使用质量力驱使流体流动。
如图4所示，传统的波尔兹曼网格法采用流入和流出压降来驱使流体流动，其在流体流动过程中会形成一定的压力差分布。然而，通过流入和流出压降来驱使流体流动，其流体流动达到平衡需要较长的时间，且使用循环边界条件及压力差模拟计算岩心渗透率的准确性影响较大。
可以理解的，该流体流动达到平衡是指该模拟区域内流速不随时间变化。
可以理解的，采用该质量力而非流入和流出压降来驱使流体流动，其能够使流体流动快速达到平衡。如图5所示，本发明的各网格均设置质量力来驱使流体流动，因此流体流动能够快速达到平衡，从而缩短了计算渗透率所需的时间。
该质量力的范围为0.001-0.01N。
S60、根据处理后的该岩心多孔介质的物理模型，计算该模拟区域的密度分布函数。
该计算该模拟区域内的密度分布函数的步骤包括，将该模拟区域内的每一网格密度分布函数进行计算并叠加，以得到该模拟区域内的流体流速，并用于计算该岩心的渗透率。
进一步的，该计算该模拟区域内的流速的步骤还包括，将该模拟区域内的密度分布函数进行流速计算∑fiei/∑fi并叠加，以得到该模拟区域内的流速，并用于计算该岩心的渗透率。
可以理解的，该模拟区域内的第i网格的密度分布函数是通过相邻网格的密度分布函数计算得到，可简单表示为fi＝f(fi-1,fi+1)，其中，n≧i≧1。
在本实施例中，在流动方向上，使用该模拟区域内的第n-1网格(或第2网格)的密度分布函数计算该镜像区域内的未知的第n+1网格(或第0网格)的密度分布函数，垂直该流体流动的方向上，使用已知的该模拟区域内的第n-1网格(或第2网格)的密度分布函数代替该镜像区域内的未知的第n+1网格(或第0网格)的密度分布函数，其中，该模拟区域内的第n-1(或第2网格)网格与该镜像区域内的第n+1(或第0网格)网格呈镜像分布。
可以理解的，该镜像区域内的第0网格与该模拟区域内的第1网格相邻，该镜像区域内的第n+1网格与该模拟区域内的第n网格相邻。
S70、根据该质量力和该密度分布函数模拟该岩心渗透率的计算公式，根据该计算公式求得该岩心渗透率。
该岩心渗透率的计算公式为：其中，k表示为渗透率，F表示为质量力，fi表示为密度分布函数，ei表示为速度函数，τ为弛豫时间，Δx表示为扫描分辨率。
可以理解的，该速度函数ei表示为LBM的速度函数。
本发明采用数字岩心模拟计算岩心渗透率的方法，通过使用镜像边界作为波尔兹曼网格法的边界条件，以获得若干模型对称的虚拟模拟区域，并根据该模拟区域内的密度分布函数和质量力模拟计算，以求得岩心渗透率。本发明采用镜像边界为边界条件，其能够保持流体流动始终畅通。此外，采用质量力驱动流体流动，能够使流体流动快速达到平衡。因此，本发明的采用数字岩心模拟计算岩心渗透率的方法不仅能够提高岩心渗透率计算的准确率，还能够缩短模拟计算的时间。此外，本发明可以在岩心无损的情况下，基于Micro-CT扫描可以准确分析出岩心的绝对渗透率，并动态观察流体的流动过程。因此，可以广泛运用用到煤层气、致密油等非常规油气的勘探开发。
下面通过具体的实施例及上述附图对本发明作进一步的详细描述，以下实施例可以使本专业技术人员更全面的理解本发明，但不以任何方式限制本发明。
在对比例1-2及本发明的实施例中，采用Micro CT扫描该岩心样品，其中，该Micro-CT扫描的分辨率为0.5微米，该岩心样品的有效物理体积1毫米立方，孔隙度为11.8％。
对比例1
采用传统的实验室气渗透率测量法测量，当注入气达到稳态后，测量流入/流出压力和流量，并使用达西方程计算渗透率，其测得岩心渗透率为121mD。
对比例2
使用传统的波尔兹曼网格方法，使用循环边界作为LBM的边界条件，并使用压力差为0.4kg/m3驱动流体流动，并使用岩心渗透率计算公式：其中，k表示为渗透率，ΔP表示为压力差，fi表示为密度分布函数，ei表示为速度函数，τ为弛豫时间，Δx表示为扫描分辨率。使用传统的波尔兹曼网格方法，其运算岩心渗透率的时间约为7个小时，计算岩心渗透率为89mD。
实施例
采用本发明采用数字岩心模拟计算岩心渗透率的方法，使用镜像边界边界作为LBM的边界条件，并使用质量力为0.01N驱动流体流动，并使用该岩心渗透率的计算公式：其中，k表示为渗透率，F表示为质量力，fi表示为密度分布函数，ei表示为速度函数，τ为弛豫时间，Δx表示为扫描分辨率。本发明采用数字岩心模拟计算岩心渗透率的方法，其运算岩心渗透率的时间时间约为3个小时，计算岩心渗透率为132mD。
性能测试
将对比文件1-2及实施例模拟计算该岩心样品的渗透率、运算该岩心样品的渗透率的时间进行测试，并比采用传统LBM模拟计算的渗透率及与的与实验室实测的渗透率进行比较，具体的测试结果如表一所示。
表一
性能参数 对比例1 对比例2 实施例
运算时间(h) —— 7 3
渗透率(mD) 121 89 132
误差(％) —— 26.4 9.1
如表一所示，对比例2计算的岩心渗透率为89mD，远小于对比例1测得的岩心渗透率，其相对误差为26.4％。然本发明实施例1计算的岩心渗透率为132mD，略高于对比例1测得的岩心渗透率，其相对误差为9.1％。因此，测试结果说明利用本发明的采用数字岩心模拟计算岩心渗透率的方法计算岩心渗透率，其计算结果更接近实际测量的岩心渗透率。相较于对比例2，本发明的特征岩心渗透率模拟计算结果的误差可减小至10％以内，且计算速度提高1倍以上。
本发明采用数字岩心模拟计算岩心渗透率的方法，通过使用镜像边界作为波尔兹曼网格法的边界条件，以获得若干与模拟区域具有相同构型的虚拟的模拟区域，并根据真实的模拟区域内的密度分布函数和质量力模拟计算，以求得岩心渗透率。不仅能够提高岩心渗透率测量的准确度，还缩短了模拟计算的时间。此外，与现有的数字岩心模拟计算方法相比，在不改变被扫描岩心物理体积的前提下，本发明采用的数字岩心模拟可扩大表征体元(Representative Elementary Volume，REV)体积。
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，以上实施方式仅是用于解释权利要求书。然本发明的保护范围并不局限于说明书。任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可轻易想到的变化或者替换，都包含在本发明的保护范围之内。

Claims (10)

1.一种采用数字岩心模拟计算岩心渗透率的方法，其包括如下步骤：
a.采用波尔兹曼网格法模拟岩心多孔介质的物理模型；
b.对该岩心多孔介质的物理模型进行物理场初始化及边界处理，该边界处理采用镜像边界作为波尔兹曼网格法的边界条件，以获取模拟区域及与该模拟区域具有相同构型的若干镜像区域；
c.在该模拟区域内各网格均设置质量力驱使流体流动，以使该流体在该模拟区域内快速达到平衡；
d.根据处理后的该岩心多孔介质的物理模型，计算该模拟区域内的密度分布函数；及
e.根据该质量力和该密度分布函数模拟计算该岩心渗透率。
2.如权利要求1所述的采用数字岩心模拟计算岩心渗透率的方法，其特征在于：在模拟岩心多孔介质的物理模型之前还包括扫描岩心样品，并获取岩心多孔介质的横截面的二维图像，以及将该岩心多孔介质的横截面的二维图像转换为波尔兹曼网格参数，并根据该波尔兹曼网格参数建立该岩心多孔介质的物理模型。
3.如权利要求1所述的采用数字岩心模拟计算岩心渗透率的方法，其特征在于：该岩心渗透率的计算公式为：其中，k表示为渗透率，F表示为质量力，fi表示为密度分布函数，ei表示为速度函数，τ为弛豫时间，Δx表示为扫描分辨率。
4.如权利要求3所述的采用数字岩心模拟计算岩心渗透率的方法，其特征在于：该物理场初始化包括对该密度函数、该质量力和该弛豫时间进行初始化。
5.如权利要求1所述的采用数字岩心模拟计算岩心渗透率的方法，其特征在于：该质量力的范围为0.001-0.01N。
6.如权利要求1所述的采用数字岩心模拟计算岩心渗透率的方法，其特征在于：该模拟区域内设有n个网格，每一网格具有不同的密度分布函数，将每一网格的密度分布函数进行计算并叠加，以得到该模拟区域内的密度分布函数。
7.如权利要求6所述的采用数字岩心模拟计算岩心渗透率的方法，其特征在于：该模拟区域内的第i网格的密度分布函数是通过相邻网格的密度分布函数计算得到，并表示为fi＝f(fi-1,fi+1)，其中，n大于等于i，i大于等于1。
8.如权利要求7所述的采用数字岩心模拟计算岩心渗透率的方法，其特征在于：在流动方向上，使用该模拟区域内的第n-1网格的密度分布函数计算该镜像区域内的未知的第n+1网格的密度分布函数，垂直该流体流动的方向上，使用已知的该模拟区域内的第n-1网格的密度分布函数代替该镜像区域内的未知的第n+1网格的密度分布函数，其中，该模拟区域内的第n-1网格与该镜像区域内的第n+1网格呈镜像分布。
9.如权利要求6所述的采用数字岩心模拟计算岩心渗透率的方法，其特征在于：还包括将该模拟区域内的每一网格的密度分布函数进行流速计算并叠加，以得到该模拟区域内的流速，并用于该岩心渗透率的计算。
10.如权利要求1所述的采用数字岩心模拟计算岩心渗透率的方法，其特征在于：与该模拟区域边界相邻的镜像区域和该模拟区域呈镜像分布。
CN201610531428.8A 2016-07-07 2016-07-07 一种采用数字岩心模拟计算岩心渗透率的方法 Active CN106202695B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610531428.8A CN106202695B (zh) 2016-07-07 2016-07-07 一种采用数字岩心模拟计算岩心渗透率的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610531428.8A CN106202695B (zh) 2016-07-07 2016-07-07 一种采用数字岩心模拟计算岩心渗透率的方法

Publications (2)

Publication Number Publication Date
CN106202695A CN106202695A (zh) 2016-12-07
CN106202695B true CN106202695B (zh) 2018-04-20

Family

ID=57472455

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610531428.8A Active CN106202695B (zh) 2016-07-07 2016-07-07 一种采用数字岩心模拟计算岩心渗透率的方法

Country Status (1)

Country Link
CN (1) CN106202695B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022053991A1 (en) * 2020-09-10 2022-03-17 Khalifa University of Science and Technology Morphology decoder to predict heterogeneous rock permeability with machine learning guided 3d vision

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106960070B (zh) * 2016-12-28 2020-02-21 山东科技大学 一种基于有限元-离散元ct重构煤体的渗流模拟方法
CN108267390B (zh) * 2016-12-30 2020-07-10 中国石油天然气股份有限公司 一种含纳米孔隙储层的气体渗透率确定方法
CN109426687A (zh) * 2017-08-29 2019-03-05 中国石油化工股份有限公司 一种油藏模拟网格转换方法及系统
MX2020004533A (es) * 2017-11-06 2020-08-20 Univ Khalifa Science & Technology Metodo y sistema para determinar la permeabilidad de un medio poroso.
CN108597018A (zh) * 2018-04-28 2018-09-28 清能艾科（深圳）能源技术有限公司 对土壤实现数值模拟的方法和装置
CN108876898B (zh) * 2018-04-28 2022-06-03 清能艾科（深圳）能源技术有限公司 实现原油饱和度预测的方法和装置、机器设备
CN109211666B (zh) * 2018-08-31 2019-12-03 山东科技大学 基于ct扫描的预测应力加载条件下煤体渗透率的方法
CN110135311A (zh) * 2019-05-06 2019-08-16 重庆科技学院 一种基于三维岩心扫描图像的孔隙与孔喉识别系统及方法
CN110441209B (zh) * 2019-08-13 2020-12-11 中国石油大学(北京) 一种基于致密储层数字岩心计算岩石渗透率的方法
CN111624147B (zh) * 2020-04-16 2023-04-07 中国石油天然气股份有限公司 岩心的相对渗透率测定方法及装置
CN112082917B (zh) * 2020-08-03 2021-04-30 西南石油大学 一种基于动态网络模拟的气水非稳态两相渗流模拟方法
CN117494597A (zh) * 2023-10-31 2024-02-02 中国石油大学(北京) 一种三维数字岩心渗透率的确定方法、装置及设备

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8081802B2 (en) * 2008-11-29 2011-12-20 Ingrain, Inc. Method for determining permeability of rock formation using computer tomograpic images thereof
CN101717625A (zh) * 2009-11-23 2010-06-02 大庆石油学院 一种适用于低渗透率油藏的粘土稳定剂

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022053991A1 (en) * 2020-09-10 2022-03-17 Khalifa University of Science and Technology Morphology decoder to predict heterogeneous rock permeability with machine learning guided 3d vision

Also Published As

Publication number Publication date
CN106202695A (zh) 2016-12-07

Similar Documents

Publication Publication Date Title
CN106202695B (zh) 一种采用数字岩心模拟计算岩心渗透率的方法
US10830713B2 (en) System and methods for computing physical properties of materials using imaging data
Zhao et al. Pore structure characterization of coal by synchrotron radiation nano-CT
CN104619952B (zh) 具可靠多相渗透性决定数字岩石分析系统及方法
AU2011258594B2 (en) Method for obtaining consistent and integrated physical properties of porous media
CA2911247C (en) Digital core sensitivity analysis
CN106837315B (zh) 裂缝性碳酸盐岩基质与裂缝耦合作用表征方法
CN106716114A (zh) 多孔介质分析系统和方法
CN105510204A (zh) 一种基于ct图像的渗透率预测方法
Shikhov et al. Evaluation of capillary pressure methods via digital rock simulations
WO2014104909A1 (ru) Способ получения характеристической трехмерной модели образца пористого материала для исследования свойств проницаемости
CN106093083A (zh) 一种采用数字岩心模拟建立岩电关系的方法
CN104737037A (zh) 使用多点模拟的储层建模
US20160306074A1 (en) Tuning digital core analysis to laboratory results
Hajizadeh et al. An algorithm for 3D pore space reconstruction from a 2D image using sequential simulation and gradual deformation with the probability perturbation sampler
Keogh et al. Data capture for multiscale modelling of the Lourinha Formation, Lusitanian Basin, Portugal: an outcrop analogue for the Statfjord Group, Norwegian North Sea
Ji et al. A multiscale reconstructing method for shale based on SEM image and experiment data
CN113167713B (zh) 数字表征岩石渗透率的方法
Jouini et al. Simulation of elastic properties in carbonates
Kang et al. Construction of complex digital rock physics based on full convolution network
US11119025B2 (en) Method for characterizing the porosity of rock
CN111894568A (zh) 裂缝性碳酸盐岩油藏饱和度模型数字岩心分析方法
Naraghi et al. 3D reconstruction of porous media from a 2D section and comparisons of transport and elastic properties
Yao et al. Upscaling of carbonate rocks from micropore scale to core scale
Zhao et al. Fractal characterization and petrophysical analysis of 3D dynamic digital rocks of sandstone

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20170516

Address after: 518057 Shenzhen Institute of technology, Nanshan District high tech Industrial Park, Guangdong,, Tsinghua University, A302

Applicant after: Shenzhen Institute of Stinghua University

Address before: 5 B block 31F, 518040 deep Kang village, 125 overseas Chinese village, Chengdong, Guangdong, Shenzhen, Futian District

Applicant before: Long Wei

Applicant before: Zan Cheng

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20170821

Address after: Nanshan District Guangdong streets of Shenzhen science and technology of Guangdong Province in 518054 southern Shenzhen Research Institute of Tsinghua University, A304-1

Applicant after: Cleanergy Aike (Shenzhen) Energy Technology Co. Ltd

Address before: 518057 Shenzhen Institute of technology, Nanshan District high tech Industrial Park, Guangdong,, Tsinghua University, A302

Applicant before: Shenzhen Institute of Stinghua University

GR01 Patent grant
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20161207

Assignee: Shenzhen Institute of Stinghua University

Assignor: Cleanergy Aike (Shenzhen) Energy Technology Co. Ltd

Contract record no.: 2018440020080

Denomination of invention: Method for calculating core permeability by adopting digital core simulation

Granted publication date: 20180420

License type: Exclusive License

Record date: 20181015