劉曉峰
(中國瑞林工程技術(shù)股份有限公司,江西南昌 330038)
尾礦庫是礦山企業(yè)生產(chǎn)中用于貯存選礦后尾礦的重要設(shè)施,也是潛伏巨大安全隱患的危險(xiǎn)源。尾礦庫的安全穩(wěn)定不僅關(guān)系到礦山的正常生產(chǎn)運(yùn)行,還關(guān)系到下游人民群眾的生命及財(cái)產(chǎn)安全。 國內(nèi)外尾礦庫潰壩事故累見不鮮[1]:2008 年山西襄汾新塔礦業(yè)“9·8”潰壩事故造成 277 人死亡;2010 年廣東信宜紫金礦業(yè)“9·21” 潰壩事故造成 22 人死亡;2014年陜西商洛場“10·4”潰壩事故造成 5 人死亡;2019 年巴西淡水河公司“1·25”潰壩事故造成235 人死亡。
導(dǎo)致尾礦庫潰壩事故發(fā)生的原因較多[2],且通常是由多種致災(zāi)因素疊加后共同作用的結(jié)果。 由于尾礦庫潰壩事故影響較大, 國家對尾礦庫的運(yùn)行安全越來越重視,因此特別要求對尾礦庫,尤其是“頭頂庫”進(jìn)行尾礦庫潰壩風(fēng)險(xiǎn)評估,研究制定治理和預(yù)防事故的措施,科學(xué)制定居民搬遷計(jì)劃。尾礦庫潰壩模擬能為尾礦庫應(yīng)急預(yù)案提供依據(jù), 指導(dǎo)確定安全緊急疏散路線, 對于尾礦庫災(zāi)害預(yù)防及下游居民生命財(cái)產(chǎn)安全保障有著重要的意義[3]。
尾礦庫潰壩的致災(zāi)機(jī)理[4]主要包括滲透破壞、地震液化、壩體滑坡及洪水漫頂?shù)取?1)滲透破壞。 滲透破壞主要是由于尾礦庫庫內(nèi)水位的上升,導(dǎo)致尾礦庫干灘長度減小,尾礦壩壩內(nèi)浸潤線升高,最終誘發(fā)滲流破壞。2)地震液化。地震液化主要是由于堆積壩尾砂粒度較細(xì)或尾礦欠密實(shí),在地震作用下,超靜孔隙水壓力迅速增加,導(dǎo)致尾砂液化,進(jìn)而形成裂隙滲流通道,造成壩體失穩(wěn)破壞等。 3)壩體滑坡。 壩體滑坡主要因壩體坡度較陡、 筑壩材料黏聚力及抗剪強(qiáng)度較低,進(jìn)而引起壩體局部或整體失穩(wěn)破壞。4)洪水漫頂。 洪水漫頂主要是因突發(fā)超標(biāo)準(zhǔn)暴雨、洪水,或排洪設(shè)施失效等,導(dǎo)致尾礦庫庫內(nèi)水位持續(xù)上升,致使子壩擋水,且洪水最終漫過子壩,進(jìn)一步誘發(fā)了壩面拉溝、壩體失穩(wěn)坍塌及尾礦水外泄等危害。
滲透破壞、地震液化及壩體滑坡這3 類致災(zāi)機(jī)理傾向于尾礦壩穩(wěn)定分析范疇,在設(shè)計(jì)過程中可通過降滲、抗液化及壩坡放緩等措施控制。而洪水漫頂則是由于出現(xiàn)了突發(fā)暴雨,且排洪系統(tǒng)堵塞失效而造成尾礦庫事故。相對前3 類潰壩方式,洪水漫頂導(dǎo)致的尾礦庫潰壩事故對下游影響范圍更廣,危害更嚴(yán)重。
尾礦壩漫頂潰決過程及其機(jī)理[5]可以大致描述為:洪水漫過尾礦壩,漫頂水流在重力作用下不斷沖刷尾礦壩。在尾礦庫下部位置,漫頂水流率先達(dá)到尾砂起動流速, 尾礦壩壩身的尾砂開始與漫頂水流一起形成挾砂的尾礦水進(jìn)一步?jīng)_刷壩體下游。 尾礦壩壩體經(jīng)過不斷地沖刷和掏深,形成了坡度很陡的深溝。 當(dāng)深溝兩側(cè)尾砂體不能維持重力作用下的穩(wěn)定時(shí),壩體發(fā)生垮塌,使得潰口不斷變寬變深,并逐漸向上游延伸,致使更多的尾礦水流出。尾礦水流速在崎嶇地形和沿程阻力的作用下逐漸降低, 尾砂的水挾砂能力也逐步降低,最終尾砂沿程沉積。
目前對于尾礦庫潰壩模擬主要有物理模型試驗(yàn)[6]和數(shù)值模擬[7]兩種方法:1)物理模型實(shí)驗(yàn)。 物理模型試驗(yàn)是將地形地貌等按比尺縮小制作, 根據(jù)水平向及豎向比尺縮小比例是否相同, 可分為正態(tài)模型和變態(tài)模型,然后再根據(jù)相似原理,選擇與實(shí)際尾砂相似條件較好的模型砂, 對潰壩過程及影響進(jìn)行實(shí)時(shí)展現(xiàn)。2)數(shù)值模擬。數(shù)值模擬則主要通過電子計(jì)算機(jī)制作好尾礦庫庫區(qū)以及其下游周邊地形的模型,再結(jié)合數(shù)值計(jì)算軟件, 選擇好數(shù)學(xué)模型, 設(shè)置模型參數(shù)、邊界條件、初設(shè)條件及監(jiān)測點(diǎn),最終通過圖像顯示潰壩過程及潰壩后對下游影響范圍及程度。 相比而言,物理模型試驗(yàn)搭建耗時(shí)較長,重復(fù)性較差,造價(jià)高,而數(shù)值模擬簡便靈活,便于邊界條件的調(diào)整,重復(fù)性較好,相對經(jīng)濟(jì)。
本文以某“頭頂庫”為例,采用CFD 軟件數(shù)值模擬分析尾礦庫漫頂潰壩后的淹沒范圍及影響程度。
根據(jù)潰壩過程的時(shí)間長短,可分為瞬間潰壩與逐漸潰壩。 由于尾礦壩壩體破壞不像剛性壩體那樣瞬間潰決,且其潰口是逐步發(fā)展擴(kuò)大的,故尾礦壩漫頂潰壩一般為逐漸潰壩??紤]到尾礦壩與水庫中的土壩相對較為相似,且尾礦庫漫頂潰壩的水流與水庫潰決也頗為相似,因此可用工程類比法,借用水利學(xué)中的潰壩公式[8]。為反應(yīng)尾礦庫逐漸潰壩過程,在制作尾礦庫模型時(shí),事先在尾礦壩上設(shè)置一缺口,其寬度根據(jù)黃河水利委員提出的堤壩潰口寬度經(jīng)驗(yàn)公式設(shè)置:
式中:b 為潰口寬度,m;W 為水庫總庫容,m3;B 為主壩長度,m;H 為壩高,m;k 為壩體材料經(jīng)驗(yàn)系數(shù)(黏土類取0.65,壤土取1.30)。
為體現(xiàn)尾礦庫潰壩后最不利情況, 初始潰口位置設(shè)置于尾礦壩最大斷面處。
潰壩數(shù)值模擬應(yīng)用成果見圖1~圖8。
圖1 尾礦庫及下游構(gòu)筑物模型
圖2 下游村莊監(jiān)測點(diǎn)
圖3 潰壩10 min 后流速云圖
圖4 潰壩20 min 后流速云圖
圖5 潰壩1 h 后流速云圖
圖6 潰壩2 h 后流速云圖
圖7 尾砂淤積厚度云圖
圖8 下游村莊尾砂厚度監(jiān)測值
通過潰壩模擬結(jié)果可知:1)模擬實(shí)驗(yàn)開始后,潰壩水流在勢能作用下沿預(yù)設(shè)潰口發(fā)生漫頂, 漫頂水流沖刷壩體下游側(cè)壁面,形成沖溝。 在水流漫頂10 min 后,沖溝正式形成,此時(shí)沖溝內(nèi)水流流速達(dá)8 m/s。2)隨著水流繼續(xù)沖刷壩體,沖溝迅速發(fā)育,發(fā)展成潰口, 下泄流量劇增, 流速增大, 在水流漫頂20 min后,最大水流流速達(dá)9 m/s,前端水舌已達(dá)谷口開闊地。 此處地勢開闊,水流速度下降,并造成淤積。 同時(shí),庫內(nèi)水面下降,潰口繼續(xù)擴(kuò)寬,下泄流量增大,導(dǎo)致下游窄谷中水位上升。3)漫頂潰壩發(fā)生1 h 后,水流開始對壩體下游村莊建筑造成影響。 此時(shí)潰口已經(jīng)擴(kuò)展到最大,庫內(nèi)水面大幅度降低,下泄流量相應(yīng)減小,尾砂沿程沉積于壩體下游溝谷中。4)庫內(nèi)水位逐步干涸,潰口停止發(fā)展,下泄流量持續(xù)減小,直至斷流。2 h 后,潰口下游基本無徑流,只有少部分低洼地帶還有少量水流殘存。5)在潰壩發(fā)生后,尾砂隨水流下泄,其中大部分尾砂都沉積于下游的溝谷之中,僅小部分尾砂隨水流移動至村莊低洼區(qū)沉積。 最終村莊內(nèi)僅幾棟建基面較低的建筑會受尾砂沉積影響(僅 probe3 和 probe4 監(jiān)測點(diǎn)附近受影響,probe1 和probe2 監(jiān)測點(diǎn)區(qū)未受影響),最大深度約1.8 m。
1)通過尾礦庫潰壩數(shù)值模擬,可預(yù)測潰壩后對下游的影響范圍及程度。 據(jù)此可為尾礦庫失事提前預(yù)警、人員疏散及相關(guān)安全預(yù)防措施等提供重要的依據(jù)。
2)利用流體力學(xué)軟件對尾礦庫潰壩進(jìn)行數(shù)值模擬是較新的一個(gè)研究方向,為使計(jì)算結(jié)果更接近真實(shí)工況,應(yīng)結(jié)合實(shí)際潰壩案例,給模型賦予準(zhǔn)確的參數(shù)。
3)為盡量避免或減輕潰壩造成的危害,應(yīng)嚴(yán)格按照相關(guān)要求,加強(qiáng)尾礦庫生產(chǎn)運(yùn)行管理;同時(shí)可引導(dǎo)找出萬一潰壩發(fā)生其潰口的位置, 即在筑壩過程中, 靠厚實(shí)山體一側(cè)較其它段的尾礦壩灘頂高程可稍低些。