CN110736400A - underwater drilling blasting vibration velocity calculation method considering rock internal structure - Google Patents

underwater drilling blasting vibration velocity calculation method considering rock internal structure Download PDF

Info

Publication number
CN110736400A
CN110736400A CN201911008520.6A CN201911008520A CN110736400A CN 110736400 A CN110736400 A CN 110736400A CN 201911008520 A CN201911008520 A CN 201911008520A CN 110736400 A CN110736400 A CN 110736400A
Authority
CN
China
Prior art keywords
formula
blasting
vibration
calculating
vibration speed
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911008520.6A
Other languages
Chinese (zh)
Other versions
CN110736400B (en
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.)
CHANGJIANG CHONGQING WATERWAY ENGINEERING BUREAU
Original Assignee
CHANGJIANG CHONGQING WATERWAY ENGINEERING BUREAU
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 CHANGJIANG CHONGQING WATERWAY ENGINEERING BUREAU filed Critical CHANGJIANG CHONGQING WATERWAY ENGINEERING BUREAU
Priority to CN201911008520.6A priority Critical patent/CN110736400B/en
Publication of CN110736400A publication Critical patent/CN110736400A/en
Application granted granted Critical
Publication of CN110736400B publication Critical patent/CN110736400B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F42—AMMUNITION; BLASTING
    • F42D—BLASTING
    • F42D3/00—Particular applications of blasting techniques
    • F42D3/04—Particular applications of blasting techniques for rock blasting
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B7/00—Special methods or apparatus for drilling
    • E21B7/007—Drilling by use of explosives
    • 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12—Simultaneous equations, e.g. systems of linear equations
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Algebra (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Software Systems (AREA)
  • Computing Systems (AREA)
  • Geochemistry & Mineralogy (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Operations Research (AREA)
  • Earth Drilling (AREA)

Abstract

In the embodiment of the application, firstly, blasting parameters for performing underwater blasting operation on a rock structure are obtained, wherein the blasting parameters comprise single-section maximum explosive quantity, a blasting center distance and an elevation difference, secondly, a -th peak value vibration speed is obtained according to the blasting parameters and a preset -th formula, and a second formula is established according to the -th peak value vibration speed, wherein the second formula meets normal distribution.

Description

underwater drilling blasting vibration velocity calculation method considering rock internal structure
Technical Field
The application relates to the technical field of underwater drilling blasting vibration velocity calculation, in particular to underwater drilling blasting vibration velocity calculation methods and devices considering internal rock structures, electronic equipment and computer readable storage media.
Background
With the vigorous development of shipping industry in China, -cycle development of underwater drilling blasting is achieved, and the peak vibration speed of the underwater drilling blasting is the most key evaluation index in response characteristic and dynamic stability evaluation of adjacent buildings (structures) under the action of blasting seismic waves.
However, the inventor researches and discovers that in the prior art, the research on the blasting vibration speed still stays at the stage of considering controllable variables such as the distance between the centers of explosion, the single-section maximum dosage, the height difference and the like, the influence of the discontinuity and complexity of the rock mass structure on the blasting vibration speed is not considered, and the accuracy of the calculation of the blasting vibration speed is not high.
Disclosure of Invention
In view of the above, the present application aims to provide methods and apparatuses for calculating vibration velocity of underwater drilling blasting considering internal rock structure, an electronic device and a computer readable storage medium, so as to solve the problems in the prior art.
In order to achieve the above purpose, the embodiment of the present application adopts the following technical solutions:
A method for calculating the vibration velocity of underwater drilling blasting by considering the internal structure of rock comprises the following steps:
acquiring blasting parameters for underwater blasting operation on a rock mass structure, wherein the blasting parameters comprise single-section maximum explosive quantity, an explosive center distance and an elevation difference;
calculating to obtain a peak value vibration speed according to the blasting parameters and a preset th formula, and establishing a second formula according to the peak value vibration speed, wherein the second formula meets normal distribution;
and calculating to obtain the blasting vibration speed for underwater blasting operation based on the second formula and a preset second peak value vibration speed.
In a preferred option of the embodiment of the present application, the step of establishing the second formula according to the th peak vibration velocity includes:
taking the peak vibration speed as the mathematical expectation of a second formula to be established, and calculating the variance of the second formula according to the peak vibration speed and preset vibration speed monitoring data;
the second formula is established based on the mathematical expectation and the variance.
In a preferred selection of the embodiment of the present application, the second formula includes:
Figure BDA0002243480020000021
wherein, F (V)S) Indicating the blasting vibration speed, V, at which the underwater blasting operation is carried outSRepresenting a preset second peak vibration speedDegree, VuRepresenting a threshold value of vibration speed, σ2Represents variance and μ represents mathematical expectation.
In a preferred selection of the embodiment of the present application, the -th formula includes:
Figure BDA0002243480020000022
wherein, VMRepresenting peak vibration speed, Q representing single-section maximum drug quantity, R representing a burst center distance, H representing an elevation difference, K and α representing parameters related to a propagation medium, and β representing an elevation difference factor.
In a preferred option of the embodiment of the present application, the data management method further includes:
and converting the second formula into a third formula meeting standard normal distribution, and calculating to obtain the single-section maximum dosage threshold according to the third formula and a preset probability.
In a preferred option of the embodiment of the present application, the step of converting the second formula into a third formula satisfying a standard normal distribution includes:
and performing linear transformation processing on the second formula to obtain a standard normal distribution formula, and obtaining the third formula according to the standard normal distribution formula.
In a preferred selection of the embodiment of the present application, the third formula includes:
Figure BDA0002243480020000031
wherein Q ispRepresenting a single segment maximum dose threshold, VuRepresenting a vibration speed threshold value, b representing a third peak vibration speed at a preset probability, mu representing a mathematical expectation, R representing a burst length, and H representing a height difference.
The embodiment of the application also provides underwater drilling blasting vibration velocity calculation devices considering rock internal structure, which comprises:
the device comprises a parameter acquisition module, a parameter calculation module and a parameter calculation module, wherein the parameter acquisition module is used for acquiring blasting parameters for carrying out underwater blasting operation on a rock mass structure, and the blasting parameters comprise single-section maximum explosive quantity, an explosive center distance and an elevation difference;
an peak vibration velocity calculating module, configured to calculate a peak vibration velocity according to the blasting parameters and a preset formula, and establish a second formula according to the peak vibration velocity, where the second formula satisfies normal distribution;
and the blasting vibration speed calculation module is used for calculating the blasting vibration speed for carrying out underwater blasting operation based on the second formula and a preset second peak value vibration speed.
The embodiment of the application also provides electronic devices, which comprise a memory and a processor, wherein the processor is used for executing an executable computer program stored in the memory to realize the underwater drilling blasting vibration velocity calculation method considering the rock internal structure.
The embodiment of the present application further provides computer-readable storage media, on which a computer program is stored, which when executed, implements the steps of the above method for calculating the vibration velocity of underwater drilling blasting considering the internal rock structure.
According to the underwater drilling blasting vibration velocity calculation method and device, the electronic equipment and the computer readable storage medium considering the rock internal structure, after the second formula meeting normal distribution is established according to the peak value vibration velocity, the change of the actual underwater drilling blasting vibration velocity is simulated by simulating the underwater drilling blasting vibration velocity fluctuating near the peak value vibration velocity through the second formula, so that the problem that the accuracy of blasting vibration velocity calculation is not high due to the fact that the underwater drilling blasting vibration velocity fluctuates near the peak value vibration velocity due to discontinuity and complexity of the rock structure in the prior art is solved, and the accuracy of blasting vibration velocity calculation is further improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained from the drawings without inventive effort.
Fig. 1 is a block diagram of an electronic device according to an embodiment of the present disclosure.
Fig. 2 is a schematic flow chart of methods for calculating vibration velocity of underwater drilling blasting considering internal rock structures according to embodiments of the present application.
Fig. 3 is a schematic diagram of a distance to a center of burst and a height difference provided in an embodiment of the present application.
Fig. 4 is a schematic flowchart of step S120 according to an embodiment of the present application.
Fig. 5 is a functional diagram of a standard normal distribution provided in the embodiment of the present application.
Fig. 6 is a structural block diagram of an underwater drilling blasting vibration velocity calculation device considering the internal rock structure according to an embodiment of the present application.
The figure comprises 10 electronic equipment, 12 memory, 14 processor, 100 underwater drilling blasting vibration velocity calculation device considering rock internal structure, 110 parameter acquisition module, 120- peak vibration velocity calculation module and 130 blasting vibration velocity calculation module.
Detailed Description
To further clarify the objects, aspects and advantages of the embodiments of the present application, reference will now be made in detail to the present embodiments of the application illustrated in the accompanying drawings, which form a part hereof, and to show by way of illustration, and not by way of limitation, some embodiments of the application .
Thus, the following detailed description of the embodiments of the present application, presented in the accompanying drawings, is not intended to limit the scope of the claimed application, but is merely representative of selected embodiments of the application. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
It should be noted that like reference numerals and letters refer to like items in the following figures, and thus once a item is defined in figures, it need not be further defined and explained by in subsequent figures.
As shown in fig. 1, the present embodiment provides electronic devices 10, wherein the electronic device 10 may include a memory 12, a processor 14, and a device 100 for calculating the vibration velocity of underwater drilling blasting considering the internal structure of rock.
The underwater drilling burst vibration velocity calculating device 100 considering the rock internal structure comprises at least software functional modules which can be stored in the memory 12 in the form of software or firmware (firmware), and the processor 14 is used for executing executable computer programs stored in the memory 12, for example, the software functional modules and the computer programs and the like which are included in the underwater drilling burst vibration velocity calculating device 100 considering the rock internal structure, so as to realize the underwater drilling burst vibration velocity calculating method considering the rock internal structure.
The Memory 12 may be, but is not limited to, a Random Access Memory (RAM), a Read Only Memory (ROM), a Programmable Read-Only Memory (PROM), an Erasable Read-Only Memory (EPROM), an electrically Erasable Read-Only Memory (EEPROM), and the like.
The Processor 14 may be an type integrated circuit Chip having signal Processing capability, and the Processor 14 may be a general-purpose Processor including a Central Processing Unit (CPU), a Network Processor (NP), a System on Chip (SoC), and the like.
It will be appreciated that the configuration shown in FIG. 1 is merely illustrative and that the electronic device 10 may include more or fewer components than shown in FIG. 1 or may have a different configuration than shown in FIG. 1.
With reference to fig. 2, the embodiment of the present application further provides methods for calculating the vibration velocity of underwater drilling blasting considering the internal rock structure, which are applicable to the above electronic device 10, wherein the method steps defined by the flow related to the method for calculating the vibration velocity of underwater drilling blasting considering the internal rock structure can be implemented by the electronic device 10, and the specific flow shown in fig. 2 will be described in detail below.
And step S110, acquiring blasting parameters for underwater blasting operation on the rock mass structure.
In the embodiment of the present application, the electronic device 10 may obtain the blasting parameter through a detection device connected to the electronic device 10. Wherein the blasting parameters comprise single-section maximum dosage, a blasting center distance and a height difference.
And S120, calculating to obtain a peak value vibration speed according to the blasting parameters and a preset th formula, and establishing a second formula according to the peak value vibration speed.
In the embodiment of the present application, after the blasting parameters are obtained in step S110, the blasting parameters may be substituted into a preset th formula to calculate a th peak vibration speed, and a second formula is established according to the th peak vibration speed.
And step S130, calculating to obtain the blasting vibration speed for underwater blasting operation based on the second formula and a preset second peak value vibration speed.
In this embodiment of the application, after the second formula is established in step S120, the preset second peak vibration speed may be substituted into the second formula to calculate the blasting vibration speed.
Wherein the electronic device 10 may obtain the second peak vibration speed through a detection device connected to the electronic device 10. That is, the second peak vibration speed refers to the detected underwater drilling blasting vibration speed in the historical data, and the current blasting vibration speed can be calculated according to the blasting vibration speed in the historical data.
By the method, after the second formula meeting the normal distribution is established according to the th peak value vibration speed, the second formula simulates the change of the actual underwater drilling blasting vibration speed by simulating the underwater drilling blasting vibration speed fluctuating near the th peak value vibration speed, so that the problem that the accuracy of calculating the blasting vibration speed is not high due to the fluctuation of the underwater drilling blasting vibration speed near the th peak value vibration speed caused by the discontinuity and complexity of the rock structure in the prior art is solved, and the accuracy of calculating the blasting vibration speed is improved.
In step S110, the underwater drilling blasting is a method of drilling a hole in a damaged object such as an underwater rock or a structure, and blasting the hole with an explosive charge. The method can fully utilize the explosive force and the brisance of the explosive to break the rock, so the blasting effect is good, the unit explosive consumption is small, the generated impact is small, and the method becomes a main blasting method of the modern underwater engineering.
In detail, the single-section maximum explosive amount refers to explosive amount, and referring to fig. 3, point a is a blasting point for placing explosives, point B is a detection point for placing detection equipment, the distance R between the blasting center and the point a is the distance between the point B and the blasting center, and the height difference H refers to the difference between the heights of A, B points, that is, the difference between two leveling planes passing through the two points.
For step S120, it should be noted that the specific formula of the preset th formula is not limited, and may be set according to the actual application requirement.
For example, in alternative examples, the formula may include:
Figure BDA0002243480020000071
wherein, VMRepresenting peak vibration speed, Q representing single-section maximum drug quantity, R representing a burst center distance, H representing an elevation difference, K and α representing parameters related to a propagation medium, and β representing an elevation difference factor.
, distance proportional medicine quantity p can be introduced for determining related parameters K, α and βr＝Q1/3Dosage p with/R and high proportionk＝Q1/3H, the following equation is obtained:
Vm＝Kpr αpk β； (1-2)
by performing the linearization process on the above equation (1-2), the following equation can be obtained:
lnV＝αlnpr+βlnpk+lnK； (1-3)
y may be set to lnV, a0＝lnK，a1＝α，a2＝β，x1＝lnpr，x2＝lnpkTo further simplify the above equation (1-3), the following equation is obtained:
y＝a0+a1x1+a2x2； (1-4)
wherein, the formula (1-4) is a binary linear regression equation requiring multiple linear regression analysis, and the formula (1-4) can be subjected to regression analysis according to least square estimation principleMAnd measured value VS(blasting vibration monitoring speed in underwater drilling blasting site) is minimum, and group a is corresponded0 *，a1 *And a2 *The value of (c).
VMAnd VSQ for total deviation therebetweeneAs a representation, the equations (1-5) can be referred to, where i represents the number of sets of field explosion monitoring data. According to the multivariate function extremum solving method, the equation set (1-6) can be established, and a can be obtained by solving the equation set (1-6)0 *，a1 *And a2 *The value of (c).
Figure BDA0002243480020000081
Figure BDA0002243480020000082
The solution can be found in that α ═ a1 *，β＝a2 *，Therefore, the th formula can be represented by formulas (1-7).
Figure BDA0002243480020000083
As can be seen from equations (1-7), VMAnd Q1/3R and Q1/3the/H is a functional relationship. However, the rock soil is a geologic body formed over a long geological age, and the interior of the geologic body contains a large number of cracks, the cracks are quite complex, discontinuity and complexity are caused in the interior of the rock, and therefore the actual peak value V of particle vibration caused by underwater drilling explosionSOften the randomness is relatively large. VSAnd Q1/3R and Q1/3The relationship of/H is not a mathematical "functional relationship", but a correlation of random variables. Here VSIs a random variable, Q1/3R and Q1/3the/H is a controllable variable, also called an independent variable. Albeit VSAnd Q1/3R and Q1/3The relationship of/H is not mathematically a "functional relationship", but studies have shown that Q1/3R and Q1/3the/H is a key factor influencing the vibration speed of the underwater drilling blasting, so that the formula can reflect the predicted underwater drilling blasting vibration speed value of the real blasting vibration speed to the extent of .
Optionally, the specific steps of establishing the second formula according to the th peak vibration speed are not limited, and may be set according to actual application requirements.
For example, in alternative examples, in conjunction with fig. 4, step S120 may include step S121 and step S122.
And step S121, taking the th peak vibration speed as the mathematical expectation of a second formula to be established, and calculating the variance of the second formula according to the th peak vibration speed and preset vibration speed monitoring data.
Step S122, establishing the second formula according to the mathematical expectation and the variance.
In step S121, it should be noted that, since it is mathematically assumed that the random variable approximately follows a normal distribution, V can be assumedSAlso approximately follows a normal distribution with a mean value of VMVariance is VSAnd VMTotal deviation of (a).
Due to VSApproximately obey a normal distribution, which can be written as VS～N(μ,σ2). Where mu and sigma2The value of (C) can be calculated according to the equations (2-1) and (2-2).
Figure BDA0002243480020000091
Figure BDA0002243480020000092
And n is the group number of the underwater drilling blasting monitoring data.
The probability density and distribution function of a random variable following a normal distribution can be obtained from the definition of the random variable following a normal distribution, see equations (2-3) and (2-4). according to the definition of the random variable, let X be random variables, X be any real number, and the function f (X) ═ P { X ≦ X } is referred to as the distribution function of the random variable X, where ∞<x<∞。VSApproximately obey a normal distribution, then F (V)S)＝P{VS≤VuIs a random variable VSThe distribution function of (2).
Figure BDA0002243480020000101
Figure BDA0002243480020000102
For step S122, it should be noted that the second formula may include formula (2-4), where F (V)S) Is shown inBlasting vibration velocity, V, of underwater blasting operationsSRepresenting a preset second peak vibration speed, VuRepresenting a threshold value of vibration speed, σ2Represents variance and μ represents mathematical expectation.
In detail, VSAnd the second peak value vibration speed is preset and can be obtained through preset vibration speed monitoring data. VuThe peak vibration speed which can be borne by an underwater drilling and blasting adjacent building (structure) before the building is damaged after bearing blasting load can be represented, and the values are shown in table 1.
TABLE 1 vibration speed threshold
Figure BDA0002243480020000103
Figure BDA0002243480020000111
, in this embodiment, after obtaining the second formula, the probability that the adjacent building can be protected by blasting at a certain time can be found based on the probability (the probability that the adjacent building can be protected by blasting) that the corresponding maximum explosive quantity Q of blasting single segment is PpThe second formula is an -like normal distribution underwater drilling blasting vibration velocity calculation formula considering the internal structure of the rock, and the formula needs to be converted into a standard normal distribution formula to solve the problem that the maximum explosive quantity of a blasting single section is Q when the probability that the blasting adjacent building structure can be protected at a certain time is Pp。
That is, the underwater drilling blasting vibration velocity calculation method considering the internal rock structure may further include the sub-steps of:
and converting the second formula into a third formula meeting standard normal distribution, and calculating to obtain the single-section maximum dosage threshold according to the third formula and a preset probability.
It should be noted that, the specific step of converting the second formula into the third formula satisfying the standard normal distribution is not limited, and may be set according to the actual application requirement.
For example, in alternative examples, the following sub-steps may be included:
and performing linear transformation processing on the second formula to obtain a standard normal distribution formula, and obtaining the third formula according to the standard normal distribution formula.
In conjunction with fig. 5, when μ is 0 and σ is larger than the threshold value2When the random variable X is referred to as "1", it follows a standard normal distribution, and is denoted as "X to N (0, 1"), and the probability density and the distribution function thereof can be expressed by equations (2-5) and (2-6), respectively.
Figure BDA0002243480020000112
Figure BDA0002243480020000121
It is known that Φ (-x) — 1 Φ (x), the size of Φ (x) can be directly obtained from table 2, and -like normal distribution can be converted into standard normal distribution by linear transformation.
TABLE 2 Standard Normal distribution probability values
If X to N (mu, sigma)2) If Z is (X- μ)/σ — N (0,1), the standard normal distribution formula may be:
Figure BDA0002243480020000123
where b ═ V may be provideduMu)/sigma, when phi (b) is P, a third peak vibration speed b corresponding to the preset probability P can be obtained through table 2, and then a corresponding single-section maximum medicine quantity threshold value Q is obtained through calculation of a formula (2-7)p。
Figure BDA0002243480020000131
Wherein the third formula may include formula (2-7), QpRepresenting a single segment maximum dose threshold, VuRepresenting a vibration speed threshold value, b representing a third peak vibration speed at a preset probability, mu representing a mathematical expectation, R representing a burst length, and H representing a height difference.
, in this embodiment, after obtaining the second formula, the probability P of damage to the adjacent building structure can be obtained by the single-stage maximum drug quantity, that is, the peak vibration velocity V can be obtained by the formula and the single-stage maximum drug quantityMFurther, the variance σ can be obtained2And a mathematical expectation μ, the corresponding blasting vibration velocity can be obtained by the second formula, so that the corresponding probability P is obtained by table 2.
With reference to fig. 6, underwater drilling blasting vibration velocity calculation apparatuses 100 considering rock internal structures are further provided, and may be applied to the electronic device 10 described above, where the underwater drilling blasting vibration velocity calculation apparatus 100 considering rock internal structures may include a parameter acquisition module 110, a peak vibration velocity calculation module 120, and a blasting vibration velocity calculation module 130.
The parameter obtaining module 110 is configured to obtain blasting parameters for performing underwater blasting operation on a rock mass structure, where the blasting parameters include single-stage maximum explosive quantity, an explosive center distance, and a height difference. In this embodiment, the parameter obtaining module 110 may be configured to execute step S110 shown in fig. 2, and for the relevant content of the parameter obtaining module 110, reference may be made to the foregoing detailed description of step S110.
The peak vibration velocity calculating module 120 is configured to calculate a peak vibration velocity according to the blasting parameters and a preset formula, and establish a second formula according to the peak vibration velocity, where the second formula satisfies a normal distribution, in this embodiment, the peak vibration velocity calculating module 120 may be configured to execute step S120 shown in fig. 2, and for the peak vibration velocity calculating module 120, the detailed description of step S120 may be referred to.
The blasting vibration velocity calculation module 130 is configured to calculate a blasting vibration velocity for performing underwater blasting operation based on the second formula and a preset second peak vibration velocity. In this embodiment, the blasting vibration speed calculation module 130 may be configured to execute step S130 shown in fig. 2, and for the relevant content of the blasting vibration speed calculation module 130, reference may be made to the foregoing detailed description of step S130.
In the embodiment of the present application, there are also computer-readable storage media, which are provided corresponding to the above-mentioned underwater drilling blasting vibration velocity calculation method considering the rock internal structure, and a computer program is stored in the computer-readable storage media, and the computer program is executed to execute the steps of the above-mentioned underwater drilling blasting vibration velocity calculation method considering the rock internal structure.
The steps executed when the computer program runs are not described in again, and reference may be made to the above explanation of the underwater drilling blasting vibration velocity calculation method considering the internal rock structure.
To sum up, after the underwater drilling blasting vibration velocity calculation methods and apparatuses, the electronic device, and the computer-readable storage medium considering the rock internal structure provided by the embodiments of the present application are established, according to the peak vibration velocity, the second formula satisfying the normal distribution is established, and the second formula simulates the actual underwater drilling blasting vibration velocity variation by simulating the underwater drilling blasting vibration velocity fluctuating near the peak vibration velocity, so as to solve the problem in the prior art that the accuracy of the calculation of the blasting vibration velocity is not high due to the fluctuation of the underwater drilling blasting vibration velocity near the peak vibration velocity due to the discontinuity and complexity of the rock structure, and further improve the accuracy of the calculation of the blasting vibration velocity.
The above description is only a preferred embodiment of the present application and is not intended to limit the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims (10)

1, underwater drilling blasting vibration velocity calculation method considering rock internal structure, which is characterized by comprising:
acquiring blasting parameters for underwater blasting operation on a rock mass structure, wherein the blasting parameters comprise single-section maximum explosive quantity, an explosive center distance and an elevation difference;
calculating to obtain a peak value vibration speed according to the blasting parameters and a preset th formula, and establishing a second formula according to the peak value vibration speed, wherein the second formula meets normal distribution;
and calculating to obtain the blasting vibration speed for underwater blasting operation based on the second formula and a preset second peak value vibration speed.
2. The method for calculating the vibration velocity of underwater drilling blasting considering the internal structure of rock as claimed in claim 1, wherein said step of establishing the second formula according to the th peak vibration velocity comprises:
taking the peak vibration speed as the mathematical expectation of a second formula to be established, and calculating the variance of the second formula according to the peak vibration speed and preset vibration speed monitoring data;
the second formula is established based on the mathematical expectation and the variance.
3. The underwater drilling blasting vibration velocity calculation method considering the internal rock structure of claim 2, wherein the second formula includes:
Figure FDA0002243480010000011
wherein, F (V)S) Indicating the blasting vibration speed, V, at which the underwater blasting operation is carried outSRepresenting a preset second peak vibration speed, VuRepresenting a threshold value of vibration speed, σ2Represents variance and μ represents mathematical expectation.
4. The method of calculating a vibro velocity of underwater drilling blasting considering internal rock formations according to any of of claim 1-3, wherein the formula comprises:
Figure FDA0002243480010000021
wherein, VMRepresenting peak vibration speed, Q representing single-section maximum drug quantity, R representing a burst center distance, H representing an elevation difference, K and α representing parameters related to a propagation medium, and β representing an elevation difference factor.
5. The method for calculating the vibration velocity of underwater drilling blasting considering internal rock structures according to any of claim 1 to 3, further comprising:
and converting the second formula into a third formula meeting standard normal distribution, and calculating to obtain the single-section maximum dosage threshold according to the third formula and a preset probability.
6. The method for calculating underwater drilling blasting vibration velocity in consideration of internal rock formation according to claim 5, wherein the step of converting the second formula into a third formula satisfying a normal distribution comprises:
and performing linear transformation processing on the second formula to obtain a standard normal distribution formula, and obtaining the third formula according to the standard normal distribution formula.
7. The underwater drilling blasting vibration velocity calculation method considering the internal rock structure of claim 5, wherein the third formula comprises:
Figure FDA0002243480010000022
wherein Q ispRepresenting a single segment maximum dose threshold, VuRepresenting a vibration speed threshold, b representing the second under a preset probabilityThree peak vibration velocities, μ represents the mathematical expectation, R represents the centre of burst distance, and H represents the elevation difference.
8, A device for calculating vibration velocity of underwater drilling blasting considering internal rock structure, comprising:
the device comprises a parameter acquisition module, a parameter calculation module and a parameter calculation module, wherein the parameter acquisition module is used for acquiring blasting parameters for carrying out underwater blasting operation on a rock mass structure, and the blasting parameters comprise single-section maximum explosive quantity, an explosive center distance and an elevation difference;
an peak vibration velocity calculating module, configured to calculate a peak vibration velocity according to the blasting parameters and a preset formula, and establish a second formula according to the peak vibration velocity, where the second formula satisfies normal distribution;
and the blasting vibration speed calculation module is used for calculating the blasting vibration speed for carrying out underwater blasting operation based on the second formula and a preset second peak value vibration speed.
An electronic device of , comprising a memory and a processor, wherein the processor is configured to execute an executable computer program stored in the memory to implement the method of calculating the vibration velocity of underwater drilling blasting considering internal rock structure of any of claims 1-7.
10, computer-readable storage medium, characterized in that, a computer program is stored thereon, which when executed, implements the steps of the method for calculating the vibration velocity of underwater drilling blasting considering the internal rock formation according to any of claims 1-7.
CN201911008520.6A 2019-10-22 2019-10-22 Underwater drilling blasting vibration velocity calculation method considering internal rock structure Active CN110736400B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911008520.6A CN110736400B (en) 2019-10-22 2019-10-22 Underwater drilling blasting vibration velocity calculation method considering internal rock structure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911008520.6A CN110736400B (en) 2019-10-22 2019-10-22 Underwater drilling blasting vibration velocity calculation method considering internal rock structure

Publications (2)

Publication Number Publication Date
CN110736400A true CN110736400A (en) 2020-01-31
CN110736400B CN110736400B (en) 2022-09-13

Family

ID=69270893

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911008520.6A Active CN110736400B (en) 2019-10-22 2019-10-22 Underwater drilling blasting vibration velocity calculation method considering internal rock structure

Country Status (1)

Country Link
CN (1) CN110736400B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113607017A (en) * 2021-08-09 2021-11-05 核工业南京建设集团有限公司 Method for predicting underwater blasting explosive loading based on test blasting data
CN114781040A (en) * 2022-05-10 2022-07-22 青岛理工大学 Building blasting vibration safety evaluation method, device, equipment and medium
CN116242209A (en) * 2022-12-14 2023-06-09 青岛理工大学 Vibration damping rate calculation method, system and equipment for vibration damping hole and readable storage medium
CN116258285A (en) * 2023-05-16 2023-06-13 青岛理工大学 Porous small-clear-distance tunnel blasting vibration speed prediction method, device, equipment and medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2006108369A (en) * 2006-03-16 2007-09-27 Государственное образовательное учреждение высшего профессионального образовани "Тихоокеанский государственный университет" (RU) METHOD FOR FORMING CHARGES WITH AIR SPACES
CN101608885A (en) * 2009-07-08 2009-12-23 长江重庆航道工程局 Underwater drilling and blasting construction method of mountain river
CN202141517U (en) * 2011-03-11 2012-02-08 中国人民解放军理工大学工程兵工程学院 Blasting vibration forecasting system
CN106289371A (en) * 2016-07-22 2017-01-04 长江重庆航道工程局 The safe determinator of underwater demolition and method
CN109506532A (en) * 2018-10-25 2019-03-22 中国地质大学（武汉） The determination method of underwater drilling and blasting explosive specific charge

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2006108369A (en) * 2006-03-16 2007-09-27 Государственное образовательное учреждение высшего профессионального образовани "Тихоокеанский государственный университет" (RU) METHOD FOR FORMING CHARGES WITH AIR SPACES
CN101608885A (en) * 2009-07-08 2009-12-23 长江重庆航道工程局 Underwater drilling and blasting construction method of mountain river
CN202141517U (en) * 2011-03-11 2012-02-08 中国人民解放军理工大学工程兵工程学院 Blasting vibration forecasting system
CN106289371A (en) * 2016-07-22 2017-01-04 长江重庆航道工程局 The safe determinator of underwater demolition and method
CN109506532A (en) * 2018-10-25 2019-03-22 中国地质大学（武汉） The determination method of underwater drilling and blasting explosive specific charge

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《中国水利发电工程》编审委员会: "《中国水力发电工程 施工卷》", 31 August 2000, 中国电力出版社 *
张小军等: "基于正态分布的爆破振动评价与安全药量计算", 《爆炸与冲击》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113607017A (en) * 2021-08-09 2021-11-05 核工业南京建设集团有限公司 Method for predicting underwater blasting explosive loading based on test blasting data
CN114781040A (en) * 2022-05-10 2022-07-22 青岛理工大学 Building blasting vibration safety evaluation method, device, equipment and medium
CN114781040B (en) * 2022-05-10 2023-07-07 青岛理工大学 Building blasting vibration safety evaluation method, device, equipment and medium
CN116242209A (en) * 2022-12-14 2023-06-09 青岛理工大学 Vibration damping rate calculation method, system and equipment for vibration damping hole and readable storage medium
CN116242209B (en) * 2022-12-14 2023-08-18 青岛理工大学 Vibration damping rate calculation method, system and equipment for vibration damping hole and readable storage medium
CN116258285A (en) * 2023-05-16 2023-06-13 青岛理工大学 Porous small-clear-distance tunnel blasting vibration speed prediction method, device, equipment and medium
CN116258285B (en) * 2023-05-16 2023-08-25 青岛理工大学 Porous small-clear-distance tunnel blasting vibration speed prediction method, device, equipment and medium

Also Published As

Publication number Publication date
CN110736400B (en) 2022-09-13

Similar Documents

Publication Publication Date Title
CN110736400A (en) underwater drilling blasting vibration velocity calculation method considering rock internal structure
US20200089733A1 (en) Sensor placement method for reducing uncertainty of structural modal identification
CN107092759B (en) Dam displacement monitoring point optimal arrangement method based on gravity dam foundation parameter inversion
CN108509754B (en) High-density packaged integrated circuit bonding wire touch risk assessment method
CN112946783A (en) Method, device and equipment for determining hydrate saturation
CN106768581B (en) A kind of horizontal stress coefficient of variation measurement method and device
Feng et al. Climate network stability measures of El Niño variability
Li et al. Modal identification of civil structures via covariance-driven stochastic subspace method
Huang et al. Wavelet‐based approach of time series model for modal identification of a bridge with incomplete input
CN108228959A (en) Using the method for Random censorship estimating system virtual condition and using its wave filter
CN111080020B (en) Robustness evaluation method and device for drilling arrangement scheme
Choi et al. Efficient calculation of statistical moments for structural health monitoring
TWI636276B (en) Method of determining earthquake with artificial intelligence and earthquake detecting system
CN108256261B (en) Stability evaluation method and device for offshore floating structure
Habibullah et al. On a newly developed estimator for more accurate modeling with an application to civil engineering
Saghir et al. On the Performance of Phase‐I Bivariate Dispersion Charts to Non‐Normality
CN111723858A (en) Material pile quality estimation method and device based on SVM, electronic equipment and medium
CN113901552A (en) Method and device for establishing weir body model
JP2022148889A (en) Damage degree estimation method of building, damage degree estimation system of building, damage degree learning method of building, damage degree learning system of building, and program
Pietruszczak et al. Analysis of Selected Dynamic Properties of Fractional Order Accelerometers for Application in Telematic Equipment
CN117235461B (en) Roadway support grade determining method and device, electronic equipment and storage medium
CN111366972B (en) Goaf power transmission line monitoring method and device
Sun et al. Extrapolation reconstruction of wind pressure fields on the claddings of high-rise buildings
Bu et al. Second order hierarchical partial least squares regression-polynomial chaos expansion for global sensitivity and reliability analyses of high-dimensional models
JP6778628B2 (en) Underground structure estimator

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