SciELO - Scientific Electronic Library Online

 
vol.22 número5Enhancing the electrochemical performance of next-generation 5 V class Li-ion batteries using LiPF6/LiBOB mixed salt electrolyteGraphene nano particles to improve lightning dissipation on transmission lines índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Journal of applied research and technology

versión On-line ISSN 2448-6736versión impresa ISSN 1665-6423

J. appl. res. technol vol.22 no.5 Ciudad de México oct. 2024  Epub 11-Ago-2025

https://doi.org/10.22201/icat.24486736e.2024.22.5.2394 

Articles

Study of recovery factor and CO2 storage behavior in a reservoir using numerical simulation

aUniversidad Estatal Península de Santa Elena, Santa Elena, Ecuador


Abstract:

In this research, numerical simulation was used to analyze the recovery factor during the CO2 displacement and storage phases in an Ecuadorian field. CO2 continuous injection and WAG methods were used. The methodology included field data collection, creation of a static reservoir model, configuration of the PVT model, construction of a dynamic model in GEM-CMG, and evaluation of CO2 injection scenarios. After a screening analysis, the optimal scenario was identified and the behavior of the recovery factor over time was studied. The simulation results indicate that the recovery factor with continuous injection was 28%, while with the WAG method, it reached 31%. It was concluded that the most efficient method was the WAG. Regarding storage, it was observed that continuous injection managed to store 83% of the CO2 in the reservoir. This is essential since it helps mitigate global warming by reducing CO2 emissions into the atmosphere.

Keywords: Recovery factor; CO2 storage; CMG; WAG

1. Introduction

In Ecuador, oil production constitutes one of the main economic sources (Ordoñez Moran & Nuñez Lapo, 2019). However, after several years of production, the energy of some fields has decreased (Rivadeneira & Baby, 2004). To recover this energy, enhanced recovery (EOR) methods can be used, such as CO2 injection.

The interaction between CO2 and oil in the reservoir leads to a phenomenon known as miscibility. This means that CO2 dissolves in the oil, resulting in a reduction in the viscosity of the oil, making it easier for the oil to move towards producing wells (Abdullah & Hasan, 2021). Miscibility in this context implies that CO2 and oil form a homogeneous phase, which maximizes the contact between both fluids and increases recovery efficiency (Verma, 2015 This miscibility can be generated by first contact or by multiple contacts (Castro, 2013).

The oil fields in Ecuador have a production history of more than 50 years, which classifies them mostly as mature fields (Chipantashi-Aneloa, 2021). The implementation of the CO2 injection technique as an Enhanced Oil Recovery (EOR) method is essential in these fields since it contributes to maintaining oil production and recovering reserves that are still exploitable (Farajzadeh et al., 2022). Furthermore, this technique allows a significant amount of carbon dioxide to be stored in the deposit, reducing its release into the atmosphere, and mitigating the effects of environmental pollution caused by this gas (Godec et al., 2011).

1.1. Generalities of CO2 injection

CO2 is injected into the oil reservoir to provide energy. CO2 dissolves in the oil, reducing its viscosity and moving to the producing wells. Currently, improved hydrocarbon recovery is the most viable economic option in the CO2 capture and storage process because it increases production and improves the cost and income balance in said processes (Alquicira Balderas, 2018).

When CO2 is injected into the reservoir, a physical and chemical interaction between the reservoir rock and the existing hydrocarbon fluid arises (Moghadasi et al., 2018). These interactions explain the mechanism of recovery+ of the remaining oil. These mechanisms are classified into:

  • Increase in oil volume.

  • Decrease in the density of oil and water.

  • Decrease in oil viscosity.

  • Decrease in interfacial tension between rock and fluid.

A detailed description of the methods used in this research, as well as the criteria on which their choice was based, is now provided.

1.2. Continuous CO 2 injection method

In this method, CO2 is injected continuously until the values of the produced gas-oil ratio reach levels so high that the project is not financially viable (Santamaría, 2018). The multi-contact process, condensation and combination mechanisms are used to generate miscibility between gas and oil (Safaei et al., 2023).

1.3. Alternate water and CO 2 injection method (WAG)

WAG injection is a method that is applied in almost all gas injection processes, it prevents the formation of CO2 channels through the oil phase, increasing the macroscopic displacement efficiency (Santamaría, 2018). Figure 1 presents several gas injection techniques, the total volume of gas to be injected and the water-gas relationship and frequency; each of these variables must be calculated for the WAG process.

Figure 1 WAG alternating water gas flooding methods, (Ameli et al., 2023). 

Furthermore, WAG injection is controlled by Valladares and Steven (2018):

  • Ratio of the volume of CO2 injected over the volume of water injected.

  • Pothole size.

  • Pothole number.

  • Advantages of WAG injection.

  • Reduce the relative permeability to CO2 and therefore its mobility.

  • Covers the CO2 requirement in each instance.

1.4. Selection mechanisms for each method

The selection criteria for CO2 injection methods are not only determined by one or two reservoir characteristics, but rather the best approach is determined based on viscosity, permeability, API, reservoir temperature and reservoir depth (Ameli et al., 2023).

Table 1 and Table 2 present the appropriate selection standards for continuous injection and WAG projects respectively.

Table 1 Selection criteria according to the properties of the deposit, (Ameli et al., 2023

Property Characteristic
Reservoir depth Deeper wells imply higher drilling costs and operating expenses because the reservoir temperature increases with depth.
Temperature Critical temperature of CO2 is 31°C, where CO2 operates as vapor under these conditions. At temperatures above this level, the density of CO2 rises with pressure.
Pressure Various degrees of miscibility can be achieved at intermediate and high pressure.
Porosity and permeability Some reservoirs have varying levels of porosity and commonly range between 11% and 30%.

Table 2 Selection criteria for WAG projects, (Khoshsima et al., 2023). 

Parameters Values
Reservoir Thickness (ft) <100
Temperature (°F) -
Average permeability (MD) <100
Previous production method Water injection
Type of training -
Depth -
Fluid API 29-45
Viscosity (cP) <2
Viscosity ratio sep-30

2. Methodology

The information required to perform the simulation in static and dynamic conditions is obtained from various data sources. In relation to information processing, the GMSH simulator will be implemented to generate the mesh and the CMG software to define the model under dynamic conditions and execute the CO2 injection methods. In addition, the GEM, WinProp and Results sections will be implemented (CMG GEM, n.d.).

Next, we proceed to describe in detail the process carried out in each of the sections used.

2.1. Static reservoir modeling

To generate the geological unit, a portion of the sand T under study is digitized, taking into consideration the location of the wells it contains. This process is carried out using the GMSH software. It is important to highlight that, to carry out this procedure, you must have an isopach map of the reservoir to be simulated.

After importing the image into GMSH, the "points" tool is used to generate values in a Cartesian plane for both "x" and "y". In this way, points are added to build all the layers immersed in the sand (Figure 2).

Figure 2 History of points generated in GMSH. 

And for the characterization of the reservoir fluids, the chromatography analysis of the fluids coming from the sand under study was implemented, using the WinProp-CMG (Figure 3).

Figure 3 WinProp component login. 

After having estimated the necessary calculations and laboratory tests of the reservoir fluid composition in the WinProp section, it is required to export the PVT model.

-Initialization of the CMG simulator.

When starting the CMG program, the GEM option is chosen, a widely recognized tool worldwide for modeling compositional, chemical, and unconventional reservoirs. Where the system of units is selected and the start date of the model is set.

-Description of the reservoir, creation of the grid.

To create the reservoir mesh (Figure 4), the mesh was configured in spatial units "i, j, k", these being 60, 106 and 8, respectively. Each cell "i, j" has a length of 100 feet. It is worth mentioning that the top of the formation is 9,735 feet deep and has a thickness of 50 feet.

Figure 4 Contour map with Orthogonal Corner Point grid. Own preparation, based on the CMG simulator. 

-Rock properties.

The general characteristics of the rock are detailed in the following data sheet (Figure 5).

Figure 5 Rock properties input window. 

-Components.

In this section, the PVT data corresponding to the fluid is entered (Table 3).

Table 3 Fluid properties. 

PROPERTIES VALUE
Oil volumetric factor (BY/BN) 1,48
Reservoir temperature (°F) 196
Reservoir pressure (Psi) 4226
°API 30,8
GOR (SCF/STB) 180
Density of oil (gr/cm3) 0,82
Bubble pressure (Psi) 835

-Rock-fluid interaction.

In this section the data for the creation of the relative permeability curves of the oil-water and oil-gas systems are entered (Table 4 and Table 5).

Table 4 Relative permeabilities of the oil-water system. 

Sw Krw Kro
0.228 0 0.7943
0.248 0.0015 0.7322
0.315 0.0036 0.6526
0.3423 0.0058 0.5456
0.4328 0.0185 0.3632
0.5156 0.0419 0.2261
0.5945 0.0854 0.0853
0.6358 0.1142 0.0632
0.6645 0.1525 0.052
0.7472 0.256 0.034
0.7924 0.327 0

Table 5 Relative permeabilities of the oil-gas system. 

Sg Krg Kro
0.0005 0 0.7943
0.06 0 0.6534
0.0858 0.002 0.5493
0.1842 0.013 0.4325
0.2683 0.0353 0.2864
0.3676 0.052 0.184
0.4895 0.156 0.08532
0.5532 0.283 0.038
0.6428 0.3623 0.0124
0.6634 0.3856 0
0.8335 0.999 0

-Initial reservoir conditions.

This section details the values of reservoir pressure, reference depth (average depth of the layer of interest) and the level of contact between water and oil (see Figure 6).

Figure 6 Initial reservoir parameters input. 

2.2. Dynamic modeling of the sand to be studied.

-Wells and recurrent data.

Initially, the seven production wells and one water injection well existing in the field will be built. The producing wells are called: SSF-032, SSF-013, SSF-017, SSF-028, SSF-037, SSF-039 and SSF-047, while the injection well will be identified as SSF-W-I.

In simulation, two main factors are taken into account for producing wells. The first is the maximum fluid flow rate at the surface, set at 150 barrels per day. The second fundamental factor is the minimum bottomhole pressure (BHP), set at 50 Psi.

The conditions set for the water injection well are as follows: Maximum injection rate of 400 barrels per day and a maximum bottom hole pressure (BHP) of 100 Psi.

Construction of injection methods

For this work, "February 1, 2020" is considered as the starting date for the execution of the injection methods to be evaluated; continuous injection of CO2 and alternating injection of water and CO2 (WAG). This process is executed for 8 years, that is, until "February 1, 2028". Furthermore, in both cases the same volume of CO2 will be injected, which is equivalent to 15% of the pore volume, i.e., 621690000 ft3.

-Continuous injection method

Two new CO2 injection wells were created in the continuous CO2 injection model. The new wells will be named "SSF-CO2-I1" and "SSF-CO2-I2". The date of creation of the new wells is "February 1, 2020".

The operational parameters established for the CO2 injector wells are as follows: The largest amount of carbon dioxide that can be injected per day is 42600 cubic feet, and the maximum downhole pressure reaches 1500 psi.

-WAG injection method.

As in the previous method, the location of the wells will be the same. The specific conditions considered for this method are presented in Table 6.

Table 6 Process operating conditions. 

WAG Injection
Year of start of WAG injection 2020
Plug size 6% VP
WAG Relationship 0.5
Total volume of gas injected 1 VP
Frequency of cycles 1 year
Water injection rate per day at each well 6100 bbl
Gas injection rate per day in each well 68100 ft3

Finally, the reservoir mesh with the wells created is shown (Figure 7). In this way, we proceed to start the simulation of our methods, with the purpose of analyzing the behavior of the reservoir during the established time.

Figure 7 Reservoir grid with 7 producing wells, 1 water injector well and 2 CO2 injector wells. Own preparation, based on the CMG simulator. 

3. Results

3.1. Recovery factor

In 2020, the recovery factor is less than 15%. However, as of 01/06/2022, there has been an increase in the recovery factor due to the implementation of CO2 injection. In the continuous CO2 injection method, a recovery factor of 28% was achieved. On the other hand, through the application of the WAG method, a significant recovery factor was achieved, reaching 31% as shown in Figure 8.

Figure 8 Recovery Factors of the Methods. 

3.2. CO 2 storage

In Figure 9 and Figure 10, it can be seen that the moles of CO2 trapped in the reservoir begin to increase significantly from 2021 onwards. This increase is due to the start of the injection this year.

Figure 9 CO2 Storage from Continuous Injection. 

Figure 10 CO2 Storage from WAG Injection. 

With the continuous CO2 injection method, a decrease is observed until 08/09/2022 where it reaches a value of 4.6x10e6 moles of CO2. Subsequently, the value of the moles of CO2 trapped in the reservoir begins to increase continuously, reaching a value of 2.3x10e7 moles in 2028.

Using the WAG method in Figure 10, it can be seen that the volume of trapped gas increases and decreases very sharply due to the CO2 and water injection cycles. At the end of the injection process in 2028, an amount of 9.0x10e6 moles is stored in the reservoir.

In this way, it is shown that the displacement or sweeping of oil is directly related to the amount of CO2 that is stored in the reservoir.

4. Conclusions

CO2 injection was implemented with the purpose of taking advantage of the associated gas that is burned in the burners in the field. This technique seeks to increase the productivity of wells by applying methods of continuous injection of CO2 and WAG at high pressures. The main objective is to achieve high miscibility between CO2 and oil, resulting in viscosity reduction, thus improving oil mobility, and increasing reservoir sweeping efficiency. These processes led to increased displacement of residual oil and ultimately more efficient oil recovery.

Based on the simulation carried out, it was determined that the recovery factors of the injection methods are as follows: in continuous injection it is 28%, while in WAG injection it amounts to 31%. In this way, it was possible to demonstrate that the most efficient method for the deposit under study is WAG. In addition, a cumulative production of 2.54 MMbbl was obtained for continuous injection during the 8 years of simulation and for WAG injection a cumulative production of 2.72 MMbbl was reached, which justifies the increase in the recovery factor until 2028.

The CO2 storage obtained until 2028, by the WAG injection method is 1.23x10e7 moles of CO2, which represents 36% and the remaining 64% corresponds to the recovered on the surface and with the continuous injection method of 5.73x10e7 moles of CO2, which represents 83% stored and 17% recovered on the surface. It should be noted that the continuous injection method proved to be more effective in terms of storage, managing to store a greater amount of CO2 compared to the WAG method.

Acknowledgements

We thank God for the opportunity He has given us to develop and execute this work.

References

Abdullah, N., & Hasan, N. (2021). Effects of miscible CO2 injection on production recovery. Journal of Petroleum Exploration and Production. https://doi.org/10.1007/s13202-021-01223-0 [ Links ]

Alquicira Balderas, M. A. (2018). Potencial del CO2 como método de recuperación mejorada de hidrocarburos en México. Repositorio Facultad de Ingeniería. http://www.ptolomeo.unam.mx:8080/xmlui/handle/132.248.52.100/16038Links ]

Ameli, F., Moghaddam, S., & Shahmarvand, S. (2023). Introduction to gas flooding technique: fundamentals and applications. Gas Injection Methods: A Volume in Enhanced Oil Recovery Series, 1-43. https://doi.org/10.1016/B978-0-12-822302-4.00008-9 [ Links ]

Castro Castañeda, L. F. (2013).Inyección de CO2 como método de recuperación mejorada[Doctoral dissertation, Tesis para obtener el título de Ingeniero Petrolero. Facultad de Ingeniería del la Universidad Autónoma de México. Ciudad Universitaria, México, D.F]. https://repositorio.unam.mx/contenidos/254921Links ]

Chipantashi-Aneloa, E. P. (2021). Potencial del fracturamiento hidráulico en el Ecuador: Estudio del caso Campo Guanta-Dureno .FIGEMPA: Investigación Y Desarrollo,12(2), 37-51. https://doi.org/10.29166/revfig.v12i2.2435 [ Links ]

Christian Santamaría Guzmán. (2018). Recuperación Mejorada de Petróleo Mediante la Inyección de Dióxido de Carbono Antropogénico. https://petroquimex.com/recuperacion-mejorada-de-petroleo-mediante-la-inyeccion-de-dioxido-de-carbono-antropogenico/Links ]

Santamaría, G., C. (2019). Recuperación Mejorada de Petróleo Mediante la Inyección de Dióxido de Carbono Antropogénico.La revista de la industría Energética. https://petroquimex.com/recuperacion-mejorada-de-petroleo-mediante-la-inyeccion-de-dioxido-de-carbono-antropogenico/Links ]

CMG GEM | Compositional and Unconventional Reservoir Simulation. (n.d.). https://www.cmgl.ca/gemLinks ]

Valladares, H., & Steven, J. (2018).Optimización de un proceso de inyección WAG mediante simulación numérica para su aplicación en yacimientos de Venezuela(Doctoral dissertation). http://hdl.handle.net/10872/18051Links ]

Farajzadeh, R., Glasbergen, G., Karpan, V., Mjeni, R., Boersma, D. M., Eftekhari, A. A., Casquera Garcia, A., & Bruining, J. (2022). Improved oil recovery techniques and their role in energy efficiency and reducing CO2 footprint of oil production. Journal of Cleaner Production, 369, 133308. https://doi.org/10.1016/J.JCLEPRO.2022.133308 [ Links ]

Godec, M., Kuuskraa, V., Van Leeuwen, T., Melzer, L. S., & Wildgust, N. (2011). CO2 storage in depleted oil fields: The worldwide potential for carbon dioxide enhanced oil recovery. Energy Procedia, 4, 2162-2169. https://doi.org/10.1016/J.EGYPRO.2011.02.102 [ Links ]

Khoshsima, A., Sedighi, M., & Mohammadi, M. (2023). Enhanced oil recovery by wáter alternating gas injection. Gas Injection Methods: A Volume in Enhanced Oil Recovery Series, 295-316. https://doi.org/10.1016/B978-0-12-822302-4.00001-6 [ Links ]

Moghadasi, R., Rostami, A., & Hemmati-Sarapardeh, A. (2018). Chapter Three - Enhanced Oil Recovery Using CO2 . Fundamentals of Enhanced Oil and Gas Recovery from Conventional and Unconventional Reservoirs, 61-99. https://doi.org/10.1016/B978-0-12-813027-8.00003-5 [ Links ]

Ordoñez Moran, A., & Nuñez Lapo, M. D. L. A. (2019). Análisis del sector petrolero y su incidencia en los aspectos socioeconómicos del Ecuador en el período 1980-2012.Revista Universidad y Sociedad,11(5), 405-409. [ Links ]

Rivadeneira, M., & Baby, P. (2004). Características geológicas generales de los principales campos petroleros de petroproducción. La Cuenca Oriente: Geología y Petróleo, 229-295. https://doi.org/10.4000/BOOKS.IFEA.3020 [ Links ]

Safaei, A., Bagheri, M., Kazemzadeh, Y., & Riazi, M. (2023). Carbon dioxide miscible flooding. Gas Injection Methods: A Volume in Enhanced Oil Recovery Series, 133-169. https://doi.org/10.1016/B978-0-12-822302-4.00003-X [ Links ]

Verma, M. K. (2015).Fundamentals of Carbon Dioxide-Enhanced Oil Recovery (CO 2-EOR)-A Supporting Document of the Assessment Methodology for Hydrocarbon Recovery Using CO 2-EOR Associated with Carbon Sequestration. (No. 2015-1071). US Geological Survey. https://doi.org/10.3133/ofr20151071 [ Links ]

Funding

The authors received no specific funding for this work.

Received: December 06, 2023; Accepted: July 23, 2024; Revised: October 31, 2024

Corresponding author. E-mail address: douglas.tomalasuarez@upse.edu.ec (D. Tomalá).

Conflict of interest

The author has no conflict of interest to declare.

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License