Salaheddine Chabab
email: salaheddine.chabab@univ-pau.fr

José Lara Cruz

Marie Poulain
email: marie.poulain@univ-pau.fr

Marion Ducousso
email: marion.ducousso@univ-pau.fr

François Contamine
email: francois.contamine@univ-pau.fr

Jean Paul

Pierre Cézac
email: pierre.cezac@univ-pau.fr

Thermodynamic Modeling of Mutual Solubilities in Gas-Laden Brines Systems Containing CO 2 , CH 4 , N 2 , O 2 , H 2 , H 2 O, NaCl, CaCl 2 , and KCl: Application to Degassing in Geothermal Processes

Keywords: CO 2 emissions, geothermal, degassing, gas solubility, bubble-point pressure

With the growing interest in geothermal energy as a renewable and sustainable energy source, nowadays engineers and researchers are facing technological and environmental challenges during geothermal wells' operation or energy recovery improvement by optimizing surface installations. One of the major problems encountered is the degassing of geothermal brines which are often loaded with dissolved gases, resulting in technical problems (scale formation, corrosion, reduced process efficiency, etc.) and environmental problems through the possible emission of greenhouse gases (CO 2 , CH 4 and water vapor) into the atmosphere. In this work, a method to predict, from readily available information such as temperature and GLR, the bubble point pressure of geothermal fluids as well as the GHG emission rate depending on the surface conditions is presented. This method is based on an extended version of the Soreide and Whitson model with new parameters optimized on the solubility data of several gases (CO 2 , CH 4 , N 2 , O 2 and H 2 ) in brine (NaCl + CaCl 2 + KCl). The developed approach has been successfully used for the prediction of water content of different gases and their solubilities in different types of brines over a wide temperature and pressure range, and has been applied for the prediction of bubble point pressure and GHG emissions by comparing the results with available industrial data of geothermal power plants including the Upper Rhine Graben sites.

Introduction

The development of renewable energy is one of the key strategies for contributing to the decrease of greenhouse gas (GHG) emissions [START_REF] Poponi | Energy Technology Perspectives 2016: Towards Sustainable Urban Energy Systems[END_REF]. Geothermal energy, used for heat and electric power generation [START_REF] Dickson | What Is Geothermal Energy; Istituto di Geoscienze E Georisorse[END_REF], appears to be a sustainable option to facilitate this energy transition. Depending on the depth of the reservoir containing the geothermal fluid, two applications are distinguished: shallow geothermal and deep geothermal [START_REF] Pan | Establishment of enhanced geothermal energy utilization plans: Barriers and strategies[END_REF][START_REF] Haehnlein | International legal status of the use of shallow geothermal energy[END_REF]. In some deep geothermal reservoirs, ground permeability is too low to properly extract the fluid from its well. Thus, the geothermal resource heat is retrieved with an Enhanced Geothermal System (EGS) [START_REF] Pan | Establishment of enhanced geothermal energy utilization plans: Barriers and strategies[END_REF], which boosts rock porosity to extract the proper amount of fluid for heat recovery [START_REF] Lu | A global review of enhanced geothermal system (EGS)[END_REF].

Currently, EGS projects are ongoing in Australia, France, Germany, the United Kingdom and the United States. However, their development stages are different from one country to another. While some of them are at plant construction and well drilling phases, such as those in Australia and the United Kingdom, others are commercial-scale plants, such as those in French territory [START_REF] Lu | A global review of enhanced geothermal system (EGS)[END_REF][START_REF] Pratiwi | Life-cycle climate-change impact assessment of enhanced geothermal system plants in the Upper Rhine Valley[END_REF][START_REF] Curtis | Geothermal energy use, country update for United Kingdom[END_REF].

France is a leading country in EGS development, notably thanks to the Soultz-sous-Fôrets project, which started in 1987 [START_REF] Lu | A global review of enhanced geothermal system (EGS)[END_REF]. The Soultz-sous-Fôrets site, at the Upper Rhine Graben, produces electricity (12 GWh/year) from geothermal fluids [START_REF] Boissavy | Geothermal energy use, country update for France[END_REF], and it was the first commercial-scale EGS power plant in the world [START_REF] Lu | A global review of enhanced geothermal system (EGS)[END_REF]. Deep geothermal on the Upper Rhine Graben, with the EGS technology, is also carried out at Rittershoffen, where heat is produced (180 GWh/year) for industrial needs [START_REF] Pratiwi | Life-cycle climate-change impact assessment of enhanced geothermal system plants in the Upper Rhine Valley[END_REF][START_REF] Mouchot | Deep geothermal plants operation in Upper Rhine Graben: Lessons learned[END_REF]. EGS projects for heat and electricity production are currently under development on the region, at Illkirch-Graffenstaden and Eckbolsheim [START_REF] Pratiwi | Life-cycle climate-change impact assessment of enhanced geothermal system plants in the Upper Rhine Valley[END_REF][START_REF] Boissavy | Geothermal energy use, country update for France[END_REF].

Usually, geothermal fluids are pumped from their wells to the surface, where they pass through heat exchangers to recover their energy. Nevertheless, geothermal fluids are often hot brines containing dissolved gases, and when they are sent to the surface, their pressure decreases and the dissolved gases evolve from their liquid phase [START_REF] Dickson | What Is Geothermal Energy; Istituto di Geoscienze E Georisorse[END_REF]. The pressure under which this degassing occurs is the bubble point pressure [START_REF] Yearsley | Common characteristics of pumped wells in geothermal power projects[END_REF][START_REF] Haizlip | The impact of high non condensable gas concentrations on well performance-Kizildere Geothermal Reservoir, Turkey[END_REF].

On the one hand, the presence of a gas phase in the geothermal fluid sent to the surface reduces the process efficiency. Firstly, heat transfer efficiency in the heat exchangers declines if the fluid entering these devices is a gas-liquid mixture rather than a single-phase liquid. Moreover, depending on the geothermal brine composition, degassing leads to changes on its pH. Then, salt scaling (e.g., calcite-CaCO 3 ) could arise in the heat exchangers [START_REF] Mouchot | Deep geothermal plants operation in Upper Rhine Graben: Lessons learned[END_REF] or in the production well [START_REF] Haizlip | The impact of high non condensable gas concentrations on well performance-Kizildere Geothermal Reservoir, Turkey[END_REF][START_REF] Yildirim | High total inorganic carbon concentration dependent carbonate scaling and mitigation system in moderate to high enthalpy geothermal fields in Turkey[END_REF], impairing the process efficiency. On the other hand, carbon dioxide (CO 2 ) and methane (CH 4 ), both GHGs, are often found in geothermal fluids. Thus, their degassing can also cause undesired GHG emissions [START_REF] Fridriksson | Greenhouse gas emissions from geothermal power production[END_REF][START_REF] Júlíusson | Tackling the challenge of H 2 S emissions[END_REF][START_REF] Mclean | Greenhouse gas emissions from New Zealand geothermal power generation in context[END_REF], when the geothermal resource is at the surface.

