Abstract
Time-dependent rheological response of geomaterials to thermal treatment is a crucial issue in geothermal energy utilization and deep mineral mining. This response, however, has not yet been fully considered in the existing rheological constitutive models for sandstones. In order to experimentally investigate such responses and establish the associated rheological constitutive model, this study considers the sandstone specimens which have been thermally treated under different temperatures. The triaxial rheological test in conjunction with the scanning electron microscope is employed in the investigation to observe the mechanically and macro-/micromorphologically rheological response. Investigation results show that the thermal treatment induces microcracks and microdefects, and subsequently, they propagate during the creep. As a consequence, the heterogeneous deformation occurs, and macrocracks are present, leading to the irregular fluctuation and mutation in strain over time. A higher temperature contributes to a more severe structure damage and in turn reduces the intactness of sandstones and elevates the rheological response. The investigation allows successful establishment of a three-dimensional constitutive equation considering the instantaneous elastic response to thermal treatment. Based on the equation, a nonlinear visco-elastoplastic rheological constitutive model is developed for sandstones. Comparison with three existing rheological models shows that the model developed in this study could well represent the rheological process of the thermally treated sandstones.
1. Introduction
Thermal treatment plays an important role in the mechanical behaviors of rocks [1, 2]. It has also widely been reported that the engineering practices, such as the geothermal energy utilization and deep mineral mining, usually involve the rheological issue of thermally treated rocks [3, 4]. Because deep rocks are usually conditioned under high stress and long-term seepage pressure coupled with high temperature, it is of great interest to understand the rheological behavior of deep rocks in response to high-temperature thermal treatment [5, 6].
Recent advances in understanding the response of rocks or rock-like media to high-temperature thermal treatment include studies by Gautam et al. [7], Fan et al. [8], Yang et al. [9], Ersoy et al. [10], and Peng and Yang [11]. Among them, Gautam et al. [7] experimentally revealed that the elastic modulus of sandstones drops by 78% due to thermal treatment when the temperature rises from room temperature (25°C) to 650°C. Fan et al. [8] experimentally found that the temperature rise causes the decrease in wave velocity and Young’s modulus of granites for the reason of thermal treatment. Yang et al. [9] observed the thermal damage of granites result from high temperature and showed the complicated response of multiple mechanical parameters (including strength, static elastic modulus, static Poisson’s ratio, peak axial strain, dynamic elastic modulus, and dynamic Poisson’s ratio) to the thermal treatment. Specifically, as temperature increases, the first two parameters rise, followed by their falling, while a falling in the static Poisson’s ratio is followed by its sharp rise, a continuing increase in the peak axial strain, and no significant change in the dynamic Poisson’s ratio. As compared to the previous experiments, Ersoy et al. [10] performed the uniaxial compressive strength test of basic volcanic rocks under a wider range of temperature (200 to 1,000°C) and then revealed that the increasing temperature induces the thermal damage and in turn causes the strength loss. A wider range of rock types (sedimentary, metamorphic, and volcanic) was investigated in Peng and Yang [11], in which the continuing decrease in wave velocity, compressive strength, and Young’s modulus is observed as temperature increases, despite an exception occurrence at the three parameters of the quartz sandstones (increase initially and then decrease) due to high porosity. Numerous previous attentions have been paid to the response of various mechanical behaviors to thermal treatment, yet particular concerns are seldom placed on the rheological response. Investigations of the rheological response to thermal treatment, such as Belmokhtar et al. [5], are rare.
On the contrary, an increasing number of constitutive models have been developed to describe rock rheological behaviors in the past years. Chen et al. [12] reported a damage mechanism-based rheological model for granites with a consideration of thermal treatment. Xiong et al. [13] presented a unified thermo-elasto-viscoplastic rheological model for soft rocks in the critical state. Zhao et al. [14] developed a nonlinear elasto-viscoplastic rheological model with consideration of various deformation components. In the model, a parallel combination of Hooke and plastic slide bodies is incorporated, in conjunction with a Hooke body, a Kelvin body, and a generalized Bingham body. On the basis of variable-order fractional derivatives and continuum damage mechanics, Tang et al. [15] gave a four-element rheological model for salt rocks. Although continuing efforts have been devoted to the improvement of rock rheological constitutive models, the existing rheological constitutive models considering the rheological response to thermal treatment, such as the model of Chen et al. [12] for granites and the model of Yang et al. [16] for sandstones, are rare.
This work adds to the limited body of (1) investigation on the potential rheological response of sandstones to thermal treatment and (2) development of a rheological constitutive model for sandstones with proper consideration of such response. Initially, the sandstone specimens for the investigation are collected from the Majiagou rock slope located in the Three Gorges Reservoir area, in central China. Then, the specimens, separately, are thermally treated under different temperatures, prior to the triaxial rheological test for investigating the mechanically rheological response. In the test, the scanning electron microscope is employed to additionally investigate the micromorphologically rheological response. Based on the investigation of rheological response to thermal treatment, a three-dimensional constitutive equation is established by incorporating such response. A nonlinear visco-elastoplastic rheological constitutive model adopting this established equation is then developed for sandstones. At last, the estimates from the developed model are compared with that from the existing models and the laboratory observations to evaluate the performance of the developed model.
2. Laboratory Investigation of Rheological Response of Sandstones to Thermal Treatment
2.1. Materials and Methods
2.1.1. Preparation of Specimens and Description of Testing Instrument
The rock specimens are steel-gray fine sandstones. Their physical parameters are listed in Table 1. The cylindrical specimens to be tested were prepared with dimensions of Φ = 50 mm and H = 100 mm.
The testing instrument, namely, triaxial rheometer, composed of the control system, oil source, axial compression system, confining pressure system, and seepage system along with sensors (shown in Figure 1), was employed. More details about the rheometer are available in Xu et al. [17]. In addition to the triaxial rheometer for observing the mechanically rheological response, a Quanta200 environmental scanning electron microscope is employed during the test to observe the micromorphologically rheological response.

2.1.2. Procedures of the Test
The following steps were adopted during the test:(1)Select the rock specimens with close longitudinal wave velocities for the triaxial rheological test.(2)Thermally treat the specimens under temperatures of 25°C (room temperature), 200°C, 400°C, and 600°C, separately, for 6 hours (h), prior to cooling down to the room temperature (shown in Figure 2). Due to time consumption of test, one specimen alone is prepared for each temperature above.(3)Saturate the specimens for 24 h with the aid of the vacuum pump to ensure the full saturation before the triaxial rheological test.(4)Conduct the triaxial rheological test using the triaxial rheometer. According to the field geostress, the confining pressure adopted in the triaxial rheological test is set at 6 MPa and the seepage pressure is set at 4 MPa. Deviatoric stress (the difference between axial pressure and confining pressure) loading is then imposed by grade. The loading duration under all levels of deviatoric stress is limited to be no less than 72 h, and the rate of deformation is limited to be no more than 0.001 mm/24 h. When the deformation becomes stable under the current level of deviatoric stress, then the next level is imposed. The test stops once the rock specimen is damaged. Note that, for sufficient contact between the specimen and rheometer, prior to graded loading, the preloading has been performed.

(a)

(b)

(c)

(d)
2.2. Investigation Results
2.2.1. Macromorphologically Rheological Response to Thermal Treatment
Figure 3 compares the macromorphology of 25, 200, 400, and 600°C thermally treated specimens after rheological failure. The 25°C thermally treated specimen after rheological failure presents an axial failure plane and some cracks. The cracks are primarily tensile ones, supplemented by shear ones.

(a)

(b)

(c)

(d)
In comparison, the 200°C and 400°C thermally treated specimens after rheological failure presents a pair of X-shaped failure planes and some cracks. The cracks are primarily shear ones, supplemented by tensile ones. And the 600°C thermally treated specimen after rheological failure is crushed, and it presents an oblique failure plane. Numerous shear cracks are observed, and a large amount of powder is produced.
The comparison shows that the sandstones exhibit significant macromorphological rheological response to thermal treatment. A higher temperature contributes to a more severe structure damage and in turn reduces the intactness of sandstones.
2.2.2. Micromorphologically Rheological Response to Thermal Treatment
Figure 4 compares the microscopic morphology of failure planes of 25, 200, 400, and 600°C thermally treated specimens. The failure plane associated with 25°C is step-shaped, with each step being neat and smooth (shown in Figure 4(a)). The heights of steps vary from a few µm to 30 µm. In contrast, the failure planes associated with 200, 400, and 600°C are rough and irregularly shaped (shown in Figures 4(b)–4(d)). Loose structure, small sized particles, and holes are revealed.

