一种缝洞型油藏注水开发效果评价方法.pdf

本发明提供了一种缝洞型油藏注水开发效果评价方法，该方法包括：将待评价缝洞单元划分成多个注采井组；计算各注采井组的井控系数，结合预定井控系数指标确定各井组的井控效果区间；计算各注采井组的水驱采收率及吨水换油率，结合预定开发效果指标确定各井组的开发效果区间；计算各注采井组的累积含水率及存水率，根据累积含水率判断是否存在水窜，结合预定存水率指标确定存水效果区间；计算各注采井组的无因次弹性产量比、单储压降及人工水驱指数，结合预定天然能量指标确定各井组的天然能量区间，结合预定人工水驱指标确定水驱效果区间。本发明全面反映了注采井组的注水效果以及达到该效果的原因，使提出的调整措施更有基础依据和针对性。

1.  一种缝洞型油藏注水开发效果评价方法，其特征在于，所述的评价方法包括：
根据地质资料、示踪剂检测资料以及生产动态响应，将待评价缝洞单元划分成多个注采井组，其中各所述注采井组分别包括多个生产井；
将每个注采井组的生产井划分为钻遇溶洞井及未钻遇溶洞井，计算钻遇溶洞井的井控系数A_D及未钻遇溶洞井的井控系数B_D，根据预定井控系数指标、所述钻遇溶洞井的井控系数A_D及未钻遇溶洞井的井控系数B_D确定所述注采井组的井控效果区间；其中，所述钻遇溶洞井的井控系数为：
AD=Σi=1nv1ωimax(VeiVti)]]>
其中，ωi=VtiΣi=1nv1Vti×100%;]]>
以及，所述未钻遇溶洞井的井控系数为：
SD=Σm=1nwΣj=1nv2ωmjfmjΘ(x)nw]]>
其中，ωmj=VtmjΣj=1nv2Vtmj,]]>Θ(x)=1dvmj&GreaterEqual;dwm0dvmj<dwm;]]>
计算各所述注采井组的水驱采收率及吨水换油率，根据预定开发效果指标及所述水驱采收率及吨水换油率确定所述注采井组的开发效果区间；
根据生产动态资料，分别计算各所述注采井组的累积含水率及存水率，根据所述累积含水率判断是否存在水窜，并根据预定存水率指标及所述存水率确定存水效果区间；
根据所述生产动态资料，分别计算各所述注采井组的无因次弹性产量比、单储压降以及人工水驱指数，根据预定天然能量指标、所述无因次弹性产量比及单储压降确定所述注采井组的天然能量区间，并根据预定人工水驱指标及所述人工水驱指数确定水驱效果区间；
其中，ω_i为所述注采井组中钻遇的第i个溶洞体积占所述注采井组中所有钻遇溶洞总体积的百分比；V_ei为生产井与钻遇的第i溶洞接触点下方溶洞体积；V_ti为钻遇 的第i个溶洞的总体积；n_v1为钻遇溶洞的数目；ω_mj为所述注采井组中第m口生产井所对应的未钻遇的第j个溶洞体积占所述第m口生产井控制的所有未钻遇溶洞总体积的百分比；f_mj为第m口生产井与其所对应的未钻遇的第j个溶洞之间的连通系数，反映生产井与溶洞空间连通程度的大小，与两者的距离和裂缝渗透率有关；Θ(x)为海维塞德函数，决定于溶洞与对应生产井的相对高低位置关系；n_w为所述注采井组中未钻遇溶洞生产井的数目；n_v2为所述注采井组中未被钻遇溶洞的数目；V_tmj为第m口生产井所对应的未钻遇的第j个溶洞的体积；d_vmj为第m口生产井所对应的未钻遇的第j个溶洞的深度；d_wm为第m口生产井的射孔深度。

2.  根据权利要求1所述的评价方法，其特征在于，根据地质资料、示踪剂检测资料以及生产动态响应，将待评价缝洞单元划分成多个注采井组，其中各所述注采井组分别包括多个生产井，包括：
根据所述待评价缝洞单元的地质资料，确定所述待评价缝洞单元内的缝洞网络，按照所述缝洞网络中地下河、廊道及断层的走向将所述待评价缝洞单元初步划分成多个注采井组；
从每个所述注采井组中选取至少一个生产井，向所选取生产井中投放示踪剂，检测所选取生产井的周边生产井中是否存在示踪剂以确定所选取生产井与其周边生产井的连通关系，根据所述连通关系对初步划分的注采井组进行二次划分；
对于不存在示踪剂的周边生产井，根据包含含水率及压力的生产动态响应，对不存在示踪剂的周边生产井进行注采井组划分。

3.  根据权利要求1所述的评价方法，其特征在于，根据预定井控系数指标、所述钻遇溶洞井的井控系数A_D及未钻遇溶洞井的井控系数B_D确定所述注采井组的井控效果区间，包括：
根据所述钻遇溶洞井的井控系数A_D及未钻遇溶洞井的井控系数B_D，采用体积加权法确定各所述注采井组的总井控系数λ_w：
λw=Σi=1nv1VtiNtAD+Σm=1nwΣj=1nv2VtmjNtSD;]]>
根据预定井控系数指标及各所述注采井组的总井控系数λ_w确定所述注采井组的井控效果区间；
其中，N_t为所述注采井组内所有储集空间的体积。

4.  根据权利要求1所述的评价方法，其特征在于，计算各所述注采井组的水驱采收率及吨水换油率，包括：
计算各所述注采井组的水驱采收率ER：
ER=Np-NpiN×100%;]]>
计算各所述注采井组的吨水换油率N_p：
Ω=NpBoWpBw×100%;]]>
其中，N_p为所述注采井组的累产油；N_pi为所述注采井组注水前的累产油；N为所述注采井组的地质储量；B_o为原油的体积系数；W_p为所述注采井组的累产水；B_w为水的体积系数。

5.  根据权利要求1所述的评价方法，其特征在于，计算各所述注采井组的累积含水率及存水率，包括：
计算各所述注采井组的累积含水率f_wc：
fwc=WpBwNpBo+WpBw;]]>
计算各所述注采井组的存水率η：
η=Wi-Wp-WeWi,]]>
其中，W_e＝N_pB_o+W_pB_w-NB_oiC_tΔp-W_iB_w；
其中，W_i为所述注采井组的累积注水量；W_e为所述注采井组的漏失量；W_e为漏失量；N为地质储量；B_oi为原始地层压力下原油的体积系数；C_t为综合压缩系数；Δp为压力下降量。

6.  根据权利要求5所述的评价方法，其特征在于，根据所述累积含水率判断是否存在水窜，包括：
根据所述累积含水率计算单位时间内各所述注采井组的累积含水率变化Δf_wc；
求取所述待评价缝洞单元内所有注采井组的累积含水率变化的均值；
将所述各注采井组累积含水率变化Δf_wc与所述累积含水率变化的均值进行比较，判断各所述注采井组是否存在水窜。

7.  根据权利要求1所述的评价方法，其特征在于，所述计算各所述注采井组的无因次弹性产量比、单储压降及人工水驱指数，包括：
计算各所述注采井组的无因次弹性产量比D_pr：
Dpr=ΔpN100Np=Bo100BoiCt-We-WpBw100BoiCtNp;]]>
计算各所述注采井组的单储压降N_pr：
Npr=NpBo+WpBwNBoiCtΔp=1+WeNBoiCtΔp;]]>
计算各所述注采井组的人工水驱指数ψ：
ψ=WiBw-WpBw-W1BwNpBo;]]>
其中，Δp为所述注采井组的压力下降量。

一种缝洞型油藏注水开发效果评价方法
技术领域
本发明涉及缝洞型油藏注水开发技术领域，尤其涉及一种缝洞型油藏注水开发效果评价方法。
背景技术
目前有关油藏注水开发技术存在诸多弊端，例如，塔河缝洞型碳酸盐岩油藏大部分区块实施了注水替油和缝洞单元注水开发，取得了较好的效果，但目前采收率只有13％-15％，部分注采单元含水上升严重，稳油控水难度逐渐加大，需要继续深化开展注水工作及采取适当的调整措施，实现老井增产增效，提高油藏采收率。因此，需要深入评价注水效果，分析注水过程中存在的问题及其原因，从而制定合理的注水开发技术对策。
目前有关油藏注水开发效果评价主要集中于常规的砂岩油藏，现有技术中，汪子昊针对典型砂岩油藏采用行业标准对比和同类油田开发指标对比的方法，结合单因素指标的评价研究了注水开发以后的储量变化趋势及储量动用程度、产量变化趋势及产量变化规律、地层压力保持水平、油田综合含水变化规律、油田注入水利用效率等。李兴训从传统的油藏工程方法到现代模糊数学理论及数值模拟等方面阐述了如何评价砂岩油藏水驱开发效果。罗银富结合油藏工程、系统工程、模糊数学、层次分析原理等方法，引入启动压力梯度的概念，研究了低渗透砂岩油藏的水驱效果评价方法。在缝洞型油藏注水开发技术研究方面，公开号为CN101942984A的中国专利申请提出了一种碳酸盐岩缝洞型油藏注水替油开采方法，通过试注、补充能量、关井、采油等步骤实现缝洞型油藏高效开发。马旭杰通过数值模拟和注水实践总结，提出了“低注高采、试注判断连通实时调整、单元整体温和注水、井组差异管理”的缝洞型油藏注水开发模式。谷雨在分析注水替油基本原理的基础上，从岩溶地貌特征、储集体类型、完井位置、水体能量及含水特征、原油物性、溶洞定容性等六个方面总结出注水替油的选井原则。
但是，由于缝洞型油藏不同于常规的砂岩油藏，储集体空间离散分布，不同部位 的缝洞尺寸、形状等存在较大的差异，传统的砂岩油藏注水开发效果方法对缝洞型油藏不适用。而有关缝洞型油藏注水开发研究又主要集中于基于经验的注水方式研究，缺乏注水效果评价方面的研究，并且注水效果评价的研究且存在较大困难，具体表现为以下几个方面：(1)缝洞单元内部储集体分布、注采关系复杂，难以有效确定或者划分注采井组，并将注采井组作为基本研究对象；(2)不同注采井组储集体内部连通关系复杂、缝洞井网匹配关系各样、天然能量来源不同；(3)不同注采井组生产动态个性化强、可比性差；(4)现有的注水开发决策制定仅依靠过往经验，缺乏有效地理论指导；(5)目前没有针对缝洞型油藏注水开发效果评价的技术方法体系，导致无法开展注水调整，使后续调整措施的制定缺乏指导。
发明内容
为了解决上述技术问题，本发明提供了一种缝洞型油藏注水开发效果评价方法，所述的评价方法包括：
根据地质资料、示踪剂检测资料以及生产动态响应，将待评价缝洞单元划分成多个注采井组，其中各所述注采井组分别包括多个生产井；
将每个注采井组的生产井划分为钻遇溶洞井及未钻遇溶洞井，计算钻遇溶洞井的井控系数A_D及未钻遇溶洞井的井控系数B_D，根据预定井控系数指标、所述钻遇溶洞井的井控系数A_D及未钻遇溶洞井的井控系数B_D确定所述注采井组的井控效果区间；其中，所述钻遇溶洞井的井控系数为：
AD=Σi=1nv1ωimax(VeiVti)]]>
其中，ωi=VtiΣi=1nv1Vti×100%;]]>
以及，所述未钻遇溶洞井的井控系数为：
BD=Σm=1nwΣj=1nv2ωmjfmjΘ(x)nw]]>
其中，ωmj=VtmjΣj=1nv2Vtmj,Θ(x)=1dvmj&GreaterEqual;dwm0dvmj<dwm;]]>
计算各所述注采井组的水驱采收率及吨水换油率，根据预定开发效果指标及所述水驱采收率及吨水换油率确定所述注采井组的开发效果区间；
根据生产动态资料，分别计算各所述注采井组的累积含水率及存水率，根据所述累积含水率判断是否存在水窜，并根据预定存水率指标及所述存水率确定存水效果区间；
根据所述生产动态资料，分别计算各所述注采井组的无因次弹性产量比、单储压降以及人工水驱指数，根据预定天然能量指标、所述无因次弹性产量比及单储压降确定所述注采井组的天然能量区间，并根据预定人工水驱指标及所述人工水驱指数确定水驱效果区间；
其中，ω_i为所述注采井组中钻遇的第i个溶洞体积占所述注采井组中所有钻遇溶洞总体积的百分比；V_ei为生产井与钻遇的第i溶洞接触点下方溶洞体积；V_ti为钻遇的第i个溶洞的总体积；n_v1为钻遇溶洞的数目；ω_mj为所述注采井组中第m口生产井所对应的未钻遇的第j个溶洞体积占所述第m口生产井控制的所有未钻遇溶洞总体积的百分比；f_mj为第m口生产井与其所对应的未钻遇的第j个溶洞之间的连通系数，反映生产井与溶洞空间连通程度的大小，与两者的距离和裂缝渗透率有关；Θ(x)为海维塞德函数，决定于溶洞与对应生产井的相对高低位置关系；n_w为所述注采井组中未钻遇溶洞生产井的数目；n_v2为所述注采井组中未被钻遇溶洞的数目；V_tmj为第m口生产井所对应的未钻遇的第j个溶洞的体积；d_vmj为第m口生产井所对应的未钻遇的第j个溶洞的深度；d_wm为第m口生产井的射孔深度。
在一实施例中，根据地质资料、示踪剂检测资料以及生产动态响应，将待评价缝洞单元划分成多个注采井组，其中各所述注采井组分别包括多个生产井，包括：
根据所述待评价缝洞单元的地质资料，确定所述待评价缝洞单元内的缝洞网络，按照所述缝洞网络中地下河、廊道及断层的走向将所述待评价缝洞单元初步划分成多个注采井组；
从每个所述注采井组中选取至少一个生产井，向所选取生产井中投放示踪剂，检测所选取生产井的周边生产井中是否存在示踪剂以确定所选取生产井与其周边生产井的连通关系，根据所述连通关系对初步划分的注采井组进行二次划分；
对于不存在示踪剂的周边生产井，根据包含含水率及压力的生产动态响应，对不存在示踪剂的周边生产井进行注采井组划分。
在一实施例中，根据预定井控系数指标、所述钻遇溶洞井的井控系数A_D及未钻遇溶洞井的井控系数B_D确定所述注采井组的井控效果区间，包括：
根据所述钻遇溶洞井的井控系数A_D及未钻遇溶洞井的井控系数B_D，采用体积加权法确定各所述注采井组的总井控系数λ_w：
λw=Σi=1nv1VtiNtAD+Σm=1nwΣj=1nv2VtmjNtBD;]]>
根据预定井控系数指标及各所述注采井组的总井控系数λ_w确定所述注采井组的井控效果区间；
其中，N_t为所述注采井组内所有储集空间的体积。
在一实施例中，计算各所述注采井组的水驱采收率及吨水换油率，包括：
计算各所述注采井组的水驱采收率ER：
ER=Np-NpiN×100%;]]>
计算各所述注采井组的吨水换油率N_p：
Ω=NpBoWpBw×100%;]]>
其中，N_p为所述注采井组的累产油；N_pi为所述注采井组注水前的累产油；N为所述注采井组的地质储量；B_o为原油的体积系数；W_p为所述注采井组的累产水；B_w为水的体积系数。
在一实施例中，计算各所述注采井组的累积含水率及存水率，包括：
计算各所述注采井组的累积含水率f_wc：
fwc=WpBwNpBo+WpBw;]]>
计算各所述注采井组的存水率η：
η=Wi-Wp-WeWi,]]>
其中，W_e＝N_pB_o+W_pB_w-NB_oiC_tΔp-W_iB_w；
其中，W_i为所述注采井组的累积注水量；W_e为所述注采井组的漏失量；W_e为漏失量；N为地质储量；B_oi为原始地层压力下原油的体积系数；C_t为综合压缩系数；Δp为压力下降量。
在一实施例中，根据所述累积含水率判断是否存在水窜，包括：
根据所述累积含水率计算单位时间内各所述注采井组的累积含水率变化Δf_wc；
求取所述待评价缝洞单元内所有注采井组的累积含水率变化的均值；
将所述各注采井组累积含水率变化Δf_wc与所述累积含水率变化的均值进行比较，判断各所述注采井组是否存在水窜。
在一实施例中，所述计算各所述注采井组的无因次弹性产量比、单储压降及人工水驱指数，包括：
计算各所述注采井组的无因次弹性产量比D_pr：
Dpr=ΔpN100Np=Bo100BoiCt-We-WpBw100BoiCtNp;]]>
计算各所述注采井组的单储压降N_pr：
Npr=NpBo+WpBwNBoiCtΔp=1+WeNBoiCtΔp;]]>
计算各所述注采井组的人工水驱指数ψ：
ψ=WiBw-WpBw-W1BwNpBo;]]>
其中，Δp为所述注采井组的压力下降量。
本发明注采井组的划分可信度大大增加，引入井控系数评价现有井网与缝洞网络的配置关系这一不可控因素的好坏，并利用井控系数对注采井组进行分类，使具有可比性的注采井组进行后期注水效果对比。水窜指标和能量指标作为评价注入水利用率和注采井组能量特征的指标，可为评价注采井组注水效果提供依据；效果指标则实现了对注水效果好坏的定量评价；这些评价指标组成一个非常全面、完整的缝洞油藏注水开发评价体系；井控指标、水窜指标、能量指标和效果指标组成有机的整体，全面了反映注采井组的自身特征及其注水开发效果，使提出的调整措施更有基础依据和针对性；使用该评价方法评价的缝洞型油藏注采井组注水开发效果与现场实际开发效果吻合度大大增加。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例缝洞型油藏注水开发效果评价方法的流程图；
图2为本发明实施例将待评价缝洞单元划分为多个注采井组的流程图；
图3为本发明实施例基于储层地质结构初次划分的注采井组分布示意图；
图4为本发明实施例根据示踪剂监测结果确定井间连通性后调整的注采井组分布示意图；
图5为本发明实施例中生产井W1、W2及W3的含水率变化示意图；
图6为本发明实施例中向生产井W4注水后W5中的油压变化示意图；
图7为本发明实施例中向生产井W4停止注水后W6的产液量变化示意图；
图8为本发明实施例中向生产井W20注水后W19的产液量变化示意图；
图9为本发明实施例中W21开井生产后W19的产液量变化示意图；
图10为本发明实施例根据储层地质结构、示踪剂检测结果及生产动态响应综合因素划分的注采井组分布示意图；
图11为本发明实施例生产井钻遇溶洞情况下井控系数计算示意图；
图12为本发明实施例生产井未钻遇溶洞情况下井控系数计算示意图；
图13(A)～(F)为本发明实施例六个注采井组的累积含水率变化曲线。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明提供了一种缝洞型油藏注水开发效果评价方法，该方法具体包括以下步骤，如图1所示。
S1、根据地质资料、示踪剂检测资料以及生产动态响应，将待评价缝洞单元划分成多个注采井组，其中各注采井组中分别包括多个生产井。
S2、将每个注采井组的生产井划分为钻遇溶洞井及未钻遇溶洞井，计算钻遇溶洞井的井控系数A_D及未钻遇溶洞井的井控系数B_D，根据预定井控系数指标、A_D及B_D确定上述各注采井组的井控效果区间；其中，上述钻遇溶洞井的井控系数为：
AD=Σi=1nv1ωimax(VeiVti)---(1)]]>
其中，ωi=VtiΣi=1nv1Vti×100%---(2)]]>
以及，上述未钻遇溶洞井的井控系数为：
BD=Σm=1nwΣj=1nv2ωmjfmjΘ(x)nw---(3)]]>
其中，ωmj=VtmjΣj=1nv2Vtmj---(4)]]>
Θ(x)=1dvmj&GreaterEqual;dwm0dvmj<dwm---(5)]]>
S3、计算各注采井组的水驱采收率及吨水换油率，根据预定开发效果指标及水驱采收率及吨水换油率确定注采井组的开发效果区间；
S4、根据生产动态资料，分别计算各注采井组的累积含水率及存水率，根据上述累积含水率判断是否存在水窜，并根据预定存水率指标及上述存水率确定存水效果区间；
S5、根据生产动态资料，分别计算各注采井组的无因次弹性产量比、单储压降以及人工水驱指数，根据预定天然能量指标、上述无因次弹性产量比及单储压降确定各注采井组的天然能量区间，并根据预定人工水驱指标及上述人工水驱指数确定水驱效果区间；
其中，ω_i为每个注采井组中钻遇的第i个溶洞体积占每个注采井组中所有钻遇溶洞总体积的百分比；V_ei为生产井与钻遇的第i溶洞接触点下方溶洞体积；V_ti为钻遇的第i个溶洞的总体积；n_v1为钻遇溶洞的数目；ω_mj为每个注采井组中第m口生产井所对应的未钻遇的第j个溶洞体积占上述第m口生产井控制的所有未钻遇溶洞总体积的百分比；f_mj为第m口生产井与其所对应的未钻遇的第j个溶洞之间的连通系数，反映生产井与溶洞空间连通程度的大小，与两者的距离和裂缝渗透率有关；Θ(x)为海维塞德函数，决定于溶洞与对应生产井的相对高低位置关系；n_w为每个注采井组中未钻遇溶洞生产井的数目；n_v2为每个注采井组中未被钻遇溶洞的数目；V_tmj为第m口生产井所对应的未钻遇的第j个溶洞的体积；d_vmj为第m口生产井所对应的未钻遇的第j个溶洞的深度；d_wm为第m口生产井的射孔深度。
其中，i为每个注采井组中所有被钻遇溶洞的序号；j为每个注采井组中所有未被 钻遇溶洞的序号；m为上述注采井组中未钻遇溶洞井的序号。
当生产井钻遇溶洞时，V_ei为生产井与钻遇的第i溶洞接触点下方溶洞体积.。例如，图11中的生产井1钻遇了三个溶洞2，生产井1与三个溶洞2的接触点下方的的体积为V_e1、V_e2及V_e3(即图11中虚线下方部分)。如图12所示，当生产井未钻遇溶洞时，生产井3分别通过裂缝5与溶洞2、溶洞3及溶洞5相连通，生产井4分别通过裂缝6与溶洞1、溶洞3、裂缝4相连通。
本实施例中，S2～S5的描述顺序并非用于限定本发明的缝洞型油藏注水开发效果评价方法中各个效果评价的先后顺序，具体实施时，可以根据需要进行调整。
图2为本发明实施例将待评价缝洞单元划分为多个注采井组的流程图，具体包括以下步骤：
S11、根据待评价缝洞单元的地质资料，确定待评价缝洞单元内的缝洞网络，按照缝洞网络中地下河、廊道及断层的走向将待评价缝洞单元初步划分成多个注采井组；
S12、从每个注采井组中选取至少一个生产井，向所选取生产井中投放示踪剂，检测所选取生产井的周边生产井中是否存在示踪剂以确定所选取生产井与其周边生产井的连通关系，根据上述连通关系对初步划分的注采井组进行二次划分；
S13、对于不存在示踪剂的周边生产井，根据包含含水率及压力的生产动态响应，对不存在示踪剂的周边生产井进行注采井组划分。
本发明采用S11～S13对待评价缝洞单元进行注采井组划分，使注采井组的划分可信度大大增加。
将待评价缝洞单元划分成若干个注采井组后，可以根据S2求取各注采井组的井控系数，对各注采井组内生产井与缝洞网络的配置关系进行评价。具体实施时：
首先，将每个注采井组的生产井划分为钻遇溶洞井及未钻遇溶洞井。
其次，根据公式(1)～(5)计算钻遇溶洞井的井控系数A_D及未钻遇溶洞井的井控系数B_D。
再次，根据预定井控系数指标、上述钻遇溶洞井的井控系数A_D及未钻遇溶洞井的井控系数B_D确定各注采井组的井控效果区间。具体实施时，根据上述钻遇溶洞井的井控系数A_D及未钻遇溶洞井的井控系数B_D，利用下述公式(6)确定各注采井组的总井控系数λ_w，根据预定井控系数指标及各所述注采井组的总井控系数λ_w确定各 注采井组的井控效果区间。
λw=Σi=1nv1VtiNtAD+Σm=1nwΣj=1nv2VtmjNtBD---(6)]]>
在上述公式(6)中，N_t为各注采井组内所有储集空间的体积。
上述注采井组的划分涵盖了地质信息、示踪剂监测资料及生产动态响应特征，使注采井组的划分可信度大大增加。
S3具体实施时，可以分别利用下述公式(7)、公式(8)计算各注采井组的水驱采收率ER及吨水换油N_p，并根据预定开发效果指标及所述水驱采收率ER及吨水换油率N_p确定各注采井组的开发效果区间。
ER=Np-NpiN×100%---(7)]]>
Ω=NpBoWpBw×100%---(8)]]>
在上述公式(7)、公式(8)中，N_p为各注采井组的累产油；N_pi为各注采井组注水前的累产油；N为各注采井组的地质储量；B_o为原油的体积系数；W_p为各注采井组的累产水；B_w为水的体积系数。
利用水驱采收率和吨水换油率这两个效果指标实现了对注采井组注水效果好坏的定量评价。
S4具体实施时，可以根据各注采井组的生产动态资料，分别利用公式(9)、公式(10)计算各注采井组的累积含水率f_wc及存水率η，并根据累积含水率f_wc判断是否存在水窜，根据预定存水率指标即存水率η确定各注采井组的存水效果。具体实施时，需要先计算单位时间内各注采井组的累积含水率变化Δf_wc，再求取上述待评价缝洞单元内所有注采井组的累积含水率变化Δf_wc的均值，最后将所述各注采井组累积含水率变化Δf_wc与上述累积含水率变化的均值进行比较，判断各注采井组是否存在水窜。
fwc=WpBwNpBo+WpBw---(9)]]>
η=Wi-Wp-WeWi---(10)]]>
在计算存水率η时，充分考虑注入水直接窜进的无效循环现象和注入水的漏失现 象，采用物质平衡法，利用公式(11)计算注入水漏失：
W_e＝N_pB_o+W_pB_w-NB_oiC_tΔp-W_iB_w   (11)
在上述公式(9)～(11)中，W_i为各注采井组的累积注水量；W_e为各注采井组的漏失量；W_e为各注采井组的漏失量；N为各注采井组的地质储量；B_oi为原始地层压力下原油的体积系数；C_t为综合压缩系数；Δp为各注采井组的压力下降量。
S5具体实施时，可以分别利用公式(12)～(14)计算各注采井组的无因次弹性产量比D_pr、单储压降N_pr及人工水驱指数ψ以确定各注采井组的天然能量区间以及水驱效果区间。
Dpr=ΔpN100Np=Bo100BoiCt-We-WpBw100BoiCtNp---(12)]]>
Npr=NpBo+WpBwNBoiCtΔp=1+WeNBoiCtΔp---(13)]]>
ψ=WiBw-WpBw-W1BwNpBo---(14)]]>
S4、S5中获得的水窜指标、能量指标作为评价各注采井组注水利用率和能量特征的指标，为评价注采井组的注水效果提供了依据。
下面以塔河油田中某一缝洞单元为例，对本发明实施例的缝洞型油藏注水开发效果评价方法进行说明。
在上述缝洞单元内，溶洞类型分为地下河、竖井、廊道各孤立溶洞。其中，地下河具有河流特征，是规模最大的岩溶储集系统，廊道和裂缝是全局分散溶洞和竖井的主要沟通通道。本发明实施例中的缝洞单元内共钻有28口生产井，分别用W1～W28表示。
为评价该缝洞单元的注水开发效果，可将该缝洞单元划分成若干个个注采井组，分别对各注采井组的注水开发效果进行评价，即可获知该缝洞单元的注水开发效果。
首先参照S1，根据该缝洞单元的地质资料、示踪剂检测资料以及生产动态响应，将上述缝洞单元划分成多个注采井组。
具体实施时，结合S11，先根据地质资料确定该缝洞单元内的缝洞网络，按照地下河、廊道和断层(即大裂缝)的走向及分布特征将该缝洞单元初步划分成五个注采井组：注采井组Ⅰ、注采井组Ⅱ、注采井组Ⅲ、注采井组Ⅳ及注采井组Ⅴ，其分布示意图如图3所示。其中，上述五个注采井组中及其所包括的生产井见表1。
表1 按照缝洞单元的地质资料初步划分的注采井组