CN113217109A - Waveform completion method for rockburst mine microseismic monitoring system - Google Patents

Waveform completion method for rockburst mine microseismic monitoring system Download PDF

Info

Publication number
CN113217109A
CN113217109A CN202110677344.6A CN202110677344A CN113217109A CN 113217109 A CN113217109 A CN 113217109A CN 202110677344 A CN202110677344 A CN 202110677344A CN 113217109 A CN113217109 A CN 113217109A
Authority
CN
China
Prior art keywords
time
waveform
vibration speed
vibration
monitoring system
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
CN202110677344.6A
Other languages
Chinese (zh)
Other versions
CN113217109B (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.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
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 China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN202110677344.6A priority Critical patent/CN113217109B/en
Publication of CN113217109A publication Critical patent/CN113217109A/en
Application granted granted Critical
Publication of CN113217109B publication Critical patent/CN113217109B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH 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
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Mining

Abstract

A method for completing a waveform of a microseismic monitoring system of a rock burst mine comprises the steps of firstly, obtaining a waveform missing interval; selecting the time t corresponding to the first trough before the missing waveform intervalmAnd the first acquisition point t before the missing waveform intervalMData point between and calculate the rate of change a of the vibration speed between two pointsmSelecting the time t corresponding to the first trough after the missing waveform intervalnAnd the first acquisition point t after the missing waveform intervalNData point between and calculate the rate of change a of the vibration speed between two pointsn(ii) a For tmTo tMAcceleration a betweenmFitting with t to obtain a functional relation am(t) analogously to obtain an(t)；am(t) and an(t) the integral over time t is the vibration velocity VmaxTo obtain tXAnd utilize am(t)、an(t) complementing the waveform missing part. The method can improve the application range and the application precision of the microseismic monitoring system and improve the microseismic monitoring and early warning level of the rock burst mine.

Description

Waveform completion method for rockburst mine microseismic monitoring system
Technical Field
The invention belongs to the field of coal rock dynamic disaster monitoring, and particularly relates to a waveform completion method for a micro-seismic monitoring system, in particular to a waveform completion method for a micro-seismic monitoring system of a rock burst mine.
Background
Rock burst is a typical mine dynamic phenomenon and has great harmfulness. Statistics show that 135 rock burst mines are mined in 13 provinces and autonomous regions in China at present, the rock burst mines tend to increase year by year, and rock burst accidents are more serious. The microseism monitoring method is a main means for monitoring dynamic phenomena of rock burst mines, the microseism technology is known as one of the most effective and most development potential monitoring methods for predicting and forecasting coal rock dynamic disasters, particularly coal mine dynamic disasters, and the microseism monitoring method is widely applied to rock burst mines in China.
The microseism sensor is generally a magnetoelectric speed sensor, and the working principle of the microseism sensor is that the speed quantity is converted into an electric signal by utilizing electromagnetic induction, and the microseism sensor is matched with a signal acquisition station, a signal recorder and the like of a microseism monitoring system to work and is used for measuring the coal wall vibration speed. The microseismic sensor has higher precision, such as an SOS microseismic monitoring system imported from Poland,the vibration speed range is generally 1 × 10-6～6.2×10-3m/s, the mine earthquake with the minimum energy of 100J can be monitored. For microseismic sensors, high precision means that the measuring range is small, so that most of rock burst mines are faced with the condition that the measuring range of the microseismic sensors is too small when a large-energy mine earthquake event occurs, and the monitored waveforms are incomplete, namely, the waveforms are lost.
At present, a corresponding method is not available in the field of microseismic monitoring to deal with the situation that the waveform acquired by a microseismic sensor is incomplete, so that the coal wall vibration speed is inaccurate, the mine earthquake energy is not accurately calculated, and the requirement of high-precision microseismic monitoring of a rock burst mine cannot be met.
Disclosure of Invention
The invention aims to provide a waveform completion method for a microseismic monitoring system of a rock burst mine, which can improve the application range and the application precision of the microseismic monitoring system and improve the microseismic monitoring and early warning level of the rock burst mine.
In order to achieve the purpose, the invention provides a waveform completion method of a rock burst mine microseismic monitoring system, which is based on the principle of least square method and comprises the following steps:
(1) obtaining waveform missing intervals according to waveform files acquired by a microseismic monitoring system and the maximum measuring range of a microseismic sensor, and naming each interval as C according to time sequence1、C2、C3、……、Cn；
In time order, C1The time and the vibration speed corresponding to the first wave trough are respectively tm1、vm1；C1The time and the vibration speed corresponding to the first wave trough are respectively tn1、vn1；C2The time and the vibration speed corresponding to the first wave trough are respectively tm2、vm2；C2The time and the vibration speed corresponding to the first wave trough are respectively tn2、vn2；C3The time and the vibration speed corresponding to the first wave trough are respectively tm3、vm3；C3The time and the vibration speed corresponding to the first wave trough are respectively tn3、vn3；……；CnThe time and the vibration speed corresponding to the first wave trough are respectively tmn、vmn；CnThe time and the vibration speed corresponding to the first wave trough are respectively tnn、vnn；
C1The time and the vibration speed of the previous acquisition point of the waveform loss are respectively tM1、vM1；C1The time and the vibration speed of the first acquisition point after the waveform is lost are respectively tN1、vN1；C2The time and the vibration speed of the previous acquisition point of the waveform loss are respectively tM2、vM2；C2The time and the vibration speed of the first acquisition point after the waveform is lost are respectively tN2、vN2；C3The time and the vibration speed of the previous acquisition point of the waveform loss are respectively tM3、vM3；C3The time and the vibration speed of the first acquisition point after the waveform is lost are respectively tN3、vN3；……；CnThe time and the vibration speed of the previous acquisition point of the waveform loss are respectively tMn、vMn；CnThe time and the vibration speed of the first acquisition point after the waveform is lost are respectively tNn、vNn；
(2) If the sampling frequency of the microseismic sensor is known to be f, the sampling interval delta t is 1/f; selecting tm1、tM1T isN1、tn1T ism2、tM2T isN2、tn2T ism3、tM3T isN3、tn3… …, tmn、tMnT isNn、tnnThe vibration acceleration a between two adjacent sampling points is obtained by calculation, and a is delta v/delta t according to tm1To tM1Calculating the intermediate data point to obtain tm1To tM1The value of the intermediate acceleration a with the time t; calculating to obtain t by the same methodN1、tn1T ism2、tM2T isN2、tn2T ism3、tM3T isN3、tn3… …, tmn、tMnT isNn、tnnThe value of the respective acceleration a with the time t;
(3) based on the principle of least square method, using the pairs tm1To tM1Acceleration a betweenm1Fitting with time t to obtain am1Functional relation a with tm1(t); by the same token obtain an1(t)、am2(t)、an2(t)、am3(t)、an3(t)、……、amn(t)、ann(t)；
(4) The principle that the peak vibration velocity after the completion of the waveform at the same time point is equal, namely am1(t) and an1(t) the peak velocities found are the same, function am1(t)、an1(t) the integral over time t is the peak vibration velocity V after the completion of the waveformmax1Is provided with Vmax1Corresponds to a time tx1To a, am1(t)、an1(t) the two functions are respectively integrated along the time axis to obtain an equation:
Figure BDA0003121429220000021
solve to obtain tx1(ii) a Will tx1Carry-in type
Figure BDA0003121429220000022
Or formula
Figure BDA0003121429220000023
Namely obtain
Figure BDA0003121429220000024
Simultaneously by the functional formula am1(t)、an1(t) completing the waveform missing interval;
in the same way, a is obtainedm2(t) and an2(t) between am3(t) and an3(t) between … …, amn(t) and ann(t) each of (V) corresponds tomax2，tx2)、(Vmax3，tx3)、……、(Vmaxn，txn) Simultaneously by functional formula am2(t)、an2(t)、am3(t)、an3(t)、……、amn(t)、ann(t) complementing the waveform-missing region.
Further, the fitting in step (3) is a piecewise fitting.
Preferably, the microseismic sensor in the step (1) is a microseismic sensor in an SOS microseismic monitoring system, and the sampling frequency f in the step (2) is 500 Hz.
The invention provides a waveform completion method based on the least square method, the principle is simple, the operability is strong, the calculation process involved in the method is easy to realize by programming, and the method is convenient to popularize and use; the method can fill the technical blank in the field of micro-seismic monitoring in the aspect of waveform completion, and can improve the use precision and the use range of a micro-seismic monitoring system; the method has important effects on accurately evaluating the mine earthquake risk and improving the monitoring and early warning efficiency of dynamic disasters such as rock burst and the like.
Drawings
FIG. 1 is a waveform of the microseismic system monitoring of the first embodiment, wherein C represents a waveform missing interval (t)m，vm)、(tn，vn) Respectively representing the time and the vibration speed of the first trough before and after the waveform missing interval;
FIG. 2 is a schematic diagram of the completion of the missing waveform interval in the first embodiment, VmaxIndicates the peak vibration velocity, t, after completion of the missing waveform sectionxRepresents VmaxA corresponding time; (t)M，vM)、(tN，vN) Respectively representing the time and the vibration speed of a first acquisition point before the waveform missing interval and after the waveform missing interval;
FIG. 3 is a fitting curve and a functional relation of the acceleration polynomial before the waveform is missing according to the first embodiment;
FIG. 4 is a fitting curve and a functional relation of the acceleration polynomial after the waveform is absent according to the first embodiment;
FIG. 5 is a microseismic system monitoring waveform of the second embodiment, C1、C2Indicates the waveform-missing region (t)m1，vm1)、(tn1，vn1) Respectively represent C1Wave form deletion interval front and C1Time and vibration speed at the first trough after the waveform missing interval; (t)m2，vm2)、(tn2，vn2) Respectively represent C2Wave form deletion interval front and C2Time and vibration speed at the first trough after the waveform missing interval;
FIG. 6 is a schematic diagram of the completion of the missing waveform interval according to the second embodiment, Vmax1、Vmax2Indicates the peak vibration velocity, t, after completion of the missing waveform sectionx1、tx2Represents Vmax1、Vmax2Respectively corresponding time; (t)M1，vM1)、(tN1，vN1) Respectively represent C1Wave form before deletion, C1Time and vibration speed of the first acquisition point after the waveform is lost; (t)M2，vM2)、(tN2，vN2) Respectively represent C2Wave form before deletion, C2Time and vibration speed of the first acquisition point after the waveform is lost;
FIG. 7 is C of the second embodiment1Fitting a curve and a function relation of an acceleration polynomial before waveform deletion;
FIG. 8 is C of the second embodiment1Linearly fitting a straight line and a function relation by the acceleration after the waveform is lost;
FIG. 9 shows a view of a second embodiment C2Fitting a curve and a function relation of an acceleration polynomial before waveform deletion;
FIG. 10 shows a view of a second embodiment C2And fitting a curve and a function relation of the acceleration polynomial after the waveform is absent.
Detailed Description
The invention is further described with reference to the following figures and specific examples.
Example one
As shown in figure 1, the invention discloses a waveform completion method of a microseismic monitoring system of a rock burst mine, which is based on the principle of least square method and comprises the following steps:
(1) waveform file and SOS micro-seismic monitor according to micro-seismic monitoring systemAcquiring a waveform missing interval at the maximum measuring range of a microseismic sensor in a measuring system, wherein C is the missing interval; according to the time sequence, the time and the vibration speed corresponding to the first wave trough before C are respectively tm、vm(ii) a The time and the vibration speed corresponding to the first trough after C are respectively tn、vn(ii) a The time and the vibration speed of the previous acquisition point of the C waveform lack are respectively tM、vM(ii) a The time and the vibration speed of the first acquisition point after the C waveform is lost are respectively tN、vN；
(2) The sampling frequency f of a microseismic sensor in the SOS microseismic monitoring system is 500Hz, and the sampling interval delta t is 0.002 s; as shown in fig. 2, t is selectedm、tMAnd tN、tnThe vibration acceleration a between two adjacent sampling points is obtained through calculation, and a is delta v/delta t;
according to tmTo tMInter data points and formula am＝(vm+1-vm) T is calculated asmTo tMThe value of the inter-acceleration a with time t is:
[(2.722,0.010488),(2.724,0.005776),(2.726,0.0076),(2.728,0.006992),(2.73,0.007144),(2.732,0.00684),(2.734,0.035872),(2.736,0.096064),(2.738,0.113848),(2.74,0.123576),(2.742,0.120688)]；
according to tNTo tnInter data points and formula an＝(vn+1-vn) T is calculated asnTo tNThe value of the inter-acceleration a with time t is:
[(2.766,-0.1022959915),(2.768,-0.1097440035),(2.77,-0.112631996808953),(2.772,-0.106551999124349),(2.774,-0.0948480010265485),(2.776,-0.071591988671571)]；
(3) due to tmTo tMAt time taWhere significant segmentation occurs, a segmentation fit is selected; based on the principle of least square method, for taTo tMAcceleration a betweenmFitting with time t to obtain amFunctional relation a with tm(t), the results are shown in FIG. 3; same principle and process as above for tnTo tNAcceleration a betweennFitting with time t to obtain anFunctional relation a with tn(t), the results are shown in FIG. 4;
(4) the principle that the peak vibration velocity after the completion of the waveform at the same time point is equal, namely am(t) and an(t) the peak velocity V obtainedmaxAre identical, function am(t)、an(t) the integral over time t is the peak vibration velocity V after the completion of the waveformmaxIs provided with VmaxCorresponds to a time txTo a, am(t)、an(t) the two functions are respectively integrated along the time axis to obtain an equation:
Figure BDA0003121429220000041
solve to obtain tx(ii) a Will txCarry-in type
Figure BDA0003121429220000042
Or formula
Figure BDA0003121429220000043
Can obtain
Figure BDA0003121429220000051
Simultaneously by the functional formula am(t)、an(t) complementing the waveform-missing region.
Example two
As shown in fig. 5, the waveform completion method for the microseismic monitoring system of the rock burst mine, which is disclosed by the invention, is based on the principle of the least square method and comprises the following steps:
(1) obtaining a waveform missing interval according to a waveform file of the microseismic monitoring system and the maximum measuring range of the SOS microseismic monitoring system, C1、C2Is a deletion interval;
in time order, C1The time and the vibration speed corresponding to the first wave trough are respectively tm1、vm1；C1The time and the vibration speed corresponding to the first wave trough are respectively tn1、vn1；C2First wave trough pairThe time and vibration speed are tm2、vm2；C2The time and the vibration speed corresponding to the first wave trough are respectively tn2、vn2；
C1The time and the vibration speed of the previous acquisition point of the waveform loss are respectively tM1、vM1；C1The time and the vibration speed of the first acquisition point after the waveform is lost are respectively tN1、vN1；C2The time and the vibration speed of the previous acquisition point of the waveform loss are respectively tM2、vM2；C2The time and the vibration speed of the first acquisition point after the waveform is lost are respectively tN2、vN2；
(2) The general sampling frequency f of the SOS microseismic monitoring system is 500Hz, and the sampling interval delta t is 0.002 s; as shown in FIG. 6, t is selectedm1、tM1T isN1、tn1T ism2、tM2T isN2、tn2The data points in between are obtained by obtaining the vibration acceleration a between two adjacent sampling points, wherein a is delta v/delta t;
according to tm1To tM1Inter data points and formula am1＝(vm1+1-vm1) T is calculated asm1To tM1The value of the inter-acceleration a with time t is:
[(4.148,0.00836),(4.15,0.008208),(4.152,0.008056),(4.154,0.008056),(4.156,0.007448),(4.158,0.0076),(4.16,0.007752),(4.162,0.007752),(4.164,0.0076),(4.166,0.00836),(4.168,0.008664),(4.17,0.008512),(4.172,0.008512),(4.174,0.00836),(4.176,0.008208),(4.178,0.007904),(4.18,0.007904),(4.182,0.007752),(4.184,0.0076),(4.186,0.007448),(4.188,0.007296),(4.19,0.00684),(4.192,0.007144),(4.194,0.00684),(4.196,0.006688),(4.198,0.006536),(4.2,0.006384),(4.202,0.006536),(4.204,0.00608),(4.206,0.012616),(4.208,0.019456),(4.21,0.019456),(4.212,0.019608),(4.214,0.018696),(4.216,0.01748),(4.218,0.016416),(4.22,0.01596),(4.222,0.016416),(4.224,0.017176),(4.226,0.018392),(4.228,0.019608),(4.23,0.020368),(4.232,0.020672),(4.234,0.020368),(4.236,0.020064),(4.238,0.020064),(4.24,0.021128),(4.242,0.022344)]；
according to tn1To tN1Inter data points and formula an1＝(vn1+1-vn1) T is calculated asn1To tN1The value of the inter-acceleration a with time t is:
[(4.29,-0.010336),(4.292,-0.010184),(4.294,-0.00988),(4.296,-0.00988),(4.298,-0.00988),(4.3,-0.009424),(4.302,-0.009424),(4.304,-0.009272),(4.306,-0.00912),(4.308,-0.008968),(4.31,-0.008968),(4.312,-0.008816),(4.314,-0.008512),(4.316,-0.015048),(4.318,-0.025688),(4.32,-0.014592),(4.322,-0.015808),(4.324,-0.014592),(4.326,-0.013984),(4.328,-0.013832),(4.33,-0.013984),(4.332,-0.014288),(4.334,-0.014896),(4.336,-0.0152),(4.338,-0.014896),(4.34,-0.01444),(4.342,-0.013832),(4.344,-0.013528),(4.346,-0.013072),(4.348,-0.013376),(4.35,-0.013832),(4.352,-0.014288),(4.354,-0.014896),(4.356,-0.0152),(4.358,-0.0152),(4.36,-0.015048),(4.362,-0.015352),(4.364,-0.015504),(4.366,-0.016264),(4.368,-0.017328),(4.37,-0.018696),(4.372,-0.019456),(4.374,-0.020064),(4.376,-0.020368),(4.378,-0.00532),(4.38,-0.000141)]；
according to tm2To tM2Inter data points and formula am2＝(vm2+1-vm2) T is calculated asm2To tM2The value of the inter-acceleration a with time t is:
[(4.42,0.010032),(4.422,0.010336),(4.424,0.010944),(4.426,0.010944),(4.428,0.010488),(4.43,0.010488),(4.432,0.010336),(4.434,0.010336),(4.436,0.009728),(4.438,0.009728),(4.44,0.009728),(4.442,0.009576),(4.444,0.009272),(4.446,0.00912),(4.448,0.008968),(4.45,0.008816),(4.452,0.008664),(4.454,0.008664),(4.456,0.008208),(4.458,0.008208),(4.46,0.008056),(4.462,0.008208),(4.464,0.007904),(4.466,0.0076),(4.468,0.007448),(4.47,0.007448),(4.472,0.007296),(4.474,0.007144),(4.476,0.006992),(4.478,0.006688),(4.48,0.00684),(4.482,0.00684),(4.484,0.006384),(4.486,0.006384),(4.488,0.006232),(4.49,0.006232),(4.492,0.005928),(4.494,0.005928),(4.496,0.005776),(4.498,0.005624),(4.5,0.005472),(4.502,0.006992),(4.504,0.034808),(4.506,0.044232),(4.508,0.041344),(4.51,0.040584),(4.512,0.03952),(4.514,0.000106)]；
according to tN2To tn2Inter data points and formula an2＝(vn2+1-vn2) T is calculated asN2To tn2The value of the inter-acceleration a with time t is:
[(4.568,-0.01064),(4.57,-0.01064),(4.572,-0.010336),(4.574,-0.010184),(4.576,-0.010184),(4.578,-0.00988),(4.58,-0.00988),(4.582,-0.009576),(4.584,-0.009576),(4.586,-0.009272),(4.588,-0.009272),(4.59,-0.00912),(4.592,-0.008968),(4.594,-0.008968),(4.596,-0.008664),(4.598,-0.008664),(4.6,-0.008512),(4.602,-0.008208),(4.604,-0.008208),(4.606,-0.008056),(4.608,-0.007904),(4.61,-0.007904),(4.612,-0.007752),(4.614,-0.0076),(4.616,-0.007448),(4.618,-0.007448),(4.62,-0.007144),(4.622,-0.006992),(4.624,-0.006992),(4.626,-0.00684),(4.628,-0.007904),(4.63,-0.033896),(4.632,-0.033744),(4.634,-0.019608),(4.636,-0.0228),(4.638,-0.022192),(4.64,-0.02204),(4.642,-0.021736),(4.644,-0.021584),(4.646,-0.020824),(4.648,-0.020368),(4.65,-0.019912),(4.652,-0.019152),(4.654,-0.018392),(4.656,-0.01748),(4.658,-0.016568),(4.66,-0.015352),(4.662,-0.002584),(4.664,-0.000133)]；
(3) due to tm1To tM1At time ta1Where significant segmentation occurs, a segmentation fit is selected; based on the principle of least square method, for ta1To tM1Acceleration a betweenm1Fitting with time t to obtain am1Functional relation a with tm1(t), the results are shown in FIG. 7;
same principle and process as above for tb1To tN1Acceleration a betweenn1Fitting with time t to obtain an1Functional relation a with tn1(t), the results are shown in FIG. 8;
same principle and process as above for ta2To tM2Acceleration a betweenm2Fitting with time t to obtain am2Functional relation a with tm2(t), the results are shown in FIG. 9;
same principle and process as above for tb2To tN2Acceleration a betweenn2Fitting with time t to obtain an2Functional relation a with tn2(t), the results are shown in FIG. 10;
(4) the principle that the peak vibration velocity after the completion of the waveform at the same time point is equal, namely am1(t) and an1(t) the peak velocity V obtainedmax1Are identical, function am1(t)、an1(t) the integral over time t is the peak vibration velocity V after the completion of the waveformmax1Is provided with Vmax1Corresponds to a time tx1To a, am1(t)、an1(t) the two functions are integrated along the time axis respectively to obtain the equation:
Figure BDA0003121429220000071
solving the available tx1(ii) a Will tx1Carry-in type
Figure BDA0003121429220000072
Or formula
Figure BDA0003121429220000073
Can obtain
Figure BDA0003121429220000074
Simultaneously by the functional formula am1(t)、an1(t) completing the waveform of the missing section;
the principle that the peak vibration velocity after the completion of the waveform at the same time point is equal, namely am2(t) and an2(t) the peak velocity V obtainedmax2Are identical, function am2(t)、an2(t) the integral over time t is the peak vibration velocity V after the completion of the waveformmax2Is provided with Vmax2Corresponds to a time tx2To a, am2(t)、an2(t) the two functions are integrated along the time axis respectively to obtain the equation:
Figure BDA0003121429220000075
solving the available tx2(ii) a Will tx2Carry-in type
Figure BDA0003121429220000076
Or formula
Figure BDA0003121429220000077
Can obtain
Figure BDA0003121429220000078
Simultaneously by the functional formula am2(t)、an2(t) filling the waveform of the missing section.
The first embodiment and the second embodiment are respectively a waveform complementing method in which the waveform missing interval is one and two, and when there are a plurality of waveform missing intervals, the waveform complementing method for the plurality of waveform missing intervals is the same as that in the second embodiment.

