결정질 암반에 위치한 가상의 방사성폐기물처분장에 대한 정확한 안전성 평가를 수행하기 위해서는 다공암반으로 대표되는 공학방벽 및 결정질 암반으로 대표되는 자연방벽을 모두 고려한 매질에서의 물질 거동 특성을 정확하게 모사하는 것이 중요하다. 본 논문에서는 단열-다공암반내 유동 및 물질전달의 정확한 해석을 위한 연결망 구축방법에 대하여 서술하였다. 실제 단열암반을 사실적으로 모사하기 위하여 단열 물성자료는 확률밀도함수를 이용하여 생성하였다. 2차원의 사각격자로 모사된 단열암반과 육면체로 모사된 다공암반간 교차선의 원활한 검색을 위하여 단열암반과 연결전 다공암반의 교차면에 가상의 단열암반을 추가적으로 도입하였다. 전체유동경로를 구성하기 저하여 생성된 단열들 간의 교차선을 효율적으로 검색할 수 있는 방법 및 단열암반 및 단열-다공암반간의 연결도를 신속히 확인하는 알고리즘을 제안하였다. 이러한 방법들은 추후 방사성폐기물처분장에서의 단열-다공 암반을 통한 핵종 이동 특성을 모사할 수 있는 수치코드 개발에 많은 도움을 줄 것으로 사료된다.

본 연구에서는 3차원 이산 균열망 수치모형을 이용하여 균열망을 구성하는 균열요소의 길이분포가 유체 흐름 특성이 미치는 영향에 대해 수치적으 로 분석하였다. 균열요소의 길이분포의 생성을 위해 절단멱분포법칙을 적용하였으며, 지수 β₁을 1.0에서부터 6.0까지 변화시키면서 유체 흐름 모의를 수행하였다. 모의결과 지수 β₁이 증가함에 따라 균열요소들의 길이분포는 점차적으로 작아지며, 이로 인해 균열망의 투수성에 영향을 미치는 균열요소들 간의 연결성은 취약해지는 것으로 나타났다. 각각의 지수 β₁에 대해 균열요소 각각에서 계산된 유량분포를 분석하였을 때 β₁= 1.0에 서의 평균유량이   6.0에 비해 약 447배 크게 산정되었으며, 균열망의 유출경계에서 계산된 유량의 경우 β₁일 때가 6.0에 비해 약 6,440배 크게 산정되었다.

균열 암반 매질에서의 지하수 흐름과 오염물질 이송에 대한 수치모의 실험이 hydromechanic 모형과 추계적 그리고 이산적 3차원 균열망 모형에 바탕을 둔 비정상상태 흐름 수치 모형을 이용하여 수행되었다. 오염물질 이송에 대한 수치모의 실험에서 random walk의 일종인 particle following 알고리즘이 사용되었다. 이 연구의 목적은 지하 깊은 곳에 위치한 Hot dry rock에서의 지열 개발을 위해 프랑스 Soultz sous Fo