CN112461728B - 一种全岩心渗透率及孔喉半径分布的确定方法及系统 - Google Patents
一种全岩心渗透率及孔喉半径分布的确定方法及系统 Download PDFInfo
- Publication number
- CN112461728B CN112461728B CN202011261788.3A CN202011261788A CN112461728B CN 112461728 B CN112461728 B CN 112461728B CN 202011261788 A CN202011261788 A CN 202011261788A CN 112461728 B CN112461728 B CN 112461728B
- Authority
- CN
- China
- Prior art keywords
- core
- permeability
- interpolation
- section
- pore throat
- 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.)
- Expired - Fee Related
Links
- 230000035699 permeability Effects 0.000 title claims abstract description 130
- 239000011148 porous material Substances 0.000 title claims abstract description 63
- 238000009826 distribution Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 43
- 239000011435 rock Substances 0.000 claims abstract description 49
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000005259 measurement Methods 0.000 claims description 13
- 239000011229 interlayer Substances 0.000 claims description 8
- 239000002356 single layer Substances 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 238000005520 cutting process Methods 0.000 claims description 5
- 238000003860 storage Methods 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000006872 improvement Effects 0.000 abstract description 2
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 9
- 229910052753 mercury Inorganic materials 0.000 description 9
- 239000011541 reaction mixture Substances 0.000 description 5
- 238000002591 computed tomography Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 238000002834 transmittance Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 230000007096 poisonous effect Effects 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 231100000331 toxic Toxicity 0.000 description 1
- 230000002588 toxic effect Effects 0.000 description 1
Images
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/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- 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
-
- 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
- G01N2015/0846—Investigating permeability, pore-volume, or surface area of porous materials by use of radiation, e.g. transmitted or reflected light
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种全岩心渗透率及孔喉半径分布的确定方法及系统，其技术方案为：对已知全岩心孔隙度的待测岩心进行剖面渗透率测量，采用反距离加权插值法在剖面内和剖面间插值，得到岩心的插值三维体；计算各点对应的岩心孔喉半径，并根据统计结果得到岩心孔喉半径分布；同时，根据渗透率值的统计结果能够计算出全岩心渗透率值。本发明操作简单、人为因素影响小、安全，节约成本，提高了效率。
Description
技术领域
本发明涉及石油开发技术领域，尤其涉及一种全岩心渗透率及孔喉半径分布的确定方法及系统。
背景技术
现在常有确定岩心孔喉半径分布的方法有两种，即CT扫描法和压汞法。CT扫描法又叫层析成像法，是发射X射线对岩心作旋转扫描，在每个位置可采集到一组一维的投影数据，再结合旋转运动，就可得到许多方向上的投影数据；综合这些投影数据，经过迭代运算就可以得到X射线衰减系数断面分布图，最后能得到岩心孔喉结构分布。
水银注入到被抽真空的岩心内需要克服岩石孔隙对其的毛管阻力，因此注入水银的每一点压力就代表一个相应孔隙大小的毛管压力，在这个压力下进入孔隙系统的水银量就代表这个相应的孔喉大小在这个系统中所联通的孔隙体积，然后根据压力、水银注入量绘制曲线即得到岩石毛管压力曲线。
发明人发现，上述两种方法均存在一定的缺陷，其中，CT扫描法的测量方法复杂，且费用较高；压汞法测量需要有毒的汞做介质，具有一定风险性，测量周期长，且对岩样具有永久破坏性。
发明内容
针对现有技术存在的不足，本发明的目的是提供一种全岩心渗透率及孔喉半径分布的确定方法及系统，操作简单、人为因素影响小、安全，节约成本，提高了效率。
为了实现上述目的，本发明是通过如下的技术方案来实现：
第一方面，本发明的实施例提供了一种全岩心渗透率及孔喉半径分布的确定方法，对已知全岩心孔隙度的待测岩心进行剖面渗透率测量，采用反距离加权插值法在剖面内和剖面间插值，得到岩心的插值三维体；计算各点对应的岩心孔喉半径，并根据统计结果得到岩心孔喉半径分布；同时，根据渗透率值的统计结果能够计算出全岩心渗透率值。
作为进一步的实现方式，包括以下步骤：
选取已知全岩心孔隙度的待测量岩心，将其沿纵向等间距切割成多个岩心剖面，以获得相应个数的切面；
将相邻两个切面上同一位置测得的岩心渗透率值求均值后即为岩心剖面对应的渗透率值；
以中间剖面建立岩心剖面坐标系，计算单层剖面纵向插值和层间剖面横向插值；
对所测得的渗透率和插值得到的渗透率值进行异常剔除，计算各点对应的岩心孔喉半径值；
对正常范围内的数据进行处理，得到岩心孔喉半径分布及全岩心渗透率值。
作为进一步的实现方式，将待测量岩心沿纵向等间距切割成三个岩心剖面；利用岩石扫描系统测量切面的渗透率。
作为进一步的实现方式，以中间岩心剖面的一个顶角作为坐标原点建立三维坐标系，沿岩心剖面长度方向为x轴、宽度方向为y轴，以垂直于岩心剖面方向为z轴。
作为进一步的实现方式，单层剖面纵向插值计算方法为：首先计算渗透率测量点到渗透率插值点的距离，之后利用任意相邻四个点进行反距离加权插值法计算插值点的渗透率值。
作为进一步的实现方式，层间剖面横向插值计算方法为：首先计算渗透率测量点到渗透率插值点的距离，之后利用相邻两剖面的八个点进行反距离加权插值法计算插值点的渗透率值。
作为进一步的实现方式，按照孔喉半径计算公式计算出各点对应的岩心孔喉半径值，然后对正常范围内的数据值和个数进行统计得到岩心孔喉半径的分布；对正常范围内渗透率的数据值和个数进行统计，然后利用加权平均求出其全岩心渗透率值。
第二方面，本发明实施例还提供了全岩心渗透率及孔喉半径分布的确定系统，包括：
剖面渗透率测量模块，被配置为对已知全岩心孔隙度的待测岩心进行剖面渗透率测量；
插值三维体获取模块，被配置为采用反距离加权插值法在剖面内和剖面间插值，得到岩心的插值三维体；
岩心孔喉半径分布获取模块，被配置为计算各点对应的岩心孔喉半径，并根据统计结果得到岩心孔喉半径分布；
全岩心渗透率值计算模块，被配置为根据渗透率值的统计结果计算全岩心渗透率值。
第三方面，本发明实施例还提供了一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的全岩心渗透率及孔喉半径分布的确定方法。
第四方面，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现所述的全岩心渗透率及孔喉半径分布的确定方法。
上述本发明的实施例的有益效果如下：
本发明的一个或多个实施方式对已知全岩心孔隙度的待测岩心进行剖面渗透率测量，采用反距离加权插值法在剖面内和剖面间插值，得到一个岩心的插值三维体，然后根据公式计算出各点对应的岩心孔喉半径，根据统计结果得到岩心孔喉半径分布；根据渗透率值的统计和个数还可计算出全岩心渗透率值；
相比于压汞法，解决了测量过程采用有毒汞做介质、时间周期长等问题，避免了产生安全风险；相比于CT扫描法，该方法操作简单，人为因素小，不仅节约了成本还提高了效率。
附图说明
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。
图1是本发明根据一个或多个实施方式的流程图；
图2是本发明根据一个或多个实施方式的全岩心剖面切割方式示意图；
图3是本发明根据一个或多个实施方式的岩心剖面渗透率测量点分布图；
图4是本发明根据一个或多个实施方式的岩心三维坐标系建立示意图；
图5是本发明根据一个或多个实施方式的单层剖面四点插值示意图；
图6是本发明根据一个或多个实施方式的层间剖面八点插值示意图。
具体实施方式
应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合；
为了方便叙述，本申请中如果出现“上”、“下”、“左”“右”字样，仅表示与附图本身的上、下、左、右方向一致，并不对结构起限定作用,仅仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的设备或元件必须具有特定的方位，以特定的方位构造和操作，因此不能理解为对本申请的限制。
岩心孔喉半径分布：沟通孔隙之间狭窄通道的大小在岩石上的分布；
渗透率：渗透率是指在一定压差下，岩石允许流体通过的能力。
实施例一：
本发明的实施例提供了一种全岩心渗透率及孔喉半径分布的确定方法，如图1所示，对已知全岩心孔隙度的待测岩心进行剖面渗透率测量，采用反距离加权插值法在剖面内和剖面间插值，得到岩心的插值三维体；计算各点对应的岩心孔喉半径，并根据统计结果得到岩心孔喉半径分布；同时，根据渗透率值的统计结果能够计算出全岩心渗透率值。
具体的，包括以下步骤：
步骤一：岩心剖面制取：
选取已知全岩心孔隙度的待测量岩心，如图2所示，将其沿纵向等间距切割成三个岩心剖面。
步骤二：岩心剖面测量：
本实施例采用美国AUTOSCAN-Ⅱ型岩石扫描系统进行渗透率测量，其定位精度1mm，渗透率测定范围0.1mD—3000mD。
如图3所示，待测量岩心切割后会有六个切面，将相邻两个切面上同一位置测得的岩心渗透率值求均值后即为剖面上对应的渗透率值，得到三个渗透率剖面。
步骤三：建立岩心剖面坐标系：
在本实施例中，切割形成矩形的岩心剖面，以中间岩心剖面左后方(以图4所示方向为参考)的顶角为坐标原点建立三维坐标系，沿岩心剖面长度方向为x轴，沿岩心剖面宽度方向为y轴，三个剖面左后方顶角在沿z轴分布。
步骤四：反距离加权插值：
(1)单层剖面纵向插值如图5所示：
设岩心剖面上有n个点,岩心剖面数为a＝1,2,3，坐标为
渗透率测点值为
其中i＝1,2,3…,n,待插值点坐标Q(xQ,yQ,zQ),插值点渗透率值HQ。
渗透率测量点
到渗透率插值点Q(xQ,yQ,zQ)的距离：
式(1)中，
表示剖面a上第i个点渗透率测量点横坐标；
表示剖面a上第i个点渗透率测量点纵坐标；xQ表示渗透率插值点横坐标；yQ表示渗透率插值点二维纵坐标；zQ表示渗透率插值点三维纵坐标；diQ表示渗透率插值点到已知点的距离。
利用相邻四个点进行反距离加权插值法计算插值点Q(xQ,yQ,zQ)的渗透率值：
式(2)中，
表示渗透率测量值，mD；diQ表示测量点到插值点之间的距离，m；HQ表示插值的透率值，mD。
根据插值点Q横纵坐标的变化可多次计算插值点渗透率的值。
(2)层间剖面横向插值如图6所示：
渗透率测量点
到渗透率插值点Q(xQ,yQ,zQ)的距离：
式(3)中，
表示剖面a上第i个点渗透率测量点横坐标；
表示剖面a上第i个点渗透率测量点纵坐标；
表示剖面a上第i个点渗透率测量点纵坐标；xQ表示渗透率插值点横坐标；yQ表示渗透率插值点纵坐标；zQ表示渗透率插值点横坐标；d′iQ表示渗透率插值点到剖面已知点的距离。
利用相邻两剖面的八个点进行反距离加权插值法计算插值点Q(xQ,yQ,zQ)的渗透率值HQ：
其中，(i≠3n)同时a≤2根据插值点Q横纵坐标的变化可多次计算插值点渗透率的值。
式(4)中，
表示渗透率测量值，mD；d′iQ表示剖面一测量点到插值点之间的距离，m；d″iQ表示剖面二测量点到插值点之间的距离，m；HQ表示插值的透率值，mD。
上述只对剖面一二进行插值，剖面二三之间原理相同，此处不再赘述。
步骤五：孔喉半径分布统计与全岩心渗透率值计算：
对所测得的渗透率和插值得到的渗透率值进行正常保留异常剔除，按照孔喉半径计算公式计算出各点对应的岩心孔喉半径值，然后对正常范围内的数据值和个数进行统计得到岩心孔喉半径的分布：
式(5)中，r表示孔喉平均半径，μm；K表示渗透率值；Ψ表示已知全岩心孔隙度，％。
对正常范围内渗透率的数据值和个数进行统计，然后利用加权平均求出其全岩心渗透率值：
式(6)中，
表示全岩心渗透率值，mD；kn表示渗透率值，mD；jn表示各渗透率值个数。
本实施例利用剖面渗透率插值方法研究孔喉半径分布，不仅操作简单、人为因素影响小、安全、还节约了成本，提高了效率。
实施例二：
本发明实施例还提供了全岩心渗透率及孔喉半径分布的确定系统，包括：
剖面渗透率测量模块，被配置为对已知全岩心孔隙度的待测岩心进行剖面渗透率测量；
插值三维体获取模块，被配置为采用反距离加权插值法在剖面内和剖面间插值，得到岩心的插值三维体；
岩心孔喉半径分布获取模块，被配置为计算各点对应的岩心孔喉半径，并根据统计结果得到岩心孔喉半径分布；
全岩心渗透率值计算模块，被配置为根据渗透率值的统计结果计算全岩心渗透率值。
实施例三：
本实施例还提供了一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现实施例一所述的全岩心渗透率及孔喉半径分布的确定方法。
实施例四：
本实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现实施例一所述的全岩心渗透率及孔喉半径分布的确定方法。
以上实施例二-四中涉及的各步骤与方法实施例一相对应，具体实施方式可参见实施例一的相关说明部分。术语“计算机可读存储介质”应该理解为包括一个或多个指令集的单个介质或多个介质；还应当被理解为包括任何介质，所述任何介质能够存储、编码或承载用于由处理器执行的指令集并使处理器执行本发明中的任一方法。
本领域技术人员应该明白，上述本发明的各模块或各步骤可以用通用的计算机装置来实现，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。本发明不限制于任何特定的硬件和软件的结合。
以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。
Claims (5)
1.一种全岩心渗透率及孔喉半径分布的确定方法，其特征在于，
对已知全岩心孔隙度的待测岩心进行剖面渗透率测量，采用反距离加权插值法在剖面内和剖面间插值，得到岩心的插值三维体；计算各点对应的岩心孔喉半径，并根据统计结果得到岩心孔喉半径分布；同时，根据渗透率值的统计结果能够计算出全岩心渗透率值；包括以下步骤：
选取已知全岩心孔隙度的待测量岩心，将其沿纵向等间距切割成多个岩心剖面，以获得相应个数的切面；以中间岩心剖面的一个顶角作为坐标原点建立三维坐标系，沿岩心剖面长度方向为x轴、宽度方向为y轴，以垂直于岩心剖面方向为z轴；
将相邻两个切面上同一位置测得的岩心渗透率值求均值后即为岩心剖面对应的渗透率值；
以中间剖面建立岩心剖面坐标系，计算单层剖面纵向插值和层间剖面横向插值；所述单层剖面纵向插值计算方法为：首先计算渗透率测量点到渗透率插值点的距离，之后利用任意相邻四个点进行反距离加权插值法计算插值点的渗透率值；所述层间剖面横向插值计算方法为：首先计算渗透率测量点到渗透率插值点的距离，之后利用相邻两剖面的八个点进行反距离加权插值法计算插值点的渗透率值；
对所测得的渗透率和插值得到的渗透率值进行异常剔除，按照孔喉半径计算公式计算出各点对应的岩心孔喉半径值，然后对正常范围内的数据值和个数进行统计得到岩心孔喉半径的分布；对正常范围内渗透率的数据值和个数进行统计，然后利用加权平均求出其全岩心渗透率值；
对正常范围内的数据进行处理，得到岩心孔喉半径分布及全岩心渗透率值。
2.根据权利要求1所述的一种全岩心渗透率及孔喉半径分布的确定方法，其特征在于，将待测量岩心沿纵向等间距切割成三个岩心剖面；利用岩石扫描系统测量切面的渗透率。
3.一种全岩心渗透率及孔喉半径分布的确定系统，其特征在于，包括：
剖面渗透率测量模块，被配置为对已知全岩心孔隙度的待测岩心进行剖面渗透率测量；
插值三维体获取模块，被配置为采用反距离加权插值法在剖面内和剖面间插值，得到岩心的插值三维体；
岩心孔喉半径分布获取模块，被配置为计算各点对应的岩心孔喉半径，并根据统计结果得到岩心孔喉半径分布；
全岩心渗透率值计算模块，被配置为根据渗透率值的统计结果计算全岩心渗透率值；
包括以下步骤：
选取已知全岩心孔隙度的待测量岩心，将其沿纵向等间距切割成多个岩心剖面，以获得相应个数的切面；以中间岩心剖面的一个顶角作为坐标原点建立三维坐标系，沿岩心剖面长度方向为x轴、宽度方向为y轴，以垂直于岩心剖面方向为z轴；
将相邻两个切面上同一位置测得的岩心渗透率值求均值后即为岩心剖面对应的渗透率值；
以中间剖面建立岩心剖面坐标系，计算单层剖面纵向插值和层间剖面横向插值；所述单层剖面纵向插值计算方法为：首先计算渗透率测量点到渗透率插值点的距离，之后利用任意相邻四个点进行反距离加权插值法计算插值点的渗透率值；所述层间剖面横向插值计算方法为：首先计算渗透率测量点到渗透率插值点的距离，之后利用相邻两剖面的八个点进行反距离加权插值法计算插值点的渗透率值；
对所测得的渗透率和插值得到的渗透率值进行异常剔除，按照孔喉半径计算公式计算出各点对应的岩心孔喉半径值，然后对正常范围内的数据值和个数进行统计得到岩心孔喉半径的分布；对正常范围内渗透率的数据值和个数进行统计，然后利用加权平均求出其全岩心渗透率值；
对正常范围内的数据进行处理，得到岩心孔喉半径分布及全岩心渗透率值。
4.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1-2任一项所述的全岩心渗透率及孔喉半径分布的确定方法。
5.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1-2任一项所述的全岩心渗透率及孔喉半径分布的确定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011261788.3A CN112461728B (zh) | 2020-11-12 | 2020-11-12 | 一种全岩心渗透率及孔喉半径分布的确定方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011261788.3A CN112461728B (zh) | 2020-11-12 | 2020-11-12 | 一种全岩心渗透率及孔喉半径分布的确定方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112461728A CN112461728A (zh) | 2021-03-09 |
| CN112461728B true CN112461728B (zh) | 2021-07-02 |
Family
ID=74825603
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011261788.3A Expired - Fee Related CN112461728B (zh) | 2020-11-12 | 2020-11-12 | 一种全岩心渗透率及孔喉半径分布的确定方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112461728B (zh) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103926267A (zh) * | 2014-04-28 | 2014-07-16 | 西安石油大学 | 一种定量评价应力敏感过程中孔喉变化程度的方法 |
| CN104809277A (zh) * | 2015-04-14 | 2015-07-29 | 中国石油天然气股份有限公司 | 一种超低渗致密储层水平井地质建模方法 |
| CN105184867A (zh) * | 2015-09-29 | 2015-12-23 | 中国矿业大学 | 基于规则格网和角点网格技术的三维地层建模方法 |
| CN105279790A (zh) * | 2014-06-13 | 2016-01-27 | 中国石油化工股份有限公司 | 裂缝网络三维数字岩心建模方法 |
| CN105974092A (zh) * | 2016-07-08 | 2016-09-28 | 重庆科技学院 | 一种致密储层孔喉全尺度表征及分析方法 |
| CN107219564A (zh) * | 2017-05-17 | 2017-09-29 | 中国石油天然气股份有限公司 | 一种处理垂向连井剖面的方法及装置 |
| CN110298105A (zh) * | 2019-06-26 | 2019-10-01 | 大连理工大学 | 饱和多孔介质大变形分析的ccpdi-impm方法 |
| KR102034868B1 (ko) * | 2018-12-29 | 2019-10-22 | 한국건설기술연구원 | 광역지하수 거동 해석 결과를 활용한 토류벽 벽체 및 지반변형 해석의 연계를 위한 수치해석 기법 |
| CN110593842A (zh) * | 2019-10-22 | 2019-12-20 | 西南石油大学 | 实验确定页岩储层水力压裂自支撑裂缝导流能力的方法 |
| CN110618158A (zh) * | 2019-10-28 | 2019-12-27 | 中石化石油工程技术服务有限公司 | 利用核磁共振信息构建岩心毛管压力曲线的方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110322427B (zh) * | 2019-04-02 | 2023-05-09 | 中国矿业大学 | 一种基于三维成像技术的岩土体渗透率计算方法 |
| CN111706396B (zh) * | 2020-06-12 | 2021-10-22 | 鞍钢集团矿业有限公司 | 一种采空区安全预警监测方法 |
-
2020
- 2020-11-12 CN CN202011261788.3A patent/CN112461728B/zh not_active Expired - Fee Related
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103926267A (zh) * | 2014-04-28 | 2014-07-16 | 西安石油大学 | 一种定量评价应力敏感过程中孔喉变化程度的方法 |
| CN105279790A (zh) * | 2014-06-13 | 2016-01-27 | 中国石油化工股份有限公司 | 裂缝网络三维数字岩心建模方法 |
| CN104809277A (zh) * | 2015-04-14 | 2015-07-29 | 中国石油天然气股份有限公司 | 一种超低渗致密储层水平井地质建模方法 |
| CN105184867A (zh) * | 2015-09-29 | 2015-12-23 | 中国矿业大学 | 基于规则格网和角点网格技术的三维地层建模方法 |
| CN105974092A (zh) * | 2016-07-08 | 2016-09-28 | 重庆科技学院 | 一种致密储层孔喉全尺度表征及分析方法 |
| CN107219564A (zh) * | 2017-05-17 | 2017-09-29 | 中国石油天然气股份有限公司 | 一种处理垂向连井剖面的方法及装置 |
| KR102034868B1 (ko) * | 2018-12-29 | 2019-10-22 | 한국건설기술연구원 | 광역지하수 거동 해석 결과를 활용한 토류벽 벽체 및 지반변형 해석의 연계를 위한 수치해석 기법 |
| CN110298105A (zh) * | 2019-06-26 | 2019-10-01 | 大连理工大学 | 饱和多孔介质大变形分析的ccpdi-impm方法 |
| CN110593842A (zh) * | 2019-10-22 | 2019-12-20 | 西南石油大学 | 实验确定页岩储层水力压裂自支撑裂缝导流能力的方法 |
| CN110618158A (zh) * | 2019-10-28 | 2019-12-27 | 中石化石油工程技术服务有限公司 | 利用核磁共振信息构建岩心毛管压力曲线的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 基于反距离加权插值法的产量劈分新方法;王立等;《断块油气田》;20180930;第25卷(第5期);第617-621页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112461728A (zh) | 2021-03-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hutchinson | Interpolation of rainfall data with thin plate smoothing splines. Part I: Two dimensional smoothing of data with short range correlation | |
| EP2153681B1 (en) | Method and arrangement for improved radio network planning, simulation and analyzing in telecommunications | |
| Martins et al. | Percolation on two-and three-dimensional lattices | |
| CN108717727A (zh) | 一种三维编织复合材料预制体结构的识别与建模方法 | |
| AU2011289788A1 (en) | Reservoir upscaling method with preserved transmissibility | |
| CN112461728B (zh) | 一种全岩心渗透率及孔喉半径分布的确定方法及系统 | |
| CN106019339B (zh) | 高精度gps定位点获取方法及系统 | |
| CN104183021B (zh) | 一种利用可移动空间网格精简点云数据的方法 | |
| Mertens et al. | Percolation in finite matching lattices | |
| CN108009281B (zh) | 基于矩形方格的渲染方法、装置及计算机可读存储介质 | |
| KR100965838B1 (ko) | 항공용 라이다로부터 추출된 빌딩폴리곤의 내재적기하정규화 방법 | |
| CN117310789B (zh) | 探测器通道响应线性校正方法、装置、设备和存储介质 | |
| CN117237435A (zh) | 肿瘤预后效果评估的方法、装置、电子设备和存储介质 | |
| Hinrichsen et al. | A fast algorithm for estimating large-scale permeabilities of correlated anisotropic media | |
| CN108226002B (zh) | 基于刻度电成像次生孔隙度谱截止值获取次生孔隙度的方法 | |
| CN107607129A (zh) | 数据校准方法、装置及电子设备 | |
| CN116029558A (zh) | 尾矿库安全风险预警方法、装置及电子设备 | |
| CN112184900A (zh) | 高程数据的确定方法、装置及存储介质 | |
| CN102682199B (zh) | 基于二维索引计算投影系数构造系统矩阵的简易方法 | |
| CN110120067A (zh) | 岩石自然结构面磨损体积确定方法、装置和计算机设备 | |
| CN112947858B (zh) | 一种raid 5校验值的更新方法、装置和介质 | |
| CN116050231A (zh) | 基于图形算法的离散元孔隙率计算与插值优化方法及系统 | |
| CN116049330A (zh) | 独立坐标系建立方法、装置、计算机设备、存储介质 | |
| CN115550190A (zh) | 网络拓扑图的确定方法、装置、设备及存储介质 | |
| CN110211199B (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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Wang Xingjian Inventor after: Liao Wanping Inventor after: Li Qingwu Inventor after: Wang Chongming Inventor after: Fan Caiwei Inventor before: Wang Xingjian Inventor before: Liao Wanping Inventor before: Li Qingwu Inventor before: Wang Chongming Inventor before: Fan Caiwei Inventor before: Gao Xilong |
|
| CB03 | Change of inventor or designer information | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210702 |