CN112130209A - 岩溶储层预测方法及装置 - Google Patents

Abstract

本发明提供了一种岩溶储层预测方法及装置，方法包括：对目的层地震数据进行多子波分解、分频扫描重构确定分频段子波扫描重构剖面；判断其是否符合实际钻井情况；符合，则利用其中反应岩溶储层的频带数据进行岩溶储层预测；不符合，则判断分频段子波扫描重构剖面是否反应岩溶响应特征；确定反应岩溶响应特征，则对岩溶响应频段的子波数据进行数据体重构，进行岩溶储层预测；否则，则进行波形分解重构确定岩溶响应最佳分量重构数据体，进行岩溶储层预测。本发明考虑了地震子波向下传播过程中形状、频率和能量的变化，能够有效克服常规单一子波的局限性，实现了准确识别岩溶储层和有效预测其含油气性的目的,并与实际钻井情况相吻合。

Description

岩溶储层预测方法及装置

技术领域

本发明涉及地球物理勘探技术，特别是关于油气勘探技术，具体的讲是一种岩溶储层预测方法及装置。

背景技术

由于地下岩层岩性差异导致的较大波阻抗差在目的层段形成强波峰反射，降低地震的纵向分辨率，掩盖了岩溶储层的地震响应特征，干扰了储层的有效预测。

现有技术中的常规储层预测的技术方法效果不明显，此外强反射层与岩溶储层地震响应存在一定的重合，降低了常规含油气检测方法可靠性。并且由于龙潭组的页岩和茅口组顶部灰岩的岩性差异形成的强波峰反射屏蔽了下覆岩溶的地震响应特征，受制于单一子波地震勘探的局限性，实际的地震记录与地下真实的地震响应存在一定的差异,导致实际生产中储层预测精度低。

发明内容

为提高岩溶储层预测的可靠性，对岩溶储层进行预测，本发明实施例提供了一种岩溶储层预测方法，方法包括：

对目的层地震数据进行多子波分解确定分频段子波数据；

对所述分频段子波数据进行分频扫描重构确定分频段子波扫描重构剖面；

根据实际钻井数据判断所述分频段子波扫描重构剖面是否符合实际钻井情况；

确定所述分频段子波扫描重构剖面符合实际钻井情况，则利用所述分频段子波数据中反应岩溶储层的频带数据进行岩溶储层预测；

判断所述分频段子波扫描重构剖面不符合实际钻井情况，则判断所述分频段子波扫描重构剖面是否反应岩溶响应特征；

确定反应岩溶响应特征，则对岩溶响应频段的子波数据进行数据体重构，利用重构后的岩溶响应频段的子波数据进行岩溶储层预测；

确定不反应岩溶响应特征，则对分频扫描重构后的子波数据进行波形分解重构确定岩溶响应最佳分量重构数据体，利用确定的岩溶响应最佳分量重构数据体进行岩溶储层预测。

本发明实施例中，所述的对目的层地震数据进行多子波分解确定分频段子波数据包括：

对目的层地震数据进行雷克子波、Paul子波或Morlet子波分解确定分频段子波数据。

本发明实施例中，所述的根据实际钻井数据判断所述分频段子波扫描重构剖面是否符合实际钻井情况包括：

根据目的层工区的实际钻井数据建立地质模型；

利用建立的地质模型进行正演模拟确定地震响应特征；

根据所述地震响应特征判断分频段子波扫描重构剖面是否符合实际钻井情况。

本发明实施例中，所述的对目的层地震数据进行多子波分解确定分频段子波数据包括：

对目的层地震数据进行频谱分析确定多子波分解的频带范围；

根据确定的频带范围进行多子波分解确定分频段子波数据。

同时，本发明还提供一种岩溶储层预测装置，包括：

子波分解模块，用于对目的层地震数据进行多子波分解确定分频段子波数据；

分频重构模块，用于对所述分频段子波数据进行分频扫描重构确定分频段子波扫描重构剖面；

分频判断模块，用于根据实际钻井数据判断所述分频段子波扫描重构剖面是否符合实际钻井情况；

