Disclosure of Invention
Therefore, the invention aims to solve the technical problems that the mine dynamic disaster monitoring and early warning method in the prior art cannot accurately judge and identify the position and the danger level of a local dangerous area, has poor early warning effect, and cannot pertinently guide a mine to perform disaster protection, and provides the mine disaster multi-parameter local dangerous area judgment and early warning method which can accurately judge and identify the position and the danger level of the local dangerous area, has good early warning effect, can pertinently guide the mine to perform disaster protection, and can meet the field requirements of various mines.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows:
a mine disaster multi-parameter local danger area judgment and early warning method comprises the following steps:
step 1, spatial positioning of mine disaster monitoring points
In the same coordinate system, carrying out space positioning on each monitoring point of each monitoring and early warning system;
step 2, obtaining the initial risk index I of each monitoring pointSingle point initialization(0≤ISingle point initialization≤1);
Step 3, retrieving the conditions of other monitoring points around each monitoring point
Obtaining the condition that other monitoring points are searched within the sphere area with the radius r by taking each monitoring point as a central measuring point;
step 4, correcting initial danger index of monitoring point
According to the range of the sphere area with the radius r of each central measuring point, the initial risk index I of the monitoring point is retrievedSingle point initializationInitial hazard index I for each central stationSingle point initializationCorrecting to obtain corrected danger index I of each central measuring pointSingle point correction;
Step 5, obtaining the position and the danger level of the local danger area
Selecting a modified risk index ISingle point correctionThe central measuring point with the maximum value, the sphere area with the radius r of the central measuring point is judged as a local danger area, and the corrected danger index I of the central measuring pointSingle point correctionThe risk level of the local risk zone.
Preferably, in the step 4, a correction coefficient k is assigned to each retrieved monitoring pointiInitial risk index I of the retrieved monitoring pointSingle point initializationThe higher the correction factor k it is assigned toiThe larger the value of (A);
modified hazard index I for each center measurement pointSingle point correction=ISingle point initialization*k1*k2*……*ki
Wherein, I in the above formulaSingle point initializationIs the initial danger index of the central measuring point, I is the number of the monitoring points to be retrieved, and when I is obtained by calculationSingle point correctionWhen the ratio is more than 1, take ISingle point correctionEqual to 1.
Preferably, the correction coefficient kiThe values of (A) are as follows:
preferably, step 2 is to provide corresponding monitoring data danger judgment threshold values for the monitoring data of different monitoring and early warning systems, and perform danger level assessment on each monitoring point according to the monitoring data danger judgment threshold values to obtain normalized initial danger index I of each monitoring pointSingle point initialization。
Preferably, the initial risk index I of each monitoring point in the step 3Single point initializationAnd obtaining the conditions of other monitoring points searched in the sphere area range with each monitoring point as a central measuring point and the radius of r from large to small.
Preferably, the radius r in step 3 is 50 meters.
Preferably, each monitoring and early warning system in the step 1 comprises a coal seam stress monitoring system, a microseismic monitoring system, a drilling cutting method monitoring system and a hydrological monitoring system.
Compared with the prior art, the technical scheme of the invention has the following advantages:
(1) the mine disaster multi-parameter local dangerous area judgment and early warning method provided by the invention is based on the existing single-parameter mine disaster monitoring method, can accurately judge and recognize the position and the dangerous grade of the local dangerous area by utilizing the monitoring data obtained by the existing various monitoring and early warning systems through operation and analysis, has good early warning effect, can pertinently guide the mine to carry out disaster protection, and can meet the field requirements of various mines.
(2) According to the mine disaster multi-parameter local dangerous area judgment and early warning method provided by the invention, the initial dangerous index of each monitoring point is corrected, so that the accuracy and reliability of judging and recognizing the position and the dangerous level of the local dangerous area are greatly improved.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1, the invention is a preferred embodiment of a multi-parameter local danger area judgment and early warning method for mine disasters. The mine disaster multi-parameter local danger area identification early warning method comprises the following steps:
step 1, spatial positioning of mine disaster monitoring points
And carrying out space positioning on each monitoring point of each monitoring and early warning system in the same coordinate system. Specifically, because the existing most mine disaster monitoring and early warning systems are not configured with coordinate information of monitoring points, the early warning method of the invention firstly needs to perform spatial positioning on each monitoring point of each monitoring and early warning system in the same coordinate system (for example, a 1980 national geodetic coordinate system commonly used in mines) according to the installation position (or the occurrence position of a microseismic event) of each monitoring point.
And (3) carrying out the following coordinate positioning on each monitoring point of each monitoring and early warning system:
note: the types of mine disaster monitoring data and the number of monitoring points in the above table can be expanded, which is not limited by the invention.
The invention utilizes various existing monitoring and early warning systems to carry out operation analysis, and the common monitoring and early warning systems comprise a coal seam stress monitoring system, a microseism monitoring system, a drilling cutting method monitoring system, a hydrological monitoring system and the like. The present embodiment only takes the partial monitoring data type as an example.
For example, the coordinates of two measuring points of stress monitoring "stress measuring point 1", "stress measuring point 2", microseismic monitoring "microseismic event 1", "microseismic event 2", drilling cutting method monitoring "construction hole 1" and "construction hole 2" are as follows:
step 2, obtaining the initial risk index I of each monitoring pointSingle point initialization
Providing corresponding monitoring data danger judgment threshold values aiming at monitoring data of different monitoring and early warning systems, and carrying out danger grade evaluation on each monitoring point according to the monitoring data danger judgment threshold values to obtain normalized initial danger index I of each monitoring pointSingle point initialization(0≤ISingle point initialization≤1)。
Because the monitoring data are more in types and the types can be expanded, the following description only takes the coal bed stress monitoring, the microseismic monitoring and the drilling cutting method monitoring as examples:
1. initial danger index (I) of coal body stress monitoring pointStress single point initiation)
(1) Calculating stress value index (I) of single measuring pointY)
Stress value index (I)Y) Giving four-level values according to the real-time stress value of the measuring point, wherein the depth hole value intervals are different, and the value markThe parameters are quasi adjustable and take the following values:
grade
|
Blue color
|
Yellow colour
|
Orange colour
|
Red colour
|
Default (MPa) shallow hole
|
[0,8)
|
[8,10)
|
[10,13)
|
[13,+∞)
|
Default value (MPa) deep hole
|
[0,9)
|
[9,12)
|
[12,15)
|
[15,+∞)
|
Stress value index (I)Y)
|
0
|
0.33
|
0.66
|
1 |
Note: the parameter values and parameter intervals in the above table can be reset according to practical application, which is not limited in the present invention.
For example, the stress values of the two shallow hole stress measuring points "stress measuring point 1" and "stress measuring point 2" are 10Mpa and 11Mpa respectively, the stress value index of "stress measuring point 1" is 0.66, and the stress value index of "stress measuring point 2" is 0.33.
(2) Calculating stress amplitude index (I) of single measuring pointI)
Increase value index (I)I) The standard is used for representing the increase amplitude between the current stress value and the initial installation value, the larger the increase amplitude is, the higher the represented danger level is, and the values are as follows:
the growth coefficient is [ (current stress value-initial stress value)/initial stress value ]. 100%
Grade
|
Blue color
|
Yellow colour
|
Orange colour
|
Red colour
|
Coefficient of growth
|
[-∞,50%)
|
[50%,80%)
|
[80%,100%)
|
[100,+∞)
|
Stress increase value index (I)I)
|
0
|
0.33
|
0.66
|
1 |
Note: the parameter values and parameter intervals in the above table can be reset according to practical application, which is not limited in the present invention.
For example, the initial stress values of the two shallow hole stress measuring points "stress measuring point 1" and "stress measuring point 2" are both 7Mpa, the stress increase coefficients of the two measuring points are 42.9% and 57.1%, the amplitude index of the "stress measuring point 1" is 0, and the amplitude index of the "stress measuring point 2" is 0.33.
(3) Calculating the stress acceleration value index (I) of the single measuring pointG)
Single-point stress acceleration index (I)G) The method is a standard for representing the rapid increase degree of the stress value in a short time, the faster the increase speed is, the higher the represented danger level is, and the calculation process is as follows:
① first, in order to avoid the influence of oil filling operation at the ground stress measuring point, it is judged whether the current stress is larger than the critical value YCritical increase of(default is 6MPa, adjustable parameter) when being less than YCritical increase ofIn time, the measurement point does not perform the acceleration index calculation.
② calculating the growth rate
Suppose that the current time is t2 and the stress value is Yt2T1 at some time before, and Y as the stress valuet1The speed increasing rate:
rate of increase of 100% (Y)t2-Yt1)/(Yt1(t2-t1))
In the formula: the default time interval is 1 hour (i.e., 1 hour from t2-t 1), which may be adjusted
③ acceleration value index (I)G) The default values and calculation method are as follows:
grade
|
Blue color
|
Yellow colour
|
Orange colour
|
Red colour
|
Rate of growth
|
[-∞,10%)
|
[10%,15%)
|
[15%,30%)
|
[30%,+∞)
|
Index of acceleration value (I)G)
|
0
|
0.33
|
0.66
|
1 |
Note: the parameter values and parameter intervals in the above table can be reset according to practical application, which is not limited in the present invention.
For example, when the default values of t2-t1 are 1 hour, the stress values of the two shallow hole stress measuring points "stress measuring point 1" and "stress measuring point 2" before 1 hour are respectively 9.5Mpa and 9Mpa, the growth speeds of the two measuring points are respectively 5.3% and 22.2%, the acceleration value index of "stress measuring point 1" is 0, and the amplification value index of "stress measuring point 2" is 0.66.
(4) Calculating the initial danger index (I) of the coal body stress monitoring pointStress single point initiation)
Calculating the initial danger index (I) of the coal body stress monitoring pointStress single point initiation) Calculating three indexes of stress value, amplification and acceleration of the stress single-measuring point in real time, and defining weight coefficients K of the three indexes according to the region display characteristicsStress1~KStress 3(∑KStress i1, adjustable parameter), the calculation process is as follows:
Istress single point initiation=KStress 1*IY+KStress 2*II+KStress 3*IG
In the formula: i isYSingle point stress value index, KStress 1Default value 0.5; i isI-a single-point stress increase value indicator,
Kstress 2Default value 0.2; i isGSingle point stress acceleration value index, KStress 3Default value 0.3
For example, KStress 1~KStress 3Taking default values, and initial danger indexes (I) of two shallow hole stress measuring points ' stress measuring point 1 ' and ' stress measuring point 2Stress single point initiation):
IStress single point initiation 1=0.5*0.66+0.2*0+0.3*0=0.33
IStress single point initiation 2=0.5*0.33+0.2*0.33+0.3*0.66=0.43
2. Microseismic event monitoring point initial risk index (I)Microseismic single point initiation)
The early warning of a local dangerous area has certain requirements on the timeliness of microseism monitoring data, microseism events rolling for 24 hours (adjustable parameters) at the current moment are taken as basic data participating in local early warning, and the danger index (I) of each microseism event monitoring point is determined according to the energy of each microseism event in the current time periodMicroseismic single point initiation) The value method is as follows: (energy less than 10)3J microseismic events do not participate in the calculation)
Index/grade
|
Blue color
|
Yellow colour
|
Orange colour
|
Red colour
|
Energy interval/J
|
[103,104)
|
[104,105)
|
[105,106)
|
[106,+∞)
|
Microseismic event monitoring point hazard index (I)Microseismic single point initiation)
|
0
|
0.33
|
0.66
|
1 |
Note: the parameter values and parameter intervals in the above table can be reset according to practical application, which is not limited in the present invention.
For example, a "microseismic event 1" energy of 1.2 × 10 within 24 hours4J, its initial hazard index IMicroseismic single point initiation 1The table above is taken to be 0.33 (yellow) and the "microseismic event 2" energy occurring within 24 hours is 1.2 × 103J, the hazard index I thereofMicroseismic single point initiation 2The table is taken 0 (blue) for the control.
3. Initial danger index (I) of drill cuttings monitoring pointsSingle point initiation of drill cuttings)
Taking the drilling inspection result which rolls for 24 hours (adjustable parameters) at the current moment as basic data participating in local early warning, and judging the initial danger index (I) of each drilling cutting monitoring point according to whether the drilling inspection data has dynamic display and exceeds the standard or notSingle point initiation of drill cuttings) The values of (A) are as follows:
note: the parameter values and parameter intervals in the above table can be reset according to practical application, which is not limited in the present invention.
For example, within 24 hours, the situation that the power is shown and the coal powder is not overproof occurs in a construction hole 1 of a certain chip method monitoring point, and the initial danger index I of the monitoring point is judgedSingle point initiation of drill cuttings 1Is 0.66. The situation that 'no power display and no standard exceeding of coal powder' occur in 'construction hole 2' of another chip method monitoring point is judged, and the initial danger index I of the monitoring point is judgedSingle point initiation of cuttings 2Is 0.
Step 3, retrieving the conditions of other monitoring points around each monitoring point
And obtaining the condition that other monitoring points are searched in the sphere area range with the radius r by taking each monitoring point as a central measuring point.
In the specific retrieval process, the initial risk index I of each monitoring point can be usedSingle point initializationAnd searching in the order from big to small to prevent missing detection, or preferably searching by using monitoring points with larger initial danger indexes after sequencing. The radius r is an empirical value, and in this embodiment, the radius r takes the value of 50 meters.
For example, the initial risk index for each monitoring point is ranked as:
the stress measuring point 2, the drilling cuttings construction hole 1, the stress measuring point 1, the microseismic measuring point 2 and the drilling cuttings construction hole 2.
According to the sequencing result, each monitoring point is sequentially taken as a central measuring point, and the condition results of other monitoring points are retrieved as follows:
step 4, correcting initial danger index of monitoring point
According to the range of the sphere area with the radius r of each central measuring point, the initial risk index I of the monitoring point is retrievedSingle point initializationInitial hazard index I for each central stationSingle point initializationCorrecting to obtain corrected danger index I of each central measuring pointSingle point correction。
Specifically, a correction coefficient k is assigned to each retrieved monitoring pointiInitial risk index I of the retrieved monitoring pointSingle point initializationThe higher the correction factor k it is assigned toiThe larger the value of (A); the correction coefficient kiThe values of (A) are as follows:
note: the parameter values and parameter intervals in the above table can be reset according to practical application, which is not limited in the present invention.
Modified hazard index I for each center measurement pointSingle point correction=ISingle point initialization*k1*k2*……*ki
Wherein, I in the above formulaSingle point initializationIs the initial danger index of the central measuring point, I is the number of the monitoring points to be retrieved, and when I is obtained by calculationSingle point correctionWhen the ratio is more than 1, take ISingle point correctionEqual to 1.
For example, the corrected danger index I of each center pointSingle point correctionThe results were:
step 5, obtaining the position and the danger level of the local danger area
According to the corrected danger index I of each central measuring pointSingle point correctionAnd judging the local danger area and the danger level.
Specifically, a modified risk index I is selectedSingle point correctionThe central measuring point with the maximum value, the sphere area with the radius r of the central measuring point is judged as a local danger area, and the corrected danger index I of the central measuring pointSingle point correctionThe risk level of the local risk zone. For a clearer and more intuitive determination, the corrected risk index I of all central measuring points can be determined firstSingle point correctionThe risk levels of the local dangerous areas are judged according to the following table in a descending order.
Hazard class
|
Blue color
|
Yellow colour
|
Orange colour
|
Red colour
|
Modified hazard index ISingle point correction |
[0,0.3)
|
[0.3,0.5)
|
[0.5,0.8)
|
[0.8,1] |
Note: the parameter values and parameter intervals in the above table can be reset according to practical application, which is not limited in the present invention.
For example, the modified risk index I for each monitoring pointSingle point correctionThe results were ranked as follows:
the sphere area within the radius range of 50m of the central measuring point stress measuring point 1 with the maximum corrected danger index is judged as a local danger area from the table, the danger index of the local danger area is 0.74, and the sphere area has an orange danger level.
In the present invention, the "adjustable parameter" means that each parameter and parameter range shown in the embodiment of the present invention can be reset according to the actual application, and the present invention does not limit these specific values and ranges.
In other embodiments, the type of the monitored data according to which the warning method of the present invention is based can be increased, decreased, or re-formulated according to the actual application.
It should be understood that the above examples are only for clarity of illustration and are not intended to limit the embodiments. Other variations and modifications will be apparent to persons skilled in the art in light of the above description. And are neither required nor exhaustive of all embodiments. And obvious variations or modifications therefrom are within the scope of the invention.