Road slope data monitoring method under complex geological conditions

Abstract

The invention provides a road slope data monitoring method under a complex geological condition, and relates to the technical field of road slope monitoring, and the method comprises the steps: carrying out the temperature compensation of a change rate, analyzing the thermal expansion characteristic of a material, calculating the deformation caused by the temperature change, correcting the monitored side length, and generating a temperature compensation side length; according to the temperature compensation side length and the surrounding rock rheological equation, correcting the two adjacent side lengths through the strain difference of the adjacent sides so as to obtain the final side length; combining the final side length with the reference interior angle to generate a topological structure; based on the topological structure, constructing a three-dimensional deformation index matrix containing a diagonal displacement difference, an adjacent side curvature ratio and an interior angle variation coefficient; and inputting the index matrix, the support structure strain space-time matrix and the rock mass fragmentation tensor into a space-time attention graph neural network together, and outputting a deformation evolution trend including a diagonal displacement difference, an adjacent side curvature ratio and an interior angle variation coefficient. According to the invention, high-precision real-time monitoring can be realized.

Description

Technical Field

[0001] The present invention relates to the technical field of road slope monitoring, and particularly to a method for monitoring road slope data under complex geological conditions. Background Art

[0002] In the construction of road projects in complex geological regions, slope stability monitoring is a core technical problem to ensure project safety. Especially in rock mass fracture zones, high in-situ stress soft rocks, and regions with significant temperature alternation, some traditional monitoring methods face the following technical bottlenecks: For example, existing monitoring technologies rely on the layout mode of discrete sensors (such as single-point displacement gauges and inclinometers), which is essentially a data acquisition logic based on "point-like observations".

[0003] Such methods have the following technical limitations: Single-point sensors can only reflect the local deformation state of their installation positions and cannot construct a continuous displacement field distribution model on the slope surface. For example, in a rock mass fracture zone, potential slip surfaces often exhibit heterogeneous expansion characteristics, and the deformation gradient information between discrete data points is missing, resulting in the inability to accurately reconstruct the spatial morphology of the slip surface through interpolation algorithms.

[0004] For example, in progressive failure, dynamic characteristics such as the phase difference and direction synergy of displacement rates at different positions are key indicators for judging the connectivity of slip surfaces. However, due to the lack of real-time cross-correlation calculation of multi-node data in existing technologies, it is difficult to capture such co-deformation laws; the sparse layout of discrete sensors is likely to cause monitoring blind spots. Especially at the intersection of structural planes or the interface of weak interlayers, local strain concentration may occur prior to macroscopic displacement, and traditional methods are not easy to detect such precursors of hidden damage due to the lack of distributed strain sensing capabilities. Summary of the Invention

[0005] The technical problem to be solved by the present invention is to provide a method for monitoring road slope data under complex geological conditions, which can achieve high-precision real-time monitoring.

[0006] To solve the above technical problem, the technical solution of the present invention is as follows: A method for monitoring road slope data under complex geological conditions, the method comprising: Measuring the lengths and interior angles of each side of a quadrilateral in real time to generate a reference topology; Based on the benchmark topology and the real-time distances between the vertices of the quadrilateral and the tunnel face, calculate the weight of each vertex; when the weight of a vertex is greater than the threshold, monitor the change rate of the corresponding edge in real time; perform temperature compensation on the change rate, analyze the thermal expansion characteristics of the material, calculate the deformation caused by temperature change, and correct the monitored side length to generate the temperature-compensated side length; according to the temperature-compensated side length and the surrounding rock rheological equation, correct the two adjacent side lengths through the strain difference between adjacent sides to obtain the final side length; combine the final side length with the reference interior angle to generate the topological structure. Based on the topological structure, construct a three-dimensional deformation index matrix including the diagonal displacement difference, the adjacent side curvature ratio, and the interior angle coefficient of variation; input the index matrix, the support structure strain spatio-temporal matrix, and the rock mass fragmentation tensor into the spatio-temporal attention graph neural network, and output the deformation evolution trend including the diagonal displacement difference, the adjacent side curvature ratio, and the interior angle coefficient of variation. Execute multi-level early warning decisions according to the component parameters of the deformation evolution trend.

[0007] The above scheme of the present invention has at least the following beneficial effects: By constructing the benchmark topology through real-time measurement of the side lengths and interior angles of the quadrilateral unit, it breaks through the discreteness limitation of traditional single-point monitoring and realizes the systematic description of the "planar deformation" of the slope. For example, when the joint surface inside the slope undergoes dislocation, the quadrilateral vertex weight calculation mechanism can automatically identify the key deformation areas near the tunnel face (such as the vertex of the free face triggered by the weight threshold), and track the length change rate of the corresponding edge in real time, which can improve the ability to identify the local landslide risk of the slope. Based on the thermal expansion characteristics of the material, establish a temperature-deformation coupling model, calculate the false deformation caused by temperature change in real time, eliminate the interference of environmental factors on the monitoring data, and improve the side length monitoring accuracy from ±2mm to ±0.5mm; combined with the Burgers surrounding rock rheological equation, use the strain difference between adjacent sides to dynamically adjust the influence of long-term creep, and solve the defect that the traditional method does not consider the time-dependent deformation of the rock mass. For example, in the mudstone-sandstone interlayer slope, it can accurately capture the progressive influence of the cohesion reduction caused by groundwater softening during the rainy season on the slope deformation.

[0008] Construct a topological deformation index including the diagonal displacement difference (resolution 0.1mm), the adjacent side curvature ratio (accuracy 0.01), and the interior angle coefficient of variation (sensitivity 0.5°), and combine multi-source parameters such as the support structure strain (such as the change rate of anchor cable tension 1%) and the rock mass fragmentation tensor (the resolution of the volume ratio of fragmented blocks 5%) to form a full-element monitoring matrix.

[0009] Based on the component parameters of the deformation evolution trend (such as the displacement difference rate 0.5mm / h, the curvature ratio growth rate 10% / h), realize multi-level responses from yellow warning (local deformation) to red warning (overall instability), which can reduce the false alarm rate and missed alarm rate. Description of the Drawings

[0010] Figure 1 It is a schematic flow chart of a method for monitoring road slope data under complex geological conditions provided by an embodiment of the present invention. Specific embodiments

[0011] Hereinafter, exemplary embodiments of the present disclosure will be described in more detail with reference to the accompanying drawings. Although the exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided so that the present disclosure can be more thoroughly understood and the scope of the present disclosure can be completely conveyed to those skilled in the art.

[0012] As Figure 1 shown, an embodiment of the present invention proposes a method for monitoring road slope data under complex geological conditions, and the method includes the following steps: Measure the lengths and interior angles of each side of the quadrilateral in real time to generate a reference topology. Based on the reference topology and the real-time distances between each vertex of the quadrilateral and the tunnel face, calculate the weight of each vertex; when the weight of the vertex is greater than the threshold, monitor the change rate of the corresponding side in real time; perform temperature compensation on the change rate, analyze the material thermal expansion characteristics, calculate the deformation amount caused by temperature change, and correct the monitored side length to generate temperature-compensated side lengths; according to the temperature-compensated side lengths and the surrounding rock rheological equation, correct the lengths of two adjacent sides through the strain difference between adjacent sides to obtain the final side lengths; combine the final side lengths with the reference interior angles to generate a topological structure. Based on the topological structure, construct a three-dimensional deformation index matrix including the diagonal displacement difference, the adjacent side curvature ratio, and the interior angle coefficient of variation; input the index matrix, the support structure strain spatio-temporal matrix, and the rock mass fragmentation tensor into the spatio-temporal attention graph neural network, and output the deformation evolution trend including the diagonal displacement difference, the adjacent side curvature ratio, and the interior angle coefficient of variation. Execute a multi-level early warning decision according to the component parameters of the deformation evolution trend.

[0013] In the embodiments of the present invention, a reference topology is constructed by measuring the side lengths and interior angles of quadrilateral units in real time, breaking through the discreteness limitation of traditional single-point monitoring and realizing the systematic description of the "surface deformation" of the slope. For example, when the joint plane inside the slope undergoes dislocation, the quadrilateral vertex weight calculation mechanism can automatically identify the key deformation areas near the tunnel face (such as the free face vertices triggered by the weight threshold), and track the length change rate of the corresponding sides in real time, which can improve the ability to identify the local landslide risk of the slope. Based on the thermal expansion characteristics of materials, a temperature-deformation coupling model is established to calculate the false deformation caused by temperature changes in real time, eliminating the interference of environmental factors on the monitoring data and improving the side length monitoring accuracy from ±2 mm to ±0.5 mm; combined with the Burgers surrounding rock rheological equation, the long-term creep effect is dynamically adjusted by using the strain difference between adjacent sides, solving the defect that the traditional method does not consider the time-dependent deformation of the rock mass. For example, in a mudstone-sandstone interlayer slope, the progressive effect of the cohesion reduction caused by groundwater softening during the rainy season on the slope deformation can be accurately captured. A topological deformation index including the diagonal displacement difference (resolution 0.1 mm), the adjacent side curvature ratio (accuracy 0.01), and the interior angle coefficient of variation (sensitivity 0.5°) is constructed, and multi-source parameters such as the strain of the support structure (such as the change rate of the anchor cable tension 1%) and the rock mass fragmentation tensor (the resolution of the volume ratio of fragmented blocks 5%) are combined to form a full-element monitoring matrix. Based on the component parameters of the deformation evolution trend (such as the displacement difference rate 0.5 mm / h and the curvature ratio growth rate 10% / h), a multi-level response from yellow warning (local deformation) to red warning (overall instability) can be realized, reducing the false alarm rate and the missed alarm rate.

[0014] In a preferred embodiment of the present invention, the reference topology is as follows: At the selected positions of the tunnel section, the crown settlement point P1, the left convergence point P2, the right convergence point P3, and the invert uplift point P4 are arranged. The initial lengths of the sides of the quadrilateral and the interior angles between adjacent sides are measured in real time by a three-dimensional laser rangefinder, and an initial quadrilateral geometric parameter database is established; Through the biaxial inclinometer array between P2 and P3, the change amount of the inclination angle of the side walls on both sides is collected, and the three-dimensional section scanner is synchronously triggered to perform point cloud scanning on the section where the quadrilateral is located. The scanning data is spatially registered with the laser ranging result to obtain calibration data; According to the calibration data, the lengths of the sides and the interior angles of the initial quadrilateral are dynamically corrected to generate a reference topology including the spatial coordinates of the side lines, the interior angle values, and the section curvature.

[0015] In the embodiments of the present invention, the acquisition of the initial geometric parameters, the layout of the key points, and the collection of the basic data are realized as follows: Install prisms or reflective targets at four characteristic positions on the vault (P1), left sidewall (P2), right sidewall (P3), and invert (P4) of the tunnel cross-section to form a quadrilateral monitoring unit (P1 - P2 - P3 - P4). These four points cover the top, both sides, and the bottom of the tunnel cross-section, constituting a closed geometric figure that can reflect the overall deformation of the cross-section.

[0016] Use a three-dimensional laser rangefinder (accuracy ±0.3 mm) to measure the initial distances of the 6 sides of the quadrilateral (P1P2, P1P3, P1P4, P2P3, P2P4, P3P4), and record the interior angles of adjacent sides (such as ∠P1P2P3, ∠P2P3P4, etc.) through an angle measurement module (accuracy ±5″), and establish an initial geometric parameter database containing side lengths (L 1 -L 6 ), interior angles (θ 1 -θ 4 ).

[0017] Choose a quadrilateral instead of a triangle because a quadrilateral has a higher degree of freedom in deformation sensitivity (a triangle is a rigid structure and it is difficult to capture in-plane shear deformation), and the four vertices cover the main deformation-sensitive areas of the tunnel cross-section (vault settlement, sidewall convergence, invert heave).

[0018] The four vertices correspond to the "top - side - bottom" characteristic positions of the tunnel cross-section. The constructed quadrilateral can synchronously monitor multi-mode deformations such as settlement, convergence, and heave, avoiding the one-sidedness of single-point monitoring. The laser ranging and angle measurement technologies provide initial geometric parameters with millimeter-level accuracy.

[0019] The fusion of inclination monitoring and point cloud scanning is achieved as follows: Install biaxial inclinometers (accuracy ±0.01°) at intervals of 0.5 - 1 m along the longitudinal direction of the tunnel at the sidewall position between the left convergence point P2 and the right convergence point P3 to form an inclination monitoring array, and continuously collect the change amounts (△α, △β) of the transverse inclination angle (reflecting horizontal convergence) and the longitudinal inclination angle (reflecting twist along the tunnel axis) of the sidewall.

[0020] Point cloud scanning and spatial registration: When the change amount of the inclination exceeds the threshold (such as △α > 0.05° or △β > 0.03°), trigger the three-dimensional cross-section scanner (resolution 0.5 mm) installed on the tunnel sidewall to collect point cloud data of the cross-section where the quadrilateral is located, and obtain thousands of dense three-dimensional coordinate points (x, y, z).

[0021] Spatially register the point cloud data with the prism coordinates of the laser ranging through the Iterative Closest Point (ICP) algorithm: Using the four vertices P1 - P4 as reference control points, adjust the point cloud coordinate system to align the scanned point cloud with the geometric framework of the laser ranging, eliminate the equipment installation error and coordinate system deviation, and generate calibrated three - dimensional coordinate data.

[0022] The inclinometer captures the minute tilt deformation of the side wall in real - time (for example, a 0.01° tilt corresponds to a 1.7 - mm displacement at a 10 - m distance), triggering high - precision scanning and avoiding the energy consumption and data redundancy of continuous full - time scanning; the point cloud scanning provides high - density cross - section profile data (such as the convexity and concavity of the side wall and the local settlement of the vault), making up for the deficiency of the laser ranging that only monitors discrete points. The integration of the two improves the geometric parameter calibration accuracy by 40% (for example, the curve fitting error of the side wall is reduced from ±2 mm to ±1.2 mm). Triggering scanning based on the threshold of the inclination change amount realizes on - demand monitoring. While ensuring data integrity, the data acquisition frequency is reduced by 60% (the conventional scanning frequency is 1 time per hour, and the triggered scanning is on average 3 - 5 times per day).

[0023] Using the calibrated point cloud data, fit the spatial straight lines of each side of the quadrilateral by the least - squares method (such as the three - dimensional straight - line equation of the P1P2 side), calculate the actual side length (considering the influence of the tunnel surface, rather than the horizontal projection distance), and correct the initial length of the laser ranging. i y i z i ), calculate the actual interior angle between adjacent sides by vector cross - product, and eliminate the angle deviation in the initial measurement (such as the angle error caused by the instrument centering error).

[0024] Perform cross - section profile fitting (such as cubic spline curve) on the calibrated point cloud data, calculate the local curvature radii (R = 1 / k, where k is the curvature) of the vault (near P1), side wall (near P2 / P3), and invert arch (near P4), reflecting the change in the cross - section shape (such as the increase in curvature caused by the settlement of the vault).

[0025] Integrate the corrected side length (L i ’), interior angle (θ i ’), three - dimensional vertex coordinates (x i y i z i ), and cross - section curvature (R 1 - R 3 ) into a reference topological database, with each parameter attached with a timestamp (accurate to seconds), supporting subsequent spatio - temporal analysis of deformation.

[0026] Upgraded from two-dimensional plane geometry (initial measurement) to three-dimensional space topology, considering the true curved surface form of the tunnel section (such as the curvature characteristics of horseshoe-shaped and circular sections), avoiding the misjudgment of deformation caused by the plane assumption (for example, the length of the inclined side can differ from the horizontal projection by 3 - 5 mm / m in three-dimensional space). Introducing section curvature as an additional parameter can sensitively identify local buckling deformation (such as the sudden change in curvature caused by the concave inward of the side wall), forming a complement with side length and interior angle: The change in side length reflects "tension / compression", the change in interior angle reflects "shear dislocation", and the change in curvature reflects "bending instability". The combination of the three can comprehensively describe the deformation mode of the section (for example, extrusion deformation is manifested as an increase in curvature + a decrease in interior angle, and shear deformation is manifested as an increase in the difference in the lengths of opposite sides). After each calibration, a new reference topology is generated, automatically covering the old data, adapting to the stress adjustment of the surrounding rock during the tunnel construction stage (such as the stable convergence deformation period after the initial support), and ensuring the real-time matching of the reference parameters with the current state of the surrounding rock.

[0027] In a preferred embodiment of the present invention, based on the reference topology and the real-time distances between the vertices of the quadrilateral and the tunnel face, the weight of each vertex is calculated, including: Using a three-dimensional laser rangefinder to continuously measure the three-dimensional Euclidean distances between the vertices P1, P2, P3, and P4 of the quadrilateral and the tunnel face at a sampling frequency of 5 times per second, generating the original dynamic distance sequence of each vertex; performing sliding window average filtering on the original sequence with a window size of 10 sampling points to generate the smoothed dynamic distance time series data, including: Install a three-dimensional laser ranging base station at the center position of the tunnel face (tunnel excavation face), with a built-in high-precision inertial navigation system (positioning accuracy ±2 mm), emit laser beams to the four vertices (P1 - P4) of the quadrilateral, and calculate the three-dimensional Euclidean distances between the vertices and the tunnel face through the TOF (time of flight) principle.

[0028] Continuously collect data at a sampling frequency of 5 Hz to generate an original distance sequence containing timestamps (for example, the distance value of point P1 is recorded every 0.2 seconds to form the original distance sequence).

[0029] Filter the original sequence using a 10-point sliding average window (time span of 2 seconds): for each moment, calculate the average value of the first 5 points to the current point, and eliminate sudden noises (such as abnormal jump values caused by equipment vibration and dust reflection) to generate the smoothed distance time series data.

[0030] The 5Hz sampling rate can capture the rapid changes in the vertex distance during the tunneling face advancement (such as the sudden change in distance during blasting excavation). Compared with the traditional 1Hz sampling, the response speed to dynamic deformation is increased by 4 times. The sliding window filtering effectively suppresses high-frequency noise (such as the ±1mm random fluctuation caused by environmental vibration), reducing the standard deviation of the distance data from ±0.8mm to ±0.3mm. The distance data of all vertices are synchronously collected and filtered to ensure the alignment of timestamps (error < 10ms), avoiding the weight calculation deviation caused by asynchronous measurement.

[0031] Based on the attenuation characteristics of rock mass stress wave propagation, a negative exponential mapping relationship between vertex weight and real-time distance is established, and the stress attenuation curve is fitted to determine the weight function parameters. The smoothed dynamic distance time series data is input into the weight function parameters to calculate the weight values of each vertex in real time, including: Through on-site tests (such as burying stress sensors under similar surrounding rock conditions), measure the attenuation law of stress wave with distance caused by tunneling face blasting or excavation, and find that the stress amplitude and the distance conform to the negative exponential attenuation relationship: , where is the initial stress at the tunneling face, is the attenuation coefficient, which is determined by fitting the on-site data through the least squares method. For example, in mudstone , and in sandstone .

[0032] Based on the principle that "the closer to the tunneling face, the greater the disturbance of the vertex by excavation and the higher the weight should be", the vertex weight is defined as being positively correlated with the stress amplitude, that is, a negative exponential weight function is constructed: , where and are normalization constants, and the range is [0, 1]. For example, when , , and when , .

[0033] In tunnel sections with different surrounding rock grades (such as grade I to grade V), by comparing the correlation between vertex deformation and weight value (such as using Pearson correlation coefficient), optimize the attenuation coefficient and the normalization constants (A, B). For example, set a larger k value in the fractured surrounding rock (grade V) to make the weight of the close-range vertices increase faster and highlight its sensitivity.

[0034] The weight function directly maps the propagation law of rock mass excavation disturbance, solves the problem that the traditional weight setting (such as empirical assignment) lacks a theoretical basis, and makes the weight value positively correlated with the actual deformation risk (for example, the weight of the vertex 1m away from the tunnel face is 0.9, which drops to 0.4 at 3m, and approaches 0 outside 5m); by calibrating parameters on site, the weight attenuation rate is dynamically adjusted in hard rock (slow attenuation) and soft rock (fast attenuation) to ensure the universality of the weight model. The non-linear characteristic (exponential attenuation) of the weight function amplifies the importance of the vertices at close range. For example, when the tunnel face advances 2m in front of point P2, the weight of P2 suddenly rises from 0.6 to 0.85, triggering high-frequency monitoring of the sides where P2 is located (such as P1P2, P2P3) (the acquisition frequency of the change rate is increased from 1 time per second to 5 times per second), realizing the intelligent allocation of monitoring resources.

[0035] The smoothed distance data , (i = 1, 2, 3, 4 corresponding to the four vertices) are input into the calibrated weight function in real time to calculate the current weight value of each vertex: ; For different vertices, due to geometric position differences, independently calibrated parameters may be adopted. For example, for the crown P1 and the invert P4, different k values may be set due to different stress wave propagation paths.

[0036] represents the weight value of the i-th vertex (i = 1, 2, 3, 4 corresponding to P1 to P4) at time t, reflecting the sensitivity of this vertex to the excavation of the tunnel face (the greater the weight, the more significant the vertex is affected by stress disturbance and belongs to the key deformation area).

[0037] represents the smoothed three-dimensional space Euclidean distance (unit: meter or millimeter) between vertex i and the tunnel face at time t, which is obtained by performing moving window (window size 10) average filtering on the original distance sequence (sampled 5 times per second) to eliminate high-frequency noise (such as instrument vibration, environmental interference). The closer to the tunnel face ( the smaller), the greater the weight , reflecting the characteristic that "the vertices near the tunnel face are more directly affected by excavation stress disturbance" (in line with the law of rock mass stress wave attenuation with distance).

[0038] represents the amplitude parameter of the weight function of vertex i, which determines the maximum possible value of the weight (when is close to 0, the value of approaches +B i ), and different vertices may set different A i due to geometric position differences (such as different stress propagation paths from the tunnel face for the crown P1 and the invert P4).

[0039] k i represents the stress attenuation coefficient of vertex i, reflecting the rate at which the weight decays as the distance increases (the larger k i , the faster the weight decreases as the distance increases). For example, for the sidewall vertices (P2, P3), since they are close to the free face of the tunnel face excavation, k i can be set to a relatively large value to make its weight more sensitive to distance changes.

[0040] B i represents the offset of the weight function of vertex i, which is used to adjust the reference value of the weight (for example, when the distance is extremely far, the weight approaches B i , to avoid missing the judgment of key areas due to the weight becoming zero).

[0041] Weight threshold linkage mechanism: Preset a weight threshold (such as 0.7). When the weight of any vertex > 0.7, the system automatically triggers real-time monitoring of the change rate of the adjacent edges of this vertex (for example, the adjacent edges of P2 are P1P2, P2P3, P2P4). Calculate the relative change amount of the side length within a unit time (such as △L / L0×100% / min), and synchronize the monitoring data to the edge computing node for real-time analysis. △L represents the absolute change amount of the side length within a unit time (such as the number of millimeters increased / decreased by a certain side within 1 minute); L0 represents the reference length of this side (the initial side length from the reference topology, used as a reference benchmark for the deformation amount). The relative change rate (△L / L0) can better reflect the deformation characteristics of the slope material than the absolute change amount (for example, the allowable strain thresholds of slopes with different materials are different). Combining the time dimension (per minute) can track the deformation rate in real time to determine whether to trigger an early warning.

[0042] Update the weight value once per second to reflect in real time the distance changes caused by the advancement of the tunnel face or the displacement of the surrounding rock (for example, when the tunnel face advances 5m per day, the vertex weight decays dynamically according to an exponential law), avoiding the defect that a fixed weight cannot adapt to working condition changes; the threshold trigger mechanism realizes "monitoring on demand": only when the vertex enters the high-weight area (close to the tunnel face), start high-frequency monitoring of the edge change rate. Compared with full-time high-frequency monitoring, it can reduce the sensor energy consumption and data transmission pressure by 40%. By focusing the weight on the vertices near the free face (for example, the weights of the sidewall vertices P2 / P3 within 3m behind the tunnel face are often higher than 0.7), the monitoring sensitivity of the system to the change rate of the side length in this area is increased by 3 times (the sampling frequency is increased from 1 time / minute to 1 time / second). For example, in a certain tunnel project, it successfully identified 2 hours in advance the sudden increase in the weight of the sidewall vertex caused by the tunnel face blasting (from 0.6 to 0.85), and captured the abnormal tensile change rate of the adjacent edge P2P3 (0.5% / min), avoiding the risk of missing the judgment of local landslides.

[0043] In a preferred embodiment of the present invention, when the weight of a vertex is greater than a threshold, the change rate of the corresponding edge is monitored in real time; temperature compensation is performed on the change rate, and the deformation caused by the temperature change is calculated by analyzing the thermal expansion characteristics of the material, including: Normalize the weight value of each vertex to the range of 0-1. When the vertex weight exceeds the preset deformation sensitivity threshold, mark the vertex as a high-risk node and activate the real-time change rate monitoring module of the edges connected to it. At the same time, bind the weight data to the spatial coordinates in the baseline topology to generate a weighted dynamic topology node information table, including: The real-time weight values ​​of each vertex are uniformly scaled to the range of 0-1 through linear transformation to eliminate the problem of incomparable weight values ​​of different vertices due to parameter differences (for example, the initial weight range of P2 is 0.3-0.9, and that of P4 is 0.1-0.6, and both are in the range of 0-1 after normalization). A deformation sensitivity threshold (such as 0.7) is preset. When the normalized weight of a vertex exceeds the threshold, it is automatically marked as a high-risk node (for example, the side wall vertices P2 / P3 within 3m behind the face often trigger this mark), and the real-time monitoring module of all the edges connected to it (for example, P2 connects the three edges of P1P2, P2P3, and P2P4) is activated.

[0044] Dynamic topology information binding: The weight value of the high-risk node is bound to its three-dimensional spatial coordinates (x, y, z) in the baseline topology and the edge information (such as connecting vertices, baseline edge length, initial internal angle) to generate a weighted dynamic topology node information table (for example, the table fields include: node ID, weight value, x coordinate, y coordinate, connecting edge 1-ID, connecting edge 2-ID, etc.), providing a spatial reference for subsequent edge positioning and deformation analysis.

[0045] Through the weight threshold, the vertices that are significantly disturbed by excavation (such as nodes near the free surface) are screened out, and high-frequency monitoring is started only for the edges connected to them. Compared with indiscriminate monitoring of all edges, this can reduce 60% of invalid data collection and improve the monitoring efficiency of key areas. The dynamic topology table associates the abstract weight value with the specific spatial location.

[0046] According to the marked high-risk vertices, locate the quadrilateral edges connected to them; start the high-frequency monitoring mode of the three-dimensional laser rangefinder for the target edge, continuously collect the original data of the edge length at a sampling frequency of 1Hz, and use the difference calculation between two adjacent sampling points to generate the original length change rate sequence; remove outliers from the original change rate data, and the removal rule is: if the single point change rate exceeds 3 times the standard deviation of the historical mean value during the same period, it is marked as noise and linear interpolation repair is performed to obtain the processed change rate data, including: For the edges connected to high-risk nodes (such as the P1P2 edge connected to P2), start the 1Hz high-frequency monitoring mode of the 3D laser rangefinder (the conventional mode is once per minute), and continuously collect the original data of the side length (such as recording the side length values L(t1), L(t2),... once per second).

[0047] Calculate the original length change rate by the adjacent sampling point difference method: For example, the side length difference between the current moment t and the previous moment t-1 is △L = L(t) - L(t-1), and the change rate is expressed as △L / L0×100% / second (L0 is the reference side length of this edge), generating an original change rate sequence (such as deformation rate data including positive and negative values).

[0048] Outlier rejection and data repair: Establish a 3-fold standard deviation rejection rule: Calculate the mean μ and standard deviation σ of the change rate in the historical same period (such as the past 1 hour). If the change rate of a single point exceeds μ±3σ, it is determined as noise (such as jump values caused by short-term occlusion of the device or sudden vibration).

[0049] Repair the noise points by the linear interpolation method (such as replacing with the average value of the previous valid point and the next valid point). For example, if the 5th sampling point is an outlier, fill it with the average value of the change rates of the 4th and 6th points (△L4+△L6) / 2 to obtain the smoothed change rate data.

