[논문]수치해석을 통한 춘천 산사태지역 토석류 거동 분석

최정해

doi:10.9720/kseg.2017.1.59

[국내논문] 수치해석을 통한 춘천 산사태지역 토석류 거동 분석
Analysis of Debris Flow of Chun-cheon Landslide Area using Numerical Methods 원문보기

지질공학 = The journal of engineering geology, v.27 no.1, 2017년, pp.59 - 66

초록
AI-Helper

최근들어 발생하는 강우는 지금까지의 강우특성과 비교하여 여름철에 단기간 집중적으로 내리는 형태로 바뀌고 있다. 2011년 강원도 춘천지역에서는 이러한 집중강우로 인해 많은 산사태가 발생하였으며, 이로 인해 재산피해와 인명피해를 야기시켰다. 본 연구의 목적은 강우에 의한 산사태 물질의 이동에 대한 특성을 분석하기 위해서 현장 자료와 수치해석 자료를 활용하여 시간에 따른 특성의 변화 및 강우량에 따른 특성의 변화를 분석하는 것이다. 본 연구에서는 항복응력 값의 변화를 통해서 강우량을 재현하였으며, 이를 통해서 어떠한 거동 특성을 보이는지 해석하였다. 그 결과 많은 강우에 의해 항복응력 값이 낮아지는 경우가 빠른 속도의 이동형태를 보이고 동시에 넓은 지역에 걸쳐 확산되는 특성을 확인하였으며, 이러한 결과는 현장에서 확인된 토석류의 흔적과 매우 유사하다는 것을 확인하였다. 본 연구를 통해서 산사태 발생 당시의 토석류 최대 속도 및 최대 두께 등을 파악하는 것이 가능하다.

Abstract ▼ AI-Helper

The characteristic of recent rainfall pattern in Korea is concentrated in summer season and it is very different compare with former characteristic. In 2011, there was heavy rainfall in Chuncheon city of northern part of Korea. Because of rainfall in short time, many landslides were occurred in narrow area and many people were killed by these landslides at that time. The purpose of this study is to calculate run-out distance of debris flow and analyze the movement properties of debris flow according to the elapsed time using numerical analysis method at that time. The debris 2D program, which is developed by prof. Liu in National Taiwan University, was used in this study. Run-out distance of debris flow was calculated under different yield strength conditions which were controlled by rainfall amount. The results reveal that absolute maximum velocity of the debris flow is about 8.1 m/s and maximum depth of debris flow is about 7 m when debris flow was occurred. The run-out distance after 500 sec is about 300 m from end of the valley. It is very well similar with actual debris flow run-out distance. From these results, we can presume the maximum velocity and depth of debris flow at that time.

주제어

AI 본문요약
AI-Helper

* AI 자동 식별 결과로 적합하지 않은 문장이 있을 수 있으니, 이용에 유의하시기 바랍니다.

문제 정의

2011년 7월에는 집중된 하계기간의 집중호우로 인해 강원도 춘천시 일대를 비롯하여 경기도 여주 및 이천 등지에서 산발적인 산사태가 발생하였으며 이로 인해 가옥 및 인명피해가 발생하였다. 본 연구는 강원도 춘천시 신북읍에서 발생한 산사태에 대해서 산사태 발생 시의 토석류의 흐름에 대해 분석을 실시하고 흐름에 영향을 줄 수 있는 항복응력(yield stress)이 강우에 의해서 감소된다는 특성을 이용하여 항복응력의 변화를 이용하여 강우에 의해 확산되는 토석류의 범위를 계산하였다. 산사태가 발생한 곳은 2011년 7월 26일부터 27일 양일에 걸친 집중호우로 인해 토석류가 발생하여 건물 및 인명피해가 발생한 지역이다.
산사태의 규모는 약 220 m 정도의 연장을 가지며 정상부 하단에서 발생하여 하부의 도로까지 진행한 것으로 보여지며 토석류의 이동 경로에 있던 가옥에 피해를 입힌 것으로 파악된다. 본 연구에서는 국립대만대학교 코페이 류(Kofei Liu) 교수가 개발한 Debris-2D 프로그램을 활용하여 인위적인 산사태 물질의 부피를 가정하여 토석류의 흐름과 확산 범위 그리고 시간대별 최대 이동 속도를 계산하여 당시 발생한 토석류에 대한 특성을 해석하였다.

