一种基于导航金字塔的张量投票处理方法技术

一种基于导航金字塔的张量投票处理方法包括以下步骤：首先，获取地震属性切片像，将其转换为灰度图像后进行导航金字塔处理，增强裂缝特征；然后将导航金字塔处理后的图像进行张量编码，也就是将数据以一个二阶对称非负的实数矩阵表示；将编码后的图像进行张量投票，其投票包括稀疏投票和稠密投票，实现裂缝等相关特征的检测；投票完成后进行张量分解，得到裂缝及溶洞边界特征图，同时在此基础上进一步投票，通过张量分析得到交叉点特征图，最终将裂缝及溶洞边界特征图和交叉点特征图叠置可视化。本发明专利技术主要特点是有效抑制噪声，增强裂缝的线性特征，提高检测的准确性，多尺度、多方向的刻画缝洞，突出裂缝之间和缝洞之间的交点。

【技术实现步骤摘要】
一种基于导航金字塔的张量投票处理方法
本专利技术属于地球物理勘探
，涉及一种基于导航金字塔的张量投票处理方法，有效进行裂缝及溶洞边界特征提取。
技术介绍
随着数学、物理学、计算机技术等学科的飞速发展，物探方法以及设备仪器也有了新的发展基础，特别是数字图像处理技术在工程物探中的应用范围逐步扩大，应用在裂缝识别提取的方法也层出不穷。目前对裂缝识别提取的方法有阈值分割法、边缘检测方法等。阈值分割法一般只适用于背景灰度一致、光照均匀和对比度较高的图像处理，而一幅地震属性切片图像中的边缘具有复杂且类型各异的特点，这些边缘所包含的信息是未知的且很多由噪声产生，难以精细刻画裂缝及溶洞边界。近年来，基于边缘检测的裂缝识别方法不断改进，解决了该方法对噪声过于敏感的缺点，为消除伪边缘并且能够准确清晰的识别裂缝，再此基础上也引用了数学形态学处理以及分形理论，此方法在近年来越来越趋向成熟，对裂缝检测具有重要的借鉴价值，但对于展现了裂缝的空间展布形态千变万化的地下复杂构造的地震属性切片图像来说，该方法对边缘特征的增强效果十分有限，特别是断裂两侧出现模糊现象时，无法识别边界内部梯度变化不大的裂缝，边界精准度达不到要求。上述单项技术都无法多尺度多方向对裂缝溶洞进行精细刻画，在实际生产中，如何利用地震属性切片数据准确的识别提取地下各尺度裂缝的空间展布以及缝洞之间的连通性可视化，并是地震储层描述综合研究亟待解决的问题。
技术实现思路
为解决上述问题，本专利技术的目的在于提供一种基于导航金字塔的张量投票处理方法，该方法能够利用叠后地震属性数据增强边界特征，多尺度多方向刻画裂缝，有效提取裂缝及溶洞边界特征。为达到以上技术目的，本专利技术提供以下技术方案。通过导航金字塔将地震图像进行多尺度多方向分解重构后进行张量投票，可同时处理多种结构特征，有效去除噪声，将视觉认知以计算机可视化展现出来，精细刻画断裂，该方法依次包括以下步骤：(1)获取地震属性切片数据，其转换为灰度图像；(2)对灰度图像进行导航金字塔处理，增强边界特征；(3)将导航金字塔处理后的图像中每个像素点编码为张量T，即将数据以一个二阶对称非负的实数矩阵表示。(4)张量投票。选择合适的投票域，将编码后的图像进行稀疏张量投票，将稀疏投票后的结果作为输入进行稠密投票。(5)投票完成后进行叠加分析，计算特征显著性。(6)步骤(5)中计算的λ1-λ2的值对应到图像中作为每个像素点值形成裂缝及边界概率图，在此基础上取裂缝及边界的极大值点，重复步骤(4)和步骤(5)，得到λ2的值对应到图像中作为像素点值形成交点概率图。(7)将裂缝及溶洞边界特征图和交叉点特征图叠置可视化。所述步骤(1)中，需要将地震属性切片数据导出为RGB图像，然后根据公式(1)将RGB图像转换为灰度图像。Gray＝0.29900*R+0.58700*G+0.11400*B(1)所述步骤(2)中的导航金字塔处理主要包括三个步骤：首先是利用极坐标滤波器进行多尺度分解，极坐标滤波器公式如下：其中，Ka和Kb分别为滤波器锥形区域的起始和截止波数。其次是使用方向可控滤波器进行多方向分解，方向可控滤波器具有任意旋转功能，它在每一层都可产生多个方向可控的子带集合，任何一个方向都可使用若干个基滤波器的线性组合而成，它对应的滤波基在频域上相互之间有一定重叠的滤波器。最后寻找最优方向进行插值重构，为得到精确的插值函数，在进行插值重构时需要准确知道所需基滤波器的数量，才能得到每个基滤波器的权重。所述步骤(3)中，是将导航金字塔处理后的图像数据编码为张量T，当图像数据不包含方向信息时，像素编码为球张量Ixy表示该位置像素值，其特征值λ1＝λ2＝1；当图像数据包含方向信息时，像素编码为棒张量其特征值分别为λ1＝1，λ2＝0。编码后的图像为张量图I。所述步骤(4)中的张量投票过程如下：投票之前先确定基本投票域，假设在二维笛卡尔坐标系中，投票者在原点O处，接受者位于点P处，投票点O处的张量主方向为Y轴方向。为获得在点P处投票的方向，过点O和P，作与X轴相切的圆，该密切圆在投票点O的径向向量CO平行于投票张量的主方向，由于通过密切圆的传播路线能很好地保持了曲率，则在接受点P的径向向量CP应平行于O向P所投选票的张量的主方向,这是因为通过密切圆的传播路径保持了曲率。其棒形投票域衰减函数为其球型投票域衰减函数为其中，s为OP的弧长，k为弧线的曲率，σ为投票尺度，C控制衰减程度，是投票尺度σ的函数，棒张量投票公式如下：其中，θ为OP与x轴正方向的夹角，DF为衰减函数，控制不同投票位置和方向上张量投票的强度。球张量投票公式如下：其中，Rθ′表示将张量Vstick旋转到方向上的旋转矩阵，即为Vstick最大的特征值所对应的特征向量，θ′表示x轴的旋转角。所述步骤(5)中的叠加分析，是将步骤(4)中各位置接收到的张量选票结果进行叠加，如下式TP＝T0+∑Vstick(p)+∑Vball(p)(8)将投票后的张量分解棒张量成分和球张量成分，即其中，λ1和λ2为张量T的特征值，为对应的特征向量，为棒张量分量，为球张量分量，λ2为球张量分量的显著性大小，λ1-λ2为棒张量分量的显著性大小。若λ1-λ2越大，则表示该点位于曲线上的可能性越大，将其对应到图像中作为每个像素点的值，则得到一幅表征裂缝出现在某个位置可能性的概率图。所述步骤(6)中，取步骤(5)中所得的概率图的局部极大值点，将其作为输入，重复步骤(4)和步骤(5)，得到球张量分量的显著性图，即交点图。所述步骤(7)中，裂缝及溶洞边界特征图和交叉点特征图叠置实现可视化。附图说明图1为本专利技术的一种基于导航金字塔的张量投票处理方法的流程框图；图2为本专利技术的导航金字塔分解重构步骤图；图3为本专利技术的张量投票过程示意图；图4为本专利技术的张量投票方法的测试效果示意图具体实施方式下面结合附图和实施例，对本专利技术的实施方式作进一步详细描述。如图1所示，本专利技术提供了一种基于导航金字塔的张量投票处理方法，包括以下步骤：(1)获取地震属性切片数据，其转换为灰度图像；(2)对灰度图像进行导航金字塔处理，增强边界特征；(3)将导航金字塔处理后的图像中每个像素点编码为张量T，即将数据以一个二阶对称非负的实数矩阵表示。(4)张量投票。选择合适的投票域，将编码后的图像进行稀疏张量投票，将稀疏投票后的结果作为输入进行稠密投票。(5)投票完成后进行叠加分析，计算特征显著性。(6)步骤(5)中计算的λ1-λ2的值对应到图像中作为每个像素点值形成裂缝及边界概率图，在此基础上取裂缝及边界的极大值点，重复步骤(4)和步骤(5)，得到λ2的值对应到图像中作本文档来自技高网...

【技术保护点】
1.一种基于导航金字塔的张量投票处理方法，其特征在于，包括以下步骤：/n(1)获取地震属性切片数据，其转换为灰度图像；/n(2)对灰度图像进行导航金字塔处理，增强边界特征；/n(3)将导航金字塔处理后的图像中每个像素点编码为张量T，即将数据以一个二阶对称非负的实数矩阵表示；/n(4)张量投票，选择合适的投票域，将编码后的图像进行稀疏张量投票，将稀疏投票后的结果作为输入进行稠密投票；/n(5)投票完成后进行叠加分析，计算特征显著性；/n(6)步骤(5)中计算的λ

【技术特征摘要】
1.一种基于导航金字塔的张量投票处理方法，其特征在于，包括以下步骤：
(1)获取地震属性切片数据，其转换为灰度图像；
(2)对灰度图像进行导航金字塔处理，增强边界特征；
(3)将导航金字塔处理后的图像中每个像素点编码为张量T，即将数据以一个二阶对称非负的实数矩阵表示；
(4)张量投票，选择合适的投票域，将编码后的图像进行稀疏张量投票，将稀疏投票后的结果作为输入进行稠密投票；
(5)投票完成后进行叠加分析，计算特征显著性；
(6)步骤(5)中计算的λ1-λ2的值对应到图像中作为每个像素点值形成裂缝及边界概率图，在此基础上取裂缝及边界的极大值点，重复步骤(4)和步骤(5)，得到λ2的值对应到图像中作为像素点值形成交点概率图；
(7)将裂缝及溶洞边界特征图和交叉点特征图叠置可视化。


2.如权利要求1所述的一种基于导航金字塔的张量投票处理方法，其特征在于，所述步骤(1)中，地震属性切片数据通过iloop软件获取并导出为RGB图像，并按照下式将图像转换为灰度图像Gray＝0.29900*R+0.58700*G+0.11400*B。


3.如权利要求1所述的一种基于导航金字塔的张量投票处理方法，其特征在于，所述步骤(2)中，导航金字塔处理包括多尺度、多方向分解和图像重构，需要对分解的级数和方向的数量进行设置，多尺度分解是将图像降采样且在频率域内根据下式进行：



其中，L(Ka,Kb,K)是极坐标滤波器函数，Ka和Kb分别为滤波器锥形区域的起始和截止波数；
多方向分解是利用二维高斯函数进行多方向分解，通过0°和9...

【专利技术属性】
技术研发人员：黄旭日，崔晓庆，徐云贵，胡叶正，曹卫平，唐静，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：四川;51