# Calculation of Elastic Modulus for Fractured Rock Mass Using Dimensional Analysis Coupled with Numerical Simulation

Authors: Fengyu Ren; Jing Zhang; Zhihua Ouyang; Hao Hu
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

Underground mining activities make the fractures in the natural rock mass develop randomly. The elastic modulus of the fractured rock massEm is changed with the redistribution and development of the fractures. An equivalent model of fractured rock mass is structured to represent the hydraulic conductivity and the rock mass strain because of the continuum theory. Dimensional analysis is very useful to build relationship of the parameters in complex physical phenomena. Based on the engineering phenomenon of groundwater flowing into the goaf along the fracture in the rock mass, a fuzzy expression among parameters such as the parameter Em, the seepage flow Q, and the exposed area of the goaf A is obtained using dimensionless analysis. To calculate the parameter Em, the fuzzy relationship is then characterized by Darcy’s law and numerical simulation. Under the scripting environment of Python, an automated program to realize the numerical simulation of all scenarios is established, which also provides convenience for drawing the dimensionless flux charts. The results show that the parameter Em can be calculated by the dimensional analysis coupled with numerical simulation. In addition, the parameter Em decreases with the increase of the parameter Q, and the integrity of rock mass is also worse. Finally, a mine example is used to verify the feasibility of the method.


---

## Body

## 1. Introduction

The elastic modulus of rock mass is one of the most important mechanical parameters in geotechnical engineering and has important guiding significance for engineering practice [1]. Considering the cutting effect of the structural plane in the rock mass, numerous fractures are randomly formed inside the rock mass. Therefore, the study of the elastic modulus of fractured rock mass is more suitable for engineering practice than the complete rock mass. Most scholars carried out many works to obtain more practical elastic modulus, such as empirical estimation method [2–8], acoustic wave test method [9, 10], and in situ test method [11]. However, the empirical estimation method is subject to subjective factors, which leads to the estimation results vary with different engineers. The acoustic wave test method and in situ test method are costly and are obviously interfered with by environmental factors. However, it is worth noting that the above methods indicate that the fractures in the rock mass have a great influence on the elastic modulus, so we attempt to study the fractures of the rock mass to calculate the elastic modulus.The fractures in natural rock mass are chaotic and disorderly [12]. Most scholars [13–16] research the mechanical behavior of fractured rock mass by establishing an equivalent continuous medium model, that is, to convert discrete fractured rock mass into continuous medium by equivalent hydraulic conductivity methods. The equivalent continuum model simplifies the study of the mechanical behavior of fractured rock mass and provides convenience for studying the elastic modulus of fractured rock. On the other hand, the natural rock mass is generally saturated with groundwater, which produces seepage within the rock mass [17, 18]. The seepage characteristics are determined by the fractures. Considering that the seepage flow in the exposed area of the goaf can be measured, we can establish the connection between the elastic modulus and the seepage flow in the goaf by utilizing the fractures of rock mass.In engineering, the elastic modulus and fractures of the rock mass change due to the influence of mining, and this change involves many physical quantities, not only related to the physical and mechanical parameters of rock mass and the mining parameters but also to the initial state of the fractures. To find the relationship between the elastic modulus and the seepage flow, we attempt to utilize dimensional analysis coupled with numerical simulation, say, finite element method (FEM). Dimensional analysis [19–21] is widely used to investigate hypotheses about complex physical phenomena. The Buckingham Pi theorem [22–25] in dimensional analysis provides a way to form a function of dimensionless variables. As a result, the stress and strain at any circumstances can be easily estimated using FEM following the pattern given in the dimensionless function.In the following sections, based on the engineering phenomenon of groundwater flow into the goaf, we regard the fractured rock mass as an equivalent continuous medium model and derive the expressions of hydraulic conductivity and rock mass strain under three-dimensional conditions. Then, we integrate the physical quantities of this physical phenomenon using the Buckingham Pi theorem and establish the functional relation between the elastic modulus and the seepage flow. Considering that the functional relationship contains the rock mass strain, we use the numerical simulation to solve the strain value in the rock mass and the indirect coupling method to calculate the hydraulic conductivity and the corresponding seepage flow. With the help of the scripting language Python, the FEM input files for all scenarios formed automatically can be executed. Consequently, we can draw the dimensionless flux charts about the elastic modulus conveniently; that is, we can obtain the elastic modulus of fractured rock mass under any working condition. Finally, we take a mine example to verify the accuracy of this method.


## 2. Equivalent Model for Fractured Rock Mass

The equivalent model of fractured rock mass (as shown in Figure1) approximately treats the rock mass object as a continuous medium and assumes that the normal direction of the fractured rock group corresponds to the direction of coordinate axis [26]. Generally, the normal directions of the surfaces in the above three groups of fractures are arbitrary and not perpendicular to each other because of the anisotropy of rock mass. However, Warren and Root [27] proposed an equivalent model for fractured rock mass that transformed the complex fracture system to a three-dimensional state with three groups of mutually orthogonal fracture systems, which greatly simplifies the study of the hydraulic conductivity. In this way, we can analyze the mechanical behavior of fractured rock mass by using the continuum mechanics method and the seepage theory from a macro perspective, regardless of the geometry structure of a particular fracture.Figure 1
The equivalent model for fractured rock mass, where the parametersB and S are the parallel fractures of aperture and spacing, respectively.Regardless of the deformation of the material, the hydraulic conductivity of a set of parallel fractures can be defined as [28, 29](1)K=gB312μkS,where g is the gravitational acceleration and μk is the kinematic viscosity of water. The parameters B and S are the parallel fractures of aperture and spacing, respectively. Equation (1) expresses the hydraulic conductivity in the absence of deformation and for a prescribed initial fracture aperture.Under the condition of orthotropic, the hydraulic conductivity of three groups of basic fractures that the normal direction is consistent with the direction of coordinate vector can be expressed as(2)Kij=Kij1+Kij2+Kij3,i,j=1,2,3,where Kij is the tensor for hydraulic conductivity and the numbers 1, 2, and 3 refer to the coordinates x, y, and z in Figure 1, respectively. Kij1 can be expanded as(3)Kij1=K111K121K131K211K221K231K311K321K331.Obviously, under the condition of orthotropic, wheni ≠ j and Kij1=0, and when i = j = 1 and Kij1=0, then(4)Kij1=0000K221000K331.And equation (2) can be expressed as(5)Kij=K112+K113000K221+K223000K331+K332.Assuming that the aperture of fractures in the three groups of fractures is expressed byB1, B2, and B3, respectively, and the spacing of fractures in the three groups of fractures is expressed by S1, S2, and S3, respectively. Then,(6)K11=K112+K113=gB2312μkS2+gB3312μkS3,K22=K221+K333=gB1312μkS1+gB3312μkS3,K33=K331+K332=gB1312μkS1+gB2312μkS2.According to equation (6), the hydraulic conductivity is positively proportional to the cube of the aperture of fractures, indicating that the aperture of fractures plays a leading role in the change of the hydraulic conductivity. Regardless of the influence of the change of the spacing of fractures, we label the variables of the aperture of fractures as Δub1, Δub2, and Δub3. Assuming that the initial aperture and spacing of the three groups of fractures in the fractured rock mass are the same (namely, B1 = B2 = B3 = B; S1 = S2 = S3 = S), the fractured rock mass equivalent hydraulic conductivity can be expressed as(7)K11=g12μkSB+Δub23+B+Δub33,K22=g12μkSB+Δub13+B+Δub33,K33=g12μkSB+Δub13+B+Δub23.As shown in Figure2, we select the model in Figure 1 to project onto the x-z plane. Taking a set of fractures in the rock mass as the analysis object, we regard the fractured rock mass as an equivalent continuous rock mass.Figure 2
The equivalent continuous model atx-z plane.From Hooke’s law, for a linear rock matrix,(8)Δεs=ΔσE,where Δεs is the strain in the direction of the applied stress Δσ and E is the elastic modulus of the rock matrix.Therefore, between two fractures in the rock mass at separationS, the displacement Δus is as follows:(9)Δus=SΔεs=SΔσE.Defining the deformation across a single fracture asΔub, then [26](10)Δub=Δσkn,where kn is the normal fracture stiffness.Total displacementΔum is obtained from the sum of component displacements as(11)Δum=Δus+Δub=S+BΔσEm,where Em is the elastic modulus of the fractured rock mass.Substituting equations (9) and (10) into equation (11) gives(12)Em=EknS+BSkn+E.Rearranging equation (11) yields(13)Rm=EmE=S+BknSkn+E,where Rm is the elastic modulus reduction ratio [26].Rearranging equations (10) and (9),(14)Δus=SknΔubE.Substituting equation (14) into equation (11) gives(15)Δub=EΔumSkn+E=ES+BΔεSkn+E,where Δε is the strain for equivalent continuous rock mass and Δε=Δσ/Em.Substituting equation (13) into equation (15) and rearranging the equation yield(16)Δub=1−RmS+BΔε.Generalizing equation (16) to the three-dimensional situation,(17)Δubi=1−RmS+BΔεi,i=1,2,3.Substituting equation (17) into equation (7) gives(18)K11=g12μkSB+1−RmS+BΔε23+B+1−RmS+BΔε33,K22=g12μkSB+1−RmS+BΔε13+B+1−RmS+BΔε33,K33=g12μkSB+1−RmS+BΔε13+B+1−RmS+BΔε23.Equation (18) expresses the functional relationship between the hydraulic conductivity and the strain. Meanwhile, the relationship shows that the strain has a great influence on the hydraulic conductivity. When the rock mass is excavated, the groundwater will flow into the excavation area. Darcy’s law [30] suggested that the relationship between the seepage flow and the hydraulic conductivity in the anisotropic medium can be expressed as(19)Qi=−KijJ,where Qi is the seepage flow. J is the hydraulic gradient.As to the orthogonally anisotropic medium, expanding the equation (19) can obtain(20)Q1Q2Q3T=K11000K22000K33J1J2J3.Therefore, as to the mined excavations, the total seepage flow of the goaf can be expressed as(21)Q=Q12+Q22+Q32.


## 3. Buckingham Pi Theorem

The Buckingham Pi theorem can describe physical processes involvingn physical quantities (x1, x2,..., xn), and the physical processes can be expressed as functional relations of the physical quantities involved as follows:(22)fx1,x2,…,xn=0.Among then physical quantities involved in the physical process, the n−m nonbasic physical quantities can be represented by the composite quantity of m basic physical quantities and converted into n−m dimensionless parameters, then equation (22) can be expressed as [22, 31](23)fπ1,π2,…,πn−m=0.The important fact to notice is that the new relation involves fewer variables than the original relation, which simplifies the theoretical analysis and experimental design alike.The development of fractures in rock mass determines the initial hydraulic conductivity and the elastic modulus of the fractured rock mass. Considering that the stress state of the surrounding rock disturbed by the mining activity, a certain displacement can be generated in the normal direction of the fracture surface, which caused the hydraulic conductivity to change. In addition, a large number of goafs of different sizes are formed by the excavation of the mine, and the exposed area of the goaf is directly related to the seepage flow, indicating that there is a fuzzy relationship between seepage flow and elastic modulus of the fractured rock mass. However, the fuzzy relationship involves many parameters (as shown in Figure3), such as the exposed area of goaf, the depth of goaf, the bulk weight of rock mass, and the aperture and spacing of fractures.Figure 3
Geometry diagram of physical parameters, wherel, w, and h are the goaf span, width, and height, respectively; H is the depth of goaf.Considering that the fuzzy relationship between the elastic modulus of the fractured rock mass and the seepage flow involves many parameters, and so far there is no clear equation to express the relationship, we utilize the Buckingham Pi theorem to derive the relationship as follows:(24)Q=fB,S,H,A,h,γ,g,Em,where Q is the seepage flow, H is the depth of the goaf, A is the exposed area of the goaf, h is the height of the goaf, γ is the bulk weight of rock mass, and g is the gravitational acceleration. Table 1 summarizes all physical quantities and corresponding dimensions in the fuzzy relationship.Table 1
Units and dimensions of physical quantities.
Physical quantitiesSymbolUnitDimensionsSeepage flow in the goafQm3/sL3T-1Parallel fractures of apertureBmLParallel fractures of spacingSmLDepth of the goafHmLExposed area of the goafAm2L2Height of the goafhmLBulk weight of rock massγkg/(m2s2)ML-2T-2Gravitational accelerationgm/s2LT-2Elastic modulus of the fractured rock massEmkg/(ms2)ML-1T-2The physical quantitiesQ, h, and Em are selected as the basic physical quantities, and the remaining physical quantities are regarded as the quantities to be determined. Then,(25)π1=Qα1hβ1Emγ1Sη1,π2=Qα2hβ2Emγ2Bη2,π3=Qα3hβ3Emγ3Aη3,π4=Qα4hβ4Emγ4γη4,π5=Qα5hβ5Emγ5gη5,π6=Qα6hβ6Emγ6Hη6,where αi, βi, γi, and ηi are the parameters to be determined, and i = 1, 2, 3, 4, 5, and 6.Taking equation (25) as an example, we can obtain(26)MLT0=L3T−1α1Lβ1ML−1T−2γ1Lη1.From equation (26), we can obtain αi = γi = 0 and ηi = −βi; then,(27)π1=hSβ1.Similarly,π2, π3, π4, π5, and π6 can be expressed as(28)π2=hBβ2,π3=h2Aβ3,π4=Emγhγ4,π5=h5gQ2η5,π6=hHβ6.Therefore, the dimensionless relationship can be expressed as(29)gh5Q2η5=fhSβ1,hBβ2,h2Aβ3,Emγhγ4,hHβ6.The values ofη5, β1, β2, β3, γ4, and β6 in equation (29) can change the axis values of the dimensionless flux charts but cannot change the final evaluation results. To make the calculation more convenient, the values of the parameters need to be clear. Considering that the parameter Q is one of the key parameters for calculating the elastic modulus for fractured rock mass, we can take η5 = 1/2. Referring to the expression of the hydraulic conductivity equation (1) and the dimension of B and S, the expression B/S can be a dimensionless expression. Besides, the magnitude of the parameters h and B is quite different; therefore, we can take β1 = 1 and β2 = −1 to obtain the dimensionless expression B/S. Generally, the value of parameter A is greater than the parameter h2; we take β3 = −1 to make the dimensionless flux chart clearer. Given that the parameter h can be eliminated by using the dimensionless expression Em/γh to multiply the dimensionless expression h/H, we take γ4 = 1 and β6 = 1. Finally, we can obtain(30)gh5Q=fhS,Bh,Ah2,Emγh,hH.Considering that the multiplication and division of dimensionless quantities have no effect on the final analysis results, equation (30) can be expressed as(31)Qgh5=fBS,Ah2,EmγH.Poisson’s ratioν can reflect the lateral deformation of the material. As one of the important factors in the calculation of elastic mechanics, it has a certain effect on the seepage flow and elastic modulus. Considering that its dimension is 1, equation (31) can be expressed as(32)Qgh5=fBS,Ah2,EmγH,ν,or(33)Qgh5=fBS,Ah2,EmνγH.


## 4. Numerical Simulation

Zhengxing Iron Mine located in Daye city, Hubei province, China, is taken in this work as the engineering background of numerical analysis. The surrounding rock mass in the mine is relatively broken and mainly consists of skarn, which is a medium-basically stable rock mass. The mining method is the sublevel open stope mining method. The mine uses top-down sequential mining, and the mined levels are −120 m, −160 m, −200 m, −240 m, −280 m, and −320 m. Multiple goafs of varying sizes are formed due to the unreasonable mining planning. The maximum exposed area of the goaf is 10000 m2, and the minimum exposed area is 2400 m2. The maximum height of the goaf is 50 m, and the minimum is 10 m. There are many developed local fault zone and derived fractures in the mine because of the unreasonable support. Considering that the in situ stress gradually increases with the mining depth, the confining pressure of the goaf is more and more complicated. In addition, the increase of mining depth accelerates the development of fractures in rock mass.To realize the automatic execution for numerical simulation of all scenarios, we use Python as the scripting environment for modeling and finite element calculation. Gmsh [32–34], a 3D modeling software, can conduct finite element mesh division. In addition, the output of this software can be .tmpl file format for Python software to call, so that we can use Python’s loop function to change the model parameters in the .tmpl file, so as to achieve automatic model establishment and mesh division. SfePy [35, 36] (Simple Finite Element in Python) is an open-source finite element analysis software written in Python, which uses the Galerkin method to solve. At present, the SfePy has mature models in various aspects such as acoustics, fluid mechanics, and elastic mechanics. Use Python to call the SfePy module and load the model created by Gmsh into SfePy for setting boundary condition and finite element calculation. Finally, the dimensionless diagram of elastic modulus can be obtained by using the scientific calculation module of Python to substitute the finite element calculation results into the function relations determined by the Buckingham Pi theorem. The flowchart of numerical simulation is shown in Figure 4.Figure 4
Flowchart of numerical simulation.Generally, numerical simulation simplifies the actual size of irregular goaf for finite element calculation. The modeling and meshing based on Gmsh are shown in Figure5. The size of the rock mass is 1000 m∗500 m∗1100 m, and the size of the goaf is an unfixed value. We use the parameters of l, w, and h in Figure 3 to control the size of the goaf. The density of the mesh is increased near the goaf.Figure 5
Modeling and meshing based on Gmsh.When the model is loaded into the SfePy, considering that the finite element calculation in this work involves the coupling calculation of elasticity and fluid mechanics, we set two boundary conditions: (1) stress boundary condition: the vertical stress is calculated according to the gravity stress, and the horizontal stress is calculated according to the Poisson effect. The lateral side of the rock mass restricts the horizontal displacement, and the bottom of the model is a fixed constraint. The Earth surface and the goaf are free. (2) Stress boundary condition: the model is a saturated rock mass medium, and a fixed water head value is applied to the boundary around the rock mass. The water head value is the distance from the underground water surface to the goaf. The Earth surface is a free boundary and the pore water pressure is zero.The physical and mechanical parameters of Zhengxing Iron Mine are shown in Table2. The parameters of γ, g, υ, and μk are fixed values, which are taken from field investigation and laboratory experiment. Then, as shown in equation (39), each dimensionless expression contains two variable parameters, for example, the expressionEm/νγH, where Em and H are the variables. However, the initial value of parameter Em for finite element calculation can take the experience value of 50 GPa for Zhengxing Iron Mine (as the parameter E shown in Table 2). The value ranges of the parameters B and h can be obtained by the mining status of Zhengwang Iron Mine. The value ranges of the parameters B and S measured based on the field survey are 1 × 10−4∼5 × 10−3 and 0.05∼10.0, respectively. To keep the value range of the dimensionless expression B/S unchanged, the value of S can be fixed as 1.0 m; then the value range of B is 1 × 10−5∼1 × 10−1 m. The parameter Rm is related to the parameters B and S [26], where we can set the reference value of the parameter to 0.8.Table 2
The physical and mechanical parameters of Zhengxing Iron Mine.
Physical quantitiesInitial valueFinal valueParallel fractures of apertureB (m)1 × 10−4∼5 × 10−31 × 10−5∼1 × 10−1Depth of the goafH (m)120∼3201 × 102∼1 × 103Span of the goafl (m)—3.0∼160.0Exposed area of the goafA (m2)2.4 × 103∼1 × 1042.4 × 103∼1 × 104Height of the goafh (m)10.0∼50.010.0∼50.0Parallel fractures of spacingS (m)0.05∼10.01.0Elastic modulus of the rockE (GPa)50.050.0Width of the goafw (m)20.020.0Bulk weight of rock massγ (kN/m3)27.027.0Gravitational accelerationg (m/s2)9.89.8Poisson's ratioυ0.250.25Kinematic viscosity of water (20°C)μk (m2/s)1.01 × 10−61.01 × 10−6Elastic modulus reduction ratioRm—0.8The value ranges of the parametersH and A in Zhengxing Iron Mine are shown in Table 2. The parameter A can be expressed as(34)A=2lw+lh+wh.Once the parametersw and h are determined, the value of A is determined by the span of the goaf l. Therefore, the value range of l is 3.0∼160.0 m. The final value of all physical parameters is shown in Table 2.When the seepage flow into the goaf is calculated, the strain induced by the mining process is first determined by the finite element model. Then, the seepage flowQ can be obtained by equation (21). The results of numerical simulation are then substituted into equation (33) to obtain the dimensionless parameters. As the size of the goaf changes, we must process a large amount of data. Consequently, we establish an automated program to handle the onerous task. This procedure is described in the flowchart of Figure 6, where the parameters B, H, and l are variables, and the parameter B controls the development of fractures in different rock masses, the parameter H controls the gravity stress, and the parameter l controls the exposed area of the goaf.Figure 6
Flowchart of finite element model.The calculation of the seepage flow is combined with theoretical analysis and numerical simulation. The mathematical model of saturated continuous medium is established by Gmsh, and the finite element numerical calculation results are obtained by SfePy. Using the value of element strain in the numerical simulation results, the hydraulic conductivity can be obtained by equation (18). Combining the hydraulic conductivity and the value of pore water pressure in the numerical simulation results, the seepage flow can be calculated by equation (21). The hydraulic gradient J can be obtained by Darcy’s law [30]:(35)J=ΔHwL,where ΔHw is the water head drop of nodal element sections and L is the spacing of nodal element sections.The water headHw can be expressed as(36)Hw=pγw+hw,where p is the pore water pressure, γw is the bulk density of water, and hw is the position head.


## 5. Results

Based on the dimensional analysis coupled with numerical simulation, the dimensionless flux charts of the fractured rock mass parameters obtained using the flow chart shown in Figure6 are shown in Figure 7. The results show that the fuzzy relationship between the parameter Q and the parameter Em can be expressed indirectly. In addition, as shown in Figure 7, the parameter Q increases with the increase in the ratio of B/S at the same height of the goaf h. In other words, the worse the integrity of rock mass, the greater the parameter Q, which is consistent with the engineering practice.Figure 7
Dimensionless flux charts of the fractured rock mass parameters: (a)B/S = 1.0 × 10−5. (b) B/S = 5.0 × 10−5. (c) B/S = 1.0 × 10−4. (d) B/S = 5.0 × 10−4. (e) B/S = 1.0 × 10−3. (f) B/S = 5.0 × 10−3. (g) B/S = 1.0 × 10−2. (h) B/S = 1.0 × 10−1.
(a)(b)(c)(d)(e)(f)(g)(h)Certainly, the relationship among the parameters ofQ, A, and Em can be obtained under a certain working condition. Assuming that the height of the goaf h is a constant, we can obtain the relationship between the parameter Q and the parameter A when the value of dimensionless expression Em/νγH is given. As shown in Figure 8, when h = 20 m, the parameter Q increases with the exposed area of the goaf A, which is in line with the actual situation. Similarly, the relationship between the parameter Q and the parameter Em is shown in Figure 9. In the same working condition, the parameter Q decreases with the increase of the parameter Em. In fact, the seepage flow Q is measured easily in the field. Therefore, we can use the measured parameter Q to estimate the parameter Em under specific working conditions.Figure 8
Diagram of the relationship betweenQ and A.Figure 9
Diagram of the relationship betweenQ and Em.Certainly, the relationship among the parameters ofQ, A, and Em can be obtained under a certain working condition. Assuming that the height of the goaf h is a constant, we can obtain the relationship between the parameter Q and the parameter A when the value of dimensionless expression Em/νγH is given. As shown in Figure 8, when h = 20 m, the parameter Q increases with the exposed area of the goaf A, which is in line with the actual situation. Similarly, the relationship between the parameter Q and the parameter Em is shown in Figure 9. In the same working condition, the parameter Q decreases with the increase of the parameter Em. In fact, the seepage flow Q is measured easily in the field. Therefore, we can use the measured parameter Q to estimate the parameter Em under specific working conditions.Figure7 contains the relationship among the physical quantities related to the fractured rock mass parameters of Zhengxing Iron Mine under all possible scenarios. Then, we can obtain Em under a certain working condition by Figure 7. Taking the mining level of −200 m in the mine as an example, the hydraulic conductivity K is 1.65 × 10−8 m/s with the help of field investigation and laboratory experiments. In addition, the average density of the rock mass in the mining area is 10 fractures per meter by measuring the fracture density on-site, and the average spacing of the fractures is calculated as S = 0.10 m. Combined with the parameter values in Table 1, we can obtain(37)B=12μkSKg1/3=12×1.01×10−6×0.1×1.65×10−89.81/3=1.27×10−5m.Selecting a goaf at the mining level as the research object, we can measure that the exposed area of the goafA is 5675 m2, the height of the goaf h is 20 m, and the seepage flow of the goaf Q is 16.74 m3/d. Therefore,(38)Ah2=5675202=14.19,Qgh5=16.7424×3600×9.8×205=3.46×10−8.According toB/S = 1.27 × 10−4, we can utilize Figure 7(c) to calculate the elastic modulus of fractured rock mass Em. As shown in Figure 10, the characteristic point G (14.19, 3.46 × 10−8) represents the coordinate point of the exposed area and seepage flow of the goaf. The interval [4.09, 4.14] of the right coordinate axis is divided into 5 equal parts to accurately obtain the corresponding elastic modulus of fractured rock mass Em at G point. Making a corresponding curve through the point G, as shown in the black dashed line in Figure 10, we can obtain that the value of log10Em/νγHis 4.14. The trend of the curve through the point G is similar to that of the red curve in Figure 10. Considering that the elevation of Earth surface is +50 m, we can obtain the depth of the goaf H is 250 m. Finally, we calculate that the parameter Em is equal to 23.29 GPa.Figure 10
The characteristic point and curve at the mining level of −200 m.Most scholars carried out many works to obtain more practical elastic modulus. The empirical estimation method establishes the calculation formula of rock mass elastic modulus by combining engineering experience and theoretical analysis. It has the advantages of simplicity and practicality.Coon and Merritt [2] and Gardner [3] proposed the expression of the relationship between the ratio Em/E and the rock mass quality index RQD, as follows:(39)EmE=0.0231RQD−1.32.Bieniawski [4] suggested that the elastic modulus of rock mass is correlated with the rock mass quality index RMR. When RMR > 50, Em can be expressed as(40)Em=2RMR−100.When RMR ≤ 50, Serafim and Pereira [5] suggested that Em can be expressed as(41)Em=10RMR−10/40.Based on the database of rock mass deformation modulus measurements, Hoek and Diederichs [6] proposed that the elastic modulus of rock mass Em can be represented as(42)Em=E0.02+1−2D1+e60+15D−GSI/11,where GSI is the geological strength index. D is a factor that depends upon the degree of disturbance to which the rock mass has been subjected to blast damage and stress relaxation.The index parameters of rock mass in the empirical estimation method introduced in this work are shown in Table3. According to the index parameters of the Zhengxing Iron Mine, the rock mass elastic modulus Em corresponding to each empirical formula can be calculated. As shown in Table 3, the value of Em calculated by the empirical formula is close to the calculation result of the method in this work, which indicates that it is feasible to use the dimensional analysis coupled with numerical simulation to calculate the elastic modulus of the fractured rock mass. In addition, when the working conditions are clear, we can directly obtain the value of Em from the dimensionless flux charts according to the seepage flow in the goaf, which facilitates the calculation of the elastic modulus on-site.Table 3
The index parameters of rock mass in the field and the calculated value ofEm.
ParametersRQDRMRDGSIEm (GPa)Coon and Merritt [2] and Gardner [3]Bieniawski [4]Hoek and Diederichs [6]Value78630.176923.1026.022.20


