河系

本研究預計以polyline 之水系 GIS 圖資，計算單位網格內之河系總長 度。

△地下水觀測井

圖4.2-15、濁水溪流域河系圖 4.3 因子分數圖層處理

4.3.1 土地利用

土地利用原始分類高達 103 種，本研究將其歸納為 7 種土地利用類 別。圖4.3-1 是由圖 4.2-1 經過網格切割後，依照各種計劃土地分類根據表 3.4-2 給予分數，再將該分類分數所佔網格面積比例加權計算各網格分數，

最後再搭配表3.4-1 以補注潛勢因子間權重加權計算所得。

圖4.3-1、加權後土地利用潛勢分數分布圖 4.3.2 表層土壤種類

圖4.3-2 是表層 0~150cm 土壤經網格切割後依照各種土壤根據表 3.4-2 給予四種深度分層土壤種類分數並依該土壤種類分類所佔網格面積比例 加權計算各網格分數，再將各網格四層分數平均，接著再搭配表 3.4-1 以 補注潛勢因子間權重加權計算所得之表層土壤種類補注潛勢分數分布 圖，對照圖4.2-2 可看出粒徑較大的土壤分布區域其潛勢分數較高；反之，

的影響，故無法像較深層之土壤分布可以看出地下水儲存之區域。表4.3-2 為本研究之分類代碼與原始圖層分類之對照，其中因為原始圖層在濁水溪 流域無土壤資料，導致分數過低，經評估後決定將該區域直接給定最高分。

表4.3-2、表層土壤種類計劃分類表

分類代碼 原始分類 細項

1 粗砂土 粗砂土，砂土

2 細砂土 細砂土，壤質粗砂土，壤質砂土 3 壤質細砂土 壤質細砂土，粗砂質壤土，砂質壤土

，細砂質壤土

4 極細砂土 極細砂土，壤質極細砂土，極細砂質壤土

5 坋土 坋土，坋質壤土

6 壤土 壤土

7 砂質粘壤土 砂質粘壤土

8 粘壤土 粘壤土，坋質粘壤土 9 坋質粘土 坋質粘土

10 粘土 粘土

11 石礫 石礫

圖4.3-2、加權後表層土壤種類潛勢分數分布圖

4.3.3 降雨與地下水位變化相關性

圖 4.3-3 是以距離倒數平方法內插出 1 公里見方網格中心座標之相關 性值，依照表3.4-2 分類後再搭配表 3.4-1 以補注潛勢因子間權重加權計算 所得。

圖4.3-3、加權後降雨與地下水位變化相關性潛勢分數分布圖

4.3.4 平均年降雨量

圖 4.3-4 是以距離倒數平方法內插出 1 公里見方網格中心座標之 1993~2008 平均年雨量值後，再依照表 3.4-2 分類後再搭配表 3.4-1 以補注 潛勢因子間權重加權計算所得。由第 3.1.2 節之平均降雨量等值圖可知靠 近扇頂海拔較高處雨量豐沛，降雨量由東往西遞減，因此本因子影響地下 水補注潛勢之分數亦是由東往西遞減。

圖 4.3-4、1993~2008 年加權後平均年降雨量潛勢分數分布圖

4.3.5 單位蓄水量變化

圖 4.3-5 是以距離倒數平方法將點位之單位蓄水量變化資料內插出 1 公里見方網格中心座標之蓄水量變化後，再依照表 3.4-2 分類後再搭配表 3.4-1 以補注潛勢因子間權重加權計算所得。發現右下角靠近斗六丘陵處 降雨量大且降雨與地下水位變化相關性高，因此蓄水量變化較大。

圖4.3-5、加權後單位蓄水量變化潛勢分數分布圖

4.3.6 透水係數

圖 4.3-6 是以距離倒數平方法內插出 1 公里見方網格中心座標之含水 層一(F1)透水係數後，再依照表 3.4-2 分類給分後，再搭配表 3.4-1 以補注 潛勢因子間權重加權計算所得的透水係數潛勢分數分布圖。由圖 4.3-7 可 知扇頂區域之透水係數相對較大。

圖4.3-6、濁水溪沖積扇含水層一加權後透水係數潛勢分數分布圖

4.3.7 河系密度

圖 4.3-7 是以 SuperGIS 計算單位網格內河系之總長度，依照表 3.4-2 分類給分後，再搭配表3.4-1 以補注潛勢因子間權重加權計算所得。

圖4.3-7、加權後河系密度潛勢分數分布圖

4.4 補注潛勢評估

在前述4.3 節完成濁水溪沖積扇各補注潛勢因子之分數計算，並且將 各因子分數分別根據因子間權重百分比加權，最後加總即可得到整個區域 內各網格之補注潛勢總分，再搭配表 4.4-1 對各網格進行補注潛勢總分分 級，分級成果如圖 4.4-1 所示。將圖 4.4-1 套疊鄉鎮邊界圖後，由圖 4.4-2 可看出 50 分以上中高補注潛勢分布的鄉鎮有彰化縣的社頭鄉、田中鎮、

