CN106703883B - A kind of personalized method for determining Water Inrush From Working-faces danger classes - Google Patents
A kind of personalized method for determining Water Inrush From Working-faces danger classes Download PDFInfo
- Publication number
- CN106703883B CN106703883B CN201611243211.3A CN201611243211A CN106703883B CN 106703883 B CN106703883 B CN 106703883B CN 201611243211 A CN201611243211 A CN 201611243211A CN 106703883 B CN106703883 B CN 106703883B
- Authority
- CN
- China
- Prior art keywords
- water inrush
- floor
- water
- risk level
- inrush
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; 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
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Mining & Mineral Resources (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Economics (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Geology (AREA)
- Educational Administration (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明公开了一种个性化确定采煤工作面底板突水危险等级的方法,属于煤炭安全开采领域。该方法综合考虑采煤工作面底板水压、隔水层厚度、煤层底板矿压破坏现象、隔水段岩层岩性组合、底板承压地下水原始导升发育现象等因素,建立这些影响因素与突水危险等级的隐式函数依赖关系,通过插值分析该隐式函数的突变点,建立评价底板突水危险等级的综合评价指标体系。该方法全面考虑影响煤炭开采底板突水威胁的各种因素,建立各影响因素与突水系数的函数依赖关系,解决了传统方法只考虑水压力与隔水层厚度、或无法建立影响因素与突水系数函数依赖关系的难题,为突水系数计算提供了新方法。
The invention discloses a method for individually determining the danger level of water inrush from the floor of a coal mining face, belonging to the field of coal safety mining. This method comprehensively considers the factors such as the water pressure of the coal mining face floor, the thickness of the water-resisting layer, the failure phenomenon of the coal seam floor pressure, the lithology combination of the rock strata in the water-resisting section, and the development of the original conduction of groundwater under pressure on the floor, and establishes the relationship between these influencing factors and the outburst. The implicit function dependence of the water inrush hazard level is analyzed through interpolation to establish a comprehensive evaluation index system for evaluating the water inrush hazard level of the floor. This method comprehensively considers various factors that affect the threat of water inrush from the coal mining floor, and establishes the functional dependence relationship between each influencing factor and the water inrush coefficient. The problem of functional dependence of water coefficient provides a new method for calculating water inrush coefficient.
Description
技术领域technical field
本发明属于煤矿安全开采技术领域,具体涉及一种个性化确定采煤工作面底板突水危险等级的方法。The invention belongs to the technical field of coal mine safety mining, and in particular relates to a method for individually determining the risk level of water inrush from the floor of a coal mining working face.
背景技术Background technique
目前评价煤矿安全开采底板突水危险性的方法主要采用突水系数法。突水系数计算公式T=p/M是我国1964年焦作水文地质大会战期间借鉴匈牙利工程师韦格弗伦斯对隔水厚度(T=M/p)概念提出的,公式中p为煤层底板含水层水压,M为煤层底板隔水层厚度。该公式提出后,以焦作、峰峰、井陉、邯郸、肥城、淄博等大水矿区底板突水案例以及基础数据为依据,计算得出T≤0.06MPa/m为安全,0.06MPa/m<T≤MPa/m为较安全,0.1MPa/m<T为危险。突水系数反映了水文地质学中地下水渗流最基本的规律。At present, the water inrush coefficient method is mainly used to evaluate the risk of floor water inrush in coal mine safety mining. The formula for calculating the water inrush coefficient T=p/M was proposed by the Hungarian engineer Wegflens on the concept of water-repelling thickness (T=M/p) during the 1964 Jiaozuo Hydrogeological Conference in my country. In the formula, p is the water content of the coal seam floor layer hydraulic pressure, and M is the thickness of the coal seam floor water-resisting layer. After the formula was put forward, based on the floor water inrush cases and basic data in Jiaozuo, Fengfeng, Jingxing, Handan, Feicheng, Zibo and other Dashui mining areas, it was calculated that T≤0.06MPa/m is safe, and 0.06MPa/m<T ≤MPa/m is relatively safe, and 0.1MPa/m<T is dangerous. Water inrush coefficient reflects the most basic law of groundwater seepage in hydrogeology.
在以后相当长的历史时期,突水系数对于指导煤矿安全开采起到了重要的作用。但是随着浅层煤层资源开采殆尽,煤炭开采逐步转向深部,奥灰水的水压逐步升高,在实际开采中突水系数远远大于《煤矿防治水规定》的0.1MPa/m。后期虽然有人考虑煤层底板破坏程度、隔水段岩层岩性组合及其它影响因素,提出了各种突水系数的计算公式,但所提出的邻接系数仍然采用1964年提出的邻接值。显然,这是不妥当的。In a rather long historical period, the water inrush coefficient played an important role in guiding the safe mining of coal mines. However, as the shallow coal seam resources are exhausted, coal mining gradually shifts to the deep, and the hydraulic pressure of the Austrian ash water gradually increases. In actual mining, the water inrush coefficient is far greater than 0.1MPa/m in the "Regulations on Water Prevention and Control in Coal Mine". In the later period, although some people put forward various calculation formulas for the water inrush coefficient considering the damage degree of the coal seam floor, the lithological combination of the rock formation in the water-resisting section, and other influencing factors, the adjacency coefficient proposed in 1964 still adopts the adjacency value proposed in 1964. Obviously, this is inappropriate.
目前的突水系数计算方法存在如下问题:The current calculation method of water inrush coefficient has the following problems:
1.影响突水危险性的因素考虑不全面1. Factors affecting the risk of water inrush are not fully considered
突水系数计算公式仅考虑了水压与隔水层厚度,没有考虑煤层底板矿压破坏程度、隔水段岩层岩性组合、底板承压地下水导升高度、含水层富水性指标等有关因素。The calculation formula of water inrush coefficient only considers the water pressure and the thickness of the water-resisting layer, but does not take into account the damage degree of the coal seam floor, the lithology combination of the rock formation in the water-resisting section, the groundwater conduction height of the floor under pressure, and the water-rich index of the aquifer.
2.没有采用机器学习方法建立突水危险性与影响因素间的函数依赖关系2. No machine learning method was used to establish the functional dependence between the risk of water inrush and its influencing factors
早期的突水系数计算公式考虑因素少,计算公式简单。后期的计算公式考虑因素多,在采集样本较少的情况下,无法拟合出突水系数与影响因素间的函数依赖关系。The calculation formula of water inrush coefficient in the early stage has few factors to consider, and the calculation formula is simple. There are many factors to be considered in the later calculation formula, and the functional dependence between the water inrush coefficient and the influencing factors cannot be fitted in the case of a small number of samples.
3.突水系数给出的仅仅是突不突水的危险性,而与突水量无关3. The water inrush coefficient only gives the risk of water inrush, but has nothing to do with the amount of water inrush
从其计算公式T=p/M可以看出,目前的突水系数仅考虑煤层底板含水层水压,煤层底板隔水层厚度等因素,与实际突水量没有关系。在实际开采过程中如果底板含水层富水性较差,即使是突水,因突水量较小,也不会对煤炭开采工作面带来较大危害。It can be seen from its calculation formula T=p/M that the current water inrush coefficient only considers the water pressure of the coal seam floor aquifer, the thickness of the coal seam floor water-resisting layer and other factors, and has nothing to do with the actual water inrush. In the actual mining process, if the water-richness of the floor aquifer is poor, even if the water inrush is small, it will not bring great harm to the coal mining face.
发明内容Contents of the invention
针对现有技术存在的上述问题,本发明提出了一种个性化确定采煤工作面底板突水危险等级的方法。Aiming at the above-mentioned problems existing in the prior art, the present invention proposes a method for individually determining the risk level of water inrush from the floor of the coal mining face.
本发明所采用的技术解决方案是:The technical solution adopted in the present invention is:
一种个性化确定采煤工作面底板突水危险等级的方法,按以下步骤进行:A method for individually determining the risk level of water inrush from the floor of a coal mining face is carried out according to the following steps:
步骤一:收集煤矿生产过程中采煤工作面综合数据,分析最大突水量与多种影响因素的相关性,选择相应种类影响因素作为影响底板突水危险等级的自变量;Step 1: Collect the comprehensive data of the coal mining face in the coal mine production process, analyze the correlation between the maximum water inrush volume and various influencing factors, and select the corresponding type of influencing factors as the independent variable affecting the risk level of the floor water inrush;
步骤二:根据具体矿山综合排水能力与最大突水量数据,个性化划分底板突水危险等级;Step 2: According to the comprehensive drainage capacity and maximum water inrush data of specific mines, individualize the risk level of floor water inrush;
步骤三:将自变量的向量组合记作x,底板突水危险等级记作y,建立自变量与底板突水危险等级的隐式函数依赖关系;Step 3: Denote the vector combination of independent variables as x and the risk level of water inrush from the floor as y, and establish the implicit functional dependence between the independent variable and the risk level of water inrush from the floor;
步骤四:将新采集到的自变量的向量组合x代入步骤三中建立的隐式函数依赖关系,计算得到底板突水危险等级y,然后依据步骤二确定底板突水危险等级。Step 4: Substituting the vector combination x of the newly collected independent variables into the implicit functional dependency established in step 3 to calculate the floor water inrush risk level y, and then determine the floor water inrush risk level according to step 2.
上述方法还包括以下步骤:The above method also includes the following steps:
步骤五:利用步骤三建立的隐式函数依赖关系,对各自变量细分插值,在各底板突水危险等级范围内,反向演算求取各自变量的极大、极小值,建立各自变量与底板突水危险等级的综合评价指标体系;Step 5: Use the implicit functional dependence relationship established in Step 3 to subdivide and interpolate the respective variables. Within the scope of the water inrush hazard level of each floor, perform reverse calculation to obtain the maximum and minimum values of the respective variables, and establish the relationship between the respective variables and A comprehensive evaluation index system for the risk level of floor water inrush;
步骤六:根据步骤五建立的各自变量与底板突水危险等级的综合评价指标体系,现场工作人员对新采集到的自变量数据,即可通过查询综合评价指标体系确定底板突水危险等级。Step 6: According to the comprehensive evaluation index system of the respective variables and the risk level of floor water inrush established in step 5, the on-site staff can determine the risk level of floor water inrush by querying the comprehensive evaluation index system for the newly collected independent variable data.
优选的,步骤一中:所述采煤工作面综合数据包括最大突水量、突水系数、煤层底板矿压破坏程度、隔水段岩层岩性组合、底板承压地下水导升高度及含水层富水性指标等;分析最大突水量分别与突水系数、煤层底板矿压破坏程度、隔水段岩层岩性组合、底板承压地下水导升高度及含水层富水性指标等多种影响因素的相关性,并根据相关性大小选择相应种类影响因素如煤层底板矿压破坏程度、底板承压地下水导升高度等作为影响底板突水危险等级的自变量。Preferably, in step 1: the comprehensive data of the coal mining face include the maximum water inrush volume, water inrush coefficient, damage degree of coal seam floor pressure, rock formation lithology combination of water-resisting section, groundwater conduction height of floor under pressure and aquifer richness. Water index, etc.; analyze the correlation between the maximum water inrush volume and various influencing factors such as water inrush coefficient, coal seam floor pressure damage degree, rock formation lithology combination of water-resisting section, groundwater conduction height of floor under pressure, and water-rich index of aquifer , and according to the size of the correlation, the corresponding types of influencing factors such as the damage degree of the coal seam floor pressure, the height of groundwater under pressure on the floor, etc. are selected as independent variables that affect the risk level of floor water inrush.
优选的,步骤二中:底板突水危险等级划分为安全、较危险、危险三个等级;最大突水量小于或等于具体矿山综合排水能力三分之二定为安全等级1;最大突水量大于具体矿山综合排水能力三分之二且小于或等于具体矿山综合排水能力定为较危险等级2;最大突水量大于具体矿山综合排水能力定为危险等级3。Preferably, in step 2: the floor water inrush hazard level is divided into three levels: safe, relatively dangerous, and dangerous; the maximum water inrush volume is less than or equal to two-thirds of the comprehensive drainage capacity of the specific mine, which is defined as safety level 1; the maximum water inrush volume is greater than the specific mine. Two-thirds of the mine's comprehensive drainage capacity and less than or equal to the specific mine's comprehensive drainage capacity are defined as relatively dangerous level 2; the maximum water inrush is greater than the specific mine's comprehensive drainage capacity is defined as dangerous level 3.
优选的,步骤三中:根据具体矿山已知n个观测样本(x1,y1),(x2,y2)……(xn,yn)在若干函数{f(x,ω)}中求一个最优函数f(x,ω0),对未知依赖关系进行估计,使式R(ω)=∫L(y,f(x,ω))dF(x,y)所示的期望风险最小。Preferably, in Step 3: According to the specific mine, n observation samples (x 1 , y 1 ), (x 2 , y 2 )...(x n , y n ) are known in several functions {f(x,ω) } to find an optimal function f(x,ω 0 ), and estimate the unknown dependencies, so that the formula R(ω)=∫L(y,f(x,ω))dF(x,y) shows Expect minimal risk.
上述步骤三中,建立的自变量与底板突水危险等级的隐式函数依赖关系,借助机器学习理论实现,具体地:In the third step above, the implicit functional dependence between the established independent variable and the risk level of water inrush on the floor is realized with the help of machine learning theory, specifically:
(1){f(x,ω)}为预测函数集,ω称为广义参数,L(y,f(x,ω))为损失函数;(1) {f(x,ω)} is the prediction function set, ω is called the generalized parameter, and L(y,f(x,ω)) is the loss function;
(2)损失函数采用L(y,f(x,ω))=(y-f(x,ω))2;(2) The loss function adopts L(y,f(x,ω))=(yf(x,ω)) 2 ;
(3)在训练过程中可以利用的信息只有样本数据,因此期望风险(3) The information that can be used in the training process is only sample data, so the expected risk
R(ω)=∫L(y,f(x,ω))dF(x,y)无法计算,采用对其进行估算;R(ω)=∫L(y,f(x,ω))dF(x,y) cannot be calculated, use estimate it;
(4)在计算过程中,训练采用回归预测的危险等级是连续值,采用四舍五入的方法取整。(4) In the calculation process, the risk level predicted by regression training is a continuous value, which is rounded to an integer.
优选的,步骤四中:自变量的向量组合x在施工勘探、采煤工作过程中实际获取;对于未进行物理勘探的地段,通过三维矿山数字模型进行插值预测,先预测、后调整。Preferably, in step 4: the vector combination x of the independent variables is actually obtained during the construction exploration and coal mining; for the section without physical exploration, the interpolation prediction is performed through the three-dimensional mine digital model, and the prediction is performed first and then adjusted.
优选的,步骤五中具体地:Preferably, in step five, specifically:
(1)分析自变量的向量组合x的各个分量的极大、极小值;(1) Analyze the maximum and minimum values of each component of the vector combination x of the independent variable;
(2)以极大值减去极小值之差的百分之一为单位,细分插值;(2) Take one percent of the difference between the maximum value minus the minimum value as the unit, and subdivide and interpolate;
(3)求取各插值点的底板突水危险等级,此时的危险等级不进行四舍五入处理;(3) Obtain the water inrush hazard level of the bottom plate at each interpolation point, and the hazard level at this time will not be rounded;
(4)求突水危险等级在区间[0,1.5)、[1.5,2.5)、[2.5,3.5)时,各自变量的极大、极小值,建立各自变量与底板突水危险等级的综合评价指标体系。(4) When the water inrush risk level is in the interval [0, 1.5), [1.5, 2.5), [2.5, 3.5), the maximum and minimum values of the respective variables are established, and the synthesis of the respective variables and the floor water inrush risk level is established Evaluation System.
本发明的有益技术效果是:The beneficial technical effect of the present invention is:
与现有技术相比,本发明综合考虑具体矿山排水能力、影响突水量的各种因素,建立突水危险评级与影响因素间的隐式函数依赖关系、采用参数区间插值确定指标框架,解决了突水数据小样本建模、隐式函数应用困难等技术难题,使得突水危险等级评价更加科学合理,符合矿山生产实际。Compared with the prior art, the present invention comprehensively considers the specific mine drainage capacity and various factors affecting the water inrush volume, establishes the implicit functional dependence relationship between the water inrush risk rating and the influencing factors, and adopts parameter interval interpolation to determine the index framework, which solves the problem of Technical problems such as small sample modeling of water inrush data and difficulties in the application of implicit functions make the evaluation of water inrush hazard levels more scientific and reasonable, which is in line with the actual mine production.
附图说明Description of drawings
图1为本发明一种个性化确定采煤工作面底板突水危险等级方法的流程图。Fig. 1 is a flow chart of a method for individually determining the risk level of water inrush from the floor of a coal mining face according to the present invention.
具体实施方式Detailed ways
本发明提供了一种个性化确定采煤工作面底板突水危险等级的方法。该方法综合考虑采煤工作面底板水压、隔水层厚度、煤层底板矿压破坏现象、隔水段岩层岩性组合、底板承压地下水原始导升发育现象等因素,建立这些影响因素与突水危险等级的隐式函数依赖关系,通过插值分析该隐式函数的突变点,建立评价底板突水危险等级的综合评价指标体系。该方法全面考虑影响煤炭开采底板突水威胁的各种因素,建立各影响因素与突水系数的函数依赖关系,解决了传统方法只考虑水压力与隔水层厚度或无法建立影响因素与突水系数函数依赖关系的难题,为突水系数计算提供了新方法。The invention provides a method for individually determining the risk level of water inrush from the floor of a coal mining face. This method comprehensively considers the factors such as the water pressure of the coal mining face floor, the thickness of the water-resisting layer, the damage phenomenon of the coal seam floor pressure, the lithology combination of the rock strata in the water-resisting section, and the development of the original conduction of groundwater under pressure on the floor, and establishes the relationship between these influencing factors and the outburst. The implicit function dependence of the water inrush hazard level is analyzed through interpolation to establish a comprehensive evaluation index system for evaluating the water inrush hazard level of the floor. This method fully considers various factors that affect the threat of water inrush from the coal mining floor, and establishes the functional dependence relationship between each influencing factor and the water inrush coefficient, which solves the problem that the traditional method only considers the water pressure and the thickness of the water-resisting layer or cannot establish the relationship between the influencing factors and the water inrush. The difficult problem of coefficient function dependence provides a new method for the calculation of water inrush coefficient.
下面结合附图以及具体实施方式对本发明作进一步详细说明:Below in conjunction with accompanying drawing and specific embodiment the present invention is described in further detail:
如图1所示,一种用于个性化确定采煤工作面底板突水危险等级的方法,按照如下步骤进行:As shown in Figure 1, a method for individually determining the risk level of water inrush from the floor of a coal mining face follows the steps below:
步骤一:收集煤矿生产过程中采煤工作面最大突水量、突水系数、煤层底板矿压破坏程度、隔水段岩层岩性组合、底板承压地下水导升高度、含水层富水性指标等综合数据,分析最大突水量与突水系数、煤层底板矿压破坏程度等其它影响因素的相关性,选择相关性较高的影响因素作为影响最大突水量的自变量。Step 1: Collect the maximum water inrush volume of the coal mining face during the coal mine production process, the water inrush coefficient, the damage degree of the coal seam floor pressure, the lithology combination of the rock strata in the water-resisting section, the groundwater conduction height of the floor under pressure, and the water-rich index of the aquifer, etc. According to the data, the correlation between the maximum water inrush and other influencing factors such as the water inrush coefficient and the damage degree of coal seam floor pressure is analyzed, and the factors with higher correlation are selected as the independent variables affecting the maximum water inrush.
步骤二:根据具体矿山综合排水能力与最大突水量数据,个性化划分各组数据的底板突水危险等级。一般情况下划分为安全、较危险、危险三个等级。最大突水量小于或等于具体矿山综合排水能力三分之二定为安全等级1;最大突水量大于具体矿山综合排水能力三分之二且小于具体矿山综合排水能力定为较危险等级2;最大突水量大于或等于具体矿山综合排水能力定为危险等级3。Step 2: According to the comprehensive drainage capacity and maximum water inrush data of specific mines, the floor water inrush risk level of each group of data is individually divided. Generally, it is divided into three levels: safe, relatively dangerous, and dangerous. The maximum water inrush is less than or equal to two-thirds of the comprehensive drainage capacity of the specific mine as safety level 1; If the water volume is greater than or equal to the comprehensive drainage capacity of the specific mine, it is classified as hazard level 3.
步骤三:将自变量的向量组合记作x,底板突水危险等级记作y。根据具体矿山已知n个观测样本(x1,y1),(x2,y2)……(xn,yn)在若干函数{f(x,ω)}中求一个最优函数f(x,ω0),对未知依赖关系进行估计,使式R(ω)=∫L(y,f(x,ω))dF(x,y)所示的期望风险最小。Step 3: Denote the vector combination of independent variables as x, and denote the water inrush hazard level of the floor as y. According to the specific mine known n observation samples (x 1 ,y 1 ),(x 2 ,y 2 )...(x n ,y n ) find an optimal function among several functions {f(x,ω)} f(x,ω 0 ), estimate the unknown dependencies, and minimize the expected risk shown by the formula R(ω)=∫L(y,f(x,ω))dF(x,y).
步骤三建立的各影响因素与底板突水危险等级的函数依赖关系为隐式的函数依赖关系,可以借助支持向量机理论等机器学习理论实现。The functional dependence between each influencing factor established in step 3 and the risk level of floor water inrush is an implicit functional dependence, which can be realized with the help of machine learning theories such as support vector machine theory.
具体地:specifically:
(1){f(x,ω)}为预测函数集,ω称为广义参数,L(y,f(x,ω))为损失函数。(1) {f(x,ω)} is the prediction function set, ω is called the generalized parameter, and L(y,f(x,ω)) is the loss function.
(2)损失函数采用L(y,f(x,ω))=(y-f(x,ω))2。(2) The loss function adopts L(y,f(x,ω))=(yf(x,ω)) 2 .
(3)在训练过程中可以利用的信息只有样本数据,因此期望风险R(ω)=∫L(y,f(x,ω))dF(x,y)无法计算,采用对其进行估算。(3) The information that can be used in the training process is only sample data, so the expected risk R(ω)=∫L(y,f(x,ω))dF(x,y) cannot be calculated, using Estimate it.
(4)在计算过程中,训练采用回归预测的危险等级是连续值,采用四舍五入的方法取整。(4) In the calculation process, the risk level predicted by regression training is a continuous value, which is rounded to an integer.
步骤四:根据步骤三建立的隐式函数依赖关系,对新采集到影响因素的数据向量即自变量的向量组合x,即可代入隐式函数计算其底板突水危险等级。底板突水危险等级影响因素可以在施工勘探、采煤工作过程中实际获取。对于未进行物理勘探的地段,可通过三维矿山数字模型进行插值预测,先预测、后调整。Step 4: According to the implicit function dependency relationship established in Step 3, the data vector of the newly collected influencing factors, that is, the vector combination x of independent variables, can be substituted into the implicit function to calculate its floor water inrush risk level. The factors affecting the risk level of floor water inrush can be actually obtained in the process of construction exploration and coal mining. For the sections that have not been physically explored, interpolation prediction can be performed through the 3D mine digital model, and the prediction is made first and then adjusted.
步骤五:为方便现场实际应用,利用步骤三建立的隐式函数依赖关系,对各自变量细分插值,在各底板突水危险等级范围内,反向演算求取各自变量的极大、极小值,建立各自变量与底板突水危险等级的综合评价指标体系。Step 5: In order to facilitate the actual application on site, use the implicit function dependence relationship established in Step 3 to subdivide and interpolate the respective variables, and perform reverse calculation to obtain the maximum and minimum values of the respective variables within the range of the water inrush hazard level of each floor value, and establish a comprehensive evaluation index system for the respective variables and the risk level of floor water inrush.
具体地:specifically:
(1)分析底板突水危险等级影响因素(自变量的向量组合x)的各个分量的极大、极小值。(1) Analyze the maximum and minimum values of each component of the factors affecting the risk level of floor water inrush (the vector combination x of independent variables).
(2)以极大值减去极小值之差的百分之一为单位,细分插值。(2) Subdividing and interpolating in units of one percent of the difference between the maximum value minus the minimum value.
(3)求取各插值点的底板突水危险等级,此时的危险等级不进行四舍五入处理。(3) Obtain the water inrush hazard level of the bottom plate at each interpolation point, and the hazard level at this time will not be rounded.
(4)求突水危险等级在区间[0,1.5)、[1.5,2.5)、[2.5,3.5)时,各自变量的极大、极小值,建立各自变量与底板突水危险等级的综合评价指标体系。(4) When the water inrush risk level is in the interval [0, 1.5), [1.5, 2.5), [2.5, 3.5), the maximum and minimum values of the respective variables are established, and the synthesis of the respective variables and the floor water inrush risk level is established Evaluation System.
步骤六:根据步骤5建立的各自变量与底板突水危险等级的综合评价指标体系,现场工作人员对新采集到影响因素的数据,即可通过查询综合评价指标体系确定底板突水危险等级。Step 6: According to the comprehensive evaluation index system of the respective variables and the risk level of floor water inrush established in step 5, the on-site staff can determine the risk level of floor water inrush by querying the comprehensive evaluation index system for the newly collected data of influencing factors.
本发明在步骤一中提出突水危险性的关键在于突水量,而不仅仅是底板水压力与隔水层厚度的比值。并在步骤一中综合考虑了影响突水量的突水系数、煤层底板矿压破坏程度、隔水段岩层岩性组合、底板承压地下水导升高度、含水层富水性指标等因素。The present invention proposes that the key to the risk of water inrush in step one is the water inrush volume, not just the ratio of the water pressure on the bottom plate to the thickness of the water-resisting layer. And in the first step, factors such as the water inrush coefficient affecting the water inrush volume, the damage degree of the coal seam floor pressure, the lithology combination of the rock formation in the water-resisting section, the groundwater conduction height of the floor under pressure, and the water-rich index of the aquifer are considered comprehensively.
在步骤二中根据具体矿山排水能力个性化划分矿山采煤工作面突水危险等级,而不是统一划分所有矿山突水危险等级。In the second step, according to the specific mine drainage capacity, the water inrush risk level of the coal mining face is divided individually, rather than all mine water inrush risk levels are uniformly divided.
在无法给出突水影响因素与最大突水量显式函数依赖关系的条件下,在步骤三中采用隐式函数依赖关系描述最大突水量与各影响因素的函数依赖关系,该函数依赖关系要求:对于数据样本(x1,y1),(x2,y2),(x3,y3)……(xn,yn),在若干函数{f(x,ω)}中求一个最优函数f(x,ω0),对未知依赖关系进行估计,使得下式所示的期望风险最小。Under the condition that the explicit functional dependence relationship between the influencing factors of water inrush and the maximum water inrush volume cannot be given, the implicit functional dependence relationship is used to describe the functional dependence relationship between the maximum water inrush volume and each influencing factor in step 3. The functional dependence relationship requires: For data samples (x 1 ,y 1 ),(x 2 ,y 2 ),(x 3 ,y 3 )...(x n ,y n ), find one of several functions {f(x,ω)} The optimal function f(x,ω 0 ) estimates the unknown dependencies so that the expected risk shown in the following formula is the smallest.
R(ω)=∫L(y,f(x,ω))dF(x,y)R(ω)=∫L(y,f(x,ω))dF(x,y)
在步骤五中对各影响因素细分插值,在各底板突水危险等级范围内,反向演算求取各自变量的极大、极小值,建立各自变量与底板突水危险等级的综合评价指标体系。In step 5, subdivide and interpolate each influencing factor, within the scope of each floor water inrush risk level, reverse calculation to obtain the maximum and minimum values of the respective variables, and establish a comprehensive evaluation index for each variable and the floor water inrush risk level system.
本发明提出了一种用于个性化确定采煤工作面底板突水危险等级的方法,与现有技术相比,本发明综合考虑具体矿山排水能力、影响突水量的各种因素,建立突水危险评级与影响因素间的隐式函数依赖关系、采用参数区间插值确定指标框架,解决了突水数据小样本建模、隐式函数应用困难等技术难题,使得突水危险等级评价更加科学合理,符合矿山生产实际。The present invention proposes a method for individually determining the risk level of water inrush on the floor of a coal mining face. Compared with the prior art, the present invention comprehensively considers the specific mine drainage capacity and various factors that affect the amount of water inrush, and establishes the water inrush risk level. The implicit functional dependency between risk rating and influencing factors, and the use of parameter interval interpolation to determine the index framework solve technical problems such as small sample modeling of water inrush data and difficulties in the application of implicit functions, making the evaluation of water inrush risk levels more scientific and reasonable. In line with the actual production of mines.
上述方式中未述及的部分采取或借鉴已有技术即可实现。Parts not mentioned in the above methods can be realized by adopting or referring to existing technologies.
上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。The above descriptions are 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 scope of the present invention shall also belong to the scope of the present invention. protected range.
Claims (7)
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611243211.3A CN106703883B (en) | 2016-12-29 | 2016-12-29 | A kind of personalized method for determining Water Inrush From Working-faces danger classes |
| PCT/CN2017/108619 WO2018121035A1 (en) | 2016-12-29 | 2017-10-31 | Customized method for determining coal mining face floor water inrush risk level |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611243211.3A CN106703883B (en) | 2016-12-29 | 2016-12-29 | A kind of personalized method for determining Water Inrush From Working-faces danger classes |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106703883A CN106703883A (en) | 2017-05-24 |
| CN106703883B true CN106703883B (en) | 2018-03-13 |
Family
ID=58906146
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611243211.3A Active CN106703883B (en) | 2016-12-29 | 2016-12-29 | A kind of personalized method for determining Water Inrush From Working-faces danger classes |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN106703883B (en) |
| WO (1) | WO2018121035A1 (en) |
Families Citing this family (48)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106703883B (en) * | 2016-12-29 | 2018-03-13 | 山东科技大学 | A kind of personalized method for determining Water Inrush From Working-faces danger classes |
| CN107165626A (en) * | 2017-06-30 | 2017-09-15 | 徐州市耐力高分子科技有限公司 | A kind of coal-face floods prevention method with cranny development country rock top plate |
| CN107237644B (en) * | 2017-06-30 | 2018-11-06 | 湖南科技大学 | Tunnel inverted arch three-dimensional gushing water destroys the determination method of critical hydraulic pressure and critical thickness |
| CN107784437B (en) * | 2017-10-16 | 2021-09-28 | 重庆大学 | Stress concentration-based coal roadway driving face outburst danger identification method |
| CN109359374A (en) * | 2018-10-10 | 2019-02-19 | 山东科技大学 | Two-level fuzzy comprehensive evaluation method for water inrush evaluation of coal seam floor |
| CN109800955A (en) * | 2018-12-24 | 2019-05-24 | 永城煤电控股集团有限公司 | Coal seam bottom water bursting hazard assessment calculation method |
| CN109948268B (en) * | 2019-01-21 | 2023-10-17 | 安迈智能(北京)矿山科技股份有限公司 | An automatic identification method for comprehensive risks on the working face |
| CN109598102B (en) * | 2019-02-01 | 2022-05-27 | 内蒙古科技大学 | Method, device, equipment and medium for predicting the danger degree of rock burst in coal mine |
| CN110532872A (en) * | 2019-07-24 | 2019-12-03 | 宁德市公路局 | A kind of landslide hierarchy system and method based on convolution supporting vector neural network |
| CN110533224B (en) * | 2019-08-06 | 2020-08-28 | 山东科技大学 | Oil shale continuous exploration drilling position optimization method |
| CN110552741B (en) * | 2019-09-09 | 2021-01-19 | 中煤科工集团西安研究院有限公司 | Coal face bottom plate water inrush comprehensive monitoring and early warning system and method |
| CN111127234B (en) * | 2019-10-11 | 2024-01-19 | 重庆大学 | A method and device for determining the first mining layer in outburst coal seam group mining |
| CN111260216B (en) * | 2020-01-15 | 2023-06-30 | 山东大学 | A comprehensive evaluation and prevention method for seepage field of underground water-sealed oil storage depot during operation period |
| CN111222254A (en) * | 2020-01-22 | 2020-06-02 | 西安科技大学 | Method and system for classifying rockburst hazard level of working face based on stress superposition method |
| CN111239840B (en) * | 2020-02-25 | 2022-05-24 | 华北科技学院 | Floor water inrush early warning method based on high-density electrical method |
| CN111563653A (en) * | 2020-04-03 | 2020-08-21 | 山东大学 | Early warning construction method for water-rich broken stratum of underground engineering |
| CN111652490A (en) * | 2020-05-28 | 2020-09-11 | 山东科技大学 | A new method for prediction of structural water inrush in deep mines based on fracture mechanics theory |
| CN111652509B (en) * | 2020-06-03 | 2024-02-13 | 华北科技学院 | Method for classifying and distinguishing water inrush risk of Taiyuan limestone aquifer based on multiple variables |
| CN111562285A (en) * | 2020-06-03 | 2020-08-21 | 安徽大学 | Mine water inrush source identification method and system based on big data and deep learning |
| CN111814322B (en) * | 2020-06-24 | 2023-11-10 | 应急管理部信息研究院 | Working face overlying strata damage height determining method based on half-plane theory |
| CN114329680B (en) * | 2020-10-09 | 2024-04-16 | 神华神东煤炭集团有限责任公司 | A stability evaluation method for underground reservoir pillar dam in mining area and its application |
| CN112668873B (en) * | 2020-12-25 | 2024-03-05 | 中国矿业大学 | Mine safety situation analysis and prediction and early warning methods |
| CN112906280A (en) * | 2021-03-11 | 2021-06-04 | 宁夏安普安全技术咨询有限公司 | Mathematical model establishing method for safety evaluation and risk prediction |
| CN112966949B (en) * | 2021-03-15 | 2022-06-10 | 北京市市政工程研究院 | Tunnel construction risk assessment method and device and storage medium |
| CN113294143B (en) * | 2021-04-16 | 2023-09-26 | 中煤能源研究院有限责任公司 | Method for evaluating investigation treatment effect of advanced ground area of limestone water damage of coal seam floor |
| CN113449414B (en) * | 2021-06-07 | 2023-03-28 | 西安科技大学 | Three-layer structure-based bottom plate slippage failure depth calculation method |
| CN113449415B (en) * | 2021-06-07 | 2023-02-24 | 西安科技大学 | A Calculation Method Based on Sliding Damage Depth of Double-layer Structure Floor |
| CN113516414A (en) * | 2021-08-09 | 2021-10-19 | 江苏徐矿能源股份有限公司 | Method for determining rock burst danger level |
| CN113565490B (en) * | 2021-08-31 | 2023-08-08 | 中煤科工集团重庆研究院有限公司 | A method for early warning of microseismic water damage |
| CN113623004B (en) * | 2021-08-31 | 2024-02-13 | 中煤科工集团重庆研究院有限公司 | Judgment method for water damage early warning |
| CN113914928B (en) * | 2021-09-06 | 2024-06-25 | 中煤科工开采研究院有限公司 | Method for dividing support areas of fully mechanized coal mining face support and accurately supporting coal mine fully mechanized coal mining face support |
| CN114087022B (en) * | 2021-10-28 | 2023-11-28 | 山东科技大学 | Coal seam floor variable parameter water inrush channel early warning system and water inrush risk judging method |
| CN114412567B (en) * | 2021-12-08 | 2023-03-14 | 中国矿业大学 | Early warning method for in-situ water-retention coal mining on water with bearing pressure of bottom plate limestone |
| CN114219260A (en) * | 2021-12-08 | 2022-03-22 | 郑州煤矿机械集团股份有限公司 | A method for evaluating coal mining adaptability in thin coal seam working face |
| CN114151133A (en) * | 2021-12-14 | 2022-03-08 | 中煤科工开采研究院有限公司 | Grading early warning method for incoming pressure of top plate of fully mechanized mining face |
| CN114548676A (en) * | 2022-01-18 | 2022-05-27 | 中交第二航务工程局有限公司 | Tunnel granite fault water burst risk level prediction method |
| CN114659898B (en) * | 2022-02-25 | 2025-03-28 | 中国矿业大学 | A method for evaluating the bearing capacity of intact floor in pressure mining |
| CN114626720A (en) * | 2022-03-16 | 2022-06-14 | 华北科技学院(中国煤矿安全技术培训中心) | Method for constructing coal spontaneous combustion risk evaluation model under goaf inerting condition |
| CN115062836A (en) * | 2022-06-10 | 2022-09-16 | 徐州矿务集团有限公司 | Method and system for predicting water inrush pressure rack accident under complex geological condition |
| CN115114476B (en) * | 2022-07-26 | 2022-11-15 | 汶上义桥煤矿有限责任公司 | Image processing-based monitoring video storage method for coal washing transmission equipment |
| CN116797020A (en) * | 2023-05-24 | 2023-09-22 | 中国矿业大学 | A microseismic early warning method for coal mine roof detachment water inrush considering the evolution of rock structure |
| CN117057601B (en) * | 2023-08-02 | 2024-01-30 | 中国安全生产科学研究院 | Non-coal mine safety monitoring and early warning system based on Internet of Things |
| CN118154046B (en) * | 2024-05-10 | 2024-07-23 | 太原向明智控科技有限公司 | Top plate pressure grade dividing method |
| CN118191967B (en) * | 2024-05-14 | 2024-08-06 | 中煤科工西安研究院(集团)有限公司 | Intelligent early warning system and method for full-space three-dimensional monitoring of water damage risk of coal seam roof |
| CN119624103B (en) * | 2024-11-21 | 2025-09-05 | 北京科技大学 | A comprehensive evaluation and prediction method for typical construction risks of deep shafts |
| CN119784149B (en) * | 2024-12-13 | 2025-09-23 | 新汶矿业集团有限责任公司 | Mine water inrush disaster prevention and control evaluation method based on fault water intelligent identification and early warning |
| CN120273778B (en) * | 2025-03-21 | 2025-11-07 | 陕西金源招贤矿业有限公司 | Method for forecasting water burst position of separation layer carrying silt of coal mining working face |
| CN120124334B (en) * | 2025-05-15 | 2025-12-16 | 北京交通大学 | A Dynamic Variable-Weight Risk Assessment and Active Control Method for Water and Sediment Inrush |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1719635A1 (en) * | 1988-12-19 | 1992-03-15 | Н.И.Никифоров, А.М.Оллыкайнен и И.Г.Новиков | Method to protect mine workings against influx of surface water |
| CN101699451A (en) * | 2009-05-08 | 2010-04-28 | 中国矿业大学(北京) | Novel practical method frangibility index method for evaluating seam floor water inrush |
| CN101894189B (en) * | 2010-07-14 | 2011-04-20 | 中国矿业大学(北京) | New method for evaluating coal seam bottom water bursting |
| CN102194056B (en) * | 2011-05-05 | 2012-03-21 | 中国矿业大学(北京) | BN-GIS (Bayesian Network-Geographic Information System) method for evaluating and predicting water inrush danger of coal-seam roof and floor |
| CN103049645B (en) * | 2012-11-28 | 2015-12-02 | 山东科技大学 | A kind of coal seam floor water-inrush risk evaluation method |
| CN103279809B (en) * | 2013-06-09 | 2017-02-08 | 山东科技大学 | Method for predicting and evaluating water-inrush from seam floor based on bidirectional impact of indexes |
| CN104156560A (en) * | 2014-07-12 | 2014-11-19 | 中国矿业大学 | Multi-level coal mine water inrush prediction method based on SaE-ELM (self-adaptive evolutionary extreme learning machine) |
| CN104502995A (en) * | 2014-12-15 | 2015-04-08 | 中国矿业大学 | A Ts-q method for risk assessment of floor water inrush in deep coal seam mining |
| CN104766242A (en) * | 2015-03-25 | 2015-07-08 | 山东科技大学 | Method for evaluating dangerousness of water inrush from coal floor |
| CN105069689B (en) * | 2015-08-21 | 2017-03-29 | 山东科技大学 | Based on the coal seam floor water-inrush risk evaluation method that grey correlation is combined with FDAHP |
| CN106703883B (en) * | 2016-12-29 | 2018-03-13 | 山东科技大学 | A kind of personalized method for determining Water Inrush From Working-faces danger classes |
-
2016
- 2016-12-29 CN CN201611243211.3A patent/CN106703883B/en active Active
-
2017
- 2017-10-31 WO PCT/CN2017/108619 patent/WO2018121035A1/en not_active Ceased
Also Published As
| Publication number | Publication date |
|---|---|
| CN106703883A (en) | 2017-05-24 |
| WO2018121035A1 (en) | 2018-07-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106703883B (en) | A kind of personalized method for determining Water Inrush From Working-faces danger classes | |
| Zhao et al. | Groundwater impact of open cut coal mine and an assessment methodology: A case study in NSW | |
| Batycky et al. | Revisiting reservoir flood-surveillance methods using streamlines | |
| CN111706317B (en) | Method for determining distribution condition of residual oil in hypotonic reservoir in encryption adjustment area | |
| CN104847379B (en) | A method for strengthening the bottom plate of working face by grouting at the top of Austrian ash | |
| CN104695950A (en) | Prediction method for volcanic rock oil reservoir productivity | |
| CN104699954A (en) | Calibration method for oil recovery rate of karst-fissure-cavity carbonatite reservoir | |
| CN116227710B (en) | Prediction method and system for ecological water level variation under coal mining in ecologically fragile areas | |
| CN111695303A (en) | Method for evaluating water filling strength of coal seam roof sandstone aquifer | |
| CN107391778B (en) | An Analytical Calculation Method of Seepage Flow in Circular Tunnel | |
| CN116025331A (en) | A method for identifying and quantifying invalid water | |
| CN109211745B (en) | Method for recovering evolution process of organic matter-rich shale porosity | |
| CN117454050A (en) | Method for analyzing and evaluating grouting reconstruction of coal seam based on fuzzy analytic hierarchy process | |
| CN103926633B (en) | The quantitative forecasting technique of seat earth coal measures water-resisting layer protrusion-dispelling performance | |
| CN110378010A (en) | A kind of Remaining Oil Distribution quantitative analysis method | |
| CN117853256B (en) | Method and system for zoning water-saving coal mining in mining areas based on post-mining groundwater level recovery time | |
| Adams et al. | A physically based model of rebound in South Crofty tin mine, Cornwall | |
| Wilson et al. | Wairau river-Wairau aquifer interaction | |
| Agarwal et al. | A streamline-based method for assisted history matching applied to an Arabian Gulf field | |
| CN106651612B (en) | A method for evaluating the status quo of oilfield development | |
| CN104573399B (en) | The dynamic water grouting amount prediction of a kind of mine water inrush and slip casting effect evaluation methodology | |
| CN115796339A (en) | Method and device for predicting water inflow of mine working face, electronic equipment and storage medium | |
| Tian et al. | Monitoring and predicting CO2 flooding using material balance equations | |
| CN103984042B (en) | Limestone of mid Ordovician paleocrust of weathering water isolating Forecasting Methodology | |
| CN113027399A (en) | Method for obtaining water flooding curve of high-water-content block based on micro-flow simulation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20170524 Assignee: Changzhou Haike Zhidong Technology Co.,Ltd. Assignor: SHANDONG University OF SCIENCE AND TECHNOLOGY Contract record no.: X2025980038503 Denomination of invention: A method for personalizing the determination of water inrush risk level of coal mining face floor Granted publication date: 20180313 License type: Common License Record date: 20251126 |