A safe blasting construction method

Abstract

The application discloses a safe blasting construction method, comprising the following steps: S1, partitioning a blasting area 5; S2, setting a blast hole parameter and a blasting parameter according to an engineering environment of a corresponding area; S3, establishing a resonance time optimization configuration model, and solving the resonance time optimization configuration model based on a genetic algorithm to obtain minimum time for resonance of a blasting hole; and S4, setting a hole delay blasting time of the corresponding area according to the obtained minimum resonance time of each blasting hole. The resonance time optimization configuration model is established, the resonance time optimization configuration model is solved based on the genetic algorithm to obtain the minimum time for resonance of the blasting hole, under the condition of determining a to-be-solved individual, the hole delay blasting time of each area is set to be greater than the minimum resonance time obtained through the genetic algorithm, so that the total blasting time is reduced under the premise that resonance is not generated in each area during blasting, and the noise pollution time to the engineering periphery is reduced.

Description

Technical Field

[0001] This invention relates to the field of engineering construction, and specifically to a safe blasting construction method. Background Technology

[0002] A cofferdam is a temporary retaining structure built during the construction of hydraulic engineering projects. Its main purpose is to prevent water and soil from entering the construction site of structures, facilitating drainage, excavation of foundation pits, and construction of buildings within the cofferdam. Cofferdams are generally dismantled after use, but sometimes they become part of permanent engineering projects.

[0003] Cofferdam demolition typically employs blasting, which carries risks such as rock collapse, impacting slope stability, and damaging underground structures. Especially when the demolition area is large and multiple detonation points are required, there is a risk of resonance during blasting, potentially damaging underground works and surrounding buildings. To mitigate resonance, inter-hole delayed blasting is generally used. Currently, the inter-hole delayed blasting time is determined based on engineering experience, which may result in prolonged blasting time and extended noise pollution around the project site. Summary of the Invention

[0004] The technical problem to be solved by this invention is that the inter-hole delay blasting time determined in current multi-point blasting demolition projects is based on empirical values, which leads to long blasting times and prolonged noise pollution around the project. The purpose is to provide a safe blasting construction method that reduces the total blasting time and the noise pollution time around the project without causing resonance in each area during blasting.

[0005] This invention is achieved through the following technical solution:

[0006] A safe blasting construction method includes the following steps:

[0007] S1 divides the blasting area into multiple zones;

[0008] S2, set the borehole parameters and blasting parameters according to the engineering environment of the corresponding area;

[0009] S3. Establish a resonance time optimization configuration model, and use a genetic algorithm to solve the resonance time optimization configuration model to obtain the minimum time for the blast hole to generate resonance;

[0010] S4, set the inter-hole delay blasting time for the corresponding area based on the minimum resonance time generated by each blasting hole.

[0011] The beneficial effects of this invention are as follows: First, the blasting area is divided into zones, and then the borehole parameters and blasting parameters are set according to the engineering environment of the corresponding zone. This facilitates the accurate setting of the corresponding borehole parameters and blasting parameters to achieve the ideal demolition effect. Furthermore, a resonance time optimization configuration model is established, and the minimum time for the blasting hole to generate resonance is obtained by solving the resonance time optimization configuration model based on a genetic algorithm. Then, the individual values ​​of each individual to be solved based on the minimum time for the blasting hole to generate resonance are obtained. Under the premise of the determined individuals to be solved, the inter-hole delay blasting time of each zone is set to be greater than the minimum resonance time obtained by the genetic algorithm. This reduces the total blasting time and the noise pollution time to the surrounding area of ​​the project without generating resonance in each zone during blasting.

[0012] In some embodiments, the blasting area in step S1 is divided into a pre-splitting blasting area, a buffer blasting area, and a main blasting area. The pre-splitting blasting area is located at the boundary between the cofferdam demolition section and the cofferdam retention section. The buffer blasting area is located at the cofferdam slope. The main blasting area is all areas of the blasting area other than the pre-splitting blasting area and the buffer blasting area. Dividing the blasting area into pre-splitting blasting area, buffer blasting area, and main blasting area facilitates adaptation to different blasting environments, achieves ideal demolition results, and reduces demolition costs.

