歐孝奪，秦金喜，全守岳，江 杰*，楊杉楠

(1.廣西大學(xué)土木建筑工程學(xué)院，南寧 530004；2.廣西大學(xué)工程防災(zāi)與結(jié)構(gòu)安全教育部重點(diǎn)實(shí)驗(yàn)室，南寧 53004；3.廣西金屬尾礦安全防控工程技術(shù)研究中心，南寧 530004；4.中鐵二院工程集團(tuán)有限責(zé)任公司南寧勘察設(shè)計(jì)研究院，南寧 530012)

隨著中國(guó)經(jīng)濟(jì)快速增長(zhǎng)，金屬工業(yè)產(chǎn)品的需求逐年上升，而開采量的增加導(dǎo)致排泥庫數(shù)量和庫容大幅增加。排泥庫作為金屬冶煉廢棄物堆存設(shè)施，是金屬非金屬礦山的重大危險(xiǎn)源。特別在巖溶地區(qū)修建排泥庫時(shí)，地下溶洞交錯(cuò)連接，常因底部溶洞未探查清楚，致使排泥庫在長(zhǎng)周期運(yùn)行后或遇暴雨、地震等其他自然因素作用下發(fā)生泄漏事故。排泥庫泥漿富含重金屬離子，一旦發(fā)生泄漏事故，庫內(nèi)的廢渣、廢水就會(huì)污染地下水，甚至淹沒下游的農(nóng)田村莊，不僅破壞生態(tài)環(huán)境，同時(shí)給人民的生命財(cái)產(chǎn)帶來巨大的安全隱患[1-3]。據(jù)宋園生[4]統(tǒng)計(jì)，溶洞塌陷占造成排泥庫泥漿泄漏事故直接原因的48.1%，因而開展巖溶區(qū)排泥庫深層泥漿的泄漏機(jī)理研究非常有必要，為后期尾礦庫防治工作提供理論基礎(chǔ)。

為研究流體的泄漏機(jī)理，國(guó)內(nèi)外學(xué)者進(jìn)行了大量有益的研究。Ger-manopoulos等[5]、Lin[6]、Alonso等[7]和中國(guó)城鎮(zhèn)供水協(xié)會(huì)[8]分別建立了不同情況下的管網(wǎng)泄漏量計(jì)算模型，在一定條件下為不同區(qū)域的流體泄漏提供了計(jì)算方法，但并未針對(duì)巖溶區(qū)尾礦庫區(qū)泥漿流體的泄漏提供相應(yīng)的計(jì)算模型。王勇智[9]認(rèn)為泥石流這類流體的濃度有一定的變化范圍，其相對(duì)黏滯系數(shù)與固體體積濃度、顆粒大小均有關(guān)系。徐黎明等[10]研究得出影響泥石流流速的因素包括泥沙粒徑、泥漿深度、泥漿密度等，并利用此參數(shù)建立BP(back propagation)神經(jīng)網(wǎng)絡(luò)。費(fèi)祥俊等[11]研究了泥石流中不同顆粒大小的運(yùn)動(dòng)軌跡等參數(shù)，證明固液兩相流的運(yùn)動(dòng)形式受顆粒大小影響，并確定了兩相流的分界粒徑。沈壽長(zhǎng)[12]研究了應(yīng)該采取何種計(jì)算模型才能準(zhǔn)確反映泥石流固液相的相互作用，并提出了相應(yīng)的本構(gòu)關(guān)系方程。泥石流等流體相關(guān)性質(zhì)雖同尾礦泥漿相似，但目前國(guó)內(nèi)外學(xué)者并未針對(duì)尾礦泥漿深層泄漏機(jī)理進(jìn)行相應(yīng)研究。宋園生[4]分析了處于巖溶區(qū)的排泥庫面臨溶洞塌陷時(shí)采取的險(xiǎn)情處置方式，給出了巖溶地區(qū)排泥庫在建設(shè)和運(yùn)行過程中的建議并得到相應(yīng)的應(yīng)用。但現(xiàn)目前，國(guó)內(nèi)外針對(duì)排泥庫的泥漿的探測(cè)研究鮮有報(bào)道，而針對(duì)大壩、水庫等處于巖溶地區(qū)中的設(shè)施泥漿泄漏通常采用地球物理法[13]、綜合示蹤法[14-16]等進(jìn)行檢測(cè)，雖方法效果顯著但較為復(fù)雜。對(duì)于流體泄漏，學(xué)者主要集中研究泥石流潰壩及水庫泄漏或僅針對(duì)喀斯特地區(qū)尾礦庫泄漏所存在的災(zāi)情處置手段等方面，而針對(duì)巖溶區(qū)排泥庫底部泥漿泄漏規(guī)律的研究仍處于空白。

為解決巖溶區(qū)排泥庫底部泄漏的快速探測(cè)問題，自主研發(fā)了一種探測(cè)裝置[17]，利用該探測(cè)裝置所測(cè)的試驗(yàn)數(shù)據(jù)運(yùn)用灰色關(guān)聯(lián)建立了泥漿參數(shù)計(jì)算模型，并用ANSYS進(jìn)行泄漏孔直徑、泥漿含水量、動(dòng)力黏度對(duì)泄漏點(diǎn)的影響半徑數(shù)值模擬研究[18]。該研究創(chuàng)新性提出新型排泥庫泄漏探測(cè)方式，建立了泥漿對(duì)探測(cè)器拉力的多元線性回歸計(jì)算模型，對(duì)快速準(zhǔn)確查找排泥庫泄漏點(diǎn)位置有極大幫助。通過數(shù)值模擬計(jì)算了深層泥漿的不同物理性質(zhì)對(duì)影響半徑的影響作用。本文對(duì)于排泥庫深層泄漏點(diǎn)快速探測(cè)及明確泄漏區(qū)域泥漿流動(dòng)狀態(tài)都具有重要的理論指導(dǎo)意義。

1 試驗(yàn)設(shè)計(jì)

試驗(yàn)主要測(cè)量泥漿探測(cè)裝置所受滲流力與泥漿流速、泥漿含水量和泥漿密度之間的關(guān)系，以期能尋找規(guī)律并建立相關(guān)理論。

1.1 試驗(yàn)裝置

試驗(yàn)裝置主要由測(cè)量裝置、流動(dòng)裝置、探測(cè)裝置、儲(chǔ)存裝置等組成，如圖1所示。

圖1 試驗(yàn)裝置模型圖

1.1.1 泥漿滲流力測(cè)量裝置

泥漿滲流力采用WH-A04L型數(shù)顯拉力計(jì)進(jìn)行測(cè)量，測(cè)量精度為0.001 kg。數(shù)顯拉力計(jì)下端連接自制的泥漿探測(cè)器。泥漿探測(cè)器采用聚氯乙烯塑料塊加工制成，為增大泥漿探測(cè)器和流動(dòng)泥漿的摩擦阻力，使泥漿對(duì)探測(cè)器的拉力作用更加敏感，將探測(cè)器做成扇貝形狀，并在探測(cè)器表面進(jìn)行粗糙條紋處理，探測(cè)器的內(nèi)部進(jìn)行機(jī)械掏空，可通過填充不同重量的填充物改變其密度，使之沉到不同的泥漿層中，探測(cè)器如圖2所示。試驗(yàn)中為使泥漿流速計(jì)算誤差減小，使之更接近出口處泥漿流速，探測(cè)器放置在聚氯乙烯(PVC)管底部靠近孔口處。

圖2 可調(diào)重量探測(cè)器

1.1.2 泥漿流動(dòng)裝置

泥漿流動(dòng)裝置由內(nèi)徑150 mm，長(zhǎng)度1 000 mm的PVC白色塑料管組成。將塑料管豎直放置，在管底使用塑料蓋進(jìn)行封底。為得到不同的泥漿流速，對(duì)蓋子進(jìn)行不同直徑大小的鉆孔，孔口直徑分別設(shè)為10、15、20、25、30、35 mm。因排泥庫泥漿隨著深度的增加，含水量有很大區(qū)別，研究表明排泥庫泥漿靜置抽干上部清水后的表層泥漿含水量基本維持在120%，排泥庫發(fā)生泄漏的泥漿一般為新排入排泥庫的渾濁泥漿，由于遠(yuǎn)距離泵送功量需要，吹填入庫的泥漿需保證一定的流動(dòng)性，使所吹填泥漿含水量比靜置的泥漿要高很多，根據(jù)此情況對(duì)泥漿進(jìn)行分組，分別配置含水量為100%、110%、120%、130%、140%、150%、160%、170%、180%的泥漿。