分频预测模块，用于确定所述分频段子波扫描重构剖面符合实际钻井情况，则利用所述分频段子波数据中反应岩溶储层的频带数据进行岩溶储层预测；

岩溶响应特征判断模块，用于判断所述分频段子波扫描重构剖面不符合实际钻井情况，则判断所述分频段子波扫描重构剖面是否反应岩溶响应特征；

岩溶响应子波重构模块，用于确定反应岩溶响应特征，则对岩溶响应频段的子波数据进行数据体重构，利用重构后的岩溶响应频段的子波数据进行岩溶储层预测；

波形分解重构模块，用于确定不反应岩溶响应特征，则对分频扫描重构后的子波数据进行波形分解重构确定岩溶响应最佳分量重构数据体，利用确定的溶响应最佳分量重构数据体进行岩溶储层预测。

本发明实施例中，所述的子波分解模块对目的层地震数据进行多子波分解确定分频段子波数据包括：

对目的层地震数据进行雷克子波、Paul子波或Morlet子波分解确定分频段子波数据。

本发明实施例中，所述的分频判断模块包括：

建模单元，用于根据目的层工区的实际钻井数据建立地质模型；

响应特征确定单元，用于利用建立的地质模型进行正演模拟确定地震响应特征；

响应判断单元，用于根据所述地震响应特征判断分频段子波扫描重构剖面是否符合实际钻井情况。

本发明实施例中，所述的子波分解模块包括：

频谱分析单元，用于对目的层地震数据进行频谱分析确定多子波分解的频带范围；

分解单元，用于根据确定的频带范围进行多子波分解确定分频段子波数据。

同时，本发明还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述方法。

同时，本发明还提供一种计算机可读存储介质，计算机可读存储介质存储有执行上述方法的计算机程序。

本发明利用多子波地震道分解与重构技术，本发明考虑了地震子波向下传播过程中形状、频率和能量的变化，能够有效克服常规单一子波的局限性，能够推断出地震波在传播过程中经过不同地质体的地震响应特征，将本发明的方案用于某地区三维地震储层预测与含油气性检测。采用该技术处理后的地震数据在排除了强反射层的干扰,纵向分辨率得到大大提高,实现了准确识别岩溶储层和有效预测其含油气性的目的,并与实际钻井情况相吻合。

为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。

附图说明

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

图1为本发明公开的岩溶储层预测方法的流程图；

图2为本发明实施例中公开的进行岩溶储层预测方法的流程图；

图3为本发明实施例中建立的地质模型；

图4为本发明实施例中建立的地质模型；

图5为本发明实施例中建立的地质模型；

图6为本实施例中茅口组岩溶储层模型正演剖面；

图7为本实施例中茅口组岩溶储层模型正演剖面；

图8为本实施例中茅口组岩溶储层模型正演剖面；

图9本发明实施例中原始地震剖面；

图10为本发明实施例中30-40hz子波重构后剖面；

图11为本发明实施例中30-40hz子波重构+波形分解第3、4分量剖面；

图12为本发明实施例中目的层相同时窗的原始地震数据根振幅；

图13为本发明实施例中子波重构后的均方根振幅；

图14为本发明公开的岩溶储层预测装置的框图；

图15为本发明实施例公开的电子设备的示意图。

具体实施方式

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

如图1所示，为本发明公开的岩溶储层预测方法的流程图，其包括：

步骤S101，对目的层地震数据进行多子波分解确定分频段子波数据；

步骤S102，对所述分频段子波数据进行分频扫描重构确定分频段子波扫描重构剖面；

步骤S103，根据实际钻井数据判断所述分频段子波扫描重构剖面是否符合实际钻井情况；

步骤S104，确定所述分频段子波扫描重构剖面符合实际钻井情况，则利用所述分频段子波数据中反应岩溶储层的频带数据进行岩溶储层预测；

步骤S105，判断所述分频段子波扫描重构剖面不符合实际钻井情况，则判断所述分频段子波扫描重构剖面是否反应岩溶响应特征；

步骤S106，确定反应岩溶响应特征，则对岩溶响应频段的子波数据进行数据体重构，利用重构后的岩溶响应频段的子波数据进行岩溶储层预测；

步骤S107，确定不反应岩溶响应特征，则对分频扫描重构后的子波数据进行波形分解重构确定岩溶响应最佳分量重构数据体，利用确定的岩溶响应最佳分量重构数据体进行岩溶储层预测。

现有技术中，由于龙潭组的页岩和茅口组顶部灰岩的岩性差异形成的强波峰反射屏蔽了下覆岩溶的地震响应特征，受制于单一子波地震勘探的局限性，实际的地震记录与地下真实的地震响应存在一定的差异,导致实际生产中储层预测精度低。本发明利用多子波地震道分解与重构技术,本发明考虑了地震子波向下传播过程中形状、频率和能量的变化，能够有效克服常规单一子波的局限性，能够推断出地震波在传播过程中经过不同地质体的地震响应特征，将该方法用于某地区三维地震储层预测与含油气性检测。采用频率域子波重构+波形分解技术处理后的地震数据排除了强反射层的干扰,纵向分辨率得到大大提高,实现了准确识别岩溶储层和有效预测其含油气性的目的,并与实际钻井情况相吻合。

本发明实施例提出的一种利用多子波分解与重构，综合地震、测井、钻井识别岩溶储层的方法，具体包括以下内容：

1、利用井约束模型正演分析，得出不同地质条件下不同岩溶的的地震响应特征；

2、应用地质层位约束的沿层分频子波重构+波形分解技术，能有效识别岩溶地震特征，不同岩溶具有不同的子波频段，高频段子波重构可以识别强屏蔽下的薄储层；

3、油气具有高频衰减特征，基于子波重构技术可以检测含油气特征，其可靠性较高。

本实施例基于多子波理论，融合地震、钻井、测井、地质信息，通过不同频带子波重构筛选出岩溶储层的响应频段，在该基础上再利用波形重构，能够有效去除龙潭组的页岩和茅口组顶部灰岩的波组抗差形成的强波峰反射，恢复强屏蔽层下岩溶的地震的响应特征。该方法对落实强反射层屏蔽下的岩溶储层及其分布范围起到重要的作用，具有非常广阔的应用前景。实践证明，该方法流程简单、实用性强，成本低，预测精度高等优点，适用于精度更高的岩性圈闭勘探。

现结合本发明的示例性实施例对本发明技术方案做详细说明，所述实施例的示例在图2中示出。以下将通过参照图2的流程步骤来说明本发明的技术方案，以便解释本发明的技术方案。图2所示步骤具体包括：

目的层频谱分析；

多子波分解；

分解后数据分频段扫描；

目的层井点处子波扫描分析；

判断岩溶响应与强屏蔽层响应频段是否重合；

确定不重合，则执行岩溶响应频段的子波重构；进行重构数据体属性分析；重构数据体含油气检测；确定岩溶储层的空间范围。

确定重合，则执行重构数据体波形分解；找出岩溶响应最佳分量；然后进行重构数据体属性分析；重构数据体含油气检测；确定岩溶储层的空间范围。

其中，图2中并未表示判断分频段子波扫描重构剖面是否符合实际钻井情况，符合则利用所述分频段子波数据中反应岩溶储层的频带数据进行岩溶储层预测，不符合则进行图2中的进行目的层井点处子波扫描分析及后续步骤。

本发明利用多子波地震道分解与重构技术,本发明考虑了地震子波向下传播过程中形状、频率和能量的变化，能够有效克服常规单一子波的局限性，能够推断出地震波在传播过程中经过不同地质体的地震响应特征，将该方法用于某地区三维地震储层预测与含油气性检测。对于分频段子波扫描重构剖面符合实际钻井情况则利用所述分频段子波数据中反应岩溶储层的频带数据进行岩溶储层预测，对于不符合分频段子波扫描重构剖面不符合实际钻井情况，则进行波形分解重构，即采用频率域子波重构+波形分解技术处理后的地震数据排除了强反射层的干扰,纵向分辨率得到大大提高,实现了准确识别岩溶储层和有效预测其含油气性的目的,并与实际钻井情况相吻合。

