Electronic detonator network synchronous detonation control system and method based on topological optimization

Abstract

The invention relates to the technical field of detonation control, in particular to an electronic detonator network synchronous detonation control system and method based on topological optimization. The method comprises the steps that information of a to-be-blasted area where electronic detonator network topological nodes are located is collected, and inter-hole delay of the to-be-blasted area and the optimal explosive loading amount of blast holes are determined based on blasting vibration prior data; geological analysis calculation is conducted on the geological parameter data to obtain a stability index of the to-be-blasted area, and an explosion performance index is obtained through calculation based on the optimal explosive loading amount of the blast hole and the initial explosion source depth of the to-be-blasted area; whether the number of topological nodes of the to-be-blasted area is adjusted or not is judged based on the explosion performance index and the stability index, so that a detonation control scheme is generated; and the detonation control scheme is executed to execute synchronous detonation control of the electronic detonator, the detonation effect can be optimized, the synchronous detonation accuracy of the electronic detonator can be improved, it is ensured that the detonation effect meets the expectation, and meanwhile construction safety is guaranteed.

Description

Technical Field

[0001] This invention relates to the field of detonation control technology, specifically to a synchronous detonation control system and method for electronic detonator networks based on topology optimization. Background Technology

[0002] The core advantage of electronic detonators lies in their use of electronic control modules to replace traditional chemical delay elements. Through programmable technology, they achieve millisecond-level delay accuracy and integrate automatic detection, fault diagnosis, and password command detonation functions. They are widely used in mining, tunneling, urban demolition, and underwater blasting.

[0003] Traditional electronic detonator initiation systems mostly adopt a wired networking method, with all detonators mounted on a busbar via long leads and controlled uniformly by an initiator. However, in large-scale mine or urban tunnel blasting, hundreds to thousands of electronic detonators need to be detonated simultaneously. The accuracy of synchronous detonation and the blasting effect are required to be high. Traditional electronic detonator initiation systems often do not take into account the charge amount of each electronic detonator and the inter-hole delay in the area to be blasted, resulting in poor blasting effect and low accuracy of synchronous detonation.

[0004] Therefore, there is an urgent need for a topology-optimized electronic detonator network synchronous initiation control system and method to solve the above problems. Summary of the Invention

[0005] The purpose of this invention is to provide a topology-optimized electronic detonator network synchronous detonation control system and method, which aims to solve the problems of poor blasting effect and low synchronous detonation accuracy in traditional electronic detonator detonation systems.

[0006] The objective of this invention can be achieved through the following technical solutions: On the one hand, a topology-optimized electronic detonator network synchronous initiation control system, the system includes: The blasting information acquisition unit is used to treat each blast hole in the area to be blasted as a topology node, construct an electronic detonator network based on the electronic detonators of the topology nodes, and acquire information about the area to be blasted where the topology nodes of the electronic detonator network are located. The information about the area to be blasted includes geological parameter data and prior data of blasting vibration. The blasting parameter determination unit is used to determine the inter-hole delay and the optimal charge amount of the borehole in the area to be blasted based on the prior data of blasting vibration. The data calculation unit is used to perform geological analysis and calculation on geological parameter data to obtain the stability index of the area to be blasted, and to calculate the explosive performance index based on the optimal charge amount of the blast hole and the initial blast source depth of the area to be blasted. The detonation control scheme generation unit is used to determine whether to adjust the number of topological nodes in the area to be blasted based on the explosion performance index and stability index. If not, the number of topological nodes is not adjusted and the detonation control scheme is generated with the initial topological nodes. If yes, the number of topological nodes is adjusted to generate the detonation control scheme. The synchronous detonation execution unit is used to execute the detonation control scheme to perform synchronous detonation control of the electronic detonator.

[0007] Furthermore, determining the inter-hole delay and optimal charge amount for the borehole in the area to be blasted based on prior data of blasting vibration specifically includes the following process: A cut-out area with the same environment as the area to be blasted was designed, and a single-hole, single-free-face blasting experiment was conducted within the cut-out area. Ground vibration curves for each hole were collected, with a one-to-one correspondence between the topological nodes of the cut-out area and the area to be blasted. The ground vibration curves were fitted using a Fourier series form. By analyzing the combined vibration velocity under different inter-hole delays for each borehole charge, the optimal charge amount and inter-hole delay for each borehole were obtained. Step 1: Perform multi-hole vibration superposition calculation on the ground vibration curve of a single hole according to the different time delays between candidate holes to obtain all synthetic waveform functions. Based on all synthetic waveform functions, obtain the upper limit of the maximum synthetic vibration velocity value. and lower limit ; Step 2, for any borehole charge Based on the limit value of the maximum synthesized vibration velocity obtained in step one This yields the set of maximum vibration velocity lower limits corresponding to all candidate charge quantities. Then, it identifies the values ​​in the set that are lower than the safe vibration velocity standard. The maximum value; Step 3: Select the set of maximum vibration velocity lower limit values ​​that are less than the safe vibration velocity standard. The maximum value of the selected charge amount is denoted as the optimal charge amount of the borehole, and the inter-hole delay corresponding to the optimal charge amount of the borehole is the inter-hole delay.

[0008] Furthermore, the collection of geological parameter data specifically includes the following processes: A three-dimensional coordinate system was established by selecting the seismic wave detector located at the center of the area to be blasted as the origin, and the position coordinates of each seismic wave detector were determined: , , , , , ,in, Given the distances of other seismic wave detectors to the central topology node, calculate the acoustic path from each seismic wave detector to the non-origin seismic wave detectors, and calculate the difference between each acoustic path. , 6. Calculation and Estimated time delay of sound waves According to the difference Delay estimate Determine the initial elevation angle of the target seismic wave detector installation location. and azimuth And based on the initial pitch angle and azimuth Determine the installation location coordinates of the target seismic wave detector, install the target seismic wave detector based on the installation location coordinates, and collect geological parameter data based on the target seismic wave detector.

[0009] Furthermore, the stability index of the area to be blasted is obtained through geological analysis and calculation of geological parameter data, specifically including the following processes: Based on geological parameter data, the root mean square amplitude, average amplitude, and average instantaneous frequency of seismic waves in the area to be blasted are obtained. Among them, the root mean square amplitude is used to identify amplitude anomalies or describe sequence stratigraphy and track seismic anomalies in the stratigraphy; the average amplitude is used to identify amplitude anomalies or sequence stratigraphy characteristics and identify lithology and lithological changes; the average instantaneous frequency is used to detect amplitude absorption anomalies and track changes in related frequency absorption characteristics caused by faults, lithology, or stratigraphy. Based on the root mean square amplitude, average amplitude, and average instantaneous frequency, the corresponding correlation coefficients are looked up in the correlation coefficient table to obtain the root mean square amplitude correlation coefficient. Correlation coefficient of average amplitude Correlation coefficient with average instantaneous frequency The correlation coefficient table stores the correlation coefficients corresponding to the ranges of the root mean square amplitude, average amplitude, and average instantaneous frequency. The root mean square amplitude correlation coefficient Correlation coefficient of average amplitude Correlation coefficient with average instantaneous frequency After normalization, the stability index of the area to be blasted is obtained by inputting it into the correlation formula. The correlation formula is as follows: in, , and This is the weight, and its value is set by the system.

[0010] Furthermore, the calculation of the explosive performance index based on the optimal charge amount of the borehole and the initial blast source depth of the area to be blasted specifically includes the following process: Determine the chemical reaction kinetics and specific thermal decomposition kinetic parameters of electronic detonator explosives: reaction order n, activation energy. Pre-exponential frequency factor and the heat of reaction Q per unit mass; Obtain the optimal charge M for the borehole, the number of topological nodes K, and the initial blast source depth of the area to be blasted. ; Obtaining the universal molar gas constant The heat released per unit mass of reactant Thermal conductivity λ, density of reactants d, the radius of the bottom surface of the reactant; c, the specific heat capacity of the reactant at constant volume. After normalization, the thermal conductivity parameter W was calculated as follows: ; After normalization, the critical explosion temperature was calculated based on the thermal conductivity parameter W. : ; After normalization, based on the critical explosion temperature The thermal explosion performance index was calculated. : .

[0011] Furthermore, determining whether to adjust the initial blast depth based on the explosion performance index and stability index specifically includes the following process: Determine whether the explosive performance index is greater than the stability index. If so, the initial blast source depth is deemed to meet the blasting conditions, and no adjustment is made to the initial blast source depth. If not, the initial blast source depth is deemed not to meet the blasting conditions, and adjustment is required.

