• Journal of the Korean Association of Geographic Information Studies
• /
• v.11 no.3
• /
• pp.151-165
• /
• 2008

Ground slope and aspect are important parameters for physical deterministic water balance models like BROOK90 or hydrological models which attempt to calculate evapotranspiration, snowmelt, and net radiation. This study constructs a Digital Elevation Model(DEM) and examines how DEM resolution can change the average ground slope and aspect of a river basin and attempts to evaluate the effects on simulation results of BROOK90, a physical deterministic water balance model. The study area is Byungcheon river basin in Korea. DEM has been constructed using a 1:25,000 digital map with the methods of TIN and Topo To Raster. The total of 20 DEMs with 10m~100m resolution have been constructed, with a 10m interval. It was found that the higher the DEM resolution, the steeper the average ground slope value of the Byungcheon river basin. In turn, the direct solar radiation of a hilly area in the model increased the evapotranspiration and reduced the stream runoff in the Byungcheon river basin. On the other hand, a lower DEM resolution tends to move the average aspect from southeast to south in the Byungcheon river basin. Accordingly, it was found that stream runoff was reduced and evapotranspiration increased.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.469-469
• /
• 2023

고랭지밭은 집중호우시 토사가 유실되고 고농도의 흙탕물이 하천으로 유입되면서 하류지역의 수질 악화로 사회적·환경적 문제를 초래하고 있다. 이에, 환경부는 비점오염원관리지역을 지정하고 비점오염저감시설을 보급하고 있으나, 저감 노력에도 불구하고 집중호우시 수질 악화는 지속적으로 나타나고 있다. 따라서 발생원 관리 노력이 더욱 많이 필요한 실정인데, 최근 강원도와 환경부에서는 토사유실에 취약한 고랭지밭의 지형적 조건을 개선하여 토양유실을 저감하는 고랭지밭 경사도 완화 사업을 추진하고 있다. 고랭지밭 경사도 완화는 급경사지 및 경사장의 경작지를 계단식으로 조성하여 비점유출을 최소화시키는 발생원 관리방안이다. 그러나, 고랭지밭 경사도 완화에 따른 정량화된 비점유출 저감효과 분석에 관한 연구는 부족한 실정이다. 따라서 본 연구에서는 홍천군 광원리와 창촌리에 조성된 경사도 완화 경작지를 대상으로 강우유출수 모니터링을 수행하였으며, 비점저감효과 정량화를 위해 인근의 경사지 밭에서도 강우유출수 모니터링을 수행하여 대조구로 분석하였다. 강우유출수 모니터링은 2020년 8월에 총 3회(광원리 2회, 창촌리 1회) 수행되었으며, 분석 결과 광원리의 경사도 완화 경작지 평균 수질농도는 탁도 124.4 NTU, SS 111.5 mg/L, TOC 4.6 mg/L로 나타났으며, 대조구의 평균 수질농도는 탁도 1,741.7 NTU, SS 673.3 mg/L, TOC 30.6 mg/L로 나타났다. 광원리의 경사도 완화 경작지는 대조구 대비 수질항목별 83.4 ~ 92.9%의 비점저감 효과가 있는 것으로 나타났다. 창촌리의 경사도 완화 경작지의 평균 수질농도는 탁도 598.1 NTU, SS 414.5 mg/L, TOC 8.5 mg/L로 나타났으며, 대조구 평균 수질농도는 탁도3,487.3 NTU, SS 3,081.2 mg/L, TOC 40.3 mg/L로 나타났다. 창촌리의 경사도 완화 경작지는 대조구 대비 수질항목별 78.9 ~ 86.5%의 수질저감효과가 있는 것으로 나타났다. 경작지의 경사도 완화로 밭에서 발생하는 비점오염 발생을 SS는 최대 86.5%, TOC는 최대 84.9% 줄일 수 있는 것으로 나타났다. 하지만 본 연구 결과는 단년간 모니터링을 통해 도출된 결과이므로 정량화된 비점오염저감효과 도출을 위해서 다양한 강우조건 등을 고려한 지속적인 모니터링이 필요할 것으로 판단되며, 향후 연구에서는 2022년도에 인제군 북면 월학리 신규 조성된 경사도 완화 경작지를 대상으로 강우유출수 모니터링을 수행할 계획이다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.21-21
• /
• 2022

