李姝景 劉曉峰 武立功 肖利興 姚 池 姜清輝1

（1.南昌大學(xué)建筑工程學(xué)院，江西南昌330031；2.江西省尾礦庫工程安全重點(diǎn)實(shí)驗(yàn)室，江西南昌330031；3.青州市應(yīng)急管理局，山東濰坊262500；4.中國(guó)瑞林工程技術(shù)股份有限公司，江西南昌330031）

尾礦庫一旦潰壩，壩內(nèi)的尾礦砂、水、泥將以泥石流的形式流出，會(huì)嚴(yán)重威脅下游居民的生命財(cái)產(chǎn)安全，因此對(duì)尾礦庫下游進(jìn)行有效防控很有必要［1-4］。

目前，關(guān)于尾礦庫潰壩特征及演進(jìn)的研究方法主要有物理模型試驗(yàn)和數(shù)值模擬兩種［5］。陶冬良［6］以鉬礦尾礦庫潰壩事件為工程背景，利用模型試驗(yàn)法，明確了潰壩事故發(fā)生時(shí)的安全撤離時(shí)間以及尾礦庫潰壩形態(tài)和潰壩洪水狀態(tài)的影響因素。張力霆等［7-8］基于某尾礦庫工程，運(yùn)用模型試驗(yàn)法和數(shù)值模擬法，研究了尾礦庫的潰壩機(jī)理及過程，總結(jié)了潰壩破壞模式。王儀心等［9］以山西襄汾尾礦庫潰壩事故為例，介紹了尾礦庫潰壩安全風(fēng)險(xiǎn)評(píng)價(jià)方法，并定量評(píng)價(jià)了其風(fēng)險(xiǎn)程度。為研究尾礦庫逐漸潰壩時(shí)潰口的變化形態(tài)以及下游淹沒范圍，李火坤等［10］利用FLOW-3D軟件建立了潰壩數(shù)值模型，分析了各位置淹沒深度及變化情況，驗(yàn)證了模擬的可行性。李鋼等［11］按一定的相似比尺建立了物理試驗(yàn)?zāi)Ｐ停治隽宋驳V庫潰壩后泥石流的演進(jìn)規(guī)律以及沖擊力，并將模型試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果進(jìn)行了對(duì)比，驗(yàn)證了所建模型的可靠性。王學(xué)良等［12］分別利用Massflow和PFC3D兩種方法模擬了尾礦庫潰壩后水砂的運(yùn)動(dòng)特征，并進(jìn)行了對(duì)比分析。王光進(jìn)等［13］利用模型試驗(yàn)方法進(jìn)行了尾礦庫漫頂潰壩分析，并探究了潰壩模式。廖威林等［14］基于ANSYS CFX軟件對(duì)尾砂的流變特性進(jìn)行了研究，結(jié)果表明：尾砂顆粒的粗細(xì)、濃度均影響剪切應(yīng)力大小。敬小非等［15］以秧田箐尾礦庫為研究對(duì)象，通過改變潰壩時(shí)的潰口形態(tài)，分析了下泄泥沙的流動(dòng)特性，認(rèn)為下泄泥沙的沖擊力、演進(jìn)規(guī)律以及淹沒范圍均受到潰口形態(tài)的影響。陳俊等［16］利用FLOW-3D軟件建立了潰壩數(shù)值模型，分析了潰壩后庫內(nèi)泥沙對(duì)下游產(chǎn)生的影響。陳小玉［17］研究了尾礦泥沙流變參數(shù)的計(jì)算方法，并分析了尾礦庫發(fā)生潰壩后下泄泥沙的演進(jìn)規(guī)律。

近年來，學(xué)術(shù)界對(duì)尾礦庫潰壩泥沙的成因及內(nèi)力的研究較多，對(duì)下游防護(hù)措施的研究較少［18］?，F(xiàn)有的尾礦庫潰壩特征研究主要集中在潰壩砂流的演進(jìn)特征、潰壩形式等方面，對(duì)提前設(shè)置防護(hù)措施以減少下游居民生命財(cái)產(chǎn)損失等方面的研究涉及較少。采取合理的防護(hù)措施對(duì)于尾礦庫工程防災(zāi)意義重大，尤其是下游開闊地形范圍廣、影響大的尾礦庫工程，相關(guān)防護(hù)的研究更為迫切。本研究以具有開闊下游區(qū)域的永平銅礦尾礦庫為例，建立下游設(shè)置不同形式攔擋壩的物理模型，分析不同形式攔擋壩對(duì)潰壩砂流的影響，并基于FLOW-3D軟件對(duì)該模型進(jìn)行三維數(shù)值模擬，計(jì)算并分析潰壩后下游淹沒范圍以及不同形式攔擋壩的攔擋效果，為尾礦庫潰壩泥石流防控研究提供參考。

1 模型試驗(yàn)

1.1 工程背景

江西省永平銅礦5座初期壩中1#壩的下游部分區(qū)域存在大量的居民區(qū)與農(nóng)田，如圖1所示。該類現(xiàn)象也普遍存在于國(guó)內(nèi)其他尾礦庫中，如湖北省楊家灣尾礦庫。此類尾礦庫發(fā)生潰壩時(shí)，會(huì)威脅到下游居民的生命財(cái)產(chǎn)安全。為研究開闊地形下的潰壩砂流演進(jìn)過程并確定下游最優(yōu)防護(hù)措施，本研究通過在實(shí)驗(yàn)室內(nèi)搭建尾礦庫潰壩試驗(yàn)三維模型，對(duì)下游尾砂淤積厚度以及潰壩砂流的影響范圍進(jìn)行分析，探究砂流影響范圍以及下游泥沙淤積厚度的變化規(guī)律，在此基礎(chǔ)上，對(duì)下游防護(hù)壩的布置形式進(jìn)行優(yōu)選。

1.2 試驗(yàn)方案

1.2.1 材料選取

為了保證尾砂的化學(xué)性質(zhì)相似，試驗(yàn)所用的模型砂取自永平燕倉尾礦庫。選取的尾砂粒徑均小于0.5 mm，顆粒級(jí)配曲線對(duì)比如圖2所示。尾砂具體物理性質(zhì)及力學(xué)參數(shù)取值見表1。

1.2.2 試驗(yàn)裝置

