CN118154046B - Top plate pressure grade dividing method - Google Patents
Top plate pressure grade dividing method Download PDFInfo
- Publication number
- CN118154046B CN118154046B CN202410571181.7A CN202410571181A CN118154046B CN 118154046 B CN118154046 B CN 118154046B CN 202410571181 A CN202410571181 A CN 202410571181A CN 118154046 B CN118154046 B CN 118154046B
- Authority
- CN
- China
- Prior art keywords
- pressure
- real
- time
- index
- calculating
- 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
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000005065 mining Methods 0.000 claims abstract description 18
- 239000003245 coal Substances 0.000 claims abstract description 11
- 230000008859 change Effects 0.000 claims abstract description 8
- 238000007781 pre-processing Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 2
- 230000006870 function Effects 0.000 description 16
- 238000004458 analytical method Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 239000008358 core component Substances 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/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)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Agronomy & Crop Science (AREA)
- Primary Health Care (AREA)
- Mining & Mineral Resources (AREA)
- Animal Husbandry (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
The invention provides a top plate pressure grade dividing method, and belongs to the field of coal mine working face pressure prediction; the problems that the working resistance data of the hydraulic supports at different positions are low in utilization rate and the pressure change of the top plate cannot be reflected correctly are solved; the method comprises the following steps: preprocessing the mine pressure historical time sequence data collected by the hydraulic support, and constructing four indexes: the method comprises the steps of constructing a global real-time early warning function, a real-time average risk function, a real-time maximum risk function and a real-time range function, constructing a comprehensive index according to four indexes, determining weight parameters and dividing thresholds of the indexes, calculating the comprehensive index of the roof pressure, and grading according to the pressure thresholds; the method is suitable for dividing the roof pressure of the fully mechanized mining face, has the advantages of comprehensiveness, systematicness, simplicity, high efficiency and the like compared with the single index of the traditional mining pressure division, has guiding significance for preventing, controlling and managing mining pressure disasters, and can ensure safe and efficient mining of the face.
Description
Technical Field
The invention provides a top plate pressure grade dividing method, and belongs to the technical field of coal mine working face pressure prediction.
Background
In the coal exploitation process, due to severe production environment and complex change, various disasters such as working face caving, roof fall, pressing frame and the like are caused by the influence of factors such as complex ground stress, exploitation dynamic pressure and the like, so that the safety exploitation of coal and the safety of personnel and equipment are seriously influenced, and a large amount of economic loss is caused.
In recent years, the degree of automation and the degree of informatization of fully mechanized coal mining working surfaces are gradually improved, and the safety and the production efficiency of the working surfaces are driven to be partially improved. More and more fully-mechanized mining faces are provided with electrohydraulic control hydraulic supports, and the change of working resistance of the hydraulic supports reflects the mechanical characteristics of overlying strata, so that the perception prediction of the mining pressures of the fully-mechanized mining faces based on massive working resistance data acquired by the hydraulic supports is an effective means for realizing early warning and early response of the periodical pressure of the working faces, and has an important role in optimizing the control quality of surrounding rocks and the early warning of accidents of the roof of the working faces.
However, in the face of massive data collected by a hydraulic support of an underground fully-mechanized mining face, the analysis of data by a general coal mine is only a preliminary stage, the deeper analysis and mining of the data are not achieved, the pressure classification of the data is simply and artificially divided into several types, the indexes are single and simple, the instantaneous and global conditions of pressure change cannot be reflected, and a unified system frame is not established.
How the hydraulic support working resistance time sequence data at different positions can be effectively analyzed is a significant research.
Disclosure of Invention
The invention provides a top plate pressure grade dividing method for solving the problems that working resistance data of hydraulic supports at different positions are low in utilization rate and the pressure change of a top plate cannot be reflected correctly.
In order to solve the technical problems, the invention adopts the following technical scheme: a top plate pressure rating method comprising the steps of:
Step one, acquiring mine pressure historical data of each hydraulic support of a fully mechanized mining face of a coal mine;
step two, preprocessing the collected time series data;
step three, determining weight parameters and different top plate pressure grade thresholds in various index formulas;
Step four, considering the column pressure values of n columns in the m x delta t time period before the t moment, and assuming that the column pressure value of the ith column in the j x delta t time period before the t moment is x i,j (t), wherein delta t is a sampling interval, calculating a real-time early warning function of the ith column The expression is as follows:
；
wherein: ， i=1, 2, … …, n, k=1, 2, … …, n, j=1, 2, … …, m, α, β are hyper-parameters;
step five, calculating a first index: global real-time early warning function of n upright posts The expression is as follows:
；
in the above formula: As a piecewise function when When the utility model is in the standing state,When =1When the utility model is not in the standing state,=0；
Step six, calculating a second index: real-time average risk function of ith columnThe expression is as follows:
；
in the above formula: ，，、 C is a super parameter;
Step seven, calculating a third index: real-time maximum risk function of ith column The expression is as follows:
；
in the above formula: x max is the maximum column pressure;
step eight, calculating a fourth index: real-time pole difference function of ith column The expression is as follows:
；
in the above formula: X min is the minimum column pressure;
Step nine, calculating a comprehensive index W (t), wherein the expression is as follows:
；
And step ten, grading the roof pressure according to the calculation result of the comprehensive index and the set grade threshold value.
In the first step, the historical data of the mine pressure is obtained through a pressure sensor, the pressure sensor is arranged in the lower cavity of the hydraulic support, the change of the liquid pressure in the lower cavity of the upright post is monitored in real time, and the historical data of the mine pressure is a plurality of pressure sensor numerical matrixes read by a plurality of hydraulic supports according to time sequence.
The pretreatment in the second step comprises the following steps: and (3) checking the pure randomness and the stationarity of the mine pressure historical data, correcting the sampling interval, processing the data missing value, and processing the data redundancy and abnormal noise value.
The weight distribution calculation of the four indexes is carried out by adopting a comprehensive index method, and the calculation of the historical pressure time series data is carried out by adopting the comprehensive index method to generate a roof pressure division model.
Compared with the prior art, the invention has the following beneficial effects: the top plate pressure grade classification method provided by the invention is calculated based on the comprehensive index method, can quickly, efficiently and systematically obtain the top plate pressure grade classification of the fully-mechanized mining face comprehensively, and has a wide application range. The method has the specific advantages that:
(1) The hydraulic support has certain universality and strong operability, provides effective basis for grading the roof pressure of the fully mechanized mining face, has guiding significance for hydraulic support selection, mine pressure disaster prevention and control and management, and is beneficial to safe and efficient production of coal mines.
(2) The invention comprehensively considers the continuity and the relativity of the pressure data of a plurality of hydraulic supports within a period of time, and integrates the time series data of the pressures of the hydraulic supports, compared with other methods which only consider the single instantaneous pressure data of a single hydraulic support, the analysis is more comprehensive and reasonable.
(3) The invention comprehensively considers statistical factors such as global and local, real-time and time period, average risk and maximum risk, and the like, and defines main indexes such as global real-time early warning function, real-time average risk function, real-time maximum risk function, real-time range function, and the like.
(4) The invention adopts a comprehensive index method for index value balance aiming at a plurality of main indexes, the comprehensive index method is the most representative, and is an evaluation model established according to related theories such as mathematical statistics and the like on the basis of analyzing the top plate pressure data, and the grade of top plate pressure division can be determined according to the model.
Drawings
The invention is further described below with reference to the accompanying drawings:
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
As shown in fig. 1, the invention provides a roof pressure grade dividing method, which adopts a data driving method to establish a roof pressure dividing model of a fully-mechanized mining working face, so that the roof pressure dividing model becomes an index and a frame which can reflect the relation among a hydraulic support, the hydraulic support, surrounding rocks, the mine pressure law and the like.
The top plate pressure grade classification method is realized based on a comprehensive index method, the method can collect the historical data of the pressure sensor of the hydraulic support, perform data preprocessing, determine weight parameters and grade thresholds in various index formulas, and respectively calculate four indexes: the method comprises the following steps of:
Step one, acquiring mining pressure historical data of each hydraulic support of a fully mechanized mining face of a coal mine, wherein the mining pressure historical data are time sequence data; the mine pressure historical data are obtained through pressure sensors, the pressure sensors in the coal mine are core components of the hydraulic support, the pressure sensors are arranged in the lower cavity of the hydraulic support, and the change of the liquid pressure in the lower cavity of the upright post is monitored in real time to know the interaction condition of the hydraulic support and the top plate.
Step two, preprocessing the collected time series data;
The specific pretreatment flow and links comprise: and (3) checking the pure randomness and the stationarity of the mine pressure historical data, correcting the sampling interval, processing the data missing value, removing redundancy and abnormal noise value of the data, and the like.
Step three, determining weight parameters and different top plate pressure grade thresholds in various index formulas;
For example: the rating threshold may be set to [0,0.2,0.4,0.6,0.8,1].
Step four, considering the column pressure values of n columns in the m x delta t time period before the t moment, and assuming that the column pressure value of the ith column in the j x delta t time period before the t moment is x i,j (t), wherein delta t is a sampling interval, calculating a real-time early warning function of the ith columnThe expression is as follows:
；
wherein: ， i=1, 2, … …, n, k=1, 2, … …, n, j=1, 2, … …, m, α, β are hyper-parameters.
Step five, calculating a first index: global real-time early warning function of n upright postsThe expression is as follows:
；
in the above formula: As a piecewise function when When the utility model is in the standing state,When =1When the utility model is not in the standing state,=0。
Step six, calculating a second index: real-time average risk function of ith columnThe expression is as follows:
；
in the above formula: ，，、 And c is a super parameter.
Step seven, calculating a third index: real-time maximum risk function of ith columnThe expression is as follows:
；
in the above formula: X max is the maximum column pressure.
Step eight, calculating a fourth index: real-time pole difference function of ith columnThe expression is as follows:
；
in the above formula: x min is the minimum column pressure.
Step nine, calculating a comprehensive index W (t), wherein the expression is as follows:
。
And step ten, grading the roof pressure according to the calculation result of the comprehensive index and the set grade threshold value.
For example: when the comprehensive index is less than or equal to 0.2, the pressure is not generated; when the comprehensive index is 0.2 and is less than or equal to 0.4, the pressure is in a weak pressure state; when the comprehensive index is less than or equal to 0.4 and less than or equal to 0.6, the pressure is in a medium pressure state; when the comprehensive index is less than or equal to 0.6 and less than or equal to 0.8, the safety state is displayed; and when the comprehensive index is 0.8 and less than or equal to 1.0, the pressure is in a strong pressure state.
The object of the top plate pressure dividing method based on the comprehensive index method adopted by the invention is the pressure sensors of a plurality of hydraulic supports, namely, the data format of the pressure sensors is a plurality of pressure sensor numerical matrixes read according to time sequence. And constructing a comprehensive index method for index value balance calculation, and calculating historical pressure time series data through the comprehensive index method to generate a roof pressure division model. The weight parameters and the level threshold in various index formulas need to be determined in advance.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some or all of the technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit of the invention.
Claims (4)
1. A top plate pressure grade dividing method is characterized in that: the method comprises the following steps:
Step one, acquiring mine pressure historical data of each hydraulic support of a fully mechanized mining face of a coal mine;
step two, preprocessing the collected time series data;
step three, determining weight parameters and different top plate pressure grade thresholds in various index formulas;
Step four, considering the column pressure values of n columns in the m x delta t time period before the t moment, and assuming that the column pressure value of the ith column in the j x delta t time period before the t moment is x i,j (t), wherein delta t is a sampling interval, calculating a real-time early warning function of the ith column The expression is as follows:
；
wherein: ， i=1, 2, … …, n, k=1, 2, … …, n, j=1, 2, … …, m, α, β are hyper-parameters;
step five, calculating a first index: global real-time early warning function of n upright posts The expression is as follows:
；
in the above formula: As a piecewise function when When the utility model is in the standing state,When =1When the utility model is not in the standing state,=0；
Step six, calculating a second index: real-time average risk function of ith columnThe expression is as follows:
；
in the above formula: ，，、 C is a super parameter;
Step seven, calculating a third index: real-time maximum risk function of ith column The expression is as follows:
；
in the above formula: x max is the maximum column pressure;
step eight, calculating a fourth index: real-time pole difference function of ith column The expression is as follows:
；
in the above formula: X min is the minimum column pressure;
Step nine, calculating a comprehensive index W (t), wherein the expression is as follows:
；
And step ten, grading the roof pressure according to the calculation result of the comprehensive index and the set grade threshold value.
2. A roof pressure rating method as set forth in claim 1 wherein: in the first step, the historical data of the mine pressure is obtained through a pressure sensor, the pressure sensor is arranged in the lower cavity of the hydraulic support, the change of the liquid pressure in the lower cavity of the upright post is monitored in real time, and the historical data of the mine pressure is a plurality of pressure sensor numerical matrixes read by a plurality of hydraulic supports according to time sequence.
3. A roof pressure rating method as set forth in claim 2 wherein: the pretreatment in the second step comprises the following steps: and (3) checking the pure randomness and the stationarity of the mine pressure historical data, correcting the sampling interval, processing the data missing value, and processing the data redundancy and abnormal noise value.
4. A roof pressure rating method according to any one of claims 1 to 3, wherein: the weight distribution calculation of the four indexes is carried out by adopting a comprehensive index method, and the calculation of the historical pressure time series data is carried out by adopting the comprehensive index method to generate a roof pressure division model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410571181.7A CN118154046B (en) | 2024-05-10 | 2024-05-10 | Top plate pressure grade dividing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410571181.7A CN118154046B (en) | 2024-05-10 | 2024-05-10 | Top plate pressure grade dividing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118154046A CN118154046A (en) | 2024-06-07 |
| CN118154046B true CN118154046B (en) | 2024-07-23 |
Family
ID=91297004
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410571181.7A Active CN118154046B (en) | 2024-05-10 | 2024-05-10 | Top plate pressure grade dividing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118154046B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104252686A (en) * | 2014-08-15 | 2014-12-31 | 国家电网公司 | Determination method for power grid safety aggregative indicators |
| CN106703883A (en) * | 2016-12-29 | 2017-05-24 | 山东科技大学 | Method for determining floor water inrush danger level of coal mining working faces in personalized manner |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106960108B (en) * | 2017-04-07 | 2020-04-24 | 中国环境科学研究院 | Reservoir upstream incoming water pressure analysis method based on Bayesian network |
| CN107559045B (en) * | 2017-08-24 | 2019-08-13 | 中国矿业大学 | Underground coal mine Face Pressure Early-Warning System and method |
| CN111256754B (en) * | 2020-01-19 | 2021-08-10 | 河海大学 | Concrete dam long-term operation safety early warning method |
| CN112145231A (en) * | 2020-08-20 | 2020-12-29 | 中煤科工开采研究院有限公司 | Early warning method for working condition of hydraulic support |
| CN113673119B (en) * | 2021-09-07 | 2023-08-18 | 重庆大学 | Coal mine rock burst danger dynamic and static coupling evaluation method based on Bayesian method |
| CN114151133A (en) * | 2021-12-14 | 2022-03-08 | 中煤科工开采研究院有限公司 | Grading early warning method for incoming pressure of top plate of fully mechanized mining face |
| CN115083134B (en) * | 2022-07-21 | 2022-12-02 | 中国长江三峡集团有限公司 | Hydraulic engineering dam flood discharge dissolved gas supersaturation early warning regulation and control system |
| CN117035197B (en) * | 2023-08-25 | 2024-06-04 | 成都理工大学 | Intelligent lost circulation prediction method with minimized cost |
-
2024
- 2024-05-10 CN CN202410571181.7A patent/CN118154046B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104252686A (en) * | 2014-08-15 | 2014-12-31 | 国家电网公司 | Determination method for power grid safety aggregative indicators |
| CN106703883A (en) * | 2016-12-29 | 2017-05-24 | 山东科技大学 | Method for determining floor water inrush danger level of coal mining working faces in personalized manner |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118154046A (en) | 2024-06-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111256754B (en) | Concrete dam long-term operation safety early warning method | |
| CN107977749B (en) | Prediction analysis method for oil field development index | |
| CN103233741B (en) | Method for quantitatively evaluating surrounding rock stability of ultra-large underground cavern group in construction period | |
| CN103292762B (en) | Displacement monitoring method for judging dam stability | |
| CN111563524A (en) | Multi-station fusion system operation situation abnormity monitoring and alarm combining method | |
| CN111042143A (en) | Foundation pit engineering early warning method and system based on analysis of large amount of monitoring data | |
| CN109918822A (en) | The shield driving parameter error calculation method and system returned based on random forest | |
| CN112145231A (en) | Early warning method for working condition of hydraulic support | |
| CN113780375B (en) | Digital twinning-based virtual-real interaction wind power plant wind power monitoring system and application | |
| CN114357691A (en) | Power facility geological foundation deformation safety assessment method | |
| CN112101811A (en) | Water supply network pipe explosion risk prediction method and system | |
| CN117454114B (en) | Subway tunnel tunneling blasting vibration safety monitoring device based on multi-point location distribution | |
| CN112861350A (en) | Temperature overheating defect early warning method for stator winding of water-cooled steam turbine generator | |
| CN118154046B (en) | Top plate pressure grade dividing method | |
| CN111476420A (en) | Method for optimizing early warning index of microseism monitoring rock burst | |
| CN113298426B (en) | Knowledge graph driven dam safety evaluation weight dynamic drafting method and system | |
| CN117271675A (en) | TBM tunneling rock mass digability identification method and device based on big data mining | |
| CN116304636A (en) | Elevator top-rushing accident dynamic prediction method and system based on fault tree | |
| CN115859205A (en) | Hydropower station risk early warning method and system based on space-time internet of things | |
| CN113158562A (en) | TBM rock machine mapping construction method and system based on physical constraint and data mining | |
| CN113586157B (en) | Extraction working face salient dangerous area rapid division method based on Kriging interpolation | |
| CN118798002A (en) | Intelligent water conservancy early warning system and method based on digital twinning | |
| KR102188096B1 (en) | Analasis method for evevt of sensor | |
| CN117474165A (en) | Ore pressure prediction method based on support vector machine | |
| CN118799104A (en) | Intelligent monitoring method and system for dynamic and supporting quality of top plate |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |