井震联合油藏储层预测方法、装置及电子设备

技术领域

本发明属于地震勘探资料解释技术领域，更具体地，涉及一种井震联合油藏储层预测方法、装置及电子设备。

背景技术

地震勘探资料解释工作是利用地震数据解释软件，以人机联作的方法进行地震资料解释。以偏移地震数据体为基础，技术人员通过认识和经验对工区内构造规律和波组特征的与解释软件相结合，全面开展地质构造、油气情况等解释工作。

储层的参数预测是寻找油气储量的基础工作之一，在油藏的开发和设计过程中占十分重要的地位。储层表征的目的是尽可能充分地把地质信息，如研究区域背景信息，地震数据，井数据和其他类型数据综合起来，并尽可能真实地刻画选定地区的地质变量的空间分布。其中，井间预测技术是石油地质统计学的核心部分。目前主要采用的是以井数据为基础来预测砂岩厚度，现有方法有反距离法，普通克里格法，指示克里格法等，也有储层表征的条件模拟技术，如序贯高斯模拟，序贯指示模拟等。

随着勘探开发的持续进行，井位的不断部署，工区的资料逐步丰富。而现有的普通克里格和序贯高斯模拟并没有完全利用丰富的资料，因此存在预测结果不够准确的问题。

发明内容

有鉴于此，本发明实施例提供了一种井震联合油藏储层预测方法、装置及电子设备，至少解决现有技术中预测结果不够准确的问题。

第一方面，本发明实施例提供了一种井震联合油藏储层预测方法，包括：

选取地震属性并对选取的地震属性进行标准化处理；

基于标准化处理的地震属性得到地震属性变异函数平面图；

基于所述地震属性变异函数平面图建立标准化处理的地震属性与油藏储层参数的关系；

基于所述关系对油藏储层进行预测。

可选的，所述选取地震属性并对选取的地震属性进行标准化处理，包括：

将获取的地震属性信息和井点信息结合；

基于对砂岩厚度敏感地震属性信息，建立地震属性信息与测井数据对应的关系，所述测井数据为经过测线的井点信息；

基于地震属性信息与测井数据对应的关系建立地震属性信息与测井数据交会图；

对于所述交会图进行相关性分析，选取对储层参数相关的地震属性信息。

可选的，所述选取地震属性并对选取的地震属性进行标准化处理，包括：

定义交点处的相交测线地震属性的偏差；

基于交点处地震属性均值得到相交测线地震属性与所述交点处地震属性均值的偏差；

基于所述地震属性的偏差和地震属性均值的偏差对地震属性进行标准化处理。

可选的，所述基于对砂岩厚度敏感地震属性信息中，所述对砂岩厚度敏感地震属性信息至少包括：平均瞬时频率和均方根振幅。

可选的，基于所述地震属性变异函数平面图建立标准化处理的地震属性与油藏储层参数的关系，包括：

基于地震属性中测线上得测点搜索附近的井点，建立测线属性与井数据的关系；

基于所述测线属性与井数据的关系得到测点的局部累积分布函数。

可选的，所述基于所述关系对油藏储层进行预测，包括：

指定网络单元模拟要访问的随机路径；

搜索网络单元中网格点的附近点，并求解克里格方程组，得到该网格点的局部累积分布函数；

基于所述局部累积分布函数对油藏储层参数进行预测。

可选的，所述搜索网络单元中网格点的附近点，并求解克里格方程组，得到该网格点的局部累积分布函数，包括：

求解克里格方程，得到每个测点对网格点的权重，从而得到该网格点的局部累积分布函数。

可选的，所述基于所述局部累积分布函数对油藏储层参数进行预测，包括：

从局部条件分布函数中随机抽取一个值作为该网格点的模拟值，并将得到的模拟值加入到原始井数据中，直到网格节点被全部模拟。

第二方面，本发明实施例还提供了一种井震联合油藏储层预测装置，包括：

选取模块，用于选取地震属性并对选取的地震属性进行标准化处理；

变异函数模块，用于基于标准化处理的地震属性得到地震属性变异函数平面图；

关系建立模块，用于基于所述地震属性变异函数平面图建立标准化处理的地震属性与油藏储层参数的关系；

预测模块，用于基于所述关系对油藏储层进行预测。

第三方面，本发明实施例还提供了一种电子设备，该电子设备包括：

存储器，存储有可执行指令；

处理器，所述处理器运行所述存储器中的所述可执行指令，以实现第一方面任一项所述的井震联合油藏储层预测方法。

本发明通过选取地震属性，并基于选取的地震属性得到地震属性变异函数平面图，从而建立标准化处理的地震属性与油藏储层参数的关系，基于建立的标准化处理的地震属性与油藏储层参数的关系实现对油藏储层的预测，将地震属性融入到储层参数的预测中，达到提高预测结果准确性的目的。

本发明的其它特征和优点将在随后具体实施方式部分予以详细说明。

附图说明

通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。