試驗(yàn)?zāi)Ｐ褪怯芍茲{系統(tǒng)、貯漿系統(tǒng)、排漿系統(tǒng)以及數(shù)據(jù)測(cè)量系統(tǒng)組成。其中，貯漿系統(tǒng)尺寸為80 cm×50 cm×100 cm（長(zhǎng)×寬×高），系統(tǒng)上方和前方均為開放的，在水槽內(nèi)設(shè)簡(jiǎn)易的斜坡裝置，使槽內(nèi)后方略高于潰口處，以便減少水槽內(nèi)尾砂淤積，確保能下泄充足的尾砂量，形成尾礦庫庫區(qū)模型。在貯漿系統(tǒng)前方開口處插入塑料擋板，便于儲(chǔ)存足夠的尾砂，達(dá)到瞬間潰壩的目的。為防止槽內(nèi)泥沙壓力過大造成擋板變形，在擋板前嵌入兩根橫向鋼筋來固定擋板，由此構(gòu)成排漿系統(tǒng)。模型試驗(yàn)裝置如圖3所示。

1.2.3 攔擋壩確定

為探究不同形式攔擋壩對(duì)潰壩下泄泥沙流的影響，從而確定最優(yōu)攔擋角度，試驗(yàn)采用了3種形式的壩體，夾角分別為 180°（工況1）、135°（工況 2）、90°（工況3），并且將壩體均設(shè)在距潰口1.5 m處。

1.2.4 試驗(yàn)過程

參照《尾礦庫手冊(cè)》［19］以及課題組相關(guān)研究經(jīng)驗(yàn)，確定試驗(yàn)所需水砂混合物體積為10 m3，水砂濃度比為1∶1，通過制漿系統(tǒng)將尾礦砂和水按照試驗(yàn)濃度比攪拌均勻，完成潰壩泥沙的初步制備，通過泄沙口將泥沙排放入貯沙池以模擬高勢(shì)能待潰頭頂尾礦庫，迅速打開排漿裝置，高勢(shì)能泥沙向下游開闊區(qū)流出，潰壩砂流演進(jìn)開始。同時(shí)，利用高速攝像機(jī)記錄砂流在開口區(qū)的演進(jìn)情況，通過數(shù)據(jù)量測(cè)系統(tǒng)記錄下游尾砂淤積厚度、尾砂影響范圍等數(shù)據(jù)。

1.3 試驗(yàn)結(jié)果及分析

1.3.1 不同攔擋壩形式下的下游泥沙淤積厚度特征

泥沙淤積厚度是研究泥沙演進(jìn)規(guī)律的一個(gè)重要參數(shù)，在本次試驗(yàn)中，通過分析泥沙淤積厚度可以直接體現(xiàn)攔擋壩對(duì)于泥沙淤積的攔截效果［18］。試驗(yàn)對(duì)泥沙淤積厚度、淹沒范圍進(jìn)行了精確測(cè)量，以潰口中心為起始點(diǎn)，橫向間隔10 cm、縱向間隔20 cm設(shè)置監(jiān)測(cè)點(diǎn)，并用MATLAB編制程序繪制出泥沙淤積厚度等高線云圖，進(jìn)一步分析3種壩型對(duì)泥沙演進(jìn)的影響。

為分析尾礦庫潰壩后尾砂下泄后的泥沙淤積情況，試驗(yàn)選取了兩個(gè)具有代表性的監(jiān)測(cè)區(qū)域，分別是潰口前20 cm×40 cm區(qū)域和攔擋壩前區(qū)域，如圖4所示。

3種工況下泥沙淤積厚度云圖如圖5所示。

分析圖5（a）可知：潰口前監(jiān)測(cè)區(qū)域的平均泥沙淤積厚度為2.5 cm，壩前泥沙平均淤積厚度為2 cm。

分析圖5（b）可知：潰口前監(jiān)測(cè)區(qū)域的平均泥沙淤積厚度為1.7 cm，壩前泥沙平均厚度為1.4 cm，相對(duì)于工況1，淹沒范圍內(nèi)各位置的淤積厚度均減少了。原因主要為：①下游設(shè)置的攔擋壩壩體可以使流動(dòng)的泥沙發(fā)生回流，泥沙從潰口下泄出來時(shí)會(huì)帶有較大的沖擊力，當(dāng)壩體夾角為180°時(shí)，整個(gè)攔擋壩會(huì)同時(shí)受到下泄泥沙的碰撞，此時(shí)由于攔擋壩的截流攔擋作用，使得泥沙發(fā)生較大范圍回流。由于泥沙碰撞壩體時(shí)間基本一致，所以壩前泥沙淤積厚度較深且均勻。壩體夾角為135°時(shí)，折角的出現(xiàn)使下泄泥沙碰撞到壩體各個(gè)位置存在時(shí)間差，以致壩前泥沙淤積厚度分布不均勻。當(dāng)泥沙碰撞折角兩側(cè)的壩體時(shí)，增大了泥沙運(yùn)動(dòng)位移，減少了壩前泥沙淤積量，由此也使下泄泥沙中的動(dòng)能減小，運(yùn)動(dòng)速度減小，回流范圍減小。②攔擋壩具有導(dǎo)流作用，壩體夾角為135°的攔擋壩的導(dǎo)流作用比夾角為180°的攔擋壩明顯，此時(shí)下泄泥沙在碰撞攔擋壩后一部分會(huì)沿著攔擋壩的方向流動(dòng)，使得壩前泥沙淤積厚度逐漸減小。對(duì)比以上兩種工況可知，減小攔擋壩的壩體夾角有利于減緩下泄泥沙對(duì)攔擋壩的沖擊力以及壩前淤積厚度，有利于壩體穩(wěn)定［20］。

分析圖5（c）可知：潰口前監(jiān)測(cè)區(qū)域泥沙平均淤積厚度為1.5 cm，壩前泥沙平均淤積厚度為1 cm，相對(duì)于工況2，各位置的淤積厚度均略有減小。進(jìn)一步分析發(fā)現(xiàn)，泥沙運(yùn)動(dòng)到135°壩體折角時(shí)沖擊力較90°工況更大，所以回流的范圍也比壩體夾角為90°的工況影響大。原因是，當(dāng)下游設(shè)置的攔擋壩的壩體夾角為90°時(shí)，壩體夾角過小，使得沖擊力不夠大，并且當(dāng)泥沙運(yùn)動(dòng)到兩側(cè)壩體時(shí)，運(yùn)動(dòng)速度已經(jīng)趨近于零，此時(shí)壩體僅是起著攔擋以及使泥沙沿壩體順勢(shì)流動(dòng)的作用，經(jīng)過一段時(shí)間后，會(huì)使壩前泥沙堆積過多，不利于壩體穩(wěn)定；當(dāng)下游設(shè)置的攔擋壩的壩體夾角為135°時(shí)，下泄泥沙到達(dá)壩體折角時(shí)的沖擊力比工況3大，雖然沖擊力未達(dá)到能帶動(dòng)大量泥沙發(fā)生回流的程度，但也可使泥沙因回流沿壩體向周圍擴(kuò)散，可以有效地下泄泥沙排到兩側(cè)［20］。

根據(jù)上述分析可知，因尾礦庫發(fā)生瞬潰時(shí)泥沙的流動(dòng)特性，淤積深度最深的區(qū)域?yàn)闈⒖谡胺健檫M(jìn)一步分析潰口前方泥沙淤積厚度的變化規(guī)律，在模型試驗(yàn)中增加了8個(gè)監(jiān)測(cè)點(diǎn)，如圖3所示。監(jiān)測(cè)點(diǎn)在不同夾角攔擋壩下的泥沙淤積厚度變化特征如圖6所示。分析該圖可知：攔擋壩的壩體夾角越小，壩前泥沙淤積厚度越小。當(dāng)下游設(shè)置壩體夾角為180°的攔擋壩時(shí)，壩前50 cm內(nèi)的泥沙淤積厚度均勻下降，而壩體夾角為135°和90°的攔擋壩，壩前10 cm處的泥沙淤積厚度會(huì)突增，這是由泥沙的回流作用以及回流后的剩余沖擊力所致。

1.3.2 不同攔擋壩形式下的下游泥沙淹沒范圍特征

淹沒范圍可直觀地反映攔擋壩的防護(hù)效果［21］。試驗(yàn)中提取高速相機(jī)拍攝的照片，并進(jìn)行圖像校正，通過計(jì)算分析系統(tǒng)獲取了各壩型攔擋作用下的尾砂淹沒范圍和淹沒面積，分別如圖7、表2所示。分析可知：在壩體夾角為135°的工況下，下游淹沒范圍最大，是由于攔擋壩既具有橫向?qū)Я髯饔?，又具有縱向?qū)Я髯饔?，兼具?dǎo)流和攔擋的作用，而其他兩種工況下攔擋壩的導(dǎo)流作用只有單向。

1.3.3 物理模型試驗(yàn)結(jié)果

綜上分析可知：夾角為180°的攔擋壩具有攔擋作用，下泄泥沙帶有較大的勢(shì)能，碰撞壩體后發(fā)生回流，增大了壩前泥沙淤積厚度，這需要壩體具有良好的抗沖擊能力和承載力。夾角為135°的攔擋壩具有攔擋和導(dǎo)流雙重作用，壩體上的折角減小了泥沙對(duì)壩體的沖擊力，阻擋了壩前泥沙的回流范圍。該壩型的雙向?qū)Я髯饔檬沟孟滦鼓嗌硵U(kuò)大了運(yùn)動(dòng)區(qū)域，能夠使下泄泥沙沿壩體排到周圍，減小了泥沙淤積厚度，降低了對(duì)壩體的危害，同時(shí)也保護(hù)了下游居民安全。夾角為90°的攔擋壩也具有攔擋和導(dǎo)流作用，但縱向?qū)Я髯饔么笥跀r擋和橫向?qū)Я髯饔?，以致泥沙的縱向流動(dòng)范圍較明顯，并且此工況下的下泄泥沙對(duì)攔擋壩的沖擊力比工況2小，所以當(dāng)泥沙觸碰壩體后減小了回流范圍，大部分泥沙都沿著壩體流向下游，但由于速度較小，泥沙長(zhǎng)時(shí)間的流出會(huì)導(dǎo)致壩前淤積較大，不利于壩體穩(wěn)定?？梢?，在下游設(shè)置壩體夾角為135°左右的攔擋壩，效果較好。

2 數(shù)值模擬及結(jié)果分析

2.1 數(shù)值模擬原理

本研究采用流體計(jì)算軟件FLOW-3D對(duì)該尾礦庫在設(shè)置不同攔擋壩時(shí)的潰壩砂流演進(jìn)進(jìn)行模擬，選擇FLOW-3D軟件中的泥沙顆粒流模塊，并結(jié)合FAVOR技術(shù)劃分網(wǎng)格，計(jì)算泥沙的淤積厚度和移動(dòng)值。

在FLOW-3D軟件中，非黏性顆粒流中的剪切應(yīng)力由3種剪切力組成：①因顆粒碰撞產(chǎn)生的剪切應(yīng)力τi；②顆粒的附加黏性剪切應(yīng)力τv；③流體剪切應(yīng)力τf。本研究通過文獻(xiàn)［22］中的公式來計(jì)算非黏性顆粒流中的剪切應(yīng)力τ：

式中，ρs為顆粒密度，kg/m3；ρf為流體密度，kg/m3；D為球體直徑，m；λ為直徑與最小間隙之比，λ=D/S；dudy為混合物速度梯度；μf為流體的動(dòng)態(tài)黏滯系數(shù)；ΔR為Couette流間隙。

2.2 數(shù)值模擬模型建立及參數(shù)設(shè)定

本研究利用Civil-3D軟件繪制幾何模型，然后導(dǎo)入FLOW-3D軟件完成幾何模型構(gòu)建，如圖8所示，模型計(jì)算參數(shù)取值見表3。

2.3 數(shù)值模擬結(jié)果及分析

2.3.1 下游泥沙淤積厚度特征對(duì)比

3種工況下的泥沙淤積厚度數(shù)值模擬結(jié)果如圖9所示。分析圖9可知：攔擋壩的壩體夾角為180°時(shí)，潰口前監(jiān)測(cè)區(qū)域的泥沙平均淤積厚度為2.3 cm，壩前泥沙淤積厚度為2 cm；攔擋壩的壩體夾角為135°時(shí)，潰口前監(jiān)測(cè)區(qū)域泥沙平均淤積厚度為2.4 cm，壩前泥沙淤積厚度為1.7 cm；攔擋壩的壩體夾角為90°時(shí)，潰壩后潰口前監(jiān)測(cè)區(qū)域泥沙平均淤積厚度為2 cm，壩前泥沙淤積厚度為1.5 cm。

模型試驗(yàn)與數(shù)值模擬的各項(xiàng)數(shù)據(jù)對(duì)比如表4所示。對(duì)比發(fā)現(xiàn)：數(shù)值模擬得出的下游泥沙淤積厚度云圖、各位置泥沙淤積厚度與模型試驗(yàn)結(jié)果相近。

3種壩型的壩前指定某個(gè)點(diǎn)的泥沙淤積厚度隨時(shí)間的變化特征如圖10所示。分析可知：最后趨于穩(wěn)定的壩前泥沙淤積厚度與模型試驗(yàn)數(shù)據(jù)相似，并且壩前泥沙淤積厚度隨著攔擋壩壩體夾角的減小逐漸減小。

2.3.2 下游泥沙淹沒范圍特征對(duì)比

下游泥沙淹沒范圍是影響攔擋壩選擇的重要因素之一，對(duì)于下游居民的生命安全有重要影響，因此有必要準(zhǔn)確掌握不同壩型對(duì)下游的淹沒范圍。圖11為數(shù)值模擬得到的各壩型的下游淹沒范圍，表5為模型試驗(yàn)與數(shù)值模擬得到的下游淹沒面積數(shù)據(jù)。綜合分析表明：數(shù)值模擬與模型試驗(yàn)得到的下游淹沒面積相近。

3 最優(yōu)防護(hù)措施確定

對(duì)比物理模型試驗(yàn)以及數(shù)值模擬的結(jié)果可知，對(duì)于該類型尾礦庫，壩體夾角為135°左右的攔擋壩可以起到較好的攔擋效果。135°左右的壩體起著攔擋和導(dǎo)流的雙重作用，壩體上的折角可減緩下泄泥沙對(duì)壩體的沖擊力，該壩型的雙向?qū)Я髯饔迷龃罅讼滦鼓嗌车牧鲃?dòng)范圍，降低了下游泥沙淤積厚度，能夠使得泥沙有效地排到周圍，減小了對(duì)下游居民的危害。

為進(jìn)一步確定最優(yōu)攔擋壩的壩體夾角，利用數(shù)值模擬方法對(duì)壩體夾角為126°、128°、130°以及136°的攔擋壩的防護(hù)效果進(jìn)行了分析。結(jié)果表明：由于該范圍內(nèi)壩體夾角變化幅度較小，下游淹沒范圍整體上變化不大；壩前淤積厚度峰值出現(xiàn)的時(shí)間以及大小也基本一致，如圖12所示。隨著下泄泥沙勢(shì)能減弱，壩前泥沙淤積厚度趨于穩(wěn)定，當(dāng)攔擋壩壩體夾角改變時(shí)，壩前泥沙淤積厚度發(fā)生了明顯變化。由圖12可知：在這4種工況下，壩前泥沙淤積厚度的最終穩(wěn)定值的變化特征為隨著壩體夾角減小而減小，原因在于壩體夾角越大，泥沙觸碰壩體時(shí)的勢(shì)能越大，泥沙在壩前發(fā)生回流的次數(shù)越多，使得壩前淤積厚度增大；攔擋壩的壩體夾角為126°時(shí)，壩前泥沙淤積厚度最終穩(wěn)定值的變化特征則恰恰相反，原因在于適當(dāng)減小壩體夾角，不僅不會(huì)改變攔擋壩的攔擋和雙向?qū)Я鞯碾p重作用，還會(huì)增強(qiáng)壩體的導(dǎo)流作用，從而增大了泥沙的流動(dòng)性；但攔擋壩的壩體夾角縮減過多，下泄泥沙運(yùn)動(dòng)距離增大，便會(huì)降低泥沙對(duì)壩體的沖擊力，減小壩前泥沙流動(dòng)性，導(dǎo)致壩體前泥沙淤泥過多，不利于壩體穩(wěn)定。綜上所述：壩體夾角為128°的攔擋壩為最優(yōu)防護(hù)措施。

4 結(jié) 論

針對(duì)下游為開闊平原地形的尾礦庫，通過模型試驗(yàn)與數(shù)值模擬兩種方法，研究了在不同下游攔擋壩布置形式下潰壩后各位置的泥沙淤積厚度以及淹沒范圍，并對(duì)下游防護(hù)攔擋壩的布置形式進(jìn)行了優(yōu)選，得出以下結(jié)論：

（1）攔擋壩具有攔擋和導(dǎo)流的作用，其中，導(dǎo)流作用分為橫向和縱向。當(dāng)泥沙觸碰到壩體時(shí)會(huì)產(chǎn)生兩種運(yùn)動(dòng)方式：①在壩體截流攔擋的作用下，大部分泥沙被截在壩前，小部分泥沙越過壩體繼續(xù)運(yùn)動(dòng)；②泥沙在攔擋壩的攔擋與導(dǎo)流雙重作用下，流體運(yùn)動(dòng)方向發(fā)生改變。

（2）不同壩型的攔擋壩所發(fā)揮的作用不同。對(duì)于壩體夾角較大的攔擋壩，壩體攔擋作用大于導(dǎo)流作用；對(duì)于壩體夾角為135°左右的攔擋壩，壩體的攔擋和導(dǎo)流作用均發(fā)揮了明顯作用；對(duì)于壩體夾角較小的攔擋壩，壩體的導(dǎo)流作用大于攔擋作用。

（3）壩體夾角為128°的攔擋壩為最優(yōu)防護(hù)措施。一方面在于夾角為128°的壩體起著攔擋和導(dǎo)流的雙重作用，壩體夾角能減緩下泄泥沙對(duì)壩體的沖擊力，同時(shí)此種壩型具有雙向?qū)Я髯饔?，增大了下泄泥沙的流?dòng)范圍，降低了下游泥沙淤積厚度，使得泥沙有效地排到周圍，減小了對(duì)下游居民的危害；另一方面，相對(duì)于135°的攔擋壩來說，適當(dāng)減小壩體夾角有利于減小泥沙對(duì)壩體的沖擊，減少壩前泥沙淤積，有利于壩體穩(wěn)定。