Application and Study on Relationship Between Relative Permeability Ratio and Water Saturation Curve
Li Jun^{1, 3, *}, Kang Kai^{1}, Peng Xiaodong^{2}, Lu Yan^{2}, Jiang Lifu^{1}
^{1}CNOOC Ltd. Tianjin, Tianjin, China
^{2}CNOOC Ltd. Zhanjiang, Zhanjiang, Guangdong, China
^{3}Engaged in Reservoir Engineering and Reserve Evaluation, CNOOC Ltd. Tianjin, Tianjin, China
Email address:
To cite this article:
Li Jun, Kang Kai, Peng Xiaodong, Lu Yan, Jiang Lifu. Application and Study on Relationship Between Relative Permeability Ratio and Water Saturation Curve. International Journal of Oil, Gas and Coal Engineering. Vol. 4, No. 3, 2016, pp. 24-29. doi: 10.11648/j.ogce.20160403.11
Received: March 29, 2016; Accepted: June 22, 2016; Published: June 28, 2016
Abstract: In the process of oil field development, the law of water cut is usually described by exponent model, but during the early water breakthrough stage and ultra-high water cut stage, relative permeability ratio of oil to water (K_{ro}/K_{rw}) and water saturation (S_{w}) is no longer linear relations in semilog coordinate, so it has an important influence on water cut increasing rate, theoretical production index, oil displacement efficiency, relationship between recovery and water cut etc. Based on relation curve character of relative permeability ratio of oil to water (K_{ro}/K_{rw}) and water saturation (S_{w}), this paper proposes a new function model in order to reduce prediction error, such as production index and water cut increasing rate. According to nonlinear properties of the model, Newton algorithm is supposed to be used for optimization solution. Finally, take relative permeability data of P Oilfield as an example to testify, and results show that fitting precision of the model is quite high and it possess nicer applicability. On this basis, the relationship between water cut and recovery, relationship between water cut increasing rate and recovery are derived. By comparing water cut and water cut increasing rate with recovery degree, it can be used to guide the development as the foundation of geology and reservoir of P Oilfield.
Keywords: Relative Permeability Ratio, Water Saturation, Optimization Solution, Newton Algorithm
1. Introduction
Evaluate water production law of oil production well is a key point in oilfield development. The precision of water cut forecast will affect reservoir engineer’s right decision. Water cut percentage f_{w} is a significant evaluation parameter to describe water production law, without considering the influence of gravity and capillary pressure, according to the flow equation, it is calculated as follows:
(1)
The water cut percentage f_{w} of ground which is converted by volume coefficient is calculated as follows:
(2)
Equation (2) shows that viscosity ratio of oil to water can be considered as a constant value for a specific reservoir, so relative permeability ratio determination is important to water cut. How to define functional relation between K_{ro}/K_{rw} and S_{w} will decide forecast precision of water cut in oilfield development.
2. Exponent Model Establishment
In semi-log coordinate relation curve of K_{ro}/K_{rw} and S_{w} has two curved ends and one straight line in the middle [1-15], which can be illustrated in Figure 1. Select the middle straight line of the curve and calculated by empirical equation (3) [1, 2, 3] below:
(3)
Ren [4] proposed characterizing relationship of and S_{w} by equation (4), established a new method to describe the various forms of water flooding curve with equation; Miao and Liu [5] et al. characterized 3 kind of the relationship between water cut and recovery as well as the type of water cut increasing rate and recovery; Gao and Li [6, 7] et al. proposed using two straight lines method for data matching separately in medium water cut period and high water cut period, calculated the recoverable reserves of Xingbei development area of Daqing oil field and verified the cumulative oil production and water cut. Song and Bing Shao [8, 9] et al. also made enormous effort on model of K_{ro}/K_{rw} and S_{w}, derived a new water flooding characteristic curve which is more applicable for high water-cut oil field by the average water saturation equation which proposed respectively by Buckley-leverett and Welge.
(4)
However the middle straight line can be only applied to medium water cut period, which is not suitable to relative low or high water saturation period. It will reduce accuracy of water production forecast. Two straight lines method take data matching separately in medium water cut stage and high water cut stage into consideration and improved prediction accuracy. But unified functional relation is not defined and it cannot match data in low water saturation period. The expression of power function of relative permeability which define functional relation between K_{ro}/K_{rw} and S_{w} improve the fitting precision, but it could not solve the problem radically. While in actual oilfield production, it is in low or high water saturation period in a long time. So according to the curve characteristic (two curved ends and straight line in the middle), a new exponential model as equation (5) is established to fit relation of K_{ro}/K_{rw} and S_{w} in this paper:
(5)
3. Exponent Model Solution
Equation (5) has strong nonlinearity, to solving parameters like a, b, c and d, function F(X) is proposed as follows:
(6)
Thereinto,
Apparently F(X) has continuous second-order partial derivative, with positive definite Hesse matrix. Therefore use optimal principle and solving parameters by Newton algorithm [9].
Solution steps are as follows:
(1) Select iteration initial point , assume k=0; calculate gradient matrix ；
(2) Calculate Hesse matrix and its inverse matrix ；
(3) Calculate and ；
(4) Test whether the condition is satisfied: if so， is optimal solution；otherwise, set k=k+1 and repeat step (2).
4. Case Demonstration and Comparison
Take relative permeability data of P oilfield as an example, which is illustrated in Table 1 and Figure 2, this paper studied relation between water saturation and relative permeability ratio of the area.
S_{w} | K_{rw} | K_{ro} |
0.22 | 0 | 1 |
0.24 | 1.40E-05 | 0.907 |
0.259 | 0.00011177 | 0.819 |
0.279 | 0.00037723 | 0.736 |
0.299 | 0.00089417 | 0.659 |
0.318 | 0.00174643 | 0.587 |
0.338 | 0.00301783 | 0.52 |
0.357 | 0.00479221 | 0.458 |
0.377 | 0.00715338 | 0.401 |
0.397 | 0.01018519 | 0.348 |
0.416 | 0.01397145 | 0.3 |
0.436 | 0.018596 | 0.257 |
0.456 | 0.02414266 | 0.217 |
0.475 | 0.03069527 | 0.181 |
0.495 | 0.03833765 | 0.15 |
0.514 | 0.04715364 | 0.121 |
0.534 | 0.05722705 | 0.097 |
0.554 | 0.06864172 | 0.076 |
0.573 | 0.08148148 | 0.057 |
0.593 | 0.09583016 | 0.042 |
0.613 | 0.11177158 | 0.03 |
0.632 | 0.12938958 | 0.02 |
0.652 | 0.14876797 | 0.012 |
0.671 | 0.1699906 | 0.007 |
0.691 | 0.19314129 | 0.003 |
0.711 | 0.21830387 | 0.001 |
0.73 | 0.24556216 | 0 |
0.75 | 0.275 | 0 |
A calculation program is written for data matching of relative permeability data for three models (, , ) respectively. Matching results are shown in Table 2 ~ Table 3 and Figure 3~Figure 5.
Figure 3 shows that model is suitable to data matching in middle straight line, while fitting accuracy is low in low or high water saturation period. Especially in high water saturation period, the fitting error can reach 394%; Figure 4 indicates that comparing to model , fitting accuracy of model has been improved, but still low; however model has higher fitting accuracy and relative lower error. It shows that can describe the relationship between relative permeability ratio of oil to water () and water saturation () more accurately.
5. Application
5.1. Relationship Between Water Cut and Recovery
According to the fitting result of equation (4), the water cut percentage f_{w} can be calculated as follows:
(7)
The relationship between recovery () and water saturation () is as follows:
(8)
In this case, we set to be A, to be B, equation (7) is converted to equation (8) as follows:
(9)
5.2. Relationship Between Water Cut Increasing Rate and Recovery
Water cut increasing rate is defined as an increase value of water cut while 1% OGIP is produced, it is calculated as follows:
(10)
Water cut increasing rate () is calculated as follows by formula (9):
(11)
We can plot the relationship between water cut increasing rate () and recovery ().
5.3. Case Analysis
In regional structure, P Oilfield lies in the east of Bohai Bay Basin and northeast end of middle section of Bonan Low bulge, and is developed on Tancheng-Lujiang fault zone. Bonan low bulge is next to Bozhong and Bodong depression in the north, with southwest-southeast surrounded by Yellow River Estuary and Miaoxi depression, hence it refers to one of the most favorable oil-gas enrichment areas in Bohai Sea, and owns superior geological conditions of oil-gas accumulation. Cenozoic and Mesozoic strata are mainly revealed in drilling of Area 1/3/8/9 of P Oilfield, which falls into Neogene Pingyuan Formation, Minghua Formation, Guantao Formation and Palaeogene Dongying Formation from top down according to lithology, electric characteristics, paleontological analysis and regional stratum comparison. Main Oil Saturation formation develops in Lower Minghua Formation and Guantao Formation of Neogene. Since Neogene, this area enters peneplanation period and forms a set of sedimentary system centered on fluvial facies. Guantao Formation is braided river deposit, and Lower Minghua Formation belongs to meandering stream deposit, with reservoir lithology referring to terrigenous clastic rock of fluvial facies deposit. According to the principle of "cycle comparison, and step control" and in combination with drilling of producing well, stratum thickness, development of oil reservoir and relationship with seismic reflection horizon, Oil Saturation target intervals are divided into 13 oil formations (respectively named as L00, L10, L20, L30, L40, L50, L60, L70, L80, L90, L100, L110 and L120), in which, main Oil Saturation target intervals of Lower Minghua Formation are divided into 5 oil formations (L00-L40 oil formations) and those of Guantao Formation are divided into 8 oil formations (L50-L120 oil formations). L50-L100 oil formations is main oil producing in P Oilfield.
According to the actual oil-water parameters in P Oilfield (Table 4), it can be concluded that the relationship between water cut and recovery degree and water cut increasing rate and recovery degree (Figure 6~Figure 7).
μ_{w} | 0.5 |
μ_{o} | 30 |
B_{w} | 1 |
B_{o} | 1.1177 |
(1) Relationship between water cut and recovery
Figure 6 shows that the error between actual water cut increasing curve and theoretical curve is bigger. The theoretical water cut rising curve increases as "S" form, the actual water cut increasing curve increase as the early stage of convex type.
(2) Relationship between water cut increasing rate and recovery
Figure 7 shows that the error between actual water cut increasing rate curve and theoretical curve is still worse. The theoretical water cut increasing rate curve is usually increased first and decreased then, the actual curve shows a high water cut increasing rate in the early stage, 2008 and 2011, the water cut increasing rate is negative because of the new platform (B/C/D) put into operation as well as the oil spill risk lead the oilfield to shut down. At present, the water cut increasing rate of Bohai P Oilfield is 5.3%.
Figure 6~Figure 7 show that P Oilfield has the poor development effect. According to the analysis, the main reason include:(1) Large sections with one-set series of strata development and severe interlayer interference. Oil wells mainly produce the Guantao Formation (8 oil formations and 35 sub-layers), with a vertical span of the production horizon 300-350m, average well control 130×10^{4} m^{3}. There is severe interlayer interference in the production process. Moreover, according to the oil well production data, with the increase of production thickness, the production capacity of oil wells shows a downward trend. The interlayer interference coefficient for separate and commingled production of the Upper and Lower Guantao Formation is 40%. The interlayer contradiction in large sections of commingled production is obvious. (2) Imbalanced water absorption of injection wells, with severe breakthrough of injected water. The water injection time in P Oilfield is quite long, and the monolayer breakthrough of injected water becomes more serious. The water absorption profile of water injection wells shows that the imbalance of water absorption in sub-layers causes the imbalance of the produced fluid profile of production wells, as well as the large difference in water cut among all layers and in water flooding effects. According to the problems which affect the development of P Oilfield, we can improve the effect of water flooding by means of profile control and optimization of water injection in short-term, make comprehensive adjustment for long term goal.
6. Conclusion
(1) In low or high water saturation period, traditional exponent model presents low prediction accuracy. This paper proposed a new exponent model, which can establish function relation between relative permeability ratio and water saturation in a high accuracy.
(2) According to new model, the relationship between water cut and recovery, relationship between water cut of increasing and recovery are derived.
(3) By comparing the water cut and recovery degree, water cut increasing rate and recovery degree, it can be seen that the development of P Oilfield is poor. At the same time, this paper puts forward some corresponding measures and suggestions.
References