Research Article  Open Access
A New Equivalent Statistical Damage Constitutive Model on Rock Block Mixed Up with Fluid Inclusions
Abstract
So far, there are few studies concerning the effect of closed “fluid inclusions” on the macroscopic constitutive relation of deep rock. Fluidmatrix element (FME) is defined based on rock element in statistical damage model. The properties of FME are related to the size of inclusions, fluid properties, and pore pressure. Using FME, the equivalent elastic modulus of rock block containing fluid inclusions is obtained with Eshelby inclusion theory and the double MT homogenization method. The new statistical damage model of rock is established on the equivalent elastic modulus. Besides, the porosity and confining pressure are important influencing factors of the model. The model reflects the initial damage (void and fluid inclusion) and the macroscopic deformation law of rock, which is an improvement of the traditional statistical damage model. Additionally, the model can not only be consistent with the rock damage experiment date and threeaxis compression experiment date of rock containing pore water but also describe the lockedin stress experiment in rocklike material. It is a new fundamental study of the constitutive relation of lockedin stress in deep rock mass.
1. Introduction
As the research of rock mechanics gradually develops to deep and complicated geological conditions, the traditional theory of rock mechanics has been continuously improved and perfected. The application of CT technology in rock mechanics helps people have a new understanding of microscopic pore structure of rock, thus establishing a microscopic system of rock [1–4].
Pores are classified into connected pores and closed pores. Actually, the theoretical and experimental research on the connected pores in rock is relatively mature [5–8], and the theory about it has developed to fluid interaction with all forms of subsurface materials, whether the materials are unconsolidated or crystalline [9]. From the perspective of engineering applications, GassmannBiot equation describes the relationship between rock physical properties and pore fluid characteristics under the conditions of low frequency [10]. Besides, it is an important theoretical basis of rock fluid replacement or seismic wave detection in oil and gas engineering [11–15]. Compared with the connected pore problems, theoretical and experimental investigations about closed pores are relatively few. Closed pore is an important carrier of lockedin stress. In 1979, the concept of “lockedin” stress was first proposed by Zongji Chen. Based on the microscopic rheological analysis of rock, Chen found that the tectonic loading and thermal loading lead to nonuniform stress field, and some of them are retained in the form of lockedin stress. Additionally, lockedin stress is considered as an important cause of engineering disaster [16]. Barrows considered that the internal shape of viscoelastic Earth material possesses an internal equilibrium pressure over extended periods of geologic time. Internal equilibrium pressure should be regarded as an intrinsic physical property of Earth materials [17]. For a long period of time, the lockedin stress hypothesis has only a few quantitative studies and developments, but the situation has changed in recent years. Wang analyzed the complex environment of diagenesis and the influence of geological processes on the inhomogeneity and discontinuity of rock based on the geological characteristics of rocks [18]. It is also an important source of lockedin stress. Using the hypothesis of lockedin stress, Yue proposed and demonstrated that highpressure fluid inclusions are concrete, measurable, and computational stress inclusions. Besides, they are prevalent in geological rocks and minerals. The effect of fluid inclusions on surrounding rock is calculated. Yue believed that small, confined, and compacted fluid inclusions are common volume stress of rock burst, surrounding rock rupture, and large deformation of roadways in deep rock excavation chambers [19–21].
In general, the research on the constitutive relations of rock fluid inclusions under deep high pressure is very few. Therefore, this study attempts to establish a concise mathematical and physical model that reflects the constitutive relationship between lockedin stress and rock block. The relationship between microcosmic characteristics and macroscopic rules has good applicability, which also lays a good foundation for the mechanical problems of deep rock mass containing closed fluid inclusions.
2. Related Work
The statistical damage model is a good mathematical tool for establishing the constitutive relationship between the microscopic and macroscopic properties of rocks. “Rock element” is the basic element of the rock statistical damage model and the “rock element” satisfies the DP yield criterion. Rock elements gradually yield under external stress, and their strength intensity conforms to Weibull distributions. Moreover, the process can well reflect the stressstrain law of the rock under external stress.
According to the Lemaitre strain equivalence hypothesis, the elastic modulus of rock after damage is as follows [22, 23]:where is modulus of elasticity before rock damage, is modulus of elasticity after rock damage, and is damage variable.
The stressstrain curve is consistent with the law of linear elasticity before material damage. The macroscopic constitutive relation considering cumulative damage of rock can be expressed as
The solution to is as follows.
According to the study of Wengui Cao, the yield judgment of rock element is appropriate with DruckerPrager criterion, and the strength intensity of rock elements conforms to Weibull distribution [24, 25]where is the parameter related to cohesion and internal friction angle ; is the spherical tensor invariant; is the deviator tensor invariant; , .
Assuming that the strength of rock element obeys the Weibull distribution, its probability density function is where and in (4) are Weibull distribution parameters:
According to (2), the constitutive relation of rock is as follows:where is confining pressure of rock and is the elasticity modulus of rock matrix. and can be obtained with value method of multivariate function. It uses the mathematical nature that the derivative of the highest point of the curve is 0:whereThe parameters , in (7) need to be gotten from fitting methods by experiment date and conform to the linear rule.
However, the damage defined in the traditional statistical model essentially reflects the damage caused by the yielding of the nondestructive matrix due to the external stress. Besides, the initial damage of the “fluid inclusions,” cracks, and holes existing in the rock before the external loading are not reflected in the model. Obviously, these are important factors influencing the initial elastic modulus of the rock. To establish the relationship between fluid inclusions and macroscopic constitutive relations of rocks, the initial factors of these rocks must be clearly described.
To solve this problem, it is necessary to establish the “rock element” which can reflect the constitutive relation of the rock part representatively, so fluid inclusionmatrix element (FME) is defined [26]. Theoretical basis of FME is RVE, which originates from rock element. Representative volume element (RVE) is an effective method to study such problems, which has been widely used in theoretical and numerical simulation of voids, fractured composites, and so on [27–31]. FME differs from rock element because a typical element of the same size contains a fluid inclusion, and the size is 3 to 5 times of the diameter around the inclusion. There are two types of pores in FME. One is connected pore, which indicates that the pore is connected with each other or connected with the cracks in the matrix. This kind of problem includes the pore which has fluid exchange with the outside, whether it has obvious boundary with the matrix or not. It covers the scope of fluid solid coupling, fluid replacement, and other engineering problems. Another kind is the closed pore, which has obvious boundary with matrix or no external exchange of fluid (gas) under certain confining pressure. This kind of FME contains rock pressure existing in the formation of geology, and it is an important carrier of the lockedin stress. Special statement is that pores mixed up with fluid inside FME are collectively referred to as “fluid inclusions” in this study.
3. Model
FME is mixed up with a fluid inclusion and rock matrix around a finite range, as shown in Figure 1 (0 is matrix; 1 is inclusion).
(a)
(b)
Figure 1 shows the difference between FME and rock element.
3.1. Fluid Inclusion Volume Modulus of FME
Pore fluid volume modulus of FME is not only related to porosity but also associated with bulk modulus of fluid and bulk modulus of skeleton . The relationship between and fluid inclusion compressibility is as follows:
can be expressed as [32, 33]where is the volume of FME and its rate of volume change is consistent with the overall rate of volume of the rock block . Therefore,
The compression factor of rock matrix can be defined aswhere is the stress acting on the rock skeleton. From (12) and (13),
From the point of the whole point of FME, the change of fluid inclusion volume is equal to the volume of fluid. So,and external stress on the rock has the relationship between fluid inclusion pressure and skeleton pressure ,where is the ratio of the rock block. Based on (16),
According to the theory of solid mechanics, the compression coefficient of the skeleton under the condition of elastic deformation can be expressed aswhere is Poisson’s ratio of the rock matrix. In general, is not obtained experimentally, while the overall modulus of matrix elasticity () is experimentally obtained [32]
In engineering applications, the contents of water, oil, and gas in the fluid inclusions are often different. When there are many kinds of fluid, the effect of fluid components on the fluid inclusion modulus should be considered. In the study of fluid replacement process, Murphy obtained the calculation method of fluid bulk modulus in fluid inclusion under different saturation [34]where is saturation. is a fluid volume modulus of fluid inclusion, and is the second fluid bulk modulus of fluid inclusion. In engineering applications, common reference for determination of fluid elastic parameters is presented in Table 1.

According to Table 1, different fluid or saturation has an important influence on the bulk modulus of elasticity in fluid inclusions. Based on the measured values of Table 1, the volume modulus of elasticity in fluid inclusions includes gas water condition 1.5~3.5 GPa, oil condition 0.9~2 GPa, and water condition 2~4 GPa.
Table 1 is an experiential value used in oil and gas engineering, which is convenient and quick to use. In most cases, there exists a situation that fluid inclusion pressure rather than fluid elastic modulus on rock stressstrain relation is the main consideration. Related researches demonstrate that the fluid inclusion pressure has a positive correlation with the pore compressibility, and the relation between the fluid inclusion volume modulus and the pressure can be depicted as follows [38]:where is the pressure increment relative to the original pore pressure. is the ratio of the initial volume to the actual volume of the fluid inclusions in the FME, and is the correction factor under different geological conditions or environment.
Study on compression coefficient or bulk modulus of FME is important for the equivalent elastic modulus based on Eshelby equivalent inclusion theory. It makes the fluid inclusion have practical meaning and physical meaning of certain form, which is different from the traditional solid mechanic problems or traditional inclusion damage model from this point. In addition, it is also a very important reason why this model is suitable to analyze lockedin stress.
3.2. Constitutive Relation of FME
The constitutive relationship of the FME containing different pore types is described by a uniform “equivalent elastic modulus,” which is based on the Eshelby equivalent inclusion theory and MT homogenization method. In 1957, Eshelby studied the law of stress around a phase in an infinite homogeneous matrix [39]. Since then, Eshelby equivalent inclusion theory and related analytical solutions have been developed into various forms with wide application in the study of composite and discontinuous medium mechanics [40–47].
As shown in Figure 2, according to the Eshelby inclusion theory, it is firstly assumed that the fluid inclusion in the FME is homogeneous with the matrix and free from external stress. If the inclusion has a strain due to some physical action (thermal expansion, etc.) in unconstrained condition (equivalent to the inclusion gotten from the matrix, as presented in Figure 1), it is called the intrinsic strain. The matrix will produce a bound strain on the inclusion due to the existence of intrinsic strain. Eshelby has proved that the relation between intrinsic strain and bound strain is as follows [41]:where is a constant tensor related to the shape of inclusions and the Poisson ratio of the matrix; it is called the Eshelby tensor and (25) can be expressed as
If it is a spherical inclusion in FME, can be expressed aswhere is matrix Poisson ratio and other tensor projects are 0.
In fact, the pore fluid inclusions in the FME are heterogeneous with the rock matrix and are subjected to external stress, which can be found in Figure 3. The strain produced by a heterogeneous inclusion with an external stress can be divided into the following three sections:(1)Compared with homogeneous material, the perturbation strain tensor produced by the substrate on the heterogeneous inclusions under the action of external stress is .(2)Confinement strain tensor of matrix to inclusion is .(3)Common strain tensor is , and it is produced by homogeneous inclusions and matrix due to external stress, conforming to the generalized Hooke’s law.
As is shown in Figure 3, the average stress tensor in inclusion is as follows:where is total stress tensor of inclusion; is perturbation stress tensor of matrix on inclusion; is the restraint stress tensor of matrix on inclusion; and is elastic modulus tensor of inclusion.
Average stress tensor of matrix in FME iswhere is total inclusion stress tensor; is elastic constants of matrix.
With Eshelby inclusion theory, (27) can be written asAccording to the MoriTanaka equivalent stress analysis method, the overall mean stress remains unchanged, which can be expressed as [46] where is volume ratio of pore in FME. It depends on 3~5 times of fluid inclusion diameter [26].
It can be obtained from (25), (27), and (30) thatwhere is unit tensor; from (25), (29), and (31), and its simple form iswhere ,
Under the action of an external stress tensor , an average strain tensor is generated, which can be obtained in accordance with Hooke’s law:where is equivalent elastic modulus tensor of FME.
3.3. Statistical Damage Model Based on Double MT Homogenization Theory
The stressstrain relationship of the “fluid inclusionsmatrix part” and the “matrix part containing no initial damage” has been obtained. The composition of the rock block is shown in Figure 4.
The stressstrain relationship between the two parts can be unified by the MT homogenization method [44]. It should be noted that the porosity of rock is and the ratio of the fluid inclusions to the FME mentioned above is . The method of double MT homogenization method is that the equivalent elastic modulus of FME is considered as an inclusion, and then the equivalent elastic modulus is solved in a wider range of rock matrix with FMEs. Therefore, the equivalent elastic modulus of rock block after double MT homogenization method iswhere is equivalent elastic modulus tensor of rock matrix based on double MT homogenization method.
According to (30), (31), (36), and (37), the statistical damage model of the initial fluid inclusions of the rock is expressed aswhere are stress tensors and is strain tensor.
In this way, the equivalent elastic modulus of FME and the constitutive equation of rock are combined, and the equivalent elastic modulus of FME is closely related to the fluid properties and pressure of fluid inclusions in FME. Then, the relationship between stressstrain curve of rock and fluid inclusions is established.
Numerous studies have demonstrated that rock porosity changes under different confining pressures should be considered. For example, in the study of rock pore seepage process, it was found that the pore porosity will change with the confining pressure. Additionally, the correlation study shows that the porosity changes with the different confining pressure change can be expressed as [48, 49]where is initial porosity of rock sample and reflects the whole process of rock pore changes. The value of is generally , which is consistent with the actual situation and can be determined according to the actual change of porosity in the experimental data.
It can be obtained from (32) thatIf the influence of fluid inclusion pressure of the initial FME on the partial elastic modulus is considered, where reflects the effect of pressure on tensor.
According to (41), the effect of porosity on the full stressstrain curve is discussed as follows.
(1) Confining pressure is 6.9 MPa, , = 90 GPa, and = 2 GPa, internal friction angle of rock remains unchanged, and changes with confining pressure.
Figure 5 shows the results. Compared with the situation that the effect of confining pressure on porosity is not considered, the effect of confining pressure on porosity can make the curve more faithfully reflect the strain softening characteristics of rock, and the slope of postpeak curve is steeper. Besides, an important feature of the model lies in the establishment of the relationship between the microscopic and macroscopic rocks. We will confirm whether the model’s macroscopic curve is consistent with the true experimental trend of stressstrain curves of rock. Combined with the experimental data of paper [50], the experimental and model data with different confining pressures are compared in later words.
(2) The initial porosity is not mentioned in the reference. When confining pressure is 6.9 MPa, the porosity = 0.35 remains unchanged and the initial porosity = 0.35 changes with confining pressure. Additionally, the porosity = 0.15 remains unchanged and the initial porosity = 0.15 changes with confining pressure, = 2 GPa, . Results are shown in Figure 6.
(3) When confining pressure is 13.8 MPa, the porosity = 0.35 remains unchanged and the initial porosity = 0.35 changes with confining pressure. The porosity = 0.15 remains unchanged and the initial porosity = 0.15 changes with confining pressure, = 2 GPa, . Results are shown in Figure 7.
Since the porosity in the original references was not considered or given, two sets of numerical values were selected for simulation. Through comparing Figures 6 and 7, when the confining pressure is 6.9 MPa, the data of = 0.25 is similar to the experimental data. When the confining pressure is 13.8 MPa, the data of = 0.25 is similar to the experimental data. Because initial porosities in initial references are not given for different samples, the initial porosity can be assumed for different confining pressures. The model of porosity changing with confining pressure can not only guarantee the accuracy of the variation law before the stress peak but also simulate the law of postpeak variation. It shows the rationality of considering fluid inclusion with confining pressure changing. Besides, there exists the possibility of explaining the influence of FME on the mechanical properties of rock.
Compared with the experimental data of the stressstrain curve, it is reasonable and effective in combining the Eshelby equivalent inclusion theory and the method of statistical mathematics to establish the model. The model has practicality and physical meaning. However, the figure also shows that there is a certain gap between the model curves beyond the yield limit and the development trend of the test curve and that the rock strain softening part is not well modeled. The theoretical establishment of this part is a direction that needs to be conducted in the future.
Meanwhile, on the basis of verification rationality, this model should not only verify the general rock compression experiment, but also validate the rock containing different types of fluid inclusions, especially for the constitutive model of rock mixed up with closed fluid inclusions or lockedin stress. Actually, experimental and verification work in this area is relatively few but significant.
4. Validation
Experimental verification includes two types of situations: one is the rock which contains fluidconnected pores and the other is to simulate the rock containing closed pores (lockedin stress) in rocklike materials.
4.1. Comparison with Experimental Results of Rock Mixed Up with Pore Water (Experimental Verification of the Connected Pore in FME)
The influence of the fluid inclusion on the total stressstrain curve of rock is rarely involved in the traditional statistical model of rock damage. In recent years, there are more studies about the influence of fluid inclusion on the overall characteristics of rock. Therefore, in this paper, the rationality of the model is validated by the threeaxis experiment data of rock containing pore water. If the influence of pore pressure variation on the total stress and strain of rock can be obtained, the applicability of the model with the connected pore in FME will be validated.
Through comparing the model with the sample data of rock in [51, 52], the unknown parameter is determined by the experimental data at the pore pressure of 0. At this time, the equivalent elastic modulus of the fluid inclusion is 0, and the independent variable is few. Based on the two sets of data, the parameters are recognized and corrected, and the unknown parameters of the original text are inversed. Rock parameters are as follows: = 14 GPa, , internal friction angle , initial porosity = 0.2, .
Fluid inclusion pressure of FME is 4 MPa and confining pressure is 20 MPa.
Fluid inclusion pressure of FME is 12 MPa and confining pressure is 20 MPa.
Figures 8 and 9 show that the model data is in good agreement with the experimental data of rock containing pore water. Before reaching the yield limit of the rock, the model basically agrees with the experimental data of rock containing pore water. Additionally, the model can simulate the law of full stress and strain and reflect the influence of pore water in FME on the total stress and strain of rock. That is to say, the relationship between the different pore water or the fluid state in the pore and the macroscopic characteristics of the rock can be analyzed using the model. Therefore, the model has practicality in analyzing rock mechanics problems with connected pore with an advantage in back analysis of pore properties through macroscopic phenomena.
4.2. Comparison with the Test Results of the LockedIn Stress in Rock (Experimental Validation of the Closed Pore in FME)
If the fluid state in FME can be reasonably reflected, the model not only can be used to analyze the rock block containing connected pores, but also has good applicability to the rock containing lockedin stress. However, the experimental researches about lockedin stress are few, and the experimental simulation itself is also confronted with a lot of difficulties. The change of equivalent fluid inclusion modulus of FME has mature application in fluid replacement in oil and gas engineering, while the experimental date about the constitutive relation of rock in this process is few. Consequently, in this part, the lockedin stress in rock is validated by simulating the closed pores with different equivalent modulus of elasticity in FME.
The experiment is designed that nitrile rubber inside rock matrix (rocklike material) can produce different expansion pressure at different temperature control conditions to simulate the different lockedin stress in rock. The equivalent elastic modulus of different inclusions in FME can be achieved from the above analysis. Besides, the matrix material consists of Portland cement and fine quartz sand in a certain proportion. The parameters are as follows: = 2.5 GPa, = 0.2 GPa, , internal friction angle , = 50 MPa, and . The experimental data of pore rubber materials were 2% and 4% respectively, and the lockedin stress was 0.8 MPa and 1.2 MPa, respectively. The transformation relationship between date of experimental situation and theoretical transformation value is presented in Table 2.

Figure 10 shows the result of the situation that porosity is 2% and equivalent elastic modulus of fluid inclusion is 0.5 GPa. Figure 11 shows the result of the situation that porosity is 4% and equivalent elastic modulus of fluid inclusion is 0.2 GPa.
The lockedin stress exerts significant influence on the experiment. The data in Figures 10 and 11 shows that the model curve reflects the change law before and after the peak in the actual experiment. Additionally, the random error control is within a certain range. However, it is undeniable that there is a certain difference between the model curve and the final part of the test date. Therefore, there are factors affecting the curve changes which have not been fully taken into consideration.
Therefore, from the experimental verification results of this part, it can be seen that the influence of lockedin stress in the deep rock mass on the macroscopic constitutive relationship should not be ignored. Highpressure and highenergy fluid inclusions in deep rock mass are the important reasons for such phenomena as rock burst and large deformation. Based on the model, we can analyze the fluid inclusions in the deep rock mass to predict or analyze some rock phenomenon. Additionally, the model provides a new idea for the solution of deep rock mechanics problems.
5. Conclusions
With the development of deep rock mechanics, more and more attention has been paid to the lockedin stress in rocks. Highpressure fluid inclusions are important parts of the lockedin stress in deep rock mass. Nevertheless, there are few studies on the mechanical properties of rock mixed up with closed pores. Based on RVE, this paper unifies the analysis methods of closed pore and connected pore, and FME is defined. Based on FME, the relationship between the microscopic fluid inclusions and the macroscopic mechanical properties of the rock is established using the statistical mathematical method.
The main works of this paper are as follows: on the basis of the traditional damage model, firstly, the constitutive relation of FME is built on Eshelby equivalent inclusions theory, and the stressstrain relation of FME is then determined. Secondly, the influence of FME on initial modulus of elasticity of rock mass is analyzed and the relationship between FME and the whole stress and strain of the matrix is established based on the Weibull distribution. The effect of fluid inclusions in FME on the stressstrain curve of rock is analyzed. Compared with the general rock compression experimental data, the model is proved to be reasonable. On this basis, the model is compared with the threeaxis confining pressure experimental date and lockedin stress simulation experiment of rock mass. Moreover, the model results are in good agreement with experimental date.
The main conclusions are as follows:
(1) FME is the medium of establishing the fluid inclusion and macroscopic stressstrain curve of rock. Combined with the principle of Eshelby equivalent inclusion, the FME can reflect the fluid state of the fluid inclusion and other factors with the concept of equivalent elastic modulus. FME is a basis of rock statistical damage model. It is an improvement of the traditional statistical damage model which cannot effectively reflect the initial fluid inclusion, crack, and other damage.
(2) The stress of the same strain changes with the increasing elastic modulus of the FME, but it does not influence the variation trend before and after the peak. However, the increase of porosity affects the gradient of the curve behind the peak, and the stressstrain curve becomes slow with the increase of porosity. The assumption that porosity varies with confining pressure can reflect the law of postpeak variation, which is greatly consistent with the experimental data.
(3) Compared with the different experimental data, the model based on FME and the mathematical statistics model is reasonable. Based on FME, it can reflect the influence of initial porosity or initial damage in a wider sense on rock mechanical properties, whether it is connected pore or closed pore in FME. As long as the initial damage can be reflected in FME, the influence of initial damage can be obtained from the model.
The innovations of the model are as follows:
(1) Using the Eshelby equivalent inclusion theory and MT uniform theory, the initial elastic modulus of rock is analyzed. Besides, it is an improvement of the traditional damage model, which cannot reflect the initial hole and damage of the rock, making the new statistical damage model more practical and reasonable.
(2) Based on FME, the relationship between rock microcosmic and macroscopic is established. The macroscopic constitutive relation of rock is analyzed by fluid inclusions in FME. Few studies have been conducted from this point of view.
(3) The model is employed to analyze the lockedin stress in deep rock mass, and the lockedin stress simulation experiment is carried out, which is groundbreaking for analyzing geological phenomena through the point of lockedin stress. Briefly, one important reason why we call it a new model is the effect of the lockedin stress on the rock. Therefore, the mathematical statistics model built on it is reasonable and effective.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors acknowledge the support of Jiangsu Natural Science Foundation of China (BK20141067).
References
 M. Vervoot, M. Wevers, R. Swennen et al., “Recent advances of Xray CT and its application for rock material,” in Proceedings of the Xray CT for Geomaterials; Soils, Concrete, Rocks, 2004. View at: Google Scholar
 L. Du, X. Zhang, J. Niu, X. Wang, and G. Feng, “The seismic CT method in measuring rock bodies,” Applied Geophysics, vol. 3, no. 3, pp. 192–195, 2006. View at: Publisher Site  Google Scholar
 D. Crandall, G. Bromhal, and D. H. Smith, “Conversion of a microCT scanned rock fracture into a useful model,” in Proceedings of the 2009 ASME Fluids Engineering Division Summer Conference, FEDSM2009, pp. 2101–2108, USA, August 2009. View at: Publisher Site  Google Scholar
 A. Sato and T. Mukunoki, “Fundamental Knowledge of Xray CT and Its Application to Rock Engineering,” Journal of MMIJ, vol. 131, no. 6, pp. 277–284, 2015. View at: Publisher Site  Google Scholar
 G. Molli, G. Cortecci, L. Vaselli et al., “Fault zone structure and fluidrock interaction of a high angle normal fault in Carrara marble (NW Tuscany, Italy),” Journal of Structural Geology, vol. 32, no. 9, pp. 1334–1348, 2010. View at: Publisher Site  Google Scholar
 B. E. Treeby and B. T. Cox, “kWave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields,” Journal of Biomedical Optics, vol. 15, no. 2, Article ID 021314, 2010. View at: Publisher Site  Google Scholar
 R. Klemd, “Highpressure fluidrock interaction as the cause of bulk rock compositional changes in blueschists and eclogites,” Egu General Assembly, vol. 14, no. 3980, 2012. View at: Google Scholar
 X. Zhuang, R. Huang, C. Liang, and T. Rabczuk, “A coupled thermohydromechanical model of jointed hard rock for compressed air energy storage,” Mathematical Problems in Engineering, vol. 2014, Article ID 179169, 11 pages, 2014. View at: Publisher Site  Google Scholar
 C. I. Steefel and K. Mäher, “Fluidrock interaction: A reactive transport approach,” Reviews in Mineralogy and Geochemistry, vol. 70, no. 1, pp. 485–532, 2009. View at: Publisher Site  Google Scholar
 F. Gassmann, “Elastic waves through a packing of spheres,” Geophysics, vol. 16, no. 4, pp. 673–685, 1951. View at: Publisher Site  Google Scholar
 B. H. Russell, K. Hedlin, F. J. Hilterman, and L. R. Lines, “Fluidproperty discrimination with AVO: A BiotGassmann perspective,” Geophysics, vol. 68, no. 1, pp. 29–39, 2003. View at: Publisher Site  Google Scholar
 J. M. Carcione, B. Gurevich, and F. Cavallini, “A generalized BiotGassmann model for the acoustic properties of shaley sandstones,” Geophysical Prospecting, vol. 48, no. 3, pp. 539–557, 2000. View at: Publisher Site  Google Scholar
 Z. Hou, C. Guo, and J. Wang, “Using Gassmann Equation Predict Marine Sediment Porosity,” Earth Science, 2016. View at: Google Scholar
 J.Y. Li, “Gas reservoir identification by seismic AVO attributes on fluid substitution,” Applied Geophysics, vol. 9, no. 2, pp. 139–148, 2012. View at: Publisher Site  Google Scholar
 N. Tisato, C. Madonna, and S. Boutareaud, “A new instrumentation to measure seismic waves attenuation,” in Proceedings of the AGU Fall Meeting, 2010, AGU Fall Meeting Abstracts. View at: Google Scholar
 T. K. Tan, “On the locked in stress, creep and dilataion of rocks and the constitutive equations,” Chinese Journal of Rock Mechanics Engineering, 1991. View at: Google Scholar
 L. Barrows, “Lockedin stress in viscoelastic Earth materials,” Eos, Transactions, American Geophysical Union, vol. 78, no. 24, p. 248, 1997. View at: Publisher Site  Google Scholar
 S. Wang, “Geological nature of rock and its deduction for rock mechanics,” Chinese Journal of Rock Mechanics and Engineering, vol. 28, no. 3, pp. 433–450, 2009. View at: Google Scholar
 Z. Q. Yue, “Gas inclusion and their expansion power as foundation of rock locked in stress hypothesis,” Journal of Engineering Geology, vol. 22, no. 4, pp. 739–756, 2014. View at: Google Scholar
 Z. Q. Yue, “On nature of earthquakes with cause of compressed methane gas expansion and migration in crustal rocks,” in Proceedings of the 5th Biot Conference on Poromechanics, BIOT 2013, pp. 507–516, Austria, July 2013. View at: Publisher Site  Google Scholar
 Z. Q. Yue, “On cause hypotheses of earthquakes with external tectonic plate and/or internal dense gas loadings,” Acta Mechanica, vol. 225, no. 45, pp. 1447–1469, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 J. Lemaitre, “A continuous damage mechanics model for ductile fracture,” Journal of Engineering Materials and Technology, vol. 107, no. 1, pp. 83–89, 1985. View at: Publisher Site  Google Scholar
 J. A. Lemaitre, “How to use damage mechanics,” Nuclear Engineering and Design, vol. 80, no. 2, pp. 233–245, 1984. View at: Publisher Site  Google Scholar
 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
 W.G. Cao and X. Li, “A new discussion on damage softening statistical constitutive model for rocks and method for determining its parameters,” Rock and Soil Mechanics, vol. 29, no. 11, pp. 2952–2956, 2008. View at: Google Scholar
 H. Xu, X. Chen, L. Dong, and Y. Yang, “Equivalent Elastic Constants of Rock Containing Pore Fluid Inclusions,” Journal of Basic Science and Engineering, vol. 25, no. 2, pp. 369–381, 2017. View at: Publisher Site  Google Scholar
 Z. F. Khisaeva and M. OstojaStarzewski, “On the size of RVE in finite elasticity of random composites,” Journal of Elasticity, vol. 85, no. 2, pp. 153–173, 2006. View at: Publisher Site  Google Scholar
 S. I. Ranganathan and M. OstojaStarzewski, “Scaling function, anisotropy and the size of RVE in elastic random polycrystals,” Journal of the Mechanics and Physics of Solids, vol. 56, no. 9, pp. 2773–2791, 2008. View at: Publisher Site  Google Scholar
 M. OstojaStarzewski, X. Du, Z. F. Khisaeva, and W. Li, “Comparisons of the size of the representative volume element in elastic, plastic, thermoelastic, and permeable random microstructures,” International Journal for Multiscale Computational Engineering, vol. 5, no. 2, pp. 73–82, 2007. View at: Publisher Site  Google Scholar
 A. Ramazani, K. Mukherjee, H. Quade, U. Prahl, and W. Bleck, “Correlation between 2D and 3D flow curve modelling of DP steels using a microstructurebased RVE approach,” Materials Science and Engineering: A Structural Materials: Properties, Microstructure and Processing, vol. 560, pp. 129–139, 2013. View at: Publisher Site  Google Scholar
 R. Pucha, “RVEbased design and analysis tools for composite materials with nanofillers,” Journal of Composite Materials, 2013. View at: Google Scholar
 M. B. Wold, L. D. Connell, and S. K. Choi, “The role of spatial variability in coal seam parameters on gas outburst behaviour during coal mining,” International Journal of Coal Geology, vol. 75, no. 1, pp. 1–14, 2008. View at: Publisher Site  Google Scholar
 I. L. Chuanliang and S. Zhu, “Compressibility coefficients of bulk volume and fluidsolid of rocks,” Lithologic Reservoirs, vol. 27, no. 2, pp. 1–5, 2015. View at: Google Scholar
 W. Murphy, A. Reischer, and H. Kai, “Modulus decomposition of compressional and shear velocities in sand bodies,” Geophysics, vol. 58, no. 2, pp. 227–239, 1993. View at: Publisher Site  Google Scholar
 Q. G. Liu and T. D. Tan, “Identifying lithological characteristics and the type of pore fluid in line of elastic moduli,” Geophysical Prospecting for Petroleum, vol. 32, no. 2, pp. 88–96, 1993. View at: Google Scholar
 X. W. Qiao and S. W. Yan, “Evaluation of the characteristics of pore fluids in permeable formations using acoustic logging date,” WLT, vol. 21, no. 3, pp. 61–66, 1997. View at: Google Scholar
 K. Huang, Z. Q. Xu, and H. X. Yang, “The relationship between P and S wave velocity ratio, modulus of elasticity and fluids contained in rocks,” XJPG, vol. 19, no. 5, pp. 20–22+85, 1998. View at: Google Scholar
 Q. Liu, R. Y. Lu, and X. Z. Li, “The Effect of the Pore Water Pressure on the Rock Mechanics Parameters,” Acta Geoscientia Sinica, vol. 29, no. 5, pp. 660–664, 2008. View at: Google Scholar
 J. D. Eshelby, “The determination of the elastic field of an ellipsoidal inclusion, and related problems,” Proceedings of the Royal Society London A: Mathematical, Physical and Engineering Sciences, vol. 241, pp. 376–396, 1957. View at: Publisher Site  Google Scholar  MathSciNet
 G. J. Rodin, “Eshelby's inclusion problem for polygons and polyhedra,” Journal of the Mechanics and Physics of Solids, vol. 44, no. 12, pp. 1977–1995, 1996. View at: Publisher Site  Google Scholar
 C. Q. Ru, “Eshelby's problem for twodimensional piezoelectric inclusions of arbitrary shape,” Proceedings of the Royal Society A Mathematical, Physical and Engineering Sciences, vol. 456, no. 1997, pp. 1051–1068, 2000. View at: Publisher Site  Google Scholar
 W. Zou, “Limitation of average eshelby tensor and its application in analysis of ellipse approximation,” Acta Mechanica Solida Sinica, vol. 24, no. 2, pp. 176–184, 2011. View at: Publisher Site  Google Scholar
 M. Kawashita and H. Nozaki, “Eshelby tensor of a polygonal inclusion and its special properties,” Journal of Elasticity: The Physical and Mathematical Science of Solids, vol. 64, no. 1, pp. 71–84, 2001. View at: Google Scholar  MathSciNet
 T. Mori and K. Tanaka, “Average stress in matrix and average elastic energy of materials with misfitting inclusions,” Acta Metallurgica et Materialia, vol. 21, no. 5, pp. 571–574, 1973. View at: Publisher Site  Google Scholar
 A. Giraud, Q. V. Huynh, D. Hoxha, and D. Kondo, “Application of results on Eshelby tensor to the determination of effective poroelastic properties of anisotropic rockslike composites,” International Journal of Solids and Structures, vol. 44, no. 1112, pp. 3756–3772, 2007. View at: Publisher Site  Google Scholar
 M. Hu, G.Y. Xu, and S.B. Hu, “Study of equivalent elastic modulus of sand gravel soil with Eshelby tensor and MoriTanaka equivalent method,” Rock and Soil Mechanics, vol. 34, no. 5, pp. 1437–1448, 2013. View at: Google Scholar
 H. Ma, C. Yan, and Q.H. Qin, “Eigenstrain boundary integral equations with local Eshelby matrix for stress analysis of ellipsoidal particles,” Mathematical Problems in Engineering, vol. 2014, Article ID 947205, 10 pages, 2014. View at: Publisher Site  Google Scholar
 Z. Chao, H. Wang, W. Xu, L. Yang, and K. Zhao, “Variation of permeability and porosity of sandstones with different degrees of saturation under stresses,” Chinese Journal of Rock Mechanics and Engineering, vol. 36, no. 3, pp. 665–680, 2017. View at: Publisher Site  Google Scholar
 H. Wang, W. Xu, J. Zuo, J. Shao, and C. Jia, “Evolution law of the permeability and porosity for lowpermeability rock based on gas permeability test,” Journal of Hydraulic Engineering, vol. 46, no. 2, pp. 208–216, 2015. View at: Publisher Site  Google Scholar
 H. W. Zhou, P. Wu, and R. Q. Yang, “Damage model of jointed rock mass,” Recent Advances in Rock Mechanics, 1989. View at: Google Scholar
 H. W. Yang, Study on Coupling Mechanism of Rock and Pore Water under Cyclic Loading, University of Chongqing, 2011.
 W. Wei, Z. Tian, and Q. Zhu, “Study of statistical damage constitutive model for rock considering pore water pressure,” Chinese Journal of Rock Mechanics Engineering, 2015. View at: Google Scholar
Copyright
Copyright © 2018 Xiao Chen 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.