Article information

Abstract

Trans Abstract

1. 서 론

최근 지구온난화의 영향으로 태풍의 강도가 증가하여 폭풍해일에 의한 해안도시의 범람피해도 증가하고 있다. 태풍에 의해 폭풍해일이 발생하면 해수면이 상승하여 해안지역을 직접 범람시키는 것이 일반적이지만, 태풍에 동반된 강우에 의한 우수가 해수면 상승으로 바다로 원활하게 배제되지 못해 범람이 가중되기도 한다. 우리나라의 경우, 2003년 태풍 Maemi로 인해 132명의 인명피해, 약 5조원의 재산피해가 발생한 것으로 파악되었다 (Kim et al., 2010). 이로 인해 연안 재해의 주요한 인자인 폭풍해일의 위험성에 대한 인식이 높아짐과 동시에, 이에 대한 연구도 활발히 수행되고 있다. 우리나라에 가장 큰 피해를 입힌 태풍 Maemi에 대한 연구로는 Kang(2004), Kang et al.(2004), Choi et al.(2004), Kawai et al. (2004, 2005), Hur et al.(2006a, 2006b), Kim et al.(2007), Chun et al.(2008)이 있다. 이러한 연구들은 태풍에 의한 폭풍해일만을 고려하여 범람을 수치모의 하였으므로 태풍에 동반되는 강우에 의한 침수피해 증가는 고려하지 않았다. 일반적으로 태풍은 강우를 동반하므로 폭풍해일에 의한 해안 도시범람을 수치모의 하는데 있어 강우의 영향이 반드시 고려되어야 한다. 또한 강우에 의한 우수 배제를 위해 우수관망이 설치되어 있는 것이 일반적이므로 우수관망 시스템이 수치모의에 고려되어야 타당한 결과를 얻을 수 있을 것으로 판단된다. 이와 같은 강우를 고려한 폭풍해일의 국내 연구로는 Park et al.(2009)이 있다. 그러나 Park et al.(2009)은 폭풍해일을 먼저 계산하고, 그 결과를 순차적으로 하천과 지상부 범람 수치모의에 이용하여 해일-하천-강우-우수관망의 동시 거동에 따른 영향을 배제한 것으로 판단된다.

본 연구에서는 폭풍해일에 의한 해안 도시범람을 수치모의 하는 데 있어 강우의 영향을 고려하기 위해 강우강도를 고려하여 우수관망 수치모형을 개발하고, 이를 폭풍해일 수치모형과 결합하였으며, 이를 단순화된 지형에 적용하여 누적강우량, 우수의 배수 및 해수의 상승에 따른 역류현상 등 결합 수치모형의 기본적인 기능을 검토하였다.

2. 지배방정식

2.2 우수관망 수치모형

우수관망 수치모형은 Yoon et al.(1997a, 1997b)이 개발한 관수로에서의 비압축성 부정류 해석 모형을 이용하였다. 이 수치모형은 관수로에 적용되는 모형이므로 유량이 작아 개수로로 흐르는 강우 초기에는 적용할 수 없으나, 강우가 어느정도 진행되어 본격적으로 흐름이 발생되면 곧 관수로로 전환되므로 실용적인 문제에 적용하기에는 무리가 없을 것으로 판단된다. 배수암거 부정류해석 연구는 비정상 흐름(unsteady flow)으로 Fig. 1과 같이 여러 개의 수조와 수조사이를 연결하는 관로로 구성된 배수계통에 대해 관로의 면적이 일정한 경우 다음의 비정상 에너지방정식 (4)와 연속방정식 (5)을 적용하여 해석할 수 있다 (Yoon et al., 1997a).

Fig. 1.

Conceptual diagram of sewer and pipe conduit.

여기서, ζ는 기준면으로부터 수면까지의 높이(m), v는 관로의 단면평균유속(m/s), Q는 관로의 유량(m³/s)이며, A_c는 수조의 단면적(m²), A는 관로의 단면적(m²), l은 인접수조간의관로길이(m), d는 관로의 직경(m), ΔQ는 관정으로부터의 순유출유량(m³/s), g는 중력가속도(m/s²), f는 관로의 마찰계수, k는 관로의 미소손실계수이다. ΔH는 인접 관정간의 속도수두를 제외한 수두차로서 다음과 같다.