[0012] Furthermore, generating the detonation control scheme without adjusting the number of topology nodes and using the initial topology nodes specifically includes the following process: A detonation control scheme is generated based on the initial set number of topology nodes, initial detonation source depth, inter-hole delay of topology nodes, and detonation sequence of topology nodes.

[0013] Furthermore, adjusting the number of topology nodes to generate the detonation control scheme specifically includes the following processes: Increase the number of topology nodes until the explosion performance index is greater than the stability index, and generate an initiation control scheme based on the adjusted number of topology nodes, initial blast source depth, inter-hole delay of topology nodes, and initiation sequence of topology nodes.

[0014] On the other hand, a synchronous detonation control method for electronic detonator networks based on topology optimization is provided, which includes: Each blast hole in the area to be blasted is taken as a topological node. An electronic detonator network is constructed based on the electronic detonators of the topological nodes. Information on the area to be blasted where the topological nodes of the electronic detonator network are located is collected. The information on the area to be blasted includes geological parameter data and prior data of blasting vibration. The inter-hole delay and optimal charge amount of the blast hole in the area to be blasted are determined based on prior data of blasting vibration. The stability index of the area to be blasted is obtained by geological analysis and calculation of geological parameter data, and the explosive performance index is obtained by calculation based on the optimal charge amount of the blast hole and the initial blast source depth of the area to be blasted. Based on the explosion performance index and stability index, determine whether to adjust the number of topology nodes in the area to be blasted. If not, do not adjust the number of topology nodes and generate an initiation control scheme with the initial topology nodes. If yes, adjust the number of topology nodes to generate an initiation control scheme. The detonation control scheme is implemented to perform synchronous detonation control of the electronic detonator.

[0015] Compared to existing solutions, the beneficial effects achieved by this invention are: On the one hand, the present invention determines the inter-hole delay and the optimal charge amount of the blast hole in the area to be blasted based on the prior data of blasting vibration. By analyzing historical blasting data and actual field conditions, it can accurately calculate the optimal detonation time interval and charge amount of each blast hole, so that the explosion energy can be more rationally distributed in time and space, optimize the blasting effect and improve the synchronous detonation accuracy of the electronic detonator. On the other hand, this invention collects geological parameter data of the area to be blasted and performs geological analysis and calculations to obtain a stability index for the area. Different geological conditions result in significant differences in the response to blasting, and the stability index can accurately reflect the geological characteristics of the area to be blasted. Adjusting blasting parameters based on the stability index, such as appropriately reducing the charge amount and adjusting the delay in geologically unstable areas, can effectively control the damage range of the blasting to the surrounding rock mass, ensuring that the blasting effect meets expectations while guaranteeing construction safety. Attached Figure Description

[0016] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments recorded in this invention. For those skilled in the art, other drawings can be obtained based on these drawings.

[0017] Figure 1 This is a system block diagram of a topology-optimized electronic detonator network synchronous initiation control system according to an embodiment of the present invention; Figure 2 This is a flowchart of a synchronous initiation control system for electronic detonator networks based on topology optimization, according to an embodiment of the present invention. Figure 3 This is a flowchart illustrating the process of a synchronous initiation control method for electronic detonator networks based on topology optimization, according to an embodiment of the present invention. Figure 4 This is a structural block diagram of an electronic device according to an embodiment of the present invention. Detailed Implementation

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

[0019] Furthermore, the described features, structures, or characteristics can be combined in any suitable manner in one or more exemplary embodiments. Numerous specific details are provided in the following description to give a full understanding of exemplary embodiments of this disclosure. However, those skilled in the art will recognize that the technical solutions of this disclosure can be practiced with one or more of the specific details omitted, or other methods, components, steps, etc., can be employed. In other instances, well-known structures, methods, implementations, or operations are not shown or described in detail to avoid obscuring various aspects of this disclosure.

[0020] This embodiment provides a topology-optimized electronic detonator network synchronous initiation control system. Figure 1 This is a system block diagram of a topology-optimized electronic detonator network synchronous initiation control system according to an embodiment of the present invention, as shown below. Figure 1 As shown, the system includes: The blasting information acquisition unit is used to treat each blast hole in the area to be blasted as a topology node, construct an electronic detonator network based on the electronic detonators of the topology nodes, and acquire information about the area to be blasted where the topology nodes of the electronic detonator network are located. The information about the area to be blasted includes geological parameter data and prior data of blasting vibration. The blasting parameter determination unit is used to determine the inter-hole delay and the optimal charge amount of the borehole in the area to be blasted based on the prior data of blasting vibration. The data calculation unit is used to perform geological analysis and calculation on geological parameter data to obtain the stability index of the area to be blasted, and to calculate the explosive performance index based on the optimal charge amount of the blast hole and the initial blast source depth of the area to be blasted. The detonation control scheme generation unit is used to determine whether to adjust the number of topological nodes in the area to be blasted based on the explosion performance index and stability index. If not, the number of topological nodes is not adjusted and the detonation control scheme is generated with the initial topological nodes. If yes, the number of topological nodes is adjusted to generate the detonation control scheme. The synchronous detonation execution unit is used to execute the detonation control scheme to perform synchronous detonation control of the electronic detonator.

[0021] In summary, this invention treats each borehole in the area to be blasted as a topological node, constructs an electronic detonator network based on these nodes, collects information about the area to be blasted where each topological node is located, and determines the inter-hole delay and optimal charge amount for each borehole based on prior blasting vibration data. Geological analysis is performed on the geological parameter data to obtain the stability index of the area to be blasted. An explosive performance index is calculated based on the optimal charge amount for each borehole and the initial blast source depth of the area to be blasted. Based on the explosive performance index and the stability index, it is determined whether the number of topological nodes in the area to be blasted needs to be adjusted; if not, the topology is not adjusted. The number of nodes is determined, and an initial topology node is used to generate an initiation control scheme. If necessary, the number of topology nodes is adjusted to generate an initiation control scheme. The initiation control scheme is executed to perform synchronous initiation control of the electronic detonators. This allows for precise calculation of the optimal initiation time interval and charge amount for each borehole, enabling a more rational distribution of explosive energy in time and space, optimizing blasting effects, and improving the accuracy of synchronous initiation of the electronic detonators. Blasting parameters are adjusted according to the stability index. For example, in geologically unstable areas, appropriately reducing the charge amount and adjusting the delay can effectively control the damage range of the blasting to the surrounding rock mass, ensuring that the blasting effect meets expectations while guaranteeing construction safety.

[0022] In some embodiments, determining the inter-hole delay and optimal charge amount for the borehole in the area to be blasted based on prior data of blasting vibration specifically includes the following process: Design a field cut area with the same environment as the area to be blasted, and conduct a single-hole single-free-face blasting experiment in the field cut area to collect the ground vibration curve of the single hole. The single hole in the field cut area corresponds one-to-one with the topological node of the area to be blasted. The ground vibration curves were fitted using Fourier series. By analyzing the combined vibration velocities under different inter-hole delays for each borehole charge, the optimal charge amount and inter-hole delay for each borehole were obtained. Step 1: Perform multi-hole vibration superposition calculation on the ground vibration curve of a single hole according to the different time delays between candidate holes to obtain all synthetic waveform functions. Based on all synthetic waveform functions, obtain the upper limit of the maximum synthetic vibration velocity value. and lower limit The method involves multi-hole vibration superposition calculation of the single-hole ground vibration curve based on different inter-hole delays, yielding all synthetic waveform functions. In blasting engineering, vibrations generated by multiple borehole detonations superimpose on the ground, forming complex vibration waveforms. To accurately predict vibration under different inter-hole delays, multi-hole vibration superposition calculation is needed for the single-hole ground vibration curve under different inter-hole delays. The basic principle is based on the linear superposition principle of vibration signals, meaning that the displacement, velocity, or acceleration values ​​of multiple vibration signals at the same moment can be simply added to obtain the value of the synthetic vibration signal at that moment.

[0023] Step 2, for any borehole charge Based on the limit value of the maximum synthesized vibration velocity obtained in step one This yields the set of maximum vibration velocity lower limits corresponding to all candidate charge quantities. Then, it identifies the values ​​in the set that are lower than the safe vibration velocity standard. The maximum value; Step 3: Select the set of maximum vibration velocity lower limit values ​​that are less than the safe vibration velocity standard. The maximum value of the selected charge amount is denoted as the optimal charge amount of the borehole, and the inter-hole delay corresponding to the optimal charge amount of the borehole is the inter-hole delay.