[0050] The 1Hz sampling frequency can capture rapid deformations caused by instantaneous loads such as blasting excavation and rainfall infiltration (such as a sudden change in side length of 0.5mm within 1 second). Compared with the conventional low-frequency monitoring, the response speed to sudden deformations is increased by 60 times, avoiding missed judgments of short-duration high-risk deformations. By rejecting outliers through statistical rules (such as ±5mm jumps caused by dust interference), the standard deviation of the change rate data is reduced from ±0.8% / min to ±0.2% / min, providing clean input data for subsequent temperature compensation and avoiding misjudgments caused by noise.

[0051] Deploy a distributed temperature sensor array on the surface and inside of the rock mass on both sides of the target edge, collect the temperature field data of the entire section of the monitored edge, and generate a temperature spatio-temporal distribution matrix; call the pre-stored rock mass material thermal expansion coefficient library, and match the corresponding thermal expansion coefficient according to the material type of the current monitored edge; calculate the predicted deformation amount caused by temperature fluctuations point by point based on the thermal expansion coefficient, the original reference side length, and the real-time temperature difference, where the real-time temperature difference is the difference between the current temperature and the initial temperature at the time of generating the reference topology, including: Optical fiber temperature sensors or thermocouples are arranged at intervals of 0.5 - 1 m on the surface (such as the sidewall concrete support layer) and inside (such as buried in boreholes) of the rock mass on both sides of the target edge, forming a temperature field acquisition network covering the entire length of the monitored edge (for example, 11 sensors are arranged on a 5 m long edge, 1 at each end, and 1 every 0.5 m in the middle), and temperature data is collected in real time (accuracy ±0.5 °C), generating a temperature spatio-temporal distribution matrix containing position coordinates and timestamps (such as a three-dimensional array: [edge position, time, temperature]).

[0052] Call the pre-stored thermal expansion coefficient library of rock mass materials (including thermal expansion coefficients of various materials such as sandstone, mudstone, and concrete, such as sandstone and concrete), and match the corresponding thermal expansion coefficient α according to the material type of the area where the monitored edge is located (such as determined through geological exploration reports or on-site sampling).

[0053] Calculate the real-time temperature difference △T = T_current - T_reference (T_reference is the initial temperature when the reference topology is generated, usually taking the ambient temperature at the start of monitoring). Based on the thermal expansion formula △L_thermal = α × L0 × △T, calculate the temperature-induced deformation amount at the corresponding position of each sensor point by point (for example, when the temperature difference at a certain point is 10 °C and the reference side length is 10 m, the thermal deformation is 10 × 10 −6 ×10000 mm × 10 = 1 mm).

[0054] The distributed sensor array covers the surface and inside of the monitored edge, solving the problem that single-point temperature measurement cannot reflect the overall temperature difference of the edge (such as when the temperature difference at both ends of the edge is 3 °C, single-point measurement will cause a 20% compensation error), improving the temperature field resolution to 0.5 m / point and the compensation accuracy by 30%. Through material-specific thermal expansion coefficient matching (instead of a unified default value), precise compensation is carried out for different rock mass / support materials (such as the thermal expansion difference between steel supports and concrete linings), improving the side length monitoring accuracy from ±2 mm to ±0.5 mm (for example, when the day-night temperature difference is 15 °C, the thermal deformation of a 10 m steel string sensor is about 1.74 mm, and this error can be completely eliminated through compensation), ensuring that the rate-of-change data only reflects the true mechanical deformation rather than temperature pseudo-variation.

[0055] In a preferred embodiment of the present invention, the monitored side length is corrected to generate a temperature-compensated side length, including:[[]] Extract the sequential temperature data collected by the temperature sensors and the corresponding theoretical temperature deformation amounts; Based on the thermal expansion characteristics of the rock mass material, perform a sliding time window integration on the theoretical temperature deformation amounts in time series, accumulate the instantaneous deformation amounts within each time slice, and generate a temperature compensation cumulative value curve; Obtain the original side length data measured by the three-dimensional laser rangefinder in real time, and align it with the temperature compensation cumulative value according to the timestamp; Perform compensation superposition calculation: Temperature-compensated side length = Original side length + Cumulative compensation value, where the sign of the compensation value is determined according to the temperature rise and fall direction. When the temperature rises and causes expansion, the compensation value is negative, and vice versa.

[0056] In the embodiments of the present invention, by accumulating the theoretical deformation (such as thermal expansion / contraction) caused by temperature changes, the interference of ambient temperature on side length monitoring is removed in real time, and the side length monitoring accuracy is improved from ±2 mm to ±0.5 mm. For example, in a scenario with a 20°C day-night temperature difference, the thermal deformation of a 10-m long steel sensor can reach 2.32 mm (steel thermal expansion coefficient 11.6×10 -6 / °C). Through compensation, this error can be completely eliminated, ensuring that the monitoring data only reflects the true mechanical deformation of the surrounding rock (such as joint dislocation, rock mass creep, etc.). By using a sliding time window integration (such as accumulating instantaneous deformation in 1-minute time slices), the continuous change of temperature over time (such as non-linear expansion during the heating process) can be captured, avoiding the defect of traditional single-point compensation (only calculating with the current temperature difference) that ignores the temperature change process. For example, in a scenario with a continuous temperature rise for 3 hours, the expansion amount per minute and second can be accurately accumulated to generate a continuous compensation value curve, making the compensation effect more consistent with the actual deformation process. By aligning the original side length data with the temperature compensation cumulative value through timestamps, it is ensured that the two strictly correspond in the space-time dimension (error < 10 ms), solving the compensation deviation problem caused by asynchronous acquisition of different sensors (laser rangefinder and temperature sensor). For example, when the laser rangefinder collects the side length at 10:00:05, the temperature cumulative compensation value at the same moment can be accurately matched, avoiding over-compensation or under-compensation caused by time misalignment (such as mistakenly using the compensation value at 10:00:00 for the data at 10:00:05, which may introduce an error of ±0.2 mm). When the temperature rises, the rock mass / sensor expands, and the measured side length will "falsely become longer". The true mechanical deformation needs to be restored by subtracting the expansion amount (the compensation value is negative); when the temperature drops, the contraction causes the measured side length to "falsely become shorter", and it needs to be corrected by adding the contraction amount (the compensation value is positive), avoiding the reverse correction caused by confusion in the compensation direction (such as mistakenly setting the temperature rise compensation value as positive), ensuring that the corrected side length completely eliminates temperature pseudo-deformation.

[0057] The side length data after temperature compensation is used as an intermediate result and directly input into the surrounding rock rheological equation (such as the Burgers model) for aging deformation correction, avoiding the interference of temperature noise on creep analysis (such as when not compensated by traditional methods, thermal expansion may be misjudged as plastic deformation of the rock mass, resulting in a calculation deviation of rheological parameters exceeding 40%). At the same time, the pure side length data improves the calculation accuracy of three-dimensional deformation indexes (such as diagonal displacement difference, adjacent side curvature ratio), provides high-quality input for the spatio-temporal attention GNN model, and finally improves the early warning accuracy from 65% to 92%.

[0058] In a preferred embodiment of the present invention, according to the temperature-compensated side lengths and the surrounding rock rheological equation, the lengths of two adjacent sides are corrected by the strain difference between the adjacent sides to obtain the final side lengths; the final side lengths are combined with the reference interior angles to generate a topological structure, including: The axial strain distribution data of two adjacent sides are collected in real time by using a distributed optical fiber strain sensor array, and the average strain difference between the two sides is calculated based on the collected strain distribution data, including: on the rock mass surface or within the support structure of two adjacent sides of a quadrilateral (such as P1P2 and P2P3), distributed optical fiber strain sensors (such as BOTDA / BOTDR) are arranged at intervals of 0.2 - 0.5 m along the axial direction, and the axial strain data of each measurement point are collected in real time (accuracy ±5 με) to form a strain distribution curve (such as the strain sequence of side P1P2 is ε1, ε2,..., εn); the arithmetic mean of the strain data of each side is taken to obtain the average strain value of the side (such as εg1 corresponds to side P1P2, and εg2 corresponds to side P2P3), and the difference △ε = εg1 - εg2 is calculated to reflect the deformation difference between adjacent sides caused by rock mass rheology.

[0059] Based on the average strain difference, the aging deformation parameters in the surrounding rock rheological equation, including the viscoelastic coefficient and the creep rate, are called, and in combination with the time difference between the current monitoring time and the reference topology generation time, the co-deformation ratio of adjacent sides caused by the rheological effect is calculated, including: the viscoelastic coefficient of the current stratum (such as the spring stiffness k1, k2 and the viscosity η1, η2 of the dashpots in the Burgers model) and the creep rate β (reflecting the strain growth trend per unit time, calibrated by in-situ creep tests, such as β = 0.01% / h for mudstone) are retrieved from a pre-established surrounding rock rheological parameter library.

[0060] The time difference △t (accurate to hours) between the current monitoring time and the reference topology generation time is calculated, and based on the time dependence of the rheological equation, the co-deformation ratio γ of adjacent sides under the rheological effect is deduced, such as , indicating the proportion of co-deformation of adjacent sides caused by rheology with the increase of time. When △t = 0, γ = 0 (no time effect, co-deformation is 0); as △t increases, γ approaches 1 (long-term creep causes the deformation of adjacent sides to tend to be coordinated), which conforms to the law of progressive development of viscoelastic deformation with time in the Burgers model.

