Prediction and early warning method and system for tunnel blasting vibration

Abstract

The invention provides a prediction and early warning method and system for tunnel blasting vibration. The method comprises the following steps: acquiring a basic prediction model and actual explosive parameters, actual blasting parameters and actual geological parameters of each target blasting area; respectively calculating an explosive parameter correction item, a blasting parameter correction item and a geological parameter correction item of each target blasting area by utilizing the actual explosive parameters, the actual blasting parameters and the actual geological parameters of each target blasting area; according to the basic prediction model and the explosive parameter correction item, the blasting parameter correction item and the geological parameter correction item of each target blasting area, respectively calculating a vibration prediction value of each target blasting area; calculating a comprehensive vibration predicted value through the vibration predicted value of each target blasting area; and carrying out vibration early warning by using the comprehensive vibration predicted value. According to the method, for a collaborative blasting scene, the comprehensive vibration prediction value is calculated through the vibration prediction value of each target blasting area, and the vibration condition can be reflected more accurately.

Description

Technical Field

[0001] The present invention relates to the technical field of tunnel blasting, and particularly to a method and system for predicting and warning tunnel blasting vibration. Background Art

[0002] With the continuous growth of the economy, the demand for infrastructure such as railways, highways, and subways has also been greatly developed. During the construction of these infrastructures, tunnel blasting is often required. During tunnel blasting, a target blasting area can be selected in the rock mass of the tunnel chamber or the surrounding rock of the tunnel, and one or more blast holes are set in the target blasting area, and then explosives are filled in the blast holes, and then the explosives are detonated to carry out the blasting.

[0003] It should be noted that in some long tunnel or large tunnel projects, in order to improve the construction efficiency, the tunnel may also be divided into multiple sections along the length direction or the cross-section direction, and then the target blasting areas are respectively selected in each section for simultaneous blasting, that is, multi-section collaborative blasting. However, this kind of collaborative blasting is likely to cause strong vibration in the surrounding rock of the tunnel and the surrounding area of the tunnel, and may affect the stability of the surrounding rock of the tunnel and the buildings in the surrounding area. Therefore, during tunnel blasting, especially during multi-section collaborative blasting, it is often necessary to predict the vibration degree caused by tunnel blasting, and then be able to issue necessary warnings based on the prediction results of the vibration degree.

[0004] Patent CN108489601A discloses a method for monitoring and controlling blasting vibration when a tunnel passes closely under underground pipelines, which includes finding a formation similar to the formation where the pipeline is located at a certain distance from the pipeline to install sensors to monitor the vibration velocity generated when the seismic wave formed after tunnel blasting propagates to the sensors, so as to achieve vibration detection.

[0005] Patent CN118965982A discloses a method and system for predicting the impact of tunnel blasting construction on existing buildings, which includes obtaining a tunnel blasting data sample set and normalizing the tunnel blasting data sample set. The tunnel blasting data samples include tunnel blasting impact parameters and corresponding blasting vibration velocities; using an optimization strategy to improve the fish swarm search algorithm, and combining a backpropagation neural network and the improved fish swarm search algorithm to construct a tunnel blasting vibration velocity prediction model; training and validating the tunnel blasting vibration velocity prediction model through the tunnel blasting data sample set, and then obtaining the predicted blasting vibration velocity according to the tunnel blasting impact parameters of the tunnel to be blasted; based on the predicted blasting vibration velocity, completing the prediction of the impact of tunnel blasting construction on existing buildings.

[0006] Obviously, in the methods disclosed in the above-mentioned patents CN108489601A and CN118965982A, the vibration prediction mainly targets single-target blasting areas. However, for the scenario of multi-section collaborative blasting, the accuracy of vibration prediction by these methods is usually low. Summary of the Invention

[0007] The present invention provides a method and system for predicting and warning tunnel blasting vibrations to solve the problem of low accuracy of vibration prediction for collaborative blasting in the prior art.

[0008] On the one hand, the present invention provides a method for predicting and warning tunnel blasting vibrations. The tunnel is divided into multiple sections along the length direction or cross-section direction, and target blasting areas are respectively selected from at least two of the multiple sections for collaborative blasting. The method includes: Obtain the basic prediction model and the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area; Use the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area to calculate the explosive parameter correction term, blasting parameter correction term, and geological parameter correction term of each target blasting area respectively; Calculate the vibration prediction value of each target blasting area through the basic prediction model and the explosive parameter correction term, blasting parameter correction term, and geological parameter correction term of each target blasting area; Calculate the comprehensive vibration prediction value through the vibration prediction values of each target blasting area; Use the comprehensive vibration prediction value for vibration warning.

[0009] Preferably, the explosive parameter correction term of each target blasting area is calculated by the following method: ; Where: i is the number of the target blasting area; D i is the actual detonation velocity of the explosive in the i-th target blasting area; Q ei is the actual heat of detonation of the explosive in the i-th target blasting area; ρ i is the actual density of the explosive in the i-th target blasting area; D 0 is the detonation velocity of the standard explosive; Q e0 is the heat of detonation of the standard explosive; ρ 0 is the density of the standard explosive; β 1i is the explosive parameter correction term of the i-th target blasting area.

[0010] Preferably, the blasting parameter correction term of each target blasting area is calculated by the following method: ; where: i is the number of the target blasting area; m i is the initiation method correlation parameter of the i-th target blasting area; c i is the actual hole spacing between blast holes in the i-th target blasting area; b i is the actual row spacing between blast holes in the i-th target blasting area; L i is the actual hole depth of the blast holes in the i-th target blasting area; c 0 is the hole spacing of the standard blast hole; b 0 is the row spacing of the standard blast hole; L 0 is the hole depth of the standard blast hole; β 2i is the blasting parameter correction term of the i-th target blasting area.

[0011] Preferably, the geological parameter correction terms of each target blasting area are calculated by the following method: ; where: i is the number of the target blasting area; p i is the geological structure correction coefficient of the i-th target blasting area; E i is the actual elastic modulus of the rock and soil mass in the i-th target blasting area; u i is the actual Poisson's ratio of the rock and soil mass in the i-th target blasting area; ρ ri is the actual density of the rock and soil mass in the i-th target blasting area; E 0 is the elastic modulus of the standard rock and soil mass; u 0 is the Poisson's ratio of the standard rock and soil mass; ρ r0 is the density of the standard rock and soil mass; β 3i is the geological parameter correction term of the i-th target blasting area.

