수질관리모델링 MIDTERM EXAM

1. 수질관리모델링MIDTERM EXAM 환경공학과 20051444 김태희

2. 1.study and explain the water quality model QUAL2E. 1)Theory (Physical, Chemical, and Biological Processes) Qual2e는 1차원 정상상태(steady state)와 1차원 가동적상태 (Dynamic state) 모의 가능한 1차원 수질예측모형으로서 하천구간 내 특정지점에서의 수질현상을 모의한다. 하천의 구간은 총 8개의 형태로 설정이 가능하며 모형의 매개변수는 수동보정을 주로 활용. (1) 수리계수 수리(유량)계수는 하천의 유량, 유속, 수심자료를 바탕으로 입력되 며, 실측자료 혹은 모형(ex: HEC-2)모의자료를 활용하나 일반적으 로 모의 결과를 활용한다.

3. 1.study and explain the water quality model QUAL2E. (2) 수질계수 BOD와 DO의 관계결정계수 조류와 영양염류 관계결정계수 일반적으로 실측자료를 활용하나 자료획득이 어려운 경우 매뉴얼 상의 값을 적용한 후 보정해 나가는 방법을 이용한다. 하천 시스템 내에서 물질이동 및 변환은 유체의 흐름에 의하여 농 도 부하가 이동되는 유속에 의한 유송과 유체흐름의 난류성과 물질 농도 차이에 의해 발생하는 확산, 화학 및 생물학적 반응에 의한 물 질의 증감, 시스템 내의 오염 부하량 유입 또는 유출에 의한 부하의 증감에 의해 좌우된다.이와 같은 과정에 대하여 물질이동방정식 을 세우면 다음 식(1.1)과 같다.

4. 1.study and explain the water quality model QUAL2E. (1.1) 여기서, C는 물질의 농도를 나타내고, A는 유로 단면적, D는 종확산계수, U는 평균유속, S는 물질의 유입부하량 혹은 유출량, V는 단위 구간의 수체적, t는 시간, X는 거리, R는 반 응에 의한 물질의 증감량을 나타낸다. QUAL2E는 에너지 보전의 원리를 이용하며 열전달을 고려한 에너 지 방정식을 이용하여 수온을 계산한다. 지배식은 식(1.2)와 같으 며 물질 이동 방정식과 흡사한 형태를 가진다. (1.2) 여기서, T는 수온이고, ρ는 물의 밀도, c는 비열이다

5. 1.study and explain the water quality model QUAL2E. 2) Numerical Algorithm QUAL2E 유한차분법 알고리즘 (1.1) 이 모형에서 사용한 수치해석 기법은 유한차분법으로서 요소간의 수질 농도를 선형적으로 가정하는 방법을 사용하고 있다.또한, 요 소 수질농도가 미지의 다른 요소 수질농도들의 함수로 표현되는 음 해법으로 차분 방법식이 세워지도록 하고 있다.식(1.1)의 물질이동 방정식을 1차반응을 가정한 선형 방정식에 대입하여 다시 표현하 면 식(1.2)와 같다. (1.2)

6. 1.study and explain the water quality model QUAL2E. 여기서, kC는 해석 수질항목의 반응에 의한 증감을 나타내고, P는 해석 수질항목의 반응을 제외한 반응을 나타낸다. 지배식(1.2)를 음해법을 이용한 유한차분식으로 표현하면 식(1.3) 과 같다. (1.3) n+1 시각의 수질항목 C에 대하여 정리하면 식(1.4)와 같다. (1.4)

7. 1.study and explain the water quality model QUAL2E. 여기서, δXi : 요소간의 길이 유한차분 방정식은 3개의 절점간의 관계로 표시된 삼각행렬로써 가우스 소거법에 의하여 해를 구하게 된다.

8. 1.study and explain the water quality model QUAL2E. 3) Run the example data

9. 1.study and explain the water quality model QUAL2E

