The application of mathematical models of the transport of chemical substances in the remediation of consequences of the uranium mining
|
|
- Geoffrey Fisher
- 5 years ago
- Views:
Transcription
1 The application of mathematical models of the transport of chemical substances in the remediation of consequences of the uranium mining Michal Beneš, Martin Stýblo Department of Mathematics, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University, Prague Trojanova 13, Prague Czech Republic Phone: Fax: JiříMaryška, JiříMužák Department of Control Engineering, Faculty of Mechatronics, Technical University of Liberec Hálkova 6, Liberec Czech Republic Phone: /268 Fax: Abstract A large contamination of the Cenomanian aquifer in the area of the Stráž deposit in Northern Bohemia due to chemical mining of uranium is one of the most important environmental problems in the Czech Republic. The paper deals with mathematical modelling of transport of dissolved chemical substances in this area. The model consists of computation of porous media flow, convection-diffusion transport of chemicals, and simulation of selected chemical reactions. Several finite element formulations (i.e. primal, mixed-hybrid and mixed) of flow and transport problems are considered. For spatial discretization, trilateral prismatic elements are used to describe the structure of stratified layers. Results of numerical experiments and practical computations including simulation of remediation scenarios in the area of leaching fields and modelling of pollution leakage from waste ponds are discussed. 1
2 2 1 Situation in the Stráž deposit Since 1968, more than 4 million tons of sulphuric acid (H 2 SO 4 ), 300 thousand tons of nitric acid (HNO 3 ), 120 thousand tons of ammonia (NH 3 ) and smaller quantities of other chemicals have been injected into the Stráž deposit. A part of the contaminants (about 0.5 %) is situated in the upper Turonian aquifer, which is a source of drinking water for the surrounding region. This contamination is caused by the escape of solutions either from surface technologies or from damaged casings. The wells of older construction had no doublecasing in the Turonian aquifer. In the Cenomanian aquifer, two different areas can be defined: area of leaching fields and area of solution excursion towards the deep mine. The lower part of the geological profile in the leaching fields area with uranium ores (washout horizon, friable sandstones) is filled with solutions of total salinity g/l. High concentrations and large vertical extents of contamination (in some places the full thickness of the fucoid sandstones, i.e. 40 meters) can be found in the area of older leaching fields where extremely high dosages of sulphuric acid were used. Similar solutions also are in more permeable parts of the profile, in the neighbourhood of the leaching fields ( meters). In the area of the solution escape, the thickness of contaminated horizon is decreasing from the hydraulic barrier to the drainage system of the deep mine. The purpose of this drainage system is to hold back acid water contaminated from the leaching process. The salinity of solutions varies from 1 to 10 g/l; ph is usually from 1.8 to 3.5. In the whole Cenomanian aquifer, there are more then 180 millions m 3 of contaminated water. For details, see [1]. 2 Mathematical models of transport The most general mathematical model under consideration consists of porous-media flow problem, transport of chemical substances with diffusion and dispersion, and thermodynamical modelling of chemical processes. Porous-media flow problem is described by the Darcy s law and the continuity equation u = A p ;.u = q in Ω, (1) where p denotes piezometric head, u is filtration velocity, A is symmetric positive definite second rank tensor of hydraulic permeability of the porous medium, q represents density of potential liquid sources in the medium. The boundary conditions are given by p = p D on Ω 1, n.(a p) =0on Ω 2 (impermeable part of the boundary) and n.(a p)+σp = σp D on Ω 3, where n is unit vector of outer normal (see [2]). The transport of chemical species dissolved in the technological solution is described by the following law l i = D(u) l c ; [D] ij = D m δ ij + α T. u.δ ij +(α L α T ) u iu j, (2) u where l i is the flux of l-th substance, l c denotes its concentration, and D is the diffusivity-dispersivity tensor. D m denotes the coefficient of molecular diffusion, α T and α L the transversal and the longitudinal dispersivity, respectively. The equation of mass balance for l-th substance is l c t + u c+ li+ l l cq l r( l c, i c,...) = l c q + ; l=1,..., L. (3) where the first term describes storage of the l-th substance, second term describes convection, and third term reflects the influence of diffusion and dispersion. The term l cq describes the influence of sink of the l-th species and the term l c q + determines source of l-th species. Last term on the left hand side equation, l r( l c, i c,...), describes chemical processes. The boundary conditions are l c = 0 on the part of boundary Ω +, where u.n 0and l i= 0 on the remaining part Ω (see [3]). The reaction term l r( l c, i c,...) describes the influence of all reactions where the l-th species is present. In general, it can be rewritten as l r( l c, i c,...) =,l r( l c, i c,...)+ +,l r( l c, i c,...) (4)
3 3 where,l r expresses the decrease of the l-th species by its reaction with other species, and,l r expresses the increase of the l-th species by reactions of other species (see [4]). The chemical situation in the mining field is very complex. Therefore, the diffusion-convection problem is solved on some simple testing examples whereas the real situation is simulated as follows: the stationary porous-media flow problem is solved; the convective transport over finite volumes is computed; the diffusion and dispersion is computed; the chemical situation is computed on each volume separately. The following types of reactions are considered: fast reactions in the technological solution described by the Guldberg-Waag law schematically written as: k G( i c)=0 (5) Such equilibrium relations represent a coupling (constraint) which should be considered when solving dynamical relations as shown below. kinetic reactions with the rock (dissolving); kinetic reactions of precipitation; sorption on the surface of the rock. Kinetic relations describing previous three cases can be written in the following form ([5]): d l c dt = ν l,i k i,l ( l c) ν l,i ( i c) ν i,l i I l + ν(i,j) l k ij( i c) ν ij ( j c) ν ji. (6) i,j l Here, I l denotes the set of indeces of species reacting with the l-th species, ν l,i is stechiometric factor of l-th species if reacting with the i-th species, and k i,l is the reaction rate. The coefficient νl,i is the modified stechiometric factor. 3 Mathematical models and numerical methods Considering several different requirements of the application, the following mathematical problems are solved: Steady-state porous media flow with free boundary is modelled by primal FEM. By an iteration process, the grid is deformed according to the phreatic surface (see [2]). For transport modelling, the mixed-hybrid FEM model of porous media is computed on a grid obtained in the previous step in order to get a conservative approximation of the flow field (see [6]). The system of convection-reaction equations is solved by an explicit finite-volume upwind method using the flow field computed in the previous step (see [7]). The chemical kinetics is computed on each element separately. The discretization of the problem is given by physical and chemical characteristics of the geological region. There are used trilateral prismatic elements with vertical faces and general nonparallel basis (cnf. Figure 1). The five-face elements and the simplices are used to model closures of the sedimented layers. These elements describe character of stratified sedimented layers and complicated horizontal profile with many wells. In case of real problems, it is necessary to consider large area of several km 2. The thickness of cenomanian sandstones is several tens of meters. Consequently, real elements have considerable size disproportions in horizontal and vertical nodal distances. The influence of these disproportions on the spectrum of the resulting system of linear equations and on the rate of convergence of numerical methods is studied.
4 4 x 6 x 4 n e 2 x 5 n e n e n e x 2 z x 1 n e 1 x 3 y x Figure 1: The RT 0 1,h and H1 h elements. x 4 x 5 t c 2 x 12 6 c23 c 13 x 1 t 1 x 3 x 2 The following approximation spaces are used: The Raviart-Thomas space RT 1,h 0 of linear vector valued functions for mixed-hybrid formulation. The finite space Hh 1 for primal and mixed formulation (6 nodal degrees of freedom). The Raviart-Thomas space of vector valued functions was generated on the following type of elements for mixed-hybrid formulation vj e := kj e [ x1 αj1 e, x 2 αj2, e βjx e 3 αj3 e ]. Parameters of these functions are determined by the conditions n e j vi e ds := δ ij, i,j =1,..., 5 f e j where fj e denotes j-th face of the element e n e j is its unit normal (with respect to the element e ). The Raviart-Thomas approximation of the porous-media flow problem yields the system of linear algebraic equations with symmetric indefinite matrix. The finite space Hh,ψ 1 (Ω) was generated for primal and mixed formulation. Any element is decomposed by intersections of diagonals of vertical faces and centres of gravity of bases into 11 subelements (simplices). On these simplices, there are defined linear functions in such a way, that base function on whole element is continous and piece-wise linear. The approximation of primal formulation leads to a system of linear algebraic equations with sparse positively definite matrix (for details, see [8], [9]). 4 Modelling of the remediation process In 1996, the uranium mining has been stopped, and the major interest of the Czech government is in the removal of consequences of former mining activity - in the optimal remediation scenario. For this purpose, the model is used to produce short- and long-term predictions of behaviour of whole area under different outer conditions. The transport equation is discretized in time with timestep that corresponds to the size of time interval in which the whole problem is solved. For the space discretization of the transport problem the same mesh and the same spaces of scalar and vector valued functions are used. The solution of chemical species transport problem in the convective field of filtration velocities leads to the large nonsymmetric sparse system of equations. The transport equation is solved on each timelevel by superimposing the influence of convection solved in the first stage and the influence of diffusion and dispersion. Obtained results of concentration distribution create an input for thermodynamical model of chemical processes. This procedure exhibits certain time stability, which is very important for calculations in long-time intervals (several decades). One of such predictions is shown in the Figure 2. In the first phase, the problem of pumping-out concentrated solutions is solved. Evolution of pressure in the Cenomanian aquifer is solved by a local regional hydraulic model (see [1]), which is a two-layer model with grid covering an area about 300 km 2. This area is mostly limited by tectonic lines. The horizont is splitted vertically into 9-13 layers. The grid consists of 1400 plane elements and spatial elements. Filtration parameters are determined from hydrogeological by experiment, and their
5 Figure 2: The situation of the water level in the Cenomanian aquifer. The situation after starting the evaporation station (in 1996) with less injection into the barrier, and finally the situation after flooding the mine. 5
6 6 vertical distribution is given by data of the Geotechnological Information System (GTIS). The initial conditions (the material concentrations in solution) in the leaching field area are given by monitoring wells. A vertical distribution of the contaminants is controlled by induction logging and the results are stored in the GTIS. Chemical processes are simulated by two different methods. A decay of caolinit by the sulphuric acid, which causes decrease of salinity and increase of aluminium contents in solution, is solved by the thermodynamical kinetic model. This model was developed from the balanced thermodynamical model for solutions and rock of the Stráž deposit. The thermodynamical system of the balanced model consists of 18 main components of solution, 32 minerals, 110 chemical reactions between the solution components and 44 reactions between minerals and solution. The main components of the solution may occur in 184 different forms. In the first stage, the kinetic model was reduced to 5 components and 11 equations. The model solves a leaching of Al from the rock and its precipitation in the form of alunit or aluminium hydroxid. The concentration of H 2 SO 4 and the total amount of sulphur ion varies at the same time. The kinetics of uranium leaching is solved by kinetic functions given by the experiment. The kinetics is simplified in comparison with the detailed models for uranium mining optimization, because it is required to simulate the situation in the whole volume for a period of several decades with reasonable time steps (see [10]). 5 Modelling of the influences of the Stráž wa-ste pond on the groundwater The Stráž wastepondisabout180ha (1.8 km 2 ). It is divided into two parts: the 1st and the 2nd stage. The first stage was built in It is filled-up to height of about 25 m above the former surface. In its upper part, there is a water surface, in which waste muds are floated. The second stage was built in 1988, when the future development of uranium production was planned. After a short term use of this second stage, production was significantly limited. At present, this part of the waste pond forms a lake with an area of about 90 ha (0.9 km 2 ). The lake contains water contaminated by waste from the Stráž uranium mill. The water surface in this lake is approximately on the same level as groundwater. The main contamination is 3 to 5 g/l of SO 4 ions. Underlying rocks of the pit are created by permeable coniacian sandstones with a filtration coefficient k f = m/s. The boundary conditions are determined by the level of water surface in the pit and in the Ploučnice river in the South (Dirichlet boundery conditions). The part of the boundary formed by impermeable rocks is described by homogenous Neumann boundary condition. In addition, the rain dotation during the season, and the dotation from surrounding hills in the North are considered (Newton boundary condition). Details can be found in [1]. 6 Conclusion The models of transport of contaminants in the area of the Stráž deposit are very important for monitoring and prediction of the environmental situation and help in decision of an appropriate remediation scenario for the mine. The real situation is simulated by a multi-stage model including porous-media flow, convective and diffusive transport, and chemical reactions. Such models will be improved by solving some of these stages simultaneously. As the description of situation near waste ponds and in subsurface layer is required, the models will be adapted for this purpose, too. Acknowledgement. 205/96/0921. The reserch is supported by the Czech Grant Agency under the project No. References [1] P. Mareček, J. Novák, and V. Wasserbauer. Mathematical modelling of the influences of the Stráž mill waste pond on the groundwater. In IAEA Technical Committee Meeting on Computer Application in Uranium Exploration and Production Case History and Current Status, Vienna, November 1994.
7 Figure 3: The time development of the concentrations of TDS (total dissolved solids) near the bottom of the friable sandstones before beginning evaporation in 1995, and a cross-sections through the deposit. 7
8 8 [2] J. Maryška. Mathematical models of the underground water flow and chemical transport problem. PhD thesis, FNSPE CTU Prague, [3] J. Maryška and J. Mužák. Mathematical modelling of the transport of chemical species in the contaminated underground water. Short communication at Congres milieux poreux, St. Etienne, [4] M. Stýblo, J. Maryška, and J. Mužák. Mathematical modelling of groundwater flow and transport of chemicals. In Proceedings of the 9th Summer School Software and Algorithms of Numerical Mathematics, pages , Železná Ruda, [5] R. Brdička, M. Kalousek, A. Schütz. Introduction into Physical Chemistry. SNTL, Prague [6] J. Maryška and J. Mužák. Hybrid mixed model of the transport of the chemical substances. In Proceedings of Numerical Modelling in Continuum Mechanics, volume Part II., pages , [7] J. Maryška, J. Mužák, and M. Stýblo. Mixed-hybrid formulation of the transport of chemicals in the contaminated underground water. In S. Míka, M. Míková, and Brandner M., editors, Lecture Notes of IMAMM 95, pages , Srní, June Univ. of West Bohemia in Plzeň. [8] J. Maryška, M. Rozložník,andM.Tůma. Mixed-hybrid finite element approximation of the potential fluid problem. J. Comput. Appl. Math., 63: , [9] J. Maryška, M. Rozložník, andm. Tůma. The potential fluid flow problem and the convergence rate of the minimal residual method. Submitted to Num. Lin. Alg. Appl., [10] M. Stýblo, J. Maryška, J. Mužák, J. Skokan, M. Tůma, J. Drkošková, M. Rozložník, Z. Strakoš, J. Novák, P. Mareček, and V. Wasserbauer. Mathematical modelling of porous media flow and transport of chemicals in underground water. In WORKSHOP 96, pages Czech Technical University in Prague and Technical University in Brno, 1996.
NUMERICAL SOLUTION OF TWO-REGION ADVECTION-DISPERSION TRANSPORT AND COMPARISON WITH ANALYTICAL SOLUTION ON EXAMPLE PROBLEMS
Proceedings of ALGORITMY 2002 Conference on Scientific Computing, pp. 130 137 NUMERICAL SOLUTION OF TWO-REGION ADVECTION-DISPERSION TRANSPORT AND COMPARISON WITH ANALYTICAL SOLUTION ON EXAMPLE PROBLEMS
More informationMODEL OF GROUNDWATER FLOW IN FRACTURED ENVIRONMENT
Proceedings of ALGORITMY 2002 Conference on Scientific Computing, pp. 138 145 MODEL OF GROUNDWATER FLOW IN FRACTURED ENVIRONMENT JIŘÍ MARYŠKA, OTTO SEVERÝN AND MARTIN VOHRALÍK Abstract. A stochastic discrete
More informationEVALUATION OF CRITICAL FRACTURE SKIN POROSITY FOR CONTAMINANT MIGRATION IN FRACTURED FORMATIONS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationParallel Implementation of BDDC for Mixed-Hybrid Formulation of Flow in Porous Media
Parallel Implementation of BDDC for Mixed-Hybrid Formulation of Flow in Porous Media Jakub Šístek1 joint work with Jan Březina 2 and Bedřich Sousedík 3 1 Institute of Mathematics of the AS CR Nečas Center
More informationSafety assessment for disposal of hazardous waste in abandoned underground mines
Safety assessment for disposal of hazardous waste in abandoned underground mines A. Peratta & V. Popov Wessex Institute of Technology, Southampton, UK Abstract Disposal of hazardous chemical waste in abandoned
More informationEXAMPLE PROBLEMS. 1. Example 1 - Column Infiltration
EXAMPLE PROBLEMS The module UNSATCHEM is developed from the variably saturated solute transport model HYDRUS-1D [Šimůnek et al., 1997], and thus the water flow and solute transport parts of the model have
More informationRATE OF FLUID FLOW THROUGH POROUS MEDIA
RATE OF FLUID FLOW THROUGH POROUS MEDIA Submitted by Xu Ming Xin Kiong Min Yi Kimberly Yip Juen Chen Nicole A project presented to the Singapore Mathematical Society Essay Competition 2013 1 Abstract Fluid
More informationSaltwater injection into a fractured aquifer: A density-coupled mass-transport model
Saltwater injection into a fractured aquifer: A density-coupled mass-transport model Junfeng Luo 1, Martina aus der Beek 2, Joachim Plümacher 2, Sven Seifert 1, Bertram Monninkhoff 1 1 DHI-WASY GmbH, Volmerstr.
More informationThree-dimensional Modelling of Reactive Solutes Transport in Porous Media
151 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 41, 214 Guest Editors: Simonetta Palmas, Michele Mascia, Annalisa Vacca Copyright 214, AIDIC Servizi S.r.l., ISBN 978-88-9568-32-7; ISSN 2283-9216
More informationKarst Topography In order to understand karst topography we must first have a basic understanding of the water cycle, the formation of limestone (carb
Karst Topography The formation of caves and other associated features in limestone bedrock is called karst topography. Limestone, a sedimentary rock made mostly of the mineral calcite and small amounts
More informationRADIONUCLIDE DIFFUSION IN GEOLOGICAL MEDIA
GEOPHYSICS RADIONUCLIDE DIFFUSION IN GEOLOGICAL MEDIA C. BUCUR 1, M. OLTEANU 1, M. PAVELESCU 2 1 Institute for Nuclear Research, Pitesti, Romania, crina.bucur@scn.ro 2 Academy of Scientists Bucharest,
More informationNumerical simulation of fracture flow with a mixedhybrid FEM stochastic discrete fracture network model
Numerical simulation of fracture flow with a mixedhybrid FEM stochastic discrete fracture network model Jiří Maryška (jiri.maryska@vslib.cz), Otto Severýn (otto.severyn@vslib.cz) and Martin Vohralík, (vohralik@km1.fjfi.cvut.cz)
More informationLithosphere: (Rocky Sphere) Solid, rocky, outer layer of the Earth. Includes the crust and part of the upper mantle. Lithosphere
Lithosphere: (Rocky Sphere) Solid, rocky, outer layer of the Earth. Includes the crust and part of the upper mantle. Lithosphere Permafrost Permafrost Ground that is at a temperature of 0 or below for
More informationSolving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions
International Workshop on MeshFree Methods 3 1 Solving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions Leopold Vrankar (1), Goran Turk () and Franc Runovc (3) Abstract:
More informationChapter 6 The lithosphere and the hydrosphere
Chapter 6 The lithosphere and the hydrosphere The lithosphere shell of the Earth, consists of crust and upper mantle contains minerals and rocks 1. Minerals solid, inorganic substances with defined composition
More informationEARTH SCIENCE 11 CHAPTER 9 NOTES WATER MOVING UNDERGROUND
EARTH SCIENCE 11 CHAPTER 9 NOTES WATER MOVING UNDERGROUND KEY A. ALL THE WORLD S WATER 1. The world has an estimated one and one third thousand million cubic kilometers of water. 2. The earth s water is
More informationGroundwater. (x 1000 km 3 /y) Oceans Cover >70% of Surface. Groundwater and the. Hydrologic Cycle
Chapter 17 Oceans Cover >70% of Surface Groundwater and the Hydrologic Cycle Vasey s Paradise, GCNP Oceans are only 0.025% of Mass Groundwater Groundwater is liquid water that lies in the subsurface in
More information11/22/2010. Groundwater in Unconsolidated Deposits. Alluvial (fluvial) deposits. - consist of gravel, sand, silt and clay
Groundwater in Unconsolidated Deposits Alluvial (fluvial) deposits - consist of gravel, sand, silt and clay - laid down by physical processes in rivers and flood plains - major sources for water supplies
More informationChemical Hydrogeology
Physical hydrogeology: study of movement and occurrence of groundwater Chemical hydrogeology: study of chemical constituents in groundwater Chemical Hydrogeology Relevant courses General geochemistry [Donahoe]
More informationSoils, Hydrogeology, and Aquifer Properties. Philip B. Bedient 2006 Rice University
Soils, Hydrogeology, and Aquifer Properties Philip B. Bedient 2006 Rice University Charbeneau, 2000. Basin Hydrologic Cycle Global Water Supply Distribution 3% of earth s water is fresh - 97% oceans 1%
More informationEssentials of Geology, 11e
Essentials of Geology, 11e Groundwater Chapter 10 Instructor Jennifer Barson Spokane Falls Community College Geology 101 Stanley Hatfield Southwestern Illinois Co Jennifer Cole Northeastern University
More informationA posteriori error estimates, stopping criteria, and adaptivity for multiphase compositional Darcy flows in porous media
A posteriori error estimates, stopping criteria, and adaptivity for multiphase compositional Darcy flows in porous media D. A. Di Pietro, M. Vohraĺık, and S. Yousef Université Montpellier 2 Marseille,
More informationNumerical Solution of the Two-Dimensional Time-Dependent Transport Equation. Khaled Ismail Hamza 1 EXTENDED ABSTRACT
Second International Conference on Saltwater Intrusion and Coastal Aquifers Monitoring, Modeling, and Management. Mérida, México, March 3-April 2 Numerical Solution of the Two-Dimensional Time-Dependent
More informationSection 1: Earth s Interior and Plate Tectonics Section 2: Earthquakes and Volcanoes Section 3: Minerals and Rocks Section 4: Weathering and Erosion
Section 1: Earth s Interior and Plate Tectonics Section 2: Earthquakes and Volcanoes Section 3: Minerals and Rocks Section 4: Weathering and Erosion Key Terms Crust Mantle Core Lithosphere Plate Tectonics
More informationGroundwater. (x 1000 km 3 /y) Reservoirs. Oceans Cover >70% of Surface. Groundwater and the. Hydrologic Cycle
Chapter 13 Oceans Cover >70% of Surface Groundwater and the Hydrologic Cycle Oceans are only 0.025% of Mass Groundwater Groundwater is liquid water that lies in the subsurface in fractures in rocks and
More informationWhat is a water table? What is an aquifer? What is the difference between a spring and a well?
CHAPTER 11 3 Water Underground SECTION The Flow of Fresh Water BEFORE YOU READ After you read this section, you should be able to answer these questions: What is a water table? What is an aquifer? What
More informationProf. Stephen A. Nelson EENS 111. Groundwater
Page 1 of 8 Prof. Stephen A. Nelson EENS 111 Tulane University Physical Geology This page last updated on 20-Oct-2003 is water that exists in the pore spaces and fractures in rock and sediment beneath
More informationENVIRONMENTAL EFFECTS OF GROUNDWATER WITHDRAWAL IN SOUTH NYÍRSÉG
PhD thesis ENVIRONMENTAL EFFECTS OF GROUNDWATER WITHDRAWAL IN SOUTH NYÍRSÉG János Szanyi Szeged, 2004 ENVIRONMENTAL EFFECTS OF GROUNDWATER WITHDRAWAL IN SOUTH NYÍRSÉG Preliminaries, the aims of the dissertation
More informationMathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts 1 and Stephan Matthai 2 3rd Febr
HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts and Stephan Matthai Mathematics Research Report No. MRR 003{96, Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL
More informationOrigin and Evolution of Formation Waters in the West-Central Part of the Alberta Basin
Page No. 004-1 Origin and Evolution of Formation Waters in the West-Central Part of the Alberta Basin Karsten Michael* University of Alberta, 1-26 ESB, Edmonton, AB T6G 2E3 karsten@ualberta.ca and Stefan
More informationA BOOKLET ON. T Rangasamy, A R Leach and A P Cook. Facilitating safety and health research in the South African mining industry
A BOOKLET ON THE HYDRAULIC DESIGN OF COAL BARRIER PILLARS T Rangasamy, A R Leach and A P Cook Facilitating safety and health research in the South African mining industry A BOOKLET ON THE HYDRAULIC DESIGN
More information1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #5: Groundwater Flow Patterns. Local Flow System. Intermediate Flow System
1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #5: Groundwater Flow Patterns c Local Flow System 10,000 feet Intermediate Flow System Regional Flow System 20,000 feet Hydrologic section
More informationNinth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673
Ninth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673 A NEW NUMERICAL APPROACH FOR THE SOLUTION OF CONTAMINANT TRANSPORT EQUATION Mohamed El-Gamel Department of Mathematical
More informationDNAPL migration through interbedded clay-sand sequences
Groundwater Quality: Natural and Enhanced Restoration of Groundwater Pollution (Proceedings ofthe Groundwater Quality 2001 Conference held al Sheffield. UK. June 2001). IAHS Publ. no. 275. 2002. 455 DNAPL
More informationImpact of the Danube River on the groundwater dynamics in the Kozloduy Lowland
GEOLOGICA BALCANICA, 46 (2), Sofia, Nov. 2017, pp. 33 39. Impact of the Danube River on the groundwater dynamics in the Kozloduy Lowland Peter Gerginov Geological Institute, Bulgarian Academy of Sciences,
More informationWATER ON AND UNDER GROUND. Objectives. The Hydrologic Cycle
WATER ON AND UNDER GROUND Objectives Define and describe the hydrologic cycle. Identify the basic characteristics of streams. Define drainage basin. Describe how floods occur and what factors may make
More informationThermohaline Flow and Reactive Solute Transport in Fractured Rock
Thermohaline Flow and Reactive Solute Transport in Fractured Rock Thomas GRAF Institute of Fluid Mechanics in Civil Engineering Alberta Geological Survey, Edmonton 2008-2009 Basin Hydrogeologist Sasdf
More informationBuilding ground level
TMA4195 MATHEMATICAL MODELLING PROJECT 212: AQUIFER THERMAL ENERGY STORAGE 1. Introduction In the project we will study a so-called Aquifer Thermal Energy Storage (ATES) system with the aim of climitizing
More informationLand subsidence due to groundwater withdrawal in Hanoi, Vietnam
Land Subsidence (Proceedings of the Fifth International Symposium on Land Subsidence, The Hague, October 1995). 1AHS Publ. no. 234, 1995. 55 Land subsidence due to groundwater withdrawal in Hanoi, Vietnam
More informationModeling Reactive Flows in Porous Media
Modeling Reactive Flows in Porous Media Peter Lichtner (lead PI), Los Alamos National Laboratory Glenn Hammond, Pacific Northwest National Laboratory Richard Tran Mills, Oak Ridge National Laboratory NCCS
More information1. Resistivity of rocks
RESISTIVITY 1) Resistivity of rocks 2) General principles of resistivity surveying 3) Field procedures, interpretation and examples 4) Summary and conclusions INDUCED POLARIZATION 1) General principles
More informationAngel International School - Manipay 1 st Term Examination November, 2015
Grade 08 Angel International School - Manipay 1 st Term Examination November, 2015 Chemistry Duration: 2 Hours Index No:- Underline the correct answer. 1. Which of the following is a fossil fuel? a) Natural
More informationABSTRACT GOVERNING EQUATIONS
A three dimensional finite element model for fluid flow and transport in confined or unconfined aquifer F. Jacob*, J.M. Crolef, P. Lesaint*, J. Mania* "Laboratoire de calcul scientifique, ^Laboratoire
More informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationGroundwater Rebound in the South Yorkshire Coalfield: A review of initial modelling
Groundwater Rebound in the South Yorkshire Coalfield: A review of initial modelling 1 S. P. Burke, 1 H. A.B. Potter and 2 A. Jarvis 1 Environment Agency: Science Group, Olton Court, Olton,Solihull UK 2
More information1. INTRODUCTION TO CFD SPRING 2019
1. INTRODUCTION TO CFD SPRING 2019 1.1 What is computational fluid dynamics? 1.2 Basic principles of CFD 1.3 Stages in a CFD simulation 1.4 Fluid-flow equations 1.5 The main discretisation methods Appendices
More informationColloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.
Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
More informationBefore writing the discipline examinations, candidates must have passed, or have been exempted from, the Basic Studies Examinations.
INTRODUCTION The Canadian Engineering Qualifications Board of Engineers Canada issues the Examination Syllabus that includes a continually increasing number of engineering disciplines. Each discipline
More informationChapter 2 Theory. 2.1 Continuum Mechanics of Porous Media Porous Medium Model
Chapter 2 Theory In this chapter we briefly glance at basic concepts of porous medium theory (Sect. 2.1.1) and thermal processes of multiphase media (Sect. 2.1.2). We will study the mathematical description
More informationKRONDORF MINERAL WATER
Water character Project Krondorf the restoration of traditional production of mineral water that preserves the source as well as the environment, utilising centuries-old procedures and focusing on water
More information' International Institute for Land Reclamation and Improvement. 2 Groundwater Investigations. N.A. de Ridder'? 2.1 Introduction. 2.
2 Groundwater Investigations N.A. de Ridder'? 2.1 Introduction Successful drainage depends largely on a proper diagnosis of the causes of the excess water. For this diagnosis, one must consider: climate,
More informationAlgebraic flux correction and its application to convection-dominated flow. Matthias Möller
Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany
More information1978 MSc Geology, South Dakota School of Mines and Technology Introduction to Groundwater Vistas (Groundwater Modeling) Albuquerque, New Mexico
Clyde Yancey, BSc, MSc, PGeo Vice President of Exploration Education 1978 MSc Geology, South Dakota School of Mines and Technology 1975 BA Geology Trinity University 1999 Introduction to Groundwater Vistas
More informationRT3D Rate-Limited Sorption Reaction
GMS TUTORIALS RT3D Rate-Limited Sorption Reaction This tutorial illustrates the steps involved in using GMS and RT3D to model sorption reactions under mass-transfer limited conditions. The flow model used
More informationEarth s Structure and Natural Processes Practice Test
Name: Earth s Structure and Natural Processes Practice Test Section: Directions: For each of the questions or incomplete statements below, choose the best of the answer choices given and write your answer
More informationSTUDY GUIDE FOR CONTENT MASTERY. Movement and Storage of Groundwater
Groundwater SECTION 10.1 Movement and Storage of Groundwater In your textbook, read about the hydrosphere, precipitation and groundwater, and groundwater storage. Use the following terms to complete the
More informationIn Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationSurface Processes Focus on Mass Wasting (Chapter 10)
Surface Processes Focus on Mass Wasting (Chapter 10) 1. What is the distinction between weathering, mass wasting, and erosion? 2. What is the controlling force in mass wasting? What force provides resistance?
More informationGEOCHEMISTRY, GROUNDWATER AND POLLUTION,
GEOCHEMISTRY, GROUNDWATER AND POLLUTION, 2 ND EDITION C.A.J. APPELO Hydrochemical Consultant, Amsterdam, the Netherlands D. POSTMA Environment & Resources DTU, Technical University of Denmark, Kgs. Lyngby,
More informationLecture Note for Open Channel Hydraulics
Chapter -one Introduction to Open Channel Hydraulics 1.1 Definitions Simply stated, Open channel flow is a flow of liquid in a conduit with free space. Open channel flow is particularly applied to understand
More informationTEACHER BACKGROUND KNOWLEDGE. Surface Processes: Weathering and Erosion
TEACHER BACKGROUND KNOWLEDGE Surface Processes: Weathering and Erosion Core Concepts Precipitation (caused by the water cycle) and wind cause rocks to be broken into smaller pieces in the process called
More information1. INTRODUCTION TO CFD SPRING 2018
1. INTRODUCTION TO CFD SPRING 018 1.1 What is computational fluid dynamics? 1. Basic principles of CFD 1.3 Stages in a CFD simulation 1.4 Fluid-flow equations 1.5 The main discretisation methods Appendices
More informationDeep Borehole Disposal Performance Assessment and Criteria for Site Selection
Deep Borehole Disposal Performance Assessment and Criteria for Site Selection Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department
More informationCode-to-Code Benchmarking of the PORFLOW and GoldSim Contaminant Transport Models using a Simple 1-D Domain
Code-to-Code Benchmarking of the PORFLOW and GoldSim Contaminant Transport Models using a Simple 1-D Domain - 11191 Robert A. Hiergesell and Glenn A. Taylor Savannah River National Laboratory SRNS Bldg.
More informationEOC Study Guide Honors
Name Date Science. SC.912.E.5.1 Cite evidence used to develop and verify the scientific theory of the Big Bang (also known as the Big Bang Theory) of the origin of the universe. 1. What is the approximate
More informationGroundwater Modeling for Flow Systems with Complex Geological and Hydrogeological Conditions
Available online at www.sciencedirect.com Procedia Earth and Planetary Science 3 ( 2011 ) 23 28 2011 Xi an International Conference on Fine Geological Exploration and Groundwater & Gas Hazards Control
More informationUranium Mining Activities in Los Gigantes, Argentina: Possible Case Study Site
Uranium Mining Activities in Los Gigantes, Argentina: Possible Case Study Site Dr. Daniel S. Cicerone Soil and Water Chemistry Division Department of Environment, Chemistry Management National Atomic Energy
More informationMATHEMATICAL MODELING AND NUMERICAL SOLUTION OF 3D ATMOSPHERIC BOUNDARY LAYER
, Vol, Pt, Special Issue Proceedings of International Conference RDAMM 585 MATHEMATICAL MODELING AND NUMERICAL SOLUTION OF D ATMOSPHERIC BOUNDARY LAYER L. Beneš, K. Kozel Department of Technical Mathematics,
More informationUnderground Storage & Disposal - The Salt Concept. Thomas Brasser - GRS
Underground Storage & Disposal - The Salt Concept Thomas Brasser - GRS Who is GRS ( Plant & Reactor Safety Ltd. ) Non-profit, independent expert and research organisation Assess and improve safety of technical
More informationJ. Environ. Res. Develop. Journal of Environmental Research And Development Vol. 8 No. 1, July-September 2013
SENSITIVITY ANALYSIS ON INPUT PARAMETERS OF 1-D GROUNDWATER FLOW GOVERNING EQUATION : SOLVED BY FINITE-DIFFERENCE AND THOMAS ALGORITHM IN MICROSOFT EXCEL Goh E.G.* 1 and Noborio K. 2 1. Department of Engineering
More informationAvailable online at ScienceDirect. Energy Procedia 114 (2017 )
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 114 (2017 ) 2772 2780 13th International Conference on Greenhouse Gas Control Technologies, GHGT-13, 14-18 November 2016, Lausanne,
More informationMonday 2 June 2014 Morning
Monday 2 June 2014 Morning A2 GCE GEOLOGY F794/01 Environmental Geology *3076922876* Candidates answer on the Question Paper. OCR supplied materials: None Other materials required: Electronic calculator
More informationChapter 13. Groundwater
Chapter 13 Groundwater Introduction Groundwater is all subsurface water that completely fills the pores and other open spaces in rocks, sediments, and soil. Groundwater is responsible for forming beautiful
More informationEvaluation of the hydraulic gradient at an island for low-level nuclear waste disposal
A New Focus on Groundwater Seawater Interactions (Proceedings of Symposium HS1001 at IUGG2007, Perugia, July 2007). IAHS Publ. 312, 2007. 237 Evaluation of the hydraulic gradient at an island for low-level
More informationNONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS
June - July, 5 Melbourne, Australia 9 7B- NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS Werner M.J. Lazeroms () Linné FLOW Centre, Department of Mechanics SE-44
More informationThe Outlook on Potential Uranium ISL Mining at Nyota Deposit (Tanzania)
The Outlook on Potential Uranium ISL Mining at Nyota Deposit (Tanzania) А. Boytsov, S. Stander, V. Martynenko URAM 2014, IAEA, Vienna, Austria June 2014 Outline of Presentation 1. Uranium One Introduction
More informationModeling and Measurement of Thermal Process in Experimental Borehole in Matlab&Simulink and Comsol Multiphysics
Modeling and Measurement of Thermal Process in Experimental Borehole in Matlab&Simulink and Comsol Multiphysics STEPAN OZANA, RADOVAN HAJOVSKY, MARTIN PIES, BLANKA FILIPOVA Department of Cybernetics and
More informationADVENTURES IN WATER DEVELOPED BY LOUISVILLE WATER COMPANY
ADVENTURES IN WATER DEVELOPED BY LOUISVILLE WATER COMPANY Tunneling for Water explains the science behind a first-of-its-kind project in the world! Louisville Water Company is the first water utility
More informationOn the validation study devoted to stratified atmospheric flow over an isolated hill
On the validation study devoted to stratified atmospheric flow over an isolated hill Sládek I. 2/, Kozel K. 1/, Jaňour Z. 2/ 1/ U1211, Faculty of Mechanical Engineering, Czech Technical University in Prague.
More informationSite Investigation and Landfill Construction I
Site Investigation and Landfill Construction I Gernot Döberl Vienna University of Technology Institute for Water Quality, Resources and Waste Management Contents Site Investigation Base Liners Base Drainage
More informationMapping Wall Rock Chemistry on the Ultimate Pit Surface of an Open Pit Mine Utilizing GIS
Mapping Wall Rock Chemistry on the Ultimate Pit Surface of an Open Pit Mine Utilizing GIS Patsy Moran, PhD Tetra Tech, Golden, U.S.A. David Richers, PhD Tetra Tech, Golden, U.S.A. Abstract Knowledge of
More informationGeophysical Monitoring Researches for CO 2 Geological Storage. Shinsuke NAKAO Geological Survey of Japan, AIST
Geophysical Monitoring Researches for CO 2 Geological Storage Shinsuke NAKAO Geological Survey of Japan, AIST 1 2 OUTLINE Introduction Multi-lateral Background, Objectives Geophysical Monitoring Researches
More informationDifferentiation of chloride source using stable chlorine isotopes
Differentiation of chloride source using stable chlorine isotopes RemTech 2009 Banff, Alberta October 14-16 2009 Dr. Alec Blyth and Tom Anthony 1 1 currently City of Calgary Background: DOW Chemical Canada
More informationMathematical model of Baltic artesian basin
Mathematical model of Baltic artesian basin Juris Sennikovs, Janis Virbulis, and Uldis Bethers Laboratory for Mathematical Modelling of Environmental and Technological Processes UNIVERSITY OF LATVIA Contents
More informationBasin-scale Modeling of CO 2 Sequestration in the Illinois Basin Status Report
Basin-scale Modeling of CO 2 Sequestration in the Illinois Basin Status Report Edward Mehnert, James Damico, Scott Frailey, Hannes Leetaru, Yu-Feng Lin, Roland Okwen Illinois State Geological Survey, Prairie
More informationTodd Arbogast. Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and Sciences (ICES)
CONSERVATIVE CHARACTERISTIC METHODS FOR LINEAR TRANSPORT PROBLEMS Todd Arbogast Department of Mathematics and, (ICES) The University of Texas at Austin Chieh-Sen (Jason) Huang Department of Applied Mathematics
More information12 10 8 6 4 2 0 40-50 50-60 60-70 70-80 80-90 90-100 Fresh Water What we will cover The Hydrologic Cycle River systems Floods Groundwater Caves and Karst Topography Hot springs Distribution of water in
More informationYear 6. Geography. Revision
Year 6 Geography Revision November 2017 Rivers and World knowledge How the water cycle works and the meaning of the terms evaporation, condensation, precipitation, transpiration, surface run-off, groundwater
More informationSAFETY ASSESSMENT CODES FOR THE NEAR-SURFACE DISPOSAL OF LOW AND INTERMEDIATE-LEVEL RADIOACTIVE WASTE WITH THE COMPARTMENT MODEL: SAGE AND VR-KHNP
SAFETY ASSESSMENT CODES FOR THE NEAR-SURFACE DISPOSAL OF LOW AND INTERMEDIATE-LEVEL RADIOACTIVE WASTE WITH THE COMPARTMENT MODEL: SAGE AND VR-KHNP J. B. Park, J. W. Park, C. L. Kim, M. J. Song Korea Hydro
More informationSustainable Energy Science and Engineering Center GEOTHERMAL ENERGY. Sustainable Energy Sources. Source:
Sustainable Energy Sources GEOTHERMAL ENERGY Earth s Temperature Profile GEOTHERMAL ENERGY Plate Tectonics Earth's crust is broken into huge plates that move apart or push together at about the rate our
More informationKey Stage 3 - Volcano Fracking
After the meeting Come out of your role. Write your own summary of the dilemmas facing the council, and recommend whether or not they should allow fracking to take place. Pupil worksheet is short for hydraulic
More informationPRELIMINARY REPORT. Evaluations of to what extent CO 2 accumulations in the Utsira formations are possible to quantify by seismic by August 1999.
SINTEF Petroleumsforskning AS SINTEF Petroleum Research N-7465 Trondheim, Norway Telephone: +47 73 59 11 Fax: +477359112(aut.) Enterprise no.: NO 936 882 331 MVA TITLE PRELIMINARY REPORT Evaluations of
More informationCapabilities of TOUGH Codes for Modeling Geologic Sequestration and Leakage of CO 2
Capabilities of TOUGH Codes for Modeling Geologic Sequestration and Leakage of CO 2 Karsten Pruess Earth Sciences Division Lawrence Berkeley National Laboratory Presented at Workshop on Leakage Modeling
More informationILLINOIS GROUNDWATER - BEDROCK
ILLINOIS GROUNDWATER - BEDROCK Northern Illinois Fresh Water Water Supply Shallow Bedrock Aquifers < 500 depth Deep Bedrock Aquifers > 500 depth Southern Illinois Saline water Oil Fracking Coal Waste injection
More informationDEPARTMENT OF GEOSCIENCES
DEPARTMENT OF GEOSCIENCES Office in Natural Resources Building, Room 322 (970) 491-7826 warnercnr.colostate.edu/geosciences-home (http:// warnercnr.colostate.edu/geosciences-home) Richard Aster, Department
More informationInverse Modelling for Flow and Transport in Porous Media
Inverse Modelling for Flow and Transport in Porous Media Mauro Giudici 1 Dipartimento di Scienze della Terra, Sezione di Geofisica, Università degli Studi di Milano, Milano, Italy Lecture given at the
More informationWESTCARB Phase I Results Review
WESTCARB Phase I Results Review Arizona Geologic Characterization Errol Montgomery Principal Errol L. Montgomery & Associates (520) 881-4912 emontgomery@elmontgomery.com November 9, 2005 Site Characterization
More informationTracers and Isotopes in Urban Hydrology
Tracers and Isotopes in Urban Hydrology What is a hydrologic tracer? Any substance that can be used for tracking water movement is a tracer An ideal tracer behaves exactly as the traced material behaves
More information