Claims (3)

1. A waveform completion method for a rock burst mine microseismic monitoring system is characterized by being based on the principle of least square method and comprising the following steps:
(1) obtaining waveform missing intervals according to waveform files acquired by a microseismic monitoring system and the maximum measuring range of a microseismic sensor, and naming each interval as C according to time sequence1、C2、C3、……、Cn；
In time order, C1The time and the vibration speed corresponding to the first wave trough are respectively tm1、vm1；C1The time and the vibration speed corresponding to the first wave trough are respectively tn1、vn1；C2The time and the vibration speed corresponding to the first wave trough are respectively tm2、vm2；C2The time and the vibration speed corresponding to the first wave trough are respectively tn2、vn2；C3The time and the vibration speed corresponding to the first wave trough are respectively tm3、vm3；C3The time and the vibration speed corresponding to the first wave trough are respectively tn3、vn3；……；CnThe time and the vibration speed corresponding to the first wave trough are respectively tmn、vmn；CnThe time and the vibration speed corresponding to the first wave trough are respectively tnn、vnn；
C1The time and the vibration speed of the previous acquisition point of the waveform loss are respectively tM1、vM1；C1The time and the vibration speed of the first acquisition point after the waveform is lost are respectively tN1、vN1；C2The time and the vibration speed of the previous acquisition point of the waveform loss are respectively tM2、vM2；C2The time and the vibration speed of the first acquisition point after the waveform is lost are respectively tN2、vN2；C3The time and the vibration speed of the previous acquisition point of the waveform loss are respectively tM3、vM3；C3The time and the vibration speed of the first acquisition point after the waveform is lost are respectively tN3、vN3；……；CnThe time and the vibration speed of the previous acquisition point of the waveform loss are respectively tMn、vMn；CnThe time and the vibration speed of the first acquisition point after the waveform is lost are respectively tNn、vNn；
(2) If the sampling frequency of the microseismic sensor is known to be f, the sampling interval delta t is 1/f; selecting tm1、tM1T isN1、tn1T ism2、tM2T isN2、tn2T ism3、tM3T isN3、tn3… …, tmn、tMnT isNn、tnnThe vibration acceleration a between two adjacent sampling points is obtained by calculation, and a is delta v/delta t according to tm1To tM1Calculating the intermediate data point to obtain tm1To tM1The value of the intermediate acceleration a with the time t; calculating to obtain t by the same methodN1、tn1T ism2、tM2T isN2、tn2T ism3、tM3T isN3、tn3… …, tmn、tMnBetween，tNn、tnnThe value of the respective acceleration a with the time t;
(3) based on the principle of least square method, using the pairs tm1To tM1Acceleration a betweenm1Fitting with time t to obtain am1Functional relation a with tm1(t); by the same token obtain an1(t)、am2(t)、an2(t)、am3(t)、an3(t)、……、amn(t)、ann(t)；
(4) The principle that the peak vibration velocity after the completion of the waveform at the same time point is equal, namely am1(t) and an1(t) the peak velocities found are the same, function am1(t)、an1(t) the integral over time t is the peak vibration velocity V after the completion of the waveformmax1Is provided with Vmax1Corresponds to a time tx1To a, am1(t)、an1(t) the two functions are respectively integrated along the time axis to obtain an equation:
Figure FDA0003121429210000021
solve to obtain tx1(ii) a Will tx1Carry-in type
Figure FDA0003121429210000022
Or formula
Figure FDA0003121429210000023
Namely obtain
Figure FDA0003121429210000024
Simultaneously by the functional formula am1(t)、an1(t) completing the waveform missing interval;
in the same way, a is obtainedm2(t) and an2(t) between am3(t) and an3(t) between … …, amn(t) and ann(t) each of (V) corresponds tomax2，tx2)、(Vmax3，tx3)、……、(Vmaxn，txn) Simultaneously by functional formula am2(t)、an2(t)、am3(t)、an3(t)、……、amn(t)、ann(t) complementing the waveform-missing region.
2. The method for completing the waveform of the microseismic monitoring system for rock burst mine according to claim 1, wherein the fitting in the step (3) is a piecewise fitting.
3. The method for completing the waveform of the microseismic monitoring system of the rock burst mine according to claim 1 or 2, wherein the microseismic sensor in the step (1) is the microseismic sensor in the SOS microseismic monitoring system, and the sampling frequency f in the step (2) is 500 Hz.
CN202110677344.6A 2021-06-18 2021-06-18 Waveform completion method for rockburst mine microseismic monitoring system Active CN113217109B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110677344.6A CN113217109B (en) 2021-06-18 2021-06-18 Waveform completion method for rockburst mine microseismic monitoring system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110677344.6A CN113217109B (en) 2021-06-18 2021-06-18 Waveform completion method for rockburst mine microseismic monitoring system

Publications (2)

Publication Number Publication Date
CN113217109A true CN113217109A (en) 2021-08-06
CN113217109B CN113217109B (en) 2022-06-10

Family

ID=77080487

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110677344.6A Active CN113217109B (en) 2021-06-18 2021-06-18 Waveform completion method for rockburst mine microseismic monitoring system

Country Status (1)

Country Link
CN (1) CN113217109B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114895353A (en) * 2022-05-27 2022-08-12 中国矿业大学 Time service alignment method for data collected by monitoring unit of well-ground integrated microseismic monitoring system

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5598508A (en) * 1991-10-18 1997-01-28 Goldman; Julian M. Real-time waveform analysis using artificial neural networks
CN105118443A (en) * 2015-06-25 2015-12-02 华南师范大学 Electrophoresis display device driving method and electrophoresis display device
CN107505839A (en) * 2017-07-21 2017-12-22 中国地质大学（武汉） A kind of synchronous waveform method and system of virtual instrument
CN107705848A (en) * 2017-11-07 2018-02-16 杭州八脉科技有限公司 A kind of method and system for recommending coordinating program according to user health situation
CN107967240A (en) * 2017-11-23 2018-04-27 大连理工大学 Shatter-proof time-histories acceleration optimization algorithm based on artificial ground shaking generating process
CN110568486A (en) * 2019-09-17 2019-12-13 电子科技大学 Seismic signal completion method based on synchronous sparse low-rank tensor completion model
CN111953380A (en) * 2020-07-02 2020-11-17 中国人民解放军战略支援部队信息工程大学 Non-periodic long code direct sequence spread spectrum signal time delay estimation method and system based on norm fitting

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5598508A (en) * 1991-10-18 1997-01-28 Goldman; Julian M. Real-time waveform analysis using artificial neural networks
CN105118443A (en) * 2015-06-25 2015-12-02 华南师范大学 Electrophoresis display device driving method and electrophoresis display device
CN107505839A (en) * 2017-07-21 2017-12-22 中国地质大学（武汉） A kind of synchronous waveform method and system of virtual instrument
CN107705848A (en) * 2017-11-07 2018-02-16 杭州八脉科技有限公司 A kind of method and system for recommending coordinating program according to user health situation
CN107967240A (en) * 2017-11-23 2018-04-27 大连理工大学 Shatter-proof time-histories acceleration optimization algorithm based on artificial ground shaking generating process
CN110568486A (en) * 2019-09-17 2019-12-13 电子科技大学 Seismic signal completion method based on synchronous sparse low-rank tensor completion model
CN111953380A (en) * 2020-07-02 2020-11-17 中国人民解放军战略支援部队信息工程大学 Non-periodic long code direct sequence spread spectrum signal time delay estimation method and system based on norm fitting

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
徐凯等: "基于随机森林方法的地震插值方法研究", 《石油科学通报》 *
王立夫等: "基于一致性稀疏表示的地震信号补全算法", 《计算机应用》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114895353A (en) * 2022-05-27 2022-08-12 中国矿业大学 Time service alignment method for data collected by monitoring unit of well-ground integrated microseismic monitoring system
CN114895353B (en) * 2022-05-27 2023-03-10 中国矿业大学 Time service alignment method for data collected by monitoring unit of well-ground integrated microseismic monitoring system
US11719841B1 (en) 2022-05-27 2023-08-08 China University Of Mining And Technology Timing alignment method for data acquired by monitoring units of borehole-surface micro-seismic monitoring system

