歐孝奪 蘇 建 鐘一和 吳光航 江 杰

（1.廣西大學(xué)土木建筑工程學(xué)院，廣西南寧530004；2.工程防災(zāi)與結(jié)構(gòu)安全重點(diǎn)實(shí)驗(yàn)室，廣西南寧530004；3.廣西金屬尾礦安全防控工程技術(shù)研究中心，廣西南寧530004）

礦石選洗所剩尾礦泥漿經(jīng)泵送排入尾礦庫中堆儲(chǔ)，具有含水率高、細(xì)粒含量多等特點(diǎn)［1］。由于缺乏有效的排水通道，尾礦漿經(jīng)十幾年時(shí)間仍未能完全固結(jié)，尾礦庫內(nèi)積水嚴(yán)重，極易引發(fā)漫頂、潰壩及庫底泥漿泄漏等事故。通過增設(shè)“逐級(jí)濕堆尾礦排水裝置”［2］為排水體，可有效加快尾礦泥漿橫縱向排水固結(jié)。該排水裝置以排水性能良好且能有效過濾細(xì)顆粒尾礦的土工排水材料作為豎向排水通道。排水材料下部固定連接至庫底排水管，上部蜷曲于可隨尾礦堆排同步升高的浮體中。隨著尾礦泥漿的逐級(jí)堆排，排水材料自浮體中隨之拉伸，過濕性尾礦上部積水及孔隙水隨著這一排水通道持續(xù)性向下排水，最終從庫底排水管集中排出。

尾礦的排水固結(jié)與尾礦類型、顆粒粒徑、沉積時(shí)間、預(yù)壓力施加時(shí)間等因素密切相關(guān)［3-5］，沉降量與沉降時(shí)間基本呈線性正相關(guān)關(guān)系，且濃度越高，相關(guān)系數(shù)越?。?］。尾礦泥漿的逐級(jí)排放工藝使尾礦土在飽和—非飽和狀態(tài)中循環(huán)［7］，并在沉積層中形成交錯(cuò)層理［8］。巫尚蔚［9］將尾礦的一維自重沉積過程分為沉降階段和固結(jié)階段。但在增設(shè)排水體后，尾礦排水固結(jié)過程及滲流計(jì)算模型更為復(fù)雜，既有研究難以有效描述其在逐級(jí)堆排條件下排水固結(jié)過程的特性變化規(guī)律。

通過等比例模型試驗(yàn)?zāi)M逐級(jí)堆排沉積條件下尾礦泥漿的排水過程，對(duì)該條件下尾礦橫縱向固結(jié)效果進(jìn)行了分析，并將尾礦排水固結(jié)分為5個(gè)階段；進(jìn)一步建立尾礦理想孔隙滲流模型，計(jì)算得到排水體最大影響半徑時(shí)序曲線，為排水體作用下逐級(jí)堆排尾礦的排水固結(jié)提供參考。

1 排水固結(jié)模型試驗(yàn)

1.1 尾礦排水固結(jié)物理模型

從單次堆排過程來看，尾礦泥漿通過泵送排放至尾礦庫內(nèi)，隨后懸浮泥漿通過豎向排水體排出自由水及孔隙水，土顆粒在自重作用下逐步固結(jié)，最終轉(zhuǎn)變?yōu)榫邆浣Y(jié)構(gòu)強(qiáng)度的沉積土。而由于尾礦庫采用間隔式的逐級(jí)堆排，故單次泥漿排水固結(jié)至一定程度后，下一次泥漿又會(huì)重新補(bǔ)充水分，使其固結(jié)狀態(tài)發(fā)生改變。根據(jù)逐尾礦逐級(jí)堆排特點(diǎn)，建立其排水固結(jié)物理模型如圖1所示。（1）分級(jí)填入既定厚度尾礦泥漿，排水體持續(xù)作業(yè)；（2）泥漿僅在自重作用下完成排水固結(jié)，故排水口在排水體下方；

（3）表層蒸發(fā)面及下端重力排水口為排水面，圓柱其余面為不排水面。

1.2 模型試驗(yàn)設(shè)計(jì)

（1）模型設(shè)計(jì)。根據(jù)上述物理模型設(shè)計(jì)制作等比例試驗(yàn)裝置（如圖2所示）。試驗(yàn)箱為高1.8 m，直徑2.0 m的不銹鋼圓筒，底部中心設(shè)置砂濾層排水口。排水口連接土工復(fù)合排水材料為豎向排水體，排水材料由表面濾網(wǎng)及骨架芯材2部分組成，上端固定在安裝支架中，使排水體保持豎直。

（2）尾礦泥漿制備及堆排。所選尾礦為黏粒含量較大的鋁土尾礦，礦泥取自中國(guó)鋁業(yè)廣西分公司2#尾礦庫，經(jīng)室內(nèi)試驗(yàn)測(cè)試其基本性質(zhì)如表1所示。

?

尾礦土樣經(jīng)過“曬干—初步粉碎—充分浸泡—攪拌”過程后制成20%濃度分散均勻的泥漿。堆排速率和每級(jí)堆排高度根據(jù)近2 a中國(guó)鋁業(yè)廣西分公司2#尾礦庫實(shí)際情況，日排泥高度為0.07 m/d，堆排21 d后試驗(yàn)箱達(dá)到飽和。最終泥面高度達(dá)到1.475 m，水面高度達(dá)到1.747 m（如圖3所示）。試驗(yàn)共持續(xù)125 d，其中前期堆排21 d，停排后繼續(xù)測(cè)試104 d。

2 雙向排水固結(jié)分析

2.1 孔隙水壓力變化分析

試驗(yàn)布設(shè)9個(gè)孔隙水壓力計(jì)（1#～9#），各位置孔隙水壓力變化時(shí)程曲線如圖3所示?？紫端畨毫η€存在一定波動(dòng)，但總體上呈現(xiàn)明顯階段性。

（1）尾礦堆積高度每增加0.5 m，孔隙水壓力平均減小約5 kPa；同一高度的各位置孔隙水壓力隨著與排水通道的距離的減小而減小，但減幅逐漸增大。

（2）堆排期孔隙水壓力呈波浪式逐級(jí)遞增，停排期孔隙水壓力呈兩段式逐漸減小，在約第35 d后下降速度明顯放緩?？紫端畨毫ψ兓瘯r(shí)間拐點(diǎn)與上層覆水排干時(shí)間吻合，故認(rèn)為停排前期（第21～35 d）尾礦主要是在外加水頭下產(chǎn)生滲流排出自由水，而停排后期（第35 d后）尾礦主要為固結(jié)過程中產(chǎn)生的滲流，尾礦壓縮后所含的孔隙水被排出，滲透性降低，孔隙水壓力值減幅變小。

2.2 沉降變化分析

試驗(yàn)共布設(shè)沉降標(biāo)9個(gè)（1#～9#），從堆排8 d后（堆填高度0.5 m布設(shè)第一層）開始測(cè)量，各位置沉降時(shí)程變化曲線如圖4所示。

從圖4來看，雙向沉降變化與孔隙水壓力變化并不同步，土體應(yīng)力變化滯后于應(yīng)變變化。

（1）對(duì)于不同高度的沉降測(cè)試點(diǎn)，下層土的沉降很快達(dá)到穩(wěn)定，且單層變形量更大。最終底層（1#～3#）平均單層沉降變形152.6 mm，遠(yuǎn)大于上層（7#～9#）的平均單層沉降變形22.7 mm；對(duì)同一高度的各沉降值點(diǎn)，隨著滲流路徑增加沉降量越來越接近，且該規(guī)律隨深度增加越發(fā)明顯，底層距離排水材料0.66 m與0.98 m位置沉降量幾乎相等。

（2）底層及中間層沉降測(cè)試點(diǎn)（1#～6#）的沉降曲線可分為2個(gè)階段：堆排期（21 d前）逐級(jí)堆載與排水同時(shí)進(jìn)行，沉降降幅較大；停排期（21 d后）僅通過排水固結(jié)壓縮尾礦體，沉降明顯放緩；尾礦表層沉降測(cè)試點(diǎn)（7#～9#）自停排后開始安裝，沉降拐點(diǎn)與底層及中間層不同步，除向下排水外其表面還受到蒸發(fā)作用長(zhǎng)期影響，試驗(yàn)后期沉降增幅明顯高于底層及中層。

2.3 逐級(jí)堆排尾礦排水過程分析

尾礦在豎向排水體作用下，由完全擾動(dòng)的水土混合懸浮泥漿轉(zhuǎn)化為土水分離的結(jié)構(gòu)性沉積土，該沉降及自重固結(jié)過程亦是泥土系統(tǒng)（懸浮態(tài)泥和沉積態(tài)土共存的系統(tǒng)）內(nèi)顆粒間從無接觸應(yīng)力到產(chǎn)生有效應(yīng)力的過程。Imai［10-11］認(rèn)為固結(jié)過程主要存在3個(gè)階段：絮凝階段、沉降階段及固結(jié)階段，張楠［12］在Imai研究基礎(chǔ)上進(jìn)一步提出對(duì)于“泥—土系統(tǒng)”最終形成沉積土的演化過程。

通過模型試驗(yàn)?zāi)M分級(jí)堆排的簡(jiǎn)化過程，基于上述對(duì)不同位置孔隙水壓力、沉降時(shí)序變化分析，可將尾礦的沉積固結(jié)過程劃分為5個(gè)階段。

（1）沉降絮凝階段。尾礦漿排放初期，由于水力沖擊、攪動(dòng)及重力作用，尾礦顆粒產(chǎn)生劇烈運(yùn)動(dòng)，此時(shí)液面渾濁。但由于土顆粒密度大于水的密度，呈現(xiàn)下沉運(yùn)動(dòng)并表現(xiàn)出集合沉降的傾向。當(dāng)土顆粒間電解質(zhì)濃度達(dá)到臨界絮凝濃度時(shí)，土顆粒在布朗運(yùn)動(dòng)、不等速沉降、以及水流紊動(dòng)引起的碰撞作用下會(huì)迅速產(chǎn)生絮凝［13］。尾礦漿上下呈現(xiàn)分選沉降層和絮凝沉降層（如圖5所示）。

（2）絮凝壓縮階段。沉降絮凝階段后，絮凝沉降層的尾礦顆粒大多結(jié)成絮團(tuán)并沉底（如圖6所示）。隨著排水的進(jìn)行，絮團(tuán)在底部受重力和布朗運(yùn)動(dòng)的雙重作用，開始相互緊密接觸，絮團(tuán)間的毛細(xì)水被擠出，膠體顆粒間的電磁力不斷增加，同時(shí)也伴隨著孔隙的形成［14］。最終，絮團(tuán)全部轉(zhuǎn)換為絮體，沉降區(qū)消失。

（3）絮體壓縮固結(jié)階段。隨著絮凝壓縮的進(jìn)行，上部壓力使得絮體距離逐漸減小，絮體區(qū)的下部開始固結(jié)形成結(jié)構(gòu)性沉積土（如圖7所示），絮體間的電磁力連接開始轉(zhuǎn)變?yōu)橥令w粒間有效應(yīng)力［15-16］。排水通道附近的沉積土在水頭差的作用下開始通過排水體濾膜排水。隨著滲流的進(jìn)行，尾礦土的孔隙不斷減小，固結(jié)程度逐漸提高。

（4）覆水自重固結(jié)階段。絮體在上覆壓力及自重作用下全部壓縮并轉(zhuǎn)化為結(jié)構(gòu)性沉積土（如圖8所示）。尾礦泥漿完成水土分離，上層覆水通過豎向排水體快速排出，超靜孔隙水壓力逐漸消散，尾礦體積不斷壓縮，土顆粒間有效應(yīng)力逐漸增大。由于排水速率快，該階段尾礦的固結(jié)速度最高。

（5）自重壓縮固結(jié)階段。上層覆水排干后（試驗(yàn)中為35 d后），尾礦進(jìn)入自重壓縮固結(jié)階段。隨著水位逐漸降低到尾礦泥表面以下位置，滲流場(chǎng)轉(zhuǎn)變?yōu)橐酝令w粒中孔隙水滲流為主要滲流形式，滲流速率及土體壓縮速率減緩。隨著排水過程的持續(xù)，尾礦上部呈現(xiàn)非飽和土區(qū)并逐漸向下加深。最后在自重作用下土體壓縮至穩(wěn)定，孔隙水壓力及沉降幾乎不變化，自重固結(jié)完成（如圖9所示）。

3 滲流模型及最大影響半徑

3.1 尾礦理想孔隙滲流模型建立

水頭差使飽和尾礦泥內(nèi)水體在尾礦孔隙內(nèi)流動(dòng)，水體所受土顆粒浮力的反力與水滲流所受阻力互為作用力與反作用力。沿滲流流線方向取一微分單元體，建立理想孔隙滲流模型（如圖10所示）。

從圖10模型中取獨(dú)立圓柱形微分體（下稱微分體）進(jìn)行受力分析，設(shè)滲流起始端a處水頭為h1、面積為A1；滲流終端b處水頭為h2、面積為A2；孔隙流體長(zhǎng)度為L(zhǎng)。作用在微分體的截面a、b處的邊界壓力（孔隙水壓力）分別為：

式中，γw為水容重；其中PWa方向與滲流同向，PWb方向與滲流反向。

水體自身重力沿ab方向的分力為：

式中，γw為水容重；n為土的孔隙率；l為微分體長(zhǎng)度；A為微分體截面積；α為滲流方向與豎直面夾角。

水體所受土顆粒浮力的反力等于與土顆粒同體積水重，其在ab向上的分力為：

水滲流時(shí)微分體土顆粒對(duì)其施加總阻力為：

式中，T為單位水滲流受到的阻力，其方向與滲流反向。

微分體受力如圖11所示，建立其平衡方程：

引入單位尾礦上的滲透力Gcr，其與阻力T互為作用力與反作用力。同時(shí)，將孔隙水壓力u1引入平衡方程式（6），得：

移項(xiàng)整理后，可得孔隙水壓力u1為：

由式（8）可知，滲流起始位置的孔隙水壓力值u1與水位h2、孔隙滲流面積A1及A2、孔隙率n及單位尾礦上的滲透力Gcr各變量呈非線性關(guān)系，隨各變量離散性變化而不斷改變。

3.2 孔隙水壓力變化機(jī)制分析

尾礦泥漿沉降絮凝初期，懸浮狀顆粒間不產(chǎn)生接觸應(yīng)力，僅受自身重力及水浮托力γwVS作用并緩慢下沉。沉降過程中，顆粒間由懸浮狀態(tài)或邊—邊接觸逐漸過渡為面—邊、面—面形式接觸，同時(shí)伴隨著大孔隙和中等孔隙發(fā)育；隨著排水固結(jié)的持續(xù)進(jìn)行，土顆粒開始接觸，有效應(yīng)力不斷增加，尾礦顆粒逐漸由蜂窩狀結(jié)構(gòu)向更為致密的絮凝狀結(jié)構(gòu)轉(zhuǎn)化，土的骨架初步形成，骨架中孔隙構(gòu)成穩(wěn)定的滲流通道。

取排水過程的某一時(shí)段進(jìn)行分析，假設(shè)該時(shí)段的初始孔隙水壓力為u10，結(jié)束時(shí)的孔隙水壓力為u11。根據(jù)式（8）可得孔隙水壓力分別為：

式中符號(hào)與式（8）相同。

為方便討論，現(xiàn)做出以下假設(shè)：

（1）滲流路徑近似為水平向，即 cosα=0；

（2）所研究的滲流孔隙路徑始于各孔隙水壓力計(jì)，終于排水管壁，則式（9）、（10）可簡(jiǎn)化為：

由于排水體處于整個(gè)尾礦泥中心且上下貫通，孔隙水滲流進(jìn)入豎直排水體時(shí)將由水平運(yùn)動(dòng)變?yōu)榇瓜蜻\(yùn)動(dòng)，故可用井損折減系數(shù)β替代試驗(yàn)中水流偏轉(zhuǎn)產(chǎn)生的水頭損失，則滲流末端的水頭高度可表示為：

式中，hw為豎向排水體水位高程?？芍紫端畨毫ψ兓禐椋?/p>

由式（14）可見，尾礦孔隙及水位差變化使得孔隙水壓力逐漸消散。

尾礦泥漿堆排初期，豎向排水通道內(nèi)水位較低，孔隙的兩頭形成水位差，從而使孔隙水沿著孔隙發(fā)生滲流。隨著排水過程的進(jìn)行以及上覆荷載的持續(xù)作用，孔隙不斷被壓縮，等效截面積A減小，加上孔隙間水頭hw也同步減小，尾礦滲流越發(fā)困難，孔隙水壓力逐漸趨于穩(wěn)定。

3.3 最大影響半徑確定

最大影響半徑邊界處不發(fā)生滲流，即孔隙水壓力變化梯度為0，可表達(dá)為：

式中，Δu為孔隙水壓力消散量；l為該位置與豎向排水通道的距離。

將式（14）Δu代入式（15）得：

由2.2節(jié)可知，式（16）中孔隙面積A及水位高度均會(huì)隨固結(jié)過程而逐漸減小，所以最大影響半徑亦是隨時(shí)間而變化的函數(shù)。根據(jù)最大影響半徑定義［17-18］，隨著滲流長(zhǎng)度的增加，尾礦孔隙水壓力下降值逐漸變小。當(dāng)孔隙水壓力下降值不隨距離變化時(shí)，該位置與排水體距離即為排水體最大影響半徑。為了解尾礦固結(jié)過程中最大影響半徑變化特征，以天為單位對(duì)尾礦各個(gè)位置孔隙水壓力下降值進(jìn)行統(tǒng)計(jì)分析?？紤]到底部尾礦受不透水邊界條件影響，滲流條件與上部不一致，故只對(duì)0.5 m及1.0 m堆積高度的4#～9#孔隙水壓力計(jì)進(jìn)行分析。根據(jù)之前測(cè)試數(shù)據(jù)，自第35 d后各層孔隙水壓力變化變緩，即影響半徑變化趨于穩(wěn)定，故僅對(duì)堆排期第11～21 d及停排期的第22～35 d進(jìn)行統(tǒng)計(jì)。

通過對(duì)試驗(yàn)箱半徑內(nèi)的不同位置孔隙水壓力下降值統(tǒng)計(jì)分析，將各斷面孔隙水壓力下降值曲線進(jìn)行擬合預(yù)測(cè)以確定其最大影響半徑值。選用不同函數(shù)擬合后發(fā)現(xiàn)，當(dāng)選用指數(shù)函數(shù)（y=ae-ex）擬合時(shí)R2值最高，各擬合曲線的R2在0.841 8～0.995之間，曲線的擬合值與實(shí)際值相關(guān)性較好。故選用指數(shù)函數(shù)擬合每天0.5 m及1.0 m堆積高度斷面孔隙水壓力下降值曲線，得到最大影響半徑，統(tǒng)計(jì)后繪于圖12。

由圖12可見，最大影響半徑隨著堆排時(shí)間增加，總體上呈逐漸減小趨勢(shì)。但由于孔隙水壓力受溫度、大氣壓等因素影響，加上各曲線擬合R2值不一所形成的可靠性差異，最大影響半徑時(shí)序變化曲線存在一定波動(dòng)性。堆排第1 d 0.5 m及1.0 m堆積高度排水體影響半徑分別達(dá)到8.34 m及5.46 m，下部最大影響半徑明顯大于上部；之后最大影響半徑迅速衰減，并在第4 d出現(xiàn)明顯拐點(diǎn)，減小速率變緩；堆排完成15 d后，0.5 m及1.0 m堆積高度排水體影響半徑分別穩(wěn)定在約2.2 m及1.7 m。這印證了排水體最大影響半徑變化主要是受尾礦固結(jié)程度的影響。固結(jié)前期尾礦孔隙較大，排水體可以通過較長(zhǎng)的滲流通道作用；而隨著上覆荷載的增加，尾礦的固結(jié)程度不斷提高，骨架被逐漸壓縮，滲流孔隙通道變得狹窄或是被滲流過程土顆粒阻斷，滲流難度增加，滲流場(chǎng)范圍逐漸減小。

4 結(jié)論

（1）從縱向上看，孔隙水壓力值、沉降變化與堆積深度成正比，但二者增幅并不相同；從徑向上看，同一高度的孔隙水壓力與沉降值隨滲流路徑增加降幅均逐漸減小，且隨深度增加該規(guī)律越發(fā)顯著。

（2）停排后（第21 d后）上覆荷載逐漸減小，使下層沉降增幅明顯減??；而孔隙水壓力則是在停排后呈兩段式減小，在上層自由水排干后（第35 d），孔隙水壓力降幅明顯減小。說明孔隙水壓力與沉降在土體排水固結(jié)過程中變化并不完全同步，應(yīng)力變化滯后于應(yīng)變變化。

（3）根據(jù)沉降及孔隙水壓力變化曲線可進(jìn)一步將逐級(jí)堆排尾礦排水過程分為沉降絮凝階段、絮凝壓縮階段、絮體壓縮固結(jié)階段、覆水自重固結(jié)階段及自重壓縮固結(jié)階段5個(gè)階段。

（4）建立了尾礦理想孔隙滲流模型，孔隙水壓力與水位、孔隙滲流面積、孔隙率及單位尾礦上的滲透力呈非線性相關(guān)，并基于孔隙水壓力數(shù)據(jù)擬合測(cè)算各時(shí)期不同高度排水體最大影響半徑。隨著尾礦的固結(jié)程度不斷提高，滲流場(chǎng)范圍逐漸減小，影響半徑在約4 d內(nèi)迅速衰減，第35 d后0.5 m及1.0 m堆積高度影響半徑分別降低至約1.7 m及2.2 m。