CN104594854A - 底水油藏水平井变参数射孔控水完井方法及装置 - Google Patents

Abstract

本发明提供了一种底水油藏水平井变参数射孔控水完井方法及装置，涉及水平井控水完井技术领域。方法包括：设置控水单元划分策略和射孔单元划分策略；生成水平井控水井段含水剖面预测模型，确定使各控水单元同时达到目标极限含水率的各控水单元产液量；确定在一预设的射孔枪弹条件下使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数，确定使射孔密度取整误差达到最小的各控水单元的射孔枪弹组合；从而确定各射孔单元的射孔密度；进行底水油藏水平井变参数射孔控水完井操作。本发明能够解决现有技术中的部分水平井段的过早见水或含水过高导致油井含水急剧上升，过早达到极限含水率，而使得未见水或低含水井段得不到有效动用的问题。

Description

底水油藏水平井变参数射孔控水完井方法及装置

技术领域

本发明涉及水平井控水完井技术领域，尤其涉及一种底水油藏水平井变参数射孔控水完井方法及装置。

背景技术

在油气开采技术中，水平井由于具有生产井段长、产层控制面积大、生产压差小等优点，广泛应用于底水油藏开发中。受井筒压力损失和储层非均质等因素影响，水平井开发底水油藏过程中存在局部底水过早突破，油井含水上升过快，储层原油动用不均衡等问题。如何实现稳油控水，最大程度地保证储层均衡动用，已成为底水油藏水平井开发的重要课题之一。

当前国内外应用较广的水平井控水完井技术有：分段变参数控水完井技术、中心管控水完井技术、流入控制装置(Inflow Control Devices，简称ICD)控水完井技术等。其中，中心管控水完井的控水原理是通过改变井筒内流体的流动方向，减轻井筒压力损失对产液不均衡的影响，从而延缓底水脊进。而分段变参数控水完井和ICD控水完井的控水原理则是通过调整完井参数或节流装置参数，使水平井达到均衡的产液剖面或者均衡的见水时间，从而延长无水采油期。

当前的控水完井技术的共同之处在于均以延长见水时间(即无水采油期)为目标，仅适用于未见水的水平井。然而，目前我国大多数油田已进入中高含水开发阶段，水平井产液速率普遍较高，无水采油期普遍较短，部分井段甚至一投产即见水，整个井段见水不均匀问题不可避免。这种部分井段的过早见水或含水过高问题将导致油井含水急剧上升，过早达到极限含水率，而使得未见水或低含水井段得不到有效动用。

发明内容

本发明的实施例提供一种底水油藏水平井变参数射孔控水完井方法及装置，以解决现有技术中的部分水平井段的过早见水或含水过高问题导致的油井含水急剧上升，过早达到极限含水率，而使得未见水或低含水井段得不到有效动用的问题。

为达到上述目的，本发明采用如下技术方案：

一种底水油藏水平井变参数射孔控水完井方法，包括：

设置控水单元划分策略和射孔单元划分策略；

根据所述控水单元划分策略，生成水平井控水井段含水剖面预测模型；

根据所述水平井控水井段含水剖面预测模型，确定使各控水单元同时达到目标极限含水率的各控水单元产液量；

根据所述射孔单元划分策略和所述各控水单元产液量，确定在一预设的射孔枪弹条件下使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数；

根据所述各射孔单元的射孔总表皮系数，确定使射孔密度取整误差达到最小的各控水单元的射孔枪弹组合；

根据所述各射孔单元的射孔总表皮系数和各控水单元的射孔枪弹组合的参数，确定各射孔单元的射孔密度；

根据所述控水单元划分策略、射孔单元划分策略、各控水单元的射孔枪弹组合和各射孔单元的射孔密度，进行水平井变参数射孔控水完井操作。

具体的，所述设置控水单元划分策略和射孔单元划分策略，包括：

设置控水单元的数量和各控水单元的长度；

设置射孔单元的数量和各射孔单元长度。

上述的水平井控水井段含水剖面预测模型，通过一公式一表示为：

$\left\{\begin{array}{c}{\mathrm{WC}}_i\left(t\right)=\frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}}/(1+\frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}})\\ \frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}}=F_{c,i}\×\frac{h_{w,i}}{d_i-h_{w,i}}\\ F_{c,i}=\frac{K_{\mathrm{rw},i}^{\′}{\mathrm{\μ}}_{o,i}}{K_{\mathrm{ro},i}^{\′}{\mathrm{\μ}}_{w,i}}\ln {\left(\frac{X_e{Y}_e{K}_{v,i}}{3d_{i}L{K}_{h,i}}\right)}^{0.5}\\ h_{w,i}=h_{\mathrm{wBT},i}+\frac{\underset{j=1}{\mathrm{\Σ}}{B}_{o,i}{q}_{o,i}^j\mathrm{\Δ}{t}_j}{X_e{Y}_e{\mathrm{\φ}}_i(1-S_{\mathrm{wc},j}-S_{\mathrm{or},i})}\\ h_{\mathrm{wBT},i}=\frac{B_{o,i}{q}_{t,i}{t}_{\mathrm{BT},i}}{X_e{Y}_e{\mathrm{\φ}}_i(1-S_{\mathrm{wc},i}-S_{\mathrm{or},i})}\\ t=t_{\mathrm{BT},i}+\underset{j=1}{\mathrm{\Σ}}\mathrm{\Δ}{t}_j\\ t_{\mathrm{BTi}}=\frac{4{\mathrm{\φ}}_i\mathrm{\Δ}{L}_i{h}_i^2}{B_{o,i}{q}_{t,i}}[1-\cot \frac{\mathrm{\π}{d}_i}{2h_i}\cos \frac{\mathrm{\π}{d}_i}{2h_i}\ln (\sec \frac{\mathrm{\π}{d}_i}{2h_i}+\tan \frac{\mathrm{\π}{d}_i}{2h_i})]\end{array}\right.(i=1,...,m)$

其中，WC_i(t)为第i个控水单元在时间t的含水率；q_w,i为第i个控水单元的产水量；q_o,i为第i个控水单元的产油量；h_w,i为第i个控水单元的平均水侵厚度；K′_rw,i为第i个控水单元的水相端点相对渗透率；K′_ro,i为第i个控水单元的油相端点相对渗透率；μ_w,i为第i个控水单元的水相粘度；μ_o,i为第i个控水单元的油相粘度；K_v,i为第i个控水单元的垂向渗透率；K_h,i为第i个控水单元的水平渗透率；φ_i为第i个控水单元的孔隙度；X_e为泄油宽度；Y_e为泄油长度；L为控水井段长度；S_wc,i为第i个控水单元的束缚水饱和度；S_or,i为第i个控水单元的残余油饱和度；为第i个控水单元在底水突破后第j个时间步的产油量；Δt_j为底水突破后第j个时间步的步长；q_t,i为第i个控水单元在底水突破前的产液量；h_wBT,i为第i个控水单元在底水突破时的平均水侵厚度；t_BT,i为第i个控水单元的底水突破时间；B_w,i为第i个控水单元的水相体积系数；B_o,i为第i个控水单元的油相体积系数；F_c,i为修正系数；h_i为第i个控水单元的拟油柱高度；d_i为第i个控水单元的拟避水高度。

另外，在根据所述控水单元划分策略，生成水平井控水井段含水剖面预测模型之后，包括：

根据预先设置的h_i和d_i的初始值以及所述公式一确定所述t_BT,i和h_wBT,i；

根据所述公式一确定水平井底水突破后的各时间步的含水率WC_i，在所述WC_i等于一初始含水率WC_0,i时，确定所述h_w,i；

根据所述公式一确定所述公式一的第k个迭代步的拟油柱高度和拟避水高度其中， $h_i^k=h_{0,i}+h_{w,i};d_i^k=d_{0,i}+h_{w,i};$

当一值小于一预先设置的绝对误差阈值时，确定所述h_i的值为确定所述d_i的值为

具体的，所述根据所述水平井控水井段含水剖面预测模型，确定使各控水单元同时达到目标极限含水率的各控水单元产液量，包括：

根据所述公式一确定一公式二：

$\left\{\begin{array}{c}{\mathrm{WC}}_1(q_{t,1},t)={\mathrm{WC}}_{\mathrm{limit}}\\ ...\\ {\mathrm{WC}}_i(Q_{t,i},t)={\mathrm{WC}}_{\mathrm{limit}}\\ ...\\ {\mathrm{WC}}_m(q_{t,m},t)={\mathrm{WC}}_{\mathrm{limit}}\end{array}\right.(i=1,...,m)$

其中，WC_limit为所述目标极限含水率；WC_i(q_t,i,t)为第i个控水单元的含水率；

根据所述公式二和一预先设置的公式三确定所述各控水单元的产液量q_t,i；其中，所述公式三为：

$\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{q}_{t,i}=Q_i(i=1,...,m)$

其中Q_t为油井的产液量。

此外，所述根据所述射孔单元划分策略和所述各控水单元产液量，确定在一预设的射孔枪弹条件下使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数，包括：

通过一公式四和公式五确定所述各射孔单元的射孔总表皮系数S_pdc,i(i＝1,2,…,n)；

其中，所述公式四为：

$\left\{\begin{array}{c}\frac{{\mathrm{\μ}}_{o,i+1}{B}_{o,i+1}{q}_{t,i+1}}{2\mathrm{\π}{K}_{h,i+1}\mathrm{\Δl}}(\ln \frac{4{\mathrm{\β}}_{i+1}{h}_{i+1}}{\mathrm{\π}{r}_w}+\ln \mathrm{tg}\frac{\mathrm{\π}{d}_{i+1}}{2h_{i+1}}+S_{\mathrm{pdc},i+1})\\ =\frac{{\mathrm{\μ}}_{o,i}{B}_{o,i}{q}_{t,i}}{2\mathrm{\π}{K}_{h,i}\mathrm{\Δl}}(\ln \frac{4{\mathrm{\β}}_i{h}_i}{\mathrm{\π}{r}_w}+\ln \mathrm{tg}\frac{\mathrm{\π}{d}_i}{2h_i}+S_{\mathrm{pdc},i})-\mathrm{\Δ}{P}_{i,i+1}\\ \mathrm{\Δ}{P}_{i,i+1}=-\frac{{\mathrm{\ρ}}_o\mathrm{f\Δl}}{2D}[\frac{A}{3q_{t,i+1}}(v_{i,i+1}^3-v_{i+1}^3)+\frac{A}{3q_{t,i}}(v_i^3-v_{i,i+1}^3)]\\ v_i=(\underset{j=i}{\overset{n}{\mathrm{\Σ}}}{q}_{t,j}-0.5q_{t,i})/A\\ v_{i,i+1}=\underset{j=i+1}{\overset{n}{\mathrm{\Σ}}}{q}_{t,j}/A\end{array}\right.(i=\mathrm{1,2},...,n-1)$

其中，K_v,i为第i个射孔单元的垂向渗透率；K_h,i为第i个射孔单元的水平渗透率；S_pdc,i为第i个射孔单元的射孔总表皮系数；v_i为第i个射孔单元中点处的截面平均流速；v_i,i+1为第i和i+1个射孔单元交点处的截面平均流速；ρ_o为油相密度；f为管壁摩擦系数；r_w为井筒半径；D为套管内径；A为套管的横截面积；μ_o,i为第i个射孔单元的油相粘度；ΔP_i,i+1为第i和i+1个射孔单元之间的压降；

所述公式五为：

Min(S_pdc,i)＝Min(S_k,min)

其中，S_k,min为预设的第k个射孔枪弹组合的最小射孔总表皮系数。

此外，所述根据所述各射孔单元的射孔总表皮系数，确定使射孔密度取整误差达到最小的各控水单元的射孔枪弹组合，包括：

根据各射孔单元的射孔总表皮系数，确定第i个控水单元的最小射孔总表皮系数 $S_{\mathrm{pdc},\min }^i,(i=1,...,m);$

确定各和S_k,min的差值中的最小非负数对应的射孔枪弹组合为第i个控水单元的射孔枪弹组合。

此外，所述根据所述各射孔单元的射孔总表皮系数和各控水单元的射孔枪弹组合的参数，确定各射孔单元的射孔密度，包括：

根据一公式六确定所述各射孔单元的射孔密度：

$\left\{\begin{array}{c}{S}_{\mathrm{pdc},i}=\frac{1}{K_{\mathrm{dr},i}}{S}_{p,i}+S_{d,i}+S_{\mathrm{cz},i}{l}_{p,i}<r_{d,i}\\ S_{\mathrm{pdc},i}={\hat{S}}_{p,i}+\frac{1}{n_{s,i}{l}_{p,i}}(\frac{1}{K_{\mathrm{cr},i}}-1)\ln \frac{r_{\mathrm{cz},i}}{r_{p,i}},l_{p,i}>r_{d,i}\end{array}\right.$

其中，K_dr,i为第i个射孔单元的污染程度系数；K_cr,i为第i个射孔单元的压实程度系数；S_p,i为第i个射孔单元的射孔几何表皮系数；S_d,i为第i个射孔单元的钻井污染表皮系数；S_cz,i为第i个射孔单元的射孔压实表皮系数；r_d,i为第i个射孔单元的污染带半径；r_cz,i为第i个射孔单元的压实带半径；l_p,i为第i个射孔单元的射孔深度；r_p,i为第i个射孔单元的孔眼半径；n_s,i为第i个射孔单元的射孔密度；为采用有效孔深l′_p，i＝l_p,i-(1-K_dr,i)(r_d,i-r_w)和有效井筒半径r′_w＝r_w+(1-K_dr,i)(r_d,i-r_w)计算的第i个射孔单元的射孔几何表皮系数。

一种底水油藏水平井变参数射孔控水完井装置，包括：

策略设置单元，用于设置控水单元划分策略和射孔单元划分策略；

模型生成单元，用于根据所述控水单元划分策略，生成水平井控水井段含水剖面预测模型；

控水单元产液量确定单元，用于根据所述水平井控水井段含水剖面预测模型，确定使各控水单元同时达到目标极限含水率的各控水单元产液量；

射孔总表皮系数确定单元，用于根据所述射孔单元划分策略和所述各控水单元产液量，确定在一预设的射孔枪弹条件下使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数；

射孔枪弹组合确定单元，用于根据所述各射孔单元的射孔总表皮系数，确定使射孔密度取整误差达到最小的各控水单元的射孔枪弹组合；

射孔密度确定单元，用于根据所述各射孔单元的射孔总表皮系数和各控水单元的射孔枪弹组合的参数，确定各射孔单元的射孔密度；

水平井变参数射孔控水完井操作单元，用于根据所述控水单元划分策略、射孔单元划分策略、各控水单元的射孔枪弹组合和各射孔单元的射孔密度，进行水平井变参数射孔控水完井操作。

另外，该策略设置单元，具体用于：

设置控水单元的数量和各控水单元的长度，以及设置射孔单元的数量和各射孔单元长度。

此外，该模型生成单元生成的水平井控水井段含水剖面预测模型，通过一公式一表示为：

$\left\{\begin{array}{c}{\mathrm{WC}}_i\left(t\right)=\frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}}/(1+\frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}})\\ \frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}}=F_{c,i}\&times;\frac{h_{w,i}}{d_i-h_{w,i}}\\ F_{c,i}=\frac{K_{\mathrm{rw},i}^{\&prime;}{\mathrm{\&mu;}}_{o,i}}{K_{\mathrm{ro},i}^{\&prime;}{\mathrm{\&mu;}}_{w,i}}\ln {\left(\frac{X_e{Y}_e{K}_{v,i}}{3d_{i}L{K}_{h,i}}\right)}^{0.5}\\ h_{w,i}=h_{\mathrm{wBT},i}+\frac{\underset{j=1}{\mathrm{\&Sigma;}}{B}_{o,i}{q}_{o,i}^j\mathrm{\&Delta;}{t}_j}{X_e{Y}_e{\mathrm{\&phi;}}_i(1-S_{\mathrm{wc},j}-S_{\mathrm{or},i})}\\ h_{\mathrm{wBT},i}=\frac{B_{o,i}{q}_{t,i}{t}_{\mathrm{BT},i}}{X_e{Y}_e{\mathrm{\&phi;}}_i(1-S_{\mathrm{wc},i}-S_{\mathrm{or},i})}\\ t=t_{\mathrm{BT},i}+\underset{j=1}{\mathrm{\&Sigma;}}\mathrm{\&Delta;}{t}_j\\ t_{\mathrm{BTi}}=\frac{4{\mathrm{\&phi;}}_i\mathrm{\&Delta;}{L}_i{h}_i^2}{B_{o,i}{q}_{t,i}}[1-\cot \frac{\mathrm{\&pi;}{d}_i}{2h_i}\cos \frac{\mathrm{\&pi;}{d}_i}{2h_i}\ln (\sec \frac{\mathrm{\&pi;}{d}_i}{2h_i}+\tan \frac{\mathrm{\&pi;}{d}_i}{2h_i})]\end{array}\right.(i=1,...,m)$

其中，WC_i(t)为第i个控水单元在时间t的含水率；q_w,i为第i个控水单元的产水量；q_o,i为第i个控水单元的产油量；h_w,i为第i个控水单元的平均水侵厚度；K′_rw,i为第i个控水单元的水相端点相对渗透率；K′_ro,i为第i个控水单元的油相端点相对渗透率；μ_w,i为第i个控水单元的水相粘度；μ_o,i为第i个控水单元的油相粘度；K_v,i为第i个控水单元的垂向渗透率；K_h,i为第i个控水单元的水平渗透率；φ_i为第i个控水单元的孔隙度；X_e为泄油宽度；Y_e为泄油长度；L为控水井段长度；S_wc,i为第i个控水单元的束缚水饱和度；S_or,i为第i个控水单元的残余油饱和度；为第i个控水单元在底水突破后第j个时间步的产油量；Δt_j为底水突破后第j个时间步的步长；q_t,i为第i个控水单元在底水突破前的产液量；h_wBT,i为第i个控水单元在底水突破时的平均水侵厚度；t_BT,i为第i个控水单元的底水突破时间；B_w,i为第i个控水单元的水相体积系数；B_o,i为第i个控水单元的油相体积系数；F_c,i为修正系数；h_i为第i个控水单元的拟油柱高度；d_i为第i个控水单元的拟避水高度。

此外，该底水油藏水平井变参数射孔控水完井装置，用于：

根据预先设置的h_i和d_i的初始值以及所述公式一确定所述t_BT,i和h_wBT,i；

根据所述公式一确定水平井底水突破后的各时间步的含水率WC_i，在所述WC_i等于一初始含水率WC_0,i时，确定所述h_w,i；

根据所述公式一确定所述公式一的第k个迭代步的拟油柱高度和拟避水高度其中， $h_i^k=h_{0,i}+h_{w,i};d_i^k=d_{0,i}+h_{w,i};$

当一值小于一预先设置的绝对误差阈值时，确定所述h_i的值为确定所述d_i的值为

另外，所述控水单元产液量确定单元，具体用于：

根据所述公式一确定一公式二：

$\left\{\begin{array}{c}{\mathrm{WC}}_1(q_{t,1},t)={\mathrm{WC}}_{\mathrm{limit}}\\ ...\\ {\mathrm{WC}}_i(Q_{t,i},t)={\mathrm{WC}}_{\mathrm{limit}}\\ ...\\ {\mathrm{WC}}_m(q_{t,m},t)={\mathrm{WC}}_{\mathrm{limit}}\end{array}\right.(i=1,...,m)$

其中，WC_limit为所述目标极限含水率；WC_i(q_t,i,t)为第i个控水单元的含水率；

根据所述公式二和一预先设置的公式三确定所述各控水单元的产液量q_t,i；其中，所述公式三为：

$\underset{i=1}{\overset{m}{\mathrm{\&Sigma;}}}{q}_{t,i}=Q_i(i=1,...,m)$

其中Q_t为油井的产液量。

另外，所述射孔总表皮系数确定单元，具体用于：

通过一公式四和公式五确定所述各射孔单元的射孔总表皮系数S_pdc,i(i＝1,2,…,n)；

其中，所述公式四为：

$\left\{\begin{array}{c}\frac{{\mathrm{\&mu;}}_{o,i+1}{B}_{o,i+1}{q}_{t,i+1}}{2\mathrm{\&pi;}{K}_{h,i+1}\mathrm{\&Delta;l}}(\ln \frac{4{\mathrm{\&beta;}}_{i+1}{h}_{i+1}}{\mathrm{\&pi;}{r}_w}+\ln \mathrm{tg}\frac{\mathrm{\&pi;}{d}_{i+1}}{2h_{i+1}}+S_{\mathrm{pdc},i+1})\\ =\frac{{\mathrm{\&mu;}}_{o,i}{B}_{o,i}{q}_{t,i}}{2\mathrm{\&pi;}{K}_{h,i}\mathrm{\&Delta;l}}(\ln \frac{4{\mathrm{\&beta;}}_i{h}_i}{\mathrm{\&pi;}{r}_w}+\ln \mathrm{tg}\frac{\mathrm{\&pi;}{d}_i}{2h_i}+S_{\mathrm{pdc},i})-\mathrm{\&Delta;}{P}_{i,i+1}\\ \mathrm{\&Delta;}{P}_{i,i+1}=-\frac{{\mathrm{\&rho;}}_o\mathrm{f\&Delta;l}}{2D}[\frac{A}{3q_{t,i+1}}(v_{i,i+1}^3-v_{i+1}^3)+\frac{A}{3q_{t,i}}(v_i^3-v_{i,i+1}^3)]\\ v_i=(\underset{j=i}{\overset{n}{\mathrm{\&Sigma;}}}{q}_{t,j}-0.5q_{t,i})/A\\ v_{i,i+1}=\underset{j=i+1}{\overset{n}{\mathrm{\&Sigma;}}}{q}_{t,j}/A\end{array}\right.(i=\mathrm{1,2},...,n-1)$

其中，K_v,i为第i个射孔单元的垂向渗透率；K_h,i为第i个射孔单元的水平渗透率；S_pdc,i为第i个射孔单元的射孔总表皮系数；v_i为第i个射孔单元中点处的截面平均流速；v_i,i+1为第i和i+1个射孔单元交点处的截面平均流速；ρ_o为油相密度；f为管壁摩擦系数；r_w为井筒半径；D为套管内径；A为套管的横截面积；μ_o,i为第i个射孔单元的油相粘度；ΔP_i,i+1为第i和i+1个射孔单元之间的压降；

所述公式五为：

Min(S_pdc,i)＝Min(S_k,min)

其中，S_k,min为预设的第k个射孔枪弹组合的最小射孔总表皮系数。

此外，所述射孔枪弹组合确定单元，具体用于：

根据各射孔单元的射孔总表皮系数，确定第i个控水单元的最小射孔总表皮系数 $S_{\mathrm{pdc},\min }^i,(i=1,...,m)$

确定各和S_k,min的差值中的最小非负数对应的射孔枪弹组合为第i个控水单元的射孔枪弹组合。

此外，所述射孔密度确定单元，具体用于：

根据一公式六确定所述各射孔单元的射孔密度：

$\left\{\begin{array}{c}{S}_{\mathrm{pdc},i}=\frac{1}{K_{\mathrm{dr},i}}{S}_{p,i}+S_{d,i}+S_{\mathrm{cz},i}{l}_{p,i}<r_{d,i}\\ S_{\mathrm{pdc},i}={\hat{S}}_{p,i}+\frac{1}{n_{s,i}{l}_{p,i}}(\frac{1}{K_{\mathrm{cr},i}}-1)\ln \frac{r_{\mathrm{cz},i}}{r_{p,i}},l_{p,i}>r_{d,i}\end{array}\right.$

其中，K_dr,i为第i个射孔单元的污染程度系数；K_cr,i为第i个射孔单元的压实程度系数；S_p,i为第i个射孔单元的射孔几何表皮系数；S_d,i为第i个射孔单元的钻井污染表皮系数；S_cz,i为第i个射孔单元的射孔压实表皮系数；r_d,i为第i个射孔单元的污染带半径；r_cz,i为第i个射孔单元的压实带半径；l_p,i为第i个射孔单元的射孔深度；r_p,i为第i个射孔单元的孔眼半径；n_s,i为第i个射孔单元的射孔密度；为采用有效孔深l′_p，i＝l_p,i-(1-K_dr,i)(r_d,i-r_w)和有效井筒半径r′_w＝r_w+(1-K_dr,i)(r_d,i-r_w)计算的第i个射孔单元的射孔几何表皮系数。

本发明实施例提供的底水油藏水平井变参数射孔控水完井方法及装置，通过设置控水单元划分策略和射孔单元划分策略，并根据所述控水单元划分策略，生成水平井控水井段含水剖面预测模型；而后根据所述水平井控水井段含水剖面预测模型，确定使各控水单元同时达到目标极限含水率的各控水单元产液量；从而根据所述射孔单元划分策略和所述各控水单元产液量，确定在一预设的射孔枪弹条件下使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数；从而根据所述各射孔单元的射孔总表皮系数，确定使射孔密度取整误差达到最小的各控水单元的射孔枪弹组合；从而根据所述各射孔单元的射孔总表皮系数和各控水单元的射孔枪弹组合的参数，确定各射孔单元的射孔密度，并根据所述控水单元划分策略、射孔单元划分策略、各控水单元的射孔枪弹组合和各射孔单元的射孔密度，进行水平井变参数射孔控水完井操作。实现了底水油藏早期见水水平井的稳油控水，避免了当前部分水平井段的过早见水或含水过高问题导致的油井含水急剧上升，过早达到极限含水率，而使得未见水或低含水井段得不到有效动用的问题。

附图说明

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

图1为本发明实施例提供的底水油藏水平井变参数射孔控水完井方法的流程图；

图2为本发明实施例中的水平井控水井段的水平渗透率的示意图；

图3为本发明实施例中的水平井控水井段的孔隙度的示意图；

图4为本发明实施例中的水平井控水井段的油柱高度的示意图；

图5为本发明实施例中的水平井控水井段的含水率分布的示意图；

图6为本发明实施例中使各控水单元同时达到目标极限含水率的各控水单元产液量的示意图；

图7为本发明实施例中各射孔单元的射孔总表皮系数的示意图；

图8为本发明实施例中各射孔单元的射孔密度的示意图；

图9为本发明实施例提供的底水油藏水平井变参数射孔控水完井装置的结构示意图。

具体实施方式

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

如图1所示，本发明实施例提供的一种底水油藏水平井变参数射孔控水完井方法，包括：

步骤101、设置控水单元划分策略和射孔单元划分策略。

步骤102、根据所述控水单元划分策略，生成水平井控水井段含水剖面预测模型。

步骤103、根据所述水平井控水井段含水剖面预测模型，确定使各控水单元同时达到目标极限含水率的各控水单元产液量。

步骤104、根据所述射孔单元划分策略和所述各控水单元产液量，确定在一预设的射孔枪弹条件下使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数。

步骤105、根据所述各射孔单元的射孔总表皮系数，确定使射孔密度取整误差达到最小的各控水单元的射孔枪弹组合。

步骤106、根据所述各射孔单元的射孔总表皮系数和各控水单元的射孔枪弹组合的参数，确定各射孔单元的射孔密度。

步骤107、根据所述控水单元划分策略、射孔单元划分策略、各控水单元的射孔枪弹组合和各射孔单元的射孔密度，进行水平井变参数射孔控水完井操作。

本发明实施例提供的底水油藏水平井变参数射孔控水完井方法，通过设置控水单元划分策略和射孔单元划分策略，并根据所述控水单元划分策略，生成水平井控水井段含水剖面预测模型；而后根据所述水平井控水井段含水剖面预测模型，确定使各控水单元同时达到目标极限含水率的各控水单元产液量；从而根据所述射孔单元划分策略和所述各控水单元产液量，确定在一预设的射孔枪弹条件下使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数；从而根据所述各射孔单元的射孔总表皮系数，确定使射孔密度取整误差达到最小的各控水单元的射孔枪弹组合；从而根据所述各射孔单元的射孔总表皮系数和各控水单元的射孔枪弹组合的参数，确定各射孔单元的射孔密度，根据所述控水单元划分策略、射孔单元划分策略、各控水单元的射孔枪弹组合和各射孔单元的射孔密度，进行水平井变参数射孔控水完井操作。实现了底水油藏早期见水水平井的稳油控水，避免了当前部分水平井段的过早见水或含水过高问题导致的油井含水急剧上升，过早达到极限含水率，而使得未见水或低含水井段得不到有效动用的问题。

上述的步骤101中的设置控水单元划分策略和射孔单元划分策略，可以通过如下方式实现：

设置控水单元的数量和各控水单元的长度。例如可以根据储层、流体物性参数，含水率、油柱高度和避水高度，采用聚类分析方法对控水井段进行粗分段，确定控水单元划分策略，即控水单元个数为m和长度为ΔL_i(i＝1,...,m)。

设置射孔单元的数量和各射孔单元长度。例如可以根据射孔枪及其接头长度等对控水井段进行细分段，确定射孔单元划分策略，即射孔单元个数为n和长度为Δl，其中，射孔单元长度Δl为单个射孔枪及其接头长度之和的整数倍。

上述步骤102中的水平井控水井段含水剖面预测模型，通过一公式一表示为：

$\left\{\begin{array}{c}{\mathrm{WC}}_i\left(t\right)=\frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}}/(1+\frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}})\\ \frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}}=F_{c,i}\&times;\frac{h_{w,i}}{d_i-h_{w,i}}\\ F_{c,i}=\frac{K_{\mathrm{rw},i}^{\&prime;}{\mathrm{\&mu;}}_{o,i}}{K_{\mathrm{ro},i}^{\&prime;}{\mathrm{\&mu;}}_{w,i}}\ln {\left(\frac{X_e{Y}_e{K}_{v,i}}{3d_{i}L{K}_{h,i}}\right)}^{0.5}\\ h_{w,i}=h_{\mathrm{wBT},i}+\frac{\underset{j=1}{\mathrm{\&Sigma;}}{B}_{o,i}{q}_{o,i}^j\mathrm{\&Delta;}{t}_j}{X_e{Y}_e{\mathrm{\&phi;}}_i(1-S_{\mathrm{wc},j}-S_{\mathrm{or},i})}\\ h_{\mathrm{wBT},i}=\frac{B_{o,i}{q}_{t,i}{t}_{\mathrm{BT},i}}{X_e{Y}_e{\mathrm{\&phi;}}_i(1-S_{\mathrm{wc},i}-S_{\mathrm{or},i})}\\ t=t_{\mathrm{BT},i}+\underset{j=1}{\mathrm{\&Sigma;}}\mathrm{\&Delta;}{t}_j\\ t_{\mathrm{BTi}}=\frac{4{\mathrm{\&phi;}}_i\mathrm{\&Delta;}{L}_i{h}_i^2}{B_{o,i}{q}_{t,i}}[1-\cot \frac{\mathrm{\&pi;}{d}_i}{2h_i}\cos \frac{\mathrm{\&pi;}{d}_i}{2h_i}\ln (\sec \frac{\mathrm{\&pi;}{d}_i}{2h_i}+\tan \frac{\mathrm{\&pi;}{d}_i}{2h_i})]\end{array}\right.(i=1,...,m)$

其中，WC_i(t)为第i个控水单元在时间t的含水率；q_w,i为第i个控水单元的产水量，单位为m³/s；q_o,i为第i个控水单元的产油量，单位为m³/s；h_w,i为第i个控水单元的平均水侵厚度；K′_rw,i为第i个控水单元的水相端点相对渗透率，单位为m²；K′_ro,i为第i个控水单元的油相端点相对渗透率，单位为m²；μ_w,i为第i个控水单元的水相粘度，单位为mPa·s；μ_o,i为第i个控水单元的油相粘度；K_v,i为第i个控水单元的垂向渗透率；K_h,i为第i个控水单元的水平渗透率；φ_i为第i个控水单元的孔隙度；X_e为泄油宽度；Y_e为泄油长度；L为控水井段长度；S_wc,i为第i个控水单元的束缚水饱和度；S_or,i为第i个控水单元的残余油饱和度；为第i个控水单元在底水突破后第j个时间步的产油量，单位为m³/s；Δt_j为底水突破后第j个时间步的步长，单位为秒；q_t,i为第i个控水单元在底水突破前的产液量；h_wBT,i为第i个控水单元在底水突破时的平均水侵厚度；t_BT,i为第i个控水单元的底水突破时间；B_w,i为第i个控水单元的水相体积系数；B_o,i为第i个控水单元的油相体积系数；F_c,i为修正系数；h_i为第i个控水单元的拟油柱高度；d_i为第i个控水单元的拟避水高度。

另外，在步骤102之后，需要对上述公式一进行求解：

根据预先设置的h_i和d_i的初始值以及所述公式一确定所述t_BT,i和h_wBT,i。

根据所述公式一确定水平井底水突破后的各时间步的含水率WC_i，在所述WC_i等于一初始含水率WC_0,i时，确定所述h_w,i。

根据所述公式一确定所述公式一的第k个迭代步的拟油柱高度和拟避水高度其中， $h_i^k=h_{0,i}+h_{w,i};d_i^k=d_{0,i}+h_{w,i}.$

当一值小于一预先设置的绝对误差阈值时，确定所述h_i的值为确定所述d_i的值为

另外，上述的步骤103中的根据所述水平井控水井段含水剖面预测模型，确定使各控水单元同时达到目标极限含水率的各控水单元产液量，可以通过如下方式实现：

根据所述公式一确定一公式二：

$\left\{\begin{array}{c}{\mathrm{WC}}_1(q_{t,1},t)={\mathrm{WC}}_{\mathrm{limit}}\\ ...\\ {\mathrm{WC}}_i(Q_{t,i},t)={\mathrm{WC}}_{\mathrm{limit}}\\ ...\\ {\mathrm{WC}}_m(q_{t,m},t)={\mathrm{WC}}_{\mathrm{limit}}\end{array}\right.(i=1,...,m)$

其中，WC_limit为所述目标极限含水率；WC_i(q_t,i,t)为第i个控水单元的含水率。

根据所述公式二和一预先设置的公式三确定所述各控水单元的产液量q_t,i；其中，所述公式三为：

$\underset{i=1}{\overset{m}{\mathrm{\&Sigma;}}}{q}_{t,i}=Q_i(i=1,...,m)$

其中Q_t为油井的产液量。

此外，上述步骤104中的根据所述射孔单元划分策略和所述各控水单元产液量，确定在一预设的射孔枪弹条件下使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数，可以通过如下方式实现：

通过一公式四和公式五确定所述各射孔单元的射孔总表皮系数S_pdc,i(i＝1,2,…,n)。

其中，所述公式四为：

$\left\{\begin{array}{c}\frac{{\mathrm{\&mu;}}_{o,i+1}{B}_{o,i+1}{q}_{t,i+1}}{2\mathrm{\&pi;}{K}_{h,i+1}\mathrm{\&Delta;l}}(\ln \frac{4{\mathrm{\&beta;}}_{i+1}{h}_{i+1}}{\mathrm{\&pi;}{r}_w}+\ln \mathrm{tg}\frac{\mathrm{\&pi;}{d}_{i+1}}{2h_{i+1}}+S_{\mathrm{pdc},i+1})\\ =\frac{{\mathrm{\&mu;}}_{o,i}{B}_{o,i}{q}_{t,i}}{2\mathrm{\&pi;}{K}_{h,i}\mathrm{\&Delta;l}}(\ln \frac{4{\mathrm{\&beta;}}_i{h}_i}{\mathrm{\&pi;}{r}_w}+\ln \mathrm{tg}\frac{\mathrm{\&pi;}{d}_i}{2h_i}+S_{\mathrm{pdc},i})-\mathrm{\&Delta;}{P}_{i,i+1}\\ \mathrm{\&Delta;}{P}_{i,i+1}=-\frac{{\mathrm{\&rho;}}_o\mathrm{f\&Delta;l}}{2D}[\frac{A}{3q_{t,i+1}}(v_{i,i+1}^3-v_{i+1}^3)+\frac{A}{3q_{t,i}}(v_i^3-v_{i,i+1}^3)]\\ v_i=(\underset{j=i}{\overset{n}{\mathrm{\&Sigma;}}}{q}_{t,j}-0.5q_{t,i})/A\\ v_{i,i+1}=\underset{j=i+1}{\overset{n}{\mathrm{\&Sigma;}}}{q}_{t,j}/A\end{array}\right.(i=\mathrm{1,2},...,n-1)$

其中，K_v,i为第i个射孔单元的垂向渗透率；K_h,i为第i个射孔单元的水平渗透率；S_pdc,i为第i个射孔单元的射孔总表皮系数；v_i为第i个射孔单元中点处的截面平均流速；v_i,i+1为第i和i+1个射孔单元交点处的截面平均流速；ρ_o为油相密度；f为管壁摩擦系数；r_w为井筒半径；D为套管内径；A为套管的横截面积；μ_o,i为第i个射孔单元的油相粘度；ΔP_i,i+1为第i和i+1个射孔单元之间的压降。

所述公式五为：

Min(S_pdc,i)＝Min(S_k,min)

其中，S_k,min为预设的第k个射孔枪弹组合的最小射孔总表皮系数。例如，a种射孔枪型与b种射孔弹型经随机组合后可形成a×b种射孔枪弹组合，各种射孔枪弹组合的最小射孔总表皮为S_k,min，k＝1,2,…,a×b，表皮系数约束条件为上述公式五。

此外，上述的步骤105中的根据所述各射孔单元的射孔总表皮系数，确定使射孔密度取整误差达到最小的各控水单元的射孔枪弹组合，可以通过如下方式实现：

根据各射孔单元的射孔总表皮系数，确定第i个控水单元的最小射孔总表皮系数 $S_{\mathrm{pdc},\min }^i,(i=1,..,m).$

确定各和S_k,min的差值中的最小非负数对应的射孔枪弹组合为第i个控水单元的射孔枪弹组合。

此外，上述步骤106中的根据所述各射孔单元的射孔总表皮系数和各控水单元的射孔枪弹组合的参数，确定各射孔单元的射孔密度，可以通过如下方式实现：

根据一公式六确定所述各射孔单元的射孔密度：

$\left\{\begin{array}{c}{S}_{\mathrm{pdc},i}=\frac{1}{K_{\mathrm{dr},i}}{S}_{p,i}+S_{d,i}+S_{\mathrm{cz},i}{l}_{p,i}<r_{d,i}\\ S_{\mathrm{pdc},i}={\hat{S}}_{p,i}+\frac{1}{n_{s,i}{l}_{p,i}}(\frac{1}{K_{\mathrm{cr},i}}-1)\ln \frac{r_{\mathrm{cz},i}}{r_{p,i}},l_{p,i}>r_{d,i}\end{array}\right.$

其中，K_dr,i为第i个射孔单元的污染程度系数；K_cr,i为第i个射孔单元的压实程度系数；S_p,i为第i个射孔单元的射孔几何表皮系数；S_d,i为第i个射孔单元的钻井污染表皮系数；S_cz,i为第i个射孔单元的射孔压实表皮系数；r_d,i为第i个射孔单元的污染带半径；r_cz,i为第i个射孔单元的压实带半径；l_p,i为第i个射孔单元的射孔深度；r_p,i为第i个射孔单元的孔眼半径；n_s,i为第i个射孔单元的射孔密度；为采用有效孔深l′_p，i＝l_p,i-(1-K_dr,i)(r_d,i-r_w)和有效井筒半径r′_w＝r_w+(1-K_dr,i)(r_d,i-r_w)计算的第i个射孔单元的射孔几何表皮系数。

各射孔单元的射孔密度可通过求解得到的n_s,i，i＝1,...,n，四舍五入得到。

下面以某底水油藏一口早期见水的水平井为例，说明本发明的底水油藏水平井变参数射孔控水完井方法。

该水平井控水井段长度400m，井筒半径0.15m，套管内径121mm，泄油长度400m，泄油宽度500m，避水高度6m。其井段水平渗透率、孔隙度、油柱高度和含水率分布分别如图2-5所示。地层原油粘度为1.1mPa·s，地层水粘度为0.4mPa·s，原油体积系数1.3，水体积系数1，残余油饱和度0.3，束缚水饱和度0.35，水相端点相对渗透率0.4，油相端点相对渗透率0.72，垂向渗透率与水平渗透率之比为0.5。现有一种102射孔枪A和四种102射孔弹B1、B2、B3、B4可供选择。其中射孔枪长度3m，射孔枪接头长度0.3m，最大孔密为40孔每米。弹型A1的孔深221.6mm、孔径17mm，弹型A2的孔深371.5mm、孔径9.9mm，弹型A3的孔深483mm、孔径10.92mm，弹型A4的孔深643.3mm、孔径13.4mm。对该井进行变参数射孔控水设计的具体实施方式为：控水单元划分策略如图2-5所示，整个控水井段共划分7个控水单元，长度分别为30m、70m、90m、50m、70m、50m、40m。射孔单元划分策略为：单元长度10m、个数40个。使各控水单元同时达到目标极限含水率的各控水单元产液量如图6所示。使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数如图7所示。而使孔密取整误差达到最小的各控水单元射孔枪弹组合和各射孔单元的射孔密度如图8所示。

对应于图1所述的底水油藏水平井变参数射孔控水完井方法的实施例，本发明实施例提供了一种底水油藏水平井变参数射孔控水完井装置，如图9所示，包括：

策略设置单元21，可以设置控水单元划分策略和射孔单元划分策略。

模型生成单元22，可以根据所述控水单元划分策略，生成水平井控水井段含水剖面预测模型。

控水单元产液量确定单元23，可以根据所述水平井控水井段含水剖面预测模型，确定使各控水单元同时达到目标极限含水率的各控水单元产液量。

射孔总表皮系数确定单元24，可以根据所述射孔单元划分策略和所述各控水单元产液量，确定在一预设的射孔枪弹条件下使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数。

射孔枪弹组合确定单元25，可以根据所述各射孔单元的射孔总表皮系数，确定使射孔密度取整误差达到最小的各控水单元的射孔枪弹组合。

射孔密度确定单元26，可以根据所述各射孔单元的射孔总表皮系数和各控水单元的射孔枪弹组合的参数，确定各射孔单元的射孔密度。

水平井变参数射孔控水完井操作单元27，可以根据所述控水单元划分策略、射孔单元划分策略、各控水单元的射孔枪弹组合和各射孔单元的射孔密度，进行水平井变参数射孔控水完井操作。

另外，该策略设置单元21，具体可以设置控水单元的数量和各控水单元的长度，以及设置射孔单元的数量和各射孔单元长度。

此外，该模型生成单元22生成的水平井控水井段含水剖面预测模型，可以通过一公式一表示为：

$\left\{\begin{array}{c}{\mathrm{WC}}_i\left(t\right)=\frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}}/(1+\frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}})\\ \frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}}=F_{c,i}\&times;\frac{h_{w,i}}{d_i-h_{w,i}}\\ F_{c,i}=\frac{K_{\mathrm{rw},i}^{\&prime;}{\mathrm{\&mu;}}_{o,i}}{K_{\mathrm{ro},i}^{\&prime;}{\mathrm{\&mu;}}_{w,i}}\ln {\left(\frac{X_e{Y}_e{K}_{v,i}}{3d_{i}L{K}_{h,i}}\right)}^{0.5}\\ h_{w,i}=h_{\mathrm{wBT},i}+\frac{\underset{j=1}{\mathrm{\&Sigma;}}{B}_{o,i}{q}_{o,i}^j\mathrm{\&Delta;}{t}_j}{X_e{Y}_e{\mathrm{\&phi;}}_i(1-S_{\mathrm{wc},j}-S_{\mathrm{or},i})}\\ h_{\mathrm{wBT},i}=\frac{B_{o,i}{q}_{t,i}{t}_{\mathrm{BT},i}}{X_e{Y}_e{\mathrm{\&phi;}}_i(1-S_{\mathrm{wc},i}-S_{\mathrm{or},i})}\\ t=t_{\mathrm{BT},i}+\underset{j=1}{\mathrm{\&Sigma;}}\mathrm{\&Delta;}{t}_j\\ t_{\mathrm{BTi}}=\frac{4{\mathrm{\&phi;}}_i\mathrm{\&Delta;}{L}_i{h}_i^2}{B_{o,i}{q}_{t,i}}[1-\cot \frac{\mathrm{\&pi;}{d}_i}{2h_i}\cos \frac{\mathrm{\&pi;}{d}_i}{2h_i}\ln (\sec \frac{\mathrm{\&pi;}{d}_i}{2h_i}+\tan \frac{\mathrm{\&pi;}{d}_i}{2h_i})]\end{array}\right.(i=1,...,m)$

其中，WC_i(t)为第i个控水单元在时间t的含水率；q_w,i为第i个控水单元的产水量；q_o,i为第i个控水单元的产油量；h_w,i为第i个控水单元的平均水侵厚度；K′_rw,i为第i个控水单元的水相端点相对渗透率；K′_ro,i为第i个控水单元的油相端点相对渗透率；μ_w,i为第i个控水单元的水相粘度；μ_o,i为第i个控水单元的油相粘度；K_v,i为第i个控水单元的垂向渗透率；K_h,i为第i个控水单元的水平渗透率；φ_i为第i个控水单元的孔隙度；X_e为泄油宽度；Y_e为泄油长度；L为控水井段长度；S_wc,i为第i个控水单元的束缚水饱和度；S_or,i为第i个控水单元的残余油饱和度；为第i个控水单元在底水突破后第j个时间步的产油量；Δt_j为底水突破后第j个时间步的步长；q_t,i为第i个控水单元在底水突破前的产液量；h_wBT,i为第i个控水单元在底水突破时的平均水侵厚度；t_BT,i为第i个控水单元的底水突破时间；B_w,i为第i个控水单元的水相体积系数；B_o,i为第i个控水单元的油相体积系数；F_c,i为修正系数；h_i为第i个控水单元的拟油柱高度；d_i为第i个控水单元的拟避水高度。

此外，该底水油藏水平井变参数射孔控水完井装置，可以根据预先设置的h_i和d_i的初始值以及所述公式一确定所述t_BT,i和h_wBT,i。根据所述公式一确定水平井底水突破后的各时间步的含水率WC_i，在所述WC_i等于一初始含水率WC_0,i时，确定所述h_w,i。根据所述公式一确定所述公式一的第k个迭代步的拟油柱高度和拟避水高度其中， $h_i^k=h_{0,i}+h_{w,i};d_i^k=d_{0,i}+h_{w,i}.$ 当一 $|h_i^k-h_i^{k-1}|$ 值小于一预先设置的绝对误差阈值时，确定所述h_i的值为确定所述d_i的值为

另外，所述控水单元产液量确定单元23，具体可以根据所述公式一确定一公式二：

