CN101770038A - Intelligent positioning method of mine microquake sources - Google Patents

Intelligent positioning method of mine microquake sources Download PDF

Info

Publication number
CN101770038A
CN101770038A CN201010100527A CN201010100527A CN101770038A CN 101770038 A CN101770038 A CN 101770038A CN 201010100527 A CN201010100527 A CN 201010100527A CN 201010100527 A CN201010100527 A CN 201010100527A CN 101770038 A CN101770038 A CN 101770038A
Authority
CN
China
Prior art keywords
microquake
microquake sources
monitoring
sensor
mine
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
Application number
CN201010100527A
Other languages
Chinese (zh)
Other versions
CN101770038B (en
Inventor
陈炳瑞
冯夏庭
徐速超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan Institute of Rock and Soil Mechanics of CAS
Original Assignee
Wuhan Institute of Rock and Soil Mechanics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan Institute of Rock and Soil Mechanics of CAS filed Critical Wuhan Institute of Rock and Soil Mechanics of CAS
Priority to CN2010101005273A priority Critical patent/CN101770038B/en
Publication of CN101770038A publication Critical patent/CN101770038A/en
Application granted granted Critical
Publication of CN101770038B publication Critical patent/CN101770038B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses an intelligent positioning method of mine microquake sources, which is a microquake resource intelligent positioning system which comprehensively considers various influencing factors and targets, and integrates with sensors distribution-noise wave filtration-microquake source positioning analysis-three-dimensional display in to a whole. The method realizes the real-time analysis and the forecast of the incubation, the development and the occurrence of mine underground geological disaster through a self-developed procedure, solves the defects that the traditional sensors distribution method is not systematic, the noise wave filtration is incomplete, the microquake wave speed model is not exact and the positioning method is easily diffused, and has the characteristics of friendly operation interface, good noise wave filtration performance, exact, convenient and fast microquake source positioning analysis, visual and pictorial result display, and wide application range. The method has wide application value in the fields of mineral engineering, water conservancy and electric power engineering, petroleum engineering, soil engineering and underground engineering and the like.

Description

Intelligent positioning method of mine microquake sources
Technical field
The invention belongs to the mine microquake monitoring technical field, more specifically relate to a kind of intelligent positioning method of mine microquake sources, this method can be widely used in Mineral Engineering, Hydraulic and Hydro-Power Engineering, petroleum engineering, Geotechnical Engineering and underground works.
Background technology
Along with the rapid growth of Chinese national economy, a main developing direction of deep mining having become China mining.The increasing of the mining degree of depth, the corresponding increase of terrestrial stress, the deterioration of mining conditions, phenomenon showed increased such as roof fall, wall caving appear when causing exploiting, even power ground such as rock burst occurs and press disaster, cause huge lives and properties and economic loss for country, bargh and the people, seriously restricted the sustainable development of national economy and bargh.Statistics show that along with the increase of mining depth, the dynamic disaster that is popular in mines such as mountain copper mine, the plumbous zinc ore of all mouths, Tongling Shizishan Copper Mine, the exquisite gold mine in Shandong, Xiangxi Gold Mine, ore deposit, wax gourd mountain obviously is increase trend, and supporting is overhauled engineering and also greatly increased; Compare China Mine Geological complex structure, mine accident pilosity with the main coal producer in the world, average 1,000,000 tons of mortality ratio are about 2-6, being about 100 times and about 30 times of South Africa of the U.S., about year death toll 6000 people, is 3 times of other producing coal country death toll summations of the whole world.Therefore, the monitoring of rock pressure and down-hole geologic hazard has become the difficult problem that China's mine safety needs to be resolved hurrily with forecast.Studies show that, no matter dynamic disasters such as mountain, right and wrong colliery rock burst, ore deposit shake, still the coal in colliery is given prominence to disasters such as (or gushing out) and mine floor gushing water with gas, all is the result of rock rupture process unstabilitys such as the micro rupture germinating brought out of the stress field disturbance in the mining process, development, perforation.Monitoring shows no matter be which kind of mine power disaster, in most cases, before dynamic disaster occurs, micro rupture (microseismic activity) omen is arranged all.Therefore, the monitoring microseismic activity, hazard forecasting is extremely important for mining geology to obtain the precursor information that disaster arrives.
For monitoring tool and equipment, carried out going deep into systematic research both at home and abroad, and for the research of microseismic monitoring sensor method for arranging, many indexs of noise technology for eliminating, microquake sources high-precision intelligent location technology, still need and will further improve, existing both at home and abroad method mainly has the following disadvantages:
1. sensor arrangement method
Transducer arrangements not only influences the monitoring of microseismic signals, and the uniqueness of locating speed, precision and the positioning result of different microseism location algorithms is also had in various degree influence.Reasonably the transducer arrangements scheme not only can monitor more effective microseismic signals more broadly, and can make location algorithm determine source location and origin time fast and accurately.It mainly is rule of thumb that moment sensor is arranged, the sensor monitors that the different staff of experience arranges to microseismic signals often differ greatly, often can not make sensor farthest monitor effective microseismic signals, be difficult to also guarantee that sensor forms an optimum array, cause microseism locating speed and precision to be subjected in various degree influence, the ore deposit shakes accurate prediction and also is restricted to a great extent.Therefore, be necessary very much transducer arrangements is optimized.
2. noise filtering technique
On-the-spot micro seismic monitoring, can noise eliminating be very important technical measures, be one of the determinative that successfully predict of rock burst.The present method of canceling noise single indexs (such as threshold value, average frequency or ring number etc.) that adopt more, but for many noises be interweaved (Construction traffic, roofbolter, air compressor machine, TBM operating noise etc.), and under the comparatively approaching condition of the actual fracturing features of noise signal feature and country rock, single index is the influence of elimination environmental noise preferably often.Therefore, be necessary to study many indexs filtering technique, further improve the recognition capability of useful signal.
3. microquake sources method for positioning analyzing
Microquake sources localization method forefathers have done big quantity research, in the practical application, and according to unusual two classes that are divided into of the parameter that participates in finding the solution, the one, the known speed model is found the solution the classical localization method of origin time and microquake sources position; The one, the combination method localization method that microquake sources position, origin time and rate pattern are found the solution together.The former, in seismic field, mining engineering, be most widely used, rate pattern is the maximum not enough of this method to forbidding, though forefathers have done many researchs to rate pattern, but because rock material is complicated, heterogeneous, contain a large amount of cracks, joint and little discontinuity surface, and position, size and the trend of prior very difficult definite these bad bodies, also be difficult to draw a clear the boundary between these bad bodies, reasonable given velocity of wave model remains difficulty in advance, and this has influenced the stability and the bearing accuracy of location algorithm to a great extent; The latter has solved rate pattern preferably and has given inaccurate problem, has improved the microquake sources bearing accuracy on largely, but microquake sources position, origin time and these parameters of medium velocity interrelated, has brought the positioning result problem of unstable again.These location algorithms mainly adopt least square method to find the solution in addition, and this class methods disadvantage is easily to disperse in the solution procedure.Although proposed various improving one's methods in order to improve the stability of solution scholars, as decomposition of singular matrix method, damped least square method etc., various improved methods also all belong to the category of linear location, have always solved this problem and have brought that problem.Therefore, be necessary to look for another way and explore new location algorithm and method for solving.
For this reason, gathering microseismic signals, when carrying out the microseism location, must carry out reasonably optimizing, making it monitor more effective information as far as possible, forming an optimum sensor array sensing station; Rate pattern is carried out correct identification, make its energy correct calculation travel-time of microseismic signals in medium; Many noises interleaved signal is carried out filtering, make it at utmost obtain useful signal; The focus localization method is improved, avoid the wild effect of separating in the microseism location as far as possible, locate position and time that microseism takes place fast and accurately, the result of conveniently omnibearing display analysis and forecast is for the prediction of mining geology disaster provides information more accurately and reliably.Microseism intelligence method for positioning analyzing proposes in order to address these problems, this method not only can monitor more microseismic signals, effectively eliminate the influence of environmental noise, accurately locate the position of mine microquake sources fast, but also can pass through rock interior Cracks Evolution real-time positioning research Cracks Evolution rule and mechanism, be the important means of indoor (or on-the-spot) experimental study rock failure mechanism of rock rule and pattern.
Summary of the invention
The objective of the invention is to be to provide a kind of microquake sources intelligent locating method, solve the problems referred to above and deficiency that existing method exists, further improve microquake sources space-time bearing accuracy, improved the mine disaster prediction precision.This method has realized that transducer arrangements Automatic Optimal, down-hole disaster preparation process real-time analysis, disaster occurrence positions and time high accuracy prediction and transducer arrangements scheme, monitoring analysis and forecast result's dynamic 3 D intuitively shows.The real-time monitoring that this method can be used for not only that the mine is underground, overall process is bred, develops, takes place in geologic hazard in the strip mining transformation process, also can be used for indoor (or on-the-spot) experimental study of hard rock failure mechanism in the loading procedure, failure mode and damage evolution rule, also be with a wide range of applications for fields such as Tunnel Engineering, Hydraulic and Hydro-Power Engineering, exploitation of oil-gas field, nuclear waste disposals simultaneously.
In order to realize above-mentioned purpose, the present invention adopts following technical measures:
1) based on numerical analysis and particle group optimizing method (Particle Swarm Optimization, PSO) microquake sources monitoring sensor layout optimization technology; 2) many indexs of microseismic signals intelligent filter technology; 3) microquake sources layering intelligence location technology; 4) microseismic analysis and the dynamic three-dimensional display technology that predicts the outcome.
A kind of microquake sources intelligent locating method the steps include:
A, based on numerical analysis and PSO optimization method microquake sources monitoring sensor layout optimization technology:
1, at first, according to the actual mining technology in mine (the foundation block chambering method is for example arranged), (for example horizontal tectonics stress is main to geologic condition, rock mass fracture growth), set up numerical simulator, the movable advanced line number value of mining is analyzed, held the zone that the mining geology disaster might occur in the recovery process on the whole, as the key area of microseismic monitoring sensor layout optimization.
2, in the recovery activity scope that Mine Safety in Production is concerned about,, determine N (for example N=100) microquake sources (three-dimensional coordinate and origin time) in conjunction with the numerical analysis result.
3, according to the distribution of heterogeneity ore deposit, mine rock wave velocity test, goaf wave velocity test and heterogeneity rock stratum, determine microseism velocity of wave propagation model and span.
4, then,, determine the scope of transducer arrangements, and in scope, utilize mixed congruence method to produce X group (for example X=16) sensing station (three-dimensional coordinate) at random, every group of M (for example M=30) sensing station according to the actual mining technology in mine.
5, to every group of sensor, judge that the focus microseismic signals propagates into tomography and the goaf that sensor experiences, and calculate the distance that microseismic signals is propagated therein, and select the proper speed model automatically, sensor monitors is then to calculate N * M (100 * 30=3000) according to formula (a).
t ij = T i + L ij V - - - ( a )
Wherein, t IjBe the time that j wave detector receives the signal that i microquake sources send, T iBe i the microquake sources origin time of earthquake, V is the equivalent velocity of wave that ripple is propagated in medium, L IjBe the distance of i microquake sources to j wave detector, computing formula is as follows:
L ij = ( x j - x i ) 2 + ( y i - y j ) 2 + ( z j - z i ) 2
In the formula, (x i, y i, z i) an i microquake sources position coordinates, (x j, y j, z j) be k wave detector position coordinates.
6, when calculating by formula (b) carry out random perturbation, Virtual Monitoring is then to obtain N * M (100 * 30=3000).
t V=(1+3c(-1) xa 0)t c (b)
Wherein, t cFor Virtual Monitoring then, t cWhen monitoring for calculating, a 0Be the number between the 0-0.05, c is the random number between 0-1, and x=(int) 3c rounds 3 times random numbers.
7, considering that transducer arrangements cost and optical cable cabling are easily under the prerequisite, so that when monitoring and the accumulation residual sum of squares (RSS) when calculating is minimum is objective function, utilize population (Particle SwarmOptimization, PSO) colony intelligence method, according to M (M=30) sensor monitors signal, N (N=100) microquake sources positioned.
8, if bearing accuracy meets the demands, and sensor can monitor microseismic signals preferably, and the position is reasonable, finishes transducer arrangements optimization; Otherwise, carry out next step.
9, if do not find the right sensors position, then utilize the PSO operation, in the scope that sensor can be arranged, produce X group (X=16) new sensing station, return (4) step of this section, new one group sensing station is carried out quality judge.
Whole optimizing process as shown in Figure 1.
B, many indexs of microseismic signals intelligent filter technology:
At first set up various types of noise databases, noise signal feature (as rising time, tale, peak counting, energy, amplitude duration etc.) is carried out analytic induction, make up neural network learning and test sample book by test; Then, utilize genetic algorithm (Genetic algorithm, GA) structure, weights and the threshold value of optimization neural network, and utilize superpower non-linear mapping capability of artificial neural network and self-learning capability neural network training, set up the mapping relations between different noise signal features and the noise types; At last, Input Monitor Connector to microseismic signals it is carried out filtering, obtain effective microseismic signals.Neural network filter construction synoptic diagram as shown in Figure 2.
C, microquake sources intelligence layering method for positioning analyzing:
1, at first utilize the intelligent filter technology that sensor array is monitored the processing such as dry, filtering that disappear of many noise signals, farthest the accurate recording microseismic signals is then chosen and is used for the microseism ripple type of focus location.
2, initialization PSO parameter, microquake sources scope and velocity of wave scope.
3, utilize mixing congruence random algorithm initialization population position (source location and velocity of wave) and particle flying speed (iteration step length).
4, according to time difference positioning principle, calculate the adaptive value of particle by formula (c), and judge whether to satisfy flight number of times and the bearing accuracy of setting in advance, satisfy, carry out (6) step of this section; Otherwise, carry out (5) step of this section;
f = Σ k = 1 n [ T k M - T k C ] 2 - - - ( c )
In the formula, T k MAnd T k CBe respectively k sensor monitors and when calculating.
5, according to formula (d) and (e) the more position of new particle (microquake sources) and flying speed (iteration step length), returned for (3) step;
V id=wV id+c 1r 1(P id-X id)+c 2r 2(P gd-X id) (d)
X id=X id+V id (e)
Wherein, w is an inertia weight; c 1And c 2Be the study factor of non-negative constant; r 1And r 2Be the random number between [0,1]; D=1,2 ..., D; With
Figure GSA00000005848600061
Figure GSA00000005848600062
Be respectively the optimum focal shock parameter that i particle position, the optimum focal shock parameter that searches up to now and whole population search up to now.
6, microquake sources coordinate and the rate pattern substitution formula (f) that recognizes obtained microquake sources origin time t;
t = Σ k = 1 n ( w k - L k V ) n - - - ( f )
7, judge that microseismic signals propagates into sensor and whether passes big tomography, dead zone, reject this signal, return (2) step of this section and reorientate if pass; Otherwise it is that correct locating with microquake sources is accurately that microseismic signals is chosen, and finishes the microquake sources location.
Whole microquake sources positioning analysis process as shown in Figure 3.
D, dimension display technologies:
The integrated AutoCAD core dynamic link library exploitation under the VC++ environment of this part forms, mainly be that monitoring result, analysis result and engineering three-dimensional stereo model image are shown, by the monitoring information database, geometric model information database and analysis result stored data base are formed, three databases are used for storing characteristics such as concrete physical dimension, position and the measuring point placement information at the various shape informations of microseism monitoring of equipment and feature, each position of Geotechnical Engineering respectively and utilize the result of above-mentioned intellectual technology Treatment Analysis, the three connects each other, information mutual communication.
The present invention compared with prior art has the following advantages and good effect:
1) sensor intelligent optimizing method for disposing, considered of the constraint of actual mining technology to sensing station, microseismic signals energy attenuation and tomography, goaf are to the influence of microseism wave velocity, under the certain situation of transducer sensitivity and monitoring range, take all factors into consideration transducer arrangements cost and microquake sources bearing accuracy, with the population intelligent optimization method for finding the solution means, may search for optimization in the position global space to sensor, reached and obtain manyly, more comprehensively and the purpose of more complete effective microseism information, Fig. 4 is the example that domestic certain metal mine is used; 2) method of utilizing genetic algorithm and neural network to combine, adopt many indexs to carry out filter analyses to the microseismic signals under many noise jamming, solved the difficult problem that single index filtering technique is difficult to handle the microseismic signals that many noises are interweaved, Fig. 5 and Fig. 6 are filter effect figure.3) microquake sources intelligent locating method, on the basis of handling (disappear dry, filtering etc.) and transducer arrangements optimization to ripple, consider of the influence of the heterogeneity of tomography, goaf and rock material to velocity of wave, based on time difference positioning principle, adopt population colony intelligence method to microseism wave velocity model and focal shock parameter (focus three-dimensional coordinate and origin time) joint inversion, the shortcoming of easily having dispersed when both having solved the classic method location, improved bearing accuracy again, Fig. 7 is the design sketch that this method is used in certain Tunnel Engineering; 4) optimize 3-D display as a result, integrated AutoCAD core dynamic link library exploitation forms under the VC++ environment, can be visual in image, aspect 3-D display transducer arrangements and microquake sources intelligence positioning result quickly.
Description of drawings
Fig. 1 is a kind of transducer arrangements optimisation technique process flow diagram;
Fig. 2 filters the synoptic diagram of making an uproar for the many indexs of artificial neural network;
Fig. 3 is a microquake sources intelligence layering positioning flow synoptic diagram;
Fig. 4 is certain stage casing, mine transducer arrangements floor map;
Fig. 5 is acoustie emission event rate before the filtering;
Fig. 6 is acoustie emission event rate after the filtering;
The design sketch of Fig. 7 for using in the microquake sources intelligence method for positioning analyzing engineering;
Fig. 8 is a numerical analysis integral grid illustraton of model;
Fig. 9 is the mining area cut-away view;
Figure 10 is a 3-D display software operation function menu.
Embodiment
A kind of microquake sources intelligent locating method (being applied as example explanation the specific embodiment of the present invention in certain mine with the present invention) describes in further detail the present invention below in conjunction with accompanying drawing, and a kind of microquake sources intelligent locating method the steps include:
1. based on numerical analysis and PSO optimization method microquake sources monitoring sensor layout optimization technology
Based on the zone that Primary Evaluation rock pressure and geologic hazard might occur, be target so that sensor monitors is wider, monitoring accuracy is higher, in conjunction with the microquake sources location algorithm, consider the influence of mining technology, utilize the PSO technology, sensing station is optimized, detailed process is as follows:
(1) at first, determine mechanics analysis model and parameter thereof according to shop experiment and field monitoring, according to geologic condition (big tomography, structural plane), consider to have foundation block chambering method mining technology characteristics, set up numerical analysis model, as Fig. 8 and as shown in Figure 9, with the energy release rate is evaluation index, numerical analysis and evaluation are carried out in the mining activity, hold the evolution rule of stress field in the recovery process on the whole, the zone that Primary Evaluation rock pressure and geologic hazard might occur is as the key area of microseismic monitoring sensor layout optimization.
(2) determine the monitoring range of 600m * 600m * 550m according to the actual recovery activity in mine,, determine 800 microquake sources in conjunction with the numerical analysis result.
(3), determine that microseism velocity of wave propagation model is V (v according to the distribution of heterogeneity ore deposit, mine rock wave velocity test, goaf wave velocity test and heterogeneity rock stratum 1, v 2, v 3, v 4) and velocity of wave span 4500m.s -1--6500m.s -1
(4) the scope 600m of transducer arrangements * 600m * 550m utilizes mixed congruence method to produce 20 groups of sensing stations at random, judges whether tomography and goaf are arranged between sensor and the microquake sources, selects the proper speed model, according to formula (a)
Figure GSA00000005848600081
Calculate 24000 sensor monitors then;
(b) t by formula when (5) calculating V=(1+3c (1) xa 0) t cCarry out random perturbation, obtain 24000 Virtual Monitoring then.
According to then virtual, microquake sources is positioned analysis, considering site operation easily under the prerequisite, make when monitoring and the accumulation residual sum of squares (RSS) when calculating less than 10 -4, finish to optimize, optimize as a result design sketch as shown in Figure 4.
2. many indexs of microseismic signals intelligent filter technology
1), determines that Construction traffic, electrical equipment noise, roofbolter, air compressor machine, fine bar such as knock at the feature of noise, set up noise database by site test; 2) be the neural network input with 12 characteristic quantities such as rise time, tale, peak counting, energy, amplitude duration, 5 kinds of noise types are that neural network output makes up 50 of neural network learnings and test sample book 10; 3) utilize structure, weights and the threshold value of genetic algorithm optimization neural network, determine that the neural network optimum structure is 12-34-8-1; 4) set up mapping relations between different noise signal features and the noise types by training; 5) Input Monitor Connector to microseismic signals it is carried out de-noising, filtering, obtain effective microseismic signals.Filter effect as shown in Figure 5 and Figure 6.
3. microquake sources intelligence layering method for positioning analyzing
(1) at first utilize the intelligent filter technology that sensor array is monitored the processing such as dry 1, filtering 2 that disappear of many noise signals, farthest the accurate recording microseismic signals 3, choose effective microseismic signals 5 and carry out the microquake sources location.
(2) then,, obtain useful signal more accurately to then revising again and picking up 6;
(3) initialization PS0 parameter 7: study factor c 1=c 2=2, population size N Pop=100, w 0=1, termination condition ε 0=1.0 * 10 -10Least square method termination condition ε 0=1.0 * 10 -10The microquake sources scope is that 600m * 600m * 550m and velocity of wave scope are 4500m.s -1-6500m.s -1
(4) utilize mixing congruence random algorithm initialization population position (source location and velocity of wave 8) and particle flying speed (iteration step length).
(5) according to formula (c) Calculate the adaptive value 9 of particle, and judge whether to satisfy flight number of times and the bearing accuracy of setting in advance 10, satisfy, carry out (6) step of this section; Otherwise, carry out (5) step of this section;
(6) according to formula (d) V Id=wV Id+ c 1r 1(P Id-X Id)+c 2r 2(P Gd-X Id) and (e) X Id=X Id+ V IdMore the position of new particle (microquake sources) and flying speed (iteration step length) 11 returned for (3) step;
(7) with the microquake sources coordinate and the rate pattern substitution formula (f) that recognize
Figure GSA00000005848600092
Obtain microquake sources origin time t12;
(8) judge that microseismic signals propagates into sensor and whether passes big tomography, dead zone 14, reject this signal, return (2) step of this section and reorientate if pass; Otherwise it is that correct and microquake sources location are accurately 15 that microseismic signals is chosen, and finishes microquake sources location 16.
Locating effect as shown in Figure 7.
4. dimension display technologies
According to independently developed real software (master menu as shown in figure 11), can directly show the characteristic such as concrete physical dimension, position and measuring point placement information at the various shape informations of micro seismic monitoring and feature, each each position of Geotechnical Engineering, mine efficiently comprehensively and utilize the result of above-mentioned intellectual technology Treatment Analysis.

Claims (1)

1. a microquake sources intelligent locating method the steps include:
A, arrange based on numerical analysis and PSO method microquake sources monitoring sensor:
(1), according to the actual mining technology in mine, geologic condition is set up numerical simulator, and the movable advanced line number value of mining is analyzed, and holds the zone that the mining geology disaster occurs in the recovery process on the whole, as the microseismic monitoring sensor layout area;
(2) in the mining scope of activities,, determine N microquake sources in conjunction with the numerical analysis result;
(3) according to the distribution of heterogeneity ore deposit, mine rock wave velocity test, goaf wave velocity test and heterogeneity rock stratum, determine microseism velocity of wave propagation model and span;
(4) then,, determine the scope of transducer arrangements, and in scope, utilize mixed congruence method to produce X group sensing station, every group of M sensing station at random according to the actual mining technology in mine;
(5) to every group of sensor, judge that the focus microseismic signals propagates into tomography and the goaf that sensor experiences, and calculate the distance that microseismic signals is propagated therein, select rate pattern automatically, calculate according to formula meter (a) and obtain N * M sensor monitors then:
t ij = T i + L ij V - - - ( a )
Wherein, t IjBe the time that j wave detector receives the signal that i microquake sources send, T iBe i the microquake sources origin time of earthquake, V is the equivalent velocity of wave that ripple is propagated in medium, L IjBe the distance of i microquake sources, be calculated as follows to j wave detector:
L ij = ( x j - x i ) 2 + ( y j - y j ) 2 + ( z j - z i ) 2
In the formula, (x i, y i, z i) an i microquake sources position coordinates, (x j, y j, z j) be k wave detector position coordinates;
(6) when calculating by formula (b) carry out random perturbation, obtain N * M Virtual Monitoring then:
t V=(1+3c(-1) xa 0)t c (b)
Wherein, t cFor Virtual Monitoring then, t cWhen monitoring for calculating, a 0Be the number between the 0-0.05, c is the random number between 0-1, and x=(int) 3c rounds 3 times random numbers;
(7) under the prerequisite of transducer arrangements, the accumulation residual sum of squares (RSS) when making when monitoring and calculating is minimum to be to utilize population colony intelligence method by objective function, according to M sensor monitors signal, N microquake sources is positioned;
(8) bearing accuracy meets the demands, and sensor monitors is to microseismic signals, and the position is reasonable, finishes transducer arrangements;
(9) do not find sensing station, utilize the PSO operation, in the scope of transducer arrangements, produce the new sensing station of X group, return (4) step of this section, new one group sensing station is judged;
B, many indexs of microseismic signals intelligent filter technology, at first set up various types of noise databases by test, to the noise signal feature: rise time, tale, peak counting, energy, amplitude duration are carried out analytic induction, make up neural network learning and test sample book; Then, utilize structure, weights and the threshold value of genetic algorithm optimization neural network, utilize superpower non-linear mapping capability of artificial neural network and self-learning capability neural network training, set up the mapping relations between different noise signal features and the noise types; At last, Input Monitor Connector to microseismic signals it is carried out filtering, obtain effective microseismic signals;
C, microquake sources intelligence layering location:
(1) utilize intelligent filter that sensor array is monitored many noise signals dry, Filtering Processing that disappears, farthest the accurate recording microseismic signals is then chosen and is used for the microseism ripple type of focus location;
(2) initialization PSO parameter, microquake sources scope and velocity of wave scope;
(3) utilize mixing congruence random algorithm initialization population position and particle flying speed;
(4) according to time difference positioning principle, by the value of formula (c) calculating particle, judge flight number of times and the bearing accuracy set in advance, satisfy, carry out this section (6) and go on foot;
f = Σ k = 1 n [ T k M - T k C ] 2 - - - ( c )
In the formula, T k MAnd T k CBe respectively k sensor monitors and when calculating;
(5) microquake sources coordinate and the rate pattern substitution formula (f) that recognizes obtained microquake sources origin time t;
t = Σ k = 1 n ( w k - L k V ) n - - - ( f )
(6) judge that microseismic signals propagates into sensor and passes tomography, dead zone, pass and reject this signal, return (2) step of this section and reorientate;
D, dimension display technologies, this part integrated AutoCAD dynamic link library under the VC++ environment is developed to, be that monitoring result, analysis result and engineering three-dimensional stereo model image are shown, by the monitoring information database, geometric model information database and analysis result stored data base are formed, three databases are used for storing various shape informations and feature, concrete physical dimension, position and the measuring point placement information characteristic at each position of Geotechnical Engineering and the result who utilizes above-mentioned Intelligent treatment to analyze of microseism monitoring of equipment respectively, the three connects each other, information mutual communication.
CN2010101005273A 2010-01-22 2010-01-22 Intelligent positioning method of mine microquake sources Expired - Fee Related CN101770038B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101005273A CN101770038B (en) 2010-01-22 2010-01-22 Intelligent positioning method of mine microquake sources

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101005273A CN101770038B (en) 2010-01-22 2010-01-22 Intelligent positioning method of mine microquake sources

Publications (2)

Publication Number Publication Date
CN101770038A true CN101770038A (en) 2010-07-07
CN101770038B CN101770038B (en) 2012-08-22

Family

ID=42503016

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101005273A Expired - Fee Related CN101770038B (en) 2010-01-22 2010-01-22 Intelligent positioning method of mine microquake sources

Country Status (1)

Country Link
CN (1) CN101770038B (en)

Cited By (50)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102129063A (en) * 2010-12-23 2011-07-20 中南大学 Method for positioning micro seismic source or acoustic emission source
CN102262220A (en) * 2011-04-28 2011-11-30 中南大学 Positioning method based on non-linear fitting micro-seismic source or acoustic emission source
CN102778668A (en) * 2012-07-23 2012-11-14 中煤科工集团西安研究院 Method for quickly and precisely positioning passive mine focus
CN103410565A (en) * 2013-03-14 2013-11-27 天地科技股份有限公司 Monitoring system and early warning method for rock burst multi-parameter process
CN103576188A (en) * 2012-07-26 2014-02-12 中国石油化工股份有限公司 Seismic source location method eliminating influences of velocity errors
CN103605151A (en) * 2013-11-20 2014-02-26 中北大学 Distributed group wave shallow-layer slight shock positioning method based on phase measuring
CN103744055A (en) * 2013-12-04 2014-04-23 桂林电子科技大学 Mining area illegal mining prevention monitoring and positioning method and device thereof
CN103953392A (en) * 2014-05-07 2014-07-30 中国科学院武汉岩土力学研究所 Method for distinguishing position of rockburst risk on deep buried tunnel section
CN103984026A (en) * 2014-05-30 2014-08-13 江苏三恒科技股份有限公司 Micro-vibration wave shape automatic labeling system
CN104360391A (en) * 2014-12-02 2015-02-18 武汉科技大学 Micro seismic source positioning method based on beam bunching array waveform
CN104406681A (en) * 2014-11-21 2015-03-11 中国矿业大学 Testing method for determining microquake wave velocity in real time
CN104502964A (en) * 2014-12-19 2015-04-08 桂林电子科技大学 Method for obtaining microearthquake wave velocity based on space geometry relationship
CN104536035A (en) * 2015-01-16 2015-04-22 淮南矿业(集团)有限责任公司 Method for obtaining position of centrum of coal measure stratum
CN104834004A (en) * 2015-04-13 2015-08-12 中南大学 Mine slight shock and blasting signal identification method based on waveform slope before and after peak value
CN105093314A (en) * 2015-07-10 2015-11-25 中联煤层气有限责任公司 Method for measuring and determining micro-seismic focus
CN105607040A (en) * 2015-09-07 2016-05-25 中国神华能源股份有限公司 Mining area illegal mining prevention monitoring and positioning method and system
CN105785436A (en) * 2016-03-17 2016-07-20 北京矿冶研究总院 Mining micro-seismic monitoring method
CN106680867A (en) * 2016-11-17 2017-05-17 大连理工大学 Dynamic parameter method for accurate positioning of micro-seismic event
CN106990435A (en) * 2017-06-07 2017-07-28 中煤科工集团西安研究院有限公司 It is a kind of to weaken the microseism localization method and device for relying on first break pickup precision
CN107884822A (en) * 2017-11-13 2018-04-06 北京矿冶研究总院 Method for improving positioning precision of mining micro-seismic source
CN107918279A (en) * 2017-11-20 2018-04-17 上海交通大学 A kind of TBM vibration-reducing control methods based on particle swarm optimization algorithm PSO
CN108896397A (en) * 2018-07-17 2018-11-27 西南大学 Roof greening charge of surety evaluation method based on On Microseismic Monitoring Technique
CN108919358A (en) * 2018-04-23 2018-11-30 中国矿业大学 A kind of mine quake disaster differentiates and signal reconfiguring method
CN108931816A (en) * 2018-08-17 2018-12-04 山东省科学院激光研究所 A kind of seismic source location method and device
CN109033607A (en) * 2018-07-19 2018-12-18 山东科技大学 A kind of optimization method of microseism seismic source location parameter
CN109061723A (en) * 2018-05-18 2018-12-21 中国科学院武汉岩土力学研究所 A kind of the microquake sources high-precision locating method and system of tunnel rock burst preparation process
CN109557588A (en) * 2018-11-16 2019-04-02 徐州工程学院 A kind of underground coal mine two dimension mine shake Velocity Inversion dimension reduction method
CN109779635A (en) * 2019-02-02 2019-05-21 韩少鹏 A kind of tunnel Engineering safe excavation method
CN109958475A (en) * 2019-04-28 2019-07-02 贵州大学 A kind of microseismic sensors structure and recognition methods identifying subtle disturbing signal
CN110018062A (en) * 2019-05-07 2019-07-16 中国科学院武汉岩土力学研究所 Rock structural face failure by shear location positioning method in a kind of direct shear test
CN110333535A (en) * 2019-04-03 2019-10-15 中国科学院武汉岩土力学研究所 A kind of scene rockmass anisotropy velocity of wave field measurement method in situ
CN110414675A (en) * 2019-09-02 2019-11-05 中北大学 A kind of underground shallow layer seismic source location method based on deep learning
CN110457758A (en) * 2019-07-16 2019-11-15 江西理工大学 Prediction technique, device, system and the storage medium in Instability of Rock Body stage
CN110609321A (en) * 2019-09-24 2019-12-24 中国科学院武汉岩土力学研究所 Micro seismic source positioning method based on speed model database
CN111405469A (en) * 2020-03-24 2020-07-10 辽宁大学 Mine earthquake monitoring system based on mobile phone mobile sensing network and crowd-sourcing positioning method
CN112505699A (en) * 2020-11-26 2021-03-16 中国矿业大学 Method for inverting underground goaf position parameters by fusing InSAR and PSO
CN112904414A (en) * 2021-01-19 2021-06-04 中南大学 Earth sound event positioning method and instability disaster early warning method thereof, earth sound perception instrument, monitoring system and readable storage medium thereof
CN112946735A (en) * 2021-02-23 2021-06-11 石家庄铁道大学 Rockfall impact positioning method and device based on micro-seismic monitoring system
CN113050158A (en) * 2021-03-19 2021-06-29 中国科学院武汉岩土力学研究所 Analysis method, device and equipment for near-field microseismic signal waveform and storage medium
CN113740899A (en) * 2021-10-20 2021-12-03 辽宁工程技术大学 Coal mine stope seismic source monitoring and positioning system based on wireless transmission
CN113803067A (en) * 2021-08-13 2021-12-17 山东省煤田地质规划勘察研究院 Local rock burst prevention and control device for coal mine
CN113917562A (en) * 2021-09-27 2022-01-11 中国科学院武汉岩土力学研究所 Macro-micro structure representation and three-dimensional space construction method and device for deep-buried soft interlayer
CN114047546A (en) * 2021-11-18 2022-02-15 辽宁大学 Crowd-sourcing spiral mine earthquake positioning method based on three-dimensional spatial joint arrangement of sensors
CN114371503A (en) * 2021-12-10 2022-04-19 煤炭科学技术研究院有限公司 Seismic source positioning method and device, electronic equipment and storage medium
CN114545500A (en) * 2022-01-29 2022-05-27 煤炭科学研究总院有限公司 Method and device for determining wave velocity distribution information
CN114563826A (en) * 2022-01-25 2022-05-31 中国矿业大学 Microseism sparse table network positioning method based on deep learning fusion drive
US11628445B2 (en) 2019-04-09 2023-04-18 Jiangxi University Of Science And Technology Material crushing cavity structure and method for designing a multi-stage nested material crushing cavity structure
CN117991349A (en) * 2024-04-07 2024-05-07 吉林大学 Microseism positioning method based on improved ant lion optimization algorithm
US11989655B2 (en) 2020-10-26 2024-05-21 Jiangxi University Of Science And Technology Prediction method, device and system for rock mass instability stages
CN118311646A (en) * 2024-06-07 2024-07-09 山东省地质矿产勘查开发局第三地质大队(山东省第三地质矿产勘查院、山东省海洋地质勘查院) Data analysis method for geological environment monitoring

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110488350B (en) * 2019-09-20 2021-10-29 西南石油大学 Seismic inversion big data generation method based on convolutional neural network

Cited By (77)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102129063A (en) * 2010-12-23 2011-07-20 中南大学 Method for positioning micro seismic source or acoustic emission source
CN102129063B (en) * 2010-12-23 2012-10-10 中南大学 Method for positioning micro seismic source or acoustic emission source
CN102262220A (en) * 2011-04-28 2011-11-30 中南大学 Positioning method based on non-linear fitting micro-seismic source or acoustic emission source
CN102262220B (en) * 2011-04-28 2013-07-17 中南大学 Positioning method of micro-seismic source or acoustic emission source based on non-linear fitting
CN102778668A (en) * 2012-07-23 2012-11-14 中煤科工集团西安研究院 Method for quickly and precisely positioning passive mine focus
CN103576188A (en) * 2012-07-26 2014-02-12 中国石油化工股份有限公司 Seismic source location method eliminating influences of velocity errors
CN103576188B (en) * 2012-07-26 2016-09-14 中国石油化工股份有限公司 A kind of seismic source location method of release rate error impact
CN103410565A (en) * 2013-03-14 2013-11-27 天地科技股份有限公司 Monitoring system and early warning method for rock burst multi-parameter process
CN103410565B (en) * 2013-03-14 2016-08-10 天地科技股份有限公司 Bump many reference amounts process monitoring system and method for early warning
CN103605151A (en) * 2013-11-20 2014-02-26 中北大学 Distributed group wave shallow-layer slight shock positioning method based on phase measuring
CN103744055A (en) * 2013-12-04 2014-04-23 桂林电子科技大学 Mining area illegal mining prevention monitoring and positioning method and device thereof
CN103744055B (en) * 2013-12-04 2016-08-17 桂林电子科技大学 A kind of mining area anti-illegal mining monitoring and positioning method and equipment thereof
CN103953392B (en) * 2014-05-07 2015-12-02 中国科学院武汉岩土力学研究所 Rockburst risk position method of discrimination on deep tunnel section
CN103953392A (en) * 2014-05-07 2014-07-30 中国科学院武汉岩土力学研究所 Method for distinguishing position of rockburst risk on deep buried tunnel section
CN103984026B (en) * 2014-05-30 2017-08-22 江苏三恒科技股份有限公司 Microseism waveform automatic marking system
CN103984026A (en) * 2014-05-30 2014-08-13 江苏三恒科技股份有限公司 Micro-vibration wave shape automatic labeling system
CN104406681B (en) * 2014-11-21 2015-11-18 中国矿业大学 A kind of method of testing determining microseism velocity of wave in real time
CN104406681A (en) * 2014-11-21 2015-03-11 中国矿业大学 Testing method for determining microquake wave velocity in real time
CN104360391A (en) * 2014-12-02 2015-02-18 武汉科技大学 Micro seismic source positioning method based on beam bunching array waveform
CN104360391B (en) * 2014-12-02 2017-02-22 武汉科技大学 Micro seismic source positioning method based on beam bunching array waveform
CN104502964A (en) * 2014-12-19 2015-04-08 桂林电子科技大学 Method for obtaining microearthquake wave velocity based on space geometry relationship
CN104536035A (en) * 2015-01-16 2015-04-22 淮南矿业(集团)有限责任公司 Method for obtaining position of centrum of coal measure stratum
CN104536035B (en) * 2015-01-16 2017-03-22 淮南矿业(集团)有限责任公司 Method for obtaining position of centrum of coal measure stratum
CN104834004A (en) * 2015-04-13 2015-08-12 中南大学 Mine slight shock and blasting signal identification method based on waveform slope before and after peak value
CN104834004B (en) * 2015-04-13 2017-04-05 中南大学 Mine microquake based on pre-and post-peaking waveform slope and blast signal recognition methodss
CN105093314B (en) * 2015-07-10 2017-09-22 中联煤层气有限责任公司 A kind of method for determining microseism focus
CN105093314A (en) * 2015-07-10 2015-11-25 中联煤层气有限责任公司 Method for measuring and determining micro-seismic focus
CN105607040A (en) * 2015-09-07 2016-05-25 中国神华能源股份有限公司 Mining area illegal mining prevention monitoring and positioning method and system
CN105785436A (en) * 2016-03-17 2016-07-20 北京矿冶研究总院 Mining micro-seismic monitoring method
CN105785436B (en) * 2016-03-17 2018-08-14 北京矿冶研究总院 Mining micro-seismic monitoring method
CN106680867A (en) * 2016-11-17 2017-05-17 大连理工大学 Dynamic parameter method for accurate positioning of micro-seismic event
CN106990435A (en) * 2017-06-07 2017-07-28 中煤科工集团西安研究院有限公司 It is a kind of to weaken the microseism localization method and device for relying on first break pickup precision
CN107884822B (en) * 2017-11-13 2019-09-27 北京矿冶研究总院 Method for improving positioning precision of mining micro-seismic source
CN107884822A (en) * 2017-11-13 2018-04-06 北京矿冶研究总院 Method for improving positioning precision of mining micro-seismic source
CN107918279A (en) * 2017-11-20 2018-04-17 上海交通大学 A kind of TBM vibration-reducing control methods based on particle swarm optimization algorithm PSO
CN108919358A (en) * 2018-04-23 2018-11-30 中国矿业大学 A kind of mine quake disaster differentiates and signal reconfiguring method
CN109061723A (en) * 2018-05-18 2018-12-21 中国科学院武汉岩土力学研究所 A kind of the microquake sources high-precision locating method and system of tunnel rock burst preparation process
CN109061723B (en) * 2018-05-18 2020-07-10 中国科学院武汉岩土力学研究所 High-precision positioning method and system for micro seismic source in tunnel rock burst inoculation process
CN108896397B (en) * 2018-07-17 2021-04-27 西南大学 Roof greening safety load evaluation method based on microseismic monitoring technology
CN108896397A (en) * 2018-07-17 2018-11-27 西南大学 Roof greening charge of surety evaluation method based on On Microseismic Monitoring Technique
CN109033607A (en) * 2018-07-19 2018-12-18 山东科技大学 A kind of optimization method of microseism seismic source location parameter
CN108931816B (en) * 2018-08-17 2020-01-21 山东省科学院激光研究所 Seismic source positioning method and device
CN108931816A (en) * 2018-08-17 2018-12-04 山东省科学院激光研究所 A kind of seismic source location method and device
CN109557588B (en) * 2018-11-16 2020-08-28 徐州工程学院 Coal mine underground two-dimensional mine seismic wave velocity inversion dimension reduction method
CN109557588A (en) * 2018-11-16 2019-04-02 徐州工程学院 A kind of underground coal mine two dimension mine shake Velocity Inversion dimension reduction method
CN109779635A (en) * 2019-02-02 2019-05-21 韩少鹏 A kind of tunnel Engineering safe excavation method
CN110333535A (en) * 2019-04-03 2019-10-15 中国科学院武汉岩土力学研究所 A kind of scene rockmass anisotropy velocity of wave field measurement method in situ
US11628445B2 (en) 2019-04-09 2023-04-18 Jiangxi University Of Science And Technology Material crushing cavity structure and method for designing a multi-stage nested material crushing cavity structure
CN109958475A (en) * 2019-04-28 2019-07-02 贵州大学 A kind of microseismic sensors structure and recognition methods identifying subtle disturbing signal
CN110018062B (en) * 2019-05-07 2020-05-08 中国科学院武汉岩土力学研究所 Method for positioning shearing failure position of rock structural surface in direct shear test
CN110018062A (en) * 2019-05-07 2019-07-16 中国科学院武汉岩土力学研究所 Rock structural face failure by shear location positioning method in a kind of direct shear test
CN110457758A (en) * 2019-07-16 2019-11-15 江西理工大学 Prediction technique, device, system and the storage medium in Instability of Rock Body stage
CN110414675A (en) * 2019-09-02 2019-11-05 中北大学 A kind of underground shallow layer seismic source location method based on deep learning
CN110414675B (en) * 2019-09-02 2022-05-27 中北大学 Underground shallow seismic source positioning method based on deep learning
CN110609321B (en) * 2019-09-24 2020-11-13 中国科学院武汉岩土力学研究所 Micro seismic source positioning method based on speed model database
CN110609321A (en) * 2019-09-24 2019-12-24 中国科学院武汉岩土力学研究所 Micro seismic source positioning method based on speed model database
CN111405469A (en) * 2020-03-24 2020-07-10 辽宁大学 Mine earthquake monitoring system based on mobile phone mobile sensing network and crowd-sourcing positioning method
US11989655B2 (en) 2020-10-26 2024-05-21 Jiangxi University Of Science And Technology Prediction method, device and system for rock mass instability stages
CN112505699A (en) * 2020-11-26 2021-03-16 中国矿业大学 Method for inverting underground goaf position parameters by fusing InSAR and PSO
CN112904414A (en) * 2021-01-19 2021-06-04 中南大学 Earth sound event positioning method and instability disaster early warning method thereof, earth sound perception instrument, monitoring system and readable storage medium thereof
WO2022156582A1 (en) * 2021-01-19 2022-07-28 中南大学 Earthquake sound event positioning method and instability disaster warning method based on same, earthquake sound sensor, monitoring system, and readable storage medium
CN112946735A (en) * 2021-02-23 2021-06-11 石家庄铁道大学 Rockfall impact positioning method and device based on micro-seismic monitoring system
CN113050158A (en) * 2021-03-19 2021-06-29 中国科学院武汉岩土力学研究所 Analysis method, device and equipment for near-field microseismic signal waveform and storage medium
CN113803067A (en) * 2021-08-13 2021-12-17 山东省煤田地质规划勘察研究院 Local rock burst prevention and control device for coal mine
CN113803067B (en) * 2021-08-13 2024-01-23 山东省煤田地质规划勘察研究院 Colliery local rock burst prevention and cure device
CN113917562B (en) * 2021-09-27 2023-02-28 中国科学院武汉岩土力学研究所 Macro-microscopic structure representation and three-dimensional space construction method and device for deep-buried soft interlayer
CN113917562A (en) * 2021-09-27 2022-01-11 中国科学院武汉岩土力学研究所 Macro-micro structure representation and three-dimensional space construction method and device for deep-buried soft interlayer
CN113740899A (en) * 2021-10-20 2021-12-03 辽宁工程技术大学 Coal mine stope seismic source monitoring and positioning system based on wireless transmission
CN114047546A (en) * 2021-11-18 2022-02-15 辽宁大学 Crowd-sourcing spiral mine earthquake positioning method based on three-dimensional spatial joint arrangement of sensors
CN114371503A (en) * 2021-12-10 2022-04-19 煤炭科学技术研究院有限公司 Seismic source positioning method and device, electronic equipment and storage medium
CN114371503B (en) * 2021-12-10 2023-08-29 煤炭科学技术研究院有限公司 Method and device for positioning seismic source, electronic equipment and storage medium
CN114563826A (en) * 2022-01-25 2022-05-31 中国矿业大学 Microseism sparse table network positioning method based on deep learning fusion drive
CN114563826B (en) * 2022-01-25 2023-03-03 中国矿业大学 Microseismic sparse table network positioning method based on deep learning fusion drive
CN114545500A (en) * 2022-01-29 2022-05-27 煤炭科学研究总院有限公司 Method and device for determining wave velocity distribution information
CN117991349A (en) * 2024-04-07 2024-05-07 吉林大学 Microseism positioning method based on improved ant lion optimization algorithm
CN118311646A (en) * 2024-06-07 2024-07-09 山东省地质矿产勘查开发局第三地质大队(山东省第三地质矿产勘查院、山东省海洋地质勘查院) Data analysis method for geological environment monitoring
CN118311646B (en) * 2024-06-07 2024-08-20 山东省地质矿产勘查开发局第三地质大队(山东省第三地质矿产勘查院、山东省海洋地质勘查院) Data analysis method for geological environment monitoring

Also Published As

Publication number Publication date
CN101770038B (en) 2012-08-22

Similar Documents

Publication Publication Date Title
CN101770038B (en) Intelligent positioning method of mine microquake sources
CN103306722B (en) Micro-seismic multi-dimensional information comprehensive region detection and evaluation method for impact danger region
CN105785471A (en) Impact danger evaluation method of mine pre-exploiting coal seam
CN111090709A (en) Big data geological analysis method for sandstone-type uranium ore mineralization prediction
CN116591777B (en) Multi-field multi-source information fusion rock burst intelligent monitoring and early warning device and method
CN111339486A (en) Deep foundation pit blasting vibration velocity risk level big data evaluation method
CN111062544A (en) Prediction method for uranium mineralization distant scenic region
CN103606019B (en) Mine goaf overlying stratum sedimentation dynamic prediction method based on time-space relationship
Wang et al. An auto-detection network to provide an automated real-time early warning of rock engineering hazards using microseismic monitoring
CN115185015B (en) Deep lithium beryllium ore investigation method
CN107798189A (en) Accurate dynamic outburst prevention method based on transparent space geophysical
CN105487117A (en) Three-dimensional earthquake observation system optimization method and apparatus
CN107942383A (en) Roof sandstone watery grade prediction technique
CN108197421B (en) Quantitative evaluation method for beneficial zone of joint development of dense gas and coal bed gas
CN114943149A (en) Method for calculating volume of rock mass damaged by rock burst in tunnel
Liu et al. Study of roof water inrush forecasting based on EM-FAHP two-factor model
Zhang Exploration on coal mining-induced rockburst prediction using Internet of things and deep neural network
Wu et al. Combination of seismic attributes using clustering and neural networks to identify environments with sandstone-type uranium mineralization
Liang et al. The influence factors of the stability of tailings dam based on multi-source information fusion method
CN110717618A (en) Submarine hydrothermal sulfide resource evaluation and prediction method based on multi-index comprehensive elements
Jianping et al. A 3-D Prediction Method for Blind Orebody Based on 3-D Visualization Model and Its Application
Cao et al. Construction and Application of “Active Prediction-Passive Warning” Joint Impact Ground Pressure Resilience Prevention System: Take the Kuan Gou Coal Mine as an Example
CN114415237B (en) Sandstone-type uranium ore control fracture identification method and system based on three-dimensional seismic data
US12123995B1 (en) Intelligent monitoring and early warning device and method for rock burst based on multi-field and multi-source information fusion
Fan et al. Application of Seismic Channel Wave Technology on Small Structure Exploration in Coal Mine

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120822

Termination date: 20190122