Next Article in Journal
A New Model for Estimating the Diffuse Fraction of Solar Irradiance for Photovoltaic System Simulations
Next Article in Special Issue
Study of the Effect of Clay Particles on Low Salinity Water Injection in Sandstone Reservoirs
Previous Article in Journal
Sectoral Electricity Consumption and Economic Growth: The Time Difference Case of China, 2006–2015
Previous Article in Special Issue
Investigation of the Flow Characteristics of Methane Hydrate Slurries with Low Flow Rates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Production Characteristics with Different Superimposed Modes Using Variogram: A Case Study of a Super-Giant Carbonate Reservoir in the Middle East

1
Research Institute of Petroleum Exploration and Development, PetroChina, Beijing 100083, China
2
School of Civil and Environmental Engineering, University of Science and Technology Beijing, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Submission received: 30 August 2016 / Revised: 15 January 2017 / Accepted: 6 February 2017 / Published: 18 February 2017
(This article belongs to the Special Issue Oil and Gas Engineering)

Abstract

:
Heterogeneity of permeability is an important factor affecting the production of a carbonate reservoir. How to correctly characterize the heterogeneity of permeability has become a key issue for carbonate reservoir development. In this study, the reservoirs were categorized into four superimposed modes based on the actual logging data from a super-giant heterogeneous carbonate reservoir in the Middle East. A modified permeability formula in terms of the variogram method was presented to reflect the heterogeneity of the reservoirs. Furthermore, the models of oil production and water cut were established and the analytical solutions were obtained. The calculation results show that the present model can predict the productivity of wells with different heterogeneous layers more accurately and rapidly. The larger the varigoram value, the stronger the heterogeneity of the reservoirs, and the faster the decline of production owing to a quicker reduction of formation pressure. With the increase in variogram value, the relative permeability of the oil phase is smaller and the water phase larger, and the water cut becomes larger. This study has provided a quick and reasonable prediction model for heterogeneous reservoir.

1. Introduction

There are two factors to influence oil recovery, which are geological conditions and development plans [1]. The former are inherent characteristics of reservoirs, such as heterogeneity, viscosity, wettability, capillary pressure, gravity, and so on [2,3]. The latter is related to the development design and technological measures, such as well spacing density, water injection, etc. [4]. The heterogeneity of vertical permeability is one of the most important factors to influence the reservoir oil recovery.
Carbonate rock is one of the most important oil and gas reservoir formations in the world, which reflects about half of the global reserves. The output has reached more than 60% of the total output. Compared with the conventional sandstone reservoir, the degree of exploration and development for carbonate reservoir is lower. Carbonate reservoirs have typical characteristics of dual media and serious permeability heterogeneity, which leads to poor development effects.
At present, there are two aspects to study permeability heterogeneity. On the one hand, the permeability heterogeneity is characterized by permeability variation coefficient and ratio of vertical permeability to horizontal permeability [5], and then the effects of heterogeneity are analyzed by theoretical calculations utilizing modified average permeability [6,7]. However, these methods can not characterize the randomness of the regional variables and represents the distribution of the strata heterogeneity. On the other hand, numerical simulation is used to study heterogeneity directly with different permeability layers [8,9]. Although the geological model established for numerical simulation is closer to the real reservoir, the simulation processes are complicated and time consuming [10].
The variogram is a method of moment estimation proposed by Motheron (1965), and it is a specific research tool in geological statistics [11], which can not only express the spatial structure with the regional variables, or characterize the randomness of the regional variables, but also reflects the change degree of the regional variables in a certain distance range and direction [12]. This method can reflect the change degree of the regional variables. In the process of the development of carbonate reservoirs, variogram analyses could be conducted to characterize the heterogeneity of their properties at the outcrop scale [13].
In this paper we, firstly, classified the heterogeneous reservoirs into four superimposed modes using the actual logging data from a super-giant carbonate reservoir in the Middle East. Then a new permeability model considering the heterogeneity was established with the variogram method. In addition, we verified our model with actual production data and the results show that the present model can predict the productivity of wells with different heterogeneous layers more accurately and rapidly. Finally, the analysis of variogram values was carried out to determine the relationship with oil production, relative permeability and water cut in a heterogeneous carbonate reservoir. The results could provide new insights and theoretical bases for the development of heterogeneous reservoirs.

2. Study Area

The oilfield studied in this paper is located in the north of the massif terrace of Arabia Peninsula. The geological strata is relatively stable, and there is almost no fault. The main rock lithology is carbonate and sandstone. Carbonate reservoirs usually have the characteristics of large lithologic variation, multiple reservoir types, and strong heterogeneity [14,15]. The oilfield studied in this paper is a typical carbonate reservoir with a heterogeneity of vertical permeability [16], and there is a great difference of geological structure and superimposed mode in different areas. According to the different values of permeability, we can divide formation into five types, which are tight type, low type, I type, II type, and III type. Furthermore the four superimposed modes were established in terms of the different ratios of permeability shown in Figure 1. The details of the parameters can be found in Table 1.
Figure 2 is the permeability distribution of four kinds of superimposed mode. We can conclude that the distribution of the superimposed modes A and D is near linear, and the distribution of the superimposed modes B and C tend toward a normal distribution. Thus, the superimposed modes B and C are more heterogeneous than the superimposed modes A and D [17].

3. Methods

The development of a carbonate oilfield has two exploitation stages. In the first stage, the driving force of the oil production is the original pressure and energy of the formation. The formation pressure gradually decreases with the development of the oil reservoir, which will lead to a significant decline of production. We can enhance oil recovery through the method of water-flooding. The oilfield studied in this paper began to use water-flooding in 2010.

3.1. Elastic Exploitation

The formation pressure gradually decreases with the development of the oil reservoir, and the change of formation pressure with time is satisfied as follows [18]:
{ 2 p r 2 + 1 r p r = 1 χ p t χ = K ϕ μ C t
where K is the permeability of the reservoir, m2; ϕ is the porosity, %; μ is the viscosity of fluid, Pa·s; and C t is the compressibility, MPa−1.
Define:
u = r 2 4 χ t
where u is dimensionless.
When the pressure of inner or outer boundary is constant, Equation (1) can be derived as follows:
{ u d 2 p d u 2 + d p d u ( 1 + u ) = 0 ( 0 < u < ) p ( u = 0 ) = p w p ( u = 450 ) = p e
where p w is the pressure under the conditions of u = 0 , and p e is the pressure under the conditions of u = 450 .
We can get the pressure distribution equation by solving Equation (2).
p ( u ) = p e Ei ( 1 , u ) ( 0.0365384234 p w 0.0365384234 p e )
Then
p ( r , t ) = p e Ei ( 1 , r 2 4 χ t ) ( 0.0365384234 p w 0.0365384234 p e )
We have analyzed the heterogeneity of the reservoir using variogram [19].
γ ( h d ) = 1 2 N ( h d ) i = 1 N ( h d ) [ Z ( k i ) Z ( k i + h d ) ] 2
K c = K ¯ ± 2 γ ( h d )
where K ¯ = i = 1 n h i k i / h t , h d is the lag distance, N ( h ) is the variation capacity, h i is the thickness of layer i , k i is the permeability of layer i , and h t is the total thickness of the reservoir; K c is the equivalent permeability.
According to the Darcy's law:
Q = K c A μ d p d r
Then we can get the production near wellbore:
Q = K c A μ p ( r w + Δ r , t ) p w Δ r
where r w is the shaft radius, A is the cross-section area of reservoir and Q is the flow rate.

3.2. Water-Flooding

The oil production decreased gradually with the decrease of formation pressure. In order to stabilize the oil production, we use the way of water-flooding to supplement the formation energy, and the flow resistance of each layer changes with the flood front. Take the case of superimposed mode A; we have analyzed the flow resistance of oil-water two phase flow in the process of water flooding, and obtained the water injection rate of each small layer [20].
R i = 1 2 π h i , j k i , j f i , 1 ln r i , 1 r w + + 1 2 π h i , j k i , j f i , n ln r e r n
where f i , n = k r o i , n μ o + k r w i , n μ w , n is the slug fluidity, i is the number of the laye and R i is resistance of layer i . This present model is the oil-water two phase flow model, so j = 1, 2.
The water injection rate of each layer can be expressed as:
Q i = p e p w R i = 1 R i 1 R 1 + 1 R 2 + 1 R i Q w
where Q w is total water injection rate.
In this paper, we divide the formation into five layers, and the flow resistance of each layer is as follows:
R 1 = 1 2 π h 1 k 1 f 1 , 1 ln r 1 , 1 r w + 1 2 π h 1 k 1 f 1 , 2 ln r e r 1 , 1
R 2 = 1 2 π h 2 k 2 f 2 , 1 ln r 2 , 1 r w + 1 2 π h 2 k 2 f 2 , 2 ln r e r 2 , 1
R 3 = 1 2 π h 3 k 3 f 3 , 1 ln r 3 , 1 r w + 1 2 π h 3 k 3 f 3 , 2 ln r e r 3 , 1
R 4 = 1 2 π h 4 k 4 f 4 , 1 ln r 4 , 1 r w + 1 2 π h 4 k 4 f 4 , 2 ln r e r 4 , 1
R 5 = 1 2 π h 5 k 5 f 5 , 1 ln r 5 , 1 r w + 1 2 π h 5 k 5 f 5 , 2 ln r e r 5 , 1
The relative permeability model in this paper is as follows:
K r o = k r o max k r o max × ( s o max s o ) n ( s o max s o r ) n
K r w = k r w max × ( s w s w r ) ( s w max s w r ) m
where s o + s w = 1 , s o is the oil saturation, s o max is the maximum oil saturation, s o r is the residual oil saturation, s w is the water saturation, s w max is the maximum water saturation, s w r is the irreducible water saturation, K r o is the oil relative permeability, and K r w is the water relative permeability, and m and n are characteristic constants.
f w (water cut ) can be expressed as follows considering the influence of capillary pressure and gravity:
f w = 1 + K r o μ r o 1 v t ( p c x Δ ρ g sin α ) 1 + μ w μ o K r o K r w

3.3. The Procedure of Simulation

Figure 3 is the simulation procedure. We can see the whole simulation process clearly through the flowchart. Firstly, the formation parameters are input into the model, then the variogram value of the whole reservoir is calculated. By comparing the parameters of well logging and seismic data, the error analysis is carried out. Then the optimal equivalent permeability is obtained through repeated calculation. Finally, the optimal model is used in the fitting and prediction. Table 2 is the parameter list for simulations.

4. Results and Discussion

In this paper, we compared the actual production data of a super-giant carbonate reservoir in the Middle East with the present model calculation, and the results showed that the present model could predict the productivity of wells in heterogeneous reservoirs more accurately and rapidly. Then, we analyzed the influence of the variogram value on the oil productivity, the relative permeability, and the water cut. Table 3 shows the production data of four wells with different superimposed modes.

4.1. The Comparison of Actual Production and the Present Model Calculation

Figure 4 displays the comparison of actual production data and the present model with different superimposed modes. The result shows that the present model can reflect the law of real production decline and predict the productivity of heterogeneous reservoirs easily and accurately. Due to more heterogeneity, the production declines of superimposed modes B and C are larger than the superimposed modes A and D under the same thickness.

4.2. Effect of Variogram Value on Production

Figure 5 shows the relationship between production and time with different variogram values. Figure 5 shows that the larger the varigoram value, the stronger the heterogeneity of reservoirs, and the faster decline of production. The stronger heterogeneity of the reservoir would lead to a reduction of the formation pressure more quickly.

4.3. Effect of Variogram Value on Water Cut

Figure 6 displays the relationship between the water cut and times with different variogram values. As shown in Figure 6, the water cut becomes larger with the increase in the variogram value. This means that the heterogeneity of the reservoir would be stronger with the increase of variogram value, which leads to faster water movement in high permeability layer than low permeability layer.

4.4. Effect of Variogram Value on Relative Permeability

Figure 7 shows relative permeability curves with different variogram values. The results show that the relative permeability of the oil phase is larger and the water phase is smaller with the decrease of the variogram value. The reason why is that the smaller variogram value is, the weaker the heterogeneity of reservoir, and then the water-free oil production period is longer, and the water cut rises slower.

4.5. The Advantages and Limitations of the Present Method

In this paper, the 2D variogram method is used to characterize the heterogeneity of carbonate reservoirs. Firstly, the 2D model can rapidly calculate the variation of the permeability in single direction, which can describe the heterogeneity of reservoir qualitatively. In addition, 2D model can provide more targeted understanding of the influence of the heterogeneity on development in single direction. However, the 2D model also has some limitations compared to 3D model. The impact of the reservoir development is various, and the 2D model cannot analyze the impact of all the factors, which leads to the inconsistency between theoretical analysis and engineering practice. We mainly study the effect of vertical heterogeneity on reservoir productivity in this paper, so the 2D model is more efficient.

5. Conclusions

The reservoirs were categorized into four superimposed modes using the actual logging data from a super-giant heterogeneous carbonate reservoir in the Middle East. A modified permeability formula in terms of the variogram method was presented to reflect the heterogeneity of reservoir. Based on the modified permeability formula, the models of oil production and water cut were established and the analytical solutions were obtained.
The results show that the present model can predict the productivity of a well with different heterogeneous layers more accurately and rapidly. Due to greater heterogeneity, the production declines if superimposed modes B and C are larger than the superimposed modes A and D under the same thickness. The larger the varigoram value, the stronger heterogeneity of reservoirs, the faster the decline of production owing to a quicker reduction of formation pressure. With the increase of the variogram value, the relative permeability of the oil phase is smaller and the water phase larger, and the water cut becomes larger.
This study has provided a quick and reasonable prediction model for heterogeneous reservoir development.

Acknowledgements

We gratefully acknowledge the PetroChina Science and Technology Major Project under Grant 2011E-2501 and the National Nature Science Foundations of China under Grant 51404024 for financial support.

Author Contributions

Chenji Wei proposed this topic and provided the actual data of oilfield. Hongqing Song and Yong Li supervised the work. Qi Zhang, Benbiao Song and Jiulong Wang were in charge of model establishment and analytical solution. Jiulong Wang organized the manuscript. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

K permeability of the reservoir
ϕ porosity
μ viscosity
C t compressibility
p w pressure under the conditions of u = 0
p e the pressure under the conditions of u = 450
h d lag distance
N ( h ) variation capacity
h i thickness of layer i
k i permeability of layer i
h t total thickness of the reservoir
K c equivalent permeability
A cross section area of reservoir
Q flow rate
r w shaft radius
Q w total water injection rate.
R i resistance of layer i
s o oil saturation
s o max maximum oil saturation
s o r residual oil saturation
s w water saturation
s w max maximum water saturation
s w r irreducible water saturation
K r o oil relative permeability
K r w water relative permeability
m characteristic constant
n characteristic constant
f w water cut

References

  1. Dang, T.Q.C.; Chen, Z.; Nguyen, T.B.N.; Bae, W. The potential of enhanced oil recovery by micellar/polymer flooding in heterogeneous reservoirs. Energy Sources Part A Recovery Util. Environ. Eff. 2014, 36, 1540–1554. [Google Scholar] [CrossRef]
  2. Hou, J.R.; Liu, Z.C.; Dong, M.Z.; Yue, X.A.; Yang, J.Z. Effect of viscosity of alkaline/surfactant/polymer (ASP) solution on enhanced oil recovery in heterogeneous reservoirs. J. Can. Pet. Technol. 2006, 45. [Google Scholar] [CrossRef]
  3. Ehrenberg, S.N. Whole core versus plugs: Scale dependence of porosity and permeability measurements in platform carbonates. AAPG Bull. 2007, 91, 835–846. [Google Scholar] [CrossRef]
  4. Yousef, A.A.; Al-Saleh, S.; Al-Jawfi, M.S. Improved/enhanced oil recovery from carbonate reservoirs by tuning injection water salinity and ionic content. In Proceedings of the Society of Petroleum Engineers (SPE) Improved Oil Recovery Symposium, Tulsa, OK, USA, 4–18 April 2012.
  5. Verma, V.K.; Datta, S. Flow in a channel filled by heterogeneous porous medium with a linear permeability variation. Spec. Top. Rev. Porous Media 2012, 3, 201–208. [Google Scholar] [CrossRef]
  6. Hulea, I.N.; Nicholls, C.A. Carbonate rock characterization and modeling: Capillary pressure and permeability in multimodal rocks—A look beyond sample specific heterogeneity. AAPG Bull. 2012, 96, 1627–1642. [Google Scholar] [CrossRef]
  7. Tidwell, V.C.; Wilson, J.L. Heterogeneity, permeability patterns, and permeability upscaling: Physical characterization of a block of massillon sandstone exhibiting nested scales of heterogeneity. SPE Reserv. Eval. Eng. 2000, 3, 283–291. [Google Scholar] [CrossRef]
  8. Chao, H.C.; Rajaram, H.; Illangasekare, T. Intermediate-scale experiments and numerical simulations of transport under radial flow in a two-dimensional heterogeneous porous medium. Water Resour. Res. 2000, 36, 2869–2884. [Google Scholar] [CrossRef]
  9. Correia, M.G.; Maschio, C.; Schiozer, D.J. Integration of multiscale carbonate reservoir heterogeneities in reservoir simulation. J. Pet. Sci. Eng. 2015, 131, 34–50. [Google Scholar] [CrossRef]
  10. Janković, I.; Fiori, A.; Dagan, G. Flow and transport in highly heterogeneous formations: 3. Numerical simulations and comparison with theoretical results. Water Resour. Res. 2003, 39. [Google Scholar] [CrossRef]
  11. Zhou, F.; Yao, G.; Wang, J. Effects of variogram characteristics of coal permeability on CBM production: A case study in Southeast Qinshui Basin, China. Energy Explor. Exploit. 2014, 32, 263–279. [Google Scholar] [CrossRef]
  12. Guadagnini, A.; Riva, M.; Franzetti, S. Flow to a well in a randomly heterogeneous aquifer with truncated multiscale variogram. Dev. Water Sci. 2002, 47, 1299–1306. [Google Scholar]
  13. Matonti, C.; Guglielmi, Y.; Viseur, S.; Bruna, P.O.; Borgomano, J.; Dahl, C.; Marié, L. Heterogeneities and diagenetic control on the spatial distribution of carbonate rocks acoustic properties at the outcrop scale. Tectonophysics 2014, 638, 94–111. [Google Scholar] [CrossRef]
  14. Al-Ameri, T.K.; Al-Khafaji, A.J.; Zumberge, J. Petroleum system analysis of the Mishrif reservoir in the Ratawi, Zubair, North and South Rumaila oil fields, Southern Iraq. GeoArabia 2009, 14, 91–108. [Google Scholar]
  15. Ziauddin, M.E.; Bize, E. The effect of pore scale heterogeneities on carbonate stimulation treatments. In Proceedings of the Society of Petroleum Engineers (SPE) Middle East Oil and Gas Show and Conference, Manama, Bahrain, 11–14 March 2007.
  16. Awadeesian, A.M.; Al-Jawed, S.N.A.; Saleh, A.H. Reservoir-scale sequence stratigraphy of Mishrif carbonates and implication to water injection strategy North Rumaila field case. Arab. J. Geosci. 2015, 8, 7025–7040. [Google Scholar] [CrossRef]
  17. Zech, A.; Zehner, B.; Kolditz, O.; Attinger, S. Impact of heterogeneous permeability distribution on the groundwater flow systems of a small sedimentary basin. J. Hydrol. 2016, 532, 90–101. [Google Scholar] [CrossRef]
  18. Kong, X.Y. The application of primitive function and Green function. In Advanced Seepage Mechanics; Press of the University of Science and Technology of China: Hefei, China, 1999. Chapter 4. (In Chinese) [Google Scholar]
  19. Ungaro, F.; Zasada, I.; Piorr, A. Mapping landscape services, spatial synergies and trade-offs. A case study using variogram models and geostatistical simulations in an agrarian landscape in North-East Germany. Ecol. Indic. 2014, 46, 367–378. [Google Scholar] [CrossRef]
  20. Song, H.Q.; Yao, Z.W.; Ming, W.; Kai, S.Y.; Tu, H.O.C. Complicated controllable displacement of equivalent filtrational resistance with several slugs. J. Univ. Sci. Technol. Beijing 2009, 175, 926–932. [Google Scholar]
Figure 1. The schematic of four superimposed modes: (a) superimposed mode A; (b) superimposed mode B; (c) superimposed mode C; (d) superimposed mode D.
Figure 1. The schematic of four superimposed modes: (a) superimposed mode A; (b) superimposed mode B; (c) superimposed mode C; (d) superimposed mode D.
Energies 10 00250 g001aEnergies 10 00250 g001b
Figure 2. Schematic of permeability distribution of different superimposed modes.
Figure 2. Schematic of permeability distribution of different superimposed modes.
Energies 10 00250 g002
Figure 3. The flow chart of simulation.
Figure 3. The flow chart of simulation.
Energies 10 00250 g003
Figure 4. Comparison of actual production data and the present model with different superimposed modes: (a) superimposed mode A; (b) superimposed mode B; (c) superimposed mode C; (d) superimposed mode D.
Figure 4. Comparison of actual production data and the present model with different superimposed modes: (a) superimposed mode A; (b) superimposed mode B; (c) superimposed mode C; (d) superimposed mode D.
Energies 10 00250 g004
Figure 5. Production with different variogram values.
Figure 5. Production with different variogram values.
Energies 10 00250 g005
Figure 6. Water cuts with different variogram values.
Figure 6. Water cuts with different variogram values.
Energies 10 00250 g006
Figure 7. Relative permeability with different variogram values.
Figure 7. Relative permeability with different variogram values.
Energies 10 00250 g007
Table 1. The permeability distribution of the different superimposed modes.
Table 1. The permeability distribution of the different superimposed modes.
Reservoir TypesPermeability (mD)Percentage (%)
Superimposed Mode ASuperimposed Mode BSuperimposed Mode CSuperimposed Mode D
I Type100–150552053
II Type70–10025252315
III Type40–708304430
Low10–404101840
Tight5–108151012
Table 2. Parameter list for simulations.
Table 2. Parameter list for simulations.
ParameterValueParameterValue
n (dimensionless)2.3 p e (MPa)6
m (dimensionless)1.4 p e (MPa)26
ϕ 0.2 r e (m)500
μ o (Pa·s)2.7 × 10−3 r w (m)0.1
μ w (Pa·s)1 × 10−3 h t (m)100
C t (MPa−1)1 × 10−7
Table 3. The production data of four superimposed modes/wells.
Table 3. The production data of four superimposed modes/wells.
ModelSuperimposed Mode A(t/d)Superimposed Mode B(t/d)Superimposed Mode C(t/d)Superimposed Mode D(t/d)
Time (d)
01000.00714.29553.00382.71
120857.14622.14511.57335.29
240714.29642.86465.14295.57
360636.00601.43437.86271.86
480571.43579.00424.86248.86
600500.00547.57393.57239.57
720442.86488.43378.71234.71
840422.57471.43356.57191.71
960404.43424.00304.71147.57
1080382.43396.29295.00147.43
1200352.43376.14276.43138.29
1320322.57368.71271.43133.29
1440285.71361.86281.14138.29
1560165.29330.00266.71133.29
1680163.29313.29244.29124.43

Share and Cite

MDPI and ACS Style

Wei, C.; Song, H.; Li, Y.; Zhang, Q.; Song, B.; Wang, J. Production Characteristics with Different Superimposed Modes Using Variogram: A Case Study of a Super-Giant Carbonate Reservoir in the Middle East. Energies 2017, 10, 250. https://0-doi-org.brum.beds.ac.uk/10.3390/en10020250

AMA Style

Wei C, Song H, Li Y, Zhang Q, Song B, Wang J. Production Characteristics with Different Superimposed Modes Using Variogram: A Case Study of a Super-Giant Carbonate Reservoir in the Middle East. Energies. 2017; 10(2):250. https://0-doi-org.brum.beds.ac.uk/10.3390/en10020250

Chicago/Turabian Style

Wei, Chenji, Hongqing Song, Yong Li, Qi Zhang, Benbiao Song, and Jiulong Wang. 2017. "Production Characteristics with Different Superimposed Modes Using Variogram: A Case Study of a Super-Giant Carbonate Reservoir in the Middle East" Energies 10, no. 2: 250. https://0-doi-org.brum.beds.ac.uk/10.3390/en10020250

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop