CN112699553A - 一种岩爆烈度等级智能预测系统方法 - Google Patents

一种岩爆烈度等级智能预测系统方法 Download PDF

Info

Publication number
CN112699553A
CN112699553A CN202011595138.2A CN202011595138A CN112699553A CN 112699553 A CN112699553 A CN 112699553A CN 202011595138 A CN202011595138 A CN 202011595138A CN 112699553 A CN112699553 A CN 112699553A
Authority
CN
China
Prior art keywords
rock burst
index
weight
prediction
entropy
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.)
Granted
Application number
CN202011595138.2A
Other languages
English (en)
Other versions
CN112699553B (zh
Inventor
李克钢
李明亮
秦庆词
赵冰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN202011595138.2A priority Critical patent/CN112699553B/zh
Publication of CN112699553A publication Critical patent/CN112699553A/zh
Application granted granted Critical
Publication of CN112699553B publication Critical patent/CN112699553B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Development Economics (AREA)
  • Geometry (AREA)
  • Game Theory and Decision Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开一种岩爆烈度等级智能预测系统方法，属于地下工程及岩土工程技术领域。本发明所述方法包括选取岩爆烈度等级预测指标；构建岩爆烈度等级预测指标体系；确定岩爆烈度等级的具体分级标准；搜集国内外典型的岩爆工程案例数据；通过德尔菲法和改进熵权法，确定岩爆烈度等级预测指标的主客观权重值；采用偏好系数法，确定岩爆烈度等级预测指标的最优化组合权重值；运用多维正态云模型理论和python编程手段，确定岩爆烈度等级的综合隶属度；根据最大隶属度原则确定岩爆烈度等级；基于python软件自主开发一种岩爆烈度等级预测系统；该岩爆预测方法拥有更好的准确性和适用性，并且，可以大量节省运算时间以及计算步骤简洁明了。

Description

一种岩爆烈度等级智能预测系统方法
技术领域
本发明涉及涉及一种岩爆烈度等级智能预测系统方法，属于深部地下工程及岩土工程技术领域。
背景技术
岩爆作为矿山、隧道、水电站等深部岩土工程施工中常遇到的一种地质灾害，常常表现为岩石极速飞弾、围岩突然抛掷、爆裂和剥离等岩石动力学灾害现象。随着岩土工程领域不断地向深部进军，遇到的岩爆现象日益增多，因此，岩爆灾害的发生严重威胁地下施工作业人员以及地下设备的安全，直接影响企业施工进度，如遇到强烈岩爆等级，则会毁坏大面积工程。因此，岩爆烈度等级预测对岩土工程领域施工程度以及作业人员至关重要。
随着国力的不断提升，众多学者对岩爆烈度等级预测研究做了大量工作，主要从不同的非线性理论、不同的数值分析、数值仿真技术和一些监测设备着手，同时取得了一些可观的效果，但是目前基于对岩爆机理尚不清楚，考虑到现场施工监测技术成本过高和现有的研究方法对岩爆烈度等级预测准确率较低等问题。
发明内容
针对现有岩爆烈度等级预测方法的不足，本发明提出一种岩爆烈度等级智能预测系统方法，旨在准确、客观和简便地预测地下工程及岩土工程中的岩爆烈度等级。
本发明采用的技术方案是：一种岩爆烈度等级智能预测系统方法，包括以下步骤：
步骤一：分析影响岩爆烈度等级预测指标；
步骤二：构建岩爆烈度等级预测指标体系；
步骤三：确定岩爆烈度等级的具体分级标准；
步骤四：搜集典型的岩爆工程数据；
步骤五：基于改进熵权法结合步骤四中的岩爆工程数据，确定岩爆烈度等级预测指标的客观权重值；
步骤六：基于德尔菲法，确定岩爆烈度等级预测指标的主观权重值；
步骤七：基于偏好系数法，确定岩爆烈度等级预测指标的最优化组合权重值；
步骤八：基于多维云模型原理结合步骤二和步骤三，确定岩爆烈度等级预测指标对应的多维云模型的数字特征：期望、熵、超熵；
步骤九：基于多维云模型的云正向发生器，结合步骤八，确定岩爆烈度等级预测多维云模型；
步骤十：基于多维云模型的云正向发生器结合步骤四和步骤七，确定步骤四中的岩爆烈度等级的综合隶属度；
步骤十一：基于最大隶属度原则确定步骤四中的岩爆烈度等级；
步骤十二：基于python软件开发一种岩爆烈度等级预测系统。
具体地，所述步骤一中，所选取的岩爆烈度等级预测指标共有三个：依次为弹性变形能指数Wet、岩石脆性系数σc/σt和岩体应力系数σθ/σc；步骤二中，选取弹性变形能指数Wet、岩体应力系数σθ/σc和岩石脆性系数σc/σt，建立岩爆烈度等级预测指标体系；步骤三中，岩爆烈度等级依次为：无岩爆(I)、轻微岩爆(II)、中等岩爆(III)和强烈岩爆(IV)四个等级。
具体地，所述步骤五中改进熵权法过程为：
(1)极值法原始数据矩阵Am×n标准化处理：Am×n＝(aij)m×n(i＝1,2,...,n)
式中：m为具体的岩爆案例工程数据个数；n为岩爆烈度等级预测指标个数；aij为原始矩阵的元素；
(2)计算熵hi：
Figure BDA0002867990800000021
其中，
Figure BDA0002867990800000022
计算熵权wi：
Figure BDA0002867990800000023
式中：hi为熵；fij为特征比重；rij为归一化后矩阵的元素；Wi为未改进的熵权法的权重；k为系数，
Figure BDA0002867990800000024
(3)基于秩比法的改进熵权法，step1：分别将各指标对岩爆影响程度由大到小排序赋值zi，最大指标赋值为L，其次分别赋值为L-1,...,影响最小指标赋值为1，则赋值后指标i的偏好校正系数Zi：
Figure BDA0002867990800000025
ri为各指标对岩爆影响程度由大到小排序赋值；step2：基于最小信息熵原理，把偏好校正系数与熵权法确定的指标权重作耦合处理，即得到改进的熵权重值，公式如下：
Figure BDA0002867990800000031
step3：用拉格朗日乘子法计算下面优化问题得到最终的计算公式：
Figure BDA0002867990800000032
式中：WZi为改进后熵权法的权重值。
具体地，所述步骤六中德尔菲法的过程为：
(1)选择咨询专家成员，设计评价意见咨询表格，进行专家咨询，并且回收咨询数据；
(2)整理回收数据，进行权重的计算，其公式如下：
Figure BDA0002867990800000033
VI＝Dj/Ej
式中：E为指标评价重要性均值；Ej为第j个指标指标评价重要性均值；VI为指标重要性评价变异系数；Dj为指标重要性评价的标准差；s为专家咨询人数；aj表示第j位专家的评分值，取值范围为1-10分，1表示极为不重要，10分表示极为重要，其他重要性评分处于极为不重要到极为重要之间，方差可以表示对同一指标不同专家评价分数的离散程度；
对各个评价指标的权重进行归一化处理，公式如下：
Figure BDA0002867990800000034
式中：Wj为德尔菲法权重值。
具体地，所述步骤七采用德尔菲法确定的主观权重与改进熵权法确定的客观权重进行加权综合，为消除主客观因素对于权重精度的影响，引入主、客观权重偏好系数，由于岩爆烈度等级预测指标之间具有一定的相关性问题，引入改进的最优化组合赋权法，计算各指标的最优组合权重，使之得到的改进最优化组合权重系数更准确；
计算偏好系数：
Figure BDA0002867990800000041
式中：α、β为客、主观权重偏好系数；α1为客观权重待定系数，β1为主观权重待定系数；WZi为改进熵权法的权重值；Wj为德尔菲法权重值；Aij为对j个指标下的第i个评价对象的评价值；
计算各指标的最优组合权重：
W优＝αWZi+βWj
式中：W优为改进最优化组合权重。
具体地，所述步骤八和步骤九中多维云模型的构建是确定预测指标岩爆烈度等级分级标准，确定预测指标各等级云特征：期望：Ex1，Ex2，Ex3，超熵：He1，He2，He3，熵：En1，En2，En3，然后基于多维正向云发生器生成多维云模型；多维云模型的期望计算公式为：Ex＝(Bmin+Bmax)/2，多维云模型的超熵的确定可以根据各个预测指标的最大范围，取一个常数h，使He≤0.5；多维云模型的熵计算公式为：En＝(Exmax)/3，式中：Bmin，Bmax分别对应等级标准的上、下限边界值，对于单边界值的情况，如(-∞，Bmax]或者(Bmin，+∞]，根据分析数据实测值的上、下限确定缺省边界；Exmax为某个因子对应的最高级别的期望值。
具体地，所述步骤十中综合隶属度的计算方法如下：
Figure BDA0002867990800000042
式中：u为综合隶属度；x(x1,x2,…,xp)为各评价指标的一次随机取值；Exj为指标的期望；En′j为超熵；j为指标个数；xp为指标取值；En'2为方差的正态分布x～P(Ex，En'2)；其中，En'又是服从以En为期望、He2为方差的正态分布En'～P(En，He2)。注：设U{x1,x2,…,xp}是一个用精确数值表示的p维定量论域，C是U{x1,x2,…,xp}上的定性概念，若定量值x∈U,且X(x1,x2,…,xp)是定性概念C的一次随机实现，若X(x1,x2,…,xp)满足：X(x1,x2,…,xp)～N(Ex(Ex1,Ex2,…,Exp))，(En′(En2′，…,Enp′)2),其中En′(En2′，…,Enp′)～N(En(En1,En2,…,Enp)，(He(He1,He2,…,Hep))2)。
具体地，所述步骤十二中岩爆烈度等级智能预测系统的创建方法如下：基于python编程处理结合提出的一种岩爆烈度等级预测方法，建立一种岩爆烈度等级预测系统，该系统由用户界面、操作界面和数据处理界面组成，用户界面由岩爆图片、用户名称和密码组成；操作界面由3个岩爆烈度等级预测指标、待预测的具体典型岩爆案例工程名称和收集的具体典型岩爆案例工程3个预测指标数据组成；数据处理界面由实际岩爆案例工程岩爆烈度等级、典型岩爆案例工程名称和开始运行按钮以及保存数据案例组成。
本发明的有益效果是：
1、本发明提出一种改进最优化组合赋权-多维正态云模型岩爆列度等级预测方法，综合考虑影响岩爆预测内外因素，选取岩性条件与应力条件下的弹性变形能指数Wet、岩石脆性系数σc/σt和岩体应力系数σθ/σc等3个指标构成岩爆烈度等级分级预测指标体系，本发明引入主客观偏好系数的方法，弥补主客观因素给权重计算精度带来的不足，提高权重计算的精度；
2、本发明采用改进熵权法，旨在降低收集的典型岩爆案例工程中极端数据的影响，采用改进熵权法获得的客观权重值准确度较高；
3、本发明以德尔菲法确定主观权重和改进熵权法确定客观权重，运用改进最优化组合赋权法获得组合权重，弥补单一赋权的不足，充分考虑主客观因素影响和解决样本数据较少导致权重值不准确等问题，结合多维正态云模型原理建立岩爆预测模型，基于python工具直接得到综合隶属度，解决岩爆预测中评价指标随机性、模糊性和计算综合隶属度繁琐等问题，较客观反映岩爆正态分布的特征。
附图说明
图1为岩爆指标评价体系图；
图2为第Ⅰ等级的多维正态云模型图；
图3为第Ⅱ等级的多维正态云模型图；
图4为第Ⅲ等级的多维正态云模型图；
图5为第Ⅳ等级的多维正态云模型图；
图6为改进最优化组合赋权-多维正态云模型的流程图；
图7为岩爆烈度等级预测系统界面。
具体实施方式
旨在较清楚表述使本发明实施例的计算原理、工作过程和方法优点，下面结合附图对本发明技术方案进行详细说明。
实施例1：如图1-7所示，本发明的岩爆烈度等级预测方法按如下步骤进行：
步骤一：分析影响岩爆烈度等级预测指标。
指标的选取会对岩爆预测准确度带来一定的影响，目前岩爆预测指标主要为岩石单轴抗压强度σc、岩石单轴抗拉强度σt、岩体完整性系数Kv、弹性能量指数Wet、硐室最大切向应力σθ和岩石脆性系数(σc/σt)等一些列指标，本发明综合考虑搜集的典型岩爆案例工程数据都具有相同指标特点，并且对选取的指标要服从易测取、可量化计算、较常见等原则，由于岩爆多数发生在水电站、隧道和地下隧硐室等工程中，因此，其岩爆多数发生在应力较集中区域的岩体上，同时，在众多岩爆案例工程中，其断面形式以张拉破坏为主，岩爆亦主要发生在硬岩上，并且目前对岩爆机理的研究多数基于能量理论，而能够较合理反映岩体储存和释放能量的能力为岩爆倾向性指数Wet，基于此，综合上述分析，本发明选取岩体应力系数(σθ/σc)、岩石脆性系数(σc/σt)和弹性变形能系数(Wet)作为岩爆预测指标。
步骤二：构建岩爆烈度等级预测指标体系。
岩爆的发生受众多因素影响，诱发岩爆因素大致分为以下几种：一是岩石本身具有岩爆倾向性和较大的弹脆性等特点；二是岩体处于“三高一扰动”等环境的外部条件；三是岩爆的发生与自身岩性有关，这些诱发岩爆的因素大多具有复杂性、突发性和不确定性。岩爆预测与多种指标体系和岩爆众多判据有关，将岩爆烈度等级分为无岩爆(Ⅰ)、轻微岩爆(Ⅱ)、中等岩爆(Ⅲ)和强烈岩爆(Ⅳ)等4种级别。
步骤三：确定岩爆烈度等级的具体分级标准。
岩爆预测受多种因素影响，主要表现为内因和外因两方面，内外因素分别为岩性条件和应力条件，考虑诱发岩爆产生内外因素，选取弹性变形能指数Wet、岩体应力系数σθ/σc和岩石脆性系数σc/σt，建立岩爆指标评价体系，见图1，参考史秀志在《未确知测度模型在岩爆烈度分级预测中的应用》和李绍红在《基于权重融合和云模型的岩爆倾向性预测研究》等研究成果，整理岩爆烈度等级预测指标的分级标准，如表1所示。
表1岩爆烈度等级预测指标及分级标准
Figure BDA0002867990800000071
步骤四：搜集国内外典型的岩爆工程案例数据。
为验证本发明提出的岩爆烈度等级预测方法的准确性，参考Sajjad Afraeia在《Developing intelligent classification models for rock burst prediction afterrecognizing significant predictor variables,Section 1:Literature review anddata preprocessing procedure》和周健在《Long-term prediction model of rockburstin underground openings usingheuristic algorithms and support vectormachines》中的研究成果，选取国内外若干典型岩爆工程案例实测数据，如表2所示。
表2国内外典型岩爆案例实测数据
Figure BDA0002867990800000072
Figure BDA0002867990800000081
步骤五：基于改进熵权法结合步骤四中的岩爆案例工程数据，确定岩爆烈度等级预测指标的客观权重值，具体步骤如下：
(1)极值法原始数据矩阵Am×n标准化处理：
Am×n＝(aij)m×n(i＝1,2,…,n)
式中：m为具体的岩爆案例工程数据个数；n为岩爆烈度等级预测指标个数；aij为原始矩阵的元素。
归一化后得到了标准化的矩阵Rm×n＝(rij)m×n(i＝1,2,…,n)
标准化公式为
Figure BDA0002867990800000091
式中：rij为归一化后矩阵的元素；aij为原始矩阵的元素。
为预防在后续的计算中出现为0的情况，对评价指标标准化后的值均加上0.1，对熵权法进行改进，使其应用更加广泛。
(2)计算熵hi：
Figure BDA0002867990800000092
其中，
Figure BDA0002867990800000093
式中：hi为熵；fij为归一化后矩阵的元素；rij为归一化后矩阵的元素；n为具体的岩爆案例工程数据个数；k为系数，
Figure BDA0002867990800000094
计算熵权wi：
Figure BDA0002867990800000095
式中：Wi为熵权法(未改进)的权重；hi为熵；m为岩爆烈度等级预测指标个数。
(3)改进熵权法
基于秩比法的改进熵权法，本发明在综合考虑指标对岩爆预测影响的基础上，分别将各指标对岩爆影响程度由大到小排序赋值ri，对岩爆影响越大，则该指标的性能越好。最大指标赋值为L，其次分别赋值为L-1,...,影响最小指标赋值为1，则赋值后指标i的偏好校正系数Zi：
Figure BDA0002867990800000096
式中：Zi为偏好校正系数；ri为各指标对岩爆影响程度由大到小排序赋值；n为最大指标赋值；ri为各指标对岩爆影响程度由大到小排序赋值
基于最小信息熵原理，把偏好校正系数与熵权法确定的指标权重作耦合处理，即得到改进的熵权重值。计算公式如下：
Figure BDA0002867990800000101
式中：Wi为熵权法(未改进)的权重；m为岩爆烈度等级预测指标个数；
针对上述优化问题，基于用拉格朗日乘子法，公式如下：
Figure BDA0002867990800000102
式中：WZi为改进后熵权法的权重；Wi为熵权法(未改进)的权重；
岩爆烈度等级各预测指标客观权重如表3所示。
表3改进熵权法权重计算结果
评价因子 σ<sub>c</sub>/σ<sub>t</sub> σ<sub>θ</sub>/σ<sub>c</sub> Wet
W<sub>i</sub> 0.1889 0.3120 0.4992
WZ<sub>i</sub> 0.1810 0.4029 0.4161
步骤六：基于德尔菲法，确定岩爆烈度等级预测指标的主观权重值，具体步骤如下：
(1)选择咨询专家成员(考虑权威性与专业性)，设计评价意见咨询表格，进行专家咨询，并且回收咨询数据，如表4所示。
表4
指标 专家1 专家2 专家3 专家4
σ<sub>θ</sub>/σ<sub>c</sub> 9 9 9 8
σ<sub>c</sub>/σ<sub>t</sub> 6 7 7 5
Wet 8 7 7 8
(2)整理回收数据，进行权重的计算，其公式如下：
Figure BDA0002867990800000103
VI＝Dj/Ej
式中：E为指标评价重要性均值；Ej为第j个指标指标评价重要性均值；VI为指标重要性评价变异系数；Dj为指标重要性评价的标准差；s为专家咨询人数；aj表示第j位专家的评分值，取值范围为1-10分，1表示极为不重要，10分表示极为要，其他重要性评分处于极为不重要到极为重要之间，方差可以表示对同一指标不同专家评价分数的离散程度。
(3)对各个评价指标的权重进行归一化处理，公式如下：
Figure BDA0002867990800000111
式中：Wj为德尔菲法权重值；。
岩爆烈度等级各预测指标主观权重如表5所示。
表5德尔菲法各指标权重计算结果
Figure BDA0002867990800000112
步骤七：基于偏好系数法，确定岩爆烈度等级预测指标的最优化组合权重值。
综合考虑评价指标的重要性和准确性，本发明采用德尔菲法确定的主观权重与改进熵权法确定的客观权重进行加权综合，为消除主客观因素对于权重精度的影响，引入主、客观权重偏好系数，考虑到岩爆烈度等级预测指标之间具有一定的相关性问题，引入改进的最优化组合赋权法，计算各指标的最优组合权重。
计算偏好系数：
Figure BDA0002867990800000113
式中：α、β为客、主观权重偏好系数；α1为客观权重待定系数，β1为主观权重待定系数；WZi为改进熵权法的权重值；Wj为德尔菲法权重值；Aij为对j个指标下的第i个评价对象的评价值。
计算各指标的最优组合权重：
W优＝αWZi+βWj
式中：W优为改进最优化组合权重。
客、主观权重偏好系数如表6所示。
表6主客观权重偏好系数计算结果
表6主客观权重偏好系数计算结果
α<sub>1</sub>＝0.6863 β<sub>1</sub>＝0.7273
α＝0.4855 β＝0.5145
岩爆烈度等级各预测指标改进最优化组合权重如表7所示。
表7改进最优组合权重计算结果
方法 σ<sub>c</sub>/σ<sub>t</sub> σ<sub>θ</sub>/σ<sub>c</sub> Wet
改进最优化组合权重 0.2308 0.3957 0.3735
步骤八：基于多维云模型原理结合步骤二和步骤三，确定岩爆烈度等级预测指标对应的多维云模型的数字特征(期望、熵、超熵)。注：设U{x1,x2,…,xp}是一个用精确数值表示的p维定量论域，C是U{x1,x2,…,xp}上的定性概念，若定量值x∈U,且X(x1,x2,…,xp)是定性概念C的一次随机实现，若X(x1,x2,…,xp)满足：X(x1,x2,…,xp)～N(Ex(Ex1,Ex2,…,Exp))，(En′(En2′，…,Enp′)2),其中En′(En2′，…,Enp′)～N(En(En1,En2,…,Enp)，(He(He1,He2,…,Hep))2)。
云模型用期望Ex、熵En和超熵He三个数字特征来整体表征定性概念，期望Ex是云滴在论域空间分布中的数学期望；熵En是定性概念的不确定性度量，一方面反映云滴的离散程度，另一方面决定论域空间中可被概念接受的云滴的确定度；超熵He是熵的不确定度量。
多维云模型的期望计算公式：
Ex＝(Bmin+Bmax)/2
式中：Bmin，Bmax分别对应等级标准的上、下限边界值，对于单边界值的情况，如(-∞，Bmax]或者(Bmin，+∞]，可以根据分析数据实测值的上、下限确定缺省边界。
多维云模型的熵计算公式：
En＝(Exmax)/3
式中：Exmax为某个岩爆烈度等级预测指标对应的最高级别的期望值；En为多维云模型的熵。
多维云模型的超熵可以根据各个岩爆烈度等级预测指标的最大范围，取一个常数h，使He≤0.5，若He＞0.5，则云滴与云滴的距离过大，云滴太过于分散，本发明超熵为0.05、0.01。
根据步骤三的岩爆预测指标及分级标准，利用实施例中步骤八里面的公式分别确定各指标各等级的期望Ex，熵En，并根据期望选取合适的超熵He，结果如表8所示。
表8各指标岩爆多维正态云模型的数字特征计算结果
Figure BDA0002867990800000131
步骤九：基于多维云模型的云正向发生器结合步骤八，确定岩爆烈度等级预测多维云模型。
(1)生成以Ex(Ex1，Ex2，…，Exp)为期望值，En(En1，En2，…，Enp)为标准差的正态随机数X(x1，x2，…，xp)；
(2)生成以En(En1，En2，…，Enp)为期望值，He(He1，He2，…，Hep)为标准差的正态随机数En′(En′1，En′2，…，En′p)；
(3)计算隶属度,那么X(x1，x2，…，xp)，u(x(x1，x2，...，xp))为多维正态云的一个云滴。
Figure BDA0002867990800000132
(4)重复步骤(1)～(3)，直到产生N个云滴，进而生成云模型。
以岩体应力系数σθ/σc和岩石脆性系数σc/σt作为岩爆烈度等级评价指标，表8各指标岩爆多维正态云模型的数字特征计算结果，采用多维云模型云正向发生器分别建立4个岩爆烈度等级依次为Ⅰ、Ⅱ、Ⅲ、Ⅳ等级的二维正态云模型，如图2～5所示。
步骤十：基于多维云模型的云正向发生器结合步骤四和步骤七，确定步骤四中的岩爆烈度等级的综合隶属度。
综合隶属度的计算方法如下：
将实施例中表7各指标岩爆多维正态云模型的数字特征与实施例中表7中各指标最优化组合权重带入实施例中隶属度公式，建立三个岩爆级别的多维云模型生成器。
步骤十一：基于最大隶属度原则确定步骤四中的岩爆烈度等级。
基于实施例中表2岩爆案例实测数据，带入建立的多维云模型生成器中，得到各岩爆等级的综合隶属度，计算出各岩爆等级的综合隶属度后，根据最大隶属度原则即可预测岩爆等级，本发明列出收集的岩爆案例工程预测结果，具体预测结果如表9所示。
表9为本发明通过改进最优化组合赋权-多维云模型计算出去收集的岩爆案例工程的综合隶属度
表9
Figure BDA0002867990800000141
Figure BDA0002867990800000151
Figure BDA0002867990800000161
注：*表示预测结果与实际不符(一共错4组，分别为第2、12、63和82组)；该模型的正确率为95.6％。
将本发明所提出岩爆烈度等级预测方法应用于锦屏二级水电站引水隧洞中，预测结果与实际等级吻合度较好，如表10所示。
表10锦屏二级水电站引水隧洞岩爆预测结果
Figure BDA0002867990800000162
Figure BDA0002867990800000171
注：*表示预测结果与实际不符
步骤十二：基于python软件开发一种岩爆烈度等级预测系统。
本发明基于提出的一种岩爆烈度等级预测方法，运用python编程处理，建立一种岩爆烈度等级预测系统，如图7所示，该系统可以同时预测多种岩爆工程案例发生的岩爆烈度等级。据图7的岩爆烈度等级预测系统界面，只需要输入岩爆预测指标σθ/σc、σc/σt、Wet对应的数值，就可以计算出岩爆的烈度等级，方便于矿山、隧道、水电站等岩土工程应用。
该岩爆烈度等级预测方法拥有更好的准确性和适用性，可以为矿山、隧道、水电站等岩土工程岩爆烈度等级预测问题提供一种新途径。
上面结合附图对本发明的具体实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。

