# Analytical Model of Waterflood Sweep Efficiency in Vertical Heterogeneous Reservoirs under Constant Pressure

Authors: Lisha Zhao; Li Li; Zhongbao Wu; Chenshuo Zhang
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

An analytical model has been developed for quantitative evaluation of vertical sweep efficiency based on heterogeneous multilayer reservoirs. By applying the Buckley-Leverett displacement mechanism, a theoretical relationship is deduced to describe dynamic changes of the front of water injection, water saturation of producing well, and swept volume during waterflooding under the condition of constant pressure, which substitutes for the condition of constant rate in the traditional way. Then, this method of calculating sweep efficiency is applied from single layer to multilayers, which can be used to accurately calculate the sweep efficiency of heterogeneous reservoirs and evaluate the degree of waterflooding in multilayer reservoirs. In the case study, the water frontal position, water cut, volumetric sweep efficiency, and oil recovery are compared between commingled injection and zonal injection by applying the derived equations. The results are verified by numerical simulators, respectively. It is shown that zonal injection works better than commingled injection in respect of sweep efficiency and oil recovery and has a longer period of water free production.


---

## Body

## 1. Introduction

Heterogeneity is a common problem encountered in oil-bearing formations. The most significant property that affects waterflooding performance is the matrix permeability, and its variation in the vertical direction causes displacing fluid to advance faster in zones of higher permeability and results in earlier breakthrough in such layers. This phenomenon negatively affects volumetric sweep efficiency and leads to low ultimate oil recovery [1–5]. Therefore, it is considerably important to establish a quantitative characterization method for calculating sweep efficiency; by this method, we can determine the criteria of dividing and reorganizing layer series and relieve interlayer interference in stratified waterflooding reservoir development.Vertical sweep efficiency is the fraction of vertical cross-sectional area of the reservoir between injection and production wells that is swept by water at a given time. It can be used to estimate and predict the unswept fraction of the reservoir by water injection and the additional oil recovery potential.Buckley and Leverett [6] presented the fractional flow equation based on mass conservation equation. Welge [7] proposed that waterflood sweep efficiency could be obtained by water cut and water saturation curves. To multilayer reservoirs, Stiles [8] presented the first model for waterflood calculations in stratified reservoirs and assumed velocities in different layers to be proportional to their absolute permeabilities. Later models for noncommunicating layers without crossflow [9–11] and models for communicating layers with complete crossflow [12, 13] were carried out. El-Khatib [14, 15] investigated the effect of crossflow on the performance of stratified reservoirs and presented a closed form analytical solution for communicating stratified systems with log-normal permeability distributions. Zhou et al. [16, 17] developed a model of linear nonpiston waterflood and calculated the sweep efficiency by production compartment. El-Khatib [18, 19] considered the gravity effect of the dip angle to sweep efficiency and Snyder and Ramey [20] improved on previous work while changing the property of layers in a logical and consistent manner. Prince [21] investigated the effect of mobility ratio and the Dykstra-Parsons permeability variation coefficient (VDP) on the performance.Firstly, many of the models mentioned above assumed piston-like displacement when predicting waterflooding performance in stratified reservoirs, which was not consistent with the frontal advance theory. Secondly, the abovementioned study of flooding process was usually based on the condition of constant rate, which provided a constant volume of water injected into each layer by production compartment according to each layer’s property. However, for multilayered reservoir development, the production pressure between injection and production well is usually stable, while the rate of liquid produced in each layer is changeable. So, in this model, the reservoir is divided into a number of layers. And each layer is considered as non-piston-like displacement on the condition of constant pressure. A mathematical model is developed for calculating waterflood sweep efficiency by extending the Buckley-Leverett displacement mechanism.


## 2. Mathematical Model Building and Assumptions

Figure1 is a schematic representation of the stratified system. The production pressure between injection and production well is equal to a constant value: ΔP=Pwf-Pinj, where Pinj is the bottom pressure of injection well, MPa; Pwf is the bottom pressure of producing well, MPa; ΔP is the pressure drop from the injection well to producing well, MPa.Figure 1
Schematic representation of stratified system.The following assumptions are made:(1)
The system is divided into a number of homogeneous layers; each has a uniform thickness and constant permeability.(2)
The system is linear and horizontal and of constant thickness.(3)
The flow is isothermal and incompressible and obeys Darcy’s law.(4)
The displacement is non-piston-like with a two-phase region where both oil and water exist.(5)
Capillary and gravity forces are negligible.(6)
The system is noncommunicating with no crossflow allowed between adjacent layers.(7)
The relative permeability characteristics are the same for all layers.(8)
The initial fluid saturation is uniform at the irreducible water saturation.(9)
The porosity is assumed to be constant in all layers.


