梁遠(yuǎn)驥,方前程,劉奇
(1.廣東省冶金建筑設(shè)計(jì)研究院有限公司, 廣東 廣州 510080;2.黃淮學(xué)院建筑工程學(xué)院, 河南 駐馬店市 463000;3.長(zhǎng)沙礦山研究院有限責(zé)任公司, 湖南 長(zhǎng)沙 410012;4.國(guó)家金屬采礦工程技術(shù)研究中心, 湖南 長(zhǎng)沙 410012)
尾礦庫(kù)是礦山企業(yè)不可或缺的基礎(chǔ)設(shè)施,但同時(shí)也是一種潛在的人造重大危險(xiǎn)源[1]。隨著礦山機(jī)械化水平的提升,礦產(chǎn)資源開采規(guī)模不斷擴(kuò)大,尾礦庫(kù)數(shù)量也在不斷攀升。我國(guó)尾礦庫(kù)呈現(xiàn)總量大、中小型庫(kù)占比高、安全等級(jí)系數(shù)普遍偏低、上游式筑壩多等突出特點(diǎn),并且尾礦庫(kù)下游一般為人口、經(jīng)濟(jì)密集區(qū),尾礦庫(kù)壩體一旦失穩(wěn),會(huì)對(duì)地下水體和土壤環(huán)境造成嚴(yán)重的污染和破壞,造成重大人員傷亡和不可估量的經(jīng)濟(jì)損失[2]。降雨是影響尾礦壩穩(wěn)定性的關(guān)鍵因素[3],因此,科學(xué)合理地開展降雨條件下的尾礦壩穩(wěn)定性評(píng)估具有重大意義。
目前,眾多專家和學(xué)者針對(duì)降雨對(duì)尾礦庫(kù)壩體穩(wěn)定性的影響開展了大量研究。何錦龍等[4]采用大型數(shù)值模擬軟件,研究了正常水位及20 a、50 a、100 a一遇日最大降雨量4種工況下尾礦壩穩(wěn)定性狀況,研究顯示降雨量與壩體穩(wěn)定性狀態(tài)呈線性相關(guān)關(guān)系;陳生水等[5]以滲流理論為基礎(chǔ),以連續(xù)降雨為研究條件,建立了尾礦壩體失穩(wěn)潰壩理論模型,結(jié)果顯示降雨入滲引起壩體浸潤(rùn)線抬升是壩體失穩(wěn)的主要誘因,為尾礦庫(kù)壩體穩(wěn)定性評(píng)估提供了理論依據(jù);隨著數(shù)學(xué)理論模型研究不斷深入,包括集對(duì)分析理論[6]、可變模糊集[7]、物元-可拓學(xué)[8]、因素空間-未確知理論[9]等數(shù)學(xué)評(píng)價(jià)模型逐漸應(yīng)用到尾礦庫(kù)壩體穩(wěn)定性評(píng)估中,取得了較為滿意的應(yīng)用效果。經(jīng)過對(duì)比發(fā)現(xiàn),上述研究成果均以確定的靜態(tài)權(quán)重評(píng)估不同降雨條件下尾礦壩的穩(wěn)定性,忽略了降雨條件改變導(dǎo)致的其他評(píng)價(jià)指標(biāo)對(duì)尾礦壩穩(wěn)定性貢獻(xiàn)度的變化。實(shí)際上,尾礦壩穩(wěn)定性影響因素眾多且關(guān)系復(fù)雜,特別是在降雨條件影響下,其他風(fēng)險(xiǎn)因素的發(fā)生程度和概率也會(huì)受到影響。因此,以往的評(píng)價(jià)模型評(píng)估結(jié)果具有一定片面性,甚至可能出現(xiàn)評(píng)估結(jié)果與實(shí)際情況不符,需要確定不同降雨條件下的動(dòng)態(tài)指標(biāo)權(quán)重,評(píng)估結(jié)果才更符合真實(shí)狀況。
鑒于此,本文在構(gòu)建尾礦庫(kù)壩體穩(wěn)定性綜合評(píng)價(jià)指標(biāo)體系基礎(chǔ)上,引入屬性識(shí)別理論,基于G1法建立了不同降雨條件下的動(dòng)態(tài)評(píng)價(jià)模型。以某尾礦庫(kù)為例,驗(yàn)證評(píng)估模型的有效性,以期進(jìn)一步為降雨條件下尾礦壩穩(wěn)定性評(píng)價(jià)提供理論參考。
尾礦庫(kù)壩體穩(wěn)定性影響因素眾多,且指標(biāo)之間存在復(fù)雜的關(guān)聯(lián)性,因此,在指標(biāo)選取時(shí),要兼顧適用、實(shí)用、全面、客觀、系統(tǒng)、可測(cè)算等原則。在參考大量文獻(xiàn)[6-9]、查閱行業(yè)規(guī)范、實(shí)地調(diào)研和咨詢相關(guān)專家的基礎(chǔ)上,從降雨條件、工程因素、管理因素和地質(zhì)條件4個(gè)方面,構(gòu)建尾礦壩穩(wěn)定性綜合評(píng)價(jià)指標(biāo)體系[10-11]。將尾礦壩穩(wěn)定性劃分為5個(gè)風(fēng)險(xiǎn)等級(jí)。指標(biāo)體系及具體分級(jí)標(biāo)準(zhǔn)見表1。
表1 尾礦壩穩(wěn)定性評(píng)價(jià)指標(biāo)分級(jí)標(biāo)準(zhǔn)
2.1 G1法計(jì)算指標(biāo)主觀權(quán)重
郭亞軍在傳統(tǒng)層次分析法的基礎(chǔ)上,優(yōu)化改進(jìn)了判斷矩陣構(gòu)建方式,提出了G1法,該方法是由專家先根據(jù)指標(biāo)相對(duì)重要性程度對(duì)指標(biāo)進(jìn)行排序,并給定各指標(biāo)之間相對(duì)重要性比值來確定指標(biāo)權(quán)重的主觀賦權(quán)法。具體計(jì)算過程如下[12]。
(1)建立指標(biāo)序列集合。為綜合分析待評(píng)對(duì)象,假設(shè)選取了n個(gè)評(píng)價(jià)指標(biāo),由專家結(jié)合主觀經(jīng)驗(yàn),在對(duì)指標(biāo)進(jìn)行兩兩相對(duì)重要性分析比較之后,在所有指標(biāo)中挑選出最重要指標(biāo),記為X1,其指標(biāo)權(quán)重記為α1,依次在剩余指標(biāo)中繼續(xù)挑選,直到所有指標(biāo)挑選完成,最終形成了按指標(biāo)相對(duì)重要性程度排序的指標(biāo)序列(X1,X2,…,Xn),對(duì)應(yīng)的指標(biāo)權(quán)重可以表示為(α1,α2,…,αn)。
(2)確定指標(biāo)相對(duì)重要性比值。相對(duì)重要性比值r j體現(xiàn)了兩兩指標(biāo)之間對(duì)評(píng)估結(jié)果的重要性程度,具體取值相對(duì)重要性情況描述見表2。計(jì)算公式為:
表2 相對(duì)重要性比值
(3)計(jì)算第n個(gè)指標(biāo)權(quán)重。第n個(gè)指標(biāo)權(quán)重計(jì)算公式為:
(4)計(jì)算所有指標(biāo)權(quán)重。根據(jù)相對(duì)重要性比值rj,確定其他指標(biāo)權(quán)重。計(jì)算公式為:
2.2 屬性識(shí)別模型
屬性識(shí)別模型基本原理是:假設(shè)有n個(gè)待評(píng)對(duì)象,選取了m個(gè)指標(biāo)I j(j=1,2,…,m),劃分了p個(gè)評(píng)價(jià)等級(jí)Ck(k=1,2,…,p),按照單指標(biāo)分級(jí)屬性特點(diǎn),構(gòu)建不同等級(jí)Ck下的單指標(biāo)屬性測(cè)度函數(shù);根據(jù)評(píng)估對(duì)象x i(i=1,2,…,n)的第j個(gè)指標(biāo)的賦值t j,結(jié)合所建立的單指標(biāo)測(cè)度函數(shù),可以計(jì)算在不同等級(jí)C k下的屬性測(cè)度(稱為單指標(biāo)測(cè)度μxj k);將單指標(biāo)測(cè)度μxj k與對(duì)應(yīng)的指標(biāo)權(quán)重相乘,計(jì)算確定綜合屬性測(cè)度,并基于屬性識(shí)別準(zhǔn)則綜合評(píng)估風(fēng)險(xiǎn)等級(jí)。屬性識(shí)別模型具有計(jì)算簡(jiǎn)便、可操作性強(qiáng)等突出優(yōu)勢(shì),目前,已經(jīng)廣泛應(yīng)用于圍巖穩(wěn)定性、煤與瓦斯突出風(fēng)險(xiǎn)評(píng)估等方面。
2.2.1 單指標(biāo)屬性測(cè)度分析
根據(jù)指標(biāo)屬性不同,建立不同的單指標(biāo)測(cè)度函數(shù),基于指標(biāo)x ij的測(cè)量值t j,計(jì)算單指標(biāo)屬性測(cè)度μxj k,單指標(biāo)等級(jí)劃分判據(jù)見表3。
表3 單指標(biāo)分級(jí)判據(jù)
對(duì)于表3中a j k,滿足a j0＜a j1＜…＜a j k,或者a j0＞a j1＞…＞a jk,令:
式中,k=1,2,…,p。
當(dāng)a j0＜a j1＜…＜a j k時(shí),構(gòu)造單指標(biāo)屬性測(cè)度函數(shù)μxjk(t),見式(6)～式(8)。
式中,k=1,2,…,p-1;j=1,2,…,m。
當(dāng)a j0＞a j1＞…＞a jk時(shí),構(gòu)造單指標(biāo)屬性測(cè)度函數(shù)μxj k(t),見式(9)～式(11)。
2.2.2 多指標(biāo)綜合屬性測(cè)度
考慮到待評(píng)對(duì)象往往受多方面因素共同影響,需要結(jié)合指標(biāo)權(quán)重,計(jì)算綜合屬性測(cè)度μxk,計(jì)算公式為:
式中,w j為指標(biāo)權(quán)重。
2.2.3 屬性識(shí)別分析
在置信度識(shí)別準(zhǔn)則下,結(jié)合綜合屬性測(cè)度μxk計(jì)算結(jié)果,綜合分析待評(píng)對(duì)象狀況。計(jì)算公式為:
式中,λ為置信度,取值區(qū)間為[0,1],一般取0.7;k0為對(duì)應(yīng)的評(píng)價(jià)等級(jí)。
以廣東省云浮市某尾礦庫(kù)主壩為例,采用本文所建立的G1法-屬性識(shí)別模型評(píng)價(jià)尾礦壩的穩(wěn)定性。該尾礦庫(kù)是1986年建成,1988年開始投入使用,由主壩、1#副壩、2#副壩、3#副壩和4#副壩構(gòu)成,屬于傍山型尾礦庫(kù)。2010年,為滿足生產(chǎn)實(shí)際需要,對(duì)尾礦庫(kù)后期運(yùn)行進(jìn)行重新設(shè)計(jì),設(shè)計(jì)最終堆積壩壩頂標(biāo)高為85.0 m,主壩總壩高為57 m,總庫(kù)容為1858.6×104m3,有效庫(kù)容為1604.7×104m3,根據(jù)《尾礦設(shè)施設(shè)計(jì)規(guī)范》(GB 50863—2013),屬三等尾礦庫(kù)。為貫徹國(guó)家節(jié)能環(huán)保減排理念和最大限度回收利用礦物資源的要求,該企業(yè)決定回收尾礦庫(kù)尾礦中的有價(jià)元素。目前,該尾礦庫(kù)仍在回收利用中。該庫(kù)區(qū)不存在巖溶、采空區(qū)、泥石流等不良地質(zhì)。
為做好此次評(píng)估,組織了5位長(zhǎng)期從事尾礦壩體穩(wěn)定性評(píng)估相關(guān)研究工作、具有高級(jí)工程師職稱的專家組成臨時(shí)評(píng)估專家組,對(duì)各指標(biāo)進(jìn)行賦值。其中,對(duì)于某些定量指標(biāo),根據(jù)設(shè)計(jì)圖紙、現(xiàn)場(chǎng)測(cè)量情況確定;對(duì)于某些定性指標(biāo),則由專家根據(jù)指標(biāo)分級(jí)標(biāo)準(zhǔn),經(jīng)過充分討論形成一致性意見后確定。
3.1 基于G1法計(jì)算指標(biāo)權(quán)重
按照權(quán)重計(jì)算原則,需要先確定不同降雨條件下工程、管理、地質(zhì)準(zhǔn)則層指標(biāo)權(quán)重,然后確定4個(gè)準(zhǔn)則層下的各指標(biāo)層權(quán)重,最后,將準(zhǔn)則層和目標(biāo)層指標(biāo)權(quán)重相乘,得到最終目標(biāo)權(quán)重。以降雨工況小雨為例,專家給出的準(zhǔn)則層指標(biāo)序列為工程＞管理＞地質(zhì)＞降雨,對(duì)應(yīng)的相對(duì)重要性比值分別為r1=1.4、r2=1.6、r3=1.1,根據(jù)式(1)至式(3),得到α4=0.159,α3=0.174,α2=0.278,α1=0.389。同理,不同降雨條件下各準(zhǔn)則層及準(zhǔn)則層下各指標(biāo)權(quán)重具體計(jì)算結(jié)果見表4。最終各指標(biāo)目標(biāo)權(quán)重分配情況見表5。
表4 各層次指標(biāo)權(quán)重計(jì)算結(jié)果
表5 不同降雨條件下指標(biāo)權(quán)重分配結(jié)果
為更加直觀反映不同降雨條件下各指標(biāo)權(quán)重的動(dòng)態(tài)變化情況,繪制了如圖1所示的權(quán)重變化趨勢(shì)圖。由圖1可知,隨著降雨強(qiáng)度的增大,各指標(biāo)權(quán)重呈現(xiàn)出不同程度的非線性起伏變化特征,即各指標(biāo)對(duì)尾礦壩穩(wěn)定性的影響程度隨著降雨條件的變化而改變,采用動(dòng)態(tài)權(quán)重進(jìn)行屬性識(shí)別綜合評(píng)價(jià),其結(jié)果更符合實(shí)際。
圖1 各降雨條件下指標(biāo)權(quán)重動(dòng)態(tài)變化情況
3.2 屬性識(shí)別計(jì)算
根據(jù)式(4)至式(11),構(gòu)建單指標(biāo)屬性測(cè)度函數(shù)。將專家指標(biāo)賦值結(jié)果代入不同等級(jí)下的單指標(biāo)屬性測(cè)度函數(shù)中,計(jì)算各等級(jí)下的單指標(biāo)屬性測(cè)度值,根據(jù)式(12),基于置信度識(shí)別原理,計(jì)算不同降雨條件下的綜合屬性測(cè)度,綜合分析評(píng)估尾礦壩體穩(wěn)定性狀態(tài),具體計(jì)算結(jié)果見表6。
由表6可知,隨著降雨量增大,尾礦壩危險(xiǎn)性也隨著增加。主要原因是隨著降雨量的變化,對(duì)應(yīng)的危險(xiǎn)等級(jí)升高,且隨著降雨量的增加,其他風(fēng)險(xiǎn)因素對(duì)尾礦壩穩(wěn)定性影響程度也發(fā)生了變化,危險(xiǎn)等級(jí)最終表現(xiàn)為逐步升級(jí);而可拓學(xué)等其他文獻(xiàn)中,不同降雨工況下各指標(biāo)權(quán)重保持恒定,難以客觀反映指標(biāo)對(duì)穩(wěn)定性狀態(tài)影響程度,本文評(píng)價(jià)結(jié)果與專家討論意見基本一致,表明所建立的綜合評(píng)價(jià)模型評(píng)估結(jié)果更為精準(zhǔn),也符合實(shí)際情況。
表6 綜合評(píng)價(jià)結(jié)果
(1)不同降雨條件對(duì)壩體穩(wěn)定性的影響具有差異性,為此引入了G1法,分別計(jì)算了5種降雨條件下的動(dòng)態(tài)指標(biāo)權(quán)重,能夠有效反映降雨條件下各指標(biāo)對(duì)尾礦壩穩(wěn)定性的影響,可以提高綜合評(píng)估結(jié)果的合理性。
(2)考慮到降雨條件下尾礦壩穩(wěn)定性評(píng)估存在復(fù)雜性,引入了屬性識(shí)別理論,建立綜合評(píng)價(jià)模型,以廣東省云浮市某尾礦庫(kù)為研究背景,分析評(píng)估了5種降雨條件下尾礦壩穩(wěn)定性等級(jí),結(jié)果顯示數(shù)學(xué)評(píng)價(jià)模型與專家實(shí)際調(diào)查討論結(jié)論相一致,表明所建立的基于G1法-屬性識(shí)別綜合模型在尾礦庫(kù)壩體穩(wěn)定性評(píng)價(jià)方面具有一定適應(yīng)性,可以為評(píng)估工作提供一條新思路。
(3)降雨對(duì)尾礦壩穩(wěn)定性影響的相關(guān)研究成果和理論還不完善,為更加真實(shí)全面客觀分析尾礦庫(kù)壩體穩(wěn)定性狀態(tài),需要進(jìn)一步完善指標(biāo)體系,優(yōu)化分級(jí)標(biāo)準(zhǔn),提高模型適應(yīng)性。
猜你喜歡尾礦庫(kù)測(cè)度尾礦尾礦庫(kù)空間信息提取與形變監(jiān)測(cè)應(yīng)用衛(wèi)星應(yīng)用(2023年1期)2023-02-21某鐵礦山尾礦庫(kù)回采工藝設(shè)計(jì)實(shí)例昆鋼科技(2022年2期)2022-07-08三個(gè)數(shù)字集生成的自相似測(cè)度的乘積譜數(shù)學(xué)物理學(xué)報(bào)(2022年3期)2022-05-25R1上莫朗測(cè)度關(guān)于幾何平均誤差的最優(yōu)Vornoi分劃數(shù)學(xué)物理學(xué)報(bào)(2022年2期)2022-04-26基于視覺識(shí)別的浮選尾礦在線測(cè)灰儀的應(yīng)用選煤技術(shù)(2021年3期)2021-10-13長(zhǎng)期運(yùn)行尾礦庫(kù)的排滲系統(tǒng)滲透特性的差異化反演分析有色金屬（礦山部分）(2021年4期)2021-08-30《固體礦產(chǎn)尾礦分類》等3項(xiàng)行業(yè)標(biāo)準(zhǔn)于2021年6月1日起實(shí)施河北地質(zhì)(2021年2期)2021-08-21非等熵Chaplygin氣體測(cè)度值解存在性數(shù)學(xué)物理學(xué)報(bào)(2020年4期)2020-09-07Cookie-Cutter集上的Gibbs測(cè)度數(shù)學(xué)年刊A輯(中文版)(2020年2期)2020-07-25某金礦重選尾礦回收金、銅的工藝研究礦產(chǎn)綜合利用(2020年1期)2020-07-24