제안 방법

본 연구에서는 강우상태에 따라 토석류를 구성하는 물질의 전단응력이 약화되고 이런 결과로 인해 항복응력(yield stress)이 달라진다는 조건으로 항복응력을 위 표와 같이 각각 400 Pa, 700 Pa, 1300 Pa로 구분하여 분석을 실시하였다. 분석결과를 이용하여 토석류가 발생하고 80분 후까지의 토석류의 속도변화에 대한 분석을 실시하고 또한 발생이후 약 10초간의 속도변화를 각 경우에 따라 분석을 실시하였다.
본 연구에서는 대략적인 조사를 통해서 현장에 분포하는 토심을 평균 1.5 m로 산정하였으며, DEM 자료를 활용하여 토석류의 근원물질에 대한 영역을 설정하였다(Fig. 3). 이를 바탕으로 식 (10)을 이용하여 사면에 분포하는 전체 토석류의 평균심도는 물이 충분히 포화되었다는 가정 하에 2.
앞서 설명한 바와 같이 강우상태에 따라 지표에 존재하는 토사의 전단강도가 저하된다는 가정 하의 분석에서 정의된 토석류의 항복응력을 조절하여 프로그램을 실행하였으며, 이로 인한 토석류의 확산 범위를 시간대별로 분석하였다. 분석 시간은 최초 토석류가 발생하고 난 이 후 5초 후, 20초 후, 200초 후로 나누어 3개의 경과시간에 따른 토석류의 확산범위를 계산하였다. 그 결과는 Fig.
본 연구에서는 강우상태에 따라 토석류를 구성하는 물질의 전단응력이 약화되고 이런 결과로 인해 항복응력(yield stress)이 달라진다는 조건으로 항복응력을 위 표와 같이 각각 400 Pa, 700 Pa, 1300 Pa로 구분하여 분석을 실시하였다. 분석결과를 이용하여 토석류가 발생하고 80분 후까지의 토석류의 속도변화에 대한 분석을 실시하고 또한 발생이후 약 10초간의 속도변화를 각 경우에 따라 분석을 실시하였다. 그 결과는 Fig.
산사태가 발생한 시점에 대한 정확한 해석은 당시의 기록을 바탕으로 재현 할 수 있으며, 현장에 남아있는 흔적들을 바탕으로 본 연구에서는 수치해석 프로그램을 활용하여 발생 시의 토석류의 흐름에 대한 해석을 실시하였다.
앞서 설명한 바와 같이 강우상태에 따라 지표에 존재하는 토사의 전단강도가 저하된다는 가정 하의 분석에서 정의된 토석류의 항복응력을 조절하여 프로그램을 실행하였으며, 이로 인한 토석류의 확산 범위를 시간대별로 분석하였다. 분석 시간은 최초 토석류가 발생하고 난 이 후 5초 후, 20초 후, 200초 후로 나누어 3개의 경과시간에 따른 토석류의 확산범위를 계산하였다.
)이다. 이를 계산하기 위해서는 Fig. 2와 같이 거친면을 가지는 판을 준비하고 시료를 위치시킨 다음 또 다른 판을 준비하여 수직된 방향으로 천천히 밀어 시료의 수직높이 변화를 측정한다. 이러한 작업을 평균 10회 정도 실시하여 식 (11)을 사용하여 항복응력 (τ₀)을 계산한다.
49 m로 산정하였다. 이를 기본으로 DEM 좌표에서 부피를 계산하고 토석류의 흐름에 대한 Debris2D 프로그램을 실행하였다. 기본적으로 프로그램은 토석류의 속도가 0.
그러나 이번 분석에서는 같은 시간동안 입력인자의 값에 따라 흐름의 양상을 분석하는 것이 목적이기 때문에 시간을 동일하게 주어 계산이 일정시간에 멈추도록 하였다. 입력조건은 강우에 따라 변화를 보이는 항복응력을 다르게 계산하여 분석을 실시하였다. 항복응력의 조건 및 시간 설정은 Table 1과 같다.

대상 데이터

본 연구지역은 용두리 편마암 복합체가 가장 광범하게 분포하고 주로 자류석편마암(미그마타이트질), 호상 흑운모편마암 및 흑운모편암으로 구성된다. 본 연구지역에 분포하는 자류석편마암은 뚜렷한 엽리면이 없고 자류석, 흑운모, 녹니석 등이 포함된 편마암이다. 호상 흑운모편마암은 연구지역에 가장 넓게 분포하고 있으며 기저를 이루는 자류석편마암과 큰 차이를 보이지 않는다.
변성퇴적암류는 이 지역의 대부분을 차지하며 하부로부터 용두리 편마암 복합체, 구봉산층군, 의암층군으로 구분되며, 용두리 편마암 복합체와 구봉산층군 사이의 관계는 미상이다. 본 연구지역은 용두리 편마암 복합체가 가장 광범하게 분포하고 주로 자류석편마암(미그마타이트질), 호상 흑운모편마암 및 흑운모편암으로 구성된다. 본 연구지역에 분포하는 자류석편마암은 뚜렷한 엽리면이 없고 자류석, 흑운모, 녹니석 등이 포함된 편마암이다.

이론/모형

사용한 구성방정식은 매우 얇은 물의 막을 형성하고 있다는 가정을 한 질량 및 운동보존의 법칙을 기초로 한다. 프로그램에 사용한 좌표계는 일반적인 Cartesian 좌표를 활용하였으며, 전체의 토석류내에 존재하는 변형률 비(strain rate)와 전단응력(shear stress)을 표현하기 위해서 Julien and Lan (1991)이 제안한 일반화된 3차원 구성관계(constitutive relation)를 활용하였다. 이 관계는 식 (1), (2)로 표현된다.

성능/효과

5에 나타나있다. Fig. 4에서 알 수 있듯이 전체적인 토석류의 흐름은 항복응력이 낮을수록 토석류의 흐름속도는 높게 나타나는 것을 확인하였다. 이는 점착력이 약해지면서 전단응력이 낮아지고 이로 인해 점착력과 전단응력의 영향을 받는 항복응력이 낮아져 발생하는 결과이다.
이는 점착력이 약해지면서 전단응력이 낮아지고 이로 인해 점착력과 전단응력의 영향을 받는 항복응력이 낮아져 발생하는 결과이다. Fig. 5를 통해서 각 경우에 대하여 최대 토석류 흐름 속도를 비교 분석한 결과 항복응력이 400 Pa인 경우는 최대 흐름 속도가 8.10 m/s이며 700 Pa인 경우는 최대 흐름 속도가 8.04 m/s, 1300 Pa인 경우는 최대 흐름 속도가 7.92 m/s인 것으로 나타났다.
이렇게 얕아지는 것은 계곡의 종점부까지 진행한 토석류가 넓은 지형을 만나 확산되면서 더 이상 두꺼운 퇴적물을 만들지 않기 때문이다. 그리고 약 2000초 이후부터 보이는 퇴적두께는 퇴적물이 도로를 지나 강과 만나는 지점에 쌓이면서 나타나는 것으로 확인되었다(Fig. 9~11). 이 분석을 통해서 토석류가 발생한 당시의 최대 순간 퇴적두께는 약 6.
이는 앞서 설명한 바와 같이 충분한 강우에 의해서 토사의 전단강도가 저하되면 이로 인해 항복응력이 저하되고 그 결과로 토석류 물질은 보다 넓은 지역으로 확산된다는 것을 확인하였다. 또한 이 결과를 통해서 좌측의 토석류는 토석류 발생 약 5초 경과 후 시설물까지 확산이 되는 것을 알 수 있으며, 우측의 토석류 경우에는 발생이후 약 20초 정도가 경과된 후에 시설물에 미치는 것을 알 수 있다. 이는 토석류가 발생한 당시 토석류가 시간에 따라 어떻게 진행되었으며, 시설물까지 도달한 시간을 대략적으로 유추하는데 큰 도움을 줄 수 있다.
6에 나타나 있다. 본 계산결과 항복응력이 낮은 경우에는 확산되는 범위는 넓게 분포하고 있으며, 이와 반대로 항복응력이 높은 경우에는 확산되는 범위가 낮은 경우에 비해서 좁게 분포한다는 것을 확인하였다. Fig.
전단강도의 저하로 인해 흙이 가지는 전단강도에 대한 항복응력 값이 저하되고 이로 인해 결국은 토석류의 흐름에 영향을 미치게 된다. 앞서 분석한 바와 같이 항복응력이 낮은 경우는 높은 경우에 비해 빨리 이동하는 특성을 확인하였으며, 또한 넓은 지역에 걸쳐 확산된다는 것을 확인하였다. 이는 강우량이 많아지고 토층사이로의 침투가 빨라 흙이 충분이 포화된다는 가정 하에 토석류 물질은 매우 빠른 속도로 이동하고 넓은 범위로 분포하게 된다.
앞서 확산범위해석에서 강우량에 의해 항복응력 값이 작아지면 토석류의 흐름은 빨라지고 넓게 분포하며, 이 경우 좌측의 토석류 물질이 시설물에 도착하는 시간이 약 발생이후 5초 경과 후에 도착하는 것으로 나타났으며, 우측 토석류 물질은 약 20초 이후에 도착하는 것으로 나타났다.
9~11). 이 분석을 통해서 토석류가 발생한 당시의 최대 순간 퇴적두께는 약 6.5 m~7.5 m 정도를 보였을 것으로 추정한다. 이는 현장에서 확인한 토석류의 흔적과도 매우 유사한 결과를 보인다.
6에 나타난 바와 같이 항복응력이 400 Pa인 경우에 비해 항복응력이 1300 Pa인 경우는 토석류가 좁은 형태를 보이며 강 까지 흘러가는 것을 확인 할 수 있다. 이는 앞서 설명한 바와 같이 충분한 강우에 의해서 토사의 전단강도가 저하되면 이로 인해 항복응력이 저하되고 그 결과로 토석류 물질은 보다 넓은 지역으로 확산된다는 것을 확인하였다. 또한 이 결과를 통해서 좌측의 토석류는 토석류 발생 약 5초 경과 후 시설물까지 확산이 되는 것을 알 수 있으며, 우측의 토석류 경우에는 발생이후 약 20초 정도가 경과된 후에 시설물에 미치는 것을 알 수 있다.
앞서 실시한 수치해석의 결과는 현장에서 확인되는 토석류의 흔적과 매우 유사한 결과를 보인다. 최대 퇴적 높이뿐만 아니라 최종적으로 확산된 토석류의 범위 또한 실제 현장의 흔적들과 매우 유사하다는 것을 확인하였다. 토석류의 특성상 실시간 해석이 매우 어렵다는 점을 고려하였을 때, 이러한 수치해석을 통한 토석류 발생 시의 특성을 해석하는 것은 토석류를 이해하는 데 매우 큰 도움을 줄 것으로 생각한다.

