CN108830140A - 一种基于电成像测井分形维数的火山岩岩性识别方法 - Google Patents
一种基于电成像测井分形维数的火山岩岩性识别方法 Download PDFInfo
- Publication number
- CN108830140A CN108830140A CN201810398397.2A CN201810398397A CN108830140A CN 108830140 A CN108830140 A CN 108830140A CN 201810398397 A CN201810398397 A CN 201810398397A CN 108830140 A CN108830140 A CN 108830140A
- Authority
- CN
- China
- Prior art keywords
- fractal dimension
- lithology
- imaging logging
- electric imaging
- identification methods
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/30—Noise filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/50—Extraction of image or video features by performing operations within image blocks; by using histograms, e.g. histogram of oriented gradients [HoG]; by summing image-intensity values; Projection analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- General Engineering & Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于电成像测井分形维数的火山岩岩性识别方法,该方法主要包括导入取心井段成像测井数据、测井数据滤波与信号增强、阈值分割及骨架信息提取,确定不同岩性的分维数范围和未取心井段的分维数计算,以及导入未取心井段成像测井数据并计算分维数进行岩性识别五部分。本发明通过对电成像测井的数据处理和分形维数计算等,可实现未取心井段岩性的自动识别,在一定程度上克服了传统的肉眼观察法带来的不确定性,提高岩性识别精度和速率。同时,本发明还能为其它岩类的测井识别带来新思路。
Description
技术领域
本发明属于地球物理测井和地质勘探领域,具体地说是涉及一种基于电成像测井分形维数的火山岩岩性识别方法。
背景技术
火山岩的岩性识别是火山岩储层测井评价和油气资源探测的重要基础。然而,受到多种因素的影响,火山岩的非均质性强,导致它们的测井响应特征十分复杂,识别难度极大。
前人基于常规测井资料的交会图分析和常规测井资料的分形维数计算等对火山岩的岩性识别进行了许多有益探索,并取得了部分成果。受制于常规测井的纵向分辨率,它们的识别精度低,很难真正把地下复杂的岩性变化体现出来。随着测井技术的进步,电成像测井通过多组极板和电极,将电流传送到地层,进而得到井周范围内的电流变化情况,将其用图像处理技术显示出来,能展示地层的岩性和储集空间等信息,具有分辨率高,直观成像等优点。
然而,现有的方法和思路都是通过简单的肉眼识别方法,建立不同岩性的典型电成像测井信息库,采用定性的方法来对岩性进行判断,精度低,主观性强,且十分耗时。
发明内容
基于上述技术问题,本发明提供一种基于电成像测井分形维数的火山岩岩性识别方法。
本发明所采用的技术解决方案是:
一种基于电成像测井分形维数的火山岩岩性识别方法,包括以下步骤:
a导入取心井段的成像测井数据,对成像测井数据进行预处理,消除噪声;
a1确定取心井段火山岩的岩性类型,应用Matlab软件,分别读入它们所对应深度段的电成像测井数据;
a2先采用中值滤波算法消除整体图像中的大部分噪声,再采用形态学滤波的开运算进一步将细小的噪声消除,突出岩石骨架和孔隙信息;
a3对去噪后的电成像测井灰度图像进行处理,得到灰度直方图,进而确定阈值,将骨架和孔隙区分;
b采用网络覆盖法求取骨架的分维数;
使用半径为δ的圆或边长为δ的正方形去覆盖图像,采用公式(1)计算分维数:
式中:D为分维数;a为常数,N(δ)为半径为δ时所覆盖图像的圆的个数;
c根据步骤a1所确定的火山岩岩性类型,得到每个深度段的分维数,进而确定不同岩性的分维数分布范围,建立岩性识别标准;
d导入未取心井段的成像测井数据,重复步骤a和步骤b,得到它们的分维数;
e根据步骤c所确定的不同岩性分维数分布范围,确定步骤d未取心井段的成像测井数据所对应的岩性。
步骤a1中:根据肉眼观察岩石薄片确定取心井段火山岩的岩性类型。
步骤a3中:得到灰度直方图后,根据最后一个峰的左波谷确定阈值。
步骤a3中:将骨架灰度置为255,将孔隙灰度置为0。
步骤a3中:所述阈值为75。
步骤b中:计算的分维数D的最大值和最小值分别为2和1。
本发明的有益技术效果是:
本发明通过对电成像测井的数据处理和分形维数计算等,可实现未取心井段岩性的自动识别,在一定程度上克服了传统的肉眼观察法带来的不确定性,提高岩性识别精度和速率。同时,本发明还能为其它岩类的测井识别带来新思路。
附图说明
下面结合附图与具体实施方式对本发明作进一步的说明:
图1为本发明提供的一种基于电成像测井分形维数的火山岩岩性识别方法技术流程图;
图2是导入的流纹岩原始电成像测井图像;
图3是对图2进行中值滤波后的流纹岩电成像测井图像;
图4是对图3进行开启运算后的流纹岩电成像测井图像;
图5是对图4进行阈值分割后的流纹岩电成像测井图像;
图6是对已取心井段确定岩性的电成像测井分维数的统计图。
具体实施方式
本发明构建了一个新的电成像测井定量处理方法,具体地说是提供一种基于电成像测井分形维数的火山岩岩性识别方法,该方法通过对电成像测井的数据处理和分形维数计算等,以实现岩性的自动识别,提高识别精度和速率。
本发明充分发挥了电成像测井分辨率高、能对岩层纹理进行精细刻画的特点,又克服了传统基于电成像测井定性识别岩性的不确定性,为火山岩储层的岩性识别提供了新思路。
一种基于电成像测井分形维数的火山岩岩性识别方法,包括以下步骤:
a导入取心井段的电成像测井数据,对电成像测井数据进行预处理,消除噪声;
a1根据肉眼观察、岩石薄片等确定取心井段火山岩的岩性类型,应用Matlab软件平台,分别读入它们所对应深度段的电成像测井数据;
a2先采用中值滤波算法消除电成像测井数据整体图像中的大部分噪声,再采用形态学滤波的开运算进一步将细小的噪声消除,突出岩石骨架和孔隙信息;
a3对去噪后的电成像测井灰度图像进行处理,得到图像的灰度直方图,根据最后一个峰的左波谷确定阈值,进而将骨架和孔隙区分出来,将骨架的灰度置为255,将孔隙的灰度置为0。
b采用网络覆盖法求取骨架的分维数。
使用半径为δ的圆去覆盖图像,采用公式(1)计算分维数:
式中:D为分维数;a为常数,N(δ)为半径为δ时所覆盖图像的圆的个数。
c根据步骤a1所确定的火山岩岩性类型,得到每个深度段所对应的分维数,进而确定不同岩性的分维数分布范围,建立岩性识别标准。
d导入未取心井段的成像测井数据,重复步骤a和步骤b,得到它们的分形维数。
e根据步骤c所确定的不同岩性分维数分布范围,确定步骤d未取心井段的成像测井数据所对应的岩性。
上述步骤中:孔隙和骨架的阈值一般为75。计算的分维数D的最大值和最小值分别为2和1。
下面结合附图对本发明作进一步说明。
图1是本发明一种基于电成像测井分形维数的火山岩岩性识别方法的技术流程图,主要包括导入取心井段成像测井数据、测井数据滤波与信号增强、阈值分割及骨架信息提取,确定不同岩性的分维数范围和未取心井段的分维数计算,以及导入未取心井段成像测井数据并计算分维数进行岩性识别五部分,这五部分必不可少,且顺序不可颠倒。
图2是流纹岩的原始电成像测井图。
图3是应用本发明所提供的方法对图2所述的流纹岩成像测井图进行中值滤波处理得到的电成像测井图像,从图中可知,图像中的大部分噪声得到有效压制。
图4是应用本发明所提供的方法对图3所述的流纹岩成像测井图进行开启运算处理得到的电成像测井图像,从图中可知,图像中的细节噪声得到有效压制。
图5是应用本发明所提供的方法对图4所述的流纹岩成像测井图进行阈值分割处理得到的电成像测井图像,从图中可知,孔隙和骨架的信息得到有效分离。
图6是应用本发明所提供的方法对取心井确定岩性的所有成像测井图计算得到的分维数分布图,从图中可知,不同岩心的分维数有明显区别。流纹岩和英安岩的分维数范围为1.1至1.3,安山岩和玄武岩的分维数范围为1.3至1.5,凝灰岩的分维数范围为1.5至1.7,角砾岩的分维数范围为1.7至1.9。
Claims (6)
1.一种基于电成像测井分形维数的火山岩岩性识别方法,其特征在于包括以下步骤:
a导入取心井段的成像测井数据,对成像测井数据进行预处理,消除噪声;
a1确定取心井段火山岩的岩性类型,应用Matlab软件,分别读入它们所对应深度段的电成像测井数据;
a2先采用中值滤波算法消除整体图像中的大部分噪声,再采用形态学滤波的开运算进一步将细小的噪声消除,突出岩石骨架和孔隙信息;
a3对去噪后的电成像测井灰度图像进行处理,得到灰度直方图,进而确定阈值,将骨架和孔隙区分;
b采用网络覆盖法求取骨架的分维数;
使用半径为δ的圆去覆盖图像,采用公式(1)计算分维数:
式中:D为分维数;a为常数,N(δ)为半径为δ时所覆盖图像的圆的个数;
c根据步骤a1所确定的火山岩岩性类型,得到每个深度段的分维数,进而确定不同岩性的分维数分布范围,建立岩性识别标准;
d导入未取心井段的成像测井数据,重复步骤a和步骤b,得到它们的分维数;
e根据步骤c所确定的不同岩性分维数分布范围,确定步骤d未取心井段的成像测井数据所对应的岩性。
2.根据权利要求1所述的一种基于电成像测井分形维数的火山岩岩性识别方法,其特征在于,步骤a1中:根据肉眼观察岩石薄片确定取心井段火山岩的岩性类型。
3.根据权利要求1所述的一种基于电成像测井分形维数的火山岩岩性识别方法,其特征在于,步骤a3中:得到灰度直方图后,根据最后一个峰的左波谷确定阈值。
4.根据权利要求1所述的一种基于电成像测井分形维数的火山岩岩性识别方法,其特征在于,步骤a3中:将骨架灰度置为255,将孔隙灰度置为0。
5.根据权利要求1所述的一种基于电成像测井分形维数的火山岩岩性识别方法,其特征在于,步骤a3中:所述阈值为75。
6.根据权利要求1所述的一种基于电成像测井分形维数的火山岩岩性识别方法,其特征在于,步骤b中:计算的分维数D的最大值和最小值分别为2和1。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810398397.2A CN108830140B (zh) | 2018-04-28 | 2018-04-28 | 一种基于电成像测井分形维数的火山岩岩性识别方法 |
PCT/CN2018/125979 WO2019205725A1 (zh) | 2018-04-28 | 2018-12-30 | 一种基于电成像测井分形维数的火山岩岩性识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810398397.2A CN108830140B (zh) | 2018-04-28 | 2018-04-28 | 一种基于电成像测井分形维数的火山岩岩性识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108830140A true CN108830140A (zh) | 2018-11-16 |
CN108830140B CN108830140B (zh) | 2020-06-16 |
Family
ID=64155770
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810398397.2A Expired - Fee Related CN108830140B (zh) | 2018-04-28 | 2018-04-28 | 一种基于电成像测井分形维数的火山岩岩性识别方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN108830140B (zh) |
WO (1) | WO2019205725A1 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019205725A1 (zh) * | 2018-04-28 | 2019-10-31 | 中国石油大学(华东) | 一种基于电成像测井分形维数的火山岩岩性识别方法 |
CN111058829A (zh) * | 2019-12-05 | 2020-04-24 | 中国矿业大学 | 基于图像处理的岩层分析方法 |
CN112213465A (zh) * | 2020-09-29 | 2021-01-12 | 西安石油大学 | 一种基于电成像测井恢复古熔岩流方位及火山机构原始形态的方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112801213A (zh) * | 2021-03-09 | 2021-05-14 | 成都理工大学 | 一种基于人工智能的岩石岩性快速分类识别方法和装置 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001020662A (ja) * | 1999-07-07 | 2001-01-23 | Fujita Corp | 安定性評価方法 |
CN101950359A (zh) * | 2010-10-08 | 2011-01-19 | 郝红卫 | 一种岩石种类的识别方法 |
CN102053268A (zh) * | 2010-11-18 | 2011-05-11 | 中国石油天然气股份有限公司 | 一种蚀变火山岩岩性识别方法及装置 |
CN102071928A (zh) * | 2009-11-25 | 2011-05-25 | 中国石油天然气股份有限公司 | 一种三维空间火山岩岩性识别方法 |
CN103325117A (zh) * | 2013-06-17 | 2013-09-25 | 中国石油天然气集团公司 | 一种基于matlab的岩心图像处理方法及系统 |
CN103353462A (zh) * | 2013-06-17 | 2013-10-16 | 中国石油大学(华东) | 一种基于核磁共振成像的岩石非均质性定量评价方法 |
US20150120197A1 (en) * | 2013-10-30 | 2015-04-30 | Chevron U.S.A. Inc. | Method of Transforming Reservoir Properties to a Seismic Attribute for Hydrocarbon and Lithology Identification |
CN105717551A (zh) * | 2014-12-05 | 2016-06-29 | 核工业北京地质研究院 | 一种火山岩型铀矿盲矿空间定位方法 |
CN106707340A (zh) * | 2016-12-13 | 2017-05-24 | 中国石油天然气股份有限公司大港油田分公司 | 一种火山岩岩相的预测方法 |
CN107092036A (zh) * | 2017-04-10 | 2017-08-25 | 中国科学院大学 | 一种基于储层真电阻率反演的火山岩地层流体识别方法及系统 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7174254B2 (en) * | 2004-05-14 | 2007-02-06 | Leroy Ellis | Mud gas isotope logging interpretative process utilizing mixing lines in oil and gas drilling operations |
CN104574346A (zh) * | 2013-10-23 | 2015-04-29 | 核工业北京地质研究院 | 一种光学遥感图像分解算法 |
CN103775071B (zh) * | 2014-01-15 | 2016-08-24 | 四川大学 | 采动煤岩体裂隙演化的测量方法 |
CN108830140B (zh) * | 2018-04-28 | 2020-06-16 | 中国石油大学(华东) | 一种基于电成像测井分形维数的火山岩岩性识别方法 |
-
2018
- 2018-04-28 CN CN201810398397.2A patent/CN108830140B/zh not_active Expired - Fee Related
- 2018-12-30 WO PCT/CN2018/125979 patent/WO2019205725A1/zh active Application Filing
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001020662A (ja) * | 1999-07-07 | 2001-01-23 | Fujita Corp | 安定性評価方法 |
CN102071928A (zh) * | 2009-11-25 | 2011-05-25 | 中国石油天然气股份有限公司 | 一种三维空间火山岩岩性识别方法 |
CN101950359A (zh) * | 2010-10-08 | 2011-01-19 | 郝红卫 | 一种岩石种类的识别方法 |
CN102053268A (zh) * | 2010-11-18 | 2011-05-11 | 中国石油天然气股份有限公司 | 一种蚀变火山岩岩性识别方法及装置 |
CN103325117A (zh) * | 2013-06-17 | 2013-09-25 | 中国石油天然气集团公司 | 一种基于matlab的岩心图像处理方法及系统 |
CN103353462A (zh) * | 2013-06-17 | 2013-10-16 | 中国石油大学(华东) | 一种基于核磁共振成像的岩石非均质性定量评价方法 |
US20150120197A1 (en) * | 2013-10-30 | 2015-04-30 | Chevron U.S.A. Inc. | Method of Transforming Reservoir Properties to a Seismic Attribute for Hydrocarbon and Lithology Identification |
CN105717551A (zh) * | 2014-12-05 | 2016-06-29 | 核工业北京地质研究院 | 一种火山岩型铀矿盲矿空间定位方法 |
CN106707340A (zh) * | 2016-12-13 | 2017-05-24 | 中国石油天然气股份有限公司大港油田分公司 | 一种火山岩岩相的预测方法 |
CN107092036A (zh) * | 2017-04-10 | 2017-08-25 | 中国科学院大学 | 一种基于储层真电阻率反演的火山岩地层流体识别方法及系统 |
Non-Patent Citations (6)
Title |
---|
ZHOU W.等: "Fractal analysis of fracture damage propagation of rock", 《TSINGHUA SCIENCE AND TECHNOLOGY》 * |
张强: "基于模拟退火的数字岩心重构技术及导电性数值模拟研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
张旭 等: "分形理论在火山岩地层抗钻特性评价中的应用", 《武汉大学学报(理学版)》 * |
牟丹: "辽河盆地中基性火成岩测井岩性识别方法研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
胡刚: "基于成像测井的火山岩结构分析", 《科技信息》 * |
陆诗阔 等: "变质岩储层岩性及裂缝测井识别方法研究进展", 《特种油气藏》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019205725A1 (zh) * | 2018-04-28 | 2019-10-31 | 中国石油大学(华东) | 一种基于电成像测井分形维数的火山岩岩性识别方法 |
CN111058829A (zh) * | 2019-12-05 | 2020-04-24 | 中国矿业大学 | 基于图像处理的岩层分析方法 |
CN111058829B (zh) * | 2019-12-05 | 2021-06-25 | 中国矿业大学 | 基于图像处理的岩层分析方法 |
CN112213465A (zh) * | 2020-09-29 | 2021-01-12 | 西安石油大学 | 一种基于电成像测井恢复古熔岩流方位及火山机构原始形态的方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2019205725A1 (zh) | 2019-10-31 |
CN108830140B (zh) | 2020-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108830140A (zh) | 一种基于电成像测井分形维数的火山岩岩性识别方法 | |
CN103592690B (zh) | 基于电成像测井孔隙度谱信息自动识别储层裂缝的方法 | |
CN102262041B (zh) | 一种基于多谱孔隙结构分析的饱和度确定方法 | |
CN110823962B (zh) | 一种滑坡体的三维成像方法及系统 | |
CN103135131A (zh) | 一种针对裂缝性储层预测的解释装置 | |
CN105986811A (zh) | 基于常规测井资料五维度岩性判识方法 | |
CN105114067A (zh) | 一种岩性测井相的方法 | |
CN106014397B (zh) | 一种成像测井缝洞储层参数提取方法 | |
CN108072915A (zh) | 一种识别碳酸盐岩颗粒滩相的方法 | |
CN110363299A (zh) | 面向露头岩层分层的空间案例推理方法 | |
Wang et al. | Fragmentation calculation method for blast muck piles in open-pit copper mines based on three-dimensional laser point cloud data | |
CN115390155A (zh) | 针对非均质储层的测井解释方法、装置、电子设备及介质 | |
CN108035710B (zh) | 基于数据挖掘划分深层岩石地质相的方法 | |
CN111626377B (zh) | 一种岩相识别方法、装置、设备及存储介质 | |
CN106468171A (zh) | 多信息耦合的递进式砂砾岩岩相判别方法 | |
CN104977602A (zh) | 一种地震数据采集施工的控制方法及装置 | |
CN116201535B (zh) | 油气藏目标井标志地层自动划分方法、装置及设备 | |
CN116466408A (zh) | 一种基于航磁数据的人工神经网络超基性岩识别方法 | |
CN111708085A (zh) | 一种基于波形分离的孔洞探测增强方法和装置 | |
CN108181665B (zh) | 裂缝的确定方法和装置 | |
CN105137495A (zh) | 一种油气检测的方法及系统 | |
CN109598049B (zh) | 钻孔岩石裂隙发育程度及区域岩石裂隙发育规律的方法 | |
CN105953766A (zh) | 基于岩芯图像的砾石磨圆度测量方法 | |
CN105824042A (zh) | 基于照明能量最优的最大纵向距设计方法 | |
CN111561266B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200616 |
|
CF01 | Termination of patent right due to non-payment of annual fee |