J. Ocean Eng. Technol. Search

CLOSE


J. Ocean Eng. Technol. > Volume 32(3); 2018 > Article
부산신항 건설이 진우도 주변 해저지형 변화에 미치는 영향 - 수치실험을 통한 정성적 평가

Abstract

In this study, a qualitative evaluation of the morphological changes in the seabed surrounding Jinwoo-Island due to the construction of the new Busan port were determined through a numerical simulation. Various scenarios for the discharge of the Nakdong river estuary dam and construction stage of the new Busan port were established and utilized for an indirect and qualitative investigation through simulation using the numerical model implemented in this study. It was concluded through a qualitative study that the morphological changes in the seabed surrounding Jinwoo-Island were typical estuary seabed changes due to the discharge of the Nakdong river estuary dam and waves from the open sea. The effects from the construction of the new Busan port were relatively small.

1. 서 론

낙동강 하구 울타리섬 주변해역의 지형변화는 수리역학적 환경(하구둑 방류, 조석, 파랑 등)에 의하여 활발하게 진행되어 왔으며 여러 연구(Kim and Ha, 2001; Kim at al., 2005; Jang and Kim, 2006; Oh et al., 2010; Jeon et al., 2010)에서 확인되었다. 낙동강 하구 울타리섬의 지형변화에 영향을 미치는 인자 중 낙동강 하구둑 방류는 인위적인 요인으로 분류된다. 즉, 하계 태풍이후 토사농도가 큰 인위적인 대량 방류수와 함께 하구둑 외해로 이송된 토사가 하구둑 인근에 퇴적되고, 다시 외해파랑 및 조류에 의하여 이송되어 울타리섬 주변지역의 지형변화에 영향을 미치게 되는 역학적 특성을 확인하였다(Kim and Hong, 2008).
부산항 신항만 건설은 자연 상태의 지형을 변화시키기 때문에 이로 인한 수리역학적 환경 변화(흐름의 양상과 크기, 파고 등)가 불가피하게 발생된다. 이러한 수리역학적인 변동은 퇴적물의 이동현상을 통하여 기존의 평형을 이루던 자연지형 및 저질 상태로부터 새로운 평형 상태를 이루게 하는 원인이 된다. 따라서 신항만 건설로 인하여 인근 해역인 낙동강 하구 울타리섬 영역에서 국지적인 침식 또는 퇴적 현상이 예상된다.
부산항 신항만 건설 이후 진우도 협수로의 지형변화로 소형 선박의 항행에 불편을 겪었다는 민원이 꾸준히 현재까지 제기되어 오고 있다. Yonn et al.(2005)은 진우도 주변의 단기간(3개월)의 퇴적환경 변화를 진우도 주변 해안선 변화를 통하여 확인하였으며, Kim et al.(2008)은 2008년 지형(낙동강 하구둑부터 부산신항을 포함하는 영역)을 적용하여 연간 진우도 주변의 지형 변화를 수치시뮬레이션을 통하여 조사하였으며, 협수로의 지형 변화 양상을 확인하였다.
진우도를 포함한 낙동강 하구지역은 장기 및 단기적인 지형 변화 폭이 큰 지역으로, 장기지형변화를 정량적으로 정확하게 평가하기 위해서는 정확하고 다양한 입력 자료, 즉 조사 과업기간 내의 수심 분포 및 변화, 과업대상해역의 외력인 바람, 파랑(평상시, 이상시), 조류 이외에도 해저질의 분포 및 물리적 특성, 지층분포 및 최상층두께, 낙동강으로부터 유입되는 유량, 토사 성분 및 농도의 시계열자료 등이 필요하다. 그러나 충분한 자료 확보가 이루어지지 않은 상태로 정확한 정량적 평가는 무리라 판단하여 Hong(2010a)은 정성적 평가를 실시한 바 있다. 즉, 부산항 신항의 건설로 인해 발생하는 조석 및 파랑에 의한 흐름장(해저토사의 이송외력)의 변화로 인한 지형변화 차이만을 평가한 것으로서, 대상해역을 진우도 주변에 한정하였으며, 미비 자료는 경험치나 주변자료로부터 유추, 적용하여 가능한 한 실제적인 장기지형 변화를 나타내고자 하였다. 즉, 부산항 신항 건설 전과 2008년의 수심도에 각각 동일한 인자, 즉 동일한 파랑 및 조류, 동일한 해저질 분포 및 물리적 특성, 동일한 지층의 최상층 두께 및 지질, 동일한 유입유량, 토사성분 및 농도자료를 적용하여 각각 1년간의 지형변화를 시뮬레이션하였으며, 진우도 인근해역의 지형변화는 신항 건설로 인한 조류 등의 흐름 변화에 기인했다기보다는 외해파랑에 의해 발생하는 과정 중에 있는 것으로서, 전형적인 외해파랑에 의한 지형변화 중에 있는 것으로 판단하였다
상기의 정성적 평가의 문제는 진우도 주변 협수로 지형변화의 주된 원인이 외해파랑이라는 것을 추론한 것으로서, 이미 기술하였듯이 부산항 신항 건설 전과 2008년의 수심도에 동일한 해양환경을 적용한 것으로 10여년에 걸친 부산항 신항 건설기간을 고려하지 않았다. 그러므로, 본 연구의 목적은 건설 이전부터 2016년까지의 단계별 수심자료를 최대한 적용하여 부산항 신항 건설 단계별로 진우도 주변해역의 지형변화를 정성적으로 좀 더 정확하게 조사하였으며, 여러 영향 인자 중 낙동강 하구둑 방류 영향과의 상대적 비교를 위하여 방류유무에 따른 영향도 조사하였다.