[0012] Preferably, the vibration prediction values of each target blasting area are calculated respectively through the basic prediction model and the explosive parameter correction terms, blasting parameter correction terms and geological parameter correction terms of each target blasting area, specifically including calculating the vibration prediction values of each target blasting area through the following formula: V final-i = V i × β 1i × β 2i ×β 3i ; where: i is the number of the target blasting area; V final-i is the vibration prediction value of the i-th target blasting area obtained by calculation; β 1i is the explosive parameter correction term of the i-th target blasting area; β 2i is the blasting parameter correction term of the i-th target blasting area; β 3i is the geological parameter correction term of the i-th target blasting area; Vi is the influence amount of the blasting vibration velocity of the i-th target blasting area on the vibration prediction location, calculated by the basic prediction model.

[0013] Preferably, the comprehensive vibration prediction value is calculated through the vibration prediction values of each target blasting area, specifically including calculating the comprehensive vibration prediction value through the following formula: ; where V total is the calculated comprehensive vibration prediction value; V final-i is the vibration prediction value of the i-th target blasting area; n is the total number of target blasting areas; ω i is the angular frequency of the vibration wave of the i-th target blasting area; is the initial phase of the vibration wave of the i-th target blasting area; t is the time interval between the moment of vibration detection and the moment of collaborative blasting.

[0014] Preferably, the comprehensive vibration prediction value is calculated through the vibration prediction values of each target blasting area, specifically including calculating the comprehensive vibration prediction value through the following formula: ; where V total’ is the calculated comprehensive vibration prediction value; V final-i is the vibration prediction value of the i-th target blasting area; n is the total number of target blasting areas; γ is a reduction coefficient, where the value of γ is greater than 0 and less than 1.

[0015] Preferably, the basic prediction model is specifically the Sadovsky formula.

[0016] Preferably, vibration early warning is performed using the comprehensive vibration prediction value, specifically including: judging whether the comprehensive vibration prediction value is greater than or equal to a threshold value; performing vibration early warning when the comprehensive vibration prediction value is greater than or equal to the threshold value.

[0017] On the other hand, the present invention provides a prediction and early warning system for tunnel blasting vibration. The tunnel is divided into multiple sections along the length direction or the cross-section direction, and target blasting areas are respectively selected from at least two of the multiple sections for collaborative blasting. The system includes: An acquisition unit for acquiring a basic prediction model and the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area; A first calculation unit, configured to calculate an explosive parameter correction term, a blasting parameter correction term, and a geological parameter correction term for each target blasting area respectively by using the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area; A second calculation unit, configured to calculate a vibration prediction value for each target blasting area respectively through the basic prediction model and the explosive parameter correction term, the blasting parameter correction term, and the geological parameter correction term of each target blasting area; A third calculation unit, configured to calculate a comprehensive vibration prediction value through the vibration prediction values of each target blasting area; An early warning unit, configured to perform vibration early warning by using the comprehensive vibration prediction value.

[0018] By using the tunnel blasting vibration prediction and early warning method provided in the embodiments of the present application, the tunnel is divided into multiple sections along the length direction or the cross-section direction, and target blasting areas are respectively selected from at least two sections among the multiple sections for collaborative blasting. Wherein, the method includes obtaining a basic prediction model and the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area, and then calculating an explosive parameter correction term, a blasting parameter correction term, and a geological parameter correction term for each target blasting area respectively by using the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area, and then calculating a vibration prediction value for each target blasting area respectively through the basic prediction model and the explosive parameter correction term, the blasting parameter correction term, and the geological parameter correction term of each target blasting area, and then calculating a comprehensive vibration prediction value through the vibration prediction values of each target blasting area, and then performing vibration early warning by using the comprehensive vibration prediction value. For the scenario of collaborative blasting of multiple target blasting areas during the tunnel construction process, the comprehensive vibration prediction value calculated through the vibration prediction values of each target blasting area can more accurately reflect the vibration condition of the vibration prediction location, so the problems in the prior art can be solved. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following will briefly introduce the drawings required for the description of the embodiments or the prior art. Obviously, the drawings in the following description are some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0020] Figure 1 It is a schematic flowchart of a tunnel blasting vibration prediction and early warning method provided by the present invention; Figure 2 It is a structural block diagram of a tunnel blasting vibration prediction and early warning system provided by the present invention; Figure 3It is a schematic structural diagram of the electronic device provided by the present invention. Specific Embodiments

[0021] To make the objectives, technical solutions and advantages of the present invention clearer, the technical solutions in the present invention will be clearly and completely described below with reference to the accompanying drawings in the present invention. Apparently, the described embodiments are some but not all of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments in the present invention without making creative efforts shall fall within the protection scope of the present invention.

[0022] As mentioned above, in some long tunnel or large tunnel projects, the construction efficiency of the tunnel can often be improved through cooperative blasting. However, cooperative blasting is likely to cause strong vibrations in the tunnel surrounding rock and the surrounding area of the tunnel, thereby affecting the stability of the tunnel surrounding rock and the buildings in the surrounding area of the tunnel. Therefore, during the cooperative blasting process, it is necessary to predict and warn of the vibration level caused by the tunnel blasting. Currently, the vibration prediction methods for tunnel blasting mainly target the vibration prediction during blasting in a single target blasting area. However, for the scenario of multi-section cooperative blasting, the accuracy of vibration prediction by these methods is usually low.

[0023] In view of this, the embodiments of the present application provide a method and system for predicting and warning of tunnel blasting vibration, which can be used to perform more accurate vibration prediction for the scenario of multi-section cooperative blasting. For the convenience of understanding, the embodiments of the present application can be described as a whole here. For example, the technical solutions provided by the embodiments of the present application can be applied to software products. For example, the method provided by the embodiments of the present application can be designed as an application program (i.e., software product), and then installed on an electronic device such as a user terminal or a server to apply the method provided by the embodiments of the present application.