## 6. Conclusions

The aperture and spacing of fractures, which are closely related to the hydraulic conductivity, have the most direct influence on the elastic modulus of the fractured rock massEm. The proposed equivalent model is useful to analyze the mechanical behavior of fractured rock mass and conceptually represents the relationship between the hydraulic conductivity and the strain. Considering that the joint stiffness parameters are elusive, the hydraulic conductivity can be obtained conveniently by introducing a modulus reduction ratio Rm.The dimensionless flux charts developed in this work integrate most physical quantities involved in the physical process of groundwater flow into goaf. With the help of dimensional analysis coupled with numerical simulation, the elastic modulus of the fractured rock massEm can be calculated. When the working conditions are specific, the parameter Em can be obtained by querying the relationship diagram of the seepage flow Q and the parameter Em. Taking the Zhengxing Iron Mine as an example, the elastic modulus for the fractured rock mass calculated by the method in this paper is close to the calculation result of the existing empirical formula, which shows that it is feasible to use the method to calculate the parameter Em. In addition, the Python introduced in this work, which is a useful tool to handle the data generated under different working conditions, provides an automated program to realize the automatic execution for numerical simulation of all scenarios.In conclusion, the method of the dimensional analysis coupled with the numerical simulation throws new light on the reasonable calculation of the parameterEm. Further field measurements and observations contribute to the straightforward practical use of the method.



---

*Source: 2803837-2021-02-13.xml*