NL2029589B1 - Identification Method of Regional Impact Risk Grade and Region Based on Inversion of Source Parameters - Google Patents

Identification Method of Regional Impact Risk Grade and Region Based on Inversion of Source Parameters Download PDF

Info

Publication number
NL2029589B1
NL2029589B1 NL2029589A NL2029589A NL2029589B1 NL 2029589 B1 NL2029589 B1 NL 2029589B1 NL 2029589 A NL2029589 A NL 2029589A NL 2029589 A NL2029589 A NL 2029589A NL 2029589 B1 NL2029589 B1 NL 2029589B1
Authority
NL
Netherlands
Prior art keywords
source
target area
impact
area
source parameters
Prior art date
Application number
NL2029589A
Other languages
Dutch (nl)
Inventor
Wang Changbin
Ma Zhifeng
Liu Yaoqi
Guo Wenhao
Shangguan Kefeng
Wang Songwei
Bai Xianqi
Xue Chengchun
Cao Anye
Peng Yujie
Wen Yingyuan
Original Assignee
Huating Coal Ind Group Co Ltd
Univ China Mining
Xuzhou Wushuo Information Tech Co Ltd
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 Huating Coal Ind Group Co Ltd, Univ China Mining, Xuzhou Wushuo Information Tech Co Ltd filed Critical Huating Coal Ind Group Co Ltd
Priority to NL2029589A priority Critical patent/NL2029589B1/en
Application granted granted Critical
Publication of NL2029589B1 publication Critical patent/NL2029589B1/en

Links

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
    • G01V1/288—Event detection in seismic signals, e.g. microseismics
    • 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
    • G01V1/282—Application of seismic models, synthetic seismograms
    • 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
    • G01V1/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/65—Source localisation, e.g. faults, hypocenters or reservoirs

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a method for identifying regional impact risk grade and region based on inversion of source parameters, which comprises the following steps: establishing a spatial coordinate system of a target area and screening mine 5 earthquake data; establishing the index system of hazard grade assessment; calculating the evaluation index value D; identifying the level of impact risk according to the index value D; determining the scope of impact risk and formulate preventive measures. According to the method; the earthquake source parameters are comprehensively utilized to predict the impact hazard level and the local area; and the 10 obtained results are convenient for field implementation.

Description

Identification Method of Regional Impact Risk Grade and Region Based on
Inversion of Source Parameters
TECHNICAL FIELD
The invention relates to a method for judging the danger level and the local area of rock burst, in particular to a method for judging the regional danger level and the local area of rock burst based on inversion of seismic source parameters for underground microseismic monitoring in coal mines.
BACKGROUND
China's coal mining depth is increasing rapidly at a rate of 10-25m per year, which leads to a sharp increase in the stress level of the original rock around the mining space. In addition, the complexity of coal seam occurrence makes the coal rock dynamic disaster represented by rock burst one of the key factors hindering efficient coal mining. The prevention and control of rock burst, especially the accurate detection of impact dangerous areas, is one of the worldwide problems. Geophysical methods, represented by microseismic monitoring, can monitor the impact dangerous areas in real time and in a large scale, and the information can be transmitted in real time, thus realizing diversified and centralized processing with other early warning information, which has broad prospects in prevention and control of rock burst risk and prediction and early warning.
At present, scholars at home and abroad have put forward many microseismic early warning indicators, such as energy, frequency, time-space diffusion, apparent stress, apparent volume, spatial concentration, etc. However, these indicators realize early warning through abnormal changes of index values, such as "abnormal changes of index values" and "abnormal high or low values", which are only judged qualitatively, resulting in low prediction efficiency. In addition, active and passive CT inversion based on microseisms, which has developed rapidly in recent years, has achieved good results in the detection of impact dangerous areas, but its long inversion period, manual shooting, wiring (which leads to increased cost and labor intensity) and mutual interference with production have become one of the main problems restricting its development. At the same time, in some typical rock-burst mining areas in Poland and China, it has been observed that low-energy rock-burst (energy less than 10*]) often induces rock-burst events, while high-energy rock-burst (energy greater than 10*J) does not induce rock-burst events, which is one of the main problems that lead to an increase in the false alarm rate of monitoring and early warning. Therefore, it is of great significance to fully mine the source information and put forward a method to identify the impact risk level and region by integrating various sources.
SUMMARY
In order to overcome the shortcomings of the prior art, the invention provides a method for judging the regional impact risk grade and the local area based on the inversion of the seismic source parameters, which is a regional impact risk detection method which comprehensively utilizes the seismic source parameter information, quantitatively evaluates the impact risk degree and judges the local area of the impact risk.
To achieve the above purpose, the present invention provides a scheme including the following steps:
Sl, establishing a target area spatial coordinate system and screening mine earthquake data;
S2, establishing an index system for evaluating the impact risk level of the target area, which comprises a strength index Q and a dispersion index H;
S3, calculating the evaluation index value of the impact hazard level of the target pC area: H
S4, judging the impact risk level of the target area according to the evaluation index value D of the impact risk level of the target area, wherein the greater the value
D is, the higher the risk level is;
S5, judging the impact hazard range according to the nearest distance between the focal center and the roadway and the relationship between the focal center and the focal radius, and formulating preventive measures according to the hazard grade and the range of the target area.
Preferably, the steps of establishing the spatial coordinate system of the target area and screening the mine earthquake data in step SI are:
S1.1, selecting a source with energy greater than 100J as a target area, establishing a source bank and analyzing;
S1.2, according to the spatial position of the target area, selecting a reference point 0 (x. Vo Zo) as an origin to establish a spatial rectangular coordinate system Z;
S1.3. according to the relative position with reference point O, all the focal coordinates in the seismic source bank are converted from world coordinates to relative coordinates, and the coordinates of any focal point are expressed as (572).
S1.4 according to the spatial relative position of the roadway in the target area and the reference point o, the geometric expression of the roadway in the coordinate system Z is listed: ”? n P| where “7/0 is the coordinate component of the coordinate (x Yo 2) of any point at the beginning and end of the roadway in the coordinate system Z, and m,n; p’ is the spatial direction vector of the roadway.
Preferably, the method for calculating the intensity index Q in step S2 is: ‚ | Vise
Q,=2. wh, W/'= Vv Vico i=1 | e J > Vi i in which, Q is the Q value of the intensity index of the j-th source in the source bank, "i is the prediction weight of the i-th source parameter in the source bank,
VP is the ith source parameter value of the j-th source in the source bank, where 1, 2 and 3 represent seismic moment, radiation energy and apparent stress respectively,
€ is the average value of the i-th source parameter of the j-th source in the source bank, and n is the total number of source parameter types used.
Preferably, the method for calculating the dispersion index H in step S2 is as follows: n
U= >" wyuy U>l
H _ ea ‘ U=1U<1 in which "7 is the prediction weight of correlation relation of source parameters, “i is the dispersion degree between pairwise corresponding source parameters,
Uy = 7 ds © which satisfies the relation Uus where x and y are a set of corresponding source parameter types, g(x) is a set of fitted correlation relation of source parameters, Um is a set of dispersion limit between source parameters, and n is the logarithm of source parameters used.
Preferably, the source parameters include seismic moment, radiation energy, apparent stress and source radius, which are calculated by Brune model.
Preferably, when solving the source parameters, the particle swarm optimization algorithm is used to solve the zero frequency limit Ay and the corner frequency Jo, and finally the source parameters are determined.
Preferably, in step S5, the specific scheme for formulating preventive measures according to the statistical regional danger level and the local area is as follows: when the impact danger level of the corresponding area is A, the roadway impact danger area does not need to carry out pressure relief work, but only needs to strengthen the detection of the target area; when the impact hazard level of the corresponding area is B, the roadway impact hazard area needs to be relieved and the detection of the target area should be strengthened; when the impact hazard level of the corresponding area is C, the personnel near the dangerous area shall be evacuated immediately, and the pressure relief work shall 5 be carried out after a certain period of time; when the dangerous areas of roadway judged by different danger levels coincide, the higher one shall prevail.
The method disclosed by the invention has the beneficial effects that the coal mine microseismic data are deeply excavated, so that the impact danger area and the danger grade can be quickly and accurately predicted, and relevant measures are taken in advance according to the prediction result, so that the impact danger in the area is reduced, and the operation safety of personnel is guaranteed; according to microseismic data, the danger zone is calculated in real time, which meets the needs of dynamic early warning of mine impact danger and meets the requirements of intelligent development of coal mines. The earthquake source parameters are comprehensively used to predict the impact hazard level and the local area, and the obtained results are convenient for on-site implementation, and the method established by the invention also has the characteristics of clear physical meaning and suitability for programming to realize intelligence.
BRIEF DESCRIPTION OF THE FIGURES
In order to explain the embodiments of the present invention or the technical scheme 1n the prior art more clearly, the following will briefly introduce the drawings used in the embodiments. Obviously, the drawings in the following description are only some embodiments of the present invention, and for ordinary technicians in the field, other drawings can be obtained according to these drawings without paying creative labor.
Figure 1 is a waveform record of each microseismic station corresponding to a source in the source library of the present invention;
Figure 2 is a plan view of the mining project in the target area of the present invention; in which ©) is working face, (2) is transport chute, (3) is return air chute;
Figure 3 is a spatial coordinate system of the target area of the present invention;
Figure 4 is a graph showing the statistical relationship between source parameters of the present invention;
Figure 5 is a schematic diagram of the layout of microseismic stations in the target area of the present invention; in which ©) is risk grade A, ©) is risk grade B, ©D) isrisk grade C, © is headentry, © is working face, {9 is tailentry, 1 is radius of hypocenter.
DESCRIPTION OF THE INVENTION
The following will clearly and completely describe the technical scheme in the embodiments of the present invention with reference to the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in the field without creative labor belong to the scope of protection of the present invention.
In order to make the above objects, features and advantages of the present invention more obvious and easy to understand, the present invention will be further explained in detail with reference to the drawings and specific embodiments.
Several microseismic events monitored during mining in 250105-1 working face of Yanbei Coal Mine, Gansu Province are selected for calculation. The waveform records of microseismic stations are shown in Figure 1. The plan of mining engineering in the target area is shown in Figure 2. (1) Establish the spatial model and spatial coordinate system of the target area as shown in Figure 3, and screen the mine earthquake data as follows: a, select the sources with energy greater than 100J in the target area to establish an analysis source bank; b, according to the spatial position of the target area, select the reference point 0 (x. Vo Zo) as the origin to establish a spatial rectangular coordinate system Z;
c‚ according to the relative position with reference point o, all the focal coordinates in the seismic source bank are converted from world coordinates to relative coordinates, and the coordinates of any seismic source can be expressed as (x,,V,Z;). d, according to the spatial relative position of the roadway in the target area and the reference point o, the geometric expression of the roadway in the coordinate system Z can be listed: n P where Xr Vr? is the coordinate component of the coordinate (x, ’ 2) of any point at the beginning and end of the roadway in the coordinate system Z, and m,n, D is the spatial direction vector of the roadway. (2) Establish an index system for evaluating the impact risk level of the target area, including strength index Q and dispersion index H; the calculation method of strength index Q is: ‚ I Vive — WJ j_ 072 ww, W, ie Vi cel in which, Q is the intensity index Q value of the jth source in the source library,
Wi is the prediction weight of the ith source parameter in the source library, and vi is the ith source parameter value of the jth source in the source library, where 1, 2 and 3 represent the seismic moment, radiation energy and apparent stress respectively, ei is the average value of the i-th source parameter of the j-th source in the source library, and n is the total number of source parameter types used.
The calculation method of dispersion index H is:
A
U= > wy UU >1 i=1,/=1
H={ ©
U=1,U<1 . . Wy . ee . . . in which, is the prediction weight of correlation relation of source parameters,
U; . oo ¥ is the dispersion degree between pairwise corresponding source parameters, yg, (x) : WT which satisfies the relation “mx where x and y are a group of gi (x) : corresponding source parameter types, is a group of fitted correlation relation u Ce of source parameters, Ymax is a group of dispersion limit between source parameters, and n is the logarithm of source parameters used. When there is historical training . p. w,, . : 7 data, the weights ”: and "7 are determined by using the prediction efficiency scoring method, and when there is no historical training data, the parameter values ; W,
Wi and / are equal.
In the process of calculating the dispersion index H, when calculating the . . U, - . . . u ‚ . . dispersion and dispersion limit “wax between source parameters, it is necessary to obtain the statistical relationship between a certain pair of source parameters as shown in Figure 4. According to this statistical relationship, the dispersion limit max . . ou can be obtained, and the dispersion # can be further calculated.
The seismic source parameters include seismic moment, radiation energy and apparent stress, which are calculated by Brune model:
Ampv FA _ 0 ee M, = 0 seismic moment: < in which, © is rock density, V is propagation speed of P wave or S wave, ” is . . A, ne . distance between source and station, ° is zero frequency limit of source, and Q is radiation factor of P wave or S wave (P wave is 0.52, S wave is 0.63),
U ' £ ! Us : 6 ’
E, =4mpv r* We) | i 0) di + 4mpv r? Vl ( dt radiation energy: Ph 5 in which is rock density, ° and VS are propagation velocity of P wave or S a U Ui wave respectively, r is distance between source and station, and S$ are source ee U, U; ij: : radiation terms respectively, | / n) and | ) are average values of Up and Us on source sphere respectively, and 4 and ft are starting point and ending point of P wave and S wave in time domain waveform respectively.
E. 0,=T ~
M, apparent stress: in which 7 is the shear modulus of rock mass at the earthquake source, E, is a M, . _. radiation energy of vibration wave and 1s seismic moment.
Use Brune model to calculate the focal radius: 234, Of I
F. “nr in which “¢ indicates the corner frequency of the seismic source
Af, and "Fis the P wave velocity.
Further, the source parameter value is finally determined according to the following formula: ¥=¢ ‚in which Pis the average value of some source mechanical parameters of the source; "is the focal mechanics parameter of each station and # is the total number of stations participating in the calculation.
Particle Swarm Optimization (PSO) is used to solve the zero frequency limit A, and corner frequency h. and the residual error between the theoretical value and the actual value in the iterative process is calculated by the following formula (when the residual error meets the conditions, the iteration is stopped):
EE
GETS
MES, in which “ indicates the residual error between the i-th theoretical source spectrum and the actual source spectrum; k is the theoretical source spectrum when i-th Jo and A, are adopted; Si is the actual source spectrum. (3) Calculate the evaluation index value of the impact risk level of the target area . ) according to the following formula: D= 2 : (4) Judging the impact risk level of the target area according to the evaluation index value D of the impact risk level of the target area, and the greater the value D, the higher the risk level. (5) According to the distance between the fracture radius of the source and the roadway, judge the boundary of impact danger, count the regional danger level and boundary, and formulate preventive measures.
The specific scheme is as follows: when D <05, the danger level is A, and the corresponding area 1s weak impact danger level. At this time, there is no need to carry out pressure relief work in the roadway impact danger area, only need to strengthen the monitoring of the target area; when 05D <0.75, the hazard level is B, and the corresponding area is of medium impact hazard level. At this time, pressure relief work should be carried out in the hazardous area, and monitoring should be strengthened.
When 220.75 the hazard level is C, and the corresponding area is of strong impact hazard level. At this time, people near the hazardous area should be evacuated immediately, and pressure relief work should be carried out after a certain period of time as shown in figure 5.
When the dangerous areas of roadway judged by different danger levels coincide, the higher one shall prevail.
A total of 27 seismic sources in the screened target area are shown in Table 1:
Table 1
Zero . a. .
Corner Seismic Radiation | Apparent | Radius of
Hypocenter frequency em EE frre mn A 6 /2160 |2.85E-06 |3.67E+12 | 1.ISE+0S |2.49E+02 [38.48 8 21.76 | 5.05E06 |7.21E+12 | 3.93E+04 | 127E+02 [3820 9 [2862 [109B05 | L65E+13 | 143E+05 | 5.55E+01 [29.04
The above embodiments only describe the preferred mode of the invention, but do not limit the scope of the invention. On the premise of not departing from the design spirit of the invention, various modifications and improvements made by ordinary technicians in the field to the technical scheme of the invention shall fall within the protection scope determined by the claims of the invention.
PREFERRED EMBODIMENTS
1. A method for judging the regional impact risk grade and the local area based on the inversion of the seismic source parameters is characterized by comprising the following steps:
S1, establishing a target area spatial coordinate system and screening mine earthquake data;
S2, establishing an index system for evaluating the impact risk level of the target area, which comprises a strength index Q and a dispersion index H;
S3, calculating the evaluation index value of the impact hazard level of the target pC area: H
S4, judging the impact risk level of the target area according to the evaluation index value D of the impact risk level of the target area, wherein the greater the value
D is, the higher the risk level is;
SS, judging the impact hazard range according to the nearest distance between the focal center and the roadway and the relationship between the focal center and the focal radius, and formulating preventive measures according to the hazard grade and the range of the target area. 2. The method for judging the regional impact risk grade and the local area based on the inversion of the seismic source parameters according to embodiment 1 is characterized in that the steps of establishing the spatial coordinate system of the target area and screening the mine earthquake data in S1 are as follows:
S1.1, selecting a source with energy greater than 100J as a target area, establishing a source bank and analyzing;
S1.2, according to the spatial position of the target area, selecting a reference point ACY as an origin to establish a spatial rectangular coordinate system Z,
S1.3. according to the relative position with reference point O, all the focal coordinates in the seismic source bank are converted from world coordinates to relative coordinates, and the coordinates of any focal point are expressed as (5.32),
S1.4 according to the spatial relative position of the roadway in the target area and the reference point o, the geometric expression of the roadway in the coordinate
X-x, JJ, ZZ, ic 1 om np Xp V2, ; : system Z is listed: , where ‘7’? + is the coordinate (x ‚Vz ) oo component of the coordinate * 7’ // of any point at the beginning and end of the roadway in the coordinate system Z, and m,n, p is the spatial direction vector of the roadway. 3. The method for judging the regional impact risk grade and the local area based on the inversion of the seismic source parameters according to embodiment 1 is characterized in that the method for calculating the intensity index Q in step S2 is as follow: 1 Vire! fn | 3
Q;7), WW W/=i Vi Co
Coe — JV <¢ i=l 7 > Vi i € in which, Q is the Q value of the intensity index of the j-th source in the source bank, ": is the prediction weight of the i-th source parameter in the source bank,
Vi! is the ith source parameter value of the j-th source in the source bank, where 1, 2 and 3 represent seismic moment, radiation energy and apparent stress respectively, ¢/ is the average value of the i-th source parameter of the j-th source in the source bank, and n 1s the total number of source parameter types used. 4. The method for judging the regional impact risk grade and the local area based on the inversion of the seismic source parameters according to embodiment 1 is characterized in that the method for calculating the dispersion index H in step S2 is as follow:
U= > wy; U>1 i=l j=1
H — i#f + among them, Wi is the prediction weight of correlation relation of source parameters, “i is the dispersion degree between pairwise corresponding source parameters, which satisfies the relation Umax where x and y are a set of corresponding source parameter types, g(x) is a set of fitted correlation relation of source parameters, Una js a set of dispersion limit between source parameters, and n is the logarithm of source parameters used. 5. The method for judging the regional impact risk grade and the local area based on the inversion of the seismic source parameters according to embodiment 3 is characterized in that the source parameters include seismic moment, radiation energy, apparent stress and source radius, which are calculated by Brune model. 6. The method for judging the regional impact risk grade and the local area based on the inversion of the seismic source parameters according to embodiment 3 is characterized in that when solving the source parameters, the particle swarm optimization algorithm is used to solve the zero frequency limit A and the corner frequency h. and finally the source parameters are determined. 7. The method for judging the regional impact risk grade and the local area based on the inversion of the seismic source parameters according to embodiment 1 is characterized in that in step S5, the specific scheme for formulating preventive measures according to the statistical regional danger level and the local area is as follows: when the impact danger level of the corresponding area is A, the roadway impact danger area does not need to carry out pressure relief work, but only needs to strengthen the detection of the target area. when the impact hazard level of the corresponding area is B, the roadway impact hazard area needs to be relieved and the detection of the target area should be strengthened; when the impact hazard level of the corresponding area is C, the personnel near the dangerous area shall be evacuated immediately, and the pressure relief work shall be carried out after a certain period of time; when the dangerous areas of roadway judged by different danger levels coincide, the higher one shall prevail.

Claims (7)

CONCLUSIESCONCLUSIONS 1. Een methode om het regionale impact risiconiveau en het lokale gebied te beoordelen op basis van de inversie van de seismische bronparameters, met als kenmerk dat de methode volgende stappen omvat: sl, een ruimtelijk coördinatensysteem voor het doelgebied opzetten en gegevens over aardbevingen in mijnen screenen; s2, een indexsysteem voor het evalueren van het impactrisiconiveau van het doelgebied opzetten, dat een sterkte-index Q en een dispersie-index H omvat; s3, de evaluatie-indexwaarde van het impactgevarenniveau van het doelgebied p-Y berekenen: H. s4, het impactrisiconiveau van het doelgebied volgens de evaluatie-indexwaarde D van het impactrisiconiveau van het doelgebied beoordelen, waarbij hoe groter de waarde D is, hoe hoger het risiconiveau is; s5, volgens de kortste afstand tussen het epicentrum en de tunnel en de relatie met de epicentrumradius het impactgevarengebied bepalen, en preventieve maatregelen volgens het gevarenniveau en het gebied van het doelgebied formuleren.1. A method to assess the regional impact risk level and the local area based on the inversion of the seismic source parameters, characterized in that the method includes the following steps: sl, setting up a spatial coordinate system for the target area and mining earthquake data screen; s2, set up an index system for evaluating the impact risk level of the target area, which includes a strength index Q and a dispersion index H; s3, calculate the evaluation index value of the impact hazard level of the target area p-Y: H. s4, evaluate the impact risk level of the target area according to the evaluation index value D of the impact risk level of the target area, where the larger the value D, the higher the risk level is; s5, according to the shortest distance between the epicenter and the tunnel and the relationship with the epicenter radius determine the impact hazard area, and formulate preventive measures according to the hazard level and the area of the target area. 2. De methode om het regionale impact risiconiveau en het lokale gebied te beoordelen op basis van de inversie van de seismische bronparameters volgens conclusie 1, met als kenmerk dat de stappen voor het vaststellen van het ruimtelijke coördinatensysteem van het doelgebied en het screenen van de mijnaardbevingsgegevens in sl zijn als volgt:2. The method of assessing the regional impact risk level and the local area based on the inversion of the seismic source parameters according to claim 1, characterized in that the steps of determining the spatial coordinate system of the target area and screening the mine earthquake data in sl are as follows: sl.1, een bron met een energie groter dan 100J als doelgebied selecteren, een brongegevens databank opzetten en analyseren;sl.1, select a source with an energy greater than 100J as the target area, set up and analyze a source data database; sl.2, volgens de ruimtelijke positie van het doelgebied, een referentiepunt 0%: 3920) selecteren als een oorsprong om een ruimtelijk rechthoekig coördinatensysteem Z vast te stellen;sl.2, according to the spatial position of the target area, select a reference point 0%: 3920) as an origin to establish a spatial rectangular coordinate system Z; s1.3. volgens de relatieve positie ten opzichte van het referentiepunt O, worden de coördinaten van alle seismische bronnen in de seismische brongegevensbank geconverteerd van wereldcoördinaten naar relatieve coördinaten, en de coördinaten me x,,V,,Z van elke seismische bron worden uitgedrukt als ( oY) ).s1.3. according to the relative position to the reference point O, the coordinates of all seismic sources in the seismic source database are converted from world coordinates to relative coordinates, and the coordinates me x,,V,,Z of each seismic source are expressed as ( oY) ). s1.4 volgens de ruimtelijke relatieve positie van de tunnel in het doelgebied en het referentiepunt o, wordt de geometrische uitdrukking van de tunnel in het xox, Yo, ITEs1.4 according to the spatial relative position of the tunnel in the target area and the reference point o, the geometric expression of the tunnel becomes xox, Yo, ITE X.,V.,Z, . coordinatensysteem Z weergegeven: ”! n P waar “1 is de ml (7) coördinaatcomponent van de coördinaat * °°’ 77 van elk punt aan het begin en einde van de tunnel in het coördinatensysteem Z, en, n;p is de ruimtelijke richtingsvector van de tunnel.X.,V.,Z, . coordinate system Z displayed: ”! n P where “1 is the ml(7) coordinate component of the coordinate * °°’ 77 of each point at the beginning and end of the tunnel in the Z coordinate system, and, n;p is the spatial direction vector of the tunnel. 3. De methode om het regionale impact risiconiveau en het lokale gebied te beoordelen op basis van de inversie van de seismische bronparameters volgens conclusie 1, met als kenmerk dat de methode voor het berekenen van de intensiteitsindex Q in stap s2 als volgt is: 1 ‚VJ > e/ 7 W/= J f V Q,= w‚W;/ dm ‚Vi < ef ij > 7 i=l €; waarin, Q de Q-waarde is van de intensiteitsindex van de j* bron in de bronbank, "i is het voorspellingsgewicht van de 7% bronparameter in de bronbank, Vi' is de i% bronparameterwaarde van de j* bron in de bronbank, waarbij 1, 2 en 3 respectievelijk seismisch moment, stralingsenergie en schijnbare spanning3. The method of assessing the regional impact risk level and local area based on the inversion of the seismic source parameters according to claim 1, characterized in that the method of calculating the intensity index Q in step s2 is as follows: 1. VJ > e/ 7 W/= J f V Q,= w‚W;/ dm ‚Vi < ef ij > 7 i=l €; where, Q is the Q value of the intensity index of the j* source in the source bank, "i is the prediction weight of the 7% source parameter in the source bank, Vi' is the i% source parameter value of the j* source in the source bank, where 1, 2 and 3 are seismic moment, radiant energy and apparent voltage respectively . Jo. . de . vertegenwoordigen, © is de gemiddelde waarde van de i bronparameter van de /® bron in de bronbank, en n is het totale aantal gebruikte bronparametertypen.. yo. . the . represent, © is the average value of the i resource parameter of the /® resource in the resource bank, and n is the total number of resource parameter types used. 4. De methode om het regionale impact risiconiveau en het lokale gebied te beoordelen op basis van de inversie van de seismische bronparameters volgens conclusie 1, met als kenmerk dat de methode voor het berekenen van de dispersie-index H in stap s2 is als volgt:The method of assessing the regional impact risk level and local area based on the inversion of the seismic source parameters according to claim 1, characterized in that the method of calculating the dispersion index H in step s2 is as follows: U= > wy; U>1 i=l j=1 H=y % U=LU<1U= > wy; U>1 i=l j=1 H=y % U=LU<1 Ww. . . . . daartussen, 7 1s het voorspellingsgewicht van de correlatie relatie van b My id di d .. derend ronparameters, is de spreidingsgraad tussen paarsgewijs corresponderende v—g;(x) Uy == bronparameters, die voldoet aan de relatie Umax waarbij x en y een reeks Sj) corresponderende bronparametertypen zijn, 7 is een set van aangepaste . . 4 . vq. . correlatie relaties van bronparameters, Ymax is een reeks spreidingslimieten tussen bronparameters, en n is de logaritme van de aangenomen bronparameters.WW. . . . . between them, 7 1s is the prediction weight of the correlation relation of b My id di d .. dend ron parameters, is the degree of dispersion between pairwise corresponding v—g;(x) Uy == source parameters, which satisfies the relation Umax where x and y have a set Sj) are corresponding source parameter types, 7 is a set of custom . . 4 . vq. . correlation relations of source parameters, Ymax is a set of spread limits between source parameters, and n is the logarithm of the assumed source parameters. 5. De methode om het regionale impact risiconiveau en het lokale gebied te beoordelen op basis van de inversie van de seismische bronparameters volgens conclusie 3, met als kenmerk dat de bronparameters omvatten seismisch moment, stralingsenergie, schijnbare spanning en bronstraal, die worden berekend door het Brune-model.The method of assessing the regional impact risk level and local area based on the inversion of the seismic source parameters according to claim 3, characterized in that the source parameters include seismic moment, radiant energy, apparent voltage and source radius, which are calculated by the Brune model. 6. De methode om het regionale impact risicoklasse en het lokale gebied te beoordelen op basis van de inversie van de seismische bronparameters volgens conclusie 3, met als kenmerk dat bij het oplossen van de bronparameters het . Co . nc A deeltjeszwermoptimalisatie-algoritme wordt toegepast om de nulfrequentielimiet op te lossen en de hoekfrequentie 4 en tot slot worden de bronparameters bepaald.The method of assessing the regional impact risk class and the local area based on the inversion of the seismic source parameters according to claim 3, characterized in that in solving the source parameters, the . Co . nc A particle swarm optimization algorithm is applied to solve for the zero frequency limit and the angular frequency 4 and finally the source parameters are determined. 7. De methode om het regionale impact risiconiveau en het lokale gebied te beoordelen op basis van de inversie van de seismische bronparameters volgens conclusie 1, met als kenmerk dat in stap s5, het specifieke schema voor het formuleren van preventieve maatregelen volgens het statistische regionale gevarenniveau en het lokale gebied als volgt is:7. The method of assessing the regional impact risk level and the local area based on the inversion of the seismic source parameters according to claim 1, characterized in that in step s5, the specific scheme for formulating preventive measures according to the statistical regional hazard level and the local area is as follows: wanneer het impactgevaarniveau van het overeenkomstige gebied A is, hoeft het gevarengebied van de tunnel geen drukontlastingswerkzaamheden uit te voeren, maar hoeft het alleen de detectie van het doelgebied te versterken;when the impact hazard level of the corresponding area is A, the tunnel hazard area does not need to carry out pressure relief work, but only needs to strengthen the detection of the target area; wanneer het impactgevaarniveau van het overeenkomstige gebied B is, moet het gevarengebied van de tunnel wel worden ontlast en tegelijkertijd moet de detectie van het doelgebied worden versterkt;when the impact hazard level of the corresponding area is B, the danger area of the tunnel must be relieved and at the same time the detection of the target area must be strengthened; wanneer het impactgevaarniveau van het overeenkomstige gebied C is, moet het personeel in de buurt van het gevaarlijke gebied onmiddellijk worden geëvacueerd en moeten de drukontlastingswerkzaamheden na een bepaalde tijdsperiode worden uitgevoerd;when the impact hazard level of the corresponding area is C, the personnel near the hazardous area should be evacuated immediately, and the pressure relief work should be carried out after a certain period of time; wanneer de gevaarlijke gebieden van de tunnel worden beoordeeld op basis van verschillende gevarenniveaus en de resultaten elkaar overlappen, dan wordt het gevarenniveau van dit gebied onderworpen aan de hoogst vastgestelde waarde.when the hazardous areas of the tunnel are assessed according to different hazard levels and the results overlap, the hazard level of this area is subject to the highest determined value.
NL2029589A 2021-11-02 2021-11-02 Identification Method of Regional Impact Risk Grade and Region Based on Inversion of Source Parameters NL2029589B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
NL2029589A NL2029589B1 (en) 2021-11-02 2021-11-02 Identification Method of Regional Impact Risk Grade and Region Based on Inversion of Source Parameters

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
NL2029589A NL2029589B1 (en) 2021-11-02 2021-11-02 Identification Method of Regional Impact Risk Grade and Region Based on Inversion of Source Parameters

Publications (1)

Publication Number Publication Date
NL2029589B1 true NL2029589B1 (en) 2023-06-01

Family

ID=86548173

Family Applications (1)

Application Number Title Priority Date Filing Date
NL2029589A NL2029589B1 (en) 2021-11-02 2021-11-02 Identification Method of Regional Impact Risk Grade and Region Based on Inversion of Source Parameters

Country Status (1)

Country Link
NL (1) NL2029589B1 (en)

Similar Documents

Publication Publication Date Title
US20230152479A1 (en) Physics-based and data-driven integrated method for rock burst hazard assessment
Wang et al. Risk Assessment of Water Inrush in Karst Tunnels Based on the Efficacy Coefficient Method.
Kabiesz et al. Application of rule-based models for seismic hazard prediction in coal mines.
Zhang et al. Risk identification for coal and gas outburst in underground coal mines: A critical review and future directions
Bakhtavar et al. Toward predicting blast-induced flyrock: a hybrid dimensional analysis fuzzy inference system
Wang et al. AdaBoost-driven multi-parameter real-time warning of rock burst risk in coal mines
Hao et al. Quantification of margins and uncertainties for the risk of water inrush in a karst tunnel: representations of epistemic uncertainty with probability
CN115577844A (en) Multi-parameter advanced prediction method for coal mine rock burst
CN116044501A (en) Advanced geological forecast dynamic monitoring and early warning system and method
Xiao et al. Hazard degree identification and coupling analysis of the influencing factors on goafs
Xu et al. Risk assessment of coal mine gas explosion based on cloud integrated similarity and fuzzy DEMATEL
CN112379425A (en) Regional impact hazard level and/30043territory identification method based on seismic source parameter inversion
Zhou et al. Evaluation of rock burst intensity based on annular grey target decision-making model with variable weight
NL2029589B1 (en) Identification Method of Regional Impact Risk Grade and Region Based on Inversion of Source Parameters
Chen et al. A new evaluation method for slope stability based on TOPSIS and MCS
Liu et al. Cluster analysis of moment tensor solutions and its application to rockburst risk assessment in underground coal mines
Zou et al. Intelligent Control of Smooth Blasting Quality in Rock Tunnels Using BP‐ANN, ENN, and ANFIS
Wang et al. Risk assessment and online forewarning of oil & gas storage and transportation facilities based on data mining
Zhang et al. Risk assessment of coal and gas outburst in driving face based on finite interval cloud model
Li et al. Research status and development trends of mine floor water inrush grade prediction
Xue et al. A method to predict rockburst using temporal trend test and its application
CN116070907A (en) Karst collapse susceptibility assessment method and system based on analytic hierarchy process
CN113374525A (en) Coal and gas outburst danger area identification comprehensive early warning method based on multi-parameter data fusion
Elmo et al. Can new technologies shake the empirical foundations of rock engineering?
Li et al. Risk Assessment of Building Foundation Pit Construction Based on Fuzzy Hierarchical Comprehensive Evaluation Method