一种储层特征分析方法、存储介质及电子设备与流程

1.本发明涉及地震勘探技术领域，特别地涉及一种基于多属性融合的储层特征分析方法、存储介质及电子设备。


背景技术：

2.随着地震勘探技术的进步和地震属性分析技术的迅速发展，地震属性在石油勘探与开发各个环节中发挥的作用逐渐彰显。地震属性是指从地震数据中导出的关于几何学、运动学、动力学及统计特性的特殊度量，是反映地震波的特征参数。每一种地震属性只是储层一个或几个特征参数的地球物理响应，并不能反映储层整个性质特征，一种地震属性只适用于储层某一特征的预测，比如，振幅属性能较好地识别岩性、溶洞、不整合等现象，相位属性在刻画地质体边界方面效果较好，频率属性对储层厚度、含油气性较为敏感，因此说明了任何一种地震属性预测都具有特殊性与针对性，从另一方面也揭示了单一的地震属性预测具有一定的局限性。
3.为了较好地解决单一属性的局限性问题，最大程度地克服使用单一属性的盲目性，研究人员通常采用多属性分析与融合技术，即从繁多的地震属性中优选出有效的属性信息，同时参与储层预测研究，充分运用地震数据具有较好横向连续性的特征，优选出与储层特征有关的优化地震属性及属性组合，刻画储层的平面展布特征，有效提高地震储层预测精度。
4.目前的三维体地震多属性分析与融合技术，大多基于经验系数和机器学习算法，经验系数在新的探区常常需要摸索较长的时间，耗时耗力，而基于神经网络、支持向量机等的机器学习算法由于缺乏地质约束和先验信息的控制，预测结果常常具有较多的多解性和不确定性。
5.因此，本发明基于上述问题，提供一种基于多属性融合的储层特征分析方法、存储介质及电子设备。


技术实现要素：