[0013] In some embodiments, the borehole parameters in step S2 include borehole inclination angle and borehole depth. The borehole inclination angle in the pre-splitting blasting zone is the same as the corresponding slope ratio, and the borehole depth is shallower than the corresponding slope step depth. The borehole inclination angle in the buffer blasting zone is the same as the corresponding slope ratio, and the borehole depth is shallower than the corresponding slope step depth. The boreholes in the main blasting zone are perpendicular to the surrounding rock, and the borehole depth is shallower than the corresponding slope step depth. By setting the borehole inclination angle and borehole depth for the pre-splitting blasting zone, buffer blasting zone, and main blasting zone respectively, different blasting environments can be accommodated.

[0014] In some embodiments, the borehole parameters in step S2 further include the borehole spacing W, which is obtained by a first calculation formula, the first calculation formula being:

[0015] W = (20~40)·d;

[0016] Where W represents the minimum resistance line (1.4m to 2.8m) and d represents the diameter of the explosive charge. The spacing between boreholes in different blasting zones is obtained by setting the first calculation formula, which facilitates better demolition results during blasting.

[0017] In some embodiments, the blasting hole parameters in step S2 further include the borehole spacing α, which is obtained by a second calculation formula, the second calculation formula being:

[0018] α = (7~12)·D;

[0019] Where α is the borehole spacing and D is the borehole diameter. By setting a second calculation formula, the borehole spacing in different blasting zones can be obtained, which facilitates better blasting demolition results.

[0020] In some embodiments, the blasting parameters in step S2 include the explosion point vibration velocity V, the maximum single-explosive charge Q, the shortest distance R from the center of the explosive charge to the building (structure), the terrain coefficient K between the blast point and the protected object, and the geological condition coefficient s between the blast point and the protected object, satisfying a first equation. The first equation is:

[0021] V = K(Q) 1 / 3 / R) s By setting the first equation, the vibration velocity V at the explosion point, the maximum single-explosive charge Q, and the shortest distance R from the center of the explosive charge to the building (structure) can satisfy the first equation, thus achieving the explosion effect while ensuring explosion safety.

[0022] In some embodiments, the shortest distance R from the center of the medicine package to the building (structure) is obtained by a third calculation formula, which is:

[0023]

[0024] Where K is a coefficient of 1 to 1.2 related to the charging path and the degree of blasting, and Q is the maximum single-explosive charge. A third calculation formula is used to obtain the shortest distance from the center of the explosive charge to a building (structure) to ensure the safety of buildings around the blast point.

[0025] In some embodiments, the blasting parameters of step S2 further include the linear charge density Δline and the ultimate compressive strength R of the surrounding rock, which satisfy the second equation. 压 The second equation is:

[0026] △ 线 =0.42·[R 压 ] 0.5 ·α 0.6 ;

[0027] The borehole spacing α is used to ensure the effectiveness of demolition at the blast points by setting a second equation.

[0028] In some embodiments, the specific process of establishing a resonant time optimization configuration model based on a genetic algorithm is as follows:

[0029] S3.1 Set the number of iterations of the genetic algorithm and the number of individuals to be solved in each generation of the population, wherein each individual to be solved includes the borehole spacing α, the borehole row spacing W, the explosion source frequency ω, the explosion source amplitude A, and the explosion source vibration velocity S.

[0030] S3.2. Establish a resonance time optimization configuration model, which includes an objective function and constraints, wherein the objective function is:

[0031] min(T S );

[0032] Among them, T s The time it takes for the explosion point to resonate;

[0033] S3.3. Encode the borehole spacing α, borehole spacing W, explosion source frequency ω, explosion source amplitude A, and explosion source vibration velocity S to obtain the encoded values ​​of the individual to be solved. Based on the constraints, randomly generate an initial population of several individuals to be solved, and let this initial population be the parent population. Let the individuals to be solved in the parent population be the parent individuals. The encoded value of each parent individual includes the borehole spacing α, borehole spacing W, explosion source frequency ω, explosion source amplitude A, and explosion source vibration velocity S.

[0034] S3.4. Input the encoded values ​​of each individual in the parent population into the simulation system to perform blasting simulation and obtain the actual resonance time T at the blasting point. s The actual resonance time T s This refers to the simulation results of the explosion simulation.

[0035] S3.5. Calculate the objective function value of each individual in the parent population based on each simulation result, and then calculate and sort the fitness value of each individual in the parent population based on the fitness function.

[0036] S3.6. Save the top M parent individuals with the largest fitness values ​​in the parent population. Select parent individuals from all parent individuals other than the top M parent individuals with the largest fitness values ​​by roulette wheel and perform crossover and mutation operations to obtain offspring individuals. Calculate the fitness values ​​of the offspring individuals after crossover and mutation and sort them. Reinsert the offspring individuals into the parent population according to their fitness values. Select a set number of individuals to be solved to form a new parent population, and then return to S4.

[0037] S3.7 Repeat S3.4-S3.6 until the number of iterations is reached or the resonance time difference between two adjacent iterations is less than the set threshold.

[0038] In some embodiments, the fitness function is

[0039] Compared with the prior art, the present invention has the following advantages and beneficial effects:

[0040] 1. Divide the blasting area into zones, and then set the borehole parameters and blasting parameters according to the engineering environment of the corresponding zone. This will enable accurate setting of the corresponding borehole parameters and blasting parameters, and achieve the desired demolition effect while ensuring blasting safety.

[0041] 2. Establish a resonance time optimization configuration model, and use a genetic algorithm to solve the resonance time optimization configuration model to obtain the minimum time for resonance to occur in the blasting hole. Then, obtain the individual values ​​of each individual to be solved based on the minimum time for resonance to occur in the blasting hole. Under the condition of the determined individuals to be solved, the inter-hole delay blasting time of each region can be set to be greater than the minimum resonance time obtained by the genetic algorithm. This reduces the total blasting time and the noise pollution time to the surrounding area of ​​the project without resonance occurring in each region during blasting. Attached Figure Description

[0042] To more clearly illustrate the technical solutions of the exemplary embodiments of the present invention, the accompanying drawings used in the embodiments will be briefly described below. It should be understood that the following drawings only show some embodiments of the present invention and should not be considered as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort. In the drawings:

[0043] Figure 1 This is a diagram showing the borehole structure corresponding to each blasting zone in this invention.

[0044] Main blasting hole 1, auxiliary hole 2, buffer blasting hole 3, pre-splitting blasting hole 4, blasting zone 5. Detailed Implementation

[0045] To make the objectives, technical solutions, and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the embodiments and accompanying drawings. The illustrative embodiments and descriptions of the present invention are only used to explain the present invention and are not intended to limit the present invention.

[0046] Throughout this specification, references to "an embodiment," "an example," or "an example" mean that a particular feature, structure, or characteristic described in connection with that embodiment or example is included in at least one embodiment of the present invention. Therefore, the phrases "an embodiment," "an example," "an example," or "an example" appearing in various places throughout the specification do not necessarily refer to the same embodiment or example. Furthermore, specific features, structures, or characteristics can be combined in one or more embodiments or examples in any suitable combination and / or sub-combination. Moreover, those skilled in the art will understand that the illustrations provided herein are for illustrative purposes and are not necessarily drawn to scale. The term "and / or" as used herein includes any and all combinations of one or more of the associated listed items.

[0047] In the description of this invention, the terms "front", "rear", "left", "right", "up", "down", "vertical", "horizontal", "high", "low", "inner", and "outer" indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings. They are only for the convenience of describing this invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limiting the scope of protection of this invention.

[0048] The terms "first," "second," etc., used in this invention are merely for clarity of description and are not intended to limit any order or emphasize importance. Furthermore, the term "connection" as used herein, unless otherwise specified, can refer to a direct connection or an indirect connection via other components.

[0049] Example

[0050] A safe blasting construction method is provided, comprising the following steps:

[0051] S1, divide the blasting zone 5 into sections;

[0052] S2, set the borehole parameters and blasting parameters according to the engineering environment of the corresponding area;

[0053] S3. Establish a resonance time optimization configuration model, and use a genetic algorithm to solve the resonance time optimization configuration model to obtain the minimum time for the blast hole to generate resonance;

[0054] S4. Based on the minimum resonance time obtained from each blasting hole, the inter-hole delay blasting time for the corresponding area is set. First, the blasting area 5 is divided into zones, and then the blasting hole parameters and blasting parameters are set according to the engineering environment of the corresponding area. This facilitates the accurate setting of the corresponding blasting hole parameters and blasting parameters to obtain the ideal demolition effect. A resonance time optimization configuration model is also established, and the minimum time for the blasting hole to generate resonance is obtained by solving the resonance time optimization configuration model based on the genetic algorithm. Then, the individual values ​​of each to be solved based on the minimum time for the blasting hole to generate resonance are obtained. Under the condition of the determined individuals to be solved, the inter-hole delay blasting time of each area is set to be greater than the minimum resonance time obtained by the genetic algorithm. This reduces the total blasting time and the noise pollution time to the surrounding area of ​​the project without generating resonance during blasting.

[0055] Specifically, the blasting area in step S1 is divided into a pre-splitting blasting area, a buffer blasting area, and a main blasting area. The pre-splitting blasting area is located at the boundary between the cofferdam being demolished and the cofferdam being retained. The buffer blasting area is located on the cofferdam slope. The main blasting area includes all areas of the blasting area except for the pre-splitting and buffer blasting areas. Dividing the blasting area into pre-splitting, buffer, and main blasting areas facilitates adaptation to different blasting environments, achieving ideal demolition results and reducing demolition costs.

[0056] See Figure 1 The borehole parameters in step S2 include the borehole inclination angle and borehole depth. In the pre-splitting blasting zone, the inclination angle of the pre-splitting blasting hole 4 is the same as the corresponding slope ratio, and the depth of the pre-splitting blasting hole 4 is shallower than the corresponding slope step depth. In the buffer blasting zone, the inclination angle of the buffer blasting hole 3 is the same as the corresponding slope ratio, and the depth of the buffer blasting hole 3 is shallower than the corresponding slope step depth. In the main blasting zone, the main blasting hole 1 is perpendicular to the surrounding rock, and the depth of the main blasting hole 1 is shallower than the corresponding slope step depth. By setting the borehole inclination angle and borehole depth separately for the pre-splitting blasting zone, buffer blasting zone, and main blasting zone, different blasting environments can be accommodated.

[0057] See Figure 1 The borehole parameters in step S2 also include the borehole spacing W, which is obtained by a first calculation formula, namely:

[0058] W = (20~40)·d;

[0059] Where W represents the minimum resistance line (1.4m to 2.8m) and d represents the diameter of the explosive charge. The spacing between boreholes in different blasting zones is obtained by setting the first calculation formula, which facilitates better demolition results during blasting.

[0060] The blasting hole parameters in step S2 also include the hole spacing α, which is obtained by a second calculation formula, namely:

[0061] α = (7~12)·D;

[0062] Where α is the borehole spacing and D is the borehole diameter. By setting a second calculation formula, the borehole spacing in different blasting zones can be obtained, which facilitates better blasting demolition results.

[0063] The blasting parameters in step S2 include the explosion point vibration velocity V, the maximum single-explosive charge Q, the shortest distance R from the center of the explosive charge to the building (structure), the terrain coefficient K between the blast point and the protected object, and the geological condition coefficient s between the blast point and the protected object, satisfying the first equation:

[0064] V = K(Q)1 / 3 / R) s By setting the first equation, the vibration velocity V at the explosion point, the maximum single-explosive charge Q, and the shortest distance R from the center of the explosive charge to the building (structure) can all satisfy the first equation, thus ensuring the explosion effect while ensuring explosion safety. The maximum single-explosive charge Q is determined according to the site design requirements.

[0065] The shortest distance R from the center of the medicine package to the building (structure) is obtained by a third calculation formula, which is:

[0066]

[0067] Where K is a coefficient of 1 to 1.2 related to the charging path and blasting intensity, and Q is the maximum single-explosive charge. A third calculation formula is used to obtain the shortest distance from the center of the explosive charge to a building or structure to ensure the safety of buildings surrounding the blast point. The maximum single-explosive charge Q is determined according to the site design requirements.

[0068] Specifically, the blasting parameters in step S2 also include the linear charge density Δline and the ultimate compressive strength R of the surrounding rock, which satisfy the second equation. 压 The second equation is:

[0069] △ 线 =0.42·[R 压 ] 0.5 ·α 0.6 ;

[0070] The borehole spacing α is used to ensure the effectiveness of demolition at the blast points by setting a second equation.

[0071] Specifically, the process of establishing a resonance time optimization configuration model based on a genetic algorithm is as follows: S3.1, setting the number of iterations of the genetic algorithm and the number of individuals to be solved in each generation of the population, wherein each individual to be solved includes the borehole spacing α, the borehole row spacing W, the explosion source frequency ω, the explosion source amplitude A, and the explosion source vibration velocity S.

[0072] Since the borehole spacing α, the borehole row spacing W, the explosion source frequency ω, the explosion source amplitude A, and the explosion source vibration velocity S can all affect the time for resonance to occur during an explosion, determining the borehole spacing, borehole row spacing, explosion source frequency, explosion source amplitude, and explosion source vibration velocity is an important step when the blasting area is determined.

[0073] S3.2. Establish a resonance time optimization configuration model, which includes an objective function and constraints, wherein the objective function is:

[0074] min(T S );

[0075] Among them, T s The time it takes for the explosion point to resonate;

[0076] S3.3. Encode the borehole spacing α, borehole spacing W, explosion source frequency ω, explosion source amplitude A, and explosion source vibration velocity S to obtain the encoded values ​​of the individual to be solved. Based on the constraints, randomly generate an initial population of several individuals to be solved, and let this initial population be the parent population. Let the individuals to be solved in the parent population be the parent individuals. The encoded value of each parent individual includes the borehole spacing α, borehole spacing W, explosion source frequency ω, explosion source amplitude A, and explosion source vibration velocity S.

[0077] S3.4. Input the encoded values ​​of each individual in the parent population into the simulation system to perform blasting simulation and obtain the actual resonance time T at the blasting point. s The actual resonance time T s This refers to the simulation results of the explosion simulation; the coded values ​​of each individual are obtained based on the design range of borehole parameters and blasting parameters.

[0078] S3.5. Calculate the objective function value of each individual in the parent population based on each simulation result, and then calculate and sort the fitness value of each individual in the parent population based on the fitness function.

[0079] S3.6. Save the top M parent individuals with the largest fitness values ​​in the parent population. Select parent individuals from all parent individuals other than the top M parent individuals with the largest fitness values ​​by roulette wheel and perform crossover and mutation operations to obtain offspring individuals. Calculate the fitness values ​​of the offspring individuals after crossover and mutation and sort them. Reinsert the offspring individuals into the parent population according to their fitness values. Select a set number of individuals to be solved to form a new parent population, and then return to S4.

[0080] S3.7 Repeat S3.4-S3.6 until the number of iterations is reached or the resonance time difference between two adjacent iterations is less than the set threshold.

[0081] The threshold is set between 1 and 2 seconds. Simulation stops when the resonance time difference between two adjacent iterations is less than this threshold. This allows the determination of the minimum resonance time generated by the blasting during the simulation. The values ​​corresponding to the individual parameters to be solved (including borehole spacing α, borehole spacing W, explosion source frequency ω, explosion source amplitude A, and explosion source vibration velocity S) based on the minimum resonance time can then be determined. This allows for the determination of the values ​​for each individual parameter (facilitating the setting of corresponding borehole and blasting parameters based on the determined individual parameter values). Under the condition of the determined individual parameter values ​​(including borehole spacing α, borehole spacing W, explosion source frequency ω, explosion source amplitude A, and explosion source vibration velocity S), the inter-hole delay blasting time for each region is set to be greater than the minimum resonance time obtained through the genetic algorithm (i.e., obtaining the inter-hole delay blasting time for the corresponding blasting region). This reduces the total blasting time and noise pollution time around the project without causing resonance in each region during blasting.

[0082] Wherein, the fitness function is

[0083] That is, if The smaller the value, the larger the fitness function F(x), and the more suitable the encoding value of the corresponding parent individual is.

