Tenth Gulf Water Conference, 22-24 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 inouiri@yahoo.fr
Outline Problematic, General conceptual model and objective, Problem formulation and resolution methodology, ALL_WATER_gw presentation, Case study application. 28/06/2012 2
Problematic How much water has to be taken from each well? 28/06/2012 3
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
Objective The main objective of this work was to develop a multi-objective optimization tool for groundwater management. 28/06/2012 5
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
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
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
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
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/2012 10
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/2012 11
ALL_WATER_gw User interface 28/06/2012 12
Read in the WEAP Area Read in the Linkage and the MODFLOW model Summary of the inputs 28/06/2012 13
Objectives and constraints GA parameters Objectives weights Stopping Progression parameters 28/06/2012 14
Statistics of the optimal solutions Choose a well Proposed optimal solutions Evaluate and save optimal solutions 28/06/2012 15
Case study application 9 groundwater nodes, 16 demand sites, 24 Transmission links, Period: 1983 2015, Monthly time step, Optimization 2013-2015, 28/06/2012 16
WEAP-MODFLOW model for Zeuss Koutine 28/06/2012 17
Thousand Cubic Meter 1 500 1 400 1 300 1 200 1 100 1 000 900 800 700 600 500 400 300 200 100 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 1982 1983 1984 1985 1986 Sep Nov 1987 1988 Jan Mar May 1990 1991 1992 Jul Sep Nov 1993 1994 1995 Jan Mar May 1997 1998 1999 Jul Sep Nov 2000 2001 2002 Jan Mar May 2004 2005 2006 Jul Sep Nov 2007 2008 2009 Jan Mar May 2011 2012 2013 Jul Sep 2014 2015 Figure. Monthly water demand for domestic and touristic demand sites nodes. 28/06/2012 18
Simulation results MODFLOW Cell Head Layer: 1, Row : 24, Column: 25 Reference changement climatique 50 40 30 20 10 0 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 1982 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 28/06/2012 19
Meter 50 44 48 46 42 4038 36 42 40 38 36 34 MODFLOW Cell Head Year: 2015, Scenario: Reference, Layer: 1, month: December 0 32 30 28 26 24 22 20 18 16 14 12 10 8 Column 35 30 25 20 15 10 5 Meter 37.5 to 40.0 35.0 to 37.5 32.5 to 35.0 30.0 to 32.5 27.5 to 30.0 25.0 to 27.5 22.5 to 25.0 20.0 to 22.5 17.5 to 20.0 15.0 to 17.5 12.5 to 15.0 10.0 to 12.5 7.5 to 10.0 5.0 to 7.5 2.5 to 5.0 0.0 to 2.5 34 32 30 28 26 24 Row 22 20 18 16 14 12 10 8 6 4 28/06/2012 20
Optimization results Figure. Overview of the DEMAND_FRACT.csv file produced by ALL_WATER_gw. 28/06/2012 21
U.S. Dollar 2.1 2.0 1.9 1.8 1.7 1.6 1.5 1.4 1.3 1.2 1.1 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 Jan 2010 Mar May 2010 2010 Jul 2010 Sep Nov 2010 2010 Jan 2011 Mar May 2011 2011 Jul 2011 Sep Nov 2011 2011 Jan 2012 Mar May 2012 2012 Average Cost of Water (per m^3) All months (12) Jul 2012 Sep Nov 2012 2012 Jan 2013 Mar May 2013 2013 Jul 2013 Sep Nov 2013 2013 Jan 2014 Mar 2014 May 2014 Jul 2014 Sep Nov 2014 2014 Jan 2015 Mar May 2015 2015 Jul 2015 Sep Nov 2015 2015 OPTIMIZATION Reference Figure. Optimized and Reference Average Cost of Water (KWh per m3). 28/06/2012 22
54 52 50 48 46 44 42 40 38 36 34 32 30 28 26 24 22 20 18 16 14 12 10 8 6 4 2 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 1982 1982 1983 1984 1985 1986 1986 1987 1988 1989 1990 1991 1991 1992 1993 1994 1995 1996 1996 1997 1998 1999 2000 2001 2001 2002 2003 2004 2005 2006 2006 2007 2008 2009 2010 2011 2011 2012 2013 2014 2015 OPTIMIZATION Reference 52 50 48 46 44 42 40 38 36 34 32 30 28 26 24 22 20 18 16 14 12 10 8 6 4 2 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 1982 1982 1983 1984 1985 1986 1986 1987 1988 1989 1990 1991 1991 1992 1993 1994 1995 1996 1996 1997 1998 1999 2000 2001 2001 2002 2003 2004 2005 2006 2006 2007 2008 2009 2010 2011 2011 2012 2013 2014 2015 OPTIMIZATION Reference Figure. Optimized and Reference cells heads of two piezometers. 28/06/2012 23
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/2012 24
شكرا Merci Thank You Danke schön 28/06/2012 25