崩塌地河道土砂堆積形貌變化之數值模擬

Numerical Modelling of Geomorphology of Landslide Induced Sediment Deposition in Incisional Channel

高嘉廷(CHIA-TING KAO)

指導教授 ： 游景雲

國立臺灣大學/工學院/土木工程學研究所/碩士(2019年)

https://doi.org/10.6342/NTU201903259

摘要

崩塌為集水區帶來大量泥砂，也是河川重要的土砂的來源。過去雖然有學者研究其河道發展，然而，如何量化輸送到下游的崩塌產生土砂量仍然是一個挑戰。本研究致力於量化受水流侵蝕的崩塌體的形貌變化和河川演化機制。我們使用去藕合二維模型來研究它的崩塌地貌受河川侵蝕演化過程。
本研究在水理方面，選用淺水波方程式來模擬水流流況。同時,利用流場求出河川輸砂率,在輸砂方面,則使用Exner方程式來計算泥沙傳輸質量的守恆，進而看出河川形貌的演化。此外，兩個方程都是以水平二維基礎以可侵蝕的河岸及沒有二次流假設下，選取的非線性項的控制方程，然而，解析方法無法解決此問題。因此，本研究採用有限體積法進行數值模擬。初步由流場結果發現，侵蝕最嚴重的區域為崩塌的下緣處，最為可能會發生二次崩塌的地方，而在進一步的河川演化過程，可以發現在本研究之假設條件下，結果發現於攻擊岸下游處為淤積處，主要控制原因為流場分布，原有的崩積體也由半圓形變得接近半橢圓,總體來說,河川的發展愈來愈趨近趨近於曲流。

關鍵字

淺水波方程 ； Exner方程式 ； 有限體積法 ； 崩塌堆積 ； LHLL ； 去耦合模型

並列摘要

Landslides produce enormous volumes of sediment in mountainous watersheds. However, quantifying landslide-derived sediment transported downstream remains a challenge. This study makes an effort in determining the amount of landslide-derived sediment that eroded by the water and the ratio of sediment yield of this problem. We try to use a decoupled two-dimensional model to investigate the evolution of landslide aggregation influenced by the channel flow.
In terms of the hydrodynamics model, we choose the shallow water equations to simulate the flow condition. Meanwhile, the flow field is used to determine the sediment transport rate of the river. In the aspect of sediment transport, the Exner equation is used to estimate the conservation of sediment transport and investigate channel evolution. Besides, both of the two equations under a horizontal 2D straight, non-erodible banks, and no secondary flow channel can be derived as our governing equation which has the nonlinear terms, and the analytical method may not work to solve the two equations under the 2D condition. As a result, this study adopts the finite volume method to proceed with this work. According to the preliminary results of the flow field, it is found that the most severely eroded area is at the lower edge of the landslides aggregation, and the place where the second collapse is most likely to occur. In the process of further river evolution, it can be found that under our assumptions, the results suggest that the deposition is in the downstream part of the attack shore and channel evolution is dominated by the flow field distribution. Also, the original landslides aggregation transforms a semi-elliptical into a semi-circle. In summary, the river develops the meandering stream gradually.

並列關鍵字

Shallow water equations ； Exner equation ； Finite volume method ； Landslide aggregation ； decouple model

參考文獻

1.Alcrudo, F., & Garcia-Navarro, P. (1993). A high-resolution Godunov-type scheme in finite volumes for the 2D shallow-water equations. International Journal for Numerical Methods in Fluids, 16(6), 489-505. doi:10.1002/fld.1650160604

2.Berthon, C., Cordier, S., Delestre, O., & Le, M. H. (2012). An analytical solution of the shallow water system coupled to the Exner equation. Comptes Rendus Mathematique, 350(3), 183-186. doi:https://doi.org/10.1016/j.crma.2012.01.007

3.Bradford, S. F., & Sanders, B. F. J. J. o. H. E. (2002). Finite-volume model for shallow-water flooding of arbitrary topography. 128(3), 289-298.

4.Brufau, P., Vázquez‐Cendón, M., & García‐Navarro, P. J. I. J. f. N. M. i. F. (2002). A numerical model for the flooding and drying of irregular domains. 39(3), 247-275.

5.Cantelli, A., Wong, M., Parker, G., & Paola, C. (2007). Numerical model linking bed and bank evolution of incisional channel created by dam removal. Water Resources Research, 43(7).