阅读说明：本技术 一种基于试爆数据预测水下爆破装药量的方法 (Method for predicting underwater blasting explosive loading based on test blasting data ) 是由 顾云 孙飞 李飞 刘迪 刘勤杰 徐静 于 2021-08-09 设计创作，主要内容包括：本发明公开了一种基于试爆数据预测水下爆破装药量的方法,首先通过试爆数据代入公式拟合出待定参数K、β-(1)和β-(2)的数值；再根据水下爆破近邻域受保护建(构)筑物特征,参考《爆破安全规程》(GB6722-2014)查出各频率范围许可爆破震动速度V；再由公式推算出爆破齐发炸药量Q-(max)；由于最大许可爆破震动速度V与爆破震动主频相关联,从而得到Q-( max)’。本发明预测方法简单,通过实验测试以观察爆破效果及爆破危害,并通过试爆结果推算施工方案的爆破震动危害,让试爆效果观察数据与爆破施工方案预期效果之间具有的科学关联,有据可依。(The invention discloses a method for predicting underwater blasting explosive loading quantity based on test blasting data K 、 β 1 And β 2 the value of (d); then according to the characteristics of the protected building (structure) in the near field of underwater blasting, the allowable blasting vibration speed in each frequency range is found out by referring to blasting safety regulations (GB6722-2014) V (ii) a Then calculating the explosive quantity of blasting explosive by formula Q max (ii) a Due to maximum allowable blasting vibration speed V Is associated with the blast vibration dominant frequency to obtain Q max ’。 The prediction method is simple, the blasting effect and the blasting hazard are observed through experimental tests, the blasting vibration hazard of the construction scheme is calculated through the trial blasting result, and scientific association between the trial blasting effect observation data and the expected effect of the blasting construction scheme is provided.)

1. A method for predicting underwater blasting explosive loading based on test blasting data is characterized by comprising the following steps:

step 1, firstly carrying out test blasting, and carrying out blasting detection on the underwater blasting operation adjacent domain embankment protection object by using a blasting vibration meter so as to obtain vibration intensity V, blasting explosive quantity Q, spatial distance R between a monitoring point and a blasting point and blasting vibration water area propagation distance R_D；

Step 2, comparing the data V, Q, R and R obtained in the step 1_DSubstituting the parameters into the following prediction model formula (1), and obtaining undetermined parameters K and beta in the prediction model formula (1) by adopting linear least square regression fitting₁And beta₂The numerical value of (A):

wherein K represents a site coefficient, beta₁Expressed is the attenuation coefficient, beta₂The attenuation coefficient caused by the change of the soil-water medium is shown;

step 3, obtaining the parameters K and beta in the step 2₁And beta₂Substituted into the following calculation model equation (2):

obtaining the relation between the charge Q and V, and obtaining the maximum value V of V in the protection category according to the table look-up_maxWill V_maxSubstituting into formula (2) to obtain maximum blasting explosive quantity Q_max。

2. The method for predicting underwater blasting explosive loading based on test blasting data as claimed in claim 1, characterized by comprising the following steps:

step one, because the maximum allowable blasting vibration speed V is related to the blasting vibration dominant frequency, Q obtained in the step 3 is used_maxSubstituting the following formula (3) to obtain f_max:

Step two, since f is known_maxFrom the table lookup, V in the protection class can be known_max', will V_max' substitution into formula (2) to obtain Q_max’。

Technical Field

The invention relates to a method for predicting underwater blasting explosive loading based on test blasting data, and belongs to the field of engineering blasting.

Background

When the explosive explodes in a medium, a part of energy is necessarily converted into seismic waves and rapidly spreads from an explosion source to a surrounding medium, so that the surface and the building (structure) are vibrated, and the building (structure) which is not a blasting target nearby can be damaged to different degrees, and the phenomenon is called blasting vibration effect. With the rapid development of national economic construction, blasting operation is widely applied, and the harm brought by corresponding blasting vibration effect draws high attention of people.

However, in engineering practice, the evaluation of the vibration hazard that may be brought by a blasting engineering technician to a blasting construction scheme is mainly determined by the following two measures: (1) an empirical formula derived from Newton mechanics by a scholars Sadow-fusi is citedPerforming accounting, wherein the coefficient K and the attenuation index alpha related to the topographic and geological conditions from the explosion point to the protected object in the accounting formula are in a large range, and the selection of parameter values during accounting completely depends on the personal experience of technicians and lacks of clear basis, so that the rationality of the accounting result is difficult to guarantee; (2) the experimental test is carried out through the trial explosion to observe the blasting effect and the blasting hazard, but the blast vibration hazard of the construction scheme is calculated through the trial explosion result and completely depends on personal experience, and scientific association between the trial explosion effect observation data and the expected effect of the blasting construction scheme is lacked, so that the scientific and reasonable design is difficult to guarantee according to the trial explosion test data fed back to the blasting construction scheme.

In addition, the two measures are widely applied to the engineering blasting industry due to simplicity and easiness in operation, the application effect in the blasting demolition field is good, and the great confusion is brought to the underwater blasting field. In particular, the blasting operation often faces a dilemma that is difficult to implement due to the loss and dispute caused by the blasting vibration hazard, and the normal engineering construction is greatly influenced.

Disclosure of Invention

The purpose of the invention is as follows: in order to overcome the defects in the prior art, the invention provides a method for predicting underwater blasting explosive loading based on test blasting data.

The technical scheme is as follows: in order to solve the technical problem, the method for predicting underwater blasting explosive loading quantity based on test blasting data comprises the steps of firstly carrying out test blasting to obtain vibration data of the test blasting; the vibration data comprises blasting vibration intensity and blasting vibration frequency; then combining with a test explosion scheme, reversely deducing explosion vibration propagation behavior and an attenuation rule based on an analysis method of the explosion vibration propagation behavior so as to obtain an explosion vibration intensity and an explosion vibration dominant frequency attenuation formula between an explosion point and a monitoring point, and specifically comprising the following steps of:

step 1, firstly carrying out test blasting, and carrying out blasting detection on a building protection object of an underwater blasting near-neighborhood embankment building (structure) by using a blasting vibration meter so as to obtain vibration intensity V, blasting explosive quantity Q, spatial distance R between a monitoring point and a blasting point and blasting vibration water area propagation distance R_D；

Step 2, mixing the data V, Q, R obtained in step 1 with R_DSubstituting the parameters into the following prediction model formula (1), and obtaining undetermined parameters K and beta in the prediction model formula (1) by adopting linear least square regression fitting₁And beta₂The numerical value of (A):

wherein K represents a site coefficient, beta₁Expressed is the attenuation coefficient, beta₂The attenuation coefficient caused by the change of the soil-water medium is shown;

step 3, according to the characteristics of the protected object of the bank building/structure in the near field of underwater blasting, consulting the maximum value V of V in the protection category in blasting safety regulations (GB6722-2014)_maxWill V_maxSubstituting into formula (2) to obtain the maximumLarge blasting explosive quantity Q_maxAnd guiding the engineering blasting:

further, step one, since the maximum allowable blasting vibration speed V is related to the blasting vibration dominant frequency, Q obtained in step 3 is obtained_maxSubstituting the following formula (3) to obtain f_max:

Step two, since f is known_maxFrom the table lookup, V in the protection class can be known_max', will V_max' substitution into formula (2) to obtain Q_max’。

Comparison Q_maxAnd Q_max' selecting the minimum value as the maximum blasting explosive quantity of the blasting operation.

The maximum blasting explosive quantity is obtained in the steps 1-3, and the more accurate maximum blasting explosive quantity can be further calculated through further calculation in the steps one and two.

And (3) deriving the dominant frequency f of the underwater blasting vibration transmitted to the protected object of the building (structure) on the bank of the adjacent domain from the blasting simultaneous firing quantity Q calculated by the allowable blasting vibration speed. By consulting the blasting safety regulations (GB6722-2014), if the dominant frequency value f and the allowable blasting vibration speed V correspond to the same frequency range in the blasting safety regulations (GB6722-2014), the allowable blasting vibration speed V is reasonably selected; correspondingly, the calculated blasting simultaneous firing quantity Q is credible, and a blasting parameter scheme design logic is completed.

As can be seen from the above, the expressions (1), (2) and (3) form a closed-loop computation logic. Firstly, fitting undetermined parameters K and beta by substituting trial explosion data into a formula (1)₁And beta₂The value of (d); then, according to the characteristics of the protected building (structure) in the near field of underwater blasting, referring to blasting safety regulations (GB6722-2014)Blasting vibration speed V is allowed in each frequency range; calculating the explosive quantity Q of the blasting explosive according to the formula (2); further, substituting the blasting explosive quantity Q into the formula (3) to derive the corresponding blasting dominant frequency f; if the dominant frequency f and the allowable blasting vibration speed V correspond to the same frequency range in blasting safety regulations (GB6722-2014), a blasting parameter scheme method is completed. If the explosive quantity is required to be further accurately blasted, the main frequency f is used_maxCalculating to obtain Q_max', comparison Q_maxAnd Q_max' selecting the minimum value as the maximum blasting explosive quantity of the blasting operation. Step one and step two can be seen to further check the calculation of step 1 to step 3 for more accuracy.

