CN101556703B - 基于连续切片图像的网络模型建立方法 - Google Patents
基于连续切片图像的网络模型建立方法 Download PDFInfo
- Publication number
- CN101556703B CN101556703B CN2009100156042A CN200910015604A CN101556703B CN 101556703 B CN101556703 B CN 101556703B CN 2009100156042 A CN2009100156042 A CN 2009100156042A CN 200910015604 A CN200910015604 A CN 200910015604A CN 101556703 B CN101556703 B CN 101556703B
- Authority
- CN
- China
- Prior art keywords
- venturi
- image
- network model
- hole
- pore
- 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
Images
Abstract
本发明公开了一种基于连续切片图像的网络模型建立方法,属于图像处理技术领域,涉及到对地下原油、天然气、地下水等流体所在的多孔介质(岩石)表征的网络模型的建立方法,特征是对微观模型内部一组连续的二维切片图像进行图像处理,进而分割出岩石颗粒和孔隙空间,在对孔隙空间进行三维图像重建的基础上,提取孔隙和喉道大小分布及拓扑信息,最终建立网络模型,实现了图像信息向三维空间中孔隙和喉道大小分布以及拓扑信息的转化。本发明是基于真实多孔介质的切片图像进行网络模型的建立,针对性强、实用性高,可有效地用于真实岩石孔隙空间及其渗流特性的描述和数字化管理。
Description
技术领域
本发明涉及图像处理技术领域,特别是对地下原油、天然气、地下水等流体所在的多孔介质(岩石)表征的网络模型的建立方法。
背景技术
在油田开发过程中,常常需要借助网络模拟手段,从微观角度来研究油水运动规律。网络模拟优势在于能够方便、灵活的构造出岩石各种类型的孔隙结构特征和润湿特征,较好地表征多孔介质的微观静态特征,快速高效地模拟不同静态特征岩石的微观流动参数。同时,随着数字化岩心技术的发展,将逐步实现岩心数据在计算机中的存储,利用微观模拟模型可以模拟岩心中的孔隙/喉道大小分布、孔隙的网络拓扑结构、相对渗透率曲线,毛管压力曲线等。微观渗流实验研究近十多年来发展较快。目前常用的微观研究手段包括:层析成像仪(CT)、核磁共振成像仪(NMRI或MRI)和显微物理模拟。其中,CT的工作原理是利用X射线在被测物上衰减程度不同检测其内部结构。在渗流研究中,它可以直观地反映流体在孔隙介质中的渗流状态,扫描多个切片便能重建三维特征。要实现微观实验和模拟的结合,需要解决的关键技术在于在准确表征三维微观孔喉结构。其核心技术是建立三维多孔介质的数据场以及孔隙和喉道大小分布以及拓扑信息的提取。在建立三维多孔介质的数据场方面,目前,建立三维多孔介质数据场的方法主要有:①根据CT切片等二维图像的统计信息,得到岩石颗粒、孔隙尺寸的分布,进而通过随机模拟的方法建立三维数据场;②直接在三维空间中通过高斯随机模拟的方法建立三维数据场。目前获取孔隙和喉道大小分布以及拓扑信息的方法主要靠人为地给定,这种方法缺点是和实际岩石中的孔隙空间分布差异较大。
发明内容
本发明的目的是要提供一种基于连续切片图像的网络模型建立方法,对地下原油、天然气、地下水等流体所在的多孔介质(岩石)进行表征。
本发明实现其目的所采取的技术方案是在对连续的二维切片图像进行分割图像处理的基础上,进行三维重建,再通过计算机细化算法得到真实岩石孔隙空间分布的孔隙和喉道大小分布以及拓扑信息,最终建立网络模型,具体过程如下:
(1)通过CT扫描手段获取微观模型内部一组连续的二维切片图像。
(2)确定研究的矩形区域,对图像进行亮度调整、对比度调整、锐化操作,进而在图像重建方向上利用距离加权平均方法进行图像插值。
(3)利用指示克里金方法提取出图像中的岩石部分和孔隙空间部分。首先以岩心实测孔隙度参数为标准,与图像中灰度值频率分布曲线相结合确定基础阈值。然后再给定两个阈值T0和T1,它们的取值范围分别为基础阈值的75%~85%和120%~130%。将需要识别的图像按以下规则确定:对于灰度值小于T0的部分就标记为孔隙空间;对于灰度值大于T1的部分就标记岩石颗粒;对于灰度值介于T0和T1之间的部分,就根据指示变量计算该处条件概率的线性概率,分割出相应的孔隙空间和岩石颗粒部分。分割图像完成后,进一步通过计算机三维图像重建得到孔隙空间的三维分布图像。
(4)对上述步骤得到的孔隙空间部分利用细化算法得到节点和路径,并进一步结合形态学、拓扑学原理对节点和路径数据,根据体素的空间分布计算节点的位置,与每个孔隙相连的路径的数目、节点的燃烧数、路径的长度、节点和路径的数目,最终得到喉道以及孔隙的几何参数和拓扑结构。
(5)利用计算机图像处理中的细化算法进行细化操作,根据Poiseuille定律、网络模型中孔隙空间大小不变原理、孔隙之间喉道的唯一性、孔隙和喉道形状因子等价这些原则对提取的孔隙和喉道大小分布以及拓扑结构进行修正,具体步骤包括:
a.网络模型中每个喉道对应到唯一的一组和它相连的孔隙修正:如果两个孔隙之间可能会存在多个喉道相连,那么需要根据Poiseuille定律,考虑形状因子的影响,对喉道的面积、形状、长度进行等效计算。其中,孔隙和喉道形状因子G定义为它们的截面面积A和周长P的函数,即
G=A/P2
b.网络模型中喉道的长度修正:通过细化算法得到的喉道长度包括了一部分孔隙半径的长度,而网络模型中喉道的长度不包含孔隙半径的长度,因而需要根据其连接关系消除孔隙半径的影响。
c.喉道形状的修正:由于网络模型喉道形状只能是三角形、圆形和正方形,因而需要根据形状因子进行近似处理。
d.孔隙形状和半径的修正:孔隙的形状因子主要根据其连接的喉道形状进行面积加权平均计算,然后按照类似修正喉道形状的方法进行近似。进一步根据喉道所占空间大小以及孔隙的形状计算孔隙半径。
(6)依照网络模型规则对孔隙和喉道大小分布以及拓扑信息进行修正后,得到以节点和路径这样的拓扑结构来表示的孔隙空间分布信息。其中,节点表示了孔隙所处的位置,而路径则表示了连接相邻两个孔隙之间的流动通道,路径上面积最小的截面即为喉道截面,与每个孔隙相连的路径数目即为孔隙配位数,节点的燃烧数即为孔隙半径,路径长度即为喉道长度,最终建立起网络模型。
本发明利用真实多孔介质的切片图像,通过计算机图像处理技术提取了孔隙空间骨架、喉道信息以及孔隙信息,建立了相应的网络模型,实现了图像信息向三维空间中孔隙/喉道大小分布以及拓扑信息的转化,克服了传统方法的缺陷。
附图说明
图1为实验用岩心模型扫描切片图(512×512像素)。
图2为实验用岩心模型图像重建区域图(221×221像素)。
图3为CT图像预处理以后的效果图。
图4(a)为原始CT切片A。
图4(b)为原始CT切片B。
图4(c)为插值后的图像效果图。
图5为用指示克里金方法分割以后的效果图。
图6为岩心孔隙部分三维图像。
图7(a)为细化算法步骤中三维图像体素标记以后结果示意图。
图7(b)为孔隙空间部分体素燃烧过程示意图。
图8为细化以后得到的节点和路径分布图。
图9为最终得到的网络模型图。
图10为网络模型提取流程图。
具体实施方式
以下通过对垦东70-1井岩心模型构建三维网络模型的实例并结合附图来详细说明本发明的内容及实现原理:
1.获取连续的二维切片图像
垦东70-1井岩心模型在进行CT切片扫描时,对未饱和流体的空岩心进行扫描。岩心扫描总长度确定为1.25mm,共扫描切片100张,每两张CT切片之间的间隔为0.0125mm。扫描切片如附图1所示。
2.确定模型区域、进行预处理和图像插值
由于CT实验扫描得到的岩心图片中的岩心区域是圆形,而进行网络模拟的模型一般是一个规则的长方体或者立方体,所以有必要在一系列连续的CT图像中取出一组定位相同的矩形区域。实际CT扫描图像(图像大小都为512×512像素)截取出的矩形区域如附图2所示。经过亮度调整、对比度调整、锐化处理操作后的结果如附图3所示。
垦东70-1井岩心模型的CT二维切片图像的分辨率为8.04微米/像素,而每两张切片之间的间隔为12.5微米。已知相邻的两张连续的切片A、B如附图4(a)、(b)所示,为保证三维空间各个方向上分辨率相同,对这两张切片图像中对应的像素点灰度值利用距离加权平均方法进行插值处理,得到的插值图像如4(c)所示。对100张CT切片图像进行插值,插值以后得到153张图像。
3.利用指示克里金方法进行分割以后,得到孔隙空间三维图像
根据岩心实测孔隙度(孔隙度为34.0%)和二维图像中灰度值频率分布曲线确定基础阈值为93,分别选择T0为75、T1为120,得到图像的分割效果如附图5所示。分割图像中孔隙空间部分对应的图像灰度值都为1,对应图中的白色部分;而岩石颗粒部分对应的图像灰度值都为0,对应图中的黑色部分。进一步通过计算机三维图像重建得到孔隙空间的三维分布图像,如附图6所示。
4.孔隙/喉道几何信息以及拓扑信息的提取
利用细化算法进行孔隙/喉道大小分布以及拓扑信息的提取,细化算法的步骤为:
①将研究区域以外的所有体素标记为-1,将所有岩石颗粒对应的体素标记为0(即该层体素的燃烧数为0),而所有孔隙空间对应的体素标记为1,如附图7(a)所示。
②赋值k=1,沿着燃烧方向,将所有紧连着具有燃烧数为k-1体素的孔隙空间对应的体素(尚未标记燃烧数)标记为k,对于每个这样的体素,记录其相邻的燃烧数为k-1的体素的相对位置,如附图7(b)所示。
③赋值k=k+1,重复步骤②,直到燃烧结束。
最终通过细化得到孔隙空间中的节点和路径分布,如附图8所示。其中,节点表示了孔隙所处的位置,而路径则表示了连接相邻两个孔隙之间的流动通道,在路径上存在的面积最小截面即为喉道截面。与每个孔隙相连的路径数目即为孔隙配位数,节点的燃烧数即为孔隙半径,路径长度即为喉道长度。
5.网络模型的建立
依照网络模型规则对孔隙和喉道大小分布以及拓扑信息进行修正,具体包括:
①网络模型中每个喉道对应到唯一的一组和它相连的孔隙修正:如果两个孔隙之间可能会存在多个喉道相连,那么需要根据Poiseuille定律,考虑形状因子的影响,对喉道的面积、形状、长度进行等效计算;
②网络模型中喉道的长度修正:通过细化算法得到的喉道长度包括了一部分孔隙半径的长度,而网络模型中喉道的长度不包含孔隙半径的长度,因而需要根据其连接关系消除孔隙半径的影响;
③喉道形状的修正:由于网络模型喉道形状只能是三角形、圆形和正方形,因而需要根据形状因子进行近似处理;
④孔隙形状和半径的修正:孔隙的形状因子主要根据其连接的喉道形状进行面积加权平均计算,然后按照类似修正喉道形状的方法进行近似。进一步根据喉道所占空间大小以及孔隙的形状计算孔隙半径。
最终得到的网络模型如附图9所示,图中球体代表了孔隙,线段代表了连接孔隙之间的喉道。附图9主要显示出孔隙和喉道的拓扑关系,事实上,每个孔隙和喉道都具有不同的形状和大小。
表1建立的网络模型基本参数
附图10为本发明网络模型提取流程图。与其他方法相比,本发明实现了CT三维图像信息向孔隙/喉道大小分布以及拓扑信息的转化,建立与真实岩石相吻合的网络模型,可有效地用于真实岩石孔隙空间及其渗流特性的描述和数字化管理。
Claims (1)
1.一种基于连续切片图像的网络模型建立方法,其特征在于:对微观模型内部一组连续的二维切片图像进行图像处理,进而分割出岩石颗粒和孔隙空间,在对孔隙空间进行三维图像重建的基础上,提取孔隙和喉道大小分布以及拓扑信息,最终建立网络模型,按如下过程实现:
(1)通过CT扫描手段获取微观模型内部一组连续的二维切片图像;
(2)确定研究的矩形区域,对图像进行亮度调整、对比度调整、锐化操作,进而在图像重建方向上利用距离加权平均方法进行图像插值;
(3)利用指示克里金方法提取出图像中的岩石部分和孔隙空间部分,首先以岩心实测孔隙度参数为标准,与图像中灰度值频率分布曲线相结合确定基础阈值,然后再给定两个阈值T0和T1,它们的取值范围分别为基础阈值的75%~85%和120%~130%;将需要识别的图像按以下规则确定:对于灰度值小于T0的部分就标记为孔隙空间;对于灰度值大于T1的部分就标记岩石颗粒;对于灰度值介于T0和T1之间的部分,就根据指示变量计算该处条件概率的线性概率,分割出相应的孔隙空间和岩石颗粒部分,分割图像完成后,进一步通过计算机三维图像重建得到孔隙空间的三维分布图像;
(4)对上述步骤得到的孔隙空间部分利用细化算法得到节点和路径,并进一步结合形态学、拓扑学原理对节点和路径数据,根据体素的空间分布计算节点的位置,与每个孔隙相连的路径的数目、节点的燃烧数、路径的长度、节点和路径的数目,得到喉道以及孔隙的几何参数和拓扑结构;
(5)利用计算机图像处理中的细化算法进行细化操作,根据Poiseuille定律、网络模型中孔隙空间大小不变原理、孔隙之间喉道的唯一性、孔隙和喉道形状因子等价这些原则对提取的孔隙和喉道大小分布以及拓扑结构进行修正,具体步骤包括:
a.网络模型中每个喉道对应到唯一的一组和它相连的孔隙修正:如果两个孔隙之间存在多个喉道相连,那么需要根据Poiseuille定律,考虑形状因子的影响,对喉道的面积、形状、长度进行等效计算,其中,孔隙和喉道形状因子G定义为它们的截面面积A和周长P的函数,即
G=A/P2
b.网络模型中喉道的长度修正:通过细化算法得到的喉道长度包括了一部分孔隙半径的长度,而网络模型中喉道的长度不包含孔隙半径的长度,因而需要根据其连接关系消除孔隙半径的影响;
c.喉道形状的修正:由于网络模型喉道形状只能是三角形、圆形和正方形,因而需要根据形状因子进行近似处理;
d.孔隙形状和半径的修正:孔隙的形状因子根据其连接的喉道形状进行面积加权平均计算,然后按照类似修正喉道形状的方法进行近似,进一步根据喉道所占空间大小以及孔隙的形状计算孔隙半径;
(6)依照网络模型规则对孔隙和喉道大小分布以及拓扑信息进行修正后,得到以节点和路径这样的拓扑结构来表示的孔隙空间分布信息,其中,节点表示了孔隙所处的位置,而路径则表示了连接相邻两个孔隙之间的流动通道,路径上面积最小的截面即为喉道截面,与每个孔隙相连的路径数目即为孔隙配位数,节点的燃烧数即为孔隙半径,路径长度即为喉道长度,最终建立起网络模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100156042A CN101556703B (zh) | 2009-05-16 | 2009-05-16 | 基于连续切片图像的网络模型建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100156042A CN101556703B (zh) | 2009-05-16 | 2009-05-16 | 基于连续切片图像的网络模型建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101556703A CN101556703A (zh) | 2009-10-14 |
CN101556703B true CN101556703B (zh) | 2011-01-26 |
Family
ID=41174808
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100156042A Expired - Fee Related CN101556703B (zh) | 2009-05-16 | 2009-05-16 | 基于连续切片图像的网络模型建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101556703B (zh) |
Families Citing this family (50)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101706966B (zh) * | 2009-11-06 | 2011-11-16 | 上海第二工业大学 | 基于二维图像和多点统计方法的多孔介质三维重构方法 |
CN102129548B (zh) * | 2011-03-15 | 2013-01-09 | 中国石油大学(华东) | 多孔介质中流动泡沫结构图像实时采集与识别方法 |
CN102222359B (zh) * | 2011-05-24 | 2012-09-26 | 中国石油天然气股份有限公司 | 一种岩心三维孔隙结构重构方法 |
US9183326B2 (en) * | 2011-07-12 | 2015-11-10 | Ingrain, Inc. | Method for simulating fractional multi-phase/multi-component flow through porous media |
CN102628354B (zh) * | 2012-04-12 | 2014-08-13 | 东北石油大学 | 孔隙微米级油水分布识别量化方法 |
US9080946B2 (en) * | 2012-06-15 | 2015-07-14 | Ingrain, Inc. | Digital rock analysis systems and methods with multiphase flow REV determination |
US10514372B2 (en) | 2012-08-23 | 2019-12-24 | Halliburton Energy Services, Inc. | Digital rock analysis systems and methods that reliably predict a porosity-permeability trend |
CN104076046B (zh) * | 2013-03-28 | 2016-08-31 | 中国石油化工股份有限公司 | 多孔介质中剩余油微观分布图像采集与定量表征方法 |
CN103325118A (zh) * | 2013-06-26 | 2013-09-25 | 中国石油大学(北京) | 一种获取碳酸盐岩岩心孔洞特征参数的方法及装置 |
US9311720B2 (en) * | 2013-11-26 | 2016-04-12 | Saudi Arabian Oil Company | Automated saw cut correction for 3D core digital modeling from computerized tomography scanner (CTS) images |
CN104778678B (zh) * | 2014-10-09 | 2016-04-06 | 中国石油大学(华东) | 一种考虑孔喉末端的孔隙喉道识别方法 |
CN105654525B (zh) * | 2014-11-10 | 2018-06-29 | 中国石油化工股份有限公司 | 反映真实岩心连通性的2d模型构建方法 |
CN104535475B (zh) * | 2015-01-08 | 2018-02-06 | 中国石油天然气股份有限公司 | 碳酸盐岩微观结构的确定方法及装置 |
CN104794709B (zh) * | 2015-04-10 | 2018-03-09 | 四川大学 | 一种三维岩心图像孔隙和喉道的分割方法 |
CN104751473B (zh) * | 2015-04-10 | 2017-10-27 | 中国石油天然气股份有限公司 | 碳酸岩多尺度孔洞特征的确定方法及装置 |
CN104933760B (zh) * | 2015-06-18 | 2019-01-08 | 中国地质大学(北京) | 一种重构土壤ct图片三维重建及土壤孔隙搜索方法 |
CN106368688A (zh) * | 2015-07-21 | 2017-02-01 | 中国石油化工股份有限公司 | 研究岩石受热破裂条件和裂缝分布的实验方法 |
CN106716114B (zh) * | 2015-09-11 | 2020-06-09 | 数岩科技(厦门)股份有限公司 | 多孔介质分析系统和方法 |
CN105427383B (zh) * | 2015-11-23 | 2017-04-05 | 中国石油大学(华东) | 一种考虑凹凸性的岩石孔隙网络模型的孔喉截面构造方法 |
CN106780744B (zh) * | 2016-12-27 | 2020-04-10 | 中国石油天然气集团公司 | 利用不同分辨率ct图像构建多尺度三维数字岩心的方法 |
CN106918484A (zh) * | 2017-03-28 | 2017-07-04 | 武汉瑞福宁科技有限公司 | 一种基于组织切片的三维模型的构建方法 |
CN108303360B (zh) * | 2017-07-31 | 2020-12-04 | 中国矿业大学 | 一种煤岩三维孔隙网络结构参数表征方法 |
CN111615625A (zh) * | 2017-11-06 | 2020-09-01 | 哈里发科学技术大学 | 用于确定多孔介质的渗透率的方法和系统 |
CN107939387B (zh) * | 2017-11-14 | 2021-04-09 | 中国石油大学(北京) | 一种微观岩石网络模型制作方法 |
CN110320137A (zh) * | 2018-03-29 | 2019-10-11 | 中国石油化工股份有限公司 | 一种基于数字岩心的多尺度融合方法 |
CN108763650B (zh) * | 2018-04-28 | 2022-03-11 | 湘潭大学 | 一种覆岩采动裂隙网络模型构建方法 |
CN108961388B (zh) * | 2018-06-05 | 2023-03-21 | 哈尔滨工业大学深圳研究生院 | 微观结构三维建模方法、装置、设备和计算机存储介质 |
CN108734763B (zh) * | 2018-06-15 | 2022-07-05 | 重庆大学 | 微装配系统的微装配空间的数字化方法及系统 |
CN108491677A (zh) * | 2018-07-04 | 2018-09-04 | 河海大学 | 基于改进最大球法的微观孔隙模型的孔隙特征统计方法 |
CN109063348B (zh) * | 2018-08-09 | 2022-01-04 | 中国石油天然气股份有限公司 | 基于孔隙喉道网络模型的驱替模拟方法及装置 |
CN109887073B (zh) * | 2018-08-16 | 2022-12-09 | 清能艾科(深圳)能源技术有限公司 | 岩心的三维数字模型构建方法及装置 |
CN109187305A (zh) * | 2018-08-24 | 2019-01-11 | 中国地质大学(武汉) | 一种页岩结构表征方法和装置 |
CN109001236A (zh) * | 2018-09-03 | 2018-12-14 | 东南大学 | 一种基于纳米ct技术表征水泥浆体中孔隙率分布的方法 |
CN109360220A (zh) * | 2018-10-25 | 2019-02-19 | 南昌工程学院 | 一种岩心断层ct序列图像的批处理系统 |
CN109596618B (zh) * | 2018-12-26 | 2021-02-26 | 太原理工大学 | 一种基于断面轮廓序列的微观多相结构三维建模测量方法 |
CN109993786B (zh) * | 2019-03-08 | 2021-05-18 | 中国石油大学(北京) | 迂曲度获取方法、装置、设备以及存储介质 |
CN111829887B (zh) * | 2019-04-22 | 2023-07-04 | 中国石油化工股份有限公司 | 一种基于高压压汞的岩石压裂模拟实验方法 |
CN110263362B (zh) * | 2019-04-25 | 2022-11-29 | 南京大学(苏州)高新技术研究院 | 基于孔隙密度流的岩土体离散元流固耦合数值模拟方法 |
CN110701636B (zh) * | 2019-05-20 | 2020-10-27 | 岳阳县湖顶智能科技有限公司 | 热力发电热源管理方法 |
CN110349159B (zh) * | 2019-06-21 | 2022-10-25 | 浙江大学宁波理工学院 | 基于权重能量自适应分布的三维形状分割方法及系统 |
CN110632653A (zh) * | 2019-08-16 | 2019-12-31 | 中国石油天然气股份有限公司 | 盐构造分析方法、装置及系统 |
CN110441220A (zh) | 2019-09-02 | 2019-11-12 | 西南石油大学 | 微生物碳酸盐岩的孔隙表征方法、装置及服务器 |
CN113155693B (zh) * | 2020-01-07 | 2024-04-02 | 中国石油化工股份有限公司 | 一种孔喉连接关系判别方法、系统、电子设备及存储介质 |
CN111784839B (zh) * | 2020-04-17 | 2023-02-07 | 中国科学院力学研究所 | 一种rev连通孔隙空间的构建方法、装置 |
CN113176187A (zh) * | 2021-03-23 | 2021-07-27 | 中国地质大学(武汉) | 一种表征岩石孔隙网络连通性的方法 |
CN113177949B (zh) * | 2021-04-16 | 2023-09-01 | 中南大学 | 一种大尺寸岩石颗粒特征识别方法及装置 |
CN114092489B (zh) * | 2021-11-02 | 2023-08-29 | 清华大学 | 多孔介质渗流通道提取、模型训练方法、装置及设备 |
CN114565658A (zh) * | 2022-01-14 | 2022-05-31 | 武汉理工大学 | 基于ct技术的孔隙尺寸计算方法及装置 |
CN116977999B (zh) * | 2023-08-03 | 2024-03-08 | 北京建工环境修复股份有限公司 | 一种基于机器视觉的岩芯智能识别方法、系统及存储介质 |
CN117611485B (zh) * | 2024-01-24 | 2024-04-02 | 西南石油大学 | 一种基于时空图神经网络的三维岩心渗透率预测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5913052A (en) * | 1997-01-24 | 1999-06-15 | Lucent Technologies Inc. | System and method for debugging digital signal processor software with an architectural view and general purpose computer employing the same |
CN1545043A (zh) * | 2003-11-10 | 2004-11-10 | 魏 茹 | 可视化电子图书馆 |
CN1746896A (zh) * | 2005-10-14 | 2006-03-15 | 清华大学 | 基于vrml模型的有约束协同装配工艺规划的仿真方法 |
-
2009
- 2009-05-16 CN CN2009100156042A patent/CN101556703B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5913052A (en) * | 1997-01-24 | 1999-06-15 | Lucent Technologies Inc. | System and method for debugging digital signal processor software with an architectural view and general purpose computer employing the same |
CN1545043A (zh) * | 2003-11-10 | 2004-11-10 | 魏 茹 | 可视化电子图书馆 |
CN1746896A (zh) * | 2005-10-14 | 2006-03-15 | 清华大学 | 基于vrml模型的有约束协同装配工艺规划的仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101556703A (zh) | 2009-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101556703B (zh) | 基于连续切片图像的网络模型建立方法 | |
CN105115874B (zh) | 多源信息融合的多组分三维数字岩心构建方法 | |
CN106780744B (zh) | 利用不同分辨率ct图像构建多尺度三维数字岩心的方法 | |
CN105193445B (zh) | 一种超声探头标定体模、超声探头标定系统及其标定方法 | |
CN105279790B (zh) | 裂缝网络三维数字岩心建模方法 | |
KR101392804B1 (ko) | 지상 라이다에서 취득된 포인트 클라우드 기반의 실내공간 3차원 모델 추출 방법 및 이를 구현하기 위한 프로그램이 기록된 기록매체 | |
CN105261068A (zh) | 基于Micro-CT技术的储层岩心三维实体模型重构方法 | |
CN105487121A (zh) | 基于ct扫描图像与电成像图像融合构建多尺度数字岩心方法 | |
CN110135311A (zh) | 一种基于三维岩心扫描图像的孔隙与孔喉识别系统及方法 | |
CN110320137A (zh) | 一种基于数字岩心的多尺度融合方法 | |
Yilmaz et al. | Comparison of data reduction algorithms for Li DAR‐derived digital terrain model generalisation | |
CN102737542B (zh) | 一种顾及多重约束条件的水深注记自动选取方法 | |
CN105447274A (zh) | 一种利用面向对象分类技术对中等分辨率遥感图像进行滨海湿地制图的方法 | |
CN109345625A (zh) | 岩心图像自适应分区三维重建方法 | |
CN104809756A (zh) | 基于X-ray CT图像的沥青混合料空隙空间结构重构方法 | |
CN113029899B (zh) | 一种基于显微图像处理的砂岩渗透率计算方法 | |
CN101477533B (zh) | 地理格网下不易获取地理要素空间渐变信息的数字制图方法 | |
CN109064556B (zh) | 一种面向isr的地貌高精度仿真建模系统 | |
CN112147680A (zh) | 岩石内部微裂缝提取和刻画方法及系统 | |
CN108038903B (zh) | 用于构建岩心模型的三维数字模型生成方法 | |
CN116756836B (zh) | 一种隧道超欠挖体积计算方法、电子设备及存储介质 | |
CN105654525B (zh) | 反映真实岩心连通性的2d模型构建方法 | |
Geisthövel et al. | Automated Swiss-style relief shading and rock hachuring | |
Romero et al. | Analytical applications of fine-scale terrestrial lidar at the imperial Inca site of Caranqui, northern highland Ecuador | |
CN110222368A (zh) | 一种利用二维切片计算岩心三维孔隙度和渗透率的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110126 Termination date: 20110516 |