Simulating the Filtration Effects of Cement-Grout in Fractured Porous Media with the 3 D Unified Pipe-Network Method

In grouting process, filtration is the retention and adsorption of cement-grout particles in a porous/fractured medium. Filtration partly/even completely blocks the transportation channels in the medium, greatly decreasing its permeability. Taking into account filtration effects is essential for accurately estimating the grout penetration region. In this paper, the 3D unified pipe-network method (UPM) is adopted for simulating 3D grout penetration process in a fractured porous medium, considering filtration effects. The grout is assumed to exhibit two-phase flow, and the filtration effects depend on not only the concentration and rheology of the grout but also the porosity and permeability of the fractured porous medium. By comparing the model with the experimental results, we firstly verify the proposed numerical model. Then sensitivity analysis is conducted, showing the influences of grout injection pressures, the water?cement ratios of grout (W/C) and the grout injection rates on filtration effect. Finally, the grout filtration process in a complex 3D fractured network is simulated, indicating that the size of the grout penetration region is limited due to filtration. Record Type: Published Article Submitted To: LAPSE (Living Archive for Process Systems Engineering) Citation (overall record, always the latest version): LAPSE:2019.0533 Citation (this specific file, latest version): LAPSE:2019.0533-1 Citation (this specific file, this version): LAPSE:2019.0533-1v1 DOI of Published Version: https://doi.org/10.3390/pr7010046 License: Creative Commons Attribution 4.0 International (CC BY 4.0) Powered by TCPDF (www.tcpdf.org)


Introduction
Cement grouting is a widely used method for sealing underground water and improving the mechanical characteristics of surrounding media, such as fractured rocks, in geotechnical engineering.Over the past few decades, many on-site experiments have been conducted to assess the grouting quality [1][2][3][4][5][6] and to improve the strength and penetration distance.Some researchers remarked that cement particles can accumulate and deposit on the surfaces of pores/fractures [7][8][9], thereby blocking transport channels within the porous/fractured medium and terminating any further penetration of the grout.This phenomenon is called filtration, an irreversible process that significantly reduces the permeability of the porous/fractured medium [10][11][12].Hence, assessments of the distance over which grout has penetrated will be greatly overestimated if the effects of filtration are ignored.
Researchers have experimentally investigated the influences of filtration on grouting with one-dimensional column injection tests considering constant grout injection pressures, injection rates, or cement grout particle concentrations.These researchers found that the cessation of cement grout penetration depends not only on the physical characteristics of the grout [13][14][15], such as the concentration of the cement grout particles and rheology of the grout, but also on the properties of the porous/fractured medium [16,17], such as the pore shape and fracture distribution.Furthermore, the permeability of the grouted medium decreases substantially in penetration zones [18][19][20].Nevertheless, while the results of experimental investigations make it possible to verify proposed analytical and numerical models, their main disadvantage is that the influencing factors are very difficult to control as planned [21].
Many numerical and analytical methods, including phenomenological models [22], continuum models [23], trajectory analysis models [24], stochastic models [25], and network models [26][27][28][29][30], have been proposed to analyse the effects of filtration in porous but not fractured media.With these methods, some important conclusions have been drawn; for example, Saada et al. [19] and Maghous et al. [31] found that small variations in the porosity can lead to high variations in the permeability, as the injection of grout towards a void can be stopped by cement particle filtration, after which no significant reduction in the porosity will occur.Researchers have also proposed single [32] and multi-phase [19] models that are all applicable, but they have shown different precisions.
Moreover, research has been presented on the simulation of grout flow considering filtration effects in a fractured, porous medium [13,33].Generally, fractures exist naturally within rock-like materials; thus, the initiation and propagation of such fracture need to be modelled by sophisticated numerical tools, such as the eXtended Finite Element Method (XFEM) [34][35][36], the strong discontinuity embedded approach [37][38][39][40][41][42], the mixed-mode FEM [43][44][45][46], and the phase-field model [47][48][49].However, when focusing only on mass/energy transport within existing fractures and not their mechanical evolution, much simpler models can be used instead of these sophisticated techniques.The unified pipe-network method (UPM) is one of these simplified methods; it uses equivalent pipe networks to simulate the mass/energy transport processes within a porous/fractured medium.With the UPM, 2D and 3D mass/energy transport problems can be transformed into equivalent 1D problems, similar to the lattice element method (LEM) [50][51][52].The fractures and porous medium can then be discretized using the same pipe-network system [53,54], which represents a significant advantage that provides a higher computing efficiency and greater numerical stability than other methods.
In this work, the grout filtration process within a 3D porous and fractured rock mass is simulated with a 3D UPM model.The grout flow is assumed to be an incompressible liquid [13,15,27].In comparison with our former work [55], which focused on grout flow with a considerable time-dependent viscosity, we focus on the blocking effects of cement particles correlated with the microscopic structures of fractured, porous media.The remainder of this paper is organized as follows.In Section 2, the 3D two-phase filtration model is clearly introduced, and both the fluid flow and the mass transportation equations are given.In Section 3, the discretized equations of the filtration model based on the 3D UPM are established, thereby covering the interactions between the grout flow and the rock matrix and fractures.In Section 4, the developed numerical model is verified through a comparison with the experimental results given in [56] in tandem with a sensitivity analysis of the filtration process, including the grout injection pressure, grout injection rate and grout water-cement ratio.In Section 5, grout penetration in a porous and fractured rock is simulated considering different fracture inclination angles and random fracture networks.In Section 6, concluding remarks are provided.

Mathematical Model
We consider two-phase (fluid and solid particles) flow to study the filtration processes in a fractured and porous rock.This model is based on the interaction of Darcy's flow and mass transportation.The cement-based grout is considered to be a Newtonian fluid.The saturated single-phase grout flow in both the fractures and rock matrix obeys Darcy's law and can be calculated based on a mass conservation equation: where τ is a sign expressing the different media in the model domain, for example, τ = m expressing a matrix and τ = f expressing a fracture, where φ m is the porosity of the rock matrix and φ f is the porosity of a fracture; ρ is the grout density; q is the source term; and u is the flow velocity vector, which can be expressed as in which µ is the grout fluid viscosity, P is the grout pressure, g is the gravitational acceleration and k τ is the rock matrix and intrinsic fracture permeability tensor.The aperture of the fractures in the current simulation is assumed to be smooth, and the laminar flows of the grout in the fractures can be described by a cubic law; thus, the intrinsic fracture permeability can be calculated as k f = a 2 /12, where a is the fracture aperture.
For analyzing the transport of cement grout particles through a porous medium, the advection-dispersion equation is used with an additional term (∂σ/∂t), which is defined as the rate of suspended particle mass transfer between the grout and the pore skeleton.Although the hydrodynamic dispersion coefficient is typically non-zero for multiphase flows in porous media [57], researchers have noted that the hydrodynamic dispersion tensors of porous media and fractures can be ignored since the effect of hydrodynamic dispersion on the dilution of grout is generally small [10,19,[58][59][60].Hence, the following relation is obtained for a porous medium: where C is the concentration of transported particles, D m is the hydrodynamic dispersion tensor of the porous medium, λ is the constant deposition coefficient and σ is the mass of grout particles being deposited per initial unit pore volume.In the above equation, the inhibition of particle deposition by inter-particle forces is assumed to be minimal [61], and the filtration process is modeled through a kinetic equation for mass transfer and represented by a first-order relationship [10,18,56].
Similarly, when assuming that filtration cannot proceed after cement grout particles have migrated into the fractures, the following relation is obtained for fractures: The filtration of cement grout particles is enhanced by a decrease in the volume of the void space.When the grout is injected into the porous medium, the porous medium becomes plugged, and the porosity decreases with space and time.The porosity is thus assumed to decrease in proportion to the mass of filtered particles as follows: where φ m 0 is the initial porosity and ρ s is the rock density.The intrinsic permeability of the porous rock decreases with the clogging of the pore space.The experimental work performed by [19] proved that the permeability is dramatically influenced by the porosity, which means that very low variations in the porosity may induce large variations in the permeability.However, the well-known Kozeny-Carman equation expresses that the change in the permeability is consistent with the change in the porosity.Thus, a hyperbolic function of the (filtration-induced) porosity change is adopted in this paper to express the relationship between the porosity and permeability: where k m 0 is the initial intrinsic permeability and b is a negative scalar.The concentration of cement grout particles has a considerable effect on its rheological properties [62].Except for the porosity φ τ and intrinsic permeability k τ , the viscosity µ and the grout density ρ vary with time because of the cement solid filtration.The effective viscosity and density of the grout can be expressed as where µ w is the viscosity of water, ρ w is the density of water, δ is the volume fraction of cement particles, which can be calculated as δ = C/ρ s , and ρ c is the density of the cement particles, where ρ c = ρ s .

The Discretization Model for the 3D Unified Pipe-Network Method
In the framework of the 3D UPM, the equivalent hydraulic parameters and diffusion parameters of the fracture and matrix pipes are derived based on discrete unstructured elements.

Matrix Pipe-Network Model
The 3D porous medium is discretized using tetrahedral elements in this model, where each tetrahedron has six edges, and these elements are regarded as the flow pipes of the matrix, as shown in Figure 1a.When assuming that the edges of the tetrahedral elements are equivalent pipes, the mass/energy transport problem in a fractured and porous medium is transformed into an equivalent transport problem in the pipe network system.To obtain the equivalent parameters of the flow matrix pipes, the center of the sphere circumscribed around the tetrahedron (point o) and the circumcenters of the four faces of the tetrahedron (points cn, (n = 1, 2, 3, 4)) are determined as division points (see Figure 1b).Linking the point o to the points c n and to the midpoint of their joint edge in sequence encloses a partition face of the tetrahedron.The tetrahedron will form six partition faces in total, each of which is perpendicular to its respective intersecting edges.The tetrahedron is thus divided into four parts by these six partition faces, and the total pore volume of each part is assigned to its corresponding node, see Figure 1a.The equations for grout flow and cement grout particle transport through the control volume of a node are derived by formal integration methods, such as the finite volume method (FVM), expressed as follows: where V is the nodal control volume and A is the nodal control area generating the 3D Voronoi element.The pressures and concentrations within a tetrahedral element can be written using the following interpolations: where N k is the linear shape function as in the FEM and P k and C k denote the pressures and grout concentrations, respectively.Here, we would like to emphasize that hydro-mechanical effects are ignored in the UPM; these effects require higher-order shape functions for the displacement field and lower-order shape functions for the pressure/concentration fields, such as the hybrid formulations used in [63][64][65], to avoid numerical instabilities [66,67].When using the same shape functions for the pressure/concentration fields in our examples, we do not encounter any numerical problems.
The grout fluid flow Q m ij in pipe ij is equal to the flow through the area oc1 f c2, see Figure 1a, and can be calculated as where n oc1 f c2 is the unit normal vector, A oc1 f c2 is the area of the face oc1 f c2, P i and P j are the pressures at nodes i and j, respectively, and K m ij is the equivalent conductance coefficient of pipe ij, which can be written as [68][69][70] where l ij is the length of pipe ij.Furthermore, fractures in the 3D domain are represented as 2D disks.
The equivalent hydraulic parameters of the pipes are derived from unstructured triangular elements, as shown in Figure 1b.The conductance coefficient K f ij for fracture pipe ij in the 3D domain is derived as [55] The grout flow in each node obeys Darcy's law.Thus, within the framework of the UPM, for each node i, the governing equation for single-phase saturated flow is where φ m i is the porosity of node i, V i is the control volume of node i, the subscript n i denotes the total number of connected pipes, and Q s i is the source term of node i.
The effective diffusion coefficient in a 3D porous medium can be derived using the same discretization method used for the aforementioned equivalent conductance coefficient.The mass of the grout particles ṁm ij transported in pipe ij is also equal to the mass being transported through the area oc1 f c2 and can be calculated as where C i and C j are the concentrations at nodes i and j, respectively.The grout particles being transported through each node obey the nodal laws.For each node i, the governing equation for mass transport is In the 3D fractured and porous rock, the fractures and rock matrix are reconstructed using 1D pipes.The pipes with the start and end nodes belonging to both the matrix and fracture pipes are defined as interface pipes, which means that both the fracture and matrix pipes will possess the same nodes at the interface.The nodal control volumes and the equivalent conductance coefficient and diffusion coefficients of the interface pipes are calculated by the superposition of the respective values of the fracture pipes and matrix pipes.In addition, the pore pressures and particle concentrations are assumed to be continuous at the interface without interaction between fractures and the matrix.Finally, the governing equations for each node are written as

Experimental Verification
Bouchelaghem [56] conducted a series of experiments on the flow and transport of microscopic to fine cement in sand using 1D column tests.These experiments use a well-graded and medium to coarse sand; thus, the porous medium can be considered to be homogeneous.A diagram of the grout injection experiment is given in Figure 2a.Following Bouchelaghem's experiments, the current numerical simulation uses a column with dimensions of 35 × 35 × 820 mm for a comparison with Bouchelaghem's results.The coupled model of fluid flow and mass transport is used to simulate the grout filtration process.The parameters are listed in Table 1.In this model, the grout injection rate and grout concentration are constant at the inlet boundary, while the pressure is a constant at the outlet boundary and is defined as zero.The initial conditions at the four other transverse boundaries are such that n • ∇P = 0 and n • ∇C = 0. Figure 3 shows the pore pressures at four elevations at different times under a constant injection velocity from both the experimental and the simulation results.The pore pressure evidently increases with time, and the pore pressures closer to the injection point are larger than those farther from the injection point.The pressure changes slowly far away from the injection point.The numerical results are in good agreement with the experimental results for the four elevations and are consistent with the numerical results obtained by Kim et al. [59] in comparison with the same experimental results.In summary, the key parameter causing the differences between the numerical simulations and experiments is the filtration rate, λ.

Sensitivity Analysis of Grout Filtration
The grout injection pressure, injection rate and water-cement ratio of grout are three significant engineering parameters influencing the grouting effect, and all three are usually considered before engineering projects are initiated.In this sensitivity study, the grout filtration process with different grouting operations is modeled.The model is a column with dimensions of 5 × 2 × 2 m, as shown in Figure 4.The meshes are generated by self-developed mesh techniques [71], and the total number of grids in this model is 103,792.The transient flow pressure, mass concentration, porosity and permeability of rock are obtained and compared in time and space.The results are compared by analyzing the points' values along line AB, as shown in Figure 4, which is defined as the central line of the simulation model.In the current simulation, grout is injected into a cuboid domain under a constant pressure and constant concentration at the inlet boundary, while the pressure at the outlet boundary is zero.The parameters used in this case study are listed in Table 2. Figure 5 compares the variations in the pore pressure, particle concentration, rock porosity and permeability under two conditions, namely, (i) filtration and (ii) non-filtration, for the purpose of studying the influence of filtration on the grout penetration process.Under the first condition, the grout penetration parameters and rock characteristic parameters are extracted at time steps of t = 800 s, t = 1600 s, t = 2400 s, t = 3200 s and t = 4000 s, whereas, under the second condition, they are extracted at a time step of t = 4000 s.When ignoring the filtration effects, the porosity and permeability of the rock matrix are constant, and the pore pressure and particle concentration decrease linearly.However, in reality, the filtration makes these parameters vary nonlinearly.When considering the filtration effects, Figure 5a shows that the pore pressure decreases quickly with increases in the penetration length and injection time, indicating that it is necessary to increase the injection pressure with time to ensure the designed grout penetration size.A remarkable decrease in the cement concentration is observed, and the position at which this decrease occurs moves away from the injection position with time, as shown in Figure 5b.Furthermore, the concentration of cement particles also decreases quickly with time, which means that, under the same injection concentration, the grout particles move more slowly with time.Figure 5c,d show that the porosity and permeability decrease with time and increase with the distance from the injection position, indicating that the variation rate of the permeability is larger than that of the porosity.

Effects of the Grout Injection Pressure
To analyze the influence of the grout injection pressure on the filtration, the injection pressure is varied from 30 kPa to 60 kPa, 100 kPa, 150 kPa and 200 kPa with a total injection time of 4000 s.As shown in Figure 6a, the pore pressure decreases quickly with increasing injection pressure.Figure 6b shows that the concentration decreases with the injection pressure.Although the injection pressures are different, the porosity and permeability of the rock matrix at the inlet boundary are the same, as shown in Figure 6c,d.This indicates that the injection pressure cannot influence the minimum porosity.The porosity and permeability decrease with increasing injection pressure at the same position, which explains the fact that the increase in injection pressure will increase the effect of pore sealing.Thus, in the same position, it is more difficult to inject the grout into the rock under the same constant pressure.

Effects of the Water-Cement Ratio
Figure 7 shows the influence of the water-cement ratio on the filtration.In the simulation, the water-cement ratio is varied from 0.5 to 0.8, 1.0, 2.0 and 4.0.The simulation time is 4000 s.The change in the water-cement ratio is defined as the variation in the initial concentration of injected cement grout particles in the current model.As shown in Figure 7a, the pore pressure decreases with decreasing water-cement ratio.The variation in the pore pressure is almost linear when the grout is relatively dilute.There is a notable breaking point in Figure 7b when the water-cement ratio is lower, expressing that the penetration length of the grout will decrease with decreasing water-cement ratio.The grout with lower water-cement ratio shows more severe filtration effects, which are also indicated by the results of porosities and permeabilities (see Figure 7c,d

Effects of the Grout Injection Rate
In this simulation model, the grout is injected into the rock at a constant rate.The injection rate is chosen as 0.0001, 0.0005, 0.001 , 0.005 and 0.01 m/s.When the injection rate is constant, although filtration has occurred, the pore pressure decreases almost linearly with the penetration length (Figure 8a).In addition, the pressure drops quickly with increasing injection rate.Furthermore, a small variation in the rate produces a significantly large variation in the pressure.Figure 8b shows that the concentration decreases with decreasing injection velocity, and Figure 8c,d show that the porosity and permeability increases with decreasing injection velocity.These observations show that the cement particles can move a longer distance and form a plug over a larger region when the grout velocity is higher.In addition, when the injection rate is low, the filtration needs substantially more time to occur.Thus, increasing the injection velocity can increase the influence of the filtration.

Filtration in 3D Fractured Rock
The physical and mechanical properties of the fractures are important geologic features, therein influencing the penetration and filtration of the grout.In this section, the physical parameters of the fracture are varied by changing the inclination angle (α), the dip direction (β) and the aperture.The effects of the physical properties of a fracture on the grout penetration considering filtration are tested and analyzed.Furthermore, the fracture network is generated, and the grout filtration process in a porous medium with a random distribution of fractures is simulated.When the cement particle concentration of grout is lower, the grout becomes relatively dilute and can be regarded as water.Thus, in the current model, the grout penetration range is described as the zone in which the concentration is greater than 450 kg/m 3 .

Effects of the Physical Properties of Fractures
In this simulation, the grout is injected into the rock with a constant injection pressure, and the parameters are as listed in Table 2.The simulation time is set as 4000 s.
The results of grout flow in a single fracture with different inclination angles are given in Figure 9.The fracture extends through the rock.The fracture inclination angle (α) varies from 0 • to 15 • , 30 • , 45 • , 60 • and 90 • .The dip direction of the fracture is 0 • , and the aperture is 1 mm.As shown in Figure 9, the existence of a fracture has considerable influences on the filtration process.The grout flows mainly along the fracture, and it flows quickly when approaching the fracture because the conductance of the fracture is much larger than that of the matrix.The inclination angle of the fracture changes the region of grout penetration in the domain, and the distribution of cement particles in the voids changes.In contrast, the filtration effect is inconspicuous in the fracture, as the porosity and permeability of the fracture do not change with the cement concentration or time.
Figure 10 shows the grout penetration distance in the porous medium under two scenarios: (i) considering the filtration effects and (ii) ignoring the filtration effects.Comparing Figure 10a with Figure 10b, the filtration has an obvious and remarkable influence on the grout penetration length.Filtration will stop the grout from flowing, resulting in a maximum grout penetration distance.Figure 11 shows the grout penetration in a porous medium with a single fracture in different dip directions.The dip direction of the fracture is defined to be 30     The permeability of the fracture is related to the aperture of the fracture.Figure 12 shows the influence of the fracture aperture on the grout penetration.In this simulation, the inclination angle and dip direction of the fracture are 90 • and 60 • , respectively.The aperture of the fracture is varied from 0.01 mm to 0.1 mm, 1 mm and 10 mm.In the cases where the permeability of the fracture is substantially larger than that of the rock matrix, the grout flows quickly along the fracture and completely fills the fracture.The matrix pores along the fracture wall become completely clogged by the cement particles.The fracture with a large aperture acts as a path for grout flow, whereas the larger aperture cannot influence the final grout penetration region, as shown in Figure 12a,b because, although the grout penetration into a fracture with a larger aperture can spread quickly, the filtration process in the matrix porosity stops the grout from further penetrating.When the aperture of the fracture is small, the permeability of the fracture is almost the same or less than that of the rock matrix.The velocity of the grout penetration into the fracture is equal to that into the rock matrix; therefore, the fracture cannot redirect the flow of the grout (see Figure 12c,d).Hence, a fracture with a smaller aperture has a minimal effect on the grout filtration process.

Effects of the Fracture Network
The size, dip angle, shape and distribution of fractures in a rock mass are very complicated and diverse.In the current simulation model, a fracture network using a cube model with a size of 5 × 5× 5 m is generated for exploring the grout penetration in a complicated fractured andporous rock.The aperture of the fracture is set as 0.1 mm.As shown in Figure 13a, two sets of fractures (with ten in each group) are generated randomly.The center coordinate of the fractures follows a normal distribution, and the widths and lengths of the rectangular fractures follow logarithmic distributions.The Fisher coefficient is 22.The randomly distributed parameters for the fracture geometries in this numerical model are listed in Table 3.All fractures are meshed with triangles (shown in Figure 13b), and the rock matrix is meshed with a tetrahedron (shown in Figure 13c) based on the conforming mesh techniques proposed by Wang [71].The grout is injected into the rock mass with an injection borehole located in the center of the model.The grout flow is modeled and analyzed under conditions with and without fracture networks at a constant injection pressure.The boundary conditions for the top and bottom faces are n • ∇P = 0 and n • ∇C = 0.The outlet boundaries are considered open boundaries.The pressure is defined as zero, and the concentration gradients are also defined as zero in the normal directions of these boundaries.The constant injection pressure (50 kPa) and the constant cement concentration (1000 kg/m 3 ) are set in the injection borehole.
Figures 14 and 15 show the grout flow processes influenced by filtration at four time steps.As shown in Figure 14, in the case without fractures, the grout penetration into the matrix is uniform near the injection borehole and forms a spherical grout area because the distribution of the porosity is homogeneous.The grout flow spreads outward from the borehole with time.However, the grout flow almost stops after 2000 s, as it is affected by filtration.When fractures exist in the domain, the range of the grout penetration and the grout direction change significantly, as shown in Figure 15.In the case with fractures, the grout initially penetrates along the fractures upon reaching them, and the grout spread speed is higher in the fractures than in the rock matrix.The grout fills the entire fracture and then flows along the direction of the fracture walls.The grout flow in the matrix stops when the calculation time exceeds 2000 s, while it continues to penetrate into nearby fractures (see Figure 15c,d).The fracture is less influenced by filtration.Thus, a larger pressure or more injection boreholes are needed for the grout to reach the expected point in such a domain.

Conclusions
In this study, we propose a numerical 3D UPM to simulate grout flow in fractured and porous rocks considering the effects of filtration.A two-phase model is used for the grout flow.A hyperbolic function is used to determine the relationship between the porosity and permeability that is consistent with experimental results [19].The proposed model is validated by comparing the numerically obtained results with the experimental results.Then, sensitivity analyses in consideration of the injection pressure, the water-cement ratio of the grout and the injection rate are conducted, thereby demonstrating the reliability of the proposed numerical method.Finally, the influences of the inclination angle, the dip direction, the aperture of the fracture and the complex fracture network on the filtration process are analyzed.
The simulation results indicate that filtration will stop the grout from penetrating farther into the porous and fractured medium.A high grouting pressure, a high water-cement ratio and a high grout injection rate can each increase the retention of cement particles while prolonging the penetration process.Furthermore, the existence of fractures in the domain has considerable influences on the grout flow.The inclination angle and dip direction of the fracture can redirect the flow of the grout and the grout penetration region can change because the permeability of the fracture is much larger than that of the rock matrix; thus, the grout flows quickly along the fracture and completely fills the fracture.When the grout flows into a randomly distributed fracture network, the fractures play a prominent role and influence the flow path of the grout.Although the grout stops penetrating into the rock matrix, it continues spreading into the fractures near it and extends along the fracture edges.

Figure 1 .
Figure 1.Discrete elements within fractured and porous media.

Figure 3 .
Figure 3.Comparison between the UPM model and the results obtained by Bouchelaghem [56].

Figure 5 .
Figure 5. Variations in the pore pressure, concentration, porosity and permeability with the penetration length.

Figure 6 .
Figure 6.Influences of the injection pressure on the pore pressure, concentration, porosity and permeability.
Figure7shows the influence of the water-cement ratio on the filtration.In the simulation, the water-cement ratio is varied from 0.5 to 0.8, 1.0, 2.0 and 4.0.The simulation time is 4000 s.The change in the water-cement ratio is defined as the variation in the initial concentration of injected cement grout particles in the current model.As shown in Figure7a, the pore pressure decreases with decreasing water-cement ratio.The variation in the pore pressure is almost linear when the grout is relatively dilute.There is a notable breaking point in Figure7bwhen the water-cement ratio is lower, expressing that the penetration length of the grout will decrease with decreasing water-cement ratio.The grout with lower water-cement ratio shows more severe filtration effects, which are also indicated by the results of porosities and permeabilities (see Figure7c,d)

Figure 7 .
Figure 7. Influences of the water-cement ratio on the pore pressure, concentration, porosity and permeability.

Figure 8 .
Figure 8. Influences of the injection velocity on the pore pressure, concentration, porosity and permeability.
• , 45 • , 60 • and 75 • .The inclination angle is 90 • , and the aperture of the fracture is 1 mm.The trend characterizing the grout flow is similar to that in the fractures at different inclination angles.A fracture with a different dip direction will redirect the flow of grout.

igure 9 .
Grout penetration with a single fracture at different inclination angles in the porous medium considering the filtration process.
(b) non-filtration with none fracture (a) filtration with none fracture

Figure 10 .
Figure 10.Grout penetration in a porous medium.

igure 11 .
Grout penetration with a single fracture with different dip directions in the porous medium considering the filtration process.

Figure 12 .
Figure 12.Grout penetration with a single fracture under different apertures in the porous medium considering the filtration process.

Table 3 .
Parameters adopted for the generation of the fracture models.Group Fracture Number Mean Length (m) St dev Dip Angle (Degree) Dip Direction (Degree) dev denotes the standard deviation.

( a )
Fracture network in the 3D model (b) Conforming triangulation of fractures (c) Tetrahedrization of the domain

Figure 13 .
Figure 13.Distribution of fractures in the porous, fractured rock mass.

Figure 14 .
Figure 14.Grout penetration processes at different time steps with a constant injection pressure (P = 50 kPa) in a homogeneous, porous rock considering filtration.

Figure 15 .
Figure 15.Grout penetration processes at different time steps with a constant injection pressure (P = 50 kPa) in a fractured and porous rock considering filtration.

Table 1 .
Parameters for the experimental verification simulation.

Table 2 .
Parameters for the sensitivity analysis.