Introduction To Groundwater Modelling: C. P. Kumar
Introduction To Groundwater Modelling: C. P. Kumar
Introduction To Groundwater Modelling: C. P. Kumar
C. P. Kumar
Scientist F
National Institute of Hydrology
Roorkee 247667 (Uttaranchal)
India
Email: cpkumar@yahoo.com
Webpage: http://www.angelfire.com/nh/cpkumar/
Presentation Outline
Groundwater in Hydrologic Cycle
Why Groundwater Modelling is needed?
Mathematical Models
Modelling Protocol
Model Design
Calibration and Validation
Groundwater Flow Models
Groundwater Modelling Resources
Soil
Moisture
Ground water
Groundwater
Baseflow
Linked to SW systems
Sustains flows
in streams
Groundwater Concerns?
pollution
groundwater mining
subsidence
Groundwater
The total volume that may be withdrawn annually from the aquifer.
dh
q
=
dx
K
Integrated
q
h0 dh = K
h
Poten
tiome
tri
Surfa c
ce
q
K
ho x
0
h(x,y,z,t)?
qx
dx h h0 =
K
x h(x)
x
h( x ) = h0
qx
K
Groundwater Modeling
The only effective way to test effects of
Groundwater flow
calculate both heads and flow
MODELING PROCESS
TYPES OF MODELS
CONCEPTUAL MODEL QUALITATIVE DESCRIPTION OF SYSTEM
"a cartoon of the system in your mind"
MATHEMATICAL MODEL MATHEMATICAL DESCRIPTION OF
SYSTEM
SIMPLE - ANALYTICAL (provides a continuous solution over the
model domain)
COMPLEX - NUMERICAL (provides a discrete solution - i.e. values are
calculated at only a few points)
ANALOG MODEL e.g. ELECTRICAL CURRENT FLOW through a
circuit board with resistors to represent hydraulic conductivity and
capacitors to represent storage coefficient
PHYSICAL MODEL e.g. SAND TANK which poses scaling problems
Mathematical Models
Mathematical model:
simulates ground-water flow and/or
solute fate and transport indirectly by
means of a set of governing equations
thought to represent the physical
processes that occur in the system.
(Anderson and Woessner, 1992)
z
x
y
1. Consider flux (q) through REV
2. OUT IN = - Storage
3. Combine with: q = -K grad h
q = - Ss (h t)
q = - K grad h
(Darcys Law)
div (K grad h) = Ss (h t)
(Ss = S / z)
h
( Kx ) +
( Ky ) +
( Kz ) = 0
x
x
y
z
y
z
with a source/sink term
h
( Kx ) +
( Ky ) +
( Kz ) = R *
x
x
y
y
z
z
h
h
h
( Kx ) +
( Ky ) +
( K z ) = Ss
R*
x
x
y
y
z
z
t
Specific Storage
Ss = V / (x y z h)
h
b
S=V/Ah
S = Ss b
Confined aquifer
Unconfined aquifer
Storativity
Specific yield
Figures taken from Hornberger et al. (1998)
General 3D equation
h
h
( Kx ) +
( Ky ) +
( K z ) = Ss
R*
x
x
y
y
z
z
t
2D confined:
2D unconfined:
h
h
(T x ) +
(T y ) = S
R
x
x
y
y
t
h
h
( hK x ) +
( hK y ) = S
R
x
x
y
y
t
Numerical Methods
hAll numerical methods involve
representing the flow domain by a
limited number of discrete points called
nodes.
hA set of equations are then derived to
relate the nodal values of the
dependent variable such that they
satisfy the governing PDE, either
approximately or exactly.
Numerical Solutions
Discrete solution of head at selected nodal points.
Involves numerical solution of a set of algebraic
equations.
Finite difference models (e.g., MODFLOW)
Finite element models (e.g., SUTRA)
K1
K2
K3
K4
Hybrid
Analytic Element Method (AEM)
Involves superposition of analytic solutions. Heads are
calculated in continuous space using a computer to do
the mathematics involved in superposition.
The AE Method was introduced by Otto Strack.
A general purpose code, GFLOW, was developed by
Stracks student Henk Haitjema, who also wrote a
textbook on the AE Method: Analytic Element Modeling
of Groundwater Flow, Academic Press, 1995.
Currently the method is limited to steady-state,
two-dimensional, horizontal flow.
Modelling Protocol
What is a model?
Analytical
Numerical
Modelling Protocol
Model Overkill?
ETHICS
Conceptual
Model
Everything should be made as simple as possible, but not simpler. Albert
Einstein
Model Design
Steady/Unsteady..1,
Steady/Unsteady..1,2,2,or
or3-D;
3-D;
Heterogeneous/Isotropic..Instantaneous/Continuous
Heterogeneous/Isotropic..Instantaneous/Continuous
Calibration
Show that Model can reproduce fieldmeasured heads and flow (concentrations if
contaminant transport)
Results in parameter data set that best
represents field-measured conditions.
Model Verification
Prediction
Presentation of Modelling
Design and Results
Effective Communication of
Modeling Effort
Postaudit
Model Redesign
NUMERICAL MODELING
DISCRETIZE
Write equations of GW Flow between each node
Darcy's Law
Conservation of Mass
Define
Material Properties
Boundary Conditions
Initial Conditions
Stresses
Calibrate
Steady State
Transient
Validate
Sensitivity
Predictions
Similar Process for Transport Modeling only Concentration and Flux is unknown
NUMERICAL MODELING
Model Design
MODELs NEED
Geometry
Material Properties (K, S, T, e, R, etc.)
Boundary Conditions (Head, Flux, Concentration etc.)
Stress - changing boundary condition
Model Design
Conceptual Model
Selection of Computer Code
Model Geometry
Grid
Boundary array
Model Parameters
Boundary Conditions
Initial Conditions
Stresses
Concept Development
Developing a conceptual model is the initial
and most important part of every modelling
effort. It requires thorough understanding
of hydrogeology, hydrology and dynamics
of groundwater flow.
Conceptual Model
A descriptive representation
of a groundwater system that
incorporates an interpretation of the
geological & hydrological conditions.
Generally includes information about
the water budget. May include
information on water chemistry.
Model Geometry
Model geometry defines the size and the
shape of the model. It consists of model
boundaries, both external and internal, and
model grid.
Boundaries
Physical boundaries are well defined
geologic and hydrologic features that
permanently influence the pattern of
groundwater flow (faults, geologic units,
contact with surface water etc.)
Boundaries
Hydraulic boundaries are derived from the
groundwater flow net and therefore
artificial boundaries set by the model
designer. They can be no flow boundaries
represented by chosen stream lines, or
boundaries with known hydraulic head
represented by equipotential lines.
HYDRAULIC BOUNDARIES
A streamline (flowline) is also a
hydraulic boundary because by
definition, flow is ALWAYS
parallel to a streamflow. It can
also be said that flow NEVER
crosses a streamline; therefore it
is similar to an IMPERMEABLE
(no flow) boundary
BUT
Stress can change the flow
pattern and shift the position of
streamlines; therefore care must
be taken when using a
streamline as the outer boundary
of a model.
SPECIFIED HEAD OR
CONSTANT HEAD BOUNDARY
h = constant
q is determined by the model.
And may be +ve or ve according
to the hydraulic gradient developed
Boundary Types
Specified Head/Concentration: a special case of constant head (ABC, EFG)
Constant Head /Concentration: could replace (ABC, EFG)
Specified Flux: could be recharge across (CD)
No Flow (Streamline): a special case of specified flux (HI)
Head Dependent Flux: could replace (ABC, EFG)
Free Surface: water-table, phreatic surface (CD)
Seepage Face: pressure = atmospheric at ground surface (DE)
Known Flux
Recharge
Evapotranspiration
Wells
Stream
No Flow boundaries
Initial Conditions
Values of the hydraulic head for each active
and constant-head cell in the model. They
must be higher than the elevation of the cell
bottom.
For transient simulation, heads to resemble
closely actual heads (realistic).
For steady state, only hydraulic heads in
constant head-cell must be realistic.
Model Parameters
Time
Space (layer top and bottom)
Hydrogeologic characteristics
(hydraulic conductivity, transmissivity,
storage parameters and effective porosity)
Time
Time parameters are specified when
modelling transient (time dependent)
conditions. They include time unit, length
and number of time steps.
Length of stress periods is not relevant for
steady state simulations
Grid
In Finite Difference model, the grid is
formed by two sets of parallel lines that are
orthogonal. The blocks formed by these
lines are called cells. In the centre of each
cell is the node the point at which the
model calculates hydraulic head. This type
of grid is called block-centered grid.
Grid
Grid mesh can be uniform or custom, a
uniform grid is better choice when
Evenly distributed aquifer characteristics data
The entire flow field is equally important
Number of cells and size is not an issue
Grid
Grid mesh can be custom when
There is less or no data for certain areas
There is specific interest in one or more smaller
areas
MODEL GRIDS
Grids
h It is generally agreed that from a practical
point-of-view the differences between grid
types are minor and unimportant.
h USGS MODFLOW employs a body-centred grid.
Effective porosity
Required to calculate velocity, used mainly
in solute transport models
Calibration Targets
associated error
calibration
value
+/0.80 m
20.24 m
Calibration to Fluxes
When recharge rate (R) is a calibration
parameter, calibrating to fluxes can help in
estimating K and/or R.
q = KI
H1
H2
Calibration - Remarks
Calibrations are non-unique.
A good calibration does not ensure that
the model will make good predictions.
You can never have enough field data.
Modelers need to maintain a healthy skepticism
about their results.
Need for an uncertainty analysis to accompany
calibration results and predictions.
Sensitivity analysis
Stochastic simulation
Modeling Chronology
1960s Flow models are great!
1970s Contaminant transport models are great!
1975
MODFLOW
USGS code
Finite Difference Model
MODFLOW 88
MODFLOW 96
MODFLOW 2000
MODFLOW
(Three-Dimensional Finite-Difference Ground-Water Flow
Model)
MT3D
(A Modular 3D Solute Transport Model)
FEFLOW
(Finite Element Subsurface Flow System)
FEFLOW is a finite-element package for simulating 3D and 2D
fluid density-coupled flow, contaminant mass (salinity) and
heat transport in the subsurface.
HST3D
(3-D Heat and Solute Transport Model)
The Heat and Solute Transport Model HST3D simulates
ground-water flow and associated heat and solute transport in
three dimensions.
SEAWAT
(Three-Dimensional Variable-Density Ground-Water Flow)
The SEAWAT program was developed to simulate threedimensional, variable- density, transient ground-water flow
in porous media.
SUTRA
(2-D Saturated/Unsaturated Transport Model)
SWIM
(Soil water infiltration and movement model)
VISUAL HELP
(Modeling Environment for Evaluating and Optimizing
Landfill Designs)
Visual MODFLOW
(Integrated Modeling Environment for MODFLOW and
MT3D)
: http://groups.yahoo.com/group/gwmodel/
: gwmodel@yahoogroups.com
: gwmodel-subscribe@yahoogroups.com
: gwmodel-unsubscribe@yahoogroups.com
: gwmodel-owner@yahoogroups.com
: http://in.groups.yahoo.com/group/visual-modflow/
: visual-modflow@yahoogroups.co.in
: visual-modflow-subscribe@yahoogroups.co.in
: visual-modflow-unsubscribe@yahoogroups.co.in
: visual-modflow-owner@yahoogroups.co.in
THANKS
HAPPY MODELLING