[0061] According to the co-deformation ratio of adjacent sides, the change in the length of adjacent sides caused by the rheological effect is deduced, including: According to the cooperative deformation ratio γ and the average strain difference △ε, the additional length change △Lrheology of each edge due to the rheological effect is calculated. For example, if the reference length of the adjacent edge P1P2 is L0, its rheological deformation is △Lrheology1=γ×△ε×L0, reflecting the proportional cooperative change of the edge length caused by the long-term creep of the rock mass (if the strain difference is positive, it means that the tensile deformation of the P1P2 edge is greater than that of the P2P3 edge, and the rheological correction needs to adjust the difference between the two).

[0062] The temperature compensated side length is taken as input, combined with the change in the length of the adjacent side, and the rheological correction model is used to reversely correct the length of the two adjacent sides to obtain the corrected side length, including: taking the temperature compensated side length Ltemperature compensation as input, combined with the length change △Lrheology caused by rheology, and adjusting the side length through reverse correction logic: if the rheology causes the side A to lengthen and the side B to shorten, then △Lrheology is subtracted from the side A and △Lrheology is added to the side B (or vice versa), so that the deformation difference of the adjacent sides conforms to the synergistic characteristics of the rock rheology, and the corrected final side length L is obtained (such as Lfinal1=Ltemperature compensation1-△Lrheology1, Lfinal2=Ltemperature compensation2+△Lrheology2).

[0063] Based on the corrected side length and reference internal angle data, with P1 as the coordinate origin, the spatial coordinates of P2, P3, and P4 are calculated vertex by vertex; the curvature of the quadrilateral section is fitted using the least squares method, and the curvature radius parameters are updated; the corrected side length, updated internal angle, and curvature radius parameters are integrated to generate a dynamic topological structure containing geometric parameters and curvature information, including: Taking the vault settlement point P1 in the reference topology as the coordinate origin (0, 0, 0), based on the corrected side lengths (L final 1, L final 2, etc.) and reference internal angles (such as ∠P1P2P3), the three-dimensional space coordinates of P2, P3, and P4 are calculated vertex by vertex through the geometric analysis method.

[0064] For example, given the side lengths of P1P2 and P1P3, as well as the interior angle of ∠P1, the coordinates of point P2 (x2, y2, 0) are calculated using the cosine theorem, and then the z coordinate of point P4 (reflecting the uplift of the invert) is determined through three-dimensional distance constraints.

[0065] The curvature fitting and parameter updating process is as follows: The coordinates of each vertex of the quadrilateral and the cross-section point cloud data (from the point cloud scan of the benchmark topology) are imported into the fitting algorithm, and the least squares method is used to fit the quadratic curve. The curvature radius R of the vault, side wall, and invert is calculated (for example, R decreases due to the sinking of the vault, reflecting the aggravation of the bending deformation), and the curvature parameters are updated (for example, R of the vault changes from 10m to 9.5m, indicating that the local bending degree has increased).

[0066] Dynamic topology integration is implemented as follows: Integrate the corrected side lengths (L_final1 - L_final6), updated interior angles (considering the small angle changes caused by rheology), vertex coordinates (x, y, z), and radius of curvature (R_crown, R_sidewall, R_invert) into a dynamic topological structure, with each parameter attached with a timestamp and weight information (e.g., high-risk edges are marked red and low-risk ones are blue) to form a three-dimensional geometric model that can be updated in real time.

[0067] In the embodiments of the present invention, the traditional method only focuses on immediate deformation and ignores time-dependent effects such as rock mass creep and relaxation (e.g., the creep strain of mudstone can reach 0.3% within 30 days). By introducing the time dimension through the rheological equation, the long-term deformation trend can be accurately captured, avoiding misjudging creep as a stable state. The internal joint dislocation or interlayer slip of the rock mass is identified by the strain difference between adjacent sides (e.g., the strain difference is enlarged due to the different creep rates of sandstone and mudstone in interbedded rock masses). The corrected side lengths are closer to the true mechanical deformation, enabling the progressive instability process of the slope to be quantitatively tracked (e.g., in a tunnel case, the accelerated deformation caused by mudstone softening was discovered 15 days in advance through strain difference correction). Integrate geometric deformations (side lengths, interior angles), spatial positions (coordinates), and morphological features (curvature) to construct a full-element topological model including "points - lines - planes". Compared with traditional single-point monitoring, it can describe the complete mode of "areal deformation" of the slope (e.g., the abnormal curvature combination of crown settlement accompanied by inward concavity of the sidewall). The dynamic update of parameters such as the radius of curvature provides real-time data for subsequent calculation of deformation indices (such as the curvature ratio of adjacent sides), supporting the spatio-temporal attention GNN model to capture non-linear deformation features (e.g., the continuous decreasing trend before the curvature mutation), increasing the early warning accuracy by 27% (from 65% to 92%). First, eliminate environmental noise through temperature compensation, and then eliminate the interference of time-dependent deformation through rheological correction to form a closed loop of "data purification - physical mechanism modeling - dynamic correction", improving the side length monitoring accuracy from ±2mm to ±0.5mm and being applicable to long-term monitoring.

[0068] In a preferred embodiment of the present invention, construct a three-dimensional deformation index matrix including diagonal displacement difference, adjacent side curvature ratio, and interior angle coefficient of variation, including: Extract the current lengths of the two diagonals of the quadrilateral from the dynamic topological structure and compare them with the initial diagonal lengths in the reference topology; calculate the absolute displacement difference and the displacement difference change rate of each diagonal to generate a displacement difference feature vector of the two diagonals; Based on the updated cross-section curvature parameters in the dynamic topology, obtain the radius of curvature and the curvature ratio of adjacent sides, normalize the curvature ratio, map the ratio to the 0 - 1 interval, and generate a curvature ratio feature vector; Read the current angle values of the four interior angles of the quadrilateral from the dynamic topology, calculate the interior angle standard deviation and the interior angle coefficient of variation, and perform a sliding window statistic on the coefficient of variation to generate a smoothed coefficient of variation sequence; Align the diagonal displacement difference eigenvector, curvature ratio eigenvector, and interior angle coefficient of variation sequence according to the timestamp; construct a three-dimensional deformation index matrix; Perform data standardization processing on the three-dimensional deformation index matrix to generate a standardized deformation index matrix; the dimension of the three-dimensional deformation index matrix is [timestamp × deformation feature], where each row contains: the first column: the rate of change of the displacement difference between two diagonals; the second column: the normalized value of the curvature ratio of adjacent sides; the third column: the coefficient of variation of the interior angle.

[0069] In the embodiment of the present invention, the implementation process of constructing the three-dimensional deformation index matrix: Obtain the real-time lengths of the two diagonals (such as P1P3, P2P4) of the current quadrilateral from the dynamic topological structure, and at the same time retrieve the initial lengths (reference values) of these two diagonals in the reference topology; calculate the absolute displacement difference of each diagonal (the current length minus the reference length), and further calculate the rate of change of the displacement difference per unit time (such as the change amount of the displacement difference per hour), forming a two-dimensional vector containing the rate of change of the two diagonals (for example, [0.3mm / h, -0.1mm / h]), as the diagonal displacement difference eigenvector. The diagonal connects the opposite vertices of the quadrilateral, and its length change sensitively reflects the in-plane shear deformation (such as joint dislocation causing one diagonal to elongate and the other to shorten), and the rate of change captures the deformation acceleration trend.

[0070] Extract the cross-sectional curvature radii (such as the vault curvature R1 and the sidewall curvature R2) corresponding to two adjacent sides (such as P1P2 and P2P3) in the dynamic topology, and calculate the ratio of the two (R1 / R2 or R2 / R1, determined according to the position of the side). Use a normalization method (such as Min-Max normalization) to map the curvature ratio to the 0-1 interval (for example, the original ratio of 5.0 is mapped to 1.0, and the ratio of 1.0 is mapped to 0.5), eliminate the influence of dimensions, and generate a one-dimensional curvature ratio eigenvector (such as [0.8] indicating that the current curvature ratio is in the high-value interval). The curvature ratio reflects the difference in the bending degree of adjacent regions, and after normalization, it is convenient for horizontal comparison of different cross-sections and different times to identify local convex and concave deformations (such as the sudden increase in the curvature ratio caused by the inward concavity of the sidewall).

[0071] Read the real-time angle values of the four interior angles (such as ∠P1, ∠P2, ∠P3, ∠P4) of the quadrilateral in the dynamic topology, calculate the standard deviation of these four angles (reflecting the degree of angle dispersion), and then divide by the average angle to obtain the coefficient of variation of the interior angle (dimensionless, characterizing the uniformity of the angle distribution). Apply a sliding window smoothing process (such as a window size of 10 time points) to the coefficient of variation sequence to eliminate sudden noise (such as angle jumps caused by instrument jitter), and generate a smoothed coefficient of variation sequence (such as outputting a smoothed value every 10 minutes). The coefficient of variation of the interior angle sensitively captures angle anomalies caused by joint dislocation or rock mass fragmentation (such as a sudden decrease of 2° in a certain interior angle, resulting in a sudden increase in the coefficient of variation), and the smoothing process improves data stability.

[0072] Strictly align the diagonal displacement difference eigenvector (2D), curvature ratio eigenvector (1D), and internal angle coefficient of variation sequence (1D) according to the timestamp (error < 1 second) to ensure that multi-feature data at the same moment corresponds to the same deformation state. Construct a three-dimensional matrix in the dimension of "timestamp × deformation feature", and each row contains three core indicators: The change rate of the displacement difference between two diagonals (2 columns); The normalized curvature ratio (1 column); The smoothed internal angle coefficient of variation (1 column).

[0073] Adopt the Z-Score standardization method to zero the mean and normalize the standard deviation of the data in each column of the matrix, eliminating the dimension difference (such as the unification of the displacement difference unit mm / h and the dimensionless coefficient of variation), and generating a standardized deformation index matrix. Time alignment ensures the spatio-temporal consistency of multi-features, and standardization provides a unified input format for the subsequent neural network model, improving the training efficiency and prediction accuracy.

[0074] The diagonal displacement difference directly reflects the in-plane shear deformation (such as the increase in the diagonal length difference caused by the internal dislocation of the landslide body). Compared with single-point displacement monitoring, it can capture the shear characteristics of "one pair of vertices approaching and the other pair moving away", and the identification efficiency is increased by 60%. Quantify the difference in the bending degree of adjacent sides (such as when the crown subsides and the side walls concave inward, the curvature ratio decreases significantly), and it is sensitive to local buckling instability (such as the abnormal curvature before the shotcrete cracks), making up for the defect that traditional monitoring does not pay attention to "shape changes". Through the change in angular uniformity (such as the increase in the coefficient of variation caused by the mutation of an internal angle), the sliding of the joint surface or the fragmentation of the rock mass can be identified early (such as the increase in the dispersion of internal angles caused by blasting vibration), and the sensitivity reaches 0.5°.

[0075] Normalization and smoothing processing can eliminate the influence of dimensions (such as the unification of the curvature radius unit m and the displacement difference unit mm), filter out high-frequency noise (such as accidental fluctuations caused by environmental vibration), and make the index sequence more conform to the true deformation trend (such as the gradual increase in the coefficient of variation during the creep of mudstone can be clearly captured); time alignment and standardization can provide high-quality input for the spatio-temporal attention GNN model (avoiding mis-correlation of features caused by time misalignment), enabling the model to effectively learn the spatio-temporal correlation of deformation indicators (such as the continuous 3-hour decrease in the curvature ratio and the increase in the internal angle coefficient of variation at a certain section indicating potential instability).

[0076] In a preferred embodiment of the present invention, input the index matrix, the spatio-temporal matrix of the support structure strain, and the rock mass fragmentation tensor into the spatio-temporal attention graph neural network, and output the deformation evolution trend including the diagonal displacement difference, the curvature ratio of adjacent sides, and the internal angle coefficient of variation, including: Align the standardized deformation index matrix, the support structure strain spatio-temporal matrix constructed from the strain time-series data of the support structure collected by fiber optic sensors, and the rock mass fragmentation tensor calculated based on the microseismic event energy distribution and rock mass fracture density according to the time stamp to obtain time-aligned data; Based on the time-aligned data, perform dimensionality reduction on the rock mass fragmentation tensor, extract its principal component feature vectors, and make its column dimension match the deformation index matrix; Using the principal component feature vectors, together with the deformation index matrix and the support structure strain spatio-temporal matrix, splice the three types of data features along the time axis to construct a fused spatio-temporal tensor; According to the fused spatio-temporal tensor, use the multi-head attention mechanism to calculate the spatio-temporal dependence relationships between the three types of data respectively; According to the spatio-temporal dependence relationships, multiply the fused spatio-temporal tensor by the attention weight matrix, and aggregate the spatial correlation features of the quadrilateral vertices and their support nodes through a graph convolutional network; Based on the spatial correlation features, use a gated recurrent unit to learn the long-term evolution trend of the deformation index; Input the long-term evolution trend into a fully connected prediction module, and output the deformation evolution trend parameters within the next 3 hours, including: the increment of the diagonal displacement difference change rate, the change gradient of the adjacent side curvature ratio, and the acceleration of the interior angle coefficient of variation.

[0077] In the embodiment of the present invention, the standardized deformation index matrix generated from dynamic topology (including diagonal displacement difference, curvature ratio, and interior angle coefficient of variation); the support structure strain spatio-temporal matrix collected in real time by fiber Bragg grating sensors (such as the strain time-series data of anchor cables and bolts, with a resolution of 1% and the time stamp accurate to seconds); the rock mass fragmentation tensor based on the microseismic monitoring system (calculated through microseismic event energy distribution and fracture density, reflecting the integrity of the rock mass, with an initial dimension of 30×30×30), match the time stamps of the three types of data with millisecond-level accuracy, and eliminate abnormal data points with a time misalignment exceeding 10 ms to ensure the one-to-one correspondence of the deformation index, support response, and rock mass state at the same moment.

[0078] Apply principal component analysis (PCA) to the high-dimensional rock mass fragmentation tensor (such as 30×30×30), calculate its covariance matrix and extract the first 3 principal components (cumulative variance contribution rate > 95%), and compress the tensor into a 3-dimensional principal component feature vector (such as [PC1, PC2, PC3]). Adjust the dimension of the feature vector to make its number of columns consistent with the deformation index matrix (3 columns: displacement difference change rate, curvature ratio, coefficient of variation) for subsequent splicing (such as expanding to 3 columns: PC1 corresponding to displacement difference, PC2 corresponding to curvature ratio, PC3 corresponding to coefficient of variation), which can reduce the data complexity, retain the main features of the rock mass fragmentation (such as fracture development direction, fracture zone distribution), and avoid overfitting of model training caused by high-dimensional data.

[0079] Concatenate the three types of data column by column along the time axis: Columns 1 - 3: Normalized deformation indices (rate of change of displacement difference, normalized value of curvature ratio, coefficient of variation of interior angle); Columns 4 - 6: Strain parameters of the support structure (such as rate of change of anchor cable tension, growth rate of bolt axial force); Columns 7 - 9: Principal component features of the rock mass fragmentation degree after dimensionality reduction.

[0080] Form a fused spatio - temporal tensor (dimension: timestamp × 9 features). For example, the data at a certain moment is [0.5mm / h, 0.8, 0.25, 1.2%, -0.5%, 0.8%, 0.6, -0.3, 0.9], which contains three types of features: deformation, support, and rock mass state.

[0081] Convert heterogeneous data (geometric deformation, structural response, rock mass properties) into a feature space of a unified dimension.

[0082] Apply a multi - head attention layer (such as 8 heads) to the fused spatio - temporal tensor. Each head independently calculates the attention weights between the three types of data: Time dimension: Analyze the dependence relationship between the features at the current moment and the features in the past 12 hours (such as the correlation between the current curvature ratio and the displacement difference 3 hours ago); Spatial dimension: Calculate the spatial correlation between the vertices of the quadrilateral (such as the correlation between the strain at point P2 and the fragmentation degree at point P3).

[0083] Generate an attention weight matrix to quantify the interaction intensity of different features, different moments, and different spatial positions (such as the higher the weight value, the closer the correlation).

[0084] The graph convolutional network aggregates spatially correlated features. The implementation process is as follows: Construct a graph structure: Use the 4 vertices of the quadrilateral and 2 support nodes (such as the anchor cable anchorage point, bolt end point) as nodes, and use the edges (vertex connection relationship, support - rock mass contact relationship) as edges. Define node features (coordinates, weights) and edge features (edge length, strain difference). Through graph convolutional operations (such as GCN layers), each node aggregates the features of adjacent nodes (such as node P2 fuses the strain and fragmentation degree information of P1, P3, and the support nodes), generating high - order features containing spatial correlations (such as a comprehensive index reflecting the coordinated deformation of the vertex - support). Utilize the geometric topology structure of the slope (quadrilateral element) to explicitly model the spatial proximity effect (such as the deformation of adjacent sides will affect each other), and solve the defect that traditional neural networks ignore spatial position relationships.

[0085] The gated recurrent unit learns the long - term evolution trend. The implementation process is as follows: Input the spatial correlation features output by the graph convolution into the GRU layer, and selectively memorize historical information through the gating mechanism: update gate, which determines how much of the past deformation trend to retain (e.g., a decreasing trend in the curvature ratio lasting for 12 hours will be strongly memorized); reset gate, which forgets short-term noise (e.g., fluctuations in the coefficient of variation caused by accidental vibrations); output long-term evolution feature vectors to capture the time dependence of deformation indicators (e.g., the rate of change of the displacement difference caused by mudstone creep gradually increases over time). By processing non-stationary time series (such as the non-linear trend in the deformation acceleration stage), compared with traditional LSTM, it captures long-term dependencies with fewer parameters and improves the training efficiency.

[0086] Input the long-term features output by the GRU into the fully connected prediction module, map them to the target space through a multi-layer perceptron (MLP), and output the deformation evolution trend parameters for the next 3 hours: the increment of the rate of change of the diagonal displacement difference (e.g., it is expected to increase by 0.2 mm / h in the next 1 hour); the change gradient of the curvature ratio of adjacent sides (e.g., the rate of decrease of the curvature ratio accelerates by 0.05 per hour); the acceleration of the coefficient of variation of the interior angle (e.g., the second derivative of the coefficient of variation, reflecting the acceleration degree of angle anomalies).

[0087] Deformation indicators (such as the curvature ratio) reflect geometric deformation, support strain (such as anchor cable tension) reflects structural response, and fragmentation tensor (such as fracture density) characterizes the properties of rock masses. The joint input of these three types of data avoids the one-sidedness of a single indicator (e.g., only looking at displacement may miss the hidden risk caused by the failure of the support). The multi-head attention mechanism identifies the temporal causal relationship of deformation indicators (e.g., when the rate of change of the displacement difference continuously rises for 3 hours, the curvature ratio usually becomes abnormal accordingly). Compared with traditional statistical models, the ability to capture non-linear trends is improved by 40%. The graph convolution explicitly utilizes the quadrilateral topology of the slope to quantify the spatial synergy effect between vertices (e.g., the deformation at point P2 will affect the curvature at point P3 through the edge P2P3), solves the defect of "unstructured input data" in traditional neural networks, and improves the utilization rate of spatially related features by 60%.

[0088] In a preferred embodiment of the present invention, multi-level early warning decisions are executed according to the component parameters of the deformation evolution trend, including: Extract the increment of the rate of change of the diagonal displacement difference, the change gradient of the curvature ratio of adjacent sides, and the acceleration of the coefficient of variation of the interior angle in the deformation evolution trend; based on the statistical correlation analysis of historical deformation data and support structure failure cases, set the grading threshold intervals for the three parameters, including the yellow early warning threshold interval, the orange early warning threshold interval, and the red early warning threshold interval; When any parameter first enters the yellow early warning threshold interval, increase the three-dimensional laser ranging frequency of the corresponding monitored side; associate the early warning information with the coordinates of high-risk nodes in the reference topology to generate a primary early warning signal with spatial positioning; When two or more parameters simultaneously enter the orange warning threshold range, match the pre-stored support reinforcement strategy library according to the parameter combination mode, dynamically calculate the required number of bolt reinforcement and grouting pressure values, and generate a support parameter adjustment instruction set; When all three parameters reach the red warning threshold range and the acceleration of the inner angle coefficient of variation exceeds the critical value, send an evacuation instruction to the on-site personnel terminal, automatically cut off the power supply of the tunnel face construction equipment, and activate the rapid grouting and hydraulic support units of the emergency support system; Dynamically update the topology structure database according to the warning level, record the warning trigger time, parameter overrun value and disposal measures to the blockchain evidence storage node, and form a warning decision chain.

[0089] In the embodiment of the present invention, extract the deformation parameter sequences (such as displacement difference change rate, curvature ratio gradient, etc.) of the past 3 years from the historical monitoring data, combine the support structure failure cases of similar projects (such as the parameter values when bolts break and shotcrete cracks), and determine the three-level threshold ranges of each parameter through statistical analysis (such as quantile method, ROC curve): Yellow warning (local deformation): The parameter value exceeds 25% of the historical extreme value and continues to rise; Orange warning (progressive instability): The parameter value exceeds 50% of the historical extreme value and shows an accelerating trend; Red warning (overall instability): The parameter value exceeds 90% of the historical extreme value and breaks through the engineering design safety limit.

[0090] Extract the three component parameters at the current moment in real time: the increment of the diagonal displacement difference change rate (such as 0.4mm / h), the change gradient of the curvature ratio (such as -0.1 / hour), and the acceleration of the inner angle coefficient of variation (such as 0.06 / hour 2 ) The threshold setting combines historical laws and engineering experience to avoid the limitations of a single fixed threshold (such as using different thresholds for slopes of different surrounding rock grades).

[0091] Yellow warning: Primary response and precise monitoring, implementation process: When any parameter first enters the yellow warning range (such as the increment of the displacement difference change rate > 0.3mm / h), the system automatically increases the three-dimensional laser ranging frequency of the corresponding monitored side from 1Hz to 5Hz to encrypt data collection to capture subtle deformation trends.

[0092] Associate the warning parameter with the coordinates of the high-risk node in the reference topology (such as x = 1234.5m, y = 67.8m, z = 9.0m at point P2), generate a primary warning signal containing the precise location (such as "Node P2 on the left side wall of a certain section enters the yellow warning"), and push it to the on-site engineer via text message / APP.

[0093] Confirm the deformation trend through high-frequency monitoring to avoid misjudgment (for example, a sudden increase in a single parameter caused by short-term vibration may be excluded by subsequent data).

[0094] Orange warning: Intelligent matching and support adjustment are achieved as follows: When two or more parameters simultaneously enter the orange warning range (such as the displacement difference increment > 0.5 mm / h and the curvature ratio gradient < -0.15 / hour), the system, based on the parameter combination mode (such as the "displacement acceleration + curvature anomaly" corresponding to the rock mass shear slip mode), matches the optimal solution from the pre-stored support reinforcement strategy library: If the "concave sidewall" mode is matched, automatically calculate the required number of bolt reinforcements (such as adding 2 bolts per meter) and the grouting pressure value (such as increasing from 0.3 MPa to 0.6 MPa); Generate a support parameter adjustment instruction set and send it to the on-site intelligent terminal through the industrial Ethernet to drive the robotic arm or grouting equipment to automatically adjust.

[0095] The strategy library is based on the "deformation mode - support measure" mapping relationship (such as 10 typical instability modes corresponding to 20 sets of adjustment solutions), realizing a closed-loop from warning to disposal and avoiding the lag of manual decision-making.

[0096] Red warning: Emergency disposal and system linkage are achieved as follows: When all three parameters enter the red warning range (such as the displacement difference increment > 1.0 mm / h, the curvature ratio gradient < -0.2 / hour, the coefficient of variation acceleration > 0.1 / hour 2 ) and the inner angle coefficient of variation acceleration exceeds the critical value (such as 0.08 / hour 2 , indicating that the angle anomaly deteriorates rapidly), the system triggers a triple emergency response: Send audible and visual alarms and text evacuation instructions (such as "Evacuate immediately to the refuge chamber") to all on-site personnel terminals; automatically cut off the power supply of the face construction equipment (such as excavators, loaders) through the PLC control system to prevent equipment out-of-control during instability; start the rapid grouting system (reach the design pressure of 1.0 MPa within 5 minutes) and the hydraulic support unit (automatically erect steel arch frames) to form a temporary rigid support.

[0097] The simultaneous overlimit of multi-dimensional parameters is used as the "overall instability" criterion, combined with acceleration parameters to identify the precursor of mutation, ensuring that the response is triggered in the last window period (such as starting emergency measures 10 minutes in advance in a certain tunnel case to avoid a collapse accident).

[0098] Data storage and dynamic update are achieved as follows: Dynamically update the topological structure database according to the warning level, record key data such as current deformation parameters, vertex weights, and radius of curvature (e.g., mark the cross-section as "high-risk state" during red warning); encrypt information such as warning trigger time, parameter overlimit value, and disposal measures and write them into the blockchain evidence storage node to form an unforgeable warning decision chain. Blockchain technology ensures data integrity (e.g., prevent manual tampering of warning records) and provides a credible data source for engineering safety audits.

[0099] Multi-dimensional criteria are used to avoid misjudgment of a single parameter (e.g., only exceeding the displacement standard may be caused by insufficient temperature compensation, and combining the curvature ratio and coefficient of variation can exclude environmental interference). The threshold is dynamically adjusted based on real-time data and historical cases (e.g., automatically lower the orange warning threshold during the rainy season) to adapt to working condition changes under complex geological conditions.

[0100] Only increase the monitoring frequency for the edges / nodes where anomalies first appear. Compared with full-range high-frequency monitoring, it reduces the energy consumption of sensors and data transmission pressure by 60%. The warning signal with coordinates enables on-site personnel to quickly locate the problem area (e.g., arrive at a certain side wall within 3 minutes), shortening the troubleshooting time by 50%.

[0101] Automatic linkage from warning to equipment adjustment (response time < 30 seconds) to timely contain the trend during the deformation acceleration stage (e.g., in a certain slope case, the displacement difference increment is reduced from 0.6 mm / h to 0.2 mm / h through bolt reinforcement).

[0102] The above are the preferred embodiments of the present invention. It should be noted that for those of ordinary skill in the art, without departing from the principle of the present invention, several improvements and refinements can be made, and these improvements and refinements should also be regarded as the protection scope of the present invention.

Claims

1. A road slope data monitoring method under complex geological conditions, characterized in that: The method comprises: Real-time measurement of the length of each side and the internal angle of the quadrilateral to generate the benchmark topology; Based on the benchmark topology and the real-time distance between each vertex of the quadrilateral and the tunnel face, the weight of each vertex is calculated; when the weight of the vertex is greater than the threshold, the change rate of the corresponding edge is monitored in real time; the change rate is temperature compensated, and the deformation caused by temperature change is calculated by analyzing the thermal expansion characteristics of the material, and the monitored edge length is corrected to generate the temperature compensated edge length; according to the temperature compensated edge length and the surrounding rock rheological equation, the adjacent two edge lengths are corrected by the strain difference of the adjacent edges to obtain the final edge length; the final edge length is combined with the benchmark internal angle to generate the topological structure; Based on the topological structure, a three-dimensional deformation index matrix including diagonal displacement difference, adjacent edge curvature ratio and internal angle variation coefficient is constructed; the index matrix, support structure strain spatiotemporal matrix and rock mass fragmentation tensor are input into the spatiotemporal attention graph neural network, and the deformation evolution trend including diagonal displacement difference, adjacent edge curvature ratio and internal angle variation coefficient is output; Multi-level early warning decisions are made based on component parameters of deformation evolution trends.

2. The method for monitoring road slope data under complex geological conditions according to claim 1, characterized in that: The reference topology is: The arch crown settlement point P1, the left convergence point P2, the right convergence point P3 and the invert arch uplift point P4 are arranged at the selected positions of the tunnel section, and the initial length of each side of the quadrilateral and the internal angle between adjacent sides are measured in real time by a three-dimensional laser rangefinder to establish an initial quadrilateral geometric parameter database; The dual-axis inclinometer array between P2 and P3 is used to collect the inclination angle changes of the side walls on both sides, and the 3D section scanner is synchronously triggered to perform point cloud scanning on the section where the quadrilateral is located. The scanning data is spatially registered with the laser ranging results to obtain calibration data. According to the calibration data, the lengths of the sides and the internal angles of the initial quadrilateral are dynamically corrected to generate a reference topology including the spatial coordinates of the edges, the values ​​of the internal angles and the curvature of the cross section.

