The model needs to define and vertical and horizontal hydraulic conductivity for every element in the model mesh. These properties might be expected to vary with lithology or geology. An analysis of 577 hydraulic conductivity measurements from the found little to no correlation between lithology and hydraulic conductivity, nor between geologic layer and hydraulic conductivity (see Section 2.1.3 in companion product 2.1-2.2 for the Hunter subregion ()). However, there was a weak correlation with depth (see Figure 26 in ). Companion submethodology M07 for groundwater modelling () proposes the use of a simple parameterisation of hydrostratigraphy in (BA) groundwater model layers that treats the layers as homogeneous, but varies hydraulic properties with depth. In the absence of a good basis for varying hydraulic properties by lithology or geology, the parameterisation of the Hunter subregion groundwater model adopts this approach.
As described in Section 184.108.40.206.3.1, every point in the model domain can be defined in terms of a layer number (N = 1 to 11), a lithology (L = 1 to 7) and a depth, d, below the model topography. The horizontal hydraulic conductivity, Kh, the vertical hydraulic conductivity, Kv, and the porosity, Φ, are assumed to be of the form:
where a is the exponential decay coefficient which varies each parameter with depth.
Reference horizontal conductivities, porosities and decay parameters are defined in Table 6, and were obtained by appropriately averaging data from Figure 26 and Figure 31 of . It is not appropriate to directly use the raw conductivity data presented in Figure 26 in the groundwater model. This is because the measured data pertain to samples on the scale of centimetres (for lab measurements) to a few tens of metres (in-situ measurements). The regional groundwater model has a best resolution of 500 m, so some ‘upscaling’ of the raw data is necessary. Consider the problem of prescribing a suitable conductivity to a 500 m region of the groundwater model. Typically, there will be some regions of low conductivity within that region but there will also be regions of high conductivity, and water will flow preferentially through those highly conductive regions, almost entirely bypassing the regions of low conductivity. Therefore, the 500 m region will have an effective conductivity that is closer to the high conductivities shown in Figure 26. Of course, given the scatter of raw data in Figure 26, the conductivity of our 500 m region will be drawn from a of conductivity values, but the mean of this probability distribution will be the (arithmetic) mean of the raw data. A classic review of this subject may be found in . This arithmetic mean is shown in Figure 26 as a function of depth, and leads to the value shown in Table 6.
Less data are available to constrain vertical conductivity. typically occur as roughly horizontal layers of low conductivity in stratified geological systems such as the Hunter subregion. Such aquitards have little effect on horizontal water flow, since it flows relatively quickly through the surrounding , but have a greater effect on vertical flow, since the water must pass through the aquitard. This means that groundwater models typically use a vertical conductivity that is a small multiple of horizontal conductivity. A base multiplying factor of 0.1 is chosen in the Hunter groundwater model (and this is varied in the analysis).
The alluvium layer (L = 7) is assumed to be more porous and conduct water more rapidly than the interburden layers (L = 1 to 6), reflecting its unconsolidated nature. Its thickness is represented as one finite element, which means at any point the conductivity of the alluvium is a single number and does not vary with depth. However because the thickness of the alluvium varies across the subregion, and conductivity is evaluated at the centroid of a 3D finite element, the effect of the decay parameter means a 20 m thick alluvium will have conductivity of 0.78 m/day, while a 10 m thick alluvium will have conductivity of 0.88 m/day. These differences are negligible. Thus the decay for the alluvium layer is unimportant since it is so thin in comparison with the overall model thickness.
Horizontal and vertical conductivities are assumed to decay at the same rate with depth.
Table 6 Reference hydraulic property values and decay parameters
In many groundwater models, the decay parameters ah, av, aΦ, and the functions Kh(L,N), Kv(L,N) and Φ(L,N) would be calibrated by observation with water heads. Here, however, these are all varied in the uncertainty analysis as described in Section 220.127.116.11.2
Product Finalisation date
- 18.104.22.168 Methods
- 22.214.171.124 Review of existing models
- 126.96.36.199 Model development
- 188.8.131.52 Boundary and initial conditions
- 184.108.40.206 Implementation of the coal resource development pathway
- 220.127.116.11 Parameterisation
- 18.104.22.168 Observations and predictions
- 22.214.171.124 Uncertainty analysis
- 126.96.36.199 Limitations and conclusions
- Currency of scientific results
- Contributors to the Technical Programme
- About this technical product