Method for calculating subsection blasting parameters in deep hole step hole of underground cavern
By conducting in-situ single-hole segmented blasting tests and function fitting, the time delay interval between holes within the hole was optimized, solving the problems of large vibration and low borehole utilization in deep-hole bench blasting. This achieved precise blasting vibration control and reduced explosive consumption per unit.
Patent Information
- Application Number
- CN202510997082.X
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-07-18
- Publication Date
- 2025-11-07
AI Technical Summary
Deep-hole bench blasting presents challenges such as large charge per hole, significant vibration during multi-hole blasting, and low borehole utilization. Existing technologies struggle to achieve vibration synthesis calculations and precise control.
Measured single-segment waveforms were obtained through in-situ single-hole segmented blasting tests. Function fitting and waveform reconstruction were performed to determine the optimal time intervals within and between holes. The segmented blasting parameters within the hole were then optimized using the Anderson linear superposition method.
It achieves precise control of deep-hole bench blasting vibration, reduces blasting vibration peak value, improves borehole utilization, and alleviates the problems of damage to surrounding rock and high explosive consumption caused by group-controlled vibration.
Smart Images
Figure CN120907388A_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the technical field of underground cavern excavation engineering, in particular to a kind of underground cavern deep hole bench hole inside segmented blasting parameter calculation method. BACKGROUND
[0002] As the core means of efficient rock breaking, the demand for fine control of blasting technology is increasingly prominent. In deep hole bench blasting, there are problems such as large single-hole charge, large multi-hole blasting vibration, and low utilization rate of blast holes. The in-hole segmented blasting technology is used to solve these problems. In the scenario of in-hole segmented blasting, there is a complex coupling relationship between the charge structure, charge, delay interval and synthesized vibration. Therefore, it is of great significance to establish an in-hole segmented blasting vibration synthesis method to achieve accurate prediction and control of blasting vibration effect.
[0003] For deep hole bench blasting vibration reduction, there are two aspects: one is the optimization of charge structure; the other is to select the optimal electronic detonator delay interval for vibration synthesis. In actual engineering, the two vibration reduction technologies are combined to form in-hole millisecond blasting technology. However, the current problem to be solved is the in-hole segmented multi-hole blasting vibration synthesis calculation method for deep hole bench blasting. The in-hole segmented vibration synthesis calculation method and the inter-hole vibration synthesis calculation method need to be combined together, and reasonable in-hole delay interval and inter-hole delay interval can achieve the effect of double vibration reduction, greatly reducing the vibration caused by electronic detonator blasting and improving the safety of blasting. However, there are few related application researches, and there are technical problems such as excessive group control vibration damage to surrounding rock and high explosive specific energy consumption. SUMMARY
[0004] To solve the above technical problems existing in the prior art, the present application provides an underground cavern deep hole bench hole inside segmented blasting parameter calculation method. The technical solution is as follows: The present application provides an underground cavern deep hole bench hole inside segmented blasting parameter calculation method, which comprises: based on in-situ single-hole in-hole segmented blasting test, obtaining the measured single-segment waveform generated by in-hole blasting of the target underground cavern deep hole bench; the measured single-segment waveform includes measured hole mouth single-segment waveform and measured hole bottom single-segment waveform; function fitting is performed on the measured single-segment waveform to obtain a single-segment blasting waveform function; the single-segment blasting waveform function includes hole mouth single-segment blasting waveform function and hole bottom single-segment blasting waveform function; based on the relationship between vibration velocity and charge, the single-segment blasting waveform function is reconstructed to obtain reconstructed single-segment waveforms under different charges; the reconstructed single-segment waveforms include hole mouth reconstructed single-segment waveforms and hole bottom reconstructed single-segment waveforms; based on different in-hole delay intervals, the hole mouth reconstructed single-segment waveforms and the hole bottom reconstructed single-segment waveforms are superimposed to obtain complete single-hole waveforms, and the optimal in-hole delay interval is determined; based on different inter-hole delay intervals, the complete single-hole waveforms are superimposed to obtain multi-hole waveforms, and the optimal inter-hole delay interval is determined.
[0005] Optionally, based on the in-situ single-hole in-hole segmented blasting test, the measured single-segment waveform generated by the in-hole blasting of the target underground cavern deep hole step is obtained, including: setting an in-situ single-hole in-hole segmented blasting test at the target underground cavern deep hole step, and setting a delay interval for the blasting of the orifice section and the bottom section in the test hole; based on the blasting vibration monitor, the measured single-segment waveform generated by the in-hole blasting of the target underground cavern deep hole step is measured and obtained.
[0006] Optionally, the measured single-segment waveform is functionally fitted to obtain a single-segment blasting waveform function, including: continuously fitting the measured single-segment waveform, and using Fourier transform for functional fitting to obtain a single-segment blasting waveform function.
[0007] Optionally, the single-segment blasting waveform function is waveform reconstructed based on the relationship between vibration velocity and charge amount to obtain a reconstructed single-segment waveform under different charge amounts, including: determining the relationship between vibration velocity and charge amount based on the Sadovskiy formula; calculating different vibration velocity and charge amount proportionality coefficients based on the relationship between vibration velocity and charge amount; reconstructing the single-segment blasting waveform function based on the proportionality coefficients to obtain a reconstructed single-segment waveform under different charge amounts.
[0008] Optionally, the orifice reconstructed single-segment waveform includes:
[0009] The bottom reconstructed single-segment waveform includes:
[0010] Wherein,
[0011]
[0012]
[0013] In the formula, n is the proportionality coefficient, v1 and v2 are different blasting vibration velocities; Q1 and Q2 are the charge amounts corresponding to v1 and v2 respectively; R is the blast center distance; K and a are coefficients and attenuation exponents related to the terrain and geological conditions between the blast source and the protected object; f1(t) and f2(t) are the orifice single-segment blasting waveform function and the bottom single-segment blasting waveform function respectively.
[0014] Optionally, based on the complete single-hole waveforms obtained by superimposing the orifice reconstructed single-segment waveforms and the bottom reconstructed single-segment waveforms at different intra-hole delay time intervals, the optimal intra-hole delay time interval is determined, comprising: superimposing the orifice reconstructed single-segment waveforms and the bottom reconstructed single-segment waveforms at different intra-hole delay time intervals to obtain complete single-hole waveforms of a single blast hole; taking the intra-hole delay time interval corresponding to the complete single-hole waveform with the minimum peak value of blasting vibration as the optimal intra-hole delay time interval; wherein the complete single-hole waveform comprises:
[0015] wherein f hole (t) is the complete single-hole waveform, f 1Q (t) and f 2Q (t) are the orifice reconstructed single-segment waveform and the bottom reconstructed single-segment waveform respectively, and Δt0 is the intra-hole delay time interval.
[0016] Optionally, based on the multi-hole waveforms obtained by superimposing the complete single-hole waveforms at different inter-hole delay time intervals, the optimal inter-hole delay time interval is determined, comprising: superimposing the complete single-hole waveforms based on the Anderson linear superposition method to obtain multi-hole intra-hole segmented blasting waveform functions corresponding to different inter-hole delay time intervals; taking the inter-hole delay time interval corresponding to the multi-hole intra-hole segmented blasting waveform function with the minimum peak value of blasting vibration as the optimal inter-hole delay time interval; wherein the multi-hole intra-hole segmented blasting waveform function comprises:
[0017] wherein F holes (t, Δt1) is the multi-hole intra-hole segmented blasting waveform function, Δt1 is the inter-hole delay time interval, i is the number of superimposed waveforms, and m is the number of superimposed waveforms.
[0018] The underground cavern deep hole bench hole intra-segmented blasting parameter calculation method provided by the embodiments of the present application can avoid problems such as excessive multi-hole vibration and low utilization rate of blast holes in deep hole bench blasting, has important significance and application value for accurate control of underground engineering deep hole bench blasting vibration, and alleviates the technical problems of excessive group control vibration damaging surrounding rock and high explosive specific energy in the prior art. BRIEF DESCRIPTION OF DRAWINGS
[0019] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed in the embodiment description will be briefly introduced as follows. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor on the basis of these drawings.
[0020] Figure 1is a flow chart of a deep hole bench hole internal segmented blasting parameter calculation method of a underground cavern provided by an embodiment of the present application. Figure 2 is a waveform analysis block diagram of a deep hole bench hole internal segmented blasting parameter calculation method of a underground cavern provided by an embodiment of the present application. Figure 3 is a hole internal segmented millisecond blasting vibration curve synthesis principle schematic diagram provided by an embodiment of the present application. Figure 4 is a multi-hole millisecond blasting vibration curve synthesis principle schematic diagram provided by an embodiment of the present application. Figure 5 is a cross section schematic diagram of a super-large cross section cavern group provided by an embodiment of the present application. Figure 6 is a space position schematic diagram of a blasting source and a nearest building provided by an embodiment of the present application. Figure 7 is a test hole setting scheme schematic diagram provided by an embodiment of the present application. Figure 8 is a blasting scheme schematic diagram provided by an embodiment of the present application. DETAILED DESCRIPTION
[0021] The technical solutions in the present application will be described below with reference to the drawings.
[0022] In the embodiments of the present application, the words such as "example", "for example" are used to represent an example, illustration or description. Any embodiment or design scheme described as "example" in the present application should not be interpreted as more preferred or more advantageous than other embodiments or design schemes. Rather, the word "example" is intended to present the concept in a specific manner. In addition, in the embodiments of the present application, the meaning expressed by "and / or" can be both, or can be one of the two.
[0023] To make the technical problems, technical solutions and advantages of the present application clearer, the following will be described in detail with reference to the drawings and specific embodiments.
[0024] Figure 1 is a flow chart of a deep hole bench hole internal segmented blasting parameter calculation method of a underground cavern provided by an embodiment of the present application. As shown in Figure 1 , the method specifically includes the following steps: Step S102, based on an in-situ single-hole hole internal segmented blasting test, obtaining a measured single-segment waveform generated by hole blasting of a target deep hole bench of a underground cavern; the measured single-segment waveform includes a measured hole mouth single-segment waveform and a measured hole bottom single-segment waveform.
[0025] Step S104, function fitting is performed on the measured single-stage waveform to obtain a single-stage blasting waveform function; the single-stage blasting waveform function includes an orifice single-stage blasting waveform function and a hole bottom single-stage blasting waveform function.
[0026] Step S106, based on the relationship between the vibration velocity and the charge amount, the single-stage blasting waveform function is reconstructed to obtain a reconstructed single-stage waveform under different charge amounts; the reconstructed single-stage waveform includes an orifice reconstructed single-stage waveform and a hole bottom reconstructed single-stage waveform.
[0027] Step S108, based on the complete single-hole waveform obtained by superimposing the orifice reconstructed single-stage waveform and the hole bottom reconstructed single-stage waveform at different hole internal delay time intervals, an optimal hole internal delay time interval is determined.
[0028] Step S110, based on the multi-hole waveform obtained by superimposing the complete single-hole waveform at different hole-to-hole delay time intervals, an optimal hole-to-hole delay time interval is determined.
[0029] Figure 2 is a waveform analysis block diagram of a method for calculating underground cavern deep hole bench step hole internal segmented blasting parameters according to an embodiment of the present application. As shown in Figure 2 , step S102 further includes the following steps: Step S1021, an in-situ single-hole internal segmented blasting test is set at a target underground cavern deep hole bench, and a delay time interval for blasting of an orifice section and a hole bottom section in the test hole is set; Step S1022, based on a blasting vibration monitor, a measured single-stage waveform generated by in-hole blasting of the target underground cavern deep hole bench is measured and obtained. The blasting vibration monitor is placed at a protection object.
[0030] Specifically, as shown in Figure 2 , step S104 includes: continuousizing the measured single-stage waveform, and using Fourier transform to perform function fitting to obtain a single-stage blasting waveform function.
[0031] Specifically, the mathematical form of the Fourier transform fitting function is as follows:
[0032]
[0033] In the formula, g ( t ) is a single-stage blasting waveform function; t is time; n is a series; ω is a fundamental frequency; a 0、 a j 、 b j is a fitting coefficient; N is the number of sampling points;t 0 is the single-segment waveform wavelength time; k is the half point; t k is the sampling point abscissa; g t k is the sampling point ordinate.
[0034] Then the Fourier transform fitting function is extended in the time domain to obtain a single-segment blast waveform function:
[0035] In the formula, f(t) is a single-segment blast waveform function.
[0036] Specifically, step S106 further includes the following steps: Step S1061, determining the relationship between vibration velocity and charge based on the Sadovskiy formula; wherein the Sadovskiy formula includes:
[0037] In the formula, v is the blast vibration velocity, Q is the charge corresponding to v, R is the blast center distance; K and a are coefficients and attenuation exponents related to the terrain and geological conditions between the blast source and the protected object.
[0038] Step S1062, calculating different vibration velocity and charge proportionality coefficients based on the relationship between vibration velocity and charge; wherein the proportionality coefficient n includes:
[0039] In the formula, v1 and v2 are different blast vibration velocities; Q1 and Q2 are the charges corresponding to v1 and v2, respectively.
[0040] Step S1063, reconstructing the single-segment blast waveform function based on the proportionality coefficient to obtain a reconstructed single-segment waveform under different charges.
[0041] Specifically, the orifice reconstructed single-segment waveform includes:
[0042] The bottom of the hole reconstructed single-segment waveform includes:
[0043] In the formula, f1(t) and f2(t) are the orifice single-segment blast waveform function and the bottom of the hole single-segment blast waveform function, respectively.
[0044] After obtaining the reconstructed single-segment waveform under different charges, the actual use charge can be calculated according to the test charge, for example, Q2 is the test charge, Q1 is the actual use charge, and the actual waveform can be calculated according to the reconstructed single-segment waveform.
[0045] Figure 3 is a hole in the segmented differential blasting vibration curve synthesis principle diagram provided by the embodiment of the application. As shown in Figure 3 , the step S108 further includes the following steps: Step S1081, superimposing the orifice reconstructed single-section waveform and the hole bottom reconstructed single-section waveform according to different hole delay time intervals to obtain a complete single-hole waveform of a single blast hole; Step S1082, determining the hole delay time interval corresponding to the complete single-hole waveform with the minimum blasting vibration peak value as the optimal hole delay time interval; wherein, the complete single-hole waveform includes:
[0046] In the formula, f hole (t) is a complete single-hole waveform, f 1Q (t) and f 2Q (t) are the orifice reconstructed single-section waveform and the hole bottom reconstructed single-section waveform respectively, and Δt0 is the hole delay time interval.
[0047] Figure 4 is a multi-hole differential blasting vibration curve synthesis principle diagram provided by the embodiment of the application. As shown in Figure 4 , the step S110 further includes the following steps: Step S1101, superimposing the complete single-hole waveform based on the Anderson linear superposition method to obtain a multi-hole hole-in-segmented blasting waveform function corresponding to different hole delay time intervals; Specifically, based on the Anderson linear superposition theory, the single-hole waveform is superimposed and calculated to obtain the multi-hole waveform function as follows:
[0048] In the formula, Δt is the delay time interval; F(t, Δt) is the predicted waveform function obtained by superposition; i is the number of superimposed waveforms, and m is the number of superimposed waveforms.
[0049] Based on the above formula, the complete single-hole waveform is superimposed:
[0050] In the formula, F holes (t, Δt1) is a multi-hole hole-in-segmented blasting waveform function, and Δt1 is a hole delay time interval.
[0051] Step S1102, determining the hole delay time interval corresponding to the multi-hole hole-in-segmented blasting waveform function with the minimum blasting vibration peak value as the optimal hole delay time interval.
[0052] The application will be described in detail below according to specific implementation methods in combination with examples and drawings. The following examples are used to illustrate the application but do not limit the use range of the application.
[0053] The engineering relied on by the application is a certain under-construction super-large cross-section cavern group, the cross-section size of the main cavern is 30m (high) x 22m (wide), the total excavation area is 616m 2 , which belongs to a super-large cross-section, as shown in Figure 5 . The cavern group passes through a village from below and dense civil buildings from above, the spatial position of the blast source and the closest building is shown in Figure 6 , the closest building to the blast source is 331m, in the early pre-feasibility study stage, the safety threshold of the peak value of the blasting vibration at the civil building is 0.35cm / s, therefore, strict controlled blasting is adopted for the excavation of the cavern group, especially the middle layer vertical bench blasting, which has large excavation volume, large explosive quantity and large vibration control difficulty. The application is a method proposed for the problem, and the specific implementation method comprises the following steps: (1) Obtain the measured single-segment waveforms at the hole mouth and the hole bottom through in-situ single-hole in-hole segmented blasting test: First, design a verification test, on the basis of the conventional design scheme of the bench blasting, set the test hole, as shown in Figure 7 , the hole mouth of the last blast hole and the previous blast hole and the delay time of the in-hole blasting of the last blast hole are set as 200ms, the TC-4850 blasting vibration monitor is used to monitor the peak value of the blasting vibration in real time, and the measured single-segment waveforms at the hole mouth and the hole bottom are obtained.
[0054] (2) Fit the measured single-segment waveforms to determine the single-segment waveform function expression of the hole mouth segment and the hole bottom segment: Use the software MATLAB to fit the sampling data of the measuring points, and obtain the continuous hole mouth segment f 1 t and the hole bottom segment f 2 t , the two waveforms obtained in the test stage are used as the basic waveforms for calculating the in-hole segmented single-segment waveforms in the field application.
[0055] (3) Reconstruct the measured single-segment waveforms to obtain single-segment waveforms under different explosive quantities: In order to obtain the single-segment waveforms of the hole mouth segment and the hole bottom segment in the field application, the fitted hole mouth segment and hole bottom segment waveforms are reconstructed, the field application blasting scheme is shown in Figure 8 , the hole mouth segment explosive quantity in the verification test is 15.3kg, the hole bottom segment explosive quantity is 10.8kg, the hole mouth segment explosive quantity in the field application is 12.6kg, and the hole bottom segment explosive quantity is 11.7kg, according to the reconstructed single-segment waveforms under different explosive quantities, the proportional coefficient n 1 of the hole mouth explosive quantity and the vibration is (10.8 / 11.7) α / 3The ratio of the amount of explosive at the bottom of the hole to the vibration. n 2 is (12.6 / 15.3) α / 3 Based on the field-measured vibration data, the results were obtained by fitting the data using the Sadovsky formula. α The value is 1.478. Therefore, the single-segment waveform of the orifice section is... n 1 f 1( t The single-segment waveform of the bottom section of the hole is as follows: n 2 f 2( t ).
[0056] (4) By superimposing the reconstructed waveforms of the orifice opening and bottom segments according to different orifice delay intervals, a complete single-orifice waveform is obtained, and the optimal orifice delay interval is determined: The two waveforms are arranged according to the delay time Δ t 0 superposition forms a single-hole waveform Δ t The superimposed single-hole waveform under different delay intervals was calculated by taking integer values from 0 to 20ms. After comparing the peak vibration velocity of the single-hole waveform corresponding to each delay interval, the peak vibration velocity reached the minimum value of 0.11cm / s when the delay interval was 6ms. Therefore, the optimal delay interval in the hole for this blasting was 6ms.
[0057] (5) Superimpose the single-hole waveform with different inter-hole delay intervals to obtain the multi-hole waveform, and determine the optimal inter-hole delay interval: After obtaining the optimal intra-hole delay interval, the delay interval Δ of the inter-hole electronic detonator is further calculated. t 1. Based on the synthesis formula of the segmented blasting waveform function in a multi-hole borehole, the peak vibration value of a 32-hole multi-hole blasting is: ,in Δ t 1. Values were taken sequentially from 0 to 200 ms, with a step size of 10 ms. After calculating and comparing the peak vibration velocity of the multi-hole waveform corresponding to each delay interval, the peak vibration velocity reached its minimum value of 0.25 cm / s at a delay interval of 50 ms. Therefore, the optimal inter-hole delay interval for this blasting operation is 50 ms. In summary, for field application, the designed intra-hole delay interval is 6 ms, and the inter-hole delay interval is 50 ms.
[0058] The above is a specific implementation example of the method of the present invention. The present invention is not limited to specific instances. Any calculation using the same method should be considered as the same type of method of the present invention.
[0059] The above merely illustrates the specific embodiments of the present application, but the protection scope of the present application is not limited thereto, any person skilled in the art can easily think of the changes or replacements within the technical range disclosed by the present application, which should be covered in the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.
Claims
1. A method for calculating parameters of segmented blasting in a deep hole bench hole of an underground cavern, characterized in that, The method comprises: Based on the in-situ single-hole in-hole segmented blasting test, the measured single-segment waveform generated by the in-hole blasting of the target underground cavern deep hole step is obtained; the measured single-segment waveform comprises a measured hole mouth single-segment waveform and a measured hole bottom single-segment waveform; The measured single-segment waveform is functionally fitted to obtain a single-segment blasting waveform function; the single-segment blasting waveform function comprises a hole mouth single-segment blasting waveform function and a hole bottom single-segment blasting waveform function; Based on the relationship between vibration velocity and charge quantity, the single-segment blasting waveform function is waveform reconstructed to obtain reconstructed single-segment waveforms under different charge quantities; the reconstructed single-segment waveforms comprise hole mouth reconstructed single-segment waveforms and hole bottom reconstructed single-segment waveforms; Based on different in-hole delay time intervals, the hole mouth reconstructed single-segment waveforms and the hole bottom reconstructed single-segment waveforms are superimposed to obtain complete single-hole waveforms, and optimal in-hole delay time intervals are determined; Based on different inter-hole delay time intervals, the complete single-hole waveforms are superimposed to obtain multi-hole waveforms, and optimal inter-hole delay time intervals are determined.
2. The method of claim 1, wherein, Based on the in-situ single-hole in-hole segmented blasting test, the measured single-segment waveform generated by the in-hole blasting of the target underground cavern deep hole step is obtained, comprising: An in-situ single-hole in-hole segmented blasting test is arranged at the target underground cavern deep hole step, and a delay time interval for hole mouth segment and hole bottom segment blasting in the test hole is arranged; Based on a blasting vibration monitor, the measured single-segment waveform generated by the in-hole blasting of the target underground cavern deep hole step is measured and obtained.
3. The method of claim 1, wherein, The measured single-segment waveform is functionally fitted to obtain a single-segment blasting waveform function, comprising: the measured single-segment waveform is continuously fitted, and Fourier transform is used for function fitting to obtain a single-segment blasting waveform function.
4. The method of claim 1, wherein, Based on the relationship between vibration velocity and charge quantity, the single-segment blasting waveform function is waveform reconstructed to obtain reconstructed single-segment waveforms under different charge quantities, comprising: Based on the Sadovnikov formula, the relationship between vibration velocity and charge quantity is determined; Based on the relationship between vibration velocity and charge quantity, different vibration velocity and charge quantity proportionality coefficients are calculated; Based on the proportionality coefficient, the single-segment blasting waveform function is reconstructed to obtain reconstructed single-segment waveforms under different charge quantities.
5. The method of claim 4, wherein, The hole mouth reconstructed single-segment waveforms comprise: The hole bottom reconstructed single-segment waveforms comprise: Wherein, In the formula, n is the proportionality coefficient, v1 and v2 are different blasting vibration velocities; Q1 and Q2 are the charge quantities corresponding to v1 and v2 respectively; R is the blast center distance; K and a are the coefficients and attenuation exponents related to the terrain and geological conditions between the blast source and the protected object, f1(t) and f2(t) are the hole mouth single-segment blasting waveform function and the hole bottom single-segment blasting waveform function respectively.
6. The method of claim 1, wherein, Based on different in-hole delay time intervals, the hole mouth reconstructed single-segment waveforms and the hole bottom reconstructed single-segment waveforms are superimposed to obtain complete single-hole waveforms, and optimal in-hole delay time intervals are determined, comprising: The hole mouth reconstructed single-segment waveforms and the hole bottom reconstructed single-segment waveforms are superimposed according to different in-hole delay time intervals to obtain complete single-hole waveforms of a single blast hole; The in-hole delay time interval corresponding to the complete single-hole waveform with the minimum blasting vibration peak value is determined as the optimal in-hole delay time interval; wherein the complete single-hole waveform comprises: where f hole (t) is the complete single-hole waveform, f 1Q (t) and f 2Q (t) are the orifice-reconstructed single-segment waveform and the bottom-reconstructed single-segment waveform, respectively, and Δt0is the time delay in the bore.
7. The method of claim 6, wherein, Based on different inter-hole delay time intervals, the complete single-hole waveforms are superimposed to obtain multi-hole waveforms, and optimal inter-hole delay time intervals are determined, comprising: Superimposing the complete single-hole waveforms based on an Anderson linear superposition method to obtain multi-hole in-hole segmented blasting waveform functions corresponding to different inter-hole delay time intervals; The inter-hole delay time interval corresponding to the multi-hole in-hole segmented blasting waveform function with the minimum blasting vibration peak value is determined as the optimal inter-hole delay time interval; wherein the multi-hole in-hole segmented blasting waveform function comprises: In the formula, F holes (t, Δti) is the multi-hole in-hole segmented blasting waveform function, Δti is the inter-hole delay interval, i is the number of superimposed waveforms, and m is the number of superimposed waveforms.