• 방사성폐기물학회지
• /
• 제21권1호
• /
• pp.23-41
• /
• 2023

Coupled thermo-hydraulic-mechanical (THM) processes are essential for the long-term performance of deep geological disposal of high-level radioactive waste. In this study, a numerical sensitivity analysis was performed to analyze the effect of rock properties on THM responses after the execution of the heater test at the Kamaishi mine in Japan. The TOUGH-FLAC simulator was applied for the numerical simulation assuming a continuum model for coupled THM analysis. The rock properties included in the sensitivity study were the Young's modulus, permeability, thermal conductivity, and thermal expansion coefficients of crystalline rock, rock salt, and clay. The responses, i.e., temperature, water content, displacement, and stress, were measured at monitoring points in the buffer and near-field rock mass during the simulations. The thermal conductivity had an overarching impact on THM responses. The influence of Young's modulus was evident in the mechanical behavior, whereas that of permeability was noticed through the change in the temperature and water content. The difference in the THM responses of the three rock type models implies the importance of the appropriate characterization of rock mass properties with regard to the performance assessment of the deep geological disposal of high-level radioactive waste.

• Structural Engineering and Mechanics
• /
• 제66권4호
• /
• pp.505-513
• /
• 2018

Steam flooding is widely used in heavy oil reservoir with coupling effects among the formation temperature change, fluid flow and solid deformation. The effective stress, porosity and permeability in this process can be affected by the multi-physical coupling of thermal, hydraulic and mechanical processes (THM), resulting in a complex interaction of geomechanical effects and multiphase flow in the porous media. Quantification of the state of deformation and stress in the reservoir is therefore essential for the correct prediction of reservoir efficiency and productivity. This paper presents a coupled fluid flow, thermal and geomechanical model employing a program (MATLAB interface code), which was developed to couple conventional reservoir (ECLIPSE) and geomechanical (ABAQUS) simulators for coupled THM processes in multiphase reservoir modeling. In each simulation cycle, time dependent reservoir pressure and temperature 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 and new porosity and permeability are obtained using volumetric strains for the next analysis step. Finally, the proposed approach is illustrated on a complex coupled problem related to steam flooding in an oil reservoir. The reservoir coupled study showed that permeability and porosity increase during the injection scenario and increasing rate around injection wells exceed those of other similar comparable cases. Also, during injection, the uplift occurred very fast just above the injection wells resulting in plastic deformation.

• 터널과지하공간
• /
• 제32권6호
• /
• pp.411-434
• /
• 2022

사용후핵연료의 심층처분 개념에서 근계영역 암반은 열-수리-역학적 복합거동을 하게 되는 것으로 잘 알려져있다. 이러한 복합거동 과정에서 암석의 여러 물성들은 변화하는데, 이러한 물성변화를 합리적으로 반영하는 경우 고준위방사성폐기물 처분장의 장기안정성의 평가를 위해 활용되는 해석 및 현장시험의 신뢰도를 향상시킬 수 있다. 이를 위해 본 기술보고에서는 암석의 열-수리-역학적 간접복합거동에 관한 국내외 연구사례를 조사하고 분석하였다. 특히, 간접복합거동의 대표적인 사례 중 지하수에 의한 포화 및 온도 증가에 따른 암석의 여러 물성 변화, 응력 변화에 의한 투수계수 변화를 중점적으로 조사·요약하였다.

• 터널과지하공간
• /
• 제30권4호
• /
• pp.359-381
• /
• 2020

벤토나이트 완충재에서의 열-수리-역학적 복합거동을 예측하기 위해 TOUGH2-MP/FLAC3D 시뮬레이터를 기반으로 개발된 Barcelona basic 모델(BBM) 해석모듈의 현장 적용성을 검토하고자 국제공동연구 DECOVALEX-2019 Task D에 참여하여 스위스 Grimsel Test Site의 현장시험(full-scale engineered barriers experiment, FEBEX) 모델링을 수행하고 현장시험에서 계측된 히터 파워, 온도, 상대습도, 응력, 포화도, 함수율 그리고 건조밀도를 계산 값과 비교하였다. 수치해석을 이용하여 시간에 따른 히터 파워와 온도 변화는 전반적으로 잘 재현되었지만, 히터 1과 히터 2에서의 파워 차이를 계산할 수는 없었으며 이를 개선하기 위해서는 FEBEX 터널 주변에 분포하는 황반암과 시험장치 및 벤토나이트 블록의 설치 공정을 반영할 필요가 있을 것으로 판단된다. 상대습도 변화와 분포 역시 전반적으로 잘 모사되었으나, 수치해석에서 히터 부근에서의 재포화과정이 상대적으로 빠르게 진행된 것으로 보아 수리모델에 대한 일부 수정이 필요할 것으로 보인다. 현장시험에서는 벤토나이트 완충재와 암반 사이에 틈이 존재하지만 수치해석에서는 완벽하게 접촉하고 있는 것으로 가정하였기 때문에 운영 초기의 응력 변화는 다소 차이를 보였지만, 전반적으로 유사한 경향을 보이는 것으로 나타났다. 해체 이후 측정한 포화도, 함수율, 그리고 건조밀도의 분포 역시 전반적으로 잘 재현되었지만, 건조밀도가 터널 중심과 히터부근에서 조금 크게 계산되어 벤토나이트 블록의 투수계수가 상대적으로 작은 값으로 반영되어 포화도와 함수율이 작게 계산된 것으로 보이며, 이를 개선하기 위해서는 건조밀도에 따른 투수계수 모델에 일부 수정이 필요할 것으로 판단된다. 본 연구의 결과를 토대로 수치모델을 수정하고 추가적인 연구를 수행한다면, 보다 나은 해석 결과와 벤토나이트 완충재에서의 THM 복합거동을 좀 더 현실적으로 예측할 수 있을 것으로 판단된다.

• 방사성폐기물학회지
• /
• 제4권2호
• /
• pp.117-131
• /
• 2006

한국의 고준위폐기물 기준 처분 시스템의 공학적 방벽에서의 T-H-M(Thermo-Hydro-Mechanical) 거동 실증을 위한 KENTEX(KAERI Engineering-scale T-H-M Experiment for Engineered Barrier System)실험 장치를 대상으로 열-수리-역학 연동현상 해석을 하여 온도, 포화도 및 응력의 변화를 예측하였다. 그리고 이들 변수와 열-수리-역학의 연동현상에 사용된 세물성법칙인 탄성물성법칙, 공극탄성 물성법칙 및 공극탄성-소성 물성법칙과의 관계를 분석하였다. 열-수리-역학 연동현상을 계산하는 데는 상용 유한요소 코드인 ABAQUS를 사용하였다. 열 계산에서 벤토나이트 내 온도는 히터 가열 후 초기에는 급격히 증가하다가 얼마의 시간이 경과한 후에는 거의 일정한 값에 도달하였다. 이 도달시간은 약 37.5일로 반경방향의 모든 지점(H=0.68m 일때)에서 정상상태에 도달한 것을 알 수 있었다. 즉, 히터와 벤토나이트 경계면에서는 $90^{\circ}C$, 벤토나이트와 외부 셀 경계면에서는 약 $70^{\circ}C$를 유지하였다. 열-수리-역학 연동현상 계산에서 시간에 따른 벤토나이트 포화도는 탄성 물성법칙, 공극탄성 물성법칙 및 공극탄성-소성 물성법칙의 세 경우 모두 거의 차이가 없었다. 열-수리-역학 계산 결과와 수리-역학 계산 결과의 비교에서 온도의 증가는 탄성 물성법칙 및 공극탄성 물성법칙 각각에 대해 시간이 경과함에 따라 포화도가 증가함을 초래해 포화가 빨리 진행됨을 알 수 있었다. 특히 히터에 가까운 쪽에서는물이 침투하고 있는 쪽 보다 포화도 증가가 큰 것으로 나타나 벤토나이트가 물로 포화되기 전의초기상태가 온도의 영향을 많이 받는 것을 알 수 있었다. 또한 응력은 세 물성 법칙 모두 시간의 경과에 따라 증가하는 경향을 보이나 탄성 물성법칙의 경우가 다른 두 경우보다 현저한 변화를 보이는데 이는 변형율이 탄성한계를 넘어서도 계속 작용하여 공극비 변화를 고려한 다른 두 물성법칙과 차이가 있음을 나타내고 있다. 그러나 공극탄성 물성법칙 및 공극탄성-소성 물성법칙의 경우에 열-수리-역학 계산 결과와 수리-역학 계산 결과를 비교하면 시간이 경과함에 따라 응력은 증가하지만 온도의 변화에 따른 서로의 응력의 차이는 작은 것을 알 수 있다. 즉 온도변화의 영향보다는 시간에 따른 포화도 변화의 영향이 더 큰 것으로 생각된다. 따라서 벤토나이트의 열-수리-역학 연동현상 해석에서 벤토나이트는 온도의 증가로 포화가 빨라지고, 포화도 증가는 응력을 증가시키는 결과를 보이므로 공극비, 열팽창 및 팽윤압 등의 영향을 받고 있는 것으로 이해된다. 그래서 벤토나이트의 열-수리-역학 연동현상 해석에서 벤토나이트는 공극비, 열팽창 및 팽윤압 등의 영향을 받으므로 탄성과 소성을 동시에 고려할 수 있는 물성법칙을 선택하는 것이 바람직하다.