尾礦壩風險評估中指標相關(guān)性權(quán)重的確定

史冬梅 楊風暴 王肖霞

(中北大學信息與通信工程學院，山西 太原 030051)

尾礦壩風險評估中指標相關(guān)性權(quán)重的確定

史冬梅 楊風暴 王肖霞

(中北大學信息與通信工程學院，山西 太原 030051)

針對現(xiàn)有尾礦壩風險評估方法未考慮監(jiān)測指標間的相關(guān)性造成的評估過程中計算重復(fù)的問題,考慮監(jiān)測指標間的相關(guān)性并將其轉(zhuǎn)化為權(quán)重，使相關(guān)性的指標中占主導(dǎo)地位的指標權(quán)具有較小的權(quán)值，反之亦然，從而降低重復(fù)計算問題。采用綜合賦權(quán)法確定權(quán)重的方法為：根據(jù)監(jiān)測指標間相互影響程度不同的特點，建立影響矩陣，由于該矩陣所需專家信息具有模糊性，故采用梯形模糊數(shù)描述專家信息，計算主觀權(quán)重；根據(jù)各指標監(jiān)測數(shù)據(jù)具有隨機、模糊等不確定性的特點，利用可能性分布對其進行表征，進而采用相關(guān)系數(shù)賦權(quán)法計算得到各指標客觀權(quán)重；通過均值法綜合主、客觀權(quán)重。實例表明，考慮了監(jiān)測指標相關(guān)性的權(quán)重確定法能減少風險評估中指標的重復(fù)計算，綜合主客觀權(quán)重，為尾礦壩各指標賦權(quán)提供了理論依據(jù)。

尾礦壩 風險評估 相關(guān)性 影響矩陣 可能性相關(guān)系數(shù)賦權(quán)法

尾礦壩的風險評估是目前尾礦研究方面的重點，現(xiàn)有風險評估模型使用的前提條件是指標間應(yīng)具有獨立性，是非相關(guān)的[1]，即指標之間沒有相互影響的關(guān)系，一個指標的變化不會導(dǎo)致另一個指標的變化。然而實際上，指標獨立是難以實現(xiàn)的，如庫水位升高將會導(dǎo)致干灘長度變短以及浸潤線的升高、浸潤線升高導(dǎo)致壩體不穩(wěn)定，壩體位移量增大等[2-3]。由于指標相關(guān)性大，會帶來評估過程中重復(fù)計算的問題。因此，需要盡可能地降低指標相關(guān)性對風險評估結(jié)果的影響。2個指標相關(guān)性越大，則一個指標可以被另一個指標解釋得越多，在指標體系中作用也越小，即指標的權(quán)重越小。因此，本研究通過考慮監(jiān)測指標相關(guān)性的權(quán)值確定方法解決評估過程中的重復(fù)計算問題。

目前尾礦壩風險評估中權(quán)重的確定方法主要有主觀賦權(quán)法中的層次分析法[4]、Delphi法[5]等、客觀賦權(quán)法中的離差最大化法[6]、熵權(quán)法[7]等，以及主客觀結(jié)合的綜合賦權(quán)法[8]。主觀賦權(quán)法主觀隨意性大，客觀賦權(quán)法所得結(jié)果可能與事實相悖，綜合賦權(quán)法能夠克服這二者的缺點，因此得到了廣泛的應(yīng)用。但是，這些權(quán)重確定方法既沒有考慮指標相關(guān)性的影響，也沒有考慮監(jiān)測信息的隨機、模糊等不確定性對權(quán)重結(jié)果的影響。

因此，利用能夠考慮專家信息模糊性以及指標間相互影響程度不同的梯形模糊數(shù)影響矩陣來計算主觀權(quán)重；在對尾礦壩各指標監(jiān)測信息分析的基礎(chǔ)上，利用可以表征多種不確定性的可能性分布描述各指標信息，進而利用相關(guān)系數(shù)賦權(quán)法[9]計算客觀權(quán)重；綜合主客觀權(quán)重，得到考慮相關(guān)性的監(jiān)測指標權(quán)重，為尾礦壩風險評估的進一步研究提供參考信息。

1 監(jiān)測指標相關(guān)性權(quán)重確定方法

1.1 主觀權(quán)重的確定

確定主觀權(quán)重最常用的方法為層次分析法，該方法構(gòu)建的判斷矩陣是用來衡量指標重要性程度的，但是并沒有考慮指標間相關(guān)性的影響。為解決這個問題，本研究構(gòu)建了指標間的影響矩陣[10]，并且考慮到專家信息的模糊性，引入梯形模糊數(shù)表征專家信息，使得構(gòu)建的影響矩陣更合理，更準確。利用該矩陣計算尾礦壩各監(jiān)測指標主觀權(quán)重的具體過程如下。

1.1.1 構(gòu)建影響矩陣

尾礦壩4個監(jiān)測指標：庫水位(u1)、浸潤線(u2)、壩體位移(u3)、干灘長度(u4)的影響矩陣為

其中，aij表示指標ui對指標uj的影響程度。值得注意的是，aij不是指標的相關(guān)度。因此，一般情況下，aij≠aji，這也是不同于層次分析法的地方。該矩陣能夠考慮到2個指標相互影響程度的不同。

上述矩陣請尾礦壩相關(guān)專家給出，通常采用表1所示標度對指標間影響程度定量衡量。但是專家打分本質(zhì)是一個模糊過程，采用確定數(shù)值會對分析結(jié)果的精確性產(chǎn)生影響。因此，以表1所示標度法為基礎(chǔ)，利用梯形模糊數(shù)表征專家信息。

表1 標度法及其含義

設(shè)專家給出的指標ui對uj影響程度的梯形模糊數(shù)為

aij=(lij,mij,nij,sij),

lij≤mij≤nij≤sij,

其隸屬函數(shù)μaij(x)為

(1)

該梯形模糊數(shù)的期望值

E(aij)=bEL(aij)+(1-b)EU(aij),

(2)

其中，

為梯形模糊數(shù)的左期望值；

為梯形模糊數(shù)的右期望值;b為樂觀系數(shù)，表明專家對所給結(jié)果的態(tài)度。如果0≤b≤0.5表明專家持偏向悲觀的態(tài)度；如果0.5

1.1.2 確定指標權(quán)重

基于上述所得影響矩陣，計算尾礦壩各監(jiān)測指標主觀權(quán)重過程如下。

(1)計算指標ui的模糊影響值

(2)計算ui的期望值E(ui)，通常取b=0.5。

(4)

1.2 客觀權(quán)重的確定

由于尾礦壩壩體結(jié)構(gòu)的復(fù)雜性和動態(tài)性，傳感器本身監(jiān)測精度與范圍的局限性，自然環(huán)境影響的強干擾性等，使得尾礦壩的監(jiān)測信息具有模糊、隨機、不完全等不確定性，通常的相關(guān)性權(quán)重計算方法并沒有考慮這些不確定性。由于可能性分布能夠處理多種不確定性，故利用可能性分布表征監(jiān)測信息，進而利用相關(guān)系數(shù)賦權(quán)法確定客觀權(quán)重。具體過程如下。

(1)利用樣本監(jiān)測數(shù)據(jù)的可能性分布構(gòu)造法，將尾礦壩各指標監(jiān)測信息轉(zhuǎn)化為三角可能性分布。

(5)

其中，π(x)為指標uj的可能性分布，m為監(jiān)測數(shù)據(jù)x的均值，m1和m2分別為以均值為界限的2組監(jiān)測數(shù)據(jù)的均值。

(2)計算指標間的相關(guān)系數(shù)

(6)

由此可構(gòu)建出各監(jiān)測指標的相關(guān)系數(shù)矩陣

(7)

其中，rij(i,j=1,2,…n)為ui與uj的相關(guān)系數(shù)，rij=1表示本身相關(guān)程度最大。該矩陣為一個對稱矩陣。

(3)計算指標uj與其他指標相關(guān)程度的均值

(8)

(4)如果一個指標與另一指標相關(guān)性越大，則其在指標體系中所占比重越小，即權(quán)重越小。此處認為相關(guān)性與權(quán)重有如下倒數(shù)關(guān)系

qj=1/δj,

(9)

(5)計算指標uj的客觀權(quán)重(歸一化處理)

(10)

計算出尾礦壩各監(jiān)測指標的客觀權(quán)重W″

W″=(W″u1,W″u2,W″u3,W″u4,W″u5).

1.3 綜合賦權(quán)法

將基于影響矩陣獲得的主觀權(quán)重與可能性相關(guān)系數(shù)賦權(quán)法所得的客觀權(quán)重利用式(10)融合計算，得到權(quán)重集合

W=(Wu1,Wu2,Wu3,Wu4,Wu5)

即為尾礦壩各監(jiān)測指標的權(quán)重?？紤]專家經(jīng)驗，同時也考慮客觀因素，能夠比較全面地反映監(jiān)測指標的相對重要程度的綜合權(quán)重為

W=βW′+(1-β)W″,

(11)

其中，β為主客觀權(quán)重的比例因子，反映人們對權(quán)重結(jié)果偏主觀一些還是客觀一些，0<β≤1。通常認為主客觀結(jié)果占同等重要的地位，故取β=1/2。

2 實例分析

以山西省某尾礦壩為例來說明庫水位、浸潤線、壩體位移、干灘長度4個指標的權(quán)重確定過程。

專家給出的尾礦壩信息影響矩陣為

根據(jù)該影響矩陣計算出主觀權(quán)重為

W″u1=0.316 0,W″u2=0.255 8,

W″u3=0.189 2,W″u4=0.239 1.

將尾礦壩各指標監(jiān)測信息利用式(5)轉(zhuǎn)化為對應(yīng)的可能性分布，如圖1所示。

根據(jù)上述可能性分布，利用式(6)計算指標間相關(guān)系數(shù)，得相關(guān)系數(shù)為

通過計算可得客觀權(quán)重：

W″u1=0.280 0,W″u2=0.248 0,

W″u3=0.235 4,W″u4=0.236 4.

熵權(quán)法所得客觀權(quán)重：

wu1=0.251 2,wu2=0.251 2,

wu3=0.247 7,wu4=0.250 0.

將2組客觀權(quán)重對比分析可知，可能性相關(guān)系數(shù)賦權(quán)法相比于熵權(quán)法，能夠?qū)⒕哂邢嚓P(guān)性的指標通過權(quán)重的方式進行調(diào)整，將在相關(guān)性中占主導(dǎo)地位的指標賦予較小權(quán)值，從而更加合理地確定指標權(quán)重，消除由于指標重復(fù)對風險評估結(jié)果的影響。

利用式(10)，得到綜合主客觀權(quán)重：

Wu1=0.298 1,Wu2=0.251 9,

Wu3=0.212 3,Wu4=0.237 8.

圖1 各指標監(jiān)測信息可能性分布

3 結(jié) 語

本研究在確定尾礦壩監(jiān)測指標權(quán)重時，考慮了指標間的相關(guān)性問題，所采用的方法為主客觀結(jié)合的綜合賦權(quán)法。主觀賦權(quán)法通過梯形模糊數(shù)影響矩陣得到，不僅考慮了每對監(jiān)測指標的相互影響程度的不同，還用梯形模糊數(shù)表征了專家信息的模糊性；客觀賦權(quán)法利用可能性相關(guān)系數(shù)賦權(quán)法計算得到，該方法利用可能性分布表征監(jiān)測信息多種不確定性，同時將指標間的相關(guān)性以定量的方式衡量，轉(zhuǎn)化為權(quán)重，從而解決風險評估中指標重復(fù)計算的問題。

[1] Hai Yanju，Gui Qinggao，Qiu Xiangli，et al.Analysis of safety and risk factors about the tailings dam in the metal mines[J].Advanced Materials Research，2012，594/597:299-302.

[2] 戴 繪,秦衛(wèi)星,張岳安.強降雨作用下山谷型尾礦壩浸潤線演化規(guī)律[J].金屬礦山，2013(7):149-152. Dai Hui，Qin Weixing，Zhang Yuean.Evolution law of phreatic line in valley-type tailings dam under heavy rainfall[J].Metal Mine，2013(7):149-152.

[3] 劉海明,曹 凈,楊春和.國內(nèi)外尾礦壩事故致災(zāi)因素分析[J].金屬礦山，2013(2):126-129. Liu Haiming，Cao Jing，Yang Chunhe.Analysis of factors for tailings dams accidents causing disaster at home and abroad[J].Metal Mine，2013(2):126-129.

[4] Yan Shanyu，Dong Yang，Pan Ke.Application of multi-level extensible method in risk assessment of tailings dams failure[J].Journal of Northeastern University，2013，34:80-83.

[5] Mei Guodong，Wu Zongzhi.Study on social risk evaluation index system for tailings pond dam-break based on the vulnerability theory[J].Advanced Materials Research:Advances in Industrial and Civil Engineering，2012,594/597:2301-2308.

[6] 張影秋,吳瑞明,張 揚.基于離差最大化與灰色關(guān)聯(lián)分析的危險源辨識方法[J].科學技術(shù)與工程，2010，10(25):6358-6362. Zhang Yingqiu，Wu Ruiming，Zhang Yang.The method of hazard identification based on maximizing deviations and grey relational analysis[J].Science Technology and Engineering，2010，10(25):6358-6362.

[7] 袁 維，白 冰，李小春，等.透鏡體對尾礦壩安全性影響的參數(shù)敏感性分析[J].中南大學學報：自然科學版，2013，44(3):1174-1183. Yuan Wei，Bai Bing，Li Xiaochun，et al.Parameters sensitivity analysis of lenticles' impacting on tailings dam safety[J].Journal of Central South University：Science and Technology Edition，2013，44(3):1174-1183.

[8] 王肖霞，楊風暴,史冬梅.可能性熵權(quán)在尾礦壩風險評估中的應(yīng)用[J].科學技術(shù)與工程，2013,13(20):5766-5771. Wang Xiaoxia，Yang Fengbao，Shi Dongmei.Application of possibility entropy weight in security risk assessment of tailing dam[J].Science Technology and Engineering，2013,13(20):5766-5771.

[9] 呂明捷，杜 云，榮 超，等.基于相關(guān)系數(shù)定權(quán)的集對分析法在湖泊富營養(yǎng)化評價中的應(yīng)用[J].南水北調(diào)與水利科技，2011，9(1):96-98. Lu Mingjie，Du Yun，Rong Chao，et al.Set pair analysis method in lake's eutrophication assessment based on correlation weight method[J].South to North Water Transfers and Water Science & Technology，2011,9(1):96-98.

(責任編輯 徐志宏)

Determining of Correlation Weights of Index in the Tailing Dam Risk Assessment

Shi Dongmei Yang Fengbao Wang Xiaoxia

(College of Information and Communications Engineering，North University of China，Taiyuan 030051，China)

To solve the problem that the present risk assessment method of tailings dam does not consider the correlation between monitoring indexes resulting in double calculation during the evaluation process，the correlation between the monitoring indexes is considered and converted them into weights，so that the dominant index had the smaller weight，and vice versa.So the problem of double counting can be lowered.An integrated weighting method was used to determine the weights:according to the characteristics that monitoring indexes has different influence on each other，an influencing matrix is established.As the experts information required in the matrix is fuzzy，the trapezoidal fuzzy number is used to describe experts information and to calculate the subjective weights；With features of fuzziness and randomness，the monitoring data are represented by possibility distributions，then the correlation coefficient weighting method was used to calculate the objective weight of every index；Subjective weight and objective weight are integrated by the mean method.Example results show that:the weight calculation method that considers the coefficient of monitoring indexes can reduce duplication in risk assessment，and integrate subjective and objective weights.This method provides the theoretical basis for weighting of each index for tailings dam.

Tailings dam，Risk assessment，Correlation，Influence matrix，Possibility of correlation coefficient weighting method

2014-09-22

山西省高等學校留學回國人員科研項目(編號：201110)，山西省研究生優(yōu)秀創(chuàng)新項目(編號：20123095)。

史冬梅(1987—)，女，碩士研究生。

TD77

A

1001-1250(2014)-11-143-04