水平井产油量的预测方法和装置与流程

本公开涉及石油开采领域，尤其涉及一种水平井产量的预测方法和装置。

背景技术：

水平井是指在石油开采中井身沿着水平方向在油层中钻进一定长度的井。对水平井产油量的预测是水平井油藏工程中的重要研究内容，也是制定水平井合理工作制度、预测水平井的生产动态以及确定水平井的开发方式的重要依据。因此，对水平井产油量的预测具有重要的理论和实际意义。

通过水平井开采石油时，通常是在水平井的驱动范围中开采石油，该驱动范围是指水平井可以开采到石油的范围。目前认为水平井的驱动范围中只有矩形部分内产油，因此通常是通过该矩形部分对水平井产油量进行预测。

然而，实际中水平井驱动范围中并不是只有矩形部内产油，因此，通过上述方法对水平井产油量进行预测的预测结果不精确。

技术实现要素：

本公开提供一种水平井产油量的预测方法和装置。可以解决对水平井产油量进行预测的预测结果不精确的问题。

第一方面，提供一种水平井产油量的预测方法，包括：

确定待预测的水平井的驱动范围中矩形部分的水平长度、所述驱动范围边界与所述水平井之间的压力差、泄油半径、所述水平井的井筒半径以及所述水平井与所述驱动范围下边界之间的距离，所述驱动范围包括所述矩形部分和两端部分；

确定所述水平井所在的油层的渗透率、所述油层的粘度、所述油层的体积系数以及所述油层的厚度；

根据所述水平长度、所述压力差、所述泄油半径、所述井筒半径、所述水平井与所述驱动范围下边界之间的距离、所述油层的渗透率、所述油层的粘度、所述油层的体积系数以及所述油层的厚度，确定所述矩形部分对应的第一产油量；

根据所述压力差、所述泄油半径、所述井筒半径、所述油层的渗透率、所述油层的粘度、所述油层的体积系数以及所述油层的厚度，确定所述两端部分对应的第二产油量；

确定所述水平井的渗透率变异系数和流度比，并根据所述渗透率变异系数和所述流度比确定所述油层的波及系数，所述波及系数用于描述所述油层的非均质性对所述水平井产油量的影响程度，所述波及系数与所述影响程度呈正比；

根据所述第一产油量、所述第二产油量和所述波及系数，确定所述水平井对应的第三产油量。

在一个可能的实现方式中，所述根据所述水平长度、所述压力差、所述泄油半径、所述井筒半径、所述水平井与所述驱动范围下边界之间的距离、所述油层的渗透率、所述油层的粘度、所述油层的体积系数以及所述油层的厚度，确定所述矩形部分对应的第一产油量，包括：

根据所述水平长度、所述压力差、所述泄油半径、所述井筒半径、所述水平井与所述驱动范围下边界之间的距离、所述油层的渗透率、所述油层的粘度、所述油层的体积系数以及所述油层的厚度，通过如下公式一确定所述第一产油量：

公式一：

其中，所述q1为所述第一产油量，所述l为所述水平长度，所述δp为所述压力差，所述re为所述泄油半径，所述rw为所述井筒半径，所述zw为所述水平井与所述驱动范围下边界之间的距离，所述k为所述油层的渗透率，所述μ0为所述油层的粘度，所述b为所述油层的体积系数，所述h为所述油层的厚度，所述ch为双曲余弦函数，所述cos为双曲余弦函数，所述ln为对数函数。

在一个可能的实现方式中，所述根据所述压力差、所述泄油半径、所述井筒半径、所述油层的渗透率、所述油层的粘度、所述油层的体积系数以及所述油层的厚度，确定所述两端部分对应的第二产油量，包括：

根据所述压力差、所述泄油半径、所述井筒半径、所述油层的渗透率、所述油层的粘度、所述油层的体积系数以及所述油层的厚度，通过如下公式二确定所述第二产油量：

公式二：

其中，所述q2为所述第二产油量，所述δp为所述压力差，所述re为所述泄油半径，所述rw为所述井筒半径，所述k为所述油层的渗透率，所述μ0为所述油层的粘度，所述b为所述油层的体积系数，所述h为所述油层的厚度，所述ln为对数函数。

在一个可能的实现方式中，所述根据所述渗透率变异系数和所述流度比确定所述油层的波及系数，包括：

根据所述渗透率变异系数和所述流度比，通过如下公式三确定所述波及系数：

公式三：

其中，所述c为所述波及系数，所述γ为所述渗透率变异系数，所述m为所述流度比。

在一个可能的实现方式中，所述根据所述第一产油量、所述第二产油量和所述波及系数，确定所述水平井对应的第三产油量，包括：

确定所述第一产油量和所述第二产油量之和，得到第四产油量；

确定所述第四产油量与所述波及系数的乘积，得到所述第三产油量。

第二方面，提供一种水平井产油量的预测装置，所述装置包括：

第一确定模块，用于确定待预测的水平井的驱动范围中矩形部分的水平长度、所述驱动范围边界与所述水平井之间的压力差、泄油半径、所述水平井的井筒半径以及所述水平井与所述驱动范围下边界之间的距离，所述驱动范围包括所述矩形部分和两端部分；

第二确定模块，用于确定所述水平井所在的油层的渗透率、所述油层的粘度、所述油层的体积系数以及所述油层的厚度；

第三确定模块，用于根据所述水平长度、所述压力差、所述泄油半径、所述井筒半径、所述水平井与所述驱动范围下边界之间的距离、所述油层的渗透率、所述油层的粘度、所述油层的体积系数以及所述油层的厚度，确定所述矩形部分对应的第一产油量；

第四确定模块，用于根据所述压力差、所述泄油半径、所述井筒半径、所述油层的渗透率、所述油层的粘度、所述油层的体积系数以及所述油层的厚度，确定所述两端部分对应的第二产油量；

第五确定模块，用于确定所述水平井的渗透率变异系数和流度比，并根据所述渗透率变异系数和所述流度比确定所述油层的波及系数，所述波及系数用于描述所述油层的非均质性对所述水平井产油量的影响程度，所述波及系数与所述影响程度呈正比；

第六确定模块，用于根据所述第一产油量、所述第二产油量和所述波及系数，确定所述水平井对应的第三产油量。

在一个可能的实现方式中，所述第三确定模块包括：

第一确定子模块，用于根据所述水平长度、所述压力差、所述泄油半径、所述井筒半径、所述水平井与所述驱动范围下边界之间的距离、所述油层的渗透率、所述油层的粘度、所述油层的体积系数以及所述油层的厚度，通过如下公式一确定所述第一产油量：

公式一：

其中，所述q1为所述第一产油量，所述l为所述水平长度，所述δp为所述压力差，所述re为所述泄油半径，所述rw为所述井筒半径，所述zw为所述水平井与所述驱动范围下边界之间的距离，所述k为所述油层的渗透率，所述μ0为所述油层的粘度，所述b为所述油层的体积系数，所述h为所述油层的厚度，所述ch为双曲余弦函数，所述cos为双曲余弦函数，所述ln为对数函数。

在一个可能的实现方式中，所述第四确定模块包括：

第二确定子模块，用于根据所述压力差、所述泄油半径、所述井筒半径、所述油层的渗透率、所述油层的粘度、所述油层的体积系数以及所述油层的厚度，通过如下公式二确定所述第二产油量：

公式二：

其中，所述q2为所述第二产油量，所述δp为所述压力差，所述re为所述泄油半径，所述rw为所述井筒半径，所述k为所述油层的渗透率，所述μ0为所述油层的粘度，所述b为所述油层的体积系数，所述h为所述油层的厚度，所述ln为对数函数。

在一个可能的实现方式中，所述第五确定模块包括：

第三确定子模块，用于根据所述渗透率变异系数和所述流度比，通过如下公式三确定所述波及系数：

公式三：

其中，所述c为所述波及系数，所述γ为所述渗透率变异系数，所述m为所述流度比。

在一个可能的实现方式中，所述第六确定模块包括：

第四确定子模块，用于确定所述第一产油量和所述第二产油量之和，得到第四产油量；

第五确定子模块，用于确定所述第四产油量与所述波及系数的乘积，得到所述第三产油量。

第三方面，提供一种水平井产油量的预测装置，所述装置包括：

处理器；

用于存储处理器可执行指令的存储器；

其中，所述处理器被配置为执行上述第一方面所述的任一项方法的步骤。

第四方面，提供一种计算机可读存储介质，所述计算机可读存储介质上存储有指令，所述指令被处理器执行时实现上述第一方面所述的任一项方法的步骤。

第五方面，提供了一种包含指令的计算机程序产品，当其在计算机上运行时，使得计算机执行上述第一方面任一项所述方法的步骤。

本公开的实施例提供的技术方案至少可以包括以下有益效果：

在本公开实施例中，确定待预测的水平井的驱动范围中矩形部分的水平长度、驱动范围边界与水平井之间的压力差、泄油半径、水平井的井筒半径、水平井与驱动范围下边界之间的距离、水平井所在的油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度之后，可以确定矩形范围中矩形部分对应的第一产油量以及两端部分对应的第二产油量；然后在确定油层的波及系数后，可以根据矩形部分的第一产油量、两端部分的第二产油量和油层的波及系数，确定水平井对应的第三产油量，提高了对水平井产油量进行预测的精确度。

应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本公开。

附图说明

此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本发明的实施例，并与说明书一起用于解释本发明的原理。

图1为水平井的驱动范围中两端部分的流场强度低于矩形部分的流场强度的示意图。

图2是根据一示例性实施例示出的一种水平井产油量的预测方法的流程图。

图3是根据一示例性实施例示出的一种水平井产油量的预测方法的流程图。

图4为根据一示例性实施例示出的一种水平井驱动范围的示意图。

图5为根据一示例性实施例示出的一种水平井驱动范围中两端部分的示意图。

图6是根据一示例性实施例示出的一种水平井产油量的预测装置的框图。

图7是根据一示例性实施例示出的一种终端的结构框图。

具体实施方式

这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反，它们仅是与本发明的一些方面相一致的装置和方法的例子。

在对本公开实施例进行详细的解释说明之前，先对本公开实施例的应用场景进行介绍：

水平井以增大油气层的裸露面积进而提高开采石油的效率的优势在石油的开采中受到了广泛地使用。水平井的井斜角达到或接近90°，井身沿着水平方向钻进一定的长度，适用于薄的油层或裂缝性油层。因此对水平井产油量的预测是目前研究者争先研究的课题。

自20世纪40年代起，国内外学者对水平井产油量的预测作了大量的研究，提出了一系列水平井产油量的预测方法，比较常用的有苏联的borisov(鲍里索夫)、法国的giger(基格)、美国joshi(乔西)以及我国的陈元千等人研究的对水平井产油量的预测的计算方法。这些方法是将水平井的驱动范围定义为椭圆形，并利用与椭圆形面积等值的拟圆形的驱动半径和产油量等的物理概念推导出来的。这些方法认为水平井在驱动范围内的流场是均匀分布的。而实际的研究结果表明，水平井的驱动范围是由一个矩形部分和两端部分的两个半圆形构成，而且两端部分的流场强度明显低于矩形部分的流场强度，如图1所示，图1为水平井的驱动范围中两端部分的流场强度低于矩形部分的流场强度的示意图。因此依据上述四种方法所预测的水平井产油量与实际水平井产油量相比明显偏高。

1996年，我国的窦宏恩利用镜像反映原理推导出了水平井产油量计算方法，这种方法从流场分布的角度讲，仅仅考虑了水平井驱动范围内矩形部分的产量，而没有考虑水平井两端部分的产量。因此计算的水平井产油量与实际水平井产油量相比明显偏低。而且窦宏恩公式仅适用于水平井位于油层中部位置的情况，适用范围窄。

图2是根据一示例性实施例示出的一种水平井产油量的预测方法的流程图，如图2所示，应用于终端，包括以下步骤。

在步骤201中，确定待预测的水平井的驱动范围中矩形部分的水平长度、该驱动范围边界与该水平井之间的压力差、泄油半径、该水平井的井筒半径以及该水平井与该驱动范围下边界之间的距离，该驱动范围包括该矩形部分和两端部分。

在步骤202中，确定该水平井所在的油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度。

在步骤203中，根据该水平长度、该压力差、该泄油半径、该井筒半径、该水平井与该驱动范围下边界之间的距离、该油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度，确定该矩形部分对应的第一产油量。

在步骤204中，根据该压力差、该泄油半径、该井筒半径、该油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度，确定该两端部分对应的第二产油量。

在步骤205中，确定该水平井的渗透率变异系数和流度比，并根据该渗透率变异系数和该流度比确定该油层的波及系数，该波及系数用于描述该油层的非均质性对该水平井产油量的影响程度，该波及系数与该影响程度呈正比。

在步骤206中，根据该第一产油量、该第二产油量和该波及系数，确定该水平井对应的第三产油量。

在一个可能的实现方式中，根据该水平长度、该压力差、该泄油半径、该井筒半径、该水平井与该驱动范围下边界之间的距离、该油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度，确定该矩形部分对应的第一产油量，包括：

根据该水平长度、该压力差、该泄油半径、该井筒半径、该水平井与该驱动范围下边界之间的距离、该油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度，通过如下公式一确定该第一产油量：

公式一：

其中，该q1为该第一产油量，该l为该水平长度，该δp为该压力差，该re为该泄油半径，该rw为该井筒半径，该zw为该水平井与该驱动范围下边界之间的距离，该k为该油层的渗透率，该μ0为该油层的粘度，该b为该油层的体积系数，该h为该油层的厚度，该ch为双曲余弦函数，该cos为双曲余弦函数，该ln为对数函数。

在一个可能的实现方式中，根据该压力差、该泄油半径、该井筒半径、该油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度，确定该两端部分对应的第二产油量，包括：

根据该压力差、该泄油半径、该井筒半径、该油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度，通过如下公式二确定该第二产油量：

公式二：

其中，该q2为该第二产油量，该δp为该压力差，该re为该泄油半径，该rw为该井筒半径，该k为该油层的渗透率，该μ0为该油层的粘度，该b为该油层的体积系数，该h为该油层的厚度，该ln为对数函数。

在一个可能的实现方式中，根据该渗透率变异系数和该流度比确定该油层的波及系数，包括：

根据该渗透率变异系数和该流度比，通过如下公式三确定该波及系数：

公式三：

其中，该c为该波及系数，该γ为该渗透率变异系数，该m为该流度比。

在一个可能的实现方式中，根据该第一产油量、该第二产油量和该波及系数，确定该水平井对应的第三产油量，包括：

确定该第一产油量和该第二产油量之和，得到第四产油量；

确定该第四产油量与该波及系数的乘积，得到该第三产油量。

