Issue
Sci. Tech. Energ. Transition
Volume 78, 2023
Selected Papers from First European Conference on Gas Hydrates (ECGH), 2022
Article Number 8
Number of page(s) 11
DOI https://doi.org/10.2516/stet/2023005
Published online 04 April 2023

© The Author(s), published by EDP Sciences, 2023

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

CO2 injection in geological formations is a key aspect of the Carbon capture and Storage (CCS) strategy. The captured CO2 will be injected into depleted oil and gas reservoirs and/or saline aquifers. The injection and potential CO2 expansion can be associated with a large Joule–Thomson effect and may in turn pose a risk of hydrate or ice formation in the formation due to the cooling effect [1]. In addition, CO2 can be injected/stored through Carbonated Water Injection technology [2], where CO2 from offshore power production is captured from flue gas and compressed/injected into the geological formation as dissolved gas in the injection water. Hydrate could be formed in the injection water with dissolved CO2.

Although different thermodynamic models have been developed to predict the carbon dioxide hydrate stability in aqueous NaCl solutions, high-accuracy experimental data are always necessary to validate or improve these predictive tools. Experimental data for carbon dioxide hydrates in equilibrium with sodium chloride solutions have been measured and reported by various authors in different hydrate regions. Table 1 gives a list of these data, reporting the temperature range and source of the experimental data. Available experimental data in the literature for CO2 in NaCl aqueous solutions have mostly been measured at low pressures below the saturation lines. As seen in the table, data for high salinity are also limited.

Table 1

Literature data for hydrate equilibria of CO2 in presence of NaCl solutions (VLwH: Vapour – Liquid (aqueous) – Hydrate / LLwH: Liquid – Liquid (aqueous) – Hydrate).

In this study, experimental investigations were conducted to measure the hydrate stability zone of CO2 in NaCl solutions over a wide range of conditions, in particular, to fill the gaps above the liquid locus and at high salinity. Moreover, experimental data on the hydrate equilibrium temperature of aqueous phases undersaturated in CO2 with different CO2 saturations and different aqueous phase compositions were carried out in order to validate and complement the available literature data (Tab. 2).

Table 2

Literature data of CO2 solubility in water and brines at hydrate–liquid water (HLw) equilibrium.

The Cubic-Plus-Association (CPA-EoS) Equation of State combined with the solid solution theory of van der Waals and Platteeuw [27] as developed by Parrish and Prausnitz [28] was employed to model the fluid and hydrate phase equilibria as previously described by Chapoy et al. [2932]. The predictions of the thermodynamic model were compared with the newly experimentally measured properties.

2 Experimental section

2.1 Materials

The CO2 used was supplied by Air Products and has a stated purity of 99.995%. The NaCl and CaCl2 were both supplied by Sigma–Aldrich and have a stated purity of 99.5%. Distilled water was used in all tests. Details of the used fluids are given in Table 3.

Table 3

Composition of the chemical used in this work.

2.2 Experimental apparatus

The tests were conducted using a fixed-volume, temperature-controlled rocking cell as described below. The experimental set-up is comprised of an equilibrium cell (Fig. 1), cryostat, rocking/pivot mechanism, and temperature/pressure recording equipment controlled by a PC. The equilibrium cell is a 300 mL titanium cylindrical pressure vessel with a mixing ball, mounted on a horizontal pivot with an associated stand for pneumatically controlled rocking through 180°. The rocking of the cell, and the subsequent movement of the mixing ball within it, ensures adequate mixing of the cell fluids.

thumbnail Fig. 1

Schematic of the setup used to measure the hydrate dissociation condition

The rig has a working temperature range of 203–323 K, with a maximum operating pressure of 69 MPa. Circulating coolant from a cryostat within a jacket surrounding the cell controls the cell temperature. The cryostat can maintain the cell temperature stability to within better than 0.05 K. To achieve good temperature stability, the jacket is insulated with polystyrene board while connecting pipe work is covered with plastic foam.

The temperature is measured and monitored using a PRT (Platinum Resistance Thermometer) located within the cooling jacket of the cell. The cell temperature is measured with an uncertainty of 0.1 K. A strain gauge pressure transducer with an uncertainty of 0.03 MPa is used to monitor pressure. Temperature and Pressure are monitored and recorded by a PC through a RS 232 serial port. The cryostat can be monitored and controlled via an interface connected to a serial port on the computer. The temperature probe and pressure transducer are regularly calibrated using NAMAS certified temperature probe and a Budenberg dead weight tester, respectively.

2.3 Experimental methods

Hydrate dissociation point measurements for saturated systems were conducted using a reliable isochoric step-heating method as described in our previous papers [33, 34].