图1示出了本发明的一个实施例的砂岩厚度变异函数各向异性平面示意图；

图2示出了本发明的一个实施例的井震联合油藏储层预测方法的流程图；

图3a示出了本发明的一个实施例的标准化前地震属性示意图；

图3b示出了本发明的一个实施例的标准化后地震属性示意图；

图4a示出了本发明的一个实施例的标准化前测线交点地震属性交会图；

图4b示出了本发明的一个实施例的标准化后测线交点地震属性交会图；

图5示出了本发明的一个实施例的主次方向均方根振幅变异函数值示意图；

图6示出了本发明的一个实施例的均方根振幅的理论变异函数平面示意图；

图7示出了本发明的一个实施例的平均瞬时频率约束井震联合预测砂岩厚度示意图；

图8示出了本发明的一个实施例的均方根振幅约束约束井震联合预测砂岩厚度示意图；

图9示出了本发明的一个实施例的两种属性共同约束井震联合预测砂岩厚度示意图，两种属性即平均瞬时频率和均方根振幅。

具体实施方式

下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。

一种井震联合油藏储层预测方法，包括：

选取地震属性并对选取的地震属性进行标准化处理；

基于标准化处理的地震属性得到地震属性变异函数平面图；

基于所述地震属性变异函数平面图建立标准化处理的地震属性与油藏储层参数的关系；

基于所述关系对油藏储层进行预测。

可选的，所述选取地震属性并对选取的地震属性进行标准化处理，包括：

将获取的地震属性信息和井点信息结合；

基于对砂岩厚度敏感地震属性信息，建立地震属性信息与测井数据对应的关系，所述测井数据为经过测线的井点信息；

基于地震属性信息与测井数据对应的关系建立地震属性信息与测井数据交会图；

对于所述交会图进行相关性分析，选取对储层参数相关的地震属性信息。

可选的，所述选取地震属性并对选取的地震属性进行标准化处理，包括：

定义交点处的相交测线地震属性的偏差；

基于交点处地震属性均值得到相交测线地震属性与所述交点处地震属性均值的偏差；

基于所述地震属性的偏差和地震属性均值的偏差对地震属性进行标准化处理。

可选的，所述基于对砂岩厚度敏感地震属性信息中，所述对砂岩厚度敏感地震属性信息至少包括：平均瞬时频率和均方根振幅。

可选的，基于所述地震属性变异函数平面图建立标准化处理的地震属性与油藏储层参数的关系，包括：

基于地震属性中测线上得测点搜索附近的井点，建立测线属性与井数据的关系；

基于所述测线属性与井数据的关系得到测点的局部累积分布函数。

可选的，所述基于所述关系对油藏储层进行预测，包括：

指定网络单元模拟要访问的随机路径；

搜索网络单元中网格点的附近点，并求解克里格方程组，得到该网格点的局部累积分布函数；

基于所述局部累积分布函数对油藏储层参数进行预测。

可选的，所述搜索网络单元中网格点的附近点，并求解克里格方程组，得到该网格点的局部累积分布函数，包括：

求解克里格方程，得到每个测点对网格点的权重，从而得到该网格点的局部累积分布函数。

可选的，所述基于所述局部累积分布函数对油藏储层参数进行预测，包括：

从局部条件分布函数中随机抽取一个值作为该网格点的模拟值，并将得到的模拟值加入到原始井数据中，直到网格节点被全部模拟。

实施例一：

本实施例多源信息空间结构特征分析技术在地震解释中应用的基本思路：

基于克里格法，结合序贯思想和随机建模思想，提出井震联合预测油藏储层参数。这种方法的主要思想是将地震信息按一定的权重，以局部累积分布函数的方式融入到储层参数的预测中，并在一定风险下给出研究区域的储层参数的预测结果，而各个模型结果之间的差别就反映了由于缺乏足够的资料等原因而引起不确定性。

(1)地震属性选取及标准化：

将地震属性信息和井点信息结合，即把地震属性与砂岩厚度之间的相关性、地震属性和砂岩厚度之间的空间结构特征结合。标准化的目的是为了消除其量纲、大小和变化范围等方面的差异性。

假设地震测线总数为L条，其在空间上有N个交点，每一个测线交点都对应有相交的两条测线的地震属性，其中的第i个交点可表示为

该定义的偏差值能反应标准化后地震属性的优劣，如果偏差值越小，那么说明测线之间地震属性的可比性就越好。对标准化参数

对于交点

/>

根据计算得到的偏差对相交测线的标准化处理的地震属性进行调整：

其中，α和β为调整系数，在计算过程中可以随机调整该参数的幅度。重复上述过程进行迭代计算，可以得到最终的标准化处理参数

对砂岩厚度比较敏感的地震属性有平均瞬时频率(FreQ)和均方根振幅(RMS)等。对于经过测线的井点数据，那么直接可以通过建立并地震属性值与测井数据对应的关系。通过建立地震属性与测井数据交会图，进行相关性分析，即可选取对储层参数影响较大的地震属性。

(2)地震属性变异函数平面图

分析储层参数的实验变异函数平面图，可以找到储层参数变异的主方向和次方向，地震属性与储层参数的变异方向是一致的，但是变程有所不同，可通过地震测线的变差函数图寻找该地震属性的主方向和次方向的变程。然后进行坐标变换，就可以得到地震数据的变异函数平面图，如图1所示。

u方向变程最大，令其为a

(3)地震属性与储层参数关系建立

地震属性与井数据的关系可用概率来表现。以均方根振幅为例，对于未经过测线的井点数据，测线上的测点搜索最近的井点，找到其对应的储层参数值，若超过一定的距离没有找到井点，则寻找下一个侧线上的测点。那么，就可以得到均方根振幅-砂岩厚度的分布图。对于一个均方根振幅小区间，都存在多个砂岩厚度的井点数。于是，可以计算出每个均方根振幅的区间对应的储层参数的条件累积分布函数。当某测线上的测点的地震属性落在某个区间，就可知道该测点的条件累积分布函数。

(4)井震联合预测

随机产生一个经过所有网格节点的随机路径，对于待估计的网格节点，搜索附近的地震测线的测点，每个测点根据地震属性和测井数据的关系，可以得到每个测点的局部累积分布函数，通过求解克里格方程，就可以知道每个测点对待估点的权重。那么，可以计算出待估网格节点的局部累积分布，从局部条件分布函数中随机抽取一个值作为该点的模拟值，并将这个新模拟的值加入到原始井数据中。重复上述步骤，直到网格节点都被模拟完为止，便完成了对研究区井震联合模拟的一个实现。当然，随机路径不同，产生等概率的实现大致趋势相同，局部细节可能存在差异。

具体步骤如图2所示。其实施过程中的样如图3a至图9所示。

实施例二：

一种井震联合油藏储层预测装置，包括：

选取模块，用于选取地震属性并对选取的地震属性进行标准化处理；

变异函数模块，用于基于标准化处理的地震属性得到地震属性变异函数平面图；

关系建立模块，用于基于所述地震属性变异函数平面图建立标准化处理的地震属性与油藏储层参数的关系；

预测模块，用于基于所述关系对油藏储层进行预测。

可选的，所述选取地震属性并对选取的地震属性进行标准化处理，包括：

将获取的地震属性信息和井点信息结合；

基于对砂岩厚度敏感地震属性信息，建立地震属性信息与测井数据对应的关系，所述测井数据为经过测线的井点信息；

基于地震属性信息与测井数据对应的关系建立地震属性信息与测井数据交会图；

对于所述交会图进行相关性分析，选取对储层参数相关的地震属性信息。

可选的，所述选取地震属性并对选取的地震属性进行标准化处理，包括：

定义交点处的相交测线地震属性的偏差；

基于交点处地震属性均值得到相交测线地震属性与所述交点处地震属性均值的偏差；

基于所述地震属性的偏差和地震属性均值的偏差对地震属性进行标准化处理。

可选的，所述基于对砂岩厚度敏感地震属性信息中，所述对砂岩厚度敏感地震属性信息至少包括：平均瞬时频率和均方根振幅。

可选的，基于所述地震属性变异函数平面图建立标准化处理的地震属性与油藏储层参数的关系，包括：

基于地震属性中测线上得测点搜索附近的井点，建立测线属性与井数据的关系；

基于所述测线属性与井数据的关系得到测点的局部累积分布函数。

可选的，所述基于所述关系对油藏储层进行预测，包括：

指定网络单元模拟要访问的随机路径；

搜索网络单元中网格点的附近点，并求解克里格方程组，得到该网格点的局部累积分布函数；

基于所述局部累积分布函数对油藏储层参数进行预测。

可选的，所述搜索网络单元中网格点的附近点，并求解克里格方程组，得到该网格点的局部累积分布函数，包括：

求解克里格方程，得到每个测点对网格点的权重，从而得到该网格点的局部累积分布函数。

可选的，所述基于所述局部累积分布函数对油藏储层参数进行预测，包括：

从局部条件分布函数中随机抽取一个值作为该网格点的模拟值，并将得到的模拟值加入到原始井数据中，直到网格节点被全部模拟。

实施例三：

本发明实施例提供一种电子设备包括存储器和处理器,

存储器，存储有可执行指令；

处理器，处理器运行存储器中的可执行指令，以实现井震联合油藏储层预测方法。

该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(RAM)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(ROM)、硬盘、闪存等。

该处理器可以是中央处理单元(CPU)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本发明的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。

本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本发明的保护范围之内。

有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。

实施例四：

本发明实施例提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现井震联合油藏储层预测方法。

根据本发明实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本发明各实施例方法的全部或部分步骤。

上述计算机可读存储介质包括但不限于：光存储介质(例如：CD－ROM和DVD)、磁光存储介质(例如：MO)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置ROM的媒体(例如：ROM盒)。

以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。