在本公开实施例中，确定待预测的水平井的驱动范围中矩形部分的水平长度、驱动范围边界与水平井之间的压力差、泄油半径、水平井的井筒半径、水平井与驱动范围下边界之间的距离、水平井所在的油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度之后，可以确定矩形范围中矩形部分对应的第一产油量以及两端部分对应的第二产油量；然后在确定油层的波及系数后，可以根据矩形部分的第一产油量、两端部分的第二产油量和油层的波及系数，确定水平井对应的第三产油量，提高了对水平井产油量进行预测的精确度。

上述所有可选技术方案，均可按照任意结合形成本公开的可选实施例，本公开实施例对此不再一一赘述。

图3是根据一示例性实施例示出的一种水平井产油量的预测方法的流程图，本实施例将对图2所示的实施例进行展开说明。如图3所示，应用于终端，该方法包括以下步骤。

在步骤301中，终端确定待预测的水平井的驱动范围中矩形部分的水平长度、驱动范围边界与水平井之间的压力差、泄油半径、水平井的井筒半径以及水平井与驱动范围下边界之间的距离。

终端在确定待预测的水平井的驱动范围中矩形部分的水平长度、驱动范围边界与水平井之间的压力差、泄油半径、水平井的井筒半径以及水平井与驱动范围下边界之间的距离时，在一种可能的实现方式中，终端可以接收用户输入的矩形部分的水平长度、驱动范围边界与水平井之间的压力差、泄油半径、水平井的井筒半径以及水平井与驱动范围下边界之间的距离。在另一种可能的实现方式中，水平井的配置信息中可能包括矩形部分的水平长度、驱动范围边界与水平井之间的压力差、泄油半径、水平井的井筒半径以及水平井与驱动范围下边界之间的距离，因此，终端还可以从水平井的配置信息中获取上述这些参数。

在本公开实施例中，待预测的水平井的驱动范围包括矩形部分和两端部分这两个部分。如图4所示，图4为水平井驱动范围的示意图，该图的方向为垂直于地面且平行于油层剖面的方向。如图5所示，图5为水平井驱动范围中两端部分的示意图，该图的方向为平行于水平井的截面的方向。

需要说明的是，驱动范围中两端部分中每个端部分的形状都为一个半球形，因此，在图4中以垂直于地面且平行于油层剖面的方向观测两端部分时，两端部分呈现半圆形。还需要说明的是，图5中是将驱动范围中的两端部分合并成一个球形，并以平行于水平井的截面的方向观测的。

其中，矩形部分的水平长度与水平井的长度相同，因此，通过对水平井长度的测量可以得到矩形部分的水平长度，该矩形部分的水平长度可以为314.1m。图4中对矩形部分的水平长度进行了注明。驱动范围边界与水平井之间的压力差为驱动范围的边界处的压力值减去水平井所在位置处的压力值，该压力差可以为3mpa。泄油半径是指水平井在平行于地面的水平方向上开采石油的半径，该泄油半径可以为200m，图5中标注了该泄油半径。水平井的井筒半径为水平井的截面半径，该井筒半径可以为0.065m，图5中标注了井筒半径。

需要说明的是，在本公开实施例中，水平井驱动范围中矩形部分的上下边界与该水平井所在油层的上下边界重合，因此，水平井与驱动范围下边界之间的距离也是水平井与油层下边界之间的距离，该水平井与驱动范围下边界之间的距离可以为7m。

在步骤302中，终端确定水平井所在的油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度。

由于水平井所在的油层会影响水平井的产油量，因此，本公开实施例还可以通过终端来确定水平井所在的油层的一些参数，这些参数包括油层的渗透率、油层的粘度、油层的体积系数以及油层的厚度。其中，油层的渗透率可以为600×10^-3μm²，油层的粘度可以为28mpa·s，油层的体积系数可以为1.083，油层的厚度可以为10m。

终端在确定水平井所在的油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度时，在一种可能的实现方式中，终端可以接收用户输入的水平井所在的油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度。在另一种可能的实现方式中，水平井的配置信息中可能包括水平井所在的油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度，因此，终端还可以从水平井的配置信息中获取上述这些参数。

需要说明的是，由于水平井驱动范围中矩形部分的上下边界与该水平井所在油层的上下边界重合，因此，将油层厚度减去水平井与油层下边界之间的距离即为水平井与油层上边界之间的距离，也即水平井与驱动范围上边界之间的距离。也即是，本公开实施例中既存在水平井与驱动范围上边界之间的距离，也存在水平井与驱动范围下边界之间的距离，且这两个距离可以相同，也可以不同。当这两个距离相同时，水平井则处于油层剖面方向上油层的中间位置；当这两个距离相同时，水平井则处于油层剖面方向上油层的非中间位置。因此，本公开实施例中水平井可以位于油层剖面方向上油层的任意一个位置处，适用范围更广。

在步骤303中，终端根据水平长度、压力差、泄油半径、井筒半径、水平井与驱动范围下边界之间的距离、油层的渗透率、油层的粘度、油层的体积系数以及油层的厚度，确定矩形部分对应的第一产油量。

由于水平井的驱动范围中包括矩形部分，也即水平井在矩形部分内也有产油量，因此可以通过步骤301和步骤302中的参数，即水平长度、压力差、泄油半径、井筒半径、水平井与驱动范围下边界之间的距离、油层的渗透率、油层的粘度、油层的体积系数以及油层的厚度，来确定矩形部分对应的第一产油量。

在一种可能的实现方式中，终端可以根据水平长度、压力差、泄油半径、井筒半径、水平井与驱动范围下边界之间的距离、油层的渗透率、油层的粘度、油层的体积系数以及油层的厚度，通过如下公式一确定第一产油量：

公式一：

其中，q1为第一产油量，l为水平长度，δp为压力差，re为泄油半径，rw为井筒半径，zw为水平井与驱动范围下边界之间的距离，k为油层的渗透率，μ0为油层的粘度，b为油层的体积系数，h为油层的厚度，ch为双曲余弦函数，cos为双曲余弦函数，ln为对数函数。

需要说明的是，公式一的确定过程中可以应用镜像反映理论和势的叠加原理，该镜像反映理论即是将位于油层边界附近井的问题转化为无限地层多井同时作用的问题的一个理论。

另外，当矩形部分的水平长度为314.1m，压力差为3mpa，泄油半径为200m，井筒半径为0.065m，水平井与驱动范围下边界之间的距离为7m，油层的渗透率可以600×10^-3μm²，油层的粘度为28mpa·s，油层的体积系数为1.083，油层的厚度为10m时，通过公式一可以计算得到矩形部分对应的第一产油量为154.72m³/d。

在步骤304中，终端根据压力差、泄油半径、井筒半径、油层的渗透率、油层的粘度、油层的体积系数以及油层的厚度，确定两端部分对应的第二产油量。

由于水平井的驱动范围中还包括两端部分，也即水平井在两端部分内也有产油量，因此也可以通过步骤301和步骤302中的参数，即压力差、泄油半径、井筒半径、油层的渗透率、油层的粘度、油层的体积系数以及油层的厚度，来确定两端部分对应的第二产油量。

在一种可能的实现方式中，终端可以根据压力差、泄油半径、井筒半径、油层的渗透率、油层的粘度、油层的体积系数以及油层的厚度，通过如下公式二确定第二产油量：

公式二：

其中，q2为第二产油量，δp为压力差，re为泄油半径，rw为井筒半径，k为油层的渗透率，μ0为油层的粘度，b为油层的体积系数，h为油层的厚度，ln为对数函数。

需要说明的是，两端部分中石油的流动形式分为平面径向流和球面向心流两种。依旧如图5所示，中石油从两端部分的边界向水平井井筒圆心中流入的过程中，存在一个过渡边界，石油从两端部分的边界流向过渡边界的流动过程中是以平面径向流的流动形式流动的，石油从过渡边界向水平井井筒圆心的流动过程中是以球面向心流的流动形式流动的。其中，可以将过渡边界与水平井井筒圆心之间的距离称为过渡半径，图5中标注了过渡半径，该过渡半径在本公开实施例中可以取1.5h，h为油层的厚度。

其中，平面径向流是指在同一平面上进行流动，例如，一个圆柱形的横截面是个圆形，从该横截面的边界向该横截面的中心流动的形式即为平面径向流。球面向心流是指从球面的各个方向上向球心流动，例如，一个球体中从球面上向球心流动的形式即为球面向心流。

另外，公式二的确定过程中可以应用势的叠加原理。并且，当矩形部分的水平长度为314.1m，压力差为3mpa，泄油半径为200m，井筒半径为0.065m，水平井与驱动范围下边界之间的距离为7m，油层的渗透率可以600×10^-3μm²，油层的粘度为28mpa·s，油层的体积系数为1.083，油层的厚度为10m时，通过公式二可以计算得到两端部分对应的第二产油量为14.66m³/d。

在步骤305中，终端确定水平井的渗透率变异系数和流度比，并根据渗透率变异系数和流度比确定油层的波及系数，该波及系数用于描述油层的非均质性对水平井产油量的影响程度，该波及系数与影响程度呈正比。

在实际应用当中，油层的非均质性是普遍存在的，并且油层的非均质性的存在会严重影响水平井的产油量。水平井的长度越长，油层的非均质性对水平井产油量的影响就越明显。因此，本公开实施例引入了波及系数，该波及系数即用于描述油层的非均质性对水平井产油量的影响程度；并且，该波及系数与影响程度呈正比，即波及系数越大，油层的非均质性对水平井产油量的影响越大，波及系数越小，油层的非均质性对水平井产油量的影响越小。

在一种可能的实现方式中，本步骤可以通过以下步骤(1)至(3)实现，包括：

(1)：终端确定水平井的渗透率变异系数。

第一种实现方式，终端获取多个样品水平井的渗透率值，根据每个样品水平井的渗透率值，确定渗透率平均值，根据样品水平井的数量、每个样品水平井的渗透率值和该渗透率平均值，通过以下公式四，确定水平井的渗透率变异系数。