충적하천의 하상과 제방은 흐름에 의해 쉽게 침식되는 사립자들로 구성되어 있다. 흐름이 더욱 강해지면 커다란 입자들도 움직이게 되고 더욱 커지면 결국 모든 입자들이 움직인다. 이렇게 흐름에 의해 하상과 제방을 구성하는 사립자들이 움직이기 시작하는 상태를 한계운동이라 하며, 이때의 흐름조건을 한계조건이라 한다. 한계조건은 이동상 수리학의 시작으로, 이때부터 흐름은 물과 유사의 혼합이라는 이상류가 되어 순수한 물의 흐름보다 더 복잡하게 된다. 개수로는 관수로와 달리 반드시 자유수면이 있으며, 따라서 물과 공기와의 마찰은 상대적으로 매우 작으므로 개수로의 전단응력 분포는 관수로와 달리 근본적으로 비대칭이다. 따라서 전단응력은 수로 바닥이나 측벽에서만 작용하게 된다. 한계조건은 이러한 추상적인 의미에 앞서 바닥이 침식되지 않는 하도나 수로의 설계 등에 기본적인 자료가 된다. 개수로에서 경계면의 전단응력을 힘으로 표시하는 것을 통상 소류력이라 하며, 개수로 경계면에서 전단응력의 분포도; 이른바 단면의 평균 전단응력의 개념을 도입하여 해석하고 설계기준으로 제시되고 있다. 본 연구에서는 자연하천과 유사한 조건의 건설기술연구원 하천실험센터의 급경사수로에서 연구를 진행하였으며, 기존 연구를 바탕으로 제작된 소류력 측정장치를 이용하였다. 하천의 설계나 평가에 적용되는 평균 소류력 개념은 복잡한 난류흐름에서 평가지표로써 대표하기 힘들기 때문에 유사 하천환경의 바닥에서 발생하는 소류력을 직접 측정하고자 연구를 진행하였고, 연구결과를 기존의 소류력 산정방법과 비교하였다. 본 연구에 사용된 장치는 실규모 실험을 위해 제작되었으며, 실규모 적용성 검토를 위해 실험실에서 충분한 검증실험을 거친 후 실규모 실험에 적용하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2018.05a
• /
• pp.346-350
• /
• 2018

우리나라는 전 국토의 70%가 산지이고 하천경사가 다른 나라에 비해 상대적으로 급하여 홍수 관리에 매우 불리한 조건을 가지고 있으며, 특히 홍수기간의 집중호우 및 돌발홍수는 인명과 재산의 막대한 피해를 입히고 있다. 최근은 기후변화로 인하여 극심한 홍수, 가뭄 등 재해의 발생빈도가 증가하는 추세로 기후변화의 영향을 최소화할 수 있는 수재해 방재관리가 필요한 상황이다. 중 대하천의 경우에는 비교적 수재해 방재관리가 잘 이루어지고 있으나, 소하천(일부 중하천 포함)의 경우에는 취약한 구조를 보이고 있다. 특히 홍수기간(7월~9월)의 인명과 재산의 피해는 주로 소하천 위주로 발생하고 있으며, 사전 사후의 체계적인 대응이 이루어지지 못하고 있다. 수재해 방재관리를 위해서는 일차적으로 수문자료의 획득에 있으며, 그 이후 해당 유역에 적합한 수재해 대응을 위한 체계적인 방법론과 방재시스템 개발 운영이 수반되어야 안전한 방재관리를 할 수 있다. 따라서 수재해 방재관리 체계를 구축하기 위해서는 중 소규모 유역 단위를 대상으로 지속적이고 신뢰성 있는 자료의 획득과 축적이 중요하므로 중 소규모 유역 단위의 대표성 있는 시험유역의 운영은 매우 의미가 있다고 볼 수 있다. 본 논문에서는 한국건설기술연구원에서 운영하는 설마천 유역(유역면적 $8.48km^2$, 유로경사 2.15%, 경기도 파주시 적성면 소재)의 신뢰성 높은 2017년 관측자료를 이용하여 강우특성, 유출특성, 증발산량과 지하수 함양량 산정 등 수문특성을 분석하였으며, 과거 관측결과와 비교하였다. 강우특성분석으로는 호우사상 분리, 주요 호우사상 분석, 지속기간별 최대강우량, 시간분포 등이 있다. 2017년은 2016년보다 최대 강우지속기간은 크게, 평균 강우지속기간은 다소 작게 나타나며, 최대 강우강도와 평균 강우강도는 크게 나타나는 호우의 특징을 보인다. 이러한 호우특징은 연 유출률을 증가시킨 원인으로 보인다. 2017년의 하천유출률은 64.8%로 과거 2010년~2016년의 평균 유출률 68.4%와 차이를 보이지 않으나, 2016년의 57.4% 보다는 많이 증가된 유출률을 보인다. 2017년의 증발산량(328.8mm)과 지하수 함양량(24.8mm)은 2016년의 증발산량(458.5mm)과 지하수 함양량(32.8mm)보다 작게 나타나는 것으로 분석되었다. 이와 같이 산정된 수문자료는 수재해 방재를 위한 기초자료로 매우 유용하게 활용될 수 있으며, 방재관리를 위한 의사결정 과정에 중요한 역할을 할 수 있으므로 지속적인 시험유역의 운영은 매우 필요하다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.442-446
• /
• 2017

우리나라는 전 국토의 70%가 산지이고 하천경사가 다른 나라에 비해 상대적으로 급하여 홍수 관리에 매우 불리한 조건을 가지고 있으며, 특히 홍수기간의 집중호우 및 돌발홍수는 인명과 재산의 막대한 피해를 입히고 있다. 최근은 기후변화로 인하여 극심한 홍수, 가뭄 등 재해의 발생빈도가 증가하는 추세로 기후변화의 영향을 최소화할 수 있는 수재해 방재관리가 필요한 상황이다. 중 대하천의 경우에는 비교적 수재해 방재관리가 잘 이루어지고 있으나, 소하천(일부 중하천 포함)의 경우에는 취약한 구조를 보이고 있다. 특히 홍수기간(7월~9월)의 인명과 재산의 피해는 주로 소하천 위주로 발생하고 있으며, 사전 사후의 체계적인 대응이 이루어지지 못하고 있다. 수재해 방재관리를 위해서는 일차적으로 수문자료의 획득에 있으며, 그 이후 해당유역에 적합한 수재해 대응을 위한 체계적인 방법론과 방재시스템 개발 운영이 수반되어야 안전한 방재관리를 할 수 있다. 따라서 수재해 방재관리 체계를 구축하기 위해서는 중 소규모 유역 단위를 대상으로 지속적이고 신뢰성 있는 자료의 획득과 축적이 중요하므로 중 소규모 유역 단위의 대표성 있는 시험유역의 운영은 매우 의미가 있다고 볼 수 있다. 본 논문에서는 한국건설기술연구원에서 운영하는 설마천 유역(유역면적 $8.48km^2$, 유로경사 2.15%, 경기도 파주시 적성면 소재)의 신뢰성 높은 2016년 관측자료를 이용하여 강우특성, 유출특성, 증발산량과 지하수 함양량 산정 등 수문특성을 분석하였으며, 과거 관측결과와 비교하였다. 강우특성 분석으로는 호우사상 분리, 주요 호우사상 분석, 지속기간별 최대강우량, 시간분포 등이 있다. 2016년은 2015년보다 최대 강우지속기간과 평균 강우지속기간, 최대 강우강도와 평균 강우강도 모두 작게 나타나는 호우의 특징을 보이나, 50mm 이상의 호우사상수의 증가에 기인하여 연 총강우량 증가로 유출률이 증가하는 유출특성을 보였다. 2016년의 하천유출률은 57.4%로 과거 2010년~2015년의 평균 유출률 60.7%와 차이를 보이지 않으나, 2015년의 36.2% 보다는 매우 많이 증가된 유출률을 보이며, 2016년의 증발산량과 지하수 함양량은 2015년과 비슷한 값을 보이는 것으로 분석되었다. 이와 같이 산정된 수문자료는 수재해 방재를 위한 기초자료로 매우 유용하게 활용되며, 방재관리를 위한 의사결정 과정에 중요한 역할을 할 수 있기 때문에 지속적인 시험유역의 운영은 매우 필요하다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.886-886
• /
• 2012

우리나라 하천의 대부분은 산지에서 발원하며 전 국토의 약 70%가 산지하천 유역에 포함된다. 최근 기후변화로 인해 여름철 집중호우가 증가하고 있는 상황에서 강우의 예측이 어렵고 경사가 급한 산지하천 유역의 피해가 가중되고 있으며 산지하천의 강우를 정량적으로 파악하고 상시 모니터링 할 수 있는 체계의 구축이 요구된다. 한국건설기술연구원(2011)에서는 산지 하천유역 모니터링 시스템을 구축하여 재해위험지역의 현장관측시스템과 레이더강우를 기반으로 하는 강우유출 시스템을 연계운영하고 있다. 본 연구에서는 하천유역 모니터링 시스템을 통해 수집되고 있는 강원도 인제군 가리산리의 관측강우량을 이용해 산지하천유역의 강우특성을 분석하고 산지유역의 강우추정을 위한 레이더 자료의 활용성을 제시하였다. 대상유역인 가리산천 유역을 대상으로 작성된 티센 면적평균 강우량과 기상레이더를 이용한 레이더 강우량에서 가리산리 관측시스템 위치의 픽셀을 추출한 후 각각의 방법으로 추정된 강우량이 관측값과 어떤 차이를 갖고 있는가를 비교하였다. 또한, 모니터링 사이트 주위의 AWS를 이용해 레이더 강우를 보정한 후 동일한 방법으로 관측강우 위치의 셀 강우를 비교하여 레이더 강우의 보정 효과를 제시하였다. 연구결과 600m 이상의 고지대에 위치한 현장관측시스템의 강우는 고도가 낮은 인근 강우관측소와 큰 차이를 나타냈으며 티센면적 평균 강우의 경우 산지하천의 강우특성을 반영하기에 한계가 있는 것으로 판단되었다. 레이더 강우 역시 실제 관측강우량에 비해 과소추정되며 대상유역 주변의 AWS를 이용해 보정한 레이더 보정강우를 활용시 현장관측시스템의 강우가 가장 유사한 결과가 도출되었다. 본 연구를 통해 산지하천 유역의 강우특성을 파악하기 위해서는 지상관측소와 레이더 자료를 병행하여 활용하는 것이 필요하며 산지하천유역의 강우를 효과적으로 모니터링 하기 위해서는 고도에 따른 관측망의 구성이 필요할 것으로 판단되었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2011.05a
• /
• pp.329-329
• /
• 2011

제방은 하천의 범람으로부터 제내지의 인명, 가옥, 재산, 각종 시설 등을 보호하고 유수의 원활한 소통을 유지하기 위해 하천을 따라 축조한 구조물로서 침식 피해를 방지하기 위해 호안을 설치하도록 규정하고 있다. 호안은 만곡부, 하천 급경사, 지형의 간섭효과 등에 의한 유수의 침입으로부터 제방 및 하안을 직접 보호하기 위해 제방 앞비탈에 설치되는 구조물을 의미한다. 따라서 호안의 안정은 곧 제체의 안정과 직결되는 요소라 할 수 있다. 호안의 안정성은 역학적인 측면에서 외력과 저항력의 크기에 따라 평가할 수 있다. 국내 실무에서 적용되고 있는 호안설계 방법은 하천설계기준 해설(2009)에 제시되어 있으며, 경험과 이론의 양면을 고려하여 설계를 수행하도록 하고 있다. 이로 인해 호안 안정성에 대한 역학적인 검토 방법의 한계 때문에 과대 또는 과소 설계의 우려가 있다. 따라서 호안설계시 유속 및 소류력의 크기에 따른 정량적 평가기법이 요구되는 상황이다. 본 연구에서는 호안의 정량적인 설계기법 개발을 위하여 수리실험을 통해 유속증가로 인한 호안재료의 이탈현상을 재현하였다. 실험수로는 폭 2 m, 길이 25 m, 수로경사 1/300의 직선수로 형태이며 우안의 1:2 경사면에 다양한 크기의 사석과 블록을 차례로 설치하여 실험을 수행하였다. 수리실험은 다른 하상변동 요인이 제한된 고정상으로 수행되었으며 정상류 흐름조건에서 유량을 변화시키며 유속변화 등 흐름현상에 의한 호안재료 이탈을 관찰하였다. 실험결과를 바탕으로 호안설계시 1차원 접근유속을 대표유속으로 적용하는 방법의 특성을 파악하고, 호안재료의 이탈유속과 흐름특성간의 상관관계를 분석하였다. 또한, 호안설계시 입력자료로 사용될 수 있는 단면평균유속과 2차원적 국부유속과의 상관관계를 파악하여 호안의 정량적 설계기법에 활용될 제원결정 경험식을 제안하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.437-441
• /
• 2017