[0084] The specific embodiments described above further illustrate the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above description is only a specific embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A safe blasting construction method, characterized in that, Includes the following steps: S1 divides the blasting area into multiple zones; S2, set the borehole parameters and blasting parameters according to the engineering environment of the corresponding area; S3. Establish a resonance time optimization configuration model, and use a genetic algorithm to solve the resonance time optimization configuration model to obtain the minimum time for the blast hole to generate resonance; S4, set the inter-hole delay blasting time for the corresponding area based on the minimum resonance time generated by each blasting hole; The borehole parameters in step S2 include the borehole spacing W, which is obtained by a first calculation formula, which is: W = (20~40)·d; Where W is the minimum resistance line of 1.4m to 2.8m, and d is the diameter of the explosive charge; the blasting hole parameters in step S2 also include the borehole spacing α, which is obtained by a second calculation formula, the second calculation formula being: α＝(7～12)·D； Where α is the borehole spacing and D is the borehole diameter; The blasting parameters in step S2 include the vibration velocity V at the blast point, the maximum single-explosive charge Q, the shortest distance R from the center of the explosive charge to the building, the topographic coefficient K between the blast point and the protected object, and the geological condition coefficient s between the blast point and the protected object, satisfying the first equation: V=K(Q 1 / 3 / R) s ; The shortest distance R from the center of the medicine pack to the building is obtained by a third calculation formula, which is: Where K is a coefficient of 1 to 1.2 related to the charging route and the degree of blasting, and Q is the maximum single-shot charge; The blasting parameters in step S2 also include the linear charge density Δline and the ultimate compressive strength R of the surrounding rock, which satisfy the second equation. 压 The second equation is: △ 线 ＝0.42·[R 压 ] 0.5 ·α 0.6 ； Where α is the borehole spacing.

2. The safe blasting construction method according to claim 1, characterized in that, The blasting area in step S1 is divided into a pre-splitting blasting area, a buffer blasting area, and a main blasting area. The pre-splitting blasting area is located at the boundary between the demolished part and the retained part of the cofferdam. The buffer blasting area is located at the slope of the cofferdam. The main blasting area is the area other than the pre-splitting blasting area and the buffer blasting area.

3. The safe blasting construction method according to claim 2, characterized in that, The borehole parameters in step S2 include the borehole inclination angle and borehole depth. The borehole inclination angle in the pre-splitting blasting zone is the same as the corresponding slope ratio, and the borehole depth is shallower than the corresponding slope step depth. The borehole inclination angle in the buffer blasting zone is the same as the corresponding slope ratio, and the borehole depth is shallower than the corresponding slope step depth. The boreholes in the main blasting zone are perpendicular to the surrounding rock, and the borehole depth is shallower than the corresponding slope step depth.

4. The safe blasting construction method according to claim 1, characterized in that, The specific process of establishing a resonance time optimization configuration model based on genetic algorithm is as follows: S3.1 Set the number of iterations of the genetic algorithm and the number of individuals to be solved in each generation of the population, wherein each individual to be solved includes the borehole spacing α, the borehole row spacing W, the explosion source frequency ω, the explosion source amplitude A, and the explosion source vibration velocity S. S3.

2. Establish a resonance time optimization configuration model, which includes an objective function and constraints, wherein the objective function is: min(T S ); Among them, T s The time it takes for the explosion point to resonate; S3.

3. Encode the borehole spacing α, borehole spacing W, explosion source frequency ω, explosion source amplitude A, and explosion source vibration velocity S to obtain the encoded values ​​of the individual to be solved. Based on the constraints, randomly generate an initial population of several individuals to be solved, and let this initial population be the parent population. Let the individuals to be solved in the parent population be the parent individuals. The encoded value of each parent individual includes the borehole spacing α, borehole spacing W, explosion source frequency ω, explosion source amplitude A, and explosion source vibration velocity S. S3.

4. Input the encoded values ​​of each individual in the parent population into the simulation system to perform blasting simulation and obtain the actual resonance time T at the blasting point. s The actual resonance time T s This refers to the simulation results of the explosion simulation. S3.

5. Calculate the objective function value of each individual in the parent population based on each simulation result, and then calculate and sort the fitness value of each individual in the parent population based on the fitness function. S3.

6. Save the top M parent individuals with the largest fitness values ​​in the parent population. Select parent individuals from all parent individuals other than the top M parent individuals with the largest fitness values ​​by roulette wheel and perform crossover and mutation operations to obtain offspring individuals. Calculate the fitness values ​​of the offspring individuals after crossover and mutation and sort them. Reinsert the offspring individuals into the parent population according to their fitness values. Select a set number of individuals to be solved to form a new parent population, and then return to S4. S3.7 Repeat S3.4-S3.6 until the number of iterations is reached or the resonance time difference between two adjacent iterations is less than the set threshold.

5. The safe blasting construction method according to claim 4, characterized in that, The fitness function is:

Images