2.1 자료

연구해역에 대한 수심과 조석 및 조류자료, 수온 및 염분자료, 하천유량 등의 관측자료를 수집하여 모델의 입력자료와 보정 및 검정을 위한 자료로 이용하였다(Fig. 1). 모델에 입력한 수심은 국립해양조사원에서 간행한 해도와 시화호에서 측량한 수심자료를 이용하였다(Fig. 2와 Fig. 3). 모델 입력자료의 기준년도는 2010년으로 하여 가용한 자료를 수집하여 사용하였다.

Fig. 1.

Computational grid system of model domain.

Fig. 2.

Bathymetry of the study area in Gyeong-Gi Bay. Squares and circles represent observation points of tide and tidal current respectively.

Fig. 3.

Detailed depth contour and observation points of temperature and salinity around Shihwa Lake.

해수유동모델의 개방경계와 보정 및 검정을 위한 관측자료는 인근 해역의 조석 4개 정점(PT1, PT2, PT3, T1)과 조류 4개 정점(PC1, PC2, PC3, PC4) 자료를 이용하였다(Fig. 2와 Table 1).

Table 1.

Observed tide and tidal current stations

조석은 인천(PT1), 안산(PT2), 평택(PT3)의 고정 조위관측소에서 2010년에 관측한 자료와 시화방조제(T1, 2002.4~2003.4)에서 기존에 관측한 자료를 이용하였다. 조류는 시화방조제 전면 수로에서 2009년 12월과 2010년 1월에 3개점(PC1, PC2, PC3)에서 관측한 중층유속 자료와 송도신도시 인근 해역에서 기존에 관측한 층별유속 자료(PC4, 2004.12~2005.1)를 이용하였다.

수온 및 염분자료는 시화호 해양환경개선사업(Ministry of Land, Transport and Maritime Affairs, 2010)에서 관측한 시화호 내측자료(Fig. 3)와 국가해양환경측정망 자료 중 인천, 아산, 가로림 연안 자료를 이용하여 초기조건을 구성하였다. 개방경계의 수온과 염분은 연안정지관측 자료 중 소청도와 선미도의 일 관측자료와 정선해양관측 자료 중 307 line의 3번 정점의 연직자료를 이용하여 연간 일 자료를 적용하였다(National Fisheries Research and Development Institute, http://portal.nfrdi.re.kr). 또한, 하천경계는 한강, 임진강, 예성강과 시화호 내측 21개 소하천을 고려하였는데, 한강의 유량은 행주대교, 임진강의 유량은 적성지역의 2010년 자료를 각각 적용하였다. 예성강의 유량은 임진강의 유량에 유역면적비를 적용하여 환산하였으며, 유입되는 모든 하천의 염분은 0psu로 입력하였다. 한강, 임진강, 예성강의 담수 유입시 적용한 수온자료는 환경부의 월별 하천수 수질현황 자료를 이용하였다.

시화호 내측으로 유입되는 하천 자료는 유역모델로 산정된 21개 지점의 유량을 입력하였다. 시화호 내측으로 유입되는 각 하천은 Table 2에 정리하였고, 하천의 위치와 유역모델을 통해 산정한 유량은 Fig. 3과 Fig. 4에 각각 제시하였다.

Table 2.

Fresh water discharge points flows into Shihwa Lake

Fig. 4.

Fresh water discharge into Shihwa Lake.

표층경계 조건으로 입력하기 위해 인천기상대와 수원기상대에서 2010년의 매 시간 기상자료 및 바람자료를 수집, 분석하였다. 기온, 운량, 기압, 상대습도, 강우량, 태양복사에너지 등은 인천기상대와 수원기상대의 자료를 평균하여 모델에 입력하였으며, 풍향, 풍속 등의 바람자료는 인천기상대의 자료를 사용하였다(Fig. 5).

Fig. 5.

Time series of meteorological data.

2.2 연구 방법

시화호 해역의 해수유동특성을 재현하기 위해 사용한 모델은 EFDC(Environmental Fluid Dynamics Code)로 모델의 설명은 Hamrick(1992)에 자세하게 제시되어 있다. EFDC를 활용한 연구결과로는 해수와 담수의 혼합과 성층현상 모의(Kang et al., 2011; Kim et al., 2008), 새만금 방조제 완공 후 유동 및 수질모의(Suh et al., 2010), 광양만의 수질모의(Park et al., 2005), 남해안의 조위 및 조류 모의(Kwon and Kang, 2007), 한강하구 및 염하수로의 순 수송량 모의(Yoon and Woo, 2012) 등이 있다.

모델의 계산영역은 시화호를 포함하여 동서방향으로 139.3km, 남북방향으로 150.2km로 설정하였고 수평방향으로는 80~2,400m의 직교곡선격자체계(orthogonal curvilinear grid system)를 채택하여 해안선과 지형을 모델에 재현하였다(Fig. 1). 연직방향으로는 5개의 sigma층으로 설정하였고, 모델의 초기조건은 해수의 흐름이 없는 상태(cold-start)로 하였으며 모델의 외해경계조건은 NAOJ(National Astronomical Observatory of Japan)의 NAO.99Jb(Matsumoto et al., 2000)로부터 유도된 조석의 주요 4개 분조(M₂, S₂, K₁, O₁)를 합성한 조위를 경계면에 부여하였다.

수치모델 실험안은 조력발전소가 가동되지 않은 상황인 2010년 관측자료를 기준으로 시화호 해역 현황(Case N)과 조력발전소 가동상황(Case A)으로 구성하였다(Table 3). 모델실험은 2010년 4월부터 11월까지 8개월간 계산하였다.

Table 3.

Model simulation cases

조력발전소 비가동시에는 시화호 남서측에 위치한 기존 갑문을 통해 해수가 유통되고, 조력발전소 가동시에는 기존갑문과 시화방조제 중앙에 건설된 조력발전소의 수문과 수차를 통해 해수가 유통된다.

시화호 기존 갑문은 8대로 구성되어 있고, 조력발전소는 10대의 수차와 8대의 수문으로 각각 구성되어 있다. 시화조력발전소의 운영은 낙조시 수문, 수차, 기존 갑문을 통해 배수하고 창조시 수차를 통해 낙차를 이용한 단류식 발전 방식으로 운영된다. 시화호 내측의 유지관리 수위는 평균해면(M.S.L.) 하 1.0m이며, 시화방조제 기존 갑문과 시화조력발전소 운영계획에 따른 유량특성식은 아래에 제시하였다(Korea Water Resources Corporation, 2005). 유량특성식에 따른 통과유량(수차, 수문 및 기존갑문)을 시화방조제 내측과 외측의 수위차 조건에 따라 모델 내부에서 방조제 내측 또는 외측 방향으로 이동하도록 적용하였다.

(1) 시화방조제 기존 갑문의 유량 특성식(수문 1대당) 시화방조제 기존 갑문의 유량특성식은 다음과 같다.

Q = C × B × H₂ × (2g∆H)^0.5

여기서,

Q: 해수 유입·출입량(m³/s)

C: 유량계수(0.8), B: 갑문폭(12 m)

H₂: 외수위(시화호측 바닥 기준, m)

g: 중력가속도(m/s²)

ΔH: 수두차(ㅣ내수위−외수위ㅣ, m)이다.

(2) 시화조력 수차의 유량 특성식(수차 1대당)

시화조력 수차의 발전시 및 방류시의 터빈 특성식을 제시하였다.

가. 발전시 터빈 특성식

Q=0                  (0.00<ΔH<1.07)

Q=150.433×ΔH+113.248179        (1.07≤ΔH<1.95)

Q=18.1077×ΔH+371.3606         (1.95≤ΔH<5.69)

Q=1767.8881-346.5186×ΔH+20.9571×ΔH²(5.69≤ΔH<7.43)

나. 방류시 터빈 특성식(수차 1대당)

Q=800×ΔH                                  (0.00<ΔH<0.05)

Q=600×ΔH+10                                  (0.05≤ΔH<0.10)

Q=250.0×ΔH+45                                  (0.10≤ΔH<0.30)

Q=59.56939963+231.1052218×ΔH−94.5018523×ΔH²+26.29460961×ΔH³−2.754795024×ΔH⁴ (0.30≤ΔH<2.50)

(3) 시화조력 수문의 유량 특성식(수문 1대당) 시화조력 수문의 유량특성식은 다음과 같다.

Q=839×ΔH^0.3855

3. 수치모델의 검정

모델의 검정을 위하여 조위 4개 정점(PT1, PT2, PT3, T1)과 조류 4개 정점(PC1, PC2, PC3, PC4), 그리고 시화호 내측과 외측을 포함한 시화방조제 인근 18개 정점(Fig. 3)에서 계절별로 관측된 수온 및 염분에 대한 모델 결과가 관측자료의 변동 특성을 잘 재현하도록 반복실험을 수행하였다.

조석의 검정은 관측된 조석의 조화상수와 계산 결과의 조화상수를 비교하였으며 모델의 정확도를 정량적으로 평가하기 위해 절대상대오차(Absolute Relative Error)와 평균오차(Mean Error)를 산정하였다.

각 정점에 대한 조석검정 결과, M₂분조, S₂분조, K₁분조, O₁분조의 진폭에 대한 절대상대오차(ARE)의 범위는 0.1~1.9%로 4개 정점에서 모두 98% 정도의 정확성을 보였다(Table 4).

Table 4.

Comparison of tidal constituents between observed and model results

조류검정은 각 정점별로 모델결과와 관측결과의 조류조화 상수로부터 산정된 U(동서성분), V(남북성분)와 합성유속을 비교하는 방법으로 수행하였다. 비교결과, 각 정점별 유속에 대한 절대상대오차(ARE)의 범위는 0.9~10.5%로 90% 정도의 정확성을 나타냈다(Table 5).

Table 5.

Comparison of tidal currents between observed and model results(unit: cm/s)

한편, 수온과 염분의 검정은 시화호 해양환경개선사업에서 관측한 시화호 내측 조사자료와 시화방조제 외측의 국가해양 환경측정망자료(inc정점)를 이용하여 계절변동 특성을 재현하였다(Fig. 3).

시화호 해양환경개선사업에서는 호 내측 12개 정점에서 2010년 4월, 6월, 7월, 8월, 10월의 5회에 걸쳐서 공간적인 표층조사가 수행되었고, 이 중 3개 정점(정점 5, 9, 12)에서는 표층과 저층에서 조사가 수행되었다. 특히 정점 9에서는 공간조사 이외에도 3월, 4월, 6월, 8월, 10월, 11월에 시기별로 3~10일간 3~6시간 간격으로 표층과 저층에서 수온, 염분조사가 수행되었다. 시화호 외측의 국가해양환경측정망 자료는 연 4회(2월, 5월, 8월, 11월)에 걸쳐 각 시기별로 1회씩 측정된 자료이다. 따라서 수온과 염분에 대한 모델 검정은 각 정점의 계절변동을 비교하는 방법으로 수행하였다.

수온과 염분의 계절별 관측 결과(Ministry of Land, Transport and Maritime Affairs, 2010)에 의하면 하계에 수온상승과 담수 유입량 증가로 인해 염분이 낮아지며, 공간적으로 수온은 정점별로 차이가 크지 않지만 염분의 경우에는 담수 유입의 영향이 큰 호 상류에서 낮게 나타난다(Table 6).

Table 6.

Observed surface temperature and salinity at each station

2010년 4월부터 11월까지 수온과 염분에 대한 모델 계산 결과를 계절별 연속측정 자료(정점 9)와 비교한 결과를 Fig. 6에 제시하였고, 시화호 내측에서 연 5회 공간 조사 자료(12개 정점)와 비교한 결과를 Fig. 7과 Fig. 8에 수온과 염분으로 각각 구분하여 제시하였다. 또한 시화호 외측 6개 정점의 관측 자료와 비교한 수온과 염분의 검정결과를 Fig. 9와 Fig. 10에 제시하였다.

Fig. 6.

Comparison of temperature and salinity between observed and model results at station 9.

Fig. 7.

Comparison of temperature between observed and model results at stations inside Shihwa Lake.

Fig. 8.

Comparison of salinity between observed and model results at stations inside Shihwa Lake.

Fig. 9.

Comparison of temperature between observed and model results at stations outside Shihwa Lake.

Fig. 10.

Comparison of salinity between observed and model results at stations outside Shihwa Lake.

정점 9에서 연속측정한 자료의 경우 8월 29일~30일에 염분이 큰 폭으로 낮아지는 현상과 모델 결과가 편차를 나타내고 있다. 관측자료에 나타난 낮은 염분은 기상자료 분석결과, 8월 29일의 강우의 영향으로 판단되며 모델에 적용한 강우조건은 인천기상대와 수원기상대의 평균 강우자료를 이용함에 따라 기상대와 시화호의 위치에 따른 편차로 인해 국지적인 강우의 영향을 모델이 잘 반영하지 못한 것으로 사료된다. 시화호 내측의 정점별 모델의 염분 시계열을 보면, 상류로 갈수록 수심이 얕아지고 유입되는 담수의 영향을 받아 염분의 변동폭이 크게 나타난다(Fig. 8). 모델결과를 관측자료와 비교 결과, 각 정점에서 8개월간 계산된 수온과 염분의 시계열은 일부 정점에서 차이를 보이지만 전반적으로 관측자료의 연변동 추세를 잘 재현한 것으로 판단된다.

이와 같은 검정과정을 통해 구축된 모델을 이용하여 시화조력발전소 운영 전과 운영 후 상황에 대한 시화호의 해수특성 변화를 분석하였다.

4.1 현황

(1) 조력발전소 비가동 상황

조력발전소가 가동되지 않는 상황(Case N)에서 기존 갑문을 통한 창조시와 낙조시의 유향과 유속을 표층과 저층으로 구분하여 Fig. 11과 Fig. 12에 각각 제시하였고, 최강유속 분포를 Fig. 13에 제시하였다. 여기서 최강유속은 대조·중조·소조기를 포함하는 모델 계산기간 동안 각 격자점에서 계산된 유속의 최대치를 나타낸 것이다.

Fig. 11.

Inflow vector plots at the surface(A) and bottom(B) without tidal power plant operation(Case N).

Fig. 12.

Outflow vector plots at the surface(A) and bottom(B) without tidal power plant operation(Case N).

Fig. 13.

Maximum current speed at the surface(A) and bottom(B) without tidal power plant operation(Case N).

계산 결과, 기존 갑문을 통해서 창조시에 유입되고 낙조시에 유출되는 것을 볼 수 있고, 시화호 외측 방조제 전면에서는 창조시 수로를 따라 동북동 방향으로 진행하고 낙조시에는 이와 반대 방향의 흐름이 나타난다. 시화호 내측의 표층 유속은 시화호 기존 갑문을 통해 해수가 유통됨으로써 갑문과 상류부 수로 부근을 제외하고는 10cm/s 이하의 약한 유속분포를 보이는 반면 시화호 갑문 인근에서는 100cm/s 이상으로 나타나며, 시화호 외측 방조제 전면 수로에서는 강한 조류의 영향으로 40~150cm/s의 유속분포를 나타낸다.

(2) 조력발전소 가동 상황

조력발전소 가동상황(Case A)에 대한 실험결과를 분석하였다. 시화 조력발전소는 8대의 수문과 10대의 수차로 구성되어 있으며, 조력발전시 수문과 수차 및 시화방조제의 기존갑문을 통해 유입·유출되는 해수에 의한 내측과 외측의 수위 및 조위 변화는 Fig. 14에 조력발전소 비가동 상황과 함께 제시하였다. 시화호 내측의 수위는 유지관리수위 평균해면(M.S.L.) 하 1.0m를 유지하고 있으며, 기존 갑문을 통해 해수가 유통되는 조력발전소 비가동 상황보다 조력발전소 가동시 수위가 큰 폭으로 변화한다.

Fig. 14.

Water level variations: Upper panel represent without tidal power plant operation and lower panel represent with tidal power plant operation.

조력발전소를 운영할 때 수차를 통한 호내측으로 해수 유입시(발전시)와 수차와 수문 및 기존갑문을 통한 유출시(배수시)에 대한 유속과 유향을 표층과 저층으로 구분하여 제시하였으며(Fig. 15와 Fig. 16), 조력발전시 최강유속의 분포를 Fig. 17에 제시하였다. 계산 결과, 발전소 가동시 해수의 유입과 유출에 따라 북서-남동방향으로 빠른 유속 분포를 나타낸다. 조력발전소 외측으로 배수시 최대 유속이 200cm/s 이상이며 시화호 내측도 발전소 가동에 따른 수위 변화폭이 증가함에 따라 유속이 증가한다.

Fig. 15.

Inflow vector plots at the surface(A) and bottom(B) with tidal power plant operation(Case A).

Fig. 16.

Outflow vector plots at the surface(A) and bottom(B) with tidal power plant operation(Case A).

Fig. 17.

Maximum current speed at the surface(A) and bottom(B) with tidal power plant operation(Case A).

4.2 결과 분석

시화호 해역의 조석, 조류, 수온, 염분에 대한 검정 과정을 거친 수치실험 결과로부터 조력발전소 가동으로 인한 해수특성의 변화를 분석하였다. 담수의 영향을 분석하기 위하여 여름철(7월~9월) 상황에 대한 결과를 이용하였다.

(1) 밀도 수평분포 변화

여름철(7월~9월)의 수온, 염분 및 밀도에 대한 평균값과 시계열자료로 부터 표층과 저층에 대한 표준편차(Standard deviation)를 구하였다.

여기서, STD는 표준편차, N은 총 자료 수, 는 N개 자료에 대한 평균이다.

시화호 내측의 밀도분포는 조력발전소 가동 전에는 표층의 경우 호 내측에서 상류로 갈수록 표준편차가 크게 나타나는데 이는 염분 표준편차에서 볼 수 있듯이 하천에서 유입되는 담수의 영향을 많이 받기 때문인 것으로 보인다. 한편, 조력발전소 가동시에는 호 내측 중앙부까지 유입되는 해수에 의해 표준편차가 작아지는 경향을 보이나 수심이 얕은 상류에서는 배수에 따른 수위 하강으로 인해 유입되는 담수의 단위 부피당 영향이 증가하기 때문에 표준편차가 소폭 커지는 경향을 보인다(Fig. 18과 Fig. 19).

Fig. 18.

Standard deviation of surface temperature, salinity and sigma-t(σ_t) in summer : Left panels represent without tidal power plant operation and right panels represent with tidal power plant operation.

Fig. 19.

Standard deviation of bottom temperature, salinity and sigma-t(σ_t) in summer : Left panels represent without tidal power plant operation and right panels represent with tidal power plant operation.

이러한 결과는 시화호 내측의 밀도분포에 영향을 미치는 주요 인자가 하천을 통해 유입되는 담수와 조력발전소 가동에 의해 유입·유출되는 해수에 의한 염분의 변화임을 알 수 있는데, 이러한 양상은 잔차류 분포에서도 나타난다. Fig. 20과Fig. 21에 잔차류에 대한 streamline과 흐름의 방향을 함께 제시하였다.

Fig. 20.

Streamlines of residual currents at the surface(A) and bottom(B) without tidal power plant operation.

Fig. 21.

Streamlines of residual currents at the surface(A) and bottom(B) with tidal power plant operation.

조력발전소 가동 전 표층의 잔차류는 호 내측으로 유입되는 담수의 영향으로 인해 하류 방향으로 호 중앙부까지 이동하는 반면에 저층은 시화호 기존 갑문으로 유입되는 해수가 표층과 반대로 상류 방향으로 이동하는 2층 순환구조(two-layer circulation)를 보인다. 기존 갑문을 통한 해수유통과 상류부에서 내려오는 하천수의 영향으로 인해 방조제 내측의 잔차류는 조력발전소를 기준으로 서측에서는 시계방향, 동측에서는 반시계방향의 흐름을 보인다.

한편, 조력발전소 가동시의 잔차류는 조력발전소를 통해 해수유통량이 증가함에 따라 염분이 증가하고, 배수시 수위하강에 따라 하류방향으로 담수의 영향이 증가하면서 시화호 중앙과 상류에서 2층 구조가 상대적으로 뚜렷해지는 경향을 보인다. 조력발전소 가동시 표층 잔차류가 외해방향으로 우세한 것은 발전소 운영시 배수로 인한 유속의 증가 때문이다. 조력발전소를 통한 해수유통의 영향으로 방조제 내측의 잔차류는 조력발전소를 기준으로 서측에서 시계방향, 동측에서 반시계방향의 흐름을 보임으로써, 조력발전소 가동 전과 유사한 경향을 보이나 발전시 수차를 통한 유입효과가 크게 나타난다.

(2) 밀도 수직구조 변화

조력발전소 가동으로 인한 하계의 평균밀도 수직구조 변화를 파악하기 위해, 시화호 내측의 6개 정점(St. 5, 7, 9, 12, 13, 15)에서 모델을 통해 계산된 수온과 염분 및 밀도의 수직 분포를 조력발전소 가동전과 가동후 상황으로 구분하여 도시하였다(Fig. 22). 하계의 평균 수직구조는 호 상류로 갈수록 표층에서 고온, 저염의 특성을 나타내고 밀도의 수직구배가 큰 반면, 하류로 갈수록 수직적으로 균질한 형태를 나타낸다.

Fig. 22.

Averaged vertical profile of temperature, salinity and sigma-t(σ_t) at 6 stations(St. 5, 7, 9, 12, 13, 15) : circle symbol(●) represent without tidal power plant operation(Case N) and x symbol(x) represent with tidal power plant operation(Case A).

조력발전소 가동시 수온은 호 상류 일부구간에서 상승하지만 하류에서는 하강하는 것으로 나타나는 반면, 염분의 경우에는 호 전체적으로 증가하는 것으로 나타난다. 따라서 조력발전소 가동시 밀도는 수층에서 전체적으로 증가하며, 호 하류에서는 수직구배가 상대적으로 약화되는 것으로 나타났다.

(3) 밀도 수직단면 변화

조력발전소 가동시 해수의 유입양상을 살펴보기 위해 조력발전소를 통과하는 단면(Line 1)과 기존갑문을 통과하는 단면(Line 2)을 설정하여 조력발전소 가동으로 인한 수직 방향의 수온과 염분분포 및 밀도분포를 분석하였다(Fig. 23). 조력발전소를 통과하는 Line 1의 경우 조력발전소 가동시 외해수의 유입으로 인해 호 내측의 평균수온이 가동 전보다 소폭 하강하고 평균염분은 증가한다. 고염수가 저층으로 유입되는 2층 순환구조(two-layer circulation)가 상대적으로 뚜렷해지며, 이에 따라 평균밀도도 증가하는 경향을 보인다(Fig. 24). 기존 갑문을 통과하는 Line 2의 경우 Line 1과 마찬가지로 조력발전소 가동시 호 내측의 수온은 소폭 하강하는 반면 평균염분과 평균밀도는 증가하는 경향을 보인다(Fig. 25).

Fig. 23.

Location of lines for cross-sectional analysis.

Fig. 24.

Averaged vertical distribution of temperature, salinity and sigma-t(σ_t) along line 1 in summer: Left panels represent without tidal power plant operation and right panels represent with tidal power plant operation.

Fig. 25.

Averaged vertical distribution of temperature, salinity and sigma-t(σ_t) along line 2 in summer : Left panels represent without tidal power plant operation and right panels represent with tidal power plant operation.

(4) 구역별 수괴 특성

조력발전소 가동으로 인한 시화호 내측의 수괴특성을 파악하기 위해 시화호를 6개의 구역으로 구분하였고(Fig. 26), 각 구역에 대한 하계의 모델 격자별 평균수온과 평균염분을 산정하여 T-S Diagram을 표층과 저층으로 구분하여 작성하였다(Fig. 27).

Fig. 26.

A segmentation map.

Fig. 27.

Averaged T-S diagrams for water mass analysis in summer.

조력발전소 가동시 시화방조제에 인접한 Zone A와 Zone B 그리고 Zone C에서는 해수유통량의 증가로 인한 염분 상승으로 전체적으로 밀도가 증가하고, 시화호 중앙부인 Zone D와 Zone E에서는 조력발전소 가동 전·후의 밀도 변화가 크지는 않지만 배수시 상류에서 유입되는 담수의 영향에 의해 각 구역 내에서 밀도가 큰 폭으로 감소하는 구역이 나타났다. 한편 담수의 영향을 가장 크게 받는 상류 지역 Zone F에서는 조력발전소 가동 전·후의 표층 밀도분포는 비슷하지만 저층의 밀도분포는 조력발전소 배수시 수위하강으로 인해 단위부피당 담수의 영향이 저층까지 작용하여 밀도가 큰 폭으로 감소하는 것으로 나타났다.

(5) 구역별 성층 변화

조력발전소 가동으로 인한 시화호 내측의 6개 구역(Fig. 26)에 대한 밀도성층의 변화를 분석하였다. 밀도성층의 변화는 성층강도 계산을 통해 분석하였는데 여러 연구에 사용된 수직 적분한 potential energy anomaly(이하 PEA) 개념을 도입하여 하구 및 연안의 성층분석을 수행하였다(Choi et al., 1998; Lee and Cho, 1998; Simpson and Bowers, 1981; Simpson et al., 1990; Rippeth et al. 2001; Lie et al., 2002; Palacios-Hernandez et al., 2002). PEA는 수층을 완전 혼합시키는데 필요한 단위 체적당의 일인데 수층의 현재 위치에너지와 이 수층이 수직적으로 완전히 혼합되었을 경우의 위치에너지의 차이를 의미한다(Kim et al., 2006).

여기서, ρ(z)는 총 수심 h인 수층의 수직밀도분포이고, g는 중력가속도, 는 수심평균 밀도이다. PEA가 0이면 완전혼합된 것이고 성층이 강할수록 값이 커지는 반면, 수직적으로 밀도가 역전된 경우는 음의 값이 된다.

여름철(7월~9월)에 6개 구역에 대한 평균 PEA를 분석한 결과(Table 7), 조력발전소 가동 전에 시화호의 PEA는 구역별로 6.8~37.0J·m⁻³으로 시화호 상류로 갈수록 성층강도가 커지는 것으로 나타났다. 조력발전소 가동 시에는 PEA가 구역별로 1.5~20.1J·m⁻³으로 상류로 갈수록 성층강도가 커지는 경향은 유사하나 조력발전소 가동 전에 비해 성층이 약해지는 것으로 나타났다. 이는 조력발전소 가동으로 인해 해수의 유입·유출이 활발해지고 수위의 변화가 커짐에 따라 혼합효과가 증가하기 때문인 것으로 사료된다.

Table 7.

Spatially averaged potential energy anomaly(J·m^-3)