CN111784136A - A Dynamic Early Warning Method of Impact Risk Based on Analytic Hierarchy Process and Fuzzy Mathematics - Google Patents

A Dynamic Early Warning Method of Impact Risk Based on Analytic Hierarchy Process and Fuzzy Mathematics Download PDF

Info

Publication number
CN111784136A
CN111784136A CN202010577804.3A CN202010577804A CN111784136A CN 111784136 A CN111784136 A CN 111784136A CN 202010577804 A CN202010577804 A CN 202010577804A CN 111784136 A CN111784136 A CN 111784136A
Authority
CN
China
Prior art keywords
early warning
fuzzy
impact
index
early
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
Application number
CN202010577804.3A
Other languages
Chinese (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.)
China Coal Research Institute CCRI
Original Assignee
China Coal Research Institute CCRI
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 China Coal Research Institute CCRI filed Critical China Coal Research Institute CCRI
Priority to CN202010577804.3A priority Critical patent/CN111784136A/en
Publication of CN111784136A publication Critical patent/CN111784136A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0635Risk analysis of enterprise or organisation activities
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Animal Husbandry (AREA)
  • Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Agronomy & Crop Science (AREA)
  • Alarm Systems (AREA)

Abstract

本发明公开一种基于层次分析和模糊数学的冲击危险性动态预警方法,该预警方法包括以下步骤:步骤一、确定预警对象和预警对象适用的冲击危险性影响因子,即预警指标;步骤二、确定各预警指标的冲击危险性分级,组成冲击危险性预警指标体系;步骤三、基于APH(层次分析)理论,确定各预警指标的权重,组成模糊权重向量A;步骤四、基于模糊数学理论,建立各指标的单因素模糊预警集,组成因素预警模糊矩阵R;步骤五、计算预警对象的模糊综合评判向量B=A×R;步骤六、根据最大隶属度原则,确定预警对象的冲击危险等级。本发明综合考虑了多种冲击地压预警指标,采用预警模型对预警对象进行冲击危险性预警,预警更加科学,预警效率更高,结果更加可靠。

Figure 202010577804

The invention discloses a dynamic early warning method for shock risk based on AHP and fuzzy mathematics. The early warning method includes the following steps: Step 1: Determine the early warning object and the shock risk influencing factor applicable to the early warning object, that is, the early warning index; Step 2, Determine the impact risk classification of each early warning index, and form a shock risk early warning index system; step 3, based on the APH (analytic hierarchy process) theory, determine the weight of each early warning index, and form a fuzzy weight vector A; step 4, based on fuzzy mathematics theory, Establish a single-factor fuzzy early-warning set of each index, which is composed of a factor early-warning fuzzy matrix R; Step 5: Calculate the fuzzy comprehensive evaluation vector B=A×R of the early-warning object; Step 6. Determine the impact risk level of the early-warning object according to the principle of maximum membership degree . The present invention comprehensively considers a variety of rockburst early warning indicators, adopts the early warning model to carry out the early warning of the shock danger to the early warning object, the early warning is more scientific, the early warning efficiency is higher, and the result is more reliable.

Figure 202010577804

Description

一种基于层次分析和模糊数学的冲击危险性动态预警方法A Dynamic Early Warning Method of Impact Risk Based on Analytic Hierarchy Process and Fuzzy Mathematics

技术领域technical field

本发明涉及煤矿安全技术领域,特别是涉及一种基于层次分析和模糊数学的冲击危险性动态预警方法。The invention relates to the technical field of coal mine safety, in particular to a dynamic early warning method for impact risk based on AHP and fuzzy mathematics.

背景技术Background technique

随着我国煤矿开采深度的增加以及开采强度的增大,煤矿面临冲击地压的形势越来越严峻,而冲击地压评价及实时监测预警是冲击地压防治的有效措施。同时,2018年发布的《防治煤矿冲击地压细则》中规定:开采冲击地压煤层时,必须采取冲击地压危险性预测、监测预警、防范治理、效果检验、安全防护等综合性防治措施;冲击地压矿井必须建立区域与局部相结合的冲击危险性监测制度,区域监测应当覆盖矿井采掘区域,局部监测应当覆盖冲击地压危险区,区域监测可采用微震监测法等,局部监测可采用钻屑法、应力监测法、电磁辐射法等。但目前我国冲击危险性预警大多采用单一指标或结合简单模型进行多指标预警,预警的效率及准确率已满足不了煤矿的现实需求。With the increase of mining depth and mining intensity of coal mines in my country, the situation of coal mines facing rockburst is becoming more and more severe, and rockburst evaluation and real-time monitoring and early warning are effective measures to prevent rockburst. At the same time, the "Detailed Rules for the Prevention and Control of Rock Burst in Coal Mine" issued in 2018 stipulates that when mining rock burst coal seams, comprehensive prevention measures such as rock burst risk prediction, monitoring and early warning, prevention and treatment, effect inspection, and safety protection must be taken; Rockburst mines must establish a regional and local impact risk monitoring system. Regional monitoring should cover the mining area of the mine, and local monitoring should cover the rockburst danger zone. Regional monitoring can use microseismic monitoring, and local monitoring can use drilling. Chip method, stress monitoring method, electromagnetic radiation method, etc. However, at present, most of my country's shock risk early warning adopts a single indicator or a combination of simple models for multi-indicator early warning, and the efficiency and accuracy of early warning cannot meet the practical needs of coal mines.

发明内容SUMMARY OF THE INVENTION

本发明的目的在于提供一种基于层次分析和模糊数学的冲击危险性动态预警方法,该方法将与冲击地压密切相关的多种参量进行综合考虑,并结合数学模型及判别准则进行实时预警,预警效率更高、结果更加可靠,对指导矿井开展冲击地压防治工作实用性更强。The purpose of the present invention is to provide a dynamic early warning method of shock risk based on AHP and fuzzy mathematics. The early warning is more efficient, the results are more reliable, and it is more practical for guiding mines to carry out rockburst prevention and control.

本发明所采用的技术解决方案是:The technical solution adopted by the present invention is:

一种基于层次分析和模糊数学的冲击危险性动态预警方法,该预警方法包括以下步骤:A dynamic early warning method for impact risk based on AHP and fuzzy mathematics, the early warning method includes the following steps:

步骤一、确定预警对象和预警对象适用的冲击危险性预警指标Step 1. Determine the early warning object and the shock risk early warning index applicable to the early warning object

根据预警对象选取其适用的冲击危险性预警指标集合U,用ui表示各预警指标,则U=(u1,u2,u3,u4…ui,…un),其中n为预警指标数量;Select the applicable shock risk early warning index set U according to the early warning object, and use u i to represent each early warning index, then U=(u 1 , u 2 , u 3 , u 4 ... u i , ... u n ), where n is Number of early warning indicators;

步骤二、对各预警指标的冲击危险性分级,组成冲击危险性预警指标体系Step 2: Classify the impact risk of each early warning index to form a shock risk early warning index system

根据预警等级分类,对各预警指标的冲击危险性分级;通过建立每个预警指标的分级集合vi,得到预警对象的冲击危险性预警指标体系V;According to the classification of early warning levels, the impact risk of each early warning index is classified; by establishing the classification set v i of each early warning index, the impact risk early warning index system V of the early warning object is obtained;

步骤三、确定各预警指标的权重,组成模糊权重向量Step 3: Determine the weight of each early warning indicator to form a fuzzy weight vector

根据APH层次分析法和1~9标度表,确定出各预警指标对预警对象冲击危险的贡献程度,进而得到各个指标组成的模糊权重向量A;According to the APH analytic hierarchy process and the 1-9 scale table, the contribution degree of each early warning index to the impact risk of the early warning object is determined, and then the fuzzy weight vector A composed of each index is obtained;

步骤四、建立各预警指标的单因素模糊预警集,组成因素预警模糊矩阵Step 4. Establish a single-factor fuzzy early-warning set of each early-warning index, and form a factor early-warning fuzzy matrix

针对某个预警指标的监测结果,结合预警隶属度函数,可确定该预警指标的全部隶属度的集合;进而依次可得到相应于每个预警指标的单因素模糊预警集;将第i个单因素模糊预警集Ri为第i行,组成因素预警模糊矩阵R;According to the monitoring results of a certain early warning indicator, combined with the early warning membership function, the set of all membership degrees of the early warning indicator can be determined; then the single factor fuzzy early warning set corresponding to each early warning indicator can be obtained in turn; The fuzzy early warning set R i is the i-th row, and the composition factor early warning fuzzy matrix R;

步骤五、计算预警对象的模糊综合评判向量Step 5. Calculate the fuzzy comprehensive evaluation vector of the early warning object

将反映各预警指标对冲击地压影响程度的模糊权重向量A与全部指标隶属度合成的因素预警模糊矩阵R,进行运算,从而得出每个冲击地压预警的模糊综合评判向量B=A×R=(b1,b2,b3,b4);Calculate the fuzzy weight vector A that reflects the influence degree of each early warning index on rock burst and the factor early warning fuzzy matrix R, which is composed of the membership degrees of all indicators, so as to obtain the fuzzy comprehensive evaluation vector of each rock burst early warning B=A× R=(b 1 , b 2 , b 3 , b 4 );

步骤六、确定预警对象的冲击危险等级Step 6: Determine the impact hazard level of the warning object

步骤五中(b1,b2,b3,b4)分别对应于(无冲击危险,弱冲击危险,中等冲击危险,强冲击危险)四个评判结果的隶属度;根据最大隶属度原则,就可以得出预警对象的冲击危险等级。In step 5 (b 1 , b 2 , b 3 , b 4 ) respectively correspond to the membership degrees of the four evaluation results (no impact risk, weak impact risk, medium impact risk, and strong impact risk); according to the principle of maximum membership degree, The impact risk level of the warning object can be obtained.

上述步骤二中,各预警指标的冲击危险性可分为5级,分别为无、弱、中等、强和不安全。步骤二中各预警指标的冲击危险性分级数量与步骤六中最终预警对象的冲击危险等级划分数量不一致,但并不影响该方法的应用效果。In the above step 2, the impact risk of each early warning indicator can be divided into 5 levels, namely none, weak, medium, strong and unsafe. The number of impact risk classifications of each early warning index in step 2 is inconsistent with the number of impact risk classifications of the final early warning object in step 6, but it does not affect the application effect of the method.

上述步骤三中,通过将各预警指标ui两两对比并对应1~9标度表进行赋值,可得到矩阵,针对矩阵求特征向量,可得到模糊权重向量A=(a1,a2,…ai,…an),归一化处理后为(a’1,a’2,…a’i,…a’n)。In the above-mentioned step 3, a matrix can be obtained by comparing the early warning indicators u i pairwise and assigning values corresponding to the scale tables 1 to 9, and the eigenvectors of the matrix can be obtained to obtain a fuzzy weight vector A=(a 1 , a 2 , ...a i , ...a n ), normalized to (a' 1 , a' 2 , ...a' i , ...a' n ).

上述步骤四中,所述隶属度函数如公式(1)所示:In the above step 4, the membership function is shown in formula (1):

Figure BDA0002549941610000021
Figure BDA0002549941610000021

Figure BDA0002549941610000022
Figure BDA0002549941610000022

Figure BDA0002549941610000023
Figure BDA0002549941610000023

Figure BDA0002549941610000024
Figure BDA0002549941610000024

式(1)中,x为测试值,S1、S2、S3和S4分别对应预警指标的四个分级阈值。In formula (1), x is the test value, and S 1 , S 2 , S 3 and S 4 correspond to the four grading thresholds of the early warning indicators, respectively.

上述步骤六中,所述最大隶属度原则是指,通过对比b1~b4大小,若最大值为bi,i取值1、2、3、4,则预警对象为第i级冲击危险性;冲击危险等级分为四类:冲击危险等级属于1级,为无冲击危险;冲击危险等级属于2级,为弱冲击危险;冲击危险等级属于3级,为中等冲击危险;冲击危险等级属于4级,为强冲击危险。In the above-mentioned step 6, the principle of the maximum membership degree means that, by comparing the sizes of b 1 to b 4 , if the maximum value is b i , and the value of i is 1, 2, 3, or 4, then the warning object is the ith level impact danger. The impact hazard level is divided into four categories: the impact hazard level belongs to level 1, which means no impact hazard; the impact hazard level belongs to level 2, which is a weak impact hazard; the impact hazard level belongs to level 3, which is a medium impact hazard; the impact hazard level belongs to Level 4 is a strong shock hazard.

本发明的有益技术效果是:The beneficial technical effects of the present invention are:

本发明综合考虑了多种冲击地压预警指标,根据不同工况,通过数学模型赋予不同的权重值,然后采用预警模型对预警对象进行冲击危险性预警,预警更加科学,预警效率更高,结果更加可靠。另外,该预警方法适用于工程中实时监测的预警装备,即冲击危险性预警结果是实时地、动态的。本发明预警方法可用于指导煤矿开展冲击地压防治工作,对于保障煤矿安全生产具有重要的现实意义。The invention comprehensively considers a variety of rockburst early warning indicators, assigns different weight values through a mathematical model according to different working conditions, and then uses the early warning model to carry out early warning on the impact risk of the early warning object, the early warning is more scientific, and the early warning efficiency is higher. more reliable. In addition, this early warning method is suitable for early warning equipment for real-time monitoring in engineering, that is, the impact risk early warning results are real-time and dynamic. The early warning method of the invention can be used to guide coal mines to carry out rockburst prevention and control work, and has important practical significance for ensuring the safe production of coal mines.

附图说明Description of drawings

下面结合附图和具体实施方式对本发明作进一步说明:The present invention will be further described below in conjunction with the accompanying drawings and specific embodiments:

图1为本发明基于层次分析和模糊数学的冲击危险性动态预警方法的流程图;Fig. 1 is the flow chart of the dynamic early warning method of impact risk based on AHP and fuzzy mathematics of the present invention;

图2为本发明具体应用实例中所建立冲击危险性预警指标体系的示意图。FIG. 2 is a schematic diagram of a shock risk early warning index system established in a specific application example of the present invention.

具体实施方式Detailed ways

如图1所示,一种基于层次分析和模糊数学的冲击危险性动态预警方法,该预警方法适用于煤矿,具体包括以下步骤:As shown in Figure 1, a dynamic early warning method of shock risk based on AHP and fuzzy mathematics, the early warning method is suitable for coal mines, and includes the following steps:

步骤一、确定预警对象和预警对象适用的冲击危险性影响因子,即预警指标Step 1. Determine the early warning object and the impact risk factor applicable to the early warning object, that is, the early warning index

预警对象可以是一个范围,包括矿井、采区(盘区或水平)、回采工作面和巷道,也可以是一段巷道、一个监测点。以监测点为例,根据预警对象选取其适用的冲击危险性预警指标集合U,用ui表示各预警指标,则U=(u1,u2,u3,u4…ui,…un),其中n为预警指标数量,包括“采动应力”、“巷道位移”、“顶板离层量”、“锚杆工作载荷”、“锚索工作载荷”和“支架阻力”等。The object of early warning can be a range, including mine, mining area (panel or level), mining face and roadway, or a section of roadway and a monitoring point. Take the monitoring point as an example, select the applicable shock risk early warning index set U according to the early warning object, and use u i to represent each early warning index, then U=(u 1 , u 2 , u 3 , u 4 ... u i , ... u n ), where n is the number of early warning indicators, including "mining stress", "roadway displacement", "roof separation amount", "anchor bolt working load", "anchor cable working load" and "support resistance", etc.

步骤二、对各预警指标的冲击危险性分级,组成冲击危险性预警指标体系Step 2: Classify the impact risk of each early warning index to form a shock risk early warning index system

根据预警等级分类,对各预警指标的冲击危险性分级,共分为5级,分别为无、弱、中等、强和不安全。以采动应力σ为例,预警指标分级集合v1=(无,弱,中等,强,不安全)=(σ≤1.5σ正常,1.5σ正常<σ≤3σ正常,3σ正常<σ≤4.5σ正常,4.5σ正常<σ≤6σ正常,σ>6σ正常)。通过建立每个预警指标的分级集合vi,得到预警对象如监测点的冲击危险性预警指标体系V。According to the classification of early warning levels, the impact risk of each early warning index is classified into 5 levels, namely none, weak, medium, strong and unsafe. Taking mining stress σ as an example, the classification set of early warning indicators v 1 = (none, weak, moderate, strong, unsafe) = (σ≤1.5σ normal , 1.5σ normal <σ≤3σ normal , 3σ normal <σ≤4.5 σ is normal , 4.5σ is normal <σ≤6σ is normal , σ>6σ is normal ). By establishing the graded set vi of each early warning index, the early warning index system V of shock risk of early warning objects such as monitoring points is obtained.

步骤三、确定各预警指标的权重,组成模糊权重向量Step 3: Determine the weight of each early warning indicator to form a fuzzy weight vector

权重表征每个预警指标对监测点冲击危险的贡献程度,如采动应力、巷道位移等分别对监测点冲击危险性影响所占的权重。The weight represents the contribution of each early warning index to the impact risk of the monitoring point, such as the weight of mining stress, roadway displacement, etc., respectively, on the impact risk of the monitoring point.

根据以往采动应力、巷道位移等监测数据与冲击显现的对应关系,再结合APH层次分析法和1~9标度表,确定出各预警指标对监测点冲击危险的贡献程度,进而得到各个指标组成的模糊权重向量A。According to the corresponding relationship between the previous mining stress, roadway displacement and other monitoring data and the impact, combined with the APH analytic hierarchy process and the 1-9 scale table, the contribution of each early warning index to the impact risk of the monitoring point is determined, and then each index is obtained. Composed of fuzzy weight vector A.

步骤四、建立各预警指标的单因素模糊预警集,组成因素预警模糊矩阵Step 4. Establish a single-factor fuzzy early-warning set of each early-warning index, and form a factor early-warning fuzzy matrix

应用模糊数学,确定每一个预警指标隶属于预警等级集合中不同预警等级的程度,称为隶属度,以rij表示。针对某个预警指标的监测结果(实际测试结果),结合预警隶属度函数,可确定该预警指标的全部隶属度的集合,即为单因素预警矩阵;进而依次可得到相应于每个预警指标的单因素模糊预警矩阵或者说预警集。根据不同的预警指标,从上至下按每行排列,将第i个单因素模糊预警集Ri为第i行,组成因素预警模糊矩阵R。如该因素预警模糊矩阵R可为n行4列。Apply fuzzy mathematics to determine the degree to which each early warning index belongs to different early warning levels in the early warning level set, which is called membership degree, and is represented by r ij . According to the monitoring results (actual test results) of a certain early warning indicator, combined with the early warning membership function, the set of all membership degrees of the early warning indicator can be determined, which is a single-factor early warning matrix; and then the corresponding early warning indicators can be obtained in turn. Single factor fuzzy early warning matrix or early warning set. According to different early warning indicators, arrange each row from top to bottom, take the i-th single-factor fuzzy early-warning set R i as the i-th row, and make up the factor early-warning fuzzy matrix R. For example, the factor warning fuzzy matrix R can be n rows and 4 columns.

步骤五、计算预警对象的模糊综合评判向量Step 5. Calculate the fuzzy comprehensive evaluation vector of the early warning object

模糊综合预警考虑所有因素对监测点的影响,将反映各预警指标对冲击地压影响程度的模糊权重向量A与全部指标隶属度合成的因素预警模糊矩阵R,进行运算,从而得出每个冲击地压预警的模糊综合评判向量B=A×R=(b1,b2,b3,b4)。The fuzzy comprehensive early warning considers the influence of all factors on the monitoring points, and calculates the fuzzy weight vector A that reflects the influence degree of each early warning index on rockburst and the factor early warning fuzzy matrix R synthesized by the membership degrees of all indicators, so as to obtain each impact Fuzzy comprehensive evaluation vector of ground pressure warning B=A×R=(b 1 , b 2 , b 3 , b 4 ).

步骤六、确定预警对象的冲击危险等级Step 6: Determine the impact hazard level of the warning object

步骤五中(b1,b2,b3,b4)分别对应于(无冲击危险,弱冲击危险,中等冲击危险,强冲击危险)四个评判结果的隶属度。根据最大隶属度原则,就可以得出监测点一段时间的监测数据所反映的冲击危险等级。In step 5, (b 1 , b 2 , b 3 , b 4 ) respectively correspond to the membership degrees of the four judgment results (no shock risk, weak shock risk, medium shock risk, and strong shock risk). According to the principle of maximum membership degree, the impact risk level reflected by the monitoring data of the monitoring point for a period of time can be obtained.

上述步骤三中,1~9标度表见下表1。In the above-mentioned step 3, the scale table of 1 to 9 is shown in Table 1 below.

表1Table 1

标度Scaling 表示含义express meaning 11 元素i和元素j同等重要Element i and element j are equally important 33 元素i比元素j稍微重要Element i is slightly more important than element j 55 元素i比元素j明显重要Element i is significantly more important than element j 77 元素i比元素j强烈重要Element i is strongly more important than element j 99 元素i比元素j极端重要Element i is extremely important than element j 2,4,6,82,4,6,8 介于上述标度之间between the above scales 倒数reciprocal 表示a<sub>ij</sub>与a<sub>ij</sub>互为倒数Indicates that a<sub>ij</sub> and a<sub>ij</sub> are reciprocals of each other

通过将各预警指标ui两两对比并对应表1进行赋值,可得到矩阵G,针对矩阵G求特征向量,可得到模糊权重矩阵A=(a1,a2,…ai,…an),归一化处理后为(a’1,a’2,…a’i,…a’n)。By comparing the early warning indicators u i pairwise and assigning them to Table 1, the matrix G can be obtained, and the eigenvectors of the matrix G can be obtained to obtain the fuzzy weight matrix A=(a 1 , a 2 ,...a i ,...a n ), normalized to (a' 1 , a' 2 ,...a' i ,...a' n ).

步骤四中,隶属度rij通过隶属度函数计算求得,隶属度函数如公式(1)所示。In step 4, the membership degree r ij is obtained by calculating the membership degree function, and the membership degree function is shown in formula (1).

Figure BDA0002549941610000051
Figure BDA0002549941610000051

