수치해석 기법의 강건성은 다양한 모델링 조건에서 계산 성능이 유지되는 것을 의미하며, 새로운 해석기법 또는 수치코드는 벤치마크 테스트를 통해 강건성이 평가될 필요가 있다. TOUGH-FLAC 모델링 기법은 국내외적으로 이산화탄소 지중저장, 사용후핵연료 지층처분, 지열 개발 등 다양한 분야에 적용되었으며, 실험 계측자료, 다른 수치코드들과의 결과 비교를 통해 모델링 유효성이 분석되었다. 본 연구에서는 해석해를 갖는 열-수리-역학적 복합거동 문제를 토대로 TOUGH-FLAC 기법의 벤치마크 테스트를 수행하였다. 적용된 해석해는 완전히 포화된 지반에 점열원 작용 시 주변매질의 온도, 간극수압, 역학적 거동과 관계되며, 해석해와 수치모사 결과를 비교하여 TOUGH-FLAC 기법의 강건성이 평가되었다. 또한, 열-수리-역학 해석의 연계항, 유체 상변화, 시간증분이 복합거동 계산에 미치는 영향을 조사하였다.
고준위방사성폐기물 처분장의 장기 안전성 확보를 위해서는 공학적방벽 및 천연방벽 내에서 발생하는 복잡한 열-수리-역학-화학적(THMC) 복합거동 해석에 대한 이해가 필수적이다. 특히 고준위방사성폐기물에서 발생하는 열로 인해 암반 및 완충재 내의 지하수에서 압력 증가 및 상변화가 발생하게 되며, 지하수의 유입으로 인해 공학적방벽 내 포화도가 변화하게 된다. 또한 포화도의 변화는 완충재 내에서의 열전달 및 다상 유동 특성에 영향을 미치게 된다. 따라서 복합거동 특성의 복잡성으로 인해 수치해석은 처분시스템에서의 THMC 복합거동 평가와 예측 및 안전성 평가에 있어 강점을 지니고 있으며, DECOVALEX 국제공동연구는 THMC 복합거동에 대한 이해도 증진 및 해석기법 검증을 목적으로 1992년부터 시작되었다. 국내에서는 2008년부터 한국원자력연구원이 지속적으로 참여하여 연구를 수행하고 있으며, 본 기술보고에서는 현재 진행 중인 DECOVALEX-2023의 주요 연구내용을 국내 암반 및 지반공학자들에게 소개하였다.
DECOVALEX 프로젝트는 고준위방사성폐기물 처분 시스템에서 발생하는 복잡한 열-수리-역학-화학적(THMC) 복합거동에 대해 수치해석을 통해 보다 더 깊이 이해하기 위해 수행되고 있는 대표적 국제공동 연구이다. 현재 DECOVALEX-2023이 7개 task를 기반으로 진행 중이며, 이 중 Task C는 Mont-Terri 지하연구시설에서 수행된 실규모 정치(FE) 시험을 대상으로 처분시스템 내 THM 복합거동을 모사하는 것을 목표로 하고 있다. 본 연구에서는 자체개발 해석 코드인 OGS-FLAC을 활용하여 수치해석 연구를 수행하였다. 수치모델에서는 FE 시험과 같이 일정 출력의 히터를 수평으로 위치시켰으며, 주어진 계측지점에서 압력 분포, 온도 변화, 역학적 변형을 계측하였다. 완충재 내부로 유입되는 유체 흐름은 완충재의 흡입력으로 인해 발생하였으며, 주변 영역에서는 열 팽창 및 열 압력이 지배적으로 작용함을 확인하였다. 해석 결과는 향후 타 참여 그룹 및 실험 결과와 비교 검증을 수행할 계획이다.
고준위방사성폐기물의 처분터널 및 처분공 간격을 결정하고 처분시스템의 성능을 평가하기 위해서는 열-수리-역학적인 복합 거동 변화에 대한 이해가 반드시 필요하고 이를 반영하여 해석해야만 한다. 하지만 한국형 기준 처분시스템에서의 처분터널 및 처분공 간격을 결정하기 위해 수행된 기존의 연구들은 이러한 복합거동 특성을 반영하지 않고 열 해석 결과만을 근거로 처분시스템을 설계하였다. 따라서 본 연구에서는 열-수리-역학적인 복합거동 특성을 반영하여 한국형 기준 처분시스템의 성능을 TOUGH2-MP/FLAC3D를 이용하여 평가하였다. 고준위방사성폐기물이 처분된 이후 방사성 붕괴열에 의해 처분시스템의 온도는 급격히 증가하다가 붕괴열의 감소로 온도는 서서히 감소하였으며, 해석 기간 1,000년 동안 벤토나이트 완충재의 최고 온도는 설계 기준인 $100^{\circ}C$ 이하로 유지되는 것으로 나타났다. 처분용기와 벤토나이트 완충재의 계면에서의 최고 온도는 약 3.21년이 지난 시점에 용기의 중간 지점에서 약 $96.2^{\circ}C$로 나타났으며, 암반에서의 최고 온도는 폐쇄 후 약 17년이 지난 시점에서 약 $68.2^{\circ}C$로 계산되었다. 처분용기 부근 벤토나이트 완충재는 처분 초기에 온도 변화에 따른 건조현상이 발생하여 포화도가 감소하지만, 시간이 지남에 따라 주변 암반으로부터의 지하수 유입에 의해 포화도가 증가하는 것으로 계산되었다. 이후, 벤토나이트 완충재 및 뒷채움재 모두 약 266년 이후 완전히 포화되는 것으로 계산되었다. 처분시스템에서의 온도 변화에 따른 열응력 그리고 벤토나이트 완충재 및 뒷채움재의 팽윤압으로 인한 응력 변화가 처분장 주변 암반에 미치는 영향을 평가하고자 수치해석에서 계산된 응력을 스폴링 강도(spalling strength)와 Mohr-coulomb 파괴 기준식과 비교하였다. 계산 결과 일축압축강도와 스폴링 강도에 도달하지 않는 것으로 나타나 처분시스템이 스폴링에 의한 파괴는 나타나지 않을 것으로 판단되며, Mohr-coulomb 파괴 기준 역시 충족하는 것으로 나타났다. 본 연구에서 사용된 수치해석 코드와 방법론은 다양한 조건에서의 한국형 기준 처분시스템에 대한 성능평가뿐만 아니라, 복층 처분시스템에 대한 설계와 성능평가에 활용될 수 있을 것으로 판단된다.
본 연구에서는 천부의 암반 공동에 대용량 고온의 열에너지를 저장하는 경우 주변 암반에 야기되는 열-수리-역학적 연계거동을 살펴보고, 이에 지하수위와 암반 투수계수 등 수리적 조건이 미치는 영향을 검토하였다. 해석대상을 투수계수가 비교적 낮은 수준($10^{-17}m^2$)인 결정질 암반으로 가정할 때 열에너지 저장으로 인한 암반 거동에 지하수가 미치는 영향은 크지 않을 것으로 예측되었다. 저장 공동이 지하수위 하부에 위치하는 경우의 온도, 주응력, 변위 분포 등은 저장공동이 불포화대에 위치하는 경우와 거의 동일하게 나타났다. 암반내 열전달 특성은 암반의 투수계수에 매우 큰 영향을 받는 것으로 나타났다. 암반의 투수계수를 $10^{-15}m^2$ 이하로 가정한 경우 열전달은 주로 암반에 전도에 의한 것으로 판단할 수 있었으나, 투수계수를 $10^{-12}m^2$으로 가정하는 경우 지하수 대류에 의한 상향 열유동이 뚜렷이 관찰되었다. 암반 투수계수의 크기에 따라 열수의 대류나 비등으로 인한 상변화 등 복합적인 유동 특성을 나타났으며, 온도, 압력, 포화도 분포가 상이하게 발달하였다.
DECOVALEX-2023 Task C에서는 6개국 9개 참여 기관들이 스위스 Mont Terri 지하처분연구시설에 서 수행된 FE 실험을 대상으로 열-수리-역학 복합거동 모사를 위한 해석코드 개발을 수행하고 있다. 현재 현장시험 결과와 비교 분석을 위한 Step 1이 진행되고 있으며, 본 연구진은 OGS-FLAC 해석 시뮬레이터를 활용하여 일련의 해석을 진행하였다. 해석 결과 히터 가열에 따른 온도 상승이 잘 구현되었고, 상 변화에 따른 완충재 내 포화도 변화를 관측할 수 있었다. 반면 완충재 흡입력의 과대평가로 완충재 내 상대습도, 온도 변화 및 Opalinus 점토암 내 압력 변화가 현장 결과와 다소간 차이를 나타내는 것을 확인하였다. 이를 통해 완충재 흡입력이 처분시스템 해석 시 유동 해석 결과에 지배적인 영향을 미침을 확인할수 있었으며, 향후 지보재 및 초기 수압 모사 개선을 통해 향상된 결과를 도출하고자 한다. 또한, Opalinus 점토암의 열, 수리, 역학적 이방성이 잘 구현되었으며 해석 결과를 통해 OGS-FLAC 시뮬레이터의 처분시스템 해석 적용성을 확인하였다.
Barcelona Basic Model(BBM)은 응력의 변화에 따른 부피변화뿐만 아니라 흡입력의 변화에 따른 팽윤거동을 설명할 수 있으며, 흡입력 변화에 따른 점착력과 선행압밀응력의 변화와 온도변화에 따른 선행압밀응력의 변화를 고려할 수 있다. 따라서, 고준위방사성폐기물 처분시스템에서 공학적방벽재로 고려되고 있는 벤토나이트 완충재의 열-수리-역학적 복합거동을 예측 및 분석하는 것에 많이 활용되고 있다. 그러나 우리나라의 암반 및 지반 공학자들에게 잘 알려져 있지 않기 때문에 BBM을 소개하고자 한다. BBM은 불포화 토질의 역학적 거동을 모사하기 위해 Modified Cam Clay(MCC) 모델을 확장하여 만들어 졌기 때문에 본 고에서는 먼저 MCC 모델을 간략하게 소개하고, 열-탄소성 모델인 BBM을 상세히 소개하였다.
암반내 불연속면이나 지반-구조물의 접촉면은 열-수리-역학적으로 연계된 거동 특성을 보이므로, 온전한 지배방정식에 근거한 경계면 요소의 개발이 필요하다. 본 논문은 경계면 요소에 대한 힘평형 방정식, 유체의 연속방정식 그리고 에너지 평형 방정식을 유도하였다. 그리고 경계면 요소에 대한 탄소성 역학 모델의 강성행렬을 제시하였다. 개발된 유한요소는 2차원 조건에서 변위는 6절점, 수압과 온도는 4절점을 사용한다. 단층내 유체 주입에 대한 완전연계된 THM 해석은 단층내의 유효응력 감소와 주위 암반의 온도 수축에 의한 주입압의 복합적인 변화을 모델링 할수 있었다. 하지만, 열적 현상을 무시한 HM해석은 수리-역학적 변수를 과다하게 산정하였다.
본 연구에서는 입자기반 개별요소모델(grain-based distinct element model, GBDEM)을 이용하여 결정질 암석 내 포함된 균열의 열-역학적 거동을 평가할 수 있는 수치해석기법을 제시하고 열에 의한 균열의 미끄러짐 거동을 해석하였다. 이는 DECOVALEX-2023 프로젝트 Task G의 일환으로 수행된 벤치마크 모델링 연구로, Task G는 결정질 암반 내 균열의 열-수리-역학적 복합거동을 해석하기 위한 수치해석기법을 개발하는 데에 목표가 있다. 여기에서는 Voronoi diagram을 이용하여 다면체 개별입자의 집합체로서 해석모델을 생성하고, 입자 및 입자간 접촉에서 발생하는 열-역학적 거동을 개별요소프로그램인 3DEC을 통해 해석하였다. 암석 시험편의 탄성거동을 재현하기 위하여 등가연속체 개념을 적용하여 입자와 접촉의 미시물성을 산정하였으며, 균열에 상응하는 접촉에는 Coulomb slip model을 부여하여 인장강도와 전단강도를 갖는 불연속면을 모사하였다. 경계응력과 열응력에 의한 균열의 거동을 수치적으로 모델링하였으며, 경계조건에 따라 균열의 미끄러짐이 발생하는 열-역학적 메커니즘을 정량적으로 분석하였다. 해석 결과, 본 연구에서 제시한 해석모델이 암석 내 열팽창과 열응력의 증가, 균열 응력과 변위, 경계조건의 영향 등을 합리적으로 재현하고 있음을 확인하였다. 본 연구의 해석모델은 Task G에 참여하는 국외 연구팀들과의 의견 교류와 워크숍을 통해 지속적으로 개선하는 한편, 향후 실내실험에 적용하여 타당성을 검증할 예정이다.
최근 대두되고 있는 지하공간 활용의 필요성 증대에 따라 지하공간의 안정적인 활용을 위해서는 수리-역학적 복합거동에 대한 이해가 필수적이다. 본 연구에서는 향후 국내 다양한 현장 및 지중 조건을 고려한 연구를 수행하는데 도구가 될 수 있는 수치해석을 위한 시뮬레이터를 개발하였다. 이를 위해 지중의 유체역학과 열역학, 역학을 해석하는 과학 소프트웨어 중 하나인 OpenGeoSys와 지반의 여러 역학적 문제들을 해석하는데 전문성을 지닌 FLAC3D를 연계하기로 하였다. 연동해석 기법에 대한 조사를 통해 OpenGeoSys를 주(master)로, FLAC3D를 종속(slave)으로 하여 연계 해석하는 파일 기반의 순차적 연계 방식으로 알고리즘을 개발하였다. 또한, 개발된 시뮬레이터의 검증을 위해 포화된 상태에서의 단상 유체 거동과 관련한 Terzaghi의 압밀 문제를 벤치마크 모델로 선정하여 비교 검증을 수행하였으며, 해석적 해와 수치해석의 결과가 잘 일치하는 것을 확인하였다.
본 웹사이트에 게시된 이메일 주소가 전자우편 수집 프로그램이나
그 밖의 기술적 장치를 이용하여 무단으로 수집되는 것을 거부하며,
이를 위반시 정보통신망법에 의해 형사 처벌됨을 유념하시기 바랍니다.
[게시일 2004년 10월 1일]
이용약관
제 1 장 총칙
제 1 조 (목적)
이 이용약관은 KoreaScience 홈페이지(이하 “당 사이트”)에서 제공하는 인터넷 서비스(이하 '서비스')의 가입조건 및 이용에 관한 제반 사항과 기타 필요한 사항을 구체적으로 규정함을 목적으로 합니다.
제 2 조 (용어의 정의)
① "이용자"라 함은 당 사이트에 접속하여 이 약관에 따라 당 사이트가 제공하는 서비스를 받는 회원 및 비회원을
말합니다.
② "회원"이라 함은 서비스를 이용하기 위하여 당 사이트에 개인정보를 제공하여 아이디(ID)와 비밀번호를 부여
받은 자를 말합니다.
③ "회원 아이디(ID)"라 함은 회원의 식별 및 서비스 이용을 위하여 자신이 선정한 문자 및 숫자의 조합을
말합니다.
④ "비밀번호(패스워드)"라 함은 회원이 자신의 비밀보호를 위하여 선정한 문자 및 숫자의 조합을 말합니다.
제 3 조 (이용약관의 효력 및 변경)
① 이 약관은 당 사이트에 게시하거나 기타의 방법으로 회원에게 공지함으로써 효력이 발생합니다.
② 당 사이트는 이 약관을 개정할 경우에 적용일자 및 개정사유를 명시하여 현행 약관과 함께 당 사이트의
초기화면에 그 적용일자 7일 이전부터 적용일자 전일까지 공지합니다. 다만, 회원에게 불리하게 약관내용을
변경하는 경우에는 최소한 30일 이상의 사전 유예기간을 두고 공지합니다. 이 경우 당 사이트는 개정 전
내용과 개정 후 내용을 명확하게 비교하여 이용자가 알기 쉽도록 표시합니다.
제 4 조(약관 외 준칙)
① 이 약관은 당 사이트가 제공하는 서비스에 관한 이용안내와 함께 적용됩니다.
② 이 약관에 명시되지 아니한 사항은 관계법령의 규정이 적용됩니다.
제 2 장 이용계약의 체결
제 5 조 (이용계약의 성립 등)
① 이용계약은 이용고객이 당 사이트가 정한 약관에 「동의합니다」를 선택하고, 당 사이트가 정한
온라인신청양식을 작성하여 서비스 이용을 신청한 후, 당 사이트가 이를 승낙함으로써 성립합니다.
② 제1항의 승낙은 당 사이트가 제공하는 과학기술정보검색, 맞춤정보, 서지정보 등 다른 서비스의 이용승낙을
포함합니다.
제 6 조 (회원가입)
서비스를 이용하고자 하는 고객은 당 사이트에서 정한 회원가입양식에 개인정보를 기재하여 가입을 하여야 합니다.
제 7 조 (개인정보의 보호 및 사용)
당 사이트는 관계법령이 정하는 바에 따라 회원 등록정보를 포함한 회원의 개인정보를 보호하기 위해 노력합니다. 회원 개인정보의 보호 및 사용에 대해서는 관련법령 및 당 사이트의 개인정보 보호정책이 적용됩니다.
제 8 조 (이용 신청의 승낙과 제한)
① 당 사이트는 제6조의 규정에 의한 이용신청고객에 대하여 서비스 이용을 승낙합니다.
② 당 사이트는 아래사항에 해당하는 경우에 대해서 승낙하지 아니 합니다.
- 이용계약 신청서의 내용을 허위로 기재한 경우
- 기타 규정한 제반사항을 위반하며 신청하는 경우
제 9 조 (회원 ID 부여 및 변경 등)
① 당 사이트는 이용고객에 대하여 약관에 정하는 바에 따라 자신이 선정한 회원 ID를 부여합니다.
② 회원 ID는 원칙적으로 변경이 불가하며 부득이한 사유로 인하여 변경 하고자 하는 경우에는 해당 ID를
해지하고 재가입해야 합니다.
③ 기타 회원 개인정보 관리 및 변경 등에 관한 사항은 서비스별 안내에 정하는 바에 의합니다.
제 3 장 계약 당사자의 의무
제 10 조 (KISTI의 의무)
① 당 사이트는 이용고객이 희망한 서비스 제공 개시일에 특별한 사정이 없는 한 서비스를 이용할 수 있도록
하여야 합니다.
② 당 사이트는 개인정보 보호를 위해 보안시스템을 구축하며 개인정보 보호정책을 공시하고 준수합니다.
③ 당 사이트는 회원으로부터 제기되는 의견이나 불만이 정당하다고 객관적으로 인정될 경우에는 적절한 절차를
거쳐 즉시 처리하여야 합니다. 다만, 즉시 처리가 곤란한 경우는 회원에게 그 사유와 처리일정을 통보하여야
합니다.
제 11 조 (회원의 의무)
① 이용자는 회원가입 신청 또는 회원정보 변경 시 실명으로 모든 사항을 사실에 근거하여 작성하여야 하며,
허위 또는 타인의 정보를 등록할 경우 일체의 권리를 주장할 수 없습니다.
② 당 사이트가 관계법령 및 개인정보 보호정책에 의거하여 그 책임을 지는 경우를 제외하고 회원에게 부여된
ID의 비밀번호 관리소홀, 부정사용에 의하여 발생하는 모든 결과에 대한 책임은 회원에게 있습니다.
③ 회원은 당 사이트 및 제 3자의 지적 재산권을 침해해서는 안 됩니다.
제 4 장 서비스의 이용
제 12 조 (서비스 이용 시간)
① 서비스 이용은 당 사이트의 업무상 또는 기술상 특별한 지장이 없는 한 연중무휴, 1일 24시간 운영을
원칙으로 합니다. 단, 당 사이트는 시스템 정기점검, 증설 및 교체를 위해 당 사이트가 정한 날이나 시간에
서비스를 일시 중단할 수 있으며, 예정되어 있는 작업으로 인한 서비스 일시중단은 당 사이트 홈페이지를
통해 사전에 공지합니다.
② 당 사이트는 서비스를 특정범위로 분할하여 각 범위별로 이용가능시간을 별도로 지정할 수 있습니다. 다만
이 경우 그 내용을 공지합니다.
제 13 조 (홈페이지 저작권)
① NDSL에서 제공하는 모든 저작물의 저작권은 원저작자에게 있으며, KISTI는 복제/배포/전송권을 확보하고
있습니다.
② NDSL에서 제공하는 콘텐츠를 상업적 및 기타 영리목적으로 복제/배포/전송할 경우 사전에 KISTI의 허락을
받아야 합니다.
③ NDSL에서 제공하는 콘텐츠를 보도, 비평, 교육, 연구 등을 위하여 정당한 범위 안에서 공정한 관행에
합치되게 인용할 수 있습니다.
④ NDSL에서 제공하는 콘텐츠를 무단 복제, 전송, 배포 기타 저작권법에 위반되는 방법으로 이용할 경우
저작권법 제136조에 따라 5년 이하의 징역 또는 5천만 원 이하의 벌금에 처해질 수 있습니다.
제 14 조 (유료서비스)
① 당 사이트 및 협력기관이 정한 유료서비스(원문복사 등)는 별도로 정해진 바에 따르며, 변경사항은 시행 전에
당 사이트 홈페이지를 통하여 회원에게 공지합니다.
② 유료서비스를 이용하려는 회원은 정해진 요금체계에 따라 요금을 납부해야 합니다.
제 5 장 계약 해지 및 이용 제한
제 15 조 (계약 해지)
회원이 이용계약을 해지하고자 하는 때에는 [가입해지] 메뉴를 이용해 직접 해지해야 합니다.
제 16 조 (서비스 이용제한)
① 당 사이트는 회원이 서비스 이용내용에 있어서 본 약관 제 11조 내용을 위반하거나, 다음 각 호에 해당하는
경우 서비스 이용을 제한할 수 있습니다.
- 2년 이상 서비스를 이용한 적이 없는 경우
- 기타 정상적인 서비스 운영에 방해가 될 경우
② 상기 이용제한 규정에 따라 서비스를 이용하는 회원에게 서비스 이용에 대하여 별도 공지 없이 서비스 이용의
일시정지, 이용계약 해지 할 수 있습니다.
제 17 조 (전자우편주소 수집 금지)
회원은 전자우편주소 추출기 등을 이용하여 전자우편주소를 수집 또는 제3자에게 제공할 수 없습니다.
제 6 장 손해배상 및 기타사항
제 18 조 (손해배상)
당 사이트는 무료로 제공되는 서비스와 관련하여 회원에게 어떠한 손해가 발생하더라도 당 사이트가 고의 또는 과실로 인한 손해발생을 제외하고는 이에 대하여 책임을 부담하지 아니합니다.
제 19 조 (관할 법원)
서비스 이용으로 발생한 분쟁에 대해 소송이 제기되는 경우 민사 소송법상의 관할 법원에 제기합니다.
[부 칙]
1. (시행일) 이 약관은 2016년 9월 5일부터 적용되며, 종전 약관은 본 약관으로 대체되며, 개정된 약관의 적용일 이전 가입자도 개정된 약관의 적용을 받습니다.