尾礦壩地震液化潰壩研究進(jìn)展

張文輝， 潘建平

(1.江西銅業(yè)集團(tuán)公司 德興銅礦， 江西 德興市 334224；2.江西理工大學(xué) 建筑與測(cè)繪工程學(xué)院， 江西 贛州市 341000)

0 引 言

隨著我國工業(yè)化對(duì)礦產(chǎn)資源需求的日益增加，開采出來的礦石經(jīng)破碎、磨礦、分選等多道工序才能選出精礦，在此過程中排出大量尾礦。我國尾礦年產(chǎn)達(dá)3億t以上，除少數(shù)被利用外，大部分以尾礦壩攔截儲(chǔ)存，現(xiàn)有尾礦庫約8540座，累積堆存尾礦超80億t[1￣2]。

尾礦壩是一種特殊的土工構(gòu)筑物，其筑壩材料、施工工藝、建筑工期、服務(wù)功能、靜動(dòng)力特性、壩體結(jié)構(gòu)形式、服務(wù)期和環(huán)境安全要求均不同于一般的水工土石壩。上游法尾礦壩后繼斷面較小及占地較少，且工藝簡單，造價(jià)低，在我國大多數(shù)尾礦壩中得到采用。然而，尾礦顆粒較細(xì)，比重較大，親水性弱，飽和疏松的尾礦是非常敏感的不穩(wěn)定結(jié)構(gòu)，而上游法尾礦壩容易形成復(fù)雜、混合的壩體結(jié)構(gòu)，許多壩體的排水系統(tǒng)不完善或易淤堵，致使壩體內(nèi)的浸潤線抬高甚至從壩面溢出，大部分壩體處于飽和狀態(tài)且比較疏松，容易引起壩體產(chǎn)生滲透破壞或滑坡、垮塌，在動(dòng)荷載條件下強(qiáng)度低，動(dòng)剪應(yīng)力比變化范圍小，地震時(shí)更易發(fā)生液化破壞[3￣4]。

在地震作用下，傾斜土體內(nèi)或其下部土層中如果發(fā)生液化，整個(gè)土體就可能流動(dòng)或向沒有支撐的一邊作側(cè)向移動(dòng)[5]，最終導(dǎo)致壩體發(fā)生流滑破壞。如1978年日本的Mochikoshi尾礦壩發(fā)生地震液化破壞，其中1號(hào)尾礦壩在地震結(jié)束時(shí)發(fā)生破壞，而2號(hào)尾礦壩則在震后約24 h垮塌，近80000 m3的尾礦水合物瞬間流向下游[6]。1985年智利Veta de Agua尾礦壩地震液化潰壩，向下游釋放了280000 m3尾礦，覆蓋長達(dá)5 km區(qū)域(見圖1)[7]。1994年加州地震中，Tapo Canyon尾礦壩發(fā)生了流滑破壞，也給當(dāng)?shù)貛砭薮蟮纳鷳B(tài)環(huán)境災(zāi)難[8]。我國1976年唐山地震中，大石河尾礦壩沉積灘發(fā)生大量噴水冒砂，且有明顯朝庫區(qū)滑移現(xiàn)象，在外坡第二級(jí)子壩上出現(xiàn)平行于壩軸線的裂縫，而在遠(yuǎn)離水邊線的沉積灘上也出現(xiàn)了眾多的微、小裂縫，后來實(shí)地勘察和計(jì)算分析都表明壩體確實(shí)曾有液化和局部滑移發(fā)生[9]。

試驗(yàn)研究已揭示了地震中砂土變形、孔壓上升、飽和砂土液化的特征，但是關(guān)于砂土液化的機(jī)理、強(qiáng)度降低以及液化砂土的性質(zhì)仍不是很清楚。另一方面，以位移為準(zhǔn)則來衡量巖土建筑物的抗震性能在近年來逐漸成為一種發(fā)展趨勢(shì)，這種方法考慮的是建筑物在地基液化中發(fā)生的位移量，用改善壩體的結(jié)構(gòu)和力學(xué)性能來進(jìn)行抗震設(shè)計(jì)[10]，從而提高壩體的抗震能力。

1 尾礦靜力和動(dòng)力特性

材料特性試驗(yàn)是研究尾礦壩抗震性能的重要環(huán)節(jié)之一，目前關(guān)于尾礦材料特性研究主要分為兩大類，一類是靜載試驗(yàn)，另一類是動(dòng)載試驗(yàn)。

Qiu等[11]通過對(duì)4種尾礦進(jìn)行試驗(yàn)，在半對(duì)數(shù)坐標(biāo)系內(nèi)，顯示有效應(yīng)力隨孔隙率的增大而降低，透水性隨孔隙率的增大而增大，近似線性關(guān)系；Fourie等[12]分析細(xì)粒含量對(duì)強(qiáng)度穩(wěn)態(tài)線位置高低的影響，表明較大細(xì)粒含量的尾礦強(qiáng)度穩(wěn)態(tài)線均高于較低含量者，意味著其需要較大的相對(duì)密度才不發(fā)生剪縮性，避免液化的發(fā)生。

阮元成等[13]對(duì)一種鐵尾礦的靜、動(dòng)強(qiáng)度和動(dòng)力變形特性進(jìn)行試驗(yàn)，結(jié)果表明尾礦料不僅在往返加載條件下動(dòng)強(qiáng)度低，動(dòng)剪應(yīng)力比變化范圍小，容易發(fā)生液化和破壞，而且，當(dāng)尾礦的密度小于某一臨界值時(shí)，靜力條件下也會(huì)發(fā)生流滑而進(jìn)入破壞狀態(tài)。張超等[14]根據(jù)銅尾礦動(dòng)力特性試驗(yàn)分析，指出當(dāng)細(xì)顆粒含量占到總量的35 %時(shí)其抗液化性能最佳。參考建筑抗震設(shè)計(jì)規(guī)范，并結(jié)合現(xiàn)場(chǎng)的標(biāo)準(zhǔn)貫入試驗(yàn)成果，提出適用于尾礦細(xì)粒含量對(duì)標(biāo)準(zhǔn)貫入擊數(shù)的修正式。陳存禮等[15]在不同固結(jié)狀態(tài)下對(duì)銅尾礦砂進(jìn)行動(dòng)三軸試驗(yàn)，認(rèn)為均壓固結(jié)和偏壓固結(jié)時(shí)孔壓增長曲線形態(tài)不同。均壓固結(jié)時(shí)，可用修正后的Seed應(yīng)力孔壓模型表達(dá)式來描述；偏壓固結(jié)時(shí)，可用指數(shù)函數(shù)來模擬。尹光志等[16￣17]利用尾礦細(xì)微觀力學(xué)與形變觀測(cè)試驗(yàn)裝置，進(jìn)行尾礦細(xì)觀觀測(cè)試驗(yàn)，結(jié)果表明受荷載作用，充水條件下各層尾礦顆粒沉降位移顯著大于未充水條件下相同位置顆粒的沉降位移。潘建平等[18]通過配置三組不同級(jí)配尾砂試樣，在不同圍壓下分別進(jìn)行不同孔隙比的固結(jié)不排水剪切試驗(yàn)，指出低圍壓時(shí)，各級(jí)配下飽和松散尾砂的應(yīng)力-應(yīng)變曲線呈軟化型，平均粒徑越小的尾砂，其峰值強(qiáng)度與穩(wěn)態(tài)強(qiáng)度越小，有必要加強(qiáng)細(xì)粒尾砂堆筑體液化勢(shì)的評(píng)估，以防液化潰壩災(zāi)難的發(fā)生。

2 尾礦壩模型試驗(yàn)

物理模型試驗(yàn)的基本思想是采用不同比尺，按照工作規(guī)律相似，對(duì)事物原型的現(xiàn)象或本質(zhì)進(jìn)行定性研究，再現(xiàn)或超前預(yù)測(cè)原型中已經(jīng)出現(xiàn)或可能出現(xiàn)的一些情況，以驗(yàn)證和解決實(shí)際工程問題。近年來隨著激振技術(shù)的不斷改進(jìn)和提高，振動(dòng)臺(tái)模型試驗(yàn)和離心機(jī)模型試驗(yàn)己逐步成為巖土工程研究領(lǐng)域的重要手段。

Zeng等[19]通過離心機(jī)模型試驗(yàn)，指出下游法修筑尾礦壩較上游法具有較高的抗震能力；盡管底部尾礦經(jīng)多年沉降固結(jié)，但細(xì)尾礦含水量仍高于液限值，說明尾礦具有較低強(qiáng)度。

陳建斌等[20]以某粉煤灰壩為研究對(duì)象，采用大型振動(dòng)臺(tái)作為試驗(yàn)手段，分別對(duì)不同高程灰壩的模型進(jìn)行了動(dòng)力特性試驗(yàn)研究。尹光志等[21]進(jìn)行細(xì)粒尾礦堆積壩加筋加固模型試驗(yàn)，通過研究獲知是否加筋的不同破壞模式。鄧濤等[22]通過室內(nèi)堆積壩模型試驗(yàn)，獲得壩體堆積過程中尾礦顆粒在庫內(nèi)沉積分布特征、浸潤線的變化規(guī)律。敬小非等[23]進(jìn)行尾礦壩在洪水情況下發(fā)生垮塌破壞的模型試驗(yàn)。結(jié)果表明壩體浸潤線變化存在滯后性；洪水導(dǎo)致壩坡尾礦砂所受的滲透力、上浮力、重力以及孔隙水壓力增大，削弱了壩體材料的黏聚力，并加大了自身荷載，降低了壩體的穩(wěn)定性。張興凱等[24]基于非恒定水流泥沙非平衡非飽和沖刷機(jī)理，根據(jù)模型相似理論和潰決侵蝕模型原理，模擬尾礦庫洪水漫頂潰壩過程，建立尾礦庫漫頂潰壩演化模型。劉磊等[25]以洛陽市欒川縣水露溝尾礦庫為工程背景，利用尾礦庫模型設(shè)計(jì)及試驗(yàn)方法建立了尾礦庫物理模型試驗(yàn)。在根據(jù)模型試驗(yàn)獲得的尾礦庫漫頂潰壩物理圖形基礎(chǔ)上，引入非平衡輸沙理論及潰口發(fā)展變化過程，建立了尾礦庫漫頂潰壩洪水預(yù)測(cè)數(shù)學(xué)模型，并利用試驗(yàn)資料對(duì)比計(jì)算模型結(jié)果。

3 尾礦壩地震液化和穩(wěn)定性分析

Seid等[26]用商用軟件模擬Mochikoshi尾礦壩的地震反應(yīng)，并對(duì)密實(shí)、反壓和穩(wěn)定柱三種抗震措施進(jìn)行簡單分析；Prodromos等[27]用有限元軟件對(duì)上游法、中線法和下游法三種壩型進(jìn)行動(dòng)力分析，指出場(chǎng)地條件、材料非線性和地震動(dòng)特征對(duì)尾礦壩動(dòng)力特性和內(nèi)部破壞都有影響。

柳厚祥等[28]采用有效應(yīng)力法，探討了高尾礦壩地震過程中和地震后孔隙水壓力的產(chǎn)生、擴(kuò)散和消散規(guī)律及其加速度、動(dòng)應(yīng)力和孔隙水壓力的響應(yīng)值。王飛躍等[29]將模糊可靠度理論應(yīng)用到尾礦壩的穩(wěn)定研究中，既考慮破壞事件的模糊性，又考慮主要變量和參數(shù)的模糊性，將模糊性和隨機(jī)性相結(jié)合，研究尾礦壩的模糊隨機(jī)可靠度分析方法?？讘椌┑萚30]用商用軟件分析尾礦壩地震液化流動(dòng)變形，發(fā)現(xiàn)壩體在不同條件下發(fā)生液化流動(dòng)變形的一些規(guī)律。結(jié)果表明增加灘長、降低水位線的高度，將有助于減小壩體液化流動(dòng)變形，增強(qiáng)壩體的穩(wěn)定性；壩基土發(fā)生液化流動(dòng)將導(dǎo)致尾礦壩整體發(fā)生流滑破壞。蔡嗣經(jīng)等[31￣32]對(duì)尾礦庫地震潰壩機(jī)理進(jìn)行了初步分析，建立評(píng)估尾礦庫抗震能力和地震潰壩可能性的指標(biāo)體系，并采用數(shù)值方法模擬得出尾礦壩易發(fā)生液化的區(qū)域主要在灘面水邊線部位、壩體中下部超孔隙水壓力較大部位、尾礦壩下游坡浸潤線溢出等部位。董隴軍等[33]針對(duì)尾礦壩地震穩(wěn)定的極限平衡分析問題，提出一種計(jì)算尾礦壩安全系數(shù)的新方法。結(jié)果表明運(yùn)用盲數(shù)運(yùn)算法可以計(jì)算穩(wěn)定性系數(shù)在不同取值區(qū)間內(nèi)的可信度，克服了傳統(tǒng)方法描述過于絕對(duì)化的問題。潘建平等以測(cè)試數(shù)據(jù)的統(tǒng)計(jì)分析結(jié)果為基礎(chǔ)，應(yīng)用一次二階矩法建立尾礦壩地震液化分析可靠度模型，可靠度理論在液化分析中能更好地考慮計(jì)算參數(shù)的變異性，不僅能判斷液化的發(fā)生與否，還能給出液化發(fā)生的概率。

4 尾礦壩地震液化潰壩研究方向

綜合上述研究現(xiàn)狀，國內(nèi)外針對(duì)尾礦壩的研究主要以靜力為主，動(dòng)力方面的研究相對(duì)較少，尾礦液化后力學(xué)特性、尾礦壩地震液化細(xì)觀破壞機(jī)制及流滑運(yùn)移特征方面的研究則更少。因此，探明尾礦液化過程中的細(xì)微觀結(jié)構(gòu)變形特征，對(duì)于研究尾礦壩的地震液化破壞機(jī)制具有重要的實(shí)際意義。

隨著剪應(yīng)變的發(fā)展，液化后尾礦呈現(xiàn)出流態(tài)化運(yùn)動(dòng)特征，其運(yùn)動(dòng)、變形的性質(zhì)與黏性流體相似，不能再視其為固態(tài)。為此，基于流體力學(xué)理論，從流動(dòng)破壞角度分析尾礦液化后的變形特性是很有必要的。針對(duì)這方面的問題開展研究，可以為尾礦壩抗震設(shè)計(jì)和防災(zāi)減災(zāi)決策提供理論依據(jù)。

參考文獻(xiàn)：

[1]趙曉晨.尾礦庫：人造“泥石流”災(zāi)害源[J].百科知識(shí), 2009.

[2]謝旭陽, 王云海, 張興凱, 等. 我國尾礦庫數(shù)據(jù)庫的建立[J].中國安全生產(chǎn)科學(xué)技術(shù),2008, 4(2):53￣56.

[3]阮元成, 郭 新.飽和尾礦料靜、動(dòng)強(qiáng)度特性的試驗(yàn)研究[J].水利學(xué)報(bào),2004(1):67￣73.

[4]魏 勇,許開立,鄭 欣.尾礦壩潰壩事故致因分析[J].中國礦業(yè),2009,18 (6):98 ￣99.

[5]陳文化.液化流滑變形模型及初步實(shí)驗(yàn)[J].自然災(zāi)害學(xué)報(bào),2004,13(3):75￣80. [6]Ishihara K. Post-earthquake failure of a tailings dam due to liquefaction of the pond deposit[C]//International conference on case histories in geotechnical engineering,Stolouis,Geotechnical Engineering, 1984(3):1129￣1143.

[7]Chronology of major tailings dam failures. http://www.wise-uranium.org/mdaf.html.

[8]Harder J L F, Stewart J P. Failure of Tapo Canyon tailings dam[J].Journal of Performance of Constructed Facilities,1996,10(3):109￣114.

[9]王余慶,王治平,等.大石河尾礦壩1976年唐山大地震震害及有關(guān)強(qiáng)震觀測(cè)記錄[J].工業(yè)建筑,1994,24(7): 38￣42. [10] Gallage C P K, Towhata l, Nishimura S. Laboratory investigation on rate-dependent properties of sand underloading low confining effective stress[J]. Soils and Foundations,2005,45(4):43￣60. [11] Qiu Y X, Sego D C. Laboratory properties of mine tailings[J].Canadian Geotechnical Journal,2001,38(1):183￣190.

[12] Fourie A B,Papageorgiou G.Defining an appropriate steady state line for Merriespruit gold tailings[J].CanadianGeotechnicalJournal,2001,38(4):695￣706.

[13]阮元成,郭 新.飽和尾礦料動(dòng)力變形特性的試驗(yàn)研究[J].水利學(xué)報(bào),2003,34(4):24￣29.

[14]張 超,楊春和.細(xì)粒含量對(duì)尾礦材料液化特性的影響[J].巖土力學(xué),2006,27(7):1133￣1137.

[15]陳存禮,何軍芳,胡再強(qiáng),等.動(dòng)荷作用下飽和尾礦砂的孔壓和殘余應(yīng)變演化特性[J].巖石力學(xué)與工程學(xué)報(bào),2006,25(z2):4034￣4039.

[16]尹光志,等.尾礦細(xì)微觀力學(xué)與變形觀測(cè)試驗(yàn)裝置的研制與應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào),2011,30(5):926￣934.

[17]尹光志,張千貴,魏作安,等.孔隙水運(yùn)移特性及對(duì)尾礦細(xì)觀結(jié)構(gòu)作用機(jī)制試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào)，2012,31(1):71￣79.

[18]潘建平,劉湘平,王宇鴿.不同顆粒級(jí)配尾砂穩(wěn)態(tài)強(qiáng)度特性試驗(yàn)研究[J].中國安全生產(chǎn)科學(xué)技術(shù),2014,10(6):39￣44.

[19] Zeng X W, Wu J,Rohlf R. Laboratory properties of mine tailings[J].Geotechnical News,1998,6:29￣32.

[20]陳建斌,周立運(yùn).粉煤灰壩動(dòng)力特性試驗(yàn)研究[J].巖土力學(xué),2005,26(3):437￣440.

[21]尹光志,魏作安,等.細(xì)粒尾礦堆壩加筋加固模型試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào)2005,24(6):1030￣1034.

[22]鄧 濤,萬 玲,魏作安.溫莊尾礦庫堆壩模型試驗(yàn)及壩體穩(wěn)定性分析[J].巖土力學(xué),2011,32(12):3647￣3652.

[23]敬小非,尹光志,魏作安,等.尾礦壩垮塌機(jī)制與潰壩模式試驗(yàn)研究[J].巖土力學(xué),2011,32(5):1377￣1384.

[24]張興凱,孫恩吉,李仲學(xué).尾礦庫洪水漫頂潰壩演化規(guī)律試驗(yàn)研究[J].中國安全科學(xué)學(xué)報(bào),2011,21(7):118￣124.

[25]劉 磊,張紅武,鐘德鈺,等.尾礦庫漫頂潰壩模型研究[J].水利學(xué)報(bào),2014,45(6):675￣681.

[26]Seid K M,Byrne P M. Embankment dams and earthquakes[J].Hydropower & Dams,2004,2:96￣102.

[27] Prodromos N P,Yiannis T.Stability of tailings dams under static and seismic loading[J]. Canadian Geotechnical Journal,2008,45(5):663￣675.

[28]柳厚祥,廖 雪,李 寧,等.高尾礦壩的有效應(yīng)力地震反應(yīng)分析.振動(dòng)與沖擊,2008,27(1):65￣70.

[29]王飛躍,等.尾礦壩穩(wěn)定性分析的模糊隨機(jī)可靠度模型及應(yīng)用[J].巖土工程學(xué)報(bào),2008,30(11):1600￣1605.

[30]孔憲京,潘建平,鄒德高.尾礦壩液化流動(dòng)變形分析[J].大連理工大學(xué)學(xué)報(bào),2008,48(4):0541￣0545.

[31]張 棟,吳宗之,等.尾礦庫地震液化機(jī)理及抗震評(píng)估指標(biāo)體系研究[J].中國安全生產(chǎn)科學(xué)技術(shù),2010,6(6):17￣22.

[32]蔡嗣經(jīng),張 棟,何 理.地震中尾礦庫液化失穩(wěn)機(jī)理及數(shù)值模擬研究[J].有色金屬科學(xué)與工程,2011,2(2):1￣6.

[33]董隴軍,等.尾礦壩地震穩(wěn)定性分析的區(qū)間模型及應(yīng)用[J].中南大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,42(1):164￣169.