절리의 방향분포가 이차원 DFN 시스템의 수리적 특성에 미치는 영향

Abstract

본 연구는 절리의 방향분포가 절리성 암반의 수리적 특성에 미치는 영향을 평가하기 위하여 등가파이프 연결구조에 기반을 둔 이차원 DFN(discrete fracture network) 유체유동 해석 프로그램 코드를 개발하고 수치실험을 수행하였다. 수치실험에 사용된 이차원 DFN 시스템은 경계효과를 고려하여 <TEX>$32m{\times}32m$</TEX> 영역의 중앙에서 <TEX>$20m{\times}20m$</TEX> 크기의 DFN 블록이다. 두 절리군을 사용하여 절리의 빈도와 길이분포를 고정하고 절리군 선주향의 평균과 표준편차를 달리하며 추계론적으로 생성한 총 15개의 DFN 블록에 대하여 매 <TEX>$30^{\circ}$</TEX> 간격으로 회전하면서 총 12 방향으로 구현한 총 180개 DFN 블록에서 블록수리전도도가 산정되었다. 또한, 각각의 DFN 블록에서 이론적 블록수리전도도, 주 수리전도도텐서 및 평균블록수리전도도를 추정하여 비교분석한 결과, 절리군의 평균 교차각이 작을수록 절리 방향분포의 변동성이 이차원 DFN 시스템의 등가연속체 취급 가능성 및 블록수리전도 특성에 더욱 영향을 미치는 것으로 평가되었다. A program code was developed to calculate block hydraulic conductivity of the 2-D DFN(discrete fracture network) system based on equivalent pipe network, and implemented to examine the effect of joint orientation distribution on the hydraulic characteristics of fractured rock masses through numerical experiments. A rock block of size <TEX>$32m{\times}32m$</TEX> was used to generate the DFN systems using two joint sets with fixed input parameters of joint frequency and gamma distributed joint size, and various normal distributed joint trend. DFN blocks of size <TEX>$20m{\times}20m$</TEX> were selected from center of the <TEX>$32m{\times}32m$</TEX> blocks to avoid boundary effect. Twelve fluid flow directions were chosen every <TEX>$30^{\circ}$</TEX> starting at <TEX>$0^{\circ}$</TEX>. The directional block conductivity including the theoretical block conductivity, principal conductivity tensor and average block conductivity were estimated for generated 180 2-D DFN blocks. The effect of joint orientation distribution on block hydraulic conductivity and chance for the equivalent continuum behavior of the 2-D DFN system were found to increase with the decrease of mean intersection angle of the two joint sets. The effect of variability of joint orientation on block hydraulic conductivity could not be ignored for the DFN having low intersection angle between two joint sets.