2. 모델 개요

본 연구의 모델은 Hong et al.(2008)에 의해 개발된 해수유동모델, Hong(2010b)에 의해 개발된 이류⋅확산 및 퇴적물 이송 모델에 파랑모델인 SWAN(Simulating waves nearshore)에 적용된 근해풍파 모델을 결합한 것이다.

2.1 해수 유동

해수유동 모델은 수심평균된 2차원 비선형 지배방정식을 2차항 정밀도를 가진 음해 ADI유한차분법(Implicit ADI finite difference scheme of second order accuracy)을 이용하여 계산하여, 해수유동 현상을 재현하는 모델이다. 지배방정식은 식(1)-(3)과 같다.
(1)
ζt+px+qy
(2)
pt+x(p2h)+x(pqh)+ghζx+gpp2+q2C2h21ρω[x(hτxx)+x(hτxy)]ΩqfVVx+hρωx(Pα)=0
(3)
qt+y(q2h)+x(pqh)+ghζy+gqp2+q2C2h21ρω[y(hτyy)+x(hτxy)]ΩpfVVy+hρωx(Pα)=0
상기식에서,
h(x, y, t) : 수심 [m]
ζ(x, y, t) : 해수면변화 [m]
p, q(x, y, t) : [m3/s/m] = (uh, vh) 방향 유량
(u, v) = 수심평균된 xy방향 유속
C(x, y): chezy저항계수 [m1/2/s]
C= M⋅h1/6 M : Manning No.
g: 중력가속도 [m/s2]
f(v): Wind 마찰계수
V, Vx, Vy(x, y, t) : x, y 풍속성분 [m/s]
Ω(x, y): Coriolis parameter (s-1)
Pa(x, y): 대기압 [kg/m/s2]
ρω: Density of water [kg/m3]
τxx, τxy, τyy: 유효전단응력 성분

2.2 이류 확산 및 퇴적물 이송

상기의 유동모델에 추가하여 용해되어 있거나 부유되어 있는 물질의 이류 및 확산을 예측하기 위하여 다음의 식 (4)와 같은 이류-확산식을 지배방정식으로 한다. 여기서, 식 (4)는 물질의 이류항과 확산항 뿐만 아니라, 물질의 용출 및 흡입에 대한 항 또한 포함하고 있다.
(4)
t(hC)+x(uhC)+y(vhC)+z(whC)=x(hDxCx)+(hDyCy)+(hDzCz)FhC+S
여기서, C : 혼합농도
u, v, w : x, y, z 방향 속도 성분 [m/s]
h: 수심 [m]
Dx, Dy, Dz: x, y, z 방향에서 확산계수 [m2/s]
F: 선형감소계수 [l/s]
S: Qs(Cs-C)
Qs: 용출 또는 용입량 [m3/s/m2]
Cs: 용출 또는 용입 유량에서의 혼합농도
위의 식 (4)에서 uv, w는 각각 x, yz 방향으로의 유속으로 유동결과로부터 그 값을 얻는다. 확산계수 Dx, Dy, Dz 등은 등방난류확산과 수심평균화에 따른 분산을 포함한다.
식 (4)를 수치해석적으로 풀어내기 위하여 3차 양해유한차분법인 Ultimate 방식(Leonard, 1991)을 사용하였다. 이 수치해석 기법은 Quickest 방법으로부터 변형된 것으로 오염물질의 이동 및 확산을 예측하는 데 많이 사용된다.
비점착성 퇴적물의 침강속도는 Van Rijn(1993)의 방법에 따라 계산된다. 즉, 부유물 또는 퇴적물의 직경에 따라 식 (5)-(7)이 적용된다.
(5)
ws,0()=(s(l)1)gds()218v,65μm<ds100μm
(6)
ws,0()=10νds[(1+0.01(s(l)1)gds()3v2)0.51],100μm<ds100μm
(7)
ws,0()=1.1[(s()1)gds()]0.5,1000μmds
상기 식에서,
s(ℓ) 퇴적물 (ℓ) 상대밀도 ρs(ℓ)/ρW
ds(ℓ)퇴적물 직경 (ℓ)
v 물의 동점성계수 [m2/s]
소류사량의 이송크기와 방향을 포함하기 위한 시뮬레이션을 위하여 Van Rijn(1993)이 개발한 근사법이 적용되어진다(식 (8)-(11)).
(8)
|Sb|=0.006ρswsd50()M0.5Me0.7
여기서, Sb = 소류사량 [kg/m/s]
M= Wave와 Current에 의한 퇴적물 유동수
Me= 초과된 퇴적물 유동수
(9)
M=νeff2(s1)gd50
(10)
Me=(νeffνcr)2(s1)gd50
(11)
νeff=νR2+Uon2
여기서, vcr는 초기운동을 위한 한계수심평균속도(m/s)이며, 은 같은 수심평균속도의 크기로 바닥층에 속도 또는 가상로그 속도 분포로 부터 계산된다.

2.3 근해 풍파 파랑

