Abstract

This research experimentally studied the effects of various fracture roughness (characterized by the fractal dimension ) and normal stress (normal loads ) applied to fracture on ultrafine cement grout nonlinear flow behavior through rough-walled plexiglass fractured sample. A high-precision and effective sealing self-made apparatus was developed to perform the stress-dependent grout flow tests on the plexiglass sample containing rough-walled fracture (fracture apertures of arbitrary variation were created by high-strength springs and normal loads according to design requirements). The real-time data acquisition equipment and high-precision self-made electronic balance were developed to collect the real-time grouting pressure and volumetric flow rate , respectively. At each , the grouting pressure ranged from 0 to 0.9 MPa, and the normal loads varied from 1124.3 to 1467.8 N. The experimental results show that (i) the Forchheimer equation was fitted very well to the results of grout nonlinear flow through rough-walled fractures. Besides, both nonlinear coefficient () and linear coefficient () in Forchheimer’s equation increased with increase of and , and the larger the was, the larger the amplitude was. (ii) For normalized transmissivity, with the increase of Re, the decline of the curves mainly went through three stages: viscous regime, weak inertia regime, and finally strong inertia regime. For a certain , as the normal load increased, the curves generally shifted downward, which shows good agreement with the single-phase flow test results conducted by Zimmerman. Moreover, with the increase of , the Forchheimer coefficient decreased. However, within smaller , decreased gradually with increasing and eventually approached constant values. (iii) At a given , increased with increasing .

1. Introduction

Fractured rock masses (Figure 1(a)) are widely distributed in underground engineering due to geological action or excavation disturbance, which may bring serious safety hazard to the engineering stability [14]. Grouting is a useful way for controlling groundwater inrush and improving mechanical properties of fractured rocks in underground mining and engineering [59]. In order to acquire a better understanding of the grouting mechanism of fractured rock, lots of experimental studies have been performed in recent decades. For example, Sui et al. [10] experimentally carried out the study on the sealing effect of permeation grouting into specimen fractures and grout propagation patterns under water flow conditions, which obtained the influence of different factors on the sealing effect, including the initial water flow speed, fracture aperture, grout take, and gel time. Lee et al. [11] studied the reinforcement effect of cement-grouted jointed rock masses. The results showed that after grouting, the stiffnesses of the filled joints were increased up to 6 times compared with those of the ungrouted joints, meaning the mechanical properties of the jointed rock masses after grouting were significantly improved than before. Funehag and Fransson [12] investigated a large number of experimental studies and field measurements on the penetration grouting into fractured natural rocks, which obtained the grout diffusion radius. In addition to experimental research, Kim et al. [13] conducted a series of numerical analysis by using UDEC to simulate the grout diffuse through smooth parallel surfaces based on the Bingham fluid model. They found that the grout penetration length was overestimated without considering the time-dependent hardening of grout material. Hässler et al. [14] proposed a model to simulate grout propagation in fractured rock considering time-dependent properties of the grout.

However, the previous studies on grouting mechanism mainly focused on the grouting reinforcement or seepage prevention effect for the fractured rocks after grouting (e.g., mechanical properties, sealing effect of grouted fractured rocks, and grout penetration length and diffusion radius). There were less studies on the flow behavior of cement grout through the rock fracture during the grout flow process before grout hardening. Similarly, for water, the flow behavior through rock fractures is governed by the well-known Navier-Stokes (NS) equations and the mass conservation equation written as follows [15, 16]: where stands for the velocity vector and mean the components of velocity vector; , , and refer to the fluid density, dynamic viscosity, and total hydraulic pressure, respectively. Nevertheless, the theory of flow through rough surface of rock fractures has not been fully developed due to the fracture roughness and fluid nonlinear flow characteristics [17]. The existence of the complex nonlinear partial differential equations and irregular rock fracture geometry makes it very difficult to solve the NS equation through real rough-walled fractures [18, 19]. For simplifying the calculation, some flow equations with assumed conditions were presented, among which the best-known equation is the cubic law that neglects the inertial terms of fluid flow through the fractures [20]. In the cubic law, the fluid flow through two smooth parallel plates was thought to be the viscous flow, and the total volumetric flow rate was linearly dependent on the grouting pressure gradient : where and are the fracture width and the hydraulic aperture, respectively. This simplified formula only applied to the laminar flow through the smooth fracture surface. However, the natural fracture is usually rough and often subjected to field stresses directly related to the fracture aperture and, hence, the fluid flow behavior through the fractures [21]. Fracture roughness has an effect on fluid flow through fractures, and the void space between opposing surfaces (fracture aperture) decreases due to the increasing normal loads [2225]. Therefore, the characteristic of the natural rough-walled fractures under stress condition makes the actual flow rate seriously inconsistent with the results calculated using the cubic law mentioned above (equation (2)) and may cause the nonlinear flow behavior through the rough-walled rock fracture. There were many studies on the nonlinear flow in rough-walled fractures. Chen et al. [26] experimentally studied the flow through rough-walled fractures and made great contribution to the indicators representing fracture geometric characteristics. Chen et al. [27] studied the effect of geometric characteristics of deformable rough fractures under confining stresses on the non-Darcy flow behaviors. They found that the Forchheimer equation can be used to describe the non-Darcy flow in rough fractures well. Zou et al. [28] simulated the nonlinear flow in 3D rough-walled rock fractures by solving the Navier-Stokes equations and investigated the influence of shear-caused aperture changes and flow conditions (inertial term) on nonlinear flow behavior, which observed that channeling flowed along the preferential paths, transverse flowed around the contact spots, and eddy flowed behind contact spots with increasing Re. Moreover, they also studied the impacts of wall surface roughness on fluid flow through rock fractures and obtained that the high-frequency secondary roughness had a major effect on the dynamic evolution of eddy flow regions in the fracture flow field, except for the Reynolds number (Re) [29]. Li et al. [30] experimentally investigated the hydromechanical behavior of rock joints by using a parallel-plate model containing contact areas and artificial fractures. They presented some empirical relations to evaluate the influences of contact area and surface roughness on the fluid flow behavior through rock fractures. From these studies, we found that the formula accepted by most researchers for nonlinear flow through rough-walled fractures is Forchheimer’s equation [31]: where and are coefficients, representing the pressure drop components caused by nonlinear and linear effects, respectively. These terms are affected by the fracture geometry and hydraulic gradient/pressure. As for the cement grout, although the grout is also a liquid before hardening, the grout flow behavior through the rock fracture is more complex than the water: due to the existence of the cement in grout, the physical properties of grout such as viscosity and density are significantly different from those of water.

In this paper, in order to study the grout flow behavior through the rough-walled fractures, a high-precision and effective sealing self-made apparatus was developed to perform stress-dependent grout flow tests on the plate specimens with fractures made of high-transparency plexiglass and high-strength springs. The plexiglass samples containing rough-walled fractures (various fracture apertures can be created by high-strength springs and normal loads) with different (the rough morphology of the single fracture was characterized by using the fractal dimension ) were prepared. The rough morphology of the fracture was cut using a high-precision laser cutting machine. A series of tests of cement grout flow through the rough-walled fractures were performed with respect to different grouting pressures ranging from 0 to 0.9 MPa, different normal loads varying from 1124.3 to 1467.8 N, and different fracture fractal dimension (1.1, 1.2, 1.3, 1.4, and 1.5). Based on the test results, the cement grout nonlinear flow behavior through rough-walled rock fractures with different and normal loads was analyzed.

2. Grout Flow Testing System

A high-precision and effective sealing self-made apparatus was developed to study the stress-dependent grout flow through fractured rocks, as shown in Figures 2 and 3. This test apparatus mainly contains three parts.

2.1. A Platform for Grout Flow through Fractured Rock Sample

The platform consists of the thick bottom steel plate ( mm) which can be used for the placement of the fractured rock sample (Figures 2(a) and 2(b)). The normal loading device ( mm) connected to the hydraulic jack provided the normal load to the fracture containing high-strength springs (Figure 2(b)). Meanwhile, the grouting pipes (5 mm in internal diameter and 9 mm in length) passing through the rectangular holes belonging to another loading device providing tangential force were inserted into the fracture and tightly sealed by using CR4305 chloroprene rubber [3234] which ensured the sealing effect of the fracture and prevented the grout overflowing from the entrance of the fracture (Figure 4(b)). A high-strength and transparent plexiglass cover plate which was pressed under two adjustable bolt cover plates as shown in Figures 2(a) and 2(b) was placed on the fractured sample. Before each test, the plexiglass rough-walled fractured samples ( mm in size) with different (1.1, 1.2, 1.3, 1.4, and 1.5; see Figure 4(c)) were placed on the platform, and the horizontal and vertical loads were applied to the fractured sample boundaries.

2.2. Grout Supply and Real-Time Data and Image Acquisition System

