• 터널과지하공간
• /
• 제13권6호
• /
• pp.465-475
• /
• 2003

본 논문은 지하공동이 건설되는 지역을 대상으로 공동내 지하수 유입량을 계산하기 위하여 영국의 AEA Technology에서 개발한 NAPSAC과 NAMMU를 이용하여 수치 모델링을 실시하였다. 이것은 3차원의 파쇄된 매질에서 지하수 유동에 대하여 아주 유효한 유한요소 소프트웨어 패키지이다. 모델링에 사용된 입력자료는 지표조사, 시추조사, 단일공 수리시험 자료 등을 이용하여 산출하였다. 공동내 지하수 유입량을 보다 정확히 산출하기 위하여 이방성 수리전도도를 고려하였으며 이방성 수리전도도값은 연구지역에 분포하는 균열체계의 확률적 분포특성이 반영된 균열망에서 역산되었다. 연속체 개념의 수치모델링을 통하여 평수기, 갈수기와 풍수기시의 지하 공동내 유입량, 유동시간, 유동경로의 변화를 계산하였다.

• 지질공학
• /
• 제26권1호
• /
• pp.87-99
• /
• 2016

본 연구에서는 단열을 포함하는 불포화 암반에서의 지하수 유동 흐름을 예측하기 위한 2차원 수치 모델링을 수행하였다. 특히, 불연속 단열망 모델링(Dscrete Fractured Network, DFN)을 통하여 도출된 투수계수의 분포를 나타내는 k-field를 모델링 입력변수로 적용하였다. 투수계수 차이는 불포화대에서 지하수의 이동속도 차이를 야기시키는 중요한 인자로 적용되었다. 연구지역의 지표 토양층으로부터 지하 대수층까지의 불포화대 깊이를 적용하여 지하수 유동 모델링의 초기조건을 실제와 유사하게 설정하였다. 강우의 지표 침투율은 인공구조물과 자연 토양층의 침투율 차이를 적용하여 실제 불포화 암반을 거쳐서 포화대까지 지하수의 이동흐름을 해석하고 시각화하였다. 특히, 오염물질의 이동 시작점이 될 수 있는 인공구조물의 하부에 모니터링 지점을 설정하여 실제 오염물질원이 지하수에 용해되어 불포화대를 이동할 때의 경로를 예측 하였으며 중력방향인 직하부로 이동하는 것을 확인되었다.

• Geomechanics and Engineering
• /
• 제16권1호
• /
• pp.21-34
• /
• 2018

Discontinuities considerably affect the mechanical and hydraulic properties of rock mass. These properties of the rock mass are influenced by the geometry of the discontinuities to a great extent. This paper aims to render an account of the geometrical parameters of several discontinuity sets related to the surrounding rock mass of Rudbar Lorestan Pumped Storage Power Plant powerhouse cavern making use of the linear and areal (circular and rectangular) sampling methods. Taking into consideration quite a large quantity of scanline and the window samplings used in this research, it was realized that the areal sampling methods are more time consuming and cost-effective than the linear methods. Having corrected the biases of the geometrical properties of the discontinuities, density (areal and volumetric) as well as the linear, areal and volumetric intensity accompanied by the other properties related to four sets of discontinuities were computed. There is an acceptable difference among the mean trace lengths measured using two linear and areal methods for the two joint sets. A 3D discrete fracture network generation code (3DFAM) has been developed to model the fracture network based on the mapped data. The code has been validated on the basis of numerous geometrical characteristics computed by use of the linear, areal sampling methods and volumetric method. Results of the linear sampling method have significant variations. So, the areal and volumetric methods are more efficient than the linear method and they are more appropriate for validation of 3D DFN (Discrete Fracture Network) codes.

• 한국수자원학회논문집
• /
• 제31권3호
• /
• pp.347-360
• /
• 1998

균열 암반 매질에서의 지하수 흐름과 오염물질 이송에 대한 수치모의 실험이 hydromechanic 모형과 추계적 그리고 이산적 3차원 균열망 모형에 바탕을 둔 비정상상태 흐름 수치 모형을 이용하여 수행되었다. 오염물질 이송에 대한 수치모의 실험에서 random walk의 일종인 particle following 알고리즘이 사용되었다. 이 연구의 목적은 지하 깊은 곳에 위치한 Hot dry rock에서의 지열 개발을 위해 프랑스 Soultz sous Foret 지역에 설치된 두개의 깊은 착정인 GPK1과 GPK2 사이에서의 tracer test 반응을 1995년에 실행된 유체순환 현장 실험으로 부터 얻어진 자료를 이용하여 예측하는 것이다. 모의 실험 결과 비반응입자(nonreactive particles)에 대한 평균 이송시간은 두 착정 사이에서 약 5일이었다.

• 한국방사성폐기물학회:학술대회논문집
• /
• 한국방사성폐기물학회 2009년도 학술논문요약집
• /
• pp.187-187
• /
• 2009

Opening of fractures induced by shear dilation or normal deformation can be a significant source of fracture permeability change in fractured rock, which is important for the performance assessment of geological repositories for spent nuclear fuel. As the repository generates heat and later cools the fluid-carrying ability of the rocks becomes a dynamic variable during the lifespan of the repository. Heating causes expansion of the rock close to the repository and, at the same time, contraction close to the surface. During the cooling phase of the repository, the opposite takes place. Heating and cooling together with the, virgin stress can induce shear dilation of fractures and deformation zones and change the flow field around the repository. The objectives of this work are to examine the contribution of thermal stress to the shear slip of fracture in mid- and far-field around a KBS-3 type of repository and to investigate the effect of evolution of stress on the rock mass permeability. In the first part of this study, zones of fracture shear slip were examined by conducting a three-dimensional, thermo-mechanical analysis of a spent fuel repository model in the size of 2 km $\times$ 2 km $\times$ 800 m. Stress evolutions of importance for fracture shear slip are: (1) comparatively high horizontal compressive thermal stress at the repository level, (2) generation of vertical tensile thermal stress right above the repository, (3) horizontal tensile stress near the surface, which can induce tensile failure, and generation of shear stresses at the comers of the repository. In the second part of the study, fracture data from Forsmark, Sweden is used to establish fracture network models (DFN). Stress paths obtained from the thermo-mechanical analysis were used as boundary conditions in DFN-DEM (Discrete Element Method) analysis of six DFN models at the repository level. Increases of permeability up to a factor of four were observed during thermal loading history and shear dilation of fractures was not recovered after cooling of the repository. An understanding of the stress path and potential areas of slip induced shear dilation and related permeability changes during the lifetime of a repository for spent nuclear fuel is of utmost importance for analysing long-term safety. The result of this study will assist in identifying critical areas around a repository where fracture shear slip is likely to develop. The presentation also includes a brief introduction to the ongoing site investigation on two candidate sites for geological repository in Sweden.

• 터널과지하공간
• /
• 제17권2호
• /
• pp.109-118
• /
• 2007

굴착손상영역(EDZ)은 굴착으로 인해 현지 암반이 역학적으로 손상을 입게 되어 응력상태, 변위상태, 암반의 안정성, 지하수의 흐름상태 등에 변화가 일어나는 영역을 의미한다. EDZ의 역학적 특성과 관련한 많은 연구들이 수행되었지만, EDZ에서의 지하수 유동 특성에 관한 연구는 아직 부족한 수준이다. 본 연구에서는 굴착으로 인한 수리-역학 상호작용(coupling)에 의해 굴착면 주변의 수리적 간극값이 변하는 영역을 수리적 굴착 thstkdduddudr이라 정의하고, 이를 3차원 분리단열망(DFN)에 적용시켜 보았다. 이를 통해 수리적 간극변화가 3차원 불연속 망에서의 전반적인 지하수 유동에 미치는 영향을 파악하였다. 또한 3차원 DFN 지하수 유동 해석 시 주로 이용되는 수두 조건과 유량 조건의 적용성을 고찰하였다. 해석 결과 수리-역학적 상호거동에 의해 발생하는 굴착면 주변의 수리적 간극변화는 터널 내부로 유입되는 유량에 큰 영향을 미치는 것으로 나타났다. 또한 DFN 해석 시 다양한 경계조건에 따른 상이한 결과를 토대로 보다 합리적인 경계조건 설정에 대한 방향을 제시하였다. 마지막으로 실제 현장에서 수리해석을 실시한 자료를 바탕으로 수리적 간극 값의 변화를 고려할 때와 고려하지 않을 때의 유입유량 차이를 통해, 3차원 지하공동의 지하수 유동해석 시 수리적 간극 값의 변화를 고려하는 것이 보다 더 보수적인 결과를 나타내는 것을 확인하였다.

• 지질공학
• /
• 제26권3호
• /
• pp.325-330
• /
• 2016

