^{1}

^{*}

^{1}

^{1}

^{1}

^{1}

This paper proposed a method of injection-production system adjustment to solve the problem that the water flooding effect was restricted because of the horizontal and vertical contradictions during the development process of fault block reservoirs. Considering the heterogeneity of reservoir, the Buckley-Leverett water flooding theory was applied to establish the relationship between the recovery and cumulative water injection. In order to achieve the goal of vertically balanced recovery of each section, the calculation method of vertical sectional injection allocation was proposed. The planar triangular seepage unit was assumed and sweep coefficients of different oil-water distribution patterns were characterized using multi-flow tube method. In order to balance and maximize the plane sweep coefficient, the calculation method of plane production system optimization was obtained. Then the injection-production system stereoscopic adjustment method based on equilibrium displacement was proposed with vertical sectional injection allocation and plane production system optimization. This method was applied to injection and production adjustment of BZ oilfield in southern Bohai. The effect of water control and oil increase was obvious. This method can greatly improve the effect of water flooding of offshore fault block reservoirs with the adjustment of injection-production system.

The complex fault-block reservoirs in the southern Bohai Sea are characterized by rapid changes in planar physical properties, multiple vertical layers, strong vertical and planar heterogeneity. Irregular well pattern with directional well is used to develop this type of reservoirs. Due to the influence of reservoir properties, fluid properties and well pattern, the vertical and horizontal contradictions always lead to the disequilibrium displacement in the water flooding development process of oilfields. And the disequilibrium displacement leads to the rapid increase of water cut and decline in production. Equilibrium displacement to development of reservoir can achieve maximum oil recovery and economic benefits. Separated layer injection adjustment and plane production system optimization are important ways to improve the effect of water flooding in reservoir development.

