Creep Rupture and Permeability Evolution in High Temperature Heat-Treated Sandstone Containing Pre-Existing Twin Flaws

: Utilizing underground coal gasiﬁcation cavities for carbon capture and sequestration provides a potentially economic and sustainable solution to a vexing environmental and energy problem. The thermal inﬂuence on creep properties and long-term permeability evolution around the underground gasiﬁcation chamber is a key issue in UCG-CCS operation in containing fugitive emissions. We complete multi-step loading and unloading creep tests with permeability measurement at conﬁning stresses of 30 MPa on pre-cracked sandstone specimens thermally heat-treated to 250, 500, 750 and 1000 ◦ C. Observations indicate a critical threshold temperature of 500 ◦ C required to initiate thermally-induced cracks with subsequent strength reduction occurring at 750 ◦ C. Comparison of histories of creep, visco-elastic and visco-plastic strains highlight the existence of a strain jump at a certain deviatoric stress level—where the intervening rock bridge between the twin starter-cracks is eliminated. As the deviatoric stress level increases, the visco-plastic strains make up an important composition of total creep strain, especially for specimens pre-treated at higher temperatures, and the development of the visco-plastic strain leads to the time-dependent failure of the rock. The thermal pre-treatment produces thermal cracks with their closure resulting in increased instantaneous elastic strains and instantaneous plastic strains. With increasing stress ratio, the steady-state creep rates increase slowly before the failure stress ratio but rise suddenly over the ﬁnal stress ratio to failure. However, the pre-treatment temperature has no clear and apparent inﬂuence on steady creep strain rates. Rock specimens subject to higher pre-treatment temperatures exhibit higher permeabilities. The pre-existing cracks close under compression with a coplanar shear crack propagating from the starter-cracks and ultimately linking these formerly separate cracks. In addition, it is clear that the specimens pre-treated at higher temperatures accommodate greater damage.


Introduction
Underground coal gasification (UCG) is a technology that converts hydrocarbons into syngas in-situ [1,2].Crucially, UCG combined with carbon capture and sequestration (CCS) may have a special promise [3].While UCG-CCS has potential promise, some key environmental and engineering issues require resolution.We focus on key rock mechanics problems including excessive subsidence and fugitive gases-each critically impacted by the response of strength and permeability to the imposed high temperatures.Similar to long-wall mining, excavation of coal from depth may cause subsidence, especially over long-term operations and under high crustal stress.It ought to be possible to inject CO 2 into former UCG voids which are currently at a depth of more than 800 m [4].The UCG process is intensely exothermic with the temperatures occasionally exceeding 900 • C in the combustion zone.Prior observations indicate that high temperatures significantly transform the physical and chemical properties and resulting mechanical behavior of rocks (e.g., heating and quenching the first caprock may result in cracking or other behavior that compromises the integrity of the storage [5]).In addition, the high permeability of the UCG goaf for geological CCS projects, CO 2 storage in artificially developed UCG goaf with high permeability appears to be an attractive proposition [4], provided storage security of the caprocks can be guaranteed.
Heat transfer mechanisms in rocks are complex and nonlinear [6][7][8][9][10][11]. UCG production boreholes typically experience extreme temperatures of the order of ~900 • C which significantly exceeds the critical temperature (Tc~400-500 • C) that induces changes in the mechanical and permeability behavior of typical sandstones [11].Temperatures from ambient to 400 • C correspond to vaporization then escape of adhered water, combined water and structural water.Between 400-600 • C, the minerals in the sandstone thermally react and transform, resulting in porosity increase, reduction in permeability and diffusivity, and changes in heat capacity [12].Reservoir rheology is also sensitive to temperature [6,13,14].An increase in temperature accelerates creep strains and consequently reduces the time-tofailure [15].The total axial strain increases with increasing temperature, with instantaneous elastic strain and instantaneous plastic strain increasing slightly as the temperature increases from 25 to 700 • C but then increasing substantially as the temperature reaches 1000 • C [16].However, the response of fractured rocks is usually different from that of intact rock-since stress concentrations are present at the ends of pre-existing cracks-promoting the initiation of new cracks and the reduction in rock mass strength [17,18].Previous investigations [19][20][21][22] addressed crack initiation and coalescence for jointed rock masses under short-term loading with key mechanisms of creep contributing to time-dependent deformation and crack propagation [23,24].Investigation of time-dependent response of jointed rock [21,[25][26][27][28] have defined key failure mechanisms but few characterizations study the comprehensive impact of thermal effects on time-dependent interaction of thermal cracking across the scales.
The high-temperature treatment of the rock mass during underground coal gasification leads to the formation of fractures and changes in the fluid transmission characteristics of the strata, which is a direct result of the combustion process [29].Thermal cracks result and increase permeability as the temperature increases from 400-600 • C [10].Prior observations have promoted our understanding of permeability evolution due to shortterm heating and loading demonstrating a linear correlation with volumetric strain [24,30].However, the timescales of permeability evolution of the goaf in UCG-CCS operations should also be noted.During creep loading the permeability will typically decrease before increasing with the onset of failure at late deformation stages [28,30].Some observations reveal that the temperature effects are on timescales of permeability evolution and that during the multi-step loading and unloading cycles process the permeability first decreases slightly as temperature increases from 25 to 300 • C and then increases with increasing temperature [28].Heating rate significantly affects thermal damage as evident in that permeability is exponentially correlated with the number of AE events [24].
Consequences of the combined presence of pre-existing fissures, high confining pressures and high temperatures on the evolution of strength and permeability for UCG-CCS operations has not been previously considered.Correspondingly, the following focuses on defining the influence of thermal damage on creep properties and permeability evolution of pre-cracked specimens of red sandstone from micro-scale to laboratory scale.The primary focus is to identify creep parameters and overall permeability evolution in the pre-cracked rock in response to different pre-treatment temperatures in defining characteristics of the evolution of transport properties and mechanical strength.

Sandstone Specimens and Heating Procedure
The red sandstone used for these experiments was collected from Rizhao City in Shandong Province, China.The dry density is 2402 kg/m 3 , and effective porosity is ~6.26-6.55%[28].The sandstone mainly comprises quartz, feldspar, dolomite, hematite and clay minerals.The tested specimens were cylindrical with 50 mm in diameter and 100 mm in height.Then, two fissures of near-identical length (2a = 12 mm) and width (d = 2 mm) were cut in the specimens by diamond wire saw.The inclination angle of two coplanar cracks was 45 • with an intervening rock bridge of 2b = 12 mm.The samples were dried for ~24 h and then heated to 250, 500, 750 and 1000 • C at a rate of 5 • C/min.The specimens were held at their target temperature for two hours before cooling naturally to room temperature (Figure 1).The color of the specimens changed as a result of the heat-treatment changing from drab-red (250 • C), progressively redder (500 and 750 • C) before becoming brown at 1000 • C. Figure 2 shows SEM images of the red sandstone after different pre-treatment temperatures.At 250 • C (Figure 2a), minor micro-cracks appear, and these develop into boundary cracks at the edges of grains during thermal cooling and differential expansion/contraction (Figure 2b,c).In addition, some transgranular cracks appear across grains.When the temperature reaches 1000 • C (Figure 2d) the density of boundary cracks increases and the crack apertures between grains also increases.The number of transgranular cracks also increases significantly.The results indicate a threshold temperature to produce thermally induced cracks in this study is 500 • C.

Testing Procedure
Conventional triaxial compression and multi-step loading-and unloading-cycle creep tests were completed under a confining pressure of 30 MPa (Figure 3a).The test system consists of three servo-hydraulic control units applying: axial, confining pressure and pore pressure loads.The maximum axial loading capacity of testing system is 800 kN, and the maximum confining pressure and pore pressure can reach 60 MPa.The axial deformation of the specimen was measured with a pair of linear variable displacement transducers (LVDTs), and the circumferential deformation was measured with a circumferential extensometer (Figure 3b).The conventional triaxial compression tests were first completed to obtain the peak strength of the double pre-cracked sandstone specimens after exposure different temperature treatments.The test procedure is as follows: (i) at a loading rate of 4 MPa/min, 30 MPa hydrostatic pressure was applied to the specimen; (ii) then the specimen was axial loaded at a rate of 0.04 mm/min until it reached the residual strength.The stress-strain curves and triaxial compressive strengths (TCS) were thereby obtained.
Multi-step loading-and unloading-cycle creep tests were subsequently completed to study the creep characteristics of the pre-cracked sandstone.Stress levels for the creep tests were set according to the TCS.The test procedure was as follows: (i) the specimens were subjected to 30 MPa hydrostatic pressure at a rate of 4 MPa/min; (ii) a pre-determined first stress level was applied under pressure-controlled condition of 5 MPa/min; (iii) the specimen was then allowed to deform for approximately 48 h under constant stress; (iv) the axial deviatoric stress was then unloaded to 0 MPa under pressure-controlled conditions of 5 MPa/min; (v) these unloading conditions were retained for 24 h.Next, these loading and unloading cycles to creep were repeated until creep failure occurred.During the test, the axial and lateral strains of the pre-cracked red sandstone specimens were recorded continuously.The basic physical parameters of the specimens are listed in Table 1 with each of the deviatoric stress levels shown in Table 2.
With inert nitrogen as the flow medium, the permeability of pre-cracked specimens was measured by steady-state flow method during the creep test.The gas pressures at upstream and downstream were 3 MPa and 0.1 MPa, respectively.The permeability was calculated from Equation (1).
where Q represents volumetric gas flow rate (m 3 /s), µ represents the viscosity of the nitrogen (Pa•s); L and A represent the length (m) and the cross-sectional area (m 2 ) of the rock specimen, respectively; P u and P d represent the inlet and outlet gas pressures (MPa), respectively.We observe that as the temperature increases the elastic modulus decreases and the Poisson's ratio overall increases (Table 1).

Creep Test Results
The axial deviatoric stress-strain curves of the pre-cracked sandstone specimens under conventional triaxial compression and cyclic creep compression are presented in Figure 4.The peak strength obtained from the cyclic creep compression is lower than that obtained from conventional compression for the different pre-treated temperatures.The creep curves plateau at each stress level, representing creep deformations under each constant stress.Each loading produced plastic deformation, with the plastic deformation increasing with an increase in the deviatoric stress.In addition, as with the increases in stress levels, the length of the stress-plateau also gradually increases, except a certain threshold stress level where coalescence of two coplanar cracks produced a larger plastic deformation, such as Figure 4a-c.A stress drop occurs after reaching the first peak strength due to the closure of the pre-existing cracks.This is because the ensemble closure of the pre-existing cracks and newly formed fractures increases the load capacity of the specimen.The resulting second peak strength was always larger than the first peak strength due to this closure of the pre-existing cracks.It should be noted that this double peak strength phenomenon is different from that for intact specimens that show only a single peak strength because intact specimens without pre-existing cracks the plastic hardening is not obvious [28,30].Figure 5 illustrates the variations in the axial and lateral strains over time for the pre-cracked specimens for different pre-treatment temperatures subject to multi-level creep loading and unloading.Instantaneous axial and lateral strains are produced as the deviatoric stress is applied with the elastic strains recovered as the deviatoric stress is removed.The increasing number of loading and unloading cycles will lead to more unrecoverable axial deformation, especially at a certain level of deviatoric stress that causes the closure of pre-existing cracks and produces large irrecoverable deformations, such as the bule dotted oval in Figure 5a-c.When the target deviatoric stress was reached, the deviatoric stress remained constant for 48 h.When the stress level is relatively low, deceleration deformation and steady creep deformation are the main components of axial deformation.As the deviatoric stress increases, the duration of the primary creep tends to increase, and when the deviatoric stress exceeds a threshold stress level, a fracture propagated between the two colinear cracks, producing a large plastic deformation.For instance, the specimens pre-treated to 250 and 500 • C (Figure 5a,b) fractured during the 3rd stress level while the specimen exposed to 750 • C (Figure 5c) fractured during only the 2nd stress level.Whereas primary, secondary and tertiary creep occurred under creep failure stress levels, as shown in the red dotted oval.The unloading deviatoric stress was 0 MPa since the elastic hysteretic behavior of rock maintained unloading conditions for 24 h.

Strain Separation
During creep experiments, the specific strains are separated by cyclic loading and unloading, including components of both recoverable strain and irrecoverable strain.The recoverable strain consists of the instantaneous elastic strain (ε me ) and visco-elastic strain (ε ve ), whereas the irrecoverable strain is composed of the instantaneous plastic strain (ε mp ) and visco-plastic strain (ε vp ), as shown in Figure 6.Moreover, total strain (ε) can also be divided into instantaneous (ε m ) and creep (ε c ) strains, which can be expressed as The instantaneous strain ε m is the measured strain when the targeted deviatoric stress level is reached and the creep strain ε c refers to the increase of strain measured with time at the stage of constant loading.The instantaneous strain ε m is composed of two parts, i.e., the instantaneous elastic strain ε me and the instantaneous plastic strain ε mp : The creep strain ε c also consists of two parts, the recoverable visco-elastic strain ε ve and the unrecoverable visco-plastic strain ε vp : The duration of the creep loading is often longer than the duration of unloadingfor instance, in this present study, the duration of the loading stage is 48 h, while the duration of unloading is 24 h.Hence, no visco-elastic recovery data are available from unloading for only 24 h.To obtain the visco-plastic strains under different stress levels, the visco-elastic strains based on experimental results were fitted by Equation (5).Then, the difference between the creep strain and the fitted visco-elastic strain may be evaluated as the visco-plastic strain under the loading conditions [16].
where A, B and m are fitting parameters (m > 0).Due to accelerated creep failure, the visco-elastic strain and visco-plastic strain in the third creep cannot be separated.and the tertiary creep should be represented by coupling visco-plasticity to damage [31][32][33].In this study, the deformation after failure of the specimen is considered irreversible and the entire creep strain is taken as the visco-plastic strain.The creep strain histories of the thermally-treated pre-cracked specimens are separated into visco-elastic and visco-plastic strains according to the prior separation method.The variation of the total creep strains, visco-elastic strains and visco-plastic strains with time for four specimens are plotted in Figure 7 (the subfigure in dashed box is the enlarged version of the dashed part in main figure).For the specimen pre-treated at 250 • C (Figure 7a), there are 4 levels of deviatoric stress in total.The creep and visco-plastic strains can be characterized by primary creep and steady-state creep during the 1st and 2nd deviatoric stress levels.At the 3rd deviatoric stress level, as a result of crack growth, the creep strain and visco-plastic strain rates gradually increases until 95.597 h when the rock bridge between the two cracks fractures.At that point, the strain rate decreases, and begins steadystate creep at 95.61 h.At the 4th deviatoric stress level, because the ultimate failure of the specimen, the creep strain cannot be separated-therefore, it is considered that the creep strain is equal to the visco-plastic strain.The creep strain sequence can be characterized into primary, steady-state and tertiary creep.The visco-elastic strain can be characterized as primary and steady-state creep over the 1st stress interval and as primary creep from the 2nd deviatoric stress level.It should be noted that, the creep strain curves of the 3rd and 4th deviatoric stress levels are different: the creep rate under the 3rd deviatoric stress level first increases, and then decreases; however, the creep rate under 4th deviatoric stress level is opposite and is characterized by a typical creep curve.Figure 12 plots D against the different stress ratios (SRs).As SR increases, D also increases-the higher the stress level the greater the damage.This is especially true when the SR increases from 0.7 to 0.8 for the specimens pre-treated at 250 and 500 °C, and as SR increases from 0.6 to 0.7 for the specimen pre-treated at 750 °C where D increases sharply due to the closure of pre-existing cracks.For identical SRs, the specimen with the higher temperature pre-treatment usually accumulates the greater damage.For example, under a stress ratio of 0.6, the D of the specimen subjected to 1000 °C is 0.56 and is larger than The separated strains for the specimen pre-heated to 500 • C are shown in Figure 7b.The results are similar to that for the 250 • C heated specimen as the creep and visco-plastic strains show both primary and steady-state creep at the 1st to 2nd deviatoric stress levels with a strain jump at the 3rd deviatoric stress level.At the 4th deviatoric stress level, the creep and visco-plastic strains are larger than those at the 1st and 2nd deviatoric stress level but smaller than at the 3rd deviatoric stress level that exhibits primary and steady-state creep.As the deviatoric stress is increased, the primary creep, steady-state creep and tertiary creep clearly exhibit at the 5th deviatoric stress level when the specimen finally fails during tertiary creep.
Figure 7c shows the separated strains for the 750 • C treatment.The creep strain curve and the visco-plastic strain curve are characterized by primary and steady-state creep at the 1st deviatoric stress levels but with very small strain rates of steady-state creep.The visco-elastic strains are lower than the visco-plastic strains at the 1st deviatoric stress level.At the 2nd deviatoric stress level, in the initial loading stage, the creep and visco-plastic strain increases rapidly, but the rate of increase decreases until 47.5 h when the creep and visco-plastic strain jump from 0.85 × 10 −3 , 0.83 × 10 −3 to 10.37 × 10 −3 , 10.33 × 10 −3 , respectively, but the visco-elastic strain only increases from 0.027 × 10 −3 to 0.042 × 10 −3 , before establishing a steady-state.The specimen ultimately fails at the 4th deviatoric stress level after 0.11 h.
The separated strains for the 1000 • C specimen are different from the previously noted specimen behaviors (Figure 7d).First, the creep and visco-plastic strains are larger than those for specimens pre-heated to 250, 500 and 750 • C. Secondly, as the creep progresses, the creep and visco-plastic strains increase smoothly and progressively with no strain jump occurring until final failure.This results from the damage produced by the intense thermal cracking-during creep loading, closure of these cracks consumes most of the strain energy with insufficient remainder to fracture the rock.
Comparison of histories of creep, visco-elastic and visco-plastic strains highlight the existence of a strain jump at a certain deviatoric stress level-where the intervening rock bridge between the two starter-cracks fractures.However, the joining of these preexisting cracks does not result in complete failure, but merely produces a larger irreversible visco-plastic strain with reduced ultimate strength.As the deviatoric stress level further increases, the visco-plastic strain accounts for an important proportion of total creep strain, especially for the high pre-treatment temperature specimens, with the development of the visco-plastic strain leading to the time-dependent failure of the rock.

Pre-Treatment Temperature Effects
These specimens exhibit dual peak strengths in compression as a result of the twin coplanar flaws (Figure 4), as illustrated by the corresponding axial strain versus pretreatment temperature in Figure 8. Figure 8a shows the two peak strengths obtained from both the triaxial and creep tests both increase as the temperature increases from 250-750 • C, but then decreases substantially as temperatures reach 1000 • C. The principal reason for this is the vaporization and escape of adhered water, combined water and structural water between temperatures of 250-750 • C, which leads to an increase in the coefficient of internal friction [11] and increase in strength.When the temperature reaches 1000 • C, boundary and transgranular cracks develop and intergranular cracks widen-again resulting in a reduction in strength.In addition to its influence on ultimate strength, thermal treatment also has a influence on the deformation characteristics.Figure 8b shows changes in the peak strains with pre-treatment temperature.The axial strain at both peak strengths increase with pre-treatment temperature as it increases from 250 to 1000 • C. To better characterize the effects of temperature on different separated strains, those strains are plotted against stress ratio for different pre-treatment temperatures Figures 9  and 10. Figure 9 shows the variation in the instantaneous elastic strains and the instantaneous plastic strains with stress ratio for the four pre-treatment temperatures.It is apparent from Figure 9a that the instantaneous elastic strains of the pre-cracked specimens increase linearly with an increase in the stress ratio.The instantaneous elastic strains increase with pre-treatment temperature, for example, the difference in value a between 250 • C and 500 • C is less than the difference in value b between 500 • C and 750 • C, and both are significantly less than the difference in value c between 750 • C and 1000 • C. Figure 9b illustrates the variation of instantaneous plastic strains with stress ratio showing that between 250~750 • C there is no obvious change with the initial stress ratio.However, when reaching the stress level for fracture, the instantaneous strain increases significantly due to the closure of pre-existing cracks.The instantaneous plastic strain of the specimen under a treatment temperature of 1000 • C is relatively large due to the closure of thermally-induced cracks with the instantaneous plastic strain increasing with an increase in the stress ratio.From the above analysis, the instantaneous elastic deformation is usually larger than the instantaneous plastic deformation, but the closure of pre-existing cracks and the fracture of rock bridges produces a large plastic deformation that presents unconventional response.The instantaneous plastic strain exceeds any instantaneous elastic strain when closure and fracture occur.The thermal pre-treatment produces thermal cracks with their closure resulting in increased instantaneous elastic strains and instantaneous plastic strains.
Figure 10 shows the evolution of the visco-elastic and visco-plastic strains with stress ratio.The visco-elastic strains increase nonlinearly with an increase in the stress ratio and exhibit a complex relationship with temperature.The visco-elastic strains of the precracked specimen after pre-treatment at 250 and 500 • C are almost equal, but smaller than that of the specimen pre-treated at 750 • C.However, as the pre-treatment temperature increases to 1000 • C, the visco-elastic strain is less than that for the specimen pre-treated at 750 • C because a larger irreversible plastic deformation is produced by the closure of thermally-induced cracks of the specimen pre-treated at 1000 • C. The visco-plastic strain apparent in the 1000 • C specimen (Figure 10b) also supports this-the visco-plastic strain of this hot specimen is larger than that for lower temperatures, except at certain stress ratios where the closure of pre-cracks and fracturing of the rock bridge causes a significant plastic deformation.
The creep strain rates during steady creep can be obtained by taking the derivative of steady creep curves, considered as of constant gradient.Figure 11 shows the steady-state creep strain rate versus stress ratio for all tested specimens.The results show that the steady-state creep rates increase slowly with stress ratio before the failure stress ratio but rise suddenly over the final stress ratio.The steady-state creep strain rate over the last stage is usually several orders of magnitude larger than that of the previous stages-requiring the use of logarithmic rates.The pre-heated temperature has no clear nor apparent influence on steady creep strain rates.The degree of material deterioration under changing stress states can be quantitatively described by the damage variable (D) [34].Usually, D can be calculated using the elastic modulus, ultrasonic wave velocity, density and indexed to energy adsorption, strain and acoustic emission [35].However, in cyclic loading-unloading creep tests, instantaneous elastic and visco-elastic strains are recovered after unloading, but instantaneous plastic and visco-plastic strains are not.Hence, the damage variable may be calculated using the ratio of non-elastic strain (instant plastic strain and visco-plastic strain) to total strains.Finally, D can be calculated from the instantaneous elastic and visco-elastic strains as [28]: where ε e and ε are the elastic and the total strains, respectively; ε me and ε cve represent the instantaneous elastic and the visco-elastic strains, respectively; ε mp and ε cvp represent the instantaneous plastic and the visco-plastic strains, respectively; and ε m and ε c represent the instantaneous and the creep strains, respectively.
Figure 12 plots D against the different stress ratios (SRs).As SR increases, D also increases-the higher the stress level the greater the damage.This is especially true when the SR increases from 0.7 to 0.8 for the specimens pre-treated at 250 and 500 • C, and as SR increases from 0.6 to 0.7 for the specimen pre-treated at 750 • C where D increases sharply due to the closure of pre-existing cracks.For identical SRs, the specimen with the higher temperature pre-treatment usually accumulates the greater damage.For example, under a stress ratio of 0.6, the D of the specimen subjected to 1000 • C is 0.56 and is larger than the 0.31 for the 750 • C, and much larger than the 0.11 and 0.13 for the 250 and 500 • C specimens.

Permeability Evolution
The evolution of volumetric creep and gas permeability for pre-treatment temperatures of 250 and 500 • C is presented in Figure 13.Equation ( 1) is adopted to calculate the coefficient of permeability.Permeability is closely related to volumetric strain during the creep process.During cyclic loading and unloading creep tests (Figure 13a) the permeability of sandstone remains near constant at low deviatoric stress (first two stress) before slightly increasing with unloading.The volumetric strain during the 3rd deviatoric stress level sharply increases and then decreases, resulting in the permeability also first increasing before subsequently decreasing.This results from the closure of pre-existing cracks and the production of some new cracks, although the newly produced cracks subsequently close under compression.At the final deviatoric stress level, the volume of the sandstone specimen initially decreases but then dilates before failure-the permeability shows an inverse trend to that expected from the volume strain signal.When the specimen fails, macro-cracks are generated and the permeability increases.Figure 13b shows the relationship between volumetric strain and permeability for the specimen after pre-heating to 500 • C. The permeability of the pre-cracked sandstone decreases with loading and increases with unloading of the deviatoric stress (Figure 13b).This results from the re-opening of pre-existing pores and fissures.With an increase in the deviatoric stress from 293 MPa to 247 MPa, the permeability due to each loading and unloading creep cycle becomes progressively lower-this can be explained by the presumed material strengthening of the sandstone under cyclic loading [30].Over the final deviatoric stress level, the large macroscopic cracks linking the starter cracks have already been formed and this precipitates structural failure of the specimen-the permeability then increases significantly to 45 × 10 −18 m 2 .Rock specimens undergoing hotter pre-treatments exhibit higher permeabilities.The average permeability of specimens pre-heated to 250 • C before failure is 5 × 10 −18 m 2 which is almost an order of magnitude lower than that for 500 • C (45 × 10 −18 m 2 ).The higher temperature pre-treatment produces more thermal cracks and provides more connected channels for gas transport.

Failure Modes
Failure modes are presented in Figure 14 showing the presence of a macro shear crack as indicative of failure with some tensile cracks induced by the shear sliding.The pre-existing cracks are closed under compression with a coplanar shear crack linking these formerly separate cracks.The pre-existing cracks exert a considerable influence on the creep failure mode in each specimen.In addition, it is clear that the specimens pre-treated at higher temperatures accommodate greater damage.For example, the shear fracture for specimens heated to 750 and 1000 • C is clearly wider than that for specimens pre-treated to 250 and 500 • C. In addition, the specimens undergoing higher temperature pre-treatments (Figure 14c,d) accumulated more severe damage than those at lower temperatures (Figure 14a,b)-apparent in the loss of some spalling along the main shear crack (Figure 14c,d).
We choose five stages to describe crack evolution in the loading and unloading creep tests (Figure 15).Due to the interaction of thermal cracks (produced by high temperature) and pre-existing cracks, the crack evolution is more complex than that in isothermal untreated intact specimens.Stage a represents the initial state with the specimen containing two coplanar pre-existing cracks dipping at 45 • relative to the direction of the maximum principal stress.We conclude from Figure 2 that the specimen contains some thermal cracks and micropores.Stage b corresponds to the first level of loading.Under compression, the micropores and some thermal cracks close and produce a small number of cracks, but the stress is insufficient to cause wholesale fracture.Stage c is when the second unloading is completed.As the loading stress is removed, the closed micropores and cracks partially recover, elastically.Stage d produces a collapse as a stress concentration at the ends of the pre-existing cracks causes a new shear crack to connect the two pre-existing cracks.A large axial deformation results due to the closure of the pre-existing cracks.Stage e is the ultimate failure of the specimen.A macro-crack forms colinear to the pre-existing cracks, linking them, and which provides a connected channel for gas flow that yields a corresponding jump in permeability.

Figure 1 .
Figure 1.Pre-cracked red sandstone specimens after exposure to different pre-treatment temperature.

Figure 11 .
Figure 11.Relationships between stress ratio with steady-state strain ratio and damage after different temperature treatment.

Figure 8 .
Figure 8.The temperature effects on strength and axial strain.(a) Strength; (b) axial strain.

Figure 9 .
Figure 9. Relationships between stress ratio with instantaneous elastic strain and plastic strain after different temperature treatment.(a) Relationship between stress ratio and instantaneous elastic strain; (b) Relationship between stress ratio and instantaneous plastic strain.

Figure 10 .
Figure 10.Relationships between stress ratio with visco-elastic strain and visco-plastic strain after different temperature treatment.(a) Visco-elastic strain; (b) visco-plastic strain.

Figure 11 .
Figure 11.Relationships between stress ratio with steady-state strain ratio and damage after different temperature treatment.

Figure 12 .
Figure 12.Relationships between stress ratio and steady-state strain ratio and damage after different temperature treatment.

Figure 15 .
Figure 15.A sketch of crack evolution of pre-cracked sandstone in (a) loading and (b) unloading creep test.

Table 1 .
Basic physical properties of sandstone specimens.

Table 2 .
Applied stress levels in creep tests.