甘肅省南部泥石流運動規律的分析與研究
——以甘肅舟曲三眼峪溝為例*

舒和平，韓　拓，齊　識，馬金珠

(蘭州大學 資源環境學院西部環境教育部重點實驗室，甘肅 蘭州 730000)

?

甘肅省南部泥石流運動規律的分析與研究
——以甘肅舟曲三眼峪溝為例*

舒和平，韓拓，齊識，馬金珠

(蘭州大學 資源環境學院西部環境教育部重點實驗室，甘肅 蘭州 730000)

摘要:以相似性原理為基礎，選擇恰當的比尺，在甘肅省隴南市武都區野外構建泥石流物理模型，重現了舟曲“8·8”特大泥石流災害發生過程。通過物理實驗模擬，測得泥石流運動的平均流速為0.741 7 m/s，沖擊力為1.99 N，最大堆積厚度為6.1 cm，堆積面積達3.64 m2,并勾勒出堆積形態。選取經驗公式，將實測值與經驗公式計算結果相比較，驗證了物理模型的合理性，進一步分析了泥石流的運動規律。根據泥石流形成的要素提取危險度評價因子,結果表明三眼峪為高度危險泥石流溝，這為泥石流災害防治工程提供可靠的科學依據。

關鍵詞:三眼峪；泥石流；模型；危險度；沖擊力；流速；甘肅舟曲

泥石流是一種廣泛分布于高寒、溫暖濕潤山區及干旱半干旱山區的常見自然現象，常發生于山區溝谷和坡地，是山區溝谷和坡地地貌的外在營力，是山地環境退化、生態平衡失調的產物，具有流量大、容重大、流動速度快、歷時短暫和破壞性大的特點[1-4]。我國是世界上遭受泥石流災害最嚴重的國家之一，尤其是甘肅省南部位于新生代亞洲印度板塊碰撞邊界，新構造運動頻繁；加之該區生態脆弱、降雨頻繁集中，且多暴雨，因此該區成為我國地質災害高發區之一。三眼峪溝位于甘肅省南部舟曲縣，自1823年以來,該溝就曾爆發過約11次大型泥石流, 以2010年的“8.8”舟曲特大泥石流災害最為嚴重[5]。因此，針對具有代表性的舟曲泥石流，研究其發生機理、運動及堆積規律，具有很強的現實意義,可為泥石流災害防治工程提供可靠的科學依據。

為了研究泥石流的運動規律，國內外學者做過不少的模擬實驗[6-9]。日本泥沙所高橋保教授做了水槽實驗,并首次在實驗中測量了縱向速度分布[10]；王光謙等[11]在云南東川蔣家溝泥石流觀測研究站現場進行了粘性泥石流現場模擬實驗，構建出了泥石流堆積模型，并進行了不同堆積坡度的模擬，驗證出了泥石流團模型的準確性；劉麗和陳洪凱[12]首次提出了“通過實測泥石流運動速速，利用頻率分析，建立泥石流流速譜、沖擊力譜等重要參數”，進而準確分析泥石流運動規律；陳中學等[13]在室內通過顆粒級配篩選出9種不同量級的土體，并設計出邊坡模型，模擬了在雨強為85 mm/s下的降雨實驗。但這些模型適用范圍較窄，對泥石流的運動過程研究不夠全面[14-16]。因此，選取具有代表性的舟曲泥石流溝，建立物理模型，進行運動過程模擬，實測流速、沖擊力及沖淤范圍等運動因素，另外將理論經驗公式計算結果與實測結果進行對比分析，并驗證模型的合理性，分析其運動規律，為當地相關政府部門在防治減災過程中提供建議，從而盡量減少泥石流災害對人民生產生活的危害和影響。

1研究區概況

舟曲縣位于甘肅省南部甘南藏族自治州，地理坐標是105°51′39″～104°45′31″E，33°13′06″～34°01′00″N。它地處青藏高原東緣，西秦嶺西翼與岷山山脈交匯地區，屬構造、侵蝕山地, 舟曲縣河流屬長江流域嘉陵江水系，主要發育河流有白龍江、拱壩河等。舟曲縣屬北亞熱帶向北溫帶的過渡區，受大氣環流和地形影響，具有垂直氣候分帶明顯和干濕季分明兩大特點。區內多年平均降雨量為435.8 mm，日最大降雨量為63.3 mm，小時最大降雨量為47.0 mm，降雨以連陰雨和暴雨形式出現，降水主要集中在6-9月。氣溫變化較小，晝夜溫差不大，多年平均氣溫12.9 ℃，最熱7月平均氣溫23.0 ℃[17]。

舟曲縣城依白龍江而建，夾于南北兩山之間，依山傍河，地域狹窄，主城區建于三眼峪溝、羅家峪溝、寨子溝、硝水溝等泥石流堆積扇上。三眼峪溝流域位于舟曲縣城北側，屬白龍江一級支流，流域面積24.1 km2(圖1)。區域上該處為西秦嶺南帶印支期冒地槽褶皺帶西段白龍江復式背斜北翼，整個流域處在兩條北西走向的區域性斷裂帶所夾的斷塊之間，與區域性斷裂走向一致的兩條正斷層、一條逆斷層分別橫穿溝谷的中段和前段[18]。斷裂兩側巖層破碎，褶曲強烈。流域內出露的前第四紀地層以二疊系為主，分布于流域中、上游，巖性主要為厚層塊狀灰巖、中到薄層硅質條帶灰巖、白云質灰巖、大理巖化灰巖及白云巖化鮞狀灰巖，巖層產狀192°～210°∠68°～ 72°。流域下游，即出山口一帶，出露中泥盆統古道嶺組，巖性為炭質板巖、千枚巖夾薄層灰巖及砂巖，巖層產狀234°∠54°。流域內第四系廣泛分布，其成因類型復雜，主要有風積黃土、殘坡積、沖洪積及崩塌、滑坡堆積物等[19]。

圖1　研究區位置圖

2實驗方法

通過調查泥石流災情和災害特征，獲取泥石流運動規律的基本要素，總結以往發生泥石流災害發生事件和技術方法，統計分析后，建立了泥石流運動規律研究技術路線，最后計算出泥石流災害發生危險程度，具體技術路線如圖2所示。

圖2　技術路線圖

(1)野外構建物理模型，重現舟曲“8·8”特大泥石流災害發生過程；

(2)通過物理實驗模擬，得到流速、流量、沖擊力、堆積形態、堆積面積、堆積厚度等多組數據，并進行統計分析；

(3)選取經驗公式，計算流速、流量、沖擊力；

(4)通過對比實測實驗數據與經驗計算數據，驗證物理模型的合理性。

(5)針對物理模型模擬時泥石流運動形成的堆積扇，分析其堆積面積、堆積形態、堆積厚度；

(6)根據泥石流形成的要素提取危險度評價因子,劃分危險度等級。

3實驗設計

3.1模型構建

模型建立是實驗的關鍵，既要保證微地形和地物同原型相似，又要確保模型實驗能準確反映原型的演變情況，因此在嚴格相似有困難的地方，就要合理進行簡化和理想化。按原型區實測地形圖和野外調查資料，并參照錄像和照片，進行模型制作。因此在甘肅省隴南市武都區典型泥石流溝段河壩溝建立了物理模擬實驗平臺，實驗平臺的平面示意圖及全景圖如圖3和圖4所示。

表1　評價因子表

表2　泥石流危險評級因子分級、評分表

(1)泥石流補給箱，為四棱錐體，邊長為120 cm，底部設置正方形供給孔，邊長為 20 cm，實驗時，泥石流的補給由閘板開關控制。

(2)泥石流流通槽，形態呈矩形，其內寬 40 cm，內高 40 cm，有效流動長度為500 cm。

(3)泥石流堆積板，為平面鐵板，其寬 200 cm，長 300 cm，以斜口方式將其與流通槽以平整形式銜接，模擬時可隨意調節其坡度。

通過記錄堆積扇在鋼質米格板上的分布，勾勒出堆積形態，以數格子的形式推算堆積面積，并測量出扇緣、扇中、扇尾堆積厚度。可進一步分析一次泥石流過程后的堆積形態及堆積量，為以后危險度的劃分提供了依據。

圖3　實驗平臺平面圖

圖4　實驗全景圖

3.2經驗公式

在泥石流運動規律分析、危險性評價和泥石流預防過程中，平均流速、流量以及沖擊力的研究是十分重要的[20-22]。國內外學者經過多年的數學建模，總結出了許多計算平均流速和沖擊力的經驗公式[23]。因此在研究三眼峪溝運動規律過程中，選取了通過顆粒級配、溝床比降、過水斷面面積等水力學參數的計算方法來計算[24-25]，該方法所需參量簡便易得，計算簡單。它不僅使我們對泥石流災害有更深入的認識，并且對于以后的防災減災也有重大的現實意義[26]。

泥石流斷面平均流速[25]：

(1)

泥石流流量：

Q=V×S。

(2)

泥石流沖擊力[27]：

F=Adρdv1vd。

(3)

式中：γ為水力半徑(m)；S為過水斷面面積(m2)；i為溝床比降(%)；g為重力加速度9.81(m/s2)；R50為泥沙顆粒中百分比小于50%的顆粒粒徑(mm)；R10為泥沙顆粒中百分比小于10%的顆粒粒徑(mm)；C為泥石流的容重(t/m3)；t為泥石流持續時間(s)；Ad為石塊與被撞物體接觸面積；ρd為石塊密度；v1為縱波波速；vd為石塊運動速度(視為與泥石流體速度相等)。

3.3危險度評價

危險度的劃分采用了簡單直觀的模型賦值法，這種方法有三種模型：分級評分模型、分級賦值模型和函數賦值模型。由于研究的區域為三眼峪，因此選用了符合單溝泥石流危險度評價模型中的分級評分模型，其具體劃步驟如下所示。

(1) 根據實地情況選取泥石流危險度評價因子。我們共選取了7個評價因子，分別是：流域面積p1，主溝比降p2，相對高差p3，主溝長p4，松散固體物質儲量p5，植被覆蓋率p6，24h最大降雨量p7,各參量數值如表1所示。

(2) 依據一定的原理方法和測度(如關聯度)，確定各危險度評價因子的權重和權數。

(3) 根據一定分級評分方法確定單溝泥石流危險度查找表(如表2所示), 然后計算各危險因子的分數。

(4) 計算各危險因子在危險度查找表中的評分總和, 即為單溝泥石流危險度, 其危險度W計算公式可表達為：

(4)

式中:Pi為第i個評價因子在危險度查找表中的評分；n為評價因子總個數。

(5)根據表3的危險度劃分區間確定泥石流溝危險等級。

表3　危險度劃分區間表

4實驗結果

4.1參數律定

實驗模型目的是為泥石流災害整治工程提供可靠的設計依據。首先將三眼峪溝道原型概化為實物模型，并在模型上重演與原型條件相似的流體運動。再將實驗模型成果應用到原型上，因而要求模型與原型之間應滿足一定的相似條件。此次實驗要求實驗流體與原型流體的水動力學特征值相似，原型與模型流體的容重相等且所含固體顆粒組成相似并遵守模擬邊界條件，模型流體與原型流體的弗勞德數相等。實驗采用正態動床模型，根據實驗目的要求和模型場地大小條件，模型和原型之間各物理量相似比尺如表4所示。

表4　各物理量相似比尺

4.2流速和流量計算

4.2.1經驗計算值

根據舟曲縣三眼峪溝相關地形參數(5)和水動力參數(6)，然后利用經驗公式進行流速計算。

將相關地形參數及水動力學參數分別帶入式(1)和式(2)，即得到斷面平均流速v和流量Q的理論計算值，即：v=8.025 m/s；Q= 410.48 m3/s。

4.2.2實驗值

