臺灣博碩士論文加值系統

由於台灣地理及氣象條件特殊，降雨型態分布不均及地形陡峭等因素，許多降雨無法被有效的利用即流入大海中，且台灣部分區域降雨量貧乏，使養殖業、農業用水需求量增加，超抽地下水，導致地下水位降低、海水入侵及地層下陷等問題。如何有效的執行地下水補注成為重要的議題。地下水補注為非飽和土壤滲流的問題，故本研究使用有限元素軟體ABAQUS來模擬人工補注井對一區域之蓄水層進行地下水補注之分析，將討論土壤不同的相對密度、補注井之深度及van Genuchten model參數α與n對地表高程變化及補注量之影響。由分析結果顯示蓄水層土壤於不同相對密度(30%、60%及90%)時，以相對密度30%的補注量最大，且蓄水層完成地下水補注所需之時間愈短。對地表隆起情形而言，蓄水層在愈疏鬆的情況下隆起量愈大，反之則愈小。在補注井深度影響方面，補注井深度愈深者，蓄水層完成地下水補注所需之時間愈短，且地表隆起量愈大，但對補注量則無明顯差異。在VG model之參數影響方面，當α值愈小時，補注量愈多，地表隆起量愈大，但蓄水層完成地下水補注所需之時間愈長；當n值愈小時，補注量愈多，地表隆起量愈大，對蓄水層完成補注所需的時間愈久。對於地表隆起量而言，當參數n減少40%時所產生之隆起量比參數α減少40%時來的明顯。

The distribution of rainfall is not uniform in Taiwan due to geographical and climatic condition. It happens that a lot of underused rainwater flows into the sea. When the consumption of water for aquaculture and agriculture increases, groundwater can be ultra-pumping. It will be the reason for reducing groundwater level, invading seawater, and soil settlement. Therefore, how to perform groundwater recharge effectively is an important issue. It is a coupled mechanical-seepage flow analysis in unsaturated soil for groundwater recharge. ABAQUS software was used to simulate a recharge well in a layer of aquifer. Parametric studies of recharge effectiveness were evaluated in this study. The influencing parameters include relative density of aquifer, the depth of recharge well and the values of unsaturated soil parameter α and n of van Genuchten model. The effects of these parameters on the amount of recharge and ground surface uplift were studied. The result shows that the amount of recharge increases and the duration of recharge decreases when relative density of aquifer reduces. The amount of uplift of ground surface increases due to the decrease of relative density. The duration of recharge is shorter and the amount of uplift of ground surface is large when the depth of recharge well is deeper. In unsaturated soil parameter of van Genuchten model, the amount of recharge and ground surface uplift are larger, the duration of recharge is longer when the values of α or n are small. The effect of the value of n on ground surface uplift is more than that of the value of α.

目錄
摘要..............................................I
Abstract........................................II
致謝............................................III
目錄.............................................IV
圖目錄........................................... VI
表目錄........................................... XI
第一章 緒論.......................................1
1.1 研究動機......................................1
1.2 研究目的......................................2
1.3 研究方法......................................3
1.4 研究內容......................................3
第二章 文獻回顧....................................6
2.1 非飽和土壤之定義...............................6
2.2 土壤水分特徵曲線...............................6
2.3 地下水補注....................................8
第三章 數值分析方法與驗證...........................12
3.1 ABAQUS......................................12
3.1.1 前處理(Preprocessing)...................... 13
3.1.2 模擬計算(Simulation)....................... 13
3.1.3 後處理(Post processing).................... 13
3.2 數值分析方法之驗證.............................14
3.2.1 砂箱模型試驗與ABAQUS數值模擬分析...............14
3.2.2 砂箱數值分析結果.............................14
3.3 補注井之數值分析方法............................15
3.3.1 參數設定及說明...............................15
3.3.2 時間增量....................................19
3.3.3 重力場及初始條件設定..........................19
3.3.4 元素種類及網格大小設定.........................21
3.3.5 補注井邊界條件之設定...........................21
第四章 數值分析結果與討論.............................40
4.1 不同相對密度之影響...............................40
4.2 不同補注井深度之影響.............................42
4.3 van Genuchten Model參數α與n之影響...............43
第五章 結論與建議....................................79
5.1 結論...........................................79
5.2 建議...........................................80
參考文獻............................................81

圖目錄
圖1-1 研究流程圖.....................................5
圖2-1 非飽和土壤組成示意圖.............................11
圖2-2 土壤水分特徵曲線各區段說明圖......................11
圖3-1 ABAQUS分析流程................................26
圖3-2 砂箱滲流試驗...................................27
圖3-3 砂箱滲流試驗示意圖..............................27
圖3-4 砂箱試驗尺寸示意圖..............................28
圖3-5 標準砂之水分特徵曲線圖...........................28
圖3-6(a) 砂箱滲流試驗(30sec，滲流速度).................29
圖3-6(b) 砂箱滲流試驗(30sec，飽和度)...................29
圖3-7(a) 砂箱滲流試驗(60sec，滲流速度).................30
圖3-7(b) 砂箱滲流試驗(60sec，飽和度)...................30
圖3-8(a) 砂箱滲流試驗(120sec，滲流速度).................31
圖3-8(b) 砂箱滲流試驗(120sec，飽和度)..................31
圖3-9(a) 砂箱滲流試驗(240sec，滲流速度)................32
圖3-9(b) 砂箱滲流試驗(240sec，飽和度)..................32
圖3-10(a) 砂箱滲流試驗(480sec，滲流速度)...............33
圖3-10(b) 砂箱滲流試驗(480sec，飽和度).................33
圖3-11(a) 砂箱滲流試驗(600sec，滲流速度)...............34
圖3-11(b) 砂箱滲流試驗(600sec，飽和度).................34
圖3-12(a) 砂箱滲流試驗(720sec，滲流速度)...............35
圖3-12(b) 砂箱滲流試驗(720sec，飽和度).................35
圖3-13(a) 砂箱滲流試驗(900sec，滲流速度)...............36
圖3-13(b) 砂箱滲流試驗(900sec，飽和度).................36
圖3-14 補注井尺寸示意圖...............................37
圖3-15 補注井邊界條件示意圖............................37
圖3-16 補注井補注位置示意圖............................38
圖3-17 不同參數α之水分特徵曲線圖(Dr=60%，n=2.9).........38
圖3-18 不同參數n之水分特徵曲線圖(Dr=60%，α=2.2).........39
圖4-1 不同相對密度之總補注量比較圖(α=2.2，n=2.9).........47
圖4-2 不同相對密度之飽和度比較圖(α=2.2，n=2.9)...........47
圖4-3 不同相對密度之地表垂直位移比較圖(12hr).............48
圖4-4 不同相對密度之地表垂直位移比較圖(24hr).............48
圖4-5 不同相對密度之地表垂直位移比較圖(10day)............49
圖4-6 不同相對密度之地表垂直位移比較圖(30day)............49
圖4-7 不同相對密度之地表垂直位移比較圖(60day)............50
圖4-8 不同相對密度之地表垂直位移比較圖(120day)...........50
圖4-9 不同補注井深度之總補注量比較圖(Dr=60%，α=2.2，n=2.9)...51
圖4-10 不同補注井深度之飽和度比較圖(Dr=60%，α=2.2，n=2.9)...51
圖4-11 不同補注井深度之地表垂直位移比較圖(12hr).............52
圖4-12 不同補注井深度之地表垂直位移比較圖(24hr).............52
圖4-13 不同補注井深度之地表垂直位移比較圖(10day)............53
圖4-14 不同補注井深度之地表垂直位移比較圖(30day)............53
圖4-15 不同補注井深度之地表垂直位移比較圖(60day)............54
圖4-16 不同補注井深度之地表垂直位移比較圖(120day)...........54
圖4-17 不同參數α之總補注量比較圖(Dr=60%，n=2.9)............55
圖4-18 不同參數α之飽和度比較圖(Dr=60%，n=2.9).............55
圖4-19 不同參數α之地表垂直位移比較圖(12hr).................56
圖4-20 不同參數α之地表垂直位移比較圖(24hr).................56
圖4-21 不同參數α之地表垂直位移比較圖(10day)................57
圖4-22 不同參數α之地表垂直位移比較圖(30day)................57
圖4-23 不同參數α之地表垂直位移比較圖(60day)................58
圖4-24 不同參數α之地表垂直位移比較圖(120day)...............58
圖4-25 不同時間下之地表垂直位移圖(參數α=1.32)..............59
圖4-26 不同時間下之地表垂直位移圖(參數α=3.08)..............59
圖4-27 不同參數n之總補注量比較圖(Dr=60%，α=2.2)...........60
圖4-28 不同參數n之飽和度比較圖(Dr=60%，α=2.2).............60
圖4-29 不同參數n之地表垂直位移比較圖(12hr).................61
圖4-30 不同參數n之地表垂直位移比較圖(24hr).................61
圖4-31 不同參數n之地表垂直位移比較圖(10day)................62
圖4-32 不同參數n之地表垂直位移比較圖(30day)................62
圖4-33 不同參數n之地表垂直位移比較圖(60day)................63
圖4-34 不同參數n之地表垂直位移比較圖(120day)...............63
圖4-35 不同時間下之地表垂直位移圖(參數n=1.74)..............64
圖4-36 不同時間下之地表垂直位移圖(參數n=4.06)..............64
圖4-37 相對密度30%之孔隙水壓分布圖(a)初始狀態(b)24小時(c)10天...65
圖4-37 相對密度30%之孔隙水壓分布圖(d)30天(e)60天(f)120天......66
圖4-38 相對密度60%之孔隙水壓分布圖(a)初始狀態(b)24小時(c)10天...67
圖4-38 相對密度60%之孔隙水壓分布圖(d)30天(e)60天(f)120天......68
圖4-39 相對密度90%之孔隙水壓分布圖(a)初始狀態(b)24小時(c)10天...69
圖4-39 相對密度90%之孔隙水壓分布圖(d)30天(e)60天(f)120天......70
圖4-40 補注井深度11m之孔隙水壓分布圖(a)初始狀態(b)24小時(c)10天...71
圖4-40 補注井深度11m之孔隙水壓分布圖(d)30天(e)60天(f)120天.....72
圖4-41 補注井深度13m之孔隙水壓分布圖(a)初始狀態(b)24小時(c)10天...73
圖4-41 補注井深度13m之孔隙水壓分布圖(d)30天(e)60天(f)120天.....74
圖4-42 參數α=1.32之孔隙水壓分布圖(a)初始狀態(b)60天(c)120天....75
圖4-43 參數α=3.08之孔隙水壓分布圖(a)初始狀態(b)60天(c)120天....76
圖4-44 參數n=1.74之孔隙水壓分布圖(a)初始狀態(b)60天(c)120天....77
圖4-45 參數n=4.06之孔隙水壓分布圖(a)初始狀態(b)60天(c)120天....78
 
表目錄
表3-1 常用的單位系統.......................................23
表3-2 材料參數設定(砂箱)................................... 23
表3-3 材料參數設定(補注井)..................................24
表3-4 不同條件下之參數說明.................................. 25
表4-1 各參數分析結果一覽表.................................. 46

[1]Fredlund, D. G., and Rahardjo, H., Soil Mechanics for Unsaturated Soils, Wiley-Interscience Publication, 1993.
[2]van Genuchten, M. T., “A closed-form equation for predicting the hydraulic conductivity of unsaturated soils,” Soil Sci. Soc. Am. J., 44, pp.892-898, 1980.
[3]Fredlund, D. G., and Xing, A., “Equations for the soil-water characteristic curve,” Canadian Geotechnical Journal, 31(3), pp. 521-532, 1994.
[4]Harrison, B. A. and Blight, G. E., “The use of indicator tests to estimate the drying leg of the soil-water characteristic curve,” Asian Conference on Unsaturation Soil, Balkema, Rotterdam, 2000.
[5]Sillers, W. S., Fredlund, D. G., and Zakerzadeh, N. “Mathematical attributes of some soil-water characteristic curve models,” Geotech. Geol. Eng., 19(3-4), pp.243-283, 2001.
[6]Das, B. M., Principles of Geotechnical Engineering, Sixth Edition, PWS, Boston, Mass, 1985.
[7]Itasca consulting group, inc. fast lagrangian analysis of continua, version 5.0, user’s manual, Minneapolis, Minnesota, USA, 2005.
[8]Das, B. M., Advanced Soil Mechanics, Second Edition, Taylor and Francis, 1997.
[9]Meyerhof, G. G. “Compaction of sands and the bearing capacity of piles,” JSMFD. ASCE, 85 (SM6), pp.1-29, 1959.
[10]Skempton, A. W., “Standard penetration test procedures and the effects in sands of overburden pressure, Relative Density, Particle Size, Ageing and Over Consolidation,” Geotechnique 36(3), pp.425-447, 1986.
[11]Schmertmann, J. H., “Static Cone to Compute Settlement Over Sand,” Journal of the Soil Mechanics and Foundations Division, American Society of Civil Engineers, Vol. 96, No. SM3, pp.1011-1043, 1970.
[12]Santamarina, J. C. and Cho, G. C., “Determination of critical state parameters in sandy soils－simple procedure,” Geotechnical Testing Journal, GTJODJ, Vol. 24, No. 2, June, pp.185-192, 2001.
[13]劉建榮，「van Genuchten 土壤特徵曲線參數對濕峰模擬與暫態補注量之影響」，碩士論文，私立逢甲大學土木及水利工程研究所，台中，2001。
[14]范嘉程、黃俊龍，「頁岩風化不飽和土壤之土壤水分特徵曲線」，中華水土保持學報，43(3)，第197-205頁，2012。
[15]蕭祝融，「以補注井延緩地層下陷之研究-本省沿海地區之個案分析評估」，碩士論文，私立逢甲大學土木及水利工程研究所，台中，1987。
[16]楊文絮，「田間地下水井補注之研究」，碩士論文，私立逢甲大學土木及水利工程研究所，台中，1988。
[17]詹子輝，「應用補注井群改善海水入侵之可行性研究」，碩士論文，國立台灣大學農業工程學系研究所，台北，1998。
[18]倪春發、許少華，「人工池補注地下水水庫之效益初步模擬與影響因子」，第三屆地下水資源及水質保護研討會論文集，1999。
[19]王銘熙，「應用人工補注法改善海水入侵現象之探討」，碩士論文，私立中原大學土木工程學系，桃園，2000。
[20]翁琮哲，「地下水補注功效之探討」，碩士論文，私立中原大學土木工程研究所，桃園，2001。
[21]莊志威，「地下水庫截水牆分析探討」，碩士論文，私立中原大學土木工程研究所，桃園，2010。
[22]劉致翔，「地下水補注對沿海含水層海水入侵整治之數值分析」，碩士論文，國立成功大學水利及海洋工程學系碩博士班，台南，2011。
[23]戴宏育，「沿海地區海水入侵之數值研究-以屏東平原為例」，碩士論文，國立成功大學水利及海洋工程學系碩博士班，台南，2013。
[24]林建翰，「人工湖補注地下水之力學分析」，碩士論文，私立朝陽科技大學營建工程系碩士班，台中，2014。
[25]費康、張建偉，「ABAQUS在岩土工程中的應用」，中國水利水電出版社，北京市，2010。
[26]曹金風、石亦平，「ABAQUS有限元分析常見問題解答」，機械工業出版社，北京市，2009。
[27]Joseph, J. A. and Jacques H., “Use of the falling-head method to assess permeability of freshly mixed cementitious-based materials,” Civ. Eng., ASCE, 25(5), pp.580-588., 2013.