[0024] For example, the software product can be installed on the electronic devices of the user terminal, including mobile phones, tablets, computers and other electronic devices, and then the method provided by the embodiments of the present application can be executed; the software product can also be installed on the electronic devices of the server side, including servers or server clusters, and then the method provided by the embodiments of the present application can be executed.

[0025] In practical applications, the method can also be applied as a hardware product. For example, a dedicated hardware device can be designed to implement the method provided by the embodiments of the present application.

[0026] Here, a holistic description of the technical solution of the method provided in the embodiments of the present application can be given first. This method is for the scenario of collaborative blasting in tunnels, and is used for vibration prediction and early warning of collaborative blasting. Specifically, the tunnel can be divided into multiple sections along the length direction or the cross-sectional direction, and target blasting areas are respectively selected from at least two of these multiple sections for collaborative blasting. For example, the tunnel can be divided into multiple sections along the length direction, and then a section is respectively selected from near the exit and the entrance of the tunnel (that is, both ends of the tunnel), and target blasting areas are respectively selected from these two sections for collaborative blasting to improve the engineering construction efficiency of the tunnel. However, this kind of collaborative blasting is obviously likely to cause strong vibrations. Therefore, the method provided in the embodiments of the present application can be used for vibration prediction and early warning.

[0027] Figure 1 This is a schematic flowchart of a method for predicting and early warning tunnel blasting vibrations provided by the embodiments of the present invention. This method can be executed by an electronic device at the user end or by an electronic device at the server end. Here, taking the server (that is, the electronic device at the server end) executing this method as an example, this method will be described. As Figure 1 shown, this method may include: Step S11: Obtain a basic prediction model and the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area.

[0028] In practical applications, for the convenience of subsequent description, each target blasting area can be numbered here. For example, if n target blasting areas are used for collaborative blasting, each target blasting area can be numbered 1, 2, 3,..., n in sequence. After that, i can be used to represent the number of each target blasting area, that is, i can be any value from 1 to n.

[0029] Among them, the actual explosive parameters specifically include the actual detonation velocity D of the explosive, the actual heat of detonation Q e and the actual density ρ of the explosive. These three key parameters can have a greater impact on the vibration situation during explosive detonation. Therefore, among the actual explosive parameters of each target blasting area obtained in this step S11, the actual detonation velocity, the actual heat of detonation, and the actual density of the explosive of each target blasting area can be included. For example, for the i-th target blasting area, its actual explosive parameters can specifically include the actual detonation velocity D of the explosive in the i-th target blasting area i , the actual heat of detonation Q ei and the actual density ρ of the explosive i . In practical applications, these three key parameters can be obtained through actual detection of the explosive.

[0030] In practical applications, the initiation method, as well as the actual hole depth of the blast holes in the target blasting area, the actual hole spacing between the blast holes, and the actual row spacing, will all affect the vibration level during the collaborative blasting process. Therefore, the actual blasting parameters in step S11 may include the actual hole spacing between the blast holes, the actual row spacing, and the actual hole depth in the target blasting area. For example, for the i-th target blasting area, its actual blasting parameters may include the actual hole spacing c i between the blast holes, the actual row spacing b i and the actual hole depth L i in the i-th target blasting area, as well as the initiation method for the i-th target blasting area. The initiation method for the i-th target blasting area is usually represented in a coded manner. For example, 1 indicates that the initiation method for the i-th target blasting area is hole-by-hole initiation, and 2 indicates that the initiation method for the i-th target blasting area is hole-by-hole initiation with row-by-row initiation.

[0031] In addition, the geological structure, the actual elastic modulus of the rock and soil mass, the actual Poisson's ratio, and the actual density in the geological parameters will also affect the vibration level during the collaborative blasting process. Therefore, the actual geological parameters in step S11 may include the geological structure, the actual elastic modulus, the actual Poisson's ratio, and the actual density of each target blasting area. For example, for the i-th target blasting area, the actual geological parameters of the i-th target blasting area may include the actual elastic modulus E i of the rock and soil mass, the actual Poisson's ratio u i and the actual density ρ ri in the i-th target blasting area, as well as the geological structure of the i-th target blasting area.

[0032] In practical applications, a rock and soil mass sample in the i-th target blasting area can be collected, and then the rock and soil mass sample can be tested to obtain the actual elastic modulus E i the actual Poisson's ratio u i and the actual density ρ ri . And the geological structure of the i-th target blasting area can be obtained by querying relevant geological data or conducting on-site surveys.

[0033] It should be noted that in the field of engineering blasting, the Sadovskii formula is usually used to predict the blasting vibration velocity (i.e., for blasting vibration prediction). However, considering that this application is for the scenario of collaborative blasting, the Sadovskii formula is difficult to directly apply to accurately predict the vibration velocity. But considering the wide application of the Sadovskii formula and the certain accuracy of the prediction results, in the scenario of this application, the Sadovskii formula is used as the basic prediction model.

[0034] Among them, the Sadovskii formula is V = K × (Q 1 / 3 / R)α In the Sadovskii formula, V is the predicted blasting vibration velocity; Q is the total charge of each blast hole in the blasted area (the target blasting area in this application); R is the distance between the blasted area and the vibration prediction location. In this application, in the case of collaborative blasting, the vibration at the vibration prediction location is predicted; α is the vibration attenuation coefficient related to the vibration transmission medium; K is the coefficient related to the vibration transmission medium and blasting conditions. Correspondingly, in this application, for the i-th target blasting area, its basic prediction model can be expressed as V i =K × (Q i 1 / 3 / R i ) α , where Q i is the total charge of each blast hole in the i-th target blasting area; R i is the distance between the i-th target blasting area and the vibration prediction location; V i is the influence amount of the blasting vibration velocity of the i-th target blasting area on the vibration prediction location calculated by this basic prediction model.

[0035] In practical applications, the vibration attenuation coefficient α and the coefficient K can be determined in advance. Specifically, the geological type of the tunnel can be determined first, and then the vibration attenuation coefficient α and the coefficient K can be obtained by querying relevant industry specifications. For example, for a geological type similar to granite geology, by querying industry specifications, the K value can be obtained in the range of 150 - 200, and the α value can be obtained in the range of 1.5 - 1.8; of course, in practical applications, other methods can also be used to determine the vibration attenuation coefficient α and the coefficient K in advance, and this is not limited here.

[0036] It should be emphasized that during the process of collaborative blasting, due to the differences in the positions of each target blasting area, different blasting strategies can be adopted, resulting in certain differences in their respective actual explosive parameters, actual blasting parameters, and actual geological parameters.

[0037] For the specific implementation method of this step S11, the user can input the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area through a user operation interface on a user-side electronic device such as a mobile phone or a computer. In this way, the user-side electronic device can send the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area to the server, so that the server can obtain the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area. Of course, the server can obtain this basic prediction model, for example, directly from the model library.

[0038] Step S12: Calculate the explosive parameter correction term, blasting parameter correction term, and geological parameter correction term for each target blasting area respectively by using the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area.

[0039] After obtaining the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area through the above Step S11, in this Step S12, the explosive parameter correction term, blasting parameter correction term, and geological parameter correction term for each target blasting area can be calculated respectively by using the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area. Specifically, for example, for the explosive parameter correction term of each target blasting area, substitute the actual explosive parameters of this target blasting area into Formula One to calculate the explosive parameter correction term of each target blasting area: Formula One In this Formula One, i is the number of the target blasting area; β 1i is the explosive parameter correction term of the i-th target blasting area; D i is the actual detonation velocity of the explosive in the i-th target blasting area; Q ei is the actual heat of detonation of the explosive in the i-th target blasting area; ρ i is the actual density of the explosive in the i-th target blasting area; D 0 is the detonation velocity of the standard explosive; Q e0 is the heat of detonation of the standard explosive; ρ 0 is the density of the standard explosive. Among them, regarding how to select the standard explosive, relevant industry manuals can be consulted, such as the "Blasting Engineer's Handbook" or the "Blasting Safety Regulations", etc., so as to determine the standard explosive according to the relevant requirements in the industry manual, and then obtain the detonation velocity D 0 , heat of detonation Q e0 and density ρ 0 of this standard explosive from the relevant database.

[0040] Therefore, after obtaining the actual explosive parameters of the i-th target blasting area through the above Step S11, that is, the actual detonation velocity D i , actual heat of detonation Q ei and actual density ρ i of the explosive, the actual detonation velocity D i , actual heat of detonation Q ei and actual density ρ i of the explosive in the i-th target blasting area can be substituted into this Formula One, so as to calculate and obtain the explosive parameter correction term β 1i of the i-th target blasting area.

[0041] For the blasting parameter correction terms for each target blasting area, the blasting parameter correction terms for each target blasting area can be calculated through Formula 2: Formula 2 In this Formula 2: i is the number of the target blasting area; m i is the initiation method correlation parameter for the i-th target blasting area, and the value of this m i is related to the initiation method of the i-th target blasting area. For example, if the initiation method of the i-th target blasting area is hole-by-hole initiation, the value of m i can be taken between 0.8 and 1.2. If the initiation method of the i-th target blasting area is row-by-row initiation, the value of m i can be taken between 0.9 and 1.1. Therefore, the value of this m i can be obtained through the initiation method of the i-th target blasting area; c i is the actual hole spacing between blast holes in the i-th target blasting area; b i is the actual row spacing between blast holes in the i-th target blasting area; L i is the actual hole depth of the blast holes in the i-th target blasting area; c 0 is the hole spacing of the standard blast holes; b 0 is the row spacing of the standard blast holes; L 0 is the hole depth of the standard blast holes; β 2i is the blasting parameter correction term for the i-th target blasting area.

[0042] Among them, for the hole spacing c 0 , row spacing b 0 , and hole depth L 0 of the standard blast holes, they can be determined according to the design specifications in blasting engineering, such as the "Safety Regulations for Blasting". This "Safety Regulations for Blasting" gives the recommended ranges of hole spacing c 0 , row spacing b 0 , hole depth L 0 and other blast hole parameters for different types of blasting projects. Appropriate values can be selected from this recommended range in combination with actual construction experience.

[0043] In actual applications like this, after obtaining the actual blasting parameters of the i-th target blasting area through the above step S11, that is, the actual hole spacing c i , actual row spacing b i , actual hole depth L i and the initiation method, the corresponding initiation method correlation parameter m i can be obtained first through the initiation method, and then the actual hole spacing c i , actual row spacing b i , actual hole depth Li and the detonation method-related parameter m i , substitute it into the second formula, so as to calculate the blasting parameter correction term β of the i-th target blasting area 2i .

[0044] For the geological parameter correction terms of each target blasting area, the geological parameter correction terms of each target blasting area can be calculated through the third formula: The third formula In the third formula: i is the number of the target blasting area; p i is the geological structure correction coefficient of the i-th target blasting area. Specifically, the geological structure of the i-th target blasting area can be obtained, and then the geological structure correction coefficient p can be determined according to the geological structure i , for example, the geological structure correction coefficient p can be determined by combining the development degree and spacing of the structural planes in the geological structure of the i-th target blasting area with engineering experience i , generally speaking, the lower the development degree of the structural plane, the closer the value of p i is to 1. When the spacing of the structural planes is relatively dense, the value of p i can be between 0.6 and 0.8; E i is the actual elastic modulus of the rock and soil mass in the i-th target blasting area; u i is the actual Poisson's ratio of the rock and soil mass in the i-th target blasting area; ρ ri is the actual density of the rock and soil mass in the i-th target blasting area; E 0 is the elastic modulus of the standard rock and soil mass; u 0 is the Poisson's ratio of the standard rock and soil mass; ρ r0 is the density of the standard rock and soil mass; β 3i is the geological parameter correction term of the i-th target blasting area.

