History, Geology, In Situ Stress Pattern, Gas Content and Permeability of Coal Seam Gas Basins in Australia: A Review

: Coal seam gas (CSG), also known as coalbed methane (CBM), is an important source of gas supply to the liquefied natural gas (LNG) exporting facilities in eastern Australia and to the Australian domestic market. In late 2018, Australia became the largest exporter of LNG in the world. 29% of the country’s LNG nameplate capacity is in three east coast facilities that are supplied primarily by coal seam gas. Six geological basins including Bowen, Sydney, Gunnedah, Surat, Cooper and Gloucester host the majority of CSG resources in Australia. The Bowen and Surat basins contain an estimated 40Tcf of CSG whereas other basins contain relatively minor accumulations. In the Cooper Basin of South Australia, thick and laterally extensive Permian deep coal seams (>2 km) are currently underdeveloped resources. Since 2013, gas production exclusively from deep coal seams has been tested as a single add ‐ on fracture stimulation in vertical well completions across the Cooper Basin. The rates and reserves achieved since 2013 demonstrate a robust statistical distribution (>130 hydraulic fracture stages), the mean of which, is economically viable. The geological characteristics including coal rank, thickness and hydrogeology as well as the present ‐ day stress pattern create favourable conditions for CSG production. Detailed analyses of high ‐ resolution borehole image log data reveal that there are major perturbations in maximum horizontal stress (S Hmax ) orientation, both spatially and with depth in Australian CSG basins, which is critical in hydraulic fracture stimulation and geomechanical modelling. Within a basin, significant variability in gas content and permeability may be observed with depth. The major reasons for such variabilities are coal rank, sealing capacity of overlying formations, measurement methods, thermal effects of magmatic intrusions, geological structures and stress regime. Field studies in Australia show permeability may enhance throughout depletion in CSG fields and the functional form of permeability versus reservoir pressure is exponential, consistent with observations in North American CSG fields.


Introduction
Coal seam gas is an unconventional gas resource and is an important source of gas supply in countries with large coal reserves such as Australia.While coal seam gas development in the United States has declined primarily due to the shale revolution, it is a thriving industry in Australia.The demand for clean energy in East Asia has resulted in substantial growth of the CSG industry in Australia.In 2018, Australia became the world's largest exporter of LNG: 29% of the nameplate capacity being from Australia's three east coast plants supplied primarily from coal seam gas resources [1]. Figure 1 shows the major CSG basins, wells and associated infrastructure including oil and gas pipelines in Australia.The CSG industry is now in a new phase where the aim is to increase the export and domestic supplies and that requires significant ongoing drilling and completion activities.There are also new exploration and drilling campaigns, particularly in Queensland.In South Australia, operators are attempting to better understand unique reservoir behaviour of deep coal seams in the Cooper Basin and unlock gas resources by optimizing hydraulic fracturing and drilling multi-fractured horizontal wells [2][3][4][5][6][7].This paper reviews the geology, in situ stress pattern, gas content and permeability of major coal seam gas basins in Australia through a review of various sources including published articles, personal communications, government reports and datasets.This paper is built on previous review papers such as a review on coalbed methane by Moore [8], relative permeability in coal by Zhang, Feng [9], modelling permeability by Pan and Connell [10] and production data analysis by Clarkson [11].An early review paper by Miyazaki [12] discussed exploration and development of coal seam gas resources in Australia.Later, a review paper published by Towler, Firouzi [13] reviewed various aspects of CSG from history and drilling to decommissioning and regulatory compliances in Queensland.In contrast, this paper integrates coal geology and in situ stress pattern and investigates variations in gas content and permeability, both spatially and with depth, in major CSG basins in Australia.Furthermore, permeability change throughout depletion is discussed in the Bowen Basin.This allows more accurate evaluation of reserves and field development scenarios in Australian CSG basins.The other important aspect of this review paper is the investigation of deep coal seams of the Cooper Basin in South Australia as a potential gas resource.
This paper classifies the history of CSG in Australia into six key phases including pre-1975, early exploration, early production, political intervention, CSG to LNG and post-2016.The history of development of deep coal seams in the Cooper Basin is reviewed separately.The geological framework of the major basins is also discussed.The in situ stress analysis and maps presented are based on information collected from CSG reservoirs and analysed by Rajabi, Tingay [14].The information about in situ stresses can have a significant contribution to future development of CSG fields by improving hydraulic fracturing designs and mitigating wellbore stability issues.Finally, reservoir properties including gas content and permeability are discussed in the Bowen, Surat, Sydney, Gunnedah and Cooper basins and field examples are presented.

History of Coal Seam Gas in Australia
Figure 2 compares CSG production in Australia with CSG production in the US.It shows that the US production peaks in 2008 at 1966 BCF while during the same year, Australia production is only 150 BCF.The US production declines afterward and in 2016, Australia production surpasses the US production.The annual CSG production (see Figure 2) and number of wells drilled each year (see Figure 3) demonstrate six key phases in history of CSG in Australia.Although some CSG drilling has occurred in nearly all Australian mainland states, it is dominated by Queensland and New South Wales.The six phases were dominantly controlled by commerciality and overcoming technical difficulties while at the same time, working closely with local communities.

Phase 1: Pre-1975
Since first being discovered by Europeans in 1797 at Coalcliff in the Illawarra of New South Wales, coal mining has been a major industry for Australia.In 2017, Australia was the world's largest exporter of coal with 32% of world trade coming from approximately 30 major mine operations, primarily in the states of Queensland and New South Wales [17].The exported coal is primarily Permian in age and is generally high-quality bituminous coal, including both thermal and coking coal [18].
In the past, methane, desorbed from coal seams, has been seen a major hazard with gas explosions in underground coal mines causing major loss of life [19,20].To overcome this concern, many of the underground, as well as high wall open-cut mines in Australia are degasified by drilling directional in-seam degasification wells [19,21].Over time, some of this gas has been used for power generation or even domestic use.For example, Miyazaki [12] documents what was perhaps the first coal seam gas well in Australia being drilled under Sydney Harbour in the Balmain Colliery during 1944.

Phase 2: 1975-1996 (Early Exploration)
The history of modern coal seam gas exploration dates from 1976 when Houston Oil and Minerals drilled two wells targeting the Permian coals of the Bowen Basin in Queensland [22].Other USA-based oil and gas companies entered the scene during the 1980′s and 1990′s hoping to leverage their experience, primarily from the San Juan Basin which had just started commercial development after initial exploration having started as early as 1950 [23].This phase focused on the Permian coals of the Bowen Basin in Queensland and the Sydney Basin in New South Wales.The cost of drilling and completion prevented economic production during this phase.By 1995, approximately 280 wells had been drilled, 55% in Queensland and 45% in New South Wales.

Phase 3: 1996-2005 (Early Production)
CSG was first produced commercially from the Permian coals of the Bowen Basin at Dawson River, near the township of Moura, in 1996.In Queensland, this development was quickly followed by other Permian coals at Fairview in 1998 and Peat/Scotia in 2002 [24].In New South Wales, the Camden Gas Project with its Rosalind Park Gas Plant in Menangle commenced production in 2001.This project currently produces 5% of New South Wales gas from 144 wells.One interesting insight is that these early production projects are all on the crests of major anticlines where permeabilities are improved by the structural setting.Several of the best production wells in these early projects demonstrate what could be "free gas" production, similar to conventional gas extraction, requiring little to no dewatering.During this phase, the potential of the Jurassic coals (Walloon Coal Measures) of the Surat Basin in Queensland was tested by some early exploration wells around Dalby and Chinchilla.1280 wells had been drilled in total by 2005 with 75% in Queensland and 25% in New South Wales.

Phase 4-2005 to 2011 (Political Intervention)
Recognizing the decreasing supply of conventional gas to the east coast market, and the possible impact of CSG, the Queensland Government introduced the Electricity (13% Gas Scheme) Amendment Bill 2004 in January 2005.The bill encouraged development of new gas supply by requiring electricity retailers to supply at least 13% of their electricity from gas-fired or renewable energy sources.This requirement was increased to 15% in 2010 and 18% by 2020 [25].This bill acted as a catalyst to increase commercial production from the Surat Basin with projects commencing in 2006 at Tipton West, Kogan and Berwyndale and ramped up rapidly.CSG production in Queensland surpassed conventional gas production in the state in 2007 and the period finished in 2011 with Surat Basin CSG production surpassing Bowen Basin CSG production.The Global Financial Crisis of 2008 had a slight impact on drilling and production but was only a minor influence.
A very different story emerged in New South Wales during this phase.Exploration had commenced in the Gunnedah, Clarence-Morton and Gloucester basins of northern New South Wales but no projects moved into commercial production.The exploration licenses in Clarence-Morton and Gloucester basins were bought back by the New South Wales government during this Phase and Phase 5 for approximately $27 million to stop CSG exploration and production in those regions.By 2011 over 4600 CSG wells had been drilled in total with 80% being in Queensland and only 20% in New South Wales.

Phase 5: 2011-2016 (CSG to LNG)
With increasing additions to reserves, CSG operators began searching for a large market for their gas.This occurred at a time when East Asian nations were searching for gas to fuel their domestic growth plans.These two forces led to the unfolding of the idea of converting CSG into LNG for export.Final investment decisions (FID) were made for these projects in 2010-2012.All three projects, which each consisted of a two-train LNG liquefaction plant, intensified drilling campaigns and gathering and compression projects, were brought to fruition during this phase (Table 1).The two-train LNG liquefaction plants were all constructed next door to each other on Curtis Island in Gladstone Harbour.This phase marks a major change and maturation of the Australian gas industry landscape with two of the projects being operated by local companies (Santos and Origin)-a first for these Australian companies in operating large-scale LNG projects [13].As well as the massive infrastructure builds, this period of ramp-up to fill the LNG facilities led to the drilling and completion of some 6000 wells in Queensland.The low level of activity in New South Wales continued in this period with only 100 wells drilled.Hence, Queensland drilled 98% of all wells in this phase.Unlike conventional LNG projects where the development wells are drilled and simply produce through the life of the project, in CSG to LNG projects the gas stream is filled by continuous drilling of development wells.The industry has now entered the post LNG start-up phase where the requirement to continue drilling and completing thousands of wells to keep the LNG (export) and domestic gas markets serviced will be the key focus.The linking of LNG price to the international oil market has also given rise to an increase in the domestic gas price.These economic factors have facilitated the expansion of the CSG industry in Queensland and it is expected to continue for the life of the LNG projects and beyond.This means new development projects are continually coming on stream and are likely to do so for many years to come.However, as of 2019 two of the three LNG facilities were operating well below capacity and EnergyQuest has suggested that this is in part because of the projects original viability being based on estimates of Proved and Probable (2P) reserves rather than Proven (1P) reserves [26].2P is the standard practice for making investment decisions by the CSG industry.One of the basic assumptions at the time was high permeability discovered around the Undulla Nose in the Fairview and Spring Gully fields would be ubiquitous in Surat and Bowen basins.
However, this assumption turned out to be overly optimistic.The various companies reduced many of the original 2P reserves, on which the LNG projects were based, but by this time, all three LNG facilities were almost completed.Figure 4 shows the 2P reserves in the Surat and Bowen Basins since 2005, and the 2P reserves estimates are at the highest in 2014 to 2016 and then they decline afterward.This is partly because the companies reduced many of their original 2P reserves that were overly optimistic.The opposite occurred in New South Wales, where political and environmental activism has led to no new CSG projects coming on stream.The sole New South Wales CSG project at Camden is now moving into its wind down phase and is likely to cease production in 2023.In September 2020, the Independent Planning Commission (IPC) conditionally approved the Narrabri Gas Project.It is estimated that the Narrabri Gas Project has the potential to supply half of the New South Wales's natural gas demand [28].Geoscience Australia [18] currently estimates Australia's CSG reserves (2P) at 43 Tcf and identifies an additional 32 Tcf of contingent resources.92.5% of the reserves lie within the Surat and Bowen basins of Queensland.

Development of Deep Coal Seam Gas in the Cooper Basin, South Australia
Early assessment of deep coal seams in the Cooper Basin was conducted by analysis of gas shows in mud log records of wells intersecting deep coal seams [29], and the first desorption canister test was conducted in Dorodillo-2 in 1998 [30].In 2007, for the first time, Patchawarra coal seams were hydraulically fractured in Moomba 77 and a gas flow rate of ~0.10 mmscf/day (2832  /day) was reported.Figure 5 shows some of dedicated unconventional gas play wells in the Cooper Basin that were drilled to test productivity of deep coal seams and compares those with deepest producing coals in the United States and China.The deepest producing coal in the United States is from the Cameo Coal zone in the Piceance Basin (8136-8142 ft.) with a reported production rate of 0.15 mmscf/day [31].Hydraulically fractured dry Permian coals of the Cooper Basin show economic flowrates greater than those recorded for Cameo Coal [3].In 2013, Roswell-1 was drilled and hydraulically fractured in the VC50 coal (Patchawarra coal) using cross-linked gel fracturing fluid and a higher proppant concentration than used in previous trials.This conductive fracture design resulted in a peak flow rate of 0.40 mmscf/day.The flow results after this design change were encouraging and kick-started a year-on-year add-frac campaign in the Cooper Basin.This program involves targeting >5 m thick coals for fracture stimulation in as many conventional gas development wells across the Cooper Basin as possible.Repeat production logging as a means of monitoring of production performance of deep coal seams has been critical to demonstrating long term (>5 years) production, from these stimulated coal seams.To date, more than 130 add-fracs have been placed in deep coal seams of the Cooper Basin in gas development wells with overall success but variable productivity and is currently an economically viable reserve for target development since 2017 [32].
Washington-1, an exploration well drilled in 2015, was hydraulically fractured in the Toolachee, Epsilon and Patchawarra coal seams.Production logging demonstrated that gas flowed from all zones, with the deepest at 3491m (11,199 ft.); and, as such, is currently the deepest producing coal in the world [3].Two further dedicated vertical wells were drilled in 2018, as a follow-up to Washington-1.Casimir-1 intersected nine coals meeting the thickness criteria but is to be fracture stimulated.The second well in the 2018 campaign, Coonatie-25, accessed both coal and conventional tight sand.Eight coal seams were hydraulically fractured with the VC50 (mid-Patchawarra coal) delivering the highest rate on test of 0.48 mmscf/day.Coonatie-25 has been on production since March 2018.
The Southern Cooper Basin Gas Project, located 100 km south of Moomba, targets deep coal seams in the Patchawarra Formation at the depth of the 2000 m.In late 2017, the Klebb pilot that includes four wells was claimed as a technical success, but the wells did not achieve commercial gas rates [33].In 2018, Jaws-1, a horizontal multi-fractured well, was planned to drill in the Vu Upper Coal, however the well encountered wellbore stability issues while drilling and hence plugged and side-tracked and drilling continued below the Vu Upper coal for the purpose of better stability.Indirect vertical fracture completion (IVFC) technique was utilized to complete the well in the Vu Upper coal.Seven vertical fractures were placed in the Jaws-1, six of which initiated outside of the target formation [2].In December 2020, Jaws-1 was shut down because gas production did not reach to commercial rates and the Klebb pilot project will be ceased in 2021 [34].

Geological Framework of Coal Deposits and Coal Seam Gas Resources in Australia
Major accumulations of coal in Australia are present in Palaeozoic, Mesozoic and Cenozoic strata [35,36].Along with coal seam gas, some of these coals have long been mined as fuel for power plants and steel production.Geological factors such as coal rank and hydrogeology are important in making favourable conditions for gas production [8,[37][38][39][40][41][42][43][44][45].These, in turn, are controlled by the geologic history of each basin, including burial, uplift, tectonic deformation, and variations in heat flow through time.In this section, six basins with CSG resources are described.

Bowen Basin
The Permian-Triassic-Bowen Basin in Queensland is one of the world's major coal basins [46][47][48].It is the northern sector of the larger Sydney-Gunnedah-Bowen Basin, which is a composite structural basin that stretches from southern New South Wales to central Queensland (see Figure 1) [49].The basin has up to 10 km of sedimentary fill and has a complex tectonic and structural history [50][51][52].It began as a back arc basin in the early Permian and was followed by a phase of thermal subsidence in the early late Permian [53].In the late Permian, a foreland basin developed together with extensive coal accumulations [54,55].The basin was later subject to a major compressive deformation and uplift episode in the Middle-Late Triassic and intruded by igneous bodies in the Cretaceous, which influenced the preservation, distribution, rank and gas content of coals [49].The basin has significant CSG resources in addition to large reserves of metallurgical and thermal coal [13,22,40,52,53,[56][57][58].CSG is produced in the Moranbah and Moura regions (northern Bowen Basin), in the Injune area (western Bowen Basin) and in the Wandoan area (eastern Bowen Basin) [24,40,53,59,60].Three stratigraphic units are major targets for CSG production in the Bowen Basin [53] including the Moranbah Coal Measures/German Creek Formation at the Moranbah Gas Project; the Baralaba Coal Measures at the Moura Gas Project, and the Bandanna Formation at the Fairview and Spring Gully Projects [13,46,57,61,62].A summary of these stratigraphic units is provided in Table 2. [22,40,49,53,57,63,64].

Stratigraphic
Unit The most significant CSG project in the Bowen Basin is the Fairview Field, associated with the delta-top/fluvial upper Permian Bandanna Formation [53].The field is located on an anticline associated with the Comet Ridge, which enhances the permeability of the coal beds to >50 mD [53].In this locality, the Bandanna Formation is between 500-880 m and contains three coal seams that have a total net thickness of 10 m [65].
The thick and laterally continuous nature of the coal seams have contributed to the viability of CSG plays in the Bowen Basin.These coals have a wide range of ranks, vitrinite content and gas content but vitrinite reflectance is generally >0.9% and vitrinite content averages ~70% [12,[66][67][68].Variations in rank are attributed to the basin's complex tectonic and burial history, with gas sweet spots associated with large anticlinal structures, including the Burunga Anticline, where there is an increase in the permeability of coals (up to 100 mD) [40].The majority of CSG in the Bowen Basin is predominantly of a thermogenic origin (associated with deep burial depths) with a small biogenic component [13].

Sydney-Gunnedah Basins
The Late Carboniferous to Triassic Sydney-Gunnedah basins of New South Wales is the southern part of the Sydney-Gunnedah-Bowen Basin (see Figure 1) [36].The Sydney and Gunnedah basins share a similar complex tectonic history and contain up to 6 km of sedimentary infill [69].The basins evolved in several phases from the late Carboniferous to the Middle Triassic [70].The origins of the basins are uncertain with crustal extension [71] and volcanic rift models [72] proposed.A foreland basin phase associated with prevalent volcanic activity then ensued during the Permian [70].The basin was subject to a phase of compressional tectonics in the Late Permian [70].In the late Triassic, sedimentation ended in the basin due to an episode of major compressive deformation [73].Further episodes of deformation occurred in the Late Cretaceous and the Tertiary [70].Similar to the Bowen Basin, the Sydney and Gunnedah basins host significant reserves of metallurgical and thermal coal [74] but CSG activity occurs only on a limited scale in the region.CSG wells are predominately located in the Camden Field in the southern Sydney Basin [75] with some located on prospects in the Hunter Coalfield in the northern Sydney Basin [76].The Narrabri Gas Project in the Gunnedah Basin is thought to be highly prospective [77].Two stratigraphic units are targets for CSG production in the Sydney-Gunnedah basins and are summarised in Table 3.The southern Sydney Basin has a cumulative reservoir thickness (>20 m of coal) and thin interburden between productive horizons (<200 m) [70].Gas consists of up to 95% methane and is of thermogenic and biogenic origin but on a local scale has been influenced by magmatism [83].There is a lack of productive aquifers in the area and a well-developed cleat system enhances permeability (1-30 mD) [70].
In the Gunnedah Basin, cumulative coal thickness is greater than 80 m [70].Igneous intrusions play a significant part in the viability of CSG plays.The intrusions take the form of Late Triassic-Early Cretaceous alkali basalt and trachyte associations, and a Tertiary flood basalt association [84].Sills and dykes typically have a detrimental effect, such as dispelling coal seam gases, reducing permeability, and replacing methane with carbon dioxide [77,81].However, intrusions may also locally increase coal rank into the optimal thermogenic gas generation window as well as acting as impermeable seals to gas flow [77].Seams adjacent to the intrusions have vitrinite reflectance values as high as 11% Rv max [77].

Surat Basin
The Triassic-Cretaceous Surat Basin overlies the southern part of the Bowen Basin and contains Australia's largest CSG resource [85].The Surat Basin is a large intracratonic basin with up to 2.5 km of sedimentary infill (see Figure 1) [53].Three mechanisms have been proposed to explain the tectonic evolution of the basin including thermal subsidence [86], dynamic platform tilting [87] and intraplate rifting [88].The basin was uplifted during the Late Jurassic and again in the mid-Cretaceous, leading to the erosion of the eastern part of the basin [89,90].The major CSG producing interval is the Middle to Upper Jurassic Walloon Coal Measures of the northern Surat Basin.Walloon Coal Measures are divided into five distinct formations including the sandstone dominated Durabilla Formation, Taroom Coal Measures, relatively coal barren Tangalooma sandstone, Upper Juandah Coal Measures and Lower Juandah Coal Measures [91][92][93][94].The depositional environment of the Walloon Coal Measures has been investigated in the literature and different environments such as fluvio-lacustrine, composite meandering to anastomosing fluvial system and waterlogged alluvial plain have been suggested [85,91,[95][96][97][98].
Coal beds are thin: between 1 cm and 30 cm thick which have limited lateral continuity with individual beds rarely traced in the subsurface for more than 5 km [99].Vitrinite content is between 62 and 90%, liptinite content is between 8 and 33% and inertinite content is ~1% [100].The ash content varies substantially with a minimum of 5.1% and moisture contents up to 20.7 wt% [100].The majority of gas production is of a biogenic origin [101].

Cooper Basin
The intracratonic Late Carboniferous-Triassic Cooper Basin in South Australia and Queensland is Australia's largest onshore oil and gas province (see Figure 1) and has been producing conventional resources for over 50 years [102].Cooling of heat flow following granite emplacement in the Carboniferous led to rapid basin subsidence creating the Cooper Basin and later infilled by 2.5 km of sediments [103].The basin experienced three episodes of uplift during the Permian caused by plate margin stresses in the east [104].An intraplate east-west compressional regime brought an end to sedimentation in the basin during the Triassic [105].The basin is overlain by the Jurassic-Cretaceous Eromanga Basin and the Cenozoic Lake Eyre Basin [103].The stratigraphic unit with the most CSG potential in the basin is the Lower Permian Patchawarra Formation which was deposited in a fluvio-lacustrine setting [104,106].The unit is 250-500 m thick with coal seams between 0.3 and 60 m in thickness [103,106].Some seams cover an area of up to 1500 km 2 [106].The coals consist predominately of inertinite (45-86%) [103].Vitrinite content of the coals is mostly between 3% and 20% with liptinite content mostly between 1% and 15% [107].Coals are bituminous in rank (0.8% to 1.8%Ro) with low moisture contents (average 3.33 wt%) and an average ash content of 3.66% [108].The VC50 (5-23 m thick), Vm3 (19.5-24 m thick) and the Vu (17-45 m thick) coals are the most prospective for potential CSG activities in the basin [109].The VC50 coal is the thickest in the basin (5-23 m), with high lateral continuity (>8km) and significant micro-porosity (4% to 9%) [4,109].Other coals with CSG potential in the basin are in the Permian Epsilon (seams up to 29 m thick) and Daralingie formations (seams up to 3 m thick) [108].

Gloucester Basin
The early Permian Gloucester Basin is a small remnant inter-arc basin located in New South Wales (see Figure 1) [110].The basin has a complex geological history starting with early to late Permian extension [111].A late Permian episode of compression led to basin inversion [111].The Waukivory Pilot Project (now closed) and the Gloucester Gas Project (which reached the development phase before being cancelled) were the only CSG projects to have commenced in the basin [111].
The Lower Permian Gloucester Coal Measures, deposited in a deltaic environment, is the primary target for CSG production [111].As described by Hughes [110] and Ward, Bocking [112], the formation is up to 2 km thick with coal seams between 5 and 30 m in thickness.The coals are characterised by considerable lateral splitting over short distances [113].The coals are composed of between 45-50% vitrinite with high amounts of inertinite [114].

In situ Stress Pattern in the Australian Coal Seam Gas Basins
The in situ stress in Australian CSG reservoirs is complex and highly variable and hence it plays an important role in gas production [81,[115][116][117][118][119][120].In this section, the in situ stress pattern of the Australian coal seam gas basins are reviewed followed by their applications in coal seam gas exploration and production.
Rajabi, Tingay [14] have developed the in situ stress map of Australia with significant stress information from CSG basins.The regional pattern of SHmax orientation in the Australian continent can be explained by four regional trends.As indicated in Figure 6, the SHmax orientation in south-eastern Australia is NW-SE and rotates to approximately ENE-WSW in most of eastern Australia (including coal seam gas basins).The SHmax in the northern and north-eastern Australia, including northern part of the Bowen Basin, is NE-SW.A predominant E-W SHmax orientation exists in many parts of Western Australia and South Australia (including the Cooper-Eromanga Basins) [14].Several studies have been conducted on the in situ stress pattern of eastern Australia with particular emphasis on the CSG basins [115][116][117][119][120][121].Statistical analyses of stress data resulted in the classification of nine stress provinces in the Australian CSG basins (Table 4) [14,122].As shown in Figure 7 and Table 4, the mean SHmax orientation in northern part of the Bowen Basin is 032°, which is consistent with regional SHmax orientation in most part of the northern Australia.Southern part of the Bowen Basin and the entire Surat Basin in Queensland have a regional SHmax orientation of 065°.The mean SHmax orientation in New South Wales CSG basins has a regional NE-SW trend but there is significant variability similar to other basins across Australia.For example, the mean SHmax orientation in the Clarence-Moreton Basin is 069° which rotates to 057° and 034° in the Gunnedah and Gloucester basins, respectively.The Sydney Basin in New South Wales, which is the southern part of the Sydney-Gunnedah-Bowen Basin, shows two regional SHmax trend including 041° and 065° in the southern and northern part of the basin, respectively.In contrast with all CSG basins in eastern Australia that show a regional trend of NE-SW orientation for the SHmax, in the Cooper, Eromanga and Galilee basins in central Australia are characterised by a regional E-W orientation (Figure 7 and Table 4).
In addition to stress variability from basin to basin, large stress perturbations, at basin scale, have been reported in each of the CSG basins.These SHmax variabilities can be seen on the maps (Figure 7) and also are represented by high standard deviations in mean SHmax orientation of each stress provinces (Table 4).Basin scale rotation of stresses in CSG basins has been explained by major geological structures and complex tectonic deformation of each basin.This complex tectonic history in each basin has resulted in extensive faulting, folding and juxtaposition of different formations [53], which are responsible for basin scale rotation of stress.The stress variations in the Bowen, Surat, Gunnedah, Clarence-Moreton, Sydney and Gloucester basins are because of the presence of geological structures such as faults, lithological contrasts and basement structures that are documented in the literature (see Figure 7) [14,81,115,116,119,120,122].
Analysis of high-resolution wellbore image log data in CSG wells have provided the opportunity to further investigate the in situ stress pattern at small scales (i.e., wellbore scale; 1-1000 m).The localised stress variations that occur in the presence of small-scale geological features and localised elastic rock property contrasts play a critical role in subsurface fluid flow in coal seams [81,117,123].In the detailed analysis of borehole image logs in the Gunnedah Basin, Rajabi, Tingay [121] observed that geological structures, including faults, fracture and magmatic intrusions, strongly perturbed stresses in the region.Following this study, Salmachi, Rajabi [81] investigated the effect of magmatic intrusions on reservoir characteristics of the Hoskissons coal seams in the Gunnedah Basin by the integration of borehole image log interpretation and drill stem test data analysis.In one of the studied wells (Lake Goran 1) with the thickest intrusion (126 m), despite the presence of open fractures, very low permeability was measured.In situ stress analysis in Lake Goran 1 shows large SHmax perturbation with depth (standard deviation of ~40°).It was interpreted that large stress perturbation along this well may have impaired fracture connectivity resulting in low permeability [81].In situ stress pattern of major coal seam gas basins in eastern and central Australia and the smoothed stress pattern (orange lines) showing the regional orientation of maximum horizontal stress [14].Different symbols show the method of measurement for stress data including focal mechanism solution of earthquakes, borehole breakouts, drilling induced fractures, over coring, and hydraulic fracturing measurements (see [14] for the details of each method).Different colours indicate different tectonic stress regime (NF, normal; SS, strike-slip; TF, thrust; U, unknown).The length of the lines for each stress indicator illustrates the quality of stress orientation based on the World Stress Map ranking criteria [124].[14].A-C refers to available stress data records in each basin based on the World Stress Map quality ranking criteria [124].

Province
Number In addition to the SHmax orientation, the analyses of stress magnitudes have also shown changes in stress regime spatially and with depth [14,115,120,122,125].For example, shallow over-coring and hydraulic fracturing measurements in the Sydney Basin in New South Wales and Bowen Basin in Queensland show a prevailing reverse stress regime with some changes to the strike-slip regime [14,126].Recent analysis of wellbore data in CSG basins in the Bowen-Surat [120] and Cooper [125] basins have highlighted the changes of stress magnitude with depth.This variability of stress magnitudes with depth has also been predicted by 3D geomechanical-numerical model of Australia.There is a prevailing reverse stress regime at shallower depths of the coal seam gas basins of Australia that changes to the strike-slip regime and in some locations to normal stress regime [122].

Applications of In Situ Stresses in Coal Seam Gas Exploration and Production
It has been demonstrated that the in situ stresses significantly control the fluid flow pattern in the fractured reservoirs because those fractures which oriented parallel (or at ~30°) to the intermediate or maximum principal stress are more susceptible to tensile or shear failure [127].In the Andersonian stress regime classification, the SHmax is always the intermediate or maximum principal stress and, therefore, it is generally observed that fractures that are sub-parallel to the orientation of SHmax typically facilitate fluid flow and hydrocarbon extraction [127].Coal seam gas reservoirs generally have low matrix permeability and coal cleats play a critical role in permeability and fluid flow in these reservoirs [96,128,129].Therefore, the variable in situ stress pattern can significantly control the fluid flow within the cleat system.Mukherjee, Rajabi [117], investigated the interaction of subsurface fractures, orientation of SHmax and permeability data in 33 CSG wells penetrated into the Walloon Coal Measures, in the eastern part of the Surat Basin, Queensland.The results revealed that there are high SHmax perturbations, from well to well, due to geological structures affecting permeability.Interestingly, high permeability zones are located where SHmax rotates from the regional trend [117].
The stability of wells is also influenced by the geomechanical parameters such as orientation and magnitudes of in situ stresses.Drilling deviated and horizontal wells is challenging where stress pattern is complex [119,130].The state of in situ stresses has been also explained as an influencing parameter in hydraulic fracture design in unconventional reservoirs such as coal seam gas and shale gas reservoirs [119][120][121]129,131].This issue is particularly important for CSG basins in Australia because of significant stress perturbations in the major coal seam gas basins.These variable stress patterns have caused complex fracturing stimulations during hydraulic fracturing tests in the CSG basins [116,120].For example, Flottmann, Brooke-Barnett [116], investigated the relationship between in situ stress, natural fractures and propagation of hydraulic fractures in the Surat Basin.In particular, they reported significant variations in both stress orientation and stress regime, which resulted in the propagation of complex hydraulic fractures.Analysis of micro-seismic and tilt-meter results by Flottmann, Brooke-Barnett [116] revealed planar hydraulic fracture propagation in areas with high differential stresses (in reverse stress regime).In areas with low differential stresses, both stresses and preexisting natural fractures control the propagation of hydraulic fractures resulting in more complex hydraulic fracture patterns.Moreover, in another study conducted by Rahman, Suarez [118], many unsuccessful hydraulic fracturing treatments, due to the presence of complex in situ stress fields, in the Australian CSG fields have been discussed.
The variable stress state has numerous implications for CSG exploration and production in eastern Australia [37,132].Therefore, it is necessary to carry out a comprehensive analysis of present-day stress to determine the regional and local-scale stress perturbations for further decisions such as horizontal drilling and hydraulic fracture design.

Gas content and Permeability in Major Australian CSG Basins
Gas content, gas saturation and permeability are critical reservoir attributes that influence the economics of any CSG project [8].In Australia, the procedure used to measure gas content follows the Australian Standard (AS 3980:2016) "determination of gas content of coal and carbonaceous material-direct desorption method" (historical editions: AS 3980-1999, AS 3980-1991).Total gas content is the summation of three components; lost gas ( ), desorbed gas ( ) and residual gas ( ).Lost gas volume depends on several factors including how quickly a sample is retrieved and placed within the desorption canister, coal diffusion coefficient and diffusion distances (cleat spacing), drilling fluid and free gas volume [20,133].Gas content is heterogeneous spatially and even vertically within an individual coal seam in Australia.This is consistent with observations in North American coal seams reported by Pashin [134], Pashin [135] and Scott [44].The primary reasons for these variabilities are interrelated factors including coal rank and depth, thermal effect of magmatic intrusions, hydrogeology and sealing capability of overlying formations [44,77,81,136,137].
High gas production may not be achieved even from extremely high gas content coals if permeability is very low [44].Hence, it is important to investigate gas content and permeability mutually.A network of connecting fractures known as cleats as well as natural fractures determine coal permeability [96,[138][139][140][141].Depositional environment and tectonics influence permeability distribution.Coal rank, composition, ash content and mineralization are also important controls on permeability [142][143][144].Coal permeability is dynamic and can change throughout depletion because of cumulative impact of two competing factors, compaction and matrix shrinkage [145][146][147][148][149][150][151][152].Magnitude and functional form of permeability change are important in CSG reservoir engineering, reserve estimation and production forecasting.In addition, permeability is an important factor in choosing a suitable completion method [153].In this section, gas content and permeability of major Australian CSG basins are discussed.

Bowen Basin
Figure 8 shows gas desorption data and coal rank in 28 cored CSG wells in the Bowen Basin.Coals are sub-bituminous to anthracite in rank and there is a weak correlation between gas content and depth.Gas content of low rank coals including sub-bituminous and high volatile bituminous does not exceed 10 m 3 /ton (daf) (353 scf/ton) in the Bowen Basin (with a few exceptions).Low volatile bituminous and anthracite are the dominant ranks and coals could be heat affected (see Figure 8).Analysis of gas composition in these core wells reveals that 86% of gas samples have methane concentration (mole percentage) higher than 95%.
Figure 8. Gas content versus depth measured in 28 CSG core wells in the Bowen Basin.Coal rank is also shown, and heat affected coals are highlighted (data sourced from Queensland Government [154]).
Gas contents of the Bandanna Coals in the Fairview Field are shown separately in Figure 9 (gas desorption data are obtained from 27 observation wells drilled across the field).Although the general trend of increase in gas content with depth is evident, there is variability in gas content within the field.One of the reasons for abnormally low gas contents, reported in some observation wells, could be linked to reservoir depletion.High permeability (>100 mD) and extensive lateral continuity of the Bandanna Coals can result in excellent wellbore connectivity across the field.Hence, gas content measurements in the observation wells, drilled recently among existing wells, could be lower than initial gas content values because of depletion by offset wells.

Bandanna Coal
In the Bowen Basin, coal permeability is low in Moranbah and Moura areas (<10 mD) and hence, in seam drilling technique is vital to promote gas production.On the other hand, permeability in Fairview and Spring Gully is high and cavity completion has been a successful completion technique [153].The cavity completions in the Fairview are reported to be more successful in more depleted areas and this could be related to permeability enhancement throughout depletion [153].Several studies in the literature that reported magnitude and/or functional form of permeability change throughout depletion in the Bowen Basin are listed in Table 5. Permeability increase has been reported in the Bandanna Coals and the functional form of increase is exponential [155] similar to Fruitland Coals of the San Juan Basin, the United States, reported by Gierhart, Clarkson [156].Yarmohammadtooski, Salmachi [65] and later Clarkson and Salmachi [160] derived functional form of permeability change (see Figure 10) for a CSG well in the Fairview Field by analysis of production data above and below desorption pressure.Figure 10 shows there is a slight decline in permeability above desorption pressure (during the depressurization phase) because of stress dependency of permeability.It is followed by an increase and then permeability remains unchanged.Production data analysis works well when there are strong constraints on flowing bottom-hole pressures and reservoir pressures with the latter being more important.When such constraints exist, trends of permeability change with reservoir pressures are very comparable to pressure transient analysis methods.Excellent results (exponential trend) were reported when pressure observation wells, offset to producers, are used to constrain reservoir pressure [161,162].Change in permeability can be investigated by time-lapse pressure transient analysis that includes multiple in situ permeability measurements during depletion.Mazumder, Scott [157] investigated permeability change in 37 horizontal producers (single lateral or chevron geometries) in the Moranbah gas field by performing time-lapse PTA on pressure recovery data during shut-in periods.Although a small number of deviant wells reported, the general trend shows permeability increase by 2 to 50 times as reservoir pressure declines [157].
Salmachi and Barkla [155] replicated pressure build up tests in three CSG wells in the Fairview Field to construct a time-lapse PTA approach to study magnitude and functional form of permeability change in the Bandanna Coals.Two of the study wells (Well B and C) followed similar trends in permeability increase and Well A, intersected lower porosity coals, exhibited permeability enhancement by more than one order of magnitude (see Figure 11).Because of limited number of permeability measurements, functional form of change was not determined using the time-lapse PTA.However, rate transient analysis (RTA) of Well B indicates that functional form of permeability enhancement is exponential and a permeability modulus of −0.00678  (0.00098335 ) was obtained.The permeability values obtained by the RTA technique are in good agreement with results of pressure transient analysis (See Figure 12).Functional form of permeability change versus reservoir pressure depends on many factors including cleat compressibility, coal geomechanical properties and matrix shrinkage [7,45,146,149,151,155,163,164].Coal failure and fine production, when occurs, can result in decrease in permeability and subsequently affect the shape of permeability change curve [165][166][167][168][169][170][171][172][173].The impact of coal failure on permeability change in Australian CSG basins is less investigated and requires further research.

Sydney-Gunnedah Basins
Thomson, Thomson [174] divided the Sydney Basin into four distinctive zones to classify gas content and composition.From surface to the depth of 100 m, carbon dioxide (in very low gas content) is the dominant coal seam gas (Zone 1).In Zone 2 (100 m to 250 m) coal seams are rich in methane with biogenic signature.From 250 m to 600 m (Zone 3), gas is a mixture of thermogenic and biogenic methane and carbon dioxide with magmatic origin.In Zone 4 (>600 m) thermogenic methane is the dominant gas with some other hydrocarbons (e.g., ethane and propane) [174].The gas content in the southern Sydney Basin ranges from less than 1 m 3 /ton to 23 m 3 /ton (daf) with methane and carbon dioxide as dominant components [175].
In an integrated study conducted by Faiz, Stalker [175] information including gas content and composition, coal rank, gas production rate and geological structures were analysed in the Camden area, southern Sydney Basin.Although there is minor variation in coal composition and rank between high and low CSG producers, gases from highly productive wells are dry and have strong biogenic signature in comparison with ethanerich gasses that are produced from low producing wells.It was concluded that meteoric waters that contain bacteria have reached to high permeable areas, bacteria have metabolized wet gases, and secondary biogenic methane has been generated.Because of this, coals in high permeable areas have higher gas saturation (gas contents closer to sorption capacity) resulting in CSG wells to produce gas more quickly [175].
The target formations in the Camden Gas Project are Bulli and Balgownie coal seams with an average thickness of 4 and 2 m respectively.Permeability in the Bulli coal seam is in the range of 1-30 mD and gas content is in the range of 6-12 m 3 /ton [113,176].In Balgownie coal seam, permeability is less than 1 mD [176].
In the Gunnedah Basin, significant exploration activities including gas content measurements for drilled cores, seismic surveys, drill stem tests and pilot gas production have been executed for the Narrabri Gas Project [113].Figure 13 shows the relationship between gas content and depth in the Gunnedah Basin.Salmachi, Rajabi [81] showed that, in some of the studied wells across the Gunnedah Basin, magmatic intrusions had a constructive impact on reservoir characteristics of the coal seams by increasing coal rank and gas content.It is also observed that gas composition is primarily methane in coals even at immediate contact with intrusions [81].These findings are comparable with those reported by Gurba and Ward [136] and Gurba and Weber [77].Salmachi, Rajabi [81] integrated borehole image logs and DST data to study presence, openness and connectivity of cleats/natural fractures in the Hoskissons Coals in the Gunnedah Basin.The results from 14 wells across the basin indicate that coal permeability can vary from 1091 mD down to almost zero and permeability decreases with depth.It was also observed that natural fractures other than cleats play a significant role in fluid flow.Figure 14 provides a summary of interpretation of cleats and natural fractures in six of the studied wells in the Gunnedah Basin.High permeable coal intervals have well developed cleat and natural fracture systems (see Glasserton-1 and Dewhurst-19 in Figure 14).Borehole image log interpretation in Kahlua-1 indicates that although coal is poorly cleated, the open fractures within the coal interval result in permeability of 376 mD.This highlights the important contribution of natural fractures to coal permeability.

Surat Basin
Coals in the Surat Basin are less thermally matured (0.42-0.64%Ro) compared to coals in the Bowen Basin (0.6-3.5%Ro) and hence they generally have lower gas content [24].Figure 15 shows gas content versus depth in the Surat Basin based on gas desorption data collected from 86 wells [154].This dataset shows that gas content is lower than 10 m 3 /ton (daf) (353 scf/ton), except for a few measurements (see Figure 15).This observation is similar to that of the Bowen Basin where low rank coals show gas contents lower than 10 m 3 /ton (daf) (353 scf/ton).The gas content generally increases with depth, but data are scattered.The coals are sub-bituminous and high volatile bituminous in rank.Methane concentration (mole percentage) higher than 95% were recorded in the majority (96%) of the gas samples.

Sub-bituminous High Volatile Bituminous
Hamilton, Esterle [101] investigated spatial trends in gas content versus depth in the eastern part of the Surat Basin using core data from 69 CSG wells.Within the studied wells, three trends for gas content were identified; gas content either (1) increase; (2) increase and then decrease; or (3) decrease with depth.The majority of the studied wells fell in the category of 'increase and then decrease' and the inflection point is around the Tangalooma sandstone [101].This is consistent with the fact that gas content in the Lower Juandah Coal Measures is generally higher than gas content in the Upper Juandah and Taroom Coal Measures [177].
Coal seams of the Walloon Coal Measures are thin and does not continue laterally similar to those in the Bowen Basin [24,99].Although thickness and gas content in the Surat Basin are lower, the key to success in the Surat Basin has been the higher permeability.Permeability of the Walloon Coal Measures ranges from 5-5000 mD [178].Structural deformation is an important mechanism responsible for enhancement of coal permeability in the Surat Basin [177].Undulla nose and Kogan anticline are two important structural features that have higher permeability and host important CSG fields in the Surat Basin.Coal composition is another factor that controls fracture density within coal seams and, hence, permeability [94,179].Core analysis in the Surat Basin indicates that cleat and fracture networks are largely related to bright vitrinite-rich coal layers.Lower Juandah Coal Measures that contain higher vitrinite compared to Upper Juandah Coal Measures and Taroom Coal Measures, have higher permeability [177].
Figure 16 shows in situ permeability versus depth for the Taroom Coal Measures and Upper and Lower Juandah Coal Measures measured by pressure transient analysis method.The general trend shows decline with depth however, there is significant variability in permeability.For example, in the Lower Juandah Coal Measures, at the depth of ~310 m, permeability values of 0.05 mD and 860 mD have been reported in two different wells (see Figure 16).The major reasons behind such differences include geological structures, coal maturity, permeability measurement methods, stress regime, fracture/cleat density, localized shear zones and cleat mineralization [81,117,[180][181][182][183].In situ permeability measurement in a CSG well provides an average value for permeability for all intersected coal seams in the tested interval.Since up to 45 individual thin coal seams can be intersected in the Walloon Subgroup [177], it is difficult to determine vertical heterogeneity of the system.

Permeability (mD) Taroom CM Upper Juandah Lower Juandah
The vertical heterogeneity and lateral discontinuity of the coal seams in the Surat Basin pose a great challenge on production history matching to obtain reservoir properties [184].Mazumder, Jiang [178] performed production history matching on 14 wells in the Surat Basin and reported permeability values in the range of 65-310 mD for the wells completed in the Taroom Coals, 90-275 mD for the wells completed in the Juandah Coals and 75-320 mD for the wells completed commingled across both Taroom and Juandah coals.The commingled production in the Surat Basin can result in coal seams with lower reservoir pressures to contribute to gas production later when high-pressure seams are depleted.This may result in increase in gas production rate after a period of decline [178].

Cooper Basin
Thick and laterally extensive Permian deep coal seams (>2 km) in the Cooper Basin, South Australia, are considered huge, underdeveloped resources.Patchawarra Formation accommodates a large bulk of deep coal seams and is the primary target for development.
Figure 17 shows the isolith map of the Patchawarra Coal showing net coal thickness across the Cooper Basin [7].The thickness of an individual coal seam in the Patchawarra formation can reach to 45 m (148 ft.) [5,29].
A number of reported gas contents for the upper Patchawarra Formation, including those tested in Dorodillo-2, are very high.This is likely an overestimation because they were conventionally cored and core retrieval process results in extended lost gas time.Hence, lost gas volume estimates are higher than actual resulting in overestimation of gas contents.Diamond and Schatzel [133] suggested that lost gas volume can be reasonably estimated when lost gas time is less than one hour.In Dorodillo-2, lost gas time was 6 h and 19 min [185] that is much longer than the limit of 1 h suggested by Diamond and Schatzel [133].Gas contents of core samples in deep coal wells that are conventionally cored should be treated very carefully because extended lost gas time increases the error in estimation method.More recent mixed isotherm analysis of three core samples of varying coal rank, suggests the total measured gas is likely within the range of 11.04-19.26m 3 /ton (daf) (390-680 scf/ton) for Patchawarra deep coals (Bindah-3; Kirralee-2; and Bobs Well-2).Total gas content in deep coal seams can be measured more accurately using pressure coring techniques.In the southern part of the Cooper Basin, Patchawarra coal seams are thick, shallower, and higher in permeability but lower in gas content.Vu upper coal seam interested in Le Chiffre-1 and Klebb-1 is 35 m thick.Figure 18 shows gas contents in Le Chiffre-1 that has 461 m of cores that were taken from 1464 m to 2044 m using the wireline retrievable coring [189].The lost gas time for samples is generally less than 1 h and gas contents range from 1.36-6.81m 3 /ton (daf).The side-wall pressure cores obtained later from Klebb-5 (located in proximity to Le Chiffre-1 and Klebb-1) indicated that average gas content (methane) in the Klebb area is 5.5 m 3 /ton (194 scf/ton) [190].Vm3, Vu upper and lower coal seams in the Klebb pilot have low gas saturation.Hence, a significant depressurization phase is needed to achieve commercial production rates.This is a barrier for development of deep coal seams in the southern part of the Cooper Basin.Permeability of deep coal seams is in the range of micro-to nano-darcy [29] hence, hydraulic fracturing is necessary to achieve economical gas rate.The Cooper Basin deep coal "add-on" frac campaign involves placing hydraulic fractures in thick deep coal seams in gas development wells.This program is considered an overall economic success with variable productivity [3].These development wells mostly target tight sandstone reservoirs in the Cooper Basin and the excess gas from deep coal seams adds further contribution to overall production.Average gas production from single-stage fracture stimulated coals is 0.2 mmscf/day and that can reach to 1 mmscf/day in some wells [3].
Dunlop, Salmachi [5], utilized a comprehensive data set from a hydraulically fractured well, completed in deep Patchawarra coal seams of the Cooper Basin, and demonstrated that hydraulic fracture conductivity can improve throughout the life of the well.The increase in fracture conductivity could be attributed to the cumulative impact Le Chiffre-1 of a number of factors including coal matrix shrinkage, recovery of residual hydraulic fracture fluid and erosion of the fracture faces due to proppant embedment.An earlier study on this wellbore by Salmachi, Dunlop [7] concluded that depletion is limited to the stimulated reservoir area adjacent to hydraulic fractures.The excessive gas desorption from near wellbore area can potentially promote matrix shrinkage effect.The matrix shrinkage can enlarge hydraulic fracture aperture and improve fracture conductivity.Further research is needed to understand the mechanisms contributing to increase in hydraulic fracture conductivity in deep coal wells.

Summary
In late 2018, Australia became the largest exporter of LNG in the world and 29% of the country's LNG nameplate capacity is in three facilities in east coast that are primarily fed by coal seam gas.The current estimate for the Australia's CSG reserves is 43 Tcf from which 92.5% lie within Surat and Bowen basins of Queensland.The 13% Gas Scheme and the idea of converting CSG to LNG for exports were two critical decisions that boosted commercial CSG production in Queensland.In the Bowen Basin, the three stratigraphic units that are major targets for CSG production include Moranbah Coal Measures/German Creek Formation, the Baralaba Coal Measures and the Bandanna Formation.In the Surat Basin, the major producing CSG interval is the Jurassic Walloon Coal Measures.Coals in the Surat Basin are low in rank and gas content is generally lower than 10 m 3 /ton (daf).The key to success in CSG development in the Surat Basin is high permeability that can range from 5-5000 mD.In the Bowen Basin, coals are higher in rank (low volatile bituminous and anthracite) and gas contents are generally higher than those in the Surat Basin are.The significant variability in gas content and permeability in Australian CSG basins are attributed to geological structures, coal rank, stress regime, measurement methods, magmatism and sealing capacity of overlying formations.Field studies indicate that coal permeability in Australian fields can vary during depletion and the functional form of permeability change with reservoir pressure is exponential, consistent with producing CSG fields in North America.
Detailed analyses of high-resolution borehole image logs in CSG wells indicate that there are major perturbations in SHmax orientation, both spatially and with depth in Australian CSG basins.The regional orientation of SHmax in the northern part of the Bowen Basin is NNE-SSW.In the Clarence-Moreton, Surat and southern part of the Bowen basins, the regional SHmax orientation is NE-SW to ENE-WSW.The Cooper Basin in central Australia has a significant E-W SHmax orientation.In the Gunnedah Basin, the geological structures, including faults, fracture and magmatism have strongly perturbed stresses in the region.Analysis of micro-seismic and tilt-meter results indicate that in areas with high differential stresses (in reverse stress regime) planar hydraulic fractures propagate.In areas with low differential stresses, both stresses and pre-exiting natural fractures control hydraulic fracture propagation resulting in more complex fracture patterns.
In 2007, for the first time, Patchawarra coal seams were hydraulically fractured in Moomba 77 and a gas flowrate of ~0.1 mmscf/day was reported.Extended lost gas times in deep coal wells that are cored conventionally pose a challenge to accurately measure gas contents.Pressure core can provide a more accurate estimation of gas content in deep coal seams than conventional desorption canisters do.Since 2013, gas production exclusively from deep coal seams has been tested as a single add-on fracture stimulation in vertical well completions across the Cooper Basin.
Informed Consent Statement: Not applicable.
Data Availability Statement: Majority of data presented in this study are openly available in government websites.Please refer to the reference section to find the address for each dataset.Restrictions may apply to the availability of some data presented in this study.

Figure 1 .
Figure 1.Location of major coal seam gas basins in eastern Australia.Wellbores (black dots) and associated infrastructure including oil and gas pipelines are also shown.

 Phase 1 :
Figure 2 compares CSG production in Australia with CSG production in the US.It shows that the US production peaks in 2008 at 1966 BCF while during the same year, Australia production is only 150 BCF.The US production declines afterward and in 2016, Australia production surpasses the US production.The annual CSG production (see Figure 2) and number of wells drilled each year (see Figure 3) demonstrate six key phases in history of CSG in Australia.Although some CSG drilling has occurred in nearly all Australian mainland states, it is dominated by Queensland and New South Wales.The six phases were dominantly controlled by commerciality and overcoming technical difficulties while at the same time, working closely with local communities. Phase 1: Pre-1975  Phase 2: 1975-1996 (Early Exploration)  Phase 3: 1996-2005 (Early Production)  Phase 4: 2005-2011 (Political Intervention)  Phase 5: 2011-2016 (CSG to LNG)  Phase 6: Post-2016

Figure 3 .
Figure 3. Number of CSG wells drilled per year in Australia since the 1980s [15].

Figure 5 .
Figure 5. Dedicated deep coal wells drilled/completed in the Cooper Basin, South Australia.Moomba 77 is the first dedicated well completed in the Patchawarra coal in 2007 by the hydraulic fracturing.

Figure 6 .
Figure 6.Stress provinces across the Australian continent [14].The orange lines show the smoothed pattern of maximum horizontal stress (SHmax) that demonstrate the variability of SHmax across the Australian continent.

Figure 7 .
Figure 7.In situ stress pattern of major coal seam gas basins in eastern and central Australia and the smoothed stress pattern (orange lines) showing the regional orientation of maximum horizontal stress[14].Different symbols show the method of measurement for stress data including focal mechanism solution of earthquakes, borehole breakouts, drilling induced fractures, over coring, and hydraulic fracturing measurements (see[14] for the details of each method).Different colours indicate different tectonic stress regime (NF, normal; SS, strike-slip; TF, thrust; U, unknown).The length of the lines for each stress indicator illustrates the quality of stress orientation based on the World Stress Map ranking criteria[124].

Figure 9 .
Figure 9.Gas content versus depth for the Bandanna Coals of the Bowen Basin.Variability in gas content is observed across the field (data sourced from Queensland Government[154]).

Figure 10 .
Figure 10.Change in permeability as a function of average reservoir pressure above and below the desorption pressure [160].There is a slight decrease in permeability above desorption pressure due to stress-dependency of permeability followed by a sharp increase in a short period due to matrix shrinkage effect and then permeability remains unchanged (This figure was originally published in Journal of Natural Gas Science and Engineering Volume 40 Pages 51-60).

Figure 11 .
Figure 11.Change of effective permeability as a function of reservoir pressure in the Fairview field, Bowen Basin [155] (This figure was originally published in The APPEA Journal Vol 59 No 1).

Figure 12 .
Figure 12.Integration of rate transient analysis and pressure transient analysis to obtain the functional form of permeability change [155] (This figure was originally published in The APPEA Journal Vol 59 No 1).(a) Comparison of field data with tank type model when permeability is constant, (b) bottom-hole flowing pressures and average reservoir pressures, (c) permeability function obtained using RTA, (d) the semi-log plot of permeability ratio as a function of average reservoir pressure is a straight-line (exponential increase).

Figure 13 .
Figure 13.Gas content versus depth for the Hosskisons coal in the Gunnedah Basin.Significant variability in gas content is observed [81] (This figure was originally published in The International Journal of Coal Geology Volume 165 Pages 278-289).

Figure 14 .
Figure 14.Cleat and fracture systems characterization by integrated borehole image log interpretation and DST data analysis for the Hoskissons Coals in six wells across the Gunnedah Basin [81] (This figure was originally published in The International Journal of Coal Geology Volume 165 Pages 278-289).

Figure 15 .
Figure 15.Gas content versus depth in the Surat Basin based on gas desorption data collected from 86 wells [154].

Figure 17 .
Figure 17.(a) Isolith map of the Patchawarra Coal showing net coal thickness in the Patchawarra Formation across the Cooper Basin, (b) the stratigraphic column of the Cooper Basin [186-188] (This figure was originally published in The AAPG Bulletin Vol 103 Issue 1 Pages 91-107).

Figure 18 .
Figure 18.Gas content of deep coal seams in Le Chiffre-1 drilled in the Weena Trough, southern part of the Cooper Basin [189].

Table 1 .
CSG to LNG projects in Australia.

Table 5 .
Some permeability enhancement studies in the Bowen basin(Australia)and the San Juan Basin in the United States.