Claims (8)

1.一种岩爆烈度等级智能预测系统方法，其特征在于：包括以下步骤：
步骤一：分析影响岩爆烈度等级预测指标；
步骤二：构建岩爆烈度等级预测指标体系；
步骤三：确定岩爆烈度等级的具体分级标准；
步骤四：搜集典型的岩爆工程数据；
步骤五：基于改进熵权法结合步骤四中的岩爆工程数据，确定岩爆烈度等级预测指标的客观权重值；
步骤六：基于德尔菲法，确定岩爆烈度等级预测指标的主观权重值；
步骤七：基于偏好系数法，确定岩爆烈度等级预测指标的最优化组合权重值；
步骤八：基于多维云模型原理结合步骤二和步骤三，确定岩爆烈度等级预测指标对应的多维云模型的数字特征：期望、熵、超熵；
步骤九：基于多维云模型的云正向发生器，结合步骤八，确定岩爆烈度等级预测多维云模型；
步骤十：基于多维云模型的云正向发生器结合步骤四和步骤七，确定步骤四中的岩爆烈度等级的综合隶属度；
步骤十一：基于最大隶属度原则确定步骤四中的岩爆烈度等级；
步骤十二：基于python软件开发一种岩爆烈度等级预测系统。
2.根据权利要求1所述的一种岩爆烈度等级智能预测系统方法，其特征在于：所述步骤一中，所选取的岩爆烈度等级预测指标共有三个：依次为弹性变形能指数Wet、岩石脆性系数σc/σt和岩体应力系数σθ/σc；步骤二中，选取弹性变形能指数Wet、岩体应力系数σθ/σc和岩石脆性系数σc/σt，建立岩爆烈度等级预测指标体系；步骤三中，岩爆烈度等级依次为：无岩爆(I)、轻微岩爆(II)、中等岩爆(III)和强烈岩爆(IV)四个等级。
3.根据权利要求2所述的一种岩爆烈度等级智能预测系统方法，其特征在于：所述步骤五中改进熵权法过程为：
(1)极值法原始数据矩阵Am×n标准化处理：Am×n＝(aij)m×n(i＝1,2,…,n)
式中：m为具体的岩爆案例工程数据个数；n为岩爆烈度等级预测指标个数；aij为原始矩阵的元素；
(2)计算熵hi：
Figure FDA0002867990790000021
其中，
Figure FDA0002867990790000022
计算熵权wi：
Figure FDA0002867990790000023
式中：hi为熵；fij为特征比重；rij为归一化后矩阵的元素；Wi为未改进的熵权法的权重；k为系数，
Figure FDA0002867990790000024
(3)基于秩比法的改进熵权法，step1：分别将各指标对岩爆影响程度由大到小排序赋值zi，最大指标赋值为L，其次分别赋值为L-1,...,影响最小指标赋值为1，则赋值后指标i的偏好校正系数Zi：
Figure FDA0002867990790000025
ri为各指标对岩爆影响程度由大到小排序赋值；step2：基于最小信息熵原理，把偏好校正系数与熵权法确定的指标权重作耦合处理，即得到改进的熵权重值，公式如下：
Figure FDA0002867990790000026
st:
Figure FDA0002867990790000027
step3：用拉格朗日乘子法计算下面优化问题得到最终的计算公式：
Figure FDA0002867990790000028
式中：WZi为改进后熵权法的权重值。
4.根据权利要求2所述的一种岩爆烈度等级智能预测系统方法，其特征在于：所述步骤六中德尔菲法的过程为：
(1)选择咨询专家成员，设计评价意见咨询表格，进行专家咨询，并且回收咨询数据；
(2)整理回收数据，进行权重的计算，其公式如下：
Figure FDA0002867990790000029
VI＝Dj/Ej
式中：E为指标评价重要性均值；Ej为第j个指标指标评价重要性均值；VI为指标重要性评价变异系数；Dj为指标重要性评价的标准差；s为专家咨询人数；aj表示第j位专家的评分值，取值范围为1-10分，1表示极为不重要，10分表示极为重要，其他重要性评分处于极为不重要到极为重要之间，方差可以表示对同一指标不同专家评价分数的离散程度；
对各个评价指标的权重进行归一化处理，公式如下：
Figure FDA0002867990790000031
式中：Wj为德尔菲法权重值。
5.根据权利要求3所述的一种岩爆烈度等级智能预测系统方法，其特征在于：所述步骤七采用德尔菲法确定的主观权重与改进熵权法确定的客观权重进行加权综合，为消除主客观因素对于权重精度的影响，引入主、客观权重偏好系数，由于岩爆烈度等级预测指标之间具有一定的相关性问题，引入改进的最优化组合赋权法，计算各指标的最优组合权重，使之得到的改进最优化组合权重系数更准确；
计算偏好系数：
Figure FDA0002867990790000032
式中：α、β为客、主观权重偏好系数；α1为客观权重待定系数，β1为主观权重待定系数；WZi为改进熵权法的权重值；Wj为德尔菲法权重值；Aij为对j个指标下的第i个评价对象的评价值；
计算各指标的最优组合权重：
W优＝αWZi+βWj
式中：W优为改进最优化组合权重。
6.根据权利要求2所述的一种岩爆烈度等级智能系统预测方法，其特征在于：所述步骤八和步骤九中多维云模型的构建是确定预测指标岩爆烈度等级分级标准，确定预测指标各等级云特征：期望：Ex1，Ex2，Ex3，超熵：He1，He2，He3，熵：En1，En2，En3，然后基于多维正向云发生器生成多维云模型；多维云模型的期望计算公式为：Ex＝(Bmin+Bmax)/2，多维云模型的超熵的确定可以根据各个预测指标的最大范围，取一个常数h，使He≤0.5；多维云模型的熵计算公式为：En＝(Exmax)/3，式中：Bmin，Bmax分别对应等级标准的上、下限边界值，对于单边界值的情况，如(-∞，Bmax]或者(Bmin，+∞]，根据分析数据实测值的上、下限确定缺省边界；Exmax为某个因子对应的最高级别的期望值。
7.根据权利要求6所述的一种岩爆烈度等级智能预测系统方法，其特征在于：所述步骤十中综合隶属度的计算方法如下：
Figure FDA0002867990790000041
式中：u为综合隶属度；x(x1，x2，…，xp)为各评价指标的一次随机取值；Exj为指标的期望；En′j为超熵；j为指标个数；xp为指标取值；En'2为方差的正态分布x～P(Ex，En'2)；其中，En'又是服从以En为期望、He2为方差的正态分布En'～P(En，He2)；注：设U{x1,x2,…,xp}是一个用精确数值表示的p维定量论域，C是U{x1,x2,…,xp}上的定性概念，若定量值x∈U,且X(x1,x2,…,xp)是定性概念C的一次随机实现，若X(x1,x2,…,xp)满足：X(x1,x2,…,xp)～N(Ex(Ex1,Ex2,…,Exp))，(En′(En2′，…,Enp′)2),其中En′(En2′,…,Enp′)～N(En(En1,En2,…,Enp)，(He(He1,He2,…,Hep))2)。
8.根据权利要求2所述的一种岩爆烈度等级智能预测系统方法，其特征在于：所述步骤十二中岩爆烈度等级智能预测系统的创建方法如下：基于python编程处理结合提出的一种岩爆烈度等级预测方法，建立一种岩爆烈度等级预测系统，该系统由用户界面、操作界面和数据处理界面组成，用户界面由岩爆图片、用户名称和密码组成；操作界面由3个岩爆烈度等级预测指标、待预测的具体典型岩爆案例工程名称和收集的具体典型岩爆案例工程3个预测指标数据组成；数据处理界面由实际岩爆案例工程岩爆烈度等级、典型岩爆案例工程名称和开始运行按钮以及保存数据案例组成。
CN202011595138.2A 2020-12-29 2020-12-29 一种岩爆烈度等级智能预测系统方法 Active CN112699553B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011595138.2A CN112699553B (zh) 2020-12-29 2020-12-29 一种岩爆烈度等级智能预测系统方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011595138.2A CN112699553B (zh) 2020-12-29 2020-12-29 一种岩爆烈度等级智能预测系统方法

