Table of Contents Author Guidelines Submit a Manuscript
Shock and Vibration
Volume 2015 (2015), Article ID 286097, 8 pages
Research Article

Study on Shale’s Dynamic Damage Constitutive Model Based on Statistical Distribution

1State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu 610500, China
2School of Geoscience and Technology, Southwest Petroleum University, Chengdu 610500, China

Received 11 January 2015; Accepted 2 April 2015

Academic Editor: Mohammad Elahinia

Copyright © 2015 Jianjun Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The dynamic constitutive model of shale is basic for shale gas reservoir reforming. In order to investigate the dynamic mechanism of shale, a new dynamic damage constitutive model of shale under uniaxial impact load was established based on the statistical damage theory and the laboratory test results of deformation and the damage characteristics under the action of SHPB impact. Compared with the theoretical results, the model can describe shale’s mechanical attributes and reveal the fracture damage mechanism as well. The results will provide theoretical basis for hydraulic fracturing on shale and other dynamic reforming technics.

1. Introduction

Shale gas is an important unconventional gas resource. Because of the characteristics of low porosity and low permeability, shale gas reservoir must be reformed using hydraulic fracturing or other technics. In order to reveal the law of shale’s mechanical property and fracture damage mechanism, a reasonable dynamic constitutive model is necessary for effective fracturing design needs. But few works are directly related to this content. In recent years, due to the increasing exploitation of shale gas, the study of shale’s mechanical behavior has become the hot issue in the academic fields and engineering areas. On one hand, the research results can provide basic theoretical foundation for reservoir reforming; on the other hand, they can also provide experimental reference and basic parameters for numerical simulation of hydraulic fracturing.

There are many works that have been done on the dynamic constitutive model of rock. For example, Lindholm et al. (1974) proposed an overstress model [1], but the physical parameters of the model are not clear. Sigenori et al. [2] (1977) and Yu [3] further developed the overstress model and obtained the correction of the overstress model, which can describe the dynamic mechanical behavior of rock material very well; Krajcinovic and Silva [4] and Dougill et al. [5] applied damage mechanics in rock and concrete, and the continuous damage process of rock and concrete was analyzed; in addition, the damage constitutive model is established by using the concept of fracture surface. Chen and Feng [6], Wenxue et al. [7], and Shang and Ning [8] established the constitutive relations equation of rock materials by the Lagrange analysis of experimental results under the blasting load. Some results have been achieved on dynamic constitutive model of rock, but the established constitutive model was mostly only for the ordinary rocks whose research backgrounds are engineering blasting, earthquake, and so on and rarely for shale whose mechanic characteristics are more complex than conventional rock in the research background of hydraulic fracturing.

In this paper, the dynamic damage constitutive model of shale under loading is established based on the theory of damage mechanics. At the same time, by considering the relationship between the mechanics parameters and the strain rate, the correction of the model has been put forward. Moreover, compared with the experimental results, the constitutive model is proved to be reasonable.

2. Establishment of Dynamic Damage Constitutive Model of Shale

The large amount of defects within shale makes the quantity of damaged parts increase and the area expand when the shale is under a loading situation, and the increased damage causes the effective stress on a microscale to increase. At the same time, due to the increasing effective stress, more defects are produced. Thus, the strength of the whole rock is reduced. Under the condition of continuous load, the changes of the damage in mechanical properties can be seen as a successive process approximately. In order to analyze the process conveniently, the assumptions are shown as follows [9].(a)Ignore the impact of inertial effect on constitutive relation within strain rate’s scope of 10~102/s.(b)Shale is treated as isotropic viscoelastic on a macroscale. The viscoelastic damage model of shale is shown in Figure 1. In addition, the damage attributes are isotropic as well.(c)Shale obeys Hooke’s law and the unit of shale is viscoelasticity before damage.(d)In a resistivity model, when it is under static load conditions, the viscous body is invalid, while, under the condition of dynamic load, the viscous body plays a role in the process of impact compression.(e)The strength of elastic body obeys the law of probability statistics.On the basis of the above assumptions, the functions of dynamic damage constitutive have been built and, in the functions, the strength of shale obeys Weibull distribution, and the probability density function is given bywhere is the distribution function of shale’s strength, that is, the failure probabilities of shale when the strength of the rock is ; is distribution variable; and are the parameters of Weibull distribution, where represents the macroscopic average strength of rock and reflects the brittleness of rock. Generally speaking, both of them reflect the mechanical properties of rock material.

Figure 1: The viscoelastic damage model of shale.
2.1. The Damage Evolution Equation

The damage of shale is caused by the rising number of internal holes under the action of load. The accumulation of the damage forms the fracture damage of materials. Shale’s breakage that occurs may be regarded as a destruction of microunit. At the same time, the breakage causes the effective area which can bear the load to decrease. On a certain level of load, the failure area is , and the damage variable is defined as the ratio of failure area and total area ; that is, ; at any interval , the failure area is , and when the load is added to , the failure area can be shown asThe damage evolution equation is

2.2. The Dynamic Damage Constitutive Equation

In 1984, the strain equivalence principle was put forward by the French famous scholar Lemaitre [10], which was defined as the strain occurring to damage materials under the stress that was equivalent to the materials in the same condition under the effective stress. And shale’s basic formula of damage constitutive relation can be established bywhere is the nominal stress vector, is the effective stress vector, is shale’s elastic matrix, is the strain matrix, and is the damage variable. The effective stress is based on the assumption that stress cannot transfer through the damage, but, by the previous research results, it can be known that damages like cracks, holes, and joints can transfer part of the stress in the process of compression. Based on the above reasons, the damage correction factor which ranges from 0 to 1 is introduced in the constitutive model. Therefore, formula (4) can be modified as In this paper, shale’s damage criterion is set as random distribution variable and is given byBecause the Drucker-Prager yield criterion has the characteristics of simple form and suitability for rock medium [11], it is determined to express shale’s yield strength, which iswhere is the function of shale’s friction angle and cohesion, and the related equation follows as , with the value of ranging from 0 to ; is the first stress invariant; is the second deviatoric stress invariant; and the formulas, respectively, are expressed as follows:By Hooke’s law, it can be gotten thatBecause of the uniaxial compression in this paper, it can be known thatSo, by formulas (8)–(10), it can be obtained thatAssuming that the shale is under dynamic loads in the -direction, the viscosity of the body in -direction will play a role. And it can be seen that the strain of elastic body is equal to the viscous body , and the sum of the two is the apparent stress value in -direction, which is expressed by the following formula:So, the constitutive model of elastomeric is For viscous part, because there is no damage characteristic, its constitutive relation isBecause of the uniaxial compression, there isSubstitute formulas (13), (14), and (15) into formula (12) and the constitutive relation of shale parallel body is obtained asFormula (16) is shale’s dynamic damage constitutive model based on Weibull distribution.

2.3. Determination of Parameters in Constitutive Model

In the shale’s dynamic damage constitutive model, six parameters, , , , , , and , need to be determined [12]. can be measured by the slope of the initial rise stage of the shale’s stress-strain curve; the value of ranges from 0 to ; and that of ranges from 0.1 to 0.5; the strain and strain rate are measured by experiments, while , , and can be determined by fitting the shale’s stress-strain curves of impact shock experiment.

For the experiment of uniaxial compression of shale, the axial stress is set as and the corresponding strain is set as , and then the shale’s dynamic statistical damage constitutive model can be established byThe formula above can be reduced to Taking the logarithm on both sides of formula (18), there is Taking the logarithm on both sides of formula (20), there isand setEquation (21) can be simplified asBy linear regression, and can be obtained, and then according to is obtained.

The solving process above is to obtain and by using curve fitting method. Then, substitute the values into formula (3) and formula (16), and the shale’s damage evolution equation and statistical damage constitutive model under impact load can be obtained.

Through the SHPB impact test for shale, the obtained test results of the shale’s stress-strain curve are shown in Figure 2.

Figure 2: The stress-strain curves under different strain rates.

Using the above method, in the shale’s dynamic statistical damage constitutive model, parameters and under different strain rates can be obtained, as shown in Table 1.

Table 1: Parameters of the shale’s dynamic statistical damage constitutive model.
2.4. The Correction of Constitutive Model
2.4.1. The Effect of Weibull Distribution Parameters and on the Constitutive Model

In order to study the effect of Weibull distribution parameters and on constitutive model, a univariate method is selected [13]. First of all, fix , and study the changes in the constitutive curve when changes. Then, fix , and study the changes of constitutive curve accompanying the value of . According to the experimental results of 3# test specimen the analysis results are shown in Figures 3 and 4.

Figure 3: The effect of parameter on constitutive curve.
Figure 4: The effect of parameter on constitutive curve.

From Figures 3 and 4, the following rules can be concluded.(1)The shale’s stress-strain curve peak increases with the increase of distribution parameters and .(2)The effect of distribution parameters and on the constitutive curve after the peak is larger than before the peak.

2.4.2. The Effect of Damage Correction Factor on Constitutive Model Curve

According to the experimental result of 7# test specimen, the effect of damage correction factor on constitutive model curve is studied and the result is shown in Figure 5.

Figure 5: The effect of damage correction factor on constitutive curve.

Figure 5 shows that, with the decrease of , the shale’s peak strength increases, the ductility increases after the peak, and the residual strength increases.

2.4.3. The Correction Process of the Constitutive Model

According to the Weibull distribution parameters and obtained above, the strain rate is set as horizontal axis, and and areas are set as vertical axis, respectively. The curves changing with the strain rate are obtained, and the results are shown in Figure 6.