10. 1.study and explain the water quality model QUAL2E

11. 1.study and explain the water quality model QUAL2E

12. 2. Study and explain the water quality model WASP6. 1) Theory (Physical, Chemical, and Biological Processes) WASP6(The Water Quality Analysis Simulation Program-6) 모형 -1981년EPA의 Di Toro 등에 의하여 개발. -사용자가 다양한 오염물 관리 계획을위해 자연 현상과 인위적인 오염원에 의한 수질 변동을 해석하고 예측할 수 있게 해줌. -WASP6 모형은 2개의 독립된 프로그램인 DYNHYD 모형(수체의 수리적 특성에 대한 모형)과 WASP5모형(수체에서 오염물질의 이 동과 반응을 예측)으로 구성. -부영양화와 BOD, DO 등 전통적인 오염원을 대상으로 하는 EUTRO 모형과 유기화합물, 금속 등의 독성 물질을 대상으로 하 는 TOXI 모형의부 프로그램으로 구성.

13. 2. Study and explain the water quality model WASP6. →화학적 생물학적 과정 -TOXI는 화학물질이 겪을 수 있는 대부분의 물리, 화학, 생물학적 반응을 다룸. -흡착, 이온화, 휘발, 가수분해, 광분해, 산화, 생분해와 여기서 언급되지 않은 여분의 반응도 사용자의 정의에 따라 다룰 수 있다. -가장 간단한 형태로 화학물질에 수체와 침전층 내에서의 이동과 화학물질의 전환, 고체로의 흡착을 적용. 이를 위해 1차 부패 상수(first-order decay constants)와 평형 분 배계수(equilibrium partition coefficients)를 사용할 수 있는데, 조 금 복잡하게는 2차 반응과 비선형 흡착 등온선이나 1차 흡착․탈착 속도 상수를 이용할 수도 있다.

14. 2. Study and explain the water quality model WASP6. 가. 전환 휘발, 생분해, 가수분해, 산화, 광분해 등의 반응에 대해 각각의 일차반응 속도상수를 입력하면, 총 반응은 다음 식에 의해 각 개별 반응의 합에 기초하게 된다. (1) 여기서, Kki = 반응 k에 대한 화학물질 ⅰ의 1차 전환 상수, day-1 생분해 -미생물에 의한 전환, 분해는 박테리아의 효소 시스템 안에서 화합물을 분 해하는것. -생분해는 유기 화학물질에 대한 생물체의 효소 공격이라는 넓고 복잡한 과정을포함. 박테리아와 적은수의 조류는 지표수의 생분해에서 중요한역할.

15. 2. Study and explain the water quality model WASP6. -독성 화학물질을 분해하는 박테리아 개체수의 성장 동역학은 잘 알려져 있지 않았음. 경쟁 기질과 그 유기물을 분해하는 다른 박테리아의 존재, 화 학물질에의 적응가능성 혹은 공대사는 개체수 변화를 정량화하기에 어렵. -결과적으로 독성 화학물질모델은 직접적으로 박테리아를 모델링하기보 다는 일정한 생물학적 활성을 가정. 때때로 다른 수계에서 측정된 1차 생 분해 속도 상수가 직접 쓰이기도 함. (2) 여기서, Ki = 화학물질 ⅰ의 segment j에서의 1차 생분해 속도 상수, day-1 TOXI에서는 1차 생분해 속도 상수 대신 수체 혹은 저니에서의 반감기를 이용 (3) 여기서, T1/2 = 반감기, days 만일 이 속도 상수들이 유사한 조건에서 측정되었다면 이 일차적 접근은 좀 더 복잡한 접근법과 비슷한 수준으로 정확할 것이다

16. 2. Study and explain the water quality model WASP6. 나. 평형 흡착 흡착은 일반적으로 다른 반응에 비해 빠르게 일어나므로 평형이 가 정. 적당한 농도(10-5M 미만이거나 용해도의 절반 미만)에서 흡착 평형은 용존 상태의화학물질의 농도에 비례. (4) 여기서, KPS = 분배 계수(partition coefficient) 평형상태에서 상간의 분포는 분배 계수 KPS에 의해 조절. 각 상에 서 화학물질의 총 질량은 KPS와 존재하는 고체상(DOC상 포함)의 양에 의해 조절. 식(4)는 순간적인 평형과 가역성이 가정될 때 사용 가능.

17. 2. Study and explain the water quality model WASP6. Water column과 benthic segments 안의 용존 화학물질은 고형물 입자와 용존유기 탄소와 반응해 용존, DOC에 흡착, 고형물에 흡착 (3가지 종류의 고형물)된 고체상의 5가지 상을 형성. 이 반응은 수 체의 단위 부피에 대하여 기술. (5) (6) 여기서, n은 공극률이다. 정반응은 흡착이고 역반응은 탈착. 이들 반응은 일반적으로 빠르고 국부평형으로 고려. Cw와 Cs, CB는 평형 분배 계수 Kps0와 KpB(L/kg)에 의해 조절.

18. 2. Study and explain the water quality model WASP6. (7) (8) 이 식은 고형물과 DOC에 흡착될 수 있는 공간이 충분하다는 가정 하에 선형 Freundlich 등온선으로 주어진다. (9) (10)

19. 2. Study and explain the water quality model WASP6. 화학물질의 총 농도는 다섯 가지 상 농도의 합이다. (11) 위 두 식을 인수분해하고 용존 부분 fD로 나타내면 다음과 같다. (12) 유사하게 고형물와 DOC에 흡착된 부분은 다음과 같이 나타낼 수 있다. (13) (14)

20. 2. Study and explain the water quality model WASP6. 이들 부분은 분배 계수, 내부적으로 계산된 공극률, 예측된 고형물 농도, 정해진 DOC의 농도로부터 모의 수행 동안 결정. 주어진 총농 도와 다섯 상의 분율, 용존, 흡착, 생물흡착 농도는 다음과 같이 결 정된다. (15) (16) (17) 이들 다섯 가지 농도는 mg/L 단위이고 각 상 내에서의 농도로 표현 될 수 있다. (18) (19) (20) 이들 농도의 단위는 각각 mg/Lw, mg/kgs, mg/kgB이다.

21. 2. Study and explain the water quality model WASP6. →물리학적 과정 고형물의 운송은 수계에서 잠재적으로 아주 중요. 과량의 고형물 은 수질에 직접 영향을 준다. 고형물의 이동은 화학물질의 이동과 반응에도 영향을 미친다. 많은 화학물질들은 고형물에 강하게 흡착 되어 침강(settling), 재부상(resuspension), 매몰(burial)된다. 흡착은 또한 화학물질의 이송(transfer)과 전환율 (transformationrate)에도 영향을 준다. 고형물의 이동속도와 농도 는 대부분의 독성물질 연구에서 추정되어야 한다. WASP6에서는 구간별로 각 고형물에 대한 물질수지를 계산하는데, 이는 benthic segment와 water segment 모두에 적용된다.

22. 2. Study and explain the water quality model WASP6. (1) Water column transport 수체 내에서 고형물과 고체상 화학물질은 하부의 water segment로 침강하고 bed segment 표면에 퇴적(deposit). WASP6에서는 운송 영역 3,4,5에서 속도와 표면적으로 침강, 매몰, 재부상 율을 표현. 입자의 운반속도는 고형물과 입자상 화학물질의 유량을 얻기 위해 단면적과 곱해진다. 침강속도는 부유 입자의 크기와 관련해 Stoke의 침강속도의 범위 에 있어야 한다. (21) 여기서, Vs = Stokes velocity, m/day g = 중력가속도, 981 cm/sec2 μ = 동점성도, 0.01 poise(g/cm3-sec) at 20℃ ρp, ρw = 입자와 물의 밀도, g/cm3 dp = 입자의 유효직경, mm

23. 2. Study and explain the water quality model WASP6. (2) Benthic Exchange 순 침전과 재부상은 고형물과 입자상 화학물질의 benthicexchange 를 유발한다. (22) 여기서, WBs = net sediment flux rate, g/day S = 고형물의 농도, g/m3 wD = 침전(deposition) 속도, m/day wR = 재부상(scour) 속도, m/day Aij = benthic segment의 표면적, m2 ⅰ = benthic segment j = water segment

24. 2. Study and explain the water quality model WASP6. 유체역학이나 수질 예측 프로그램의 기본적인 원리는 질량 보존이 다. 대상이되는 수체의 부피와 수질을 결정하는 물질들은 시간과 공간의 변화에 따라 물질수지식을 이용해 설명된다. 이때 수체 내 용존 물질에 대한 물질 수지식은 직접적인 경로를 통하거나 확산에 의해 발생하는 부하, 이송․확산에 의한 물질 운반, 물리․화학․생물학 적 전환을 통해 유입되거나 유출되는 모든 물질을 설명해야 한다. 그림 1와 같은 유체 내에서의 물질 수지식은 다음 식(23)과 같이 표현할 수있다.

25. 2. Study and explain the water quality model WASP6. (23) 여기서,C = 수체 내 오염물의 농도, mg/L 혹은 g/m3 t = 시간, days UX, UY, UZ = X, Y, Z 방향 이송 속도, m/day EX, EY, EZ = X, Y, Z 방향 확산 계수, m2/day SL = 직접․확산 부하량, g/m3-day SB = 경계 부하량, g/m3-day SK = total kinetic transformation rate; positive is source, negative is sink, g/m3-day

26. 2. Study and explain the water quality model WASP6. 식(23)은 수직, 수평 방향의 완전 혼합을 가정하여 1차원으로 유도 될 수 있다. (24) 여기서,A = 단면적, m2 위 식(24)는 크게 세 가지 항으로 나누어 볼 수 있는데, 각각 운송, 부하, 전환 항을 나타낸다. 거리 및 시간에 대하여 물질의 농도가 연속함수이므로, 1차원으로 가정된 물질수지식을 Taylor series로 전개하면, 다음과 같다. (25) (26)

27. 2. Study and explain the water quality model WASP6. 고차 항을 무시하고, 식(25)에서 식(26)를 빼서 정리하면, 다음과 같다. (27) 우변의 첫 번째 항과 두 번째 항을 동일한 방법으로 정리하면, (28) 이므로, (29) (30) (31) 이다.

28. 2. Study and explain the water quality model WASP6. 그림 2 수질해석을 위한 모식도 (여기서, x는 위치, j는 구획번호, L은 길이) 식(29) 및 식(31)를 그림2을 참조하여 식(24)에 대입하고 구획 j에 대하여정리하면, 다음과 같다. (32)

29. 2. Study and explain the water quality model WASP6. 양변에 Lj를 곱하면, 다음과 같다. (33) 여기서, Rj,j+1 및 Rj-1,j 는 확산계수를 포함하는 항으로 이며, Qj,j+1은 구획 j에서 구획 j+1로의 흐름을 나타낸다. Cj,j+1은 구획 j와 구획 j+1의 경계면 농도를 의미하며, 가증인수 ν(0∼1)를 도입하면, Cj,j+1 = νCj+1 + (1-ν)Cj로 표현되므로, 구획 j에 구획 ⅰ가 연결되어 있고, 경계면을 ij라 하여 물질수지방정식을 정리 하 면, 다음과 같다. (34)

30. 2. Study and explain the water quality model WASP6. 3) Run the example data

31. 2. Study and explain the water quality model WASP6.

32. 2. Study and explain the water quality model WASP6.

