Watershed area

Analysis on inundation characteristics by compound external forces in coastal areas

연안 지역의 복합 외력에 의한 침수 특성 분석

Taeuk KangaDongkyun SunbSangho Leec*
강 태욱a선 동균b이 상호c*
aResearch Professor, Disaster Prevention Research Institute, Pukyong National University, Busan, KoreabResearcher, Disaster Prevention Research Institute, Pukyong National University, Busan, KoreacProfessor, Department of Civil Engineering, Pukyong National University, Busan, Korea
a부경대학교 방재연구소 전임연구교수b부경대학교 방재연구소 연구원c부경대학교 공과대학 토목공학과 교수*Corresponding Author

ABSTRACT

연안 지역은 강우, 조위, 월파 등 여러가지 외력에 의해 침수가 발생될 수 있다. 이에 이 연구에서는 연안 지역에서 발생될 수 있는 단일 및 복합 외력에 의한 지역별 침수 특성을 분석하였다. 연구에서 고려한 외력은 강우와 폭풍 해일에 의한 조위 및 월파이고, 분석 대상지역은 남해안 및 서해안의 4개 지역이다. 유역의 강우-유출 및 2차원 지표면 침수 분석에는 XP-SWMM이 사용되었고, 폭풍 해일에 의한 외력인 조위 및 월파량 산정에는 ADCSWAN (ADCIRC와 UnSWAN) 모형과 FLOW-3D 모형이 각각 활용되었다. 단일 외력을 이용한 분석 결과, 대부분의 연안 지역에서는 강우에 의한 침수 영향보다 폭풍 해일에 의한 침수 영향이 크게 나타났다. 복합 외력에 의한 침수 분석 결과는 대체로 단일 외력에 의한 침수 모의 결과를 중첩시켜 나타낸 결과와 유사하였다. 다만, 특정 지역에서는 복합 외력을 고려함에 따라 단일 외력만을 고려한 침수모의에서 나타나지 않았던 새로운 침수 영역이 발생하기도 하였다. 이러한 지역의 침수 피해 저감을 위해서는 복합 외력을 고려한 분석이 요구되는 것으로 판단되었다.키워드연안 지역 침수 분석 강우 폭풍 해일 복합 외력

The various external forces can cause inundation in coastal areas. This study is to analyze regional characteristics caused by single or compound external forces that can occur in coastal areas. Storm surge (tide level and wave overtopping) and rainfall were considered as the external forces in this study. The inundation analysis were applied to four coastal areas, located on the west and south coast in Republic of Korea. XP-SWMM was used to simulate rainfall-runoff phenomena and 2D ground surface inundation for watershed. A coupled model of ADCIRC and SWAN (ADCSWAN) was used to analyze tide level by storm surge and the FLOW-3D model was used to estimate wave overtopping. As a result of using a single external force, the inundation influence due to storm surge in most of the coastal areas was greater than rainfall. The results of using compound external forces were quite similar to those combined using one external force independently. However, a case of considering compound external forces sometimes created new inundation areas that didn’t appear when considering only a single external force. The analysis considering compound external forces was required to reduce inundation damage in these areas.KeywordsCoastal area Inundation analysis Rainfall Storm surge Compound external forces

MAIN

1. 서 론

우리나라는 반도에 위치하여 삼면이 바다로 둘러싸여 있는 지리적 특성을 가지고 있다. 이에 따라 해양 산업을 중심으로 부산, 인천, 울산 등 대규모의 광역도시가 발달하였을 뿐만 아니라, 창원, 포항, 군산, 목포, 여수 등의 중․소규모 도시들도 발달되어 있다. 또한, 최근에는 연안 지역이 바다를 전망으로 하는 입지 조건을 가지고 있어 개발 선호도가 높고, 이에 따라 부산시 해운대의 마린시티, 엘시티와 같은 주거 및 상업시설의 개발이 지속되고 있다(Kang et al., 2019b).

한편, 최근 기후변화에 따른 지구 온난화 현상으로 평균 해수면이 상승하고, 해수면 온도도 상승하면서 태풍 및 강우의 강도가 커지고 있어 전 세계적으로 자연 재해로 인한 피해가 증가하고 있다(Kim et al., 2016). 실제로 2020년에는 최장기간의 장마가 발생하여 부산, 울산은 물론, 전국에서 50명의 인명 피해와 3,489세대의 이재민이 발생하였다1). 특히, 연안 지역은 강우, 만조 시 해수면 상승, 폭풍 해일(storm surge)에 의한 월파(wave overtopping) 등 복합적인 외력(compound external forces)에 의해 침수될 수 있다(Lee et al., 2020). 일례로, 2016년 태풍 차바 시 부산시 해운대구의 마린시티는 강우와 폭풍 해일에 의한 월파가 발생함에 따라 대규모 침수를 유발하였다(Kang et al., 2019b). 또한, 2020년 7월 23일에 부산에서는 시간당 81.6 mm의 집중호우와 약최고고조위를 상회하는 만조가 동시에 발생하였고, 이로 인해 감조 하천인 동천의 수위가 크게 상승하여 하천이 범람하였다(KSCE, 2021).

