何曉英，陳洪凱，唐紅梅，朱繡竹

(重慶交通大學(xué) 巖土工程研究所，重慶 400074)

0 引 言

龍頭是泥石流運(yùn)動(dòng)的一種特有屬性，龍頭的出現(xiàn)增加了泥石流對(duì)建筑物的破壞能力，預(yù)測(cè)泥石流龍頭的沖擊力是泥石流防治工程設(shè)計(jì)中一個(gè)重要的課題[1]。目前，關(guān)于泥石流龍頭的研究大多集中于泥石流龍頭的形態(tài)描述與形成機(jī)制的探索，如T.Ilstad，等[2]認(rèn)為泥石流龍頭后部一定區(qū)域可代表典型的黏塑性流，具有較高的孔隙壓力；陳春光，等[3]針對(duì)入?yún)R主河的泥石流龍頭運(yùn)動(dòng)機(jī)理，基于野外觀測(cè)資料和模型實(shí)驗(yàn)建立了潛入式交匯的泥石流龍頭運(yùn)動(dòng)方程，發(fā)現(xiàn)混合流黏滯系數(shù)和繞流系數(shù)對(duì)龍頭運(yùn)動(dòng)速度影響較大；王兆印[4]通過室內(nèi)泥石流實(shí)驗(yàn)研究泥石流龍頭的形成和發(fā)展，發(fā)現(xiàn)混石流龍頭隆起高度與龍頭卵石粒徑成正比，并建立了龍頭運(yùn)動(dòng)的能量理論解釋泥石流的形成與運(yùn)動(dòng)機(jī)理；胡凱衡，等[5]描述了蔣家溝黏性泥石流溝道里的龍頭-龍尾陣次交替的陣性波狀運(yùn)動(dòng)過程，將泥石流波狀運(yùn)動(dòng)分為緩波、滾動(dòng)短波、滾動(dòng)長(zhǎng)波等3類；陳洪凱，等[6]將泥石流龍頭的形成過程可概化為泥石流勻速流動(dòng)、差異流動(dòng)、初始龍頭形成和間歇龍頭形成4個(gè)階段，建立泥石流龍頭壓脹機(jī)理解釋了泥石流陣流、逆序結(jié)構(gòu)和顆粒浮動(dòng)等現(xiàn)象；Z.F.Haza，等[7]采用模型實(shí)驗(yàn)分析水下泥石流龍頭運(yùn)動(dòng)，得到泥石流龍頭流速在0.164～0.439 m/s之間浮動(dòng)，其最小與最大剪應(yīng)力分別為4.47 與127.72 Pa；T.Ilstad，等[8]通過對(duì)水下泥石流的室內(nèi)模型實(shí)驗(yàn)發(fā)現(xiàn)，對(duì)于黏性泥石流，其龍頭存在滑水現(xiàn)象，并討論了在不同重度，不同泥漿組成條件下泥石流龍頭的形態(tài)和運(yùn)動(dòng)規(guī)律；S.T.Lancaster，等[9]認(rèn)為泥石流龍頭部位攜帶有大量樹枝、碎石等物質(zhì)，造成泥石流龍頭運(yùn)動(dòng)速度加劇，侵蝕深度增加；魏鴻[10]通過穩(wěn)定均勻顆粒水石兩相流龍頭沖擊壩體的水槽實(shí)驗(yàn)分析了沖擊荷載峰值的構(gòu)成情況，認(rèn)為泥石流沖擊力與泥石流流速成正比；R.M.Iverson[11]認(rèn)為泥石流頭部聚集大量的粗大顆粒主要是由于分選較差的泥石流體與溝床之間的摩擦作用增強(qiáng)所致。為分析泥石流龍頭沖擊力與泥石流固液兩相介質(zhì)之間的關(guān)系，筆者修筑較大尺度的泥石流溝實(shí)驗(yàn)?zāi)Ｐ?，通過對(duì)泥石流漿體黏度、固相粒徑及固相比等因素的科學(xué)組合，實(shí)施系列泥石流龍頭沖擊特性模型實(shí)驗(yàn)，分析泥石流龍頭沖擊荷載曲線，揭示實(shí)驗(yàn)條件下黏性泥石流龍頭沖擊特性。

1 泥石流龍頭運(yùn)動(dòng)機(jī)理

