Mine intelligent management and control optimization method based on digital twinning

Abstract

The invention relates to the technical field of mine management and control, and particularly discloses an intelligent mine management and control optimization method based on digital twinning. By monitoring and analyzing the geological condition, the microseismic influence, the slope form and the support instability of the monitoring area where the target mine slope is located corresponding to the current monitoring time period and calculating the deformation degree evaluation coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring time period, the deformation condition of the target mine slope is visually displayed. The problems of complexity, time consumption, labor consumption and the like of periodic patrol and measurement of mine management personnel are effectively solved, subjective errors of mine slope deformation analysis are reduced to a certain extent, the accuracy, reliability and rationality of a mine slope deformation analysis result are ensured, and meanwhile, the mining safety of the mine personnel is effectively ensured.

Description

Technical Field

[0001] The present invention relates to the technical field of mine management and control, and specifically to a method for optimizing intelligent mine management and control based on digital twins. Background Art

[0002] With the continuous expansion of open-pit mining, the problem of slope stability has become increasingly prominent. Once the slope is deformed and is not detected and handled in time, it may cause serious safety accidents, such as large-scale landslides and mud-rock flows, which will not only cause huge damage to mining equipment and facilities, but also threaten the lives of on-site workers. At the same time, it will also cause irreversible damage to the surrounding environment, posing a serious threat to mine safety production. Therefore, it is particularly important to effectively monitor and analyze the deformation of mine slopes.

[0003] Traditional methods of monitoring the deformation of mine slopes are relatively simple and lack of systematization, mostly relying on irregular manual inspections and simple measuring tools. This method has obvious deficiencies in the comprehensiveness, accuracy and timeliness of data collection. Due to the inability to grasp the subtle changes in the slope in a timely manner, it is impossible to effectively identify potential dangers in the early stages of slope deformation, and the best opportunity for prevention and control is missed. Moreover, there is a lack of in-depth and efficient analysis methods for the collected data, making it difficult to extract valuable information from complex data to accurately assess the degree and trend of slope deformation. Summary of the invention

[0004] The purpose of the present invention is to provide a mine intelligent management and control optimization method based on digital twins to solve the problems raised in the above background.

[0005] The purpose of the present invention can be achieved through the following technical solutions: a method for intelligent mine management and control optimization based on digital twins, including: a data acquisition unit, a digital model construction unit, a geological monitoring and analysis unit, a microseismic monitoring and analysis unit, a slope morphology analysis unit, a support monitoring and analysis unit, a slope deformation analysis unit, a deformation anomaly assessment unit, an early warning terminal and a system repository.

[0006] The data acquisition unit is used to collect three-dimensional image data, geological data, microseismic data, slope morphology data, and support data of each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and obtain the three-dimensional image data, geological data, microseismic data, slope morphology data, and support data of each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located;

[0007] The digital model building unit is used to process the data of each monitoring sub-time period in the monitoring area where the target mine slope is located corresponding to the current monitoring period, and generate a digital surface model through a DSM generation algorithm;

[0008] The geological monitoring and analysis module is used to monitor and analyze the geological data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and obtain the geological impact assessment coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period;

[0009] The microseismic monitoring and analysis unit is used to monitor and analyze the microseismic data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and obtain the microseismic impact assessment coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period;

[0010] The slope morphology analysis unit is used to monitor and analyze the slope morphology data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, thereby analyzing and obtaining the slope morphology evaluation coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period;

[0011] The support monitoring and analysis module is used to monitor and analyze the support data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and obtain the support safety assessment index of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located;

[0012] The slope deformation analysis unit is used to analyze the deformation evaluation index of the monitoring area where the target mine slope is located corresponding to the current monitoring period, and obtain the deformation degree evaluation coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period;

[0013] The deformation anomaly assessment unit is used to analyze the deformation anomaly level of the monitoring area where the target mine slope is located corresponding to the current monitoring period, and obtain the deformation anomaly level of the monitoring area where the target mine slope is located corresponding to the current monitoring period;

[0014] The early warning terminal is used to obtain early warning information corresponding to the current monitoring period in the monitoring area where the target mine slope is located, so as to support corresponding auxiliary decision-making.

[0015] Preferably, the geological data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located is monitored and analyzed, and the specific execution steps are as follows:

[0016] Divide the monitoring area where the target mine slope is located into various soil monitoring areas according to the set area division method, and obtain various soil monitoring areas corresponding to the monitoring area where the target mine slope is located;

[0017] The length, width and depth of each soil crack in each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period are obtained, and the sum of them is calculated to obtain the total length, total width and total depth of soil cracks in each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period, which are recorded as i represents the number of each soil monitoring area, i=1,2,....,n, j represents the number of each monitoring sub-period, j=1,2,....,m;

[0018] According to the formula Calculate the soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period L', D', H' respectively represent the set total length of reference soil cracks, the total width of reference soil cracks, and the total depth of reference soil cracks; ΔL, ΔD, ΔH respectively represent the set allowable total length difference of soil cracks, the total width difference of soil cracks, and the total depth difference of soil cracks; a1, a2, a3 respectively represent the weight factors corresponding to the set total length of soil cracks, the total width of soil cracks, and the total depth of soil cracks; β1 represents the correction factor corresponding to the set soil crack state compliance index;

[0019] The soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period is calculated on average, and the average soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to the current monitoring period is obtained, which is recorded as λ i , the average soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to the current monitoring period is subtracted from the set soil crack state compliance index to obtain the soil crack state compliance index difference, and the maximum soil crack state compliance index difference is extracted from the soil crack state compliance index differences, recorded as Δλ max ;

[0020] According to the formula Calculate the average soil stability η of the monitoring area of ​​the target mine slope corresponding to the current monitoring period, Δλ max 'represents the reference average soil crack state compliance index and the reference maximum soil crack state compliance index difference respectively, a4 and a5 represent the weight factors corresponding to the set average soil crack state compliance index and the maximum soil crack state compliance index difference respectively, β2 represents the correction factor corresponding to the set soil average stability;

[0021] From the soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located, the soil crack state compliance index corresponding to the first monitoring sub-period and the soil crack state compliance index corresponding to the last monitoring sub-period are extracted respectively, and recorded as At the same time, obtain the interval length between the first monitoring sub-period and the last monitoring sub-period, recorded as ΔT;

[0022] According to the formula Calculate the soil change rate ξ of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to the current monitoring period i ;

[0023] According to the formula Calculate the soil status evaluation index φ of the monitoring area of ​​the target mine slope corresponding to the current monitoring period 土 ,η',ξ i 'represent the set reference soil average stability and reference average soil change rate, a6 and a7 represent the weight factors corresponding to the set soil average stability and average soil change rate, respectively;

[0024] By installing a steel ruler water level gauge in the observation well, the groundwater level of each monitoring sub-period in the current monitoring period corresponding to the monitoring area of ​​the target mine slope is obtained, and the maximum and minimum water levels are extracted from them. At the same time, the interval length of the maximum and minimum water levels corresponding to the monitoring interval time period is obtained, and the groundwater level change rate of the monitoring area of ​​the target mine slope corresponding to the current monitoring period is obtained by analysis;

[0025] Obtain the groundwater seepage velocity of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, extract the maximum seepage velocity and the minimum seepage velocity, and obtain the interval length of the maximum seepage velocity and the minimum seepage velocity corresponding to the monitoring interval time period, thereby analyzing the change rate of groundwater seepage velocity in the monitoring area where the target mine slope is located corresponding to the current monitoring period;

[0026] The groundwater level change rate and groundwater seepage velocity change rate of the target mine slope monitoring area corresponding to the current monitoring period are comprehensively analyzed to obtain the groundwater impact index φ of the target mine slope monitoring area corresponding to the current monitoring period. 水 ;

[0027] According to the formula Calculate the geological condition assessment coefficient ω of the monitoring area where the target mine slope is located corresponding to the current monitoring period 地质 , b3 and b4 represent the weight factors corresponding to the set soil state assessment index and groundwater impact index respectively.

[0028] Preferably, the microseismic data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located is monitored and analyzed, and the specific execution steps are as follows:

[0029] Obtain the microseismic waveform diagram of each microseismic event in each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and extract the amplitude, frequency and duration from it, thereby obtaining the amplitude, frequency and duration of the microseismic waveform of each microseismic event in each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and obtain the microseismic waveform coefficient of each microseismic event in each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located through analysis;

[0030] The number of microseismic events is used as the horizontal coordinate, and the microseismic waveform coefficient is used as the vertical coordinate. Based on this, a two-dimensional coordinate system of microseismic waveform coefficients is constructed, and the microseismic waveform coefficients obtained from each microseismic event are plotted on the two-dimensional coordinate system of microseismic waveform coefficients by plotting points.

[0031] A microseismic waveform coefficient reference line is set on the two-dimensional coordinate system of the microseismic waveform coefficient, and the points above the microseismic waveform coefficient reference line are marked as abnormal points, and the points on and below the microseismic waveform coefficient reference line are marked as normal points. The number of normal points and the number of abnormal points in the monitoring area where the target mine slope is located corresponding to each monitoring sub-time period in the current monitoring period are statistically obtained. The microseismic anomaly degree of each monitoring sub-time period in the monitoring area where the target mine slope is located corresponding to the current monitoring period is obtained through analysis.

[0032] Obtain the magnitude and released energy of each microseismic event in each monitoring sub-time period in the monitoring area where the target mine slope is located, which are recorded as The mean values ​​are calculated respectively to obtain the mean magnitude and mean released energy of the microseismic events in each monitoring sub-time period in the monitoring area of ​​the target mine slope corresponding to the current monitoring period, which are recorded as M i 、E i , according to the formula Calculate the intensity of microseismicity in each monitoring sub-time period in the monitoring area of ​​the target mine slope corresponding to the current monitoring period M' and E' represent the reference mean magnitude and reference mean released energy, respectively; c6, c7, c8, and c9 represent the weight factors corresponding to the set magnitude standard deviation, released energy standard deviation, mean magnitude, and mean released energy, respectively; α3 represents the correction factor corresponding to the set microseismic intensity;

[0033] According to the formula Calculate the microseismic impact assessment coefficient ω of the monitoring area where the target mine slope is located corresponding to the current monitoring period 微震, f1 and f2 represent the coefficient factors corresponding to the set microseismic stability and microseismic intensity, respectively.

[0034] Preferably, the slope morphology data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located is monitored and analyzed, and the specific execution steps are as follows:

[0035] The slope morphology monitoring points are arranged in the monitoring area where the target mine slope is located according to the set method, and the slope morphology monitoring points corresponding to the monitoring area where the target mine slope is located are obtained;

[0036] The horizontal displacement and vertical displacement of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period are obtained, and the horizontal displacement corresponding to the first monitoring sub-period and the horizontal displacement corresponding to the last monitoring sub-period are respectively extracted from the horizontal displacement of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period, and the vertical displacement corresponding to the first monitoring sub-period and the vertical displacement corresponding to the last monitoring sub-period are respectively extracted from the vertical displacement of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period, and the interval length of the time period corresponding to the first monitoring sub-period and the last monitoring sub-period is obtained at the same time, thereby analyzing and obtaining the horizontal displacement change rate and vertical displacement change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period;

[0037] The average value of the horizontal displacement change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period is calculated to obtain the average horizontal displacement change rate of the monitoring area where the target mine slope is located corresponding to the current monitoring period. The average value of the vertical displacement change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period is calculated to obtain the average vertical displacement change rate of the monitoring area where the target mine slope is located corresponding to the current monitoring period;

[0038] The displacement change degree θ of the monitoring area where the target mine slope is located corresponding to the current monitoring period is obtained through comprehensive analysis 位移 ;

[0039] Obtain the slope of each monitoring sub-period in the current monitoring period of each slope morphology monitoring point in the monitoring area where the target mine slope is located, and extract the slope corresponding to the first monitoring sub-period and the slope corresponding to the last monitoring sub-period from the slope of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period, and obtain the interval length of the time period corresponding to the first monitoring sub-period and the last monitoring sub-period, and obtain the slope change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period through analysis;

[0040] The slope change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period is averaged to obtain the average slope change rate of the monitoring area where the target mine slope is located corresponding to the current monitoring period. At the same time, the maximum slope change rate is extracted from the slope change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period to obtain the maximum slope change rate of the monitoring area where the target mine slope is located corresponding to the current monitoring period;

[0041] The slope change degree θ of the target mine slope in the monitoring area corresponding to the current monitoring period is obtained through comprehensive analysis. 坡度 ;

[0042] According to the formula Calculate the slope morphology assessment coefficient ω of the target mine slope in the monitoring area corresponding to the current monitoring period 形态 , w1 and w2 represent the weight factors corresponding to the set reference displacement change degree and slope change degree respectively.

[0043] The support monitoring and analysis module is used to monitor and analyze the support data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and obtain the support safety assessment index of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located;

[0044] Preferably, the support data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located is monitored and analyzed, and the specific execution steps are as follows:

[0045] The supporting structure area corresponding to the monitoring area where the target mine slope is located is obtained, and the apparent image of the supporting structure area in the monitoring area where the target mine slope is located in the last monitoring sub-period in the current monitoring period is collected, and the abnormal area in the supporting structure area is identified from the image, wherein the abnormal area includes a fracture area, a detached area, a corroded area, and a damaged area, thereby extracting the number of abnormal areas and the area of ​​each abnormal area, and summing up the area of ​​each abnormal area to obtain the total area of ​​the abnormal area, which is used as the apparent abnormal area number and the apparent abnormal area total area of ​​the supporting structure area in the monitoring area where the target mine slope is located corresponding to the current monitoring period, thereby analyzing and obtaining the apparent defect degree γ of the supporting structure area in the monitoring area where the target mine slope is located corresponding to the current monitoring period 表观 ;

[0046] Obtain the ultrasonic echo signal waveform of the last monitoring sub-period in the current monitoring period of the supporting structure area of ​​the monitoring area where the target mine slope is located, and extract the amplitude therefrom as the amplitude of the ultrasonic echo signal of the supporting structure area of ​​the monitoring area where the target mine slope is located in the current monitoring period, denoted as A, and compare it with the reference ultrasonic echo signal waveform of the supporting structure area stored in the system repository to obtain the ultrasonic echo signal overlap waveform length of the last monitoring sub-period in the current monitoring period of the supporting structure area of ​​the monitoring area where the target mine slope is located, as the ultrasonic echo signal overlap waveform length of the supporting structure area of ​​the monitoring area where the target mine slope is located in the current monitoring period, denoted as L;

[0047] At the same time, the values ​​of amplitude and waveform length in the reference ultrasonic echo signal waveform diagram of the support structure area stored in the system library are obtained, which are recorded as A' and L 波形 ';

[0048] According to the formula Calculate the internal defect degree γ of the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the current monitoring period 内部 , x3 and x4 represent the weight factors corresponding to the number of apparent abnormal regions and the total area of ​​apparent abnormal regions of the set pair, respectively, and μ2 represents the set apparent defect degree correction factor;

[0049] The stress and strain detection points are arranged in the support structure area to obtain the stress and strain detection points in the support structure area corresponding to the monitoring area where the target mine slope is located;

[0050] The stress and strain of each stress-strain detection point in the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the last monitoring sub-period in the current monitoring period are obtained, and the strain-stress curve is drawn with strain as the horizontal coordinate and stress as the vertical coordinate. Thus, the strain-stress curve of each stress-strain detection point in the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the last monitoring sub-period in the current monitoring period is obtained, which is used as the strain-stress curve of each stress-strain detection point in the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the current monitoring period, and the stress extreme value and strain extreme value are extracted from it, and the maximum stress extreme value is subtracted from the reference yield strength set by the support structure stored in the system repository to obtain the stress reserve value, and the strain extreme value is subtracted from the system repository to obtain the stress reserve value. The strain margin value is obtained by subtracting the maximum allowable strain set by the support structure stored in the system storage, thereby obtaining the stress reserve value and strain margin value of each stress-strain detection point in the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the last monitoring sub-period in the current monitoring period, as the stress reserve value and strain margin value of each stress-strain detection point in the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the current monitoring period, and the minimum stress reserve value and the minimum strain margin value are extracted from them, and the minimum stress reserve value and the minimum strain margin value of the support structure area in the monitoring area where the target mine slope is located corresponding to the current monitoring period are obtained, thereby analyzing and obtaining the support instability assessment coefficient ω of the monitoring area where the target mine slope is located corresponding to the current monitoring period 支护 .

[0051] The slope deformation analysis unit is used to analyze the deformation evaluation index of the monitoring area where the target mine slope is located corresponding to the current monitoring period, and obtain the deformation degree evaluation coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period;

[0052] Preferably, the deformation assessment index of the monitoring area where the target mine slope is located corresponding to the current monitoring period is analyzed, and the specific analysis method is as follows:

[0053] According to the formula Calculate the deformation degree assessment coefficient ψ of the monitoring area where the target mine slope is located corresponding to the current monitoring period, z1, z2, z3, z4 are respectively represented as the assessment factors corresponding to the set geological condition assessment coefficient, microseismic impact assessment coefficient, slope morphology assessment coefficient, and support instability assessment coefficient;

[0054] The deformation degree assessment index of the monitoring area where the target mine slope is located corresponding to the current monitoring period is compared with the set deformation degree assessment index threshold. If the deformation degree assessment index of the monitoring area where the target mine slope is located corresponding to the current monitoring period is greater than the set deformation degree assessment index threshold, it is determined that the deformation of the monitoring area where the target mine slope is located corresponding to the current monitoring period is normal; otherwise, it is determined that the deformation of the monitoring area where the target mine slope is located corresponding to the current monitoring period is abnormal.

[0055] Preferably, the deformation anomaly level of the monitoring area where the target mine slope is located corresponding to the current monitoring period is analyzed, and the specific analysis method is as follows:

[0056] If the monitoring area where the target mine slope is located is abnormally deformed during the current monitoring period, the deformation degree assessment index of the monitoring area where the target mine slope is located corresponding to the current monitoring period is subtracted from the set deformation degree assessment index threshold to obtain the deformation degree assessment index difference of the monitoring area where the target mine slope is located corresponding to the current monitoring period, and the deformation degree assessment index difference of the monitoring area where the target mine slope is located corresponding to the current monitoring period is matched with the deformation degree assessment index difference corresponding to each set deformation abnormality level to obtain the deformation abnormality level of the monitoring area where the target mine slope is located corresponding to the current monitoring period;

[0057] The deformation anomaly level of the monitoring area where the target mine slope is located corresponding to the current monitoring period constitutes the early warning information of the monitoring area where the target mine slope is located corresponding to the current monitoring period.

[0058] Beneficial effects of the present invention:

[0059] The present invention monitors and analyzes the geological conditions, microseismic impacts, slope morphology, and support instability of the monitoring area where the target mine slope is located corresponding to the current monitoring period, and calculates the deformation degree assessment coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period. It intuitively displays the deformation of the target mine slope, effectively solves the problems of tediousness, time and labor consumption of periodic inspections and measurements of mine management personnel, reduces the subjective error of mine slope deformation analysis to a certain extent, ensures the accuracy, reliability and rationality of the mine slope deformation analysis results, and effectively ensures the mining safety of mine personnel.

[0060] The present invention monitors and analyzes soil stability, soil change rate, groundwater level change and groundwater seepage velocity change, thereby analyzing the geological conditions corresponding to the monitoring area where the target mine slope is located, thereby realizing multi-dimensional analysis of the geological conditions of the mine slope, thereby ensuring the credibility and accuracy of the geological condition analysis results of the monitoring area where the target mine slope is located, and providing strong data support for the subsequent determination and analysis of the deformation of the target mine slope.

[0061] The present invention performs deformation anomaly level analysis when the deformation degree assessment coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period is greater than a threshold value, and performs corresponding early warning operations according to the deformation anomaly level. This not only improves the real-time and timeliness of mine managers' awareness of slope deformation risks, but also improves the efficiency of mine managers in handling slope deformation, and reduces the safety hazards caused by slope deformation. BRIEF DESCRIPTION OF THE DRAWINGS

[0062] The present invention will be further described below in conjunction with the accompanying drawings.

[0063] Figure 1 It is the overall system block diagram of the present invention.

[0064] Figure 2 It is a flow chart of the present invention. DETAILED DESCRIPTION

[0065] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0066] like Figure 1 and Figure 2 As shown, the present invention is a method for intelligent mine management and control optimization based on digital twins, including: a data acquisition unit, a digital model construction unit, a geological monitoring and analysis unit, a microseismic monitoring and analysis unit, a slope morphology analysis unit, a support monitoring and analysis unit, a slope deformation analysis unit, a deformation anomaly assessment unit, an early warning terminal and a system repository.

[0067] A data acquisition unit is used to collect three-dimensional image data, geological data, microseismic data, slope morphology data, and support data of each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and obtain the three-dimensional image data, geological data, microseismic data, slope morphology data, and support data of each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located;

[0068] It should be noted that high-resolution oblique photography equipment such as drones is used to perform oblique photography of the monitoring area where the target mine slope is located according to the set shooting angle, flight altitude, and shooting interval, thereby collecting three-dimensional image data of the monitoring area where the target mine slope is located corresponding to each monitoring sub-time period in the current monitoring period;

[0069] The three-dimensional information and deformation data of the slope are obtained through other sensors (such as lidar, GNSS, geological radar, etc.), where the deformation data includes geological data, microseismic data, slope morphology data, and support data.

[0070] A digital model building unit is used to process the data of each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and generate a digital surface model through a DSM generation algorithm;

[0071] It should be noted that data processing includes data cleaning, data correction, and data fusion;

[0072] The geological monitoring and analysis module is used to monitor and analyze the geological data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located. The specific execution steps are as follows:

[0073] Divide the monitoring area where the target mine slope is located into various soil monitoring areas according to the set area division method, and obtain various soil monitoring areas corresponding to the monitoring area where the target mine slope is located;

[0074] The length, width and depth of each soil crack in each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period are obtained, and the sum of them is calculated to obtain the total length, total width and total depth of soil cracks in each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period, which are recorded as i represents the number of each soil monitoring area, i=1,2,....,n, j represents the number of each monitoring sub-period, j=1,2,....,m;

[0075] According to the formula Calculate the soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period L', D', H' respectively represent the set total length of reference soil cracks, the total width of reference soil cracks, and the total depth of reference soil cracks; ΔL, ΔD, ΔH respectively represent the set allowable total length difference of soil cracks, the total width difference of soil cracks, and the total depth difference of soil cracks; a1, a2, a3 respectively represent the weight factors corresponding to the set total length of soil cracks, the total width of soil cracks, and the total depth of soil cracks; β1 represents the correction factor corresponding to the set soil crack state compliance index;

[0076] The soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period is calculated on average, and the average soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to the current monitoring period is obtained, which is recorded as λ i , the average soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to the current monitoring period is subtracted from the set soil crack state compliance index to obtain the soil crack state compliance index difference, and the maximum soil crack state compliance index difference is extracted from the soil crack state compliance index differences, recorded as Δλ max ;

[0077] According to the formula Calculate the average soil stability η of the monitoring area of ​​the target mine slope corresponding to the current monitoring period, Δλ max 'represents the reference average soil crack state compliance index and the reference maximum soil crack state compliance index difference respectively, a4 and a5 represent the weight factors corresponding to the set average soil crack state compliance index and the maximum soil crack state compliance index difference respectively, β2 represents the correction factor corresponding to the set soil average stability;

[0078] From the soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located, the soil crack state compliance index corresponding to the first monitoring sub-period and the soil crack state compliance index corresponding to the last monitoring sub-period are extracted respectively, and recorded as At the same time, obtain the interval length between the first monitoring sub-period and the last monitoring sub-period, recorded as ΔT;

[0079] According to the formula Calculate the soil change rate ξ of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to the current monitoring period i ;

[0080] According to the formula Calculate the soil status evaluation index φ of the monitoring area of ​​the target mine slope corresponding to the current monitoring period 土 ,η',ξ i 'represent the set reference soil average stability and reference average soil change rate, a6 and a7 represent the weight factors corresponding to the set soil average stability and average soil change rate, respectively;

[0081] The groundwater level of the monitoring area of ​​the target mine slope corresponding to each monitoring period in the current monitoring period is obtained by installing a steel ruler water level gauge in the observation well, and the maximum and minimum water levels are extracted from them. Thus, the maximum and minimum water levels of the monitoring area of ​​the target mine slope corresponding to the current monitoring period are obtained, which are recorded as sh max ,sh min , and at the same time obtain the interval length of the monitoring interval time period corresponding to the maximum water level and the minimum water level, recorded as T1;

[0082] According to the formula Calculate the groundwater level change rate κ of the monitoring area where the target mine slope is located corresponding to the current monitoring period 水位 ;

[0083] Obtain the groundwater seepage velocity of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, extract the maximum seepage velocity and the minimum seepage velocity from them, and thus obtain the maximum seepage velocity and the minimum seepage velocity of the monitoring area where the target mine slope is located corresponding to the current monitoring period, which are recorded as sv max 、sv min , and at the same time, obtain the interval length of the monitoring interval time period corresponding to the maximum seepage velocity and the minimum seepage velocity, recorded as T2;

[0084] According to the formula Calculate the groundwater seepage velocity change rate κ of the target mine slope monitoring area corresponding to the current monitoring period 水速 ;

[0085] According to the formula Calculate the groundwater impact index φ of the monitoring area of ​​the target mine slope corresponding to the current monitoring period 水 , κ 水位 ',κ 水速 'represent the set reference groundwater level change rate and reference groundwater seepage velocity change rate, b1 and b2 represent the weight factors corresponding to the set groundwater level change rate and groundwater seepage velocity change rate, respectively;

[0086] According to the formula Calculate the geological condition assessment coefficient ω of the monitoring area where the target mine slope is located corresponding to the current monitoring period 地质 , b3 and b4 represent the weight factors corresponding to the set soil state assessment index and groundwater impact index respectively.

[0087] In a specific embodiment, the present invention monitors and analyzes soil stability, soil change rate, groundwater level change, and groundwater seepage velocity change, thereby analyzing the geological conditions corresponding to the monitoring area where the target mine slope is located, thereby realizing multi-dimensional analysis of the geological conditions of the mine slope, thereby ensuring the credibility and accuracy of the geological condition analysis results of the monitoring area where the target mine slope is located, and providing strong data support for the subsequent determination and analysis of the deformation of the target mine slope.

[0088] The microseismic monitoring and analysis unit is used to monitor and analyze the microseismic data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located. The specific execution steps are as follows:

[0089] The microseismic waveforms of each microseismic event in each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located are obtained, and the amplitude, frequency and duration are extracted from them. The amplitude, frequency and duration of the microseismic waveforms of each microseismic event in each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located are obtained, which are recorded as f represents the number of each microseismic event, f = 1, 2, ..., g;

[0090] According to the formula Calculate the microseismic waveform coefficients of each microseismic event in each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located c1, c2, and c3 represent the weight factors corresponding to the set amplitude, frequency, and duration, respectively, and α1 represents the correction factor corresponding to the set microseismic waveform coefficient;

[0091] The number of microseismic events is used as the horizontal coordinate, and the microseismic waveform coefficient is used as the vertical coordinate. Based on this, a two-dimensional coordinate system of microseismic waveform coefficients is constructed, and the microseismic waveform coefficients obtained from each microseismic event are plotted on the two-dimensional coordinate system of microseismic waveform coefficients by plotting points.

