[发明专利]一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法

说明书

本发明涉及一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法，主要包括读取成像测井图片中图例像素点的RGB值，根据像素点的个数，插值生成一条归一化的电阻率曲线，使每一组不同的RGB值对应一个特定的电阻率值；逐点读取成像测井图像像素点的RGB值，用最小二乘法计算成像测井图像上每一个像素点的RGB值与图例中各像素点RGB值离差的平方和，离差的平方和越小代表这两个点的RGB值越接近；对成像测井图片上的每一个像素点都寻找图例中与其RGB值最接近的像素点，将图例上该点对应的电阻率值赋给成像测井图片中的对应点；对成像测井图像的电阻率值逐行求取算数平均值，使成像测井图像转变为一条一维的电阻率曲线；对这条电阻率曲线选用小波包分解，运用分解得到的低频成分判断层中岩相的韵律变化。

技术领域：

本发明涉及一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法。

背景技术：

微电导率成像测井曲线反映了地层微电导率的相对变化，这种变化可以转换成彩色图像或灰度图像，从而将井中测遇的地层序列转换成地层成像测井资料，浅色代表低电导率地层，深色代表高电导率地层，这样就可以精细地反映砂体内部的岩石类型、岩石结构、沉积构造、沉积韵律、沉积能量变化。根据成像测井图像上所反映出的特征建立解释模式及典型模板，判断沉积微相的方法在砂砾岩地区有着良好的效果，但由于成像测井图像自身精度的限制这种方法在粉细砂岩地区应用效果不佳，因此可以尝试针对粉细砂岩地区的沉积相研究提出一种新的研究方法。

小波变换的基函数随着尺度j的减小，其时频窗口宽度也减小，而相应的频域窗口宽度增大，也就是说，相应小波基函数的频域窗口随尺度减小而增大。所谓小波包，简单的说就是一个函数族。由他们构造出L²(R)的规范正交基库。从此库中可以选出L²(R)的许多组规范正交基，即小波包是对小波变换概念的推广。小波包变换实际上是将信号通过一系列的带通滤波器，对高频和低频频带做二进制划分，最后整个频带都被划分为均匀的频带；通过小波包分解，可以得到有尺度函数组成的子空间和小波函数组成的子空间(图1)。它克服了小波变换对高频信号频率分辨率低的缺陷小波包分析可以实现对信号更加精细的分析。小波包的空间分解为：

其一般表达式为：

这里，可以写为

当k＝0和m＝0时，子空间还原为其正交基就是小波基，即

发明内容：

本发明涉及一种利用成像测井图像识别粉砂岩与泥岩交替地层中岩相韵律变化的方法，主要包括读取成像测井图片中图例像素点的RGB值，根据像素点的个数，插值生成一条归一化的电阻率曲线，使每一组不同的RGB值对应一个特定的电阻率值；逐点读取成像测井图像像素点的RGB值，用最小二乘法计算成像测井图像上每一个像素点的RGB值与图例中各像素点RGB值离差的平方和，离差的平方和越小代表这两个点的RGB值越接近；对成像测井图片上的每一个像素点都寻找图例中与其RGB值最接近的像素点，将图例上该点对应的电阻率值赋给成像测井图片中的对应点；对成像测井图像的电阻率值逐行求取算数平均值，使成像测井图像转变为一条一维的电阻率曲线；对这条电阻率曲线选用小波包分解，运用分解得到的低频成分判断层中岩相的韵律变化。

附图说明

图1.小波包分解原理

图2.成像测井图像一维与处理后的成像测井一维数组对比

图3.处理后的成像测井一维数组小波包分解结果

具体实施方式：