350 likes | 619 Views
Applications of GW model , from mathematic equations to colorful graphics. Introduction to GW model GUI ( eg.,PMWIN ) Applications on lessons Real world poblems. Proccessing MODFLOW (PMWIN). Werapol Bejranonda and Daniela Petzke
E N D
Applications of GW model,frommathematicequations to colorfulgraphics • Introduction to GW model GUI (eg.,PMWIN) • Applications on lessons • Real worldpoblems Proccessing MODFLOW (PMWIN) WerapolBejranondaand Daniela Petzke Geohydraulik und Ingenieurhydrologie, Uni Kassel, 17.06.2010
Introduction to GW model GUI Why GW model GUI ??
Introduction to GW model GUI- Whatis GW model GUI ? • GUI (graphical-user interfaces) Preprocessor and postprocessor graphical-user interfaces for preparing GW MODEL input data and viewing model output • Preprocessor and postprocessor software • Argus ONE • Groundwater Modeling System • Groundwater Vistas • ModIME • Processing Modflow • Visual Modflow • … etc.
Introduction to GW model GUI - Whyweneed GUI ? Files required in MODFLOW Sample of MODFLOW input (boundary-fileformat) 2 39 29 1 1 FREE 0 1 1 1(20I3) -1 6. DATA IBOUND Array for Layer 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 …… ......., Sample of visulization in PMWIN (contourlines) http://en.wikipedia.org/wiki/MODFLOW
Introduction to GW model GUI - Whatisthe PMWIN ? • PMWIN a pre-/post-processing of GW simulation systems supported various GW models (i.e., MODFLOW, MT3D, MT3DMS, PMPATH, PEST2000, UCODE, …) • Pre-processing Prepare/Converse/Create MODFLOW inputdatavisual. –> numerical method • Post-processing Analyze/Coverse/Presentoutput to Graphics and otherattachedmodels, like PMPATH and MT3D etc. numerical method –> Visual.
Introduction to GW model GUI - How PMWIN work on MODFLOW? conceptual model visualized results output data 30 30 9 8.993021 8.986081 8.98039 8.974634 8.96138 8.921746 8.858334 8.811367 8.792134 8.783539 8.777983 8.772 8.766218 8.761964 8.757025 8.745215 8.709628 8.652272 8.60986 8.592547 8.582397 8.575535 8.569161 8.539918 8.485279 8.391317 8.262108 8.146635 8 9 8.991026 8.984963 8.981595 8.975132 8.963697 8.928112 8.86082 8.809385 8.79072 8.782465 8.777608 8.771564 8.766438 8.761987 8.755163 8.742675 8.710963 8.64966 8.601792 8.592494 8.584826 8.576152 8.5663 8.539207 8.48043 8.391795 8.26943 8.162137 8 9 8.992068 8.984685 8.982265 8.97536 8.960163 8.926689 8.861733 8.808471 8.790564 8.783772 8.775874 8.769582 8.765231 8.760021 8.754039 8.746641 8.713333 8.644117 8.59821 8.59292 8.587028 8.578007 8.561831 8.527905 8.46121 8.34574 8.227065 8.171682 8 9 8.991956 8.986297 8.982117 8.973289 8.957443 8.914145 8.848021 8.807922 8.792113 8.786644 8.776596 8.767121 8.763245 8.759284 8.753279 8.743539 8.718549 8.652662 8.598843 8.591493 8.584188 8.575192 re-analyze MODFLOW engages real world map visualized model grid-cells design input data 2 39 29 1 1 FREE 0 1 1 1(20I3) -1 6. DATA IBOUND Array for Layer 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 …… ......., data & properties
Applications on lessons Apply the groundwater into visualized mode
Application on lessons : Standard Sample- PMWIN tutorials • Official Book 3D-groundwater modeling with PMWIN: a simulation system for modeling groundwater flow and transport processes (Wen-Hsing Chiang and Wolfgang Kinzelbach, 2005) # read free in google book and uni-ebook • Online tutorial / step-by-step http://www.cs.amherst.edu/~dac/mellon/tutorial.pdf • Sample files in folder C:\Simcore\PMWIN 5.3\Examples
Application on lessons- What is the applications? Application of MODFLOW • Schematic flow system • Effects of pumping on shallow and deep water levels (drawdown) • Natural/Artificial recharge to ground water • Surface-water and groundwater interactions • Sources of water to wells • etc. “Improving Traditional Conjunctive Use Management With Surface Waterand Groundwater Dynamic Interaction Modelling” Bejranonda W., M. Koch and S. Koontanakulvong (2009)
Application on lessons- PMWIN tutorials : Study model In folder # C:\Simcore\PMWIN 5.3\Examples\misc\misc2
Applications on lessons Real world application : A study case of Pump station design
Application on lessons : Real world- Real world : Study case Effect of Pump-station(s) at East bank of flood plain in North of Thailand Impervious rock Lake Impervious rock High Permeability soil River West bank East bank Impervious rock
Application on lessons : Real world- Map & Grid (1) Definegridsize (meshsize) 2 layars 20 Grid x 30 Grid (width x height) 200 x 300 m. (width x height)
Application on lessons : Real world- Map & Grid (2) Insertthemap as background (Grid->meshsize -> options -> Maps -> Raster graphic) Environmentaloptions -> Coordinatesystem
Application on lessons : Real world- Map & Grid (3) Finite difference grid system Real world map with hydrologic system = +
Application on lessons : Real world- Grid refinement Prepare for the wells Grid divided into 10 small columns 10 small rows
Application on lessons : Real world- Grid properties Based on finite difference method • Size of X,Y,Z dimension Cell dimension (Finite difference grid) • Aquifer type Set equations for each layer • Elevations Define grid/layer size (z dimension)
Application on lessons : Real world- Layer type Aquifer type 1.Confined 2.Unconfined 3.combination of both (conversable) pressured aquifer Layer 1 : Unconfined Layer 2 : Confined
Application on lessons : Real world- Boundary condition no flow constant head 0 for inactive cell 1 for computed head -1 for constant head
Application on lessons : Real world- Aquifer properties Aquifer’selevations • Top L1 = 15 • Top L2 = 2 • Bottom L1 = 3 • Bottom L2 = 0 Thickness of aquifer L1 = 12 m. L2 = 2 m. +15.00 Layer 1 Thickness of aquifer +3.00 +2.00 Layer 2 +0.00
Application on lessons : Real world- Aquifer properties • Initial head • - initial value for simulation • Hydraulic conductivity • properties of porous media • Vertical leakance • -bottom material property • Effective porosity
Application on lessons : Real world- Head & Elevation Hydraulic head • Initial Head L1 = 5 • Initial Head L2 = 3 • Head boundary = 10 average river stage = +10.00 m. +15.00 Layer 1 +10.00 +5.00 +3.00 +3.00 +2.00 Layer 2 +0.00
Application on lessons : Real world- Aquifer properties Layer 1 Layer 2
Application on lessons : Real world- Sources & Sinks • Such as : • Rivers • Streams • Horizontal Flow Barriers • Drains • Seepage Face • Wells
Application on lessons : Real world- Sources & Sinks : Recharge Lake +15.00 Layer 1 +10.00 +5.00 +3.00 +3.00 +2.00 Layer 2 +0.00
Application on lessons : Real world- Sources & Sinks : Wells Lake Pump +15.00 Layer 1 +10.00 +5.00 +3.00 +3.00 +2.00 Layer 2 +0.00
Application on lessons : Real world- Engage MODFLOW simulation
Application on lessons : Real world- Graphical contour lines • Results extractor -> read -> save (ASCII Matrix : head.dat) • Back to main program • Presentation -> Value -> Matrix -> Load (head.dat) • Options -> Environment -> Contours • Trick +visible / Click +Restore default • Click Level HEADER (to adjust the interval max/min)
Application on lessons : Real world- Compare 2 scenarios Case A: 1 pump x 0.001 m3/s Pump Lake +15.00 Layer 1 +10.00 +5.00 +3.00 +3.00 +2.00 Layer 2 +0.00 Case B: 2 pumps x 0.0005 m3/s Lake Pumps +15.00 Layer 1 +10.00 +5.00 +3.00 +3.00 +2.00 Layer 2 +0.00
Application on lessons : Real world- Model results : Water level Case A: 1 pump x 0.001 m3/s Case B: 2 pumps x 0.0005 m3/s
Application on lessons : Real world- Model results : Water Budget
Application on lessons : Real world- Model results : Water Table
Application on lessons : Real world- What is exactly in the real world? • Unknown parameter / Lack of field value • Heterogeneous / Uncertainty • Calibration & Verification • Degree of freedom? • etc., …. Field data Observed data Data interpretation • Adjustment ? • Which parameter • How to • Uncertainty Model calibration
Fragen & Antworten Darcy's law : “Les Fontaines Publiques de la Ville de Dijon”, 1856 ….BIOMOC (DOS/DG/SGI/Sun) Version 1.0, 1999/03/10CFP (Win) Version 1.0, 2008/04/14GSFLOW (Linux/Win) Version 1.1.1, 2010/02/12GWM-2000 and GWM-2005 (Win) Version 1.1.3 and 1.2.1, 2010/01/26HST3D (Unix/Linux/Win) Version 2.2.13, 2006/05/26HYDROTHERM (Linux/Win) Version 3.1.0, 2008/03/01INFIL3.0 (Win) Version 1.0, 2008/06/19JUPITER API (DOS/Win) Version 1.4.0, 2010/04/06MF2K-FMP and MF2005-FMP2 (Win) Version 1.0.00 and 1.0.00, 2009/10/28 and 2006/05/19MF2K-GWT (DOS/Unix) Version 1.9.8, 2007/10/22MMA (Sun/DOS) Version 1.100, 2009/11/25MODFE (DG/Sun) Version 1.2, 1998/07/08MODFLOW, MODFLOW-2005, and related programsMODOPTIM (Win) Version 1.0, 2005/05/05OPR-PPR (Win/DOS) Version 1.00, 2007/08/13PHAST (Linux/Win) Version 1.5.1, 2008/03/24PHREEQC (Mac/Linux/Win) Version 2.15.0, 2008/02/05SEAWAT Version 4.00.03, 2009/09/22SHARP (DOS/DG/SGI/Sun) Version 1.1, 09/03/2004STLK1 and STWT1 (Win) Version 1.0, 2009/12/04STRMDEPL08 Version 1.1, 2010/01/26SUTRA and related programsSUTRA-MS (DOS) Version 2004.1, 2005/05/09TopoDrive and ParticleFlow (Sun/Win) Version 1.0, 2002/02/13UCODE_2005 (DOS/UNIX) Version 1.019, 2010/04/06 What’s next ? , 2010…