6.本发明的目的在于：针对上述问题，本发明提供了一种基于多属性融合的储层特征分析方法、存储介质及电子设备，解决了现有的三维体地震多属性分析方法存在耗时耗力、预测结果具有多解性和不确定性的问题，达到可以兼顾储层沿层切片的地质规律，又能在数据驱动的基础上迅速找到最优解的目的。
7.本发明采用的技术方案如下：
8.为实现上述目的，第一方面，本发明提供一种储层特征分析方法，所述方法包括：
9.获取目标储层的若干个三维地震属性数据体，以及每个三维地震属性数据体所对应的沿层时间切片；
10.对于目标储层的每个层位，利用沿层时间切片方程将所有沿层时间切片融合为一个综合沿层切片，以此获得目标储层多个层位的多个综合沿层切片；
11.利用最小二乘法对所述多个综合沿层切片进行反演，得到若干个最优权重系数，利用所述若干个最优权重系数对所述若干个三维地震属性数据体进行加权处理得到三维地震属性融合体；
12.利用三维地震属性融合体对目标储层的整体特征进行分析。
13.根据本发明的实施例，可选的，上述储层特征分析方法中，获取目标储层的若干个三维地震属性数据体，以及每个三维地震属性数据体所对应的沿层时间切片包括：
14.根据地质约束确定目标储层的n个地震属性；
15.利用叠前地震道集和/或叠后地震属性获取与n个地震属性对应的n个三维地震属性数据体p1,p2,
…
,pn；
16.对所述n个三维地震属性数据体p1,p2,
…
,pn分别进行反演，以得到与所述n个三维地震属性数据体p1,p2,
…
,pn对应的n个沿层时间切片。
17.根据本发明的实施例，可选的，上述储层特征分析方法中，所述地质约束包括目标储层的测井数据和地质认识分析，所述地质认识分析包括目标储层的地质成因分析。
18.根据本发明的实施例，可选的，上述储层特征分析方法中，所述沿层时间切片方程为：
19.m＝ω1s1 ω2s2
…
ωnsn，
20.其中，n表示地震属性的个数，ω1,ω2,
…
,ωn表示n个地震属性在沿层时间切片中的权重系数，s1,s2,
…
,sn表示沿层时间切片数据，m表示综合沿层切片。
21.根据本发明的实施例，可选的，上述储层特征分析方法中，在利用所述沿层时间切片方程将所有沿层时间切片融合为一个综合沿层切片时，所述n个地震属性在沿层时间切片中的权重系数ω1,ω2,
…
,ωn根据不同地震属性对于目标储层发育的贡献值进行调整。
22.根据本发明的实施例，可选的，上述储层特征分析方法中，利用下式获得目标储层多个层位的多个综合沿层切片：
[0023][0024]
其中，i表示目标储层的层位数，ω
i,n
表示第n个地震属性在第i层地层的第n个沿层时间切片中的权重系数，mi表示第i层地层的综合沿层切片。
[0025]
根据本发明的实施例，可选的，上述储层特征分析方法中，利用最小二乘法对所述多个综合沿层切片进行反演得到若干个最优权重系数时的计算式为：
[0026][0027]
其中，n表示地震属性的个数，表示n个地震属性在三维地震属性融合体中的最优权重系数，i表示目标储层的层位数，ω
i,n
表示第n个地震属性在第i层地层的第n个沿层时间切片中的权重系数。
[0028]
根据本发明的实施例，可选的，上述储层特征分析方法中，利用所述若干个最优权重系数对所述若干个三维地震属性数据体进行加权处理得到三维地震属性融合体时的计算式为：
[0029][0030]
其中，z表示三维地震属性融合体，p1,p2,
…
,pn表示n个三维地震属性数据体。
[0031]
第二方面，本发明提供了一种存储介质，其上存储有计算机程序，所述计算机程序可被一个或多个处理器执行，以实现如上述的储层特征分析方法。
[0032]
第三方面，本发明提供了一种电子设备，所述电子设备包括存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，执行如上述的储层特征分析方法。
[0033]
与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：
[0034]
1.本发明提供的一种储层特征分析方法、存储介质及电子设备，通过对于目标储层的每个层位，利用沿层时间切片方程将获取的所有沿层时间切片融合为多个综合沿层切片，再利用最小二乘法对多个综合沿层切片进行反演得到最优权重系数，再结合获取的若干个三维地震属性数据体进行加权处理得到三维地震属性融合体，最后再利用三维地震属性融合体对目标储层的整体特征进行分析；通过对目标储层多个沿层切片的地质趋势和目标储层的地层趋势进行分析和控制，开展基于最小二乘算法的三维地震属性融合体的参数求取，针对地震预测中的多解性提供了确定解和唯一解，有效提高地震储层特征分析的精度和准确度，从而提高地震勘探的效率。
[0035]
2.本发明得到的三维地震属性融合体可以兼顾地层沿层切片中反映的地质规律和预测目标，又能在地震数据的驱动下形成三维空间融合体数据，迅速找到最优解，为后续的油藏地质建模、油藏精细开发等研究提供快捷可靠的三维地震成果数据体。
[0036]
3.本发明利用最小二乘法对所述多个综合沿层切片进行反演，可以简便地求得得到若干个最优权重系数，充分考虑到了不同地层之间的相关性，防止反演结果陷入局部最优解，具有更高的稳健性和可靠性。
附图说明
[0037]
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
[0038]
在下文中将基于实施例并参考附图来对本发明进行更详细的描述。
[0039]
图1为本发明实施例二提供的一种储层特征分析方法的流程示意图。
[0040]
图2为本发明实施例三提供的一种储层特征分析方法中horizon-1层的5个沿层时间切片图。
[0041]
图3为本发明实施例三提供的一种储层特征分析方法中horizon-1层的综合沿层切片图。
[0042]
图4为本发明实施例三提供的一种储层特征分析方法中horizon-2层的综合沿层
切片图。
[0043]
图5为本发明实施例三提供的一种储层特征分析方法中horizon-3层的综合沿层切片图。
[0044]
图6为本发明实施例三提供的一种储层特征分析方法中horizon-4层的综合沿层切片图。
[0045]
图7为本发明实施例三提供的一种储层特征分析方法得到的三维地震属性融合体。
[0046]
在附图中，相同的部件使用相同的附图标记，附图并未按照实际的比例绘制。
具体实施方式
[0047]
以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达到相应技术效果的实现过程能充分理解并据以实施。本发明实施例以及实施例中的各个特征，在不相冲突前提下可以相互结合，所形成的技术方案均在本发明的保护范围之内。
[0048]
实施例一
[0049]
本发明实施例提供了一种可应用于电子设备的基于多属性融合的储层特征分析方法，所述方法应用于所述电子设备时执行以下步骤：
[0050]
步骤1：获取目标储层的若干个三维地震属性数据体，以及每个三维地震属性数据体所对应的沿层时间切片；
[0051]
步骤1.1：确定目标储层的n个地震属性；
[0052]
步骤1.2：利用叠前地震道集获取与n个地震属性对应的n个三维地震属性数据体p1,p2,
…
,pn；
[0053]
步骤1.3：对所述n个三维地震属性数据体p1,p2,
…
,pn分别进行反演，以得到与所述n个三维地震属性数据体p1,p2,
…
,pn对应的n个沿层时间切片；
[0054]
步骤2：对于目标储层的每个层位，利用沿层时间切片方程将所有沿层时间切片融合为一个综合沿层切片，以此获得目标储层多个层位的多个综合沿层切片；
[0055]
步骤2.1：对于目标储层的每个层位，利用沿层时间切片方程将该层位的n个沿层时间切片融合为一个综合沿层切片m，所述沿层时间切片方程为：
[0056]
m＝ω1s1 ω2s2
…
ωnsn，
[0057]
其中，n表示地震属性的个数，ω1,ω2,
…
,ωn表示n个地震属性在沿层时间切片中的权重系数，s1,s2,
…
,sn表示沿层时间切片数据，m表示综合沿层切片；
[0058]
步骤2.2：根据步骤2.1利用下式得到目标储层的i个层位的i个综合沿层切片m1,m2,
…
,mi：
[0059][0060]
其中，i表示目标储层发育区的层位数，ω
i,n
表示第n个地震属性在第i层地层的第n个沿层时间切片中的权重系数，mi表示第i层地层的综合沿层切片；
[0061]
步骤3：利用最小二乘法对所述多个综合沿层切片进行反演，得到若干个最优权重系数，利用所述若干个最优权重系数对所述若干个三维地震属性数据体进行加权处理得到三维地震属性融合体；
[0062]
步骤3.1：利用最小二乘法对i个综合沿层切片m1,m2,
…
,mi进行反演得到n个最优权重系数：
[0063][0064]
其中，n表示地震属性的个数，表示n个地震属性在三维地震属性融合体中的最优权重系数，i表示目标储层发育区的层位数，ω
i,n
表示第n个地震属性在第i层地层的第n个沿层时间切片中的权重系数，采用最小二乘法充分考虑到了不同地层之间的相关性，防止反演结果陷入局部最优解，具有更高的稳健性和可靠性；
[0065]
步骤3.2：利用所述步骤3.1得到的n个最优权重系数和所述步骤1.2得到的n个三维地震属性数据体p1,p2,
…
,pn进行加权处理得到三维地震属性融合体z：
[0066][0067]
步骤4：利用三维地震属性融合体对目标储层的整体特征进行分析。
[0068]
本实施例提供的基于多属性融合的储层特征分析方法，通过对于目标储层的每个层位，利用沿层时间切片方程将获取的所有沿层时间切片融合为多个综合沿层切片，再利用最小二乘法对多个综合沿层切片进行反演得到最优权重系数，再结合获取的若干个三维地震属性数据体进行加权处理得到三维地震属性融合体，最后再利用三维地震属性融合体对目标储层的整体特征进行分析；通过对目标储层多个沿层切片的地质趋势和目标储层的地层趋势进行分析和控制，开展基于最小二乘算法的三维地震属性融合体的参数求取，针对地震预测中的多解性提供了确定解和唯一解，有效提高地震储层特征分析的精度和准确度。
[0069]
实施例二
[0070]
请参阅图1，本发明实施例在实施例一的基础上提供了一种可应用于电子设备的基于多属性融合的储层特征分析方法。
[0071]
进一步地，获取目标储层的若干个三维地震属性数据体，以及每个三维地震属性数据体所对应的沿层时间切片包括：
[0072]
根据地质约束确定目标储层的n个地震属性；
[0073]
利用叠后地震属性获取与n个地震属性对应的n个三维地震属性数据体p1,p2,
…
,pn；
[0074]
对所述n个三维地震属性数据体p1,p2,
…
,pn分别进行反演，以得到与所述n个三维地震属性数据体p1,p2,
…
,pn对应的n个沿层时间切片；
[0075]
其中，所述地质约束包括目标储层的测井数据和地质认识分析，所述地质认识分析包括目标储层的地质成因分析。
[0076]
进一步地，在利用所述沿层时间切片方程将所有沿层时间切片融合为一个综合沿
层切片时，所述n个地震属性在沿层时间切片中的权重系数ω1,ω2,
…
,ωn根据不同地震属性对于目标储层发育的贡献值进行调整。
[0077]
具体地，所述方法应用于所述电子设备时具体执行以下步骤：
[0078]
步骤1：获取目标储层的若干个三维地震属性数据体，以及每个三维地震属性数据体所对应的沿层时间切片；
[0079]
步骤1.1：根据地质约束确定目标储层的n个地震属性；
[0080]
其中，所述地质约束包括目标储层的测井数据和地质认识分析，所述地质认识分析包括目标储层的地质成因分析；
[0081]
步骤1.2：利用叠后地震属性获取与n个地震属性对应的n个三维地震属性数据体p1,p2,
…
,pn；
[0082]
步骤1.3：对所述n个三维地震属性数据体p1,p2,
…
,pn分别进行反演，以得到与所述n个三维地震属性数据体p1,p2,
…
,pn对应的n个沿层时间切片；
[0083]
具体地，结合地质约束沿层向下开时窗，以开展层位属性面的提取，获得对应的n个沿层时间切片；
[0084]
步骤2：对于目标储层的每个层位，利用沿层时间切片方程将所有沿层时间切片融合为一个综合沿层切片，以此获得目标储层多个层位的多个综合沿层切片；
[0085]
步骤2.1：对于目标储层的每个层位，利用沿层时间切片方程将该层位的n个沿层时间切片融合为一个综合沿层切片m，所述沿层时间切片方程为：
[0086]
m＝ω1s1 ω2s2
…
ωnsn，
[0087]
其中，n表示地震属性的个数，ω1,ω2,
…
,ωn表示n个地震属性在沿层时间切片中的权重系数，s1,s2,
…
,sn表示沿层时间切片数据，m表示综合沿层切片；
[0088]
在上述过程中，所述n个地震属性在沿层时间切片中的权重系数ω1,ω2,
…
,ωn根据不同地震属性对于目标储层发育的贡献值进行调整；
[0089]
具体地，以所述n个沿层时间切片为输入，输入到所述沿层时间切片方程中，并结合地质约束和这n个地震属性对于储层发育的贡献值不断调整输入的这n个沿层时间切片的权重系数，调整到最能反映储层特征的n个沿层时间切片的权重系数后，输出得到既能反映地层由于区域应力或者断层破裂引起的大、中、小尺度的断裂-裂缝，又能反映其中由于岩性等引起的裂缝的综合沿层切片m；
[0090]
步骤2.2：根据步骤2.1利用下式得到目标储层的i个层位的i个综合沿层切片m1,m2,
…
,mi：
[0091][0092]
其中，i表示目标储层发育区的层位数，ω
i,n
表示第n个地震属性在第i层地层的第n个沿层时间切片中的权重系数，mi表示第i层地层的综合沿层切片；
[0093]
步骤3：利用最小二乘法对所述多个综合沿层切片进行反演，得到若干个最优权重系数，利用所述若干个最优权重系数对所述若干个三维地震属性数据体进行加权处理得到三维地震属性融合体；
[0094]
步骤3.1：利用最小二乘法对i个综合沿层切片m1,m2,
…
,mi进行反演得到n个最优权重系数：
[0095][0096]
其中，n表示地震属性的个数，表示n个地震属性在三维地震属性融合体中的最优权重系数，i表示目标储层发育区的层位数，ω
i,n
表示第n个地震属性在第i层地层的第n个沿层时间切片中的权重系数，采用最小二乘法充分考虑到了不同地层之间的相关性，防止反演结果陷入局部最优解，具有更高的稳健性和可靠性；
[0097]
步骤3.2：利用所述步骤3.1得到的n个最优权重系数和所述步骤1.2得到的n个三维地震属性数据体p1,p2,
…
,pn进行加权处理得到三维地震属性融合体z：
[0098][0099]
步骤4：利用三维地震属性融合体对目标储层的整体特征进行分析。
[0100]
本实施例提供的基于多属性融合的储层特征分析方法，通过对于目标储层的每个层位，利用沿层时间切片方程将获取的所有沿层时间切片融合为多个综合沿层切片，再利用最小二乘法对多个综合沿层切片进行反演得到最优权重系数，再结合获取的若干个三维地震属性数据体进行加权处理得到三维地震属性融合体，最后再利用三维地震属性融合体对目标储层的整体特征进行分析；得到的三维地震属性融合体可以兼顾地层沿层切片中反映的地质规律和预测目标，又能在地震数据的驱动下形成三维空间融合体数据，迅速找到最优解，为后续的油藏地质建模、油藏精细开发等研究提供快捷可靠的三维地震成果数据体。
[0101]
实施例三
[0102]
请参阅图1至图7，本发明最佳实施例基于实施例二提供了一种在地震三维解释中的属性降维技术，具体针对工区的储层类型，提供了一种可应用于电子设备的基于多属性融合的储层特征分析方法，如图1所示，所述方法应用于所述电子设备时执行以下步骤：
[0103]
步骤1：获取目标储层的若干个三维地震属性数据体，以及每个三维地震属性数据体所对应的沿层时间切片；
[0104]
步骤1.1：根据地质约束确定目标储层的n个地震属性；
[0105]
其中，所述地质约束包括目标储层的测井数据和地质认识分析，所述地质认识分析包括目标储层的地质成因分析；
[0106]
本实施例中，以实际中的工区为例，发育的为裂缝性储层，裂缝性储层是指以裂缝为主要储集空间、渗流通道的储集层，根据该裂缝性储层的测井数据和裂缝的地质成因分析可知，倾角属性、泥质含量、afe属性、叠前各向异性属性和曲率属性这5种地震属性对最终的裂缝发育都有不同程度的贡献，因此确定本实施例中的地震属性包括倾角属性、泥质含量、afe属性、叠前各向异性裂缝发育强度体属性和曲率属性，n＝5；
[0107]
步骤1.2：利用叠后地震属性获取与n个地震属性对应的n个三维地震属性数据体p1,p2,
…
,pn；
[0108]
本实施例中，利用叠后地震属性的沿层切片数据获取到5个三维地震属性数据体p1,p2,p3,p4,p5；
[0109]
步骤1.3：对所述n个三维地震属性数据体p1,p2,
…
,pn分别进行反演，以得到与所述n个三维地震属性数据体p1,p2,
…
,pn对应的n个沿层时间切片；
[0110]
具体地，结合地质约束沿层向下开时窗，以开展层位属性面的提取，获得对应的n个沿层时间切片；
[0111]
本实施例中，以地层horizon-1为例，结合前期的测井数据和裂缝的地质成因分析成果，沿层向下开12ms的时窗获得了如图2所示的5个沿层时间切片，由图可知，5个沿层时间切片的整体趋势一致，且对储层最终的发育都有一定程度的贡献，这与测井数据和裂缝的地质成因分析是相符合的；
[0112]
步骤2：对于目标储层的每个层位，利用沿层时间切片方程将所有沿层时间切片融合为一个综合沿层切片，以此获得目标储层多个层位的多个综合沿层切片；
[0113]
步骤2.1：对于目标储层的每个层位，利用沿层时间切片方程将该层位的n个沿层时间切片融合为一个综合沿层切片m，所述沿层时间切片方程为：
[0114]
m＝ω1s1 ω2s2
…
ωnsn，
[0115]
其中，n表示地震属性的个数，ω1,ω2,
…
,ωn表示n个地震属性在沿层时间切片中的权重系数，s1,s2,
…
,sn表示沿层时间切片数据，m表示综合沿层切片；
[0116]
在上述过程中，所述n个地震属性在沿层时间切片中的权重系数ω1,ω2,
…
,ωn根据不同地震属性对于目标储层发育的贡献值进行调整；
[0117]
具体地，以所述n个沿层时间切片为输入，输入到所述沿层时间切片方程中，并结合地质约束和这n个地震属性对于储层发育的贡献值不断调整输入的这n个沿层时间切片的权重系数，调整到最能反映储层特征的n个沿层时间切片的权重系数后，输出得到既能反映地层由于区域应力或者断层破裂引起的大、中、小尺度的断裂-裂缝，又能反映其中由于岩性等引起的裂缝的综合沿层切片m；
[0118]
本实施例中，5个地震属性对于储层发育的贡献值也有不同，例如，afe属性反映的是大尺度的断裂-裂缝特征，叠前各向异性裂缝发育强度体属性反映的是小尺度各向异性裂缝发育带，曲率属性反映的是地震三维体沿层弯曲程序的裂缝发育；
[0119]
本实施例中，先融合该裂缝性储层的horizon-1层，以如图2所示的horizon-1层的5个沿层时间切片为输入，输入到沿层时间切片方程：
[0120]
m1＝ω1s1 ω2s2 ω3s3 ω4s4 ω5s5，
[0121]
得到如图3所示的horizon-1层的综合沿层切片m1；
[0122]
步骤2.2：根据所述步骤2.1利用下式得到目标储层的i个层位的i个综合沿层切片m1,m2,
…
,mi：
[0123][0124]
其中，i表示目标储层发育区的层位数，ω
i,n
表示第n个地震属性在第i层地层的第n个沿层时间切片中的权重系数，mi表示第i层地层的综合沿层切片；
[0125]
本实施例中，整个裂缝性储层发育区为4个层位horizon-1、horizon-2、horizon-3、horizon-4，因此，根据步骤2.1利用下式得到了如图3至6所示的4个综合沿层切片m1,m2,m3,m4，具体计算式为：
[0126][0127]
步骤3：利用最小二乘法对所述多个综合沿层切片进行反演，得到若干个最优权重系数，利用所述若干个最优权重系数对所述若干个三维地震属性数据体进行加权处理得到三维地震属性融合体；
[0128]
步骤3.1：利用最小二乘法对i个综合沿层切片m1,m2,
…
,mi进行反演得到n个最优权重系数：
[0129][0130]
其中，n表示地震属性的个数，表示n个地震属性在三维地震属性融合体中的最优权重系数，i表示目标储层发育区的层位数，ω
i,n
表示第n个地震属性在第i层地层的第n个沿层时间切片中的权重系数；最小二乘法又称最小平方法，利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小，采用最小二乘法充分考虑到了不同地层之间的相关性，防止反演结果陷入局部最优解，具有更高的稳健性和可靠性；
[0131]
本实施例中，利用最小二乘法对4个综合沿层切片m1,m2,m3,m4进行反演，得到5个最优权重系数，具体计算式为：
[0132][0133]
使用for循环求解上式，得到5个最优权重系数充分考虑到了不同地层之间的相关性，防止反演结果陷入局部最优解，因此采用最小二乘算法，具有更高的稳健性和可靠性；
[0134]
步骤3.2：利用所述步骤3.1得到的n个最优权重系数和所述步骤1.2得到的n个三维地震属性数据体p1,p2,
…
,pn进行加权处理，得到三维地震属性融合体z：
[0135][0136]
本实施例中，利用所述步骤3.1得到的5个最优权重系数以及所
述步骤1.1得到的5个三维地震属性数据体p1,p2,p3,p4,p5进行加权处理，得到三维地震属性融合体z：
[0137][0138]
计算上式，得到如图7所示的三维地震属性融合体z；
[0139]
步骤4：利用三维地震属性融合体z对目标储层的整体特征进行分析；
[0140]
本实施例中得到的三维地震属性融合体z综合考虑了该裂缝性储层沿层的地质趋势和测井分析结果，兼顾了各个地震属性对于储层的物理意义和地质意义，针对地震预测中的多解性提供了确定解和唯一解，有效提高地震储层特征分析的精度和准确度，从而提高地震勘探的效率。该三维地震属性融合体z对于后续的油藏地质建模、油藏精细开发、新探井的勘探开发等研究提供了一种快捷可靠的三维地震成果数据体。
[0141]
本实施例提供的基于多属性融合的储层特征分析方法，通过对于目标储层的每个层位，利用沿层时间切片方程将获取的所有沿层时间切片融合为多个综合沿层切片，再利用最小二乘法对多个综合沿层切片进行反演得到最优权重系数，再结合获取的若干个三维地震属性数据体进行加权处理得到三维地震属性融合体，最后再利用三维地震属性融合体对目标储层的整体特征进行分析；通过对目标储层多个沿层切片的地质趋势和目标储层的地层趋势进行分析和控制，开展基于最小二乘算法的三维地震属性融合体的参数求取，针对地震预测中的多解性提供了确定解和唯一解，有效提高地震储层特征分析的精度和准确度，从而提高地震勘探的效率；得到的三维地震属性融合体可以兼顾地层沿层切片中反映的地质规律和预测目标，又能在地震数据的驱动下形成三维空间融合体数据，迅速找到最优解，为后续的油藏地质建模、油藏精细开发等研究提供快捷可靠的三维地震成果数据体；能够解决现有的三维体地震多属性分析方法存在耗时耗力、预测结果具有多解性和不确定性的问题，达到可以兼顾储层沿层切片的地质规律，又能在数据驱动的基础上迅速找到最优解的目的；利用最小二乘法对所述多个综合沿层切片进行反演，可以简便地求得得到若干个最优权重系数，充分考虑到了不同地层之间的相关性，防止反演结果陷入局部最优解，具有更高的稳健性和可靠性。
[0142]
实施例四
[0143]
本发明实施例提供一种计算机可读存储介质，如闪存、硬盘、多媒体卡、卡型存储器(例如，sd或dx存储器等)、随机访问存储器(ram)、静态随机访问存储器(sram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、可编程只读存储器(prom)、磁性存储器、磁盘、光盘、服务器、app应用商城等等，其上存储有计算机程序，所述计算机程序被处理器执行时可以实现如下方法步骤：
[0144]
步骤1：获取目标储层的若干个三维地震属性数据体，以及每个三维地震属性数据体所对应的沿层时间切片；
[0145]
步骤1.1：根据地质约束确定目标储层的n个地震属性；
[0146]
其中，所述地质约束包括目标储层的测井数据和地质认识分析，所述地质认识分析包括目标储层的地质成因分析；
[0147]
步骤1.2：利用叠前地震道集或叠后地震属性获取与n个地震属性对应的n个三维地震属性数据体p1,p2,
…
,pn；
[0148]
步骤1.3：对所述n个三维地震属性数据体p1,p2,
…
,pn分别进行反演，以得到与所
述n个三维地震属性数据体p1,p2,
…
,pn对应的n个沿层时间切片；
[0149]
具体地，结合地质约束沿层向下开时窗，以开展层位属性面的提取，获得对应的n个沿层时间切片；
[0150]
步骤2：对于目标储层的每个层位，利用沿层时间切片方程将所有沿层时间切片融合为一个综合沿层切片，以此获得目标储层多个层位的多个综合沿层切片；
[0151]
步骤2.1：对于目标储层的每个层位，利用沿层时间切片方程将所述n个沿层时间切片融合为一个综合沿层切片m，所述沿层时间切片方程为：
[0152]
m＝ω1s1 ω2s2
…
ωnsn，
[0153]
其中，n表示地震属性的个数，ω1,ω2,
…
,ωn表示n个地震属性在沿层时间切片中的权重系数，s1,s2,
…
,sn表示沿层时间切片数据，m表示综合沿层切片；
[0154]
在上述过程中，所述n个地震属性在沿层时间切片中的权重系数ω1,ω2,
…
,ωn根据不同地震属性对于目标储层发育的贡献值进行调整；
[0155]
具体地，以所述n个沿层时间切片为输入，输入到所述沿层时间切片方程中，并结合地质约束和这n个地震属性对于储层发育的贡献值不断调整输入的这n个沿层时间切片的权重系数，调整到最能反映储层特征的n个沿层时间切片的权重系数后，输出得到既能反映地层由于区域应力或者断层破裂引起的大、中、小尺度的断裂-裂缝，又能反映其中由于岩性等引起的裂缝的综合沿层切片m；
[0156]
步骤2.2：根据步骤2.1利用下式得到目标储层的i个层位的i个综合沿层切片m1,m2,
…
,mi：
[0157][0158]
其中，i表示目标储层发育区的层位数，ω
i,n
表示第n个地震属性在第i层地层的第n个沿层时间切片中的权重系数，mi表示第i层地层的综合沿层切片；
[0159]
步骤3：利用最小二乘法对所述多个综合沿层切片进行反演，得到若干个最优权重系数，利用所述若干个最优权重系数对所述若干个三维地震属性数据体进行加权处理得到三维地震属性融合体；
[0160]
步骤3.1：利用最小二乘法对i个综合沿层切片m1,m2,
…
,mi进行反演得到n个最优权重系数：
[0161][0162]
其中，n表示地震属性的个数，表示n个地震属性在三维地震属性融合体中的最优权重系数，i表示目标储层发育区的层位数，ω
i,n
表示第n个地震属性在第i层地层的第n个沿层时间切片中的权重系数；
[0163]
使用for循环求解上式，得到最优权重系数
[0164]
步骤3.2：利用所述步骤3.1得到的n个最优权重系数和所述步骤1.2得到的n个三维地震属性数据体p1,p2,
…
,pn进行加权处理得到三维地震属性融合体z：
[0165][0166]
步骤4：利用三维地震属性融合体z对目标储层的整体特征进行分析。
[0167]
上述方法步骤的具体实施例过程可参见实施例三，本实施例在此不再重复赘述。
[0168]
实施例五
[0169]
本实施例提供了一种电子设备，该电子设备可以是手机、电脑或平板电脑等，包括存储器和处理器，所述存储器上存储有计算器程序，该计算机程序被处理器执行时实现一种基于多属性融合的储层特征分析方法。可以理解，电子设备还可以包括，多媒体组件，输入/输出(i/o)接口，以及通信组件。
[0170]
其中，处理器用于执行一种基于多属性融合的储层特征分析方法中的全部或部分步骤。存储器用于存储各种类型的数据，这些数据例如可以包括电子设备中的任何应用程序或方法的指令，以及应用程序相关的数据。
[0171]
所述处理器可以是专用集成电路(application specific integrated circuit，简称asic)、数字信号处理器(digital signal processor，简称dsp)、数字信号处理设备(digital signal processing device，简称dspd)、可编程逻辑器件(programmable logic device，简称pld)、现场可编程门阵列(field programmable gate array，简称fpga)、控制器、微控制器、微处理器或其他电子元件实现，用于执行一种基于多属性融合的储层特征分析方法。
[0172]
所述存储器可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(static random access memory，简称sram)，电可擦除可编程只读存储器(electrically erasable programmable read-only memory，简称eeprom)，可擦除可编程只读存储器(erasable programmable read-only memory，简称eprom)，可编程只读存储器(programmable read-only memory，简称prom)，只读存储器(read-only memory，简称rom)，磁存储器，快闪存储器，磁盘或光盘。
[0173]
多媒体组件可以包括屏幕和音频组件，所述屏幕可以是触摸屏，音频组件用于输出和/或输入音频信号。例如，音频组件可以包括一个麦克风，麦克风用于接收外部音频信号。所接收的音频信号可以被进一步存储在存储器或通过通信组件发送。音频组件还包括至少一个扬声器，用于输出音频信号。
[0174]
i/o接口为处理器和其他接口模块之间提供接口，上述其他接口模块可以是键盘，鼠标，按钮等。这些按钮可以是虚拟按钮或者实体按钮。
[0175]
通信组件用于该电子设备与其他设备之间进行有线或无线通信。无线通信，例如wi-fi，蓝牙，近场通信(near field communication，简称nfc)，2g、3g或4g，或它们中的一种或几种的组合，因此相应的该通信组件405可以包括：wi-fi模块，蓝牙模块，nfc模块。
[0176]
进一步地，上述基于多属性融合的储层特征分析方法包括以下步骤：
[0177]
步骤1：获取目标储层的若干个三维地震属性数据体，以及每个三维地震属性数据体所对应的沿层时间切片；
[0178]
步骤1.1：根据地质约束确定目标储层的n个地震属性；
[0179]
其中，所述地质约束包括目标储层的测井数据和地质认识分析，所述地质认识分
析包括目标储层的地质成因分析；
[0180]
步骤1.2：利用叠前地震道集或叠后地震属性获取与n个地震属性对应的n个三维地震属性数据体p1,p2,
…
,pn；
[0181]
步骤1.3：对所述n个三维地震属性数据体p1,p2,
…
,pn分别进行反演，以得到与所述n个三维地震属性数据体p1,p2,
…
,pn对应的n个沿层时间切片；
[0182]
具体地，结合地质约束沿层向下开时窗，以开展层位属性面的提取，获得对应的n个沿层时间切片；
[0183]
步骤2：对于目标储层的每个层位，利用沿层时间切片方程将所有沿层时间切片融合为一个综合沿层切片，以此获得目标储层多个层位的多个综合沿层切片；
[0184]
步骤2.1：对于目标储层的每个层位，利用沿层时间切片方程将所述n个沿层时间切片融合为一个综合沿层切片m，所述沿层时间切片方程为：
[0185]
m＝ω1s1 ω2s2
…
ωnsn，
[0186]
其中，n表示地震属性的个数，ω1,ω2,
…
,ωn表示n个地震属性在沿层时间切片中的权重系数，s1,s2,
…
,sn表示沿层时间切片数据，m表示综合沿层切片；
[0187]
在上述过程中，所述n个地震属性在沿层时间切片中的权重系数ω1,ω2,
…
,ωn根据不同地震属性对于目标储层发育的贡献值进行调整；
[0188]
具体地，以所述n个沿层时间切片为输入，输入到所述沿层时间切片方程中，并结合地质约束和这n个地震属性对于储层发育的贡献值不断调整输入的这n个沿层时间切片的权重系数，调整到最能反映储层特征的n个沿层时间切片的权重系数后，输出得到既能反映地层由于区域应力或者断层破裂引起的大、中、小尺度的断裂-裂缝，又能反映其中由于岩性等引起的裂缝的综合沿层切片m；
[0189]
步骤2.2：根据步骤2.1利用下式得到目标储层的i个层位的i个综合沿层切片m1,m2,
…
,mi：
[0190][0191]
其中，i表示目标储层发育区的层位数，ω
i,n
表示第n个地震属性在第i层地层的第n个沿层时间切片中的权重系数，mi表示第i层地层的综合沿层切片；
[0192]
步骤3：利用最小二乘法对所述多个综合沿层切片进行反演，得到若干个最优权重系数，利用所述若干个最优权重系数对所述若干个三维地震属性数据体进行加权处理得到三维地震属性融合体；
[0193]
步骤3.1：利用最小二乘法对i个综合沿层切片m1,m2,
…
,mi进行反演得到n个最优权重系数：
[0194]
[0195]
其中，n表示地震属性的个数，表示n个地震属性在三维地震属性融合体中的最优权重系数，i表示目标储层发育区的层位数，ω
i,n
表示第n个地震属性在第i层地层的第n个沿层时间切片中的权重系数；
[0196]
使用for循环求解上式，得到最优权重系数
[0197]
步骤3.2：利用所述步骤3.1得到的n个最优权重系数和所述步骤1.2得到的n个三维地震属性数据体p1,p2,
…
,pn进行加权处理得到三维地震属性融合体z：
[0198][0199]
步骤4：利用三维地震属性融合体z对目标储层的整体特征进行分析。
[0200]
上述方法步骤的具体实施例过程可参见实施例三，本实施例在此不再重复赘述。
[0201]
综上，本发明提供了一种基于多属性融合的储层特征分析方法、存储介质及电子设备，通过对于目标储层的每个层位，利用沿层时间切片方程将获取的所有沿层时间切片融合为多个综合沿层切片，再利用最小二乘法对多个综合沿层切片进行反演得到最优权重系数，再结合获取的若干个三维地震属性数据体进行加权处理得到三维地震属性融合体，最后再利用三维地震属性融合体对目标储层的整体特征进行分析。
[0202]
本发明通过对目标储层多个沿层切片的地质趋势和目标储层的地层趋势进行分析和控制，开展基于最小二乘算法的三维地震属性融合体的参数求取，针对地震预测中的多解性提供了确定解和唯一解，有效提高地震储层特征分析的精度和准确度，从而提高地震勘探的效率；本发明得到的三维地震属性融合体可以兼顾地层沿层切片中反映的地质规律和预测目标，又能在地震数据的驱动下形成三维空间融合体数据，迅速找到最优解，为后续的油藏地质建模、油藏精细开发等研究提供快捷可靠的三维地震成果数据体。
[0203]
需要说明的是，由于说明书附图不得着色和涂改，所以本发明附图中部分区别明显的地方比较难以显示，若有必要，可提供彩色图片。
[0204]
在本发明实施例所提供的几个实施例中，应该理解到，所揭露的系统和方法，也可以通过其它的方式实现。以上所描述的系统和方法实施例仅仅是示例性的。
[0205]
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0206]
虽然本发明所揭露的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。