具体实施方式

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

本文实施例所述的方法和装置应用于页岩油水平井的重复压裂选井中，由于待选井水平井组中水平井的数量庞大，且对每个井进行重复压裂计算的工作量大，因此将待选井水平井组中的所有水平井均进行重复压裂计算，然后根据计算结果得出能够进行重复压裂的水平井的方法将带来巨大的工作量，影响水平井的生产。本文实施例所述的方法在待选井水平井组中优选出能够进行重复压裂计算的水平井，指导工作人员对优选出的水平井进行重复压裂计算，减小了计算量，并且考虑了井间窜扰对选井的影响，提高了选井精度。

如图1所示为本文实施例一种页岩油水平井重复压裂选井装置的结构示意图，在本图中描述了页岩油水平井重复压裂选井装置的结构，具体包括被窜扰影响产量计算单元101、水平井综合指标计算单元102、选井单元103：

被窜扰影响产量计算单元101获取待选井水平井组中各水平井产生窜扰前的产油水平以及窜扰期间的产油量，计算各水平井的被窜扰影响产量；

水平井综合指标计算单元102根据包括被窜扰影响产量计算单元101得到的被窜扰影响产量的所述各水平井的生产指标，计算各水平井包括被窜扰程度反向百分位数的综合指标；

选井单元103根据水平井综合指标计算单元102得到的各水平井的综合指标以及设定的阈值选出能够进行重复压裂计算的水平井。

如图2所示为文实施例一种页岩油水平井重复压裂选井方法的流程图，在本图中描述了重复压裂选井的过程，该方法包括：

步骤201：获取待选井水平井组中各水平井产生窜扰前的产油水平以及窜扰期间的产油量，计算各水平井的被窜扰影响产量；

步骤202：根据包括所述被窜扰影响产量的所述各水平井的生产指标，计算各水平井包括被窜扰程度反向百分位数的综合指标；

步骤203：根据所述各水平井的综合指标以及设定的阈值选出能够进行重复压裂计算的水平井。

在本文实施例中，待选井水平井组中包括多个水平井。所述生产指标还包括各水平井的短期累产、中期累产和长期累产以及水平段长度，在本文实施例中，各水平井的进行重复压裂的周期为2至5年，所述短期累产为水平井上一次进行重复压裂后短期的累计产油量，作为优选地，本文实施例中所述的短期累产为水平井上一次进行重复压裂后30天的累计产油量。所述中期累产与所述长期累产的含义与所述短期累产相似，作为优选地，本文实施例中所述的中期累产为水平井上一次进行重复压裂后90天的累计产油量，所述长期累产为水平井上一次进行重复压裂后360天的累计产油量。所述水平段长度为水平井着陆点到完钻井深的长度，在水平段中开采石油，由于各水平井的开采地理位置以及储油层位置不同，因此各水平井的水平段长度也不同，分别获取各水平井的水平段长度用于重复压裂选井。所述阈值可以根据人工经验得到，例如，综合指标大于90的水平井有推荐重复压裂井，可以直接进行重复压裂，无需再进行重复压裂计算，综合指标在80-90之间的水平井为建议考虑重复压裂井，能够进行重复压裂计算进一步确定重复压裂的结果。

通过本文实施例的方法，以被窜扰影响产量为指标来表征页岩油水平井的被窜扰程度，并计算待选井水平井组中各水平井包括被窜扰程度反向百分位数的综合指标，根据综合指标和设定的阈值选出能够进行重复压裂计算的水平井，考虑了井间窜扰对选井的影响，提高了选井的精度。本文所述方法能够在待选井水平井组中选出储层品质较好、初次改造不足以及被窜扰影响不严重的水平井作为能够进行重复压裂计算的水平井，以便于后续对选出的水平井进行进一步的重复压裂计算，减小了重复压裂计算的工作量。

根据本文的一个实施例，为了进一步减小计算量，提高选井效率，步骤201获取待选井水平井组中各水平井产生窜扰前的产油水平以及窜扰期间的产油量，计算各水平井的被窜扰影响产量之前，还包括，剔除所述待选井水平井组中井筒损坏到无法进行重复压裂的水平井。

在本步骤中，若某个水平井由于地理位置或储油层分布等因素，进行上一次重复压裂后，该水平井的井筒已经损坏到无法再次进行重复压压裂，因此仍计算该水平井的被窜扰影响产量并计算其综合指标没有意义，且会增加计算量。因此通过人工经验等方式对待选井水平井组中各水平井的井筒进行观察分析，剔除井筒损坏到无法进行重复压裂的水平井。

