Heavy metals released from kerogen to produced water during oil/gas extraction have caused major environmental concerns. To curtail water usage and production in an operation and to use the same process for carbon sequestration, supercritical CO2 (scCO2) has been suggested as a fracking fluid or an oil/gas recovery agent. It has been shown previously that injection of scCO2 into a reservoir may cause several chemical and physical changes to the reservoir properties including pore surface wettability, gas sorption capacity, and transport properties. Using molecular dynamics simulations, we here demonstrate that injection of scCO2 might lead to desorption of physically adsorbed metals from kerogen structures. This process on one hand may impact the quality of produced water. On the other hand, it may enhance metal recovery if this process is used for in-situ extraction of critical metals from shale or other organic carbon-rich formations such as coal.
Swelling clay hydration/dehydration is important to many environmental and industrial processes. Experimental studies usually probe equilibrium hydration states in an averaged manner and thus cannot capture the fast water transport and structural change in interlayers during hydration/dehydration. Using molecular simulations and thermogravimetric analyses, we observe a two-stage dehydration process. The first stage is controlled by evaporation at the edges: water molecules near hydrophobic sites and the first few water molecules of the hydration shell of cations move fast to particle edges for evaporation. The second stage is controlled by slow desorption of the last 1-2 water molecules from the cations and slow transport through the interlayers. The two-stage dehydration is strongly coupled with interlayer collapse and the coordination number changes of cations, all of which depend on layer charge distribution. This mechanistic interpretation of clay dehydration can be key to the coupled chemomechanical behavior in natural/engineered barriers.
Using molecular dynamics simulations, we investigate the molecular scale origin of crystal face selectivity when one gibbsite particle attaches to another in water. A comparison of the free energy per unit surface area of particle–particle attachment indicates that particle attachment through edge surfaces, where the edge surfaces are either (1 0 0) or (1 1 0) crystal faces, is more energetically favorable compared to attachment between two basal surfaces (i.e., (0 0 1) crystal faces) or between the basal surface of one particle and the edge surface of another. This result suggests that gibbsite crystals with low basal/edge surface area ratio will preferentially attach through edge surfaces, potentially helping the crystals grow laterally. However, for larger gibbsite particles (high basal/edge surface area ratio) the total free energy, not normalized by surface area, of particle attachment through the basal surfaces is lower (more negative) than attachment through the edge surfaces, indicating that larger gibbsite particles will preferentially aggregate through basal surface attachments. The short-range electrostatic interactions including the interparticle hydrogen bonds from surface –OH groups drive particle attachment, and the dominant contribution to the free energy minimum is enthalpic rather than entropic. However, the enthalpy of basal-edge attachment is significantly offset by the entropy leading to a higher free energy (less negative) compared to that of basal-basal attachment. Study of the free energy for a few imperfect attachments of two particles indicates a higher free energy (i.e., less negative, less stable), compared to a perfect attachment
Upon extraction/injection of a large quantity of gas from/into a subsurface system in shale gas production or carbon sequestration, the gas pressure varies remarkably, which may significantly change the wettability of porous media involved. Mechanistic understanding of such changes is critical for designing and optimizing a related subsurface engineering process. Using molecular dynamics simulations, we have calculated the contact angle of a water droplet on various solid surfaces (kerogen, pyrophyllite, calcite, gibbsite, and montmorillonite) as a function of CO2 or CH4 gas pressure up to 200 atm at a temperature of 300 K. The calculation reveals a complex behavior of surface wettability alteration by gas pressure variation depending on surface chemistry and structure, and molecular interactions of fluid molecules with surfaces. As the CO2 gas pressure increases, a partially hydrophilic kerogen surface becomes highly hydrophobic, while a calcite surface becomes more hydrophilic. Considering kerogen and calcite being the major components of a shale formation, we postulate that the wettability alteration of a solid surface induced by a gas pressure change may play an important role in fluid flows in shale gas production and geological carbon sequestration.
The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heatgenerating nuclear waste disposition in deep-seated clay/shale/argillaceous rock. International collaboration activities such as heater tests, continuous field data monitoring, and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum and discrete) to tackle issues related to flow and transport at various scales of the host-rock, its interactions with barrier materials, and EBS design concept.
This interim report is an update of ongoing experimental and modeling work on bentonite material described in Jové Colón et al. (2019, 2020) from past international collaboration activities. As noted in Jové Colón et al. (2020), work on international repository science activities such as FEBEX-DP and DECOVALEX19 is either no longer continuing by the international partners. Nevertheless, research activities on the collected sample materials and field data are still ongoing. Descriptions of these underground research laboratory (URL) R&D activities are described elsewhere (Birkholzer et al. 2019; Jové Colón et al. 2020) but will be explained here when needed. The current reports recent reactive-transport modeling on the leaching of sedimentary rock.
Water flow in nanometer or sub-nanometer hydrophilic channels bears special importance in diverse fields of science and engineering. However, the nature of such water flow remains elusive. Here, we report our molecular-modeling results on water flow in a sub-nanometer clay interlayer between two montmorillonite layers. We show that a fast advective flow can be induced by evaporation at one end of the interlayer channel, that is, a large suction pressure created by evaporation (∼818 MPa) is able to drive the fast water flow through the channel (∼0.88 m/s for a 46 Å-long channel). Scaled up for the pressure gradient to a 2 μm particle, the velocity of water is estimated to be about 95 μm/s, indicating that water can quickly flow through a μm-sized clay particle within seconds. The prediction seems to be confirmed by our thermogravimetric analysis of bentonite hydration and dehydration processes, which indicates that water transport at the early stage of the dehydration is a fast advective process, followed by a slow diffusion process. The possible occurrence of a fast advective water flow in clay interlayers prompts us to reassess water transport in a broad set of natural and engineered systems such as clay swelling/shrinking, moisture transport in soils, water uptake by plants, water imbibition/release in unconventional hydrocarbon reservoirs, and cap rock integrity of supercritical CO2 storage.
In shale gas production, gas composition may vary over time. To understand this phenomenon, we use molecular dynamics simulations to study the permeation of CH4, C2H6 and their mixture from a source container through a pyrophyllite nanopore driven by a pressure gradient. For a pure gas, the flow rate of CH4 is always higher than that of C2H6, regardless of pore size. For a 1:1 C2H6: CH4 mixture, however, C2H6:CH4 flow rate ratio is higher than the compositional ratio in the container (i.e., 1:1) when the pore size is smaller than ~1.8 nm. The selective transport is caused by the competitive adsorption of C2H6 over CH4 in the nanopore. The selectivity is also determined by the interplay between the surface diffusion of the adsorbed molecules and the viscous flow in the center of the pore, and it diminishes as the viscous flow becomes to dominate the surface diffusion when the pore size becomes larger than 1.8 nm. Our work shows that compositional differentiation of shale gas in production is a consequence of nanopore confinement and therefore a key characteristic of an unconventional reservoir. The related compositional information can potentially be used for monitoring the status of a production well such as its recovery rate.
A molecular-scale understanding of the transition between hydration states in clay minerals remains a challenging problem because of the very fast stepwise swelling process observed from X-ray diffraction (XRD) experiments. XRD profile modeling assumes the coexistence of multiple hydration states in a clay sample to fit the experimental XRD pattern obtained under humid conditions. While XRD profile modeling provides a macroscopic understanding of the heterogeneous hydration structure of clay minerals, a microscopic model of the transition between hydration states is still missing. Here, for the first time, we use molecular dynamics simulation to investigate the transition states between a dry interlayer, one-layer hydrate, and two-layer hydrate. We find that the hydrogen bonds that form across the interlayer at the clay particle edge make an important contribution to the energy barrier to interlayer hydration, especially for initial hydration.