6.1 結論
1. 本研究整合條件模擬與分佈型水平衡模式，並應用蒙地卡羅模擬 的理念，藉由模擬濁水溪沖積扇 1000 組合透水係數所對應之 1000 組抽水量，完成對濁水溪沖積扇抽水量推估不確定性之研究。分 析成果顯示平均而言，年補注量第一含水層(F1)為 5.15 億噸、年 抽水量第二及第三含水層(F2 與 F3)分別為 4.56 億噸及 2 億噸，
而標準差各層分別為 0.64、0.61 與 0.43 億噸，此結果顯示透水係 數之不確定性對抽水量推估確實在各層造成了 12%、13%及 22%
的影響。
2. 條件模擬部份，本研究應用 GSLIB 之 SASIM 程式針對濁水溪沖 積扇各分層分別產生出 100 組透水係數之網格資料，並經驗證顯 示，各分層之各組合成網格資料的半變異元值皆符合各分層觀測 資料之半變異元，顯示條件模擬所合成之網格資料與觀測資料之 半變異元一致。
3. 前述條件模擬所得之各分層合成透水係數資料顯示，各分層透水 係數之平均值分別為第一層 49.75、第二層 47.48 與第三層
48.81(m/day)，及標準差分別為 36.41、16.68 及 25.97(m/day)。
4. 傳統模擬模式（如 MODFLOW）為給定參數及其他抽水或補注量
62
值而計算出水位值，為符合本研究蒙地卡羅模擬，需模擬多組透 水係數所對應之多組抽水量的需求，本研究另行以水平衡條件為 基礎，發展一在給定水文地質參數與水位值下，計算補注量之模 擬模式，應用此模式可對多組透水係數場快速計算出其相對之多 組抽水量。
5. 本研究將抽水量之標準差求得其空間分佈圖，可發現在 F1 標準 差分佈顯示值較大者集中於扇頂，F2 之標準差分佈則較鬆散並無 明顯集中的情況，F3 標準差分佈則顯示較大值集中於扇央。
6.2 建議
1. 本研究目前仍為穩態模擬，建議將來可進行暫態模擬之研究，除 了抽水量空間及總量之不確定性外，可進一步評估抽水量時間分 布上之不確定。
2. 本研究實際案例為濁水溪沖積扇，未來可以將本研究應用於不同 的地區並比較其相關的半變異元。
63
參考文獻
1. Avoufirassi , M. and M. A. Marino, ”Cokrining of Aquifer Transmissivity from Field Measurements of Transmissivity and SpecificCapacity”, Mathematical Gelolgy, Vol.16, No.1, pp.19-34.
（1984）
2. Deutsch, C. V. & Journel , A. G. ,GSLIB: Geostatistical Software Library and User’s Guide,Oxford University Press , New York
3. 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.
4. Pete Bettinger, David Graetz, Kevin Boston, John Sessions and WoodamChung, 2002. Eight Heuristic Planning Techniques Applied to Three Increasingly Difficult Wildlife Planning Problems. Fennica, 36(2),561-584.
5. Warrick, A. W., D. E. Myers and D. R. Nielsen,”
GeostatisticsMethods Applied to Soil Science”, in A. Klute(ed.) Methoes of SoilAnalysis,part1.P53-80（1986）
6. Y.P. Lin,Yih-Chi Tan,Shahrokh Rouhani,”Identifying spatial characteristicsof transmissivity using simulatedannealing and kriging methods”,Environmental Geology,Vol41,pp.200-208(2001)
7. 中興工程顧問公司，濁水溪沖積扇地下水可開發潛能評估。台灣 省政府水利處委託計畫報告，第 5-1 至 5-12 頁，第 6-1 至 6-28 頁，第 7-1 至 7-16 頁，1997。
64
8. 中興工程顧問公司，濁水溪沖積扇地下水人工補注規劃。台灣省 政府水利處委託計畫報告，第 4-6 至 4-18 頁，1998。
9. 巨廷工程顧問股份有限公司、國立交通大學，多元化水資源經營 管理方案－地下水補注分析。經濟部水利署委託計畫成果報告第 2-1 至 2-24 頁，第 4-1 至 4-18 頁，2005。
10. 台灣大學水工試驗所，濁水溪沖積扇扇央扇尾平原區地表垂直向 補注量評估。經濟部水資源局委託計畫報告，1998。
11. 台灣大學土木工程研究所，濁水溪沖積扇扇頂平原地區地表垂向 補注量估計，經濟部委託計畫報告，1997。
12. 能邦科技股顧問股份有限公司，台灣地區地下水補注量估算。經 濟部水資源局委託計畫報告，第 4-1 至 4-70 頁，第 5-1 至 5-140 頁，2000。
13. 農業工程研究中心，7C 濁水溪平原地下水數學模式之評估與應用。
台灣省水利局委託計畫報告，1989。
14. 劉聰桂，濁水溪沖積扇之地下水資源－碳十四與氚定年/示踨研究。
濁水溪沖積扇地下水及水文地質研討會論文集，第 145-164 頁，
1996。
15. 經濟部中央地質調查所，台灣地區地下水觀測網第一期計畫，濁 水溪沖積扇水文地質調查研究總報告，1999。
16. 江崇榮、黃智昭、陳瑞娥，以地下水歷線分析法評估濁水溪沖積 扇之地下水收支，經濟部中央地質調查所彙刊第十九號，第 61-90 頁，2006。
65
17. 經濟部中央地質調查所，台灣地區地下水區水文地質調查及地下 水資源評估，地下水補注潛勢評估與地下水模式建置，2009。
18. 經濟部水利署，地下水資源整體營運規劃與綜合評估，2002。
19. 王雲直，應用專家系統於地下水模式參數檢定之研究─以濁水溪 沖積扇為例，國立交通大學，碩士論文，民國 99 年
20. 張弼舜，應用專家系統於穩健型地下水參數檢定模式之發展，國 立交通大學，碩士論文，民國 100 年
21. 楊深惠，應用衛星影像辨識與河道水里演算於濁水溪沖積扇地下 水數值模擬，國立交通大學，碩士論文，民國 101 年
22. 潘禎哲，改良式降低變異數分析於地下水觀測井網規劃之應用，
國立交通大學，碩士論文，民國 83 年
23. 許文鴻，因子克利金應用於地閜水質監測井網設計之研究，國立 交通大學，碩士論文，民國 93 年
24. 邱怡芬，聯克利金法於地下水觀測井網規劃之應用，國立交通大 學，碩士論文，民國 85 年
25. 陳敏華，應用模擬退火法及景觀指數於及水區土地利用最佳化規 劃，台灣大學，碩士論文，民國 96 年
26. 蕭朱杏、莊愷瑋，地理統計於醫學與環境之應用，地理統計在農 業和環境科學之應用研討會論文集，pp.93-112.，民國 90 年
66
27. 黃信彰、徐國錦，以指標地質統計法建構水文地質模式之研究，
成功大學，碩士論文，民國 91 年
28. 番成正，優勢水流對坡地安定影響之研究 - 以草嶺為例，成功大學，
碩士論文，民國 92 年
29. 劉原，地質統計方法對水文地質場模擬影響之研究，成功大學，
碩士論文，民國 98 年