Hence, it is necessary to determine the bubble point depth to decide where a pump should be placed in the well to increase the fluid pressure [START_REF] Yearsley | Common characteristics of pumped wells in geothermal power projects[END_REF] and minimize the amount of gas sent to the surface. The sites of Rittershoffen and Soultz-sous-Fôrets are pumped-well geothermal plants [START_REF] Pratiwi | Life-cycle climate-change impact assessment of enhanced geothermal system plants in the Upper Rhine Valley[END_REF][START_REF] Mouchot | Deep geothermal plants operation in Upper Rhine Graben: Lessons learned[END_REF]. In Rittershoffen, the production well contains a pump set at a depth of 460 m [START_REF] Pratiwi | Life-cycle climate-change impact assessment of enhanced geothermal system plants in the Upper Rhine Valley[END_REF], whereas the Soultz-sous-Fôrets production well presents a pump immersed at 350 m [START_REF] Mouchot | Deep geothermal plants operation in Upper Rhine Graben: Lessons learned[END_REF]. The bubble point pressure of these sites geothermal fluids is estimated to be about 2.0 MPa [START_REF] Mouchot | Deep geothermal plants operation in Upper Rhine Graben: Lessons learned[END_REF].

In some cases, despite the efforts to avoid its degassing, the geothermal resource arrives flashed at the wellhead of plants such as those in Rittershoffen and Soultz-sous-Fôrets. Thus, a vapor-liquid separator ensures that only the liquid fraction of the geothermal fluid is sent to the heat exchangers, while the gas phase, containing the GHGs, is vented to the atmosphere.

Most of the published data on gas emissions from geothermal sites focus on carbon dioxide emissions, which may be interesting information, as this gas is the most abundant in geothermal fluids [START_REF] Fridriksson | Greenhouse gas emissions from geothermal power production[END_REF]. These data concern only geothermal power plants, whose average carbon dioxide emission rate is 122 gCO 2 •kWh -1 [16]. However, this average emission rate is lower than those typically observed in fossil fuel power plants: 450 gCO 2 •kWh -1 in natural gas plants; 675 gCO 2 •kWh -1 in oil plants; and 900 gCO 2 •kWh -1 in coal plants [START_REF] Fridriksson | Greenhouse gas emissions from geothermal power production[END_REF].

Moreover, the configuration of the plant can influence its GHG emissions. The EGS sites in the Upper Rhine Graben employ process configurations where the geothermal fluid extracted from the well is reinjected in the ground [START_REF] Lu | A global review of enhanced geothermal system (EGS)[END_REF][START_REF] Pratiwi | Life-cycle climate-change impact assessment of enhanced geothermal system plants in the Upper Rhine Valley[END_REF]. Gas emissions would occur only at the plant liquid-vapor separator. These type of sites should emit low GHG quantities [START_REF] Pratiwi | Life-cycle climate-change impact assessment of enhanced geothermal system plants in the Upper Rhine Valley[END_REF][START_REF] Fridriksson | Greenhouse gas emissions from geothermal power production[END_REF][START_REF] Fridleifsson | The possible role and contribution of geothermal energy to the mitigation of climate change[END_REF]. However, in some extreme cases, such as those of geothermal power plants in Turkey, where the geothermal fluid is also reinjected in its reservoir, carbon dioxide emissions (1050 gCO 2 •kWh -1 ) [START_REF] Aksoy | Power generation from geothermal resources in Turkey[END_REF] can exceed those of some coal power plants [START_REF] Layman | Geothermal projects in Turkey: Extreme greenhouse gas emission rates comparable to or exceeding those from coal-fired plants[END_REF]. This stems from the fact that the geothermal reservoirs exploited in this country are rich in carbonate-dominated rocks [START_REF] Haizlip | Origin and impacts of high concentrations of carbon dioxide in geothermal fluids of western Turkey[END_REF].

To the best of our knowledge, there is no available data on GHG emissions arising from the degassing of geothermal fluids at the Upper Rhine Graben sites during the process of energy production. However, as this environmentally harmful phenomenon can actually occur, it would be interesting to assess it at these geothermal sites.

To estimate GHG emissions and optimize the process, knowing the breakout pressure of geothermal fluid is essential. Thermodynamic characterization of the liquid-vapor equilibrium of these fluids, which depends on downhole composition data, is necessary to predict their boiling points. Each formation brine is unique and its composition is intimately linked to the geochemistry of the place where it is recovered on the planet [START_REF] Iberl | Evaluation of corrosion of materials for application in geothermal systems in Central Europe[END_REF]. This huge diversity of fluid composition participates in the complexity of this kind of study. Moreover, it could be difficult to determine the brine's composition in a reservoir due to technical difficulties in setting up downhole sampling, related to temperature and pressure conditions. Therefore, although some articles describe techniques to perform this type of sampling [START_REF] Klyen | Sampling Techniques for Geothermal Fluids[END_REF][START_REF] Brown | Precious metals in high-temperature geothermal systems in New Zealand[END_REF][START_REF] Arnórsson | Wet-steam well discharges. I. Sampling and calculation of total discharge compositions[END_REF][START_REF] Wolff-Boenisch | Review of available fluid sampling tools and sample recovery techniques for groundwater and unconventional geothermal research as well as carbon storage in deep sedimentary aquifers[END_REF], downhole compositions are usually obtained with thermodynamics models from physical and chemical data related to the fluid recovered at wellhead [START_REF] Yildirim | High total inorganic carbon concentration dependent carbonate scaling and mitigation system in moderate to high enthalpy geothermal fields in Turkey[END_REF][START_REF] Pátzay | Modeling of scale formation and corrosion from geothermal water[END_REF][START_REF] Pátzay | Preliminary investigations of scaling and corrosion in high enthalpy geothermal wells in Hungary[END_REF][START_REF] Akin | Modeling of calcite scaling and estimation of gas breakout depth in a geothermal well by using PHREEQC[END_REF][START_REF] Wamalwa | Evaluation of factors controlling the concentration of non-condensible gases and their possible impact on the performance of wells in Olkaria, Kenya[END_REF]. Among these data, one of the main parameters used to deduce downhole composition is the Gas-Liquid Ratio (GLR), which can be obtained by flashing a sample of the geothermal fluid at atmospheric pressure and wellhead temperature. Then, the gas and liquid volumes issued from this flash operation are measured, and the GLR can be expressed in gas volume units at standard or normal conditions per liquid volume units at flash conditions [START_REF] Mouchot | Scale and Corrosion control program, Example of two geothermal plants in Operation in the Upper Rhine Graben[END_REF][START_REF] Vandenberghe | The Merksplas-Beerse geothermal well (17W265) and the Dinantian reservoir[END_REF]. Finally, the gas phase composition is analyzed by chromatography [START_REF] Sanjuan | Major geochemical characteristics of geothermal brines from the Upper Rhine Graben granitic basement with constraints on temperature and circulation[END_REF], and these analysis results are associated to GLR data to ensure a further characterization of the geothermal fluid composition.

As the bubble point pressure is used to decide on the position of a submerged pump on the production well, its miscalculation leads to problems for geothermal plant design. Thus, any tool predicting geothermal fluid bubble point pressure should be accurate enough to avoid these complications. In order to ensure the calculations' accuracy, the model parameters must be fitted from experimental data obtained for thermodynamic systems representatives of the natural geothermal resources of interest. Experimental thermodynamic characterization of gas-brine systems is therefore essential to obtain the breakout pressure. In case of the Upper Rhine Graben, the geothermal brine composition is well-known [START_REF] Sanjuan | Major geochemical characteristics of geothermal brines from the Upper Rhine Graben granitic basement with constraints on temperature and circulation[END_REF]. This fluid is a very salty formation water (TDS: Total Dissolved Species around 90 g•L -1 ) which is registered as Na-Ca-Cl type and also contains dissolved gas. The off-gas obtained from the natural brines is composed of 86.2% carbon dioxide, with smaller quantities of nitrogen (N 2 , 9.85%) and methane (2.28%). In the literature, some studies have obtained experimental solubility data of gas mixtures in pure water: CO 2 + CH 4 + H 2 O [START_REF] Dhima | Solubility of hydrocarbons and CO 2 mixtures in water under high pressure[END_REF][START_REF] Qin | Experimental measurements of vapor-liquid equilibria of the H 2 O+ CO 2 + CH 4 ternary system[END_REF][START_REF] Hayashi | Direct measurement of solubility of methane and carbon dioxide in high pressure water using gas chromatography[END_REF][START_REF] Al Ghafri | Experimental and modeling study of the phase behavior of (methane+ CO 2 + water) mixtures[END_REF] or CO 2 + N 2 + H 2 O [START_REF] Foltran | Understanding the solubility of water in carbon capture and storage mixtures: An FTIR spectroscopic study of H 2 O+ CO 2 + N 2 ternary mixtures[END_REF][START_REF] Hassanpouryouzband | Solubility of flue gas or carbon dioxide-nitrogen gas mixtures in water and aqueous solutions of salts: Experimental measurement and thermodynamic modeling[END_REF]. Nevertheless, the studied gas mixtures were not in accordance with the Upper Rhine Graben representative gas mixture. In addition, these investigations were carried out in a salt-free aqueous solution. By considering only the main gas, carbon dioxide, a recent work on the CO 2 + H 2 O + Na + Ca + K + Cl system has shown that the minor species of the brine did not impact the carbon dioxide solubility [START_REF] Poulain | Experimental measurements of carbon dioxide solubility in Na-Ca-K-Cl solutions at high temperatures and pressures up to 20 MPa[END_REF]. The Upper Rhine Graben water could therefore be synthetized by water containing only sodium chloride (NaCl) and calcium chloride (CaCl 2 ).

From this information, the aim of this paper is to use a simple modelling approach to determine the bubble point pressure and GHG emissions in the Upper Rhine Graben and other geothermal sites worldwide. To perform this modeling work, the well-known Soreide and Whitson model with new optimized parameters will be used to estimate the downhole composition from the GLR, temperature, salinity (NaCl, CaCl 2 and KCl) and dry composition of the degassed fluid. A comparison with available industrial data of geothermal power plants will also be made.

Thermodynamic Modeling

The idea of the approach proposed in this section is to link the information easily available to geothermal operators, such as GLR, temperature and salinity, with the thermodynamics of phase equilibria in order to estimate equilibrium properties that are very important for the design, operation, and monitoring of geothermal wells. These properties mainly concern the prediction of the bubble-point pressure of gas-laden brines and the esti-mation of greenhouse gas emissions, and more specifically the amount of carbon dioxide and methane emitted in the process.

From GLR to the Gas-Laden Brine Composition

To calculate the bubble-point pressure, it is necessary to know the composition of all molecular species dissolved in the geothermal fluid. To do so, two types of approaches are possible: [START_REF] Poponi | Energy Technology Perspectives 2016: Towards Sustainable Urban Energy Systems[END_REF] an iterative approach to find the liquid composition by minimizing an objective function based on the differences between the calculated properties and the input information (GLR and gas composition, for instance). However, this can lead to optimization problems by converging on the local minimum instead of the global minimum. Hence, it requires a robust global optimization algorithm; (2) a simpler non-iterative approach based on a material balance by including the phases and compounds existing in the system. This work is based on the latter approach.

First, to link the GLR to the gas composition in geothermal fluids (brines), a material balance must be established on a flash system at a fixed temperature and pressure. Under certain conditions of temperature and pressure, a mixture of overall (monophasic) composition z i and molar flow rate F may split into several phases. This is called a flash or a phase split equilibrium. The most common case is the two-phase liquid-vapor separation (Figure 1), where a liquid phase of flow rate L and composition x i and a vapor phase of flow rate V and composition y i are formed. The different compounds i are distributed by number of moles in both phases according to their volatilities. The global and partial material balances are written as follows:

F = V + L (1)
Fz i = Vy i + Lx i = Vy i + (F -V)x i (2)
By defining the molar Gas-to-Liquid Ratio (GLR) mol = V(T,P) L(T,P) and the partition coefficient

K i = y i x i
, we obtain from Equation (1):

z i = y i (GLR) mol + 1 K i (GLR) mol + 1 (3)
From Equation (3), we can easily calculate the composition z i in the single-phase geothermal brine if we have the vapor composition y i , the GLR and the K i partition coefficient. In Equation (3), the composition y i is the wet steam composition. Since in most cases one only knows the dry steam composition, a calculation must be performed to obtain the wet steam composition by including the water content y w estimated by Equation ( 6) (or by an equation of state, see below) and by renormalizing the composition on all molecular species (water and gas). In Equation ( 4), f L i and f V i are the fugacities in the liquid (L) and vapor (V) phases, a w is the activity of water in brine calculated by an activity coefficient model (e.g., Pitzer model) and P sat w is the water vapor pressure.

f L i T, P, n L i = f V i T, P, n V i (4)
The determination of the GLR is performed at low pressure (typically atmospheric), and under these conditions, with the Poynting factor and the fugacity coefficients very close to 1, one can consider:

f L w(=water) T, P, n L w = a w × P sat w and f V w T, P, n V w = y w × P (5)
therefore:

y w = a w × P sat w (T) P (6)
Concerning the GLR, most of the data found in the literature are in the form of a volume ratio (Nm 3 /m 3 ), whereas the GLR in the above material balance equations is molar (mol/mol), which therefore requires the knowledge of the densities of the liquid and the vapor phases. The following Equation ( 7) can be used to convert the GLR from volume ratio to molar ratio.

(GLR) mol = (GLR) vol ρ Vap × Mw Liq ρ Liq × Mw Vap (7)
In this work, the density of the liquid phase (brine) ρ Liq was calculated by the model proposed by Al ghafri et al., 2012 [START_REF] Al Ghafri | Densities of aqueous MgCl 2 (aq), CaCl 2 (aq), KI (aq), NaCl (aq), KCl (aq), AlCl 3 (aq), and (0.964 NaCl+ 0.136 KCl)(aq) at temperatures between (283 and 472) K, pressures up to 68.5 MPa, and molalities up to 6 mol•kg-1[END_REF] and that of the gas phase ρ Vap by the multiparametric equation of state GERG-2008 ( Kunz and Wagner 2012 [41]). Concerning the GERG-2008 model, its Fortran code for calculating the thermophysical properties of gas mixtures is freely released by NIST on github [42].

Regarding the partition coefficient K i , correlations such as those of Giggenbach 1980 [START_REF] Giggenbach | Geothermal gas equilibria[END_REF] and D' Amore and Truesdell 1988 [44] can be used, but with prudence, since they only consider the temperature and not the effect of pressure and salinity. A better option is the one we use in this work which is to calculate the partition coefficient by the equation of state that will be used for the calculation of phase equilibrium and that will be presented below.

Equation of State (EoS)