여기서, ΔP_c는 인접 관정의 공기압 차(N/m²), Δζ는 인접수조사이의 수위차(m), γ_w는 물의 단위중량(N/m³)이다. 모든 관정이 개방되어 있는 경우 관정 내부에서의 수면은 대기압의 영향에 있으므로 압력수두차 P_c는 사라진다. 식 (4)는 유속의 함수이므로 수치해석의 편의를 위해 연속방정식에서 요구하는 유량변수로 다시 나타내면 다음 식 (7)과 같이 나타낼 수 있다.

위의 식 (7)에서 사용된 계수 K는 식 (4)를 이용하여 다음과 같이 표현할 수 있다.

이를 통해 개발된 우수관망 수치모형은 Fig. 2와 같이 비선형천수방정식의 연속방정식과 운동방정식 사이에 적용하여 수치모의를 수행하도록 구성하였다.

Fig. 2.

Simulation flow chart of the coupled surge-rainfall-sewer numerical model.

3. 우수관망 수치모의

3.1 우수관망에서의 배수 및 역류현상 해석 알고리즘

지상과 우수관정에서의 수위는 각각 폭풍해일 모형과 우수관망 모형의 연속방정식을 통해 산정된다. 이후 Fig. 2의 ‘Sewer System Drainage Process’를 통해 지상과 우수관정의 수위차에 의해 지상으로부터 우수관정으로의 유입유량이 결정되고, 이에 따라 지상과 우수관정의 수위가 조정된다. 강우 또는 해일의 범람으로 인해 지상의 수위가 우수관정의 수위보다 높아진 경우 지상에 고인 물이 우수관으로 유입되어 바다로 배수된다. 반면 해수면의 높이가 지상의 수위보다 높은 경우에는 해수가 우수관을 통해 역류하게 된다. Fig. 3과 Fig. 4는 각각 우수의 배수 및 해수의 역류현상이 발생할 때 나타날 수 있는 경우들을 보여준다. Fig. 3은 우수의 배수현상에 따른 우수관정의 수위변화를 나타낸 것으로서, 지상의 우수 용기는 폭풍해일 모형의 한 격자에 해당하며, 그 격자 중앙 하부에는 관정이 연결되어 있다. 그러나 그림에서는 이해를 돕기 위해 지상의 격자와 관정을 분리하여 도시하였다. Fig. 3의 Case 1은 지상에 고인 우수의 수위가 우수관정의 수위보다 높은 경우로서 Δt시간 동안 관정으로 유입되는 우수량(오리피스 공식으로 계산)이 관정을 충분히 채울 수 있을 정도로 많이 유입되어 지상의 수위와 관정의 수위가 같아지는 경우를 보여준다. Case 2는 Δt시간 동안 유입되는 우수의 양이 작아 관정을 충분히 채울 수 없으며, 지상의 수위는 하강 후에도 바닥보다 높은 상태를 유지하는 경우를 나타낸 것이다. Case 3은 지상에 고인 우수량이 적어 지상에 고인 모든 우수가 우수관으로 유입되어도 우수관정의 수위가 크게 상승하지 않는 경우이다. Case 1~Case 3의 경우 모두 우수관으로 유입된 우수에 의해 우수관정의 수위가 조정되는 경우이다. 반면 Case 4는 지상에 우수가 존재하지 않는 경우이다. 이러한 경우 관정으로의 유입이 발생하지 않아 지상 및 우수관정 수위에 영향을 미치지 않으므로 폭풍해일 모형과 우수관망 모형에서 각각 산정된 수위를 그대로 유지하게 된다.

Fig. 3.

Drainage process for storm sewer system in the cases of surface water elevation higher than junction water elevation.

Fig. 4.

Back flow process for storm sewer system in the cases of surface water elevation lower than junction water elevation.

Fig. 4는 해수의 역류현상에 따른 우수관정의 수위변화를 나타낸 것이다. Case 5는 관정의 수위가 관정이 위치하는 지상의 수위보다 높아 해수의 역류가 발생하는 경우이다. 이러한 경우 우수관정과 지상의 수위가 다시 산정된다. 반면 Case 7의 경우 관정의 수위와 지상의 수위가 동일한 경우이다. 이러한 경우 지상과 우수관정에서 수위차가 존재하지 않기 때문에 기존에 산정된 수위를 그대로 유지하게 된다.

3.2 수치모의 대상영역 및 우수관망 배치구조

