Ruina Xu

Two‐phase flow properties of a sandstone rock for the CO2/water system: Core‐flooding experiments, and focus on impacts of mineralogical changes

The two-phase flow characterization (CO2/water) of a Triassic sandstone core from the Paris Basin, France, is reported in this paper. Absolute properties (porosity and water permeability), capillary pressure, relative permeability with hysteresis between drainage and imbibition, and residual trapping capacities have been assessed at 9 MPa pore pressure and 28°C (CO2 in liquid state) using a single core-flooding apparatus associated with magnetic resonance imaging. Different methodologies have been followed to obtain a data set of flow properties to be upscaled and used in large-scale CO2 geological storage evolution modeling tools. The measurements are consistent with the properties of well-sorted water-wet porous systems. As the mineralogical investigations showed a nonnegligible proportion of carbonates in the core, the experimental protocol was designed to observe potential impacts on flow properties of mineralogical changes. The magnetic resonance scanning and mineralogical observations indicate mineral dissolution during the experimental campaign, and the core-flooding results show an increase in porosity and water absolute permeability. The changes in two-phase flow properties appear coherent with the pore structure modifications induced by the carbonates dissolution but the changes in relative permeability could also be explained by a potential increase of the water-wet character of the core. Further investigations on the impacts of mineral changes are required with other reactive formation rocks, especially carbonate-rich ones, because the implications can be significant both for the validity of laboratory measurements and for the outcomes of in situ operations modeling.

Pore-scale numerical investigations of fluid flow in porous media using lattice Boltzmann method

Purpose – Lattice Boltzmann method (LBM) is employed to explore friction factor of single-phase fluid flow through porous media and the effects of local porous structure including geometry of grains in porous media and specific surface of porous media on two-phase flow dynamic behavior, phase distribution and relative permeability. The paper aims to discuss this issue. Design/methodology/approach – The 3D single-phase LBM model and the 2D multi-component multi-phase Shan-Chen LBM model (S-C model) are developed for fluid flow through porous media. For the solid site, the bounce back scheme is used with non-slip boundary condition. Findings – The predicted friction factor for single-phase fluid flow agrees well with experimental data and the well-known correlation. Compared with porous media with square grains, the two-phase fluids in porous media with circle grains are more connected and continuous, and consequently the relative permeability is higher. As for the factor of specific porous media surface, t...

Numerical Investigation of the Flow Dynamics and Evaporative Cooling of Water Droplets Impinging onto Heated Surfaces: An Effective Approach To Identify Spray Cooling Mechanisms

Numerical investigations of the dynamics and evaporative cooling of water droplets impinging onto heated surfaces can be used to identify spray cooling mechanisms. Droplet impingement dynamics and evaporation are simulated using the presented numerical model. Volume-of-fluid method is used in the model to track the free surface. The contact line dynamics was predicted from a dynamic contact angle model with the evaporation rate predicted by a kinetic theory model. A species transport equation was solved in the gas phase to describe the vapor convection and diffusion. The numerical model was validated by experimental data. The physical effects including the contact angle hysteresis and the thermocapillary effect are analyzed to offer guidance for future numerical models of droplet impingement cooling. The effects of various parameters including surface wettability, surface temperature, droplet velocity, droplet size, and droplet temperature were numerically studied from the standpoint of spray cooling. The numerical simulations offer profound analysis and deep insight into the spray cooling heat transfer mechanisms.

Effect of reactive surface area of minerals on mineralization trapping of CO2 in saline aquifers

The reactive surface area, an important parameter controlling mineral reactions, affects the amount of mineralization trapping of CO2 which affects the long-term CO2 storage. The effect of the reactive surface area on the mineralization trapping of CO2 was numerically simulated for CO2 storage in saline aquifers. Three kinds of minerals, including anorthite, calcite and kaolinite, are involved in the mineral reactions. This paper models the relationship between the specific surface area and the grain diameter of anorthite based on experimental data from literature (Brantley and Mellott, 2000). When the reactive surface areas of anorthite and calcite decrease from 838 to 83.8 m2/m3, the percentage of mineralization trapping of CO2 after 500 years decreases from 11.8% to 0.65%. The amount of dissolved anorthite and the amounts of precipitated kaolinite and calcite decrease significantly when the reactive surface areas of anorthite and calcite decrease from 838 to 83.8 m2/m3. Calcite is initially dissolved in the brine and then precipitates during the geochemical reactions between CO2-H2O and the minerals. Different reactive surface areas of anorthite and calcite lead to different times from dissolution to precipitation. The pH of the brine decreases with decreasing reactive surface areas of anorthite and calcite which influences the acidity of the saline aquifer. The gas saturation between the upper and lower parts of the saline aquifer increases with decreasing reactive surface areas of anorthite and calcite. The mass density distribution of brine solution shows that the CO2+brine solution region increases with decreasing reactive surface areas of anorthite and calcite.

