专利名称：一种碳酸盐岩内油水识别的方法及系统的制作方法
技术领域：
本发明关于地球物理勘探领域，特别是关于油气勘探开发领域，具体的讲是一种碳酸盐岩内油水识别的方法及系统。
背景技术：
近年来，碳酸盐岩逐渐成为油气勘探非常重要的一个新领域。我国的碳酸盐岩油气田通常埋藏很深，储层成因复杂，主要包括同生期和准同生期岩溶作用、表生岩溶作用和埋藏岩溶作用。而不同地质背景下和不同成因形成的储层类型多种多样，碳酸盐岩储层通常具有复杂的孔、洞和缝结构，非均值性强。在勘探中准确预测碳酸盐岩储层类型已属不易，而辨析碳酸盐岩储层中复杂的油气水分布关系更为困难。因为，和均值砂岩储层相比，大多数碳酸盐岩油藏没有统一的油水界面，即便是相邻很近的储集空间内充注的流体也可
能不同。因此，如何进行碳酸盐岩储层油水识别是本领域的一个技术难题。目前，油气田勘探开发中常用的烃类检测技术是利用地震反射(或折射)特征进行油气藏的判识，包括地震属性技术、AVO (Amplitude Versus Offset,振幅随偏移距的变化)技术、合成声波测井技术、纵横波速度比技术等，其中最具代表性是AVO技术。但是传统的AVO分析只与介质分界面两边的弹性参数有关。另外，传统的AVO技术在实际应用中存在不便在众多的三维地震剖面上依靠肉眼观察来寻找振幅增强或减弱现象并不完全可靠，容易造成遗漏或重复；在平面上无法成图，不便于分析地质现象在平面上的分布规律。

发明内容
为了克服现有技术存在的上述问题，本发明提出了一种碳酸盐岩内油水识别的方法及系统，通过在对地震资料做分频处理的基础上，求取近、远道地震振幅差值的方法来进行油水识别，既考虑了频率因素，也将地震剖面上定性的振幅变化现象转换为可以量化描述的振幅差值的信息，可以在平面和剖面上同时成图，更便于单个井点的油水识别和平面上地质规律的统计分析，充分发挥其钻前预测的优势。本发明的目的之一是，提供一种碳酸盐岩内油水识别的方法，所述的方法包括采集地震三维工区的地震纵波传输信号；对所述的地震纵波传输信号进行叠前处理；对叠前处理后的地震数据进行常规地震层位标定和解释；对叠前处理后的地震数据进行波阻抗反演；根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率；对所述的频率对应的地震数据进行归一化处理；根据所述的归一化处理的结果进行油水识别。本发明的目的之一是，提供了一种碳酸盐岩内油水识别的系统，所述的系统包括采集装置，用于采集地震三维工区的地震纵波传输信号；叠前处理装置，用于对所述的地震纵波传输信号进行叠前处理；层位标定装置，用于对叠前处理后的地震数据进行常规地震层位标定和解释；波阻抗反演装置，用于对叠前处理后的地震数据进行波阻抗反演；频率提取装置，用于根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率；归一化处理装置，用于对所述的频率对应的地震数据进行归一化处理；油水识别装置，用于根据所述的归一化处理的结果进行油水识别。本发明的有益效果在于，通过在对地震资料做分频处理的基础上，求取近、远道地震振幅差值的方法来进行油水识别，既考虑了频率因素，也将地震剖面上定性的振幅变化现象转换为可以量化描述的振幅差值的信息，可以在平面和剖面上同时成图，更便于单个井点的油水识别和平面上地质规律的统计分析，充分发挥其钻前预测的优势，预测结果更准确。


为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可
以根据这些附图获得其他的附图。图I为本发明实施例提供的一种碳酸盐岩内油水识别的方法的流程图；图2为本发明实施例提供的一种碳酸盐岩内油水识别的方法中步骤S102的具体流程图；图3为图I中的步骤S103的具体流程图；图4为本发明实施例提供的一种碳酸盐岩内油水识别的方法中步骤S105的具体流程图；图5为本发明实施例提供的一种碳酸盐岩内油水识别的系统的实施方式一的结构框图；图6为本发明实施例提供的一种碳酸盐岩内油水识别的系统的实施方式二的结构框图；图7为本发明实施例提供的一种碳酸盐岩内油水识别的系统的实施方式三的结构框图；图8为本发明实施例提供的一种碳酸盐岩内油水识别的系统的实施方式四的结构框图；图9为本发明实施例提供的一种碳酸盐岩内油水识别方法中全频段地震数据剖面；图10为本发明实施例提供的一种碳酸盐岩内油水识别方法中叠前反演波阻抗剖面；图11中，(a)为IOHz的近道集的振幅剖面；(b)为IOHz的远道集的振幅剖面；(c)是40Hz的近道集的振幅剖面；(d)是40Hz的远道集分量的振幅剖面；图12为40Hz的近、远道集的振幅差值剖面图；图13为40Hz的近、远道集的振幅差值平面图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。图I为本发明实施例提供的一种碳酸盐岩内油水识别的方法的流程图，由图I可知，所述的方法包括SlOl :采集地震三维工区的地震纵波传输信号。在具体的实施例中，诸如在塔里木盆地某地震三维工区通过布设信号激发和接收装置，采集地震纵波传输信号。S102 :对所述的地震纵波传输信号进行叠前处理。图2为本发明实施例提供的一种碳酸盐岩内油水识别的方法中步骤S102的具体流程图，由图2可知，步骤S102具体包括S201 :对所述的地震纵波传输信号进行常规校正，主要包括预处理、去噪、振幅补偿、反褶积和静校正等。S202 :对所述的地震纵波传输信号进行保幅叠前时间偏移处理，在具体的实施方式中，采用的主要方法是以经典的地震波绕射叠加理论为基础，引入保幅加权函数，应用Kirchhoff积分法实现保幅偏移。S103 :对叠前处理后的地震数据进行常规地震层位标定和解释，图3为该步骤的具体流程图，由图3可知，该步骤具体包括S301 :根据已知钻井资料以及所述的地震纵波传输信号制作合成记录；S302 :根据所述的合成记录，将钻井层位与地震层位进行对应；S303 :根据所述钻井层位与地震层位的对应关系，对地震层位进行解释。在具体的实施例中，将经过叠前时间偏移的地震数据导入到地震解释反演软件平台中(诸如Jason地震解释模块)，根据已知钻井资料诸如哈7-3井(Ha7_3)的声波时差制作合成记录，将测井层位和地震层位对应起来，并在地震解释模块中完成剖面上层位解释(如图9中全频段地震数据剖面)。在该实施例中解释了三个地震层位TS、T03t>和TCV2y 2，分别对应志留系柯坪塔格组底部界线、奥陶系土木休克组底部界线和奥陶系鹰山组第二段底部界线，各地层命名采用塔里木盆地通用名，其中土木休克组底部界线是不整合面，储层易发育在此界面附近。S104 :对叠前处理后的地震数据进行波阻抗反演。在具体的实施例中，诸如应用Jason软件地震反演模块，对叠前全频段地震数据进行波阻抗反演，波阻抗剖面根据波阻抗值大小进行颜色填充，低阻抗值对应于储层段(如图10中叠前反演波阻抗剖面，不同波阻抗值用不同颜色填充，从剖面上可以直观地观察到波阻抗低值区，即指示着储层发育区，在剖面上标注了储层的具体位置)，即灰色填充区，但储层内充填的流体类型尚无法确定。S105 :根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率。图4为该步骤的具体流程图，由图4可知，该步骤具体包括S401 :根据所述的波阻抗，对地震数据进行分角度部分的叠加得到近道集数据和远道集数据。该步骤可通过如下公式进行
nearNear = J^x(I)
HU I
权利要求
1.一种碳酸盐岩内油水识别的方法，其特征是，所述的方法包括 采集地震三维工区的地震纵波传输信号； 对所述的地震纵波传输信号进行叠前处理； 对叠前处理后的地震数据进行常规地震层位标定和解释； 对叠前处理后的地震数据进行波阻抗反演； 根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率； 对所述的频率对应的地震数据进行归一化处理； 根据所述的归一化处理的结果进行油水识别。
2.根据权利要求I所述的方法，其特征是，对所述的地震纵波传输信号进行叠前处理具体包括 对所述的地震纵波传输信号进行常规校正； 对所述的地震纵波传输信号进行保幅叠前时间偏移处理。
3.根据权利要求I所述的方法，其特征是，所述的对叠前处理后的地震数据进行常规地震层位标定和解释具体包括 根据已知钻井资料以及所述的地震纵波传输信号制作合成记录； 根据所述的合成记录，将钻井层位与地震层位进行对应； 根据所述钻井层位与地震层位的对应关系，对地震层位进行解释。
4.根据权利要求I所述的方法，其特征是，根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率具体包括 根据所述的波阻抗，对地震数据进行分角度的部分叠加得到近道集数据和远道集数据； 分别在所述的近道集数据以及远道集数据的目的层段内随机抽取多个地震道； 对所述的多个地震道进行频谱分析，确定频段参数； 对所述的近道集数据、远道集数据进行分频处理； 提取所述的地震数据中AVO效应明显的地震近道集数据以及远道集数据对应的频率。
5.根据权利要求4所述的方法，其特征是，根据所述的波阻抗，对地震数据进行分角度部分叠加得到近道集数据和远道集数据通过如下公式进行 nearNear = Yx(i)i=I(I) farFar = y xf/)-(2) 其中，Near为小角度范围内的道集叠加起来得到的近道集数据，Far为大角度范围内的道集叠加起来得到的远道集数据，η为大角度范围的起始角度，x(i)为地震道。
6.根据权利要求4所述的方法，其特征是，根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率还包括 将分频处理后的近道集数据、远道集数据进行平滑处理。
7.根据权利要求6所述的方法，其特征是，所述的将分频处理后的近道集数据、远道集数据进行平滑处理通过如下公式进行
8.根据权利要求7所述的方法，其特征是，对所述的频率对应的地震数据进行归一化处理具体包括 将所述的频率对应的近道集数据和远道集数据的振幅值归并到相同量纲级； 确定所述的近道集数据和远道集数据的振幅差； 根据所述的振幅差进行不同颜色的填充。
9.一种碳酸盐岩内油水识别的系统，其特征是，所述的系统包括 采集装置，用于采集地震三维工区的地震纵波传输信号； 叠前处理装置，用于对所述的地震纵波传输信号进行叠前处理； 层位标定装置，用于对叠前处理后的地震数据进行常规地震层位标定和解释； 波阻抗反演装置，用于对叠前处理后的地震数据进行波阻抗反演； 频率提取装置，用于根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率； 归一化处理装置，用于对所述的频率对应的地震数据进行归一化处理； 油水识别装置，用于根据所述的归一化处理的结果进行油水识别。
10.根据权利要求9所述的系统，其特征是，所述的叠前处理装置具体包括 校正单元，用于对所述的地震纵波传输信号进行常规校正； 保幅偏移单元，用于对所述的地震纵波传输信号进行保幅叠前时间偏移处理。
11.根据权利要求9所述的系统，其特征是，所述的层位标定装置具体包括 合成记录制作单元，用于根据已知钻井资料以及所述的地震纵波传输信号制作合成记录; 井震对应单元，用于根据所述的合成记录，将钻井层位与地震层位进行对应； 地震层位解释单元，用于根据所述钻井层位与地震层位的对应关系，对地震层位进行解释。
12.根据权利要求9所述的系统，其特征是，所述的频率提取装置具体包括 叠加单元，用于根据所述的波阻抗，对地震数据进行分角度的部分叠加得到近道集数据和远道集数据； 地震道提取单元，用于分别在所述的近道集数据以及远道集数据的目的层段内随机抽取多个地震道； 分频单元，用于对所述的多个地震道进行频谱分析，确定频段参数，对所述的近道集数据、远道集数据进行分频处理； 频率提取单元，用于提取所述的地震数据中AVO效应明显的近道集数据和远道集数据对应的频率。
13.根据权利要求12所述的系统，其特征是，所述的叠加单元通过如下公式进行
14.根据权利要求12所述的系统，其特征是，所述的频率提取装置还包括 平滑处理单元，用于将分频后的近道集数据、远道集数据进行平滑处理。
15.根据权利要求14所述的系统，其特征是，所述的平滑处理单元通过如下公式进行
16.根据权利要求15所述的系统，其特征是，所述的归一化处理装置具体包括 振幅差确定单元，用于将所述的频率对应的近道集数据和远道集数据的振幅值归并到相同量纲级，并确定所述的近道集数据和远道集数据的振幅差； 颜色填充单元，用于根据所述的振幅差进行不同颜色的填充。
全文摘要
本发明实施例提供了提供一种碳酸盐岩内油水识别的方法及系统，所述的方法包括采集地震三维工区的地震纵波传输信号；对所述的地震纵波传输信号进行叠前处理；对叠前处理后的地震数据进行常规地震层位标定和解释；对叠前处理后的地震数据进行波阻抗反演；根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率；对所述的频率对应的地震数据进行归一化处理；根据所述的归一化处理的结果进行油水识别。可以在平面和剖面上同时成图，更便于单个井点的油水识别和平面上地质规律的统计分析，充分发挥其钻前预测的优势。
文档编号G01V1/28GK102819038SQ201210265359
公开日2012年12月12日 申请日期2012年7月27日 优先权日2012年7月27日
发明者李艳东, 张研, 李明, 周波, 陶小晚, 唐俊伟, 李浩武, 马德波 申请人:中国石油天然气股份有限公司