A method for determining a time sequence of initiating a blast hole and related equipment

By acquiring the vibration waveform of the borehole, determining the recombined IMF and peak vibration velocity, selecting a reference borehole, calculating the activation coefficient, and optimizing the borehole detonation time sequence, the problem of existing technologies being unable to fully consider complex factors is solved, thus improving the blasting effect.

CN119618009BActive Publication Date: 2026-07-24BAOTOU IRON & STEEL (GROUP) CO LTD
View PDF 1 Cites 0 Cited by

Patent Information

Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
BAOTOU IRON & STEEL (GROUP) CO LTD
Filing Date
2024-10-16
Publication Date
2026-07-24

AI Technical Summary

Technical Problem

Existing technologies cannot fully consider the complex factors in the blasting process when determining the time sequence of borehole detonation, resulting in poor blasting effects.

Method used

By acquiring the vibration waveforms of multiple boreholes, the reconstructed IMF corresponding to the vibration waveform of each borehole is determined, and the vibration velocity peak is extracted from the reconstructed IMF. A reference borehole is selected, and based on the vibration velocity peak of the reference borehole and the vibration velocity peak of other boreholes, the activation coefficient is calculated, linear superposition and time-domain extension are performed, and the detonation time sequence is combined according to Anderson's principle to optimize the detonation time of the borehole.

Benefits of technology

This enabled coordinated blasting between different blast holes, improving the blasting effect.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN119618009B_ABST
    Figure CN119618009B_ABST
Patent Text Reader

Abstract

The application provides a blast hole initiation time sequence determination method and related equipment, and relates to the field of blasting vibration control. The method acquires vibration waveforms of multiple blast holes; determines reorganized intrinsic mode components (IMFs) corresponding to the vibration waveforms of the blast holes; extracts more accurate vibration velocity peaks of the blast holes from the reorganized IMFs corresponding to the vibration waveforms of the blast holes; determines target IMFs corresponding to the vibration waveforms of the blast holes based on vibration velocity peaks and reorganized IMFs corresponding to the vibration waveforms of a reference blast hole and vibration velocity peaks of other blast holes; determines multiple superposition initiation time sequences; combines each superposition initiation time sequence with the target IMFs corresponding to the vibration waveforms of the blast holes to obtain multiple synthetic IMFs; extracts multiple synthetic vibration velocity peaks from the synthetic IMFs; and takes a time sequence corresponding to the smallest synthetic vibration velocity peak as the initiation time sequence of the blast hole. The scheme can realize coordinated blasting among the blast holes, optimizes the initiation time sequence, and is conducive to improving the blasting effect.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This application belongs to the field of blasting vibration control, and in particular relates to a method for determining the detonation time sequence of blast holes and related equipment. Background Technology

[0002] With the rapid development of blasting technology and the continuous expansion of blasting scale, large-scale deep-hole blasting involves numerous boreholes distributed over a wide area. The detonation time sequence of these boreholes is not only a key means to control blasting vibrations and improve blasting efficiency and safety, but also a significant driving force for advancing blasting technology to higher levels. Therefore, determining the detonation time sequence of the boreholes is of paramount importance.

[0003] Currently, methods based on empirical formulas or experimental summaries are commonly used to determine the detonation time sequence of boreholes.

[0004] However, empirical formulas and experimental 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 and related equipment for determining the detonation time sequence of blast holes. In order to improve the blasting effect, the specific scheme is as follows:

[0006] The first aspect of this application provides a method for determining the detonation time sequence of a borehole, including:

[0007] Obtain vibration waveforms from multiple boreholes;

[0008] Based on the vibration waveform of each borehole, the reconstructed IMF corresponding to the vibration waveform of each borehole is determined, and the peak vibration velocity of each borehole is extracted from the reconstructed IMF corresponding to the vibration waveform of each borehole.

[0009] A reference borehole is determined from each of the boreholes. Based on the peak vibration velocity and the reconstructed IMF corresponding to the vibration waveform of the reference borehole, as well as the peak vibration velocity of other boreholes, the target IMF corresponding to the vibration waveform of each of the boreholes is determined. The other boreholes are the boreholes other than the reference borehole.

[0010] Multiple superimposed detonation time sequences are determined, and each superimposed detonation time sequence is combined with the target IMF corresponding to the vibration waveform of each borehole according to Anderson's principle to obtain multiple synthetic IMFs;

[0011] Multiple synthetic vibrational velocity peaks are extracted from each of the synthetic IMFs, and the time series corresponding to the smallest synthetic vibrational velocity peak is used as the detonation time series of the multiple boreholes.

[0012] In one possible implementation, determining the reconstructed IMF corresponding to the vibration waveform of each of the boreholes, based on the vibration waveform of each borehole, includes:

[0013] The vibration waveforms of each borehole are decomposed to obtain the intrinsic mode components (IMFs) corresponding to the vibration waveforms of each borehole.

[0014] The number of layers of the IMF corresponding to the vibration waveform of each borehole is obtained. Based on the IMF with the fewest layers, the vibration waveform of each borehole is re-decomposed to obtain the reconstructed IMF corresponding to the vibration waveform of each borehole.

[0015] In one possible implementation, determining the target IMF corresponding to the vibration waveform of each borehole based on the reconstructed IMF corresponding to the peak vibration velocity and vibration waveform of the reference borehole, as well as the peak vibration velocity of other boreholes, includes:

[0016] Based on the peak vibration velocity of the reference borehole and the peak vibration velocity of the other boreholes, the activation coefficients of the other boreholes are calculated.

[0017] Based on the reconstructed IMF corresponding to the vibration waveform of the reference borehole and the activation coefficients of the other boreholes, the reconstructed IMF corresponding to the vibration waveform of the other boreholes is calculated.

[0018] Linear superposition and time-domain extension are performed on the reconstructed IMFs corresponding to the vibration waveforms of each borehole to obtain the target IMFs corresponding to the vibration waveforms of each borehole.

[0019] In one possible implementation, determining multiple superimposed detonation time series includes:

[0020] Based on the propagation velocity of the vibration waveform of each of the blast holes and the target distance of each of the blast holes, the detonation time difference sequence is calculated, where the target distance is the distance between the position of the blast hole and the blasting vibration monitoring position.

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

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

[0023] The acquisition unit is used to acquire the vibration waveforms of multiple boreholes;

[0024] The first determining unit is used to determine the reconstructed IMF corresponding to the vibration waveform of each of the boreholes based on the vibration waveform of each of the boreholes, and extract the peak vibration velocity of each of the boreholes from the reconstructed IMF corresponding to the vibration waveform of each of the boreholes.

[0025] The second determining unit is used to determine a reference borehole from each of the boreholes, and to determine the target IMF corresponding to the vibration waveform of each of the boreholes based on the peak vibration velocity and the reconstructed IMF corresponding to the vibration waveform of the reference borehole, as well as the peak vibration velocity of other boreholes. The other boreholes are the boreholes in each of the boreholes excluding the reference borehole.

[0026] The third determining unit is used to determine multiple superimposed detonation time sequences, and according to Anderson's principle, combine each superimposed detonation time sequence with the target IMF corresponding to the vibration waveform of each borehole to obtain multiple synthetic IMFs.

[0027] The fourth determining unit is used to extract multiple synthetic vibration velocity peaks from each of the synthetic IMFs, and to take the time sequence corresponding to the smallest synthetic vibration velocity peak as the detonation time sequence of the multiple boreholes.

[0028] In one possible implementation, the first determining unit for determining the reconstructed IMF corresponding to the vibration waveform of each of the boreholes, based on the vibration waveform of each borehole, includes:

[0029] The decomposition subunit is used to decompose the vibration waveform of each of the boreholes to obtain the intrinsic mode components (IMFs) corresponding to the vibration waveform of each of the boreholes.

[0030] The first determining subunit is used to obtain the number of layers of the IMF corresponding to the vibration waveform of each borehole, and to re-decompose the vibration waveform of each borehole based on the IMF with the fewest layers to obtain the reconstructed IMF corresponding to the vibration waveform of each borehole.

[0031] In one possible implementation, the second determining unit, which determines the target IMF corresponding to the vibration waveform of each borehole based on the peak vibration velocity and the reconstructed IMF corresponding to the vibration waveform of the reference borehole, and the peak vibration velocity of other boreholes, includes:

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

[0033] The second calculation subunit is used to calculate the reconstructed IMF corresponding to the vibration waveform of the other boreholes based on the reconstructed IMF corresponding to the vibration waveform of the reference borehole and the activation coefficient of the other boreholes.

[0034] An extended subunit is used to linearly superimpose and time-domain extend the reconstructed IMFs corresponding to the vibration waveforms of each borehole to obtain the target IMFs corresponding to the vibration waveforms of each borehole.

[0035] A third aspect of this application provides a computer program product including computer-readable instructions that, when executed on an electronic device, cause the electronic device to implement the method for determining the detonation time sequence of a borehole as described in the first aspect or any implementation thereof.

[0036] A fourth aspect of this application provides an electronic device, including 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 borehole in the first aspect or any implementation thereof.

[0039] The fifth aspect of this application provides a computer storage medium carrying one or more computer programs that, when executed by an electronic device, enable the electronic device to perform a method for determining the detonation time sequence of a borehole as described in the first aspect or any implementation thereof.

