臺灣博碩士論文加值系統

為了能夠瞭解滯洪池在夏天易於散發惡臭的問題，以及冬天寒流來襲所造成養殖業池塘的寒害，本文乃探討水溫、溶氧(DO)濃度與生化需氧量(BOD)濃度在水中的傳遞特性以及其相互關聯性。
本研究針對上述的問題進行以下的探討：
1. 單層夏天和冬天的水溫、DO濃度與BOD濃度的分析
2. 雙層夏天和冬天的水溫、DO濃度與BOD濃度的分析
3. 夏天滯洪池增加水深，水質變化的情形
4. 冬天寒流來襲時，增加水體深度，藉此瞭解養殖池水溫變化的情形
商用軟體其實已開發，可研究相關的問題，例如wasp，但是其使用的方式是以代數的方式進行運算，但擴散方程式為一微分方程式，在進行傳輸時，會有exponential function的邊界層型態，因此需更精密的運算。
本研究僅建立一個計算架構，對於更多細節都可以在此架構之下應用求得。
本研究已求得單層水溫在任何時間下的解析解與雙層水溫在時間無窮大之下的解析解，可供數值解作為驗證。數值方法乃修改Kuo et al., 2008所建立的「垂直分層延散計算法」。
根據本研究結果共獲得四個主要結論：
1. 在夏天DO濃度的變化主要來自於BOD濃度與溫度的影響，但是在冬天DO濃度的變化主要來自於初始濃度與BOD濃度。
2. 本研究另將垂向水體虛擬分為兩層，以模擬DO與BOD在一維雙層的非均質傳輸現象。即在本案例中，設定下半層熱傳導係數與延散係數遠低於上半層，模擬結果發現，以此方法模擬DO與BOD濃度，對於DO與BOD濃度的變化是沒有影響的，所以可以推知在此條件下，主要影響DO濃度與BOD濃度的是控制方程式的源項。
3. 夏天滯洪池溶氧濃度若降低至臨界值或臨界比例時，則可以增加水深的方式，補充溶氧以及增加溶氧變好的速率。
4. 模擬冬天寒流來襲時，若能增加水深，則可以明顯減少水溫降低速率。

In order to deal with the environmental problems including the stench of retarding basin in summer and the damage of aquaculture basin caused by the cold current, we analysis the relationship between water temperature, BO (dissolved oxygen), BOD (biological oxygen demand) and vertical delaminate. Besides the diffusion effects of these physic parameters in the water were analyzed in this thesis. Based on upon statement we will produce the following analysis: 1.Regarding the water in retarding basin as one layer and analyzing the water temperature, the concentration of DO and BOD in summer and winter. 2. Regarding the water in retarding basin as two layers and analyzing the water temperature, the concentration of DO and BOD in summer and winter. 3. Observing the variation of water quality as pouring water into the retarding basin in summer. 4. Realizing the variation of water temperature in the aquaculture basin caused by pouring water during cold current. In this thesis, we have already obtained the analytical solutions of the water temperature in one layer at any time and in two layers at the infinite time. The analytical solution was put to proof the simulation solution which was modified from “computation method for dispersion with vertical layers” built by Kuo et al. at 2008. We concluded four primary conclusions. 1. The variation of the DO concentration was mainly affected by the BOD concentration and the water temperature in summer. However, the variation of the DO concentration was mainly affected by initial concentration and the BOD concentration in winter. 2. In order to simulate the 1D heterogeneous diffusion effect of the DO and BOD concentration, we assumed the water layer as two parts in vertical direction. Besides, the heat diffusion and convection coefficients in layer below are far lower than in layer upon. Based on the simulation results, we found the DO and BOD concentration will not be affected by upon assumption. The source term of governing equation is the primary parameter which affected the DO and BOD concentration in this assumption. 3. In summer, if the dissolved oxygen of the retarding basin decreased to the critical value or ratio, we could pour water into the basin to improve the amount of the dissolved oxygen. 4. During the cold current in winter, if we pour water into the aquaculture basin, the amount of decreasing velocity of water temperature will be reduced obviously.

目錄
中文摘要 I
目錄 IV
表目錄 VI
圖目錄 VII
符號說明： IX
第一章 緒論 1
1.1 研究動機與目的 1
1.2 文獻回顧 2
1.3 研究方法 3
1.4 章節介紹 3
第二章 水溫、DO濃度、BOD濃度之數學模型 5
2.1 水溫的數學模型建立 5
2.1.1 太陽輻射引致熱傳輸的熱方程式 5
2.1.2 水體熱能量傳遞的邊界方程式 6
2.1.3 水體熱能量傳遞的初始方程式 7
2.1.4 一維單層熱擴散方程式 7
2.1.5 一維雙層熱擴散方程式： 7
2.2 DO質量傳輸方程式 8
2.2.1 一維垂向溶氧質量傳輸方程式 8
2.2.2 溶氧質量傳輸的邊界方程式 8
2.2.3 DO質量傳輸方程式初始方程式 9
2.2.4 一維單層DO質量傳輸方程式 10
2.2.5 一維雙層DO質量傳輸方程式 10
2.3 BOD質量傳輸方程式 11
2.3.1 一維垂向生化需氧量質量傳輸方程式 11
2.3.2 生化需氧量質量傳輸的邊界方程式 12
2.3.3 BOD質量傳輸方程式初始方程式 12
2.3.4 一維單層BOD質量傳輸方程式 12
2.3.5 一維雙層BOD質量傳輸方程式： 13
第三章 解析解求解計算過程 15
3.1 一維單層熱擴散方程式求解過程 15
3.2 一維雙層熱擴散方程式求解過程 17
第四章 數值模式 22
4.1 數值方法： 22
4.2 模式之計算步驟： 24
4.3 數值解的驗證 27
第五章 結果分析 31
5.1 一維單層夏天的案例 31
5.2 一維單層冬天的案例 38
5.3 一維雙層夏天的案例 44
5.4 一維雙層冬天的案例 48
5.5 夏天滯洪池(以增加水深解決問題)的案例 53
5.6 冬天寒流(以增加水深解決問題)的案例 58
第六章 結論與建議 64
參考文獻 66
附錄A 水溫參數探討 68
附錄B 數學模型所遇到的困難 73
附錄C　 與 的求解過程 75
附錄D 與 的求解過程 80
附錄E 數值方法所遇到的困難 84
附錄F 垂向分層延散計算法計算過程 92


表目錄
表4. 1 本研究的數學模型對照表 22
表4. 2 Kuo et al., 2008建立的數學模型 23
表4. 3 水溫驗證所使用的參數表 27
表5. 1 本研究案例參數表 31
表5. 2 無因次時間與有因次時間的對照表 63


圖目錄
圖2. 1 一維單層熱擴散數學模型示意圖 7
圖2. 2 一維雙層熱擴散數學模型示意圖 8
圖2. 3 一維單層DO質量傳輸數學模型示意圖 10
圖2. 4 一維雙層DO質量傳輸數學模型示意圖 11
圖2. 5 一維單層BOD質量傳輸數學模型示意圖 13
圖2. 6 一維雙層BOD質量傳輸數學模型示意圖 14
圖4. 1 數值流程圖 26
圖4. 2 一維單層水溫的解析解與數值解驗證圖(t=1年) 28
圖4. 3 一維單層水溫的解析解與數值解驗證圖(t=2年) 28
圖4. 4 一維雙層時間無窮大之垂向數值水溫與解析解驗證圖 29
圖4. 5 一維單層水溫解析解與實測值驗證圖 30
圖5. 1 一維單層夏天垂向無因次溫度隨時間變化圖 33
圖5. 2 一維單層夏天垂向無因次DO濃度隨時間變化圖 33
圖5. 3 一維單層夏天垂向無因次BOD濃度隨時間變化圖 33
圖5. 4 一維單層夏天無因次垂向DO與BOD濃度比較圖 35
圖5. 5 一維單層夏天無因次垂向DO濃度與溫度關係圖 36
圖5. 6 一維單層夏天無因次垂向BOD濃度與溫度關係圖 38
圖5. 7 一維單層冬天無因次垂向溫度隨時間變化圖 39
圖5. 8 一維單層冬天無因次垂向DO濃度隨時間變化圖 40
圖5. 9 一維單層冬天無因次垂向BOD濃度隨時間變化圖 40
圖5. 10 一維單層冬天無因次垂向DO與BOD濃度關係圖 41
圖5. 11　一維單層冬天無因次垂向DO濃度與溫度關係圖 43
圖5. 12　一維單層冬天無因次垂向BOD濃度與溫度關係圖 44
圖5. 13 一維夏天無因次垂向溫度單層與雙層比較圖 46
圖5. 14 一維夏天無因次垂向DO濃度單層與雙層比較圖 47
圖5. 15 一維夏天無因次垂向BOD濃度單層與雙層比較圖 48
圖5. 16 一維冬天無因次垂向溫度單層與雙層比較圖 50
圖5. 17 一維冬天無因次垂向DO濃度單層與雙層比較圖 51
圖5. 18 一維冬天無因次垂向BOD濃度單層與雙層比較圖 52
圖5. 19 夏天滯洪池無因次垂向溫度變化量與不同水深的比較圖 55
圖5. 20 夏天滯洪池無因次垂向DO濃度變化量與不同水深的比較圖 56
圖5. 21 夏天滯洪池無因次垂向BOD濃度變化量與不同水深的比較圖 58
圖5. 22 冬天寒流來襲無因次垂向溫度與不同水深的比較圖 60
圖5. 23 冬天寒流來襲無因次垂向DO濃度與不同水深的比較圖 61
圖5. 24 冬天寒流來襲無因次垂向BOD濃度與不同水深的比較圖 63
圖A. 1 Secchi disc的照片 68
圖A. 2 參數 =0.03時一維垂向無因次溫度隨時間的變化圖 68
圖A. 3 參數 =10時一維垂向無因次溫度隨時間的變化圖 69
圖A. 4 一維垂向無因次溫度與不同參數 的關係比較圖 70
圖A. 5 參數 =25時一維垂向無因次溫度隨時間的變化圖 71
圖A. 6 參數 =125時一維垂向無因次溫度隨時間的變化圖 71
圖A. 7 一維垂向無因次溫度與不同參數 的關係比較圖 72
圖B. 1 經過一年不同下邊界條件的關係比較圖 74
圖B. 2 經過二年不同下邊界條件的關係比較圖 74
圖E. 1 一維垂向無因次DO濃度隨時間關係圖(法一, 50個格點) 85
圖E. 2 一維垂向無因次DO濃度隨時間關係圖(法一, 100個格點) 85
圖E. 3 一維垂向無因次DO濃度隨時間關係圖(法一, 200個格點) 86
圖E. 4 一維垂向無因次DO濃度不同網格數的關係圖一(法一) 86
圖E. 5 一維垂向無因次DO濃度不同網格數的關係圖二(法一) 87
圖E. 6 一維垂向無因次DO濃度不同網格數的關係圖三(法一) 87
圖E. 7 一維垂向無因次DO濃度隨時間關係圖(法三, 50個格點) 88
圖E. 8 一維垂向無因次DO濃度隨時間關係圖(法三, 100個格點) 89
圖E. 9 一維垂向無因次DO濃度隨時間關係圖(法三, 200個格點) 89
圖E. 10 一維垂向無因次DO濃度不同網格數的關係圖一(法三) 90
圖E. 11 一維垂向無因次DO濃度不同網格數的關係圖二(法三) 90
圖E. 12 一維垂向無因次DO濃度不同網格數的關係圖三(法三) 91

1.Bowie, G.L., Mills, W.B., Porcella, D.B., 1985., Rates, constants and kinetics formulations in surface water quality modeling. EPA/600/3-85/040.
2.Brown, L.C. and Barnwell, Jr., T.O., 1987. The enhanced stream water quality models QUAL2E and QUAL2E-UNCAS: Documentation and user manual. EPA-600/3-87/007, U.S. Environmental Protection Agency, Athens, GA.
3.Chapman, A.J., 1984. Heat transfer. 4th. Taipei Publications Trading company.
4.Chapra, S.C., 1997. Surface Water-Quality Modeling. McGraw-Hill Inc, New York.
5.Dake, J.M.K., Harleman, D.R.F., 1969. Thermal stratification in lake- analytical and laboratory studies. Water Resour. Res. 5, 484-495.
6.Farlow, S.J., 1993. Partial Differential Equations for Scientists and Engineers. Dover Publications, Inc.
7.Hargrave, B.T., 1969. Similarity of oxygen uptake by benthic communities. Limnol. Oceanogr., 14: 801-805
8.Henderson-Seller, B., 1984. Engineering Limnology. Pitman Publishing, London.
9.Kuo, Y.C., Huang, L.H., and Tsai, T.L., 2008., Water Resour. Res., 44, W03414, doi:10.1029/2007WR006084.
10.Ozisik, M.N., 1993. Heat conduction. 2nd. John Wiley & Sons, Inc.
11.Stefan, H.G. and Fang, X., 1994., Dissolved oxygen model for regional lake analysis. Ecol. Modeling, 71：37-68
12.Vassilis Z.A., and Soultana K.G, 2003., Simulation of water temperature and dissolved oxygen distribution in Lake Vegoritis, Greece., Ecol model., 160, 39-53.
13.Wilfried, B., 1934., Evaporation into the atmosphere., Kluwer academic publishing.
14.林憲德, 蘇瑞泉, 台灣地區月平均日射量分布之研究, 氣象學報, 第32卷第3期, 1986.
15.唐榮澤, 臺南地區日射量之分析, 氣象學報, 第25卷第2期, 1979.
16.郭遠錦, 水平分區、垂直分層之三維地下水移流及延散計算法, 國立台灣大學碩士論文, 台北, 民國93年.
17.陳秋楊, 陳伯中, 王敏昭, 鄭克聲, 德基水庫集水區水質監測與管理計畫, 經濟部水利署, 民國89年