In order to obtain the stable grouting pressure, the high-precision nitrogen regulator was installed on the nitrogen tank to achieve stable nitrogen source pressure (Figure 3(a)). Then, the stable nitrogen drove the grout to flow at a steady grouting pressure in a certain amount of time. The grouting pressure ranged from 0-0.9 MPa which provided different grouting pressure gradient required for the tests of grout flow through fractures. The real-time data acquisition equipment consists of the high-precision pressure transmitter and paperless recorder (Figure 3(a)) which were used to collect the real-time grouting pressure . The volumetric flow rate out of each fracture was derived when the quality change rates of grout flowing out of fracture were roughly the same, and the real-time grout fluid quality was obtained by the self-made high-precision electronic balance (see Figure 3(b)).

2.3. Load Supply System

The horizontal (- and -directions) loads on the fracture specimen were produced using adjustable hydraulic jacks (Figure 2(b)). The horizontal loading device connected to load jacks (pressure head) can apply normal uniform loads on the boundaries of the fractured sample. Moreover, the high-precision pressure sensors (Figure 2(b)) installed between the pressure head and the counterforce frames were used to feedback the values of the horizontal loads, thus guiding the adjustable hydraulic jacks to provide the fracture specimen with accurate normal boundary loads. Vertical loads applied on the specimen surface created by high-strength and transparent plexiglass cover plate were for balancing the vertical grouting pressure in the fracture and then further sealing the fractured sample.

3. Material Preparation

3.1. Grout Material

The ultrafine cement grouts (W/C ratios of 1.0 with (kg/m3) and (Pa)) were selected in the experiment. The ultrafine cement grouts have been widely applied in the reinforcement for underground engineering due to its more stable physical and mechanical properties, greater range of strength adjustability and better permeability, and nonpolluting nature [35]. Please note that because of the time of the grout flow through fracture in the test being very short, the grout has been in good flow state during the test. Therefore, the time-dependent hardening of grout was not considered in this paper.

3.2. Sample Material

For better conducting the visualization study of grout flow process through fracture and preset the roughness of plexiglass fractured sample which was similar to that of natural fracture, the plate fractured rock samples ( mm in size) made of high-transparency plexiglass were prepared. Fracture morphologies are mostly irregular in actual underground engineering [3639]. In order to simulate the irregular characteristics of fracture in nature, Ju et al. [22] presented a fractal model of single rough fracture based on the methodology. The rough morphology of the single fracture was characterized by using the fractal dimension which was based on the Weierstrass-Mandelbrot fractal function [40] where the parameter is a real number larger than 1, reflecting the degree of deviation of the curve from the straight line; is an arbitrary phase angle; Fractal dimension. Using the real part of the function as the fractal control function , where the function is a fractal curve that is continuous, nondifferentiable, and -dimensional. The fractal dimension satisfies where represents a constant; represents Hausdorff-Besicovitch dimension. MATLAB was used to generate the fractal curve with different fractal dimensions on the computer, and then, the fractal fractures (rough morphologies) were cut according to the function on the plexiglass mold by using the high-precision laser cutting machine (Figure 4(a)). , , , , and were selected to study the effect of fracture roughness on the grout flow behavior through rough-walled fractures, and the higher the fractal dimension, the larger the fracture roughness was (Figure 4(c)). The high-strength springs that were used to simulate the variations in the fracture aperture were the important component of the fractured sample (Figure 4(b)). On the plexiglass sample, the reasonable setting of sizes of inlet (outlet) and springs ensured the uniform distribution of normal loads to the sample boundary (see Figure 4(b)).

3.3. Sealing of the Fractured Sample

For ensuring that no air or grout leakage occurred under grouting pressure condition, sealing of the fracture in the sample was the key to the success of the grouting test. The “CR4305 chloroprene rubber” was selected in this experiment after many attempts (Figure 4(b)). The CR4305 chloroprene rubber is a compressible and waterproof material that is compressed under normal loads, which seals the entrance and exit of the fracture. A piece of transparent crystal plate of a suitable size and 2 mm thickness was attached on the fractured sample, and then, the fractured sample with the transparent crystal plate was covered with a high-strength and transparent plexiglass cover plate of a suitable size and 30 mm thickness. When the fracture was preliminarily pressed tight by the transparent plexiglass cover plate, the sealed fractured sample was placed on the platform, and then, the horizontal loading device applied horizontal loads to the sample. In the grout flow testing apparatus, the transparent plexiglass cover plate provided uniform and powerful vertical loads generated by two strong bolt cover plates to further seal the fracture and balance the vertical grouting pressure in the fracture (see Figures 2(a) and 2(b)).

4. Results and Analysis

4.1. Stress-Dependent Grout Flow Behavior through Rough-Walled Fracture

From the theory of hydraulics, is defined as the ratio of the hydraulic head difference to the flow length, calculated as follows: where here is called grout hydraulic gradient in this paper. is the grouting pressure obtained by the pressure transmitter and paperless recorder in real time. is the grout density ( (kg/m3)). is the flow length that is 0.24 m. Note that is linearly related to at the fluid inlet: . The hydraulic head at the grout outlet of the fracture was assumed to equal zero. Thus, the grouting pressure gradient ranged from 0 to 3.75 MPa/m.

From Figure 5, the relationship between and showed obvious nonlinear characteristics. The linear Darcy’s law can no longer be used to describe the nonlinear grout flow behavior through the rock fracture. The zero-intercept quadratic equation used to describe the nonlinear flow behavior through fractures was presented by Forchheimer in 1901 (equation (3)). This nonlinear flow model has been widely adopted by many scholars [24, 41, 42]. Combined with equation (7), equation (3) can be written to where and

The relationship between the grout hydraulic gradient and the volumetric flow rate of grout flow through rough-walled fracture with different subjected to various values ranging from 1124.3 to 1467.8 N was displayed in Figure 5. Referred to equation (8), the experimental data were well fitted and analyzed by the Forchheimer equation (Figure 5 and Table 1); the values of the correlation coefficient for all cases were larger than 0.99, as listed in Table 1, which shows that all of the experimental values almost agreed well with the fitting curves of zero-intercept quadratic equation. From Figure 5, at a given , with the decrease of , all of the maximum volumetric flow rate increased, meaning has an obvious effect on the grout flow through rough-walled fracture. For example, was increased from () to (), increasing by 0.68 times (); to by 0.90 times (); to by 0.69 times (); to by 1.54 times (), respectively. Meanwhile, at a given , as was increased (decreased), all of the maximum volumetric flow rate decreased, meaning also has an obvious effect on the grout flow through rough-walled fracture. For example, as ranged from 1124.3 N to 1467.8 N, the maximum volumetric flow rate decreased from to , decreasing by 95.99% (); to by 96.32% (); to , by 94.85% (); to , by 94.10% (); to , by 93.95% ().

In equation (8), the nonlinear and linear coefficients and were both calculated (see Table 1). From Figure 6, both and of all increased with increasing . Both coefficients and showed increasing trends with increasing , and the variation trend was generally consistent with the experiments performed by Yin et al. [43]. Moreover, the larger the was, the larger the amplitude was. Taking coefficients and for under varying as an example, for , 1353.3, and 1467.8 N was , , and , respectively, increasing by 7.16, 186.23, and 357.50 times over the value of for . Meanwhile, for , 1353.3, and 1467.8 N was , , and , respectively, increasing by 6.68, 21.22, and 32.31 times over the value of for .

4.2. Normalized Transmissivity

The transmissivity was used to estimate the nonlinear flow regimes for rock fracture samples in previous research works [44, 45]. However, is a constant value related to the permeability and flow area of the fractured rock in Darcy’s law. In this study, for grout nonlinear flow through rough-walled fractures with different and subjected to various normal loads , was expressed as follows:

In order to evaluate the flow behavior through the rock fractures, the Reynolds number (Re) is defined as the ratio of inertial forces to viscous forces, which is typically used to quantify the onset of nonlinear flow:

For analyzing the nonlinear flow behavior through the fracture, and were introduced in the study [45], where the Forchheimer formula (equation (3)) was rewritten: where (obtained when the flow rate was sufficiently low and the inertial force was negligibly small) refers to the intrinsic transmissivity and refers to the apparent transmissivity. is the Forchheimer coefficient (Figure 7), which is expressed as referred to the Forchheimer equation. Based on , the relations between and Re were fitted and analyzed with different and (Figure 8). From Figure 8, for all cases, when Re was small, the viscous effect was dominant in the grout flow behavior, and presented a horizontal straight line (close to 1.0), the stage of which was called the viscous stage (taking Figure 8(e) as an example, Red number 1). As Re increased, the inertial effect gradually increased, and started to bend downward. However, during this stage, the overall grout flow behavior was still dominated by the viscous effect, and the inertial effect can be ignored, the stage of which was called the weak inertial effect stage (Figure 8(e), Red number 2). When Re increased to a certain value, started to decline linearly due to the strong inertial effect. At this time, the inertial effect cannot be ignored, the stage of which was called the strong inertial effect stage (Figure 8(e), Red number 3). Meanwhile, for a certain , as the normal load increased, the curves generally shifted downward. Although there were two different fluids, this variation trend of curves showed good agreement with the single-phase flow test results conducted by Zimmerman et al. [45]. From Figure 7, with the increase of , most of the values of decreased. This decreasing trend was especially evident when the normal forces were larger (e.g., and 1467.8 N). Taking larger as an example, with the increase of from 1.1 to 1.5, decreased from 0.02018 to 0.01372, or by 32.01%, while decreased gradually with the increase of , when was smaller (such as 1124.3 N). However, when was 1238.8 N, as increased, was almost unchanged.

4.3. Critical Grout Hydraulic Gradient Analysis

The flow regime does not obey the linear Darcy’s law under higher fluid flow velocity (higher pressure gradient ) [27, 46, 47]. In order to further illustrate the non-Darcy flow mechanism of fluid flow through rough-walled fractures and quantify the degree of nonlinear flow effect, a factor was proposed by Zeng and Grigg [48]: where and were energy losses due to the inertial and viscous dissipation mechanisms in the fractures, respectively. Thus, the factor expressed the percentage of the grout hydraulic gradient drop related to the nonlinear term to the total grout hydraulic gradient. For grout flow through rock fractures, nonlinear flow effects may become more obvious when the volumetric flow velocity increased at the grout inlet. Previous studies have shown that in practical engineering, when the hydraulic gradient drop caused by the nonlinear flow effect was greater than 10% of the total gradient drop, the nonlinear term cannot be ignored [49, 50]. Therefore, the critical factor was selected to evaluate the flow regime through fractures, and the corresponding is defined as the critical hydraulic gradient . The variations in for all cases with different and were as shown in Figure 9 and Table 1. From Figure 9, at a given , as increased, for all cases showed increasing trends, which means the fracture roughness has an important impact on the flow regime through fracture. The reason for these variations may be as follows: with the increase of , the fracture was rougher, which would make it harder for the grout to flow through the fracture. Therefore, the grout hydraulic gradient required for grout flow entering into the nonlinear flow regime was larger. Taking as an example, was 128.390 (), 132.593 (), 165.246 (), and 209.498 (), respectively, increasing by 43.66%, 48.37%, 84.90%, and 134% over of 89.369 for . Please note that, for magnitude, the magnitude of in our study was not same as that in Yin et al. [43], and the difference was about an order of magnitude ( in our study was larger), which may be due to the viscosity of the cement grout being higher than that of water, grout entering the nonlinear flow regime required larger grout flow rate than that of water, so of grout would become larger accordingly.

5. Conclusions

The effects of the fracture roughness (characterized by the fractal dimension ) and normal loads on the ultrafine cement grout nonlinear flow behavior through rough-walled fracture were studied. A self-made apparatus with high precision and effective sealing characteristics was used for stress-dependent grout flow tests. The nonlinear flow behavior (-), nonlinear flow regime, nonlinear and linear coefficients ( and ), normalized transmissivity (), Forchheimer coefficient (), and critical hydraulic gradient () for grout flow through fractures were obtained. The main experimental results were as follows: (1)The correlation between grout hydraulic gradient and volume flow rate was demonstrated to be nonlinear based on the flow tests, which was fitted by the Forchheimer equation. At a given , with the decrease of , all of the maximum volumetric flow rate increased, meaning the fracture roughness has an obvious effect on the grout flow through rough-walled fracture. Meanwhile, at a given , with the increase of , all of the decreased, meaning the normal stress also has an obvious effect on the grout flow through rough-walled fracture. Moreover, both nonlinear coefficient () and linear coefficient () of all increased with increasing and , and the larger the FN was, the larger the amplitude(2)For normalized transmissivity, when Re was small, with the increase of Re, the curves presented a horizontal straight line (close to 1.0), the reason of which was that the viscous effect was dominant in the grout flow behavior through fracture; as Re continued to increase, started to bend downward. During this stage, the grout flow behavior was still dominated by the viscous effect, and the inertial effect can be ignored. When Re increased to a certain value, started to decline linearly. During this stage, the inertial effect cannot be ignored, and it became dominant in the grout flow behavior through fracture. Moreover, for a certain , as the normal load increased, the curves generally shifted downward. The curve variation trend showed good agreement with the single-phase flow test results conducted by Zimmerman. For the Forchheimer coefficient , with the increase of , decreased. Within smaller , decreased gradually with the increase of and eventually approached constant values(3)At a given , increased with the increase of , which means the fracture roughness has a significant effect on the flow regime through fracture. The magnitude of in our study was an order larger than that in Yin et al. [43], which may be because the cement grout viscosity was higher than that of water, grout entering the nonlinear flow regime needed larger grout flow rate than that of water. Therefore, of grout would be larger accordingly

In the future research, the split grouting in the deep fractured rock mass and cement grout flow behavior through 3D fractures with contacts will be the focus of our study.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Acknowledgments

This work was supported by the joint Ph.D. program of “double first rate” construction disciplines of CUMT.