[0040] Based on the above technical solution, this application provides a method and related equipment for determining the detonation time sequence of blast holes. The method acquires the locations of multiple blast holes and blasting vibration monitoring locations within the blast zone, and collects the vibration waveforms of each blast hole. Based on the vibration waveforms of each blast hole, it determines the reconstructed IMF corresponding to the vibration waveform of each blast hole, and extracts more accurate vibration velocity peak values ​​from the reconstructed IMFs corresponding to the vibration waveforms of each blast hole. It determines a reference blast hole from among the blast holes, and based on the vibration velocity peak value of the reference blast hole, the reconstructed IMF corresponding to the vibration waveform of the reference blast hole, and the vibration velocity peak values ​​of other blast holes, it determines the target IMF corresponding to the vibration waveform of each blast hole. The other blast holes are those excluding the reference blast hole. It determines multiple superimposed detonation time sequences, and according to Anderson's principle, combines each superimposed detonation time sequence with the target IMF corresponding to the vibration waveform of each blast hole to obtain multiple synthetic IMFs. It extracts multiple synthetic vibration velocity peak values ​​from each synthetic IMF, and uses the time sequence corresponding to the smallest synthetic vibration velocity peak value as the detonation time sequence of multiple blast holes. This scheme optimizes the detonation time sequence of blast holes, enabling coordinated blasting between different blast holes and improving blasting effectiveness. Attached Figure Description

[0041] Various other advantages and benefits will become apparent to those skilled in the art upon reading the following detailed description of preferred embodiments. The accompanying drawings are for illustrative purposes only and are not intended to limit the scope of this application. Furthermore, the same reference numerals denote the same parts throughout the drawings. In the drawings:

[0042] Figure 1 This is a flowchart illustrating a borehole detonation time sequence method disclosed in an embodiment of this application;

[0043] Figure 2 This is a schematic diagram of the structure of a borehole detonation time sequence device disclosed in an embodiment of this application;

[0044] Figure 3 This is a hardware structure block diagram of a borehole detonation time sequence device disclosed in an embodiment of this application. Detailed Implementation

[0045] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0046] The terms "first," "second," etc., used in the specification, claims, and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such terms are interchangeable where appropriate; this is merely a way of distinguishing objects with the same attributes in the embodiments of this application. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion, so that a process, method, system, product, or apparatus that comprises a series of elements is not necessarily limited to those elements, but may include other elements not explicitly listed or inherent to those processes, methods, products, or apparatuses.

[0047] To improve blasting effects, this application provides a method for determining a regression model of coal coke thermal properties. The method for reflecting implicit orders in trading data provided in this application is further explained in detail below with reference to the accompanying drawings and specific embodiments.

[0048] Please see the appendix Figure 1 , Figure 1 This is a flowchart illustrating a method for determining the detonation time sequence of a borehole, provided in an embodiment of this application. The method may include the following steps:

[0049] Step S101: Obtain the vibration waveforms of multiple boreholes.

[0050] In this application, sensors (such as vibration sensors or accelerometers) can be used to record the vibration waveforms generated by each borehole during the blasting process, which contain vibration signals transmitted from each borehole.

[0051] Step S102: Based on the vibration waveform of each borehole, determine the reconstructed intrinsic mode component (IMF) corresponding to the vibration waveform of each borehole, and extract the peak vibration velocity of each borehole from the reconstructed IMF corresponding to the vibration waveform of each borehole.

[0052] It should be noted that IMF (Intrinsic Mode Function) is a fundamental concept in empirical mode decomposition methods, representing the local characteristic scale of a signal. Each IMF satisfies two conditions: the number of extrema and the number of zero-crossings must be equal or differ by at most one in the entire dataset; and at any point, the average of the upper envelope formed by local maxima and the lower envelope formed by local minima is zero.

[0053] In this application, determining the reconstructed IMF corresponding to the vibration waveform of each borehole based on the vibration waveform of each borehole includes: decomposing the vibration waveform of each borehole to obtain the IMF corresponding to the vibration waveform of each borehole; obtaining the number of layers of the IMF corresponding to the vibration waveform of each borehole; and re-decomposing the vibration waveform of each borehole based on the IMF with the fewest layers to obtain the reconstructed IMF corresponding to the vibration waveform of each borehole.

[0054] Complex vibration waveforms can be decomposed into a series of relatively simple IMFs using Empirical Mode Decomposition (EMD) or its improved methods (such as ensemble EMD). These IMFs 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 differ. The IMF with the fewest layers can be selected as a benchmark to re-decompose the vibration waveforms of all boreholes. Then, the peak vibration velocities are extracted from the reconstructed IMFs of each borehole; these peaks reflect the maximum vibration intensity generated during the borehole explosion.

[0055] By obtaining the charge amount and target distance of the borehole, a functional relationship can be established between the peak vibration velocity of the borehole and the charge amount and target distance. The target distance indicates the distance between the borehole location and the vibration monitoring point. See the following formula for details:

[0056]

[0057] In the formula: V k Let $k$ be the peak velocity of the $k$-th recombined IMF, where $k \in {1, 2, ..., K}. min}, K min The minimum number of layers, Q is the charge amount in the borehole, R is the distance between the borehole location and the vibration monitoring point location, and α is the minimum number of layers. k β k These are parameters related to the properties of explosives, rock properties, and topography and geological conditions.

[0058] Step S103: Determine the reference borehole from each borehole. Based on the peak vibration velocity and the reconstructed IMF corresponding to the vibration waveform of the reference borehole, as well as the peak vibration velocity of other boreholes, determine the target IMF corresponding to the vibration waveform of each borehole. The other boreholes are the boreholes other than the reference borehole.

[0059] In this application, key vibration control locations are selected from vibration monitoring points. The borehole corresponding to this location can be used as a reference borehole. Based on the peak vibration velocity of the reference borehole and the peak vibration velocity of other boreholes, the activation coefficients of the other boreholes are calculated. Based on the reconstructed IMF corresponding to the vibration waveform of the reference borehole and the activation coefficients of other boreholes, the reconstructed IMF corresponding to the vibration waveform of the other boreholes is calculated. The reconstructed IMFs corresponding to the vibration waveforms of each borehole are linearly superimposed and time-domain extended to obtain the target IMF corresponding to the vibration waveform of each borehole.

[0060] Specifically, based on the assumption that the functional relationship between the peak vibration velocity of a borehole and the charge amount and target distance holds true at any identical moment within the duration of the vibration waveform, the ratio of the peak vibration velocity of other boreholes to that of the reference borehole is calculated and used as an activation coefficient. This activation coefficient can be viewed as a scaling or adjustment factor of the vibration waveform of other boreholes relative to the reference borehole vibration waveform. See the following formula for details:

[0061]

[0062] In the formula: A j,k (Q0,R0,Q j ,R j ) represents the activation coefficient of the k-th layer IMF corresponding to the vibration waveform of the j-th borehole, v j,k v0 represents the peak vibration velocity of other boreholes; Q0 represents the peak vibration velocity of the reference borehole; R0 represents the charge amount of the reference borehole; and Q represents the distance from the reference borehole to the key vibration control location. j R represents the charge amount for the j-th other borehole. j Let be the distance from the j-th other borehole to the critical location for vibration control.

[0063] Based on the reconstructed IMF corresponding to the vibration waveform of the reference borehole, the reconstructed IMF corresponding to the vibration waveforms of other boreholes is obtained by multiplying it by the activation coefficient of other boreholes. See the following formula for details:

[0064]

[0065] In the formula: f' j,k (t) represents the k-th layer IMF of the j-th other borehole, f 0,k (t) represents the k-th layer IMF of the reference borehole.

[0066] The IMFs of each layer of the other boreholes are combined to obtain the IMF of the other boreholes. See the following formula for details:

[0067]

[0068] In the formula: f' j (t) represents the IMF of the j-th other borehole.

[0069] By linearly superimposing and time-domain extending the reconstructed IMFs corresponding to the vibration waveforms of each borehole, the target IMF corresponding to the vibration waveform of each borehole is obtained. The specific formula can be found in the following formula:

[0070]

[0071] In the formula: f j f(t) represents the target IMF of the j-th other borehole in the full time domain, and tdur represents the duration of the reference borehole f0(t).

[0072] Step S104: Determine multiple superimposed detonation time sequences, and combine each superimposed detonation time sequence with the target IMF corresponding to the vibration waveform of each borehole according to Anderson's principle to obtain multiple synthetic IMFs.

[0073] In this application, the initiation time difference sequence can be calculated based on the distance from the reference borehole to the vibration control key position and the distance from the j-th other borehole to the vibration control key position, as well as the propagation velocity of the vibration waveform of each borehole. Multiple superimposed initiation time sequences are then calculated based on multiple preset initiation time sequences and initiation 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] In the formula: t j-0 c is the time difference between the vibration wave generated by the j-th other borehole and the vibration wave generated by the reference borehole, which travels to the critical vibration control location. v Let t' be the longitudinal wave velocity. j t represents the actual time t takes for the j-th other borehole to participate in the superposition calculation. j Let be the detonation time of the j-th other blast hole.

[0077] According to Anderson's principle, each superimposed detonation time series is combined with the target IMF corresponding to the vibration waveform of each borehole to obtain multiple composite IMFs. See the following formula for details:

[0078]

[0079] In the formula, f j (t-t'j) represents the target IMF corresponding to the vibration waveform of the j-th other borehole; F(t) represents the composite IMF.

[0080] Furthermore, by changing the detonation time sequence, the synthetic IMF corresponding to each detonation time sequence is obtained sequentially.

[0081] Step S105: Extract multiple synthetic vibration velocity peaks from each synthetic IMF, and use the time series corresponding to the smallest synthetic vibration velocity peak as the detonation time series of multiple boreholes.

[0082] In this application, multiple composite velocity peaks Vmax(t,{t) are extracted from each composite IMF. j +t j-0 The peak values ​​of the synthesized vibration velocity are the positive and negative maximum velocities in the synthesized waveform. Then, from multiple peak values ​​of the synthesized vibration velocity Vmax(t,{t... j +t j-0 Determine the minimum peak value of the combined vibrational velocity V in}). 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 peak composite vibration velocity was used as the detonation time series for multiple boreholes.

[0085] In summary, this application provides a method for determining the detonation time sequence of blast holes. This method acquires the locations of multiple blast holes and blasting vibration monitoring locations within the blast zone, and collects the vibration waveforms of each blast hole. Based on the vibration waveforms of each blast hole, a reconstructed IMF corresponding to the vibration waveform of each blast hole is determined, and a more accurate vibration velocity peak value is extracted from the reconstructed IMF corresponding to the vibration waveform of each blast hole. A reference blast hole is determined from among the blast holes. Based on the vibration velocity peak value of the reference blast hole, the reconstructed IMF corresponding to the vibration waveform of the reference blast hole, and the vibration velocity peak values ​​of other blast holes, a target IMF corresponding to the vibration waveform of each blast hole is determined. The other blast holes are those excluding the reference blast hole. Multiple superimposed detonation time sequences are determined, and according to Anderson's principle, each superimposed detonation time sequence is combined with the target IMF corresponding to the vibration waveform of each blast hole to obtain multiple synthetic IMFs. Multiple synthetic vibration velocity peak values ​​are extracted from each synthetic IMF, and the time sequence corresponding to the smallest synthetic vibration velocity peak value is used as the detonation time sequence of multiple blast holes. This scheme optimizes the detonation time sequence of blast holes, enabling coordinated blasting between different blast holes and improving blasting effectiveness.

[0086] The methods described in the above-disclosed embodiments of this application are detailed in terms of the methods. The methods of this application can be implemented by various forms of devices. Therefore, this application also discloses a device for determining the detonation time sequence of a blast hole. Specific embodiments are given below for detailed description.

[0087] Please see the appendix Figure 2 , Figure 3 This is a schematic diagram of a device for determining the detonation time sequence of a borehole, as disclosed in an embodiment of this application. The device includes:

[0088] Acquisition unit 11 is used to acquire the vibration waveforms of multiple boreholes.

[0089] The first determining unit 12 is used to determine the reconstructed IMF corresponding to the vibration waveform of each of the boreholes based on the vibration waveform of each of the boreholes, and to extract the peak vibration velocity of each of the boreholes from the reconstructed IMF corresponding to the vibration waveform of each of the boreholes.

[0090] The second determining unit 13 is used to determine a reference borehole from each of the boreholes, and to determine the target IMF corresponding to the vibration waveform of each borehole based on the peak vibration velocity and the reconstructed IMF corresponding to the vibration waveform of the reference borehole, as well as the peak vibration velocity of other boreholes, wherein the other boreholes are the boreholes other than the reference borehole.

[0091] The third determining unit 14 is used to determine multiple superimposed detonation time sequences, and according to Anderson's principle, combine each superimposed detonation time sequence with the target IMF corresponding to the vibration waveform of each borehole to obtain multiple synthetic IMFs.

[0092] The fourth determining unit 15 is used to extract multiple synthetic vibration velocity peaks from each of the synthetic IMFs, and to take the time sequence corresponding to the smallest synthetic vibration velocity peak as the detonation time sequence of the multiple boreholes.

[0093] As one possible implementation, the first determining unit 12 for determining the reconstructed IMF corresponding to the vibration waveform of each of the boreholes, based on the vibration waveform of each borehole, includes:

[0094] The decomposition subunit is used to decompose the vibration waveform of each of the boreholes to obtain the intrinsic mode components (IMFs) corresponding to the vibration waveform of each borehole.

[0095] The first determining subunit is used to obtain the number of layers of the IMF corresponding to the vibration waveform of each borehole, and to re-decompose the vibration waveform of each borehole based on the IMF with the fewest layers to obtain the reconstructed IMF corresponding to the vibration waveform of each borehole.

[0096] As one possible implementation, the second determining unit 13, which determines the target IMF corresponding to the vibration waveform of each borehole based on the reconstructed IMF corresponding to the peak vibration velocity and vibration waveform of the reference borehole, and the peak vibration velocity of other boreholes, includes:

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

[0098] The second calculation subunit is used to calculate the reconstructed IMF corresponding to the vibration waveform of the other boreholes based on the reconstructed IMF corresponding to the vibration waveform of the reference borehole and the activation coefficient of the other boreholes.

[0099] An extended subunit is used to linearly superimpose and time-domain extend the reconstructed IMFs corresponding to the vibration waveforms of each borehole to obtain the target IMFs corresponding to the vibration waveforms of each borehole.

[0100] As one possible implementation, the third determining unit 14 for determining multiple superimposed detonation time sequences includes:

[0101] The third calculation subunit is used to calculate the detonation time difference sequence based on the propagation velocity of the vibration waveform of each of the blast holes and the target distance of each of the blast holes, wherein 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 multiple superimposed detonation time sequences based on multiple preset detonation time sequences and the detonation time difference sequence.

[0103] This application also provides an electronic device in its embodiments. (See reference...) Figure 3 The diagram illustrates a structural schematic suitable for implementing the electronic device in the embodiments of this application. The electronic device in the embodiments of this application may include, but is not limited to, fixed terminals such as mobile phones, laptops, PDAs (personal digital assistants), PADs (tablet computers), desktop computers, etc. Figure 3 The electronic device shown is merely an example and should not impose any limitation on the functionality and scope of use of the embodiments of this application.

[0104] like Figure 3 As shown, the electronic device may include a processing unit (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 into a random access memory (RAM) 303. When the electronic device is powered on, the RAM 303 also stores various programs and data required for the operation of the electronic device. The processing unit 301, ROM 302, and RAM 303 are interconnected via a bus 304. An input / output (I / O) interface 303 is also connected to the bus 304.

[0105] Typically, the following devices can be connected to I / O interface 305: input devices 306 including, for example, touchscreens, touchpads, keyboards, mice, cameras, microphones, accelerometers, gyroscopes, etc.; output devices 307 including, for example, liquid crystal displays (LCDs), speakers, vibrators, etc.; storage devices 308 including, for example, memory cards, hard drives, etc.; and communication devices 309. Communication device 309 allows electronic devices to communicate wirelessly or wiredly with other devices to exchange data. Although Figure 3 Electronic devices with various devices are shown, but it should be understood that it is not required to implement or have all of the devices shown. More or fewer devices may be implemented or have alternatively.

[0106] This application also provides a computer program product including computer-readable instructions, which, when executed on an electronic device, cause the electronic device to implement any of the methods for determining the detonation time sequence of boreholes provided in this application.

[0107] This application also provides a computer-readable storage medium carrying one or more computer programs. When the one or more computer programs are executed by an electronic device, the electronic device can implement any of the methods for determining the detonation time sequence of blast holes provided in this application.

[0108] It should also be noted that the device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. In addition, in the device embodiment drawings provided in this application, the connection relationship between modules indicates that they have a communication connection, which can be implemented as one or more communication buses or signal lines.

[0109] Through the above description of the embodiments, those skilled in the art can clearly understand that this application can be implemented by means of software plus necessary general-purpose hardware, or it can be implemented by special-purpose hardware including application-specific integrated circuits, special-purpose CPUs, special-purpose memory, special-purpose components, etc. Generally, any function performed by a computer program can be easily implemented by corresponding hardware, and the specific hardware structure used to implement the same function can also be diverse, such as analog circuits, digital circuits, or special-purpose circuits. However, for this application, software program implementation is more often the preferred implementation method. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product is stored in a readable storage medium, such as a computer floppy disk, USB flash drive, mobile hard disk, ROM, RAM, magnetic disk, or optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, training equipment, or network device, etc.) to execute the methods described in the various embodiments of this application.

[0110] In the above embodiments, implementation can be achieved, in whole or in part, through software, hardware, firmware, or any combination thereof. When implemented in software, it can be implemented, in whole or in part, as 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, all or part of the processes or functions described in the embodiments of this application are generated. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions may be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, the computer instructions may be transmitted from one website, computer, training device, or data center to another website, computer, training device, or data center via wired (e.g., coaxial cable, fiber optic, digital subscriber line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) means. 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 or data center that integrates one or more available media. The available media may be magnetic media (e.g., floppy disks, hard disks, magnetic tapes), optical media (e.g., DVDs), or semiconductor media (e.g., solid-state drives (SSDs)).

Claims

1. A method for determining the detonation time sequence of a borehole, characterized in that, include: Obtain vibration waveforms from multiple boreholes; Based on the vibration waveform of each borehole, the reconstructed intrinsic mode component (IMF) corresponding to the vibration waveform of each borehole is determined, and the peak vibration velocity of each borehole is extracted from the reconstructed IMF corresponding to the vibration waveform of each borehole. A reference borehole is determined from each of the boreholes. Based on the peak vibration velocity and the reconstructed IMF corresponding to the vibration waveform of the reference borehole, as well as the peak vibration velocity of other boreholes, the target IMF corresponding to the vibration waveform of each of the boreholes is determined. The other boreholes are the boreholes other than the reference borehole. The method involves determining the target IMF corresponding to the vibration waveform of each borehole based on the peak vibration velocity and the reconstructed IMF corresponding to the vibration waveform of the reference borehole, as well as the peak vibration velocity of other boreholes. This includes: Based on the peak vibration velocity of the reference borehole and the peak vibration velocity of the other boreholes, the activation coefficients of the other boreholes are calculated. Specifically, based on the assumption that the functional relationship between the peak vibration velocity of the borehole and the charge amount and target distance holds true at any identical moment within the duration of the vibration waveform, the ratio of the peak vibration velocity of other boreholes to that of the reference borehole is calculated and used as an activation coefficient. The activation coefficient can be regarded as a scaling or adjustment factor of the vibration waveform of other boreholes relative to the vibration waveform of the reference borehole. The specific formula is as follows: In the formula: A j,k (Q0,R0,Q j ,R j ) represents the activation coefficient of the k-th layer IMF corresponding to the vibration waveform of the j-th borehole, v j,k v0 represents the peak vibration velocity of other boreholes; Q0 represents the peak vibration velocity of the reference borehole; R0 represents the charge amount of the reference borehole; and Q represents the distance from the reference borehole to the key vibration control location. j R represents the charge amount for the j-th other borehole. j Let be the distance from the j-th other borehole to the critical vibration control location; Based on the reconstructed IMF corresponding to the vibration waveform of the reference borehole and the activation coefficients of the other boreholes, the reconstructed IMF corresponding to the vibration waveform of the other boreholes is calculated. Specifically, based on the reconstructed IMF corresponding to the vibration waveform of the reference borehole, the reconstructed IMF corresponding to the vibration waveform of other boreholes is obtained by multiplying it by the activation coefficient of other boreholes. The specific formula is as follows: In the formula: f j,k (t) represents the k-th layer IMF of the j-th other borehole, f 0,k (t) represents the k-th layer IMF of the reference borehole; Linear superposition and time-domain extension are performed on the reconstructed IMFs corresponding to the vibration waveforms of each borehole to obtain the target IMFs corresponding to the vibration waveforms of each borehole. Specifically, the IMFs of each layer of the other borehole are combined to obtain the IMF of the other borehole. The specific formula is as follows: In the formula: f' j (t) represents the IMF of the j-th other borehole; By linearly superimposing and time-domain extending the reconstructed IMFs corresponding to the vibration waveforms of each borehole, the target IMF corresponding to the vibration waveform of each borehole is obtained. The specific formula is as follows: In the formula: f j (t) represents the target IMF of the j-th other borehole in the entire time domain, and tdur represents the duration of the reference borehole f0(t). Multiple superimposed detonation time sequences are determined, and each superimposed detonation time sequence is combined with the target IMF corresponding to the vibration waveform of each borehole according to Anderson's principle to obtain multiple synthetic IMFs; Multiple synthetic vibrational velocity peaks are extracted from each of the synthetic IMFs, and the time series corresponding to the smallest synthetic vibrational velocity peak is used as the detonation time series of the multiple boreholes.

2. The method for determining the detonation time sequence of a borehole according to claim 1, characterized in that, The determination of the reconstructed IMF corresponding to the vibration waveform of each of the boreholes, based on the vibration waveform of each borehole, includes: The vibration waveforms of each borehole are decomposed to obtain the IMF corresponding to the vibration waveforms of each borehole. The number of layers of the IMF corresponding to the vibration waveform of each borehole is obtained. Based on the IMF with the fewest layers, the vibration waveform of each borehole is re-decomposed to obtain the reconstructed IMF corresponding to the vibration waveform of each borehole.

3. The method for determining the detonation time sequence of a borehole according to claim 1, characterized in that, The determination of multiple superimposed detonation time sequences includes: Based on the propagation velocity of the vibration waveform of each of the blast holes and the target distance of each of the blast holes, the detonation time difference sequence is calculated, where the target distance is the distance between the position of the blast hole and the blasting vibration monitoring position. Based on multiple preset detonation time sequences and the detonation time difference sequence, multiple superimposed detonation time sequences are calculated.

4. A device for determining the detonation time sequence of a blast hole, characterized in that, include: The acquisition unit is used to acquire the vibration waveforms of multiple boreholes; The first determining unit is used to determine the reconstructed IMF corresponding to the vibration waveform of each of the boreholes based on the vibration waveform of each of the boreholes, and extract the peak vibration velocity of each of the boreholes from the reconstructed IMF corresponding to the vibration waveform of each of the boreholes. The second determining unit is used to determine a reference borehole from each of the boreholes, and to determine the target IMF corresponding to the vibration waveform of each of the boreholes based on the peak vibration velocity and the reconstructed IMF corresponding to the vibration waveform of the reference borehole, as well as the peak vibration velocity of other boreholes. The other boreholes are the boreholes in each of the boreholes excluding the reference borehole. The second determining unit, which determines the target IMF corresponding to the vibration waveform of each borehole based on the peak vibration velocity and the reconstructed IMF corresponding to the vibration waveform of the reference borehole, and the peak vibration velocity of other boreholes, includes: The first calculation subunit is used to calculate the activation coefficient of the other blast holes based on the peak vibration velocity of the reference blast hole and the peak vibration velocity of the other blast holes. Specifically, based on the assumption that the functional relationship between the peak vibration velocity of the borehole and the charge amount and target distance holds true at any identical moment within the duration of the vibration waveform, the ratio of the peak vibration velocity of other boreholes to that of the reference borehole is calculated and used as an activation coefficient. The activation coefficient can be regarded as a scaling or adjustment factor of the vibration waveform of other boreholes relative to the vibration waveform of the reference borehole. The specific formula is as follows: In the formula: A j,k (Q0,R0,Q j ,R j ) represents the activation coefficient of the k-th layer IMF corresponding to the vibration waveform of the j-th borehole, v j,k v0 represents the peak vibration velocity of other boreholes; Q0 represents the peak vibration velocity of the reference borehole; R0 represents the charge amount of the reference borehole; and Q represents the distance from the reference borehole to the key vibration control location. j R represents the charge amount for the j-th other borehole. j Let be the distance from the j-th other borehole to the critical vibration control location; The second calculation subunit is used to calculate the reconstructed IMF corresponding to the vibration waveform of the other boreholes based on the reconstructed IMF corresponding to the vibration waveform of the reference borehole and the activation coefficient of the other boreholes. Specifically, based on the reconstructed IMF corresponding to the vibration waveform of the reference borehole, the reconstructed IMF corresponding to the vibration waveform of other boreholes is obtained by multiplying it by the activation coefficient of other boreholes. The specific formula is as follows: In the formula: f j,k (t) represents the k-th layer IMF of the j-th other borehole, f 0,k (t) represents the k-th layer IMF of the reference borehole; An extended subunit is used to linearly superimpose and time-domain extend the reconstructed IMFs corresponding to the vibration waveforms of each borehole to obtain the target IMFs corresponding to the vibration waveforms of each borehole. Specifically, the IMFs of each layer of the other borehole are combined to obtain the IMF of the other borehole. The specific formula is as follows: In the formula: f' j (t) represents the IMF of the j-th other borehole; By linearly superimposing and time-domain extending the reconstructed IMFs corresponding to the vibration waveforms of each borehole, the target IMF corresponding to the vibration waveform of each borehole is obtained. The specific formula is as follows: In the formula: f j (t) represents the target IMF of the j-th other borehole in the entire time domain, and tdur represents the duration of the reference borehole f0(t). The third determining unit is used to determine multiple superimposed detonation time sequences, and according to Anderson's principle, combine each superimposed detonation time sequence with the target IMF corresponding to the vibration waveform of each borehole to obtain multiple synthetic IMFs. The fourth determining unit is used to extract multiple synthetic vibration velocity peaks from each of the synthetic IMFs, and to take the time sequence corresponding to the smallest synthetic vibration velocity peak as the detonation time sequence of the multiple boreholes.

5. The device for determining the detonation time sequence of a borehole according to claim 4, characterized in that, The first determining unit for determining the reconstructed IMF corresponding to the vibration waveform of each of the boreholes, based on the vibration waveform of each borehole, includes: The decomposition subunit is used to decompose the vibration waveform of each of the boreholes to obtain the intrinsic mode components (IMFs) corresponding to the vibration waveform of each of the boreholes. The first determining subunit is used to obtain the number of layers of the IMF corresponding to the vibration waveform of each borehole, and to re-decompose the vibration waveform of each borehole based on the IMF with the fewest layers to obtain the reconstructed IMF corresponding to the vibration waveform of each borehole.

6. A computer program product, characterized in that, It includes computer-readable instructions that, when executed on an electronic device, cause the electronic device to implement the method for determining the detonation time sequence of a borehole as described in any one of claims 1 to 3.

7. An electronic device, characterized in that, It includes 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 to enable the electronic device to implement the method for determining the detonation time sequence of the borehole as described in any one of claims 1 to 3.

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