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

단일절리에서 2상유체의 거동을 모의하기 위해 개발된 수치모형의 검증을 위해서 상대투과계수 특성식을 수치모형에 적용하여 가스와 물의 동시거동을 해석한 후, 수치모의 결과를 모형실험결과와 비교하였다. 절리면의 거칠기와 간극의 크기에 댸한 민감도 분석을 실시한 결과, 가스의 이동속도는 절리면의 거칠기와 반비례하였으며, 절리간극의 크기와는 상대투과계수 특성식의 영향으로 단상유체의 흐름에서와 같은 간극크기의 제곱에 비례하는 향상은 보이지 않았다. 수치모형의 현장

본 연구에서는 단일절리에서 2상유체 동시거동을 해석하기 위해서 2차원 유한차분 수치모형을 개발하였다. 개발된 모형은 압력에 따른 점성의 변화가 포화도에 따른 상대투과계수의 변화를 절리간극의 크기별로 고려할 수 있다. 수치기법으로는 IMPES해법을 적용하여 물과 가스의 압력변화량과 포화도를 차례로 구하였다. 개발된 수치 모형에 이용할 상대투과계수의 특성식 도출을 위해서 일곱가지 경우의 평판모형실험을 실시하였다. 실험으로부터 도출된 상대투과계수 특성곡선은

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