For undersaturated systems mixtures of CO2 and aqueous phase, composed of either distilled water or water with salts were made in a separate 600 mL piston vessel. In the first instance, a measured amount of CO2 was injected into the clean evacuated vessel. The amount of CO2 was measured using a balance with a precision of 0.01 g. The aqueous phase was then pumped into the vessel using a Quizix high-pressure pump. The volume of injected liquid was determined according to the density of the aqueous phase and the desired CO2 concentration. The weight of the aqueous phase injected into the vessel was checked using the balance. This method gave an uncertainty of ±0.05 weight % for the mixtures used in this study. The pressure of the mixture was increased by injecting nitrogen into the vessel on the opposite side of the piston to the test fluid. The pressure was increased at least 14 MPa above the expected bubble point of the mixture and was mixed for at least 30 min to ensure it was a single phase. The equilibrated mixture was then transferred to the rocking cell. During the transfer, the pressure of the mixture in the 600 mL cell and the test fluid pressure was maintained at a pressure at least 14 MPa above the expected bubble point to ensure it was always a single phase.

Once the test fluids had been loaded into the test cell the pressure and temperature were adjusted in order to make a HET measurement. The pressure was adjusted by injecting or withdrawing equilibrated fluid into the vessel. The cell was then rocked to achieve equilibrium before lowering the temperature stepwise below the expected HET temperature. The presence of hydrates was indicated by a change in the P versus T trend as shown in Figure 2. The temperature was then increased stepwise allowing sufficient time for equilibrium to be achieved at each step. The equilibrium P/T data were then plotted to determine the point at which hydrates are no longer present as indicated by a change in the P versus T plot. An example is shown in Figure 2 below.

thumbnail Fig. 2

Plot of equilibrium P/T data recorded during test xCO2 = 2.11 × 10−2 mole fraction in distilled water.

3 Thermodynamic modeling

A general phase equilibrium model based on the uniformity of component fugacities in all phases is used to predict phase equilibria, water activity, and the hydrate forming conditions. A description of the thermodynamic model and parameters can be found elsewhere [2932]. The model has been validated against carbon dioxide hydrate dissociation data [35] and gas solubility in brines at temperatures above the carbon dioxide critical temperature [36]. The model can also predict the gas solubility at low temperatures (T < 283.15 K) with an excellent agreement as seen in Figure 3.

thumbnail Fig. 3

Predicted and experimental CO2 solubility in aqueous sodium chloride. Black lines: predictions using the CPA-EoS model. : 5.52 wt% NaCl [24]; : 14.92 wt% NaCl [24]; : 22.61 wt% NaCl [24]; : 10 wt% NaCl (Heriot–Watt unpublished data); : 20 wt% NaCl (Heriot–Watt unpublished data).

4 Results and discussions

The experimental hydrate dissociations for CO2 in equilibrium with 5, 10, 15, 20, 23, and 25 wt% sodium chloride solutions are reported in Table 4 and plotted in Figures 4 and 5. In the vapor region (Fig. 4) the model is in excellent agreement with our new experimental data and most of the available literature data, with typical agreement better than 0.5 K. For the 25 wt% solutions, the problem is more complex as when the system is cooled down the system forms NaCl·2H2O before clathrate hydrates and the concentration in salt will therefore decrease in the remaining aqueous phase until it reaches the eutectic concentration 23.2 wt%, nevertheless taking into account the salting-out effect the model can still predict the dissociation conditions. In the liquid region (Fig. 5), errors between models and experimental results are slightly higher, however not greater than 1 K.

thumbnail Fig. 4

Predicted and experimental hydrate stability of the CO2 system in presence of aqueous sodium chloride solution in the VLwH region; Black lines: hydrate stability zone predicted using the CPA-EoS model using an aqueous mole fraction of 0.8; dotted blue lines: vapour–liquid locus of CO2 (no water) using the CPA-EoS model.

thumbnail Fig. 5

Predicted and experimental hydrate stability of the CO2 system in presence of aqueous sodium chloride solution in VLH and LLH regions; black lines: hydrate stability zone predicted using the CPA-EoS model using an aqueous mole fraction of 0.8; dotted blue lines: vapour–liquid locus of CO2 (no water) using the CPA-EoS model.

Table 4

Experimental CO2 hydrate dissociation conditions in the presence of sodium chloride solutions (AqFr: Aqueous Fraction).

Table 5 gives a summary of the measurements made in this work for undersaturated systems. As can be seen, in all cases increasing the pressure leads to an increase in the HET for the saturated system. All of the data for distilled water are summarized in Figures 68 showing the effect of CO2 saturation on HET for the tests with distilled water. The deviations between the measured HET and model predictions for distilled water are on average lower than 0.5 K and a maximum of 1 K for the lowest CO2 concentration.

thumbnail Fig. 6

Experimental and predicted hydrate equilibrium temperature at CO2 isopleth (xCO2 is in mol fraction) in distilled water (red lines: HLw predictions; black lines: VHLw/LHLw predictions).

thumbnail Fig. 7

Comparison of experimental and predicted solubility of CO2 in water at H–Lw and V–Lw equilibria.

thumbnail Fig. 8

Comparison of experimental and predicted solubility of CO2 in water at H–Lw equilibria at 5, 6, 10, 20, 30, and 40 MPa (lines from left to right).

Table 5

Hydrate equilibrium temperatures at HLw as function of CO2 concentration (xCO2) in distilled water and saline systems (CO2 concentration is on free-salt basis and mole fraction).

For salt solutions, the results are plotted in Figure 9 where larger deviations between experimental data and model predictions are observed. For some CO2 concentrations, the measured HET increases with increasing salinity, while the model predicts the opposite behavior. There is a need for further experimental measurements and model development for aqueous systems undersaturated with CO2.

thumbnail Fig. 9

Comparison of experimental and predicted solubility of CO2 in NaCl aqueous solutions at H–Lw equilibria (grey dashed lines predictions for distilled water at 10 MPa; green lines: predictions for 3 wt% NaCl at 10 (dashed) and 20 MPa; blue lines: predictions for 5.52 wt% NaCl at 10 (dashed), and 20 MPa; red lines: predictions for 6 wt% NaCl at 10 (dashed) and 20 MPa).

5 Conclusion

In this work, new experimental data have been reported for hydrates formed from carbon dioxide in the presence of aqueous solutions of sodium chloride over a wide range of temperatures, pressure, and concentration. The predictions of the developed model are compared against independent experimental data and the data generated in this work over a wide range of temperature, pressure, and NaCl concentration. A good agreement between predictions and experimental data is observed, demonstrating the reliability of the developed model. The hydrate stability pressure–temperature zone of dissolved CO2 in the presence of salt show that a decrease in the system pressure and/or an increase in salt concentration favors hydrate formation, as both factors reduce equilibrium gas solubility in the aqueous phase. This behavior is unlike that of the system including a gas phase, where a higher gas pressure corresponds to an increased gas solubility which favors hydrate formation.

To further validate the model, Hydrate Equilibrium (melting) Temperature (HET) of aqueous phases undersaturated with CO2, at different CO2 saturation levels and with different aqueous phase compositions were measured. The model is in relatively good agreement with the experimental results for distilled water, but larger deviations are observed for saline systems highlighting the need for further work to improve the model.

Acknowledgments

Some of this work was part of a Joint Industrial Project (JIP) conducted jointly at the Institute of GeoEnergy Engineering, Heriot-Watt University, and the CTP laboratory of MINES ParisTech between 2018 and 2021. The JIP was supported by Equinor, GALP Energia, Linde AG, Petrobras, Petronas, and TOTAL, which is gratefully acknowledged. HET measurements for undersaturated measurements were financed by Equinor.

References

Appendix

Temperature uncertainty

In this setup uncertainty of temperature measurement was calculated according to: U c ( T ) = u thermometer 2 + u calibration 2 + u repeatability 2 . $$ {U}_c(T)=\sqrt{{{u}_{\mathrm{thermometer}}}^2+{{u}_{\mathrm{calibration}}}^2+{{u}_{\mathrm{repeatability}}}^2}. $$(A.1)

The system, repeatability, and calibration uncertainties were found to be 0.1, 0.05, and 0.05 K, respectively. As a result, Uc(T) was found to be 0.12 K for all points.

Pressure uncertainty:

Similar to the temperature uncertainty, the uncertainty of pressure measurements was calculated according to the following equations: U c ( P ) = u system 2 + u calibration 2 + u repeatability 2 . $$ {U}_c(P)=\sqrt{{{u}_{\mathrm{system}}}^2+{{u}_{\mathrm{calibration}}}^2+{{u}_{\mathrm{repeatability}}}^2}. $$(A.2) Uc(P) was found to be 0.043 MPa with three digits accuracy for all the measurements.

HET for undersaturated systems

In this study, uncertainties of the measured hydrate equilibrium temperature were estimated based on the “Evaluation of Measurement Data – Guide to the Expression of Uncertainty in Measurement (GUM)” [37]. For HET measurements, considering the impact of uncertainties associated with calibration parameters, measured temperature, pressure, and sample preparation, the combined standard uncertainty of the measured density is calculated by: u ( T HET ) 2 = u ( T ) calibration 2 +   u ( T ) measurements 2 . $$ {u({T}_{\mathrm{HET}})}^2={u(T)}_{\mathrm{calibration}}^2+\enspace {u(T)}_{\mathrm{measurements}}^2. $$(A.3)

For the calibration uncertainties, the temperature probe was calibrated using a NAMAS-certified temperature probe. For the uncertainties of HET measurements, equation (A.4) was used to determine the standard uncertainties: u 2 ( T HET ) measurements = ( T HET p ) 2 u 2 ( p ) + ( T HET T ) 2 u 2 ( T ) + ( T x ) 2 u 2 ( x CO 2 ) . $$ {u}^2{({T}_{\mathrm{HET}})}_{\mathrm{measurements}}={\left(\frac{\mathrm{\partial }{T}_{{HET}}}{\mathrm{\partial }p}\right)}^2{u}^2(p)+{\left(\frac{\mathrm{\partial }{T}_{\mathrm{HET}}}{\mathrm{\partial }T}\right)}^2{u}^2(T)+{\left(\frac{\mathrm{\partial }T}{\mathrm{\partial }x}\right)}^2{u}^2\left({x}_{\mathrm{CO}2}\right). $$(A.4)

( T x ) $ \left(\frac{\mathrm{\partial }T}{\mathrm{\partial }x}\right)$ was calculated using the thermodynamic model. Considering the systematic uncertainties of measurement devices, the calibration procedures, and the reproducibility of the results, the standard uncertainties of the measured HET were found to be 0.74 K.

HET for saturated systems

For the uncertainties of HET measurements in presence of NaCl, the procedure described by Stringari et al. [38] was used, and the standard uncertainties of the measured HET were found to be 0.4 K.

All Tables

Table 1

Literature data for hydrate equilibria of CO2 in presence of NaCl solutions (VLwH: Vapour – Liquid (aqueous) – Hydrate / LLwH: Liquid – Liquid (aqueous) – Hydrate).

Table 2

Literature data of CO2 solubility in water and brines at hydrate–liquid water (HLw) equilibrium.

Table 3

Composition of the chemical used in this work.

Table 4

Experimental CO2 hydrate dissociation conditions in the presence of sodium chloride solutions (AqFr: Aqueous Fraction).

Table 5

Hydrate equilibrium temperatures at HLw as function of CO2 concentration (xCO2) in distilled water and saline systems (CO2 concentration is on free-salt basis and mole fraction).

All Figures

thumbnail Fig. 1

Schematic of the setup used to measure the hydrate dissociation condition

In the text
thumbnail Fig. 2

Plot of equilibrium P/T data recorded during test xCO2 = 2.11 × 10−2 mole fraction in distilled water.

In the text
thumbnail Fig. 3

Predicted and experimental CO2 solubility in aqueous sodium chloride. Black lines: predictions using the CPA-EoS model. : 5.52 wt% NaCl [24]; : 14.92 wt% NaCl [24]; : 22.61 wt% NaCl [24]; : 10 wt% NaCl (Heriot–Watt unpublished data); : 20 wt% NaCl (Heriot–Watt unpublished data).

In the text
thumbnail Fig. 4

Predicted and experimental hydrate stability of the CO2 system in presence of aqueous sodium chloride solution in the VLwH region; Black lines: hydrate stability zone predicted using the CPA-EoS model using an aqueous mole fraction of 0.8; dotted blue lines: vapour–liquid locus of CO2 (no water) using the CPA-EoS model.

In the text
thumbnail Fig. 5

Predicted and experimental hydrate stability of the CO2 system in presence of aqueous sodium chloride solution in VLH and LLH regions; black lines: hydrate stability zone predicted using the CPA-EoS model using an aqueous mole fraction of 0.8; dotted blue lines: vapour–liquid locus of CO2 (no water) using the CPA-EoS model.

In the text
thumbnail Fig. 6

Experimental and predicted hydrate equilibrium temperature at CO2 isopleth (xCO2 is in mol fraction) in distilled water (red lines: HLw predictions; black lines: VHLw/LHLw predictions).

In the text
thumbnail Fig. 7

Comparison of experimental and predicted solubility of CO2 in water at H–Lw and V–Lw equilibria.

In the text
thumbnail Fig. 8

Comparison of experimental and predicted solubility of CO2 in water at H–Lw equilibria at 5, 6, 10, 20, 30, and 40 MPa (lines from left to right).

In the text
thumbnail Fig. 9

Comparison of experimental and predicted solubility of CO2 in NaCl aqueous solutions at H–Lw equilibria (grey dashed lines predictions for distilled water at 10 MPa; green lines: predictions for 3 wt% NaCl at 10 (dashed) and 20 MPa; blue lines: predictions for 5.52 wt% NaCl at 10 (dashed), and 20 MPa; red lines: predictions for 6 wt% NaCl at 10 (dashed) and 20 MPa).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.