Method for determining detonation time sequence of blast hole and related equipment

Abstract

The invention provides a blasthole detonation time sequence determination method and related equipment, and relates to the field of blasting vibration control, and the method comprises the following steps: obtaining vibration waveforms of a plurality of blastholes; determining a recombination intrinsic mode component IMF corresponding to the vibration waveform of each blast hole, and extracting a more accurate vibration velocity peak value of each blast hole from the recombination IMF corresponding to the vibration waveform of each blast hole; determining a target IMF corresponding to the vibration waveform of each blast hole based on the vibration velocity peak value of the reference blast hole, the recombined IMF corresponding to the vibration waveform and the vibration velocity peak values of other blast holes; and determining a plurality of superposition detonation time sequences, combining each superposition detonation time sequence with the target IMF corresponding to the vibration waveform of each blast hole to obtain a plurality of synthetic IMFs, extracting a plurality of synthetic vibration velocity peak values from the synthetic IMFs, and taking the time sequence corresponding to the minimum synthetic vibration velocity peak value as the detonation time sequence of the blast hole. According to the scheme, coordinated blasting among the blast holes can be achieved, the detonating time sequence is optimized, and the blasting effect can be improved.

Description

Technical Field

[0001] The present application belongs to the field of blasting vibration control, and in particular, relates to a method for determining a detonation time sequence of a blasthole and related equipment. Background Art

[0002] With the rapid development of blasting technology and the continuous expansion of blasting scale, in large-scale deep hole blasting, there are a large number of blastholes and their distribution range is wide. The detonation time sequence of blastholes is not only a key means to achieve blasting vibration control and improve blasting efficiency and safety, but also an important driving force to promote the development of blasting technology to a higher level. Therefore, it is particularly important to determine the detonation time sequence of blastholes.

[0003] At present, the method based on empirical formula or experimental summary is usually used to determine the detonation time sequence of blastholes.

[0004] However, empirical formulas and experimental summary methods often cannot fully consider the various complex factors in the blasting process, and the calculation results may have large errors, resulting in poor blasting effects. Summary of the invention

[0005] In view of the above problems, this application proposes a method for determining the detonation time sequence of a blasthole and related equipment. In order to improve the blasting effect, the specific scheme is as follows:

[0006] The first aspect of the present application provides a method for determining a blast hole detonation time sequence, comprising:

[0007] Obtain vibration waveforms of multiple blastholes;

[0008] Based on the vibration waveform of each of the blastholes, determining the reorganized IMF corresponding to the vibration waveform of each of the blastholes, and extracting the vibration velocity peak value of each of the blastholes from the reorganized IMF corresponding to the vibration waveform of each of the blastholes;

[0009] Determine a reference blasthole from each of the blastholes, and determine a target IMF corresponding to the vibration waveform of each of the blastholes based on the vibration velocity peak value and the recombined IMF corresponding to the vibration waveform of the reference blasthole and the vibration velocity peak value of other blastholes, wherein the other blastholes are blastholes other than the reference blasthole among the blastholes;

[0010] Determine a plurality of superimposed detonation time sequences, and combine each of the superimposed detonation time sequences with a target IMF corresponding to the vibration waveform of each of the blastholes according to the Anderson principle to obtain a plurality of synthetic IMFs;

[0011] A plurality of synthetic vibration velocity peak values ​​are extracted from each of the synthetic IMFs, and a time series corresponding to a minimum synthetic vibration velocity peak value is used as a detonation time series of the plurality of blastholes.

[0012] In a possible implementation, determining the reorganized IMF corresponding to the vibration waveform of each blasthole based on the vibration waveform of each blasthole includes:

[0013] Decomposing the vibration waveform of each blasthole to obtain the intrinsic modal component IMF corresponding to the vibration waveform of each blasthole;

[0014] The number of layers of the IMF corresponding to the vibration waveform of each blasthole is obtained, and the vibration waveform of each blasthole is re-decomposed based on the IMF with the least number of layers to obtain the recombined IMF corresponding to the vibration waveform of each blasthole.

[0015] In a possible implementation, the reconstructed IMF corresponding to the vibration waveform of each blasthole based on the vibration velocity peak value and the vibration waveform of the reference blasthole and the vibration velocity peak value of other blastholes is determined, including:

[0016] Based on the vibration velocity peak value of the reference blasthole and the vibration velocity peak value of the other blastholes, the activation coefficients of the other blastholes are calculated;

