Blasting vibration signal coordinate correction method, electronic equipment and storage medium

Abstract

The invention belongs to the technical field of engineering blasting, and relates to a blasting vibration signal coordinate correction method, electronic equipment and a storage medium, and the method comprises the following steps: obtaining a three-component original blasting vibration monitoring signal at a measuring point during blasting of a blasting source; obtaining a relative position relation between the explosion source and the measuring point, and determining a P-wave real polarization angle; determining a P-wave polarization angle based on the original blasting vibration monitoring signal; determining an error angle from a monitoring data coordinate system to a real data coordinate system based on the P-wave real polarization angle and the P-wave polarization angle, and expanding the error angle into a three-dimensional coordinate transformation matrix; and correcting the three-component original blasting vibration monitoring signal based on the three-dimensional coordinate transformation matrix to obtain a real blasting vibration signal. According to the method, the azimuth error existing in actual monitoring can be compensated, so that the fine analysis of the blasting vibration signal is realized, the analysis error is reduced, the precision requirement of the blasting signal analysis is met, and the reliability of the blasting vibration signal analysis is remarkably improved.

Description

TECHNICAL FIELD

[0001] The present application belongs to the technical field of engineering blasting, and particularly relates to a blasting vibration signal coordinate correction method, an electronic device and a storage medium. BACKGROUND

[0002] When explosives explode, the energy released thereby is partly used for rock breaking, and the rest is dissipated in various forms, one of which is in the form of wave propagation to form blasting seismic waves. Blasting seismic waves are composed of wave components with different polarization characteristics, and can be generally divided into two categories, namely body waves (P waves, S waves) and surface waves (R waves, L waves). According to the difference in propagation speed of different waves, the fastest wave is the compression wave (P wave), followed by the shear wave (S wave) and the Rayleigh surface wave (R wave). Therefore, the P wave, as the first wave to reach the monitoring device, has clear polarization characteristics and is not disturbed by other complex wave components, and can be better utilized.

[0003] In the process of engineering blasting construction, blasting vibration, as the first negative effect, has been widely concerned by researchers, and the evaluation and hazard control of blasting vibration have been a research hotspot. Blasting vibration monitoring, as one of the most direct and effective methods for studying blasting vibration, can provide reliable data support for in-depth study of blasting vibration.

[0004] When monitoring and collecting blasting vibration signals, the monitoring device is often arranged in a safe area away from the blast area within a certain range, because the blast area within a certain range can seriously threaten the safety of personnel and equipment. The monitoring device is easily affected by the test environment, human factors and the instrument system itself when collecting signals. If the orientation of the sensor is not aligned with the true blasting vibration orientation, the three-dimensional coordinates of the collected blasting vibration signals will not be the true blasting vibration orientation, which will cause a certain deviation and affect the subsequent analysis of the vibration signals.

[0005] However, the above-mentioned influencing factors are not considered in the current analysis of blasting vibration signals, so there is an error between the research data and the actual situation, which causes the analysis result to be inaccurate and cannot meet the precision requirements of blasting signal analysis. SUMMARY

[0006] Therefore, the present application provides a blasting vibration signal coordinate correction method, an electronic device and a storage medium to solve the problems in the prior art.

[0007] The technical scheme of the present application is as follows: A blasting vibration signal coordinate correction method, comprising the following steps: acquiring three-component original blasting vibration monitoring signals at the measuring point when the explosive source explodes; Obtaining the relative position relationship between the explosion source and the measuring point, and determining the true polarization angle of the P wave; Determining the polarization angle of the P wave based on the original blasting vibration monitoring signal; Determining the error angle from the monitoring data coordinate system to the true data coordinate system based on the true polarization angle of the P wave and the polarization angle of the P wave, and expanding the error angle into a three-dimensional coordinate transformation matrix; Correcting the original three-component blasting vibration monitoring signal based on the three-dimensional coordinate transformation matrix to obtain the true blasting vibration signal.

[0008] Further, the true polarization angle of the P wave is determined based on the following formula: , wherein, θ is the true polarization angle of the P wave, H is the depth of the blast hole where the explosion source is located, d is the distance between the blast hole where the explosion source is located and the measuring point.

