Revision of the Numerical Model for the Lower Hutt Groundwater Zone APRIL 2003
Total Page:16
File Type:pdf, Size:1020Kb
APRIL 2003 Revision of the numerical model for the Lower Hutt groundwater zone APRIL 2003 Revision of the numerical model for the Lower Hutt groundwater zone Prepared for: Greater Wellington – The Regional Council ? phreatos GROUNDWATER CONSULTING Contents 1. Introduction...............................................................................................................1 2. Previous modelling .................................................................................................2 3. Data sources.............................................................................................................3 3.1 Geological data...........................................................................................................3 3.2 Harbour bathymetry...................................................................................................4 3.3 Groundwater levels....................................................................................................4 3.4 River stage data .........................................................................................................5 3.5 Groundwater usage...................................................................................................5 4. Hydrogeology...........................................................................................................5 4.1 The Lower Hutt groundwater zone (LHGZ) ...........................................................5 4.2 Hydrostratigraphy.......................................................................................................6 4.2.1 Taita alluvium..............................................................................................................6 4.2.2 Petone marine beds and melling peat....................................................................7 4.2.3 Waiwhetu artesian gravels .......................................................................................7 4.2.4 Wilford shell bed.........................................................................................................9 4.2.5 Moera gravels .............................................................................................................9 4.2.6 Deep strata..................................................................................................................9 4.2.7 Greywacke ................................................................................................................10 4.3 Hydraulic properties.................................................................................................10 4.3.1 Taita alluvium............................................................................................................10 4.3.2 Petone marine beds/melling peat..........................................................................11 4.3.3 Upper Waiwhetu gravels.........................................................................................11 4.3.4 Lower Waiwhetu aquifer.........................................................................................12 4.3.5 Wilford shell beds.....................................................................................................13 4.3.6 Moera basal gravels ................................................................................................13 4.3.7 Deep strata................................................................................................................13 4.3.8 Greywacke ................................................................................................................14 4.4 Recharge...................................................................................................................14 4.4.1 River recharge..........................................................................................................14 4.4.2 Rainfall recharge......................................................................................................15 4.5 Groundwater flows and aquifer discharge ...........................................................16 5. Numerical model design......................................................................................17 5.1 Model code................................................................................................................17 5.2 Finite difference grid design...................................................................................17 5.3 Model boundaries ....................................................................................................19 5.4 River simulation........................................................................................................19 5.4.1 River bed elevation..................................................................................................19 5.4.2 River stage................................................................................................................21 5.4.3 River bed conductance ...........................................................................................21 5.4.4 River spreadsheet....................................................................................................22 5.4.5 Hutt River south of KGB and Waiwhetu Stream.................................................22 5.5 Discharge simulation...............................................................................................23 6. Model calibration ...................................................................................................23 6.1 Procedure..................................................................................................................23 6.2 Steady state calibration...........................................................................................24 6.2.1 Input parameters......................................................................................................25 6.2.2 Steady state model calibration results..................................................................26 6.2.3 Steady state calibration verification.......................................................................27 6.2.4 Steady state model sensitivity analysis................................................................28 6.3 Transient flow calibration........................................................................................30 6.3.1 1996-1997 transient calibration.............................................................................30 6.3.2 14-year calibration...................................................................................................33 6.3.3 2000-2001 calibration..............................................................................................34 7. Summary and Conclusions ................................................................................35 References.............................................................................................................................37 REVISION OF THE NUMERICAL MODEL FOR THE LOWER HUTT GROUNDWATER ZONE 1. Introduction The Lower Hutt – Port Nicholson sedimentary basin (Figure 1) contains a regionally important groundwater resource which supplies up to 40% of the water demand for the greater Wellington region. The basin contains several confined artesian, semi-confined and unconfined gravel aquifer units which collectively constitute a layered aquifer system known as the Lower Hutt Groundwater Zone (LHGZ). Currently, major municipal abstraction takes place only from the most productive Upper Waiwhetu Artesian Gravels at an average rate of approximately 60 ML/day but peaking at 100 ML/day during the summer months. Management of the groundwater resource is reliant upon a robust evaluation of the sustainable yield in conjunction with a strategic groundwater level and water quality monitoring system to ensure that saline coastal waters do not invade the aquifers. Given increasing abstraction demands on the Lower Hutt Groundwater Zone, there is a need to review the safe yields for groundwater system, and to review management safeguards to minimise the risk of saline water intrusion. Recent advances in the geological and hydrogeological understanding of the groundwater system prompted the Wellington Regional Council to revise the existing numerical model (HAM1) built some eight years ago (Reynolds, 1993). To avoid inherent limitations of the old model, which focussed only upon the Waiwhetu Artesian Gravels and was limited by poor definition of the system recharge and discharge mechanisms, construction of an entirely new model has been required (HAM2). The HAM2 has a less complex the spatial zonation of hydraulic properties than the previous model, and a more detailed layer structure based upon revised geological interpretation of the sedimentary sequence. The model boundaries have been re-designed allowing aquifer discharge processes to be more appropriately modelled, and the simulation of recharge through the Hutt River bed has been refined. These changes have resulted in an improved calibration which has undergone sensitivity and uncertainty analysis using the automated parameter estimation routine, PEST (Watermark Computing, 1998). The objectives of re-building the Hutt Aquifer Model can be summarised as follows: · Facilitate the re-assessment and refinement of the sustainable and safe yield of the Waiwhetu Aquifer with greater confidence. · Enable an assessment of the feasibility and advantages of abstracting groundwater from lower stratigraphic levels