一种基于权重导向的地震属性保边滤波方法技术

本发明专利技术公开了一种基于权重导向的地震属性保边滤波方法。本发明专利技术中，通过结构算子初步计算属性中各点的边缘强度值，能够初步判断属性中各地质体的大致边缘位置，再利用最佳阈值法将边缘强度矩阵进行二值分割，并将分割结果作为后续滤波处理的先导矩阵，能够区分出属性中各地质体的强弱边缘，帮助后续进行滤波处理；其次将先导矩阵作为向导，对于地质体边缘处值跃变的区域进行非线性滤波进行保边处理，保留跃变趋势，在地质体内部值过渡平缓区域进行线性滤波处理，抹平值间差异，这种处理方式能够在压制噪声的前提下得到保留有地质体边缘信息的初始滤波处理结果；最后再对初始滤波结果进行误差分析，这种滞后处理方式能够对原始数据进行保真还原。始数据进行保真还原。

【技术实现步骤摘要】
一种基于权重导向的地震属性保边滤波方法


[0001]本专利技术属于地震勘探
，具体为一种基于权重导向的地震属性保边滤波方法。

技术介绍

[0002]随着地震勘探水平的不断提高，薄层储层内部的剩余油开发进入新的阶段，从薄层储层内部拾取出各地质体连续性边缘是当下亟待解决的重要问题，而基于地震平面属性进行地质体边界预测逐渐进入研究人员的视野。但是在地震平面属性中噪声是不可避免的存在，其严重影响边界信息预测的精度，而保持地质体边界信息的前提下压制平面属性数据中的噪声，对于地质体边界预测有着重要意义。
[0003]但是在实际应用中，地震平面属性地质体边缘处属性值具有跃变的趋势，在滤波处理中，此处主要影响滤波效果的是值权重；地质体内部则过渡趋势较为平缓，在滤波处理中，此处主要影响滤波效果的是空间距离权重。属性值的不同变化趋势需要不同类型的滤波方法进行滤波处理。相关研究表明，属性值过渡趋势较为平缓的区域，主要考虑空间距离权重的线性的滤波方法能够帮助减小值间差异，减小连续性噪声的影响；属性值跃变区域需要优先考虑属性值的影响，提高属性值在滤波过程中所占权重，相关非线性的滤波方法能够保留值的跃变趋势，减小非连续性噪声的影响。然而，现有的平面属性滤波技术方法并不具备这种特性，限制了地质体边界预测的精度。

技术实现思路

[0004]本专利技术的目的在于：为了解决上述提出的问题，提供一种基于权重导向的地震属性保边滤波方法。
[0005]本专利技术采用的技术方案如下：一种基于权重导向的地震属性保边滤波方法，所述基于权重导向的地震属性保边滤波方法包括以下步骤：
[0006]步骤S1：输入地震属性At，计算属性数据中各地质体的边缘强度S；
[0007]步骤S2：基于上述边缘强度数据S，使用最佳阈值法进行得到阈值T，并根据阈值T将属性数据中各点进行二值分割，得到先导矩阵Gd；
[0008]步骤S3：根据先导矩阵G，对平面属性At进行分类滤波处理，得到初始滤波数据At1；
[0009]步骤S4：基于初始滤波数据At1，根据误差常数σ进行误差分析，得到最终滤波数据At2。
[0010]在一优选的实施方式中，所述步骤S1包括以下步骤：
[0011]步骤S11：输入地震属性At并将其规约为0～255的灰度数据G，其每点的计算表达式为：
[0012][0013]式中，(x,y)横纵坐标，max(At(x,y))代表At中的最大值；
[0014]步骤S12：计算灰度数据G中各地质体的水平方向的边缘强度S
x
，其表达式为：
[0015][0016]式中，*表示卷积运算；
[0017]步骤S13：计算灰度数据G中各地质体的垂直方向的边缘强度S
y
，其表达式为：
[0018][0019]步骤S14：基于水平方向的边缘强度S
x
和垂直方向的边缘强度S
y
计算属性数据中各地质体的边缘强度S，其每点的计算表达式为：
[0020][0021]在一优选的实施方式中，所述步骤S2包括以下步骤：
[0022]步骤S21：基于属性数据中各地质体的边缘强度S，以最佳阈值分割法计算全局最优二值化阈值T；
[0023]步骤S22：根据全局最优二值化阈值T分割边缘强度数据S得到先导矩阵Gd，其每点的计算表达式为：
[0024][0025]在一优选的实施方式中，所述步骤S21包括以下步骤：
[0026]步骤S211：基于S的归一化图像直方图I，计算S的零阶累积矩A，其每个灰度值为k的计算表达式为：
[0027][0028]式中，I(i)表示值i在S中出现的次数，i,k为灰度值，i,k∈[0,255]；
[0029]步骤S212：基于S的归一化图像直方图I，计算S的一阶累积矩B，其每个灰度值为k的计算表达式为：
[0030][0031]步骤S213：基于S的归一化图像直方图I，计算k＝255时S的累积矩M(255)，其表达式为：
[0032][0033]步骤S214：基于A，B和M(255)计算S中每个灰度值的平均方差σ1，其表达式为：
[0034][0035]步骤S215：选取σ1中的最大值σ1(m)，此时的m为优选的全局最佳阈值T，其表达式为：
[0036][0037]在一优选的实施方式中，所述步骤S3包括以下步骤：
[0038]步骤S31：构建线性滤波模板L
n
；
[0039]步骤S32：构建非线性滤波模板NL
n
；
[0040]步骤S33：基于先导矩阵Gd，线性滤波模板L
n
和非线性滤波模板NL
n
对地震平面属性At进行分类滤波处理，具体为在Gd(x,y)＝1处进行非线性滤波处理，在Gd(x,y)＝0处进行线性滤波处理，得到初始滤波数据At1，At中每点的计算表达式为：
[0041][0042]式中，At(i,j)为计算中以At(x,y)为靶点，大小为n的滑动窗口矩阵。
[0043]在一优选的实施方式中，所述步骤S4包括以下步骤：
[0044]步骤S41：设定误差系数θ，针对地震平面属性数据At和初始滤波数据At1进行误差分析，并进行数据校正处理得到最终滤波结果At2，At2中每点的计算表达式为：
[0045][0046]综上所述，由于采用了上述技术方案，本专利技术的有益效果是：
[0047]本专利技术中，首先通过结构算子初步计算属性中各点的边缘强度值，能够初步判断属性中各地质体的大致边缘位置，再利用最佳阈值法将边缘强度矩阵进行二值分割，并将分割结果作为后续滤波处理的先导矩阵，能够区分出属性中各地质体的强弱边缘，帮助后续进行滤波处理；其次将先导矩阵作为向导，对于地质体边缘处值跃变的区域进行非线性滤波进行保边处理，保留跃变趋势，在地质体内部值过渡平缓区域进行线性滤波处理，抹平值间差异，这种处理方式能够在压制噪声的前提下得到保留有地质体边缘信息的初始滤波处理结果；最后再对初始滤波结果进行误差分析，这种滞后处理方式能够在一定程度上对原始数据进行保真还原。同时也能够通过检测地质体边缘并以此作为向导对属性中地质体内外信息进行区分处理，这种处理方式能够极大程度保留边缘、压制噪声，从而为后续地质体边缘信息提取、标定提供基础保障。
具体实施方式
[0048]为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0049]实施例：
[0050]一种基于权重导向的地震属性保边滤波方法，所述基于权重导向的地震属性保边滤波方法包括以下步骤：
[0051]步骤S1：输入本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于权重导向的地震属性保边滤波方法，其特征在于：所述基于权重导向的地震属性保边滤波方法包括以下步骤：步骤S1：输入地震属性At，计算属性数据中各地质体的边缘强度S；步骤S2：基于上述边缘强度数据S，使用最佳阈值法进行得到阈值T，并根据阈值T将属性数据中各点进行二值分割，得到先导矩阵Gd；步骤S3：根据先导矩阵G，对平面属性At进行分类滤波处理，得到初始滤波数据At1；步骤S4：基于初始滤波数据At1，根据误差常数σ进行误差分析，得到最终滤波数据At2。2.如权利要求1所述的一种基于权重导向的地震属性保边滤波方法，其特征在于：所述步骤S1包括以下步骤：步骤S11：输入地震属性At并将其规约为0～255的灰度数据G，其每点的计算表达式为：式中，(x,y)横纵坐标，max(At(x,y))代表At中的最大值；步骤S12：计算灰度数据G中各地质体的水平方向的边缘强度S
x
，其表达式为：式中，*表示卷积运算；步骤S13：计算灰度数据G中各地质体的垂直方向的边缘强度S
y
，其表达式为：步骤S14：基于水平方向的边缘强度S
x
和垂直方向的边缘强度S
y
计算属性数据中各地质体的边缘强度S，其每点的计算表达式为：3.如权利要求1所述的一种基于权重导向的地震属性保边滤波方法，其特征在于：所述步骤S2包括以下步骤：步骤S21：基于属性数据中各地质体的边缘强度S，以最佳阈值分割法计算全局最优二值化阈值T；步骤S22：根据全局最优二值化阈值T分割边缘强度数据S得到先导矩阵Gd，其每点的计算表达式为：4.如权利要求3所述的一种基于权重导向的地震属性保边滤波方法...

【专利技术属性】
技术研发人员：丁峰，邓琳平，张栋，
申请(专利权)人：四川中质鼎峰勘查技术有限公司，
类型：发明
国别省市：