CN110441209A - 一种基于致密储层数字岩心计算岩石渗透率的方法 - Google Patents
一种基于致密储层数字岩心计算岩石渗透率的方法 Download PDFInfo
- Publication number
- CN110441209A CN110441209A CN201910742867.7A CN201910742867A CN110441209A CN 110441209 A CN110441209 A CN 110441209A CN 201910742867 A CN201910742867 A CN 201910742867A CN 110441209 A CN110441209 A CN 110441209A
- Authority
- CN
- China
- Prior art keywords
- rock
- permeability
- micro
- compact reservoir
- 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.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000035699 permeability Effects 0.000 title claims abstract description 40
- 239000012530 fluid Substances 0.000 claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims abstract description 13
- 238000004088 simulation Methods 0.000 claims abstract description 9
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 7
- 238000004141 dimensional analysis Methods 0.000 claims abstract description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 6
- 239000011148 porous material Substances 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000002591 computed tomography Methods 0.000 claims description 3
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 3
- 229910052753 mercury Inorganic materials 0.000 claims description 3
- 229910052757 nitrogen Inorganic materials 0.000 claims description 3
- 238000001179 sorption measurement Methods 0.000 claims description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000001143 conditioned effect Effects 0.000 description 1
- 238000000280 densification Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
-
- 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
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- 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/11—Region-based segmentation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Dispersion Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biochemistry (AREA)
- Fluid Mechanics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
Abstract
本发明涉及一种基于致密储层数字岩心计算岩石渗透率的方法，包括以下步骤：(1)基于微纳米CT实验结果，建立致密储层岩心的复杂数字岩心，其中包括：岩石骨架、孔隙空间和微纳孔系统；(2)基于所建立的数字岩心，利用灰色格子Boltzmann方法模拟流体的流动；(3)流动方向的边界条件为循环边界，并且外加作用力，其他方向为封闭边界条件；(4)根据计算的流速和压力结果，结合量纲分析，得到岩石的渗透率。本发明利用灰色格子Boltzmann方法，在传统算法的基础上，综合考虑了微纳米孔隙的流动，高效准确地实现了岩石渗透率的求取，方法新颖、可推广性高。
Description
技术领域
本发明涉及岩土领域，具体的说涉及一种基于致密储层数字岩心计算岩石渗透率的方法。
背景技术
致密油藏的开发已成为石油工业增储上产的关键领域。致密储层渗透率的大小直接关系到油藏开发方案的选取，是储层重要的物性参数之一。目前室内试验测量难度大，稳态法几乎难以实施，脉冲测量方法不确定性高，而且会掺杂气体的非达西效应。基于致密储层的数字岩心来模拟计算岩石的渗透率被认为是一种有效的替代方法。但目前通常采用的微纳米CT精度不足，只能识别比较大的粒间孔隙，而对于微纳米孔隙系统无法识别，但致密岩心很多时候微纳米孔占比高达一半以上，而且微纳米孔起到了沟通粒间孔隙的作用，在此数字岩心上采用传统的格子Boltzmann方法会忽略微纳米孔对流动储存和流动的贡献，计算结果势必会致密储层复杂数字岩心的建立有很大的误差产生。
发明内容
本发明说是要解决的技术问题至少包括：现有技术中的岩石渗透率的计算会产生很大的误差的问题。
本发明提出一种基于致密储层数字岩心计算岩石渗透率的方法，包括以下步骤：
(1)致密储层复杂数字岩心的建立：基于微纳米CT实验结果，建立致密储层岩心的复杂数字岩心，其中包括：岩石骨架、孔隙空间和微纳孔系统；
(2)利用灰色格子Boltzmann方法模拟流体的流动：基于所建立的数字岩心，利用灰色格子Boltzmann方法模拟流体的流动；
(3)流动条件的确定及流速的计算：流动方向的边界条件为循环边界，并且外加作用力，其他方向为封闭边界条件；
(4)岩石的渗透率的计算：根据计算的流速和压力结果，结合量纲分析，得到岩石的渗透率。
所述步骤1)中，利用微纳米CT扫描岩样得到岩心的三维数据体，依据体素灰度值的大小分布，识别出岩石骨架、孔隙空间和微纳米孔隙系统(大约对应于灰度直方图上的三个峰)，分别标记为1、0和0.5，进行图像分割，这样就得到了致密储层的复杂数字岩心。
所述步骤2)中，在已生成的数字岩心上，基于灰色格子Boltzmann方法进行流动模拟计算。体素点为0时，流体存在，体素点为1时采用经典的half-way反弹格式，而体素点为0.5时(灰色格子)采用部分反弹格式，即乘以部分反弹率ns，当ns＝0，即为孔隙空间，不存在反弹；当ns＝1，即为岩石骨架，完全反弹；部分反弹介于两者之间。
部分反弹格式中的反弹率ns需要依据微纳米孔隙的渗透率获得。首先可通过氮气吸附或高压压汞所得到的孔隙半径进行估算得到微纳米孔隙的渗透率，转化到格子单位可以得到对应的kmicro，继而ns＝v/(2kmicro+v)。
所述步骤3)中，流动方向上采用循环边界条件，而其他方向均采用封闭边界条件，即添加一层封闭的岩石骨架。在流动方向上施加作用力，在LBGK模型中，对应的作用力表达式为：其中pg为流动方向的压力梯度；ei是格子的方向向量，fi为粒子在各个方向上的密度分布函数。
在模拟计算过程中，时间步长取1ts，格子长度取1lu，密度取1mu/lu3，松弛因子τ取1，对应的动力粘度为保证算法的精度，调节pg的大小，使得任一格子中流体流速不高于0.1lu/ts。以实际体素点尺寸，密度和粘度为基础，即可建立格子单位与SI单位之间的关系，继而实现了两者的转化。
所述步骤4中，经过一定的迭代计算，等到格子的流速不发生变化时，计算流动方向上所有格子流速的平均值，即则可得到渗透率按照量纲转化，就可得到实际岩心的渗透率。
综上，本发明所设计一种基于致密储层数字岩心计算岩石渗透率的方法，该方法利用灰色格子Boltzmann方法，在传统算法的基础上，综合考虑了微纳米孔隙的流动，高效准确地实现了岩石渗透率的求取，方法新颖、可推广性高。本发明克服了传统计算中在此数字岩心上采用传统的格子Boltzmann方法会忽略微纳米孔对流动储存和流动的贡献而产生很大的误差的问题，计算的准确度更高。
附图说明
图1是本发明实施方式的基于致密储层数字岩心计算岩石渗透率的方法的流程图；
图2是本发明所涉及致密储层复杂数字岩心的示意图。
具体实施方式
下面结合具体实施方式对本发明进行进一步的详细描述，给出的实施例仅为了阐明本发明，而不是为了限制本发明的范围。以下提供的实施例可作为本技术领域普通技术人员进行进一步改进的指南，并不以任何方式构成对本发明的限制。
如图1所示，本发明的一种基于致密储层数字岩心计算岩石渗透率的方法，其包括以下步骤：
S10、采集并制作致密储层岩样，利用微纳米CT扫描岩样，得到由体素点组成的三维数字岩心数据体。根据三维数字岩心数据体绘制体素灰度值的频率分布图，从频率分布图中可以发现，所述频率分布图中存在三个峰，分别对应为岩石骨架、孔隙空间和微纳米孔隙系统。由于微纳米孔隙系统一般低于所用CT的分辨率，因此其体素点的灰度值介于岩石骨架和孔隙空间之间。根据频率分布图，选取三个峰之间的分界点作为界限，对原三维灰度数据体进行分割，最终得到致密储层的三维数字岩心，对三类体素点(岩石骨架、孔隙空间和微纳米孔隙系统)进行标注，如图2所示，图2是本发明所涉及致密储层的复杂数字岩心的示意图，图2上的黑色白色和灰色区域分别对应岩石骨架、孔隙空间和微纳米孔隙系统，并将其标注为1、0和0.5。
S20、依据步骤S10生成的致密储层的复杂数字岩心，以体素点中心为格子中心，建立灰色格子Boltzmann模拟方法的基本网格。灰色格子Boltzmann模拟方法是在LBGK模型的基础上发展而来的，LBGK模型是格子Boltzmann方法应用最为广泛的一类模拟算法。所谓灰色格子Boltzmann就是在微纳米孔隙系统的体素点采用介于0/1的反射条件。当体素点为0时，流体存在，体素点为1时采用经典的half-way反弹格式，而体素点为0.5时(灰色格子)采用部分反弹格式，即乘以部分反弹率ns，当ns＝0，即为孔隙空间，不存在反弹；当ns＝1，即为岩石骨架，完全反弹；部分反弹介于两者之间。
S30、选择一个流动方向，在流动方向上采用循环边界条件，而其他方向均采用封闭边界条件，即添加一层封闭的岩石骨架。在流动方向上施加作用力用以等效代替压力梯度的施加，在LBGK模型中，对应的作用力表达式为：其中pg为流动方向的压力梯度；ei是格子的方向向量，fi为粒子在各个方向上的密度分布函数，Fi表示体积力，τ表示松弛时间，ρ表示流体密度，u表示粒子速度，eq表示平衡态。初始时，设定格子的密度分布为1-ns，即，孔隙空间格子的流体密度为1，岩石骨架格子的流体密度为0，而微纳米孔格子的流体密度为1-ns。在存在边界条件的影响下，即可进行迭代计算直至流动稳定。
其中，在模拟计算过程中，时间步长取1ts，格子长度取1lu，密度取1mu/lu3，松弛因子τ取1，此时，对应的动力粘度其中，v表示动力粘度，为保证算法的精度，调节pg的大小，使得任一格子中流体流速不高于0.1lu/ts。以实际体素点尺寸，密度和粘度为基础，即可建立格子单位与SI单位之间的关系，继而实现了两者的转化。
同时，部分反弹格式中的反弹率ns需要依据微纳米孔隙的渗透率获得。首先可通过氮气吸附或高压压汞所得到的孔隙半径进行估算得到微纳米孔隙的渗透率。此估算可以利用经典的Kozeny公式，即其中φ为孔隙度，r为孔喉半径，T是迂曲度。将得到的渗透率转化到格子单位可以得到对应的kmicro，继而通过公式ns＝v/(2kmicro+v)就可得到对应格子单位处的微纳米孔隙(灰色格子)的部分反射率。
S40、经过一定的迭代计算，等到格子的流速不发生变化时，计算流动方向上所有格子流速的平均值，即则可得到渗透率按照量纲转化(所述量纲转化是指不同的量纲之间的相互转化)，就可得到实际岩心的渗透率，此时实际岩心的渗透率为:k*(lR)2，其中lR为一个格子(一个体素点)所对应实际物理长度。
综上，本发明利用灰色格子Boltzmann方法，在传统算法的基础上，综合考虑了微纳米孔隙的流动，高效准确地实现了岩石渗透率的求取，方法新颖、可推广性高。
以上对本发明进行了详述。对于本领域技术人员来说，在不脱离本发明的宗旨和范围，以及无需进行不必要的实验情况下，可在等同参数、浓度和条件下，在较宽范围内实施本发明。虽然本发明给出了特殊的实施例，应该理解为，可以对本发明作进一步的改进。总之，按本发明的原理，本申请欲包括任何变更、用途或对本发明的改进，包括脱离了本申请中已公开范围，而用本领域已知的常规技术进行的改变。按以下附带的权利要求的范围，可以进行一些基本特征的应用。
Claims (8)
1.一种基于致密储层数字岩心计算岩石渗透率的方法，其特征在于，包括以下步骤：
1)致密储层复杂数字岩心的建立；
2)利用灰色格子Boltzmann方法模拟流体的流动；
3)流动条件的确定及流速的计算；
4)岩石的渗透率的计算。
2.如权利要求1所述的基于致密储层数字岩心计算岩石渗透率的方法，其特征在于，包括以下步骤：
1)基于微纳米CT实验结果，建立致密储层岩心的复杂数字岩心；
2)基于步骤1)所建立的数字岩心，利用灰色格子Boltzmann方法模拟流体的流动；
3)以流动方向的边界条件为循环边界，并且外加作用力，其他方向为封闭边界条件，进行迭代计算直至压力场稳定；
4)根据计算的流速和压力结果，结合量纲分析，得到岩石的渗透率。
3.如权利要求2所述的基于致密储层数字岩心计算岩石渗透率的方法，其特征在于：所述步骤1)中，利用微纳米CT扫描岩样得到由体素点组成的三维数据体，根据三维数字岩心数据体体素灰度值的大小分布，识别出岩石骨架、孔隙空间和微纳米孔隙系统，进行图像分割，得到致密储层的复杂数字岩心。
4.如权利要求2所述的基于致密储层数字岩心计算岩石渗透率的方法，其特征在于：所述步骤2)中，基于灰色格子Boltzmann方法，其中，体素点为0时，流体存在，体素点为1时采用经典的half-way反弹格式，而体素点为0.5时采用部分反弹格式，即乘以部分反弹率ns，当ns＝0，即为孔隙空间，不存在反弹；当ns＝1，即为岩石骨架，完全反弹；部分反弹介于两者之间。
5.如权利要求4所述的基于致密储层数字岩心计算岩石渗透率的方法，其特征在于：所述部分反弹率ns依据微纳米孔隙的渗透率获得：首先可通过氮气吸附或高压压汞所得到的孔隙半径进行估算得到微纳米孔隙的渗透率，转化到格子单位可以得到对应的kmicro，继而其中kmicro表示微纳米孔隙的渗透率，ν表示运动粘度。
6.如权利要求2所述的基于致密储层数字岩心计算岩石渗透率的方法，其特征在于，所述步骤3)中，在流动方向上施加作用力，在LBGK模型中，对应的作用力表达式为：其中pg为流动方向的压力梯度；ei是格子的方向向量，fi为粒子在各个方向上的密度分布函数，Fi表示体积力，τ表示松弛时间，ρ表示流体密度，u表示粒子速度，eq表示平衡态。
7.如权利要求6所述的基于致密储层数字岩心计算岩石渗透率的方法，其特征在于，所述步骤3)中，时间步长取1ts，格子长度取1lu，密度取1mu/lu3，松弛因子τ取1，对应的动力粘度为保证算法的精度，调节pg的大小，使得任一格子中流体流速不高于0.1lu/ts。以实际体素点尺寸，密度和粘度为基础，即可建立格子单位与SI单位之间的关系，继而实现了两者的转化。
8.如权利要求2所述的基于致密储层数字岩心计算岩石渗透率的方法，其特征在于：所述步骤4)中，经过所述步骤3)中的迭代计算至格子的流速不发生变化时，计算流动方向上所有格子流速的平均值，即则可得到渗透率量纲转化，得到实际岩心的渗透率，其中表示平均速度，ei表示格子的方向向量，fi表示粒子在各个方向上的密度分布函数，ρ表示流体密度，v表示运动粘度，pg表示流动方向的压力梯度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910742867.7A CN110441209B (zh) | 2019-08-13 | 2019-08-13 | 一种基于致密储层数字岩心计算岩石渗透率的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910742867.7A CN110441209B (zh) | 2019-08-13 | 2019-08-13 | 一种基于致密储层数字岩心计算岩石渗透率的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110441209A true CN110441209A (zh) | 2019-11-12 |
| CN110441209B CN110441209B (zh) | 2020-12-11 |
Family
ID=68434945
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910742867.7A Active CN110441209B (zh) | 2019-08-13 | 2019-08-13 | 一种基于致密储层数字岩心计算岩石渗透率的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110441209B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113029899A (zh) * | 2021-02-24 | 2021-06-25 | 西南石油大学 | 一种基于显微图像处理的砂岩渗透率计算方法 |
| CN113125325A (zh) * | 2021-04-26 | 2021-07-16 | 东北石油大学 | 一种煤岩裂隙特征表征与渗透性模拟方法 |
| CN114428040A (zh) * | 2021-09-13 | 2022-05-03 | 中国石油化工股份有限公司 | 一种页岩油储层储渗空间定量表征及参数获取方法 |
| CN117494597A (zh) * | 2023-10-31 | 2024-02-02 | 中国石油大学(北京) | 一种三维数字岩心渗透率的确定方法、装置及设备 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140343858A1 (en) * | 2013-05-16 | 2014-11-20 | Exa Corporation | Mass exchange model for relative permeability simulation |
| CN106202695A (zh) * | 2016-07-07 | 2016-12-07 | 龙威 | 一种采用数字岩心模拟计算岩心渗透率的方法 |
| CN108267390A (zh) * | 2016-12-30 | 2018-07-10 | 中国石油天然气股份有限公司 | 一种含纳米孔隙储层的气体渗透率确定方法 |
| CN108729908A (zh) * | 2018-05-21 | 2018-11-02 | 中国石油大学(华东) | 一种基于孔隙网络模型的致密油流动模拟及渗透率预测方法 |
| CN109242970A (zh) * | 2018-10-11 | 2019-01-18 | 中国科学院力学研究所 | 页岩lREV尺度数字岩心重构方法及装置 |
-
2019
- 2019-08-13 CN CN201910742867.7A patent/CN110441209B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140343858A1 (en) * | 2013-05-16 | 2014-11-20 | Exa Corporation | Mass exchange model for relative permeability simulation |
| CN106202695A (zh) * | 2016-07-07 | 2016-12-07 | 龙威 | 一种采用数字岩心模拟计算岩心渗透率的方法 |
| CN108267390A (zh) * | 2016-12-30 | 2018-07-10 | 中国石油天然气股份有限公司 | 一种含纳米孔隙储层的气体渗透率确定方法 |
| CN108729908A (zh) * | 2018-05-21 | 2018-11-02 | 中国石油大学(华东) | 一种基于孔隙网络模型的致密油流动模拟及渗透率预测方法 |
| CN109242970A (zh) * | 2018-10-11 | 2019-01-18 | 中国科学院力学研究所 | 页岩lREV尺度数字岩心重构方法及装置 |
Non-Patent Citations (4)
| Title |
|---|
| MINGQIANG CHEN 等: "Pore network modelling of fluid flow in tight formations considering boundary layer effect and media deformation", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 * |
| YAN-LONG ZHAO 等: "Lattice Boltzmann simulation of gas flow and permeability prediction in coal fracture networks", 《JOURNAL OF NATURAL GAS SCIENCE AND ENGINEERING》 * |
| 张潇丹 等: "REV尺度多孔介质格子Boltzmann方法的数学模型及应用的进展研究", 《化工进展》 * |
| 李如如 等: "部分渗透格子玻尔兹曼流动及渗透率模拟计算", 《山西大学学报》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113029899A (zh) * | 2021-02-24 | 2021-06-25 | 西南石油大学 | 一种基于显微图像处理的砂岩渗透率计算方法 |
| CN113029899B (zh) * | 2021-02-24 | 2022-07-01 | 西南石油大学 | 一种基于显微图像处理的砂岩渗透率计算方法 |
| CN113125325A (zh) * | 2021-04-26 | 2021-07-16 | 东北石油大学 | 一种煤岩裂隙特征表征与渗透性模拟方法 |
| CN114428040A (zh) * | 2021-09-13 | 2022-05-03 | 中国石油化工股份有限公司 | 一种页岩油储层储渗空间定量表征及参数获取方法 |
| CN114428040B (zh) * | 2021-09-13 | 2022-07-22 | 中国石油化工股份有限公司 | 一种页岩油储层储渗空间定量表征及参数获取方法 |
| CN117494597A (zh) * | 2023-10-31 | 2024-02-02 | 中国石油大学(北京) | 一种三维数字岩心渗透率的确定方法、装置及设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110441209B (zh) | 2020-12-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110441209A (zh) | 一种基于致密储层数字岩心计算岩石渗透率的方法 | |
| CN108763711B (zh) | 一种基于岩心扫描图像分块数值模拟的渗透率预测方法 | |
| CN107817199B (zh) | 一种致密砂岩多尺度孔隙模型的构建方法和应用 | |
| CN104619952B (zh) | 具可靠多相渗透性决定数字岩石分析系统及方法 | |
| CN111428321B (zh) | 一种基于简化数字岩心的砾岩储层孔隙网络模型建模方法 | |
| RU2639727C2 (ru) | Системы и способы цифрового анализа горных пород с определением сэо многофазного потока | |
| CN111624147A (zh) | 岩心的相对渗透率测定方法及装置 | |
| CN113468829B (zh) | 基于孔隙网络模型的非稳态非牛顿两相流体驱替模拟方法 | |
| CN108802073A (zh) | 一种基于数字岩心的岩石电性参数获取方法及装置 | |
| CN104751473B (zh) | 碳酸岩多尺度孔洞特征的确定方法及装置 | |
| CN108267390B (zh) | 一种含纳米孔隙储层的气体渗透率确定方法 | |
| CN107655908A (zh) | 一种构建数字岩心的方法及装置 | |
| CN112082917B (zh) | 一种基于动态网络模拟的气水非稳态两相渗流模拟方法 | |
| CN114283254B (zh) | 基于核磁共振数据的岩心数字化孔隙网络模型构建方法 | |
| CN110320137A (zh) | 一种基于数字岩心的多尺度融合方法 | |
| CN110322427A (zh) | 一种基于三维成像技术的岩土体渗透率计算方法 | |
| Ji et al. | A multiscale reconstructing method for shale based on SEM image and experiment data | |
| CN104374679B (zh) | 一种确定储层渗透率贡献值的方法 | |
| CN114998246A (zh) | 基于图像分层处理技术计算多孔介质绝对渗透率的方法 | |
| CN110222368A (zh) | 一种利用二维切片计算岩心三维孔隙度和渗透率的方法 | |
| Wu et al. | Fractal characteristics of low-permeability sandstone reservoirs | |
| CN110598323B (zh) | 一种渗透破坏离散元模拟方法 | |
| CN105241911A (zh) | 基于lbm模拟低场核磁共振分析流体的方法及装置 | |
| CN111179296A (zh) | 一种基于数字岩心技术研究岩石导热特性的新方法 | |
| CN110441204B (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 |