1.2 試驗(yàn)步驟

測(cè)量泥漿泄漏滲流力的試驗(yàn)步驟如下。

(1)采用常規(guī)土工密度實(shí)驗(yàn)對(duì)不同含水量的泥漿進(jìn)行密度測(cè)量，得到相應(yīng)的數(shù)值。試驗(yàn)前，將管底泄漏孔口預(yù)先用東西封住，往PVC管倒?jié)M一定含水量的泥漿。

(2)打開底部泄漏孔，同時(shí)開始計(jì)時(shí)并進(jìn)行上部注漿，保持泥漿水頭不變，進(jìn)而控制底下流速保持一致。由于流速較快，故只記錄前面5 s的流量，5 s后停止泄漏，將泄漏出來的泥漿進(jìn)行稱重，記錄5 s鐘泄漏的流量和上面的拉力讀數(shù)。

(3)重復(fù)上面的試驗(yàn)，更換不同含水量的泥漿和不同大小的泄漏孔口，得到正交數(shù)據(jù)。

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

2.1 關(guān)聯(lián)度分析

為了解泥漿的各物理力學(xué)參數(shù)與泥漿滲流力探測(cè)裝置的關(guān)聯(lián)程度，通過對(duì)多種方法進(jìn)行比較分析，選擇灰色關(guān)聯(lián)分析法從泥漿的各物理力學(xué)參數(shù)中篩選出與探測(cè)裝置所受泥漿滲流力關(guān)聯(lián)度較高的影響因子作為定量因子。

通過以上室內(nèi)土工實(shí)驗(yàn)和泄漏模型試驗(yàn)，獲取了大量的排泥庫泥漿物理力學(xué)性質(zhì)數(shù)據(jù)和泥漿泄漏流動(dòng)數(shù)據(jù)，共計(jì)54組原始數(shù)據(jù)。

原始數(shù)據(jù)處理具體分析步驟如下。

2.1.1 流速換算

由于試驗(yàn)所測(cè)數(shù)據(jù)為泥漿流量，需要將其換算為探測(cè)器位置的泥漿流速。換算公式為

(1)

式(1)中：v為泥漿流速，cm/s；m為泄漏出來的泥漿質(zhì)量，g；ρ為泥漿密度，g/cm3；t為時(shí)間，此處取t=5 s；A為內(nèi)徑150 mm的PVC圓管截面面積。

2.1.2 確定參考數(shù)列和比較數(shù)列

能反映系統(tǒng)行為特征的數(shù)據(jù)序列稱為參考數(shù)列X0，本文參考數(shù)列需要代表泥漿探測(cè)裝置的受力指標(biāo)，此處選擇泥漿探測(cè)器所受滲流力大小作為參考數(shù)列X0。影響系統(tǒng)行為的因素組成的數(shù)據(jù)序列，稱比較數(shù)列，此處選擇泥漿流動(dòng)參數(shù)和泥漿物理性質(zhì)指標(biāo)作為比較數(shù)列Xi，分別為泥漿的流速v(X1)、天然含水量ω(X2)、天然密度ρ(X3)。

2.1.3 原始數(shù)據(jù)無量綱化

由于系統(tǒng)中各因素的物理意義不同，導(dǎo)致數(shù)據(jù)的量綱不相同，不便于比較計(jì)算，難以得出正確的結(jié)論。因此在進(jìn)行灰色關(guān)聯(lián)度分析時(shí)，需要對(duì)數(shù)據(jù)進(jìn)行無量綱化處理后才可進(jìn)行關(guān)聯(lián)性分析。常用的無量綱化處理方法有三種，分別為初值化、均值化和區(qū)間值化處理。經(jīng)過分析，本次試驗(yàn)數(shù)據(jù)沒有“0”，且含水量是從100%開始計(jì)算，采用初值法進(jìn)行無量綱化處理，可以使計(jì)算簡(jiǎn)單、不易出錯(cuò)，故對(duì)排泥庫泄漏試驗(yàn)所測(cè)原始數(shù)據(jù)使用初值法進(jìn)行無量綱化處理。

2.1.4 對(duì)序列求差處理

