专利名称：基于非下采样Contourlet变换的拟声波曲线构建方法
技术领域：
本发明涉及地震勘探数据处理技术，更具体地，涉及一种基于非下采样 Contourlet变换的拟声波曲线构建方法，用于构建高精度拟声波曲线以提高地震层位标定的准确性和地震反演的分辨率。
背景技术：
如何充分利用各种测井数据曲线来提高地震反演的分辨率和精度，一直是油气勘探开发地震研究工作努力的方向。长期以来，在地震反演中采用声波测井数据曲线来进行地震地质层位标定和波阻抗反演。然而，在很多情况下，由于井筒污染、储层胶结程度和孔隙度或其它非地层岩性因素影响，声波测井数据曲线中高频信息在很多情况下不能代表岩性变化，分不清地层剖面上的岩性，测井声波不能很好地反映储层和围岩的差异，导致岩性识别困难，从而造成测井数据曲线与地震剖面匹配较差，波阻抗反演结果与钻井地质不吻合，储层预测困难。近年来，在地震地质层位标定和波阻抗反演领域聚焦于拟声波的技术，目前很多生产单位都在试图用拟声波阻抗取代常规波阻抗反演。EPS软件作为独有的特色模块，在其反演系统中加入了拟声波曲线构建的功能。目前，常见的构建拟声波曲线的方法有(1)常规的测井数据曲线校正该方法仅简单地对声波曲线上因受到环境干扰而出现的异常值进行校正，并没有针对储层特征构建新的声波曲线。(2)经验公式或统计拟合法经验公式法是采用某种经验公式进行测井数据曲线间的转换，如应用Gardner公式将密度曲线转换为声波曲线，应用Faust公式将电阻率曲线转换为声波曲线等。这类方法虽然应用广泛，但因忽视了声波曲线中地层背景速度低频信息，所以理论上存在不足。统计拟合是将声波曲线与其他测井数据曲线通过交会图的方式(将两种测井数据在平面图上交会，根据交会点的坐标定出所求参数的数值或范围)统计拟合为一个关系式进行曲线转换。该方法同样没有考虑声波曲线中地层背景速度低频信息。(3)信息统计加权法 该方法将自然伽马(Gamma)、自然电位、电阻率等各种测井数据加权叠加为声波曲线以拟合声波数据，以突出声波测井、密度测井资料响应不明显的地质目标。但它仅在声波曲线基础上进行了简单的加减运算，而且在多数情况下只利用了单一的测井数据曲线进行重构。(4)基于小波变换的重构基本原理是对多种经过标准化处理的测井数据进行小波变换后，按分辨率的大小将细节特征映射到各自的塔式结构中，基于在相同分辨率下进行信息融合的原则，在对应层不同的频率段进行区域的特征选择。通常情况下都是以声波曲线为基础进行重构。因此，一般采用声波测井数据曲线的低频以及其它测井数据曲线的高频进行重构，再根据具体数据处理的要求，可以实现有针对性的地质目标增强，以利于观察和分析地质目标。
发明内容
本发明的基本构思是对ASCII格式的测井数据进行预处理(曲线校正和归一化)后，进行多级非下采样Contourlet变换和逆变换，从而得到最终构建的拟声波曲线。主要技术内容包括=ASCII格式测井数据的读入；测井数据的预处理；多级非下采样 Contourlet变换的正反演算；拟声波曲线的显示与对比。为了实现上述目的，根据本发明的一方面，提供了一种基于非下采样Contour let 变换的拟声波曲线构建方法，包括以下步骤读取ASCII格式的测井数据；对读取的测井数据曲线进行预处理，预处理包括曲线校正、数据标准化和归一化处理；对预处理后的测井数据曲线分别进行多级非下采样Contourlet变换，并根据测井数据曲线之间的匹配度进行加权融合，以生成高频融合曲线和低频融合曲线；将经过多级非下采样Contourlet变换得到的高频融合曲线和低频融合曲线进行非下采样Contourlet逆变换以生成拟声波曲线。根据本发明的一方面，所述构建方法还包括将拟声波曲线与原始声波曲线进行对比显示。根据本发明的一方面，所述测井数据曲线是声波数据曲线和非声波数据曲线，或者仅仅是非声波数据曲线。根据本发明的一方面，非声波数据曲线是Gamma曲线、电阻率曲线和密度曲线中的一个或多个。根据本发明的一方面，在预处理过程的曲线校正中，采用加权滑动平均法对测井数据进行滤波，其中，新采样的测井数据的权值高于先前采样的测井数据的权值。在对测井曲线进行标准化的步骤中，将目标井的标准层测井数据的频率交会图或直方图分别与关键井的标准层频率的测井数据的交会图或直方图进行对比来确定测井数据曲线之间的刻度误差，其中，将目标井的标准层的测井曲线的频率交会图或频率直方图与关键井的标准层的测井曲线的频率交会图或频率直方图的峰值之间的差值用作刻度误差的校正值。根据本发明的一方面，在非下采样Contourlet变换中采用非下采样金字塔和非下采样方向滤波器组，所述非下采样金字塔和非下采样方向滤波器组是二通道的迭代非下采样滤波器组。
通过下面结合附图进行的描述，本发明的上述和其他目的和特点将会变得更加清楚，其中图1是示出根据本发明实施例的非下采样Contourlet变换的拟声波曲线构建方法的流程图；图2是示出利用声波曲线、Gamma曲线和电阻率曲线构建的拟声波曲线的示意图；图3是示出利用密度曲线、Gamma曲线和电阻率曲线构建的拟声波曲线的示意图。
具体实施方式
根据本发明的拟声波曲线构建方法主要包括=ASCII格式测井数据的读入、测井数据的预处理、多级非下采样Contourlet变换的正反算和拟声波曲线的显示与对比等技术。以下，参照附图来详细说明本发明的实施方式。图1示出了根据本发明实施例的非下采样Contourlet变换的拟声波曲线构建方法的流程图。首先，在步骤100，读取测井数据，这里，通常采用ASCII格式的测井数据。接下来，在步骤110，对读取的测井数据进行预处理。在本实施例中，所述预处理可包括曲线的校正、标准化和归一化。具体说来，首先，对读入的每条测井数据曲线分别进行平滑滤波处理，以去掉测井数据曲线的统计起伏和毛刺干扰。本发明实施例采用加权滑动平均法对数据进行滤波，即把N个数据看成一个队列，队列的长度固定为N，每进行一次新的采样，把新采样结果放入队尾，而去掉排队首的一个数据。对于不同的数据应用不同的加权。通常，越接近待估计点的数据，其权值取得越大。这样每进行一次测量就可计算出新的滤波值，则第η次采样经滤波后的输出为以下的等式(1)
权利要求
1.一种基于非下采样Contourlet变换的拟声波曲线构建方法，包括以下步骤 读取测井数据；对读取的测井数据曲线进行预处理，预处理包括曲线校正、数据标准化和归一化处理；对预处理后的测井数据曲线分别进行多级非下采样Contourlet变换，并根据测井数据曲线之间的匹配度进行加权融合，以生成高频融合曲线和低频融合曲线；将经过多级非下采样Contourlet变换得到的高频融合曲线和低频融合曲线进行非下采样Contourlet逆变换以生成拟声波曲线。
2.如权利要求1所述的方法，还包括将拟声波曲线与测井数据中的原始声波曲线进行对比显示。
3.如权利要求1所述的方法，其中，所述测井数据曲线是声波数据曲线和非声波数据曲线，或者所述测井数据曲线仅包括非声波数据曲线。
4.如权利要求3所述的方法，其中，非声波数据曲线是Gamma曲线、电阻率曲线和密度曲线中的一个或多个。
5.如权利要求1所述的方法，其中，在预处理过程的曲线校正中，采用加权滑动平均法对测井数据进行滤波，其中，新采样的测井数据的权值高于先前采样的测井数据的权值。
6.如权利要求1所述的方法，其中，在对测井数据曲线进行标准化的步骤中，将目标井的标准层的测井数据的频率交会图或直方图分别与关键井的标准层的测井数据的频率交会图或直方图进行对比来确定测井数据曲线之间的刻度误差，其中，将目标井的标准层的测井数据的频率交会图或频率直方图与关键井的标准层的测井数据的频率交会图或频率直方图的峰值之间的差值用作刻度误差的校正值。
7.如权利要求1所述的方法，其中，在非下采样Contourlet变换中采用非下采样金字塔和非下采样方向滤波器组，所述非下采样金字塔和非下采样方向滤波器组是二通道的迭代非下采样滤波器组。
8.如权利要求1所述的方法，其中，测井数据是ASCII格式的数据。
全文摘要
提供了一种基于非下采样Contourlet变换的拟声波曲线构建方法，包括以下步骤读取测井数据；对读取的测井数据曲线进行预处理，预处理包括曲线校正、数据标准化和归一化处理；对预处理后的测井数据曲线分别进行多级非下采样Contourlet变换，并根据测井数据曲线之间的匹配度进行加权融合，以生成高频融合曲线和低频融合曲线；将经过多级非下采样Contourlet变换得到的高频融合曲线和低频融合曲线进行非下采样Contourlet逆变换以生成拟声波曲线。
文档编号G01V1/40GK102353985SQ20111014574
公开日2012年2月15日 申请日期2011年6月1日 优先权日2011年6月1日
发明者刘丽红, 张洞君, 彭真明, 邹文, 陶正喜, 陶韬, 黄东山 申请人:中国石油集团川庆钻探工程有限公司地球物理勘探公司