1 / 27

SHUVALOVO area First step GWFM

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

lilac
Télécharger la présentation

SHUVALOVO area First step GWFM

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. 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

  2. 1. Problem Definition 2. Conceptual Model 3. MathematicalModel 4. Numerical Formulation 5. Sensitivity Analysis 6. Calibration 7. Model Application Modelling phases procedure

  3. 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 • SOIL-GIS 2008 • Definition of Modflow and MT3D input: grid cell dimensions, boundary conditions, layers geometry, permeability, recharge, withdrawals, porosity, pollutant concentration, ecc … • 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 PHASE 1 PHASE 2 PHASE 3 PHASE 4 PHASE 5-6 PHASE 7 Modelling phases procedure

  4. GG1 cross-section Conceptual model All layers are continuous in the model !!!

  5. 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 Numerical formulation Top of prequaternary is the bottom of 8 layer

  6. Numerical formulation Tops Elevation Table

  7. Numerical formulation

  8. 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 Tops correction

  9. In this example layer 8 bottom (below/green) intersects layer 8 top (above/grey) in three areas Tops correction

  10. 3D Surfer representation 9 surfaces elevation

  11. 3D Surfer representation 9 surfaces elevation

  12. 3D Surfer representation 9 surfaces elevation

  13. Tables GEODB Database Shapefiles Raster Geological and Hydrogeological interpretation Interpolation GridConvert ESRI grid for SOIL-GIS 3D Surfer grid Modelling methodology

  14. Constant head Lakes as constant head Rivers Lakes as constant head Rivers GW aquifer Boundary conditions (layers 1-2-3)

  15. ? Constant head Constant head IM aquifer Boundary conditions (layers 5-6-7)

  16. Grid cells of 10-50 meters are the ideal to study Shuvalovo area and in particular the lakes, the most important elements to modelling. Discretization domain

  17. 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. In the next future …

  18. K range ? ? Is correct? Permeability

  19. ? Dusty sand in Lithology.shp Fine sand in Boreholes.mdb Permeability

  20. Various grain size sand in Lithology.shp Medium grain size sand in Boreholes.mdb ? Permeability

  21. Piezometry Differences between GW_Water_level_max.shp in blue and Daylight_Relief.shp in black

  22. Piezometry Differences between GW_Water_level_max.shp in blue and interpolation of GW boreholes maximum water level in red

  23. Piezometry Differences between IM_Water_level_max.shp in blue and interpolation of IM boreholes minimum water level in red

  24. DEM Differences between Daylight_Relief.shp in black and interpolation of altitude points in red

  25. Static and dynamic piezometric surface (checked with target control) Water exchange estimate between GW aquifer and IM aquifer according to reciprocal K values Flow expected results

  26. 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

  27. THANKS for your precious attention !!!

More Related