本实施例的实施步骤如下：

(1)目的层频谱分析，从叠前时间偏移剖面分别抽取线号、道号间隔5km的剖面进行频谱分析，确定多子波数据分解的频带范围。

由于工区激发接收条件和地下地质情况的不同，不同区域岩溶储层的频带宽度不一样，本实施例中通过频谱分析确定多子波数据分解的频带范围，一方面确保进行多子波分解时不漏掉有效信号，另一方面确定了频带范围能够提高分解的效率，避免过宽的频带分解节约生产成本，并且为保证分解数据包含地震数据体所有有效信息，在频谱分析后的频带可以适当拓宽频带分解。

(2)多子波分解；本实施例中选取三种子波实验，分解结果与实际的地震数据对比，选取雷克子波作为分解的子波。为了提高分解效率，设置分解时窗和数据分解范围。沿层分频扫描重构，根据先宽后窄的思想，扫描增量由大到小，精细确定有效反射的频带宽度。

本实施例中，通过三种子波试验能够筛选出与实际的地震子波最匹配的一系列子波序列分解，这样多子波分解的结果与原始地震数据匹配度更高，后期数据重构、波形分解的结果更可靠，应用重构后的数据体来识别岩溶储层和含油气预测也更可靠。

本实施例中，子波实验中选取了Ricker(雷克)子波、Paul子波、Morlet子波，对本领域技术人员可知，任何能够筛选出与实际的地震子波最匹配的子波均涵盖在本申请的方案中，并不以本申请公开的内容为限。

具体的选用哪种子波试验根据情况而定，本发明实施例中，Ricker就能够满足生产的需求，而在特殊的的地质条件下，地震数据体的子波与雷克子波有较大的区别，需要采用其他的子波序列更能够反映用于分解的地震数据内部子波的变化。

本发明实施例中采用的是软件提供的三种子波，也可以用户自定义子波。通常情况下选取雷克子波分解，因为波波形简单、收敛速度迅速、与实际的地震信号相匹配，特殊情况下，当Ricker序列与数据体内地震子波区别较大时，分别将实验的不同子波序列分解的数据体进行数据重构，重构的数据体与原始地震数据对比，结果相差最小的就选择为分解的子波。

本发明提供的方法中，将实验的不同子波序列分解的数据体进行数据重构得到分频段子波扫描重构剖面，本发明提供的方法中，对于确定的分频段子波扫描重构剖面，判断分频段子波扫描重构剖面是否符合实际钻井情况，符合则利用所述分频段子波数据中反应岩溶储层的频带数据进行岩溶储层预测；

本实施例中，利用已知钻井情况进行分频扫描重构的剖面分析以及验证重构的剖面是否吻合具体包括：

(3)根据工区内钻、测井资料以及实际的地层结构，并考虑茅口组岩溶储层发育特征建立地质模型，如图3-图5所示；图3为地层结构；图4为茅口组顶部发育岩溶储层；图5去除上覆地层结构；其中，图3仅考虑了地层结构，未考虑茅口组发育岩溶储层；图4是在图3的基础上考虑了距离茅口组顶20-40m发育厚度为10-15m的岩溶型储层；图5为去除上覆地层结构的岩溶型储层模型；

图6-图8为本实施例中茅口组岩溶储层模型正演剖面，其中，图6为地层结构；图7为茅口组顶部发育岩溶储层；图8为去除上覆地层结构。

利用模型正演确定不同地质模型下的地震响应特征，将不同模型的正演结果与分频段子波扫描重构的剖面结果对比验证重构的剖面是否与实钻吻合。

本实施例中分频段子波扫描重构剖面不符合实际钻井情况，则对不同频段扫描重构后的数据与模型正演的剖面对比，剔除强反射层频段范围进行子波重构，如果分频段重构的结果能够消除强反射的影响，与模型正演中图8结果吻合，直接采用该段频率子波进行重构。

如果不能压制强反射影响突出岩溶特征，在分频段子波重构的基础上开展波形分解，波形分解是对多子波分解后的数据体将具有某一相似特征的地震波形统计分类，具体的，本实施例中筛选出岩溶响应的波形分量及对筛选出的波形分量后进行波形分解重构包括：

