• Geomechanics and Engineering
• /
• v.10 no.4
• /
• pp.423-436
• /
• 2016

Currently, there is a great interest in the coupling between multiphase fluid flow and geomechanical effects in hydrocarbon reservoirs and surrounding rocks. The ideal solution for this coupled problem is to introduce the geomechanical effects through the stress analysis solution and implement an algorithm, which assures that the equations governing the flow and stress analyses are obeyed in each time step. This paper deals with the implementation of a program (FORTRAN90 interface code), which was developed to couple conventional reservoir (ECLIPSE) and geomechanical (ABAQUS) simulators, using a partial coupling algorithm. The explicit coupled hydro-mechanical behavior of Iranian field during depletion and $CO_2$ injection is studied using the soils consolidation procedure available in ABAQUS. Time dependent reservoir pressure fields obtained from three dimensional compositional reservoir models were transferred into finite element reservoir geomechanical models in ABAQUS as multi-phase flow in deforming reservoirs cannot be performed within ABAQUS. The FEM analysis of the reservoir showed no sign of plastic strain under production and $CO_2$ injection scenarios in any part of the reservoir and the stress paths do not show a critical behavior.

• Tunnel and Underground Space
• /
• v.33 no.4
• /
• pp.265-280
• /
• 2023

Water leakage in tunnels is a defect that can affect tunnel stability and the ground movement by changing the stress and pore water pressure of the surrounding ground. Long-term or large-scale water leaks may lead to damage of tunnel structure and the surrounding environment, such as tunnel lining instability and ground surface settlement. The present study numerically investigated the effects of water leakage on the structural stability of a tunnel and the ground behavior. The tunnel was assumed to be under undrained conditions for preventing the inflow of the surrounding water and leaks occurred in the concrete lining after completion of the tunnel construction. A coupled hydro-mechanical analysis using a TOUGH-FLAC simulator developed in Python was conducted for assessing the leakage induced-behavior of the tunnel structure and ground under different conditions of the amount and location of water leak. Additionally, the effect of hydro-mechanical coupling terms on the results of coupled response was investigated and discussed.

• Tunnel and Underground Space
• /
• v.28 no.5
• /
• pp.400-425
• /
• 2018

This study presents the research results and current status of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to simulate the coupled hydro-mechanical behavior of fault, including slip or reactivation, induced by water injection. The first research step of Task B is a benchmark simulation which is designed for the modelling teams to familiarize themselves with the problem and to set up their own codes to reproduce the hydro-mechanical coupling between the fault hydraulic transmissivity and the mechanically-induced displacement. We reproduced the coupled hydro-mechanical process of fault slip 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. A methodology to formulate the hydro-mechanical coupling relations of two different hydraulic aperture models and link the solid element of TOUGH2 and the interface element of FLAC3D was suggested. In addition, we developed a coupling module to update the changes in geometric features (mesh) and hydrological properties of fault caused by water injection at every calculation step for TOUGH-FLAC simulator. Then, the transient responses of the fault, including elastic deformation, reactivation, progressive evolutions of pathway, pressure distribution and water injection rate, to stepwise pressurization were examined during the simulations. The results of the simulations suggest that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOLVAEX-2019 Task B and validated using the field data from fault activation experiments in a further study.

• Tunnel and Underground Space
• /
• v.34 no.2
• /
• pp.127-142
• /
• 2024

The robustness of a numerical method means that its computational performance is maintained under various modeling conditions. New numerical methods or codes need to be assessed for robustness through benchmark testing. The TOUGH-FLAC modeling approach has been applied to various fields such as subsurface carbon dioxide storage, geological disposal of spent nuclear fuel, and geothermal development both domestically and internationally, and the modeling validity has been examined by comparing the results with experimental measurements and other numerical codes. In the present study, a benchmark test of the TOUGH-FLAC approach was performed based on a coupled thermal-hydro-mechanical behavior problem with an analytical solution. The analytical solution is related to the temperature, pore water pressure, and mechanical behavior of a fully saturated porous medium that is subjected to a point heat source. The robustness of the TOUGH-FLAC approach was evaluated by comparing the analytical solution with the results of numerical simulation. Additionally, the effects of thermal-hydro-mechanical coupling terms, fluid phase change, and timestep on the computation of coupled behavior were investigated.

• Tunnel and Underground Space
• /
• v.23 no.6
• /
• pp.470-479
• /
• 2013

Hydraulic fracturing is used as a method for promoting the fluid flow in the rock and, in the energy field such as geothermal development and the development of sales gas, many studies has been actively conducted. In many cases, hydraulic fracturing is not performed in isotropic rock and especially in the case of sedimentary rocks, hydraulic fracturing is conducted in the transverse isotropic rock. The direction of the crack growth on hydraulic fracturing does not necessarily coincides with the direction of maximum principal stress in the transverse isotropic rock. Therefore, in this study, bonded particle model with hydro-mechanical coupling analysis was adopted for analyzing the characteristics of hydraulic fracturing in transverse isotropic rock. In addition, experiments of hydraulic fracturing were conducted in laboratory-scale to verify the validity of numerical analysis. In this study, the crack growth and crack patterns showed significant differences depending on the viscosity of injection fluid, the angle of bedding plane and the influence of anisotropy. In the case of transverse isotropic model, the shear crack growth due to hydraulic fracturing appeared prominently.

• Tunnel and Underground Space
• /
• v.9 no.1
• /
• pp.1-11
• /
• 1999

A finite element code was developed to analyze coupled thermo-hydro-mechanical phenomena. This code is based on the finite element formulation provided by Noorishad et al. (1984) and Joint behavior was simulated Goodman's joint constitutive model. The developed code was applied for T-H-M coupling analysis for two kinds of shaft models, with a joint or without a joint respectively. For a model without a joint, temperature increased from the shaft wall to outward evidently. The radial displacement showed opposite directions of outward and inward at some distance from shaft wall. For a model with a joint, closure of joint was found due to thermal expansion. The temperature distribution along a joint showed relatively lower than that of rock matrix because of low thermal conductivity and high specific heat of water. And it could be concluded that effects of thermal flow to joint were more than that of hydraulic flow in a rock mass.

• Tunnel and Underground Space
• /
• v.27 no.3
• /
• pp.146-160
• /
• 2017

$CO_2$ sequestration for alleviating global warming is a hot issue in the world. In this study, TOUGH2 and FLAC3D were combined for analyzing the hyro-mechanical coupling behaviors expected in $CO_2$ sequestration and applied it to Sleipner project carried out in Norway. In the analysis, the influence of pore pressure on in situ stress was considered and the influence of caprock permeability on hydro-mechanical behaviors was analyzed. In the condition of constant injection rate, pressure and saturation at the injection well, liquid and gas saturation in rock, major and minor stress variations with time and distance from the injection well, and horizontal and vertical displacements after injection could be investigated. The major principal stress was quickly increased in the early stage and then slowly decreased to a stable value, which was higher than the initial value. In contrast, the minor principal stress returned to initial value after some increase in the early stage. Surface upheaval was steadily increased and it was up to 15mm in 2 years after injection. When the caprock's permeability was changed from $3e-15m^2{\sim}3e-18m^2$, it was found that the injection well pressure and surface upheaval were inversely propotional to the permeability.

• Nuclear Engineering and Technology
• /
• v.53 no.12
• /
• pp.3990-4002
• /
• 2021

CSPACE (Core meltdown, Safety and Performance Analysis CodE for nuclear power plants) for a simulation of severe accident progression in a Pressurized Water Reactor (PWR) is developed by coupling of verified system thermal hydraulic code of SPACE (Safety and Performance Analysis CodE for nuclear power plants) and core damage progression code of COMPASS (Core Meltdown Progression Accident Simulation Software). SPACE is responsible for the description of fluid state in nuclear system nodes, while COMPASS is responsible for the prediction of thermal and mechanical responses of core fuels and reactor vessel heat structures. New heat transfer models to each phase of the fluid, flow blockage, corium behavior in the lower head are added to COMPASS. Then, an interface module for the data transfer between two codes was developed to enable coupling. An implicit coupling scheme of wall heat transfer was applied to prevent fluid temperature oscillation. To validate the performance of newly developed code CSPACE, we analyzed typical severe accident scenarios for OPR1000 (Optimized Power Reactor 1000), which were initiated from large break loss of coolant accident, small break loss of coolant accident, and station black out accident. The results including thermal hydraulic behavior of RCS, core damage progression, hydrogen generation, corium behavior in the lower head, reactor vessel failure were reasonable and consistent. We demonstrate that CSPACE provides a good platform for the prediction of severe accident progression by detailed review of analysis results and a qualitative comparison with the results of previous MELCOR analysis.

• Tunnel and Underground Space
• /
• v.8 no.1
• /
• pp.17-25
• /
• 1998

This paper investigates coupled thermal, mechanical and hydraulic phenomena in deep rock mass especially for underground heat storage system. Firstly, concepts of underground heat storage were presented and coupling phenomena in this area were illustrated. In order to understand the basic mechanism of thermal, hydraulic and deformation behavior in rock cavern disturbed by thermal gradient about 10$0^{\circ}C$, various numerical experiments were conducted using several codes. The study involves the behavior of fractured rock mass including rock joint. In spite of the limitation of codes modelling fully coupled effects, these codes could be applied in analysis of underground heat storage. The heat loss in rock mass, which is a major factor in heat storage, is insignificant in all results.

• Journal of the Korean Geotechnical Society
• /
• v.27 no.3
• /
• pp.95-102
• /
• 2011

A model of solid-water-air coupling in triphasic mixture is compared to that of solid-water coupling in biphasic mixture with an application to partially saturated soils. Based on thermodynamics, the mathematical framework governing the behavior of a partially saturated soil is derived by using balance equations, and numerical implementation through drainage experiment of a sand column is carried out to validate the obtained formulations. The role of the air phase in the hydro-mechanical behavior of triphasic mixture can be analyzed from the interaction among phases and from the solid skeleton's constitutive behavior, and the three-phase model found applications in geotechnical engineering problems, such as $CO_2$ sequestration and air storage in an aquifer.