确定油藏注水井窜流量的方法及装置与流程

本申请涉及油藏工程动态分析方法，具体而言，涉及一种确定油藏注水井窜流量的方法及装置。

背景技术：

油田开发中后期，受井网形式、注采配置关系、平面非均质性等因素影响，各井组之间注采对应程度差别很大，部分井组注采极不均衡，出现平面干扰严重、水淹、水窜等现象，使注入水低效、无效循环，降低注水利用效率。

确定储层中注入水的窜流状况是制定水井调剖对策的重要依据。国外对于油藏窜流状况的研究起步较早，主要是针对由于成岩作用形成的高渗层和裂缝而发育的水驱原生窜流通道；国内的研究发展比较迅速且成熟，包括生产动态分析、试井、测井资料识别、示踪剂监测，以及基于层次分析、模糊综合评判法的数学理论体系。这些技术主要针对陆相砂岩油藏水驱形成的次生高渗通道，可以对窜流通道的平面、纵向、井间发育状况进行描述，但现有技术缺乏针对注水井的窜流状况如窜流量的定量研究方法，不能满足油藏后期精细调剖的技术要求。

技术实现要素：

本申请的主要目的在于提供一种确定油藏注水井窜流量的方法，以解决现有技术中缺乏针对注水井的窜流状况进行定量研究的问题。

为了实现上述目的，根据本申请的一个方面，提供了一种确定油藏注水井窜流量的方法，该方法包括：

获取预设时间周期内各井组中每口注水井对油井的注水量、各井组中油井的产液量、以及各注水井的总注水量；

根据所述注水量以及所述产液量计算出各井组中每口注水井对油井的贡献值；

根据所述产液量计算出各井组中油井的需求分配总水量；

根据所述贡献值以及所述需求分配总水量计算出各井组中每口注水井对油井的分配水量；

根据所述分配水量以及所述总注水量计算出各注水井的窜流量。

为了实现上述目的，根据本申请的另一方面，提供了一种确定油藏注水井窜流量的装置，该装置包括：

数据获取模块，用于获取预设时间周期内各井组中每口注水井对油井的注水量、各井组中油井的产液量、以及各注水井的总注水量；

贡献值计算模块，用于根据所述注水量以及所述产液量计算出各井组中每口注水井对油井的贡献值；

需求分配总水量计算模块，用于根据所述产液量计算出各井组中油井的需求分配总水量；

分配水量计算模块，用于根据所述贡献值以及所述需求分配总水量计算出各井组中每口注水井对油井的分配水量；

窜流量计算模块，用于根据所述分配水量以及所述总注水量计算出各注水井的窜流量。

为了实现上述目的，根据本申请的另一方面，还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述确定油藏注水井窜流量的方法中的步骤。

为了实现上述目的，根据本申请的另一方面，还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序在计算机处理器中执行时实现上述确定油藏注水井窜流量的方法中的步骤。

本申请的有益效果为：本申请利用静态地质资料和井网生产动态数据相结合的方法，可以较准确地确定水窜的注水井，并定量化计算出注水井的窜流量，可以为油藏后期精细调剖提供数据依据和技术指导。

附图说明

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

图1是本申请实施例确定油藏注水井窜流量方法的流程图；

图2是本申请实施例油藏中区块a井网示意图；

图3是本申请实施例计算注水井贡献值的流程示意图；

图4是本申请实施例计算各井组中注水井对油井的贡献值的流程示意图；

图5是本申请实施例计算各注水井窜流量的流程示意图；

图6是本申请实施例计算油井需求分配总水量的流程示意图；

图7是本申请实施例确定油藏注水井窜流量的装置的结构框图；

图8是本申请实施例贡献值计算模块的结构框图；

图9是本申请实施例关联度计算单元的结构框图；

图10是本申请实施例窜流量计算模块的结构框图；

图11是本申请实施例分配水量计算模块的结构框图；

图12是本申请实施例油藏中区块b油水井注采对应状况分析表；

图13是本申请实施例油藏中区块b注水井窜流量及窜流级别统计表；

图14是本申请实施例油藏中区块b注水井窜流状况统计表。

具体实施方式

为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分的实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。

需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。

需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。

需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“井组”又称为“注采井组”，其中每个“井组”中包括一个油井和至少一个向该油井注水的注水井。此外，一个藏注的井网中可以同时包括多个井组，每个注水井根据其对应的油井，也可以同时属于两个及以上的井组。

如图1所示，本申请实施例的确定油藏注水井窜流量方法包括步骤s101至步骤s105。

步骤s101，获取预设时间周期内各井组中每口注水井对油井的注水量、各井组中油井的产液量、以及各注水井的总注水量。由于在本领域中通常使用月注水量以及月产液量等来描述注水井和油井的生产参数，因此在本实施例中，上述预设时间周期内的注水量、产液量以及总注水量，可以分别为：各井组中每口注水井对油井的月注水量、各井组中油井的月产液量以及各注水井的月总注水量。

在本实施例中，井网中包括多个井组，每个井组包括一个油井以及至少一个注水井，本申请在计算注水井的窜流量之前，需要先获取每个井组中每口注水井对油井的月注水量、每个井组中油井的月产液量、以及该井网中所有井组所包含的每个注水井的月总注水量。

图2是本申请实施例油藏中区块a的井网示意图，如图2所示，该井网中共有3口油井分别为a、b和c，该井网中共有6口注水井分别为1、2、3、4、5和6，该井网中共有3个井组，分别为：由油井a和向油井a注水的注水井1、2和3组成的第一井组；由油井b以及向油井b注水的注水井2、3、5和6组成的第二井组；油井c以及向油井c注水的注水井3、4和5组成的第三井组。

下面结合图2对本步骤进行说明，在计算油藏中区块a的井网中注水井的窜流量之前，需要先获取：每个井组中每口注水井对油井的月注水量，例如第一井组，需要分别获取第一井组中注水井1、2和3向油井a的月注水量；每个井组中油井的月产液量，即油井a、b和c的月产液量；以及该井网中所有井组所包含的每个注水井的月总注水量，即分别获取注水井1、2、3、4、5和6的月总注水量。

步骤s102，根据所述注水量以及所述产液量计算出各井组中每口注水井对油井的贡献值。在本步骤中，计算出各井组中每口注水井对油井的贡献值，在本实施例中，该贡献值可以为贡献比。在本实施例中该贡献比的计算方法为：先根据井组中每口注水井的月注水量以及油井的月产液量，分别生成每口注水井的月注水量序列以及油井的月产液量序列，进而根据月注水量序列和月产液量序列，通过灰色关联度算法，计算出该井组中每口注水井对油井的关联度，而油井的贡献比为该油井的关联度与该井组中所有油井的关联度之和的比值。

下面结合图2对本步骤进行说明，在本步骤中计算出图2所示的井网中每个井组中每口注水井对油井的贡献值，例如上述的第一井组，分别计算出第一井组中注水井1、2和3对油井a的贡献值。

步骤s103，根据所述产液量计算出各井组中油井的需求分配总水量。在本步骤中，根据各井组中油井的产液量，分别计算出需要分配给各油井的总水量，即各油井的需求分配总水量，该需求分配总水量可以表示油井出产一定量的油液而需要分配的水量，通常可以通过分析油井的月产液量和动液面变化趋势，并根据物质平衡的基本原理而得出。在本实施例中，油井的月需求分配总水量qw'为该油井的月产液量ql与分配系数ε的乘积，公式为：

qw'＝ε×ql(1)

该分配系数ε为本领域的工作人员参照油水井各韵律层注采对应状况并结合工作经验而确定的一个常数，一般该分配系数ε的取值区间为0.8至1.2。

步骤s104，根据所述贡献值以及所述需求分配总水量计算出各井组中每口注水井对油井的分配水量。在本步骤中，根据在步骤s102中计算出的贡献值以及通过步骤s103中计算出的油井需求分配总水量，计算出各井组中每口注水井对油井的分配水量。在本实施例中，井组中每口注水井对油井的分配水量为每口注水井对油井的贡献比与油井的需求分配总水量的乘积。在本实施例中，如果一个注水井同时属于两个及以上个井组时，则会分别计算出该注水井对每个井组中的油井的分配水量。

下面结合图2对本步骤进行说明，在本步骤中计算出图2所示的井网中每个井组中每口注水井对油井的分配水量，以上述第一井组为例，根据在上述步骤s102中分别计算得出的注水井1、2和3对油井a的贡献值，以及在步骤s103中计算得出的油井a的需求分配总水量，可以分别计算出注水井1、2和3对油井a的分配水量，同理，可以计算得出第二井组中注水井2、3、5和6对油井b的分配水量，以及第三井组中注水井3、4和5对油井c的分配水量。如图2所示，由于注水井2同时属于第一井组和第二井组，即同时向油井a和油井b注水，则在本步骤中会分别计算出注水井2对油井a和油井b的分配水量。