Because of their simplicity and versatility, the most widely used thermodynamic models for determining phase diagrams are the cubic equations of state. However, among their weaknesses is the poor representation of the phase equilibria of mixtures containing water, especially the Type III phase diagrams according to the Van Konynenburg and Scott classification [START_REF] Van Konynenburg | Critical lines and phase equilibria in binary van der Waals mixtures[END_REF]. This very particular molecule (water) has a molecular structure and thermophysical properties very different from those of hydrocarbons, which completely changes the phase equilibria, and even more so in the case of the existence of dissolved salts in water. The need for a good representation of the phase equilibria of mixtures of gases with water or brine therefore requires modifications to the cubic equations of state, either by modifying some parameters (binary interaction coefficient Kij, cohesion parameter a, etc.) as in the Soreide and Whitson model [START_REF] Søreide | Peng-Robinson predictions for hydrocarbons, CO 2 , N 2 , and H 2 S with pure water and NaCl brine[END_REF] or by adding other terms to the equation of state to take into account the presence of water or saline water (CPA type models [START_REF] Chabab | Thermodynamic study of the CO 2 -H 2 O-NaCl system: Measurements of CO 2 solubility and modeling of phase equilibria using Soreide and Whitson, electrolyte CPA and SIT models[END_REF][START_REF] Olsen | Investigation of the performance of e-CPA for a wide range of properties for aqueous NaCl solutions[END_REF]).

The Soreide and Whitson (SW) EoS [START_REF] Søreide | Peng-Robinson predictions for hydrocarbons, CO 2 , N 2 , and H 2 S with pure water and NaCl brine[END_REF] was selected in this work since it is widely used for gas-brine systems and is implemented in several process modeling and reservoir simulator tools (Simulis Thermodynamics [START_REF] Baudouin | Simulis ® thermodynamics: An open framework for users and developers[END_REF], MUFITS [START_REF] Afanasyev | Compositional modeling of multicomponent gas injection into saline aquifers with the MUFITS simulator[END_REF], Eclipse 300 [START_REF] Schlumberger | PVTi and ECLIPSE 300: An Introduction to PVT Analysis and Compositional Simulation[END_REF], IHRRS [START_REF] Petitfrere | Three-phase EoS-based reservoir simulation with salinity dependent phaseequilibrium calculations[END_REF], etc.). It is an update of the well-known Peng Robinson equation of state (Equations ( 8)-( 13) by adding a NaCl salinity m NaCl (mol/kgw) dependence to some model parameters. The alpha function of water (Equation ( 14)) and the binary gas/water interaction parameters specific to each phase (aqueous k AQ ij and non-aqueous k N A ij , see Equations ( 15)-( 18)) are concerned. It makes the model "asymmetric" and adds more flexibility for the simultaneous representation of liquid and vapor phases by different parameters.

P(T, v) = RT v -b - a(T) v(v + b) + b(v -b) (8)
where P is the pressure, T the temperature, R the ideal gas constant and v the molar volume.

The energy parameter a and the molar co-volume b for a pure compound i are given by:

a i (T) = 0.457235529 R 2 T 2 c,i P c,i × α(T) (9) α(T) = 1 + 0.37464 + 1.54226ω i -0.26992ω i 2 1 - T T c,i 2 (10)
a m (T) = n c ∑ i=1 n c ∑ i=1 x i x j a i (T)a j (T) 1 -k ij ( 11
) b i = 0.0777960739 RT c,i P c,i (12)
b m = ∑ i x i b i (13)
where P c,i , T c,i and ω i are respectively the critical pressure, the critical temperature and the acentric factor of compound i.

α w (T, m s ) = 1 + 0.4530 1 - T T c,w 1 -0.0103m 1.1 NaCl + 0.0034 T T c,w -3 -1 2 (14)
a AQ m (T) =

n c ∑ i=1 n c ∑ i=1 x i x j a i (T)a j (T) 1 -k AQ ij ( 15
)
a N A m (T) = n c ∑ i=1 n c ∑ i=1 y i y j a i (T)a j (T) 1 -k N A ij (16)
The binary interaction parameters k AQ ij and k N A ij are respectively adjusted on the gas solubility data and water content data of gas + water + NaCl systems. The k N A ij parameter is generally constant or slightly dependent on temperature while the k AQ ij paramater depends on temperature and NaCl molality m NaCl . The binary gas/gas interaction parameters are those of the Peng-Robinson model and are widely available in the literature ( Whitson and Brule 2000 [START_REF] Whitson | Phase Behavior[END_REF], Lasala 2016 [START_REF] Lasala | Advanced Cubic Equations of State for Accurate Modelling of Fluid Mixtures[END_REF]. The gas/water parameters are listed in the "Results and Discussion" section.

Calculation Procedure

The procedure for calculating the bubble point pressure from the GLR, temperature and dry gas composition using an equation of state is shown in Figure 2. In the first step the dry gas composition is converted to the "true" composition (including water) and the GLR is converted to a molar ratio. Then, to estimate the composition of the geothermal fluid, the partition coefficients K i are generated by a classical PT-Flash calculation (by the Rachford-Rice method) [START_REF] Firoozabadi | Thermodynamics of Hydrocarbon Reservoirs[END_REF][START_REF] Michelsen | Thermodynamic Models: Fundamentals & Computational Aspects[END_REF][START_REF] Chabab | Thermodynamic Study of Complex Systems Containing Gas, Water, and Electrolytes: Application to Underground Gas Storage[END_REF] using the same temperature and pressure used to determine the GLR. However, the partition coefficients are very weakly dependent on the global composition z, so values of z i estimated close to those of the real case or simply values equal to 1/nc (relative to the number of compounds nc) can be used. Once the partition coefficients and the composition in the vapor phase y i are determined, the composition of the downhole fluid can be determined by Equation (3). Hence, the calculation of the bubble point could be started, knowing the required set of input data (temperature and liquid composition x i = z i ). It consists of an iterative search for the bubble point pressure, making sure each time that the vapor composition is stable, the sum of the gas molar fraction is equal to 1. The iterative part can be performed by the Newton method [START_REF] Michelsen | Thermodynamic Models: Fundamentals & Computational Aspects[END_REF] or by substitutions [START_REF] Sandler | Biochemical, and Engineering Thermodynamics[END_REF] (presented in Figure 2), which is simpler and sufficient. With this method, the bubble point pressure can be estimated at any point in the reservoir by knowing the temperature gradient, which will allow the determination of at which depth level the bubbling takes place. With this method, bubble point pressure can be estimated at any point at the surface or in the reservoir by knowing the temperature or its gradient. The GHG emission rate is directly related to the vaporized fraction β = V F which can be determined by a PT-Flash calculation at known temperature and pressure. Therefore, the GHG emission rate depends on the temperature and pressure conditions at the surface or in the tank and on the difference between the actual pressure and the bubble point pressure. In other words, the lower the pressure is below the bubble point pressure, the higher the vaporized fraction and the quantity of gas emitted. Generally, the flow rates F of geothermal fluids are given in terms of volume (e.g., m 3 /h), hence there is a need to know its density in order to be able to calculate the amount of gas emitted. The density of the liquid and vapor phases are given directly by the equation of state. This represents an advantage of using the phi-phi approach (use of an equation of state to describe both liquid and vapor phases) instead of the gamma-phi approach. The volumetric flow rate must be converted to a molar flow rate F using the molar density given by the equation of state at T, P and z i (Equation ( 3)), then multiplied by the vaporized fraction β to obtain the degassed quantity.

Results and Discussion

3.1. Mutual Solubilities of Gas-Water and Gas-Brine (Single or Mixed Salts) Systems a.

Water content in gas rich-phase

The asymmetry of the binary interaction parameters used in the Soreide and Whitson EoS permits one to overcome the lack of "flexibility" of the symmetrical (phi-phi) approach concerning the distinct representation of the liquid and vapor phases by using specific parameters for each phase. In addition to the improvement of the calculation of the liquid phase, the other advantage of this adaptation is a very good representation of the vapor phase in particular the water content in the gas-rich phase. In Figure 3, experimental literature data of water content in several gases of interest (CO 2 , CH 4 , N 2 , O 2 , and H 2 ) are compared with the SW model predictions using the optimized interaction parameters proposed in this work (see Section 3.1 d). As shown in the figure, the model accurately predicts the water content in these gases over a wide range of pressure and temperature, covering reservoir and surface conditions. The average absolute deviations between the predictions and the experimental literature data vary between 2 and 4% for the mentioned gases. CO 2 solubility in different brines (NaCl, CaCl 2 , or KCl) CO 2 is generally the major species among the gases dissolved in geothermal waters, and this is justified by its relatively high solubility compared to other gases. Therefore, in this section we focus on the characterization of its solubility under thermodynamic conditions in geothermal reservoirs. As previously mentioned, the model takes into account the salinity by including the molality of NaCl in the model parameters (alpha function of water and kij in the aqueous phase). Since sodium and chloride are the main species in natural brine, NaCl is the only salt considered in the original version of SW model. Considering NaCl as the only salt is a good first assumption to simplify the physicochemical calculations. However, in some geological reservoirs, other salts may be dissolved in the brine such as CaCl 2 , KCl, etc. These salts can exist in significant quantities and exhibit salting-out effects more or less important than NaCl for the same concentrations. Hence, it is important to consider these species for a more "realistic" representation of formation water. In this work, the binary interaction parameters used to model the CO 2 + H 2 O + NaCl system are those proposed by Chabab et al., 2019 [47] which improve very significantly the representation of the CO 2 solubility in particular in very concentrated and saturated NaCl-brines compared to the original parameters.

Regarding the modeling of CO 18)) was not modified for the other salts because: (1) there are no (or few) water content data of CO 2 in equilibrium with CaCl 2 and KCl brines; (2) the salting-out effect on water content is much less important than on gas solubility. Thus, for the salts CaCl 2 and KCl, the kij used are: Equation ( 17) for the aqueous phase with the coefficients listed in Table 1 and Equation ( 18) for the non-aqueous phase. In Figure 4, the experimental literature data of CO 2 solubilities in the different salts (NaCl, CaCl 2 and KCl) are compared to the calculations generated by the SW model using the proposed optimized parameters. As shown in the figure, the model captures well the salting-out effect of the different salts on the solubility of CO 2 under the different conditions of temperature, pressure and molality of salts. The model reproduces the solubility in NaCl and CaCl 2 solutions perfectly, but a little less well in KCl solution. This could be due to several factors (data consistency, unsuitable kij expression, optimization problem (local minimum), etc.), but is probably due to the consistency of the literature data, knowing that there is a mismatch between the data of Hou et al., 2013 [START_REF] Hou | Phase equilibria of (CO 2 + H 2 O+ NaCl) and (CO 2 + H 2 O+ KCl): Measurements and modeling[END_REF] and those of Zhao et al., 2015 [START_REF] Zhao | Experimental studies and modeling of CO 2 solubility in high temperature aqueous CaCl 2 , MgCl 2 , Na 2 SO 4 , and KCl solutions[END_REF], and that the model reproduces very well the data of Hou et al. (see Table 2). 2.97 a : using the optimized parameters proposed in this work (Table 1).

