臺灣博碩士論文加值系統

本研究目的為建立一區域性地下水流模式之率定方法。首先，先對地下水區域進行主成分分析並找出其時空間變數，並以其為基準增設新輔助井，再以經驗正交函數法配合地下水流模式之蓄水量變化歷線與模擬誤差歷線，快速且準確地掌握各補注量與水文地質參數之時間與空間分布資訊，並針對地下水系統之地面水補注量以及水文地質參數進行修正，最後將該方法應用於名竹盆地。
本研究欲建立一地下水流數值模式率定方法，首先，先對原觀測井進行主成分分析，得到該區域之時空間變數分布，接著找出特徵值為0之交線，並於交點上新增輔助井。之後，設定目標函數使模擬水位與觀測水位之均方根誤差(RMSE值)最小，而本研究中選定之決策變數為水平水力傳導係數、垂直水力傳導係數以及各補注量之時空分布；限制式主要分為三條：(1)各補注量於率定過程中須符合質量守恆；(2)地下水位之模擬結果須符合地下水流控制方程式；(3)水平水力傳導係數與垂直水力傳導係數有一合理範圍之限制。求解流程首先為設定決策變數初始值，輸入地下水流模式進行地下水位之模擬，並計算其目標函數判斷是否達成停止條件，若否則計算蓄水量誤差歷線，並利用經驗正交函數分析地下水流模式之蓄水量模擬誤差歷線，計算決策變數之修正量，以進行各補注量與水文地質參數之修正，完成一次迭代過程，經過數次迭代求解達到停止條件後便完成地下水流參數優選模式之率定，獲得各淨補注量與水文地質參數之最佳時空分布。
本研究將所建立之參數率定模式應用於名竹盆地地下水流數值模式，其地下水流模式之模擬年限為2012年1月至2012年12月，以日為時間單位進行模擬，求解過程中同時針對兩層含水層中之水平水力傳導係數分區、垂直水力傳導係數分區以及第一含水層中之降雨補注量、邊界補注量、河水補注量時空分布進行率定計算修正量。率定結果發現，於模式迭代初期，其地下水位模擬誤差RMSE值下降幅度最大，於幾次迭代之後其地下水位模擬誤差RMSE值幾乎呈現穩定狀態不再下降而達到停止條件，完成模式率定。率定完成之水平水力傳導係數與垂直水力傳導係數皆在合理範圍內，其地下水流數值模式能夠準確模擬含水層(一)之水位變化趨勢與豐枯水期地下水位值；含水層(二)之地下水位模擬結果亦能大致抓到水位變化之趨勢，且在新增輔助井後，與無輔助井之模擬結果比較後可發現，地下水位模擬情形改善許多。顯示本研究方法能夠快速且準確的掌握地下水補注量與水文地質參數之時間與空間分布資訊，並反饋回地下水流數值模式中，以獲得準確且良好之地下水流數值模式。

This study is aimed to develop a regional groundwater numerical model calibration method. First, use principal component analysis (PCA) on the groundwater section to find out its temporal-spatial variable,and use it as a reference to create new assistance well. Then, applies empirical orthogonal function (EOF) with the change hydrograph of groundwater storage and simulated error hydrograph of groundwater level to quickly and accurately catch and calibrate the temporal-spatial distribuation of water recharge and hydrogeological parameters. The established method was applied on the groundwater system of Ming Chu Basin.
This study is aimed to develop a groundwater numerical model calibration method. First, use principal component analysis (PCA) on the groundwater section to get its temporal-spatial variable distribution, and finding the line of eigenvalue=0, and create the new assistance well on it. After, setting the objective function is minimizing the the root mean square error (RMSE) of the simulated and observed groundwater level. The decision variables are horizontal hydraulic conductivity, vertical hydraulic conductivity and water recharge. There are three constraints of the optimization model: (1) the water recharge of groundwater system in every iteration of calibrating process must obey the mass balance; (2) the simulated groundwater level must follow the governing equation of groundwater flow; (3) the value of horizontal hydraulic conductivity and vertical hydraulic conductivity are restricted to a reasonable limits. The process of the optimization model sets the initial value of decision variables first, and inputs the variables to groundwater model. Thus, the groundwater level can be simulated and the objective function will be estimated. If the objective function doesn’t satisfy the stop condition, the simulated error hydrograph of groundwater level will be calculated and analyzed with EOF. Then, the modified decision variables is calculater according to the simulated error hydrograph of groundwater level end the result of EOF analysis. From iterations, the optimal temporal-spatial distribuation of surface water recharge and hydrogeological parameters can be obtain.
This study applied the model on the calibration of the groundwater system in Ming Chu Basin. The simulated period is from January 2012 to December 2012 daily. The decision variables were selected in this study are horizontal hydraulic conductivity, vertical hydraulic conductivity of two acqufiers and rain water recharge, river water recharge and boundary water recharge of hydraulic conductivity in first acquifer. The result show that the RMSE is decreased dramatically in early iteration of the calibration and become smoothly after several iterations. The calibrated hydraulic conductivity and vertical leakence are in reasonable limits. The simulated groundwater level can reflect the approximately trendance in all acquifer and can capture the peak of the observed value in first acquifer.
Hence, the established method of this study can effectively and accurately calibrate temporal-spatial distribution of surface water recharge and hydrogeological parameters.