[0045] For the elastic modulus E 0 , Poisson's ratio u 0 and density ρ r0 of the standard rock and soil mass in the third formula, they are obtained by querying and calculating the "Code for Geotechnical Investigation". Among them, for various different types of rock and soil, the "Code for Geotechnical Investigation" gives the empirical value ranges of mechanical parameters. Through the mechanical parameters of these various different types of rock and soil, the elastic modulus E 0 , Poisson's ratio u 0 and density ρ r0 can be estimated. Of course, the elastic modulus E 0 , Poisson's ratio u 0 and density ρ r0 can also be obtained by actual measurement.

[0046] Therefore, after obtaining the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area through the above-mentioned step S11, in this step S12, the explosive parameter correction term, blasting parameter correction term, and geological parameter correction term of each target blasting area can be calculated respectively by using the above-mentioned formulas 1-3.

[0047] Step S13: Calculate the vibration prediction values of each target blasting area respectively through the basic prediction model and the explosive parameter correction term, blasting parameter correction term, and geological parameter correction term of each target blasting area.

[0048] After obtaining the explosive parameter correction term, blasting parameter correction term, and geological parameter correction term of each target blasting area through the above-mentioned step S12, the vibration prediction values of each target blasting area can be further calculated through this step S13. Specifically, the vibration prediction values of each target blasting area can be calculated respectively by using the following formula 4: V final-i = V i × β 1i × β 2i ×β 3i Formula 4 In this formula 4, i is the number of the target blasting area; V final-i is the vibration prediction value of the i-th target blasting area calculated; β 1i is the explosive parameter correction term of the i-th target blasting area; β 2i is the blasting parameter correction term of the i-th target blasting area; β 3i is the geological parameter correction term of the i-th target blasting area; Vi is the influence amount of the blasting vibration velocity of the i-th target blasting area on the vibration prediction location calculated through the basic prediction model. Obviously, the vibration prediction value V final-i of the i-th target blasting area reflects the influence amount of the blasting vibration velocity of the i-th target blasting area on the vibration prediction location.

[0049] Among them, the values of β 1i , β 2i and β 3i can be calculated respectively by using the above-mentioned formulas 1-3. In the above-mentioned step S11, it is mentioned that V i =K × (Q i 1 / 3 / R i ) α , so the total charge Q i of each blast hole in the i-th target blasting area, and the distance R i between the i-th target blasting area and the vibration prediction location can be substituted into the formula Vi =K × (Q i 1 / 3 / R i ) α so as to calculate the value of this V i , and then substitute the values of β 1i , β 2i and β 3i as well as the value of V i into Formula Four, and the vibration prediction value V final-i of the i-th target blasting area can be calculated.

[0050] Step S14: Calculate the comprehensive vibration prediction value based on the vibration prediction values of each target blasting area.

[0051] After calculating the vibration prediction values of each target blasting area through the above Step S13, since the vibration prediction values of each target blasting area can only reflect the influence amount of the corresponding target blasting area on the blasting vibration velocity at the vibration prediction location. For example, the vibration prediction value V final-i of the i-th target blasting area mentioned above reflects the influence amount of the i-th target blasting area on the blasting vibration velocity at the vibration prediction location. Therefore, in this Step S14, it is necessary to further calculate the comprehensive vibration prediction value based on the vibration prediction values of each target blasting area. Among them, the comprehensive vibration prediction value reflects the total influence amount of the coordinated blasting of each target blasting area on the blasting vibration velocity at the vibration prediction location.

[0052] In practical applications, the specific implementation method for this Step S14 can be to use the following Formula Five to calculate the comprehensive vibration prediction value: Formula Five In this Formula Five, V total is the calculated comprehensive vibration prediction value; V final-i is the vibration prediction value of the i-th target blasting area; n is the total number of target blasting areas; ω i is the angular frequency of the vibration wave of the i-th target blasting area; is the initial phase of the vibration wave of the i-th target blasting area; t is the time interval between the moment of vibration detection and the moment of coordinated blasting.

[0053] In this Formula Five, the vibration prediction values V final-i of each target blasting area are vectorially superposed to calculate the comprehensive vibration prediction value V total . The comprehensive vibration prediction value V total obtained by this calculation methodThe accuracy is relatively high. However, in Formula Five, it is necessary to determine the angular frequency ω of the vibration wave in each target blasting area, as well as the initial phase of the vibration wave in each target blasting area , such as the angular frequency ω of the vibration wave in the i-th target blasting area i and the initial phase . Therefore, it is necessary to pre-set sensors in the target blasting area and the vibration prediction location to detect and calculate the angular frequency ω i and the initial phase . Of course, the angular frequency ω i and the initial phase can also be pre-estimated, so the implementation cost is relatively high

[0054] In practical applications, in order to simplify the calculation method, Formula Six or Formula Seven shown below can also be used to calculate the comprehensive vibration prediction value Formula Six Formula Seven In this Formula Six, V total-peak is the calculated comprehensive vibration prediction value. Obviously, the comprehensive vibration prediction value V calculated by this Formula Six total-peak is the direct superposition (rather than vector superposition) of the vibration prediction values of each target blasting area. Therefore, this comprehensive vibration prediction value V total-peak reflects the peak value of the vibration at the vibration prediction location. Obviously, the actual vibration situation at the vibration prediction location will be less than this comprehensive vibration prediction value V total-peak . Therefore, in this application, this comprehensive vibration prediction value V total-peak can also be used for vibration detection. For example, in practical applications, for the case of tunnel collaborative blasting, this comprehensive vibration prediction value V can be calculated by this Formula Six total-peak to predict the most extreme vibration situation that may occur, so as to give an early warning

[0055] In this Formula Seven, V total’ is the calculated comprehensive vibration prediction value; γ is a reduction coefficient, where the value of γ is greater than 0 and less than 1. Obviously, the comprehensive vibration prediction value V calculated by this Formula Seven total’ is V calculated by Formula Six total-peak multiplied by the reduction coefficient γ. Specifically, in engineering practice, V total-peak reflects the peak value of the vibration at the vibration prediction location. However, when multiple target blasting areas are collaboratively blasting, due to the differences in the angular frequency and initial phase of the vibration waves generated during the blasting of each target blasting area, the actual vibration situation will generally be less than this V total-peak . In this way, in this Formula Seven, by multiplying V total-peak by the reduction coefficient γ to obtain Vtotal’ , which can more accurately reflect the vibration conditions at the vibration prediction location.