연안 지역의 복합 외력을 고려한 침수 분석에 관한 사례로서, 우선 강우와 조위를 고려한 연구 사례는 다음과 같다. Han et al. (2014)은 XP-SWMM을 이용하여 창원시 배수 구역을 대상으로 침수 모의를 수행하였는데, 연안 도시의 침수 모의에는 조위의 영향을 반드시 고려해야 함을 제시하였다. Choi et al. (2018a)은 경남 사천시 선구동 일대에 대하여 초과 강우 및 해수면 상승 시나리오를 조합하여 침수 분석을 수행하였다. Choi et al. (2018b)은 XP-SWMM을 이용하여 여수시 연등천 및 여수시청 지역에 대하여 강우 시나리오와 해수위 상승 시나리오를 고려한 복합 원인에 의한 침수 모의를 수행하여 홍수예경보 기준표를 작성하였다. 한편, 강우, 조위, 월파를 고려한 연구 사례로서, Song et al. (2017)은 부산시 해운대구 수영만 일원에 대하여 XP-SWMM으로 월파량의 적용 유무에 따른 침수 면적을 비교하였다. Suh and Kim (2018)은 부산시 마린시티 지역을 대상으로 태풍 차바 때 EurOtop의 경험식을 ADSWAN에 적용하여 월파량을 반영하였다. Chen et al. (2017)은 TELEMAC-2D 및 SWMM을 기반으로 한 극한 강우, 월파 및 조위를 고려하여 중국 해안 원자력 발전소의 침수를 예측하고 분석하기 위한 결합 모델을 개발한 바 있다. 한편, Lee et al. (2020)은 수리‧수문학 분야와 해양공학 분야에서 사용되는 물리 모형의 기술적 연계를 통해 연안 지역의 침수 모의의 재현성을 높였다.

상기의 연구들은 공통적으로 연안 지역에 대하여 복합 외력을 고려했을 때 발생되는 침수 현상의 재현 또는 예측을 목적으로 수행되었다. 이 연구는 이와 차별하여 복합 외력을 고려하는 경우 나타날 수 있는 연안 지역의 침수 특성 분석을 목적으로 수행되었다. 이를 위해 단일 외력을 독립적으로 고려했을 때 발생되는 침수 양상과 동시에 고려하는 경우의 침수 현상을 비교, 분석하였다. 복합 외력에 의한 지역적 침수 특성 분석은 우리나라 남해안과 서해안에 위치한 4개 지역에 대하여 적용되었다.

1) 장연제, 47일째 이어진 긴 장마, 50명 인명피해… 9년만에 최대, 동아닷컴, 2020년 8월 9일 수정, 2021년 3월 4일 접속, https://www.donga.com/news/article/all/20200809/102369692/2

2. 연구 방법

2.1 연안 지역의 침수 영향 인자

연안 지역의 침수는 크게 세 가지의 메카니즘으로 발생될 수 있다. 우선, 연안 지역은 바다와 인접하고 있기 때문에 그 영향을 직접적으로 받는다. Kim (2018)에 의하면, 연안 지역의 침수는 폭풍 해일에 의해 상승한 조위와 월파로 인해 발생될 수 있다(Table 1). 특히, 경상남도의 창원과 통영, 인천광역시의 소래포구 어시장 등 남해안 및 서해안 지역의 일부는 백중사리, 슈퍼문(super moon) 등 만조 시 조위의 상승으로 인한 침수가 발생하는 지역이 존재한다(Kang et al., 2019a). 두 번째는 강우에 의한 내수 침수 발생이다. ME (2011)에서는 도시 지역의 우수 관거를 10 ~ 30년 빈도로 계획하도록 지정하고 있고, 펌프 시설은 30 ~ 50년 빈도의 홍수를 배수시킬 수 있도록 정하고 있다. 하지만 최근에는 기후변화의 영향으로 도시 지역 배수시설의 설계 빈도를 초과하는 강우가 빈번하게 나타나고 있다. 실제로 2016년의 태풍 차바 시 울산 기상관측소에 관측된 시간 최대 강우량은 106.0 mm로서, 이는 300년 빈도 이상의 강우량에 해당하였다(Kang et al., 2019a). 따라서 배수시설의 설계 빈도 이상의 강우는 연안 도시 지역의 침수를 유발할 수 있다. 세 번째, 하천이 인접한 연안 도시에서는 하천의 범람으로 인해 침수가 발생할 수 있다. 하천의 경우, 기본계획이 수립되기는 하지만, 설계 빈도를 상회하는 강우의 발생, 제방, 수문 등 홍수 방어시설의 기능 저하, 예산 등의 문제로 하천기본계획 이행의 지연 등에 의해 범람할 가능성이 존재한다.

Table 1.

Type of natural hazard damage in coastal areas (Kim, 2018)

ItemRisk factor
Facilities damage∙ Breaking of coastal facilities by wave
– Breakwater, revetment, lighters wharf etc.
∙ Local scouring at the toe of the structures by wave
∙ Road collapse by wave overtopping
Inundation damage∙ Inundation damage by wave overtopping
∙ Inundation of coastal lowlands by storm surge
Erosion damage∙ Backshore erosion due to high swell waves
∙ Shoreline changes caused by construction of coastal erosion control structure
∙ Sediment transport due to the construction of artificial structures

상기의 내용을 종합하면, 연안 지역은 조위 및 월파에 의한 침수, 강우에 의한 내수 침수, 하천 범람에 의한 침수로 구분될 수 있다. 이 연구에서는 폭풍 해일에 의한 조위 상승 및 월파와 강우를 연안 지역의 침수 유발 외력으로 고려하였다. 하천 범람의 경우, 상대적으로 사례가 희소하여 제외하였다.

2.2 복합 외력을 고려한 침수 모의 방법

이 연구에서는 조위 및 월파와 강우를 연안 지역의 침수 발생에 관한 외력 조건으로 고려하였다. 따라서 해당 외력 조건을 고려하여 침수 분석을 수행할 수 있어야 한다. 이와 관련하여 Lee et al. (2020)은 Fig. 1과 같이 수리‧수문 및 해양공학 분야에서 사용되는 물리 기반 모형의 연계를 통해 조위, 월파, 강우를 고려한 침수 분석 방법을 제시하였고, 이 연구에서는 해당 방법을 이용하였다.

/media/sites/kwra/2021-054-07/N0200540702/images/kwra_54_07_02_F1.jpg
Fig. 1.

Connection among the models for inundation analysis in coastal areas (Lee et al., 2020)