步骤s105，根据所述分配水量以及所述总注水量计算出各注水井的窜流量。在本申请中，注水井的窜流量为该注水井的总注水量与总分配水量的差。在本实施例中，在计算某注水井的窜流量时，需要先得到该注水井的总分配水量，当该油井同时属于两个及以上个井组时，即该油井同时向两个及以上个油井注水时，需要根据在上述步骤s104中计算得出的该油井分别对每个油井的分配水量，求和得出该油井的总分配水量，进而将该注水井的总注水量与总分配水量求差，得到该注水井的窜流量。在本步骤中，如果一个注水井只属于一个井组时，则该注水井的总分配水量就为该注水井对该井组的油井的分配水量。

下面结合图2对本步骤进行说明，在本步骤中，在计算注水井的窜流量时，需要先得到该注水井的总分配水量，如图2所示的注水井3，该注水井3同时向油井a、b和c注水，通过在上述步骤s104中计算得出的注水井3分别对油井a、b和c的分配水量，加和得到注水井3的总分配水量，进而将注水井3的总注水量与总分配水量求差，得到该注水井3的窜流量。如图2所示的注水井1，由于该注水井1只属于第一井组，即只为油井a注水，则注水井1的总分配水量就为注水井1对油井a的分配水量。

从以上的描述中，可以看出，本申请利用静态地质资料和井网生产动态数据相结合的方法，可以较准确地确定水窜的注水井，并定量化计算出注水井的窜流量，可以为油藏后期精细调剖提供数据依据和技术指导。

如图3所示，上述步骤s102，根据所述注水量以及所述产液量计算出各井组中每口注水井对油井的贡献值，具体包括步骤s201至步骤s202。

步骤s201，根据所述注水量以及所述产液量计算出各井组中每口注水井对油井的关联度，在本申请的实施例中，可以通过灰色关联度算法，计算出井组中每口注水井对油井的灰色关联度。

步骤s202，根据所述关联度计算出各井组中每口注水井对油井的贡献值，在本步骤中，根据井组中每口注水井对油井的灰色关联度计算出每口注水井对油井的贡献值。在本实施例中，该贡献值为贡献比，油井的贡献比为该注水井的灰色关联度与该井组中所有注水井的灰色关联度之和的比值。

图4为本申请实施例的根据所述注水量以及所述产液量计算出各井组中每口注水井对油井的贡献值的方法流程图，如图4所示，该方法包括步骤s301至步骤s304。

步骤s301，根据所述注水量生成各井组中每口注水井的周期注水量序列。在本步骤中，根据在上述步骤s101中得到的各井组中每口注水井对油井的月注水量，生成各井组中每口注水井的周期注水量序列，作为后续计算的自变量序列。在本实施例中，某井组中每口注水井的周期注水量序列可以为xi’＝(xi’(1)，xi’(2)，…，xi’(n))^t，其中i＝1,2,…,m，m为该井组中注水井的数量，xi’表示该井组中第i口注水井的周期注水量序列，n为周期注水量序列中的数据的个数。

步骤s302，根据所述产液量生成各井组中油井的周期产液量序列。在本步骤中，根据在上述步骤s101中得到的各井组中油井的月产液量，生成各井组中油井的周期产液量序列，作为后续计算的因变量序列。在本实施例中，某井组中油井的周期产液量序列可以为x0’＝(x0’(1),x0’(2),…，x0’(n))^t，其中n为周期注水量序列中的数据的个数。

步骤s303，根据所述周期注水量序列以及周期产液量序列，通过灰色关联度算法，计算出各井组中每口注水井对油井的关联度。在本申请中，通过灰色关联度的算法来计算每个井组中注水井对油井的关联度。在本实施例中，根据灰色关联度分析原理，首先将在上述步骤s301中得到的周期注水量序列和在上述步骤s302中得到的周期注水量序列作为分析的基础序列，以n×(m+1)矩阵形式表示为：

在公式(2)中：x0’＝(x0’(1),x0’(2),…，x0’(n))^t表示某井组中由油井月产液量构成的因变量序列；xi’＝(xi’(1)，xi’(2)，…，xi’(n))^t表示某井组中由注水井月注水量构成的自变量序列，其中i＝1,2,…,m，m为该井组中注水井的数量，xi’表示该井组中第i口注水井的周期注水量序列，n为序列中的数据的个数。