본 연구에서 개발된 우수관망 수치모형의 강우에 의한 침수 과정이 잘 작동하는 지를 확인하기 위해 Fig. 5와 같이 하나의 단순한 영역을 설정하였다. 격자간격 Δx와 Δy는 10 m로 동일하게 구성하였으며 Δt는 0.1 s로 하였다. 그림에서 x축과 y축의 숫자는 격자번호를 나타내며, 2점간의 실제거리는 사각형 box내에 표시하였다. 지상영역과 바다 사이에 짙은 회색으로 표시된 부분은 방조제를 나타낸 것으로서 지상 영역에 쌓인 우수가 우수관망 시스템을 통하지 않고 바다로 유출되는 것을 방지하기 위해 설정한 것이다. 또한 영역 안에 존재하는 물이 밖으로 새어나가지 않도록 외곽 경계에 벽을 설치하였다.

Fig. 5.

Storm sewer system designed for numerical simulation test.

우수관망 시스템은 지상영역에서 강우가 관정을 통해 신속하게 배수될 수 있도록 우수관을 Fig. 5와 같이 배치하였다. 우수관정의 단면적은 모두 2.0 m²이며, 우수관로의 직경은 종방향의 주관로인 P1∼P12와 P14~P26가 1.0 m, 바다에 접한 2개의 관로인 P13과 P27은 1.5 m, 나머지 지류 관로는 0.6 m이다. 그림에서 J는 우수관정(juction), P는 우수관(pipe)을 나타낸다. 각 알파벳 뒤에 부여된 숫자는 해당 관정이나 관로의 번호이다.

4. 수치모의 결과

본 연구에서 구축된 수치모형의 타당성을 확인하기 위해 Fig. 5에 보인 단순화된 지형을 이용하여 강우강도의 적용, 우수의 배수, 해수의 역류에 대한 수치모의를 수행하여 결과를 검토하였다. Fig. 6은 우수관망 시스템이 없는 경우에 대해 80 mm/hr의 강우강도를 적용하였을 때의 결과이다. 이때 강우는 지상영역에만 적용되었고, 지반, 해수면 및 방조제의 높이는 각각 2.0 m, 0.0 m, 3.2 m로 하였으며, 지상영역의 초기 침수 깊이는 0.0 m로서 건조상태이다. Fig. 6(a)는 강우강도 80mm/hr를 적용하였을 때 지상영역의 수위 변화를 나타낸 것이다. 그림을 보면 지상영역의 수위가 시간에 따라 일정하게 증가하는 것을 확인할 수 있다. 또한 적용된 누적강우량과 지상영역에 축적되는 우수량을 비교하기 위해 Fig. 6(b)와 같이 나타내었다. 그림을 통해 누적강우량과 지상영역에 고이는 누적우수량이 일치하는 것을 확인함으로써 강우강도가 정상적으로 적용되는 것을 알 수 있다.

Fig. 6.

Numerical simulation of rainfall without storm sewer system.

Fig. 7은 우수관망 시스템의 배수과정에 대한 검토 결과이다. 이때의 지반고, 해수면 및 방조제의 높이는 각각 2.0 m, 0.0 m, 3.2 m로 하였고, 강우는 없는 것으로 하였다. 지상영역의 초기 우수 깊이는 1.0 m로 하였으며, 각 관정의 초기 수위는 지상영역의 초기 수위와 동일한 3.0 m이다. 해수위는 매우 넓은 바다로 가정하여 0.0 m로 고정하였다. Fig. 7(a)는 우수가 바다로 배수될 때 시간에 따른 우수의 수위 변화를 나타낸 것이다. 그림을 보면 지상영역의 우수가 시간에 따라 감소하는 것을 확인할 수 있다. 관정이 방조제에 가까울수록 관정 위의 지상부 수심이 주변보다 크게 하강하는 것은 바다에 가까운 관정일수록 배수 유량이 크기 때문이다. 또한 우수관망 시스템을 통해 배수되는 누적 우수량과 지상영역에서 감소하는 누적 우수량을 비교하여 Fig. 7(b)와 같이 나타내었다. 관망을 통해 바다로 배수되는 누적 우수량은 Fig. 5에서 2개 관망의 바다 쪽 끝에 위치하는 관로인 P13과 P27의 통과 수량을 누적한 것이다. 그림을 통해 우수의 누적 배수량과 지상영역에서 감소되는 누적 우수량이 일치하는 것을 확인함으로써 우수 배수가 정상적으로 수행되는 것을 알 수 있다. 또한 600분 이후 배수된 누적 우수량이 일정해지며, 이는 지상영역의 우수가 모두 배수되어 더 이상 배수가 이루어지지 않음을 뜻한다.

