一种基于波形特征箱体的储层构型边界预测方法及系统

技术领域

本发明涉及油气田开发技术领域，尤其涉及一种基于波形特征箱体的储层构型边界预测方法及系统。

背景技术

复合砂体是指多期单砂体叠置发育组合形成的复合体，内部单砂体之间在纵向、横向上彼此接触或被泥岩隔、夹层分隔。复合砂体整体被认为是地下储层，但伴随着油气田开发程度加深，储层精细描述面临着在识别储层的基础上需要进一步厘清复合砂体内部单砂体的叠置关系，以此研究储层连通性，指导井位部署与油藏配产。因此，预测单砂体之间的接触位置与接触关系(即构型边界)的意义非常重大，但现有技术主要利用地震属性预测复合砂体，针对更精细的单砂体构型边界无成熟预测方法，多利用多属性融合、聚类等来变通实现，但该做法缺乏理论依据，效果整体一般。

发明内容

本发明提供一种基于波形特征箱体的储层构型边界预测方法及系统，用以解决现有技术无法准确地预测储层构型边界的缺陷。

本发明提供一种基于波形特征箱体的储层构型边界预测方法，包括：

修正储层顶解释层面和储层底解释层面，使储层顶解释层面位于地震数据反射波谷位置，储层底解释层面位于地震数据反射波峰位置；

根据储层顶解释层面和储层底解释层面确定时窗范围，得到每个地震道的储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值；

根据储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值，得到构型系数；

根据储层段最大能量值、储层段地震波形偏度值以及构型系数，绘制每个地震道的波形特征箱体，以预测储层构型边界。

根据本发明提供的一种基于波形特征箱体的储层构型边界预测方法，所述修正储层顶解释层面和储层底解释层面，使储层顶解释层面位于地震数据反射波谷位置，储层底解释层面位于地震数据反射波峰位置，包括：

根据油田储层地质数据，构建复合砂体纵波阻抗模型；

根据复合砂体纵波阻抗模型，得到振幅反射能量、地震子波以及地震反射系数；

根据振幅反射能量、地震子波以及地震反射系数进行褶积运算，得到地震反射剖面；

基于地震反射剖面，修正储层顶解释层面和储层底解释层面，使储层顶解释层面位于地震数据反射波谷位置，储层底解释层面位于地震数据反射波峰位置。

根据本发明提供的一种基于波形特征箱体的储层构型边界预测方法，所述根据储层顶解释层面和储层底解释层面确定时窗范围，得到每个地震道的储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值，包括：

根据储层段地震振幅值，通过第一表达式得到每个地震道的储层段最大能量值，其中第一表达式为：E

根据储层段地震振幅值，通过第二表达式得到每个地震道的储层段地震振幅值的中值，其中第二表达式为：m

根据储层段地震振幅值，通过第三表达式得到每个地震道的储层段地震振幅值的均值，其中第三表达式为：a

根据储层段地震振幅值，通过第四表达式得到每个地震道的储层段地震波形偏度值，其中第四表达式为：

γ

根据本发明提供的一种基于波形特征箱体的储层构型边界预测方法，所述根据储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值，得到构型系数，具体为：

根据储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值，通过第五表达式得到构型系数，其中第五表达式为：

c

根据本发明提供的一种基于波形特征箱体的储层构型边界预测方法，得到对构型系数具有影响力的计算阈值的方法包括：

获取样本数据；

根据样本数据，得到复合砂体纵波阻抗模型的波形特征箱体以及实际砂体构型边界；

通过调节对构型系数具有影响力的计算阈值，使复合砂体纵波阻抗模型的波形特征箱体与实际砂体构型边界的吻合度达到第一预设范围，得到最终计算阈值。

根据本发明提供的一种基于波形特征箱体的储层构型边界预测方法，通过第六表达式得到复合砂体纵波阻抗模型的波形特征箱体与实际砂体构型边界的吻合度，其中第六表达式为：

acc＝n/N·100％

acc表示复合砂体纵波阻抗模型的波形特征箱体与实际砂体构型边界的吻合度，n表示复合砂体纵波阻抗模型的波形特征箱体反映的砂体构型边界与实际砂体构型边界一致的地震道数，N表示复合砂体纵波阻抗模型的总地震道数。

根据本发明提供的一种基于波形特征箱体的储层构型边界预测方法，所述根据储层段最大能量值、储层段地震波形偏度值以及构型系数，绘制每个地震道的波形特征箱体，以预测储层构型边界，包括：

根据储层段最大能量值，得到波形特征箱体的高度；

根据储层段地震波形偏度值，得到波形特征箱体的宽度；

根据构型系数，得到波形特征箱体的颜色；

根据波形特征箱体的高度、宽度、颜色，绘制每个地震道的波形特征箱体，并将每个地震道的波形特征箱体按照预设间隔排列于地震反射剖面；

观察波形特征箱体之间的变化，将变化程度达到第二预设范围的位置作为预测的储层构型边界，预测的储层构型边界两侧的波形特征箱体之间的变化反映单砂体之间的接触关系的变化。

本发明还提供一种基于波形特征箱体的储层构型边界预测系统，包括：

修正模块，用于：修正储层顶解释层面和储层底解释层面，使储层顶解释层面位于地震数据反射波谷位置，储层底解释层面位于地震数据反射波峰位置；

数值计算模块，用于：根据储层顶解释层面和储层底解释层面确定时窗范围，得到每个地震道的储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值；

构型系数得到模块，用于：根据储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值，得到构型系数；

储层构型边界预测模块，用于：根据储层段最大能量值、储层段地震波形偏度值以及构型系数，绘制每个地震道的波形特征箱体，以预测储层构型边界。

本发明还提供一种电子设备，包括处理器和存储有计算机程序的存储器，其特征在于，所述处理器执行所述计算机程序时实现上述任一种所述的基于波形特征箱体的储层构型边界预测方法。

本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述任一种所述的基于波形特征箱体的储层构型边界预测方法。

本发明还提供一种计算机程序产品，所述计算机程序产品包括计算机程序，计算机程序可存储在非暂态计算机可读存储介质上，所述计算机程序被处理器执行时，计算机能够执行上述任一种所述的基于波形特征箱体的储层构型边界预测方法。

地震反射波形包含丰富的振幅、频率、相位等信息，地震波形的变化能够反映纵向上单砂体构型样式的变化，本发明提供的一种基于波形特征箱体的储层构型边界预测方法及系统，将储层顶解释层面修正于地震数据反射波谷位置，储层底解释层面修正于地震数据反射波峰位置，根据储层段最大能量值、储层段地震波形偏度值以及构型系数绘制每个地震道的波形特征箱体，利用地震波形特征变化信息直接、准确地预测储层构型边界，提供储层精细描述，优化油气田开发的速度和质量，为油气田开发带来重要价值。

附图说明

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

图1为本发明提供的一种基于波形特征箱体的储层构型边界预测方法的流程示意图之一。

图2(a)示出复合砂体纵波阻抗模型，图2(b)示出地震数据、储层顶解释层面储层底解释层面、以及波形特征箱体图。

图3示出渤海A油田沿S1砂体井轨迹地震反射剖面及通过本发明提供的一种基于波形特征箱体的储层构型边界预测方法得到的波形特征箱体图。

图4示出渤海A油田沿S1砂体井轨迹地震反射剖面及测井随钻解释。

图5为本发明提供的一种基于波形特征箱体的储层构型边界预测系统的结构示意图。

图6为本发明提供的电子设备的结构示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例，它们不应该理解成对本发明的限制。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。在本发明的描述中，需要理解的是，所用到的术语仅仅是用于描述的目的，而不能理解为指示或暗示相对重要性。

下面结合图1-图6描述本发明提供的基于波形特征箱体的储层构型边界预测方法、系统、设备及介质。

图1是本发明提供的基于波形特征箱体的储层构型边界预测方法的流程示意图。

参照图1，本发明提供的一种基于波形特征箱体的储层构型边界预测方法，可以包括：

步骤S110、修正储层顶解释层面和储层底解释层面，使储层顶解释层面位于地震数据反射波谷位置，储层底解释层面位于地震数据反射波峰位置；

步骤S120、根据储层顶解释层面和储层底解释层面确定时窗范围，得到每个地震道的储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值；

步骤S130、根据储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值，得到构型系数；

步骤S140、根据储层段最大能量值、储层段地震波形偏度值以及构型系数，绘制每个地震道的波形特征箱体，以预测储层构型边界。

需要说明的是，本发明提供的基于波形特征箱体的储层构型边界预测方法的执行主体可以是任何符合技术要求的终端侧设备，例如基于波形特征箱体的储层构型边界预测装置等。

地震反射波形包含丰富的振幅、频率、相位等信息，地震波形的变化能够反映纵向上单砂体构型样式的变化，本发明提供的基于波形特征箱体的储层构型边界预测方法及系统基于此利用地震波形特征变化信息预测储层构型边界，能够得到更加准确及精细的储层构型边界描述，有助于油气田开发。

在一种实施例中，步骤S110可以包括：

根据油田储层地质数据，构建复合砂体纵波阻抗模型；

根据复合砂体纵波阻抗模型，得到振幅反射能量、地震子波以及地震反射系数；

根据振幅反射能量、地震子波以及地震反射系数进行褶积运算，得到地震反射剖面；

基于地震反射剖面，修正储层顶解释层面和储层底解释层面，使储层顶解释层面位于地震数据反射波谷位置，储层底解释层面位于地震数据反射波峰位置。

具体的，以渤海A油田为例，基于渤海A油田储层地质条件，确定该油田砂岩速度约2450m/s，泥岩速度约2650m/s，砂岩密度约2.10g/cm3，泥岩密度约2.25g/cm3，地震子波反射能量约为5000，复合砂体厚度约15m，砂体宽厚比约50:1，根据上述油田储层地质数据，构建复合砂体纵波阻抗模型(参见图2)，得到振幅反射能量、地震子波以及地震反射系数，再根据褶积表达式得到地震反射剖面。

褶积表达式为：

d＝s·w*r (1)

褶积表达式中，d表示地震反射剖面，s表示振幅反射能量，取s＝5000，w表示地震子波，采用30Hz主频的Ricker子波(该子波是勘探地球物理领域最为常用的理论子波)，r表示地震反射系数，可根据上述构建的复合砂体纵波阻抗模型结合式(2)得到。

r

式(2)中，i表示第i采样点，r

X(n)＝Σh(m)×x(n-m) (3)

式(3)中，n代表序列，m为偏移量，X(n)表示h(m)与x(n-m)的褶积。

以地震反射剖面为基础，解释储层顶解释层面和储层底解释层面，可以通过采用现有软件中的波峰、波谷层位自动追踪解释功能修正储层顶解释层面和储层底解释层面，确保储层顶解释层面位于地震数据反射波谷位置(即反射值为负，同时振幅最大的位置)，储层底解释层面位于地震数据反射波峰(即反射值为正，同时振幅最大的位置)。

在一种实施例中，步骤S120可以包括：

根据储层段地震振幅值，通过第一表达式得到每个地震道的储层段最大能量值，其中第一表达式为：E

根据储层段地震振幅值，通过第二表达式得到每个地震道的储层段地震振幅值的中值，其中第二表达式为：m

根据储层段地震振幅值，通过第三表达式得到每个地震道的储层段地震振幅值的均值，其中第三表达式为：a

根据储层段地震振幅值，通过第四表达式得到每个地震道的储层段地震波形偏度值，其中第四表达式为：

式(4)中，γ

在一种实施例中，步骤S130可以具体为：

根据储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值，通过第五表达式得到构型系数，其中第五表达式为：

式(5)中，c

具体的，可以通过以下步骤得到对构型系数具有影响力的计算阈值：

获取样本数据；

根据样本数据，得到复合砂体纵波阻抗模型的波形特征箱体以及实际砂体构型边界；

通过调节对构型系数具有影响力的计算阈值，使复合砂体纵波阻抗模型的波形特征箱体与实际砂体构型边界的吻合度达到第一预设范围，得到最终计算阈值。

同样以渤海A油田为例，以渤海A油田的油田储层地质数据作为样本数据，可以得到其复合砂体纵波阻抗模型，继而得到复合砂体纵波阻抗模型的波形特征箱体，然后可以根据渤海A油田的实钻井得到实际砂体构型边界，在调整对构型系数具有影响力的计算阈值时，当复合砂体纵波阻抗模型的波形特征箱体与实际砂体构型边界的吻合度达到第一预设范围时，即可得到最终计算阈值。

需要说明的是，本实施例中，第一预设范围为80％以上，即，通过调节对构型系数具有影响力的计算阈值，当复合砂体纵波阻抗模型的波形特征箱体与实际砂体构型边界的吻合度达到80％时，以此时的计算阈值作为最终计算阈值。

可以通过第六表达式得到复合砂体纵波阻抗模型的波形特征箱体与实际砂体构型边界的吻合度，其中第六表达式为：

acc＝n/N·100％ (6)

式(6)中，acc表示复合砂体纵波阻抗模型的波形特征箱体与实际砂体构型边界的吻合度，n表示复合砂体纵波阻抗模型的波形特征箱体反映的砂体构型边界与实际砂体构型边界一致的地震道数，N表示复合砂体纵波阻抗模型的总地震道数。

在一种实施例中，步骤S140可以包括：

根据储层段最大能量值，得到波形特征箱体的高度；

根据储层段地震波形偏度值，得到波形特征箱体的宽度；

根据构型系数，得到波形特征箱体的颜色；

根据波形特征箱体的高度、宽度、颜色，绘制每个地震道的波形特征箱体，并将每个地震道的波形特征箱体按照预设间隔排列于地震反射剖面；

观察波形特征箱体之间的变化，将变化程度达到第二预设范围的位置作为预测的储层构型边界，预测的储层构型边界两侧的波形特征箱体之间的变化反映单砂体之间的接触关系的变化。

需要说明的是，每个波形特征箱体包括高度、宽度、颜色三个维度的数据，在得到储层段最大能量值、储层段地震波形偏度值以及构型系数后，需要先对储层段最大能量值和储层段地震波形偏度值进行横向平滑处理，再以平滑后的储层段最大能量值、储层段地震波形偏度值以及构型系数作为绘制波形特征箱体的数据基础。

根据储层段最大能量值绘制波形特征箱体的高度，具体方法：地震反射剖面中每个地震道均具有一个储层段最大能量值，将整个地震反射剖面的所有储层段最大能量值归一化为0至1数值范围内的数值，记归一化后第k地震道的储层段最大能量值为F

根据储层段地震波形偏度值绘制波形特征箱体的宽度，具体方法：地震反射剖面中每个地震道均具有一个储层段地震波形偏度值，将整个地震反射剖面的所有储层段地震波形偏度值归一化为0至1数值范围内的数值，记归一化后第k地震道的波形偏度为Q

根据构型系数绘制波形特征箱体的颜色，具体做法为：地震反射剖面中每个地震道均具有一个构型系数，在绘制第k波形特征箱体时，构型系数c

按照上述方法绘制好每个地震道的波形特征箱体后，将波形特征箱体按照预设间隔(本实施例中，取每隔1个地震道绘制一个波形特征箱体，可根据实际地震反射剖面道间距适当调整)绘制在地震反射剖面上，即完成在地震反射剖面上绘制每道的波形特征箱体。