由于原始变量序列具有不同的量纲或数量级，为了保证分析结果的可靠性，进一步对变量序列进行无量纲化处理得到如下矩阵：

在公式(3)中：x0＝(x0(1),x0(2),…，x0(n))^t表示无量纲化因变量序列；

xi＝(xi(1)，xi(2)，…，xi(n))^t，i＝1,2,…,m表示无量纲化自变量序列。

上述无量纲化采用归一化方法：

在公式(4)中：i＝0,1,2…m，k＝1,2…n。

在公式(3)中第一列(因变量序列)与其余各列(自变量序列)对应元素的绝对差值，形成如下绝对差值矩阵：

其中：

δi(k)＝|x0(k)-xi(k)|(6)

在公式(6)中：i＝1,2…m，k＝1,2…n。

将δ(max)、δ(max)分别定义为绝对差值矩阵中的最大数和最小数：

根据公式(7)和(8)对绝对差值阵中数据作如下变换得到关联系数：

由公式(9)得到关联系数矩阵：

在公式(9)中：ρ为分辨系数，表示关联系数间的差异，根据式(10)中数据分布情况多在0.1～0.5之间取值，一组参数对应一个ρ值，一般取0.1，以提高因素之间的关联度差异；ξi(k)为关联系数，是不超过1的正数，δi(k)越小，ξi(k)越大，它反映自变量序列xi与因变量序列x0的关联程度。

对式(10)中第i列的n个关联系数求平均得到灰色关联度：

进而根据公式(11)求得该井组注水井的灰色关联度序列：

γ＝{γ1,γ2,…,γm}(12)

在公式(12)中，γ为灰色关联度，m为该井组中注水井的数量，γm为该井组中第m个注水井对井组中油井的灰色关联度。

步骤s304，根据所述关联度计算出各井组中每口注水井对油井的贡献值。在本步骤中，可以根据在上述步骤s303中得到的各井组注水井的灰色关联度序列，计算出各井组中每口注水井对油井的贡献值。井组中每口注水井对油井的贡献值序列为：

g＝{g1,g2,…,gm}(13)

其中，贡献值序列(13)为根据上述灰色关联度序列(12)计算得出，具体计算公式为：

在上述公式(14)和(15)中的ε为在上述公式(1)中提到的分配系数，取值区间为0.8至1.2。

在公式(14)中，gi为井组中第i个注水井对油井的贡献值。

在通过公式(13)至(15)得到各井组中每口注水井对油井的贡献值之后，上述步骤s104，根据所述贡献值以及所述需求分配总水量计算出各井组中每口注水井对油井的分配水量，就可以用如下公式来表示：

qwi＝gi·ql(16)

其中，gi为井组中第i口注水井对油井的贡献值，ql为该井组油井的需求分配总水量，qwi即为井组中第i口注水井对油井的分配水量。

根据公式(16)可以得出上述井组中第i口注水井向与其对应的所有油井的分配水量的数据序列qwi,j，其中j为与上述第i口注水井对应的第j口油井。

在通过公式(16)得到各井组中每口注水井对油井的分配水量之后，进而上述步骤s105，根据所述分配水量以及所述总注水量计算出各注水井的窜流量，就可以用如下公式来表示：

在公式(17)中：qrwi为第i口水井窜流量；qawi为第i口注水井总注水量；k为与该注水井对应的油井数量。以此可以计算出各井组中每口注水井的窜流量。

如图5所示，上述步骤s105，根据所述分配水量以及所述总注水量计算出各注水井的窜流量，具体包括步骤s401和步骤s402。

步骤s401，根据所述分配水量计算出各注水井的总分配水量。在本实施例中，在计算某注水井的窜流量时，需要先得到该注水井的总分配水量，当该油井同时属于两个及以上个井组时，即该注水井同时向两个及以上个油井注水时，需要根据在上述步骤s104中计算得出的该注水井分别对每个油井的分配水量，求和得出该注水井的总分配水量。如果一个注水井只属于一个井组时，则该注水井的总分配水量就为该注水井对该井组的油井的分配水量。在得到注水井的总分配水量之后，进而通过步骤s402，计算所述总注水量与所述总分配水量之差，得到各注水井的窜流量。

如图6所示，上述步骤s103，根据所述产液量计算出各井组中油井的需求分配总水量，具体包括步骤s501和步骤s502。

步骤s501，获取预设的各井组中油井的分配系数。在本步骤中，由于油井的需求分配总水量qw'为该油井的产液量ql与分配系数ε的乘积，因此在计算各井组中油井的需求分配总水量时，需要先获取设定的各井组中油井的分配系数ε。该分配系数ε为本领域的工作人员参照油水井各韵律层注采对应状况并结合工作经验而确定的一个常数，一般该分配系数ε的取值区间为0.8至1.2。进而在步骤s502中,通过计算各井组中油井的产液量与分配系数的乘积，得到各井组中油井的需求分配总水量。

在本申请的一可选实施例中，在上述步骤s403中计算出各井组中注水井的关联度序列γ＝{γ1,γ2,…,γm}之后，还根据预设的对应关系，确定了井组中注水井对油井的关联性级别，其中关联性级别由强到弱依次分为a：0.75≤γ＜1.00，b：0.50≤γ＜0.75，c：0.25≤γ＜0.50，d：0.00＜γ＜0.25。

在本申请的一可选实施例中，在通过上述步骤s404中的公式(17)计算出注水井窜流量qrwi之后，还根据窜流量与窜流级别的对应关系确定出各注水井的窜流级别，其中根据生产数据统计规律，将注水井窜流程度由强到弱划分为四个级别，a：5≤qrwi，b：3≤qrwi＜5，c：1≤qrwi＜3，d：0≤qrwi＜1，(其中qrwi单位为10⁴m³)。

在本申请的一可选实施例中，还根据计算出的注水井的窜流量以及注水井的总注水量，来计算出该注水井的窜流比例，其中窜流比例为窜流量与总注水量的比值。

图12为本申请实施例油藏中区块b油水井注采对应状况分析表，如图12所示，油藏区块b中的油井p1142同时对应四口注水井，其中注水井p1150与该油井的关联性级别最高；而油井p1170只对应注水井p1171，而注水井p1171对该油井的关联性级别达到a级。

图13为本申请实施例油藏中区块b注水井窜流量及窜流级别统计表，如图13所示，在油藏区块b中注水井p1145和p1132的窜流量和窜流级别相对较高，同时其窜流比例也相对较高，需要有针对性的采取注水井调剖措施。

图14为本申请实施例油藏中区块b注水井窜流状况统计表，从图13中可以看出，油藏区块a的窜流状况比较严重，13口注水井中窜流井有8口，其中a、b、c、d级窜流井分别为2口、6口、0口、0口，分别占目前总水井数的15.4％、46.2％、0％、0％，需要有针对性的采取注水井调剖措施。

需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。

基于同一发明构思，本申请实施例还提供了一种确定油藏注水井窜流量的装置，可以用于实现上述实施例所描述的方法，如下面的实施例所述。由于确定油藏注水井窜流量的装置解决问题的原理与确定油藏注水井窜流量的方法相似，因此确定油藏注水井窜流量的装置的实施可以参见确定油藏注水井窜流量的方法的实施例，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

如图7所示，该装置包括：数据获取模块1、贡献值计算模块2、需求分配总水量计算模块3、分配水量计算模块4以及窜流量计算模块5。

数据获取模块1，用于获取预设时间周期内各井组中每口注水井对油井的注水量、各井组中油井的产液量、以及各注水井的总注水量。由于在本领域中通常使用月注水量以及月产液量等来描述注水井和油井的生产参数，因此在本申请的实施例中，上述预设时间周期内的注水量、产液量以及总注水量，分别为：各井组中每口注水井对油井的月注水量、各井组中油井的月产液量以及各注水井的月总注水量。本申请在计算注水井的窜流量之前，需要先通过数据获取模块1获取每个井组中每口注水井对油井的月注水量、每个井组中油井的月产液量、以及该井网中所有井组所包含的每个注水井的月总注水量。

贡献值计算模块2，用于根据所述注水量以及所述产液量计算出各井组中每口注水井对油井的贡献值。在本申请的可选实施例中，贡献值的计算方法为：先根据井组中每口注水井的月注水量以及油井的月产液量，分别生成每口注水井的月注水量序列以及油井的月产液量序列，进而根据月注水量序列和月产液量序列，通过灰色关联度算法，计算出该井组中每口注水井对油井的关联度，而油井的贡献值为该油井的关联度与该井组中所有油井的关联度之和的比值。

需求分配总水量计算模块3，用于根据所述产液量计算出各井组中油井的需求分配总水量。在本申请实施例中，油井的月需求分配总水量为该油井的月产液量与分配系数的乘积，其中分配系数为本领域的工作人员参照油水井各韵律层注采对应状况并结合工作经验而确定的一个常数，一般该分配系数ε的取值区间为0.8至1.2。

分配水量计算模块4，用于根据所述贡献值以及所述需求分配总水量计算出各井组中每口注水井对油井的分配水量。在本申请实施例中，井组中每口注水井对油井的分配水量为每口注水井对油井的贡献比与油井的需求分配总水量的乘积。在本申请实施例中，如果一个注水井同时属于两个及以上个井组时，则会分别计算出该注水井对每个井组中的油井的分配水量。

窜流量计算模块5，用于根据所述分配水量以及所述总注水量计算出各注水井的窜流量。在本申请中，注水井的窜流量为该注水井的总注水量与总分配水量的差。在本申请实施例中，在计算某注水井的窜流量时，需要先得到该注水井的总分配水量，当该油井同时属于两个及以上个井组时，即该油井同时向两个及以上个油井注水时，需要根据该油井分别对每个油井的分配水量，求和得出该油井的总分配水量，进而将该注水井的总注水量与总分配水量求差，得到该注水井的窜流量。在本申请实施例中，如果一个注水井只属于一个井组时，则该注水井的总分配水量就为该注水井对该井组的油井的分配水量。

如图8所示，上述贡献值计算模块包括2包括：关联度计算单元6和算法单元7。

关联度计算单元6，用于根据所述注水量以及所述产液量计算出各井组中每口注水井对油井的关联度。在本申请的实施例中，关联度计算单元6通过灰色关联度算法，计算出井组中每口注水井对油井的灰色关联度。

算法单元7，用于根据所述关联度计算出各井组中每口注水井对油井的贡献值。在本申请实施例中，算法单元7根据井组中每口注水井对油井的灰色关联度计算出每口注水井对油井的贡献值。在本申请实施例中，该贡献值为贡献比，油井的贡献比为该注水井的灰色关联度与该井组中所有注水井的灰色关联度之和的比值。

如图9所示，上述关联度计算单元6包括：周期注水量序列生成单元601、周期产液量序列生成单元602和灰色关联度计算单元603。

周期注水量序列生成单元601，用于根据所述注水量生成各井组中每口注水井的周期注水量序列。

周期产液量序列生成单元602，用于根据所述产液量生成各井组中油井的周期产液量序列；

灰色关联度计算单元603，用于根据所述周期注水量序列以及周期产液量序列，通过灰色关联度算法，计算出各井组中每口注水井对油井的关联度。在本申请的实施例中，灰色关联度计算单元603可以采用上述步骤s303中的计算公式，来计算出各井组中每口注水井对油井的关联度。

如图10所示，上述窜流量计算模块5包括：总分配水量计算单元501和计算单元502。

总分配水量计算单元501，用于根据所述分配水量计算出各注水井的总分配水量。在本申请实施例中，在计算某注水井的窜流量时，需要先得到该注水井的总分配水量，当该注水井同时属于两个及以上个井组时，即该注水井同时向两个及以上个油井注水时，总分配水量计算单元501获取该注水井分别对每个油井的分配水量，求和得出该油井的总分配水量。如果一个注水井只属于一个井组时，则该注水井的总分配水量就为该注水井对该井组的油井的分配水量。

计算单元502，用于通过计算所述总注水量与所述总分配水量之差，得到各注水井的窜流量。

如图11所示，上述分配水量计算模块4包括：分配系数获取单元401和计算单元402。

分配系数获取单元401，用于获取预设的各井组中油井的分配系数。在本申请中，由于油井的需求分配总水量为该油井的产液量与分配系数的乘积，因此在计算各井组中油井的需求分配总水量时，需要先获取设定的各井组中油井的分配系数。该分配系数为本领域的工作人员参照油水井各韵律层注采对应状况并结合工作经验而确定的一个常数，一般该分配系数的取值区间为0.8至1.2。

计算单元402，用于通过计算所述产液量与所述分配系数的乘积，得到各井组中油井的需求分配总水量。

本申请的另一方面，还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述确定油藏注水井窜流量的方法中的步骤。

本申请的另一方面，还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序在计算机处理器中执行时实现上述确定油藏注水井窜流量的方法中的步骤。

显然，本领域的技术人员应该明白，上述的本申请的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本申请不限制于任何特定的硬件和软件结合。

以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。