CN116242209B - Vibration damping rate calculation method, system and equipment for vibration damping hole and readable storage medium - Google Patents

Vibration damping rate calculation method, system and equipment for vibration damping hole and readable storage medium Download PDF

Info

Publication number
CN116242209B
CN116242209B CN202211601008.4A CN202211601008A CN116242209B CN 116242209 B CN116242209 B CN 116242209B CN 202211601008 A CN202211601008 A CN 202211601008A CN 116242209 B CN116242209 B CN 116242209B
Authority
CN
China
Prior art keywords
detection point
information
distance
explosion
blasting
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202211601008.4A
Other languages
Chinese (zh)
Other versions
CN116242209A (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.)
Linyi Urban Construction Group Co ltd
Qingdao No1 Municipal Engineering Co ltd
Shandong Expressway Bridge Group Co ltd
Qingdao University of Technology
First Construction Co Ltd of China Construction Eighth Engineering Division Co Ltd
Shandong Luqiao Group Co Ltd
Original Assignee
Linyi Urban Construction Group Co ltd
Qingdao No1 Municipal Engineering Co ltd
Shandong Expressway Bridge Group Co ltd
Qingdao University of Technology
First Construction Co Ltd of China Construction Eighth Engineering Division Co Ltd
Shandong Luqiao 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 Linyi Urban Construction Group Co ltd, Qingdao No1 Municipal Engineering Co ltd, Shandong Expressway Bridge Group Co ltd, Qingdao University of Technology, First Construction Co Ltd of China Construction Eighth Engineering Division Co Ltd, Shandong Luqiao Group Co Ltd filed Critical Linyi Urban Construction Group Co ltd
Priority to CN202211601008.4A priority Critical patent/CN116242209B/en
Publication of CN116242209A publication Critical patent/CN116242209A/en
Application granted granted Critical
Publication of CN116242209B publication Critical patent/CN116242209B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F42AMMUNITION; BLASTING
    • F42DBLASTING
    • F42D1/00Blasting methods or apparatus, e.g. loading or tamping
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F42AMMUNITION; BLASTING
    • F42DBLASTING
    • F42D3/00Particular applications of blasting techniques
    • F42D3/04Particular applications of blasting techniques for rock blasting
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Algebra (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a vibration reduction rate calculation method, a system, equipment and a readable storage medium of a vibration reduction hole, and relates to the technical field of vibration reduction hole calculation; fitting and calculating according to the first explosive distance, the second information and the third information to obtain a first fitting coefficient and a second fitting coefficient; calculating according to the second explosive distance, the second information, the first fitting coefficient and the second fitting coefficient to obtain fourth information, wherein the fourth information comprises the theoretical vibration speed of the explosion of the second detection point under the first explosion; and calculating according to the third information and the fourth information to obtain fifth information, wherein the fifth information is the first blasting vibration reduction rate corresponding to the first height. The method can directly determine the vibration reduction effect of the vibration reduction hole, thereby playing a role in predicting the explosion vibration speed after vibration reduction and playing a certain role in protecting underground pipelines, building structures and the like possibly existing in front of the face.

Description

Vibration damping rate calculation method, system and equipment for vibration damping hole and readable storage medium
Technical Field
The invention relates to the technical field of calculation of vibration damping holes, in particular to a vibration damping rate calculation method, a vibration damping rate calculation system, vibration damping rate calculation equipment and a readable storage medium of the vibration damping holes.
Background
In the existing engineering practice and related research, the vibration damping rate research of the blasting vibration damping hole actually comprises the combined damping action of surrounding rock and the vibration damping hole on a propagation path, and the existing method cannot know the vibration speed damping condition of the blasting vibration wave after passing through the vibration damping hole, so that vibration damage can be caused to a building structure, an underground pipeline and the like existing in front of a tunnel face. Therefore, it is needed to provide a precise calculation method of the vibration damping rate of the vibration damping hole to directly determine the vibration damping effect of the vibration damping hole, so as to play a role in predicting the explosion vibration speed after vibration damping, and also can be used for assisting in designing the explosion of the vibration damping hole, thereby playing a certain role in protecting underground pipelines, building structures and the like existing in front of the face.
Disclosure of Invention
The technical problems to be solved by the application are as follows: provided are a vibration damping rate calculation method, system, device and readable storage medium for a vibration damping hole, which solve the above problems.
In order to achieve the above object, the embodiment of the present application provides the following technical solutions:
in a first aspect, an embodiment of the present application provides a method for calculating a damping rate of a damping hole, where the method includes: acquiring first information, second information and third information, wherein the first information comprises a first core distance and a second core distance under the first blasting, the first core distance is a core distance from the core of the first blasting to a first detection point, the second core distance is a core distance from the core of the first blasting to a second detection point, the difference between the first distance of the first detection point and the second distance of the second detection point is a wind well diameter, the first detection point and the second detection point are both arranged on the wind well wall, the height of the first detection point and the height of the second detection point are both the first height, the first distance is a horizontal distance from the core of the first blasting to the first detection point, the second distance is a horizontal distance from the core of the first blasting to the second detection point, the second information comprises an initiating explosive amount, and the third information comprises a first blasting vibration speed acquired by the first detection point and a first blasting vibration speed acquired by the second detection point; according to the first explosive center distance, the second information and the first explosive vibration speed acquired by the first detection point, fitting and calculating to obtain a first fitting coefficient and a second fitting coefficient; calculating according to the second explosive distance, the second information, the first fitting coefficient and the second fitting coefficient to obtain fourth information, wherein the fourth information comprises the theoretical vibration speed of blasting of a second detection point under the first blasting; and calculating according to the first blasting vibration speed and fourth information acquired by the second detection point to obtain fifth information, wherein the fifth information is the first blasting vibration reduction rate corresponding to the first height.
In a second aspect, an embodiment of the present application provides a vibration damping rate calculation system for a vibration damping hole, the system including: the device comprises an acquisition module, a control module and a control module, wherein the acquisition module is used for acquiring first information, second information and third information, the first information comprises a first core distance and a second core distance under the first blasting, the first core distance is a core distance from the first core to the first detection point, the second core distance is a core distance from the first core to the second detection point, the difference between the first distance of the first detection point and the second distance of the second detection point is a wind shaft diameter, the first detection point and the second detection point are both arranged on the wind shaft wall, the height of the first detection point and the height of the second detection point are both the first height, the first distance is a horizontal distance from the core to the first detection point, the second distance is a horizontal distance from the core to the second detection point, the second information comprises a detonation quantity, and the third information comprises a first detonation velocity acquired by the first detection point and a first detonation velocity acquired by the first detection point; the first calculation module is used for obtaining a first fitting coefficient and a second fitting coefficient according to the first explosive distance, the second information and the first explosion vibration speed acquired by the first detection point; the second calculation module is used for calculating according to the second explosive distance, the second information, the first fitting coefficient and the second fitting coefficient to obtain fourth information, wherein the fourth information comprises the theoretical vibration speed of blasting of a second detection point under the first blasting; the first processing module is used for calculating according to the first blasting vibration speed and the fourth information acquired by the second detection point to obtain fifth information, wherein the fifth information is the first blasting vibration reduction rate corresponding to the first height.
In a third aspect, an embodiment of the present application provides a vibration damping rate calculation apparatus for a vibration damping orifice, the apparatus including a memory and a processor. The memory is used for storing a computer program; and the processor is used for realizing the step of the vibration damping rate calculating method of the vibration damping hole when executing the computer program.
In a fourth aspect, an embodiment of the present application provides a readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the damping ratio calculation method of a damping hole described above.
The beneficial effects of the application are as follows: according to the method, blasting vibration detection points are arranged in surrounding rocks at two sides of a wind well excavated above a tunnel, the wind well is regarded as a vibration damping hole, a Sarkowski formula in the surrounding rocks at two sides is obtained through fitting, and the vibration damping rate of the vibration damping hole is determined based on the attenuation of the vibration speeds of the detection points after passing through the vibration damping hole. In actual construction, wind wells with different diameters (regarded as vibration damping holes) can be excavated respectively to obtain vibration speed attenuation rates of the vibration damping holes with different diameters, and guidance is provided for blasting-like vibration damping engineering.
Drawings
FIG. 1 is a schematic flow chart of a method for calculating vibration damping rate of a vibration damping hole according to an embodiment of the present application;
FIG. 2 is a schematic view of a first height structure of the present invention;
FIG. 3 is a schematic view of the tangential distance structure from the center of the explosion to the axis of the wind well in the first explosion according to the present invention;
FIG. 4 is a schematic view of the structure of the present invention with the preset height set equal to the first height;
FIG. 5 is a schematic diagram of a vibration damping rate calculation device with vibration damping holes according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of a vibration damping rate calculation system of a vibration damping hole according to an embodiment of the present invention;
FIG. 7 is a schematic diagram of a second module according to the present invention;
the marks in the figure:
1. a core of the first explosion; 2. a secondary blasting core; 3. a formation; 4. ground surface; 5. a wind shaft; 6. a tunnel; 901. an acquisition module; 902. a first computing module; 903. a second computing module; 904. a first processing module; 905. a second processing module; 906. a third processing module; 907. a fourth processing module; 908. a fifth processing module; 909. a sixth processing module; 9011. a first acquisition unit; 9012. a first calculation unit; 9013. a second acquisition unit; 9014. a second calculation unit; 9015. a third acquisition unit; 9016. a third calculation unit; 9051. a fourth acquisition unit; 9052. a ninth calculation unit; 9053. a tenth calculation unit; 9054. an eleventh calculation unit; 9055. a twelfth calculation unit; 9060. a tenth acquisition unit; 9061. a thirteenth calculation unit; 9062. a fourteenth calculation unit; 9063. a fifteenth calculation unit; 9070. an eleventh acquisition unit; 9071. a sixteenth calculation unit; 9072. a seventeenth calculation unit; 9073. an eighteenth calculation unit; 9074. a nineteenth calculation unit; 9081. a twentieth calculation unit; 9082. a twenty-first calculation unit; 9083. a twenty-second calculation unit; 9091. a twelfth acquisition unit; 9092. a twenty-third calculation unit; 9093. a twenty-fourth calculation unit; 9094. a twenty-fifth calculation unit; 9095. a twenty-sixth calculation unit; 90511. a fifth acquisition unit; 90512. a fourth calculation unit; 90513. a sixth acquisition unit; 90514. a fifth calculation unit; 90515. a seventh acquisition unit; 90516. a sixth calculation unit; 90517. an eighth acquisition unit; 90518. a seventh calculation unit; 90519. a ninth acquisition unit; 90520. an eighth calculation unit; 800. damping rate calculation means for the damping hole; 801. a processor; 802. a memory; 803. a multimedia component; 804. an I/O interface; 805. a communication component.
Detailed Description
Example 1:
the embodiment provides a vibration damping rate calculating method of a vibration damping hole, referring to fig. 1, the method is shown to include step S1, step S2, step S3 and step S4.
Step S1, acquiring first information, second information and third information, wherein the first information comprises a first core distance and a second core distance under the first blasting, the first core distance is a core distance from the first blasting to the first detection point, the second core distance is a core distance from the first blasting to the second detection point, the difference between the first distance of the first detection point and the second distance of the second detection point is a wind shaft diameter, the first detection point and the second detection point are both arranged on the wind shaft wall, the height of the first detection point and the height of the second detection point are both the first height, the first distance is a horizontal distance from the core of the first blasting to the first detection point, the second distance is a horizontal distance from the core of the first blasting to the second detection point, the second information comprises a detonation quantity, and the third information comprises a first detonation velocity acquired by the first detection point and a first detonation velocity acquired by the second detection point;
In step S1, in the step of obtaining the first information, if the step of calculating the first explosive distance is S11, specific steps of calculating S11 include:
step S111: acquiring the horizontal distance from the explosion center of the first explosion to a first detection point, the tangential distance from the explosion center of the first explosion to the axis of the wind well and the height value of the first height;
step S112: calculating according to the horizontal distance from the explosion center of the first explosion to the first detection point, the tangential distance from the explosion center of the first explosion to the axis of the wind well and the height value of the first height to obtain a first explosion center distance;
wherein the height value of the first height is set asLocation of thenThe first explosive distance is calculated as follows:
in the above formula, x represents the horizontal distance from the explosion center of the first explosion to the first detection point, y represents the tangential distance from the explosion center of the first explosion to the axis of the wind well, h represents the height value of the first height,representing the positionFirst explosive distance (position)The range of values of (1) to (n). The positions of x, y and h can be seen in fig. 2 and 3.
If the second explosive distance is calculated in the step S12, the specific calculation in the step S12 includes:
step S121: acquiring the horizontal distance from the explosion center of the first explosion to a first detection point and the diameter of the air shaft;
Step S122: according to the horizontal distance from the explosion center of the first explosion to the first detection point and the diameter calculation of the air shaft, the horizontal distance from the explosion center of the first explosion to the second detection point is obtained;
the calculation formula of the horizontal distance from the explosion center of the first explosion to the second detection point is as follows:
in the above formula, m is the horizontal distance from the explosion center of the first explosion to the second detection point,representing the horizontal distance of the burst core of the first burst to the first detection point,indicating the diameter of the wind shaft. Wherein the positions of x and D can be seen in fig. 2 and 3.
Step S123: the tangential distance from the explosion center of the first explosion to the axis of the wind well and the height value of the first height are obtained;
step S124: calculating according to the horizontal distance from the explosion center of the first explosion to the second detection point, the tangential distance from the explosion center of the first explosion to the axis of the wind well and the height value of the first height to obtain a second explosion center distance;
in step S124, the height value of the first height is set to beLocation of thenThe calculation formula of the second explosive center distance is as follows:
in the above formula, m is the horizontal distance from the explosion center of the first explosion to the second detection point, y is the tangential distance from the explosion center of the first explosion to the axis of the wind well, h represents the height value of the first height, Representing the positionThe second explosive distance (position)The range of values of (1) to (n). The positions of y and h can be seen in fig. 2 and 3.
It is understood that the detonation amount is determined according to engineering implementation requirements; and the third information is acquired through a blasting vibration monitor.
S2, fitting calculation is carried out according to the first explosive distance, the second information and the first explosive vibration speed acquired by the first detection point, so as to obtain a first fitting coefficient and a second fitting coefficient;
in the fitting calculation, the first fitting coefficient and the second fitting coefficient can be calculated through a Sargassy formula, wherein the Sargassy formula is as follows:
in the above-mentioned method, the step of,the blasting vibration velocity acquired at the position i of the first detection point under the first blasting is represented,representing the positionThe first explosive distance (the range of the position i is 1 to n), Q is the explosive quantity, K is a first fitting coefficient,and is the second fitting coefficient.
Step S3, calculating according to the second explosive distance, the second information, the first fitting coefficient and the second fitting coefficient to obtain fourth information, wherein the fourth information comprises the theoretical vibration speed of the explosion of a second detection point under the first explosion;
In step S3, the calculation formula is:
in the above-mentioned method, the step of,the theoretical vibration of blasting, where the second detection point is located at the position j under the first blasting, K is a first fitting coefficient,as a result of the second fitting coefficient,the second explosive distance of the position j is represented (the value range of the position j is 1 to n), and Q is the explosive quantity.
And S4, calculating according to the first blasting vibration speed and fourth information acquired by the second detection point to obtain fifth information, wherein the fifth information is the first blasting vibration reduction rate corresponding to the first height.
In step S4, the calculation formula is:
in the above-mentioned method, the step of,the theoretical vibration velocity of the blast at the position j of the second detection point under the first blast is shown,the explosion vibration speed acquired at the position j of the second detection point under the first explosion is represented,the first shot damping rate corresponding to the first height at the first shot is shown.
After step S4, to further ensure the accuracy of the first high vibration reduction rate, steps S5 to S9 are further included:
s5, obtaining sixth information and seventh information, wherein the sixth information comprises a third explosion center distance and a fourth explosion center distance under the second explosion, the third explosion center distance is the explosion center distance from the explosion center of the second explosion to the first detection point, the fourth explosion center distance is the explosion center distance from the explosion center of the second explosion to the second detection point, and the seventh information comprises a second explosion vibration speed acquired by the first detection point and a second explosion vibration speed acquired by the second detection point;
In step S5, in the step of obtaining the sixth information, if the step of calculating the third explosive distance is S51, specific steps of calculating S51 include:
step S511: acquiring the horizontal distance from the explosion center of the first explosion to the first detection point and the single explosion cycle footage;
step S512: calculating according to the horizontal distance from the explosion center of the first explosion to the first detection point and the single explosion cyclic footage to obtain the horizontal distance from the explosion center of the second explosion to the first detection point;
step S513, obtaining a tangential distance from the explosion center of the first explosion to the axis of the wind well and a height value of the first height;
step S514, calculating according to the tangential distance from the explosion center of the first explosion to the axis of the wind well, the height value of the first height and the horizontal distance from the explosion center of the second explosion to the first detection point to obtain a third explosion center distance;
the calculation formula of the third explosive center distance is as follows:
in the above formula, x represents the horizontal distance from the explosion core of the first explosion to the first detection point, L represents the single explosion cycle footage, y represents the tangential distance from the explosion core of the first explosion to the axis of the wind well, and h represents the first heightThe height value of the degree is set to,representing the positionThird explosive distance (position)The range of values of (1) to (n). The position of x, y, h, L can be seen in fig. 2 and 3, among others.
If the fourth explosive distance is calculated in step S52, the specific calculation steps in step S52 include:
step S521: acquiring the horizontal distance from the explosion center of the first explosion to a first detection point and the diameter of the air shaft;
step S522: according to the horizontal distance from the explosion center of the first explosion to the first detection point and the diameter calculation of the air shaft, the horizontal distance from the explosion center of the first explosion to the second detection point is obtained;
the calculation formula of the horizontal distance from the explosion center of the first explosion to the second detection point is as follows:
in the above formula, m is the horizontal distance from the explosion center of the first explosion to the second detection point, x is the horizontal distance from the explosion center of the first explosion to the first detection point, and D is the diameter of the wind well. Wherein the positions of x and D can be seen in fig. 2 and 3.
Step S523, obtaining a single blasting circulation footage;
step S524, calculating according to the horizontal distance from the explosion core of the first explosion to the second detection point and the single explosion cyclic footage to obtain the horizontal distance from the explosion core of the second explosion to the second detection point;
step S525, obtaining a tangential distance from a blasting center of the first blasting to an axis of the wind well and a height value of the first height;
and S526, calculating according to the horizontal distance from the explosion center of the second explosion to the second detection point, the tangential distance from the explosion center of the first explosion to the axis of the wind well and the height value of the first height, and obtaining a fourth explosion center distance.
In step S526, the calculation formula is:
in the above formula, m is the horizontal distance from the explosion center of the first explosion to the second detection point, L represents the single explosion cycle footage, y is the tangential distance from the explosion center of the first explosion to the axis of the wind well, h represents the height value of the first height,the fourth explosive center distance of the position j is represented (the value range of the position j is 1 to n). The positions of y, h and L can be seen in fig. 2 and 3.
It is understood that the detonation amount is determined according to engineering implementation requirements; the seventh information is obtained through a blasting vibration monitor.
S6, calculating according to the third explosive distance, the second information and the second blasting vibration speed acquired by the first detection point, and updating a first fitting coefficient and a second fitting coefficient according to a calculation result;
step S7, according to the fourth explosive distance, the second information, the first fitting coefficient and the second fitting coefficient, eighth information is obtained, wherein the eighth information comprises the theoretical vibration speed of blasting of a second detection point under the second blasting;
in step S7, the calculation formula is:
in the above-mentioned method, the step of,indicating that the second detection point is located at the position under the second blastingj the theoretical vibration velocity of the blasting is, As a first coefficient of fit,as a result of the second fitting coefficient,the fourth explosive distance of the position j is shown (the value range of the position j is 1 to n), and Q is the explosive quantity.
Step S8, calculating according to the second blasting vibration speed acquired by the second detection point and the eighth information to obtain ninth information, wherein the ninth information is a second blasting vibration reduction rate corresponding to the first height;
in the above-mentioned method, the step of,the theoretical vibration velocity of the second detection point at position j in the second blasting is shown,the blasting vibration speed acquired by the second detection point at the position j under the second blasting is represented,the second shot damping rate corresponding to the first height at the second shot is shown.
And step S9, calculating according to the fifth information and the ninth information to obtain tenth information, wherein the tenth information is the comprehensive vibration reduction rate of the first height.
In step S9, the calculation formula is:
representing the first shot damping rate corresponding to the first altitude at the first shot,representing the second shot damping rate corresponding to the first altitude at the second shot,representing the integrated damping rate of the first altitude.
When the first height is taken as a reference, a first explosion vibration reduction rate of a preset height is obtained through calculation, in the step S1, the first information further comprises a fifth explosion center distance and a sixth explosion center distance under the first explosion, the fifth explosion center distance is an explosion center distance from the explosion center of the first explosion to a third detection point, the sixth explosion center distance is an explosion center distance from the explosion center of the first explosion to a fourth detection point, the difference between the third distance of the third detection point and the fourth distance of the fourth detection point is a wind well diameter, wherein the third detection point and the fourth detection point are both arranged on a wind well wall, the height of the third detection point and the height of the fourth detection point are both preset heights, the preset heights are the multiple relation of the first height, the third distance is a horizontal distance from the explosion center of the first explosion to the third detection point, and the fourth distance is a horizontal distance from the explosion center of the first explosion to the fourth detection point;
The preset height is a multiple relation of the first height, and the multiple relation is set to be an integral multiple for conveniently calculating the vibration reduction rate of the uniform stratum; in order to facilitate calculation of the vibration reduction rate of the uneven stratum, the multiple relation is set to be non-integral multiple, and at the moment, the vibration reduction rate calculation of the complex stratum is guaranteed.
When the first explosion vibration reduction rate of the preset height is calculated, after step S9, steps S10 to S13 are further included:
s10, obtaining eleventh information, wherein the eleventh information comprises a first blasting vibration speed acquired by a third detection point and a first blasting vibration speed acquired by a fourth detection point;
step S11, fitting calculation is carried out according to the fifth explosive distance, the second information and the first explosion vibration speed acquired by the third detection point, so as to obtain a third fitting coefficient and a fourth fitting coefficient;
step S12, according to the sixth explosive distance, the second information, the third fitting coefficient and the fourth fitting coefficient, twelfth information is obtained through calculation, and the twelfth information comprises the theoretical vibration speed of blasting of a fourth detection point under the first blasting;
and S13, calculating according to the first blasting vibration speed acquired by the fourth detection point and the twelfth information to obtain thirteenth information, wherein the thirteenth information is the first blasting vibration reduction rate corresponding to the preset height.
In step S10 to step S13, the calculation principle is the same as that of step S1 to step S4, and will not be described in detail here.
In order to further ensure the accuracy of the preset high vibration reduction rate, in step S5, the sixth information further includes a seventh core distance and an eighth core distance under the second blasting, where the seventh core distance is a core distance from the core of the second blasting to the third detection point, and the eighth core distance is a core distance from the core of the second blasting to the fourth detection point;
after the first blasting vibration reduction rate corresponding to the preset height is calculated, the method further comprises the steps of S14-S18:
s14, acquiring fourteenth information, wherein the fourteenth information comprises a second blasting vibration speed acquired by a third detection point and a second blasting vibration speed acquired by a fourth detection point;
s15, calculating according to the seventh explosion center distance, the second information and the second explosion vibration speed acquired by the third detection point, and updating a third fitting coefficient and a fourth fitting coefficient according to a calculation result;
s16, calculating according to the eighth explosive distance, the second information, the third fitting coefficient and the fourth fitting coefficient to obtain fifteenth information, wherein the fifteenth information comprises the theoretical vibration speed of blasting at a fourth detection point in the second blasting;
Step S17, according to the second blasting vibration velocity and fifteenth information acquired by the fourth detection point, sixteenth information is obtained through calculation, wherein the sixteenth information is a second blasting vibration damping rate corresponding to a preset height;
and step S18, calculating according to the thirteenth information and the sixteenth information to obtain seventeenth information, wherein the seventeenth information is the comprehensive vibration reduction rate of a preset height.
In steps S14 to S18, the calculation principle is the same as that of steps S5 to S9, and will not be described in detail here.
When the comprehensive vibration reduction rate of the stratum is researched from the angles of the first blasting and the second blasting after the comprehensive vibration reduction rate of the preset height is calculated, the step S18 further comprises the steps S19 to S21:
step S19, calculating according to the first blasting vibration reduction rate corresponding to the first height and the first blasting vibration reduction rate corresponding to the preset height to obtain a first blasting average vibration reduction rate;
step S20, calculating according to the second blasting vibration reduction rate corresponding to the first height and the second blasting vibration reduction rate corresponding to the preset height to obtain a second blasting average vibration reduction rate;
step S21: and calculating according to the first blasting average vibration reduction rate and the second blasting average vibration reduction rate to obtain the comprehensive vibration reduction rate.
When the preset heights are multiple groups, and the detection points of the multiple preset heights are equidistantly arranged, the calculation formula of the fifth explosive center distance is as follows:
in the above formula, x represents the horizontal distance from the explosion center of the first explosion to the first detection point, y represents the tangential distance from the explosion center of the first explosion to the axis of the wind well, h represents the height value of the first height, h 2 The height value of the equidistant measuring point on the same side is represented, i represents the measuring point position (positionThe value of (2) is in the range of 1 to n),the fifth explosive center distance of the position i is represented (the value range of the position i is 1 to n). Wherein x, h, and h 2 The position can be seen in fig. 4.
The calculation formula of the sixth explosive core distance is as follows:
in the above formula, m is the horizontal distance from the explosion center of the first explosion to the second detection point, y is the tangential distance from the explosion center of the first explosion to the axis of the wind well, h represents the height value of the first height, h 2 The height value of the equidistant measuring point on the same side is represented, j represents the measuring point position of the fourth measuring point (the value range of the position j is 1 to n),representing the positionThe sixth explosive core distance (the value range of the position j is 1 to n). Wherein h and h 2 The position can be seen in fig. 4.
The calculation formula of the seventh explosion center distance is as follows:
in the above formula, x represents the horizontal distance from the explosion center of the first explosion to the first detection point, L represents the single explosion cycle footage, y represents the tangential distance from the explosion center of the first explosion to the axis of the wind well, h represents the height value of the first height, and h 2 The height value of the equidistant measuring point on the same side is represented, i represents the measuring point position of the third measuring point (the range of the position i is 1 to n),the seventh explosive center distance of the position i is represented (the value range of the position i is 1 to n). Wherein x, h and h 2 The position can be seen in fig. 4.
The calculation formula of the eighth explosive center distance is as follows:
in the above formula, m is the horizontal distance from the explosion center of the first explosion to the second detection point, L represents the single explosion cycle footage, y represents the tangential distance from the explosion center of the first explosion to the axis of the wind well, h represents the height value of the first height, and h 2 The height value of the equidistant measuring point on the same side is represented, j represents the measuring point position of the fourth measuring point (the value range of the position j is 1 to n),the eighth explosive center distance of the position j is represented (the value range of the position j is 1 to n). Wherein x, h and h 2 The position can be seen in fig. 4.
When the preset heights are multiple groups and the detection points of the multiple preset heights are equidistantly arranged, in the fitting calculation process, the point values of multiple groups of third fitting coefficients and the point values of multiple groups of fourth fitting coefficients can be obtained through a Sarkowski formula, and then the point values of the multiple groups of third fitting coefficients and the point values of the multiple groups of fourth fitting coefficients are fitted through a least square method, so that the third fitting coefficients and the fourth fitting coefficients are obtained. If there is a layer A The average vibration reduction rate of the first blasting of the stratum A is as follows:
in the above-mentioned method, the step of,the average damping rate for the first shot of formation a is indicated,the damping rate at the first burst detection point k is shown, and k is the detection point.
Similarly, the average vibration reduction rate of the second blasting of the stratum A is obtained as follows:
in the above-mentioned method, the step of,the average vibration reduction rate for the second shot of formation a is shown,the damping rate at the detection point k in the second explosion is shown, and k is the detection point.
The damping rate of formation a is:
in the above-mentioned method, the step of,the average damping rate for the first shot of formation a is indicated,the average vibration reduction rate for the second shot of formation a is shown,the damping rate of formation a is shown.
As specific engineering examples:
TABLE 1
TABLE 2
The damping rate can be obtained:
in order to correct the first burst damping rate corresponding to the first height, after step S21, steps S22 to S26 are further included:
step S22, eighteenth information and nineteenth information are acquired, the eighteenth information comprises a ninth explosion center distance and a tenth explosion center distance under the first explosion, the ninth explosion center distance is from an explosion center of the first explosion to an explosion center distance of a fifth detection point, the tenth explosion center distance is from an explosion center of the first explosion to an explosion center distance of a sixth detection point, the fifth detection point is arranged on the same side of the first detection point, the fifth detection point is determined through a preset horizontal distance with the first detection point, the sixth detection point is arranged on the same side of the second detection point, the sixth detection point is determined through a preset horizontal distance with the second detection point, the height of the fifth detection point and the height of the sixth detection point are both the first height, and the nineteenth information comprises a first explosion vibration speed acquired by the fifth detection point and a first explosion vibration speed acquired by the sixth detection point;
The preset horizontal distance is 0cm to 20cm so as to reduce damage to the stratum.
S23, calculating according to the ninth explosive distance, the second information and the first explosion vibration speed acquired by the fifth detection point, and updating a first fitting coefficient and a second fitting coefficient according to a calculation result;
step S24: according to the tenth explosion center distance, the second information, the first fitting coefficient and the second fitting coefficient, obtaining twentieth information, wherein the twentieth information comprises the theoretical vibration speed of the explosion of a sixth detection point under the first explosion;
step S25: calculating according to the first blasting vibration speed and the twentieth information acquired by the sixth detection point to obtain a first blasting vibration reduction rate correction value corresponding to the first height;
step S26: and calculating according to the first explosion vibration reduction rate correction value corresponding to the first height and fifth information to obtain the first explosion vibration reduction rate of the first height correction.
Example 2:
as shown in fig. 6, the present embodiment provides a vibration damping rate calculation system of a vibration damping hole, which includes an acquisition module 901, a first calculation module 902, a second calculation module 903, and a first processing module 904, wherein,
The acquiring module 901 is configured to acquire first information, second information and third information, where the first information includes a first shot distance and a second shot distance under the first shot, the first shot distance is a shot distance from the first shot to the first detection point, the second shot distance is a shot distance from the first shot to the second detection point, a difference between the first distance of the first detection point and the second distance of the second detection point is a wind shaft diameter, the first detection point and the second detection point are both disposed on a wind shaft wall, the first detection point and the second detection point are both at a first height, the first distance is a horizontal distance from the first shot to the first detection point, the second distance is a horizontal distance from the first shot to the second detection point, the second information includes a detonation velocity, and the third information includes a first detonation velocity acquired by the first detection point and a first detonation velocity acquired by the second detection point;
the first calculation module 902 is configured to perform fitting calculation according to the first heart rate distance, the second information, and the first blasting vibration velocity acquired by the first detection point, to obtain a first fitting coefficient and a second fitting coefficient;
The second calculation module 903 is configured to calculate, according to the second center of burst distance, the second information, the first fitting coefficient, and the second fitting coefficient, obtain fourth information, where the fourth information includes a theoretical vibration velocity of blasting at a second detection point during the first blasting;
the first processing module 904 is configured to calculate, according to the first blasting vibration speed and the fourth information acquired by the second detection point, obtain fifth information, where the fifth information is a first blasting vibration reduction rate corresponding to the first height.
In one embodiment of the present disclosure, the acquisition module 901 includes a first acquisition unit 9011 and a first calculation unit 9012, wherein,
a first obtaining unit 9011, configured to obtain a horizontal distance from a center of explosion of the first explosion to a first detection point, a tangential distance from the center of explosion of the first explosion to an axis of the wind well, and a height value of the first height;
the first calculating unit 9012 is configured to calculate, according to the horizontal distance from the explosion center of the first explosion to the first detection point, the tangential distance from the explosion center of the first explosion to the axis of the wind well, and the height value of the first height, obtain a first explosion center distance;
in one embodiment of the present disclosure, after the first computing unit 9012, a second acquiring unit 9013, a second computing unit 9014, a third acquiring unit 9015, and a third computing unit 9016 are further included, wherein,
A second obtaining unit 9013, configured to obtain a horizontal distance from the explosion center of the first explosion to the first detection point and a diameter of the wind shaft;
a second calculating unit 9014, configured to calculate, according to the horizontal distance from the explosion center of the first explosion to the first detection point and the diameter of the air shaft, obtain the horizontal distance from the explosion center of the first explosion to the second detection point;
a third obtaining unit 9015, configured to obtain a tangential distance from the explosion center of the first explosion to the axis of the wind well and a height value of the first height;
a third calculating unit 9016, configured to calculate, according to the horizontal distance from the first explosion center to the second detection point, the tangential distance from the first explosion center to the axis of the wind well, and the height value of the first height, to obtain a second explosion center distance;
in a specific embodiment of the present disclosure, after the first processing module 904, a second processing module 905 is further included, as shown in fig. 7, where the second processing module 905 includes a fourth acquiring unit 9051, a ninth calculating unit 9052, a tenth calculating unit 9053, an eleventh calculating unit 9054, and a twelfth calculating unit 9055, where,
a fourth obtaining unit 9051, configured to obtain sixth information and seventh information, where the sixth information includes a third shot distance and a fourth shot distance under the second shot, the third shot distance is a shot distance from a shot of the second shot to the first detection point, the fourth shot distance is a shot distance from a shot of the second shot to the second detection point, and the seventh information includes a second shot vibration velocity acquired by the first detection point and a second shot vibration velocity acquired by the second detection point;
A ninth calculating unit 9052, configured to calculate, according to the third center distance, the second information, and the second blasting vibration velocity acquired by the first detection point, and update a first fitting coefficient and a second fitting coefficient according to a calculation result;
a tenth calculation unit 9053, configured to calculate, according to the fourth heart distance, the second information, the first fitting coefficient, and the second fitting coefficient, obtain eighth information, where the eighth information includes a theoretical vibration velocity of blasting at a second detection point during a second blasting;
an eleventh calculating unit 9054, configured to calculate, according to the second blasting vibration velocity acquired by the second detection point and the eighth information, to obtain ninth information, where the ninth information is a second blasting vibration damping rate corresponding to the first height;
and a twelfth calculating unit 9055, configured to calculate tenth information according to the fifth information and the ninth information, where the tenth information is the integrated damping rate of the first height.
In one embodiment of the present disclosure, the fourth acquiring unit 9051 includes a fifth acquiring unit 90511, a fourth calculating unit 90512, a sixth acquiring unit 90513, and a fifth calculating unit 90514, wherein,
A fifth acquiring unit 90511, configured to acquire a horizontal distance from the explosion center of the first explosion to the first detection point and a single explosion cycle footage;
a fourth calculating unit 90512, configured to calculate, according to the horizontal distance from the center of the first shot to the first detection point and the single shot cyclic footage, obtain the horizontal distance from the center of the second shot to the first detection point;
a sixth acquiring unit 90513, configured to acquire a tangential distance from a center of explosion of the first blasting to an axis of the wind shaft and a height value of the first height;
a fifth calculating unit 90514, configured to calculate, according to the tangential distance from the explosion center of the first explosion to the axis of the wind well, the height value of the first height, and the horizontal distance from the explosion center of the second explosion to the first detection point, a third explosion center distance;
in one embodiment of the present disclosure, after the fifth computing unit 90514, a seventh acquiring unit 90515, a sixth computing unit 90516, an eighth acquiring unit 90517, a seventh computing unit 90518, a ninth acquiring unit 90519, and an eighth computing unit 90520 are further included, wherein,
a seventh acquiring unit 90515, configured to acquire a horizontal distance from the explosion center of the first explosion to the first detection point and a diameter of the blast well;
A sixth calculating unit 90516, configured to calculate, according to the horizontal distance from the explosion center of the first explosion to the first detection point and the diameter of the air shaft, obtain the horizontal distance from the explosion center of the first explosion to the second detection point;
an eighth acquiring unit 90517 for acquiring a single blasting cycle footage;
a seventh calculating unit 90518, configured to calculate, according to the horizontal distance from the center of the first shot to the second detection point and the single shot cyclic footage, obtain the horizontal distance from the center of the second shot to the second detection point;
a ninth acquiring unit 90519, configured to acquire a tangential distance from a blasting center of the first blasting to an axis of the wind shaft and a height value of the first height;
and an eighth calculating unit 90520, configured to calculate, according to the horizontal distance from the explosion center of the second explosion to the second detection point, the tangential distance from the explosion center of the first explosion to the axis of the wind well, and the height value of the first height, obtain a fourth explosion center distance.
When calculating the first burst damping rate of the preset height, in one embodiment of the present disclosure, after the second processing module 905, a third processing module 906 is further included, where the third processing module 906 includes a tenth acquisition unit 9060, a thirteenth calculation unit 9061, a fourteenth calculation unit 9062 and a fifteenth calculation unit 9063, where,
A tenth acquiring unit 9060, configured to acquire eleventh information, where the eleventh information includes a first blasting vibration speed acquired by the third detecting point and a first blasting vibration speed acquired by the fourth detecting point;
a thirteenth calculation unit 9061, configured to obtain a third fitting coefficient and a fourth fitting coefficient according to the fifth heart distance, the second information, and the first blasting vibration velocity acquired by the third detection point;
a fourteenth calculation unit 9062, configured to calculate, according to the sixth explosive distance, the second information, the third fitting coefficient, and the fourth fitting coefficient, obtain twelfth information, where the twelfth information includes a theoretical vibration velocity of blasting at a fourth detection point under the first blasting;
a fifteenth calculating unit 9063, configured to calculate, according to the first blasting vibration velocity acquired by the fourth detection point and the twelfth information, thirteenth information, where the thirteenth information is a first blasting vibration damping rate corresponding to a preset height.
In a specific embodiment of the present disclosure, after the third processing module 906, a fourth processing module 907 is further included, the fourth processing module 907 including an eleventh obtaining unit 9070, a sixteenth calculating unit 9071, a seventeenth calculating unit 9072, an eighteenth calculating unit 9073, and a nineteenth calculating unit 9074, wherein,
An eleventh acquiring unit 9070, configured to acquire fourteenth information, where the fourteenth information includes a second blasting vibration speed acquired by the third detecting point and a second blasting vibration speed acquired by the fourth detecting point;
a sixteenth calculating unit 9071, configured to calculate, according to the seventh heart rate distance, the second information, and the second blasting vibration velocity acquired by the third detection point, and update a third fitting coefficient and a fourth fitting coefficient according to a calculation result;
a seventeenth calculating unit 9072, configured to calculate, according to the eighth heart distance, the second information, the third fitting coefficient, and the fourth fitting coefficient, obtain fifteenth information, where the fifteenth information includes a theoretical vibration velocity of blasting at a fourth detection point during the second blasting;
an eighteenth calculating unit 9073, configured to calculate, according to the second blasting vibration velocity and the fifteenth information acquired by the fourth detecting point, obtain sixteenth information, where the sixteenth information is a second blasting vibration damping rate corresponding to a preset height;
a nineteenth calculating unit 9074, configured to calculate, according to the thirteenth information and the sixteenth information, seventeenth information, where the seventeenth information is a comprehensive vibration reduction rate of a preset height.
In one embodiment of the present disclosure, after the fourth processing module 907, a fifth processing module 908 is further included, where the fifth processing module 908 includes a twentieth computing unit 9081, a twenty-first computing unit 9082, and a twenty-second computing unit 9083, where,
a twentieth calculation unit 9081, configured to calculate, according to the first blasting vibration reduction rate corresponding to the first height and the first blasting vibration reduction rate corresponding to the preset height, a first blasting average vibration reduction rate;
a twenty-first calculation unit 9082, configured to calculate, according to the second blasting vibration reduction rate corresponding to the first height and the second blasting vibration reduction rate corresponding to the preset height, a second blasting average vibration reduction rate;
and a twenty-second calculation unit 9083, configured to calculate, according to the first blasting average damping rate and the second blasting average damping rate, a comprehensive damping rate.
In a specific embodiment of the present disclosure, after the fifth processing module 908, a sixth processing module 909 is further included, where the sixth processing module 909 includes a twelfth acquisition unit 9091, a twenty-third calculation unit 9092, a twenty-fourth calculation unit 9093, a twenty-fifth calculation unit 9094, and a twenty-sixth calculation unit 9095, where,
A twelfth obtaining unit 9091, configured to obtain eighteenth information and nineteenth information, where the eighteenth information includes a ninth shot center distance and a tenth shot center distance under the first shot, the ninth shot center distance is a shot center distance from the first shot to a fifth detection point, the tenth shot center distance is a shot center distance from the first shot to a sixth detection point, the fifth detection point is disposed on the same side as the first detection point, the fifth detection point is determined by a preset horizontal distance from the first detection point, the sixth detection point is disposed on the same side as the second detection point, the sixth detection point is determined by a preset horizontal distance from the second detection point, and the height of the fifth detection point and the height of the sixth detection point are both the first height, and the nineteenth information includes a first shot vibration speed collected by the fifth detection point and a first shot vibration speed collected by the sixth detection point;
a twenty-third calculation unit 9092, configured to calculate, according to the ninth explosive distance, the second information, and the first explosion vibration velocity acquired by the fifth detection point, and update a first fitting coefficient and a second fitting coefficient according to a calculation result;
a twenty-fourth calculating unit 9093, configured to calculate, according to the tenth explosion center distance, the second information, the first fitting coefficient, and the second fitting coefficient, obtain twentieth information, where the twentieth information includes a theoretical vibration velocity of the sixth detection point during the first explosion;
A twenty-fifth calculating unit 9094, configured to calculate, according to the first blasting vibration velocity acquired by the sixth detection point and the twentieth information, to obtain a first blasting vibration damping rate correction value corresponding to the first height;
and a twenty-sixth calculating unit 9095, configured to calculate, according to the first explosion vibration damping rate correction value corresponding to the first height and the fifth information, to obtain a first explosion vibration damping rate with the first height corrected.
Example 3:
corresponding to the above method embodiment, there is also provided a vibration damping rate calculating apparatus of a vibration damping hole in the present embodiment, and the vibration damping rate calculating apparatus of a vibration damping hole described below and the vibration damping rate calculating method of a vibration damping hole described above may be referred to in correspondence with each other.
FIG. 5 is a block diagram of a damper ratio calculation device 800 of a damper orifice shown in accordance with an exemplary embodiment. As shown in fig. 5, the damping rate calculation device 800 of the damping hole may include: a processor 801, a memory 802. The damping rate computing device 800 of the damping orifice may also include one or more of a multimedia component 803, an i/O interface 804, and a communication component 805.
Wherein the processor 801 is configured to control the overall operation of the damping-rate calculation apparatus 800 for the damper orifice to perform all or part of the steps in the damping-rate calculation method for the damper orifice described above. The memory 802 is used to store various types of data to support the operation of the damper rate computing device 800 at the damper orifice, which may include, for example, instructions for any application or method operating on the damper rate computing device 800 at the damper orifice, as well as application related data such as contact data, messages, pictures, audio, video, and the like. The Memory 802 may be implemented by any type or combination of volatile or non-volatile Memory devices, such as static random access Memory (Static Random Access Memory, SRAM for short), electrically erasable programmable Read-Only Memory (Electrically Erasable Programmable Read-Only Memory, EEPROM for short), erasable programmable Read-Only Memory (Erasable Programmable Read-Only Memory, EPROM for short), programmable Read-Only Memory (Programmable Read-Only Memory, PROM for short), read-Only Memory (ROM for short), magnetic Memory, flash Memory, magnetic disk, or optical disk. The multimedia component 803 may include a screen and an audio component. Wherein the screen may be, for example, a touch screen, the audio component being for outputting and/or inputting audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signals may be further stored in the memory 802 or transmitted through the communication component 805. The audio assembly further comprises at least one speaker for outputting audio signals. The I/O interface 804 provides an interface between the processor 801 and other interface modules, which may be a keyboard, mouse, buttons, etc. These buttons may be virtual buttons or physical buttons. The communication component 805 is configured to perform wired or wireless communication between the vibration damping rate computing device 800 of the vibration damping hole and other devices. Wireless communication, such as Wi-Fi, bluetooth, near field communication (Near FieldCommunication, NFC for short), 2G, 3G or 4G, or a combination of one or more thereof, the respective communication component 805 may thus comprise: wi-Fi module, bluetooth module, NFC module.
In an exemplary embodiment, the damper ratio calculation device 800 of the damper orifice may be implemented by one or more application specific integrated circuits (Application Specific Integrated Circuit, abbreviated as ASIC), digital signal processor (DigitalSignal Processor, abbreviated as DSP), digital signal processing device (Digital Signal Processing Device, abbreviated as DSPD), programmable logic device (Programmable Logic Device, abbreviated as PLD), field programmable gate array (Field Programmable Gate Array, abbreviated as FPGA), controller, microcontroller, microprocessor, or other electronic component for performing the damper ratio calculation method of the damper orifice described above.
In another exemplary embodiment, there is also provided a computer readable storage medium including program instructions which, when executed by a processor, implement the steps of the damping rate calculation method of a damping hole described above. For example, the computer readable storage medium may be the memory 802 described above including program instructions executable by the processor 801 of the damping-hole damping-rate calculation device 800 to perform the damping-hole damping-rate calculation method described above.
Example 4:
corresponding to the above method embodiment, there is also provided a readable storage medium in this embodiment, and a readable storage medium described below and a method for calculating a damping rate of a damping hole described above may be referred to correspondingly to each other.
A readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the damping ratio calculation method of the damping hole of the above method embodiment.
The readable storage medium may be a usb disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), a magnetic disk, or an optical disk, and the like.

Claims (10)

1. The method for calculating the damping rate of the damping hole is characterized by comprising the following steps:
acquiring first information, second information and third information, wherein the first information comprises a first core distance and a second core distance under the first blasting, the first core distance is a core distance from the core of the first blasting to a first detection point, the second core distance is a core distance from the core of the first blasting to a second detection point, the difference between the first distance of the first detection point and the second distance of the second detection point is a wind well diameter, the first detection point and the second detection point are both arranged on the wind well wall, the height of the first detection point and the height of the second detection point are both the first height, the first distance is a horizontal distance from the core of the first blasting to the first detection point, the second distance is a horizontal distance from the core of the first blasting to the second detection point, the second information comprises an initiating explosive amount, and the third information comprises a first blasting vibration speed acquired by the first detection point and a first blasting vibration speed acquired by the second detection point;
According to the first explosive center distance, the second information and the first explosive vibration speed acquired by the first detection point, fitting and calculating to obtain a first fitting coefficient and a second fitting coefficient;
calculating according to the second explosive distance, the second information, the first fitting coefficient and the second fitting coefficient to obtain fourth information, wherein the fourth information comprises the theoretical vibration speed of blasting of a second detection point under the first blasting;
and calculating according to the first blasting vibration speed and fourth information acquired by the second detection point to obtain fifth information, wherein the fifth information is the first blasting vibration reduction rate corresponding to the first height.
2. The vibration reduction rate calculation method of a vibration reduction hole according to claim 1, characterized in that: according to the first blasting vibration speed acquired by the second detection point and the fourth information, calculating to obtain fifth information, and then, including:
obtaining sixth information and seventh information, wherein the sixth information comprises a third explosion center distance and a fourth explosion center distance under the second explosion, the third explosion center distance is the explosion center distance from the explosion center of the second explosion to the first detection point, the fourth explosion center distance is the explosion center distance from the explosion center of the second explosion to the second detection point, and the seventh information comprises a second explosion vibration speed acquired by the first detection point and a second explosion vibration speed acquired by the second detection point;
Calculating according to the third explosive distance, the second information and the second blasting vibration speed acquired by the first detection point, and updating a first fitting coefficient and a second fitting coefficient according to a calculation result;
calculating according to the fourth explosive distance, the second information, the first fitting coefficient and the second fitting coefficient to obtain eighth information, wherein the eighth information comprises the theoretical vibration speed of blasting of a second detection point under the second blasting;
according to the second blasting vibration speed acquired by the second detection point and the eighth information, obtaining ninth information, wherein the ninth information is a second blasting vibration reduction rate corresponding to the first height;
and calculating according to the fifth information and the ninth information to obtain tenth information, wherein the tenth information is the comprehensive vibration reduction rate of the first height.
3. The vibration reduction rate calculation method of a vibration reduction hole according to claim 1, characterized in that: the first information further comprises a fifth core distance and a sixth core distance under the first blasting, wherein the fifth core distance is a core distance from the core of the first blasting to a third detection point, the sixth core distance is a core distance from the core of the first blasting to a fourth detection point, the difference between the third distance of the third detection point and the fourth distance of the fourth detection point is a wind shaft diameter, the third detection point and the fourth detection point are both arranged on the wall of the wind shaft, the height of the third detection point and the height of the fourth detection point are both preset heights, the preset heights are a multiple relation of the first height, the third distance is a horizontal distance from the core of the first blasting to the third detection point, and the fourth distance is a horizontal distance from the core of the first blasting to the fourth detection point;
After the fifth information is obtained by calculating according to the first blasting vibration speed and the fourth information acquired by the second detection point, the method further comprises the following steps:
acquiring eleventh information, wherein the eleventh information comprises a first blasting vibration speed acquired by a third detection point and a first blasting vibration speed acquired by a fourth detection point;
according to the fifth explosive center distance, the second information and the first explosion vibration speed acquired by the third detection point, fitting and calculating to obtain a third fitting coefficient and a fourth fitting coefficient;
according to the sixth explosive distance, the second information, the third fitting coefficient and the fourth fitting coefficient, twelfth information is obtained through calculation, and the twelfth information comprises the theoretical vibration speed of blasting of a fourth detection point under the first blasting;
and calculating according to the first blasting vibration speed acquired by the fourth detection point and the twelfth information to obtain thirteenth information, wherein the thirteenth information is the first blasting vibration reduction rate corresponding to a preset height.
4. A vibration reduction ratio calculation method of a vibration reduction hole according to any one of claims 1 to 3, characterized in that: according to the first blasting vibration speed and fourth information acquired by the second detection point, calculating to obtain fifth information, wherein the fifth information comprises the following steps:
The method comprises the steps of obtaining eighteenth information and nineteenth information, wherein the eighteenth information comprises a ninth explosion center distance and a tenth explosion center distance under first explosion, the ninth explosion center distance is from an explosion center of the first explosion to a tenth explosion center distance from the explosion center of the first explosion to a sixth explosion center distance, the tenth explosion center distance is from the explosion center of the first explosion to the explosion center of a sixth explosion center distance, the fifth explosion center distance is arranged on the same side of a first detection point, the fifth detection point is determined through a preset horizontal distance from the first detection point, the sixth detection point is arranged on the same side of a second detection point, the sixth detection point is determined through a preset horizontal distance from the second detection point, the height of the fifth detection point and the height of the sixth detection point are both the first height, and the nineteenth information comprises a first explosion vibration speed acquired by the fifth detection point and a first explosion vibration speed acquired by the sixth detection point;
calculating according to the ninth explosive distance, the second information and the first explosion vibration speed acquired by the fifth detection point, and updating a first fitting coefficient and a second fitting coefficient according to a calculation result;
according to the tenth explosion center distance, the second information, the first fitting coefficient and the second fitting coefficient, obtaining twentieth information, wherein the twentieth information comprises the theoretical vibration speed of the explosion of a sixth detection point under the first explosion;
Calculating according to the first blasting vibration speed and the twentieth information acquired by the sixth detection point to obtain a first blasting vibration reduction rate correction value corresponding to the first height;
and calculating according to the first explosion vibration reduction rate correction value corresponding to the first height and fifth information to obtain the first explosion vibration reduction rate of the first height correction.
5. A vibration damping rate calculation system for a vibration damping hole, comprising:
the device comprises an acquisition module, a control module and a control module, wherein the acquisition module is used for acquiring first information, second information and third information, the first information comprises a first core distance and a second core distance under the first blasting, the first core distance is a core distance from the first core to the first detection point, the second core distance is a core distance from the first core to the second detection point, the difference between the first distance of the first detection point and the second distance of the second detection point is a wind shaft diameter, the first detection point and the second detection point are both arranged on the wind shaft wall, the height of the first detection point and the height of the second detection point are both the first height, the first distance is a horizontal distance from the core to the first detection point, the second distance is a horizontal distance from the core to the second detection point, the second information comprises a detonation quantity, and the third information comprises a first detonation velocity acquired by the first detection point and a first detonation velocity acquired by the first detection point;
The first calculation module is used for obtaining a first fitting coefficient and a second fitting coefficient according to the first explosive distance, the second information and the first explosion vibration speed acquired by the first detection point;
the second calculation module is used for calculating according to the second explosive distance, the second information, the first fitting coefficient and the second fitting coefficient to obtain fourth information, wherein the fourth information comprises the theoretical vibration speed of blasting of a second detection point under the first blasting;
the first processing module is used for calculating according to the first blasting vibration speed and the fourth information acquired by the second detection point to obtain fifth information, wherein the fifth information is the first blasting vibration reduction rate corresponding to the first height.
6. The vibration reduction ratio calculation system of a vibration reduction hole according to claim 5, further comprising a second processing module after the first processing module, the second processing module comprising:
a fourth obtaining unit, configured to obtain sixth information and seventh information, where the sixth information includes a third shot distance and a fourth shot distance under the second shot, the third shot distance is a shot distance from a shot of the second shot to the first detection point, the fourth shot distance is a shot distance from a shot of the second shot to the second detection point, and the seventh information includes a second shot vibration velocity acquired by the first detection point and a second shot vibration velocity acquired by the second detection point;
A ninth calculation unit, configured to calculate, according to the third explosive distance, the second information, and the second blasting vibration velocity acquired by the first detection point, and update a first fitting coefficient and a second fitting coefficient according to a calculation result;
a tenth calculation unit, configured to calculate, according to the fourth explosive distance, the second information, the first fitting coefficient, and the second fitting coefficient, obtain eighth information, where the eighth information includes an explosion theoretical vibration velocity of the second detection point under the second explosion;
an eleventh calculation unit, configured to calculate, according to the second blasting vibration velocity acquired by the second detection point and the eighth information, to obtain ninth information, where the ninth information is a second blasting vibration damping rate corresponding to the first height;
and a twelfth calculation unit, configured to calculate, according to the fifth information and the ninth information, tenth information, where the tenth information is a comprehensive vibration reduction rate of the first height.
7. The vibration reduction rate calculation system of a vibration reduction hole according to claim 5, further comprising a third processing module after the first processing module, wherein the first information further comprises a fifth heart distance and a sixth heart distance under the first explosion, the fifth heart distance is a heart distance from a heart of the first explosion to a third detection point, the sixth heart distance is a heart distance from the heart of the first explosion to a fourth detection point, a difference between the third distance of the third detection point and the fourth distance of the fourth detection point is a wind well diameter, the third detection point and the fourth detection point are both arranged on a wind well wall, the height of the third detection point and the height of the fourth detection point are both preset heights, the preset heights are a multiple relation of the first height, the third distance is a horizontal distance from the heart of the first explosion to the third detection point, and the fourth distance is a horizontal distance from the heart of the first explosion to the fourth detection point; the third processing module includes:
A tenth acquisition unit, configured to acquire eleventh information, where the eleventh information includes a first blasting vibration speed acquired by the third detection point and a first blasting vibration speed acquired by the fourth detection point;
a thirteenth calculation unit, configured to obtain a third fitting coefficient and a fourth fitting coefficient according to the fifth heart distance, the second information, and the first blasting vibration velocity acquired by the third detection point;
a fourteenth calculation unit, configured to calculate, according to the sixth explosive center distance, the second information, the third fitting coefficient, and the fourth fitting coefficient, obtain twelfth information, where the twelfth information includes a theoretical vibration velocity of blasting at a fourth detection point during the first blasting;
the fifteenth calculation unit is used for calculating according to the first blasting vibration speed acquired by the fourth detection point and the twelfth information to obtain thirteenth information, wherein the thirteenth information is the first blasting vibration reduction rate corresponding to a preset height.
8. The vibration reduction ratio calculation system of a vibration reduction hole according to claim 5, further comprising a sixth processing module after the first processing module, the sixth processing module comprising:
A twelfth obtaining unit, configured to obtain eighteenth information and nineteenth information, where the eighteenth information includes a ninth shot distance and a tenth shot distance under the first shot, the ninth shot distance is a shot distance from a shot of the first shot to a fifth detection point, the tenth shot distance is a shot distance from a shot of the first shot to a sixth detection point, the fifth detection point is disposed on the same side of the first detection point, the fifth detection point is determined by a preset horizontal distance from the first detection point, the sixth detection point is disposed on the same side of the second detection point, and the sixth detection point is determined by a preset horizontal distance from the second detection point, both the height of the fifth detection point and the height of the sixth detection point are the first height, and the nineteenth information includes a first shot vibration speed collected by the fifth detection point and a first shot vibration speed collected by the sixth detection point;
a twenty-third calculation unit, configured to calculate, according to the ninth explosive distance, the second information, and the first blasting vibration velocity acquired by the fifth detection point, and update a first fitting coefficient and a second fitting coefficient according to a calculation result;
a twenty-fourth calculation unit, configured to calculate, according to the tenth explosion center distance, the second information, the first fitting coefficient, and the second fitting coefficient, obtain twentieth information, where the twentieth information includes a theoretical vibration velocity of the sixth detection point during the first explosion;
A twenty-fifth calculation unit, configured to calculate, according to the first blasting vibration speed collected by the sixth detection point and the twentieth information, to obtain a first blasting vibration damping rate correction value corresponding to the first height;
and a twenty-sixth calculation unit, configured to calculate, according to the first explosion vibration reduction rate correction value corresponding to the first height and the fifth information, to obtain a first explosion vibration reduction rate with the first height corrected.
9. Damping rate calculation apparatus of damping hole, characterized by comprising:
a memory for storing a computer program;
a processor for implementing the steps of the vibration damping ratio calculation method of the vibration damping hole according to any one of claims 1 to 4 when executing the computer program.
10. A readable storage medium, wherein a computer program is stored on the readable storage medium, which when executed by a processor, implements the steps of the damping ratio calculation method of a damping hole according to any one of claims 1 to 4.
CN202211601008.4A 2022-12-14 2022-12-14 Vibration damping rate calculation method, system and equipment for vibration damping hole and readable storage medium Active CN116242209B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211601008.4A CN116242209B (en) 2022-12-14 2022-12-14 Vibration damping rate calculation method, system and equipment for vibration damping hole and readable storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211601008.4A CN116242209B (en) 2022-12-14 2022-12-14 Vibration damping rate calculation method, system and equipment for vibration damping hole and readable storage medium

