Permeability is one of the most fundamental reservoir rock properties required for modeling hydrocarbon production. However, shale permeability is not yet fully understood because of the high temperature of shale reservoirs. The third thermal stress that is caused by temperature change will decrease the permeability of shale. In this work, a theoretical model has been derived to describe the permeability of shale considering the third thermal stress; the principles of thermodynamics and the mechanics of elasticity have been employed to develop this model. The elastic modulus parameters of the shale were measured, along with Poisson’s ratio, as required. Lastly, the permeability of shale was tested by transient pulse-decay. Isothermal flow experiments were carried out at 303, 313, 323, and 333 K to assess the effects of shale expansion and deformation on shale permeability caused by the third thermal stress. The permeability of shale samples, as predicted by the model, was found to agree well with experimental observations. The model may provide useful descriptions of the gas flow in shale. The correction accuracy of the permeability was found to increase at lower permeability. However, the development of completely predictive models for shale permeability will require additional experimental data and further testing.

1. Introduction

Shale gas is an unconventional but promising gas resource that has been used with significant success in the world. Permeability is one of the most fundamental properties of any reservoir rock, and it is required for modeling hydrocarbon production. Heller et al. [1] measured the permeability of the Barneet and Eagle Ford shales at low pressure (6.9 MPa) and reported the influence of confining pressure and pore pressure on permeability. Javadpour et al. [2] obtained nanoscale images of shale porosity using an atomic force microscope and then presented a model for gas flow in the nanopores of mudrocks based on the theory of Knudsen diffusion and slip flow. Kwon et al. [3] investigated the effect of effective stress on shale permeability using Wilcox shale. Chalmers et al. [4] measured the permeability of gas shale from the Liard Basin and assessed the effects of origin and distribution, total organic carbon (TOC) content of minerals, porosity, and effective stress on the permeability of the matrix. Mahnaz et al. [5] used helium as a medium to measure the permeability of shale and assessed the Klinkenberg effect on the permeability of gas shale. The permeability of rock change by the wetting-drying cycles and freeze-thaw-fatigue loading was reported by Song et al. [6, 7].

The thickness of gas shales varies widely. For example, a study of five shale gas systems in the USA indicates a range in shale thicknesses of 2 to 700 m, and the maximum reported thickness of Chinese gas shales is 925 m [810]. Therefore, geothermal differences associated with depth will change the permeability of the shale. Zhang et al. [11] measured the permeability of Carrara marble at different temperatures. They found that marble penetration was significantly increased when the temperature rose to 600-700 K. Li and Xian [12] carried out coal permeability experiments under different stress and temperature conditions using methane and helium. Liang et al. [13] measured the permeability of coal samples from Bang Xin Zu. They concluded that there is an exponential relationship between the permeability of coal and the temperature. Feng et al. [14] measured the permeability and deformation of the high-rank coal while increasing the temperature from 298 to 873 K.

Current studies mainly focus on the temperature effects on the permeability of rock. However, temperature effects on the permeability properties of rocks containing complicated components such as shale are comparatively unknown. Thermal stress can be divided into three categories [15]: (1) the first thermal stress, which is produced by the external deformation of homogeneous materials by heating; (2) the second thermal stress, which is produced as a result of the nonuniform temperature distribution within the same object that cannot freely deform; and (3) the third thermal stress, which occurs in materials comprised of multiple components with different linear expansion coefficients that produce thermal stress due to differences in expansion coefficients as the temperature changes. The thermal stress produced in a shale under the impacts of temperature generally consists of the third thermal stress. The reason for this is that shale accumulations in China are generally composed of organic matter and minerals (quartz, calcite, and clay minerals), and the linear expansion coefficients of these minerals do not have the same magnitude [16, 17].

In this paper, a model is proposed that is based on the principles of thermodynamics and mechanics of elasticity, to describe the temperature impacts of shale permeability under the third thermal stress conditions. The correction factor is tested by experimental measurements on the permeability of shale at 303, 313, 323, and 333 K.

2. Model

2.1. Calculation of the Third Thermal Stress

Since shale is composed of different mineral components, the linear expansion coefficient of shale is calculated considering the linear expansion coefficients of the constituent minerals and organic matter:

The thermal stress caused by temperature is expressed as follows [18]: where is the percentage of the total volume of the various materials, is the linear expansion coefficient of a particular material, the volume expansion coefficient is , is the elasticity modulus, is the Poisson’s ratio, and is the change in temperature.

Thermal stress is a vector and its direction points from the region of largest expansion quantities to the region of smallest expansion, as shown in Figure 1.

The direction of the thermal stress in bulk shale points from the organic matter to the minerals. Assuming that the minerals are fully mixed with the organic matter, the value of the thermal stress in small units can be calculated by the model of Figure 1. The thermal stress is expressed as follows: where is the linear expansion coefficient of mixture material and is the linear expansion coefficient of small materials.

2.2. Balance Differential Equation for the Thermal Stress Field

Temperature change in shale formations can be determined by the depth of the shale. A shale formation at the same depth can be considered to have isothermal conditions because the temperature would not change at the same depth of shale burial. Additionally, shale can be assumed to be an isotropic solid for calculation simplicity [19]. Therefore, a planar coordinate system is used to establish differential governing equations for the temperature and stress fields (Figure 2). The temperature is determined by , where is the geothermal gradient, K/100 m.

The equilibrium, geometric, and physical equations in Cartesian coordinates are as follows [18]:

The equilibrium differential equation, using the displacement as the basic unknown function, is obtained through a transformation of Equation (4).

The solution to Equation (5) is obtained by introducing the thermal elastic strain potential function (,).

A set of solutions to the equation of state of constant temperature, homogeneous plus, is the general solution of differential equations of equilibrium equations, boundary conditions to obtain simultaneous, displacement solution for determining: where is the vertical crustal stress, is the horizontal stress, and is the buried depth.

2.3. Correction Coefficient for Thermal-Hydro-Mechanics-Coupled Permeability

The shale strain in the and axes (horizontal plane) are the same, because the temperature of the shale is equal at the same depth. So, the shale volume strain can be formulated as

The empirical formula can be used to calculate the permeability of the homogeneous medium under the fluid solid heat coupled conditions, where is the porosity, and is the specific surface area, which is the ratio of value of pore surface area and the whole inner volume.

By introducing a modification coefficient to consider the third thermal stress, the equation can be obtained, where is the permeability without considering the third thermal stress.

The shale porosity is not affected by the third thermal stress, so . For the linear elastic strain of a homogeneous medium, an approximation of the internal surface area change is made, which is expressed as .

The above equations are simplified by omitting higher order terms because the linear elastic deformation is small:

External factors that affect the correction coefficient of the shale permeability include the temperature difference, geothermal gradient, vertical stress, and horizontal stress. These factors are also related to the medium’s elastic modulus, Poisson’s ratio, coefficient of expansion, and the mass ratio of its mineral components.

Equation (9) shows that permeability decreases as temperature increases. Permeability decreases largely as a result of increasing the temperature when the difference in the expansion coefficient of shale components is higher and the organic matter content is higher.

2.4. The Permeability Test

It is very difficult to measure shale permeability using conventional steady-state measuring methods because the shale permeability is very low (10-3 to 10-6 mD). So, an alternative method such as the transient pulse-decay (TPD) method [20] or the beam-bending method [21] was proposed to measure the shale permeability. In contrast with conventional steady-state permeability test methods, the pulse-decay method does not need to record rock outlet velocity. It has the advantage of high measurement accuracy and efficiency, but it has three disadvantages: (1) difficult sample preparation, (2) high sealing of the experimental device, and (3) a requirement for high-precision pressure sensors. Measurement precision with the beam-bending method is higher than with the conventional steady-state permeability test method. A permeability range of 0.009 to 400 nD can be measured by this method, but it also has three disadvantages: (1) complex mathematical models, (2) difficult sample cutting, and (3) small samples with poor representation.

In this work, the transient pulse-decay experiment was used to measure the permeability of shale samples affected by the temperature.

3. Experiment and Methods

3.1. Samples

The shale samples used in the laboratory experiments were collected from the lower Silurian Longmaxi Formation in Yibin City, Sichuan Province. The total organic carbon (TOC) content of the samples is 7.88%, and the vitrinite reflectance () of the samples is 2.85%. These values are in the range of most favorable conditions for the occurrence of shale gas () [22]. Shale samples were prepared with rock specimen blocks were obtained (Figure 3). Subsequently, the specimens were put into the oven and baked for 24 h at a temperature of 353 K. Then, the specimens were cooled and sealed in plastic bags.

The mineralogical compositions of shale samples were tested using X-ray diffraction. The X-ray diffraction data were obtained using a Siemens D5000 diffractometer, using CuKα radiation. The samples were scanned from 2 to 75° 2θ, with a step time of 2 s per 0.02° step. The results are shown in Table 1.

The triaxial compression strength, elastic modulus, and Poisson’s ratio of the samples were tested using MTS815. The results are shown in Table 2.

3.2. Apparatus

Figure 4 presents a schematic diagram of the permeability measurement apparatus. This apparatus consisted of a temperature control system, an experimental system, a data acquisition system, an auxiliary system, and a triaxial pressure chambers.

The maximum oil field thermostat system can be heated to 473 K, with a temperature control accuracy of 0.1 K. The servo loading system provides a maximum 3000 kN axial load to provide a maximum 200 MPa of confining pressure.

The experimental system consisted of a sample cell, two reference cell (5 mL), two pressure transducer (Rosemount, American, model 3051, range: 0 to 13.79 MPa, accuracy: 8 kPa), and a pressure difference sensor (Rosemount, American, model 3051, range: 0 to 0.68 MPa, accuracy: 0.68 kPa). The pore pressure is provided by the ISCO pump, with a maximum available gas pressure of about 5000PSI.

3.3. Procedure

Prior to measuring the permeability of shale samples, the apparatus was checked for pressure leaks. The helium was injected into the reference cell. If the values of pressure collected from the pressure transducer decreased less than the accuracy of the pressure transducer (8 kPa) in 24 h, then it was no leaks in the system [23].