Water Droplet Spreading and Wicking on Nanostructured Surfaces

Phase-change heat transfer on nanostructured surfaces is an efficient cooling method for high heat flux devices due to its superior wettability. Liquid droplet spreading and wicking effect then dominate the heat transfer. Therefore, this study investigates the flow behavior after a droplet touches a nanostructured surface focusing on the ZnO nanowire surface with three different nanowire sizes and two array types (regular and irregular). The spreading diameter and the wicking diameter are measured against time. The results show that the average spreading and wicking velocities on a regular nanostructured surface are both smaller than those on an irregular nanostructured surface and that the nanowire size affects the liquid spreading and capillary wicking.

Effect of Mineral Dissolution/Precipitation and CO2 Exsolution on CO2 transport in Geological Carbon Storage

Geological carbon sequestration (GCS) in deep saline aquifers is an effective means for storing carbon dioxide to address global climate change. As the time after injection increases, the safety of storage increases as the CO2 transforms from a separate phase to CO2(aq) and HCO3- by dissolution and then to carbonates by mineral dissolution. However, subsequent depressurization could lead to dissolved CO2(aq) escaping from the formation water and creating a new separate phase which may reduce the GCS system safety. The mineral dissolution and the CO2 exsolution and mineral precipitation during depressurization change the morphology, porosity, and permeability of the porous rock medium, which then affects the two-phase flow of the CO2 and formation water. A better understanding of these effects on the CO2-water two-phase flow will improve predictions of the long-term CO2 storage reliability, especially the impact of depressurization on the long-term stability. In this Account, we summarize our recent work on the effect of CO2 exsolution and mineral dissolution/precipitation on CO2 transport in GCS reservoirs. We place emphasis on understanding the behavior and transformation of the carbon components in the reservoir, including CO2(sc/g), CO2(aq), HCO3-, and carbonate minerals (calcite and dolomite), highlight their transport and mobility by coupled geochemical and two-phase flow processes, and consider the implications of these transport mechanisms on estimates of the long-term safety of GCS. We describe experimental and numerical pore- and core-scale methods used in our lab in conjunction with industrial and international partners to investigate these effects. Experimental results show how mineral dissolution affects permeability, capillary pressure, and relative permeability, which are important phenomena affecting the input parameters for reservoir flow modeling. The porosity and the absolute permeability increase when CO2 dissolved water is continuously injected through the core. The MRI results indicate dissolution of the carbonates during the experiments since the porosity has been increased after the core-flooding experiments. The mineral dissolution changes the pore structure by enlarging the throat diameters and decreasing the pore specific surface areas, resulting in lower CO2/water capillary pressures and changes in the relative permeability. When the reservoir pressure decreases, the CO2 exsolution occurs due to the reduction of solubility. The CO2 bubbles preferentially grow toward the larger pores instead of toward the throats or the finer pores during the depressurization. After exsolution, the exsolved CO2 phase shows low mobility due to the highly dispersed pore-scale morphology, and the well dispersed small bubbles tend to merge without interface contact driven by the Ostwald ripening mechanism. During depressurization, the dissolved carbonate could also precipitate as a result of increasing pH. There is increasing formation water flow resistance and low mobility of the CO2 in the presence of CO2 exsolution and carbonate precipitation. These effects produce a self-sealing mechanism that may reduce unfavorable CO2 migration even in the presence of sudden reservoir depressurization.

Numerical Investigation of Jet Impingement Cooling of a Flat Plate with Carbon Dioxide at Supercritical Pressures

