南崁溪因流經桃園市、蘆竹鄉及大園鄉等人口稠密地區，工商業發展迅速，造成各河段的污染程度日益嚴重，因此桃園市政府環境保護局發布桃園市南崁溪流域廢（污）水重金屬銅排放總量管制方式，劃定南崁溪流域為重金屬銅總量管制區，並指定總量管制區之特定承受水體為南崁溪。
本研究應用Water Quality Analysis Simulation Program (WASP)建立南崁溪水質模式，利用平均絕對百分誤差法(Mean Absolute Percentage Error, MAPE)進行評估。
研究結果顯示率定MAPE 18%，驗證 MAPE 27%，符合優良及合理預測，將Q75流量輸入至模式進行模擬，實際污染總量4月為45.3 kg/day、8月為45.5 kg/day，河川涵容能力為7.81 kg/day，並根據重金屬銅達成率統計後得知南崁溪流域達成率低於50%之污染熱點測站，以達成地面水體分類及水質標準（0.03 mg/L）為目標，利用情境分析來模擬熱點測站污染源削減後之河段水質變化，假設7種情境，分析南崁溪削減污染總量，情境一為污染源全數削減50%之污染量，其削減總量4月與8月皆為24.9 kg/day；情境二為污染源全數削減65%之污染量，其削減總量4月為29.5 kg/day，8月為29.6 kg/day；情境三為污染源全數削減70%之污染量，其削減總量4月為31.7 kg/day，8月為32 kg/day；情境四為僅針對東門溪、高速公路橋下排水、瓦窯溝排水污染源削減50%之污染量，削減總量4月與8月皆為14.7 kg/day；情境五為僅針對東門溪、高速公路橋下排水、瓦窯溝排水污染源削減60%之污染量，削減總量4月與8月皆為17.7 kg/day；情境六為僅針對東門溪、高速公路橋下排水、瓦窯溝排水污染源削減65%之污染量，削減總量4月與8月皆為19.1 kg/day；情境七為僅針對東門溪、高速公路橋下排水、瓦窯溝排水污染源削減70%之污染量，削減總量4月與8月皆為20.6 kg/day。
綜整以上七種情境之結果，其削減總量介於14.72-31.73 kg/day（4月）及14.72-31.8 kg/day（8月）之間，而在有限資源及影響範圍較小的概念下，本研究建議僅針對部分污染源進行削減，如情境七之狀況較為合適，且各河段之濃度較為平均。

Nankang River flows through the areas with high population density such as Taoyuan City, Luzhu Township and Dayuan Township. With the rapid development of industry and commerce, the pollution levels of various river sections have become increasingly serious. Therefore, the Department of Environmental Protection, Taoyuan has released the Total Maximum Daily Load (TMDL) Method of heavy metal copper in waste (sewage) water of Nankang River basin, delineated Nankang River basin of heavy metal copper control area, and designated it as the specifi endorsed water.
In this study, the Water Quality Analysis Simulation Program (WASP) was used to establish the water quality model of Nankang River and evaluated by Mean Absolute Percentage Error (MAPE).
The results showed that the MAPE of calibration and validation are 18% and 27%, both of them are good and reasonable prediction. Put the Q75 flow rate into the model for simulation, the actual total pollution of April and August are 45.3 kg/day and 45.5 kg/day, and the capacity of the river is 7.81 kg/day. With calculating of the achieving Rate of heavy metal copper, we can find out the pollution stations whose achieving Rate is than 50% as the hot spot stations in Nankang River. In order to reach the standard of surface water body classification and water quality (0.03 mg/L), the situation analysis was used to simulate the change of water quality after the reduction of pollution sources in hot spot stations. Here are seven situations are assumed to analyze the total amount of pollution reduction in Nankang River. In the first situation, the total amount of pollution reduced by 50% from the sources of pollution is 24.9 kg/day in April and August; in the second situation, the total amount of pollution reduced by 65% from the total amount of pollution sources are 29.5 kg/day in August and 29.6 kg/day in August; Situation 3 is the total amount of pollution that is reduced by 70% from the source of pollution and the total reduction are 31.7 kg/day in April and 32 kg/day in August. In the situation 4, it is only for Highway bridge drainage and Wayaogou drainage pollution sources to reduce the amount of pollution by 50%, the total reduction in April and August are 14.7 kg/day; Situation 5 is only for the East Gate Creek, Highway Bridge drainage, Wayaogou drainage pollution sources The reduction of pollution by 60% and the total reduction by 17.7 kg/day in April and August respectively. Situation 6 is the reduction of pollution by 65% only for Dongmen Creek and expressway bridge drainage, In April and August, both of them are 19.1 kg/day. Situation7 is the reduction of 70% of the pollution from the drainage of Dongmen Creek and the highway bridges and drainage sources and the total reduction was 20.6 kg/day.
Integrate the results of the seven situations, the total amount of reduction of April and August are 14.72-31.73 kg/day and 14.72-31.80 kg/day. Under the concept of limited resources and smaller impact area, in this study, aimming for the specific pollution sources to reduce is the suggested strategy such as the situation 7,furthermore,the concentration in the various sections of the river is more close to the average.

目錄
目錄	I
表目錄	III
圖目錄	V
第一章 前言	1
1-1 研究緣起	1
1-2 研究目的	4
1-3 研究流程	4
第二章 文獻回顧	6
2-1 研究區域	6
2-1-1 南崁溪流域河川整治計畫	6
2-2 重金屬銅的特性與危害	8
2-3 水質模式	10
2-3-1 國內常用水質模式	11
2-3-2 模式評選	13
2-3-3 模式應用	13
第三章 研究區域與研究方法	16
3-1 研究區域	16
3-1-1南崁溪模擬河段介紹	16
3-1-2 環保署河川監測資料	18
3-2 研究方法	19
3-2-1 WASP模式建置	20
3-2-2 模擬結果判定指標	32
第四章 結果與討論	34
4-1 台灣地區河川重金屬水質標準達成率統計	34
4-1-1 各項重金屬水質標準達成率變化趨勢	34
4-1-2 南崁溪流域重金屬銅水質標準達成率變化趨勢	43
4-2 WASP水質模式模擬南崁溪流域重金屬	47
4-2-1模式率定及驗證	47
4-2-2率定及驗證結果	48
4-2-3河川涵容能力	49
4-2-4污染削減策略	53
4-2-5污染削減策略適用評估	69
第五章 結論與建議	71
5-1 結論	71
5-2 建議	71
參考文獻	72
 
表目錄
表 2 1 南崁溪流域相關計劃水質研究整理	7
表 2 2 台灣重金屬來源	9
表 2 3 各項標準之管制重金屬項目與限值	9
表 2 4水質模式之發展	10
表 2 5台灣地區常使用水質模式表	12
表 2 6國內研究統整	14
表 2 7國外研究統整	15
表 3 1 106年8月桃園市各區現住戶口統計表	18
表 3 2 南崁溪流域水體分類	19
表 3 3 WASP模式簡介	21
表 3 4 WASP模擬重金屬於水中流佈情形	25
表 3 5 支流排水資料	27
表 3 6 南崁溪流域各斷面之幾河特性表	31
表 3 7 南崁溪WASP 模式參數表	32
表 3 8平均誤差百分比判定指標表	33
表 4 1 2006-2015年各年度河川重金屬水質項目達成率統計表(單位：%)	34
表 4 2 地面水體分類及水質標準	35
表 4 3 2010-2015年各項重金屬達成率皆為100%之河川計有22條	35
表 4 4 2006-2015年54條河川重金屬水質項目達成率統計（單位：%）	37
表 4 5 105年度重金屬銅達成率低於85%之三大流域	43
表 4 6 桃園市南崁溪2010至2016年各測站之銅達成率(單位:%)	45
表 4 7桃園市南崁溪2010-2016年各測站相對位置	45
表 4 8南崁溪水質測站實測值	47
表 4 9南崁溪流域各斷面Q75整理表	50
表 4 10銅實際污染量	55
表 4 11污染源全數削減50%之污染量	57
表 4 12污染源全數削減65%之污染量	59
表 4 13污染源全數削減70%之污染量	61
表 4 14僅針對東門溪、高速公路橋下及瓦窯溝排水污染源削減50%之污染量	63
表 4 15僅針對東門溪、高速公路橋下排水及瓦窯溝排水污染源削減60%污染量	65
表 4 16僅針對東門溪、高速公路橋下排水及瓦窯溝排水污染源削減65%污染量	67
表 4 17僅針對東門溪、高速公路橋下排水及瓦窯溝排水污染源削減70%污染量	69

 
圖目錄
圖 1 1 我國11條重點污染整治河川分布	2
圖 1 2 桃園市南崁溪流域重金屬銅總量管制區	3
圖1 3 研究流程圖	5
圖 3 1南崁溪流域圖	17
圖 3 2 水質模式模擬流程圖	20
圖 3 3南崁溪流域主要支流排水幹渠魚骨圖	28
圖 3 4南崁溪流域水質測站位置圖	29
圖 3 5 南崁溪水質模式模擬網格圖	30
圖 4 1 桃園市南崁溪2006-2015年銅達成率變化趨勢	44
圖 4 2 桃園市南崁溪2010-2016年各測站銅達成率合鬚圖	45
圖 4 3 南崁溪流域水質測站位置圖	46
圖 4 4水體銅濃度模擬結果（105年4月7日）	48
圖 4 5水體銅濃度模擬結果（105年8月2日）	49
圖 4 6南崁溪流域斷面劃分與水質測站位置圖	52
圖 4 7 南崁溪4月Q75銅實際污染量模擬情形	54
圖 4 8南崁溪8月Q75銅實際污染量模擬情形	54
圖 4 9 污染源全數削減50%之銅濃度曲線（採用4月支排流實測值）	56
圖 4 10污染源全數削減50%之銅濃度曲線（採用8月支排流實測值）	56
圖 4 11 污染源全數削減65%之銅濃度曲線（採用4月支排流實測值）	58
圖 4 12污染源全數削減65%之銅濃度曲線（採用8月支排流實測值）	58
圖 4 13 污染源全數削減70%之銅濃度曲線（採用4月支排流實測值）	60
圖 4 14污染源全數削減70%之銅濃度曲線（採用8月支排流實測值）	60
圖 4 15 僅針對東門溪、高速公路橋下排水及瓦窯溝排水污染源削減50%之銅濃度曲線(採用4月支排流實測值)	62
圖 4 16 僅針對東門溪、高速公路橋下排水及瓦窯溝排水污染源削減50%之銅濃度曲線(採用8月支排流實測值)	62
圖 4 17 僅針對東門溪、高速公路橋下排水及瓦窯溝排水污染源削減60%之銅濃度曲線(採用4月支排流實測值)	64
圖 4 18 僅針對東門溪、高速公路橋下排水及瓦窯溝排水污染源削減60%之銅濃度曲線(採用8月支排流實測值)	64
圖 4 19僅針對東門溪、高速公路橋下排水及瓦窯溝排水污染源削減65%之銅濃度曲線(採用4月支排流實測值)	66
圖 4 20僅針對東門溪、高速公路橋下排水及瓦窯溝排水污染源削減65%之銅濃度曲線(採用8月支排流實測值)	66
圖 4 21僅針對東門溪、高速公路橋下排水及瓦窯溝排水污染源削減70%之銅濃度曲線(採用4月支排流實測值)	68
圖 4 22僅針對東門溪、高速公路橋下排水及瓦窯溝排水污染源削減70%之銅濃度曲線(採用8月支排流實測值)	68
圖 4 23 原污染情境與情境七模擬情境比較圖（採用4月支排流實測值）	70
圖 4 24原污染情境與情境七模擬情境比較圖（採用8月支排流實測值）	70

1.Allison, J.D., and Allison, T.L.. Partition Coefficients for Metals in Surface Water, Soil, and Waste.  U.S. EPA, Athens, GA. EPA/600/R-05-074 , 2005.
2.Di Toro,D.M.,Connolly,J.P.,“Mathematical Models of Water Quality in Large Lakes,Part 2:Lake Erie”,U.S. EPA-600/3-80-065,1980.
3.O’connor,D.J., and W.E. Dobbins, “The mechanism of reaeration innatural streams,”Transactions of the American Society of Civil Engineers, Vol.123, pp.59-72, 1956.
4.USEPA .Water Quality Analysis Simulation Program (WASP) User’s Manual. ,2008.
5.Streeter,H.W. and E.B.Phelp,“A Study of Pollution and Natural Purification of the Ohio  River”,Public Health Bulletin, NO.146,P.1-P.75,1925.
6.Sen Yu, Li He , Hongwei Lu, “An environmental fairness based optimisation model for the decision-support of joint control over the water quantity and quality of a river basin” , Renewable Energy, North China Electric Power University, China,2016.
7.陳國洲，新店溪流域污水廠配置對淡水河流域水質變化影響研究，淡江大學水資源及環境工程學系碩士論文，2001。
8.行政院環保署，流域別水污染總量管制暨配套措施之研究，2005。
9.周子翔，氣候變遷影響下濕地規劃之探討-以淡水河重要流域為例，國立台北科技大學建築系建築與都市設計碩士論文，2016。
10.黃聖授，高屏溪涵容能力之評估，國立中山大學環境工程碩士論文，2001。
11.張鈞凱，大漢溪流域污染調查與削減策略評估，國立台灣大學環境工程學系碩士論文，2005。
12.邱智慧，徵收水汙染防治費對河川水質之影響，國立台灣大學環境工程學研究所碩士論文，2006。
13.黃珮瑜，QUAL2K模式評估人工溼地改善大漢溪水質之研究，淡江大學水資源及環境工程學系碩士論文，2011。
14.陳宣佑，淡水河流域河川水質監測站檢討之研究，逢甲大學水利工程與資源保育學系碩士論文，2015。
15.經濟部水利署水利規劃試驗所，淡水河流域因應氣候變遷防洪及土砂研究計畫(1/2)，2012。
16.行政院環境保護署，環境影響評估河川水質評估模式技術規範，2011。
17.侯嘉琦，利用WASP推估河川涵容能力及研擬管理策略 -以屏東縣武洛溪為例，國立屏東科技大學環境工程與科學系碩士論文，2016。
18.陳伸安，二維水理棲地模式運用於南崁溪生態規劃之研究，國立中央大學土木工程研究所碩士論文，2006。
19.古沼格，運用生態工程改善河川水質並提昇河川自淨能力－以南崁溪為例，國立中央大學土木工程研究所碩士論文，2008。
20.行政院環境保護署，105年度河川污染整治策略評估及總量管制整合執行計畫（北區）期末報告，2016。
21.林永倫，南崁溪上游安但汙染源之調查及整治規劃，東南科技大學防災科技研究所碩士論文，2011。
22.行政院環境保護署，推動水區劃定計畫（北區），2012。
23.桃園縣政府環境保護局，桃園縣都會型河川污染整治及河道生態環境改善之評估規劃及細部設計計畫，2012。
24.行政院環境保護署，105年度河川污染整治策略評估及總量管制整合執行計畫（北區）期末報告，2016。
25.行政院環境保護署，重點河川污染整治策略評析及執行計畫期末報告，2015。
26.行政院環境保護署，重點河川污染整治策略評估與行動整合計畫期末報告，2014。
27.行政院環境保護署，重點河川污染整治推動行動整合計畫，2013。
28.張瑞奉，南崁溪及老街溪水質趨勢分析，國立中央大學環境工程研究所，2017。
29.行政院環境保護署，底泥污染來源及傳輸模式調查計畫－以重點河川為例，2013。
30.行政院環境保護署，102年河川水庫底泥污染傳輸模式應用及研析計畫，2014。
31.行政院環境保護署，底泥整治評估技術與對策支援專案計畫，2014。
32.行政院環境保護署，河川水庫底泥調查技術及整體環境危害評估計畫，2016。
33.曹婉如，高雄愛河沉積物重金屬元素分佈研究，國立臺灣海洋大學海洋環境資訊學系，2014。
34.李榮華，重金屬污染土壤之物質流分析與管理架構研擬－以桃園市污染場址為例，國立臺北大學公共事務學院自然資源與環境管理研究所，2015。
35.行政院環境保護署，底泥污染傳輸模式技術應用示範計畫，2016。
36.行政院環境保護署，運用河川水質模式評估開發行為對地面水體水質之影響，2017。