Figure BDA0002549941610000052
Figure BDA0002549941610000052

Figure BDA0002549941610000053
Figure BDA0002549941610000053

式(1)中,x为测试值,S1、S2、S3和S4分别对应预警指标的四个分级阈值。In formula (1), x is the test value, and S 1 , S 2 , S 3 and S 4 correspond to the four grading thresholds of the early warning indicators, respectively.

步骤六中,所述最大隶属度原则是指,通过对比b1~b4大小,若最大值bi,i取值为1、2、3或4,则预警对象为第i级冲击危险性。冲击危险等级分为四类;冲击危险等级属于1级,为无冲击危险;冲击危险等级属于2级,为弱冲击危险;冲击危险等级属于3级,为中等冲击危险;冲击危险等级属于4级,为强冲击危险。In step 6, the maximum membership degree principle means that, by comparing the sizes of b 1 to b 4 , if the maximum value b i , i takes a value of 1, 2, 3 or 4, the warning object is the ith level impact risk. . The impact hazard level is divided into four categories; the impact hazard level belongs to level 1, which means no impact hazard; the impact hazard level belongs to level 2, which is a weak impact hazard; the impact hazard level belongs to level 3, which is a medium impact hazard; the impact hazard level belongs to level 4 , for a strong shock hazard.

下面以河南某煤矿为具体应用实例,进一步介绍本发明基于层次分析和模糊数学的冲击危险性动态预警方法。The following takes a coal mine in Henan as a specific application example to further introduce the dynamic early warning method of shock risk based on AHP and fuzzy mathematics.

背景:河南某煤矿为典型的冲击地压矿井,在开采2-3煤过程中曾发生多次冲击地压事故。该矿13200工作面上巷掘进过程中安装了应力计和位移计进行实时监测,现采用该方法对该巷道掘进过程中的冲击危险性进行实时预警。Background: A coal mine in Henan is a typical rockburst mine. There have been many rockburst accidents during the mining of 2-3 coal. The stress meter and displacement meter are installed for real-time monitoring during the excavation of the roadway on the 13200 working face of the mine, and this method is now used to provide real-time early warning of the impact risk during the excavation of the roadway.

(1)确定预警对象和预警对象适用的冲击危险性预警指标:根据应力计和位移计安装位置,确定预警对象13200上巷距入口500~600m段,确定冲击危险性预警指标为应力梯度▽σ和位移梯度▽ε。(1) Determine the early warning object and the impact risk early warning index applicable to the early warning object: According to the installation position of the stress meter and displacement meter, determine the 500-600m section of the early warning object 13200 upper lane from the entrance, and determine the impact risk early warning index as the stress gradient ▽σ and the displacement gradient ▽ε.

(2)确定各预警指标的冲击危险性分级,组成冲击危险性预警指标体系,如图2所示。(2) Determine the impact risk classification of each early warning index, and form a shock risk early warning index system, as shown in Figure 2.

(3)确定各预警指标的权重,组成模糊权重向量:根据各预警指标之间的重要程度对比,结合APH理论和1~9标度表,可得到模糊权重矩阵A=(a1’,a2’)=(0.1667,0.8333)。(3) Determine the weight of each early warning index and form a fuzzy weight vector: according to the comparison of the importance of each early warning index, combined with the APH theory and the 1-9 scale table, the fuzzy weight matrix A = (a 1 ', a 2 ')=(0.1667, 0.8333).

(4)建立各指标的单因素模糊预警集,组成因素预警模糊矩阵。其中各预警指标的分级阈值参考图2,计算得到因素预警模糊矩阵R为(4) Establish a single-factor fuzzy early-warning set of each index, and form a factor early-warning fuzzy matrix. The grading thresholds of each early warning index refer to Figure 2, and the factor early warning fuzzy matrix R is calculated as

Figure BDA0002549941610000061
Figure BDA0002549941610000061

(5)计算预警对象的模糊综合评判向量:B=A×R=(0.9603,0.0397,0,0)。(5) Calculate the fuzzy comprehensive evaluation vector of the warning object: B=A×R=(0.9603, 0.0397, 0, 0).

(6)确定预警对象的冲击危险等级:由最大隶属度原则得出,该测点的模糊综合评判指数为0.9603,冲击危险等级为无。(6) Determine the impact risk level of the early warning object: According to the principle of maximum membership degree, the fuzzy comprehensive evaluation index of this measuring point is 0.9603, and the impact risk level is none.

根据现场钻屑量监测结果及动力显现情况可知,监测时间段内钻屑量未超标,且现场无动力显现,可见该监测点为无冲击危险,与基于层次分析和模糊数学的冲击危险性动态预警方法所得结果一致。According to the monitoring results of on-site drilling cuttings volume and the dynamic display situation, it can be seen that the drilling cuttings volume did not exceed the standard during the monitoring time period, and there was no dynamic display on the site. It can be seen that the monitoring point has no impact risk, which is different from the dynamic impact risk based on AHP and fuzzy mathematics. The results of the early warning method were consistent.

上述方式中未述及的部分采取或借鉴已有技术即可实现。The parts not mentioned in the above manner can be realized by adopting or learning from the existing technology.

当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。Of course, the above description is not intended to limit the present invention, and the present invention is not limited to the above examples. Changes, modifications, additions or substitutions made by those skilled in the art within the essential scope of the present invention should also belong to the present invention. the scope of protection of the invention.

Claims (4)

1. A dynamic early warning method for impact danger based on hierarchical analysis and fuzzy mathematics is characterized by comprising the following steps:
step one, determining an early-warning object and an impact risk early-warning index applicable to the early-warning object
Selecting an applicable impact risk early warning index set U according to the early warning object, and using UiIf each warning indicator is represented, U is equal to (U)1,u2,u3,u4…ui,…un) Wherein n is the number of early warning indicators;
step two, classifying the impact dangerousness of each early warning index to form an impact danger early warning index system
Classifying according to the early warning grades, and grading the impact risks of the early warning indexes; by establishing a hierarchical set v of each pre-warning indicatoriObtaining an impact risk early warning index system V of the early warning object;
step three, determining the weight of each early warning index to form a fuzzy weight vector
Determining the contribution degree of each early warning index to the impact risk of the early warning object according to an APH (advanced persistent threat) analytic hierarchy process and a 1-9 standard table, and further obtaining a fuzzy weight vector A formed by each index;
step four, establishing a single-factor fuzzy early warning set of each early warning index to form a factor early warning fuzzy matrix
Aiming at the monitoring result of a certain early warning index, a set of all membership degrees of the early warning index can be determined by combining an early warning membership function; and then a single-factor fuzzy early warning set corresponding to each early warning index can be obtained in sequence; the ith single-factor fuzzy early warning set RiForming a factor early warning fuzzy matrix R for the ith row;
step five, calculating a fuzzy comprehensive judgment vector of the early warning object
Calculating a factor early warning fuzzy matrix R synthesized by a fuzzy weight vector A reflecting the influence degree of each early warning index on the rock burst and all index membership degrees, thereby obtaining a fuzzy comprehensive judgment vector B of each rock burst early warning, wherein the fuzzy comprehensive judgment vector B is A × R (B)1,b2,b3,b4);
Step six, determining the impact danger level of the early warning object
In step five (b)1,b2,b3,b4) Respectively corresponding to the membership degrees of the four judgment results (no impact risk, weak impact risk, medium impact risk and strong impact risk); according to the maximum membership principle, the impact danger level of the early warning object can be obtained.
2. The dynamic early warning method for impact risk based on hierarchical analysis and fuzzy mathematics as claimed in claim 1, wherein: in the third step, each early warning index u is usediPairwise comparison is carried out and the value is assigned corresponding to a 1-9 scale table, a matrix can be obtained, the characteristic vector is solved aiming at the matrix, and a fuzzy weight vector A (a) is obtained1,a2,…ai,…an) And is (a ') after normalization treatment'1,a’2,…a’i,…a’n)。
3. The dynamic early warning method for impact risk based on hierarchical analysis and fuzzy mathematics as claimed in claim 1, wherein: in the fourth step, the membership function is shown as formula (1):
Figure FDA0002549941600000021
in the formula (1), x is a test value, S1、S2、S3And S4And the four grading threshold values respectively correspond to the early warning indexes.
4. The dynamic early warning method for impact risk based on hierarchical analysis and fuzzy mathematics as claimed in claim 1, wherein: in the sixth step, the principle of maximum membership degree means that b is compared1~b4Size, if maximum value is biIf i takes values of 1, 2, 3 and 4, the early warning object is the ith-level impact risk; impact hazard grade fractionIs classified into four types: the impact danger level belongs to level 1 and is no impact danger; the impact risk grade belongs to grade 2 and is a weak impact risk; the impact risk rating is class 3, which is a medium impact risk; the impact hazard classification belongs to class 4 and is a strong impact hazard.
CN202010577804.3A 2020-06-22 2020-06-22 A Dynamic Early Warning Method of Impact Risk Based on Analytic Hierarchy Process and Fuzzy Mathematics Pending CN111784136A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010577804.3A CN111784136A (en) 2020-06-22 2020-06-22 A Dynamic Early Warning Method of Impact Risk Based on Analytic Hierarchy Process and Fuzzy Mathematics

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010577804.3A CN111784136A (en) 2020-06-22 2020-06-22 A Dynamic Early Warning Method of Impact Risk Based on Analytic Hierarchy Process and Fuzzy Mathematics

Publications (1)

Publication Number Publication Date
CN111784136A true CN111784136A (en) 2020-10-16

Family

ID=72757645

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010577804.3A Pending CN111784136A (en) 2020-06-22 2020-06-22 A Dynamic Early Warning Method of Impact Risk Based on Analytic Hierarchy Process and Fuzzy Mathematics

Country Status (1)

Country Link
CN (1) CN111784136A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115496342A (en) * 2022-09-05 2022-12-20 煤炭科学技术研究院有限公司 A method and device for early warning of rock burst based on subjective and objective dynamic weights
CN115860582A (en) * 2023-02-28 2023-03-28 山东科技大学 Intelligent impact risk early warning method based on self-adaptive lifting algorithm
CN116227982A (en) * 2022-12-30 2023-06-06 中国矿业大学(北京) Quantification method and device for pollution degree of coal dust

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张宁博 等: "基于一孔多点式应力与位移监测系统的掘进巷道冲击危险性评价", 煤炭学报, pages 140 - 149 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115496342A (en) * 2022-09-05 2022-12-20 煤炭科学技术研究院有限公司 A method and device for early warning of rock burst based on subjective and objective dynamic weights
CN116227982A (en) * 2022-12-30 2023-06-06 中国矿业大学(北京) Quantification method and device for pollution degree of coal dust
CN116227982B (en) * 2022-12-30 2023-10-31 中国矿业大学(北京) Quantification method and device for pollution degree of coal dust
CN115860582A (en) * 2023-02-28 2023-03-28 山东科技大学 Intelligent impact risk early warning method based on self-adaptive lifting algorithm

Similar Documents

Publication Publication Date Title
CN109034582B (en) Tunnel crossing fault water inrush and mud inrush risk evaluation method based on cloud model and combined weighting
CN105785471B (en) A kind of shock hazard evaluation method of the pre- working seam of mine
CN113222347B (en) Method for evaluating gray system of safety risk of open-air blasting
CN111784136A (en) A Dynamic Early Warning Method of Impact Risk Based on Analytic Hierarchy Process and Fuzzy Mathematics
CN113570226A (en) Evaluation method of probability grade of water inrush disaster in tunnel in fault fracture zone
CN110610285A (en) A risk classification evaluation method for gobs in underground metal mines
CN106296475B (en) Multi-type evidence synthesis assessment method for tunnel and underground engineering water inrush disaster
CN110458367B (en) A method for evaluating the geological hazard of coal mine working face gas
CN112085921A (en) Landslide comprehensive monitoring and early warning method based on displacement and power multi-parameter
CN106407545A (en) Soft soil foundation ditch safety evaluation method based on fuzzy comprehensive judgment method
WO2014166244A1 (en) Method for evaluating water inrush weakness of coal seam floor
CN104133985A (en) Evaluation method of geological disaster risk attribute interval of tunnel and underground construction
CN110378551A (en) A kind of vcehicular tunnel facility military service method of evaluating performance based on big data
CN109359374A (en) Two-level fuzzy comprehensive evaluation method for water inrush evaluation of coal seam floor
Wang et al. Risk Assessment of Water Inrush in Karst Tunnels Based on the Efficacy Coefficient Method.
CN112100851A (en) Method for evaluating tunnel water inrush disaster risk based on set pair analysis
CN115994398A (en) Method for evaluating collapse risk of deep-buried granite tunnel
CN105046080A (en) Rock mass quality evaluation method
CN115660487A (en) A Comprehensive Quality Evaluation Method for Tunnel Smooth Blasting
CN113216946B (en) A detection method for borehole wall caving and falling blocks based on cuttings and element logging data
CN115146971A (en) Rock burst mine distinguishing and grading method based on equivalent depth
CN115526389A (en) A geological risk prediction method for coal mine dynamic disasters based on sedimentary facies
CN106840843A (en) A kind of Rock burst proneness based on bayesian theory differentiates and method for early warning
CN107784437B (en) Stress concentration-based coal roadway driving face outburst danger identification method
CN106991521A (en) A kind of tunnels and underground engineering slip casting effect grade evaluation method and system

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: 20201016