33. 3. Study and explain the hydraulic model DYNHYD5. 1) Theory (Physical, Chemical, and Biological Processes) WASP5는 수리모델인 DYNHYD5와 부영양화모델인 EUTRO5, 독성물질 모델인 TOXI5로 구성되어 있는 데 그 중 WASP5내의 수리 모형인 DYNHYD5는 Potomac 연안 모형이었던 DYNHYD2를 확장한 모형이다. 이들은 서로 독자적으로 사용이 가능하고 연합하여서도 사용 가능. DYNHYD5 모형은 변하기 쉬운 조수 감만의 주기, 바람 및 불안정한 물질 의 유입을 가상하는 간단한 유체역학 모형. DYNHYD는 수질 프로그램, WASP이 수질 내의 오염물질에 대한 운동 그 리고 상호 작용을 가상하는 동안 물의 운동을 가상하는 유체역학 프로그램. 이 경우에 DYNHYD5가 동수역학적 계산결과를 WASP5로 넘겨주게 되는 데, 다른 동수역학적 프로그램들에 의한 결과도 WASP5 모형과 연계 운영 될 수 있다.

34. 3. Study and explain the hydraulic model DYNHYD5. DYNHYD5 모형은 절점과 수로로서 구성된 계산격자망에 대하여 바람에 의한 응력, 마찰력, 중력의 영향을 고려한 1차원 연속방정 식, 운동방정식을 지배방정식으로 하며, 수로에서는 운동방정식을 이용하여 유속과 유량을 계산하고, 절점에서는 연속방정식을 이용 하여 수위를 계산한다. 수리모델링 결과인 유속 및 수심은 수질모델의 입력자료로 활용된 다. 이러한 유체역학 모형은 에너지와 질량에 대한 보전 법칙에서부터 시작된다. DYNHYD5모델을 사용하고자 하는 경우에는 각 소구간의 수심과 수두 높이 등의 소구간 특성을 정확하게 파악하는 것이 무엇보다도 중요하며, 정확한 자료가 입력된 경우에 모형의 실행이 가능하다.

35. 3. Study and explain the hydraulic model DYNHYD5. 2) Numerical Algorithm

36. 3. Study and explain the hydraulic model DYNHYD5.

37. 3. Study and explain the hydraulic model DYNHYD5.

38. 3. Study and explain the hydraulic model DYNHYD5.

39. 3. Study and explain the hydraulic model DYNHYD5. 연속방정식 The equation of continuity is given by: A = cross-sectional area, m2 Q = flow, m3/sec For rectangular channels of constant width B (refer to Figure 7): B = width, m • H= water surface elevation (head), m = rate of water surface elevational change with respect to time, m/sec • = rate of water volume change with respect to distance per • unit width, m/sec

40. 3. Study and explain the hydraulic model DYNHYD5. 3) Run the example data DYNHYD5와 EUTRO5을 이용하여 용담하류 금강의 수질 모의 실험과 그 결과 DYNHYD5입력자료

41. 3. Study and explain the hydraulic model DYNHYD5. EUTRO5입력자료

42. 3. Study and explain the hydraulic model DYNHYD5. 용담하류 금강의 수질 모의 실험 결과

43. 3. Study and explain the hydraulic model DYNHYD5. DYNHYD5에서 보여지는 값들

44. 3. Study and explain the hydraulic model DYNHYD5. Wasp과 DYNHYD5의 구조

45. 4. Study and explain the Watershed and Sewer Management model SWMM. • Theory (Physical, Chemical, and Biological Processes) SWMM 모형은 도시지역의 수량 및 수질을 모델링하기 위하여 개 발된 모형이다. 유역내 강우, 지표유출, 수질변화에 관계된 모든 부 분의 해석을 수행할수 있게 개발되었으며, Extran의 모듈을 사용하 여 Backwater, 배수관 망에서의 수리 유동상황을 역동적으로 해 석할 수 있다. 1969년에서 1971년에 걸쳐 미국 EPA에서 개발된 SWMM은 도시 지역의 지표유출 해석에 광범위하게 사용되었으며, 1994년에 4.3 Version으로 개량 되었다. 유역내의 우수관로, 합류식관로, 자연배 수로 등에서의 유동, 수위, 오염물 농도 등을 해석한다.

46. 4. Study and explain the Watershed and Sewer Management model SWMM. SWMM은 실행, 유출, 이동, Extran, 저류 및 처리의 5개의 부모형 으로 구성된다. 실행모형은 모델링 전체를 조절하는 기능을 수행한 다. 즉, 부모형 간 자료를 구성하고, 입출력을 지시한다. 유출모형 은 강우에 의해 유역내 에서 발생하는 지표유출과 오염도 변화를 추적하여, 수문곡선과 오염도 곡 선을 출력한다. 유출모형의 계산 결과를 다른 모형에서 이용하여 모델링을 수행한다. 이동모형은 우기 및 건기의 하수시스템으로의 침입수량을 계산 한다. Extran 모형은 배수관망에서의 유량 및 수심 등의 유동 상황을 추적 한다. 관로내의 유량을 해석하기 위하여 연속방정식과 운동방정식을 연 계하여 계산한다. 저류 및 처리 모형은 수리, 수질에 대한 저류 및 처 리장치의 효율을 평가한다.

47. 4. Study and explain the Watershed and Sewer Management model SWMM. →물리학적 과정 강우와 용설에 의해 유출이 발생, 비선형 저류방정식에 의해 지표면 유출을 묘, 침투량 산정에는 Horton 혹은 Green-Apmt방정식을 사용, 지표수 수질은 건기 동안 오염물질의 선형/비선형축적과 감소 인자에 의한 오염 물질의 유실을 기본으로 함. 건기시 유량 및 수질계산은 인구밀도와 같은 통계적인 인자와 Shield 기준에 의한 관로내 SS의 퇴적을 기본으 로 함

48. 4. Study and explain the Watershed and Sewer Management model SWMM. →화학적 과정 - Runoffbolck에서는 10개의수질, transportmodul에서의 3개의 수질모의가 가능 - 다른 모형과의 연계–수질모형인 WASP와 DYNHTD모형, 독립적 이거나 전체의SWMM모형과 HEC,STORM,QUAL-Ⅱ모형 - 모형의 정확성–유출은 약간의 보정에 의해 비교적 정확한 결과를 얻을 수 있 다. 수질은 폭 넓은 보정을 필요로 하며 실측 자료가 없을 시 매우 부정확함. →모형의 구조 Blosk 형태로 구성– 5개의 serviceblock, 4개의 computational block 각각의 block은 특정한 기능을 가지며 block 특정한 기능을 가지며 다른 block에서 사용.

49. 4. Study and explain the Watershed and Sewer Management model SWMM. 일반적으로 executive block과 하나 혹은 두개의computationnal block과 함께 실행

50. 4. Study and explain the Watershed and Sewer Management model SWMM. -RunoffblocK *본 모형은 swmm모델의 가장 중요한 모형 *강우 및 온도 또는 사용자 가정의 한우량 분포로부터 기상자료를 계산한 후 눈융해 모의실험기능과 함께 비선형 저류 접근법을 사용 하여 강우-유출과정을 모의실험 *최종적으로 유출모형은 하수관로 유입부에서의 수분곡선 및 오염 도 곡선을 출력한다. *본 모형은 몇분–몇년의 가간에 대한 모델링을 수행한다. *몇 주 미만의 모의실험은 단일사상방식, 장기모의 실험은 연속방 식으로 수행된다. 융설을 제외하고는, 모든 계산은 두 가지 경우에 대해 동일하게 수행 *유출부모형에서의 수질분석 과정은 옵션의 변화를 통해 지표 유출 성분 부하의 산출을 포함