[0056] Among them, the value of the reduction coefficient γ needs to be set according to the specific conditions of the tunnel, including factors such as the geological conditions around the tunnel and the spacing between the target blasting areas, and then determined in combination with relevant industry standards. Specifically, if the foundation around the tunnel is a rock foundation and the spacing between each target blasting area is small, the value of the reduction coefficient γ is relatively small and can be taken between 0.5 and 0.7; if the foundation around the tunnel is a soil foundation, the value of the reduction coefficient γ is relatively large and can be taken between 0.7 and 0.9.

[0057] In addition, obviously, the V calculated in Formula VII total’ compared with the V calculated in Formula V total in terms of, since this V total is obtained by vector superposition calculation, its accuracy will be higher, but the implementation cost is relatively higher. And the way to calculate V through Formula VII total’ although the accuracy is relatively low, the implementation cost is also relatively low. Therefore, in practical applications, according to actual needs, for example, when particularly accurate vibration prediction is required, the method of Formula V can be used to calculate the comprehensive vibration prediction value; when cost control is required, the methods of Formula VI or Formula VII can be used to calculate the comprehensive vibration prediction value.

[0058] In addition, for the scenario of multi-terminal interaction, for example, the user inputs the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area through the user operation interface on user-side electronic devices such as mobile phones and computers, and then in this step S11, the server can obtain the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area, and the basic prediction model. At this time, after the server calculates the comprehensive vibration prediction value through this step S14, the method can further include that the server feeds back the calculated comprehensive vibration prediction value to the user-side electronic device and displays it on the user operation interface of the user-side electronic device, so as to facilitate the user to check the comprehensive vibration prediction value in a timely manner.

[0059] Step S15: Use the comprehensive vibration prediction value for vibration early warning.

[0060] After obtaining the comprehensive vibration prediction value through the above-mentioned step S14, the comprehensive vibration prediction value can be used for vibration early warning. For example, a threshold can be preset in advance, and the comprehensive vibration prediction value is compared with the threshold. Specifically, it can be determined whether the comprehensive vibration prediction value is greater than or equal to the threshold. When the comprehensive vibration prediction value is greater than or equal to the threshold, vibration early warning is carried out, such as emitting an alarm sound or an alarm indicator light, etc. On the contrary, when the comprehensive vibration prediction value is less than the threshold, vibration early warning may not be carried out. Among them, for the setting method of the threshold, for example, it can be set according to the building safety standard of the vibration prediction location, or it can also be set to a certain estimated value according to experience.

[0061] Adopting the tunnel blasting vibration prediction and early warning method provided by the embodiments of the present application, the tunnel is divided into multiple sections along the length direction or the cross-section direction, and target blasting areas are respectively selected from at least two of the multiple sections for collaborative blasting. Among them, the method includes obtaining a basic prediction model and the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area, and then using the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area to calculate the explosive parameter correction term, blasting parameter correction term, and geological parameter correction term of each target blasting area respectively. Then, through the basic prediction model and the explosive parameter correction term, blasting parameter correction term, and geological parameter correction term of each target blasting area, the vibration prediction value of each target blasting area is calculated respectively. Then, through the vibration prediction values of each target blasting area, the comprehensive vibration prediction value is calculated. Then, the comprehensive vibration prediction value is used for vibration early warning. This method is aimed at the scenario of collaborative blasting of multiple target blasting areas during the tunnel construction process. The comprehensive vibration prediction value calculated through the vibration prediction values of each target blasting area can more accurately reflect the vibration situation of the vibration prediction location, so it can solve the problems in the prior art.

[0062] Based on the same inventive concept as the tunnel blasting vibration prediction and early warning method provided by the embodiments of the present application, the embodiments of the present application can also provide a tunnel blasting vibration prediction and early warning system. Among them, the tunnel is divided into multiple sections along the length direction or the cross-section direction, and target blasting areas are respectively selected from at least two of the multiple sections for collaborative blasting. As Figure 2 shown in the specific structural schematic diagram of the system 20, the system 20 includes: an acquisition unit 201, a first calculation unit 202, a second calculation unit 203, a third calculation unit 204, and an early warning unit 205, where: The acquisition unit 201 is used to acquire a basic prediction model and the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area; The first calculation unit 202 is configured to calculate the explosive parameter correction term, the blasting parameter correction term, and the geological parameter correction term for each target blasting area respectively by using the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area. The second calculation unit 203 is configured to calculate the vibration prediction value for each target blasting area respectively through the basic prediction model and the explosive parameter correction term, the blasting parameter correction term, and the geological parameter correction term of each target blasting area. The third calculation unit 204 is configured to calculate the comprehensive vibration prediction value through the vibration prediction values of each target blasting area. The early warning unit 205 is configured to perform vibration early warning by using the comprehensive vibration prediction value.

[0063] By adopting the system 20 provided in the embodiment of the present application, since the system 20 adopts the same inventive concept as the method provided in the embodiment of the present application, when the method can solve the problems in the prior art, the system 20 can also solve the problems in the prior art, and details are not described herein again.

[0064] Wherein, the first calculation unit 202 can specifically calculate the explosive parameter correction term for each target blasting area in the following manner: ; Where: i is the number of the target blasting area; D i is the actual detonation velocity of the explosive in the i-th target blasting area; Q ei is the actual heat of detonation of the explosive in the i-th target blasting area; ρ i is the actual density of the explosive in the i-th target blasting area; D 0 is the detonation velocity of the standard explosive; Q e0 is the heat of detonation of the standard explosive; ρ 0 is the density of the standard explosive; β 1i is the explosive parameter correction term for the i-th target blasting area.

[0065] Wherein, the first calculation unit 202 can specifically calculate the blasting parameter correction term for each target blasting area in the following manner: ; Where: i is the number of the target blasting area; m i is the initiation method related parameter of the i-th target blasting area; c i is the actual hole spacing between blast holes in the i-th target blasting area; b i is the actual row spacing between blast holes in the i-th target blasting area; L i is the actual hole depth of the blast hole in the i-th target blasting area; c 0 is the hole spacing of the standard blast hole; b0 is the row spacing of standard blast holes; L 0 is the hole depth of standard blast holes; β 2i is the blasting parameter correction term for the i-th target blasting area.

[0066] Among them, the first calculation unit 202 can specifically calculate the geological parameter correction terms of each target blasting area in the following way: ; where: i is the number of the target blasting area; p i is the geological structure correction coefficient of the i-th target blasting area; E i is the actual elastic modulus of the rock and soil mass in the i-th target blasting area; u i is the actual Poisson's ratio of the rock and soil mass in the i-th target blasting area; ρ ri is the actual density of the rock and soil mass in the i-th target blasting area; E 0 is the elastic modulus of the standard rock and soil mass; u 0 is the Poisson's ratio of the standard rock and soil mass; ρ r0 is the density of the standard rock and soil mass; β 3i is the geological parameter correction term of the i-th target blasting area.

[0067] Among them, by using the basic prediction model and the explosive parameter correction terms, blasting parameter correction terms and geological parameter correction terms of each target blasting area, the vibration prediction values of each target blasting area can be calculated respectively, which can specifically include calculating the vibration prediction values of each target blasting area through the following formula: V final-i = V i × β 1i × β 2i ×β 3i ; where: i is the number of the target blasting area; V final-i is the vibration prediction value of the i-th target blasting area calculated; β 1i is the explosive parameter correction term of the i-th target blasting area; β 2i is the blasting parameter correction term of the i-th target blasting area; β 3i is the geological parameter correction term of the i-th target blasting area; Vi is the influence amount of the blasting vibration velocity of the i-th target blasting area on the vibration prediction location calculated by the basic prediction model.

[0068] Among them, by using the vibration prediction values of each target blasting area, the comprehensive vibration prediction value can be calculated, which can specifically include calculating the comprehensive vibration prediction value through the following formula: ; where, Vtotal is the calculated comprehensive vibration prediction value; V final-i is the vibration prediction value of the i-th target blasting area; n is the total number of target blasting areas; ω i is the angular frequency of the vibration wave in the i-th target blasting area; is the initial phase of the vibration wave in the i-th target blasting area; t is the time when vibration detection is carried out, which is the time interval between the time of vibration detection and the time of collaborative blasting.

[0069] Among them, calculating the comprehensive vibration prediction value through the vibration prediction values of each target blasting area may specifically include calculating the comprehensive vibration prediction value through the following formula: ; Among them, V total’ is the calculated comprehensive vibration prediction value; V final-i is the vibration prediction value of the i-th target blasting area; n is the total number of target blasting areas; γ is a reduction coefficient, where the value of γ is greater than 0 and less than 1.

[0070] Among them, the basic prediction model may specifically be the Sadovsky formula.

[0071] Among them, using the comprehensive vibration prediction value for vibration early warning may specifically include: judging whether the comprehensive vibration prediction value is greater than or equal to the threshold; in the case where the comprehensive vibration prediction value is greater than or equal to the threshold, carrying out vibration early warning.

[0072] Figure 3 Illustrates a schematic physical structure diagram of an electronic device, such as Figure 3As shown in the figure, the electronic device may include: a processor 310, a communications interface 320, a memory 330, and a communication bus 340. Among them, the processor 310, the communications interface 320, and the memory 330 complete mutual communication through the communication bus 340. The processor 310 may call the logical instructions in the memory 330 to execute the method for predicting and warning tunnel blasting vibration provided in the embodiments of the present application, including obtaining a basic prediction model and the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area, and then using the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area to calculate the explosive parameter correction term, blasting parameter correction term, and geological parameter correction term of each target blasting area respectively. Then, through the basic prediction model and the explosive parameter correction term, blasting parameter correction term, and geological parameter correction term of each target blasting area, calculate the vibration prediction value of each target blasting area respectively. Then, calculate the comprehensive vibration prediction value through the vibration prediction values of each target blasting area, and then use the comprehensive vibration prediction value for vibration warning.

[0073] Among them, in practical applications, the electronic device may be an electronic device at the user end or an electronic device at the server end.

[0074] Obviously, since the processor 310 can call the logical instructions in the memory 330 to execute the method provided in the embodiments of the present application, this method is for the scenario of collaborative blasting in multiple target blasting areas during the tunnel construction process. The comprehensive vibration prediction value calculated through the vibration prediction values of each target blasting area can more accurately reflect the vibration situation at the vibration prediction location, so it can solve the problems in the prior art.

[0075] In addition, when the logical instructions in the above-mentioned memory 330 are implemented in the form of software functional units and sold or used as an independent product, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art or a part of this technical solution can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), magnetic disks, or optical discs that can store program codes.

[0076] On the other hand, the present invention also provides a computer program product, which includes a computer program. The computer program can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the prediction and early warning method for tunnel blasting vibration provided by the embodiments of the present application, including obtaining a basic prediction model and the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area, and then using the actual explosive parameters, actual blasting parameters, and actual geological parameters of each target blasting area to calculate the explosive parameter correction term, blasting parameter correction term, and geological parameter correction term of each target blasting area respectively. Then, through the basic prediction model and the explosive parameter correction term, blasting parameter correction term, and geological parameter correction term of each target blasting area, calculate the vibration prediction value of each target blasting area respectively. Then, calculate the comprehensive vibration prediction value through the vibration prediction values of each target blasting area, and then use the comprehensive vibration prediction value for vibration early warning.

[0077] Obviously, since the computer can execute the method provided by the embodiments of the present application when the computer program is executed by the processor, and this method is for the scenario of collaborative blasting in multiple target blasting areas during the tunnel construction process, the comprehensive vibration prediction value calculated through the vibration prediction values of each target blasting area can more accurately reflect the vibration situation at the vibration prediction location. Therefore, it can solve the problems in the prior art.

[0078] On another aspect, the present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, it is the method provided by the embodiments of the present application.

[0079] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separated. The components shown as units may or may not be physical units, that is, they may be located in one place or distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. Those of ordinary skill in the art can understand and implement it without creative labor.

[0080] Through the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on such an understanding, the essence of the above technical solution, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to enable a computer device (which can be a personal computer, server, or network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0081] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements for some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for predicting and warning tunnel blasting vibration, characterized in that: The tunnel is divided into a plurality of sections along the length direction or the cross-section direction, and target blasting areas are respectively selected from at least two of the plurality of sections for coordinated blasting. The method comprises: Obtain the basic prediction model and actual explosive parameters, actual blasting parameters and actual geological parameters of each target blasting area; Using actual explosive parameters, actual blasting parameters and actual geological parameters of each target blasting area, respectively calculate the explosive parameter correction item, blasting parameter correction item and geological parameter correction item of each target blasting area; Calculate the vibration prediction value of each target blasting area respectively by using the basic prediction model and the explosive parameter correction item, blasting parameter correction item and geological parameter correction item of each target blasting area; Calculate the comprehensive vibration prediction value through the vibration prediction values ​​of each target blasting area; The comprehensive vibration prediction value is used to provide vibration warning.

2. The method according to claim 1, characterized in that The explosive parameter correction items for each target blasting area are calculated by the following method: ； Where: i is the number of the target blasting area; D i is the actual detonation velocity of the explosive in the i-th target blasting area; Q ei is the actual explosion heat of the explosive in the i-th target blasting area; ρ i is the actual density of explosives in the i-th target blasting area; D0 is the detonation velocity of standard explosives; Q e0 is the explosion heat of standard explosive; ρ0 is the density of standard explosive; β 1i is the explosive parameter correction item for the i-th target blasting area.

3. The method according to claim 1, characterized in that The blasting parameter correction items for each target blasting area are calculated by the following method: ； Where: i is the number of the target blasting area; m i is the associated parameter of the detonation mode of the i-th target blasting area; c i is the actual blasthole spacing between blastholes in the i-th target blasting area; b i is the actual row spacing between blastholes in the i-th target blasting area; L i is the actual depth of the blasthole in the i-th target blasting area; c0 is the blasthole spacing of the standard blasthole; b0 is the row spacing of the standard blasthole; L0 is the depth of the standard blasthole; β 2i is the blasting parameter correction item of the i-th target blasting area.

4. The method according to claim 1, characterized in that: The geological parameter correction items for each target blasting area are calculated by the following method: ； Where: i is the number of the target blasting area; p i is the geological structure correction coefficient of the i-th target blasting area; E i is the actual elastic modulus of the rock and soil in the i-th target blasting area; u i is the actual Poisson's ratio of the rock and soil in the i-th target blasting area; ρ ri is the actual density of the rock and soil in the i-th target blasting area; E0 is the elastic modulus of the standard rock and soil; u0 is the Poisson's ratio of the standard rock and soil; ρ r0 is the density of standard rock and soil mass; 3i is the geological parameter correction item of the i-th target blasting area.

5. The method according to claim 1, characterized in that The vibration prediction value of each target blasting area is calculated respectively by using the basic prediction model and the explosive parameter correction item, blasting parameter correction item and geological parameter correction item of each target blasting area, specifically including calculating the vibration prediction value of each target blasting area respectively by the following formula: V final-i = V i × b 1i × b 2i ×b 3i ; Where: i is the number of the target blasting area; V final-i is the vibration prediction value of the i-th target blasting area calculated; β 1i is the explosive parameter correction term for the i-th target blasting area; β 2i is the blasting parameter correction term of the i-th target blasting area; β 3i is the geological parameter correction term of the i-th target blasting area; V i It is the influence of the blasting vibration velocity of the i-th target blasting area on the vibration prediction location, calculated by the basic prediction model.

6. The method according to claim 1, characterized in that The comprehensive vibration prediction value is calculated by using the vibration prediction values ​​of each target blasting area, specifically including calculating the comprehensive vibration prediction value by the following formula: ； Among them, V total is the calculated comprehensive vibration prediction value; V final-i is the vibration prediction value of the i-th target blasting area; n is the total number of target blasting areas; ω i is the angular frequency of the vibration wave in the ith target blasting area; is the initial phase of the vibration wave in the i-th target blasting area; t is the time interval between the moment of vibration detection and the moment of coordinated blasting.

7. The method according to claim 1, characterized in that The comprehensive vibration prediction value is calculated by using the vibration prediction values ​​of each target blasting area, specifically including calculating the comprehensive vibration prediction value by the following formula: ； Among them, V total’ is the calculated comprehensive vibration prediction value; V final-i is the vibration prediction value of the i-th target blasting area; n is the total number of target blasting areas; γ is the reduction coefficient, where the value of γ is greater than 0 and less than 1.

8. The method according to claim 1, characterized in that: The basic prediction model is specifically the Sadovsky formula.

9. The method according to claim 1, characterized in that: Using the comprehensive vibration prediction value to perform vibration early warning specifically includes: Determining whether the comprehensive vibration prediction value is greater than or equal to a threshold; When the comprehensive vibration prediction value is greater than or equal to the threshold, a vibration warning is issued.

10. A tunnel blasting vibration prediction and warning system, characterized in that: The tunnel is divided into a plurality of sections along the length direction or the cross-section direction, and target blasting areas are respectively selected from at least two of the plurality of sections for coordinated blasting. The system comprises: An acquisition unit, used to acquire the basic prediction model and actual explosive parameters, actual blasting parameters and actual geological parameters of each target blasting area; The first calculation unit is used to calculate the explosive parameter correction item, the blasting parameter correction item and the geological parameter correction item of each target blasting area respectively by using the actual explosive parameter, the actual blasting parameter and the actual geological parameter of each target blasting area; A second calculation unit is used to calculate the vibration prediction value of each target blasting area respectively by using the basic prediction model and the explosive parameter correction item, blasting parameter correction item and geological parameter correction item of each target blasting area; A third calculation unit is used to calculate a comprehensive vibration prediction value through the vibration prediction values ​​of each target blasting area; The early warning unit is used to use the comprehensive vibration prediction value to provide a vibration early warning.

Images