우선, 태풍에 의해 발생되는 폭풍 해일의 영향을 분석하기 위해서는 태풍에 의해 발생되는 기압 강하, 해상풍, 진행 속도 등을 고려하여 해수면의 변화 양상 및 조석-해일-파랑을 충분히 재현 가능해야 한다. 이 연구에서는 국내․외에서 검증 및 공인된 폭풍 해일 모형인 ADCIRC 모형과 파랑 모형인 UnSWAN이 결합된 ADCSWAN (coupled model of ADCIRC and UnSWAN)을 이용하였다. 정수압 가정의 ADCSWAN은 월파량 산정에 단순 경험식을 적용하는 단점이 있지만 넓은 영역을 모의할 수 있고, FLOW-3D는 해안선의 경계를 고해상도로 재현이 가능하다. 이에 연구에서는 먼 바다 영역에 대해서는 ADCSWAN을 이용하여 분석하였고, 연안 주변의 바다 영역과 월파량 산정에 대해서는 FLOW-3D 모형을 이용하였다. 한편, 연안 지역의 침수 모의를 위해서는 유역에서 발생하는 강우-유출 현상과 우수 관거 등의 배수 체계에 대한 분석이 가능해야 한다. 또한, 배수 체계로부터 범람한 물이 지표면을 따라 흘러가는 현상을 해석할 수 있어야 하고, 바다의 조위 및 월파량을 경계조건으로 반영할 수 있어야 한다. 이 연구에서는 이러한 현상을 모의할 수 있고, 도시 침수 모의에 활용도가 높은 XP-SWMM을 이용하였다.

2.3 침수 분석 대상지역

연구의 대상지역은 조위 및 월파에 의한 침수와 강우에 의한 내수 침수의 영향이 복합적으로 발생할 수 있는 남해안과 서해안에 위치한 4개 지역이다. Table 2는 침수 분석 대상지역을 정리하여 나타낸 표이고, Fig. 2는 각 지역의 유역 경계를 나타낸 그림이다.

Table 2.

Target region for inundation analysis

ClassificationAdministrative districtTarget regionArea
(km2)
Main cause of inundationPump
facility
Number of
major outfall
The south
coast
Haundae-gu, BusanMarine City area0.53Wave overtopping9
Haundae-gu, BusanCentum City area4.76Poor interior drainage at high tide level12
The west
coast
GunsanJungang-dong area0.79Poor interior drainage at high tide level23
BoryeongOcheon Port area0.41High tide level5

/media/sites/kwra/2021-054-07/N0200540702/images/kwra_54_07_02_F2.jpg
Fig. 2.

Watershed area

남해안의 분석 대상지역 중 부산시 해운대구의 마린시티는 바다 조망을 중심으로 조성된 주거지 및 상업시설 중심의 개발지역이다. 마린시티는 2016년 태풍 차바 및 2018년 태풍 콩레이 등 태풍 내습 시 월파에 의한 해수 월류로 인해 도로 및 상가 일부가 침수를 겪은 지역이다. 부산시 해운대구의 센텀시티는 과거 수영만 매립지였던 곳에 조성된 주거지 및 상업시설 중심의 신도시 지역이다. 센텀시티 유역의 북쪽은 해발고도 El. 634 m의 장산이 위치하는 등 산지 특성도 가지고 있어 상대적으로 유역 면적이 넓고, 배수시설의 규모도 크고 복잡하다. 하지만 수영강 하구의 저지대 지역에 위치함에 따라 강우 시 내수 배제가 불량하고, 특히 만조 시 침수가 잦은 지역이다.

서해안 분석 대상지역 중 전라북도 군산시의 중앙동 일원은 군산시 내항 내측에 조성된 구도시로서, 금강 및 경포천 하구에 위치하는 저지대이다. 이에 따라 군산시 풍수해저감종합계획에서는 해당 지역을 3개의 영역으로 구분하여 내수재해 위험지구(영동지구, 중동지구, 경암지구)로 지정하였고, 이 연구에서는 해당 지역을 모두 고려하였다. 한편, 군산시 중앙동 일원은 특히, 만조 시 내수 배제가 매우 불량하여 2개의 펌프시설이 운영되고 있다. 충청남도 보령시의 오천면에 위치한 오천항은 배후의 산지를 포함한 소규모 유역에 위치한다. 서해안의 특성에 따라 조석 간만의 차가 크고, 특히 태풍 내습 시 폭풍 해일에 의한 침수가 잦은 지역이다. 산지의 강우-유출수는 복개된 2개의 수로를 통해 바다로 배제되고, 상가들이 위치한 연안 주변 지역에는 강우-유출수 배제를 위한 3개의 배수 체계가 구성되어 있다.

3. 연구 결과

3.1 침수 모의 모형 구축

XP-SWMM을 이용하여 분석 대상지역별 침수 모의 모형을 구축하였다. 적절한 침수 분석 수행을 위해 지역별 수치지형도, 도시 공간 정보 시스템(urban information system, UIS), 하수 관망도 등의 수치 자료와 현장 조사를 통해 유역의 배수 체계를 구성하였다. 그리고 2차원 침수 분석을 위해 무인 드론 및 육상 라이다(LiDAR) 측량을 수행하여 평면해상도가 1 m 이하인 고해상도 수치지형모형(digital terrain model, DTM)을 구성하였고, 침수 모의 격자를 생성하였다.

Fig. 3은 XP-SWMM의 상세 구축 사례로서 부산시 마린시티 배수 유역에 대한 소유역 및 관거 분할 등을 통해 구성한 배수 체계와 고해상도 측량 결과를 이용하여 구성한 수치표면모형(digital surface model, DSM)을 나타낸다. Fig. 4는 각 대상지역에 대해 XP-SWMM을 이용하여 구축한 침수 모의 모형을 나타낸다. 침수 분석을 위해서는 침수 모의 영역에 대한 설정이 필요한데, 다수의 사전 모의를 통해 유역 내에서 침수가 발생되는 지역을 검토하여 결정하였다.

/media/sites/kwra/2021-054-07/N0200540702/images/kwra_54_07_02_F3.jpg
Fig. 3.

Analysis of watershed drainage system and high-resolution survey for Marine City

/media/sites/kwra/2021-054-07/N0200540702/images/kwra_54_07_02_F4.jpg
Fig. 4.