首先，将输入的地震数据体运用奇异值分解的方式得到不同的解，即不同的特征波形，每个解对应一个特征值。对所有解按其特征值平方和排列，即按能量大小排序。如第一分量代表了地震中能量最强的波形，能量次之的为第二分量，…，依此类推。每一分量代表具有相似特征的地震响应，具有相似特征的地震响应能够反映一定特征的地质信息，比如地质界面、异常体等。所有的分量中，第一分量能量最强，通常为该段数据体贡献最大地质信息的地震响应，第二分量反映次一级的地质特征，…，依此类推。通过波形分解得到不同波形的一系列不同分量的结果，选取其中的分量或者分量组合与模型正演的结果对比，筛选出能够代表岩溶响应的波形分量或分量组合重构数据体。

(4)频率重构的基础上，沿拾取的层位设置时窗，该时窗包含岩溶储层的情况下尽量小，通常根据储层与拾取的强反射层距离大小而定，在该时窗内进行波形分解，波形分解按不同的波形得到一系列不同的结果，按能量大小排序为第一分量、第二分量、第三分量…，能够反映不同的地质信息，从中筛选出某一分量或某几个分量组合进行重构，将重构结果与正演模拟的结果对比，选出与正演结果吻合度最高的分量或分量组合重构，代表岩溶储层的地震响应。

例如，四川盆地川南地区某工区经钻井证实茅口组岩溶储层被上覆的强反射层上二叠底界屏蔽，实际的钻井与地震剖面不符合，利用基于子波重构的技术能够有效识别岩溶储层。图9-图11为四川南部某工区原始地震剖面和多子波重构前后的剖面对比，其中，图9为原始剖面；图10为30-40hz子波重构后剖面；图11为30-40hz子波重构+波形分解第3、4分量剖面。

井A茅三段钻遇岩溶储层获高产气，井B未钻遇岩溶储层。从图9可以看出原始地震剖面上过井A、井B处茅口组顶部均为强波峰反射；图10是30-40hz子波重构剖面，井点B处消除了上二叠底界强波峰影响为空白反射，井点A处上二叠底界强波峰影响部分得到压制，残留部分与岩溶储层子波频率域重合的地震响应。图11为基于30-40hz子波重构后沿上二叠底界上下15ms时窗的波形分解后第3和第4分量重构剖面，井点A处为弱波峰反射，与模型正演结果一致，压制了上二叠底界的屏蔽作用，突出了茅口组岩溶储层的地震响应。本发明实施例能够有效识别强反射层下的岩溶储层地震响应。

(5)首先，将工区内数据体重构，在地震属性中优选出最能反映岩溶储层特征的属性沿层提取该地震属性获得其平面分布图，确定岩溶储层发育有利区。然后在该基础上对重构数据体开展以频谱、频谱衰减分析为主的含油气性检测，获取含油气范围有利区，结合之前岩溶储层预测有利区预测平面图最终确定含油气的范围。

例如，四川盆地川南地区某工区传统的地震属性分析技术与应用多子波分解与重构后的地震属性分析的一个实例。图12、图13为目的层相同时窗的原始地震数据和应用地震子波重构技术的均方根振幅属性平面图，图12为原始地震数据均方根振幅；图13为子波重构后的均方根振幅。

本实施例中，井A、D、E均钻遇岩溶储层，井B、井C未钻遇岩溶储层。

图12为传统的地震属性分析技术，并未对地震数据本身做任何处理开展的地震属性储层预测工作，可以看出因受到上覆强反射层影响中五口井均位于强振幅区域，并不能确定岩溶储层的发育有利区；图13应用多子波分解与重构技术后的地震属性分析储层预测，从平面图可以看到，该方法消除了上覆地层因影响，其中井A、D、E位于强振幅区域，井B、井C位于弱振幅区，与实钻吻合，能够用于预测强反射层下岩溶储层发育的有利区。