二水鄉、北斗鎮、溪州鄉、和雲林縣的莿桐鄉、林內鄉、斗六市、西螺鎮、

古坑鄉以及南投縣的竹山鎮、名間鄉等地區。

本研究將補注潛勢總分 50 分以上的網格劃分出來，即可概略推估補 注區界線。圖4.4-3 為本研究以 50 分以上為補注區界線與工研院資料推斷 之界線對比圖，黃線為本研究推估之補注區界線，綠線為以工研院(1996) 地電阻測深法及觀測井資料所推測之補注區界線，由圖可發現兩界線相當 接近。

本研究推估之補注區約涵蓋500 個網格(面積約 500 平方公里)，大致 占濁水溪沖積扇面積之23%，而以工研院資料推估之補注區約涵蓋 450 個 網格(面積約 450 平方公里) ，大致占濁水溪沖積扇面積之 21%，。圖 4.4-4 為兩者補注區範圍對比圖，由圖可看出本研究涵蓋範圍在上方較廣而下方 較狹小。主要差異為本研究上方有涵蓋到社頭鄉，但工研院並無；工研院

研究涵蓋虎尾鎮三分之一、斗南鎮六分之一以及完整涵蓋古坑鄉。雖然兩 者有些許差異，但整體來說相當接近，故驗證了本研究方法之可行性。

表4.4-1、地下水補注潛勢分數分級表

地下水補注潛勢分數 分級

大於80 很高(very high) 介於70 到 80 之間 高(high)

介於60 到 70 之間 中高(high moderate) 介於50 到 60 之間 中(moderate)

介於40 到 50 之間 中低(low moderate) 介於30 到 40 之間 低(low)

小於30 很低(very low)

圖4.4-1、地下水補注潛勢總分分級圖

圖4.4-2、地下水補注潛勢總分分級以及鄉鎮邊界圖

圖4.4-3、50 分以上補注區界線與地電阻探測法研判之界線對比圖

圖4.4-4、本研究與地電阻探測法研判之補注區範圍對比圖

第五章 結論與建議

5.1 結論

1. 本研究整合 DRASTIC 之因子分析方法以及 Shaban 等人之因子間因果 關係分析方式，再以地理資訊系統(GIS)為圖層處理平台，發展出一系 統化的地下水補注潛勢分析流程，並將其應用於濁水溪沖積扇補注潛 勢評估。

2. 本研究目前針對濁水溪沖積扇考量之地下水補注潛勢因子有「土地利 用」、「表層土壤種類」、「降雨與地下水位變化相關性」、「平均 年降雨量」、「單位蓄水量變化」、「透水係數」以及「河系密度」

等七項補注潛勢因子。

3. 分析結果顯示「降雨與地下水位變化相關性」、「平均年降雨量」、

「單位蓄水量變化」、「透水係數」等因子圖層皆反應出高補注區位 於扇頂區域，並且有由扇頂向扇尾遞減之趨勢。

4. 在同樣補注水量之下，非受壓含水層由於土層蓄水能力較受壓含水層 高，即比出水量大於儲水係數，前者之水位變化將較後者小，因此為 反映系統所受外來補注量的大小，本研究以「單位蓄水量變化」代表 之，分析結果發現其值較高者位於扇頂區域以及斗六丘陵地帶。

5. 若以補注潛勢總分 50 分為界，推估濁水溪沖積扇扇頂補注區之範圍，

莿桐、林內及南投縣竹山等鄉鎮，與工研院能源與資源研究所於 1996

參考文獻

1. Bou Kheir R., A. Shaban, M-C. Girard, M. Khawlie, and C.

