• Tunnel and Underground Space
• /
• v.29 no.6
• /
• pp.468-479
• /
• 2019

It is essential to comprehend coupled hydro-mechanical behavior to utilize subsurface for the recent demand for underground space usage. In this study, we developed a new simulator for numerical simulation as a tool for researching to consider the various domestic field and subsurface conditions. To develop the new module, we combined OpenGeoSys, one of the scientific software package that handles fluid mechanics (H), thermodynamics (T), and rock and soil mechanics (M) in the subsurface with FLAC3D, one of the commercial software for geotechnical engineering problems reinforced. In this simulator development, we design OpenGeoSys as a master and FLAC3D as a slave via a file-based sequential coupling. We have chosen Terzaghi's consolidation problem related to single-phase fluid flow at a saturated condition as a benchmark model to verify the proposed module. The comparative results between the analytical solution and numerical analysis showed a good agreement.

• Journal of the Korean Geotechnical Society
• /
• v.31 no.9
• /
• pp.5-15
• /
• 2015

In this study, rainfall-induced slope stability and coupling effect are investigated using hydro-mechanical finite element model. This model is developed by formulating constitutive and coupled balance equations and is verified by comparing the numerical results with field matric suction. The homogeneous soil layer (soil column) and soil slope are modeled by this model, and the results of variation in matric suction, mean effective stress, porosity, displacement, factor of safety are compared with those of staggered analysis. It is found that the vertical and horizontal displacement from coupling analysis considering change in porosity is larger than that of staggered analysis. The displacement and matric suction from coupling analysis by rainfall infiltration can affect slope instability, which shows a progressive failure behavior. The lowest factor of safety is observed under short-term rainfall. This results confirm the fact that coupling analysis is needed to design soil slope under severe rain condition.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.67-67
• /
• 2016

최근 기후온난화로 인한 이상기후와 지진의 발생가능성으로 인해 직렬 혹은 병렬로 위치한 2개 이상의 댐(저수지)들의 연속 붕괴 가능성도 점점 커지고 있어 연속 댐 붕괴에 대한 비상대처계획도 수립에 대한 관심도 증대되고 있다. 국내에서는 댐(저수지) 붕괴로 인한 극한홍수해석이나 비상대처계획 수립시 국내에서는 DAMBRK, FLDWAV, HEC-RAS와 같은 1차원 수리동역학 모형이 주로 사용되어지고 있다. 하지만 1차원 모형은 흐름을 하나의 방향으로만 한정하여 해석하고, 각각의 적용 횡단면에서 동일 수위를 가지다는 가정으로 홍수범람해석을 수행하므로 정확성뿐만 아니라 실제 적용성에서도 한계를 가질 수밖에 없다. 댐(저수지) 붕괴로 인한 홍수범람해석에서 2차원 이상의 고차원 모형은 앞서 언급한 1차원 모형에 적용된 비현실적인 가정을 포함하지 않으므로 더욱 정교하고 신뢰성 있는 결과를 얻을 수 있다. 하지만, 홍수범람해석을 위한 상용 프로그램들은 단일 댐(저수지) 붕괴의 적용에는 큰 어려움이 없으나, 지형단면을 초기에 한번만 고려하는 문제로 인하여 연속 댐(저수지) 붕괴의 고려에는 한계를 가진다. 따라서 본 연구에서는 댐(저수지)의 연속 붕괴를 해석할 수 있는 2차원 수리동역학 모형을 개발하고자 한다. 각 댐(저수지)의 붕괴 전과 후의 지형 단면을 여러 번 반영할 수 있는 모형을 개발함으로써 시간차를 두고 붕괴되는 댐(저수지)의 연속 붕괴를 해석할 수 있도록 하였다. 개발모형은 2002년 태풍 루사로 인해 실제로 붕괴된 저수지에 적용되었으며, 이 저수지들은 붕괴시 시간차이를 두고 붕괴가 이루어져 개발모형의 적용 유역으로 선택하였다. 저수지의 연속붕괴 모델링을 위해 지형자료로는 저수지 단면, 댐 제체 및 여수로 붕괴 단면, 하천 단면 그리고 홍수터 지형 반영을 위한 수치지도, 경계조건으로는 저수지로의 유입유량과 하류단 조위조건이 고려되었다. 그리고 태풍 루사 당시의 기록적인 강우를 반영하기 위해 연구유역 인근에 관측된 강우를 모형에서 하천 및 홍수터에서 고려할 수 있도록 하였다.

• Tunnel and Underground Space
• /
• v.33 no.1
• /
• pp.57-69
• /
• 2023

In this study, we simulated induced seismicity in the Groningen natural gas reservoir using 2D hydro-mechanical coupled discrete element modelling (DEM). The code used is PFC2D (Particle Flow Code 2D), a commercial software developed by Itasca, and in order to apply to this study we further developed 1)initialization of inhomogeneous reservoir pressure distribution, 2)a non-linear pressure-time history boundary condition, 3)local stress field monitoring logic. We generated a 2D reservoir model with a size of 40 × 50 km² and a complex fault system, and simulated years of pressure depletion with a time range between 1960 and 2020. We simulated fault system failure induced by pressure depletion and reproduced the spatiotemporal distribution of induced seismicity and assessed its failure mechanism. Also, we estimated the ground subsidence distribution and confirmed its similarity to the field measurements in the Groningen region. Through this study, we confirm the feasibility of the presented 2D hydro-mechanical coupled DEM in simulating the deformation of a complex fault system by hydro-mechanical coupled processes.

• 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.

• Tunnel and Underground Space
• /
• v.32 no.3
• /
• pp.203-218
• /
• 2022

A hydro-mechanical study was performed to analyze the relationship between the magnitude of injection-induced seismicity and shut-in. In hydraulic analysis, the suspension of fluid injection makes the pore pressure gradient smaller while the pore pressure at the pressure front can reach the critical value for several hours after shut-in, which leads to the additional slip with wider area than during injection. The hydro-mechanical numerical analysis was performed to model the simplified fault system, and simulated the largest magnitude earthquake during shut-in stage. The effect of the abrupt suspension of fluid injection on the large magnitude earthquake was investigated in comparison with the continuous injection. In addition to the pore pressure distribution, it was found that the geometry of multiple faults and the stress redistribution are also important in evaluating the magnitude of the induced seismicity.

• Tunnel and Underground Space
• /
• v.18 no.3
• /
• pp.175-184
• /
• 2008

In this study, state-of-the-art of $CO_2$ geological sequestration as a method of greenhouse gas reduction was reviewed. Thermal-Hydraulic-Mechanically(THM) coupled simulation technology and its application to a stability analysis of geological formation due to $CO_2$ injection as well as a leakage path analysis were investigated and introduced.

• Journal of the KSME
• /
• v.50 no.4
• /
• pp.46-49
• /
• 2010

수리동역학 코드를 사용하여 폭발성형 관통자의 생성과정을 해석하고 다양한 해석기법을 개발하였다. 수치해석 결과 섬광 X선 장비를 사용한 정치시험 결과와 비교하여 해석 결과의 신뢰성을 확인하였다. 폭발성형 관통자의 관통성능 증대를 위하여 다양한 라이너 모델에 대한 수치해석을 수행하고, 그 결과를 비교하였다.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2008.10a
• /
• pp.377-389
• /
• 2008

• Tunnel and Underground Space
• /
• v.32 no.6
• /
• pp.518-529
• /
• 2022

In Task C of the DECOVALEX-2023 project, nine institutes from six nations are developing their numerical codes to simulate thermo-hydro-mechanical coupled behavior for the FE experiment performed at Mont Terri underground rock laboratory, Switzerland. Currently, Step 1 for comparing the simulation results to field data is the ongoing stage, and we used the OGS-FLAC simulator for a series of numerical simulations. As a result, temperature increase depending on the heating hysteresis was well simulated, and saturation variation in the bentonite depending on phase change was observed. However, due to the suction overestimation, relative humidity and temperature change in the bentonite and the pressure variation in the Opalinus clay showed a difference compared to the field data. From the observation, it is confirmed that the effect of the bentonite capillary pressure is dominant to the flow analysis in the disposal system. We further plan to draw improved results considering tunnel support material and accurate initial water pressure distribution. Additionally, the thermal, hydrological, and mechanical anisotropy of the Opalinus clay was well simulated. From the simulation results, we confirmed the applicability of the OGS-FLAC simulator in the disposal system analysis.