Simulation model for inundation analysis by target region using XP-SWMM

한편, 이 연구에서는 월파량 및 조위의 산정 과정과 침수 모의 모형의 보정에 관한 내용 등은 다루지 않았다. 관련된 내용은 선행 연구인 Kang et al. (2019b)와 Lee et al. (2020)을 참조할 수 있다.

3.2 침수 모의 설정

3.2.1 분석 방법

복합 외력에 의한 침수 영향을 검토하기 위해서는 외력 조건에 대한 빈도와 지속기간의 설정이 필요하다. 이 연구에서는 재해 현상이 충분히 나타날 수 있도록 강우와 조위 및 월파의 빈도를 모두 100년으로 설정하였다. 이때, 조위와 월파량의 산정에는 만조(약최고고조위) 시, 100년 빈도에 해당하는 태풍 내습에 따른 폭풍 해일의 발생 조건을 고려하였다.

지역별 강우 발생 특성과 유역 특성을 고려하기 위해 MOIS (2017)의 방재성능목표 기준에 따라 임계 지속기간을 결정하여 대상지역별 강우의 지속기간으로 설정하였다. 이때, 강우의 시간 분포는 MLTM (2011)의 Huff 3분위를 이용하였다. 그리고 조위와 월파의 경우, 일반적인 폭풍 해일의 지속기간을 고려하여 5시간으로 결정하였다. 한편, 침수 모의를 위한 계산 시간 간격, 2차원 모의 격자 등의 입력자료는 분석 대상지역의 유역 규모와 침수 분석 대상 영역을 고려하여 결정하였다. 참고로 침수 분석에 사용된 수치지형모형은 1 m 급의 고해상도로 구성되었지만, 2차원 침수 모의 격자의 크기는 지역별로 3 ~ 4 m이다. 이는 연구에서 사용된 XP-SWMM의 격자 수(100,000개) 제약에 따른 설정이나, Sun (2021)은 민감도 분석을 통해 2차원 침수 분석을 위한 적정 격자 크기를 3 ~ 4.5 m로 제시한 바 있다.

Table 3은 이 연구에서 설정한 침수 모의 조건과 분석 방법을 정리하여 나타낸 표이다.

Table 3.

Simulation condition and method

ClassificationTarget regionSimulation conditionSimulation method
RainfallStorm surgeSimulation time interval2D
grid size
Return
period
DurationTemporal
distribution
Return
period
DurationWatershed
routing
Channel
routing
2D
inundation
The south coastMarine City area100 yr1 hr3rd quartile
of Huff’s
method
1005 hr5 min10 sec1 sec3 m
Centum City area1 hr1005 min10 sec1 sec4 m
The west coastJungang-dong area2 hr1005 min10 sec1 sec3.5 m
Ocheon Port area1 hr1001 min10 sec1 sec3 m

3.2.2 복합 재해의 동시 고려

이 연구의 대상지역들은 모두 소규모의 해안가 도시지역이고, 이러한 지역에 대한 강우의 임계지속기간은 1시간 ~ 2시간이나, 이 연구에서 분석한 폭풍 해일의 지속기간은 5시간으로 강우의 지속기간과 폭풍 해일의 지속기간이 상이하다. 이에 이 연구에서는 서로 다른 지속기간을 가진 강우와 폭풍 해일 또는 조위를 고려하기 위해 강우의 중심과 폭풍 해일의 중심이 동일한 시간에 위치하도록 설정하였다(Fig. 5).

XP-SWMM은 폭풍 해일이 지속되는 5시간 전체를 모의하도록 설정하였고, 폭풍 해일이 가장 큰 시점에 강우의 중심이 위치하도록 강우 발생 시기를 결정하였다. 다만, 부산 마린시티의 경우, 폭풍 해일에 의한 피해가 주로 월파에 의해 발생되므로 강우의 중심과 월파의 중심을 일치시켰고(Fig. 5(a)), 상대적으로 조위의 영향이 큰 3개 지역은 강우의 중심과 조위의 중심을 맞추었다. Fig. 5(b)는 군산시 중앙동 지역의 복합 외력에 의한 침수 분석에 사용된 강우와 조위의 조합이다.

한편, 100년 빈도의 확률강우량만을 고려한 침수 분석에서는 유역 유출부의 경계조건으로 우수 관거의 설계 조건을 고려하여 약최고고조위가 일정하게 유지되도록 설정하였다.

/media/sites/kwra/2021-054-07/N0200540702/images/kwra_54_07_02_F5.jpg
Fig. 5.

Consideration of external force conditions with different durations

3.2.3 XP-SWMM의 월파량 고려

XP-SWMM에 ADCSWAN 및 FLOW-3D 모형에 의해 산정된 월파량을 입력하기 위해 해안가 지역에 절점을 생성하여 월파 현상을 구현하였다. XP-SWMM에서 월파량을 입력하기 위한 절점의 위치는 FLOW-3D 모형에서 월파량을 산정한 격자의 중심 위치이다.

Fig. 6(a)는 마린시티 지역에 대한 월파량 입력 지점을 나타낸 것으로서, 유역 경계 주변에 동일 간격으로 원으로 표시한 지점들이 해당된다. Fig. 6(b)는 XP-SWMM에 월파량 입력 지점들을 반영하고, 하나의 절점에 월파량 시계열을 입력한 화면을 나타낸다.

/media/sites/kwra/2021-054-07/N0200540702/images/kwra_54_07_02_F6.jpg
Fig. 6.

Considering wave overtopping on XP-SWMM

3.3 침수 모의 결과

3.3.1 단일 외력에 의한 침수 모의 결과