Fig. 7.

Numerical simulation of drainage process through storm sewer system for the case of an open sea.

Fig. 8은 강우가 없을 때 해수면의 상승에 따른 해수의 역류현상을 수치모의 하였을 때의 결과이다. 이 경우 지반고와 방조제의 높이는 각각 2.0 m와 3.2 m로 배수과정 검토를 위한 수치모의 조건과 동일하게 하였고, 지상영역의 초기 침수 깊이는 0.0 m로서 초기에는 지상영역이 건조한 상태에서 시작하였다. 또한 지상영역에 위치하는 관정에서의 초기 수위는 평균해면고인 0.0 m로 하였으며, 해수면의 높이는 3.0로 고정하였다. Fig. 8(a)는 해수의 역류가 발생하였을 때 시간에 따라 지상영역에 쌓이는 해수의 수위 변화를 나타낸 것이다. 그림을 보면 해수가 우수관망을 통해 유입됨으로써 지상 영역에서의 수위가 증가하는 것을 확인할 수 있다. 또한 우수관망 시스템을 통해 역류한 해수의 양과 지상부에 쌓이는 해수의 증가량을 Fig. 8(b)와 같이 나타내었다. 그림을 통해 역류한 해수량과 지상영역에 쌓인 해수의 증가량이 일치하는 것을 확인함으로써 해수의 역류과정 또한 정상적으로 작동되는 것을 알 수 있다. 900분 이후 우수관망을 통해 역류하는 해수의 양이 일정해진 바, 이는 지상영역의 수위와 해수면의 높이가 동일해져 더 이상 해수가 역류하지 않아 이러한 결과가 나타난 것이라 판단된다.

Fig. 8.

Numerical simulation of back flow from the sea for the case of an open sea.

Fig. 9는 Fig. 7에서 검토한 우수의 배수과정에서 Fig. 5의 바다 경계를 따라 불투수벽을 설정하여 바다의 면적이 제한 될 때 지상영역의 우수배수가 바다의 해수면을 상승시키는 경우에 대한 수치모의 결과를 보여준다. 이때의 지반고, 초기 해수면, 방조제의 높이 및 지상영역의 초기 수위는 우수의 배수과정을 수치모의 하였을 때의 조건과 동일하게 하였다. Fig. 9(a)는 우수의 배수가 진행함에 따른 지상영역 및 해수면의 수위변화를 나타낸 것이다. 그림을 보면 지상영역의 초기 우수량이 시간에 따라 감소함과 동시에 해수면의 높이가 상승하는 것을 확인 할 수 있다. 또한 지상영역의 우수의 감소량, 우수관망시스템을 통해 배수되는 우수량 및 해수의 증가량을 비교하여 Fig. 9(b)와 같이 나타내었다. 그림을 통해 3가지 수량이 동일한 것으로 나타나 우수의 배수에 따른 해수면 상승 역시 정상적으로 수행되고 있음을 알 수 있다. 배수가 시작된 지 300분 이후에는 지상영역과 바다의 수위가 동일해져 더 이상의 배수가 발생하지 않는 것으로 나타났다.

Fig. 9.

Numerical simulation of drainage process through storm sewer system for the case of a closed sea.