[0009] Further, the polarization angle of the P wave is determined based on the original blasting vibration monitoring signal, including the following steps: Numerically integrating the original blasting vibration monitoring signal to obtain the blasting vibration displacement time history curve; Drawing the particle motion trajectory vector diagram of the XOY horizontal plane and the particle motion trajectory vector diagram of the ROZ vertical plane based on the blasting vibration displacement time history curve; Determining the polarization angle of the P wave based on the particle motion trajectory vector diagram of the XOY horizontal plane and the particle motion trajectory vector diagram of the ROZ vertical plane.

[0010] Further, the numerical integration of the original blasting vibration monitoring signal is performed to obtain the blasting vibration displacement time history curve, including the following steps: Vector synthesis of the original blasting vibration monitoring signals in the X and Y directions to obtain the blasting vibration signal in the R direction; Numerically integrating the blasting vibration signal in the R direction to obtain the blasting vibration displacement time history curve.

[0011] Further, the particle motion trajectory vector diagram of the XOY horizontal plane and the particle motion trajectory vector diagram of the ROZ vertical plane are drawn based on the blasting vibration displacement time history curve, including the following steps: In the horizontal XOY rectangular coordinate system, the X and Y direction displacement coordinates of the particle are sequentially drawn in time sequence based on the blasting vibration displacement time history curve to obtain the particle motion trajectory vector diagram of the XOY horizontal plane; In the vertical ROZ rectangular coordinate system, the R and Z direction displacement coordinates of the particle are sequentially drawn in time sequence based on the blasting vibration displacement time history curve to obtain the particle motion trajectory vector diagram of the ROZ vertical plane.

[0012] Further, the P-wave polarization angle is determined based on the particle motion trajectory vector diagram in XOY horizontal plane and the particle motion trajectory vector diagram in ROZ vertical plane, comprising the following steps: the polarization angle of the P-wave in the particle motion trajectory vector diagram in XOY horizontal plane is obtained as a first P-wave polarization angle α ; the polarization angle of the P-wave in the particle motion trajectory vector diagram in ROZ vertical plane is obtained as a second P-wave polarization angle θ ’.

[0013] Further, the error angle of the monitoring data coordinate system to the real data coordinate system is determined, and the error angle is expanded to a three-dimensional coordinate transformation matrix, comprising the following steps: in XOY horizontal plane, the error angle of the monitoring data coordinate system to the real data coordinate system is equal to the first P-wave polarization angle α； in ROZ vertical plane, the error angle of the monitoring data coordinate system to the real data coordinate system is obtained by the following formula: ; wherein, is a vertical error angle, θ is a real P-wave polarization angle, is the second P-wave polarization angle; the three-dimensional coordinate transformation matrix is determined by the following formula: wherein, T is a three-dimensional coordinate transformation matrix, is a vertical error angle, α is a first P-wave polarization angle, is a transformation matrix of rotation around Y axis, is a transformation matrix of rotation around Z axis.

[0014] Further, the real blasting vibration signal is obtained by the following formula: wherein, , , V’ ( t )=[ X’ ( t ), Y’ ( t ), Z’ ( t ] T , T is a three-dimensional coordinate transformation matrix, is a vertical error angle, αis the first P-wave polarization angle, V t ) is the original blasting vibration monitoring signal, X t Y t Z t V’ t X’ t Y’ t Z’ t

[0015] The electronic device comprises a storage, a processor and a computer program stored in the storage and capable of running on the processor, and the processor executes the computer program to realize the blasting vibration signal coordinate correction method.

[0016] The computer readable storage medium has a computer program stored thereon, and the computer program is executed by a processor to realize the blasting vibration signal coordinate correction method.

[0017] Compared with the prior art, the blasting vibration signal coordinate correction method, the electronic device and the storage medium provided by the present application can compensate for the azimuth error existing in actual monitoring, realize fine analysis of the blasting vibration signal, reduce analysis error, meet the precision requirement of blasting signal analysis, and have the following advantages: (1) The present application utilizes the characteristics that the P-wave velocity in the blasting seismic wave is the fastest, the first wave usually reaches the monitoring instrument and is not affected by other waves, discriminates the polarization direction by using the polarization analysis method, accurately determines the real blasting vibration data coordinate azimuth by using matrix coordinate transformation, realizes accurate correction of the X, Y and Z directions, significantly reduces the data analysis deviation caused by the azimuth error, and improves the scientificity and reliability of the blasting vibration signal analysis.

