Optimizing Geothermal Energy Extraction in CO 2 Plume Geothermal Systems

: CPG (CO 2 Plume Geothermal) has recently emerged as a promising technology to combine the extraction of geothermal energy with underground CO 2 storage, thus achieving double positive results. The idea is to inject CO 2 in its supercritical form to create a plume and replace the reservoir brine, which is continuously circulated to transfer heat from the reservoir to surface facilities. Apart from the positive aspects of this technology, including the reduced energy to inject and lift the working ﬂuid as well as its enhanced mobility in the reservoir and reduced environmental footprint, there are also negative issues that must be handled by adequately studying the geological ﬁeld/reservoir and appropriately designing the production system. In this work, we present a ﬁnite volume numerical simulation that can study a geothermal reservoir from its geological origin to the dynamic simulation of CO 2 injection and estimate the geothermal energy extraction. It is shown that the system performance is strongly related to the selected schedule, and optimizing it in conjunction with the related cost is of the utmost importance for the Final Investment Decision to be taken and for the viability of such multipurpose projects under a sustainable future.


Introduction
Geothermal power plants utilize the temperature disparity between the Earth's heated subsurface rock and the cooler surface to produce electricity.These systems convey thermal energy from beneath the ground to the surface using a working fluid, which subsequently undergoes a partial transformation of thermal energy to electricity within a power plant.The cooled working fluid is typically reinjected into a subsurface reservoir to establish hydraulic sustainability.Traditional geothermal energy technologies employ hot brine as the working fluid, and, as a result, power plants are often situated in regions with active tectonic or volcanic activity where the temperature gradient is exceptionally high [1].Such fields are characterized as high enthalpy ones and are utilized for power generation.In Greece, for example, they are directly related to the well-known subduction of the African lithospheric plate beneath the Aegean microplate and the subsequent development of the South Aegean Active Volcanic Arc (SAAVA).High enthalpy geothermal fields in Greece appear in the Cyclades group of Islands to be relevant to the aforementioned arc volcanic activity (e.g., the "proven" Milos and Nisyros fields) and in the wider area surrounding Lesvos island (older location of the arc that currently moved southwards) [2], while within the same area, shallow-depth, low-temperature geothermal fields occur as well, e.g., that of SW Santorini Island, which is characterized as "probable" [3].Nonetheless, such distinctive thermal resources are limited in terms of their size and location, prompting the need for innovative technologies to exploit the vast amount of thermal energy present in the Earth's crust.
Enhanced Geothermal Systems (EGS) have been suggested as a means to expand access to geothermal resources by artificially constructing reservoirs in areas that lack suitable conditions [4].These systems typically involve the injection of cold pressurized water for the hydraulic fracturing of a subsurface formation.The injected water receives heat and is pumped back to the surface.Carbon dioxide (CO 2 ) has been proposed as an alternative working fluid due to its non-toxic, abundant, and non-flammability properties [5].However, the challenges linked to using CO 2 in EGS, such as its reduced density at reservoir conditions and environmental impact in the case of a leakage, necessitate further research and development.Addressing the risks associated with employing CO 2 in EGS, CO 2 Plume Geothermal (CPG) energy systems have been recently introduced and even developed [6][7][8].These systems use CO 2 as the primary subsurface working fluid in naturally permeable sedimentary basins or EGS, generating a large-scale CO 2 plume.Furthermore, a buoyancy-driven thermosiphon can be established by exploiting CO 2 density variations between the injection and production wells.By performing this, the requirement for expensive pumping, which is commonly associated with conventional hydrothermal setups, can be eliminated.
Sedimentary basins appropriate for CPG systems are found globally [9].These basins often comprise aquifers with high salinity, rendering them unsuitable for drinking or industrial water resources.They may also encompass partially depleted oil and gas fields utilized for Enhanced Oil Recovery (EOR) operations [10].In a CPG system, the buoyant CO 2 must be contained by very low-permeability or fully impervious caprock beds that cover the permeable reservoir.Additionally, CPG systems can be integrated with CO 2 capture and storage (CCS) sites [11][12][13], enabling the simultaneous generation of electricity and heat while sequestering CO 2 .This integration ensures reservoir stability by alleviating overpressurization in relation to standalone CCS operations, which could provoke humaninduced seismicity and CO 2 leakage.Moreover, a combined CPG-CCS system can boost the profitability of CCS, thus promoting global initiatives to combat climate change.
In this work, a deep, high-permeable, deep saline aquifer is selected for the development of the joint CPG-CCS system.The dynamics of injection involve two-dimensional variables.Namely, these include the mobility ratio and the gravity number.Via construction, the gravity number controls the strength of buoyancy (i.e., with or without buoyancy) and depends on petrophysical parameters.On the other hand, the mobility ratio depends on the CO 2 and brine fluid properties as well as the capillary forces between the two fluids and the wettability of the rock.As such, it is understood that the success of the injection depends on the flow regime that the injection scenario falls into.Nevertheless, various analytical solutions, verified by numerical simulation, have been proposed to handle related issues such as cap rock uplift [14], plume pressure buildup [15], and flow regime analysis [16].However, these solutions do not provide any information about the porous formation under inclination or under the action of thermal gradients, thus justifying the use of numerical modeling.

Model Description
The aquifer under consideration possesses distinctive characteristics.Situated at a depth of 3.5 km, it lacks hydraulic connectivity with adjacent formations both laterally and vertically, making it an isolated system with Neumann no-flow boundary conditions applied.Initially, it is slightly underpressurized with an average pressure of 3800 psi, which is slightly below what would be projected based on the water and rock formation gradient.Due to magmatic influences, the initial temperature is set at 500 • F. It exhibits a substantial inclination of 3 degrees, with an average thickness of 12.7 m and an areal extent reaching 80 square kilometers.To maintain simplicity in our analysis, we assumed the reservoir to be both homogeneous and isotropic in the xy plane.The elevated temperature within the aquifer exceeded what would be anticipated from a typical thermal gradient, suggesting an underlying magmatic influence.Thanks to its high and isotropic permeability, pressure is nearly uniformly distributed across all cells.Monitoring pressure within the aquifer ensures that the fracturing limit, which could result in unintended rock fractures, is not surpassed.Our calculations are based on the numerical solution of the mass, momentum, and energy differential equations governing fluid flow within the porous medium.
The selected aquifer for this simulation mirrors PVT properties, which are commonly found in the SAAVA area.The simulator utilized in this study [17] is a 3D black oil model that integrates fluid flow and heat transfer.Given our dual focus on CCS and geothermal engineering, technical units are consistently employed throughout the text, with conversion factors detailed in Table 1 to aid reader comprehension.The characteristics of the aquifer are presented in Table 2.

Schedule
Since CO 2 is lighter than water even in its supercritical phase, the inclination of the system exhibited in Figure 1 can be effectively utilized.The proposed methodology involves extracting water from the bottom of the aquifer to create space for injecting CO 2 at the top for two years (phase A).The density disparity between the two fluids generates a CO 2 plume that expands rapidly within the upper layer, allowing it to reach deeper production wells.However, due to the aquifer inclination and the limited vertical permeability compared to the horizontal one (kz = 0.1 × kx), the migration of the CO 2 plume experiences a slight deceleration.The Voidage Replacement Ratio (VRR), defined as the ratio between the volume of the injected fluid and the volume of the produced fluid, exhibits variations during the process due to the influence of the downhole temperature on the density of the injected fluid and the occurrence of CO2 breakthrough toward the latter stages.On average, the VRR maintains a value of 91%, resulting in a minor decrease in pressure.exhibits variations during the process due to the influence of the downhole temperature on the density of the injected fluid and the occurrence of CO2 breakthrough toward the latter stages.On average, the VRR maintains a value of 91%, resulting in a minor decrease in pressure.Once a substantial amount of water has been produced and the simultaneous injection of CO2 has taken place, Wells 3 and 4, which were formerly used to withdraw brine, are closed to prevent the concurrent production of CO2 after CO2 breakthrough has occurred.Wells 1 and 2 continue operating as before (i.e., injecting CO2) but with an increased injection rate for another thirteen years (phase B) to model the arrival of additional emitters.This phase is characterized by a significant pressure buildup within the aquifer.It is crucial to monitor the maximum pressure to ensure it remains below the safe limit of the fracture pressure of the caprock.After the maximum quantity of CO2 has been securely stored, the CPG phase is initiated.Well 3 is converted into a production well, while Well 4 maintains its injection activity at a slightly higher rate than the production rate of Well 3.This adjustment is made to account for the presence of a relatively small-scale CO2 emitter, allowing for its safe storage throughout the CPG phase without risking uncontrolled pressure escalation.Consequently, a functional CPG system is set in motion, the sustainability of which is evaluated for the next ten years (phase C).The control scheme is expressed in detail in Tables 3 and 4.  Once a substantial amount of water has been produced and the simultaneous injection of CO 2 has taken place, Wells 3 and 4, which were formerly used to withdraw brine, are closed to prevent the concurrent production of CO 2 after CO 2 breakthrough has occurred.Wells 1 and 2 continue operating as before (i.e., injecting CO 2 ) but with an increased injection rate for another thirteen years (phase B) to model the arrival of additional emitters.This phase is characterized by a significant pressure buildup within the aquifer.It is crucial to monitor the maximum pressure to ensure it remains below the safe limit of the fracture pressure of the caprock.After the maximum quantity of CO 2 has been securely stored, the CPG phase is initiated.Well 3 is converted into a production well, while Well 4 maintains its injection activity at a slightly higher rate than the production rate of Well 3.This adjustment is made to account for the presence of a relatively small-scale CO 2 emitter, allowing for its safe storage throughout the CPG phase without risking uncontrolled pressure escalation.Consequently, a functional CPG system is set in motion, the sustainability of which is evaluated for the next ten years (phase C).The control scheme is expressed in detail in Tables 3 and 4.

Discussion/Results
The simulations are conducted using the Reveal simulator, a software developed by Petroleum Experts (IPM) [17].Throughout the simulations, the average pressure within the reservoir remains safely below the fracturing pressure.It is important to note that, in contrast to oil production, where maximizing profits often drives greedy optimization schedules, CCS is inherently unprofitable.Therefore, there is no necessity to adopt a greedy schedule that could push the aquifer's pressure close to its fracture limit.
The most important result obtained is the total mass of CO 2 sequestered, as depicted in Figure 2, and the geothermal power output that can be extracted from CO 2 produced before it is recycled, as shown in Equation ( 1) where h wh prod , h wh inj are the specific enthalpy of CO 2 at the wellhead conditions in the production (Well_3) and injection (Well_4) wells, respectively, and .m is the mass flow rate.

Discussion/Results
The simulations are conducted using the Reveal simulator, a software developed by Petroleum Experts (IPM) [17].Throughout the simulations, the average pressure within the reservoir remains safely below the fracturing pressure.It is important to note that, in contrast to oil production, where maximizing profits often drives greedy optimization schedules, CCS is inherently unprofitable.Therefore, there is no necessity to adopt a greedy schedule that could push the aquifer's pressure close to its fracture limit.
The most important result obtained is the total mass of CO2 sequestered, as depicted in Figure 2, and the geothermal power output that can be extracted from CO2 produced before it is recycled, as shown in Equation ( 1) where ℎ  ℎ , ℎ  ℎ are the specific enthalpy of CO2 at the wellhead conditions in the production (Well_3) and injection (Well_4) wells, respectively, and ̇ is the mass flow rate.After 25 years of control, the total mass of CO2 that is sequestered within the reservoir is estimated at 36.3 MMtonnes.To calculate the power output, the produced CO2 mass rate is directly obtained from the simulator and equals ̇ = 6.58Mtonnes/day.CoolProps [18] software was utilized to estimate CO2-specific enthalpy by inputting the temperature and pressure values, as shown in Figure 3.The average specific to the production well is calculated at 238.7 BTU/lbm and at the injection well at 138.6 BTU/lbm.Subsequently, power output can be calculated in a straightforward manner using Equation ( 1) and the mass flow rate from Table 4 after appropriate unit conversions.It is noteworthy that up to After 25 years of control, the total mass of CO 2 that is sequestered within the reservoir is estimated at 36.3 MMtonnes.To calculate the power output, the produced CO 2 mass rate is directly obtained from the simulator and equals .m = 6.58Mtonnes/day.CoolProps [18] software was utilized to estimate CO 2 -specific enthalpy by inputting the temperature and pressure values, as shown in Figure 3.The average specific to the production well is calculated at 238.7 BTU/lbm and at the injection well at 138.6 BTU/lbm.Subsequently, power output can be calculated in a straightforward manner using Equation (1) and the mass flow rate from Table 4 after appropriate unit conversions.It is noteworthy that up to this point, Well_3 (producer) was an injection well for over a decade.Subsequently, the bottomhole temperature is initially only 160 • F and rapidly converges the desired range due to the rock heating it and hot fluids arriving at the bottomhole.The temperature observed during the later stages is considered representative of enthalpy calculations.Additionally, due to the ongoing injection of CO 2 at a slightly higher rate than the production rate, the pressure gradually increases during phase 3. The average pressure value is employed as the representative for enthalpy estimation purposes.
this point, Well_3 (producer) was an injection well for over a decade.Subsequently, the bottomhole temperature is initially only 160 °F and rapidly converges the desired range due to the rock heating it and hot fluids arriving at the bottomhole.The temperature observed during the later stages is considered representative of enthalpy calculations.Additionally, due to the ongoing injection of CO2 at a slightly higher rate than the production rate, the pressure gradually increases during phase 3. The average pressure value is employed as the representative for enthalpy estimation purposes.High enthalpy wells can be considered almost isenthalpic, implying that fluid enthalpy remains constant along the well as long as CO2 remains in a supercritical state.The power output can be calculated straightforwardly and is determined to be E = 17.7 MW.To further increase the power output, the injection and production rates can be increased as needed.There is no need for global concern regarding the average pressure increase, as the fluid is injected and produced simultaneously at similar rates.Additionally, there is no local concern regarding pressure buildup due to the high and isotropic permeability of the reservoir.When considering the increasing recycling rates, the main consideration is the time it takes for the returned fluid to reach the reservoir's temperature.Ideally, it is preferable for the fluid to reach the aquifer's temperature before being produced, as the density of supercritical CO2 decreases when an increase in temperature and geothermal energy retrieval is maximized.However, even with the moderate flow rates simulated in this study, the temperature of the produced fluid ultimately converges to a value significantly lower than the temperature of the aquifer.Nevertheless, there are many control options to be exploited in order to optimize the production temperature, such as moving the wells further apart, utilizing deep saline aquifers with higher magmatic activity, or even halting the CPG for a few months in order for the fluid to reach the aquifer's temperature, which actually cannot be considered a practical option for an industrial application.

Conclusions
In conclusion, the integration of the CPG-CCS joint system presents a groundbreaking approach that synergizes energy production with carbon-negative emissions.This study uniquely focused on a deep saline aquifer situated in a region characterized by pronounced subsurface magmatic activity.The outcomes unequivocally demonstrate the remarkable success in sequestering over 36 million tonnes of CO2 over a 25-year period.Additionally, the establishment of a geothermal system effectively harnesses the sequestered CO2 for energy generation.This innovative system offers multifaceted advantages.Firstly, it significantly contributes to carbon negativity by securely storing a substantial volume of CO2.Moreover, it facilitates energy extraction through geothermal means, augmenting the overall energy output.This unique amalgamation of carbon sequestration High enthalpy wells can be considered almost isenthalpic, implying that fluid enthalpy remains constant along the well as long as CO 2 remains in a supercritical state.The power output can be calculated straightforwardly and is determined to be E = 17.7 MW.To further increase the power output, the injection and production rates can be increased as needed.There is no need for global concern regarding the average pressure increase, as the fluid is injected and produced simultaneously at similar rates.Additionally, there is no local concern regarding pressure buildup due to the high and isotropic permeability of the reservoir.When considering the increasing recycling rates, the main consideration is the time it takes for the returned fluid to reach the reservoir's temperature.Ideally, it is preferable for the fluid to reach the aquifer's temperature before being produced, as the density of supercritical CO 2 decreases when an increase in temperature and geothermal energy retrieval is maximized.However, even with the moderate flow rates simulated in this study, the temperature of the produced fluid ultimately converges to a value significantly lower than the temperature of the aquifer.Nevertheless, there are many control options to be exploited in order to optimize the production temperature, such as moving the wells further apart, utilizing deep saline aquifers with higher magmatic activity, or even halting the CPG for a few months in order for the fluid to reach the aquifer's temperature, which actually cannot be considered a practical option for an industrial application.

Conclusions
In conclusion, the integration of the CPG-CCS joint system presents a groundbreaking approach that synergizes energy production with carbon-negative emissions.This study uniquely focused on a deep saline aquifer situated in a region characterized by pronounced subsurface magmatic activity.The outcomes unequivocally demonstrate the remarkable success in sequestering over 36 million tonnes of CO 2 over a 25-year period.Additionally, the establishment of a geothermal system effectively harnesses the sequestered CO 2 for energy generation.This innovative system offers multifaceted advantages.Firstly, it significantly contributes to carbon negativity by securely storing a substantial volume of CO 2 .Moreover, it facilitates energy extraction through geothermal means, augmenting the overall energy output.This unique amalgamation of carbon sequestration and geothermal energy underscores the potential for a system that not only yields surplus energy but also mitigates carbon emissions.These compelling findings underscore the feasibility, importance, and novelty of the CPG-CCS joint system, positioning it as a sustainable and impactful solution to address carbon emissions and fulfill our growing energy needs.Continued research and development in this realm hold immense promise for driving meaningful progress in the fight against climate change.

Figure 2 .
Figure 2. Cumulative mass of CO2 sequestered in the aquifer (10 10 lbm) vs. time (years).At phase A, 308 Mtonnes of CO2 are injected each year, and a percentage of them is produced after the breakthrough.At phase B, 2600 Mtonnes are injected and sequestered each year; therefore, the cumulative mass stored is much higher, represented by a considerable slope increase.Finally, in phase C, CPG is initiated where 2500 Mtonnes are injected and simultaneously 2400 Mtonnes are produced, resulting in a net sequestration rate of 100 Mtonnes per year.

Figure 2 .
Figure 2. Cumulative mass of CO 2 sequestered in the aquifer (10 10 lbm) vs. time (years).At phase A, 308 Mtonnes of CO 2 are injected each year, and a percentage of them is produced after the breakthrough.At phase B, 2600 Mtonnes are injected and sequestered each year; therefore, the cumulative mass stored is much higher, represented by a considerable slope increase.Finally, in phase C, CPG is initiated where 2500 Mtonnes are injected and simultaneously 2400 Mtonnes are produced, resulting in a net sequestration rate of 100 Mtonnes per year.

Figure 3 .
Figure 3. Bottomhole pressures and temperatures for Well_3 and Well_4 during phase C.

Figure 3 .
Figure 3. Bottomhole pressures and temperatures for Well_3 and Well_4 during phase C.

Table 3 .
Well schedule production (P) and injection (I) rate.

Table 3 .
Well schedule production (P) and injection (I) rate.

Table 4 .
Well schedule production (P) and injection (I) mass rate.
Notes: (P) denotes the production rate, (I) denotes the injection rate.Well_1 and Well_2 are dedicated to production.Well_3 is dedicated to injection in phases A and B and switches to production in phase C. Well_4 is dedicated to injection in all phases.
Notes: (P) denotes the production rate, (I) denotes the injection rate.Well_1 and Well_2 are dedicated to production.Well_3 is dedicated to injection in phases A and B and switches to production in phase C. Well_4 is dedicated to injection in all phases.