CN113504172A - 一种基于数字岩心的渗透性能评价方法及装置 - Google Patents
一种基于数字岩心的渗透性能评价方法及装置 Download PDFInfo
- Publication number
- CN113504172A CN113504172A CN202111058732.2A CN202111058732A CN113504172A CN 113504172 A CN113504172 A CN 113504172A CN 202111058732 A CN202111058732 A CN 202111058732A CN 113504172 A CN113504172 A CN 113504172A
- Authority
- CN
- China
- Prior art keywords
- evaluation
- permeability
- parameters
- parameter
- dimensionless
- 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.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 243
- 230000035699 permeability Effects 0.000 title claims abstract description 167
- 239000011435 rock Substances 0.000 claims abstract description 62
- 238000000034 method Methods 0.000 claims abstract description 35
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 10
- 238000009826 distribution Methods 0.000 claims description 7
- 238000011161 development Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 abstract description 7
- 230000007547 defect Effects 0.000 abstract description 2
- 239000011148 porous material Substances 0.000 description 28
- 238000005516 engineering process Methods 0.000 description 8
- 238000002591 computed tomography Methods 0.000 description 7
- 239000012530 fluid Substances 0.000 description 5
- 238000000691 measurement method Methods 0.000 description 5
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010835 comparative analysis Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000005325 percolation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 238000002834 transmittance Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004624 confocal microscopy Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011496 digital image analysis Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000001000 micrograph Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000004660 morphological change Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000036314 physical performance Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000003313 weakening 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
-
- 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
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/901—Indexing; Data structures therefor; Storage structures
- G06F16/9017—Indexing; Data structures therefor; Storage structures using directory or table look-up
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/903—Querying
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/906—Clustering; Classification
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- General Engineering & Computer Science (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Dispersion Chemistry (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Radiology & Medical Imaging (AREA)
- Computational Linguistics (AREA)
- Software Systems (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开一种基于数字岩心的渗透性能评价方法及装置。所述方法包括根据不同岩性类型，选取与渗透率相关性强的数字岩心定量参数；将评价参数中有量纲的参数转变为无量纲，形成无量纲量表，并为不同的评价参数设置不同权重系数，形成权重表；当评价某一岩石样本时，根据权重表和无量纲量表计算CT渗透率评价因子；根据CT渗透率评价因子系统查找预先创建的CT评价类别解释模板数据表，得到对应的评价类别。本发明结合各项参数创建CT渗透率评价因子，实现量化综合评价岩心孔喉连通性及渗流级别好坏程度，解决仅使用某一项而不是结合多项参数得出结论的评价方法不足，实现不同岩心样品之间渗透性能对比分析。
Description
技术领域
本发明涉及数字岩心微观孔隙结构扫描处理领域，尤其涉及一种基于数字岩心的渗透性能评价方法及装置。
背景技术
数字岩心，是岩心经过X射线扫描后的图像，通过计算机反演、建模、重构获取模拟岩心，其核心建模方法可分为两大类：物理实验方法和数值重建方法。物理实验方法均借助高倍光学显微镜、扫描电镜或CT成像仪等高精度仪器获取岩心的平面图像，之后对平面图像进行三维重建即可得到数字岩心。
CT全称是Computed Tomography，是计算机断层扫描技术，CT数字岩心应用在石油勘探开发地质储层岩心分析方面，它是岩心在X射线扫描下形成的数字图像分析与应用。其物理原理是基于射线与物质的相互作用，岩心内部不同密度与X射线响应程度不同，岩石骨架密度高X射线透过率低，岩心孔隙密度小X射线透过率高，通过对样品进行360°全方位扫描，获取一系列二维投影图像并进行三维重构，获得岩样三维模型，构建数字岩心。由于数字岩心能够完整、真实地反应岩心的微观结构，因此许多学者利用数字岩心技术进行渗流特征的研究，包括有利用整合法建立了碳酸盐岩双孔隙网络模型，对碳酸盐岩不同尺度孔隙的渗流规律进行了研究，有基于多孔介质的分形理论，对岩心孔隙度和渗透率进行了预测，证明了数字岩心在预测岩心物性方面的准确性，还有利用二维扫描电镜图像重构数字岩心，构建了基于数字岩心技术的页岩气藏与致密气藏表观渗透率计算模型，并对表观渗透率的影响因素进行了分析。然而关于利用数字岩心技术对储层渗透性的综合评价及分类方面的研究较少。
目前，岩心孔隙结构测量方法有许多种，分为间接测定法、直接测量法和数字岩心法。数字岩心法分为序列切片成像法、激光扫描共聚焦显微镜法和X射线CT扫描法。CT数字岩心孔隙类别评价方法属于X射线CT扫描法，在X射线CT扫描法中，各项检测参数常作为单一物性性能的指标参与评价，各单一参数（喉道半径、迂曲度、配位数）只能一定程度反映岩心渗透能力的优劣，但无法获取一个统一的参数进行岩心渗透性的评价，单一数据存在“孤岛”现象，没有统一的量化定义。
此外，CT数字岩心分析中，没有建立有效的渗透性评价分类方法，实现不同样品的对比分析。CT数字岩心属于新的储层评价技术，石油天然气行业标准中（SY/T 6285-2011，以下简称行业标准）各种类型划分依据间接测量法和直接测量法所得，因采用技术不同提供参数项也不同，CT数字岩心产生的参数无现行的CT行业岩心评价标准，其分析评价也只能按行业标准参照执行，因油气行业标准建立所用实验手段与CT扫描物理原理有较大差别，所以使用油气行业标准评价CT扫描岩心数据存在较大偏差，无法快速准确地对储层渗透性分类及评价。
发明内容
本发明提供了一种基于数字岩心的渗透性能评价方法，包括：
根据不同岩性类型，选取与渗透率相关性强的数字岩心定量参数，作为不同岩性类型的评价参数；
将评价参数中有量纲的参数转变为无量纲，形成无量纲量表，并为不同的评价参数设置不同权重系数，形成权重表，以赋予被评价对象不同侧面重要程度的定量分配；
当评价某一岩石样本时，根据权重表和无量纲量表计算CT渗透率评价因子；
根据CT渗透率评价因子系统查找预先创建的CT评价类别解释模板数据表，得到对应的评价类别，实现单颗岩心样品储层渗透性综合评价。
如上所述的基于数字岩心的渗透性能评价方法，其中，针对常规砂岩，采用平均喉道半径、平均配位数、迂曲度作为评价参数；针对微裂缝发育的致密型岩性，增添关于微裂缝的参数，具体采用平均喉道半径、平均配位数、迂曲度、平均微裂缝宽度、微裂缝贡献率作为评价参数。
如上所述的基于数字岩心的渗透性能评价方法，其中，计算CT渗透率评价因子具体包括如下子步骤：
当评价某一岩石样本时，扫描得到该样本的CT参数表；
根据CT参数表中的参数类别查找权重表和无量纲量表，得到该参数对应的权重数和无量纲量，根据查找到的权重数和无量纲量计算CT渗透率评价因子。
如上所述的基于数字岩心的渗透性能评价方法，其中，计算CT渗透率评价因子，具体包括：
为岩石样本的CT参数表中的每一个参数查询权重表和无量纲量表，得到每个参数对应的权重值和无量纲量；
运用下式计算CT渗透率评价因子：
其中，
为每颗岩石样品的CT渗透率评价因子，范围在0~1之间；
为第n个评价参数对应的权重值，
为第n个评价参数对应的无量纲量，n的取值为1到N，N为所有评价参数的总数。
如上所述的基于数字岩心的渗透性能评价方法，其中，预先为CT渗透率评价因子与评价类别创建相对应的CT评价类别解释模板数据表，所述CT评价类别解释模板数据表是由CT渗透率评价因子和渗透率评价标准级别相关关系所建立的，根据CT渗透率评价因子查找所述CT评价类别解释模板数据表得到对应的评价类别。
本发明还提供一种基于数字岩心的渗透性能评价装置，包括：
评价参数选定模块，用于根据不同岩性类型，选取与渗透率相关性强的数字岩心定量参数，作为不同岩性类型的评价参数；
权重表构建模块，用于为不同的评价参数设置不同权重系数，形成权重表，以赋予被评价对象不同侧面重要程度的定量分配；
无量纲量表构建模块，用于将评价参数中有量纲的参数转变为无量纲，形成无量纲量表；
CT渗透率评价因子计算模块，用于当评价某一岩石样本时，根据权重表和无量纲量表计算CT渗透率评价因子；
数字岩心储层渗透性类别评价模块，用于根据CT渗透率评价因子查找预先创建的CT评价类别解释模板数据表，得到对应的数字岩心储层渗透率评价类别，实现单颗岩心样品储层渗透性评价。
如上所述的基于数字岩心的渗透性能评价装置，其中，CT渗透率评价因子系数计算模块，具体用于当评价某一岩石样本时，扫描得到该样本的CT参数表；根据CT参数表中的参数类别查找权重表和无量纲量表，得到该参数对应的权重数和无量纲量，根据查找到的权重数和无量纲量计算CT渗透率评价因子。
如上所述的基于数字岩心的渗透性能评价装置，其中，CT渗透率评价因子系数计算模块，具体用于：
为岩石样本的CT参数表中的每一个参数查询权重表和无量纲量表，得到每个参数对应的权重值和无量纲量；运用下式计算CT渗透率评价因子：
其中，
为每颗岩石样品的CT渗透率评价因子，范围在0~1之间；
为第n个评价参数对应的权重值，
为第n个评价参数对应的无量纲量，n的取值为1到N，N为所有评价参数的总数。
如上所述的基于数字岩心的渗透性能评价装置，其中，数字岩心储层渗透性类别评价模块，还用于预先为CT渗透率评价因子与评价类别创建相对应的CT评价类别解释模板数据表，所述CT评价类别解释模板数据表是由CT渗透率评价因子和渗透率评价标准级别相关关系所建立的，根据CT渗透率评价因子查找所述CT评价类别解释模板数据表得到对应的评价类别。
本发明还提供一种计算机可读存储介质，其特征在于，包括至少一个存储器和至少一个处理器；
存储器用于存储一个或多个程序指令；
处理器，用于运行一个或多个程序指令，用以执行上述任一项所述的一种基于数字岩心的渗透性能评价方法。
本发明实现的有益效果如下：本发明结合各项参数创建CT渗透率评价因子，即反映每颗样品微观渗透率的评价系数，实现量化综合评价岩心孔喉连通性及渗流级别好坏程度，解决仅使用某一项而不是结合多项参数得出结论的评价方法不足，实现不同岩心样品之间渗透性能对比分析，填补行业内CT扫描数字岩心技术综合评价岩心渗透性方法的空白。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，还可以根据这些附图获得其他的附图。
图1是本发明实施例一提供的一种基于数字岩心的渗透性能评价方法流程图；
图2是计算CT渗透率评价因子的具体操作流程图；
图3是CT渗透率评价因子与实测渗透率数据相关性图版示例图。
具体实施方式
下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在介绍本发明的数字岩心微观孔隙结构类别评价方法前，先对岩石的各属性参数进行解释说明，围绕岩石所具有的喉道的几何形状、大小、分布、相互连通情况、迂曲度特征，以及孔隙与喉道间的配置关系等，来反映储层中渗流通道的形态特征；
优选采用5微米、8微米精度扫描分析岩心，产生了孔隙度、有效孔隙度、平均孔隙半径、平均喉道半径、孔喉比、平均配位数、迂曲度、裂缝宽度、孔隙表面积、孔隙比表面积、高密度矿物含量等众多参数，选定上述若干数据作为评价参数；例如，本发明选用平均喉道半径、平均配位数、迂曲度作为评价参数，以下对各参数意义描述如下：
喉道半径：岩石的孔隙系统由孔隙和喉道两部分组成。孔隙为系统中的膨大部分，连通孔隙的细小部分称为喉道。孔隙是流体赋存于岩石中的基本储集空间，而喉道则是控制流体在岩石中渗流的重要通道。流体在自然界复杂的孔隙系统中流动时，都要经历一系列的交替着的孔隙和喉道。喉道半径是以能够通过孔隙喉道的最大球体半径来衡量的，单位是微米(μm)。若喉道半径大，孔隙空间的连通性好，液体在孔隙系统中的渗流能力就强。
配位数：配位数是每个孔隙所连通喉道的个数。配位数的大小不仅对流体的渗流起重要的作用，而且对于油气的采出、驱油效率也有重要的控制作用。
迂曲度：迂曲度是刻画喉道沿流体流动方向上形态变化的参数，体现了渗流通道的复杂程度，迂曲度最小值为1，代表理想化的直线型喉道，迂曲度越大代表渗流通道越复杂，越不利于渗流；迂曲度越接近1代表渗流通道越平直，越有利于渗流。
基于上述岩石属性参数，本发明以下实施例详细描述数字岩心渗透性的评价方法。
实施例一
如图1所示，本发明实施例一提供一种基于数字岩心的渗透性能评价方法，包括：
步骤110、根据不同岩性类型，选取与渗透率相关性强的数字岩心定量参数，作为不同岩性类型的评价参数；
具体地，根据所评价的不同种类的岩石岩性微观孔隙结构差别及其产生的参数数值特点选取与渗透率相关性强的数字岩心定量参数：例如针对常规砂岩，优选采用平均喉道半径、平均配位数、迂曲度等作为评价参数；针对微裂缝发育的致密型岩性，要增添关于微裂缝的参数，优选采用平均喉道半径、平均配位数、迂曲度、平均微裂缝宽度、微裂缝贡献率等作为评价参数。
步骤120、将评价参数中有量纲的参数转变为无量纲，形成无量纲量表，并为不同的评价参数设置不同权重系数，形成权重表，以赋予被评价对象不同侧面重要程度的定量分配；
本发明优选利用灰色关联法衡量不同评价参数之间的关联关系，为不同评价参数设置对应的权重系数，对于某一岩性类型的所有评价参数的权重之和为1，灰色关联法是根据不同的评价参数之间发展趋势的相似或相异程度，作为衡量参数间关联程度的一种方法；
灰色关联法的具体构造方法：
为了从数据信息的内部结构上分析被评判事物与其影响因素之间的关系，必须用某种数量指标定量反映被评判事件的性质。这种按一定顺序排列的数量指标，称为关联分析的母序列，记为
；子序列是决定或影响被评判事物性质的各子因素数据的有序排列，考虑主因素的m个子因素，则有子
；
由于不同储层微观孔隙结构特征参数的量纲、数值均存在较大差异，每种类型中不同评价参数的物理量有不同的数据范围和单位，为了使各参数具有可比性，首先需要采用极差化法对各参数进行无量纲化处理，即
，
为第i个评价参数的无量纲量，
为第i个评价参数的物理量，
为第i个评价参数的极大值，
为第i个评价参数的极小值；
然后计算同一观测时刻子因素与母因素之间的绝对差值，记为：
同一观测时刻各子因素与母因素之间的绝对差值的最大值为：
母序列与子序列的关联系数为：
其中，
为关联系数；
为分辨系数，
，分辨系数作用为削弱最大绝对差数值太大而失真的影响，提高关联系数之间的差异显著性；
各子因素对母因素之间的关联度为：
根据关联度，经归一化处理可得到权系数为：
以下以常规砂岩为例，针对常规砂岩，利用灰色关联法设定各参数间的权重系数，创建如表1所示：
表1 常规砂岩权重表
将表1中的各评价参数转变为如表2所示的无量纲量表：
表2 无量纲量表
步骤130、当评价某一岩石样本时，根据权重表和无量纲量表计算CT渗透率评价因子；
本发明实施例中，如图2所示，计算CT渗透率评价因子具体包括如下子步骤：
步骤210、当评价某一岩石样本时，扫描得到该样本的CT参数表；
例如，扫描某一岩石样本的CT分析主要参数如下表3所示：
表3 某岩石样本CT分析主要参数表
步骤220、根据CT参数表中的参数类别查找权重表和无量纲量表，得到该参数对应的权重数和无量纲量，根据查找到的权重数和无量纲量计算CT渗透率评价因子；
其中，计算CT渗透率评价因子，具体包括：
①为岩石样本的CT参数表中的每一个参数查询权重表和无量纲量表，得到每个参数对应的权重值和无量纲量；
当选择岩石样本的平均喉道半径参数项，表3中平均喉道半径值为14.07μm，查表1“常规砂岩权重表”确定参数名称为平均喉道半径所对应的权重数M1为0.342，对应的无量纲量为0.268219833；
当选择岩石样本的平均配位数参数项，表3中平均配位数值为0.4843，查表1“常规砂岩权重表”确定参数名称为平均配位数对应的权重数M2为0.416，对应的无量纲量为0.129715084；
当选择岩石样本的迂曲度参数项，表3中迂曲度值为2.97，查表1“常规砂岩权重表”确定参数名称为迂曲度对应的权重数M3为0.242，对应的无量纲量为0.820261438。
②在上述查表得到权重值和对应的无量纲量之后，运用下式计算CT渗透率评价因子：
其中，
为每颗岩石样品的CT渗透率评价因子，范围在0~1之间；
为第n个评价参数对应的权重值，
为第n个评价参数对应的无量纲量，n的取值为1到N，N为所有评价参数的总数；
本发明中常规砂岩的CT渗透率评价因子：
返回参见图1，步骤140、根据CT渗透率评价因子系统查找预先创建的CT评价类别解释模板数据表，得到对应的评价类别，实现单颗岩心样品储层渗透性综合评价；
预先为CT渗透率评价因子与评价类别创建相对应的CT评价类别解释模板数据表，所述CT评价类别解释模板数据表是由CT渗透率评价因子和渗透率评价标准级别相关关系所建立的，根据CT渗透率评价因子查找所述CT评价类别解释模板数据表得到对应的评价类别；
例如，如下表4所示的CT渗透率评价因子数据：
表4 CT渗透率评价级别解释模板数据表
表4中建立了四个区间一类、二类、三类、四类CT类别解释级别，实现对样品综合评价，其中使用数字岩心技术评价渗透性中：一类对应油气行业高渗及以上储层（>500mD）、二类对应中渗储层（50mD~500mD）、三类对应低渗及特低渗储层（1mD~50mD）、四类对应超低渗储层（<1mD）；
在步骤130示例中计算得到CT渗透率评价因子为0.3442，查询“CT评价类别解释模板数据表”，得到本示例岩石样品CT渗透性评价类别为四类。
实施例二
本发明实施例二提供一种基于数字岩心的渗透性能评价装置，包括：
评价参数选定模块，用于根据不同岩性类型，选取与渗透率相关性强的数字岩心定量参数，作为不同岩性类型的评价参数；
具体地，针对常规砂岩，采用平均喉道半径、平均配位数、迂曲度作为评价参数；针对微裂缝发育的致密型岩性，增添关于微裂缝的参数，具体采用平均喉道半径、平均配位数、迂曲度、平均微裂缝宽度、微裂缝贡献率作为评价参数。
权重表构建模块，用于为不同的评价参数设置不同权重系数，形成权重表，以赋予被评价对象不同侧面重要程度的定量分配；
无量纲量表构建模块，用于将评价参数中有量纲的参数转变为无量纲，形成无量纲量表；
CT渗透率评价因子计算模块，用于当评价某一岩石样本时，根据权重表和无量纲量表计算CT渗透率评价因子；
其中，CT渗透率评价因子系数计算模块，具体用于当评价某一岩石样本时，扫描得到该样本的CT参数表；根据CT参数表中的参数类别查找权重表和无量纲量表，得到该参数对应的权重数和无量纲量，根据查找到的权重数和无量纲量计算CT渗透率评价因子。
具体地，CT渗透率评价因子系数计算模块，具体用于：
为岩石样本的CT参数表中的每一个参数查询权重表和无量纲量表，得到每个参数对应的权重值和无量纲量；运用下式计算CT渗透率评价因子：
其中，
为每颗岩石样品的CT渗透率评价因子，范围在0~1之间；
为第n个评价参数对应的权重值，
为第n个评价参数对应的无量纲量，n的取值为1到N，N为所有评价参数的总数。
数字岩心储层渗透性类别评价模块，用于根据CT渗透率评价因子查找预先创建的CT评价类别解释模板数据表，得到对应的数字岩心储层渗透率评价类别，实现单颗岩心样品储层渗透性评价。
其中，数字岩心储层渗透性类别评价模块，还用于预先为CT渗透率评价因子与评价类别创建相对应的CT评价类别解释模板数据表，所述CT评价类别解释模板数据表是由CT渗透率评价因子和渗透率评价标准级别相关关系所建立的，根据CT渗透率评价因子查找所述CT评价类别解释模板数据表得到对应的评价类别。
以下为本发明技术方案在实际应用中进行数字岩心渗透性类别评价的实例：
在某油田应用中，取样点数据26颗进行岩石渗透性评价分类。通过CT扫描参数分析得到一系列岩石微观孔喉结构定量参数，其中包括计算CT渗透率评价因子所需的三个参数：平均喉道半径、平均配位数、迂曲度，利用极差化方法计算各样品点对应的无量纲量详见下表5：
表5 样品点渗透性评价分类表
根据表5数据建立如图3所示的CT渗透率评价因子与实测渗透率数据相关性图版，证明本发明提出的评价因子参数与岩心渗透性相关性强，在评价渗透性方面具有很好的效果。将本发明对渗透性的评价类别与常规评价渗透率级别相比对，其中一类对应高渗符合率为100%；二类对应中渗符合率为100%，三类对应低渗及特低渗符合率为75%；四类对应超低渗符合率为73.33%。数据分析证明，该方法评价岩心渗透性整体符合率较高，能实现岩心快速渗透性评价。
与上述实施例对应的，本发明实施例提供一种计算机可读存储介质，包括：至少一个存储器和至少一个处理器；
存储器用于存储一个或多个程序指令；
处理器，用于运行一个或多个程序指令，用以执行一种基于数字岩心的渗透性能评价方法。
本发明所公开的实施例提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机程序指令，当所述计算机程序指令在计算机上运行时，使得计算机执行一种基于数字岩心的渗透性能评价方法。
在本发明实施例中，处理器可以是一种集成电路芯片，具有信号的处理能力。处理器可以是通用处理器、数字信号处理器(Digital Signal Processor，简称DSP)、专用集成电路(Application Specific Integrated Circuit，简称ASIC)、现场可编程门阵列(FieldProgrammable Gate Array，简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。
可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。处理器读取存储介质中的信息，结合其硬件完成上述方法的步骤。
存储介质可以是存储器，例如可以是易失性存储器或非易失性存储器，或可包括易失性和非易失性存储器两者。
其中，非易失性存储器可以是只读存储器(Read-Only Memory，简称ROM)、可编程只读存储器(Programmable ROM，简称PROM)、可擦除可编程只读存储器(Erasable PROM，简称EPROM)、电可擦除可编程只读存储器(Electrically EPROM，简称EEPROM)或闪存。
易失性存储器可以是随机存取存储器(Random Access Memory，简称RAM)，其用作外部高速缓存。通过示例性但不是限制性说明，许多形式的RAM可用，例如静态随机存取存储器(Static RAM，简称SRAM)、动态随机存取存储器(Dynamic RAM，简称DRAM)、同步动态随机存取存储器(Synchronous DRAM，简称SDRAM)、双倍数据速率同步动态随机存取存储器(Double Data RateSDRAM，简称DDRSDRAM)、增强型同步动态随机存取存储器(EnhancedSDRAM，简称ESDRAM)、同步连接动态随机存取存储器(Synchlink DRAM，简称SLDRAM)和直接内存总线随机存取存储器(DirectRambus RAM，简称DRRAM)。
本发明实施例描述的存储介质旨在包括但不限于这些和任意其它适合类型的存储器。
本领域技术人员应该可以意识到，在上述一个或多个示例中，本发明所描述的功能可以用硬件与软件组合来实现。当应用软件时，可以将相应功能存储在计算机可读介质中或者作为计算机可读介质上的一个或多个指令或代码进行传输。计算机可读介质包括计算机存储介质和通信介质，其中通信介质包括便于从一个地方向另一个地方传送计算机程序的任何介质。存储介质可以是通用或专用计算机能够存取的任何可用介质。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的技术方案的基础之上，所做的任何修改、等同替换、改进等，均应包括在本发明的保护范围之内。
Claims (10)
1.一种基于数字岩心的渗透性能评价方法，其特征在于，包括：
根据不同岩性类型，选取与渗透率相关性强的数字岩心定量参数，作为不同岩性类型的评价参数；
将评价参数中有量纲的参数转变为无量纲，形成无量纲量表，并为不同的评价参数设置不同权重系数，形成权重表，以赋予被评价对象不同侧面重要程度的定量分配；
当评价某一岩石样本时，根据权重表和无量纲量表计算CT渗透率评价因子；
根据CT渗透率评价因子系统查找预先创建的CT评价类别解释模板数据表，得到对应的评价类别，实现单颗岩心样品储层渗透性综合评价。
2.如权利要求1所述的基于数字岩心的渗透性能评价方法，其特征在于，针对常规砂岩，采用平均喉道半径、平均配位数、迂曲度作为评价参数；针对微裂缝发育的致密型岩性，增添关于微裂缝的参数，具体采用平均喉道半径、平均配位数、迂曲度、平均微裂缝宽度、微裂缝贡献率作为评价参数。
3.如权利要求1所述的基于数字岩心的渗透性能评价方法，其特征在于，计算CT渗透率评价因子具体包括如下子步骤：
当评价某一岩石样本时，扫描得到该样本的CT参数表；
根据CT参数表中的参数类别查找权重表和无量纲量表，得到该参数对应的权重数和无量纲量，根据查找到的权重数和无量纲量计算CT渗透率评价因子。
4.如权利要求3所述的基于数字岩心的渗透性能评价方法，其特征在于，计算CT渗透率评价因子，具体包括：
为岩石样本的CT参数表中的每一个参数查询权重表和无量纲量表，得到每个参数对应的权重值和无量纲量；
运用下式计算CT渗透率评价因子：
其中， 
为每颗岩石样品的CT渗透率评价因子，范围在0~1之间；
为第n个评价参数对应的权重值，
为第n个评价参数对应的无量纲量，n的取值为1到N，N为所有评价参数的总数。
5.如权利要求1所述的基于数字岩心的渗透性能评价方法，其特征在于，预先为CT渗透率评价因子与评价类别创建相对应的CT评价类别解释模板数据表，所述CT评价类别解释模板数据表是由CT渗透率评价因子和渗透率评价标准级别相关关系所建立的，根据CT渗透率评价因子查找所述CT评价类别解释模板数据表得到对应的评价类别。
6.一种基于数字岩心的渗透性能评价装置，其特征在于，包括：
评价参数选定模块，用于根据不同岩性类型，选取与渗透率相关性强的数字岩心定量参数，作为不同岩性类型的评价参数；
权重表构建模块，用于为不同的评价参数设置不同权重系数，形成权重表，以赋予被评价对象不同侧面重要程度的定量分配；
无量纲量表构建模块，用于将评价参数中有量纲的参数转变为无量纲，形成无量纲量表；
CT渗透率评价因子计算模块，用于当评价某一岩石样本时，根据权重表和无量纲量表计算CT渗透率评价因子；
数字岩心储层渗透性类别评价模块，用于根据CT渗透率评价因子查找预先创建的CT评价类别解释模板数据表，得到对应的数字岩心储层渗透率评价类别，实现单颗岩心样品储层渗透性评价。
7.如权利要求6所述的基于数字岩心的渗透性能评价装置，其特征在于，CT渗透率评价因子系数计算模块，具体用于当评价某一岩石样本时，扫描得到该样本的CT参数表；根据CT参数表中的参数类别查找权重表和无量纲量表，得到该参数对应的权重数和无量纲量，根据查找到的权重数和无量纲量计算CT渗透率评价因子。
8.如权利要求7所述的基于数字岩心的渗透性能评价装置，其特征在于，CT渗透率评价因子系数计算模块，具体用于：
为岩石样本的CT参数表中的每一个参数查询权重表和无量纲量表，得到每个参数对应的权重值和无量纲量；运用下式计算CT渗透率评价因子：
其中，
为每颗岩石样品的CT渗透率评价因子，范围在0~1之间；
为第n个评价参数对应的权重值，
为第n个评价参数对应的无量纲量，n的取值为1到N，N为所有评价参数的总数。
9.如权利要求6所述的基于数字岩心的渗透性能评价装置，其特征在于，数字岩心储层渗透性类别评价模块，还用于预先为CT渗透率评价因子与评价类别创建相对应的CT评价类别解释模板数据表，所述CT评价类别解释模板数据表是由CT渗透率评价因子和渗透率评价标准级别相关关系所建立的，根据CT渗透率评价因子查找所述CT评价类别解释模板数据表得到对应的评价类别。
10.一种计算机可读存储介质，其特征在于，包括至少一个存储器和至少一个处理器；
存储器用于存储一个或多个程序指令；
处理器，用于运行一个或多个程序指令，用以执行如权利要求1-5任一项所述的一种基于数字岩心的渗透性能评价方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111058732.2A CN113504172A (zh) | 2021-09-10 | 2021-09-10 | 一种基于数字岩心的渗透性能评价方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111058732.2A CN113504172A (zh) | 2021-09-10 | 2021-09-10 | 一种基于数字岩心的渗透性能评价方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113504172A true CN113504172A (zh) | 2021-10-15 |
Family
ID=78016620
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111058732.2A Pending CN113504172A (zh) | 2021-09-10 | 2021-09-10 | 一种基于数字岩心的渗透性能评价方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113504172A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116431567A (zh) * | 2023-06-12 | 2023-07-14 | 煤炭科学研究总院有限公司 | 基于数字岩石力学的煤岩数据采集、管理与评估方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6792354B1 (en) * | 2000-11-13 | 2004-09-14 | O'meara, Jr. Daniel J. | Method for determining reservoir fluid volumes, fluid contacts, compartmentalization, and permeability in geological subsurface models |
| CN110887772A (zh) * | 2018-09-11 | 2020-03-17 | 中国石油化工股份有限公司 | 一种碳酸盐岩储层渗透率识别方法、系统及装置 |
| CN111537419A (zh) * | 2020-04-28 | 2020-08-14 | 北京润泽创新科技有限公司 | 一种数字岩心微观孔隙结构类别评价方法及装置 |
| CN111598440A (zh) * | 2020-05-14 | 2020-08-28 | 中国海洋石油集团有限公司 | 多角度驱动的复杂介质储层渗透率定量评价方法及系统 |
| WO2020215524A1 (zh) * | 2019-04-24 | 2020-10-29 | 山东科技大学 | 一种基于分形理论及ct扫描的煤体孔渗参数预测方法 |
| CN112855109A (zh) * | 2020-12-31 | 2021-05-28 | 西南石油大学 | 基于灰色关联法和层次分析法的压裂酸化选井选层方法 |
| CN113158490A (zh) * | 2021-05-13 | 2021-07-23 | 中国石油大学(华东) | 致密砂岩储层渗透率计算模型的建立方法 |
| CN113189122A (zh) * | 2021-05-13 | 2021-07-30 | 中海石油(中国)有限公司海南分公司 | 一种射孔损伤室内综合评价方法 |
-
2021
- 2021-09-10 CN CN202111058732.2A patent/CN113504172A/zh active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6792354B1 (en) * | 2000-11-13 | 2004-09-14 | O'meara, Jr. Daniel J. | Method for determining reservoir fluid volumes, fluid contacts, compartmentalization, and permeability in geological subsurface models |
| CN110887772A (zh) * | 2018-09-11 | 2020-03-17 | 中国石油化工股份有限公司 | 一种碳酸盐岩储层渗透率识别方法、系统及装置 |
| WO2020215524A1 (zh) * | 2019-04-24 | 2020-10-29 | 山东科技大学 | 一种基于分形理论及ct扫描的煤体孔渗参数预测方法 |
| CN111537419A (zh) * | 2020-04-28 | 2020-08-14 | 北京润泽创新科技有限公司 | 一种数字岩心微观孔隙结构类别评价方法及装置 |
| CN111598440A (zh) * | 2020-05-14 | 2020-08-28 | 中国海洋石油集团有限公司 | 多角度驱动的复杂介质储层渗透率定量评价方法及系统 |
| CN112855109A (zh) * | 2020-12-31 | 2021-05-28 | 西南石油大学 | 基于灰色关联法和层次分析法的压裂酸化选井选层方法 |
| CN113158490A (zh) * | 2021-05-13 | 2021-07-23 | 中国石油大学(华东) | 致密砂岩储层渗透率计算模型的建立方法 |
| CN113189122A (zh) * | 2021-05-13 | 2021-07-30 | 中海石油(中国)有限公司海南分公司 | 一种射孔损伤室内综合评价方法 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116431567A (zh) * | 2023-06-12 | 2023-07-14 | 煤炭科学研究总院有限公司 | 基于数字岩石力学的煤岩数据采集、管理与评估方法 |
| CN116431567B (zh) * | 2023-06-12 | 2023-09-29 | 煤炭科学研究总院有限公司 | 基于数字岩石力学的煤岩数据采集、管理与评估方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2011258594B2 (en) | Method for obtaining consistent and integrated physical properties of porous media | |
| Peth et al. | Three‐dimensional quantification of intra‐aggregate pore‐space features using synchrotron‐radiation‐based microtomography | |
| US20150331145A1 (en) | Method for producing a three-dimensional characteristic model of a porous material sample for analysis of permeability characteristics | |
| CN113552146B (zh) | 一种基于数字岩心技术的储层评价方法及装置 | |
| Moon et al. | Statistical inference over persistent homology predicts fluid flow in porous media | |
| Claes et al. | A three-dimensional classification for mathematical pore shape description in complex carbonate reservoir rocks | |
| CN111537419A (zh) | 一种数字岩心微观孔隙结构类别评价方法及装置 | |
| CN113504172A (zh) | 一种基于数字岩心的渗透性能评价方法及装置 | |
| CN113435066B (zh) | 一种基于数字岩心技术的测井解释储层评价方法 | |
| Pavičić et al. | Geometric and fractal characterization of pore systems in the Upper Triassic dolomites based on image processing techniques (example from Žumberak Mts, NW Croatia) | |
| Zhao et al. | Research on micro-pore structure and 3D visual characterization of inter-salt shale based on X-CT imaging digital core technology | |
| Eriavbe et al. | Machine learning application to permeability prediction using log & core measurements: A realistic workflow application for reservoir characterization | |
| Hou et al. | Characterization and evaluation of carbonate reservoir pore structure based on machine learning | |
| Xie et al. | Investigation of the Oriented Structure Characteristics of Shale Using Fractal and Structural Entropy Theory | |
| Stadtműller et al. | Characterization of the carbonate formation fracture system based on well logging data and results of laboratory measurements | |
| Tian et al. | Pore microstructure and multifractal characterization of lacustrine oil-prone shale using high-resolution SEM: a case sample from natural qingshankou shale | |
| CN115905917A (zh) | 综合静动态数据构建海域低渗气层分类曲线的方法 | |
| CN113047830B (zh) | 储层的流动单元识别方法、装置、计算机及存储介质 | |
| Schembre-McCabe et al. | A framework to validate digital rock technology | |
| Liu et al. | Quantitative characterization of tight rock microstructure of digital core | |
| CN112231881A (zh) | 基质渗透率计算方法及系统 | |
| Zhang et al. | Shale Oil–Water Two-Phase Flow Simulation Based on Pore Network Modeling | |
| Zhao et al. | Quantitative Classification and Prediction of Diagenetic Facies in Tight Gas Sandstone Reservoirs via Unsupervised and Supervised Machine Learning Models: Ledong Area, Yinggehai Basin | |
| Krakowska-Madejska et al. | Development of a Permeability Formula for Tight and Shale Gas Reservoirs Based on Advanced High-Precision Lab Measurement Techniques | |
| Caja et al. | Digital Rock Physics in Cuttings Using High-Resolution Thin Section Scan Images |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211015 |
|
| RJ01 | Rejection of invention patent application after publication |