Title

分離元素法山崩模型建置及震滑機制行為初探

Translated Titles

Landslide model set up and simulation using discrete element method

Authors

翁晨桓

Key Words

分離元素法 ; 數值地形模型 ; 草嶺山崩 ; 集集地震 ; Discrete element method ; Digital terrain model ; Tsaoling landslide ; Chi-Chi earthquake

PublicationName

臺北科技大學土木與防災研究所學位論文

Volume or Term/Year and Month of Publication

2011年

Academic Degree Category

碩士

Advisor

張國楨

Content Language

繁體中文

Chinese Abstract

本研究以1999年集集地震所引發之草嶺山崩事件為例,運用數值地形模型(DTM)資訊,結合顆粒流分離元素之PFC 3D程式分析山崩之行為及機制。為決定山崩塊體及山崩影響範圍,分別運用災前、災後DTM來分析,並用以規範崩塌區之上、下部地形面,定義其間所夾空間即為崩塌之區塊。首先,本研究嘗試數種不同方法來建置崩塌區塊體及模型顆粒叢集等數值模擬之準備,包含NO SHADOW指令法、霣降法、三角柱法、重力排除法、顆粒汰選法等,比較各方法之優缺點及模型可行性等。本研究結果指出顆粒汰選法為最佳方法,其方法為依模擬岩石試體製作流程預先製作一含蓋整個崩塌塊體之顆粒匯集區,再將顆粒予以排除或汰選,配合施予顆粒元素間鍵結作用,進而建立一接近實際外幾何之形態、厚度之崩塌區塊,滑動面亦由DTM資料來建立,符合真實地形特徵之地表模型。 本研究考慮草嶺強震站之地震加速度資料,將之轉換為給定模型之地表地動速度,並考慮岩層中存在孔隙水壓所造成之影響,另考量岩層因受震產生高速滑移,並引致岩層摩擦係數弱化之影響。考量以上因素後,最後以簡化方式一併列入模擬之不同參數及相關模擬條件。 由模擬結果數據資料建置之DTM與真實DTM比對顯示岩層弱化後摩擦係數為0.087時可得最接近真實之堆積型態,其堆積區高程差平均值為-3.2m,樣本標準差為17.5。山崩發生之最大滑動距離為3662 m,最大落距高度為632 m,最高速度為87.2 m/s,本速度峰值之時間序列和地震訊號所得結果吻合。 研究成果顯示以三維離散元素顆粒流數值方法可模擬山崩發生之完整動態歷程,相較於一般常用之極限平衡法及有限元素法,除可求得破壞面及應力應變得知穩定性外,該方法還可模擬山崩發生之大變形行為,包含裂縫、運動軌跡、速度變化、滑動距離,乃至三維的堆積型態,進一步掌握災害發生之影響範圍。研究成果指示了本模擬之成效,以及震滑相關現象之推論。

English Abstract

This research takes Tsaoling landslide caused by Chi Chi Earthquake as an example, adopts from DTM data and also integrates the PFC 3D program of Particle flow discrete element method to analyze the sling mechanism and the landslide behaviors. To determine the affected area of landslide and to constraint the slid block, we use pre-disaster DTM and post-event DTM to build two topographic surfaces (upper topographic surface and lower topographic surface), thus the slid block between upper topographic surface and lower topographic surface could be defined. First of all, in this study several different methods are developed and tried to build the slidblock and the particle clusters to confine the slidblock, including NO SHADOW command Method, Pluviation Method, Triangular Prism Method, Gravity Elimination Method、Particle Detection Method, etc. The above-mentioned methods are compared according to the availability and feasibility. The Particle Detection Method is composed of several procedures: Firstly, prepare one dense and cubic particle clusters and then transpose to the mass center of the slid block; Secondary, the particle cluster should be rescaled in order to cover the entire slid block; Third, detect the particlesif the position is situated within the slid blocks confined by the upper and lower surface of the DTM before and after landslide, otherwise, the particle is been eliminated from the cluster. After the procedure, all there maining particles are bonded according to the material properties and served as the model for further landslide simulation. Therefore an almost realistic geometry concerning the thickness and shapes of slid block could be established, and the sliding surface also could be built from DTM data to match the real terrain features of topography. This study considers earthquake ground excitation by the strong motion acceleration data gather from the CHY080 seismostation situated just nearby the sliding area. Furthermore the existence of pore water pressure within the rock mass and the effect of friction coefficient attenuation caused by rock mass shaking and sliding are also integrated in this study. Compare the topography results from the simulation data with the real DTM, it says that the gained accumulation patterns of the simulation result could be mostly similar to real situation (DTM after landslide) at the friction coefficient of 0.087. The average difference between the DTM and PFC results is -3.2m and the standard deviation is about 17.5. When landslides occur, the largest sliding distance is about 3662 m, the maximum drop height is about 632m and the highest sliding speed is 87.2 m/s. The timing of maximum sling velocity coincides with seismograph of CHY080. The results shows that use the numerical method of three-dimensional discrete element particle flow can simulate the complete landslide process dynamically. Comparing with the limit equilibrium method and finite element method which are the most common to be used, the discrete element particle flow provides not only the states of instability but also the dynamical process of sliding. The result from this study states the effect of this simulation and states the inference of the related earthquake-sliding phenomena as well.

Topic Category 工程學院 > 土木與防災研究所
工程學 > 土木與建築工程
工程學 > 市政與環境工程
Reference
  1. 【1】 中興工程顧問股份有限公司,「草嶺崩坍之處理對策評估」,經濟部水利處,1990。
    連結:
  2. 【4】 李宏輝,「砂岩力學行為之微觀機制-以個別元素法探討」,博士論文,國立台灣大學土木工程學系,台北,2008。
    連結:
  3. 【5】 李後彥,「數值地形模型於山崩及後續地形演化之應用-以草嶺山崩為例」,碩士論文,國立台北科技大學土木與防災研究所,台北,2009。
    連結:
  4. 【10】 林立欣,「離散元素模擬之平行策略」,碩士論文,國立台灣大學土木工程學研究所,台北,2005。
    連結:
  5. 【11】 林昇鴻,「卵形顆粒介質運動之數值模擬」,碩士論文,國立中央大學土木工程研究所,桃園,2006。
    連結:
  6. 【12】 林郁修,「分離元素法於岩石貫切破壞試驗之模擬分析」,碩士論文,國立台北科技大學土木與防災研究所,台北,2007。
    連結:
  7. 【14】 唐昭榮,「臺灣遽變式山崩傳送與堆積之顆粒流離散元素模擬」,博士論文,國立臺灣大學理學院地質科學研究所,台北,2010。
    連結:
  8. 【15】 高贈智,「集集地震引致台灣中部山區邊坡臨界滑移量之分析」,碩士論文,國立台灣大學土木工程學研究所,台北,2004。
    連結:
  9. 【16】 張來福,「車籠埔斷層周圍岩石力學特性之初探」,碩士論文,國立台灣大學土木工程學研究所,台北,2008。
    連結:
  10. 【17】 張家偉,「應用分離元素法探討紅菜坪地區地滑演化」,碩士論文,國立台灣大學土木工程學系,台北,2007。
    連結:
  11. 【18】 張家詮,「分離元素法於擬脆性岩材微觀破裂機制之初探」,碩士論文,國立台北科技大學土木與防災研究所,台北,2007。
    連結:
  12. 【19】 連俊凱,「具節理強度岩體因落石作用撞擊碎裂與堆積行為」,碩士論文,國立台灣大學土木工程學系,台北,2010。
    連結:
  13. 【20】 陳正旺,「車籠埔斷層周圍岩石力學特性之初探」,碩士論文,國立台灣大學土木工程學系,台北,2005。
    連結:
  14. 【21】 陳春華,「非連續體變形分析法模擬集集地震引致草嶺邊坡崩塌模擬之研究」,碩士論文,國立成功大學土木工程研究所,台南,2008。
    連結:
  15. 【22】 彭健豪,「九份二山地滑區滑動歷程與積行為之研究」,碩士論文,國立台灣大學土木工程學系,台北,2008。
    連結:
  16. 【26】 蕭安佑,「洪氾及山崩模擬-以金山及小林村為例」,碩士論文,國立台北科技大學土木與防災研究所,台北,2010。
    連結:
  17. 【27】 簡士堯,「低凝聚力材料趾部受侵蝕崩塌滑動模式研究」,碩士論文,國立台灣大學土木工程學系,台北,2009。
    連結:
  18. 【28】 Chen, R.-F., K.-J. Chang, J. Angelier, Y.-C. Chan, B. Deffontaines, C.-T. Lee, and M.-L. Lin (2006), Topographical changes revealed by high-resolution airbone LiDAR data: The 1999 Tsaoling landslide induced by the Chi-Chi earthquake, Eng. Geol., 88, 160-172.
    連結:
  19. 【29】 Chen, R.-F., Y.-C. Chan, J. Angelier, J.-C. Hu, C. Huang, K.-J. Chang, and T.-Y. Shih (2005), Large earthquake-triggered landslides and mountain belt erosion: The Tsaoling case, Taiwan, Comptes Rendus Géoscience, Acad. des Sci. Paris, 337, 1164-1172.
    連結:
  20. 【32】 Hung, J.-J. (2000), Chi-Chi earthquake induced landslides in Taiwan, Earthquake Eng. Eng. Seismol., 2 (2), 25-33.
    連結:
  21. 【34】 Kuo, C.-Y., Y.-C. Tai, F. Bouchut, A. Mangeney, M. Pelanti, R.-F. Chen, and K.-J. Chang (2009), Simulation of Tsaoling landslide, Taiwan, based on Saint Venant equations over general topography, Eng. Geol., 104, 181-189.
    連結:
  22. 【35】 Potyondy, D.O., and P.A. Cundall (2004), A bonded-particle model for rock, Int. J. Rock Mech. Min. Sci., 41, 1239-1364.
    連結:
  23. 【2】 何維信,航空攝影測量學,大中國圖書公司,台北,1995。
  24. 【3】 呂榮華,「以單顆粒勁度探討卵礫石層之抗剪強度」,碩士論文,國立中興大學土木工程研究所,台中,1999。
  25. 【6】 李憲忠,「九二一集集地震三維震源過程與震波傳遞分析」,博士論文,國立中央大學地球物理研究所,桃園,2003。
  26. 【7】 李錫堤、洪如江、林銘郎、蔡龍珆,「草嶺崩塌地工程地質調查與穩定性評估」,中興工程顧問社,專題研究報告,1993。
  27. 【8】 李錫堤、林銘郎、吳禮浩、鄭俊昇等人,「草嶺大崩山區的地質調查及歷次大崩山滑動面的決定」,1994 年岩盤工程研討會,1994。
  28. 【9】 李錫堤、洪如江、林銘郎、蔡龍珆等人,「草嶺崩塌地工程地質調查與穏定性評估-專題研究報告」,經濟部水資局統一規劃委員會,1993。
  29. 【13】 洪如江、李錫堤、林美聆、林銘郎、鄭富書、陳正與,「天塹可以飛渡、崩山足以斷流(草嶺順向坡滑動) 」,地工技術第77 期,2000。
  30. 【23】 黃鑑水、何信昌、劉桓吉等人,「台灣中部草嶺地區之地質與山崩」,經濟部地質調所彙刊,第二號,1983。
  31. 【24】 葉怡成,測量學-21世紀觀點,東華書局,2006,台北。
  32. 【25】 廖豐益,「以數值分析探討粒狀材料之破壞行為」,碩士論文,國立台灣工業技術學院營建工程技術研究所,台北,1994。
  33. 【30】 Cundall, P. A. and O. D. L. Strack, “A Discrete Numerical Model for
  34. 【31】 Cundall, P.A. (1971), A computer model for simulating progressive large scale movements in blocky rock systems. Proc. of the Sym, Int. Soc. Rock Mech., 1 No. II-8.
  35. Granular Assemblies,” Geotechnique, Vol. 29, No. 1, pp. 47-65 (1979).
  36. 【33】 Itasca (2006), User’s Manual, PFC 3D, Itasca Consulting Group, Inc., Minneapolis, Minnesota.
  37. 【36】 Togo, T., T. Shimamoto, S.L. Ma, and T. Hirose (2009), High-velocity friction of faults: A review and implication for landslide studies. The Next Generation of Research on Earthquske-induces landslides, An Int. Conf. Comm. 10th Anni. of the Chi-Chi Earthquake, 205-216.
  38. 【37】 USGS (2004), Landslide Types and Processes, Complied by Highland, L., USGS Fact Sheet 2004-3072.
Times Cited
  1. 周書玄(2013)。離散元素法初探材料微、巨觀參數之關係。臺北科技大學土木與防災研究所學位論文。2013。1-96。 
  2. 張賀翔(2012)。應用分離元素法探討順向坡失穩歷程-以國道三號崩塌事件為例。臺北科技大學土木與防災研究所學位論文。2012。1-99。 
  3. 林承翰(2014)。應用斜坡單元及分離元素法探討大型崩塌之演育。臺灣大學土木工程學研究所學位論文。2014。1-112。 
  4. 鄭憲聰(2014)。應用分離元素法於台二線邊坡災害模擬。臺北科技大學土木與防災研究所學位論文。2014。1-48。