CN102819038B - 一种碳酸盐岩内油水识别的方法及系统 - Google Patents

Abstract

本发明实施例提供了一种碳酸盐岩内油水识别的方法及系统，所述的方法包括：采集地震三维工区的地震纵波传输信号；对所述的地震纵波传输信号进行叠前处理；对叠前处理后的地震数据进行常规地震层位标定和解释；对叠前处理后的地震数据进行波阻抗反演；根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率；对所述的频率对应的地震数据进行归一化处理；根据所述的归一化处理的结果进行油水识别。可以在平面和剖面上同时成图，更便于单个井点的油水识别和平面上地质规律的统计分析，充分发挥其钻前预测的优势。

Description

一种碳酸盐岩内油水识别的方法及系统

技术领域

本发明关于地球物理勘探领域，特别是关于油气勘探开发领域，具体的讲是一种碳酸盐岩内油水识别的方法及系统。

背景技术

近年来，碳酸盐岩逐渐成为油气勘探非常重要的一个新领域。我国的碳酸盐岩油气田通常埋藏很深，储层成因复杂，主要包括同生期和准同生期岩溶作用、表生岩溶作用和埋藏岩溶作用。而不同地质背景下和不同成因形成的储层类型多种多样，碳酸盐岩储层通常具有复杂的孔、洞和缝结构，非均值性强。在勘探中准确预测碳酸盐岩储层类型已属不易，而辨析碳酸盐岩储层中复杂的油气水分布关系更为困难。因为，和均值砂岩储层相比，大多数碳酸盐岩油藏没有统一的油水界面，即便是相邻很近的储集空间内充注的流体也可能不同。因此，如何进行碳酸盐岩储层油水识别是本领域的一个技术难题。

目前，油气田勘探开发中常用的烃类检测技术是利用地震反射(或折射)特征进行油气藏的判识，包括地震属性技术、AVO（Amplitude Versus Offset,振幅随偏移距的变化）技术、合成声波测井技术、纵横波速度比技术等，其中最具代表性是AVO技术。但是传统的AVO分析只与介质分界面两边的弹性参数有关。另外，传统的AVO技术在实际应用中存在不便：在众多的三维地震剖面上依靠肉眼观察来寻找振幅增强或减弱现象并不完全可靠，容易造成遗漏或重复；在平面上无法成图，不便于分析地质现象在平面上的分布规律。

发明内容

为了克服现有技术存在的上述问题，本发明提出了一种碳酸盐岩内油水识别的方法及系统，通过在对地震资料做分频处理的基础上，求取近、远道地震振幅差值的方法来进行油水识别，既考虑了频率因素，也将地震剖面上定性的振幅变化现象转换为可以量化描述的振幅差值的信息，可以在平面和剖面上同时成图，更便于单个井点的油水识别和平面上地质规律的统计分析，充分发挥其钻前预测的优势。

本发明的目的之一是，提供一种碳酸盐岩内油水识别的方法，所述的方法包括：采集地震三维工区的地震纵波传输信号；对所述的地震纵波传输信号进行叠前处理；对叠前处理后的地震数据进行常规地震层位标定和解释；对叠前处理后的地震数据进行波阻抗反演；根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率；对所述的频率对应的地震数据进行归一化处理；根据所述的归一化处理的结果进行油水识别。

本发明的目的之一是，提供了一种碳酸盐岩内油水识别的系统，所述的系统包括：采集装置，用于采集地震三维工区的地震纵波传输信号；叠前处理装置，用于对所述的地震纵波传输信号进行叠前处理；层位标定装置，用于对叠前处理后的地震数据进行常规地震层位标定和解释；波阻抗反演装置，用于对叠前处理后的地震数据进行波阻抗反演；频率提取装置，用于根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率；归一化处理装置，用于对所述的频率对应的地震数据进行归一化处理；油水识别装置，用于根据所述的归一化处理的结果进行油水识别。

