ORIGINAL RESEARCH article

Front. Earth Sci., 30 January 2023

Sec. Environmental Informatics and Remote Sensing

Volume 10 - 2022 | https://doi.org/10.3389/feart.2022.1100196

Permeation grouting mechanism of viscous time-varying fluid considering diffusion path

  • 1. College of Civil Engineering and Architecture, Shandong University of Science and Technology, Qingdao, China

  • 2. Geotechnical and Structural Engineering Research Center, Shandong University, Jinan, China

  • 3. School of Transportation and Civil Engineering, Shandong Jiaotong University, Jinan, China

  • 4. Qingdao West Coast Rail Transit Co., Ltd., Qingdao, China

Abstract

Due to the hidden nature of the grouting project, it has not been possible to get a clear picture of the slurry transport law, and most of the existing permeation grouting theories have ignored the tortuous characteristics of the pore channels of porous media. Based on the fractal theory of porous media and the equation of seepage motion of viscous time-varying Bingham fluid, a theoretical model of permeation grouting with viscous time-varying Bingham fluid is established, and viscosity tests are carried out for C-S slurry to obtain the time-varying equation of viscosity of C-S slurry. A set of permeation grouting simulation experiment device was designed, and the permeation grouting simulation experiment of C-S slurry was carried out to obtain the variation curve of grouting pressure with time under different injection media and different grouting rates. The results show that the grouting pressure considering the diffusion path of the slurry is 0.9–1.1 times of the experimental results, the grouting pressure obtained from the experimental is 0.7–1.2 times of the theoretical value considering the diffusion path, and the calculated value considering the diffusion path agrees well with the experimental value. The theoretical model can provide scientific guidance for field grouting construction.

1 Introduction

Grouting is to inject some solidified slurry into the cracks or pores of rock and soil mass through grouting equipment, and to achieve the purpose of seepage prevention and leakage stoppage of rock and soil mass through permeation, filling, compaction and splitting. Among them, permeation grouting is widely used in porous medium strata, such as medium, coarse sand layers and gravel layers, etc. These strata have large pore sizes, and slurry mostly diffuses in the form of permeation (). Because the gel time of quick-setting slurry is controllable, it ranges from tens of seconds to several minutes (), the gel time can be effectively controlled, and the quick-setting slurry can be used to quickly and effectively block and reinforce the sand layer in view of the water inrush and other disasters in the construction process. At present, there is a lack of scientific guidance for the application of penetration grouting of quick-setting slurry, and the relevant grouting parameters are still determined by engineering experience.

At present, in the theoretical research of permeation grouting in sand layer, Maag () derives the theoretical formula of permeation grouting with Newtonian fluid based on Darcy’s law and the initial viscosity of slurry; Ma () gave the theoretical formula of spherical diffusion of Newtonian fluid; Yang () considered four factors: grouting pressure, grouting time, permeability coefficient and water-cement ratio, and studied the permeation diffusion mechanism of slurry in sand layer through grouting simulation test. Yang () derived the calculation formulas of cylindrical and spherical diffusion radius based on the rheological equation of Bingham fluid; Yang (; ) established theoretical models of Bingham fluid and power-law fluid permeation grouting respectively; Yang () established a theoretical model of cylindrical diffusion of power-law fluid based on the rheological equation of power-law fluid. With the development of permeation grouting theory, the viscosity parameters of viscous time-varying fluid were gradually introduced, and Li () established the diffusion model of C-S slurry in a single slab crack; Yang () derived the formula for calculating the permeation diffusion radius of Bingham fluid based on the equation of time-dependent viscosity of slurry; Zhou () obtained the calculation formula of diffusion radius of viscous time-varying fluid through theoretical derivation; Based on the rheological equation of Bingham fluid, Zhang () established the theoretical model of permeation grouting diffusion of quick-setting slurry; Ye (; ) and Liu () studied the permeation diffusion mechanism of viscous time-varying slurry behind the shield wall. In addition, some scholars have studied the theory of permeation grouting diffusion from the aspect of percolation filtration effect, such as Axelsson M (), Jong-Sun Kim () and so on. Feng () and Li () have studied the influence of percolation filtration effect on slurry diffusion process.

To sum up, at present, some progress has been made in the research of permeation grouting theory, but most of the above theoretical studies have not considered the tortuous effect of the pore channel of the injected sand body. Because of the pore structure characteristics of the injected sand body itself, the diffusion path of the slurry is not a simple straight line, and it often migrates and spreads along the channels connected with the pores of the sand layer. Therefore, the above theory has low applicability to the permeation grouting engineering of the sand layer. Zhou (), Chen (), Yang () and Zhang () discussed the theoretical models of the diffusion path of Newton, Bingham and power-law slurries in porous medium respectively, but the research on the migration and diffusion law of viscous time-varying slurries in porous medium was not deep enough. In this paper, based on fractal theory of porous medium, a theoretical model of permeation grouting with viscous time-varying fluid is established.

2 Theoretical model of Bingham fluid permeation grouting considering diffusion path of porous medium

2.1 Basic assumptions of grouting theory

To study the permeation diffusion mechanism of Bingham fluid in injected porous medium, the following assumptions were made:

  • 1. The slurry is isotropic, homogeneous and incompressible Bingham fluid.

  • 2. The injected porous medium is homogeneous and isotropic, which meets the conditions of injectability.

  • 3. The grouting rate is constant, the seepage of slurry is laminar flow, and the flow pattern of slurry remains unchanged during the seepage process.

  • 4. The influence of gravity is not considered in the grouting process.

2.2 Fractal theory of porous medium

Due to the fractal characteristics of porous medium, the seepage path of slurry in porous medium often presents a “tortuosity” effect. Figure 1 shows a schematic diagram of the seepage diffusion path of slurry in porous medium, l denotes the actual route length of the seepage channel and l0 denotes the straight line length of the seepage channel. It can be seen from Figure 1 that the seepage path of slurry in porous medium is not all straight lines, but flows along the curved channels connected by pores, showing a typical “tortuosity” effect.

FIGURE 1

The seepage channel in porous medium conforms to the fractal scale law, and the fractal form of seepage channel () can be specifically expressed as follows

Derivation of Eq. 1 gives:Where: Dt is the fractal dimension of tortuosity of pore channel; b is the radius of the pore channel.Tortuous is usually defined as ()

The relationship between average tortuosity and porosity φ can be expressed as ()

The tortuosity fractal dimension of pore channels in porous medium () can be expressed asWhere: is the average tortuosity of pore channels of porous medium; is the average radius of pore channels in porous medium.

According to references (), the average radius of pore channels in porous medium can be expressed as

The relationship between and porosity φ of porous medium can be expressed by the following formula:

The maximum radius of porous medium can be expressed by the following formula:Where: Df is the fractal dimension of pore channels; bmin and bmax represent the minimum radius and maximum radius of porous medium pore channels, respectively; d+ is generally taken as 24; r is the average radius of particles in porous medium.

According to fractal theory, the fractal dimension of pores can be expressed as ()

Reference () gives the structural parameters of porous medium:

The expressions for the fractal dimension of the tortuosity of the pore channels in porous media can be obtained by combining (Eqs 410).

2.3 Temporal and spatial distribution equation of grouting pressure

According to advance mech mechanics of fluids in porous medium (), the average velocity of Bingham fluid in a single seepage channel isWhere: is the average velocity in the seepage channel; b is the radius of seepage channel; is the pressure gradient in the direction of slurry seepage; τ0 is the initial yield stress of the slurry.

The seepage velocity v at any point in the injected porous medium, the average velocity of seepage channel at that point and the porosity of the porous medium satisfy ():

Radius b of seepage channel of injected porous medium and permeability k satisfy ():

Under the condition of constant grouting rate, the grouting flow rate meets

The pressure gradient along the diffusion direction of the Bingham fluid is obtained by combining (Eqs 1215)

Substituting Eqs 2, 14 into Eq. 16 to obtainWhen there is boundary condition , , are groundwater pressure.

Grouting pressure at grouting port position can be expressed asWhere: q is the slurry flow rate; μ is the viscosity of slurry; S is the cross-sectional area of slurry diffusion; k is the permeability of the injected medium; τ0 is the initial yield stress of slurry; φ is the porosity of the injected medium; Dt is the tortuosity fractal dimension of seepage channel; p0 is the groundwater pressure; pc is the pressure at the grouting mouth.

3 Permeation grouting simulation experiment of sand body

3.1 Experiment device

In order to study the penetration diffusion mechanism of C-S slurry in sand body, a set of visual penetration grouting simulation experiment device has been designed, as shown in Figure 2. The experiment device consists of transparent plexiglass tube, pressure sensor, flow sensor, double-liquid mixer and manual grouting pump.

FIGURE 2

3.1.1 Permeation grouting tube

In order to obtain the dynamic process of slurry diffusion clearly and intuitively, the permeation grouting tube is made of PMMA (plexiglass), which can withstand the grouting pressure of 10 MPa. The permeation grouting tube is arranged vertically with an inner diameter of 10 cm and a length of 2 m. The inner part of the transparent plexiglas tube is filled with river sand as the injected medium with fixing device at both ends. The fixing device only allows the slurry to pass through and does not allow the injected medium to pass through, so as to prevent the whole movement of river sand during grouting.

3.1.2 Data recorder

In order to monitor the grouting rate and grouting pressure in the process of grouting in real-time, flow sensor and pressure sensor are set on the grouting tubeline. The pressure sensor has a measuring range of 0–10 MPa and a measuring accuracy of 0.02 MPa. The flow sensor has a measuring range of 0–30 L/min and a measuring accuracy of 0.3 L/min.

3.1.3 Manual grouting pump

The cement slurry and silicate slurry used in the experiment are pumped by manual grouting pump. The manual grouting pump can control the grouting pressure at 0–10 MPa and the pumping flow range is 0–10 L/min.

3.2 Viscosity time-dependent test of C-S slurry

The C-S slurry is a two-liquid slurry material made of a mixture of cement slurry and silicate slurry, the gel time of the slurry can be controlled within a few tens of seconds to a few minutes. The cement used for the test was 32.5R ordinary portland cement, the quality of the cement complied with the standard “Common Portland Cement” GB175-2007, the silicate modulus M=3.0, the concentration Be`=38, the water/cement ratio W/C is 1:1 and the mixing volume ratio is 1:1.

SV-10/SV-100 sine wave vibration viscometer (see Figure 3) is selected as the test instrument with a measuring range of 0.3–10,000 mPa s/1–100 Pa s. The viscosity and temperature curve of C-S slurry can be measured directly with time.

FIGURE 3

The curves of the viscosity and temperature of the C-S slurry with time are shown in Figure 4.

FIGURE 4

It can be seen from Figure 4 that although the temperature of slurry increases after mixing, the change range is within 1°C, so the slurry temperature can be considered to be basically unchanged. The initial viscosity growth of C-S slurry is small. With the reaction time, the slurry viscosity increases rapidly. The slurry viscosity curve can be fitted by power function form :Where: is the apparent viscosity of C-S slurry; t is the time after mixing the C-S slurry.

3.3 Experiment plan

Three kinds of river sand with different particle size gradations are selected as the injected medium and the grouting experiment environment is anhydrous. The permeability coefficient and porosity of the injected medium are determined according to the constant head permeability test and “Standard for geotechnical testing method” ().

The physical and mechanical parameters of the three types of river sands injected are shown in Table 1.

TABLE 1

Number of injected mediaParticle size range/mmPorosity φ (%)Permeability k/m2
10.5–136.42.218 × 10−9
21–237.45.857 × 10−9
32–538.51.512 × 10−8

The physical and mechanical parameters of injected medium.

During this experiment the grouting rate was kept constant and designed at 1, 2 and 3 L/min. 5 sets of experiment conditions were designed according to the porosity of the injected media and the grouting rate, as shown in Table 2.

TABLE 2

Condition numberGrouting rate/(L/min)Porosity of injected medium φ (%)Grouting time/s
1137.460
2237.460
3337.450
4238.560
5236.450

Experiment conditions.

It is inevitable that there will be certain fluctuation of grouting rate when using manual grouting pump in the experiment. Therefore, it is considered that the fluctuation of grouting rate is constant within 10%.

4 Analysis of experimental results

With the growth of time, the slurry pressure is increasing, and the growth rate of slurry pressure also has a significant increase. The analysis suggests that the C-S slurry is at a low viscosity stage due to the low viscosity of the slurry mix at the beginning, the slurry viscosity grows slowly and the slurry flow diffusion resistance is at a low range; as time grows, the C-S slurry viscosity increases rapidly and the slurry flow diffusion resistance rises rapidly, making the grouting pressure increase as well. The variation curve of grouting pressure is shown in Figure 5.

FIGURE 5

The grouting pressure is positively related to the grouting rate, while the grouting pressure is negatively related to the porosity of the injected medium.

5 Comparison and analysis of theoretical model and experimental results

During the experiment the diffusion front on the wall of the Plexiglas tube and the diffusion front in the river sand are basically consistent, the slurry is diffused in a one-dimensional form, the slurry diffusion distance can be directly obtained by measuring the position of the slurry front on the wall of the tube.

For one-dimensional diffusion penetration grouting, the diffusion distance and grouting time of grout meet:

The relationship between grouting pressure and diffusion distance in one-dimensional permeation diffusion process can be obtained by the simultaneous Eqs 19, 21.

Simultaneous (Eqs 16, 21) obtainable:Whereas Eq. 22 is the temporal and spatial distribution equation of grouting pressure considering diffusion path and Eq. 23 is the temporal and spatial distribution equation of grouting pressure without considering seepage path.

Again, by substituting into Eqs 22, 23 the equation for grouting pressure versus time can be obtained, without going into too much detail here. The relevant calculation parameters for the theoretical model are shown in Table 3.

TABLE 3

Condition numberq (L/min)k (10−9m2)t (s)S (m2)P0 (MPa)
115.8570.374600.007850
225.8570.37460
335.8570.37450
4215.120.38560
522.2180.37450

Calculation parameters of theoretical model.

5.1 Analysis of the variation of grouting pressure with time

By substituting the relevant parameters into the variation of grouting pressure versus time equation, the variation of grouting pressure versus time curve can be obtained, and the specific calculation results are compared with the experimental data as shown in Figure 6.

FIGURE 6

As can be seen from Figure 6, at the beginning of the grouting period, the grouting pressure is small and the growth rate is also small; as the grouting proceeds, there is a significant increase in the growth rate of the grouting pressure. There is a certain difference between the calculated results considering the slurry diffusion path compared to the model without considering the slurry diffusion path, and the theoretical difference increases with time. When the final design grouting time is reached, the calculated results considering the slurry diffusion path are 2.9–3.6 times higher than those without considering the slurry diffusion path, and the grouting pressure considering the slurry diffusion path is 0.9–1.1 times higher than the experimental results, and the experimental values are in good agreement with the theoretical.

5.2 Analysis of relationship between grouting pressure and diffusion distance

The calculated parameters are substituted into Eqs 22, 23 to obtain the variation curve of grouting pressure versus diffusion distance, and the specific calculated results are compared with the test data as shown in Figure 7.

FIGURE 7

As can be seen from Figure 7, the grouting pressure has a significant effect on the slurry diffusion distance when the grouting pressure is at a low level; when the grouting pressure is at a high level, the slurry diffusion distance is significantly reduced by the grouting pressure. The theoretical model considering the diffusion path of slurry has a faster growth rate, while the theoretical model without the diffusion path has a slower growth rate, and the difference in slurry pressure increases as the diffusion distance increases. When the same diffusion distance as the test value is reached, the grouting pressure value obtained from the experimental is 2.4–4.1 times the value without considering the diffusion path, and the grouting pressure value obtained from the experimental is 0.7–1.2 times the theoretical value with considering the diffusion path, which is within the acceptable error range.

In summary, during penetration grouting, the grouting pressure values obtained from the experiment differed significantly from the calculated values without considering the diffusion path, while the calculated values considering the diffusion path agreed better with the experimental values, so it was necessary to consider the slurry diffusion path.

6 Conclusion

  • (1) Based on the fractal theory of porous media and the equation of seepage motion of viscous time-varying Bingham fluid, a theoretical model of permeation grouting with viscous time-varying Bingham fluid is established.

  • (2) The viscosity of the C-S slurry was tested and the variation curve of the viscosity of the C-S slurry with time was obtained. A set of permeation grouting simulation experiment device was designed, and the permeation grouting simulation experiment of C-S slurry was carried out to obtain the variation curve of grouting pressure with time under different injection media and different grouting rates.

  • (3) When the designed grouting time is reached, the calculation results considering the diffusion path are 2.9–3.6 times the results without the diffusion path, and the grouting pressure considering the slurry diffusion path is 0.9–1.1 times the experimental results.

  • (4) For the same diffusion distance conditions as the experimental, the grouting pressure values obtained from the experimental are 2.4–4.1 times the values obtained without considering the diffusion path, and the grouting pressure obtained from the experimental is 0.7–1.2 times the theoretical value considering the diffusion path.

Statements

Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author contributions

HW is mainly responsible for the overall writing of the full text. LD and QZ are responsible for data processing and analysis. ZL is responsible for the innovative ideas of the article. PZ are responsible for the test operation. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 52109131); the Natural Science Foundation of Shandong Province (Grant No. ZR2020QE290); and the National Natural Science Foundation of China (Grant No. 52179120 and 52171267).

Conflict of interest

Author PZ was employed by Qingdao West Coast Rail Transit Co., Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Correction note

This article has been corrected with minor changes. These changes do not impact the scientific content of the article.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Summary

Keywords

permeation grouting, fractal theory, porous medium, diffusion path, viscous time-varying Bingham fluid, C-S slurry

Citation

Wang H, Dong L, Zhang Q, Li Z and Zhang P (2023) Permeation grouting mechanism of viscous time-varying fluid considering diffusion path. Front. Earth Sci. 10:1100196. doi: 10.3389/feart.2022.1100196

Received

16 November 2022

Accepted

12 December 2022

Published

30 January 2023

Corrected

24 July 2026

Volume

10 - 2022

Edited by

Shengli Yang, China University of Mining and Technology, Beijing, China

Reviewed by

Qin Yan, Taiyuan University of Technology, China

Tianqi Jiang, Kyushu University, Japan

Updates

Copyright

*Correspondence: Zhipeng Li,

This article was submitted to Environmental Informatics and Remote Sensing, a section of the journal Frontiers in Earth Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics