Sensitivity of Reservoir and Operational Parameters on the Energy Extraction Performance of Combined CO 2 -EGR – CPG Systems

: There is a potential for synergy effects in utilizing CO 2 for both enhanced gas recovery (EGR) and geothermal energy extraction (CO 2 -plume geothermal, CPG) from natural gas reservoirs. In this study, we carried out reservoir simulations using TOUGH2 to evaluate the sensitivity of natural gas recovery, pressure buildup, and geothermal power generation performance of the combined CO 2 -EGR – CPG system to key reservoir and operational parameters. The reservoir parameters included horizontal permeability, permeability anisotropy, reservoir temperature, and pore-size-distribution index; while the operational parameters included wellbore diameter and ambient surface temperature. Using an example of a natural gas reservoir model, we also investigated the effects of different strategies of transitioning from the CO 2 -EGR stage to the CPG stage on the energy-recovery performance metrics and on the two-phase fluid-flow regime in the production well. The simulation results showed that overlapping the CO 2 -EGR and CPG stages, and having a relatively brief period of CO 2 injection, but no production (which we called the CO 2 -plume establishment stage) achieved the best overall energy (natural gas and geothermal) recovery performance. Permeability anisotropy and reservoir temperature were the parameters that the natural gas recovery performance of the combined system was most sensitive to. The geothermal power generation performance was most sensitive to the reservoir temperature and the production wellbore diameter. The results of this study pave the way for future CPG-based geothermal power-generation optimization studies. For a CO 2 -EGR – CPG project, the results can be a guide in terms of the required accuracy of the reservoir parameters during exploration and data acquisition. of establishment stage, the minimum gas (i.e., mainly CO saturation in the reservoir pore space surrounding the production well inside production well These results imply that the amount of water in the bottom-hole of the production well is 2.86% vol. (case 1-A), 5.11% vol. (case 1-B), 9.88% vol. (case 2-A), and 20.84% vol. (case 2-B).


Introduction
The effects of global warming on the environment (including sea-level rise, extreme weather conditions etc.) are on the rise due to anthropogenic emissions of greenhouse gases, namely carbon dioxide (CO2). The development and use of clean, low-carbon, energy-efficient technologies and renewable energy sources are ways of reducing such emissions. To this effect, the technology known as carbon capture and storage (CCS), which features the capture of CO2 from flue gases of power plants (and other large CO2-emitting industrial plants) that are then stored in suitable, carefully selected geological formations, is one of the technologies that can contribute to achieving net-zero CO2 emissions [1][2][3].
(i) The combined system exhibits less pore-water influence (due to the presence of residual natural gas), potentially reducing the resistance to flow of the CO2 in the reservoir. It may also reduce the extent of unwanted geochemical CO2-rock-brine interaction, such as salt precipitation; (ii) Additional natural gas and geothermal energy are extracted for power generation, which leads to an increase in the gas field's total amount of producible energy; (iii) The natural gas-based power generation would likely be operated with CCS, providing the CO2 for the EGR and CPG operations; (iv) Economic (cost-saving) benefits are achieved by using/sharing already-existing multidisciplinary datasets (on reservoir parameters) and infrastructure (surface facilities, wells etc.). Hence, investment costs are significantly reduced; (v) The combined system extends the useful lifetime of the gas reservoir, recovering otherwise stranded assets, such as wells, offshore platforms etc., thereby postponing the expensive decommissioning phase of the wells and abandonment stages of the gas field.
Different stages were identified for energy (and power) generation from the combined system (see also ref. [45]). The first stage is the conventional natural gas recovery (CNGR) stage. Here, natural gas is produced by the primary-recovery drive. Heat is extracted from the natural gas at the surface and converted to power using an organic Rankine cycle (ORC) or a CO2-based Rankine cycle (CRC). The CNGR stage is followed by the CO2-EGR stage when the fluid pressure or natural gas is depleted and the remaining natural gas cannot be economically recovered by just the natural primary drive. CO2 is injected into the reservoir to recover the remaining natural gas and reduce the residual methane content in the reservoir. The associated heat in the produced and remaining natural gas is extracted and converted to power as described above. At some point during the CO2-EGR stage, CO2 breakthrough occurs in the production well. Hence, the installation of a methane-CO2 separator is required at the surface.
A transition period (TP) exists at the end of the CO2-EGR stage. It is the period when the mass fraction of CO2 in the produced fluid is greater than 10% and less than 90%. We assume that the separation of the produced mixed fluid is necessary. The separation process of CO2 and methane (CH4) is cost-and energy-intensive, requiring higher energy input than the geothermal power generated during the CO2-EGR stage. The longer this transition period lasts, the higher the parasitic power required to separate the mixed fluid at the surface.
Ref. [45] described a CO2-plume establishment (PE) stage, where the production well is shut-in and CO2 injection is continued to charge the reservoir with CO2 and establish a CO2 plume between the injection and the production wells. The PE stage can come before or after the CO2-EGR stage (the latter is the case in the simulations shown later), depending on the operational strategy.
After these stages, the CPG stage commences and all produced fluid (mostly CO2) is sent to the CO2 turbine for direct-CPG [41,43,44] power generation. The CO2 leaving the turbine is further cooled and reinjected into the reservoir. When the reservoir heat is depleted, after decades as shown in refs. [41,44], the CPG-based injection and production wells are shut down and the injected CO2 is permanently stored in the natural gas reservoir.
In this study, we used a similar anticlinal natural gas reservoir model, as described in ref. [45], to carry out reservoir simulations in TOUGH2 [52,53], aimed at expanding the previous study on the combined CO2-EGR-CPG system to: (i) Accommodate lessons learned from ref. [54], including that the bottom-hole production flowrate directly influences how much water enters the production well and that the two-phase (water/CO2) flow regime in the production well is an important design parameter; (ii) Determine the effect of the residual CH4 content of the reservoir and the effect of the CO2-plume-establishment stage on the performance of the combined system, as quantified by the amount of natural gas recovery, fluid-pressure buildup, and electric power generation; (iii) Determine the sensitivity of the performance metrics mentioned in (i) and (ii) to various important reservoir and operational (non-reservoir) parameters.
Hence, this work provides a preliminary guide for an effective implementation strategy of the combined CO2-EGR-CPG system. The sensitivity study provides information necessary for the optimization of such a system.

Methodology
In this section, we describe the natural gas model, numerical simulation, and performance metrics used for this study. We report the changes we made to the original natural gas reservoir model (presented in ref. [45]) and the implementation strategy/concept of the combined system for effective natural gas and geothermal energy extraction from the natural gas reservoir. We describe the numerical model we used to simulate the three main stages (CNGR, CO2-EGR+TP, and CPG stages) associated with the combined system. Finally, we explain how the sensitivity study was conducted.

Reservoir Model
Using the reservoir properties from some examples of hot/deep natural gas fields around the world, we set up a similar natural gas reservoir to that of our previous study (see ref. [45]), in terms of model geometry, well configuration, boundary conditions, and relative permeability functions. However, the fluid properties of the current model were updated, such that the salt concentration changed from zero in the previous model to 150,000 ppm in the current model. Furthermore, the y-axis dimension of the full model was extended from the previous 3 km to 4.5 km, so that the updated model had equal dimensions on all sides. We also introduced a horizontal-to-vertical permeability anisotropy of kh/kv = 2 in the base-case model. We conducted our simulations using the same reservoir simulator, TOUGH2 with the EOS7C module, developed for simulating gas and brine flow as well as heat transport in (natural gas) reservoirs [52,53]. The rock and fluid properties of the new model, as well as the initial conditions, are summarized in Table 1. Note that the model comprised of two fluid phases: brine and gas. The gas phase may contain CO2 and/or CH4.

. Reservoir Simulation Schemes
Based on the model of the combined CO2-EGR-CPG described above, the numerical simulations were carried out in three main stages: a. The CNGR stage with a base-case duration of 25 years at a production flowrate of 4 kg/s/well; b. The CO2-EGR stage with all simulated cases lasting for 1 year with a high injectionproduction flowrate ratio [45,51] (i.e., injection flowrate of 30 kg/s/well and production flowrate of 4 kg/s/well). The high injection-production flowrate ratio is beneficial for achieving a good CO2-EGR performance as well as a short duration for establishing an adequate CO2-plume reservoir [51]. After the 1 year period of high injectionproduction flow-rate ratio, the production rate was increased to 30 kg/s/well (equal to the injection flowrate). The CO2-EGR stage ended when the mass fraction of CO2 (XCO2) at the production-well region reached 90%. In this study, we considered some cases that were associated with 1.5 years of the CO2-plume establishment (PE) stage after the 1 year CO2-EGR stage. After the PE stage, the transition period, which involves CO2-CH4 separation, continued till 90% CO2 mass fraction (in the gas phase) was reached at the production well region; c. Finally, the CPG stage, where the base-case CPG-stage duration was 30 years at a circulation flowrate of 30 kg/s/well.

Performance Metrics
In this section, we present how we calculated the metrics used to measure the energyextraction performance for the different cases considered in this study. These metrics included: (a) The CO2 saturation (in the reservoir and in the production well), and the corresponding flow regime established (at the bottom-hole section of the production well) at the time of the highest water saturation around the production-well inlet region of the reservoir (ref. [45]). This performance metric was only applicable for the reservoir parameters, and it can be used to determine the importance of the PE stage for the combined system to achieve an annular flow regime (dominant CO2 flow) in the production well. The method to calculate this metric, using the gas saturation in the well and flowrate, can be found in ref. [54]; (b) The natural gas recovery performance (NGRP), which incorporates the amount of natural gas recovered (in terms of natural gas recovery factor) and the duration of the CNGR and CO2-EGR stages (including the transition period, TP). This implies that the duration of natural gas recovery also includes the time of CH4-CO2 separation at the land surface. The natural gas recovery factor during the CNGR and CO2-EGR stages can be measured as a percentage of the original gas in place (OGIP). These respective factors are calculated as: where −CNGR is the volume of gas produced during the CNGR stage, and OGIP is the volume of gas initially in place, where −EGR is the volume of gas produced during the CO2-EGR stage, Ultimate recovery factor, UR (%) = CNGR + EGR The volumes of the produced gas during the CNGR and CO2-EGR stages can be obtained from the TOUGH2 simulation output files. The natural gas recovery performance was measured using a natural gas recovery index (RI), which is calculated as the ultimate recovery factor [%] divided by the sum of the durations of the CNGR stage, CNGR [year] and CO2-EGR+TP stage, EGR [year] (Equation (4)).
The recovery index provides a way to select the best reservoir parameters and strategies that favor natural gas recovery and a shorter duration of the CO2-CH4 separation. The shorter this transition period is, the higher the energy efficiency of the combined system. Hence, a high value of RI is favorable; (c) The pressure buildup at the injection wells was one of the performance metrics considered for the combined system. It can compromise the integrity of the caprock overlying the natural gas reservoir. The non-dimensional pressure buildup (PBU) metric, used in this study, compares the maximum pressure in the injection well, wb−max , to the initial reservoir pressure, 0 , as shown in Equation (5).
where is the injectivity index, calculated as − is the average of the maximum pressures in the injection grid cells of the two injection wells. The perforation layer thickness, , is 20 m, the well radius, , is 0.07 m, and the effective gridblock radius, = √ ⁄ . Here, A is the grid block area, which is 100 m 2 . The is the CO2 mass flowrate [kg/s], is the relative permeability [-], and and are the density [kg/m 3 ] and the dynamic viscosity [Pa‧s] of the gas (i.e., CO2) phase, respectively; (d) The average net geothermal electricity (measured in gigawatt-hours [GWeh]) generated using the produced natural gas (via the organic/CO2-based Rankine cycle) during the CNGR, CNGR , and EGR, EGR , stages and from the produced CO2 (via the direct CO2 turbomachinery) during the direct-CPG stage, CPG . The average net power generated was calculated using the output wellhead temperature and pressure results obtained from the wellbore heat transfer model described in ref. [45]. The power system models applied in this study for the indirect and direct CO2 turbomachinery power systems were extensively described in ref. [43] and ref. [45].

Residual Methane Content and CO2-Plume Establishment Stage
In our simulations, the residual CH4 content was accounted for in the simulations using two main cases: case 1 (base case) and case 2, which considered the CNGR stage for 25 years and 26 years, respectively ( Table 2). This implies that case 2 had a lower residual CH4 content than case 1 (Figure 1). Case 2 was equivalent to the example case presented in ref. [45]. Case 2 is an example of a depleted reservoir, while case 1 is an example of a partially depleted natural gas reservoir because it had a higher proportion of the residual CH4 in the reservoir after the CNGR stage. These two cases were further classified into four sub-cases to represent cases with and without the 1.5 years of the CO2-plume establishment (PE) stage. Cases 1-A and 2-A considered the 1.5 years of the PE stage, whereas cases 1-B and 2-B did not consider the PE stage (Table 2). This enabled us to investigate the effects of the PE stage on the production performance of the combined system in terms of ensuring annular flow in the production well, and in reducing the duration of the "expensive" transition period. The corresponding effects of these four sub-cases on the minimum CO2 saturation in the well, natural gas recovery (during the CNGR and CO2-EGR stages-including the transition period), the geothermal energy generation (during the three stages), and performance metrics were also evaluated and discussed in this study (in Section 3.1).
(a) (b) Figure 1. Quarter model showing the gas (i.e., CH4) saturation in the reservoir pore space after the CNGR stage for the partially depleted (a) and depleted (b) natural gas reservoir model. Note that the part of the model with gas (i.e., CH4) saturation less than 15% (i.e., high brine saturation) has been blanked.

. Reservoir and Operational Parameters
We varied some key reservoir and operational (non-reservoir) parameters in a fixed range of ±5%, as compared to the base-case conditions (presented in Table 1), to investigate the sensitivity of the performance metrics of the combined system described in Section 2.2 to these parameters. The reservoir parameters (and their respective base-case values) that were studied include permeability (100 mD), anisotropy kh/kv (2.0 [-]), relative permeability (van Genuchten parameter of 0.77 [-]), and initial temperature (120 °C). The operational parameters (and their respective base-case values) that were studied include well diameter (0.14 m) and average ambient surface temperature (15 °C). These base-case values corresponded to the values chosen for the reservoir example presented in ref. [45], except the reservoir initial temperature (150 °C was chosen in our previous study).
The simulations for each of the reservoir and operational parameter spaces considered in this study were run for a duration of 25 years for the CNGR stage and without considering the PE stage. The simulation results (high and low bounds) of the perfor-mance metrics for each of the respective parameter spaces were compared to the simulation results of the base-case example presented as case 1-B (i.e., partially depleted natural gas reservoir without the PE stage).
The comparison indicated the sensitivity on the performance metrics (m) for a 5% change (both high and low) in the value of the parameter (X) from the base-case value considered in this study. Mathematically, we represent this sensitivity value, σ, as where 0 and 0 are the base-case performance metric result and the value of the parameter for the base case, respectively. The subscripts "high" and "low" denote the respective upper and lower bounds of the 5% deviation from the base-case value. We chose ( high − low )/ 0 = 0.1 (i.e., 5% both ways from the base case). This implies that σ is a function of the performance metric results, , and Equation (8) becomes: It should be noted that the non-reservoir (operational) parameters considered in this study do not affect the natural gas recovery performance and the pressure-buildup metrics because these performance metrics are only dependent on the reservoir conditions. In addition, the sensitivity values obtained for the minimum CO2 saturation in the well did not show any significant changes with the 5% deviation from the base case.

Effects of Residual CH4 Content and CO2-Plume Establishment Stage on the Performance Metrics
When the water saturation around the production well inlet was the highest during the CO2-EGR+TP stage, we considered whether the CO2-plume establishment (PE) stage was required to ensure the desired [49] annular flow regime in the bottom-hole region of the production well. Recall that "A" denotes that the PE stage exists, whereas "B" denotes that the PE stage does not exist. Figure 2 shows that, at the time of highest water saturation at the bottom-hole, case 1-A had the highest CO2 saturation in the reservoir and in the production well. This implies that the PE stage helped reduce the amount of water entering the production well. The results in Figure 3 show that the lower the residual CH4 content after the CNGR stage, the higher the probability that the production well may experience slug/churn flow. However, for relatively low-diameter production wells (14 cm and 21 cm), the desired [49] annular flow was sustained even when the residual CH4 content was low (cases 2-A and 2-B). For a production-well diameter of 33 cm, case 2-B (depleted reservoir with no CO2 plume establishment option) exhibited slug/churn flow at the bottom of the production well. Hence, if the diameter of the production well is relatively large, it may be necessary to include the CO2 plume establishment (PE) stage. As discussed in ref. [54], the optimal diameter of the production well needs to be determined not only to achieve an annular flow regime, but also to maintain minimal pressure and heat losses in the production well.
We also investigated the effect of residual CH4 content (in terms of partially depleted or depleted natural gas) on the natural gas recovery performance. The natural gas recovered during CNGR was greater for the depleted reservoir scenario (see Table 3 and Figure  A1a in Appendix A). The natural gas recovery during the CO2-EGR stage was higher for the cases without CO2-plume establishment (case 1-B and case 2-B). This is because of the relatively high CO2-CH4 mixing rate that is associated with the CO2-plume establishment stage, which may slightly reduce natural gas recovery during the CO2-EGR+TP stage. Table 3. Simulation results of the percentages of the original gas in place (OGIP), recovered during the CNGR and CO2-EGR+TP stages, and associated natural gas recovery index (RI). From Table 3, we observed a decrease in the number of years for the expensive separation period in the CO2-EGR+TP stage for the two cases that considered the PE stage, which led to a higher recovery index (RI). For the partially depleted natural gas reservoir example cases, the expensive separation duration reduced when the PE stage was considered (i.e., case 1-A), and this led to better performance (higher RI) for case 1-A than for case 1-B. Table 3 also shows that the calculated RI was lower for the depleted natural gas reservoir example cases, and that there was only a slight increase in RI between cases 2-A and 2-B. Hence, the simulation results showed that the effect of the PE stage on the natural gas recovery index was more pronounced when the reservoir was partially depleted of natural gas. However, if the natural gas reservoir was depleted, the PE stage was important to establish a CO2 connection between the injectors and producers. During the CNGR stage, more geothermal energy was generated from the natural gas for the depleted reservoir cases (case 2) than the partially depleted reservoir cases (case 1), shown in Table 4 and Figure A2, as the amount of natural gas recovered from the depleted natural gas reservoir was greater than that recovered from the partially depleted natural gas reservoir. The opposite was the case for the CO2-EGR+TP stage, during which more natural gas, and hence more associated geothermal energy, was produced from the partially depleted reservoir (see Table 3). During the CPG stage, case 1-A had the highest value of geothermal electricity generated (Table 4). This was because it had the shortest duration of the CO2-EGR+TP stage (1.50 years), which ensured an earlier start, and longer duration, of the CPG stage.

Cases
Of all four cases considered, the combined CO2-EGR-CPG system in a partially depleted natural gas reservoir, which includes a PE stage (case 1-A), yielded the best results in terms of: (i) the desired [49] annular flow regime in the production well, (ii) the natural gas recovery performance (including having the shortest duration of the CO2-EGR+TP stage), and (iii) the average total net electricity generation performance. Hence, it appears advantageous to plan the combined CO2-EGR-CPG system (including the PE stage) to commence before the natural gas reservoir is completely depleted. , where G−well is the gas-phase saturation in the well, L−well is the water-phase saturation in the well, and denotes density. The plotted symbols show the lowest gas-phase mass fraction points of the four investigated cases ( Table 2) for three production well diameters, D. The lines indicate the minimum flowrate required to achieve the annular flow regime in the production well. Hence, conditions above each line indicate the annular flow regime, while conditions below each line indicate the undesired slug/churn flow regime in the production well.

Natural Gas Recovery Performance and Maximum Fluid-Pressure Buildup
Recall that the natural gas recovery performance (NGRP) metric is only dependent on the changes in the reservoir parameters and not on non-reservoir parameters. As described in Section 2.2, the NGRP was determined in this study by the value of the natural gas recovery index (RI), which can be calculated using Equation (4). Figure A1b-e show the volume of natural gas recovered (in %OGIP) during the CNGR and CO2-EGR stages for the different reservoir parameters considered in this study. Figure 4 shows the sensitivity for the reservoir parameters considered in this study on the natural gas recovery performance for the CNGR and the CO2-EGR+TP stages (when natural gas was being produced and methane was separated from CO2). The parameters plotted on the upper part of the horizontal logarithmic axis (Figure 4) indicated those with a positive sensitivity value (σ), and those parameters that were plotted on the lower part of the axis have a negative σ. A positive (+) value of σ shows that there was an increase in the NGRP when the value of the parameter increased, whereas the opposite was the case for a negative (-) value of σ.
The sensitivity results, plotted in Figure 4, show that permeability anisotropy and reservoir temperature were the most sensitive parameters (both ~0.2) to the NGRP, with each having a different sign of the sensitivity value. We observed that, as the permeability anisotropy increases, the NGRP significantly decreased (i.e., negative σ value). This was the case as the production well was perforated only over the topmost layer interval (20 m out of the entire 100 m reservoir thickness), and a decrease in the vertical permeability (i.e., increase in permeability anisotropy) led to less upward flow of the natural gas towards the production well perforation. Hence, the decrease in NGRP. However, an increase in reservoir temperature would lead to a significant increase in the NGRP, as the kinematic viscosity of the gas is reduced at higher temperatures, increasing the gas mobility. The horizontal permeability parameter was the least sensitive parameter for the NGRP, with a sensitivity value of +0.0045, which was two orders of magnitude less than those of the most sensitive parameters. . Sensitivity values (plotted on a logarithmic axis) for the natural gas recovery performance (NGRP) metric for the reservoir parameters considered in this study for the CNGR and the CO2-EGR+TP stages. The positive (upper) part signifies that an increase in the parameter values will improve the NGRP. The negative (lower) part identifies parameters that, when being increased, will result in a decrease in the NGRP.
The sensitivity results for the maximum fluid-pressure buildup, plotted in Figure 5, showed that the van Genuchten relative permeability parameter was the most sensitive parameter, followed by horizontal permeability and reservoir temperature. The sensitivity to the van Genuchten parameter was because the maximum fluid-pressure buildup happened during the beginning of the CO2-EGR+TP stage when two-phase flow processes near the injection wells are important. An increase in horizontal permeability obviously leads to a decrease in pressure-buildup intensity, hence a high (negative) sensitivity value was observed for this parameter. The sensitivity to reservoir temperature was a result of the constant CO2 injection mass flowrate used. If a constant volume flowrate was applied instead, the fluid-pressure buildup would be significantly less sensitive to reservoir temperature. Figure 5 also shows that the fluid-pressure-buildup metric was not sensitive to permeability anisotropy changes. Figure 5. Sensitivity values (plotted on a logarithmic axis) for the fluid-pressure buildup metric for the reservoir parameters considered in this study. The positive (upper) part signifies that an increase in the parameter values will increase the fluid-pressure buildup. The negative (lower) part identifies parameters that, when being increased, will result in a decrease in the fluid-pressure buildup.

Geothermal Energy (Electricity) Generation Performance during the CNGR, CO2-EGR+TP, and CPG Stages
During the CNGR stage, the average net power and electricity generated by the basecase example (i.e., case 1-B, described in Section 3.1) were 0.042 MWe and 9.16 GWeh over 25 years. The time series plots of the simulation results, showing the net power generated for the different parameters considered in this study, can be found in the Appendix A ( Figures A3-A8). Figure 6 shows that the most sensitive parameter that influenced the net electrical geothermal energy generated during the CNGR stage was the reservoir temperature (σ = +5.0), followed by the mean surface ambient temperature (σ = -0.5). It is no surprise that the net energy generated increased with an increase in reservoir temperature, but decreased with an increase in the mean surface ambient temperature. The permeability anisotropy parameter was the least sensitive parameter for this performance metric during the CNGR stage. The sensitivity value (-0.0015) for the permeability anisotropy differed by 3 orders of magnitude from that of the reservoir temperature parameter. Figure 6 also shows that, during the CNGR stage, changes in the other parameters considered in this study had little effect on the geothermal energy generation performance, with their sensitivity values of σ < ±0.1. Figure 6. Sensitivity values (plotted on a logarithmic axis) for the geothermal electricity generation performance metric for the reservoir and operational parameters considered in this study. This is calculated for the three main stages of the combined system. The positive (upper) part signifies that an increase in the parameter values will increase the geothermal power output. The negative (lower) part identifies parameters that, when being increased, will result in a decrease in the geothermal power output.
During the CO2-EGR+TP stage, the average net power and the electricity generated for the base-case example were 0.053 MWe and 1.78 GWeh (over 3.82 years). From Figure  6, we observed that the most sensitive parameter for the CO2-EGR+TP stage was the production wellbore diameter, with σ = +1.50. This is because during this stage, the CO2 mass fraction increases, leading to a higher density of the fluid flowing in the well and causing higher pressure losses due to friction as the well diameter decreases. This reduces the amount of energy extractable from the produced/separated methane. The least sensitive parameter was the mean ambient surface temperature, with σ = +0.028.
During the CPG stage (when the CO2 mass fraction was >90% and only the direct CO2-thermosiphon power system model was used), the average net power and the geothermal energy generated for the base case example were 0.739 MWe and 173.44 GWeh (over 27.18 years). This was about 17.6 and 13.9 times higher than the average net power generated during the CNGR and CO2-EGR+TP stages, respectively. This is due to the use of the direct CPG power system [40] and the increase in production mass flowrate from 4 kg/s to 30 kg/s applied during this stage.
As expected, the most sensitive parameter for this CPG stage was reservoir temperature, with a sensitivity value of σ > +10.0 ( Figure 6). The wellbore diameter and the mean surface ambient temperature also showed high sensitivity values because, as the wellbore diameter increases, the fluid pressure and the heat losses decrease, while lower mean surface ambient temperatures favor higher net energy extraction and power generation rates at the land surface. The least sensitive parameter during this CPG stage was the relative permeability with a σ that was more than 2 orders of magnitude lower than that for the reservoir temperature. The relatively sensitivity of horizontal permeability was due to the fact that a constant CO2 flowrate was used in this study. If, instead, the fluid-pressure difference between the injection and the production wells were kept constant, permeability would be a key parameter.
In this study, we observed that reservoir temperature was the most important parameter for the combined CO2-EGR-CPG system, especially during the CNGR and CPG stages. The second most important parameter was the production-wellbore diameter, which can be adjusted, i.e., engineered. Interestingly, the performance metrics, except for pressure buildup, studied here were not particularly sensitive to the van Genuchten relative-permeability parameter for pore-size distribution.

Conclusions
We present a reservoir simulation study of a combined CO2 enhanced gas recovery (EGR)-CO2 plume geothermal (CPG) system. We investigated how the residual CH4 content and the CO2-plume establishment (PE) stage affected the performance of the CO2-EGR-CPG system. This performance was quantified by assessing the amount of natural gas recovered, the amount of geothermal energy generated, and whether an annular flow regime could be achieved in the production well when a specified mass flowrate is used. In addition, the sensitivity of the system performance to key reservoir and operational parameters was investigated. Based on the simulation results, our key results and findings were: 1. The PE stage was important to establish the desired annular flow regime near the bottom of the production well, especially when large-diameter wells were used. However, it was possible to achieve an annular flow regime by using wells with smaller diameters; 2. Commencing with the injection of CO2 into the natural gas reservoir before it was completely depleted and including the PE stage increased the natural gas recovery performance of the combined CO2-EGR-CPG system; 3. Of the four reservoir parameters (permeability anisotropy, horizontal permeability, relative permeability, and reservoir temperature) considered in this study, permeability anisotropy and reservoir temperature were the parameters that most strongly affected the natural gas recovery performance. The fluid pressure buildup at the injection wells was most sensitive to changes in the van Genuchten relative permeability and horizontal-permeability parameters; 4. The sensitivity results revealed that, for a given CO2 flowrate, the reservoir temperature was the parameter the geothermal power generation performance was most sensitive to. The production wellbore diameter was the second most sensitive parameter.
Changes in the van Genuchten relative permeability parameter for pore-size distribution did not have a significant influence on the geothermal power generation performance of the combined CO2-EGR-CPG system.