CN109239775A - Mineral resources are by illegal mining tracking positioning method - Google Patents
Mineral resources are by illegal mining tracking positioning method Download PDFInfo
- Publication number
- CN109239775A CN109239775A CN201811089429.7A CN201811089429A CN109239775A CN 109239775 A CN109239775 A CN 109239775A CN 201811089429 A CN201811089429 A CN 201811089429A CN 109239775 A CN109239775 A CN 109239775A
- Authority
- CN
- China
- Prior art keywords
- explosion events
- illegal mining
- explosion
- event
- events
- 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.)
- Granted
Links
- 238000005065 mining Methods 0.000 title claims abstract description 111
- 238000000034 method Methods 0.000 title claims abstract description 33
- 229910052500 inorganic mineral Inorganic materials 0.000 title claims abstract description 14
- 239000011707 mineral Substances 0.000 title claims abstract description 14
- 238000004880 explosion Methods 0.000 claims abstract description 131
- 238000004458 analytical method Methods 0.000 claims abstract description 47
- 238000012544 monitoring process Methods 0.000 claims abstract description 41
- 230000000694 effects Effects 0.000 claims abstract description 19
- 230000033001 locomotion Effects 0.000 claims abstract description 19
- 230000002776 aggregation Effects 0.000 claims abstract description 17
- 238000004220 aggregation Methods 0.000 claims abstract description 17
- 238000004519 manufacturing process Methods 0.000 claims abstract description 16
- 238000007621 cluster analysis Methods 0.000 claims abstract description 9
- 239000011435 rock Substances 0.000 claims description 35
- 230000035939 shock Effects 0.000 claims description 13
- 238000001228 spectrum Methods 0.000 claims description 11
- 239000004744 fabric Substances 0.000 claims description 9
- 230000004044 response Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000009825 accumulation Methods 0.000 claims description 3
- 238000000151 deposition Methods 0.000 claims 1
- 238000012216 screening Methods 0.000 abstract description 3
- 238000007405 data analysis Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 5
- 241000208340 Araliaceae Species 0.000 description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 description 3
- 238000005422 blasting Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 235000008434 ginseng Nutrition 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000009412 basement excavation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000009172 bursting Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000010358 mechanical oscillation Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005180 public health Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 238000005070 sampling Methods 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
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
A kind of mineral resources of the present invention are by illegal mining tracking positioning method, including Microseismic monitoring system is installed in target area, pass through the microseismic event and explosion events in Microseismic monitoring system monitoring objective region.It selects the fixed cycle, such as daily, carry out the identification and screening of illegal mining explosion events, rejects microseismic event and mine normally produces the explosion events of generation, and then obtain illegal mining activity explosion events.Time-space serial analysis, space cluster analysis are carried out for the illegal mining explosion events drawn a circle to approve.Finally according to data analysis result, it is described in detail in conjunction with the engineering model or ore body model at scene, the three-dimensional coordinate of aggregation reference point of dated illegal mining explosion events and the spatial relationship between production operation and movable period, the motion profile of illegal mining work surface both can be obtained by multiple analysis and observation, and then real-time tracing and positioning can have been carried out to illegal mining activity.
Description
Technical field
The present invention relates to minerals resource safety fields more particularly to a kind of mineral resources by illegal mining tracking positioning method.
Background technique
Mineral resources are one of most important physical resources of national society's economic life, but there are mineral resources for various regions
Illegal illegal mining situation, cause many serious safety accidents.Mineral resources illegal mining activity is brought to the daily production in mining area
Serious threat, such as leaving a large amount of unknown goaf will cause frequent mining area Ground pressure activity, surface subsidence and roof fall accident frequency
Take place frequently life, in turn results in safety accident, the waste of resource and the destruction of environment.Due to the movable concealment of illegal mining, it is difficult to obtain
The specific region of illegal mining activity and illegal mining track.
It is therefore desirable to study a kind of technical method for grasping illegal mining active characteristics information.
Summary of the invention
(1) goal of the invention
The technical problem to be solved by the present invention is to illegal mining active characteristics information can not be grasped by illegal mining in face of mineral resources
The case where, propose a kind of method for realizing the information collection of illegal mining active characteristics and data process&analysis.
(2) technical solution
In order to achieve the above object, the main technical schemes that the present invention uses include:
A kind of mineral resources are by illegal mining tracking positioning method, comprising:
Step S1, mining Microseismic monitoring system is installed in target area, and in real-time monitoring and record target area
Rock vibration event；
Step S2, classify to the rock vibration event of record, obtain explosion events collection and illegal mining explosion events collection；
Step S3, Time-space serial analysis, space cluster analysis are carried out to illegal mining explosion events collection；
Step S4, the analysis result combining target regional model based on step S3 obtains the aggregation ginseng of illegal mining explosion events
Examination point three-dimensional coordinate, movable period and the spatial relationship between normal mining production operation；
Step S5, step S1-S4 is repeated, motion profile, fltting speed, shot cycle and the excavation of illegal mining operation are obtained
Amount is realized and carries out real-time tracing positioning and monitoring to illegal mining activity.
The step S1 further includes,
After mining Microseismic monitoring system is installed in target area, coordinate method is demarcated by fixed point big gun, repeatedly calibration system
System velocity of wave, carries out system accuracy debugging, until the Microseismic monitoring system positioning accuracy is less than 6m.
Preferably, the mining Microseismic monitoring system has the sample rate of 10KHz or more, the dynamic model greater than 115dB
Enclose, between sensor response frequency range 0.5Hz~6000Hz, timing tracking accuracy it is small less than 1 microsecond, state event location error
In 6m.
The step S2 includes:
Step S2a, signature analysis is carried out to the rock vibration event of record, obtains explosion events collection；
Step S2d, Characteristics of Source Parameters is carried out to explosion events collection and compares analysis, rejected mine and normally produce generation
Explosion events obtain illegal mining explosion events collection.
The step S2a at least one of includes the following steps:
Wave f orm analysis is carried out to each rock vibration event of record, if it exists the wave character of vibration event
Continuously to there is multiple wave crest, and S wave characteristic is unobvious, then the vibration event is denoted as explosion events；
Spectrum sigtral response is carried out to each rock vibration event of record, if it exists the frequency spectrum high frequency of vibration event
Part accounting is much larger than low frequency part, and corner frequency is greater than 500Hz, then the vibration event is denoted as explosion events；
The model analysis of cloth dragon is carried out to each rock vibration event of record, vibration event meets Bu Longmo if it exists
Type, then the vibration event information is microseismic event information, which is rejected；
Audio files is converted to the waveform of each rock vibration event of record, carries out event sound signature analysis,
The sound characteristic of vibration event meets explosion events sound characteristic if it exists, then the vibration event is denoted as explosion events.
The step S3 includes:
Step S3a, the focal shock parameter of illegal mining explosion events is obtained；
Step S3b, Time-space serial analysis is carried out to illegal mining explosion events collection based on focal shock parameter, obtains illegal mining explosion thing
The scene aggregation situation of part and motion frequency on a timeline change；
Step S3c, the distance between the blast source for calculating illegal mining explosion events two-by-two, by the hair of each explosion events
Raw position carries out clustering as clustering object, to explosion events, obtains the accumulation point of explosion events.
The step S3c includes:
Step 3c1, the location information for the illegal mining explosion events that obtaining step 2 obtains, calculate two-by-two between blast source away from
From；
Step 3c2, using the generation position of explosion events as clustering object, all blast sources constitute a class by itself, and obtain explosion
The 3 d space coordinate (x, y, z) of position, to cluster element；
Cluster element calculation formula include:
xik=(xi1,xi2,xi3),
Wherein, i=1,2,3 ... n, n are blast source number；
It step 3c3, will be apart from the smallest two blast source xikAnd xjkAnd one kind is done, calculate xikAnd xjkThe distance between
dij,
Calculation formula includes:
Step 3c4, calculate new class and remaining it is all kinds of at a distance from, then by nearest two classes merging, if class obtains number still at this time
Greater than 1, then above-mentioned steps are repeated, until all explosion events are classified as 1 class.
The step S5 includes:
Step S5a, distance threshold is set, step S1-S4 is repeated, obtains repeatedly band-like fortune of the cluster event on space-time
Dynamic direction is the movable space motion path of illegal mining；
Step S5b, by each cluster result and event period of right time, shot cycle is obtained；
Step S5c, apparent volume is added up by focal shock parameter and Modulus of Elasticity of Rock Mass parameter calculates, obtain excavated volume.
Preferably, the distance threshold is 6m.
(3) beneficial effect
The beneficial effects of the present invention are: the engineering model at scene or ore body model is combined to be described in detail, indicates and steal
The three-dimensional coordinate of aggregation reference point and the spatial relationship between production operation and the movable period for adopting explosion events, pass through
The motion profile of illegal mining work surface had both can be obtained in multiple analysis and observation, and then it is fixed tracking in real time can be carried out to illegal mining activity
Position provides for bargh's unit and governance illegal mining activity protection mineral resources and effectively instructs foundation.
Microseismic monitoring system is installed in target area, passes through the microseism thing in Microseismic monitoring system monitoring objective region
Part and explosion events.Select the fixed cycle, such as daily, carry out illegal mining explosion events identification and screening, reject microseismic event and
Mine normally produces the explosion events of generation, and then obtains illegal mining activity explosion events.For the illegal mining explosion drawn a circle to approve
Event carries out Time-space serial analysis, space cluster analysis.Finally according to data analysis result, in conjunction with scene engineering model or
Ore body model is described in detail, indicate illegal mining explosion events aggregation reference point three-dimensional coordinate and with production operation it
Between spatial relationship and the movable period, the motion profile of illegal mining work surface both can be obtained by multiple analysis and observation, in turn
Real-time tracing and positioning and quantitative statistics can be carried out to illegal mining activity, be that bargh's unit and governance illegal mining activity are protected
Shield mineral resources, which provide, effectively instructs foundation.
Detailed description of the invention
Fig. 1 is high-precision of embodiment of the present invention Microseismic monitoring system topological diagram；
Fig. 2 is vibration signal waveforms of embodiment of the present invention feature and audio files classification schematic diagram；
Fig. 3 is the aggregation situation schematic diagram of explosion events of the embodiment of the present invention；
Fig. 4 is that motion frequency variation is shown on the time shaft for the explosion events that the analysis of Time-space serial of the embodiment of the present invention obtains
It is intended to；
Fig. 5 is the cluster result schematic diagram that space cluster analysis of the embodiment of the present invention obtains illegal mining explosion events.
Specific embodiment
In order to preferably explain the present invention, in order to understand, below by specific embodiment, the present invention is made detailed
Description.
A kind of mineral resources of the present embodiment are by illegal mining tracking positioning method, and specifically, method includes the following steps:
Step S1, mining Microseismic monitoring system is installed in target area, and in real-time monitoring and record target area
Rock vibration event.
A set of mining Microseismic monitoring system, the rock vibration thing in real-time monitoring target area are installed in target area
Part, and the rock vibration event for triggering monitoring system is stored into local server.
Preferably, it is desirable that the set system has the sample rate of 10KHz or more, the dynamic range greater than 115dB, sensor
Between response frequency range 0.5Hz~6000Hz, timing tracking accuracy less than 1 microsecond, state event location error be less than 6m.
Sample rate refers to system sampling frequency, requires according to signal acquisition, and sample rate is higher than 2 times of signal maximum frequency
More than, sample rate is higher, and signal resolution ability is stronger；
Dynamic range refers to the ratio of variable signal maxima and minima, and the bigger expression system performance of dynamic range is more
It is good；
Timing tracking accuracy refers to that every equipment time synchronization error, synchronous error are the smaller the better in system；
Position error refers to the gap between system positioning coordinate position and practical focus coordinate position, position error
It is required that being less than seismic source rupture scale.
A set of satisfactory high-precision Microseismic monitoring system, system architecture diagram such as Fig. 1 institute are installed in target area
To show, sensor installation is installed by the way of bored grouting, and drilling depth requires 3-5 times that is greater than tunnel diameter, so that
It is few that sensor receives rock vibrational waveform noise, effectively ensures waveform quality；It acquires substation and configures battery, prevent disconnected
System leads to not carry out information collection and transmission due to that can not power in electric situation；In the time service of earth's surface GP configuring S or BDS
Server guarantees that time synchronization error is less than 1 microsecond between each equipment of system using network ptp time synchronization protocol.
Preferably, step S1 further includes after installing mining Microseismic monitoring system in target area, passing through fixed point big gun mark
Position fixing method, calibration system velocity of wave, carries out system accuracy debugging, until the Microseismic monitoring system positioning accurate repeatedly
Degree is less than 6m.
System is system accuracy according still further to carry out system debug, the performance indicator mainly debugged after finishing, and is led to
Fixed point big gun calibration coordinate method is crossed, repeatedly calibration system velocity of wave, finally obtaining one can accurately reflect monitoring region rock mass characteristics
Elastic wave propagation speed, and guarantee state event location precision be less than 6m.
Step S2, classify to the rock vibration event of record, obtain explosion events collection and illegal mining explosion events collection；
Wave f orm analysis, spectrum sigtral response, cloth imperial (Brune) model analysis, event are carried out for trigger event
Sound characteristic analysis has microseismic event, mechanical oscillation event and explosion events so that trigger event be classified, and collects mine
Production resettlement record in mountain rejects the explosion thing that mine normally produces generation by comparing the Characteristics of Source Parameters of explosion events
Part, to obtain the explosion events of improper production.
Explosion events refer to that fractured rock generates strong rock vibration in turn by the way of explosive blasting for exploitated ore
It moves to be determined as explosion events by Microseismic monitoring system；
Microseismic event refers to rock in a kind of breakoff phenomenon for being squeezed with generating under shearing force, this rock
Breakoff phenomenon can also generate elastic wave and be detected by Microseismic monitoring system, but its wave character and Blast waveform feature exist
Significant otherness；
Illegal mining explosion events refer to the explosion events in addition to explosion events are normally exploited in mine.
Step S2a, signature analysis is carried out to the rock vibration event of record, obtains explosion events collection.
The recognition methods of explosion events and microseismic event has: being used by the event to Microseismic monitoring system trigger collection
Wave character identification, spectrum analysis, imperial (Brune) model analysis of cloth, event sound signature analysis identify explosion events class
Type.
The step S2a at least one of includes the following steps:
Wave f orm analysis is carried out to each rock vibration event of record, if it exists the wave character of vibration event
Continuously to there is multiple wave crest, and S wave characteristic is unobvious, then the vibration event is denoted as explosion events；
It, can be to the rock of record by the rock vibration event to record by wave f orm analysis such as P wave and S wave analysis
Body vibration event is classified, and Blast waveform feature has and multiple wave crest continuously occurs, and Blast waveform S wave characteristic is unknown
It is aobvious.
Spectrum sigtral response is carried out to each rock vibration event of record, if it exists the frequency spectrum high frequency of vibration event
Part accounting is much larger than low frequency part, and corner frequency is greater than 500Hz, then the vibration event is denoted as explosion events；
By carrying out Fourier's variation to event waveforms, explosion events spectrum signature, explosion events frequency spectrum high frequency can get
Part accounting is more, and corner frequency is high, generally higher than 500Hz.
The model analysis of cloth dragon is carried out to each rock vibration event of record, vibration event meets Bu Longmo if it exists
Type, then the vibration event information is microseismic event information, which is rejected；
By judging according to whether spectrum curve meets cloth dragon model, the spectrum curve of general mine microquake event is accorded with
Imperial (Brune) model of cloth is closed, explosion events do not meet the model；
Cloth dragon (Brune) model is the Source Model proposed in the 1970's, which destroys broken to shear
Bad, entire destroyed area has continuous sliding velocity, has apparent variation characteristic in the spectrum analysis of waveform, this
It is characterized in for the waveform of non-microseismic event.
Audio files is converted to the waveform of each rock vibration event of record, carries out event sound signature analysis,
The sound characteristic of vibration event meets explosion events sound characteristic if it exists, then the vibration event is denoted as explosion events；
System monitoring to event can be converted to audio files by waveform, explosion events are distinguished by sound characteristic.
Event sound signature analysis identifies explosion events type, as shown in Figure 2.
Step S2d, Characteristics of Source Parameters is carried out to explosion events collection and compares analysis, rejected mine and normally produce generation
Explosion events obtain illegal mining explosion events collection.
Improper operation explosion events are determined according to explosion events focal shock parameter, to finally identify that illegal mining operation is quick-fried
Broken event as shown in figure 3, obtaining monitoring data from mine installation Microseismic monitoring system, and carries out identification to monitoring data and picks
It removes, draws a circle to approve the illegal mining explosion events in this period.
Alternatively, mining production generally has the fixed explosion period, by rejecting in the period by time piecewise analysis
Explosion events, then remaining explosion events are regarded as illegal mining explosion events.
It selects the fixed cycle, such as daily, carry out the identification and screening of illegal mining explosion events, rejects microseismic event and mine
The explosion events that normal production generates, and then obtain illegal mining activity explosion events.
Step S3, Time-space serial analysis, space cluster analysis are carried out to illegal mining explosion events collection.
Improper production resettlement event is the illegal mining bursting work activity monitored in region, for what is drawn a circle to approve
Illegal mining explosion events carry out Characteristics of Source Parameters analysis, Time-space serial analysis, space cluster analysis.
Step S3a, the focal shock parameter of illegal mining explosion events is obtained.
Explosion events focal shock parameter refers to that explosion events source physics feature describes, general to be retouched with " when, sky, by force "
State a shock event.
Step S3b, Time-space serial analysis is carried out to illegal mining explosion events collection based on focal shock parameter, obtains illegal mining explosion thing
The scene aggregation situation of part and motion frequency on a timeline change.
Illegal mining explosion events carry out Time-space serial analysis and refer to the scene aggregation situation of explosion events and in time shaft
On motion frequency variation.
Step S3c, the distance between the blast source for calculating illegal mining explosion events two-by-two, by the hair of each explosion events
Raw position carries out clustering as clustering object, to explosion events, obtains the accumulation point of explosion events.
Space cluster analysis is one of main means of Spatial Data Mining and Knowledge Discovery, is widely used to geography
The numerous areas such as, geology, meteorology, cartography, astronomy and public health.Use for reference earthquake spatial distribution often
The characteristics of non-spherical aggregation is presented, the correlation experience of application space clustering, for mine blasting event data, using poly-
Class method handles the irregular aggregation shape of explosion.
The step S3b includes:
Step 3c1, the location information for the illegal mining explosion events that obtaining step 2 obtains, calculate two-by-two between blast source away from
From.
The principle of space cluster analysis method is to calculate the distance between blast source two-by-two, by comparing, selection distance symbol
Close the threshold value of setting.
Step 3c2, using the generation position of explosion events as clustering object, all blast sources constitute a class by itself, and obtain explosion
The 3 d space coordinate (x, y, z) of position, to cluster element.
Using the generation position of explosion events as clustering object, cluster 3 d space coordinate that element is explosion position (x,
Y, z), following formula calculates:
xik=(xi1,xi2,xi3) (i=1,2,3 ... n)
It step 3c3, will be apart from the smallest two blast source xikAnd xjkAnd one kind is done, calculate xikAnd xjkThe distance between
dij。
When starting to calculate, each blast source constitutes a class by itself (event number n then shares n class at this time).Calculate each blast source it
Between distance, measure their close and distant degree with the distance between they, by two nearest blast sources and do one kind.It adopts
Use dijTo indicate xikAnd xjkThe distance between, calculation formula is as follows:
Step 3c4, calculate new class and remaining it is all kinds of at a distance from, then by nearest two classes merging, if class obtains number still at this time
Greater than 1, then above-mentioned steps are repeated, until all explosion events are classified as 1 class.
Clustered using knearest neighbour method, that is, calculate new class with remaining it is all kinds of at a distance from, then nearest two classes are closed
And if class obtains number and is still greater than 1 at this time, above-mentioned steps are repeated, until all explosion events are classified as 1 class.Wherein class
The event of Gq and class Gp, drkDistance definition are as follows:
drk=min { dpk,dqk}
By monitoring that explosion events carry out clustering to Microseismic monitoring system, its aggregation place is found out.
Characteristics of Source Parameters analysis is carried out for the illegal mining explosion events drawn a circle to approve, Time-space serial is analyzed, space is poly-
Alanysis, as shown in Figure 4 and Figure 5.66 illegal mining explosion locating events are shared in the analytical cycle, these events are mainly distributed
In squat lower wall, gulf lower wall and two sides.It is main to concentrate near driving middle section.From the time of generation, morning 2:00~
4:00 point, morning 10:00~15:00 point, evening 20:00~21:00 point or so.
The illegal mining explosion events space distribution situation out of this analytical cycle is it is found that illegal mining clusters of events degree one on the whole
As, consider that aggregation scale there will be certain redundancy relative to systematic error, selects with 6m (empirical value) as cluster cell ruler
Degree, parameter threshold are 2 (empirical values)；
Step S4, the analysis result combining target regional model based on step S3 obtains the aggregation ginseng of illegal mining explosion events
Examination point three-dimensional coordinate, movable period and the spatial relationship between normal mining production operation.
Result is analyzed based on step 3 and the engineering model at scene or ore body model is combined to be described in detail, and is indicated and is stolen
The three-dimensional coordinate of aggregation reference point and the spatial relationship between production operation and the movable period of explosion events are adopted, then
The movable work surface position of illegal mining is grasped.
It is described in detail in conjunction with the engineering model or ore body model at scene, indicates the aggregation ginseng of illegal mining explosion events
It is the three-dimensional coordinate of examination point and the spatial relationship between production operation and movable period, as shown in the table:
As a specific example of step S4, according to the analysis of step 3 as a result, combining the engineering model or mine at scene
Body Model is described in detail, and indicates the three-dimensional coordinate of the aggregation reference point of illegal mining explosion events and between production operation
Spatial relationship and the movable period, it is as shown in the table:
Step S5, step S1-S4 is repeated, motion profile, fltting speed, shot cycle and the excavation of illegal mining operation are obtained
Amount is realized and carries out real-time tracing positioning and monitoring to illegal mining activity.
The movement rail that illegal mining work surface is pushed ahead can be obtained by the way that step S1, step S2, step S3 process is repeated several times
Mark can carry out real-time tracing to the movable fltting speed of illegal mining, shot cycle, blasting energy by quantization statistical analysis
Positioning and monitoring provide for illegal mining activity improvement and effectively reliably instruct foundation.
Step S5a, distance threshold is set, step S1-S4 is repeated, obtains repeatedly band-like fortune of the cluster event on space-time
Dynamic direction is the movable space motion path of illegal mining.
By the way that reasonable distance threshold is arranged, it is preferable that the distance threshold is 6m, i.e. the focus of illegal mining explosion is maximum
Rupture size will cluster each time by repeated multiple times progress illegal mining explosion events clustering and be considered that illegal mining explosion is primary
Propulsion activity is tunneled, multiple cluster event will appear the band-like direction of motion on space-time, and zonal motion direction is living for illegal mining
Dynamic space motion path.
Step S5b, by each cluster result and event period of right time, shot cycle is obtained.
Shot cycle is determined by each cluster result and event period of right time.
Step S5c, apparent volume is added up by focal shock parameter and Modulus of Elasticity of Rock Mass parameter calculates, obtain excavated volume.
According to explosion hypocenter released energy calculate illegal mining excavated volume, excavated volume can be added up by focal shock parameter apparent volume and
Modulus of Elasticity of Rock Mass parameter is calculated.
It is difficult to accurately carry out illegal mining activity trajectory description when monitoring data deficiency, needs repeatedly to be monitored number for a long time
According to analysis.
In order to which better tracing and positioning illegal mining activity needs to be equipped with high-precision Microseismic monitoring system and detailed mine is quick-fried
The data that break a record are carried out referring to analysis.
Finally, bargh can carry out site inspection confirmation according to illegal mining moving position, while reinforcing the production in the region
The safety of safety management, support personnel and equipment.
It should be clear that the invention is not limited to specific configuration described above and shown in figure and processing.
For brevity, it is omitted here the detailed description to known method.In the above-described embodiments, it has been described and illustrated several
Specific step is as example.But method process of the invention is not limited to described and illustrated specific steps, ability
The technical staff in domain can be variously modified, modification and addition, or change the step it after understanding spirit of the invention
Between sequence.
Finally, it should be noted that above-described embodiments are merely to illustrate the technical scheme, rather than to it
Limitation；Although the present invention is described in detail referring to the foregoing embodiments, those skilled in the art should manage
Solution: it can still modify to technical solution documented by previous embodiment, or special to part of or whole technologies
Sign is equivalently replaced；And these modifications or substitutions, various embodiments of the present invention skill that it does not separate the essence of the corresponding technical solution
The range of art scheme.
Claims (9)
1. a kind of mineral resources are by illegal mining tracking positioning method, which is characterized in that the described method includes:
Step S1, mining Microseismic monitoring system, and the rock mass in real-time monitoring and record target area are installed in target area
Vibration event；
Step S2, classify to the rock vibration event of record, obtain explosion events collection and illegal mining explosion events collection；
Step S3, Time-space serial analysis, space cluster analysis are carried out to illegal mining explosion events collection；
Step S4, the analysis result combining target regional model based on step S3 obtains the aggregation reference point of illegal mining explosion events
Three-dimensional coordinate, movable period and the spatial relationship between normal mining production operation；
Step S5, step S1-S4 is repeated, motion profile, fltting speed, shot cycle and the excavated volume of illegal mining operation are obtained, it is real
Real-time tracing positioning and monitoring now are carried out to illegal mining activity.
2. the method according to claim 1, wherein the step S1 further includes,
After mining Microseismic monitoring system is installed in target area, coordinate method is demarcated by fixed point big gun, repeatedly calibration system wave
Speed carries out system accuracy debugging, until the Microseismic monitoring system positioning accuracy is less than 6m.
3. the method according to claim 1, wherein
The mining Microseismic monitoring system has sample rate, the dynamic range greater than 115dB, the sensor response of 10KHz or more
Between frequency range 0.5Hz~6000Hz, timing tracking accuracy less than 1 microsecond, state event location error be less than 6m.
4. the method according to claim 1, wherein the step S2 includes:
Step S2a, signature analysis is carried out to the rock vibration event of record, obtains explosion events collection；
Step S2d, Characteristics of Source Parameters is carried out to explosion events collection and compares analysis, reject the explosion thing that mine normally produces generation
Part obtains illegal mining explosion events collection.
5. according to the method described in claim 3, it is characterized in that, the step S2a at least one of includes the following steps:
Wave f orm analysis is carried out to each rock vibration event of record, the wave character of vibration event is continuous if it exists
There is multiple wave crest, and S wave characteristic is unobvious, then the vibration event is denoted as explosion events；
Spectrum sigtral response is carried out to each rock vibration event of record, the frequency spectrum high frequency section of vibration event accounts for if it exists
Than being much larger than low frequency part, corner frequency is greater than 500Hz, then the vibration event is denoted as explosion events；
The model analysis of cloth dragon is carried out to each rock vibration event of record, vibration event meets cloth dragon model if it exists, then
The vibration event information is microseismic event information, which is rejected；
Audio files is converted to the waveform of each rock vibration event of record, event sound signature analysis is carried out, if depositing
Meet explosion events sound characteristic in the sound characteristic of vibration event, then the vibration event is denoted as explosion events.
6. the method according to claim 1, wherein the step S3 includes:
Step S3a, the focal shock parameter of illegal mining explosion events is obtained；
Step S3b, Time-space serial analysis is carried out to illegal mining explosion events collection based on focal shock parameter, obtains the hair of illegal mining explosion events
Radix Rehmanniae point assembles situation and motion frequency on a timeline changes；
Step S3c, the distance between the blast source for calculating illegal mining explosion events two-by-two, by the generation position of each explosion events
As clustering object, clustering is carried out to explosion events, obtains the accumulation point of explosion events.
7. the method according to claim 1, wherein the step S3c includes:
Step 3c1, the location information for the illegal mining explosion events that obtaining step 2 obtains calculates the distance between blast source two-by-two；
Step 3c2, using the generation position of explosion events as clustering object, all blast sources constitute a class by itself, and obtain explosion position
3 d space coordinate (x, y, z), for cluster element；
Cluster element calculation formula include:
xik=(xi1,xi2,xi3),
Wherein, i=1,2,3 ... n, n are blast source number；
It step 3c3, will be apart from the smallest two blast source xikAnd xjkAnd one kind is done, calculate xikAnd xjkThe distance between dij,
Calculation formula includes:
Step 3c4, calculate new class and remaining it is all kinds of at a distance from, then by nearest two classes merging, if class obtains number and is still greater than at this time
1, then above-mentioned steps are repeated, until all explosion events are classified as 1 class.
8. the method according to claim 1, wherein the step S5 includes:
Step S5a, distance threshold is set, step S1-S4 is repeated, obtains repeatedly zonal motion side of the cluster event on space-time
To for the movable space motion path of illegal mining；
Step S5b, by each cluster result and event period of right time, shot cycle is obtained；
Step S5c, apparent volume is added up by focal shock parameter and Modulus of Elasticity of Rock Mass parameter calculates, obtain excavated volume.
9. according to the method described in claim 8, it is characterized in that,
The distance threshold is 6m.
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 true CN109239775A (en) | 2019-01-18 |
| CN109239775B 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) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111009087A (en) * | 2020-01-08 | 2020-04-14 | 桂林赛普电子科技有限公司 | Mining area anti-border-crossing illegal mining monitoring and early warning method and system |
| CN112700612A (en) * | 2020-12-04 | 2021-04-23 | 徐州勇强自动化设备有限公司 | Mine anti-theft mining monitoring equipment and operation system thereof |
| CN113009558A (en) * | 2021-02-23 | 2021-06-22 | 武汉理工大学 | Micro-seismic monitoring system for preventing illegal mining of surface ore body |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100238765A1 (en) * | 2009-03-20 | 2010-09-23 | Grechka Vladimir | Reservoir characterization from multicomponent microseismic data |
| CN102262220A (en) * | 2011-04-28 | 2011-11-30 | 中南大学 | Positioning method based on non-linear fitting micro-seismic source or acoustic emission source |
| CN104715508A (en) * | 2015-02-13 | 2015-06-17 | 中南大学 | Grid modeling method and device |
| CN105607040A (en) * | 2015-09-07 | 2016-05-25 | 中国神华能源股份有限公司 | Mining area illegal mining prevention monitoring and positioning method and system |
| CN106680867A (en) * | 2016-11-17 | 2017-05-17 | 大连理工大学 | Dynamic parameter method for accurate positioning of micro-seismic event |
| CN107478725A (en) * | 2017-08-31 | 2017-12-15 | 北京市政建设集团有限责任公司 | Rock partition wall method for estimating stability is pressed from both sides in a kind of vcehicular tunnel |
-
2018
- 2018-09-18 CN CN201811089429.7A patent/CN109239775B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100238765A1 (en) * | 2009-03-20 | 2010-09-23 | Grechka Vladimir | Reservoir characterization from multicomponent microseismic data |
| CN102262220A (en) * | 2011-04-28 | 2011-11-30 | 中南大学 | Positioning method based on non-linear fitting micro-seismic source or acoustic emission source |
| CN104715508A (en) * | 2015-02-13 | 2015-06-17 | 中南大学 | Grid modeling method and device |
| CN105607040A (en) * | 2015-09-07 | 2016-05-25 | 中国神华能源股份有限公司 | Mining area illegal mining prevention monitoring and positioning method and system |
| CN106680867A (en) * | 2016-11-17 | 2017-05-17 | 大连理工大学 | Dynamic parameter method for accurate positioning of micro-seismic event |
| CN107478725A (en) * | 2017-08-31 | 2017-12-15 | 北京市政建设集团有限责任公司 | Rock partition wall method for estimating stability is pressed from both sides in a kind of vcehicular tunnel |
Non-Patent Citations (1)
| Title |
|---|
| 罗一忠: "大面积采空区失稳的重大危险源辨识", 《中国优秀博硕士学位论文全文数据库 (博士)工程科技Ⅰ辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111009087A (en) * | 2020-01-08 | 2020-04-14 | 桂林赛普电子科技有限公司 | Mining area anti-border-crossing illegal mining monitoring and early warning method and system |
| CN112700612A (en) * | 2020-12-04 | 2021-04-23 | 徐州勇强自动化设备有限公司 | Mine anti-theft mining monitoring equipment and operation system thereof |
| CN113009558A (en) * | 2021-02-23 | 2021-06-22 | 武汉理工大学 | Micro-seismic monitoring system for preventing illegal mining of surface ore body |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109239775B (en) | 2020-08-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Xu et al. | Microseismic monitoring and stability evaluation for the large scale underground caverns at the Houziyan hydropower station in Southwest China | |
| Feng et al. | Sectional velocity model for microseismic source location in tunnels | |
| CN109239775A (en) | Mineral resources are by illegal mining tracking positioning method | |
| CN103544791B (en) | Based on the underground system for monitoring intrusion of seismic event | |
| CN103306722B (en) | Micro-seismic multi-dimensional information comprehensive region detection and evaluation method for impact danger region | |
| CN107727737A (en) | One kind excavates the unknown geological structure activation recognition methods of underground engineering under Unloading Effect | |
| CN102913285B (en) | Mine roof hazard warning method | |
| CN109595036A (en) | A kind of method for early warning of Mine Ground Pressure Disaster | |
| CN105807312A (en) | Coal mine roof rock mass vertical zone determination method based on micro-seismic monitoring | |
| CN106199700A (en) | A kind of underground water seal oil storage cave depot micro seismic monitoring method and system | |
| CN103953392A (en) | Method for distinguishing position of rockburst risk on deep buried tunnel section | |
| CN103336297A (en) | Micro-rupture vector scanning method | |
| CN106443761A (en) | Full-frequency band and wide-frequency domain shake monitoring system for mine earthquakes and rock bursts | |
| Feng et al. | A highly accurate method of locating microseismic events associated with rockburst development processes in tunnels | |
| CN105607040A (en) | Mining area illegal mining prevention monitoring and positioning method and system | |
| Occhiena et al. | A microseismic-based procedure for the detection of rock slope instabilities | |
| CN114814939A (en) | Coal mine micro-seismic table net monitoring effect evaluation method | |
| CN113189644A (en) | Microseismic source positioning method and system | |
| CN110376643B (en) | Micro-seismic effect data processing method for jet grouting pile diameter detection | |
| Eberhardt et al. | Landslide monitoring: The role of investigative monitoring to improve understanding and early warning of failure | |
| CN111405469B (en) | Mine earthquake monitoring system based on mobile phone mobile sensing network and crowd-sourcing positioning method | |
| CN103941281A (en) | Distributed mine earthquake detection method and device | |
| CN107515419A (en) | The method of estimation and device of rock stability | |
| CN112230275A (en) | Seismic waveform identification method and device and electronic equipment | |
| CN107728218B (en) | A kind of homologous electromagnetic signal method of discrimination of coal petrography rupture |
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 | ||
| 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 |
|
| CP03 | Change of name, title or address | ||
| 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 |