CN101738402A - 一种岩样二维ct图像分析系统 - Google Patents
一种岩样二维ct图像分析系统 Download PDFInfo
- Publication number
- CN101738402A CN101738402A CN200810226277A CN200810226277A CN101738402A CN 101738402 A CN101738402 A CN 101738402A CN 200810226277 A CN200810226277 A CN 200810226277A CN 200810226277 A CN200810226277 A CN 200810226277A CN 101738402 A CN101738402 A CN 101738402A
- Authority
- CN
- China
- Prior art keywords
- image
- rock sample
- interest
- gray
- region
- 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 49
- 238000010191 image analysis Methods 0.000 title abstract description 3
- 238000004451 qualitative analysis Methods 0.000 claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims description 26
- 230000006870 function Effects 0.000 claims description 13
- 208000013201 Stress fracture Diseases 0.000 claims description 11
- 238000000034 method Methods 0.000 claims description 7
- 238000011161 development Methods 0.000 claims description 6
- 239000003086 colorant Substances 0.000 claims description 4
- 230000005311 nuclear magnetism Effects 0.000 claims description 4
- 238000012113 quantitative test Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000004744 fabric Substances 0.000 claims description 3
- 235000019994 cava Nutrition 0.000 claims description 2
- 238000004445 quantitative analysis Methods 0.000 abstract description 7
- 238000004141 dimensional analysis Methods 0.000 abstract 2
- 239000011159 matrix material Substances 0.000 abstract 1
- 239000000725 suspension Substances 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 239000007789 gas Substances 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 239000003345 natural gas Substances 0.000 description 2
- 239000003209 petroleum derivative Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000010244 region-of-interest analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及一种岩样二维CT图像分析系统;图像读取模块打开从岩样核磁CT仪器中采集的tiff格式图片,从内存中取出对应的图片数据显示在主窗口中;显示控制模块对已选取的图像进行灰度输出控制,提供直观的图像对岩芯进行定性分析;一维分析模块在岩心图像上显示水平和垂直两条线,显示鼠标悬停点的灰度值,直方图的形式显示在主窗口右侧,实现对岩样特征的定量分析;二维分析模块首先创建一个感兴趣区,兴趣区是一个与岩心图像等大的布尔类型矩阵,设为1的点被认为是兴趣区;对兴趣区内的点进行分析,对面孔率、分选系数等参数进行计算,从而得到任意形状的裂缝、微裂缝、孔洞以和非均质性特点的定量结果;实现对岩样定量、定性分析。
Description
技术领域
本发明涉及一种油气储层岩样二维CT图像的定性和定量分析系统。
现有技术
针对石油天然气储层岩石样品中不同程度地发育有裂缝、微裂缝、孔洞和具有非均质性特点,对岩样进行定性和定量的分析能够科学的对岩样做出评价,有利于油气田开发。
目前国内没有专用于岩样二维CT图像分析的软件,只能通过直接观察岩样二维CT图像的方式对岩样作定性分析,不能进行定量分析。医学上应用的CT分析软件不针对岩石样品,不能够对岩石的非均质性进行分析,因此使用功能受到很大局限,不适用于岩石样品分析。
发明内容
本发明的目的是根据石油天然气储层岩石样品中不同程度地发育有裂缝、微裂缝、孔洞和具有非均质性特点,提供一种岩样二维CT图像分析的系统。通过这套系统,提高定性分析的准确性,实现岩样二维CT图像定量分析,达到对岩样内裂缝、微裂缝、孔洞的发育特征、发育程度和岩石非均质性特征进行定量分析的目的。
本发明的岩样二维CT图像分析系统是:通过计算机技术,对从岩样核磁CT仪器中采集的tiff格式图片进行定性和定量分析。本岩样二维CT图像分析系统分为图像读取、显示控制、一维分析和二维分析四大部分:
(1)图像读取模块允许同时打开多幅从岩样核磁CT仪器中采集的tiff格式图片,对其进行解码之后装入内存,同时计算每幅图的缩略图,并显示在缩略图列表中;点选缩略图时,根据鼠标点击位置计算出选择的图像,并从内存中取出对应的图片数据显示在主窗口中;
(2)显示控制模块主要有两个函数,其中伪彩色处理函数能够将灰度图转化为彩色图;灰度显示上下限控制函数能够使上限以上的像素全部显示为白色,下限以下的像素显示为黑色,上下限之间灰度从0-255均匀分布;显示模块的目的是增强视觉显示效果,提供更直观的图像以突出岩石样品中不同程度发育的裂缝、微裂缝、孔洞及非均质性等特点,实现对岩心进行定性分析;
(3)一维分析模块在岩心图像上显示水平和垂直两条直线,并使用鼠标移动它们的交点,同时显示鼠标悬停点的灰度值;通过提取两条线上的图像数据,将它们以直方图的形式显示在主窗口右侧,或者输出到文本中;一维分析模块用于实现对岩样内裂缝、微裂缝、孔洞等的发育特征、发育程度以及岩石非均质性特征等进行定量分析;
(4)二维分析模块运行时,首先需要创建一个感兴趣区,这一过程分两步:第一步是画定区域,有三个形状绘制模块,用于实现使用鼠标在岩心图像上绘制圆形、矩形和任意形状的功能;其中圆形和矩形的尺寸位置可以随时修改,任意形状的边界计算采用了一种逐行扫描以获得选区的算法;另外还有一个模块提供直接输入灰度范围的功能,在选择区域之后,第二步是根据用户的决定将区域添加到兴趣区或从兴趣区内删除,兴趣区是一个与岩心图像等大的布尔类型矩阵,设为1的点被认为是兴趣区,兴趣区创建以后,二维分析模块将对兴趣区内的点进行分析,首先是统计点数并绘出直方图,之后对面孔率、分选系数参数进行计算,从而达到定量研究任意形状的裂缝、微裂缝、孔洞以及非均质性特点的目的。
除了用于打开图像的图像读取模块是每次运行必操作的模块之外,其他三个模块都是各自独立的。它们负责完成不同的任务,共同为岩样二维CT图像分析服务。
发明的效果
该系统在多个油气田储层岩石裂缝、微裂缝、孔洞等的发育特征分析中得到了实际应用,如应用于大庆深层天然气特殊岩性火山岩、砾岩储层岩石孔隙结构特征的分析,取得了很好的应用效果。
附图说明
图1岩样二维CT图像分析系统总流程图。
图2岩样二维CT图像分析系统图像分析处理模块流程图。
图3岩样二维CT图像分析系统兴趣区分析模块流程图。
具体实施实方式
下面根据附图对本发明提供的一种岩样二维CT图像分析系统结构加以叙述。
系统运行并初始化之后,如果用户点击打开图片按钮,那么会弹出一个打开图片对话框,帮助用户选择要打开的图片;如果用户单击新建按钮,则清空已经打开的图片;如果用户单击缩略图,则将该图片设置为活动图,并显示在主窗口中;如果用户单击退出按钮则结束程序。在已设置活动图的基础上,系统允许用户对其进行分析和处理,如图1所示。
点击打开/关闭灰度条按钮会显示或者隐藏灰度条对话框,该对话框用于控制灰度显示上下限,控制方程为:
式中gray为像素点的灰度,top和down分别为灰度条上设置的最大和最小灰度。
点击黑白/彩色按钮能够切换图像色彩(灰度/伪彩色)显示方式,伪彩色算法如下:
式中函数r(),g(),b()分别用于求该像素点红、绿、蓝三个通道的色值,gray为该像素点的灰度值。
用户点击投影线按钮会切换显示/隐藏水平和垂直两条投影线。当用户在主窗口中单击鼠标左键拖动的时候,会改变这两条线交点的位置,并重新计算两条线所经过点的灰度值,这些灰度值被以直方图的形式显示在主窗口右侧。单击保存按钮会打开一个保存文件对话框,用于保存这些数值。
单击兴趣区按钮会打开/关闭感兴趣区(二维分析)模块。
以上这部分功能的流程见图2。
用户在兴趣区模块中创建兴趣区之后,自动统计不同的灰度值在兴趣区中所包含的点数,依此做直方图。并找出最大最小灰度值,作为直方图坐标值的范围。
面孔率、分选系数、均质系数和平均值的计算方法如下:
(1)面孔率(%)=面孔频率/总频率;
(3)均质系数a:a=(∑riai)/rmax;
(4)平均值Gaver为几何平均值。
这套方法的软件在Windows XP以上版本的操作系统下运行,并在中石油廊坊分院渗流所承担的多项研究工作中得到具体的实施。例如在大庆深层天然气特殊岩性火山岩、砾岩储层岩石孔隙结构特征的分析工作中,应用这套方法分析了170个岩心的二维CT岩石图像,给出了每个岩心的图像灰度频率分布图及图像特征描述,不仅对气孔溶洞、裂缝、微裂缝、高密度条带等岩石特征进行了定性描述,还给出了面孔率、图像灰度平均值、分选系数、均质系数等定量分析结果,取得了很好的应用效果。
Claims (1)
1.一种岩样二维CT图像分析系统,包括核磁CT仪器和计算机系统,其特征在于:计算机系统还包括有图像读取、显示控制、一维分析和二维分析四大部分:
(1)图像读取模块允许同时打开多幅从岩样核磁CT仪器中采集的tiff格式图片,对其进行解码之后装入内存,同时计算每幅图的缩略图,并显示在缩略图列表中,点选缩略图时会根据鼠标点击位置计算出选择的是哪一幅图,之后从内存中取出对应的图片数据显示在主窗口中;
(2)显示控制模块是对已选取的图像进行灰度输出控制,包括有两个函数,其中伪彩色处理函数能够将灰度图转化为彩色图,同时调整彩色输出范围,突出岩石样品中不同程度发育的裂缝、微裂缝、孔洞和非均质性特点;灰度显示上下限控制函数能够使上限以上的像素全部显示为白色,下限以下的像素显示为黑色,上下限之间灰度丛0-255均匀分布;提供直观的图像对岩芯进行定性分析;
(3)一维分析模块在岩心图像上显示水平和垂直两条线,显示鼠标悬停点的灰度值,鼠标拖拽时显示水平和垂直投影线灰度分布直方图,查看和保存直方图数据,通过提取两条线上的图像数据,将它们以直方图的形式显示在主窗口右侧,或者输出到文本中实现对岩样内裂缝、微裂缝、孔洞的发育特征、发育程度和岩石非均质性特征进行定量分析;
(4)二维分析模块运行时,首先创建一个感兴趣区,这一过程分两步:第一步是画定区域;包括三个形状绘制模块,提供使用鼠标在岩心图像上绘制圆形、矩形和任意形状的功能,其中圆形和矩形的尺寸位置可以修改,任意形状的边界计算采用一种逐行扫描以获得选区的算法;另外还包括一个提供直接选择灰度范围功能的模块;在选择区域之后,第二步是根据用户的决定将区域添加到兴趣区或从兴趣区内删除;兴趣区是一个与岩心图像等大的布尔类型矩阵,设为1的点被认为是兴趣区;兴趣区创建之后,二维分析模块将对兴趣区内的点进行分析,首先是统计点数并绘出直方图,之后对面孔率、分选系数参数进行计算,从而得到任意形状的裂缝、微裂缝、孔洞以和非均质性特点的定量结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102262770A CN101738402B (zh) | 2008-11-12 | 2008-11-12 | 一种岩样二维ct图像分析系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102262770A CN101738402B (zh) | 2008-11-12 | 2008-11-12 | 一种岩样二维ct图像分析系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101738402A true CN101738402A (zh) | 2010-06-16 |
CN101738402B CN101738402B (zh) | 2011-08-03 |
Family
ID=42462150
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008102262770A Active CN101738402B (zh) | 2008-11-12 | 2008-11-12 | 一种岩样二维ct图像分析系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101738402B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183450A (zh) * | 2011-04-20 | 2011-09-14 | 东北石油大学 | 储层岩心微观孔隙结构原子力显微镜的表征方法 |
CN102230875A (zh) * | 2011-06-17 | 2011-11-02 | 东南大学 | 基于工业x射线断层照相的水泥基材料孔隙率分布测试方法及装置 |
CN103091342A (zh) * | 2011-10-31 | 2013-05-08 | 中国石油化工股份有限公司 | 一种对岩芯样品进行ct扫描分析处理的方法 |
CN103325118A (zh) * | 2013-06-26 | 2013-09-25 | 中国石油大学(北京) | 一种获取碳酸盐岩岩心孔洞特征参数的方法及装置 |
CN103573251A (zh) * | 2012-08-06 | 2014-02-12 | 中国石油化工股份有限公司 | 大尺寸火山岩水力裂缝起裂与扩展ct扫描裂缝监测方法 |
CN103698803A (zh) * | 2012-09-27 | 2014-04-02 | 中国石油天然气股份有限公司 | 一种岩石孔隙结构表征方法及装置 |
CN104299271A (zh) * | 2013-07-18 | 2015-01-21 | 达索系统公司 | 确定三维建模对象组件的分解图的分解路径的计算机实现方法 |
CN106204515A (zh) * | 2015-05-05 | 2016-12-07 | 中国石油化工股份有限公司 | 一种岩心扫描图像的分析及展示系统 |
CN106529022A (zh) * | 2016-11-09 | 2017-03-22 | 中国地质大学(北京) | 一种储层分选性精确划分方法 |
CN107064187A (zh) * | 2017-04-21 | 2017-08-18 | 中国石油天然气股份有限公司 | 一种烃源岩有机碳含量测定方法及装置 |
CN107449707A (zh) * | 2017-07-03 | 2017-12-08 | 中国石油天然气股份有限公司 | 页岩储层中不同尺度孔隙定量的三维表征确定方法和装置 |
CN110572626A (zh) * | 2019-09-03 | 2019-12-13 | 云南白药集团健康产品有限公司 | 一种图像处理方法及设备 |
CN116930244A (zh) * | 2023-09-19 | 2023-10-24 | 东北石油大学三亚海洋油气研究院 | 一种稠油裂缝的核磁图版绘制方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102866080A (zh) * | 2012-09-10 | 2013-01-09 | 中国十七冶集团有限公司 | 基于图像处理分析的护堤抛石密度测定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101086533B (zh) * | 2007-07-06 | 2010-04-14 | 福州华虹智能科技开发有限公司 | 基于嵌入式Linux的浅层地震勘探仪器软件体系结构的装置 |
-
2008
- 2008-11-12 CN CN2008102262770A patent/CN101738402B/zh active Active
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183450B (zh) * | 2011-04-20 | 2012-07-25 | 东北石油大学 | 储层岩心微观孔隙结构原子力显微镜的表征方法 |
CN102183450A (zh) * | 2011-04-20 | 2011-09-14 | 东北石油大学 | 储层岩心微观孔隙结构原子力显微镜的表征方法 |
CN102230875A (zh) * | 2011-06-17 | 2011-11-02 | 东南大学 | 基于工业x射线断层照相的水泥基材料孔隙率分布测试方法及装置 |
CN102230875B (zh) * | 2011-06-17 | 2013-06-05 | 东南大学 | 基于工业x射线断层照相的水泥基材料孔隙率分布测试方法及装置 |
CN103091342A (zh) * | 2011-10-31 | 2013-05-08 | 中国石油化工股份有限公司 | 一种对岩芯样品进行ct扫描分析处理的方法 |
CN103091342B (zh) * | 2011-10-31 | 2015-02-11 | 中国石油化工股份有限公司 | 一种对岩芯样品进行ct扫描分析处理的方法 |
CN103573251B (zh) * | 2012-08-06 | 2016-04-27 | 中国石油化工股份有限公司 | 大尺寸火山岩水力裂缝起裂与扩展ct扫描裂缝监测方法 |
CN103573251A (zh) * | 2012-08-06 | 2014-02-12 | 中国石油化工股份有限公司 | 大尺寸火山岩水力裂缝起裂与扩展ct扫描裂缝监测方法 |
CN103698803A (zh) * | 2012-09-27 | 2014-04-02 | 中国石油天然气股份有限公司 | 一种岩石孔隙结构表征方法及装置 |
CN103698803B (zh) * | 2012-09-27 | 2017-02-08 | 中国石油天然气股份有限公司 | 一种岩石孔隙结构表征方法及装置 |
CN103325118A (zh) * | 2013-06-26 | 2013-09-25 | 中国石油大学(北京) | 一种获取碳酸盐岩岩心孔洞特征参数的方法及装置 |
CN104299271A (zh) * | 2013-07-18 | 2015-01-21 | 达索系统公司 | 确定三维建模对象组件的分解图的分解路径的计算机实现方法 |
CN106204515A (zh) * | 2015-05-05 | 2016-12-07 | 中国石油化工股份有限公司 | 一种岩心扫描图像的分析及展示系统 |
CN106529022A (zh) * | 2016-11-09 | 2017-03-22 | 中国地质大学(北京) | 一种储层分选性精确划分方法 |
CN106529022B (zh) * | 2016-11-09 | 2019-09-06 | 中国石油大学(华东) | 一种储层分选性精确划分方法 |
CN107064187A (zh) * | 2017-04-21 | 2017-08-18 | 中国石油天然气股份有限公司 | 一种烃源岩有机碳含量测定方法及装置 |
CN107449707A (zh) * | 2017-07-03 | 2017-12-08 | 中国石油天然气股份有限公司 | 页岩储层中不同尺度孔隙定量的三维表征确定方法和装置 |
CN107449707B (zh) * | 2017-07-03 | 2020-01-07 | 中国石油天然气股份有限公司 | 页岩储层中不同尺度孔隙定量的三维表征确定方法和装置 |
CN110572626A (zh) * | 2019-09-03 | 2019-12-13 | 云南白药集团健康产品有限公司 | 一种图像处理方法及设备 |
CN110572626B (zh) * | 2019-09-03 | 2021-05-28 | 云南白药集团健康产品有限公司 | 一种图像处理方法及设备 |
CN116930244A (zh) * | 2023-09-19 | 2023-10-24 | 东北石油大学三亚海洋油气研究院 | 一种稠油裂缝的核磁图版绘制方法 |
CN116930244B (zh) * | 2023-09-19 | 2023-12-01 | 东北石油大学三亚海洋油气研究院 | 一种稠油裂缝的核磁图版绘制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101738402B (zh) | 2011-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101738402B (zh) | 一种岩样二维ct图像分析系统 | |
Seethepalli et al. | RhizoVision Explorer: Open-source software for root image analysis and measurement standardization | |
Kokalj et al. | Why not a single image? Combining visualizations to facilitate fieldwork and on-screen mapping | |
Rocchini et al. | Satellite remote sensing to monitor species diversity: Potential and pitfalls | |
Zhang et al. | Adobe photoshop quantification (PSQ) rather than point-counting: A rapid and precise method for quantifying rock textural data and porosities | |
JP5624881B2 (ja) | 位置特有のマルチパラメータデータセットを分析し比較するグラフィックユーザインターフェイス | |
US7308139B2 (en) | Method, system, and apparatus for color representation of seismic data and associated measurements | |
US10753848B2 (en) | Efficient contours and gating | |
Magnini et al. | Theory and practice for an object-based approach in archaeological remote sensing | |
CN114897781B (zh) | 基于改进R-UNet深度学习的透水混凝土孔隙自动识别方法 | |
CN106484990B (zh) | 一种发动机试验数据三维瀑布图重建、显示及分析方法 | |
CN106777707A (zh) | 一种利用改进的蜘蛛网图进行测井岩性定量识别的方法 | |
CN110889834A (zh) | 一种基于云计算的公路隧道围岩分级方法 | |
Cooper et al. | Expediting and standardizing stone artifact refitting using a computerized suitability model | |
Hájek | Process-based approach to automated classification of forest structures using medium format digital aerial photos and ancillary GIS information | |
CN106204515A (zh) | 一种岩心扫描图像的分析及展示系统 | |
WO2018207885A1 (ja) | 解析装置、地層年代推定装置、解析方法、地層年代推定方法、およびプログラム | |
CN111177607B (zh) | 一种基于b/s架构的岩心图像处理和应用方法 | |
CN115438115A (zh) | 油田综合研究大数据样本标注方法 | |
Pham et al. | Interactive visualization of spatial and temporal patterns of diversity and abundance in ecological data | |
CN1277394A (zh) | 一种显微生物医学图像识别系统及其分析方法 | |
van Wijk et al. | Methods for environment: Productivity trade-off analysis in agricultural systems | |
CN112119425A (zh) | 用于图像控制和颗粒分析仪图像显示的方法和装置 | |
CN118097093B (zh) | 基于病理大模型在数字病理切片数据集上的以图搜图系统 | |
CN117763213B (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 |