$\left\{\begin{array}{c}{\mathrm{WC}}_1(q_{t,1},t)={\mathrm{WC}}_{\mathrm{limit}}\\ ...\\ {\mathrm{WC}}_i(Q_{t,i},t)={\mathrm{WC}}_{\mathrm{limit}}\\ ...\\ {\mathrm{WC}}_m(q_{t,m},t)={\mathrm{WC}}_{\mathrm{limit}}\end{array}\right.(i=1,...,m)$

其中，WC_limit为所述目标极限含水率；WC_i(q_t,i,t)为第i个控水单元的含水率。

根据所述公式二和一预先设置的公式三确定所述各控水单元的产液量q_t,i；其中，所述公式三为：

$\underset{i=1}{\overset{m}{\mathrm{\&Sigma;}}}{q}_{t,i}=Q_i(i=1,...,m)$

其中Q_t为油井的产液量。

另外，所述射孔总表皮系数确定单元24，可以通过一公式四和公式五确定所述各射孔单元的射孔总表皮系数S_pdc,i(i＝1,2,…,n)。

其中，所述公式四为：

$\left\{\begin{array}{c}\frac{{\mathrm{\&mu;}}_{o,i+1}{B}_{o,i+1}{q}_{t,i+1}}{2\mathrm{\&pi;}{K}_{h,i+1}\mathrm{\&Delta;l}}(\ln \frac{4{\mathrm{\&beta;}}_{i+1}{h}_{i+1}}{\mathrm{\&pi;}{r}_w}+\ln \mathrm{tg}\frac{\mathrm{\&pi;}{d}_{i+1}}{2h_{i+1}}+S_{\mathrm{pdc},i+1})\\ =\frac{{\mathrm{\&mu;}}_{o,i}{B}_{o,i}{q}_{t,i}}{2\mathrm{\&pi;}{K}_{h,i}\mathrm{\&Delta;l}}(\ln \frac{4{\mathrm{\&beta;}}_i{h}_i}{\mathrm{\&pi;}{r}_w}+\ln \mathrm{tg}\frac{\mathrm{\&pi;}{d}_i}{2h_i}+S_{\mathrm{pdc},i})-\mathrm{\&Delta;}{P}_{i,i+1}\\ \mathrm{\&Delta;}{P}_{i,i+1}=-\frac{{\mathrm{\&rho;}}_o\mathrm{f\&Delta;l}}{2D}[\frac{A}{3q_{t,i+1}}(v_{i,i+1}^3-v_{i+1}^3)+\frac{A}{3q_{t,i}}(v_i^3-v_{i,i+1}^3)]\\ v_i=(\underset{j=i}{\overset{n}{\mathrm{\&Sigma;}}}{q}_{t,j}-0.5q_{t,i})/A\\ v_{i,i+1}=\underset{j=i+1}{\overset{n}{\mathrm{\&Sigma;}}}{q}_{t,j}/A\end{array}\right.(i=\mathrm{1,2},...,n-1)$

其中，K_v,i为第i个射孔单元的垂向渗透率；K_h,i为第i个射孔单元的水平渗透率；S_pdc,i为第i个射孔单元的射孔总表皮系数；v_i为第i个射孔单元中点处的截面平均流速；v_i,i+1为第i和i+1个射孔单元交点处的截面平均流速；ρ_o为油相密度；f为管壁摩擦系数；r_w为井筒半径；D为套管内径；A为套管的横截面积；μ_o,i为第i个射孔单元的油相粘度；ΔP_i,i+1为第i和i+1个射孔单元之间的压降。

所述公式五为：

Min(S_pdc,i)＝Min(S_k,min)

其中，S_k,min为预设的第k个射孔枪弹组合的最小射孔总表皮系数。

此外，所述射孔枪弹组合确定单元25，可以根据各射孔单元的射孔总表皮系数，确定第i个控水单元的最小射孔总表皮系数确定各和S_k,min的差值中的最小非负数对应的射孔枪弹组合为第i个控水单元的射孔枪弹组合。

此外，所述射孔密度确定单元26，可以根据一公式六确定所述各射孔单元的射孔密度：

$\left\{\begin{array}{c}{S}_{\mathrm{pdc},i}=\frac{1}{K_{\mathrm{dr},i}}{S}_{p,i}+S_{d,i}+S_{\mathrm{cz},i}{l}_{p,i}<r_{d,i}\\ S_{\mathrm{pdc},i}={\hat{S}}_{p,i}+\frac{1}{n_{s,i}{l}_{p,i}}(\frac{1}{K_{\mathrm{cr},i}}-1)\ln \frac{r_{\mathrm{cz},i}}{r_{p,i}},l_{p,i}>r_{d,i}\end{array}\right.$

其中，K_dr,i为第i个射孔单元的污染程度系数；K_cr,i为第i个射孔单元的压实程度系数；S_p,i为第i个射孔单元的射孔几何表皮系数；S_d,i为第i个射孔单元的钻井污染表皮系数；S_cz,i为第i个射孔单元的射孔压实表皮系数；r_d,i为第i个射孔单元的污染带半径；r_cz,i为第i个射孔单元的压实带半径；l_p,i为第i个射孔单元的射孔深度；r_p,i为第i个射孔单元的孔眼半径；n_s,i为第i个射孔单元的射孔密度；为采用有效孔深l′_p，i＝l_p,i-(1-K_dr,i)(r_d,i-r_w)和有效井筒半径r′_w＝r_w+(1-K_dr,i)(r_d,i-r_w)计算的第i个射孔单元的射孔几何表皮系数。

值得说明的是，本发明实施例提供的底水油藏水平井变参数射孔控水完井装置的具体实现方式可以参见图1中的方法实施例，此处不再赘述。

本发明实施例提供的底水油藏水平井变参数射孔控水完井装置，通过设置控水单元划分策略和射孔单元划分策略，并根据所述控水单元划分策略，生成水平井控水井段含水剖面预测模型；而后根据所述水平井控水井段含水剖面预测模型，确定使各控水单元同时达到目标极限含水率的各控水单元产液量；从而根据所述射孔单元划分策略和所述各控水单元产液量，确定在一预设的射孔枪弹条件下使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数；从而根据所述各射孔单元的射孔总表皮系数，确定使射孔密度取整误差达到最小的各控水单元的射孔枪弹组合；从而根据所述各射孔单元的射孔总表皮系数和各控水单元的射孔枪弹组合的参数，确定各射孔单元的射孔密度，并根据所述控水单元划分策略、射孔单元划分策略、各控水单元的射孔枪弹组合和各射孔单元的射孔密度，进行水平井变参数射孔控水完井操作。实现了底水油藏早期见水水平井的稳油控水，避免了当前部分水平井段的过早见水或含水过高问题导致的油井含水急剧上升，过早达到极限含水率，而使得未见水或低含水井段得不到有效动用的问题。

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

Claims (16)

1.一种底水油藏水平井变参数射孔控水完井方法，其特征在于，包括：

设置控水单元划分策略和射孔单元划分策略；

根据所述控水单元划分策略，生成水平井控水井段含水剖面预测模型；

根据所述水平井控水井段含水剖面预测模型，确定使各控水单元同时达到目标极限含水率的各控水单元产液量；

根据所述射孔单元划分策略和所述各控水单元产液量，确定在一预设的射孔枪弹条件下使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数；

根据所述各射孔单元的射孔总表皮系数，确定使射孔密度取整误差达到最小的各控水单元的射孔枪弹组合；

根据所述各射孔单元的射孔总表皮系数和各控水单元的射孔枪弹组合的参数，确定各射孔单元的射孔密度；

根据所述控水单元划分策略、射孔单元划分策略、各控水单元的射孔枪弹组合和各射孔单元的射孔密度，进行水平井变参数射孔控水完井操作。

2.根据权利要求1所述的底水油藏水平井变参数射孔控水完井方法，其特征在于，所述设置控水单元划分策略和射孔单元划分策略，包括：

设置控水单元的数量和各控水单元的长度；

设置射孔单元的数量和各射孔单元长度。

3.根据权利要求2所述的底水油藏水平井变参数射孔控水完井方法，其特征在于，所述水平井控水井段含水剖面预测模型，通过一公式一表示为：

$\left\{\begin{array}{c}{\mathrm{WC}}_i\left(t\right)=\frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}}/(1+\frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}})\\ \frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}}=F_{c,i}\&times;\frac{h_{w,i}}{d_i-h_{w,i}}\\ F_{c,i}=\frac{K_{\mathrm{rw},i}^{\&prime;}{\mathrm{\&mu;}}_{o,i}}{K_{\mathrm{ro},i}^{\&prime;}{\mathrm{\&mu;}}_{w,i}}\ln {\left(\frac{X_e{Y}_e{K}_{v,i}}{{3d}_{i}L{K}_{h,i}}\right)}^{0.5}\\ h_{w,i}=h_{\mathrm{wBT},i}+\frac{\underset{j=1}{\mathrm{\&Sigma;}}{B}_{o,i}{q}_{o,i}^j{\mathrm{\&Delta;t}}_j}{X_e{Y}_e{\mathrm{\&phi;}}_i(1-S_{\mathrm{wc},i}-S_{\mathrm{or},i})}\\ h_{\mathrm{wBT},i}=\frac{B_{o,i}{q}_{t,i}{t}_{\mathrm{BT},i}}{X_e{Y}_e{\mathrm{\&phi;}}_i(1-S_{\mathrm{wc},i}-S_{\mathrm{or},i})}\\ t=t_{\mathrm{BT},i}+\underset{j=1}{\mathrm{\&Sigma;}}{\mathrm{\&Delta;t}}_j\\ t_{\mathrm{Bt},i}=\frac{{4\mathrm{\&phi;}}_i{\mathrm{\&Delta;L}}_i{h_i}^2}{B_{o,i}{q}_{t,i}}[1-\cot \frac{{\mathrm{\&pi;d}}_i}{2h_i}\cos \frac{{\mathrm{\&pi;d}}_i}{2h_i}\ln (\sec \frac{{\mathrm{\&pi;d}}_i}{2h_i}+\tan \frac{{\mathrm{\&pi;d}}_i}{2h_i})]\end{array}\right.(i=1,...,m)$

其中，WC_i(t)为第i个控水单元在时间t的含水率；q_w,i为第i个控水单元的产水量；q_o,i为第i个控水单元的产油量；h_w,i为第i个控水单元的平均水侵厚度；K′_rw,i为第i个控水单元的水相端点相对渗透率；K′_ro,i为第i个控水单元的油相端点相对渗透率；μ_w,i为第i个控水单元的水相粘度；μ_o,i为第i个控水单元的油相粘度；K_v,i为第i个控水单元的垂向渗透率；K_h,i为第i个控水单元的水平渗透率；φ_i为第i个控水单元的孔隙度；X_e为泄油宽度；Y_e为泄油长度；L为控水井段长度；S_wc,i为第i个控水单元的束缚水饱和度；S_or,i为第i个控水单元的残余油饱和度；为第i个控水单元在底水突破后第j个时间步的产油量；Δt_j为底水突破后第j个时间步的步长；q_t,i为第i个控水单元在底水突破前的产液量；h_wBT,i为第i个控水单元在底水突破时的平均水侵厚度；t_BT,i为第i个控水单元的底水突破时间；B_w,i为第i个控水单元的水相体积系数；B_o,i为第i个控水单元的油相体积系数；F_c,i为修正系数；h_i为第i个控水单元的拟油柱高度；d_i为第i个控水单元的拟避水高度。