假設(shè)泥石流為固液兩相流體，將其等效為漿體(水和細(xì)顆粒泥沙摻混而成)和固體(泥石流體中較粗的泥沙顆粒)，如圖1。

圖1 泥石流兩相等效模型Fig.1 Two-phase equivalent model of debris flow

假定泥石流運(yùn)動(dòng)為一維兩相流運(yùn)動(dòng)體系(圖2)，泥石流運(yùn)動(dòng)中，泥石流龍頭受到自身重力、溝底摩擦力、溝床坡度等因素的影響，則溝床上的泥石流體的受力情況如下。

泥石流龍頭所受推力F為：

F=Gsinθ+Gmcosθ

(1)

G=gh2[α(ρs-ρf)+ρf]

(2)

(3)

式中：θ為溝床坡度，(°)；G，Gm分別為泥石流龍頭自重應(yīng)力與后續(xù)流推力，Pa；h為泥石流龍頭高度，m；ρs，ρf分別為泥石流固相與液相的密度，kg/m3；α，α1分別為泥石流龍頭濃度與后續(xù)流濃度。

泥石流龍頭所受阻力FR：

FR=(μsGs+μfGf)cosθ

(4)

式中：μs，μf分別為固相顆粒與液相漿體與泥石流溝床相互作用的摩擦系數(shù)。

圖2 泥石流龍頭受力Fig.2 Force of debris flow heads

同時(shí)，泥石流龍頭在運(yùn)動(dòng)過程中固相顆粒聚集，大量實(shí)驗(yàn)資料表明泥石流固相運(yùn)動(dòng)速度遠(yuǎn)小于液相漿體的運(yùn)動(dòng)速度[12]，液相漿體對(duì)固相的運(yùn)動(dòng)具有拖曳作用，其拖曳速度為(vf-vs)，泥石流固液兩相的差異性流動(dòng)越明顯，泥石流龍頭更易出現(xiàn)[13]。單位泥石流體受力情況如圖3。

圖3 泥石流龍頭單元體物相相互作用Fig.3 Interacting force between solid phase and fluid phase in element debris flow heads

Ff—漿體之間的黏性應(yīng)力；Fsf，F(xiàn)fs—漿體與固體之間的拖曳

力及摩阻力(二者大小相同，方向相反)；Fs—相顆粒間摩擦力

泥石流液相漿體的黏性應(yīng)力根據(jù)Binhanm的流變方程計(jì)算：

(5)

式中：τB為漿體的Binhanm極限屈服應(yīng)力，Pa；μ為漿體的黏度系數(shù)，Pa·s；ρf為泥石流液相的密度，kg/m3；L為泥石流體中流層之間的混摻長(zhǎng)度，即液相漿體微團(tuán)因脈動(dòng)流速的作用而移動(dòng)的距離，m；dvf/dy為泥石流漿體豎向流速梯度。

根據(jù)以上分析可知，泥石流運(yùn)動(dòng)過程與泥石流固相與液相介質(zhì)的性質(zhì)(液相漿體黏度、固液兩相介質(zhì)密度)泥石流固相與液相介質(zhì)之間的相互作用(液相拖曳力，固相摩阻力、固相摩擦力)及泥石流固相與液相介質(zhì)與泥石流溝床的摩擦力有關(guān)。在這些作用的影響下，泥石流在運(yùn)動(dòng)過程中，泥石流體內(nèi)部固液兩相介質(zhì)相互作用，易形成壓脹核；泥石流體外部，由于固液兩相介質(zhì)與溝床的摩擦系數(shù)不同，固相顆粒所受摩擦力大于液相所受摩擦力，固液兩相介質(zhì)之間流速產(chǎn)生差異，形成差異性流動(dòng)，固相顆粒聚集，泥石流端部產(chǎn)生壅高，龍頭形成。

2 實(shí)驗(yàn)?zāi)Ｐ驮O(shè)計(jì)

2.1 實(shí)驗(yàn)?zāi)Ｐ?/h3>

以四川涼山州平川泥石流為原型，考慮尺寸相似與地形相似，修筑泥石流實(shí)驗(yàn)?zāi)Ｐ?。泥石流主溝槽長(zhǎng)15.0 m，平均溝槽寬度0.5 m，平均深度0.6 m，高差3.6 m；設(shè)置兩個(gè)彎道，坡度分別為27°，10°。溝槽頂端為1.5 m3的水源箱，出口放置體積為0.26 m3的木質(zhì)結(jié)構(gòu)碎石漏斗。泥石流溝口設(shè)置動(dòng)態(tài)應(yīng)力傳感器，接受泥石流沖擊信號(hào)(圖4、圖5)。

圖4 實(shí)驗(yàn)?zāi)Ｐ驮O(shè)計(jì)平面(單位：m)Fig.4 Designing plane figure of the testing model

圖5 泥石流運(yùn)動(dòng)與沖擊特性實(shí)驗(yàn)?zāi)Ｐ虵ig.5 Testing model of the debris flow movement and impacting properties

2.2 實(shí)驗(yàn)工況與材料

選取5組泥石流漿體黏度η(單位：Pa·s)(0，0.13，0.37，0.72，0.93)、5組固相比r(0，0.02， 0.04，0.08，0.16)，及4組顆粒粒徑d(cm)(河沙0.075～0.2，瓜米石0.3～0.5，碎石0.8～1.0、1.3～1.5)進(jìn)行正交組合，則實(shí)驗(yàn)工況為：Cηrd(η= 1，2，3，4，5；r= 1，2，3，4，5；d= 1，2，3，4)，共85組。

泥石流漿體黏度采用淀粉醚、熟膠粉與水配比而成，采用NDJ-1旋轉(zhuǎn)黏度計(jì)測(cè)定配制后的泥石流漿體，使其達(dá)到預(yù)期黏度要求。固相比為固相顆粒體積與泥石流體的體積之比。

2.3 觀測(cè)儀器及實(shí)驗(yàn)過程

動(dòng)態(tài)應(yīng)力傳感器，線性精度為5‰；DH5922動(dòng)態(tài)信號(hào)測(cè)試系統(tǒng)，采樣頻率設(shè)定為1000 Hz；Fastcamultimal024型高精度攝像儀，測(cè)試精度400 fpt。實(shí)驗(yàn)過程如下：蓄水池蓄滿水，加入相應(yīng)配比的淀粉醚與熟膠粉，攪拌均勻直至反應(yīng)完全；按固相比將不同粒徑的固相顆粒倒入碎石漏斗；打開水源箱閥門與碎石漏斗閥門，控制固相顆粒加入速度，確保其均勻分布于泥石流漿體中；高速攝像機(jī)連續(xù)記錄泥石流運(yùn)動(dòng)過程，傳感器接收泥石流沖擊，動(dòng)態(tài)信號(hào)測(cè)試系統(tǒng)連續(xù)記錄泥石流沖擊荷載變化過程。

3 實(shí)驗(yàn)結(jié)果分析

3.1 泥石流龍頭沖擊力

根據(jù)攝像機(jī)拍攝到的圖像，甄別出泥石流龍頭形態(tài)，從動(dòng)態(tài)信號(hào)測(cè)試系統(tǒng)中讀取泥石流龍頭部分沖擊力峰值、平均值及沖擊作用時(shí)間，將得到的數(shù)據(jù)與泥石流整體沖擊力進(jìn)行對(duì)比，結(jié)果如表1。

表1 泥石流龍頭沖擊力 Table 1 Impacting force of debris flow heads

分析表1實(shí)驗(yàn)數(shù)據(jù)知：

1)隨著漿體黏度的增大，泥石流龍頭作用越不明顯，當(dāng)η≥0.72 Pa·s時(shí)，泥石流沖擊過程中不再出現(xiàn)明顯的龍頭。

2)泥石流龍頭沖擊作用時(shí)間隨著液相漿體黏度的增大而減小。當(dāng)η= 0 Pa·s時(shí)，泥石流龍頭沖擊時(shí)間平均值為7.4 s；當(dāng)η= 0.13 Pa·s時(shí)，泥石流龍頭沖擊時(shí)間平均值為2.6 s；當(dāng)η= 0.37 Pa·s時(shí)，泥石流龍頭沖擊時(shí)間平均值為1.6 s。

3)在漿體黏度相同的條件下，泥石流龍頭作用時(shí)間隨著固相比增大而減小。

4)泥石流龍頭峰值在低黏度與低固相比組合條件下為沖擊全程峰值(如η= 0.13 Pa·s，r= 0.02，0.04時(shí)、η= 0.37 Pa·s ，r= 0.02時(shí))，但隨著黏度與固相比的增大，龍頭部分有峰值但不是全程峰值。

3.2 泥石流龍頭沖擊形態(tài)

分析實(shí)驗(yàn)中高速攝像儀與DH5922動(dòng)態(tài)測(cè)試系統(tǒng)獲得的85種工況下泥石流沖擊圖像與沖擊數(shù)字信號(hào)，并結(jié)合表1的龍頭沖擊力數(shù)據(jù)，可將泥石流龍頭形態(tài)分為3種形式。

3.2.1 高速高強(qiáng)度龍頭形態(tài)

該類龍頭高度較大，呈飛濺狀，龍頭速度最大，而后迅速向尾部減??；龍頭在接觸傳感器的短時(shí)間內(nèi)出現(xiàn)較龍身及尾部都高的壓力峰值，常出現(xiàn)泥石流龍頭飛離龍身，出現(xiàn)短暫斷流現(xiàn)象(圖6)。該類形態(tài)一般出現(xiàn)在低黏度(η= 0，0.13 Pa·s時(shí))、低固相比實(shí)驗(yàn)工況下，例如C222及C224工況(圖7)。

圖6 高速高強(qiáng)度龍頭形態(tài)Fig.6 High speed-high intensity heads configuration

圖7 C222，C224沖擊荷載Fig.7 Impacting force figure of work condition C222 and C224

3.2.2 過渡型龍頭形態(tài)

該類龍頭高度較龍身仍較高，但趨于平順，沒有飛濺的泥石流體脫離泥石流整體的運(yùn)動(dòng)，龍頭速度仍然較大，但是與龍身流速差異減小(圖8)。

圖8 過渡型龍頭形態(tài)Fig.8 Transitional type heads configuration

龍頭部位的沖擊力有明顯峰值，但不是泥石流整體沖擊力的峰值。該類形態(tài)一般出現(xiàn)在低黏度高固相比(η= 0，0.13 Pa·s，r>0.08時(shí)，如C241工況)、中黏度低固相比(η= 0.37 Pa·s，r<0.16時(shí)，如C322)工況下(圖9)。

圖9 C241，C322工況沖擊荷載Fig.9 Impacting force figure of work condition C241 and C322

3.2.3 低速低強(qiáng)度形態(tài)

該類龍頭低矮，流體形態(tài)平順，流速較小，泥石流紊動(dòng)較弱，不存在紊動(dòng)虛擬高度。龍頭至龍尾光滑平順(圖10)。

圖10 低速低強(qiáng)度形態(tài)Fig.10 Low speed-low intensity heads configuration.

龍頭部位的沖擊力有一定紊動(dòng)但不存在明顯峰值。該類形態(tài)是本次實(shí)驗(yàn)的大部分情況，出現(xiàn)在中黏度高固相比(η= 0.37 Pa·s，r> 0.08時(shí)，如C333)及高黏度(η= 0.72 ，0.93 Pa·s，如C542)工況下，且黏度與固相比越大，龍頭形態(tài)越低矮，紊動(dòng)越不明顯，泥石流龍頭沖擊力呈平滑上升趨勢(shì)(圖11)。

圖11 C333，C542工況沖擊荷載Fig.11 Impacting force figure of work condition C333 and C542

4 結(jié) 論

將泥石流等效為兩相流體系，分析泥石流龍頭內(nèi)部固液兩相介質(zhì)之間的相互作用探討了泥石流龍頭形成機(jī)理。并以平川泥石流溝為原型自行研制泥石流運(yùn)動(dòng)與沖擊特性實(shí)驗(yàn)?zāi)Ｐ?，探索了不同黏度、不同固相比、不同固相顆粒粒徑組合條件下泥石流龍頭沖擊特性，實(shí)施了85組泥石流沖擊實(shí)驗(yàn)，獲取了5 084 500多個(gè)測(cè)試數(shù)據(jù)。分析得到各工況條件下泥石流龍頭沖擊形態(tài)及沖擊荷載變化規(guī)律：

