CN109708550B - Blind gun identification method based on blasting vibration signal detection - Google Patents

Blind gun identification method based on blasting vibration signal detection Download PDF

Info

Publication number
CN109708550B
CN109708550B CN201910022593.4A CN201910022593A CN109708550B CN 109708550 B CN109708550 B CN 109708550B CN 201910022593 A CN201910022593 A CN 201910022593A CN 109708550 B CN109708550 B CN 109708550B
Authority
CN
China
Prior art keywords
blasting
blasting vibration
vibration signal
blind
component
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.)
Active
Application number
CN201910022593.4A
Other languages
Chinese (zh)
Other versions
CN109708550A (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.)
Zhejiang Institute of Hydraulics and Estuary
Zhejiang Guangchuan Engineering Consulting Co Ltd
Original Assignee
Zhejiang Institute of Hydraulics and Estuary
Zhejiang Guangchuan Engineering Consulting 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 Zhejiang Institute of Hydraulics and Estuary, Zhejiang Guangchuan Engineering Consulting Co Ltd filed Critical Zhejiang Institute of Hydraulics and Estuary
Priority to CN201910022593.4A priority Critical patent/CN109708550B/en
Publication of CN109708550A publication Critical patent/CN109708550A/en
Application granted granted Critical
Publication of CN109708550B publication Critical patent/CN109708550B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention relates to an engineering blasting technology, in particular to a blind shot identification method based on blasting vibration signal detection. The invention relates to a blind shot identification method based on blasting vibration signal detection, which comprises the following steps: (1) obtaining and recording an original blasting vibration signal f (t); (2) analyzing the original blasting vibration signal f (t) obtained in the step (1) by empirical mode decomposition to obtain an IMF component of the blasting vibration signal and obtain a principal component of the IMF component; (3) acquiring an envelope function of a principal vibration component; (4) confirming whether a blind shot is generated; (5) and determining the position of the blind shot. The invention provides a quick, economic and reliable blind shot identification method based on blasting vibration signal detection, so as to avoid safety accidents caused by inaccurate and timely blind shot elimination.

Description

