#### 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.

#### 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

Figure 1 is a schematic representation of the stratified system. The production pressure between injection and production well is equal to a constant value: , where is the bottom pressure of injection well, MPa; is the bottom pressure of producing well, MPa; is the pressure drop from the injection well to producing well, MPa.

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 from to , where is the water saturation of displacement front, f, and 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.

**(a) Before water breakthrough**

**(b) At the time of breakthrough**

**(c) After water breakthrough**

Based on the material balance equation and Buckley-Leverett’s frontal advance theory, the total volume of water injected in a given time equalswhere is the flow rate of oil, m^{3}/s; is the width of each layer, m; is the height of each layer, m; is the porosity of matrix, f; is the displacement front position, m; is the initial position of displacement front, m; is the fractional flow of water, f.

The calculation of sweep efficiency in a single layer is given bywhere is the length of model, m. Before water breakthrough, average saturation is determined by fractional flow curve () [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:where is no longer a constant and should be determined by the saturation of the producing well . Consequently, it is necessary to calculate the waterflood front location before breakthrough and the saturation of the outlet face 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 position between injection well and the front location is given bywhere is the absolute permeability, *µ*m^{2}; is the water relative permeability, f; is the oil relative permeability, f; is the viscosity of water, mPas; is the viscosity of oil, mPas.

Based on the Buckley-Leverett frontal displacement theory, the water cut is expressed asSubstituting (5) into (4) and integrating, we obtain the pressure drop between the positions of :

The location of the injection well is treated as the position of , and the frontal displacement function can be written asEquation (7) can be deformed asFrom the position of to , the water saturation changes from to . Substituting (8) into (6), we getwhere is given by

##### 3.2. Time of Water Breakthrough

The pressure drop from the injection well to the displacement front location is defined as . Since the average saturation is a constant before water breakthrough, setting to value and substituting value into (9) yieldSubstituting (7) into (11), we getwhere 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:In the pure oil region, there is only one phase, and the pressure drop of this region isSubstituting (11) and (12) into (14), we getwhere is the ratio of oil viscosity and water viscosity, f. Taking the derivative of (13) and substituting it into (12), we obtainCombining (15) with (16) and integrating and rearranging the equation, the accumulated volume of water from to is obtained:Substituting (7) into (17), the position of displacement front is given by

When the front reaches the outlet face, equals ; substituting it into (18) and rearranging the equation, the time of water breakthrough can be obtained:

##### 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 iswhere 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 to .

Substituting (7) into (13) and taking the derivative, we getSolving (20) and (21) simultaneously, we getMaking the integral of (22), where is given byGetting (8) deformed and substituting it into (24),

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) Layer

(1)The layer The time of water breakthrough in layer can be calculated according to (19), and the sweep volume of layer is calculated by where is as in Welge’s equation [7]: (2)The other layers, for example, layer Substituting into (18), we can get the front position of layer . The sweep volume of layer is(3)Total sweep efficiency

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

(1)The layer It is the same as in Section 4.1(1).(2)Layers with breakthrough before layer , for example, layer The time of breakthrough in layer is calculated by (19). The saturation and water cut of layers that got breakthrough before can be calculated by After water breakthrough, Welge’s equation reforms as follows: Substituting (31) into (26), the sweep volume of layer can be obtained.(3)Layers without breakthrough, for example, layer This is the same as in Section 4.1(2).(4)Total sweep efficiency

#### 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 Table 1, and the relative permeability curves are plotted in Figure 3. By analyzing the data of logging, five layers are classified of the main layer .

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. is the breakthrough time of the most permeable layer (layer 1). 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 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 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.

**(a) Results of equation calculation**

**(b) Results of numerical simulation**

**(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.

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.

#### Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

#### Acknowledgments

This work was supported by China National Petroleum Co., Ltd., Major Science and Technology Projects (2012E-3304).