本发明在于对用于储层预测和含油气检测的地震数据做了提前处理，从输入源上提高了数据的精度和可靠性，采用本发明的技术得到的地震数据能够准确的识别岩溶储层的地震响应，因此，相比传统的地震属性预测和含油气性检测直接在原始的地震处理后的成果数据上预测，其本身提供的地震数据体上都不能很好的识别出岩溶储层的响应，必然会导致预测精度的下降。本发明基于地震子波理论，从用于预测的数据体出发，重构出新的数据体，综合多种信息将地震数据上的岩溶响应识别出来，相对与井资料应用地震数据能够更好的反映岩溶储层的空间变化的优势，达到预测岩溶储层有利区的目的，这样做能够大大提高预测精度。

本发明通过分频扫描和分波段的进行波形分解重构确定岩溶储层的子波频段响应段，采用精细分区频段扫描分析，多井验证空间上解决岩溶地震响应精度问题。合理筛选波形分量表征岩溶地震响应特征。利用钻遇岩溶的钻、测井资料，结合该区地质认识，正演岩溶的地震模型，与波形分量对比分析，确定岩溶响应的波形分量，提高预测精度。

同时，本发明还提供一种岩溶储层预测装置，如图14所示包括：

子波分解模块401，用于对目的层地震数据进行多子波分解确定分频段子波数据；

分频重构模块402，用于对所述分频段子波数据进行分频扫描重构确定分频段子波扫描重构剖面；

分频判断模块403，用于根据实际钻井数据判断所述分频段子波扫描重构剖面是否符合实际钻井情况；

分频预测模块404，用于确定所述分频段子波扫描重构剖面符合实际钻井情况，则利用所述分频段子波数据中反应岩溶储层的频带数据进行岩溶储层预测；

岩溶响应特征判断模块405，用于判断所述分频段子波扫描重构剖面不符合实际钻井情况，则判断所述分频段子波扫描重构剖面是否反应岩溶响应特征；

岩溶响应子波重构模块406，用于确定反应岩溶响应特征，则对岩溶响应频段的子波数据进行数据体重构，利用重构后的岩溶响应频段的子波数据进行岩溶储层预测；

波形分解重构模块407，用于确定不反应岩溶响应特征，则对分频扫描重构后的子波数据进行波形分解重构确定岩溶响应最佳分量重构数据体，利用确定的溶响应最佳分量重构数据体进行岩溶储层预测。

本发明实施例中，所述的子波分解模块对目的层地震数据进行多子波分解确定分频段子波数据包括：

对目的层地震数据进行雷克子波、Paul子波或Morlet子波分解确定分频段子波数据。

本发明实施例中，所述的分频判断模块包括：

建模单元，用于根据目的层工区的实际钻井数据建立地质模型；

响应特征确定单元，用于利用建立的地质模型进行正演模拟确定地震响应特征；

响应判断单元，用于根据所述地震响应特征判断分频段子波扫描重构剖面是否符合实际钻井情况。

本发明实施例中，所述的子波分解模块包括：

频谱分析单元，用于对目的层地震数据进行频谱分析确定多子波分解的频带范围；

分解单元，用于根据确定的频带范围进行多子波分解确定分频段子波数据。

通过前面实施例的描述，对本领域技术人员而言，可清楚获知本申请岩溶储层预测装置的实施方式，在此不再赘述。

本实施例还提供一种电子设备，该电子设备可以是台式计算机、平板电脑及移动终端等，本实施例不限于此。在本实施例中，该电子设备可以参照前述方法及装置的实施例，其内容被合并于此，重复之处不再赘述。

图15为本发明实施例的电子设备600的系统构成的示意框图。如图15所示，该电子设备600可以包括中央处理器100和存储器140；存储器140耦合到中央处理器100。值得注意的是，该图是示例性的；还可以使用其他类型的结构，来补充或代替该结构，以实现电信功能或其他功能。

一实施例中，岩溶储层预测方法功能可以被集成到中央处理器100中。其中，中央处理器100可以被配置为进行如下控制：

对目的层地震数据进行多子波分解确定分频段子波数据；

对所述分频段子波数据进行分频扫描重构确定分频段子波扫描重构剖面；

根据实际钻井数据判断所述分频段子波扫描重构剖面是否符合实际钻井情况；

