Two hydrogeologic models explore reflux dolomitization using two outcrop datasets at different scales to constrain transient boundary conditions and heterogeneous petrophysical properties. A platform-scale petrophysical model of the Permian San Andres Formation was built from outcrop and subsurface data following a reservoir modeling approach that preserves outcrop heterogeneity and incorporates a sequence stratigraphic framework. This model was used as input for hydrogeological simulations of hypersaline fluid flow and solute transport during the accumulation and compaction of the platform. Boundary conditions change over time, as relative sealevel fluctuations drive sedimentation, depositional environment migration, topographic gradients, and location, size and salinity of the brine source. The potential volume and distribution of dolomite formed is inferred by a magnesium mass-balance. The composite result of reflux events at various orders of stratigraphic hierarchy is a complex dolomite pattern that resembles that observed on San Andres outcrops. Dolostone bodies across the platform may be generated by different combinations of favorable conditions, including proximity to the brine source, zones of higher permeability, permeability contrasts, and latent reflux. A meter-scale reactive transport model of the Albian Upper Glen Rose Formation simulates deposition of three high-frequency cycles punctuated by three brine reflux events. The simulator determines flow, solute and reactive transport along the flow paths, revealing the spatial and temporal distribution of calcite dissolution, and precipitation of dolomite and sulfate. The model recreates fully and partially dolomitized cycles within the time and lithological constrains on Glen Rose outcrops. Our results show that the distribution of dolomite within a high-frequency cycle may be the net result of intercycle processes, whereby dolomitizing fluids sourced from younger cycles flow across stratigraphically significant boundaries. We also show that variations in dolomite abundance and the unfulfilled dolomitization potential control the contemporaneous propagation of multiple dolomite fronts and the coalescence of discrete dolomite bodies. Results show that reflux is an effective and efficient mechanism to dolomitize carbonate formations that progresses simultaneously with sediment accumulation. Dolomitization is the cumulative result of many short-lived reflux events, sourced in different locations and times, and amalgamation of successive dolostone bodies. This model contrasts with previous studies that approached dolomitization of a carbonate platform as a discrete reflux event and current interpretations that relate dolomite bodies to their most immediate stratigraphic surfaces. / text

No description available.

Seawater-fresh groundwater natural equilibrium conditions in the aquifers of the Sel&ccedil / uk sub-basin may be disturbed by the discharge occurred during the past 30 years in order to supply water for irrigation and domestic purpose usage. Two dimensional density dependent cross sectional saturated flow and solute transport simulations were carried out to determine whether seawater intrusion has occurred in the Sel&ccedil / uk sub-basin due to the imposed discharge or the determined salt-water in the western section of the aquifers represents natural interface equilibrium conditions. The numerical simulation model was calibrated using field measurements. The results of the simulations suggest that the seawater intrusion has occurred in the study area since the pumping activity increased in the region. Five scenarios were simulated for future predictions: (1) Present recharge and discharge conditions are maintained, (2) Discharge increases at a rate of the municipality need as a result of increasing population, (3) Discharge decreases by 12% from the present value, (4) Discharge decreases by 25% from the present value, and (5) No discharge occurs or recharge amount equals to the present discharge amount condition supplied to the aquifer. These scenario results suggest that seawater intrusion in the Sel&ccedil / uk sub-basin would progress in the following years unless proper management measures are taken into consideration. KEYWORDS: seawater intrusion, density dependent groundwater flow, solute transport, contamination, Sel&ccedil / uk sub-basin.