公式四：

其中，ki为第i个样品水平井的渗透率值，为渗透率平均值，n为样品水平井的数量。

第二种实现方式，终端可以在对水平井产油量进行预测之前，按照第一种实现方式获取并存储水平井的渗透率变异系数。在对水平井产油量进行预测时，直接获取已存储的水平井的渗透率变异系数，不需要重复确定水平井的渗透率变异系数，提高了处理效率。

(2)：终端确定水平井的流度比。

终端确定水平井在见水时油水前缘后的平均含水饱和度对应的水的相对渗透率、油水前缘前的平均含水饱和度对应的油的相对渗透率、地层中水的粘度和油层的粘度；根据见水时油水前缘后的平均含水饱和度对应的水的相对渗透率、油水前缘前的平均含水饱和度对应的油的相对渗透率、底层中水的粘度和油层的粘度，通过以下公式五，确定水平井的流度比。

公式五：

其中，krw为见水时油水前缘后的平均含水饱和度对应的水的相对渗透率，kr0为油水前缘前的平均含水饱和度对应的油的相对渗透率，μw为底层中水的粘度，μ0为油层的粘度。

(3)：终端根据该渗透率变异系数和流度比，通过如下公式三确定波及系数：

公式三：

其中，c为波及系数，γ为渗透率变异系数，m为流度比。

在步骤306中，终端确定第一产油量和第二产油量之和，得到第四产油量。

终端在确定第一产油量和第二产油量之后，可以将第一产油量与第二产油量相加，得到第四产油量，该第四产油量为在油层是均质性的情况下水平井的产油量。

在一种可能的实现方式中，终端可以通过如下公式六确定第一产油量和第二产油量之和，得到第四产油量：

公式六：q4＝q1+q2

其中，q4为第四产油量。

需要说明的是，当矩形部分对应的第一产油量为154.72m³/d，两端部分对应的第二产油量为14.66m³/d时，第四产油量为169.38m³/d。

另外，可以先执行步骤305再执行步骤306，也可以先执行步骤306再执行步骤305。也即是，可以先确定水平井的渗透率变异系数和流度比，并根据渗透率变异系数和流度比确定油层的波及系数，再确定第一产油量和第二产油量之和，得到第四产油量的动作；也可以先确定第一产油量和第二产油量之和，得到第四产油量，然后确定水平井的渗透率变异系数和流度比，并根据渗透率变异系数和流度比确定油层的波及系数，本公开实施例对此不做限定。

在步骤307中，终端确定第四产油量与波及系数的乘积，得到第三产油量。

在确定油层的波及系数以及在油层是均质性的情况下水平井的第四产油量后，可以将第四产油量与波及系数相乘，得到第四产油量与波及系数的乘积，该乘积即为水平井对应的第三产油量。

在一种可能的实现方式中，终端可以通过如下公式七确定第四产油量与波及系数的乘积，得到第三产油量：

所述第一产油量满足：

公式七：q3＝q4×c

其中，q3为第三产油量。

另外，当第四产油量为169.38m³/d，波及系数为70.5％时，得到第三产油量为119.42m³/d。该119.42m³/d是通过矩形部分的水平长度为314.1m，压力差为3mpa，泄油半径为200m，井筒半径为0.065m，水平井与驱动范围下边界之间的距离为7m，油层的渗透率可以600×10^-3μm²，油层的粘度为28mpa·s，油层的体积系数为1.083，油层的厚度为10m等计算得到的。而在实际生产中，该待预测的水平井的时机产油量是115.32m³/d，因此，通过本公开实施例的方法得到的水平井的第三产油量与实际产油量的误差为3.55％，该误差小于5％。

在本公开实施例中，确定待预测的水平井的驱动范围中矩形部分的水平长度、驱动范围边界与水平井之间的压力差、泄油半径、水平井的井筒半径、水平井与驱动范围下边界之间的距离、水平井所在的油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度之后，可以确定矩形范围中矩形部分对应的第一产油量以及两端部分对应的第二产油量；然后在确定油层的波及系数后，可以根据矩形部分的第一产油量、两端部分的第二产油量和油层的波及系数，确定水平井对应的第三产油量，提高了对水平井产油量进行预测的精确度。

图6是根据一示例性实施例示出的一种水平井产油量的预测装置600的框图。如图6所示，该装置包括第一确定模块601，第二确定模块602、第三确定模块603、第四确定模块604、第五确定模块605、第六确定模块606。

第一确定模块601，用于确定待预测的水平井的驱动范围中矩形部分的水平长度、该驱动范围边界与该水平井之间的压力差、泄油半径、该水平井的井筒半径以及该水平井与该驱动范围下边界之间的距离，该驱动范围包括该矩形部分和两端部分

第二确定模块602，用于确定该水平井所在的油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度；

第三确定模块603，用于根据该水平长度、该压力差、该泄油半径、该井筒半径、该水平井与该驱动范围下边界之间的距离、该油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度，确定该矩形部分对应的第一产油量；

第四确定模块604，用于根据该压力差、该泄油半径、该井筒半径、该油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度，确定该两端部分对应的第二产油量；

第五确定模块605，用于确定该水平井的渗透率变异系数和流度比，并根据该渗透率变异系数和该流度比确定该油层的波及系数，该波及系数用于描述该油层的非均质性对该水平井产油量的影响程度，该波及系数与该影响程度呈正比；

第六确定模块606，用于根据该第一产油量、该第二产油量和该波及系数，确定该水平井对应的第三产油量。

在一个可能的实现方式中，该第三确定模块603包括：

第一确定子模块，用于根据该水平长度、该压力差、该泄油半径、该井筒半径、该水平井与该驱动范围下边界之间的距离、该油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度，通过如下公式一确定该第一产油量：

公式一：

其中，该q1为该第一产油量，该l为该水平长度，该δp为该压力差，该re为该泄油半径，该rw为该井筒半径，该zw为该水平井与该驱动范围下边界之间的距离，该k为该油层的渗透率，该μ0为该油层的粘度，该b为该油层的体积系数，该h为该油层的厚度，该ch为双曲余弦函数，该cos为双曲余弦函数，该ln为对数函数。

在一个可能的实现方式中，该第四确定模块604包括：

第二确定子模块，用于根据该压力差、该泄油半径、该井筒半径、该油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度，通过如下公式二确定该第二产油量：

公式二：

其中，该q2为该第二产油量，该δp为该压力差，该re为该泄油半径，该rw为该井筒半径，该k为该油层的渗透率，该μ0为该油层的粘度，该b为该油层的体积系数，该h为该油层的厚度，该ln为对数函数。

在一个可能的实现方式中，该第五确定模块605包括：

第三确定子模块，用于根据该渗透率变异系数和该流度比，通过如下公式三确定该波及系数：

公式三：

其中，该c为该波及系数，该γ为该渗透率变异系数，该m为该流度比。

在一个可能的实现方式中，该第六确定模块606包括：

第四确定子模块，用于确定该第一产油量和该第二产油量之和，得到第四产油量；

第五确定子模块，用于确定该第四产油量与该波及系数的乘积，得到该第三产油量。

在本公开实施例中，确定待预测的水平井的驱动范围中矩形部分的水平长度、驱动范围边界与水平井之间的压力差、泄油半径、水平井的井筒半径、水平井与驱动范围下边界之间的距离、水平井所在的油层的渗透率、该油层的粘度、该油层的体积系数以及该油层的厚度之后，可以确定矩形范围中矩形部分对应的第一产油量以及两端部分对应的第二产油量；然后在确定油层的波及系数后，可以根据矩形部分的第一产油量、两端部分的第二产油量和油层的波及系数，确定水平井对应的第三产油量，提高了对水平井产油量进行预测的精确度。

关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。

图7示出了本发明一个示例性实施例提供的终端700的结构框图。该终端700可以是：智能手机、平板电脑、笔记本电脑或台式电脑。终端700还可能被称为用户设备、便携式终端、膝上型终端、台式终端等其他名称。

通常，终端700包括有：处理器701和存储器702。

