Two-phase flow visualization under reservoir conditions for highly heterogeneous conglomerate rock: A core-scale study for geologic carbon storage

Geologic storage of carbon dioxide (CO2) is considered a viable strategy for significantly reducing anthropogenic CO2 emissions into the atmosphere; however, understanding the flow mechanisms in various geological formations is essential for safe storage using this technique. This study presents, for the first time, a two-phase (CO2 and brine) flow visualization under reservoir conditions (10 MPa, 50 °C) for a highly heterogeneous conglomerate core obtained from a real CO2 storage site. Rock heterogeneity and the porosity variation characteristics were evaluated using X-ray computed tomography (CT). Multiphase flow tests with an in-situ imaging technology revealed three distinct CO2 saturation distributions (from homogeneous to non-uniform) dependent on compositional complexity. Dense discontinuity networks within clasts provided well-connected pathways for CO2 flow, potentially helping to reduce overpressure. Two flow tests, one under capillary-dominated conditions and the other in a transition regime between the capillary and viscous limits, indicated that greater injection rates (potential causes of reservoir overpressure) could be significantly reduced without substantially altering the total stored CO2 mass. Finally, the capillary storage capacity of the reservoir was calculated. Capacity ranged between 0.5 and 4.5%, depending on the initial CO2 saturation.

Atmospherically released carbon dioxide (CO 2 ) is considered to be a major deriver behind climate change, and as such geologic CO 2 storage is considered a key technology in climate change mitigation strategies 1,2 . Both industry and research communities are currently evaluating the safety and feasibility of long-term CO 2 sequestration, and a number of pilot-and demonstration-scale projects have been conducted as a part of this effort to test, monitor, and verify technologies in various subsurface geological environments [3][4][5] .
Previous research related to geologic CO 2 storage has until now focused on evaluating sandstone formations, as their relatively high porosity and permeability suggest greater economic viability than other rock formations. For this reason, the reservoir lithology of most pilot-scale (e.g., Frio, Nagaoka, Ketzin and Otway) and demonstration-or commercial-scale projects (e.g., MGSC Decatur, Sleipner, Snøhvit, In Sala and Gorgon) has been sandstone 3,5 . The average reservoir porosity has been 5-35% and reservoir permeabilities have ranged from as low as 5 mD (In Salah) to 5,000 mD (Sleipner).
One exception was the MRSCP Gaylord project (a Midwest regional carbon sequestration partnership), where the targeted storage formation consisted of dolomite. This was characterized by interbedded, laminated algal dolomudstone, represented by a mean porosity and permeability of 13% and 22.6 mD, respectively 6 . In an effort to assess the geologic storage potentials, basalt formations have also been considered as alternative storage formations 7,8 .
Although a number of CO 2 storage projects are being conducted worldwide, the need for additional safe storage verification experiments using various geological formations cannot be overemphasized. The fact that only limited data has been collected during the injection and post-injection phases demonstrates the necessity of additional field assessments of the processes leading to plume stabilization and long term trapping 9 . Followed by on-going international efforts, the first onshore pilot-scale (~10,000 ton) CO 2 storage project is being conducted in South Korea at the Janggi sedimentary formation, located on a south-eastern portion of the Korean Peninsula ( Fig. 1(a)). One distinctive feature differing from other pilot-scale projects is that the targeted formation is a conglomerate, consisting of gravel and cobbles in a silty sand matrix. Figure 1(b) shows a thickness-permeability cross-plot of values selected from various existing storage projects. This plot reveals the injectivity at Janggi falls between 0.1 and 1 darcy-meter, similar to the estimated injectivities at the Nagaoka and In Salah sites.
In addition to field-demonstration studies, well-controlled laboratory experiments are essential for obtaining a fundamental understanding of physical 10,11 and chemical processes 12,13 occurring in potential storage formations, as well as predicting their CO 2 trapping capacities 14 . A number of studies have conducted reactive transport experiments to understand the complex processes of CO 2 -water-mineral interactions and their coupling effects on the CO 2 -water multiphase flow at in-situ reservoirs conditions [15][16][17][18] . These experimental data can be utilized to validate the numerical studies and build-up confidence in the modeling tools. For example, Smith et al. 18 used detailed experimental data to constrain three-dimensional reactive transport models in order to describe and predict the evolution of pore-space and permeability during geologic CO 2 storage in carbonate reservoirs.
Core-flooding is a technique that has been used to conduct various laboratory experiments on core samples under conditions closely imitating those of subsurface environments 19 . Due to this advantage, core-flooding tests have been utilized for characterizing the dynamic behaviour of CO 2 under reservoir conditions. Such tests have focused on resolving various problems related to important factors controlling the behaviours of two-phase fluids, including permeability heterogeneity and multiphase flow [20][21][22] , measuring relative permeabilities [23][24][25] , interfacial tension 26,27 , and capillary pressures [28][29][30] . Additionally, the core-flooding experiments have been applied to assessing residually trapped CO 2 amounts under reservoir conditions 23,25,[31][32][33] . Still other research topics have included multiphase fluid behaviour in fractured porous media 34,35 , salt-precipitation dynamics 36,37 and the non-equilibrium dissolution of CO 2 in a heterogeneous core 32,38 .
Despite these numerous core-scale studies, there has been little investigation on how CO 2 migrates in a highly heterogeneous porous media. This study presents, for the first time, a two-phase flow visualization under reservoir conditions for a highly heterogeneous conglomerate core that has been obtained from the real CO 2 storage site. The three-dimensional (3D) distribution of clasts within the matrix was determined utilizing X-ray computed tomography (CT). Additionally, real-time X-ray scanning techniques (conducted during the CO 2 and brine flow tests) captured the dynamic distribution of CO 2 saturation in the conglomerate core. Multiphase transport simulations were also conducted to assess the effects of heterogeneity on pressure build-up as well as spatial variation of CO 2 saturation. Finally, based on the initial-residual (IR) characteristic curve, the storage potential of the conglomerate formation is discussed.

