Unlock the Potentials to Further Improve CO 2 Storage and Utilization with Supercritical CO 2 Emulsions When Applying CO 2 -Philic Surfactants

: Supercritical CO 2 (ScCO clearly demonstrate the potentials of CO 2 emulsion on CO 2 storage and oil recovery over conventional tertiary productions. Relative to traditional aqueous soluble surfactant emulsion, the novel surfactant emulsion contributes to higher injectivity, CO 2 storage capability, oil recovery and energy utilization efﬁciency. The CIDS could further reduce water injection cost and energy consumption. The ﬁndings here reveal the potentials of further improving CO 2 storage and utilization when applying ScCO 2 -philic surfactant emulsion, to compromise both environmental and economic concerns. inch-thick rubber sleeve to prevent leakage and axial absolute and differential transducers pressure variations and propagation. differential transducers 1 (2 inches long), (4 (4 long), respectively. Two BPRs are connected in series to the outlet of the at and 1100 psi respectively, to realize stepwise pressure reduction and prevent severe pressure shock during CO 2 phase conversion. The employed are from outcrops, by and under 110 ◦ for CO 2 diffusion and penetration, with multiple aluminium then heat 12

As mentioned, even though CO 2 emulsion with ASS has been implemented in multiple field applications, the performances could be dramatically improved by the employment of novel CPS. Up to now, however, there has been a lack of systematic comparison studies to demonstrate the superiorities of CPS caused by its dual phase partition capacity, which is addressed here. In current communication, comprehensive comparisons are carried out experimentally and numerically with a commercially available ASS (CD-1045) and a proprietary nonionic CPS whose aqueous stability, partition coefficient between ScCO 2 /brine and solubility in ScCO 2 have been determined in prior publications [61]. Individual surfactant affinity to the oleic phase is indicated through a phase behavior test. To gain some insight on the electro-kinetic behaviors of selected surfactants, static adsorption measurements are carried out on Silurian dolomite outcrop. CO 2 storage and retention in porous media are indicated by corefloodings with both Berea sandstone and Silurian dolomite under a twophase condition (W/ScCO 2 ) to evaluate individual emulsion performances, characterized by incubation times, emulsion propagation, pressure drop magnitude and liquid desaturations. Harsher conditions are employed to reveal the incremental oil production over CO 2 flooding and co-injection of brine/ScCO 2 through fractured carbonate. Moreover, field scale simulations are employed to predict the superiorities of ScCO 2 foam with CPS under a gravity field. Multiple tertiary injections act as baselines (ScCO 2 and WAG) to illustrate the necessity of conformance control when mitigating gravity segregation. SAG injections are carried out to compare two types of CO 2 emulsions, which are assessed by phase saturation, surfactant/pressure distribution, CO 2 storage/retention factor, oil production and the CO 2 utilization ratio. The studies here attempt to meet the gap between academic research and field implementation of this novel technology, which could benefit both environmental concerns via the final sink of CO 2 even in the form of gas hydrates [62,63], as well as industrial revenue through substantial increased oil recovery for waste valorization towards a circular economy [64].

Materials
Synthetical brine with 3 wt% pure NaCl (99.99%) is used to be consistent with the conditions employed for solubility and partition coefficient measurements [61], as too do the system pressure (1500 psi) and temperature (35 • C). The viscosity and density of Wasson crude oil are 7.1 cp (@35 • C) and 0.86 g/cm 3 respectively. Chaser CD 1045 (anionic) is commercially available without modifications. A novel proprietary nonionic hydrocarbon CPS (2-ethyl-1-hexanol, 2EH-PO 5 -EO 15 ) is selected, whose solubility in ScCO 2 and partition coefficient between brine/ScCO 2 are shown in the Appendix A, as previously determined [61]. Both Silurian dolomite and Berea sandstone outcrops are employed, and properties in details are displayed below (Table 1). Industrial-grade liquid CO 2 , delivered in cylinders with 900 psig initial pressure, is compressed to 1500 psi.

Experimental Apparatus and Procedures 2.2.1. Static Adsorption
In current studies, except for surfactant adsorption, there are no other geochemical reactions of ScCO 2 with carbonate accounted. Static surfactant adsorption on Silurian dolomite outcrop is conducted at 35 • C and pH of 7.2. With two sieves of 40 and 100 mesh, fine particles ground from the rock sample, whose diameters span 0.15 mm to 0.425 mm, are used. In every test, 10 g of prepared fine particles are mixed with 20 g aqueous surfactant solution. The composed sample is rested for 24 h after being gently stirred for 1 h, and centrifuged for 5 min to separate the particle and solution. Then, high-performance liquid chromatography (HPLC) is employed to analyze the surfactant concentration in the supernatant. Even though it is known that static adsorption may overestimate the surfactant retention over a dynamic test [65], the studies here attempt to pay more attention on comparative evaluation rather than individual performance.

Phase Behavior Test
It is important to understand the preferences to interfaces of employed surfactants with presence of the third phase (oleic phase), owing to their amphiphilic characteristics. Foam stability could be impaired, potentially caused by the mitigation of the surfactant from the G/W interface to the oleic phase due to preferential partition [66], which may also generate an additional microemulsion phase. The composition and structure of microemulsion and macroemulsion are determined by aqueous salinity, temperature and other factors [67]. Meanwhile, it is known that presence of microemulsion or macroemulsion may individually facilitate the increase of the desaturation capillary number through the reduction of interfacial tension [68] or an increase of viscous force [69]. However, those "side effects" would confound the observations' interpretations, and it is also uncertain that synergism or antagonism would appear with the addition of foam. Hence, the affinity to the O/W interface of a current foaming agent is indicated through a phase behavior test. Here, mother brine fluid is diluted to reach eight different salinities in 2 wt% increment. Thereafter, the surfactant is added respectively to achieve 0.2 wt% solutions. The opening of the 5 cc capacity pipette is sealed after mixing 2 cc surfactant solution with 1 cc oil. The fluid is mixed by rotating the pipette gently and frequently, and then the remainder in a convection oven under 35 • C for 10 days to reach equilibrium.

Two-Phase Coreflooding without Artificial Fracture
A schematic graph of coreflooding setup is shown in Figure 1, which comprises a fluid injection system, oven, core holder, pressure transducers, back pressure regulator (BPR) and effluent collector. Brine and surfactant are injected by a TELEDYNE ISCO Model 500D syringe pump, while compressed ScCO 2 is displaced by DI water through a high-pressure accumulator separated by a piston and connected to a QX-6000 Quizix pump. A phoenix Hassler-type core holder for a 1 ft-long & 2-inch diameter core is mounted vertically, and fluids are injected from the top. Overburden pressure is exerted by hydraulic oil assisted by the 0.25 inch-thick rubber sleeve to prevent leakage and ensure the axial flow. Two absolute and three differential transducers (±0.1 psi) align the side of the core holder to indicate pressure variations and fluid propagation. Three differential transducers correspond to the three sections arraying downward, designated as Section 1 (2 inches long), 2 (4 inches long) and 3 (4 inches long), respectively. Two BPRs are connected in series to the outlet of the core holder, set at 1500 psi and 1100 psi respectively, to realize stepwise pressure reduction and prevent severe pressure shock during CO 2 phase conversion. The employed cores are drilled from outcrops, followed by cleaning and drying under 110 • C for 24 h. To prevent CO 2 diffusion and penetration, cores are wrapped with multiple layers of aluminium foil, and then a teflon heat shrink tube. After loading into the core holder, it is vacuumed for 12 h before porosity and permeability measurement with brine, which are listed in Table 1 for details. Cores are not reused among different corefloodings.
Sustain. Chem. 2021, 2, FOR PEER REVIEW mounted vertically, and fluids are injected from the top. Overburden pressure is exe by hydraulic oil assisted by the 0.25 inch-thick rubber sleeve to prevent leakage and sure the axial flow. Two absolute and three differential transducers (±0.1 psi) align side of the core holder to indicate pressure variations and fluid propagation. Three d ential transducers correspond to the three sections arraying downward, designate Section 1 (2 inches long), 2 (4 inches long) and 3 (4 inches long), respectively. Two B are connected in series to the outlet of the core holder, set at 1500 psi and 1100 psi res tively, to realize stepwise pressure reduction and prevent severe pressure shock du CO2 phase conversion. The employed cores are drilled from outcrops, followed by c ing and drying under 110 °C for 24 h. To prevent CO2 diffusion and penetration, core wrapped with multiple layers of aluminium foil, and then a teflon heat shrink tube. A loading into the core holder, it is vacuumed for 12 h before porosity and permeab measurement with brine, which are listed in Table 1 for details. Cores are not re among different corefloodings. In-situ generation of CO2 emulsion is employed without a pre-generator thro simultaneous injection of ScCO2 and surfactant solutions on the initially fully brine-s rated cores. A 0.2 wt% surfactant concentration is used along all of the experiments a dissolved surfactant in the aqueous phase regardless of the surfactant nature. Con fluid injection rates (CO2, 0.1 cc/min, 1.4 ft/d; liquid, 0.3 cc/min, 4.2 ft/d) are maintain achieve 75% foam quality (injection gas fraction) at reservoir conditions (1500 psi, 35 Pressure drops are recorded in real time, and the average water saturation inside the is derived from the mass balance calculations of effluents.  In-situ generation of CO 2 emulsion is employed without a pre-generator through simultaneous injection of ScCO 2 and surfactant solutions on the initially fully brine-saturated cores. A 0.2 wt% surfactant concentration is used along all of the experiments with a dissolved surfactant in the aqueous phase regardless of the surfactant nature. Constant fluid injection rates (CO 2 , 0.1 cc/min, 1.4 ft/d; liquid, 0.3 cc/min, 4.2 ft/d) are maintained to achieve 75% foam quality (injection gas fraction) at reservoir conditions (1500 psi, 35 • C). Pressure drops are recorded in real time, and the average water saturation inside the core is derived from the mass balance calculations of effluents. The same coreflooding setup is used as above ( Figure 1). The same procedures are employed to prepare the unfractured core, and matrix porosity & permeability are measured in the same manner (Table 1). Thereafter, the core is fractured artificially with a band saw into two halves along the longitudinal axis, followed by drying out in an oven for 24 h, which creates a smooth surface with minimal roughness [51]. Then, three precut Teflon film strips with 75-micron thickness are placed onto the surface to build an artificial fracture whose contributes negligibly (<0.1%) on total pore volume (PV) (Figure 2a). The split halves are composed back and fixed with tape ( Figure 2b). Similarly, brine saturation in 1 cc/min stops after 10 PV injections under 2000 psi confining pressure. The permeability to brine of the fractured composite core is measured and derived with Darcy's law, as listed in Table 1. 2 PV of crude oil is injected under constant pressure (400 psi), and the initial oil saturations before each following tested strategy are also documented in Table 1. The same coreflooding setup is used as above ( Figure 1). The same procedures are employed to prepare the unfractured core, and matrix porosity & permeability are measured in the same manner (Table.1). Thereafter, the core is fractured artificially with a band saw into two halves along the longitudinal axis, followed by drying out in an oven for 24 h, which creates a smooth surface with minimal roughness [51]. Then, three precut Teflon film strips with 75-micron thickness are placed onto the surface to build an artificial fracture whose contributes negligibly (<0.1%) on total pore volume (PV) (Figure 2a). The split halves are composed back and fixed with tape ( Figure 2b). Similarly, brine saturation in 1 cc/min stops after 10 PV injections under 2000 psi confining pressure. The permeability to brine of the fractured composite core is measured and derived with Darcy's law, as listed in Table 1. 2 PV of crude oil is injected under constant pressure (400 psi), and the initial oil saturations before each following tested strategy are also documented in Table 1. Pure ScCO2 injection is carried out at 0.8 cc/min (2.333 ft/D) on composite fractured core as a baseline, until oil cut in effluent is close to 0, to indicate the impacts of such extreme heterogeneity even though minimal miscible pressure is satisfied. The simultaneous injection of brine and ScCO2 at 0.2 and 0.6 cc/min respectively at reservoir conditions (1500 psi and 35 °C) to maintain 75% injection foam quality behave as the other baseline. Then, CO2 emulsion corefloodings are conducted similar to the above W&G co-injection, except for an additional 0.2 wt% surfactant in the brine. The pre-saturation of the surfactant solution is not employed to resemble the practical conditions in the field. Outcrop cores are not reused, and properties such as initial oil saturation are listed in Table 1. Pressure drops across the core and oil production with time are recorded.

Simulation Deck
As mentioned, foam behaviors are sensitive to in-situ phase velocities, which implies the radial flow in a cylindrical coordinator would be more appropriate for foam modeling, and better resemble the practical scenario. A 2D sector (theta 15°, porosity 20%) simulation deck is built, as shown in Figure 3, which is homogeneous but anisotropic with respect to permeability, except for the outermost grid column. The permeabilities in radial and vertical directions of the main body are set as 200 and 400 md, respectively. In the outermost column, the vertical producer locates holds permeability as high as 10,000 Darcy in every direction to mimic an open boundary condition and prevent artificial backflow [70]. An injector is positioned in innermost column grids vertically. The sector is 440 ft in a radial direction, which is divided into 100 grids unevenly. The grid size is 3 ft for the first 30 grids, and increases to 5 ft for the remains. It is as thick as 100 ft in the Z direction, which is cut to 20 layers evenly. Initial oil saturation is 0.4 without the presence of a gaseous phase and any surfactants. Oil properties have been tuned with CMG/WINPROP to ensure the achievement of multiple contact miscible (MMP~1082 psi) under initial reservoir conditions (1500 psi, 35 °C). Measured partition coefficients between brine/ScCO2 of selected CPS (2EH-PO5-EO15) is fed into STARS directly. Pure ScCO 2 injection is carried out at 0.8 cc/min (2.333 ft/D) on composite fractured core as a baseline, until oil cut in effluent is close to 0, to indicate the impacts of such extreme heterogeneity even though minimal miscible pressure is satisfied. The simultaneous injection of brine and ScCO 2 at 0.2 and 0.6 cc/min respectively at reservoir conditions (1500 psi and 35 • C) to maintain 75% injection foam quality behave as the other baseline. Then, CO 2 emulsion corefloodings are conducted similar to the above W&G co-injection, except for an additional 0.2 wt% surfactant in the brine. The pre-saturation of the surfactant solution is not employed to resemble the practical conditions in the field. Outcrop cores are not reused, and properties such as initial oil saturation are listed in Table 1. Pressure drops across the core and oil production with time are recorded.

Simulation Deck
As mentioned, foam behaviors are sensitive to in-situ phase velocities, which implies the radial flow in a cylindrical coordinator would be more appropriate for foam modeling, and better resemble the practical scenario. A 2D sector (theta 15 • , porosity 20%) simulation deck is built, as shown in Figure 3, which is homogeneous but anisotropic with respect to permeability, except for the outermost grid column. The permeabilities in radial and vertical directions of the main body are set as 200 and 400 md, respectively. In the outermost column, the vertical producer locates holds permeability as high as 10,000 Darcy in every direction to mimic an open boundary condition and prevent artificial backflow [70]. An injector is positioned in innermost column grids vertically. The sector is 440 ft in a radial direction, which is divided into 100 grids unevenly. The grid size is 3 ft for the first 30 grids, and increases to 5 ft for the remains. It is as thick as 100 ft in the Z direction, which is cut to 20 layers evenly. Initial oil saturation is 0.4 without the presence of a gaseous phase and any surfactants. Oil properties have been tuned with CMG/WINPROP to ensure the achievement of multiple contact miscible (MMP~1082 psi) under initial reservoir conditions (1500 psi, 35 • C). Measured partition coefficients between brine/ScCO 2 of selected CPS (2EH-PO 5 -EO 15 ) is fed into STARS directly. Sustain

Foam Model
Different modeling techniques and numerical models were proposed, classified as implicit texture models [71], pore network [72], fractional flow theory [73] and population balance equations [74]. Here, in CMG/STARS, an empirical/implicit texture foam model is employed to ascribe foam effects to gas relative permeability reduction, represented by a dimensionless factor, FM, which is the product of multiple factors, where fmmob denotes as reference mobility reduction factor and F1 to F7 stand for various effects from an application point of view. To meet the purpose of demonstrations and comparisons but also prevent the unnecessary complexity, only the effects of surfactant concentration (F1) and dry out (F7) are employed here, which are defined in the following, where Cs denotes as in-situ surfactant concentration; fmsurf is a critical concentration above which foam strength would not increase further; epsurf indicates the dependence of gas mobility on surfactant concentration; fmdry is a critical water saturation at which foam would coalescence significantly and epdry controls the shape of the curve when approaching fmdry. Here, fmsurf is set as the same with surfactant concentration in an injection aqueous stream, and a typical value of foam mobility reduction in the field [75] is chosen for fmmob (Table 2). To prevent unrealistic gas trapping, fmdry is slightly higher than irreducible water saturation (Table 2).  Table 3 summarizes the tested injection schemes which are under a constant rate constraint. Continuous pure ScCO2 injection is carried out as a baseline, in rates of 420,000 scf/d to maintain the same voidage displacement with the following alternating injections. WAG injection acts as another baseline, with rates at 45 bbl/d and 315,000 scf/d respec-

Foam Model
Different modeling techniques and numerical models were proposed, classified as implicit texture models [71], pore network [72], fractional flow theory [73] and population balance equations [74]. Here, in CMG/STARS, an empirical/implicit texture foam model is employed to ascribe foam effects to gas relative permeability reduction, represented by a dimensionless factor, FM, which is the product of multiple factors, where fmmob denotes as reference mobility reduction factor and F 1 to F 7 stand for various effects from an application point of view. To meet the purpose of demonstrations and comparisons but also prevent the unnecessary complexity, only the effects of surfactant concentration (F 1 ) and dry out (F 7 ) are employed here, which are defined in the following, where C s denotes as in-situ surfactant concentration; fmsurf is a critical concentration above which foam strength would not increase further; epsurf indicates the dependence of gas mobility on surfactant concentration; fmdry is a critical water saturation at which foam would coalescence significantly and epdry controls the shape of the curve when approaching fmdry. Here, fmsurf is set as the same with surfactant concentration in an injection aqueous stream, and a typical value of foam mobility reduction in the field [75] is chosen for fmmob (Table 2). To prevent unrealistic gas trapping, fmdry is slightly higher than irreducible water saturation ( Table 2).  Table 3 summarizes the tested injection schemes which are under a constant rate constraint. Continuous pure ScCO 2 injection is carried out as a baseline, in rates of 420,000 scf/d to maintain the same voidage displacement with the following alternating injections. WAG injection acts as another baseline, with rates at 45 bbl/d and 315,000 scf/d Sustain. Chem. 2021, 2 134 respectively, to achieve a 75% injection gas volume fraction at reservoir conditions, when maintaining the slug ratio of brine to ScCO 2 at 1:1 and the unit slug size at 180 days. With corresponding law, the compressibility factor of ScCO 2 is derived with its known critical pressure and temperature. Accordingly, CO 2 volumes could be converted between downhole and surface (14.7 psi, 20 • C). As mentioned above, 420,000 scf maintains a similar reservoir volume to the summation of 45 bbl and 315,000 scf. Meanwhile, the implementation time for alternating injection is doubled (2880 days) relative to the continuous injection to maintain the same amount of total fluids injection (brine and ScCO 2 ). Lastly, the novel strategy, CIDS, is studied with reduced surfactant concentration (molar fraction) in ScCO 2 phase to remain similar among surfactant injections. Here, the results at the end of injection will be evaluated through the observed gas production rate, well bottom hole pressure (BHP), gas saturation, surfactant distribution, pressure distribution, CO 2 storage (the difference between injection & production CO 2 ), CO 2 retention factor (CO 2 storage/CO 2 injection × 100%), oil production and CO 2 utilization ratio (oil production/CO 2 injection). Multiple factors would impose significant impacts on surfactant performances on the rock surface, including surface nature, fluid/rock interactions and fluid electrolyte composition/strength [76]. The HPLC analyses indicate that the static adsorption of CD-1045 and 2EH-PO 5 -EO 15 on Silurian dolomite are 1.1 and 0.25 mg/g, respectively. The more than four-times-higher affinity to dolomite rock of employed ASS than selected CPS could be ascribed to the intrinsic molecular structure differences. Ionic surfactant behaviors are mainly affected by the electrostatic attraction/repulsion between the charged head group and rock surface [12]. Dolomite surface is positively charged in neutral pH brine (~7.2), with an iso-potential point of around pH 8.2 [77], which results in strong attraction with negatively charged anionic CD-1045. On the contrary, nonionic surfactant adsorption on dolomite depends on the silica and clay content in the mineral, and electrostatic force is not the dominant factor [78]. Therefore, the adsorption of 2EH-PO 5 -EO 15 is mainly driven by hydrogen bonding between oxygen in the ethoxy groups (EO) of surfactant and the hydroxyl groups on the mineral surface, which is weaker than the electrostatic force. In addition, the relatively large number of EO groups in 2EH-PO 5 -EO 15 enhances its hydrophilicity, which in turn reduces its adsorption [65].

Phase Behavior Test
It is important to ensure the foam/CO 2 emulsion generation and stability are not adversely impacted by the preferential partition of employed surfactants to the W/O interface, which may result in surfactant loss and the generation of stable viscous W/O emulsion or microemulsion [79,80]. It is also crucial to prevent misinterpretations by foam alone if microemulsion [68] or macroemulsion [69] is indeed present. Figure 4 illustrates the phase behaviors observations of selected surfactant solutions with Wasson crude oil under variant salinities. Neither of these display the typical Winsor type phase behaviors. It is observed that CD-1045 produces small amount of light brown separate phase (Figure 4a), which is slightly more viscous than original oil. It is similar to unstable liquid crystal which may be caused by a certain high molecular component in oil, such as asphaltene. In general, the detrimental impacts to foam stability should not be severe, which is supported by the literatures that state CD-1045 is a proven good foaming agent [81]. Meanwhile, the novel CPS demonstrates negligible affinity to employed Wasson oil, with an almost indiscernible inter phase between oil and brine. Similar conclusions were drawn by Chen et al. [82], with another nonionic CPS whose partition coefficient between W/O was very small. Hence, the observations here eliminate the concerns of impairment to foam stability by surfactant loss or interference to observation interpretation by the potential presence of additional phases, which implies the poor affinities of both selected ASS and CPS to Wasson crude oil as foaming agents.
Sustain. Chem. 2021, 2, FOR PEER REVIEW 9 It is observed that CD-1045 produces small amount of light brown separate phase ( Figure  4a), which is slightly more viscous than original oil. It is similar to unstable liquid crystal which may be caused by a certain high molecular component in oil, such as asphaltene.
In general, the detrimental impacts to foam stability should not be severe, which is supported by the literatures that state CD-1045 is a proven good foaming agent [81]. Meanwhile, the novel CPS demonstrates negligible affinity to employed Wasson oil, with an almost indiscernible inter phase between oil and brine. Similar conclusions were drawn by Chen et al. [82], with another nonionic CPS whose partition coefficient between W/O was very small. Hence, the observations here eliminate the concerns of impairment to foam stability by surfactant loss or interference to observation interpretation by the potential presence of additional phases, which implies the poor affinities of both selected ASS and CPS to Wasson crude oil as foaming agents.

Two-Phase Flooding with Unfractured Berea Sandstone
Sandstone is normally negatively charged with neutral pH brine (~7.2) with an isopotential point of around pH 5 [83], which results in similar adsorption levels for both anionic and nonionic surfactants. Therefore, the comparisons between two types of CO2 emulsions would be ascribed to their foam behaviors only. Figure 5 displays the sectional pressure drops across the cores of ASS ( Figure 5a) and CPS (Figure 5b), which indicate not only the propagation, but also the apparent viscosity of CO2 emulsions. After satisfying minimum pressure gradient for foam generation [83], strong ASS CO2 emulsion appears at three total injection PV (TIPV, liquid & CO2) at Section 1, while a much shorter TIPV is required for CPS emulsion (~1 TIPV). This earlier response is mainly attributed to the dual partition of CPS, which has one of the most supreme superiorities over conventional ASS. Before the strong foam generation, ScCO2 mobility (~0.058 cp, 708 kg/m 3 and water-wet @ 1500 psi and 35 °C) was much higher than brine, resulting in the faster spreading of CPS [84], which could promote the apparent emulsion propagation. Further propagating downward, CPS emulsion reaches Sections 2 and 3 at 5 and 12 TIPV, respectively (Figure 5b), which are still earlier than responses from ASS emulsion at 8 and 15 TIPV, respectively. However, those similar subsequent observations should be more attributed to the better stabilizing ability of CPS to emulsion film rather than surfactant spreading, as addressed above. It has been claimed that foam would not affect the liquid relative permeability and propagation directly [85]. Therefore, either type of employed surfactant would already spread to the end of core after such amount of TIPV with current 75% injection foam quality. It is important to note that owing to its thermodynamically instability, foam actually "propagates" in the manner of breaking and reforming [86]. In other words, foam or CO2 emulsion could not stay intact as a homogenous and isotropic phase with unchanged properties to propagate, as conventional polymer solution does. Instead, it has to regenerate downstream after strong foam has been generated upstream. Hence, the observed earlier responses in Section 2 and 3 from CPS emulsion indeed imply its better stabilizing ability to film. It has been reported that the CPS employed here could reduce the ScCO2/water IFT to 5.6 mN/m (@ 2000 psi, 24 °C, 0.01 wt%) [87], relative to 9.5 mN/m, which CD-1045 could achieved under similar conditions [88]. Lower IFT would explicitly

Two-Phase Flooding with Unfractured Berea Sandstone
Sandstone is normally negatively charged with neutral pH brine (~7.2) with an isopotential point of around pH 5 [83], which results in similar adsorption levels for both anionic and nonionic surfactants. Therefore, the comparisons between two types of CO 2 emulsions would be ascribed to their foam behaviors only. Figure 5 displays the sectional pressure drops across the cores of ASS ( Figure 5a) and CPS (Figure 5b), which indicate not only the propagation, but also the apparent viscosity of CO 2 emulsions. After satisfying minimum pressure gradient for foam generation [83], strong ASS CO 2 emulsion appears at three total injection PV (TIPV, liquid & CO 2 ) at Section 1, while a much shorter TIPV is required for CPS emulsion (~1 TIPV). This earlier response is mainly attributed to the dual partition of CPS, which has one of the most supreme superiorities over conventional ASS. Before the strong foam generation, ScCO 2 mobility (~0.058 cp, 708 kg/m 3 and water-wet @ 1500 psi and 35 • C) was much higher than brine, resulting in the faster spreading of CPS [84], which could promote the apparent emulsion propagation. Further propagating downward, CPS emulsion reaches Sections 2 and 3 at 5 and 12 TIPV, respectively (Figure 5b), which are still earlier than responses from ASS emulsion at 8 and 15 TIPV, respectively. However, those similar subsequent observations should be more attributed to the better stabilizing ability of CPS to emulsion film rather than surfactant spreading, as addressed above. It has been claimed that foam would not affect the liquid relative permeability and propagation directly [85]. Therefore, either type of employed surfactant would already spread to the end of core after such amount of TIPV with current 75% injection foam quality. It is important to note that owing to its thermodynamically instability, foam actually "propagates" in the manner of breaking and reforming [86]. In other words, foam or CO 2 emulsion could not stay intact as a homogenous and isotropic phase with unchanged properties to propagate, as conventional polymer solution does. Instead, it has to regenerate downstream after strong foam has been generated upstream. Hence, the observed earlier responses in Sections 2 and 3 from CPS emulsion indeed imply its better stabilizing ability to film. It has been reported that the CPS employed here could reduce the ScCO 2 /water IFT to 5.6 mN/m (@ 2000 psi, 24 • C, 0.01 wt%) [87], relative to 9.5 mN/m, which CD-1045 could achieved under similar conditions [88]. Lower IFT would explicitly reduce the difference of capillary pressures between the plateau border and film, which in turn weakens the capillary suction that is the main driving mechanism of film coalescence [89]. This superior film stabilization ability of CPS significantly reduces the time needed in every section approaching the plateau values (Figure 5a,b). In addition, prominently higher apparent emulsion strength, indicated by more than doubled total pressure drops (Figure 5a,b), could also be attributed to the more robust emulsion film stabilized by CPS. In turn, the combinations of more rapid propagation, faster foam strength buildup and higher apparent emulsion viscosity result in a much steeper water desaturation rate and lower eventual water saturation (Figure 5c, CPS 0.2 VS ASS 0.46), which clearly manifests the superiorities of CPS ScCO 2 emulsion over the conventional one in fluid displacement efficiency and CO 2 storage.
Sustain. Chem. 2021, 2, FOR PEER REVIEW 10 reduce the difference of capillary pressures between the plateau border and film, which in turn weakens the capillary suction that is the main driving mechanism of film coalescence [89]. This superior film stabilization ability of CPS significantly reduces the time needed in every section approaching the plateau values (Figure 5a,b). In addition, prominently higher apparent emulsion strength, indicated by more than doubled total pressure drops (Figure 5a,b), could also be attributed to the more robust emulsion film stabilized by CPS. In turn, the combinations of more rapid propagation, faster foam strength buildup and higher apparent emulsion viscosity result in a much steeper water desaturation rate and lower eventual water saturation (Figure 5c, CPS 0.2 VS ASS 0.46), which clearly manifests the superiorities of CPS ScCO2 emulsion over the conventional one in fluid displacement efficiency and CO2 storage.

Two-Phase Flooding with Unfractured Silurian Dolomite
When replacing consolidate porous media with Silurian dolomite, similar relative behaviors are observed (Figure 6), except for the Section 1, which may be caused by the lithology of employed dolomite demonstrating a stronger inlet effect [90]. It takes around 2 TIPV for ASS emulsion to demonstrate a discernible pressure response in Section 2, while only as less as 0.5 TIPV is needed for CPS emulsion. Furthermore, once strong emulsion has been generated, it takes much less TIPV for CPS to buildup in each section and propagate downstream (Figure 6a,b, 1.3 TIPV VS 5 TIPV). In particular, with Silurian dolomite, ASS emulsion could hardly reach comparable strength in both sections, which partially results in lower total emulsion strength relative to CPS emulsion. Accordingly, much clearer contrast is displayed on the average water saturation in the cores (Figure 6c). The desaturation rate of CPS emulsion is very high, and CO2 storage reaches more than 90%

Two-Phase Flooding with Unfractured Silurian Dolomite
When replacing consolidate porous media with Silurian dolomite, similar relative behaviors are observed (Figure 6), except for the Section 1, which may be caused by the lithology of employed dolomite demonstrating a stronger inlet effect [90]. It takes around 2 TIPV for ASS emulsion to demonstrate a discernible pressure response in Section 2, while only as less as 0.5 TIPV is needed for CPS emulsion. Furthermore, once strong emulsion has been generated, it takes much less TIPV for CPS to buildup in each section and propagate downstream (Figure 6a,b, 1.3 TIPV VS 5 TIPV). In particular, with Silurian dolomite, ASS emulsion could hardly reach comparable strength in both sections, which partially results in lower total emulsion strength relative to CPS emulsion. Accordingly, much clearer contrast is displayed on the average water saturation in the cores (Figure 6c). The desaturation rate of CPS emulsion is very high, and CO 2 storage reaches more than 90% in 2.8 TIPV, while much gentler stepwise behaviors are shown by ASS emulsion. Moreover, the eventual water saturation is reduced to as low as 0.03 by CPS emulsion which is five times less than the conventional one. The observations here could be ascribed to significantly reduced CPS adsorption on the dolomite rock surface, in addition to the superiorities of CPS emulsion as observed on Berea sandstone.
Sustain. Chem. 2021, 2, FOR PEER REVIEW 11 in 2.8 TIPV, while much gentler stepwise behaviors are shown by ASS emulsion. Moreover, the eventual water saturation is reduced to as low as 0.03 by CPS emulsion which is five times less than the conventional one. The observations here could be ascribed to significantly reduced CPS adsorption on the dolomite rock surface, in addition to the superiorities of CPS emulsion as observed on Berea sandstone.

Three-Phase CO2 & W/G Floodings with Fractured Silurian Dolomite
The permeability of an open space between parallel plates-that is, the configuration of a fracture in current composition system-is proportional to the square of a fracture aperture. Thereby, 75 micron corresponds to 475 Darcy, which increases the composite core permeability to around 1 Darcy with 150 mD matrix permeability. Thus, the fracture behaves as a severe theft zone, with tremendously high conductivity. In turn, typical pure ScCO2 flooding produces very low oil recovery (24%), corresponding to negligible pressure drop (Figure 7), even though MMP is currently achieved. In such a short core without mobility control, with the presence of a high conductivity fracture along central axis, advective flow in longitudinal direction is dominant over transverse flow in radial direction, not to mention the molecular diffusion, where multiple contact miscibility would play a very limited role [91]. In contrast, the co-injection of brine and ScCO2 slightly improves the performances, but to a limited extent. Injection fluids still break through early at 0.3 TIPV, and level off at 35% oil recovery, accompanied with pressure drops in a similar order of magnitude (Figure 7). The effects of additional brine phase injection are two-fold. On one hand, it was suggested that the simultaneous injection could provide better mobility control than continuous single-phase injection (brine of gas) to reduce oil capillary entrapment via the multiple phase relative permeability effect [92]; nevertheless, the en-

Three-Phase CO 2 & W/G Floodings with Fractured Silurian Dolomite
The permeability of an open space between parallel plates-that is, the configuration of a fracture in current composition system-is proportional to the square of a fracture aperture. Thereby, 75 micron corresponds to 475 Darcy, which increases the composite core permeability to around 1 Darcy with 150 mD matrix permeability. Thus, the fracture behaves as a severe theft zone, with tremendously high conductivity. In turn, typical pure ScCO 2 flooding produces very low oil recovery (24%), corresponding to negligible pressure drop (Figure 7), even though MMP is currently achieved. In such a short core without mobility control, with the presence of a high conductivity fracture along central axis, advective flow in longitudinal direction is dominant over transverse flow in radial direction, not to mention the molecular diffusion, where multiple contact miscibility would play a very limited role [91]. In contrast, the co-injection of brine and ScCO 2 slightly improves the performances, but to a limited extent. Injection fluids still break through early at 0.3 TIPV, and level off at 35% oil recovery, accompanied with pressure drops in a similar order of magnitude (Figure 7). The effects of additional brine phase injection are two-fold. On one hand, it was suggested that the simultaneous injection could provide better mobility control than continuous single-phase injection (brine of gas) to reduce oil capillary entrapment via the multiple phase relative permeability effect [92]; nevertheless, the enhancement may be more valid for small-scale heterogeneity. On the other hand, a significant amount of brine may reduce CO 2 -oil contact for a miscible process through the water shielding effect [31]. Therefore, it seems like the current scenario is still dominated by the lack of sufficient mobility control, and only the regions near fracture are swept. hancement may be more valid for small-scale heterogeneity. On the other hand, a sign cant amount of brine may reduce CO2-oil contact for a miscible process through the w shielding effect [31]. Therefore, it seems like the current scenario is still dominated by lack of sufficient mobility control, and only the regions near fracture are swept.

Three-Phase ScCO2 Emulsion Floodings with Fractured Silurian Dolomite
Relative to the above baselines, the performance improvements with the presenc ASS are prominent. As shown in Figure 8, even though a similar 30% oil is recover indicating the sweeping of near fracture regions, further continuous oil production is served, implying the improvement of mobility control and desaturation in a radial di tion, which results in an eventual 54% oil recovery with more than a forty-times incre of the pressure drop. The enhancement of liquid production also suggests the mar improvement of CO2 storage. As mentioned, the presence of a surfactant could stabi the thin film, called lamella, which separates the dispersing bubbles. Due to this, disp ing phase (gas or ScCO2) mobility could be reduced via increased apparent viscosity, ing to shear stresses between porous media and film [93] and through reduced effec relative permeability due to dispersing phase trapping [85]. It has been proven that C 1045 was an effective CO2 foaming agent under variable reservoir conditions [81]. In c trast, when replacing stabilizer agent with CPS, more than 10 times the pressure d across the composite core is observed, and oil recovery is further enhanced to 70%, wh nearly reaches the levels in unfractured systems. The high conductivity fracture co further promote the CPS spreading rate when it transports with high-mobility ScCO2. superior performances here of CPS emulsion in such harsh reservoir conditions could ascribed to not only the better film stabilization ability, as mentioned above, but also extraordinary tolerance to the oleic phase, as illustrated in prior test, since oil saturat in a fracture is much higher than that of a matrix, when it behaves as a flow path.

Three-Phase ScCO 2 Emulsion Floodings with Fractured Silurian Dolomite
Relative to the above baselines, the performance improvements with the presence of ASS are prominent. As shown in Figure 8, even though a similar 30% oil is recovered, indicating the sweeping of near fracture regions, further continuous oil production is observed, implying the improvement of mobility control and desaturation in a radial direction, which results in an eventual 54% oil recovery with more than a forty-times increase of the pressure drop. The enhancement of liquid production also suggests the marked improvement of CO 2 storage. As mentioned, the presence of a surfactant could stabilize the thin film, called lamella, which separates the dispersing bubbles. Due to this, dispersing phase (gas or ScCO 2 ) mobility could be reduced via increased apparent viscosity, owing to shear stresses between porous media and film [93] and through reduced effective relative permeability due to dispersing phase trapping [85]. It has been proven that CD-1045 was an effective CO 2 foaming agent under variable reservoir conditions [81]. In contrast, when replacing stabilizer agent with CPS, more than 10 times the pressure drop across the composite core is observed, and oil recovery is further enhanced to 70%, which nearly reaches the levels in unfractured systems. The high conductivity fracture could further promote the CPS spreading rate when it transports with high-mobility ScCO 2 . The superior performances here of CPS emulsion in such harsh reservoir conditions could be ascribed to not only the better film stabilization ability, as mentioned above, but also the extraordinary tolerance to the oleic phase, as illustrated in prior test, since oil saturation in a fracture is much higher than that of a matrix, when it behaves as a flow path.
shielding effect [31]. Therefore, it seems like the current scenario is still dominated by lack of sufficient mobility control, and only the regions near fracture are swept.

Three-Phase ScCO2 Emulsion Floodings with Fractured Silurian Dolomite
Relative to the above baselines, the performance improvements with the presen ASS are prominent. As shown in Figure 8, even though a similar 30% oil is recove indicating the sweeping of near fracture regions, further continuous oil production i served, implying the improvement of mobility control and desaturation in a radial d tion, which results in an eventual 54% oil recovery with more than a forty-times incr of the pressure drop. The enhancement of liquid production also suggests the ma improvement of CO2 storage. As mentioned, the presence of a surfactant could stab the thin film, called lamella, which separates the dispersing bubbles. Due to this, disp ing phase (gas or ScCO2) mobility could be reduced via increased apparent viscosity, ing to shear stresses between porous media and film [93] and through reduced effe relative permeability due to dispersing phase trapping [85]. It has been proven that 1045 was an effective CO2 foaming agent under variable reservoir conditions [81]. In trast, when replacing stabilizer agent with CPS, more than 10 times the pressure across the composite core is observed, and oil recovery is further enhanced to 70%, w nearly reaches the levels in unfractured systems. The high conductivity fracture c further promote the CPS spreading rate when it transports with high-mobility ScCO2 superior performances here of CPS emulsion in such harsh reservoir conditions coul ascribed to not only the better film stabilization ability, as mentioned above, but also extraordinary tolerance to the oleic phase, as illustrated in prior test, since oil satura in a fracture is much higher than that of a matrix, when it behaves as a flow path.

Simulations
Five injection scenarios are tested to illustrate the necessary of mobility control and superiors of CPS emulsion in practical reservoir conditions with the presence of a gravity field. Figure 9 demonstrates the gas production rates of variable scenarios, in which the ScCO 2 times are highlighted (upper right sub-graph). ScCO 2 BTs rapidly in 55 days for pure ScCO 2 injection, and approaches the plateau quickly, which indicates the necessity of mobility control. With the presence of a WAG injection, it indeed delays the gas BT until 240 days. Densities of brine and ScCO 2 are conspicuously distinct, which could result in three flow zones in the reservoir, owing to the countercurrent flow of water and gas [94], i.e., mixing zone (near wellbore comprising both water and gas), override zone (mainly consisting of gas) and underride zone (mainly consisting of water). Therefore, the downwards flow of brine impedes the migration of ScCO 2 upwards to some extent. When adding ASS in water stream, ScCO 2 BT is further delayed to 282 days with a prominently reduced gas production rate, which is caused by the obviously reduced gas mobility than WAG [95,96]. With novel CPS, even better performances are observed that ScCO 2 BTs until 295 days with another~10% gas production rate decrease. Owing to its dual phase partition capacity, CPS would be carried by injected ScCO 2 to the traditional override zone, and better mitigate the issue of gravity segregation. When applying the most novel strategy, CIDS, relative to pure ScCO 2 injection, it did reduce the gas production rate by 25% after ScCO 2 BT, even though the BT time is just slightly delayed to 60 days. CIDS could generate ScCO 2 emulsion with in-situ brine, and further analysis would give more details in the following.

Simulations
Five injection scenarios are tested to illustrate the necessary of mobility con superiors of CPS emulsion in practical reservoir conditions with the presence of a field. Figure 9 demonstrates the gas production rates of variable scenarios, in w ScCO2 times are highlighted (upper right sub-graph). ScCO2 BTs rapidly in 55 pure ScCO2 injection, and approaches the plateau quickly, which indicates the n of mobility control. With the presence of a WAG injection, it indeed delays th until 240 days. Densities of brine and ScCO2 are conspicuously distinct, which cou in three flow zones in the reservoir, owing to the countercurrent flow of water [94], i.e., mixing zone (near wellbore comprising both water and gas), overr (mainly consisting of gas) and underride zone (mainly consisting of water). There downwards flow of brine impedes the migration of ScCO2 upwards to some exten adding ASS in water stream, ScCO2 BT is further delayed to 282 days with a prom reduced gas production rate, which is caused by the obviously reduced gas mobi WAG [95,96]. With novel CPS, even better performances are observed that ScCO til 295 days with another ~10% gas production rate decrease. Owing to its dual ph tition capacity, CPS would be carried by injected ScCO2 to the traditional overri and better mitigate the issue of gravity segregation. When applying the most no egy, CIDS, relative to pure ScCO2 injection, it did reduce the gas production rate after ScCO2 BT, even though the BT time is just slightly delayed to 60 days. CID generate ScCO2 emulsion with in-situ brine, and further analysis would give mor in the following. More direct illustrations are shown in Figure 10 of gas saturation distributio end of injections. As aforementioned, pure ScCO2 injection (Figure 10a) and ev flooding (Figure 10b) cause a lack of appropriate mobility control to the gas ph severe gravity segregation would promote low viscosity and cause the density C to mitigate to the override zone at a few top layers, and leave significant upswe by miscible flooding; the presence of a surfactant, especially CPS (Figure 10d), m those adverse impacts remarkably. In particular, CPS emulsion not only extend zone pronouncedly deeper into reservoir, as displayed by ASS emulsion (Figure also pulls the override zone downwards, which is ascribed to the better gas mob trol in the override zone. More striking effectiveness is provided by CIDS injectio no area is unswept by ScCO2 (Figure 10e), which is expected to give excellent C More direct illustrations are shown in Figure 10 of gas saturation distributions at the end of injections. As aforementioned, pure ScCO 2 injection ( Figure 10a) and even WAG flooding (Figure 10b) cause a lack of appropriate mobility control to the gas phase, and severe gravity segregation would promote low viscosity and cause the density CO 2 phase to mitigate to the override zone at a few top layers, and leave significant upswept areas by miscible flooding; the presence of a surfactant, especially CPS (Figure 10d), mitigates those adverse impacts remarkably. In particular, CPS emulsion not only extends mixing zone pronouncedly deeper into reservoir, as displayed by ASS emulsion (Figure 10c), but also pulls the override zone downwards, which is ascribed to the better gas mobility control in the override zone. More striking effectiveness is provided by CIDS injection where no area is unswept by ScCO 2 (Figure 10e), which is expected to give excellent CO 2 final retention. As shown in Figure 11a, CPS emulsion stores more than 63% injected CO 2 than emulsion using ASS, while another 28% improvement is achieved by novel CIDS injection. However, it is noted that cumulative CO 2 injections are different among tested scenarios, thereby the CO 2 retention factor (CO 2 storage/injection) could have a better quantity for storage efficiency. Figure 11b indicates that alternating injection with CPS provides the supreme CO 2 retention efficiency, which is 65% higher than conventional ASS emulsion, while CIDS performs comparably, as just 4% lower.
Sustain. Chem. 2021, 2, FOR PEER REVIEW 14 retention. As shown in Figure 11a, CPS emulsion stores more than 63% injected CO2 than emulsion using ASS, while another 28% improvement is achieved by novel CIDS injection. However, it is noted that cumulative CO2 injections are different among tested scenarios, thereby the CO2 retention factor (CO2 storage/injection) could have a better quantity for storage efficiency. Figure 11b indicates that alternating injection with CPS provides the supreme CO2 retention efficiency, which is 65% higher than conventional ASS emulsion, while CIDS performs comparably, as just 4% lower. (a) (b) Figure 11. Gas retention analysis among tested injection scenarios, (a) CO2 storage; and (b) CO2 retention factor.
The surfactant concentration distributions ( Figure 12) provide better interpretations of the above observations, and reveal the essential causes of CPS superiorities. Owing to water slumping and only the aqueous solubility of ASS, a thin topmost layer behaves as a surfactant vacuum zone lack of CO2 mobility control, and most ASS concentrates near wellbore. On the contrary, CPS with alternating injection could be carried by injected Sustain. Chem. 2021, 2, FOR PEER REVIEW 14 retention. As shown in Figure 11a, CPS emulsion stores more than 63% injected CO2 than emulsion using ASS, while another 28% improvement is achieved by novel CIDS injection. However, it is noted that cumulative CO2 injections are different among tested scenarios, thereby the CO2 retention factor (CO2 storage/injection) could have a better quantity for storage efficiency. Figure 11b indicates that alternating injection with CPS provides the supreme CO2 retention efficiency, which is 65% higher than conventional ASS emulsion, while CIDS performs comparably, as just 4% lower. (a) (b) Figure 11. Gas retention analysis among tested injection scenarios, (a) CO2 storage; and (b) CO2 retention factor.
The surfactant concentration distributions ( Figure 12) provide better interpretations of the above observations, and reveal the essential causes of CPS superiorities. Owing to water slumping and only the aqueous solubility of ASS, a thin topmost layer behaves as a surfactant vacuum zone lack of CO2 mobility control, and most ASS concentrates near wellbore. On the contrary, CPS with alternating injection could be carried by injected Figure 11. Gas retention analysis among tested injection scenarios, (a) CO 2 storage; and (b) CO 2 retention factor.
The surfactant concentration distributions ( Figure 12) provide better interpretations of the above observations, and reveal the essential causes of CPS superiorities. Owing to water slumping and only the aqueous solubility of ASS, a thin topmost layer behaves as a surfactant vacuum zone lack of CO 2 mobility control, and most ASS concentrates near wellbore. On the contrary, CPS with alternating injection could be carried by injected ScCO 2 and be present in those top layers to slow down the gas mobility and facilitate the CO 2 diversion downwards. A much more prominent contrast is shown by CIDS, which displays almost uniform high concentration in the whole swept region. Moreover, it is interesting to note that in both scenarios employing CPS, the surfactant vacuum zones instead appear near the wellbore vertically, which could be attributed to the liquid dry out effect and the dual phase partition of CPS. The aqueous phase has marginal solubility in CO 2 phase but is nonnegligible under the scenario with a significant amount of CO 2 injection. Here, water saturation near wellbore would be reduced markedly, resulting in a noteworthy increase of surfactant concentration in the aqueous phase. Due to this, the preferential partition of CPS into an injected ScCO 2 stream will be driven by thermodynamic equilibrium under the constant partition coefficient. Consequently, foam would be weakened, caused by a decrease of local surfactant concentration. Meanwhile, conspicuously reduced water saturation may invoke the so-called foam dry out effect, resulting in the further impairment of foam strength. Here, the dynamic partition, in conjunction with water dry-out, promotes the appearance of the CPS vacuum zone near the wellbore, which potentially improves the injectivity, as shown below. The superiorities of CPS emulsion are also manifested by the economic analysis for oil recovery and CO 2 utilization efficiency (oil production/CO 2 injection), as shown in Figure 13. Alternating injections of CPS emulsion doubles the oil production over WAG, and produces 30% more than the scenario using ASS (Figure 13a). More strikingly, CIDS injection could provide another 13% additional oil recovery. The variation of CO 2 utilization ratio is consistent with oil production and monotonically increasing among tested scenarios, which corresponds with the increase of industrial revenue (Figure 13b). In addition, without liquid phase injection, the operation cost will be further reduced to benefit the balance sheet.
Sustain. Chem. 2021, 2, FOR PEER REVIEW 15 ScCO2 and be present in those top layers to slow down the gas mobility and facilitate the CO2 diversion downwards. A much more prominent contrast is shown by CIDS, which displays almost uniform high concentration in the whole swept region. Moreover, it is interesting to note that in both scenarios employing CPS, the surfactant vacuum zones instead appear near the wellbore vertically, which could be attributed to the liquid dry out effect and the dual phase partition of CPS. The aqueous phase has marginal solubility in CO2 phase but is nonnegligible under the scenario with a significant amount of CO2 injection. Here, water saturation near wellbore would be reduced markedly, resulting in a noteworthy increase of surfactant concentration in the aqueous phase. Due to this, the preferential partition of CPS into an injected ScCO2 stream will be driven by thermodynamic equilibrium under the constant partition coefficient. Consequently, foam would be weakened, caused by a decrease of local surfactant concentration. Meanwhile, conspicuously reduced water saturation may invoke the so-called foam dry out effect, resulting in the further impairment of foam strength. Here, the dynamic partition, in conjunction with water dry-out, promotes the appearance of the CPS vacuum zone near the wellbore, which potentially improves the injectivity, as shown below. The superiorities of CPS emulsion are also manifested by the economic analysis for oil recovery and CO2 utilization efficiency (oil production/CO2 injection), as shown in Figure 13. Alternating injections of CPS emulsion doubles the oil production over WAG, and produces 30% more than the scenario using ASS (Figure 13a). More strikingly, CIDS injection could provide another 13% additional oil recovery. The variation of CO2 utilization ratio is consistent with oil production and monotonically increasing among tested scenarios, which corresponds with the increase of industrial revenue (Figure 13b). In addition, without liquid phase injection, the operation cost will be further reduced to benefit the balance sheet.  Sustain. Chem. 2021, 2, FOR PEER REVIEW 15 ScCO2 and be present in those top layers to slow down the gas mobility and facilitate the CO2 diversion downwards. A much more prominent contrast is shown by CIDS, which displays almost uniform high concentration in the whole swept region. Moreover, it is interesting to note that in both scenarios employing CPS, the surfactant vacuum zones instead appear near the wellbore vertically, which could be attributed to the liquid dry out effect and the dual phase partition of CPS. The aqueous phase has marginal solubility in CO2 phase but is nonnegligible under the scenario with a significant amount of CO2 injection. Here, water saturation near wellbore would be reduced markedly, resulting in a noteworthy increase of surfactant concentration in the aqueous phase. Due to this, the preferential partition of CPS into an injected ScCO2 stream will be driven by thermodynamic equilibrium under the constant partition coefficient. Consequently, foam would be weakened, caused by a decrease of local surfactant concentration. Meanwhile, conspicuously reduced water saturation may invoke the so-called foam dry out effect, resulting in the further impairment of foam strength. Here, the dynamic partition, in conjunction with water dry-out, promotes the appearance of the CPS vacuum zone near the wellbore, which potentially improves the injectivity, as shown below. The superiorities of CPS emulsion are also manifested by the economic analysis for oil recovery and CO2 utilization efficiency (oil production/CO2 injection), as shown in Figure 13. Alternating injections of CPS emulsion doubles the oil production over WAG, and produces 30% more than the scenario using ASS (Figure 13a). More strikingly, CIDS injection could provide another 13% additional oil recovery. The variation of CO2 utilization ratio is consistent with oil production and monotonically increasing among tested scenarios, which corresponds with the increase of industrial revenue (Figure 13b). In addition, without liquid phase injection, the operation cost will be further reduced to benefit the balance sheet.  Injectivity is always an important concern during the field implementation of foam. Figure 14 displays the pressure distributions and contour lines at end of injection time frames, which provides some clues. It is interesting to note that with the variation of scenarios, not only the magnitude of pressures but also the pressure distributions, including the contour lines, are varying. Owing to the lack of mobility control, pure ScoCO 2 and WAG floodings show almost horizontal pressure contour lines, and high pressures distribute at the bottom (Figure 14a,b), corresponding to the phase separations in the reservoir. When applying the surfactants, pressure contour lines become more vertical-like morpha, indicating piston-like displacement and improvement of mobility control, but with different configurations. The pressure distribution of ASS emulsion is very uneven (Figure 14c), where high pressure zones compact in a few feet, and dissipate quickly away from wellbore. On the contrary, the scenarios using CPS demonstrate much more even pressure gradients with lower injection pressures (Figure 14d,e), corresponding to the surfactant vacuum zones near wellbore, as addressed above. The tilt extent of contour lines and pressure gradient distribution imply the energy utilization efficiency, and are relevant to the project cost. Furthermore, the maximum readings from above pressure distributions (Figure 15a) indicate the well injectivities, and correspond to the bottom hole pressures under the constant rate injection mode (Figure 15b). For alternating injection, CPS emulsion could reduce the max BHP as much as 250 psi, or equivalent to as high as a 76% injectivity increase, if using 1500 psi as average reservoir pressure in the far end. Moreover, when considering CIDS, one may have to compromise the remarkable injectivity reduction with improved oil production and CO 2 utilization ratio, relative to alternating injection, as addressed above.
Sustain. Chem. 2021, 2, FOR PEER REVIEW 16 Injectivity is always an important concern during the field implementation of foam. Figure 14 displays the pressure distributions and contour lines at end of injection time frames, which provides some clues. It is interesting to note that with the variation of scenarios, not only the magnitude of pressures but also the pressure distributions, including the contour lines, are varying. Owing to the lack of mobility control, pure ScoCO2 and WAG floodings show almost horizontal pressure contour lines, and high pressures distribute at the bottom (Figure 14a,b), corresponding to the phase separations in the reservoir. When applying the surfactants, pressure contour lines become more vertical-like morpha, indicating piston-like displacement and improvement of mobility control, but with different configurations. The pressure distribution of ASS emulsion is very uneven (Figure 14c), where high pressure zones compact in a few feet, and dissipate quickly away from wellbore. On the contrary, the scenarios using CPS demonstrate much more even pressure gradients with lower injection pressures (Figure 14d,e), corresponding to the surfactant vacuum zones near wellbore, as addressed above. The tilt extent of contour lines and pressure gradient distribution imply the energy utilization efficiency, and are relevant to the project cost. Furthermore, the maximum readings from above pressure distributions (Figure 15a) indicate the well injectivities, and correspond to the bottom hole pressures under the constant rate injection mode (Figure 15b). For alternating injection, CPS emulsion could reduce the max BHP as much as 250 psi, or equivalent to as high as a 76% injectivity increase, if using 1500 psi as average reservoir pressure in the far end. Moreover, when considering CIDS, one may have to compromise the remarkable injectivity reduction with improved oil production and CO2 utilization ratio, relative to alternating injection, as addressed above.

Conclusions
Here, systematic comparison studies have been carried out between aqueous and CO2 soluble surfactants as CO2 emulsion stabilizers. Laboratory experiments are conducted in homogeneous and fractured outcrops under both two-and three-phase conditions facilitated by auxiliary tests, which indicate the superiorities of this CO2-philic agent. Field scale simulations further predict its ascendancies over traditional aqueous soluble agents in practical conditions. The findings here would benefit both greenhouse control and storage, as well as enhancing oil recovery, which are summarized as the following:

•
The employed anionic aqueous soluble surfactant displays more than four times higher adsorption on Silurian dolomite than selected CO2-philic one, which is ascribed to stronger electrostatic attraction owing to opposite charges. The behaviors of nonionic CPS are dominated by hydrogen bonding, which is weaker than the electrostatic force.

•
Both ASS (CD-1045) and CPS (2EH-PO5-EO15) show poor affinity to Wasson oil, particularly the latter. In brief, the observations here eliminate the concern of surfactant loss to W/O interface, and both surfactants are mainly used for foaming purposes.

•
In the two-phase unfractured Berea sandstone system, CPS emulsion propagates much faster, which is attributed to its dual phase partition capacity before strong foam generation upstream, but to its superior film stabilization ability when it further propagates downstream. The better emulsion stability is also supported by the much faster pressure buildup and higher plateau pressure drops, which in turn facilitate CPS emulsion to achieve higher fluid displacement efficiency and CO2 storage. • Similar relative behaviors are observed in the two-phase unfractured Silurian dolomite system, with even sharper contrast in the displacement rate and CO2 storage. CPS emulsion could reduce brine saturation to less than 0.1 in 2.8 TIPV, and CO2 saturation could be as high as 97% eventually.

•
In the three-phase system with fractured Silurian dolomite, either continuous ScCO2 injection or co-injection of brine and ScCO2 provides relatively poor performances with negligible pressure drops, and oil recoveries are 24% and 35% respectively, even though MMP is achieved. It seems like here it is dominated by the lack of sufficient mobility control, and only the regions near the fracture are swept.

•
The selected ASS (CD-1045) indeed improves oil recovery to 54%, with more than a forty times increase of the pressure drop. When replacing the stabilizer agent with CPS, it is further enhanced to 70%, with another over 10 time pressure drop increase, which also indicates the far better CO2 storage.

Conclusions
Here, systematic comparison studies have been carried out between aqueous and CO 2 soluble surfactants as CO 2 emulsion stabilizers. Laboratory experiments are conducted in homogeneous and fractured outcrops under both two-and three-phase conditions facilitated by auxiliary tests, which indicate the superiorities of this CO 2 -philic agent. Field scale simulations further predict its ascendancies over traditional aqueous soluble agents in practical conditions. The findings here would benefit both greenhouse control and storage, as well as enhancing oil recovery, which are summarized as the following:

•
The employed anionic aqueous soluble surfactant displays more than four times higher adsorption on Silurian dolomite than selected CO 2 -philic one, which is ascribed to stronger electrostatic attraction owing to opposite charges. The behaviors of nonionic CPS are dominated by hydrogen bonding, which is weaker than the electrostatic force.

•
Both ASS (CD-1045) and CPS (2EH-PO 5 -EO 15 ) show poor affinity to Wasson oil, particularly the latter. In brief, the observations here eliminate the concern of surfactant loss to W/O interface, and both surfactants are mainly used for foaming purposes.

•
In the two-phase unfractured Berea sandstone system, CPS emulsion propagates much faster, which is attributed to its dual phase partition capacity before strong foam generation upstream, but to its superior film stabilization ability when it further propagates downstream. The better emulsion stability is also supported by the much faster pressure buildup and higher plateau pressure drops, which in turn facilitate CPS emulsion to achieve higher fluid displacement efficiency and CO 2 storage. • Similar relative behaviors are observed in the two-phase unfractured Silurian dolomite system, with even sharper contrast in the displacement rate and CO 2 storage. CPS emulsion could reduce brine saturation to less than 0.1 in 2.8 TIPV, and CO 2 saturation could be as high as 97% eventually.

•
In the three-phase system with fractured Silurian dolomite, either continuous ScCO 2 injection or co-injection of brine and ScCO 2 provides relatively poor performances with negligible pressure drops, and oil recoveries are 24% and 35% respectively, even though MMP is achieved. It seems like here it is dominated by the lack of sufficient mobility control, and only the regions near the fracture are swept.

•
The selected ASS (CD-1045) indeed improves oil recovery to 54%, with more than a forty times increase of the pressure drop. When replacing the stabilizer agent with CPS, it is further enhanced to 70%, with another over 10 time pressure drop increase, which also indicates the far better CO 2 storage.

Simulations
• Relative to ScoCO 2 and WAG injections, it is necessary to employ a surfactant for proper mobility control. CPS emulsion could further slow down the gas BT and reduce the CO 2 production rate, relative to ASS emulsion. In particular, CPS emulsion not only extends the mixing zone pronouncedly deeper into the reservoir, but also pulls the override zone downwards to provide better conformance control there. CIDS injection could even sweep the entire reservoir area. Accordingly, CPS emulsion stores more than 63% CO 2 than ASS emulsion, while another 28% improvement is achieved by novel CIDS injection. Meanwhile, CPS emulsion is also 65% more efficient on CO 2 retention than a conventional one. • Conventional ASS emulsion would leave an unswept surfactant vacuum zone on the topmost layer, owing to water slumping, while the dual phase partition capacity of CPS endues its possibility to present in those layers to realize the mobility control and facilitate the CO 2 diversion downwards. In turn, the alternating injection of CPS emulsion doubles the oil production than WAG, and produces 30% more than ASS emulsion, while CIDS injection could provide another 13% more oil recovery. The variation of the CO 2 utilization ratio is consistent with oil production among tested cases.

•
The pressure distribution of ASS emulsion is uneven and the pressure gradient is sharp near the wellbore, which implies low energy utilization efficiency. On the contrary, CPS emulsion provides a much more even pressure gradient. In addition, there exists a CPS vacuum zone near the wellbore, which is caused by the combination of ScCO 2 extraction to water, the dynamic partition of CPS into injected ScCO 2 stream and the foam dry-out effect. Thereby, CPS emulsion demonstrates significantly lower BHP and much higher injectivity. Critical water saturation for foam collapsing fmmob Reference mobility reduction factor fmsurf Critical surfactant concentration for "strong" foam (mole fraction) F 1 Effect of surfactant concentration for foam F 7 Effect of foam dry out by waster saturation FM Foam mobility reduction factor K G/L partition coefficient