論文口試委員審定書 #
誌謝 I
中文摘要 II
ABSTRACT IV
目錄 VI
圖目錄 VIII
表目錄 X
第一章 緒論 1
1.1 研究緣起 1
1.2 研究目的 2
1.3 研究流程 2
第二章 文獻回顧 4
2.1 地下水系統流出流入量辨識 4
2.2 地下水流數值模式 7
2.3 主成分分析 9
2.4 經驗正交函數 10
第三章 方法建立 12
3.1 地下水流數值模式率定優選模式 12
3.1.1 優選模式之建立 12
3.1.2 水平衡分析與地下水流數值模式 16
3.1.3 優選模式之求解 26
3.2 以經驗正交函數改進優選模式之求解 28
3.2.1 以EOF分析蓄水變化量歷線與蓄水量誤差歷線 29
3.2.2 合理假設與計算決策變數修正量公式 30
3.2.3 地下水補注量之初始設定與誤差修正限制式 36
3.3 利用主成分分析新設輔助井於增進參數率定精度 39
3.3.1 新設輔助井網之設計流程 39
3.3.2 利用主成分分析決定輔助井位置 40
第四章 方法應用 43
4.1 研究區域概述 43
4.2 含水層概念分層與邊界辨識 50
4.3 地下水系統流出流入量辨識 53
4.4 地下水流數值模式建立 60
4.4.1 模擬範圍界定與邊界條件設定 60
4.4.2 地層分層與數值網格設定 61
4.4.3 水文地質參數分區與設定 62
4.4.4 河流演算套件設定 66
4.4.5 抽水量設定 69
4.4.6 地下水系統補注量分區與初始設定 69
4.5 地下水流數值模式率定結果 71
第五章 結論與建議 89
5.1 結論 89
5.2 建議 90
參考文獻 91

1.經濟部中央地質調查所，「濁水溪中游地區地下水補注調查與評估」，1996。
2.經濟部中央地質調查所，「台灣地區地下水觀測網第一期計畫濁水溪沖積扇水文地質調查研究總報告」，經濟部水資源局委託計畫報告，1999。
3.江崇榮、陳瑞娥、賴慈華、黃智昭，「濁水溪沖積扇地下水區之補注水源評估」，第二屆資源工程研討會，第86-91頁，2005。
4.江崇榮、黃智昭、陳瑞娥，「以地下水歷線分析法評估濁水溪沖積扇之地下水收支」，經濟部中央地質調查所彙刊，第十九號，第61-89頁，2006。
5.中興工程顧問股份有限公司，「濁水溪沖積扇地面地下水聯合運用管理模式建立與機制評估」，經濟部水利署中區水資源局委託計畫報告，2007。
6.謝壎煌、陳忠偉、葉信富、李振誥，「應用河道水位變化評估新虎尾溪地下水補注量之研究」，農業工程學報，第53卷第2期，第50-60頁，2007。
7.經濟部水利署，「雲林內陸地層下陷地區地下水抽水機制之探討與評估期末報告」，2011。
8.陳聖傑，「地下水水位之主成分分析-以濁水溪沖積扇為例」，國立成功大學資源工程學系碩士論文，2012。
9.王韋勳，「名竹盆地地下水流數值模式之建立與應用」，國立台灣大學土木工程學系碩士論文，2012。
10.林聖婷，「濁水溪沖積扇補注量與抽水量空間分佈模式建立」，國立台灣大學土木工程學系碩士論文，2012。
11.經濟部水利署，「經濟部水利署水文水資源資料管理供應系統」網站， 2012，
http：//gweb.wra.gov.tw/wrweb/。
12.經濟部中央地質調查所，「地下水補注地質敏感區畫定計畫書」，2014。
13.徐年盛、黃建霖、劉宏仁、李冠緯，「地下水抽補量與參數時空分布率定之方法建立與數值實驗驗證」，農業工程學報，第61卷第3期，第1-28頁，2015。
14.行政院農業委員會，「提升灌溉用水效率策略規劃與優化灌溉管理決策輔助機制」，2015。
15.許富建，「區域性地下水系統水流模式率定方法建立與應用-以濁水溪沖積扇為例」，國立台灣大學土木工程學系碩士論文，2016。
16.Arnold, J. G., Muttiah, R. S., Srinivasan, R., and Allen, P. M. (2000). Regional estimation of base flow and groundwater recharge in the Upper Mississippi river basin. Journal of Hydrology, 227(1-4), 21-40.
17.Becker, L., and Yeh, W. W. G. (1972), Identification of parameters in unsteady open channel flows, Water Resources Research, 8(4), 956-965.
18.Dong, D. L., Wu, Q., Zhang, R., Song, Y. X., Chen, S. K., Li, P., Liu, S. Q., Bi, C. C., Lu, Z. Q., and Huang, S. L. (2007). Environmental Characteristics of Groundwater: an Application of PCA to Water Chemistry Analysis in Yulin. Journal of China University of Mining and Technology, 17(1), 73-77.
19.Gerhart, J. M. (1986). Groundwater recharge and its effects on nitrate concentration beneath a manured field site in Pennsylvania. Groundwater,24(4),483-489.
20.Hannachi, A., Jolliffe, I. T., and Stephenson, D. B. (2007). Empirical orthogonal functions and related techniques in atmospheric science: A review. International Journal of Climatology, 27(9), 1119-1152.
21.Healy, R. W., and Cook, P. G. (2002). Using groundwater levels to estimate recharge. Hydrogeology Journal, 10(1), 91-109.
22.Hsu, N. S., Chiang, C. J., Wang, C. H., Liu, C. W., Huang, C. L., and Liu, H. J. (2013). Estimation of pumpage and recharge in alluvial fan topography under multiple irrigation practices. Journal of Hydrology, 479, 35-50.
23.Igboekwe, M. U., Gurunadha Rao, V. V. S., and Okwueze, E. E. (2008). Groundwater flow modelling of Kwa Ibo River watershed, southeastern Nigeria. Hydrological Processes, 22(10), 1523-1531.
24.Lee, C. H., Yeh, H. F., and Chen, J. F. (2008). Estimation of groundwater recharge using the soil moisture budget method and the base-flow model. Environmental Geology, 54(8), 1787-1797.
25.Martinez-Santos, P., and Martinez-Alfaro, P. E. (2010). Estimating groundwater withdrawals in areas of intensive agricultural pumping in central Spain. Agricultural Water Management, 98(1), 172-181.
26.Meinzer, O. E. (1933). Geophysical interpretation of ground-water levels. Transactions-American Geophysical Union,14,36-37.
27.Moon, S. K., Woo, N. C., and Lee, K. S. (2004). Statistical analysis of hydrographs and water-table fluctuation to estimate groundwater recharge. Journal of Hydrology, 292(1-4), 198-209.
28.Munoz, B., Lesser, V. M., and Ramsey, F. L. (2008). Design-based empirical orthogonal function model for environmental monitoring data analysis. Environmetrics, 19(8), 805-817.
29.Shorieh, A., Porel, G., Razack, M. (2015). Assessment of Groundwater Quality in the Dogger Aquifer of Poitiers, Poitou-Charentes Region, France. Journal of Water Resource and Protection, 7(3), 171-182.
30.Takounjou, A. F., Gurunadha Rao, V. V. S., Ngoupayou, J. N., Nkamdjou, L. S., and Ekodeck, G. E. (2009). Groundwater flow modelling in the upper Anga’a river watershed, Yaounde, Cameroon. African Journal of Environmental Science and Technology, 3(10), 341-352.
31.Yu, H. L., and Lin, Y. C. (2015). Analysis of space-time non-stationary patterns of rainfall-groundwater interactions by integrating empirical orthogonal function and cross wavelet transform methods. Journal of Hydrology, 525, 585-597.
32.Yu, H. L., and Chu, H. J. (2010). Understanding space-time patterns of groundwater system by empirical orthogonal functions: A case study in the Choshui River alluvial fan, Taiwan. Journal of Hydrology, 381(3-4), 239-247.