Results and Discussion
Heterogeneity characterization. A core-plug obtained from the Janggi conglomerate formation was used to investigate the spatial distribution of clasts within a silty sand matrix. Figure 2(a) presents two cross-sectional CT images taken along the longitudinal axis at different angles. In addition, Fig. 2(b) shows the binary image generated from the 8-bit CT images, where the boundaries between the clasts and matrix are distinguished more clearly and discontinuities within the clasts are highlighted. The clast fractioning profiles along the longitudinal core-axis are shown in Fig. 2(c). Nine cross-sectional images, taken along the white-dashed lines appearing in   Fig. 2(d). In the upstream region (30 mm) of the core, the constituents are relatively homogeneous and only a small fraction of clasts were preserved (f < 0.3). In the midstream portion (70 mm), both medium size of clasts and matrix existed together, whereas larger clast sizes were present in the downstream region (100 mm). Figure 1(e) shows the three-dimensional distribution of clasts in the core-plug.
The pore size distribution were analysed for both the matrix and clasts ( Fig. 3(a)). The clast regime was characterized by a skewed uni-modal pore-size distribution dominated by micropores (<1 μm diameter). The silty sand matrix regime demonstrated a bi-modal pore size distribution dominated by micropores (<1 μm diameter)  and mesopores (>2 μm diameter) ( Fig. 3(a)). The bi-modal pore size distribution characteristics of the matrix were also reflected in the shape of the capillary pressure (P c ) curve, appearing as a superposition of two uni-modal pore size distribution characteristic curves. The P c (S w = 0.5) of the clasts was more than two orders of magnitude greater than that of the matrix ( Fig. 3(b)).