根据本文的一个实施例，步骤201中计算各水平井的被窜扰影响产量的公式为(1)：

IP_Affected_Reverse＝IP×d-IP_Affected (1)

其中，IP_Affected_Reverse表示被窜扰影响产量，IP表示产生窜扰前的产油水平，d表示窜扰的总时间，IP_Affected表示窜扰期间的产油量。

在本步骤中，可以通过页岩油水平井的生产曲线得出该水平井的产生窜扰的时间、窜扰的总时间以及产生窜扰前的产油水平、窜扰期间的产油量等。

具体地，页岩油水平井的生产曲线可以如图7所示，其中包括日产液曲线、日产油曲线、含水量曲线等，当某水平井产生窜扰时，该井的日产油量会急速下降，日产液量和含水量急速上升，对日产液曲线、日产油曲线、含水量曲线分别进行积分，求得各曲线急速下降部分子曲线和急速上升部分子曲线的拐点，对日产液曲线、日产油曲线、含水量曲线急速下降部分子曲线的拐点取平均值作为产生窜扰的时间，对日产液曲线、日产油曲线、含水量曲线急速上升部分子曲线的拐点取平均值作为窜扰结束的时间，进而得到窜扰的总时间，并通过日产油曲线计算产生窜扰前的产油水平和窜扰期间的产油量。

根据本文一个实施例，步骤202所述的综合指标还包括储层品质百分位数以及初次改造百分位数。

在本步骤中，所述储层品质百分位数表示水平井的储油层品质，可以根据产油量等生产指标求得。所述初次改造效果表示水平井经过重复压裂后的产油效果，不限于第一次对水平井进行压裂，每进行一次重复压裂后都可以求得压裂后的产油效果作为所述初次改造效果。

根据本文的一个实施例，步骤202根据包括所述被窜扰影响产量的所述各水平井的生产指标，计算各水平井包括被窜扰程度反向百分位数的综合指标进一步包括，根据所述生产指标计算所述各水平井单位水平段的中期产能以及各水平井的初次改造效果；根据所述单位水平段的中期产能计算所述储层品质百分位数；根据所述初次改造效果计算所述初次改造效果百分位数；根据所述被窜扰影响产量计算所述被窜扰程度反向百分位数；根据所述各水平井的所述储层品质百分位数、所述初次改造百分位数、以及所述被窜扰程度反向百分位数计算所述各水平井的综合指标。

在本步骤中，计算所述水平井单位水平段的中期产能的公式为(2)：

Norm_IP90＝IP90/SLL (2)

其中，Norm_IP90表示所述单位水平段的中期产能，IP90表示所述中期累产，SLL表示所述水平段长度。作为优选地，本文实施例中所述中期累产为重复压裂后的90天累产。

计算所述初次改造效果的公式为(3)：

IP_first＝IP30/IP360 (3)

其中，IP_first表示所述初次改造效果，IP30表示所述短期累产，IP360表示所述长期累产。作为优选地，本文实施例中所述短期累产为重复压裂后的30天累产，所述长期累产为重复压裂后的360天累产。

根据本文的一个实施例，根据所述被窜扰影响产量计算所述被窜扰程度反向百分位数进一步包括，将所述各水平井的被窜扰影响产量进行升序排序得到被窜扰影响产量的排序表，根据所述水平井在所述被窜扰影响产量的排序表中的位置以及所述待选井水平井组中水平井的数量计算所述水平井的被窜扰程度反向百分位数。

在本步骤中，对各水平井的被窜扰影响产量进行升序排序，再根据被窜扰影响产量的排序表和待选井水平井组中水平井的数量计算每个水平井的被窜扰程度反向百分位数，得到了各水平井的被窜扰影响产量在待选井水平井组中的精确排名。与现有技术中的选井方法相比，在选井过程中加入被窜扰影响产量，考虑了井间窜扰对选井的影响，提高了选井精度。

此外，还对各水平井的单位水平段的中期产能进行升序排序得到单位水平段的中期产能的排序表，再根据单位水平段的中期产能的排序表和待选井水平井组中水平井的数量计算每个水平井的储层品质百分位数，得到了各水平井的储层品质在待选井水平井组中的精确排名。因此，在选井过程中，除了考虑井间窜扰的影响之外，还考虑了水平井的储油层品质，进一步提高选井的精度。

此外，还对各水平井的初次改造效果进行升序排序得到初次改造效果的排序表，再根据初次改造效果的排序表和待选井水平井组中水平井的数量计算每个水平井的初次改造效果百分位数，得到了各水平井的初次改造效果在待选井水平井组中的精确排名。因此，在选井过程中，除了考虑井间窜扰的影响和储油层品质之外，还考虑了水平井的初次改造效果，进一步提高选井的精度。

因此，在选井过程中加入被窜扰影响产量，考虑了井间窜扰对选井的影响，此外，还加入水平井的储油层品质和初次改造效果，综合选出被窜扰影响不严重、储油层品质较好、初次改造效果不足的水平井作为能够进行重复压裂计算的水平井，提高了选井的精度，减小了对待选井水平井组进行重复压裂计算的工作量。

根据本文的一个实施例，计算所述被窜扰程度反向百分位数的公式为(4)：

其中，k表示所述水平井在所述被窜扰影响产量的排序表中的位置，IP_Affected_Reverse_Percentile(z(k))表示该水平井的被窜扰程度反向百分位数，m表示所述待选井水平井组中水平井的数量。

计算所述储层品质百分位数的公式为(5)：

其中，i表示所述水平井在所述单位水平段的中期产能的排序表中的位置，Norm_IP90_Percentile(x(i))表示该水平井的储层品质百分位数，m表示所述待选井水平井组中水平井的数量。

计算所述初次改造效果百分位数的公式为(6)：

其中，j表示所述水平井在所述初次改造效果的排序表中的位置，IP_first_Percentile(y(j))表示该水平井的初次改造效果百分位数，m表示所述待选井水平井组中水平井的数量。

根据本文的一个实施例，根据所述各水平井的所述储层品质百分位数、所述初次改造百分位数、以及所述被窜扰程度反向百分位数计算所述各水平井的综合指标进一步包括，对所述储层品质百分位数、所述初次改造百分位数、以及所述被窜扰程度反向百分位数进行加权并取平均值，得到所述综合指标。

在本步骤中，所述权值可以根据该水平井的生产指标计算得到，例如，获取各水平井的历史生产数据，分别计算各水平井的储层品质百分位数、初次改造百分位数、以及被窜扰程度反向百分位数，通过主成分分析法确定权重，然后计算各水平井的综合指标。也可以取各水平井的储层品质百分位数、初次改造百分位数、以及被窜扰程度反向百分位数的均值作为自身的综合指标，其公式为(7)：

其中，C表示一个所述水平井的综合指标，NormIP90_Percentile表示该水平井的所述储层品质百分位数，IP_first_Percentile表示该水平井的所述初次改造百分位数，IP_Affected_Reverse_Percentile表示该水平井的所述被窜扰程度反向百分位数。

如图3所示为本文实施例一种页岩油水平井重复压裂选井装置的详细结构图，在本图中具体描述了页岩油水平井重复压裂选井装置的详细结构，包括被窜扰影响产量计算单元301、水平井综合指标计算单元302、选井单元303。

根据本文的一个实施例，被窜扰影响产量计算单元301进一步包括初筛模块3011，剔除待选井水平井组中井筒损坏到无法进行重复压裂的水平井，进一步减小计算量，提高选井效率。可以通过人工经验等方式对待选井水平井组中各水平井的井筒进行观察分析，剔除井筒损坏到无法进行重复压裂的水平井。

根据本文的一个实施例，被窜扰影响产量计算单元301进一步包括被窜扰时间计算模块3012，根据页岩油水平井的生产曲线计算产生窜扰的时间以及窜扰的总时间。具体地，页岩油水平井的生产曲线可以如图7所示，当某水平井产生窜扰时，该井的日产油量会急速下降，日产液量和含水量急速上升，可以对日产液曲线、日产油曲线、含水量曲线分别进行积分，求得日产液曲线、日产油曲线、含水量曲线急速下降部分子曲线和急速上升部分子曲线的拐点，对日产液曲线、日产油曲线、含水量曲线急速下降部分子曲线的拐点取平均值作为产生窜扰的时间，对各曲线急速上升部分子曲线的拐点取平均值作为窜扰结束的时间，进而得到窜扰的总时间，以便于被窜扰影响产量计算模块3014根据窜扰的总时间计算被窜扰影响产量。

根据本文的一个实施例，被窜扰影响产量计算单元301进一步包括产量获取模块3013，通过水平井的生产曲线获取水平井的30天累产、90天累产、360天累产以及产生窜扰前的平均天产、窜扰期间的产油量等，以便于被窜扰影响产量计算模块3014根据产生窜扰前的平均天产、窜扰期间的产油量计算被窜扰影响产量，以及水平井综合指标计算单元302计算水平井的综合指标。

根据本文的一个实施例，被窜扰影响产量计算单元301进一步包括被窜扰影响产量计算模块3014，根据被窜扰时间计算模块3012得到的被窜扰总时间以及产量获取模块3013获取到的产生窜扰前的平均天产、窜扰期间的产油量利用本说明书中公式(1)计算被窜扰影响产量。

根据本文的一个实施例，水平井综合指标计算单元302进一步包括单位水平段的中期产能计算模块3021，根据产量获取模块3013获取到水平井的90天累产以及该水平井的水平段长度，利用本说明书公式(2)计算水平井单位水平段的中期产能。

根据本文的一个实施例，水平井综合指标计算单元302进一步包括初次改造效果计算模块3022，根据产量获取模块3013获取到水平井的30天累产以及360天累产利用本说明书公式(3)计算水平井的初次改造效果。

根据本文的一个实施例，水平井综合指标计算单元302进一步包括储层品质百分位数计算模块3023，首先对单位水平段的中期产能计算模块3021计算的各水平井的中期产能进行升序排序，得到单位水平段的中期产能的排序表，然后根据单位水平段的中期产能的排序表以及待选井水平井组中水平井的数量利用本说明书公式(5)计算各水平井的储层品质百分位数。

根据本文的一个实施例，水平井综合指标计算单元302进一步包括初次改造效果百分位数计算模块3024，首先对初次改造效果计算模块3022计算的各水平井的初次改造效果进行升序排序，得到初次改造效果的排序表，然后根据初次改造效果的排序表以及待选井水平井组中水平井的数量利用本说明书公式(6)计算各水平井的初次改造效果百分位数。

根据本文的一个实施例，水平井综合指标计算单元302进一步包括被窜扰程度反向百分位数计算模块3025，首先对被窜扰影响产量计算单元301得到的各水平井的被窜扰影响产量进行升序排序，得到被窜扰影响的产量的排序表，然后根据被窜扰影响的产量的排序表以及待选井水平井组中水平井的数量利用本说明书公式(4)计算各水平井的被窜扰程度反向百分位数。

根据本文的一个实施例，水平井综合指标计算单元302进一步包括综合指标计算模块3026，对储层品质百分位数计算模块3023得到的储层品质百分位数、初次改造效果百分位数计算模块3024得到的初次改造效果百分位数以及被窜扰程度反向百分位数计算模块3025得到的被窜扰程度反向百分位数进行加权并取平均值，得到所述综合指标。可以通过产量获取模块3013获取各水平井的历史生产数据，然后储层品质百分位数计算模块3023、初次改造效果百分位数计算模块3024和被窜扰程度反向百分位数计算模块3025分别计算各水平井的储层品质百分位数、初次改造百分位数、以及被窜扰程度反向百分位数，通过主成分分析法确定权重，然后计算各水平井的综合指标。也可以取各水平井的储层品质百分位数、初次改造百分位数、以及被窜扰程度反向百分位数利用本说明书公式(7)求得均值作为综合指标。

根据本文的一个实施例，选井单元303进一步包括阈值设定模块3031，根据人工经验得到阈值，例如，综合指标大于90的水平井有推荐重复压裂井，可以直接进行重复压裂，无需再进行重复压裂计算，综合指标在80-90之间的水平井为建议考虑重复压裂井，能够进行重复压裂计算进一步确定重复压裂的结果。

根据本文的一个实施例，选井单元303进一步包括综合指标对比模块3032，对比综合指标计算模块3026得到的待选井水平井组中各水平井的综合指标以及阈值设定模块3031设定的阈值，选出能够进行重复压裂计算的水平井。

如图4所示为本文实施例利用被窜扰影响产量进行重复压裂选井的流程图，本文实施例所述的待选井水平井组中包括多个水平井，本图中描述了利用被窜扰影响产量进行重复压裂选井的流程。具体包括：

步骤401：分析待选井水平井组中各水平井的生产曲线。

在本步骤中，所述水平井的生产曲线如图7所示，其中包括日产液曲线、日产油曲线、含水量曲线，可以通过生产曲线获得该水平井的短期累产、长期累产、中期累产以及计算被窜扰影响产量。

