If
a solution is expressed in terms of head rather than drawdown, then the initial
heads must be defined. Numerical modeling is conducted in terms of head and allows
us to define complex initial conditions. There are some important concepts to
keep in mind: Now lets
explore the transient character of the MODFLOW problem that we have been using
as a class example. By undertaking a transient modeling exercise involving variable
recharge and pumping.
1. Initial material properties and boundary conditions
must be consistent with the initial heads. If you start with initial heads that
are contoured from field measurements and you do not apply a stress to the system,
the heads will adjust to the properties and boundaries, so you are inadvertently
introducing a stress by defining inconsistent values for starting conditions.
The most common way to deal with this problem is to calculate a pre-stress steady
state solution for use as initial heads.
2. If the field system being simulated is not in equilibrium, an earlier equilibrium
condition can be identified and defined as a starting point. All subsequent stresses
must be simulated from the time when equilibrium prevailed until time when the
initial conditions are needed is reached, then the early stresses must continue
along with the new stress of interest if the early stresses continue in the field.
A classic example is this hydrograph of a well on Long Island displaying the impact
of sewering which started in the mid 50's and was completed in mid-60's.
The upper dashed line represents a long term steady state condition. There are
daily and seasonal variations but the average values are relatively constant in
time. The lower dashed line represents the new long term steady state condition
that prevails after the full impact of sewering has been realized.
The water levels decline as a result of sewering because recharge that used to
enter the ground water via septic tanks is now routed away by the sewers.
Suppose you wanted to simulate the impact of a proposed municipal water supply
well that would go on line in 1962. In order to obtain a consistent set of starting
values for initial conditions in 1962, it would be necessary to simulate the early
equilibrium that prevailed until 1955, add the stress of the gradually reducing
recharge in space through time until reaching 1957, then continue simulating the
reductions in recharge along with the addition of the pumping well.
It is vital that you develop hydrographs of an area and come
to understand these transitions before undertaking an elaborate modeling project.
Another item to note while viewing this diagram is this.
If you had water level data from a number of wells prior to 1955, the average
of those values would reflect the average water levels during the early period
for one set of boundary conditions and stresses. However if data from 1950 through
1975 were mixed together and averaged, the result will NOT reflect a set of boundary
conditions and stresses that occurred together in the field.
3. What can we do if we cannot use a steady state initial condition because our
problem is dependent on a short term response during a particular time of year?
In this case there may be a quasi-cyclic equilibrium in the area, that is a general
repetition of conditions that occur year after year and we may attempt to represent
that situation, then start our stress at the appropriate point in the cycle.
Draw diagrams of the variations and the timing of your stress and coordinate this
with a diagram of model stress periods and time steps to help yourself get everything
timed correctly in your model.
4. What if we do not have sufficient data to establish an acceptable steady initial
condition to commence our cyclic equilibrium? In this case, we may be able to
start with a rather arbitrary initial condition but simulate the cycle long enough
such that we simulate the same values at the same times in subsequent cycles.
5. It may be that there is enough information in your transient data to estimate
initial conditions. This can be judged by evaluating the sensitivity of initial
conditions to the available data and the correlation of parameters defining those
conditions to other parameters that are estimated.
6. Finally, it is useful to note that the further, in time, the simulation is
from the initial conditions the less influence those initial conditions have on
the simulated values.
7. Recall the procedure for discretizing
time as discussed in Unit 7.
Click here to undertake the transient exercise
GO "BACK" to the transient main page