1)當(dāng)黏度大于0.72 Pa·s時(shí)，泥石流運(yùn)動(dòng)中不產(chǎn)生明顯的龍頭，泥石流龍頭作用時(shí)間隨著漿體黏度與固相顆粒固相比的增加而減??；

2)龍頭沖擊力峰值在低黏度與低固相比組合條件下為沖擊全程峰值但隨著黏度與固相比的增大，龍頭沖擊力峰值減小，不是全程峰值，據(jù)此將泥石流龍頭形態(tài)分為高速低強(qiáng)度龍頭形態(tài)、過渡型龍頭形態(tài)及低速高強(qiáng)度形態(tài)。

[1] 沈壽長(zhǎng),李良勛,魏鴻.稀性泥石流龍頭形成機(jī)制的實(shí)驗(yàn)研究[J].中國鐵道科學(xué),1996,17(3):42-49.Sheng Shouchang,Li Liangxun,Wei Hong.Experiments on the surge of uncoherent debris flow [J].China Railway Science,1996,17(3):42-49.

[2] Ilstad T,Elverhφi A, Issler D,et al.Sub-aqueous debris flow behavior and its dependence on the sand/clay ratio:a laboratory study using particle tracking [J].Marine Geology,2004,213:4l5-438.

[3] 陳春光,姚令侃,楊慶華.入?yún)R主河的泥石流龍頭運(yùn)動(dòng)機(jī)理研究[J].水利學(xué)報(bào),2004(1):11-16.Chen Chunguang,Yao Lingkan,Yang Qinghua.Movement of debris flow head in main channel at the confluence [J].Journal of Hydraulic Engineering,2004 (1):11-16.

[4] 王兆印.泥石流龍頭運(yùn)動(dòng)的實(shí)驗(yàn)研究及能量理論[J].水利學(xué)報(bào),2001(3):18-26. Wang Zhaoyin.Experimental study oil debris flow head and the energy theory [J].Journal of Hydraulic Engineering,2001(3):18-26.

[5] 胡凱衡,韋方強(qiáng),李泳,等.泥石流的陣性波狀運(yùn)動(dòng)[J].山地學(xué)報(bào),2004 (11):707-712.Hu Kaiheng,Wei Fangqiang,Li Yong,et al.Characteristics of debris-flow surge [J].Journal of Mountain Research,2004 (11):707-712.

[6] 陳洪凱,杜榕桓,唐紅梅,等.泥石流龍頭壓脹機(jī)理探析[J].重慶交通大學(xué)學(xué)報(bào):自然科學(xué)版,2008,27(5):790-793. Chen Hongkai,Du Ronghuan,Tang Hongmei,et al.Study on blasting mechanism of high pressure core in debris flow [J].Journal of Chongqing Jiaotong University:Natural Science,2008,27(5):790-793.

[7] Haza Z F,Harahap I S H,Dakssa L M.Experimental studies of the flow-front and drag forces exerted by sub-aqueous mudflow on inclined base [J].Nat Hazards,2013,3:643-651.

[8] Ilstada T,Elverhbia A,Isslera D,et al.Sub-aqueous debris flow behavior and its dependence on the sand/clay ratio: a laboratory study using particle tracking [J].Marine Geology,2004,213:415-438.

[9] Lancaster S T,Hayes S K.Effects of wood on debris flow run-out in small mountain watersheds [J].Water Resources Research,2003,39 (6):1168-1188.

[10] 魏鴻.泥石流龍頭對(duì)壩體沖擊力的實(shí)驗(yàn)研究[J].中國鐵道科學(xué),1996,17(3):50-62.Wei Hong.Experimental study on impact force of debris flow heads [J].China Railway Science,1996,17(3):50-62.

[11] Iverson R M.Dynamic pore-pressure fluctuations in rapidly shearing granular materials [J].Advancement of Science A,1989,246:796-799.

[12] 錢寧.高含沙水流運(yùn)動(dòng)[M].北京:清華大學(xué)出版社,1989:156-163.Qian Ning.The Hyper-concentrated Sediment Flow [M].Beijing:Tsinghua University Press,1989:156-163.

[13] Chen Hongkai,Tang Hongmei,Chen Yeying,et al.Research on method to calculate velocities of solid phase and liquid phase in debris flow [J].Applied Mathematics and Mechanics:English Edition,2006,27(3):357-364.