然后观察波形特征箱体之间的变化，将变化程度达到第二预设范围的位置作为预测的储层构型边界(本实施例将相邻波形特征箱体的颜色发生直接变化(例如黑变绿、红变绿、绿变黑等等)或相邻波形特征箱体的胖瘦具有较为明显变化(例如左边波形特征箱体的宽度只有右边波形特征箱体的宽度的十分之一等等)的位置作为预测的储层构型边界)，预测的储层构型边界两侧的波形特征箱体之间的变化反映单砂体之间的接触关系的变化，接触关系例如有两期砂体前下后上叠置、两期砂体前上后下叠置、单期砂体无叠置等等，可以根据实际情况判断波形特征箱体之间的变化对构型边界的预测是否具有参考性。

本发明提取渤海A油田沿着S1复合砂体开发井轨迹的地震反射剖面，将储层顶解释层面和储层底解释层面分别修正于地震数据反射波峰位置和地震数据反射波谷位置，根据储层段最大能量值、储层段地震波形偏度值以及构型系数，在沿井轨迹地震反射剖面中绘制每个地震道的波形特征箱体(参见图3)，变化程度达到第二预设范围的位置(例如相邻波形特征箱体的颜色发生直接变化的位置)作为预测的储层构型边界，图4示出钻遇该砂体的开发井随钻岩性解释，预测结果与随钻岩性解释结果一致，证实了该井钻遇该砂体的构型边界，验证了通过本发明提供的基于波形特征箱体的储层构型边界预测方法得到的预测结果的合理性和准确性。

地震反射波形包含丰富的振幅、频率、相位等信息，地震波形的变化能够反映纵向上单砂体构型样式的变化，本发明提供的一种基于波形特征箱体的储层构型边界预测方法及系统，将储层顶解释层面修正于地震数据反射波谷位置，储层底解释层面修正于地震数据反射波峰位置，根据储层段最大能量值、储层段地震波形偏度值以及构型系数绘制每个地震道的波形特征箱体，利用地震波形特征变化信息直接、准确地预测储层构型边界，提供储层精细描述，优化油气田开发的速度和质量，为油气田开发带来重要价值。

下面对本发明提供的基于波形特征箱体的储层构型边界预测系统进行描述，下文描述的基于波形特征箱体的储层构型边界预测系统与上文描述的基于波形特征箱体的储层构型边界预测方法可相互对应参照。

参照图5，本发明提供的一种基于波形特征箱体的储层构型边界预测系统，可以包括：

修正模块510，用于：修正储层顶解释层面和储层底解释层面，使储层顶解释层面位于地震数据反射波谷位置，储层底解释层面位于地震数据反射波峰位置；

数值计算模块520，用于：根据储层顶解释层面和储层底解释层面确定时窗范围，得到每个地震道的储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值；

构型系数得到模块530，用于：根据储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值，得到构型系数；

储层构型边界预测模块540，用于：根据储层段最大能量值、储层段地震波形偏度值以及构型系数，绘制每个地震道的波形特征箱体，以预测储层构型边界。

在一种实施例中，修正模块可以包括：

模型构建子模块，用于：根据油田储层地质数据，构建复合砂体纵波阻抗模型；

数据得到子模块，用于：根据复合砂体纵波阻抗模型，得到振幅反射能量、地震子波以及地震反射系数；

褶积子模块，用于：根据振幅反射能量、地震子波以及地震反射系数进行褶积运算，得到地震反射剖面；

修正子模块，用于：基于地震反射剖面，修正储层顶解释层面和储层底解释层面，使储层顶解释层面位于地震数据反射波谷位置，储层底解释层面位于地震数据反射波峰位置。

在一种实施例中，数值计算模块可以包括：

第一计算子模块，用于：根据储层段地震振幅值，通过第一表达式得到每个地震道的储层段最大能量值，其中第一表达式为：E

第二计算子模块，用于：根据储层段地震振幅值，通过第二表达式得到每个地震道的储层段地震振幅值的中值，其中第二表达式为：m

第三计算子模块，用于：根据储层段地震振幅值，通过第三表达式得到每个地震道的储层段地震振幅值的均值，其中第三表达式为：a

第四计算子模块，用于：根据储层段地震振幅值，通过第四表达式得到每个地震道的储层段地震波形偏度值，其中第四表达式为：

γ

在一种实施例中，构型系数得到模块具体用于：

根据储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值，通过第五表达式得到构型系数，其中第五表达式为：

c

需要说明的是，得到对构型系数具有影响力的计算阈值的方法包括：

获取样本数据；

根据样本数据，得到复合砂体纵波阻抗模型的波形特征箱体以及实际砂体构型边界；

通过调节对构型系数具有影响力的计算阈值，使复合砂体纵波阻抗模型的波形特征箱体与实际砂体构型边界的吻合度达到第一预设范围，得到最终计算阈值。

需要说明的是，通过第六表达式得到复合砂体纵波阻抗模型的波形特征箱体与实际砂体构型边界的吻合度，其中第六表达式为：

acc＝n/N·100％

acc表示复合砂体纵波阻抗模型的波形特征箱体与实际砂体构型边界的吻合度，n表示复合砂体纵波阻抗模型的波形特征箱体反映的砂体构型边界与实际砂体构型边界一致的地震道数，N表示复合砂体纵波阻抗模型的总地震道数。

在一种实施例中，储层构型边界预测模块可以包括：

第一绘制子模块，用于：根据储层段最大能量值，得到波形特征箱体的高度；

第二绘制子模块，用于：根据储层段地震波形偏度值，得到波形特征箱体的宽度；

第三绘制子模块，用于：根据构型系数，得到波形特征箱体的颜色；

第四绘制子模块，用于：根据波形特征箱体的高度、宽度、颜色，绘制每个地震道的波形特征箱体，并将每个地震道的波形特征箱体按照预设间隔排列于地震反射剖面；

预测子模块，用于：观察波形特征箱体之间的变化，将变化程度达到第二预设范围的位置作为预测的储层构型边界，预测的储层构型边界两侧的波形特征箱体之间的变化反映单砂体之间的接触关系的变化。

图6示例了一种电子设备的实体结构示意图，如图6所示，该电子设备可以包括：处理器(processor)810、通信接口(Communications Interface)820、存储器(memory)830和通信总线840，其中，处理器810，通信接口820，存储器830通过通信总线840完成相互间的通信。处理器810可以调用存储器830中的逻辑指令，以执行基于波形特征箱体的储层构型边界预测方法，该方法包括：

修正储层顶解释层面和储层底解释层面，使储层顶解释层面位于地震数据反射波谷位置，储层底解释层面位于地震数据反射波峰位置；

根据储层顶解释层面和储层底解释层面确定时窗范围，得到每个地震道的储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值；

根据储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值，得到构型系数；

根据储层段最大能量值、储层段地震波形偏度值以及构型系数，绘制每个地震道的波形特征箱体，以预测储层构型边界。

此外，上述的存储器830中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

另一方面，本发明还提供一种计算机程序产品，所述计算机程序产品包括计算机程序，计算机程序可存储在非暂态计算机可读存储介质上，所述计算机程序被处理器执行时，计算机能够执行上述各方法所提供的基于波形特征箱体的储层构型边界预测方法，该方法包括：

修正储层顶解释层面和储层底解释层面，使储层顶解释层面位于地震数据反射波谷位置，储层底解释层面位于地震数据反射波峰位置；

根据储层顶解释层面和储层底解释层面确定时窗范围，得到每个地震道的储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值；

根据储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值，得到构型系数；

根据储层段最大能量值、储层段地震波形偏度值以及构型系数，绘制每个地震道的波形特征箱体，以预测储层构型边界。

又一方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各方法提供的基于波形特征箱体的储层构型边界预测方法，该方法包括：

修正储层顶解释层面和储层底解释层面，使储层顶解释层面位于地震数据反射波谷位置，储层底解释层面位于地震数据反射波峰位置；

根据储层顶解释层面和储层底解释层面确定时窗范围，得到每个地震道的储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值；

根据储层段最大能量值、储层段地震振幅值以及储层段地震波形偏度值，得到构型系数；

根据储层段最大能量值、储层段地震波形偏度值以及构型系数，绘制每个地震道的波形特征箱体，以预测储层构型边界。

以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。

通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。

最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。