사면 내의 지하수 유동과 사면의 안정성에 대한 강수 영향의 완전 연동된 수리지질역학적 수치 해석

A Fully Coupled Hydrogeomechanical Numerical Analysis of Rainfall Impacts on Groundwater Flow in Slopes and Slope Stability

사면 내의 지하수 유동과 사면의 안정성에 대한 강수의 영향을 평가하기 위하여 수리지질역학적 수치 모델이 제시되었다. 이 수치 모델은 변형성 지질매체 내에서의 포화-불포화 지하수 유동을 설명하는 완전 연동된 간극탄성론적 지배 방정식들과 Galerkin 유한요소법에 근거하여 개발되었다. 이렇게 개발된 완전 연동된 수리지질역학적 수치 모델은 다양한 강수량 조건 하에 있는 불포화 사면에 대한 일련의 수치모의실험에 적용되었다. 이러한 수치모의실험 결과들은 강수량이 증가할수록 사면의 전반적인 수리역학적 안정성이 저하되며 잠재적인 파괴가 사면 전단부에서부터 시작되어 사면 정상부 쪽으로 팽창됨을 보여주고 있다. 강수량이 증가함에 따라 수리지질학적으로는 압력수두와 전체 수리수두가 증가한다. 그 결과 지하수면이 상승하고 불포화대가 감소하며 삼출면이 사면 전단부로부터 사면 정상부쪽으로 팽창하고 이러한 삼출면에서의 지하수 유동 속도도 증가하게 된다. 한편 강수량이 증가함에 따라 지질역학적으로는 사면 전단부를 향해 수평 변위는 증가하지만 수직 변위는 감소하게 된다. 이는 강수량 증가에 따른 지하수면의 상승에 수반하는 부력에 기인하는 것으로 해석된다. 그 결과 사면의 전반적인 변형이 사면 전단부를 향해 심화되고 전단파괴 안전율이 1 이하인 불안전한 지역이 사면 전단부에서 두꺼워지면서 사면 전단부에서부터 사면 정상부 쪽으로 팽창하게 된다. 또한 수치모의실험 결과들은 이러한 잠재적인 전단파괴면과 지표면 사이의 지반에서는 잠재적인 인장파괴가 발생할 수 있음을 보여주고 있다.

A hydrogeomechanical numerical model is presented to evaluate rainfall impacts on groundwater flow in slopes and slope stability. This numerical model is developed based on the fully coupled poroelastic governing equations for groundwater flow in deforming variably saturated geologic media and the Galerkin finite element method. A series of numerical experiments using the model developed are then applied to an unsaturated slope under various rainfall rates. The numerical simulation results show that the overall hydromechanical slope stability deteriorates, and the potential failure nay initiate from the slope toe and propagate toward the slope crest as the rainfall rate increases. From the viewpoint of hydrogeology, the pressure head and hence the total hydraulic head increase as the rainfall rate increases. As a result, the groundwater table rises, the unsaturated zone reduces, the seepage face expands from the slope toe toward the slope crest, and the groundwater flow velocity increases along the seepage face. From the viewpoint of geomechanics, the horizontal displacement increases, and the vertical displacement decreases toward the slope toe as the rainfall rate increases. This may result from the buoyancy effect associated with the groundwater table rise as the rainfall rate increases. As a result, the overall deformation intensifies toward the slope toe, and the unstable zone, in which the factor of safety against shear failure is less than 1, becomes thicker near the slope toe and propagates from the slope toe toward the slope crest. The numerical simulation results also suggest that the potential tension failure is likely to occur within the slope between the potential shear failure surface and the ground surface.