CN116432815A - Tunnel blasting earth surface vibration waveform prediction method - Google Patents

Tunnel blasting earth surface vibration waveform prediction method Download PDF

Info

Publication number
CN116432815A
CN116432815A CN202310146284.4A CN202310146284A CN116432815A CN 116432815 A CN116432815 A CN 116432815A CN 202310146284 A CN202310146284 A CN 202310146284A CN 116432815 A CN116432815 A CN 116432815A
Authority
CN
China
Prior art keywords
blasting
vibration
tunnel
explosive
wave
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.)
Pending
Application number
CN202310146284.4A
Other languages
Chinese (zh)
Inventor
戴斌
王志博
沈韬
张学民
韩淼
万国强
马亮亮
曾晓辉
龚涛
王树英
周贤舜
张燕勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming Metro Construction Management Co ltd
Central South University
China Railway No 3 Engineering Group Co Ltd
China Railway Development Investment Group Co Ltd
Fifth Engineering Co Ltd of China Railway No 3 Engineering Group Co Ltd
Original Assignee
Kunming Metro Construction Management Co ltd
Central South University
China Railway Development Investment Group Co Ltd
Fifth Engineering Co Ltd of China Railway No 3 Engineering Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Kunming Metro Construction Management Co ltd, Central South University, China Railway Development Investment Group Co Ltd, Fifth Engineering Co Ltd of China Railway No 3 Engineering Group Co Ltd filed Critical Kunming Metro Construction Management Co ltd
Priority to CN202310146284.4A priority Critical patent/CN116432815A/en
Publication of CN116432815A publication Critical patent/CN116432815A/en
Pending legal-status Critical Current

Links

Images

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
    • 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
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • 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)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • General Engineering & Computer Science (AREA)
  • Marketing (AREA)
  • Data Mining & Analysis (AREA)
  • Tourism & Hospitality (AREA)
  • Mathematical Physics (AREA)
  • General Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Quality & Reliability (AREA)
  • Computational Mathematics (AREA)
  • Game Theory and Decision Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Development Economics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Drilling And Exploitation, And Mining Machines And Methods (AREA)

Abstract

The invention provides a tunnel blasting earth surface vibration waveform prediction method, which comprises the following steps: firstly, deducing a short column explosive blasting vibration waveform function according to a Heelan theory; then constructing a vibration calculation form suitable for tunnel cut tunneling blasting based on the short column explosive package vibration waveform function, and calculating a tunnel cut hole blasting ground surface vibration waveform function according to a superposition principle; simplifying and correcting the group hole blasting waveform function formula to obtain a tunnel group hole blasting surface vibration waveform function; and finally, comparing the obtained data with field actual measurement data to determine specific parameters in a formula, thereby obtaining a tunnel blasting vibration earth surface waveform function curve suitable for the field. The method can comprehensively and accurately predict the waveform time course curve of the whole process of the vibration of the surface of the tunnel blasting, but not only single blasting peak value prediction, can effectively promote the field application of the vibration waveform prediction of the surface of the tunnel blasting, and has important significance for optimizing blasting parameters and ensuring the safety of various facilities and personnel around the blasting area.

Description

