鑑於台灣工業快速發展與都市人口群聚，以及工業產生之高污染廢水不當排入河川，且台灣地區河川特性為豐枯水期流量差異明顯，上游坡度大，枯水季節之基流量（Base Flow）甚小，平均流速甚大，致再曝氣能力強，但河川實際上可涵容之污染總量，如生化需氧量(Biochemical Oxygen Demand, BOD)等皆不大，對於區域環境造成相當大之污染負荷，進而超過河川涵容能力(Natural Assimilative Capacity, NAC)之臨界值，並妨害水體正常用途。
本研究選用QUAL2E模式以朴子溪流域之水理、水文特性與歷年監測及實測水質資料進行模擬，並結合LINDO最佳化模式以河川水體使用分類所需考慮之條件與限制參數，設定河川流域之涵容能力為一個目標函數並使其最大化，進而推估朴子溪流域各河段之可允許污染量(Allowable Pollution Loading, APL)和未來水質改善策略。台灣大多數之汙染源都集中於廣域之河川流域，放流水標準之管制下，仍難達成水質環境基準之要求，為求有效改善必需對影響水質汙染負荷量做削減之規定。因此本研究之成果對建立河川廢水管制制度之施行可發揮最佳之貢獻，未來亦可為各工業區廢水管制提供參考策略。

In this paper, the performance of water quality modeling is investigated in the analysis of the natural assimilative capacity (NAC) of a river system. The new data available, along with the appropriate integrated model, allow for the analysis of the NAC of the entire river. The integrated model developed here emphasized the implementation of the QUAL2E simulation model and the appropriate allocation of waste loads (optimization model) in a subtropical area.
Simulations are calibrated and validated according to the data from the Water Quality and Flow Rate Monitoring Project using the QUAL2E model. On the basis of the optimal allocation model, the results of the NAC and allowable pollution loading in the Putzu River basin are discussed. The study provides a framework of the integrated model that can be used as a guideline to determine the NAC of subtropical rivers.

1.Barnwell, T. O., Jr., Brown, L. C., and Whittemore, R. C. (2004). Importance of field data in stream water quality modeling using QUAL2E-UNCAS. Journal of Environmental Engineering, 130, 643–647.
2.Bishop, A. B., and Grenney, W. J. (1976). Coupled optimization–simulation water quality model. Journal of the Environmental Engineering Division, 102, 1071–1086.
3.Bjorndal, E., Bjorndal, M., Pardalos, P. M., and Ronnqvist, M.(Eds.). (2010). Energy, natural resources and environmentaleconomics, New York: Springer.
4.Bowie, G. L., W. B. Mills, D. B. Porcella, C. L. Campbell, J. R. Pagenkopf, G. L. Rupp, Johnson, K. M., Chan, P. W. H., Gherini, S. A. and C. E. Chamberlin. (1985). Rates, constants, and kinetics formulations in surface water quality modeling. EPA/600/3-85/040. U.S. Environmental Protection Agency, Environmental Research Laboratory, Athens, GA. 455pp.
5.Brown, L. C., and Barnwell, T. O., Jr. (1987). The enhanced stream water quality models QUAL2E and QUAL2E-UNCAS. Documentation and users manual. USEPA Report EPA/600/3-87/007. Athens: Environmental Research Laboratory, Office of Research and Development.
6.Burn, D. H. (1989). Water-quality management through combined simulation–optimization approach. Journal of Environmental Engineering, 115, 1011–1024.
7.Chaudhury, R. R., Sobrinho, J. A. H., Wright, R. M., and Sreenivas, M. (1998). Dissolved oxygen modeling of the Blackstone River (Northeastern United States). Water Research, 32, 2400–2412.
8.Cho, J. H., Ahn, K. H., Chung, W. J., and Gwon, E. M. (2003).Waste load allocation for water quality management of a heavily polluted river using linear programming. Water Science and Technology, 48, 185–190.
9.Dobbins, W. E. (1964). BOD and oxygen relationships in streams. J.Sanit. Eng. Div., Am. Soc. Civ. Eng., 90(3), 53–78.
10.Kao, J. J., and Bau, S. F. (1996). Risk analysis for flow duration curve based seasonal discharge management programs. Water Research, 30, 1369–1376.
11.Lee, C. S., and Chang, S. P. (2005). Interactive fuzzy optimization for an economic and environmental balance in a river system. Water Research, 39, 221–231.
12.Lee, C. S., Chang, S. P., and Wang, C. M. (2002). Using hand-held ystem in water quality modeling of river. 8th Mainland China–Taiwan Environmental Protection Conference (pp. 741–746). Hsin-Chu.
13.Leopold, L. B., and Maddock, T. (1953). The hydraulic geometry of stream channels and some physiographic implications. U.S. Geological Survey Professional Paper 252. Washington, D.C.
14.Linfield, C. and Thomas, O. B. Jr., (1985). Computer Program Documentation for the Enhanced Stream Water Quality Model QUAL2E, Department of Civil Engineering Tufts University and Environmental Research Laboratory.
15.Masch, F. D. and Associates, (1971) “Tidal Hydrodynamic and Salinity Model for San Antonio and Matagorda Bays”, Texas, A Report to Texas Water Development Board, Austin, Texas.
16.Mujumdar, P. P., and Vemula, V. R. S. (2004). Fuzzy water load allocation model: Simulation–optimization approach. Journal of Computing in Civil Engineering, 18, 120–131.
17.Pardalos, P. M., and Resende, M. G. C. (Eds.). (2002). Handbook of applied optimization. New York: Oxford University Press.
18.Park, C. C. (1977). World-wide variations in hydraulic geometry exponents of stream channel: An analysis and some observations. Journal of Hydrology, 33, 133–146.
19..QUAL2E Windows Interface User’s Guide,United States Environmental Protection Agency,August 1995.
20.Schrange, L. (1991). User’s manual for linear, integer, and quadratic programming with LINDO, Release 5.0. San Francisco: Scientific Press.
21.TEPA (1993). Pollution control planning in the Pei-Kung River and Putzu River Basins–Putzu River. Taipei: Environmental Protection Administration of Taiwan.
22.TEPA (2002). The evaluation of water-quality improvement using surface aerators in the Putzu River. Taipei: Environmental Protection Administration of Taiwan.
23. U. S. Environmental Protection Agency, http://www.epa.gov/ .
24.U. S. Environmental Protection Agency, (1991). Guidance for Water-Quality-based Decisions: The TMDL Process: Office of Water Washington, DC 20460, EPA 440/4-91-001.
25.Yang, C. C., Chen, C. S. and Lee, C. S. (2011) Comprehensive River Water Quality Management by Simulation and Optimization Models. Environmental Model Asscess, 16, 283-294.
26.台灣省水利局 (1998)，朴子溪(含牛稠溪)治理規劃檢討報告。
27.台灣省自來水公司 (2000)，水質二年報。
28.行政院環保署(2001)，「朴子溪流域排水幹渠水質水量監測及特性分析」。
29.行政院環保署 (1993)，北港溪及朴子溪流域污染整治規劃。行政院環保署，台北市。
30.行政院環保署 (2002)，朴子溪流域水質改善曝氣設施設置評估規劃。行政院環保署，台北市。
31.行政院環保署(2003)，九十一年度“中港溪流域污染整治實施計劃”技術轉移及教育訓練會議，「QUAL2E 水質模式建立」。
32.高正忠 (1995)，國立交通大學環境系統分析講義。
33.國立台灣大學土木研究所 (1998)，高屏溪、東港溪污染量調查及整治規劃之高屏溪總量管制－ 期中報告，行政院環境保護署計劃，EPA-87-G105-03-03。
34.許銘熙、郭義雄、郭振泰、柳文成 (1998a)，淡水河感潮段垂直二維水理與水質動態傳輸模式(一)。行政院國科會專題計畫研究成果報告。
35.許銘熙、郭振泰、郭義雄、柳文成 (1998b)，淡水河系潮流、河口環流與鹽分佈之研究(一)。國立臺灣大學水工試驗所研究報告239號，台北市。
36.許銘熙、郭振泰、郭義雄、柳文成 (1999)，淡水河系潮流、河口環流與鹽分佈之研究(二)。國立臺灣大學水工試驗所研究報告273號，台北市。
37.溫清光 (1989)，高屏溪流域河川水質規劃—高屏溪水資源污染調查規劃。
38.劉思明 (1988)， QUAL2E 在河川涵容能力求法上之應用，國立成功大學環境工程研究所碩士論文。
39.嘉義縣環保局 (2001)，朴子溪污染來源申報資料。
40.嘉義縣環境保護局 (2001)，朴子溪流域污染整治實施計畫，嘉義縣環保局。
41.嘉義縣環境保護局 (2002)，嘉義縣河川水體水質及地下水水質監測工作計畫，嘉義縣環保局。
42.嘉義縣環保局 (2004)，嘉義縣河川水質改善評估與規劃工作計畫，嘉義縣環保局。
43.蕭培元 (1999)，淡水河流域污水廠最佳配置研究，淡江大學水資源及環境工程學系碩士論文。