The effect of in-situ or redistributed stress on solute transport in fractured rocks is one of the major concerns for many subsurface engineering problems. However, it remains poorly understood due to the difficulties in experiments and numerical modeling. The main aim of this thesis is to systematically investigate the influences of stress on solute transport in fractured rocks, at scales of single fractures and fracture networks, respectively. For a single fracture embedded in a porous rock matrix, a closed-form solution was derived for modeling the coupled stress-flow-transport processes without considering damage on the fracture surfaces. Afterwards, a retardation coefficient model was developed to consider the influences of damage of the fracture surfaces during shear processes on the solute sorption. Integrated with particle mechanics models, a numerical procedure was proposed to investigate the effects of gouge generation and microcrack development in the damaged zones of fracture on the solute retardation in single fractures. The results show that fracture aperture changes have a significant influence on the solute concentration distribution and residence time. Under compression, the decreasing matrix porosity can slightly increase the solute concentration. The shear process can increase the solute retardation coefficient by offering more sorption surfaces in the fracture due to gouge generation, microcracking and gouge crushing. To study the stress effects on solute transport in fracture systems, a hybrid approach combing the discrete element method for stress-flow simulations and a particle tracking algorithm for solute transport was developed for two-dimensional irregular discrete fracture network models. Advection, hydrodynamic dispersion and matrix diffusion in single fractures were considered. The particle migration paths were tracked first by following the flowing fluid (advection), and then the hydrodynamic dispersion and matrix diffusion were considered using statistic methods. The numerical results show an important impact of stress on the solute transport, by changing the solute residence time, distribution and travel paths. The equivalent dispersion coefficient is scale dependent in an asymptotic or exponential form without stress applied or under isotropic compression conditions. Matrix diffusion plays a dominant role in solute transport when the hydraulic gradient is small. Outstanding issues and main scientific achievements are also discussed. / QC 20111011

Dans cette thèse, nous avons étudié expérimentalement, numériquement et théoriquement la sédimentation d’un soluté en cellule de Hele-Shaw, modèle analogue de milieu poreux, avec un écoulement zonal horizontal. Une cellule, alimentée par deux grands réservoirs assurant un écoulement zonal uniforme, a été réalisée au laboratoire. La vitesse a été mesurée par PIV. L’injection d’une solution saline dans cet écoulement permet de visualiser la propagation du panache de pollution, et de mesurer le champ de concentration par LIF. Trois configurations typiques ont été observées : panache sans digitation (« stable »), avec une seule digitation (« faiblement instable »), et avec de nombreuses digitations (« fortement instable »). L’influence du contraste de densité, de la vitesse zonale, et du débit d’injection sur la stabilité du panache a été analysée. Ces résultats ont été comparés à des simulations numériques sous COMSOL Multiphysics®. Les champs de concentration simulés sont en accord quantitatif avec les champs expérimentaux. Ces simulations ont de plus permis d’étudier l’influence de l’anisotropie du tenseur de dispersion sur la stabilité du panache. Une étude théorique linéaire, effectuée à partir d’un écoulement de base unidirectionnel et quasi-stationnaire, est venue compléter ces résultats. Elle montre que le panache est toujours convectivement instable pour les valeurs des paramètres considérées, et permet de retrouver et de quantifier les principaux résultats mis en évidence lors des études expérimentale et numérique, comme l’effet stabilisant de la dispersivité longitudinale ou l’effet déstabilisant du contraste de densité / This thesis is an experimental and numerical investigation of solute transport and sedimentation in a Hele-Shaw cell with zonal flow. A vertical cell joining two large tanks has been built to provide a uniform horizontal flow, the velocity of which has been measured by using PIV. The concentration field of a solute injected in this flow has been measured by using LIF. Three typical plume configurations have been observed: without digitations (“stable”), with a single digitation (“weakly unstable”), and with various digitations appearing at the lower interface (“unstable plume”). The influence of density contrast, zonal flow and solute flow rate on the appearance of the various configurations has been investigated. These results have been compared to numerical simulations using the COMSOL Multiphysics® software. The numerical concentration fields are in qualitative agreement with the experimental ones. In addition, these simulations allowed us to investigate the effect of the anisotropy of the dispersion tensor on the occurrence of unstable configurations. A linear stability analysis based on a uniform parallel flow superposed to a quasi-steady horizontal layer of solute has been performed. It shows that the plume is convectively unstable for the parameters considered here, and confirms the various behaviors observed in both the experimental and the numerical analyses, like the stabilizing effect of the longitudinal dispersivity or the destabilizing effect of the solute concentration