处理器701可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器701可以采用dsp(digitalsignalprocessing，数字信号处理)、fpga(field－programmablegatearray，现场可编程门阵列)、pla(programmablelogicarray，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器701也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称cpu(centralprocessingunit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器701可以在集成有gpu(graphicsprocessingunit，图像处理器)，gpu用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器701还可以包括ai(artificialintelligence，人工智能)处理器，该ai处理器用于处理有关机器学习的计算操作。

存储器702可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器702还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器702中的非暂态的计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器701所执行以实现本申请中方法实施例提供的水平井产油量的预测方法。

在一些实施例中，终端700还可选包括有：外围设备接口703和至少一个外围设备。处理器701、存储器702和外围设备接口703之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口703相连。具体地，外围设备包括：射频电路704、触摸显示屏705、摄像头706、音频电路707、定位组件708和电源709中的至少一种。

外围设备接口703可被用于将i/o(input/output，输入/输出)相关的至少一个外围设备连接到处理器701和存储器702。在一些实施例中，处理器701、存储器702和外围设备接口703被集成在同一芯片或电路板上；在一些其他实施例中，处理器701、存储器702和外围设备接口703中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。

射频电路704用于接收和发射rf(radiofrequency，射频)信号，也称电磁信号。射频电路704通过电磁信号与通信网络以及其他通信设备进行通信。射频电路704将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路704包括：天线系统、rf收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路704可以通过至少一种无线通信协议来与其它终端进行通信。该无线通信协议包括但不限于：城域网、各代移动通信网络(2g、3g、4g及5g)、无线局域网和/或wifi(wirelessfidelity，无线保真)网络。在一些实施例中，射频电路704还可以包括nfc(nearfieldcommunication，近距离无线通信)有关的电路，本申请对此不加以限定。

显示屏705用于显示ui(userinterface，用户界面)。该ui可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏705是触摸显示屏时，显示屏705还具有采集在显示屏705的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器701进行处理。此时，显示屏705还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏705可以为一个，设置终端700的前面板；在另一些实施例中，显示屏705可以为至少两个，分别设置在终端700的不同表面或呈折叠设计；在再一些实施例中，显示屏705可以是柔性显示屏，设置在终端700的弯曲表面上或折叠面上。甚至，显示屏705还可以设置成非矩形的不规则图形，也即异形屏。显示屏705可以采用lcd(liquidcrystaldisplay，液晶显示屏)、oled(organiclight-emittingdiode,有机发光二极管)等材质制备。

摄像头组件706用于采集图像或视频。可选地，摄像头组件706包括前置摄像头和后置摄像头。通常，前置摄像头设置在终端的前面板，后置摄像头设置在终端的背面。在一些实施例中，后置摄像头为至少两个，分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种，以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及vr(virtualreality，虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中，摄像头组件706还可以包括闪光灯。闪光灯可以是单色温闪光灯，也可以是双色温闪光灯。双色温闪光灯是指暖光闪光灯和冷光闪光灯的组合，可以用于不同色温下的光线补偿。

音频电路707可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波，并将声波转换为电信号输入至处理器701进行处理，或者输入至射频电路704以实现语音通信。出于立体声采集或降噪的目的，麦克风可以为多个，分别设置在终端700的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器701或射频电路704的电信号转换为声波。扬声器可以是传统的薄膜扬声器，也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时，不仅可以将电信号转换为人类可听见的声波，也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中，音频电路707还可以包括耳机插孔。

定位组件708用于定位终端700的当前地理位置，以实现导航或lbs(locationbasedservice，基于位置的服务)。定位组件708可以是基于美国的gps(globalpositioningsystem，全球定位系统)、中国的北斗系统、俄罗斯的格雷纳斯系统或欧盟的伽利略系统的定位组件。

电源709用于为终端700中的各个组件进行供电。电源709可以是交流电、直流电、一次性电池或可充电电池。当电源709包括可充电电池时，该可充电电池可以支持有线充电或无线充电。该可充电电池还可以用于支持快充技术。

在一些实施例中，终端700还包括有一个或多个传感器710。该一个或多个传感器710包括但不限于：加速度传感器711、陀螺仪传感器712、压力传感器713、指纹传感器714、光学传感器715以及接近传感器716。

加速度传感器711可以检测以终端700建立的坐标系的三个坐标轴上的加速度大小。比如，加速度传感器711可以用于检测重力加速度在三个坐标轴上的分量。处理器701可以根据加速度传感器711采集的重力加速度信号，控制触摸显示屏705以横向视图或纵向视图进行用户界面的显示。加速度传感器711还可以用于游戏或者用户的运动数据的采集。

陀螺仪传感器712可以检测终端700的机体方向及转动角度，陀螺仪传感器712可以与加速度传感器711协同采集用户对终端700的3d动作。处理器701根据陀螺仪传感器712采集的数据，可以实现如下功能：动作感应(比如根据用户的倾斜操作来改变ui)、拍摄时的图像稳定、游戏控制以及惯性导航。

压力传感器713可以设置在终端700的侧边框和/或触摸显示屏705的下层。当压力传感器713设置在终端700的侧边框时，可以检测用户对终端700的握持信号，由处理器701根据压力传感器713采集的握持信号进行左右手识别或快捷操作。当压力传感器713设置在触摸显示屏705的下层时，由处理器701根据用户对触摸显示屏705的压力操作，实现对ui界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。

指纹传感器714用于采集用户的指纹，由处理器701根据指纹传感器714采集到的指纹识别用户的身份，或者，由指纹传感器714根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时，由处理器701授权该用户执行相关的敏感操作，该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器714可以被设置终端700的正面、背面或侧面。当终端700上设置有物理按键或厂商logo时，指纹传感器714可以与物理按键或厂商logo集成在一起。

光学传感器715用于采集环境光强度。在一个实施例中，处理器701可以根据光学传感器715采集的环境光强度，控制触摸显示屏705的显示亮度。具体地，当环境光强度较高时，调高触摸显示屏705的显示亮度；当环境光强度较低时，调低触摸显示屏705的显示亮度。在另一个实施例中，处理器701还可以根据光学传感器715采集的环境光强度，动态调整摄像头组件706的拍摄参数。

接近传感器716，也称距离传感器，通常设置在终端700的前面板。接近传感器716用于采集用户与终端700的正面之间的距离。在一个实施例中，当接近传感器716检测到用户与终端700的正面之间的距离逐渐变小时，由处理器701控制触摸显示屏705从亮屏状态切换为息屏状态；当接近传感器716检测到用户与终端700的正面之间的距离逐渐变大时，由处理器701控制触摸显示屏705从息屏状态切换为亮屏状态。

本领域技术人员可以理解，图7中示出的结构并不构成对终端700的限定，可以包括比图示更多或更少的组件，或者组合某些组件，或者采用不同的组件布置。

本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本发明的其它实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本公开未公开的本技术领域中的公知常识或惯用技术手段。

应当理解的是，本发明并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。