현재 경주 중저준위방사성폐기물 처분장(이하 '경주방폐장')에서는 2단계 표층처분시설이 계획중에 있으며, 포화대에 위치한 1단계 처분시설과는 달리 불포화대 상부에 위치하게 된다. 단열을 포함하는 불포화대의 특성상 지하수 및 용질의 대부분이 단열을 통해 이동할 것으로 예상된다. 따라서 불포화 암반 매질에 대한 정밀한 해석을 위하여 단열망 연속체와 암반 매질 연속체를 구분하여 해석하는 TOUGH2 전산코드의 meshmaker 모듈의 MINC 기법을 활용하였다. TOUGH2 MINC 기법의 기존 국내 연구 사례가 미미하여 본 연구에서는 MINC를 이용한 mesh 구성 방법에 대한 절차를 개발하였으며, 단열 연속체와 암반매질 연속체의 k-field를 생성하였다. 이와 같이 생성된 도메인은 향후 이중 연속체를 기반으로 경주방폐장의 지하수 유동 및 오염물질 이동 등에 활용될 뿐만 아니라 단열이 발달한 암반에서 단열-암반매질 연결성을 고려한 단열망 유동 특성을 분석하는데 참고가 될 것으로 기대한다.

• 방사성폐기물학회지
• /
• 제6권4호
• /
• pp.265-282
• /
• 2008

중 저준위 방사성폐기물 처분 부지의 부지특성조사 결과를 이용하여 처분장 부지의 지하수 유동체계를 이해하기 위한 수치 모델링을 수행하였다. 부지의 투수성 단열대 및 암반단열의 분포 특성에 근거하여 단열망 모델을 구축하고, 이를 이용하여 생성된 10개의 수리전도도장을 지하수 유동 모델링에 반영한 추계론적 Hybrid-EPM 방법으로 수치 모델을 구성하였다. 10회의 지하수 유동 모델링 결과, 처분 부지의 지하수두 및 지하수 흐름은 지표 근처에서 지형적인 요소에 크게 지배를 받는 것으로 나타나며 처분장 심도에서는 주변에 존재하는 투수성이 높은 단열대에 의해 영향을 받음을 확인할 수 있었다. 특히, 처분 시설 건설 중 사일로 주변 지역에서 수위 강하가 크게 발생하는 것으로 분석되었다. 처분 시설 폐쇄한 후 지하수위는 1년 이내에 급속히 회복되며, 대략 2년이 지난 후 완전히 회복 될 것으로 분석되었다.

• 한국터널공학회:학술대회논문집
• /
• 한국터널공학회 2005년도 학술발표회 논문집
• /
• pp.348-358
• /
• 2005

Although the evaluation of the mechanical properties and behavior of discontinuous rock masses is very important for the design of underground openings, it has always been considered the most difficult problem. One of the difficulties in describing the rock mass behavior is assigning the appropriate constitutive model. This limitation may be overcome with the progress in discrete element software such as PFC, which does not need the user to prescribe a constitutive model for rock mass. Instead, the micro-scale properties of the intact rock and joints are defined and the macro-scale response results from those properties and the geometry of the problem. In this paper, a $30m{\times}30m{\times}30m$ jointed rock mass of road tunnel site was analyzed. A discrete fracture network was developed from the joint geometry obtained from core logging and surface survey. Using the discontinuities geometry from the DFN model, PFC simulations were carried out, starting with the intact rock and systematically adding the joints and the stress-strain response was recorded for each case. With the stress-strain response curves, the mechanical properties of discontinuous rock masses were determined and compared to the results of empirical methods such as RMR, Q and GSI. The values of Young's modulus, Poisson's ratio and peak strength are almost similar from PFC model and Empirical methods. As expected, the presence of joints had a pronounced effect on mechanical properties of the rock mass. More importantly, the mechanical response of the PFC model was not determined by a user specified constitutive model.

• 한국방사성폐기물학회:학술대회논문집
• /
• 한국방사성폐기물학회 2005년도 춘계 학술대회
• /
• pp.231-238
• /
• 2005

본 연구에서는 유동저항이론을 기초로 하여 연속다공체, 분리단열망 및 연속다공체-분리단열망 공존암반과 같은 3가지 암반을 대상으로 암반 특성에 따른 지하수 유동저항 개념 모델링 및 관계식을 제안하였다. 정상상태조건에서 밀도변동은 고려하지 않았으며 유한 체적법을 이용하였다. 각종 물성치는 블록 중심에서 정의되고, flux는 블록면에서 정의되는 staggered 격자 체계하에서 모든 블록에 대해 Darcy 법칙이 적용되었다. 접촉면에서의 투수계수는 인접면 중심에서 정의된 물성치의 조화평균값을 사용하였다. 유동저항개념을 이용하여 인접한 블록간의 상대압력차와 flux의 관계를 표현하였다. 개개의 단열에서의 유동은 다공암반에서 이용된 방정식과 동일한 형태의 2차원 방정식으로 모사되었다. 본 논문에서 제안된 모델은 추후 다양한 암반 특성별 유동 모사 기법을 개발하는데 많은 기여를 할 것으로 기대된다.