Un dispositif en sol nu comprenant trois parcelles (90 m²) et six cases lysimétriques (4 m²) instrumentées a permis d’étudier l’influence d’apports répétés d’une boue de station d’épuration (BOUE) et de la même boue compostée avec des déchets verts (DVB) sur (i) les modifications des propriétés physico-chimiques et hydriques du sol ; (ii) le transport des ions bromures ; (iii) la dynamique de minéralisation du carbone et de l’azote organiques dans le sol. Finalement, le calage des paramètres hydrodynamiques et de décomposition de la matière organique dans le sol, propres à chaque horizon des profils de sols instrumentés, a permis une description satisfaisante de la dynamique hydrique et du transport des ions bromures et nitrates observés expérimentalement. Ces résultats aideront à l’adaptation de scénarios agronomiques sur le long terme pour répondre plus efficacement aux problématiques liées aux nouvelles exigences de recyclage des déchets et de protection de l’environnement. / An instrumented experimental bare soil site including three plots (90 m²) and six lysimeters (4 m²) has enabled to study the influence of repeated sewage sludge (SS) applications, as well as the same sewage sludge composted with green waste (SGW), on (i) the soil physicochemical and hydric properties ; (ii) the transport of bromide ions ; (iii) the organic carbon and nitrogen mineralization dynamics in soil. Finally, fitting the governing hydraulic parameters and the organic matter decomposition parameters in each soil horizon allowed a satisfying description of the experimentally observed water and bromide and nitrate ions transport in soil. These results will help conceive long-term agronomical scenarios aimed at providing more efficient answers to the new requirements of environmental protection and waste recycling policies.

