RESUMEN
Transport networks, such as vasculature or river networks, provide key functions in organisms and the environment. They usually contain loops whose significance for the stability and robustness of the network is well documented. However, the dynamics of their formation is usually not considered. Such structures often grow in response to the gradient of an external field. During evolution, extending branches compete for the available flux of the field, which leads to effective repulsion between them and screening of the shorter ones. Yet, in remarkably diverse processes, from unstable fluid flows to the canal system of jellyfish, loops suddenly form near the breakthrough when the longest branch reaches the boundary of the system. We provide a physical explanation for this universal behavior. Using a 1D model, we explain that the appearance of effective attractive forces results from the field drop inside the leading finger as it approaches the outlet. Furthermore, we numerically study the interactions between two fingers, including screening in the system and its disappearance near the breakthrough. Finally, we perform simulations of the temporal evolution of the fingers to show how revival and attraction to the longest finger leads to dynamic loop formation. We compare the simulations to the experiments and find that the dynamics of the shorter finger are well reproduced. Our results demonstrate that reconnection is a prevalent phenomenon in systems driven by diffusive fluxes, occurring both when the ratio of the mobility inside the growing structure to the mobility outside is low and near the breakthrough.
RESUMEN
Salt solubility is generally determined under isotropic stress conditions. Yet, in the context of salt weathering of porous media, mechanical constraints on the in-pore growth of salt crystals are likely to be orientation-dependent, resulting in an anisotropic stress state on the crystal. In this paper, we determine by molecular simulation the solubility of NaCl in water when the crystal is subjected to anisotropic stress. Such anisotropy causes the chemical potential of the crystal to be orientation-dependent, and proper thermodynamic formulation requires describing the chemical potential as a tensor. The solute and crystal chemical potentials are computed from free energy calculations using Hamiltonian thermodynamic integration, and the usual condition of solubility is reformulated to account for the tensorial nature of the crystal chemical potential. We investigate in detail how the uniaxial compression of the crystal affects its solubility. The molecular simulation results led to revisiting the Correns law under anisotropic stress. Regarding the solute, the non-ideal behavior of the liquid phase is captured using Pitzer's ion interaction approach up to high concentrations of interest for in-pore crystallization and beyond the concentrations addressed in the existing literature. Regarding NaCl crystals, the validity of the generalized Gibbs-Duhem equation for a tensorial chemical potential is carefully verified, and it is found that crystallization progresses almost orthogonally to the crystal surface even under high shear stresses. Comparing uniaxial and isotropic compression highlights the major differences in solubility caused by stress anisotropy, and the revisited Correns law offers an appropriate framework to capture this phenomenon.
RESUMEN
Aqueous geochemistry datasets from regional groundwater monitoring programs can be a major asset for environmental baseline assessment (EBA) in regions with development of natural gases from unconventional hydrocarbon resources. However, they usually do not include crucial parameters for EBA in areas of shale gas development such as methane concentrations. A logistic regression (LR) model was developed to predict the probability of methane occurrence in aquifers in Alberta (Canada). The model was calibrated and tested using geochemistry data including methane concentrations from two groundwater monitoring programs. The LR model correctly predicts methane occurrence in 89.8% (n = 234 samples) and 88.1% (n = 532 samples) of groundwater samples from each monitoring program. Methane concentrations strongly depend on the occurrence of electron donors such as sulfate and to a lesser extent on well depth and the total dissolved solids of groundwater. The model was then applied to a province-wide public health groundwater monitoring program (n = 52,849 samples) providing aqueous geochemistry data but no methane concentrations. This approach allowed the prediction of methane occurrence in regions where no groundwater gas data are available, thereby increasing the resolution of EBA in areas of shale gas development by using basic hydrochemical parameters measured in high-density groundwater monitoring programs.
Asunto(s)
Agua Subterránea , Contaminantes Químicos del Agua , Alberta , Monitoreo del Ambiente , Metano , Gas Natural , Yacimiento de Petróleo y GasRESUMEN
Due to increasing concerns over the potential impact of shale gas and coalbed methane (CBM) development on groundwater resources, it has become necessary to develop reliable tools to detect any potential pollution associated with hydrocarbon exploitation from unconventional reservoirs. One of the key concepts for such monitoring approaches is the establishment of a geochemical baseline of the considered groundwater systems. However, the detection of methane is not enough to assess potential impact from CBM and shale gas exploitation since methane in low concentrations has been found to be naturally ubiquitous in many groundwater systems. The objective of this study was to determine the methane sources, the extent of potential methane oxidation, and gas-water-rock-interactions in shallow aquifers by integrating chemical and isotopic monitoring data of dissolved gases and aqueous species into a geochemical PHREEQC model. Using data from a regional groundwater observation network in Alberta (Canada), the model was designed to describe the evolution of the concentrations of methane, sulfate and dissolved inorganic carbon (DIC) as well as their isotopic compositions (δ34SSO4, δ13CCH4 and δ13CDIC) in groundwater subjected to different scenarios of migration, oxidation and in situ generation of methane. Model results show that methane migration and subsequent methane oxidation in anaerobic environments can strongly affect its concentration and isotopic fingerprint and potentially compromise the accurate identification of the methane source. For example elevated δ13CCH4 values can be the result of oxidation of microbial methane and may be misinterpreted as methane of thermogenic origin. Hence, quantification of the extent of methane oxidation is essential for determining the origin of methane in groundwater. The application of this model to aquifers in Alberta shows that some cases of elevated δ13CCH4 values were due to methane oxidation resulting in pseudo-thermogenic isotopic fingerprints of methane. The model indicated no contamination of shallow aquifers by deep thermogenic methane from conventional and unconventional hydrocarbon reservoirs under baseline conditions. The developed geochemical and multi-isotopic model describing the sources and fate of methane in groundwater is a promising tool for groundwater assessment purposes in areas with shale gas and coalbed methane development.