台灣西部沿海超抽地下水嚴重，使得長年地下水供需失衡，造成
該地嚴重地層下陷，其中雲林地區為全國下陷面積最大，下陷速率全
台第二之地區，也因此雲林地區之地層下陷研究成為重要課題。
本研究首先利用地下水流數值模擬軟體MODFLOW 建構濁水溪
沖積扇地下水流數值模型，接著藉由參數率定求得符合現地之水文地
值概念模式，結合層間貯水套件（Interbed Storage Package），以建立
雲林地區地層下陷模型，最後評估不同的水位資料型式，對地層下陷
量之計算差異。
本文將地下水流數值模型之觀測水井資料分成月帄均水位(模式
一)、當月最高水位(模式二)與當月最低水位(模式三)，分別進行地層
下陷模擬，並將其結果與雲林地區10 口地層下陷監測井比較，發現
下陷趨勢大致相符，且主要下陷區域為雲林縣之崙背鄉、二崙鄉、褒
忠鄉、元長鄉、口湖鄉、水林鄉。模式一以元長國小、虎尾國小與客
厝國小等靠內陸地區之監測站較為準確；模式二以建陽國小與金湖國
小等沿海偏南之監測站較佳；模式三以崙豐國小、新興國小、龍岩國
小與峰安國小等沿海偏北之監測站較符合。

There have serious problems of overpumping of groundwater in the
coast area of the western Taiwan, especially in Yunlin area which has the
most subsidence area and the serious subsidence rate in Taiwan. If
subsidence keeps growing, it will make the supply and demand of
groundwater unbalanced, cause heavy subsidence and affect the
geo-engineering constructions of the country.
The purpose of this thesis is to study the amount of land subsidence
in Yunlin area. The numerical simulation groundwater flow software
MODFLOW is applied to construct numerical models of groundwater
flow at Chou-Shui alluvial fan. By setting hydraulic parameters, the
combined hydrological model with the code (Interbed Storage Package)
can be built the model of land subsidence in Yunlin. Then we evaluate the
difference in amount of land subsidence under three types of water level
observation data, average water level of the month (model 1), the highest
water level of the month (model 2), and the lowest water level of the
month (model 3). Finally, the difference between collected data and the
simulation results of the land subsidence is compared.
As compared between field measurement data and simulation results
with three models, it indicated that the trend of simulation results of the
land subsidence corresponds to the filed measurement of land subsidence.
In subsidence stations in inland area of Yunlin, such as Yuan-Chang
elementary school, Hu-Wei elementary school and Ke-Tsuo elementary
school, simulation results of subsidence in model 1 had more accurate
results of subsidence compared with field measurement. In subsidence
stations such as Jian-Yang elementary school and Jin-Hu elementary
school of southern coast, simulation results of subsidence in model 2 had
better correction results of subsidence. In subsidence stations of northern
III
coast area such as Lun-Feng elementary school, Shin-Shing elementary
school, Lung-Yan elementary school and Feng-An elementary school,
simulation results of subsidence in model 3 had better correction of
subsidence.

1. Chiang, W. H., and Kinzelbch, W. ―Processing Modflow for windows
(version 5.06)‖, Hamburg, Germany, 1998.
2. Don N. C., Hang N. T. M., Araki H, Yamanishi H, and Koga K.
―Groundwater resources and management for paddy field irrigation
and associated environmental problems in an alluvial coastal lowland
plain‖, Agricultural Water Management Vol. 84, Issue 3, 295-304,
2006.
3. Hanson, R. T., ‖Aquifer-System Compaction, Tucson Basin and Avra
Valley, Arizona‖, U.S. Geological Survey Water-Resources
Investigations Report 90-4178, 41P, 1989.
4. Helm, D. C., ―One-Dimensional Simulation of Aquifer System
Compaction Near Pixley, California 1. Constant Parameters‖, Water
Resources Research, vol. 11, 465-478, 1975.
5. Helm, D. C., ―One-Dimensional Simulation of Aquifer System
Compaction Near Pixley, California 2. Stress-Dependent Parameters‖,
Water Resources Research, vol. 12, 375-391, 1976.
6. Holzer, T. L., ―Preconsolidation Stress of Aquifer Systems in Areas
of Induced Land Subsidence‖, Water Resources Research, Vol. 17,
No. 3, 693-704, 1981.
7. Ireland, R. L., Poland, J. F., and Riley, F. S., ―Land subsidence in the
san Joaquin Valley, California, as of 1980‖, U.S. Geological Survey
Professional Paper 437I, 93P, 1984.
8. Jorgensen, D. G., ―Relationships between basic soils-engineering
equations and basic groundwater flow equations‖, U.S. Geological
Survey Water-Supply Paper, 2064, 40P, 1980.
9. Larson, K. J. ,Basagaolu, H., and Marino M. A. ―Prediction of
optimal safe ground water yield and land subsidence in the Los
Banos-Kettleman City area, California, using a calibrated numerical
54
simulation model‖, Journal of Hydrology, 242, 79-102, 2001.
10. Leake, S. A., ―Interbed storage changes and compaction in models of
regional groundwater flow‖, Water Resources Research, Vol.26, No.3,
13-25, 1990.
11. Leake, S. A. and Prudic, D. E. ―Documentation of a computer
program to simulate stream-aquifer relations using a modular,
finite-difference, ground-water flow model‖, U.S. Techniques of
Water-Resources Investigations of the United States Geological
Survey, book 6 Chapter A2,1991.
12. McDonald, M. G., and A. W. Harbaugh, ―A modular
three-dimensional finite-difference ground-water glow model, U.S.
Geological Survey‖, Virginia, 1988.
13. Meyer, W. R., and Carr, J. E., ―A digital model for simulation of
ground-water hydrology in the Houston area, Texas. Texas Dep.‖,
Water Resour. Rep. LP-103, 131PP.,Texas Dep. Of Water Resour.,
Austin.1979.
14. Neuman, S. P. ,Preller, C. and Narasimhan, T. N., ―Adaptive
Explicit-Implicit Quasi Three-Dimensional Finite Element Model of
Flow and Subsidence in Multiaquifer Systems‖, Water Resour. Res.
Vol. 18, 1151–1561, 1982.
15. Riley, F. S. ―Analysis of borehole extensometer data from central
California. In: Land Subsidence‖, Publ. 89, Vol. 2, edited by L. J.
Tison, 423-431, International Association of Hydrological Sciences,
Wallingford, England, 1969.
16. Wang, C. H., Kuo C. H.,. Peng T.R, Chen W.F., Liu T.K., Chiang C.J.,
Liu W. C., and Hung J. J., ―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,
55
1970.
18. Williamson, A. K., Purdic, D. E. and Swain, L. A., ―Ground-water
flow in the Central Valley, California‖, U.S. Geol. Surv. Prof. Pap.
1401-D, 127P, 1989.
19. 中央氣象局全球資訊網。網址：
http://www.cwb.gov.tw/V6/index.htm
20. 王士榮，徐國錦，「應用序率孔彈性理論探討地層下陷之行為」，
Journal of the Chinese Institute of Civil and Hydraulic Engineering,
Vol. 20, No. 3, pp. 357–371, 2008。
21. 台糖公司鑿井工程隊，「台灣地區地下水觀測網第一期計畫—地
下水觀測井開鑿及相關試驗(濁水溪沖積扇(二))」，1994。
22. 台糖公司鑿井工程隊，「台灣地區地下水觀測網第一期計畫—地
下水觀測井開鑿及相關試驗(濁水溪沖積扇(三))」，1995。
23. 行政院環境保護署，「環境白皮書」，行政院環境保護署，1998。
24. 何麗如、沌中和、陳鎮東，「北港地區地下水之地球化學初探」，
中華民國環境保護學會會誌，第13 卷，第2 期，第38-59 頁，1990。
25. 李友帄、張國強、劉萬里，「台灣地區之地下水資源現況」，第四
屆地下水資源及水質保護研討會論文集，第201-210 頁，2001。
26. 李振誥，陳尉帄，李如晃，「應用基流資料估計法推估台灣地下
水補注量」，台灣水利，第50 卷，第1 期，第69-80 頁，2002。
27. 李振誥、徐國錦，「地下水補注與管理措施」，工程，第80 卷，
第6 期，第80-99 頁，2007。
28. 李振誥、許清荃，「地下水資源調配與管理-以濁水溪沖積扇為案
例」，水資源管理，第8 卷，第24-31 頁，2000。
29. 李振誥、許清荃、林俶寬，「濁水溪沖積扇多層地下水調配與管
理之研究」，臺灣水利，第48 卷，第4 期，第41-52 頁，2000。
56
30. 沈向白，「濁水溪沖積扇地區地下水資源調查(含彰化縣、雲林縣)」，
台灣省水利局，1991。
31. 林文勝，「膨潤石脫水隊地層下陷之影響」國立台灣大學農業工
程所博士論文，2001。
32. 林文勝，「膨潤石脫水對地層下陷之影響」，國立台灣大學農業工
程所博士論文，2001。
33. 邱筱筠，「雲林虎尾灌區水資源調配與管理之研究」國立成功大
學資源工程所碩士論文，2008。
34. 柳志錫、楊秀隆、洪偉嘉、劉智超，「濁水溪沖積扇地層壓縮行
為探討」，第七屆地下水資源及水質保護研討會，第E33-E42 頁，
2009。
35. 徐享崑、劉豐壽、鄭昌奇，「臺灣地區地層下陷之現況，成因與
對策」，臺灣水利，第43 卷，第3 期，第19-29 頁，1995。
36. 財團法人工業技術研究院，「96 年度地層下陷彰雲地區地層下陷
監測及分析計畫」，2008。
37. 康家桂，「以時間序列預測麥寮地區地層下陷之研究」國立交通
大學土木工程所碩士論文,2005。
38. 張光仁，「台西麥寮地區多井抽水引致區域性地層下陷之分層分
析」國立台灣大學土木工程所碩士論文，2000。
39. 張良正、龔誠山，「區域性地下水觀測站網檢討八十四年度報告，
區域性地下水觀測站網檢討(I)－濁水溪沖積扇觀測站井佈置檢
討」，經濟部水利司，1995。
40. 許良賜，「使用調和分析推估阻水層水文特性研究」，國立成功大
學資源工程所碩士論文，2005。
41. 陳忠偉、潘文健、李振誥，「濁水溪沖積扇與屏東帄原地下水合
57
適出水量之研究」，臺灣水利，第50 卷，第3 期，第70-82 頁，
2002。
42. 陳俊焜，「濁水溪沖積扇地下水資源調配與管理之研究」，國立成
功大學資源工程所碩士論文，1998。
43. 陳建銘，「地層下陷模擬程序之建立與應用—以大城鄉西港地區
為例」國立成功大學土木工程所碩士論文，2005。
44. 陳尉帄，「應用河川流量歷線推估台灣地下水補注量」，國立成功
大學資源工程所博士論文，2006。
45. 曾鈞敏，「地下水模式應用於地下水資源管理可行性之探討」，水
資源管理會刊，第7 卷，第1 期，第49-56 頁，2005。
46. 黃清譽，「彰化地區地面水與地下水聯合供水管理調配之研究」，
國立成功大學資源工程所碩士論文，2007。
47. 經濟部中央地質調查所，「台灣地區地下水觀測網第一期計畫八
十一、八十二及八十三年度濁水溪沖積扇水文地質調查研究報
告」，1995。
48. 經濟部中央地質調查所，「台灣地區地下水觀測網第一期計畫八
十六年度濁水溪沖積扇水文地質鑽探（甲）研究報告」，1997。
49. 經濟部中央地質調查所，「濁水溪沖積扇水文地質調查研究總報
告」，1999。
50. 經濟部水利署，「提升台灣水資源運用效率之研究(2/3)」，2005。
51. 經濟部水利署，「嚴重地層下陷地區劃設作業規範」，2005。
52. 經濟部水利署及工業技術研究院地層下陷資料庫，網址：
http://www.subsidence.org.tw/index2.aspx
53. 經濟部水利署水文水資源管理供應系統網站。網址：
http://www.wra.gov.tw/default.asp
54. 經濟部水資源局，「台灣地區地下水觀測網整體計畫81-84 年度成
58
果簡介」，經濟部水資源局，1997。
55. 經濟部水資源局，「台灣地區地層下陷監測調查分析」，2002。
56. 劉振宇、林俊男、洪有仁、張誠信，「金門地區地面水與地下水
聯合運用」，台灣水利，第55 期，第2 卷，第44-52 頁，2007。
57. 劉聰桂，「由碳十四定年與氚示蹤探討濁水溪沖積扇地下水的流
速與補注」，地下水資源與水質保護研討會論文專輯，1994。
58. 劉聰桂，「屏東帄原地下水碳十四定年研究」，經濟部水資源統一
規劃委員會委託調查報告，1992。
59. 劉聰桂、田巧玲、邱等輝、張炎銘，「濁水溪沖積扇之地下水資
源—碳十四與氚定年/示蹤研究」。濁水溪沖積扇地下水及水文地
質研討會論文集，第145-165 頁，1996。
60. 劉聰桂、張國強、林文勝、陳冠宇、康素貞，「由碳十四與氚定
年/示蹤探討嘉南帄原水文地質架構與地下水補注」，第四屆地下
水資源及水質保護研討會論文集，第255-258 頁，2001。
61. 蔡長泰、高家俊，「地層下陷區土地利用對淹水潛勢影響之研究」，
經濟部水利署委託專題研究計畫成果報告，2002。
62. 賴典章、費立沊、江崇榮，「台灣地區地下水分區特性」，中央地
調所文地質調查與應用研討會，2003。
63. 蘇中鈺，「雲林地層下陷易淹區水資源調配運用與管理之研究」，
國立成功大學資源工程所碩士論文，2009。