Comments
Description
Transcript
Document
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