Multiphase flow experiments.
This experimental study aimed to characterize the spatio-temporal distribution of CO 2 saturation build-up during the multiphase flow tests. Two different tests were conducted: (i) CO 2 injection into a brine-saturated core, simulating the CO 2 injection period (drainage) conditions; and (ii) brine injection after the CO 2 injection test, reproducing the post-injection period (imbibition) conditions. CO 2 saturation distribution snapshots revealed the propagation of CO 2 saturation throughout the conglomerate core with increases in pore volume (PV) (Supplementary Figs S1 and S2). Figure 4(a) presents the vertically averaged CO 2 saturations from the five different locations (30, 60, 70, 110, and 120 mm) delineated in Fig. 2(a). Although the averaging is smoothing out the fluctuation of CO 2 saturation related to small-scale heterogeneity, it clearly presents the increasing tendency within the core sample. All profiles showed logarithmic increases with the PVs. The CO 2 saturation profile at 30 mm increased slowly when approaching the asymptote during the flow tests; however, the other profile (such as at 110 mm) approached their asymptotes more rapidly, implying that CO 2 saturation was quickly stabilized throughout the experiment. It is interesting to note that CO 2 saturation at the downstream side of the core was greater than that at the upstream side in early stages (<1.0 PV). The spatial inversion of CO 2 saturation distribution was attributed to the heterogeneous pore space along the conglomerate core, including the presence of complex discontinuities in clasts ( Fig. 2(b)). On the upstream side where the matrix was dominant, CO 2 saturation was low due to the presence of relatively large sized pores, which accommodated the injected CO 2 . Equal amounts of CO 2 flowing through the downstream side, dominated by clasts, showed elevated CO 2 saturation because clast discontinuities provided only small pore-spaces.
To perform a detailed assessment of the spatial distribution of CO 2 saturation, three subdomains (20 × 20 mm) were selected, one each located in the upstream (domain I), midstream (domain II), and downstream (domain III) regions as shown in Fig. 4(b). These domains primarily composed of sand, small clasts, and large clasts, respectively. Domain I (primarily sand) was characterized by a relatively homogeneous CO 2 saturation distribution over time (Fig. 4(c)) with a small (sub-core) scale heterogeneous CO 2 saturation observed due to the effect of the bi-modal pore-size distribution preserved in the matrix (Fig. 3(a)). Domain II, which was composed of irregularly spaced small clast fragments, demonstrated heterogeneous CO 2 saturation caused by irregularly spaced small clast fragments. Low CO 2 saturations resulted from the presence of clasts, whereas high CO 2 saturations were observed in high porosity sand. Highly heterogeneous CO 2 saturation was observed in domain III, caused by a dense network of discontinuities within the large clasts where multiple discontinuities with relatively large apertures (110 mm) were located, as seen in Fig. 4(c). This discontinuity network provided well-connected pathways for CO 2 flow, resulting in a rapid increase in CO 2 saturation.
The vertically averaged CO 2 saturation was also plotted for the brine injection tests (Fig. 4(d)). As discussed in Oh et al. 32 , the dynamic transition of CO 2 trapping could be assessed during the imbibition test. Three regimes of CO 2 trapping mechanisms were characterized: (i) immediate mobile CO 2 displacement by injected brine (R1), (ii) immobile CO 2 preservation as residual trapping (R2), and (iii) the gradual dissolution of residually trapped CO 2 into the fresh brine (R3). The dynamic transition of CO 2 trapping at 30 mm displayed initial decrease in CO 2 saturation from 0.35 to 0.27 during the early stage (0.01-0.1 PV) primarily resulted from the mobile CO 2 displacement (Regime 1). In Regime 2, the CO 2 saturation was relatively consistent and stabilized at 0.27 between 0.1 and 1.0 PV, which implies the preservation of immobile CO 2 (residually trapped). This is possible because the brine that enters the core will initially dissolve a small amount of CO 2 at the inlet of the core and becomes fully-saturated preventing it from dissolving any more CO 2 . However, as more fresh brine was injected into the core, the residually trapped CO 2 dissolved into fresh brine (Regime 3). Figure 4(e) shows CO 2 saturation snapshots from the three different domains during the brine injection test. In domain I where the matrix dominated, CO 2 was distributed uniformly at the end of drainage test. Brine displacement tests induced similar uniform distributions in residually trapped CO 2 , which eventually dissolved into the brine. However, in domains II and III, a larger fraction of pores was initially filled with CO 2 , resulting in more snap-off and trapping as brine invaded the pore space. Under these conditions, CO 2 dissolution was highlighted only within the matrix and discontinuities.
Multiphase transport simulations. The above experiments were conducted with fluid injection at the left boundary, but CO 2 injection from the opposite boundary would cause difference in the distribution of CO 2 saturation due to nature in highly heterogeneous conglomerate core. We used TOUGH2 39 simulator to model multiphase and multicomponent fluid transport to assess directional effect of fluid flow on pressure build-up as well as CO 2 saturation distribution along the core sample. A two-dimensional model was constructed based on the CT image of the core sample to reflect the distribution of clasts. Then, two different cases were simulated (case 1: left-to-right flow; case 2: reverse flow). Figure 5 presents simulation results of pressure (a and f) and CO 2 saturation map (b and g), profiles for CO 2 pressure (c and h) and saturation (d and i) along the center of the model domain, and pressure build-up (e and j) for the two different cases. Both cases showed a significant pressure drop between points A and B, which lie on a boundary of a large sized clast. The heterogeneous pressure distribution resulted in different CO 2 saturation within the model for the two cases. For example, case 1 showed relatively high CO 2 saturation (S g~0 .5) at point A and low CO 2 saturation (S g~0 .3) at point B, while case 2 showed relatively high CO 2 saturation (S g~0 .5) at point B. This implies that the distribution of CO 2 saturation is dependent on the flow direction. The pressure difference between the inlet and outlet (ΔP) showed ~150 kPa during the early stage of injection and stabilized at ~30 kPa for both cases (Fig. 5(e and j)).

CO 2 mass at different capillary numbers.
A dimensionless capillary number (N c ), representing the relative contributions of viscous versus capillary forces, was set to reflect the reservoir conditions. Two CO 2 injection tests (q co 2 = 1.0 and 0.1 ml/min) were performed in this study, producing associated N c values of ~10 −1 and ~10 −2 , respectively. The high N c case was observed in a transition regime between the capillary and viscous limits, whereas the low N c case was under capillary dominated flow conditions. The transient pressure difference (ΔP) values across the core for the two different cases are shown in Fig. 6. The ΔP build-up under a high injection rate (N c ~ 10 −1 ) reached 2,300 kPa in early injection stages before stabilizing to ~700 kPa; however, the low injection rate (N c ~ 10 −2 ) displayed a small buildup of ΔP build-up (~150 kPa) in early injection stages before ΔP stabilized to ~65 kPa for the remainder of the multiphase flow experiment. Figure 7 displays the CO 2 saturation distribution maps for the two injection rates (q co2 = 1.0 and 0.1 ml/min) and the CO 2 saturation contrast (ΔSco 2 = S q=1.0 −S q=0.1 ) maps between the high and low injection rate cases. The  measured ΔSco 2 displayed a strong dependency on core heterogeneity, showing negative values on the upstream side (implying lower CO 2 saturation with a greater injection rate (S q=1.0 < S q=0.1 )). The high injection rate created excessed pressure on CO 2 (nonwetting phase), enabling CO 2 to transport across the high P c region on the midstream side. This had a negative impact on the build-up of CO 2 mass on the upstream side. At low injection rates, however, additional CO 2 occupied the pore-spaces on the upstream side because CO 2 transport across the high P c region on the midstream side was hindered. The CO 2 mass within the core was calculated for both N c cases to be 2.24 and 2.18 g, respectively, once a steady-state condition was reached. This demonstrated a 10× reduction in injection rate caused only a 2% decrease in CO 2 mass, suggesting that injection rates, which may cause reservoir overpressure when high, can be significantly reduced without substantially altering the total CO 2 mass. Storage capacity. During the post-injection stage, capillary forces at the pore scale are responsible for breaking up the CO 2 phase into the form of blobs or ganglia, becoming effectively immobile 40,41 . The capillary trapping potential was evaluated in these conditions based on the performed multiphase flow experiments. The initial-residual (IR) characteristic curves have commonly been used to describe the residual trapping potentials of rock, with numerous developed models describing this relationship. This study applied the Land model 42 , one of the earliest and most widely used trapping models. Figure 8(a) presents the relationship between the initial (Sco 2,i ) and residual saturations (Sco 2,r ), where Sco 2,i refers to the initial CO 2 saturation prior to brine injection testing whereas Sco 2,r refers to the residually trapped CO 2 saturation. For analysis, highly heterogeneous conglomerate core was divided into 50 equal subdomains. Then, based on the brine injection testing, the residual CO 2 saturation achieved at 0.1 PV (R2) was determined (Fig. 4(d)), and subsequently, the amount of capillary trapped CO 2 was assessed using the Land model given by . Here, C is the dimensionless constant known as the Land coefficient. In this study, the calculated C was calculated as 0.9, showing a monotonic increase in the residual saturation as a function of Sco 2 , i . This was in agreement with numerous previous works investigating water-wet porous medias.
The capillary-trapping capacity, defined as ϕSco 2,r (ϕ is porosity) is a rescaling of the traditional IR saturation plot 43,44 . This quantity is of great importance in geologic CO 2 storage applications because it states how much CO 2 can be stored securely per unit volume of a rock. Figure 8(b) presents the relationship between Sco 2,i and ϕSco 2,r at different porosities (blue, red and orange circle symbols). The IR relationships of other core samples (Berea and Pohang sandstones) are plotted alongside for comparison. Core sample details are described in the Supplementary Table S1. The trapping capacities of both Berea and Pohang sandstones showed more or less linear increases with Sco 2,i , whereas the ϕSco 2,r values for heterogeneous conglomerate cores were widely spread and dependent on the size of porosity. Although the overall capacity of the conglomerate core was smaller than those of the Berea or Pohang sandstones, the conglomerate high porosity region displayed a similar capacity to the sandstones (e.g., ϕSco 2,r ~4.5% at Sco 2,i = 0.3).

Conclusions
Although sandstone formations are considered to be the most appropriate for CO 2 sequestration, other types of geologic formations are being investigated for this purpose. For example, the target rock formation in MRSCP Gaylord project was dolomite 45 , and the Midale project (an extension of the Weyburn project) demonstrated the possibility for CO 2 containment in a fractured reservoir 5 . Basaltic formations have also been considered as alternative storage formations 7,8 . In South Korea, a conglomerate formation is being tested for geologic CO 2 storage. Figure 8. (a) A trapping curve relating initial CO 2 saturation (Sco 2,i ) to residual CO 2 saturation (Sco 2,r ). The Sco 2,i refers to the initial CO 2 saturation prior to imbibition test while the Sco 2,r refers to the residually trapped CO 2 saturation. Red circles represent the relationship between the Sco 2,i and Sco 2,r within equally divided subdomains; error bars represent one standard deviation. (b) Trapping capacity (ϕSco 2,r ) versus initial CO 2 saturation (Sco 2,i ); colors represent different ranges of porosity. This study reported a thorough core-scale investigation of a highly heterogeneous conglomerate formation. The rock heterogeneity was investigated with the use of X-ray CT and MICP analysis revealing that the targeted storage formation was composed of a silty sand matrix and irregularly spaced clasts of various sizes. Laboratory multiphase flow tests with flow visualization under reservoir conditions provided important information on CO 2 migration behaviours as well as the storage capacity of the highly heterogeneous conglomerate. Variations in clast spatial configurations gave rise to different storage features. Dense networks of discontinuities within clasts were observed to provide well-connected CO 2 flow pathways, causing a rapid increase in CO 2 saturation, and potentially helping to reduce overpressure. Two flow tests, one in the capillary-dominated condition and the other in a transition regime between the capillary and viscous limits, revealed that a 10× reduction in injection rate caused only a 2% decrease in the CO 2 mass. This implies the heterogeneity positively influenced pressure management because high injection rates, which may cause overpressure in the reservoirs, could be significantly reduced without substantially altering the total stored CO 2 mass. We note that although chemical reactions are not considered in this study, the complex process of CO 2 -water-mineral interactions may play important roles on multiphase flow in subsurface reservoirs and trapping mechanisms.
Among the four potential trapping mechanisms, capillary trapping is a rapid and effective mechanism for rendering an injected fluid immobile and reducing the need to ensure caprock integrity. The capillary-trapping capacity of the conglomerate core investigated in this study was calculated using imbibition test results and the IR characteristic curve. While the Berea and Pohang sandstone cores showed approximately linear increases in capacity with initial CO 2 saturation, the Janggi conglomerate displayed significant variations due to its highly heterogeneous porosity distribution. The conglomerate capacity ranged between 0.5 and 4.5%, depending on the initial saturation (0-0.4): for an initial saturation of 0.4, the conglomerate capacity ranged between 3-4%. Despite the fact that the trapping analyses showed lower overall capacities relative to homogeneous sandstone, the high porosity region of the conglomerate core displayed a similar capacity to sandstones, implying conglomerate formations could be an alternative for CO 2 storage formation. CO 2 injection and comprehensive monitoring in conglomerate formations will validate the feasibility of CO 2 storage in the highly heterogeneous conglomerate formations.

