Grout diffusion model in porous media considering the variation in viscosity with time

Abstract

Grouting engineering is widely used in water plugging for geotechnical engineering. However, grout is usually treated as a Newton fluid and the viscosity is considered unchangeable over time during the grouting design process. This study proposes a grout diffusion model for porous media that considers the variation in viscosity with time. The flow equation is derived for a single smooth tubule. Then, the microequation of Bingham grout flowing in porous media is obtained. Finally, an assembled diffusion model of spheres and cylinders for grouting using a perforated pipe is proposed. A numerical simulation method is used to verify the grout diffusion mode. The research results can guide grouting design and practical grouting engineering in water plugging and reinforcement.

Keywords

Time variation in viscosity Bingham grouts porous media assembled model grout diffusion

Introduction

During the construction process of tunnels and underground engineering, many disasters of water inrush have occurred in China, for example, in the Yesanguan tunnel,^1,2 the Yuanliangshan tunnel,³ and the Qiyueshan tunnel.^4–7 These disasters induce underground water and a crisis in the living environment near the project. Due to increasing viscosity with time, grouting has been widely used in water plugging and the reinforcement of tunnels.^8–10 Through grouting, the permeability of rock and soil can be dramatically decreased, and the physical and mechanical behaviors can be significantly improved. This technique has been widely applied in geotechnical engineering to prevent underground water inflow and protect the environment.

New technologies and materials for grouting have been developed to a large extent. However, the growth in grouting mechanisms has been very slow, and the related theory substantially lags behind the requirements of grouting design and construction. There are four reasons for the slow development of grouting theory. First, the injected medium is complicated; second, grouts behave with various characteristics and changeable viscosity over time. Some non-Newtonian fluids show a time-dependent change in viscosity and nonlinear stress–strain behavior, in which the longer the fluid undergoes shear stress, the lower its viscosity. Third, the interaction among grouts, injected media, and underground water has not been clearly understood. Finally, grouting engineering is a concealed project that is difficult to detect.

There are many kinds of grouting from different aspects, such as the aim, target, and grout material. According to the theory of grout diffusion, grouting can be divided into many types: seepage grouting,¹¹ split grouting,¹² filling grouting,¹³ jet grouting,¹⁴ and so on.¹⁵ Comparatively, the theory of seepage grouting is relatively mature and is represented by Maag’s spherical, Raffle’s cylindrical diffusion theory, Baker’s formula, Wallner’s formula, and G. Lombardi’s formula from which some grout diffusion formulas have been developed.^16–23 Cement materials and grout injection models are two critical parts in grout diffusion study. Some studies discussed different cement materials, such as Portland cement and silica fume, and compared the performances of various grout injection models.^24,25

Due to its low price and maturity, cement slurry has been widely applied. During practical cement grouting, perforated pipes and cement-based slurries with water-to-cement ratios between 0.8 and 1.0 have been extensively employed. In these studies, the time features and the corresponding control methods have been considered.^26–36 A starting pressure is required to initiate flow without considering the underground water pressure. Therefore, the cement-based slurry belongs to a Bingham fluid instead of a Newton fluid. Therefore, the cement slurry should be treated as a Bingham fluid during the grouting design and construction process. However, currently, we regard cement grouts as Newton fluids.

With the exception of grout properties, the diffusion rule should also be implemented. As it is known, in grouting design, spherical and cylindrical grout diffusion formulas have been used. However, based on the practical pre-grouting of tunnel faces and the results of indoor tests, the shape of cement paste is not exactly a sphere or a cylinder. We also need to study the grout diffusion model to improve the theory of grouting.

Based on the above discussion, an assembled diffusion model with Bingham grout jetted into a porous medium and a variable viscosity in time is proposed in this article. A porous medium is a medium that contains pores, and the pores are typically filled with a fluid. In addition, a skeletal medium is usually a solid with a foam-like structure. During this study, the injected medium is simplified to an isotropic and homogeneous medium, and underground water is treated as a source of pressure. The dilution and mixing of the water and grout are not considered in this study. The diffusion distance of cement grout under different grouting pressures and void ratios is analyzed. The detailed derivation process and the simulated results are described. The spherical and cylindrical diffusion equation is used to certify the results, which show that the spherical and cylindrical formula is a special example. Finally, the results of the numerical computational method also agree with the diffusion mode.

Diffusion mechanism of Bingham grouts

The diffusion mechanism of Bingham grouts is based on the rheology equation. The flow equation has been derived for a single smooth tubule. The rheological equation of Bingham fluids is represented by

τ = τ_{s} + η_{p} \cdot γ

(1)

where $τ$ and $τ_{s}$ represent the shearing stress and yield stress, respectively; $η_{p}$ is the plastic viscosity; $γ$ is shear rate ( $γ = - (dv / dr)$ ); and v and r are the flow velocity and radius, respectively.

The law for the variation in viscosity with time is expressed as

η_{p} (t) = η_{p 0} e^{kt}

(2)

where $η_{p} (t)$ and $η_{p 0}$ are the plastic viscosity at time t and the initial time, respectively.

From equations (1) and (2), the rheological equation of a Bingham fluid considering the time-dependent viscosity is obtained at a given time t, which is as follows

τ (t) = τ_{s} (t) - η_{p 0} e^{kt} \cdot \frac{dv}{dr}

(3)

At a given time t, selecting a representative elementary volume (REV; Figure 1) and using the mechanical equilibrium analysis method, equation (4) can be established

(p + dp) π r^{2} - p π r^{2} - 2 π r \cdot τ (t) \cdot dl = 0

(4)

where p, l, and r represent the pressure, the grouting radius at a given time, and the radius, respectively. Equation (4) can be simplified as follows

τ (t) = - \frac{r}{2} \frac{dp}{dl}

(5)

Substituting equation (5) into equation (3), the following equation can be obtained

\frac{dv}{dr} = \frac{e^{- kt}}{η_{p 0}} (- \frac{r}{2} \cdot \frac{dp}{dl} + τ_{s} (t))

(6)

At the same time, considering $τ = τ_{s} (t)$ and $r = r_{p}$ , equation (5) can be expressed as

r_{p} = \frac{- 2 τ_{s} (t)}{dp / dl}

(7)

Figure 1.

Flow diagram of the Bingham fluid.

A non-slip boundary condition is used at the pipe’s inner surface, namely, $v = 0$ when $r = r_{0}$ . Considering the boundary condition and integrating equation (6), the following expression can be obtained

v = \frac{e^{- kt}}{η_{p 0}} [\frac{1}{4} \cdot \frac{dp}{dl} (r_{0}^{2} - r^{2}) + τ_{s} (t) (r_{0} - r)]

(8)

Based on the flow pattern of the Bingham fluid in the thin pipe, as shown in Figure 1, the velocity of the Bingham fluid with slow migration (while considering the variation in viscosity with time) can be described as follows

{\begin{matrix} v = \frac{e^{- kt}}{η_{p 0}} [\frac{1}{4} \cdot \frac{dp}{dl} ({r_{0}}^{2} - {r_{p}}^{2}) + τ_{s} (t) (r_{0} - r_{p})] (0 < r < r_{p}) \\ v = \frac{e^{- kt}}{η_{p 0}} [\frac{1}{4} \cdot \frac{dp}{dl} ({r_{0}}^{2} - r^{2}) + τ_{s} (t) (r_{0} - r)] (r_{p} < r < r_{0}) \end{matrix}

(9)

The quantity of flow passing through the pipe’s given section per unit time is the sum of the flow in the piston section and the shear slippage section. The flow quantity per unit time can be represented as follows

Q = π r_{p}^{2} v_{p} + \int_{r_{p}}^{r_{0}} 2 π rvdr

(10)

Substituting equation (9) into equation (10), the flow quantity Q per unit time can be expressed as equation (11)

Q = - \frac{π e^{- kt}}{η_{p 0}} (\frac{1}{8} \frac{dp}{dl} r_{0}^{4} - \frac{1}{8} \frac{dp}{dl} r_{p}^{4} + \frac{1}{3} τ_{s} (t) r_{0}^{3} - \frac{1}{3} τ_{s} (t) r_{p}^{3})

(11)

Substituting equation (7) into equation (11), the flow quantity Q per unit time can be reached as follows

Q = - \frac{π e^{- kt}}{η_{p 0}} \cdot \frac{dp}{dl} \cdot r_{0}^{4} [\frac{2}{3} {(\frac{τ_{s} (t) / r_{0}}{dp / dl})}^{4} + \frac{1}{3} \frac{τ_{s} (t) / r_{0}}{dp / dl} + \frac{1}{8}]

(12)

Therefore, the average velocity at a given time t in the pipe can be expressed as

\vec{v} = \frac{Q}{π r_{0}^{2}} = - \frac{e^{- kt}}{η_{p 0}} \cdot \frac{dp}{dl} \cdot r_{0}^{2} [\frac{2}{3} {(\frac{τ_{s} (t) / r_{0}}{dp / dl})}^{4} + \frac{1}{3} \frac{τ_{s} (t) / r_{0}}{dp / dl} + \frac{1}{8}]

(13)

If the Bingham fluid can flow in the pipe, then the primary shear stress $τ_{s} (t)$ should be overcome first. Thus, the Bingham fluid initiates the pressure gradient. From the initial conditions, at the initial time, the average velocity equals zero ( $\vec{v} = 0$ ), so equation (13) can be simplified as follows

\frac{τ_{s} (t) / r_{0}}{dp / dl} = - \frac{1}{2}

(14)

Equation (14) can also be expressed as

- \frac{dp}{dl} = \frac{2 τ_{s} (t)}{r_{0}}

(15)

Defining

\frac{2 τ_{s} (t)}{r_{0}} = λ (t)

where $λ (t)$ is the gradient of the start-up pressure for the Bingham slurry with the viscosity dependent on time.

Assuming that the grouted media is a porous isotropic media, and its porosity is $ϕ$ , so the seepage speed in the grouted media can be represented as follows

\vec{V} = ϕ \vec{v} = - ϕ \frac{e^{- kt}}{η_{p 0}} \cdot \frac{dp}{dl} \cdot r_{0}^{2} [\frac{2}{3} {(\frac{τ_{s} (t) / r_{0}}{dp / dl})}^{4} + \frac{1}{3} \frac{τ_{s} (t) / r_{0}}{dp / dl} + \frac{1}{8}]

(16)

To facilitate usage, two variables, K and $β$ , are adopted as follows

\begin{matrix} K = \frac{ϕ r_{0}^{2}}{8 μ} \\ β = \frac{η_{p 0}}{μ} \end{matrix}

Thus, equation (16) can be transformed to the following format

\vec{V} = - e^{- kt} \cdot \frac{K}{β} \cdot \frac{dp}{dl} {(1 + \frac{λ (t)}{dp / dl})}^{2} \cdot {2 - \frac{4}{3} [(1 + \frac{λ (t)}{dp / dl}) - \frac{1}{4} {(1 + \frac{λ (t)}{dp / dl})}^{2}]}

(17)

Equation (17) is the equation of Bingham grouts considering the viscosity variation with time in isotropic porous media.

Combined diffusion model of cement slurry flowing in porous media

For grouting in tunnels or other underground engineering, a perforated pipe is widely applied to avoid the collapse of a grouting hole and to reinforce the surrounding rock. During the whole grouting process, the driven force of grout flow is unbalanced. The grout flow is in the same direction as the driven force, so the cement slurry should fill the pipe first and then diffuse into the voids or fractures of the injected media. Therefore, there are two stages for the whole grouting process. In the first stage, the grout fills the pipe. The diameter of the pipe is approximately 42 mm, and its volume is small. The pipe also has a small resistance because of its hollow structure. Therefore, the first stage of the grouting process is very short. In the second stage, the grout is injected from the pipe into the surround media. Compared to the second stage, the first stage is very short in time and its influences are relatively small, so the first stage is not discussed in this article.

By investigating the injecting cement paste of the grouting practice and indoor tests, when using a perforated pipe, grout diffusion occurs along the cylindrical part of the pipe and in the spherical part of the pipe at the end point. The diffusion scheme is shown in Figures 2 and 3. The blue and black arrows represent the grouting pressure and underground water pressure, respectively.

Figure 2.

Profile of the plugging.

Figure 3.

Top view of the plugging.

In this research, several assumptions are adopted. First, the grouted medium is treated as an isotropic medium; second, the grout is assumed to be an isotropic media with a variable viscosity in time; third, slurry permeates into the grouted medium after filling the pipe; fourth, the migration of slurry in the injected medium belongs to laminar conditions, which conforms to Darcy law; and finally, the effect of gravity is ignored. The theoretical permeation model is shown in Figures 2 and 3, where $p_{1}$ is the external water pressure, $l_{0}$ is the radius of the grouting pipe, l is the diffusion radius of grouting at a given time t, and m is the effective length of the grouting pipe. In this diffusion pattern, the grout quantity and the average velocity can be represented as follows

\vec{V} = \frac{Q}{A} = \frac{Q}{2 π lm + 2 π l^{2}}

(18)

Combining equations (17) and (18), equation (19) can be obtained by integration

\begin{array}{l} - e^{- k t} \cdot \frac{K}{β} \cdot \frac{d p}{d l} {(1 + \frac{λ (t)}{d p / d l})}^{2} \cdot \\ {2 - \frac{4}{3} [(1 + \frac{λ (t)}{d p / d l}) - \frac{1}{4} {(1 + \frac{λ (t)}{d p / d l})}^{2}]} = \frac{Q}{2 π l m + 2 π l^{2}} \end{array}

(19)

From equations (17)–(19), we can obtain the relationship between the grouting pressure p, the grouting volume Q, and the radius of grouting l

- e^{- kt} \cdot \frac{K}{β} \cdot \frac{dp}{dl} [\frac{1}{3} {(\frac{λ (t)}{dp / dl})}^{4} + \frac{4}{3} \frac{λ (t)}{dp / dl} + 1] = \frac{Q}{2 π lm + 2 π l^{2}}

(20)

In the beginning of the grouting period, the injection pressure is much higher than the yield pressure of the Bingham fluid. This parameter is shown as $λ (t) < < J$ , where J is the pressure gradient, $J = dp / dl$ . Thus, equation (20) can be simplified as

- e^{- kt} \cdot \frac{K}{β} \cdot \frac{dp}{dl} (\frac{4}{3} \frac{λ (t)}{dp / dl} + 1) = \frac{Q}{2 π lm + 2 π l^{2}}

(21)

From equation (21), we can obtain the grout diffusion governing equation as follows. Equation (21) includes several variables, such as the injection pressure, the diffusion radius, and the grout property

\frac{dp}{dl} = \frac{Q β e^{kt}}{(2 π lm + 2 π l^{2}) \cdot K} + \frac{4}{3} λ (t)

(22)

For the kind of slurry that does not consider time-dependent viscosity, equation (22) can be simplified as follows

\frac{dp}{dl} = \frac{Q β}{(2 π lm + 2 π l^{2}) \cdot K} + \frac{4}{3}

(23)

Considering the boundary and initial conditions, at the beginning of grouting, the time is $t_{0}$ , the pressure is $p_{0}$ and the grouting radius is $l_{0}$ . The grouting radius and the pressure are $l_{1}$ and $p_{1}$ , respectively. In addition, $t_{1}$ indicates the random time.

Considering the boundary and initial conditions, the integration of equation (23) yields

\begin{matrix} Δ p = p_{1} - p_{0} \\ = \frac{Q (t_{1}) β e^{k t_{1}}}{2 π Km} \cdot \ln (- \frac{l_{1}}{m + l_{1}}) + \frac{4}{3} λ (t_{1}) l_{1} \\ - \frac{Q (t_{0}) β e^{k t_{0}}}{2 π Km} \cdot \ln (- \frac{l_{0}}{m + l_{0}}) - \frac{4}{3} λ (t_{0}) l_{0} \end{matrix}

(24)

where $Q (t)$ is the grout flow rate from the initial time to time t.

During the whole grouting process and based on this diffusion pattern, the total grout quantity injected into the isotropic medium can be expressed as follows

Q' = \int Q (t) dt = V_{0} \cdot ϕ = (π l^{2} m + \frac{2}{3} π l^{3}) \cdot ϕ = ϕ π l^{2} (m + \frac{2}{3} l)

(25)

By simplification, at time t, equation (25) can be transformed as follows

Q (t) = \frac{ϕ π l^{2}}{t} (m + \frac{2}{3} l)

(26)

For the Bingham grouts without considering the viscosity time variation, equation (24) can be simplified as follows

Δ p = \frac{ϕ π l_{1}^{2} (m + \frac{2}{3} l_{1}) β e^{kt}}{2 π Kmt} \ln \frac{l_{1} (m + l_{0})}{l_{0} (m + l_{1})} + \frac{4}{3} λ (l_{1} - l_{0})

(27)

For Newtonian grouts without considering the viscosity time variation, equation (24) can be simplified as follows

Δ p = \frac{ϕ l_{1}^{2} (m + \frac{2}{3} l_{1}) β}{2 Kmt} \ln \frac{l_{1} (m + l_{0})}{l_{0} (m + l_{1})}

(28)

Verification of the grout diffusion model

Based on classical theory, the grout diffusion model proposed in this article is verified by a numerical simulation method.

Classical theory

1. When m tends to infinity, equation (28) can be simplified as equation (29)

Δ p = \frac{l_{1}^{2} β ϕ \ln \frac{l_{1}}{l_{0}}}{2 Kt}

(29)

The classical theory of grout column diffusion is as follows

l_{1} = \sqrt{\frac{2 Kt Δ p}{ϕ β \ln \frac{l_{1}}{l_{0}}}}

(30)

where $l_{1}$ is the radius of diffusion (cm), K is the permeability coefficient of soil (cm/s), $Δ p$ is the grouting pressure represented as the head height (cm), $l_{0}$ is the radius of the grouting pipe (cm), t is the grouting time (s), the porosity of soil is $ϕ$ , and $β$ represents the ratio of viscosity of slurry and that of water. Thus, the cylindrical diffusion equation is just a special case of this equation.

2. When m tends to $r_{0}$ , equation (28) can be simplified as equation (31)

Δ p = \ln \frac{r_{1} (r_{0} + r_{0})}{r_{0} (r_{0} + r_{1})} \cdot \frac{ϕ r_{1}^{2} β (r_{0} + \frac{2}{3} r_{1})}{2 K r_{0} t}

(31)

For the reason of $r_{1} > > r_{0}$ , equation (31) can be transformed as follows

Δ p = \frac{β ϕ l^{3}}{3 Kt l_{0}}

This equation is Maag’s formula³⁷ of a Newtonian fluid that seeps into the porous media. Therefore, Maag’s formula is also a special example of this study’s results.

Verification by numerical simulation

Geometric model

The various parameters in the geometric model are selected by repeated experiments. The length of the injected model square is 10 m. The injected pipe is located at the center. The depth and diameter of the pipe are 3 m and 10 cm, respectively, as shown in Figure 4.

Figure 4.

The geometric model.

Introduction of the simulation

Taking the coarse sand of the yellow river as the injected medium, the porosity is 39% and the bulk density is 1480 kg/m³. The pressure of the surrounding underground water is 0.2 MPa. The cement grout with a water cement ratio of 0.8 belongs to the Bingham fluid, and the density and dynamic viscosity are 1350 kg/m³ and 20 MPa·s, respectively. The environment temperature is 20°C.

The flow rate of the grout poured into the injected media by the pipe is 1 m³/s, and the grouting time is 2700 s. The entire outer boundary is treated as a continual boundary, and the grout gravity is neglected.

Computational results

To obtain the diffusion model of Bingham fluid flow in porous media, computational fluid dynamics (CFD) software is used, and some results are obtained. CFD software is a tool with amazing flexibility, accuracy, and breadth of application.

The diffusion mode is generally fit to a cylindrical–spherical shape. The three-dimensional results are shown in Figure 5. The vertical and longitudinal views are shown in Figures 6 and 7, respectively. The diffusion range of the grout changes with time, as shown in Figure 8. At the beginning, the grout radius increases quickly, and then the velocity decreases with time, but the shape of the region remains as cylindrical–semispherical. This effect is because the grout quickly extends with a small radius at the beginning. With continuous injection of the grout, the grout slowly extends. Thus, the grout radius quickly increases. Due to the same pressure on the grout from all directions, the region still maintains the same shape.

Figure 5.

3D form of grout diffusion.

Figure 6.

Vertical view of the grout pressure.

Figure 7.

Longitudinal view of the grout pressure.

Figure 8.

Grout region enlarged with time.

By analyzing the results of the numerical simulation, the velocity of grouting can be obtained. Using the integration method, the volume flow rate is reached, as shown in Figure 9. The results of the numerical simulation and theoretical formula are very similar, that is, the largest difference is 4.87%, the smallest difference is 0.56%, and the average difference is 2.03%. In addition, 93.3% of the result points exhibit a difference of less than 4%. As time develops, the variance decreases. The difference may be caused during the determination of the values of m and l from Figures 8 and 3. The correctness of the combined transport pattern of the slurry and the diffusion theoretical formula is verified, as shown in Figure 9. The numerical simulation method is used to verify the model validity, and the simulation result of the diffusion shape is similar to the cylindrical–semispherical mode proposed by this work.

Figure 9.

Change in mean volume flow rate with time.

Conclusion

This study is based on the rheology equation and considers the variation in viscosity with time. The flow equation is derived for a single smooth tubule. Then, the microequation of Bingham grouts flowing in the porous media is obtained. The main conclusions are listed as follows:

By field investigation and survey, the grout diffusion mode of cylindrical–semispherical is proposed, which fits grouting using a pipe with a many-hole uniform layout.

Based on the diffusion mode, the formula of Bingham grouts considering the viscosity changing with time has been derived. The relationship among the grouting pressure, diffusion radius, grouting time, and grouting parameters is determined.

Classical theory is used to verify the mode validity, and the results show that the classical formula of the cylinder and sphere mode is one of the special cases for this study.

The numerical simulation method is used to verify the mode validity, and the simulation result of the diffusion shape is similar to the cylindrical–semispherical mode proposed by this work.

The evaluation indicators of the results are minimal in this article. More evaluation indicators will be introduced to scientifically explain the results in the future. Further studies will provide additional contributions to grout diffusion in composite structures. Furthermore, the grout fills the pipe during the first stage of the grouting process. However, the first stage is not discussed in this article. Further investigations will consider the first stage of the grouting process.

Footnotes

Handling Editor: Gang Chen

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by the National Natural Science Foundation of China (nos 51108387, 51678495, and 51578463).

ORCID iD

Qingdong Wu

References

Bing

Risk analysis of karst disaster during constructing Yesanguan tunnel of Yiwan railway and study on its countermeasures. China Safe Sci J 2006; 16: 4–9.

Zhang

Sun

GQ.

Technology for treating burst port of filling solution cavity with high-pressure and rich water of Yesanguan tunnel on Yichang–Wanzhou railway. J Rail Eng Soc, 2010; 27: 81–86.

Liu

Wang

SR.

Study on karst waterburst mechanism and prevention countermeasures in Yuanliangshan tunnel. Rock Soil Mech 2006; 27: 228–208.

Sun

et al . Study on grouting techniques for preventing water inflows in karst tunnels. Modern Tunnel Technol 2015; 5: 178–183.

Shi

et al . Control of water inrush in typical karst tunnels in three gorges reservoir area and its application. Chinese J Rock Mech Eng 2014; 33: 1887–1896.

Sun

Zhang

Meng

et al . Primary resonance analysis and vibration suppression for the harmonically excited nonlinear suspension system using a pair of symmetric viscoelastic buffers. Nonlinear Dyn 2018; 94: 1243–1265.

Zhou

et al . An attribute synthetic evaluation system for risk assessment of floor water inrush in coal mines. Mine Water Environ 2015; 34: 288–294.

Varol

Dalgıç

Grouting applications in the Istanbul metro, turkey. Tunnel Underground Space Technol 2006; 21: 602–612.

Andjelkovic

Lazarevic

Nedovic

et al . Application of the pressure grouting in the hydraulic tunnels. Tunnel Underground Space Technol 2013; 37: 165–179.

10.

Chen

Using grouting of shield tunneling to reduce settlements of overlying tunnels: case study in Shenzhen metro construction. J Construct Eng Manage 2012; 138: 574–584.

11.

Tirupati

Experimental and numerical investigations of chemical grouting in heterogeneous porous media. PhD Thesis, University of Windsor, ON, Canada, 2005.

12.

Vliet

MJ.

Finite element modeling of fracture grouting by means of discrete fracture. Delft: Delft University of Technology, 2002.

13.

Fransson

Tsang

Rutqvist

et al . A new parameter to assess hydromechanical effects in single-hole hydraulic testing and grouting. Int J Rock Mech Min Sci 2007; 44: 1011–1021.

14.

Nikbakhtan

Ahangari

Rahmani

Estimation of jet grouting parameters in Shahriar dam, Iran. Int J Mini Sci Technol 2010; 20: 472–477.

15.

Zhang

et al . Research on mechanism of grout diffusion of dynamic grouting and plugging method in water inrush of underground engineering. Chinese J Rock Mech Eng 2011; 30: 2377–2396.

16.

Fleming

WGK

. The improvement of pile performance by base grouting. Civil Eng 1993; 97: 88–93.

17.

Kong

Sun

et al . A bi-level programming for bus lane network design. Transport Res Part C Emerg Technol 2015; 55: 310–327.

18.

Yao

et al . Transit network design based on travel time reliability. Transport Res Part C Emerg Technol 2014; 43: 233–248.

19.

Lombardi

. The role of cohesion in cement grouting of rock. In: Proceedings of the 15th international congress on large dams, Lausanne, 24–28 June 1985, pp.235–260.

20.

Yao

et al . An improved particle swarm optimization for carton heterogeneous vehicle routing problem with a collection depot. Ann Operat Res 2016; 242: 303–320.

21.

Yao

Chen

Zhang

et al . Allocation method for transit lines considering the user equilibrium for operators. Transport Res Part C Emerg Technol. Epub ahead of print 25 September 2018. DOI: 10.1016/j.trc.2018 .09.019

22.

Huang

Rong

Analysis of driving-water permeation grouting of saturated sand seam. Rock Soil Mech 2009; 30: 2016–2020.

23.

Wang

Shan

et al . Prediction of bus travel time using random forests based on near neighbors. Comput Aided Civil Infrastruct Eng 2017; 33: 333–350.

24.

Assaad

Daou

Cementitious grouts with adapted rheological properties for injection by vacuum techniques. Cement Concrete Res 2014; 59: 43–54.

25.

Assaad

Daou

Harb

Influence of thixotropy on performance of grouts placed using vacuum injection techniques. ACI Mater J 2015; 112: 189–198.

26.

Wen-Jun

Spreading model of grouting in rock mass fissures based on time-dependent behavior of viscosity of cement-based grouts. Chinese J Rock Mech Eng 2005; 24: 2709–2714.

27.

Yao

Chen

Cao

et al . Short-term traffic speed prediction for an urban corridor. Computer Aided Civil Infrastruct Eng 2017; 32: 154–169.

28.

Zhang

Wang

et al . Analysis of the displacement increment induced by removing temporary linings and corresponding countermeasures. Tunnelling and Underground Space Technology 2018; 73: 236–243.

29.

Jiang

Sheu

Peng

et al . Hinterland patterns of China Railway (CR) express in China under the belt and road initiative: a preliminary analysis. Transport Res Part E Logistic Transport Rev 2018; 119: 189–201.

30.

Zhang

Wang

Sun

et al . The sightseeing bus schedule optimization under park and ride system in tourist attractions. Ann Operat Res. Epub ahead of print 7 November 2016. DOI: 10.1007/s10479-016-2364-4

31.

Zhang

Wang

et al . Study on self-tuning tyre friction control for developing main-servo loop integrated chassis control system. IEEE Access 2017; 5: 6649–6660.

32.

Chin

Lin

JY.

Experimental validation of open-frame ROV model for virtual reality simulation and control. J Marine Sci Technol 2018; 23: 267–287.

33.

Bin

Wenxuan

Zhen

et al . A heuristic algorithm based on leaf photosynthate transport. Simulation 2018; 94: 593–607.

34.

Xiong

Zhu

Zhang

Energy recovery strategy numerical simulation for dual axle drive pure electric vehicle based on motor loss model and big data calculation. Complexity 2018; 2018: 4071743.

35.

Lin

Cheng

SC.

Block diagonal dominant remotely operated vehicle model simulation using decentralized model predictive control. Adv Mech Eng 2017; 9: 1–24.

36.

Cheng

Atmodihardjo

Woo

et al . Remote temperature monitoring device using a multiple patients-coordinator set design approach. Robomech J 2015; 2: 4.

37.

Maag

Methode zur feldmäßigen Bestimmung der Wasserdurchlässigkeit. Straße und Verkehr 1941; 19: 335–338.