石滿生

（江西天億礦業(yè)有限公司，江西 貴溪 335424）

1 引言

傾 大多數(shù)尾礦壩依山而建，下游不乏人口密集區(qū)、重要建構(gòu)筑物、交通要道和河流等，一旦潰壩，給下游人民生命和財產(chǎn)安全帶來嚴(yán)重危害，且污染下游生態(tài)環(huán)境［1,2]。因此，尾礦庫潰壩災(zāi)害成為礦山防災(zāi)減災(zāi)的重要難題，深入研究尾砂流滑運移特征及對下游區(qū)域建構(gòu)物的影響范圍，具有重要的現(xiàn)實意義［3,4]。李彥娥等［5]采用單溝泥石流經(jīng)驗預(yù)測公式對尾礦庫潰壩后泥石流覆蓋范圍及影響程度進行定量分析評價?；矢P龍等［6]應(yīng)用二維潰壩水流數(shù)值模型和尾礦壩潰決砂流數(shù)值模型，分別模擬某尾礦壩潰決洪水、砂流運動狀況。劉磊等［7]對尾礦壩潰口及壩體崩塌變化過程、下泄洪水演化過程等進行觀測。郭天勇等［8]建立潰壩滑坡體滑移距離的計算公式，分析擴散角、空氣阻力、殘余強度、撞擊動壓力和攔擋壩等對滑移距離的影響。

尾礦庫潰壩作為一種災(zāi)害事故具有瞬間爆發(fā)性和很強的破壞性，要想得到真實可靠的潰壩現(xiàn)場監(jiān)測數(shù)據(jù)是非常困難的，并且需要巨大的人力和財力的投入［9]。本文針對尾礦庫潰壩事故特點，把由尾砂和水組成的固液混合物視為均質(zhì)、連續(xù)、不可壓縮的非牛頓流體，建立潰壩尾砂流滑數(shù)值模型，采用FLUENT軟件模擬潰壩尾砂的流滑運移過程，探究潰決尾砂流滑影響范圍，為尾砂潰壩災(zāi)害的預(yù)報與災(zāi)害的防治提供理論支撐與技術(shù)指導(dǎo)。

2 工程簡況

該尾礦庫位于選廠北側(cè)長約1600m，寬約400m的條帶狀溝谷中，植被茂密，清泉長流。最終堆積壩頂高程為310m，初期壩軸線底部最低高程為235m，計算尾礦壩最大壩高Hmax=75m，尾礦庫的總庫容為247.75×104m3。根據(jù)《尾礦設(shè)施設(shè)計規(guī)范》（GB 50863-2013）［10]，該尾礦庫為三等庫，庫內(nèi)主要水工構(gòu)筑物為3級，次要構(gòu)筑物為5級，臨時構(gòu)筑物為5級。具體地形分布見圖1。尾礦庫下游有公路、河流，防護山脊下有居民區(qū)。

圖1 尾礦庫地形

3 FLUENT計算原理

尾礦庫潰決后，庫內(nèi)尾水夾帶尾礦形成泥漿一起沖向庫區(qū)下游，對下游建構(gòu)物造成巨大沖擊。把潰壩所形成的泥漿體看作是空氣和泥砂的兩相流體，可用歐拉-歐拉模型來描述其運動。在歐拉-歐拉模型中，不同相被認為互相貫穿的連續(xù)介質(zhì)。體積率是空間和時間的連續(xù)函數(shù)，各相的體積率之和等于1.0。從各相的守恒方程可以推導(dǎo)出一組方程，這些方程對于所有的相都具有類似的形式［11]。從實驗得到的數(shù)據(jù)可以建立一些特定的關(guān)系，從而能使上述方程封閉，另外，對于小顆粒流（Granular Flows），則可以通過應(yīng)用分子運動論的理論使方程封閉。在歐拉-歐拉模型中，常見的模型有三種，分別為：流體體積模型（VOF）、混合物模型及歐拉模型。

經(jīng)過驗證，選用流體體積模型（即VOF模型）對尾礦壩潰壩水砂流的運動過程進行數(shù)值模擬比較合適［9]。水氣界面的追蹤可通過求解相體積分?jǐn)?shù)F的運輸方程完成

4 潰壩數(shù)值模擬

4.1 計算模型

根據(jù)模擬現(xiàn)場壩體幾何尺寸，確定尾礦庫潰壩的計算模型。為了方便模型的建立、降低硬件資源消耗和縮短計算時程，將實際情況進行了適當(dāng)簡化，如庫內(nèi)溝谷上大下小斷面，折算成一個寬度的斷面，下游坡度和落差按實際地形確定，所建模型重點考慮潰壩水砂流對下游左側(cè)民宅的影響。壩體外坡為1∶4，總壩高為75m，庫容量為250萬m3，模型如圖2所示。

圖2 計算模型

假設(shè)尾礦庫潰壩是在最不利的情況下發(fā)生的，也就是瞬間全潰并且潰決水砂量最大的情況。經(jīng)估算并結(jié)合現(xiàn)場情況，綜合確定模型壩體部分的尺寸，本文分析以壩高75m情況下潰決水砂流的流動過程模擬，一種是未采取防范措施，另一種是在下游228m平臺處修3m高混凝土擋土墻，通過模擬觀察潰壩水砂流的淹沒過程和壓力變化情況。

根據(jù)類似工程有限元數(shù)值分析的經(jīng)驗，結(jié)合現(xiàn)場實際情況，確定邊界條件。由于該尾礦庫建在山谷中，則計算區(qū)域的下游斷面為開邊界，上游、兩岸和地面為閉邊界。

4.2 材料特性定義

在FLUENT中對多相流的設(shè)置中需要一些材料參數(shù)，這些參數(shù)是進行模擬的基礎(chǔ)和成敗的關(guān)鍵。根據(jù)現(xiàn)場取樣尾礦的物理力學(xué)性質(zhì)測試的結(jié)果，以礦漿濃度為40%，下游溝谷按實際坡度和落差，潰決形式為瞬間全潰為基礎(chǔ)，參照國內(nèi)外有關(guān)尾礦漿和泥石流等相關(guān)資料，綜合分析后，確定出各計算參數(shù)，見表1。

表1 模型計算參數(shù)

4.3 結(jié)果分析

建立幾何模型后，設(shè)立邊界條件，導(dǎo)入到FLUENT程序中，檢查網(wǎng)格質(zhì)量，輸入計算參數(shù)，進行計算。

從圖3可知，10～20s時間段下游淹沒區(qū)增長最快，這主要是初期壩下游有一個7m的臺階落差，使得此段水砂流速度很大。由于下游左側(cè)突出防護小山體的阻擋作用，使得洪峰最高中心往遠離民宅方向發(fā)生了偏移，這大大減小了下游民宅的防洪壓力。在36s時刻，淹沒民宅區(qū)范圍最大，漫流水砂遇到民宅阻擋，水位抬到5m多高，對房子形成一定的沖擊壓力，整個影響時間大約在15s左右。在44～46s時間段，下游河邊的淹沒區(qū)最大，然后隨著上游水位的下降和下游的泄流，淹沒區(qū)緩慢減少。

圖3 不同時刻潰壩泥砂流的運行情況

針對上述模擬結(jié)果，經(jīng)過實地調(diào)研并與業(yè)主充分討論，本著安全、經(jīng)濟與實用的原則，擬采取在228m平臺處靠近民宅一側(cè)修建一3m高混凝土擋墻，以預(yù)防水砂流漫到居民區(qū)，起到保護民宅的作用。對比圖3（d）與圖4可知，修建擋土墻后，在同一時刻（36s），潰壩尾砂對民宅的影響大大減小。

從圖5可知，在潰壩泥砂作用下，增設(shè)擋土墻的總壓應(yīng)力先增大后減小，在20s時最大，此時主要是上游尾礦從臺階上沖下，對下游平臺有較大的沖擊力所致，所以擋土墻設(shè)置處的防護措施應(yīng)嚴(yán)格到位。隨著下游水砂的下泄，總壓應(yīng)力總體上越來越小。

圖4 36s時潰壩泥砂流的淹沒情況（修筑擋墻后）

圖5 擋土墻總壓應(yīng)力隨時間的變化

5 結(jié) 論

（1）為了減小潰壩發(fā)生的概率，汛期應(yīng)加強對排洪系統(tǒng)、壩體的巡查，確保排洪系統(tǒng)完好暢通及壩體完好。

（2）事先對壩體潰壩數(shù)值進行研究分析，做好應(yīng)急預(yù)案及防護措施，確保下游人民生命財產(chǎn)安全、減少對下游流域生態(tài)環(huán)境帶來的影響。

（3）采用FLUENT軟件模擬尾礦庫潰壩后泥砂流滑運移特征及影響范圍，進一步豐富了尾礦庫潰壩災(zāi)害風(fēng)險防范的手段，為業(yè)主在尾礦庫的安全管理及風(fēng)險防范提供參考。

（4）本案由于下游坡度較陡，河谷狹窄，未采取任何措施時，尾礦庫潰壩水砂流會漫到居民區(qū)，對下游左側(cè)民宅存在一定影響。采取修建1處3m高的混凝土擋土墻，預(yù)防尾砂流對民宅的影響，從分析結(jié)果來看，取得較好防護效果。