基于敏感测井曲线的高分辨率反演方法与流程

1.本发明涉及地球物理反演领域，尤其涉及一种基于敏感测井曲线的高分辨率反演方法。


背景技术：

2.地震反演(seismic inversion)是利用地表观测地震资料，以已知地质规律和钻井、测井资料为约束，对地下岩层空间结构和物理性质进行成像(求解)的过程。在地震反演过程中，通过地面上的一系列点上，利用人工激发地震波，地震波向地下传播，当遇到波阻抗界面(即上下地层波阻抗不相等面)时，在波阻抗界面上地震波产生反射现象，地震波传播方向发生改变，地震波开始向上传播，在地面上的一系列接收点上安置着接收器，接收向上传播的地震波。
3.目前，常采用基于褶积模型的叠后声波阻抗反演方法。然而该方法只利用了声波测井信息和密度测井信息，并没有使用电性和放射性等测井资料，因此，直接利用这种声波阻抗进行反演难以区分地下储层和非储层，从而无法有效识别岩性，影响了对储层的认识。


技术实现要素：

4.本发明的目的在于提供一种基于敏感测井曲线的高分辨率反演方法，旨在解决目前常规阻抗反演方法对岩性识别不清、储层厚度预测不准、分辨率较低的问题。
5.第一方面，本发明实施例提供了一种基于敏感测井曲线的高分辨率反演方法，包括：
6.获取井曲线数据；
7.对井曲线数据进行预处理，获得标准化曲线数据；
8.分析并筛选标准化曲线数据，得到敏感曲线数据；
9.对敏感曲线数据进行曲线重构，形成拟声波曲线数据；
10.对拟声波曲线数据进行地震反演。
11.可选的，曲线重构的步骤包括：提取敏感曲线数据，敏感曲线数据包括对应低、中、高频段的多条敏感曲线；计算各频段内敏感曲线的匹配程度；比较匹配程度的大小，若匹配程度的数值大，则将相应敏感曲线映射至原频段的声波曲线数据中，形成拟声波曲线数据。
12.可选的，匹配程度的计算公式为：
[0013][0014]
其中，表示第i个频段第j条敏感曲线，表示第i个频段第m条敏感曲线，j、m都属于敏感曲线集a，cor表示敏感曲线之间的相关程度，x n表示目的层段内，p
i，j
表示第i个频段第j条敏感曲线整体的匹配程度。
[0015]
可选的，地震反演的步骤包括：标定精细储层；根据精细储层，约束拟声波曲线数据以构建地质框架模型；结合地质框架模型进行地震反演。
[0016]
可选的，标定精细储层的方法包括：选择较强的地震反射界面作为标准层；利用标准子波进行合成地震记录初始标定；根据大套地层对应的波组关系提取、统计井旁道地震数据体中的子波并进行二次标定；通过慢度曲线监控调整内部反射波组的对应关系。
[0017]
可选的，还包括：在进行地震反演之前，根据建立的地质框架模型，约束井阻抗的动态变化范围并进行内插外推，以确定约束条件；在进行地震反演时，根据约束条件，对地质框架模型进行稀疏约束脉冲反演。
[0018]
可选的，包括对井曲线数据进行预处理。
[0019]
可选的，预处理的步骤包括：异常值处理、曲线重采样、环境校正、基线校正、归一化、标准化。
[0020]
可选的，对井曲线数据进行预处理，包括：对扩径井段采集的声波测井曲线数据和密度测井曲线数据进行环境校正；和/或，对自然电位测井曲线数据进行基线校正。
[0021]
可选的，对经预处理后的井曲线数据进行标准化处理。
[0022]
本发明的有益效果在于：通过获取对岩性较为敏感的敏感曲线数据，再将敏感曲线数据加入到原声波曲线数据中，以获得新的拟声波曲线数据，再对新获得的拟声波曲线数据进行地震反演，以提高岩性的识别度和储层厚度的准确度，从而获得能反映储层的三维阻抗体，进而便于对储层进行预测。
附图说明
[0023]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0024]
图1为研究区内声波阻抗区分岩性直方图；
[0025]
图2示出了本发明实施例1提供的一种基于敏感测井曲线的高分辨率反演方法流程图；
[0026]
图3a-3b示出了本发明实施例1中的曲线预处理前后对比图；
[0027]
图4a-4b示出了本发明实施例1中的电阻率与中子测井交会图、电阻率与声波阻抗交会图；
[0028]
图5a-5c示出了本发明实施例1中的敏感曲线数据与原始声波曲线的融合示意图；
[0029]
图6示出了本发明实施例1中的拟声波曲线数据区分岩性直方图；
[0030]
图7示出了本发明实施例1中的拟声波曲线数据地震反演结果图。
具体实施方式
[0031]
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
[0032]
图1为该研究区内声波阻抗区分岩性直方图。由图1可以看出目的层段声波阻抗对储层区分不敏感，叠置严重，因此，提供一种基于敏感测井曲线的高分辨率反演方法，通过筛选出敏感测井曲线，将反映岩性的敏感曲线信息加入到声波测井曲线中形成新的拟声波曲线数据，并对其开展地震反演，从而获得能反映储层的三维阻抗体，进而进行储层预测。
[0033]
本发明一实施例提供了一种基于敏感测井曲线的高分辨率反演方法，图2为本发明实施例1提供的一种基于敏感测井曲线的高分辨率反演方法流程图，请参考图2，该基于敏感测井曲线的高分辨率反演方法，包括：
[0034]
s01：获取井曲线数据；
[0035]
s02：判断所述井曲线数据是否异常，以获得标准化曲线数据；
[0036]
s03：分析并筛选标准化曲线数据，得到敏感曲线数据；
[0037]
s04：对敏感曲线数据进行曲线重构，形成拟声波曲线数据；
[0038]
s05：对拟声波曲线数据进行地震反演。
[0039]
参考图3至6详细说明本实施例提供的基于敏感测井曲线的高分辨率反演方法。
[0040]
s01，获取井曲线数据。在本实施例中，待研究区内存在至少三口井，井曲线数据需要预先针对研究区内的井分别进行收集。井曲线数据包括对应低、中、高各频段的声波测井曲线数据、密度测井曲线数据、自然电位测井曲线数据、电阻率测井曲线数据等。
[0041]
s02，判断所述井曲线数据是否异常，以获得标准化曲线数据。
[0042]
由于获取的井曲线数据并未经过处理，其中可能包含存在异常的曲线，因此需要判断井曲线数据是否异常，以保证各曲线的合理性。在必要情况下，还需要在判断井曲线数据是否存在异常后，对井曲线数据进行预处理，预处理的步骤包括：异常值处理、曲线重采样、环境校正、基线校正、归一化、标准化。值得注意的是，对于扩径的井段，为了保证后续阻抗反演的正确进行，需要对声波测井曲线数据和密度测井曲线数据进行环境校正；另外，对自然电位测井曲线数据进行基线校正。
[0043]
另外，还需要在此基础上对目的层段的测井曲线进行标准化处理，包括：选取标准层，标准层应选取分布范围稳定、具有一定厚度、曲线特征明显的层段；将标准层上各测井曲线的标准方差和均值归一化，以获得偏移量；将偏移量叠加至经上述校正后的井曲线数据中，从而对其进行标准化处理。需要说明的是，对正常井曲线数据进行标准化处理的过程可参照对经校正后的异常井曲线的标准化处理过程，此处不再赘述。将经标准化处理后的异常井曲线数据和正常井曲线数据汇总形成标准化曲线数据。
[0044]
参考图3a-图3b，以电阻率测井曲线为例进行的曲线预处理前后对比，首先对电阻率取对数，去除异常值，再在此基础上进行标准化处理，图3a为电阻率测井曲线未经预处理图，图3b为电阻率测井曲线经预处理图，由对比可知，经预处理后的电阻率曲线的曲线值范围更为集中。
[0045]
s03：分析、筛选标准化曲线数据，得到敏感曲线数据。
[0046]
在本实施例中，分析并筛选标准化曲线数据的步骤包括：获取测井解释结果；提取标准化曲线数据中的电性、放射性等井曲线数据；根据所述测井解释结果将所述电性、放射性测井曲线数据进行直方图或交会图分析；根据分析结果，筛选出对岩性敏感的敏感曲线数据。需要说明的是，第一处理结果需要预先采用常规反演方法获得。另外，测井解释结果为预先通过常规反演方法获得的地质结果。
[0047]
在实际分析过程中，以将电阻率曲线和中子测井曲线进行交会图分析、将电阻率曲线与声波阻抗曲线进行交会图分析为例，仅对比可知，相比与声波阻抗曲线，电阻率与中子测井对储层的区分度更大，如图4a-图4b所示，由此可见选择电性、放射性等井曲线数据进行直方图或交会图分析，可以得到对岩性更为敏感的敏感曲线数据。
[0048]
s04：对敏感曲线数据进行曲线重构，形成拟声波曲线数据。
[0049]
由于声波曲线不能有效识别研究对象，因此还需要通过一定数学手段重构出一条能反映储层地球物理特征的拟声波曲线，用于地震反演，重构方法采用的是基于小波变换的重构技术，具体包括：获取敏感曲线数据，敏感曲线数据包括对应低、中、高频段的多条敏感曲线；计算各频段内敏感曲线的匹配程度；比较匹配程度的大小，若匹配程度的数值大，则将相应敏感曲线映射至原频段的声波曲线数据中，形成拟声波曲线数据。
[0050]
在本实施例中，在比较匹配程度时，需要多条敏感曲线的匹配程度进行比较，以确定匹配程度数值大的敏感曲线，并将该敏感曲线映射至原频段声波曲线数据中。在其他实施例中，在比较匹配程度时，还可以设定一阈值，将多条敏感曲线的匹配程度的大小与阈值进行比较，当大于该阈值时，可将该将该敏感曲线映射至原频段声波曲线数据中。
[0051]
在本实施例中，匹配程度的计算公式为：
[0052][0053]
其中，表示第i个频段第j条敏感曲线，表示第i个频段第m条敏感曲线，j、m都属于敏感曲线集a，cor表示敏感曲线之间的相关程度，x n表示目的层段内，p
i，j
表示第i个频段第j条敏感曲线整体的匹配程度。
[0054]
应当注意，在进行曲线重构过程中，为便于针对地质目标进行储层反演，选取声波曲线数据进行曲线重构，以作为拟声波曲线数据的低频部分；选取敏感曲线数据进行曲线重构，以作为拟声波曲线的中高频部分。具体而言，以采用电阻率曲线作为敏感测井曲线为例，将电阻率曲线中的高频段与声波曲线的低频段进行融合，如图5a-图5c所示，图中实线为原始曲线，虚线为滤波曲线，图5a为声波曲线，图5b为电阻率曲线，图5c为新得到的拟声波曲线。由此得到的拟声波曲线与密度乘积得到的新阻抗曲线对储层区分更好，叠置区域明显减少，如图6所示。
[0055]
s05：对拟声波曲线数据进行地震反演。
[0056]
具体地，利用步骤s04得到的拟声波曲线开展高分辨率地震反演，具体步骤为：标定精细储层；根据精细储层，约束拟声波曲线数据以构建地质框架模型；结合地质框架模型进行地震反演。
[0057]
首先，需要先标定精细储层，以便于后续较为准确的构建地质框架模型，标定精细储层具体包括：先选择较强的地震反射界面作为标准层；利用标准子波进行合成地震记录初始标定；根据大套地层对应的波组关系提取统计井旁道地震数据体中的子波，并进行二次标定；通过慢度曲线监控调整内部反射波组的对应关系。需要说明的是，标准子波可以为雷克子波、带通子波等；另外，在进行调整内部反射波组对应关系时，可利用井和地震重新提取子波以进行精细调整。
[0058]
其次，构建地质框架模型。在构建地质框架模型时，需要根据预先标定的精细储层，约束拟声波曲线数据在每个地层内进行内插，形成一个平滑、闭合的地质框架模型。
[0059]
最后，进行地震反演，在地震反演之前，还需要根据建立的地质框架模型，约束井阻抗的动态变化范围并进行内插外推，以确定约束条件，从而便于在地震反演时，根据约束条件，对地质框架模型进行稀疏约束脉冲反演。需要说明的是，通过在进行内插外推，可以
确定各个待反演地震道的各采样点处的阻抗动态范围，从而便于为反演提供准确的约束条件。
[0060]
在利用上述方法进行反演之后，得到新的声波阻抗反演结果，该反演结果能较好的区分储层与非储层，以便于对储层进行准确预测，参考图7，由图中可以看出，反演结果基本能区分储层与非储层，效果较好，最终也证明了该方法的有效性。
[0061]
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。