CN-122016278-A - Downward filling mining method filling body lateral cantilever safety monitoring method based on dynamic response

CN122016278ACN 122016278 ACN122016278 ACN 122016278ACN-122016278-A

Abstract

The invention discloses a dynamic response-based downward filling mining method filling body lateral cantilever safety monitoring method, which relates to the technical field of mine safety engineering and structural health monitoring, and is characterized in that vibration data are collected through a three-way acceleration sensor by taking regular micro-vibration in a roadway environment as a natural excitation source, steady-state vibration segments are screened through standardization, covariance matrix feature decomposition and minimum projection criteria, main vibration azimuth sequences are separated, natural frequencies are extracted, a quantitative mapping model based on the natural frequencies and the safety coefficients is combined with the current natural frequencies to calculate corresponding safety coefficients, and stability of the cantilever is evaluated and support decision is assisted based on the safety coefficients. The invention realizes non-invasive continuous on-line monitoring, provides quantitative basis for supporting decision-making, avoids resource waste and potential safety hazard, and assists intelligent safety control of mines.

Inventors

LI JICHUAN
ZHANG XIAOYONG
WANG CHUANGYE
WU XIANGYE
YANG ZHANFENG
SHI MING
WANG JIAN
YIN HUIFANG

Assignees

内蒙古科技大学

Dates

Publication Date: 20260512
Application Date: 20260203

Claims (8)

1. A downward filling mining method filling body lateral cantilever safety monitoring method based on dynamic response, which is characterized by comprising the following steps: S1, arranging a three-way vibration acceleration sensor on the surface of a cantilever structure of a roadway lateral filling body, collecting three-way vibration acceleration time sequence data, dividing the three-way vibration acceleration time sequence data into a plurality of time periods, and constructing a data matrix of the three-way vibration acceleration time sequence data of each time period; S2, carrying out standardization processing on the data matrix to obtain a covariance matrix, carrying out matrix eigenvalue decomposition on the covariance matrix to obtain eigenvalues and corresponding eigenvectors, and constructing an eigenvalue matrix based on the eigenvectors; s3, screening column vectors in the eigenvalue matrix of each time period based on a minimum projection criterion to determine a steady-state vibration time period index; s4, repeating the step S2, and extracting a feature vector corresponding to the maximum feature value as a first main direction feature vector, and extracting a feature vector corresponding to the minimum feature value as a second main direction feature vector based on the reconstructed triaxial vibration acceleration time sequence data; s5, respectively projecting the reconstructed triaxial vibration acceleration data onto the first main direction feature vector and the second main direction feature vector to obtain vibration time sequence data of the main vibration direction of the cantilever body and reference vibration time sequence data of surrounding rock; s6, respectively calculating frequency spectrums of the vibration time sequence data of the main vibration direction of the cantilever body and the reference vibration time sequence data of the surrounding rock, and calculating a amplitude-frequency response function based on the frequency spectrums of the vibration time sequence data and the reference vibration time sequence data of the surrounding rock; And S7, identifying the frequency corresponding to the amplitude-frequency response function peak value, taking the frequency as the current natural frequency of the cantilever body, calculating the corresponding safety coefficient by combining the current natural frequency based on a mapping relation function of the pre-calibrated natural frequency of the cantilever body and the safety coefficient, and evaluating the stability of the cantilever body and assisting in supporting decision based on the safety coefficient.
2. The method for monitoring the lateral cantilever safety of the filling body by the downward filling mining method based on the dynamic response is characterized in that a three-way vibration acceleration sensor is arranged on the surface of a cantilever structure formed by the lateral filling body of a roadway in S1, normal micro-vibration existing in a mine roadway environment is used as an excitation source, three-axis vibration acceleration time sequence data of the cantilever structure are collected in real time, the three-way vibration acceleration time sequence data are divided into a plurality of time periods, and a data matrix of the three-way vibration acceleration time sequence data of each time period is constructed.
3. The method for monitoring the safety of a lateral cantilever of a filling body by a downward filling mining method based on a dynamic response according to claim 1, wherein the construction of the eigenvalue matrix based on the eigenvector comprises equally dividing three-way vibration acceleration time series data into Segment of the first pair Repeating step S2 to obtain feature vector Constructing a characteristic value matrix Wherein, the method comprises the steps of, Characteristic value vectors of the ith three-way vibration sequence are respectively as follows 、 And (3) with 。
4. A method for monitoring the lateral cantilever safety of a filling body by a downward filling mining method based on a dynamic response according to claim 3, wherein the step S3 specifically comprises: S3.1, calculate Sum of absolute values of correlation coefficients of each column vector and other column vectors: ; In the formula, Calculating operators for the cross-correlation coefficients; Is a matrix Is the first of (2) A column vector; Is a kronecker function; s3.2, initializing an orthogonal complement space operator, and setting To correspond to Is recorded as a column vector of (1) The column space formed by the sheet is recorded as Initialized orthogonal complement space operator: Wherein, the method comprises the steps of, A space projection operator is orthogonally complemented; Is a unit matrix; S3.3 pair Projection norm calculation of residual candidate column vector ; S3.4, select Updating : ; S3.5 repeating S3.3-S3.4 until Sharing in common Column vector and construct index set ; S3.6, based on index set Selecting three-dimensional vibration acceleration sequence segments, and performing end-to-end reconstruction to obtain reconstructed three-dimensional vibration acceleration time sequence data.
5. The method for monitoring the lateral cantilever safety of a filling body by a downward filling mining method based on power response according to claim 4, wherein S4 specifically comprises reconstructing a matrix based on reconstructed time series data Extracting the main vibration azimuth vector of the cantilever body through S2 And lateral filling surrounding rock reference vibration vector Wherein, the method comprises the steps of, The feature vector corresponding to the maximum feature value, And the feature vector corresponding to the minimum feature value.
6. The method for monitoring the safety of a lateral cantilever of a filling body by a downward filling mining method based on power response according to claim 5, wherein S5 comprises the following steps of calculating a data sequence of vibration time sequence data of a main vibration direction of the cantilever body Data sequence of surrounding rock reference vibration time sequence data 。
7. The method for monitoring the safety of a lateral cantilever of a filling body by a downward filling mining method based on power response according to claim 5, wherein S6 specifically comprises a data sequence of vibration time series data of the main vibration direction of the cantilever body Data sequence of surrounding rock reference vibration time sequence data Respectively calculating frequency spectrums, and calculating amplitude-frequency response functions based on the frequency spectrum data Wherein, the method comprises the steps of, And (3) with The main vibration azimuth spectrum of the cantilever body and the surrounding rock reference spectrum are respectively obtained.
8. The method for monitoring the safety of the lateral cantilever body of the filling body of the downward filling mining method based on the dynamic response according to claim 1, wherein the mapping relation function based on the pre-calibrated natural frequency of the cantilever body and the safety coefficient specifically comprises the following steps: ; Wherein, the For the natural frequency of the cantilever body calibrated in advance, As a safety factor, the safety factor of the device, Is a mapping relation function.