把上述所得數據用相似性比尺處理，得到實際參量、模型參量及比尺對照表(表7)。

流通區流速及泥位測量采用LDM41型激光測距儀觀測。激光測距傳感器采用相位比較原理來測定出傳感器激光發射點與物源之間的距離。實驗所選用的激光測距儀為10 Hz，即每秒記錄10個距離h0，將儀器端口到流通槽的距離設為H，則得到泥位Δh=H-h0。將流通槽底部寬度設為a，則斷面平均流量q為a×h，進而單位時間斷面平均流量為a×h,即斷面瞬時流速v為40×a×h。其中，已知H=1.4 m，a=0.4m。具體流速和流量如表8所示。

表5　三眼峪相關地形參數

表6　三眼峪溝水動力學參數

表7　實際參量、模型參量及比尺對照表

表8　流速和流量表

4.3沖擊力

4.3.1經驗計算值

基于表2和表3的基礎數據，并參考相關參數：由遙感地質圖可知，研究區以灰巖為主，其密度約為2.74 g /cm3；通過查閱“巖性密度縱波波速測試表”得到，該區縱波波速為4 614 m/s；石塊與被撞物體的接觸面積按堆積扇上最大石塊的中、短徑所在平面面積的10%計算[28]；經文獻查閱及實地調查，該區泥石流體流速約為5.4 m/s。將相關參數代入式(3)中，得到沖擊力：F=19 784 t 。

4.3.2實驗值

實驗中，沖擊力的測量采用光柵傳感分析儀，SA-1系列光柵傳感分析儀是采用掃描激光技術的光纖光柵傳感解調分析儀器。在流通槽上、中、下部3個斷面分別布設圓餅狀壓力傳感器FX-FBG-P-0004、FX-FBG-P-0005、FX-FBG-P-0006，再將壓力傳感器與調節儀相連，最后連入電腦。當泥石流體流過傳感器時會產生一定的應力，力以數據的形式呈現在電腦上,具體實驗數據如表9所示。

圖5　應力變化趨勢圖

三個測點的應力總體呈現下降趨勢。其中測點1應力變化平緩；測點2應力先變小再變大，最后減小；測點3應力變化幅度也相對較小。經過計算可知：實驗沖擊力為1.99 N，換算成質量為0.204 kg，而實際沖擊力為19 784 t，進而算得相似性比尺約為1:96 000 000，符合相似性標準。

表9　實測沖擊力數據統計表

4.4堆積結果

鋼質米格板由150的小正方形組成，每小格邊長20 cm，通過堆積形態圖分析可知：堆積物所占格數約為91格，總格數為150格,占鋼板總面積的60.67%。堆積形態呈梨形。最高泥位出現在縱向第5格和第6格之間、橫向第2格交匯處，最大泥位深度為6.1 cm。具體堆積形態如圖6所示。

圖6　實測堆積形態分布圖

將圖6中最上面第一條橫線設為橫0線，第二條橫線設為橫1線，以此類推，最下面一條橫向線則為橫15線；同理，將左數第一條縱向線設為縱0線，則右端最后一條縱向線為縱10線。其交點以坐標形式記錄，則部分交點泥位深度表如表10所示。

表10　堆積區堆積深度統計表　cm

根據表10可以得出如下幾個結論。

(1) 縱5線為泥石流堆積體的軸線，最大堆積厚度就出現在軸線尾部，且在軸線后部堆積厚度明顯高于左右兩側。

(2) 在扇尾，堆積厚度明顯高于扇中及扇緣。

(3) 在扇中，左側堆積厚度略高于右側；在扇緣，右側堆積厚度則高于左側。

(4) 軸線在扇緣堆積厚度較淺，且低于左右兩側。

4.5危險度

5結論

通過對舟曲三眼峪溝泥石流運動相關參數進行分析，依據相似性原理，選擇恰當的比尺，在隴南市武都區野外構建了泥石流模型，然后進行模擬實驗，測出流速、沖擊力，并將模型實測數據與經驗公式計算值對比分析，通過堆積面積、堆積厚度、堆積形態劃分危險度。得出的主要結論如下。

(1) 舟曲泥石流形成區面積達22.64km2，占流域總面積的94%。形成方式有支溝雨洪匯流型和主溝沖蝕型兩類，為“8.8”舟曲特大泥石流災害提供了豐富的物源；流通區位于大小眼峪交匯處至峪門口溝段,占流域總面積的2.4%，徑流通道對泥石流的形成及流速、流量、沖淤情況等起著控制作用；三眼峪口至白龍江之間為泥石流的堆積區，占流域總面積的1.7%。

(2) 通過公式計算出流速為8.025m/s、流量410.48m3/s、沖擊力為19 784t。經過相似比尺計算，在誤差允許范圍內與實測值吻合，驗證了模型的合理性。

(3) 構建物理模型，概化后完全遵循相似性原理，對泥石流模型進行調節，旨在保證模型實驗的合理性和科學性，從而再現“8·8”特大災害性泥石流發生時的情形；經過實驗測試得到平均流速為0.741 7m/s，沖擊力為1.99N，最大堆積厚度為6.1cm。

(4) 對泥石流體以及堆積扇的堆積面積、堆積形態、堆積厚度進行了詳細的分析研究，根據泥石流形成要素提取危險度評價因子,計算得到危險度值為71，確定三眼峪溝為重度危險泥石流溝。

參考文獻：

[1]王繼康.泥石流防治工程技術[M].北京:中國鐵道出版社,1996.

[2]陳洪凱,唐紅梅,馬永泰,等.公路泥石流研究及治理[M].北京:人民交通出版社,2004.

[3]任非凡.G212線隴南段泥石流危險性區劃研究[D].蘭州:蘭州大學,2006.

[4]康志成, 李焯芬,馬藹乃,等.中國泥石流研究[M].北京:科學出版社,2004.

[5]胡桂勝,陳寧生,鄧明楓,等.甘肅舟曲三眼峪溝泥石流粗大顆粒沖擊力特征分析[J].地球與環境,2011,39(4):478-484.

[6]唐興邦,章書成.泥石流研究[J].中國科學院院刊,1992,31(2):119-123.

[7]HutterK,SvendsenB,RickenmannD.Debrisflowmodeling:Areview[J].ContinuumMechanics&Thermodynamics,1994,8(1):1-35.

[8]GordonEG.Criticalflowconstrainflowhydraulicsinmobile-bedstreams:anewhypothesis[J].WaterResourcesResearch,1997,33(2):349-358.

[9]MergiliM,FellinW,MoreirasSM,etal.SimulationofdebrisflowsintheCentralAndesbasedonOpenSourceGIS:possibilities,limitationsandparametersensitivity[J].NaturalHazards,2011, 61(3):1051-1081.

[10]Iverson R M, Vallance J W. New views of granular mass flows[J]. Geology, 2001, 29(29):115-118.

[11]王光謙,邵頌東,費祥俊. 泥石流模擬:I-模型[J].泥沙研究,1998(3):7-13.

[12]劉麗,陳洪凱.泥石流運動力學研究現狀及趨勢[J].重慶交通大學學報(自然科學版),2010,29(2):233-239.

[13]陳中學,汪稔,胡明鑒,等.黏土顆粒含量對蔣家溝泥石流啟動影響分析[J].巖土力 學,2010,31(7):2197-2011.

[14]魏新功,王振國,包紅霞.降水原因造成的舟曲縣地質災害分析[J].甘肅科技,2008,24(21):84-88.

[15]趙成.“8.8”舟曲暴雨泥石流的成災模式[J].西北地質,2011,44(3):63-70.

[16]韓國營.甘肅省白龍江流域苔類植物的研究[D].濟南:山東師范大學,2009.

[17]舒和平,劉東飛,顧春杰,等. 中小尺度區域泥石流災害風險評價[J].山地學報,2015,32(6):754-760.

[18]陜西省地質局區測隊.武都幅區域地質圖及說明書(1:20萬)[R].西安:陜西省地質局區測隊,1970.

[19]胡向德,黎志恒,魏潔,等. 舟曲縣三眼峪溝特大型泥石流的形成和運動特征[J].水文地質工程地質,2011,38(4):82-87.

[20]馬東濤,祁龍.三眼峪溝泥石流災害及其綜合治理[J].水土保持通報,1997,17(4):26-31.

[21]張成檢.舟曲縣自龍江流域地質環境及地質災害分布特征[J].甘肅水利水電技術，2010,46(12):26-28.

[22]胡凱衡,葛永剛,崔鵬,等.對甘肅舟曲特大泥石流災害的初步認識[J].山地學報,2010,28(5):628-634.

[23]黎志恒,朱立峰,胡向德,等.三眼峪特大泥石流溝流域分區特征[J].西北地質,2011,44(3):38-43.

[24]牛最榮.建設水文應急監測隊伍提高水文應急監測能力-舟曲特大山洪泥石流水文應急監測思考[J].甘肅水利水電技術,2010,46(12):8-14.

[25]余斌.粘性泥石流的平均運動速度研究[J].地球科學進展,2008,23(5):524-531.

[26]陳洪凱,唐紅梅,陳野鷹.公路泥石流力學[M].北京:科學出版社,2007.

[27]章書成,OldrichHungr,OlavSlaymaker.泥石流中巨石沖擊力的計算[M].北京：科學出版社,1996:67-72.

[28]余斌,楊永紅,蘇永超,等.甘肅省舟曲8.7特大泥石流調查研究[J].工程地質學報,2010,18(4):437-443.

*收稿日期：2016-01-28修回日期：2016-03-07

基金項目：國家科技支撐計劃項目(2011BAK12B05)

第一作者簡介：舒和平(1989-),男,瑤族,湖南溆浦縣人,博士研究生,主要從事水文過程與泥石流災害方面的研究. E-mail: shuhp15@lzu.edu.cn

中圖分類號：S157.8; X43

文獻標志碼：A

文章編號：1000-811X(2016)03-0119-07

doi:10.3969/j.issn.1000-811X.2016.03.020

Analysis and Research on the Law of Motion of Debris Flows in Southern Gansu Province——A Case Study of Sanyanyu Gully in Zhouqu of Gansu Province

SHU Heping, HAN Tuo, QI Shi and MA Jinzhu

(KeyLaboratoryofWesternChina’sEnvironmentSystemsoftheMinitryofEducation,CollegeofEarthandEnvironmentalSciences,LanzhouUniversity,Lanzhou73000,China)

Abstract:Based on the similarity principle, we chose the appropriate scale to construct a debris flow model in Wudu district of Longnan City in the southern Gansu Province, it can rebuilt the “8.8” catastrophic debris flow in Zhouqu. Through physical simulation experiment, we can know the average velocity was 0.7417 m/s, the impact force is 1.99N, the thickness of the accumulation is 6.1 cm and the accumulation area is 3.64 m2. At the same time, we have also draw accumulation form. In addition, empirical formulas were selected to calculate debris flows dynamic values which are compared with the experimental data, and they can verify the accuracy of the model. Meanwhile, the motion law of debris flows is further analyzed. According to the elements of debris flow formation, hazard assessment factors are extracted. The result showed that Sanyanyu gully is a high hazard debris flow gully. This result can provide reliable scientific basis for debris flow disaster prevention project.

Key words:Sanyanyu gully; debris flow; model; hazard; impact force; flow velocity; Zhouqu of Gansu Province

舒和平，韓拓，齊識，等. 甘肅省南部泥石流運動規律的分析與研究——以甘肅舟曲三眼峪溝為例[J].災害學, 2016,31(3):119-124,137.[SHU Heping, HAN Tuo, QI Shi,et al.Analysis and Research on the Law of Motion of Debris Flows in Southern Gansu Province——A Case Study of Sanyanyu Gully in Zhouqu of Gansu Province[J].Journal of Catastrophology, 2016,31(3)：119-124,137.]