CN116258285B - Porous small-clear-distance tunnel blasting vibration speed prediction method, device, equipment and medium - Google Patents
Porous small-clear-distance tunnel blasting vibration speed prediction method, device, equipment and medium Download PDFInfo
- Publication number
- CN116258285B CN116258285B CN202310545335.0A CN202310545335A CN116258285B CN 116258285 B CN116258285 B CN 116258285B CN 202310545335 A CN202310545335 A CN 202310545335A CN 116258285 B CN116258285 B CN 116258285B
- Authority
- CN
- China
- Prior art keywords
- tunnel
- vibration speed
- vibration
- blasting
- speed data
- 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.)
- Active
Links
- 238000005422 blasting Methods 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims abstract description 98
- 230000009467 reduction Effects 0.000 claims abstract description 91
- 238000012544 monitoring process Methods 0.000 claims abstract description 56
- 230000005641 tunneling Effects 0.000 claims abstract description 56
- 238000004364 calculation method Methods 0.000 claims abstract description 40
- 230000008569 process Effects 0.000 claims description 24
- 238000001514 detection method Methods 0.000 claims description 19
- 239000011435 rock Substances 0.000 claims description 15
- 238000004880 explosion Methods 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 7
- 230000009172 bursting Effects 0.000 claims description 4
- 239000002360 explosive Substances 0.000 claims description 3
- 238000010276 construction Methods 0.000 abstract description 6
- 238000012360 testing method Methods 0.000 description 16
- 238000004891 communication Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 8
- 230000005236 sound signal Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012545 processing Methods 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
- 230000002411 adverse Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000009412 basement excavation 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
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/08—Construction
-
- 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)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Entrepreneurship & Innovation (AREA)
- Algebra (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Primary Health Care (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The application provides a method, a device, equipment and a medium for predicting the blasting vibration speed of a porous small-clear-distance tunnel, which relate to the technical field of tunnel blasting and comprise a first tunnel, a second tunnel and a third tunnel, wherein the method comprises the steps of acquiring first vibration speed data of the first tunnel during the second tunnel tunneling, and calculating according to the first vibration speed data to obtain a blasting vibration speed prediction formula and a first reduction coefficient; according to the first vibration speed data and the second vibration speed data, third vibration speed data of the first tunnel and fourth vibration speed data of the second tunnel during third tunneling are obtained through calculation, and the method is used for solving the problems that in the prior art, the acquisition of the blasting vibration speed of the porous small-clearance tunnel is mainly based on-site monitoring, a related theoretical calculation method is lacked, and the blasting construction scheme cannot be optimized in advance.
Description
Technical Field
The application relates to the technical field of tunnel blasting, in particular to a method, a device, equipment and a medium for predicting the blasting vibration speed of a porous small-clear-distance tunnel.
Background
With the development of tunnels and underground engineering, the idea of multi-hole (three-hole and more) small clear-distance tunnels is increasingly applied to practical engineering. When a tunnel is excavated by using a drilling and blasting method, blasting vibration generated by construction of a backward tunnel often has adverse effects on primary supports and secondary liners of a preceding tunnel, and the primary supports and the secondary liners ensure safety mainly by controlling the maximum vibration speed not to exceed a safety allowable vibration speed standard. At present, the acquisition of the blasting vibration speed of the porous small-clear-distance tunnel is mainly based on-site monitoring, and a related theoretical calculation method is lacked, so that the blasting construction scheme is difficult to optimize in advance.
Disclosure of Invention
The application aims to provide a method, a device, equipment and a medium for predicting the blasting vibration speed of a porous small-clear-distance tunnel so as to solve the problems. In order to achieve the above purpose, the technical scheme adopted by the application is as follows:
the application provides a method for predicting blasting vibration speed of a porous small-clear-distance tunnel, which comprises a first tunnel which is excavated, a second tunnel and a third tunnel which are positioned at two sides of the first tunnel, wherein the first tunnel and the second tunnel are tunnels to be excavated;
acquiring first vibration speed data of a first tunnel during tunneling of a second tunnel, and calculating according to the first vibration speed data to obtain a blasting vibration speed prediction formula and a first reduction coefficient, wherein the first reduction coefficient is a reduction coefficient when blasting vibration waves pass through the first tunnel;
calculating to obtain second vibration speed data of a second tunnel according to the blasting vibration speed prediction formula and the blasting center distance during tunneling of the second tunnel;
calculating a second reduction coefficient according to the first vibration speed data and the second vibration speed data, wherein the second reduction coefficient is a reduction coefficient when blasting vibration waves pass through surrounding rock between the first tunnel and the second tunnel;
and calculating third vibration speed data of the first tunnel and fourth vibration speed data of the second tunnel during third tunneling according to the first and second reduction coefficients and the explosion vibration speed prediction formula.
In a second aspect, the application also provides a method for predicting the blasting vibration speed of the porous small-clear-distance tunnel, which comprises a first tunnel which is excavated, a second tunnel and a third tunnel which are positioned at two sides of the first tunnel, wherein the first tunnel and the second tunnel are tunnels to be excavated;
the acquisition module is used for: the method comprises the steps of obtaining first vibration speed data of a first tunnel during tunneling of a second tunnel, and calculating according to the first vibration speed data to obtain a blasting vibration speed prediction formula and a first reduction coefficient, wherein the first reduction coefficient is a reduction coefficient when blasting vibration waves pass through the first tunnel;
a first calculation module: the second vibration speed data of the second tunnel is obtained through calculation according to the explosion vibration speed prediction formula and the explosion center distance during the second tunnel tunneling;
a second calculation module: the method comprises the steps of calculating a second reduction coefficient according to first vibration speed data and second vibration speed data, wherein the second reduction coefficient is a reduction coefficient when blasting vibration waves pass through surrounding rock between a first tunnel and a second tunnel;
a third calculation module: and the third vibration speed data of the first tunnel and the fourth vibration speed data of the second tunnel are obtained through calculation according to the first and second reduction coefficients and the explosion vibration speed prediction formula.
In a third aspect, the present application further provides a device for predicting a blasting vibration velocity of a porous small-clear-distance tunnel, including:
a memory for storing a computer program;
and the processor is used for realizing the step of the porous small-clear-distance tunnel blasting vibration speed prediction method when executing the computer program.
In a fourth aspect, the present application also provides a computer readable storage medium, on which a computer program is stored, which when executed by a processor, implements the steps of the method for predicting a blasting vibration velocity based on a porous small clear distance tunnel.
The beneficial effects of the application are as follows:
aiming at the problem that the blasting vibration speed of the porous small-clear-distance tunnel is difficult to predict, the application provides a small-clear-distance tunnel blasting vibration speed calculation method by researching the reduction coefficient of blasting seismic waves in different media. By using the vibration speed calculation method, when the back tunnel is in blasting construction, the blasting vibration speed of the front tunnel can be predicted according to the method so as to judge whether the lining structure of the front tunnel is damaged or not, and then the blasting scheme is adjusted and optimized, so that accurate guidance is provided for construction, and the safety of the construction is ensured.
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 predicting the blasting vibration velocity of a porous small-clear-distance tunnel according to an embodiment of the application;
FIG. 2 is a schematic diagram of the positional relationship of a porous small clear-distance tunnel according to an embodiment of the present application;
FIG. 3 is a schematic diagram of the location of a monitoring point in an embodiment of the present application;
FIG. 4 is a schematic diagram illustrating a position of a second predicted point according to an embodiment of the present application;
FIG. 5 is a schematic diagram of the locations of a first tunnel and a second tunnel;
FIG. 6 is a second schematic diagram of the position of the point E to be solved according to the embodiment of the present application;
FIG. 7 is a schematic structural diagram of a device for predicting the blasting vibration velocity of a porous small-clear-distance tunnel according to an embodiment of the application;
fig. 8 is a schematic structural diagram of a porous small-clear-distance tunnel blasting vibration velocity prediction device according to an embodiment of the application.
The marks in the figure:
800. porous small clear distance tunnel blasting vibration speed prediction equipment; 801. a processor; 802. a memory; 803. a multimedia component; 804. an I/O interface; 805. a communication component.
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 porous small-clear-distance tunnel blasting vibration speed prediction method.
Referring to fig. 1 and 2, fig. 2 shows that the method includes an excavated first tunnel, a second tunnel and a third tunnel to be excavated, wherein the second tunnel is positioned at the right side of the first tunnel, the third tunnel is positioned at the left side of the first tunnel, and the diameters of the first tunnel, the second tunnel and the third tunnel are the same;
s1, acquiring first vibration speed data of a first tunnel during tunneling of a second tunnel, and calculating a blasting vibration speed prediction formula and a first reduction coefficient according to the first vibration speed data, wherein the first reduction coefficient is a reduction coefficient when blasting vibration waves pass through the first tunnel;
specifically, the step S1 includes:
s11, symmetrically arranging a plurality of monitoring points on two side walls of a first tunnel along the tunneling direction, and forming a pair of monitoring points by a group of symmetrical monitoring points;
specifically, as shown in fig. 3, k monitoring points are arranged at the excavated part of the right side wall of the first tunnel, which are respectively，，…The method comprises the steps of carrying out a first treatment on the surface of the K monitoring points are arranged on the excavated part of the left side wall of the first tunnel, which are respectively，，…The method comprises the steps of carrying out a first treatment on the surface of the It should be noted that the monitoring points on the left and right sides should be arranged at the same mileage and the same elevation. Preferably, a three-way blasting vibration meter NUBOX-6016 is adopted for monitoring, and the vibration meter is tightly attached to the tunnel excavation rock wall.
In the field of tunnel blasting, monitoring points can only be arranged in an excavated tunnel, and only the excavated tunnel needs to monitor or predict the blasting vibration speed so as to judge whether damage can be caused to the excavated tunnel.
S12, when the monitoring points are all arranged on the primary support of the side wall, obtaining the vibration speed of all the monitoring points in the second tunneling process、、…；、、…；
S13, fitting a Sargassy formula according to the vibration velocity of the right monitoring point to obtain a blasting vibration velocity prediction formula;
specifically, the second tunnel is positioned on the right side of the first tunnel, so that the vibration speed of the monitoring point on the right side of the first tunnel、、…Fitting the Sarcopsis formula to obtain a blasting vibration velocity prediction formula of the stratum where the tunnel is located:
（1）。
in the method, in the process of the application,the predicted vibration velocity is indicated to be the same,representing coefficients related to the conditions of the blastfield,represented as the seismic wave attenuation coefficient,indicating the explosive quantity.
S14, calculating a first vibration speed ratio of each pair of monitoring points according to the vibration speed of each pair of monitoring points;
s15, calculating average values of all the first vibration ratio to obtain a first reduction coefficient;
（2）。
in the method, in the process of the application,representing the first reduction coefficient.
Based on the above embodiment, the method further includes:
s2, calculating second vibration speed data of the second tunnel according to the explosion vibration speed prediction formula and the explosion center distance during tunneling of the second tunnel, wherein the explosion vibration speed prediction formula is also applicable to the third tunnel and the second tunnel because the third tunnel and the second tunnel are the same as the stratum where the first tunnel is located;
specifically, the step S2 includes:
s21, selecting a plurality of first prediction points on the side wall, close to the first tunnel, of the second tunnel along the tunneling direction;
specifically, as shown in fig. 3, m first predicted points are selected on the side wall of the portion to be excavated of the second tunnel, where the first predicted points are respectively:，，…the positions of the first predicted point and the monitoring point are in one-to-one correspondence
S22, acquiring the core explosion distances of all the first predicted points during the second tunneling;
s23, obtaining the vibration speed of each first predicted point in the second tunneling process according to a blasting vibration speed prediction formula and a blasting center distance calculation, and forming second vibration speed data by the vibration speeds of all the first predicted points、……。
Based on the above embodiment, the method further includes:
s3, calculating a second reduction coefficient according to the first vibration speed data and the second vibration speed data, wherein the second reduction coefficient is a reduction coefficient when blasting vibration waves pass through surrounding rocks between the first tunnel and the second tunnel;
specifically, the step S3 includes:
s31, forming a plurality of pairs of detection points by the detection points close to the second tunnel and the first detection points close to the first tunnel;
s32, calculating a second vibration velocity ratio of each pair of detection points according to the vibration velocity of each pair of detection points;
s33, calculating the average value of all second vibration ratio values to obtain a second reduction coefficient:
（3）。
in the method, in the process of the application,representing a second reduction factor.
Based on the above embodiment, the method further includes:
s4, calculating third vibration speed data of the first tunnel and fourth vibration speed data of the second tunnel during third tunneling according to the first and second reduction coefficients and the blasting vibration speed prediction formula:
specifically, the third vibration speed data of the first tunnel includes a vibration speed of the left side of the first tunnelAnd right vibration velocity：
（4）。
In the method, in the process of the application,the explosion center distance of any point on the left side wall of the first tunnel is shown when the third tunnel is exploded;
（5）。
specifically, the fourth vibration speed data of the second tunnel includes a vibration speed of the left side of the second tunnelAnd right vibration velocity：
（6）。
（7）。
Because the diameters of the three tunnels are the sameThe method is suitable for three tunnels.
The third vibration speed data of the first tunnel and the fourth vibration speed data of the second tunnel may be vibration speeds of any positions on the side wall.
Example 2:
the embodiment provides a porous small-clear-distance tunnel blasting vibration speed prediction method.
Referring to fig. 4, the method includes a first tunnel which is excavated, a second tunnel to be excavated, a third tunnel and a fourth tunnel, wherein the second tunnel is positioned on the right side of the first tunnel, the third tunnel is positioned on the left side of the first tunnel, the fourth tunnel is positioned on the right side of the second tunnel, the diameters of the first tunnel, the second tunnel, the third tunnel and the fourth tunnel are the same, and the method includes the following steps:
s1, acquiring first vibration speed data of a first tunnel during tunneling of a second tunnel, and calculating a blasting vibration speed prediction formula and a first reduction coefficient according to the first vibration speed data, wherein the first reduction coefficient is a reduction coefficient when blasting vibration waves pass through the first tunnel;
s2, calculating to obtain second vibration speed data of a second tunnel according to the blasting vibration speed prediction formula and the blasting center distance during tunneling of the second tunnel;
s3, calculating a second reduction coefficient according to the first vibration speed data and the second vibration speed data, wherein the second reduction coefficient is a reduction coefficient when blasting vibration waves pass through surrounding rocks between the first tunnel and the second tunnel;
s4, calculating third vibration speed data of the first tunnel and fourth vibration speed data of the second tunnel during tunneling of the third tunnel according to the first reduction coefficient and the second reduction coefficient.
S5, selecting n second prediction points on the side wall, close to the first tunnel, of the third tunnel along the tunneling direction, wherein the n second prediction points are respectively:，，…；
the second prediction points are positioned on the side wall of the third tunnel which is not excavated, and the positions of the second prediction points and the positions of the monitoring points are in one-to-one correspondence;
s6, acquiring the explosion center distances of all the second predicted points in the third tunneling process, calculating the vibration speed of each second predicted point according to an explosion vibration speed prediction formula and the explosion center distance, and calculating the vibration speed of each second predicted point according to the explosion vibration speed prediction formula、……；
S7, forming a plurality of pairs of predicted points by the monitored points close to the third tunnel and the second predicted points;
s8, calculating a third vibration speed ratio of each pair of predicted points according to the vibration speed of each pair of predicted points;
s9, calculating average values of all third vibration ratios to obtain a third reduction coefficient, wherein the third reduction coefficient is a reduction coefficient when blasting vibration waves pass through surrounding rocks between the first tunnel and the third tunnel;
（8）。
in the method, in the process of the application,representing the third reduction factor.
S10, calculating fifth vibration speed data of a third tunnel during tunneling of the fourth tunnel according to the first reduction coefficient, the second reduction coefficient and the third reduction coefficient;
specifically, the fifth vibration speed data of the third tunnel includes a left vibration speedAnd right side vibration velocity：
（9）。
In the method, in the process of the application,the bursting distance of any point on the right side wall of the second tunnel is represented when the fourth tunnel bursts;
（10）。
similarly, the vibration speed data of the first tunnel and the second tunnel can be based onAndcalculating to obtain;
referring to fig. 6, based on the above embodiment, the method further includes, when the line type of the tunnels is irregular, the thickness of the middle rocks formed between the tunnels is a variable, and at this time, calculating a coefficient λ related to the formation property according to a series of vibration ratios and the corresponding thickness x of the middle rocks on site:
（11）。
after the coefficient lambda related to the stratum property is determined, the reduction coefficient corresponding to the thickness of any surrounding rock between tunnels can be obtained.
In FIG. 6, the thickness of the intermediate rock formed between the second tunnel and the first tunnel is varied byHas been calculated in example 1 and can therefore be based onCalculating a coefficient related to the formation property between the second tunnel and the first tunnel。
Example 3: the embodiment provides a porous small-clear-distance tunnel blasting vibration speed prediction method.
The method comprises an excavated first tunnel, a second tunnel and a third tunnel which are positioned at two sides of the first tunnel and are to be excavated, wherein the second tunnel is positioned at the right side of the first tunnel, the third tunnel is positioned at the left side of the first tunnel, and the diameters of the first tunnel, the second tunnel and the third tunnel are the same;
s1, acquiring first vibration speed data of a first tunnel during tunneling of a second tunnel, and calculating a blasting vibration speed prediction formula and a first reduction coefficient according to the first vibration speed data, wherein the first reduction coefficient is a reduction coefficient when blasting vibration waves pass through the first tunnel;
s11, symmetrically arranging a plurality of monitoring points on two side walls of a first tunnel along the tunneling direction, and forming a pair of monitoring points by a group of symmetrical monitoring points;
s12, when the monitoring points are arranged on the two liners of the first tunnel, the monitoring points are the first testing pointsObtaining vibration speed of a first test point of a second tunnel in one tunneling process；
S13, acquiring a first length of primary tunneling of a second tunnelArranging a second test point at a position with a first length from the first test pointThe second test point is arranged on a primary support of the first tunnel;
s14, acquiring vibration speed of a second test point during secondary tunneling of a second tunnel；
S15, obtaining an attenuation factor according to the ratio of the vibration speed of the first test point to the vibration speed of the second test point:
（12）。
in the method, in the process of the application,representing the attenuation factor of the second liner.
S2, calculating to obtain second vibration speed data of a second tunnel according to the blasting vibration speed prediction formula and the blasting center distance during tunneling of the second tunnel;
s3, calculating a second reduction coefficient according to the first vibration speed data and the second vibration speed data, wherein the second reduction coefficient is a reduction coefficient when blasting vibration waves pass through surrounding rocks between the first tunnel and the second tunnel;
s4, as shown in FIG. 6, calculating to obtain third vibration speed data of the first tunnel and fourth vibration speed data of the second tunnel during the third tunnel blasting according to the first reduction coefficient, the second reduction coefficient and the attenuation factor.
The third vibration speed data comprises the vibration speed of the left side of the first tunnelAnd right vibration velocity：
（13）。
（14）。
Specifically, the fourth vibration velocity data includes a vibration velocity of the left side of the second tunnelAnd right vibration velocity：
（15）。
（16）。
Representing the thickness of the strata in the first tunnel and the second tunnel.
Example 4:
as shown in fig. 7, the embodiment provides a device for predicting the blasting vibration speed of a porous small-clear-distance tunnel, which comprises a first tunnel which is excavated, a second tunnel and a third tunnel which are positioned at two sides of the first tunnel, wherein the first tunnel and the second tunnel are tunnels to be excavated;
the acquisition module is used for: the method comprises the steps of obtaining first vibration speed data of a first tunnel during tunneling of a second tunnel, and calculating according to the first vibration speed data to obtain a blasting vibration speed prediction formula and a first reduction coefficient, wherein the first reduction coefficient is a reduction coefficient when blasting vibration waves pass through the first tunnel;
a first calculation module: the second vibration speed data of the second tunnel is obtained through calculation according to the explosion vibration speed prediction formula and the explosion center distance during the second tunnel tunneling;
a second calculation module: the method comprises the steps of calculating a second reduction coefficient according to first vibration speed data and second vibration speed data, wherein the second reduction coefficient is a reduction coefficient when blasting vibration waves pass through surrounding rock between a first tunnel and a second tunnel;
a third calculation module: and the third vibration speed data of the first tunnel and the fourth vibration speed data of the second tunnel are obtained through calculation according to the first and second reduction coefficients and the explosion vibration speed prediction formula.
Based on the above embodiments, the acquisition module includes:
an arrangement unit: a plurality of monitoring points are symmetrically arranged on the two side walls of the first tunnel along the tunneling direction, a pair of monitoring points is formed by a group of symmetrical monitoring points;
a first acquisition unit: when all the monitoring points are arranged on the primary support of the side wall, acquiring the vibration speed of all the monitoring points in the second tunneling process;
fitting unit: the method comprises the steps of fitting a Sargassy formula according to the vibration speed of a right monitoring point to obtain a blasting vibration speed prediction formula;
a first calculation unit: the first vibration speed ratio of each pair of monitoring points is obtained through calculation according to the vibration speed of each pair of monitoring points;
a second calculation unit: and calculating the average value of all the first vibration ratio to obtain a first reduction coefficient.
Based on the above embodiments, the first calculation module includes:
the selecting unit: the method comprises the steps that a plurality of first prediction points are selected on the side wall, close to a first tunnel, of a second tunnel along the tunneling direction, wherein the positions of the first prediction points and the positions of the monitoring points are in one-to-one correspondence;
a second acquisition unit: the method comprises the steps of acquiring the core bursting distance of all first predicted points during the second tunneling;
a third calculation unit: and the vibration speed of each first predicted point in the second tunneling process is calculated according to the blasting vibration speed prediction formula, and the vibration speeds of all the first predicted points form second vibration speed data.
Based on the above embodiment, the second calculation module includes:
fourth acquisition unit: the system comprises a plurality of pairs of detection points, wherein the pairs of detection points consist of a monitoring point close to a second tunnel and a first prediction point close to a first tunnel;
a fourth calculation unit: the second vibration speed ratio of each pair of detection points is obtained through calculation according to the vibration speed of each pair of detection points;
a fifth calculation unit: and calculating the average value of all the second vibration ratio to obtain a second reduction coefficient.
Based on the above embodiment, the method further comprises a fourth tunnel to be excavated, wherein the fourth tunnel is located at the other side of the second tunnel:
a second selection unit: the method comprises the steps that a plurality of second prediction points are selected on the side wall, close to a first tunnel, of a third tunnel along the tunneling direction, and the positions of the second prediction points and the positions of the monitoring points are in one-to-one correspondence;
fifth acquisition unit: the method comprises the steps of obtaining the explosion center distance of all second predicted points in the third tunneling process, and obtaining the vibration speed of each second predicted point according to an explosion vibration speed prediction formula and the explosion center distance;
sixth acquisition unit: the system comprises a plurality of pairs of prediction points, wherein the pairs of prediction points are formed by a monitoring point close to a third tunnel and a second prediction point;
a fifth calculation unit: the third vibration speed ratio of each pair of predicted points is obtained by calculation according to the vibration speed of each pair of predicted points;
a sixth calculation unit: the method comprises the steps of calculating the average value of all third vibration ratios to obtain a third reduction coefficient, wherein the third reduction coefficient is a reduction coefficient when blasting vibration waves pass through surrounding rock between a first tunnel and a third tunnel;
seventh calculation unit: and the fifth vibration speed data of the third tunnel during the fourth tunneling is calculated according to the first reduction coefficient, the second reduction coefficient and the third reduction coefficient.
Based on the above embodiment, when the monitoring point is arranged on the two liners of the first tunnel, the monitoring point is made to be the first test point, and further includes:
seventh acquisition unit: the method comprises the steps of obtaining the vibration speed of a first test point of a second tunnel in one tunneling process;
eighth acquisition unit: the method comprises the steps of acquiring a first length of primary tunneling of a second tunnel, and arranging a second test point at a position which is away from the first test point by the first length, wherein the second test point is arranged on a primary support of the first tunnel;
a ninth acquisition unit: the vibration speed of the second test point is obtained when the second tunnel is driven for the second time;
an eighth calculation unit: the damping factor is obtained from the ratio of the vibration speed of the first test point to the vibration speed of the second test point;
a ninth calculation unit: and the third vibration speed data of the first tunnel and the fourth vibration speed data of the second tunnel are obtained through calculation according to the first reduction coefficient, the second reduction coefficient and the attenuation factor.
It should be noted that, regarding the apparatus in the above embodiments, the specific manner in which the respective modules perform the operations has been described in detail in the embodiments regarding the method, and will not be described in detail herein.
Example 5:
corresponding to the above method embodiment, a porous small clear distance tunnel blasting vibration speed prediction device is further provided in this embodiment, and a porous small clear distance tunnel blasting vibration speed prediction device described below and a porous small clear distance tunnel blasting vibration speed prediction method described above can be referred to correspondingly.
Fig. 8 is a block diagram illustrating a porous small clear distance tunnel blasting vibration velocity prediction apparatus 800 according to an exemplary embodiment. As shown in fig. 8, the porous small clear distance tunnel blasting vibration velocity prediction apparatus 800 may include: a processor 801, a memory 802. The porous small clear distance tunnel blasting vibration velocity prediction apparatus 800 may further comprise 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 apparatus 800 for predicting a vibration velocity of a blasting in a small clear distance tunnel, so as to complete all or part of the steps in the method for predicting a vibration velocity of a blasting in a small clear distance tunnel. The memory 802 is used to store various types of data to support the operation of the multi-hole small clear distance tunnel blasting vibration prediction device 800, which may include, for example, instructions for any application or method operating on the multi-hole small clear distance tunnel blasting vibration prediction device 800, as well as application-related data, such as contact data, messaging, 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 porous small clear distance tunnel blasting vibration velocity prediction 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 multi-hole small clear tunnel blasting vibration velocity prediction apparatus 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 apparatus (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 multi-hole small clear tunnel blasting vibration velocity prediction method described above.
In another exemplary embodiment, a computer readable storage medium is also provided, comprising program instructions which, when executed by a processor, implement the steps of the porous small clear tunnel blasting vibration velocity prediction 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 multi-hole small clear distance tunnel blasting vibration velocity prediction apparatus 800 to perform the multi-hole small clear distance tunnel blasting vibration velocity prediction method described above.
Example 6:
corresponding to the above method embodiment, a computer readable storage medium is also provided in this embodiment, and a computer readable storage medium described below and a method for predicting a blasting vibration velocity of a porous small-clear tunnel described above may be referred to correspondingly with each other.
A computer readable storage medium, on which a computer program is stored, which when executed by a processor implements the steps of the method for predicting the vibration velocity of a multi-hole small clear-distance tunnel blasting of the above method embodiment.
The computer 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, which may store program codes.
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.
Claims (6)
1. The method is characterized by comprising a first tunnel which is excavated, a second tunnel and a third tunnel which are positioned at two sides of the first tunnel, wherein the second tunnel and the third tunnel are tunnels to be excavated;
acquiring first vibration speed data of a first tunnel during second tunneling, and calculating according to the first vibration speed data to obtain a blasting vibration speed prediction formula and a first reduction coefficient, wherein the first reduction coefficient is a reduction coefficient when blasting vibration waves pass through the first tunnel, and comprises the following steps:
symmetrically arranging a plurality of monitoring points on two side walls of the first tunnel along the tunneling direction, wherein a pair of monitoring points are formed by a group of symmetrical monitoring points;
when all the monitoring points are arranged on the primary support of the side wall, the vibration speed of all the monitoring points in the second tunneling process is obtained;
fitting a Sarkowski formula according to the vibration speed of the right monitoring point to obtain a blasting vibration speed prediction formula;
calculating a first vibration speed ratio of each pair of monitoring points according to the vibration speed of each pair of monitoring points;
calculating the average of all the first vibration ratios to obtain a first reduction coefficient;
calculating to obtain second vibration speed data of a second tunnel according to the blasting vibration speed prediction formula and the blasting center distance during tunneling of the second tunnel;
according to the first vibration speed data and the second vibration speed data, a second reduction coefficient is obtained through calculation, wherein the second reduction coefficient is a reduction coefficient when blasting vibration waves pass through surrounding rock between the first tunnel and the second tunnel, and the method comprises the following steps:
a plurality of pairs of detection points are formed by the detection points close to the second tunnel and the first detection points close to the first tunnel;
calculating a second vibration velocity ratio of each pair of detection points according to the vibration velocity of each pair of detection points;
calculating the average value of all the second vibration ratio to obtain a second reduction coefficient;
according to the first reduction coefficient, the second reduction coefficient and the blasting vibration speed prediction formula, third vibration speed data of the first tunnel and fourth vibration speed data of the second tunnel during third tunneling are obtained through calculation;
the third vibration speed data of the first tunnel comprises the vibration speed of the left side of the first tunnelAnd right vibration +.>：
；
In the method, in the process of the application,the explosion center distance of any point of the left side wall of the first tunnel when the third tunnel is exploded is represented by +.>Representing coefficients related to the conditions of the blasting field, < +.>Expressed as seismic attenuation coefficient>The blasting explosive quantity is represented;
;
the fourth vibration speed data of the second tunnel comprises the vibration speed of the left side of the second tunnelAnd right vibration +.>：
;
；
In the method, in the process of the application,representing the first reduction coefficient, ">Representing a second reduction factor.
2. The method for predicting the explosion vibration speed of the porous small clear distance tunnel according to claim 1, wherein the calculating the second vibration speed data of the second tunnel according to the explosion vibration speed prediction formula and the explosion center distance during the second tunnel tunneling comprises the following steps:
selecting a plurality of first predicted points on the side wall, close to the first tunnel, of the second tunnel along the tunneling direction, wherein the positions of the first predicted points and the positions of the monitoring points are in one-to-one correspondence;
the core bursting distance of all the first predicted points in the second tunneling process is obtained;
and calculating the vibration speed of each first predicted point in the second tunneling process according to the blasting vibration speed prediction formula, and forming second vibration speed data by the vibration speeds of all the first predicted points.
3. The device is characterized by comprising a first tunnel which is excavated, a second tunnel and a third tunnel which are positioned at two sides of the first tunnel, wherein the second tunnel and the third tunnel are tunnels to be excavated;
the acquisition module is used for: the method is used for acquiring first vibration speed data of a first tunnel during tunneling of a second tunnel, calculating a blasting vibration speed prediction formula and a first reduction coefficient according to the first vibration speed data, wherein the first reduction coefficient is a reduction coefficient when blasting vibration waves pass through the first tunnel, and comprises the following steps:
symmetrically arranging a plurality of monitoring points on two side walls of the first tunnel along the tunneling direction, wherein a pair of monitoring points are formed by a group of symmetrical monitoring points;
when all the monitoring points are arranged on the primary support of the side wall, the vibration speed of all the monitoring points in the second tunneling process is obtained;
fitting a Sarkowski formula according to the vibration speed of the right monitoring point to obtain a blasting vibration speed prediction formula;
calculating a first vibration speed ratio of each pair of monitoring points according to the vibration speed of each pair of monitoring points;
calculating the average of all the first vibration ratios to obtain a first reduction coefficient;
a first calculation module: the second vibration speed data of the second tunnel is obtained through calculation according to the explosion vibration speed prediction formula and the explosion center distance during the second tunnel tunneling;
a second calculation module: the method is used for calculating a second reduction coefficient according to the first vibration speed data and the second vibration speed data, wherein the second reduction coefficient is a reduction coefficient when blasting vibration waves pass through surrounding rock between the first tunnel and the second tunnel, and comprises the following steps:
a plurality of pairs of detection points are formed by the detection points close to the second tunnel and the first detection points close to the first tunnel;
calculating a second vibration velocity ratio of each pair of detection points according to the vibration velocity of each pair of detection points;
a third calculation module: the method comprises the steps of calculating third vibration speed data of a first tunnel and fourth vibration speed data of a second tunnel during third tunneling according to a first reduction coefficient, a second reduction coefficient and a blasting vibration speed prediction formula;
the third vibration speed data of the first tunnel comprises the vibration speed of the left side of the first tunnelAnd right vibration +.>：
；
In the method, in the process of the application,the explosion center distance of any point of the left side wall of the first tunnel when the third tunnel is exploded is represented by +.>Representing coefficients related to the conditions of the blasting field, < +.>Expressed as seismic attenuation coefficient>The blasting explosive quantity is represented;
;
the fourth vibration speed data of the second tunnel comprises the vibration speed of the left side of the second tunnelAnd right vibration +.>：
;
；
In the method, in the process of the application,representing the first reduction coefficient, ">Representing a second reduction factor.
4. The porous small clear distance tunnel blasting vibration velocity prediction device according to claim 3, wherein the first calculation module comprises:
the selecting unit: the method comprises the steps that a plurality of first prediction points are selected on the side wall, close to a first tunnel, of a second tunnel along the tunneling direction, wherein the positions of the first prediction points and the positions of the monitoring points are in one-to-one correspondence;
a second acquisition unit: the method comprises the steps of acquiring the core bursting distance of all first predicted points during the second tunneling;
a third calculation unit: and the vibration speed of each first predicted point in the second tunneling process is calculated according to the blasting vibration speed prediction formula, and the vibration speeds of all the first predicted points form second vibration speed data.
5. A porous small clear distance tunnel blasting vibration velocity prediction apparatus, comprising:
a memory for storing a computer program;
a processor for implementing the steps of the porous small clear distance tunnel blasting vibration velocity prediction method according to any one of claims 1 to 2 when executing the computer program.
6. A computer readable storage medium, characterized in that the computer readable storage medium has stored thereon a computer program which, when executed by a processor, implements the steps of the porous small clear tunnel blasting vibration velocity prediction method according to any one of claims 1 to 2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310545335.0A CN116258285B (en) | 2023-05-16 | 2023-05-16 | Porous small-clear-distance tunnel blasting vibration speed prediction method, device, equipment and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310545335.0A CN116258285B (en) | 2023-05-16 | 2023-05-16 | Porous small-clear-distance tunnel blasting vibration speed prediction method, device, equipment and medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116258285A CN116258285A (en) | 2023-06-13 |
| CN116258285B true CN116258285B (en) | 2023-08-25 |
Family
ID=86686523
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310545335.0A Active CN116258285B (en) | 2023-05-16 | 2023-05-16 | Porous small-clear-distance tunnel blasting vibration speed prediction method, device, equipment and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116258285B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003075247A (en) * | 2001-08-31 | 2003-03-12 | Nishimatsu Constr Co Ltd | Method for measuring elastic wave speed in face in tunnel pit |
| CN110736400A (en) * | 2019-10-22 | 2020-01-31 | 长江重庆航道工程局 | underwater drilling blasting vibration velocity calculation method considering rock internal structure |
| CN111426243A (en) * | 2020-01-20 | 2020-07-17 | 重庆中环建设有限公司 | Blasting parameter determination method based on differential vibration synthesis of different blast holes in near zone of blasting source |
| CN115371791A (en) * | 2022-07-19 | 2022-11-22 | 青岛理工大学 | Underground pipeline vibration velocity measuring method and system |
| CN115408929A (en) * | 2022-07-22 | 2022-11-29 | 北京交通大学 | Data processing method and device for tunnel driving construction prediction |
| CN115688428A (en) * | 2022-11-01 | 2023-02-03 | 青岛理工大学 | Method, device and equipment for calculating critical value of explosive quantity of blasting near zone supporting structure |
| CN115935753A (en) * | 2022-12-26 | 2023-04-07 | 西南交通大学 | Method for simulating and analyzing influence of small-clear-distance tunnel blasting on adjacent tunnel |
-
2023
- 2023-05-16 CN CN202310545335.0A patent/CN116258285B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003075247A (en) * | 2001-08-31 | 2003-03-12 | Nishimatsu Constr Co Ltd | Method for measuring elastic wave speed in face in tunnel pit |
| CN110736400A (en) * | 2019-10-22 | 2020-01-31 | 长江重庆航道工程局 | underwater drilling blasting vibration velocity calculation method considering rock internal structure |
| CN111426243A (en) * | 2020-01-20 | 2020-07-17 | 重庆中环建设有限公司 | Blasting parameter determination method based on differential vibration synthesis of different blast holes in near zone of blasting source |
| CN115371791A (en) * | 2022-07-19 | 2022-11-22 | 青岛理工大学 | Underground pipeline vibration velocity measuring method and system |
| CN115408929A (en) * | 2022-07-22 | 2022-11-29 | 北京交通大学 | Data processing method and device for tunnel driving construction prediction |
| CN115688428A (en) * | 2022-11-01 | 2023-02-03 | 青岛理工大学 | Method, device and equipment for calculating critical value of explosive quantity of blasting near zone supporting structure |
| CN115935753A (en) * | 2022-12-26 | 2023-04-07 | 西南交通大学 | Method for simulating and analyzing influence of small-clear-distance tunnel blasting on adjacent tunnel |
Non-Patent Citations (1)
| Title |
|---|
| 复线隧道开挖爆破振动监测分析与振速预测;周建军;《科学之友》;摘要，第66页左栏第一段-第69页右栏最后一段 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116258285A (en) | 2023-06-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cardu et al. | Analysis of predictor equations for determining the blast-induced vibration in rock blasting | |
| CN114781040B (en) | Building blasting vibration safety evaluation method, device, equipment and medium | |
| Shen et al. | FRACOD modeling of rock fracturing and permeability change in excavation-damaged zones | |
| Yu et al. | Analysis for microseismic energy of immediate rockbursts in deep tunnels with different excavation methods | |
| Shadabfar et al. | Probabilistic modeling of excavation-induced damage depth around rock-excavated tunnels | |
| Wu et al. | Parameter calculation of the initiating circuit with mixed use of nonel detonators and electronic detonators in tunnel controlled-blasting | |
| CN109444463A (en) | A kind of ore-rock moving monitoring system and method | |
| CN116258285B (en) | Porous small-clear-distance tunnel blasting vibration speed prediction method, device, equipment and medium | |
| JP2018100823A (en) | Blasting method | |
| CN115600368A (en) | Deep hole blasting optimization method, device and equipment and readable storage medium | |
| CN107152964A (en) | Blasting vibration monitoring method for deeply-buried underground water supply pipeline in tunnel step method construction | |
| Li et al. | Blasting vibration safety criterion analysis with equivalent elastic boundary: based on accurate loading model | |
| Wang et al. | Experimental Study of Blast‐Induced Vibration Characteristics Based on the Delay‐Time Errors of Detonator | |
| CN116595608B (en) | Blasting and settlement pipeline safety control evaluation method, device, equipment and medium | |
| KR20220153338A (en) | Blasting management system for analysis of vibration and fragmentation caused by blasting | |
| CN114964469B (en) | Method and system for testing blasting vibration of asymmetric high-stress tunnel | |
| CN115688428A (en) | Method, device and equipment for calculating critical value of explosive quantity of blasting near zone supporting structure | |
| CN116243401A (en) | Submarine tunnel gushing water disaster multi-element information monitoring method and system | |
| Behera et al. | A 3D nearfield vibration model for simultaneous blasting of multiple holes in the sedimentary rock formation. | |
| Proprenter et al. | Risk Management in tunneling-a joint approach of all Involved | |
| Li et al. | Modifying Sadov’s Empirical Formula for Group-Hole Blasting in Open-Pit Mine Conditions | |
| CN116242209B (en) | Vibration damping rate calculation method, system and equipment for vibration damping hole and readable storage medium | |
| CN116307045B (en) | Method, system, equipment and medium for calculating vibration speed of building under tunnel blasting | |
| Stephansson et al. | FRACOD modeling of rock fracturing and permeability change in excavation damaged zones | |
| Man et al. | A new excavation technology of blasting combined with TBM |
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 |