본 연구에서 적용한 근해 풍파 모델은 네덜란드에서 개발되어 공개된 SWAN으로 근해지역에서 단주기 Short crested 파랑들의 전파, 생성과 소멸을 산정하는 모델이다. 이 모델은 수심변화에 의한 굴절과 천수 효과, 그리고 바닥마찰과 파랑의 붕괴에 인한 에너지 분산을 포함한다. 또한 이 모델은 파랑과 해류의 상호작용에 의한 효과도 포함하고 고정된, 매개변수의 방향이 분리된 모델이다. 해류의 영향을 참조하기 위해서, 이 모델에서의 기본 방정식들은 스펙트럼 파동 밀도에 관한 보존방정식을 사용한다. 보존방정식의 매개변수화는 독립변수로서 파동 스펙트럼의 영번째와 첫번째 모멘트를 사용함으로써 얻어진다.
이것은 다음과 같은 두 개의 편미분 방정식으로 식(12)-(13)과 같다.
(12)
(cgxm0)x+(cgym0)y+(cθm0)θ=T0
(13)
(cgxm1)x+(cgym1)y+(cθm1)θ=T1
여기서, m0(x, y, θ) : 파랑작용 스펙트럼의 0차 모멘트
m1(x, y, θ): 파랑작용 스펙트럼의 1차 모멘트
cgxcgy : 각각 군속도 cgx,y방향의 성분들
cθ: θ 방향의 전파속도
xy : 직교 좌표계
θ : 파랑전파의 방향
T0T1 : Source
모멘트 mnθ는 다음과 같다.
(14)
mn(θ)=0ωnA(ω,θ)dω
여기서, ω는 절대 주파수이고 는 스펙트럼 파장 작용 밀도이다. 전파속도 cgx, cgycθ는 선형파 이론에 사용하여 구해진다. 기본방정식의 좌측항은 굴절과 천수 효과를 고려한 것이다. T0T1은 바닥마찰과 쇄파에 기인한 국부적인 바람의 발생과 에너지 손실을 고려한 것이다. 이런 현상들에 대한 흐름의 영향들도 내포되어 있다.
기본방정식들은 Eulerian 유한차분법을 이용함으로써 풀며 파동 스펙트럼의 영 번째와 첫 번째 모멘트는 다수의 분리된 방향들에 대한 사각형의 격자에서 계산되어진다. Once- through marching procedure는 파랑전파의 주된 방향에 대해 적용한다. 모델로부터 유의파고, 평균파주기, 평균파향, 방향표준편차, 잉여응력과 같은 적분된 파랑 매개변수가 얻어지며 방향을 가진 파랑에너지의 분산 형태의 스펙트럼 출력 데이터도 얻을 수도 있다.

3. 정성적 평가 비교 방안 수립

부산항 신항 건설 단계별 진우도 주변 지형변화 및 낙동강 하구둑 방류 영향과의 상대적 비교를 위하여 정성적 평가방안을 수립하였다.
우선, 부산항 신항 건설 단계별에 따는 진우도 주변 지형변화에 미치는 영향을 조사하기 위하여 주된 인자를 건설 이전부터 2016년까지의 단계별 수심자료로 보았다. 그러나 확보할 수 있는 단계별 수심자료는 Fig. 1에 나타낸 A영역에 한정되어 있으므로 B영역의 해저지형은 일정 시점(신항 건설 전)에서의 수심으로 고정하고 A영역의 해저지형 변화만 반영하여 시뮬레이션을 수행하였다. A영역의 해저 지형은 Table 1의 주요 5단계를 고려(신항 건설 전, 투기장 조성, 북컨 조성, 준설 단계, 현상태)하며, 외해 파랑 및 조류 조건, 홍수기 하천 유입 토사 및 공사 중 부유토사는 가능한 한 많은 자료를 확보하여 각 단계별로 동일한 입력 자료를 사용하였다. 장기 지형변화 시뮬레이션을 수행한 후, 각 단계별 진우도 주변 해저지형의 침⋅퇴적 특성을 정량적으로 상대 비교함으로써 부산항 신항 건설이 진우도 주변 협수로 지형에 미치는 영향을 검토하고자 하였다.
앞서 언급한 동일한 모델을 적용하여 하구둑 방류만을 제거하고 수행한 결과와 하구둑 방류를 포함하여 얻은 결과를 정량적으로 상대 비교함으로써 하구둑 방류의 영향도 조사하였다. 이 경우에는 부산항 신항 건설 전의 수심도는 1980년 발행 해도를 기준으로 수립하였고, 2008년 수심도는 2006년 발행 수치 해도를 기본으로 2008년까지의 변화를 반영하여 수립한 것으로 그림 Fig. 2에 나타내었다. Fig. 2(a)Fig. 2(b)를 비교하면 낙동강 하구지역의 지형변화가 활발히 일어나고 있음을 알 수 있다. 즉, 그림의 눌차만 해역은 부산항 신항 건설로 인해 협수로 폭이 급격히 감소하였으며, 진우도는 북동쪽 일부 해안 형상이 변하며 외해방향으로 진우도 전체가 전진하는 양상을 보여준다. 또한, 신자도 A, B구역은 병합되어 신자도가 하나의 섬으로 이어지며 장자도, 대마등 및 백합등은 형상 및 크기가 변하고, 백합등 남쪽에 도요등이 새로 생성된다.
이외에도, 주목할 만한 현상으로 Fig. 2(a)의 낙동강 하구둑 근처 C구역의 모래톱이 Fig. 2(b)에서 보듯이 A와 B구역으로 분리되어 지류를 형성하게 되어 낙동강 하구 주변의 하구둑 방류에 의한 유속장이 건설 전에 비하여 차이가 나게 됨을 추측할 수 있다.