[0018] (2) The traditional blasting site data monitoring usually needs workers to accurately install and deploy the monitoring instrument, but due to the complex actual construction site environment, it is difficult to ensure the accurate alignment of the monitoring instrument. The present application does not need to accurately deploy the instrument by workers on site, and even if the equipment has obvious installation inclination or azimuth deviation, the accurate alignment of the data can still be realized through post-correction, which improves the universality and easy operability of the method in complex site environment.

[0019] ​​​​​​​​​​​​​​​(3) If the azimuth error of the monitoring data cannot be accurately corrected, it is easy to cause deviation of the blasting vibration risk evaluation, which may have an adverse effect on safety control. The implementation of the present application can effectively ensure the accuracy of the identification and separation of the components of the blasting seismic wave in subsequent analysis, and further effectively analyze the component characteristics and propagation law of the blasting seismic wave, so as to make the blasting vibration evaluation more accurate, further improve the prediction accuracy of the blasting vibration hazard and the scientificity of the safety control, improve the engineering safety level, has strong practicability, and is worth promoting. BRIEF DESCRIPTION OF DRAWINGS

[0020] Figure 1 A flow chart of a blasting vibration signal coordinate correction method in an embodiment of the present application.

[0021] Figure 2 A schematic diagram of the arrangement of the blast hole and the measuring point in which the blast source is located in an embodiment of the present application.

[0022] Figure 3 A schematic diagram of the arrangement of the XOY horizontal plane and the ROZ vertical plane of the blast hole and the measuring point in which the blast source is located in an embodiment of the present application.

[0023] Figure 4 X-direction, Y-direction and Z-direction original blasting vibration monitoring signals obtained by a vibration sensor in an embodiment of the present application.

[0024] Figure 5 X-direction, Y-direction and Z-direction blasting vibration displacement time history curves in an embodiment of the present application.

[0025] Figure 6 Original blasting vibration monitoring signals and blasting vibration displacement time history curves in the horizontal radial (R-direction) in an embodiment of the present application.

[0026] Figure 7 Particle motion trajectory vector diagrams of the XOY horizontal plane and the ROZ vertical plane in an embodiment of the present application. DETAILED DESCRIPTION

[0027] The present application provides a blasting vibration signal coordinate correction method, an electronic device and a storage medium to solve the above problems. In order to enable those skilled in the art to better understand the technical solutions of the present application and implement them, the technical solutions in the present application will be described in detail below with reference to the drawings.

[0028] In the description of the present application, it needs to be understood that the orientation or positional relationship indicated by the terms "center", "longitudinal", "transverse", "length", "width", "thickness", "upper", "lower", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", "axial", "radial", "circumferential", and the like is based on the orientation or positional relationship shown in the drawings, and is only for the purpose of facilitating the description of the present application and simplifying the description, and does not indicate or imply that the device or element referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore cannot be understood as a limitation of the present application.

[0029] In addition, it needs to be further pointed out that, in the description of the embodiments of the present application, unless otherwise specified, " / " means or, for example, A / B can mean A or B: "and / or" in the text only describes the relationship between the associated objects, which means that there can be three relationships, for example, A and / or B, which can mean that A exists alone, A and B exist together, and B exists alone, and in addition, in the description of the embodiments of the present application, "multiple" means two or more than two.

[0030] The following terms "first", "second", "third", "fourth" are only for descriptive purposes, and cannot be understood as implying or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined with "first", "second", "third", "fourth" can explicitly or implicitly include one or more of the features, and in the description of the present application, unless otherwise specified, the meaning of "multiple" is at least two, such as two, three, etc., unless otherwise explicitly specified.

[0031] In the present application, unless otherwise explicitly specified and limited, the terms "mounting", "connection", "connection", "fixing" and the like should be understood in a broad sense, for example, it can be fixedly connected, or it can be detachably connected, or it can be integrated; it can be mechanically connected, or it can be electrically connected; it can be directly connected, or it can be indirectly connected through an intermediate medium; it can be the internal communication of two elements or the interaction relationship between two elements, unless otherwise explicitly limited. For those skilled in the art, the specific meaning of the above terms in the present application can be understood according to the specific circumstances.

[0032] In the present application, unless otherwise explicitly specified and limited, the first feature is "on" or "under" the second feature, which can be direct contact between the first and second features, or indirect contact between the first and second features through an intermediate medium. Moreover, the first feature "above", "above" and "above" the second feature can be directly above or obliquely above the first feature, or only indicate that the horizontal height of the first feature is higher than that of the second feature. The first feature "below", "below" and "below" the second feature can be directly below or obliquely below the first feature, or only indicate that the horizontal height of the first feature is less than that of the second feature.

