一种基于Contourlet变换的地震多属性融合方法技术

本发明专利技术属于地震资料油气预测、油气资源勘探及信号处理应用领域，提供了一种基于Contourlet变换的地震多属性融合方法，通过对信号的多分辨率、多方向性、局部化特性及各向异性的表示，提高描述属性信息的有效性，达到充分利用现有的各种属性信息，提高地震目标标定的准确性，从而能够更准确地进行储层预测及流体识别。

【技术实现步骤摘要】

本专利技术属于地震资料油气预测、油气资源勘探及信号处理应用领域，是一种利用Contourlet变换实现地震多属性融合的方法，采用本专利技术方法得到的地震属性融合结果信息量丰富、分辨率高，且消除了冗余信息，从而降低了单一地震属性油气预测的多解性，为降低油气勘探开发风险做出科学合理的指导。
技术介绍
近年来，地震属性技术发展迅速，已广泛应用于地层层序分析、油藏特征描述以及油藏动态检测等各个领域，成为了油藏地球物理的核心部分。地震属性也从早期的振幅属性发展到目前常用的数百种。因此，利用地震属性来研究地下储层的物性，乃至流体识别具有很大的前景。然而，每一种地震属性都只对某些地质特征敏感，利用各种地震属性信息解释复杂的地质问题时，由于地质因素、观测条件、测量精度、处理方式和主观因素等限制因 素，使用单一地震属性来预测储层会产生严重的多解性，缺乏可信度。并且多种地震属性信息参数按各自的方法原理和特征变化进行解释，有可能产生相互矛盾的结果。因此，开展地震多属性融合技术研究有利于降低多解性，能够更准确地进行地下储层的预测及流体识另O，为降低油气勘探开发风险做出科学合理的指导，是目前国外许多石油公司及科研机构的研究重点之地震多属性融合就是利用计算机自动实现两种或两种以上地震属性进行综合的有效手段，是地震属性分析、综合解释必备手段。纵观国内外现状，多属性综合或融合的基本理论和方法大多借用了多源信息融合理论及模式。目前，常用的融合手段和方法有(I)多属性或解释结果的人工叠合方法该方法属于低级融合，精度不高，需要人工配准。(2)基于属性切片的图像融合方法该方法为一种较好的计算机自动融合方法，主要思想来源于多光谱遥感图像处理技术，目前主要采用RGB-HIS颜色信息的2D信息融合。该方法的主要难度和影响精度的因素是属性源的配准问题。(3)基于地震属性数据的融合该方法是在属性数据上(1D、2D或3D)进行，融合效果和精度显然优于已经成像后基于图像的融合。目前，利用小波多分辨分析可以实现地震多属性融合，根据Mallat给出的多尺度分解与重构算法，将多种经过标准化处理后的地震属性在不同尺度上进行分解，分解到粗尺度上的信号称为平滑信号(低频分量)，分解到细尺度上的信号称为细节信号(高频分量)，其细节特征按分辨率的大小映射到各自的塔式结构中，本着在相同分辨率下进行信息融合的原则，在对应层不同的频率段进行融合。通常情况下，低频分量采用加权平均融合规贝U，对高频分量则采用区域方差协方差融合或加权融合准则，最终得到融合后的小波金字塔。对融合所得的小波金字塔进行小波逆变换，所得到的重构信号即为地震属性融合结果。基于小波多分辨分析的地震多属性融合方法，可以很好的描述信号的空域和时域的局部化特征，能够高效地对一维分段连续信号进行分析，但小波变换在一维所具有的特性并不能简单地推广到二维或更高维。由一维小波张成的可分离小波只具有有限的方向，是各向同性的，无法精确地表达信号边缘的方向，也不能“最优”表示含“线”或“面”奇异的二维信息，缺乏平移不变性和方向敏感性，使小波变换的应用受到了一定的限制。而事实上具有线或面奇异的函数在高维空间中非常普遍。例如，自然物体光滑边界使得自然信号的不连续性往往体现为光滑曲线上的奇异性，而并不仅仅是点奇异。可见，在高维情况下，小波变换并不能充分利用数据本身特有的几何特征，并不是最优的或者说“最稀疏”的函数表示算法。
技术实现思路
本专利技术的目的在于针对
技术介绍
中小波变换存在的问题，研究设计了，通过对信号的多分辨率、多方向性、局部化特性及各向异性的表示，提高描述属性信息的有效性，达到充分利用现有的各种属性信息，提高地 震目标标定的准确性，从而能够更准确地进行储层预测及流体识别。本专利技术解决小波变换存在的问题所采用的技术方案如下—种基于Contourlet变换的地震多属性融合方法,包括以下步骤①读入同一地质体的两种或多种二维地震属性数据。如均方根振幅属性、相干体切片属性等。②将步骤①中读入的地震属性数据进行预处理，包括地震属性的异常点剔除和标准化处理。目前，多信息储层预测可用的地震属性数据量大，在对众多属性融合之前，必须对地震属性进行预处理，主要有以下两个原因一是由于在采集或者提取属性数据的时候会出现一些提取错误的值，即异常点，这些离群异常值的存在对后面的标准化处理及属性融合处理有着致命的影响。二是由于不同地震属性的单位、量纲以及数值大小、变化范围不相同，某些属性特征往往会淹没在区域背景上，不能直接使用原始数据进行地质属性融合。如果，直接使用原始数据，就会突出绝对值大的属性，而压制绝对值小的属性，为避免这种不合理现象，需要将各个属性值换算到统一的规范尺度之下，即标准化处理(归一化处理)。首先，对读入的每种地震属性分别进行去噪处理，即去掉地震属性的异常点和毛刺。通常属性数据的异常点会远远大于或者小于平均值，可能会超出属性范围，本专利技术由公式(I)来判断异常点，满足该公式的点即为异常点，并用邻值代替法来处理异常点，用前面一个数据代替异常点。/-/ > (Xiax + Zmin )12(I)式中f、/、fmax和fmin分别为原始属性值、平均值、最大数据和最小数据。然后，对剔除异常点的地震属性数据进行标准化处理，即归一化处理。本专利技术采用极差正规化进行归一化处理，将属性的每个数据减去该属性的最小值，然后除以该属性的极差。标准化后的数据都在(Ti的范围内，最大值为1，最小值为O。具体变换公式为权利要求1.，其特征在于包括以下步骤 ①读入同一地质体的两种或多种二维地震属性数据； ②将步骤①中读入的地震属性数据进行预处理，包括地震属性的异常点剔除和标准化处理； ③将预处理后的地震属性数据分别进行多级Contourlet变换，并根据地震属性反应信息情况进行加权融合，得到融合低频分量和融合的不同尺度不同方向上的高频分量； ④进行Contourlet逆变换，即将步骤③中得到的融合低频分量和融合的不同尺度不同方向上的高频分量进行Contourlet重构，从而得到最后的地震属性融合结果。2.根据权利要求I所述的，其特征在于所述步骤②中的地震属性的异常点剔除和标准化处理包括以下步骤 a、由公式(I)来判断异常点，满足该公式的点即为异常点，并用邻值代替法来处理异常点，用前面一个数据代替异常点；3.根据权利要求I所述的，其特征在于所述步骤③中进行多级Contourlet变换得到多分辨率、多方向的Contourlet系数,然后根据地震属性反应的目标情况设定相应的权值，将各种地震属性的高频分量和低频分量分别进行加权融合，得到变换域的高频融合分量和低频融合分量。4.根据权利要求I所述的，其特征在于所述步骤③中的多级Contourlet变换利用拉普拉斯金字塔变换和方向滤波器组实现，首先将信号进行LP多尺度分解，然后将经过LP分解产生的高频子带输入DFB，DFB对LP分解后的高频子带进行方向分解，变换在LP分解后的低频子带上迭代进行，并逐渐将点奇异连成线形结构，从而捕获信号的轮廓，最终，Contourlet变换将信号分解到多尺度多方向的Contourlet域内，形成Contourlet系数,从而用于加权融合处理。5.根据权利要求4所述的，其特征在于所述拉普本文档来自技高网...

【技术保护点】
一种基于Contourlet变换的地震多属性融合方法，其特征在于包括以下步骤：①读入同一地质体的两种或多种二维地震属性数据；②将步骤①中读入的地震属性数据进行预处理，包括地震属性的异常点剔除和标准化处理；③将预处理后的地震属性数据分别进行多级Contourlet变换，并根据地震属性反应信息情况进行加权融合，得到融合低频分量和融合的不同尺度不同方向上的高频分量；④进行Contourlet逆变换，即将步骤③中得到的融合低频分量和融合的不同尺度不同方向上的高频分量进行Contourlet重构，从而得到最后的地震属性融合结果。

【技术特征摘要】

【专利技术属性】
技术研发人员：彭真明，谢成芳，李全忠，赵敏，黄振星，张伟，王圣川，张萍，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：