The vibration excitation of the underwater blasting on the structure near the bank has two significant differences compared with the vibration excitation of the geotechnical blasting on the structure near the bank: firstly, because the water medium in a fluid state can not bear shear load and can not transmit transverse waves (also called shear waves), the transverse wave components are naturally filtered out by the vibration excitation of the underwater explosion on the structure close to the embankment, and the transverse wave components are mainly shown as longitudinal wave excitation; secondly, due to the difference of wave impedance of the water medium and the rock-soil medium, longitudinal waves from underwater blasting are incident, reflected and transmitted on a water-soil interface, and the specific incident, reflected and transmitted wave intensities are related to the wave impedance of the water and rock-soil medium. In general, compared with rock-soil blasting, the vibration excitation behavior of the underwater blasting on the structure near the bank shows obvious difference on the propagation behavior of waves on an aqueous medium and an interface between the aqueous medium and the rock-soil.

Has the advantages that: (1) due to the difference of media, the method aims at the vibration excitation of the near-neighborhood structure of the embankment caused by the underwater explosion, realizes that the predicted vibration excitation behavior of the near-neighborhood structure of the embankment caused by the underwater explosion cannot directly refer to a rock-soil explosion vibration check formula widely adopted in the current engineering practice due to the load response behavior and the wave impedance difference of water and soil media, needs to take account that the water media cannot transmit transverse waves, and the incident, reflected and transmitted intensities of the longitudinal waves at a water-soil interface are redistributed. Therefore, an operator is supplemented on the basis of a rock-soil blasting vibration prediction formula through dimensional analysis and least square regression fitting, the operator is shown as the third item on the right side of the formula (1) to represent the difference of water and soil media, and the difference is integrated into the prediction formula, so that the underwater blasting vibration correspondingly improves the dam near-neighborhood structure vibration excitation prediction formula; (2) according to the invention, through further checking calculation, more accurate blasting explosive quantity is obtained, so that safer guidance is provided for underwater engineering blasting; (3) the prediction method is simple, the blasting effect and the blasting hazard are observed through experimental tests, the blasting vibration hazard of the construction scheme is calculated through the trial blasting result, and scientific association between the trial blasting effect observation data and the expected effect of the blasting construction scheme is provided and can be relied on; (4) the method is specially improved aiming at underwater blasting, and the excitation behavior of the underwater blasting on the vibration of the structure in the vicinity of the embankment is obtained through theoretical analysis, so that a foundation is provided for subsequent research.

Drawings

FIG. 1 is a logic block diagram of a blast construction scheme predicted based on test blasting data in accordance with the present invention.

Detailed Description

The present invention will be further described with reference to the accompanying drawings.

Examples

The invention discloses a method for predicting underwater blasting explosive loading quantity based on test blasting data, which comprises the steps of firstly carrying out test blasting to obtain test blasting vibration data; the vibration data comprises blasting vibration strength and blasting vibration frequency; then combining with a test explosion scheme, reversely deducing explosion vibration propagation behavior and an attenuation rule based on an analysis method of the explosion vibration propagation behavior so as to obtain an explosion vibration intensity and an explosion vibration dominant frequency attenuation formula between an explosion point and a monitoring point, and specifically comprising the following steps of:

step 1, firstly carrying out test blasting, and carrying out blasting detection on a building protection object of an underwater blasting near-neighborhood embankment building (structure) by using a blasting vibration meter so as to obtain vibration intensity V, blasting explosive quantity Q, spatial distance R between a monitoring point and a blasting point and blasting vibration water area propagation distance R_D；

Step 2, mixing the data V, Q, R obtained in step 1 with R_DSubstituting into the following prediction model equation (1) by using linear least square methodObtaining undetermined parameters K and beta in the prediction model formula (1) by regression fitting₁And beta₂The numerical value of (A):

wherein K represents a site coefficient, beta₁Expressed is the attenuation coefficient, beta₂The attenuation coefficient caused by the change of the soil-water medium is shown;

step 3, according to the characteristics of the protected object of the bank building/structure in the near field of underwater blasting, consulting the maximum value V of V in the protection category in blasting safety regulations (GB6722-2014)_maxWill V_maxSubstituting into formula (2) to obtain maximum blasting explosive quantity Q_maxAnd guiding the engineering blasting:

step one, because the maximum allowable blasting vibration speed V is related to the blasting vibration dominant frequency, Q obtained in the step 3 is used_maxSubstituting the following formula (3) to obtain f_max:

Step two, since f is known_maxFrom the table lookup, V in the protection class can be known_max', will V_max' substitution into formula (2) to obtain Q_max’。

Comparison Q_maxAnd Q_max' selecting the minimum value as the maximum blasting explosive quantity of the blasting operation.

The maximum blasting explosive quantity is obtained in the steps 1-3, and the more accurate maximum blasting explosive quantity can be further calculated through further calculation in the steps one and two.

By referring to the safety regulations for blasting (GB6722-2014), as shown in table 2, if the dominant frequency value f and the allowable blasting vibration velocity V correspond to the same frequency range in table 2, it is determined that the allowable blasting vibration velocity V is reasonably selected.

TABLE 2 blasting vibration safety allowance criteria

The prediction model formula (1) is a functional relation of underwater blasting vibration strength derived by screening physical quantities related to blasting vibration dominant frequency and applying dimensional analysis in consideration of water and soil medium change characteristics accompanied by propagation of underwater blasting vibration on the bank and the adjacent areas thereof, and the specific derivation process is as follows:

when the propagation intensity attenuation rule of underwater blasting vibration on the bank and the neighboring area is researched, the influence of the change characteristics of the water and soil media on the propagation intensity of the blasting vibration is considered because the vibration propagation media relate to water and soil. The invention adopts a dimension analysis method to deduce a vibration intensity attenuation formula. The independent variables affecting the blasting vibration strength are given as follows: explosive quantity Q, explosive center distance R and distance R from measuring point to embankment_DThe density ρ of the medium, the propagation speed c of the shock wave in the medium and the detonation time t. The blasting vibration displacement mu, the vibration speed V, the vibration acceleration a and the vibration dominant frequency f are used as dependent variables, and the following functional expressions can be obtained according to Buchkinson's theorem (Π theorem) along with the variables affected by the dependent variables:

V＝f(Q,R,c,μ,a,f,ρ,R_D,t) (5)

from the formula (5), the total number of physical quantities for the analysis problem is n ═ 10. According to pi theorem, the basic variables are Q, R and c, so the dimension number of the basic variables is 3, the dimension numbers of the derived quantity and the dependent variable are 7, pi and pi₁～Π₇Representing dimensionless variables, then:

the substitution of formula (6) for formula (5) has:

if a dimensionless variable that does not contribute to the dependent variable Π can be noted as constant 1, then:

for the attributes of the blasting object field medium, both the medium density ρ and the propagation velocity c of the vibration wave can be recorded as constants, which can be further simplified by the formula (8):

considering blasting vibration and explosive quantity Q, explosive source distance R and distance R between measuring point and embankment_DCan be written as:

order toThen there are:

lnV₀＝α₁+β′₁lnQ-3β′₁lnR (11)

wherein alpha is₁+β′₁lnQ shows the combined effect of the medium condition and explosive quantity of the blasting object field on the velocity of the vibration particle; -beta'₁lnR represents the attenuation relation between the vibration particle speed and the explosion source distance R; beta'₁The damping coefficient is shown and mainly reflects the influence of the field medium condition on the velocity of the vibration particles.

Let lnK₁＝α₁，β₁＝-3β′₁Then, there are:

wherein the content of the first and second substances,the formula of Sadawski under open blasting conditions is shown.

A substitution of formula (12) for formula (10) has:

let lnK₂＝-ε₁，β₂＝-η₁Under the condition of considering the water area factor, the calculation formula of the velocity of the blasting vibration particle by introducing the distance ratio coefficient is as follows:

wherein, K₁The field coefficient of the blasting object is shown; k₂The landform influence coefficient is represented; beta is a₁The attenuation coefficient is shown; beta is a₂Expressed versus distance coefficient;

in the formula (14), K may be replaced by a coefficient K₁,K₂The method comprises the following steps:

in step 2, according to the actually measured blasting vibration data, regression fitting is carried out on undetermined parameters in the prediction model formula by adopting a linear least square method, and the specific fitting process is as follows:

taking logarithm on two sides of equal sign of formula (2) comprises:

let y equal lnV, beta₀＝lnK，x₂＝ln(R_DR), then formula (16) can be:

y＝β₀+β₁x₁+β₂x₂ (17)

to determine the parameter beta to be determined₀,β₁,β₂The optimum estimated value of (a) is required to make the vibration measured data y_iFitting value to equation (17)Minimum sum of squared residuals, i.e. solve the objective function:

from the theorem of the extreme value of the multivariate function, f is a minimum requirement:

then, a non-homogeneous linear equation set solution β can be obtained₀，β₁，β₂. Then, the variable substitution can be used to obtainThen the specific underwater blasting explosive loading expression is as follows:

wherein, the calculation model formula in step 3 is:

in the formula (21), the blasting explosive quantity Q can be calculated according to the maximum allowable vibration speed value V of the prediction point, so as to guide engineering blasting.

The foregoing illustrates and describes the principles, general features, and advantages of the present invention. It should be understood by those skilled in the art that the above embodiments do not limit the present invention in any way, and all technical solutions obtained by using equivalent alternatives or equivalent variations fall within the scope of the present invention.