林俊男臺灣大學：生物環境系統工程學研究所董仁Tung, JenJenTung2007-11-272018-06-292007-11-272018-06-292004http://ntur.lib.ntu.edu.tw//handle/246246/56061地下水流動是一個複雜的物理現象，由於數值方法及電子計算機的快速發展，多數學者以數學模型來模擬地下水流。以往在二維非拘限地下水流動，大都以水平（X-Y面）作為討論的對象，其控制方程式為非線性之偏微分方程式，而垂直斷面又涉及非飽和孔隙介質流之物理現象，使控制方程式更複雜，參數更多且不易確立。本文之目的，為減少非飽如孔隙介質流參數之推估，而不需引用非飽和層的理論。將單位時間內抽水量視為含水層排水量及孔隙壓縮效應所產生排水量二者之和，以有限元素法將非拘限含水層中之自由水面作為模擬區域的變動邊界，針對非拘限含水層之飽和部份，以有限元素法推估自由水面。 控制方程式應用有限元素法離散後，可獲得一線性聯立方程式。於穩態時以水頭值當作其自由水面高程進行疊代，直到二者差值絕對值小於允許誤差值E即為模擬自由水面高程。於抽水實驗時，本研究首先藉由調整孔隙壓縮產生之排水量，即為模式中之抽水量，使觀測井對應節點之洩降與實驗觀測井洩降之差值小於給定之允許誤差值EA。進而計算含水層中之有效孔隙率使孔隙壓縮產生之排水量與含水層洩降之排水量二者之總和等於實際抽水量。而本文模擬結果與利用雙對數標準曲線比對法及半對數分析法所求得之有效孔隙率，亦在同一個數階之內，顯示本模式之可靠性及可行性。Due to the rapid development of numerical analysis and computer, scholars employ mathematical model to simulate the groundwater flow, a complex physical phenomenon. It is used to discuss the flow in a two-dimensional unconfined aquifer in the horizontal (x-y) plane. The governing equation is a nonlinear partial differential one, and the vertical section involves the phenomenon of unsaturated porous media flow which has a more complex phenomenon and more uncertain parameters. Therefore, the purpose of the article is to reduce the estimate of parameters in an unsaturated porous media flow, and not to quote some theories of an unconfined aquifer. We calculate the total volume of water released from storage in the aquifer per unit time due to the elastic storage and the actual drainage of water from the aquifer, take the free surface as a moving boundary in a simulated region, and estimate the groundwater stage in the saturated region of unconfined aquifer by using the Finite Element Method. A linear system can be obtained by solving the governing equation using the Finite Element Method. We use iteration method and take head as groundwater stage until the absolute value of the difference is under the allowed value. This study makes sum of the drainage in this mode and drawdown cone equal to the total pumping rate by adjusting the effective porosity after making the absolute value of the difference between the stage of the observation well and the corresponding node under the allowed value EA by adjusting the drainage in this mathematical model. By comparing with the values of effective porosity obtained from normal logarithm curves and semilog curves analysis method, the result is in a same order, which confirms the reliability and feasibility of the method propose in this paper.目 錄 中文摘要……………………………………………I 英文摘要……………………………………………II 目錄…………………………………………………III 圖目錄………………………………………………IV 表目錄………………………………………………V 符號說明……………………………………………VI 第一章 緒論……………………………………1 第二章 文獻回顧………………………………4 第三章 模式理論………………………………7 3-1 地下水流理論……………………………7 3-1-1 控制方程式…………………………7 3-1-2 初始條件及邊界條件………………13 3-2 數值方法…………………………………14 3-3 系統解析…………………………………16 3-3-1 控制方程式之離散…………………16 3-3-2 空間近似……………………………17 3-3-3 時間系統解析………………………20 3-3-4 移動邊界建立………………………21 第四章 模式演算………………………………22 第五章 資料分析與結論建議…………………33 5-1 資料分析…………………………………33 5-2 結論………………………………………35 5-3 建議………………………………………36 參考文獻……………………………………………37 附錄A………………………………………………40 附錄B………………………………………………45 圖 目 錄 圖 1-1 研究步驟流程圖………………………3 圖 3-1 控制體積圖……………………………7 圖 4-1 Case 1之元素分割圖…………………23 圖 4-2 Case 1流程圖…………………………25 圖 4-3 Case 2之元素分割圖…………………28 圖 4-4 Case 2計算邊界流量流程圖…………30 圖 4-5 Case 2計算有效孔隙率流程圖………31 圖 5-1 洩降圓錐水量與時間關係圖…………33 表 目 錄 表 2-1 Lawrenceville市抽水試驗雙對數標準曲線 與半對數直接分析法之比較………………………………6 表 4-1 Case 1之Dirichlet邊界值…………………24 表 4-2 Case 1之Neumann邊界值……………………24 表 4-3 Case 1允許誤差值0.3之節點水頭值………26 表 4-4 Case 1之節點水頭值………………………26 表 4-5 Case 2 之相關含水層參數…………………27 表 4-6 Case 2 之Dirichlet邊界值…………………29 表 4-7 Case 2 之Neumann邊界值……………………29 表 4-8 Case 2 之洩降水位…………………………29 表 4-9 Case 2 之累積抽水量………………………30 表 5-1 有效孔隙率比較……………………………34458224 bytesapplication/pdfen-US有限元素法移動邊界孔隙壓縮效應Unconfined aquifersFinite element MethodMoving boundarythesishttp://ntur.lib.ntu.edu.tw/bitstream/246246/56061/1/ntu-93-R87622013-1.pdf