지점빈도해석과 지역빈도해석

왜 배우는가: 홍수량 빈도해석은 실측 홍수 자료로 설계홍수량을 정하는 가장 직접적인 방법이며, 자료가 충분한지에 따라 지점해석과 지역해석 중 무엇을 쓸지가 달라진다. 두 방법의 개념과 차이를 알아야 수문 설계와 기술사 답안의 절차 문제를 정확히 풀 수 있다. 시험에서는 절차와 적용 조건을 묻는 형태로 반복 출제된다.

출제 이력: 제138회 2교시 2번, 제125회 3교시 6번

1. 한 줄 정의와 핵심 키워드

지점빈도해석은 대상 지점의 연최대홍수량 자료를 확률분포에 맞춰 재현기간별 설계홍수량을 구하는 방법이고, 지역빈도해석은 동질한 여러 관측소 자료를 통합해 설계홍수량을 추정하는 방법이다.

핵심 키워드: 지점빈도해석, 지역빈도해석, 연최대홍수량(AMS), 재현기간, 초과확률, 동질지역, 지수홍수량, 지역 성장곡선

2. 개념 풀이

쉽게 말하면, 지점빈도해석은 "이 지점의 과거 홍수 기록으로 얼마나 큰 홍수가 몇 년에 한 번 오는지 계산하는 것"이다. 지역빈도해석은 "기록이 짧은 지점은 비슷한 주변 유역 기록을 빌려서 계산하는 것"이다.

정식 정의로는, 지점빈도해석은 지점빈도해석(at-site frequency analysis)으로 한 지점의 연최대홍수량 계열에 확률분포를 적합하여 재현기간 T에 대응하는 홍수량을 구하는 방법이다. 재현기간은 같은 크기 이상의 홍수가 평균적으로 한 번 발생하는 간격이다.

재현기간(return period)은 T = \dfrac{1}{P} 로 정의하며, P는 1년 동안 초과되는 확률(P = 초과확률)이다. 예를 들어 T=100년이면 초과확률은 P=0.01, 즉 매년 1%다. 재현기간과 초과확률은 서로 역수 관계이므로 답안에서는 두 값을 함께 쓴다.

지역빈도해석(regional frequency analysis)은 동질성이 확인된 관측소들의 자료를 지수홍수량(index flood)으로 무차원화해 통합하고, 그 결과로 지역 성장곡선(regional growth curve) G_T를 구하는 방법이다. 대상 지점의 설계홍수량은 지수홍수량 I와 성장계수 G_T의 곱으로 구한다.

$Q_T = I \times G_T$

여기서 Q_T는 재현기간 T의 설계홍수량(m³/s), I는 대상 지점의 지수홍수량(m³/s), G_T는 재현기간 T의 성장계수(무차원)다.

동질지역(homogeneous region)은 기후·지형·유출 특성이 비슷해서 같은 확률분포를 따른다고 볼 수 있는 관측소 집합이다. 지역빈도해석의 성패는 이 동질성 판단에 달려 있다.

3. 등장 배경과 필요성

홍수량 자료는 설계의 출발점이지만 관측 기간이 짧은 경우가 많다. 우리나라는 실측 홍수량 자료가 적다는 점이 자주 지적되며, 이 때문에 빈도해석의 신뢰도가 떨어질 수 있다. 재현기간이 긴 설계홍수량을 짧은 기록으로 외삽하면 표본오차가 커지므로, 여러 지점의 정보를 모으는 방법이 필요해졌다.

지역빈도해석은 이 문제를 해결하려고 발전했다. 자료가 없는 미계측 유역에서도 유사한 관측소의 정보를 빌려 설계홍수량을 추정할 수 있다. 미국 지질조사소(USGS)가 미계측 유역의 홍수 분위수를 예측하려고 지역 빈도해석 방식을 개발한 것이 대표적인 배경이다.

또한 최근에는 토지이용 변화와 기후변화로 홍수 자료의 정상성(stationarity) 가정이 흔들리고 있다. 그래서 비정상 빈도해석과 지역해석의 결합이 새로운 쟁점으로 떠오르고 있으며, 두 방법의 기본 절차를 정확히 이해하는 것이 그 기반이 된다.

4. 구성요소·원리·절차

지점빈도해석과 지역빈도해석은 같은 확률론적 틀을 쓰지만 자료를 누구에게서 얻는가가 다르다. 지점빈도해석(At-site Frequency Analysis)은 대상 지점의 연최대홍수량 자료만 쓰고, 지역빈도해석(Regional Frequency Analysis)은 동질지역 안의 여러 관측소 자료를 통합한다. 두 방법의 구성요소를 비교하면 다음과 같다.

표를 보면 두 방법의 차이는 분포를 한 지점에서 맞추는가, 여러 지점을 묶어 공통 형태(성장곡선)를 쓰는가에 있다. 지역해석은 자료 부족을 보완하지만 동질성 가정이 깨지면 결과가 크게 틀어진다.

원리는 다음과 같다. 지점해석에서는 표본 계열이 독립이고 동일한 분포를 따른다고 보고, 재현기간 T의 설계홍수량 Q_T를 분포의 분위수로 구한다. 재현기간과 초과확률의 관계는 다음과 같다.

$P(Q \ge Q_T) = \frac{1}{T}$

여기서 T는 재현기간(년), Q_T는 재현기간 T의 설계홍수량(m³/s)이다. 지역해석에서는 지수홍수량 I_i(예: 평균 연최대홍수량, m³/s)와 무차원 성장곡선 G_T를 곱해 대상 지점의 설계홍수량을 구한다.

$Q_{T,i} = I_i \times G_T$

여기서 Q_{T,i}는 지점 i의 재현기간 T 설계홍수량(m³/s), I_i는 지점 i의 지수홍수량(m³/s), G_T는 재현기간 T의 성장계수(무차원)이다. 두 식 모두 단위가 m³/s로 유지되는지 확인하면 된다.

지점해석은 자료의 개별성을, 지역해석은 자료의 공통성을 활용한다

지점빈도해석 절차

• 연최대홍수량 계열을 정리하고 이상치, 독립성, 정상성을 점검한다.
• 확률분포(GEV, Gumbel 등)를 선정하고 모멘트법, 최우도법, L-모멘트법으로 매개변수를 추정한다.
• 확률도와 적합도 검정으로 분포의 적합성을 판단한다.
• 재현기간 T에 해당하는 설계홍수량을 산정한다.

지역빈도해석 절차

• 유역특성과 유출특성이 유사한 관측소를 모아 동질지역을 설정한다.
• 불일치도(discordancy)와 이질성 검정으로 지역 내 지점의 동질성을 확인한다.
• 각 지점의 지수홍수량으로 관측값을 무차원화하고 통합한다.
• L-모멘트 등으로 지역 성장곡선 G_T를 추정한다.
• 대상 지점의 I_i를 산정하고 Q_{T,i} = I_i \times G_T로 설계홍수량을 구한다.

{"type":"flow","title":"지역빈도해석 절차","nodes":[{"id":"a","label":"관측소 자료","sub":"연최대홍수량 계열"},{"id":"b","label":"동질지역 설정","sub":"유역·유출 특성 유사"},{"id":"c","label":"동질성 검정","sub":"불일치도 검토"},{"id":"d","label":"지수홍수량 표준화","sub":"무차원 관측값"},{"id":"e","label":"지역 성장곡선 G_T","sub":"L-모멘트 추정"},{"id":"f","label":"설계홍수량 산정","sub":"I_i × G_T"}],"edges":[{"from":"a","to":"b","label":""},{"from":"b","to":"c","label":""},{"from":"c","to":"b","label":"재설정"},{"from":"c","to":"d","label":""},{"from":"d","to":"e","label":""},{"from":"e","to":"f","label":""}]}

5. 관련 법령·기준·표준

홍수량 빈도해석의 기준은 하천설계기준(KDS 51 00 00)의 설계수문량 항목에 있으며, 그 아래 홍수량 빈도해석 절에서 지점·지역해석의 적용 절차를 정한다. 설계홍수량의 명칭 체계는 기본홍수량과 계획홍수량으로 나뉜다. 기본홍수량은 댐과 유역변경 등 인위적 개발을 무시한 자연상태 홍수량의 확률홍수량이다. 계획홍수량은 개발을 고려해 하도가 부담하는 홍수량이다. 기준값과 적용 조건은 다음 표와 같다.

최근에는 기후변화와 토지이용 변화로 인한 비정상성(Nonstationarity) 검토가 요구되는 추세이며, 비정상성을 반영한 지역빈도해석 연구가 진행되고 있다.