A regional numerical model of the Galilee Basin – here referred to as the Galilee Basin hydrogeological (GBH) model – has been developed utilising data and interpretation undertaken for the by the Assessment team. The model was jointly funded by the Commonwealth Department of the Environment’s Office of Water Science and the Queensland Government. A consultant – HydroSimulations – was engaged to produce a calibrated and stress-tested transient numerical groundwater flow model capable of simulating the of proposed coal mining developments. The GBH model was built over a 2.5 month period. The following summaries are based on a detailed description of the GBH model, provided in .
Input data for the model was compiled as part of this assessment for the Galilee subregion. Companion product 2.1-2.2 for the Galilee subregion () provides detail on data used as input to the GBH. These data include: the geological model of the Galilee subregion, potentiometric surfaces for in the subregion, and volumes, among other . The GBH model was built utilising the new MODFLOW-USG platform and AlgoMesh, which uses flexible meshes to allow refinement around key features while maintaining a reasonable cell count for computational efficiency. The GBH model consists of an unstructured grid (except for mines) consisting of 1.93 million cells ranging in size from 150 m to 10 km and covers an area of 299,400 km2. A cell size of 165 m to 250 m is used around modelled mines. Modelling was carried out in three stages, resulting in steady-state, transient and predictive simulations. Mines simulated in the transient model are: China First, Kevin’s Corner, Carmichael, South Galilee, China Stone, Hyde Park and Alpha. These are mines that comprise the modelled , as outlined in companion product 2.3 for the Galilee subregion ().
The GBH model consists of 13 layers representing the major hydrostratigraphic units present within the subregion (Table 6). In many cases, individual model layers represent multiple stratigraphic units that are grouped based on similar hydraulic properties.
Table 6 Model layer thickness and corresponding hydrostratigraphic groupings (after Turvey et al., 2015)
Significant vertical hydraulic gradients exist through the upper Permian coal measures. For this reason the upper Permian coal measures were modelled as three layers (layers 10, 11 and 12). From top to bottom these sub-units are informally designated BC1, BC2 and BC3. The partitioning was done on the basis of approximately similar groundwater pressures existing in each sub-unit (Table 7). The upper Permian coal measures are designated as a partial aquifer because of low hydraulic conductivity (variable Kh values, but less than 0.2 m/day).
Table 7 Stratigraphic units and marker beds used to define subdivision of the upper Permian coal measures
Major rivers within the subregion are represented in the model using the MODFLOW River package and are assumed to be connected to the . Rivers modelled include: Carmichael River, Belyando River, Native Companion Creek, Alpha Creek, Aramac Creek, Thomson River, Flinders River, Diamantina River, Barcoo River, Warrego River, and Blackwater Creek.
Groundwater from 121 are simulated in the model using the MODFLOW River package (head-dependent boundary conditions) by setting the stage and bed elevation of the river cells to be equal, allowing one-way flow (discharge from groundwater) only. The River package allows for the specification of spring – discharge observation points where the fluxes from these river cells are functions of the calculated groundwater level and assigned conductance. Major springs represented include Doongmabulla, Mellaluka, Edgbaston, Coreena and Corinda springs. Based on limited data available from the Queensland springs database, the major springs were assigned flows greater than 100 m3/day; all other springs modelled were limited to less than 100 m3/day.
Long-term mean recharge rates used in the model are based on chloride mass balance estimates reported in companion product 2.1-2.2 for the Galilee subregion (). Estimated recharge rates vary from a low of 0.1 mm/year, for the Joe Joe Group, to a high of 5 mm/year for Cenozoic alluvial deposits.
It is considered that no lateral groundwater or outflow occurs along the eastern margin of the Galilee subregion. Consequently, no boundary condition is assigned to the eastern edge of all model layers. In contrast, general-head boundaries are applied along the western boundary of the model to the aquifer model layers 2, 4 and 6, representing the regional hydraulic gradient and communication with aquifers outside of the model domain.
Existing groundwater is represented in the model using the MODFLOW Well package. For steady-state simulations the mean recorded artesian discharge between 1983 and 2012 was applied to be consistent with the data used for calibration. For the historical transient model the actual annual discharge rate values were used, while predictive model runs used a constant flow rate of 895 ML/day, which represents estimated total daily extraction for the year 2012.
noted that discharge estimates provided may be inaccurate and that for some bores there is about the source aquifer. As a result it is possible that well stresses are misapplied in the model. To compensate, variable bore discharge rates have been applied during model simulation to prevent incorrectly assigned discharge rates from drying out cells in low conductivity layers, causing model instability.
Mine for the seven mines included in the predictive model is simulated using the MODFLOW Drain package. The package integrates the proposed mine plan location, timing, depth and method of extraction. Open-cut mine pits are assumed to be the full thickness of the relevant target coal seams. Drain cells are progressively added to the transient model to simulate extraction of mine panels over time.
Longwall mine induced fracturing is simulated using the MODFLOW-USG Time-Variant Materials package. Hydraulic parameters are changed with time in the and overlying fracture zones during extraction of each longwall panel.
18.104.22.168.6.1 Steady-state calibration
Steady-state calibration was assessed against water levels corrected for temperature and salinity to ensure comparability with MODFLOW predictions which do not take account of density variation.
Calibration statistics for the steady-state simulation are shown in Table 8. Predictions within ±10 m of target levels were distributed evenly across the model domain, however a tendency remains to under predict water levels in the west and over predict water levels in the east. An area of significant over prediction (>50 m) occurs within layers 11 and 12 at the south-east of the model. Target residuals at the southern group of mines (Kevin’s Corner, Alpha, China First and South Galilee) are relatively good, with an overprediction of approximately 13 m being the most significant discrepancy. However, at Carmichael, particularly north of the Carmichael River, significant model overpredictions of up to 38 m occur within the coal seam layers. No data were available for calibration at either China Stone or Hyde Park mines at the time of model development.
Table 8 Steady-state calibration statistics (after Turvey et al., 2015)
includes a comparison of simulated potentiometric surface contours derived from the steady-state calibration and contours interpreted from observed water levels for companion product 2.1-2.2 for the Galilee subregion (). found that there is general agreement between the pattern of the Eromanga Basin modelled surfaces and interpreted surfaces. A comparison of model layers 1 and 2 and interpreted surfaces for the Cenozoic and Winton and Mackunda shows differences in elevation generally less than 20 m and similar inferred flow directions. A comparison between the model layer 4 surface and the interpreted surface for the Hooray Sandstone showed modelled elevations 10 to 30 m below the interpreted surface and similar inferred flow directions. These properties are the same for the comparison between the model layer 6 surface and the interpreted surface for the Hutton Sandstone.
The discrepancies between the modelled Galilee Basin units (model layers 8, 10, 11 and 12) and interpreted surfaces are more variable with a general pattern of under prediction of modelled water levels over most of the model domain and over prediction of modelled water levels toward the eastern boundary of the subregion where outcrop. However, inferred flow directions are generally similar for both modelled and interpreted surfaces.
The general underestimation of modelled water levels in the location of some is likely to account for predicted spring flows that are less than observed flows. The total modelled spring flow under estimates reported flow by around one order of magnitude.
The water balance for the steady-state simulation indicates that just under half of the to is from direct rainfall recharge, with an approximately equal contribution from river leakage and along model boundaries accounting for the remaining portion of inflow. Losses from the system in order of magnitude are: evapotranspiration; groundwater ; to rivers and springs; followed by groundwater outflow along model boundaries. All but two rivers in the model (Thomson and Diamantina rivers) act as net losing systems overall.
22.214.171.124.6.2 Transient calibration
Manual transient calibration was carried out for the period January 1998 to December 2012. Due to the limited data available for transient calibration, steady-state Kh and Kv values were not varied and only variation in storage parameters (Ss and Sy) was allowed. High specific storage (Ss) values were required to prevent excessive due to over the calibration period. For the transient simulation, rates were varied using a series of factors for the period 1983 to 2011 derived from the Australian Water Resources Assessment landscape model (AWRA-L) in order to scale steady-state recharge values to derive a transient recharge . Recharge estimates for the period 2012 to 2014 were derived using the Penman-Grindley recharge model ().
Calibration statistics for the transient simulation are shown in Table 9.
Table 9 Transient model calibration statistics (after Turvey et al., 2015)
The transient model mass balance indicates that the majority of water that enters the model does so via a combination of river leakage and rainfall recharge. Regional groundwater flow into the model, via general-head boundaries, is approximately double that of outflow. Over the calibration period there was a net loss in storage suggesting a depleting resource.
Two predictive model runs were simulated and run out to the year 2220 in order to allow for sufficient simulation of the water level recovery phase.
The first predictive model run for the GBH model represents a baseline which consists of existing use without groundwater due to coal mining (). For , the is defined as a future that includes all coal mines and coal seam gas (CSG) fields that are commercially producing as of December 2012 (see Section 2.3.4 in companion product 2.3 () for further detail). As of December 2012, there were no operational coal mines or CSG fields in the (). Hence, the first predictive run of the GBH model represents the coal resource development baseline defined for this assessment.
The second predictive model run includes the seven coal mines that comprise the modelled (see Section 126.96.36.199 in companion product 2.3 () for detail) and represents the cumulative of these potential coal mining operations. Additional stresses to those imposed on the baseline predictive model run include: a progressive excavation and associated of the seven open-cut and longwall mines; progressive development of fracture zones above underground mines and associated changes in ; and progressive placement of fill in open-cut mines to end of mine life using modified hydraulic parameters and .
The GBH model is a complex regional model produced in a very short time period and within a limited budget. Upon completion, a number of limitations and issues have been identified by the contractors – HydroSimulations – and reviewers from Queensland Government, Office of Water Science and relevant project and discipline leads from the Bioregional Assessment Technical Programme. The limitations noted by various reviewers include:
- The ability of the model to simulate to the unconfined Cenozoic – model layer 1 – is limited. Predictive model runs did not show any impact to model layer 1, particularly in the vicinity of open-cut mine pits. It would be expected that and of Cenozoic units would occur adjacent to mine pits which penetrate and expose the Cenozoic units. suggested that certain options within the MODFLOW-USG code do not allow the model to represent negative pressures (i.e. dry cells representing unsaturated conditions). This means that an artificial head is maintained at the base of each layer, which may result in under prediction of the magnitude of drawdown within each cell. This will have had particular implication for impacts predicted in model layer 1, essentially resulting in water being held at the base of model layer 1 for the entire model simulation and resulting in negligible reported drawdown within this layer.
- Combining several hydrogeological units into a single model layer means that accurate representation of real-world groundwater levels for all members of the GAB and Galilee Basin units cannot be produced. This is a particular issue for the upper Permian coal measure layers (model layers 10, 11 and 12), where bulk hydraulic conductivity parameters applied to the model representing coal seam layers are higher than those of the interburden layers (that are not included). Additionally, no vertical hydraulic conductivity variation exists within the bulked layers. As discussed in Section 188.8.131.52.2.4, combining hydrogeological layers together can be justified when conservatively estimating drawdowns. The lumping of hydrogeological layers, however, can result in biased parameter estimates if local groundwater level observations are used to constrain the aggregated hydraulic parameters. Another potential issue is that mine dewatering rates will be over-estimated as the aggregated unit will be dewatered in the model, rather than just the mined interval. Such an overestimate can be considered conservative and therefore justified. It will contribute to discrepancies between estimated pumping rates between models, further complicating comparisons between model results.
- The representation of losing and gaining streams in the GBH model is suited to predict the hydrological change due to coal resource development. However, in a deterministic aquifer simulator such as the GBH model, incorrectly specified – groundwater boundary conditions is a large source of . This may bias parameter estimates and therefore bias predictions. Debate exists about the representation of the Thomson, Alice, Carmichael and Belyando rivers as losing streams. Initial assessments done for companion product 2.1-2.2 for the Galilee subregion () suggests these streams could be predominantly gaining. Section 184.108.40.206.1 identifies the large potential of improving the characterisation of the connection status of the river network to reduce predictive uncertainty.
- Queensland Government reviewers noted issues with some of the input data, in particular the boundary definitions for some model layers (e.g. the Rewan Group) and the accuracy of some flow rates. There is a possibility that over of groundwater from the model due to erroneous input well flow rate data has resulted in calibrated hydraulic conductivities and storage properties being too high. Anomalously high hydraulic conductivities and storage properties are likely to result in increased mine and greater resultant drawdowns.
- Limited data were available for calibration; most data collected were of single-point observations over a period of 30 years. Therefore, reliability of the data is unknown as long-term trends were not able to be analysed. Section 220.127.116.11.1 highlights the need for additional field data collection.
- Limited formal analysis has been carried out due to the short model development time frame and agreed scope for the project. It is therefore difficult to clearly indicate which parameters are controlling calibration, inflows and predicted impacts.
- The GBH model input data does not include more recent groundwater data that became available after the release of the groundwater modelling report for China Stone Coal Project. Extra data in would refine the results of the steady-state and transient models in the vicinity of China Stone and Carmichael coal projects.
Product Finalisation date
- 18.104.22.168 Methods
- 22.214.171.124 Review of existing models
- 126.96.36.199.1 Alpha and Kevin's Corner model review
- 188.8.131.52.2 Carmichael model review
- 184.108.40.206.3 China First model review
- 220.127.116.11.4 China Stone model review
- 18.104.22.168.5 South Galilee model review
- 22.214.171.124.6 Galilee Basin hydrogeological model review
- 126.96.36.199.7 Suitability of existing groundwater models
- 188.8.131.52 Model development
- 184.108.40.206 Boundary and initial conditions
- 220.127.116.11 Implementation of the coal resource development pathway
- 18.104.22.168 Parameterisation
- 22.214.171.124 Observations and predictions
- 126.96.36.199 Uncertainty analysis
- 188.8.131.52 Limitations and conclusions
- Currency of scientific results
- Contributors to the Technical Programme
- About this technical product