[0024] In some embodiments, the acquisition of geological parameter data specifically includes the following process: A three-dimensional coordinate system was established by selecting the seismic wave detector located at the center of the area to be blasted as the origin, and the position coordinates of each seismic wave detector were determined: , , , , , ,in, Given the distances of other seismic wave detectors to the central topology node, calculate the acoustic path from each seismic wave detector to the non-origin seismic wave detectors, and calculate the difference between each acoustic path. , 6. Calculation and Estimated time delay of sound waves According to the difference Delay estimate Determine the initial elevation angle of the target seismic wave detector installation location. and azimuth And based on the initial pitch angle and azimuth Determine the installation location coordinates of the target seismic wave detector, install the target seismic wave detector based on the installation location coordinates, and collect geological parameter data based on the target seismic wave detector.

[0025] Among them, the cross-correlation calculation formula between sound source signals is used to calculate... and Inter-delay estimate : ; in, This represents the expression for the delayed estimate. For Weighting function It is the cross-power spectrum of the signal. These are parameters for frequency domain analysis.

[0026] Among them, based on the difference Delay estimate Determine the initial pitch angle of the target sound source and azimuth Specifically, the process includes the following: Assuming the target seismic wave detector is installed at the following location: Calculate the distance between the target seismic wave detector and the origin seismic wave detector. : ; in, Indicates the speed of sound; A set of equations is established based on the locations of each seismic wave detector: ; Based on spherical coordinate expressions, equations and The expression yields the initial pitch angle. and azimuth : ; .

[0027] In some embodiments, Figure 2 This is a flowchart of a topology-optimized electronic detonator network synchronous initiation control system according to an embodiment of the present invention, as shown below. Figure 2 As shown, the process of obtaining the stability index of the area to be blasted through geological parameter data through geological analysis and calculation includes the following steps: Step S201: Obtain the root mean square amplitude, average amplitude, and average instantaneous frequency of seismic waves in the area to be blasted based on geological parameter data; Among them, the root mean square amplitude is used to identify amplitude anomalies or describe sequence stratigraphy and track seismic anomalies in stratigraphy; the mean amplitude is used to identify amplitude anomalies or sequence stratigraphy features and identify lithology and lithological changes; and the mean instantaneous frequency is used to detect amplitude absorption anomalies and track changes in related frequency absorption features caused by faults, lithology or stratigraphy. Step S202: Based on the root mean square amplitude, average amplitude, and average instantaneous frequency, look up the corresponding correlation coefficients in the correlation coefficient table to obtain the root mean square amplitude correlation coefficient. Correlation coefficient of average amplitude Correlation coefficient with average instantaneous frequency ; The correlation coefficient table stores the correlation coefficients corresponding to the value ranges of the root mean square amplitude, average amplitude, and average instantaneous frequency; this correlation coefficient table is derived from historical blasting data.

[0028] Step S203: Calculate the root mean square amplitude correlation coefficient. Correlation coefficient of average amplitude Correlation coefficient with average instantaneous frequency After normalization, the stability index of the area to be blasted is obtained by inputting it into the correlation formula. .

[0029] The correlation formula is as follows: ; in, , and This is the weight, and its value is set by the system; the selectable values ​​are... 0.2, 0.3, : 0.5.

[0030] In some embodiments, the calculation of the explosive performance index based on the optimal charge amount of the borehole and the initial blast source depth of the area to be blasted specifically includes the following process: Determine the chemical reaction kinetics and specific thermal decomposition kinetic parameters of electronic detonator explosives: reaction order n, activation energy. Pre-exponential frequency factor and the heat of reaction Q per unit mass; Obtain the optimal charge M for the borehole, the number of topological nodes K, and the initial blast source depth of the area to be blasted. ; Obtaining the universal molar gas constant The heat released per unit mass of reactant Thermal conductivity λ, density of reactants d, the radius of the bottom surface of the reactant; c, the specific heat capacity of the reactant at constant volume. After normalization, the thermal conductivity parameter W was calculated as follows: ; After normalization, the critical explosion temperature was calculated based on the thermal conductivity parameter W. : ; After normalization, based on the critical explosion temperature The thermal explosion performance index was calculated. : .

[0031] In some embodiments, determining whether to adjust the initial blast source depth based on the explosion performance index and stability index specifically includes the following process: Determine whether the explosive performance index is greater than the stability index. If so, the initial blast source depth is deemed to meet the blasting conditions, and no adjustment is made to the initial blast source depth. If not, the initial blast source depth is deemed not to meet the blasting conditions, and adjustment is required.

[0032] Furthermore, generating a detonation control scheme without adjusting the number of topology nodes and using the initial topology nodes specifically includes: generating a detonation control scheme based on the initially set number of topology nodes, initial blast source depth, inter-hole delay of topology nodes, and detonation sequence of topology nodes. Generating a detonation control scheme by adjusting the number of topology nodes specifically includes: increasing the number of topology nodes until the explosive performance index is greater than the stability index, and generating a detonation control scheme based on the adjusted number of topology nodes, initial blast source depth, inter-hole delay of topology nodes, and detonation sequence of topology nodes.

[0033] In some embodiments, the present invention also provides a method for synchronous initiation control of electronic detonator networks based on topology optimization. Figure 3 This is a flowchart illustrating a synchronous initiation control method for electronic detonator networks based on topology optimization, according to an embodiment of the present invention. Figure 3 As shown, the method includes the following steps: Step S301: Treat each blast hole in the area to be blasted as a topology node, construct an electronic detonator network based on the electronic detonators of the topology nodes, and collect information on the area to be blasted where the topology nodes of the electronic detonator network are located. The information on the area to be blasted includes geological parameter data and prior data of blasting vibration. Step S302: Determine the inter-hole delay and optimal charge amount of the borehole in the area to be blasted based on the prior data of blasting vibration; Step S303: Perform geological analysis and calculation on the geological parameter data to obtain the stability index of the area to be blasted, and calculate the explosive performance index based on the optimal charge amount of the blast hole and the initial blast source depth of the area to be blasted. Step S304: Based on the explosion performance index and stability index, determine whether to adjust the number of topology nodes in the area to be blasted. If not, do not adjust the number of topology nodes and generate an initiation control scheme with the initial topology nodes. If yes, adjust the number of topology nodes to generate an initiation control scheme. Step S305: Execute the detonation control scheme to perform synchronous detonation control of the electronic detonator.

[0034] This embodiment also provides an electronic device. Figure 4 This is a structural block diagram of an electronic device according to an embodiment of the present invention, such as... Figure 4 As shown, the electronic device includes a memory 401 and a processor 402. The memory 401 stores a computer program. When the computer program is executed by the processor 402, the processor 402 executes the electronic detonator network synchronous initiation control system and method based on topology optimization as described in any of the above embodiments.

[0035] The memory 401 may be an electronic memory such as flash memory, EEPROM (Electrically Erasable Programmable Read-Only Memory), EPROM, hard disk, or ROM. The memory 401 has storage space 403 for program code 413 for performing any of the method steps described above. For example, the storage space 403 for program code may include individual program codes 413 for implementing the various steps in the methods described above. This program code can be read from or written to one or more computer program products. These computer program products include program code carriers such as hard disks, CDs, memory cards, or floppy disks. The program code may be compressed, for example, in a suitable form. When run by a computing processing device, this code causes the computing processing device to perform the various steps in the methods described above. This program code can be read from or written to one or more computer program products. These computer program products include program code carriers such as hard disks, CDs, memory cards, or floppy disks. The program code may be compressed, for example, in a suitable form. When this code is run by a computing device, it causes the computing device to perform the various steps in the topology-optimized electronic detonator network synchronous initiation control system described above.

[0036] The above embodiments can be implemented, in whole or in part, by software, hardware, firmware, or any other combination thereof. When implemented using software, the above embodiments can be implemented, in whole or in part, as a computer program product. The computer program product includes one or more computer instructions or computer programs. When the computer instructions or computer programs are loaded or executed on a computer, all or part of the processes or functions described in the embodiments of this application are generated. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions can be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, the computer instructions can be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via wired (e.g., infrared, wireless, microwave, etc.) means. The computer-readable storage medium can be any available medium that a computer can access or a data storage device such as a server or data center that includes one or more sets of available media. The available medium can be a magnetic medium (e.g., floppy disk, hard disk, magnetic tape), an optical medium (e.g., DVD), or a semiconductor medium. A semiconductor medium can be a solid-state drive.

[0037] Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0038] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0039] In the several embodiments provided in this application, it should be understood that the disclosed systems, apparatuses, and methods can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative; for instance, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some interfaces; the indirect coupling or communication connection between apparatuses or units may be electrical, mechanical, or other forms.

