The bottom water reservoir of x oilfield had the characteristics of high crude oil viscosity, low oil column height and sufficient bottom water energy. The development mode of single sand body of horizontal well was adopted, which caused water cut rise rapidly in horizontal production wells and bottom water coning irresistible. When the oil well entered the ultra-high water cut stage, the sweeping volume tended to be stable, and the main way to improve oil recovery rate was to improve oil displacement efficiency. At present, x Oilfield Group was in the stage of ultra-high water cut and 68% of production wells had water cut more than 90%. The dynamic characteristics of production wells showed that the water cut of horizontal wells in the bottom water reservoir rised rapidly, the recoverable reserves were mainly produced after high water cut period, and the cumulative oil production in the stage of extra-high water cut accounts for 60% - 70% of the recoverable reserves. Therefore, the development characteristics of the bottom water reservoir showed that the extra-high water cut stage was an important development stage of this reservoir type.
At present, the problem was that the common empirical formula of recovery rate in oilfield was mainly obtained through statistical analysis, and all calculation parameters were applicable within a certain range (Permadi & Jayadi, 2010) . However, the viscosity, permeability and well pattern density of crude oil in x oilfield were not applicable to the empirical formula, so there was no reliable conclusion on recovery rate (Kuchuk, 1988) . In the actual development process of the oilfield, the local scour multiple index was as high as tens of thousands (Gao, Jiang, Wang, et al., 2016; Zheng, 1993) . With the increase of pore volume injection multiple, the oil displacement efficiency in the swept area could be further improved, due to large differences of residual oil saturation acquired by laboratory test and actual condition. Under the condition of high scour, the end point of residual oil saturation must shift to the right, so the relative permeability curve needs to be further modified, and the final understanding of water drive efficiency needs further theoretical study.
2. Theoretical Model Establishment and Analysis of Bottom Water Reservoir
2.1. Quantitative Characterization of Bottom Water Ridge
The flow of horizontal wells in bottom water reservoir was shown in Figure 1. The hypothesis was: 1) Oil-water two-phase fluid seepage in reservoir; 2) Uniform reservoir thickness; 3) Reservoir fluid percolation conforms to Darcy’s law; 4) Reservoir rocks were anisotropic and heterogeneous; 5) Regardless of capillary pressure; 6) Regardless of gravity (Fan, 1993; Zheng, Xu, & Chen, 2013) .
According to the basic formula of the balance principle of statics:
The relationship between the height of water ridge and the radius of water ridge could be obtained through the corresponding differentiation and integration:
Figure 1. Bottom water ridge of horizontal well in bottom-water reservoir.
By integrating, the following formula was proposed:
Through the above formula, the height under different water ridge radius could be calculated, as well as the distribution and remaining oil quantity.
2.2. Vector Characteristics and Calculation Method of Directional Reservoir Permeability
Research and production practice showed that permeability had obvious directionality, so permeability was vector (Xin, 2011; Li, Zhou, Xiong, et al., 2015) . The directional difference of permeability will directly affect the well arrangement, well spacing and the technical design of various development measures such as artificial fracturing. In view of the above important role of vector permeability, based on previous studies, this paper introduced the reasons for the vector characteristics of rock permeability, and analyzed some understandings that were easy to be confused when understanding and applying the vector characteristics of permeability (Liu & Hu, 2011, Chen & Tao, 1997; Bing, 2012) .
As shown in Figure 2, OB was the azimuth reference line (the azimuth of OB was 0˚, generally taking the due east direction as the azimuth reference line); in the derivation process, the X and Y axes of the coordinate system were consistent with the direction of the extreme permeability, and the extreme permeability in the x-axis direction was Kmx. The extreme permeability in y direction was Kmy; was the displacement pressure gradient in n direction; A was the area of the vertical seepage section perpendicular to ; was displacement of fluids. The displacement effect of x component and y component on fluid. The above principle was called equivalent displacement principle. If was used to express the flow of fluid through section a in X direction under the action of . The flow rate of fluid through section a in Y direction under the action of was expressed by , According to the equivalent displacement principle, Flow Q through seepage section a under the action of, Should be
Figure 2. Derivation of permeability calculation model.
Derivation of quantitative calculation model of vector permeability. Using Kn to express permeability in n direction; Fluid seepage velocity in the direction of Vn: μ was the viscosity of the fluid, according to Darcy’s law.
Then the flow through seepage section A is:
As shown in Figure 3, Use Ax to express the effective seepage area of section A in z direction (that was the seepage area perpendicular to x direction); was the azimuth of direction n; represents the azimuth of the extreme permeability Kmx, that is, the azimuth of the x-axis. Because the two sides of the angle and the two sides of the angle were perpendicular to each other, so be equal to , so:
Similarly, the effective seepage area of interface a in Y direction (that was, the seepage area perpendicular to Y direction was
The component in the X direction was:
Under the action of , the seepage velocity of fluid passing through section a in X direction was:
You could get it by combining
The above formula was the quantitative calculation model of anisotropic permeability.
Where Kn was the permeability in n direction; Kmx was the extreme permeability in the x-axis direction; Kmy was the extreme permeability in the y-axis direction; was the azimuth of direction n; was the azimuth of the extreme permeability Kmx.
2.3. Distribution Pattern and Grid Division of Stream in Water Ridge Area of Bottom Water Reservoir
The gas oil ratio was low in X Oilfield, the water ridge profile of horizontal well could be simplified as oil-water two-phase flow chart, as shown in Figure 3. After the production well entered the ultra-high water cut period, the sweep area was stable. The flow pipe model could be used to simulate the oil-water two-phase flow in the sweep area through the quantitative characterization of water ridge shape (Sun, Zhou, Hu, et al., 2018, Yan, Li, Yin, et al., 2009; Chen, Bai, Lu, et al., 2015) .
The flow tube model assumed that the immiscible displacement process still followed the same streamline when the homogeneous fluid in the porous medium was stable. According to the flow line of single-phase fluid flow in the early stage, the simulation could be divided into multiple flow tubes. The linear displacement model was used to simulate the immiscible displacement in the flow tube. The dynamic of this region was simulated by combining the dynamic of all the flow tubes that made up this region in the same time. For each flow tube, the
Figure 3. Mesh generation of water coning in bottom-water reservoir.
simulation mainly included the solution of forward propulsion equation, the solution of total flow, and the solution of the relationship between total flow and the timing.
According to the two-dimensional seepage theory, the streamline distribution of stable condition was obtained. Assume that the streamline was constant, the displacement path of bottom water was obtained. In the process of displacement, the displacement space of bottom water was expanding along the horizontal direction.
In the flow tube model established by Higgins and Leighton, the flow tube was divided into N equal volume grids. For each production step, displacement process of simulation by advancing saturation of water drive front, until the water breakthrough occurs in the n grid. Finally, by combining the result of different flow tubes in the same injection time value, we could get the total dynamic state of all flow tubes in the same injection production well.
The solution steps of Higgins-Leighton method were as follows: the streamline and flow pipe distribution was determined in single-phase flow, for each flow tube, cumulative water injection multiple was given; Calculating the saturation distribution and average apparent viscosity distribution corresponding to each flow tube; Calculating the total flow volume in the flow pipe; Calculating the corresponding time value; Summarizing the results of each flow pipe in the same injection time value for getting the total flow dynamics. According to the streamline distribution of the bottom water area, equal volume grid was the principle of dividing.
2.4. Realization of Bottom Water Drive Model Calculation Program
The model was calculated by MATLAB software, this software was a high level programming language widely used in the field of engineering calculation and numerical analysis. It had strong numerical calculation function, which was easy to learn and easy to operate. The software was divided into one core solver and four modules, which included that data processing module, equal volume mesh module, water drive front saturation tracking module, post-processing module, data input and processing module.
3. Genetic Algorithm Coupled with Bottom Water Flooding Theory Model
3.1. The Form of the Permeability Curve Expressed
Usually, the phase permeability curve was expressed as follows:
where, Sw, water saturation; Swi, irreducible water saturation; Sor, residual oil saturation; , the initial maximum oil relative permeability; constant Cw, Co, refer to oil phase index and water phase index, . The phase permeability curve could be well fitted according to formula (13).
3.2. The Method of Genetic Algorithm (GA)
Genetic algorithm was a computational model simulating the natural selection and genetic mechanism of Darwinian evolution. It was a method to search the optimal solution by simulating the natural evolution (Yin, Zhao, Dong, et al., 2012) . Taking small bubbles of different colors and sizes as examples, in order to find bubbles of specific colors and sizes, it was necessary to optimize the initial number of small bubbles (survival and elimination), and then produce a new group of small bubbles after genetic calculation until the optimal solution was found.
In this paper, the genetic algorithm was coupled with the bottom water drive model. Different characteristic parameter of relative permeability curve was generated through initial population, and theoretical curves were calculated by the theoretical model between water cut and recovery degree. The optimal solution with minimum of forecast error was determined compared with the actual production data, as shown in Figure 4.
3.3. Calculation Results of Genetic Algorithm Coupled with Theoretical Model of Bottom Aquifer Drive
The residual oil saturation was 0.05 by the optimized infiltration curve, as shown in Figure 5, which was lower than the original oil saturation. This indicates that the oil displacement efficiency in the affected area can be greatly improved depending on the high PV in the ultra-high water cut stage of the bottom water reservoir. Further analysis shows that the optimized curve reflects the stronger hydrophilic property of the reservoir, as a result of isopermeability point shifts to the right.
Figure 4. Flow chart of the research coupling genetic algorithm and bottom-water drive mode.
Figure 5. Calculation result of optimized relative permeability curve based on genetic algorithm.
4. Method Application and Recognition
4.1. Distribution Characteristics of Micro Remaining Oil in Bottom Water Reservoir
Through the bottom water drive model, the distribution rule of remaining oil under the condition of constant pressure difference/constant flow rate can be realized. From Figure 6, it can be seen that the bottom water wave and area within the production range of oil layer present the “inverted V” shape, and the production well had the fastest water flow in the vertical direction, which was the dominant mainstream line, and its flow was far greater than the edge area.
It could be seen from the figure that there were obvious differences in oil displacement efficiency in the affected area. The vertical direction of production wells was reflected as “dominant channel”, and the scour multiple was higher than 1000 times. It could be seen from the figure that under the high porosity volume injection multiple, the oil displacement efficiency was the highest. Through this model, the quantitative characterization of the remaining oil spatial difference could be realized.
4.2. Effective Production Radius of Water Drive in Bottom Water Reservoir
Based on the quasi streamline method, the bottom water drive model of oil-water two-phase homogeneous horizontal well was established, which could realize the quantitative characterization of bottom water wave and remaining oil in the area, under the condition of constant pressure difference and constant flow rate, as well as the difference analysis of parameters such as injection multiple of pore volume and displacement efficiency in the affected area. On this basis, the concept of effective production radius of bottom water drive reservoir was put forward for the first time. Minimum distance of interference between wells was obtained. In addition, the theoretical template of effective production radius of horizontal well under different fluid properties and water avoiding height was given, as shown in Figure 7, which further enriches the development boundary of
Figure 6. Residual oil distribution of bottom water drive model.
Figure 7. The recoverable reserves of single well under different production radius with crude oil viscosity 142 mPa·s.
bottom water reservoir and provides guidance for further development and adjustment of oilfield.
1) Combined with Genetic algorithm and the flow tube simulation of the bottom water, the phase permeability curve was calculated from the production data of the ultra-high water cut stage for the bottom water reservoir.
2) The residual oil saturation by the optimized infiltration curve was lower than the original oil saturation, indicating that the oil displacement efficiency in the affected area can be greatly improved depending on the high PV in the ultra-high water cut stage of the bottom water reservoir.
3) Through the bottom water drive model, the distribution rule of remaining oil under the condition can be obtained, and the theoretical template of effective production radius of horizontal well under different fluid properties and water avoiding height was given, in order to provide guidance for further development and adjustment of oilfield.
 Chen, B., Bai, Z. H., Lv, J. W. et al. (2015). Study on Relationship between Sedimentary Microfacies and Water Flooding Pore Volume as Well as Remaining Oil: Taking Lower V Oil Layers in Shuanghe Oilfield as an Example. Journal of Xi’an Shiyou University (Natural Science Edition), 30, 53-58.
 Kuchuk, F. J. (1988). Pressure Transient Analysis for Horizontal Wells. In SPE Annual Technical Conference and Exhibition (SPE-18300-MS). Houston: Society of Petroleum Engineers. https://doi.org/10.2118/18300-MS
 Li, L. F., Zhou, F. X., Xiong, J. H. et al. (2015). Quantitative Evaluation Method of Water Flooded Degree of Oil Reservoirs Developed by Horizontal Wells and Its Application. Journal of Xi’an Shiyou University (Natural Science Edition), 5, 59-63.
 Permadi, P., & Jayadi, T. (2010). An Improved Water Coning Calculation for Horizontal well. In SPE Russian Oil and Gas Conference and Exhibition (SPE-133162-RU). Houston: Society of Petroleum Engineers.