CN110610285A - 一种地下金属矿采空区危险度分级评价方法 - Google Patents
一种地下金属矿采空区危险度分级评价方法 Download PDFInfo
- Publication number
- CN110610285A CN110610285A CN201910659265.5A CN201910659265A CN110610285A CN 110610285 A CN110610285 A CN 110610285A CN 201910659265 A CN201910659265 A CN 201910659265A CN 110610285 A CN110610285 A CN 110610285A
- Authority
- CN
- China
- Prior art keywords
- index
- evaluation
- goaf
- risk
- weight
- 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.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 128
- 239000002184 metal Substances 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 88
- 239000011159 matrix material Substances 0.000 claims description 33
- 238000004364 calculation method Methods 0.000 claims description 29
- 239000013598 vector Substances 0.000 claims description 21
- 239000011435 rock Substances 0.000 claims description 16
- 238000010606 normalization Methods 0.000 claims description 15
- 238000001514 detection method Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000009933 burial Methods 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 150000007524 organic acids Chemical class 0.000 claims description 3
- 238000012854 evaluation process Methods 0.000 abstract 1
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000005314 correlation function Methods 0.000 description 2
- 238000013210 evaluation model Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 231100000817 safety factor Toxicity 0.000 description 2
- 238000012163 sequencing technique Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
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/0635—Risk analysis of enterprise or organisation activities
-
- 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)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Primary Health Care (AREA)
- Animal Husbandry (AREA)
- General Health & Medical Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Marine Sciences & Fisheries (AREA)
- Agronomy & Crop Science (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种地下金属矿采空区危险度分级评价方法,该方法首先选取16个指标构建了地下金属矿采空区危险度分级评价指标体系,将危险度分为4个等级并确定了分级标准;然后采用AHP、熵权法和博弈论组合赋权法确定各指标的综合权重;并建立采空区危险度评价的物元可拓模型;将各评价指标的权重值代入物元可拓模型中,得到隶属于各危险度等级的各指标关联度和综合关联度值;最后根据最大关联度原则确定出危险度等级。本发明能够全面、客观、准确地评价采空区危险度等级,评价过程清晰、评价结果客观,可为地下金属矿采空区危险度分级评价提供一种新方法。
Description
技术领域
本发明属于地下矿山安全领域,具体涉及一种地下金属矿采空区危险度分级评价方法。
背景技术
采空区是地下金属矿山开采的重大危险源之一,我国许多大型地下金属矿山,都存在着大量的采空区。大量采空区的存在使得矿山开采条件恶化和隐患资源量增加,极易引发井下采空区的大面积冒落、地表塌陷和巷道变形等重大灾害,造成严重的人员伤亡和设备损坏。因此,在实施采空区精细探测的基础上,开展采空区危险度分级评价,对于制定采空区处理措施,保障矿山安全生产具有重要现实意义。
在已公开的成果中,层次分析法、模糊数学法、数值计算法、BP神经网络法、灰色理论等方法已被应用于采空区的危险度评价中。然而由于每个采空区的赋存位置、空间形态、体积大小、工程地质条件等情况均存在差异,涉及的因素也很多,给采空区的危险度分级评价工作带来了难度。此外,各种评价方法也有自身的局限性,导致评价结果不能准确、客观、全面的反映出采空区的危险度等级。例如采用层次分析法确定权重时,容易受主观因素影响;采用熵权法确定权重时,又过于依赖客观信息而忽略指标本身的属性;BP神经网络法易陷入局部极小点和收敛速度慢;可拓学评价中往往只使用一种方法来确定权重,导致评价结果与实际情况不符。因此,采用某一种评价方法难以对采空区危险度等级进行准确、快速的评价。
鉴于此,本发明以采空区探测结果为基础,建立采空区危险度评价指标体系,然后采用层次分析法、熵权法以及基于博弈论的组合赋权法来确定评价指标的组合权重,同时结合物元可拓理论,建立基于组合赋权的物元可拓评价模型,从而可以获得准确、合理的采空区危险度评价结果,进而提出了一种地下矿采空区危险度分级评价方法。
发明内容
本申请旨在至少解决现有技术中存在的技术问题之一。为此,本发明的目的之一在于提供一种基于组合赋权的物元可拓评价模型的地下金属矿采空区危险度分级评价方法,该方法能够客观、准确、快速地评价地下金属矿采空区的危险度等级。
为解决上述技术问题,本申请采用如下技术方案:
一种地下金属矿采空区危险度分级评价方法,包括以下步骤:
步骤1:对地下金属矿采空区进行现场探测,并收集采空区相关的工程地质资料;
步骤2:选取n个评价指标,建立采空区危险度分级评价指标体系,并确定危险度等级的分级标准;
步骤3:运用层次分析法(AHP)计算出各评价指标的主观权重Wa,并进行总体一致性检验;
采用熵权法计算采空区危险度各评价指标的信息熵Ei,得到指标的客观权重Wb;
根据AHP、熵权法得到的权重,采用博弈论组合赋权法,计算出各评价指标的综合权重Wc。
步骤4:采用可拓学理论,建立采空区危险度评价的物元可拓模型;
步骤5:将组合赋权法得到的各指标综合权重值代入至物元可拓模型中,对各项指标进行关联度运算,得到关于分级等级的关联度计算结果;
步骤6:根据最大关联度原则,判定采空区的危险度等级。
进一步的,所述步骤2中,选定采空区类型、地质构造、岩体结构、围岩支护、顶板特征、地表特征、水文因素、相邻空区分布、工程布置共9个定性指标,以及采空区规模、埋深、高跨比、暴露面积、矿柱安全系数、岩石抗压强度和岩石质量标准共7个定量指标,合计16个指标建立采空区危险度分级评价指标体系,并将危险度等级从低至高分为差(Ⅰ级)、中(Ⅱ级)、良(Ⅲ级)和优(Ⅳ级)四个等级。
进一步的,所述步骤3中,根据AHP法计算指标的主观权重,并进行总体一致性检验的过程具体包括以下步骤:
(1)建立层次结构模型。以采空区危险度为目标层,以各指标为指标层,建立指标体系的层次结构模型;
(2)构造判断矩阵。采用基于Satty提出的1~9标度法,对各指标对采空区危险度的相对重要性进行比较,用数值表示出来,并写成矩阵形式,即构成判断矩阵X:
X=(xij)m×n (1)
式(1)中,xij表示指标xi与指标xj重要性比较的取值,i,j=1,2,...,n;
(3)计算判断矩阵特征向量W和判断矩阵的最大特征根λmax:
W=(w1,w2,ΛΛ,wn)T (2)
(4)对判断矩阵X进行一致性检验。计算判断矩阵的一致性指标和一致性比率:
式中,C.I.为一致性指标,C.R.为一致性比率,R.I为随机一致性指标,n为判断矩阵的阶数。
当时,即认为判断矩阵具有满意的一致性。
(5)层次总排序及一致性检验。参照步骤(2)-(4),计算出指标层对于采空区危险度相对重要性的排序,即权重Wa,并进行总排序的一致性检验。
进一步的,所述步骤4中,采用熵权法计算采空区危险度各评价指标的客观权重的过程包括以下步骤:
(1)根据初始数据获得共m个评价样本(采空区),n个评价指标,构建初始数据矩阵R:
R=(rij)m×n,i=1,2,...,m;j=1,2,...,n; (6)
式(6)中,rij为第j个指标下第i个采空区的评价值。
(2)对指标进行数据无量纲和归一化处理:
对于正向指标(数值越高越优),归一化公式为:
对于负向的指标(数值越低越优),归一化公式为:
(3)根据归一化处理后得到的数据r′ij,计算出第j指标下第i个采空区占该指标的比重:
(4)计算第j项指标的信息熵Ei,计算公式为:
(5)根据公式(11)计算出各指标的权重值,得到指标熵权值的排序向量Wb
进一步的,所述步骤5中,采用博弈论组合赋权法计算采空区危险度各评价指标的综合权重的过程包括以下步骤:
(1)根据L种方法对同一指标获得的L个指标权重值,构造权重向量uk=[uk1,uk2,...ukn],k=1,2,...,L,将L个权重向量任意线性组合记为:
式中:αk为线性组合系数,αk>0;u为可能的权重向量集合。
L在这里是指采用博弈论组合赋权方法时,可能采用了多种(L种)方法得到了同一指标的L个权重值。本专利中L为2,即:采用了AHP和熵权法2种方法得到了同一个指标的2个权重值,即Wa和Wb)。
(2)为了使可能向量集u和各个uk的离差值达到最小,对组合系数αk进行优化,进而得到uk的最优解uk *,计算公式如下:
计算可得:则对应的线性方程组为:
(3)对式(14)进行求解,得到组合系数[α1,α2,…,αL],然后对其进行归一化处理,计算公式如下:
(4)求出采空区危险度各评价指标的综合权重,计算公式如下:
进一步的,所述步骤6中,构建采空区危险度评价的物元可拓模型的过程包括以下步骤:
(1)物元模型
假设待评价的采空区评价等级为H,评价指标为Ci,各项指标的特征向量值为Vi,以有序三元组R={H,Ci,Vi}构建待评价采空区危险度的物元,如下式所示:
(2)构建评价对象的经典域物元模型,表达式如下:
式中:Rj为第j个同征物元;Hj为评价对象的第j个危险度等级;Ci为评价对象的第i个评价指标,i=1,2,…,n;Vji为危险度等级j第i个评价指标Ci标准化后的经典域;aji和bji分别为Vji的取值界限。
(3)构建评价对象的节域物元模型,计算表达式如下:
式中:HP为评价对象危险度等级的全体;VPi为P关于评价指标Ci所取的值域,即HP节域;api和bpi分别为VPi的取值界限。
(4)构建评价对象的现状物元,表达式如下:
式中:Hk为评价对象;Vki为评价对象的第i个评价指标Ci的标准化取值。
(5)建立物元等级的关联度函数及关联度。关联函数表示的是指标的数值与某一评价等级取值范围的程度,经典域物元可拓模型的关联函数公式为:
式中:|Vij|=|bij-aij|;ρ(Voi,Vji)表示点Voi距离区间Vji的距离;ρ(Voi,Vpi)表示点Voi距离区间VPi的距离。
(6)计算待评物元的综合关联度。将各个评价指标与等级之间的关联度进行加权求和,即可求得待评价对象属于各危险度等级的关联度,计算式如下:
式中,Kj(R)为加权关联度值,Kj(Voi)为危险度等级j的关联度,wi为各评价指标的权重。
如果待评价对象的相关物元R属于等级jo,则用下式计算:
令,则可得到等级变量特征值的计算表达式如下:
式中,j*为采空区危险度的等级变量特征值,依据j*可知其偏向相邻等级的程度。
与现有技术相比,本发明的有益效果在于:
(1)本发明基于采空区的实测结果和工程地质调查,综合考虑采空区稳定性的影响因素,选取16个定性和定量指标作为评价指标,建立了采空区危险度等级判别的多指标评价体系,克服了因单一指标不能很好反映采空区危险度而导致评价结果准确性低的缺点;同时应用AHP、熵权法以及博弈论组合赋权法求得各评价指标的综合权重,有效避免了权重求取过程中人为主观因素的影响,确保了本发明方法预测结果的可靠度;
(2)本发明使评价结果直观可靠,运用物元可拓模型可以将定量指标和定性指标学合理地关联到了一起,并将危险度等级这一定性概念转换为最大关联度,实现了定性与定量的转换,结果直观、准确,便于实际应用分析;
(3)本发明具有很好的继承性,只要获得采空区的9个定性指标和7个定量指标值,就可迅速、准确的对采空区安全度等级进行评价,在实际工程中具有一定的应用意义。
附图说明
图1为本发明的流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种地下金属矿采空区危险度分级评价方法,其特征在于,包括以下步骤:
步骤1:对地下金属矿采空区进行现场探测,并收集采空区相关的工程地质资料;
步骤2:选取n个评价指标,建立采空区危险度分级评价指标体系,并确定危险度等级的分级标准;
步骤3:运用AHP、熵权法和博弈论组合赋权法确定出各指标的综合权重;
步骤4:建立采空区危险度评价的物元可拓模型;
步骤5:将各指标权重值代入物元可拓模型,对各项指标进行关联度运算,得到关于分级等级的关联度计算结果;
步骤6:根据最大关联度原则,判定采空区的危险度等级。
具体的,在步骤3中,首先运用层次分析法(AHP)计算出各评价指标的主观权重Wa,并进行总体一致性检验;采用熵权法计算采空区危险度各评价指标的信息熵Ei,得到指标的客观权重Wb;根据AHP、熵权法得到的权重,采用博弈论组合赋权法,计算出各评价指标的综合权重Wc。
所述步骤3中,根据AHP法计算指标的主观权重,并进行总体一致性检验的具体过程包括以下步骤:
(1)以采空区危险度为目标层,以各指标为指标层,建立指标体系的层次结构模型;
(2)采用基于Satty提出的1~9标度法,对各指标对采空区危险度的相对重要性进行比较,用数值表示出来,并写成矩阵形式,即构成判断矩阵X:
X=(xij)m×n (1)
式(1)中,xij表示指标xi与指标xj重要性比较的取值,i,j=1,2,...,n;
(3)计算判断矩阵特征向量W和判断矩阵的最大特征根λmax:
W=(w1,w2,ΛΛ,wn)T (2)
(4)对判断矩阵X进行一致性检验:计算判断矩阵的一致性指标和一致性比率:
式中,C.I.为一致性指标,C.R.为一致性比率,R.I为随机一致性指标,n为判断矩阵的阶数。
当时,即认为判断矩阵具有满意的一致性;
(5)层次总排序及一致性检验,参照步骤(2)-(4),计算出指标层对于采空区危险度相对重要性的排序,即权重Wa,并进行总排序的一致性检验。
所述步骤3中,采用熵权法计算采空区危险度各评价指标的客观权重的过程包括以下步骤:
(1)根据初始数据获得共m个评价(采空区),n个评价指标,构建初始数据矩阵R:
R=(rij)m×n,i=1,2,...,m;j=1,2,...,n; (6)
式(6)中,rij为第j个指标下第i个采空区的评价值;
(2)对指标进行数据无量纲和归一化处理:
对于数值越高越优的正向指标,归一化公式为:
对于数值越低越优的负向指标,归一化公式为:
(3)根据归一化处理后得到的数据r′ij,计算出第j指标下第i个采空区占该指标的比重:
(4)计算第j项指标的信息熵Ei,计算公式为:
(5)根据公式(11)计算出各指标的权重值,得到指标熵权值的排序向量Wb
所述步骤3中,采用博弈论组合赋权法计算采空区危险度各评价指标的综合权重的过程包括以下步骤:
(1)根据L种方法对同一指标获得的L个指标权重值,构造权重向量uk=[uk1,uk2,…ukn],k=1,2,...,L,将L个权重向量任意线性组合记为:
式中:αk为线性组合系数,αk>0;u为可能的权重向量集合;
(2)为了使可能向量集u和各个uk的离差值达到最小,对组合系数αk进行优化,进而得到uk的最优解uk *,计算公式如下:
计算可得:则对应的线性方程组为:
(3)对式(14)进行求解,得到组合系数[α1,α2,…,αL],然后对其进行归一化处理,计算公式如下:
(4)求出采空区危险度各评价指标的综合权重,计算公式如下:
在步骤4中,物元可拓模型包括经典域物元、节点物元和现状物元,构建采空区危险度评价的物元可拓模型的过程包括以下步骤:
(1)物元模型
假设待评价的采空区评价等级为H,评价指标为Ci,各项指标的特征向量值为Vi,以有序三元组R={H,Ci,Vi}构建待评价采空区危险度的物元,如下式所示:
(2)构建评价对象的经典域物元模型,表达式如下:
式中:Rj为第j个同征物元;Hj为评价对象的第j个危险度等级;Ci为评价对象的第i个评价指标,i=1,2,…,n;Vji为危险度等级j第i个评价指标Ci标准化后的经典域;aji和bji分别为Vji的取值界限;
(3)构建评价对象的节域物元模型,计算表达式如下:
式中:HP为评价对象危险度等级的全体;VPi为P关于评价指标Ci所取的值域,即HP节域;api和bpi分别为VPi的取值界限;
(4)构建评价对象的现状物元,表达式如下:
式中:Hk为评价对象;Vki为评价对象的第i个评价指标Ci的标准化取值;
(5)建立物元等级的关联度函数及关联度,关联函数表示的是指标的数值与某一评价等级取值范围的程度,经典域物元可拓模型的关联函数公式为:
式中:|Vij|=|bij-aij|;ρ(Voi,Vji)表示点Voi距离区间Vji的距离;ρ(Voi,Vpi)表示点Voi距离区间VPi的距离;
(6)计算待评物元的综合关联度;将各个评价指标与等级之间的关联度进行加权求和,即可求得待评价对象属于各危险度等级的关联度,计算式如下:
式中,Kj(R)为加权关联度值,Kj(Voi)为危险度等级j的关联度,wi为各评价指标的权重;
如果待评价对象的相关物元R属于等级jo,则用下式计算:
令,则可得到等级变量特征值的计算表达式如下:
式中,j*为采空区危险度的等级变量特征值,依据j*可知其偏向相邻等级的程度。
以下将结合具体实施例对本发明作详细说明:
实施例1
步骤1:对地下金属矿采空区进行现场探测,并收集采空区相关的工程地质资料。
步骤2:建立采空区危险度分级评价指标体系,并确定危险度等级的分级标准;
综合考虑地下金属矿采空区的实际情况,选定采空区类型(X1)、地质构造(X2)、岩体结构(X3)、围岩支护(X4)、顶板特征(X5)、地表特征(X6)、水文因素(X7)、相邻空区分布(X8)、工程布置(X9)、采空区规模(X10)、埋深/m(X11)、高跨比(X12)、暴露面积/m2(X13)、矿柱安全系数(X14)、岩石抗压强度/Mpa(X15)、岩石质量标准(X16)共16项指标建立采空区危险度分级评价指标体系。其中,采空区类型、地质构造、岩体结构、围岩支护、顶板特征、地表特征、水文因素、相邻空区分布、工程布置9个指标为定性指标;采空区规模、埋深、高跨比、暴露面积、矿柱安全系数、岩石抗压强度和岩石质量标准7个指标为定量指标;同时,将危险度等级从低至高划分为四个等级:差(Ⅰ级)、中(Ⅱ级)、良(Ⅲ级)和优(Ⅳ级)。
在确定评价指标体系的基础上,结合相关研究成果,分别获得了危险度等级和各定性、定量评价指标的关系,如表1,表2所示。
表1定量指标分级标准
表2指标分级标准及量化标准赋值
步骤3:运用AHP、熵权法和博弈论组合赋权法确定出各指标的综合权重;
选取一个采空区为实例进行说明,采用AHP、熵权法和博弈论组合赋权法确定采空区评价指标的综合权重,其步骤如下:
(1)构造判断矩阵:将16个评价指标,构建两两比较判断矩阵,如下形式:
(2)根据本发明内容中的步骤3,利用编程技术计算出16个评价指标总排序Wa,并行一致性检验;
(3)根据本发明内容中的步骤3,计算出16个评价指标的客观权重Wb。
(4)根据本发明内容中的步骤3,采用博弈论组合赋权法,计算出各评价指标的综合权重Wc,见表3所示。
表3 AHP-熵权法确定指标权重
步骤4:根据本发明内容中的步骤4,分别确定出实例采空区危险度评价的经典域物元R0、节域物元Rp、现状物元Rk。
步骤5:根据本发明内容中的步骤5,将指标权重值代入物元可拓模型进行关联度运算,得到16项指标关于I、II、III、IV等级的关联度计算结果,如表4所示。
表4指标特征及级别关联度
步骤6:根据本发明内容中的步骤6,按照最大关联度原则,判定采空区的危险度等级,得出待评价采空区关于不同危险度等级的关联度,见表5。
表5采空区危险度分级评价结果
根据最大关联度原则,该采空区危险度属于Ⅰ(差)级。该评价结果与采空区的实际情况基本一致。
上述实施例仅仅是清楚地说明本发明所作的举例,而非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里也无需也无法对所有的实施例予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明的保护范围之中。
Claims (8)
1.一种地下金属矿采空区危险度分级评价方法,其特征在于,包括以下步骤:
步骤1:对地下金属矿采空区进行现场探测,并收集采空区相关的工程地质资料;
步骤2:选取n个评价指标,建立采空区危险度分级评价指标体系,并确定危险度等级的分级标准;
步骤3:运用AHP、熵权法和博弈论组合赋权法确定出各指标的综合权重;
步骤4:建立采空区危险度评价的物元可拓模型;
步骤5:将各指标权重值代入物元可拓模型,对各项指标进行关联度运算,得到关于分级等级的关联度计算结果;
步骤6:根据最大关联度原则,判定采空区的危险度等级。
2.根据权利1所述的评价方法,其特征在于:步骤2中选取采空区类型、地质构造、岩体结构、围岩支护、顶板特征、地表特征、水文因素、相邻空区分布、工程布置共9个定性指标,以及采空区规模、埋深、高跨比、暴露面积、矿柱安全系数、岩石抗压强度和岩石质量标准共7个定量指标,合计16个指标建立采空区危险度分级评价指标体系,并将危险度等级从低至高分为差(Ⅰ级)、中(Ⅱ级)、良(Ⅲ级)和优(Ⅳ级)四个等级。
3.根据权利2所述的评价方法,其特征在于:在步骤3中,首先运用层次分析法(AHP)计算出各评价指标的主观权重Wa,并进行总体一致性检验;采用熵权法计算采空区危险度各评价指标的信息熵Ei,得到指标的客观权重Wb;根据AHP、熵权法得到的权重,采用博弈论组合赋权法,计算出各评价指标的综合权重Wc。
4.根据权利3所述的评价方法,其特征在于:所述步骤3中,根据AHP法计算指标的主观权重,并进行总体一致性检验的具体过程包括以下步骤:
(1)以采空区危险度为目标层,以各指标为指标层,建立指标体系的层次结构模型;
(2)采用基于Satty提出的1~9标度法,对各指标对采空区危险度的相对重要性进行比较,用数值表示出来,并写成矩阵形式,即构成判断矩阵X:
X=(xij)m×n (1)
式(1)中,xij表示指标xi与指标xj重要性比较的取值,i,j=1,2,...,n;
(3)计算判断矩阵特征向量W和判断矩阵的最大特征根λmax:
W=(w1,w2,ΛΛ,wn)T (2)
(4)对判断矩阵X进行一致性检验:计算判断矩阵的一致性指标和一致性比率:
式中,C.I.为一致性指标,C.R.为一致性比率,R.I为随机一致性指标;
当时,即认为判断矩阵具有满意的一致性;
(5)层次总排序及一致性检验,参照步骤(2)-(4),计算出指标层对于采空区危险度相对重要性的排序,即主观权重Wa,并进行总排序的一致性检验。
5.根据权利3所述的评价方法,其特征在于:所述步骤3中,采用熵权法计算采空区危险度各评价指标的客观权重的过程包括以下步骤:
(1)根据初始数据获得共m个评价采空区,n个评价指标,构建初始数据矩阵R:
R=(rij)m×n,i=1,2,...,m;j=1,2,...,n; (6)
式(6)中,rij为第j个指标下第i个采空区的评价值;
(2)对指标进行数据无量纲和归一化处理:
对于数值越高越优的正向指标,归一化公式为:
对于数值越低越优的负向指标,归一化公式为:
(3)根据归一化处理后得到的数据r′ij,计算出第j指标下第i个采空区占该指标的比重:
(4)计算第j项指标的信息熵Ei,计算公式为:
(5)根据公式(11)计算出各指标的权重值,得到指标熵权值的排序向量,也即客观权重Wb
6.根据权利3所述的评价方法,其特征在于:所述步骤3中,采用博弈论组合赋权法计算采空区危险度各评价指标的综合权重的过程包括以下步骤:
(1)根据L种方法对同一指标获得的L个指标权重值,构造权重向量uk=[uk1,uk2,…ukn],k=1,2,…,L,将L个权重向量任意线性组合记为:
式中:αk为线性组合系数,αk>0;u为可能的权重向量集合;
(2)为了使可能向量集u和各个uk的离差值达到最小,对组合系数αk进行优化,进而得到uk的最优解uk *,计算公式如下:
计算可得:则对应的线性方程组为:
(3)对式(14)进行求解,得到组合系数[α1,α2,…,αL],然后对其进行归一化处理,计算公式如下:
(4)求出采空区危险度各评价指标的综合权重,计算公式如下:
7.根据权利2所述的评价方法,其特征在于:在步骤4中,物元可拓模型包括经典域物元、节点物元和现状物元。
8.根据权利7所述的评价方法,其特征在于:构建采空区危险度评价的物元可拓模型的过程包括以下步骤:
(1)物元模型
假设待评价的采空区评价等级为H,评价指标为Ci,各项指标的特征向量值为Vi,以有序三元组R={H,Ci,Vi}构建待评价采空区危险度的物元,如下式所示:
(2)构建评价对象的经典域物元模型,表达式如下:
式中:Rj为第j个同征物元;Hj为评价对象的第j个危险度等级;Ci为评价对象的第i个评价指标,i=1,2,…,n;Vji为危险度等级j第i个评价指标Ci标准化后的经典域;aji和bji分别为Vji的取值界限;
(3)构建评价对象的节域物元模型,计算表达式如下:
式中:HP为评价对象危险度等级的全体;VPi为P关于评价指标Ci所取的值域,即HP节域;api和bpi分别为VPi的取值界限;
(4)构建评价对象的现状物元,表达式如下:
式中:Hk为评价对象;Vki为评价对象的第i个评价指标Ci的标准化取值;
(5)建立物元等级的关联度函数及关联度,关联函数表示的是指标的数值与某一评价等级取值范围的程度,经典域物元可拓模型的关联函数公式为:
式中:|Vij|=|bij-aij|;ρ(Voi,Vji)表示点Voi距离区间Vji的距离;ρ(Voi,Vpi)表示点Voi距离区间VPi的距离;
(6)计算待评物元的综合关联度;将各个评价指标与等级之间的关联度进行加权求和,即可求得待评价对象属于各危险度等级的关联度,计算式如下:
式中,Kj(R)为加权关联度值,Kj(Voi)为危险度等级j的关联度,wi为各评价指标的权重;
如果待评价对象的相关物元R属于等级jo,则用下式计算:
Kj0(R)=max(Kj(R)) (25)
令,则可得到等级变量特征值的计算表达式如下:
式中,j*为采空区危险度的等级变量特征值,依据j*可知其偏向相邻等级的程度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910659265.5A CN110610285A (zh) | 2019-07-22 | 2019-07-22 | 一种地下金属矿采空区危险度分级评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910659265.5A CN110610285A (zh) | 2019-07-22 | 2019-07-22 | 一种地下金属矿采空区危险度分级评价方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110610285A true CN110610285A (zh) | 2019-12-24 |
Family
ID=68890266
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910659265.5A Pending CN110610285A (zh) | 2019-07-22 | 2019-07-22 | 一种地下金属矿采空区危险度分级评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110610285A (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111126858A (zh) * | 2019-12-26 | 2020-05-08 | 冯威 | 一种大型超市管理绩效评价系统 |
CN111178749A (zh) * | 2019-12-26 | 2020-05-19 | 冯威 | 一种儿童生长发育状况的评价系统 |
CN111523796A (zh) * | 2020-04-21 | 2020-08-11 | 西南石油大学 | 非煤隧道有害气体危害评价方法 |
CN111582706A (zh) * | 2020-04-30 | 2020-08-25 | 城云科技(中国)有限公司 | 一种基于智慧城市综合治理的数字化管理系统 |
CN111967777A (zh) * | 2020-08-19 | 2020-11-20 | 国网河南省电力公司经济技术研究院 | 一种储能促进清洁能源利用价值的综合评价方法 |
CN111967776A (zh) * | 2020-08-19 | 2020-11-20 | 国网河南省电力公司经济技术研究院 | 一种园区综合能源系统运营价值链的评估方法 |
CN112070367A (zh) * | 2020-08-20 | 2020-12-11 | 中国人民解放军海军工程大学 | 一种基于改进综合成新率的大型武器装备价值评估方法 |
CN112580078A (zh) * | 2020-12-22 | 2021-03-30 | 江苏方天电力技术有限公司 | 一种针对电力数据指纹的评估方法 |
CN112633672A (zh) * | 2020-12-18 | 2021-04-09 | 鄂尔多斯市伊化矿业资源有限责任公司 | 一种采空区大面积悬顶危险性评估系统及方法 |
CN112712484A (zh) * | 2020-12-07 | 2021-04-27 | 鞍钢集团矿业有限公司 | 一种采用模糊数学理论圈定采空区电阻率影像异常的方法 |
CN113887978A (zh) * | 2021-10-13 | 2022-01-04 | 北京工商大学 | 一种基于改进物元可拓模型的食品风险综合评价方法 |
CN117194879A (zh) * | 2023-09-08 | 2023-12-08 | 哈尔滨工业大学 | 一种针对寒区冬季河流的水环境健康多指标体系评价算法 |
CN117474329A (zh) * | 2023-10-31 | 2024-01-30 | 青岛理工大学 | 一种基于修正ahp法和可拓理论的基坑风险评估方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103310390A (zh) * | 2013-07-05 | 2013-09-18 | 武汉大学 | 一种电网安全综合评价方法 |
CN103839057A (zh) * | 2014-03-28 | 2014-06-04 | 中南大学 | 一种锑浮选工况识别方法及系统 |
CN106897555A (zh) * | 2017-02-22 | 2017-06-27 | 辽宁工程技术大学 | 一种露天煤矿气体环境的评价方法 |
CN107330590A (zh) * | 2017-06-19 | 2017-11-07 | 华北电力大学 | 一种基于熵值法和物元可拓法的核电站安全运行评价方法 |
CN109118121A (zh) * | 2018-09-13 | 2019-01-01 | 中国电力科学研究院有限公司 | 基于改进PageRank的电力市场出清方式评价方法 |
-
2019
- 2019-07-22 CN CN201910659265.5A patent/CN110610285A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103310390A (zh) * | 2013-07-05 | 2013-09-18 | 武汉大学 | 一种电网安全综合评价方法 |
CN103839057A (zh) * | 2014-03-28 | 2014-06-04 | 中南大学 | 一种锑浮选工况识别方法及系统 |
CN106897555A (zh) * | 2017-02-22 | 2017-06-27 | 辽宁工程技术大学 | 一种露天煤矿气体环境的评价方法 |
CN107330590A (zh) * | 2017-06-19 | 2017-11-07 | 华北电力大学 | 一种基于熵值法和物元可拓法的核电站安全运行评价方法 |
CN109118121A (zh) * | 2018-09-13 | 2019-01-01 | 中国电力科学研究院有限公司 | 基于改进PageRank的电力市场出清方式评价方法 |
Non-Patent Citations (2)
Title |
---|
唐朝晖: "基于泡沫图像特征和物元可拓模型的锑浮选工况识别", 《控制与决策》 * |
李山存: "金属矿山地下采空区定量风险评价方法研究", 《现代矿业》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111126858A (zh) * | 2019-12-26 | 2020-05-08 | 冯威 | 一种大型超市管理绩效评价系统 |
CN111178749A (zh) * | 2019-12-26 | 2020-05-19 | 冯威 | 一种儿童生长发育状况的评价系统 |
CN111523796A (zh) * | 2020-04-21 | 2020-08-11 | 西南石油大学 | 非煤隧道有害气体危害评价方法 |
CN111582706A (zh) * | 2020-04-30 | 2020-08-25 | 城云科技(中国)有限公司 | 一种基于智慧城市综合治理的数字化管理系统 |
CN111967777A (zh) * | 2020-08-19 | 2020-11-20 | 国网河南省电力公司经济技术研究院 | 一种储能促进清洁能源利用价值的综合评价方法 |
CN111967776A (zh) * | 2020-08-19 | 2020-11-20 | 国网河南省电力公司经济技术研究院 | 一种园区综合能源系统运营价值链的评估方法 |
CN111967777B (zh) * | 2020-08-19 | 2023-10-27 | 国网河南省电力公司经济技术研究院 | 一种储能促进清洁能源利用价值的综合评价方法 |
CN112070367A (zh) * | 2020-08-20 | 2020-12-11 | 中国人民解放军海军工程大学 | 一种基于改进综合成新率的大型武器装备价值评估方法 |
CN112712484B (zh) * | 2020-12-07 | 2023-10-17 | 鞍钢集团矿业有限公司 | 一种采用模糊数学理论圈定采空区电阻率影像异常的方法 |
CN112712484A (zh) * | 2020-12-07 | 2021-04-27 | 鞍钢集团矿业有限公司 | 一种采用模糊数学理论圈定采空区电阻率影像异常的方法 |
CN112633672A (zh) * | 2020-12-18 | 2021-04-09 | 鄂尔多斯市伊化矿业资源有限责任公司 | 一种采空区大面积悬顶危险性评估系统及方法 |
CN112580078B (zh) * | 2020-12-22 | 2022-11-08 | 江苏方天电力技术有限公司 | 一种针对电力数据指纹的评估方法 |
CN112580078A (zh) * | 2020-12-22 | 2021-03-30 | 江苏方天电力技术有限公司 | 一种针对电力数据指纹的评估方法 |
CN113887978A (zh) * | 2021-10-13 | 2022-01-04 | 北京工商大学 | 一种基于改进物元可拓模型的食品风险综合评价方法 |
CN113887978B (zh) * | 2021-10-13 | 2024-05-31 | 北京工商大学 | 一种基于改进物元可拓模型的食品风险综合评价方法 |
CN117194879A (zh) * | 2023-09-08 | 2023-12-08 | 哈尔滨工业大学 | 一种针对寒区冬季河流的水环境健康多指标体系评价算法 |
CN117474329A (zh) * | 2023-10-31 | 2024-01-30 | 青岛理工大学 | 一种基于修正ahp法和可拓理论的基坑风险评估方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110610285A (zh) | 一种地下金属矿采空区危险度分级评价方法 | |
CN109034582B (zh) | 基于云模型和组合赋权的隧道穿越断层突水突泥风险评价方法 | |
CN107943880B (zh) | 一种基于层次分析法的地质灾害易发性改进评估方法 | |
KR101642951B1 (ko) | Gis 기반 실시간 지진피해 예측 방법 | |
CN110059963A (zh) | 一种基于模糊多态贝叶斯网络的隧道风险评价方法 | |
CN103268420A (zh) | 一种岩石高边坡的危险性评价方法 | |
CN110210776A (zh) | 一种基于监控量测的隧道施工动态风险评估方法 | |
CN104133985A (zh) | 隧道及地下工程地质灾害风险属性区间评价方法 | |
CN108332696B (zh) | 滑坡监测方法选择方法 | |
CN111022119B (zh) | 一种基于模式识别的盾构隧道结构抗震安全评价方法 | |
CN114926067A (zh) | 一种基于信息量模型的层次分析地质灾害风险评价方法 | |
Wang et al. | Risk Assessment of Water Inrush in Karst Tunnels Based on the Efficacy Coefficient Method. | |
CN111242499A (zh) | 基于层次-可拓分析的既有隧道衬砌结构病害评价方法 | |
Wang et al. | A novel model of the ideal point method coupled with objective and subjective weighting method for evaluation of surrounding rock stability | |
CN114881396A (zh) | 一种基于ahp和topsis的隧道坍塌风险评估方法 | |
CN112529267A (zh) | 一种基于模糊综合评价法的超前地质预报方法及系统 | |
CN116434475A (zh) | 一种基于降雨强度与地表位移变化率的滑坡预警方法 | |
CN112907130A (zh) | 基于盾构施工参数的施工风险等级确定方法、装置及介质 | |
Hu et al. | Hazard degree identification of goafs based on scale effect of structure by RS-TOPSIS method | |
Wang et al. | Risk Assessment of Water Inrush in Karst Tunnels Based on the Ideal Point Method. | |
CN115271297A (zh) | 一种基于ahp和综合指数法的屋后切坡危险性评价方法 | |
CN105046095B (zh) | 生态敏感区公路隧道施工对地下水环境扰动的评价方法 | |
CN110499768A (zh) | 一种基于层次分析法的边坡安全稳定性等级评价方法 | |
CN117172547A (zh) | 一种基于未确知测度理论的公路隧道施工安全评价方法 | |
CN114493295B (zh) | 一种盐岩地质下隧道监测项目的确定方法及监测系统 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191224 |