[0040] The units described as separate components may or may not be physically separate. 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 units can be selected to achieve the purpose of this embodiment according to actual needs.

[0041] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A synchronous initiation control system for electronic detonator networks based on topology optimization, characterized in that, The system includes: The blasting information acquisition unit is used to treat each blast hole in the area to be blasted as a topology node, construct an electronic detonator network based on the electronic detonators of the topology nodes, and acquire information about the area to be blasted where the topology nodes of the electronic detonator network are located. The information about the area to be blasted includes geological parameter data and prior data of blasting vibration. The blasting parameter determination unit is used to determine the inter-hole delay and the optimal charge amount of the borehole in the area to be blasted based on the prior data of blasting vibration. The data calculation unit is used to perform geological analysis and calculation on geological parameter data to obtain the stability index of the area to be blasted, and to calculate the explosive performance index based on the optimal charge amount of the blast hole and the initial blast source depth of the area to be blasted. The detonation control scheme generation unit is used to determine whether to adjust the number of topological nodes in the area to be blasted based on the explosion performance index and stability index. If not, the number of topological nodes is not adjusted and the detonation control scheme is generated with the initial topological nodes. If yes, the number of topological nodes is adjusted to generate the detonation control scheme. The synchronous detonation execution unit is used to execute the detonation control scheme to perform synchronous detonation control of the electronic detonator.

2. The electronic detonator network synchronous initiation control system based on topology optimization according to claim 1, characterized in that, Determining the inter-hole delay and optimal charge amount for the blasting area based on prior data of blasting vibration specifically includes the following process: A cut-out area with the same environment as the area to be blasted was designed, and a single-hole, single-free-face blasting experiment was conducted within the cut-out area. Ground vibration curves for each hole were collected, with a one-to-one correspondence between the topological nodes of the cut-out area and the area to be blasted. The ground vibration curves were fitted using a Fourier series form. By analyzing the combined vibration velocity under different inter-hole delays for each borehole charge, the optimal charge amount and inter-hole delay for each borehole were obtained. Step 1: Perform multi-hole vibration superposition calculation on the ground vibration curve of a single hole according to the different time delays between candidate holes to obtain all synthetic waveform functions. Based on all synthetic waveform functions, obtain the upper limit of the maximum synthetic vibration velocity value. and lower limit ; Step 2, for any borehole charge Based on the limit value of the maximum synthesized vibration velocity obtained in step one This yields the set of maximum vibration velocity lower limits corresponding to all candidate charge quantities. Then, it identifies the values ​​in the set that are lower than the safe vibration velocity standard. The maximum value; Step 3: Select the set of maximum vibration velocity lower limit values ​​that are less than the safe vibration velocity standard. The maximum value of the selected charge amount is denoted as the optimal charge amount of the borehole, and the inter-hole delay corresponding to the optimal charge amount of the borehole is the inter-hole delay.

3. The electronic detonator network synchronous initiation control system based on topology optimization according to claim 1, characterized in that, Specific details of collecting geological parameter data The process includes the following: A three-dimensional coordinate system was established by selecting the seismic wave detector located at the center of the area to be blasted as the origin, and the position coordinates of each seismic wave detector were determined: , , , , , ,in, Given the distances of other seismic wave detectors to the central topology node, calculate the acoustic path from each seismic wave detector to the non-origin seismic wave detectors, and calculate the difference between each acoustic path. , 6. Calculation and Estimated time delay of sound waves According to the difference Delay estimate Determine the initial elevation angle of the target seismic wave detector installation location. and azimuth And based on the initial pitch angle and azimuth Determine the installation location coordinates of the target seismic wave detector, install the target seismic wave detector based on the installation location coordinates, and collect geological parameter data based on the target seismic wave detector.

