Leading practices for assessing the integrity of confining strata: application to mining and coal-seam gas extraction

Similar documents
10. GEOTECHNICAL EXPLORATION PROGRAM

Aquitard Characterization The Legend of Indiana s Magic Clay Layer. Juliet Port, LPG #2214 July 2014

1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #5: Groundwater Flow Patterns. Local Flow System. Intermediate Flow System

Workshop on Aquifers, Aquitards & Mining Edited by: Steve Bouzalakos, UNSW Australia and Wendy Timms, UNSW Australia

Measurement, Monitoring and Verification (MMV)

Reservoir Geomechanics and Faults

11/22/2010. Groundwater in Unconsolidated Deposits. Alluvial (fluvial) deposits. - consist of gravel, sand, silt and clay

Surface Processes Focus on Mass Wasting (Chapter 10)

Water Inflow Issues above Longwall Panels

Determining In Situ Properties of Claystone Aquitards Using Pore Pressure Responses from Grouted-in Pressure Transducers

Groundwater. (x 1000 km 3 /y) Oceans Cover >70% of Surface. Groundwater and the. Hydrologic Cycle

Pit Slope Optimization Based on Hydrogeologic Inputs

SASKATCHEWAN STRATIGRAPHY GLACIAL EXAMPLE BOULDERS IN GLACIAL DEPOSITS

Recommendations for Injection and Storage Monitoring

Finding Large Capacity Groundwater Supplies for Irrigation

HYDROGEOLOGICAL PROPERTIES OF THE UG2 PYROXENITE AQUIFERS OF THE BUSHVELD COMPLEX

Assessing the Tier 2 Trigger for Fractured Sedimentary Bedrock Sites

Lecture 15: Subsidence

Accelerated gravity testing of aquitard core permeability and implications at formation and regional scale

Current challenges at CO 2 Sites

Chapter 12 Subsurface Exploration

Applied Geophysics for Environmental Site Characterization and Remediation

Gotechnical Investigations and Sampling

Groundwater. (x 1000 km 3 /y) Reservoirs. Oceans Cover >70% of Surface. Groundwater and the. Hydrologic Cycle

4.11 Groundwater model

PREDICTION OF ACID MINE DRAINAGE POTENTIAL FROM COAL MINES

Saltwater injection into a fractured aquifer: A density-coupled mass-transport model

CNSC Review of the Long-Term Safety Case for a Deep Geologic Repository

Prof. Stephen A. Nelson EENS 111. Groundwater

Azimuthal Resistivity to Characterize Fractures in a Glacial Till. Mark Boris, University of Saskatchewan Jim Merriam, University of Saskatchewan

Vertical hydraulic conductivity of a clayey-silt aquitard: accelerated fluid flow in a centrifuge permeameter compared with in situ conditions

Boreholes. Implementation. Boring. Boreholes may be excavated by one of these methods: 1. Auger Boring 2. Wash Boring 3.

1. Rock Mechanics and mining engineering

Delineation of Zones at Risk from Groundwater Inflows at an Underground Platinum Mine in South Africa

The Interaction of Reservoir Engineering and Geomechanics (a story)

IAEA SAFETY STANDARDS Geotechnical Aspects of Site Evaluation and Foundations in NPPs, NS-G-3.6

Blanca Payàs 1, Ryan Morris 2, Alain Arnaud 1, David Albiol 1, John Bradshaw 3, Javier Duro 1. Abstract

Site Investigation and Landfill Construction I

Predicting the soil-water characteristics of mine soils

ENCE 3610 Soil Mechanics. Site Exploration and Characterisation Field Exploration Methods

REPORT ON ESTIMATES OF GROUNDWATER INFLOW INTO AND DRAWDOWN AROUND THE PROPOSED OPEN PIT SOUTHDOWN IRON ORE PROJECT. Submitted to:

4D stress sensitivity of dry rock frame moduli: constraints from geomechanical integration

Geological Mapping using Geophysics

University of Pretoria. Matthys Dippenaar and Louis van Rooy Engineering Geology and Hydrogeology, Geology Department July 2017, Livingstone, Zambia

RADIONUCLIDE DIFFUSION IN GEOLOGICAL MEDIA

Before writing the discipline examinations, candidates must have passed, or have been exempted from, the Basic Studies Examinations.

Interpretation and application of hydrogeological concepts to commercial scale brine extraction projects

How to Interpret Mining Company Drill Reports & Announcements

DATA ACQUISITION METHODS FOR GROUNDWATER INVESTIGATION AND THE SITING OF WATER SUPPLY WELLS

A Risk-based Groundwater Modelling Study for Predicting Thermal Plume Migration from SAGD Well-pads

ENVIRONMENTAL EFFECTS OF GROUNDWATER WITHDRAWAL IN SOUTH NYÍRSÉG

Storage: Deep Monitoring and Verification

Appendix Z. Groundwater Assessment Supplementary Report. Document Name i Insert Month/Year

Prediction of subsoil subsidence caused by opencast mining

Geologic Considerations of Shallow SAGD Caprock; Seal Capacity, Seal Geometry and Seal Integrity, Athabasca Oilsands, Alberta Canada

Settlement characteristics of major infrastructures in Shanghai

Manual on Subsurface Investigations National Highway Institute Publication No. FHWA NHI Federal Highway Administration Washington, DC

GPS and GIS Assisted Radar Interferometry

b Blackwell Science Environmental Geology PRINCIPLES AND PRACTICE

- Cased-hole geophysical well logs - Airborne geophysics - Seismic data - LiDAR. GeoConvention 2015: New Horizons 1

Underground Risk Management Course Marina Del Rey, California November, Geotechnical Data Reports. Greg Raines, PE

Deep Borehole Disposal Performance Assessment and Criteria for Site Selection

FUNDAMENTALS OF ENGINEERING GEOLOGY

GEOTECHNICAL ENGINEERING INVESTIGATION HANDBOOK Second Edition

ractical Geomechanics for Oil & Gas Industry

Instructional Objectives

Geotechnical Models and Data Confidence in Mining Geotechnical Design

Cone Penetration Testing in Geotechnical Practice

Chapter 14: Groundwater. Fig 14.5b

The process of determining the layers of natural soil deposits that will underlie a proposed structure and their physical properties is generally

Methods of Interpreting Ground Stress Based on Underground Stress Measurements and Numerical Modelling

Chapter 7 Permeability and Seepage

Geotechnical modelling based on geophysical logging data

Geophysical Exploration in Water Resources Assessment. John Mundell, P.E., L.P.G., P.G. Ryan Brumbaugh, L.P.G. Mundell & Associates, Inc.

Origin and Evolution of Formation Waters in the West-Central Part of the Alberta Basin

Safety assessment for disposal of hazardous waste in abandoned underground mines

Understanding the Causes of Roof Control Problems on A Longwall Face from Shield Monitoring Data - a Case Study

Lines-of-Evidence Approach to the Evaluation of Stray Gas Incidents

CE 240 Soil Mechanics & Foundations Lecture 5.2. Permeability III (Das, Ch. 6) Summary Soil Index Properties (Das, Ch. 2-6)

The UK GeoEnergy Test Bed Ceri J Vincent British Geological Survey

CO 2 storage capacity and injectivity analysis through the integrated reservoir modelling

ENGINEERING GEOLOGY AND ROCK ENGINEERING ASPECTS OF OPERATION AND CLOSURE OF KBS-3

1 Water Beneath the Surface

5 IEAGHG CCS Summer School. Geological storage of carbon dioxide (a simple solution)

Geomechanics for reservoir and beyond Examples of faults impact on fluid migration. Laurent Langhi Team Leader August 2014

25th International Conference on Ground Control in Mining

MUDLOGGING, CORING, AND CASED HOLE LOGGING BASICS COPYRIGHT. Coring Operations Basics. By the end of this lesson, you will be able to:

the Quarrying Industry Dewatering and the Quarrying Industry the Quarrying Industry

Wisconsin s Hydrogeology: an overview

Underground Storage & Disposal - The Salt Concept. Thomas Brasser - GRS

*** ***! " " ) * % )!( & ' % # $. 0 1 %./ +, - 7 : %8% 9 ) 7 / ( * 7 : %8% 9 < ;14. " > /' ;-,=. / ١

Exploration, Drilling & Production

Groundwater Modeling for Flow Systems with Complex Geological and Hydrogeological Conditions

PART I Hot Dry Rock Geothermal Energy: History and Potential of the Newest and Largest Renewable Energy Resource

Land subsidence due to groundwater withdrawal in Hanoi, Vietnam

WESTCARB Regional Partnership

The Mine Geostress Testing Methods and Design

' International Institute for Land Reclamation and Improvement. 2 Groundwater Investigations. N.A. de Ridder'? 2.1 Introduction. 2.

Evaluation of the hydraulic gradient at an island for low-level nuclear waste disposal

Site Investigations and Geotechnical Risk For Underground Construction Greg Raines, PE

Transcription:

Leading practices for assessing the integrity of confining strata: application to mining and coal-seam gas extraction Wendy Timms 1,3, 4, Ian Acworth 2,3, Adam Hartland 2,3, David Laurence 1,4 1 School of Mining Engineering, University of New South Wales, Sydney, NSW 2052, Australia, w.timms@unsw.edu.au 2 School of Civil and Environmental Engineering, University of New South Wales 3 National Centre for Groundwater Research and Training and UNSW Connected Waters Initiative 4 Australian Centre for Sustainable Mining Practices, University of New South Wales Abstract Confining strata, or aquitards that act as low permeability seepage barriers, can limit potential impacts of depressurization or dewatering and migration of contaminants associated with mining and coal seam gas (CSG) development. Innovative geophysical and hydrochemical methods are recommended in a staged approach to best practice aquitard assessment including coupled numerical modelling. For example, tests of drill core in the new NCGRT geotechnical centrifuge have enabled relatively rapid measurement of permeability or vertical hydraulic conductivity (K v ) of clayey sediments and shale (K v 10-9 to <10-12 m/s). Stratigraphic and structural data and 3D modelling are also required to assess lateral continuity, and the likelihood of preferential leakage paths. Key words: confining strata, aquitard, seepage barrier, permeability, integrity Introduction Naturally occuring seepage barriers can limit potential impacts of depressurization or dewatering and migration of contaminants. Such low permeability strata are known as confining strata or aquitards. An ideal aquitard is a deep or thick material of low permeability, that is laterally continuous and without preferential flow paths (e.g. fracture networks). Plastic clay rich material can form an effective aquitard provided that there is no loss of hydraulic integrity due to faults or geological structures, and where they are resistant to stress fractures from drilling or subsidence. Potential contaminants can be effectively contained by an aquitard where solute fluxes are controlled by chemical gradients (diffusive transport), and less controlled by flow (advective transport), and geochemical reactions act to degrade or isolate salts, trace metal or organic constituents. The volume of groundwater flow to underground mine workings or during coal seam gas (CSG) production can have implications for safe and efficient extraction and for potential impacts on shallow aquifers, or nearby rivers, lakes and wetlands. For example, tunneling and coal mining beneath large water bodies with negligible inflows is practised where geological conditions are favourable. Extraction of water for coal seam gas (CSG) production is required to reduce the hydraulic head at the coal seam sufficiently to enable gas desorption. The feasible design of water extraction systems to achieve this low hydraulic head depends on many factors including the permeability of the coal seam, and the degree of vertical disconnection or connectivity through the overburden. Proponents of mining or CSG developments in areas of perceived high risk of vertical connectivity will be subject to consent conditions that require a mining method to be demonstrated with no adverse environmental impact. Community scrutiny of these developments has reached such an intensity that the Federal Government has appointed an expert committee to ensure that both coal mining and coal seam gas proposals will have minimal impact on water resources (see http://www.environment.gov.au/epbc/about/coal-seaminterim-committee.html). The objective of this paper is to discuss leading practice methods for assessing the hydraulic integrity of aquitards for underground mining and CSG operations. A staged approach is recommended for assessing aquitard integrity, depending on the degree of risk and the proximity of high value aquifers and surface waters. Examples of innovative methods are given, although detailed information on each method is beyond the scope of this overview paper. 1

Methods Table 1 outlines a proposed approach for leading practice assessment of the integrity of confining strata (aquitards) for underground mining and coal seam gas extraction. The appropriate level of practice for a given site (Level 1, 2 or 3) would depend on the risk of significant vertical leakage. If vertical leakage is significant, depressurization may not be effective, or environmental impacts could be detected and considered to be unacceptable. Table 1 Leading practice for assessing the integrity of confining strata (aquitards) for underground mining and coal seam gas extraction. Information & assessment Level 1 Level 2 Level 3 Geophysical surveys Drilling and strata sampling Appropriate airborne/remote sensing methods Grab samples from rotary mud drilling L1 + appropriate ground surveys eg. Seismic reflection Minimally disturbed core L2 + baseline InSAR and geodetic surveys, microseismic monitoring L2 + preservation of target core in sealed plastic Geophysical logging Qualitative logs (eg. Natural gamma) L1 + quantitative logs to identify strata Groundwater pressure Aquifers only L1 + data within confining strata Hydrochemistry Aquifers only, ph, L1 + stable isotope EC, temperature and tracers and aquitard major ions porewater chemistry from core, multiple sample sites/cores Permeability Geomechanical integrity Regional structural analysis Depressurization and/or subsidence Groundwater modelling approach Integrated assessments Aquifers only (AS 2368-1990) n/a (requires cores) Small scale geological mapping Estimated cone of depression extent over time, estimated subsidence/upsidence (valley bulge) as appropriate Aquifer layers, conductance factor for vertical flow, steady state calibration Groundwater focus L1 + bore permeability tests (packer, drill stem tests etc.) Index properties (eg. plasticity), rock quality designation Large scale geological mapping & 3D geological model Modelled estimates of depressurization or subsidence/upsidence, horizontal stress conditions Model multiple realisations of heterogeneity and structures, hydraulic values constrained by direct data, transient calibration and verification etc. Multi-disciplinary approaches using established methods L2 + logs for defining strata properties and defects L2 + continuous logging L2 + multiple sampling events, organic and radioisotope tracers, assess potential leaching of solid phases L2 + core matrix permeability testing (ASTM D6527 or ASTM D7664) L2 + Extensometers, rock mass strength & tensile strength of core L2 + groundwater data in each key structural area (eg. either side of fault), assess lateral continuity & subcropping L2 + models for variations over time (eg. Pore pressure recovery rates, principal stresses) constrained by direct data and observations L2 + flow and recharge boundaries constrained by direct data L2 + coupled processes for deformable, fractured or multiphase flow as appropriate L2 + research and development of improved/more efficient methods 2

Results and Discussion The concept of leading practice is simply the best possible way of conducting activities for a given area, such as the conceptual system in Figure 1. The key hydrogeological features of this coal bearing section are illustrated, including the role of confining strata, and possible seepage pathways. As new challenges emerge and new solutions are developed, it is important that leading practice be flexible in developing solutions that match site-specific requirements (Laurence, 2011). Leading practice is therefore as much about possible approaches (Table 1) as it is about a fixed set of practices or a particular technology. For example, new technologies such as the centrifuge permeameter can provide data on matrix permeability and transport properties of confining strata, as one part of an overall assessment approach. Level 1 information and assessment approaches are typical of approaches with limited emphasis on confining strata, while Level 3 provides an indication of more robust and innovative approaches that could be deployed for a high risk site. The practices in Table 1 are not intended to be prescriptive or exclusive. A key point of these suggested leading practices is that multi-disciplinary approaches that combine a number of independent and complimentary assessments would increase the confidence in the findings. Figure 1 Groundwater conditions in a coal bearing sequence (modified after Thomas 2002) Remote sensing and geophysical surveys are typically part of early geological assessments of potential coal mining and coal seam gas extraction. Table 1 suggests that Level 1 airborne and remote sensing could be enhanced at Level 2 with appropriate ground surveys, and at Level 3 with the addition of more advanced methods if there is a high risk of impacts due to ground disturbance. For example, radar interferometry is a well-known technique that can measure the ground movement of a wide-area, and which has been widely used for monitoring subsidence due to underground mining (Ng et al., 2010). The differential interferometric SAR (DInSAR) technique has the potential to precisely observe the ground displacement with an accuracy down to a few millimetres. To detect ground motion, satellite images over time are compared, and verified with available geodetic ground surveys. Methods are being developed by a UNSW team, as part of NCGRT research, to improve numerical models by assessing aquitard integrity at field sites, and with laboratory testing. Field site characterization has included geophysical techniques (gravity, resistivity, cross-hole seismic) to assess moisture and cracking status of shallow aquitards. Drilling projects to date have confirmed that coring, although more expensive and time consuming, produces far more accurate lithology logs than rotary mud methods and grab sampling that miss fine sediment, and permeable pathways that can occur within confining strata. For example, Timms and Hendry (2003) reported that the downwards migration of brine through a fractured till was very limited by a deeper unoxidized till with no visible fractures (Figure 2a). Geophysical logging, core porewater analysis and centrifuge permeameter testing were all part of the assessment of the integrity of the low permeability barrier at this potash mine site. 3

Many techniques for assessing confining strata were developed for glacial tills in North America, including for example, direct analysis of stable isotopes on core samples (Hendry and Wassenaar, 2010). Deuterium isotopes often have distinctive end-members (enriched or depleted) that provide unique tracers of mixing, and enable more detailed vertical profiles of porewater characteristics than is possible from a limited number of multi-level piezometers (Figure 2b). A Level 3 assessment with analysis of minimally disturbed core, could for example, identify strata that form the most effective hydraulic barrier. Without such detailed information, averaging of groundwater level data collected from aquifers above and below a confining strata could be misleading (Figure 2b). The utility of multiple organic tracers for assessing potential vertical connectivity is presented by Hartland et al. (2012). Dissolved organic matter (DOM) was characterised using liquid chromatography organic carbon detection (LC-OCD), fluorescence spectroscopy, and total organic carbon (TOC) concentration and stable carbon isotopes ( 13 C), supplemented by liquid chromatography tandem mass spectrometry (LC-MS/MS) analysis of anthropogenic tracers such as atrazine and pharmaceuticals. a) b) Figure 2 a) Geophysical logging to detect base of oxidized fractured zone in glacial till at a potash mine (Timms and Hendry, 2003) b) Conceptual schematic of information provided by Level 1-3 aquitard assessment The hydraulic conductivity of aquifers and aquitards is important for assessing whether or not vertical connectivity could be significant. A conceptual cross-section indicates that seepage can occur laterally and vertically (Figure 1). Differences in hydraulic conductivity measured at matrix scale (core samples), during bore tests (eg. aquifer pump test, or drill stem tests) and for regional scale models may be attributed to preferential flow. Flow pathways could occur for instance, via conductivity faults, while other faults filled with weathered clay act as flow barriers that compartmentalise a groundwater system (Figure 1). Tests of drill core in the new NCGRT geotechnical centrifuge using ASTM D7664-10 have enabled relatively rapid measurement of vertical hydraulic conductivity (Kv) of clayey sediments and shale compared with standard test methods (Table 2). The technical specifications of the Broadbent GT18 (875 RPM max, 2 m diam.) centrifuge, including the permeameter module (550 G max for 2 4.2 kg samples) has been established in a new facility at UNSW (Timms, 2012). The centrifuge permeameter enables large permeability datasets to be obtained in reasonable time for confining strata. Recent K v data of semi-consolidated smectite varied from 10-9 to 10-10 m/s depending on stress state, at depths of up to 40 m in the Gunnedah coalfields area. By contrast, shale can be 100 to 1000 times less permeable with test values to date ranging from 10-10 to <10-12 m/s (or <0.0000001 m/day). With future development of instrumentation in the NCGRT-UNSW facility, even lower K values could be measured, along with the variation of K value with moisture or saturation. 4

Table 2 Standard test methods for measurement of hydraulic conductivity - fine grained soils during lowering of groundwater level and wetting and drying (ASTM 2010, D 7664-10). Test Method Equipment Advantages Disadvantages Testing time for silty clays Rigid wall permeameter 1 Oedometer with outflow measurement 2 Flexible wall permeameter 1 Permeameter with outflow control 2 Stress control, volume change measurements Impedance of porous stone 1-2 weeks Centrifuge permeameter Centrifuge permeameter, instrumentation Fast testing time, best for wetting and drying hysteresis Equipment requirements <1 week 1 Axis translation method, 2 with high entry porous disc Findings so far have highlighted the sensitivity of K to moisture content and small fractures. Minimally disturbed core samples from depths greater than 100 m have been tested, although load restrictions in this centrifuge cannot match in situ lithostatic stresses at such depths. Nevertheless, the permeability values are consistent with larger permeability datasets for shales (Neuzil, 1994). In contrast to aquicludes (zero flow), the flow of groundwater through aquitards is typically very slow, but could be significant for large areas over long time periods if underlying aquifers remain depressurized (Timms, 2012). Effective zero flow conditions are possible when K and hydraulic gradient reduce linear flow velocity to a point where solute transport is dominated by diffusion (associated with a solute concentration gradient). Steady state flow conditions in the centrifuge also enable reactive transport experiments to quantify the degree of contaminant retardation under in situ conditions (Figure 3). A range of other tests are also possible for no flow and transient flow conditions in the centrifuge permeameter module (Table 3). For example, retardation of trace metals and analogues of radio-nuclides were defined for a natural glacial till seepage barrier over thousands of years during 3 months of centrifuge modelling (Timms et al. 2009). Analysis of both aqueous phase (core effluent) and solid phase (core dissection after testing) enable a very large range of retardation values to be measured for realistic flow conditions (Figure 3). Table 3 Centrifuge permeameter test setups Figure 3 Steady-state permeability and contaminant transport testing for low permeability core samples in a centrifuge permeameter compared with static batch tests. Test Setup No flow Steady flow Transient flow Centrifuge permeameter measurements Compressibility, specific storage, pore water extraction Permeability, reactive transport, effective porosity Permeability with variable moisture, water retention curve 5

Few studies to date have directly assessed the geomechanical integrity of sediment and rock aquitards, particularly hydraulic responses to varying effective stresses associated with depressurization. For example, the assessment of rock mass strength and tensile strength of strata is recommended for high risk sites (Table 1), along with regional structural analysis and mapping. While faults are often inferred to act as barriers to flow, detailed groundwater pressure data could provide evidence of flow barriers within sedimentary sequences. Such hydrogeological and geomechanical information are important to better constrain 3D modelling of seepage and depressurization processes. Discussion of coupled modelling and integrated assessment of confining strata integrity is beyond the scope of this paper, as is consideration of barriers to lateral flow around mine pits and engineered seepage barriers. However, it is noted that design and construction of seepage barriers is an option for mining operations to limit subsurface flows between pits or tailings and the surrounding catchment. Combined surface water and groundwater models have been used for example as part of assessing a seepage barrier at the edge of a pit in a fractured rock (Wasko et al., 2011). Conclusions This paper has outlined examples of leading practice assessment of confining strata that are important for safe and efficient mining operations that minimise environmental impacts. While not intended to be prescriptive, or comprehensive, a staged assessment approach has been recommended, with a suite of innovative tools now available for application to proposed coal seam gas or underground mining developments that are identified as high risk or are located near high value rivers or surface waters. Acknowledgements This work was funded by the Australian Research Council and National Water Commission, as part of National Centre for Groundwater Research and Training, Program 1B. Mark Whelan is thanked for designing centrifuge instrumentation systems for testing, and Dayna McGeeney is thanked for improving core test methods. References AS (1990). Testing pumping of water wells. Australian Standard AS2368-1990. ASTM, (2010). Standard test method for measurement of hydraulic conductivity of unsaturated soils. American Society for Testing and Materials International, D 7664-10. Hartland, A, Timms, W., Khan, S and J McDonald (2012). Application of multiple organic tracers for assessing fluxes through aquitards and potential connectivity with surface waters. IAH Congress, Niagara, Canada, September 2012. Hendry M.J. and Wassenaar, L.I. (2010). Inferring Heterogeneity in Aquitards Using High-Resolution δ 2 H and δ 18 O Profiles. Groundwater 47(5): 639 645 Laurence, D. (2011). A guide to leading practice sustainable development in mining. Australian Government, Department of Resources, Energy and Tourism. Neuzil C.E., (1994). How permeable are clays and shales? Water Resources Research, 30(2):145-150 Ng, A., Ge, L, Yan,Y., Li, X, Chang, H, Zhang,k. & Rizos, C, (2010). An approach for mapping accumulated mine subsidence using small stack SAR differential interferograms in the Southern Coalfield of New South Wales, Australia', Engineering Geology, 115 (1-2)1-15 Timms, W, Hendry, J., Muise J, and Kerrich, R. (2009). Coupling Centrifuge Modeling and Laser Ablation ICP- MS to determine contaminant retardation in clays. Environmental Science and Technology. 2009, 43, 1153 1159 Timms, W.A. and Hendry, M.J., (2003). Application of centrifuge modelling to assess long term brine migration in thick clay till, Saskatchewan, Canada., Australian Institute of Mining and Metallurgy Conference Water in Mining, Brisbane, October 13-15, 2003, pp. 363-372. Timms, W. (2012). Environmental time machine - the integrity of aquitards overlying coal seams. AusIMM Bulletin, April 2012 issue, pages 79-81. Thomas, L. (2002) Coal Geology. John Wiley & Sons, West Sussex, England. Wasko C.D., W.A. Timms, and B.M. Miller (2011). Groundwater Modelling in a Highly Seasonal, Fractured Rock Environment using FEFLOW. 6