GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
Integrating Geological Information in City
Management to Prevent Environmental Risks
GeoInforM
LIFE06 TCY/ROS/000267
SHUVALOVO area
First step GWFM
Michel Fasoli
Piergiorgio Valentini, Roberta Donati
Province of Milan
Environmental Information System Department
Central Directorate Environmental Resources
direzione
centrale
risorse
ambientali
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
Modelling phases procedure
1. Problem Definition
2. Conceptual Model
3. Mathematical Model
4. Numerical Formulation
5. Sensitivity Analysis
6. Calibration
7. Model Application
direzione
centrale
risorse
ambientali
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
direzione
centrale
risorse
ambientali
Modelling phases procedure
PHASE 1
PHASE 2
•Aims identification for modelling
•Pilot area identification
•Collecting data about pilot area and regional area
•Organizing data with GIS software first of all in an ArcMap project and
then in a geodatabase with all original data e new data elaboration
•Geological and hydrogeological data interpretation
•Conceptual model definition
PHASE 3
SOIL-GIS 2008
PHASE 4
Definition of Modflow and MT3D input: grid cell dimensions, boundary
conditions, layers geometry, permeability, recharge, withdrawals,
porosity, pollutant concentration, ecc …
PHASE 5-6
PHASE 7
•Flow simulation with SOIL-GIS 2008
•Transport simulation with SOIL-GIS 2008
•Sensitivity analysis compared with K
•Calibration limits compared with target wells
Application of the model in the case study after the calibration
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
Conceptual model
GG1 cross-section
All layers are continuous in the model !!!
direzione
centrale
risorse
ambientali
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
Numerical formulation
Elevation data from:
• Altitude points
• Boreholes (from
database and
from sections)
• Water levels
• Daylight relief
• isolines
• Lakes fictitious
• points
• Top inter moraine isolines
• Prequaternary relief
isolines
• Tops correction
fictitious points
Top of prequaternary
is the bottom of 8 layer
direzione
centrale
risorse
ambientali
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
Numerical formulation
Tops Elevation Table
direzione
centrale
risorse
ambientali
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
Numerical
formulation
direzione
centrale
risorse
ambientali
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
direzione
centrale
risorse
ambientali
Tops correction
• The interpolation for some layers (with few
elevation points) gives problems of intersections
between contiguous layers (see example in the
next slide)
• In this case is necessary to put in some fictitious
points to avoid the problem
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
Tops correction
In this example layer 8 bottom (below/green) intersects layer 8 top
(above/grey) in three areas
direzione
centrale
risorse
ambientali
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
9 surfaces elevation
3D Surfer representation
direzione
centrale
risorse
ambientali
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
9 surfaces elevation
3D Surfer representation
direzione
centrale
risorse
ambientali
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
9 surfaces elevation
3D Surfer representation
direzione
centrale
risorse
ambientali
direzione
centrale
risorse
ambientali
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
Modelling methodology
GEODB
Tables
Database
Shapefiles
Raster
Geological and
Hydrogeological
interpretation
Interpolation
GridConvert
ESRI grid for
SOIL-GIS
3D Surfer grid
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
direzione
centrale
risorse
ambientali
GW aquifer Boundary conditions
(layers 1-2-3)
Constant
head
Lakes as
constant
head
Rivers
Rivers
Lakes as
constant
head
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
direzione
centrale
risorse
ambientali
IM aquifer Boundary conditions
(layers 5-6-7)
?
Constant
head
Constant
head
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
direzione
centrale
risorse
ambientali
Discretization domain
Grid cells of 10-50 meters are the ideal to
study Shuvalovo area and in particular the
lakes, the most important elements to
modelling.
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
direzione
centrale
risorse
ambientali
In the next future …
The next steps will be the definition of:
• permeability grids for each layers (we will
simulate a situation with minimum K values
and another with maximum K values);
• porosity grids estimate from K values;
• withdrawals grids (known or estimated);
• recharge grids (% infiltration);
• wells as control target.
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
direzione
centrale
risorse
ambientali
K range
?
Permeability
?
Is correct?
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
Permeability
?
Dusty sand in Lithology.shp
Fine sand in Boreholes.mdb
direzione
centrale
risorse
ambientali
direzione
centrale
risorse
ambientali
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
Permeability
Various grain size sand
in Lithology.shp
Medium grain size sand
in Boreholes.mdb
?
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
direzione
centrale
risorse
ambientali
Piezometry
Differences between GW_Water_level_max.shp in blue and
Daylight_Relief.shp in black
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
direzione
centrale
risorse
ambientali
Piezometry
Differences between GW_Water_level_max.shp in blue and interpolation of GW
boreholes maximum water level in red
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
direzione
centrale
risorse
ambientali
Piezometry
Differences between IM_Water_level_max.shp in blue and interpolation of IM
boreholes minimum water level in red
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
direzione
centrale
risorse
ambientali
DEM
Differences between Daylight_Relief.shp in black and interpolation of altitude
points in red
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
direzione
centrale
risorse
ambientali
Flow expected results
• Static and dynamic piezometric surface
(checked with target control)
• Water exchange estimate between GW
aquifer and IM aquifer according to
reciprocal K values
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
direzione
centrale
risorse
ambientali
Future actions
?
To respect the future actions of the Province of Milan
we need to know the data necessary for a complete flow and transport simulation
GeoInforM Project - LIFE06 TCY/ROS/000267
3° Steering Committee – ST. Petersburg, 17th-18th April 2008
THANKS
for your precious attention !!!
direzione
centrale
risorse
ambientali