이상의 수치모의를 통해 지상부가 수평인 지형에 대해 우수관망 수치모형의 기본적인 수량보존 성능을 검토하였다. 그러나 실제로는 지형이 경사진 곳이 대부분이므로 경사진 지형에 대한 강우 시의 우수관망 성능을 검토하였다. 이를 위해 Fig. 10에 보인 바와 같은 경사지를 대상으로 지상영역에 80 mm/hr의 강우강도를 적용하였고, 개방된 바다의 해수면의 높이는 각각 0.5 m, 1.5 m, 2.0 m, 2.5 m로 적용하여 수치모의를 수행하였다. Fig. 10(a)는 해수면의 높이에 따른 지상영역의 최종 수위 즉 최대 침수위를 비교한 것이다. 그림을 보면 대부분의 경우 지상영역의 침수심은 해수위가 높을수록 우수관로의 배수능력이 저하되어 증가하는 것으로 나타나 예상되는 자연현상이 잘 재현되었다. 특히 해수위가 2.5 m로 대부분의 지면부보다 높은 경우, 강우 초기에는 강우에 의한 침수와 바다로부터의 해수 역류 영향이 겹쳐 침수심이 빠르게 증가하다가, 지상부 수위가 해수위보다 높아진 다음부터 우수가 바다로 배수되어 배수량이 강우량과 동일해지는 시점에서 침수심은 더 이상 증가하지 않고 평형상태를 이루었다. 이때 지상부 수위는 해수위보다 0.35 m 정도 높은 상태를 유지하였다. 이러한 현상은 Fig. 10(b)에 제시한 바와 같이 누적강우량, 관로를 통한 배수-역류 수량과 지상영역에 쌓이는 물의 양을 비교함으로써 이해할 수 있다. 그림에서 강우에 의한 누적강우량이 우수의 배수-역류 수량과 지상영역에 쌓이는 물의 양을 합한 것과 일치하는 것으로 보아 강우와 우수관망 수치모형이 동시에 적용되었을 때에도 수치모의가 정상적으로 수행되는 것을 확인할 수 있다. 또한 그림에서 우수관망 시스템을 통해 배수-역류되는 물의 양을 보면 초기에는 음의 값을 나타내는 바, 이는 해수면의 높이가 지상부의 수위보다 높아 해수가 역류하는 현상이 잘 재현되고 있음을 알 수 있다. 우수관로 시스템이 작동한지 800분 이후 누적강우량과 누적배수량은 일정하게 증가하나, 지상부에 쌓인 물의양은 더 이상 변화가 없어 평형상태에 이른 다음에는 강우로 인한 우수가 모두 우수관로 시스템을 통해 바다로 지속적으로 배수되고 있음을 보여준다.

Fig. 10.

Numerical simulation of discharge process through storm sewer system for the case of rainfall over sloping ground.

5. 결 론

본 연구에서는 기존의 폭풍해일 수치모형에 강우에 의한 침수 증대 효과를 고려하기 위한 해일-강우-우수관망 결합 수치모형을 개발하였다. 폭풍해일의 연속방정식에 강우효과를 고려하였고, 동시에 우수관망 수치모형을 개발하여 폭풍해일 수치모형과 결합하였다. 수치모형을 검증하기 위해 이를 단순화된 지형을 이용하여 수치모의를 수행하였다. 먼저 우수관망 미작동시 강우에 의해 지상영역에 쌓이는 누적우수량이 총 강우량과 일치하는 지를 검토하였다. 다음으로는 강우가 없을 때 우수관망 설치에 따른 우수의 배수, 해수의 역류 및 우수의 배수에 따른 해수면의 변화에 대한 수치모의를 수행 하였다. 우수의 배수에서는 우수의 배수량과 지상영역에서의 우수의 감소량을 비교하였고, 해수의 역류에서는 해수면의 상승에 따른 해수의 역류량과 내륙지역 쌓이는 해수의 증가량을 비교하였다. 그 결과 각 조건에서 비교한 값이 일치하는 것을 확인하였다. 또한 우수의 배수로 인한 해수면의 높이가 변하는 경우에는 우수의 배수량, 내륙지역에 적용된 우수의 감소량 및 해수의 증가량을 비교한 결과 일치하는 것을 확인 하였다. 이와 같이 모든 조건에서 비교한 값이 일치한 것을 확인함으로써 적용된 강우강도 및 우수관망 수치모형의 기본 성능을 검증하였다. 최종적으로 강우강도와 우수관망 시스템을 동시에 적용하였을 때에도 정상적으로 수치모의를 수행하는지 추가적으로 검토하였다. 이때에는 누적강우량, 우수관망을 통해 배수 및 역류되는 물의 양과 내륙지역에 쌓이는 물의 양을 비교하였다. 그 결과 누적강우량, 우수관망을 통과하는 우수 및 해수의 양과 내륙지역에 쌓이는 물의 양을 더한 것과 일치하는 것을 확인함으로써 본 연구에서 구축된 우수관망 수치모형의 실효성을 검증할 수 있었다. 향후 이를 실제 지형에 적용하여 2차원 수치모의를 수행함으로써 기존의 폭풍해일에 의한 침수범람 결과보다 더욱 현실성 있는 수치모의 결과를 산정할 수 있을 것으로 기대된다.

Article information Continued