一种针对分布于煤岩间受强反射影响的砂岩储层预测方法

技术领域

本发明属于三维地震储层预测技术领域，涉及一种针对分布于煤岩间受强反射影响的砂岩储层预测方法。

背景技术

煤岩是油气藏最主要的烃源岩之一，由于煤岩岩石骨架本身具有低速度、低密度等特征，容易与围岩形成强烈的反射界面，在三维地震中通常表现为横向连续稳定的低频强反射现象。已有相关研究证实，煤岩强反射数倍于储层反射，使得煤岩之间的储层反射信息被严重干涉，储层预测无法获得准确结果。为解决这一问题，前人重点采用井震结合三维建模、定量反演预测等技术手段，并在一定范围内取得了效果。但是这些技术通常容易受井控程度、钻井分布均匀程度等诸多条件限制，而且建模或反演过程中参数设置的人为因素影响较大。除此之外，由于煤岩岩性相对较脆，钻井过程中往往伴有不同程度的井壁垮塌，使得测井资料难以获得煤岩准确岩电参数，测井资料的应用受到了极大限制，这些都为煤系地层砂岩储层预测带来困难，特别是对于分布在煤岩之间的砂岩储层预测，难度更大。

因此，如何提供一种针对分布于煤岩间受强反射影响的砂岩储层预测方法，以解决现有技术对于分布在煤岩之间的砂岩储层预测难度大、储层预测无法获得准确结果等缺陷，实已成为本领域技术人员亟待解决的技术问题。

发明内容

鉴于以上所述现有技术的缺点，本发明的目的在于提供一种针对分布于煤岩间受强反射影响的砂岩储层预测方法，用于解决现有技术对于分布在煤岩之间的砂岩储层预测难度大、储层预测无法获得准确结果的问题。

为实现上述目的及其他相关目的，本发明提供一种针对分布于煤岩间受强反射影响的砂岩储层预测方法，包括：分别统计目标区已钻井煤岩层数、各层厚度、煤间地层厚度参数，以及煤岩之间砂岩、储层的叠合厚度参数；对目标区三维地震数据进行分频处理，选取对煤岩识别最有效的数据；拟合煤岩之间地层与煤岩反射波时差，实现煤岩反射波时差平面数据从时间域到深度域的转换；分别拟合煤岩之间地层与砂岩、砂岩与储层叠合厚度之间的关系式，筛选满足误差要求的钻井，并对煤岩间地层、砂岩及储层厚度平面数据进行分步校正，实现逐级误差控制。

于本发明的一实施例中，对所述目标区三维地震数据进行分频处理，以选取对煤岩识别最有效的三维地震数据的步骤包括：对所述目标区三维地震数据进行分频处理，形成不同频段地震数据体，并分别提取已钻井点处各套煤岩不同频段的反射波振幅值，形成振幅数据集；将煤岩厚度与不同频段的反射波振幅值进行散点交汇，选取相关度最高的拟合方式所对应的数据对，形成第一拟合关系式和第一相关度，作为分频参数优选的依据。对所述对煤岩识别最有效的三维地震数据中煤岩反射波进行层位解释，对所解释的煤岩反射波时间层位进行相减并取绝对值，形成煤间反射波时差平面数据。

于本发明的一实施例中，拟合煤岩之间地层厚度与煤岩反射波时差之间的关系式，并优选出相关度最高的拟合关系式，以实现煤岩反射波时差平面数据从时间域到深度域的转换，并筛选满足误差率要求的已钻井进行校正的步骤包括：利用已筛选对煤岩识别最有效的三维地震数据，提取各井点处所述煤间反射波时差数值，并与实钻煤岩之间地层厚度进行散点交汇，选取相关度最高的拟合方式，形成第二拟合关系式及第二相关度；将所述煤间反射波时差平面数据按照第二拟合关系式进行计算，以获取校正前的深度域煤间地层厚度平面数据；将满足第一误差率要求的已钻井定义为第一校正井集合，并对校正前的深度域煤间地层厚度平面数据进行校正，形成校正后的深度域煤间地层厚度平面数据。

于本发明的一实施例中，分别拟合煤岩之间地层与砂岩、砂岩与储层的叠合厚度之间的关系式，优选出相关对最高的拟合关系式，分步筛选满足误差率要求的已钻井，进而对砂岩及储层厚度进行分步校正，实现逐级误差控制的步骤包括：将已有钻井的煤岩之间的砂岩叠合厚度与煤岩之间地层厚度进行散点交汇，选取相关度最高的拟合方式，形成第三拟合关系式及第三相关度；将校正后的深度域煤间地层厚度平面数据按照第三拟合关系式进行计算，以获取校正前砂岩叠合厚度平面数据；将满足第二误差率要求的已钻井定义为第二校正井集合，并对校正前砂岩叠合厚度平面数据进行校正，形成校正后砂岩叠合厚度平面数据。

于本发明的一实施例中，分别拟合煤岩之间地层与砂岩、砂岩与储层的叠合厚度之间的关系式，优选出相关对最高的拟合关系式，分步筛选满足误差率要求的已钻井，进而对砂岩及储层厚度进行分步校正，实现逐级误差控制的步骤还包括：将已有钻井的煤岩之间的砂岩叠合厚度与储层叠合厚度进行散点交汇，选取相关度最高的拟合方式，形成第四拟合关系式及第四相关度；对校正后砂岩叠合厚度平面数据按照第四拟合关系式进行计算，得到校正前储层叠合厚度平面数据；将满足第三误差率要求的已钻井定义为第三校正井集合，并对校正前砂岩叠合厚度平面数据进行校正，以获取煤间最终储层叠合厚度分布结果。

于本发明的一实施例中，第一误差率要求是指第一误差率小于(1-第二相关度)*100％；第二误差率要求是指第二误差率小于(1-第三相关度)*100％；第三误差率要求是指第三误差率小于(1-第四相关度)*100％。

于本发明的一实施例中，第一误差率是指指已钻井点处经拟合关系式函数计算所得到的地层厚度与实钻地层厚度差值，取绝对值后与相应实钻厚度的百分比；第二误差率是指指已钻井点处经拟合关系式函数计算所得到的砂岩叠合厚度与实钻砂岩叠合厚度差值，取绝对值后与相应实钻厚度的百分比；第三误差率是指已钻井点处经拟合关系式函数计算所得到的储层叠合厚度与实钻储层叠合厚度差值，取绝对值后与相应实钻厚度的百分比。

如上所述，本发明所述的针对分布于煤岩间受强反射影响的砂岩储层预测方法，具有以下有益效果：

本发明所述针对分布于煤岩间受强反射影响的砂岩储层预测方法通过对地震资料的优化处理及优选，明确对煤岩最敏感的地震数据，以落实煤岩之间的地层时间域厚度展布，进而建立地层与砂岩、砂岩与储层之间的最佳拟合关系式，进行递推和逐步校正。由于只利用各类岩性的厚度信息，从而避免了由于井径垮塌等造成的测井曲线不真实，直接建模或反演失真的问题，同时也有效削弱了建模或反演过程中参数设置的人为因素影响。

附图说明

图1显示为本发明的针对分布于煤岩间受强反射影响的砂岩储层预测方法于一实施例中的实施流程图。

图2显示为本发明的针对分布于煤岩间受强反射影响的砂岩储层预测方法的方法流程示意图。

元件标号说明

S11～S20 步骤

具体实施方式

以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。

实施例一

本实施例提供一种针对分布于煤岩间受强反射影响的砂岩储层预测方法，包括：

根据岩性及储层划分标准，对目标区已有钻井的煤系地层进行测井解释，分别统计各钻井的煤岩层数、各层厚度、煤间地层厚度参数，以及煤岩之间砂岩、储层的叠合厚度参数；

对所述目标区三维地震数据进行分频处理，以选取对煤岩识别最有效的三维地震数据；

拟合煤岩之间地层厚度与煤岩反射波时差之间的关系式，并优选出相关度最高的拟合关系式，以实现煤岩反射波时差平面数据从时间域到深度域的转换，并筛选满足误差率要求的已钻井进行校正；

分别拟合煤岩之间地层与砂岩、砂岩与储层的叠合厚度之间的关系式，优选出相关对最高的拟合关系式，分步筛选满足误差率要求的已钻井，进而对砂岩及储层厚度进行分步校正，实现逐级误差控制。

以下将结合图示对本实施例所提供的针对分布于煤岩间受强反射影响的砂岩储层预测方法进行详细描述。请参阅图1，显示为针对分布于煤岩间受强反射影响的砂岩储层预测方法于一实施例中的流程示意图。如图1所示，所述针对分布于煤岩间受强反射影响的砂岩储层预测方法具体包括以下几个步骤：

S11，根据岩性及储层划分标准，对目标区已有钻井的煤系地层进行测井解释，分别统计各钻井的煤岩层数、各层厚度、煤间地层厚度参数，以及煤岩之间砂岩、储层的叠合厚度参数。

所述S11中的岩性及储层划分标准需根据目标区实际沉积特征、储层特征、测录井及试油气数据制定，不具有统一标准。

以某工区为例，该地区已有完钻井20余口，均于煤系地层钻遇三套主力煤岩，由于沼泽相沉积环境相对稳定，因此煤岩横向分布稳定井间对比可靠，仅煤岩厚度、煤间地层厚度及砂岩、储层叠合厚度有所差异。根据实施例所处区域地质沉积背景，将三套煤岩编号为5号、6号及7号煤，本实施例以5号及6号煤之间储层识别为例。

例如，根据岩性及储层划分标准，对目标区已有钻井的煤岩地层进行测井解释，统计各完钻井5号及6号煤岩厚度，以及5-6号煤之间的地层厚度；统计5-6号煤之间砂岩及储层叠合厚度。

S12，对所述目标区三维地震数据进行分频处理，形成不同频段地震数据体，并分别提取已钻井点处各套煤岩不同频段的反射波振幅值，形成振幅数据集。

所述S12中的分频处理需视实际地震资料品质进行参数设置，经分频处理后的地震资料需能够将各套煤岩反射波在纵向上有所区分。

在本实施例中，所述S12中的反射波振幅值需依据井震时深标定关系而定，若煤岩标定于波峰则统一读取最大振幅值，若标定于波谷则统一读取最小振幅值。

S13，将S11中的煤岩厚度与S12中不同频段的反射波振幅值进行散点交汇，选取相关度最高的拟合方式所对应的数据对，形成第一拟合关系式和第一相关度，作为分频参数优选的依据。在本实施例中，拟合方式包括对比线性、二项式、三项式等。

例如，目标区三维地震数据有效频带6-45Hz，主频22-25Hz，由于煤岩反射主要为低频强反射特征，经分析，目标区20Hz以内地震数据分辨率较低，难以将5号及6号煤岩反射波完全分离开，读取振幅值容易产生混淆，即视20Hz以下为无效反射频段；对20-45Hz地震资料以2Hz为间隔、10Hz为长度进行分频，形成8个有效分频数据体(20-30Hz、22-32Hz...34-44Hz)。

根据井震标定结果，5号及6号煤岩均标定于反射波波谷范围内，分别提取上述8个分频数据体煤岩所对应的反射波最小振幅值，形成最小振幅数据集。

将5号及6号煤岩实钻厚度与8个分频数据所提取的最小振幅值分别进行散点交汇，并分别对比线性、二项式、三项式等拟合方式，最终确定26-36Hz所对应的e指数函数相关度最高(作为第一拟合关系式)，达到0.81(作为第一相关度)，因此，选取26-36Hz频段，作为对煤岩反应最敏感的地震数据。

S14，将S13中优选的分频段地震数据中煤岩反射波进行层位解释，对所解释的煤岩反射波时间层位进行相减并取绝对值，形成煤间反射波时差平面数据。

S15，根据S14所得煤间反射波时差平面数据，分别提取各井点处所述煤间反射波时差平面数据，并与实钻煤岩之间地层厚度进行散点交汇，选取相关度最高的拟合方式(拟合方式包括对比线性、二项式、三项式等)，形成第二拟合关系式及第二相关度。

例如，对步骤S14所确定的26-36Hz分频地震数据中5号和6号煤岩反射波的最小振幅进行层位解释，并取两套煤岩反射波层位时间差值的绝对值，形成5-6号煤之间的时差平面数据；提取各井点处煤间反射波时差数值，并与实钻5-6号煤间地层厚度进行散点拟合，优选相关度最高达到0.84的三项式拟合，作为第二相关度及第二拟合关系式。

S16，将S14中所述煤间反射波时差平面数据按照第二拟合关系式进行计算，以获取校正前的深度域煤间地层厚度平面数据；优选满足第一误差率要求的已钻井作为第一校正井点集合，对校正前的深度域煤间地层厚度平面数据进行校正，形成校正后的深度域煤间地层厚度平面数据。

第一误差率是指指已钻井点处经拟合关系式函数计算所得到的地层厚度与实钻地层厚度差值，取绝对值后与相应实钻厚度的百分比。

第一误差率要求是指第一误差率小于(1-第二相关度)*100％。

例如，按照第二拟合关系式，对煤间反射波时差平面数据进行计算，得到校正前的深度域煤间地层厚度平面数据，分别提取已钻井点处对应的地层厚度数值；并与已钻井实钻5-6号煤间地层厚度取差值的绝对值，得到每口井煤间地层厚度误差率；

由于26-36Hz地震振幅数值对煤岩厚度的最大相关度为0.84，即该分频数据对5-6号煤系地层识别准确度，因此实钻井与预测可允许的最大误差率为(1-0.84)*100％＝16％。

将实钻煤间地层厚度与预测煤间地层厚度误差率小于16％的井作为第一校正井集合，并对所得校正前的深度域煤间地层厚度平面数据进行校正。误差率大于16％的井不可参与校正，否则容易造成局部畸变。

S17，根据S11中的识别结果，将已有钻井煤岩之间的砂岩叠合厚度与煤间地层厚度进行散点交汇，选取相关度最高的拟合方式(拟合方式包括对比线性、二项式、三项式等)，形成第三拟合关系式及第三相关度；将S16中校正后的深度域煤间地层厚度平面数据按照第三拟合关系式进行计算，以获取校正前砂岩叠合厚度平面数据。

例如，建立S11中钻井实测的5-6号煤之间地层厚度以及砂岩叠合厚度散点交汇，并优选相关度最高达到0.79的二项式拟合方式，作为第三相关度及第三拟合关系式；并对S16中校正后的深度域平面数据进行二项式函数计算，得到校正前煤间砂岩叠合厚度平面数据。

S18，优选满足第二误差率要求的已钻井作为第二校正井点集合，对S17中校正前砂岩叠合厚度平面数据进行校正，形成校正后砂岩叠合厚度平面数据。

第二误差率是指已钻井点处经拟合关系式函数计算所得到的砂岩叠合厚度与实钻砂岩叠合厚度差值，取绝对值后与相应实钻砂岩叠合厚度的百分比。

第二误差率要求是指第二误差率小于(1-第三相关度)*100％。

例如，建立钻井实测的5-6号煤之间地层厚度以及砂岩叠合厚度散点交汇，并优选相关度最高达到0.79的二项式拟合方式，作为第三相关度及第三拟合关系式；并对校正后砂岩厚度平面数据进行函数计算，得到校正前砂岩叠合厚度平面数据；

优选误差率小于21％的井作为第二校正井集合，并对校正前砂岩叠合厚度平面数据进行校正，形成校正后砂岩叠合厚度平面数据。

S19，根据S11中的识别结果，将已有钻井的煤岩之间的砂岩叠合厚度与储层叠合厚度进行散点交汇，选取相关度最高的拟合方式(拟合方式包括对比线性、二项式、三项式等)，形成第四拟合关系式及第四相关度；对S18中校正后砂岩叠合厚度平面数据按照第四拟合关系式进行计算，得到校正前储层叠合厚度平面数据。

例如，建立已钻井煤岩之间的砂岩叠合厚度与储层叠合厚度的散点交汇，并优选相关度为0.89的线性拟合方式，作为第四相关度及第四拟合关系式；并对S19中的砂岩叠合厚度平面数据进行线性函数计算，得到校正前煤间储层叠合厚度平面数据。

S20，优选满足第三误差率要求的已钻井作为第三校正井点集合，对校正前储层叠合厚度平面数据进行校正，以获取煤间最终储层叠合厚度分布结果。

第三误差率是指已钻井点处经拟合关系式函数计算所得到的储层厚度与实钻储层厚度差值，取绝对值后与相应实钻厚度的百分比。

第三误差率要求是指第三误差率小于(1-第四相关度)*100％。

例如，建立已钻井煤岩之间的砂岩叠合厚度与储层叠合厚度的散点交汇，并优选相关度为0.89的线性拟合方式，作为第四相关度及第四拟合关系式；并对校正后砂岩叠合厚度平面数据进行线性函数计算，得到校正前砂岩叠合厚度平面数据。

优选误差率小于11％的井作为第三校正井集合，并对校正前储层叠合厚度平面数据进行校正，得到煤间最终储层叠合厚度分布结果。

综上所述，本发明所述分布于煤岩之间受强反射影响的砂岩储层预测方法，通过对地震资料的优化处理及优选，明确对煤岩最敏感的地震数据，以落实煤岩之间的地层时间域厚度展布，进而建立地层与砂岩、砂岩与储层之间的最佳拟合关系式，进行递推和逐步校正。由于只利用各类岩性的厚度信息，从而避免了由于井径垮塌等造成的测井曲线不真实，直接建模或反演失真的问题，同时也有效削弱了建模或反演过程中参数设置的人为因素影响。本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。