過往在邊坡穩定性分析的研究中大多是使用二維的數值分析方法以及結果做呈現，應變都以平面狀態作為基礎假設的方式，此方式做出的結果通常會以一剖面上的曲線來表示滑動位置，若想知道整體滑動面積或範圍會較難從結果得知，且與實際邊坡自然的三維崩塌情形難免有所異同。本研究比起以往邊坡穩定的數值分析更具有研究價值的地方在於，相較二維邊坡穩定分析在模擬結果的顯示上更直觀且具合理性。
本研究通過整合GMS與FLAC3D兩套數值軟體，先於建模功能強大的地下水數值軟體GMS中建出三維地層模型，透過MATLAB進行資料轉換後匯入FLAC3D中，以擷取兩套軟體的強項做結合應用的方式協助補足FLAC3D擬真建模能力偏弱的部分，而後再以FLAC3D內建語言FISH撰寫副程式來進行檔案寫入以及後續的邊坡穩定性分析，並以茂林區D048大規模崩塌區作為模擬邊坡的對象，並以鳳凰颱風期間地下水位監測數據做為參考，設定四種地下水位面進行不同情境的模擬，並以安全係數、位移、剪應變增量以及塑性區來對模擬成果進行下水位對於邊坡崩塌影響的探討。
本研究以擬真三維地層模擬結果來探討地下水與邊坡穩定性關係的成果，使模型更貼近現址狀況，並且可以成為未來學者研究的基礎，或是可供政府於災害防治時用作參考。

In this study, by integrating two sets of numerical software, GMS and FLAC3D, a three-dimensional stratigraphic model was built in the groundwater numerical software GMS with powerful modeling capabilities, and the data was converted through MATLAB and then imported into FLAC3D to capture the strengths of the two sets of software. Combined with the application, it helps to make up for the weak part of FLAC3D's immersive modeling ability, and then write the sub-program in FLAC3D's built-in language FISH for file writing and subsequent slope stability analysis. Taking the large-scale collapse area of D048 in Maolin District as the object of slope simulation, the simulation results are used to discuss the impact of water level on slope collapse. This research can be used as a basis for scholars' research in the future, or can be used as a reference for the government in disaster prevention and control.

1. Aquaveo, 2020, GMS User Manual (v10.5).
2. Brown, .B.&King P. ,1966, Automatic embankment analysis:equilibrium and instability conditions”, Geotechnique, 16 (3), pp.209-219.
3. Cheng, Y. M. & Yip, C. J. , 2007, “Three-Dimensional asymmetrical slope stability analysis extension of Bishop’s, Janbu’s, and Morgenstern–Price’s techniques”, Journal of geotechnical and geoenvironmental engineering, 133 (12), pp.1544-1555.
4. Griffiths, D. V. & Lane, P. A. ,1999, “Slope stability analysis by finite elements”, Geotechnique, 49 (3), pp.387-403.
5. Salencon, J., 1983, Calcul a la rupture et analyse limite, Paris, Presses ENPC.
6. Varnes D.J. 1978, “Slope movement types and processes, “In: Special Report 176: Landslides: Analysis and Control (Eds: Schuster, R. L. & Krizek, R. J.), Transportation and Road Research Board, National Academy of Science, Washington D. C., pp.11-33.
7. Zhang, L. L., D. G. Fredlund, M. D. Fredlund and G. W. Wilson., 2014, Modeling the Unsaturated Soil Zone in Slope Stability Analysis., Canadian Geotechnical Journal 51 (12) ,pp.1384-1398.
8. Zienkiewicz, O.C. and Cheng Y.K., 1967, The finite element method in structural and continuum mechanics: numerical solution of problems in structural and continuum mechanics, McGraw-Hill, London.
9. 行政院農業委員會水土保持局，2014，「萬山潛在大規模崩塌區調查監測計畫成果報告書」。
10. 行政院農業委員會水土保持局，2015，「萬山潛在大規模崩塌地區監測系統為運維與擴充計畫期末報告書」。
11. 行政院農業委員會水土保持局，2019，「高雄市茂林區 D048 大規模崩塌區調查及多納溫泉溪集水區治理規劃期末報告書」。
12. 邱子芸，2017，應用Scoops3D於潛在大規模崩塌影響之研究，國立台灣海洋大學河海工程學系碩士論文。
13. 吳憶珊，2018，應用三維數值模擬技術探討地下水位與邊坡崩塌之研究，國立台灣海洋大學河海工程學系碩士論文。
14. 李怡萱，2020，以三維擬真地層模擬探討降雨入滲對邊坡穩定性影響之研究，國立成功大學資源工程學系碩士論文。
15. 林芷瑩，2006，剪力強度折減法應用於地震邊坡之研究，國立成功大學資源工程學系碩士論文。
16. 周宗仁，2020，邊界元素法在海岸工程應用，http://1.34.137.146/bem_html/bem.html
17. 美國地質調查局USGS(2014),” Landslide Types and Processes”,https://pubs.usgs.gov/fs/2004/3072/fs-2004-3072.html.
18. 陳榮河、洪勇善，1999，崩塌地整治工法之介紹，地工技術雜誌，第72 期，5-12。
19. 陳志昌，2001，FLAC程式應用於土壤邊坡穩定分析，碩士論文，國立中央大學應用地質研究所。
20. 陳育民、徐鼎平，2009，FLAC/FLAC3D基礎與工程實例，中國水利水電出版社。
21. 陳容花，2015，降雨條件下入滲對不飽和土壤邊坡之影響及臨界體積含水比探討，國立台灣科技大學營建工程系碩士論文。
22. 陳容花，2015，降雨條件下入滲對不飽和土壤邊坡之影響及臨界體積含水比探討，國立台灣科技大學營建工程系碩士論文。
23. 陳欣郁，2019，以地層模擬結果進行三維度邊坡穩定分析之研究，國立成功大學資源工程學系碩士論文。
24. 曾惠蒂，2005，變形分析法應用於地震邊坡之研究，國立成功大學資源工程學系碩士論文。
25. 馮正一、譚學斌，2015，以邊坡物理模型試驗探討降雨與地下水位上升引起之崩塌，社團法人中華水土保持學會104 年度會員代表暨會員大會─大會手冊暨研討會。
26. 黃喬雅，2016，「坡地崩塌知多少？」，台灣颱風洪水研究中心颱洪科普網，http://www.ttfri.narl.org.tw/sp/。
27. 經濟部中央地質調查所，2012，莫拉克災區潛在大規模崩塌地區分析報告，中央災害防救會報。
28. 鍾明劍、陳建新、李璟芳、黃韋凱、譚志豪，2014，潛在大規模崩塌地對鄰近聚落衝擊之調查與評估方法，中華水土保持學報，第47卷第3期： 122-134。