Many studies about separated layer injection adjustment and plane production system optimization have been carried on. For the separated layer injection adjustment, methods based on static parameters and field experiences are commonly used. WU Jiawen (2006) proposed an optimization method for layered water injection project based on the distribution of remaining oil [

J. I. Bingyu (2008) studied the method of oil production calculation for areal well pattern of low-permeability reservoir with non-Darcy seepage flow [

On the basis of previous studies, using non-piston water flooding theory, this paper studies the calculation method of injection allocation for realizing balanced recovery of each section. Irregular well pattern is divided into different triangular injection-production seepage units in plane, and the optimization method of injection and production pressure difference is obtained by MTM to realize production system optimization. According to the calculation results of water injection allocation volume of each layer, separate injection of the injection wells are adjusted, and according to the calculation results of the plane injection production pressure difference, the oil nozzles or pump frequency of the production wells are adjusted, then injection and production system adjustment method based on balanced displacement is obtained by combining the two methods. The method has been applied in oil field, and the effect of oil stabilization and water control is achieved.

The purpose of sectional water injection is to alleviate the contradiction between layers and reduce the interlayer interference caused by the difference of reservoir and fluid physical properties. By adjusting the injection volume of each section, the oil recovery of each section is consistent at a certain time, then the vertical balanced displacement is realized eventually. According to the degree of recovery at different time, the cumulative injection can be calculated, and the sectional injection is obtained under constant speed injection.

If injection water doesn’t breakthrough in the jth section, according to Bukley-Leverett’s non-piston water flooding theory [

W j = ∫ 0 T Q j d t = x f j ϕ A / f ′ w ( S wf j ) (1)

where W_{j} is the accumulated water injection in the jth section, T is the cumulative water injection time, Q_{j} is the water injection volume in section, t is the injection time, X_{fj} is the oil-water front position in the jth section, ϕ is reservoir porosity, A is the seepage section area, S_{wfj} is the water saturation of oil-water front in section, f ′ w ( S wf j ) is the derivative of the water cut corresponding to the water saturation of oil-water front.

Under the condition that the amount of injection remains constant, the average injection volume can be obtained according to the following formula.

Q ¯ j = ϕ A [ L / f ′ w ( S w 2 j ) − x f j / f ′ w ( S wf j ) ] / Δ t (2)

where Q ¯ j is the allocation of water injection in the jth section, Δt is allocation of water injection time, L is the distance of injection to prodution, S_{w2j} is the water saturation at the outlet of the jth section after water breakthrough, f ′ w ( S w 2 j ) is the water content derivative corresponding to water saturation at the outlet after water breakthrough. The allocation water injection rate cannot be solved directly, but can be obtained iteratively by setting the predicted outlet saturation.

Before water breakthrough, the recovery of reservoir in the jth section can be obtained by cumulative injection water:

R j = W j / ( ϕ A L ) (3)

After water breakthrough, the recovery of reservoir in the jth section can be obtained by average water saturation:

R j = ( S ¯ w j − S wi ) / ( 1 − S wi ) = { S w 2 j + [ 1 − f w ( S w 2 j ) ] / f ′ w ( S w 2 j ) − S wi } / ( 1 − S wi ) (4)

where R_{j} is the oil recovery factor, S_{wi} is the irreducible water saturation, f_{w}(S_{w2j}) is the outlet water content in the jth section.

If injection water has breakthrough in the jth section, the cumulative injection volume required for the water saturation from S_{wj} to S_{w2j} at the outlet of the jth section can be expressed as

∫ T T + Δ t Q j d t = L ϕ A [ 1 / f ′ w ( S w 2 j ) − 1 / f ′ w ( S w j ) ] (5)

where S_{wj} is the outlet water saturation of the jth section before allocation, f ′ w ( S w j ) is the water content derivative corresponding to water saturation at the outlet after water breakthrough.

Assuming that the injection allocation of each section remains constant in Δt time, the injection allocation of the jth section can be obtained by

Q ¯ j = L ϕ A [ 1 / f ′ w ( S w 2 j ) − 1 / f ′ w ( S w j ) ] / Δ t (6)

When this method is used to segmental injection allocation, the recovery factor of each section can be obtained according to the accumulated water volume of each layer. By setting allocation time and target recovery factor, the accumulated water volume of each segment is calculated, and then the injection allocation amount of each section can be obtained.

The injection-production well pattern is divided into different triangular seepage units (1/2 of the maximum area can be affected by water injection wells). The two vertices of the triangle seepage unit are water injection well and corresponding oil production well. The triangle seepage unit is assumed to be homogeneous and equal thickness reservoir. The oil and water well are composed of flow tubes. The seepage rule in the flow tubes conforms to the oil-water two-phase piston displacement. The compressibility of formation rocks and fluids, the capillary force and gravity are neglected, and the injection-production pressure difference is constant. As shown in

α / β = α 0 / β 0 (7)

where α is the angle between streamline AO and injection-production connection line, β is the angle between streamline BO and injection-production connection line, α_{0} is the angle between sweep maximum area triangle side AC and injection-production connection line, β_{0} is the angle between sweep maximum area triangle side BC and injection-production connection line.

In the case of oil and water two-phase flow, the length of the water phase displacement of any angle at any time t is [

ξ = C ( p inj − p wf ) μ o K o − μ w K w { μ o K o L ( sin α + sin β ) sin ( α + β ) p inj − p wf − − 2 ( μ o K o − μ w K w ) p inj − p wf t + [ μ o K o L ( sin α + sin β ) sin ( α + β ) p inj − p wf ] 2 } (8)

where C is a constant(C = 0.0864), p_{inj} is the bottom hole pressure of injection well, p_{wf} is the bottom hole flow pressure of production well, μ_{o} is the viscosity of reservoir crude oil, μ_{w} is formation water viscosity, K_{o} is the oil phase permeability, K_{w} is the water phase permeability.

According to the distribution of oil and water at different times, the whole displacement process can be divided into four periods: when t is less than T ′ 1 (α = 0, the time when the oil-water front reaches the turning point O of AB), the oil-water fronts in flow tubes don’t reach the turning point of the flow tube. When t is greater than T ′ 1 and less than T ′ 2 (α = 0, the time when the oil-water front reaches point B), part of the oil-water fronts in flow tubes have reached the turning point, but not reached the production well. When t is greater than T ′ 2 and less than T ′ 3 (α = α_{0}, the time of the oil-water front reaching C point), most of the oil-water fronts in flow tubes have crossed the turning point, and some of them reach the production well. When t is greater than T ′ 3 and less than T ′ 4 (α = α_{0}, the time of the oil-water front reaching B point), the oil-water fronts in flow tubes have crossed the turning point, and reached the production well.

The area sweep coefficient of triangular seepage unit is the ratio of sweep area to unit area, and the water content is the ratio of flow rate in flow tubes that the front of water drives to the well to all flow rate in flow tubes.

f w = ∑ j = 1 n Q j / ∑ i = 1 m Q i (9)

E a = ∑ i = 1 n S i / S (10)

where f_{w} is the water cut of production well, E_{a} is the area sweep coefficient, m is the total number of flow pipes, n is the number of flow pipes that the front of the water reaches the production well, Q_{i}, Q_{j} are the flow rates of the ith and jth flow pipe, S_{i} is the area of the ith flow pipe, S is the area of ΔABC.

The area sweep coefficient and water cut of production well are calculated by judging the distribution in each flow pipe at different time. Setting the target area sweep coefficient of the adjustment time, the injection-production pressure difference of the corresponding seepage unit is inversely calculated. Through the adjustment of production system optimization, the eventual balanced displacement can be achieved.

Based on the above theory, a method of injection-production system adjustment for offshore fault-block reservoirs is established. According to cumulative water injection and water content of each section vertically, the recovery degree and water saturation at the outlet can be obtained by Equation (1) or Equation (5). The target water saturation at the outlet end can be calculated by Equation (4) after setting the target recovery degree of each section. And the average injection allocation of each section in the target time is calculated according to Equation (2) and Equation (6). On the plane, injection-production triangle seepage units are divided. And fluid properties and well pattern parameters of each unit are obtained. According to the current water content of production well in the unit, the area sweep coefficient and the oil-water front distribution are derived. The area sweep coefficient of each seepage unit is compared, and the seepage unit with smaller plane area sweep coefficient needs to increasing injection-production pressure difference to improve water flooding effect. By setting the sweep coefficient of the target area of each seepage unit, the injection-production pressure difference is calculated according to Equation (8). The production system of production wells is adjusted according to injection-production pressure difference. Equilibrium displacement can be achieved through vertical and horizontal injection-production adjustment, thus the development effect of water flooding can be improved. The calculation procedure is shown in

The injection-production adjustment of S1 block in BZ oilfield has been conducted by using this method, as shown in _{1}m^{l} I-1 layer and I-2 layer. The development mode of commingled production by directional production wells and sectional water injection by directional water injection wells. The injection and production system adjustment of well I1-P1 is shown in this paper. The well spacing from I1 well to P1 well and P2 well is respectively 590 m and 420 m, and the well spacing between P1 well and P2 well is 745 m. The reservoir and fluid parameters are shown in

Layer | Parameters | ||||
---|---|---|---|---|---|

Thickness (m) | Porosity (%) | Permeability (mD) | Viscosity of formation water (mPa∙s) | Viscosity of crude oil (mPa∙s) | |

I-1 | 6.5 | 30.0 | 1051 | 0.1 | 12.9 |

I-2 | 8.1 | 32.0 | 1983 | 0.1 | 12.5 |

The test results of production profile of well P1 before allocating shows the water cut in I-1 layer is 45.0% and the water content in I-2 layer is 60.8%. Calculated by Equation (5), the recovery factor of layer I-1 is 13.8%, the recovery factor in layer I-2 is 17.7%. The calculation result shows that due to good reservoir physical properties of I-2 layer, the injected water front advances rapidly along this layer. The recovery of layer I-1 is lower than that of layer I-2. Water injection allocation is needed to tap the remaining oil potential. The allocation time is set to 300 days, and the target recovery of the two layers is 20.0%. The average injection allocation of layer I-1 is 96 m^{3}/d and that of layer I-2 is 68 m^{3}/d, which is calculated by Equation (6).

The well group is divided into two seepage units as shown in _{1} is 35 degrees, and the angle θ_{2} is 53 degrees. The area sweep coefficient is calculated for each seepage unit respectively by applying this method. The area sweep coefficient of seepage unit (1) is 0.52. The area sweep coefficient of seepage unit (2) is 0.78. In order to achieve equilibrium displacement, the sweep area of seepage unit (1) needs to be increased. The target time is 300 days and the target area sweep coefficient is 0.80. The injection-production pressure difference in well I1-P1 direction is optimized. By increasing the production frequency of pump in well P1, the injection-production pressure difference is increased from 10.0 MPa to 11.7 MPa.

After adjustment, the effect of water control and oil increase is obvious, as shown in ^{3}, and the water cut remains stable. This method has been applied to adjust injection-production system of other well groups in this block, and the effect is remarkable. The daily oil increase is 110 m^{3}, and the annual production decline rate is only 1.5%, which improves the development effect of water flooding.

It should be figured out that absolute equilibrium displacement cannot be achieved due to the heterogeneity of actual reservoirs. In addition, the adjustment of injection-production system is a dynamic process. A relatively balanced displacement needs to be adjusted in time according to the development dynamics and combined with the corresponding technological measures.

1) Aiming at the prominent contradiction between layers in multi-layer development of offshore fault-block reservoirs, based on Buckely-Leverett non-piston water flooding theory, a sectional injection allocation method is established. By setting the unified target recovery of each section, the vertical balanced displacement can be achieved.

2) In order to maximize the plane sweep area of water flooding, injection-production well group is divided into triangle seepage units. By using multi-flow pipe method, the adjustment method of production system optimization for offshore fault-block reservoirs with irregular pattern development is established to realize horizontal balanced displacement.

3) Combining the vertical sectional injection allocation method and planar production system optimization method, the injection-production system adjustment method of fault block reservoir is established. This method is applied to guide the injection-production system adjustment in block S1 of BZ oilfield. The practical application results show that this method has obvious effect on water control and oil production increase. And the annual production decline rate of the S1 block is only 1.5%. This method has achieved good application effect, which can conduct the injection-production system adjustment of offshore fault block reservoirs.

4) It should be pointed out that the method of sectional water injection calculation proposed in this paper does not consider the influence of plane well pattern. Moreover, the hypothesis of triangular seepage unit division with planar production system optimization is limited to the reservoir with strong heterogeneity. Further study is needed to improve the above limitations.

The authors declare no conflicts of interest regarding the publication of this paper.

Yang, M., Chen, C.L., Wang, Y., Wu, X.H. and Ma, D. (2020) Study on Injection-Production System Adjustment of Fault Block Reservoir Based on Equilibrium Displacement. Open Journal of Yangtze Gas and Oil, 5, 35-45. https://doi.org/10.4236/ojogas.2020.52004

W_{j} Accumulated water injection in the jth section, m^{3}

T Cumulative water injection time, d

Q_{j} Water injection volume in section, m^{3}

t Injection time, d

X_{fj} Oil-water front position in the jth section, m

ϕ Reservoir porosity, dimensionless

A Seepage section area, m^{2}

S_{wfj} Water saturation of oil-water front in section, m

f ′ w ( S wf j ) Derivative of the water cut corresponding to the water saturation of oil-water front, dimensionless

Q ¯ j Allocation of water injection in the jth section, m^{3}

Δt Allocation of water injection time, d

L Distance of injection to prodution, m

S_{w2j} Water saturation at the outlet of the jth section after water breakthrough, dimensionless

f ′ w ( S w 2 j ) Water content derivative corresponding to water saturation at the outlet after water breakthrough, dimensionless

R Oil recovery factor, dimensionless

S_{wi} Irreducible water saturation, dimensionless

f_{w}(S_{w2j}) Outlet water content in the jth section, dimensionless

S_{wj} Outlet water saturation of the jth section before allocation, dimensionless

f ′ w ( S w j ) Water content derivative corresponding to water saturation at the outlet after water breakthrough, dimensionless

p_{inj} Bottom hole pressure of injection well, MPa

p_{wf} Bottom hole flow pressure of production well, MPa

μ_{o} Viscosity of reservoir crude oil, mPa∙s

μ_{w} Viscosity of formation water, mPa∙s

K_{o} Oil phase permeability, dimensionless

K_{w} Water phase permeability, dimensionless

f_{w} Water cut of production well, dimensionless

m Total number of flow pipes, dimensionless

n Number of flow pipes that the front of the water reaches the production well, dimensionless

Q_{i} Flow rates of the ith flow pipe, dimensionless

E_{a} Area sweep coefficient, dimensionless

S Area of triangle seepage unit, m^{2}

S_{i} Area of the ith flow pipe, m^{2}

j Serial number of layer

i Serial number of flow tube

w water phase

o oil phase