## Abstract

The interlayer interference is very serious in the process of water flooding development, especially when the reservoir adopts commingling production. The implementation of various interlayer interference mitigation measures requires that the production performance parameters and remaining oil distribution of each layer of the reservoir should be clearly defined, and the accurate production splitting of oil wells is the key. In this paper, the five-spot pattern is simplified to a single well production model of commingled production centered on oil well. The accurate production splitting results are obtained through automatic history matching of single well production performance. The comparison between the calculation results of this method and that of reservoir numerical simulation shows that the method is simple, accurate, and practical. In the field application, for the multilayer commingled production reservoir without accurate numerical simulation, this method can quickly and accurately realize the production splitting of the reservoir according to the development performance data.

## 1. Introduction

Due to the influence of reservoir physical property difference among layers, the interlayer interference is serious in the process of water flooding development when the reservoir adopts commingling production, and the difference of producing degree among layers is large. The implementation of measures to reduce the interlayer interference, such as separate layer water injection and formation reorganization, needs to make clear the production performance of each layer and the distribution of remaining oil. Accurate oil well production splitting is the basis of dynamic analysis, development effect evaluation, and remaining oil research of each layer. At present, the commonly used methods of production splitting include liquid production profile coefficient method [1–3], effective thickness method [4], static land splitting coefficient (KH) method [5–7], and dynamic splitting equation method [8–10]. These methods only use effective thickness, effective permeability, and other formation parameters to split production, and seldom consider the influence of remaining oil saturation difference between layers at different time on production splitting, which cannot truly reflect the actual production situation of each underground layer, resulting in low accuracy of production splitting results. Although the reservoir numerical simulation method can obtain more accurate liquid and oil production of each layer at different times, it needs to adjust the parameters manually with long cycle and heavy workload. Tian et al. and Gu et al. [11, 12] introduced the catastrophe method into the production splitting of multilayer combined production wells. They established a new production splitting method of multilayer combined production gas wells considering reserve characteristics, development characteristics, and geological characteristics and applied it in various oilfields [13–15]. However, this method does not fully consider the impact of development dynamics on production splitting and has some limitations.

The automatic history matching method can realize the parameter inversion and adjustment of the numerical simulation model and ensure the accuracy of the numerical simulation to the development performance prediction after effective test and correction. In this paper, the production splitting method of multilayer commingled production reservoir based on automatic history matching is established. Based on automatic history matching, the production of multilayer oil wells is further divided, and more accurate production splitting results are obtained.

## 2. Establishment and Solution of Single Well Production Mathematical Model in Multilayer Commingled Production Reservoir

### 2.1. Establishment of Mathematical Model

- (1)
Filtration equation

*μ*m

^{2}), $Kro$ and $Krw$ are the relative permeability of oil phase and water phase, $\mu o$ and $\mu w$ are the viscosity of oil phase and water phase (mPa·s), $Po$ and $Pw$ are the pressure of oil phase and water phase ($10\u20101MPa$), $So$ and $Sw$ are the saturation of oil phase and water phase, $\rho o$ and $\rho w$ are the density of oil phase and water phase (kg/m

^{3}), $\varphi $ is the formation porosity, $qo$ is the volume flow rate of oil injected or produced in the unit in unit time under standard ground conditions (m

^{3}/day), and $qw$ is the volume flow rate of water injected or produced in the unit in unit time under standard ground conditions (m

^{3}/day).

- (2)
Auxiliary equation

- (3)
Initial condition

- (4)
Boundary condition

### 2.2. The Solution of Mathematical Model

#### 2.2.1. Well Treatment Method

The treatment of the well in the mathematical model is the treatment of the boundary conditions in the model, and the production mode of constant liquid production is adopted in the mathematical model of multilayer commingled production in this paper. For the mathematical model, when a well passes through multiple layers vertically, the given liquid production at the wellhead is the sum of the production of all the layers that the well passes through. Therefore, we need to allocate the given liquid production to each layer, that is, liquid rate splitting.

*μ*m

^{2}), $Krw$ is the relative permeability of water phase, $Kro$ is the relative permeability of oil phase, $A$ is the flow cross-section area of the $k$ layer (m

^{2}), $\lambda w$ is the mobility of water, $\lambda o$ is the mobility of oil, $Sw$ is water saturation, and $i$ is the $i$

*-*th grid.

#### 2.2.2. Solution of Multilayer Commingled Production Reservoir Model

There are four unknowns ($Sw$, $So$, $Pw$, and $Po$) in the mathematical model of multilayer reservoirs, and the saturation term in the equation should be eliminated first when solving the pressure. In this paper, the implicit pressure and explicit saturation [16–18] (IMPES) method is used to establish the difference equations and the numerical method is used to solve them.

## 3. Establishment and Solution of Optimization Model for Single Well Automatic History Fitting

Based on the single well production mathematical model in multilayer commingled production reservoir, the optimization theory is introduced, the fitting parameters and fitting indexes are reasonably selected, the objective function of single well automatic history fitting is established, and the appropriate algorithm is selected to solve the problem, so as to complete the single well automatic history fitting.

### 3.1. The Establishment of Objective Function

### 3.2. Selection of Fitting Parameters

Injection-production well spacing is the most important fitting parameter. According to the actual dynamic production data of oil wells, the control reserves of single well are calculated by using water drive curve method, and then, the control area of single well is obtained, which is taken as the constraint condition of fitting injection-production well spacing.

^{4}t), $Np$ is the cumulative oil production (10

^{4}t), $a1$ is the intercept of the straight line segment (which is a constant related to rock properties or fluid properties), and $b1$ is the slope of straight line section (which is a constant related to geological conditions, well pattern arrangement, and oilfield management measures).

In the linear regression of A-type water drive curve, the occurrence time of the straight line should be considered. In the early stage of production, the first half of A-type water drive curve drawn by production data is not linear. Only when the oilfield enters the stage of comprehensive water drive development and the oil well reaches a certain water cut, the water drive curve at this stage is approximate to the inclined straight line, which has a certain representativeness.

^{4}t), $A$ is oil bearing area (m

^{2}), $H$ is average effective thickness (m), $\varphi $ is average effective porosity, $Soi$ is average original oil saturation of reservoir, $\rho osc$ is average surface crude oil density (t/m

^{3}), $Boi$ is volume coefficient of crude oil under original formation pressure (m

^{3}/m

^{3}), and $ER$ is oil recovery.

^{3}/m

^{3}).

In the automatic history fitting of single well, the interlayer heterogeneity becomes the main factor affecting the water cut of single well when the production well is producing with constant liquid rate. The interlayer heterogeneity is mainly reflected in the permeability difference of each small layer. Therefore, the permeability of each small layer is selected as one of the fitting parameters of single well history fitting.

*μ*m

^{2}) and $n$ is the total number of all layers.

### 3.3. Solution of Optimization Model

After establishing the objective function of single well automatic history fitting, it needs to choose the appropriate algorithm to solve. In this paper, the genetic algorithm [21–24] is chosen which has strong compatibility, fast solving speed, and is suitable for large-scale parallel computing.

The operation process of genetic algorithm is a typical iterative process. Combined with the single well production mathematical model and single well automatic history fitting optimization model studied in this paper, the algorithm flow is shown in Figure 1.

The termination condition is whether the maximum iteration number of the genetic algorithm is achieved or whether the stagnation iteration number of the genetic algorithm is achieved. When the adjustable parameter $M$ (Equation (22)) is in a certain range, a group of parameters $Mc$ is found to minimize the objective function $J$ (Equation (13)), and the solution $Mc$ is the optimal solution of the objective function of single well automatic history fitting.

## 4. Example Verification

### 4.1. Model Parameter Setting

The mathematical model of multilayer commingled production reservoir is established, and P1 oil well is selected as the research object. There are five small layers in the well, and the basic parameters of the physical properties of the model reservoir are set as follows: viscosity of crude oil is 15 mPa·s, viscosity of formation water is 0.5 mPa·s, crude oil density is 850 kg/m^{3}, and formation water density is 1000 kg/m^{3}. The production time and liquid production of P1 well are shown in Table 1. At the same time, the reservoir numerical simulation model is established for simulation calculation, and the accuracy of this method is verified by the reservoir numerical simulation results.

According to the production data of P1 well, the cumulative oil production and cumulative water production data of the well are used to draw the A-type water drive curve, as shown in Figure 2.

According to the linear regression of A-type water drive curve, the calculated parameters and results are shown in Table 2. The control area of single well in the table can be used as one of the limiting conditions for fitting injection production well spacing. The well spacing calculated from the single well control area in the table is taken as the initial value of fitting.

### 4.2. Calculation Results

The objective function is solved by genetic algorithm, and the relevant genetic parameters are set as follows. The number of genetic independent variables (nvars) is 6, and the upper and lower limits of genetic variables are given. Fitness function (fitnessfun) of genetic algorithm is the objective function established in this paper. The generation of genetic algebra is set to 50. StallGenLimit is set to 50.

The automatic history fitting process and fitting results of P1 well are shown in Figure 3. It can be seen from the figure that in the process of applying genetic algorithm to single well automatic history fitting, the observed data of daily oil production, daily water production, and water cut of oil wells will approach the actual production data with the process of genetic variation. It can be seen from the figure that the water cut curve of P1 oil well obtained by the model simulation calculation after fitting has little difference with the actual water cut dynamic change of oil well, which can verify the feasibility and practicability of the single well automatic history fitting method.

Figure 4 shows the daily oil production comparison between the actual model and the model after single well automatic history fitting. It can be seen from the figure that the production performance of oil wells in multilayer reservoir after single well automatic history fitting is highly consistent with the production performance of actual model, which proves that the permeability of each layer after fitting is quite close to the real permeability value of each layer.

The daily oil production of each small layer during the development period was calculated by using the fitted multilayer combined production single well flow model and compared with the numerical simulation results to verify the split production results of P1 well. Figures 5–7 show the daily oil production comparison of layer 1, layer 3, and layer 5, respectively. In the diagram, the red curve represents the time-varying daily oil production of layer calculated by numerical simulation, and the blue curve represents the daily oil production of layer after fitting. It can be seen that in the whole production and development stage, the daily oil production situation of each layer is almost consistent with the actual daily oil production situation of each layer, which proves that the splitting result of the fitted model is highly accurate and can be used as a reliable basis for clarifying the distribution of remaining oil in each layer. At the same time, it provides help for the design of production and injection allocation plan or the improvement of remaining oil tapping measures in the later stage of production and development.

## 5. Conclusions

- (1)
The mathematical model of single well production in multilayer commingled production reservoir is established, and the optimization theory is introduced. The objective function of single well automatic history fitting for multilayer coproduction reservoir is established, and the genetic algorithm is selected to solve the objective function

- (2)
A production splitting method based on automatic history fitting of single well is proposed. Based on single well production mathematical model of multilayer commingled production reservoir, automatic history fitting of single well is carried out, and the reservoir parameters of each layer near the single well are accurately reflected, which makes the production splitting result more accurate

- (3)
Through the example verification, it can get more accurate daily oil production of layer by using single well automatic history fitting for a production well in multilayer commingled production reservoir and then split the production. In the field application, this method can quickly and accurately realize the production splitting of the reservoir according to the development performance data

## Data Availability

All the data are in the paper.

## Conflicts of Interest

The authors declare no conflicts of interest.

## Acknowledgments

This work was supported by the National Natural Science Foundation of China with No. 51974343 and Project funded by China Postdoctoral Science Foundation with No. 2021M703588. The first author would like to acknowledge the College of Petroleum Engineering at China University of Petroleum (Qingdao) for providing experimental equipment.