Fig. 7은 단일 외력을 고려한 지역별 침수 모의 결과이다. 즉, Fig. 7의 왼쪽 그림들은 지역별로 100년 빈도 강우에 의한 침수 모의 결과를 나타내고, Fig. 7의 오른쪽 그림들은 만조 시 100년 빈도 폭풍 해일에 의한 침수 모의 결과이다. 대체로 강우에 의한 침수 영역은 유역 중․상류 지역의 유역 전반에 걸쳐 발생하였고, 폭풍 해일에 의한 침수 영역은 해안가 전면부에 위치하는 것을 볼 수 있다. 이는 폭풍 해일에 의한 조위 상승과 월파의 영향이 상류로 갈수록 감소하기 때문이다.

한편, 4개 지역 모두에서 공통적으로 강우에 비해 폭풍 해일에 의한 침수 영향이 상대적으로 크게 분석되었다. 이러한 결과는 연안 지역의 경우, 폭풍 해일에 대비한 침수 피해 저감 노력이 보다 중요함을 의미한다.

/media/sites/kwra/2021-054-07/N0200540702/images/kwra_54_07_02_F7.jpg
Fig. 7.

Simulation results by single external force (left: rainfall, right: storm surge)

3.3.2 복합 외력에 의한 침수 모의 결과

Fig. 8은 복합 외력을 고려한 지역별 침수 모의 결과이다. 즉, 강우 및 폭풍 해일을 동시에 고려함에 따라 발생된 침수 영역을 나타낸다. 복합 외력을 고려하는 경우, 단일 외력만을 고려한 분석 결과(Fig. 7)보다 침수 영역은 넓어졌고, 침수심은 깊어졌다.

복합 외력에 의한 침수 분석 결과는 대체로 단일 외력에 의한 침수 모의 결과를 중첩시켜 나타낸 결과와 유사하였고, 이는 일반적으로 예상할 수 있는 결과이다. 주목할만한 결과는 군산시 중앙동의 침수 분석에서 나타났다. 즉, 군산시 중앙동의 경우, 단일 외력만을 고려한 침수 모의 결과에서 나타나지 않았던 새로운 침수 영역이 발생하였다(Fig. 8(c)). 이와 관련된 상세 내용은 3.4절의 고찰에서 기술하였다.

/media/sites/kwra/2021-054-07/N0200540702/images/kwra_54_07_02_F8.jpg
Fig. 8.

Simulation results by compound external forces

3.4 결과 고찰

외력 조건별 침수의 영향을 정량적으로 비교하기 위해 침수 면적을 이용하였다. 이 연구에서는 강우만에 의해 유발된 침수 면적을 기준(기준값: 1)으로 하고, 폭풍 해일(조위+월파량)에 의한 침수 면적과 복합 외력에 의한 침수 면적의 상대적 비율로 분석하였다(Table 4).

Table 4.

Impact evaluation for inundation area by external force

ConditionMarine City, BusanCentum City, BusanJungang-dong area,
Gunsan
Ocheon Port area,
Boryeong
Inundation area
(km2)
RateInundation area
(km2)
RateInundation area
(km2)
RateInundation area
(km2)
Rate
Single
external force
Rainfall (①)0.01641.00.07591.00.04571.00.01751.0
Storm surge (②)0.03632.210.06850.900.14633.200.04122.35
Compound
external forces
Combination
(①+②)
0.05243.190.15051.980.26325.760.04732.70

분석 결과, 부산 센텀시티를 제외한 3개 지역은 모두 폭풍 해일에 의한 침수 면적이 강우에 의한 침수 면적에 비해 2.2 ~ 3.2배 넓은 것으로 분석되었다. 한편, 복합 외력에 의한 침수 면적은 마린시티와 센텀시티의 경우, 각각의 외력에 의한 침수 면적의 합과 유사하게 나타났다. 이는 각각의 외력에 의한 침수 영역이 상이하여 거의 중복되지 않음을 의미한다. 반면에, 오천항에서는 각각의 외력에 의한 침수 면적의 합이 복합 외력에 의한 면적보다 크게 나타났다. 이는 오천항의 경우, 유역면적이 작고 배수 체계가 비교적 단순하여 강우와 폭풍 해일에 의한 침수 영역이 중복되기 때문인 것으로 분석되었다(Fig. 7(d)).

군산시 중앙동 일대의 경우, 복합 외력에 의한 침수 면적이 각각의 독립적인 외력 조건에 의한 침수 면적의 합에 비해 37.1% 크게 나타났다. 이러한 현상의 원인을 분석하기 위해 복합 외력 조건에서만 나타난 우수 관거(Fig. 8(c)의 A 구간)에 대하여 종단을 검토하였다(Fig. 9). Fig. 9(a)는 강우만에 의해 분석된 우수 관거 내 흐름 종단을 나타내고, Fig. 9(b)는 폭풍 해일만에 의한 우수 관거의 종단이다. 그림을 통해 각각의 독립적인 외력 조건 하에서는 해당 구간에서 침수가 발생되지 않은 것을 볼 수 있다. 다만, 강우만을 고려하더라도 우수 관거는 만관이 된 상태를 확인할 수 있다(Fig. 9(a)). 반면에, 만관 상태에서 폭풍 해일이 함께 고려됨에 따라 해수 범람과 조위 상승에 의해 우수 배제가 불량하게 되었고, 이로 인해 침수가 유발된 것으로 분석되었다(Fig. 9(c)). 따라서 이러한 지역은 복합 외력에 대한 취약지구로 판단할 수 있고, 단일 외력의 고려만으로는 침수를 예상하기 어려운 지역임을 알 수 있다.

/media/sites/kwra/2021-054-07/N0200540702/images/kwra_54_07_02_F9.jpg
Fig. 9.

A part of drainage profiles by external force in Jungang-dong area, Gunsan

4. 결 론

이 연구에서는 외력 조건에 따른 연안 지역의 침수 특성을 분석하였다. 연구에서 고려된 외력 조건은 두 가지로서 강우와 폭풍 해일(조위와 월파)이다. 분석 대상 연안 지역으로는 남해안에 위치하는 2개 지역(부산시 해운대구의 마린시티와 센텀시티)과 서해안의 2개 지역(군산시 중앙동 일원 및 보령시 오천항)이 선정되었다.