[0017] Based on the reorganized IMF corresponding to the vibration waveform of the reference blasthole and the activation coefficients of the other blastholes, the reorganized IMF corresponding to the vibration waveform of the other blastholes is calculated;

[0018] The recombined IMF corresponding to the vibration waveform of each blasthole is linearly superimposed and expanded in time domain to obtain the target IMF corresponding to the vibration waveform of each blasthole.

[0019] In a possible implementation, determining a plurality of superimposed detonation time sequences includes:

[0020] Based on the acquired propagation speed of the vibration waveform of each blast hole and the target distance of each blast hole, a detonation time difference sequence is calculated, wherein the target distance is the distance between the position of the blast hole and the blasting vibration monitoring position;

[0021] Based on a plurality of preset detonation time sequences and the detonation time difference sequence, a plurality of superimposed detonation time sequences are calculated.

[0022] A second aspect of the present application provides a device for determining a blast hole detonation time sequence, comprising:

[0023] An acquisition unit, used for acquiring vibration waveforms of a plurality of blastholes;

[0024] A first determination unit is used to determine the reorganized IMF corresponding to the vibration waveform of each blasthole based on the vibration waveform of each blasthole, and extract the vibration velocity peak value of each blasthole from the reorganized IMF corresponding to the vibration waveform of each blasthole;

[0025] A second determination unit is used to determine a reference blasthole from each of the blastholes, and determine a target IMF corresponding to the vibration waveform of each of the blastholes based on the vibration velocity peak value and the recombined IMF corresponding to the vibration waveform of the reference blasthole and the vibration velocity peak value of other blastholes, wherein the other blastholes are blastholes other than the reference blasthole among the blastholes;

[0026] A third determination unit is used to determine a plurality of superimposed detonation time sequences, and combine each of the superimposed detonation time sequences with a target IMF corresponding to the vibration waveform of each of the blastholes according to the Anderson principle to obtain a plurality of synthetic IMFs;

[0027] The fourth determination unit is used to extract multiple synthetic vibration velocity peak values ​​from each of the synthetic IMFs, and use the time series corresponding to the minimum synthetic vibration velocity peak value as the detonation time series of the multiple blastholes.

[0028] In a possible implementation, the first determining unit for determining the reorganized IMF corresponding to the vibration waveform of each blasthole based on the vibration waveform of each blasthole includes:

[0029] A decomposition subunit, used for decomposing the vibration waveform of each blast hole to obtain an intrinsic modal component IMF corresponding to the vibration waveform of each blast hole;

[0030] The first determination subunit is used to obtain the number of IMF layers corresponding to the vibration waveform of each blasthole, and re-decompose the vibration waveform of each blasthole based on the IMF with the least number of layers to obtain the recombined IMF corresponding to the vibration waveform of each blasthole.

[0031] In a possible implementation, the second determining unit for determining the target IMF corresponding to the vibration waveform of each blasthole based on the vibration velocity peak value and the reconstructed IMF corresponding to the vibration waveform of the reference blasthole and the vibration velocity peak value of other blastholes includes:

[0032] A first calculation subunit is used to calculate the activation coefficient of the other blast holes based on the vibration velocity peak value of the reference blast hole and the vibration velocity peak value of the other blast holes;

[0033] A second calculation subunit is used to calculate the reorganized IMF corresponding to the vibration waveform of the other blastholes based on the reorganized IMF corresponding to the vibration waveform of the reference blasthole and the activation coefficient of the other blastholes;

[0034] The expansion subunit is used to linearly superpose and time-domain expand the recombined IMF corresponding to the vibration waveform of each blast hole to obtain the target IMF corresponding to the vibration waveform of each blast hole.

[0035] The third aspect of the present application provides a computer program product, including computer-readable instructions. When the computer-readable instructions are executed on an electronic device, the electronic device implements the method for determining the detonation time sequence of a blasthole according to the first aspect or any implementation of the first aspect.

[0036] A fourth aspect of the present application provides an electronic device, comprising at least one processor and a memory connected to the processor, wherein:

[0037] The memory is used to store computer programs;

[0038] The processor is used to execute the computer program so that the electronic device can implement the method for determining the detonation time sequence of the blasthole according to the first aspect or any implementation manner of the first aspect.

[0039] A fifth aspect of the present application provides a computer storage medium, which carries one or more computer programs. When the one or more computer programs are executed by an electronic device, the electronic device can implement the method for determining the detonation time sequence of a blast hole according to the first aspect or any implementation of the first aspect.

[0040] Based on the above technical scheme, the present application provides a method for determining the detonation time series of a blasthole and related equipment, the method obtains the positions of multiple blastholes in the blasting area and the blasting vibration monitoring positions, and collects the vibration waveform of each blasthole; based on the vibration waveform of each blasthole, determines the recombined IMF corresponding to the vibration waveform of each blasthole, and extracts a more accurate vibration velocity peak value of each blasthole from the recombined IMF corresponding to the vibration waveform of each blasthole; determines a reference blasthole from each blasthole, and determines the target IMF corresponding to the vibration waveform of each blasthole based on the vibration velocity peak value and the recombined IMF corresponding to the vibration waveform of the reference blasthole, as well as the vibration velocity peak values ​​of other blastholes, the other blastholes being the blastholes except the reference blasthole in each blasthole; determines multiple superimposed detonation time series, and combines each superimposed detonation time series with the target IMF corresponding to the vibration waveform of each blasthole according to the Anderson principle to obtain multiple synthetic IMFs; extracts multiple synthetic vibration velocity peaks from each synthetic IMF, and uses the time series corresponding to the smallest synthetic vibration velocity peak as the detonation time series of multiple blastholes. The scheme optimizes the blast hole detonation time sequence, which can achieve coordinated blasting among different blast holes and is beneficial to improving the blasting effect. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] Various other advantages and benefits will become apparent to those of ordinary skill in the art by reading the detailed description of the preferred embodiments below. The accompanying drawings are only for the purpose of illustrating the preferred embodiments and are not to be considered as limiting the present application. Also, the same reference symbols are used throughout the accompanying drawings to represent the same components. In the accompanying drawings:

[0042] Figure 1 A schematic flow chart of a method for detonating a blast hole in an embodiment of the present application;

[0043] Figure 2 It is a structural schematic diagram of a blast hole detonation time sequence device disclosed in an embodiment of the present application;

[0044] Figure 3 The hardware structure block diagram of a blast hole detonation time sequence device disclosed in an embodiment of the present application. DETAILED DESCRIPTION

[0045] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0046] The terms "first", "second", etc. in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and need not be used to describe a specific order or sequential order. It should be understood that the terms used in this way can be interchangeable under appropriate circumstances, which is only to describe the distinction mode adopted by the objects of the same attributes when describing in the embodiments of the present application. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, so that the process, method, system, product or equipment comprising a series of units need not be limited to those units, but may include other units that are not clearly listed or inherent to these processes, methods, products or equipment.

[0047] In order to improve the blasting effect, the present application provides a method for determining a regression model of coal coke thermal properties. The following is a further detailed description of the method for reflecting implicit orders in the trading quotation provided by the present application in combination with the accompanying drawings and specific implementation methods.

[0048] Please refer to the attached Figure 1 , Figure 1 A flow chart of a method for determining a blast hole detonation time sequence provided in an embodiment of the present application. The method may include the following steps:

[0049] Step S101, obtaining vibration waveforms of multiple blast holes.

[0050] In the present application, a sensor (such as a vibration sensor or an accelerometer) may be used to record the vibration waveform generated by each blast hole during the blasting process, and the vibration waveform includes the vibration signal transmitted from each blast hole.

[0051] Step S102: based on the vibration waveform of each blast hole, determine the recombined natural modal component IMF corresponding to the vibration waveform of each blast hole, and extract the vibration velocity peak value of each blast hole from the recombined IMF corresponding to the vibration waveform of each blast hole.

[0052] It should be noted that IMF (Intrinsic Mode Function) is a basic concept in the empirical mode decomposition method, which represents the local characteristic scale in the signal. Each IMF satisfies two conditions: in the entire data set, the number of extreme points and the number of zero crossings must be equal or differ by at most one; at any point, the average value of the upper envelope formed by the local maximum points and the lower envelope formed by the local minimum points is zero.

[0053] In the present application, based on the vibration waveform of each blasthole, the reorganized IMF corresponding to the vibration waveform of each blasthole is determined, including: decomposing the vibration waveform of each blasthole to obtain the IMF corresponding to the vibration waveform of each blasthole. Obtaining the number of layers of the IMF corresponding to the vibration waveform of each blasthole, re-decomposing the vibration waveform of each blasthole based on the IMF with the least number of layers, and obtaining the reorganized IMF corresponding to the vibration waveform of each blasthole.

[0054] The complex vibration waveform can be decomposed into a series of relatively simple IMFs through empirical mode decomposition or its improved methods (such as ensemble empirical mode decomposition), which represent the characteristics of the signal in different frequency bands. Since the vibration waveforms of different boreholes may have different complexities and frequency characteristics, the number of IMF layers obtained after EMD decomposition may also be different. The IMF with the least number of layers can be selected as the benchmark to re-decompose the vibration waveforms of all boreholes. Then, the vibration velocity peaks are extracted from the reorganized IMF of each borehole, which reflect the maximum vibration intensity generated when the borehole explodes.

[0055] By obtaining the charge amount and target distance of the blasthole, we can establish a functional relationship between the peak value of the blasthole's vibration velocity, the charge amount, and the target distance. The target distance is used to indicate the distance between the position of the blasthole and the position of the vibration monitoring point. For details, please refer to the following formula:

[0056]

[0057] Where: V k is the peak value of the vibration velocity of the kth reorganized IMF, k∈{1, 2, ..., K min}, K min is the minimum number of layers, Q is the charge of the blasthole, R is the distance between the blasthole and the vibration monitoring point, α k , β k It is a parameter related to the performance of explosives, rock properties, and terrain and geological conditions.

[0058] Step S103, determine a reference borehole from each borehole, and determine the target IMF corresponding to the vibration waveform of each borehole based on the vibration velocity peak value and the reconstructed IMF corresponding to the vibration waveform of the reference borehole and the vibration velocity peak value of other boreholes. The other boreholes are the boreholes except the reference borehole among the boreholes.

[0059] In the present application, a key position for vibration control is selected from the vibration monitoring points, and the blast hole corresponding to the position can be used as a reference blast hole, and the activation coefficients of other blast holes are calculated based on the vibration velocity peak value of the reference blast hole and the vibration velocity peak values ​​of other blast holes. Based on the reorganized IMF corresponding to the vibration waveform of the reference blast hole and the activation coefficients of other blast holes, the reorganized IMF corresponding to the vibration waveform of other blast holes is calculated. The reorganized IMF corresponding to the vibration waveform of each blast hole is linearly superimposed and time-domain expanded to obtain the target IMF corresponding to the vibration waveform of each blast hole.

[0060] Specifically, based on the assumption that the functional relationship between the peak value of the vibration velocity of the blasthole, the charge amount, and the target distance is valid at any same time within the duration range of the vibration waveform, the ratio of the peak value of the vibration velocity of other blastholes to the peak value of the vibration velocity of the benchmark blasthole is calculated and used as the activation coefficient. The activation coefficient can be regarded as a scaling or adjustment factor of the vibration waveform of other blastholes relative to the vibration waveform of the benchmark blasthole. For details, please refer to the following formula:

[0061]

[0062] Where: A j,k (Q0,R0,Q j ,R j ) is the activation coefficient of the k-th layer IMF corresponding to the vibration waveform of the j-th blasthole, v j,k is the vibration velocity peak value of other blastholes; v0 is the vibration velocity peak value of the reference blasthole, Q0 is the charge of the reference blasthole, R0 is the distance from the reference blasthole to the key position of vibration control, Q j is the charge of the jth other blasthole, R j is the distance from the jth other blasthole to the key position of vibration control.

[0063] Based on the reorganized IMF corresponding to the vibration waveform of the reference blasthole, multiply it by the activation coefficient of other blastholes to obtain the reorganized IMF corresponding to the vibration waveform of other blastholes. For details, please refer to the following formula:

[0064]

[0065] Where: f' j,k (t) is the k-th layer IMF of the j-th other blasthole, f 0,k (t) is the kth layer IMF of the benchmark blasthole.

[0066] Combine the IMFs of the jth other blastholes to get the IMF of the jth other blasthole. For details, please refer to the following formula:

[0067]

[0068] Where: f' j (t) is the IMF of the jth other blasthole.

[0069] The reorganized IMF corresponding to the vibration waveform of each blasthole is linearly superimposed and expanded in the time domain to obtain the target IMF corresponding to the vibration waveform of each blasthole. For details, please refer to the following formula:

[0070]

[0071] Where: f j (t) is the target IMF of the jth other blasthole in the full time domain, and tdur is the duration of the reference blasthole f0(t).

[0072] Step S104: determine a plurality of superimposed detonation time sequences, and combine each superimposed detonation time sequence with the target IMF corresponding to the vibration waveform of each blasthole according to the Anderson principle to obtain a plurality of synthetic IMFs.

[0073] In the present application, the detonation time difference sequence can be calculated based on the distance from the reference blasthole to the key position of vibration control and the distance from the jth other blasthole to the key position of vibration control, as well as the propagation speed of the vibration waveform of each blasthole obtained, and multiple superimposed detonation time sequences can be calculated based on multiple preset detonation time sequences and detonation time difference sequences. For details, please refer to the following formula:

[0074] t j-0 =(R j -R0) / c v

[0075] t′ j =t j +t j-0

[0076] Where: t j-0 is the time difference between the vibration wave generated by the jth other blasthole and the vibration wave generated by the reference blasthole and transmitted to the key position of vibration control, c v is the longitudinal wave velocity, t' j is the time when the jth other blasthole actually participates in the superposition calculation, t j is the detonation time of the jth other blast hole.

[0077] According to the Anderson principle, each superimposed detonation time series is combined with the target IMF corresponding to the vibration waveform of each blast hole to obtain multiple synthetic IMFs. For details, please refer to the following formula:

[0078]

[0079] In the formula, f j (t-t'j) is the target IMF corresponding to the vibration waveform of the jth other blasthole; F(t) is the synthetic IMF.

[0080] Furthermore, the detonation time sequence is changed to obtain the synthetic IMF corresponding to each detonation time sequence in turn.

[0081] Step S105: extract multiple synthetic vibration velocity peak values ​​from each synthetic IMF, and use the time series corresponding to the minimum synthetic vibration velocity peak value as the detonation time series of multiple blastholes.

[0082] In this application, multiple synthetic vibration velocity peak values ​​Vmax(t,{t j +t j-0}), the composite vibration velocity peak is the maximum positive vibration velocity and the maximum negative vibration velocity in the composite waveform. Then from multiple composite vibration velocity peaks Vmax(t,{t j +t j-0}) to determine the minimum peak value of the synthetic vibration velocity V min For details, please refer to the following formula:

[0083] V min = {V max (t,{t j +t j-0})} min

[0084] Finally, the time series corresponding to the minimum synthetic vibration velocity peak is used as the detonation time series of multiple blastholes.

[0085] In summary, the present application provides a method for determining the detonation time series of a blasthole, which obtains the positions of multiple blastholes in the blasting area and the blasting vibration monitoring positions, and collects the vibration waveform of each blasthole. Based on the vibration waveform of each blasthole, the recombined IMF corresponding to the vibration waveform of each blasthole is determined, and a more accurate vibration velocity peak value of each blasthole is extracted from the recombined IMF corresponding to the vibration waveform of each blasthole. A reference blasthole is determined from each blasthole, and the target IMF corresponding to the vibration waveform of each blasthole is determined based on the vibration velocity peak value of the reference blasthole and the recombined IMF corresponding to the vibration waveform, as well as the vibration velocity peak values ​​of other blastholes. The other blastholes are the blastholes except the reference blasthole in each blasthole. Multiple superimposed detonation time series are determined, and each superimposed detonation time series is combined with the target IMF corresponding to the vibration waveform of each blasthole according to the Anderson principle to obtain multiple synthetic IMFs. Multiple synthetic vibration velocity peaks are extracted from each synthetic IMF, and the time series corresponding to the smallest synthetic vibration velocity peak is used as the detonation time series of multiple blastholes. The scheme optimizes the blast hole detonation time sequence, which can achieve coordinated blasting among different blast holes and is beneficial to improving the blasting effect.

[0086] The method is described in detail in the embodiments disclosed in the above-mentioned application. The method of the application can be implemented by various forms of devices. Therefore, the application also discloses a device for determining the detonation time sequence of a blasthole. A specific embodiment is given below for detailed description.

[0087] Please refer to the attached Figure 2 , Figure 3 This is a schematic diagram of the structure of a device for determining the detonation time sequence of a blasthole disclosed in an embodiment of the present application, the device comprising:

[0088] The acquisition unit 11 is used to acquire the vibration waveforms of multiple blast holes.

[0089] The first determination unit 12 is used to determine the reorganized IMF corresponding to the vibration waveform of each blasthole based on the vibration waveform of each blasthole, and extract the vibration velocity peak value of each blasthole from the reorganized IMF corresponding to the vibration waveform of each blasthole.

[0090] The second determination unit 13 is used to determine a reference borehole from each of the boreholes, and determine a target IMF corresponding to the vibration waveform of each of the boreholes based on the vibration velocity peak value and the reconstructed IMF corresponding to the vibration waveform of the reference borehole and the vibration velocity peak values ​​of other boreholes. The other boreholes are the boreholes among the boreholes except the reference borehole.

[0091] The third determination unit 14 is used to determine a plurality of superimposed detonation time sequences, and combine each of the superimposed detonation time sequences with a target IMF corresponding to the vibration waveform of each of the blastholes according to the Anderson principle to obtain a plurality of synthetic IMFs.

[0092] The fourth determining unit 15 is used to extract multiple synthetic vibration velocity peak values ​​from each of the synthetic IMFs, and use the time series corresponding to the minimum synthetic vibration velocity peak value as the detonation time series of the multiple blastholes.

[0093] As an implementable embodiment, the first determining unit 12 for determining the reorganized IMF corresponding to the vibration waveform of each blasthole based on the vibration waveform of each blasthole includes:

[0094] The decomposition subunit is used to decompose the vibration waveform of each blast hole to obtain the intrinsic modal component IMF corresponding to the vibration waveform of each blast hole.

[0095] The first determination subunit is used to obtain the number of IMF layers corresponding to the vibration waveform of each blasthole, and re-decompose the vibration waveform of each blasthole based on the IMF with the least number of layers to obtain the recombined IMF corresponding to the vibration waveform of each blasthole.

[0096] As an implementable embodiment, the second determining unit 13 for determining the target IMF corresponding to the vibration waveform of each blasthole based on the vibration velocity peak value and the recombined IMF corresponding to the vibration waveform of the reference blasthole and the vibration velocity peak value of other blastholes includes:

[0097] The first calculation subunit is used to calculate the activation coefficients of the other blast holes based on the vibration velocity peak value of the reference blast hole and the vibration velocity peak values ​​of the other blast holes.

[0098] The second calculation subunit is used to calculate the reorganized IMF corresponding to the vibration waveform of the other blastholes based on the reorganized IMF corresponding to the vibration waveform of the reference blasthole and the activation coefficients of the other blastholes.

[0099] The expansion subunit is used to linearly superpose and time-domain expand the recombined IMF corresponding to the vibration waveform of each blast hole to obtain the target IMF corresponding to the vibration waveform of each blast hole.

[0100] As an implementable embodiment, the third determination unit 14 for determining a plurality of superimposed detonation time sequences includes:

[0101] The third calculation subunit is used to calculate the detonation time difference sequence based on the acquired propagation speed of the vibration waveform of each blast hole and the target distance of each blast hole, and the target distance is the distance between the position of the blast hole and the blasting vibration monitoring position.

[0102] The fourth calculation subunit is used to calculate a plurality of superimposed detonation time sequences based on a plurality of preset detonation time sequences and the detonation time difference sequence.

[0103] The present application also provides an electronic device in an embodiment. Figure 3 As shown, it shows a schematic diagram of the structure of an electronic device suitable for implementing the embodiments of the present application. The electronic device in the embodiments of the present application may include but is not limited to fixed terminals such as mobile phones, notebook computers, PDAs (personal digital assistants), PADs (tablet computers), desktop computers, etc. Figure 3 The electronic device shown is merely an example and should not bring any limitation to the functions and scope of use of the embodiments of the present application.

[0104] like Figure 3 As shown, the electronic device may include a processing device (e.g., a central processing unit, a graphics processing unit, etc.) 301, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 302 or a program loaded from a storage device 308 to a random access memory (RAM) 303. When the electronic device is powered on, various programs and data required for the operation of the electronic device are also stored in the RAM 303. The processing device 301, the ROM 302, and the RAM 303 are connected to each other via a bus 304. An input / output (I / O) interface 303 is also connected to the bus 304.

[0105] Typically, the following devices may be connected to the I / O interface 305: input devices 306 including, for example, a touch screen, a touch pad, a keyboard, a mouse, a camera, a microphone, an accelerometer, a gyroscope, etc.; output devices 307 including, for example, a liquid crystal display (LCD), a speaker, a vibrator, etc.; storage devices 308 including, for example, a memory card, a hard disk, etc.; and communication devices 309. The communication device 309 may allow the electronic device to communicate with other devices wirelessly or by wire to exchange data. Although Figure 3 An electronic device having various devices is shown, but it should be understood that it is not required to implement or possess all the devices shown. More or fewer devices may be implemented or possessed instead.

[0106] Also provided in an embodiment of the present application is a computer program product including computer-readable instructions. When the computer-readable instructions are executed on an electronic device, the electronic device implements any method for determining the detonation time sequence of a blast hole provided in an embodiment of the present application.

[0107] A computer-readable storage medium is also provided in an embodiment of the present application. The storage medium carries one or more computer programs. When the one or more computer programs are executed by an electronic device, the electronic device can implement any method for determining the detonation time sequence of a blast hole provided in an embodiment of the present application.

[0108] It should also be noted that the device embodiments described above are merely schematic, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed over multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. In addition, in the drawings of the device embodiments provided by the present application, the connection relationship between the modules indicates that there is a communication connection between them, which may be specifically implemented as one or more communication buses or signal lines.

[0109] Through the description of the above implementation mode, the technicians in the field can clearly understand that the present application can be implemented by means of software plus necessary general hardware, and of course, it can also be implemented by special hardware including special integrated circuits, special CPUs, special memories, special components, etc. In general, all functions completed by computer programs can be easily implemented by corresponding hardware, and the specific hardware structure used to implement the same function can also be various, such as analog circuits, digital circuits or special circuits. However, for the present application, software program implementation is a better implementation mode in more cases. Based on such an understanding, the technical solution of the present application is essentially or the part that contributes to the prior art can be embodied in the form of a software product, which is stored in a readable storage medium, such as a computer floppy disk, a U disk, a mobile hard disk, a ROM, a RAM, a disk or an optical disk, etc., including a number of instructions to enable a computer device (which can be a personal computer, a training device, or a network device, etc.) to execute the methods described in each embodiment of the present application.

[0110] In the above embodiments, all or part of the embodiments may be implemented by software, hardware, firmware or any combination thereof. When implemented by software, all or part of the embodiments may be implemented in the form of a computer program product.

[0111] The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, the process or function described in the embodiment of the present application is generated in whole or in part. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices. The computer instructions may be stored in a computer-readable storage medium, or transmitted from one computer-readable storage medium to another computer-readable storage medium. For example, the computer instructions may be transmitted from a website site, a computer, a training device, or a data center by wired (e.g., coaxial cable, optical fiber, digital subscriber line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) mode to another website site, computer, training device, or data center. The computer-readable storage medium may be any available medium that a computer can store or a data storage device such as a training device, a data center, etc. that includes one or more available media integrations. The available medium may be a magnetic medium, (e.g., a floppy disk, a hard disk, a tape), an optical medium (e.g., a DVD), or a semiconductor medium (e.g., a solid-state drive (SSD)), etc.

Claims

1. A method for determining the detonation time sequence of a blasthole, characterized in that: include: Obtain vibration waveforms of multiple blastholes; Based on the vibration waveform of each of the blastholes, a recombined intrinsic modal component IMF corresponding to the vibration waveform of each of the blastholes is determined, and a vibration velocity peak value of each of the blastholes is extracted from the recombined IMF corresponding to the vibration waveform of each of the blastholes; Determine a reference blasthole from each of the blastholes, and determine a target IMF corresponding to the vibration waveform of each of the blastholes based on the vibration velocity peak value and the recombined IMF corresponding to the vibration waveform of the reference blasthole and the vibration velocity peak value of other blastholes, wherein the other blastholes are blastholes other than the reference blasthole among the blastholes; Determine a plurality of superimposed detonation time sequences, and combine each of the superimposed detonation time sequences with a target IMF corresponding to the vibration waveform of each of the blastholes according to the Anderson principle to obtain a plurality of synthetic IMFs; A plurality of synthetic vibration velocity peak values ​​are extracted from each of the synthetic IMFs, and a time series corresponding to a minimum synthetic vibration velocity peak value is used as a detonation time series of the plurality of blastholes.

2. The method for determining the detonation time sequence of a blasthole according to claim 1, characterized in that: The step of determining the reorganized IMF corresponding to the vibration waveform of each blasthole based on the vibration waveform of each blasthole comprises: Decomposing the vibration waveform of each blasthole to obtain the IMF corresponding to the vibration waveform of each blasthole; The number of layers of the IMF corresponding to the vibration waveform of each blasthole is obtained, and the vibration waveform of each blasthole is re-decomposed based on the IMF with the least number of layers to obtain the recombined IMF corresponding to the vibration waveform of each blasthole.

3. The method for determining the detonation time sequence of a blasthole according to claim 1, characterized in that: The method of determining the target IMF corresponding to the vibration waveform of each blasthole based on the vibration velocity peak value and the recombined IMF corresponding to the vibration waveform of the reference blasthole and the vibration velocity peak values ​​of other blastholes includes: Based on the vibration velocity peak value of the reference blasthole and the vibration velocity peak value of the other blastholes, the activation coefficients of the other blastholes are calculated; Based on the reorganized IMF corresponding to the vibration waveform of the reference blasthole and the activation coefficients of the other blastholes, the reorganized IMF corresponding to the vibration waveform of the other blastholes is calculated; The recombined IMF corresponding to the vibration waveform of each blasthole is linearly superimposed and expanded in time domain to obtain the target IMF corresponding to the vibration waveform of each blasthole.

4. The method for determining the detonation time sequence of a blasthole according to claim 1, characterized in that: The determining of a plurality of superimposed detonation time sequences comprises: Based on the acquired propagation speed of the vibration waveform of each blast hole and the target distance of each blast hole, a detonation time difference sequence is calculated, wherein the target distance is the distance between the position of the blast hole and the blasting vibration monitoring position; Based on a plurality of preset detonation time sequences and the detonation time difference sequence, a plurality of superimposed detonation time sequences are calculated.

5. A device for determining the detonation time sequence of a blasthole, characterized in that: include: An acquisition unit, used for acquiring vibration waveforms of a plurality of blastholes; A first determination unit is used to determine the reorganized IMF corresponding to the vibration waveform of each blasthole based on the vibration waveform of each blasthole, and extract the vibration velocity peak value of each blasthole from the reorganized IMF corresponding to the vibration waveform of each blasthole; A second determination unit is used to determine a reference blasthole from each of the blastholes, and determine a target IMF corresponding to the vibration waveform of each of the blastholes based on the vibration velocity peak value and the recombined IMF corresponding to the vibration waveform of the reference blasthole and the vibration velocity peak value of other blastholes, wherein the other blastholes are blastholes other than the reference blasthole among the blastholes; A third determination unit is used to determine a plurality of superimposed detonation time sequences, and combine each of the superimposed detonation time sequences with a target IMF corresponding to the vibration waveform of each of the blastholes according to the Anderson principle to obtain a plurality of synthetic IMFs; The fourth determination unit is used to extract multiple synthetic vibration velocity peak values ​​from each of the synthetic IMFs, and use the time series corresponding to the minimum synthetic vibration velocity peak value as the detonation time series of the multiple blastholes.

6. The device for determining the detonation time sequence of a blasthole according to claim 5, characterized in that: The first determining unit for determining the reorganized IMF corresponding to the vibration waveform of each blasthole based on the vibration waveform of each blasthole includes: A decomposition subunit, used for decomposing the vibration waveform of each blast hole to obtain an intrinsic modal component IMF corresponding to the vibration waveform of each blast hole; The first determination subunit is used to obtain the number of IMF layers corresponding to the vibration waveform of each blasthole, and re-decompose the vibration waveform of each blasthole based on the IMF with the least number of layers to obtain the recombined IMF corresponding to the vibration waveform of each blasthole.

7. The device for determining the detonation time sequence of a blasthole according to claim 5, characterized in that: The second determination unit for determining the target IMF corresponding to the vibration waveform of each blasthole based on the vibration velocity peak value and the recombined IMF corresponding to the vibration waveform of the reference blasthole and the vibration velocity peak values ​​of other blastholes includes: A first calculation subunit is used to calculate the activation coefficient of the other blast holes based on the vibration velocity peak value of the reference blast hole and the vibration velocity peak value of the other blast holes; A second calculation subunit is used to calculate the reorganized IMF corresponding to the vibration waveform of the other blastholes based on the reorganized IMF corresponding to the vibration waveform of the reference blasthole and the activation coefficient of the other blastholes; The expansion subunit is used to linearly superpose and time-domain expand the recombined IMF corresponding to the vibration waveform of each blast hole to obtain the target IMF corresponding to the vibration waveform of each blast hole.

8. A computer program product, characterized in that It comprises computer-readable instructions, and when the computer-readable instructions are executed on an electronic device, the electronic device implements the method for determining the detonation time sequence of a blasthole as claimed in any one of claims 1 to 4.

9. An electronic device, characterized in that: The method comprises at least one processor and a memory connected to the processor, wherein: The memory is used to store computer programs; The processor is used to execute the computer program so that the electronic device can implement the method for determining the detonation time sequence of a blasthole as described in any one of claims 1 to 4.

10. A computer storage medium, characterized in that: The storage medium carries one or more computer programs, which, when executed by an electronic device, enable the electronic device to implement the method for determining the detonation time sequence of a blast hole as described in any one of claims 1 to 4.

Images