Publications (2)

Publication Number Publication Date
CN116242209A CN116242209A (en) 2023-06-09
CN116242209B true CN116242209B (en) 2023-08-18

Family

ID=86628533

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211601008.4A Active CN116242209B (en) 2022-12-14 2022-12-14 Vibration damping rate calculation method, system and equipment for vibration damping hole and readable storage medium

Country Status (1)

Country Link
CN (1) CN116242209B (en)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06249599A (en) * 1993-02-26 1994-09-06 Sato Kogyo Co Ltd Expert system of blasting control
JP2006090573A (en) * 2004-09-21 2006-04-06 Shimizu Corp Method of estimating blasting vibration in controlled blasting, and blasting control method
JP2014222119A (en) * 2013-05-13 2014-11-27 清水建設株式会社 Method and apparatus for detecting detonation
CN104390537A (en) * 2014-09-23 2015-03-04 武汉大学 Method for controlling damage of slope pre-splitting blasting digging based on blast vibration test
CN108592720A (en) * 2018-05-16 2018-09-28 厦门爆破工程公司 A kind of vibration damping hole distribution method reducing ground blasting vibration
CN110736400A (en) * 2019-10-22 2020-01-31 长江重庆航道工程局 underwater drilling blasting vibration velocity calculation method considering rock internal structure
CN111426243A (en) * 2020-01-20 2020-07-17 重庆中环建设有限公司 Blasting parameter determination method based on differential vibration synthesis of different blast holes in near zone of blasting source
CN111780634A (en) * 2020-07-20 2020-10-16 中国铁建投资集团有限公司 Method for reducing blasting vibration velocity
CN113343503A (en) * 2021-07-26 2021-09-03 内蒙聚力工程爆破有限公司 Blasting method and device for predicting blasting vibration damage
CN113588064A (en) * 2021-07-06 2021-11-02 中铁二十局集团第六工程有限公司 Method for measuring blasting vibration data
WO2022051890A1 (en) * 2020-09-08 2022-03-17 中煤科工集团淮北爆破技术研究院有限公司 Vibration reduction method based on blast vibration prediction technology
CN114781040A (en) * 2022-05-10 2022-07-22 青岛理工大学 Building blasting vibration safety evaluation method, device, equipment and medium

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06249599A (en) * 1993-02-26 1994-09-06 Sato Kogyo Co Ltd Expert system of blasting control
JP2006090573A (en) * 2004-09-21 2006-04-06 Shimizu Corp Method of estimating blasting vibration in controlled blasting, and blasting control method
JP2014222119A (en) * 2013-05-13 2014-11-27 清水建設株式会社 Method and apparatus for detecting detonation
CN104390537A (en) * 2014-09-23 2015-03-04 武汉大学 Method for controlling damage of slope pre-splitting blasting digging based on blast vibration test
CN108592720A (en) * 2018-05-16 2018-09-28 厦门爆破工程公司 A kind of vibration damping hole distribution method reducing ground blasting vibration
CN110736400A (en) * 2019-10-22 2020-01-31 长江重庆航道工程局 underwater drilling blasting vibration velocity calculation method considering rock internal structure
CN111426243A (en) * 2020-01-20 2020-07-17 重庆中环建设有限公司 Blasting parameter determination method based on differential vibration synthesis of different blast holes in near zone of blasting source
CN111780634A (en) * 2020-07-20 2020-10-16 中国铁建投资集团有限公司 Method for reducing blasting vibration velocity
WO2022051890A1 (en) * 2020-09-08 2022-03-17 中煤科工集团淮北爆破技术研究院有限公司 Vibration reduction method based on blast vibration prediction technology
CN113588064A (en) * 2021-07-06 2021-11-02 中铁二十局集团第六工程有限公司 Method for measuring blasting vibration data
CN113343503A (en) * 2021-07-26 2021-09-03 内蒙聚力工程爆破有限公司 Blasting method and device for predicting blasting vibration damage
CN114781040A (en) * 2022-05-10 2022-07-22 青岛理工大学 Building blasting vibration safety evaluation method, device, equipment and medium

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
小净距双洞隧道下穿建筑物爆破振速控制技术研究;李立功等;隧道建设;36(05);第592-599页 *

Also Published As

Publication number Publication date
CN116242209A (en) 2023-06-09

Similar Documents

Publication Publication Date Title
CN107884822B (en) Method for improving positioning precision of mining micro-seismic source
CN102829899B (en) Quick surrounding rock stress measuring method applicable to buried circular tunnels
CN107478523B (en) Method and system for testing blasting vibration speed of middle rock wall of small-spacing tunnel
CN115600368A (en) Deep hole blasting optimization method, device and equipment and readable storage medium
CN116242209B (en) Vibration damping rate calculation method, system and equipment for vibration damping hole and readable storage medium
JP6998014B2 (en) Blasting method
JP6420054B2 (en) Elastic wave velocity measurement method
CN115146487B (en) Deep-buried tunnel blasting parameter evaluation method
CN103389489A (en) Micro earthquake monitoring and positioning method based on highly-deviated well
CN115371791A (en) Underground pipeline vibration velocity measuring method and system
CN111751872A (en) System and method for predicting blasting vibration speed of central area of side wall of underground cavern
CN107060732A (en) The accurate measurement and positioning method of multiple spot three-dimensional coordinate in curved drilling hole
Ragam et al. Selection and evolution of MEMS accelerometer sensor for measurement of blast-induced peak particle velocity
CN104930937B (en) Circular divided blasting tunnel construction method capable of effectively controlling blasting vibration velocity
CN111784025B (en) System and method for predicting blasting vibration speed of side wall of vertical shaft based on simply supported beam model
CN115688428B (en) Method, device and equipment for calculating explosive quantity critical value of blasting near zone supporting structure
CN114723863B (en) Vibration speed cloud picture prediction method and system for open bench blasting and electronic equipment
CN107036495A (en) Hyperboloid of one sheet reaming chock blasting method from the bottom of one's heart
CN114936342A (en) Method, device, equipment and medium for calculating close-range blasting explosive quantity of protective object
CN116307045B (en) Method, system, equipment and medium for calculating vibration speed of building under tunnel blasting
CN116258285B (en) Porous small-clear-distance tunnel blasting vibration speed prediction method, device, equipment and medium
CN107576319A (en) Without pipe omnidirectional self-correction real-time displacement measurement pipe
CN107023301A (en) Model test apparatus and method for studying Tunnel Engineering blast simulated acid rain
CN106404161A (en) Calculation method for precisely positioning vibration position in tunnel construction
CN105841715A (en) High-dynamic dual-axis angular-rate gyroscope, and error compensation for zero offset and scale factor

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