CN116882213A - Method and system for calculating continuous detonation delay time of electronic detonator - Google Patents
Method and system for calculating continuous detonation delay time of electronic detonator Download PDFInfo
- Publication number
- CN116882213A CN116882213A CN202311145957.0A CN202311145957A CN116882213A CN 116882213 A CN116882213 A CN 116882213A CN 202311145957 A CN202311145957 A CN 202311145957A CN 116882213 A CN116882213 A CN 116882213A
- Authority
- CN
- China
- Prior art keywords
- delay time
- blasting
- calculating
- vibration
- parameters
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000005474 detonation Methods 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000005422 blasting Methods 0.000 claims abstract description 108
- 238000004364 calculation method Methods 0.000 claims abstract description 48
- 230000014509 gene expression Effects 0.000 claims description 72
- 230000006870 function Effects 0.000 claims description 33
- 238000012545 processing Methods 0.000 claims description 27
- 230000000977 initiatory effect Effects 0.000 claims description 26
- 239000002245 particle Substances 0.000 claims description 10
- 238000004088 simulation Methods 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 8
- 238000009933 burial Methods 0.000 claims description 7
- 239000002360 explosive Substances 0.000 claims description 5
- 238000010168 coupling process Methods 0.000 claims description 4
- 230000003993 interaction Effects 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 239000007787 solid Substances 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 239000011435 rock Substances 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 7
- 238000004891 communication Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000005236 sound signal Effects 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 238000009412 basement excavation Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- KLDZYURQCUYZBL-UHFFFAOYSA-N 2-[3-[(2-hydroxyphenyl)methylideneamino]propyliminomethyl]phenol Chemical compound OC1=CC=CC=C1C=NCCCN=CC1=CC=CC=C1O KLDZYURQCUYZBL-UHFFFAOYSA-N 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 201000001098 delayed sleep phase syndrome Diseases 0.000 description 1
- 208000033921 delayed sleep phase type circadian rhythm sleep disease Diseases 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004141 dimensional analysis Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Computational Mathematics (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Drilling And Exploitation, And Mining Machines And Methods (AREA)
Abstract
The application relates to the technical field of tunnel blasting, and relates to a method and a system for calculating continuous detonation delay time of an electronic detonator, wherein the method comprises the steps of obtaining modeling parameters and a first function, wherein the first function is used for calculating vibration speeds transmitted to a target at blastholes under different blasting delay times, and the blasting delay time comprises blasting delay time between a current blasthole and a next blasthole; establishing at least one numerical model according to the modeling parameters, wherein the numerical model is used for simulating and calculating the vibration speed of the single blasthole detonation transferred to the target; calculating vibration speeds transmitted to a target position from the blasting position under different blasting delay time by using at least one numerical model to obtain vibration speed information; fitting the first function according to the vibration speed information to obtain a vibration speed calculation formula; the application provides a calculation basis for the blasting delay time between the current blasthole and the next blasthole when the electronic detonator single-hole continuous detonation technology is adopted.
Description
Technical Field
The application relates to the technical field of tunnel blasting, in particular to a method and a system for calculating continuous detonation delay time of an electronic detonator.
Background
With the continuous development of the urban process, the construction requirements of urban subways, highways and railway tunnels are continuously increased. In the face of the complex environment of cities, the tunnel blasting excavation mostly adopts an electronic detonator single-hole continuous detonation technology to control blasting vibration, so that the influence on the environment is reduced. At present, no mature calculation method is formed for the continuous detonation delay time of the electronic detonator. If the delay interval of detonation between the blast holes is too large, although vibration can be effectively reduced, the blasting effect is poor, especially the cut hole blasting is easy to cause cut failure and the blast hole utilization rate is low; the smaller the detonation delay interval between the blast holes is, the larger the vibration is, and the damage to surrounding building (construction) is more likely to occur, so that a method for calculating the continuous detonation delay time of the electronic detonator is needed to select a proper detonation delay time when the single-hole continuous detonation technology of the electronic detonator is adopted.
Disclosure of Invention
The application aims to provide a method and a system for calculating continuous detonation delay time of an electronic detonator so as to solve the problems.
In order to achieve the above object, the embodiment of the present application provides the following technical solutions:
in one aspect, the embodiment of the application provides a method for calculating continuous detonation delay time of an electronic detonator, which comprises the following steps:
obtaining modeling parameters and a first function, wherein the modeling parameters comprise tunnel section size parameters, tunnel burial depth parameters, blast hole position parameters, blast hole loading quantity parameters and blast hole depth parameters, the first function is used for calculating the vibration speed transmitted to a target position from a blast hole under different blasting delay time, and the blasting delay time comprises the blasting delay time between a current blast hole and a next blast hole;
establishing at least one numerical model according to the modeling parameters, wherein the numerical model is used for simulating and calculating the vibration speed of single blasthole detonation transferred to a target;
calculating vibration speeds transmitted to a target position from the blasting position under different blasting delay time by using at least one numerical model to obtain vibration speed information;
fitting the first function according to the vibration speed information to obtain a vibration speed calculation formula;
and calculating the blasting delay time between the current blast hole and the next blast hole according to the vibration speed calculation formula.
In a second aspect, an embodiment of the present application provides a system for calculating a continuous initiation delay time of an electronic detonator, where the system includes:
the system comprises an acquisition module, a control module and a control module, wherein the acquisition module is used for acquiring modeling parameters and a first function, the modeling parameters comprise tunnel section size parameters, tunnel burial depth parameters, blast hole position parameters, blast hole charge quantity parameters and blast hole depth parameters, the first function is used for calculating vibration speeds transmitted to a target from a blast hole under different blasting delay time, and the blasting delay time comprises blasting delay time between a current blast hole and a next blast hole;
the establishing module is used for establishing at least one numerical model according to the modeling parameters, and the numerical model is used for simulating and calculating the vibration speed of the single blasthole detonation transferred to the target;
the simulation module is used for calculating the vibration speed transmitted to the target position from the blasting position under different blasting delay time by using at least one numerical model to obtain vibration speed information;
the first processing module is used for fitting the first function according to the vibration speed information to obtain a vibration speed calculation formula;
and the second processing module is used for calculating the blasting delay time between the current blast hole and the next blast hole according to the vibration speed calculation formula.
In a third aspect, an embodiment of the present application provides an electronic detonator continuous initiation delay time calculation device, where the device includes a memory and a processor. The memory is used for storing a computer program; and the processor is used for realizing the step of the method for calculating the continuous detonation delay time of the electronic detonator when executing the computer program.
In a fourth aspect, an embodiment of the present application provides a readable storage medium, where a computer program is stored, where the computer program, when executed by a processor, implements the steps of the method for calculating a continuous detonation delay time of an electronic detonator described above.
The beneficial effects of the application are as follows:
according to the application, the vibration speed transmitted to the target position from the blasting position under different blasting delay time is simulated by establishing at least one numerical model, so that vibration speed information is obtained, then a first function is fitted according to the vibration speed information and the blasting vibration period, so that a vibration speed calculation formula is obtained, the blasting delay time between the current blasthole and the next blasthole can be reversely deduced through the vibration speed calculation formula, the problem that the Sadawski formula in the prior art is only suitable for vibration prediction of same-segment blasting is solved, the problem of damping blasting vibration under the condition of delay time blasting between electronic detonator holes is not considered, and a calculation basis is provided for the delay time of blasting between the current blasthole and the next blasthole when the electronic detonator single-hole continuous blasting technology is adopted.
Additional features and advantages of the application will be set forth in the description which follows, and in part will be apparent from the description, or may be learned by practice of the embodiments of the application. The objectives and other advantages of the application will be realized and attained by the structure particularly pointed out in the written description and claims thereof as well as the appended drawings.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the embodiments will be briefly described below, it being understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and other related drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a schematic flow chart of a method for calculating continuous initiation delay time of an electronic detonator according to an embodiment of the application.
Fig. 2 is a schematic structural diagram of a system for calculating continuous initiation delay time of an electronic detonator according to an embodiment of the application.
Fig. 3 is a schematic structural diagram of a continuous initiation delay time calculation device for an electronic detonator according to an embodiment of the application.
Fig. 4 is a schematic diagram of a continuous initiation mode of an electronic detonator according to an embodiment of the application.
The drawing is marked: 901. an acquisition module; 902. establishing a module; 903. a simulation module; 904. a first processing module; 905. a second processing module; 9011. the system comprises a first acquisition unit, a 9012 query unit; 9013. a first calculation unit; 90131. a second calculation unit; 90132. a third calculation unit; 90133. a fourth calculation unit; 901321, fifth calculation unit; 901322, sixth computing unit; 901323, seventh calculation unit; 9021. a building unit; 9022. dividing units; 9023. a first processing unit; 9024. a second processing unit; 9041. a second acquisition unit; 9042. a third processing unit; 9043. a fourth processing unit; 800. the electronic detonator continuous detonation delay time calculating device; 801. a processor; 802. a memory; 803. a multimedia component; 804. an I/O interface; 805. a communication component; 11. an auxiliary hole; 12. and (5) cutting the slot.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present application more apparent, the technical solutions of the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application, and it is apparent that the described embodiments are some embodiments of the present application, but not all embodiments of the present application. The components of the embodiments of the present application generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the application, as presented in the figures, is not intended to limit the scope of the application, as claimed, but is merely representative of selected embodiments of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
It should be noted that: like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further definition or explanation thereof is necessary in the following figures. Meanwhile, in the description of the present application, the terms "first", "second", and the like are used only to distinguish the description, and are not to be construed as indicating or implying relative importance.
Example 1:
the embodiment provides a method for calculating continuous detonation delay time of an electronic detonator, and it can be understood that a scene, for example, a scene of blasting tunnel cutting, can be paved in the embodiment.
Referring to fig. 1, the method is shown to include steps S1, S2, S3, S4, and S5.
Step S1, obtaining modeling parameters and a first function, wherein the modeling parameters comprise tunnel section size parameters, tunnel burial depth parameters, blast hole position parameters, blast hole charge quantity parameters and blast hole depth parameters, the first function is used for calculating vibration speeds transmitted to a target position from a blast hole under different blasting delay time, and the blasting delay time comprises blasting delay time between a current blast hole and a next blast hole;
it can be understood that the tunnel section size parameter, the tunnel burial depth parameter, the blasthole position parameter, the blasthole charge parameter and the blasthole depth parameter are all actual parameters measured in an actual tunnel blasting scene.
It can be understood that the step S1 further includes a step S11, a step S12, and a step S13, where specific details are:
step S11, acquiring first information, wherein the first information comprises physical parameters affecting blasting vibration propagation;
it is understood that the physical parameters affecting the propagation of the blast vibration include the dosage, the distance between the blasts, the velocity of the wave, the detonation time, the period of the blast vibration, the particle vibration velocity, the vibration displacement, the frequency, the particle acceleration, the medium density, and the blast delay time between the current blast hole and the next blast hole.
Step S12, inquiring in a preset parameter library to obtain a dimension expression of each physical parameter;
in this step, the dimension expression of each physical parameter is obtained by searching in a preset parameter library as shown in table 1:
TABLE 1 influence the dimensions of the blasted vibration corresponding to each physical parameter propagated
And step S13, calculating according to the dimensional expression of each physical parameter to obtain a first function.
It can be understood that the step S13 further includes a step S131, a step S132, and a step S133, where specifically:
step S131, calculating the dimensionless expression of each physical parameter according to the dimensionless principles to obtain second information, wherein the second information comprises a dimensionless expression corresponding to a blasting vibration period, a dimensionless expression corresponding to a vibration wave speed, a dimensionless expression corresponding to a blasting delay time, a dimensionless expression corresponding to a particle vibration speed and a dimensionless expression corresponding to a medium density;
it is understood that the dimensionless expression corresponding to the blasting vibration period in this step is specifically:
as can be seen from table 1, among them,converting it into dimensionless form pi 1 The representation is:
in the above formula, T is detonation time, T is blasting vibration period, and the dimensionless expression corresponding to the vibration wave speed can be obtained by the same principle according to the dimensionless uniformity, wherein the method specifically comprises the following steps:
in the above formula, c is the vibration wave velocity, t is the detonation time, and R is the detonation distance.
Similarly, the dimensionless expression corresponding to the blasting delay time is represented by pi 3 The representation is specifically as follows:
in the above-mentioned method, the step of,blasting delay time.
Similarly, the dimensionless expression corresponding to the vibration speed of the particles is represented by pi 4 The representation is specifically as follows:
in the above equation, V is the particle vibration velocity.
Similarly, the dimensionless expression corresponding to the medium density is represented by pi 5 The representation is specifically as follows:
in the above formula, Q is the explosive amount,is the medium density.
Step S132, calculating according to the second information to obtain at least two dimensionless expressions;
it can be understood that the step S132 further includes a step S1321, a step S1322, and a step S1323, where specific details are:
step S1321, calculating according to a dimensionless expression corresponding to the blasting vibration period and a dimensionless expression corresponding to the blasting delay time to obtain a first dimensionless expression;
it can be understood that the dimensionless harmonic theory is applied toAnd dimensionless->Performing calculation, and setting the result as +.>I.e. a first dimensionless expression, wherein +.>The method comprises the following steps:
step S1322, calculating according to the dimensionless expression corresponding to the vibration wave speed and the dimensionless expression corresponding to the particle vibration speed to obtain a second dimensionless expression;
it is understood that the same theory is applied to dimensionless pi 2 And dimensionless pi 4 Performing calculation, and setting the result asI.e. the second dimensionless expression, gives pi 7 The method comprises the following steps:
and step S1323, calculating according to the dimensionless expression corresponding to the medium density and the first dimensionless expression to obtain a third dimensionless expression.
It will be appreciated that the dimensionless expression pi for medium density corresponds to 5 And a first dimensionless expression pi 6 Performing calculation, and setting the result asI.e. the third dimensionless expression, gives pi 8 The method comprises the following steps:
and step S133, calculating according to at least two dimensionless expressions to obtain the first function.
It will be appreciated that the second dimensionless expression is known from the principle of dimensional analysisAnd a third dimensionless expression pi 8 All have dimensionless, vibration wave velocity c and value density +.>Assuming a constant, a first function can be obtainedThe method is characterized in that when the tunnel blasting is carried out, the detonation delay time between the current blasthole and the next blasthole is not too long or too short, when the detonation delay time is too long, the blasting excavation efficiency is low, and adverse effects are easily brought to surrounding buildings when the detonation delay time is too short, so that when the electronic detonator single-hole continuous detonation technology is adopted, a proper detonation delay time is required to be selected, but a Sagnac empirical formula is adopted in the engineering of the current tunnel blasting construction, a basis is provided for the on-site blasting parameter design, the traditional Sagnac formula is only suitable for vibration prediction of simultaneous detonation, the condition that the detonation time between the current blasthole and the next blasthole is different when the electronic detonator single-hole continuous detonation technology is adopted is not considered, and the factor of different detonation time between the current blasthole and the next blasthole when the electronic detonator single-hole continuous detonation technology is adopted is comprehensively considered.
S2, establishing at least one numerical model according to the modeling parameters, wherein the numerical model is used for simulating and calculating the vibration speed of a single blasthole detonation transmitted to a target;
it can be understood that the numerical model is built based on the actually measured parameters as modeling parameters, so that the process of simulating explosion is closer to reality, and the calculation accuracy is improved.
It can be understood that the step S2 further includes a step S21, a step S22, a step S23, and a step S24, where specific details are:
s21, establishing a geometric model of single-hole continuous detonation of the tunnel electronic detonator by using the modeling information;
it can be understood that the establishment of the geometric model for single-hole continuous initiation of the tunnel electronic detonator according to the modeling information by using finite element software is a technical scheme well known to those skilled in the art, and therefore is not described herein.
S22, performing independent mapping division on the blast holes in the geometric model and setting boundary conditions on the numerical model;
it can be understood that the blast holes are individually mapped and divided by Mesh, so that each blast hole is an individual PART, and different blasting delay times are set for each blast hole.
S23, selecting material parameters which accord with field reality;
and S24, establishing a numerical model of the interaction between the rock, the explosive and the air during tunnel blasting by adopting a fluid-solid coupling method.
S3, calculating vibration speeds transmitted to a target position from the blasting position under different blasting delay time by using at least one numerical model to obtain vibration speed information;
in this step, ANSYS/LS-DYNA software is used to create a blasthole n=6, and the blastholes are numbered as blasthole 1, blasthole 2, blasthole 3, blasthole 4, blasthole 5, and blasthole 6, the blastholes 1-6 are blastholes in different rows, and the blasthole 1 is adjacent to the blasthole 2, the blasthole 2 is adjacent to the blasthole 3, the blasthole 3 is adjacent to the blasthole 4, the blasthole 4 is adjacent to the blasthole 5, the blasthole 5 is adjacent to the blasthole 6, and the blasting delay time is set to 5 different blasting delay times of 1ms, 2ms, 3ms, 4ms, 5ms, etc. Defining blast hole detonation time by a keyword of initiation_detonation, setting a model 1 to be single-hole detonation, namely, only blast hole 1 is detonated, changing blast hole detonation time on the basis of the model 1 by a model 2, setting blasting delay time of the blast hole to be 1ms, namely, blast hole 1 is detonated for 0ms, blast hole 2 is detonated for 1ms, and blast hole 3 is detonated for … … ms; meanwhile, the blasting delay time of the models 3, 4, 5 and 6 is respectively 2ms, 3ms, 4ms and 5ms, and the vibration speed transmitted to the target position from the blasting position under different blasting delay time can be simulated through the setting.
S4, fitting the first function according to the vibration speed information to obtain a vibration speed calculation formula;
it can be understood that the step S4 further includes a step S41, a step S42, and a step S43, where specific details are:
s41, performing simulation by using a numerical model to obtain a blasting vibration waveform diagram;
the blast vibration waveform in this step simulates a single hole blast for a numerical model to obtain a blast vibration waveform at a distance of 15m from the blast hole.
Step S42, determining a blasting vibration period according to the blasting vibration waveform diagram;
and step S43, fitting the first function according to the blasting vibration period and the vibration speed information to obtain a vibration speed calculation formula.
It will be appreciated that when the target site, i.e. the measuring point, is determined to be 15m from the blast hole, the following simulation data can be obtained through simulation, as shown in the following table:
the vibration velocity calculation formula obtained by fitting the first function according to the simulation data obtained by the table is specifically:
and S5, calculating blasting delay time between the current blast hole and the next blast hole according to the vibration speed calculation formula.
The optimal blasting delay time between the current blast hole and the next blast hole can be obtained by performing back-pushing according to the vibration velocity calculation formula, wherein the optimal blasting delay time is specifically as follows: according to the requirements of the site actual vibration safety control standard, the maximum surface vibration speed at 15m from the blast holes is 0.7cm/s, the explosive loading quantity of each blast hole is preset to be 0.6KG, and the vibration period at 15m from the blast holes is 10ms, so that the optimal delay time is calculated through the formulaMeanwhile, the minimum time interval of the delay time is set to be 2.53ms, and when the delay time is 2ms, the vibration speed exceeds the field actual vibration safety control standard requirement (0.7 cm/s), so that the delay time between holes can be 3ms to be used as the optimal blasting delay time according to a prediction formula and actual conditions.
The blasting mode of the blast hole is determined according to the type of the blast hole; detonating according to the blasting mode and the optimal blasting delay time, wherein when the type of the blast holes is the cut hole 12, the inner two rows of blast holes are detonated at intervals in a staggered way, and the outer two rows of blast holes are detonated at intervals in a staggered way; when the type of the blast hole is the auxiliary hole 11, after the same row of blast holes are sequentially detonated at intervals, the next row of blast holes are sequentially detonated at intervals, wherein the specific detonation mode is shown as a schematic diagram of the continuous detonation mode of the electronic detonator in fig. 4.
Example 2:
as shown in fig. 2, the present embodiment provides a system for calculating continuous initiation delay time of an electronic detonator, which includes an acquisition module 901, a setup module 902, an analog module 903, a first processing module 904, and a second processing module 905, wherein the system specifically includes:
the obtaining module 901 is configured to obtain modeling parameters and a first function, where the modeling parameters include a tunnel section size parameter, a tunnel burial depth parameter, a blasthole position parameter, a blasthole charge parameter, and a blasthole depth parameter, and the first function is configured to calculate a vibration velocity transferred from a blasthole to a target under different blasting delay times, and the blasting delay times include a blasting delay time between a current blasthole and a next blasthole;
the establishing module 902 is configured to establish at least one numerical model according to the modeling parameters, where the numerical model is used to simulate and calculate a vibration velocity of a single blasthole detonation transferred to a target;
the simulation module 903 is configured to calculate, using at least one of the numerical models, a vibration velocity transmitted from a blasting location to a target location under different blasting delay times, so as to obtain vibration velocity information;
the first processing module 904 is configured to fit the first function according to the vibration velocity information to obtain a vibration velocity calculation formula;
and a second processing module 905, configured to calculate a blasting delay time between the current blasthole and the next blasthole according to the vibration velocity calculation formula.
In a specific embodiment of the disclosure, the acquiring module 901 further includes a first acquiring unit 9011, a querying unit 9012, and a first calculating unit 9013, where the specific steps are:
a first acquiring unit 9011 for acquiring first information including a physical parameter affecting propagation of blasting vibration;
a query unit 9012, configured to query in a preset parameter library to obtain a dimension expression of each physical parameter;
a first calculating unit 9013, configured to calculate according to the dimensional expression of each physical parameter, to obtain a first function.
In a specific embodiment of the disclosure, the first computing unit 9013 further includes a second computing unit 90131, a third computing unit 90132, and a fourth computing unit 90133, where specifically:
a second calculating unit 90131, configured to calculate a dimensionless expression of each physical parameter according to a dimensionless alignment principle, to obtain second information, where the second information includes a dimensionless expression corresponding to a blasting vibration period, a dimensionless expression corresponding to a vibration wave speed, a dimensionless expression corresponding to a blasting delay time, a dimensionless expression corresponding to a particle vibration speed, and a dimensionless expression corresponding to a medium density;
a third computing unit 90132, configured to perform computation according to the second information to obtain at least two dimensionless expressions;
a fourth calculating unit 90133, configured to calculate according to at least two dimensionless expressions, to obtain the first function.
In a specific embodiment of the disclosure, the third computing unit 90132 further includes a fifth computing unit 901321, a sixth computing unit 901322, and a seventh computing unit 901323, where specifically:
a fifth calculating unit 901321, configured to calculate according to a dimensionless expression corresponding to the blasting vibration period and a dimensionless expression corresponding to the blasting delay time, to obtain a first dimensionless expression;
a sixth calculating unit 901322, configured to calculate according to the dimensionless expression corresponding to the vibration wave velocity and the dimensionless expression corresponding to the particle vibration velocity, to obtain a second dimensionless expression;
and a seventh calculating unit 901323, configured to calculate according to the dimensionless expression corresponding to the medium density and the first dimensionless expression, to obtain a third dimensionless expression.
In a specific embodiment of the disclosure, the establishing module 902 further includes an establishing unit 9021, a dividing unit 9022, a first processing unit 9023, and a second processing unit 9024, where specific details are:
the establishing unit 9021 is used for establishing a geometric model of single-hole continuous detonation of the tunnel electronic detonator by using the modeling information;
a dividing unit 9022, configured to separately map and divide the blastholes in the geometric model and set boundary conditions for the numerical model;
a first processing unit 9023, configured to select a material parameter that meets an actual field;
a second processing unit 9024 is configured to establish a numerical model of the rock-explosive-air interaction during tunnel blasting using a fluid-solid coupling method.
In a specific embodiment of the disclosure, the first processing module 904 further includes a second obtaining unit 9041, a third processing unit 9042, and a fourth processing unit 9043, where specific details are:
a second obtaining unit 9041, configured to obtain a blasting vibration waveform map through simulation by using a numerical model;
a third processing unit 9042, configured to determine a blasting vibration period according to the blasting vibration waveform map;
and a fourth processing unit 9043, configured to fit the first function according to the blasting vibration period and the vibration velocity information, to obtain a vibration velocity calculation formula.
It should be noted that, regarding the system in the above embodiment, the specific manner in which the respective modules perform the operations has been described in detail in the embodiment regarding the method, and will not be described in detail herein.
Example 3:
corresponding to the above method embodiment, an electronic detonator continuous detonation delay time calculating device is further provided in this embodiment, and an electronic detonator continuous detonation delay time calculating device described below and an electronic detonator continuous detonation delay time calculating method described above can be referred to correspondingly with each other.
Fig. 3 is a block diagram illustrating an electronic detonator sequential initiation delay time calculation device 800 in accordance with an exemplary embodiment. As shown in fig. 3, the electronic detonator continuous initiation delay time calculation device 800 may include: a processor 801, a memory 802. The electronic detonator sequential initiation delay time computing device 800 may also include one or more of a multimedia component 803, an i/O interface 804, and a communication component 805.
The processor 801 is configured to control the overall operation of the electronic detonator continuous initiation delay time computing device 800, so as to complete all or part of the steps in the electronic detonator continuous initiation delay time computing method. The memory 802 is used to store various types of data to support operation of the time delay computing device 800 at the electronic detonator continuous initiation, such data may include, for example, instructions for any application or method operating on the time delay computing device 800 at the electronic detonator continuous initiation, as well as application related data such as contact data, messages sent and received, pictures, audio, video, and the like. The Memory 802 may be implemented by any type or combination of volatile or non-volatile Memory devices, such as static random access Memory (Static Random Access Memory, SRAM for short), electrically erasable programmable Read-Only Memory (Electrically Erasable Programmable Read-Only Memory, EEPROM for short), erasable programmable Read-Only Memory (Erasable Programmable Read-Only Memory, EPROM for short), programmable Read-Only Memory (Programmable Read-Only Memory, PROM for short), read-Only Memory (ROM for short), magnetic Memory, flash Memory, magnetic disk, or optical disk. The multimedia component 803 may include a screen and an audio component. Wherein the screen may be, for example, a touch screen, the audio component being for outputting and/or inputting audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signals may be further stored in the memory 802 or transmitted through the communication component 805. The audio assembly further comprises at least one speaker for outputting audio signals. The I/O interface 804 provides an interface between the processor 801 and other interface modules, which may be a keyboard, mouse, buttons, etc. These buttons may be virtual buttons or physical buttons. The communication component 805 is configured to perform wired or wireless communication between the electronic detonator continuous initiation delay time calculation device 800 and other devices. Wireless communication, such as Wi-Fi, bluetooth, near field communication (Near FieldCommunication, NFC for short), 2G, 3G or 4G, or a combination of one or more thereof, the respective communication component 805 may thus comprise: wi-Fi module, bluetooth module, NFC module.
In an exemplary embodiment, the electronic detonator sequential initiation delay time calculation device 800 may be implemented by one or more application specific integrated circuits (Application Specific Integrated Circuit, abbreviated as ASIC), digital signal processors (DigitalSignal Processor, abbreviated as DSP), digital signal processing devices (Digital Signal Processing Device, abbreviated as DSPD), programmable logic devices (Programmable Logic Device, abbreviated as PLD), field programmable gate arrays (Field Programmable Gate Array, abbreviated as FPGA), controllers, microcontrollers, microprocessors, or other electronic components for performing the electronic detonator sequential initiation delay time calculation method described above.
In another exemplary embodiment, a computer readable storage medium is also provided that includes program instructions that, when executed by a processor, implement the steps of the electronic detonator sequential initiation delay time calculation method described above. For example, the computer readable storage medium may be the memory 802 described above including program instructions executable by the processor 801 of the electronic detonator continuous initiation delay time calculation device 800 to perform the electronic detonator continuous initiation delay time calculation method described above.
Example 4:
corresponding to the above method embodiment, a readable storage medium is also provided in this embodiment, and a readable storage medium described below and a method for calculating continuous initiation delay time of an electronic detonator described above may be referred to correspondingly.
A readable storage medium, on which a computer program is stored, which when executed by a processor, implements the steps of the method for calculating the continuous detonation delay time of an electronic detonator according to the above method embodiment.
The readable storage medium may be a usb disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), a magnetic disk, or an optical disk, and the like.
The above description is only of the preferred embodiments of the present application and is not intended to limit the present application, but various modifications and variations can be made to the present application by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the protection scope of the present application.
The foregoing is merely illustrative of the present application, and the present application is not limited thereto, and any person skilled in the art will readily recognize that variations or substitutions are within the scope of the present application. Therefore, the protection scope of the application is subject to the protection scope of the claims.
Claims (10)
1. The method for calculating the continuous detonation delay time of the electronic detonator is characterized by comprising the following steps of:
obtaining modeling parameters and a first function, wherein the modeling parameters comprise tunnel section size parameters, tunnel burial depth parameters, blast hole position parameters, blast hole loading quantity parameters and blast hole depth parameters, the first function is used for calculating the vibration speed transmitted to a target position from a blast hole under different blasting delay time, and the blasting delay time comprises the blasting delay time between a current blast hole and a next blast hole;
establishing at least one numerical model according to the modeling parameters, wherein the numerical model is used for simulating and calculating the vibration speed of single blasthole detonation transferred to a target;
calculating vibration speeds transmitted to a target position from the blasting position under different blasting delay time by using at least one numerical model to obtain vibration speed information;
fitting the first function according to the vibration speed information to obtain a vibration speed calculation formula;
and calculating the blasting delay time between the current blast hole and the next blast hole according to the vibration speed calculation formula.
2. The method for calculating the continuous detonation delay time of the electronic detonator according to claim 1, wherein the step of establishing at least one numerical model according to the modeling parameters comprises the steps of:
establishing a geometric model of single-hole continuous detonation of the tunnel electronic detonator by using the modeling information;
performing independent mapping division on the blast holes in the geometric model and setting boundary conditions on the numerical model;
selecting material parameters which accord with field reality;
and establishing a numerical model of the interaction between the rock, the explosive and the air during tunnel blasting by adopting a fluid-solid coupling method.
3. The method of claim 1, wherein obtaining the first function comprises:
acquiring first information, wherein the first information comprises physical parameters affecting blasting vibration propagation;
inquiring in a preset parameter library to obtain a dimension expression of each physical parameter;
and calculating according to the dimensional expression of each physical parameter to obtain a first function.
4. The method for calculating the continuous detonation delay time of the electronic detonator according to claim 3, wherein the calculating according to the dimensional expression of each physical parameter to obtain the first function comprises:
calculating the dimensionless expression of each physical parameter according to the dimensionless homogeneous principle to obtain second information, wherein the second information comprises a dimensionless expression corresponding to a blasting vibration period, a dimensionless expression corresponding to a vibration wave speed, a dimensionless expression corresponding to a blasting delay time, a dimensionless expression corresponding to a particle vibration speed and a dimensionless expression corresponding to a medium density;
calculating according to the second information to obtain at least two dimensionless expressions;
and calculating according to at least two dimensionless expressions to obtain the first function.
5. The continuous detonation delay time calculation system of the electronic detonator is characterized by comprising the following components:
the system comprises an acquisition module, a control module and a control module, wherein the acquisition module is used for acquiring modeling parameters and a first function, the modeling parameters comprise tunnel section size parameters, tunnel burial depth parameters, blast hole position parameters, blast hole charge quantity parameters and blast hole depth parameters, the first function is used for calculating vibration speeds transmitted to a target from a blast hole under different blasting delay time, and the blasting delay time comprises blasting delay time between a current blast hole and a next blast hole;
the establishing module is used for establishing at least one numerical model according to the modeling parameters, and the numerical model is used for simulating and calculating the vibration speed of the single blasthole detonation transferred to the target;
the simulation module is used for calculating the vibration speed transmitted to the target position from the blasting position under different blasting delay time by using at least one numerical model to obtain vibration speed information;
the first processing module is used for fitting the first function according to the vibration speed information to obtain a vibration speed calculation formula;
and the second processing module is used for calculating the blasting delay time between the current blast hole and the next blast hole according to the vibration speed calculation formula.
6. The electronic detonator sequential initiation delay time computing system of claim 5 wherein said setup module comprises:
the establishing unit is used for establishing a geometric model of single-hole continuous detonation of the tunnel electronic detonator by using the modeling information;
the dividing unit is used for carrying out independent mapping division on the blast holes in the geometric model and setting boundary conditions on the numerical model;
the first processing unit is used for selecting material parameters which accord with the field reality;
and the second processing unit is used for establishing a numerical model of the interaction between the rock, the explosive and the air during tunnel blasting by adopting a fluid-solid coupling method.
7. The electronic detonator sequential initiation delay time computing system of claim 5 wherein said acquisition module comprises:
a first acquisition unit configured to acquire first information including physical parameters affecting propagation of blasting vibration;
the query unit is used for querying in a preset parameter library to obtain a dimension expression of each physical parameter;
and the first calculation unit is used for calculating according to the dimensional expression of each physical parameter to obtain a first function.
8. The electronic detonator sequential initiation delay time computing system of claim 7 wherein said first computing unit comprises:
the second calculation unit is used for calculating the dimensionality expression of each physical parameter according to the dimensionality homogeneous principle to obtain second information, wherein the second information comprises a dimensionless expression corresponding to a blasting vibration period, a dimensionless expression corresponding to a vibration wave speed, a dimensionless expression corresponding to a blasting delay time, a dimensionless expression corresponding to a particle vibration speed and a dimensionless expression corresponding to a medium density;
the third calculation unit is used for calculating according to the second information to obtain at least two dimensionless expressions;
and the fourth calculation unit is used for calculating according to at least two dimensionless expressions to obtain the first function.
9. An electronic detonator continuous initiation delay time computing device, comprising:
a memory for storing a computer program;
a processor for implementing the steps of the method for calculating the continuous detonation delay time of the electronic detonator according to any one of claims 1 to 4 when executing the computer program.
10. A readable storage medium, characterized in that it has stored thereon a computer program which, when executed by a processor, implements the steps of the method for calculating the sequential detonation delay time of an electronic detonator according to any one of claims 1 to 4.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311145957.0A CN116882213B (en) | 2023-09-07 | 2023-09-07 | Method and system for calculating continuous detonation delay time of electronic detonator |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311145957.0A CN116882213B (en) | 2023-09-07 | 2023-09-07 | Method and system for calculating continuous detonation delay time of electronic detonator |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116882213A true CN116882213A (en) | 2023-10-13 |
CN116882213B CN116882213B (en) | 2023-11-24 |
Family
ID=88262578
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311145957.0A Active CN116882213B (en) | 2023-09-07 | 2023-09-07 | Method and system for calculating continuous detonation delay time of electronic detonator |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116882213B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20080071271A (en) * | 2007-01-30 | 2008-08-04 | 조선대학교산학협력단 | The method of prediction of blasting vibration by superposition on modeling data of single hole waveform |
WO2014036978A2 (en) * | 2012-09-10 | 2014-03-13 | 北京北方邦杰科技发展有限公司 | Detonator initiation control method, device, and initiator |
WO2014134913A1 (en) * | 2013-03-08 | 2014-09-12 | 葛洲坝易普力股份有限公司 | Detonation system having digital electronic detonator able to identify blast hole location and control method thereof |
CA2962230A1 (en) * | 2014-09-23 | 2016-03-31 | In Mu Ryu | Method of the shortest inter-hole delay blast and the blasting and delaying means |
CN109752085A (en) * | 2018-11-09 | 2019-05-14 | 北方爆破科技有限公司 | A kind of digital electric detonator hole-by-hole initiation vibration prediction method |
CN113251880A (en) * | 2021-05-25 | 2021-08-13 | 鞍钢矿业爆破有限公司 | Inter-hole time-delay blasting control vibration reduction method based on seismic wave energy attenuation |
CN115114833A (en) * | 2022-07-26 | 2022-09-27 | 武九铁路客运专线湖北有限责任公司 | Reasonable differential time control method for tunnel tunneling blasting |
-
2023
- 2023-09-07 CN CN202311145957.0A patent/CN116882213B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20080071271A (en) * | 2007-01-30 | 2008-08-04 | 조선대학교산학협력단 | The method of prediction of blasting vibration by superposition on modeling data of single hole waveform |
WO2014036978A2 (en) * | 2012-09-10 | 2014-03-13 | 北京北方邦杰科技发展有限公司 | Detonator initiation control method, device, and initiator |
WO2014134913A1 (en) * | 2013-03-08 | 2014-09-12 | 葛洲坝易普力股份有限公司 | Detonation system having digital electronic detonator able to identify blast hole location and control method thereof |
CA2962230A1 (en) * | 2014-09-23 | 2016-03-31 | In Mu Ryu | Method of the shortest inter-hole delay blast and the blasting and delaying means |
CN109752085A (en) * | 2018-11-09 | 2019-05-14 | 北方爆破科技有限公司 | A kind of digital electric detonator hole-by-hole initiation vibration prediction method |
CN113251880A (en) * | 2021-05-25 | 2021-08-13 | 鞍钢矿业爆破有限公司 | Inter-hole time-delay blasting control vibration reduction method based on seismic wave energy attenuation |
CN115114833A (en) * | 2022-07-26 | 2022-09-27 | 武九铁路客运专线湖北有限责任公司 | Reasonable differential time control method for tunnel tunneling blasting |
Non-Patent Citations (2)
Title |
---|
康永全;孟海利;薛里;孙崔源;管仁生;: "浅埋隧道爆破施工电子雷管降振延时时间的分析", 铁道建筑, no. 01 * |
张小军;汪旭光;于亚伦;杨德强;: "确定延时爆破间隔时间的新方法", 金属矿山, no. 03 * |
Also Published As
Publication number | Publication date |
---|---|
CN116882213B (en) | 2023-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109657310B (en) | Real-time dynamic virtual reality simulation method for gas explosion in coal mine | |
CN115325893B (en) | Tunnel blasting expansion excavation construction method and system | |
CN114781040B (en) | Building blasting vibration safety evaluation method, device, equipment and medium | |
CN109242144A (en) | A kind of surface blasting prediction technique, apparatus and system | |
JP6998014B2 (en) | Blasting method | |
CN116882213B (en) | Method and system for calculating continuous detonation delay time of electronic detonator | |
CN115600368A (en) | Deep hole blasting optimization method, device and equipment and readable storage medium | |
CN109752085A (en) | A kind of digital electric detonator hole-by-hole initiation vibration prediction method | |
CN113343503A (en) | Blasting method and device for predicting blasting vibration damage | |
Li et al. | Blasting vibration safety criterion analysis with equivalent elastic boundary: based on accurate loading model | |
CN113189644A (en) | Microseismic source positioning method and system | |
KR102636445B1 (en) | Blasting design system reflecting the situation of the blasting site and operation mathod of the same | |
CN114936342B (en) | Method, device, equipment and medium for calculating close-range blasting explosive quantity of protective object | |
Kara et al. | The latest generation of the electronic system enhanced safety and productivity | |
JP2023548563A (en) | Fracture density modeling system, method, and apparatus | |
CN114723863A (en) | Open bench blasting vibration speed cloud chart prediction method and system and electronic equipment | |
CN114564775A (en) | Double-scale coupling numerical simulation method suitable for building tunnel above goaf | |
CN114578428A (en) | Efficient construction method and device for velocity field | |
CN111259095A (en) | Method, device and equipment for calculating boundary of ore rock | |
Milanović et al. | Application of softwares for drilling and blasting | |
CN117610315B (en) | Tunnel intelligent blasting design system based on multiple geological information | |
Toraño et al. | Analysis of the environmental effects of ground vibrations produced by blasting in quarries | |
Ercins et al. | Performance Analysis of the Explosion Applications Realized with Electronic Ignition System at Different Times in the Same Field | |
CN118094835A (en) | Method, device, equipment, medium and product for predicting spatial distribution of detonation stack grade | |
Stanković et al. | A practical approach to the ground oscillation velocity measurement method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |