王薦霖

(成都理工大學(xué),四川 成都 610059)

泥石流災(zāi)害與人類工程活動(dòng)息息相關(guān),尤其是在中國(guó)西南地區(qū),每年因泥石流暴發(fā)造成的損失相當(dāng)巨大,這使得泥石流災(zāi)害的防治顯得尤為重要[1]。泥石流基本特征值作為泥石流災(zāi)害防治設(shè)計(jì)的重要依據(jù),采用合理的計(jì)算方法,對(duì)其結(jié)果的安全性與可靠性有著極其重要的意義。

1 工程概況

甸心洪溝位于云南省雙廊鎮(zhèn)長(zhǎng)育村。該泥石流溝谷中上游溝道較窄,縱坡降較陡,溝道形態(tài)呈現(xiàn)出典型的“V”字型。泥石流以沖刷為主,中下游溝道較寬,縱坡相對(duì)較小,以淤積沉積為主。

甸心洪溝溝長(zhǎng)約5.08 km,河源高程2 850 m,最低點(diǎn)位于洱海邊緣地帶,高程為1 964 m,最大高差886 m,溝道平均縱坡降180‰。溝口—溝頂為陡坡地形,坡度30°～40°,局部>40°,溝口以下為堆積區(qū),為緩坡地形,地形坡度一般為5°～15°,平均縱坡降68‰,流域總面積10.57 km2。

自上世紀(jì)80年代以來(lái),甸心洪溝多次發(fā)生泥石流災(zāi)害,由于近年來(lái)在堆積區(qū)上新建居民建筑及公路對(duì)原泥石流堆積區(qū)地貌改動(dòng)較大,很難調(diào)查泥石流沖出物質(zhì)堆積范圍及危害區(qū)范圍信息,通過(guò)對(duì)當(dāng)?shù)鼐用裨L問(wèn),其回憶最大一次泥石流導(dǎo)致較多房屋被沖毀,停靠路邊的拖拉機(jī)被沖入洱海,由此表明當(dāng)時(shí)泥石流規(guī)模應(yīng)達(dá)到中型以上,結(jié)合甸心洪溝泥石流溝道內(nèi)可移動(dòng)物源非常豐富,可移動(dòng)物源約14.90萬(wàn)m3,極端暴雨條件下可形成中型規(guī)模泥石流。綜上所述,甸心洪溝泥石流屬溝谷型暴雨型形成期低頻稀性泥石流。

2 泥石流基本特征值計(jì)算

由于此次甸心洪溝泥石流并無(wú)監(jiān)測(cè)資料，因此，其泥石流基本特征值的計(jì)算主要以野外調(diào)查獲取的泥痕、溝道斷面等特征為參照，并結(jié)合擬設(shè)工程的需要綜合獲取[2]。需計(jì)算獲取的泥石流基本特征值包括：泥石流流體重度及固體物質(zhì)密度、設(shè)計(jì)清水洪峰流量、泥石流流速、泥石流流量、一次泥石流過(guò)流總量、一次泥石流固體沖出物總量。

2.1 泥石流流體重度及固體密度

泥石流流體重度需采用形態(tài)調(diào)查法和現(xiàn)場(chǎng)重度配方法,并結(jié)合量化評(píng)分綜合確定[3]?，F(xiàn)場(chǎng)重度配方法采用《工程地質(zhì)手冊(cè)》(第四版)中的變種公式計(jì)算[4],如下:

式中:γc為泥石流重度(t/m3);Gc為配制泥漿重量(t);V為配制泥漿體積(m3)。

形態(tài)調(diào)查判定泥石流漿體為稠漿狀,容重為1.40～1.50 t/m3;現(xiàn)場(chǎng)重度配方法流體密度為1.44～1.49 t/m3,平均值為1.46 t/m3,由于受取樣的限制,較大碎塊石不能參與配漿試驗(yàn),所得出的試驗(yàn)值偏小,本次綜合后取泥石流容重為1.50 g/cm3?？紤]到該泥石流溝堆積物以粗顆粒為主,其原巖成分主要為灰?guī)r,根據(jù)經(jīng)驗(yàn),泥石流固體密度取2.60 t/m3。

2.2 設(shè)計(jì)清水洪峰流量

清水洪峰流量采用中國(guó)水利水電科學(xué)院提出的經(jīng)驗(yàn)公式及大理水文法兩種方法計(jì)算地表水匯水流量。

(1) 中國(guó)水利水電科學(xué)院經(jīng)驗(yàn)公式法。依據(jù)《泥石流防治工程技術(shù)》第三章中的相關(guān)公式進(jìn)行計(jì)算[5]。計(jì)算結(jié)果見(jiàn)表1。

匯流面積F≥3 km2時(shí),Qp=0.278φSF2/3;