4.根据权利要求3所述的底水油藏水平井变参数射孔控水完井方法，其特征在于，在根据所述控水单元划分策略，生成水平井控水井段含水剖面预测模型之后，包括：

根据预先设置的h_i和d_i的初始值以及所述公式一确定所述t_BT,i和h_wBT,i；

根据所述公式一确定水平井底水突破后的各时间步的含水率WC_i，在所述WC_i等于一初始含水率WC_0,i时，确定所述h_w,i；

根据所述公式一确定所述公式一的第k个迭代步的拟油柱高度和拟避水高度其中， $h_i^k=h_{0.i}+h_{w,i};d_i^k=d_{0,i}+h_{w,i};$

当一值小于一预先设置的绝对误差阈值时，确定所述h_i的值为确定所述d_i的值为

5.根据权利要求4所述的底水油藏水平井变参数射孔控水完井方法，其特征在于，所述根据所述水平井控水井段含水剖面预测模型，确定使各控水单元同时达到目标极限含水率的各控水单元产液量，包括：

根据所述公式一确定一公式二：

$\left\{\begin{array}{c}{\mathrm{WC}}_1(q_{t,1},t)={\mathrm{WC}}_{\mathrm{limit}}\\ ...\\ {\mathrm{WC}}_i(q_{t,i},t)={\mathrm{WC}}_{\mathrm{linit}}\\ ...\\ {\mathrm{WC}}_m(q_{t,m},t)={\mathrm{WC}}_{\mathrm{limit}}\end{array}\right.(i=1,...,m)$

其中，WC_limit为所述目标极限含水率；WC_i(q_t,i,t)为第i个控水单元的含水率；

根据所述公式二和一预先设置的公式三确定所述各控水单元的产液量q_t,i；其中，所述公式三为：

$\underset{i=1}{\overset{m}{\mathrm{\&Sigma;}}}{q}_{t,i}=Q_t,(i=1,...,m)$

其中Q_t为油井的产液量。

6.根据权利要求5所述的底水油藏水平井变参数射孔控水完井方法，其特征在于，所述根据所述射孔单元划分策略和所述各控水单元产液量，确定在一预设的射孔枪弹条件下使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数，包括：

通过一公式四和公式五确定所述各射孔单元的射孔总表皮系数S_pdc,i(i＝1,2,…,n)；

其中，所述公式四为：

$\left\{\begin{array}{c}\frac{{\mathrm{\&mu;}}_{o,i+1}{B}_{o.i+1}{q}_{t,i+1}}{2\mathrm{\&pi;}{K}_{h,i+1}\mathrm{\&Delta;l}}(\ln \frac{4{\mathrm{\&beta;}}_{i+1}{h}_{i+1}}{\mathrm{\&pi;}{r}_w}+\ln \mathrm{tg}\frac{{\mathrm{\&pi;d}}_{i+1}}{2h_{i+1}}+S_{\mathrm{pdc},i+1})\\ =\frac{{\mathrm{\&mu;}}_{o,i}{B}_{o,i}{q}_{t,i}}{2\mathrm{\&pi;}{K}_{h,i}\mathrm{\&Delta;l}}(\ln \frac{4{\mathrm{\&beta;}}_i{h}_i}{{\mathrm{\&pi;r}}_w}+\ln \mathrm{tg}\frac{\mathrm{\&pi;}{d}_i}{2h_i}+S_{\mathrm{pdc},i})-{\mathrm{\&Delta;P}}_{i,i+1}\\ \mathrm{\&Delta;}{P}_{i,i+1}=-\frac{{\mathrm{\&rho;}}_o\mathrm{f\&Delta;l}}{2D}[\frac{A}{3q_{t,i+1}}(v_{i,i+1}^3-v_{i+1}^3)+\frac{A}{{3q}_{t,i}}(v_i^3-v_{i,i+1}^3)]\\ v_i=(\underset{j=i}{\overset{n}{\mathrm{\&Sigma;}}}{q}_{t,j}-0.5q_{t,i})/A\\ v_{i,i+1}=\underset{j=i+1}{\overset{n}{\mathrm{\&Sigma;}}}{q}_{t,j}/A\end{array}\right.(i=\mathrm{1,2},...,n-1)$

其中，K_v,i为第i个射孔单元的垂向渗透率；K_h,i为第i个射孔单元的水平渗透率；S_pdc,i为第i个射孔单元的射孔总表皮系数；v_i为第i个射孔单元中点处的截面平均流速；v_i,i+1为第i和i+1个射孔单元交点处的截面平均流速；ρ_o为油相密度；f为管壁摩擦系数；r_w为井筒半径；D为套管内径；A为套管的横截面积；μ_o,i为第i个射孔单元的油相粘度；ΔP_i,i+1为第i和i+1个射孔单元之间的压降；

所述公式五为：

Min(S_pdc,i)＝Min(S_k,min)

其中，S_k,min为预设的第k个射孔枪弹组合的最小射孔总表皮系数。

7.根据权利要求6所述的底水油藏水平井变参数射孔控水完井方法，其特征在于，所述根据所述各射孔单元的射孔总表皮系数，确定使射孔密度取整误差达到最小的各控水单元的射孔枪弹组合，包括：

根据各射孔单元的射孔总表皮系数，确定第i个控水单元的最小射孔总表皮系数 $S_{\mathrm{pdc},\min }^i,(i=1,...,m);$

确定各和S_k,min的差值中的最小非负数对应的射孔枪弹组合为第i个控水单元的射孔枪弹组合。

8.根据权利要求7所述的底水油藏水平井变参数射孔控水完井方法，其特征在于，所述根据所述各射孔单元的射孔总表皮系数和各控水单元的射孔枪弹组合的参数，确定各射孔单元的射孔密度，包括：

根据一公式六确定所述各射孔单元的射孔密度：

$\left\{\begin{array}{c}{S}_{\mathrm{pdc},i}=\frac{1}{K_{\mathrm{dr},i}}{S}_{p,i}+S_{d,i}+S_{\mathrm{cz},i},l_{p,i}<r_{d,i}\\ S_{\mathrm{pdc},i}={\hat{S}}_{p,i}+\frac{1}{n_{s,i}{l}_{p,i}}(\frac{1}{K_{\mathrm{cr},i}}-1)\ln \frac{r_{\mathrm{cz},i}}{r_{p,i}},l_{p,i}>r_{d,i}\end{array}\right.$

其中，K_dr,i为第i个射孔单元的污染程度系数；K_cr,i为第i个射孔单元的压实程度系数；S_p,i为第i个射孔单元的射孔几何表皮系数；S_d,i为第i个射孔单元的钻井污染表皮系数；S_cz,i为第i个射孔单元的射孔压实表皮系数；r_d,i为第i个射孔单元的污染带半径；r_cz,i为第i个射孔单元的压实带半径；l_p,i为第i个射孔单元的射孔深度；r_p,i为第i个射孔单元的孔眼半径；n_s,i为第i个射孔单元的射孔密度；为采用有效孔深l_p,i′＝l_p,i-(1-K_dr,i)(r_d,i-r_w)和有效井筒半径r_w′＝r_w+(1-K_dr,i)(r_d,i-r_w)计算的第i个射孔单元的射孔几何表皮系数。

9.一种底水油藏水平井变参数射孔控水完井装置，其特征在于，包括：

策略设置单元，用于设置控水单元划分策略和射孔单元划分策略；

模型生成单元，用于根据所述控水单元划分策略，生成水平井控水井段含水剖面预测模型；

控水单元产液量确定单元，用于根据所述水平井控水井段含水剖面预测模型，确定使各控水单元同时达到目标极限含水率的各控水单元产液量；

射孔总表皮系数确定单元，用于根据所述射孔单元划分策略和所述各控水单元产液量，确定在一预设的射孔枪弹条件下使近井附加渗流压降达到最小的各射孔单元的射孔总表皮系数；

射孔枪弹组合确定单元，用于根据所述各射孔单元的射孔总表皮系数，确定使射孔密度取整误差达到最小的各控水单元的射孔枪弹组合；

射孔密度确定单元，用于根据所述各射孔单元的射孔总表皮系数和各控水单元的射孔枪弹组合的参数，确定各射孔单元的射孔密度；

水平井变参数射孔控水完井操作单元，用于根据所述控水单元划分策略、射孔单元划分策略、各控水单元的射孔枪弹组合和各射孔单元的射孔密度，进行水平井变参数射孔控水完井操作。

10.根据权利要求9所述的底水油藏水平井变参数射孔控水完井装置，其特征在于，所述策略设置单元，具体用于：

设置控水单元的数量和各控水单元的长度，以及设置射孔单元的数量和各射孔单元长度。

11.根据权利要求10所述的底水油藏水平井变参数射孔控水完井装置，其特征在于，所述模型生成单元生成的水平井控水井段含水剖面预测模型，通过一公式一表示为：

$\left\{\begin{array}{c}{\mathrm{WC}}_i\left(t\right)=\frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}}/(1+\frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}})\\ \frac{B_{w,i}{q}_{w,i}}{B_{o,i}{q}_{o,i}}=F_{c,i}\&times;\frac{h_{w,i}}{d_i-h_{w,i}}\\ F_{c,i}=\frac{K_{\mathrm{rw},i}^{\&prime;}{\mathrm{\&mu;}}_{o,i}}{K_{\mathrm{ro},i}^{\&prime;}{\mathrm{\&mu;}}_{w,i}}\ln {\left(\frac{X_e{Y}_e{K}_{v,i}}{{3d}_{i}L{K}_{h,i}}\right)}^{0.5}\\ h_{w,i}=h_{\mathrm{wBT},i}+\frac{\underset{j=1}{\mathrm{\&Sigma;}}{B}_{o,i}{q}_{o,i}^j{\mathrm{\&Delta;t}}_j}{X_e{Y}_e{\mathrm{\&phi;}}_i(1-S_{\mathrm{wc},i}-S_{\mathrm{or},i})}\\ h_{\mathrm{wBT},i}=\frac{B_{o,i}{q}_{t,i}{t}_{\mathrm{BT},i}}{X_e{Y}_e{\mathrm{\&phi;}}_i(1-S_{\mathrm{wc},i}-S_{\mathrm{or},i})}\\ t=t_{\mathrm{BT},i}+\underset{j=1}{\mathrm{\&Sigma;}}{\mathrm{\&Delta;t}}_j\\ t_{\mathrm{Bt},i}=\frac{{4\mathrm{\&phi;}}_i{\mathrm{\&Delta;L}}_i{h_i}^2}{B_{o,i}{q}_{t,i}}[1-\cot \frac{{\mathrm{\&pi;d}}_i}{2h_i}\cos \frac{{\mathrm{\&pi;d}}_i}{2h_i}\ln (\sec \frac{{\mathrm{\&pi;d}}_i}{2h_i}+\tan \frac{{\mathrm{\&pi;d}}_i}{2h_i})]\end{array}\right.(i=1,...,m)$

