On the reliability of analytical models to predict solute transport in a fracture network

Abstract. In hydrogeology, the application of reliable tracer transport model approaches is a key issue to derive the hydrodynamic properties of aquifers. Laboratory- and field-scale tracer dispersion breakthrough curves (BTC) in fractured media are notorious for exhibiting early time arrivals and late time tailing that are not captured by the classical advection–dispersion equation (ADE). These "non-Fickian" features are proven to be better explained by a mobile–immobile (MIM) approach. In this conceptualization the fractured rock system is schematized as a continuous medium in which the liquid phase is separated into flowing and stagnant regions. The present study compares the performances and reliabilities of the classical MIM and the explicit network model (ENM), taking expressly into account the network geometry for describing tracer transport behavior in a fractured sample at bench scale. Though ENM shows better fitting results than MIM, the latter remains still valid as it proves to describe the observed curves quite well. The results show that the presence of nonlinear flow plays an important role in the behavior of solute transport. First, the distribution of solute according to different pathways is not constant, but it is related to the flow rate. Second, nonlinear flow influences advection in that it leads to a delay in solute transport respect to the linear flow assumption. However, nonlinear flow is not shown to be related with dispersion. The experimental results show that in the study case the geometrical dispersion dominates the Taylor dispersion. However, the interpretation with the ENM shows a weak transitional regime from geometrical dispersion to Taylor dispersion for high flow rates. Incorporating the description of the flow paths in the analytical modeling has proven to better fit the curves and to give a more robust interpretation of the solute transport

Evidence of non-Darcy flow and non-Fickian transport in fractured media at laboratory scale

Abstract. During a risk assessment procedure as well as when dealing with cleanup and monitoring strategies, accurate predictions of solute propagation in fractured rocks are of particular importance when assessing exposure pathways through which contaminants reach receptors. Experimental data obtained under controlled conditions such as in a laboratory allow to increase the understanding of the fundamental physics of fluid flow and solute transport in fractures. In this study, laboratory hydraulic and tracer tests have been carried out on an artificially created fractured rock sample. The tests regard the analysis of the hydraulic loss and the measurement of breakthrough curves for saline tracer pulse inside a rock sample of parallelepiped shape (0.60 × 0.40 × 0.08 m). The convolution theory has been applied in order to remove the effect of the acquisition apparatus on tracer experiments. The experimental results have shown evidence of a non-Darcy relationship between flow rate and hydraulic loss that is best described by Forchheimer's law. Furthermore, in the flow experiments both inertial and viscous flow terms are not negligible. The observed experimental breakthrough curves of solute transport have been modeled by the classical one-dimensional analytical solution for the advection–dispersion equation (ADE) and the single rate mobile–immobile model (MIM). The former model does not properly fit the first arrival and the tail while the latter, which recognizes the existence of mobile and immobile domains for transport, provides a very decent fit. The carried out experiments show that there exists a pronounced mobile–immobile zone interaction that cannot be neglected and that leads to a non-equilibrium behavior of solute transport. The existence of a non-Darcian flow regime has showed to influence the velocity field in that it gives rise to a delay in solute migration with respect to the predicted value assuming linear flow. Furthermore, the presence of inertial effects enhance non-equilibrium behavior. Instead, the presence of a transitional flow regime seems not to exert influence on the behavior of dispersion. The linear-type relationship found between velocity and dispersion demonstrates that for the range of imposed flow rates and for the selected path the geometrical dispersion dominates the mixing processes along the fracture network

Laboratory experimental investigation of heat transport in fractured media

Low enthalpy geothermal energy is a renewable resource that is still underexploited nowadays in relation to its potential for development in society worldwide. Most of its applications have already been investigated, such as heating and cooling of private and public buildings, road defrosting, cooling of industrial processes, food drying systems or desalination. Geothermal power development is a long, risky and expensive process. It basically consists of successive development stages aimed at locating the resources (exploration), confirming the power generating capacity of the reservoir (confirmation) and building the power plant and associated structures (site development). Different factors intervene in influencing the length, difficulty and materials required for these phases, thereby affecting their cost. One of the major limitations related to the installation of low enthalpy geothermal power plants regards the initial development steps that are risky and the upfront capital costs that are huge. Most of the total cost of geothermal power is related to the reimbursement of invested capital and associated returns. In order to increase the optimal efficiency of installations which use groundwater as a geothermal resource, flow and heat transport dynamics in aquifers need to be well characterized. Especially in fractured rock aquifers these processes represent critical elements that are not well known. Therefore there is a tendency to oversize geothermal plants. In the literature there are very few studies on heat transport, especially on fractured media. This study is aimed at deepening the understanding of this topic through heat transport experiments in fractured networks and their interpretation. Heat transfer tests have been carried out on the experimental apparatus previously employed to perform flow and tracer transport experiments, which has been modified in order to analyze heat transport dynamics in a network of fractures. In order to model the obtained thermal breakthrough curves, the Explicit Network Model (ENM) has been used, which is based on an adaptation of Tang’s solution for the transport of the solutes in a semi-infinite single fracture embedded in a porous matrix. Parameter estimation, time moment analysis, tailing character and other dimensionless parameters have permitted a better understanding of the dynamics of heat transport and the efficiency of heat exchange between the fractures and the matrix. The results have been compared with the previous experimental studies on solute transport

Impacts of desalinated and recycled water in the Abu Dhabi surficial aquifer

In Abu Dhabi, one of the most arid regions in the world, in recent decades, desalinated water has been identified as a prime solution in solving the water demand issues. In this study, a three-dimensional coupled density-dependent flow and solute transport model was set up in order to study the effect of the artificial recharge using desalinated water and the influence of nonconventional water with a salt concentration in the range 0.1–2 g/L The results confirm that this region demands the adoption of a more rational use of irrigation water or additional usage of desalinated water and recycled water together with optimizing groundwater pumping at locations that are vulnerable to further quality degradation and depletion. The long-term storage of desalinated freshwater with a maximum radial distance of 653 m in the dune surface is ensured with the formation of the transition zone, and change in the groundwater head up to 5 km. The maximum recovery obtained by immediate recovery is 70%. The study expresses the long-term feasibility of desalinated freshwater storage and the need for further management practices in quantifying the contribution of desalinated and recycled water for agriculture activities which might have improved groundwater quality and increased hydraulic head at some locations

Experimental study of forced convection heat transport in porous media

Abstract. The present study is aimed at extending this thematic issue through heat transport experiments and their interpretation at laboratory scale. An experimental study to evaluate the dynamics of forced convection heat transfer in a thermally isolated column filled with porous medium has been carried out. The behavior of two porous media with different grain sizes and specific surfaces has been observed. The experimental data have been compared with an analytical solution for one-dimensional heat transport for local nonthermal equilibrium condition. The interpretation of the experimental data shows that the heterogeneity of the porous medium affects heat transport dynamics, causing a channeling effect which has consequences on thermal dispersion phenomena and heat transfer between fluid and solid phases, limiting the capacity to store or dissipate heat in the porous medium

Geostatistic application for the study of contaminated groundwater in broad industrial areas

A geostatistic application is proposed for the realization of a hydraulic conductivity map in an abandoned industrial area, characterized by the extension of about 300 ha, in which a characterization study is in progress. The study has made it possible to define the distribution of the grade of vulnerability of the area and, in relation to the knowledge of the centres of danger, has given further elements for the individuation of hot spots and the paths of contamination. In order to obtain a map of the hydraulic conductivity, through investigation of a direct nature, experimental values of permeability and trasmissivity have been subjected to specific experimental checks, and to calibration made through theoretic hydrogeologic approaches in order to establish the reliability of the data that have been given. The results, which are the object of verification, have been treated afterwards through a geostatistic approach in order to get a punctual valuation of the hydraulic conductivity, also in non-sampled points which are inside the examined area

Bench scale laboratory tests to analyze non-linear flow in fractured media

The knowledge of flow phenomena in fractured rocks is very important for groundwater resources management in hydrogeological engineering. brbr A critical emerging issue for fractured aquifers is the validity of the Darcian-type "local cubic law", which assumes a linear relationship between flow rate and pressure gradient to accurately describe flow patterns. brbr Experimental data obtained under controlled conditions such as in a laboratory increase our understanding of the fundamental physics of fracture flow and allow us to investigate the presence of non-linear flow inside fractures that generates a substantial deviation from Darcy's law. brbr In this study the presence of non-linear flow in a fractured rock formation has been analyzed at bench scale in laboratory tests. The effects of non-linearity in flow have been investigated by analyzing hydraulic tests on an artificially created fractured rock sample of parallelepiped (0.60 Ã\u97 0.40 Ã\u97 0.8 m) shape. brbr The volumes of water passing through different paths across the fractured sample for various hydraulic head differences have been measured, and the results of the experiments have been reported as specific flow rate vs. head gradient. The experimental results closely match the Forchheimer equation and describe a strong inertial regime. The results of the test have been interpreted by means of numerical simulations. For each pair of ports, several steady-state simulations have been carried out varying the hydraulic head difference between the inlet and outlet ports. The estimated linear and non-linear Forchheimer coefficients have been correlated to each other and respectively to the tortuosity of the flow paths. A correlation among the linear and non-linear Forchheimer coefficients is evident. Moreover, a tortuosity factor that influences flow dynamics has been determined. © Author(s) 2013

Nuovi dati dall'insula 104: i contesti della Casa dell'Iscrizione Dipinta (VII-IX d.C.)

In hydrogeology, the application of reliable tracer transport model approaches is a key issue to derive the hydrodynamic properties of aquifers. Laboratory- and field-scale tracer dispersion breakthrough curves (BTC) in fractured media are notorious for exhibiting early time arrivals and late time tailing that are not captured by the classical advection–dispersion equation (ADE). These "non-Fickian" features are proven to be better explained by a mobile–immobile (MIM) approach. In this conceptualization the fractured rock system is schematized as a continuous medium in which the liquid phase is separated into flowing and stagnant regions. The present study compares the performances and reliabilities of the classical MIM and the explicit network model (ENM), taking expressly into account the network geometry for describing tracer transport behavior in a fractured sample at bench scale. Though ENM shows better fitting results than MIM, the latter remains still valid as it proves to describe the observed curves quite well. The results show that the presence of nonlinear flow plays an important role in the behavior of solute transport. First, the distribution of solute according to different pathways is not constant, but it is related to the flow rate. Second, nonlinear flow influences advection in that it leads to a delay in solute transport respect to the linear flow assumption. However, nonlinear flow is not shown to be related with dispersion. The experimental results show that in the study case the geometrical dispersion dominates the Taylor dispersion. However, the interpretation with the ENM shows a weak transitional regime from geometrical dispersion to Taylor dispersion for high flow rates. Incorporating the description of the flow paths in the analytical modeling has proven to better fit the curves and to give a more robust interpretation of the solute transport