4. 수치모델 수립 및 시뮬레이션 시나리오 작성

4.1 모델 영역, 격자체계 및 외해 경계조건

Table 2에서 보듯이 파랑의 경우 세부역에 3가지 종류의 Dynamic nesting 격자시스템을, 조석의 경우 3가지 종류의 Dynamic nesting 격자시스템을 적용하되, 각 격자시스템간의 자료교환은 자동경계전이 방식(Dynamic nesting method)을 적용하여 시뮬레이션 시간을 단축하고 결과의 정확성을 증가시켰다. 이때, 자동경계전이(Dynamic nesting)는 경계에서의 유속장을 자동전이(Dynamic transfer) 할 수 있게 하였다. 앞서 Table 2에서 기술하였듯이 Grid 간격이 일정한 Dynamic nesting 격자시스템을 사용하였으며, 외해경계조건으로 조석모델의 경우 Hong(2010a)의 광역격자체계에서 얻은 결과로부터 산정하여 사용하였다, 즉, 조화상수를 적용하여 만든 시간이력조위에 항류에 의한 추가조위를 합산하여 얻었다. 파랑모델의 경우는 해양연구원의 Hinding casting 모델로 부터 Table 1의 각 시점의 자료를 활용하여 계산한 결과를 12시간 간격으로 선취하고, 각 지점간의 파랑조건은 보간법을 적용하여 사용하였다.

4.2 입력자료 및 시뮬레이션 시나리오

낙동강 하구둑에서 방류되는 유량은 Fig. 3Fig. 4에서 보듯이 2015년 수자원공사 자료를 사용하여, 하구둑 방류 구간 전체에 균일하게 방류하는 것으로 가정하였으며, 방류 SS농도는 Hong(2010a)에서 적용했던 토사농도를 사용하였다. 이외에도 녹산배수펌프장 및 서낙동강 지류로부터 유입되는 유량은 낙동강 하구둑 방출유량의 10%로 보았으며, 해저질은 신항내 점토와 실트 혼합지역으로, 이외 지역은 사토로 보았으며, 표층두께분포는 2005년 부산신항 남컨(2-2 단계) 진우도 전면해상 지반조사(동아지질) 자료와 Hong(2010a) 시뮬레이션 결과로부터 유추하여 세부역에 대하여 Fig. 5와 같이 수정함으로써 장기지형 변화의 지역적인 패턴을 반영하였다.
배경농도는 연결잔교 내외측 관측(2013.04-2014.02) 부유사농도 관측치중 Fig. 6(a)에서 보듯이 10ppm~60ppm사이의 농도중 1.8×10e-6 이상 Frequency를 제거하여 사용하였으며(Fig. 6(b)), 해저질 입경 중 D50 분포는 금번 연구에서 관측한자료 및 수자원공사 낙동강 하구둑 인근 해저질 자료를 바탕으로 정성적 추정하여 Fig. 7과 같이 작성하였다.
Table 3은 시뮬레이션 시나리오별 하구둑 방류유무 및 적용수심을 나타낸 것으로, 크게 낙동강 방류 유무 및 부산항 신항 건설 단계별로 나누어진다. 부산항 신항 건설 단계별의 경우는 전술한 부산항 신항 영역(Table 2의 A영역)의 지형 변화에 따른 B영역의 1년간의 지형 반응을 조사하였다.
하구둑 방류유무의 경우는 Hong(2010a)이 수립하여 적용했던 모델에 하구둑 방류만을 제거하여 하구둑 방류에 의한 지형반응만을 조사하기 위한 것으로, 부산항 신항 건설 단계별의 경우와는 전혀 다른 입력인자를 가지고 수립된 모델이다. 즉, 주요 인자중 하구둑 방류유량 적용시기가 부산항 신항 건설 단계별의 경우와는 관계가 없는 경우이다.

5. 진우도 협수로 지형에 미치는 영향조사 및 분석

5.1 유속장의 변화