질의응답

핵심어	질문	논문에서 추출한 답변
	여름철에 단기간 집중적으로 내린 강우의 예는?	최근들어 발생하는 강우는 지금까지의 강우특성과 비교하여 여름철에 단기간 집중적으로 내리는 형태로 바뀌고 있다. 2011년 강원도 춘천지역에서는 이러한 집중강우로 인해 많은 산사태가 발생하였으며, 이로 인해 재산피해와 인명피해를 야기시켰다. 본 연구의 목적은 강우에 의한 산사태 물질의 이동에 대한 특성을 분석하기 위해서 현장 자료와 수치해석 자료를 활용하여 시간에 따른 특성의 변화 및 강우량에 따른 특성의 변화를 분석하는 것이다.
	최근의 강우특성은?	최근의 강우특성은 예전의 강우특성과 비교하여 하계기간에 집중되는 양상을 보이고 있다. 특히 최근 들어 이러한 하계기간의 집중호우로 인해 많은 건물 및 인명피해를 입은 사례가 증가하고 있다.
	토석류의 흐름 및 확산범위를 계산하기 위한 입력요소로 필요한 것은?	프로그램을 사용하여 토석류의 흐름 및 확산범위를 계산하기 위해서 필요한 입력요소(input parameter)는 현장에서 조사한 토층의 깊이를 포함한 토석류의 부피를 계산할 수 있는 기초자료와 흐름을 해석하기 위한 DEM 및 항공사진 자료가 반드시 필요하다. 현장조사 및 항공사진을 통해서 산사태 발생지역에 대한 예상 건조부피(Vd)를 계산 할 수 있다.

참고문헌 (13)

Julien, P. Y. and Lan, Y., 1991, Rheology of hyperconcentrations, Journal of Hydraulic Engineering ASCE, 117(3), 346-353.