Methods
Field site and geological setting. An onshore pilot-scale geologic CO 2 storage project was launched in South Korea in 2011. The primary tasks (Phase I, 2012-2015) were the site-selection and initial characterization of potential pilot-scale storage sites. The Janggi conglomerate, composed of both clasts and silty sand matrix, was chosen as the target storage formation for the project. More recently, Phase II (2015-2017) was completed, including the drilling of seven boreholes to investigate the lithology of targeted and capping formations. The representative lithological log for this is illustrated in Fig. 1(a) 46 . The average porosity and permeability of the storage formation were 16.1% and 8.5 mD, respectively 47 . The overlying Seongdongri formation, composed of thick mudstone and massive dacitic tuffs, forms the regional sealing formation for CO 2 storage, and a monitoring well (~40 m up-gradient from the injection well) was drilled. A comprehensive monitoring system was installed including a down-hole pressure and temperature measurements, vertical seismic profile (VSP), electrical resistivity tomography (ERT), distributed temperature sensing (DTS) and distributed acoustic sensing (DAS) along the well via an optical fiber, and U-tube downhole sampler. X-ray micro CT. The dimensions of the core-plug obtained from the Janggi conglomerate formation were 130 × 48 mm (length × diameter). X-ray computed tomography (CT), a non-destructive and non-invasive method, was used to investigate the spatial distribution of clasts within the silty sand matrix. The CT images used in this study were obtained at 140 kV and 200 μA. A total of 2,000 cross-sectional images were taken at 0.065 mm intervals with a micrometer-order pixel size. Subsequently, two-dimensional (2D) images were stacked to generate 3D core sample images.
Multiphase flow tests with in-situ imaging technology. The apparatus for core-flooding experiments was designed to conduct multiphase flow tests. It was comprised of fluid (supercritical CO 2 and brine) injection, core-holder, temperature controlling, confining, and back-pressure control systems (Supplementary Fig. S5). A 2D X-ray scanning system, consisting of an X-ray tube and a detector, was integrated into the core-flooding system in order to capture the real-time distribution of CO 2 saturation distributions along the core sample. The detector board had two arrays of 128 channels each, corresponding to 256 detectors with a detector pitch of 0.4 mm. The array of the X-ray detector arrays were perpendicular to the core-axis, performing scans parallel to core-axis. The data output was in 16-bit format, ranging between 0 and 65,535. The produced 2D images were converted into CO 2 saturations using the following equation: , where Ico 2 and I brinesat are the values of grayscale intensity obtained from the background scans of the core, which were saturated with CO 2 and brine in the core, respectively. Finally, I exp was the value obtained for the multiphase condition. The obtained CO 2 saturation map may include noise mainly due to the relatively small density contrast between the two fluid phases (CO 2 and brine) compared with the contrast between any one of the fluids and the rock 48 . In the experiments, NaI was used as a dopant to enhance the X-ray attenuation so that the CO 2 saturation could be distinguished from the brine in the scanned image. Two multiphase flow tests were conducted: (i) CO 2 injection into a brine-saturated core, simulating the CO 2 injection period (drainage) conditions; and (ii) brine injection following the CO 2 injection test, reproducing post-injection period (imbibition) conditions. Downstream pressure was maintained during these two experiments at 10 MPa, and the temperature was set at 50 °C to replicate the subsurface environment at Janggi pilot site. L and H represent the length and diameter of the core, respectively; ΔP represents the pressure drop; and ΔP c was calculated from the capillary pressure curves), were ~10 −1 and ~10 −2 , for q = 1.0 and 0.1 ml/min, respectively.
Pore size distribution and capillary pressure. The pore size distributions and capillary pressure curves for the silty sand matrix and clasts in the conglomerate sample were determined using a mercury injection porosimeter (MICP). To perform MICP tests, two small subsamples (~1 cm 3 ) were taken from the conglomerate core each representing either the silty sand matrix or clasts. These subsamples were dried in a vacuum oven at 70 °C for 24 hours prior to testing. The MICP and corresponding mercury saturation measured in the mercury/air system were transformed into a capillary pressure and corresponding CO 2 saturation for a CO 2 /brine system using the following relationship: , where P c (Pa) represents the capillary pressure, σ (mN m −1 ) represents the interfacial tension (IFT) between two fluids, and θ (deg) represents the contact angle measured in the wetting phase. The IFT of σ Hg = 485 mN m −1 and σ CO 2 =32.6 mN m −1 was used in the conversion for the geologic CO 2 sequestration conditions. The contact angles for the two systems were assumed to be equal 29 .