Figure 6: The relationship between model parameters and with strain rate .

By fitting the plot, the relationship between and with strain rate is established and expressed by As shown in Figure 7, the relationship between the damage correction factor and strain rate turns out to be obviously linear and the equation of linear fitting is Substitute formulas (25) and (26) into the constitutive equation (16) and the shale’s dynamic statistical damage constitutive model of specific solutions can be obtained.

Figure 7: Relationship between damage correction factor and strain rate .
2.5. Verification of Constitutive Model

The parameters of the statistical constitutive model can be solved by fitting the stress-strain curves, and then the theoretical curves can be obtained. Compared with the experimental curves, the shale’s dynamic statistical damage constitutive model is verified to be correct and reasonable, and the comparison of theoretical curve and experimental curve is shown in Figure 8.

Figure 8: Comparison of theoretical curve and experimental curve.

From Figure 8, it can be seen that, in this paper, the shale’s dynamic damage constitutive model based on the statistical distribution can better reflect the shale’s stress-strain relationship. Theoretical curves and experimental curves have a good alignment, which illustrates that the constitutive model established in this paper is reasonable.

3. Conclusions

(1)Shale’s stress-strain curves will change in line with the model parameters, while the curves of the peak strength increase with the Weibull distribution parameters and increasing and decrease with the damage correction factor increasing.(2)Parameters and of the statistical damage constitutive model and strain rate abide by a linear relationship, while and strain rate abide by the quadratic polynomial relationship.(3)The statistical damage constitutive model is proved to be reasonable by the comparison of theoretical curves and experimental curves.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The research is financially supported by the National Science Foundation of China (Grant no. 51174170).


  1. U. S. Lindholm, L. M. Yeakley, and A. Nagy, “The dynamic strength and fracture properties of dresser basalt,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 11, no. 5, pp. 181–191, 1974. View at Publisher · View at Google Scholar · View at Scopus
  2. K. Sigenori, S. Azuhiko, and K. Minoru, “On the mechanical behavior of rocks under impulsive loading,” Bulletin of the Faculty of Engineering Hokkaido University, vol. 83, pp. 51–62, 1977. View at Google Scholar
  3. Y. Yu, Rock Dynamics, University of Science and Technology Beijing Press, Beijing, China, 1994.
  4. D. Krajcinovic and M. A. G. Silva, “Statistical aspects of the continuous damage theory,” International Journal of Solids and Structures, vol. 18, no. 7, pp. 551–562, 1982. View at Publisher · View at Google Scholar · View at Scopus
  5. J. W. Dougill, J. C. Lau, and N. J. Burt, “Toward a theoretical model for progressive failure and softening in rock, concrete and similar material,” in Proceedings of the ASCE-EMD Specialty Conference on Mechanics in Engineering, pp. 335–355, Waterloo, Canada, May 1976.
  6. Y. Chen and S. Feng, “Lagrangian analysis method and constitutive relation of cement stone on intense dynamic loading,” Chinese Academy of Railway Sciences, vol. 17, no. 3, pp. 1–11, 1996. View at Google Scholar
  7. G. Wenxue, Y. Jun, and H. Fenglei, “The constitutive relation of rock under strong impact loading,” Journal of Beijing Institute of Technology, vol. 20, no. 2, pp. 165–170, 2000. View at Google Scholar
  8. L. Shang and J.-G. Ning, “Dynamic constitutive relationship of concrete subjected to shock loading,” Engineering Mechanics, vol. 22, no. 2, pp. 116–120, 2005. View at Google Scholar · View at Scopus
  9. X. Jinyu, F. Jianshe, and L. Xiaocong, Dynamic Mechanical Properties of Rock under the Condition of Confining Pressure, Northwestern Polytechnic University Press, 2012.
  10. J. A. Lemaitre, “How to use damage mechanics,” Nuclear Engineering and Design, vol. 80, no. 2, pp. 233–245, 1984. View at Publisher · View at Google Scholar · View at Scopus
  11. J. Liu, L. Lin, R. Song, and J. Zhao, “A pore scale modeling of fluid flow in porous medium based on navier-stokes equation,” Disaster Advances, vol. 6, supplement 1, pp. 5–11, 2013. View at Google Scholar
  12. Q. He, L. Dianshu, and W. Shuren, “Experimental study of rock's constitutive relation on high strain ration,” Blasting, vol. 13, no. 2, pp. 1–7, 1996. View at Google Scholar
  13. W. Cao, M. Zhao, and C. Liu, “Study on the model and its modifying method for rock softening and damage based on Weibull random distribution,” Chinese Journal of Rock Mechanics and Engineering, vol. 23, no. 19, pp. 3226–3231, 2004. View at Google Scholar · View at Scopus