為管理雲林地區地層下陷易淹區之地下水資源，首先利用MODFLOW軟體對雲林地區地下水流進行三維數值模擬，以求得符合現地狀況之水文地質架構及水力參數。利用建構完成之三維數值模式，配合MODOFC軟體對模式產生響應矩陣並進行線性規劃，最後求得限制條件與目標函數下之最佳解，由此構成地下水資源管理模式，以達區域地下水管理之目的。
本研究在進行地下水流模擬中，水平透水係數之大小的調整對於地下水流場改變量較大，其次為補注量。在模式率定結果中，乾濕季水位變化、未知之抽水行為及邊界效應影響，扇頂、扇尾及邊界附近之井位率定結果水位偏差較大。
在管理模式中，由最佳化管理抽水方案之管理結果可得知扇頂、扇央及扇尾之最佳抽水地點分別為古坑鄉、西螺鎮及麥寮鄉。以2007年一年期管理之最佳抽水方案之最佳抽水地點為位於扇央區之古坑鄉，總抽水量為每年5134.9萬噸，扇頂區地下水位下降量為1.9至6.3公尺。
由最佳化管理補注方案之管理結果可得知扇頂、扇央及扇尾之最佳補注地點分別為林內鄉、虎尾鎮及台西鄉。以2007年一年期管理之最佳補注地點為位於扇尾區之台西鄉，總補注量為每年1351.6萬噸，地下水位上升量為0.3至1.2公尺。

The groundwater resources management in Yunlin land subsidence area has become a important issue. In this study, two stages of the study are performed. In the first stage, a regional hydro-geological model in Yunlin was constructed. Groundwater level for Yunlin was modelled using the numerical model MODFLOW. From the calibration result, the aquifer hydraulic parameters of the system can be estimated. In the second stage, the numerical model is combined with the optimal model, MODOFC, which can generate response matrix of allowable pumping or recharging water and solve in Linear Programming to obtain the best value for groundwater resources management.
By calibrating with groundwater level, we can find out that the parameter of conductivity is the most significant factor in model and recharge is less significant. The calibration result shows that unknown pumping, seasonal groundwater level variation and the boundary influence cause model deviation in top-fan, distal-fan and nearby boundary area.
In the pumping scenario of management model, the optimal pumping location at top-fan , mid-fan and distal-fan is Gukeng Township, Siluo Township and Mailiao Township, respectively. In 2007, from 1 through 12 of a period of one month, the best pumping scenario allows to an increasing 51.3 million m3 of groundwater pumping. This optimal solution is evaluated to lower the groundwater level about 1.9 m to 6.3 m at top-fan.
In the recharging scenario of management model, the optimal recharging location at top-fan , mid-fan and distal-fan is Linnei Township, Huwei Township and Taisi Township, respectively. In 2007, from 1 through 12 of a period of one month, the best recharging scenario allows to an increasing 13.5 million m3 of groundwater water recharging. Finally, this optimal solution is evaluated to raise the groundwater level about 0.3 m to 1.2 m at distal-fan.