Blind gun identification method based on blasting vibration signal detection
Technical Field
The invention relates to an engineering blasting technology, in particular to a blind shot identification method based on blasting vibration signal detection, which can be used for blind shot identification in a differential hole-by-hole detonating network.
Background
The one-time dosage of the engineering blasting is larger, and in order to achieve good blasting effect and reduce blasting vibration, the general sectional blasting is adopted and differential delay detonation is adopted. During the blasting process, blind shots may be generated during the blasting process due to some uncertain factors. Accidents are easily caused after blind shots are generated by blasting, and the traditional blind shot identification method is mainly judged by the experience of workers, so that the actual situation is often difficult to accurately reflect by the result. Therefore, the method for identifying the blind cannon is practical, rapid and reliable, and has great practical significance.
Disclosure of Invention
Aiming at the technical current situation of blind shot identification in the current blasting engineering, the invention aims to provide a quick, economic and reliable blind shot identification method so as to avoid safety accidents caused by inaccurate and timely blind shot elimination.
Aiming at the aim of the invention, the blind shot identification method based on the blasting vibration signal detection mainly comprises the following steps:
the blind shot identification method based on the blasting vibration signal detection comprises the following steps:
(1) establishing a hole-by-hole detonation type blasting network, arranging a blasting vibration tester at a position for recording blasting vibration signals, and recording original blasting vibration signals f (t);
(2) analyzing the original blasting vibration signal f (t) obtained in the step (1) by empirical mode decomposition to obtain an IMF component of the blasting vibration signal and obtain a principal component of the IMF component;
(3) performing Hilbert transformation on the principal component of the IMF component obtained in the step (2) to obtain an envelope function of the principal vibration component;
(4) analyzing the envelope function obtained in the step (3), wherein the number of each pole on the envelope function is the number of each actual initiation section of the blasting network, comparing the number with the number of theoretical sections, and if the number is equal, generating no blind shot; if the number is smaller than the theoretical segment number, generating blind shots;
(5) and (4) obtaining the actual detonation time of each section of explosive according to the envelope function obtained in the step (4), sequentially subtracting the detonation times of the adjacent sections to obtain the actual detonation delay time of the adjacent blast holes, and comparing the actual detonation delay time with the theoretical delay time of blasting design so as to determine the position of the blind shot.
The blasting vibration tester is a TC-4850 blasting vibration tester.
The specific steps of the step (2) are as follows:
(1) using cubic spline interpolation to form an upper envelope line and a lower envelope line of all local extreme points of an original blasting vibration signal f (t), wherein the upper envelope line and the lower envelope line comprise all data points of the signal;
(2) the average value of the upper envelope line and the lower envelope line obtained in the step (1) is recorded as m1(t) obtaining
Figure DEST_PATH_IMAGE001
(3) For different blasting vibration signals f (t), if h1(t) obtained in the step (2) is an IMF component or not, if h1(t) is an IMF component, the result is marked as c1(t) = h1 (t); if h1(t) is not an IMF component, taking the h1(t) signal obtained in the step (2) as an original signal, repeating the steps (1) to (2) k times to obtain the k-th screening data h1k (t):
Figure 476330DEST_PATH_IMAGE002
h1k (t) is not an IMF component, and there must be a termination condition for the screening process, determined by calculating the standard deviation SD between two successive processing results:
Figure DEST_PATH_IMAGE003
when the value of SD is 0.1-0.4, the first satisfied h1k (t) is the first order, and c1(t) = h1k (t), then c1(t) is the first-order IMF component of the signal f (t);
(4) decomposing the 1 st IMF component c obtained in the step (2)1(t) removing the original blasting vibration signal f (t) to obtain a residual time sequence r1(t):
Figure 455787DEST_PATH_IMAGE004
(5) When r is1(t) is a monotonic function or only one pole exists, empirical mode decomposition is completed, and r is1(t) recording as residual error, otherwise, recording r obtained in the step (4)1(t) repeating steps (1) through (4) as raw data to obtain other order components of the signal, such as c1(t)~ ci(t)。
The value of SD is preferably 0.2-0.3.
Preference is given to c2(t),c3(t) or c4(t) is a principal component of the IMF component. Comparing principal component components of Intrinsic Mode Function (IMF) components, and selecting IMF components capable of reflecting principal components of signal information, wherein the obtaining principle is as follows: firstly, comparing the IMF component with the original signal characteristics, such as the approximate time point of waveform change, waveform form and the like (generally, the 1 st IMF component is a noise component, and the 2 nd, 3 rd and 4 th components are main characteristic components); observing the waveform information retention degree of the selected component, wherein the later component loses the original waveform information due to more screening times; and comparing whether the effect is optimal or not.
The specific steps of the step (3) are as follows:
hilbert transform is performed on principal component c (t) of the selected IMF components,
Figure DEST_PATH_IMAGE005
wherein H [ c (t) ] is Hilbert transform function of c (t), PV is Cauchy's principal value, and c (t) is selected IMF principal component;
Figure 861623DEST_PATH_IMAGE006
wherein a (t) is the envelope function of the IMF principal component c (t).
The blind shot identification method based on the blasting vibration signal detection can effectively avoid a method for judging blind shots by manual experience, and can more effectively identify the blind shots and the positions of the blind shots.
Drawings
FIG. 1 is a schematic flow chart of a blind shot identification method based on blasting vibration signal detection according to the embodiment;
FIG. 2 is a diagram of a blasting network of the blind shot identification method based on the detection of the blasting vibration signal according to the embodiment;
fig. 3 is an original blasting vibration signal f (t) of the blind shot identification method based on the blasting vibration signal detection in the embodiment;
FIG. 4 shows the components (i.e., c) of signals IMF 1-IMF 13 of the blind shot identification method based on the detection of the blasting vibration signal in the present embodiment1(t)~ ci(t),c1(t) corresponding to IMF1, c2(t) corresponds to IMF2, and so on) and remainder r;
FIG. 5 shows the IMF5 components of the main components when the blind shot identification method based on the detection of the blasting vibration signal is used without a blind shot;
FIG. 6 shows the IMF5 component Hilbert transform modulus value when the blind shot identification method based on the blast vibration signal detection has no blind shot in the embodiment;
fig. 7 shows the Hilbert transform module values of the IMF principal component components when blind shots exist in the blind shot identification method based on the blast vibration signal detection according to the embodiment.
Detailed Description
As further described below with reference to fig. 1 to 6, a blind shot identification method based on blast vibration signal detection includes the following steps (fig. 1 is a schematic flow chart of the present application):
(1) establishing a hole-by-hole detonation type blasting network (established according to the figure 2), setting a blasting vibration tester at a position for recording a blasting vibration signal, and recording an original blasting vibration signal f (t), wherein f (t) is as shown in figure 3 (the embodiment adopts a TC-4850 blasting vibration tester);
(2) analyzing the original blasting vibration signal f (t) obtained in the step (1) by empirical mode decomposition to obtain an IMF component of the blasting vibration signal and obtain a principal component of the IMF component; the specific steps of the step (2) are as follows:
firstly, interpolating all local extreme points of an original blasting vibration signal f (t) (shown in figure 3) by using cubic splines to form an upper envelope line and a lower envelope line, wherein the upper envelope line and the lower envelope line contain all data points of the signal;
secondly, the average value of the upper envelope line and the lower envelope line obtained in the step one is recorded as m1(t) obtaining
Figure DEST_PATH_IMAGE007
Thirdly, for different blasting vibration signals f (t), if h1(t) obtained in the step (2) is an IMF component or not, if h1(t) is an IMF component, the value is recorded as c1(t) = h1 (t); if h1(t) is not an IMF component, taking the h1(t) signal obtained in the step (2) as an original signal, repeating the steps (1) to (2) k times to obtain the k-th screening data h1k (t):
Figure 319149DEST_PATH_IMAGE008
h1k (t) is not an IMF component, and there must be a termination condition for the screening process, determined by calculating the standard deviation SD between two successive processing results:
Figure DEST_PATH_IMAGE009
when the value of SD is 0.1-0.4, the first satisfied h1k (t) is the first order, and c1(t) = h1k (t), then c1(t) is the first-order IMF component of the signal f (t);
fourthly, the 1 st IMF component c resolved by the third step1(t) removing the original blasting vibration signal f (t) to obtain a residual time sequence r1(t):
Figure 195838DEST_PATH_IMAGE010
When r is1(t) is a monotonic function or exists onlyAt one extreme, empirical mode decomposition is completed, r1(t) is recorded as residual error, otherwise r obtained in the fourth step is recorded1(t) repeating steps one-four as raw data to obtain other order components of the signal, e.g. c1(t)~ ci(t) of (d). (obtaining FIG. 4)
Comparing the principal component components of the IMF components, and selecting the IMF components which can reflect the principal components of the signal information, wherein the acquisition principle is as follows: firstly, comparing IMF components with original signal characteristics, such as approximate time points of waveform change, waveform forms and the like (generally, the 1 st IMF component is a noise component, and the 2 nd to 5 th components are main characteristic components); observing the waveform information retention degree of the selected component, wherein the later component loses the original waveform information due to more screening times; and comparing whether the effect is optimal or not. (obtaining FIG. 5)
(3) Performing Hilbert transformation on the principal component of the IMF component obtained in the step (2) to obtain an envelope function of the principal vibration component; the method comprises the following specific steps:
hilbert transform is performed on principal component c (t) of the selected IMF components,
Figure DEST_PATH_IMAGE011
wherein H [ c (t) ] is Hilbert transform function of c (t), PV is Cauchy's principal value, and c (t) is selected IMF principal component;
Figure 641470DEST_PATH_IMAGE012
wherein a (t) is the envelope function of the IMF principal component c (t). (obtaining FIG. 6 or FIG. 7)
(4) Analyzing the envelope function obtained in the step (3), wherein the number of each pole on the envelope function is the number of each actual initiation section of the blasting network, comparing the number with the number of theoretical sections, and if the number is equal, generating no blind shot; if the number is smaller than the theoretical segment number, generating blind shots; (FIG. 6 shows the case of a blind gun, and FIG. 7 shows the case of a non-blind gun)
(5) And (4) obtaining the actual detonation time of each section of explosive according to the envelope function obtained in the step (4), sequentially subtracting the detonation times of the adjacent sections to obtain the actual detonation delay time of the adjacent blast holes, and comparing the actual detonation delay time with the theoretical delay time of blasting design so as to determine the position of the blind shot.
The blind shot identification mainly comprises two steps, wherein the first step is to judge whether a blind shot is generated: the method is characterized in that a single sudden change of blasting vibration signals of monitoring points on a time curve can be caused according to the fact that the detonation of each section of detonator in the differential blasting means the sudden loading of energy, the Hilbert transformation is carried out on IMF main component components on the basis, the number of each pole on a modulus diagram is the actual detonation section of a blasting network, the theoretical detonation section can be known through a blasting design network, and the theoretical detonation section are compared to judge whether blind shots are generated or not; secondly, identifying the position of the blind shot: the differential blasting delay interval can be defined as the time difference interval between the initiation moments of two adjacent sections of detonators, the actual initiation interval time between the adjacent sections can be obtained, and the actual initiation interval time is compared with the theoretical initiation interval time to judge the blind shot position.
The specific description is as follows:
FIG. 6 is a representation of IMF5 component Hilbert transform mode values (without blind shots); the time point of the abscissa corresponding to each pole on the modulus diagram is the actual detonation moment of the differential blasting detonator. As can be clearly seen from fig. 6, the Hilbert transform module value diagram of the signal principal component shows that 7 local singular points appear, and the time points of the local singular points are 55ms, 99ms, 151ms, 247ms, 343ms, 449ms and 496ms respectively. The differential blasting vibration signal shown in fig. 3 is formed by superposing 7 sections of blasting vibration waveforms. The delay interval of the differential blasting can be defined as the time difference interval between the initiation moments of two adjacent sections of detonators, and the actual initiation moments of the detonators of the obtained sections are 44ms, 52ms, 96ms, 96ms, 106 ms and 47ms respectively.
The IMF5 component Hilbert transform module values shown in FIG. 6 have 7 extreme points and no blind shot.
FIG. 7 is a principal component Hilbert transform model value (with blind shot); as is clear from fig. 6, the time points of occurrence of the 6 local singular points are 55ms, 99ms, 151ms, 247ms, 343ms and 496ms, respectively. According to the graph 5, the actual detonation moments of the detonators of the different sections are 44ms, 52ms, 96ms, 96ms and 153ms respectively.
As shown in fig. 7, the Hilbert transform module value of the principal component has 6 extreme points and generates blind shots, and it is determined that blind shots are generated and the number of the blind shots is 1; and secondly, judging the generation position of the blind shot, wherein the specific method comprises the steps of sequentially obtaining delay interval time of adjacent sections, as shown in table 1, judging that the blind shot grows at the position, namely the 6 th gun hole is not detonated, and generating the blind shot, wherein the detonation interval time of the fifth section is 153ms and is not corresponding to a designed theoretical value, but the interval time is within the sum of the theoretical interval time of the 5 th section and the theoretical interval time of the 6 th section (120 ms-430 ms). Therefore, whether the blind cannon occurs or not can be conveniently and quickly determined and the position of the blind cannon can be found through the blind cannon detection method and the blind cannon detection device.
Table 1 blasting network theory and actual delay table units: ms is
Figure 908503DEST_PATH_IMAGE014

