Groutability of cement‐based grout with consideration of viscosity and filtration phenomenon

The groutability depends on the properties of the grout, its injection processes, and on the mechanical properties of the soil formation. During the process of pouring cement‐based grouting into a porous medium, a variation with time occurs in the viscosity of grout suspension. In addition, the particle filtration phenomenon will limit the expansion of the grouted zone because cement particles are progressively stagnant within the soil matrix. In this paper, a closed‐form solution was derived by implementing the mass balance equations and the generalized phenomenological filtration law, which can be used to evaluate the deposition of cement‐based grout in the soil matrix. The closed‐form solution relevant to a particular spherical flow was modified by a step‐wise numerical calculation, considering the variable viscosity caused by a chemical reaction, and the decrease in porosity resulting from grout particle deposition in the soil pores. A series of pilot‐scale chamber injection tests was performed to verify that the developed step‐wise numerical calculation is able to evaluate the injectable volume of grout and the deposition of grout particles. The results of the chamber injection tests concurred well with that of the step‐wise numerical calculation. Based on the filtration phenomenon, a viable approach for estimating the groutability of cement‐based grout in a porous medium was also suggested, which might facilitate a new insight in the design of the grouting process. Copyright © 2009 John Wiley & Sons, Ltd.     

Mathematical and numerical filtration–advection–dispersion model of miscible grout propagation in saturated porous media

The development of a predictive model of behaviour of porous media during injection of miscible grout, taking into account convection, dilution and filtration of grout solution with interstitial water, as well as consolidation aspects, is presented. Model assumptions are reviewed and discussed first. During the establishment of the model, we insist on surface terms and their physical relevance in expressing adsorption effects. Constitutive laws such as Fick's law for diffusive mass transport, hydrodynamic dispersion tensor dealing with miscibility, are modified by taking into account filtration effects. A new surface term appears in mass balance equations as a consequence of filtration. According to the filtration laws used, an initial filtration rate is estimated on the basis of a one‐dimensional experimental campaign. The field equations are discretized by using Galerkin finite element and θ‐scheme standard method. For transport equation, Streamline Upwind Petrov Galerkin method is employed to prevent numerical oscillations. Lastly, confrontation of numerical results with laboratory experiments constitutes a first step to validate the model on a realistic basis. Copyright © 2001 John Wiley & Sons, Ltd.     

An FEM/VOF hybrid formulation for fracture grouting modelling

A numerical model using a hybrid formulation of a finite element method (FEM) coupled with the volume of fluid (VOF) technique to simulate the fracture grouting processes in soils is described. The numerical model considered the couplings of the stress distribution, with two-phase fluid flows, and the mesh element damage. The hardening of grout in soil is described by a time-dependent Young’s modulus and viscosity. Crack initiation, branching, propagation, and grout vein growth in homogeneous and heterogeneous soils can be numerically reproduced. Although the method is developed particularly for simulating fracture grouting, the processes of compaction grouting and permeation grouting can also be numerically simulated. Some grouting cases have been simulated with results similar to the experimental results. This further confirms the adequacy and the power of the numerical approach.     

Modelling of cement suspension flow in granular porous media

A theoretical model of cement suspensions flow in granular porous media considering particle filtration is presented in this paper. Two phenomenological laws have been retained for the filtration rate and the intrinsic permeability evolution. A linear evolution with respect to the volume fraction of cement in the grout has been retained for the filtration rate. The intrinsic permeability of the porous medium is looked for in the form of a hyperbolic function of the porosity change. The model depends on two phenomenological parameters only. The equations of this model are solved analytically in the one‐dimensional case. Besides, a numerical resolution based on the finite element method is also presented. It could be implemented easily in situations where no analytical solution is available. Finally, the predictions of the model are compared to the results of a grout injection test on a long column of sand. Copyright © 2005 John Wiley & Sons, Ltd.     

Modeling of a semi-real injection test in sand

This article presents a model of flow and transport with filtration in porous media which is used to analyze large-scale grouting tests. A program based on the finite element method is developed to solve the model equations; a particular attention is paid to inherent issues of transport problems. The analysis of these tests aims at providing insights on the propagation patterns associated to the injection of cement-based grouts in sand. To apprehend particular features that characterize field injections, the experiments are performed by using a tube-a-manchette and a patented grout. Finally, the role of filtration during the tests is discussed.     

裂隙注浆偏流机理及帷幕体采动效应研究综述与展望

注浆堵水是地下工程和矿山水害防控的常用技术措施,但是,目前浆液扩散机理与采动效应的研究滞后于工程实践的要求,制约了注浆工程科学有效地实施。本文跟踪裂隙岩体中注浆浆液扩散研究,受裂隙地下水“偏流效应”启发,提出了深部矿山裂隙岩体动水注浆浆液扩散偏流效应的命题,包括2个关键科学问题：裂隙岩体动水注浆浆液扩散偏流机理、偏流形成的“浆-岩”组合体采动变形破坏和渗透性变化机理。结合研究趋势,提出了跨学科综合研究方法。以典型矿山注浆防治水工程为背景,研究深部矿山水文地质结构及其采动时空演化,揭示裂隙岩体动水条件下注浆浆液扩散偏流机理及其影响因素,获得裂隙岩体内帷幕注浆形成的三维异形“浆-岩”组合体形态和力学特性,建立浆液偏流扩散数学模型,构建“浆-岩”组合体渗透性采动效应评价方法,为深部矿山水害防治和安全生产提供科学基础。     

Numerical modelling of compensation grouting above shallow tunnels

This paper describes the development of a numerical model for compensation grouting which is a useful technique for the protection of surface structures from the potentially damaging movements arising from tunnel construction. Pipes are inserted into the ground between the tunnel and the overlaying structure from an access shaft. Buildings on the surface are instrumented and movements are carefully monitored. Once the deformations exceed a certain Trigger Level, grout is injected into the ground to prevent damage. In the finite element model described here, compensation grouting is modelled by applying an internal pressure to zero‐thickness interface elements embedded in the mesh. An ‘observational algorithm’ is used, where the deformations of the surface are monitored and used to control the injection process. Example analyses of compensation grouting are given for three‐dimensional tunnel construction underneath a greenfield site. Different strategies are used to control the injection process and their effectiveness in preventing surface movement is assessed. The numerical model is shown to replicate general behaviour expected in the field and is capable of modelling the control of ground surface movements at a greenfield site. Copyright © 2005 John Wiley & Sons, Ltd.     

Effects of medium and pore‐fluid compressibility on chemical‐dissolution front instability in fluid‐saturated porous media

Theoretical analysis and computational simulations have been carried out to investigate how medium and pore‐fluid compressibility affects the chemical‐dissolution front propagation, which is associated with a fully‐coupled nonlinear problem between porosity, pore‐fluid pressure, pore‐fluid density and reactive chemical‐species transport within a deformable and fluid‐saturated porous medium. When the fully‐coupled nonlinear system is in a subcritical state, some analytical solutions have been derived for a special case, in which the ratio of the equilibrium concentration to the solid molar density of the chemical species is approaching zero. To investigate the effect of either medium compressibility or pore‐fluid compressibility on the evolutions of chemical dissolution fronts in supercritical chemical dissolution systems, numerical algorithms and procedures have been also proposed. The related theoretical and numerical results have demonstrated that: (i) not only can pore‐fluid compressibility affect the propagating speeds of chemical dissolution fronts in both subcritical and supercritical systems, but also it can affect the growth and amplitudes of irregular chemical dissolution fronts in supercritical systems; (ii) medium compressibility may have a little influence on the propagating speeds of chemical dissolution fronts, but it can have significant effects on the growth and amplitudes of irregular chemical dissolution fronts in supercritical systems; and (iii) both medium and pore‐fluid compressibility may stabilize irregular chemical‐dissolution‐fronts in supercritical chemical dissolution systems. Copyright © 2011 John Wiley & Sons, Ltd.     

Discrete fracture model for simulating waterflooding processes under fracturing conditions

In our study, we develop a model for simulating fracturing processes in a poroelastic medium. The proposed approach combines the discrete fracture model enriched with contact plane mechanics. The model captures mechanical interactions of fractures and a deformable medium, fluid, and heat transfer in fractures and in a porous medium. Both effects of poroelasticity and thermoelasticity are accounted in our model. Mass and heat conservation equations are approximated by the finite volume method, and mechanical equilibrium equations are discretized by means of the Galerkin finite element approach. Two‐dimensional grid facets between 3‐dimensional finite elements are considered as possible fracture surfaces. Most of these facets are inactive from the beginning and are activated throughout the simulation. A fracture propagation criterion, based on Irwin's approach, is verified on each nonlinear iteration. When the criterion is satisfied, additional contact elements are added into finite element and discrete fracture model formulations respectively. The proposed approach allows modeling of existing natural and artificially created fractures within one framework. The model is tested on single‐ and multiple‐phase fluid flow examples for both isothermal and thermal conditions and verified against existing semianalytical solutions. The applicability of the approach is demonstrated on an example of practical interests where a sector model of an oil reservoir is simulated with different injection and production regimes.     

A fully coupled element‐free Galerkin model for hydro‐mechanical analysis of advancement of fluid‐driven fractures in porous media

Hydraulic fracturing (HF) of underground formations has widely been used in different fields of engineering. Despite the technological advances in techniques of in situ HF, the industry uses semi‐analytical tools to design HF treatment. This is due to the complex interaction among various mechanisms involved in this process, so that for thorough simulations of HF operations a fully coupled numerical model is required. In this study, using element‐free Galerkin (EFG) mesh‐less method, a new formulation for numerical modeling of hydraulic fracture propagation in porous media is developed. This numerical approach, which is based on the simultaneous solution of equilibrium and continuity equations, considers the hydro‐mechanical coupling between the crack and its surrounding porous medium. Therefore, the developed EFG model is capable of simulating fluid leak‐off and fluid lag phenomena. To create the discrete equation system, the Galerkin technique is applied, and the essential boundary conditions are imposed via penalty method. Then, the resultant constrained integral equations are discretized in space using EFG shape functions. For temporal discretization, a fully implicit scheme is employed. The final set of algebraic equations that forms a non‐linear equation system is solved using the direct iterative procedure. Modeling of cracks is performed on the basis of linear elastic fracture mechanics, and for this purpose, the so‐called diffraction method is employed. For verification of the model, a number of problems are solved. According to the obtained results, the developed EFG computer program can successfully be applied for simulating the complex process of hydraulic fracture propagation in porous media. Copyright © 2016 John Wiley & Sons, Ltd.     

An anisotropic damage–plasticity model for saturated quasi‐brittle materials

This article is devoted to numerical modeling of anisotropic damage and plasticity in saturated quasi‐brittle materials such as rocks and concrete. The damaged materials are represented by an isotropic poroelastic matrix containing a number of families of microcracks. Based on previous works, a discrete thermodynamic approach is proposed. Each family of microcracks exhibits frictional sliding along crack surfaces as well as crack propagation. The frictional sliding is described by a Coulomb–Mohr‐type plastic criterion by taking into account the effect of fluid pressure through a generalized effective stress concept. The damage evolution is entirely controlled by and coupled with the frictional sliding. The effective elastic properties as well as Biot's coefficients of cracked porous materials are determined as functions of induced damage. The inelastic deformation due to frictional sliding is also taken into account. The procedure for the identification of the model's parameters is presented. The proposed model is finally applied to study both mechanical and poromechanical responses of a typical porous brittle rock in drained and undrained compression tests as well as in interstitial pressure controlled tests. The main features of material behaviors are well reproduced by the model. Copyright © 2012 John Wiley & Sons, Ltd.     

Theoretical and numerical analyses of chemical‐dissolution front instability in fluid‐saturated porous rocks

The chemical‐dissolution front propagation problem exists ubiquitously in many scientific and engineering fields. To solve this problem, it is necessary to deal with a coupled system between porosity, pore‐fluid pressure and reactive chemical‐species transport in fluid‐saturated porous media. Because there was confusion between the average linear velocity and the Darcy velocity in the previous study, the governing equations and related solutions of the problem are re‐derived to correct this confusion in this paper. Owing to the morphological instability of a chemical‐dissolution front, a numerical procedure, which is a combination of the finite element and finite difference methods, is also proposed to solve this problem. In order to verify the proposed numerical procedure, a set of analytical solutions has been derived for a benchmark problem under a special condition where the ratio of the equilibrium concentration to the solid molar density of the concerned chemical species is very small. Not only can the derived analytical solutions be used to verify any numerical method before it is used to solve this kind of chemical‐dissolution front propagation problem but they can also be used to understand the fundamental mechanisms behind the morphological instability of a chemical‐dissolution front during its propagation within fluid‐saturated porous media. The related numerical examples have demonstrated the usefulness and applicability of the proposed numerical procedure for dealing with the chemical‐dissolution front instability problem within a fluid‐saturated porous medium. Copyright © 2007 John Wiley & Sons, Ltd.     

帷幕灌浆扩散半径及数值模拟的研究

从连续性方程出发,考虑浆液压力对土体孔隙率的影响,分析推导了浆液在多孔介质中的渗流规律,并给出了扩散半径的简单近似计算公式;然后将二相流理论应用到注浆研究中,假设孔隙由水和浆液完全充填,且二者不相混溶,建立了浆液驱水的非稳定渗流模型。根据某大坝的实际情况,用FLAC软件中的二相流模块对大坝灌浆过程进行了模拟分析,并与推导公式相比较。由于二相流理论考虑了毛细压力的作用,模拟得到的浆液扩散速率递减得更快。模拟结果表明,浆液和水之间存在着一个过渡带,浆液的饱和度在不同时间和位置上是变化的。随着时间的延长,浆液扩散得越远,但其速率逐渐减小。扩散半径不仅与渗透系数有关,还与孔隙度有关,而且孔隙度较渗透系数对扩散半径有着更大的影响。二相流理论可以更好地模拟帷幕灌浆的浆液扩散情况。     

Surface waves in a cylindrical borehole through partially-saturated porous solid

Propagation of surface waves is discussed in a cylindrical borehole through a liquid-saturated porous solid of infinite extent. The porous medium is assumed to be a continuum consisting of a solid skeletal with connected void space occupied by a mixture of two immiscible inviscid fluids. This model also represents the partial saturation when liquid fills only a part of the pore space and gas bubbles span the remaining void space. In this isotropic medium, potential functions identify the existence of three dilatational waves coupled with a shear wave. For propagation of plane harmonic waves along the axially-symmetric borehole, these potentials decay into the porous medium. Boundary conditions are chosen to disallow the discharge of liquid into the borehole through its impervious porous walls. A dispersion equation is derived for the propagation of surface waves along the curved walls of no-liquid (all gas) borehole. A numerical example is studied to explore the existence of cylindrical waves in a particular model of the porous sandstone. True surface waves do not propagate along the walls of borehole when the supporting medium is partially saturated. Such waves propagate only beyond a certain frequency when the medium is fully-saturated porous or an elastic one. Dispersion in the velocity of pseudo surface waves is analysed through the changes in consolidation, saturation degree, capillary pressure or porosity.     

Analytical solutions for dynamic pressures of coupling fluid–porous medium–solid due to SV wave incidence

This paper presents the results of theoretical investigation on the dynamic coupling of an ideal fluid‐porous medium‐elastic half‐space system subjected to SV waves to study the effect of sediment on the seismic response of dams for reservoirs that are deposited with a significant amount of sediment after a long period of operation. The effects of the porous medium and the incident wave angle on dynamic pressures in the overlying ideal fluid are analyzed, and the reflection and transmission coefficients of the wave at the material interfaces are derived using an analytical solution in terms of displacement potentials. The numerical test of modeling shows that the dynamic pressures significantly depend on the properties of porous medium. The fully saturated porous medium reduces the response peaks slightly, while the partially saturated porous medium causes a considerable increase in the resonance peaks. Copyright © 2009 John Wiley & Sons, Ltd.     

Hydro‐mechanical modeling of cohesive crack propagation in multiphase porous media using the extended finite element method

In this paper, a numerical model is developed for the fully coupled hydro‐mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non‐wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two‐phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid skeleton is accounted. The cohesive crack model is integrated into the numerical modeling by means of which the nonlinear fracture processes occurring along the fracture process zone are simulated. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three‐phase formulation. The other variables are incorporated into the model via the experimentally determined functions, which specify the relationship between the hydraulic properties of the fracturing porous medium, that is saturation, permeability and capillary pressure. The spatial discretization is implemented by employing the extended finite element method, and the time domain discretization is performed using the generalized Newmark scheme to derive the final system of fully coupled nonlinear equations of the hydro‐mechanical problem. It is illustrated that by allowing for the interaction between various processes, that is the solid skeleton deformation, the wetting and the non‐wetting pore fluid flow and the cohesive crack propagation, the effect of the presence of the geomechanical discontinuity can be completely captured. Copyright © 2012 John Wiley & Sons, Ltd.     

考虑扩散路径的宾汉姆流体渗透注浆机制

浆液在多孔介质中的扩散路径对渗透扩散范围和注浆效果具有非常重要的影响。采用理论分析,以分形特征与宾汉姆流体在多孔介质中的渗流运动方程为基础,揭示了考虑扩散路径的宾汉姆流体渗透注浆机制,并利用团队前期开展的渗透注浆试验对其进行了验证。分析了多孔介质孔隙率、宾汉姆水泥浆液水灰比、多孔介质渗透系数、注浆压力、地下水压力等对扩散半径的影响变化规律。同时,基于Comsol Multiphysics平台,采用计算机编程技术二次开发得到了考虑扩散路径的宾汉姆流体渗透注浆机制的渗透注浆三维数值模拟程序,并以此开展了宾汉姆水泥浆液在多孔介质中渗透扩散形态效果的数值模拟。研究结果表明：与不考虑扩散路径的宾汉姆流体渗透注浆球形扩散公式获得的扩散半径理论计算值相比,采用考虑扩散路径的宾汉姆流体渗透注浆机制得到的扩散半径理论计算值更接近试验值。该研究成果可为实践注浆工程提供一定的理论支撑。     

基于流体体积法的劈裂注浆有限元分析

劈裂注浆作为一种有效的土体加固方法,其理论研究落后于工程实践。提出了基于弥散裂缝模型和流体体积法的劈裂注浆有限元分析方法,并通过ABAQUS二次开发编写劈裂注浆有限元程序。数值分析结果能与室内试验较好吻合,验证了有限元算法的合理性。在此基础上研究了注浆孔埋深和注浆流量对劈裂浆脉形状的影响。结果表明,劈裂注浆过程可分为起劈和劈裂发展两个阶段。当劈裂浆脉扩展到模型边界后,继续注浆会引起注浆压力的大幅提高和已有浆脉宽度的增加。随着注浆孔埋深的增加,浆脉分支减少,长度减小,宽度增加,劈裂浆脉形状的主要控制因素从土体参数的随机性变成大小主应力值的差异。注浆量一定的情况下,注浆速率越大,劈裂浆脉长度越短,宽度越大,注浆终压也越大。研究为劈裂注浆的工程应用提供理论支撑。     

单一裂隙动水注浆扩散模型

通过理论推导建立了单一裂隙动水注浆扩散模型,用计算机对单一裂隙动水注浆扩散模型进行编程和分析,试验验证该注浆扩散模型是合理的。研究结果表明,单一裂隙动水注浆在注浆前期浆液扩散面积呈圆形而在注浆后期浆液扩散面积呈椭圆形;浆液扩散半径与浆液黏度成反比,与裂隙开度和注浆压力成正比;水流流动有利于浆液沿顺水流方向扩散而抑制浆液沿逆水流方向扩散;水流流速增大或减小到一定程度时浆液扩散半径受水流流速的影响程度将逐渐减小;水流流动对逆水流方向扩散半径的影响程度大于顺水流方向扩散半径的影响程度;浆液扩散半径受水流流速影响最大,受裂隙开度影响次之,受注浆压力和浆液黏度的影响相对较小