Also Published As

Publication number Publication date
CN113217109B (en) 2022-06-10

Similar Documents

Publication Publication Date Title
CN105022031B (en) A kind of layered velocity localization method of region rock mass microseism focus
CN102435980B (en) Analytical solution-based acoustic emission source or micro seismic source positioning method
CN103697999B (en) A kind of heavily stressed hard rock TBM construction tunnel microseism velocity of wave real time acquiring method
CN102262220B (en) Positioning method of micro-seismic source or acoustic emission source based on non-linear fitting
CN103018788A (en) Advanced detection device and method for unfavorable geology and rock mass mechanical properties of deep and long tunnels
CN108490485B (en) Double-track tunnel microseismic monitoring system and positioning accuracy assessment method thereof
CN202975361U (en) Advance detection apparatus for unfavorable geology and rock mass mechanical property of deep and long tunnel
CN106437843B (en) coal mine bottom plate water guide channel identification method based on microseismic monitoring
CN104453879B (en) The Forecasting Methodology of pressure before boring
CN101363916A (en) High precision refraction static correction data inversion method
CN113217109B (en) Waveform completion method for rockburst mine microseismic monitoring system
CN105258662A (en) Method for measuring end face spatial displacement and angular variation of shafting engineering member on the basis of stay-supported displacement sensors
CN108226288A (en) A kind of subway tunnel railway roadbed comes to nothing quantity monitoring method
CN103235334B (en) Low drainage roadway fault-throw advanced-detection method based on seismic signals
CN102608223B (en) Digital signal analysis method based on spectral energy ratio
CN104459824A (en) Device for monitoring fracturing effect of micro-earthquakes and monitoring method with device
CN105022091B (en) The far field focus method for rapidly positioning that a kind of nothing tests the speed in advance
CN103606240A (en) Method for fire alarming through distributed type fiber optical temperature sensor system
CN109239779B (en) Testing method for loose circle of tunnel surrounding rock and surrounding rock damage grading method
CN114415116A (en) Coal mining monitoring method and device and computer equipment
CN112130204B (en) On-site testing method suitable for wave velocity of surrounding rock of deep-buried tunnel
CN102096093B (en) Method for calculating propagation speed of earthquake waves in mine region by utilizing microearthquake point as earthquake source
CN103742156A (en) Method for determining timing and modes for changing opposite advancing into unidirectional advancing before deep hard rock tunnel holing-through
CN207232405U (en) The microseism information monitoring system of rock partition wall is pressed from both sides in a kind of freeway tunnel
CN109521221B (en) Method for acquiring microwave wave velocity of hard rock tunnel constructed by drilling and blasting method in real time

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant