CN103867196B - 一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法 - Google Patents
一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法 Download PDFInfo
- Publication number
- CN103867196B CN103867196B CN201410126102.8A CN201410126102A CN103867196B CN 103867196 B CN103867196 B CN 103867196B CN 201410126102 A CN201410126102 A CN 201410126102A CN 103867196 B CN103867196 B CN 103867196B
- Authority
- CN
- China
- Prior art keywords
- pixel
- imaging logging
- resistivity
- rgb value
- rhythm
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 18
- 230000033764 rhythmic process Effects 0.000 title claims abstract description 18
- 239000004575 stone Substances 0.000 title claims abstract description 9
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 13
- 230000000452 restraining effect Effects 0.000 claims description 4
- 230000005611 electricity Effects 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 239000011435 rock Substances 0.000 description 4
- 239000004576 sand Substances 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- QGWYLXFVPIMLDO-UHFFFAOYSA-N ethyl n-[4-[benzyl(2-phenylethyl)amino]-2-(2,4,5-trimethoxyphenyl)-1h-imidazo[4,5-c]pyridin-6-yl]carbamate Chemical class N=1C(NC(=O)OCC)=CC=2NC(C=3C(=CC(OC)=C(OC)C=3)OC)=NC=2C=1N(CC=1C=CC=CC=1)CCC1=CC=CC=C1 QGWYLXFVPIMLDO-UHFFFAOYSA-N 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法，主要包括读取成像测井图片中图例像素点的RGB值，根据像素点的个数，插值生成一条归一化的电阻率曲线，使每一组不同的RGB值对应一个特定的电阻率值；逐点读取成像测井图像像素点的RGB值，用最小二乘法计算成像测井图像上每一个像素点的RGB值与图例中各像素点RGB值离差的平方和，离差的平方和越小代表这两个点的RGB值越接近；对成像测井图片上的每一个像素点都寻找图例中与其RGB值最接近的像素点，将图例上该点对应的电阻率值赋给成像测井图片中的对应点；对成像测井图像的电阻率值逐行求取算数平均值，使成像测井图像转变为一条一维的电阻率曲线；对这条电阻率曲线选用小波包分解，运用分解得到的低频成分判断层中岩相的韵律变化。
Description
技术领域：
本发明涉及一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法。
背景技术：
微电导率成像测井曲线反映了地层微电导率的相对变化，这种变化可以转换成彩色图像或灰度图像，从而将井中测遇的地层序列转换成地层成像测井资料，浅色代表低电导率地层，深色代表高电导率地层，这样就可以精细地反映砂体内部的岩石类型、岩石结构、沉积构造、沉积韵律、沉积能量变化。根据成像测井图像上所反映出的特征建立解释模式及典型模板，判断沉积微相的方法在砂砾岩地区有着良好的效果，但由于成像测井图像自身精度的限制这种方法在粉细砂岩地区应用效果不佳，因此可以尝试针对粉细砂岩地区的沉积相研究提出一种新的研究方法。
小波变换的基函数随着尺度j的减小，其时频窗口宽度也减小，而相应的频域窗口宽度增大，也就是说，相应小波基函数的频域窗口随尺度减小而增大。所谓小波包，简单的说就是一个函数族。由他们构造出L2(R)的规范正交基库。从此库中可以选出L2(R)的许多组规范正交基，即小波包是对小波变换概念的推广。小波包变换实际上是将信号通过一系列的带通滤波器，对高频和低频频带做二进制划分，最后整个频带都被划分为均匀的频带；通过小波包分解，可以得到有尺度函数组成的子空间和小波函数组成的子空间(图1)。它克服了小波变换对高频信号频率分辨率低的缺陷小波包分析可以实现对信号更加精细的分析。小波包的空间分解为：
其一般表达式为：
这里，可以写为
当k＝0和m＝0时，子空间还原为其正交基就是小波基，即
发明内容：
本发明涉及一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法，主要包括读取成像测井图片中图例像素点的RGB值，根据像素点的个数，插值生成一条归一化的电阻率曲线，使每一组不同的RGB值对应一个特定的电阻率值；逐点读取成像测井图像像素点的RGB值，用最小二乘法计算成像测井图像上每一个像素点的RGB值与图例中各像素点RGB值离差的平方和，离差的平方和越小代表这两个点的RGB值越接近；对成像测井图片上的每一个像素点都寻找图例中与其RGB值最接近的像素点，将图例上该点对应的电阻率值赋给成像测井图片中的对应点；对成像测井图像的电阻率值逐行求取算数平均值，使成像测井图像转变为一条一维的电阻率曲线；对这条电阻率曲线选用小波包分解，运用分解得到的低频成分判断层中岩相的韵律变化。
附图说明
图1.小波包分解原理
图2.成像测井图像一维与处理后的成像测井一维数组对比
图3.处理后的成像测井一维数组小波包分解结果
具体实施方式：
步骤一，读取成像测井图片中图例像素点的RGB值，根据像素点的个数，插值生成一条归一化的电阻率曲线，使每一组不同的RGB值对应一个特定的电阻率值，其特征是选择成像测井图例像素点所做组成矩阵的任意一行，共n组RGB值，用R1表示，R1＝(xRi，xGi，xBi)，(i＝1，2，…n)；根据像素点的个数，差值生成一条由n点组成的电阻率曲线，其中这一行的第一个像素点对应的电阻率值为0，最后一个像素点对用的电阻率值为1；
步骤二，逐点读取成像测井图像像素点的RGB值，R2＝(xRjk，xGjk，xBjk)，用最小二乘法计算成像测井图像上每一个像素点的RGB值与图例中各像素点RGB值离差的平方和，离差的平方和越小代表这两个点的RGB值越接近；对成像测井图片上的每一个像素点都寻找图例中与其RGB值最接近的像素点，将图例上该点对应的电阻率值赋给成像测井图片中的对应点；
步骤三，对成像测井图像的电阻率值逐行求取算数平均值，使成像测井图像转变为一条一维的电阻率曲线，其特征是求取每一行电阻率的算数平均值，使成像测井图像转变为一维数组(图1)；
步骤四，选用db小波包对曲线进行8层分解，用分解后得到的低频成分判断地层韵律型，判断原则是若电阻率从左往右逐渐增大即为正韵律地层；若电阻率从左往右逐渐减小即为反韵律地层(图2)。
Claims (4)
1.一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法，主要包括读取成像测井图片中图例像素点的RGB值，根据像素点的个数，插值生成一条归一化的电阻率曲线，使每一组不同的RGB值对应一个特定的电阻率值；逐点读取成像测井图像像素点的RGB值，用最小二乘法计算成像测井图像上每一个像素点的RGB值与图例中各像素点RGB值离差的平方和，离差的平方和越小代表这两个点的RGB值越接近；对成像测井图片上的每一个像素点都寻找图例中与其RGB值最接近的像素点，将图例上该点对应的电阻率值赋给成像测井图片中的对应点；对成像测井图像的电阻率值逐行求取算数平均值，使成像测井图像转变为一条一维的电阻率曲线；对这条电阻率曲线选用小波包分解，运用分解得到的低频成分判断层中岩相的韵律变化；
所述的读取成像测井图片中图例像素点的RGB值，根据像素点的个数，插值生成一条归一化的电阻率曲线，使每一组不同的RGB值对应一个特定的电阻率值，其特征是选择成像测井图例像素点所做组成矩阵的任意一行，共n组RGB值，用R1表示，R1＝(xRi，xGi，xBi)，(i＝1，2，…n)；根据像素点的个数，插值生成一条由n点组成的电阻率曲线，其中这一行的第一个像素点对应的电阻率值为0，最后一个像素点对应的电阻率值为1。
2.根据权利要求1所述的一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法，所述的逐点读取成像测井图像像素点的RGB值，用最小二乘法计算成像测井图像上每一个像素点的RGB值与图例中各像素点RGB值离差的平方和，离差的平方和越小代表这两个点的RGB值越接近，其特征是读取成像测井静态图像数组R2＝(xRjk，xGjk，xBjk)，
3.根据权利要求1所述的一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法，所述的对成像测井图像的电阻率值逐行求取算数平均值，使成像测井图像转变为一条一维的电阻率曲线，其特征是求取每一行电阻率的算数平均值，使成像测井图像转变为1*h列的一维数组。
4.根据权利要求1所述的一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法，所述的对这条电阻率曲线选用小波包分解，运用分解得到的低频成分判断层中岩相的韵律变化，其特征是选用db小波包对曲线进行8层分解，用分解后得到的低频成分判断地层韵律型，判断原则是若电阻率从左往右逐渐增大即为正韵律地层；若电阻率从左往右逐渐减小即为反韵律地层。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410126102.8A CN103867196B (zh) | 2014-04-01 | 2014-04-01 | 一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410126102.8A CN103867196B (zh) | 2014-04-01 | 2014-04-01 | 一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103867196A CN103867196A (zh) | 2014-06-18 |
| CN103867196B true CN103867196B (zh) | 2019-03-22 |
Family
ID=50906158
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410126102.8A Expired - Fee Related CN103867196B (zh) | 2014-04-01 | 2014-04-01 | 一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103867196B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016144202A1 (ru) * | 2015-03-10 | 2016-09-15 | Алексей Алексеевич СУХАНОВ | Способ определения зрелости сапропелевого органического вещества нефтегазоматеринских толщ |
| CN106014397B (zh) * | 2016-05-13 | 2019-06-18 | 中国石油天然气集团公司 | 一种成像测井缝洞储层参数提取方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102033247A (zh) * | 2010-10-22 | 2011-04-27 | 中国石油天然气股份有限公司 | 一种逐点刻度电成像资料计算视地层水电阻率谱及参数的方法 |
| CN102053268A (zh) * | 2010-11-18 | 2011-05-11 | 中国石油天然气股份有限公司 | 一种蚀变火山岩岩性识别方法及装置 |
| CN102418518A (zh) * | 2011-04-12 | 2012-04-18 | 北京师范大学 | 神经网络模拟交会图识别油层水淹级别的方法 |
| CN103592690A (zh) * | 2013-10-24 | 2014-02-19 | 长江大学 | 基于电成像测井孔隙度谱信息自动识别储层裂缝的方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2440591C2 (ru) * | 2008-04-10 | 2012-01-20 | Шлюмбергер Текнолоджи Б.В. | Способ получения характеристик геологической формации, пересекаемой скважиной |
| GB2503010B (en) * | 2012-06-14 | 2018-04-18 | Reeves Wireline Tech Ltd | A method of processing geological log data |
-
2014
- 2014-04-01 CN CN201410126102.8A patent/CN103867196B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102033247A (zh) * | 2010-10-22 | 2011-04-27 | 中国石油天然气股份有限公司 | 一种逐点刻度电成像资料计算视地层水电阻率谱及参数的方法 |
| CN102053268A (zh) * | 2010-11-18 | 2011-05-11 | 中国石油天然气股份有限公司 | 一种蚀变火山岩岩性识别方法及装置 |
| CN102418518A (zh) * | 2011-04-12 | 2012-04-18 | 北京师范大学 | 神经网络模拟交会图识别油层水淹级别的方法 |
| CN103592690A (zh) * | 2013-10-24 | 2014-02-19 | 长江大学 | 基于电成像测井孔隙度谱信息自动识别储层裂缝的方法 |
Non-Patent Citations (2)
| Title |
|---|
| "利用成像测井资料进行沉积微相分析";彭跃 等;《石油天然气学报》;20090831;第31卷(第4期);第379-383页 |
| "基于微电阻率扫描成像测井的沉积微相识别";孙鲁平 等;《测井技术》;20090831;第33卷(第4期);第269-273页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103867196A (zh) | 2014-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Perez-Muñoz et al. | Wavelet transform analysis for lithological characteristics identification in siliciclastic oil fields | |
| Purkis et al. | Remote sensing of geomorphology and facies patterns on a modern carbonate ramp (Arabian Gulf, Dubai, UAE) | |
| Gao et al. | Microfacies and depositional environments of the Late Ordovician Lianglitage Formation at the Tazhong Uplift in the Tarim Basin of northwest China | |
| Feranie et al. | 3D fractal dimension and flow properties in the pore structure of geological rocks | |
| Li et al. | A new method for evaluating the pore structure complexity of digital rocks based on the relative value of fractal dimension | |
| CN107179550B (zh) | 一种数据驱动的地震信号零相位反褶积方法 | |
| Barnaby et al. | Outcrop analog for mixed siliciclastic–carbonate ramp reservoirs—stratigraphic hierarchy, facies architecture, and geologic heterogeneity: Grayburg Formation, Permian Basin, USA | |
| CN103867196B (zh) | 一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法 | |
| Scott et al. | Multi-scale image-based pore space characterisation and pore network generation: Case study of a north sea sandstone reservoir | |
| Milne et al. | Spectral and wavelet analysis of gilgai patterns from air photography | |
| Theologou et al. | An MICP-based pore typing workflow–core scale to log scale | |
| Abbas et al. | Lithofacies classification of carbonate reservoirs using advanced machine learning: a case study from a Southern Iraqi oil field | |
| Wang et al. | Coherence cube enhancement based on local histogram specification | |
| Larson et al. | Machine learning classification of Austin Chalk chemofacies from high-resolution x-ray fluorescence core characterization | |
| Lehane et al. | Analytical tools for quantifying the morphology of invertebrate trace fossils | |
| Berg et al. | 2-D pore architecture characterization of a carbonate mudrock reservoir: Insights from the mid-continent “Mississippi Lime” | |
| Purkis et al. | Visualizing lateral anisotropy in modern carbonates | |
| CN107678065B (zh) | 提高地震分辨率的保构造井控空间反褶积方法和装置 | |
| Gaci et al. | Two-dimensional multifractional Brownian motion-based investigation of heterogeneities from a core image | |
| Ferreira et al. | The importance of outcrop reservoir characterization in oil-industry facies modelling workflows-a case study from the Middle Jurassic of the Maciço Calcário Estremenho, Portugal | |
| Xiao et al. | Detecting sedimentary cycles using autocorrelation of grain size | |
| Alward et al. | Integrated Workflow for Petrophysical Evaluation in Heterogeneous Carbonate Reservoirs | |
| Ottinger et al. | Integration of Subsurface Data, Modern Analog Concepts, and History Matching Feedback to Model a Thick Carbonate Reservoir, Offshore Abu Dhabi | |
| CN113671565B (zh) | 一种针对巨厚储层气藏开发的地震多尺度储层预测方法 | |
| Chalabi et al. | Spatial Distribution of Sedimentary Facies on Modern Carbonate Platforms Interpreted via Remote Sensing, Northwest Celebes Sea, Malaysia |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into 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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190322 |