1 / 118

Modelling of groundwater flow and solute transport

Mathematical models in hydrogeology. Modelling of groundwater flow and solute transport. Summer School Polish Geological Institute. Dr.ir. Gualbert Oude Essink 1. Netherlands Institute of Applied Geosciences 2. Faculty of Earth Sciences, Free University of Amsterdam. PGI: Aug. 22-24, 2002.

josie
Download Presentation

Modelling of groundwater flow and solute transport

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Mathematical models in hydrogeology Modelling of groundwater flow and solute transport Summer School Polish Geological Institute Dr.ir. Gualbert Oude Essink 1. Netherlands Institute of Applied Geosciences 2. Faculty of Earth Sciences, Free University of Amsterdam PGI: Aug. 22-24, 2002

  2. Curriculum Vitae • Delft University of Technology, Civil Engineering: till 1997 • Ph.D.-thesis: Impact of sea level rise on groundwater flow regimes • Utrecht University, Earth Sciences: till 2002 • Lectures in groundwater modelling and transport processes • Variable-density groundwater flow • Salt water intrusion and heat transport • Netherlands Institute of Applied Geosciences (NITG): till 2004 • Salinisation processes in Dutch coastal aquifers • Free University of Amsterdam, Earth Sciences: till 2004 • EC-project CRYSTECHSALIN (o.a. with J. Bear) • Crystallisation processes in porous media NITG: g.oudeessink@nitg.tno.nl VU: oudg@geo.vu.nl

  3. Programme (I) Modelling of groundwater flow and solute transport • Modelling protocol • Discretisation Partial Differential Equation (PDE) • Groundwater flow: MODFLOW • Solute transport: MOC3D PGI: Aug. 22-24, 2002

  4. Programme (II) Thursday Augustus 22, 2002: modelling protocol Darcy’s Law, steady state continuity equation Laplace-equation Taylor series development boundary conditions: Dirichlet, Neumann, Cauchy non-steady state: explicit, implicit, Crank-Nicolson • Friday Augustus 23, 2002: • groundwater flow with MODFLOW: mathematical description • short introduction Graphical User Interface PMWIN Saturday Augustus 24, 2002: solute transport with MOC3D: mathematical description numerical dispersion PGI: Aug. 22-24, 2002

  5. The Hydrological Circle

  6. Ten steps of the Modelling Protocol • 1. Problem definition • 2. Purpose definition • 3. Conceptualisation • 4. Selection computer code • 5. Model design • 6. Calibration • 7. Verification • 8. Simulation • 9. Presentation • 10. Postaudit

  7. Modelling protocol Why numerical modelling? • +: • cheaper than scale models • analysis of very complex systems is possible • a model can be used as a database • simulation of future scenarios • -: • simplification of the reality • only a tool, no purpose on itself • garbage in=garbage out: (field)data important • perfect fit measurement and simulation is suspicious

  8. Modelling protocol 3. Conceptualisation (I) • Model is only a schematisation of the reality • Which hydro(geo)logical processes are relevant? • Which processes can be neglected? • Boundary conditions • Variables and parameters: • subsoil parameters • fluxes in and out • initial conditions • geochemical data • Mathematical equations

  9. Modelling protocol 3. Concept (II): example of salt water intrusion • Relevant processes: • groundwater flow in a • heterogeneous porous • medium • solute transport • variable density flow • natural recharge • extraction of groundwater The Netherlands • Negligible processes: • heat flow • swelling of clayey aquitard • non-steady groundwater flow • Boundary conditions • no flow at bottom • flux in dune area • constant head in polder area

  10. 4. Selection computer code (II) • Groundwater • computer • codes Hemker (1994)

  11. 4. Selection computer code There are numerous good groundwater computer codes available, so don’t make your own code! • See the internet, e.g.: • U.S. Geological Survey: water.usgs.gov/nrp/gwsoftware/ • Scientific Software Group: www.scisoftware.com/ • Waterloo Hydrogeologic: www.flowpath.com/ • www.groundwatersoftware.com/newsletter.htm • Groundwater Digest: groundwater@groundwater.com

  12. Modelling protocol 5. Model design (I) • Choice grid Dx: • depends on natural variation • in the groundwater system • concept model • data collection Choice time step Dt • Conditions: • initial conditions • boundary conditions: • Dirichlet: head • Neumann: flux, e.g., no flow • Cauchy/Robin: mixed boundary condition

  13. Modelling protocol 5. Model design (II): example Geometry, subsoil parameters, boundary conditions

  14. Modelling protocol 5. Model design (III): choice time step Dt

  15. Modelling protocol 6. Calibration (I) • Fitting the groundwater model: is your model okay? • trial and error • automatic parameter estimation/inverse modelling (PEST, UCODE)

  16. Modelling protocol 6. Calibration (II): example Measured and computed freshwater heads

  17. Modelling protocol 6. Calibration (III): example • Measured and computed freshwater heads • 95 useful observation wells • maximum measured head in area: +1.83 m • minimum measured head in area: -4.98 m • average (measured -computed) = -0.04 m • average absolute|measured -computed| = 0.34 m • standard deviation = 0.46 m

  18. Modelling protocol 6. Calibration: errors during modelling protocol • Wrong model concept • resistance layer not considered • heat transport on groundwater flow not considered (r not constant) • Incomplete equations • decay term of solute transport not considered • Inaccurate parameters and variables • Errors in computer code • Numerical inaccuracies • Dt, Dt, numerical dispersion

  19. Modelling protocol 7. Verification: testing the calibrated model • 1D, 2D -->3D • stationary dataset --> non-steady state simulation ‘verification problem’: there is always a lack of data

  20. Modelling protocol 8. Simulation 9. Presentation • Simulation of scenarios • Computation time depends on: • computer speed • size model • efficiency compiler • output format

  21. Modelling protocol 10. Postaudit • Postaudit: • analysing model results after a long time • Anderson & Woessner(’92): • four postaudits from the 1960’s • Errors in model results are mainly caused by: • wrong concept • wrong scenarios

  22. Modelling scheme Anderson & Woessner (1992) Modelling protocol

  23. gives where K= hydraulic conductivity [L/T ] Discretisation PDE I. Darcy’s law (1856)

  24. In pressures: Definition piezometric head: q=Darcian specific discharge [L/T] f=piezometric head [L] p=pressure [M/LT2] k=intrinsic permeability [L2] m=dynamic viscosity [M/LT] r=density [M/L3] g=gravitaty [L/T2] ki=hydraulic conductivity [L/T] Discretisation PDE II. Darcy’s Law

  25. Gives: if then: gives Discretisation PDE III. Darcy’s Law

  26. Discretisation PDE Analogy physical processes

  27. Discretisation PDE Continuity equation: steady state Mass balance [M T-1] r=density [M L-3]

  28. + Continuity equation Flow equation (Darcy’s Law) = Groundwater flow equation If ki=constant and r=constant then: Laplace equation Discretisation PDE Steady state groundwater flow equation (PDE=Partial Differential Equation)

  29. - Discretisation PDE Taylor series development (I) Best estimate of f i+1 is based on f i

  30. + Discretisation PDE Taylor series development (II)

  31. Discretisation in x-direction (i ): Discretisation in y-direction (j ): If Dx = Dy then: ‘Fivepoint operator’ Discretisation PDE Laplace equation in 2D

  32. Discretisation PDE Fivepoint operator (I): constant head example

  33. Discretisation PDE Fivepoint operator (II): no-flow example Nodes on the edges of an element

  34. Characteristics: Hydraulic conductivity k=10 m/dag Groundwater recharge N=0 mm/dag Thickness aquifer D=50m Boundary conditions BOUNDARY 1: no-flow near the mountains BOUNDARY 2: linear from 100 m (near mountains) to 0 m (near sea) BOUNDARY 3: constant seawater level of 0 m BOUNDARY 4: no-flow Discretisation PDE Fivepoint operator (III): example

  35. Discretisation PDE Fivepoint operator (IV): example On BOUNDARY 1: On BOUNDARY 4: EXCEL example!

  36. Discretisation PDE Fivepoint operator (V): effect convergence criterion • Solution I: • Initial estimate head=50m • Convergence-criterion=1m • Number of iterations=13 • Head 1=29,17m • Head 2=25,11m • Solution II: • Initial estimate head=50m • Convergence-criterion=0.01m • Number of iterations=74 • Head 1=20,67m • Head 2=23.07m

  37. Jacobi iteration: only old values • Gauss-Siedel iteration: also two new values • Overrelaxation: Discretisation PDE Numerical schemes: iterative methods

  38. + Non-steady state continuity equation = Flow equation (Darcy’s Law) Groundwater flow equation Ss=specific storativity [1/L] W’=source-term Multiply with constant thickness D of the aquifer gives: S=elastic storage coefficient [-] T=kD= transmissivity [L2/T] non-steady state Non-steady state groundwater flow equation

  39. Explicit (‘forwards in space’): non-steady state Explicit numerical 1D solution (I) One-dimensional non-steady state groundwater flow equation: • Properties: • direct solution • can be numerical instable

  40. Groundwater flow in aquifer between two ditches: recharge=1mm/day ditch ditch Dx=10 phreatic storage=1mm/day I. II. non-steady state Explicit numerical 1D solution (II): example • Subsoil parameters: • storage=1/3 • T=kD=10 m2/day • N=0.001 m/day • Numerical parameters: two sets • timestep Dt=1 day & Dt=10 days • length step Dx=10m • NDt/m=0.003 m & 0.03 m • TDt/mDx2=0.3 & 3 EXCEL example!

  41. non-steady state Explicit numerical 1D solution (III): stability analysis Suppose:

  42. non-steady state Implicit numerical solution One-dimensional non-steady groundwater flow equation (now no source-term W): Implicit (‘backwards in space’): • Characteristics: • not a direct solution • solving with matrix Af=R • numerical always stable • needs more memory

  43. non-steady state Crank-Nicolson numerical solution One-dimensional non-steady state groundwater flow equation (now no source term W): Crank-Nicolson (‘central in space’): a=0.5 • Characteristics: • solving as implicit • stable with temporarily oscillations

  44. non-steady state Numerical schemes

  45. MODFLOW MODFLOW • a modular 3D finite-difference ground-water flow model • (M.G. McDonald & A.W. Harbaugh, from 1983 on) • USGS, ‘public domain’ • non steady state • heterogeneous porous medium • anisotropy • coupled to reactive solute transport • MOC3 (Konikow et al, 1996) • MT3D, MT3DMS (Zheng, 1990) • RT3D • easy to use due to numerous Graphical User Interfaces (GUI’s) • PMWIN, GMS, Visual Modflow, Argus One, Groundwater Vistas, etc.

  46. MODFLOW Nomenclature MODFLOW element

  47. MODFLOW MODFLOW: start with water balance of one element [i,j,k]

  48. MODFLOW Continuity equation (I) In - Out = Storage

  49. MODFLOW Continuity equation (II) In = positive

  50. MODFLOW Flow equation (Darcy’s Law) where is the conductance [L2/T]

More Related