CN111273361A - Earthquake monitoring table net special for coal mine - Google Patents
Earthquake monitoring table net special for coal mine Download PDFInfo
- Publication number
- CN111273361A CN111273361A CN202010144884.3A CN202010144884A CN111273361A CN 111273361 A CN111273361 A CN 111273361A CN 202010144884 A CN202010144884 A CN 202010144884A CN 111273361 A CN111273361 A CN 111273361A
- Authority
- CN
- China
- Prior art keywords
- time
- earthquake
- seismic
- ground
- underground
- 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.)
- Pending
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 50
- 239000003245 coal Substances 0.000 title claims abstract description 29
- 238000012806 monitoring device Methods 0.000 claims abstract description 25
- 238000012545 processing Methods 0.000 claims abstract description 14
- 238000005065 mining Methods 0.000 claims abstract description 9
- 238000009826 distribution Methods 0.000 claims abstract description 5
- 230000008859 change Effects 0.000 claims abstract description 4
- 238000000034 method Methods 0.000 claims description 18
- 230000001360 synchronised effect Effects 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 239000013307 optical fiber Substances 0.000 claims description 4
- 230000003993 interaction Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000009977 dual effect Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 230000002035 prolonged effect Effects 0.000 claims description 2
- 239000002131 composite material Substances 0.000 claims 1
- 239000011435 rock Substances 0.000 description 12
- 230000033001 locomotion Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 6
- 230000009471 action Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000013500 data storage Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 2
- 239000002243 precursor Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 208000010392 Bone Fractures Diseases 0.000 description 1
- 241000510672 Cuminum Species 0.000 description 1
- 235000007129 Cuminum cyminum Nutrition 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 206010017076 Fracture Diseases 0.000 description 1
- 208000013201 Stress fracture Diseases 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 230000002354 daily effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 210000001503 joint Anatomy 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/10—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices using induction coils
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
- E21F17/18—Special adaptations of signalling or alarm devices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/26—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
- G01V3/28—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device using induction coils
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Electromagnetism (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A special earthquake monitoring station network for coal mines comprises a ground monitoring device, an underground monitoring device, a GPS antenna, an NTP/PTP time synchronization server, a power supply, a ground monitoring server and a data processing computer; the ground monitoring device comprises a ground seismic data collector and a ground seismic sensor, and the underground monitoring device comprises an underground seismic data collector and an underground seismic sensor; the underground and ground seismic data acquisition devices are used for time service through a GPS antenna or an NTP/PTP time synchronization server connected with the underground and ground seismic data acquisition devices; the ground monitoring device and the underground monitoring device form a ground and underground special earthquake three-dimensional monitoring table network, the ground and underground combined arrangement mode is adopted, the earthquake source positioning precision is improved, the real-time online monitoring is carried out on the earthquake phenomenon in the mining area range, and the earthquake source distribution and the earthquake change trend are visually displayed.
Description
Technical Field
The invention relates to coal mine earthquake and rock burst monitoring, in particular to a special earthquake monitoring table network for a coal mine, and belongs to the technical field of mine safety and earthquake monitoring.
Background
Almost all major coal mining countries in the world are threatened by rock burst to varying degrees. The first coal mine rock burst phenomenon was reported in 1783 in the united kingdom worldwide. Since then, rock burst occurs in several dozen countries and regions such as the former soviet union, south africa, germany, united states, canada, india, and uk.
In China, rock burst occurs in the smoothing coal mine earlier than 1933. Later, with the increase of mining depth and the continuous expansion of mining range, rock burst phenomena occur in many mines in mine areas such as Beijing, Fushun, Shandong, Henan, Shanxi and the like. With the continuous increase of the mining depth, the rock burst disaster is more and more prominent, and the harm to the coal mine safety production is more and more large. Research shows that the dynamic disaster problem of coal mines such as rock burst and mine earthquake and the dynamic disaster instability problem of rock engineering are the results of instability of rock fracture processes such as initiation, development and penetration of micro fractures induced by stress field disturbance in the engineering activity process. Therefore, regardless of the kind of dynamic disaster of rock, there is mostly a precursor of micro-fracturing before the dynamic disaster occurs, and the direct cause of inducing micro-fracturing activity is an increase in stress or strain in the rock formation.
The scientific prediction and forecast of mine vibration is always the aim of cumin of students in various countries. Research in this area has made significant progress in recent years. At present, microseismic monitoring systems adopted in China mainly comprise the SOS, ARAMIS and ESG of Canada in Poland and KJ551 in China, and geophone monitoring systems mainly comprise the SAK, ARES-5/E and the like in Poland. The microseismic monitoring frequency band is generally 5-200 Hz, and the ground sound monitoring frequency band is 60-2000 Hz.
Because the impact mechanism, the precursor information and the vibration parameters are different at different stages, only one frequency band vibration monitoring device is equipped or the vibration monitoring data of a plurality of frequency bands are not subjected to seamless butt joint, so that an ideal monitoring effect is difficult to obtain. At present, mine vibration is mostly monitored by micro-vibration and earthquake sound, and seismic monitoring can realize lower-frequency, wider-area and more sensitive signal capture, so that the method has important significance for mine vibration monitoring.
Disclosure of Invention
The earthquake monitoring network special for the coal mine comprises a ground monitoring device, an underground monitoring device, a GPS antenna, an NTP/PTP time synchronization server, a power supply, a ground monitoring server and a data processing computer; the ground monitoring device comprises a ground seismic data collector and a ground seismic sensor, and the underground monitoring device comprises an underground seismic data collector and an underground seismic sensor; the underground seismic sensors are connected in parallel and are connected with an underground seismic data acquisition unit; the ground seismic sensors are connected in parallel and are connected with a ground seismic data acquisition unit; the underground and ground seismic data acquisition devices are used for time service through a GPS antenna or an NTP/PTP time synchronization server connected with the underground and ground seismic data acquisition devices; the underground and ground seismic data acquisition units are connected with a ground monitoring server and a data processing computer, and the data processing computer is provided with a station network real-time receiving software, a station network operation monitoring software, a station network data analysis software and a station network calibration processing software; the ground monitoring device and the underground monitoring device form a ground and underground special earthquake three-dimensional monitoring table network, the ground and underground combined arrangement mode is adopted, the earthquake source positioning precision is improved, the real-time online monitoring is carried out on the earthquake phenomenon in the mining area range, and the earthquake source distribution and the earthquake change trend are visually displayed.
The ground seismic sensor adopts a common seismometer, the underground seismic sensor is an explosion-proof type seismometer which accords with the coal mine safety standard, the seismometer adopts a three-component integrated electronic feedback mode, the seismometer is formed by integrally mounting three independent component velocity sensors (1 vertical UD, 1 east-west EW and 1 north-south NS), and the feedback structure expands the folding pendulum length of the seismometer through electronic feedback, so that the aim of prolonging the equivalent working period of the seismometer is fulfilled. The sensitivity of the sensor reaches 2000Vs/m, the dynamic range is better than 120dB, the observation frequency band of the short-period seismometer is 1S-50 Hz, the observation frequency band of the broadband seismometer is 60S-50 Hz, and the sensor is suitable for being installed under different ground and underground environments.
The ground seismic data acquisition unit and the underground seismic data acquisition unit are provided with high-speed analog signal input channels and can be connected with a three-component seismometer; by adopting GPS synchronous time service, on one hand, GPS clock signals output by the IRIG-B code generator can be received through optical fibers, and real-time data service is provided through the optical fibers, on the other hand, the whole system is subjected to network time service through the NTP/PTP time synchronization server, and the collector is internally provided with the NTP/PTP time service function to provide the real-time data function.
The seismic sensors and the seismic data collectors are respectively arranged on the ground and the underground to form a ground and underground special seismic three-dimensional monitoring table network, the ground and underground combined arrangement mode is adopted, the seismic source positioning precision is obviously improved, particularly, the vertical positioning depth has obvious advantages, the real-time online monitoring can be carried out on seismic phenomena such as earthquakes, mine earthquakes, micro-earthquakes and the like in a mining area range, and the seismic source distribution and the seismic variation trend are visually displayed.
The ground time service adopts GPS antenna and data acquisition unit cable conductor directly to link, and the GPS antenna of each ground station is independent, and signal reception is stable, and synchronous precision is high. The underground station data collector is internally provided with an NTP/PTP time service function, on one hand, each underground data collector is time-served through an optical cable through a ground time-service GPS antenna, the time-service precision is high, and the synchronous error is less than 10 mu s, on the other hand, the time-service is carried out through an NTP/PTP time synchronization server network, the ground GPS antenna receives time-service information, the underground data collector is time-served through the NTP/PTP time synchronization server connected in the network, and the time-service synchronous error is less than 1 ms.
The earthquake monitoring station network special for the coal mine comprises real-time receiving software, operation monitoring software, data analysis software and calibration processing software. Real-time receiving software, real-time receiving station network data, real-time displaying waveform, real-time monitoring of circuit on-off and real-time data storage; running monitoring software, monitoring the working state of the station network data processing software, and dynamically displaying the earthquake positioning result; data analysisThe software realizes earthquake waveform data man-machine interaction processing, earthquake phase identification, earthquake positioning and waveform file format conversion; and the calibration processing software is used for calculating pulse calibration and sine calibration, calculating a transfer function and detecting the working state of the seismometer. The seismic positioning method adopts grid positioning, and the grid positioning is carried out by using an objective function minimum search method, namely, a grid m (m) is searched in a full space range1,m2,m3). Different objective functions are established. Firstly, when the seismic source has the first arrival P wave and the second arrival S wave, according to the principle of the least square method, if the earthquake occurs at the position, the sum of the squares of the residual errors of the calculated travel time and the actual observed travel time of the seismic waves arriving at each station is minimum, and the seismic source position can be considered. Secondly, when only the first arrival P wave exists, according to the concept of equal time difference (EDT) surface, namely: t isi-Tj=ti(m1,m2,m3)-tj(m1,m2,m3) (ii) a That is to say the time T of a certain earthquake reaching the station i and ji、TjDifference sum grid (seismic source) (m)1,m2,m3) Theoretical travel time t to reach the two stationsi(m1,m2,m3)、tj(m1,m2,m3) The difference is equal. According to the least squares principle, this grid is the seismic source when the sum of the squared residuals of the difference between observed time and the difference between theoretical travel time is minimal.
The system adopts seismic duration magnitude MDThe magnitude of the vibration is calculated, the magnitude of the vibration is measured by recording the duration of the vibration, the magnitude of the vibration can better meet the actual situation than the magnitude of the vibration is measured by using the maximum amplitude, the duration tau is generally defined from the initial motion of the P wave, the time until the double amplitude of the wake wave is attenuated to 2 millimeters is traced backwards, and the wake wave more than 2 millimeters does not appear in 10 seconds later.
Description of the drawings:
in order to more clearly explain the technical solution of the present invention, the following description will be made of the necessary design details in the solution by means of schematic diagrams or schematic diagrams:
three-direction integrated structure of seismometer in figure 1
FIG. 2 is a schematic diagram of a feedback type
FIG. 3 is a block diagram of the working principle of the seismometer
FIG. 4 shows a structure diagram of a seismic monitoring network dedicated for coal mine
The specific implementation mode is as follows:
as shown in fig. 4, the earthquake monitoring network dedicated for coal mine comprises a ground monitoring device, an underground monitoring device, a GPS antenna, an NTP/PTP time synchronization server, a power supply, a ground monitoring server, and a data processing computer; the ground monitoring device comprises a ground seismic data collector and a ground seismic sensor, and the underground monitoring device comprises an underground seismic data collector and an underground seismic sensor; the underground seismic sensors are connected in parallel and are connected with an underground seismic data acquisition unit; the ground seismic sensors are connected in parallel and are connected with a ground seismic data acquisition unit; the underground and ground seismic data acquisition devices are used for time service through a GPS antenna or an NTP/PTP time synchronization server connected with the underground and ground seismic data acquisition devices; the underground and ground seismic data acquisition unit is connected with a ground monitoring server and a data processing computer, and the data processing computer is provided with a station network real-time receiving software, a station network operation monitoring software, a station network data analysis software and a station network calibration processing software.
The ground monitoring device and the underground monitoring device form a ground and underground special earthquake three-dimensional monitoring table network, the ground and underground combined arrangement mode is adopted, the earthquake source positioning precision is improved, the real-time online monitoring is carried out on the earthquake phenomenon in the mining area range, and the earthquake source distribution and the earthquake change trend are visually displayed.
The whole earthquake sensor adopts a three-direction integrated structure, and the sensitive axis directions of EW, NS and UD are pairwise orthogonal, so that a space measurement coordinate system is formed. The three-direction integrated structure of the seismometer is shown in the attached figure 1.
Furthermore, when the seismic sensor is installed, the geographical positions of the three components are adjusted, the EW points to the east and west, the NS points to the south and north, the UD points to the vertical direction, and the level bubble is adjusted to ensure that the seismic sensor is in the horizontal state and ensure the seismic waveform data accuracy.
Further, the whole earthquake sensor is of an electronic feedback type structure, namely, all directions are formed by additionally arranging electronic feedback circuits on a mechanical pendulum body in a certain form, the equivalent working period of the seismometer can be effectively prolonged, the mechanism of the electronic feedback type structure is shown in the attached drawing 2, a simple pendulum is taken as an example, if a certain movement A with a certain size exists on the ground, if the action of a feedback moment does not exist, the simple pendulum I should move to a position 1 under the action of the ground movement, but under the action of the feedback moment, the simple pendulum I only moves to a position 2, namely the position 2 is the output position of the feedback type simple pendulum I relative to the movement A on the ground. The output situation at this time is equivalent to the situation of a pure simple pendulum II, namely, under the action of ground motion A with the same size, the simple pendulum II and the feedback type simple pendulum I have the same rotation angle output. Obviously, the pure pendulum II has a longer physical pendulum length and a larger natural vibration period than the pure pendulum I. The folding pendulum length of the seismometer is expanded through electronic feedback, so that the purpose of extending the equivalent working period of the seismometer is achieved.
Further, the inside of the seismic sensor adopts an electromagnetic conversion principle, and a working principle block diagram of the seismometer is shown in fig. 3. When some motion (such as a mine earthquake) exists in a coal mine, the mass block in the seismic sensor generates relative motion relative to the suspended frame, so that the working coil generates relative motion relative to the magnetic field, and therefore induced voltage is generated in the working coil, the voltage is amplified by the amplifier and flows through the feedback network, and feedback torque is generated in the feedback coil in the form of current. Meanwhile, the voltage output from the amplifier passes through the compensation network and the differential driver to form double-end output of the seismometer.
Furthermore, the double ends of the seismometer output voltage signals, the seismometer is connected with a data acquisition unit for data acquisition, analog-to-digital conversion is carried out on the seismometer analog signals, and the data acquisition unit adopts general seismic data acquisition and recording equipment which has high resolution and large dynamic range, can output low-delay real-time data streams and is suitable for severe environments. The device can convert the input of a plurality of channels of analog voltage and frequency quantities into digital quantity to be output, has the functions of network and serial port data transmission, supports large-capacity data storage, and has the functions of data acquisition, recording, network data service, GPS synchronous time service and NTP/PTP time synchronous server network time service.
The seismic data acquisition device adopts GPS antenna timing, the GPS adopts IRIG-B code format synchronous timing, the IRIG-B code is a time code of one frame per second, and each pulse in the time format is called code element. The "on time" reference point for a symbol is its leading edge, the repetition rate of the symbol is referred to as the symbol rate, and the symbol rate of the B code is 100 pps. One index count for each symbol. The time interval between the leading edges of two adjacent code elements is the index counting interval, and the index counting interval of the B code is 10 ms. The index count starts at a frame reference point with a "0" and is incremented by 1 every other index count interval thereafter until the end of the frame. The index count interval of each frame of the B code is 100 until the end of the frame. The index counting interval of each frame of the B code is 100, and the index counting number is from 0 to 99. A time format frame reference mark begins. Consists of all symbols between two adjacent frame reference markers. The repetition rate of the time frame is the time frame rate, and the period of the time frame is the time frame period. The time frame rate of the B code is 1/sec, and the time frame period is 1 sec. The shorter the time frame period, the longer the information bits. The B code is 30 bits, with 10 bits (from 001 to 365 or 366) every day, 6 bits, 7 bits apart, 7 bits per second. The time sequence is second- > mi- > hour- > day. The positions are between P0 and P5. Pure binary seconds code of day time (SBS code): A. the B-format time code has 17 bits of pure binary second code of day time, 0 second in midnight, 86399 second time sequence at the maximum count, preceding low bit and succeeding high bit, besides the BCD code of year time.
Further, the seismic data acquisition unit adopts NTP/PTP time synchronization server network time service, the synchronization server network time service is suitable for a client/server program and a protocol which are operated on a computer, the program is written by a user serving as an NTP/PTP client, a server or both, under the basic condition, the NTP/PTP client sends out a time request to exchange time with the time server, and the exchange results are that the client can calculate the time delay, the compensation value of the time delay and adjust the time synchronization with the server. Typically, there are 6 exchanges in 5 to 10 minutes at the beginning of the setup. Once synchronized, a single information exchange is typically required to synchronize with the server time every 10 minutes. Redundant servers and different network paths are used to ensure accuracy of reliability.
The earthquake data collector is connected with the earthquake sensor well, acquires earthquake data information in real time, collects earthquake data of each station to the server through a special earthquake data transmission network built on the ground of the coal mine and in the pit, and the processor carries out daily data processing by accessing the data of the server.
The earthquake monitoring station network special for the coal mine comprises real-time receiving software, operation monitoring software, data analysis software and calibration processing software. Real-time receiving software, real-time receiving station network data, real-time displaying waveform, real-time monitoring of circuit on-off and real-time data storage; running monitoring software, monitoring the working state of the station network data processing software, and dynamically displaying the earthquake positioning result; the data analysis software is used for realizing the man-machine interaction processing, seismic facies recognition, seismic positioning and waveform file format conversion of seismic waveform data; and the calibration processing software is used for calculating pulse calibration and sine calibration, calculating a transfer function and detecting the working state of the seismometer.
The earthquake positioning method special for coal mine adopts grid positioning which is carried out by using an objective function minimum search method, namely, a grid m (m) is searched in a full space range1,m2,m3). Different objective functions are established.
Firstly, when the seismic source has the first arrival P wave and the second arrival S wave, according to the principle of the least square method, if the earthquake occurs at the position, the sum of the squares of the residual errors of the calculated travel time and the actual observed travel time of the seismic waves arriving at each station is minimum, and the seismic source position can be considered.
The objective function is as follows:
wherein: f (m)1,m2,m3T) is the residual square sum of the theoretical travel time and the observation travel time of each station and the center of the grid;
Obsiis the observed time of the seismic station; thoi(m1,m2,m3) The arrival time of the seismic station is calculated by a front database;
secondly, when only the first arrival P wave exists, according to the concept of equal time difference (EDT) surface, namely: t isi-Tj=ti(m1,m2,m3)-tj(m1,m2,m3) (ii) a That is to say the time T of a certain earthquake reaching the station i and ji、TjDifference sum grid (seismic source) (m)1,m2,m3) Theoretical travel time t to reach the two stationsi(m1,m2,m3)、tj(m1,m2,m3) The difference is equal; the objective function is established as follows:
where T and (m)1,m2,m3) Respectively the observed time and the theoretical travel time, T, of the earliest arriving stationiAnd ti(m1,m2,m3) Is the observed time and theoretical travel time of each station (m)1,m2,m3) Representing the position of a seismic source, wherein n is the number of stations; according to the least squares principle, this grid is the seismic source when the sum of the squared residuals of the difference between observed time and the difference between theoretical travel time is minimal.
The system adopts seismic duration magnitude MDCalculating magnitude of vibration, measuring magnitude of vibration by recording duration of vibration, better conforming to actual conditions than measuring magnitude of vibration by maximum amplitude, the duration tau being defined generally as the time from initial movement of P wave, backward pursuit until dual amplitude of wake wave decays to 2 mm, and no wake wave above 2 mm appearing in 10 seconds laterIn the range of the earthquake amplitude, tau is almost irrelevant to the earthquake centre distance, and the larger earthquake amplitude is not considered, so that the earthquake magnitude can be easily measured by a computer for a transmission table network. The method can be used as a station formal measuring and reporting MLThe method of (1). The general formula of the calculation is as follows:
MD=α0+α1logτ(s)+α2Δ
α 0, α 1, α 2 are constant in a certain range of magnitude, α 0 is-1.3-1.0, which varies with the magnification of the seismograph, α 1 is 1.7-2.6, which is mostly near 2.0 and varies with the passband of the seismograph, α 2 is very small, about ten thousandths, so it can be ignored when △ < 200 km.
Claims (10)
1. A special earthquake monitoring station network for coal mines comprises a ground monitoring device, an underground monitoring device, a GPS antenna, an NTP/PTP time synchronization server, a power supply, a ground monitoring server and a data processing computer; the ground monitoring device comprises a ground seismic data collector and a ground seismic sensor, and the underground monitoring device comprises an underground seismic data collector and an underground seismic sensor; the underground seismic sensors are connected in parallel and are connected with an underground seismic data acquisition unit; the ground seismic sensors are connected in parallel and are connected with a ground seismic data acquisition unit; the underground and ground seismic data acquisition devices are used for time service through a GPS antenna or an NTP/PTP time synchronization server connected with the underground and ground seismic data acquisition devices; the underground and ground seismic data acquisition units are connected with a ground monitoring server and a data processing computer, and the data processing computer is provided with a station network real-time receiving software, a station network operation monitoring software, a station network data analysis software and a station network calibration processing software; the method is characterized in that:
the ground monitoring device and the underground monitoring device form a ground and underground special earthquake three-dimensional monitoring table network, the ground and underground combined arrangement mode is adopted, the earthquake source positioning precision is improved, the real-time online monitoring is carried out on the earthquake phenomenon in the mining area range, and the earthquake source distribution and the earthquake change trend are visually displayed.
2. The special earthquake monitoring table network for coal mines as set forth in claim 1, which is characterized in that: the ground seismic sensor is a common seismometer, and the underground seismic sensor is an explosion-proof seismometer which is in a composite coal mine safety standard.
3. The special earthquake monitoring table network for the coal mine as claimed in claim 2, which is characterized in that: the whole seismic sensor adopts a three-direction integrated structure, and the sensitive axis directions of east-west, south-north and vertical direction are orthogonal pairwise, so that a space measurement coordinate system is formed.
4. The special earthquake monitoring table network for the coal mine as claimed in claim 3, which is characterized in that: the whole earthquake sensor is of an electronic feedback structure, namely, all directions are formed by additionally arranging electronic feedback circuits on a mechanical pendulum body in a certain form, so that the equivalent working period of the seismometer can be effectively prolonged.
5. The special earthquake monitoring table network for coal mines as set forth in claim 1, which is characterized in that: the ground seismic data acquisition unit and the underground seismic data acquisition unit are provided with high-speed analog signal input channels which are used for connecting seismic sensors; by adopting GPS synchronous time service, on one hand, GPS clock signals output from an IRIG-B code generator are received through optical fibers, and real-time data service is provided through the optical fibers, on the other hand, the whole system is subjected to network time service through an NTP/PTP time synchronization server, and the collector is internally provided with an NTP/PTP time service function to provide a real-time data function.
6. The special earthquake monitoring table network for the coal mine as claimed in claim 5, which is characterized in that: when the seismic data acquisition device adopts GPS antenna timing, the GPS adopts IRIG-B code format for synchronous timing; when the earthquake data collector adopts NTP/PTP time synchronization server network time service, the synchronization server network time service is suitable for a client/server program and a protocol which are operated on a computer, the program is compiled by users which are used as NTP/PTP clients, server terminals or both, under the basic condition, the NTP/PTP clients send out time requests to exchange time with the time server, and the exchange result is that the client can calculate the time delay, make up value of the time delay and adjust the time synchronization with the server.
7. The special earthquake monitoring table network for coal mines as set forth in claim 1, which is characterized in that: the real-time receiving software is used for receiving station network data in real time, displaying waveforms in real time, monitoring the on-off of a line in real time and storing real-time data; the monitoring software is operated to monitor the working state of the station network data processing software and dynamically display the earthquake positioning result; the data analysis software realizes earthquake waveform data man-machine interaction processing, earthquake phase identification, earthquake positioning and waveform file format conversion; the calibration processing software calculates pulse calibration and sine calibration, calculates a transfer function and detects the working state of the seismometer.
8. A method of locating seismic sources using the coal mine dedicated seismic monitoring station network of any of claims 1 to 7, characterized by: the positioning method adopts grid positioning which is carried out by using an objective function minimum search method and searches a grid m (m) in the full space range1,m2,m3) (ii) a Establishing different objective functions;
when the first arrival P wave and the second arrival S wave exist at the same time, according to the principle of the least square method, if an earthquake occurs at the position, the sum of squares of residual errors between the calculated travel time and the actual observed travel time of the earthquake waves reaching each station is minimum, and the calculated travel time and the actual observed travel time can be considered as the earthquake source position;
the objective function is as follows:
wherein: f (m)1,m2,m3T) is the residual square sum of the theoretical travel time and the observation travel time of each station and the center of the grid;
Obsiis the observed time of the seismic station; thoi(m1,m2,m3) The arrival time of the seismic station is calculated by a front database;
when only the first arrival P-wave, according to the concept of the equal time difference (EDT) plane, namely: t isi-Tj=ti(m1,m2,m3)-tj(m1,m2,m3) (ii) a That is to say the time T of a certain earthquake reaching the station i and ji、TjDifference sum grid (seismic source) (m)1,m2,m3) Theoretical travel time t to reach the two stationsi(m1,m2,m3)、tj(m1,m2,m3) The difference is equal; the objective function is established as follows:
where T and (m)1,m2,m3) Respectively the observed time and the theoretical travel time, T, of the earliest arriving stationiAnd ti(m1,m2,m3) Is the observed time and theoretical travel time of each station (m)1,m2,m3) Representing the position of a seismic source, wherein n is the number of stations; according to the least squares principle, this grid is the seismic source when the sum of the squared residuals of the difference between observed time and the difference between theoretical travel time is minimal.
9. A method for determining seismic magnitude using the coal mine dedicated seismic monitoring network of any one of claims 1 to 7, characterized by: using seismic duration magnitude MDMagnitude of the magnitude is calculated and the magnitude is measured by recording the duration of the vibration, duration τ being defined generally as the time from the onset of the P-wave, back-tracking until the dual amplitude of the wake decays to 2 mm, and no longer appearing above 2 mm in the next 10 seconds.
10. A method of determining seismic magnitude as claimed in claim 9, wherein: the calculation formula of the magnitude is as follows:
MD=α0+α1logτ(s)+α2Δ
α0、α1、α2constant within a certain range of magnitude α01.3-1.0, which varies with the magnification of the seismograph α11.7-2.6, mostly around 2.0, which varies with the passband of the seismograph α2Very small, on the order of parts per million, and therefore negligible at △ < 200 km.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010144884.3A CN111273361A (en) | 2020-03-04 | 2020-03-04 | Earthquake monitoring table net special for coal mine |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010144884.3A CN111273361A (en) | 2020-03-04 | 2020-03-04 | Earthquake monitoring table net special for coal mine |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111273361A true CN111273361A (en) | 2020-06-12 |
Family
ID=70997549
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010144884.3A Pending CN111273361A (en) | 2020-03-04 | 2020-03-04 | Earthquake monitoring table net special for coal mine |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111273361A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114198147A (en) * | 2021-11-16 | 2022-03-18 | 中国矿业大学(北京) | Coal mine rock burst multi-parameter monitoring system |
RU2818988C1 (en) * | 2023-06-19 | 2024-05-08 | Федеральное государственное бюджетное научное учреждение "Республиканский академический научно-исследовательский и проектно-конструкторский институт горной геологии, геомеханики, геофизики и маркшейдерского дела" (РАНИМИ) | Modular multichannel mine seismic station |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030132751A1 (en) * | 2002-01-15 | 2003-07-17 | Stolarczyk Larry G. | Radio-imaging of underground structures |
CN105510981A (en) * | 2015-12-04 | 2016-04-20 | 中国科学院地质与地球物理研究所 | Magnetite gob physical geography judgment method and device |
CN106483556A (en) * | 2016-10-09 | 2017-03-08 | 华北科技学院 | A kind of lasting earthquake magnitude based on Mine Earthquakes monitoring system and Richter scale conversion method |
CN106569255A (en) * | 2016-11-14 | 2017-04-19 | 中国矿业大学 | Wireless transmission-based coal mine stope earthquake source monitoring and locating system |
-
2020
- 2020-03-04 CN CN202010144884.3A patent/CN111273361A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030132751A1 (en) * | 2002-01-15 | 2003-07-17 | Stolarczyk Larry G. | Radio-imaging of underground structures |
CN105510981A (en) * | 2015-12-04 | 2016-04-20 | 中国科学院地质与地球物理研究所 | Magnetite gob physical geography judgment method and device |
CN106483556A (en) * | 2016-10-09 | 2017-03-08 | 华北科技学院 | A kind of lasting earthquake magnitude based on Mine Earthquakes monitoring system and Richter scale conversion method |
CN106569255A (en) * | 2016-11-14 | 2017-04-19 | 中国矿业大学 | Wireless transmission-based coal mine stope earthquake source monitoring and locating system |
Non-Patent Citations (2)
Title |
---|
刘逎泉: "《加速器理论》", 31 December 1990 * |
金星,等: "一种新地震定位方法研究", 《地震工程与工程振动》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114198147A (en) * | 2021-11-16 | 2022-03-18 | 中国矿业大学(北京) | Coal mine rock burst multi-parameter monitoring system |
RU2818988C1 (en) * | 2023-06-19 | 2024-05-08 | Федеральное государственное бюджетное научное учреждение "Республиканский академический научно-исследовательский и проектно-конструкторский институт горной геологии, геомеханики, геофизики и маркшейдерского дела" (РАНИМИ) | Modular multichannel mine seismic station |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bilich et al. | GPS seismology: Application to the 2002 M w 7.9 Denali fault earthquake | |
KR101179135B1 (en) | Apparatus and method for generating gps time | |
Kouba | Measuring seismic waves induced by large earthquakes with GPS | |
CN103471544B (en) | A kind of high precision displacement deformation monitoring application system based on the Big Dipper | |
US20240053431A1 (en) | Broadcast-type underwater navigation positioning system and method | |
Murray et al. | Regional global navigation satellite system networks for crustal deformation monitoring | |
CN101430384A (en) | Real-time accurate monitoring method for earth crust millimeter-level displacement | |
Shu et al. | Broadband velocities and displacements from integrated GPS and accelerometer data for high‐rate seismogeodesy | |
Butyrin et al. | Digital seismic logger “Ermak-5”. Experience of development and implementation | |
Samanes et al. | Estimating the VLF modal interference distance using the South America VLF Network (SAVNET) | |
CN111273361A (en) | Earthquake monitoring table net special for coal mine | |
CN103148845A (en) | Satellite combined compass and method for measuring and computing azimuth angle and pitch angle | |
Fang et al. | Epicenter and magnitude of large earthquake determined from high-rate GPS observations: A case study of the 2008 M 8.0 Wenchuan earthquake | |
CN103412338A (en) | Pseudo-random code electrical instrument | |
Dahmen et al. | Comparison and Combination of GNSS and Strong‐Motion Observations: A Case Study of the 2016 Mw 7.0 Kumamoto Earthquake | |
CN103941281A (en) | Distributed mine earthquake detection method and device | |
Kumar et al. | Earthquake early warning system for northern India | |
Melbourne et al. | 25‐second determination of 2019 M w 7.1 Ridgecrest earthquake coseismic deformation | |
JP2016128766A (en) | Method of computing pressure sensor output frequency, tsunami warning apparatus based on atmospheric pressure observation using the same, and tsunami warning system | |
LU504006B1 (en) | Interconnected microseismic monitoring system for mine water inrush disaster | |
Guo et al. | Rapid earthquake focal mechanism inversion using high-rate GPS velometers in sparse network | |
CN110146917A (en) | A kind of mining Microseismic monitoring system | |
Jiang et al. | Assessing Short‐Term Clock Errors and Drifts of Temporary Seismic Networks Using the Active Airgun Source in Binchuan, Yunnan | |
Doherty et al. | Timing potentials of Loran-C | |
Omerbashich | Detection and mapping of Earth body resonances with continuous GPS |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20201209 Address after: 271099 No.13, 3-3 / F, No.1 building, Lijing Garden Phase III, Daizong street, Tai'an City, Shandong Province Applicant after: SHANDONG SUCCEED MINING SAFETY ENGINEERING Co.,Ltd. Applicant after: NORTH CHINA INSTITUTE OF SCIENCE AND TECHNOLOGY Address before: 467 Xueyuan street, Sanhe City, Langfang City, Hebei Province Applicant before: NORTH CHINA INSTITUTE OF SCIENCE AND TECHNOLOGY |
|
TA01 | Transfer of patent application right |