計(jì)算參考數(shù)列與各比較序列差值的絕對(duì)值，計(jì)算公式為

Δi(k)=|X0(K)-Xi(K)|

(2)

求出差值絕對(duì)值的最大值和最小值，計(jì)算公式如式(3)、式(4)。

max[maxΔi(k)]=1.52

(3)

min[minΔi(k)]=0.00

(4)

2.1.5 求關(guān)聯(lián)系數(shù)與關(guān)聯(lián)度

用式(5)求參考數(shù)列與比較數(shù)列的關(guān)聯(lián)系數(shù)。

(5)

因?yàn)殛P(guān)聯(lián)系數(shù)是比較數(shù)列與參考數(shù)列在各個(gè)時(shí)刻的關(guān)聯(lián)程度值，所以它的數(shù)不止一個(gè)，而信息過于分散不便于進(jìn)行整體性比較。因此有必要將關(guān)聯(lián)系數(shù)集中為一個(gè)值，即求其平均值，作為比較數(shù)列與參考數(shù)列間關(guān)聯(lián)程度的數(shù)量表示。

通過上述數(shù)據(jù)處理手段進(jìn)行試驗(yàn)原始數(shù)據(jù)的處理得到如圖3所示的各因素與滲流力關(guān)聯(lián)系數(shù)。

圖3 各因素與滲流力關(guān)聯(lián)度

用式(6)求解關(guān)聯(lián)度ri：

(6)

式(6)中：ri為X0與Xi的關(guān)聯(lián)度，其值越接近1，則關(guān)聯(lián)度越好。各指標(biāo)與泥漿所受滲流力關(guān)聯(lián)度分別為r1=0.818 1、r2=0.705 2、r3=0.846 7，由此可知各關(guān)聯(lián)度大小順序?yàn)閞3>r1>r2。

由灰色關(guān)聯(lián)分析結(jié)果可知，排泥庫泥漿含水量、密度和流速三者對(duì)泥漿滲流力的關(guān)聯(lián)系數(shù)均較高，其中排泥庫泥漿含水量與泥漿探測(cè)器所受滲流力關(guān)聯(lián)度較低，而泥漿密度和泥漿探測(cè)器所受滲流力關(guān)聯(lián)度較高，相互之間的關(guān)系較為密切，均可作為計(jì)算模型考慮的因素。

2.2 泥漿泄漏參數(shù)多元線性回歸計(jì)算模型

通過運(yùn)用多元線性回歸理論對(duì)前期試驗(yàn)數(shù)據(jù)進(jìn)行分析，計(jì)算時(shí)取系數(shù)個(gè)數(shù)p=3，樣本數(shù)量取n=54，根據(jù)數(shù)值分析相關(guān)回歸計(jì)算公式計(jì)算得到：

y=-18.535 5-0.017 3x1+0.013 6x2+

10.802 3x3

(7)

式(7)中：y為泥漿對(duì)探測(cè)器的滲流力大小，N；x1為泥漿流速，：cm/s；x2為泥漿含水量，%；x3為泥漿密度，g/cm3。

這就是排泥庫泥漿深層泄漏時(shí)泥漿對(duì)探測(cè)器的滲流力與泥漿內(nèi)部參數(shù)的計(jì)算模型。

對(duì)此計(jì)算模型進(jìn)行顯著性檢驗(yàn)。運(yùn)用F分布的方差分析法計(jì)算本次數(shù)據(jù)的F=740.667，查F分布上側(cè)分位數(shù)表知F0.05(3，50)<2.84，顯而易見F>F0.05(3，50)，證明此次計(jì)算結(jié)果線性回歸顯著，結(jié)論可靠。

3 數(shù)值模擬分析

在建立泥漿探測(cè)器滲流力與泥漿參數(shù)的計(jì)算模型后，為能快速準(zhǔn)確查到泄漏點(diǎn)位置，還需確定泄漏點(diǎn)的影響半徑。

3.1 模型建立

3.1.1 基本假定

排泥庫泥漿由于在重力作用下發(fā)生沉積，其含水量隨深度發(fā)生變化并出現(xiàn)分層現(xiàn)象，模擬泥漿分層較為復(fù)雜，為簡(jiǎn)化計(jì)算，縮短計(jì)算時(shí)間，故假設(shè)模擬過程泥漿均為含水量單一的介質(zhì)。且運(yùn)用VOF模型時(shí)需控制體積必須充滿單一流體相或存在多相的聯(lián)合但其有一相作為可壓縮的理想氣體，故本模型上部排泥庫庫面以上假設(shè)為可壓縮的理想氣相。壁面處為非滑移邊界條件，熱流量假設(shè)為零熱流量，即壁面假設(shè)為絕熱壁面、光滑壁面。

3.1.2 軟件及模型選取

本次數(shù)值模擬使用ANSYS ICEM CFD作為前處理器，F(xiàn)luent 15.0作為求解器和后處理軟件。該模型模擬的兩相是空氣和泥漿，選用Fluent中的VOF模型，故只能絕對(duì)速度的壓力基瞬態(tài)求解器。模擬泥漿旋轉(zhuǎn)效應(yīng)不是特別明顯，采用常用的標(biāo)準(zhǔn)k-ε湍流模型。由于討論的泥漿流動(dòng)是瞬態(tài)的，采用更適宜瞬態(tài)計(jì)算的PISO算法，收斂標(biāo)準(zhǔn)取10-4，有效提高收斂性和縮短計(jì)算時(shí)間。

3.1.3 建模尺寸

通過對(duì)廣西百色地區(qū)尾礦庫區(qū)調(diào)研發(fā)現(xiàn)，為貼近現(xiàn)場(chǎng)實(shí)際情況，更真實(shí)模擬排泥庫泥漿泄漏情況，模型模擬尺寸為50 m×50 m×30 m。整個(gè)模型模擬的是一個(gè)底部泄漏的排泥庫，模型上部模擬的是排泥庫庫面，模型側(cè)面和底面模擬的是排泥庫壁面，圓孔模擬的是排泥庫泥漿泄漏孔，幾何模型示意如圖1所示。

3.1.4 關(guān)鍵參數(shù)及邊界條件

(1)關(guān)鍵參數(shù)

通過前期查閱相關(guān)國(guó)內(nèi)外學(xué)者對(duì)于流體泄漏研究?jī)?nèi)容以及對(duì)排泥庫現(xiàn)場(chǎng)調(diào)研發(fā)現(xiàn)，排泥庫泄漏點(diǎn)的影響半徑與泄漏孔直徑、含水量、動(dòng)力黏度、泄漏孔形狀、泥漿均勻程度、庫區(qū)是否閉庫等因素有關(guān)。但目前學(xué)者研究流體泄漏影響半徑主要有泄漏孔直徑、含水量、動(dòng)力黏度，而泥漿泄漏方面缺乏研究。因此文章借鑒流體泄漏方面研究，主要考慮泄漏孔直徑、含水量、動(dòng)力黏度對(duì)泥漿泄漏影響半徑的影響。泥漿泄漏孔直徑取1、2、3 m，含水量取120%、150%、180%，動(dòng)力黏度取1、5、10 Pa·s。

(2)邊界條件

空氣、泥漿兩相流中的氣體為理想可壓縮氣體，氣體是可壓縮的，泥漿不可壓縮，空氣隨著泥漿泄漏過程不斷下降，體積變大，密度和壓力均不變，故模型的上表面采用壓力輸入入口。模型的出口存在大氣壓力，設(shè)為壓力出口，模型壁面設(shè)為固體壁面，壁面出為非滑移邊界條件。由于Interface的連接面之間的節(jié)點(diǎn)可以不一致，選用Interface類型。

模型網(wǎng)格底部劃分如圖4所示。

圖4 網(wǎng)格底部

3.2 模擬結(jié)果分析

數(shù)值模擬對(duì)泄漏點(diǎn)主要影響因素不同大小泄漏孔口、泥漿含水量、動(dòng)力黏度的分析結(jié)果如下。

3.2.1 泥漿在不同大小泄漏孔發(fā)生泄漏時(shí)的影響半徑分析

為了更真實(shí)地模擬排泥庫泥漿的泄漏情況，本次模擬的排泥庫尺寸取50 m×50 m×30 m，泥漿的含水量設(shè)為120%，密度數(shù)值大小采用試驗(yàn)所測(cè)數(shù)值，見表1，泄漏孔直徑大小分別設(shè)為1、2、3 m。

表1 泥漿含水量與密度的關(guān)系

由于泥漿探測(cè)器能探測(cè)到的最小泥漿速度數(shù)值未知，泥漿本身具有一定的動(dòng)力黏度，流動(dòng)現(xiàn)象沒有清水那么明顯，而流速太小又不能被探測(cè)到，故設(shè)置泥漿流速0.5 m/s作為流動(dòng)半徑的界定標(biāo)準(zhǔn)。為同模型試驗(yàn)相對(duì)應(yīng)，使研究規(guī)律較明顯，選取泄漏10 s時(shí)的狀態(tài)作為研究對(duì)象。根據(jù)觀測(cè)模型試驗(yàn)狀態(tài)和現(xiàn)場(chǎng)實(shí)際結(jié)合，推測(cè)在泄漏孔泄漏時(shí)間在20、60 s…后，由于泥漿介質(zhì)相對(duì)均質(zhì)單一，泄漏孔的影響半徑會(huì)保持不變。而當(dāng)泥漿下沉到泄漏孔的影響范圍時(shí)，出現(xiàn)一些肉眼可見的現(xiàn)象：表層逐漸出現(xiàn)漩渦，由外而內(nèi)流速逐漸變大，在泄漏孔中心流速最大，下沉最多，呈明顯的漏斗狀分布。圖5是模擬排泥庫泥漿在不同泄漏孔直徑下泄漏10 s時(shí)的速度云圖。

圖5 不同直徑泄漏孔泥漿速度云圖

根據(jù)以上對(duì)不同泄漏孔直徑模擬結(jié)果，測(cè)出影響半徑如表2所示。

表2 模擬影響半徑結(jié)果

由圖5、表2可知，在泥漿高度一樣，孔口大小不一的情況下，洞口處的泄漏泥漿流速大小有所差別，洞口越大，流速越大。由于泥漿存在一定的動(dòng)力黏度，在排泥庫底部發(fā)生泄漏時(shí)，泥漿泄漏孔存在一定的影響半徑，水平影響半徑平均為泄漏孔半徑的4.61倍，豎直影響半徑平均為泄漏孔半徑的5倍，泄漏孔中心流速最大，在庫深30 m，孔口大小為1～3 m情況下，最大流速平均為22.8 m/s。由此能解釋排泥庫泥漿泄漏時(shí)，泥漿表面沒有漩渦和漏斗、表層水或泥漿流動(dòng)不明顯的原因。

根據(jù)模擬結(jié)果可推測(cè)泥漿深度較大的排泥庫泥漿流動(dòng)過程如下：

(1)泥漿表層整體下沉，表面無漩渦和漏斗。

(2)泥漿下沉到泄漏孔的影響范圍時(shí)，泥漿表層逐漸出現(xiàn)漩渦，由外而內(nèi)流速逐漸變大，在泄漏孔中心流速最大，下沉最多，呈明顯的漏斗狀分布。

(3)泥漿泄漏完畢，由于排泥庫底部泥漿密度和動(dòng)力黏度均較大，含水量較低，泥漿內(nèi)部存在相互作用力，無法完全從泄漏孔流出，在排泥庫底發(fā)生淤積，上部密度和動(dòng)力黏度較小，含水量較大的泥漿則全部從泄漏孔流出。

3.2.2 不同含水量泥漿發(fā)生泄漏時(shí)影響半徑分析

由于排泥庫泥漿從排礦口排出靜置一段時(shí)間后，在重力作用下發(fā)生沉積，泥漿含水量發(fā)生變化，上部泥漿含水量較小，底部泥漿含水量較大。使用軟件模擬泥漿分層較為復(fù)雜，故模擬過程均是采用含水量單一的介質(zhì)。為了研究不同泥漿含水量對(duì)模擬結(jié)果是否產(chǎn)生較大影響，取泄漏孔直徑為3 m的排泥庫模型，分別模擬泥漿在含水量分別為120%、150%、180%情況下的泄漏影響半徑。

由數(shù)值模擬可知，在泄漏孔口直徑為3 m，含水量分別在120%、150%、180%情況下，泄漏孔影響范圍基本相同，三種不同含水量速度云圖大致如圖6所示，只有數(shù)值上有些許差別，證明不同含水量對(duì)排泥庫泥漿泄漏影響范圍不大。

圖6 含水量泥漿速度云圖

3.2.3 不同動(dòng)力黏度的泥漿發(fā)生泄漏時(shí)影響半徑分析

由于泥漿與普通清水參數(shù)不一樣，它具有較大的動(dòng)力黏度，在泥漿內(nèi)部發(fā)生相對(duì)剪切運(yùn)動(dòng)時(shí)，會(huì)產(chǎn)生相互阻礙的力，使得泥漿產(chǎn)生流動(dòng)需要克服泥漿內(nèi)部的阻力。為了研究泥漿在不同的動(dòng)力黏度情況下是否會(huì)對(duì)泥漿泄漏影響范圍產(chǎn)生影響，對(duì)不同動(dòng)力黏度的泥漿泄漏情況進(jìn)行數(shù)值模擬。模擬過程中，取泥漿泄漏孔口大小為3 m，泥漿動(dòng)力黏度μ分別設(shè)置為1、5、10 Pa·s。

由圖7可知，在泄漏孔口直徑3 m，泥漿動(dòng)力黏度μ分別為1、5、10 Pa·s情況下，泄漏孔處泥漿最大速度依舊僅有微小變化，三者數(shù)值分析圖片只有數(shù)值上有些許差別，故只放置了一張?jiān)茍D。從得到的結(jié)果分析泥漿動(dòng)力黏度對(duì)泥漿泄漏影響范圍變化不明顯，說明泥漿動(dòng)力黏度對(duì)泥漿泄漏影響范圍影響不大。

圖7 動(dòng)力黏度泥漿速度云圖

4 結(jié)論

以廣西鋁土排泥庫為研究背景，通過室內(nèi)模型試驗(yàn)和Fluent 15.0流體計(jì)算軟件研究了排泥庫深層泥漿發(fā)生泄漏時(shí)的流動(dòng)機(jī)理，得到了以下結(jié)論。

(1)排泥庫泥漿含水量、密度和流速三者對(duì)泥漿滲流力的關(guān)聯(lián)系數(shù)均較高，均可作為排泥庫泥漿泄漏計(jì)算模型考慮的因素。

(2)運(yùn)用多元線性回歸計(jì)算理論建立了泥漿探測(cè)器所受滲流力與泥漿流速、含水量和密度之間的多元線性回歸模型。并運(yùn)用數(shù)理統(tǒng)計(jì)理論對(duì)此多元線性回歸模型進(jìn)行驗(yàn)證，證明結(jié)果顯著。

(3)通過Fluent 15.0流體計(jì)算軟件計(jì)算了在排泥庫底部發(fā)生泄漏時(shí)，泥漿泄漏孔存在一定的影響半徑，水平影響半徑平均為泄漏孔半徑的4.61倍，豎直影響半徑平均為泄漏孔半徑的5倍。并可以根據(jù)模擬結(jié)果推測(cè)厚度較大的排泥庫泥漿流動(dòng)過程。不同含水量和動(dòng)力黏度對(duì)排泥庫泥漿泄漏范圍影響不大。

(4)通過對(duì)泄漏孔半徑的影響因素分析，影響的最關(guān)鍵因素是泄漏孔半徑，而泥漿含水量和動(dòng)黏度系數(shù)則對(duì)泄漏孔影響半徑不大。隨著排泥庫底部泄漏孔半徑的增加當(dāng)泥漿發(fā)生泄漏時(shí)其影響半徑最大，會(huì)出現(xiàn)泥漿表面下沉、產(chǎn)生旋渦等現(xiàn)象。因此在進(jìn)行排泥庫選址時(shí)應(yīng)充分探測(cè)熟知該地區(qū)的地質(zhì)構(gòu)造情況，必要時(shí)對(duì)危害排泥庫安全的溶洞進(jìn)行施工處理，這在源頭上進(jìn)行防控和治理有極大的意義。

首次對(duì)排泥庫泥漿泄漏機(jī)理進(jìn)行研究，采用自制探測(cè)裝置對(duì)排泥庫泥漿深層泄漏點(diǎn)的快速查找，彌補(bǔ)了排泥庫泄漏機(jī)理研究領(lǐng)域的空白，研究成果的實(shí)際應(yīng)用將大幅減少泄漏事故的危害和損失，對(duì)排泥庫泄漏的安全防控和后期治理起到極大的指導(dǎo)意義。