부산항 신항 건설 및 낙동강 하구둑 방류가 진우도 주변 지역의 지형변동에 미치는 유속장(파랑에 의한 유동장 포함)의 변화를 조사하기 위하여 Fig. 8에서 보듯이 눌차만 협수로 인근 3개 지점 및 진우도-신자도 사이 중간 3개 지점의 연간 절대평균유속의 크기를 비교하여 Table 4에 나타내었다. Table 4에서 보듯이 신항 건설단계에 따른 유속의 변화는 눌차만 북안 협수로 일부 지역에서 약간 변화할 뿐 다른 정점에서는 거의 미미한 것으로 나타났다. 그러나 낙동강 하구둑 방류유무에 따른 유속의 변화는 신항 건설 전의 지형을 적용하였을 경우는 차이가 거의 나지 않는데 반하여 2008년 지형을 적용하였을 경우는 진우도와 신자도 사이의 협수로 지역에서는 뚜렷하게 차이가 남을 알 수 있었다. 이유는 전술한 것처럼 낙동강 하구둑 근처 모래톱이 분리되어 지류를 형성하게 됨으로써 낙동강 하구 주변의 하구둑 방류에 의한 유속장이 차이가 나는 것으로 판단되었다.
결론적으로, 진우도 주변해역에서의 부산항 신항 건설로 인한 장기 유속장의 변화는 미미하며, 낙동강 하구둑 방류유무에 의한 장기유속장의 변화 또한 울타리섬 인근 지형의 형상에 따른 지류 형성 등에 의하여 차이가 나는 것으로 판단되었다. 본 연구에서 조사 분석한 절대 연평균 유속의 변화는 장기지형의 변화에 미치는 영향의 정도를 단순하게 판단하는 기준으로 사용할 수 있으나, 국부적인 장기 침⋅퇴적 양상을 나타내지는 못한다. 그러므로 다음절에서 시뮬레이션 시나리오별 장기(연간) 침⋅퇴적 양상을 조사하였다.

5.2 부산항 신항 건설단계별 침⋅퇴적 변화

Fig. 9에 나타낸 A, B 및 C영역에서의 부산항 신항 건설단계별 침식 및 퇴적량을 Table 5에, 1단계 대비 침식 및 증감비율을 Table 6에 나타내었다. 부산항 신항 건설로 인하여 건설단계가 진행될수록 대부분 퇴적이 증가하는 경향을 보이며, 그런 경향은 A영역이 가장 크게 나타나며, 다음으로 B영역, C영역의 순서로 작게 나타났다. 즉, 부산항 신항에서 멀어질수록 영향이 줄어듦을 알 수 있었다.
영역내 지역별 침⋅퇴적 양상을 조사하기 위하여 1단계 대비 침⋅퇴적 증감 두께분포를 Fig. 10부터 Fig. 12에 나타내었다. Fig. 10는 A영역의 1단계 대비 증감 두께분포를 각 단계(2단계부터 5단계)별로 나타낸 것으로 ○표식 지역은 퇴적우세영역, □표식 영역은 침식우세영역으로 눌차만 입구 남측 해안의 협수로 지역에 주로 퇴적이 발생하는 것으로 나타났다. 반면, 눌차도와 진우도 사이는 침식이 주로 발생하며 눌차도 주변 해안에 집중되는 경향을 띠고 있다. Fig. 11은 B영역의 1단계 대비 증감 두께 분포를 각 단계(2단계부터 5단계)별로 나타낸 것으로 건설이 진행됨에 따라 거의 변화가 나타나지 않다가 5단계에 접어들어 진우도 남측 전면(○표식 지역)에 약간의 지형변화가 나타났다.
Fig. 12는 C영역의 1단계 대비 증감 두께분포를 각 단계(2단계부터 5단계)별로 나타낸 것으로 4단계부터 5단계에 와서 ○표식 지역에 퇴적이 집중되어 나타난다. 즉, 진우도 동측 및 신자도 서측 끝부분에 국부적으로 퇴적이 집중되어 나타날 뿐 다른 지역의 침⋅퇴적 양상에는 거의 변화가 없었다. 결론적으로, 부산항 신항 건설이 진우도 주변 지형에 미치는 영향은 눌차만 인근 지역에서만 나타나며, 진우도 및 신자도 주변에서는 극히 일부분(진우도 동측 및 신자도 서측 끝 해안)에 국부적으로 나타나는데 이는 앞 절에서 기술한 바와 같이 진우도 주변해역에서의 부산항 신항 건설로 인한 장기 유속장의 변화는 크지 않다는 결론과 일맥 상통하는 결과이다.

5.3 낙동강 하구둑 방류유무에 따른 침⋅퇴적 변화

낙동강 하구둑 방류유무에 따른 침⋅퇴적 변화를 조사하기 위하여 Hong(2010a)에 수립하여 적용했던 모델 결과와 동일 모델에 하구둑 방류만을 제거하여 하구둑 방류에 의한 지형반응만을 조사하였다. 2010년 수립된 모델의 경우 수심 및 지형은 부산항 신항 건설 전과 2008년 수심 및 지형을 적용하였으므로, 본 과업에서도 동일한 수심 및 지형에 하구둑 방류만을 제거하여 방류를 포함하였을 경우와 비교하였다. Fig. 13은 부산항 신항 건설 전 및 2008년 경우의 지형 및 침⋅퇴적 변화 조사를 위한 3개 영역(A, B 및 C영역)을 나타내었으며, A, B 및 C영역에서의 방류유무에 따른 침식 및 퇴적 변화 비교를 Table 7에 나타내었다.
Table 7에서 보듯이 신항 건설 전과 2008년의 경우 경향상 뚜렷한 차이는 나지 않으나, 방류로 인한 영향은 C영역에서 지배적으로 나타났다. 세부영역별 방류유무에 따른 침⋅퇴적 특성을 조사하기 위하여 Fig. 14부터 Fig. 16에 방류의 경우 미방류 대비 토사 증감 두께 분포도를 신항 건설 전 지형 및 2008년 지형의 각각의 경우에 대하여 나타내었다.
A영역에서는 방류로 인해 전반적으로 지형변화가 발생하며, 특히 2008년 지형을 적용한 경우가 좀 더 크게 발생함을 알 수 있다. B영역에서도 A영역과 유사한 경향을 보이며 A영역에 비하여 영향이 좀 더 증가함을 알 수 있다. 그러나 C영역에서는 방류로 인한 영향이 신항 건설 전 지형을 적용한 경우나 2008년 지형을 적용한 경우 둘 다 상당히 크게 미침을 알 수 있으며, 신항 건설 전은 진우도와 신자도 북쪽 영역에서, 2008년의 경우는 진우도와 신자도 사이 및 남쪽 영역에서 방류로 인한 영향이 두드러짐을 알 수 있다.
결론적으로 낙동강 하구둑 방류로 인한 영향은 눌차만 주변 지형에 비해 진우도와 신자도 주변지형에 지배적인 영향을 미치며, 파랑 및 방류 등에 의한 낙동강하구 지형의 전반적인 지형변화가 유속장의 변화를 발생하고, 이것이 방류로 인한 지형변화 양상에도 영향을 미침을 알 수 있다.

5.4 신항 건설 영향과 낙동강하구둑 방류 영향의 정성적 비교

Table 3에서 보듯이 본 과업에서 수행한 시뮬레이션은 ①울타리섬 영역은 부산항 신항 건설 전의 지형을 적용하면서 부산항 신항 영역만 건설단계별로 변화시켜 진우도 주변 지형변화를 조사하는 것과, ②과업 전체영역을 부산항 신항 건설 전과 2008년 지형을 적용하여 낙동강 하구둑 방류유무에 따른 진우도 주변지형 변화를 조사하는 것으로 ①의 경우와 ②의 경우 시뮬레이션의 주요 조건이 서로 다르므로 정량적 상대비교 대신 정성적 상대비교를 통하여 신항 건설(신항 영역은 현상태, 울타리섬 영역은 신항 건설 전)로 인한 영향과 낙동강 하구둑 방류 영향을 정성적으로 비교하였다.
Fig. 17부터 Fig. 19는 A, B 및 C영역에서 신항 건설단계별 영향의 시뮬레이션 결과 중에서 1단계 대비 5단계 증감 분포와 하구둑 방류의 영향으로 인한 증감 분포 비교를 나타낸 것이다. Fig. 17에서 보듯이 A영역에서는 신항 건설로 인한 영향이 하구둑방류에 의한 영향보다 우세하나, Fig. 18의 B영역에서는 신항 건설보다 하구둑 방류의 영향이 진우도 주변 지역에 한정되어 우세하게 나타남을 알 수 있었다.
이런 경향은 동진할수록 더 증가하여 Fig. 19에서 보듯이 하구둑 방류 영향이 상대적으로 신항 건설의 영향보다 두께 분포 변화가 상당히 많이 차이남을 알 수 있다. 즉, 진우도 주변지형의 변화에 미치는 영향은 부산항 신항 건설보다 낙동강 하구둑 방류에 의한 영향이 상대적으로 지배적이라고 판단되었다.
지금까지의 결과 분석으로부터, 진우도 주변의 지형변화는 부산항 신항 건설로 인한 영향이 적으며, 하구지역의 방류 및 외해파랑에 의한 전형적인 하구지역에서 발생하는 하구지형 변화 중에 있는 것으로 판단되었다.

6. 결 론

본 연구에서 부산항 신항 건설이 진우도 주변 지형에 미치는 영향을 조사하기 위하여 간접적 정성적인 검토를 위한 시나리오 작성 및 모델을 수립하여 시뮬레이션을 수행함으로써, 하구둑 방류유무에 따른 진우도 주변 지형 반응과 신항 건설단계별 진우도 주변 지형의 침⋅퇴적 양상을 상대적으로 비교하여 다음과 같은 결론을 얻었다.
(1) 신항 건설의 영향은 눌차만 인근 지역이 영향이 크고 진우도 주변 지역으로 동진할수록 영향이 상대적으로 적은 것으로 나타났다.
(2) 하구둑 방류의 영향은 진우도와 신자도 주변에서 크고 서진할수록 영향이 상대적으로 적은 것으로 나타났다.
(3) 신자도 동쪽 지형의 변화는 하구둑 방류로 인하여 발생하며, 낙동강 울타리섬의 지속적인 지형변화로 인한 유속장의 변화가 진우도 주변 지형의 변화를 발생시키는 것으로 판단된다. 또한, 하구둑 방류로 인한 토사공급과 이를 운반하는 파랑의 작용도 영향을 미치는 것으로 보인다.
결론적으로, 진우도 주변의 지형변화는 부산항 신항 건설로 인한 영향이 상대적으로 적으며, 하구지역의 방류 및 외해파랑에 의한 전형적인 하구지역에서 발생하는 하구지형 변화 중에 있는 것으로 판단되었다.

Fig. 1
Area for this study
joet-32-3-192f1.jpg
Fig. 2
Bathemetry before the construction of Busan New Port and 2008
joet-32-3-192f2.jpg
Fig. 3
Dischage from Nakdong estuary dam [m3/sec]
joet-32-3-192f3.jpg
Fig. 4
SS concentration of discharge from Nakdong dam
joet-32-3-192f4.jpg
Fig. 5
Seabe thickness of nested model area B [m]
joet-32-3-192f5.jpg
Fig. 6
Evaluation of background SS cocentration [mg/l]
joet-32-3-192f6.jpg
Fig. 7
Distribution of D50 grain size [mm]
joet-32-3-192f7.jpg
Fig. 8
Location for comparison of year-average current velocity
joet-32-3-192f8.jpg
Fig. 9
Area for comparison of year-erosion and deposition
joet-32-3-192f9.jpg
Fig. 10
Distribution of thickness of erosion and deposition for each construction stage to 1st construction stage at A area
joet-32-3-192f10.jpg
Fig. 11
Distribution of thickness of erosion and deposition for each construction stage to 1st construction stage at B area
joet-32-3-192f11.jpg
Fig. 12
Distribution of construction stage at C area. thickness of erosion and deposition for each construction stage to 1st
joet-32-3-192f12.jpg
Fig. 13
Area for comparison of year-erosion and deposition; (a) before construction, (b) bathemetry of 2008
joet-32-3-192f13.jpg
Fig. 14
Distribution of thickness of erosion and deposition due to discharge at A area.; (a) before construction, (b) bathemetry of 2008
joet-32-3-192f14.jpg
Fig. 15
Distribution of thickness of erosion and deposition due to discharge at B area; (a) before construction, (b) bathemetry of 2008
joet-32-3-192f15.jpg
Fig. 16
Distribution of thickness of erosion and deposition due to discharge at C area.; (a) before construction, (b) bathemetry of 2008
joet-32-3-192f16.jpg
Fig. 17
Comparison at A area between effect of port construction and effect of dischage of estuary dam.; (a) Change of thickness of erosion and deposition due to port construction (b) Change of thickness of erosion and deposition due to port construction
joet-32-3-192f17.jpg
Fig. 18
Comparison at B area between effect of port construction and effect of dischage of estuary dam.; (a) Change of thickness of erosion and deposition due to port construction (b) Change of thickness of erosion and deposition due to port construction
joet-32-3-192f18.jpg
Fig. 19
Comparison at C area between effect of port construction and effect of dischage of estuary dam.; (a) Change of thickness of erosion and deposition due to port construction (b) Change of thickness of erosion and deposition due to port construction
joet-32-3-192f19.jpg
Table 1
Five main stage for construction of Busan New Port
Stage Layout Feature
1 joet-32-3-192i1.jpg Before constructi-on
2 joet-32-3-192i2.jpg During constructi-on of dumping area of dredged soil and north dock
3 joet-32-3-192i3.jpg Starting of dredgin-g in port area
4 joet-32-3-192i4.jpg Completion of north dock and 1st dredging work
5 joet-32-3-192i5.jpg During construction of south dock
Table 2
Model area, grid system and boundary condition for wave and tide
Wave Tride
Model area joet-32-3-192i6.jpg joet-32-3-192i7.jpg
Dynamic nesting grid system A 500 m × 500 m A 100 m × 100 m
B 100 m × 100m B 50 m × 50 m
C 25 m × 25 m C 10 m × 10 m
Boundary condition From results from KORDI Hong (2010a)
Table 3
Bathemetry and dischagrge state for simulation cases
Simulation case # Bathemetry Existence of discharge of estuary dam
Busan New Port area Fencing island area
Existence of discharge of estuary dam 1-1 Before construction, bathemetry of 1996 1980 nautical chart O
1-2 2008 construction stage 2006 nautical chart O
2-1 Before construction, bathemetry of 1996 1980 nautical chart X
2-2 2008 construction stage 2006 nautical chart X

Construction stage of Busan New port 3-1 Before construction, bathemetry of 1996 1980 nautical chart O
3-2 Construction of North container and starting dumping area 1980 nautical chart O
3-3 Starting dredging in port area 1980 nautical chart O
3–4 Finishing 1st dredging and north container area 1980 nautical chart O
3–5 during south container 1980 nautical chart O
Table 4
Year-average current velocity at location of comparison of year-average current velocity for simulation cases
Year average current speed for each construction stage [m/sec] Year average current speed to existence of discharge of estuary dam [m/sec]
Bathemetry before construction of port Bathemetry of 2008

1st 2nd 3rd 4th 5th with without with without
m1 0.01 0.01 0.02 0.02 0.01 0.08 0.09 0.12 0.13
m2 0.04 0.04 0.05 0.06 0.08 0.25 0.25 0.22 0.24
m3 0.33 0.34 0.25 0.27 0.28 0.25 0.25 0.15 0.15
s1 0.21 0.21 0.21 0.21 0.20 0.18 0.18 0.37 0.49
s2 0.20 0.19 0.20 0.19 0.19 0.19 0.19 0.46 0.58
s3 0.17 0.17 0.17 0.17 0.16 0.18 0.21 0.38 0.47
Table 5
Erosion and deposition rate for construction stage
Area Erosion and deposition rate for each stage [m3]
A 1st 2nd 3rd 4th 5th
Deposition 70984 62594 51393 53893 47861
Erosion −82537 −69548 −61740 −60048 −54375
Total −11553 −6953 −10347 −6155 −6514

B 1st 2nd 3rd 4th 5th
Deposition 8170 8289 8466 9210 9795
Erosion −7156 −7202 −7395 −8114 −8582
Total 1014 1086 1071 1095 1212

C 1st 2nd 3rd 4th 5th
Deposition 24614 24981 24332 26650 27612
Erosion −25232 −25617 −24949 −27335 −28322
Total −617 −636 −617 −685 −710
Table 6
Erosion and deposition rate for each construction to 1st stage
Area Erosion and deposition percent rate to 1st stage [%]
A 2nd 3rd 4th 5th
Deposition −11.82 −27.60 −24.08 −32.57
Erosion 15.74 25.20 27.25 34.12
Total 39.81 10.44 46.72 43.62

B 2nd 3rd 4th 5th
Deposition 1.45 3.62 12.72 19.88
Erosion −0.65 −3.33 −13.39 −19.93
Total 7.13 5.61 8.00 19.52

C 2nd 3rd 4th 5th
Deposition 1.49 −1.15 8.27 12.18
Erosion −1.53 1.12 −8.33 −12.25
Total −2.99 0.05 −10.91 −15.01
Table 7
Comparion of erosion and deposition rate due to dischage of estuary dam
Area Classification Bthemetry before construction Bathemetry of 2008
with [m3] without [m3] percent to without [%] with [m3] without [m3] percent to without [%]
A deposition erosion total 88933 75367 18 −149738 −118857 −26
−109003 −94483 −15 121932 98969 23
−20070 −19117 −5 −27806 −19888 −40

B deposition erosion total 157697 123776 27 −356189 −320834 −11
−349327 −319768 −9 83139 102268 −19
−191630 −195992 2 −273050 −218566 −25

C deposition erosion total 177007 182007 −3 −244214 −201062 −21
−560977 −334950 −67 181875 62566 191
−383970 −152943 −151 −62339 −138496 55

References

Kim, SY., & Ha, JS. (2001). Sedimentary Faces and Environmental Changes of Nakdong River Estuary and Adjacent Coastal Area. Korean Journal of Fisheries and Aquatic Sciences, 34(3), 268-278.

Kim, BW., Kim, BG., & Lee, SY. (2005). Development of Mosaic Aerial Photographs for Shoreline Change Study in Nakdong Estuary. Ocean and Polar Reasearch, 27(4), 497-507.
crossref pdf
Jang, ST., & Kim, KC. (2006). Change of Oceanographic Environment in the Nakdong Estuary. The Sea, 11(1), 11-20.

Oh, CY., Park, SY., Choi, CU., & Jeon, SW. (2010). Change Detection at the Nakdong Estuary Delta Using Satellite Image and GIS. Journal of the Korean Society for Geospatial Information Science, 18(1), 21-29.

Jeon, YH., Yoon, HS., & Lee, KS. (2010). Spatiotemporal Variations of Coastal Sediment Transport at Barrier Islands in the Nakdong River Estuary. Journal of the Korean Society of Marine Environment and Safety, 16(2), 161-168.

Kim, KS., & Hong, NS. (2008). The Change of Ocean around Estuary Caused by Change of Shape of Nokdong River Estuary a Sandy Plain.. Proceedings of the 17th Conference of the Korean Society of Coastal and Ocean Engineers Mokpo: 100-103.

Kim, KS., Hong, NS., Chae, JW., & Lee, JM. (2008). Morphological Simulations near Coastal Area of Jinwoo-do.. Proceedings of 12th Conference of the Korean Association of Ocean and Technology Societies Jeju: 374-377.

Yoon, HS., Yoo, CI., Ryu, CR., & Lee, IC. (2005). Sediment Environmental Change around Jinwoo-do in Nakdong River Estuary. Proceedings of Spring Conference of the Korean Society of Marine Environment and Energy, 44-49.

Hong, NS. (2010a). The 4th Project Report of Research and Development for Ocean Environmental Hyraulics of Busan New Port - Spreading of Suspended Sediment and Morphological Change. Busan Port Construction Office, Miinistry of Ocean and Fisheries.

Hong, NS. (2010b). Investigation of Reducing Characteristics for the Spreading of Dredging Soil and the Diffusion of Contaminant by Silt Protector Curtain through Three Dimensional Numerical Model Experiment. Journal of Ocean Engineering and Technology, 24(4), 78-85.

Hong, NS., Kim, GY., & Kang, YG. (2008). Three Dimensional Numerical Model for Flow with Silt Protector. Journal of Ocean Engineering and Technology, 22(3), 1-7.

Leonard, BP. (1991). The Ultimate Conservative Difference Scheme Applied to Unsteady One-Dimensional Advection. Computer Methods in Applied Mechanics and Engineering, 88(1), 17-74.
crossref
Van Rijn, LC. (1993). Principles of Sediment Transport in Rivers, Estuaries and Coastal Seas.. AQUA Publications: The Netherlands.

TOOLS
Share :
Facebook Twitter Linked In Google+ Line it
METRICS Graph View
  • 0 Crossref
  •   Scopus
  • 10,431 View
  • 117 Download


ABOUT
BROWSE ARTICLES
ARTICLE CATEGORY

Browse all articles >

PUBLICATION ETHICS
FOR CONTRIBUTORS
Editorial Office
President Office BD Rm. 1302, 13 Jungang-daero 180beon-gil, Dong-gu, Busan 48821, Republic of Korea
Tel: +82-51-759-0656    Fax: +82-51-759-0656    E-mail: ksoehj@ksoe.or.kr                

Copyright © 2024 by The Korean Society of Ocean Engineers.

Developed in M2PI

Close layer
prev next