Description

Downward filling mining method filling body lateral cantilever safety monitoring method based on dynamic response Technical Field The invention relates to the technical field of mine safety engineering and structural health monitoring, in particular to a downward filling mining method filling body lateral cantilever safety monitoring method based on power response. Background The downward route filling mining method has become one of core technologies for safely and efficiently mining deep broken ore bodies by virtue of remarkable advantages of the downward route filling mining method in the aspects of controlling the ground pressure, improving the recovery rate of ore pillars and adapting to complex ore body conditions. In the process of stoping the approach, the filling bodies serving as surrounding rocks at two sides are separated from the lower-layer ore body, and the irregular convex parts of the local filling bodies form a cantilever structure. One end of the structure is restrained by the surrounding rock of the lateral filling body, and the other end is in a free state. The suspension section of the filling body and the joint of the matrix are easy to generate micro-cracks under the influence of blasting vibration, and the cantilever structure is separated along the crack surface along with the continuous expansion of the dominant crack, so that the crack-falling type instability is finally initiated under the action of dead weight. Because the treatment of the raised cantilever structure may interrupt the recovery process, collide with the neighboring approach blasting plan, and the forced cutting of the raised cantilever structure may cause new cracks in the filling body, the raised cantilever body should be preferably treated by support. However, excessive support tends to cause resource cost waste, and how to realize accurate decision of support time and strength on the premise of ensuring safety becomes a key contradiction of field engineering optimization. Therefore, a set of quantitative dynamic risk sensing mechanism needs to be established, and the critical point of damage evolution is identified through an easily-realized monitoring means, so that resource redundancy caused by empirical support is avoided, and meanwhile, a safety early warning blind area caused by monitoring lag is avoided. Disclosure of Invention In view of the above, the invention provides a downward filling mining method filling body lateral cantilever safety monitoring method based on dynamic response, which is used for solving the problems in the background technology. In order to achieve the above purpose, the present invention adopts the following technical scheme: A downward filling mining method filling body lateral cantilever safety monitoring method based on dynamic response, comprising: S1, arranging a three-way vibration acceleration sensor on the surface of a cantilever structure of a roadway lateral filling body, collecting three-way vibration acceleration time sequence data, dividing the three-way vibration acceleration time sequence data into a plurality of time periods, and constructing a data matrix of the three-way vibration acceleration time sequence data of each time period; S2, carrying out standardization processing on the data matrix to obtain a covariance matrix, carrying out matrix eigenvalue decomposition on the covariance matrix to obtain eigenvalues and corresponding eigenvectors, and constructing an eigenvalue matrix based on the eigenvectors; s3, screening column vectors in the eigenvalue matrix of each time period based on a minimum projection criterion to determine a steady-state vibration time period index; s4, repeating the step S2, and extracting a feature vector corresponding to the maximum feature value as a first main direction feature vector, and extracting a feature vector corresponding to the minimum feature value as a second main direction feature vector based on the reconstructed triaxial vibration acceleration time sequence data; s5, respectively projecting the reconstructed triaxial vibration acceleration data onto the first main direction feature vector and the second main direction feature vector to obtain vibration time sequence data of the main vibration direction of the cantilever body and reference vibration time sequence data of surrounding rock; s6, respectively calculating frequency spectrums of the vibration time sequence data of the main vibration direction of the cantilever body and the reference vibration time sequence data of the surrounding rock, and calculating a amplitude-frequency response function based on the frequency spectrums of the vibration time sequence data and the reference vibration time sequence data of the surrounding rock; And S7, identifying the frequency corresponding to the amplitude-frequency response function peak value, taking the frequency as the current natural frequency of the cantilever body, calculating the corresponding safety coeffi