确定所述分频段子波扫描重构剖面符合实际钻井情况，则利用所述分频段子波数据中反应岩溶储层的频带数据进行岩溶储层预测；

判断所述分频段子波扫描重构剖面不符合实际钻井情况，则判断所述分频段子波扫描重构剖面是否反应岩溶响应特征；

确定反应岩溶响应特征，则对岩溶响应频段的子波数据进行数据体重构，利用重构后的岩溶响应频段的子波数据进行岩溶储层预测；

确定不反应岩溶响应特征，则对分频扫描重构后的子波数据进行波形分解重构确定岩溶响应最佳分量重构数据体，利用确定的岩溶响应最佳分量重构数据体进行岩溶储层预测。

本发明实施例中，所述的对目的层地震数据进行多子波分解确定分频段子波数据包括：

对目的层地震数据进行雷克子波、Paul子波或Morlet子波分解确定分频段子波数据。

本发明实施例中，所述的根据实际钻井数据判断所述分频段子波扫描重构剖面是否符合实际钻井情况包括：

根据目的层工区的实际钻井数据建立地质模型；

利用建立的地质模型进行正演模拟确定地震响应特征；

根据所述地震响应特征判断分频段子波扫描重构剖面是否符合实际钻井情况。

本发明实施例中，所述的对目的层地震数据进行多子波分解确定分频段子波数据包括：

对目的层地震数据进行频谱分析确定多子波分解的频带范围；

根据确定的频带范围进行多子波分解确定分频段子波数据。

如图15所示，该电子设备600还可以包括：通信模块110、输入单元120、音频处理单元130、显示器160、电源170。值得注意的是，电子设备600也并不是必须要包括图15中所示的所有部件；此外，电子设备600还可以包括图15中没有示出的部件，可以参考现有技术。

如图15所示，中央处理器100有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器100接收输入并控制电子设备600的各个部件的操作。

其中，存储器140，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器100可执行该存储器140存储的该程序，以实现信息存储或处理等。

输入单元120向中央处理器100提供输入。该输入单元120例如为按键或触摸输入装置。电源170用于向电子设备600提供电力。显示器160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。

该存储器140可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器140还可以是某种其它类型的装置。存储器140包括缓冲存储器141(有时被称为缓冲器)。存储器140可以包括应用/功能存储部142，该应用/功能存储部142用于存储应用程序和功能程序或用于通过中央处理器100执行电子设备600的操作的流程。

存储器140还可以包括数据存储部143，该数据存储部143用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器140的驱动程序存储部144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。

通信模块110即为经由天线111发送和接收信号的发送机/接收机110。通信模块(发送机/接收机)110耦合到中央处理器100，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。

基于不同的通信技术，在同一电子设备中，可以设置有多个通信模块110，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)110还经由音频处理器130耦合到扬声器131和麦克风132，以经由扬声器131提供音频输出，并接收来自麦克风132的音频输入，从而实现通常的电信功能。音频处理器130可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器130还耦合到中央处理器100，从而使得可以通过麦克风132能够在本机上录音，且使得可以通过扬声器131来播放本机上存储的声音。

本发明实施例还提供一种计算机可读程序，其中当在电子设备中执行所述程序时，所述程序使得计算机在所述电子设备中执行如上面实施例所述的岩溶储层预测方法。

本发明实施例还提供一种存储有计算机可读程序的存储介质，其中所述计算机可读程序使得计算机在电子设备中执行上面实施例所述的岩溶储层预测。

以上参照附图描述了本发明的优选实施方式。这些实施方式的许多特征和优点根据该详细的说明书是清楚的，因此所附权利要求旨在覆盖这些实施方式的落入其真实精神和范围内的所有这些特征和优点。此外，由于本领域的技术人员容易想到很多修改和改变，因此不是要将本发明的实施方式限于所例示和描述的精确结构和操作，而是可以涵盖落入其范围内的所有合适修改和等同物。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

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

Claims (10)

1.一种岩溶储层预测方法，其特征在于，所述的方法包括：

对目的层地震数据进行多子波分解确定分频段子波数据；

对所述分频段子波数据进行分频扫描重构确定分频段子波扫描重构剖面；

根据实际钻井数据判断所述分频段子波扫描重构剖面是否符合实际钻井情况；

确定所述分频段子波扫描重构剖面符合实际钻井情况，则利用所述分频段子波数据中反应岩溶储层的频带数据进行岩溶储层预测；

确定所述分频段子波扫描重构剖面不符合实际钻井情况，则判断所述分频段子波扫描重构剖面是否反应岩溶响应特征；

确定反应岩溶响应特征，则对岩溶响应频段的子波数据进行数据体重构，利用重构后的岩溶响应频段的子波数据进行岩溶储层预测；

确定不反应岩溶响应特征，则对分频扫描重构后的子波数据进行波形分解重构确定岩溶响应最佳分量重构数据体，利用确定的岩溶响应最佳分量重构数据体进行岩溶储层预测。

2.如权利要求1所述的岩溶储层预测方法，其特征在于，所述的对目的层地震数据进行多子波分解确定分频段子波数据包括：

对目的层地震数据进行雷克子波、Paul子波或Morlet子波分解确定分频段子波数据。

3.如权利要求1所述的岩溶储层预测方法，其特征在于，所述的根据实际钻井数据判断所述分频段子波扫描重构剖面是否符合实际钻井情况包括：

根据目的层工区的实际钻井数据建立地质模型；

利用建立的地质模型进行正演模拟确定地震响应特征；

根据所述地震响应特征判断分频段子波扫描重构剖面是否符合实际钻井情况。

4.如权利要求1所述的岩溶储层预测方法，其特征在于，所述的对目的层地震数据进行多子波分解确定分频段子波数据包括：

对目的层地震数据进行频谱分析确定多子波分解的频带范围；

根据确定的频带范围进行多子波分解确定分频段子波数据。

5.一种岩溶储层预测装置，其特征在于，所述的装置包括：

子波分解模块，用于对目的层地震数据进行多子波分解确定分频段子波数据；

分频重构模块，用于对所述分频段子波数据进行分频扫描重构确定分频段子波扫描重构剖面；

分频判断模块，用于根据实际钻井数据判断所述分频段子波扫描重构剖面是否符合实际钻井情况；

分频预测模块，用于确定所述分频段子波扫描重构剖面符合实际钻井情况，则利用所述分频段子波数据中反应岩溶储层的频带数据进行岩溶储层预测；

岩溶响应特征判断模块，用于判断所述分频段子波扫描重构剖面不符合实际钻井情况，则判断所述分频段子波扫描重构剖面是否反应岩溶响应特征；

岩溶响应子波重构模块，用于确定反应岩溶响应特征，则对岩溶响应频段的子波数据进行数据体重构，利用重构后的岩溶响应频段的子波数据进行岩溶储层预测；

波形分解重构模块，用于确定不反应岩溶响应特征，则对分频扫描重构后的子波数据进行波形分解重构确定岩溶响应最佳分量重构数据体，利用确定的溶响应最佳分量重构数据体进行岩溶储层预测。

6.如权利要求5所述的岩溶储层预测装置，其特征在于，所述的子波分解模块对目的层地震数据进行多子波分解确定分频段子波数据包括：

对目的层地震数据进行雷克子波、Paul子波或Morlet子波分解确定分频段子波数据。

7.如权利要求5所述的岩溶储层预测装置，其特征在于，所述的分频判断模块包括：

建模单元，用于根据目的层工区的实际钻井数据建立地质模型；

响应特征确定单元，用于利用建立的地质模型进行正演模拟确定地震响应特征；

响应判断单元，用于根据所述地震响应特征判断分频段子波扫描重构剖面是否符合实际钻井情况。

8.如权利要求5所述的岩溶储层预测装置，其特征在于，所述的子波分解模块包括：

频谱分析单元，用于对目的层地震数据进行频谱分析确定多子波分解的频带范围；

分解单元，用于根据确定的频带范围进行多子波分解确定分频段子波数据。

9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一项所述方法。

10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至4任一项所述方法的计算机程序。

Images