Claims (4)

1. The blind shot identification method based on the blasting vibration signal detection is characterized by comprising the following steps of:
(1) establishing a hole-by-hole detonation type blasting network, arranging a blasting vibration tester at a position for recording blasting vibration signals, and recording original blasting vibration signals f (t);
(2) analyzing the original blasting vibration signal f (t) obtained in the step (1) by empirical mode decomposition to obtain an IMF component of the blasting vibration signal and obtain a principal component of the IMF component;
(3) performing Hilbert transformation on the principal component of the IMF component obtained in the step (2) to obtain an envelope function of the principal vibration component;
(4) analyzing the envelope function obtained in the step (3), wherein the number of each pole on the envelope function is the number of each actual initiation section of the blasting network, comparing the number with the number of theoretical sections, and if the number is equal, generating no blind shot; if the number is smaller than the theoretical segment number, generating blind shots;
(5) obtaining the actual detonation time of each section of explosive according to the envelope function obtained in the step (4), sequentially subtracting the detonation times of adjacent sections to obtain the actual detonation delay time of adjacent blast holes, and comparing the actual detonation delay time with the blasting design theory delay time to determine the position of the blind shot;
the specific steps of the step (2) are as follows:
(a) using cubic spline interpolation to form an upper envelope line and a lower envelope line of all local extreme points of an original blasting vibration signal f (t), wherein the upper envelope line and the lower envelope line comprise all data points of the signal;
(b) the average value of the upper envelope line and the lower envelope line obtained in the step (a) is recorded as m1(t) obtaining
Figure 915230DEST_PATH_IMAGE001
(c) H is obtained in the step (b) for different blasting vibration signals f (t)1(t) is an IMF component or not, if h1(t) is an IMF component, denoted as c1(t)=h1(t); if h1(t) is not an IMF component, in which case h is obtained in step (b)1(t) repeating the steps (a) to (b) k times by using the signal as an original signal to obtain the k-th screening data h1k(t):
Figure 295658DEST_PATH_IMAGE002
h1k(t) is an IMF component, determined by calculating the standard deviation SD between two consecutive processing results:
Figure 94987DEST_PATH_IMAGE003
when the SD value is 0.1-0.4, the first h is satisfied1k(t) is the first order, denoted c1(t)=h1k(t), then c1(t) is the first order IMF component of signal f (t);
(d) decomposing the first order IMF component c from step (b)1(t) removing the original blasting vibration signal f (t) to obtain a residual time sequence r1(t):
Figure 618372DEST_PATH_IMAGE004
(e) When r is1(t) is a monotonic function or only one pole exists, empirical mode decomposition is completed, and r is1(t) recording as residual, otherwise recording r obtained in step (d)1(t) repeating steps (a) - (d) as raw data to obtain other order components c of the signal1(t)~ ci(t)。
2. The blind shot identification method based on the blast vibration signal detection according to claim 1, characterized in that: the blasting vibration tester is a TC-4850 blasting vibration tester.
3. The blind shot identification method based on the blast vibration signal detection according to claim 1, characterized in that: the SD value is 0.2-0.3.
4. The blind shot identification method based on the blast vibration signal detection according to claim 1, characterized in that: c. C2(t) or c3(t) or c4(t) is a principal component.
CN201910022593.4A 2019-01-10 2019-01-10 Blind gun identification method based on blasting vibration signal detection Active CN109708550B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910022593.4A CN109708550B (en) 2019-01-10 2019-01-10 Blind gun identification method based on blasting vibration signal detection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910022593.4A CN109708550B (en) 2019-01-10 2019-01-10 Blind gun identification method based on blasting vibration signal detection

Publications (2)

Publication Number Publication Date
CN109708550A CN109708550A (en) 2019-05-03
CN109708550B true CN109708550B (en) 2021-06-18

Family

ID=66259942

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910022593.4A Active CN109708550B (en) 2019-01-10 2019-01-10 Blind gun identification method based on blasting vibration signal detection

Country Status (1)

Country Link
CN (1) CN109708550B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110865171B (en) * 2019-11-14 2022-04-19 北京龙德时代技术服务有限公司 Blasting safety analysis method and system based on digital noise detection
CN112611275B (en) * 2020-12-14 2022-10-18 江西理工大学 Detection method for blasting blind gun

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102997769A (en) * 2012-10-31 2013-03-27 江西理工大学 Method for detecting and identifying blind shot for blasting chambers
CN103256874A (en) * 2013-05-17 2013-08-21 四川雅化实业集团工程爆破有限公司 Misfire detecting method by means of blast vibration waves
CN103344156A (en) * 2013-07-16 2013-10-09 四川大学 Blind cannonball identification method in blasting works
CN104568024A (en) * 2015-01-21 2015-04-29 山东理工大学 Vibration type flow meter characteristic signal extraction method
CN204944329U (en) * 2015-09-17 2016-01-06 平顶山市碧源科技有限公司 A kind of colliery miso-fire search instrument

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102997769A (en) * 2012-10-31 2013-03-27 江西理工大学 Method for detecting and identifying blind shot for blasting chambers
CN103256874A (en) * 2013-05-17 2013-08-21 四川雅化实业集团工程爆破有限公司 Misfire detecting method by means of blast vibration waves
CN103344156A (en) * 2013-07-16 2013-10-09 四川大学 Blind cannonball identification method in blasting works
CN104568024A (en) * 2015-01-21 2015-04-29 山东理工大学 Vibration type flow meter characteristic signal extraction method
CN204944329U (en) * 2015-09-17 2016-01-06 平顶山市碧源科技有限公司 A kind of colliery miso-fire search instrument

Also Published As

Publication number Publication date
CN109708550A (en) 2019-05-03

Similar Documents

Publication Publication Date Title
CN109708550B (en) Blind gun identification method based on blasting vibration signal detection
CN110764063B (en) Radar signal sorting method based on combination of SDIF and PRI transformation method
CN111337888B (en) Dense decoy jamming method, computer device and computer readable storage medium
CN108132914B (en) High-confidence statistical inference and confirmation method for failure rate of limited number of product performances
CN110207974B (en) Circuit breaker fault identification method based on vibration signal time-frequency energy distribution characteristics
CN105117550A (en) Product multidimensional correlation-oriented degradation failure modeling method
CN110956004B (en) Concrete penetration depth prediction algorithm with fusion of multi-empirical formula and BP neural network model
CN107024141A (en) Sound and vibration monitoring and the defect positioning method of quality are assembled for Ramming Device
CN108226868A (en) A kind of localization method and system of shell water-drop point
CN108897354A (en) A kind of aluminium fusion process fire box temperature prediction technique based on depth confidence network
CN115166563A (en) Power battery aging state evaluation and decommissioning screening method and system
CN106324583A (en) Vector array passive sonar abnormal data online eliminating method
CN111862036B (en) Tool abnormal placement detection method and system based on double-index measurement learning
CN103487828B (en) Method for automatically detecting earthquake auxiliary channel
CN109669211B (en) Method and device for predicting seabed lithologic thin sand body reservoir, storage medium and terminal
CN111562111A (en) Engine cold state test fault diagnosis method
CN108509576B (en) Electronic map-based filtering optimization method for wrong point positions of patrol track
CN111123356B (en) Abnormal track intelligent identification method based on first arrival information
CN112541392A (en) Open bench blasting prediction method based on deep neural network
CN108960032B (en) Tristable logic stochastic resonance method
Stepanova et al. The development of a defect-rejection procedure for multiple-pass welding by the distribution of the principal parameters of acoustic-emission signals
CN111981905A (en) Method for resisting high repetition frequency laser active interference
CN117540326B (en) Construction state abnormality identification method and system for tunnel construction equipment by drilling and blasting method
CN110352334A (en) Hit detection method, strike detection device and armoring trolley
Song et al. A Radar Signal Sorting Algorithm with Adaptive PRI Box Overlap Rate

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