序率譜方法制定異質性含水層水井捕集區

以作者查詢圖書館館藏

、以作者查詢臺灣博碩士

、以作者查詢全國書目

、勘誤回報

、線上人數：7

、訪客IP：3.140.186.241

姓名

林淇平(Che-Ping Lin) 查詢紙本館藏

畢業系所

應用地質研究所

論文名稱

序率譜方法制定異質性含水層水井捕集區
(Stochastic spectral method to delineate well capture zones in heterogeneous aquifers)

相關論文

★ 延散效應對水岩交互作用反應波前的影響	★ 跨孔式注氣試驗方法推估異質性非飽和層土壤氣體流動參數
★ 現地跨孔式抽水試驗推估異質性含水層水文地質特性	★ iTOUGH2應用於實驗室尺度非飽和土壤參數之推估
★ HYDRUS-1D模式應用於入滲試驗推估非飽和土壤特性參數	★ 沿海含水層異質性對海淡水交界面影響之不確定性分析
★ 非拘限砂質海岸含水層中潮汐和沙灘坡度水文動力條件影響苯傳輸	★ 利用MODFLOW配合SUB套件推估雲林地區垂向平均長期地層下陷趨勢
★ 高雄平原地區抽水引致汙染潛勢評估	★ 利用自然電位法監測淺層土壤入滲歷程
★ 利用LiDAR點雲及影像資料決定露頭節理結合面之研究	★ 臺灣西部沿海海水入侵與地下水排出模擬分析
★ 三氯乙烯地下水污染場址整治後期傳輸行為分析¬-應用開源FreeFEM++有限元素模式架構	★ 都會地區滯洪池增設礫石樁之入滲效益模擬與分析
★ 利用數值模擬探討二氧化碳於異向性及異質性鹽水層之遷移行為	★ 三維離散裂隙網路水流與溶質傳輸模式發展

檔案

[Endnote RIS 格式]

[Bibtex 格式]

[相關文章]

[文章引用]

[完整記錄]

[館藏目錄]

[檢視]

[下載]

本電子論文使用權限為同意立即開放。
已達開放權限電子全文僅授權使用者為學術研究之目的，進行個人非營利性質之檢索、閱讀、列印。
請遵守中華民國著作權法之相關規定，切勿任意重製、散佈、改作、轉貼、播送，以免觸法。

摘要(中)

水井捕集區（well capture zone）的制定，對地下水的水源劃定精確的水井保護區(well head protection area, WHPA)非常重要。由於自然的含水層系統的水力參數包含著不同程度的異質性，這些參數的變異性將直接地影響流場的測量與水井捕集區的制定，序率方法可量化由於自然含水層複雜的異質性以及含水層參數獲得的限制所造成的不確定性。本研究中，利用一階近似法（first-order method）之非定常性譜方法（Nonstationary spectral method, NSM）與譜方法近似解（approximate spectral method, ASM）配合反向質點追蹤法（backward particle tracking），發展出利用序率模式制定異質性受壓含水層水井捕集區。在求解水井捕集區的過程中，先利用一階近似法求解地下水流速的不確定性，並利用質點追蹤的直接傳遞概念定義在特定時間內水井捕集區所形成的不確定性區間。本研究利用數個測試例量化發展的一階近似法之精確度，其中包含有限邊界地下水流場、多井系統流場、多尺度水力傳導係數以及由於源匯項所造成的非定常性流場，並考慮不同程度的異質性(ln K 小尺度擾動量之變異數 =0.1、0.5與 1.0)，推估由於含水層異質性對制定水井捕集區所造成的影響。由本研究測試例結果顯示，在相同的水文條件下，ASM與NSM所推估出的結果與蒙地卡羅法（Monte Carlo method, MCS）近似，譜方法近似解計算效率快速，並且所推估出的結果與NSM符合。由模擬的結果得知，ASM所需要的計算機計算時間僅數秒至數分鐘，而NSM需要數分鐘至數小時，MCS則需要長達數小時甚至數週的數值計算時間。ASM法的提出，包含了解析與數值兩種序率理論的優點，提供了現實中解決序率地下水模擬問題的方法。