本发明的有益效果在于，通过在对地震资料做分频处理的基础上，求取近、远道地震振幅差值的方法来进行油水识别，既考虑了频率因素，也将地震剖面上定性的振幅变化现象转换为可以量化描述的振幅差值的信息，可以在平面和剖面上同时成图，更便于单个井点的油水识别和平面上地质规律的统计分析，充分发挥其钻前预测的优势，预测结果更准确。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例提供的一种碳酸盐岩内油水识别的方法的流程图；

图2为本发明实施例提供的一种碳酸盐岩内油水识别的方法中步骤S102的具体流程图；

图3为图1中的步骤S103的具体流程图；

图4为本发明实施例提供的一种碳酸盐岩内油水识别的方法中步骤S105的具体流程图；

图5为本发明实施例提供的一种碳酸盐岩内油水识别的系统的实施方式一的结构框图；

图6为本发明实施例提供的一种碳酸盐岩内油水识别的系统的实施方式二的结构框图；

图7为本发明实施例提供的一种碳酸盐岩内油水识别的系统的实施方式三的结构框图；

图8为本发明实施例提供的一种碳酸盐岩内油水识别的系统的实施方式四的结构框图；

图9为本发明实施例提供的一种碳酸盐岩内油水识别方法中全频段地震数据剖面；

图10为本发明实施例提供的一种碳酸盐岩内油水识别方法中叠前反演波阻抗剖面；

图11中，（a）为10Hz的近道集的振幅剖面；

（b）为10Hz的远道集的振幅剖面；

（c）是40Hz的近道集的振幅剖面；

（d）是40Hz的远道集分量的振幅剖面；

图12为40Hz的近、远道集的振幅差值剖面图；

图13为40Hz的近、远道集的振幅差值平面图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

图1为本发明实施例提供的一种碳酸盐岩内油水识别的方法的流程图，由图1可知，所述的方法包括：

S101：采集地震三维工区的地震纵波传输信号。在具体的实施例中，诸如在塔里木盆地某地震三维工区通过布设信号激发和接收装置，采集地震纵波传输信号。

S102：对所述的地震纵波传输信号进行叠前处理。图2为本发明实施例提供的一种碳酸盐岩内油水识别的方法中步骤S102的具体流程图，由图2可知，步骤S102具体包括：

S201：对所述的地震纵波传输信号进行常规校正，主要包括：预处理、去噪、振幅补偿、反褶积和静校正等。

S202：对所述的地震纵波传输信号进行保幅叠前时间偏移处理，在具体的实施方式中，采用的主要方法是以经典的地震波绕射叠加理论为基础，引入保幅加权函数，应用Kirchhoff积分法实现保幅偏移。

S103：对叠前处理后的地震数据进行常规地震层位标定和解释，图3为该步骤的具体流程图，由图3可知，该步骤具体包括：

S301：根据已知钻井资料以及所述的地震纵波传输信号制作合成记录；

S302：根据所述的合成记录，将钻井层位与地震层位进行对应；

S303：根据所述钻井层位与地震层位的对应关系，对地震层位进行解释。

在具体的实施例中，将经过叠前时间偏移的地震数据导入到地震解释反演软件平台中（诸如Jason地震解释模块），根据已知钻井资料诸如哈7-3井（Ha7-3）的声波时差制作合成记录，将测井层位和地震层位对应起来，并在地震解释模块中完成剖面上层位解释（如图9中全频段地震数据剖面）。在该实施例中解释了三个地震层位：TS、TO_{3t_zkz}和TO_{1-2y_2}，分别对应志留系柯坪塔格组底部界线、奥陶系土木休克组底部界线和奥陶系鹰山组第二段底部界线，各地层命名采用塔里木盆地通用名，其中土木休克组底部界线是不整合面，储层易发育在此界面附近。