복합 외력을 고려한 연안 지역의 침수 모의를 위해서는 유역의 강우-유출 현상과 바다의 조위 및 월파량을 경계조건으로 반영할 수 있는 침수 모의 모형이 요구되는데, 이 연구에서는 XP-SWMM을 이용하였다. 한편, 조위 및 월파량 산정에는 ADCSWAN (ADCIRC와 UnSWAN) 및 FLOW-3D 모형이 이용되었다.

연안 지역별 침수 모의는 100년 빈도의 강우와 폭풍 해일을 독립적으로 고려한 경우와 복합적으로 고려한 경우를 구분하여 수행되었다. 우선, 외력을 독립적으로 고려한 결과, 대체로 폭풍 해일만 고려한 경우가 강우만 고려한 경우에 비해 침수 영향이 크게 나타났다. 따라서 연안 지역의 경우, 폭풍 해일에 의한 침수 피해 방지 계획이 상대적으로 중요한 것으로 분석되었다. 두 번째, 복합 외력에 의한 침수 분석 결과는 대체로 단일 외력에 의한 침수 모의 결과를 중첩시켜 나타낸 결과와 유사하였다. 다만, 특정 지역에서는 복합 외력을 고려함에 따라 단일 외력만을 고려한 침수 모의에서 나타나지 않았던 새로운 침수 영역이 발생하기도 하였다. 이러한 결과는 독립적인 외력 조건에서는 우수 관거가 만관 또는 그 이하의 상태가 되지만, 두 가지의 외력이 동시에 고려됨에 따라 우수 관거의 통수능 한계를 초과하여 나타났다. 이러한 지역은 복합 외력에 대한 취약지구로 판단되었고, 해당 지역의 적절한 침수 방지 대책 수립을 위해서는 복합적인 외력 조건이 고려되어야 함을 시사하였다.

현행, 자연재해저감종합계획에서는 침수와 관련된 재해 원인 지역을 내수재해, 해안재해, 하천재해 등으로 구분하고 있다. 하지만 이 연구에서 검토된 바와 같이, 연안 지역의 침수 원인은 복합적으로 나타날 뿐만 아니라, 복합 외력을 고려함에 따라 추가적으로 나타날 수 있는 침수 위험 지역도 존재한다. 따라서 기존의 획일적인 재해 원인의 구분보다는 지역의 특성에 맞는 복합적인 재해 원인을 검토할 필요가 있음을 제안한다.

Acknowledgements

본 논문은 행정안전부 극한 재난대응 기반기술 개발사업의 일환인 “해안가 복합재난 위험지역 피해저감 기술개발(연구과제번호: 2018-MOIS31-008)”의 지원으로 수행되었습니다.

References

1
Chen, X., Ji, P., Wu, Y., and Zhao, L. (2017). “Coupling simulation of overland flooding and underground network drainage in a coastal nuclear power plant.” Nuclear Engineering and Design, Vol. 325, pp. 129-134. 10.1016/j.nucengdes.2017.09.028
2
Choi, G., Song, Y., and Lee, J. (2018a). “Analysis of flood occurrence type according to complex characteristics of coastal cities.” 2018 Conference of the Korean Society of Hazard Mitigation, KOSHAM, p. 180.
3
Choi, J., Park, K., Choi, S., and Jun, H. (2018b). “A forecasting and alarm system for reducing damage from inland inundation in coastal urban areas: A case study of Yeosu City.” Journal of Korean Society of Hazard Mitigation, Vol. 18, No. 7, pp. 475-484. 10.9798/KOSHAM.2018.18.7.475
4
Han, H., Kim, Y., Kang, N., and, Kim, H.S. (2014). “Inundation analysis of a coastal urban area considering tide level.” 2014 Conference of Korean Society of Civil Engineers, KSCE, pp. 1507-1508.
5
Kang, T., Lee, S., and Sun, D. (2019a). “A technical review for reducing inundation damage to high-rise and underground-linked complex buildings in Coastal Areas (1): Proposal for analytical method.” Journal of Korean Society of Hazard Mitigation, Vol. 19, No. 5, pp. 35-43. 10.9798/KOSHAM.2019.19.5.35
6
Kang, T., Lee, S., Choi, H., and Yoon, S. (2019b). “A technical review for reducing inundation damage to high-rise and underground-linked complex buildings in coastal areas (2): Case analysis for application.” Journal of Korean Society of Hazard Mitigation, Vol. 19, No. 5, pp. 45-53. 10.9798/KOSHAM.2019.19.5.45
7
Kim, J.O., Kim, J.Y., and Lee, W.H. (2016). “Analysis on complex disaster information contents for building disaster map of coastal cities.” Journal of the Korean Association of Geographic Information Studies, Vol. 19, No. 3, pp. 43-60. 10.11108/kagis.2016.19.3.043
8
Kim, P.J. (2018). Improvement measures on the risk area designation of coastal disaster in consideration of natural hazards. Ph.D. dissertation, Chonnam National University.
9
Korean Society of Civil Engineers (KSCE) (2021). A report on the cause analysis and countermeasures establishment for Dongcheon flooding and lowland inundation. Busan/Ulsan, Gyungnam branch.
10
Lee, S., Kang, T., Sun, D., and Park, J.J. (2020). “Enhancing an analysis method of compound flooding in coastal areas by linking flow simulation models of coasts and watershed.” Sustainability, Vol. 12, No. 16, 6572. 10.3390/su12166572
11
Ministry of Environment (ME) (2011). Standard for sewerage facilities. Korea Water and Wastewater Works Association.
12
Ministry of Land, Transport and Maritime Affairs (MLTM) (2011). Improvement and complementary research for probability rainfall.
13
Ministry of the Interior and Safety (MOIS) (2017). Criteria for establishment and operation of disaster prevention performance target by region: Considering future climate change impacts.
14
Song, Y., Joo, J., Lee, J., and Park, M. (2017). “A study on estimation of inundation area in coastal urban area applying wave overtopping.” Journal of Korean Society of Hazard Mitigation, Vol. 17, No. 2, pp. 501-510. 10.9798/KOSHAM.2017.17.2.501
15
Suh, S.W., and Kim, H.J. (2018). “Simulation of wave overtopping and inundation over a dike caused by Typhoon Chaba at Marine City, Busan, Korea.” Journal of Coastal Research, Vol. 85, pp. 711-715.
16
Sun, D. (2021). Sensitivity analysis of XP-SWMM for inundation analysis in coastal area. M.Sc. Thesis, Pukyong National University.

입자 직경과 질량이 서로 다른 두 가지 변화로 인한 시뮬레이션 결과

유체 역학을 사용한 미세 유체 입자 분류

Microfluidics Particle Sorting Using Hydrodynamics

미세 유체 입자 분류는 진단, 화학적 및 생물학적 분석, 식품 및 화학 처리, 환경 평가에 적용됩니다. 미세 유체 분류 플랫폼 사용의 주요 이점은 더 작은 샘플 볼륨이 필요하다는 것입니다. 그 결과 비용과 시간이 감소합니다. 진단에서는 환자의 침습성을 감소시킵니다. 또한, 이러한 소형 플랫폼은 대규모 병렬화가 가능하므로 더 적은 시간에 더 많은 수의 입자 분류가 가능합니다.

수동 또는 능동 기술을 사용하여 정렬을 수행 할 수 있습니다. 패시브 기술은 외부 필드를 필요로하지 않으며 입자, 유동장 및 채널 구조 간의 상호 작용에만 기반합니다. 반면에 능동적 기술은 자기 또는 전기와 같은 외부 필드를 사용합니다. 아래 애니메이션은 미세 유체 플랫폼의 유체 역학을 기반으로하는 수동 분류 기술을 사용하여 세 가지 입자 종의 입자 분류를 보여줍니다.

최첨단 FlowSight에서 후 처리 된 FLOW-3D  입자 물리학 모델을 사용한 미세 유체 입자 분류 시뮬레이션

이 블로그에서는 위에 표시된 정렬 기술의이면에있는 물리학과 이러한 시뮬레이션을 수행 할 때 FLOW-3D 를 사용하는 방법에 대해 설명합니다 .

유체 역학적 분류 기술 뒤에있는 물리학

이 기술은 낮은 레이놀즈 수 영역에서 입자가 질량과 직경을 기준으로 유동장에서 특정 유선을 따른다는 원칙에 따라 작동합니다. 질량이 일정하다는 점을 감안할 때 직경이 작은 입자는 항력이 적고 직경이 큰 입자는 항력이 더 커집니다. 이로 인해 더 큰 입자가 주변의 흐름과 함께 쉽게 제거됩니다. 더 작은 입자와 그 궤적은 유체 역학적 힘의 영향을 덜받습니다.

직경이 모든 입자 종에 대해 동일하지만 밀도와 질량이 다를 때 입자 분류 중에 다른 동작을 볼 수 있습니다. 같은 직경의 입자에 대한 항력은 동일하지만 더 무거운 입자는 더 큰 관성력의 영향을 받기 때문에 감속하기가 어렵습니다. 반대로, 더 가벼운 입자는 감속하기가 더 쉽습니다. 따라서 결과적으로 가벼운 입자는 주변의 흐름과 함께 쉽게 사라지고 무거운 입자는 그 경로를 유지합니다.

입자 분류 시뮬레이션을위한 FLOW-3D

FLOW-3D 의 입자 모델을 사용하면 입자 정렬 시뮬레이션을 매우 쉽게 수행 할 수 있습니다. 모델에는 마커, 질량, 유체, 가스 또는 보이드 입자와 같은 다양한 입자 클래스를 설정할 수있는 옵션이 있습니다. 이 시뮬레이션을 위해 질량 입자가 사용되었습니다. 특정 유형의 입자 등급은 직경과 밀도에 따라 다른 종을 가질 수 있습니다. 예를 들어, 위의 애니메이션에는 질량 입자 클래스에 대해 세 가지 종이 있습니다.

질량 입자의 역학은 확산 계수, 항력 계수, 난류 슈미트 수 및 복원 계수와 같은 속성으로 제어 할 수 있습니다. 질량 입자는 열적 및 전기적 특성을 지정할 수도 있습니다. 사용자가 입자에 동시에 작용하는 여러 힘을 연구하려는 경우 이러한 속성을 완전히 활용할 수 있습니다.

결과

아래 이미지는 질량과 반경을 변경하는 두 가지 시뮬레이션 사례의 결과를 보여줍니다.

입자 직경과 질량이 서로 다른 두 가지 변화로 인한 시뮬레이션 결과
입자 직경과 질량이 서로 다른 두 가지 변화로 인한 시뮬레이션 결과

직경이 작은 입자 (왼쪽 창에 파란색) 또는 질량 (오른쪽 창에 녹색)이 적은 입자는 수축의 상단으로 이동하고 발산 지점에서 위쪽으로 발산하는 유선을 따릅니다. 직경이 더 크거나 (왼쪽 창에서 녹색) 질량이 더 큰 (오른쪽 창에서 파란색) 입자는 수축의 아래쪽으로 이동합니다. 수축을 빠져 나가면 이러한 입자는 아래쪽으로 갈라지는 유선을 따릅니다.

채널 형상, 입자 매개 변수 및 유동 특성을 기반으로하는 미세 유체 입자 분류 장치의 미세 입자 분류에 대한 정확한 수치 분석은 이러한 미세 장치의 더 나은 설계에 사용될 수 있습니다.

함께 FLOW-3D 의 강력한 입자 모델, 미세 유체 입자 정렬 시뮬레이션을 설정하기 쉽습니다. 블로그 시작 부분의 애니메이션은 서로 다른 입자 종의 깔끔한 분류와 서로 다른 콘센트에있는 각 종의 수집을 보여줍니다.

Moving Boundaries: An Eulerian Approach

Moving Boundaries: An Eulerian Approach