Abdallah, “Caract´erisation morpho-p´edologique des zones karstiques du Libansensibilit´e des sols `a l’´erosion hydrique”, S´echeresse 14:4, 2003.

2. Bromley, J., W. M. Edmunds, E. Fellman, J. Brouwer, S. R. Gaze, J.

Sudlow,and J. D. Taupin, “Estimation of rainfall inputs and direct recharge to the deep unsaturated zone of southern Niger using the chloride profile method”, Journal of Hydrology, Vol. 188-189, P.139-154, 1997.

3. Braun, Gary M., Norman S. Levine, Sheila J. Roberts and Arthur N. Samel,

“A Geographic Information Systems Methodology for the Identification of Groundwater Recharge Areas in Waukesha County, Wisconsin”, Environmental and Engineering Geoscience; v. 9; no. 3; p. 267-278, August 2003.

4. Bhuiyan, C., Ramesh P. Singh, W. A. Flugel, “Modelling of ground water recharge-potential in the hard-rock Aravalli terrain India a GIS approach”, Environ Earth Sci 59:929–938, 2009.

5. Edet, A. E., C. S. Okereke, S. C. Teme, and E. O. Esu, “Application of remotesensing data to groundwater exploration: a case study of the Cross River State,southeastern Nigeria”, Hydrogeology Journal, 6(3): 394-404, 1998.

6. Greenbaum, D., “Review of remote sensing applications to groundwater exploration in basement and regolith”, British Geological Survey Report OD 85(8):36, 1985.

7. Jaiswal, R. K., S. Mukherjee, J. Krishnamurthy, and R. Saxena, “Role of remote sensing and GIS techniques for generation of groundwater prospect zones towards rural development-an approach”, International Journal of Remote Sensing, 24(5): 993-1008, 2003.

8. Jasrotia, A.S.,R. Kumar,A. K. Saraf, “Delineation of groundwater recharge sites using integrated remote sensing and GIS in Jammu district, India”, International Journal of Remote Sensing, Volume 28 , Issue 22 , Pages 5019-5036,November 2007.

9. Krishnamurthy, J., N. Venkatesa Kumar, V. Jayaraman, and M. Manivel,

“An approach to demarcate groundwater potential zones through remote sensing and geographic information system”, International Journal of Remote Sensing, 17: 1867-1884, 1996.

10. Lebuc, C., G. Favreau, and P. Schroeter, “Long-term rise in a Sahelian water table: the Continental Terminal in southwest Niger”, Journal of Hydrology, Vol. 243, P.43-54, 2001.

11. Murthy, K. S. R., “Groundwater potential in a semi-arid region of Andhra

groundwater recharge zone in hard rock terrain”, Current Science, Vol. 87, 5, 10, 2004.

13. Mukherjee, S., “Targeting saline aquifer by remote sensing and geophysical methods in a part of Hamirpur-Kanpur, India”, Hydrology Journal, 19, 53-64, 1996.

14. National Water Well Association, “DRASTIC: A Standardized System for Evaluating Ground Water Pollution Potential Using Hydrogeologic Settings”, USEPA, May 1985.

15. Shahid, S., S. K. Nath, and J. Roy, “Groundwater potential modeling in a soft rock area using a GIS”, International Journal of Remote Sensing, 21:

1919-1924, 2000.

16. Shaban, A., M. Khawlie, R. Bou Kheir, C. A, “Assessment of road instability along a typical mountainous road using GIS and aerial photos, Lebanon – eastern Mediterranean”, Bull Eng Geol Env, 60:93-101, 2001.

17. Sener, E., A. Davraz, and M. Ozcelik, “An integration of GIS and remote sensing in groundwater investigations: A case study in Burdur, Turkey”, Hydrogeology Journal, 13: 826-834, 2005.

18. Shaban, A., M. Khawlie, and C. Abdallah, “Use of remote sensing and GIS to determine recharge potential zone: the case of Occidental Lebanon”, Hydrogeology Journal, 14: 433-443, 2006.

19. Yeh, H.F., Lee, C.H., Hsu, K.C., and Chang, P.H., “GIS for the assessment of the groundwater recharge potential zone”, Environ. Geol., Vol. 58, 185-195, 2008.

20. Kang-Tsung Chang, INTRODUCTION TO GEOGRAPHIC INFORMATION SYSTEMS, McGraw-Hill Companies, 2008.

21. 田巧玲、張炎銘、林維侃、蔡明坤，「濁水溪沖積扇之地下水水文概況」，

頁，2003。

39. 環保署全國環境水質監測網，http://wqshow.epa.gov.tw/，2009。

40. 朱子豪，「美國民間 GIS 發展動態」，國土資訊系統通訊，第十三期，

程系碩士論文，2008。

附錄 A 地理資訊系統簡介

地理資訊系統(Geographic Information System，簡稱GIS)狹義來說是電 腦系統(Computer System)，但廣義來說資訊系統(Information System)並不 一定要建置在計算機上。例如早期圖書館的圖書查詢櫃，可以按照編目卡 或作者卡來查詢書籍。

地理資訊系統(GIS)和資訊系統(Information System)主要的功能都包 含了取得(Capturing)、儲存(Storing)、搜尋(Querying)、分析(Analyzing)以 及展示(Displaying)的功能。而GIS的數據(Geographically Referenced Data) 包含了位置(Location)和屬性(Characteristic)，所以地理資訊系統可以對空 間進行分析，這和其他單純的資訊系統是不同的。本節將對地理資訊系統 的發展、地理資訊系統的桌上型軟體、以及本研究圖層資料標準化進行詳 細的說明。

 2.2.1地理資訊系統發展史

地理資訊系統的存在與發展已歷經30餘年。用戶的需要、技術的進 步、應用方法論的提高，以及有關組織機構的建立等因素，深深地影響著 地理資訊系統的發展。綜觀GIS發展，尤其是北美地區的實際情況，大致 可將地理資訊系統發展分為以下四個階段：

第一階段：

處理人口統計局資料（如美國人口調查局建立的DIME）、資源普查資料

（如加拿大統計局的GRDSR）等。許多大學研製了一些基於柵格系統的 套裝軟體，如哈佛的SYMAP、馬里蘭大學的MANS等。綜合來看，初期 地理資訊系統發展的動力來自於諸多方面，如學術探討、新技術的應用、

大量空間資料處理的生產需求等。對於這個時期地理資訊系統的發展來

大量空間資料處理的生產需求等。對於這個時期地理資訊系統的發展來