Publications (2)

Publication Number Publication Date
CN112699553A true CN112699553A (zh) 2021-04-23
CN112699553B CN112699553B (zh) 2021-11-02

Family

ID=75512018

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011595138.2A Active CN112699553B (zh) 2020-12-29 2020-12-29 一种岩爆烈度等级智能预测系统方法

Country Status (1)

Country Link
CN (1) CN112699553B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113076700A (zh) * 2021-04-27 2021-07-06 昆明理工大学 一种基于数据分析原理的svm-lda岩爆机器学习预测模型方法
CN113221471A (zh) * 2021-06-15 2021-08-06 贵州大学 基于改进可变模糊集与神经网络的岩爆等级判别方法
CN115526036A (zh) * 2022-09-19 2022-12-27 长安大学 一种岩爆规模等级的判别方法及系统
CN116843228A (zh) * 2023-07-21 2023-10-03 昆明理工大学 基于模糊集理论和组合赋权的煤层冲击倾向性评价方法
CN117110584A (zh) * 2023-10-23 2023-11-24 中铁四局集团有限公司 一种应变型隧道岩爆烈度等级预测方法
CN117272644A (zh) * 2023-09-21 2023-12-22 安徽理工大学 基于钻孔参数的岩爆倾向性预测方法及系统
CN117332240A (zh) * 2023-12-01 2024-01-02 中铁四局集团有限公司 岩爆预测模型构建方法、存储介质、岩爆预测方法及系统

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104008302A (zh) * 2014-06-09 2014-08-27 上海电力学院 一种基于组合赋权和模糊评分的配电网可靠性评估方法
CN104732070A (zh) * 2015-02-27 2015-06-24 广西大学 一种基于信息向量机的岩爆等级预测方法
RU2560008C1 (ru) * 2014-06-16 2015-08-20 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Национальный минерально-сырьевой университет "Горный" Способ крепления горных выработок на удароопасных пластах
CN106407493A (zh) * 2016-03-15 2017-02-15 中南大学 一种基于多维高斯云模型的岩爆等级评价方法
CN107748103A (zh) * 2017-09-01 2018-03-02 中国科学院武汉岩土力学研究所 一种隧道岩爆预测方法、设备、存储介质和系统
CN109523200A (zh) * 2019-01-03 2019-03-26 东北大学 基于最小二乘法的铸造车间安全评价综合权重计算方法
CN109934500A (zh) * 2019-03-15 2019-06-25 南方电网科学研究院有限责任公司 一种配电台区的低电压评价方法、装置、设备及存储介质
JP2019199708A (ja) * 2018-05-15 2019-11-21 株式会社竹中土木 トンネル切羽監視方法
CN110516907A (zh) * 2019-07-17 2019-11-29 吉林大学 一种基于ahp-熵权云模型的岩爆等级评价方法
CN110826223A (zh) * 2019-11-05 2020-02-21 山东省交通科学研究院 基于综合属性测度的岩爆危险性预测方法、系统及介质
CN111582748A (zh) * 2020-05-15 2020-08-25 国网陕西省电力公司电力科学研究院 一种基于评估指标体系和组合权重的黑启动方案评估方法
CN111598448A (zh) * 2020-05-15 2020-08-28 青岛理工大学 一种基于混凝土t型梁火灾后损伤评估方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104008302A (zh) * 2014-06-09 2014-08-27 上海电力学院 一种基于组合赋权和模糊评分的配电网可靠性评估方法
RU2560008C1 (ru) * 2014-06-16 2015-08-20 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Национальный минерально-сырьевой университет "Горный" Способ крепления горных выработок на удароопасных пластах
CN104732070A (zh) * 2015-02-27 2015-06-24 广西大学 一种基于信息向量机的岩爆等级预测方法
CN106407493A (zh) * 2016-03-15 2017-02-15 中南大学 一种基于多维高斯云模型的岩爆等级评价方法
CN107748103A (zh) * 2017-09-01 2018-03-02 中国科学院武汉岩土力学研究所 一种隧道岩爆预测方法、设备、存储介质和系统
JP2019199708A (ja) * 2018-05-15 2019-11-21 株式会社竹中土木 トンネル切羽監視方法
CN109523200A (zh) * 2019-01-03 2019-03-26 东北大学 基于最小二乘法的铸造车间安全评价综合权重计算方法
CN109934500A (zh) * 2019-03-15 2019-06-25 南方电网科学研究院有限责任公司 一种配电台区的低电压评价方法、装置、设备及存储介质
CN110516907A (zh) * 2019-07-17 2019-11-29 吉林大学 一种基于ahp-熵权云模型的岩爆等级评价方法
CN110826223A (zh) * 2019-11-05 2020-02-21 山东省交通科学研究院 基于综合属性测度的岩爆危险性预测方法、系统及介质
CN111582748A (zh) * 2020-05-15 2020-08-25 国网陕西省电力公司电力科学研究院 一种基于评估指标体系和组合权重的黑启动方案评估方法
CN111598448A (zh) * 2020-05-15 2020-08-28 青岛理工大学 一种基于混凝土t型梁火灾后损伤评估方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
孟凡奇 等: "最优组合赋权法在泥石流危险度评价中的应用", 《人民长江》 *
李明亮 等: "基于变异系数与序关分析法-多维正态云模型的岩爆预测", 《岩石力学与工程学报》 *
王正选 等: "基于改进熵权法的水资源承载力评价-以曲陆坝区为例", 《水资源与水工程学报》 *
王迎超 等: "基于正态云模型的深埋地下工程岩爆烈度分级预测研究", 《岩土力学》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113076700A (zh) * 2021-04-27 2021-07-06 昆明理工大学 一种基于数据分析原理的svm-lda岩爆机器学习预测模型方法
CN113221471A (zh) * 2021-06-15 2021-08-06 贵州大学 基于改进可变模糊集与神经网络的岩爆等级判别方法
CN115526036A (zh) * 2022-09-19 2022-12-27 长安大学 一种岩爆规模等级的判别方法及系统
CN116843228A (zh) * 2023-07-21 2023-10-03 昆明理工大学 基于模糊集理论和组合赋权的煤层冲击倾向性评价方法
CN116843228B (zh) * 2023-07-21 2024-05-14 昆明理工大学 基于模糊集理论和组合赋权的煤层冲击倾向性评价方法
CN117272644A (zh) * 2023-09-21 2023-12-22 安徽理工大学 基于钻孔参数的岩爆倾向性预测方法及系统
CN117272644B (zh) * 2023-09-21 2024-03-05 安徽理工大学 基于钻孔参数的岩爆倾向性预测方法及系统
CN117110584A (zh) * 2023-10-23 2023-11-24 中铁四局集团有限公司 一种应变型隧道岩爆烈度等级预测方法
CN117110584B (zh) * 2023-10-23 2024-03-26 中铁四局集团有限公司 一种应变型隧道岩爆烈度等级预测方法
CN117332240A (zh) * 2023-12-01 2024-01-02 中铁四局集团有限公司 岩爆预测模型构建方法、存储介质、岩爆预测方法及系统
CN117332240B (zh) * 2023-12-01 2024-04-16 中铁四局集团有限公司 岩爆预测模型构建方法、存储介质、岩爆预测方法及系统