S104：对叠前处理后的地震数据进行波阻抗反演。在具体的实施例中，诸如应用Jason软件地震反演模块，对叠前全频段地震数据进行波阻抗反演，波阻抗剖面根据波阻抗值大小进行颜色填充，低阻抗值对应于储层段（如图10中叠前反演波阻抗剖面，不同波阻抗值用不同颜色填充，从剖面上可以直观地观察到波阻抗低值区，即指示着储层发育区，在剖面上标注了储层的具体位置），即灰色填充区，但储层内充填的流体类型尚无法确定。

S105：根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率。图4为该步骤的具体流程图，由图4可知，该步骤具体包括：

S401：根据所述的波阻抗，对地震数据进行分角度部分的叠加得到近道集数据和远道集数据。该步骤可通过如下公式进行：

$\mathrm{Near}=\underset{i=1}{\overset{\mathrm{near}}{\mathrm{\Σ}}}x\left(i\right)---\left(1\right)$

$\mathrm{Far}=\underset{i=n}{\overset{\mathrm{far}}{\mathrm{\Σ}}}x\left(i\right)---\left(2\right)$

其中，Near为小角度范围内的道集叠加起来得到的近道集数据，Far为大角度范围内的道集叠加起来得到的远道集数据，n为大角度范围的起始角度，x(i)为地震道。诸如在具体的实施例中，对角度范围在[5°～15°]度之间的地震波数据进行叠加得到近道叠加数据，对[25°～35°]范围内的数据进行叠加得到远道叠加数据。

S402：分别在所述的近道集数据以及远道集数据的目的层段内随机抽取多个地震道；

S403：对所述的多个地震道进行频谱分析，确定频段参数；

S404：对所述的近道集数据、远道集数据进行分频处理，该步骤具体包括：

根据基于Wigner-Ville分布的频谱分解算法，对所述的多个地震道进行分频。在具体的实施例中，诸如采用基于Wigner-Ville分布的频谱分解算法，在基于MATLAB平台二次开发的模块PowerReservoir中来实现的。

S405：提取所述的地震数据中AVO效应明显的地震近道集数据以及远道集数据对应的频率。例如，在具体的实施例中分别抽取了10Hz的近、远道集的叠加数据体，40Hz的近、远道集的叠加数据体，从剖面上观察其振幅谱，如图11中，（a）为10Hz的近道集的振幅剖面，（b）为10Hz的远道集的振幅剖面，（c）是40Hz的近道集的振幅剖面，（d）是40Hz的远道集分量的振幅剖面，在实际操作中可以抽取地震频带范围内任意频率值的地震道集数据体。

步骤S105除了上述S401-S405之外，还包括：将分频后的多个地震到进行平滑处理。在具体的实施例中，诸如当对WV分布进行了平滑处理，平滑WV分布通过如下公式进行：

$\mathrm{SPWD}(t,f)=\underset{-\∞}{\overset{+\∞}{\∫}}\underset{-\∞}{\overset{+\∞}{\∫}}\mathrm{\φ}(u,v){\mathrm{PWD}}_x(t-u,f-v)\mathrm{dudv}---\left(3\right)$

其中，其中，φ(u，v)为二维低通滤波器，u为时间积分变量，v为频率积分变量，t为时间，f为频率。

S106：对所述的频率对应的地震数据进行归一化处理。该步骤具体包括：

（1）、将所述的频率对应的近道集数据和远道集数据的振幅值归并到相同量纲级；

（2）、确定所述的近道集数据和远道集数据的振幅差；

（3）、根据所述的振幅差进行不同颜色的填充。

在具体的实施方式中，优选AVO效应最明显的地震近、远道集所对应的频率。诸如在本实施例中，从10Hz近、远道集的叠加数据体（图11（a）、（b））振幅剖面上看，振幅差别不明显，即在该频率下的地震AVO效应不明显，不足以进行流体类型的判别；但从40Hz近、远道集的叠加数据体（图11（c）、（d））振幅剖面上看，振幅差别较明显，即在该频率下的充填不同流体介质的地震AVO效应显现出来。根据油、水特定的AVO效应表现形式，即：油层表现为远道振幅大于近道振幅，水层表现为远道振幅小于近道振幅，据此规律可区别流体类型。由此，优选40Hz的近、远道集进行振幅差值计算。

S107：根据所述的归一化处理的结果进行油水识别。求取40Hz地震近、远道集振幅差，振幅差的计算在具体的实施方式中可采用基于MATLAB平台二次开发的模块PowerReservoir来实现，根据差值大小进行颜色填充，分别在剖面（诸如图12为40Hz的近、远道集的振幅差值剖面图，用不同颜色表示振幅差值大小，黑色表示远道振幅大于近道，为油层，灰色表示远道振幅小于近道，为底水层）和平面上（诸如图13为40Hz的近、远道集的振幅差值平面图，用不同颜色表示振幅差值大小，黑色区域表示油层平面分布范围，灰色区域表示水层分布范围）显示，黑色表示油区分布，灰色表示水区分布，据此进行区域的油水识别。

图5为本发明实施例提供的一种碳酸盐岩内油水识别的系统的实施方式一的结构框图，由图5可知，所述的系统包括：

采集装置100，用于采集地震三维工区的地震纵波传输信号。在具体的实施例中，诸如在塔里木盆地某地震三维工区通过布设信号激发和接收装置，采集地震纵波传输信号。

叠前处理装置200，用于对所述的地震纵波传输信号进行叠前处理。图6为本发明实施例提供的一种碳酸盐岩内油水识别的系统的实施方式二的结构框图，由图6可知，叠前处理装置200具体包括：

校正单元201，用于对所述的地震纵波传输信号进行常规校正，主要包括：预处理、去噪、振幅补偿、反褶积和静校正等。

保幅偏移单元202，用于对所述的地震纵波传输信号进行保幅叠前时间偏移处理，在具体的实施方式中，采用的主要方法是以经典的地震波绕射叠加理论为基础，引入保幅加权函数，应用Kirchhoff积分法实现保幅偏移。

所述的系统还包括：

层位标定装置300，用于对叠前处理后的地震数据进行常规地震层位标定和解释。图7为本发明实施例提供的一种碳酸盐岩内油水识别的系统的实施方式三的结构框图，由图7可知，所述的层位标定装置300具体包括：

合成记录制作单元301，用于根据已知钻井资料以及所述的地震纵波传输信号制作合成记录；

井震对应单元302，用于根据所述的合成记录，将钻井层位与地震层位进行对应；

地震层位解释单元303，用于根据所述钻井层位与地震层位的对应关系，对地震层位进行解释。

在具体的实施例中，将经过叠前时间偏移的地震数据导入到地震解释反演软件平台中（诸如Jason地震解释模块），根据已知钻井资料诸如哈7-3井（Ha7-3）的声波时差制作合成记录，将测井层位和地震层位对应起来，并在地震解释模块中完成剖面上层位解释（如图9中全频段地震数据剖面）。在该实施例中解释了三个地震层位：TS、TO_{3t_zkz}和TO_{1-2y_2}，分别对应志留系柯坪塔格组底部界线、奥陶系土木休克组底部界线和奥陶系鹰山组第二段底部界线，各地层命名采用塔里木盆地通用名，其中土木休克组底部界线是不整合面，储层易发育在此界面附近。

所述的系统还包括：波阻抗反演装置400，用于对叠前处理后的地震数据进行波阻抗反演。在具体的实施例中，诸如应用Jason软件地震反演模块，对叠前全频段地震数据进行波阻抗反演，波阻抗剖面根据波阻抗值大小进行颜色填充，低阻抗值对应于储层段（如图10中叠前反演波阻抗剖面，不同波阻抗值用不同颜色填充，从剖面上可以直观地观察到波阻抗低值区，即指示着储层发育区，在剖面上标注了储层的具体位置），即灰色填充区，但储层内充填的流体类型尚无法确定。

频率提取装置500，用于根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率。图8为本发明实施例提供的一种碳酸盐岩内油水识别的系统的实施方式四的结构框图，由图8可知，所述的频率提取装置400具体包括：

叠加单元501，用于根据所述的波阻抗，对地震数据进行分角度部分的叠加得到近道集数据和远道集数据。该叠加单元可通过如下公式进行：

$\mathrm{Near}=\underset{i=1}{\overset{\mathrm{near}}{\mathrm{\Σ}}}x\left(i\right)---\left(1\right)$

$\mathrm{Far}=\underset{i=n}{\overset{\mathrm{far}}{\mathrm{\Σ}}}x\left(i\right)---\left(2\right)$

其中，Near为小角度范围内的道集叠加起来得到的近道集数据，Far为大角度范围内的道集叠加起来得到的远道集数据，n为大角度范围的起始角度，x(i)为地震道。诸如在具体的实施例中，对角度范围在[5°～15°]度之间的地震波数据进行叠加得到近道叠加数据，对[25°～35°]范围内的数据进行叠加得到远道叠加数据。

地震道提取单元502，用于分别在所述的近道集数据以及远道集数据的目的层段内随机抽取多个地震道；

分频单元503，用于对所述的多个地震道进行频谱分析，确定频段参数，对所述的近道集数据、远道集数据进行分频处理.

根据基于Wigner-Ville分布的频谱分解算法，对所述的多个地震道进行分频，在具体的实施例中，诸如采用基于Wigner-Ville分布的频谱分解算法，在基于MATLAB平台二次开发的模块PowerReservoir中来实现的。

频率提取单元504，用于提取所述的地震数据中AVO效应明显的近道集数据和远道集数据对应的频率。

例如，在具体的实施例中分别抽取了10Hz的近、远道集的叠加数据体，40Hz的近、远道集的叠加数据体，从剖面上观察其振幅谱，如图11中，（a）为10Hz的近道集的振幅剖面，（b）为10Hz的远道集的振幅剖面，（c）是40Hz的近道集的振幅剖面，（d）是40Hz的远道集的振幅剖面，在实际操作中可以抽取地震频带范围内任意频率值的地震道集数据体。

频率提取装置500还包括：平滑处理单元，用于将分频后的多个地震到进行平滑处理。在具体的实施例中，诸如当对WV分布进行了平滑处理，平滑WV分布通过如下公式进行：

$\mathrm{SPWD}(t,f)=\underset{-\∞}{\overset{+\∞}{\∫}}\underset{-\∞}{\overset{+\∞}{\∫}}\mathrm{\φ}(u,v){\mathrm{PWD}}_x(t-u,f-v)\mathrm{dudv}---\left(3\right)$

其中，其中，φ(u，v)为二维低通滤波器，u为时间积分变量，v为频率积分变量，t为时间，f为频率。

所述的系统还包括：归一化处理装置600，用于对所述的频率对应的地震数据进行归一化处理，主要包括：

振幅差确定单元，用于将所述的频率对应的近道集数据和远道集数据的振幅值归并到相同量纲，并确定所述的近道集数据和远道集数据的振幅差；

颜色填充单元，用于根据所述的振幅差进行不同颜色的填充。

在具体的实施方式中，优选AVO效应最明显的地震近、远道集所对应的频率。诸如在本实施例中，从10Hz近、远道集的叠加数据体（图11（a）、（b））振幅剖面上看，振幅差别不明显，即在该频率下的地震AVO效应不明显，不足以进行流体类型的判别；但从40Hz近、远道集的叠加数据体（图11（c）、（d））振幅剖面上看，振幅差别较明显，即在该频率下的充填不同流体介质的地震AVO效应显现出来。根据油、水特定的AVO效应表现形式，即油层表现为远道振幅大于近道振幅，水层表现为远道振幅小于近道振幅，据此规律可区别流体类型。由此，优选40Hz的近、远道集进行振幅差值计算。

所述的系统还包括识别装置700，用于根据所述的归一化处理的结果进行油水识别。

求取40Hz地震近、远道集振幅差，振幅差的计算在具体的实施方式中可采用基于MATLAB平台二次开发的模块PowerReservoir来实现，根据差值大小进行颜色填充，分别在剖面（诸如图12为40Hz的近、远道集的振幅差值剖面图，用不同颜色表示振幅差值大小，黑色表示远道振幅大于近道，为油层，灰色表示远道振幅小于近道，为底水层）和平面上（诸如图13为40Hz的近、远道集的振幅差值平面图，用不同颜色表示振幅差值大小，黑色区域表示油层平面分布范围，灰色区域表示水层分布范围）显示，黑色表示油区分布，灰色表示水区分布，据此进行区域的油水识别。

综上所述，本发明的有益成果是：提出了一种碳酸盐岩内油水识别的方法及系统，通过在对地震资料做分频处理的基础上，求取近、远道地震振幅差值的方法来进行油水识别，既考虑了频率因素，也将地震剖面上定性的振幅变化现象转换为可以量化描述的振幅差值的信息，可以在平面和剖面上同时成图，更便于单个井点的油水识别和平面上地质规律的统计分析，充分发挥其钻前预测的优势，预测结果更准确。

本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (16)

1.一种碳酸盐岩内油水识别的方法，其特征是，所述的方法包括：

采集地震三维工区的地震纵波传输信号；

对所述的地震纵波传输信号进行叠前处理；

对叠前处理后的地震数据进行常规地震层位标定和解释；

对叠前处理后的地震数据进行波阻抗反演；

根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率；

对所述的频率对应的地震数据进行归一化处理；

根据所述的归一化处理的结果进行油水识别。

2.根据权利要求1所述的方法，其特征是，对所述的地震纵波传输信号进行叠前处理具体包括：

对所述的地震纵波传输信号进行常规校正；

对所述的地震纵波传输信号进行保幅叠前时间偏移处理。

3.根据权利要求1所述的方法，其特征是，所述的对叠前处理后的地震数据进行常规地震层位标定和解释具体包括：

根据已知钻井资料以及所述的地震纵波传输信号制作合成记录；

根据所述的合成记录，将钻井层位与地震层位进行对应；

根据所述钻井层位与地震层位的对应关系，对地震层位进行解释。

4.根据权利要求1所述的方法，其特征是，根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率具体包括：

根据所述的波阻抗，对地震数据进行分角度的部分叠加得到近道集数据和远道集数据；

分别在所述的近道集数据以及远道集数据的目的层段内随机抽取多个地震道；

对所述的多个地震道进行频谱分析，确定频段参数；

对所述的近道集数据、远道集数据进行分频处理；

提取所述的地震数据中AVO效应明显的地震近道集数据以及远道集数据对应的频率。

5.根据权利要求4所述的方法，其特征是，根据所述的波阻抗，对地震数据进行分角度部分叠加得到近道集数据和远道集数据通过如下公式进行：

$\mathrm{Near}=\underset{i=1}{\overset{\mathrm{near}}{\mathrm{\Σ}}}x\left(i\right)---\left(1\right)$

$\mathrm{Far}=\underset{i=n}{\overset{\mathrm{far}}{\mathrm{\Σ}}}x\left(i\right)---\left(2\right)$

其中，Near为小角度范围内的道集叠加起来得到的近道集数据，Far为大角度范围内的道集叠加起来得到的远道集数据，n为大角度范围的起始角度，x(i)为地震道。

6.根据权利要求4所述的方法，其特征是，根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率还包括：

将分频处理后的近道集数据、远道集数据进行平滑处理。

7.根据权利要求6所述的方法，其特征是，所述的将分频处理后的近道集数据、远道集数据进行平滑处理通过如下公式进行：

$\mathrm{SPWD}(t,f)=\underset{-\∞}{\overset{+\∞}{\∫}}\underset{-\∞}{\overset{+\∞}{\∫}}\mathrm{\φ}(u,v){\mathrm{PWD}}_x(t-u,f-v)\mathrm{dudv}---\left(3\right)$

其中，φ(u,v)为二维低通滤波器，u为时间积分变量，v为频率积分变量，t为时间，f为频率。

8.根据权利要求7所述的方法，其特征是，对所述的频率对应的地震数据进行归一化处理具体包括：

将所述的频率对应的近道集数据和远道集数据的振幅值归并到相同量纲级；

确定所述的近道集数据和远道集数据的振幅差；

根据所述的振幅差进行不同颜色的填充。

9.一种碳酸盐岩内油水识别的系统，其特征是，所述的系统包括：

采集装置，用于采集地震三维工区的地震纵波传输信号；

叠前处理装置，用于对所述的地震纵波传输信号进行叠前处理；

层位标定装置，用于对叠前处理后的地震数据进行常规地震层位标定和解释；

波阻抗反演装置，用于对叠前处理后的地震数据进行波阻抗反演；

频率提取装置，用于根据所述的波阻抗以及所述的地震数据的振幅随偏移距的变化AVO效应提取地震数据对应的频率；

归一化处理装置，用于对所述的频率对应的地震数据进行归一化处理；

油水识别装置，用于根据所述的归一化处理的结果进行油水识别。

10.根据权利要求9所述的系统，其特征是，所述的叠前处理装置具体包括：

校正单元，用于对所述的地震纵波传输信号进行常规校正；

保幅偏移单元，用于对所述的地震纵波传输信号进行保幅叠前时间偏移处理。

11.根据权利要求9所述的系统，其特征是，所述的层位标定装置具体包括：

合成记录制作单元，用于根据已知钻井资料以及所述的地震纵波传输信号制作合成记录；

井震对应单元，用于根据所述的合成记录，将钻井层位与地震层位进行对应；

地震层位解释单元，用于根据所述钻井层位与地震层位的对应关系，对地震层位进行解释。

12.根据权利要求9所述的系统，其特征是，所述的频率提取装置具体包括：

叠加单元，用于根据所述的波阻抗，对地震数据进行分角度的部分叠加得到近道集数据和远道集数据；

地震道提取单元，用于分别在所述的近道集数据以及远道集数据的目的层段内随机抽取多个地震道；

分频单元，用于对所述的多个地震道进行频谱分析，确定频段参数，对所述的近道集数据、远道集数据进行分频处理；

频率提取单元，用于提取所述的地震数据中AVO效应明显的近道集数据和远道集数据对应的频率。

13.根据权利要求12所述的系统，其特征是，所述的叠加单元通过如下公式进行：

$\mathrm{Near}=\underset{i=1}{\overset{\mathrm{near}}{\mathrm{\Σ}}}x\left(i\right)---\left(1\right)$

$\mathrm{Far}=\underset{i=n}{\overset{\mathrm{far}}{\mathrm{\Σ}}}x\left(i\right)---\left(2\right)$

其中，Near为小角度范围内的道集叠加起来得到的近道集数据，Far为大角度范围内的道集叠加起来得到的远道集数据，n为大角度范围的起始角度，x(i)为地震道。

14.根据权利要求12所述的系统，其特征是，所述的频率提取装置还包括：

平滑处理单元，用于将分频后的近道集数据、远道集数据进行平滑处理。

15.根据权利要求14所述的系统，其特征是，所述的平滑处理单元通过如下公式进行：

$\mathrm{SPWD}(t,f)=\underset{-\∞}{\overset{+\∞}{\∫}}\underset{-\∞}{\overset{+\∞}{\∫}}\mathrm{\φ}(u,v){\mathrm{PWD}}_x(t-u,f-v)\mathrm{dudv}---\left(3\right)$

其中，φ(u，v)为二维低通滤波器，u为时间积分变量，v为频率积分变量，t为时间，f为频率。

16.根据权利要求15所述的系统，其特征是，所述的归一化处理装置具体包括：

振幅差确定单元，用于将所述的频率对应的近道集数据和远道集数据的振幅值归并到相同量纲级，并确定所述的近道集数据和远道集数据的振幅差；

颜色填充单元，用于根据所述的振幅差进行不同颜色的填充。

Images