ABSTRACT Confined round jet impingement cooling of a flat plate at constant heat flux with carbon dioxide at supercritical pressures was investigated numerically. The pressure ranged from 7.8 to 10.0 MPa, which is greater than the critical pressure of carbon dioxide, 7.38 MPa. The inlet temperature varied from 270 to 320 K and the heat flux ranged from 0.6 to 1.6 MW/m2. The shear-stress transport turbulence model was used and the numerical model was validated by comparison with experimental results for jet impingement heating with hot water at supercritical pressures. Radial conduction in the jet impingement plate was also considered. The sharp variations of the thermal-physical properties of the fluid near the pseudocritical point significantly influence heat transfer on the target wall. For a given heat flux, the high specific heat near the wall for the proper inlet temperature and pressure maximizes the average heat transfer coefficient. For a given inlet temperature, the heat transfer coefficient remains almost unchanged with increasing surface heat flux at first and then decreases rapidly as the heat flux becomes higher due to the combined effects of the thinner high specific heat layer and the smaller thermal conductivity at higher temperature.

CO2 Exsolution from CO2 Saturated Water: Core-Scale Experiments and Focus on Impacts of Pressure Variations

For CO2 sequestration and utilization in the shallow reservoirs, reservoir pressure changes are due to the injection rate changing, a leakage event, and brine withdrawal for reservoir pressure balance. The amounts of exsolved CO2 which are influenced by the pressure reduction and the subsequent secondary imbibition process have a significant effect on the stability and capacity of CO2 sequestration and utilization. In this study, exsolution behavior of the CO2 has been studied experimentally using a core flooding system in combination with NMR/MRI equipment. Three series of pressure variation profiles, including depletion followed by imbibitions without or with repressurization and repetitive depletion and repressurization/imbibition cycles, were designed to investigate the exsolution responses for these complex pressure variation profiles. We found that the exsolved CO2 phase preferentially occupies the larger pores and exhibits a uniform spatial distribution. The mobility of CO2 is low during the imbibition process, and the residual trapping ratio is extraordinarily high. During the cyclic pressure variation process, the first cycle has the largest contribution to the amount of exsolved CO2. The low CO2 mobility implies a certain degree of self-sealing during a possible reservoir depletion.

Experimental Investigation on the Behavior of Supercritical CO2 during Reservoir Depressurization

CO2 sequestration in saline aquifers is a promising way to address climate change. However, the pressure of the sequestration reservoir may decrease in practice, which induces CO2 exsolution and expansion in the reservoir. In this study, we conducted a core-scale experimental investigation on the depressurization of CO2-containing sandstone using NMR equipment. Three different series of experiments were designed to investigate the influence of the depressurization rate and the initial CO2 states on the dynamics of different trapping mechanisms. The pressure range of the depressurization was from 10.5 to 4.0 MPa, which covered the supercritical and gaseous states of the CO2 (named as CO2(sc) and CO2(g), respectively). It was found that when the aqueous phase saturated initially, the exsolution behavior strongly depended on the depressurization rate. When the CO2 and aqueous phase coexisting initially, the expansion of the CO2(sc/g) contributed to the incremental CO2 saturation in the core only when the CO2 occurred as residually trapped. It indicates that the reservoir depressurization has the possibility to convert the solubility trapping to the residual trapping phase, and/or convert the residual trapping to mobile CO2.

General slip regime permeability model for gas flow through porous media

A theoretical effective gas permeability model was developed for rarefied gas flow in porous media, which holds over the entire slip regime with the permeability derived as a function of the Knudsen number. This general slip regime model (GSR model) is derived from the pore-scale Navier-Stokes equations subject to the first-order wall slip boundary condition using the volume-averaging method. The local closure problem for the volume-averaged equations is studied analytically and numerically using a periodic sphere array geometry. The GSR model includes a rational fraction function of the Knudsen number which leads to a limit effective permeability as the Knudsen number increases. The mechanism for this behavior is the viscous fluid inner friction caused by converging-diverging flow channels in porous media. A linearization of the GSR model leads to the Klinkenberg equation for slightly rarefied gas flows. Finite element simulations show that the Klinkenberg model overestimates the effective permeability by as much as 33% when a flow approaches the transition regime. The GSR model reduces to the unified permeability model [F. Civan, “Effective correlation of apparent gas permeability in tight porous media,” Transp. Porous Media 82, 375 (2010)] for the flow in the slip regime and clarifies the physical significance of the empirical parameter b in the unified model.