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 PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 239000011159 matrix material Substances 0.000 claims abstract description 23
- 238000012544 monitoring process Methods 0.000 claims description 27
- 239000011435 rock Substances 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims 1
- 230000002085 persistent effect Effects 0.000 claims 1
- 230000035939 shock Effects 0.000 abstract description 22
- 238000011156 evaluation Methods 0.000 abstract description 9
- 239000003245 coal Substances 0.000 description 12
- 238000005065 mining Methods 0.000 description 11
- 238000006073 displacement reaction Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000002265 prevention Effects 0.000 description 5
- 238000005553 drilling Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000009412 basement excavation Methods 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 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/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
- 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
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- 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)
- 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;步骤六、根据最大隶属度原则,确定预警对象的冲击危险等级。本发明综合考虑了多种冲击地压预警指标,采用预警模型对预警对象进行冲击危险性预警,预警更加科学,预警效率更高,结果更加可靠。
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.
Description
技术领域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):
式(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
通过将各预警指标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).
式(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
(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)
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)
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 |
-
2020
- 2020-06-22 CN CN202010577804.3A patent/CN111784136A/en active Pending
Non-Patent Citations (1)
Title |
---|
张宁博 等: "基于一孔多点式应力与位移监测系统的掘进巷道冲击危险性评价", 煤炭学报, pages 140 - 149 * |
Cited By (4)
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 |