水利工程論文-泥石流泥沙體積濃度之研究.doc

水利工程論文-泥石流泥沙體積濃度之研究摘要：本研究系探討泥石流先端部和整體平均(包括先端部及后續(xù)水流)泥沙體積濃度隨著溪床坡度之變化規(guī)律，以及其縱剖面泥沙體積濃度沿著水深的分布特性。經(jīng)采用其他研究者渠槽試驗資料進行檢定及驗證后發(fā)現(xiàn)，當系數(shù)=0.1時，本研究所提出泥石流平衡泥沙體積濃度之半經(jīng)驗公式，不僅可模擬泥石流先端部泥沙體積濃度隨著溪床坡度之變化趨勢，亦可適于推估不同溪床坡度下之泥石流整體平均泥沙體積濃度。同時，本研究也應(yīng)用最大熵原理及平衡泥沙體積濃度半經(jīng)驗公式，導出了縱剖面泥沙體積濃度分布之理論公式，經(jīng)Tsubaki實測資料驗證結(jié)果顯示，本文模式可以掌握縱剖面泥沙體積濃度之分布趨勢而優(yōu)于Tsubaki及Chen兩研究者之理論公式。關(guān)鍵詞：泥石流體積濃度模式分析1前言泥石流是指發(fā)生于陡峻的溪谷或斜坡面上未固結(jié)之沉積物或經(jīng)風化之礫石、巖屑等堆積土體，受豪雨形成的地表逕流或土層地下水位上升等水文過程之影響，而失去原有安定之狀態(tài)，土砂礫石伴隨著洪流在重力作用下沿著自然坡面形成的一種高濃度集體流動現(xiàn)象，對區(qū)域環(huán)境安全和生產(chǎn)的影響甚巨，而逐漸引起人們的重視。在過去的二十余年間，已有不少學者專家投入泥石流的形成、運動、停積及治理措施等各方面課題之研究，并獲得若干成果。但受限于對泥石流形成過程及其流動特性變化規(guī)律的了解程度，致使在探討其行為機制的各項相關(guān)因子時，仍不免存在著一些盲點有待解決，例如泥石流泥沙體積濃度(或泥沙含量)的變化規(guī)律，就是一個相當顯著的例子。泥石流泥沙體積濃度系指單位體積中固體泥沙顆粒所占有的體積，屬于泥石流基本流體特性之一，其變化規(guī)律密切影響著泥石流的行為機制，而與泥石流運動、流出土方量、泛濫范圍及其成災規(guī)模等問題關(guān)系密切，為泥石流防災相當重要的指標參數(shù)之一，惟目前除了日本高橋(Takahashi)的理論模式外，相關(guān)研究成果卻是寥寥可數(shù)。根據(jù)高橋建構(gòu)泥石流平衡泥沙體積濃度理論公式所引用的理論基礎(chǔ)，可以有兩種不同的表達形式，一是以巴格諾(Bagnold)高泥沙含量固液兩相流的本構(gòu)關(guān)系為基礎(chǔ)；另一則是運用了水力學原理。茲分述如下：巴格諾1從實驗中發(fā)現(xiàn)，當固液兩相流中泥沙含量很高時，由于固體顆粒間可以透過不斷地碰撞來進行動量轉(zhuǎn)換而形成兩種不同的力量，其中沿著水流方向的動量交換會形成粒間離散剪應(yīng)力(,inter-particledispersiveshearstress)，而與水流垂直方向的動量交換則形成粒間離散應(yīng)力(Pd,interparticledispersivestress)。同時，他也提出了粒間離散剪應(yīng)力與粒間離散應(yīng)力間的本構(gòu)關(guān)系(constitutiverelationship)可表為=-Pdtan(1)式中，tan=動摩擦系數(shù)，當水流處于完全慣性區(qū)(fullyinertialrange)時，tan0.321。高橋2即利用此一關(guān)系式，經(jīng)由泥石流固液兩相的動量守恒方程導出了泥石流平衡泥沙體積濃度公式Cd=tan/(s-)(tan-tan)(2)式中，Cd=平衡泥沙體積濃度；=溪床傾斜角度；s、=固體泥沙顆粒及流體之密度。由上式得知，因不能使用于tantan0.32的溪床條件，加上泥石流常發(fā)生在溪床傾角15的陡峻地形上，使得應(yīng)用式(2)在計算泥沙體積濃度時，因合理值范圍相當有限而難以符合實況。例如，取tan=0.32、tan=0.268(=15)、s=2.6g/cm3及=1.0g/cm3代入式(2)，則平衡泥沙體積濃度Cd=3.2,該值已超出泥沙體積濃度之合理值，實際上是不存在的。因此，高橋運用水力學原理推導平衡泥沙體積濃度公式3Cd=tan/(s-)(tan-tan)(3)式中，tan=靜摩擦系數(shù)。與式(2)比較得知，上式分母中系采用靜摩擦系數(shù)tan，而非動摩擦系數(shù)，是兩方程式間的主要不同點。惟高橋經(jīng)實驗驗證顯示，式(3)的形式在模擬泥石流泥沙體積濃度變化規(guī)律上優(yōu)于式(2)，故乃建議采用式(3)作為預測穩(wěn)定均勻條件下的泥石流先端部泥沙體積濃度。不過，Chen4曾對此提出看法認為，泥石流既已達穩(wěn)定均勻流況，其流動型態(tài)自應(yīng)處于完全慣性狀態(tài)，固體泥沙顆粒間的摩擦作用理當采用動摩擦系數(shù)作為參數(shù)，較為合理，而非靜摩擦系數(shù)。筆者5則自非牛頓流體流變方程式出發(fā)，認為泥石流在高濃度、高剪應(yīng)變率時，由于顆粒間始終保持直接接觸，不同粒徑的石礫相互疊置碰撞，使顆粒間除承受來自于顆粒間彈性碰撞所產(chǎn)生的離散剪應(yīng)力作用外，還必須考慮顆粒間之摩擦阻力。這樣，泥石流先端部泥沙體積濃度可表為Cd=(3-2)tan/(s-)(2-3)tan+3tan(4)式中=摩擦阻力系數(shù)(resistancecoefficient),且01。根據(jù)上式，若不考慮粒間摩擦阻力系數(shù)，取=0時，可獲得和式(3)度表達式。歐國強6于流槽通過實驗方式研究泥石流先端部及其后續(xù)水流兩部分之整體平均泥沙體積濃度，并建立了以溪流坡度為主的泥沙體積濃度經(jīng)驗公式，即Cdt=4.3Cm(tan)1.5/1+4.3Cm(tan)1.5(5)式中Cdt=泥石流整體(包括先端部及后續(xù)流)平均泥沙體積濃度；Cm=溪床面靜止泥沙之最大泥沙體積濃度(thegrainconcentrationinvolumeinthestaticdebrisbed)。綜合以上各泥沙體積濃度公式繪出泥沙體積濃度與溪床坡度之關(guān)系曲線，如圖1所示。圖中顯示圖1泥石流泥沙體積濃度與溪床坡度之關(guān)系曲線Relationofsedimentconcentration,Cd,andchannel-bedslope,tan,inthedebrisflow1.高橋與筆者之泥石流體積濃度曲線均呈向上揚升趨勢，當溪床坡度略大于某數(shù)值時，其泥沙體積濃度即因快速增加而使之大于1.0，這與泥沙體積濃度的基本定義相互矛盾。2.歐國強泥沙體積濃度曲線系隨著溪床坡度增加而漸趨平衡，即使在較大的溪床坡度下，泥沙體積濃度亦未超出合理的范圍，惟因?qū)僖唤?jīng)驗公式，其通用性仍有待檢驗。3.不論是泥石流先端部泥沙體積濃度或其整體平均泥沙體積濃度，從單一泥石流段波的整體流動特性來看，它們之間理應(yīng)具有某種程度的共通特性，并可依相同或近似的機制建立其理論模式。惟目前系分別采用兩種截然不同的機制及公式予以計算和模擬，使得先端部和整體平均泥沙體積濃度間的變化過程無法連貫，其機制自是無法統(tǒng)一地表達。因此，虞及泥沙體積濃度對泥石流數(shù)值模擬和實際運用上的重要性，本文擬通過一定的理論過程和渠槽試驗資料，建立泥石流泥沙體積濃度之通用公式，不僅可以修正高橋公式在模擬泥石流先端部泥沙體積濃度上不盡合理的變化趨勢，亦可運用于推估泥石流的整體平均泥沙體積濃度。2理論分析考慮一無限長且厚度均勻的無粘聚性飽和堆積土層，厚度為D，溪床面傾角為。當土層表面形成水深h0及泥沙體積濃度C0的水砂混合水流時，受重力作用之影響，水流不斷地沖刷床面泥沙顆粒，使泥沙含量沿程增加，從而提高了流和力量。由于這種水流的行為機制系受沿程泥沙含量的增減所左右，不僅具有碰撞過程的應(yīng)力轉(zhuǎn)換效應(yīng)，亦與溪床表面泥沙顆粒的沖刷互為相關(guān)。因此，同時考慮巴格諾流體內(nèi)部固體泥沙顆粒間的本構(gòu)關(guān)系和水力學溪床表面的泥沙沖刷機制，將水流中泥沙顆粒間的離散剪應(yīng)力和剪阻力轉(zhuǎn)化為凈剪應(yīng)力(netshearstress)，并與溪床表面泥沙顆粒的起動應(yīng)力相結(jié)合，以決定水流對溪床表面泥沙的沖淤作用。這樣，水流沿著流動方向上的凈剪應(yīng)力，可表為粒間離散剪應(yīng)力與剪阻力之差，即r=-Pdtan=(s-)C0+gh0sin-(s-)gC0h0costan=gh0sin1+s-/(6)C0(1-tan/tan)圖2無因次凈剪應(yīng)力與泥沙體積濃度之關(guān)系Relationofdimensionlessnetshearstressandsedimentconcentration圖3表面流作用床面堆積土層模型圖Thesimplifieddiagramofflowdepositlayersonthebedsurface由上式得知，水流凈剪應(yīng)力系隨著泥沙體積濃度C0的變化而改變，且與溪床坡度和動摩擦系數(shù)相關(guān)。為此，將上式寫為無因次形式，即T*=r/ghosin=1+s-/(1-tan/tan)(7)如圖2為在給定的動摩擦系數(shù)下，不同溪床坡度時無因次凈剪應(yīng)力與泥沙體積濃度之關(guān)系圖。圖中，無因次凈剪應(yīng)力T*隨著泥沙體積濃度C0之增加，分別呈現(xiàn)上升或下降，其趨勢與溪床坡度和動摩擦系數(shù)間的大小互為相關(guān)；當tantan時，無因次凈剪應(yīng)力T*呈緩和的上升趨勢，而在tantan的情形下，無因次凈剪應(yīng)力T*呈快速的下降趨勢。但是，不論水流凈剪應(yīng)力隨著泥沙體積濃度呈上升或下降，來自于溪床表面泥沙顆粒之阻滯效應(yīng)，水流凈剪應(yīng)力和床面泥沙顆粒間的互動機轉(zhuǎn)中，總會朝著某種動態(tài)的平衡關(guān)系發(fā)展，而使床面泥沙顆粒不再被大量的沖刷或者淤積；換句話說，在適當?shù)拇裁鏃l件配合下，水流泥沙含量終將趨于穩(wěn)定，不再作激烈的增減，此時水流凈剪應(yīng)力和床面泥沙顆粒間會存在一種等勢的關(guān)系，表為方程式，則可寫為r=fcr(8)當水流和床面泥沙顆粒具有以上的關(guān)系時，水流與溪床間的泥沙交換將達到平衡，水流泥沙體積濃度趨于飽和，溪床泥沙既不沖刷也不淤積。式中fcr為溪床面泥沙顆粒起動之臨界剪應(yīng)力(criticalshearstressforincipientmovementofsedimentparticlesonthechannel-bed)，可表為單位面積溪床面上泥沙顆粒水中自重沿著流動方向上的抗沖力，即fcr=(s-)gdscos(9)式中ds=床面泥沙特征粒徑；=綜合因子，它是靜摩擦系數(shù)、泥沙顆粒形狀、級配、排列及其粘聚特性等因素綜合影響的修正系數(shù)，需由實測資料決定之。根據(jù)式(8)得知，當水流凈剪應(yīng)力與床面泥沙顆粒起動之臨界剪應(yīng)力相等(r=fcr)時，因床沙與水流間的泥沙交換已趨平衡，水流中必然存在一飽和的泥沙體積濃度Cd，使其得以穩(wěn)定均勻流條件向下游行進，故由式(6)和式(9)分別代入(8)式中，經(jīng)整理可得，恒定均勻水流條件下之平衡泥沙體積濃度Cd=tan/(s-)(tan-tan)-/tan-tands/h0(10)上式適用于一般挾沙水流及泥石流之流況。式中，泥沙粒徑與水深比值(ds/h0)為與水流之泥沙體積濃度相關(guān)，其關(guān)系式可透過一物理圖形予以求解，如圖3所示，為表面流作用床面堆積層之簡單模型圖。假設(shè)表面流從斷面0-0至斷面1-1的有限距離x中，可移動床面泥沙厚度平均達h，則由砂礫連續(xù)方程式可得C0U0h0+CmUx=CdU1h1(11)式中床面泥沙沖刷速度U=h/t，且x=(U0+U1)t/2。為簡化問題起見，假設(shè)入流速度U0與出流速度U1相接近時，即U0U1，又h1=h0+h，則上式經(jīng)整理可得h/h0=Cd-C0/Cm-Cd(12)設(shè)表面流于斷面0-0的泥沙體積濃度C00；同時，因床面泥沙移動厚度h至少必須大于或等于單一泥沙顆粒之特征粒徑ds,即h=nds,且n1，故上式可改寫為ds/h0=1/nCd/Cm-Cd(13)或者，選取圖3斷面1-1處單位底面積之柱體分析，當斷面0-0泥沙體積濃度等于零時，則柱體全部泥沙量h1Cd應(yīng)來自于床面堆積泥沙，故有h1Cd=(h+h0)Cd=hCm(14)上式經(jīng)化簡后與式(13)完全一致，這顯示式(13)在推導過程的部分假設(shè)尚稱合理。因此，將式(13)代入式(10)中，經(jīng)整理可得恒定均勻水流平衡泥