CN108710984B - 一种矿山地质环境综合评价方法及系统 - Google Patents
一种矿山地质环境综合评价方法及系统 Download PDFInfo
- Publication number
- CN108710984B CN108710984B CN201810300100.4A CN201810300100A CN108710984B CN 108710984 B CN108710984 B CN 108710984B CN 201810300100 A CN201810300100 A CN 201810300100A CN 108710984 B CN108710984 B CN 108710984B
- Authority
- CN
- China
- Prior art keywords
- evaluation index
- weight
- evaluation
- effective
- geological environment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 481
- 239000011159 matrix material Substances 0.000 claims abstract description 92
- 238000000034 method Methods 0.000 claims abstract description 51
- 238000013139 quantization Methods 0.000 claims abstract description 38
- 238000012545 processing Methods 0.000 claims abstract description 22
- 238000005065 mining Methods 0.000 claims description 28
- 230000000295 complement effect Effects 0.000 claims description 18
- 239000013598 vector Substances 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000007613 environmental effect Effects 0.000 claims description 7
- 238000010606 normalization Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 4
- 238000004451 qualitative analysis Methods 0.000 claims description 3
- 238000004445 quantitative analysis Methods 0.000 claims description 3
- 238000011161 development Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 description 4
- 238000011835 investigation Methods 0.000 description 4
- 239000011707 mineral Substances 0.000 description 4
- 230000006378 damage Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 239000002910 solid waste Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Physics & Mathematics (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种矿山地质环境综合评价方法及系统,该方法包括建立矿山地质环境评价分层结构模型;该模型包括多个评价指标;构建评价指标判断矩阵,并根据此评价指标判断矩阵,确定有效评价指标;对所有有效评价指标进行无纲量化处理得到无纲量化权值;采用层次分析法计算每个有效评价指标的主观权重;利用信息熵计算每个有效评价指标的熵权值;采用层次分析法对熵权值的一致性关系进行检验确定每个有效评价指标的客观权重;根据每个有效评价指标的主观权重、客观权重以及无纲量化权值,确定矿山地质环境的评价等级。因此,本发明通过将层次分析法和信息熵理论紧密结合,能够使矿山地质环境的评价更加合理化、综合化以及准确化。
Description
技术领域
本发明涉及矿山地质环境评价技术领域,特别涉及一种矿山地质环境综合评价方法及系统。
背景技术
面临矿产开发引发的历史遗留问题以及时代发展对矿山提出的新要求,开展矿山地质环境调查和矿山地质环境评价工作,对矿山地质环境作出充分的评价,以此科学合理的指导矿产开发工作亟不可待。虽然现在存在多种矿山地质环境评价方法,但是普遍存在以下几个方面的问题:(1)现有矿山地质环境评价方法都偏向于主观或者客观方面的单一评价,评价结果缺乏综合性和准确性;(2)现有矿山地质环境评价方法主要针对大范围多类型矿山的综合评价,而缺少对某一矿山的针对性评价;(3)指标选取方面,现有矿山地质环境评价指标的选取是针对区域性的多类型的综合评价,但就单一矿山的评价,指标体系应当更具有针对性和科学性。
发明内容
针对以上缺陷,本发明的目的是提供了一种矿山地质环境综合评价方法及系统,能够使矿山地质环境的评价更加合理化、综合化以及准确化;该方法不仅适用于整体矿山地质区域环境的评价,还适用于单一矿山地质区域环境的评价。
为实现上述目的,本发明提供了如下方案:
一种矿山地质环境综合评价方法,所述矿山地质环境综合评价方法包括:
建立矿山地质环境评价分层结构模型;所述矿山地质环境评价分层结构模型包括三层,分别为目标层、准则层、因素层;所述目标层的目标为矿山地质环境评价等级;所述准则层为与所述目标对应的多个单元,分别为矿山地质环境背景单元、矿山基本概况单元和矿山地质环境问题单元;每个所述单元包括多个评价指标;所有所述评价指标组成所述因素层;
构建每个所述单元的评价指标判断矩阵;
根据所述评价指标判断矩阵,计算所述评价指标判断矩阵对应的所述单元中的每个所述评价指标的权值,并将所述权值小于第一阈值的所述评价指标剔除,确定有效评价指标;所述有效评价指标为所述权值大于或者等于所述第一阈值的评价指标;
对所有所述有效评价指标进行无纲量化处理,得到每个所述有效评价指标的无纲量化权值;
采用层次分析法,计算每个所述有效评价指标的主观权重;
利用信息熵计算每个所述有效评价指标的熵权值;
采用层次分析法对所述熵权值的一致性关系进行检验,确定每个所述有效评价指标的客观权重;
根据每个所述有效评价指标的主观权重、客观权重以及无纲量化权值,确定所述矿山地质环境的评价等级。
可选的,建立所述矿山地质环境评价分层结构模型的原则包括科学性原则、代表性原则、简洁性原则、综合性原则、客观性原则、矿山地质环境定性与定量分析原则。
可选的,所述构建每个所述单元的评价指标判断矩阵,具体包括:
根据1~9标度法,计算每个所述单元中任意两个评价指标的重要性比值;
根据所述重要性比值,确定初步判断矩阵;
对所述初步判断矩阵进行归一化处理,得到每个所述单元的评价指标判断矩阵。
可选的,所述根据所述评价指标判断矩阵,计算所述评价指标判断矩阵对应的所述单元中的每个所述评价指标的权值,并将所述权值小于第一阈值的所述评价指标剔除,确定有效评价指标,具体包括:
根据所述评价指标判断矩阵,计算所述评价指标判断矩阵的特征向量;
计算每个所述评价指标的权值;所述权值为所述评价指标对应的所述特征向量与所有所述特征向量之和的比值;
判断所述评价指标的权值是否小于第一阈值,得到第一判断结果;
若第一判断结果表示所述评价指标的权值小于所述第一阈值,则将所述权值小于所述第一阈值的所述评价指标剔除;
若第一判断结果表示所述评价指标的权值大于或者等于所述第一阈值,则将所述权值大于或者等于所述第一阈值的评价指标保留;所述有效评价指标为所述权值大于或者等于所述第一阈值的评价指标。
可选的,在所述判断所述评价指标的权值是否小于第一阈值,得到第一判断结果之前,所述矿山地质环境综合评价方法还包括:
计算所述评价指标判断矩阵的随机一致性比率;
判断所述评价指标判断矩阵的随机一致性比率是否小于第二阈值,得到第二判断结果;
若所述第二判断结果表示所述评价指标判断矩阵的随机一致性比率小于所述第二阈值,则确定所述评价指标的权值分配合理,执行判断所述评价指标的权值是否小于第一阈值得到第一判断结果步骤;
若所述第二判断结果表示所述评价指标判断矩阵的随机一致性比率大于或者等于所述第二阈值,则确定所述评价指标的权值分配不合理,对所述评价指标的权值进行重新分配。
可选的,所述对所有所述有效评价指标进行无纲量化处理,得到每个所述有效评价指标的无纲量化权值,具体包括:
根据评价指标值与环境值的对应关系,对所述有效评价指标进行分类,确定分类结果;所述分类结果包括第一有效评价指标和第二有效评价指标;所述第一有效评价指标为评价指标值越大环境值越高的评价指标;所述第二有效评价指标为评价指标值越小环境值越高的评价指标;
根据以下公式对所述第一有效评价指标进行无纲量化处理,得到所述第一有效评价指标的无纲量化权值;所述公式为:
Wx(Bi,Xm)=[W(Bi,Xm)-Zmin(Bi,Xm)]/[Zmax(Bi,Xm)-Zmin(Bi,Xm)];
根据以下公式对所述第二有效评价指标进行无纲量化处理,得到所述第二有效评价指标的无纲量化权值;所述公式为:
Wx(Bi,Xm)=[Zmax(Bi,Xm)-W(Bi,Xm)]/[Zmax(Bi,Xm)-Zmin(Bi,Xm)];
其中:W(Bi,Xm)表示有效评价指标r的原始值,Zmax(Bi,Xm)和Zmin(Bi,Xm)分别表示有效评价指标r对应的所述原始值对应的等级划分标准中对应的区间的最大值和最小值;Bi表示第i个单元;i=1,2,3;m表示所述第i个单元中包括的所述有效评价指标的个数;r表示所述有效评价指标,r=1~m。
可选的,所述利用信息熵计算每个所述有效评价指标的熵权值,具体包括:
利用信息熵,计算所述有效评价指标的概率变量值Pir;
根据所述熵值eir和公式gir=1-eir,计算差异性系数gir;
对所述差异性系数gir进行归一化处理,得到所述有效评价指标的熵权值。
可选的,所述采用层次分析法对所述熵权值的一致性关系进行检验,确定每个所述有效评价指标的客观权重,具体包括:
采用层次分析法,构建互补判断矩阵;
计算所述互补判断矩阵的随机一致性比率;
判断所述互补判断矩阵的随机一致性比率是否小于第三阈值,得到第三判断结果;
若所述第三判断结果表示所述互补判断矩阵的随机一致性比率小于所述第三阈值,则确定所述有效评价指标的熵权值分配合理,并将所述熵权值确定为所述有效评价指标的客观权重;
若所述第三判断结果表示所述互补判断矩阵的随机一致性比率大于或者等于所述第三阈值,则确定所述有效评价指标的熵权值分配不合理,对所述有效评价指标的熵权值进行重新分配。
可选的,所述根据每个所述有效评价指标的主观权重、客观权重以及无纲量化权值,确定所述矿山地质环境的评价等级,具体包括:
根据每个所述有效评价指标的主观权重和客观权重,确定所述有效评价指标的综合权重;
根据所述有效评价指标的综合权重、所述无纲量化权值,计算每个所述单元的综合评价函数;
确定所述单元对所述目标层的单元判断矩阵、信息熵评价值;
根据所述单元判断矩阵、所述信息熵评价值,计算每个所述单元的权重;
根据每个所述单元的权重、所述综合评价函数,计算矿山地质环境综合评价函数;
根据所述矿山地质环境综合评价函数,确定所述矿山地质环境的评价等级。
本发明还提供了一种矿山地质环境综合评价系统,所述矿山地质环境综合评价系统包括:
矿山地质环境评价分层结构模型建立模块,用于建立矿山地质环境评价分层结构模型;所述矿山地质环境评价分层结构模型包括三层,分别为目标层、准则层、因素层;所述目标层的目标为矿山地质环境评价等级;所述准则层为与所述目标对应的多个单元,分别为矿山地质环境背景单元、矿山基本概况单元和矿山地质环境问题单元;每个所述单元包括多个评价指标;所有所述评价指标组成所述因素层;
评价指标判断矩阵构建模块,用于构建每个所述单元的评价指标判断矩阵;
有效评价指标确定模块,用于根据所述评价指标判断矩阵,计算所述评价指标判断矩阵对应的所述单元中的每个所述评价指标的权值,并将所述权值小于第一阈值的所述评价指标剔除,确定有效评价指标;所述有效评价指标为所述权值大于或者等于所述第一阈值的评价指标;
无纲量化权值得到模块,用于对所有所述有效评价指标进行无纲量化处理,得到每个所述有效评价指标的无纲量化权值;
主观权重计算模块,用于采用层次分析法,计算每个所述有效评价指标的主观权重;
熵权值计算模块,用于利用信息熵计算每个所述有效评价指标的熵权值;
客观权重确定模块,用于采用层次分析法对所述熵权值的一致性关系进行检验,确定每个所述有效评价指标的客观权重;
评价等级确定模块,用于根据每个所述有效评价指标的主观权重、客观权重以及无纲量化权值,确定所述矿山地质环境的评价等级。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种矿山地质环境综合评价方法及系统,该方法包括建立矿山地质环境评价分层结构模型;所述矿山地质环境评价分层结构模型包括三层,分别为目标层、准则层、因素层;所述目标层的目标为矿山地质环境评价等级;所述准则层为与所述目标对应的多个单元,分别为矿山地质环境背景单元、矿山基本概况单元和矿山地质环境问题单元;每个所述单元包括多个评价指标;所有所述评价指标组成所述因素层;构建每个所述单元的评价指标判断矩阵;根据所述评价指标判断矩阵,计算所述评价指标判断矩阵对应的所述单元中的每个所述评价指标的权值,并将所述权值小于第一阈值的所述评价指标剔除,确定有效评价指标;所述有效评价指标为所述权值大于或者等于所述第一阈值的评价指标;对所有所述有效评价指标进行无纲量化处理,得到每个所述有效评价指标的无纲量化权值;采用层次分析法,计算每个所述有效评价指标的主观权重;利用信息熵计算每个所述有效评价指标的熵权值;采用层次分析法对所述熵权值的一致性关系进行检验,确定每个所述有效评价指标的客观权重;根据每个所述有效评价指标的主观权重、客观权重以及无纲量化权值,确定所述矿山地质环境的评价等级。本发明提供的方法或者系统是将层次分析法和信息熵理论紧密结合,完整体现了矿山地质环境主观意愿和客观实际相融合的过程。尤其在评价指标权重选择衡量上,层次分析法根据主观经验对评价指标权重进行计算,利用信息熵技术对评价指标存在的客观不确定性进行数学描述,两者合理耦合构建趋于更合理化、准确化、综合化的评价指标权重体系。另外,该方法或者系统不仅对目标区的整体矿山地质环境进行有效评价,还可以针对隶属于目标层或者准则层的同一等级矿山地质环境指标进行二次评价。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例矿山地质环境综合评价方法的流程示意图;
图2为本发明矿山地质环境评价分层结构模型示意图;
图3为本发明实施例矿山地质环境综合评价系统的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供了一种矿山地质环境综合评价方法及系统,能够使矿山地质环境的评价更加合理化、综合化以及准确化;该方法不仅适用于整体矿山地质区域环境的评价,还适用于单一矿山地质区域环境的评价。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例矿山地质环境综合评价方法的流程示意图,如图1所示,本发明实施例提供的矿山地质环境综合评价方法具体包括以下几个步骤:
步骤101:建立矿山地质环境评价分层结构模型。
图2为本发明矿山地质环境评价分层结构模型示意图,如图2所示,该矿山地质环境评价分层结构模型包括三层,第一层为目标层,第二层为准则层,第三层为因素层;所述目标层的目标为矿山地质环境评价等级;所述准则层为与所述目标层对应的多个单元,分别为矿山地质环境背景单元、矿山基本概况单元和矿山地质环境问题单元;每个单元包括多个评价指标;所有评价指标组成所述因素层。
其中,矿山地质环境背景单元包括七个评价指标,分别为:水文地质、工程地质、环境地质、土地利用现状、地形地貌、岩性与地质构造以及其他人类工程影响。
所述矿山基本概况单元包括七个评价指标,分别为矿山面积、矿山类型、生产能力、采矿方式、开采深度、采空区面积以及矿区经济效益。
矿山地质环境问题单元包括六个评价指标,分别为矿山地质灾害、地下含水层破坏、地形地貌景观破坏、压占破坏土地、固体废弃物排放以及地质环境综合治理。
目标层选取:矿山地质环境调查和评价的主要目的是对矿山地质环境和矿山地质环境问题及危害进行充分调查,分析评价矿产资源开发对地质环境的影响,以此提出具体的矿山地质环境保护与治理恢复对策措施,合理科学的开发矿产资源,保护矿山地质环境,实施矿山地质环境监督管理等提供依据,按照其调查评价主要目的,矿山地质环境调查主要内容包括矿山基本情况,矿山开发引起的地质环境问题、影响及危害,矿山地质环境治理措施及效果调查等三大类,所以,在本发明的矿山地质环境指标体系中目标层为:矿山地质环境评价等级。
矿山地质环境评价分层结构模型的构建原则:科学性原则、代表性原则、简洁性原则、综合性原则、客观性原则,除此之外,一定要与矿山的实际情况进行紧密结合,充分调动指标选择和指标定量的灵活性,严格遵循矿山地质环境定性与定量分析的基本原则。
步骤102:构建每个所述单元的评价指标判断矩阵。
步骤103:根据所述评价指标判断矩阵,计算所述评价指标判断矩阵对应的所述单元中的每个所述评价指标的权值,并将所述权值小于第一阈值的所述评价指标剔除,确定有效评价指标;所述有效评价指标为所述权值大于或者等于所述第一阈值的评价指标。
步骤104:对所有所述有效评价指标进行无纲量化处理,得到每个所述有效评价指标的无纲量化权值。
步骤105:采用层次分析法,计算每个所述有效评价指标的主观权重。
步骤106:利用信息熵计算每个所述有效评价指标的熵权值。
步骤107:采用层次分析法对所述熵权值的一致性关系进行检验,确定每个所述有效评价指标的客观权重。
步骤108:根据每个所述有效评价指标的主观权重、客观权重以及无纲量化权值,确定所述矿山地质环境的评价等级。
其中,层次分析法(AHP)作为一种经典的多目标综合评价方法,在构建多类型评价指标体系和评价模型中发挥了重要作用,但评价方法偏于主观性,缺少对评价指标不确定性因素影响的考虑。信息熵技术利用模糊数学和概率论原理,充分描述了评价指标存在的不确定性,并利用数学方法对指标不确定性进行了解释和剔除。
其中,步骤102具体包括:
步骤1021:根据1~9标度法,计算每个所述单元中任意两个评价指标的重要性比值;具体为:针对准则层中各个单元实现目标和专家经验值判断,分别在各个单元内部,根据1~9标度法对每个单元Bi中包括的评价指标Xij重要性进行确定。例如两个评价指标为Xik、Xij,若这两个评价指标重要性相同,则重要性比值取1,若两个评价指标重要性不同,则重要性比值取Ckj(Cjk=1/Ckj)。
步骤1022:根据所述重要性比值,确定初步判断矩阵;所述初步判断矩阵为:C(i,n)=(Ckj)n×n。
步骤1023:对所述初步判断矩阵进行归一化处理,得到每个所述单元的评价指标判断矩阵。所述评价指标判断矩阵为:
其中:Ckj表示两个评价指标所体现的对于单元Bi的重要性比值;S(i,n)表示归一化判断矩阵(评价指标判断矩阵),其中i表示评价指标对应的准则层中的单元(i=1~3),n表示准则层中的单元所包括的评价指标的个数,k=1~n;j=1~n。
步骤103具体包括:
其中:j表示第j个评价指标,j=1~n。
③由于评价指标的权值的选择主要是基于专家主观的经验评判,所以难免会存在经验判断上的不合理性,可能会造成评价指标的权值分配存在不合理性或者脱离实际性,为了消除权值分配可能存在的不平衡性,必须对权值的分配合理性进行判断。
权值分配合理性以评价指标判断矩阵所表现的随机一致性比率CR确定,CR=CI/RI;其中CI表示评价指标判断矩阵的一致性指标,CI=(λmax-n)/(n-1),λmax表示评价指标判断矩阵的最大特征根;RI表示平均随机一致性指标值,见表1。其中:初步判断矩阵C(i,n)=(Ckj)n×n,Q为权值集合
表1平均随机一致性指标RI值
具体为:
计算评价指标判断矩阵的随机一致性比率;
判断评价指标判断矩阵的随机一致性比率是否小于第二阈值;优选的,第二阈值为0.1。
若是,确定评价指标的权值分配合理,执行判断评价指标的权值是否小于第一阈值步骤;若否,则确定评价指标的权值分配不合理,对评价指标的权值进行重新分配。
判断评价指标的权值是否小于第一阈值;优选的,第一阈值为0.1。
若是,则将权值小于所述第一阈值的所述评价指标剔除;若否,则将权值大于或者等于所述第一阈值的评价指标保留;所述有效评价指标为所述权值大于或者等于所述第一阈值的评价指标。
步骤104具体包括:
根据所选择的矿山地质环境评价指标可以看出,评价指标属性不尽相同,为了建立统一的评价体系,就必须对不同属性的评价指标进行无纲量化,建立统一无纲量化的越大环境指数越高的相对指标值{Wx(Bi,Xm)}。矿山地质环境评价指标所产生的影响值可以概括分为两类:评价指标值越大对应环境值越高的第一有效评价指标和评价指标值越小对应环境值越高的第二有效评价指标,对于两种不同的指标分类,无纲量化的基本方法是基于有效评价指标所对应的等级划分标准来对该有效评价指标进行统一无量纲。
根据以下公式对所述第一有效评价指标进行无纲量化处理,得到所述第一有效评价指标的无纲量化权值;所述公式为:
Wx(Bi,Xm)=[W(Bi,Xm)-Zmin(Bi,Xm)]/[Zmax(Bi,Xm)-Zmin(Bi,Xm)]。
根据以下公式对所述第二有效评价指标进行无纲量化处理,得到所述第二有效评价指标的无纲量化权值;所述公式为:
Wx(Bi,Xm)=[Zmax(Bi,Xm)-W(Bi,Xm)]/[Zmax(Bi,Xm)-Zmin(Bi,Xm)]。
其中:W(Bi,Xm)表示有效评价指标r的原始值,Zmax(Bi,Xm)和Zmin(Bi,Xm)分别表示有效评价指标r对应的所述原始值对应的等级划分标准中对应的区间的最大值和最小值;Bi表示第i个单元;i=1,2,3;m表示所述第i个单元中包括的所述有效评价指标的个数;r表示所述有效评价指标,r=1~m。
步骤105具体包括:
步骤106具体包括
c.计算有效评价指标的的熵权值
根据公式gir=1-eir,计算差异性系数gir,并对差异性系统gir进行归一化处理,得到有效评价指标的熵权值其中i表示有效评价指标对应的准则层中的单元(i=1~3),r表示有效评价指标,r=1~m。
步骤107具体包括:
利用信息熵确定的各有效评价指标的熵权值为独立量,为了真正反映信息熵确定的有效评价指标权值之间的一致性关系和确权的有效性,必须对初步熵权值进行一致性检验,本发明利用层次分析法对熵权值的一致性关系进行检验。
a.构建互补判断矩阵Ai
其中应当满足如下条件:0≤arl≤1;arl+alr=1;arl=0.5表示有效评价指标指标r和有效评价指标l的重要性相同,arl≥0.5表示有效评价指标指标r和有效评价指标l重要。
b.判断互补矩阵的一致性关系
基于构建的互补判断矩阵,计算互补判断矩阵的随机一致性比率CR,随机一致性比率CR<0.1,则表示该互补判断矩阵有效,各有效评价指标的熵权值分配合理,反之,则应当对各有效评价指标的熵权值进行重新分配。
c.互补判断矩阵修正
其中:CIC(ni)称为一致性指标系数,即该值系数越小则互补矩阵的一致性关系越高,以此重新分配各有效评价指标的熵权值。
客观权值体现了有效评价指标表现出的不确定性和客观效应的事实,本发明利用信息熵对有效评价指标的客观权值进行确定,即将熵权值定义为客观权重,通过一致性检验合格后的熵权值定义为最终的客观权重。
步骤108具体包括:
其中:a表示优先系数,即判断某一有效评价指标的客观权重和主观权重的优先考虑度,为了综合考虑主、客观权重的合理性,一般a选择0.5。
第二步:根据所述有效评价指标的综合权重、所述无纲量化权值,计算每个所述单元的综合评价函数,所述单元的综合评价函数:
其中:Z(i,r)值越大,表示单元Bi所对应的环境评价越高,例如矿山地质环境背景单元B1中土地利用现状评价结果值相对较高,则说明评价目标区的土地利用现状较为合理,属于矿山地质环境评价整体结果的加分项,而如果目标区的矿山地质环境问题中评价指标综合评价较低,说明该地区矿山地质灾害较严重,为该目标区存在的主要矿山地质环境问题。
第三步:构建准则层面对目标层的判断矩阵AS和信息熵评价值,从而获取针对目标层的各权重{qs(i)|i=1~3};
第四步根据每个所述单元的权重、所述综合评价函数,计算矿山地质环境综合评价函数;矿山地质环境综合评价函数为:
其中:S值越大,表示矿山地质环境越好,矿山地质环境问题越少,例如矿山地质环境问题单元B3评价指数较低,说明该地区矿山地质环境问题严重,并且矿山地质环境问题已经直接影响到了目标区整体矿山地质环境的发展。
第五步:根据所述矿山地质环境综合评价函数,确定所述矿山地质环境的评价等级。
本发明实施例提供的方法是将层次分析法和信息熵理论紧密结合,完整体现了矿山地质环境主观意愿和客观实际相融合的过程。尤其在评价指标权重选择衡量上,层次分析法根据主观经验对评价指标权重进行计算,利用信息熵技术对评价指标存在的客观不确定性进行数学描述,两者合理耦合构建趋于更合理化的评价指标权重体系。该方法不仅对目标区的整体矿山地质环境进行有效评价,还可以针对隶属于目标层或者准则层的同一等级某方面矿山地质环境指标进行二次评价。
为实现上述目的,本发明还提供了一种矿山地质环境综合评价系统。
图3为本发明实施例矿山地质环境综合评价系统的的结构示意图,如图3所示,本发明实施例提供的所述矿山地质环境综合评价系统包括:
矿山地质环境评价分层结构模型建立模块301,用于建立矿山地质环境评价分层结构模型;所述矿山地质环境评价分层结构模型包括三层,分别为目标层、准则层、因素层;所述目标层的目标为矿山地质环境评价等级;所述准则层为与所述目标对应的多个单元,分别为矿山地质环境背景单元、矿山基本概况单元和矿山地质环境问题单元;每个所述单元包括多个评价指标;所有所述评价指标组成所述因素层。
评价指标判断矩阵构建模块302,用于构建每个所述单元的评价指标判断矩阵。
有效评价指标确定模块303,用于根据所述评价指标判断矩阵,计算所述评价指标判断矩阵对应的所述单元中的每个所述评价指标的权值,并将所述权值小于第一阈值的所述评价指标剔除,确定有效评价指标;所述有效评价指标为所述权值大于或者等于所述第一阈值的评价指标。
无纲量化权值得到模块304,用于对所有所述有效评价指标进行无纲量化处理,得到每个所述有效评价指标的无纲量化权值。
主观权重计算模块305,用于采用层次分析法,计算每个所述有效评价指标的主观权重。
熵权值计算模块306,用于利用信息熵计算每个所述有效评价指标的熵权值。
客观权重确定模块307,用于采用层次分析法对所述熵权值的一致性关系进行检验,确定每个所述有效评价指标的客观权重。
评价等级确定模块308,用于根据每个所述有效评价指标的主观权重、客观权重以及无纲量化权值,确定所述矿山地质环境的评价等级。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (9)
1.一种矿山地质环境综合评价方法,其特征在于,所述矿山地质环境综合评价方法包括:
建立矿山地质环境评价分层结构模型;所述矿山地质环境评价分层结构模型包括三层,分别为目标层、准则层、因素层;所述目标层的目标为矿山地质环境评价等级;所述准则层为与所述目标对应的多个单元,分别为矿山地质环境背景单元、矿山基本概况单元和矿山地质环境问题单元;每个所述单元包括多个评价指标;所有所述评价指标组成所述因素层;
构建每个所述单元的评价指标判断矩阵;
根据所述评价指标判断矩阵,计算所述评价指标判断矩阵对应的所述单元中的每个所述评价指标的权值,并将所述权值小于第一阈值的所述评价指标剔除,确定有效评价指标;所述有效评价指标为所述权值大于或者等于所述第一阈值的评价指标;
对所有所述有效评价指标进行无纲量化处理,得到每个所述有效评价指标的无纲量化权值;
采用层次分析法,计算每个所述有效评价指标的主观权重;
利用信息熵计算每个所述有效评价指标的熵权值;
采用层次分析法对所述熵权值的一致性关系进行检验,确定每个所述有效评价指标的客观权重;
根据每个所述有效评价指标的主观权重、客观权重以及无纲量化权值,确定所述矿山地质环境的评价等级;
所述根据每个所述有效评价指标的主观权重、客观权重以及无纲量化权值,确定所述矿山地质环境的评价等级,具体包括:
根据每个所述有效评价指标的主观权重和客观权重,确定所述有效评价指标的综合权重;利用层次分析法计算的主观权重和信息熵确定的客观权重获取有效评价指标的综合权重qir,即:a表示优先系数,即判断某一有效评价指标的客观权重和主观权重的优先考虑度;
根据所述有效评价指标的综合权重、所述无纲量化权值,计算每个所述单元的综合评价函数;确定所述单元对所述目标层的单元判断矩阵、信息熵评价值;根据所述单元判断矩阵、所述信息熵评价值,计算每个所述单元的权重;根据每个所述单元的权重、所述综合评价函数,计算矿山地质环境综合评价函数;根据所述矿山地质环境综合评价函数,确定所述矿山地质环境的评价等级。
2.根据权利要求1所述的矿山地质环境综合评价方法,其特征在于,建立所述矿山地质环境评价分层结构模型的原则包括科学性原则、代表性原则、简洁性原则、综合性原则、客观性原则、矿山地质环境定性与定量分析原则。
3.根据权利要求1所述的矿山地质环境综合评价方法,其特征在于,所述构建每个所述单元的评价指标判断矩阵,具体包括:
根据1~9标度法,计算每个所述单元中任意两个评价指标的重要性比值;
根据所述重要性比值,确定初步判断矩阵;
对所述初步判断矩阵进行归一化处理,得到每个所述单元的评价指标判断矩阵。
4.根据权利要求3所述的矿山地质环境综合评价方法,其特征在于,所述根据所述评价指标判断矩阵,计算所述评价指标判断矩阵对应的所述单元中的每个所述评价指标的权值,并将所述权值小于第一阈值的所述评价指标剔除,确定有效评价指标,具体包括:
根据所述评价指标判断矩阵,计算所述评价指标判断矩阵的特征向量;
计算每个所述评价指标的权值;所述权值为所述评价指标对应的所述特征向量与所有所述特征向量之和的比值;
判断所述评价指标的权值是否小于第一阈值,得到第一判断结果;
若第一判断结果表示所述评价指标的权值小于所述第一阈值,则将所述权值小于所述第一阈值的所述评价指标剔除;
若第一判断结果表示所述评价指标的权值大于或者等于所述第一阈值,则将所述权值大于或者等于所述第一阈值的评价指标保留;所述有效评价指标为所述权值大于或者等于所述第一阈值的评价指标。
5.根据权利要求4所述的矿山地质环境综合评价方法,其特征在于,在所述判断所述评价指标的权值是否小于第一阈值,得到第一判断结果之前,所述矿山地质环境综合评价方法还包括:
计算所述评价指标判断矩阵的随机一致性比率;
判断所述评价指标判断矩阵的随机一致性比率是否小于第二阈值,得到第二判断结果;
若所述第二判断结果表示所述评价指标判断矩阵的随机一致性比率小于所述第二阈值,则确定所述评价指标的权值分配合理,执行判断所述评价指标的权值是否小于第一阈值得到第一判断结果步骤;
若所述第二判断结果表示所述评价指标判断矩阵的随机一致性比率大于或者等于所述第二阈值,则确定所述评价指标的权值分配不合理,对所述评价指标的权值进行重新分配。
6.根据权利要求1所述的矿山地质环境综合评价方法,其特征在于,所述对所有所述有效评价指标进行无纲量化处理,得到每个所述有效评价指标的无纲量化权值,具体包括:
根据评价指标值与环境值的对应关系,对所述有效评价指标进行分类,确定分类结果;所述分类结果包括第一有效评价指标和第二有效评价指标;所述第一有效评价指标为评价指标值越大环境值越高的评价指标;所述第二有效评价指标为评价指标值越小环境值越高的评价指标;
根据以下公式对所述第一有效评价指标进行无纲量化处理,得到所述第一有效评价指标的无纲量化权值;所述公式为:
Wx(Bi,Xm)=[W(Bi,Xm)-Zmin(Bi,Xm)]/[Zmax(Bi,Xm)-Zmin(Bi,Xm)];
根据以下公式对所述第二有效评价指标进行无纲量化处理,得到所述第二有效评价指标的无纲量化权值;所述公式为:
Wx(Bi,Xm)=[Zmax(Bi,Xm)-W(Bi,Xm)]/[Zmax(Bi,Xm)-Zmin(Bi,Xm)];
其中:W(Bi,Xm)表示有效评价指标r的原始值,Zmax(Bi,Xm)和Zmin(Bi,Xm)分别表示有效评价指标r对应的所述原始值对应的等级划分标准中对应的区间的最大值和最小值;Bi表示第i个单元;i=1,2,3;m表示所述第i个单元中包括的所述有效评价指标的个数;r表示所述有效评价指标,r=1~m。
8.根据权利要求1所述的矿山地质环境综合评价方法,其特征在于,所述采用层次分析法对所述熵权值的一致性关系进行检验,确定每个所述有效评价指标的客观权重,具体包括:
采用层次分析法,构建互补判断矩阵;
计算所述互补判断矩阵的随机一致性比率;
判断所述互补判断矩阵的随机一致性比率是否小于第三阈值,得到第三判断结果;
若所述第三判断结果表示所述互补判断矩阵的随机一致性比率小于所述第三阈值,则确定所述有效评价指标的熵权值分配合理,并将所述熵权值确定为所述有效评价指标的客观权重;
若所述第三判断结果表示所述互补判断矩阵的随机一致性比率大于或者等于所述第三阈值,则确定所述有效评价指标的熵权值分配不合理,对所述有效评价指标的熵权值进行重新分配。
9.一种矿山地质环境综合评价系统,其特征在于,所述矿山地质环境综合评价系统包括:
矿山地质环境评价分层结构模型建立模块,用于建立矿山地质环境评价分层结构模型;所述矿山地质环境评价分层结构模型包括三层,分别为目标层、准则层、因素层;所述目标层的目标为矿山地质环境评价等级;所述准则层为与所述目标对应的多个单元,分别为矿山地质环境背景单元、矿山基本概况单元和矿山地质环境问题单元;每个所述单元包括多个评价指标;所有所述评价指标组成所述因素层;
评价指标判断矩阵构建模块,用于构建每个所述单元的评价指标判断矩阵;
有效评价指标确定模块,用于根据所述评价指标判断矩阵,计算所述评价指标判断矩阵对应的所述单元中的每个所述评价指标的权值,并将所述权值小于第一阈值的所述评价指标剔除,确定有效评价指标;所述有效评价指标为所述权值大于或者等于所述第一阈值的评价指标;
无纲量化权值得到模块,用于对所有所述有效评价指标进行无纲量化处理,得到每个所述有效评价指标的无纲量化权值;
主观权重计算模块,用于采用层次分析法,计算每个所述有效评价指标的主观权重;
熵权值计算模块,用于利用信息熵计算每个所述有效评价指标的熵权值;
客观权重确定模块,用于采用层次分析法对所述熵权值的一致性关系进行检验,确定每个所述有效评价指标的客观权重;
评价等级确定模块,用于根据每个所述有效评价指标的主观权重、客观权重以及无纲量化权值,确定所述矿山地质环境的评价等级;
所述根据每个所述有效评价指标的主观权重、客观权重以及无纲量化权值,确定所述矿山地质环境的评价等级,具体包括:
根据每个所述有效评价指标的主观权重和客观权重,确定所述有效评价指标的综合权重;利用层次分析法计算的主观权重和信息熵确定的客观权重获取有效评价指标的综合权重qir,即:a表示优先系数,即判断某一有效评价指标的客观权重和主观权重的优先考虑度;
根据所述有效评价指标的综合权重、所述无纲量化权值,计算每个所述单元的综合评价函数;确定所述单元对所述目标层的单元判断矩阵、信息熵评价值;根据所述单元判断矩阵、所述信息熵评价值,计算每个所述单元的权重;根据每个所述单元的权重、所述综合评价函数,计算矿山地质环境综合评价函数;根据所述矿山地质环境综合评价函数,确定所述矿山地质环境的评价等级。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810300100.4A CN108710984B (zh) | 2018-04-04 | 2018-04-04 | 一种矿山地质环境综合评价方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810300100.4A CN108710984B (zh) | 2018-04-04 | 2018-04-04 | 一种矿山地质环境综合评价方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108710984A CN108710984A (zh) | 2018-10-26 |
CN108710984B true CN108710984B (zh) | 2021-06-08 |
Family
ID=63867043
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810300100.4A Active CN108710984B (zh) | 2018-04-04 | 2018-04-04 | 一种矿山地质环境综合评价方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108710984B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109858788A (zh) * | 2019-01-17 | 2019-06-07 | 北京师范大学 | 一种基于生态型地区的地质资源评价方法 |
CN109784730A (zh) * | 2019-01-17 | 2019-05-21 | 北京师范大学 | 适用于生态型地区的矿产资源的评价方法 |
CN110080818B (zh) * | 2019-06-11 | 2021-04-20 | 中国矿业大学 | 一种矿山遗留矿柱分级分类充填回收方法 |
CN110210790A (zh) * | 2019-06-18 | 2019-09-06 | 北京蓝尊科技有限公司 | 一种绿色矿山建设监管方法及系统 |
CN110728402B (zh) * | 2019-10-10 | 2023-11-03 | 辽宁工程技术大学 | 一种基于地质标本分析的区域矿产资源预测系统 |
CN111059998B (zh) * | 2019-12-31 | 2020-11-13 | 中国地质大学(北京) | 一种基于高分辨率的时序InSAR形变监测方法及系统 |
CN111260187A (zh) * | 2020-01-08 | 2020-06-09 | 长春工程学院 | 一种智能矿山地质环境信息评价系统及评价方法 |
CN111476472A (zh) * | 2020-01-14 | 2020-07-31 | 中化地质矿山总局地质研究院 | 硫铁矿山地质环境评价方法 |
CN112036732B (zh) * | 2020-08-27 | 2023-05-12 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种山区环境质量双标双控的评估方法 |
CN112001568B (zh) * | 2020-09-11 | 2022-05-20 | 新疆大学 | 高海拔高寒金属矿开采风钻作业效率影响因素的评估方法 |
CN112305634B (zh) * | 2020-10-31 | 2023-11-17 | 中国海洋石油集团有限公司 | 一种基于成因及形态分析的湖底扇含砂率预测方法 |
CN112598304A (zh) * | 2020-12-28 | 2021-04-02 | 中国能源建设集团陕西省电力设计院有限公司 | 一种基于联系熵的电力工程建设采空区评价方法 |
CN112987079A (zh) * | 2021-03-15 | 2021-06-18 | 河北省地质环境监测院(河北省地质灾害应急技术中心、河北省地质环境监测总站) | 一种矿山地质环境水土检测指标实时传输分析系统 |
CN113780719A (zh) * | 2021-07-26 | 2021-12-10 | 重庆地质矿产研究院 | 一种矿山地质环境综合评价方法 |
CN114646750B (zh) * | 2022-05-13 | 2022-10-14 | 中国地质环境监测院(自然资源部地质灾害技术指导中心) | 一种土壤环境实时监测装置 |
CN115809834B (zh) * | 2023-02-03 | 2023-04-25 | 广东智环创新环境科技有限公司 | 一种基于环境影响评价数据的生态环境分析系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105205328A (zh) * | 2015-09-28 | 2015-12-30 | 辽宁工程技术大学 | 一种基于广义分维和特征参数的地质动力环境评价方法 |
CN105469196A (zh) * | 2015-11-18 | 2016-04-06 | 山东科技大学 | 一种矿井建设项目过程后评价的综合评价方法及系统 |
CN106021873A (zh) * | 2016-05-11 | 2016-10-12 | 吉林大学 | 矿山地质环境评价方法 |
-
2018
- 2018-04-04 CN CN201810300100.4A patent/CN108710984B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105205328A (zh) * | 2015-09-28 | 2015-12-30 | 辽宁工程技术大学 | 一种基于广义分维和特征参数的地质动力环境评价方法 |
CN105469196A (zh) * | 2015-11-18 | 2016-04-06 | 山东科技大学 | 一种矿井建设项目过程后评价的综合评价方法及系统 |
CN106021873A (zh) * | 2016-05-11 | 2016-10-12 | 吉林大学 | 矿山地质环境评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108710984A (zh) | 2018-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108710984B (zh) | 一种矿山地质环境综合评价方法及系统 | |
CN109034582B (zh) | 基于云模型和组合赋权的隧道穿越断层突水突泥风险评价方法 | |
CN108280553B (zh) | 基于gis-神经网络集成的山洪灾害风险区划及预测方法 | |
Zhang et al. | Developing a cloud model based risk assessment methodology for tunnel-induced damage to existing pipelines | |
CN104346538B (zh) | 基于三种灾情因子控制的地震灾害评估方法 | |
CN106997509B (zh) | 一种不确定信息融合的应急物资分布需求预测方法 | |
CN110610285A (zh) | 一种地下金属矿采空区危险度分级评价方法 | |
CN108256141B (zh) | 一种基于Copula理论的主余震联合易损性分析方法 | |
Zhang et al. | Risk-based estimate for operational safety in complex projects under uncertainty | |
CN104951588A (zh) | 一种矿井通风系统辅助设计方法 | |
CN110020749B (zh) | 一种基于微震分形预测冲击危险等级的综合智能预测方法 | |
KR102379472B1 (ko) | 시공간 피해특성을 고려한 다종 데이터 융합 분석방법 | |
CN111242499A (zh) | 基于层次-可拓分析的既有隧道衬砌结构病害评价方法 | |
CN112101811A (zh) | 一种供水管网爆管风险预测方法及系统 | |
CN111882198A (zh) | 一种项目绩效评价方法及系统 | |
CN113887872A (zh) | 一种基于态势评估技术的交通风险分级方法 | |
CN106970986A (zh) | 基于深度学习的城市内涝影响程度挖掘方法及系统 | |
CN108805471A (zh) | 基于复合系统作用关系分析的水资源承载能力评价方法 | |
CN111476472A (zh) | 硫铁矿山地质环境评价方法 | |
CN113505978A (zh) | 一种不同形式城市社区的防灾功能评估方法和装置 | |
CN109118004A (zh) | 一种工程构筑选址适宜区预测方法 | |
Chen et al. | Genetic programming for predicting aseismic abilities of school buildings | |
CN103971300A (zh) | 基于gis空间信息的变电站站址定位方法 | |
CN115906663A (zh) | 房屋安全评估模型建立方法、评估方法、服务器及系统 | |
Xin et al. | Review of fragility analyses for major building types in China with new implications for intensity–PGA relation development |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |