Groundwater Management Optimization. Dr. Issam NOUIRI

Transcription

1 Tenth Gulf Water Conference, April 2012, Doha, Qatar Groundwater Management Optimization by ALL_WATER_gw New problem formulation & additional functionalities. Dr. Issam NOUIRI National Institute of Agronomy of Tunisia

2 Outline Problematic, General conceptual model and objective, Problem formulation and resolution methodology, ALL_WATER_gw presentation, Case study application. 28/06/2012 2

3 Problematic How much water has to be taken from each well? 28/06/2012 3

4 General Conceptual Model Demand site nodes D(d) D(1) D(d, t) & Qty(d) WEAP Transmission links F max D(se,d) & Cost(se,d) Groundwater nodes S(se) S(1) Qty (se) & Vol(se) MODFLOW Wells Max Abstraction(well) & Max Drawdown(well) 28/06/2012 4

5 Objective The main objective of this work was to develop a multi-objective optimization tool for groundwater management. 28/06/2012 5

6 Problem formulation (1/4) OBJECTIVES FUNCTIONS Objective 1 : Min f Max( Nse max FD( se, d, t) D( d, t)) / D( d, t) Satisfy DS the water demands of the sites d at every time step t. se1 d=1,, NDmax and t=δt,, Tmax Objective 2 : DWDW( w) VMaxDD ( w) Min f DD Max( ) ( Hi( w) BOTM ( w))* I( w) Minimize the maximal drawdown. 28/06/2012 6

7 Problem formulation (2/4) Objective 3 : Tmax NSEmax NDmax Minimize ( the ( unit FD ( se, dcost, t) ( VOLof ( se, water. d, t) VOS( Min f C tt se1 d 1 Tmax NSEmax NDmax ( tt se1 d 1 F D se, t)) FOL( se, d, t) CCL( se, d, t)) CCS( se, t) FOS( se, t)) ( se, d, t)) T max ( CMaxS CMaxL) 28/06/2012 7

8 Problem formulation (3/4) CONSTRAINTS F D ( se, d, t) F D( se, d) max se, d and t 1. Water supply from each groundwater node to any demand site must be lower than the maximal transmission link capacity. 2. Water Q( well abstraction, t) Q from max( each well well ) must be lower than the maximal acceptable value. DD( Cell ( well ), t) DDmax( Cell ( well )) 3. Water drawdown in each well cell must be lower than the maximal acceptable value. 4. Water quality supplied to each of the demand sites must be lower than the maximal acceptable value. 28/06/2012 8

9 Resolution methodology Fitness 1 p 1 f DS (1 MaxEQty ) p2 f DD ( 1 p 1 p 2 MaxEQag) Fitness 2 f c 1. Multi-objective Genetic Algorithm 2. PARETO optimality Concept 28/06/2012 9

10 Start year Present End year WEAP - MODFLOW DSS Historic MODFLOW model Future MODFLOW model MODFLOW simulation to compute the cell heads Abstractions optimization 28/06/

11 ALL_WATER_gw in the Framework WEAP MODFLOW DSS WEAP Model MODFLOW Historic and Future Models Optimal Abstractions WEAP-MODFLOW DSS ALL_WATER_gw MODFLOW Model User: Priorities - Constraints Options - Parameters Outputs for analyses 28/06/

12 ALL_WATER_gw User interface 28/06/

13 Read in the WEAP Area Read in the Linkage and the MODFLOW model Summary of the inputs 28/06/

14 Objectives and constraints GA parameters Objectives weights Stopping Progression parameters 28/06/

15 Statistics of the optimal solutions Choose a well Proposed optimal solutions Evaluate and save optimal solutions 28/06/

16 Case study application 9 groundwater nodes, 16 demand sites, 24 Transmission links, Period: , Monthly time step, Optimization , 28/06/

17 WEAP-MODFLOW model for Zeuss Koutine 28/06/

18 Thousand Cubic Meter Water Demand (not including loss, reuse and DSM) Scenario: Reference, All months (12) Benguerdene Jerba Medenine Other Tataouine Zarzis 0 Jan Feb Apr Jun Jul Sep Nov Jan Mar May Jul Sep Nov Jan Mar May Jul Sep Nov Jan Mar May Jul Sep Nov Jan Mar May Jul Sep Figure. Monthly water demand for domestic and touristic demand sites nodes. 28/06/

19 Simulation results MODFLOW Cell Head Layer: 1, Row : 24, Column: 25 Reference changement climatique Jan Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec Dec /06/

20 Meter MODFLOW Cell Head Year: 2015, Scenario: Reference, Layer: 1, month: December Column Meter 37.5 to to to to to to to to to to to to to to to to Row /06/

21 Optimization results Figure. Overview of the DEMAND_FRACT.csv file produced by ALL_WATER_gw. 28/06/

22 U.S. Dollar Jan 2010 Mar May Jul 2010 Sep Nov Jan 2011 Mar May Jul 2011 Sep Nov Jan 2012 Mar May Average Cost of Water (per m^3) All months (12) Jul 2012 Sep Nov Jan 2013 Mar May Jul 2013 Sep Nov Jan 2014 Mar 2014 May 2014 Jul 2014 Sep Nov Jan 2015 Mar May Jul 2015 Sep Nov OPTIMIZATION Reference Figure. Optimized and Reference Average Cost of Water (KWh per m3). 28/06/

23 MODFLOW Cell Head Layer: 1, Row : 23, Column: 28 0 Jan Oct Aug Jun Apr Feb Dec Oct Aug Jun Apr Feb Dec Oct Aug Jun Apr Feb Dec Oct Aug Jun Apr Feb Dec Oct Aug Jun Apr Feb Dec Oct Aug Jun Apr Feb Dec Oct Aug Jun Apr OPTIMIZATION Reference MODFLOW Cell Head Layer: 1, Row : 23, Column: 30 0 Jan Oct Aug Jun Apr Feb Dec Oct Aug Jun Apr Feb Dec Oct Aug Jun Apr Feb Dec Oct Aug Jun Apr Feb Dec Oct Aug Jun Apr Feb Dec Oct Aug Jun Apr Feb Dec Oct Aug Jun Apr OPTIMIZATION Reference Figure. Optimized and Reference cells heads of two piezometers. 28/06/

24 Summary & Conclusions ALL_WATER_gw can be used in the framework of WEAP-MODFLOW or as stand alone tool. It considers the main objectives and constraints for groundwater management, For real groundwater, it was able to identify optimal management solutions, ALL_WATER_gw is a free tool. 28/06/

25 شكرا Merci Thank You Danke schön 28/06/