匯流面積F<3 km2時(shí),Qp=0.278φSF。

式中:QP為頻率為P的暴雨洪峰流量;φ為產(chǎn)流系數(shù);S為設(shè)計(jì)1 h降雨強(qiáng)度(mm/h);F為各斷面匯水面積(km2)。

表1 甸心洪溝泥石流各斷面清水洪峰流量計(jì)算結(jié)果表(中國(guó)水利水電科學(xué)院經(jīng)驗(yàn)公式法)Table 1 Calculation result table of peak discharge in each section of debris flow in the meadow gully

(2) 大理水文法。計(jì)算公式參照《云南省水文手冊(cè)》(表2),計(jì)算結(jié)果見(jiàn)表3。

表2 大理水文法計(jì)算公式Table 2 Calculation formula of Dali hydrologic method

表3 甸心洪溝泥石流各斷面清水洪峰流量計(jì)算結(jié)果表(大理水文法)Table 3 Calculation result table of peak discharge in each section of debris flow in the meadow gully(Dali hydrologic method)

經(jīng)過(guò)兩種方法計(jì)算結(jié)果對(duì)比,雖然中國(guó)水利水電科學(xué)院經(jīng)驗(yàn)公式法計(jì)算結(jié)果更加保守,但大理的地方計(jì)算公式更能因地制宜地說(shuō)明該地區(qū)的實(shí)際情況,更符合該地區(qū)的實(shí)際現(xiàn)狀,故采用大理水文法計(jì)算清水洪峰流量。

2.3 泥石流流速

根據(jù)調(diào)查和泥石流重度計(jì)算,泥石流重度<1.6 t/m3,屬稀性泥石流。本次采用稀性泥石流流速計(jì)算公式(西南地區(qū))計(jì)算各斷面流速,即《泥石流災(zāi)害防治工程勘查規(guī)范》(DZ/T 0220—2016)中的公式[6]。

式中:Vc為泥石流平均流速(m/s);1/n為河床糙率系數(shù);φ為泥石流修正系數(shù);γH為泥石流固體物質(zhì)密度(t/m3);R為水力半徑,一般用平均泥深Hc(m)代替;Ic為泥石流水力坡度(‰),用溝床縱坡代替。

表4 泥石流流速計(jì)算成果表Table 4 Results table of velocity calculation of debris flow

2.4 泥石流流量

2.4.1 雨洪法

根據(jù)《泥石流災(zāi)害防治工程勘查規(guī)范》(DZ/T 0220—2016)中的公式進(jìn)行計(jì)算[6](表5)。

Qc=(1+φ)QBDc

式中:Qc為頻率為P的泥石流洪峰值流量(m3/s);QB為頻率為P的暴雨洪水設(shè)計(jì)流量(m3/s);Dc為泥石流堵塞系數(shù),根據(jù)溝谷特征查經(jīng)驗(yàn)表。

2.4.2 形態(tài)調(diào)查法

根據(jù)《泥石流災(zāi)害防治工程勘查規(guī)范》(DZ/T 0220—2016)中的公式進(jìn)行計(jì)算[6](表6)。

Qc=WcVc

式中:Qc為泥石流斷面峰值流量(m3/s);Wc為泥石流過(guò)流斷面面積(m2);Vc為泥石流斷面平均流速(m/s)。

對(duì)比兩種方法計(jì)算結(jié)果表,采用形態(tài)調(diào)查法求得的泥石流峰值流量普遍大于雨洪法(P=3.33%降雨頻率)計(jì)算求得的結(jié)果,總體上兩種方法的計(jì)算結(jié)果是相對(duì)應(yīng)的。但是,采用形態(tài)調(diào)查法計(jì)算結(jié)果沒(méi)有暴雨頻率的概念,僅能代表當(dāng)次泥石流的特征值,而雨洪法則根據(jù)現(xiàn)有溝域面積進(jìn)行計(jì)算,具有預(yù)測(cè)性質(zhì)。因此,泥石流峰值流量采用雨洪法。

2.5 一次泥石流過(guò)流總量

根據(jù)《泥石流災(zāi)害防治工程勘查規(guī)范》(DZ/T 0220—2016)中的公式進(jìn)行計(jì)算,一次泥石流總量Q(表7)。

表5 泥石流流量計(jì)算表(雨洪法)Table 5 Calculation table of discharge of debris flow(rain flood method)

表6 泥石流流量計(jì)算(形態(tài)調(diào)查法)Table 6 Calculation table of discharge of debris flow(morphologic investigation method)

Q=K×T×Qc

式中:當(dāng)流域面積F<5 km2時(shí),K=0.202;F=5～10 km2時(shí),K=0.113;F=10～100 km2時(shí),K=0.037 8;F>100 km2時(shí),K<0.025 2;Q為一次泥石流的總量(m3);T為泥石流歷時(shí)(s);Qc為泥石流的洪峰流量(m3/s)。

表7 一次泥石流過(guò)流總量計(jì)算表Table 7 Volume table of debris flow

根據(jù)計(jì)算結(jié)果,對(duì)比降雨標(biāo)準(zhǔn),上一次暴發(fā)泥石流時(shí),甸心洪溝40 min降雨強(qiáng)度頻率近30年一遇,所以該次泥石流強(qiáng)度頻率為3.30%,計(jì)算甸心洪溝泥石流沖出量約16 366.0 m3。

2.6 一次泥石流固體沖出物總量

一次泥石流固體沖出物按照《泥石流災(zāi)害防治工程勘查規(guī)范》(DT/T 0220—2016)進(jìn)行計(jì)算[6](表8)。

QH=Q(γc-γw)/(γH-γw)

式中:QH為一次泥石流沖出固體物質(zhì)總量(m3);Q為一次泥石流過(guò)程總量(m3);γc為泥石流重度(t/m3);γw為水的重度(t/m3);γH為泥石流固體物質(zhì)的重度(t/m3)。

根據(jù)計(jì)算,上一次暴發(fā)泥石流強(qiáng)度頻率按3.3%計(jì),甸心洪溝一次泥石流沖出固體物質(zhì)量約5 728.0 m3,與堆積區(qū)估算量基本一致。

表8 一次泥石流固體沖出物總量計(jì)算成果表Table 8 Results of the calculation of the total amount of solid debris flow in a debris flow

3 計(jì)算結(jié)果的設(shè)計(jì)應(yīng)用

在進(jìn)行泥石流災(zāi)害防治設(shè)計(jì)時(shí),經(jīng)專家探討,甸心洪溝泥石流采用以“攔擋”為主,并在下游輔以“排導(dǎo)、護(hù)岸”的防治思路。

在進(jìn)行谷坊壩、攔擋壩結(jié)構(gòu)設(shè)計(jì)時(shí),根據(jù)泥石流流體重度、流速等相關(guān)因素確定壩身尺寸,其壩體的庫(kù)容應(yīng)大于泥石流過(guò)流總量,且壩體溢流口的大小應(yīng)滿足泥石流的過(guò)流;在進(jìn)行下游雙邊護(hù)岸堤結(jié)構(gòu)設(shè)計(jì)時(shí),在滿足其自身結(jié)構(gòu)穩(wěn)定的前提下,同樣還需滿足泥石流的過(guò)流。

4 結(jié)語(yǔ)

(1) 通過(guò)工程實(shí)例,采用規(guī)范公式,經(jīng)相關(guān)專家審核過(guò)后,基本特征值的計(jì)算結(jié)果具有較高的安全性及可靠性,能夠?yàn)榈樾暮闇夏嗍鞯姆乐卧O(shè)計(jì)提供重要依據(jù)。依托該計(jì)算結(jié)果,可進(jìn)行可行性研究及施工圖設(shè)計(jì),能夠保護(hù)該地區(qū)居民約3 000余人,解除泥石流災(zāi)害對(duì)居民及建筑的威脅,達(dá)到明顯的減災(zāi)效益。

(2) 由于現(xiàn)場(chǎng)試驗(yàn)條件的制約、村民記憶的偏差等因素,會(huì)導(dǎo)致泥石流基本特征值計(jì)算結(jié)果有誤差。本文主要在于為泥石流基本特征值的計(jì)算提供思路、方法,在實(shí)際生產(chǎn)過(guò)程中,可用多次測(cè)量取平均值、多種方法對(duì)比綜合取值的方法,來(lái)盡量減小誤差。

[1] 李三明,張艷濤,湯江紅.泥石流特征值的計(jì)算與確定[J].資源環(huán)境與工程,2010,24(3):297-300.

[2] 張生祥,何明陽(yáng),謝光輝.地震災(zāi)區(qū)特大型泥石流基本特征值計(jì)算與分析[J].四川地質(zhì)學(xué)報(bào),2014,34(s2):36-40.

[3] Б.С.Степанов,孟河清.泥石流特征值的計(jì)算問(wèn)題[J].力學(xué)進(jìn)展,1987,17(4):98-106.

[4] 常士驃,張?zhí)K民.工程地質(zhì)手冊(cè)[M].第四版.北京:中國(guó)建筑工業(yè)出版社，2007.

[5] 王繼康.泥石流防治工程技術(shù)[M].北京:中國(guó)鐵道出版社,1996.

[6] 四川省國(guó)土資源廳.泥石流災(zāi)害防治工程勘察規(guī)范:DZ/T0220—2016[S].北京:中國(guó)標(biāo)準(zhǔn)出版社，2016.