• Tunnel and Underground Space
• /
• v.29 no.3
• /
• pp.197-213
• /
• 2019

We simulated the fault reactivation experiment conducted at 'Main Fault' intersecting the low permeability clay formations of Mont Terri Underground Research Laboratory in Switzerland using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. We formulate the hydro-mechanical coupling relation of hydraulic aperture to consider the elastic fracture opening and failure-induced dilation for reproducing the abrupt changes in injection flow rate and monitoring pressure at fracture opening pressure. A parametric study was conducted to examine the effects of in-situ stress condition and fault deformation and strength parameters and to find the optimal parameter set to reproduce the field observations. In the best matching simulation, the fracture opening pressure and variations of injection flow rate and monitoring pressure showed good agreement with field experiment results, which suggests the capability of the numerical model to reasonably capture the fracture opening and propagation process. The model overestimated the fault displacement in shear direction and the range of reactivated zone, which was attributed to the progressive shear failures along the fault at high injection pressure. In the field experiment results, however, fracture tensile opening seems the dominant mechanism affecting the hydraulic aperture increase.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.359-381
• /
• 2020

Within the framework of DECOVALEX-2019 Task D, full-scale engineered barriers experiment (FEBEX) at Grimsel Test Site was numerically simulated to investigate an applicability of implemented Barcelona basic model (BBM) into TOUGH2-MP/FLAC3D simulator, which was developed for the prediction of the coupled thermo-hydro-mechanical behavior of bentonite buffer. And the calculated heater power, temperature, relative humidity, total stress, saturation, water content and dry density were compared with in situ data monitored in the various sections. In general, the calculated heater power and temperature provided a fairly good agreement with experimental observations, however, the difference between power of heater #1 and that of heater #2 could not captured in the numerical analysis. It is necessary to consider lamprophyre with low thermal conductivity around heater #1 and non-simplified installation progresses of bentonite blocks in the tunnel for better modeling results. The evolutions and distributions of relative humidity were well reproduced, but hydraulic model needs to be modified because the re-saturation process was relatively fast near the heaters. In case of stress evolutions due to the thermal and hydraulic expansions, the computed stress was in good agreement with the data. But, the stress is slightly higher than the measured in situ data at the early stage of the operation, because gap between rock mass and bentonite blocks have not been considered in the numerical simulations. The calculated distribution of saturation, water content, and dry density along the radial distance showed good agreement with the observations after the first and final dismantling. The calculated dry density near the center of the FEBEX tunnel and heaters were overestimated compared with the observations. As a result, the saturation and water content were underestimated with the measurements. Therefore, numerical model of permeability is needed to modify for the production of better numerical results. It will be possible to produce the better analysis results and more realistically predict the coupled THM behavior in the bentonite blocks by performing the additional studies and modifying the numerical model based on the results of this study.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.11
• /
• pp.69-78
• /
• 2012

Conventional numerical analysis for rainfall-induced slope stability has been estimated by separating seepage and stress-strain behavior, respectively. Many researchers' models from commercial softwares and literatures define that partially saturated permeability is the only function of degree of saturation (or matrix suction) and then they do not consider hydraulic-mechanical characteristics for the analysis. However, in practice, the water flow processes in a deformable soil are influenced by soil skeleton movement and the pore water pressure changed due to seepage will lead to changes in stress and to deformation of a soil. The relationship between seepage and soil behavior causes a change of partially saturated permeability as well as saturated permeability with the lapse of time. Instability of partially saturated soil slopes due to infiltration would be analyzed from reduction of negative pore water pressure calculating the process of water flow based on predicted partially saturated permeability. Therefore, partially saturated permeability should be defined by the function of degree of saturation (or matric suction) and porosity. The paper presents the comparison between staggered and monolithic coupled analysis regarding seepage and stress deformation problems. As a result, the decrease in matric suction on soil slope from monolithic analysis is slower than that from staggered analysis.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.19 no.2
• /
• pp.265-281
• /
• 2017

For shallow tunnel constructions, settlement of the ground surface is a main issue. Recent technical developments in shield TBM tunneling technique have enabled a decrease in such settlements based on tunneling with ground deformation controls. For this objective, the tail void grouting is a common practice. Generally surface settlements in a soil of low permeability occur during a tunnel construction but also during a long period after completion of the tunnel. The long-term settlements occur mainly due to consolidation around the tunnel. The consolidation process is caused and determined by the tail void grouting which leads to an excess pore water pressure in the vicinity of the tunnel. Because of this, the grouting pressure has a strong effect on the long-term settlements in the shield tunneling. In order to investigate this effect, a series of coupled hydro-mechanical 3D finite element simulations have been performed. The results show that an increase in grouting pressure reduces the short-term settlements, but in many cases, it doesn't lead to a reduction of the final settlements after the completion of consolidation. Thereby, the existence of a critical grouting pressure is identified, at which the minimal settlements are expected.