1. Aguado, E. and I. Remson, “Groundwater hydraulics in aquifer management”, Journal of the Hydraulics Division, 100(1)：103-118, 1974.
2. Alley, W. M., E. Aguado, and I. Remson, “Aquifer management under transient and steady-state condition”, Water Resources Bulletin, 12(5)：963-972, 1976.
3. Atwood, D.F., and S. M. Gorelick, “Hydraulic gradient control for groundwater contaminant removal”, Journal of Hydrology, 76(1-2)：85-106, 1985.
4. Cheng, Y., C. H. Lee, Y. C. Tan and H. F. Yeh, “An optimal water allocation for an irrigation district in Pingtung country, Taiwan”, Irrig. and Drain., DOI: 10.1002/ird, 2008.
5. Deninger, R.A., “System analysis of water supply systems”, Water Resources Bulletin, 6(4)：573-579, 1970.
6. Don, N. C., N. T. M. Hang, H. Araki, H. Yamanishi, and K. Koga, “Groundwater resources management under environmental constraints in Shiroishi of Saga plain, Japan”, Environ Geol, 49: 601-609 , 2006.
7. Gong, H. , M. Li and X. Hu, “Management of groundwater in Zhengzhou City, China”, Water Research, 34(1)：57-62, 2000.
8. Gorelick, S. M., “A review of distributed parameter groundwater management modeling methods”, Water Resources Research, 19(2)：305-319, 1983.
9. Gorelick, S. M., C. I. Voss, P.E. Gill, W. Murrary, M. A. Saunnders and M.H. Wright, “Aquifer reclamation design : The use of contaminant transport simulation coupled with nonlinear programming”, Water Resources Research, 20(4)：415-427, 1984.
10. Lemoine, P. H., E. G. Reichard, and I. Remson, “An efficient response matrix method for coupling a groundwater simulator and a regional agricultural management model”, Water Resources Bulletin, 22(3)：417-423, 1986.
11. McDonald, M. G., and A. W. Harbaugh, “A modular three-dimensional finite-difference ground-water glow model, U.S. Geological Survey”, Virginia, 1988.
12. Peralta, R. C., H. Azarmnia, and S. Takahashi, “Embedding and response matrix techniques for maximizing steady-state ground-water extraction: computational comparsion”, Ground Water, 29(3)：357-364, 1991.
13. Psilovikos, A. A., “Optimization models in groundwater management, based on linear and mixed integer programming. An application to a Greek hydrogeological basin”, Physics and Chemistry of the Earth, Part B: Hydrology, Oceans and Atmosphere, 24(1-2)：139-144, 1999.
14. Takahashi, S., and R. C. Peralta, “Optimal perennial yield planning for complex nonlinear aquifers: methods and examples”, Advances in Water Resources, 18(1)：49-62, 1995.
15. Ting, C. S., “Groundwater Resources Evaluation and Management Studies for Pingtung Plain, Taiwan, phd theisis Free University, Amsterdam, ISBN-90- 9008794x, 1997.
16. Wang, C.H., C.H. Kuo, T.R. Peng, W.F. Chen, T.K. Liu, C.J. Chiang, W.C. Liu and J.J. Hung, “Isotope characteristics of Taiwan groundwater”, Western Pacific Earth Sciences, 1(4)：415-428, 2001.
17. Wattenbarger, R. A., “Maximizing seasonal withdrawals from gas storage reservoirs”, Journal of Petroleum Technology, 22(8)：948-998, 1970.
18. Wills, R. L., “A planning model for the management of groundwater quality”, Water Resources Research, 15(6)：1305-1312, 1979.
19. 丁澈士，「屏東平原地下水管理最佳化模式」，地下水資源及水質保護研討會，第265-284頁，1994。
20. 丁澈士、黃信恩，「屏東平原地下水人工補注水資源優化管理之可行性研究-以林邊溪流域為例」，臺灣水利，第51卷，第2期，第54-61頁，2003。
21. 中央氣象局全球資訊網。網址：http://www.cwb.gov.tw/V6/index.htm
22. 行政院環境保護署，「環境白皮書」，行政院環境保護署，1998。
23. 何麗如、汪中和、陳鎮東，「北港地區地下水之地球化學初探」，中華民國環境保護學會會誌，第13卷，第2期，第38-59頁，1990。
24. 吳漢雄、鄧聚龍、溫坤禮，「灰色分析入門」，高立圖書有限公司，1996。
25. 李友平、張國強、劉萬里，「台灣地區之地下水資源現況」，第四屆地下水資源及水質保護研討會論文集，第201-210頁，2001。
26. 李振誥，陳尉平，李如晃，「應用基流資料估計法推估台灣地下水補注量」，台灣水利，第50卷，第1期，第69-80頁，2002。
27. 李振誥、徐國錦，「地下水補注與管理措施」，工程，第80卷，第6期，2007。
28. 李振誥、許清荃，「地下水資源調配與管理-以濁水溪沖積扇為案例」，水資源管理，第8卷，第24-31頁，2000。
29. 李振誥、許清荃、林俶寬，「濁水溪沖積扇多層地下水調配與管理之研究」，臺灣水利，第48卷，第4期，第41-52頁，2000。
30. 沈向白，「濁水溪沖積扇地區地下水資源調查(含彰化縣、雲林縣)」，台灣省水利局，1991。
31. 季月華、朱國榮、江思珉，「地下水管理模型軟件GWM簡介及算例」，勘察科學技術，第1期，2007年。
32. 徐享崑、劉豐壽、鄭昌奇，「臺灣地區地層下陷之現況，成因與對策」，臺灣水利，第43卷，第3期，第19-29頁，1995。
33. 張良正、龔誠山，「區域性地下水觀測站網檢討八十四年度報告，區域性地下水觀測站網檢討(I)－濁水溪沖積扇觀測站井佈置檢討」，經濟部水利司，1995。
34. 陳忠偉、潘文健、李振誥，「濁水溪沖積扇與屏東平原地下水合適出水量之研究」，臺灣水利，第50卷，第3期，第70-82頁，2002。
35. 陳俊焜，「濁水溪沖積扇地下水資源調配與管理之研究」，國立成功大學資源工程所碩士論文，1998。
36. 陳美惠，「屏東隘寮圳灌區水資源調配運用及管理最佳化之研究」，國立成功大學資源工程所碩士論文，2005。
37. 陳尉平，「應用河川流量歷線推估台灣地下水補注量」，國立成功大學資源工程所博士論文，2006。
38. 陳愛光、李慈君、曹劍鋒，「地下水資源管理」，地質出版社，1991。
39. 曾鈞敏，「地下水模式應用於地下水資源管理可行性之探討」，水資源管理會刊，第7卷，第1期，第49-56頁，2005。
40. 黃清譽，「彰化地區地面水與地下水聯合供水管理調配之研究」，國立成功大學資源工程所碩士論文，2007。
41. 經濟部中央地質調查所，「台灣地區地下水觀測網第一期計畫八十一、八十二及八十三年度濁水溪沖積扇水文地質調查研究報告」，1995。
42. 經濟部中央地質調查所，「台灣地區地下水觀測網第一期計畫八十六年度濁水溪沖積扇水文地質鑽探（甲）研究報告」，1997。
43. 經濟部中央地質調查所，「濁水溪沖積扇水文地質調查研究總報告」，1999。
44. 經濟部水利署，「提升台灣水資源運用效率之研究(2/3)」，2005。
45. 經濟部水利署，「嚴重地層下陷地區劃設作業規範」，2005。
46. 經濟部水利署，「金門地下水資源調查分析(2/2)」，2007。
47. 經濟部水利署及工業技術研究院地層下陷資料庫，網址：http://www.subsidence.org.tw/index2.aspx
48. 經濟部水利署水文水資源管理供應系統網站。網址： http://www.wra.gov.tw/default.asp
49. 經濟部水資源局，「台灣地區地下水觀測網整體計畫81-84年度成果簡介」，經濟部水資源局，1997。
50. 經濟部水資源局，「台灣地區地層下陷監測調查分析」，2002。
51. 葉一隆，「斜率灰色模式與一維地下水流分析」，國立台灣大學農業工程學研究所博士論文， 2001。
52. 劉志文，「屏東沖積平原地下水資源調配與管理之研究」，國立成功大學地球科學研究所碩士論文，2000。
53. 劉振宇、林俊男、洪有仁、張誠信，「金門地區地面水與地下水聯合運用」，台灣水利，第55期，第2卷，第44-52頁，2007。
54. 劉聰桂，「由碳十四定年與氚示蹤探討濁水溪沖積扇地下水的流速與補注」，地下水資源與水質保護研討會論文專輯，1994。
55. 劉聰桂，「屏東平原地下水碳十四定年研究」，經濟部水資源統一規劃委員會委託調查報告，1992。
56. 劉聰桂、田巧玲、邱等輝、張炎銘，「濁水溪沖積扇之地下水資源—碳十四與氚定年/示蹤研究」。濁水溪沖積扇地下水及水文地質研討會論文集，第145-165頁，1996。
57. 劉聰桂、張國強、林文勝、陳冠宇、康素貞，「由碳十四與氚定年/示蹤探討嘉南平原水文地質架構與地下水補注」，第四屆地下水資源及水質保護研討會論文集，第255-258頁，2001。
58. 蔡長泰、高家俊，「地層下陷區土地利用對淹水潛勢影響之研究」，經濟部水利署委託專題研究計畫成果報告，2002。
59. 鄭文傑、倪至寬，「侷限含水層之暫態非全程貫入井抽降水試驗暨地盤水理特性研析」，臺北科技大學學報，第41-1期，2008。
60. 鄧聚龍，「灰色控制系統」，華中理工大學出版社，1997。
61. 譚義績、李峑毅、蘇明信，「台灣地區地下水管理決策系統之建議」，土木水利，第34卷，第4期，第101-111頁，2007。