CN109239775B - Method for tracking and positioning stolen mining of mineral resources - Google Patents
Method for tracking and positioning stolen mining of mineral resources Download PDFInfo
- Publication number
- CN109239775B CN109239775B CN201811089429.7A CN201811089429A CN109239775B CN 109239775 B CN109239775 B CN 109239775B CN 201811089429 A CN201811089429 A CN 201811089429A CN 109239775 B CN109239775 B CN 109239775B
- Authority
- CN
- China
- Prior art keywords
- blasting
- event
- mining
- events
- analysis
- 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
- 238000005065 mining Methods 0.000 title claims abstract description 104
- 238000000034 method Methods 0.000 title claims abstract description 26
- 229910052500 inorganic mineral Inorganic materials 0.000 title claims abstract description 15
- 239000011707 mineral Substances 0.000 title claims abstract description 15
- 238000005422 blasting Methods 0.000 claims abstract description 166
- 238000004458 analytical method Methods 0.000 claims abstract description 55
- 230000000694 effects Effects 0.000 claims abstract description 42
- 238000012544 monitoring process Methods 0.000 claims abstract description 39
- 238000004519 manufacturing process Methods 0.000 claims abstract description 23
- 238000012300 Sequence Analysis Methods 0.000 claims abstract description 12
- 239000011435 rock Substances 0.000 claims description 35
- 230000002776 aggregation Effects 0.000 claims description 17
- 238000004220 aggregation Methods 0.000 claims description 17
- 238000001228 spectrum Methods 0.000 claims description 9
- 238000009412 basement excavation Methods 0.000 claims description 8
- 238000009826 distribution Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 230000004044 response Effects 0.000 claims description 3
- 230000003203 everyday effect Effects 0.000 abstract description 3
- 238000012216 screening Methods 0.000 abstract description 3
- 238000007405 data analysis Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 6
- 230000001960 triggered effect Effects 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000006378 damage Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 238000005553 drilling Methods 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 230000005641 tunneling Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000002354 daily effect Effects 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005180 public health Effects 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for tracking and positioning stolen mining of mineral resources, which comprises the steps of installing a micro-seismic monitoring system in a target area, and monitoring micro-seismic events and blasting events in the target area through the micro-seismic monitoring system. And selecting a fixed period, such as every day, identifying and screening the illegal mining blasting events, and eliminating microseismic events and blasting events generated by normal production of mines so as to obtain the illegal mining movable blasting events. And performing space-time sequence analysis and space clustering analysis on the well-defined illegal mining blasting event. And finally, according to the data analysis result, combining with an engineering model or an ore body model on site to explain in detail, noting the three-dimensional coordinates of the gathering reference point of the illegal mining blasting event, the spatial relation between the gathering reference point and the production operation and the activity time period, and obtaining the motion trail of the illegal mining operation surface through multiple times of analysis and observation so as to track and position the illegal mining activity in real time.
Description
Technical Field
The invention relates to the field of mineral resource safety, in particular to a method for tracking and positioning stolen mining of mineral resources.
Background
Mineral resources are one of the most important material resources in national social and economic life, but illegal stealing of the mineral resources exists in various places, so that a plurality of serious safety accidents are caused. The stealing mining activity of mineral resources brings serious threats to the daily production of a mining area, and if a large amount of unknown goafs are left, the frequent occurrence of the ground pressure activity, the surface subsidence and the roof fall accidents of the mining area can be caused, so that safety accidents, the waste of resources and the damage of the environment can be caused. Due to the concealment of the stealing activities, the specific areas and the stealing tracks of the stealing activities are difficult to obtain.
Therefore, there is a need to develop a technical method for mastering the characteristic information of the illegal mining activity.
Disclosure of Invention
Objects of the invention
The invention aims to solve the technical problem that the characteristic information of the illegal mining activity cannot be mastered due to the illegal mining of mineral resources, and provides a method for realizing the characteristic information acquisition and data processing and analysis of the illegal mining activity.
(II) technical scheme
In order to achieve the purpose, the invention adopts the main technical scheme that:
a method for tracking and positioning stolen mining of mineral resources comprises the following steps:
s1, installing a mining micro-seismic monitoring system in the target area, and monitoring and recording rock mass vibration events in the target area in real time;
step S2, classifying the recorded rock mass vibration events to obtain a blasting event set and a stealing mining blasting event set;
step S3, performing space-time sequence analysis and space clustering analysis on the illegal mining blasting event set;
step S4, acquiring the three-dimensional coordinates of the gathering reference point of the illegal mining blasting event, the time period of the activity and the spatial relationship with the normal mine production operation based on the analysis result of the step S3 and the target area model;
and S5, repeating the steps S1-S4, obtaining the motion track, the propulsion speed, the blasting period and the excavation amount of the illegal mining operation, and realizing real-time tracking, positioning and monitoring of the illegal mining activities.
The step S1 may further include,
after the mining micro-seismic monitoring system is installed in the target area, the wave speed of the system is repeatedly calibrated by a fixed-point gun coordinate calibration method, and the positioning accuracy of the system is debugged until the positioning accuracy of the micro-seismic monitoring system is less than 6 m.
Preferably, the mining microseismic monitoring system has a sampling rate of more than 10KHz, a dynamic range of more than 115dB, a sensor response frequency range of 0.5Hz to 6000Hz, time synchronization precision of less than 1 microsecond and event positioning error of less than 6 m.
The step S2 includes:
s2a, performing characteristic analysis on the recorded rock mass vibration event to obtain a blasting event set;
and S2d, performing seismic source parameter characteristic comparison analysis on the blasting event set, and eliminating blasting events generated by normal production of the mine to obtain an illegal mining blasting event set.
The step S2a includes at least one of the following steps:
analyzing the waveform characteristics of each recorded rock mass vibration event, and if the waveform characteristics of the vibration event exist are that multiple wave crests continuously appear and the S wave characteristics are not obvious, marking the vibration event as an explosion event;
performing spectral feature analysis on each recorded rock mass vibration event, and if the frequency spectrum high-frequency part occupation ratio of the vibration event is far greater than that of the low-frequency part and the corner frequency is greater than 500Hz, recording the vibration event as a blasting event;
carrying out dragon distribution model analysis on each recorded rock mass vibration event, if the vibration event conforms to the dragon distribution model, taking the vibration event information as microseismic event information, and rejecting the vibration event;
and converting the recorded waveform of each rock mass vibration event into a sound file, performing event sound characteristic analysis, and recording the vibration event as a blasting event if the sound characteristics of the vibration event meet the sound characteristics of the blasting event.
The step S3 includes:
s3a, acquiring seismic source parameters of the illegal mining blasting event;
s3b, performing space-time sequence analysis on the illegal mining blasting event set based on the seismic source parameters to obtain the occurrence place aggregation condition of the illegal mining blasting event and the activity frequency change on a time axis;
and S3c, calculating the distance between blasting sources of every two illegal mining blasting events, taking the occurrence position of each blasting event as a clustering object, and carrying out clustering analysis on the blasting events to obtain the aggregation points of the blasting events.
The step S3c includes:
step 3c1, obtaining the location information of the illegal mining blasting event obtained in the step 2, and calculating the distance between every two blasting sources;
step 3c2, taking the occurrence position of the blasting event as a clustering object, wherein all blasting sources are classified into one type, and obtaining the spatial three-dimensional coordinates (x, y, z) of the blasting position as clustering elements;
the calculation formula of the clustering elements comprises:
xik＝(xi1,xi2,xi3)，
wherein, i is 1,2,3,. n, n is the number of blasting sources;
step 3c3, two blasting sources x with the minimum distanceikAnd xjkAnd doing the same, calculating xikAnd xjkA distance d therebetweenij，
The calculation formula comprises:
and 3c4, calculating the distance between the new class and the other classes, merging the two closest classes, and repeating the steps until all blasting events are classified as 1 if the number of the classes is still greater than 1.
The step S5 includes:
s5a, setting a distance threshold, repeating the steps S1-S4, and obtaining the strip-shaped motion direction of the multiple clustering events in space-time, which is the space motion track of the stealing activities;
s5b, obtaining a blasting cycle through each clustering result and event occurrence time interval;
and S5c, calculating the excavation amount through the seismic source parameter accumulated visual volume and the rock mass elastic modulus parameter.
Preferably, the distance threshold is 6 m.
(III) advantageous effects
The invention has the beneficial effects that: the method is characterized in that a detailed explanation is carried out by combining an engineering model or an ore body model on site, three-dimensional coordinates of an aggregation reference point of a illegal mining blasting event, a spatial relation between the aggregation reference point and production operation and an activity time period are noted, a motion track of an illegal mining operation surface can be obtained through multiple times of analysis and observation, then the illegal mining activity can be tracked and positioned in real time, and an effective guiding basis is provided for mine enterprise units and government administration illegal mining activity to protect mineral resources.
And installing a micro-seismic monitoring system in the target area, and monitoring a micro-seismic event and a blasting event in the target area through the micro-seismic monitoring system. And selecting a fixed period, such as every day, identifying and screening the illegal mining blasting events, and eliminating microseismic events and blasting events generated by normal production of mines so as to obtain the illegal mining movable blasting events. And performing space-time sequence analysis and space clustering analysis on the well-defined illegal mining blasting event. And finally, according to the data analysis result, combining with a field engineering model or a mineral body model to explain in detail, noting the three-dimensional coordinates of the gathering reference point of the illegal mining blasting event, the spatial relation between the gathering reference point and the production operation and the activity time period, and obtaining the motion trail of the illegal mining operation surface through multiple times of analysis and observation, thereby carrying out real-time tracking positioning and quantitative statistics on the illegal mining activity and providing an effective guide basis for mine enterprise units and governments to control the illegal mining activity and protect mineral resources.
Drawings
FIG. 1 is a topological diagram of a high-precision microseismic monitoring system according to an embodiment of the present invention;
FIG. 2 is a diagram illustrating waveform characteristics of a vibration signal and classification of sound files according to an embodiment of the present invention;
FIG. 3 is a schematic diagram illustrating the aggregation of blast events according to an embodiment of the present invention;
FIG. 4 is a diagram illustrating changes in activity frequency over a time axis of a blasting event obtained by spatio-temporal sequence analysis in accordance with an embodiment of the present invention;
fig. 5 is a diagram illustrating a clustering result of a theft mining blasting event obtained by spatial clustering analysis according to an embodiment of the present invention.
Detailed Description
For a better understanding of the present invention, reference will now be made in detail to the present invention by way of specific embodiments thereof.
The embodiment of the invention relates to a method for tracking and positioning stolen mining of mineral resources, which specifically comprises the following steps:
and S1, installing the mining micro-seismic monitoring system in the target area, and monitoring and recording the rock mass vibration event in the target area in real time.
A set of mining micro-seismic monitoring system is installed in the target area, rock mass vibration events in the target area are monitored in real time, and the rock mass vibration events triggering the monitoring system are stored in a local server.
Preferably, the system is required to have a sampling rate of more than 10KHz, a dynamic range of more than 115dB, a sensor response frequency range of 0.5Hz to 6000Hz, time synchronization accuracy of less than 1 microsecond and event positioning error of less than 6 m.
The sampling rate is the system sampling frequency, and is more than 2 times higher than the maximum frequency of a signal according to the signal acquisition requirement, and the higher the sampling rate is, the stronger the signal resolution capability is;
the dynamic range refers to the ratio of the maximum value to the minimum value of the variable signal, and the larger the dynamic range is, the better the system performance is;
the time synchronization precision refers to the time synchronization error of each device in the system, and the smaller the synchronization error is, the better the synchronization error is;
the positioning error refers to the difference distance between the system positioning coordinate position and the actual seismic source coordinate position, and the positioning error is required to be smaller than the seismic source fracture scale.
A set of high-precision microseismic monitoring system meeting the requirements is installed in a target area, the system architecture diagram is shown in figure 1, the sensor is installed in a drilling and grouting mode, the drilling depth is required to be 3-5 times larger than the diameter of a roadway, so that the noise of the vibration waveform of the rock received by the sensor is less, and the waveform quality is effectively guaranteed; the acquisition substation is provided with a battery, so that the problem that information acquisition and transmission cannot be performed due to the fact that a system cannot supply power under the condition of power failure is prevented; a time service device of a GPS or BDS is configured on the earth surface, and a network ptp time synchronization protocol is adopted to ensure that the time synchronization error between each device of the system is less than 1 microsecond.
Preferably, step S1 further includes, after the mining microseismic monitoring system is installed in the target area, repeatedly calibrating the system wave velocity by using a fixed-point gun coordinate calibration method, and performing system positioning accuracy debugging until the positioning accuracy of the microseismic monitoring system is less than 6 m.
And the system is debugged according to the fact that the system is debugged after the system is debugged, the main performance index of debugging is system positioning accuracy, the system wave speed is repeatedly calibrated by a fixed-point gun calibration coordinate method, and finally the elastic wave propagation speed capable of accurately reflecting rock mass characteristics in a monitoring area is obtained, and the event positioning accuracy is guaranteed to be smaller than 6 m.
Step S2, classifying the recorded rock mass vibration events to obtain a blasting event set and a stealing mining blasting event set;
and performing waveform characteristic analysis, spectrum characteristic analysis, Brune (Brune) model analysis and event sound characteristic analysis on the triggered events, so as to classify the triggered events, wherein the triggered events comprise microseismic events, mechanical vibration events and blasting events, mine production blasting records are collected, and blasting events generated by normal production of mines are removed by comparing the seismic source parameter characteristics of the blasting events, so that blasting events generated by abnormal production are obtained.
The blasting event is determined to be a blasting event by a microseismic monitoring system by adopting an explosive blasting mode to break rock and further generate strong rock vibration for mining ores;
the microseism event is a destruction phenomenon generated under the action of extrusion and shearing force of rocks, the destruction phenomenon of the rocks can also generate elastic waves which are detected by a microseism monitoring system, but the waveform characteristics and the blasting waveform characteristics have obvious difference;
a theft blast event refers to a blast event other than a normal mine mining blast event.
And S2a, performing characteristic analysis on the recorded rock mass vibration event to obtain a blasting event set.
The method for identifying the blasting event and the microseismic event comprises the following steps: the type of the blasting event is identified by adopting waveform characteristic identification, frequency spectrum analysis, Brune model analysis and event sound characteristic analysis for the event triggered by the microseismic monitoring system.
The step S2a includes at least one of the following steps:
analyzing the waveform characteristics of each recorded rock mass vibration event, and if the waveform characteristics of the vibration event exist are that multiple wave crests continuously appear and the S wave characteristics are not obvious, marking the vibration event as an explosion event;
the recorded rock mass vibration events can be classified by analyzing the waveform characteristics of the recorded rock mass vibration events, such as P wave analysis and S wave analysis, the blasting waveform characteristics have multiple continuous wave crests, and the blasting waveform S wave characteristics are not obvious.
Performing spectral feature analysis on each recorded rock mass vibration event, and if the frequency spectrum high-frequency part occupation ratio of the vibration event is far greater than that of the low-frequency part and the corner frequency is greater than 500Hz, recording the vibration event as a blasting event;
by performing Fourier transform on the event waveform, the spectral characteristics of the blasting event can be obtained, the high-frequency part of the blasting event spectrum has more proportion, and the corner frequency is high and is generally more than 500 Hz.
Carrying out dragon distribution model analysis on each recorded rock mass vibration event, if the vibration event conforms to the dragon distribution model, taking the vibration event information as microseismic event information, and rejecting the vibration event;
judging whether the frequency spectrum curve conforms to a bloom model or not, wherein the frequency spectrum curve of the general mine microseismic event conforms to the bloom (Brune) model, and the blasting event does not conform to the model;
the Brune model was a seismic source model proposed in the 1970 s, which assumed that the source failure was shear failure, with continuous slip velocities throughout the failure region, with significant variation characteristics in the spectral analysis of the waveform relative to the waveform of non-microseismic events.
Converting the recorded waveform of each rock mass vibration event into a sound file, performing event sound characteristic analysis, and recording the vibration event as a blasting event if the sound characteristics of the vibration event meet the sound characteristics of the blasting event;
the events monitored by the system can be converted into sound files through waveforms, and blasting events can be distinguished through sound characteristics. The event sound feature analysis identifies the blasting event type, as shown in FIG. 2.
And S2d, performing seismic source parameter characteristic comparison analysis on the blasting event set, and eliminating blasting events generated by normal production of the mine to obtain an illegal mining blasting event set.
And determining the abnormal operation blasting event according to the blasting event source parameters so as to finally identify the illegal mining operation blasting event, acquiring monitoring data from a mine installation microseismic monitoring system, identifying and removing the monitoring data and delineating the illegal mining blasting event in the period as shown in fig. 3.
Or, through time period analysis, mine production generally has a fixed blasting time period, and by eliminating blasting events in the time period, the rest blasting events can be regarded as illegal mining blasting events.
And selecting a fixed period, such as every day, identifying and screening the illegal mining blasting events, and eliminating microseismic events and blasting events generated by normal production of mines so as to obtain the illegal mining movable blasting events.
And step S3, performing space-time sequence analysis and space clustering analysis on the illegal mining blasting event set.
The abnormal production blasting event is the stealing mining blasting operation activity in the monitoring area, and the confined stealing mining blasting event is subjected to seismic source parameter characteristic analysis, space-time sequence analysis and space clustering analysis.
And S3a, acquiring the seismic source parameters of the illegal mining blasting event.
The parameters of the blasting event source refer to the physical characteristics of the blasting event source, and a primary vibration event is generally described by time, space and strength.
And S3b, performing space-time sequence analysis on the illegal mining blasting event set based on the seismic source parameters to obtain the occurrence place aggregation condition of the illegal mining blasting event and the activity frequency change on a time axis.
The time-space sequence analysis of the illegal mining blasting event refers to the aggregation of the occurrence sites of the blasting event and the change of the activity frequency on a time axis.
And S3c, calculating the distance between blasting sources of every two illegal mining blasting events, taking the occurrence position of each blasting event as a clustering object, and carrying out clustering analysis on the blasting events to obtain the aggregation points of the blasting events.
Spatial clustering analysis is one of the main means of spatial data mining and knowledge discovery, and has been widely applied to various fields such as geography, geology, meteorology, cartography, astronomy, public health and the like. By using the characteristic that natural earthquake space distribution often shows non-spherical aggregation, the related experience of space clustering analysis is applied, and a clustering method is adopted to process the aggregation shape with irregular blasting aiming at mine blasting event data.
The step S3b includes:
and 3c1, acquiring the site information of the illegal mining blasting event obtained in the step 2, and calculating the distance between every two blasting sources.
The principle of the spatial clustering analysis method is to calculate the distance between every two blasting sources, and select the distance to accord with the set threshold value through comparison.
And 3c2, taking the occurrence positions of the blasting events as clustering objects, wherein all the blasting sources form a class, and obtaining the spatial three-dimensional coordinates (x, y and z) of the blasting positions as clustering elements.
Taking the occurrence position of the blasting event as a clustering object, wherein the clustering element is a space three-dimensional coordinate (x, y, z) of the blasting event, and calculating according to the following formula:
xik＝(xi1,xi2,xi3)(i＝1,2,3,...n)
step 3c3, two blasting sources x with the minimum distanceikAnd xjkAnd doing the same, calculating xikAnd xjkA distance d therebetweenij。
At the start of the calculation, the shots are classified into one type (n types are shared at this time if the number of events is n). And calculating the distance between the blasting sources, measuring the degree of affinity and sparseness of the blasting sources according to the distance between the blasting sources, and combining the two closest blasting sources into one class. By dijTo represent xikAnd xjkThe calculation formula of the distance between the two is as follows:
and 3c4, calculating the distance between the new class and the other classes, merging the two closest classes, and repeating the steps until all blasting events are classified as 1 if the number of the classes is still greater than 1.
And (3) clustering by adopting a shortest distance method, namely calculating the distance between the new class and the other classes, merging the two closest classes, and repeating the steps until all blasting events are classified as 1 class if the number of the classes is still greater than 1. Wherein events of class Gq and class Gp, drkThe distance is defined as:
drk＝min{dpk,dqk}
and clustering and analyzing the blasting events monitored by the microseismic monitoring system to find out the gathering place of the blasting events.
And performing seismic source parameter characteristic analysis, space-time sequence analysis and spatial clustering analysis on the confined illegal mining blasting event, as shown in fig. 4 and 5. The analysis period has 66 stealing mining blasting positioning events which are mainly distributed on the footwall of the small ore body, the footwall of the large ore body and two sides. Mainly centralizes near the middle section of tunneling. From the time of occurrence, 2: 00-4: 00 am, 10: 00-15: 00 am, 20: 00-21 pm: about 00 o' clock.
According to the spatial distribution condition of the illegal mining blasting events in the analysis period, the aggregation degree of the illegal mining events is general on the whole, a certain redundancy degree of an aggregation scale relative to a system error is considered, 6m (an empirical value) is selected as a clustering unit scale, and a parameter threshold value is 2 (an empirical value);
and step S4, acquiring the three-dimensional coordinates of the gathering reference point of the illegal mining blasting event, the time period of the activity and the spatial relationship with the normal mine production operation based on the analysis result of the step S3 and the target area model.
And (3) based on the analysis result of the step (3), a detailed description is carried out by combining an engineering model or an ore body model on site, three-dimensional coordinates of an aggregation reference point of the illegal mining blasting event, a spatial relation between the aggregation reference point and production operation and an activity time period are noted, and then the operation surface position of the illegal mining activity is mastered.
The detailed description is carried out by combining with an engineering model or an ore body model on site, and the three-dimensional coordinates of the gathering reference point of the illegal mining blasting event and the space relation and the activity time period between the gathering reference point and the production operation are noted, and the following table shows that:
as a specific example of step S4, according to the analysis result of step three, a detailed description is made in conjunction with the engineering model or ore body model of the site, and the three-dimensional coordinates of the aggregation reference point of the illegal mining blasting event and the spatial relationship and the activity time period with the production operation are noted as shown in the following table:
and S5, repeating the steps S1-S4, obtaining the motion track, the propulsion speed, the blasting period and the excavation amount of the illegal mining operation, and realizing real-time tracking, positioning and monitoring of the illegal mining activities.
The forward propulsion motion trail of the illegal mining operation surface can be obtained by repeating the processes of the step S1, the step S2 and the step S3 for multiple times, and the propulsion speed, the blasting period and the blasting energy of the illegal mining activity can be tracked, positioned and monitored in real time through quantitative statistical analysis, so that an effective and reliable guide basis is provided for the treatment of the illegal mining activity.
And S5a, setting a distance threshold, repeating the steps S1-S4, and obtaining the zonal motion direction of the multiple clustering events in space-time, which is the space motion track of the stealing activities.
By setting a reasonable distance threshold value, preferably, the distance threshold value is 6m, namely the maximum breaking dimension of the seismic source of the illegal mining blasting, and by repeatedly carrying out cluster analysis on the illegal mining blasting events for multiple times, each cluster is regarded as one tunneling and propelling activity of the illegal mining blasting, the clustering events for multiple times can generate a banded motion direction in space-time, and the banded motion direction is a space motion track of the illegal mining activity.
And step S5b, obtaining the blasting cycle through each clustering result and event occurrence time interval.
And determining the blasting cycle through the clustering result and the event occurrence time period.
And S5c, calculating the excavation amount through the seismic source parameter accumulated visual volume and the rock mass elastic modulus parameter.
And calculating the illegal mining excavation amount according to the release energy of the blasting seismic source, wherein the excavation amount can be calculated through the accumulated apparent volume of the seismic source parameters and the elastic modulus parameters of the rock mass.
When the monitoring data is insufficient, the illegal mining activity track is difficult to describe accurately, and the monitoring data needs to be analyzed for a long time for many times.
In order to better track and position the illegal mining activity, a high-precision microseismic monitoring system and detailed mine blasting record data need to be equipped for reference analysis.
And finally, mine enterprises can carry out on-site investigation and confirmation according to the position of the illegal mining activity, meanwhile, the production safety management of the area is enhanced, and the safety of personnel and equipment is guaranteed.
It is to be understood that the invention is not limited to the specific arrangements and instrumentality described above and shown in the drawings. A detailed description of known methods is omitted herein for the sake of brevity. In the above embodiments, several specific steps are described and shown as examples. However, the method processes of the present invention are not limited to the specific steps described and illustrated, and those skilled in the art can make various changes, modifications and additions or change the order between the steps after comprehending the spirit of the present invention.
Finally, it should be noted that: the above-described embodiments are only for illustrating the technical solution of the present invention, and are not to be construed as limiting the same; although the present 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 solutions described in the foregoing embodiments may be modified, or some or all of the technical features may be equivalently replaced; the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (5)
1. A method for tracking and positioning stolen mining of mineral resources is characterized by comprising the following steps:
s1, installing a mining micro-seismic monitoring system in the target area, and monitoring and recording rock mass vibration events in the target area in real time;
step S2, classifying the recorded rock mass vibration events to obtain a blasting event set and a stealing mining blasting event set;
step S3, performing space-time sequence analysis and space clustering analysis on the illegal mining blasting event set;
step S4, acquiring the three-dimensional coordinates of the gathering reference point of the illegal mining blasting event, the activity time period and the spatial relation with the normal mine production operation based on the analysis result of the step S3 and the target area model;
s5, repeating the steps S1-S4, obtaining the motion track, the propulsion speed, the blasting period and the excavation amount of the illegal mining operation, and realizing real-time tracking, positioning and monitoring of the illegal mining activities;
step S2 includes:
s2a, performing characteristic analysis on the recorded rock mass vibration event to obtain a blasting event set;
s2d, performing seismic source parameter characteristic comparison analysis or passing time period analysis on the blasting event set, and eliminating blasting events generated by normal production of the mine to obtain an illegal mining blasting event set, wherein the passing time period analysis is based on that the blasting events generated by normal production have a fixed blasting time period;
the step S3 includes:
s3a, acquiring seismic source parameters of the illegal mining blasting event;
s3b, performing space-time sequence analysis on the illegal mining blasting event set based on the seismic source parameters to obtain the occurrence place aggregation condition of the illegal mining blasting event and the activity frequency change on a time axis;
s3c, calculating the distance between blasting sources of two illegal mining blasting events, taking the occurrence position of each blasting event as a clustering object, and carrying out clustering analysis on the blasting events to obtain an aggregation point of the blasting events;
the step S5 includes:
s5a, setting a distance threshold, repeating the steps S1-S4, and obtaining the zonal motion direction of the multiple clustering events in space-time, wherein the zonal motion direction is the space motion track of the stealing mining activity, and the distance threshold is 6 m;
s5b, obtaining a blasting cycle through each clustering result and event occurrence time interval;
and S5c, calculating the accumulated apparent volume of the seismic source parameters and the elastic modulus parameters of the rock mass to obtain the excavation amount.
2. The method according to claim 1, wherein the step S1 further comprises,
after the mining micro-seismic monitoring system is installed in the target area, the wave speed of the system is repeatedly calibrated by a fixed-point gun coordinate calibration method, and the positioning accuracy of the system is debugged until the positioning accuracy of the micro-seismic monitoring system is less than 6 m.
3. The method of claim 1,
the mining microseismic monitoring system has the sampling rate of more than 10KHz, the dynamic range of more than 115dB, the sensor response frequency range of 0.5Hz to 6000Hz, the time synchronization precision of less than 1 microsecond and the event positioning error of less than 6 m.
4. The method according to claim 1, wherein the step S2a comprises at least one of the following steps:
analyzing the waveform characteristics of each recorded rock mass vibration event, and if the waveform characteristics of the vibration event exist are that multiple wave crests continuously appear and the S wave characteristics are not obvious, marking the vibration event as a blasting event;
carrying out spectrum characteristic analysis on each recorded rock mass vibration event, and if the frequency spectrum high-frequency part of the vibration event is far larger than the frequency spectrum low-frequency part and the corner frequency is larger than 500Hz, recording the vibration event as a blasting event;
carrying out dragon distribution model analysis on each recorded rock mass vibration event, if the vibration event conforms to the dragon distribution model, taking the vibration event information as microseismic event information, and rejecting the vibration event;
and converting the recorded waveform of each rock mass vibration event into a sound file, performing event sound characteristic analysis, and recording the vibration event as a blasting event if the sound characteristics of the vibration event meet the sound characteristics of the blasting event.
5. The method according to claim 1, wherein the step S3c comprises:
step 3c1, obtaining the location information of the illegal mining blasting event obtained in the step 2, and calculating the distance between every two blasting sources;
step 3c2, taking the occurrence position of the blasting event as a clustering object, wherein all blasting sources are classified into one type, and obtaining the spatial three-dimensional coordinates (x, y, z) of the blasting position as clustering elements;
the calculation formula of the clustering elements comprises:
xik＝(xi1,xi2,xi3)，
wherein, i is 1,2,3,. n, n is the number of blasting sources;
step 3c3, two blasting sources x with the minimum distanceikAnd xjkAnd doing the same, calculating xikAnd xjkA distance d betweenij，
The calculation formula comprises:
and 3c4, calculating the distance between the new class and the other classes, merging the two closest classes, and repeating the steps until all blasting events are classified as 1 if the number of the classes is still greater than 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811089429.7A CN109239775B (en) | 2018-09-18 | 2018-09-18 | Method for tracking and positioning stolen mining of mineral resources |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811089429.7A CN109239775B (en) | 2018-09-18 | 2018-09-18 | Method for tracking and positioning stolen mining of mineral resources |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109239775A CN109239775A (en) | 2019-01-18 |
| CN109239775B true CN109239775B (en) | 2020-08-14 |
Family
ID=65058200
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811089429.7A Active CN109239775B (en) | 2018-09-18 | 2018-09-18 | Method for tracking and positioning stolen mining of mineral resources |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109239775B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111009087B (en) * | 2020-01-08 | 2021-06-15 | 桂林赛普电子科技有限公司 | Mining area anti-border-crossing illegal mining monitoring and early warning method and system |
| CN112700612B (en) * | 2020-12-04 | 2022-12-13 | 徐州勇强自动化设备有限公司 | Mine anti-theft mining monitoring equipment and operation system thereof |
| CN113009558B (en) * | 2021-02-23 | 2022-06-14 | 武汉理工大学 | Micro-seismic monitoring system for preventing illegal mining of surface ore body |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102262220A (en) * | 2011-04-28 | 2011-11-30 | 中南大学 | Positioning method based on non-linear fitting micro-seismic source or acoustic emission source |
| CN106680867A (en) * | 2016-11-17 | 2017-05-17 | 大连理工大学 | Dynamic parameter method for accurate positioning of micro-seismic event |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100238765A1 (en) * | 2009-03-20 | 2010-09-23 | Grechka Vladimir | Reservoir characterization from multicomponent microseismic data |
| CN104715508B (en) * | 2015-02-13 | 2018-02-27 | 中南大学 | A kind of mesh modeling method and device |
| CN105607040A (en) * | 2015-09-07 | 2016-05-25 | 中国神华能源股份有限公司 | Mining area illegal mining prevention monitoring and positioning method and system |
| CN107478725B (en) * | 2017-08-31 | 2020-05-19 | 北京市政建设集团有限责任公司 | Method for evaluating stability of rock-sandwiched partition wall in highway tunnel |
-
2018
- 2018-09-18 CN CN201811089429.7A patent/CN109239775B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102262220A (en) * | 2011-04-28 | 2011-11-30 | 中南大学 | Positioning method based on non-linear fitting micro-seismic source or acoustic emission source |
| CN106680867A (en) * | 2016-11-17 | 2017-05-17 | 大连理工大学 | Dynamic parameter method for accurate positioning of micro-seismic event |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109239775A (en) | 2019-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Frémont et al. | High precision relative locations of earthquakes at Mount St. Helens, Washington | |
| CN109239775B (en) | Method for tracking and positioning stolen mining of mineral resources | |
| CN105510959A (en) | Method for recognizing type of vibration source of tunnel, and vibration source positioning method | |
| CN101787897A (en) | System and method for predicting coal and gas outburst risk of mine in real time | |
| CN105974469A (en) | Micro-seismic monitoring early warning analysis system and method under action of rock slope excavating perturbation | |
| CN107807381B (en) | Dynamic monitoring method and device for slope instability risk based on rock mass fracture microseism wave activity rule | |
| CN106437843B (en) | coal mine bottom plate water guide channel identification method based on microseismic monitoring | |
| CN113189644B (en) | Microseismic source positioning method and system | |
| CN104266894A (en) | Mine microearthquake signal preliminary wave moment extracting method based on correlation analysis | |
| CN116591777A (en) | Multi-field multi-source information fusion rock burst intelligent monitoring and early warning device and method | |
| CN114814939A (en) | Coal mine micro-seismic table net monitoring effect evaluation method | |
| CN105607040A (en) | Mining area illegal mining prevention monitoring and positioning method and system | |
| CN114994791B (en) | Method for evaluating monitoring capability of well-ground integrated micro-seismic monitoring system | |
| Feng et al. | A highly accurate method of locating microseismic events associated with rockburst development processes in tunnels | |
| CN115220092A (en) | Microseismic statistical method for determining advanced impact danger range of working face | |
| CN111405469B (en) | Mine earthquake monitoring system based on mobile phone mobile sensing network and crowd-sourcing positioning method | |
| CN117189239B (en) | Tunnel surrounding rock damage monitoring method | |
| CN107728218B (en) | A kind of homologous electromagnetic signal method of discrimination of coal petrography rupture | |
| CN109521221A (en) | A kind of drill bursting construction hard rock tunnel microseism wave velocity of wave real time acquiring method | |
| CN113281806B (en) | Microseismic monitoring method and system for mine border-crossing and layer-crossing mining | |
| Vinoth et al. | Slope stability monitoring by quantification and behavior of microseismic events in an opencast coal mine | |
| CN114047546A (en) | Crowd-sourcing spiral mine earthquake positioning method based on three-dimensional spatial joint arrangement of sensors | |
| CN114565712A (en) | Three-dimensional synchronous modeling and ground pressure response processing method and system based on empty area scanning | |
| Verkholantsev et al. | Seismic recorder" Ermak-5" as the part of the system of seismological monitoring of Verkhnekamskoe potash deposit | |
| CN110109176B (en) | Rock burst risk assessment method under pathological condition of tunnel microseismic sensor monitoring station network |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: Tracking and locating methods for stolen mining of mineral resources Effective date of registration: 20230323 Granted publication date: 20200814 Pledgee: Agricultural Bank of China Limited Changsha Lushan Sub branch Pledgor: CHANGSHA DIGITAL MINE CO.,LTD. Registration number: Y2023430000010 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| CP03 | Change of name, title or address |
Address after: 410205 10th floor, C1 building, Lugu enterprise Plaza, 27 Wenxuan Road, Yuelu District, Changsha City, Hunan Province Patentee after: Changsha Dimai Technology Co.,Ltd. Country or region after: China Address before: 410205 10th floor, C1 building, Lugu enterprise Plaza, 27 Wenxuan Road, Yuelu District, Changsha City, Hunan Province Patentee before: CHANGSHA DIGITAL MINE CO.,LTD. Country or region before: China |
|
| CP03 | Change of name, title or address |