第36卷第8期 北京科技大学学报 Vol.36 No.8 2014年8月 Journal of University of Science and Technology Beijing Aug.2014 采矿方法优选多目标决策一致性组合权重研究及应用 谭玉叶12,宋卫东123四,李铁一》,夏鸿》 1)北京科技大学土木与环境工程学院,北京1000832)北京科技大学金属矿山高效开采与安全教育部重点实验室,北京100083 3)武钢矿业公司金山店铁和矿,武汉430080 ☒通信作者,E-mail:luobiha435@163.com 摘要以金山店铁矿充填采矿法优选为工程背景,将采矿方法优选视为多目标决策的系统工程问题来处理.综合考虑技 术、经济、资源利用、劳动效率、安全、环境等多方面因素,构建了较为全面的采矿方法优选多目标决策指标体系。分别采用层 次分析法、模糊判断矩阵及Delphi法对决策指标进行赋权,得到了三组单一权重.在此基础上,通过Kendall非参数检验及 Spearman等级相关性分析,确定了各赋权方法的相关程度,得到了较为合理的采矿方法优选多目标决策指标的组合权重,并 基于组合权重,采用模糊多目标决策确定了最优采矿方法为分段空场嗣后充填法.本文采用的组合权重赋权方法较好地解决 了多目标决策时指标权重的分配问题,提高了采矿方法优选指标权重确定的合理性及科学性. 关键词采矿方法:优选:多目标决策分析:权重 分类号TD853 Application and research on consistency combination weights for mining method optimization by multi-objective decision TAN Yu-ye2,S0 NG Wei-dong'2g,LI Tie-yi》,XIA Hong》 1)School of Civil and Environmental Engineering,University of Science and Technology Beijing,Beijing 100083,China 2)Key Laboratory of High-Efficient Mining and Safety of Metal Mines (Ministry of Education),University of Science and Technology Beijing,Beijing 100083,China 3)Jinshandian Iron Mine,Wuhan Iron Steel Group Minerals Co.Ltd.,Wuhan 430080,China Corresponding author,E-mail:luobiha435@163.com ABSTRACT Mining method selection was handled as a system engineering problem with multi-objective decision-making.First,a comprehensive multi-objective decision-making index system for mining method optimization was built in consideration of technical fea- sibility,economic benefit,resource utilization,labor productivity,safety and environment factors.Then,the analytic hierarchy process,fuzzy judgment matrix and Delphi method were used to get three groups of single weights for the decision-making index system. On this basis,a group of reasonable combination weights were determined through Kendall nonparametric test and Spearman rank corre- lation analysis.This provides a better solution for weight distribution in multi-objective decision-making,and can improve the rationali- ty and scientificity of weight-determining methods for the decision-making index system.Finally,Jinshandian iron mine was taken as a case study on mining method selection,and the optimal mining method,sublevel open stope with backfill,was determined by this com- bination weight and fuzzy multi-objective decision-making method. KEY WORDS mining methods:optimization:multi-objective decision analysis:weights 采矿方法的选择是地下金属矿床开采设计中最 益和矿产资源的回收利用程度0:而采矿方法的选 核心的内容,其选择正确与否直接影响矿山经济效 择受到多种因素的影响制约,不仅涉及技术经济、资 收稿日期:2013-08-12 基金项目:国家自然科学基金资助项目(51274023):高等学校博士学科点专项科研基金资助项目(20120006110022) DOI:10.13374/j.issn1001-053x.2014.08.019:http://journals.ustb.edu.cn
第 36 卷 第 8 期 2014 年 8 月 北京科技大学学报 Journal of University of Science and Technology Beijing Vol. 36 No. 8 Aug. 2014 采矿方法优选多目标决策一致性组合权重研究及应用 谭玉叶1,2) ,宋卫东1,2) ,李铁一3) ,夏 鸿3) 1) 北京科技大学土木与环境工程学院,北京 100083 2) 北京科技大学金属矿山高效开采与安全教育部重点实验室,北京 100083 3) 武钢矿业公司金山店铁矿,武汉 430080 通信作者,E-mail: luobiha435@ 163. com 摘 要 以金山店铁矿充填采矿法优选为工程背景,将采矿方法优选视为多目标决策的系统工程问题来处理. 综合考虑技 术、经济、资源利用、劳动效率、安全、环境等多方面因素,构建了较为全面的采矿方法优选多目标决策指标体系. 分别采用层 次分析法、模糊判断矩阵及 Delphi 法对决策指标进行赋权,得到了三组单一权重. 在此基础上,通过 Kendall 非参数检验及 Spearman 等级相关性分析,确定了各赋权方法的相关程度,得到了较为合理的采矿方法优选多目标决策指标的组合权重,并 基于组合权重,采用模糊多目标决策确定了最优采矿方法为分段空场嗣后充填法. 本文采用的组合权重赋权方法较好地解决 了多目标决策时指标权重的分配问题,提高了采矿方法优选指标权重确定的合理性及科学性. 关键词 采矿方法; 优选; 多目标决策分析; 权重 分类号 TD 853 Application and research on consistency combination weights for mining method optimization by multi-objective decision TAN Yu-ye1,2) ,SONG Wei-dong1,2) ,LI Tie-yi3) ,XIA Hong3) 1) School of Civil and Environmental Engineering,University of Science and Technology Beijing,Beijing 100083,China 2) Key Laboratory of High-Efficient Mining and Safety of Metal Mines ( Ministry of Education) ,University of Science and Technology Beijing,Beijing 100083,China 3) Jinshandian Iron Mine,Wuhan Iron & Steel Group Minerals Co. Ltd. ,Wuhan 430080,China Corresponding author,E-mail: luobiha435@ 163. com ABSTRACT Mining method selection was handled as a system engineering problem with multi-objective decision-making. First,a comprehensive multi-objective decision-making index system for mining method optimization was built in consideration of technical feasibility,economic benefit,resource utilization,labor productivity,safety and environment factors. Then,the analytic hierarchy process,fuzzy judgment matrix and Delphi method were used to get three groups of single weights for the decision-making index system. On this basis,a group of reasonable combination weights were determined through Kendall nonparametric test and Spearman rank correlation analysis. This provides a better solution for weight distribution in multi-objective decision-making,and can improve the rationality and scientificity of weight-determining methods for the decision-making index system. Finally,Jinshandian iron mine was taken as a case study on mining method selection,and the optimal mining method,sublevel open stope with backfill,was determined by this combination weight and fuzzy multi-objective decision-making method. KEY WORDS mining methods; optimization; multi-objective decision analysis; weights 收稿日期: 2013--08--12 基金项目: 国家自然科学基金资助项目( 51274023) ; 高等学校博士学科点专项科研基金资助项目( 20120006110022) DOI: 10. 13374 /j. issn1001--053x. 2014. 08. 019; http: / /journals. ustb. edu. cn 采矿方法的选择是地下金属矿床开采设计中最 核心的内容,其选择正确与否直接影响矿山经济效 益和矿产资源的回收利用程度[1]; 而采矿方法的选 择受到多种因素的影响制约,不仅涉及技术经济、资
·1116 北京科技大学学报 第36卷 源等方面,同时也要考虑采矿方法本身对环境带来 析法及熵值法属于客观赋权法,虽然避免了主观因 的污染和破坏,是一个多因素、多目标决策的复杂系 素的干扰,但应用时需要大量的样本数据,不适应于 统工程问题 采矿方法优选这种预测性决策中因.本文以金山店 对于多目标决策问题而言,决策指标的选取及 铁矿采矿方法优选为工程背景,在构建多目标决策 权重的确定具有举足轻重的地位,决策结果是否科 指标体系的基础上,采用组合赋权的思想集结层次 学合理,很大程度上取决于这两项工作的客观合理 分析法、模糊判断矩阵及德尔菲法法三种赋权方法 性.近20年来国内外采矿工作者积极寻求新的途 所得到的单一权重,以提高采矿方法优选的科学合 径,采用不同的现代数学方法对采矿方法优选的指 理性 标选取及权重确定问题进行了有益的探讨P-习.然 1工程地质概况 而,采矿方法的优选目前还没有形成较为完整的理 论体系,尚有许多不足之处,从优选指标的选取来 金山店张福山铁矿床西起罗咸益第3勘探线, 看,大多数采矿工作者都把注意力集中在技术经济 东经张福山至李道士52勘探线,东西长3500m,南 指标上,只有极少数的学者考虑到了安全、环境、资 北宽1000m,面积约3.5km2,属接触交代矽卡岩型 源等因素对采矿方法选择的影响:从指标权重的确 矿床.I号和Ⅱ号矿体是张福山矿床的两个主要矿 定看来,都是采用单一的主观权重确定方法,确定结 体,均属于急倾斜、多层、中厚矿体,具有岩石夹层 果随机性较大,缺乏科学性 多,矿体和夹层厚度变化大的特点;块状和浸染状磁 目前指标权重确定方法主要有层次分析法(an- 铁石属于中等稳固,粉状磁铁石松软破碎极不稳固 alytic hierarchy process,AHP)、模糊判断矩阵法、德 (I号矿体粉矿占27.15%,Ⅱ号矿体粉矿占 尔菲法(Delphi)、主成分分析法、熵值法等.其中, 37.94%),如表1所示.考虑到块状和浸染状磁铁 层次分析法、模糊判断矩阵及德尔菲法法属于主观 石所占比例较大,约占总开采矿量的65%以上,本 赋权法,易于应用,但存在一定的主观性,处理不妥 次研究主要针对稳定性中等稳固的块状与浸染状矿 当时往往会夸大或降低某些指标的作用:主成分分 体采用组合权重进行采矿方法选择. 表1张福山和矿旷床开采技术条件 Table 1 Technical mining condition of Zhangfushan iron ore deposit 矿体赋存条件 围岩条件 矿体 矿体厚 倾角/ 走向矿石品 上盘 下盘 矿石稳固性 度/m () 长1m位/% 名称 稳固性 名称 稳固性 0.18-83.96 50-85 块状、浸染状、粉状磁铁石: 矽卡岩、泥质黑云 石英闪长岩、 I号 2690 43.11 中等稳固 中等稳固 平均24.50平均69 中等稳固~极不稳固 母角岩、粉砂岩 大理岩 0.70-131.3150~81 块状、浸染状、粉状磁铁石:矽卡岩、泥质黑云 石英闪长岩、 Ⅱ号 1020 39.02 中等稳固 中等稳固 平均26.00 平均67 中等稳固~极不稳固 母角岩、粉砂岩 大理岩 张福山矿床以25勘探线的中央主井为界,划分 联分析,综合考虑采矿方法技术可行性、矿山经济效 为东、西两个矿区,均采用无底柱分段崩落法开采 益、矿产资源的回收利用程度、劳动生产效率、回采 西区地表早己形成塌陷坑,并采用废石进行回填处 工作面安全程度、环境保护程度等多个方面因素,通 理:东区从2007年开始出现地表塌陷,塌陷的速度 过Delphi咨询统计及反复修改,将定性与定量、成 和范围远超出了预测,严重影响到地表建筑及村庄 本与收益、效率与安全等各个指标因素集成起来,以 的安全,一度造成东区开采的完全停止.矿山从生 整体最优为目标,构建了一个多层次、多耦合的,由 产安全、保护生态环境等多方面出发,考虑改用充填 低层详细表述到高层决策、特征提取和评估的采矿 法进行后续矿石的回采 方法优选多目标决策指标体系,如图1所示 该指标体系共包括三层:目标层(A)为最优采 2 构建多目标决策指标体系 矿方法;准则层(P)为技术因素P,、经济因素P2、资 为进行全面科学的决策,首先需要建立合理的 源利用率P,、劳动生产率P,、安全因素P,和环境效 评价指标体系,通过对金山店铁矿及国内典型地下 应P。:指标层包括矿体稳固性x1、矿体倾角x,等18 金属矿山的生产流程调查及流程中各环节指标的关 个指标
北 京 科 技 大 学 学 报 第 36 卷 源等方面,同时也要考虑采矿方法本身对环境带来 的污染和破坏,是一个多因素、多目标决策的复杂系 统工程问题. 对于多目标决策问题而言,决策指标的选取及 权重的确定具有举足轻重的地位,决策结果是否科 学合理,很大程度上取决于这两项工作的客观合理 性. 近 20 年来国内外采矿工作者积极寻求新的途 径,采用不同的现代数学方法对采矿方法优选的指 标选取及权重确定问题进行了有益的探讨[2 - 5]. 然 而,采矿方法的优选目前还没有形成较为完整的理 论体系,尚有许多不足之处. 从优选指标的选取来 看,大多数采矿工作者都把注意力集中在技术经济 指标上,只有极少数的学者考虑到了安全、环境、资 源等因素对采矿方法选择的影响; 从指标权重的确 定看来,都是采用单一的主观权重确定方法,确定结 果随机性较大,缺乏科学性. 目前指标权重确定方法主要有层次分析法( analytic hierarchy process,AHP) 、模糊判断矩阵法、德 尔菲法( Delphi) 、主成分分析法、熵值法等. 其中, 层次分析法、模糊判断矩阵及德尔菲法法属于主观 赋权法,易于应用,但存在一定的主观性,处理不妥 当时往往会夸大或降低某些指标的作用; 主成分分 析法及熵值法属于客观赋权法,虽然避免了主观因 素的干扰,但应用时需要大量的样本数据,不适应于 采矿方法优选这种预测性决策中[6]. 本文以金山店 铁矿采矿方法优选为工程背景,在构建多目标决策 指标体系的基础上,采用组合赋权的思想集结层次 分析法、模糊判断矩阵及德尔菲法法三种赋权方法 所得到的单一权重,以提高采矿方法优选的科学合 理性. 1 工程地质概况 金山店张福山铁矿床西起罗咸益第 3 勘探线, 东经张福山至李道士 52 勘探线,东西长 3500 m,南 北宽 1000 m,面积约 3. 5 km2 ,属接触交代矽卡岩型 矿床. Ⅰ号和Ⅱ号矿体是张福山矿床的两个主要矿 体,均属于急倾斜、多层、中厚矿体,具有岩石夹层 多,矿体和夹层厚度变化大的特点; 块状和浸染状磁 铁石属于中等稳固,粉状磁铁石松软破碎极不稳固 ( Ⅰ 号 矿 体 粉 矿 占 27. 15% ,Ⅱ 号 矿 体 粉 矿 占 37. 94% ) ,如表 1 所示. 考虑到块状和浸染状磁铁 石所占比例较大,约占总开采矿量的 65% 以上,本 次研究主要针对稳定性中等稳固的块状与浸染状矿 体采用组合权重进行采矿方法选择. 表 1 张福山矿床开采技术条件 Table 1 Technical mining condition of Zhangfushan iron ore deposit 矿体 矿体赋存条件 围岩条件 矿体厚 度/m 倾角/ ( °) 走向 长/m 矿石品 位/% 矿石稳固性 上盘 下盘 名称 稳固性 名称 稳固性 Ⅰ号 0. 18 ~ 83. 96 平均 24. 50 50 ~ 85 平均 69 2690 43. 11 块状、浸染状、粉状磁铁石: 中等稳固 ~ 极不稳固 矽卡岩、泥质黑云 母角岩、粉砂岩 中等稳固 石英闪长岩、 大理岩 中等稳固 Ⅱ号 0. 70 ~ 131. 31 平均 26. 00 50 ~ 81 平均 67 1020 39. 02 块状、浸染状、粉状磁铁石: 中等稳固 ~ 极不稳固 矽卡岩、泥质黑云 母角岩、粉砂岩 中等稳固 石英闪长岩、 大理岩 中等稳固 张福山矿床以 25 勘探线的中央主井为界,划分 为东、西两个矿区,均采用无底柱分段崩落法开采. 西区地表早已形成塌陷坑,并采用废石进行回填处 理; 东区从 2007 年开始出现地表塌陷,塌陷的速度 和范围远超出了预测,严重影响到地表建筑及村庄 的安全,一度造成东区开采的完全停止. 矿山从生 产安全、保护生态环境等多方面出发,考虑改用充填 法进行后续矿石的回采. 2 构建多目标决策指标体系 为进行全面科学的决策,首先需要建立合理的 评价指标体系. 通过对金山店铁矿及国内典型地下 金属矿山的生产流程调查及流程中各环节指标的关 联分析,综合考虑采矿方法技术可行性、矿山经济效 益、矿产资源的回收利用程度、劳动生产效率、回采 工作面安全程度、环境保护程度等多个方面因素,通 过 Delphi 咨询统计及反复修改,将定性与定量、成 本与收益、效率与安全等各个指标因素集成起来,以 整体最优为目标,构建了一个多层次、多耦合的,由 低层详细表述到高层决策、特征提取和评估的采矿 方法优选多目标决策指标体系,如图 1 所示. 该指标体系共包括三层: 目标层( A) 为最优采 矿方法; 准则层( P) 为技术因素 P1、经济因素 P2、资 源利用率 P3、劳动生产率 P4、安全因素 P5和环境效 应 P6 ; 指标层包括矿体稳固性 x1、矿体倾角 x2等 18 个指标. · 6111 ·
第8期 谭玉叶等:采矿方法优选多目标决策一致性组合权重研究及应用 ·1117· 最优采矿方法 (A) 技术因素 经济因素 资源利用率 劳动生产率 安全因素 环境效应 (P) (P) (P) (P) (P) (P) 矿 组 采 通 作 体稳固性 体倾角 岩稳固性 采切 采矿回 石 度 成本 收率 贫化率 生产能力 采矿 活适应性 简程 织管理难度 动强度 条件 面安全性 灾害 环境污染 5 图1 采矿方法优选多目标决策指标体系 Fig.1 Multiobjective decision-making index system for mining method optimization 的内在组合关系,在矿业的许多方面得到广泛应用 3 采矿方法模糊聚类初选 并取得了良好的效果 采矿方法选择一般分两步进行,首先需要根据 参照本文作者文献采矿方法模糊聚类初选的步 矿体的开采技术条件(技术因素)进行初选,再采用 骤可,根据所研究矿体的开采技术条件(表1),初选 数值或者数学方法进行更为细化的优选.本节采用 出三种充填采矿法,分别为分层充填采矿法A,、分 模糊聚类法按照矿体稳固性x1、矿体倾角x2、矿体 段空场嗣后充填法A,和阶段空场嗣后充填法A,· 厚度x?、围岩稳固性x四个技术因素的相似程度进 结合金山店铁矿的开采技术条件,并参照国内外矿 行初选 山采用上述采矿方法开采所取得的效果,确定相关 模糊聚类分析作为一种对所研究的事物按相似 技术经济指标如表2所示. 程度或亲疏关系进行分类的方法,能够反应样本间 表2初选采矿方案技术经济指标 Table 2 Technical and economic indicators for the primary selection of mining methods 采矿法 采矿成本/(元t) 采切比/(m3…万t) 采矿回收率/% 矿石贫化率/% 矿块生产能力/(td1) A 110 620 93 7 91 A2 105 475 90 8 265 A3 8 185 88.5 9 750 4 较,从而构造层次判断矩阵。有些学者采取将所有 基于一致性程度组合赋权 元素都和某个元素比较,即n个元素只做n-1次比 前文已采用模糊聚类法按照决策指标x,、x2、x3 较.这种做法的缺点在于,任何一个判断的失误都 和x4四个技术因素进行初选,故在此只对决策指标 可导致不合理的排序的出现,而个别评价者某个判 x;~x8进行赋权.不同的采矿方法,各指标的权值 断的失误在如此繁复的比较中往往是难以避免的 会有所差别,本文是针对充填法进行采矿方法优选, 本文将所有的元素都与其他元素做比较判断,进行 因此在后文采用不同的方法进行赋权时,均基于充 2n(n-1)次不同角度的反复比较,可以提供更多更 填法进行. 全面的信息,从而导出合理的排序.根据1~9比较 4.1 决策指标层次分析赋权 标准度表回,对指标进行重要度评价,构层次判断 层次分析法能够把复杂问题中的各因素划分成 矩阵R,其中R→和R,-如表3所示 相关联的有序层次,并采用数学方法确定各个元素 4.1.2计算层次相对权重 的相对重要次序,是一种较好的权值确定方法网 根据所得到的层次判断矩阵,采用方根法求出 采用层次分析法确定指标权重,可按以下步骤进行: 其最大特征根入m所对应的特征向量W,并对W进 ①构造层次判断矩阵;②计算层次相对权重;③一致 行归一化处理,便可得到各个决策指标的相对权 性检验.详述如下 重0.所采用的公式如下: 4.1.1构造层次判断矩阵 、(RW): 采用二元对比法对同层次的相关因素进行比 入mx= (1) nW
第 8 期 谭玉叶等: 采矿方法优选多目标决策一致性组合权重研究及应用 图 1 采矿方法优选多目标决策指标体系 Fig. 1 Multi-objective decision-making index system for mining method optimization 3 采矿方法模糊聚类初选 采矿方法选择一般分两步进行,首先需要根据 矿体的开采技术条件( 技术因素) 进行初选,再采用 数值或者数学方法进行更为细化的优选. 本节采用 模糊聚类法按照矿体稳固性 x1、矿体倾角 x2、矿体 厚度 x3、围岩稳固性 x4四个技术因素的相似程度进 行初选. 模糊聚类分析作为一种对所研究的事物按相似 程度或亲疏关系进行分类的方法,能够反应样本间 的内在组合关系,在矿业的许多方面得到广泛应用 并取得了良好的效果. 参照本文作者文献采矿方法模糊聚类初选的步 骤[7],根据所研究矿体的开采技术条件( 表 1) ,初选 出三种充填采矿法,分别为分层充填采矿法 A1、分 段空场嗣后充填法 A2 和阶段空场嗣后充填法 A3 . 结合金山店铁矿的开采技术条件,并参照国内外矿 山采用上述采矿方法开采所取得的效果,确定相关 技术经济指标如表 2 所示. 表 2 初选采矿方案技术经济指标 Table 2 Technical and economic indicators for the primary selection of mining methods 采矿法 采矿成本/( 元·t - 1 ) 采切比/( m3 ·万 t - 1 ) 采矿回收率/% 矿石贫化率/% 矿块生产能力/( t·d - 1 ) A1 110 620 93 7 91 A2 105 475 90 8 265 A3 95 185 88. 5 9 750 4 基于一致性程度组合赋权 前文已采用模糊聚类法按照决策指标 x1、x2、x3 和 x4四个技术因素进行初选,故在此只对决策指标 x5 ~ x18进行赋权. 不同的采矿方法,各指标的权值 会有所差别,本文是针对充填法进行采矿方法优选, 因此在后文采用不同的方法进行赋权时,均基于充 填法进行. 4. 1 决策指标层次分析赋权 层次分析法能够把复杂问题中的各因素划分成 相关联的有序层次,并采用数学方法确定各个元素 的相对重要次序,是一种较好的权值确定方法[8]. 采用层次分析法确定指标权重,可按以下步骤进行: ①构造层次判断矩阵; ②计算层次相对权重; ③一致 性检验. 详述如下. 4. 1. 1 构造层次判断矩阵 采用二元对比法对同层次的相关因素进行比 较,从而构造层次判断矩阵. 有些学者采取将所有 元素都和某个元素比较,即 n 个元素只做 n - 1 次比 较. 这种做法的缺点在于,任何一个判断的失误都 可导致不合理的排序的出现,而个别评价者某个判 断的失误在如此繁复的比较中往往是难以避免的. 本文将所有的元素都与其他元素做比较判断,进行 2n( n - 1) 次不同角度的反复比较,可以提供更多更 全面的信息,从而导出合理的排序. 根据 1 ~ 9 比较 标准度表[9],对指标进行重要度评价,构层次判断 矩阵 R,其中 RA--P和 RP4--X如表 3 所示. 4. 1. 2 计算层次相对权重 根据所得到的层次判断矩阵,采用方根法求出 其最大特征根 λmax所对应的特征向量 W,并对 W 进 行归一化处理,便可得到各个决策指标的相对权 重[10]. 所采用的公式如下: λmax = ∑ n i = 1 ( RW) i nWi . ( 1) · 7111 ·
·1118 北京科技大学学报 第36卷 表3指标层次判断矩阵 Table 3 AHP judgment matrix A-P P2 P P Ps P-X g x10 x11 x12 13 x14 P2 1 2 3 1 g 3 4 5 3 P 12 1/2 3 x10 13 1 2 12 3 1 1/3 2 1/3 2 11 1/4 12 1 13 12 Ps 1 2 3 x12 12 2 3 4 1/4 1/3 12 1/4 13 15 13 12 1/4 13 x14 13 1 2 1/2 1 1,2,…,mj=1,2,…,n. 所得到的层次权重计算结 (2) 果如表4. 式中:W为第j个特征向量所对应的第i个元素;i= 表4层次相对权重计算结果 Table 4 Relative weight calculation results by AHP 准则 A-P 指标 P2-X 指标 P3-X 指标 P-X 指标 P5-X 指标 P6-X 0.3089 5 0.4996 0.5695 0.3639 x14 0.0709 17 0.6667 Ps 0.1774 0.0790 0.3331 x10 0.1364 x15 0.1543 X18 0.3333 P 0.1391 0.2809 0.0974 0.0803 0.4959 Ps 0.3089 0.1405 12 0.2317 x17 0.2789 0.0657 13 0.0512 x14 0.1364 4.1.3判断矩阵一致性检验 w0=(0.1543,0.0244,0.1878,0.1025,0.0506, 需要对层次判断矩阵进行一致性检验,以判断 0.0190,0.0285,0.0322,0.0071,0.0409, 上述层次相对权重是否合理,检验公式为 0.0477,0.1532,0.1299,0.0219)T CR =C1/RI, (3) 4.2决策指标模糊判断矩阵赋权 Ci=(入mx-n)/(n-l) (4) 采用模糊判断矩阵确定指标权重,可按下面步 式中:C,为一致性检验指标;R为平均随机一致性指 骤进行:①构造Fuzzy判断矩阵:②对判断矩阵 标,对n=1,2,…,9,心理学家Saaty给出了R的 进行一致性检验;③对不一致性元素进行调整:④确 值.当Cs<0.1时认为判断举证具有满意的一致 定指标模糊权重.详述如下. 性,否则需对其进行调整以使其具有满意的一致 4.2.1构造模糊判断矩阵 性).检验结果(表5)表明各层次判断矩阵均通过 根据0.1-0.9标度表,通过专家对采矿优 致性检验 选指标进行重要度模糊评价,仍采用二元对比法进 表5层次判断矩阵一致性检验结果 行2n(n-1)次比较,构造准则层(P)及指标层(X) Table 5 Consistency test results of AHP judgment matrix 各因素的模糊判断矩阵R,其中RAp和R-x如表6 判断矩阵 入m C R CR 所示. A-P 5.4814 0.1103 1.12 0.0985 4.2.2判断矩阵一致性检验 P2-X 4.0277 0.0092 0.90 0.0103 受到指标的多样性、比较的繁复性及评价者思 P:-X 3.0246 0.0123 0.58 0.0212 维的不一致等多方面的主客观因素的影响,所得到 P-X 6.0730 0.0146 1.24 0.0118 的模糊判断矩阵往往不具有完全一致性,因而需要 Pj-X 4.0725 0.0242 0.90 0.0269 对其进行一致性检验及调整, P6-X 2.0000 0 0 0 采用指标C来表示模糊判断矩阵R的不一致 整理可得一组层次权重矩阵: 性程度:
北 京 科 技 大 学 学 报 第 36 卷 表 3 指标层次判断矩阵 Table 3 AHP judgment matrix A--P P2 P3 P4 P5 P6 P4 --X x9 x10 x11 x12 x13 x14 P2 1 2 3 1 4 x9 1 3 4 2 5 3 P3 1 /2 1 2 1 /2 3 x10 1 /3 1 2 1 /2 3 1 P4 1 /3 2 1 1 /3 2 x11 1 /4 1 /2 1 1 /3 2 1 /2 P5 1 2 3 1 4 x12 1 /2 2 3 1 4 2 P6 1 /4 1 /3 1 /2 1 /4 1 x13 1 /5 1 /3 1 /2 1 /4 1 1 /3 x14 1 /3 1 2 1 /2 3 1 Wi = n ∏ n j = 1 W 槡 ij ∑ n j = 1 Wj . ( 2) 式中: Wij为第 j 个特征向量所对应的第 i 个元素; i = 1,2,…,m; j = 1,2,…,n. 所得到的层次权重计算结 果如表 4. 表 4 层次相对权重计算结果 Table 4 Relative weight calculation results by AHP 准则 A--P 指标 P2--X 指标 P3--X 指标 P4--X 指标 P5--X 指标 P6--X P2 0. 3089 x5 0. 4996 x7 0. 5695 x9 0. 3639 x14 0. 0709 x17 0. 6667 P3 0. 1774 x6 0. 0790 x8 0. 3331 x10 0. 1364 x15 0. 1543 x18 0. 3333 P4 0. 1391 x7 0. 2809 x11 0. 0974 x11 0. 0803 x16 0. 4959 P5 0. 3089 x8 0. 1405 x12 0. 2317 x17 0. 2789 P6 0. 0657 x13 0. 0512 x14 0. 1364 4. 1. 3 判断矩阵一致性检验 需要对层次判断矩阵进行一致性检验,以判断 上述层次相对权重是否合理,检验公式为 CR = CI /RI, ( 3) CI = ( λmax - n) /( n - 1) . ( 4) 式中: CI为一致性检验指标; RI为平均随机一致性指 标,对 n = 1,2,…,9,心理学家 Saaty 给出了 RI 的 值. 当 CR < 0. 1 时认为判断举证具有满意的一致 性,否则需对其进行调整以使其具有满意的一致 性[11]. 检验结果( 表5) 表明各层次判断矩阵均通过 一致性检验. 表 5 层次判断矩阵一致性检验结果 Table 5 Consistency test results of AHP judgment matrix 判断矩阵 λmax CI RI CR A--P 5. 4814 0. 1103 1. 12 0. 0985 P2--X 4. 0277 0. 0092 0. 90 0. 0103 P3--X 3. 0246 0. 0123 0. 58 0. 0212 P4--X 6. 0730 0. 0146 1. 24 0. 0118 P5--X 4. 0725 0. 0242 0. 90 0. 0269 P6--X 2. 0000 0 0 0 整理可得一组层次权重矩阵: W( 1) = ( 0. 1543,0. 0244,0. 1878,0. 1025,0. 0506, 0. 0190,0. 0285,0. 0322,0. 0071,0. 0409, 0. 0477,0. 1532,0. 1299,0. 0219) T . 4. 2 决策指标模糊判断矩阵赋权 采用模糊判断矩阵确定指标权重,可按下面步 骤进行[12]: ①构造 Fuzzy 判断矩阵; ②对判断矩阵 进行一致性检验; ③对不一致性元素进行调整; ④确 定指标模糊权重. 详述如下. 4. 2. 1 构造模糊判断矩阵 根据 0. 1 ~ 0. 9 标度表[13],通过专家对采矿优 选指标进行重要度模糊评价,仍采用二元对比法进 行 2n( n - 1) 次比较,构造准则层( P) 及指标层( X) 各因素的模糊判断矩阵 R',其中 R'A--P和 R'P4--X如表 6 所示. 4. 2. 2 判断矩阵一致性检验 受到指标的多样性、比较的繁复性及评价者思 维的不一致等多方面的主客观因素的影响,所得到 的模糊判断矩阵往往不具有完全一致性,因而需要 对其进行一致性检验及调整[14]. 采用指标 IC 来表示模糊判断矩阵 R 的不一致 性程度: · 8111 ·
第8期 谭玉叶等:采矿方法优选多目标决策一致性组合权重研究及应用 ·1119· 表6指标模糊判断矩阵 Table 6 Fuzzy judgment matrix A-P P2 P P Ps Ps P-X g x10 x11 x12 13 x14 P2 0.5 0.7 0.8 0.5 0.9 g 05 0.7 0.8 0.6 0.9 0.8 P 0.3 0.5 0.4 0.3 0.8 10 0.3 0.5 0.7 0.4 0.8 0.5 P 0.2 0.6 0.5 0.4 0.8 11 0.2 0.3 0.5 0.3 0.7 0.3 Ps 0.5 0.7 0.6 0.5 0.9 x12 0.4 0.6 0.7 0.5 0.8 0.9 0.1 0.2 0.2 0.1 0.5 x13 0.1 0.2 0.2 0.2 0.5 0.2 14 0.2 0.5 0.1 0.1 0.8 0.5 阵R的不一致性的影响程度: 表7模糊判断矩阵一致性检验结果 IC n(n-1)(n-2) 5 Table 7 Consistency check results of the fuzzy judgment matrix 判断矩阵AP P,-X P-X P-X P-X 式中,,为x,与之间的二元对比值,n为模糊判断 IC 0.09000.05000.10000.15670.0500 矩阵的阶数.当R为一致性判断矩阵时,IC=0;否 则IC>0,且IC越大,表明R的一致性程度越差. 各判断矩阵的IC值如表7所示.从检验结果 N(r)= (6) 可知,各模糊判断矩阵均存在不一致性元素,均未能 n-2 通过一致性检验.其中Rx阶数为2,无需进行一 所得判断矩阵AP和P4-X的N(r)矩阵如表 致性检验 8所示.可知,模糊判断矩阵RP4-x的不一致元素较 采用N(r)来表示二元对比值r,对模糊判断矩 多,相应的其IC值也较大 表8模糊判断矩阵N(r)计算结果 Table 8 N(r)of the fuzzy judgment matrix A-P P2 P P. P Ps P.-X x10 x11 x12 13 14 P2 0 0.10000.2000 0.06670.1000 g 0 0.07500.17500.0500 0.10000.1750 P 0.1000 0 0.1000 0.0333 0.1000 x10 0.0750 0 0.2500 0.1250 0.07500.1000 P 0.20000.1000 0 0.06670.1000 x11 0.17500.2500 0 0.15000.10000.3000 Ps 0.0667 0.0333 0.0667 0 0.0333 x12 0.0500 0.1250 0.1500 0 0.15000.3250 Po 0.1000 0.1000 0.1000 0.0333 0 13 0.1000 0.0750 0.1000 0.1500 0 0.2000 14 0.1750 0.1000 0.3000 0.3250 0.2000 0 4.2.3判断矩阵一致性调整 一致性影响最大的二元对比值进行调整,直至得到 选取N(,)最大的二元对比值,即对模糊判断 具有满意的一致性为止. 调整并检验后的判断 矩阵R不一致性影响最大的值进行调整,得R:若 矩阵A-P和P,-X如表9所示 R仍不具有满意的一致性,则继续选取对R;的不 表9调整后的模糊判断矩阵 Table 9 Fuzzy judgment matrix adjusted A-P P2 P P5 Ps P-X g x10 12 x13 14 P2 0.5 0.7 0.6 0.5 0.9 g 0.5 0.7 0.6 0.6 0.9 0.7 P3 0.3 0.5 0.4 0.3 0.7 10 0.3 0.5 0.4 0.4 0.7 0.5 P 0.4 0.6 0.5 0.4 0.8 x11 0.4 0.6 0.5 0.5 0.8 0.6 0.5 0.7 0.6 0.5 0.9 x12 0.4 0.6 0.5 0.5 0.8 0.6 0.1 0.3 0.2 0.1 0.5 x13 0.1 0.3 0.2 0.2 0.5 0.3 14 0.3 0.5 0.4 0.4 0.7 0.5
第 8 期 谭玉叶等: 采矿方法优选多目标决策一致性组合权重研究及应用 表 6 指标模糊判断矩阵 Table 6 Fuzzy judgment matrix A--P P2 P3 P4 P5 P6 P4 --X x9 x10 x11 x12 x13 x14 P2 0. 5 0. 7 0. 8 0. 5 0. 9 x9 0. 5 0. 7 0. 8 0. 6 0. 9 0. 8 P3 0. 3 0. 5 0. 4 0. 3 0. 8 x10 0. 3 0. 5 0. 7 0. 4 0. 8 0. 5 P4 0. 2 0. 6 0. 5 0. 4 0. 8 x11 0. 2 0. 3 0. 5 0. 3 0. 7 0. 3 P5 0. 5 0. 7 0. 6 0. 5 0. 9 x12 0. 4 0. 6 0. 7 0. 5 0. 8 0. 9 P6 0. 1 0. 2 0. 2 0. 1 0. 5 x13 0. 1 0. 2 0. 2 0. 2 0. 5 0. 2 x14 0. 2 0. 5 0. 1 0. 1 0. 8 0. 5 IC = ∑ n i = 1 ∑ n j = 1 j≠i ∑ n k = 1 k≠i,j rij - r( ik + rkj - ) 1 2 n( n - 1) ( n - 2) . ( 5) 式中,rij为 xi与 xj之间的二元对比值,n 为模糊判断 矩阵的阶数. 当 R 为一致性判断矩阵时,IC = 0; 否 则 IC > 0,且 IC 越大,表明 R 的一致性程度越差. 各判断矩阵的 IC 值如表 7 所示. 从检验结果 可知,各模糊判断矩阵均存在不一致性元素,均未能 通过一致性检验. 其中 R' P6-X阶数为 2,无需进行一 致性检验. 采用 N( rij) 来表示二元对比值 rij对模糊判断矩 阵 R 的不一致性的影响程度: 表 7 模糊判断矩阵一致性检验结果 Table 7 Consistency check results of the fuzzy judgment matrix 判断矩阵 A--P P2--X P3--X P4--X P5--X IC 0. 0900 0. 0500 0. 1000 0. 1567 0. 0500 N( rij) = ∑ n k = 1 k≠i,j rij - r( ik + rkj - ) 1 2 n - 2 . ( 6) 所得判断矩阵 A--P 和 P4--X 的 N( rij) 矩阵如表 8 所示. 可知,模糊判断矩阵 R' P4--X的不一致元素较 多,相应的其 IC 值也较大. 表 8 模糊判断矩阵 N( rij) 计算结果 Table 8 N( rij) of the fuzzy judgment matrix A--P P2 P3 P4 P5 P6 P4 --X x9 x10 x11 x12 x13 x14 P2 0 0. 1000 0. 2000 0. 0667 0. 1000 x9 0 0. 0750 0. 1750 0. 0500 0. 1000 0. 1750 P3 0. 1000 0 0. 1000 0. 0333 0. 1000 x10 0. 0750 0 0. 2500 0. 1250 0. 0750 0. 1000 P4 0. 2000 0. 1000 0 0. 0667 0. 1000 x11 0. 1750 0. 2500 0 0. 1500 0. 1000 0. 3000 P5 0. 0667 0. 0333 0. 0667 0 0. 0333 x12 0. 0500 0. 1250 0. 1500 0 0. 1500 0. 3250 P6 0. 1000 0. 1000 0. 1000 0. 0333 0 x13 0. 1000 0. 0750 0. 1000 0. 1500 0 0. 2000 x14 0. 1750 0. 1000 0. 3000 0. 3250 0. 2000 0 4. 2. 3 判断矩阵一致性调整 选取 N( rij) 最大的二元对比值,即对模糊判断 矩阵 R'不一致性影响最大的值进行调整,得 R' 1 ; 若 R' 1 仍不具有满意的一致性,则继续选取对 R' 1 的不 一致性影响最大的二元对比值进行调整,直至得到 具有满意的一致性为止[15]. 调整并检验后的判断 矩阵 A--P 和 P4--X 如表 9 所示. 表 9 调整后的模糊判断矩阵 Table 9 Fuzzy judgment matrix adjusted A--P P2 P3 P4 P5 P6 P4 --X x9 x10 x11 x12 x13 x14 P2 0. 5 0. 7 0. 6 0. 5 0. 9 x9 0. 5 0. 7 0. 6 0. 6 0. 9 0. 7 P3 0. 3 0. 5 0. 4 0. 3 0. 7 x10 0. 3 0. 5 0. 4 0. 4 0. 7 0. 5 P4 0. 4 0. 6 0. 5 0. 4 0. 8 x11 0. 4 0. 6 0. 5 0. 5 0. 8 0. 6 P5 0. 5 0. 7 0. 6 0. 5 0. 9 x12 0. 4 0. 6 0. 5 0. 5 0. 8 0. 6 P6 0. 1 0. 3 0. 2 0. 1 0. 5 x13 0. 1 0. 3 0. 2 0. 2 0. 5 0. 3 x14 0. 3 0. 5 0. 4 0. 4 0. 7 0. 5 · 9111 ·
·1120* 北京科技大学学报 第36卷 4.2.4确定模糊权重 4.4.2赋权结果相关性计算 由于各判断矩阵的元素个数较少(n≤10),且 采用斯皮尔曼(Spearman)等级相关系数来检验 均通过了一致性检验,可由模糊一致判断矩阵直接 赋权结果的相关性.第i种赋权法与第k种赋权 计算排序权向量W=(w1,心2,…,心n)T 法的等级相关系数记为P,则 w,=(2立-1/m-D]. (7) P=1- 6 n(n2-1)启 p0-p)2, (9) 计算得到各判断矩阵的排序权向量,整理可得 得到斯皮尔曼(Spearman)等级相关系数矩阵为 模糊权重: 「1.00000.91650.8901 w②=(0.0864,0.0352,0.1518,0.1273,0.0480, 0.91651.00000.9648 0.0336,0.0721,0.0408,0.0192,0.0880, L0.89010.96481.0000J 0.0640,0.0896,0.1056,0.0384)T P越大,则两种赋权方法的权重排序结果就越一致, 4.3决策指标Delphi咨询赋权 相关性越大.由Spearman等级相关系数矩阵可知: 邀请20位长期从事采矿方法设计、采矿工艺及 层次分析法和模糊判断矩阵的赋权结果相关性较 理论、采矿经济及矿山安全等方面研究的学者专家, 好,与Delphi咨询赋权结果相关性次之;模糊判断 发出采矿方法优选指标权重Delphi咨询表.经过统 矩阵赋权结果与Delphi咨询赋权结果相关性非常 计、整理可得一组Delphi权重: 好,相关系数高达0.9648.总体来说,三种赋权结果 w)=(0.0957,0.0299,0.1731,0.1206,0.0476, 的相关性较好,均大于0.8 0.0378,0.0654,0.0438,0.0174,0.0640, 采用如下公式计算第辶种赋权法与其余赋权法 0.0593,0.0721,0.1294,0.0438)T 的平均相关程度: 4.4决策指标组合赋权 1 P:= P (10) 以上分别采用了层次分析法、模糊判断矩阵及 g-1t-1i Delphi咨询对决策指标进行了赋权,接下来需采对 可得p1=0.9033,P2=0.9407,P3=0.9275.归一化 这三组赋权结果进行一致性检验,以选取合理组合 得B1=0.3259,B2=0.3349,B3=0.3347. 权重集结方法6 4.4.3基于一致性程度组合赋权 4.4.1赋权结果一致性检验 采用加权处理方法进行组合权重的集合,公式 当赋权方法9≥3时,一般采用非参数统计方法 如下: 中的肯德尔(Kendall)协和系数检验法来检验赋权 结果的一致性).计算公式为 W=∑BW (11) 计算得基于一致性程度的组合权重向量为 (p") W=(0.1117,0.0299,0.1707,0.1170,0.0487, K=j=1 \=1 3(n+1) 9n(n-1) n-1 (8) 0.0302,0.0556,0.0390,0.0147,0.0646, 0.0571,0.1045,0.1215,0.0348)T 式中:,为第k赋权方法所确定的权重的排序号 向量,k=1,2,…,9j=1,2,…,np,p2,…, 5 组合权重模糊多目标决策 p为1~n的一个全排列.根据前文所得的w= (W,W2,W),可得 5.1模糊隶属度矩阵的确定 对于一个多目标决策问题,由于各决策指标的属 P=(p0,p②,p)= 性不同、单位量纲不同、数量级及含义不同,在进行方 「134141092561781211 371 案决策之前,必须对决策隶属度矩阵进行规范化处理 11 4141363751981012 2 对于决策指标体系中的定量指标(如表9),分 1121412639518710134」 别采用公式cg=xg/max(xg)及cg=1-xg/max(xg) 计算可得K=2.0519.给定显著性水平a= 对收益性指标采矿回收率x,、矿块生产能力x,及消 0.O5,查肯德尔(Kendall)协和系数的临界值表得临 耗性指标采矿成本x、采切比x。、矿石贫化率x进 界值K。=0.363,即K.≤K,不能通过一致性检验 行规范化化处理.对于定性指标,如采矿工效和劳 因此,不宜采用加权平均的方法进行组合赋权,必须 动强度,采用优选关系法确定其隶属函数,并进行规 考虑各赋权结果之间相性,再进行加权处理. 范化处理呵,可得决策指标隶属度矩阵如下:
北 京 科 技 大 学 学 报 第 36 卷 4. 2. 4 确定模糊权重 由于各判断矩阵的元素个数较少( n≤10) ,且 均通过了一致性检验,可由模糊一致判断矩阵直接 计算排序权向量 W = ( w1,w2,…,wn ) T . wi = 2( ∑ n j = 1 rij - 1 ) [n( n - 1) ]. ( 7) 计算得到各判断矩阵的排序权向量,整理可得 模糊权重: W( 2) = ( 0. 0864,0. 0352,0. 1518,0. 1273,0. 0480, 0. 0336,0. 0721,0. 0408,0. 0192,0. 0880, 0. 0640,0. 0896,0. 1056,0. 0384) T . 4. 3 决策指标 Delphi 咨询赋权 邀请 20 位长期从事采矿方法设计、采矿工艺及 理论、采矿经济及矿山安全等方面研究的学者专家, 发出采矿方法优选指标权重 Delphi 咨询表. 经过统 计、整理可得一组 Delphi 权重: W( 3) = ( 0. 0957,0. 0299,0. 1731,0. 1206,0. 0476, 0. 0378,0. 0654,0. 0438,0. 0174,0. 0640, 0. 0593,0. 0721,0. 1294,0. 0438) T . 4. 4 决策指标组合赋权 以上分别采用了层次分析法、模糊判断矩阵及 Delphi 咨询对决策指标进行了赋权,接下来需采对 这三组赋权结果进行一致性检验,以选取合理组合 权重集结方法[16]. 4. 4. 1 赋权结果一致性检验 当赋权方法 q≥3 时,一般采用非参数统计方法 中的肯德尔( Kendall) 协和系数检验法来检验赋权 结果的一致性[17]. 计算公式为 K = 12∑ n j = ( 1 ∑ q k = 1 p( k) j ) 2 q 2 n( n - 1) - 3( n + 1) n - 1 . ( 8) 式中: p( k) j 为第 k 赋权方法所确定的权重的排序号 向量,k = 1,2,…,q,j = 1,2,…,n; p( k) 1 ,p( k) 2 ,…, p( k) n 为 1 ~ n 的一个全排列. 根据前文所得的 W = ( W( 1) ,W( 2) ,W( 3) ) ,可得 P = ( p( 1) ,p( 2) ,p( 3) ) = 13 4 14 10 9 2 5 6 1 7 8 12 11 3 11 4 14 13 6 3 7 5 1 9 8 10 12 2 11 2 14 12 6 3 9 5 1 8 7 10 13 4 T . 计算可得 K = 2. 0519. 给定显著性水平 α = 0. 05,查肯德尔( Kendall) 协和系数的临界值表得临 界值 Kα = 0. 363,即 Kα≤K,不能通过一致性检验. 因此,不宜采用加权平均的方法进行组合赋权,必须 考虑各赋权结果之间相性,再进行加权处理. 4. 4. 2 赋权结果相关性计算 采用斯皮尔曼( Spearman) 等级相关系数来检验 赋权结果的相关性[18]. 第 i 种赋权法与第 k 种赋权 法的等级相关系数记为 ρik,则 ρik = 1 - 6 n( n2 - 1) ∑ n j = 1 ( p( i) j - p( k) j ) 2 , ( 9) 得到斯皮尔曼( Spearman) 等级相关系数矩阵为 ρ = 1. 0000 0. 9165 0. 8901 0. 9165 1. 0000 0. 9648 0. 8901 0. 9648 1. 0000 . ρik越大,则两种赋权方法的权重排序结果就越一致, 相关性越大. 由 Spearman 等级相关系数矩阵可知: 层次分析法和模糊判断矩阵的赋权结果相关性较 好,与 Delphi 咨询赋权结果相关性次之; 模糊判断 矩阵赋权结果与 Delphi 咨询赋权结果相关性非常 好,相关系数高达 0. 9648. 总体来说,三种赋权结果 的相关性较好,均大于 0. 8. 采用如下公式计算第 i 种赋权法与其余赋权法 的平均相关程度: ρi = 1 q - 1 ∑ q k = 1,k≠i ρik . ( 10) 可得 ρ1 = 0. 9033,ρ2 = 0. 9407,ρ3 = 0. 9275. 归一化 得 β1 = 0. 3259,β2 = 0. 3349,β3 = 0. 3347. 4. 4. 3 基于一致性程度组合赋权 采用加权处理方法进行组合权重的集合,公式 如下: W = ∑ q k = 1 βkWk . ( 11) 计算得基于一致性程度的组合权重向量为 W = ( 0. 1117,0. 0299,0. 1707,0. 1170,0. 0487, 0. 0302,0. 0556,0. 0390,0. 0147,0. 0646, 0. 0571,0. 1045,0. 1215,0. 0348) T . 5 组合权重模糊多目标决策 5. 1 模糊隶属度矩阵的确定 对于一个多目标决策问题,由于各决策指标的属 性不同、单位量纲不同、数量级及含义不同,在进行方 案决策之前,必须对决策隶属度矩阵进行规范化处理. 对于决策指标体系中的定量指标( 如表 9) ,分 别采用公式 cij = xij /max( xij) 及 cij = 1 - xij /max( xij) 对收益性指标采矿回收率 x7、矿块生产能力 x9及消 耗性指标采矿成本 x5、采切比 x6、矿石贫化率 x8进 行规范化化处理. 对于定性指标,如采矿工效和劳 动强度,采用优选关系法确定其隶属函数,并进行规 范化处理[19],可得决策指标隶属度矩阵如下: · 0211 ·
第8期 谭玉叶等:采矿方法优选多目标决策一致性组合权重研究及应用 ·1121· 「0.00000.00001.00001.00000.00000.00001.00000.00000.00000.00000.00000.00001.00001.0000 C=0.33330.33330.57960.60650.26400.22450.76150.31710.38100.28570.82611.00000.63640.9000 L1.00001.00000.36790.36791.00001.00000.00001.00001.00001.00001.00000.50000.00000.0000 5.2模糊多目标决策 探索 采用模糊综合评判法进行采矿方法优选多目标 决策四.根据前文所得到的决策指标的权重矩阵 参考文献 W及隶属度矩阵C,按矩阵乘法计算权向量与隶属 Wang YG.Important impact caused by mining method choice onto 度矩阵的乘积,计算公式如下: mining engineering design and ore production.Nonferrous Mines, 1995(3):4 B=∑WR=(b,b2,…,bn). (12) (王玉国.采矿方法选择对工程设计和矿山生产的重要影响. 有色矿山,1995(3):4) 式中,b表示方案A,综合隶属度.根据方案的综合隶 Huang DZ,Hu Y Q,Chen X H,et al.Intelligent decision sup- 属度对方案集A进行排序,以综合隶属度最大的采 port system of mining method based on fuzzy expert system and 矿方法为最优采矿方法.通过计算可得到B,= neural net.Nonferrous Met,2002,54(2):104 (0.4996,0.5072,0.4482).即多目标决策结果为: (黄德锘,胡运权,陈孝华,等.基于模糊推理和神经网络的采 矿方法智能决策系统.有色金属,2002,54(2):104) 方案A,优于方案A1,方案A优于方案A3·最终确定 B] Wang X M,Xie S Q,Zhao B,et al.Optimal selection of the min- 方案A,分段空场嗣后充填法为适合金山店铁矿的 ing method based on grey relational method.Met Mine,2009(7): 最优采矿方法 10) (王新民,谢盛青,赵彬,等.基于灰关联分析法的采矿方法 6结论 优选.金属矿山,2009(7):10) 4 (1)综合考虑采矿方法技术可行性、矿山经济 Wu A X,Zhang W F,Guo L.Application of the similarity ratio value engineering method in optimization of mining method.Cent 效益、矿产资源的回收利用程度、劳动生产效率、回 South Univ Technol,2000,31(4):294 采工作面安全程度、环境保护程度等多方面因素,将 (吴爱祥,张卫锋,郭立.相似率价值工程法在采矿方法优选 定性与定量、成本与收益、效率与安全等各个指标因 中的应用.中南工业大学学报,2000,31(4):294) 素集成起来,以整体最优为目标,构建了由低层详细 5] Zhang X Y,Zhang H J.Application of fuzzy mathematics in min- 表述到高层决策的采矿方法优选多目标决策指标 ing method optimization in Jinchuan No.2 mine area.Nonferrous M,2013,64(2):80 体系 (张小义,张红军.模糊数学在金川二矿区采矿方法优选中的 (2)分别采用层次分析法、模糊判断矩阵及 应用.有色金属,2013,64(2):80) Delphi法对决策指标进行赋权,并通过一致性检验 [6 He Y T.Zou H.A dynamic empowerment method by multiple 及调整确定了较为客观的三组单一权重;基于组合 strategies for evaluation indices.J Suzhou Univ,2010,25(2): 赋权思想,通过Kendall非参数检验及Spearman等 47 级相关性分析,确定了各赋权方法相关程度,得到了 (何永太,邹海.一种多策略的评价指标动态赋权方法宿州 学院学报,2010,25(2):47) 较为合理采矿方法优选指标组合权重,解决了多目 ] Tan Y Y,Song W D,Lei Y K,et al.Synthetic judgment for min- 标决策时指标权重的分配问题,提高了采矿方法优 ing method optimization based on fuzzy cluster analysis and analyt- 选指标权重确定的合理性及科学性 ic hierarchy process.JUnie Sci Technol Beijing,2012,34(5): (3)在模糊聚类初选的基础上,分别采用合理 489 的规范化及定量化方法,通过模糊多目标决策,确定 (谭玉叶,宋卫东,雷远坤,等。基于模糊聚类及层次分析法 的采矿方法综合评判优选.北京科技大学学报,2012,34 方案A,分段空场嗣后充填法为金山店铁矿的最优 (5):489) 采矿方法 [8]Ren H G,Tan Z Y,Cai X F,et al.AHP-Fuzzy optimization of 本文所采用的组合赋权方法具有较好的逻辑性 structural parameters in sublevel openstope succedent filling meth- 及科学性,然而考虑到对权重值的分配仍然容易受 od.J Unin Sci Technol Beijing,2010,32(11):1383 到一定的主观因素影响以及对定性指标采用专家评 (任红岗,谭卓英,蔡学蜂,等。分段空场嗣后充填法采场结 分制等,优选所得到的结果可以作为采矿方法选择 构参数AHP-Fuzy优化.北京科技大学学报,2010,32(11): 1383) 时的重要参考依据,但不能作为惟一的决定因 9]Ye H W,Chang J.Choosing mining method based on fuzzy deci- 素0.同时,如何认识各个决策指标之间的内在联 sion and analytic hierarchy process.I Wuhan Unir Technol,2009, 系,提高采矿方法优选的科学性,还有待进一步研究 31(8):145
第 8 期 谭玉叶等: 采矿方法优选多目标决策一致性组合权重研究及应用 C = 0. 0000 0. 0000 1. 0000 1. 0000 0. 0000 0. 0000 1. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 1. 0000 1. 0000 0. 3333 0. 3333 0. 5796 0. 6065 0. 2640 0. 2245 0. 7615 0. 3171 0. 3810 0. 2857 0. 8261 1. 0000 0. 6364 0. 9000 1. 0000 1. 0000 0. 3679 0. 3679 1. 0000 1. 0000 0. 0000 1. 0000 1. 0000 1. 0000 1. 0000 0. 5000 0. 0000 0. 0000 . 5. 2 模糊多目标决策 采用模糊综合评判法进行采矿方法优选多目标 决策[20]. 根据前文所得到的决策指标的权重矩阵 W 及隶属度矩阵 C,按矩阵乘法计算权向量与隶属 度矩阵的乘积,计算公式如下: Bj = ∑ n i = 1 WiRij = ( b1,b2,…,bn ) . ( 12) 式中,bj表示方案 Aj综合隶属度. 根据方案的综合隶 属度对方案集 A 进行排序,以综合隶属度最大的采 矿方法为最优采矿方法. 通过计算可得到 Bj = ( 0. 4996,0. 5072,0. 4482) . 即多目标决策结果为: 方案 A2优于方案 A1,方案 A1优于方案 A3 . 最终确定 方案 A2分段空场嗣后充填法为适合金山店铁矿的 最优采矿方法. 6 结论 ( 1) 综合考虑采矿方法技术可行性、矿山经济 效益、矿产资源的回收利用程度、劳动生产效率、回 采工作面安全程度、环境保护程度等多方面因素,将 定性与定量、成本与收益、效率与安全等各个指标因 素集成起来,以整体最优为目标,构建了由低层详细 表述到高层决策的采矿方法优选多目标决策指标 体系. ( 2) 分别采用层次分析法、模糊判断矩阵及 Delphi 法对决策指标进行赋权,并通过一致性检验 及调整确定了较为客观的三组单一权重; 基于组合 赋权思想,通过 Kendall 非参数检验及 Spearman 等 级相关性分析,确定了各赋权方法相关程度,得到了 较为合理采矿方法优选指标组合权重,解决了多目 标决策时指标权重的分配问题,提高了采矿方法优 选指标权重确定的合理性及科学性. ( 3) 在模糊聚类初选的基础上,分别采用合理 的规范化及定量化方法,通过模糊多目标决策,确定 方案 A2分段空场嗣后充填法为金山店铁矿的最优 采矿方法. 本文所采用的组合赋权方法具有较好的逻辑性 及科学性,然而考虑到对权重值的分配仍然容易受 到一定的主观因素影响以及对定性指标采用专家评 分制等,优选所得到的结果可以作为采矿方法选择 时的重要参考依据,但不能作为惟一的决定因 素[20]. 同时,如何认识各个决策指标之间的内在联 系,提高采矿方法优选的科学性,还有待进一步研究 探索. 参 考 文 献 [1] Wang Y G. Important impact caused by mining method choice onto mining engineering design and ore production. Nonferrous Mines, 1995( 3) : 4 ( 王玉国. 采矿方法选择对工程设计和矿山生产的重要影响. 有色矿山,1995( 3) : 4) [2] Huang D Z,Hu Y Q,Chen X H,et al. Intelligent decision support system of mining method based on fuzzy expert system and neural net. Nonferrous Met,2002,54( 2) : 104 ( 黄德锗,胡运权,陈孝华,等. 基于模糊推理和神经网络的采 矿方法智能决策系统. 有色金属,2002,54( 2) : 104) [3] Wang X M,Xie S Q,Zhao B,et al. Optimal selection of the mining method based on grey relational method. Met Mine,2009( 7) : 10) ( 王新民,谢盛青,赵彬,等. 基于灰关联分析法的采矿方法 优选. 金属矿山,2009( 7) : 10) [4] Wu A X,Zhang W F,Guo L. Application of the similarity ratio value engineering method in optimization of mining method. J Cent South Univ Technol,2000,31( 4) : 294 ( 吴爱祥,张卫锋,郭立. 相似率价值工程法在采矿方法优选 中的应用. 中南工业大学学报,2000,31( 4) : 294) [5] Zhang X Y,Zhang H J. Application of fuzzy mathematics in mining method optimization in Jinchuan No. 2 mine area. Nonferrous Met,2013,64( 2) : 80 ( 张小义,张红军. 模糊数学在金川二矿区采矿方法优选中的 应用. 有色金属,2013,64( 2) : 80) [6] He Y T,Zou H. A dynamic empowerment method by multiple strategies for evaluation indices. J Suzhou Univ,2010,25 ( 2) : 47 ( 何永太,邹海. 一种多策略的评价指标动态赋权方法. 宿州 学院学报,2010,25( 2) : 47) [7] Tan Y Y,Song W D,Lei Y K,et al. Synthetic judgment for mining method optimization based on fuzzy cluster analysis and analytic hierarchy process. J Univ Sci Technol Beijing,2012,34( 5) : 489 ( 谭玉叶,宋卫东,雷远坤,等. 基于模糊聚类及层次分析法 的采矿方法 综 合 评 判 优 选. 北京科技大学学报,2012,34 ( 5) : 489) [8] Ren H G,Tan Z Y,Cai X F,et al. AHP-Fuzzy optimization of structural parameters in sublevel openstope succedent filling method. J Univ Sci Technol Beijing,2010,32( 11) : 1383 ( 任红岗,谭卓英,蔡学峰,等. 分段空场嗣后充填法采场结 构参数 AHP--Fuzzy 优化. 北京科技大学学报,2010,32( 11) : 1383) [9] Ye H W,Chang J. Choosing mining method based on fuzzy decision and analytic hierarchy process. J Wuhan Univ Technol,2009, 31( 8) : 145 · 1211 ·
·1122 北京科技大学学报 第36卷 (叶海旺,常剑.基于模糊决策和层次分析法的采矿方法选 Eng,2007,22(3):257 择.武汉理工大学学报,2009,31(8)145) (黎涛,朱建军,王梦光,等.模糊判断矩阵的一致性判定及 [10]Chang J E,Jiang T L.Research on the weight of coefficient 改进方法研究.系统工程学报,2007,22(3):257) through analytic hierarchy process.Wuhan Unie Technol, [16]Deng H W,Hu P L,Yang NG.Mining method choice based on 2007,29(1):153 empowerment combination TOPSIS.J Guangxi Unir Nat Sci Ed, (常建娥,蒋太立.层次分析法确定权重的研究.武汉理工 2012,37(5):990 大学学报,2007,29(1):153) (邓红卫,胡普仑,杨念哥.基于组合赋权TOPSIS的采矿方 [11]Wang X M,Zhao B,Zhang Q L.Mining method choice based on 法优选.广西大学学报:自然科学版,2012,37(5):990) AHP and fuzzy mathematics.I Cent South Unir Sci Technol, [17]Huang HT,Yang J Y,Fan Y T,et al.Deciding the combina- 2008,39(5):875 tion weights of urban land grading factors based on the correlated (王新民,赵彬,张软礼.基于层次分析和模糊数学的采矿 consistency.J Anhui Agric Sci,2009,37 (2):743 方法选择.中南大学学报:自然科学版,2008,39(5):875) (黄海涛,杨建宇,樊雅婷,等.基于相关一致性的城镇土地 [12]Liang B S.Cao DL.Fuzzy Mathematics and Its Applications. 定级因素综合权重的确定.安徽农业科学,2009,37(2): Beijing:Science Press,2007 743) (梁保松,曹殿立.模糊数学及其应用.北京:科学出版社, [18]Yang Y,Zhang X D.Research on the technology progresses to 2007) evaluate the index sign classification based on meta-evaluation.I [13]Yuan H J,Chen Y M,Zhao J F.The application of FAHP in Harbin Unir Commerce Soc Sci Ed,2007(3):77 priorities of the stope security's risk analysis.Nonferrous Met Sci (杨毅,张晓东.基于元评价的科技进步综合评价指标分类 Eng,2010,1(1):55 研究.哈尔滨商业大学学报:社会科学版,2007(3):77) (原虎军,陈玉明,赵继锋.FAPH在采场安全风险分析中权 09] Guo W B,Liu Y X,Li X S.Fuzzy clustering analysis of mining 重分配的应用.有色金属科学与工程,2010,1(1):55) induced damages of buildings./Min Saf Eng,2007,24(3): 14]Song G X,Yang D L.Methods for identifying and improving the 288 consistency of fuzzy judgment matrix.Syst Eng,2003,21 (1): (郭文兵,刘义新,李小双.采动影响下建筑物损害程度的 110 模糊聚类分析.采矿与安全工程学报,2007,24(3):288) (宋光兴,杨德礼。模糊判断矩阵的一致性检验及一致性改 [20] Gong S W,Cai M Y,Li X B.Application of AHP-FCE in opti- 进方法.系统工程,2003,21(1):110) mizing mining method.J Min Saf Eng,2009,26(1):41 [15]Li T,Zhu JJ,Wang M G,et al.Research on consistency test (龚声武,蔡明悦,李夕兵.AHP-FCE法在采矿方法优选中 and modification approach of fuzzy comparison matrix.Syst 的应用.采矿与安全工程学报,2009,26(1):41)
北 京 科 技 大 学 学 报 第 36 卷 ( 叶海旺,常剑. 基于模糊决策和层次分析法的采矿方法选 择. 武汉理工大学学报,2009,31( 8) : 145) [10] Chang J E,Jiang T L. Research on the weight of coefficient through analytic hierarchy process. J Wuhan Univ Technol, 2007,29( 1) : 153 ( 常建娥,蒋太立. 层次分析法确定权重的研究. 武汉理工 大学学报,2007,29( 1) : 153) [11] Wang X M,Zhao B,Zhang Q L. Mining method choice based on AHP and fuzzy mathematics. J Cent South Univ Sci Technol, 2008,39( 5) : 875 ( 王新民,赵彬,张钦礼. 基于层次分析和模糊数学的采矿 方法选择. 中南大学学报: 自然科学版,2008,39( 5) : 875) [12] Liang B S,Cao D L. Fuzzy Mathematics and Its Applications. Beijing: Science Press,2007 ( 梁保松,曹殿立. 模糊数学及其应用. 北京: 科学出版社, 2007) [13] Yuan H J,Chen Y M,Zhao J F. The application of FAHP in priorities of the stope security's risk analysis. Nonferrous Met Sci Eng,2010,1( 1) : 55 ( 原虎军,陈玉明,赵继锋. FAPH 在采场安全风险分析中权 重分配的应用. 有色金属科学与工程,2010,1( 1) : 55) [14] Song G X,Yang D L. Methods for identifying and improving the consistency of fuzzy judgment matrix. Syst Eng,2003,21 ( 1) : 110 ( 宋光兴,杨德礼. 模糊判断矩阵的一致性检验及一致性改 进方法. 系统工程,2003,21( 1) : 110) [15] Li T,Zhu J J,Wang M G,et al. Research on consistency test and modification approach of fuzzy comparison matrix. J Syst Eng,2007,22( 3) : 257 ( 黎涛,朱建军,王梦光,等. 模糊判断矩阵的一致性判定及 改进方法研究. 系统工程学报,2007,22( 3) : 257) [16] Deng H W,Hu P L,Yang N G. Mining method choice based on empowerment combination TOPSIS. J Guangxi Univ Nat Sci Ed, 2012,37( 5) : 990 ( 邓红卫,胡普仑,杨念哥. 基于组合赋权 TOPSIS 的采矿方 法优选. 广西大学学报: 自然科学版,2012,37( 5) : 990) [17] Huang H T,Yang J Y,Fan Y T,et al. Deciding the combination weights of urban land grading factors based on the correlated consistency. J Anhui Agric Sci,2009,37( 2) : 743 ( 黄海涛,杨建宇,樊雅婷,等. 基于相关一致性的城镇土地 定级因素综合权重的确定. 安徽农业科学,2009,37 ( 2) : 743) [18] Yang Y,Zhang X D. Research on the technology progresses to evaluate the index sign classification based on meta-evaluation. J Harbin Univ Commerce Soc Sci Ed,2007( 3) : 77 ( 杨毅,张晓东. 基于元评价的科技进步综合评价指标分类 研究. 哈尔滨商业大学学报: 社会科学版,2007( 3) : 77) [19] Guo W B,Liu Y X,Li X S. Fuzzy clustering analysis of mining induced damages of buildings. J Min Saf Eng,2007,24 ( 3) : 288 ( 郭文兵,刘义新,李小双. 采动影响下建筑物损害程度的 模糊聚类分析. 采矿与安全工程学报,2007,24( 3) : 288) [20] Gong S W,Cai M Y,Li X B. Application of AHP-FCE in optimizing mining method. J Min Saf Eng,2009,26( 1) : 41 ( 龚声武,蔡明悦,李夕兵. AHP--FCE 法在采矿方法优选中 的应用. 采矿与安全工程学报,2009,26( 1) : 41) · 2211 ·