在以往邊坡穩定分析中通常無法完全顯示與現地狀況近乎相同的起伏、層數、各層厚度及地下水位之分佈，而經常以單一或數個縱斷面的型式來表示。惟邊坡破壞之寬度可能連續至數十或數百公尺寬，其深度與各縱斷面也不近相同，因此以單位長度來計算滑動土方之數量也會出現很大的誤差；另若位處泥岩地質除其具有親水性强、透水性弱，且遇水則易軟化、塑性變形、抗風化能力弱與崩解等特性外，其受地下水位之影響也相當顯著，最重要的是它為臺灣中南部地區常見之地質，亦有坡地之癌的稱號。
本研究提出一套三維數值分析方法，過程中結合擬真地層模擬與剪力強度折減法探討不同降雨強度下之擬真三維地下水面，進而分析邊坡整體安全係數之分佈與其可能潛在滑動面與滑動量，再利用反饋分析機制求取邊坡臨界與警戒累積降雨量。而其做法是將觀測鑽孔所得土層及水位資料在GMS（Groundwater Modeling System）軟體下建立三維擬真數值地層網格，再透過MATLAB（MATrix LABoratory）軟體，撰寫轉換程式，將網格資料轉換成以連體力學為運算之FLAC3D（Fast Lagrangian Analysis of Continua）軟體所能讀取之資料，再配合FLAC3D內建程式轉換資料逕行構建成三維擬真網格模型、水位面和地層分層狀況，並給定相關參數後進行邊坡穩定性分析。
本研究最後再以一公路邊坡為案例進行分析，選取為台南市道175線25K+410泥岩滑動區，並利用臺南市政府工務局長期對該地區之地質調查﹅邊坡滑動監測﹅地下水位等資料及研究期間中央氣象局的雨量觀測資料，並結合本研究所提出之分析模式與架構，探討不同雨量下的安全係數﹅位移量﹅地下水位﹅剪應變增量及滑動趨勢的三維數據，再利用前述這些數據回歸出該地區需警戒及行動的雨量值與位移量，作為日後該地區的災害預防、預警、復建等參考。

Short-time, heavy rainfall is an important influencing factor for slope slip in sensitive areas. This research combined three dimensional pseudo-stratigraphic simulation, and shear strength reduction method to explore the pseudo -three-dimensional groundwater surface under different rainfall intensity firstly; then analyzed the overall safety factor distribution and its potential sliding surface and amount. Finally, using the feedback analysis mechanism to obtain the slope critical and alert cumulative rainfall. This research applied the GMS software to create a 3D pseudo-numeric stratigraphic grid, and then use MATLAB software to convert the grid data into FLAC3D software can be converted into a 3D realistic mesh model, water level surface and stratified layer status through built-in programs in FLAC3D, and slope stability analysis is performed after relevant parameters are given.
This research take a mudstone slope of a highway as a case and collect relative data ( i.e. the geological survey of the area, slope slip monitoring, groundwater level…),then applied the proposed model to explore the 3 dimension data of safety factor, displacement, groundwater level, shear strain and sliding trend under different rainfall. Then conduct regression analysis predicts the rainfall values and displacements in the area that require alertness and action, as a reference for future disaster prevention, early warning, and reconstruction in the area.
Three-dimensional slope stability analysis process established in this study is closer to the pseudo-stratigraphic conditions, and its numerical simulation results are also consistent with the on-site investigation results and occurrence conditions.

1. Aquaveo (2017), " GMS 10.3".
2. Aquaveo (2017), "GMS User Manual (v10.3) " .
3. Brown, C. B. & King, I. P.(1966), "Automatic embankment analysis: equilibrium and instability conditions", Geotechnique, 16 (3), pp.209-219.
4. Chen, R., H. Chen, K. Chen and H. Zhung (2009),"Simulation of a slope failure induced by rainfall infiltration", Environmental Geology 58(5): pp.943-952.
5. Cruden,D.M. and D. J.Varnes (1996),"Landslide types and processes", Landslides-Investigation and Mitigation, Special Report,247: pp.36-75.
6. Clough, R. W. & Woodward, R. J.(1967), "Analysis of embankment stresses and deformations", Journal of Soil Mechanics & Foundations Div.
7. Dawson, E. M., Roth, W. H. & Drescher, A.(1999),"Slope stability analysis by strength reduction", Geotechnique, 49 (6), pp.835-840.
8. Duncan, J.M. (1996), "State of the Art: Limit Equilibrium and Finite-Element Analysis of Slopes", Journal of Geotechnical Engineering, 122, pp.577-596.
9. Duncan, J. M. & Wright, S. G.(1980)," The accuracy of equilibrium methods of slope stability analysis", Engineering geology, 16 (1-2), pp.5-17.
10. Dingwall, J. C. & Scrivner, F. H.(1954),"Application of the elastic theory to highway embankments by use of difference equations", In Proceedings, Highway Research Board, 33, pp.476-481.
11. Griffiths, D. V. & Lane, P. A.(1999), "Slope stability analysis by finite elements", Geotechnique, 49 (3), pp.387-403.
12. Giam, S. K. & Donald, I. B.(1988),"Determination of critical slip surfaces for slopes via stress-strain calculations", In Fifth Australia-New Zealand Conference on Geomechanics, Institution of Engineers, Australia, pp.461
13. Huang, C. C., Ju, Y. J., Hwu, L. K., and Lee, J. L. (2009),"Internal soil moisture and piezometric responses to rainfall-induced shallowslope failures", Journal of Hydrology, pp.39-51.
14. Hammah, RE, Yacoub, TE, Corkum, B. & Curran (2005),"The Shear Strength Reduction Method for the Gen-eralized Hoek-Brown Criterion", In Proceedings of the 40thUS Symposium on Rock Mechanics, Alaska Rocks 2005,Anchorage, Alaska.
15. Hoek, E. and Bray, J.W. (1981),"Rock Slope Engineering",Revised 3rd Edition, The Institution of Mining and Metallurgy, London, pp. 341-351.
16. ITASCA(2018),"FLAC3D 5.00"。
17. Jinyu, L., Y. Qing and L. Maotian (2010),"Limit equilibrium analysis of the influence of matric suction on unsaturated soil slope stability", Journal of Convergence Information Technology,5(10), pp. 167-174.
18. Krahn, J.(2004), "Stability modeling with Slope/W", An Engineering Methodology, 1 (1), pp.396.
19. Lu, M. and J. Liu (2009),"FLAC 3D Numerical simulation of water content effect on thick soil slope distortion", Journal of Sanming University, 4(25).
20. Lin, M. L. and Hung, J.J. (1982),"The Influence of Moisture Content on Mechanical Properties of some Sedimentary Rocks in Taiwan", Proc. of the 7th southeast Asian Geotechnical Conference, Hong Kang.
21. Ma, T.H., C.J. Li, Z.M. Lu & Q.Y. Bao (2015),"Rainfall intensity–duration thresholds for the initiation of landslides in Zhejiang Province, China",Geomorphology, 245, pp. 193-206.
22. Marcato, G., M. Mantovani, A. Pasuto, L. Zabuski and L. Borgatti (2011),"Monitoring, numerical modelling and hazard mitigation of the Moscardo landslide (Eastern Italian Alps)" , Engineering Geology ,128, pp.95-107.
23. Manzari, M.T. and M.A. Nour(2000)," Significance of soil dilatancy in slope stability analysis" , J. Geotech. Geoenviron. Eng., 126, pp. 75-80.
24. Matsui, T. & San, K. C.(1992),"Finite element slope stability analysis by shear strength reduction technique", Soils and foundations, 32(1), pp.59-70.
25. Rahimi, A., H. Rahardjo and E. C. Leong (2011),"Effect of antecedent rainfall patterns on rainfall-induced slope failure", Journal of Geotechnical and Geoenvironmental Engineering ,137, pp.483.
26. Shukha R. and Baker R. (2003),"Mesh geometry effects on slope stability calculation by FLAC strength reduction method - Linear and non-linear failure criteria", Proc. of the 3rd Intr. Symp. On FLAC & Numerical modeling in Geomechanics. Sudbury, Ont. Canada, pp.109-116.
27. Siddharthan, R.V. & M. EL-Gamal(1998),"Permanent Rotational Deformation of Dry Cohesionless Slopes Under Seismic Excitation", Transportation Research Board of the National Academies, 1633 (361), pp.45-50.
28. Turner, A. Keith and Schuster, Robert L. (1996)," Landslides: Investigation and Mitigation",Transportation Research Board.
29. Terlien, M. T. J.(1996),"Modelling spatial and temporal variations in rainfall-triggered landslides ", Enschede, Netherlands: ITC Publ Press.
30. The MathWorks (2018)," MATLAB 9.5"。
31. Ugai, K.(1989),"A method of calculation of total safety factor of slope by elasto-plastic FEM", Soils and Foundations, 29 (2), pp.190-195.
32. Ugai, K. & Leshchinsky, D. O. V.(1995),"Three-dimensional limit equilibrium and finite element analyses: a comparison of results", Soils and foundations, 35 (4), pp.1-7.
33. Varnes, D. J. (1978),"Slope movement types and processes", Transportation Research Board Special Report, 176, pp.11-33.
34. Wright, S. G., Kulhawy, F. G. & Duncan, J. M.(1973),"Accuracy of equilibrium slope stability analysis", Journal of Soil Mechanics & Foundations Div, pp.99.
35. Zou, J. Z., Williams, D. J. & Xiong, W. L.(1995),"Search for critical slip surfaces based on finite element method", Canadian Geotechnical Journal, 32 (2), pp.233-246.
36. Zienkiewicz, O. C., Humpheson, C. & Lewis, R. W.(1975), "Associated and non-associated visco-plasticity and plasticity in soil mechanics", Geotechnique, 25 (4), pp.671-689.
37. 方世杰、李德河、林宏明、吳秉晃(2001)，「集集地震後阿里山公路邊坡之崩壞特性探討」，2001年公共工程學術研討會，頁21‐34，臺南。
38. 中央地質調查所(2008)，「五萬分之一臺灣區域地質圖」，臺灣：中央地質調查所(網路資料)。
39. 正昇工程顧問(2018)，「台南市道174線50K+500處地錨工程Ⅰ」。
40. 正昇工程顧問(2017)，「台南市道174線50K+500處地錨工程Ⅱ」。
41. 行政院農業委員會水土保持局(2018)，「山坡地土地可利用限度分類標準」，水土保持法規。
42. 向冠瑋(2017)，「利用三維數值技術分析集水井橫向集水管之控制參數及其對邊坡穩定之影響」，國立中興大學水土保持學系碩士論文。
43. 吳憶珊(2018)，「應用三維數值模擬技術探討地下水位與邊坡崩塌之研究」，國立臺灣海洋大學河海工程學系碩士論文。
44. 李思緯(2012)，「基於有限元素法之邊坡穩定分析的模型不確定性」，國立臺灣大學土木工程學研究所碩士論文。
45. 邱子芸(2017)，「應用 Scoops3D 於潛在大規模崩塌影響之研究」，國立臺灣海洋大學河海工程學系碩士論文。
46. 青山工程顧問股份有限公司(2019)，「六甲區南勢坑市道175線25K+410道路邊坡整治第一期工程規劃報告 」。
47. 青山工程顧問(2018)，「台南市道174線50K+500地滑調查及邊坡監測等評估服務工作成果報告-核定版 」。
48. 青山工程顧問股份有限公司(2017)，「市道174線50K+500地滑調查及邊坡監測成果報告 」。
49. 青山工程顧問股份有限公司(2005)，「174線50K+500整體規劃調查報告」。
50. 周南山（2005），「山區道路邊坡災害防治」，森林遊憩設施規劃設計與施工研習會，台灣。
51. 倪勝火(2013)，「台南174縣道地滑整治對策與施工探討」，國立成功大學土木工程學系碩士論文。
52. 展鴻工程公司(1999)，「174 線50K+060 段災害鑽探工程-地質鑽探報告」，交通部公路總局第五區養護工程處。
53. 財團法人中興工程顧問社(2007-2011)，「易淹水地區上游集水區地質調查及資料庫建置(第1期96年度至第3期100年度)-集水區水文地質對坡地穩定性影響之調查評估計畫-計畫成果報告」，經濟部中央地質調查所研究報告。
54. 陳欣郁(2019)，「以地層模擬結果進行三維度邊坡穩定分析之研究」，國立成功大學資源工程學系碩士論文。
55. 陳昭旭(2019)，「南勢坑大規模崩塌潛勢區邊坡行為-子計畫四：三維度邊坡穩定分析技術之研究計畫」。
56. 陳彥澄(2010)，「應用彈塑性模型模擬地震引致邊坡破壞之滑動量」，國立臺灣大學土木工程學研究所碩士論文。
57. 陳志昌(2000)，「FLAC 程式應用於土壤邊坡穩定分析」，國立中央大學應用地質研究所碩士論文。
58. 陳宏宇（1998），「山崩」，地球科學園地，第6期，第12-21頁。
59. 黃博寅、李德河、林宏明 (2017)，「降雨量與地下水位變動對南部山區邊坡穩定性影響之研究」，中華防災學刊，第九卷第二期，頁133-146。
60. 黃博寅(2016)，「降雨量與地下水位變動對南部地區邊坡穩定性影響之研究」，國立成功大學土木工程學系碩士論文。
61. 黃柏源(2014)，「降雨條件下不飽和土壤邊坡穩定性之參數探討」，國立臺北科技大學土木與防災學系碩士論文。
62. 黃亞雯(2013)，「考慮地下水位變化之擬似三維邊坡穩定模式之建立」，國立交通大學土木工程學系碩士論文。
63. 黃克紹(2005)，「考量大地材料變異性之三度空間邊坡穩定分析」，國立成功大學土木工程學系碩士論文。
64. 黃健政、李錦發、紀宗吉(2002) ，「坡地防災技術及資料庫建置」，經濟部中央地質調查所。
65. 張郁涵(2018)，「邊坡穩定分析方法比較及影響因子探討」，國立中興大學水土保持學系碩士論文。
66. 張家騰(2013)，「台南174縣道地滑整治對策與施工探討」，國立成功大學土木工程學系碩士論文。
67. 張光宗、林哲儀(2008)，「以極限平衡法分析邊坡穩定-岩坡及崩積層邊坡為例」，中華水土保持學報，第三十九卷第二期，頁155-165。
68. 張國欽(2005)，「邊坡在降雨入滲狀況下之穏定性分析與評估」，國立中興大學水土保持學系碩士論文。
69. 張元良(2004)，「以剪力強度折減法進行邊坡穩定分析之研究」，國立中興大學土木工程研究所博士論文。
70. 郭仁立(2005)，「有限元素方法於三維度土壤邊坡穩定分析之應用」，國立中興大學土木工程學系碩士論文。
71. 傅桂霖、吳瑞鵬、黃祥慶(2014)，「泥岩區地滑監測與分析之探討」，水土保持學報，第四十六卷第三期，頁1079-1092。
72. 曾致仁、李德河、吳建宏、林宏明 (2018)，「降雨引致邊坡地下水位變動及作為邊坡崩壞預警指標之研究」，中華防災學刊，第十卷第一期，頁17-29。
73. 曾國鈞(2015)，「降雨引致地下水位變動對公路邊坡穩定影響之研究」，國立成功大學土木工程學系碩士論文。
74. 曾子育(2011)，「公路邊坡破壞與颱風雨型關係之研究-以台南縣道174、175為例」，國立成功大學地球科學系碩士論文。
75. 楊政憲(2000)，「有限差分法應用於受地震力邊坡之變分穩定分析硏究」，臺灣大學土木工程學研究所碩士論文。
76. 王濤、韓煊、趙先宇、朱永生(2015)，「數值模擬方法及工程應用-深入剖析FLAC 3D 5.0」，中國建築工業出版社。
77. 廖洪鈞、廖瑞堂(1999)，「坡地社區開發安全監測手冊」。
78. 蔡宜長(1989)，「應用有限元素法探討滲流對邊坡穩定性的影響」，國立中興大學土木工程研究所碩士論文。
79. 劉小偉、劉高、諶文武、韓文峰(2003)，「降雨對邊坡變形壞影響的综合分析」，岩石力學與工程學報，第二十二卷增二，頁2715-2718。
80. 聯合大地工程顧問公司(2007)，「曾文水庫越域引水計畫：烏山嶺隧道檢討及調查規劃地質調查試驗計畫-地質初步探查試驗」，經濟部水利署水利規劃試驗所。
81. 譚志豪、鍾明劍、許世孟、蘇泰維、陳勉銘、費立沅(2011)，「山崩之地質與水文特性」，2011流域地質與坡地災害研討會論文集，第101-118頁。
82. 鍾明劍、譚志豪、陳勉銘、蘇泰維(2013)，「以定率法評估邊坡山崩臨界雨量-以南勢坑為例」，中華水土保持學報，第四十四卷第一期，頁66-77。
83. 羅佳明(2013)，「蘇花公路115.9 K邊坡調查穩定分析及運移模擬之研究」，建國科大理工期刊，第三十三卷第一期，頁1-23。
84. 蘇郁婷(2013)，「利用滲透試驗探討雙層土壤之滲流沖蝕行為」，國立臺灣大學土木工程學研究所碩士論文。
85. 蘇冠憲(2007)，「砂泥岩互層邊坡滑動機制及保護對策探討-以縣道174線公路50k+650為例」，立德管理學院資源環境系碩士論文。