[0033] It is to be understood that when an element such as a layer, region or substrate is referred to as being "on" or "connected to" another element, it can be directly on or connected to the other element or intervening elements can be present. In addition, it is to be understood that the term "connected" refers to physical or logical connections by itself or in conjunction with other elements. Also, it is to be understood that the term "vertical" "horizontal" "upper" "lower" "left" "right" and similar terms are used only for illustrative purposes and do not indicate the only orientation.

[0034] Embodiment 1 The present application provides a blasting vibration signal coordinate correction method, as shown in the method comprises the following steps: Figure 1 The method comprises the following steps: obtaining three-component original blasting vibration monitoring signals of the blasting source at the measuring point; obtaining the relative position relationship between the blasting source and the measuring point, and determining the true polarization angle of P wave; determining the polarization angle of P wave based on the original blasting vibration monitoring signals; determining the error angle from the monitoring data coordinate system to the true data coordinate system based on the true polarization angle of P wave and the polarization angle of P wave, and expanding the error angle into a three-dimensional coordinate transformation matrix; correcting the three-component original blasting vibration monitoring signals based on the three-dimensional coordinate transformation matrix to obtain the true blasting vibration signals.

[0035] Specifically, as a further refinement of the above method, the true polarization angle of P wave is determined based on the following formula: , wherein, θ is the true polarization angle of P wave, H is the depth of the blast hole where the blasting source is located, d is the distance between the blast hole where the blasting source is located and the measuring point.

[0036] Specifically, as a further refinement of the above method, the polarization angle of P wave is determined based on the original blasting vibration monitoring signals, comprising the following steps: performing numerical integral processing on the original blasting vibration monitoring signals to obtain a blasting vibration displacement time history curve; drawing a particle motion trajectory vector diagram of XOY horizontal plane and a particle motion trajectory vector diagram of ROZ vertical plane based on the blasting vibration displacement time history curve; determining the polarization angle of P wave based on the particle motion trajectory vector diagram of XOY horizontal plane and the particle motion trajectory vector diagram of ROZ vertical plane.

[0037] Specifically, as a further refinement of the above method, the numerical integral processing is performed on the original blasting vibration monitoring signals to obtain a blasting vibration displacement time history curve, comprising the following steps: vector synthesis of the original blasting vibration monitoring signals in the X direction and the Y direction to obtain a blasting vibration signal in the R direction; numerical integration of the blasting vibration signal in the R direction to obtain a blasting vibration displacement time-history curve.

[0038] Specifically, as a further refinement of the above method, based on the blasting vibration displacement time-history curve, a particle motion trajectory vector diagram in the XOY horizontal plane and a particle motion trajectory vector diagram in the ROZ vertical plane are drawn, including the following steps: In the horizontal XOY rectangular coordinate system, based on the blasting vibration displacement time-history curve, the X direction and the Y direction displacement coordinates of the particle are sequentially drawn in time sequence to obtain the particle motion trajectory vector diagram in the XOY horizontal plane; In the vertical ROZ rectangular coordinate system, based on the blasting vibration displacement time-history curve, the R direction and the Z direction displacement coordinates of the particle are sequentially drawn in time sequence to obtain the particle motion trajectory vector diagram in the ROZ vertical plane.

[0039] Specifically, as a further refinement of the above method, based on the particle motion trajectory vector diagram in the XOY horizontal plane and the particle motion trajectory vector diagram in the ROZ vertical plane, the P wave polarization angle is determined, including the following steps: obtain the polarization angle of the P wave in the particle motion trajectory vector diagram in the XOY horizontal plane as a first P wave polarization angle α ; obtain the polarization angle of the P wave in the particle motion trajectory vector diagram in the ROZ vertical plane as a second P wave polarization angle θ ’.

[0040] Specifically, as a further refinement of the above method, the error angle of the monitoring data coordinate system to the real data coordinate system is determined, and the error angle is expanded into a three-dimensional coordinate transformation matrix, including the following steps: In the XOY horizontal plane, the error angle of the monitoring data coordinate system to the real data coordinate system is equal to the first P wave polarization angle α； In the ROZ vertical plane, the error angle of the monitoring data coordinate system to the real data coordinate system is obtained by the following formula: ; wherein, is the vertical error angle, θ is the real polarization angle of the P wave, is the second P wave polarization angle; determine the three-dimensional coordinate transformation matrix according to the following formula: wherein, T is the three-dimensional coordinate transformation matrix, is the vertical error angle,α is a first P-wave polarization angle, is a transformation matrix of rotation around the Y axis, is a transformation matrix of rotation around the Z axis.

[0041] Specifically, as a further refinement of the above method, the real blasting vibration signal is obtained using the following formula: wherein, , , V’ t ) = [ X’ t ], Y’ t , Z’ t ] T , T is a three-dimensional coordinate transformation matrix, is a vertical error angle, α is a first P-wave polarization angle, V t is an original blasting vibration monitoring signal, X t , Y t and Z t represent the original velocity time history curves in the X direction, Y direction and Z direction, respectively, V’ t is a real blasting vibration signal, X’ t , Y’ t and Z’ t represent the real velocity time history curves in the X direction, Y direction and Z direction, respectively.

[0042] As a further description of the above method, with reference to the accompanying drawings, taking a single-hole blasting experiment of a certain petrochemical base as an example, the blasting-induced vibration is monitored by a vibration sensor, and the blasting vibration data azimuth correction is carried out based on the P-wave first arrival polarization direction. The refined landing implementation scheme of the blasting vibration signal coordinate correction method provided in this embodiment includes the following steps: Step S1: obtaining three-component original blasting vibration monitoring signals at the measuring point during blasting of the explosion source.

[0043] The blasting vibration monitoring instrument is set at the measuring point, and the X direction, Y direction and Z direction original blasting vibration monitoring signals generated during blasting of the explosion source are obtained by the blasting vibration monitoring instrument​​​​​​​​​​​​​ , see Figure 4 .

[0044] Step S2: Obtain the relative position relationship between the explosion source and the measuring point, and determine the true polarization angle of P wave.

[0045] The layout in actual application is shown in Figure 2 , the relevant data is obtained by measurement, and it is known that the depth of the blast hole where the explosion source is located is H = 6.0 m, and the horizontal distance between the blast hole where the explosion source is located and the measuring point is d = 10.0 m.

[0046] The true polarization angle of P wave is determined based on the following formula: , wherein θ is the true polarization angle of P wave, H is the depth of the blast hole where the explosion source is located, d is the distance between the blast hole where the explosion source is located and the measuring point.

[0047] Substitute the relevant data, .

[0048] Step S3: Numerical integral processing is performed on the original blasting vibration monitoring signal to obtain the blasting vibration displacement time history curve.

[0049] For the original blasting vibration monitoring signals in X direction, Y direction and Z direction, the blasting vibration displacement time history curves in X direction, Y direction and Z direction as shown in Figure 5 can be obtained through numerical integral calculation.

[0050] The original blasting vibration monitoring signals in X direction and Y direction are vector synthesized to obtain the blasting vibration signal in horizontal radial direction, i.e. R direction, and the blasting vibration displacement time history curve is obtained through numerical integral calculation, see the attached Figure 6 .

[0051] Step S4: Draw the particle motion trajectory vector diagram of XOY horizontal plane and ROZ vertical plane based on the blasting vibration displacement time history curve.

[0052] In the horizontal XOY rectangular coordinate system, the X direction and Y direction displacement coordinates of the particle are sequentially drawn in time sequence based on the blasting vibration displacement time history curve, and the particle motion trajectory vector diagram of XOY horizontal plane is obtained.

[0053] In the vertical ROZ rectangular coordinate system, the R direction and Z direction displacement coordinates of the particle are sequentially drawn in time sequence based on the blasting vibration displacement time history curve, and the particle motion trajectory vector diagram of ROZ vertical plane is obtained, see the attached Figure 7 .

[0054] Step S5: Determine the real polarization direction of P wave and the error angle between the monitoring data coordinate system and the real data coordinate system.

[0055] First, measure the polarization angle of P wave in particle motion trajectory vector diagram, see attached Figure 7 .

[0056] In XOY horizontal plane, the polarization angle of P wave is the first P wave polarization angle α = 5°.

[0057] In ROZ vertical plane, the polarization angle of P wave is the second P wave polarization angle θ ’= 23°, then determine the real polarization direction of P wave, see attached Figure 3 .

[0058] In XOY horizontal plane, the real polarization direction of P wave is the direction of the connecting line between the explosion source and the measuring point.

[0059] In ROZ vertical plane, the real polarization direction of P wave can be calculated according to the relative position relationship between the explosion source and the measuring point as follows: .

[0060] Finally, determine the error angle between the monitoring data coordinate system and the real data coordinate system: In XOY horizontal plane, the error angle between the monitoring data coordinate system (X’, Y’) and the real data coordinate system (X, Y) is equal to the first P wave polarization angle α，α = 5°.

[0061] In ROZ vertical plane, the error angle between the monitoring data coordinate system (R’, Z’) and the real data coordinate system (R, Z) is obtained by using the following formula: ; Wherein, is the vertical error angle, θ is the real polarization angle of P wave, is the second P wave polarization angle.

[0062] Substitute the relevant data, .

[0063] Step S6: Calculate the three-dimensional coordinate transformation matrix.

[0064] Define the coordinate transformation matrix of rotation around Z axis R Z ( α= 5), the monitoring data coordinate system (X, Y, Z) can be transformed into the horizontal real data coordinate system (X’, Y’, Z) after transformation, and then define the coordinate transformation matrix of rotation around Y axis R Y ( β= 6.3), the transformed horizontal real data coordinate system (X', Y', Z) can be transformed into the real data coordinate system (X', Y', Z').

[0065] Coordinate transformation matrix R Z ( α= 5) as follows: , Coordinate transformation matrix R Y ( β = 6.3) as follows: , Three-dimensional coordinate transformation matrix T as follows: Step S7: correcting the original blasting vibration monitoring signal by using the three-dimensional coordinate transformation matrix to obtain the real blasting vibration signal.

[0066] Correcting the original blasting vibration data T ( V t ) by using the three-dimensional coordinate transformation matrix, the real blasting vibration signal can be calculated as follows: that is: , that is: .

[0067] The method can also be implemented by relying on an electronic device, and the structure of the electronic device specifically includes a storage, a processor, and a computer program stored on the storage and executable on the processor. The processor executes the computer program to implement the blasting vibration signal coordinate correction method, the electronic device, and the storage medium as described above.

[0068] A communication interface is arranged between the storage and the processor, and the communication interface realizes signal connection and data transmission between the two. The communication interface can be one of a serial interface or a parallel interface.

[0069] The processor can be a central processing unit (CPU), or a specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application.

[0070] The method can also be implemented by relying on a computer readable storage medium, and the computer readable storage medium stores a computer program. The computer program is executed by the processor to implement the blasting vibration signal coordinate correction method, the electronic device, and the storage medium as described above.

[0071] ​The blasting vibration signal coordinate correction method, the electronic equipment and the storage medium provided by the application can accurately correct the data coordinate deviation of the blasting vibration monitoring signal caused by the deployment deviation of the monitoring instrument, significantly reduce the data analysis deviation caused by the azimuth error, improve the scientificity and reliability of the blasting vibration signal analysis, and the universality and easy operability in the complex field environment, can ensure the accuracy of the identification and separation of the blasting seismic wave components in the subsequent analysis, and then effectively analyze the component characteristics and propagation law of the blasting seismic wave, so that the blasting vibration evaluation is more accurate, and the practicability is strong, and it is worth promoting.

[0072] The above disclosure is only the preferred embodiment of the application, but the embodiments of the application are not limited thereto, and any changes that can be thought of by those skilled in the art should fall within the protection scope of the application.

Claims

1. A method of correcting a blasting vibration signal coordinate, characterized by, The method comprises the following steps: obtaining three-component original blasting vibration monitoring signals of the blasting source at the measuring point during blasting; obtaining the relative position relationship between the blasting source and the measuring point, and determining a true polarization angle of P wave; determining a P wave polarization angle based on the original blasting vibration monitoring signals; determining an error angle of a monitoring data coordinate system to a true data coordinate system based on the true polarization angle of P wave and the P wave polarization angle, and expanding the error angle into a three-dimensional coordinate transformation matrix; correcting the three-component original blasting vibration monitoring signals based on the three-dimensional coordinate transformation matrix to obtain true blasting vibration signals.

2. The method of blast vibration signal coordinate correction according to claim 1, characterized in that, The true polarization angle of P wave is determined based on the following formula: ， wherein, θ is the true polarization angle of P-wave, H is the depth of the blast hole where the explosion source is located, d is the distance between the blast hole where the explosion source is located and the measuring point.

3. The method of claim 1, wherein the method further comprises: The P wave polarization angle is determined based on the original blasting vibration monitoring signals, comprising the following steps: performing numerical integral processing on the original blasting vibration monitoring signals to obtain a blasting vibration displacement time-history curve; drawing a particle motion trajectory vector diagram of an XOY horizontal plane and a particle motion trajectory vector diagram of an ROZ vertical plane based on the blasting vibration displacement time-history curve; determining the P wave polarization angle based on the particle motion trajectory vector diagram of the XOY horizontal plane and the particle motion trajectory vector diagram of the ROZ vertical plane.

4. The method of blast vibration signal coordinate correction according to claim 3, wherein, The numerical integral processing on the original blasting vibration monitoring signals to obtain a blasting vibration displacement time-history curve comprises the following steps: performing vector synthesis on the original blasting vibration monitoring signals in the X direction and the Y direction to obtain a blasting vibration signal in the R direction; performing numerical integral on the blasting vibration signal in the R direction to obtain a blasting vibration displacement time-history curve.

5. The method of blast vibration signal coordinate correction according to claim 4, wherein, The drawing of the particle motion trajectory vector diagram of the XOY horizontal plane and the particle motion trajectory vector diagram of the ROZ vertical plane based on the blasting vibration displacement time-history curve comprises the following steps: in a horizontal XOY rectangular coordinate system, sequentially drawing X direction and Y direction displacement coordinates of a particle based on the blasting vibration displacement time-history curve in time sequence to obtain the particle motion trajectory vector diagram of the XOY horizontal plane; in a vertical ROZ rectangular coordinate system, sequentially drawing R direction and Z direction displacement coordinates of the particle based on the blasting vibration displacement time-history curve in time sequence to obtain the particle motion trajectory vector diagram of the ROZ vertical plane.

6. The method of blast vibration signal coordinate correction according to claim 5, wherein, The determination of the P wave polarization angle based on the particle motion trajectory vector diagram of the XOY horizontal plane and the particle motion trajectory vector diagram of the ROZ vertical plane comprises the following steps: Obtaining a polarization angle of a P wave in a particle motion trajectory vector diagram of an XOY horizontal plane as a first P wave polarization angle α ; acquiring a polarization angle of a P-wave in a particle motion trajectory vector diagram of a ROZ vertical plane as a second P-wave polarization angle θ ’ 7. The method of blast vibration signal coordinate correction according to claim 6, wherein, The determination of the error angle of the monitoring data coordinate system to the true data coordinate system and the expansion of the error angle into the three-dimensional coordinate transformation matrix comprise the following steps: In the XOY horizontal plane, the error angle of the monitoring data coordinate system to the real data coordinate system is equal to the first P-wave polarization angle α； in the ROZ vertical plane, the error angle of the monitoring data coordinate system to the true data coordinate system is obtained by using the following formula: ； wherein, is the vertical error angle, θ is the true polarization angle of the P-wave, is the second P-wave polarization angle; The three-dimensional coordinate transformation matrix is determined by using the following formula: wherein T is a three-dimensional coordinate transformation matrix, is a vertical error angle, α is a first P-wave polarization angle, is a transformation matrix for rotation around the Y-axis, is a transformation matrix for rotation around the Z-axis.

8. The method of blast vibration signal coordinate correction according to claim 7, wherein, The true blasting vibration signals are obtained by using the following formula: wherein , ， V’ ( t )=[ X’ ( t )， Y’ ( t )， Z’ ( t )] T ， T is a three-dimensional coordinate transformation matrix, is a vertical error angle, α is a first P-wave polarization angle, V t is an original blast vibration monitoring signal, X t , Y t and Z t represent original velocity time histories in X, Y and Z directions, respectively, V’ t is a true blast vibration signal, X’ t , Y’ t and Z’ t represent true velocity time histories in X, Y and Z directions, respectively.​​​​​​​​ 9. An electronic device, characterized by It comprises: a storage, a processor, and a computer program stored in the storage and capable of running on the processor, wherein the processor executes the computer program to implement the blasting vibration signal coordinate correction method according to claim 1.

10. A computer readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by the processor to implement the blasting vibration signal coordinate correction method according to claim 1.

Images