Investigation of the relationship between dynamic and static deformation moduli of rocks

The determination of deformation parameters of rock material is an essential part of any design in rock mechanics. The goal of this paper is to show, that there is a relationship between static and dynamic modulus of elasticity (E), modulus of rigidity (G) and bulk modulus (K). For this purpose, different data on igneous, sedimentary and metamorphic rocks, all of which are widely used as construction materials, were collected and analyzed from literature. New linear and nonlinear relationships have been proposed and results confirmed a strong correlation between static and dynamic moduli of rock species. According to rock types, for igneous rocks, the best correlation between static and dynamic modulus of elasticity (E) were nonlinear logarithmic and power ones; for sedimentary rocks were linear and for metamorphic rocks were nonlinear logarithmic and power correlation. Moreover, with respect to different published linear correlations between static modulus of elasticity (Estat) and dynamic modulus of elasticity (Edyn), an interesting correlation for rock material constants was established. It was found that the static modulus of elasticity depends on the dynamic modulus only with one parameter formula.

methods to calculate the moduli: destructive and nondestructive procedures. In the destructive one, moduli are calculated from the stress-strain curves of the rock material. This is characteristic of the modulus of elasticity. For the non-destructive one, the most common method is an ultrasonic test, measuring both longitudinal and shear wave velocities.
Measuring the longitudinal and shear wave velocities, the dynamic elastic material parameters for isotropic and ideal elastic rocks are calculated with the help of the following formula : where E dyn = dynamic modulus of elasticity, GPa, G dyn = dynamic modulus of rigidity, GPa, K dyn = dynamic bulk modulus, GPa, J dynamic ¼ dynamic Poisso's ratio, q ¼ density; g cm 3 ,V S ¼ shearvelocity; km/s, V P ¼ longitudinal infinite medium velocity; km/s. The Young's modulus obtained from the compression (destructive) test is called the static elastic modulus (E stat ). The International Society for Rock Mechanics (ISRM) suggests three standard methods for its determination (Ulusay and Hudson 2007). They are as followings: • Tangent Young's modulus E tan -at fixed percentage of ultimate stress. This is defined as the slope of a line tangent to the stress-strain curve at a fixed percentage of the ultimate strength (Fig. 1a); • Average Young's modulus E av -of the straightline part of a curve. The elastic modulus is defined as the slope of the straight-line part of the stressstrain curve for the given test (Fig. 1b); • Secant Young's modulus E sec -at a fixed percentage of ultimate stress. It is defined as the slope of the line from the origin (usually point (0; 0)) to some fixed percentage of ultimate strength, usually 50% (Fig. 1c). In this paper, the secant static modulus has been calculated following ASTM D 3148-69 (1996).
The difference between dynamic (E dyn ) and static (E stat ) Young's Modulus for rocks has been addressed widely in rock engineering (van Heerden 1977;Lama and Vutukuri 1978;Barton 2006). Ratios for E dyn /E stat are typically in the range from 1 to 2 (Eissa and Kazi 1988).
Generally, the dynamic modulus of elasticity is slightly higher than the static value (Zhang 2006;Stacey et al. 1987;Al-Shayea 2004;Ide 1936;Kolesnikov 2009;Vanheerden 1987). The discrepancies between the dynamic and static elastic moduli have been widely attributed to microcracks and pores in the rocks. Figure 2 shows the ratio of dynamic elastic modulus to static elastic modulus compiled by Stacey et al. (1987). The ratio varies between about 1 and 3.
In other research ) conducted a laboratory experiment on ten different carbonate rocks quarried in Spain and received the following diagram (Fig. 3). As shown, for the majority of cases dynamic (E dyn ) Young's Modulus is higher than static (E stat ) Young's Modulus.
The rocks consist of homogeneous limestones, low anisotropic travertines, limestones and dolostones with abundant stylolites, veins and fissures. Ultrasonic waves were measured using non-polarised Panametric transducers (1 MHz), a precise ultrasonic device consisting of signal emitting-receiving equipment and an oscilloscope (TDS 3012BTektronix)   (Fig. 4).
Types of rocks are denoted below the picture and can be described as following. Blanco Alconera ( There are several relationships between the static and dynamic Young's moduli (Belikov et al. 1970;King 1983;Eissa and Kazi 1988;McCann and Entwisle 1992;Eissa and Kazi 1988;Christaras et al. 1994;Nur and Wang 1999;Brotons et al. 2014Brotons et al. , 2016Małkowskia et al. 2018). These equations can be divided into the following groups: • Linear • Non-linear: Logarithmic (power-law) and polynomial

Linear relationships
Up to now, several linear relationships have been established between the dynamic and the static Young's modulus. The following form was used:  (Stacey et al. 1987) where a and b are material parameters. These values are summarized in Table 1.

Non-linear relationships
Some published papers use a logarithmic relationship. Analyzing the measured data, they suggest the following formula between the dynamic and the static Young's modulus: The values of E stat , and E dyn are expressed in GPa while that of q bulk in g/cm 3 . Here c and d are material constants (the published values are summarized in Table 2).
These equations can be rewritten to the following form: where the parameters are presented in Table 3.
The power-law relation (7) was suggested in other researches, too (see Table 4). E.g. according to the data of Ohen (2003). the dynamic Young's modulus is about 18 times the static Young's modulus (Peng and Zhang 2007).
From ultrasonic test data of 600 core samples in the Gulf of Mexico, Lacy (1997) obtained the following polynomial correlation for sandstones (Peng and Zhang 2007): A similar connection exists for shales and sedimentary rocks (Horsrud 2001; see Table 5). For this study, 40 samples of different types of rocks from the literature were analyzed (Lama and Vutukuri 1978). Data were classified according to rock types and the relationship between dynamic and static constants was investigated in Table 6. Also, the histogram of investigated parameters is shown in Fig. 5.
According to Fig. 6, there is a linear regression between static and dynamic modulus of elasticity for all studied rocks, giving R 2 = 0.89. More precisely, Fig. 7 shows the relationship between static and dynamic modulus of elasticity based on rock types. As it is shown, for igneous rocks there is a linear correlation that gives the value of R 2 = 0.95, for sedimentary rocks the value is R 2 = 0.90 and for        Fig. 8 illustrates the correlation between static and dynamic modulus of rigidity for all samples, giving R 2 = 0.89. Figure 9 demonstrates the relationship between static and dynamic modulus of rigidity based on rock types. As it is clear, for igneous rocks the value of R 2 = 0.96, for sedimentary rocks the value of R 2 = 0.91 and for metamorphic rocks the value of R 2 = 0.63. Figure 10 exhibits the linear correlation between static and dynamic bulk modulus for all studied rocks, giving R 2 = 0.77. Figure 11 shows the relationship between static and dynamic bulk modulus with respect to rock types. As can be seen, for igneous rock the value of R 2 = 0.77, for sedimentary rocks the value of R 2 = 0.38 and for metamorphic rocks, the value of R 2 = 0.87. The results are summarised in Tables 7, 8, 9.
Our achieved results in Fig. 12, with previously published correlations, are illustrated in Fig. 13. It can be concluded that the formulas well fit the data.
According to our analysis, with respect to different published linear correlations between static modulus of elasticity (E stat ) and dynamic modulus of elasticity (E dyn ), a relationship with high correlation (R 2 = 0.91) was observed between a and b parameters as it is shown in Fig. 14. It should be stated that data published by McCann & Entwisle (1992) for Crystalline rocks were not included in this analysis. The reason is that by applying their equation, the amount of correlation decreases from (R 2 = 0.91) to (R 2 = 0.57). It might be related to the link between the crystallized structure of rock and wave propagation.

Results and discussion
In this research, the basic geomechanical properties of different types of rocks were measured and analyzed. The results show that there is a good correlation between static and dynamic elasticity modulus, rigidity modulus and bulk modulus. Regarding the relationships between static and dynamic modulus of elasticity, the best correlation found to be nonlinear logarithmic and power regression with the value of (R 2 = 0.91). Similarly, Brotons et al. (2014Brotons et al. ( , 2016 established the nonlinear correlation between static and dynamic elastic modulus of different types of rocks with the value of R 2 = 0.99. Eissa and Kazi (1988) carried out similar research and discovered that the best correlation was to be nonlinear with the value of R 2 = 0.92. On the other hand, some other researchers adopted the linear correlation as wellfitted regression line (Belikov et al. 1970) for Granite with R 2 = 0.92, (King 1983) for Igneous and metamorphic rocks with R 2 = 0.82, (McCann and Entwisle 1992) for crystalline rocks with R 2 = 0.82, (Christaras et al. 1994) for all types of rocks with R 2 = 0.99, (Nur and Wang 1999) for all types of rocks with R 2 = 0.8).
Nevertheless, in the present study, we established the linear correlations for igneous rocks with R 2 = 0.95, for sedimentary rocks with R 2 = 0.90 and for metamorphic rocks with R 2 = 0.69. Considering the relationship between static and dynamic modulus of rigidity, the best correlation observed was nonlinear logarithmic regression with giving the value of (R 2 = 0.97); when it comes to bulk modulus, the best correlation was linear with the value of (R 2 = 0.77). According to rock types, for igneous rock, the best correlation between static and dynamic modulus of elasticity (E) was nonlinear logarithmic and power ones with the value of (R 2 = 0.96); However, King (1983) found that the best correlation was linear for igneous-metamorphic rocks with the value of (R 2 = 0.82). For sedimentary rocks, the best correlation was linear with the value of (R 2 = 0.88) and for metamorphic rocks was nonlinear logarithmic and power with the value of (R 2 = 0.93). For igneous rocks, the best correlation between static and dynamic modulus of rigidity (G) was nonlinear logarithmic with the value of (R 2 = 0.96); for sedimentary rocks was nonlinear logarithmic with the value of (R 2 = 0.98); for metamorphic rocks was nonlinear logarithmic with the value of (R 2 = 0.99). For igneous rocks, the best correlation between static and dynamic bulk modulus (K) was nonlinear logarithmic with the value of (R 2 = 0.88); for sedimentary rocks was linear with the value of (R 2 = 0.38); for metamorphic rocks was nonlinear logarithmic with (R 2 = 0.98). These values demonstrate an interesting finding that there is a higher correlation between static and dynamic constants in igneous rocks rather than sedimentary and metamorphic rocks except for bulk modulus. Also, based on an analysis of previously obtained linear correlations between static and dynamic modulus of elasticity in Table 1, there is a good logarithmic correlation between constant parameters (a, b). Interestingly enough, our achieved result as shown in Fig. 14 fits well with previously published results with high correlation R 2 = 0.91. It means the static modulus of elasticity depends on the dynamic modulus only with a one-parameter formula: where b is rock type -dependent parameter: Additionally, for a more accurate comparison, rootmean-square (v) errors between the dynamic and static elastic modulus was calculated as: where P obs j ð Þ is the observed value of a parameter in the jth sample, here is the modulus of elasticity (E), P cal j ð Þ is the calculated value of a parameter in the jth sample, j = 1, 2,…,n, is the number of tested samples. Based on our analyses we received the following correlation, Fig. 6 Relationship between static and dynamic modulus of elasticity for all rocks Fig. 7 Relationship between static and dynamic modulus of elasticity based on rock types b Fig. 5 Histogram of investigated parameters log 10 E stat ¼ 1:37 log 10 q bulk E dyn À Á À 1:32; R 2 ¼ 0:91; v ¼ 10:52 The difference between dynamic and static elastic moduli is highly associated with mineralogical differences, differences in grain/crystal size, differences in porosity. A deviation of dynamic elastic modulus from static elastic modulus can be attributed to the presence of fractures, cracks, cavities and planes of weakness and foliation (Al-Shayea 2004;Guéguen and Palciauskas 1994). In other words, as the number of discontinuities increase, the lower value of Young's modulus and the higher discrepancy between static and dynamic values are expected. The most crucial  petrographic parameter that influences the deformation behavior of rock is porosity. The trend observed in porous rocks was that the static elastic modulus was inversely proportional to porosity (Brotons et al. 2016;Garcia-del-Cura et al. 2012). However, Crystalline rocks exhibit lower values of elastic moduli due to the fact that their porous system is constituted by a dense microcrack network. That is, crystalline rocks act as non-continuous solid, while porous rocks with interparticle porosity behave as a more continuous solid due to the presence of cement and matrix between grains. The effect of mineralogy on elastic moduli of rocks have been studied by several authors and found to be much less than other factors such as porosity and crystal size (Heap and Faulkner 2008;Palchik and Hatzor 2002). Regarding the effect of grain size on the elastic modulus of nanocrystalline, with the decrease of grain size, the elastic modulus decreases (Kim and Bush 1999;Chaim 2004;Zhang and Tahmasebi 2019;Tugrul and Zarif 1999).
It is worth mentioning that more detailed material models beyond ideal elasticity give an exact relationship between the elastic and static moduli. Notably, the observed relations can be explained in a universal thermodynamic framework where internal variables  (Asszonyi et al. 2015;Berezovski and Ván 2017). The difference of dynamic and elastic moduli is natural in this framework as it is clear from the corresponding dispersion relation and related laboratory experiments (Barnaföldi et al. 2017;Ván et al. 2019) These constitutive models are based only on universal principles of thermodynamics, are independent of particular mechanisms and are successful in characterizing rheological phenomena in rocks including and beyond simple creep and relaxation. This is in accordance with the difficulty for finding a very detailed quantitative mesoscopic mechanism for the dynamics of dissipative phenomena in rocks as well.

Conclusion
Generally, the performance difference between the linear and power-law formulas was small, the linear relationship fitted well the data. Therefore, considering also the possibility of one parameter formulation, given in (11), the linear relationship between dynamic and static elastic moduli is suitable for rock engineering. It was expected, that the modulus of rigidity and bulk modulus, as original theoretical Lamé parameters, would correlate better. However, the correlation between the static and dynamic Young modulus (R 2 = 0.89) is as good as for the other ones (R 2 = 0.89 and 0.77), in spite of the fact that the Young's modulus is a composite parameter. The reason could be the uncertainty in the Poisson ration measurements and also the difference in rock types for the bulk modulus measurements is remarkable. It is worth mentioning that thermodynamic principles more detailed material models beyond ideal elasticity give detailed relationships between the elastic and static moduli in general. Particularly, the observed relations can be explained in a universal thermodynamic framework where internal variables are characterizing the structural changes in the rock.