많은 문제에서, 유체 및 고체 영역의 내부 경계가 그 안에서 이동할 수 있도록하면서 공간에 고정 된 그리드를 유지하는 것이 유리합니다. 이는 리 메싱의 필요성을 피할 수 있으므로 이러한 경계의 형태에 급격한 변화가 발생할 때마다 적절합니다. 메시 생성도 크게 단순화되었습니다.

고정 그리드 내에서 유체 인터페이스, 침전물, 응고 된 유체 및 탄성 재료의 경계 이동을 모델링하기위한 다양한 접근 방식이 표시됩니다. 유체 경계의 이동은 VOF (Volume-of-Fluid) 방법의 변형으로 수행되며, 각 계산 셀에서 유체의 양을 나타내는 양이 고정 메시를 통해 조정됩니다.

퇴적물의 침식 및 퇴적은 퇴적물 수색 모델을 사용하여 계산됩니다. 국부적 인 침식 속도는 패킹 된 퇴적물 / 유체 경계면에 존재하는 국부적 인 전단 응력을 기반으로하며, 증착은 Stokes 유동 근사치로 예측됩니다.

Emptying of gravure cell (same cell dimensions as filling case); a
three-dimensional perspective is shown. The transfer roll surface
(block at top) is moving away from the gravure roll at 0.5m/s. The
static contact of the fluid with all surfaces is 30°. The elapsed time
is 150

충진 층 경계면은 퇴적물 농도와 퇴적물의 포장 분율에 따라 달라집니다. 용융 금속은 온도가 빙점 아래로 떨어지면 굳을 수 있습니다. 응고 된 “유체”는 동결 및 용융을 유발하는 열유속의 양으로부터 결정된대로 표면이 증가하거나 수축하는 고체처럼 처리됩니다.

탄성 응력은 응고 된 재료 / 공기 인터페이스를 예측하는 VOF 방법을 사용하여 동일한 고정 그리드 내의 운동량 균형에 탄성 응력 계산을 추가하여 응고 된 영역에서 계산됩니다.

매우 일시적인 흐름 문제의 경우 유체와 공극 공간 사이 또는 두 개의 혼합 불가능한 유체 사이에있는 유체 인터페이스는 문제의 역학에 따라 자유롭게 움직여야합니다.

한 가지 해결책은 인터페이스와 함께 변형되는 메시를 만드는 것입니다. 이것은 시뮬레이션 중에 인터페이스의 형태가 거의 변경되지 않는 상황에서 잘 작동합니다. 그러나보다 일반적인 경우에는 시뮬레이션 중에 새 메시를 반복적으로 생성해야하거나 변경되지 않은 메시 내에서 자유 표면 경계를 생성하는 방법이 필요합니다. 이 작업은 후자를 제시합니다. VOF (Vol-of-fluid) 함수는 자유 표면의 위치를 추적하는 데 사용됩니다. 또한이 함수는 곡률을 계산하여 표면 장력의 영향을 예측하는 데 사용됩니다.

https://www.flow3d.com/wp-content/uploads/2014/08/Moving-Boundaries-an-Eularian-Approach.pdf

Settling Tanks

FLOW3D의 침전물 제거 및 드리프트 플럭스 모델은 침전 탱크, 이차 클라리퍼 및 기타 유사한 구조물을 따라 순환하는 입자의 정착 과정을 모델링 하는 데 매우 유용한 도구입니다. 이러한 모델을 사용하면 사용자가 입자 지름, 드래그 계수 등과 같은 개별 특성을 사용하여 모델에 여러개의 구성을 입력할 수 있습니다.

표준 출력물에는 유체 내 침전물 농도와 포장된 침전물 축적 및 가까운 기하학의 단단한 표면에 누적되는 비율이 포함됩니다. 전체적으로 시스템의 추가적인 유압 특성(예:전체적인 유압 효율과 잔류 시간)은 스칼라 추적 모델의 동시 사용을 통해 쉽게 얻을 수 있다. 이동 물체(밸브나 게이트 개방의 작용 등)를 포함하여 단순하고 효율적으로 솔리드 객체를 표현하는 FLOW-3DFAVOR™방법을 사용하면 재메싱 작업에 소요되는 시간이 최소화된 기하학적 매개 변수 공간을 매우 효율적으로 연구할 수 있습니다.

FLOW-3D에서 벽, 다공성 벽 또는 플럭스 측정 표면의 역할을 할 수 있는 완화된 객체는 이러한 유형의 시뮬레이션에 사용할 수 있는 모델링 접근 방식에 추가적인 유연성을 제공합니다.

3D 프린팅 시뮬레이션 기술 자료 모음

본 자료에서는 분말 또는 와이어를 층별로 적층 제조하는 3D프린팅 과정에 대해 3차원 수치해석 시뮬레이션이 가능한 FLOW-3D 제품과 기술에 관련된 자료를 찾아볼 수 있습니다. FLOW-3D는 주요 금속 적층 제조 공정인 파우더 베드 융접(PBF), 직접 에너지 증착(DED) 및 바인더 분사 공정에 대해 FLOW-3D는 높은 정확도의 시뮬레이션 기능을 제공합니다.

현재 3D 프린팅 및 레이저 용접 시뮬레이션에 FLOW-3D를 사용하고 있는 기관들은 포항공대, KAIST, 부산대학교, 성균관대학교 등 국내 여러 대학들과 한국생산기술연구원, 한국기계연구원, LG전자, 현재조선해양(현대중공업) 등 많은 기관에서 연구개발에 사용되고 있습니다.

FLOW-3D는 자유수면(VOF), 상변화, 열전달 등 여러 면에서 탁월한 정확도를 가지고 있는 제품으로 특히 국내외 적층 제조 분야에서 독보적인 위치를 차지하고 있습니다.
아래에 3D 프린팅 관련 기술을 연구 개발하는데 참고가 될 만한 관련 자료 및 링크를 제공합니다.

  1. Weld
  2. Additive Manufacturing
  3. DEM(Discrete Element Method)