The sample was subsequently installed in the test cell. The system was allowed to reach thermal equilibrium at the target temperature, based on heating in the oil field, ensuring that any temperature-induced expansion took place. The swelling of the shale in response to increased temperature was considered to have reached its maximum when the strain gauges showed the strain of the sample had plateaued, at which point a vacuum was applied.

We apply a confining pressure of 9 MPa, then open the relevant valve of the pulse permeability instruments, and the was used to up to pressure of about 7 MPa to make the net confining pressure on the rock sample to 2 MPa. After the pressure of the instrument’s upstream chamber, rock sample pores, and downstream chamber, it reaches thermal equilibrium according to the instrument prompts, manually lowering the pressure in the downstream chamber by about 70 k Pa to form the initial attenuation pressure pulse.

4. Results and Discussion

4.1. Strain Caused by Gas

Kumar et al. [24] measured the permeability of shale samples with helium, methane, and carbon dioxide. They found that the permeability of samples decreased by 50% using carbon dioxide and decreased by 20% using methane. It is obviously that the strain on the shale caused by gas could affect the permeability of shale. According to the report by Ao et al. [25, 26], there are two effects that induce strain in the sample caused by gas. Gas was slowly adsorbed by the sample, and so the surface potential of the sample was reduced. The volume was also increased with increases in the thickness of the surface layer while free gas penetrated the pores and cracks. In addition, the gas infiltrated into cracks that were greater than or equal to the molecular size of the gas molecules. In this study, the was used to measure the permeability. The have not obviously adsorption ability, so the adsorption-induced strain on shale could not occur. In addition, the change in effective stress caused by gas pressure does not occur instantaneously. The experiment uses the transient method to measure the permeability, so the test time is very short. Therefore, it is difficult for the gas pressure to produce effective stress changes in this experiment. To sum up, there is no strain caused by the gas.

4.2. Experimental Results of Shale Permeability

The permeability was calculated as where is the permeability, mD; is the compressibility of the pore fluid of the rock sample, dimensionless; μ is the gas viscosity, mPa⋅s; is the porosity of shale, dimensionless; is the length of shale, mm; is the slope of the pressure difference between the upstream and downstream chambers and time of the pulse permeability tester tested in the single logarithmic graph, dimensionless; is the ratio of the pore volume of the test rock sample to the volume of the upstream and downstream chambers of the pulse permeability meter, respectively; when the value of is equal to , the amount of is 1.17. The results of the shale permeability determinations are shown in Table 3.

4.3. Verification of the Permeability Prediction Model

In this work, we have established a model to predict shale permeability. Eight parameters are included in this model, i.e., temperature difference , geothermal gradient , vertical stress , horizontal stress , the elastic modulus , Poisson’s ratio , the thermal expansion coefficient , and the proportions of mineral components . was determined based on the results reported by Rahman [27], who pointed out that the TOC of shale is about 10%, whereas the mineral content of shale is about 90%. The parameters used in the model are shown in Table 4.

The correction coefficient for the temperature difference was calculated using equation (9) and reported in Table 5.

To eliminate the effect of the second thermal stress on the results, and based on the experimental conditions, the permeability at 303 K was considered the initial permeability of the model. Therefore, only the third thermal stress will affect the results. Figure 5 presents the model representations for the shale permeability.

Figure 5 compares the experimental and modeled shale permeability values. Note in the figure that the permeability predicted by the model is in good agreement with the measured data. In addition, the permeability was found to decrease as temperature increased. The model and experimental values were closer at lower permeabilities. This indicates that for a mixture of a hypotonic medium with a low permeability, the effect of the third thermal stress is higher as a result of the reduction in the permeability.

5. Conclusions

This paper presents a theoretical model to describe the impacts of temperature on shale permeability under the third thermal stress condition. This model was built on the principles of thermodynamics and mechanics of elasticity. All parameters included in this model have physical meanings, which provide a predictive basis for constraining the shale permeability.

The main constituents in the shale sample are minerals and organic matter. A shale sample will produce uneven swelling because the thermal expansion coefficients of the minerals and organic matters are different as a result of the effects of temperature change. The third thermal stress is caused by uneven swelling. Furthermore, shale permeability will change under the effects of matrix swelling, which is caused by both ground and thermal stresses. We have also modified the formula for the correction coefficient of the shale permeability. The shale permeability predicted by the model coincides with the experimentally measured permeability, which decreases as the temperature rises. Additionally, the accuracy of the permeability corrections was found to be higher when the permeability was reduced. The shale permeability decreased more rapidly when the temperature increased, particularly while the difference in the shale component expansion coefficients was larger and the organic matter content was higher.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This study was funded by the National Key Basic Research Program of China (No. 2014CB239206), the Natural Science Foundation of Chongqing (Nos. cstc2019jcyj-msxmX0507 and cstc2019jcyj-zdxmX0024), and the Scientific and Technological Research Program of Chongqing Municipal Education Commission (KJQN20200).