CN109118552B - 一种数字岩心图像的孔喉连通区域识别方法 - Google Patents
一种数字岩心图像的孔喉连通区域识别方法 Download PDFInfo
- Publication number
- CN109118552B CN109118552B CN201710491368.6A CN201710491368A CN109118552B CN 109118552 B CN109118552 B CN 109118552B CN 201710491368 A CN201710491368 A CN 201710491368A CN 109118552 B CN109118552 B CN 109118552B
- Authority
- CN
- China
- Prior art keywords
- pixel
- pore
- matrix
- adjacent
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T9/00—Image coding
- G06T9/20—Contour coding, e.g. using detection of edges
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
Abstract
本申请提出了一种数字岩心图像的孔喉连通区域识别方法,包括以下步骤:S1、构建像素矩阵,并对像素矩阵中待识别的孔喉像素点进行编号,得到编号矩阵;S2、构建相邻矩阵;S3、对所述相邻矩阵进行DM分解,得到像素编号排列向量和像素编号分类向量;S4、根据像素编号排列向量和像素编号分类向量识别数字原始岩心图像中的连通区域。本发明仅通过一次运算,即可计算得到数字岩心中所有的孔喉连通区域,计算复杂程度低,数据存储量小,运算速度快。能够正确的对数字岩心中的连通孔喉进行识别。
Description
技术领域
本发明涉及油田开发技术领域,尤其是涉及一种数字岩心图像的孔喉连通区域识别方法。
背景技术
岩心的三维数字图像由CT或扫描电镜等实验手段得到,图像中通常包含了岩心内部连通和非连通的全部孔喉空间。在进行流动模拟时,采用的数字岩心中应仅包含参与流动的孔喉空间,即与流动通道入口和出口均相连的连续孔喉空间。对于孤立的孔喉空间(死孔隙),应在流动模拟之前进行删除,以减少流动计算的工作量。
假设数字岩心中0代表骨架,1代表孔喉。目前常用有两种方法来判断孔喉的连通性。
第一种算法称为“洪水算法”,具体是:①从岩心中选取值为1的像素点a作为种子,将种子值设为2;②判断与a相邻的26个像素是否为1。如果为1,则认为该像素b代表的孔喉与种子相连;③将b作为新的种子(值设置为2);④重复②~③直到周围再也找不到值为1的点。此时,仅且全部与种子a连通的点被标记为2。
第二种算法称为“扫描线算法”,具体是:①从岩心中选取值为1的像素作为种子a,种子值标记为2;②以a为中心沿x轴向两端查找边界,得到一条通过a的扫描线La,将扫描线上的点的值标记为2。③在y方向上与La相邻的像素中找到值为1的点,作为新的种子b,重复②得到新的扫描线Lb。④重复3得到新的种子c,d…和相应的扫描线Lc,Ld…;⑤以b,c,d…作为新的种子重复①~④,直到在第③步中找不到任何一个与扫描线相邻的值为1的点。此时,仅且全部与种子a连通的点被标记为2。
在得到孔喉连通性的基础上,判断该连通区域是否与岩心流动的入口和出口相连。有以下两种方法:一是将种子布置在岩心的入口端,并判断该连通区域内是否有种子位于岩心出口端;二是直接求取该连通区域内所有点中是否同时包含位于入口和出口的点。
但是上述两种方法均存在一定问题。“洪水算法”查找连通孔喉空间有两个主要缺点:一是运算效率低,每次都要搜索种子周围的26个像素,运算量大;二是需要的堆栈空间大,由于是递归算法,压入堆栈的种子数量很大。“扫描线算法”所用到的堆栈要小于洪水算法,但是由于要查找扫描线周围的点,运算效率仍然不高。
基于计算方式的不同发展出了不同的识别方法,基于DM(Dulmage-Mendelsohn)分解的图像识别方法是将DM分解应用在数字岩心图像所对应的相邻矩阵上,通过DM分解对相邻矩阵进行分析,从而得到数字岩心图像的孔喉连通区域。DM算法是1958年由A.L.Dulmage和Nathan Mendelsohn提出的一种在图论中使用的算法,能够按照一定规律将二分图中的顶点分割为不同的集合,并且不同集合之间满足一定的偏序关系。
发明内容
针对现有技术中所存在的上述技术问题,本发明提出了一种数字岩心图像的孔喉连通区域识别方法,该方法仅通过一次运算,即可得到数字岩心中所有的孔喉连通区域,计算复杂程度低,数据存储量小,运算速度快,能够正确的对数字岩心中的连通孔喉进行识别。
本申请提出了一种数字岩心图像的孔喉连通区域识别方法,包括以下步骤:
S1、根据数字岩心图像构建像素矩阵,并对像素矩阵中待识别的孔喉像素点进行编号,得到编号矩阵;
S2、根据相邻判定规则判断编号矩阵中像素点的相邻关系,并根据所述相邻关系生成相邻矩阵;
S3、对所述相邻矩阵进行DM分解,得到像素编号排列向量和像素编号分类向量;
S4、根据像素编号排列向量和像素编号分类向量识别数字原始岩心图像中的连通区域。
在一个实施方式中,所述步骤S1具体包括:
将原始数字岩心图像中的孔喉表示为1,原始数字岩心图像中的骨架表示为0,构建像素矩阵。
将所述像素矩阵中待识别的孔喉按顺序进行编号,得到编号矩阵。
在一个实施方式中,所述步骤S2具体包括:
采用4邻域、8邻域或其它相邻判定规则判断编号矩阵中像素点的相邻关系,将相邻的两个像素点编号记录成对。4邻域就是四个方向的像元;8邻域是八个方向的邻域。
根据记录成对的像素点编号生成相邻矩阵。
在一个实施方式中,所述相邻矩阵的行为a={a1,a2,…,ak},列为b={b1,b2,…,bm},当ai所代表的像素点和bj所代表的像素点相邻时,相邻矩阵中的(ai,bj)为1,当ai所代表的像素点和bj所代表的像素点不相邻时,相邻矩阵中的(ai,bj)为0;
其中,ai∈a,bj∈b。
在一个实施方式中,所述步骤S3的具体方法为:利用DM分解算法,得到强连通分量和偏序关系,根据偏序关系获取像素编号排列向量和像素编号分类向量。
在一个实施方式中,根据所述分类向量获取所述排列向量中所对应的孔喉连通区域的像素点的首尾编号,从而获取孔喉连通区域的像素点编号。
在一个实施方式中,所述数字岩心图像中第n个孔喉连通区域的像素点编号为:
在一个实施方式中,所述相邻矩阵存储于稀疏矩阵中。
在一个实施方式中,所述相邻矩阵中的元素采用布尔类型,每个元素仅占1Byte。
与现有技术相比,本发明的优点在于,本发明仅通过一次运算,即可得到数字岩心中所有的孔喉连通区域,计算复杂程度低,数据存储量小,运算速度快。能够正确的对数字岩心中的连通孔喉进行识别。
附图说明
下面将结合附图来对本发明的优选实施例进行详细地描述。在图中:
图1显示了根据本发明的实施例所述的数字岩心图像的孔喉连通区域识别方法的流程图。
图2显示了一种原始数字岩心图像。
图3显示了根据本发明的实施例所述的数字岩心图像的孔喉连通区域识别方法对图2所示原始数字岩心图像进行识别得到的相邻矩阵示意图。
图4显示了根据本发明的实施例所述的数字岩心图像的孔喉连通区域识别方法对图2所示原始数字岩心图像进行识别得到的识别效果示意图。
图5显示了根据本发明的实施例所述的数字岩心图像的孔喉连通区域识别方法所得到的孔喉连通区域示意图。
在附图中,相同的部件使用相同的附图标记。附图并未按照实际的比例绘制。
具体实施方式
以下结合附图及实施例来详细说明本发明的实施方式,借此对本发明如何应用技术手段来解决技术问题,并达成技术效果的实现过程能充分理解并据以实施。需要说明的是,只要不构成冲突,本发明中的各个实施例以及各实施例中的各个特征可以互相结合,所形成的技术方案均在本发明的保护范围之内。
如图1所示,本发明公开了一种数字岩心图像的孔喉连通区域识别方法,采用平面图的对偶图相邻矩阵进行DM分解得到连通区域。包括以下步骤:
S1、根据数字岩心图像构建像素矩阵,并对像素矩阵中待识别的孔喉像素点进行编号,得到编号矩阵;
S2、根据相邻判定规则判断编号矩阵中像素点的相邻关系,并根据所述相邻关系生成相邻矩阵;
S3、对所述相邻矩阵进行DM分解,得到像素编号排列向量p和像素编号分类向量r;
S4、根据像素编号排列向量和像素编号分类向量识别数字原始岩心图像中的连通区域。
在一个实施方式中,所述步骤S1具体包括:
如图2所示,提供了一张原始数字岩心图像,该图像中的黑色区域表示孔喉,白色区域表示骨架。该原始数字岩心图像的原图像素为13*12,将该原始数字岩心图像中的孔喉表示为1,原始数字岩心图像中的骨架表示为0,构建像素矩阵。该像素矩阵具体如下所示:
将所述像素矩阵中待识别的孔喉按顺序进行编号,该实施方式中,自上而下沿横向依次编号,得到编号矩阵。该编号矩阵如下所示:
在一个实施方式中,所述步骤S2具体包括:
采用4邻域、8邻域或其它相邻判定规则判断编号矩阵中像素点的相邻关系,将相邻的两个像素点编号记录成对。
在本实施例中,采用4邻域相邻判定规则来判断编号矩阵中像素点的相邻关系。例如,针对编号矩阵中的像素点编号1,与其相邻的分别为右侧的像素点编号9和下方的像素点编号2,故将相邻的像素点编号记录成对,也即(1,2)和(1,9)。又如,针对编号矩阵中的像素点编号70,与其相邻的分别为上方的像素点编号69,下方的像素点编号71,左侧的像素点编号66和右侧的像素点编号74,故将该相邻的像素点编号记录成对,也即(70,69)、(70,71)、(70,66)和(70,74)。该编号矩阵中像素点的相邻关系具体如下所示:
根据上述记录成对的像素点编号生成相邻矩阵。所述相邻矩阵如图3所示,所述相邻矩阵中的每个点代表了上述的一对编号像素点。所述相邻矩阵的行为a={a1,a2,…,ak},列为b={b1,b2,…,bm},当ai所代表的像素点和bj所代表的像素点相邻时,相邻矩阵中的(ai,bj)为1,当ai所代表的像素点和bj所代表的像素点不相邻时,相邻矩阵中的(ai,bj)为0;
其中,ai∈a,bj∈b。
在一个实施方式中,所述步骤S3的具体方法为:利用DM分解算法,得到强连通分量和偏序关系,根据偏序关系获取像素编号排列向量p和像素编号分类向量r,其中,
r=(1 17 61 81)。
在一个实施方式中,所述步骤S4具体包括:
根据所述分类向量r获取所述排列向量p中所对应的孔喉连通区域的像素点的首尾编号,从而获取孔喉连通区域的像素点编号。数字岩心图像中第n个孔喉连通区域的像素点编号为:
在该实施例中,r1=1、r2=17、r3=61、r4=81,根据分类向量r将排列向量p分为三个区间:(1-16)、(17-60)和(61-80),也即该数字岩心图像有3个连通区域,根据公式(1)可以得到所述三个孔喉连通区间的像素点编号分别为:
X1={p1,p2,…,p16} (2)
X2={p17,p18,…,p60} (3)
X3={p61,p62,…,p80} (4)
具体的,也即:(1-16)表示排列向量p中第1-16个像素点,所述(17-60)表示排列向量p中第17-60个像素点,所述(61-80)表示排列向量p中第61-80个像素点。
故该数字岩心图像的连通区域的像素点编号为:
X1=[1 2 8 9 10 11 18 19 20 21 30 31 32 33 42 43];
X2=[3 4 5 6 7 12 13 14 15 16 17 22 23 24 25 26 27 28 29 34 35 36 3738 39 40 41 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60];
X3=[61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80]。
根据所述连通区域的像素点编号可以获取最终的识别效果,如图4所示,所述区域X1、区域X2和区域X3即为对图2采用本发明所述识别方法识别出的孔喉连通区域。所述区域X1对应图2数字岩心图像中的A区域,区域X2对应图2数字岩心图像中的B区域,区域X3对应图2数字岩心图像中的C区域。从图4能够看出,A区域的孔喉相连通,B区域的孔喉相连通;C区域的孔喉相连通。A、B和C区域之间不连通。
此为一种较简单的实施例,图5给出了较为复杂的数字岩心图像采用本发明所述识别方法识别出的孔喉连通区域示意图。图中颜色深度一致的区域表示相连通的孔喉,例如:a1、a2、a3和a4区域颜色深度一致,表明该四个孔喉相连通,同样的,b1、b2、b3和b4区域颜色深度一致,表明该四个孔喉相连通。可以理解的,颜色深度不一致的区域之间不连通。
得到孔喉连通性的基础上,可以通过以下两种方法判断该连通区域是否与岩心流动的入口和出口相连。一是将种子布置在岩心的入口端,并判断该连通区域内是否有种子位于岩心出口端;二是直接求取该连通区域内所有点中是否包含位于入口和出口的点。
在一个实施方式中,所述相邻矩阵存储于稀疏矩阵中。且所述相邻矩阵中的元素采用布尔类型,每个元素仅占1Byte。
与现有技术相比,本发明的优点在于,本发明仅通过一次运算,即可计算得到数字岩心中所有的孔喉连通区域,计算复杂程度低,数据存储量小,运算速度快。能够正确的对数字岩心中的连通孔喉进行识别。
以上所述仅为本发明的优选实施方式,但本发明保护范围并不局限于此,任何本领域的技术人员在本发明公开的技术范围内,可容易地进行改变或变化,而这种改变或变化都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求书的保护范围为准。
Claims (7)
1.一种数字岩心图像的孔喉连通区域识别方法,其特征在于,包括以下步骤:
S1、根据数字岩心图像构建像素矩阵,并对像素矩阵中待识别的孔喉像素点进行编号,得到编号矩阵;
S2、根据相邻判定规则判断编号矩阵中像素点的相邻关系,并根据所述相邻关系生成相邻矩阵;
S3、对所述相邻矩阵进行DM分解,得到像素编号排列向量和像素编号分类向量;
S4、根据像素编号排列向量和像素编号分类向量识别数字原始岩心图像中的连通区域;
其中,利用DM分解算法,得到强连通分量和偏序关系,根据所述偏序关系获取所述像素编号排列向量和所述像素编号分类向量;根据所述像素编号分类向量获取所述像素编号排列向量中所对应的每个孔喉连通区域的像素点的首尾编号,从而获取孔喉连通区域的像素点编号。
2.根据权利要求1所述的数字岩心图像的孔喉连通区域识别方法,其特征在于,所述步骤S1具体包括:
将原始数字岩心图像中的孔喉表示为1,原始数字岩心图像中的骨架表示为0,构建像素矩阵;
将所述像素矩阵中待识别的孔喉按顺序进行编号,得到编号矩阵。
3.根据权利要求1所述的数字岩心图像的孔喉连通区域识别方法,其特征在于,所述步骤S2具体包括:
采用4邻域或8邻域的相邻判定规则判断编号矩阵中像素点的相邻关系,将相邻的像素点编号记录成对,并根据记录成对的像素点编号生成相邻矩阵。
4.根据权利要求1-3任一项所述的数字岩心图像的孔喉连通区域识别方法,其特征在于,所述相邻矩阵的行为a={a1,a2,…,ak},列为b={b1,b2,…,bm},当ai所代表的像素点和bj所代表的像素点相邻时,相邻矩阵中的(ai,bj)为1,当ai所代表的像素点和bj所代表的像素点不相邻时,相邻矩阵中的(ai,bj)为0;
其中,ai∈a,bj∈b。
6.根据权利要求1所述的数字岩心图像的孔喉连通区域识别方法,其特征在于,所述相邻矩阵存储于稀疏矩阵中。
7.根据权利要求1所述的数字岩心图像的孔喉连通区域识别方法,其特征在于,所述相邻矩阵中的元素采用布尔类型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710491368.6A CN109118552B (zh) | 2017-06-23 | 2017-06-23 | 一种数字岩心图像的孔喉连通区域识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710491368.6A CN109118552B (zh) | 2017-06-23 | 2017-06-23 | 一种数字岩心图像的孔喉连通区域识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109118552A CN109118552A (zh) | 2019-01-01 |
| CN109118552B true CN109118552B (zh) | 2022-05-20 |
Family
ID=64733824
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710491368.6A Active CN109118552B (zh) | 2017-06-23 | 2017-06-23 | 一种数字岩心图像的孔喉连通区域识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109118552B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111738978B (zh) * | 2020-03-27 | 2023-07-25 | 中国石油化工股份有限公司 | 储层孔喉连通性的评价方法、装置、电子设备及存储介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103114848A (zh) * | 2013-01-18 | 2013-05-22 | 西南石油大学 | 一种基于岩心测量的地层裂缝空间重构方法 |
| CN104594888A (zh) * | 2014-12-03 | 2015-05-06 | 中国石油天然气股份有限公司 | 一种基于致密储层导电因素实验的油层识别方法 |
| CN106053315A (zh) * | 2016-06-06 | 2016-10-26 | 中国石油大学(华东) | 一种生物碎屑灰岩储层孔隙结构分类方法 |
| CN106199745A (zh) * | 2015-05-05 | 2016-12-07 | 中国石油化工股份有限公司 | 一种连通的孔喉的识别方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016049204A1 (en) * | 2014-09-25 | 2016-03-31 | Ingrain, Inc. | Digital pore alteration methods and systems |
-
2017
- 2017-06-23 CN CN201710491368.6A patent/CN109118552B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103114848A (zh) * | 2013-01-18 | 2013-05-22 | 西南石油大学 | 一种基于岩心测量的地层裂缝空间重构方法 |
| CN104594888A (zh) * | 2014-12-03 | 2015-05-06 | 中国石油天然气股份有限公司 | 一种基于致密储层导电因素实验的油层识别方法 |
| CN106199745A (zh) * | 2015-05-05 | 2016-12-07 | 中国石油化工股份有限公司 | 一种连通的孔喉的识别方法 |
| CN106053315A (zh) * | 2016-06-06 | 2016-10-26 | 中国石油大学(华东) | 一种生物碎屑灰岩储层孔隙结构分类方法 |
Non-Patent Citations (1)
| Title |
|---|
| 基于连通区域矩阵的文本图像识别方法;郭晓宇 等;《信息工程大学学报》;20121029;第13卷(第3期);第1-5页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109118552A (zh) | 2019-01-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Behjati et al. | Single image super-resolution based on directional variance attention network | |
| CN111738318B (zh) | 一种基于图神经网络的超大图像分类方法 | |
| CN109493346B (zh) | 一种基于多损失的胃癌病理切片图像分割方法和装置 | |
| Zhang et al. | Joint fully convolutional and graph convolutional networks for weakly-supervised segmentation of pathology images | |
| Xiao et al. | Example‐Based Colourization Via Dense Encoding Pyramids | |
| CN113379764A (zh) | 一种基于域对抗自监督学习的病理图像分割方法 | |
| JP5974589B2 (ja) | 画像処理装置およびプログラム | |
| CN103617596A (zh) | 一种基于流型转换的图像色彩风格变换方法 | |
| WO2016066042A1 (zh) | 商品图片的分割方法及其装置 | |
| CN104504745A (zh) | 一种基于图像分割和抠图的证件照生成方法 | |
| CN116310666B (zh) | 一种遥感图像自监督训练方法、装置和电子设备 | |
| CN108229576A (zh) | 一种跨倍率病理图像特征学习方法 | |
| CN110084136A (zh) | 基于超像素crf模型的上下文优化室内场景语义标注方法 | |
| CN111753714B (zh) | 基于字符分割的多方向自然场景文本检测方法 | |
| JP6389742B2 (ja) | 画像セグメンテーション方法、装置、及びプログラム | |
| CN117830848B (zh) | 基于类引导注意力机制的高光谱图像分类方法 | |
| JP2005513632A (ja) | 水源方法を使用したイメージの分割 | |
| CN109118552B (zh) | 一种数字岩心图像的孔喉连通区域识别方法 | |
| CN116310899B (zh) | 基于YOLOv5改进的目标检测方法及装置、训练方法 | |
| CN106204621A (zh) | 一种自适应分治的高分辨率图像对象精准提取方法 | |
| CN106023222A (zh) | 一种基于融合的数字抠图方法 | |
| CN116152458B (zh) | 一种基于影像的三维仿真建筑生成方法 | |
| CN111899161B (zh) | 一种超分辨率重建方法 | |
| CN107392967B (zh) | 一种基于多峰高斯分布函数的彩色图像灰度化方法 | |
| CN107464273B (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 |




















