scholarly journals Feedback mechanisms between precipitation and dissolution reactions across randomly heterogeneous conductivity fields

2021 ◽  
Vol 25 (11) ◽  
pp. 5905-5915
Author(s):  
Yaniv Edery ◽  
Martin Stolar ◽  
Giovanni Porta ◽  
Alberto Guadagnini

Abstract. Our study investigates interplays between dissolution, precipitation, and transport processes taking place across randomly heterogeneous conductivity domains and the ensuing spatial distribution of preferential pathways. We do so by relying on a collection of computational analyses of reactive transport performed in two-dimensional systems where the (natural) logarithm of conductivity is characterized by various degrees of spatial heterogeneity. Our results document that precipitation and dissolution jointly take place in the system, with the latter mainly occurring along preferential flow paths associated with the conductivity field and the former being observed at locations close to and clearly separated from these. High conductivity values associated with the preferential flow paths tend to further increase in time, giving rise to a self-sustained feedback between transport and reaction processes. The clear separation between regions where dissolution or precipitation takes place is imprinted onto the sample distributions of conductivity which tend to become visibly left skewed with time (with the appearance of a bimodal behavior at some times). The link between conductivity changes and reaction-driven processes promotes the emergence of non-Fickian effective transport features. The latter can be captured through a continuous-time random-walk model where solute travel times are approximated with a truncated power law probability distribution. The parameters of such a model shift towards values associated with increasingly high non-Fickian effective transport behavior as time progresses.

2021 ◽  
Author(s):  
Yaniv Edery ◽  
Martin Stolar ◽  
Giovanni Porta ◽  
Alberto Guadagnini

Abstract. Our study investigates interplays between dissolution, precipitation, and transport processes taking place across randomly heterogeneous conductivity domains and the ensuing spatial distribution of preferential pathways. We do so by relying on a collection of computational analyses of reactive transport performed in two-dimensional systems where the (natural) logarithm of conductivity is characterized by various degrees of spatial heterogeneity. Our results document that precipitation and dissolution jointly take place in the system, the latter mainly occurring along preferential flowpaths associated with the conductivity field, the former being observed at locations close to and clearly separated from these. High conductivity values associated with the preferential flowpaths tend to further increase in time, giving rise to a self-sustained feedback between transport and reaction processes. The clear separation between regions where dissolution or precipitation takes place is imprinted onto the sample distributions of conductivity which tend to become visibly left skewed with time (with the appearance of a bimodal behavior at some times). The link between conductivity changes and reaction-driven processes promotes the emergence of non-Fickian effective transport features. The latter can be captured through a continuous time random walk model where solute travel times are approximated with a truncated power law probability distribution. The parameters of such a model shift towards values associated with increasingly high non-Fickian effective transport behavior as time progresses.


2016 ◽  
Author(s):  
Georg Kaufmann ◽  
Franci Gabrovšek ◽  
Douchko Romanov

Abstract. Soluble rocks such as limestone, anhydrite, and gypsum are characterised by their large secondary permeability, which results from the interaction of water circulating through the rock and dissolving the soluble fracture walls. This highly selective dissolution process enlarges the fractures to voids and eventually cavities, which then carry the majority of flow through an aquifer along preferential flow paths. We employ a numerical model describing the evolution of secondary porosity in a soluble rock to study the evolution of isolated fractures in different rock types. Our main focus is three-fold: The identification of shallow versus deep flow paths and their evolution for different rock types; the effect of precipitation of the dissolved material in the fracture; and finally the complication of fracture enlargement in fractures composed of several different soluble materials. Our results show that the evolution of fractures composed of limestone and gypsum is comparable, but the evolution time scale is drastically different. For anhydrite, owing to its difference in the kinetical rate law describing the removal of soluble rock, when compared to limestone and anhydrite, the evolution is even faster. Precipitation of the dissolved rock due to changes in the hydrochemical conditions can clog fractures fairly fast, thus changing the pattern of preferential pathways in the soluble aquifer, especially with depth. Finally, limestone fracture coated with gypsum, as frequently observed in caves, will result in a substantial increase in fracture enlargement with time, thus giving these fractures a hydraulic advantage over pure limestone fractures in their competition for capturing flow.


2005 ◽  
Vol 7 ◽  
pp. 41-44
Author(s):  
Bertel Nilsson ◽  
Jens Aamand ◽  
Ole Stig Jacobsen ◽  
René K. Juhler

Recent research on Danish groundwater has focused on clarifying the fate and transport of pesticides that leach through clayey till aquitards with low matrix permeability. Previously, these aquitards were considered as protective layers against contamination of underlying groundwater aquifers due to their low permeability characteristics. However, geological heterogeneities such as fractures and macropores have been recognised as preferential flow paths within low permeable clayey till (e.g. Beven & Germann 1982). The flow velocities within these preferential flow paths can be orders of magnitude higher than in the surrounding clay matrix and pose a major risk of transport of contaminants to the underlying aquifers (e.g. Nilsson et al. 2001). Previous studies of transport in fractured clayey till have focused on fully saturated conditions (e.g. Sidle et al. 1998; McKay et al. 1999). However, seasonal fluctuations of the groundwater table typically result in unsaturated conditions in the upper few metres of the clay deposits, resulting in different flow and transport conditions. Only a few experiments have examined the influence of unsaturated conditions on flow and solute (the dissolved inorganic and organic constituents) transport in fractured clayey till. These include smallscale laboratory column experiments on undisturbed soil monoliths (e.g. Jacobsen et al. 1997; Jørgensen et al. 1998), intermediate scale lysimeters (e.g. Fomsgaard et al. 2003) and field-scale tile drain experiments (e.g. Kjær et al. 2005). The different approaches each have limitations in terms of characterising flow and transport in fractured media. Laboratory studies of solute transport in soils (intact soil columns) are not exactly representative of field conditions due to variations in spatial variability and soil structure. In contrast, field studies hardly allow quantification of fluxes and mechanisms of transport. Column and lysimeter experiments are often limited in size, and tile-drain experiments on field scale do not provide spatial resolution and often have large uncertainties in mass balance calculations. Thus, in order to represent the overall natural fracture network systems on a field scale with respect to acquiring insights into flow and transport processes, the lysimeter needs to be larger than normal lysimeter size (< 1 m3). A modified large-scale lysimeter was therefore constructed by the Geological Survey of Denmark and Greenland (GEUS) at the Avedøre experimental field site 15 km south of Copenhagen (Fig. 1). This lysimeter consisted of an isolated block (3.5 ×3.5 ×3.3 m) of unsaturated fractured clayey till with a volume sufficient to represent the overall preferential flow paths (natural fracture network) within lowpermeable clayey till at a field scale.


2002 ◽  
Vol 66 (2) ◽  
pp. 347 ◽  
Author(s):  
Zhi Wang ◽  
Jianhang Lu ◽  
Laosheng Wu ◽  
Thomas Harter ◽  
William A. Jury

2005 ◽  
Vol 69 (2) ◽  
pp. 291-300 ◽  
Author(s):  
Young-Jin Kim ◽  
Christophe J. G. Darnault ◽  
Nathan O. Bailey ◽  
J.-Yves Parlange ◽  
Tammo S. Steenhuis

SPE Journal ◽  
2015 ◽  
Vol 20 (04) ◽  
pp. 689-700 ◽  
Author(s):  
S.. Ameen ◽  
A. Dahi Taleghani

Summary Injectivity loss is a common problem in unconsolidated-sand formations. Injection of water into a poorly cemented granular medium may lead to internal erosion, and consequently formation of preferential flow paths within the medium because of channelization. Channelization in the porous medium might occur when fluid-induced stresses become locally larger than a critical threshold and small grains are dislodged and carried away; hence, porosity and permeability of the medium will evolve along the induced flow paths. Vice versa, flowback during shut-in might carry particles back to the well and cause sand accumulation inside the well, and subsequently loss of injectivity. In most cases, to maintain the injection rate, operators will increase injection pressure and pumping power. The increased injection pressure results in stress changes and possibly further changes in channel patterns around the wellbore. Experimental laboratory studies have confirmed the presence of the transition from uniform Darcy flow to a fingered-pattern flow. To predict these phenomena, a model is needed to fill this gap by predicting the formation of preferential flow paths and their evolution. A model based on the multiphase-volume-fraction concept is used to decompose porosity into mobile and immobile porosities where phases may change spatially, evolve over time, and lead to development of erosional channels depending on injection rates, viscosity, and rock properties. This model will account for both particle release and suspension deposition. By use of this model, a methodology is proposed to derive model parameters from routine injection tests by inverse analysis. The proposed model presents the characteristic behavior of unconsolidated formation during fluid injection and the possible effect of injection parameters on downhole-permeability evolution.


Biofouling ◽  
2013 ◽  
Vol 29 (9) ◽  
pp. 1069-1086 ◽  
Author(s):  
Simona Bottero ◽  
Tomas Storck ◽  
Timo J. Heimovaara ◽  
Mark C.M. van Loosdrecht ◽  
Michael V. Enzien ◽  
...  

2017 ◽  
Vol 21 (11) ◽  
pp. 5503-5515 ◽  
Author(s):  
Hiroyuki Hirashima ◽  
Francesco Avanzi ◽  
Satoru Yamaguchi

Abstract. The heterogeneous movement of liquid water through the snowpack during precipitation and snowmelt leads to complex liquid water distributions that are important for avalanche and runoff forecasting. We reproduced the formation of capillary barriers and the development of preferential flow through snow using a three-dimensional water transport model, which was then validated using laboratory experiments of liquid water infiltration into layered, initially dry snow. Three-dimensional simulations assumed the same column shape and size, grain size, snow density, and water input rate as the laboratory experiments. Model evaluation focused on the timing of water movement, thickness of the upper layer affected by ponding, water content profiles and wet snow fraction. Simulation results showed that the model reconstructs relevant features of capillary barriers, including ponding in the upper layer, preferential infiltration far from the interface, and the timing of liquid water arrival at the snow base. In contrast, the area of preferential flow paths was usually underestimated and consequently the averaged water content in areas characterized by preferential flow paths was also underestimated. Improving the representation of preferential infiltration into initially dry snow is necessary to reproduce the transition from a dry-snow-dominant condition to a wet-snow-dominant one, especially in long-period simulations.


2013 ◽  
Vol 17 (3) ◽  
pp. 947-959 ◽  
Author(s):  
D. M. Krzeminska ◽  
T. A. Bogaard ◽  
J.-P. Malet ◽  
L. P. H. van Beek

Abstract. The importance of hydrological processes for landslide activity is generally accepted. However, the relationship between precipitation, hydrological responses and movement is not straightforward. Groundwater recharge is mostly controlled by the hydrological material properties and the structure (e.g., layering, preferential flow paths such as fissures) of the unsaturated zone. In slow-moving landslides, differential displacements caused by the bedrock structure complicate the hydrological regime due to continuous opening and closing of the fissures, creating temporary preferential flow paths systems for infiltration and groundwater drainage. The consecutive opening and closing of fissure aperture control the formation of a critical pore water pressure by creating dynamic preferential flow paths for infiltration and groundwater drainage. This interaction may explain the seasonal nature of the slow-moving landslide activity, including the often observed shifts and delays in hydrological responses when compared to timing, intensity and duration of precipitation. The main objective of this study is to model the influence of fissures on the hydrological dynamics of slow-moving landslide and the dynamic feedbacks between fissures, hydrology and slope stability. For this we adapt the spatially distributed hydrological and slope stability model (STARWARS) to account for geotechnical and hydrological feedbacks, linking between hydrological response of the landside and the dynamics of the fissure network and applied the model to the hydrologically controlled Super-Sauze landslide (South French Alps).


Sign in / Sign up

Export Citation Format

Share Document