Die Modellierung von komplexen Systemen, wie dem Untergrund, ist ein Hilfsmittel zur Beschreibung der in der Realität ablaufenden Prozesse. Die Durchführung von Experimenten an einem Modell, um qualitative Aussagen über ein reales System zu erhalten, wird als Simulation bezeichnet. Dabei können vielfältige Modelle, wie z.B. physikalische und mathematische, zum Einsatz kommen. Die ungesättigte Bodenzone (vadose Zone) bezeichnet den Bereich zwischen der Landoberfläche und dem Grundwasserspiegel, innerhalb dessen der Wassergehalt geringer als bei Vollsättigung, und der Druck geringer als der Atmosphärendruck ist. Dieser Bodenbereich hat für die Landwirtschaft, Geobiologie, aerobe Abbauprozesse und Grundwasserneubildung eine große Bedeutung. Für die Nachbildung von Strömungs- und Stofftransportprozessen der ungesättigten Bodenzone existieren numerische Simulationsprogramme. Ziel der vorliegenden Arbeit ist eine umfangreiche Validierung des Programms PCSiWaPro® (entwickelt an der TU-Dresden, Institut für Abfallwirtschaft und Altlasten) für unterschiedliche Anwendungsfälle. Ein weiteres Ziel der Arbeit besteht in der Untersuchung der Anwendbarkeit des aktuellen Stands des Simulationsprogramms PCSiWaPro® auf unterschiedliche Praxisfälle bei Kombination der ungesättigten Bodenzone mit technischen Anlagen. Vier Anwendungsfälle mit unterschiedlichen Zielen wurden dafür im Rahmen dieser Arbeit untersucht: die Simulation von dezentraler Abwasserversickerung (Kleinkläranlage - KKA) anhand entsprechender Säulen- und Feldversuche, die Berechnung der Grundwasserneubildung am Beispiel von Lysimetern, der Wasserhaushalt von Erddämmen und die Modellierung von Deponieabdeckungssystemen. Die Anwendungsfälle unterscheiden sich durch den Zweck der Simulation, die Geometrie, die Größe, die festgelegten Anfangs- und Randbedingungen, die Simulationszeit, die Materialien, das Koordinatensystem sowie die Ein- und Ausgabewerte. Die Simulationsergebnisse konnten eindeutig zeigen, dass das Programm PCSiWaPro® für alle im Rahmen der vorliegenden Arbeit untersuchten Fälle, mit unterschiedlichen Strömungsregimen, Stofftransport-Parametern, Randbedingungen, Koordinatensystemen sowie Raum- und Zeitdiskretisierungen anwendbar ist. Die Simulationsergebnisse der Säulenversuche am Beispiel dezentraler Abwasserversickerung zeigten eine sehr gute Übereinstimmung zwischen gemessenen und mittels PCSiWa-Pro® berechneten Werten des Wasser- und Stoffhaushaltes (Druckhöhe, Abfluss und Stoff-konzentration) der untersuchten Bodentypen B3 (schwachschluffiger Sand), B4 (Grobsand) und B5 (mittelschluffiger Sand). Die Wurzel des mittleren quadratischen Fehlers (RMSE) betrug für die Berechnung der Druckhöhe 1,84 cm bei B5, 3,61 cm bei B3 und 1,27 cm bei B4. Die relative Abweichung betrug für die Berechnung der Druckhöhe 2,19 % bei B5, 1,3 % bei B3 und ca. 5,3 % bei B4. Die Durchführung der Sensitivitätsanalyse der für die Modellierung relevanten Parameter zeigte eine sehr hohe Sensitivität der VAN GENUCHTEN-Parameter und der gesättigten hydraulischen Leitfähigkeit des Bodens. Darüber hinaus führten die Parameter nach DIN 4220 und die mithilfe von Pedotransferfunktionen aus Siebanalysen genommenen Parameter zu unterschiedlichen Ergebnissen. Im Rahmen des am Institut für Abfallwirtschaft und Altlasten durchgeführten Projektes EGSIM wurden die Programme SENSIT und ISSOP (in Zusammenarbeit mit DUALIS GmbH IT Solution) entwickelt und zur Parameteridentifikation/-kalibrierung benutzt. Die im Rahmen dieser Arbeit erzielten Ergebnisse konnten nachweisen, unter welchen Bedingungen eine Nachklärung des vollbiologisch gereinigten Abwassers innerhalb der Bodenzone möglich ist, so dass am Ort der Beurteilung (Grundwasseroberfläche) kein unzulässiger Schadstoffeintrag erfolgt. In Bezug auf die KKA-Feldmodelle ist die Anwendung des rotationssymmetrischen Koordinatensystems als Voraussetzung der Realität besser zu entsprechen und nicht als Option zu betrachten. Darüber hinaus wurden anhand der Feldmodelle verschiedene Szenarien mit kontinuierlicher und diskontinuierlicher Versickerung sowie zwei unterschiedlich großen Einleitflächen durchgeführt. Das Programm PCSiWaPro® ist sowohl für ungesättigte als auch für variabel-gesättigte porösen Medien anwendbar. Dies wurde im Rahmen der Simulation des Wasserhaushaltes eines Erddamms nachgewiesen. Die durchschnittliche relative Abweichung zwischen gemessenen und mittels PCSiWaPro® berechneten Wasserständen des entsprechenden Beobachtungspunkts im untersuchten Dammkörper lag bei 0,08 % (entspricht 5,8 cm bei einer Müchtigkeit von ca. 70 m) und das Bestimmtheitsmaß (R2) betrug 0,987. Die Simulation des Wasserhaushaltes unterschiedlicher Deponieabdichtungssystemen mittels PCSiWaPro® zeigte im Allgemeinen ein funktionierendes Ableiten des auf Deponien anfallenden Regenwassers (auch bei Starkregenereignissen). Darüber hinaus haben die durchgeführten Bewuchs-Modelle nachweisen können, dass die Vegetation der Deponieoberflächen den Wassergehalt, durch Pflanzenwurzelentzug bzw. Evapotranspiration, reduzieren können. Die Simulationsergebnisse der durchgeführten Szenarien des Wasserhaushaltes von Lysimetern zur prognostischen Berechnung der Grundwasserneubildung mittels PCSiWaPro® konnten nachweisen, dass das Programm für die Berechnung der Grundwasserneubildungsrate für diesen Zweck anwendbar ist. Die relativen Abweichungen der be-rechneten von den gemessenen Grundwasserneubildungsraten sind auf die verwendeten Materialparameter sowie auf Vernachlässigung der möglicherweise in Lysimetern sich befin-denden Makroporen (duale Porosität) zurückzuführen. / The modelling of complex systems such as the underground is a means to describe the processes occurring in the reality. The conducting of experiments on a model to obtain qualitative evidence about a real system is referred to as a simulation. Thereby, various models (e.g. physical and mathematical models) can be used. The unsaturated zone (vadose zone) is the region between the land surface and the water table, in which the water content is less than full saturation, and the pressure is lower than the atmospheric pressure. The unsaturated zone is very significant for agriculture, geobiology, aerobic degradation processes and groundwater recharge. The processes of water flow and solute transport in the unsaturated zone can be described by means of numerical simulation programs. The aim of the present work is a comprehensive validation of the simulation program PCSiWaPro® (developed at the TU-Dresden, Institute of Waste Management and Contaminated Site Treatment) for different applications. Another aim of this work is to investigate the applicability of the current version of PCSiWaPro® for different cases of a combination between the unsaturated zone and technical facilities. Four application cases with different objectives were investigated within the present work, which are: the simulation of decentralized wastewater infiltration with corresponding column and field experiments, the computation of groundwater recharge by means of lysimeters, the water balance of earth dams and the modelling of landfill covering systems. The application cases differ from each other by the objective of the simulation, the geometry, the size, the specified initial and boundary conditions, the simulation time, the applied materials, the coordinate system, the input and output data. The simulation results clearly showed that PCSiWaPro® is applicable for all investigated cases under consideration of different flow and solute transport regimes, parameters, boundary conditions, spatial and temporal discretization, and coordinate systems. The simulation results of the experimental soil columns for the decentralized treated wastewater infiltration case showed a very good agreement between measured and computed values of water and solute balance (pressure head, flow and solute concentration) of the investigated soil types B3 (slightly silty sand), B4 (coarse sand / gravel) and B5 (medium silty sand). The root of the mean squared error (RMSE) for the computation of the pressure head was 1,84 cm at B5, 3,61 cm at B3 and 1,27 cm at B4. The relative deviation in case of pressure head computation was 2,19 % at B5, 1,3 % at B3 and 5,3 % at B4. The implementation of the sensitivity analysis of the relevant parameters for the modelling showed a very high sensitivity of the VAN GENUCHTEN parameters and the saturated hydraulic conductivity of the soil. Moreover, the parameters according to DIN 4220 led to different results than the estimated ones according to pedotransfer methods based on sieve analysis. Within the project EGSIM, which was carried out at the Institute for waste management and contaminated sites treatment in collaboration with DUALIS GmbH IT Solution, the programs SENSIT and ISSOP were developed and used for parameter identification/ calibration. The results obtained in this Work showed under which conditions is a secondary treatment of full biologically treated wastewater in the soil possible, so that no unallowable pollutants entry in the groundwater occurs. With regard to the field models of this application the implementation of the rotationally symmetric coordinate system should be considered as a condition and not as an option for a better corresponding to the reality. Furthermore, different scenarios of the field models were carried out with continuous and discontinuous infiltration, as well as under different initiation areas. PCSiWaPro® could be applied for both unsaturated and variably-saturated porous media. This could be proven by the simulation of the water balance in an earth dam. The average relative deviation between measured and simulated water levels of the corresponding observation point in the investigated dam embankment was 0,08 % (corresponding to 5,8 cm at 70 m thickness) and the coefficient of determination (R2) was 0,987. In general, the simulation of the water balance using PCSiWaPro® of different landfill covering systems showed a successful draining of the falling rainwater (even under heavy rainfall). In addition, the implemented vegetation models have proven that the vegetation of the landfill surface can reduce the water content in the landfill by evapotranspiration and water uptake by roots. The water balance simulation results of the scenarios for the computation of groundwater recharge by means of lysimeters showed that the program is applicable for this case. The relative deviation of the simulated from the measured groundwater recharge rates occur due to the implemented material parameters as well as to the neglect of macro pores effects (dual porosity).

A model is developed to describe solute transport and retention in fractured rocks. It accounts for the fact that solutes not only can diffuse directly from the flowing channel into the adjacent rock matrix composed of different geological layers but can also at first diffuse into the stagnant water zone occupied in part of the fracture and then from there into the rock matrix adjacent to it. Moreover, the effect of radioactive decay-chain has also been studied in the presence of matrix comprising different geological layers. In spite of the complexities of the system, the analytical solution obtained for the Laplace-transformed concentration at the outlet of the flowing channel can conveniently be transformed back to the time domainby use of e.g. De Hoog algorithm. This allows one to readily include it into a fracture network modelorachannelnetwork model to predictnuclide transport through channels in heterogeneous fracturedmedia consisting of an arbitrary number of rock units withpiecewise constant properties. Simulations made in this study indicate that, in addition to the intact wall rock adjacent to the flowing channel, the stagnant water zone and the rock matrix adjacent to it may also lead to a considerable retardation of solute in cases with a narrow channel. The results further suggest that it is necessary to account for decay-chain and also rock matrix comprising at least two different geological layers in safety and performance assessment of the repositories for spent nuclear fuel. The altered zone may cause a great decrease of the nuclide concentration at the outlet of the flowing channel. The radionuclide decay, when accounted for, will drastically decrease the concentration of nuclides, while neglecting radioactive ingrowth would underestimate the concentration of daughter nuclides. / <p>QC 20140224</p>

No description available.

No description available.