k AQ CO 2 ,H 2 O = T T c,CO 2 a + b T T c,CO 2 + c T T c,CO 2 m salt + m 2 salt d + e T T c,CO 2 + f (17) k N A CO 2 ,H 2 O = 0.68208385571 × 10 -3 T -2.066623464504 × 10 -2 (18)
In order to make the SW model capable of predicting the solubility of CO 2 in real brines (mixed salt systems), correlations (Duhring lines) between the molality of the salts (CaCl 2 and KCl) and their NaCl equivalents were developed in this work. This was undertaken by an iterative (Newton-Raphson) approach of finding equal CO 2 solubilities in the different salts considering NaCl as reference salt ( Markham and Kobe 1941 [76], Ratnakar et al., 2020 [START_REF] Ratnakar | On the prediction of gas solubility in brine solutions with single or mixed salts: Applications to gas injection and CO 2 capture/sequestration[END_REF]). The data resulting from the search for equal solubilities have been correlated in the following equations (Equations ( 19) and ( 20)) to make the conversion from salt molality to NaCl equivalent molality straightforward. The correspondence between salts and NaCl is independent of pressure but depends on temperature. Figure 5 shows the Duhring lines using the salt→NaCl equivalence correlations developed to estimate the NaCl equivalent of CaCl 2 and KCl at equal CO 2 solubilities for two isotherms (323 and 373 K).

m NaCl = m CaCl 2 6.4824099 - In Figure 6, the solubility of CO 2 in CaCl 2 and KCl solutions calculated by the SW model using the Kij corresponding to each salt (Table 1) were compared to the calculations using their NaCl equivalent (with the help of Equations ( 19) and ( 20)) and the NaCl coefficients for the binary interaction parameter (Table 1). The results show a perfect adequacy between the calculations by the "apparent" molality and those by the "real" Kij corresponding to each salt. Thanks to the salt→NaCl equivalence correlations proposed in this work, the SW model can be applied to the study of gas solubility in mixed-salt systems, which is more realistic than considering NaCl as the only dissolved salt in groundwater. The idea of this approach is to convert with the help of these correlations the molalities of the other salts (CaCl 2 and KCl) into equivalent NaCl molality, add them to the NaCl molality, and this total is considered in the calculations of the water alpha function and the binary interaction parameters. To test this approach, the prediction results for the solubility of CO 2 in a salt mixture (CaCl 2 + NaCl) using NaCl equivalence were compared with predictions without including this molality correction. This comparison is illustrated in Figure 7 taking as reference the experimental data recently measured by Poulain et al., 2019 [39] and Lara Cruz et al., 2021 [START_REF] Cruz | Experimental study of carbon dioxide solubility in sodium chloride and calcium chloride brines at 333.15 and 453.15 K for pressures up to 40 MPa[END_REF], and this corresponds to five isotherms between 323 and 453 K at fixed salt molality (1.2 mol/kgw NaCl + 0.2 mol/kgw CaCl 2 ) and for pressures up to 40 MPa. As shown in Figure 7, the inclusion of the salt equivalence correlation improves the prediction of CO 2 solubility in mixed-salt solution (NaCl + CaCl 2 ) and the enhancement becomes clearer and more important with increasing pressure (Figure 7a). Finally, these equivalence relations based on the salting-out effects of the different salts improves the prediction in multi-salt systems and especially when other salts exist in non-negligible quantities in the studied environments. 3. It should be noted that the calculations by the SW model are very sensitive to the critical data (critical temperature and pressure and acentric factor) used, so it is important to use the appropriate critical data provided with the binary interaction parameters.

T T c,CH 4 2 1 -4.198569 × 10 -2 m NaCl k N A CH 4 ,H 2 O = 0.494435 O 2 T c,O 2 = 154.581 K P c,O 2 = 50.43 bar ω O 2 = 0.0222 k AQ O 2 ,H 2 O = -1.167744 1 + 3.361921 × 10 -2 m NaCl 0.8 +0.466607 T T c,O 2 1 + 8.457306 × 10 -2 m NaCl 0.8 k N A O 2 ,H 2 O = 0.581650 N 2 T c,N 2 = 126.19 K P c,N 2 = 33.958 bar ω N 2 = 0.0372 k AQ N 2 ,H 2 O = -1.709096 1 + 1.792130 × 10 -2 m NaCl 0.8 +0.450487 T T c,N 2 1 + 0.066426m NaCl 0.8 k N A N 2 ,H 2 O = 0.385438 H 2 T c,H 2 = 33.145 K P c,H 2 = 12.964 bar ω H 2 = -0.219 k AQ H 2 ,H 2 O = -2.34 1 + 3.88 × 10 -3 m NaCl 0.443 + 0.166 T T c,H 2 1 + 0.049m NaCl 0.799 + -12.69 exp -0.474 T T c,H 2 k N A H 2 ,H 2 O = -0.3776 + 0.08385 T T c,H 2
In Figure 8, the prediction of the solubility of CH 4 , N 2 , O 2 , and H 2 in water and NaCl brine is compared with the literature's experimental data. As shown in this figure, the predictions are in very good agreement with the experimental data. Indeed, the use of the Søreide-Whitson model with the proposed optimized parameters allows the estimation of the solubility of these gases in brine over a wide range of temperature, pressure and salinity and with an average absolute deviation of about 3%. The difference in the studied gas solubilities can be explained by several factors such as the size of the molecules [START_REF] Battino | The O 2 /N 2 ratio gas solubility mystery[END_REF]. Indeed, large molecules need more energy to form a cavity in the aqueous phase. Moreover, this effect is important at high pressure. As such, it can explain the strong increase in the solubility of H 2 at elevated pressure (see Figure 8), which can even exceed the solubility of other gases (O 2 and CH 4 ). The contribution of hydrogen bonds (solvation) to solubility is very important and explains the high dissolution of CO 2 (which is a Lewis acid) [START_REF] Kontogeorgis | Solvation phenomena in association theories with applications to oil & gas and chemical industries[END_REF] compared to other gases. Another factor to explain the variation in solubility could be the polarizability, but it is less dominant (in water) compared to hydrogen bonding and coulombic interactions [START_REF] Campanell | On the role of solute polarizability in determining the solubilities of gases in liquids[END_REF]. Finally, the notion of solubility is influenced by several parameters which can also have counterbalances of these different dependencies. Hence, the comprehension of the difference between gas solubility could be difficult to explain directly. e.

Mixed-gas cosolubility: CO 2 + CH 4 + H 2 O and CO 2 + N 2 + H 2 O systems For mixed gas cosolubility modeling, thanks to the mixing rules, the model is easily applicable to multicomponent mixtures when binary interaction parameters are available. In this work, we are interested in gas mixtures including mainly CO 2 (70-99%) with other gases in more or less important fractions such as CH 4 and N 2 . Concerning the CO 2 + N 2 + H 2 O ternary mixture, the only solubility data available in the literature (Liu et al., 2012 [84]) are only at low temperature (up to 318 K) which have been studied for other applications such as the transport of flue gases. In Figure 9, the cosolubility data of the mixture of 90%CO 2 + 10%N 2 in water at 318 K are compared to the predictions of the SW model. In general, the predictions are good especially for N 2 but less good for CO 2 since the SW model underestimates the solubility of CO 2 at low temperatures (see Chabab et al., 2019 [47]). This is explained by the fact that the model parameters were fitted to data at higher temperatures and which are more representative of those found in geological reservoirs. Concerning the CO 2 + CH 4 + H 2 O ternary mixture, there are some experimental data available in the literature [START_REF] Dhima | Solubility of hydrocarbons and CO 2 mixtures in water under high pressure[END_REF][START_REF] Qin | Experimental measurements of vapor-liquid equilibria of the H 2 O+ CO 2 + CH 4 ternary system[END_REF] and especially in the temperature and pressure range representative of deep geological formations. In Figure 10, the cosolubility data of 72%CO 2 + 28%CH 4 in water at 376 K are compared to the predictions of the SW model. The predictions are very good for both CO 2 and CH 4 over a very wide pressure range (up to 50 MPa). It should be noted that these results were obtained without any additional parameters and without the inclusion of these data in the adjustment of the binary interaction parameters. This illustrates the advantage of the "predictive" aspect of the equations of state as opposed to the empirical correlations/models. Finally, knowing the temperature, pressure, dry gas composition and salt concentrations, the phase equilibria of the gas(s)/brine(s) systems can be predicted with good accuracy with the SW model using the interaction parameters and salt equivalence correlations proposed in this work. Moreover, any other types of phase equilibrium calculations are applicable to this model such as the search for the bubble point of gas-laden geothermal brines, which is the subject of the next section.

Application: Geothermal Fields Operations

In this part, the approach proposed in Section 2 is applied to characterize the degassing of geothermal fluids. The bubble pressure of the gas-laden brines and the GHG emission rate are estimated from information available to the geothermal operators such as the gas-to-liquid ratio (GLR), the temperature (T) in the well and the dry gas compositions of the degassed compounds.

Estimation of the Bubble-Point Pressure from the GLR

In order to validate the developed approach, this input information (GLR, T, salinity and gas composition) and their corresponding bubble-point pressures must be known. However, in the literature, these data are not widely available and generally incomplete. Two data sources [START_REF] Vandenberghe | The Merksplas-Beerse geothermal well (17W265) and the Dinantian reservoir[END_REF][START_REF] Ungemach | Handling of corrosion and scaling shortcommings in low enthalpy geothermal environments[END_REF] providing some of this necessary information have been considered for comparison and validation of the proposed predictive approach. The data from the Paris area reported by Ungemach [START_REF] Ungemach | Handling of corrosion and scaling shortcommings in low enthalpy geothermal environments[END_REF] contain more information but the salinity of the brine samples is not given, so we made the calculations at different equivalent NaCl molalities. In Table 4 and Figures 11 and12, the bubble pressure data reported from PVT analyses of samples taken from the geothermal wells are listed, illustrated and compared to those predicted by the model. Considering a NaCl molality of 0.5 mol/kgw, in Figure 11, the model reproduces well the Ungemach data with an AAD% of 5.7% for GLR between 0.23 and 0.28 and temperatures between 333 and 337 K knowing that the uncertainties of these data are not reported either. For samples taken at the same depth (see Table 4), in contrast to the Ungemach data (Paris) [START_REF] Ungemach | Handling of corrosion and scaling shortcommings in low enthalpy geothermal environments[END_REF], the Vandenberghe et al. (North Belgium) [START_REF] Vandenberghe | The Merksplas-Beerse geothermal well (17W265) and the Dinantian reservoir[END_REF] data appear to have very large uncertainty (poor repeatability) with large differences in gas composition, GLR and inhomogeneity of salt content in the analyzed brine samples (from the same site) and consequently wide variation in bubble pressure. For this reason, larger deviations between the model predictions and the Vandenberghe data are observed. However, the model captures qualitatively very well the variation of the bubble point as a function of GLR and gas composition and gives coherent values with three out of five points among these data (see Table 4 and Figure 12). [START_REF] Mouchot | First year of operation from EGS geothermal plants in Alsace, France: Scaling issues[END_REF]) from a few geothermal wells operated in the Upper Rhine Graben, the presented approach was used to predict the bubble point pressure and GHG emission rates in these sites and more precisely in the production wells GPK2 (Soultz-sous-Forêts) and GRT2 (Rittershoffen). These data are listed in Table 5 along with the predicted bubble pressure values and flow rates of the degassed fluids at different temperatures and pressures which are further detailed in Table 6. According to Mouchot et al., at the geothermal fluid temperature of 150-170 • C, a pressure of 2.3-2.5 MPa is maintained at the surface to avoid degassing and related problems such as corrosion and scale formation. This pressure, determined just by feedback, represents a pressure around the bubble point to keep the geothermal fluid in a liquid state. At a NaCl molality of 1.3-1.4 mol/kgw (obtained by conversion using the equivalence correlation proposed in this work), the model predicts bubble pressures of 2.296 and 2.59 MPa for both sites, which are very close to those reported (2.3 and 2.5 MPa) and used as lower limits in surface facilities. As a result of this excellent bubble pressure prediction, the method can also be considered for estimating the GHG rates that can be emitted and which are directly related to the vaporized fraction.

The estimates of GHG emissions expressed in grams of CO 2 equivalent per kWh of energy produced are listed in Table 6. They have been determined from the energy production of these facilities and the quantities of degassed fluids, notably CO 2 and CH 4 , considering that the global warming potential of CH 4 is 25 times higher than that of CO 2 . As shown in Table 6, by decreasing the pressure and/or increasing the temperature, the degassing and the GHG emissions become important. Compared to the average of 122 gCO 2 /kWh (see Introduction), the gas emissions are low, and may be zero if the pressure is maintained above the bubble pressure. Even at high depressurization, the emissions at these two sites remain low compared to the average emissions from geothermal sites and even lower than those from fossil fuel-based plants. This conclusion is in line with that of Lu et al. [START_REF] Lu | A global review of enhanced geothermal system (EGS)[END_REF] and Pratiwi et al. [START_REF] Pratiwi | Life-cycle climate-change impact assessment of enhanced geothermal system plants in the Upper Rhine Valley[END_REF] on the fact that the Rhine Graben sites are environmentally friendly and indeed offer low-carbon geothermal energy.

Conclusions

In the context of geothermal energy exploitation, the knowledge of the bubble pressure of geothermal fluids is very important to avoid technical and environmental problems as well as for the design, optimization and monitoring of surface facilities. In this work, a simple and efficient procedure for the prediction of phase equilibria (bubble point and flash calculation) of gas/brine systems during geothermal fluid exploitation is proposed. Being the main factors involved in the degassing problem, water content in the gas-rich phase and gas solubility have been successfully modeled in different brines at conditions covering widely the temperature, pressure and salinity range of geothermal applications. The model reproduces accurately the data reported by industrial operators and can therefore be used to estimate the equilibrium properties of any gas-laden brine from some readily available information such as temperature, GLR and dry gas composition. Finally, to fill the great lack of equilibrium data (solubility, GLR and bubble pressure) of gas/brine systems, an experimental pilot is under development, which will allow us to generate data that can be compared and used to improve the model presented in this work.

Figure 1 .

 1 Figure 1. Phase-splitting at constant temperature and pressure (PT-Flash).

Figure 2 .

 2 Figure 2. Algorithm for calculating the bubble point pressure of geothermal fluids.

Figure 3 .

 3 Figure 3. Water content in gas (CO 2 , CH 4 , N 2 , O 2 and H 2 )-rich phase: Comparison of SW EoS predictions (solid lines) with experimental literature data (symbols) [59-64].

Figure 4 .

 4 Figure 4. CO 2 solubility in different brines (NaCl, CaCl 2 , or KCl): comparison of predictions using SW EoS (solid lines) with literature data (symbols) [63,64,66-73].

3 ( 20 )Figure 5 .

 3205 Figure 5. Duhring lines of equal CO 2 solubility in different brines (NaCl as the reference salt).

  in mixed-salt brine (NaCl + CaCl 2 + H 2 O)

Figure 6 .

 6 Figure 6. Salt-dependent Kij vs. equivalent NaCl molality using SW EoS.

Figure 7 .

 7 Figure 7. Prediction of CO 2 solubility in mixed-salt brine (NaCl + CaCl 2 + H 2 O) with (solid lines) and without (dashed lines) salt molality correction. Data (symbols) from: (a) [75]; (b) [39].

  d. CH 4 , N 2 , O 2 , and H 2 solubility in water and NaCl brine Besides the CO 2 presence, other gases such as CH 4 , N 2 , O 2 and H 2 are generally dissolved in geothermal brines in varying amounts and must be considered as they significantly change the bubble pressure of gas-laden brines. In their original paper, Soreide and Whitson have also proposed binary interaction parameters for CH 4 + H 2 O and N 2 + H 2 O. However, they have not proposed parameters for the H 2 + H 2 O and O 2 + H 2 O binaries. The parameters used for these two binaries are those recently proposed in the works of Chabab et al. [78,79]. On the other hand, the coefficients of the Kij correlations proposed by Soreide and Whitson for CH 4 + H 2 O and N 2 + H 2 O have been readjusted in this work. The expressions of the Kij and their coefficients are presented in Table

Gas-H 2 O

 2 Interaction Parameters in Aqueous and Non-Aqueous Phases CH 4 T c,CH 4 = 190.56 K P c,CH 4 = 45.992 bar ω CH 4 = 0.01142 k AQ CH 4 ,H 2 O = -1.625685 1 + 8.590105 × 10 -21 m NaCl +1.114873 T T c,CH 4 1 + 1.812763 × 10 -3 m NaCl -0.169968

Figure 8 .

 8 Figure 8. CH 4 , N 2 , O 2 and H 2 solubility in water and NaCl brine at 323 ± 1 K: Comparison of SW EoS predictions (solid lines) with experimental literature data (symbols) [78-80].

Figure 9 .

 9 Figure 9. Mixed-gas cosolubility in CO 2 + N 2 + H 2 O system: Comparison of SW EoS predictions (solid lines) with experimental literature data (symbols) [84].

Figure 10 .

 10 Figure 10. Mixed-gas cosolubility in CO 2 + CH 4 + H 2 O system: Comparison of SW EoS predictions (solid lines) with experimental literature data (symbols) [34].

  2 + H 2 O + CaCl 2 and CO 2 + H 2 O + KCl systems, the following assumptions were considered to extend the SW model to these salt solutions:

	-	The coefficients of the CO 2 -H 2 O binary interaction parameter expression in the aque-
		ous phase k AQ CO 2 ,H 2 O (Equation (17)) were readjusted on the CO 2 solubility data in the mentioned brines (H 2 O + CaCl 2 and H 2 O + KCl). The readjusted Kij coefficients are
		listed in Table 1.
	-	The CO 2 -H 2 O binary interaction parameter in the non-aqueous phase k N A CO 2 ,H 2 O (Equa-tion (

-

The original alpha function of water was kept as a function of the NaCl molality and was not modified for the other salts for the following reasons: (1) To propose a salt equivalence technique considering NaCl as the reference salt (see hereafter); (2) the alpha function proposed by SW is obtained by fitting on vapor pressure data of the H 2 O + NaCl system and the difference between the vapor pressures of the brines of different salts does not have significant effects on the representation of the phase equilibria especially at the thermodynamic conditions of interest (high pressure).

Table 1 .

 1 Coefficients of the k

		AQ CO 2 ,H 2 O correlation (Equation (17)).
		NaCl	CaCl 2	KCl
	a	0.43575155	0.478902106	0.44269716
	b	-5.766906744 × 10 -2	-7.35093213 × 10 -2	-5.746192569 × 10 -2
	c	8.26464849 × 10 -3	2.05379759 × 10 -2	6.655117 × 10 -3
	d	1.29539193 × 10 -3	2.04457138 × 10 -3	-1.318885435 × 10 -3
	e	-1.6698848 × 10 -3	-3.22673398 × 10 -3	6.161936936 × 10 -4
	f	-0.47866096	-0.50876635	-0.488832169

Table 2 .

 2 Deviation between experimental data and SW model prediction of CO 2 solubility in brine.

	Data Reference	AAD% (SW Model) a

  3.808876 × 10 -3 × T -923.14 × NaCl = 0.0095690161 × T × m KCl + 0.0004886086 × T × m KCl 2 -1.1835076 × m KCl -0.16724994 × m KCl 2 -1.170219186 × 10 -5 × m KCl × T 2 -21.4278766 × 0.00467461197 m KCl × m KCl

	-	11×m CaCl 2 +0.0262706×T-6.739781 5.0858999×m CaCl 2	m CaCl 2 T	(19)

m

Table 3 .

 3 Gas-H 2 O interaction parameters in aqueous and non-aqueous phases for the SW EoS.

Table 4 .

 4 Comparison of predicted bubble pressures (Pb) with those obtained from PVT data of downhole samples[START_REF] Vandenberghe | The Merksplas-Beerse geothermal well (17W265) and the Dinantian reservoir[END_REF][START_REF] Ungemach | Handling of corrosion and scaling shortcommings in low enthalpy geothermal environments[END_REF]. .2.2. Estimation of Bubble-Point Pressure and GHG Emission Rate in the Upper Rhine Graben Geothermal Site Based on recently published data[START_REF] Mouchot | First year of operation from EGS geothermal plants in Alsace, France: Scaling issues[END_REF]

	Data Source	Depth (m bgl)	CO 2 (% mol)	CH 4 (% mol)	N 2 (% mol)	O 2 (% mol)	GLR (m 3 /m 3 )	T Exp (K)	Pb Exp (MPa)	NaCl Molality (mol/kgw)	Pb Calc (MPa)	AD% Pb
		1646	0.9367	0.0340	0.0271	0.0022	1.62	347.05	2.654	2.5	2.470	6.95
	Vandenberghe et al. (2001) [31]	1646 1740 1740	0.9225 0.8311 0.9243	0.0321 0.0627 0.0354	0.0443 0.1015 0.0323	0.0011 0.0047 0.0079	1.11 0.96 0.71	347.05 347.05 347.05	1.482 2.654 1.276	2.5 2.5 2.5	2.074 3.158 1.331	39.93 18.95 4.38
		1646	0.8795	0.0443	0.0751	0.0011	1.22	347.05	2.896	2.5	3.107	7.28
		160	0.5241	0.2113	0.2646	-	0.23	335.35	0.790	0.5	0.737	6.77
	Ungemach, P.	1818	0.5353	0.2006	0.2641	-	0.25	337.35	0.820	0.5	0.767	6.41
	(2001) [85]	1818	0.5245	0.1937	0.2818	-	0.26	337.35	0.840	0.5	0.820	2.41
		160	0.5307	0.1994	0.2699	-	0.28	333.35	0.830	0.5	0.890	7.24

Figure 11. Comparison of the predicted (lines) bubble pressures with those from literature (red symbols) [85] at different values of GLR and NaCl molality.

Figure 12. Comparison of the predicted (black symbols) bubble pressures with those from literature (red symbols) [31] at different values of GLR and at 2.5 mol/kgw of NaCl.

3

Table 5 .

 5 Prediction of bubble pressures (Pb) and GHG emission rates in the Upper Rhine Graben geothermal sites. Data from[START_REF] Mouchot | First year of operation from EGS geothermal plants in Alsace, France: Scaling issues[END_REF] [START_REF] Mouchot | First year of operation from EGS geothermal plants in Alsace, France: Scaling issues[END_REF]. b Equivalent salinity in NaCl (mNaCl = 1.2 and mCaCl 2 = 0.2 mol/kgw) calculated by the correlation developed in this work (Equation (19)). c Predicted using the developed model.

	Geothermal Site	Gas Composition at Production Well a CO 2 (% mol) a CH 4 (% mol) a N 2 (% mol) a	GLR (Nm 3 /m 3 ) at 0 • C and 1 atm a	T (K)	Pb (MPa) Exp at T a	NaCl Molality (mol/kgw) b	Pb Calc (MPa) c	P (MPa)	Geothermal Brine Flow Rate (kg/s) a	Degassed Amount (kg/h) c
		0.91	0.02	0.07	1.03	423.15	≈ 2.3	1.426	2.296	≥2.296	30	0
	Soultz-sous-Forêts	0.91 0.91	0.02 0.02	0.07 0.07	1.03 1.03	423.15 423.15	N.A. N.A.	1.426 1.426	2.296 2.296	2.100 1.900	30 30	7.7763 18.6904
		0.91	0.02	0.07	1.03	423.15	N.A.	1.426	2.296	1.5	30	62.6332
		0.91	0.02	0.07	1.2	443.15	≈ 2.5	1.3	2.59	≥2.59	72.5	0
	Rittershoffen	0.91 0.91	0.02 0.02	0.07 0.07	1.2 1.2	443.15 443.15	N.A. N.A.	1.3 1.3	2.59 2.59	2.100 1.900	72.5 72.5	84.3599 154.0723
		0.91	0.02	0.07	1.2	443.15	N.A.	1.3	2.59	1.5	72.5	425.6754

a

Table 6 .

 6 Prediction of GHG emission rates in the Upper Rhine Graben geothermal sites at different pressures and temperatures: estimated amounts and compositions of emitted gases.

					Geother-mal	Geother-mal									Degassed
	Geothermal Site	T (K)	P (MPa)	β= V F (Molar)	Brine Flow Rate F (kg/s) a	Brine Flow Rate F (mol/s) b	y H 2 O	b	y CO 2	b	y CH 4	b	y N 2	b	Amount V (mol/s)

b Degassed Amount V (kg/s) b GHG Emitted (g-eqCO 2 /kWh) b

	Soultz-sous-Forêts (12 GWh e /y)	423.15 423.15 423.15	2.10 1.90 1.50	4.01 × 10 -5 9.53 × 10 -5 3.16 × 10 -4	30 30 30	1662.561 1662.561 1662.561	0.23081 0.25370 0.31739	0.45380 0.48601 0.53799	0.04583 0.04174 0.02783	0.26956 0.21855 0.11679	0.0666364 0.1585203 0.5254182	0.0021601 0.0051918 0.0173981	6.71 15.86 48.06
	Rittershoffen (180 GWh th /y)	443.15 443.15 443.15	2.10 1.90 1.50	1.89 × 10 -4 3.47 × 10 -4 1.00 × 10 -3	72.5 72.5 72.5	4017.522 4017.522 4017.522	0.38136 0.41889 0.52317	0.43862 0.44965 0.41666	0.02994 0.02408 0.01255	0.15008 0.10738 0.04762	0.7593836 1.3927135 4.0359467	0.0234333 0.04279786 0.11824318	4.16 7.18 16.51

a Data from

[START_REF] Mouchot | First year of operation from EGS geothermal plants in Alsace, France: Scaling issues[END_REF]

[START_REF] Mouchot | First year of operation from EGS geothermal plants in Alsace, France: Scaling issues[END_REF]

. b Predicted using the developed model.

Acknowledgments: S.C. is funded through the ISITE programme E2S, supported by ANR PIA and Région Nouvelle-Aquitaine. Authors also thank two anonymous reviewers for insightful comments that improved the manuscript and Challote Cheng and the other members, for their editorial work.

Funding: This research was funded by the Région Nouvelle-Aquitaine, and the Agence Nationale de la Recherche (Projet Investissement d'Avenir (programme E2S)).

Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.

Data Availability Statement: Not applicable.

Conflicts of Interest:

The authors declare no conflict of interest.