其中，WC_i(t)为第i个控水单元在时间t的含水率；q_w,i为第i个控水单元的产水量；q_o,i为第i个控水单元的产油量；h_w,i为第i个控水单元的平均水侵厚度；K′_rw,i为第i个控水单元的水相端点相对渗透率；K′_ro,i为第i个控水单元的油相端点相对渗透率；μ_w,i为第i个控水单元的水相粘度；μ_o,i为第i个控水单元的油相粘度；K_v,i为第i个控水单元的垂向渗透率；K_h,i为第i个控水单元的水平渗透率；φ_i为第i个控水单元的孔隙度；X_e为泄油宽度；Y_e为泄油长度；L为控水井段长度；S_wc,i为第i个控水单元的束缚水饱和度；S_or,i为第i个控水单元的残余油饱和度；为第i个控水单元在底水突破后第j个时间步的产油量；Δt_j为底水突破后第j个时间步的步长；q_t,i为第i个控水单元在底水突破前的产液量；h_wBT,i为第i个控水单元在底水突破时的平均水侵厚度；t_BT,i为第i个控水单元的底水突破时间；B_w,i为第i个控水单元的水相体积系数；B_o,i为第i个控水单元的油相体积系数；F_c,i为修正系数；h_i为第i个控水单元的拟油柱高度；d_i为第i个控水单元的拟避水高度。

12.根据权利要求11所述的底水油藏水平井变参数射孔控水完井装置，其特征在于，用于：

根据预先设置的h_i和d_i的初始值以及所述公式一确定所述t_BT,i和h_wBT,i；

根据所述公式一确定水平井底水突破后的各时间步的含水率WC_i，在所述WC_i等于一初始含水率WC_0,i时，确定所述h_w,i；

根据所述公式一确定所述公式一的第k个迭代步的拟油柱高度和拟避水高度其中， $h_i^k=h_{0.i}+h_{w,i};d_i^k=d_{0,i}+h_{w,i};$

当一值小于一预先设置的绝对误差阈值时，确定所述h_i的值为确定所述d_i的值为

13.根据权利要求12所述的底水油藏水平井变参数射孔控水完井装置，其特征在于，所述控水单元产液量确定单元，具体用于：

根据所述公式一确定一公式二：

$\left\{\begin{array}{c}{\mathrm{WC}}_1(q_{t,1},t)={\mathrm{WC}}_{\mathrm{limit}}\\ ...\\ {\mathrm{WC}}_i(q_{t,i},t)={\mathrm{WC}}_{\mathrm{linit}}\\ ...\\ {\mathrm{WC}}_m(q_{t,m},t)={\mathrm{WC}}_{\mathrm{limit}}\end{array}\right.(i=1,...,m)$

其中，WC_limit为所述目标极限含水率；WC_i(q_t,i,t)为第i个控水单元的含水率；

根据所述公式二和一预先设置的公式三确定所述各控水单元的产液量q_t,i；其中，所述公式三为：

$\underset{i=1}{\overset{m}{\mathrm{\&Sigma;}}}{q}_{t,i}=Q_t,(i=1,...,m)$

其中Q_t为油井的产液量。

14.根据权利要求13所述的底水油藏水平井变参数射孔控水完井装置，其特征在于，所述射孔总表皮系数确定单元，具体用于：

通过一公式四和公式五确定所述各射孔单元的射孔总表皮系数S_pdc,i(i＝1,2,…,n)；

其中，所述公式四为：

$\left\{\begin{array}{c}\frac{{\mathrm{\&mu;}}_{o,i+1}{B}_{o.i+1}{q}_{t,i+1}}{2\mathrm{\&pi;}{K}_{h,i+1}\mathrm{\&Delta;l}}(\ln \frac{4{\mathrm{\&beta;}}_{i+1}{h}_{i+1}}{\mathrm{\&pi;}{r}_w}+\ln \mathrm{tg}\frac{{\mathrm{\&pi;d}}_{i+1}}{2h_{i+1}}+S_{\mathrm{pdc},i+1})\\ =\frac{{\mathrm{\&mu;}}_{o,i}{B}_{o,i}{q}_{t,i}}{2\mathrm{\&pi;}{K}_{h,i}\mathrm{\&Delta;l}}(\ln \frac{4{\mathrm{\&beta;}}_i{h}_i}{{\mathrm{\&pi;r}}_w}+\ln \mathrm{tg}\frac{\mathrm{\&pi;}{d}_i}{2h_i}+S_{\mathrm{pdc},i})-{\mathrm{\&Delta;P}}_{i,i+1}\\ \mathrm{\&Delta;}{P}_{i,i+1}=-\frac{{\mathrm{\&rho;}}_o\mathrm{f\&Delta;l}}{2D}[\frac{A}{3q_{t,i+1}}(v_{i,i+1}^3-v_{i+1}^3)+\frac{A}{{3q}_{t,i}}(v_i^3-v_{i,i+1}^3)]\\ v_i=(\underset{j=i}{\overset{n}{\mathrm{\&Sigma;}}}{q}_{t,j}-0.5q_{t,i})/A\\ v_{i,i+1}=\underset{j=i+1}{\overset{n}{\mathrm{\&Sigma;}}}{q}_{t,j}/A\end{array}\right.(i=\mathrm{1,2},...,n-1)$

其中，K_v,i为第i个射孔单元的垂向渗透率；K_h,i为第i个射孔单元的水平渗透率；S_pdc,i为第i个射孔单元的射孔总表皮系数；v_i为第i个射孔单元中点处的截面平均流速；v_i,i+1为第i和i+1个射孔单元交点处的截面平均流速；ρ_o为油相密度；f为管壁摩擦系数；r_w为井筒半径；D为套管内径；A为套管的横截面积；μ_o,i为第i个射孔单元的油相粘度；ΔP_i,i+1为第i和i+1个射孔单元之间的压降；

所述公式五为：

Min(S_pdc,i)＝Min(S_k,min)

其中，S_k,min为预设的第k个射孔枪弹组合的最小射孔总表皮系数。

15.根据权利要求14所述的底水油藏水平井变参数射孔控水完井装置，其特征在于，所述射孔枪弹组合确定单元，具体用于：

根据各射孔单元的射孔总表皮系数，确定第i个控水单元的最小射孔总表皮系数 $S_{\mathrm{pdc},\min }^i,(i=1,...,m);$

确定各和S_k,min的差值中的最小非负数对应的射孔枪弹组合为第i个控水单元的射孔枪弹组合。

16.根据权利要求15所述的底水油藏水平井变参数射孔控水完井装置，其特征在于，所述射孔密度确定单元，具体用于：

根据一公式六确定所述各射孔单元的射孔密度：

$\left\{\begin{array}{c}{S}_{\mathrm{pdc},i}=\frac{1}{K_{\mathrm{dr},i}}{S}_{p,i}+S_{d,i}+S_{\mathrm{cz},i},l_{p,i}<r_{d,i}\\ S_{\mathrm{pdc},i}={\hat{S}}_{p,i}+\frac{1}{n_{s,i}{l}_{p,i}}(\frac{1}{K_{\mathrm{cr},i}}-1)\ln \frac{r_{\mathrm{cz},i}}{r_{p,i}},l_{p,i}>r_{d,i}\end{array}\right.$

其中，K_dr,i为第i个射孔单元的污染程度系数；K_cr,i为第i个射孔单元的压实程度系数；S_p,i为第i个射孔单元的射孔几何表皮系数；S_d,i为第i个射孔单元的钻井污染表皮系数；S_cz,i为第i个射孔单元的射孔压实表皮系数；r_d,i为第i个射孔单元的污染带半径；r_cz,i为第i个射孔单元的压实带半径；l_p,i为第i个射孔单元的射孔深度；r_p,i为第i个射孔单元的孔眼半径；n_s,i为第i个射孔单元的射孔密度；为采用有效孔深l_p,i′＝l_p,i-(1-K_dr,i)(r_d,i-r_w)和有效井筒半径r_w′＝r_w+(1-K_dr,i)(r_d,i-r_w)计算的第i个射孔单元的射孔几何表皮系数。