Also Published As

Publication number Publication date
CN112699553B (zh) 2021-11-02

Similar Documents

Publication Publication Date Title
CN112699553B (zh) 一种岩爆烈度等级智能预测系统方法
CN111582718B (zh) 基于网络层次分析法的电缆通道火灾风险评估方法及装置
CN113792982A (zh) 基于组合赋权与模糊灰色聚类的科技服务质量评估方法和装置
Butler et al. A multiple attribute utility theory approach to ranking and selection
CN104461896B (zh) 基于可信属性的航天系统关键软件评价方法
Lahdelma et al. Stochastic multicriteria acceptability analysis (SMAA)
CN111797364B (zh) 一种基于复合云模型的滑坡多层次安全评价方法
CN105205329A (zh) 一种大坝安全综合评价方法
CN111967721A (zh) 一种综合能源系统绿色化水平评价方法及系统
CN114065223A (zh) 一种基于cvss的多维度软件安全风险评估方法
CN110889440A (zh) 基于主成分分析和bp神经网络的岩爆等级预测方法及系统
Obsie et al. Prediction of student academic performance using neural network, linear regression and support vector regression: a case study
Morrice et al. An approach to ranking and selection for multiple performance measures
CN114444910A (zh) 一种面向电力物联网的边缘网络系统健康度评估方法
CN115063056A (zh) 基于图拓扑分析改进的建造行为安全风险动态分析方法
CN110472363B (zh) 适用于高速铁路隧道的围岩变形等级预测方法及系统
Shi et al. A Grey Model for Evaluation of Information Systems Security.
CN114862267A (zh) 一种输油气管道报警管理体系的评价方法及系统
CN117745096A (zh) 一种控制网布测方案评价方法及系统
CN113779893B (zh) 基于回声状态网络的保障性仿真输出结果可信性评价方法
CN113177748B (zh) 一种输气管道震害评价方法
CN113793017A (zh) 一种配电房老化评价方法、装置、终端设备及存储介质
CN113191684A (zh) 基于变权模糊综合评判的北斗卫星合格评判方法
CN113221332B (zh) 一种基于云模型理论的海岸侵蚀脆弱性评估方法
CN112116238A (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
GR01 Patent grant
GR01 Patent grant