Experimental and theoretical investigation of the mechanisms of drying during CO2 injection into saline reservoirs

Scientific Reports volume 13, Article number: 9155 (2023) Cite this article

1 Altmetric

Metrics details

A viable CO2 storage resource must have sufficient storage capacity, reliable containment efficiency and adequate well injectivity. Deep saline formations stand out in terms of storage capacity and containment efficiency. However, formation brine dry-out and salt precipitation in the near well region could impair CO2 injectivity in deep saline reservoirs, thus reducing their potential for CO2 storage. Core-flood experiments and analytical modelling were used to investigate various mechanisms of external and internal salt precipitation. Particularly, the impact of the extension of the dry-out region on CO2 injectivity was investigated. It was found that, for high permeability rocks, injection of CO2 at relatively low injection rates could result in salt cake deposition at the injection inlet especially under high salinity conditions. It was also found that extension of the dry-out region does not have significant impact on CO2 injectivity. Although the magnitude of CO2 injectivity impairment increased more than two-fold when initial brine salinity was doubled, real-time changes in CO2 injectivity during the drying process was found to be independent of initial brine salinity. We have shown that the bundle-of-tubes model could provide useful insight into the process of brine vaporization and salt deposition in the dry-out region during CO2 injection. This work provides vital understanding of the effect of salt precipitation on CO2 injectivity.

The prerequisites for successful Carbon Capture, Utilization and Storage (CCUS) are robust containment efficiency, adequate storage volume and sufficient well injectivity, to inject large quantities of CO2 at practical flow rates1. Deep saline formations are suitable storage resources for CCUS based on their storage capacity and containment2,3,4,5,6,7. However, salt precipitation due to brine vaporization, especially in the vicinity of the wellbore, during CO2 injection could affect CO2 injectivity in deep saline formations8,9,10,11,12,13. The underlying factors of CO2 injectivity impairment induced by salt precipitation have been extensively studied and identified to include the concentration of salt in the brine, the petrophysical and petrographic properties of the rock, the drying rate, the extent of the dry-out zone, solid salt saturation in the pore spaces after drying, distribution of precipitated salts within the pores and the petrophysical properties of the reservoir rock14,15.

Salt precipitation or scaling has been a major formation damage challenge in oilfield operations since the inception of the industry. In field operations involving injection, storage and production of natural gas; various levels of injectivity impairment related directly and indirectly to salt precipitation have been encountered and reported16,17,18,19. Permeability impairment ranging between 13 and 83% and porosity reduction around 2–15% have been reported from laboratory experiments9,11,12,20,21,22,23,24. Rigorous theoretical simulations have also confirmed the reported experimental and field findings23,25,26,27,28,29,30. Cui et al., (2023) have compiled a more recent update of injectivity impairment induced by salt precipitation that have been reported by various researchers through experimental and modelling studies. Generally, porosity impairment has been lower than permeability changes as the deposition of salt in the flow pathways have more impact on permeability than porosity.

During injection of CO2 into brine-filled rock, the injected gas initially displaces mobile brine out of the rock. During this immiscible displacement stage, mass transfer between the gas and the displaced aqueous phase is minimal. After the mobile brine has been displaced, continuous injection of CO2, especially under typical field injection conditions, results in vaporization of brine, drying and salt precipitation. Generally, the dry-out zone extends into the formation with injection time after the onset of the drying process. Some experimental and numerical studies have in part examined mechanisms underlying the development of the dry-out zone qualitatively and quantitatively9,25,27,31,32,33,34,35. However, to the best of knowledge of the authors, there has not been an experimental or modelling studies that has attempted to monitor the extension of the dry-out zone in real-time and examine its impact on injectivity impairment.

Miri and Hellevang14 have identified the main underlying mechanisms that govern the drying rate and extension of the dry-out zone during salt precipitation. These factors include: (1) immiscible two-phase displacement of the resident brine by injected CO2, (2) vaporization of brine into the flowing CO2 stream, (3) capillary-driven back-flow of brine towards the injection inlet, (4) diffusion of dissolved salt in the aqueous phase, (5) gravity override of injected CO2, and (6) salt self-enhancing. These factors have also been confirmed under practical field injection conditions36. It has also been reported that very low brine evaporation rate in the drying front, may result in homogeneous distribution of precipitated salt throughout the porous medium9,10,14,37. For high vaporization rates, there is no sufficient time for the salt concentration gradient to diffuse away from the drying front, resulting in nonhomogeneous accumulation of salt10,38. The position where salt accumulation is maximum is still largely debatable. Numerical experiments conducted by Roels et al.39 suggests that precipitated salt may accumulate far from the wellbore. However, other research works6,7,9,17,22,40 report that precipitated salt accumulates near the wellbore where the fluxes and brine vaporization are the highest. However, Berntsen et al.41 identified three different drying regimes in different wellbore regions when they investigated drying and salt clogging under near-realistic radial CO2 flow conditions. This suggests that the distribution of precipitated salt is not uniform across the drying region although the exact relative distribution of precipitated salt and the governing factors are yet to be thoroughly investigated. A detailed review and assessment of CO2 injectivity impairment mechanisms has been presented by Hajiabadi et al.42.

Several analytical and numerical models have been developed to study the physics behind salt precipitation at the core and field scales27,28,43,44,45. More recently, machine learning-based modelling has also been adopted to study the mechanisms of CO2 injectivity46,47,48. The microstructure of natural porous media are very complex with tortuousness and often noncircular pore channels49,50,51. Three pore-scale models have been widely used in attempt to reconstruct representative analogs of porous structure to study fluid flow in porous media; the network models, the sphere-pack model, and the bundle-of-tubes model52. The simplest pore-scale model is the bundle-of-tubes model which is normally derived from the Hagen–Poiseuille equation53. Early traditional bundle-of-tubes models represent porous media as an assemblage of independent, uniform circular capillary tubes54,55. Bundle-of-tubes models may be interacting or non-interacting, depending on fluid communication between the individual tubes56.

In the present work, we have conducted core-flood experiments using sandstone core plugs to investigate the mechanisms of internal and external salt deposition during CO2 injection into deep saline reservoirs. Particularly, we have investigated the development of the dry-out zone in near-real-time and attempted to monitor and quantify the impact of the extension of the drying region on CO2 injectivity. We then derived a tractable bundle-of-tubes model describing the experimental observations, in attempts to model the relationship between CO2 injectivity impairment and the extent of the dry-out zone, i.e.to establish a quantitative relationship between the development of the dry-out zone and CO2 injectivity. The results have been discussed in light of established literature.

Outcrop Berea sandstone core plugs considered to be homogeneous with brine permeability in the range of 60–100 mD and porosity between 18 and 20% were used as the main reservoir rock for the study. Each core sample has a length of 20 cm and 3.81 cm diameter. These long core plugs were selected purposely to increase the residence time of CO2 during the drying process.

Liquefied CO2 with percentage purity of about 99.7%, was used as the non-wetting phase. During liquid CO2 flooding, the fluid was injected at 80 bar and 25 °C. For supercritical CO2 flooding, the injection conditions were set to 80 bar and 45 °C.

NaCl brines, salinity of 75 g/l and 150 g/l, were selected to represent low salinity (LS) and high salinity (HS) formation brines, respectively. The HS brine is expected to represent deep saline formation brine, while the LS brine was selected to test the sensitivity of brine salinity.

Figure 1 depicts schematics of the core-flooding apparatus used in the experimental work. Initially, the core was loaded into the core-holder. The Quizix pump is used to deliver brine through the connected piston cell for temporary aging to attain a set uniform temperature before it arrives at the core inlet. The ISCO CO2 pump receives liquid CO2 from the gas bottle through a pressure regulator. Either liquid or supercritical CO2 can be injected into the core depending on the set conditions. The injected fluid passes the piston cell which holds the fluid until a preset temperature is attained in the oven. A differential pressure gauge and a pressure transducer are used to monitor the pressure drop across the core and record the pore pressure. A backpressure of 80 bar was set at the outlet, during CO2 injection and the effluent fluid was safely collected in a piston cell for analysis and safe disposal.

The experimental set up used for the CO2 core-flooding experiments.

The core sample was first cleaned and dried at 65 °C for about 24 h. The core was then wrapped in shrinking Teflon sleeve before it was inserted into a rubber sleeve to prevent CO2 leakage. The carefully wrapped core sample was then mounted into the core holder. A confining pressure of about 150 bar was applied to the core in the core holder. The experimental procedure consists of the following:

Initial liquid CO2 permeability (\({K}_{i})\) of the core sample was measured.

The core was vacuum-saturated with a specific salinity brine (HS or LS).

The saturated core sample was flooded with specific pore volumes (PV) of supercritical CO2 to vaporize brine and dry the core starting at the injection inlet.

The core was taken out, inspected, and the extent of the dry-out zone measured.

Final liquid CO2 permeability (\({K}_{f}\)) of the core was the measured.

In Step 2, the vacuum-saturated core was prepared and inserted into the core-holder followed by brine injection at 1 ml/min until complete saturation was attained. In Step 3, supercritical CO2 was injected into the saturated core at 5 ml/min to vaporize brine and possibly precipitate salt into the pores. After each round of supercritical CO2 injection, the core was inspected, and the extent of the dry-out zone was measured by scanning. In Step 1 and Step 5, liquid CO2 was injected at 2 ml/min to measure permeability before and after brine vaporization, and salt precipitation.

Pressure drop profiles and pore-pressure were recorded during CO2 flooding to qualitatively study the drying process and the impact of the drying on injectivity. Pressure drop profiles might not provide thorough information of pore level events, but they give useful insight into real time changes in flow properties in the core during the drying process.

For constant injection rates (\({q}_{i}={q}_{f}\)) under linear flow conditions in a homogenous rock, we may define the relative injectivity change index, \(\beta \) from Darcy law as:

The terms \(\left(\frac{{\Delta p}_{i}}{{\Delta p}_{f}}\right)\) and \(\left(\frac{{K}_{f}}{{K}_{i}}\right)\) are evaluated and compared for consistency. The \(\beta \) provides an indirect estimation of CO2 injectivity induced by the drying process. Salt accumulation in the pore constrictions will increase the pressure drop, thus decreasing the rock permeability and CO2 injectivity.

The main objective of the modelling work is to develop a tractable physical analytical model with just enough capability to complement the experimental studies in terms of estimating injectivity impairment induced after salt precipitation without sacrificing the effectiveness of the model to capture the major underlying mechanisms. The bundle of tubes model was selected because of its ability to capture the main physical processes of transport in porous media, its computational efficiency and flexibility for modelling different pore-scale events. The main assumptions of the model include:

A fully homogeneous reservoir rock

No chemical interaction between the injection fluid and the contents of the reservoir rock

Single phase flow of injected CO2 in the reservoir during injectivity tests

Uniform deposition of precipitated salt in the dry-out region of the rock

Although the model is limited by these assumptions, it captures the basic mechanisms and therefore capable of providing acceptable estimates of the quantitative impact of salt deposition on CO2 storage in deep saline reservoirs at the core-scale.

Figure 2 shows a Berea core that was initially saturated with LS brine and flooded with about 160 PV of supercritical CO2 at an injection flow rate of 5 ml/min. It can be observed that the part of the core close to the injection inlet dried out completely after brine vaporization while the remaining section close to the outlet remained wet with immobile brine.

A Berea sandstone core after a period of drying. The core was initially saturated with LS brine after which it was flooded with about 160 PV of supercritical CO2 at injection flow rate of 5 ml/min.

By inspection, the core can be sectioned into a dry-out region and a wet region after brine vaporization and drying (Fig. 2). Salt precipitation is expected in the dry-out region after brine vaporization. The cylindrical core, radius \(R\) and length \(L\) was reconstructed with a bundle of parallel cylindrical capillary tubes with varying radii (\({r}_{1},{r}_{2},{r}_{3},\dots {r}_{N}\)) interspersed between a non-porous mass (shaded regions) which represent the rock matrix (Fig. 3). The model is sectioned into a dry-out region (\({L}_{1}\)) where salt has been precipitated into the tubes and the wet region (\({L}_{2}\)) where the pores contain immobile brine.

A Schematic of the bundle-of-tubes model. The core of length L and radius R is represented by a bundle of parallel cylindrical capillary tubes with varying radii interspersed between non-porous mass.

In the dry-out region, the precipitated solid salt is assumed to accumulate on the walls of the pore constriction, represented by the capillary tubes. As a result, the accumulated salt in the dry-out region would reduce the flow area of the tube \(r\) by \(\Delta r\), for a total of \(N\) capillary tubes in the porous medium.

We define a dimensionless parameter \({l}_{d},\) to monitor the extension of the dry-out region from the inlet to the outlet of the core of total length \(L\) given by:

Using the Hagen − Poiseuille equation, the fluid injectivity through the bundle of capillary tubes can be expressed as:

In Eq. (3), \(\mu \) is the viscosity of the injected fluid, \(Q\) is the injection flow rate and \(\Delta p\) is the total pressure drop across the bundle of tubes. The relative injectivity change index \(\beta \) can then be derived as:

In Eq. (4), \({I}_{i}\) and \({I}_{f}\) are the fluid injectivity before and after drying and salt precipitation, respectively. If the pore size distribution of the rock is known, Eq. (4) can be used to estimate the impact of precipitated salt on injectivity if \({l}_{d}\) and \({\Delta r}_{i}\) for each \({r}_{i}\) of the total \(N\) capillary tubes are known. These are the three unknown parameters in the equation that must be determined.

The pore volume of the core is approximately equal to the sum of the internal volume of all the capillary tubes. Using the definition of porosity,\(\phi ,\) the total number of capillary tubes that is representative of the pore volume of the core can be related to the porosity, the size of the core \(R\) and the average pore radius, \(\overline{{r }_{i}}\) by:

Integrating \({r}_{i}^{2}\) from zero to the maximum salt thickness \({\Delta r}_{max}\) gives:

Similarly, integrating \({r}_{i}\) over the same salt thickness yields:

Combining Eqs. (6) and (7) gives a relationship between \(\overline{{r }_{i}^{2}}\) and \(\overline{{r }_{i}}\):

where \({\overline{r} }_{i}^{2}\) is the square of the average of the capillary tube radii. Substituting Eq. (8) into Eq. (5), the total number of capillary tubes that will give a specific porosity can be estimated:

From Eq. (9), the total number of capillary tubes that represents the pore volume of the core quantitatively can be estimated, given the porosity \(\phi \) and average pore size, \(\overline{{r }_{i}}\), given the pore size distribution of the rock.

The solid salt saturation \({S}_{si}\) deposited in a capillary tube is defined by:

where \({V}_{si}\) is the volume of solid salt in the tube,\({V}_{pi}\) is the volume available in the tube for salt deposition and \({r}_{e}=({r}_{i}-{\Delta r}_{i}\)). Substituting for \({r}_{e}=({r}_{i}-{\Delta r}_{i}\)) gives:

In Eq. (11), assuming \({r}_{i}\gg {\Delta r}_{i}\) that \(\frac{\Delta {r}_{i}}{{r}_{i}^{2}}=0\) we can derive \({\Delta r}_{i}\) in terms of \({S}_{si}\):

The total mass of solid salt deposited in all the \(N\) tubes can be estimated as:

where \({\rho }_{s}\) is the density of solid salt. But the total mass of salt deposited in the core can also be expressed as:

where \({V}_{st}\) is the total volume of salt deposited in the core. Combining Eqs. (13) and (14) gives:

Substituting Eqs. (9) into (15) gives:

Since \({S}_{s}\) and \({l}_{d}\) are constant at any point in time, the uncertainty in \(\overline{{\Delta r }_{i}}\) in Eq. (16) is mainly associated with uncertainty in \(\overline{{r }_{i}}\). Therefore, Eq. (16) can be used to estimate the average thickness of the deposited solid salt on the pore walls in terms of the cumulative solid salt saturation, \({S}_{s}\) by:

If \({S}_{s}\) and \({l}_{d}\) are known, we can compute \({\Delta r}_{i}\) for each \({r}_{i}\) after salt precipitation. Using mass balance, Pruess27 derived an equation to estimate the solid salt saturation as:

In Eq. (18), \({\overline{S} }_{g,d}\) is the average gas saturation behind the dry-out front, \({\rho }_{aq}\) is the density of brine, \({X}_{s}\) is the mass fraction of salt in the brine and \({\rho }_{s}\) is the density of solid salt. Since \({\overline{S} }_{g,d}\) is also a measure of the extension of the dry-out region similar to \({l}_{d}\) in the current model, a correlation was derived for the solid salt saturation, by fitting experimental data:

First, we measured \({\rho }_{aq}\) and \({X}_{s}\) for a 100 g/l NaCl brine. Then, a Berea sandstone core-sample with known initial permeability, was saturated with the brine and about 300 PV of supercritical CO2 was flooded through the core at a rate of 5 ml/min. After every 50 PV of CO2 injections, the core was inspected by scanning to determine \({l}_{d}\) after which permeability was measured. We then calculated \(\beta \) for each \({l}_{d}\) and plotted the data. An initial \({S}_{s}\) correlation similar to Eq. (18) by replacing \({\overline{S} }_{g,d}\) with \({l}_{d}\) without the matching parameters was used to compute \(\beta \). The saturation correlation \({S}_{s}\) was then optimized to fit the initial experimental data. The optimum correlation in Eq. (19) was then used throughout the various experiments. A brine concentration outside the range used in the main experimental set was used to calibrate the correlation and to ensure repeatability.

Equation (19) is used to estimate \({S}_{s}\) in the dry-out zone. Once \({S}_{s}\) is known, \({\Delta r}_{i}\) are computed from Eq. (17) and then \(\beta \) for each \({l}_{d}\) calculated from Eq. (4).

The mechanisms of salt precipitation can be grouped into two: salt cake formation at the injection inlet and drying effects. Under some favorable conditions, salt cake may form on the surface of the core inlet during early stages of brine vaporization before the onset of drying.

To investigate mechanisms of salt cake formation, a clean Bentheimer core was vacuum-saturated with 120 g/L NaCl brine and flooded with about 50 PV of dry supercritical CO2 at a rate of 1 mL/min. The Bentheimer core was selected because of its relatively high permeability. Pressure drop profiles were recorded during the CO2 flooding. From Fig. 4A, it can be observed that no salt was formed at the core outlet. Figure 4B shows massive salt cake deposition at the core inlet. At the onset of injection, when the core is fully saturated with brine, the injected supercritical CO2 leaves brine behind the inlet due to poor sweep of the brine by the injected fluid. Salinity of the brine increases gradually due to mass transfer of water from the brine into the dry supercritical CO2. If the initial brine salinity is high enough, the brine left behind at the inlet could reach supersaturation and precipitate salt onto the inlet surface. The precipitated salt creates a saturation gradient that draws more brine into the inlet region through capillary backflow, leading to precipitation of more salts at the inlet. Thus, the salt cake formation could increase with brine salinity and inefficient brine displacement at the injection inlet.

Photographs of Bentheirmer core after CO2 was injected at 1 mL/min into the core initially saturated with 120 g/L NaCl brine. (A) No salt cake observed at the core outlet. (B) Massive salt cake found at the injection inlet.

To investigate the impact of sweep on salt cake development, supercritical CO2 injection flow rate was increased from 1 to 5 mL/min. Salt deposition at the injection inlet decreased when the injection flow rate was increased from 1 to 5 mL/min (Fig. 5). The sweep at the injection inlet improved with increasing injection flow rate, leaving less brine behind the injection inlet for salt precipitation.

Photographs showing salt cake development at the core inlet when supercritical CO2 injection rate was increased from (A) 1 mL/min to (B) 5 mL/min. Increase in CO2 injection rate decreased the amount of deposited salt.

The initial brine salinity was then reduced to 75 g/L, keeping the CO2 injection flow rate constant at 5 mL/min. The amount of salt cake precipitated at the injection inlet further decreased significantly when brine salinity was decreased (Fig. 6). At constant rate of vaporization, lowering the saturating brine salinity delays supersaturation, allowing a significant portion of the brine left behind at the injection inlet to be swept into the core which in tend reduces the amount of salt cake formed at the inlet.

Photographs showing salt cake development at the core inlet when brine salinity was decreased from (A) 120 g/L to (B) 75 g/L. Decrease in brine salinity further decreased salt deposition.

From these experimental studies, it can be observed that salt cake deposition at the injection inlet during CO2 injection into saline porous media may depend on saturating brine salinity and the brine sweep at the injection inlet.

Pore size distribution is required to model the porous medium with the bundle-of-tubes model derived in section "Modelling work". Based on a pore size distribution analysis from mercury injection on a Berea sandstone, Shi et al.57 have found an average pore radius of about 6.7 μm. Dullien and Dhawan58 reported pore constriction sizes between 0.5 and 5.0 μm and pore chamber sizes ranging from 5.0 to 50 μm in Berea sandstone. From these data, we calibrated our model to an average pore size of 6 μm using a lognormal distribution of tube radii, as shown in Fig. 7. In Fig. 7, it can be observed that tubes with radii greater than 20 μm consist of less than 5% of the total number of capillary tubes. The minimum tube radius was set to 0.38 μm.

Distribution of tube radii representing pore size distribution of Berea sandstone core. The tube radii were drawn from a lognormal distribution with average tube radius of about 6 μm which is the average pore throat size of Berea sandstone.

Natural reservoir rocks have open pore bodies connected by two or more pore throats depending on the coordination number of the type of rock. The average coordination number of Berea sandstones is between 4 and 859,60. The pore size distribution in Fig. 7 compares favorably with this range of coordination number.

A Berea sandstone core with known initial permeability was saturated with LS brine prior to the drying experiment. The core was then flooded with supercritical CO2 at 5 ml/min. After every 100 PV of CO2 injection, the core was taken out and inspected, and the advancing dry-out region was measured to estimate \({l}_{d}\). CO2 injectivity change caused by salt precipitation was also measured and \(\beta \) was computed. Figure 8 shows the impact of the advancing dry-out region, \({l}_{d}\) on CO2 injectivity impairment, \(\beta \) for the experimental and modelling studies.

Effect of the advancing dry-out region on CO2 injectivity impairment. The core was initially saturated with LS brine. CO2 Injectivity change, and extent of the dry-out region were measured after every 100 PV of CO2 injection.

Figure 8 shows the simulation data agrees favorably with the measured data. It can be observed that CO2 injectivity impairment was highest at the onset of drying. Injectivity impairment decreased to a minimum for \({l}_{d}\) of about 0.45 and then increased slightly towards the core effluent end. At the commencement of drying, injectivity impairment could be induced by salt precipitation and relative permeability effects. Brine vaporization rate is highest close to the inlet region because of the high capillary driven back-fluxes. When the brine attains supersaturation, salt could be precipitated into the pores which in turn could reduce injectivity. Also, at the onset of drying, a large portion of the rock still contains brine. Capillary backflow may also drive the aqueous phase towards the injection inlet. The presence of mobile brine in the pores, could reduce the space available for CO2 flow, thus decreasing the relative permeability of CO2, which could also have an impact on injectivity around the inlet region.

As the drying front advances into the core, brine vaporization and salt precipitation are expected to decrease as a result of the reduced fluxes. In addition, the mobile brine saturation will decrease as more brine is swept out of the core. This will increase the relative permeability to CO2. When the core is almost fully dried, brine vaporization and salt precipitation at the outlet end of the core are almost negligible. Also, the absence of mobile brine in the rock further improves the relative permeability of the rock to CO2. However, the capillary driven back fluxes and changes in the distribution of deposited salt as a result of the continuous injection of large pore volumes of CO2 after salt deposition was probably responsible for the slight increase in injectivity impairment towards the end of the rock. Salt redistribution within the pores can cause heterogeneous salt deposition which tends to impair CO2 injectivity as also observed by other researchers14,38,61.

To investigate the effect of brine salinity, the experiment was repeated with HS brine. The measured and simulated results for the LS and HS brine cases are presented in Fig. 9.

Effect of brine salinity on CO2 injectivity in relation to advancement of the dry-out region. The magnitude of CO2 injectivity impairment is seen to increase more than two-fold when brine salinity was doubled from 75 to 150 g/l.

Generally, CO2 injectivity impairment increased more than two-fold when brine salinity was doubled from 75 to 150 g/l (Fig. 9). This result compares well with the findings of Jeddizahed and Rostami62 who reported increased permeability impairment from 21 to 50%, when brine salinity was increased from 50 to 100 g/l, under similar experimental conditions. Other researchers8,25,27 have also reported direct increase in CO2 injectivity impairment proportional to increase in brine salinity.

Figure 9 also suggests that the same mechanisms are responsible for CO2 injectivity impairment as the dry-out region extents from the injection inlet regardless of the initial saturating brine salinity. Increase in brine salinity from LS to HS, increased the magnitude of CO2 injectivity impairment but had no direct impact on the development of CO2 injectivity impairment during the drying process. Also, the simulation results compare favorably with the measured data for the increased brine salinity.

Generally, solid salt saturation increases with increase in brine salinity10,20,26,44. However, it can be observed that increase in brine salinity had negligible effect on the mechanisms underlying the changes in CO2 injectivity as the dry-out zone extends into the rock. The two main mechanisms, rate of salt precipitation during the drying process and relative permeability effects are for the most part dependent on CO2 injection rate. Increase in brine salinity only increases the magnitude of salt precipitated but the rate of precipitation depends primarily on the brine vaporization rate. Also, the relative permeability effects are influenced by the amount of mobile brine present in the rock in real time.

In terms of storage space, deep saline formations are the most viable candidates for CCUS. However, salt precipitation during CO2 injection into deep saline reservoirs could impair CO2 injectivity and reduce storage potential. Insight into the drying process, the changes in the dry-out region and its impact on CO2 injectivity could improve understanding of the mechanisms of drying and salt precipitation in saline reservoirs. Core-flood experiments were conducted to investigate the mechanisms of drying during CO2 injection into saline reservoirs and the impact of the size of the advancing dry-out region on CO2 injectivity. The experimental results were used to calibrate a bundle-of-tubes model to track the size of the dry-out region and its impact on CO2 injectivity impairment induced by solid salt precipitation. Some highlights of our findings include the following:

Salt could be precipitated externally on the surface of the injection inlet in the form of filter a cake. Brine salinity and sweep efficiency of brine around the injection inlet have been identified as the main underlying parameters of external salt precipitation.

For drying and internal salt precipitation, CO2 injectivity impairment was highest at the onset of the drying process due to high brine vaporization rate and low relative permeability of CO2. Injectivity improved slightly as the drying front advanced into the middle of the core but decreased gradually towards the end of the drying process. Overall, increase in the dry-out region did not impose significant effect on the magnitude of CO2 injectivity.

Doubling the initial brine salinity from 75 to 150 g/l reduced CO2 injectivity more than two-fold, but the successive change in injectivity due to the extension of the dry-out region was independent of initial brine salinity.

Although the analytical model was largely basic, the accompanying experimental findings are very important for understanding the mechanisms of drying and salt precipitation in saline reservoirs. Insight from the current study could provide a solid basis for understanding the relationship between CO2 injectivity induced by salt deposition and the advancement of a dry-out region.

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

IEA. Technology roadmap-Carbon capture and Storage. Technol. Roadmap [Internet]. 59 (2013).

Tang, Y. et al. Experiment on CO2-brine-rock interaction during CO2 injection and storage in gas reservoirs with aquifer. Chem. Eng. J. [Internet]. 413, 127567 (2021).

Article CAS Google Scholar

Kumar, S. et al. A comprehensive review of value-added CO2 sequestration in subsurface saline aquifers. J. Nat. Gas Sci. Eng. 81, 103437 (2020).

Article CAS Google Scholar

Ringrose, P. S. et al. Storage of carbon dioxide in saline aquifers: Physicochemical processes, key constraints, and scale-up potential. Annu. Rev. Chem. Biomol. Eng. 12, 471–494 (2021).

Article CAS PubMed Google Scholar

Arora, V. et al. Separation and sequestration of CO2 in geological formations. Mater. Sci. Energy Technol. 2, 647–656 (2019).

Google Scholar

Yusof, M. A. M. et al. Experimental study of CO2 injectivity impairment in sandstone due to salt precipitation and fines migration. J. Pet Explor. Prod. Technol. 1, 1 (2022).

Google Scholar

Md Yusof, M. A. et al. Influence of brine-rock parameters on rock physical changes during CO2 sequestration in saline aquifer. Arab. J. Sci. Eng. 1, 1–15 (2021).

Google Scholar

Bacci, G., Durucan, S. & Korre, A. Experimental and numerical study of the effects of halite scaling on injectivity and seal performance during CO2 injection in saline aquifers. Energy Proc. 37, 3275–3282 (2013).

Article CAS Google Scholar

Peysson, Y., André, L. & Azaroual, M. Well injectivity during CO2 storage operations in deep saline aquifers-Part 1: Experimental investigation of drying effects, salt precipitation and capillary forces. Int. J. Greenh. Gas Control. 22, 291–300 (2014).

Article CAS Google Scholar

Miri, R. et al. New insights on the physics of salt precipitation during injection of CO2 into saline aquifers. Int. J. Greenh. Gas Control. 43, 10–21 (2015).

Article CAS Google Scholar

Muller, N. et al. CO2 injection impairment due to halite precipitation. Energy Proc. 1, 3507–3514 (2009).

Article CAS Google Scholar

Edem, D., Abba, M., & Nourian, A. et al. Experimental investigation of the extent of the impact of Halite precipitation on CO2 injection in deep saline aquifers. SPE Eur. OnePetro (2020).

Cui, G. et al. The influence of complicated fluid-rock interactions on the geothermal exploitation in the CO2 plume geothermal system. Appl. Energy. 227, 49–63 (2018).

Article CAS Google Scholar

Miri, R. & Hellevang, H. Salt precipitation during CO2 storage—A review. Int. J. Greenh. Gas Control. 51, 136–147 (2016).

Article CAS Google Scholar

Cui, G. et al. A review of salt precipitation during CO2 injection into saline aquifers and its potential impact on carbon sequestration projects in China. Fuel 334, 126615 (2023).

Article CAS Google Scholar

Jasinski, R., Sablerolle, W. & Amory, M. ETAP: scale prediction and control for the heron cluster (Society of Petroleum Engineers, 1997).

Google Scholar

Kleinitz, W., Koehler, M., & Dietzsch, G., et al. The precipitation of salt in gas producing wells. In: SPE, editor. SPE Eur Form damage Conf. The Hague: SPE. p. 1–7 (2001).

Golghanddashti, H. et al. Experimental investigation of water vaporization and its induced formation damage associated with underground gas storage. J. Porous Media. 1, 16 (2013).

Google Scholar

Place Jr., M. C. an unusual case of salt plugging in a high-pressure sour gas well. In: SPE, editor. 59th Annu Tech Conf Exhib. Houston, Texas: SPE; p. 13 (1984).

Tang, Y. et al. Experimental study of formation damage caused by complete water vaporization and salt precipitation in sandstone reservoirs. Transp. Porous Media. 107, 205–218 (2015).

Article CAS Google Scholar

Kim, M., Sell, A. & Sinton, D. Aquifer-on-a-Chip: Understanding pore-scale salt precipitation dynamics during CO 2 sequestration. Lab. Chip. 13, 2508–2518 (2013).

Article CAS PubMed Google Scholar

Bacci, G., Korre, A. & Durucan, S. Experimental investigation into salt precipitation during CO2 injection in saline aquifers. Energy Proc. 4, 4450–4456 (2011).

Article CAS Google Scholar

Ott, H., Snippe, J. & de Kloe, K. Salt precipitation due to supercritical gas injection: II Capillary transport in multi porosity rocks. Int. J. Greenh. Gas Control. 105, 103233 (2021).

Article CAS Google Scholar

Sokama-Neuyam, Y. A. & Ursin, J. R. The coupled effect of salt precipitation and fines mobilization on CO2 injectivity in sandstone. Greenh. Gases Sci. Technol. 8, 1066–1078 (2018).

Article CAS Google Scholar

André, L., Peysson, Y. & Azaroual, M. Well injectivity during CO 2 storage operations in deep saline aquifers—Part 2: Numerical simulations of drying, salt deposit mechanisms and role of capillary forces. Int. J. Greenh. Gas Control. 22, 301–312 (2014).

Article Google Scholar

Guyant, E. et al. Salt precipitation and CO2/brine flow distribution under different injection well completions. Int. J. Greenh. Gas Control. 37, 299–310 (2015).

Article CAS Google Scholar

Pruess, K. Formation dry-out from co2 injection into saline aquifers: 2 analytical model for salt precipitation. Water Resour. Res. 45, 1–6 (2009).

Google Scholar

Grimm Lima, M. et al. Permeability impairment and salt precipitation patterns during CO2 injection into single natural brine-filled fractures. Water Resour. Res. 56, 7213 (2020).

Article Google Scholar

Hajirezaie, S. et al. Numerical simulation of mineral precipitation in hydrocarbon reservoirs and wellbores. Fuel 238, 462–472 (2019).

Article CAS Google Scholar

Parvin, S. et al. Continuum scale modelling of salt precipitation in the context of CO2 storage in saline aquifers with MRST compositional. Int. J. Greenh. Gas Control. 99, 103075 (2020).

Article CAS Google Scholar

Hurter, S., Berge, J. & Labregere, D. Simulations for CO2 injection projects with compositional simulator. Offshore Eur. 1, 4–7 (2007).

Google Scholar

Abdelaal, M. & Zeidouni, M. Pressure falloff testing to characterize CO2 plume and dry-out zone during CO2 injection in saline aquifers. Int. J. Greenh. Gas Control. 103, 103160 (2020).

Article CAS Google Scholar

Norouzi, A. M. et al. Analytical solution for predicting salt precipitation during CO2 injection into saline aquifers in presence of capillary pressure. Water Resour. Res. 58, 32612 (2022).

Article ADS Google Scholar

Grimm Lima, M. M., Schädle, P., & Vogler, D., et al. A numerical model for formation dry-out during CO2 injection in fractured reservoirs using the MOOSE framework: implications for CO2-based geothermal energy extraction. World Geotherm Congr (WGC 2020+ 1). ETH Zurich, Institute of Geophysics; p. 15024 (2020).

Huijskes, T. & Schiferli, W. Implications of near-wellbore phenomena during CO2 injection in depleted gas fields. EAGE GET 2022(2022), 1–5 (2022).

Google Scholar

Talman, S. et al. Salt precipitation at an active CO2 injection site. Gas Inject. Geol. Form Relat. Top. 1, 183–199 (2020).

Google Scholar

Ott, H., Roels, S. M. & Kloe, K. Salt precipitation due to supercritical gas injection: I Capillary-driven flow in unimodal sandstone. Int. J. Greenh. Gas Control 43, 247–255 (2015).

Article CAS Google Scholar

Peysson, Y., Fleury, M. & Blázquez-Pascual, V. Drying rate measurements in convection-and diffusion-driven conditions on a Shaly sandstone using nuclear magnetic resonance. Transp. Porous Med. 90, 1001–1016 (2011).

Article CAS Google Scholar

Roels, S. M., Ott, H. & Zitha, P. L. J. μ-CT analysis and numerical simulation of drying effects of CO2 injection into brine-saturated porous media. Int. J. Greenh Gas Control. 1, 1 (2014).

Google Scholar

Pruess, K. & Muller, N. Formation dry-out from CO2 injection into saline aquifers: 1 effects of solids precipitation and their mitigation. Water Resour. Res. 45, 1–11 (2009).

Google Scholar

Berntsen, A. et al. Salt clogging during supercritical CO2 injection into a downscaled borehole model. Int. J. Greenh. Gas Control. 86, 201–210 (2019).

Article CAS Google Scholar

Hajiabadi, S. H. et al. Well Injectivity during CO2 geosequestration: A review of hydro-physical, chemical, and geomechanical effects. Energy Fuels 35, 9240–9267 (2021).

Article CAS Google Scholar

Zeidouni, M., Pooladi-Darvish, M. & Keith, D. Analytical solution to evaluate salt precipitation during CO2 injection in saline aquifers. Int. J. Greenh. Gas Control. 3, 600–611 (2009).

Article CAS Google Scholar

Azizi, E., & Cinar Y. Approximate analytical solutions for CO2 injectivity Into saline formations. SPE Reserv. Eval. Eng. [Internet]. 123–133 (2013).

Whittle, T., Park, P., & Coll, C. A semi-analytical model for the prediction of CO2 injectivity into saline aquifers or depleted hydrocarbon reservoirs. SPE Eur Energy Conf Featur 83rd EAGE Annu Conf Exhib. OnePetro (2022).

Md Yusof, M. A., Ibrahim, M. A., Mohamed, M. A., et al. Predictive modelling of CO2 injectivity impairment due to salt precipitation and fines migration during sequestration. Int Pet Technol Conf. OnePetro (2021).

Mardhatillah, M. K., Md Yusof, M. A., & Said, A. A., et al. Predictive modelling of carbon dioxide injectivity using SVR-hybrid. Offshore Technol Conf Asia. OnePetro (2022).

Yusof, M. A. M. et al. Combined impact of salt precipitation and fines migration on CO2 injectivity impairment. Int. J. Greenh. Gas. Control. 110, 103422 (2021).

Article Google Scholar

Kitron-Belinkov, M. et al. Groovy drops: Effect of groove curvature on spontaneous capillary flow. Langmuir 23, 8406–8410 (2007).

Article CAS PubMed Google Scholar

Tuller, M., Or, D. & Dudley, L. M. Adsorption and capillary condensation in porous media: Liquid retention and interfacial configurations in angular pores. Water Resour. Res. 35, 1949–1964 (1999).

Article ADS Google Scholar

Mason, G. & Morrow, N. R. Capillary behavior of a perfectly wetting liquid in irregular triangular tubes. J. Colloid Interface Sci. 141, 262–274 (1991).

Article ADS CAS Google Scholar

Idorenyin, E. & Shirif, E. A new simulation model for two-phase flow in porous media. Braz. J. Pet. Gas. 6, 1–18 (2012).

Google Scholar

Cai, J. et al. Generalized modeling of spontaneous imbibition based on hagen-poiseuille flow in tortuous capillaries with variably shaped apertures. Langmuir 30, 5142–5151 (2014).

Article CAS PubMed Google Scholar

Gates, J. I., & Lietz, W.T. Relative permeabilities of California cores by the capillary-pressure method. Drill Prod Pract. American Petroleum Institute (1950).

Fatt, I. & Dykstra, H. Relative permeability studies. J. Pet. Technol. 3, 249–256 (1951).

Article Google Scholar

Van Brakel, J. Pore space models for transport phenomena in porous media review and evaluation with special emphasis on capillary liquid transport. Powder Technol. 11, 205–236 (1975).

Article Google Scholar

Shi, J. Q., Xue, Z. & Durucan, S. Supercritical CO2 core flooding and imbibition in Berea sandstone—CT imaging and numerical simulation. Energy Proc. [Internet]. 4, 5001–5008 (2011).

Article CAS Google Scholar

Dullien, F. A. L. & Dhawan, G. K. Characterization of pore structure by a combination of quantitative photomicrography and mercury porosimetry. J. Colloid Interface Sci. 47, 337–349 (1974).

Article ADS Google Scholar

Lin, C. & Slattery, J. C. Three-dimensional, randomized, network model for two-phase flow through porous media. AICh.E J. 28, 311–324 (1982).

Article CAS Google Scholar

Yanuka, M., Dullien, F. A. L. & Elrick, D. E. Serial sectioning and digitization of porous media for two-and three-dimensional analysis and reconstruction. J. Microsc. 135, 159–168 (1984).

Article Google Scholar

Peysson, Y. Permeability alteration induced by drying of brines in porous media. Eur. Phys. J. Appl. Phys. 60, 12 (2012).

Article Google Scholar

Jeddizahed, J. & Rostami, B. Experimental investigation of injectivity alteration due to salt precipitation during CO2 sequestration in saline aquifers. Adv. Water Resour. 96, 23–33 (2016).

Article ADS CAS Google Scholar

Download references

The authors are grateful to Prof. Emmanuel A. Donkor, the Head, Department of Petroleum Engineering, Kwame Nkrumah University of Science and Technology, Ghana for his support.

Department of Petroleum Engineering, Kwame Nkrumah University of Science and Technology, PMB Kumasi, Kumasi, Ghana

Yen Adams Sokama-Neuyam, Shadrack Kofi Owusu, Victor Darkwah-Owusu, Joshua Nsiah Turkson & Adwoa Sampongmaa Otchere

Department of Petroleum Engineering, Universiti Teknologi PETRONAS, Seri Iskandar, 32610, Perak, Malaysia

Muhammad Aslam Md Yusof

Department of Energy and Petroleum Engineering, University of Stavanger, 4036, Stavanger, Norway

Jann Rune Ursin

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Y.A.S.-N. designed the work and wrote the main manuscript, M.A.M.Y. analysed the data and edited the manuscript, S.K.O. and A.S.O. conducted the experiments, V.D.-O. and J.N.T. worked on the theoretical modelling, J.R. Ursin supervised the work.

Correspondence to Yen Adams Sokama-Neuyam.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Sokama-Neuyam, Y.A., Yusof, M.A.M., Owusu, S.K. et al. Experimental and theoretical investigation of the mechanisms of drying during CO2 injection into saline reservoirs. Sci Rep 13, 9155 (2023). https://doi.org/10.1038/s41598-023-36419-3

Download citation

Received: 13 December 2022

Accepted: 03 June 2023

Published: 06 June 2023

DOI: https://doi.org/10.1038/s41598-023-36419-3

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.