(a)

(b)

(c)

(d)
This could be due to that thermal oxidation decomposition occurred under high temperature, and gas was produced. As a consequence, void was yielded within the sandstones. It could be hence inferred that a higher temperature renders a more severe rheological damage to sandstones.
2.2.3. Mechanically Rheological Response to Thermal Treatment
An example of the axial strain of thermally treated specimens over time is shown in Figure 5, accompanied with the change in instantaneous elastic modulus with the temperature used in thermal treatment. The instantaneous elastic modulus is determined using the following equation [9, 18]:where is the instantaneous elastic modulus, is the axial stress, is the confining pressure, is the instantaneous stain, and is the bulk modulus.

(a)

(b)
In Figure 5(a), the strain sometimes exhibits slight fluctuation. As revealed in Figure 5(b), the instantaneous elastic modulus is significantly subject to the temperature used in thermal treatment. A higher temperature contributes to a smaller instantaneous elastic modulus value.
Figure 6 selects the axial strain and its rate over time during the last level of deviatoric stress loading. This figure shows that the sandstone experiences the following three rheological stages ( is the axial strain, and is the second derivative of ):(i)The deceleration rheological stage ()(ii)The stable rheological stage ()(iii)The acceleration rheological stage ()

3. Developed Model considering Rheological Response to Thermal Treatment
Based on the investigation above of rheological response to thermal treatment, a three-dimensional constitutive equation is established here by incorporating such response. A nonlinear visco-elastoplastic rheological constitutive model adopting this established equation is then developed for sandstones.
3.1. Development of the Model
3.1.1. Identification of Instantaneous Elastic Response
Instantaneous elastic response of rocks to thermal treatment is normally quantified by an instantaneous elastic index, , as defined by Kachanov [19]:where is the instantaneous elastic modulus before thermal treatment and is the instantaneous elastic modulus after thermal treatment.
The ranges between 0 and 1. The investigation data in Section 2.2 give rise to the relation between and the temperature () used in thermal treatment (Figure 7). As shown in this figure, increases with . A higher value tends to result in a larger value, signifying an elevated instantaneous elastic response of sandstones to thermal treatment. Fitting the data by the least square method gives

3.1.2. Identification of Long-Term Rheology
Long-term rheological index, , as defined by Zhang et al. [20], is introduced to quantify the long-term rheology of sandstones:where is the initial rheological viscoelastic coefficient or instantaneous elastic modulus of rock, is the failure time point, and is the rheological viscoelastic coefficient at .
tends to be a stable value () under long-term seepage pressure. As stated by Zhou et al. [21], rock failure occurs in the situation where deviatoric stress is greater than the yield stress. Accordingly, equivalently denotes the time point when the deviatoric stress starts to exceed the yield stress. Thus, equation (4) is rewritten as follows [20]:where is a coefficient related to damage tolerance.
3.1.3. Consideration of Instantaneous Elastic Response to Thermal Treatment and Long-Term Rheology into the Model
From Figures 5(a) and 6, the instantaneous elastic deformation (i.e., deceleration rheological stage) of the thermally treated sandstone occurs once the axial stress is loaded, indicating that the rheological model should contain a series-wound elastic element. In the stable rheological stage, the strain under a low level of deviatoric stress tends to be constant, which suggests that the rheological model should include an adhesive element and an elastic element in parallel. In comparison, under a high level of deviatoric stress, the strain increases with time, which implies that a plastic component should be included in the rheological model. In the acceleration rheological stage, the deformation increases sharply with time, indicating nonlinearity. It suggests that the rheological model should include a nonlinear viscous element.
Based on the analysis above, a nonlinear visco-elastoplastic rheological model, as shown in Figure 8, is developed. The model is composed of (1) a Hook element, (2) a Kelvin–Voigt element, (3) a visco-elastoplastic body, and (4) a nonlinear viscous body. The instantaneous elastic index D defined by Kachanov [19] has been introduced into model so as to consider the rheological response of sandstones to thermal treatment. Also, the long-term rheological index defined by Zhang et al. [20] has been introduced to consider the long-term rheology.

The nonlinear viscous body used in the developed model is the parallel combination of a nonlinear Newton body and a plastic body. The constitutive equation for the Newton body can be expressed as follows:where is the total stress, is the total strain, and is the viscous coefficient.
Denotewhere is the coefficient of the acceleration rheological initial viscosity and is determined by the rheological curve fitting.
Substituting equation (7) to equation (6) gives
Elastic modulus, , can be expressed as
Substituting equation (5) to equation (9) gives
Similarly, viscous coefficient, (for the constitutive bodies other than the Newton body), can be expressed aswhere is the initial viscosity coefficient.
Substituting equation (5) to equation (11) gives
In the conventional triaxial test, the stress status is . The stress status and load of each level are constant after loading. Accordingly, the three-dimensional rheological equation of nonlinear visco-elastoplastic rheological model which considers the rheological response of sandstones to thermal treatment is shown as follows:where is the bulk modulus and , , and are the shear modulus, respectively, corresponding to , , and in the three-dimensional stress state. The four parameters , , , and can be determined as follows [22]:
As mentioned in Section 2.1.2, prior to graded loading, the preloading has been performed. As is known, the preloading would result in some creep. This creep resulting from the low preloading pressure acted, however, would be rather small so that the influence of such creep could be even ignored [23, 24]. For this reason, the creep due to preloading is not considered into the developed model.
3.2. Performance of the Developed Model
In this section, the estimates from the developed model are compared with that from the existing models and the laboratory observations to evaluate the performance of the developed model. The existing models selected to compare includes the following:(i)The Kelvin–Voigt model [25, 26] (Figure 9(a))(ii)The Burgers model [27, 28] (Figure 9(b))(iii)The Nishihara model [29, 30] (Figure 9(c)).

(a)

(b)

(c)
In order to obtain the parameter values of the developed model and the existing models, the rheological test data are processed through the Boltzmann superposition principle (Sone and Zoback [31]) first. Then, the data are fitted to the models, with the fitted model parameter values listed in Table 2. The developed model and the existing models adopting the respective parameters are thus established.
The estimates from the developed model and the laboratory observations are compared in Figure 10. The comparison shows that the estimate from the developed model is consistent with the sandstone triaxial rheological test result.

The estimates from the developed model and that from the selected existing models are compared in Figure 11. Figure 11 demonstrates that the developed model better reflects the rheological process of the thermally treated sandstones in comparison with the Kelvin–Voigt model, the Burgers model, and the Nishihara model.

4. Conclusions
This work adds to the limited body of (1) investigation on the potential rheological response (mechanically and macro-/micromorphologically) of sandstones to thermal treatment and (2) development of a rheological model for sandstones with proper consideration of such response.
Investigation results show that the thermal treatment induces microcracks and microdefects, and subsequently, they propagate during the creep. As a consequence, the heterogeneous deformation occurs and macrocracks are present, leading to the irregular fluctuation and mutation in strain over time. A higher temperature contributes to a more severe structure damage and in turn reduces the intactness of sandstones and elevates the rheological response.
Based on the investigation of rheological response to thermal treatment, a three-dimensional constitutive equation is established by incorporating such response. A nonlinear visco-elastoplastic rheological model adopting this established equation is then developed for sandstones. In comparison with the existing models, the developed model shows a better fit with the triaxial rheological test observations, revealing the capability to represent the rheological process of the thermally treated sandstones. The developed model can serve for numerical simulation as the constitutive model.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The research was supported by two funds, namely, the Open Fund Project of Key Laboratory for Geo-Hazards in Loess Area of China (Grant no. KLGLAMLR201505) and the Open Foundation of Jiangxi Engineering Research Center of Water Engineering Safety and Resources Efficient Utilization (Grant no. OF201602).