摘要(英)

The delineation of well capture zones is of great importance to accurately define well head protection area (WHPA) for potential groundwater resources. Because natural aquifer systems typically involve different extent of heterogeneity in aquifer parameters, such parameter variations will then directly influence the estimations of flow fields and the delineations of well capture zones. Stochastic methods can be used to quantify uncertainties that are caused by natural complexity of aquifer heterogeneity and limited capability to obtain sufficient aquifer parameters. The ultimate results of stochastic delineation of capture zones are to contribute the remediation designs and risk analysis for contaminated aquifers. In this study, two first-order methods, including nonstationary spectral method (NSM) and approximate spectral method (ASM), associated with particle backward tracking algorithm are developed to delineate stochastic capture zones in heterogeneous confined aquifer systems. The Monte Carlo simulation (MCS) is employed for validating the developed first-order methods (i.e., ASM and NSM). To obtain a stochastic capture zone, the developed first-order methods are first employed to solve for velocity uncertainties. The concept of direct propagation of uncertainties of particle tracks is then used to define the uncertainty bandwidth of the stochastic capture zone for a specified time. In this study the developed first-order methods are first assessed to quantify the accuracy of delineated well capture zones under a variety of conditions, including bounded flow domains, multiple wells flow systems, multiple hydraulic conductivity scales, and nonstationary flows caused by complex sources and sinks in the modeling areas. Additionally, different degrees of heterogeneity (i.e., the small-scale ln K variance =0.1, 0.5, and 1.0) are considered in all modeling cases to get general insight into the effect of aquifer heterogeneity on the delineations of well capture zones. The simulation results reveal that the solutions of the NSM and ASM agree reasonably with those of MCS under same hydrogeological conditions. The ASM is efficient and the ASM solutions agree well with the full version of numerical spectral method NSM. In general the computation times for the ASM solutions are in the orders of seconds to minutes based on the computational system with i7 CPU system, while the computation times for the NSM solutions are in the orders of minutes to hours. However, the computation times for the costly MCS are in the orders of hours to possibly several weeks. The proposed ASM method has conceptually taken the advantages of analytical and numerical stochastic theories and provided an opportunity to include stochastic theories in practical groundwater modeling problems.

關鍵字(中)

★ 蒙地卡羅模擬法
★ 不確定性
★ 反向質點追蹤
★ 譜方法
★ 水井捕集區
★ 序率
★ 含水層異質性
★ 一接近似法

關鍵字(英)

★ well capture zone
★ stochastic method
★ aquifer heterogeneity
★ first-order method
★ particle backward tracking
★ velocity uncertainty
★ Monte Carlo simulation

論文目次

摘要 i
ABSTRACT ii
誌謝 iv
目錄 v
圖目錄 vii
表目錄 xiv
符號說明 xv
第一章緒論 1
1.1 研究背景 1
1.2 文獻回顧 2
1.2.1 定率模式（Deterministic model） 2
1.2.2 序率模式（Stochastic model） 4
1.3 研究動機與目的 5
1.4 論文架構 7
第二章研究理論 8
2.1 研究理論概述 8
2.2 地下水水流控制方程式 8
2.3 質點追蹤法（Particle tracking method） 9
2.4 點在多邊形內演算法（Point in a polygon algorithm）求解質點均值與位移變異區間 14
2.5 蒙地卡羅模擬法（Monte Carlo simulation, MCS） 16
2.6 一階近似法（First order method） 21
2.6.1 一階近似法地下水水流與達西定律之均值與擾動方程式 21
2.6.2 非定常性譜方法（Nonstationary spectral method, NSM） 22
2.6.3 譜方法近似解（Approximate spectral method, ASM） 24
2.6.4 一階近似法求解質點均值軌跡（mean trajectory）與位移變異數（displacement variance） 25
第三章模式例建立 28
3.1 模式流場建立 28
3.2 定值均值水力傳導係數流場 29
3.3 非定常性均值水力傳導係數流場 30
3.4 多井系統非定常性均值水力傳導係數流場 31
3.5 複合性地下水流系統流場 33
第四章結果與討論 36
4.1 定常性水力傳導係數流場模擬結果 36
4.2 非定常性水力傳導係數流場模擬結果 50
4.3 多井系統非定常性水力傳導係數流場模擬結果 64
4.4 複合系統非定常性地下水流流場 78
4.5 討論 91
第五章結論與建議 93
5.1 結論 93
5.2 建議 94
參考文獻 95

參考文獻

[1] U.S. EPA. (1994) Ground Water and Wellhead Protection Handbook. EPA-625/R-94-001. Washington, DC: U.S. Environmental Protection Agency, Office of water.
[2] Paradis, D., Martel, R., Karanta, G., Lefebvre, R., Michaud, Y., Therrien, R. and Nastev, M. (2007) “Comparative study of methods for WHPA delineation”. Ground Water, Vol 45, No.2, 158-167.
[3] Bear, J. and Jacobs, M. (1965) “On the movement of water bodies injected into aquifer”, Journal of Hydrology 3 NO.1: 37-57.
[4] Todd, D.K. (1980) Groundwater Hydrology, 2nd ed. New York: John Wiley.
[5] Javandel, I. and Tsang, C.F. (1987) “Capture-zone type curves: a tool for aquifer cleanup” Ground Water, 24(5): 616-625.
[6] Zhan, H. (1998) “Analytical and numerical modeling of a double well capture zone” Journal of Mathematical Geology, 31(2): 175-193.
[7] Zhan, H. (1999) “Analytical study of capture zone to horizontal well” Journal of Hydrology 217: 46-54.
[8] Shan, C. (1999). “An analytical solution for the capture zone of two arbitrarily located wells”. Journal of Hydrology 222: 123-128.
[9] Rock, G. and Kupfersberger, H. (2002). “Numerical delineation of transient capture zones”. Journal of Hydrology 269: 134-149.
[10] Pollock, D.W. (1998). Documentation of computer programs to compute and display pathlines using results from the U.S. Geological Survey modular three dimensional finite-difference ground-water flow model. USGS Open-File Report 89-381. Reston, Virginia: Office of Ground Water.
[11] Franz, T. and Guiguer, N. (1990), “FLOWPATH, two-dimensional horizontal aquifer simulation model”. Waterloo Hydrogeologic Software, Waterloo.
[12] U.S. EPA. (1991). WHPA A modular semi-analytical model for the delineation of wellhead protection zones. Office of Groundwater Protection, Washington, DC.
[13] Kinzelbach, W., Marburger, M. and Chiang, W.H. (1992). “Determination of groundwater catchment areas in two and three spatial dimensions”, Journal of Hydrology. 134, 221-246.
[14] Schafer, D.C. (1996). “Determining 3D capture zone in homogeneous, an isotropic aquifer” Ground Water 34(4):628-369.
[15] McDonald, M.G. and Harbaugh, A.W. (1998). A modular three dimensional finite-difference ground-water flow model. USGS Techniques of Water Resources Investigations Book 6, Ch. A1, 586p.
[16] Zheng, C. and Bennett, G.D. (2002). Applied Contaminant Transport Modeling. John Wiely and Sons, Inc., New York.
[17] Shafer, J. M. (1990). GWPATH. Version 4.0 Champaign, Illinois State Water Survet Division.
[18] Bair, E. S. (1991). “Delineation of travel time-related capture areas of wells using analytical flow models and particle-tracking analysis.” Ground Water 29, no. 3: 387-397
[19] Haitjema, H., Wittman, J., Kelson, V. and Bauch, N. (1994). “WHAEM: Program documentation for the wellhead analytic element model”, EPA/600/R-94/210. Washington, DC: U.S. EPA, Office of Research and Development.
[20] Franzetti, S. and Guadagnini, A. (1996). “Probabilistic estimation of well catchments on heterogeneous aquifers”, Journal of Hydrology 174: 149-171.
[21] Bakr, M.I. and Bulter, A.P. (2004). “Worth of head data in well capture zone design: deterministic and stochastic analysis”. Journal of Hydrology 290 202-216.
[22] Bakr, M.I. and Bulter, A.P. (2005). “Nonstationary stochastic analysis in well capture zone design using first-order Taylor’s series approximation”. Water Resource Research. Vol.41, W01004.
[23] Van Leeuwen, M., Chris, B.M. te Stroet., Butler, A.P. and Tompkins J.A. (1999). “Stochastic determination of the Wierden (Netherkands) capture zones”. Ground Water 37, 8-17.
[24] Guadagnini, A. and Franzetti, S. (1999). “Time-related capture zones contaminants in randomly heterogeneous formations”, Ground Water 97(2):253-260.
[25] Kunsrmann, H. and Kinelbach, M. (2000). “Computation of stochastic wellhead protection zones by combining the first-order second-moment and Kolmogorov backward equation analysis”, Journal of Hydrology. 297(3-4):127-146.
[26] Lu, Z. and Zhang, D. (2003). “On stochastic study of well capture zones in bounded, randomly heterogeneous media”, Water Resources Research, Vol. 39, NO.4, 1100.
[27] Zhang, D. and Lu, Z. (2004). “Stochastic delineation of well capture zones”. Stochastic Environmental Research and Risk Assessment, 18(1): 39-46.
[28] Stauffer, F., Guadagnini, A., Butler, A., Franssen, H-J. H., Van De Wiel, N., Bakr, M., Riva, M. and Guadagnini, L. (2005). “Delineation of source protection zones using statistical methods”. Water Resources Management 19, 163-185.
[29] Kunstmann, H., Kastens, M. (2006) “Determination of stochastic wellhead protection zones by direct propagation of uncertainties of particle tracks”. Journal of Hydrology, Vol. 323, 215-229
[30] Riva, M. Guadagnini, A. and De Simoni, M. (2006). “Assessment of uncertainty associated with the estimation of well catchments by moment equations”. Advances in Water Resources 29: 679-691.
[31] Guadagnini, A. Sanchez-vila, X., Riva, M. and De Simoni, M. (2003). “Mean travel time of conservative solutes in randomly heterogeneous unbounded domains under mean uniform flow”, Water Resource Research. 39, 1050.
[32] Hernadez, A. F., Neuman, S.P., Guadagnini, A. and Carrera, J. (2003). “Conditioning mean steady state flow on hydraulic head and conductivity through geostatistical inversion”, Stochastic Environmental Research and Risk Assessment, 17, 329-338.
[33] Li, S.G., Liao, H.S., Ni, C.F. (2004). “Stochastic modeling of complex nonstationary groundwater systems”. Advances in Water Resources 27:1087-1104
[34] Ni, C.F. and Li, S.G. (2005). “Simple closed form formulas for predicting groundwater flow model uncertainty in complex, heterogeneous trending media”. Water Resources Research. 41(11), W11503.
[35] Ni, C.F. and Li, S.G. (2006). “Modeling groundwater velocity uncertainty in nonstationary composite porous media”. Advances in Water Resources 29(12), 1866-1875.
[36] Ni, C.F., Li, S.G. (2009). Efficient conceptual framework to quantify flow uncertainty in large-scale, highly nonstationary groundwater systems. Journal of Hydrologic Engineering (in press).
[37] 陳詹閔，邱式鴻，2004，『整合建物平面資訊和數值覆面模型自動建構建物三維模型』，第二十三屆測量學術及應用研討會論文集，第673-680頁。
[38] Gelhar, L.W. (1986). “Stochastic analysis of phreatic aquifers”, Water Resource Research Vo1. 10., NO3, pp.539-545.
[39] Gelhar, L.W. (1993). Stochastic Subsurface Hydrology. Englewood Cliffs, NJ: Prentice-Hall.
[40] Li, S.G. and McLaughlin, D.B. (1991) “A non-stationary spectral method for solving stochastic groundwater problems: unconditional analysis.” Water Resources Research. 27(7): 1598-605.
[41] Li, S.G. and McLaughlin, D.B. (1995) “Using the nonstationary spectral method to analyze floe through heterogeneous trending media”, Water Resources Research. 31(3): 541-51.

指導教授

倪春發(Chuen-Fa Ni)

審核日期

2009-7-20

推文