透過您的圖書館登入
IP:44.197.214.36
  • 學位論文

變飽和地下水流數值模式之研究

A study of variably saturated numerical model for groundwater flow

指導教授 : 李天浩

摘要


本研究以混合型理查氏方程式變飽和度地下水流方程式,建構能同時模擬飽和與未飽和地下水流,確保模式能達成局部與全域質量守恆的數值模式。首先,探討一維模式在模擬時遭逢的不收斂或不穩定問題,研發解決的方式,以及利用多層次網格細化方法,提升模式計算效率;其次,建立垂直二維並考慮與坡面正交水平方向坡寬的擬三維變飽和度地下水流模式,測試與驗證滲流面邊界條件的處理方法,改進模式解矩陣的效率,探討坡地形狀對於飽和地下水分佈與滲流面的影響。 對於乾燥土壤給予定流量邊界條件,容易發生數值不穩定問題,本研究比較測試三種入滲邊界數值處理方法:(1)以質量平衡原則迭代水頭邊界條件取代流量邊界條件的方法,(2)利用解析解作為迭代求解水頭的初估值,與(3)使用邊界半差分元素控制體積質量守恆法。研究結果發現,方法(1)在所有狀況下,都可以達到收斂與質量守恆,唯迭代計算需花費較多的時間;方法(2)必須利用土壤保水和導水試驗資料推估解析法專用的特性曲線,可以改善收斂性,唯效率與土壤特性曲線與解析法專用的特性曲線的類似程度有關;方法(3)在大部分狀況下可以收斂,也是三種方法中效率最高者,唯在極端狀況下有不能收斂的可能性。 當定流量上邊界大於土壤飽和水力傳導係數時,初期可全部入滲,當地表土層飽和後開始積水,由於外部邊界條件為給定流量邊界條件,在地表土壤飽和後無法以外部給定的定流量入滲。在此狀況下,本研究採用的邊界條件處理方法,是將定流量上邊界轉換為定水頭,在不考慮逕流下,以積水深度與入滲量的和必須等於降雨量(定流量)的方式,迭代決定積水深度(定水頭邊界條件),以及模式是否達到收斂標準。式模擬試驗證實,在重新處理過邊界條件後,不同降雨條件下,均能夠有效率且確定模式收斂的情形下完成模擬。 在空氣可以自由從側邊排除假設下的一維變飽和數值模擬問題中,地表飽和的入滲濕鋒段與地下水位接觸時,土層瞬間飽和將使土壤水壓於短時間內發生急劇變化,可能導致數值不穩定。本研究對於此問題的數值處理方法,是採用土層達到完全飽和前的入滲率,與飽和前後全域水質量差異,推估時間差分步幅的方式處理,可以滿足質量守恆的要求,並且跳過全域壓力水頭應變數大幅驟變的數值處理困難。 將理查氏方程式比擬為運移-擴散微分方程式,由運移-擴散方程式的係數u/D比值,計算佩克勒數(Peclet number),Pe,數值試驗結果顯示,當Pe值越大時,數值擴散情形越明顯。本研究利用多層次細化網格增加模式計算效率的做法,是選擇Pe值門檻,當Pe的數值大於某個增加網格的Pe值門檻時,以網格尺寸減半的比率進行網格細化、降低Pe值;(在濕鋒後端)當Pe的數值小於某個刪除網格的Pe值門檻時,刪除局部最細網格,使網格尺寸倍增減少不必要的計算。利用多層次網格模式自動產生與去除、並且相鄰網格大小比率不超過2:1和1:2範圍的設計,可以使模擬區域的Pe值均維持在一定的數值以下。數值模擬顯示,多層次網格模式與全域採用細網格模擬的含水量及壓力水頭都幾乎重合,局部網格細化所得的數值準確性與全域採用細網格的結果相當。多層次網格細化法,除了可減少數值擴散問題,改善入滲濕鋒位置與梯度的模擬外;也增加網格配置上的彈性,利用數量較少的網格,達成全域使用細網格模擬相同的計算精度,縮短計算所需時間,提升模式的執行效率;同時,亦可以適用於分層異質土壤問題。 由於多層次網格細化模式建立在質量守恆模式的基礎上,數值模擬可以達到完全的質量守恆;但是在增減網格時,因為保水曲線曲率的凹凸,可能增加或減少水量,造成質量不守恆問題。另外,由於濕鋒前端預備細化網格的配置是採用已知時刻的Pe值進行研判,因此可能在模擬下一時間步的水流時,濕鋒範圍超過細網格範圍的問題,尚待進一步解決。 擬三維數值模式之處理方式是擴充一維模式,加上沿坡面方向的坡向維度,以及與坡面正交水平方向的坡寬變化,以考量寬度變化的質量守恆式搭配達西定律重新推導得擬三維坡地地下水流方程式。在邊界的處理上則多了滲流面邊界的迭代求解過程。求解單次迭代過程中線性矩陣的方式是採用預設式共軛梯度法,能有效的減少運算資源,提高計算效率。由擬三維模式與實驗資料的驗證結果可發現,擬三維模式對於二維與楔形實驗土槽之土壤水分與滲流面模擬具有相當不錯的掌握能力。   利用擬三維數值模式模擬穩態與非穩態降雨和坡地滲流,結果顯示,坡地滲流問題必須採用能模擬未飽和層水的流動與勢能分布的變飽和度數值模式,才能正確銜接土壤-大氣接觸的邊界條件與飽和地下水位之間的物理作用機制,對於滲流面分佈的估計具有相當重要的影響。穩態數值試驗顯示,採用杜布假設、以飽和地下水流模式模擬坡地水流,無法正確表現近下游出口處的地下水位曲率和滲流面,並且可能低估實際流量。在非穩態坡地水流問題中,必須具備未飽和土壤水流數值模式,才能描述水自地表入滲到進入飽和地下水面的流動行為,和降雨停止後的退水行為;說明飽和地下水流模式對於坡地水文行為描述能力有限。   由於未飽和層水流主要在垂直方向移動,飽和地下水則橫向移動顯著,欲完整描述坡地之水文反應,坡地水流模式至少需具備垂直二維的模擬能力。擬三維模式利用準單向(沿坡地方向)流動理論與質量平衡觀念,推導出具有束縮或擴張因子的擬三維變飽和度水流控制方程式,可配合應用於卡氏座標與圓柱座標系統問題。透過合理的垂直二維與寬度函數結合,模式可應用於三維坡地之水流模擬,並可節省模式運算資源。   坡面的形狀影響地下水位在坡地土壤中的變化,向下束縮的坡面容易出現滲流面,同時在靠近坡面下游處,單寬流量和壓力水頭梯度較大;向下擴展的坡面則較少出現滲流面,在坡面下游擴展處壓力水頭梯度較小,上游處寬度較窄,單寬流量和壓力水頭梯度較大。坡面上下游水深差異亦會影響滲流面的高度,相同的坡地長度與土壤厚度,當上下游水深差異較大時,下游較容易出現滲流面,即使在寬度向下游擴展的坡面亦是。   坡面飽和水力傳導係數會影響入滲量的多寡,在降雨強度小於飽和水力傳導係數時,降雨可以全部入滲,反之可能只部分入滲。坡面飽和水力傳導係數較降雨強度小者,短延時降雨的累積入滲量不多,加上入滲速度較慢,因此形成飽和地下水位的時間也較慢,只能形成淺薄的飽和地下水層,過程中土壤有效應力均保持較大,坡地較為穩定。飽和水力傳導係數較大的坡面,降雨會快速入滲,雨停後也能夠快速的排水。在降雨強度大、降雨延時比入滲鋒抵達飽和地下水位時間長的狀況下,坡面可能出現滲流面。

並列摘要


The purpose of this study is to develop a Finite-Analytic numerical model to calculate moisture and heat fluxes across the soil-atmosphere interface and the vertical distribution of soil moisture and temperature for the bare soil in the unsaturated zone. With more accurate flux estimates, micrometeorology related and natural ground water recharge predictions are likely to be improved. The study contains the following features: (a) It utilizes the equations for moisture and heat transfer following the works of Chen (1995), and reforms the heat flux at the surface boundary condition for evaporation. (b) It accounts for the temperature effect in the metric head and hydraulic conductivity and uses the experimental results of Hopmans and Dane (1985) to calculate the metric head and the water capacity as functions of temperature. (c) It uses the method of de Vries (1975) to calculate the thermal conductivity as well as the correction factor for the average macroscopic temperature gradient. Both are decided by water content and the volumetric fraction of quartz and clay. (d) It applies a Finite-Analytic Numerical Model for efficient solution of strongly non-linear problems for abrupt changes in atmosphere and soil conditions. (e) The results are tested against field measurements, Jackson (1973) and Jackson et al. (1973). Including surface energy-balance components, water content and soil temperature. According to the results, the simulated and observational trends of the soil water content are difference at begging. With few hours ago, the effect of initial condition is reductive and than the simulated and observational trends of the soil water content is closer to closer. The simulated values of the soil temperature also can response the reasonable physical phenomenon. In the process to calculate the metric head and soil temperature, the convergence and the calculating time of metric head is harder and longer than soil temperature. Otherwise, hydraulic conductivity, soil surface roughness length, and sky cover will influence the varies of soil moisture and soil temperature. The hydraulic conductivity controls and regulates the water movement upward to the surface or downward to the groundwater table where the evaporation or infiltration occurs. The surface roughness length controls the value of mass transfer and sensible heat transfer coefficients, which affects the evaporation rate, latent heat fluxes, and sensible heat flux. The sky cover determines the quantity of solar energy, which can produce significant differences in water content, soil temperature and evaporation rate.

參考文獻


70. Vauclin, M., Vachaud, G. and Khanji, J., “Two dimensional numerical analysis of transient water transfer in saturated-unsaturated soils”. In: G. C. Vansteenkiste (Editor), Modeling and Simulation of Water Resources Systems, North-holland, Amsterdam, pp. 299-323., 1975.
9. Allen, M. B. and Murphy, C. L., “A finite element collocation method for variably saturated flows in porous media”, Numer. Methods Partial Differential Equations, 1, 229-239, 1985.
10. Allen, M. B. and Murphy, C. L., “A finite element collocation method for variably saturated flow in two space dimensions” Water Resour. Res., 22, 1537-1542, 1986.
11. Aly, I. Ei-Kadi. and Ling.G., “The Courant and Peclet number criteria for the numerical solution of Richards’ equation”, Water Resour. Res, 29(10), 3485-3494, 1993.
12. Anderson, M.G. and Burt, T.P., Process Studies in Hillslope Hydrology, John-Wiley & Sons. 1990.

被引用紀錄


林芳慈(2007)。土層瞬間飽和上揚力效應之初步研究〔碩士論文,國立臺北科技大學〕。華藝線上圖書館。https://doi.org/10.6841%2fNTUT.2007.00388
羅瑞祥(2018)。以坡地地下水模擬為目的之水文過程演算法與參數化研究〔碩士論文,國立臺灣大學〕。華藝線上圖書館。https://doi.org/10.6342%2fNTU201800458
張絜晰(2016)。前期降雨對土砂災害預警系統之影響〔碩士論文,國立臺灣大學〕。華藝線上圖書館。https://doi.org/10.6342%2fNTU201603293
洪政暘(2009)。水力緩衝區之研究及其應用〔碩士論文,國立臺灣大學〕。華藝線上圖書館。https://doi.org/10.6342%2fNTU.2009.00178
洪國展(2007)。邊坡埋設排水管對於地下水位差異分析之水文模式改進研究〔博士論文,國立臺灣大學〕。華藝線上圖書館。https://doi.org/10.6342%2fNTU.2007.01341

延伸閱讀