우리나라는 전 국토의 70%가 산지이고 하천경사가 다른 나라에 비해 상대적으로 급하여 홍수 관리에 매우 불리한 조건을 가지고 있으며, 특히 홍수기간의 집중호우 및 돌발홍수는 인명과 재산의 막대한 피해를 입히고 있다. 최근은 기후변화로 인하여 극심한 홍수, 가뭄 등 재해의 발생빈도가 증가하는 추세로 기후변화의 영향을 최소화할 수 있는 수재해 방재관리가 필요한 상황이다. 중 대하천의 경우에는 비교적 수재해 방재관리가 잘 이루어지고 있으나, 소하천(일부 중하천 포함)의 경우에는 취약한 구조를 보이고 있다. 특히 홍수기간(7월~9월)의 인명과 재산의 피해는 주로 소하천 위주로 발생하고 있으며, 사전 사후의 체계적인 대응이 이루어지지 못하고 있다. 수재해 방재관리를 위해서는 일차적으로 수문자료의 획득에 있으며, 그 이후 해당유역에 적합한 수재해 대응을 위한 체계적인 방법론과 방재시스템 개발 운영이 수반되어야 안전한 방재관리를 할 수 있다. 따라서 수재해 방재관리 체계를 구축하기 위해서는 중 소규모 유역 단위를 대상으로 지속적이고 신뢰성 있는 자료의 획득과 축적이 중요하므로 중 소규모 유역 단위의 대표성 있는 시험유역의 운영은 매우 의미가 있다고 볼 수 있다. 본 논문에서는 한국건설기술연구원에서 운영하는 차탄천 유역(유역면적 $190.64km^2$, 유로경사 0.96%, 경기도 연천군 소재)의 신뢰성 높은 2016년 관측자료를 이용하여 강우특성, 유출특성, 증발산량 등 수문특성을 분석하였으며, 과거 관측결과와 비교하였다. 강우특성 분석으로는 호우사상 분리, 주요 호우사상 분석, 지속기간별 최대강우량, 시간분포 등이 있다. 2016년은 2015년보다 최대 강우지속기간은 적게, 평균 강우지속기간은 크게, 최대 강우강도와 평균 강우강도는 모두 적게 나타나는 호우의 특징을 보이나, 설마천 유역과 같이 50mm 이상의 호우사상수의 증가에 기인하여 연 총강우량 증가로 유출률이 증가하는 유출특성을 보인다. 2016년의 하천유출률은 강우량 대비 39.6%(장진교, 유역출구)와 58.9%(보막교, 중간 소유역)로 과거 2012년~2015년의 평균 유출률 48.9%와는 다소 차이를 보인다. 2015년의 장진교 38.2%는 차이를 보이지 않으나, 보막교의 38.3%는 많은 차이를 보인다. 강우-유출특성 분석결과 강우량의 증가에 기인하여 2016년의 연간 하천유출량은 2015년보다 크게 증가 하였다. 그리고 2016년의 증발산량은 강우량 대비 장진교 50.1%, 보막교 35.4%로 2015년의 장진교 49.9%와 비슷한 값을 보이나, 2015년의 보막교는 54.9%로 약 19.5% 정도의 적은값을 보인다. 온도, 습도, 풍속, 일조시간에 영향을 받는 증발산량은 2015년보다 일평균 습도는 변화가 거의 없었으나, 일평균 풍속의 증가, 일평균 기온과 일조시간의 감소에 기인하여 적은 증발산량을 보이는 것으로 분석되었다. 이와 같이 산정된 수문자료는 수재해 방재를 위한 기초자료로 매우 유용하게 활용되므로 지속적인 시험유역의 운영은 매우 필요하다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2020.06a
• /
• pp.256-256
• /
• 2020

하천의 유출률은 하천의 경사, 연장, 강우특성 및 규모 등에 따라 크게 영향을 받으므로, 제주도의 지역별 고도별 토양 및 강우 특성과 독특한 지형에 따른 간헐하천의 유출특성을 고려할 때 관측에 의해 유출특성을 파악하는 것은 현실적으로 어렵다. 특히나 홍수시 관측되는 유량의 불확실성과 지역별·고도별 강우량의 편차를 고려한다면 강수대비 정확한 관측유출률의 산정 자체가 쉽지 않다. 이는 기존 관측된 유출률이 하천유역별로 대상기간별로 차이가 크게 나타나는 것을 통해서도 알 수 있다. 따라서 특정 지점에서의 몇 개의 강우사상에 대한 유출률이 아닌 다양한 사상과 지역적 특성을 고려한 유출특성을 파악하기 위해서는 모델링에 의한 접근법이 필요하다. 본 연구에서는 제주도 하천유역에서의 건천화 현상과 간헐적 유출특성을 모의하는데 최적화된 제주형 SWAT-K 모형을 기반으로 제주도 11개 하천유역(유역면적 30 ㎢ 이상)에 대한 유출특성을 분석하였다. 유출률을 중심으로 계절별, 지역별, 하천유역별 차이를 비교 분석함으로써 강우규모에 따른 유출 특성과 시·공간적 유출률 변화를 평가하였다. 1992~2013년 기간에 대해 분석한 결과 제주도 전역에 대한 연간 유출률은 13.3%~30.5% (평균 21.8%)로 나타났다. 하천유역별로는 천미천(31.6%), 가시천(31.1%), 화북천(29.3%), 창고천(27.4%), 의귀천(27.1%)의 유출률이 상대적으로 크게 나타났으며, 효돈천(18.9%)와 도순천(19.0%)은 상대적으로 낮은 것으로 분석되었다. 11개 하천유역에 대한 월별 평균 유출률은 강수량이 많은 8월에 가장 높은 유출률(26.1%)을 보였다.

• The Korean Journal of Ecology
• /
• v.22 no.3
• /
• pp.101-108
• /
• 1999

This survey was performed from March 1993 to March 1998, in order to clarify the relationships between water quality and topographical factor. The study sites were two reservoir basins; Kaesim and Jangchan in Iwon-myon, Okchon-gun, Chungcho'ngbukdo Province. Basin shape factors of Kaesim reservoir were at 0.030∼0.210 (mean value 0.090), those of Jangchan reservoir were at 0.217∼0.452 (mean value 0.325). The mean basin shape factor of Jangchan reservoir was 3.61 times larger than that of Kaesim reservoir because its stream width was narrower and mean stream length was shorter. In the correlation between distance from the source of stream (L) and basin area (A), Iwonchon basin was calculated as L=1.44A/sup 0.6/. Circularity ratio was 17.114 in Kaesim (22% of Kum River), and 7.444 in Jangchan. Elongation ratio was 0.357 in Kaesim, 0.636 in Jangchan and 0.282 in Kum River. Precipitation summation period of Jangchan was 1.54 times slower than that of Kaesim. Rainfall reaching time in each small basin was 337.53 min. in A'(Jangchan-ri) basin of Jangchan and 49.26 min in H (Iwon-ri) basin of Kaesim. In the relationship between watershed frequency (Df) and drainage density (Dd), the regression equation was Df=0.023Dd² in Kaesim and Df=0.189Dd² in Jangchan reservoir. As slope degree increased, DO became higher (Y/sub DO/=0.19X＋6.5927, r=0.8l), but COD(Y/sub COD/=-0.2092X＋9.7104, r=0.52) became lower. Total nitrogen was increased with the increase of basin shape factor and circularity ratio. Ratio of B/sub OD/ to COD was 1/1.2(Y/sub BOD/ = 1.2984 X/sub COD/-3.2004, r=0.9l).