4. The electronic detonator network synchronous initiation control system based on topology optimization according to claim 1, characterized in that, The process of obtaining the stability index of the area to be blasted through geological parameter data through geological analysis and calculation includes the following steps: Based on geological parameter data, the root mean square amplitude, average amplitude, and average instantaneous frequency of seismic waves in the area to be blasted are obtained. Among them, the root mean square amplitude is used to identify amplitude anomalies or describe sequence stratigraphy and track seismic anomalies in the stratigraphy; the average amplitude is used to identify amplitude anomalies or sequence stratigraphy characteristics and identify lithology and lithological changes; the average instantaneous frequency is used to detect amplitude absorption anomalies and track changes in related frequency absorption characteristics caused by faults, lithology, or stratigraphy. Based on the root mean square amplitude, average amplitude, and average instantaneous frequency, the corresponding correlation coefficients are looked up in the correlation coefficient table to obtain the root mean square amplitude correlation coefficient. Correlation coefficient of average amplitude Correlation coefficient with average instantaneous frequency The correlation coefficient table stores the correlation coefficients corresponding to the ranges of the root mean square amplitude, average amplitude, and average instantaneous frequency. The root mean square amplitude correlation coefficient Correlation coefficient of average amplitude Correlation coefficient with average instantaneous frequency After normalization, the stability index of the area to be blasted is obtained by inputting it into the correlation formula. The correlation formula is as follows: in, , and This is the weight, and its value is set by the system.

5. The electronic detonator network synchronous initiation control system based on topology optimization according to claim 1, characterized in that, The calculation of the explosive performance index based on the optimal charge amount of the borehole and the initial blast source depth of the area to be blasted includes the following process: Determine the chemical reaction kinetics and specific thermal decomposition kinetic parameters of electronic detonator explosives: reaction order n, activation energy. Pre-exponential frequency factor and the heat of reaction Q per unit mass; Obtain the optimal charge M for the borehole, the number of topological nodes K, and the initial blast source depth of the area to be blasted. ; Obtaining the universal molar gas constant The heat released per unit mass of reactant Thermal conductivity λ, density of reactants d, the radius of the bottom surface of the reactant; c, the specific heat capacity of the reactant at constant volume. After normalization, the thermal conductivity parameter W was calculated as follows: ; After normalization, the critical explosion temperature was calculated based on the thermal conductivity parameter W. : ； After normalization, based on the critical explosion temperature The thermal explosion performance index was calculated. : 。 6. The electronic detonator network synchronous initiation control system based on topology optimization according to claim 1, characterized in that, The decision to adjust the initial blast depth is based on the explosive performance index and stability index. The process includes the following: Determine whether the explosive performance index is greater than the stability index. If so, the initial blast source depth is deemed to meet the blasting conditions, and no adjustment is made to the initial blast source depth. If not, the initial blast source depth is deemed not to meet the blasting conditions, and adjustment is required.

7. The electronic detonator network synchronous initiation control system based on topology optimization according to claim 1, characterized in that, Generating a detonation control scheme without adjusting the number of topology nodes and using the initial topology nodes specifically includes the following process: A detonation control scheme is generated based on the initial set number of topology nodes, initial detonation source depth, inter-hole delay of topology nodes, and detonation sequence of topology nodes.

8. The electronic detonator network synchronous initiation control system based on topology optimization according to claim 1, characterized in that, Adjusting the number of topology nodes to generate a detonation control scheme specifically includes the following processes: Increase the number of topology nodes until the explosion performance index is greater than the stability index, and generate an initiation control scheme based on the adjusted number of topology nodes, initial blast source depth, inter-hole delay of topology nodes, and initiation sequence of topology nodes.

9. A method for synchronous initiation control of electronic detonator networks based on topology optimization, characterized in that, The method applicable to the topology-optimized electronic detonator network synchronous initiation control system according to any one of claims 1 to 8, the method comprising: Each blast hole in the area to be blasted is taken as a topological node. An electronic detonator network is constructed based on the electronic detonators of the topological nodes. Information on the area to be blasted where the topological nodes of the electronic detonator network are located is collected. The information on the area to be blasted includes geological parameter data and prior data of blasting vibration. The inter-hole delay and optimal charge amount of the blast hole in the area to be blasted are determined based on prior data of blasting vibration. The stability index of the area to be blasted is obtained by geological analysis and calculation of geological parameter data, and the explosive performance index is obtained by calculation based on the optimal charge amount of the blast hole and the initial blast source depth of the area to be blasted. Based on the explosion performance index and stability index, determine whether to adjust the number of topology nodes in the area to be blasted. If not, do not adjust the number of topology nodes and generate an initiation control scheme with the initial topology nodes. If yes, adjust the number of topology nodes to generate an initiation control scheme. The detonation control scheme is implemented to perform synchronous detonation control of the electronic detonator.

Images