상세보기
Liu, K. F. and Huang, M. C., 2006, Numerical simulation of debris flow with application on hazard area mapping, Computational Geoscience, 10(2), 221-240.

상세보기
Hanson, G. J. and Simon, A., 2001, Erodibility of cohesive streambeds in the loess area of the midwestern USA, Hydrological Processes, 15(1), 23-38.

상세보기
Hanson, G. J. and Cook, K. R., 2004, Apparatus, test procedures, and analytical methods to measure soil erodibility in situ, Applied Engineering in Agriculture, 20(4), 455-462.

상세보기
Nam, S., Petrie, J., Diplas, P., and Gutierrez, M., 2010, Effects of spatial variability on the estimation of erosion rates for cohesive riverbanks, International Conference on Fluvial Hydraulics (River Flow 2010), 8-10.
Takahashi, T., 2007, Debris flow :Mechanics, Prediction and Countermeasures, Taylor & Francis, New York, USA, 448.
Arattano, M. and Marchi, L., 2008, Systems and sensors for debris-flow monitoring and warning, Sensors 8, 2436-2452.

상세보기
Liu, K. F. and Hsu, Y. C., 2008, Study on the sensitivity of parameters relating to debris flow spread, In Proceedings of the International Conference, Debris Flows: Disasters, Risk, Forecast, Protection, Pyatigorsk, Russia.
Liu, K. F. and Lee, F. C., 1997, Experimental analysis on impact mechanism of granular flow, Chinese Journal of Mechanics, 13(1) 87-100.
Liu, K. F., Li, H. C. and Hsu, Y. C., 2009, Debris flow hazard assessment with numerical simulation, Natural Hazards, 49, 137-161.

상세보기
Liu, K. F., Wei, S. C., and Li, P. C., 2013, The influence of accumulated precipitation on debris flow hazard area, Journal of Chinese Soil and Water Conservation, 44(3), 225-233.
Liu, K. F. and Wu, Y. H., 2010, The Assessment of Debris Flow Hazard in Korea Using Debris-2D, 2010-International Symposium in Pacific Rim, Taipei, Taiwan, 820-827.
Liu, K. F. and Wu, Y. H., 2013, Introduction to Debris-2D - A Debris Flow Simulation Program, TXT-tool 3.886-1.1: In (Sassa K. eds.) ICL Landslide Teaching Tools (ISBN: 978-4-9903382-2-0), 238-246.

저자의 다른 논문 :

LOADING...

표제어: PCR

동의어: Packet Collision Rate

용어 설명 출처 목록 (6)

용어 설명: PCR은 세균 특이성이 있는 primer를 이용하여 적은 수의 세균이 있을지라도 쉽게 검출할 수 있는 유용한 방법이며, 이를 이용하여 구강 내 치면세균막이나 타액에서 직접 세균을 검출할 수 있게 되었다[8].

내보내기 구분	파일저장 인쇄 메일전송
구성항목	기본정보 상세정보 관리번호, 논문명, 저널/프로시딩명, 저자 , 발행년, 권, 호, 시작페이지, 끝페이지, 발행기관 관리번호, 논문명, 대등논문명, 저자 , 저널/프로시딩명, 발행기관, 발행년, 발행언어, 권, 호, 시작페이지, 끝페이지, ISBN, ISSN, 주제분야, 키워드, 초록(한글), 초록(영문), 저자(소속기관)
저장형식	Text(ASCII format) Excel format RefWorks Direct Export RIS format (for Reference Manager, ProCite, EndNote), Scholar's Aids, Mendeley
메일정보	받는사람 (필수) @ 보내는사람 (선택) @ 제목 내용 KISTI 검색결과 이메일 서비스
안내	총 건의 자료가 검색되었습니다. 다운받으실 자료의 인덱스를 입력하세요. (1-10,000) 검색결과의 순서대로 최대 10,000건 까지 다운로드가 가능합니다. 데이타가 많을 경우 속도가 느려질 수 있습니다.(최대 2~3분 소요) 다운로드 파일은 UTF-8 형태로 저장됩니다. 파일의 내용이 제대로 보이지 않을실 때는 웹브라우저 상단의 보기 -> 인코딩 -> 자동선택 여부를 확인하십시오. ~ Text(ASCII format) Excel format

연합인증