[0092] A microseismic waveform coefficient reference line is set on the two-dimensional coordinate system of the microseismic waveform coefficient, and the points above the microseismic waveform coefficient reference line are marked as abnormal points, and the points on and below the microseismic waveform coefficient reference line are marked as normal points. The number of normal points and the number of abnormal points in each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located are statistically obtained, which are recorded as ZN j , YN j ;

[0093] According to the formula Calculate the microseismic anomaly of the target mine slope in the monitoring area corresponding to each monitoring sub-time period in the current monitoring period e is represented as a natural constant, c4 and c5 are weight factors corresponding to the set number of normal points and the number of abnormal points, respectively, and α2 is represented as a correction factor corresponding to the set microseismic anomaly;

[0094] Obtain the magnitude and released energy of each microseismic event in each monitoring sub-time period in the monitoring area where the target mine slope is located, which are recorded as The mean values ​​are calculated respectively to obtain the mean magnitude and mean released energy of the microseismic events in each monitoring sub-time period in the monitoring area of ​​the target mine slope corresponding to the current monitoring period, which are recorded as M i 、E i , according to the formula Calculate the intensity of microseismicity in each monitoring sub-time period in the monitoring area of ​​the target mine slope corresponding to the current monitoring period M' and E' represent the reference mean magnitude and reference mean released energy, respectively; c6, c7, c8, and c9 represent the weight factors corresponding to the set magnitude standard deviation, released energy standard deviation, mean magnitude, and mean released energy, respectively; α3 represents the correction factor corresponding to the set microseismic intensity;

[0095] According to the formula Calculate the microseismic impact assessment coefficient ω of the monitoring area where the target mine slope is located corresponding to the current monitoring period 微震 , f1 and f2 represent the coefficient factors corresponding to the set microseismic stability and microseismic intensity, respectively.

[0096] In a specific embodiment, the present invention monitors and analyzes the microseismic anomaly and microseismic intensity of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, thereby analyzing the microseismic impact assessment coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period, and effectively assessing the impact of microseismic events on the mine slope.

[0097] The slope morphology analysis unit is used to monitor and analyze the slope morphology data of the target mine slope in the monitoring area corresponding to each monitoring sub-period in the current monitoring period. The specific execution steps are as follows:

[0098] The slope morphology monitoring points are arranged in the monitoring area where the target mine slope is located according to the set method, and the slope morphology monitoring points corresponding to the monitoring area where the target mine slope is located are obtained;

[0099] Obtain the horizontal displacement and vertical displacement of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period, and extract the horizontal displacement corresponding to the first monitoring sub-period and the horizontal displacement corresponding to the last monitoring sub-period from the horizontal displacement of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period, respectively, and record them as From the vertical displacements of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period, the vertical displacement corresponding to the first monitoring sub-period and the vertical displacement corresponding to the last monitoring sub-period are extracted and recorded as p represents the number of each slope morphology monitoring point, p = 1, 2, ..., v, and the interval length of the corresponding time period between the first monitoring sub-period and the last monitoring sub-period is obtained at the same time, recorded as ΔT;

[0100] It should be noted that the horizontal displacement and the vertical displacement are respectively monitored by displacement sensors installed in the horizontal direction and the vertical direction at each slope morphology monitoring point.

[0101] According to the formula Calculate the horizontal displacement change rate LB of each slope morphology monitoring point in the monitoring area of ​​the target mine slope corresponding to the current monitoring period p ;

[0102] According to the formula Calculate the vertical displacement change rate CB of each slope morphology monitoring point in the monitoring area of ​​the target mine slope corresponding to the current monitoring period p ;

[0103] The average value of the horizontal displacement change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period is calculated to obtain the average horizontal displacement change rate of the monitoring area where the target mine slope is located corresponding to the current monitoring period, recorded as LB; the average value of the vertical displacement change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period is calculated to obtain the average vertical displacement change rate of the monitoring area where the target mine slope is located corresponding to the current monitoring period, recorded as CB;

[0104] According to the formula Calculate the displacement change θ of the monitoring area where the target mine slope is located corresponding to the current monitoring period 位移, LB', CB', LP', CP' represent the set reference horizontal displacement change rate, reference vertical displacement change rate, reference horizontal displacement change rate standard deviation, reference vertical displacement change rate standard deviation, d1, d2, d3, d4 represent the set horizontal displacement change rate, vertical displacement change rate, horizontal displacement change rate standard deviation, vertical displacement change rate standard deviation corresponding weight factors, δ1 represents the set displacement change correction factor;

[0106] Obtain the slope of each monitoring sub-period in the current monitoring period of each slope morphology monitoring point in the monitoring area where the target mine slope is located, and extract the slope corresponding to the first monitoring sub-period and the slope corresponding to the last monitoring sub-period from the slope of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period, respectively, and record them as At the same time, obtain the interval length between the first monitoring sub-period and the last monitoring sub-period, recorded as ΔT;

[0107] According to the formula Calculate the slope change rate PD of each slope morphology monitoring point in the monitoring area of ​​the target mine slope corresponding to the current monitoring period p ;

[0108] The average slope change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period is calculated to obtain the average slope change rate of the monitoring area where the target mine slope is located corresponding to the current monitoring period, which is recorded as At the same time, the maximum slope change rate is extracted from the slope change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period, and the maximum slope change rate of the monitoring area where the target mine slope is located corresponding to the current monitoring period is obtained, which is recorded as PD max ;

[0109] According to the formula Calculate the slope change θ of the target mine slope in the monitoring area corresponding to the current monitoring period 坡度 , P.D. p '、 PD max 'respectively represent the set reference slope change rate, reference mean slope change rate, and reference maximum slope change rate, d8, d9, and d10 respectively represent the weight factors corresponding to the set slope change rate, mean slope change rate, and maximum slope change rate, and δ2 represents the set slope change correction factor;

[0110] According to the formula Calculate the slope morphology assessment coefficient ω of the target mine slope in the monitoring area corresponding to the current monitoring period 形态, w1 and w2 represent the weight factors corresponding to the set reference displacement change and slope change, respectively;

[0111] In a specific embodiment, the present invention monitors and analyzes the horizontal displacement changes, vertical displacement changes and slope changes in the monitoring area where the target mine slope is located corresponding to the current monitoring period, and obtains the slope morphology assessment coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period. To a certain extent, it intuitively displays the law of mine slope morphology changes, improves the adequacy of the deformation degree assessment and analysis of the monitoring area where the target mine slope is located, and at the same time ensures the mining safety of mine personnel, so that the mine slope deformation analysis achieves the expected effect.

[0112] The support monitoring and analysis module is used to monitor and analyze the support data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located. The specific execution steps are as follows:

[0113] Obtain the supporting structure area corresponding to the monitoring area where the target mine slope is located, collect the apparent image of the last monitoring sub-period in the current monitoring period of the supporting structure area in the monitoring area where the target mine slope is located, and identify the abnormal area in the supporting structure area from the image, wherein the abnormal area includes a fracture area, a detached area, a corroded area, and a damaged area, thereby extracting the number of abnormal areas and the area of ​​each abnormal area, and summing up the area of ​​each abnormal area to obtain the total area of ​​the abnormal area, thereby serving as the apparent number of abnormal areas and the total area of ​​the apparent abnormal areas corresponding to the current monitoring period of the supporting structure area in the monitoring area where the target mine slope is located, which are recorded as YN and YS respectively;

[0114] According to the formula γ 表观 =(e YN×x1+YS×x2 )×μ1 calculates the apparent defect degree γ of the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the current monitoring period 表观 , x1 and x2 represent the weight factors corresponding to the set number of apparent abnormal regions and the total area of ​​apparent abnormal regions, respectively, and μ1 represents the set apparent defect correction factor;

[0115] Obtain the ultrasonic echo signal waveform of the last monitoring sub-period in the current monitoring period of the supporting structure area of ​​the monitoring area where the target mine slope is located, and extract the amplitude therefrom as the amplitude of the ultrasonic echo signal of the supporting structure area of ​​the monitoring area where the target mine slope is located in the current monitoring period, denoted as A, and compare it with the reference ultrasonic echo signal waveform of the supporting structure area stored in the system repository to obtain the ultrasonic echo signal overlap waveform length of the last monitoring sub-period in the current monitoring period of the supporting structure area of ​​the monitoring area where the target mine slope is located, as the ultrasonic echo signal overlap waveform length of the supporting structure area of ​​the monitoring area where the target mine slope is located in the current monitoring period, denoted as L;

[0116] At the same time, the values ​​of amplitude and waveform length in the reference ultrasonic echo signal waveform diagram of the support structure area stored in the system library are obtained, which are recorded as A' and L 波形 ';

[0117] According to the formula Calculate the internal defect degree γ of the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the current monitoring period 内部 , x3 and x4 represent the weight factors corresponding to the number of apparent abnormal regions and the total area of ​​apparent abnormal regions of the set pair, respectively, and μ2 represents the set apparent defect degree correction factor;

[0118] The stress and strain detection points are arranged in the support structure area to obtain the stress and strain detection points in the support structure area corresponding to the monitoring area where the target mine slope is located;

[0119] Obtain the stress and strain of each stress-strain detection point in the support structure area in the monitoring area where the target mine slope is located, corresponding to the last monitoring sub-period in the current monitoring period, and draw the strain-stress curve with strain as the horizontal coordinate and stress as the vertical coordinate. Thus, the strain-stress curve of each stress-strain detection point in the support structure area in the monitoring area where the target mine slope is located corresponding to the last monitoring sub-period in the current monitoring period is obtained, which is used as the strain-stress curve of each stress-strain detection point in the support structure area in the monitoring area where the target mine slope is located corresponding to the current monitoring period, and extract the stress extreme value and strain extreme value from it, and make a difference between the maximum stress extreme value and the reference yield strength set by the support structure stored in the system repository. The stress reserve value is obtained, and the strain extreme value is subtracted from the maximum allowable strain set by the support structure stored in the system repository to obtain the strain margin value. Thus, the stress reserve value and strain margin value of each stress-strain detection point in the support structure area in the monitoring area where the target mine slope is located corresponding to the last monitoring sub-period in the current monitoring period are obtained, which are used as the stress reserve value and strain margin value of each stress-strain detection point in the support structure area in the monitoring area where the target mine slope is located corresponding to the current monitoring period, and the minimum stress reserve value and the minimum strain margin value are extracted from them, and the minimum stress reserve value and the minimum strain margin value of the support structure area in the monitoring area where the target mine slope is located corresponding to the current monitoring period are obtained, which are recorded as Δσ respectively. min , Δε min ;

[0120] According to the formula Calculate the support instability assessment coefficient ω of the monitoring area where the target mine slope is located corresponding to the current monitoring period 支护 ,

[0122] γ 表观 ' , γ 内部 ' , Δσ min ' , Δε min ' They respectively represent the set reference apparent defect, reference internal defect, reference minimum stress reserve value, and reference minimum strain margin value; y1, y2, y3, and y4 respectively represent the weight factors corresponding to the set apparent defect, internal defect, minimum stress reserve value, and minimum strain margin value.

[0123] In a specific embodiment, the present invention monitors and analyzes the apparent defects, internal defects, stresses and strains of the support structure area in the monitoring area where the target mine slope is located corresponding to the current monitoring period, and obtains the support instability assessment coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period, which helps to timely discover the instability of mine slope support, and makes the deformation analysis of the monitoring area where the target mine slope is located corresponding to the current monitoring period more targeted and effective.

[0124] The slope deformation analysis unit is used to analyze the deformation evaluation index of the monitoring area of ​​the target mine slope corresponding to the current monitoring period. The specific analysis method is as follows:

[0125] According to the formula Calculate the deformation degree assessment coefficient ψ of the monitoring area where the target mine slope is located corresponding to the current monitoring period, z1, z2, z3, z4 are respectively represented as the assessment factors corresponding to the set geological condition assessment coefficient, microseismic impact assessment coefficient, slope morphology assessment coefficient, and support instability assessment coefficient;

[0126] The deformation degree assessment index of the monitoring area where the target mine slope is located corresponding to the current monitoring period is compared with the set deformation degree assessment index threshold. If the deformation degree assessment index of the monitoring area where the target mine slope is located corresponding to the current monitoring period is greater than the set deformation degree assessment index threshold, it is determined that the deformation of the monitoring area where the target mine slope is located corresponding to the current monitoring period is normal; otherwise, it is determined that the deformation of the monitoring area where the target mine slope is located corresponding to the current monitoring period is abnormal.

[0127] It should be noted that the threshold of the deformation degree assessment index is dynamically adjusted according to different seasons and different climates. A smaller threshold is set when the risk of geological disasters is high, and a larger threshold is set vice versa.

[0128] The deformation anomaly assessment unit is used to analyze the deformation anomaly level of the monitoring area of ​​the target mine slope corresponding to the current monitoring period. The specific analysis method is as follows:

[0129] If the monitoring area where the target mine slope is located is abnormally deformed during the current monitoring period, the deformation degree assessment index of the monitoring area where the target mine slope is located corresponding to the current monitoring period is subtracted from the set deformation degree assessment index threshold to obtain the deformation degree assessment index difference of the monitoring area where the target mine slope is located corresponding to the current monitoring period, and the deformation degree assessment index difference of the monitoring area where the target mine slope is located corresponding to the current monitoring period is matched with the deformation degree assessment index difference corresponding to each set deformation abnormality level to obtain the deformation abnormality level of the monitoring area where the target mine slope is located corresponding to the current monitoring period;

[0130] The deformation anomaly level of the monitoring area where the target mine slope is located corresponding to the current monitoring period constitutes the early warning information of the monitoring area where the target mine slope is located corresponding to the current monitoring period.

[0131] Early warning terminal: used to obtain early warning information corresponding to the current monitoring period in the monitoring area where the target mine slope is located, so as to support corresponding auxiliary decision-making.

[0132] In a specific embodiment, the present invention performs deformation anomaly level analysis when the deformation degree assessment coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period is greater than a threshold value, and performs corresponding early warning operations according to the deformation anomaly level. This not only improves the real-time and timeliness of mine managers' awareness of slope deformation risks, but also improves the efficiency of mine managers in handling slope deformation and reduces safety hazards caused by slope deformation.

[0133] The system repository is used to store the reference ultrasonic echo signal waveform of the support structure area, store the reference yield strength set for the support structure, and store the maximum allowable strain set for the support structure.

[0134] The above contents are merely examples and explanations of the structure of the present invention. The technicians in this technical field may make various modifications or additions to the specific embodiments described or replace them in a similar manner. As long as they do not deviate from the structure of the invention or exceed the scope defined by the claims, they should all fall within the protection scope of the present invention.

Claims

1. The intelligent mine management and control optimization method based on digital twins is characterized by: The following steps are involved: Step 1: Collect 3D image data, geological data, microseismic data, slope morphology data, and support data of each monitoring sub-time period in the current monitoring period in the monitoring area where the target mine slope is located; Step 2: Process the data of each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and generate a digital surface model through a DSM generation algorithm; Step 3: monitor and analyze the geological data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and obtain the geological impact assessment coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period; Step 4: monitor and analyze the microseismic data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and obtain the microseismic impact assessment coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period; Step 5: monitor and analyze the slope morphology data of each monitoring sub-period in the current monitoring period in the monitoring area where the target mine slope is located, and obtain the slope morphology evaluation coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period; Step 6: Monitor and analyze the support data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and obtain the support safety assessment index of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located; Step 7: Analyze the deformation assessment index of the monitoring area where the target mine slope is located corresponding to the current monitoring period to obtain the deformation degree assessment coefficient of the monitoring area where the target mine slope is located corresponding to the current monitoring period; Step 8: Analyze the deformation anomaly level of the monitoring area where the target mine slope is located corresponding to the current monitoring period to obtain the deformation anomaly level of the monitoring area where the target mine slope is located corresponding to the current monitoring period; Step 9: Obtain the early warning information of the monitoring area where the target mine slope is located corresponding to the current monitoring period, so as to support the corresponding auxiliary decision.

2. The method for optimizing intelligent mine management and control based on digital twins according to claim 1 is characterized in that: The monitoring and analysis of the geological data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located is specifically performed in the following steps: Divide the monitoring area where the target mine slope is located into various soil monitoring areas according to the set area division method, and obtain various soil monitoring areas corresponding to the monitoring area where the target mine slope is located; The length, width and depth of each soil crack in each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period are obtained, and the sum of the soil cracks is calculated to obtain the total length, total width and total depth of soil cracks in each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period. The soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period is calculated by the above parameters; The soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period is processed to obtain the average soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to the current monitoring period; the soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period is processed to obtain the soil crack state compliance index difference, and the maximum soil crack state compliance index difference is extracted from the soil crack state compliance index differences, and the soil average stability of the monitoring area where the target mine slope is located corresponding to the current monitoring period is calculated by the above parameters; The soil crack state compliance index corresponding to the first and last monitoring sub-periods is extracted from the soil crack state compliance index of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period, and the interval length of the time period corresponding to the first and last monitoring sub-periods is obtained. The soil change rate of each soil monitoring area in the monitoring area where the target mine slope is located corresponding to the current monitoring period is calculated according to the above parameters. The soil state evaluation index φ of the monitoring area where the target mine slope is located corresponding to the current monitoring period is calculated by the above parameters. 土 ; The groundwater level of each monitoring sub-period in the current monitoring period in the monitoring area where the target mine slope is located is obtained by installing a steel ruler water level gauge in the observation well, and the maximum and minimum water levels are extracted from them. At the same time, the interval length of the maximum and minimum water levels corresponding to the monitoring interval time period is obtained. The groundwater level change rate of the monitoring area where the target mine slope is located corresponding to the current monitoring period is calculated by the above parameters; Obtain the groundwater seepage velocity of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, extract the maximum seepage velocity and the minimum seepage velocity, and obtain the interval length of the maximum seepage velocity and the minimum seepage velocity corresponding to the monitoring interval time period. Calculate the groundwater seepage velocity change rate of the monitoring area where the target mine slope is located corresponding to the current monitoring period through the above parameters; The groundwater level change rate and groundwater seepage velocity change rate of the target mine slope monitoring area corresponding to the current monitoring period are calculated and analyzed to obtain the groundwater impact index φ of the target mine slope monitoring area corresponding to the current monitoring period. 水 ; According to the formula Calculate the geological condition assessment coefficient ω of the monitoring area where the target mine slope is located corresponding to the current monitoring period 地质 , b3 and b4 represent the weight factors corresponding to the set soil state assessment index and groundwater impact index respectively.

3. The method for optimizing intelligent mine management and control based on digital twins according to claim 1 is characterized in that: The monitoring and analysis of the microseismic data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located is specifically performed in the following steps: Obtain the microseismic waveform diagram of each microseismic event in each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and extract the amplitude, frequency and duration from it. Obtain the microseismic waveform coefficient of each microseismic event in each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located through the above parameter calculation and analysis; The number of microseismic events is used as the horizontal coordinate, and the microseismic waveform coefficient is used as the vertical coordinate. Based on this, a two-dimensional coordinate system of microseismic waveform coefficients is constructed, and the microseismic waveform coefficients obtained from each microseismic event are plotted on the two-dimensional coordinate system of microseismic waveform coefficients by plotting points. A microseismic waveform coefficient reference line is set on the microseismic waveform coefficient two-dimensional coordinate system, and the points above the microseismic waveform coefficient reference line are marked as abnormal points, and the points on and below the microseismic waveform coefficient reference line are marked as normal points. The number of normal points and the number of abnormal points in the monitoring area where the target mine slope is located corresponding to each monitoring sub-time period in the current monitoring period are statistically obtained. The microseismic anomaly degree of each monitoring sub-time period in the monitoring area where the target mine slope is located corresponding to the current monitoring period is calculated through the above parameters. Obtain the magnitude and released energy of each microseismic event in each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located, and calculate the average value for each microseismic event in each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located. The microseismic intensity of each monitoring sub-time period in the current monitoring period corresponding to the monitoring area where the target mine slope is located is calculated through the above parameters. j represents the number of each monitoring sub-period, j = 1, 2, ..., m; According to the formula Calculate the microseismic impact assessment coefficient ω of the monitoring area where the target mine slope is located corresponding to the current monitoring period 微震 , f1 and f2 represent the coefficient factors corresponding to the set microseismic stability and microseismic intensity, respectively.

4. The method for optimizing intelligent mine management and control based on digital twins according to claim 1 is characterized in that: The monitoring and analysis of the slope morphology data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located is specifically performed in the following steps: Obtain the horizontal displacement and vertical displacement of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period, and obtain the average horizontal displacement change rate and average vertical displacement change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period through calculation and analysis; Through the above parameter calculation and analysis, the displacement change degree of the monitoring area of ​​the target mine slope corresponding to the current monitoring period is obtained. Obtain the slope of each monitoring sub-period in the current monitoring period of each slope morphology monitoring point in the monitoring area where the target mine slope is located, and extract the slope corresponding to the first monitoring sub-period and the slope corresponding to the last monitoring sub-period from the slope of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to each monitoring sub-period in the current monitoring period, and obtain the interval length of the time period corresponding to the first monitoring sub-period and the last monitoring sub-period, and obtain the slope change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period through analysis; The slope change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period is averaged to obtain the average slope change rate of the monitoring area where the target mine slope is located corresponding to the current monitoring period. At the same time, the maximum slope change rate is extracted from the slope change rate of each slope morphology monitoring point in the monitoring area where the target mine slope is located corresponding to the current monitoring period to obtain the maximum slope change rate of the monitoring area where the target mine slope is located corresponding to the current monitoring period; Through the above parameter calculation and analysis, the slope change degree of the target mine slope in the monitoring area corresponding to the current monitoring period is obtained. According to the formula Calculate the slope morphology assessment coefficient ω of the target mine slope in the monitoring area corresponding to the current monitoring period 形态 , w1 and w2 represent the weight factors corresponding to the set reference displacement change degree and slope change degree respectively.

5. The method for optimizing intelligent mine management and control based on digital twins according to claim 1 is characterized in that: The monitoring and analysis of the support data of each monitoring sub-period in the current monitoring period corresponding to the monitoring area where the target mine slope is located is specifically performed in the following steps: The supporting structure area corresponding to the monitoring area where the target mine slope is located is obtained, and the apparent image of the supporting structure area in the monitoring area where the target mine slope is located in the last monitoring sub-period in the current monitoring period is collected, and the abnormal area in the supporting structure area is identified from the image, thereby extracting the number of abnormal areas and the area of ​​each abnormal area, and summing up the area of ​​each abnormal area to obtain the total area of ​​the abnormal area, which is used as the apparent number of abnormal areas and the total area of ​​the apparent abnormal areas corresponding to the current monitoring period of the supporting structure area in the monitoring area where the target mine slope is located. The apparent defect degree γ of the supporting structure area in the monitoring area where the target mine slope is located corresponding to the current monitoring period is obtained through the above parameter calculation and analysis. 表观 ; Obtain the ultrasonic echo signal waveform of the last monitoring sub-period in the current monitoring period of the supporting structure area of ​​the monitoring area where the target mine slope is located, and extract the amplitude therefrom as the amplitude of the ultrasonic echo signal of the supporting structure area of ​​the monitoring area where the target mine slope is located in the current monitoring period, denoted as A, and compare it with the reference ultrasonic echo signal waveform of the supporting structure area stored in the system repository to obtain the ultrasonic echo signal overlap waveform length of the last monitoring sub-period in the current monitoring period of the supporting structure area of ​​the monitoring area where the target mine slope is located, as the ultrasonic echo signal overlap waveform length of the supporting structure area of ​​the monitoring area where the target mine slope is located in the current monitoring period, denoted as L; At the same time, the values ​​of amplitude and waveform length in the reference ultrasonic echo signal waveform diagram of the support structure area stored in the system library are obtained, which are recorded as A' and L 波形 '; According to the formula Calculate the internal defect degree γ of the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the current monitoring period 内部 , x3 and x4 represent the weight factors corresponding to the number of apparent abnormal regions and the total area of ​​apparent abnormal regions of the set pair, respectively, and μ2 represents the set apparent defect degree correction factor; The stress and strain detection points are arranged in the support structure area to obtain the stress and strain detection points in the support structure area corresponding to the monitoring area where the target mine slope is located; The stress and strain of each stress-strain detection point in the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the last monitoring sub-period in the current monitoring period are obtained, and the strain-stress curve is drawn with strain as the horizontal coordinate and stress as the vertical coordinate. Thus, the strain-stress curve of each stress-strain detection point in the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the last monitoring sub-period in the current monitoring period is obtained, which is used as the strain-stress curve of each stress-strain detection point in the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the current monitoring period, and the stress extreme value and strain extreme value are extracted from it, and the maximum stress extreme value is subtracted from the reference yield strength set by the support structure stored in the system repository to obtain the stress reserve value, and the strain extreme value is subtracted from the system repository to obtain the stress reserve value. The strain margin value is obtained by subtracting the maximum allowable strain set by the support structure stored in the system storage, thereby obtaining the stress reserve value and strain margin value of each stress-strain detection point in the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the last monitoring sub-period in the current monitoring period, as the stress reserve value and strain margin value of each stress-strain detection point in the support structure area of ​​the monitoring area where the target mine slope is located corresponding to the current monitoring period, and the minimum stress reserve value and the minimum strain margin value are extracted from them, and the minimum stress reserve value and the minimum strain margin value of the support structure area in the monitoring area where the target mine slope is located corresponding to the current monitoring period are obtained, thereby analyzing and obtaining the support instability assessment coefficient ω of the monitoring area where the target mine slope is located corresponding to the current monitoring period 支护 .

6. The method for optimizing intelligent mine management and control based on digital twin according to claim 1 is characterized in that: The deformation assessment index of the monitoring area where the target mine slope is located corresponding to the current monitoring period is analyzed, and the specific analysis method is as follows: According to the formula Calculate the deformation degree assessment coefficient ψ of the monitoring area where the target mine slope is located corresponding to the current monitoring period, z1, z2, z3, z4 are respectively represented as the assessment factors corresponding to the set geological condition assessment coefficient, microseismic impact assessment coefficient, slope morphology assessment coefficient, and support instability assessment coefficient; The deformation degree assessment index of the monitoring area where the target mine slope is located corresponding to the current monitoring period is compared with the set deformation degree assessment index threshold. If the deformation degree assessment index of the monitoring area where the target mine slope is located corresponding to the current monitoring period is greater than the set deformation degree assessment index threshold, it is determined that the deformation of the monitoring area where the target mine slope is located corresponding to the current monitoring period is normal; otherwise, it is determined that the deformation of the monitoring area where the target mine slope is located corresponding to the current monitoring period is abnormal.

7. The method for optimizing intelligent mine management and control based on digital twins according to claim 1 is characterized in that: The deformation anomaly level of the monitoring area of ​​the target mine slope corresponding to the current monitoring period is analyzed, and the specific analysis method is as follows: If the monitoring area where the target mine slope is located is abnormally deformed during the current monitoring period, the deformation degree assessment index of the monitoring area where the target mine slope is located corresponding to the current monitoring period is subtracted from the set deformation degree assessment index threshold to obtain the deformation degree assessment index difference of the monitoring area where the target mine slope is located corresponding to the current monitoring period, and the deformation degree assessment index difference of the monitoring area where the target mine slope is located corresponding to the current monitoring period is matched with the deformation degree assessment index difference corresponding to each set deformation abnormality level to obtain the deformation abnormality level of the monitoring area where the target mine slope is located corresponding to the current monitoring period; The deformation anomaly level of the monitoring area where the target mine slope is located corresponding to the current monitoring period constitutes the early warning information of the monitoring area where the target mine slope is located corresponding to the current monitoring period.

Images