层次分析法(AHP)在尾矿库安全运行分析中的应用.pdf
第 29 卷增刊 岩 土 力 学 Vol.29 Supp. 2008 年 11 月 Rock and Soil Mechanics Nov. 2008 收稿日期2008-07-25 作者简介 王涛, 男, 1980 年生, 硕士研究生, 助理工程师, 主要从事安全技术及工程, 岩土灾害发生机理与防治方面的研究。 E-mail kmlgwangtao 文章编号文章编号1000-7598-2008 增刊-680-07 层次分析法 (层次分析法 (AHP) 在尾矿库安全运行分析中的应用) 在尾矿库安全运行分析中的应用 王 涛 1,侯克鹏1,郭振世2,张成良1 (1. 昆明理工大学 国土资源工程学院,昆明 650093;2. 西安理工大学 水利水电学院,西安 710048) 摘 要摘 要尾矿库是用于堆存金属非金属矿山进行矿石选别后排出尾矿的场所,是一种高势能的人造泥石流源,一旦发生事故 将会严重威胁到库区周边环境和下游居民的生命及财产安全。考虑到我国小型尾矿库数量多、安全度水平低以及尾矿库发生 事故的复杂性、多样性及不确定性等特点,根据影响尾矿库安全运行的因素建立层次结构模型,运用 Delphi 法构造相关因 素的各层次判断矩阵,应用定性与定量相结合的层次分析法进行分析,最后确定出各影响因素的权重并进行总排序。结果表 明,在尾矿库运行期间排洪系统对其整体安全影响最大,其次为筑坝堆存系统、尾矿库管理与维护、回水系统和尾矿输送系 统。该研究结果可以为相关单位制定合理的安全防范措施提供必要的参考价值,从而有针对性的消除各种事故隐患,提高安 全管理水平,确保尾矿库的安全运行。 关关 键键 词词尾矿库;Delphi 法;层次分析法;安全运行;影响因素 中图分类号中图分类号TD 926.4; TV 6; X 913.4 文献标识码文献标识码A Application of analytic hierarchy process to tailings pond safety operation analysis WANG Tao1, HOU Ke-peng1, GUO Zhen-shi2, ZHANG Cheng-liang1 1. College of Land Resources Engineering, Kunming University of Science and Technology, Kunming 650093, China; 2. School of Water Resources and Hydroelectric Power, Xian University of Technology, Xian 710048, China Abstract Tailings pond is a place where is used for depositing tailings of the metal and nonmetal mine after the selecting process. It is one of the high potential energy and human-caused mud-rock flows. It will bring huge harms to the district peripheral environment and downstream inhabitants life and property when tailings pond’s paroxysmal accidents happen. In view of the quantity of small and medium-sized tailings pond, the lower safety management level, furthermore, diversified and uncertain causes during accidents, the paper has made hierarchical structure model according to the influence factors of tailings pond safety operation. Using Delphi to construct correlation factor’s judgment matrix, the combination of qualitative and quantitative analytic hierarchy process AHP is applied; and finally can be obtained in the weight of coefficient and arrayed in order. The analysis shows that the biggest influence factor is the expel-flood system. the second is the embankment system, manage or maintenance, back-water and transportation system. Therefore the related unit may depend on the sequences of the facility to have constituted respectively guard measure, which could supply guiding sense and reference value so as to eliminate each kind of accident hidden dangers and enhance the management level and ensure the tailings pond safety operation. Key words tailings pond; Delphi ; analytic hierarchy process; influencing factors; safety operation 1 引 言 尾矿库是矿业开发中必不可少的重要生产设 施, 其所处位置一般难以避开居民区和人口密集区, 属重大危险源;一旦发生事故将会引起滑坡、泥石 流等重大灾害,对库区下游人民的生命和财产造成 巨危害,对环境造成严重污染。美国克拉克大学公 害评定小组的研究表明,尾矿库事故造成的危害, 在世界 93 种事故、 灾害和辐射隐患中名列第 18 位, 它仅次于核武器爆炸、DDT、神经毒气、核辐射以 及其它 13 种灾害,而比航空失事、火灾等其它 60 种灾害还要严重[1]。 据国家安监总局统计, 截止 2007 年底,全国共有尾矿库 8 540 座;2001~2007 年, 全国共发生尾矿库事故 43 起[2],其中2001 年 3 增刊 王 涛等层次分析法(AHP)在尾矿库安全运行分析中的应用 起,2003 年 2 起,2004 年 3 起,2005 年 9 起,2006 年 12 起,2007 年 14 起。 目前,国外对尾矿库的研究主要集中在坝体安 全和环境两个方面,Kemper T[3]等对尾矿库的污染 和污染对象监测、污水治理方法和理论等方面进行 了研究。国内关于尾矿库的理论研究及工程实践都 是围绕坝体设计、稳定性、地震液化等展开的对尾 矿库的整体安全运行研究相对较少,李明[4]等人利 用 ANSYS 对尾矿库进行建模、有限元分析、计算、 仿真等,然后分析尾矿库坝体稳定性;潘建平、孔 宪京、邹德高[5]对尾矿坝地震液化稳定进行了简化 分析,提出了尾矿坝液化分析方法。 由于尾矿库的运行受多方面因素影响,且各影 响因素对尾矿库的运行有着不同的侧重,同时很多 影响因素都无法直接量化。本文旨在应用层次分析 法将这些影响因素条理化、层次化,将半定量半定 性问题定量化,从而准确、直观地反应出各影响因 素之间的相互关系,并确定出各影响因素相对重要 性的权值及次序,为相关部门进行尾矿库管理提供 参考。 2 尾矿库概论 尾矿库是指在山谷口部或洼地的周围筑坝,将 金属和非金属矿山企业进行矿石选别后的尾矿或其 他工业废渣排入其内进行沉淀的贮存场所,其运行 设施一般由尾矿输送系统、尾矿筑坝堆存系统、尾 矿排洪系统、 尾矿回水系统等组成。 根据地形不同, 常见的有以下三种类型[6]; (1)山谷型尾矿库。山谷型尾矿库是在山谷谷 口处筑坝形成的尾矿库,见图 1。它的特点是初期 坝相对较短,坝体工程量较小;后期管理和维护相 对较容易, 当后期筑坝较高时, 可获得较大的库容, 但汇水面积较大,排水设施工程量大。我国大多数 尾矿库属于该类型。 图图 1 山谷型尾矿库山谷型尾矿库 Fig.1 Tailings pond in the slope of valley (2) 平地型尾矿库。 平地型尾矿库是在平地四周 筑坝围成的尾矿库,见图 2。其特点是初期坝和后 期尾矿堆筑坝工程量大,建设费用高,维护管理较 麻烦;由于四周堆坝,库区面积越来越小。 图图 2 平地型尾矿库平地型尾矿库 Fig.2 Planar tailings pond (3)傍山型尾矿库。傍山型尾矿库是在山坡脚 下依山筑坝所围成的尾矿库,见图 3。其特点是初 期坝相对较长, 初期坝和后期尾矿筑坝工程量较大, 库容较小,管理维护相对复杂。 图图 3 傍山型尾矿库傍山型尾矿库 Fig.3 Tailings pond in the slope of hillside 3 层次分析法(AHP)及其步骤 3.1 Delphi 专家调查法专家调查法 Delphi 专家调查法, 亦称德尔菲法,是 20 世纪 40 年代美国兰德公司提出的征集专家们意见据以 判断决策的一种系统分析方法,其本质上是一种反 馈匿名函询法。这种方法简单、节省费用,能把理 论知识和实践经验的各方面专家对同一问题的意见 集中起来,适用于研究资料少、未知因素多、主要 靠主观判断和粗略估计来确定的问题,是较多地用 于长期预测和动态预测的一种重要的预测方法[8]。 经典德尔菲法是对所要预测的问题征得专家的 意见后,进行整理、归纳、统计、再匿名反馈给各 专家,然后再次征求意见,再次归纳、统计、再反 馈给各专家,直至得到稳定的意见。征求过程一般 包含四轮的征询调查,且在调查过程中包含轮与轮 间的反馈,在第四轮调查结束后,组织者依然要将 回收的调查表进行汇总整理、统计分析与预测,并 寻找出收敛程度较高的专家意见,最后将最终的调 查结果应用层次分析法进行计算。 681 岩 土 力 学 2008 年 3.2 层次分析法基本原理及计算步骤层次分析法基本原理及计算步骤 层次分析法(Analytical Hierarchy Process,简 称 AHP)是美国运筹学家、匹兹堡大学 T. L. Saaty 教授于 20 世纪 70 年代提出的一种定性和定量相结 合的系统化、 层次化分析方法。 应用 AHP 解决问题 的思路是首先,把要解决的问题分层系列化,即 根据问题的性质和要达到的目标,将问题分解为不 同的组成因素,按照因素之间的相互影响和隶属关 系将其分层类聚组合,形成一个递阶的、有序的层 次结构模型,然后,对模型中每一层次因素的相对 重要性,依据人们对客观现实的判断(通过 Delphi 专家调查法)给予定量表示,再利用数学方法确定 每一层次全部因素相对重要性次序的权值,最后, 通过综合计算各层因素相对重要性的权值,以此作 为评价和选择方案的依据。应用层次分析法解决问 题一般包括以下 4 个步骤[9] (1)建立多级递阶层次结构模型。。根据对问题 的了解和初步分析,将分析系统涉及的各要素按性 质分层排列,最简单的层次结构可分为 3 级,第 1 级是目标层(最高层) ,该级是系统要达到的目标, 一般情况下只有一个目标,如果有多个分目标,可 以在下一级设立一个分目标层。 第 2 级是准则层 (中 间层) , 该级列出衡量达到目标的各项准则, 如果某 些准则还需具体化,即做进一步解释说明,则可在 下一级再设立一个准则层。第 3 级是措施层(最低 层) , 该级排列了各种可能采取的方案或措施。 不同 层次的各要素间的关系用连线表示,如果要素间有 连线,表示二者相关,否则表示不相关。根据层次 分析法中的标度规定,每一层中的指标数一般不超 过 9 个,元素过多会给两两比较带来困难。 (2)构造判断矩阵。。任何系统分析都是以一定 的信息为基础, AHP 的信息基础主要是人们对每一 层次各要素的相对重要性给出的数值判断,所以, 判断矩阵是 AHP 的基本信息, 也是进行相对重要度 计算,进行层次单排序的依据。其方法是以上一级 的某一要素 A 作为准则,对本级的要素进行两两比 较来确定判断矩阵元素的,见表 1。 在该矩阵中,按照左上角的准则 A,将左侧列 元素中的 B1与最上一行的 B1,B2,B3,BjBn进 行比较;然后再用左侧列元素中的 B2与 B1,B2, B3BjBn进行比较,如此类推。在比较时,具体地 使用数字标度来代表一个元素针对准则超越另一个 元素的相对重要性,数字标度定义并解释了递阶层 次中每层元素针对上层准则而进行成对比较时,两 个指标哪一个更重要,重要多少,并按“1-9”标度 对重要性程度赋值,表 2 中列出了“1-9”标度的含 义。任何判断矩阵都应满足 1 ii b ; 1/ ijji bb ,1,2,..., i jn (1) 表表 1 判断矩阵判断矩阵 Table 1 Judgment matrix A B1 B2 B3 Bj Bn B1 1 b12 b13 b1j b1n B2 b21 1 b23 b2j b2n B3 b31 b32 1 b3j b3n Bi bi1 bi2 bi3 bij bin Bn bn1 bn2 bn3 bnj 1 表表 2 判断矩阵标度与含义判断矩阵标度与含义 Table 2 Judgment matrix scales and meanings 标度 含义 1 Bi比 Bj同样重要 3 Bi比 Bj稍重要 5 Bi比 Bj重要 7 Bi比 Bj重要得多 9 Bi比 Bj极为重要 2,4,6,8 Bi比 Bj相比, 重要性程度处于上述 相应两个数之间 1/3 Bi比 Bj稍次要 1/5 Bi比 Bj次要 1/7 Bi比 Bj次要得多 1/9 Bi比 Bj极为次要 (3)层次单排序及一致性检验。。在建立了判断 矩阵后,要根据判断矩阵计算本级要素相对于上一 级某一些要素来讲,本级与之有联系的要素之间相 对重要性次序的权值,即根据判断矩阵求出其最大 特征根 max λ,然后再求出所对应的特征向量 W。方 程如下 max BλWW (2) 式中W的分量( 1 W, 2 W,..., n W)就是对应于 n个要素的相对重要度,即权重系数。权重的计算 方法主要有和积法、方根法、特征根方法、对数 最小二乘法和最小二乘法等。本文采用“和积法” 进行权重计算,计算步骤如下[9] ① 将判断矩阵每一列正规化。即 1 , , 1,2,..., ij ij n kj k b bi jn b ∑ (3) ② 每一列经正规化后的判断矩阵按行相加。 即 682 增刊 王 涛等层次分析法(AHP)在尾矿库安全运行分析中的应用 1 , ,1,2,..., n iij j Wbi jn ∑ (4) ③ 对向量 122 ,,..., T W WW⎤⎡ ⎣⎦ W正规化。即 1 , 1,2,..., i n j j W Win W ∑ (5) 所得到的] T 12 , , ..., n WWW⎡ ⎣ W即为所求特征向量。 ④ 计算判断矩阵最大特征根 max λ。即 max 1 n i i i AW nW λ ∑ (6) 式中iAWAWi表示向量的第 个分量。 表示向量 AW的第i个分量。 为检验判断矩阵的一致性如何,需要计算它的 一致性指标CI,定义 max 1 n CI n λ− − (7) 显然, 当判断矩阵具有完全一致性时,0CI 。 为了检验判断矩阵是否具有满足的一致性,需要将 CI与平均随机一致性指标RI(如表3所示)进行 比较,记为CR。 CI CR RI (8) 当CR<0.10时,判断矩阵具有满意的一致性, 否则就需要对判断矩阵进行调整。 表表 3 “1-9”阶矩阵阶矩阵 RI 值值 Table 3 The values of RI 阶数/ n 1 2 3 4 5 6 7 8 9 RI 0.00 0.00 0.58 0.90 1.12 1.24 1.32 1.411.45 (4)层次总排序及一致性检验。利用同一层次 中所有层次单排序的结果, 计算针对上一层次而言, 本层次所有要素重要性的权值。层次总排序需要从 上到下逐层顺序进行,对于最高层下面的第二层, 其层次单排序即为总排序。为评价层次总排序计算 结果的一致性,需要计算与单排序类似的检验量。 CI为层次总排序一致性指标;RI为层次总排序平 均随机一致性指标;CR为层次总排序随机一致性 比例。各表达式分别为 1 n ii i CIBCI ∑ (9) 式中 i CI为与 i B对应的B层次中判断矩阵的一致 性指标;Bi为判断距阵A-B(即第一层与第二层目 标组成的判断矩阵)计算所得的向量经正规化后的 特征向量W的值。 1 n ii i RIB RI ∑ (10) 式中 i RI为与 i B对应的B层次中判断矩阵的平均 随机一致性指标。 CI CR RI (11) 同样,当CR<0.10时,认为层次总排序的计算结 果具有满意的一致性。 4 层次分析法在尾矿库安全运行分 析中的应用 4.1 工程实例工程实例 某尾矿库位于秦岭南坡山区,三面靠山,一面 筑坝,属山谷型尾矿库,库区无滑坡、泥石流、管 涌等不良地质现象,岸坡稳定,水土保持较好。 该尾矿库初期坝为堆石透水坝,坝高40.5 m, 坝顶长115 m,坝顶宽4.0 m,坝底宽157.50 m,上 游坡比11.7,下游坡比12.0。后期堆筑坝采用上 游法堆筑(用推土机或人工将滩面尾矿堆筑而成) , 设计外坡比15.0,分期逐级堆筑子坝,每期高度 3.0 m,台阶宽度8.0~11.0 m,顶宽3.5 m,外坡比为 11.5,在子坝顶布设尾矿管道进行分散放矿,尾矿 分级沉积,形成沉积滩。 该尾矿库堆存的尾矿是经五级泵站扬送至坝顶 进行分散放矿,沿途泵站建有尾矿事故池及事故返 砂泵、尾矿输送自流管桥、靠近尾矿库建有新、老 放矿流槽,流槽接坝顶的放矿主管;回水系统原采 用浮船泵站进行扬水,现采用隧洞自流回水,自流 回水系统由进水塔、 斜槽、 输水隧洞和输水管组成; 尾矿库的老排洪系统设于库的左岸,由两座排洪井 及排洪斜槽、排洪涵管、排洪隧洞组成;随着库内 水位的逐渐升高,新修建的排洪系统由竖井和排洪 隧洞组成。 4.1 建立递阶层次模型建立递阶层次模型 根据尾矿库的实际工况,首先将尾矿库安全运 行作为第1级层次 (即最高层) , 然后将所涉及的相 关影响因素进行分层,第1级层次包括5个因素, 即{} 12345 , , , ABBB BB;第2级层次包括16个因 素, 即{ 111121314 , ,, ,Bbbbb} 15 b,{} 2212223 , , Bbbb, {} 3313233 , , Bbbb, 4 B{} 4142 , bb,{} 5515253 ,,Bbbb, 如图4所示。 683 岩 土 力 学 2008 年 图图 4 尾矿库安全运行层次结构模型尾矿库安全运行层次结构模型 Fig.4 Hierarchical structure model for the tailings pond safety operation 4.3 构造判断矩阵、计算各指标值构造判断矩阵、计算各指标值 按照图4 表示的层次结构模型,通过两两比较 下一层次诸因素对上一层次某元素的相对重要性, 并赋予一定的分值(采用1-9阶标度法赋值),建 立上下层次之间的隶属关系构造判断矩阵(见表 4~9),并用“和积法”计算相应的权重,然后做 一致性检验。 表表 4 判断矩阵判断矩阵 A-B 数值表数值表 Table 4 Judgment matrix for A-B A B1 B2 B3 B4 B5W 指标 B1 1 1/2 3 3 20.25 B2 2 1 6 3 2 0.33 max 5.13λ B3 1/3 1/6 1 1 1 0.08 CI 0.033 B4 1/3 1/3 1 1 1 0.24 RI 1.12 B5 1/2 1/2 1 1 1 0.10 CR 0.029<0.10 表中各指标值计算过程如下(以下各表中的指标值 方法类同) ① 按照公式(3)计算,得到按列正规化后的 判断矩阵 0.240.200.250.330.29 0.480.400.500.330.29 0.080.070.080.110.14 0.080.130.080.110.14 0.120.200.080.110.14 ⎡⎤ ⎢⎥ ⎢⎥ ⎢⎥ ⎢⎥ ⎢⎥ ⎢⎥ ⎣⎦ (12) ② 按照公式(4) 。将各行相加,得 1 1 2 3 4 5 0.24 0.20 0.25 0.33 0.29 1.31 0.48 0.40 0.50 0.33 0.292.00 0.08 0.070.08 0.11 0.140.48 0.08 0.13 0.08 0.11 0.140.54 0.12 0.20 0.08 0.11 0.140.65 n ij j Wb W W W W ⎫ ⎪ ⎪ ⎪ ⎪ ⎬ ⎪ ⎪ ⎪ ⎪ ⎭ ∑ (13) ③ 按公式(5)将向量W [1.31 2.00 0.48 T 0.54 0.65]正规化,得 1 1.312.000.480.540.654.98 n j j W ∑ (14) 1 1 1 23 35 1.31 0.26 4.98 0.40 0.10 0.11 0.13 n j j W W W WW WW ⎫ ⎪ ⎪ ⎪ ⎬ ⎪ ⎪ ⎪ ⎭ ∑ (15) 则所求特征向量 [] T 0.26 0.40 0.10 0.11 0.13W (16) ④ 按照公式(6)计算判断矩阵的最大特征根 max λ 0.26 1 1332 2 0.40 21632 11 1 1 10.10 36 11 1 1 10.11 33 11 1 1 1 0.13 22 ⎡⎤⎡ ⎤ ⎢⎥⎢ ⎥ ⎢⎥⎢ ⎥ ⎢⎥⎢ ⎥ ⎢⎥⎢ ⎥ ⎢⎥⎢ ⎥ ⎢ ⎥⎢ ⎥ ⎢⎥⎢ ⎥ ⎢⎥⎢ ⎥ ⎢⎥⎢ ⎥ ⎢⎥⎢ ⎥ ⎢⎥⎢ ⎥ ⎣⎦ ⎣⎦ AW (17) 1 23 45 1 1 0.260.403 0.10 2 3 0.1120.131.35 2.11 0.49 0.56 0.67 AW AWAW AWAW ⎫ ⎪ ⎪ ⎪ ⎬ ⎪ ⎪ ⎪⎭ (18) 12 max 1 12 354 345 55 555 n i i i AWAWAW nWWW AWAWAW WWW λ ∑ 尾矿库安全运行 A 尾矿库筑坝堆存系统 B1 尾矿库排洪系统 B2尾矿运输系统 B3尾矿库回水系统 B4尾矿库管理与维护 B5 b11 b12 b13 b14 b15 筑 坝 方 式 及 工 艺 尾 矿 颗 粒 组 成 及 物 理 特 性 堆 积 坝 坡 比 及 侵 润 线 高 度 堆 积 坝 护 坡 及 防 冲 刷 工 程 坝 体 抗 震 能 力 尾 矿 库 调 洪 与 排 洪 能 力 排 水 隧 洞 的 坚 固 性 排 洪 设 施 安 全 性 与 可 靠 性 输 送 管 道 能 力 及 完 好 状 况 浓 缩 池 及 浓 缩 机 的 状 况 砂 泵 站 运 行 情 况 尾 矿 库 回 水 设 施 完 好 情 况 回 水 泵 站 运 行 情 况 尾 矿 库 管 理 机 构 设 置 尾 矿 库 操 作 工 的 技 术 培 训 尾 矿 库 设 施 管 理 与 维 护 b21 b22 b23b31b32 b33b41b42b51 b52 b53 684 增刊 王 涛等层次分析法(AHP)在尾矿库安全运行分析中的应用 1.352.110.49 50.2650.4050.10 0.560.67 5.13 50.1150.13 (19) ⑤ 按照公式(7) 、 (8)进行一致性检验 max 5.135 0.033 151 n CI n λ−− −− (20) 1.12RI 0.033 0.0290.10 1.12 CI CR RI (21) 可见判断矩阵具有满意的一致性。 表表 5 判断矩阵判断矩阵 B1-b 数值表数值表 Table 5 Judgment matrix for B1-b B1 b11 b22 b33 b44 b55W 指标 b11 1 2 1/2 3 1/20.25 b22 1/2 1 1/3 1 1/30.33 max 5.13λ b33 2 3 1 3 1 0.08 CI 0.028 b44 1/3 1 1/3 1 1/20.24 RI 1.12 b55 2 3 1 2 1 0.10 CR 0.025<0.10 表表 6 判断矩阵判断矩阵 B2-b 数值表数值表 Table 6 Judgment matrix for B2-b B2 b21 b22 b33 W 指标 b21 1 7 4 0.70 max 5.11λ b22 1/7 1 1/3 0.10 CI 0.028 RI 1.12 b33 1/1 3 1 0.21 CR 0.025<0.10 表表 7 判断矩阵判断矩阵 B3-b 数值表数值表 Table 7 Judgment matrix for B3-b B3 b21 b22 b33 W 指标 b31 1 2 1/2 0.31 max 3.06λ B32 1/2 1 1/2 0.20 CI 0.030 RI 0.58 b33 2 2 1 0.49 CR 0.025<0.10 表表 8 判断矩阵判断矩阵 B4-b 数值表数值表 Table 8 Judgment matrix for B4-b B4 b41 b42 W 指标 b41 1 2 0.50 b42 1 1 0.50 max 2.00λ CI 0.030 RI 0.00 CR 0<0.10 表表 9 判断矩阵判断矩阵 B5-b 数值表数值表 Table 9 Judgment matrix for B5-b B5 b21 b22 b33 W 指标 b51 1 2 1/2 0.30 max 3.02λ b52 1/2 1 1/3 0.16 CI 0.01 RI 0.58 b53 2 3 1 0.54 CR 0.017 < 0.10 4.4 层次总排序及一致性检验层次总排序及一致性检验 按照公式(9) 、 (10)分别计算CI和RI,并按 照公式(11)进行总排序的一致性检验 5 1 0.260.0280.400.025 0.100.0300.11 00.13 0.010 0.0073 0.0100.003300.0013 0.0219 ii i CIBCI ∑ (22) 5 0.26 1.120.400.58 0.100.580.11 00.13 0.580.2912 ii i RIB RI ∑ 0.23200.058000.0754 0.6566 (23) 0.0219 0.033 0.6566 CI CR RI <0.10 (24) 可见层次总排序结果同样具有满意的一致性。 表表 10 各层元素对目标层的权重各层元素对目标层的权重 Table 10 The weight of each element for target stratum B1 B2 B3 B4 B5 层次 b 0.260.400.10 0.11 0.13 权值排序 b11 0.20/ / / / 0.0527 b12 0.10/ / / / 0.02613 b13 0.31/ / / / 0.0813 b14 0.11 / / / / 0.02912 b15 0.29/ / / / 0.0754 b21 / 0.70/ / / 0.2801 b22 / 0.10/ / / 0.0409 b23 / 0.21/ / / 0.0842 b31 / / 0.31 / / 0.03111 b32 / / 0.20 / / 0.02015 b33 / / 0.49 / / 0.0498 b41 / / / 0.50 / 0.0556 b42 / / / 0.50 / 0.0556 b51 / / / / 0.30 0.03910 b52 / / / / 0.16 0.02114 b53 / / / / 0.54 0.0705 4.5 结果分析与讨论结果分析与讨论 由于在影响因素两两比较时,存在一定的主观 性,权重系数可能随着专家的侧重点不同而发生变 化,根据上述计算结果可知,按照尾矿库设施构成 将影响尾矿库安全运行的因素排序为 排洪系统 (权 重为0.40)>尾矿堆存系统(权重为0.26)>尾矿 管理与维护 (权重为0.13)>回水系统 (权重为0.11) >尾矿输送系统(权重为0.10) ;按照各构成设施系 统包括的子系统影响作用排序为尾矿库排洪和调 洪能力(权重为0.028)>排洪设施的安全性和可靠 性(权重为0.084)>尾矿坝内坡排渗及浸润线高度 (权重为0.081)>坝体抗震能力(权重为0.075)> 尾矿库设施管理与维护(权重为0.070)>回水设施 完好情况及回水泵站运行情况(权重同为0.55)> 筑坝方式及工艺(权重为0.52)>砂泵站运行情况 (权重为0.049)>排水隧洞的坚固性(权重为0.40) 685 岩 土 力 学 2008 年 >尾矿库管理机构设置(权重为0.039)>输送管道 能力及完好情况(权重为0.031)>堆积坝坡比及防 冲刷 (权重为0.029)>尾矿颗粒组成及物理特性 (权 重为0.026)>尾矿操作工的技术培训(权重为 0.021)>浓缩池及浓缩机的状况(权重为0.020) 。 综上所述,在尾矿库的生产运行期间,首先要 保障的是排洪系统的安全运行,其中排洪和调洪能 力是主要因素,因为尾矿库调洪能力或排洪能力不 足时,安全超高和库内干滩长度将不能满足要求, 容易造成溃坝事故;其次是尾矿筑坝堆存系统,其 中坝体内坡排渗及浸润线高度是主要因素,因为尾 矿库坝体浸润线过高,坝体抗滑稳定性不足,导致 坝体失稳;然后是尾矿库的管理与维护,在我国一 些小型尾矿库存在不根据正规设计施工,不按照规 程操作管理,基本没有防洪、排渗等安全设施等违 规现象,为杜绝尾矿库安全事故的发生,尾矿设施 的管理与维护显得极其重要,“三分设计, 七分管理” 充分说明这一点。 5 结 论 尾矿库是一个多系统的复杂工程,而系统中又 包含众多子系统,且各子系统之间相互作用、相互 影响,使得各变量之间的关系相对复杂。因此,根 据科学性、系统性和实用性的原则,结合尾矿库安 全运行受多因素影响的实际特点,运用不确定层次 分析法对尾矿库的安全运行进行分析研究,确定各 影响因素的权重,其分析结果与尾矿库的实际运行 情况基本相符。因此,在应用层次分析法对尾矿库 的安全运行状况分析时,将各影响因素进行排序并 确定各因素所起作用的比重,可以提高尾矿库安全 管理的针对性和有效性,同时也促使企业相关部门 在尾矿库日常管理和定期安全巡视、检查工作中有 侧重点地进行,防患于未然,真正做到“安全第一, 预防为主,综合治理” ,确保尾矿库的安全运行。 参参 考考 文文 献献 [1] 王凤江. 国外尾矿坝事故调查分析[J]. 金属矿山, 2004, 8增刊 49-52. [2] 国家安全生产监督管理总局. 2007 年年全国非煤矿山 事故分析[EB/OL]. http// www. Chinasafety.gov. cn [3] KEMPER T, SOMMER S. Estimate of heavy metal contamination in soils after a mining accident using reflectance spectroscopy[J]. Environmental Science and Technology, 2002, 3612 2 742-7 471. [4] 李明, 胡乃联, 于芳, 等. ANSYS 软件在尾矿坝稳定性 分析中的应用研究[J]. 金属矿山, 2005, 8 56-59. [5] 潘建平, 孔宪京, 邹德高. 尾矿坝地震液化稳定的简化 分析[J]. 水利学报, 2006, 3710 1 224-1 229. [6] 田文旗, 薛剑光. 尾矿库安全技术与管理[M]. 北京 煤炭工业出版社, 2006. [7] 王涛, 郭振世.尾矿库的危险有害因素辨识及安全管理 措施[J]. 黄金, 2008, 294 51-55. [8] 田军, 张朋柱, 王刊良, 等. 基于德尔菲法的专家意见 集成模型研究[J]. 系统工程理论与实践, 2004, 1 57 -62, 69. [9] 胡保生, 彭勤科. 系统工程原理与应用[M]. 北京 化 学工业出版社, 2007. 686