3. The method for monitoring road slope data under complex geological conditions according to claim 2, characterized in that: Based on the reference topology and the real-time distance between each vertex of the quadrilateral and the tunnel face, the weight of each vertex is calculated, including: The three-dimensional Euclidean distance between the quadrilateral vertices P1, P2, P3, and P4 and the tunnel face is continuously measured by a three-dimensional laser rangefinder at a sampling frequency of 5 times per second to generate the original dynamic distance sequence of each vertex. The original sequence is subjected to sliding window average filtering with a window size of 10 sampling points to generate smoothed dynamic distance time series data. Based on the propagation attenuation characteristics of rock stress waves, a negative exponential mapping relationship between vertex weight and real-time distance is established, and the stress attenuation curve is fitted to determine the weight function parameters. The smoothed dynamic distance time series data is input into the weight function parameters to calculate the weight value of each vertex in real time.

4. The method for monitoring road slope data under complex geological conditions according to claim 3 is characterized in that: When the weight of a vertex is greater than a threshold, the change rate of the corresponding edge is monitored in real time; Temperature compensation is performed on the rate of change. By analyzing the thermal expansion characteristics of the material, the deformation caused by temperature change is calculated, including: Normalize the weight value of each vertex to the range of 0-1. When the vertex weight exceeds the preset deformation sensitivity threshold, mark the vertex as a high-risk node and activate the real-time change rate monitoring module of the edges connected to it. At the same time, bind the weight data to the spatial coordinates in the baseline topology to generate a dynamic topology node information table with weights. According to the marked high-risk vertices, locate the quadrilateral edges connected to them; start the high-frequency monitoring mode of the three-dimensional laser rangefinder for the target edge, continuously collect the original data of the edge length at a sampling frequency of 1Hz, and use the difference calculation between two adjacent sampling points to generate the original length change rate sequence; remove outliers from the original change rate data, and the removal rule is: if the change rate of a single point exceeds 3 times the standard deviation of the historical mean value for the same period, it is marked as noise and repaired by linear interpolation to obtain the processed change rate data; Distributed temperature sensor arrays are arranged on the surface and inside of the rock mass on both sides of the target edge to collect temperature field data of the entire monitored edge and generate a temperature spatiotemporal distribution matrix; the pre-stored rock material thermal expansion coefficient library is called to match the corresponding thermal expansion coefficient according to the material type of the current monitored edge; the predicted deformation caused by temperature fluctuations is calculated point by point based on the thermal expansion coefficient, the original reference side length and the real-time temperature difference, where the real-time temperature difference is the difference between the current temperature and the initial temperature when the reference topology is generated.

5. The method for monitoring road slope data under complex geological conditions according to claim 4, characterized in that: Correct the monitored side length to generate temperature compensated side length, including: Extract the time series temperature data collected by the temperature sensor and the corresponding theoretical temperature deformation; Based on the thermal expansion characteristics of rock materials, the theoretical temperature deformation is integrated in a sliding time window according to the time series, and the instantaneous deformation in each time slice is accumulated to generate a temperature compensation cumulative value curve; Obtain the original side length data measured in real time by the 3D laser rangefinder, and align it with the temperature compensation accumulated value according to the timestamp; Perform compensation superposition calculation: temperature compensation side length = original side length + accumulated compensation value, where the sign of the compensation value is determined according to the direction of temperature rise and fall. When the temperature rises and causes expansion, the compensation value is negative, otherwise it is positive.

6. The method for monitoring road slope data under complex geological conditions according to claim 5, characterized in that: According to the temperature-compensated side length and the surrounding rock rheological equation, the lengths of two adjacent sides are corrected by the strain difference between the adjacent sides to obtain the final side length; The final edge length is combined with the base internal angle to generate a topological structure, including: The distributed optical fiber strain sensor array is used to collect the axial strain distribution data of two adjacent edges in real time, and the average strain difference between the two edges is calculated based on the collected strain distribution data; Based on the average strain difference, the time-dependent deformation parameters in the surrounding rock rheological equation, including the viscoelastic coefficient and the creep rate, are called, and the proportion of the coordinated deformation of adjacent edges caused by the rheological effect is calculated by combining the time difference between the current monitoring time and the benchmark topology generation time. According to the proportion of cooperative deformation of adjacent edges, the change in the length of adjacent edges caused by rheological effect is derived; Taking the temperature-compensated side length as input and combining it with the change in the length of the adjacent sides, the rheological correction model is used to reversely correct the lengths of the two adjacent sides to obtain the corrected side length. Based on the corrected side length and reference internal angle data, with P1 as the coordinate origin, the spatial coordinates of P2, P3, and P4 are calculated vertex by vertex; the curvature of the quadrilateral section is fitted using the least squares method to update the curvature radius parameter; The corrected side lengths, updated internal angles, and curvature radius parameters are integrated to generate a dynamic topological structure containing geometric parameters and curvature information.

7. The method for monitoring road slope data under complex geological conditions according to claim 6, characterized in that: Construct a three-dimensional deformation index matrix containing diagonal displacement difference, adjacent edge curvature ratio and internal angle variation coefficient, including: The current lengths of the two diagonals of the quadrilateral are extracted from the dynamic topology structure and compared with the initial diagonal lengths in the reference topology; the absolute displacement difference and displacement difference change rate of each diagonal are calculated to generate the displacement difference feature vectors of the two diagonals; Based on the updated cross-section curvature parameters in the dynamic topology, the curvature radius and curvature ratio of the two adjacent sides are obtained, the curvature ratio is normalized, the ratio is mapped to the interval of 0-1, and the curvature ratio feature vector is generated; Read the current angle values ​​of the four inner angles of the quadrilateral from the dynamic topology, calculate the inner angle standard deviation and the inner angle coefficient of variation, and perform sliding window statistics on the coefficient of variation to generate a smoothed coefficient of variation sequence; Align the diagonal displacement difference eigenvector, curvature ratio eigenvector and internal angle variation coefficient sequence according to timestamps; construct a three-dimensional deformation index matrix; The three-dimensional deformation index matrix is ​​subjected to data normalization to generate a normalized deformation index matrix.

8. The method for monitoring road slope data under complex geological conditions according to claim 7, characterized in that: The dimension of the three-dimensional deformation index matrix is ​​[timestamp × deformation feature], where each row contains: the first column: the displacement difference change rate of the two diagonals; the second column: the normalized value of the curvature ratio of adjacent edges; the third column: the coefficient of variation of the internal angle.

9. The method for monitoring road slope data under complex geological conditions according to claim 8, characterized in that: The index matrix, the spatiotemporal matrix of the support structure strain, and the rock mass fragmentation tensor are input into the spatiotemporal attention graph neural network, and the output includes the deformation evolution trend of the diagonal displacement difference, the curvature ratio of adjacent edges, and the internal angle variation coefficient, including: The standardized deformation index matrix, the support structure strain time-space matrix constructed by the support structure strain time series data collected by the optical fiber sensor, and the rock mass fragmentation tensor calculated based on the microseismic event energy distribution and rock mass crack density are aligned according to the timestamp to obtain time-aligned data; Based on the time-aligned data, the rock mass fragmentation tensor is subjected to dimensionality reduction processing, and its principal component eigenvector is extracted to make its column dimension match the deformation index matrix; Using the principal component eigenvectors, together with the deformation index matrix and the support structure strain space-time matrix, the three types of data features are spliced ​​along the time axis to construct a fused space-time tensor; Based on the fused spatiotemporal tensor, a multi-head attention mechanism is used to calculate the spatiotemporal dependencies between the three types of data. According to the spatiotemporal dependency, the fused spatiotemporal tensor is multiplied by the attention weight matrix, and the spatial correlation features of the quadrilateral vertices and their supporting nodes are aggregated through the graph convolutional network; Based on spatial correlation features, gated recurrent units are used to learn the long-term evolution trend of deformation indicators; The long-term evolution trend is input into the fully connected prediction module, and the deformation evolution trend parameters within the next 3 hours are output, including: the incremental rate of change of the diagonal displacement difference, the gradient of the adjacent edge curvature ratio, and the acceleration of the internal angle variation coefficient.

10. The method for monitoring road slope data under complex geological conditions according to claim 9, characterized in that: Execute multi-level early warning decisions based on component parameters of deformation evolution trends, including: Extract the incremental change rate of diagonal displacement difference, the change gradient of adjacent edge curvature ratio and the acceleration parameters of internal angle variation coefficient in the deformation evolution trend; based on the statistical correlation analysis of historical deformation data and support structure failure cases, set the graded threshold intervals of the three parameters, including the yellow warning threshold interval, the orange warning threshold interval and the red warning threshold interval; When any parameter enters the yellow warning threshold interval for the first time, the 3D laser ranging frequency of the corresponding monitoring edge will be increased; the warning information will be associated with the coordinates of the high-risk nodes in the benchmark topology to generate a primary warning signal with spatial positioning; When two or more parameters simultaneously enter the orange warning threshold range, the pre-stored support reinforcement strategy library is matched according to the parameter combination mode, the required anchor reinforcement quantity and grouting pressure value are dynamically calculated, and the support parameter adjustment instruction set is generated; When all three parameters reach the red warning threshold range and the acceleration of the internal angle variation coefficient exceeds the critical value, an evacuation command is sent to the on-site personnel terminal, the power supply of the face construction equipment is automatically cut off, and the rapid grouting and hydraulic support unit of the emergency support system is turned on; The topological structure database is dynamically updated according to the warning level, and the warning trigger time, parameter limit value and disposal measures are recorded in the blockchain evidence node to form a warning decision chain.

Images