Hydro-mechanical Coupling Effect Research Articles

Excavation-induced hydraulic conductivity reduction around a tunnel – Part 1: Guideline for estimate of ground water inflow rate

Field observations indicate that current engineering practice does not consistently estimate ground water flows into unlined rock excavated tunnels due to various factors that analytical solutions do not take into account. These factors include significant geological features, groundwater drawdown, inadequate estimates of hydraulic conductivity from packer tests, and stress-induced rock-mass permeability reduction in the vicinity of tunnel (lining-like zone). A key variable that is not properly accommodated in current practice, is the hydro-mechanical interaction within the joints in the surrounding rock mass. The significance of this variable is discussed in the 1st part of the paper which presents an analytical solution assessing ground water inflow rate into a tunnel using a mathematical derivation that takes into account the excavation-induced rock permeability reduction in the vicinity of a tunnel based on hydro-mechanical coupling effect. In the 2nd part of the paper, results from numerical analysis are presented to verify the proposed analytical solution for estimating ground water inflow rate into a tunnel. Further studies are currently underway to identify other key variables and their impact on the behavior of unlined tunnels and hydrological flow regime in the surrounding fractured rock mass using a distinct element method program which can fully consider hydro-mechanical coupled behavior of joints.

Observations and predictions of hydromechanical coupling effects in the Boom clay, Mol Underground Research Laboratory, Belgium

The Boom clay is considered as a candidate host rock for the Belgian disposal of radioactive waste. Thanks to the recent sinking of a new shaft to extend the Underground Research Laboratory (URL) at the Mol site, new in-situ data have been obtained. They consist in the observation of significant fracturing, evidenced during the excavation of the shaft annexes, and in correlated hydraulic perturbations in the far field (maximum of 0.2 MPa at 60-m radius). These observations cannot be explained from simple poroplastic models that predict a maximum perturbation extent of 25 m under full deconfinement. Thus, models with increasing complexity are built to explain this discrepancy. It is proposed that the low-support pressure imposed by the primary lining, combined with the long time over which this support condition held, has favoured the decompression of the clay massif through time-dependent effects, mainly skeleton viscosity. This resulted in near-field fracturing, leading to an apparent increase of the excavated radius and to a hydraulic perturbation in the far field. This interpretation still needs to be validated/invalidated with future research. The effect on the overall performance of the repository is expected to remain small due to the self-healing nature of the fractures.

Linear coupled analysis of desiccation shrinkage in a double‐layer partially saturated medium: semi‐explicit solutions

AbstractSolutions are presented for the problem of isothermal dessiccation shrinkage in a double‐layer porous partially saturated medium. The rheological model taken into account is linear poroelastic. Hence the analysis is mainly focused on hydromechanical coupling effects and contrasts of mechanical and hydraulic properties between two materials: a low thickness skin comprised between the outer boundary and the reference porous material. Three one‐dimensional ideal structures are taken into account: a wall of finite thickness (cartesian geometry), a thick cylinder and a thick sphere. The solution of the time‐dependent problem is arrived at by applying Laplace transforms to the field variables. Exact solutions are obtained in Laplace transform space using Mathematica© to solve the field equations whilst taking into account the continuity equations at the interface and the boundary conditions. The Talbot's modified algorithm has been performed to invert the Laplace transform solutions. A bibliographical and numerical study shows that this method is remarkably precise, stable and close to the analytical inversion. Results are presented using poroelastic data representative of a concrete material and involve a strong coupling effect between hydraulical and mechanical behaviours. A first approach elastic modelling of degradation process have been presented using a thin outer layer. Apart from emphasising the semi‐explicit solution utility due to accurate speed calculation, this paper deals with more complex problems than those which can be solved using purely analytical solutions. Copyright © 2001 John Wiley & Sons, Ltd.