Tunnel blasting earth surface vibration waveform prediction method
Technical Field
The invention relates to the technical field of tunnels, in particular to a tunnel blasting earth surface vibration waveform prediction method.
Background
With the continuous promotion of the urban process, urban population grows year by year, so that huge pressure is brought to urban traffic, and urban congestion is effectively relieved, and many cities are being built or subway tunnel projects penetrating through urban downtown areas are to be built. The drilling and blasting method is an economic, effective and rapid excavation method, and becomes one of the main methods of subway construction projects. However, when blasting is performed in urban areas, the blasting may pose a great risk to the building(s) and the underground pipelines, etc., due to the adjacent surrounding buildings and underground buried pipelines, etc. Therefore, a related calculation prediction study needs to be carried out on the earth surface vibration caused during tunnel blasting construction so as to avoid negative effects caused by the blasting vibration.
At present, a Sargassy formula is often adopted for the prediction of the tunnel blasting surface vibration for calculation, and the following formula is adopted. However, the Sargassy formula only considers the influence of the dosage and distance on the blasting vibration, and does not consider the influence of surrounding rock elastic modulus, longitudinal wave velocity, explosive parameters (explosive diameter, detonation velocity) and the like on the blasting vibration, and the blasting vibration velocity is very sensitive to the change of the parameters, so that a larger error exists in the blasting vibration result calculated by the Sargassy formula. Meanwhile, the Sargassy formula only reflects the change of the blasting vibration peak value under single blasting and fails to reflect the content of the frequency and duration under blasting, and the factors are also important for the research of blasting vibration.
Figure BDA0004089226650000021
The background description provided herein is for the purpose of generally presenting the context of the disclosure. Unless otherwise indicated herein, the materials described in this section are not prior art to the claims in this application and are not admitted to be prior art by inclusion in this section.
Disclosure of Invention
Aiming at the technical problems in the related art, the invention provides a tunnel blasting earth surface vibration waveform prediction method which is characterized by comprising the following steps:
s1, deducing a short-column explosive blasting vibration waveform function according to a Heelan theory, wherein the short-column explosive blasting excitation radial and axial vibration speeds are as follows:
Figure BDA0004089226650000022
wherein,,
ω＝2πf
Figure BDA0004089226650000023
the frequency f is considered to be attenuated along with the dosage and the explosive distance; k, alpha is a coefficient and an attenuation index related to the terrain, geological conditions and vibration direction of the blasting point, and Q; r is the distance from the observation point to the center of the medicine bag; v (V) p And V s Longitudinal wave velocity and transverse wave velocity of surrounding rock respectively, wherein subscripts r and z respectively represent directions;
s2, constructing an actual tunnel cut blasting vibration calculation form by using the short-column explosive-bag vibration waveform function, dividing the cut hole into a plurality of short-column explosive bags based on a superposition principle, and performing superposition calculation to obtain a tunnel cut hole blasting earth surface vibration waveform function;
s3, simplifying and correcting the group hole blasting waveform function formula to obtain a tunnel group hole blasting earth surface vibration waveform function;
and S4, verifying the group hole blasting waveform function according to the field actually measured vibration time course curve, determining related parameters in a formula, and constructing a tunnel blasting vibration surface waveform function curve suitable for the field.
Specifically, in step S1, the calculation method of the short column explosive blasting vibration waveform function is as follows:
according to Heelan theory, the displacement of P wave and S wave of the short column cavity under the action of lateral pressure load is solved as follows:
Figure BDA0004089226650000031
Figure BDA0004089226650000032
wherein,,
Figure BDA0004089226650000033
Figure BDA0004089226650000034
Figure BDA0004089226650000035
Figure BDA0004089226650000036
R p 、R s 、Z p z is as follows s R-direction displacement and Z-direction displacement of the P wave and the S wave respectively; v (V) p And V s Propagation speeds of the P wave and the S wave respectively; r is the distance from the observation point to the center of the medicine bag;
Figure BDA0004089226650000041
is the included angle between the wave propagation direction and the negative direction of the z axis; p (t) is the radial pressure exerted by the stub cavity; />
Figure BDA0004089226650000042
And->
Figure BDA0004089226650000043
Source functions for P-wave and S-wave; delta is the volume of the stub cavity; g is the shear modulus of the surrounding rock, E is the elastic modulus of the surrounding rock; ρ is the surrounding rock density; mu (mu)Is the poisson ratio of the surrounding rock.
Summing and deriving the formula (1) to obtain the excitation radial and axial vibration velocity formulas of the short column medicine bag, wherein the formulas are calculated as follows:
Figure BDA0004089226650000044
the maximum peak vibration speed of the measuring point in the radial direction and the axial direction is calculated by the formula (2), meanwhile, the vibration wave is a function gradually decreasing along with time as known from test data, in order to simplify the function form of calculating the vibration waveform, the on-site vibration waveform is considered to fluctuate in a sine wave form, and the vibration formula can be expressed as follows:
Figure BDA0004089226650000045
specifically, the calculation method of the radial pressure P (t) applied to the stub cavity is as follows:
P(t)＝P m e -at+b (4)
wherein,,
Figure BDA0004089226650000046
a. b is the detonation attenuation coefficient of the explosive, and is related to the explosive material, the property parameter and the axial uncoupled coefficient of the blast hole; p (P) m The maximum pressure under the expansion of explosive stress wave and explosive gas is considered; n is the detonation gas pressure increase multiple; ρ e Is the density of the explosive; d (D) e Is the detonation velocity of the explosive; d, d c Is the diameter of the explosive; d, d b Is the diameter of the blast hole.
Specifically, in the step S2, the calculation method of the tunnel cut hole blasting earth surface vibration waveform function is as follows:
because the short column medicine bag is uniformly and outwards diffused in a column wave form, the vibration speed of any point on the wave front in the radial r direction and the axial Z direction can be calculated by adopting the step S1; and then the three-way vibration waveform function of any measuring point position of the ground surface is obtained by converting a cylindrical coordinate system and a rectangular coordinate system, namely decomposing the r-direction vibration speed, and the formula is as follows:
X＝r·sinθ
Y＝r·cosθ
Z＝Z(5)
and finally, superposing the vibration caused by the short column explosive charges to obtain the vibration caused by the actual column explosive charges, wherein the calculation flow of the vibration of the surface of the cut hole blasting is as follows:
firstly, dividing a medicine bag with the length of l into n short column medicine bags with the same size, wherein the mass of each short column medicine bag is the same, and the axial distance from the ith short column medicine bag to a measuring point is as follows:
Figure BDA0004089226650000051
the explosive core distance from the short column medicine bag of the ith section to the measuring point is as follows:
Figure BDA0004089226650000052
then, the time for the vibration wave of the short column medicine bag in the ith section to reach the measuring point is as follows:
Figure BDA0004089226650000053
Figure BDA0004089226650000054
substituting formulas (6) - (8) into formula (3) and superposing to obtain the columnar explosive package blasting measuring point B, wherein the vibration formulas in the directions r and z are as follows:
Figure BDA0004089226650000061
Figure BDA0004089226650000062
and finally decomposing the vibration speed in the r direction along the x and y axes to obtain the horizontal and vertical vibration speeds of the measuring point B, wherein the vibration speed is as follows:
V x (t)＝V r (t)·sinθ
V y (t)＝V r (t)·cosθ (10)
wherein θ is the included angle between the connecting line of the measuring point and the medicine bag and the vertical line of the ground surface
Specifically, the simplified calculation method of the tunnel group hole blasting earth surface vibration waveform function in the step S3 is as follows:
when the dosage and the volume of blast holes are consistent and the geological conditions experienced by earthquake waves formed by explosion are the same, each single-hole blasting vibration function is only related to the blasting center distance R, the blasting center distance from the earth surface measuring point to each blast hole is far greater than the hole distance between the cut holes, the change of the blasting center distances caused by the arrangement positions of the cut holes is calculated to be smaller, the vibration velocity waveforms formed by the blast holes are negligible, the vibration velocity waveforms formed by the blast holes are consistent, the group hole blasting function is only the reproduction of single-hole blasting in time, and when the interval time of the blast holes is the same, the vibration waveform functions of the tunnel group Kong Weicha blasting vibration are simplified to be obtained:
Figure BDA0004089226650000063
wherein T is the interval time of blasting hole initiation; n is the total number of blastholes detonated once, and m is the sequence of blasthole detonating.
Specifically, the method further comprises the following steps: s5, predicting blasting vibration peak values and frequency changes of all buildings around the tunnel according to the obtained tunnel blasting earth surface vibration waveform function, and adjusting blasting parameters and schemes.
The tunnel blasting earth surface vibration waveform prediction method provided by the invention has the following beneficial effects:
(1) Consideration of more comprehensive
The tunnel blasting earth surface vibration waveform prediction method fully considers factors such as surrounding rock elastic modulus, shear modulus, explosive density, detonation velocity, diameter and the like, and can more accurately reflect the actual situation of blasting vibration.
(2) Can more truly reflect the blasting vibration signal
The tunnel blasting earth surface vibration waveform prediction method can directly predict the time course curve of blasting vibration instead of single vibration peak value prediction, can more truly evaluate the influence of peak value, frequency and the like caused by the blasting vibration, and furthest reduces the damage of the blasting vibration.
(3) The monitoring workload is small
When factors such as blasting circulation footage, surrounding rock level, distance and the like are changed, the blasting vibration function can be obtained by only adjusting related parameters in a theoretical calculation formula, and the field multiple tests are not needed.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for predicting vibration waveforms of tunnel blasting ground surfaces, which is provided by an embodiment of the invention;
FIG. 2 is a schematic diagram of the relative position of a measurement point to a tunnel;
FIG. 3 is a schematic diagram of the relative positions of the tunnel bags and the measuring points;
FIG. 4 is a schematic diagram of the model calculation of the present embodiment;
fig. 5 is a comparison of the calculation and actual measurement of tunnel group hole blasting in the present embodiment.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which are derived by a person skilled in the art based on the embodiments of the invention, fall within the scope of protection of the invention.
Example 1
Referring to fig. 1, the embodiment discloses a tunnel blasting earth surface vibration waveform prediction method, which comprises the following steps:
s1, deducing a short column explosive blasting vibration waveform function according to a Heelan theory;
according to Heelan theory, the displacement of P wave and S wave of the short column cavity under the action of lateral pressure load is solved as follows:
Figure BDA0004089226650000081
Figure BDA0004089226650000082
wherein,,
Figure BDA0004089226650000091
Figure BDA0004089226650000092
Figure BDA0004089226650000093
Figure BDA0004089226650000094
R p 、R s 、Z p z is as follows s R-direction displacement and Z-direction displacement of the P wave and the S wave respectively; v (V) p And V s Longitudinal wave velocity and transverse wave velocity of surrounding rock respectively; r is the distance from the observation point to the center of the medicine bag;
Figure BDA0004089226650000095
is the negative of the wave propagation direction and the z-axisAn included angle in the direction; p (t) is the radial pressure exerted by the stub cavity; />
Figure BDA0004089226650000096
And->
Figure BDA0004089226650000097
Source functions for P-wave and S-wave; delta is the volume of the stub cavity; g is the shear modulus of the surrounding rock, E is the elastic modulus of the surrounding rock; ρ is the surrounding rock density; mu is the surrounding rock Poisson's ratio.
Summing and deriving the formula (1) to obtain the excitation radial and axial vibration velocity formulas of the short column medicine bag, wherein the formulas are calculated as follows:
Figure BDA0004089226650000098
the maximum peak vibration speed of the measuring point in the radial direction and the axial direction can be calculated according to the formula (2), meanwhile, the vibration wave is a function gradually decreasing along with time as known from test data, in order to simplify the function form of calculating the vibration waveform, the on-site vibration waveform is considered to fluctuate according to the sine wave form, and the vibration formula can be expressed as follows:
Figure BDA0004089226650000099
wherein,,
ω＝2πf
Figure BDA0004089226650000101
the frequency f is considered to be attenuated along with the dosage and the explosive distance; k and α are coefficients and damping indexes related to the terrain, geological conditions and vibration direction of the blasting point, and in general, α is 1 to 2, k=30 to 70 is taken in the rock, k=150 to 250 is taken in the soil, and the harder the rock body, the smaller the values of k and α are. Wherein the subscripts r, z denote the r, z directions, e.g. V, respectively, in the RTZ coordinate system pr Representing the component of the longitudinal wave velocity of the surrounding rock in the r direction, and others are similar.
The calculation method of the radial pressure P (t) applied to the short column cavity is as follows:
P(t)＝P m e -at+b (4)
wherein,,
Figure BDA0004089226650000102
a. b is the detonation attenuation coefficient of the explosive, and is related to the explosive material, the property parameter and the axial uncoupled coefficient of the blast hole; p (P) m The maximum pressure under the expansion of explosive stress wave and explosive gas is considered; n is the pressure increase multiple of detonation gas, and 8-11 is taken; ρ e Is the density of the explosive; d (D) e Is the detonation velocity of the explosive; d, d c Is the diameter of the explosive; d, d b Is the diameter of the blast hole.
S2, constructing an actual tunnel cut blasting vibration calculation form by using the short-column explosive-bag vibration waveform function, dividing the cut hole into a plurality of short-column explosive bags based on a superposition principle, and performing superposition calculation to obtain a tunnel cut hole blasting earth surface vibration waveform function;
since the short column medicine bag is uniformly and outwards diffused in the form of a column wave, the vibration speed of any point on the wave front in the radial r direction and the axial Z direction can be calculated by adopting the content in the step 1, as shown in figure 1. And then the r-direction vibration speed can be decomposed through the transformation of the cylindrical coordinate system and the rectangular coordinate system, so that a three-way vibration waveform function of any measuring point position of the ground surface can be obtained, and the formula is as follows:
X＝r·sinθ
Y＝r·cosθ
Z＝Z(5)
since the detonation velocity of the columnar explosive is usually between 2000 and 7000m/s in field practice, the lateral load P (t) cannot be applied to the wall of the blast hole at the same time, in order to consider the deviation of the vibration velocity calculation caused by the detonation velocity, the columnar explosive is decomposed into a plurality of columnar explosive, the blasting vibration waveform of the columnar explosive is considered to be formed by sequential blasting of the plurality of columnar explosive, and finally the vibration caused by the plurality of columnar explosive is superimposed, namely the vibration caused by the actual columnar explosive, and the relative positions of the tunnel explosive and the measuring point are shown in fig. 2. The calculation flow of the vibration of the ground surface of the cut hole blasting is as follows:
firstly, dividing a medicine bag with the length of l into n short column medicine bags with the same size, wherein the mass of each short column medicine bag is the same, and as can be seen from fig. 2, the axial distance from the ith short column medicine bag to a measuring point is as follows:
Figure BDA0004089226650000111
the explosive core distance from the short column medicine bag of the ith section to the measuring point is as follows:
Figure BDA0004089226650000112
then, the time for the vibration wave of the short column medicine bag in the ith section to reach the measuring point is as follows:
Figure BDA0004089226650000113
Figure BDA0004089226650000114
substituting formulas (6) - (8) into formula (3) and superposing to obtain the columnar explosive package blasting measuring point B, wherein the vibration formulas in the directions r and z are as follows:
Figure BDA0004089226650000121
Figure BDA0004089226650000122
and finally decomposing the vibration speed in the r direction along the x and y axes to obtain the horizontal and vertical vibration speeds of the measuring point B, wherein the vibration speed is as follows:
V x (t)＝V r (t)·sinθ
V y (t)＝V r (t)·cosθ(10)
s3, simplifying and correcting the group hole blasting waveform function formula to obtain a tunnel group hole blasting earth surface vibration waveform function;
when the explosive quantity and the volume of the blastholes are consistent and the geological conditions experienced by the seismic waves formed by explosion are the same, each single-hole blasting vibration function is only related to the explosive distance R. Under normal conditions, the explosion distance from the earth surface measuring point to each blast hole is far greater than the hole distance between the cut holes, and the change of the explosion distance caused by the arrangement position of the cut holes is calculated to be small and negligible. Therefore, the vibration velocity waveforms formed by the blast holes are consistent, the group hole blasting function is only the reproduction of single hole blasting in time, and when the blast hole interval time is the same, the tunnel group Kong Weicha blasting vibration waveform function is simplified to be obtained as follows:
Figure BDA0004089226650000123
wherein T is the interval time of blasting hole initiation; n is the total number of blastholes detonated once, and m is the sequence of blasthole detonating.
And S4, verifying the group hole blasting waveform function according to the field actually measured vibration time course curve, determining related parameters in a formula, and constructing a tunnel blasting vibration surface waveform function curve suitable for the field.
Further, in the embodiment, the blasting vibration peak value and the frequency change of each building (structure) around the tunnel are predicted according to the obtained tunnel blasting earth surface vibration waveform function, so that blasting parameters and schemes can be adjusted in time, and the damage of tunnel blasting to the surrounding building (structure) is reduced.
Namely, the method also comprises the following steps:
s5, predicting blasting vibration peak values and frequency changes of all buildings around the tunnel according to the obtained tunnel blasting earth surface vibration waveform function, and adjusting blasting parameters and schemes.
The implementation relies on an in-out line undercut tunnel section field test of a certain urban rail transit line, double holes are arranged for simultaneous detonation, a second emulsion explosive is adopted, the explosive parameters are shown in table 1, the diameter of the blast hole is 42mm, the diameter of the explosive is 32mm, the diameter of the single Kong Yaoliang is 0.6kg, the depth of the blast hole is 1.0m, and the hole opening is plugged by stemming with the length of 0.4m. And arranging an instrument on the surface axis at a position 12.1m away from the face for blasting vibration monitoring, wherein the surrounding rock mechanical parameters are shown in Table 2.
Table 1 table of explosive parameters
Figure BDA0004089226650000131
TABLE 2 mechanical parameters table of surrounding rock
Figure BDA0004089226650000132
1): the tunnel cut package is divided, and the model calculation diagram is shown in fig. 4. The horizontal axial distance from the measuring point to the explosive core is 12.1m, namely z=12.1, the height h=10m from the measuring point to the explosive source is namely y=10m, and the measuring point is arranged in the axial direction of the tunnel, so that x=0m. According to the principle of linear superposition division, the columnar medicine bag is divided into superposition of 8 short column medicine bags.
The longitudinal wave speed and the transverse wave speed of the surrounding rock can be obtained by the Heelan theory:
Figure BDA0004089226650000141
Figure BDA0004089226650000142
taking detonation gas pressure increase multiple n=8, then obtaining the lateral pressure peak value of each short column medicine bag according to the formula 4 as follows:
Figure BDA0004089226650000143
the axial distance from each short column of medicine bag to the measuring point is obtained by the formula (6):
Figure BDA0004089226650000144
the explosive core distance of each short column of the medicine bag is obtained by the formula (7):
Figure BDA0004089226650000145
the time for each short column of the medicine package vibration wave to reach the measuring point is obtained by the formula (8):
Figure BDA0004089226650000146
Figure BDA0004089226650000147
the blasting field is IV-level surrounding rock, k=60, alpha=1.39 is taken, and the attenuation frequency of the measuring point is obtained as follows:
Figure BDA0004089226650000151
so the angular frequency ω=2pi f= 269.3
And substituting the calculation result into the formula (9) to solve the surface vibration waveform function of the single-hole blasting.
2) Substituting the calculated surface vibration waveform function of single hole blasting into the formula (11) to obtain the surface vibration waveform function of group hole blasting, and simultaneously performing simultaneous blasting on two holes, so that the differential time T=0 ms.
3) And comparing the calculated vibration waveform of the blasting earth surface of the tunnel group hole with the measured data, and drawing the vibration waveform as shown in figure 4. The graph shows that the waveform function obtained by theoretical calculation is very consistent with the actual measurement function in vibration speed peak value and vibration form, so that the accuracy of the vibration waveform function of the tunnel group hole columnar medicine bag is verified. The calculation result of the formula can more accurately predict the surface vibration waveform.
Compared with the prior art, the tunnel blasting earth surface vibration waveform prediction method provided by the embodiment has the following beneficial effects:
(1) Consideration of more comprehensive
The tunnel blasting earth surface vibration waveform prediction method provided by the embodiment fully considers factors such as surrounding rock elastic modulus, shear modulus, explosive density, explosion velocity, diameter and the like, and can more accurately reflect the actual situation of blasting vibration.
(2) Can more truly reflect the blasting vibration signal
The tunnel blasting earth surface vibration waveform prediction method provided by the embodiment can directly predict the time course curve of blasting vibration, but not single vibration peak value prediction, can more truly evaluate the influence of peak value, frequency and the like caused by the blasting vibration, and furthest reduces the damage of the blasting vibration.
(3) The monitoring workload is small
When factors such as blasting circulation footage, surrounding rock level, distance and the like are changed, the blasting vibration function can be obtained by only adjusting related parameters in a theoretical calculation formula, and the field multiple tests are not needed.
The foregoing description of the preferred embodiments of the invention is not intended to be limiting, but rather is intended to cover all modifications, equivalents, alternatives, and improvements that fall within the spirit and scope of the invention.

Claims (6)

1. The tunnel blasting earth surface vibration waveform prediction method is characterized by comprising the following steps of:
s1, deducing a short-column explosive blasting vibration waveform function according to a Heelan theory, wherein the short-column explosive blasting excitation radial and axial vibration speeds are as follows:
Figure FDA0004089226640000011
wherein,,
ω＝2πf
Figure FDA0004089226640000012
the frequency f is considered to be attenuated along with the dosage and the explosive distance; k, alpha is a coefficient and an attenuation index related to the terrain, geological conditions and vibration direction of the blasting point, and Q; r is the distance from the observation point to the center of the medicine bag; v (V) p And V s Longitudinal wave velocity and transverse wave velocity of surrounding rock respectively, wherein subscripts r and z respectively represent directions;
s2, constructing an actual tunnel cut blasting vibration calculation form by using the short-column explosive-bag vibration waveform function, dividing the cut hole into a plurality of short-column explosive bags based on a superposition principle, and performing superposition calculation to obtain a tunnel cut hole blasting earth surface vibration waveform function;
s3, simplifying and correcting the group hole blasting waveform function formula to obtain a tunnel group hole blasting earth surface vibration waveform function;
and S4, verifying the group hole blasting waveform function according to the field actually measured vibration time course curve, determining related parameters in a formula, and constructing a tunnel blasting vibration surface waveform function curve suitable for the field.
2. The tunnel blasting earth's surface vibration waveform prediction method according to claim 1, characterized by comprising the steps of: the calculation method of the short cartridge blasting vibration waveform function in the step S1 is as follows:
according to Heelan theory, the displacement of P wave and S wave of the short column cavity under the action of lateral pressure load is solved as follows:
Figure FDA0004089226640000021
Figure FDA0004089226640000022
wherein,,
Figure FDA0004089226640000023
Figure FDA0004089226640000024
Figure FDA0004089226640000025
Figure FDA0004089226640000026
R p 、R s 、Z p z is as follows s R-direction displacement and Z-direction displacement of the P wave and the S wave respectively; v (V) p And V s Propagation speeds of the P wave and the S wave respectively; r is the distance from the observation point to the center of the medicine bag;
Figure FDA0004089226640000027
is the included angle between the wave propagation direction and the negative direction of the z axis; p (t) is the radial pressure exerted by the stub cavity; />
Figure FDA0004089226640000028
And->
Figure FDA0004089226640000029
Source functions for P-wave and S-wave; delta is the volume of the stub cavity; g is the shear modulus of the surrounding rock, E is the elastic modulus of the surrounding rock; ρ is the surrounding rock density; mu is the Poisson's ratio of the surrounding rock,
summing and deriving the formula (1) to obtain the excitation radial and axial vibration velocity formulas of the short column medicine bag, wherein the formulas are calculated as follows:
Figure FDA00040892266400000210
the maximum peak vibration speed of the measuring point in the radial direction and the axial direction is calculated by the formula (2), meanwhile, the vibration wave is a function gradually decreasing along with time as known from test data, in order to simplify the function form of calculating the vibration waveform, the on-site vibration waveform is considered to fluctuate in a sine wave form, and the vibration formula can be expressed as follows:
Figure FDA0004089226640000031
3. the tunnel blasting earth' S surface vibration waveform prediction method according to claim S2, wherein: the calculation method of the radial pressure P (t) applied to the stub cavity is as follows:
P(t)＝P m e -at+b (4)
wherein,,
Figure FDA0004089226640000032
a. b is the detonation attenuation coefficient of the explosive, and is related to the explosive material, the property parameter and the axial uncoupled coefficient of the blast hole; p (P) m The maximum pressure under the expansion of explosive stress wave and explosive gas is considered; n is the detonation gas pressure increase multiple; ρ e Is the density of the explosive; d (D) e Is the detonation velocity of the explosive; d, d c Is the diameter of the explosive; d, d b Is the diameter of the blast hole.
4. The tunnel blasting earth's surface vibration waveform prediction method according to claim 1, characterized by comprising the steps of: the calculation method of the tunnel cut hole blasting earth surface vibration waveform function in the step S2 is as follows:
because the short column medicine bag is uniformly and outwards diffused in a column wave form, the vibration speed of any point on the wave front in the radial r direction and the axial Z direction can be calculated by adopting the step S1; and then the three-way vibration waveform function of any measuring point position of the ground surface is obtained by converting a cylindrical coordinate system and a rectangular coordinate system, namely decomposing the r-direction vibration speed, and the formula is as follows:
X＝r·sinθ
Y＝r·cosθ
Z＝Z (5)
and finally, superposing the vibration caused by the short column explosive charges to obtain the vibration caused by the actual column explosive charges, wherein the calculation flow of the vibration of the surface of the cut hole blasting is as follows:
firstly, dividing a medicine bag with the length of l into n short column medicine bags with the same size, wherein the mass of each short column medicine bag is the same, and the axial distance from the ith short column medicine bag to a measuring point is as follows:
Figure FDA0004089226640000041
the explosive core distance from the short column medicine bag of the ith section to the measuring point is as follows:
Figure FDA0004089226640000042
then, the time for the vibration wave of the short column medicine bag in the ith section to reach the measuring point is as follows:
Figure FDA0004089226640000043
Figure FDA0004089226640000044
substituting formulas (6) - (8) into formula (3) and superposing to obtain the columnar explosive package blasting measuring point B, wherein the vibration formulas in the directions r and z are as follows:
Figure FDA0004089226640000045
Figure FDA0004089226640000046
and finally decomposing the vibration speed in the r direction along the x and y axes to obtain the horizontal and vertical vibration speeds of the measuring point B, wherein the vibration speed is as follows:
V x (t)＝V r (t)·sinθ
V y (t)＝V r (t)·cosθ (10)
and θ is the included angle between the connecting line of the measuring point and the medicine bag and the vertical line of the ground surface.
5. The tunnel blasting earth's surface vibration waveform prediction method according to claim 1, characterized by comprising the steps of: the simplified calculation method of the tunnel group hole blasting earth surface vibration waveform function in the step S3 is as follows:
when the dosage and the volume of blast holes are consistent and the geological conditions experienced by earthquake waves formed by explosion are the same, each single-hole blasting vibration function is only related to the blasting center distance R, the blasting center distance from the earth surface measuring point to each blast hole is far greater than the hole distance between the cut holes, the change of the blasting center distances caused by the arrangement positions of the cut holes is calculated to be smaller, the vibration velocity waveforms formed by the blast holes are negligible, the vibration velocity waveforms formed by the blast holes are consistent, the group hole blasting function is only the reproduction of single-hole blasting in time, and when the interval time of the blast holes is the same, the vibration waveform functions of the tunnel group Kong Weicha blasting vibration are simplified to be obtained:
Figure FDA0004089226640000051
wherein T is the interval time of blasting hole initiation; n is the total number of blastholes detonated once, and m is the sequence of blasthole detonating.
6. The method of claim 1, the method further comprising: s5, predicting blasting vibration peak values and frequency changes of all buildings around the tunnel according to the obtained tunnel blasting earth surface vibration waveform function, and adjusting blasting parameters and schemes.
CN202310146284.4A 2023-02-22 2023-02-22 Tunnel blasting earth surface vibration waveform prediction method Pending CN116432815A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310146284.4A CN116432815A (en) 2023-02-22 2023-02-22 Tunnel blasting earth surface vibration waveform prediction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310146284.4A CN116432815A (en) 2023-02-22 2023-02-22 Tunnel blasting earth surface vibration waveform prediction method

Publications (1)

Publication Number Publication Date
CN116432815A true CN116432815A (en) 2023-07-14

Family

ID=87082183

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310146284.4A Pending CN116432815A (en) 2023-02-22 2023-02-22 Tunnel blasting earth surface vibration waveform prediction method

Country Status (1)

Country Link
CN (1) CN116432815A (en)

Similar Documents

Publication Publication Date Title
St John et al. Aseismic design of underground structures
CN107066753A (en) The Forecasting Methodology of shallow tunnel blast working earth's surface particle peak value vibration velocity
CN109470101B (en) Method for researching influence of adjacent tunnel blasting on existing tunnel based on field model
CN111241679B (en) Tunnel blasting parameter design method based on digital electronic detonator detonation
Fan et al. Advanced stability analysis of the tunnels in jointed rock mass based on TSP and DEM
CN103399342A (en) Method for forecasting transient unloading induced vibrations based on strain energy of rock masses
CN113552629A (en) Tunnel surrounding rock longitudinal wave velocity determination method and device and computer equipment
Anderson Blast monitoring: regulations, methods and control techniques
CN109612356B (en) Method for acquiring dynamic response characteristic of adjacent concrete pipeline during foundation pit excavation blasting
CN117371111B (en) TBM card machine prediction system and method based on deep neural network and numerical simulation
CN108592720B (en) Vibration damping hole arrangement method for reducing ground blasting vibration
McGarr et al. Coseismic and other short-term strain changes recorded with Sacks-Evertson strainmeters in a deep mine, South Africa
US11920472B2 (en) Reasonable millisecond time control method for excavation blasting of tunnel
Xu et al. Optimal design of microseismic monitoring networking and error analysis of seismic source location for rock slope
Khokhlov et al. Conducting industrial explosions near gas pipelines
DOWDING III Response of buildings to ground vibrations resulting from construction blasting
CN116432815A (en) Tunnel blasting earth surface vibration waveform prediction method
Yamamoto et al. Experimental and theoretical study on smooth blasting with electronic delay detonators
Jiang et al. Safety assessment of upper buried gas pipeline under blasting vibration of subway tunnel: A case study in Beijing subway line
Hua-you et al. Analysis of characteristics of compound vibration and effects to surrounding gas pipeline caused by impact and explosion
Gran et al. Dynamic response of an opening in jointed rock
Bull Soil-Structure interaction: numerical analysis and modelling
CN111365072B (en) Guide hole exploration structure and method for huge cavern of fractured rock mass
CN115030235A (en) Model test device for researching failure of pile foundation under blasting effect
CN110427688B (en) Ground stress prediction method based on actual measurement vibration

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
TA01 Transfer of patent application right

Effective date of registration: 20230811

Address after: China Railway Building, Caiyun South Road, Hangong District, Kunming City, Yunnan Province, 650200

Applicant after: CHINA RAILWAY DEVELOPMENT INVESTMENT GROUP Co.,Ltd.

Applicant after: CENTRAL SOUTH University

Applicant after: CHINA RAILWAY NO.3 ENGINEERING GROUP Co.,Ltd.

Applicant after: China Railway No.3 Engineering Group No.5 Engineering Co.,Ltd.

Applicant after: Kunming Metro Construction Management Co.,Ltd.

Address before: China Railway Building, Caiyun South Road, Chenggong District, Kunming City, Yunnan Province 650200

Applicant before: CHINA RAILWAY DEVELOPMENT INVESTMENT GROUP Co.,Ltd.

Applicant before: CENTRAL SOUTH University

Applicant before: China Railway No.3 Engineering Group No.5 Engineering Co.,Ltd.

Applicant before: Kunming Metro Construction Management Co.,Ltd.

TA01 Transfer of patent application right