步骤402：获取短期累产(IP30)。

在本步骤中，根据步骤401的各水平井的生产曲线分别获取各水平井的短期累产，即经过重复压裂后的30天累产。

步骤403：获取长期累产(IP360)。

在本步骤中，根据步骤401的各水平井的生产曲线分别获取各水平井的长期累产，即经过重复压裂后的360天累产。

步骤404：获取中期累产(IP90)。

在本步骤中，根据步骤401的各水平井的生产曲线分别获取各水平井的中期累产，即经过重复压裂后的90天累产。

步骤405：计算被窜扰影响产量(IP_Affected)。

在本步骤中，根据步骤401的各水平井的生产曲线分别计算各水平井的被窜扰影响产量(IP_Affected)。

步骤406：计算初次改造效果(IP_First)。

在本步骤中，根据步骤402获取的水平井的30天累产(IP30)以及步骤403获取的水平井的360天累产(IP360)根据本说明书公式(3)计算水平井的初次改造效果(IP_First)。

步骤407：计算单位水平段的中期产能(Norm_IP90)。

在本步骤中，根据步骤404获取的水平井的90天累产(IP90)以及该水平井的水平段长度利用本说明书公式(2)计算单位水平段的中期产能(Norm_IP90)。

步骤408：计算初次改造效果百分位数(IP_First_Percentile)。

在本步骤中，将步骤406计算的各水平井的初次改造效果(IP_First)进行升序排序得到初次改造效果的排序表，再根据初次改造效果的排序表和待选井水平井组中水平井的数量利用本说明书公式(6)计算每个水平井的初次改造效果百分位数(IP_First_Percentile)。

步骤409：计算储层品质百分位数(Norm_IP90_Percentile)。

在本步骤中，将步骤407计算的各水平井的单位水平段的中期产能(Norm_IP90)进行升序排序得到单位水平段的中期产能的排序表，再根据单位水平段的中期产能的排序表和待选井水平井组中水平井的数量利用本说明书公式(5)计算每个水平井的储层品质百分位数(Norm_IP90_Percentile)。

步骤410：计算被窜扰程度反向百分位数(IP_Affected_Reverse_Percentile)。

在本步骤中，将步骤405计算的各水平井的被窜扰影响产量(IP_Affected)进行升序排序，再根据被窜扰影响产量的排序表和待选井水平井组中水平井的数量利用本说明书公式(4)计算每个水平井的被窜扰程度反向百分位数(IP_Affected_Reverse_Percentile)。

步骤411：计算综合指标。

在本步骤中，根据步骤408得到的初次改造效果百分位数(IP_First_Percentile)、步骤409得到的储层品质百分位数(Norm_IP90_Percentile)以及步骤410得到的被窜扰程度反向百分位数(IP_Affected_Reverse_Percentile)利用本说明书公式(7)计算各水平井的综合指标。

步骤412：根据综合指标和阈值进行重复压裂选井。

在本步骤中，分别比较各水平井的综合指标和设定的阈值，选出能够进行重复压裂计算的水平井。

通过图4描述的方法，在选井过程中加入被窜扰影响产量，考虑了井间窜扰对选井的影响，此外，还加入水平井的储油层品质和初次改造效果，综合选出被窜扰影响产量最小、储油层品质最好、初次改造效果最好的水平井进行重复压裂计算，提高了选井的精度，以便于进一步对选井结果中的水平井进行重复压裂计算，得到重复压裂结果，减小了重复压裂计算的工作量。

进一步地，步骤405计算待选井水平井组中一个水平井的被窜扰影响产量(IP_Affected)的过程如图5所示，具体包括：

步骤501：分析页岩油水平井的生产曲线，得到该水平井产生窜扰的时间、窜扰的总时间以及产生窜扰前的产油水平、窜扰期间的产油量。

在本步骤中，页岩油水平井的生产曲线如图7所示，当某水平井产生窜扰时，该井的日产油量会急速下降，日产液量和含水量急速上升，可以对日产液曲线、日产油曲线、含水量曲线分别进行积分，求得日产液曲线、日产油曲线、含水量曲线急速下降部分子曲线和急速上升部分子曲线的拐点，对日产液曲线、日产油曲线、含水量曲线急速下降部分子曲线的拐点取平均值作为产生窜扰的时间，对日产液曲线、日产油曲线、含水量曲线急速上升部分子曲线的拐点取平均值作为窜扰结束的时间，进而得到窜扰的总时间，并通过日产油曲线计算产生窜扰前的产油水平和窜扰期间的产油量。

步骤502：计算该水平井被窜扰影响产量。

在本步骤中，根据步骤501得到的窜扰的总时间、产生窜扰前的产油水平以及窜扰期间的产油量利用本说明书公式(1)计算该水平井的被窜扰影响产量。

进一步地，步骤410计算待选井水平井组中一个水平井的被窜扰程度反向百分位数(IP_Affected_Reverse_Percentile)的过程如图6所示，具体包括：

步骤601：对被窜扰影响产量进行升序排序，得到被窜扰影响产量的排序表。

在本步骤中，对图4所述方法得到的待选井水平井组中各水平井的被窜扰影响产量进行升序排序。

步骤602：根据该水平井在被窜扰影响产量的排序表中的位置以及待选井水平井组中水平井的数量计算该水平井的被窜扰程度反向百分位数。

在本步骤中，根据水平井在步骤601得到的被窜扰影响产量的排序表中的位置以及待选井水平井组中水平井的数量利用本说明书公式(4)计算该水平井的被窜扰程度反向百分位数。

如图8所示为本文实施例计算机设备的结构示意图，本文中的页岩油水平井重复压裂选井装置可以为本实施例中的计算机设备，执行上述本文的方法。计算机设备802可以包括一个或多个处理设备804，诸如一个或多个中央处理单元(CPU)，每个处理单元可以实现一个或多个硬件线程。计算机设备802还可以包括任何存储资源806，其用于存储诸如代码、设置、数据等之类的任何种类的信息。非限制性的，比如，存储资源806可以包括以下任一项或多种组合：任何类型的RAM，任何类型的ROM，闪存设备，硬盘，光盘等。更一般地，任何存储资源都可以使用任何技术来存储信息。进一步地，任何存储资源可以提供信息的易失性或非易失性保留。进一步地，任何存储资源可以表示计算机设备802的固定或可移除部件。在一种情况下，当处理设备804执行被存储在任何存储资源或存储资源的组合中的相关联的指令时，计算机设备802可以执行相关联指令的任一操作。计算机设备802还包括用于与任何存储资源交互的一个或多个驱动机构808，诸如硬盘驱动机构、光盘驱动机构等。

计算机设备802还可以包括输入/输出模块810(I/O)，其用于接收各种输入(经由输入设备812)和用于提供各种输出(经由输出设备814)。一个具体输出机构可以包括呈现设备816和相关联的图形用户接口(GUI)818。在其他实施例中，还可以不包括输入/输出模块810(I/O)、输入设备812以及输出设备814，仅作为网络中的一台计算机设备。计算机设备802还可以包括一个或多个网络接口820，其用于经由一个或多个通信链路822与其他设备交换数据。一个或多个通信总线824将上文所描述的部件耦合在一起。

通信链路822可以以任何方式实现，例如，通过局域网、广域网(例如，因特网)、点对点连接等、或其任何组合。通信链路822可以包括由任何协议或协议组合支配的硬连线链路、无线链路、路由器、网关功能、名称服务器等的任何组合。

本文实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如下步骤：

获取待选井水平井组中各水平井产生窜扰前的产油水平以及窜扰期间的产油量，计算各水平井的被窜扰影响产量；

根据包括所述被窜扰影响产量的所述各水平井的生产指标，计算各水平井包括被窜扰程度反向百分位数的综合指标；

根据所述各水平井的综合指标以及设定的阈值选出能够进行重复压裂计算的水平井。

对应于图2、图4-图6中的方法，本文实施例还提供了一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该计算机程序被处理器运行时执行上述方法的步骤。

本文实施例还提供一种计算机可读指令，其中当处理器执行所述指令时，其中的程序使得处理器执行如图2、图4-图6所示的方法。

应理解，在本文的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本文实施例的实施过程构成任何限定。

还应理解，在本文实施例中，术语“和/或”仅仅是一种描述关联对象的关联关系，表示可以存在三种关系。例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。

本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本文的范围。

所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

在本文所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口、装置或单元的间接耦合或通信连接，也可以是电的，机械的或其它的形式连接。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本文实施例方案的目的。

另外，在本文各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。

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

本文中应用了具体实施例对本文的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本文的方法及其核心思想；同时，对于本领域的一般技术人员，依据本文的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本文的限制。