## 3. Calculation for Single Layer

The performance before displacement front reaches the outlet face is quite identical in each layer. After water breakthrough, injection continues, assuming that the front of displacement remains going forward, and the outlet saturation rises fromSwf to 1-Sor, where Swf is the water saturation of displacement front, f, and Sor is the residual oil saturation, f. Figure 2 illustrates the three stages during waterflooding. This process is split into three subperiods, and piecewise functions of sweep efficiency are deduced with the breakthrough time as the endpoint.Figure 2
Different stages during water/oil displacement.
(a)
Before water breakthrough (b)
At the time of breakthrough (c)
After water breakthroughBased on the material balance equation and Buckley-Leverett’s frontal advance theory, the total volume of water injected in a given time equals(1)∫0tQdt=∫x0xfϕwhSwx,t-Swcdx=ϕwhxf-x0fw′Swf,where Q is the flow rate of oil, m3/s; w is the width of each layer, m; h is the height of each layer, m; ϕ is the porosity of matrix, f; xf is the displacement front position, m; x0 is the initial position of displacement front, m; fw(Sw) is the fractional flow of water, f.The calculation of sweep efficiency in a single layer is given by(2)η=∫0xfSw-Swc/1-Sor-SwchdxhLe=S¯w-Swc1-Sor-Swc·xfLe,where Le is the length of model, m. Before water breakthrough, average saturation S¯w is determined by fractional flow curve (fw-Sw) [6], which is a constant value. So, the flooding swept volume is only determined by the position of the water frontal advance; after breakthrough, (2) deforms as follows:(3)η=S¯w-Swc1-Sor-Swc,where S¯w is no longer a constant and should be determined by the saturation of the producing well Swe. Consequently, it is necessary to calculate the waterflood front location before breakthrough and the saturation of the outlet face Swe after breakthrough. Firstly, we analyze the waterflooding performance in a single layer by studying the pressure relationship based on the constant pressure condition.


### 3.1. Pressure in Oil and Water Region

By Applying Darcy’s law for a single layer, the total velocity at a certain positionxa between injection well and the front location is given by(4)v=vw+vo=-KKrwμw+Kroμodpdx,where K is the absolute permeability, 10-3 µm2; Krw is the water relative permeability, f; Kro is the oil relative permeability, f; μw is the viscosity of water, mPa·s; μo is the viscosity of oil, mPa·s.Based on the Buckley-Leverett frontal displacement theory, the water cut is expressed as(5)fwSw=QwQw+Qo=11+Kroμw/Krwμo.Substituting (5) into (4) and integrating, we obtain the pressure drop between the positions of x=xa:(6)ΔPa=Pa-Pinj=vμwK∫0xafwSwKrwdx.The location of the injection well is treated as the position ofx=0, and the frontal displacement function can be written as(7)x=fw′Swϕwh∫0tQdt.Equation (7) can be deformed as(8)x=xffw′Swfw′Swf.From the position of x=0 to x=xa, the water saturation changes from 1-Sor to Swa. Substituting (8) into (6), we get(9)ΔPa=vμwK·xffw′Swf·FSwa,where F(Swa) is given by(10)FSwa=∫0fw′SwafwSwKrwdfw′Sw.


### 3.2. Time of Water Breakthrough

The pressure drop from the injection well to the displacement front location is defined asΔPf. Since the average saturation is a constant before water breakthrough, setting F(Swf) to value a and substituting value a into (9) yield(11)ΔPf=Pf-Pinj=vμwK·xffw′Swf·a.Substituting (7) into (11), we get(12)ΔPf=vμwaWtϕwhK,where W(t) represents the accumulated volume of water injected from the time of beginning, and it is equal to the total liquid volume produced from the oil well:(13)Wt=∫0tQdt=∫0twhvdt.In the pure oil region, there is only one phase, and the pressure drop of this region is(14)ΔPo=ΔP-ΔPf=μovKLe-xf.Substituting (11) and (12) into (14), we get(15)ΔPf=aΔPWtϕwhμrLe-μrWtfw′Swf+aWt,where μr is the ratio of oil viscosity and water viscosity, f. Taking the derivative of (13) and substituting it into (12), we obtain(16)ΔPf=μwaWtϕwh2KdWtdt.Combining (15) with (16) and integrating and rearranging the equation, the accumulated volume of water from t=0 to t=tn is obtained:(17)Wtn=μoϕwhLe2+2μwa-μofw′Swfϕwh2KΔPtn-μoϕwhLeμwa-μofw′Swf.Substituting (7) into (17), the position of displacement front xf is given by(18)xf=fw′Swfϕ·μoϕLe2+2μwa-μofw′SwfϕKΔPt-μoϕLeμwa-μofw′Swf.When the front reaches the outlet face,xf equals Le; substituting it into (18) and rearranging the equation, the time of water breakthrough can be obtained:(19)T=μwϕLea/fw′Swf2-μoϕLe22μwa-μofw′SwfϕKΔP.


### 3.3. Saturation after Water Breakthrough

After the time of water breakthrough, the producing well starts to produce water; according to (9), the pressure drop between the two wells is(20)ΔP=vμwK·xffw′Swf·FSwe,where Swe is the saturation of the producing well, f. After breakthrough, water injection continues, assuming that the front of displacement remains going forward, and the outlet saturation rises from Swf to Swe.Substituting (7) into (13) and taking the derivative, we get(21)dWt=whvdt=ϕwhLed1fw′Swe.Solving (20) and (21) simultaneously, we get(22)dt=μwLe2ϕKΔP·FSwefw′Swfd1fw′Swe.Making the integral of (22), (23)t-T=μwLe2ϕKΔPG1fw′Swe,where G(Swe) is given by(24)G1fw′Swe=∫1/fw′Swf1/fw′SwFSwefw′Swfd1fw′Swe.Getting (8) deformed and substituting it into (24),(25)t-T=μwLe2ϕKΔPGxLe·1fw′Sw.Equation (25) shows the water saturation of any position within the two wells versus injecting time after water breakthrough. Equation (23) is a critical condition of (25) since it refers to the position of the oil well.


## 3.1. Pressure in Oil and Water Region

By Applying Darcy’s law for a single layer, the total velocity at a certain positionxa between injection well and the front location is given by(4)v=vw+vo=-KKrwμw+Kroμodpdx,where K is the absolute permeability, 10-3 µm2; Krw is the water relative permeability, f; Kro is the oil relative permeability, f; μw is the viscosity of water, mPa·s; μo is the viscosity of oil, mPa·s.Based on the Buckley-Leverett frontal displacement theory, the water cut is expressed as(5)fwSw=QwQw+Qo=11+Kroμw/Krwμo.Substituting (5) into (4) and integrating, we obtain the pressure drop between the positions of x=xa:(6)ΔPa=Pa-Pinj=vμwK∫0xafwSwKrwdx.The location of the injection well is treated as the position ofx=0, and the frontal displacement function can be written as(7)x=fw′Swϕwh∫0tQdt.Equation (7) can be deformed as(8)x=xffw′Swfw′Swf.From the position of x=0 to x=xa, the water saturation changes from 1-Sor to Swa. Substituting (8) into (6), we get(9)ΔPa=vμwK·xffw′Swf·FSwa,where F(Swa) is given by(10)FSwa=∫0fw′SwafwSwKrwdfw′Sw.


## 3.2. Time of Water Breakthrough

The pressure drop from the injection well to the displacement front location is defined asΔPf. Since the average saturation is a constant before water breakthrough, setting F(Swf) to value a and substituting value a into (9) yield(11)ΔPf=Pf-Pinj=vμwK·xffw′Swf·a.Substituting (7) into (11), we get(12)ΔPf=vμwaWtϕwhK,where W(t) represents the accumulated volume of water injected from the time of beginning, and it is equal to the total liquid volume produced from the oil well:(13)Wt=∫0tQdt=∫0twhvdt.In the pure oil region, there is only one phase, and the pressure drop of this region is(14)ΔPo=ΔP-ΔPf=μovKLe-xf.Substituting (11) and (12) into (14), we get(15)ΔPf=aΔPWtϕwhμrLe-μrWtfw′Swf+aWt,where μr is the ratio of oil viscosity and water viscosity, f. Taking the derivative of (13) and substituting it into (12), we obtain(16)ΔPf=μwaWtϕwh2KdWtdt.Combining (15) with (16) and integrating and rearranging the equation, the accumulated volume of water from t=0 to t=tn is obtained:(17)Wtn=μoϕwhLe2+2μwa-μofw′Swfϕwh2KΔPtn-μoϕwhLeμwa-μofw′Swf.Substituting (7) into (17), the position of displacement front xf is given by(18)xf=fw′Swfϕ·μoϕLe2+2μwa-μofw′SwfϕKΔPt-μoϕLeμwa-μofw′Swf.When the front reaches the outlet face,xf equals Le; substituting it into (18) and rearranging the equation, the time of water breakthrough can be obtained:(19)T=μwϕLea/fw′Swf2-μoϕLe22μwa-μofw′SwfϕKΔP.


## 3.3. Saturation after Water Breakthrough

After the time of water breakthrough, the producing well starts to produce water; according to (9), the pressure drop between the two wells is(20)ΔP=vμwK·xffw′Swf·FSwe,where Swe is the saturation of the producing well, f. After breakthrough, water injection continues, assuming that the front of displacement remains going forward, and the outlet saturation rises from Swf to Swe.Substituting (7) into (13) and taking the derivative, we get(21)dWt=whvdt=ϕwhLed1fw′Swe.Solving (20) and (21) simultaneously, we get(22)dt=μwLe2ϕKΔP·FSwefw′Swfd1fw′Swe.Making the integral of (22), (23)t-T=μwLe2ϕKΔPG1fw′Swe,where G(Swe) is given by(24)G1fw′Swe=∫1/fw′Swf1/fw′SwFSwefw′Swfd1fw′Swe.Getting (8) deformed and substituting it into (24),(25)t-T=μwLe2ϕKΔPGxLe·1fw′Sw.Equation (25) shows the water saturation of any position within the two wells versus injecting time after water breakthrough. Equation (23) is a critical condition of (25) since it refers to the position of the oil well.


## 4. Calculation for Multiple Layers

### 4.1. Water Breakthrough in the First (Most Permeable) Layeri

(1)
The layeri The time of water breakthroughTi in layer i can be calculated according to (19), and the sweep volume of layer i is calculated by(26)Ai=S¯w-Swc1-Sor-Swc·whiLe, whereS¯w is as in Welge’s equation [7]: (27)S¯w=Swc+1fw′Swf. (2)
The other layers, for example, layerj SubstitutingTi into (18), we can get the front position xfj of layer j. The sweep volume of layerj is(28)Aj=S¯w-Swc1-Sor-Swc·whjxfj. (3)
Total sweep efficiency(29)η=Ai+∑j=1n-1AjhLe.


### 4.2. Water Breakthrough in the Medium (including the Least) Permeable Layeri

(1)
The layeri It is the same as in Section4.1(1). (2)
Layers with breakthrough before layeri, for example, layer j The time of breakthroughTj in layer j is calculated by (19). The saturation and water cut of layers that got breakthrough beforeTi can be calculated by(30)Ti-Tj=μwLe2ϕKjΔPG1fw′Swe. After water breakthrough, Welge’s equation reforms as follows:(31)S¯w=Swe+1-fwSwefw′Swe. Substituting (31) into (26), the sweep volume of layer j can be obtained. (3)
Layers without breakthrough, for example, layerk This is the same as in Section4.1(2). (4)
Total sweep efficiency(32)η=Ai+∑j=1mAj+∑k=1n-m-1AkhLe.


## 4.1. Water Breakthrough in the First (Most Permeable) Layeri

(1)
The layeri The time of water breakthroughTi in layer i can be calculated according to (19), and the sweep volume of layer i is calculated by(26)Ai=S¯w-Swc1-Sor-Swc·whiLe, whereS¯w is as in Welge’s equation [7]: (27)S¯w=Swc+1fw′Swf. (2)
The other layers, for example, layerj SubstitutingTi into (18), we can get the front position xfj of layer j. The sweep volume of layerj is(28)Aj=S¯w-Swc1-Sor-Swc·whjxfj. (3)
Total sweep efficiency(29)η=Ai+∑j=1n-1AjhLe.


## 4.2. Water Breakthrough in the Medium (including the Least) Permeable Layeri

(1)
The layeri It is the same as in Section4.1(1). (2)
Layers with breakthrough before layeri, for example, layer j The time of breakthroughTj in layer j is calculated by (19). The saturation and water cut of layers that got breakthrough beforeTi can be calculated by(30)Ti-Tj=μwLe2ϕKjΔPG1fw′Swe. After water breakthrough, Welge’s equation reforms as follows:(31)S¯w=Swe+1-fwSwefw′Swe. Substituting (31) into (26), the sweep volume of layer j can be obtained. (3)
Layers without breakthrough, for example, layerk This is the same as in Section4.1(2). (4)
Total sweep efficiency(32)η=Ai+∑j=1mAj+∑k=1n-m-1AkhLe.


## 5. Case Study and Verification

Take one injection-production unit of Yaerxia oilfield in Yumen for the case study; the parameters of reservoir and fluid properties are listed in Table1, and the relative permeability curves are plotted in Figure 3. By analyzing the data of logging, five layers are classified of the main layer k1g22-2.Table 1
The parameters of reservoir and fluid properties.
Variables Value Length of the model,Le (m) 200 Width of the model,w (m) 100 Thickness of each layer,h (m) 1 Pressure drop of the unit,ΔP (MPa) 10 Oil viscosity in the reservoir,μo (mPa⋅s) 5.6 Water viscosity in the reservoir,μw (mPa⋅s) 0.6 Average saturation before breakthrough,S¯w (f) 0.47 Saturation of front before breakthrough,Swf (f) 0.33 F ( S w f ) before breakthrough, a (f) 61.5 Residual oil saturation,Sor (f) 0.30 Irreducible water saturation,Swc (f) 0.32 f ′ ( S w f ) before breakthrough, (f) 6.65 Porosity of the model,ϕ (f) 0.1 Permeability of layer 1,k1 (10-3 μm2) 16.6 Permeability of layer 2,k2 (10-3 μm2) 13 Permeability of layer 3,k3 (10-3 μm2) 0.8 Permeability of layer 4,k4 (10-3 μm2) 2.4 Permeability of layer 5,k5 (10-3 μm2) 1.2Figure 3
Relative permeability curves of the case.There are two kinds of injection programs to be compared with: commingled injection and zonal injection. In commingled injection, all layers are injected and produced under the same pressure condition, while, in zonal injection, the pressure drop and the rate of water flow are different and can be adjusted according to the waterflooding performance. In this case, considering their property differences, five layers are divided into two groups in zonal injection. Put a packer between layer 2 and layer 3 and make layer 1 and layer 2 group I and layers 3~5 group II. In order to achieve a balanced frontal advance, the pressure drop of group I is reduced by half and that of group II is tripled. Use the method above to calculate water frontal position and sweep efficiency of these two injection programs, and the results are verified by a numerical simulator.


### 5.1. Water Frontal Position

Use the equations for single layer to calculate the time of water breakthrough of the five layers by (19) for both commingled injection and zonal injection. T1 is the breakthrough time of the most permeable layer (layer 1). T1 is equal to 214.1 days in the commingled injection program while it is 428.2 days in the zonal injection program because of the pressure variation. The water frontal position of all layers on T1 is calculated by (18). The corresponding sweep volume of each layer is calculated by (26) and (28). The results of equation calculation are listed in Table 2 and results of numerical simulation are listed in Table 3. We plot the water frontal position on T1 for the two injection programs (Figures 4 and 5). The calculation results are compared with numerical simulation results, where the former are plotted on the top and the latter are plotted right below. The comparisons show that the simulation results match the calculation results well, which verify the equations of single layer.Table 2
The results of analytical model in commingled injection and zonal injection.
Layers Commingled injection Zonal injection Breakthrough time,T (day) Frontal position onT1, xf (m) Sweep volume onT1, A (m3) Breakthrough time,T (day) Frontal position onT1, xf (m) Sweep volume onT1, A (m3) Layer 1 214.1 200.0 7985.7 428.1 200.0 7985.7 Layer 2 273.3 157.3 6280.6 546.7 157.3 6280.6 Layer 3 4441.8 9.8 392.3 1480.6 58.7 2342.1 Layer 4 1480.6 29.4 1174.5 493.5 173.9 6945.4 Layer 5 2961.2 14.7 588.1 987.1 87.7 3502.8Table 3
The results of numerical model in commingled injection and zonal injection.
Layers Commingled injection Zonal injection Breakthrough time,T (day) Frontal position onT1, xf (m) Sweep volume onT1, A (m3) Breakthrough time,T (day) Frontal position onT1, xf (m) Sweep volume onT1, A (m3) Layer 1 182.0 200.0 7985.7 386.0 200.0 7985.7 Layer 2 231.0 158.0 6308.7 493.0 156.0 6228.9 Layer 3 3567.0 12.0 479.1 1238.0 56.0 2236.0 Layer 4 1214.0 31.0 1237.7 441.0 172.0 6867.7 Layer 5 2502.0 16.0 638.9 874.0 86.0 3433.9Figure 4
Results comparison of water frontal position of commingled injection onT1.
(a)
Results of equation calculation (b)
Results of numerical simulationFigure 5
Results comparison of water frontal position of zonal injection onT1.
(a)
Results of equation calculation (b)
Results of numerical simulation


### 5.2. Water Cut and Sweep Efficiency

Use the equations for multiple layers to calculate water cut of the model in both of the injection programs by (5) and calculate sweep efficiency by (29) and (32), and oil recovery can be obtained. The curves versus time are shown in Figures 6–8, respectively. The calculation results are compared with numerical simulation results. In the pictures, the solid lines refer to equation calculations and dashed lines refer to numerical simulations. Generally, in oil field, the ultimate sweep efficiency is obtained when the water cut of the reservoir reaches 98%. The comparisons show that when water cut of the model comes to 98%, the sweep efficiency of commingled injection is 46.2%, while it is 58.1% for zonal injection; and the oil recovery of commingled injection is 27.8%, while it is 34.9% for zonal injection. Figures 9 and 10 illustrate the comparison of water saturation profile between commingled injection and zonal injection after 12 months and 60 months, respectively.Figure 6
Results comparison of water cut versus time in commingled injection and zonal injection.Figure 7
Results comparison of sweep efficiency versus time in commingled injection and zonal injection.Figure 8
Results comparison of oil recovery versus time in commingled injection and zonal injection.Figure 9
Comparison of water saturation between commingled injection and zonal injection after 12 months.Figure 10
Comparison of water saturation between commingled injection and zonal injection after 60 months.Since the numerical simulator has considered the effect of gravity and capillary forces, the values of water cut calculated by numerical simulation are a bit higher than the results of equation calculations, and the breakthrough time of the simulation is earlier than that of equation calculation. However, the general trends of these curves are in complete agreement, which verify the equations of multilayers.


### 5.3. Results Analysis

The comparison results of the two kinds of water injection show the following: in zonal injection, the producing degree of the less permeable layers increases apparently and achieves a much higher sweep efficiency and oil recovery compared with commingled injection. By using zonal injection program, it can restrain the monolayer breakthrough of higher permeable layers and solve interlayer conflicts. Therefore, a balanced frontal advance is achieved. In oil field development, engineers shall subdivide the layer series with zonal injection based on the reservoir’s property differences, which could significantly reduce heterogeneity and postpone the time of water breakthrough. As a consequence, better oil production is achieved.


## 5.1. Water Frontal Position

Use the equations for single layer to calculate the time of water breakthrough of the five layers by (19) for both commingled injection and zonal injection. T1 is the breakthrough time of the most permeable layer (layer 1). T1 is equal to 214.1 days in the commingled injection program while it is 428.2 days in the zonal injection program because of the pressure variation. The water frontal position of all layers on T1 is calculated by (18). The corresponding sweep volume of each layer is calculated by (26) and (28). The results of equation calculation are listed in Table 2 and results of numerical simulation are listed in Table 3. We plot the water frontal position on T1 for the two injection programs (Figures 4 and 5). The calculation results are compared with numerical simulation results, where the former are plotted on the top and the latter are plotted right below. The comparisons show that the simulation results match the calculation results well, which verify the equations of single layer.Table 2
The results of analytical model in commingled injection and zonal injection.
Layers Commingled injection Zonal injection Breakthrough time,T (day) Frontal position onT1, xf (m) Sweep volume onT1, A (m3) Breakthrough time,T (day) Frontal position onT1, xf (m) Sweep volume onT1, A (m3) Layer 1 214.1 200.0 7985.7 428.1 200.0 7985.7 Layer 2 273.3 157.3 6280.6 546.7 157.3 6280.6 Layer 3 4441.8 9.8 392.3 1480.6 58.7 2342.1 Layer 4 1480.6 29.4 1174.5 493.5 173.9 6945.4 Layer 5 2961.2 14.7 588.1 987.1 87.7 3502.8Table 3
The results of numerical model in commingled injection and zonal injection.
Layers Commingled injection Zonal injection Breakthrough time,T (day) Frontal position onT1, xf (m) Sweep volume onT1, A (m3) Breakthrough time,T (day) Frontal position onT1, xf (m) Sweep volume onT1, A (m3) Layer 1 182.0 200.0 7985.7 386.0 200.0 7985.7 Layer 2 231.0 158.0 6308.7 493.0 156.0 6228.9 Layer 3 3567.0 12.0 479.1 1238.0 56.0 2236.0 Layer 4 1214.0 31.0 1237.7 441.0 172.0 6867.7 Layer 5 2502.0 16.0 638.9 874.0 86.0 3433.9Figure 4
Results comparison of water frontal position of commingled injection onT1.
(a)
Results of equation calculation (b)
Results of numerical simulationFigure 5
Results comparison of water frontal position of zonal injection onT1.
(a)
Results of equation calculation (b)
Results of numerical simulation


## 5.2. Water Cut and Sweep Efficiency

Use the equations for multiple layers to calculate water cut of the model in both of the injection programs by (5) and calculate sweep efficiency by (29) and (32), and oil recovery can be obtained. The curves versus time are shown in Figures 6–8, respectively. The calculation results are compared with numerical simulation results. In the pictures, the solid lines refer to equation calculations and dashed lines refer to numerical simulations. Generally, in oil field, the ultimate sweep efficiency is obtained when the water cut of the reservoir reaches 98%. The comparisons show that when water cut of the model comes to 98%, the sweep efficiency of commingled injection is 46.2%, while it is 58.1% for zonal injection; and the oil recovery of commingled injection is 27.8%, while it is 34.9% for zonal injection. Figures 9 and 10 illustrate the comparison of water saturation profile between commingled injection and zonal injection after 12 months and 60 months, respectively.Figure 6
Results comparison of water cut versus time in commingled injection and zonal injection.Figure 7
Results comparison of sweep efficiency versus time in commingled injection and zonal injection.Figure 8
Results comparison of oil recovery versus time in commingled injection and zonal injection.Figure 9
Comparison of water saturation between commingled injection and zonal injection after 12 months.Figure 10
Comparison of water saturation between commingled injection and zonal injection after 60 months.Since the numerical simulator has considered the effect of gravity and capillary forces, the values of water cut calculated by numerical simulation are a bit higher than the results of equation calculations, and the breakthrough time of the simulation is earlier than that of equation calculation. However, the general trends of these curves are in complete agreement, which verify the equations of multilayers.


## 5.3. Results Analysis

The comparison results of the two kinds of water injection show the following: in zonal injection, the producing degree of the less permeable layers increases apparently and achieves a much higher sweep efficiency and oil recovery compared with commingled injection. By using zonal injection program, it can restrain the monolayer breakthrough of higher permeable layers and solve interlayer conflicts. Therefore, a balanced frontal advance is achieved. In oil field development, engineers shall subdivide the layer series with zonal injection based on the reservoir’s property differences, which could significantly reduce heterogeneity and postpone the time of water breakthrough. As a consequence, better oil production is achieved.


## 6. Conclusion

(1)
An analytical model is developed to explain the waterflooding process in the condition of constant pressure, which substitutes for the condition of constant rate in traditional methods. We deduce the theoretical relationship between the frontal position of waterflooding and water saturation of producing well versus time to describe the dynamic changes of these variables. (2)
By applying this model from one single layer to multilayers, this method gives a quantitative expression of sweep efficiency of the three stages during waterflooding, which can be used to accurately calculate the sweep efficiency of heterogeneous reservoirs and evaluate the degree of waterflooding in multilayer reservoirs. (3)
By comparing the performances of commingled injection and zonal injection in the case study, zonal injection works better than commingled injection in respect of sweep efficiency and oil recovery and has a longer period of water free production. So, it is essential to subdivide the layer series with zonal injection in order to reduce heterogeneity and postpone the time of water breakthrough.



---

*Source: 6273492-2016-11-30.xml*