第七章 結論與建議
7.2 建議
1. 由於專家系統的知識庫可持續累積，故未來遇到特殊邊界條件或特殊地 質情況，可於知識庫中增加新的規則，將能適用於這些地區，而不需重 新建置新的系統。
2. 本研究為穩態模式之參數檢定，未來可增加新規則，使之能應用在暫態 模式之參數檢定。
3. 本研究於實際案例之概念分層第一層，僅能檢定淨補注量，未來可由水 位歷線法等方法，算出實際補注量並建置於數值模式裡，將可檢定出概 念分層第一層之抽水量，及系統總抽水量。
4. 本研究假設待檢定參數完全未知，因此位於同一層之分區，初始猜值均 設為相同，如此亦能檢定出合理之參數。未來在高維度模式中，可將本 檢定結果，當成優選法之初始猜值，以得出更精確之檢定結果。
5. 本研究 K 值設定與陳韋圻(2008)之研究一致，皆使用參數分區(Zonation) 法，降低參數維度，因此本研究之專家系統未來可直接整合陳韋圻(2008) 之參數檢定規則，擴充檢定K 值之能力。
參考文獻
1. K.W. Chau, “Intelligent manipulation of calibration parameters in numerical modeling”, Advances in Environmental Research, 8, p. 467–476, 2004.
2. K.W. Chau, “Selection and calibration of numerical modeling in flow and water quality”, Environmental Modeling and Assessment, 9, p. 169–178, 2004.
3. K.W. Chau, “A review on integration of artificial intelligence into water quality modeling”, Marine Pollution Bulletin, 52, p. 726–733, 2006.
4. K.W. Chau, “A review on the integration of artificial intelligence into coastal modeling”, Journal of Environmental Management, 80, p. 47–57, 2006.
5. Sang Min Kim, Brian L. Benham, Kevin M. Brannan, Rebecca W. Zeckoski, John Doherty, 2007, , Comparision of hydrologic calibration of HSPF using automatic and manual methods, WATER RESOURCES RESEARCH, VOL.
43, W01402, doi:10.1029/2006WR004883.
6. Chau, K.W., 2003. Manipulation of numerical coastal flow and water quality models. Environmental Modelling and Software 18 (2), 99–108.
7. Chau, K.W., Albermani, F., 2002. Expert system application on preliminary design of liquid retaining structures. Expert Systems with Applications 22 (2), 169–178.
8. Chau, K.W., Albermani, F., 2003. Knowledge-based system on optimum design of liquid retaining structures with genetic algorithms. Journal of Structural Engineering, ASCE 129 (10), 1312–1321.
9. Chau, K.W., Chen, W., 2001. An example of expert system on numerical modelling system in coastal processes. Advances in Engineering Software 32 (9), 695–703.
10. Henrik Madsen, Geoffery Wilson, Hans Christian Ammentorp, 2002, Comparison of different automated strategies for calibration of rainfall run-off models.
11. Abbott, M.B., 1991. Hydroinformatics: Information Technology and the Aquatic Environment. Avebury Technical, Aldershot
12. Madsen, H., 2003. Parameter estimation in distributed hydrological catchment modelling using automatic calibration with multiple objectives.
13. Madsen, H., 2000. Automatic calibration of a conceptual rainfall-runoff model using multiple objectives. J. Hydrol. 235, pp. 276–288.
14. Hill, M.C., 1992. A computer program (MODFLOWP) for estimating parameters of a transient, three-dimensional, ground-water flow model using nonlinear regression. US Geological Survey, Open File Report 91-484, pp. 3–4 and 15.
15. Schwartz, F.W., and Zhang, H., 2003, Fundamentals of ground water, John Wiley & Sons, New York.
16. 中興工程顧問公司，濁水溪沖積扇地下水可開發潛能評估。台灣省政府 水利處委託計畫報告，第 5-1 至 5-12 頁，第 6-1 至 6-28 頁，第 7-1 至 7-16 頁，1997。
17. 中興工程顧問公司，濁水溪沖積扇地下水人工補注規劃。台灣省政府水 利處委託計畫報告，第 4-6 至 4-18 頁，1998。
18. 巨廷工程顧問股份有限公司、國立交通大學，多元化水資源經營管理方 案－地下水補注分析。經濟部水利署委託計畫成果報告第 2-1 至 2-24 頁，第 4-1 至 4-18 頁，2005。
19. 台灣大學水工試驗所，濁水溪沖積扇扇央扇尾平原區地表垂直向補注量 評估。經濟部水資源局委託計畫報告，1998。
20. 台灣大學土木工程研究所，濁水溪沖積扇扇頂平原地區地表垂向補注量 估計，經濟部委託計畫報告，1997。
21. 能邦科技股顧問股份有限公司，台灣地區地下水補注量估算。經濟部水 資源局委託計畫報告，第 4-1 至 4-70 頁，第 5-1 至 5-140 頁，2000。
22. 農業工程研究中心，7C 濁水溪平原地下水數學模式之評估與應用。台 灣省水利局委託計畫報告，1989。
23. 劉聰桂，濁水溪沖積扇之地下水資源－碳十四與氚定年/示踨研究。濁 水溪沖積扇地下水及水文地質研討會論文集，第 145-164 頁，1996。
24. 經濟部中央地質調查所，台灣地區地下水觀測網第一期計畫，濁水溪沖 積扇水文地質調查研究總報告，1999。
25. 江崇榮、黃智昭、陳瑞娥，以地下水歷線分析法評估濁水溪沖積扇之地 下水收支，經濟部中央地質調查所彙刊第十九號，第 61-90 頁，2006。
26. 經濟部中央地質調查所，台灣地區地下水區水文地質調查及地下水資源 評估，地下水補注潛勢評估與地下水模式建置，2009。
27. 經濟部水利署，地下水資源整體營運規劃與綜合評估，2002。
28. 曾憲雄，人工智慧與專家系統 : 理論.實務.應用，初版，旗標出版社，
台北市，2005。
29. 陳韋圻，應用專家系統於地下水模式自動化參數檢定之研究，國立交通 大學，碩士論文，民國 97 年。
附錄 A MODFLOW 簡介
MODFLOW為美國地質調查局(U.S.G.S.)發展之程式。該程式可解二維 及三維之地下水流問題，含水層之種類可為自由、受壓﹑半受壓含水層，
依地質特性分類可為均質、非均質及等向性、非等向性含水層。MODFLOW 程式乃利用有限差分法(Block Centered Finite Difference Approach)解水流控 制方程式，計算機數值求解方法乃採用兩種疊代技巧強制隱式法(SIP)及鬆 Conductivity)(LT−1)
h：管壓水頭(Potentiometric Head)（L)
W ： 單 位 體 積 的 體 積 流 率 (Volumetric Flux) ， 代 表 源 匯 項 (Sources/Sinks)(T−1)
Ss：孔隙介質的比儲水量(Specific Storage)（L−1） t：表時間(T )
上述式（附A.1）若結合了含水層系統邊界情況、起始條件等資訊，可
cells((i-1,j,k),(i+1,j,k),(i,j-1,k),(i,j+1,k),(i,j,k-1),(i,j,k+1))。如圖附A.1所示：
假設在列（row) 方向的 cell(j,j-1,k) 流進 cell(i,j,k) 的流量為：
( )
現在若在多加以考慮源、匯（Sources、Sinks)，則連續性方程式變為：
( ) ( ) ( )
MODFLOW中所採用的為後向差分(Backward Difference)，所以對於 cell(i,j,k)來說，若以 tm 和 tm−1 之間來代表 ∆t，則：
而附A.7式即是MODFLOW程式所解之差分式。
圖附A.1 cell(i,j-1,k)進入 cell 之地下水流