應用系統動力學於高屏溪流域水量水質改善策略之模擬分析
全文
(2) 應用系統動力學於高屏溪流域水量水質改善策略之 模擬分析 Application of System Dynamics on Accessing the Remedying Strategies of Water Quantity and Quality in Kaopin River. 研 究 生 : 李任馥. Student : Jen-Fu Li. 指 導 教 授: 張良正. Advisor : Liang-Cheng Chang. 國 立 交 通 大 學 土 木 工 程 學 系 碩 士 班 碩 士 論 文. A Thesis Submitted to Institute of Civil Engineering National Chiao Tung University in Partial Fulfillment of the Requirements for the Degree of Master of Science in Civil Engineering August 2007 Hsinchu, Taiwan, Republic of China. 中 華 民 國 九 十 六 年 八 月.
(3) 應用系統動力學於高屏溪流域水量水質改善策略之 模擬分析 學生:李任馥. 指導教授:張良正 博士 國立交通大學土木工程研究所. 摘要 工業發展和生活水準的提升,常使得用水需求及污染排放量不 斷地增加,因而造成缺水及水質污染等問題,高屏溪流域即為一典型 的例子。惟河川水質水量問題機制複雜,需要一個能綜和描述相關機 制的模式,始可進行較完善的整體評估。因此,本研究以高屏溪流域 為研究區域,應用系統動力學建立一個水量水質改善策略的綜合評估 模式,並以 VENSIM 為系統開發之軟體環境,模式的建立應用了系 統動力學的五大步驟即:問題的定義、系統的描述、因果回饋圖、模 型的建置及情境模擬等。研究結果顯示,只施行單一種水量或水質之 方案皆不能同時滿足水量和水質標準,而同時施行多種方案則可有多 種組合能同時滿足水量和水質標準,對此可行之組合,本研究進一步 計算其淨收益,以提供決策者更多的資訊做參考。研究發現若將改善 策略實施的時間點延後五年,則同時滿足水量水質標準的組合數就會 減少至兩種,此顯示高屏溪水量及水質問題之解決乃是刻不容緩。本 研究成果除了可作為高屏溪水量水質問題之改善參考外,並展示了應 用系統動力學於河川水量水質管理策略評估的可行性。. I.
(4) Application of System Dynamics on Accessing the Remedying Strategies of Water Quantity and Quality in Kaopin River Student:Jen-Fu Li. Advisor:Dr. Liang-Cheng Chang. Department of Civil Engineering National Chiao Tung University Abstract The water demand and usage in Kaopin area has increased significantly in recent years due to industrial growth and increased living standards. The increased water demand has raised the risk of water deficit. The raising water usage has also increased the amount of waste water disposal and thus causing the pollution of Kaopin river. Because water quality and water quantity are closely related, this study uses the System Dynamic Method(SD)to develop an integrated simulation model for strategically combining water quality and water quantity in the Kaopin river. System dynamic software, VENSIM, has been used to construct a simulation model. The model was developed through the procedures of problem definition, system description, causal loop feedback drawing, model implementation and scenario analysis. Investigation results demonstrate that a single strategy for remedying the water quality and quantity can not satisfy the water quality and quantity standard simultaneously, yet compound strategies can. Compound strategies have nine feasible combinations. The study calculates its net income in order to offer policymakers more options to make decision. If compound strategies are implemented until fifth year, the feasible strategies will be II.
(5) reduced to two. The results show that effective planning and management of water quality and quantity in the Kaopin river should be a priority for government. This study also demonstrated the practice of applying SD to access the remedying policies of water quality and quantity in the Kaopin river.. III.
(6) 謝誌 時光飛逝,似乎就在轉眼間,就要告別了這所美麗的校園,背起 滿滿的行囊,踏上人生另一段旅程。來這裡遇到了各個優秀的老師及 學長、同窗及學弟妹,以及那系館外那開的火紅的鳳凰花,交織成我 在交大土木所的一片片美好回憶…… 能夠有如此多有形無形的收穫,首先我要感謝在學業上啟蒙我良 多的張良正教授,指引出個人這次的研究方針,也成就了這本拙作; 以及逢甲大學的李漢鏗教授在水質領域上的教導;還要感謝楊朝仲學 長及何智超學長在研究計畫以及論文寫作上的指導,攜手帶我一起突 破各式各樣難關,陳宇文、葉明生、朱宏杰、蔡瑞彬學長不時提供專 業的建議;再來要感謝時常在研究室一起打拼到半夜的同窗好友昀 柏、君儀、浚瑋、全佑、嘉晉以及瘋狂阿給的戰友們,每回合的討論 都發出燦爛的火花令人難忘,回味再三;也感謝學弟妹婉如、怡釧、 正偉、玟豪、柏成、為善、昱維、敏威、韋圻及瀚聖一聲聲親切的問 候,讓研究室倍感溫馨。 最後,不能免俗的,要感謝最偉大的父母親和親愛的手足,以及 時常在身邊陪伴的女友僑耘,有了他們的加油打氣,才能在每回腸思 枯竭之時,產生新的信心與靈感。這次的研究不是結束而是一個新的 開始,我會帶著這樣的研究精神繼續面對未來的人生,再次由衷的感 謝以上的所有人。. IV.
(7) 章 節 目 錄. 中文摘要..................................................................................................... I 英文摘要.................................................................................................... II 謝誌........................................................................................................... IV 章節目錄.................................................................................................... V 表目錄......................................................................................................VII 圖目錄.................................................................................................... VIII 第一章 緒論...............................................................................................1 1.1 前言 ...............................................................................................1 1.2 文獻回顧 .......................................................................................2 第二章 研究方法論 ..................................................................................5 2.1 系統動力學之角色 .......................................................................5 2.2 系統動力學之原理 .......................................................................6 2.3 建模步驟 .......................................................................................8 2.4 因果回饋分析 ...............................................................................9 2.4.1 問題解決導向的因果回饋圖........................................... 10 2.4.2 水庫操作原則的因果回饋圖........................................... 11 第三章 研究範疇之界定 ........................................................................ 13 3.1 高屏溪水資源現況 ..................................................................... 13 3.2 高屏溪水質現況 ......................................................................... 14 3.3 問題與策略 ................................................................................. 17 3.3.1 水量改善的策略擬定 ...................................................... 17 3.3.2 水質改善的策略擬定 ...................................................... 18 第四章 系統的描述 ................................................................................ 19 V.
(8) 4.1 高屏溪流域基本資料 .................................................................. 19 4.2 高屏溪水量系統 .......................................................................... 21 4.3 高屏溪水質系統 .......................................................................... 27 第五章 因果回饋分析 ............................................................................ 28 第六章 系統動力模型建置 .................................................................... 31 6.1 流圖及模型建置 ......................................................................... 31 6.2 模型檢定 ..................................................................................... 44 第七章 情境模擬與分析 ........................................................................ 48 7.1 單一型策略模擬結果 ................................................................. 48 7.2 複合型策略模擬結果 ................................................................. 68 第八章 結論與建議 ................................................................................ 82 8.1 結論.............................................................................................. 82 8.2 建議.............................................................................................. 83 參考文獻................................................................................................... 84. VI.
(9) 表 目 錄. 表 3.2-1 高屏溪現行水體用途分類一覽表 .......................................... 16 表 3.2-2 地面水體分類及水質標準 ...................................................... 16 表 3.2-3 陸上水體分類及用途 .............................................................. 17 表 4.1-1 高屏溪各主支流基本資料 ...................................................... 21 表 4.1-2 高屏溪流域水系資料表 .......................................................... 21 表 4.2-1 自來水系統各水廠使用水資源與出水概況.......................... 25 表 4.2-2 高屏溪流域水質監測站與各主支流取水口關係表.............. 26 表 6.1-1 臺灣各地區污水量與用水量之比值表.................................. 34 表 6.1-2 高雄縣民國 81~90 年人口各年成長率.................................. 34 表 6.1-3 BOD 一般流達率建議值.......................................................... 35 表 6.1-4 畜牧業放流水標準 .................................................................. 38 表 6.1-5 高屏溪流域各種污染源比重表 ............................................ 39 表 6.1-6 高屏溪流域各種污染源排放表 ............................................ 40 表 7.1 單一型策略模擬結果 ................................................................. 49 表 7.2-1 策略調整之範圍 ...................................................................... 68 表 7.2-2 複合型策略模擬結果表(時刻初即加入改善策略).............. 69 表 7.2-3 複合型策略模擬結果表(第五年才加入改善策略).............. 70. VII.
(10) 圖 目 錄. 圖 2.1 系統動力學與決策支援系統建構流程圖 ................................... 5 圖 2.2 系統動力模式元件關係示意圖 ................................................... 7 圖 2.3 系統動力學建模與策略分析流程 ............................................... 9 圖 2.4-1 正向因果鍵 .............................................................................. 10 圖 2.4-2 負向因果鍵 .............................................................................. 10 圖 2.4-3 正向因果回饋環路 .................................................................. 10 圖 2.4-4 負向因果回饋環路 .................................................................. 10 圖 2.4.1 問題解決的因果回饋圖 .......................................................... 11 圖 2.4.2 水庫操作原則的因果回饋圖 .................................................. 11 圖 2.4.3 水量隨時間變化圖 .................................................................. 12 圖 4.1 高屏溪流域水系圖 .................................................................... 20 圖 4.2-1 水量系統圖 .............................................................................. 22 圖 4.2-2 甲仙堰位置圖 .......................................................................... 23 圖 4.2-3 高屏溪攔河堰位置圖 .............................................................. 24 圖 4.3 水量及水質系統圖 ..................................................................... 27 圖 5 高屏溪流域因果回饋圖 ................................................................ 30 圖 6.1-1 高屏溪流域水量流圖 .............................................................. 32 圖 6.1-2 廢水處理流程 .......................................................................... 36 圖 6.1-3 高屏溪流域水質流圖 .............................................................. 42 圖 6.1-4 高屏溪流域水量水質流圖 ...................................................... 43 圖 6.2-1 水質實測值,水質平均±一倍標準偏差................................ 45 圖 6.2-2 水質檢定結果圖 ...................................................................... 45 圖 6.2-3 水質驗證結果圖 ...................................................................... 46 VIII.
(11) 圖 6.2-4 高屏溪流域水質測站圖 .......................................................... 46 圖 6.2-5 高屏大橋水質測站資料 .......................................................... 47 圖 7.1-1 單一型策略模擬結果圖(現況) .............................................. 50 圖 7.1-2 單一型策略模擬結果圖(Q95) ................................................ 51 圖 7.1-3 單一型策略模擬結果圖(SI=1) .............................................. 52 圖 7.1-4 單一型策略模擬結果圖(Q85) ................................................ 53 圖 7.1-5 單一型策略模擬結果圖(Q65) ................................................ 54 圖 7.1-6 單一型策略模擬結果圖(Q55) ................................................ 55 圖 7.1-7 單一型策略模擬結果圖(接管速率 5%/年)......................... 56 圖 7.1-8 單一型策略模擬結果圖(猪隻削減 100%)........................... 57 圖 7.1-9 單一型策略模擬結果圖(猪隻削減 80%)............................. 58 圖 7.1-10 單一型策略模擬結果圖(猪隻削減 70%)........................... 59 圖 7.1-11 單一型策略模擬結果圖(抽取地下水 500 萬噸/月).......... 60 圖 7.1-12 單一型策略模擬結果圖(抽取地下水 800 萬噸/月).......... 61 圖 7.1-13 單一型策略模擬結果圖(抽取地下水 1000 萬噸/月)........ 62 圖 7.1-14 單一型策略模擬結果圖(抽取地下水 1200 萬噸/月)........ 63 圖 7.1-15 單一型策略模擬結果圖(坪頂淨水廠 30 萬噸/天)............ 64 圖 7.1-16 單一型策略模擬結果圖(坪頂淨水廠 60 萬噸/天)............ 65 圖 7.1-17 單一型策略模擬結果圖(南化淨水廠 30 萬噸/天)............ 66 圖 7.1-18 單一型策略模擬結果圖(南化淨水廠 60 萬噸/天)............ 67 圖 7.2-1 複合型策略模擬結果圖(時刻初即加入改善策略)(1)........ 71 圖 7.2-2 複合型策略模擬結果圖(時刻初即加入改善策略)(2)........ 72 圖 7.2-3 複合型策略模擬結果圖(時刻初即加入改善策略)(3)........ 73 圖 7.2-4 複合型策略模擬結果圖(時刻初即加入改善策略)(4)........ 74 圖 7.2-5 複合型策略模擬結果圖(時刻初即加入改善策略)(5)........ 75. IX.
(12) 圖 7.2-6 複合型策略模擬結果圖(時刻初即加入改善策略)(6)........ 76 圖 7.2-7 複合型策略模擬結果圖(時刻初即加入改善策略)(7)........ 77 圖 7.2-8 複合型策略模擬結果圖(時刻初即加入改善策略)(8)........ 78 圖 7.2-9 複合型策略模擬結果圖(時刻初即加入改善策略)(9)........ 79 圖 7.2-10 複合型策略模擬結果圖(第五年才加入改善策略)(10).... 80 圖 7.2-11 複合型策略模擬結果圖(第五年才加入改善策略)(11) .... 81. X.
(13) 第一章 緒論 1.1 前言 近年來台灣地區因為經濟的快速發展,導致用水需求不斷地成 長,使水資源的調度面臨比以往更嚴峻的考驗,因此,水資源的永續 利用對於台灣未來的經濟與社會發展有非常重大的影響,因為缺水會 影響產業的投資意願,而水的污染則會衝擊生態環境。所以,水質與 水量有效的規劃與管理應是刻不容緩的工作。此外,水質和水量的互 動關係非常密切,例如:若為了提供日漸成長的公共用水需求,而增 加攔河堰的引水量,將嚴重的影響河川的流量,而河川流量的減少則 會降低河川的涵容能力,河川涵容能力的降低,將會提高河川水污染 的程度,尤其在枯水期特別明顯。 基於前述,在規劃任何解決缺水問題的策略時,不僅必須同時考 量其對水質面的影響,也要提供其他預防或改善河川水質的配套措 施。很不幸地,在台灣水質和水量的問題通常是被分開考量,並且水 量的滿足與否比水質的好壞更被重視,進而間接忽視了環境生態的維 護,使環境及生活居住的品質不斷地惡化,此乃重水量而輕水質之策 略所導致。以高屏溪為例,高屏溪為台灣流量最大的河川亦為南部地 區最重要的河川,惟高屏地區水質不佳是長期以來的問題,而且由於 水資源有限,限制了高屏地區的經濟發展,水質問題與水資源的供給 是高屏地區急迫需要解決的問題。 有鑑於此,本研究將以水量水質整體考量的觀點,分析高屏地區 所面臨的水質水量問題,並試圖以系統思考的角度,研擬出可行而有 效的水質水量改善策略,藉以協助達成區域水資源之永續經營管理。. 1.
(14) 1.2 文獻回顧 有效的永續水資源規劃與管理,水量和水質的同時考量是必要 的,此一議題國內外相關研究如下 : L. Gabriel et al. (2000) 架構出一 個同時結合了水量模式 (MOSIM) 和水質模式 (QUAL2E-UNCAS) 的水資源管理決策支援系統,並將其運用在 Piracicaba 流域;Sharon et al. (2001) 整合水量模式 (MODSIM) 和水質模式 (HEC-5Q) ,並 將其運用在 Klamath 流域,藉由此一整合模式的模擬結果,篩選出 能改善水質並滿足水量要求的策略。李魁裕 (2001) 同時進行高屏溪 水量改善和水質改善方案的成本效益分析,特色為使用排放轉換函數 結合水質改善方案的成本效益分析,排放轉換函數是以 QUAL2E 水 質模式模擬。 上述研究均為水量模式與水質模式的搭配模擬,彼此間的連接與 整合純為資料的相互傳遞,亦非水量與水質各單元變數因果連接的方 式,因此不易瞭解水量與水質系統實際互動影響的方式,僅能就模擬 結果進行參數敏感度分析來解釋其變化趨勢。為了解決這樣的問題, 本研究將採用系統動力學來建置水質水量整合模擬模型。 傳統上進行嚴謹的水資源系統分析時,一般皆會針對問題繪製系 統圖與流程圖,以釐清問題的解決方式與處理邏輯。系統圖是將系統 本身以圖形表達,定義了問題的結構與範圍;流程圖是為達成特定目 的而實行的一連串邏輯相關的工作以圖形表現,表達了問題處理步 驟。系統圖與流程圖各自表達了解題所需的部份資訊,例如:水資源 調配的系統圖可以展現各供水與需水節點的空間相對位置及節點間 水源流動的順序,但是各節點放水原則及計算方式就得從流程圖才能 得知。同樣地,只單看流程圖雖然可以瞭解水源運用操作及演算的順 2.
(15) 序,但是無法得知流程圖各變數於實體系統的相對位置與連接型態。 系統動力學是由美國麻省理工史隆管理學院的 J.W Forrester 教 授於 1960 年代所發展的一門科學,結合了控制(Cybernetics)、系統 論(System Theory)、資訊理論(Information Theory)、決策論(Decision Theory)、電腦模擬(Computer Simulation)等理論成為管理的新方法、 新工具和新概念。系統動力學為描述、探索和分析複雜系統內部流程 與訊息、組織範圍與策略的一種嚴謹的研究方法,其可透過定量化的 系統模擬分析來進行系統結構與行為之設計。系統動力學是處理訊息 回饋系統之動態行為的一種方法論,它提供一種實驗的、定量分析的 方法,因此對於極複雜的動態、回饋且具時間滯延(Time Delay)的問 題,能提供整體、長期且較週延的解決方法。非常適合於應用在與時 間演進有密切關係之課題研究。現今水質和水量的互動的問題,已具 備了複雜的時間隨動態變化的特性,因此適合利用系統動力學來解決 此類問題。除了強調探討系統隨時間變化,系統動力學另一重要之基 (1)藉由因果關係的確認來說明系 礎為因果之關連性,該意義有三: 統之問題; (2)藉由因果關係的確認將複雜之問題作簡潔而系統化之 表示; (3)藉由變數間之因果關係來說明系統之範圍。其中,因果關 係之表示旨在說明兩個變數間之關係為正向或負向,而無數量上之意 義,意即數量之大小並不會影響到因果關係鍵之存在與否。進而若將 一系列之因果鍵串接成因果回饋關係環路,則可將因果關係發展為正 向之因果回饋關係環路或負向之因果回饋關係環路;正向之因果回饋 關係環路,係表示環路上任何變數的變動,會使系統在下一週期朝同 方向加強其變動幅度,造成自我之強化,而負向之因果回饋關係環 路,則表示任何變數的變動,會使系統在下一週期朝反方向變動幅. 3.
(16) 度,變數產生抑制變動之效果,造成自我之規律。 系統動力學應用於水資源及水質的國內外文獻如下,Simonovic (1997)將系統動力學應用於埃及的水資源策略分析上,Stave(2003) 將 系統動力學應用於拉斯維加斯的水資源管理上,XU et al.(2002)應用 水 資 源 系 統 動 力 學 (WRSD) 分 析 中 國 大 陸 黃 河 流 域 的 永 續 性 , H.C.Guo(2001) 將系統動力學應用於中國大陸洱海(Erhai)湖流域的 水質規劃與管理。李孟璁(2002) 以高雄愛河為研究對象,利用河川 水質模式,配合上游集水區的人口變化、土地使用以及降雨型態的差 異,並使用系統動力學的軟體 STELLA 來建構模式模擬現況。黃鈺 珊(2001)建構高屏溪流域水資源使用與都市系統動力模型,並以情境 分析,針對當前高屏溪流域之水資源與都市發展政策進行政策試驗, 檢視其是否達到永續發展之理念。李明益(2001)建立河川流域系統與 水管理決策的系統動力學模式。童智偉(2000)則依據驅力-狀態-反應 之永續指標架構,探討流域性環境水質之永續力。以濁水溪流域為 例,配合 VENSIM 系統動力學軟體建立環境管理決策支援系統,進 行水質評量系統及動態模擬分析。詹麗梅(2001)利用系統動力學學理 論及 STELLA 系統動力學軟體建立區域供水系統之模擬模式,並以 大基隆供水區為例,探討不同供水策略對供水系統永續性之影響。 上述系統動力學應用於水資源領域的文章並無同時考量水質和 水量的互動影響,因此本研究將以系統動力學五大步驟-問題的定 義、系統的描述、因果回饋圖、模型的建置、情境模擬,來架構水質 水量互動的問題,並利用 VENSIM 來建置系統動力模型,最後透過 情境模擬來找出能同時解決水質水量問題的策略組合。本研究的應用 區域為高屏溪流域。. 4.
(17) 第二章 研究方法論 2.1 系統動力學之角色 系統動力學建構在整個支援平台過程所扮演的角色如圖 2.1 所 示。首先針對欲研究之議題進行相關水文、地文與人文資料的蒐集, 接著藉由資料礦掘、模擬模式或優選模式等技術對相關資料進行資訊 擷取的工作。將各領域所獲得的資訊與知識,透過系統動力學進行相 互間的組織與串接。建構一個決策支援平台,最後將可進行即時的情 境模擬,模擬所產生的分析結果與相支援相關決策者如單位主管或透 過委員會進行溝通、協調與可行性之研究後,將能產生一兼具公眾利 益與環境品質的決策。. 圖 2.1 系統動力學與決策支援系統建構流程圖 5.
(18) 2.2 系統動力學之原理 一、 系統動力學之概念 系統動力學(System Dynamics)是由美國麻省理工史隆管理學 院(MIT)的 J.W Forrester 教授於 1960 年代所發展的一門科學,結 合 了 控 制 (Cybernetics) 、 系 統 理 論 (System Theory) 、 資 訊 理 論 (Information Theory)、決策論(Decision Theory)、電腦模擬(Computer Simulation)等理論成為一體的管理新方法、新工具和新概念,用以研 究具資訊回饋、決策與行動之間存在時間滯延現象等多種因素間互動 關係所形成之複雜動態行為,適用於解決非線性關係、因果循環、資 訊回饋、時間滯延等動態性複雜之問題解決上,且與其他應用於解決 一般問題的靜態、線性的工具有所不同。 根據 Jay W.Forrester 教授指出系統動力學學的特性如下: z 藉由變數間因果關係的確認以說明問題與問題內部系統的結 構,並將複雜的管理問題予以系統化的表示。 z 以因果關係為基礎,透過定義存量變數(Stock Variable)與流量 變數(FlowVariable)及其交互影響關係形成決策回饋環路。 z 導入訊息回饋的系統概念,在模式建構中將決策過程所運用 之訊息、訊息延遲及行為延誤的現象予以具體表達,並可模 擬真實系統的運作狀況。 z 可動態模擬一真實系統之訊息放大及干擾、時間延遲等現 象,建立系統的發展歷程,提供給系統分析者改變模式中之 參數組合以觀察系統行為的變化,提供管理決策之參考。 (河川流域水管理系統動力模式之發展與建立,2002) 系統動力學是處理訊息回饋系統之動態行為的一種方法論,它提. 6.
(19) 供一種實驗的、定量分析的方法,因此對於極複雜的動態、回饋且具 時間滯延(Time Delay)的問題,能提供整體、長期且較週延的解決方 法。近年來系統動力學應用日益廣泛,應用之課題包括了模擬地表水 污染、模擬水庫供水系統、生態系統族群變化、大氣化學與污染傳輸、 溫室氣體與全球暖化和現金流量的問題。由於系統動力學與時間演進 具有密切的連動關係,因此非常適合於應用在與時間演進有密切關係 之課題研究。現今流域永續經營管理的問題,已具備了所謂時間動態 複雜的特性,因此適合利用系統動力學來解決此類問題。因為系統動 力學模擬使用帶有滯延的微分方程組,能夠方便地處理非線性與時間 變化的關係,能夠做長期的、動態的分析與研究,適合大尺度、大範 圍的系統模擬。 二、 系統動力學的基本元素 系統動力模式在構成回饋系統各部門之相互關係及系統之回饋 環路結構時,一定包含兩種基本之變數,而這兩種變數是構成決策回 饋環路之必要與充分條件,其中第一種變數具有儲存之觀念,第二種 具有流動之觀念,再加上輔助其關係變化之第三種變數,構成三類主 要之系統組成變數元件,其分別為: (1)儲存(Level & Stock) 、 (2) 流動(Rate & Flow)、(3)助動(Auxiliary & Converter),另外,再 以箭頭連結表示各元件彼此間交互影響之關係,如圖 2.1 所示,其中 圖左為系統動力模式軟體 VENSIM 之示意圖,圖右則為 STELLA 之 示意圖。 Level. Stock. Rate. Flow. Auxiliary. 圖 2.2 系統動力模式元件關係示意圖 7. Converter.
(20) 積量(Stock)係表示某一系統變數在某一特定時刻的狀態,其 數值大小是累加了流入率(Inflow rate)與流出率(Outflow rate)的 淨差額所產生之結果,可說是系統過去活動結果之累積,如同水庫中 儲水量即屬於一儲存之概念;流量(Flow)則表示某種儲存變數變化 之快慢,代表著一種瞬間之行為,其數值多由儲存變數與助動變數之 交互關係來決定,如同水庫之放水量即屬於一流動之概念;輔助變數 (Auxiliary & Converter)則用來針對前兩種變數進行各變數間交互關 係之補充說明,類似一般計算方程式中之參數,如水庫中之放水規則 即屬於助動之概念。 三、 系統動力學之優點 系統動力學基於組成元件之因果關係建立系統架構,因此可透過 樹狀搜尋找到某變數之所有影響因數,有助於分析問題成因,以作為 未來解決問題之參考。 由於系統動力學將系統之組成元件與因果關係以積量、流量、輔 助變數和箭線四個元件表達,因此不同領域議題在基於四個共同組合 元件下,可容易進行整合分析。策略模擬結果與策略評估可即時以圖 形化呈現。. 2.3 建模步驟 模型建置模擬與策略分析流程如圖 2.3 所示。首先藉由「問題定 義」 、 「系統描述」 、 「因果回饋圖」以及「可能策略研擬」等建模步驟 來建置系統動力模型。在模型檢定無誤後,即可進行策略之情境模 擬,所有研擬策略的模擬結果透過「各問題代表性指標衝突分析」, 來初步篩選出可行策略。然而指標無法反映出策略施行時時間滯延的 問題,故於可行策略集合形成後,再進行策略之時間滯延效應分析, 分析的結果將提供決策者於協商過程時之重要參考。 8.
(21) 圖 2.3 系統動力學建模與策略分析流程 建模步驟的主要目的為設計模型建置之參考藍圖,其主要工作項 目包含了「問題定義」 、 「系統描述」 、 「因果回饋圖」以及「可能策略 研擬」。. 2.4 因果回饋分析 系統動力學重要基礎之一為因果之關連性,其意義有三: (1)藉 由因果關係的確認來說明系統之問題; (2)藉由因果關係的確認將複 雜之問題作簡潔而系統化之表示; (3)藉由變數間之因果關係來說明 系統之範圍。其中,因果關係之表示在說明兩個變數間之關係為正向 9.
(22) 或負向,如圖 2.4-1 及圖 2.4-2 所示,並而無數量上之意義,意即數 量之大小並不會影響到因果鍵之存在與否。進而若將一系列之因果鍵 串接成因果回饋關係環路,則可將因果關係發展為正向之因果回饋關 係環路或負向之因果回饋關係環路;正向之因果回饋關係環路,係表 示任何變數的變動,會使系統在下一週期朝同方向加強其變動幅度, 造成自我之強化,如圖 2.4-3 所示,而負向之因果回饋關係環路,則 表示任何變數的變動,會使系統在下一週期朝反方向變動幅度,變數 產生抑制變動之效果,造成自我之規律,如圖 2.4-4 所示。. -. +. B. A. B. A. 圖 2.4-1 正向因果鍵. 圖 2.4-2 負向因果鍵. + C. B. +. + C. B. +. +. A. 圖 2.4-3 正向因果回饋環路. A. 圖 2.4-4 負向因果回饋環路. 藉由因果回饋圖的繪製,我們可以清楚了解各變數的因果關係與 作用方向,並瞭解結構的基本特性。 2.4.1 問題解決導向的因果回饋圖 問題解決導向的因果回饋圖可以圖 2.4.1 為例說明,當供水量和 需水量產生差距時,表示供水不足而產生缺水量,可由缺水量計算出 年缺水指標(SI) ,策略的實施則根據實際之年 SI 和標準年 SI(=1) 之差異。. 10.
(23) 缺水指標 (SI) 差距 (缺水量). 產出 (供水量). 目標 (需水量). 狀態 (水庫水量). 中間產出 (可利用水量). 圖 2.4.1 問題解決的因果回饋圖 2.4.2 水庫操作原則的因果回饋圖 水庫操作原則的因果回饋圖如圖 2.4.2 所示,其表示水庫蓄水量 越多則可利用水量(入流量+水庫蓄水量)越多,當可利用水量越多時 水庫供水量就可以越多,但供水量越多的結果亦會導致水庫蓄水量的 減少,如此即構成了一個封閉的負向因果回饋關係環路,這表示水庫 蓄水量(累積變數)會隨著時間的前進,而趨近於某一個平衡水位。. -. 產出 (供水量) +. 狀態 (水庫水量) 入流量 +. 中間產出 (可利用水量). 圖 2.4.2 水庫操作原則的因果回饋圖. 11.
(24) 圖 2.4.3 水量隨時間變化圖. 12.
(25) 第三章 研究範疇之界定 高屏地區水質不佳是長期以來的問題,而且由於水資源有限,限 制了高屏地區的經濟發展,水質問題與水資源的供給是高屏地區急迫 需要解決的問題。. 3.1 高屏溪水資源現況 在水資源方面,以年逕流量來看,高屏溪是臺灣水資源最豐沛之 河川流域,但是由於雨量集中在豐水期(5~10 月) ,豐水期的流量為 枯水期之 8~10 倍,而高屏溪河川流短坡陡加上流域內又無大型水 庫,使得豐水期的水資源,有將近九成直接流入台灣海峽。目前高屏 溪流域內的地面水資源貯蓄設施僅澄清湖、中正湖、鳳山水庫三座, 以及以引水為主的高屏溪攔河堰,這些貯蓄設施的容量非常有限,僅 能蓄存豐水期的少部份逕流量,若要有效地貯蓄豐水期的逕流量,需 要有更多的蓄水設施。 為了因應高屏地區不斷成長的用水需求及增加南部地區的水資 源,經濟部水資源局在民國 78 年完成美濃水庫的可行性規劃,但受 到民眾及環保團體強烈的反彈及抗議,所以希望能以其他替代方案解 決高屏地區的水資源問題,例如興建瑪家水庫、南化水庫壩體加高、 曾文水庫越域引水以及開挖高屏大湖等。. 13.
(26) 3.2 高屏溪水質現況 高屏溪流域的水質問題,主要造成水質污染的原因是流域內養猪 戶排放大量未經由處理的畜牧廢水及高屏溪下游沿岸居民所排放的 家庭污水,而高屏溪的流量在乾季與濕季差異十分懸殊,在濕季豐水 期時,豐沛的水量可以稀釋排入河川中污染物的濃度,因此河川水質 尚可達到河川水體水質標準,但是在乾季枯水期時,河川流量減小, 無法稀釋大量的污染物質,而造成河川水質惡化。高屏溪的污染源主 要來自生活污水、工業廢水、畜牧廢水、垃圾滲出水及非點源污染, 以下針對各種污染源敘述如下: 一、生活污水 高屏溪流域之人口分佈不均,人口主要集中在流域下游,如屏東 市、旗山、美濃、大樹、高樹等地,而目前流域內人口密集之地區已 完成污水下水道之規劃且有部分已開始建置但目前接管率仍很低,所 以生活污水大多未經由下水道接管至污水廠處理而直接排放入河川 中,嚴重影響高屏溪之水質。目前流域內已規劃之污水下水道系統有 4 個,分別為屏東市系統、旗美系統、高樹系統、大樹系統及杉林系 統,預計將對流域內污染削減有所幫助。 二、工業廢水 高屏溪流域內之工業區主要分佈於河川下游,主要之工業區為屏 東、林園、內埔、大發、永安,高屏溪流域內之工業區,排放廢水之 工廠列管共有 867 家,其中大發及林園工業區的廢水並未排入高屏 溪,故實際將廢水直接排入高屏溪者只有 229 家,而金屬工業比例最 多,所以未來應以重金屬監測為主要工作目標。 三、畜牧廢水. 14.
(27) 高屏溪流域的畜牧業相當發達,其中以養猪業最多,高屏溪流域 養猪頭數在民國 85 年約有 200 萬頭,行政院於民國 88 年起實施養 猪離牧法,優先削減高屏堰上游養猪戶,而至目前為止,猪隻數約為 100 萬頭,雖然已削減一半的猪隻數,但因養猪廢水大多未經處理即 排入高屏溪中和養猪廢水之 BOD 為家庭污水之 BOD 的好幾十倍, 所以養猪廢水目前對高屏溪水質影響很大。 四、垃圾滲出水 台灣地區由於垃圾廠址取得不易,故多將垃圾棄置場或掩埋場設 置於河川兩岸高灘地,由於掩埋場滲出水處理不易,其滲出之污水造 成土壤、地下水及河川之污染。高屏溪流域之 24 個鄉鎮中,除了旗 山鎮及長治鄉的掩埋場距離河川比較遠外,其他地區的垃圾廠大都在 河川附近,而且都未做好滲出水的處理,故將來必須加強垃圾場滲出 水處理即遷移至行水區外。 五、非點源污染 水體污染的來源可分為點源污染與非點源污染,以高屏溪流域而 言,流域內主要非點源污染以林地及農地之農業迴歸水為主(高屏溪 涵容能力之評估,2001)。 高屏溪水體用途分類如表 3.2-1 所示,高屏溪流域各河段的水體 水質分類為旗山溪自發源地至甲仙淨水場取水口,隘寮溪發源地至三 地門橋,荖濃溪發源地至濁口溪河口、濁口溪自發源地至荖濃溪交會 口為甲類水體,荖濃溪自濁水溪交會口至旗山溪交會口、隘寮溪自三 地門橋至荖濃溪交會口及高屏溪主流自旗山溪和荖濃溪交會口為乙 類水體,高屏溪主流自牛稠溪交會口至河口為丙類水體;一般而言河 川水體分類和其用途之關係為:用於一級公共給水之河川為甲類水 15.
(28) 體,用於二級公共給水之河川為乙類水體,而用於三級公共給水之河 川為丙類水體,用於灌溉用水之河川為丁類水體,各類水體水質分類 標準及陸上水體分類及用途如表 3.2-2 及 3.2-3 所示。 表 3.2-1 高屏溪現行水體用途分類一覽表. 資料來源:水污染防治法 表 3.2-2 地面水體分類及水質標準(BOD) 水質項目\ 水體分類. 甲. 乙. 丙. 丁. 戊. 溶氧量(mg/L). >=6.5. >=5.5. >=4.5. >=2.0. >=2.0. 生化需氧量(mg/L). 1. 2. 4. /. /. 氨氮(mg/L). 0.1. 0.3. 0.3. /. /. 硝酸鹽氮(mg/L). 10. /. /. /. /. 總磷(mg/L). 0.02. 0.05. /. /. /. 資料來源:水污染防治法. 16.
(29) 表 3.2-3 陸上水體分類及用途. 資料來源:水污染防治法 一、一級公共用水:指經消毒處理即可供公共給水之水源。 二、二級公共用水:指需經混凝、沈澱、過濾、消毒等一般通用之淨水 方法處理可供公共給水之水源。 三、三級公共用水:指經活性碳吸附、離子交換、逆滲透等特殊或高度 處理可供公共給水之水源。 四、一級水產用水:在陸域地面水體, 指可供鱒魚、香魚及鱸魚培養用 水之水源; 在海域水體, 指可供嘉臘魚及紫菜類培養用水之水源。 五、二級水產用水:在陸域地面水體, 指可供鰱魚、草魚及貝類培養用 水之水源;在海域水體,指虱目魚、烏魚及龍鬚菜培養用水之水源。 六、一級工業用水:指可供製造用水之水源。 七、二級工業用水:指可供冷卻用水之水源。. 3.3 問題與策略 水質和水量的互動關係是非常密切的,舉例來說,若為了提供日 漸成長的用水需求,而增加攔河堰的引水量,將嚴重的影響河川的流 量,而河川流量的減少則會降低河川涵容能力,提高水污染的程度, 因此本研究同時討論水量改善及水質改善之方案;水量改善策略包括 了興建淨水廠與抽取地下水兩方案,水質改善策略包括了養猪離牧與 興建污水下水道系統兩方案。 3.3.1 水量改善的策略擬定 水量改善的策略擬定有興建淨水廠、抽取地下水兩方案 17.
(30) 1.興建淨水廠 本研究興建淨水廠有兩處,一是坪頂淨水廠,在坪頂淨水廠附近 再興建一個新的淨水廠來處理高屏溪豐水期的水量,而目前坪頂淨水 廠之處理水量為 130 萬噸/日,二是南化淨水廠,在南化淨水廠附近 再興建一個新的淨水廠來處理南化水庫水量以增加枯水期高雄台南 地區的供給量,而目前南化淨水廠之處理水量為 80 萬噸/日。 2.抽取地下水 本研究假定屏東平原尚有多餘地下水資源,抽取地下水以滿足各 個時期的需求量。 3.3.2 水質改善的策略擬定 目前政府所規劃的高屏溪水質改善策略主要為興建污水下水道 系統、養猪離牧政策、課徵水污染防治費、遷移行水區內的垃圾場、 興建聚落小型污水處理廠等,以下僅就養猪離牧政策及興建污水下水 道系統此二方案進行策略分析。 1.養猪離牧政策 行政院已於 87 年施行養猪離牧政策以去除流域內養猪廢水的排 放,優先削減高屏溪攔河堰上游約 51 萬頭猪隻,而至目前猪隻數尚 有 100 萬頭,所以本研究之策略將削減高屏溪之猪隻數。 2.興建污水下水道系統 高屏溪流域內人口主要集中在下游,如屏東市、旗山、美濃、大 樹、高樹等地,而本研究之研究範圍未包含屏東市,因此策略分為興 建旗美、大樹及高樹等地之污水下水道系統。. 18.
(31) 第四章 系統的描述 4.1 高屏溪流域基本資料 高屏溪流域水系圖如圖 4.1,發源於中央山脈玉山附近,向南流 經高雄、屏東兩縣,於新園鄉注入台灣海峽,東以知本溪、卑南溪及 秀姑巒溪為鄰,北以濁水溪為鄰,西以曾文溪為鄰,南接東港溪與林 邊溪,流域面積為 3,256.85 平方公里,集水面積為 812.03 平方公里, 是全台灣流域面積最大的河川。高屏溪全長 170.9 公里,主要支流為 荖濃溪呈東北西南流向,匯流於隘寮溪,之後繼續流至嶺口匯入旗山 溪後合稱為高屏溪,高屏溪有武洛溪和牛稠溪匯入,於屏東新園鄉匯 入臺灣海峽。高屏溪流經過高雄縣及屏東縣 24 個鄉鎮市,流域內人 口約為 60 萬人,提供了大高雄地區農業、工業及公共用水的水源。 高屏溪之主要支流為旗山溪、荖濃溪、隘寮溪、濁口溪及美濃溪,各 支流發源地、流域面積、長度,與高屏溪流域平均坡度及流經區域之 基本資料如表 4.1-1 及表 4.1-2 所示。. 19.
(32) 圖4.1 高屏溪流域水系圖 圖4.1 高屏溪流域水系圖 資料來源:行政院環保署網站. 20.
(33) 表4.1-1 高屏溪各主支流基本資料. 資料來源:行政院環保署網站 表4.1-2 高屏溪流域水系資料表. 資料來源:經濟部水利署網站. 4.2 高屏溪水量系統 本研究考量之高屏溪流域水源調配系統主要是供應台南和高雄 地區之需求,然其主要水源除來自高屏攔河堰外南化水庫亦是重要供 水來源,因此水源調配所考量之系統必須包含高屏溪和南化水庫才算. 21.
(34) 完整,如圖 4.2-1。. 圖 4.2-1 水量系統圖 高屏溪流域上實體的水工結構物分述如下: 一、南化水庫 南化水庫位於台南縣南化鄉曾文溪之支流後堀溪上(如圖 4.2-2) ,水庫相關設施完成於民國 82 年,於民國 83 年 7 月開始出水, 目前正進行旗山溪越域引水工程,水庫集水區面積 104 平方公里, 水庫容量 158.05 百萬立方公尺,水庫主要功能為蓄水供給嘉義台南 高雄等地區民生用水,是擔負台南高雄地區公共給水之重要水庫。 二、南化水庫旗山溪越域引水工程 南化水庫越域引水工程係在旗山溪甲仙橋上游約 500 公尺處(如 圖 4.2-2) ,興建甲仙攔河堰引旗山溪豐水期(每年 6 月至 10 月)之 水量至南化水庫,引水地點集水區面積 407.89 平方公里,型式為混 凝土垂直跌落式堰長 120 公尺,隧道長 3,054 公尺,導水引入南化水 庫,隧道最大輸水量為 30cms。甲仙堰工程已於民國 88 年 4 月完工, 南化水庫越域引水工程完工後,可增加南化水庫供水量每日約 30 萬 22.
(35) 噸,提高南化水庫使用率,並可有效改善台南、高雄之缺水及水質不 佳問題 。南化水庫若與高屏溪攔河堰聯合運用,其供水量可由 80 萬 CMD 提高至 130 萬 CMD,於豐水期時,高屏溪攔河堰計畫出水 100 萬 CMD,南化水庫出水 30 萬 CMD,枯水期時,南化水庫計畫出水 130 萬 CMD,而高屏溪攔河堰則不出水。. 圖 4.2-2 甲仙堰位置圖 資料來源:台灣南部地區跨流域水資源統籌調配 動態模式建立之研究,2001 三、高屏溪攔河堰 高屏堰之堰址位於高雄縣大樹鄉高屏鐵路橋上游 2 公里處(如圖 4.2-3) ,為一長 674 公尺固定堰用以抬高水位,可利用豐水期之水量。 高屏堰與南化水庫聯合運用,豐水期可出水 130 萬 CMD 供應台南及. 23.
(36) 高雄地區。. 圖 4.2-3 高屏溪攔河堰位置圖 資料來源:台灣南部地區跨流域水資源統籌調配 動態模式建立之研究,2001 四、淨水廠 高屏地區之自來水飲用水源包含高屏溪、東港溪及地下水三部 份,其中以高屏流域為水資源的自來水廠,包括坪頂、翁公園、澄清 湖及拷潭淨水廠等,表 4.2-1 為高雄自來水系統各水場使用水資源與 出水概況,表 4.2-2 則為高屏溪流域水質監測站與主支流去水口關係 表。 24.
(37) 表 4.2-1 自來水系統各水廠使用水資源與出水概況. 資料來源:高屏溪涵容能力之評估,2001. 25.
(38) 表 4.2-2 高屏溪流域水質監測站與各主支流取水口關係表. 資料來源:高屏溪涵容能力之評估,2001. 26.
(39) 4.3 高屏溪水質系統 本研究考量之水質系統為高屏溪流域上游至高屏大橋,如圖 4.3 虛線部分。. 高屏大橋. 圖 4.3 水量及水質系統圖. 27.
(40) 第五章 因果回饋分析 本研究考量之高屏溪水量及水質問題的因果回饋圖如圖 5 所 示,其可進一步說明如下: 一、 水量面 圖 5 中藍色線為高屏溪流域水量面因果回饋圖,本研究考量之高 屏溪流域水源調配系統主要是供應台南和高雄地區之需求,然其主要 水源來自高屏堰和南化水庫,因此高屏溪流域水量面因果回饋圖以高 屏堰和南化水庫為主說明;高屏堰-高屏堰供水量受到兩個變數的影 響,高屏堰入流量和民生用水需求,民生用水需求先由高屏堰引水至 坪頂淨水廠處理後供給(高屏堰供水),高屏堰供水量亦受到上游入 流量(高屏堰入流量)大小的影響,豐水期可引較多的水量來滿足需 求,而高屏堰供水量越多則高屏溪剩餘水量(高屏溪流量)就會越少, 高屏堰入流量為甲仙堰入流量加上側流量,也就是旗山溪之流量扣除 甲仙堰越域引水量後再加上荖濃溪之流量,民生用水需求和高屏堰供 水量之不足量為高屏堰供水量不足量;南化水庫-高屏堰供水量不足 量再由南化水庫放水至南化淨水廠處理後供給(南化水庫供水量), 南化水庫供水量受到水庫可利用水量影響,水庫可利用水量越多則南 化水庫供水量越多,水庫可利用水量受到三個變數的影響分別為南化 水庫上游入流量,甲仙堰越域引水量和水庫水量,而水庫操作原則的 因果回饋環路於南化水庫供水量. 水庫水量. 可利用水. 量,其關係為南化水庫供水量越多則水庫水量越少,為負向因果鍵, 水庫水量越少則水庫可利用水量越少,為正向因果鍵,水庫可利用水 量越少則南化水庫供水量越少,為正向因果鍵,此因果回饋環路中有 一個負向因果鍵和兩個正向因果鍵,所以此為一個負向因果回饋環 路,則表示任何變數的變動,會使系統在下一週期朝反方向變動幅 28.
(41) 度,變數產生抑制變動之效果,造成自我之規律;而民生用水需求由 高屏堰供水和南化水庫供水後仍不足則由地下水抽水供給,地下水抽 水量受到地下水可用水量影響;所以水量面之總供水量為高屏堰供水 量、南化水庫供水量和地下水抽水量之總和,若民生需求由高屏堰供 水、南化水庫供水和地下水抽水仍然不足則產生總缺水量,再由總缺 水量計算出年 SI,總缺水量到年 SI 之箭頭上多了兩條直線表示年 SI 為總模擬時刻之最終值。 二、 水質面 圖 5 中橘色線為高屏溪流域水質面因果回饋圖,河川水質(BOD) 計算為排入高屏溪總流達量(上游)加上排入高屏溪總流達量(下游) 除上高屏溪流量後得出,高屏溪總流達量(上游)包含了旗山溪、荖 濃溪和隘寮溪之污染量,排入高屏溪總流達量(下游)為高屏溪主流, 排入高屏溪總流達量(下游)為高屏溪污染排放量乘上流達率,高屏 溪污染排放量為污染產生量扣除總處理量(現況處理量+溼地處理 量),污染產生量為家庭污水、畜牧廢水和工業廢水之總和,現況處 理量為流域上現有之污水處理廠之處理量,溼地處理量為武洛溪溼地 之處理量,年平均 BOD 為總模擬時刻之各月平均值。. 29.
(42) 圖 5 高屏溪流域因果回饋圖. 30.
(43) 第六章 系統動力模型建置 6.1 流圖及模型建置 根據第三章、第四章及第五章所述之問題的定義、系統的描述及 因果回饋圖,本章節將就高屏溪流域水量和水質系統利用 venism 的 四個基本元素進行流圖與模型建置。 高屏溪流域水量流圖建置,首先必須了解高屏溪流域上有哪些實 體的水工結構物須被考量,這些實體的水工結構物將用 venism 的基 本元素「積量」(Level)來表示,例如水庫與攔河堰,而水的流動則用 「率量」(Flow)來表示,例如攔河堰引水入水庫之流量,高屏溪流域 水量系統包含高屏溪和南化水庫才算完整,高屏溪流域水質系統包含 高屏溪流域整體,相關說明分述如下: 高屏溪流域水量流圖建置如圖 6.1-1,流圖右半邊入流量 1(高屏 溪支流旗山溪)由甲仙堰越域引水至南化水庫,越域引水後之剩餘水 量和入流量 2(高屏溪支流隘寮溪)匯流為高屏堰入流量,高屏堰入 流量=入流量 1(高屏溪支流旗山溪)-越域引水+入流量 2(高屏 溪支流隘寮溪),高屏堰在豐水期時引水至坪頂淨水廠,引水至坪頂 淨水廠之剩餘水量再扣除至澄清湖淨水廠之抽水量為高屏大橋之流 量,高屏溪流量=高屏堰入流量-引水量,高屏大橋流量=高屏溪流 量-至淨水廠;流圖左半邊之可利用水量為入流量 3(南化水庫上游 入流量)加上越域引水量和南化水庫之水量,可利用水量=入流量 3 +越域引水+南化水庫,南化水庫供水量由高屏堰供水不足量(民生 用水需求和高屏堰供水量之差距)所決定,高屏堰供水不足量=民生 用水需求-高屏堰供水量,若南化水庫供水後仍不足則由地下水抽水 量供水,地表水缺水量=高屏堰供水不足量-南化水庫供水量,總供 水量(高屏堰供水量+南化水庫供水量+地下水抽水量)和民生用水 31.
(44) 需求之差距為總缺水量,總缺水量=民生用水需求-總供水量,再由 總缺水量計算出年 SI。. 圖 6.1-1 高屏溪流域水量流圖 高屏溪流域水質系統流圖及模型建置所需相關資訊一、 污染量之分析及推估 以八十七年生化需氧量(BOD5)每天流達污染量約 30.2 公噸為高 屏溪涵容能力 9.2 公噸/天之 3.3 倍,而依生化需氧量每天流達污染量 所佔比例計,家庭污水佔 49﹪為最大,其次以畜牧廢水佔全流域 41 ﹪,足見本流域水體之主要污染來源為與家庭污水畜牧廢水(高屏溪 河川水污染防治計劃可行性之檢討,1991),本研究中之污染量推估 與輸入,以家庭污水、工業廢水和畜牧廢水三類污染源為主。本研究 32.
(45) 對污染物推估之主要項目為 BOD。 1. 家庭污水 (1) 推估資料收集分析 家庭污染量推估為家庭污水量乘上家庭污水之 BOD 濃度;家庭 污水量之計算是參考每人每日之用水量,由於每人每日污水產生量與 每人每日用水量成正比關係,一般而言,污水量之估計是以用水量乘 上一適當之係數,可推估出每人每日之污水產生量。表 6.1-1 係台灣 各地區採用之用水量與污水量之比值表,比值介於 0.68~0.9,然而台 灣大部分地區採用污水量與用水量之比值為 0.8 來推估每人每日污水 量,故本研究將採用平均值 0.8 作為推估係數;此外,根據「臺灣省 自來水公司第七區管理處」資料顯示,民國 90 年高雄區自來水每人 每日平均配水量為 342 公升,而用水量成長分析則根據經濟部水資會 「台南、高雄、屏東等區自來水需水量估計」的報告, 90 年至 100 年每人每日需水量年增加值為 3.5 公升/年;根據「高雄市污水下水道 系統檢討規劃」,每人每日平均用水量在民國 120 年為 360 至 400 公 升;一般家庭污水之 BOD 濃度約介於 150-250mg/L 間,估計於民國 90 年時,排放污水之 BOD 濃度約為 200 mg/L(河川流域管理- WASP 水質模式評估愛河之整治方案,2003),在配合化糞池處理效率或下 水道接管率及污水處理廠處理效率的因子考量後,可推得最終之家庭 污染排放量,排放量再乘上合理之流達率,即為流達量,化糞池處理 效率一般約為 30%。 人口成長推估來源為高雄縣統計資料庫,高雄縣民國 81~90 年人 口各年成長率如表 6.1-2,10 年內人口成長率平均為 0.837%。. 33.
(46) 表 6.1-1 臺灣各地區污水量與用水量之比值表. 資料來源:河川流域管理- WASP 水質模式評估愛河之整治方案,2003 表 6.1-2 高雄縣民國 81~90 年人口各年成長率. 資料來源:高雄縣政府網站 (2) 污染流達率 一般而言,污染排放量從源頭到河川時污染量會有一定比例上的 減少,所以實際到達河川之污染量為污染流達量,而污染流達量與污 染排放量的比值即為流達率,流達率可以實測之方式獲得,其步驟如 下: 34.
(47) z 現地測量排水幹線排入河川前的水質和流量,推算污染流達 量。 z 估算該集污區內各類污染源之污染排放量。 z 再由污染流達量除以污染排放量為集污區之流達率。 (河川流 域管理- WASP水質模式評估愛河之整治方案,2003) 由於實測流達率非常浪費時間及人力,所以一般多參考文獻,而 本研究參考資料如表 6.1-3,最後流達率的值再經由模式的檢定來修 正。 表 6.1-3 BOD 一般流達率建議值. 資料來源:河川流域管理- WASP 水質模式評估愛河之整治方案,2003 2. 畜牧廢水 (1)推估資料收集分析 畜牧廢水的污染種類包括猪、牛、雞、鴨,其中就污染產生量來 說,猪所產生的污染量最多,所以畜牧廢水的推估主要以猪為主。一 般猪場所產生的猪糞尿廢水,主要是由猪糞尿、猪舍清洗用水、飼料 殘渣與猪隻飲用水所組成,而猪糞尿的排泄量也會因飼養方式、猪舍 地面形式與猪隻體重的變化而有所不同(猪場廢水處理規劃線上諮詢 系統,2003),除了猪隻所排放的糞水外,每清掃猪舍之沖洗水量, 亦為計算養猪廢水之重要指標;由於每一種處理設備在設計時各有不 35.
(48) 同之淨化功能,因此在要求的污水濃度也就不盡相同,為了發揮已有 的處理設備功能,因此於每日沖洗水量上必須把握正確,台灣之養猪 戶習慣以水沖洗地面,夏季並為猪隻沖涼,廢污之稀釋率約在 5~15 倍之間,每頭猪所排放的廢水量約為 30 L/day(河川流域管理- WASP 水質模式評估愛河之整治方案,2003),在水質方面,其生化需氧量 為 6700mg/L,由於農政單位的推廣與補助,稍具規模之養猪場皆有 設置污水處理設備,其設備型式以農政單位所推廣之「三段式猪廢水 處理設備」為主,此三段式為依序進行固液分離、厭氣發酵及活性污 泥處理,在 BOD 處理上,經過三階段處理後,去除率已達 90%。 (2)三段式養猪廢水處理系統 目前國內猪糞尿廢水的處理一般都採用三段式猪糞尿處理系統 ,其主要處理步驟包括固液分離、厭氣發酵及好氣處理,再加入計量 槽、調整池、初沉池、污泥回流及污泥處理,其過程簡示如圖 6.1-2 (猪場廢水處理規劃線上諮詢系統,2003):. 圖 6.1-2 廢水處理流程 資料來源:猪場廢水處理規劃線上諮詢系統,2003. 36.
(49) z 固液分離:利用固液分離機將原廢水分離成固體與液體部分,其 中固體部分可發酵做成堆肥,而液體部分則進行第二步驟之厭氣 發酵。BOD 與 SS 在此步驟分別可去除 20%及 40%。(猪場廢 水處理規劃線上諮詢系統,2003) z 厭氣發酵:利用厭氣微生物可分解有機物質之原理,以紅泥膠皮 做成厭氣發酵槽,處理固液分離後之廢水。此步驟之 BOD 與 SS 去除率都可達到 80%,此階段產生的沼氣可以作為能源。(猪場 廢水處理規劃線上諮詢系統,2003) z 好氣處理:利用曝氣機曝氣、增加溶氧,以好氣性微生物處理經 厭氣發酵後之廢水,在經沉澱後即可放流。養猪廢水之BOD與SS 去除率都超過90%以上(猪場廢水處理規劃線上諮詢系統, 2003)。 (3)放流水標準 猪隻所產生之畜牧廢水,雖然廢水量沒有工業廢水量多,但因為 有機物濃度高,另外在有機物質、懸浮固體物、重金屬等水質項目都 有很高的濃度,如不將其作妥善的處理會嚴重的污染河川水質,行政 院環境保護署為解決畜牧廢水對河川水質之衝擊於水污染防治法規 中規範畜牧業廢水之排放標準,台灣地區畜牧業放流水標準於表 6.1-4,可以看得出放流水標準日趨嚴格,例如在民國 76~81 年間畜牧 業放流水標準在大規模養猪場之生物需氧量(BOD)為 200 mg/L, 小規模養猪場為 400 mg/L,懸浮固體物(SS)分別為 300 mg/L 和 400 mg/L,而在民國 76~81 年間生物需氧量和懸浮固體物分別提高為 100mg/L 和 200mg/L 外,另外再增加化學需氧量(chemical oxygen demand, COD)為 400 mg/L;畜牧廢水以三段式猪糞尿處理法,即好 37.
(50) 氣處理、固液分離、厭氣發酵及三部份以達到畜牧業放流水標準。而 到民國 92 年,中華民國養猪協會與各地方協會及代表連署提案,要 求環保署修訂畜牧放流水標準,環保署評估養猪戶處理技術,修正化 學需氧量,亦即由現行的 250 mg/L 調整為 600 mg/L,而生化需氧量 及懸浮固體物,仍維持 80 及 150 mg/L(猪場廢水處理規劃線上諮詢 系統,2003)。 表6.1-4 畜牧業放流水標準. 資料來源:猪場廢水處理規劃線上諮詢系統,2003 3. 工業廢水 高屏溪流域內之工業區主要分佈於河川下游,主要之工業區為屏 東、林園、內埔、大發、永安,高屏溪流域內之工業區,排放廢水之 工廠列管共有 867 家,其中大發及林園工業區的廢水並未排入高屏 溪,故實際將廢水直接排入高屏溪者只有 229 家, 高屏溪流域各種 污染源比重如表 6.1-5,由表可知本研究區域內集污分區屬溪埔寮排 水至大樹、姑婆寮排水,其工業廢水比例甚小,故本研究以家庭污水 和畜牧廢水推估為主,工業廢水污染排放量(kgBOD/day)則參考表 6.1-6(行政院環保署網站)作為實際輸入值。. 38.
(51) 表 6.1-5 高屏溪流域各種污染源比重表. 資料來源:高屏溪涵容能力之評估,2001. 39.
(52) 表6.1-6 高屏溪流域各種污染源排放表. 資料來源:行政院環保署網頁 40.
(53) 高屏溪流域水質流圖建置如圖 6.1-3 ,河川水質(BOD)=排入 高屏溪總流達量(下游)+排入高屏溪總流達量(上游)/高屏大橋 流量;排入高屏溪總流達量(下游)=總流達量-溼地處理量(武洛 溪溼地),武洛溪溼地基地位置以南側堤防外水利用地,第一期面積 約 2 公頃,第二期面積約 18 公頃,分階段進行水質淨化處理程序, 第一階段處理採自然疊落、礫間接觸、生物接觸氧化及自然沈澱等方 式,以去除水中砂粒及部分懸浮固體,增加水中溶氧量,能初步削減 武洛溪水中的懸浮固體達五成,生化需氧量 (BOD)去除率也可達三 成,第二階段則以人工渠道方式,將溪水引流入三座自由水流系統的 人工溼地,經由溼地植物的過濾、淨化及吸附氮、磷後,進一步達到 淨化水質的目的,處理過後的水中懸浮固體及生化需氧量去除率,也 可分別大幅提升到六至七成以上; 水質流圖中各項目計算式如下: 總流達量=家庭污水流達量(kg)+畜牧廢水流達量(kg)+工業廢 水流達量(kg) 各污染源之流達量(kg)=各污染源之排放量(kg)× 流達率 家庭污水流達量(kg)=家庭污水排放量#1(kg)×流達率 C+家庭 污水排放量#2(kg)×流達率 D 家庭污水產生量 A=每人每日之用水量×0.8 家庭污水產生量 B=家庭污水產生量 A-化糞池處理 家庭污水產生量#1(kg)=家庭污水排放量(未處理)×家庭污水 BOD 濃度×(1-化糞池處理效率) 家庭污水排放量(未處理)=家庭污水產生量 B-流入污水廠 41.
(54) 家庭污水排放量#2(kg)=污水廠排放(已處理)×放流水標準#污水 下水道系統 畜牧廢水排放量(kg)=(畜牧廢水產生量-三階段式養猪廢水處理 系統)×放流水標準#畜牧 工業廢水排放量參考表 6.1-6(行政院環保署網站)作為實際輸入 值;高屏溪流域水量水質系統流圖如圖 6.1-4 所示。. 圖 6.1-3 高屏溪流域水質流圖. 42.
(55) 圖 6.1-4 高屏溪流域水量水質流圖. 43.
(56) 6.2 模型檢定 高屏溪水量水質系統動力模型建置完成後,必須經由模型檢定驗 證後才可實地應用,模型檢定指調整水量和水質的參數和係數,使實 測水質資料與模擬結果比較,兩者差異最小而模擬結果可以被接受之 過程;本研究水質測站選定為高屏大橋,高屏大橋位於高屏溪主流(如 圖 6.2-4) ,測站地址為高雄縣大樹鄉 1 號公路,高屏大橋水質測站資 料如圖 6.2-5,而在高屏溪流域最下游的水質測站為雙園大橋(如圖 6.2-4) ,測站地址為高雄縣大樹鄉 17 號公路,本研究未選定其為水質 檢定測站,原因為雙園大橋位置靠近出海口為感潮河段,會造成水質 ,首先, 模擬上的不準確;河川水質檢定的項目為生化需氧量(BOD5) 選定 5 年(1997、1998、2002、2003、2004)水質實測資料,計算這 5 年各月水質之平均值與一倍標準偏差如圖 6.2-1;選定與 5 年水質平 均值較接近之年份進行水質檢定,選定之年份為 2002 年,再來對 2002 年進行水量和水質的參數和係數調整,結果如圖 6.2-2,待模型檢定 完成後,為確認模型之可靠性需進一步驗證,驗證的年份為 1997、 1998、2004,結果如圖 6.2-3,驗證的標準為平均值正負一個標準偏 差,結果顯示大致上符合這個標準,所以此模式是可以被接受的。. 44.
(57) 圖 6.2-1 水質實測值,水質平均±一倍標準偏差. 圖 6.2-2 水質檢定結果圖. 45.
(58) 圖 6.2-3 水質驗證結果圖. 圖 6.2-4 高屏溪流域水質測站圖 資料來源:行政院環保署網頁 46.
(59) 圖 6.2-5 高屏大橋水質測站資料 資料來源:行政院環保署網頁. 47.
(60) 第七章 情境模擬與分析 本章節即利用第六章節所建置完成之高屏溪水量水質系統動力 模型進行策略模擬分析,其包含興建淨水廠、抽取地下水、養猪離牧 政策、興建污水下水道系統與生態基流量等單一型策略及複合型策 略,透過衝擊分析以評估其效益。本研究之策略評比標準為缺水程度 年 SI≦1,生化需氧量平均 BOD≦2 mg/l,同時滿足兩個標準方為可 接受之模擬方案,而淨成本這個項目是供決策者選出可接受之模擬方 案後之另一項參考項目。. 7.1 單一型策略模擬結果 表 7.1 為單一型策略模擬結果,目標年為民國 110 年,分別對興 建淨水廠、抽取地下水、養猪離牧政策、興建污水下水道系統與生態 基流量等進行單一策略之模擬,表格最上方為模擬時刻初(民國 91 年)之情況即在模擬時刻初的供水能力(坪頂淨水廠 110 萬噸/天, 南化淨水廠 80 萬噸/天)及污水處理情況下不做任何策略模擬至目標 年,相關模擬結果如圖 7.1-1~7.1-18,模擬結果發現年 SI 和平均 BOD 皆超出標準,各項水量及水質之單一型策略模擬結果顯示皆不能同時 滿足水量和水質標準,所以必須考量複合型策略來同時解決水量和水 質問題。. 48.
(61) 表 7.1 單一型策略模擬結果 策略. 年 SI. 平均 水質 收入 BOD 不合 格月 數. 成本. 淨收益 單位: 億元. 現況 Q75(26.18cms) (目標:民國 110). 3.8. 4.38. 128. 4.96. 0. 4.96. 生態基流量. Q95(9.61cms)(-63%). 0.7. 6.4. 128. 5.61. 0. 5.61. 3150cms(-41%). 1. 5.7. 128. 5.51. 0. 5.51. Q85(16.57cms)(-37%). 1.7. 5. 128. 5.35. 0. 5.35. Q65(35.71cms)(+36%). 5.2. 4.2. 128. 4.72. 0. 4.72. Q55(55.72cms)(+113%). 7.8. 4.1. 128. 4.42. 0. 4.42. 污水下水道. 接管速率 5%/年. 3.8. 4.36. 128. 4.96. 22.3. -17.37. 猪隻削減. 100%. 3.8. 1.1. 42. 4.96. 46. -41.04. 80%. 3.8. 1.7. 100. 4.96. 36.8. -31.84. 70%. 3.8. 2.0. 105. 4.96. 32.2. -27.24. 500(萬噸/月). 2.3. 4.38. 128. 5.2. 4.91. 0.3. 800. 1.9. 4.38. 128. 5.28. 6.48. -1.2. 1000. 1.5. 4.38. 128. 5.36. 8. -2.65. 1200. 1.2. 4.38. 128. 5.43. 9.48. -4.05. 30(萬噸/天). 3.1. 4.4. 128. 5.06. 8.7. -3.64. 60. 2.9. 4.4. 128. 5.08. 17.4. -12.31. 30. 3.68. 4.38. 128. 4.99. 8.7. -3.71. 60. 3.66. 4.38. 128. 4.99. 17.4. -12.40. 抽取地下水. 興建淨水廠 (坪頂) 興建淨水廠 (南化). 49.
(62) 圖 7.1-1 單一型策略模擬結果圖(現況). 50.
(63) 圖 7.1-2 單一型策略模擬結果圖(Q95). 51.
(64) 圖 7.1-3 單一型策略模擬結果圖(SI=1). 52.
(65) 圖 7.1-4 單一型策略模擬結果圖(Q85). 53.
(66) 圖 7.1-5 單一型策略模擬結果圖(Q65). 54.
(67) 圖 7.1-6 單一型策略模擬結果圖(Q55). 55.
(68) 圖 7.1-7 單一型策略模擬結果圖(接管速率 5%/年). 56.
(69) 圖 7.1-8 單一型策略模擬結果圖(猪隻削減 100%). 57.
(70) 圖 7.1-9 單一型策略模擬結果圖(猪隻削減 80%). 58.
(71) 圖 7.1-10 單一型策略模擬結果圖(猪隻削減 70%). 59.
(72) 圖 7.1-11 單一型策略模擬結果圖(抽取地下水 500 萬噸/月). 60.
(73) 圖 7.1-12 單一型策略模擬結果圖(抽取地下水 800 萬噸/月). 61.
(74) 圖 7.1-13 單一型策略模擬結果圖(抽取地下水 1000 萬噸/月). 62.
(75) 圖 7.1-14 單一型策略模擬結果圖(抽取地下水 1200 萬噸/月). 63.
(76) 圖 7.1-15 單一型策略模擬結果圖(坪頂淨水廠 30 萬噸/天). 64.
(77) 圖 7.1-16 單一型策略模擬結果圖(坪頂淨水廠 60 萬噸/天). 65.
(78) 圖 7.1-17 單一型策略模擬結果圖(南化淨水廠 30 萬噸/天). 66.
(79) 圖 7.1-18 單一型策略模擬結果圖(南化淨水廠 60 萬噸/天). 67.
(80) 7.2 複合型策略模擬結果 複合型策略模擬策略施行時間點分為兩種,第一種為在時刻初, 即加入改善策略,第二種為在第五年才加入改善策略,藉由 vensim 進行情境模擬包含策略施行之大小調整,策略調整之大小如表 7.2-1,在此範圍內調整策略規模大小至同時滿足水量和水質標準,第 一種情況模擬結果於表 7.2-2,表中列出多種可行性策略組合,及列 出各種可行性策略組合之淨收益,淨收益計算為收入減去成本,相關 的成本係數有自來水水價:7 元~11.5 元/噸(自來水公司) ,淨水廠處 理單位水成本:10.955 元/噸,地下水抽水固定成本:0.48 元/噸,地 下水抽水操作成本:0.9 元/噸(工研院能源資源所「地下水抽用成本 之研究」 ) ,淨水廠建造成本:17.4 億元/60 萬噸(中部地區水資源利 用整體檢討規劃);決策者可根據其需要去選定策略,第二種情況模 擬結果於表 7.2-3,表中列出兩種可行性策略組合,第二種情況因為 策略施行的時間點較晚,而在既有的範圍內可供選擇的組合數就會減 少,表示解決高屏溪之水量及水質問題是刻不容緩的,相關模擬結果 如圖 7.2-1~7.2-11。. 表 7.2-1 策略調整之範圍 MAX. MIN. 生態基流量. Q55. Q95. 猪隻削減. 100%. 0%. 抽取地下水. 1200 萬噸/月. 0. 增加淨水廠處理量. 60 萬噸/天. 0. 接管速率. 5%. 0. 68.
(81) 表 7.2-2 複合型策略模擬結果表(時刻初即加入改善策略) 策略組合. 年 SI. 平均 BOD. 水質不合 格月數. 收入. 成本. 淨收益 單位: 億元. 1. 抽取地下水(1200) +坪頂淨水廠(10)+猪隻削減(71%). 1.0. 2.0. 105. 5.47. 44.85. -39.38. 2. 抽取地下水(1100)+坪頂淨水廠(20)+猪隻削減(71%). 1.0. 2.0. 106. 5.46. 46.85. -41.39. 3. 抽取地下水(1000)+坪頂淨水廠(40)+猪隻削減(71%). 1.0. 2.0. 108. 5.45. 51.82. -46.37. 4. 抽取地下水(1000)+坪頂淨水廠(30)+坪頂淨水廠(30) +猪隻削減(72%). 1.0. 2.0. 108. 5.49. 57.95. -52.45. 5. 抽取地下水(1200)+猪隻削減(71%)+生態基流量(23.80cms). 1.0. 2.0. 105. 5.50. 41.74. -36.24. 6. 坪頂淨水廠(60)+坪頂淨水廠(60)+猪隻削減(79%) +生態基流量(18.02cms). 1.0. 2.0. 105. 5.50. 71.14. -65.63. 7. 坪頂淨水廠(30)+南化淨水廠(30)+猪隻削減(75%) +生態基流量(21.76cms). 1.0. 2.0. 105. 5.50. 56.02. -50.52. 8. 坪頂淨水廠(20)+猪隻削減(73%)+生態基流量(23.56cms)(-10%). 1.0. 2.0. 105. 5.46. 45.89. -40.43. 9. 抽取地下水(800)+坪頂淨水廠(40)+猪隻削減(74%) +生態基流量(23.56cms)(-10%). 1.0. 2.0. 106. 5.46. 51.35. -45.90. 69.
(82) 表 7.2-3 複合型策略模擬結果表(第五年才加入改善策略) 策略組合. 年 SI. 平均 BOD. 水質不合 收入 格月數. 成本. 淨收益 單位: 億元. 1. 抽取地下水(1200)+坪頂淨水廠(40)+南化淨水廠(20) +猪隻削減(97%)+生態基流量(23.98cms). 1.0. 2.0. 83. 5.52. 68.33. -62.81. 2. 抽取地下水(1200) +坪頂淨水廠(30) +南化淨水廠(10) +猪隻削減(97%)+生態基流量(23.56cms)(-10%). 1.0. 2.0. 85. 5.51. 62.67. -57.16. 70.
(83) 圖 7.2-1 複合型策略模擬結果圖(時刻初即加入改善策略)(1). 71.
(84) 圖 7.2-2 複合型策略模擬結果圖(時刻初即加入改善策略)(2). 72.
(85) 圖 7.2-3 複合型策略模擬結果圖(時刻初即加入改善策略)(3). 73.
(86) 圖 7.2-4 複合型策略模擬結果圖(時刻初即加入改善策略)(4). 74.
(87) 圖 7.2-5 複合型策略模擬結果圖(時刻初即加入改善策略)(5). 75.
(88) 圖 7.2-6 複合型策略模擬結果圖(時刻初即加入改善策略)(6). 76.
(89) 圖 7.2-7 複合型策略模擬結果圖(時刻初即加入改善策略)(7). 77.
(90) 圖 7.2-8 複合型策略模擬結果圖(時刻初即加入改善策略)(8). 78.
(91) 圖 7.2-9 複合型策略模擬結果圖(時刻初即加入改善策略)(9). 79.
(92) 圖 7.2-10 複合型策略模擬結果圖(第五年才加入改善策略)(1). 80.
(93) 圖 7.2-11 複合型策略模擬結果圖(第五年才加入改善策略)(2). 81.
(94) 第八章 結論與建議 8.1 結論 1. 本研究以系統動力學建立了一個高屏溪流域水量水質的綜合評估 模式,其中水質模擬範圍上起高屏溪流域上游下至高屏大橋,水 質計算假設模擬河段為單一河段 BOD 之質量守衡,並未隨著時間 和距離而衰減,且尚未考慮屏東市之點污染。 2. 由模擬結果得知,對只考量水量或水質單一方案之單一型策略, 其模擬結果顯示皆不能同時滿足水量和水質標準,所以必須以同 時考量水量及水質多種方案之複合型策略,來解決水量和水質問 題。 3. 本研究對各種策略組合皆計算其淨成本,複合型策略模擬結果顯 示有 9 種以上之策略組合能同時滿足水量和水質標準,其中又以 抽取地下水加上猪隻削減的策略組合之淨成本最低。 4. 若將各可行之複合型策略實施時間延後至第五年,則在給定的限 制條件下,可同時滿足水量水質要求的策略組合數將會大幅減少 至兩種,表示解決高屏溪水量及水質問題乃是刻不容緩。. 82.
數據
相關文件
(2000), “Enhancing entrepreneurial orientation research:Operationalizing and measuring a key strategic decision making process,” Journal of Management, Vol. Kenneth, Wehrung,
(1996)“Urban knowledge parks and economic and social development strategies”, Journal of Urban Planning and Development, 122(2):33-45. (1997)“Strategic behavior and
(1999), “Relationship Marketing and Data Quality Management," SAM Advanced Management Journal, Vol. (2004), “The Role of Multichannel Integration in Customer
在軟體的使用方面,使用 Simulink 來進行。Simulink 是一種分析與模擬動態
In order to ensure a sufficient quality and quantity of city water, the engineering quality of relevant works must be controlled so as to comply with the designed
Randall,D.,te al,(1997),Water Supply Planning Simulation Model Using Mixed-Integer Linear Programming, Jouenal of Water Resources Planning and
T., “The qualitative habitat evaluation index (QHEI), rationale, methods, and application,” Ohio EPA, Division of Water Quality Planning and Assessment, Ecological Assessment
T., “Criteria for establishing natural and scenic river study priorites,” Journal of Soil and Water Conservation, Nov.-Dec., pp. P., “Recreation and river type:Social-environmental