Multi-sensing safety monitoring system for pumped storage power station

Abstract

The invention relates to the technical field of pumped storage power station multi-sensor monitoring network time-space synchronous measurement, in particular to a pumped storage power station multi-sensor safety monitoring system, which comprises a sensor data acquisition module, a time-space synchronous network module, a data calibration module, a dynamic association model construction module, a multi-dimensional data fusion module and a graded early warning module. The dynamic correlation model building module simulates dynamic response of seepage pressure and surrounding rock strain through a fluid-structure interaction equation and a discrete element digital twinborn model, updates fracture permeability boundary conditions in real time, and extracts cross-sensor space-time correlation characteristics by adopting a convolutional long-short-term memory network improved by a multi-head attention mechanism. According to the method, the real-time early warning precision of abnormal events and the reliability of structural health assessment in the complex hydraulic and mechanical coupling environment are improved.

Description

Technical Field

[0001] The present invention relates to the technical field of spatio-temporal synchronous measurement of multi-sensor monitoring networks in pumped-storage power stations, and particularly to a multi-sensor safety monitoring system for pumped-storage power stations. Background Art

[0002] In existing multi-sensor safety monitoring systems for pumped-storage power stations, dynamic acquisition and analysis of key structural parameters are achieved by integrating a variety of high-precision sensors. This system uses devices such as strain sensors, displacement gauges, and piezometers to monitor key parts such as the surrounding rock mass of the underground powerhouse, penstock, dam, and water conveyance tunnel in real time, and combines temperature sensors and vibration sensors to obtain environmental variable and mechanical operation state data. Through the fusion processing of multiple data, a correlation model between structural response and load action is established, and pattern recognition algorithms are used to identify abnormal signals, thereby triggering a hierarchical early warning mechanism. The core of its technology lies in eliminating the uncertainty of single-sensor data through spatio-temporal correlation analysis, improving the reliability of monitoring results, providing a quantitative basis for structural health assessment and risk prevention and control during the whole life cycle of the power station, and effectively reducing potential safety hazards caused by material aging, hydraulic impact, or geological activities.

[0003] In existing multi-sensor safety monitoring systems for pumped-storage power stations, the spatio-temporal synchronization and correlation of sensor data are insufficient, resulting in limited identification accuracy of the coupling mechanism of abnormal events. In a complex hydro-mechanical coupling environment, sensors of different physical quantities generate timing deviations and spatial calibration errors due to differences in acquisition frequencies or transmission delays, thereby weakening the synergistic effect of data fusion. For example, the asynchronous acquisition of the deformation of the surrounding rock mass of a pumped-storage power station and the seepage data of the adjacent water conveyance tunnel masks the dynamic interaction characteristics between the two, cannot accurately reflect the correlation between hydraulic fracturing and rock mass displacement, and reduces the timeliness of early warning of structural instability precursors. Summary of the Invention

[0004] Aiming at the deficiencies of the existing technology, the present invention provides a multi-sensor safety monitoring system for pumped-storage power stations, and the present invention solves the problem of deviation in the identification of the coupling mechanism of abnormal events caused by insufficient spatio-temporal synchronization accuracy of the multi-sensor monitoring network in a pumped-storage power station under a complex hydro-mechanical coupling environment.

[0005] To solve the above technical problems, the specific technical solutions of the present invention are as follows: The multi-sensor safety monitoring system for pumped-storage power stations provided by the present invention includes: a sensor data acquisition module, a spatio-temporal synchronization network module, a data calibration module, a dynamic correlation model construction module, a multi-dimensional data fusion module, and a hierarchical early warning module; The sensor data acquisition module is used to collect stress, seepage pressure, displacement and vibration signals of key monitoring points through redundant sensor nodes, and perform noise suppression processing on the collected signals to generate preprocessed data; the spatio-temporal synchronization network module is used to receive the preprocessed data output by the sensor data acquisition module, achieve microsecond-level time alignment of sensor nodes through a hierarchical clock synchronization architecture, and dynamically adjust the data transmission timing based on cluster-based communication optimization, and output synchronized timing data; The data calibration module is used to receive the synchronized timing data output by the spatio-temporal synchronization network module, align asynchronous data through a transmission delay compensation algorithm, and combine a three-dimensional laser scanner to calibrate the global coordinate system of the sensor nodes to generate spatio-temporal aligned data; The dynamic correlation model construction module is used to input the spatio-temporal aligned data output by the data calibration module into a fluid-structure interaction model, drive a digital twin model based on the discrete element method to simulate the dynamic response of seepage pressure and surrounding rock strain, and output simulation verification data; The multi-dimensional data fusion module is used to extract the spatio-temporal correlation features of the spatio-temporal aligned data and the simulation verification data, and generate fusion decision data through a weighted fusion algorithm; The hierarchical early warning module is used to trigger multi-level early warning thresholds according to the fusion decision data, and iteratively optimize the clock synchronization parameters of the spatio-temporal synchronization network module and the fluid-structure interaction model coefficients of the dynamic correlation model construction module based on the fusion decision data.

[0006] Furthermore, for the multi-sensor safety monitoring system of the pumped storage power station of the present invention, the sensor data acquisition module includes: A redundant sensor node group, an abnormal data elimination unit and an adaptive noise suppression unit, where: The redundant sensor node group adopts a triple redundancy configuration and is deployed in the stress concentration area of the penstock. Three sensors of the same type are set at each monitoring point to form a data acquisition unit, and the three sensors synchronously collect signals of the same physical quantity; The abnormal data elimination unit is connected to the redundant sensor node group and eliminates abnormal data through a majority voting algorithm within a sliding time window. Specifically: normalize the data collected by the three sensors within the sliding time window, calculate the Euclidean distance between two-by-two data, and if the distance of a certain data exceeds the dynamic threshold from the other two data, it is determined as abnormal data and eliminated, and the mean value of the remaining data is used as the effective output; The adaptive noise suppression unit is connected to the abnormal data rejection unit, and performs noise reduction processing on the effective output signal by using an adaptive wavelet threshold filtering algorithm, including: dynamically selecting a wavelet basis function based on the signal spectrum analysis result, separating the electromagnetic interference noise frequency band through an improved SUREShrink threshold estimation, and retaining the time-domain waveform characteristics and frequency-domain harmonic characteristics when reconstructing the signal, generating preprocessed data and outputting it to the spatio-temporal synchronization network module.

[0007] Further, for the multi-sensor safety monitoring system of the pumped storage power station of the present invention, the spatio-temporal synchronization network module includes a clock synchronization sub-module and a cluster communication optimization sub-module, where: The clock synchronization sub-module is used to receive the preprocessed data output by the sensor data acquisition module, including: The master node obtains a reference clock signal from the Beidou satellite system through the satellite time synchronization protocol, and broadcasts the reference clock to all slave nodes; Each slave node establishes a time synchronization link with the master node through the precision clock synchronization protocol, and at the same time uses a linear regression model to predict the cumulative drift of the crystal oscillator of each slave node in real time, and dynamically corrects the timestamp of the sensor data according to the drift for microsecond-level time alignment between sensor nodes; The cluster communication optimization sub-module is connected to the clock synchronization sub-module, including: Dividing the monitoring area into multiple geological unit subnets according to the geological structure characteristics, and forming independent communication clusters for the sensor nodes within each subnet; Within the subnet, fixed data transmission time slots are allocated to each node through the time division multiple access mechanism. The master node monitors the channel load intensity in real time. When the data volume of the nodes within the subnet exceeds the preset threshold, a dynamic time slice adjustment algorithm is used to shorten the length of a single time slot and increase the time slot density to optimize the data transmission timing; Output the data that has been time-aligned and timing-optimized as synchronous timing data, and transmit it to the data calibration module.

[0008] Further, for the multi-sensor safety monitoring system of the pumped storage power station of the present invention, the data calibration module includes a time delay compensation sub-module and a space calibration sub-module, where: The time delay compensation sub-module is used to receive the synchronous timing data output by the spatio-temporal synchronization network module, including: Measuring the round-trip transmission link delay between the master node and the slave node by using the two-way timestamp protocol, fitting the delay fluctuation curve through the least squares method to generate a reference transmission delay amount; Predicting the random delay deviation caused by network jitter based on the Kalman filtering algorithm, constructing an interpolation function in combination with the reference transmission delay amount, and performing time-axis resampling on the asynchronously arriving sensor data for microsecond-level time alignment; The spatial calibration sub-module is connected to the time delay compensation sub-module and includes: Obtain the three-dimensional point cloud data of the monitoring area through a three-dimensional laser scanner, and extract the actual point cloud coordinates of the installation positions of the sensor nodes; Adopt the iterative closest point algorithm with normal vector constraint to register the actual point cloud coordinates with the theoretical coordinates in the design drawing, and calculate the coordinate offset of the sensor nodes; Based on the quaternion rotation matrix, convert the local coordinate systems of the sensor nodes to the global engineering coordinate system, and output the spatio-temporal alignment data to the dynamic correlation model construction module.

[0009] Furthermore, for the multi-sensor safety monitoring system of the pumped storage power station of the present invention, the dynamic correlation model construction module includes a fluid-structure interaction solution sub-module and a digital twin driving sub-module, where: The fluid-structure interaction solution sub-module is used to receive the spatio-temporal alignment data output by the data calibration module, and includes: based on the coupled partial differential equation of the seepage pressure gradient tensor and the surrounding rock strain tensor, calculate the dynamic correlation relationship between the seepage field and the stress field through an implicit iterative solver, and generate the coupled response data of seepage pressure and strain; The digital twin driving sub-module is connected to the fluid-structure interaction solution sub-module and includes: constructing a three-dimensional digital twin model of the surrounding rock fracture propagation by using the discrete element method, mapping the real-time seepage pressure data as the permeability boundary condition of the fracture surface, and simulating the non-linear relationship between permeability and strain during the fracture propagation process through the dynamic relaxation algorithm; Update the fracture network topology of the digital twin model according to the coupled response data of seepage pressure and strain, drive the dynamic evolution of the digital twin model, and output the simulation verification data to the multi-dimensional data fusion module.

[0010] Furthermore, for the multi-sensor safety monitoring system of the pumped storage power station of the present invention, the multi-dimensional data fusion module includes a spatio-temporal feature extraction sub-module and a dynamic weighted fusion sub-module, where: The spatio-temporal feature extraction sub-module is used to receive the simulation verification data output by the dynamic correlation model construction module and the spatio-temporal alignment data output by the data calibration module, and includes: Adopt a convolutional long short-term memory network improved by the multi-head attention mechanism, input the seepage pressure time series data and the displacement rate spatial distribution data, extract the local deformation features through the spatio-temporal convolutional kernel, and use the multi-head attention mechanism to capture the correlation between pressure and displacement across sensor nodes, and generate a spatio-temporal correlation feature matrix; The dynamic weighted fusion sub-module is connected to the spatio-temporal feature extraction sub-module and includes: calculating the model confidence weight according to the error rate between the simulation verification data and the measured data, and calculating the data reliability weight based on the historical stability index of the sensor nodes; Perform a normalized product operation on the model confidence weight and the data reliability weight to generate the dynamic weighted fusion coefficient of each sensor node; Perform multi-dimensional fusion on the spatio-temporal correlation feature matrix through a weighted decision function, and output the fusion decision data to the hierarchical early warning module.

[0011] Furthermore, for the multi-sensor safety monitoring system of the pumped storage power station of the present invention, the hierarchical early warning module includes an early warning trigger sub-module and a closed-loop parameter optimization sub-module, where: The early warning trigger sub-module is used to receive the fusion decision data output by the multi-dimensional data fusion module, including: Divide the three-level early warning thresholds for seepage erosion, crack expansion, and rock mass instability based on the structural safety margin model, where: The seepage erosion threshold is determined by calculating the ratio of the seepage velocity to the critical permeability through Darcy's law; The crack expansion threshold is based on the fracture mechanics J-integral model to calculate the stress intensity factor at the crack tip; The rock mass instability threshold determines the safety factor of the surrounding rock shear slip surface according to the Mohr-Coulomb criterion; Compare the fusion decision data with the three-level thresholds in real time. When the continuous over-limit duration of the data exceeds the preset tolerance, trigger the corresponding level of audible and visual alarm signals; The closed-loop parameter optimization sub-module is connected to the early warning trigger sub-module, including: using the clock synchronization error rate of the spatio-temporal synchronization network module and the root mean square of the prediction residuals of the dynamic correlation model construction module as the fitness function, and using the genetic algorithm for multi-objective optimization; Select elite individuals from the population through the tournament selection strategy. The crossover operator uses arithmetic recombination, and the mutation operator introduces Gaussian perturbation to iteratively search for the optimal combination of clock synchronization parameters and fluid-structure coupling model coefficients; Real-time feedback the optimized parameters to the clock synchronization sub-module of the spatio-temporal synchronization network module and the digital twin drive sub-module of the dynamic correlation model construction module.

[0012] Advantages of the present invention; The present invention effectively eliminates abnormal data caused by single-sensor failure or environmental interference through the triple redundancy configuration of redundant sensor node groups combined with the sliding time window majority voting algorithm, improving the reliability of multi-source heterogeneous sensor data acquisition; the spatio-temporal synchronization network module adopts a hierarchical clock synchronization architecture and a cluster-based communication optimization mechanism, predicts the crystal oscillator drift through the satellite time synchronization protocol and the linear regression model, realizes microsecond-level time alignment and dynamically adjusts the data transmission timing, and solves the multi-node communication congestion problem in complex geological environments; the data calibration module compensates for the transmission delay based on the bidirectional timestamp protocol and the Kalman filtering algorithm, and completes the global coordinate system calibration by combining the normal vector constrained iterative closest point algorithm of the three-dimensional laser scanner, eliminating the influence of spatio-temporal reference deviation on data fusion; the dynamic correlation model construction module collaboratively simulates through the fluid-structure interaction equation and the discrete element digital twin model, updates the fracture permeability boundary conditions in real time and simulates the non-linear correlation between seepage pressure and surrounding rock strain, enhancing the analytical ability of the hydraulic and mechanical coupling mechanism; the multi-dimensional data fusion module uses a spatio-temporal convolutional network improved by the multi-head attention mechanism to extract cross-sensor correlation features, combines the dynamic weighted fusion algorithm to suppress the interference of low-confidence data, and optimizes the identification accuracy of abnormal event coupling features; the hierarchical early warning module establishes a three-level threshold judgment rule based on the structural safety margin model, and realizes the closed-loop optimization iteration of the clock synchronization parameters and the model coefficients through the genetic algorithm, improving the adaptability and early warning timeliness of the system under complex working conditions. BRIEF DESCRIPTION OF THE DRAWINGS

[0013] In order to more clearly illustrate the technical solutions of the present invention, the drawings required for use in the embodiments will be briefly introduced below. Obviously, for those of ordinary skill in the art, other drawings can also be obtained based on the drawings without creative efforts.

[0014] Figure 1 FIG. is a system architecture diagram of the multi-sensor safety monitoring system for a pumped storage power station provided by an embodiment of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0015] In order to make the objectives, technical solutions, and advantages of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below in conjunction with the specific embodiments and corresponding drawings of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts belong to the scope of protection of the present invention. The technical solutions provided by each embodiment of the present invention will be described in detail below in conjunction with the drawings. In order to better understand the objectives of the present invention, the present invention will be further described in detail below.

[0016] Please refer to Figure 1, the multi-sensor safety monitoring system provided by the present invention for pumped storage power stations includes: a sensor data acquisition module, a spatio-temporal synchronization network module, a data calibration module, a dynamic correlation model construction module, a multi-dimensional data fusion module, and a hierarchical early warning module; The sensor data acquisition module is used to collect stress, seepage pressure, displacement and vibration signals of key monitoring points through redundant sensor nodes, and perform noise suppression processing on the collected signals to generate preprocessed data; The spatio-temporal synchronization network module is used to receive the preprocessed data, achieve microsecond-level time alignment of sensor nodes through a hierarchical clock synchronization architecture, and dynamically adjust the data transmission timing through cluster-based communication optimization to output synchronized timing data; The data calibration module is used to receive the synchronized timing data, align asynchronous data through a transmission delay compensation algorithm, and calibrate the global coordinate system of sensor nodes through a three-dimensional laser scanner to generate spatio-temporal aligned data; The dynamic correlation model construction module is used to input the spatio-temporal aligned data into a fluid-structure interaction model, drive a digital twin model based on the discrete element method to simulate the dynamic responses of seepage pressure and surrounding rock strain, and output simulation verification data; The multi-dimensional data fusion module is used to extract the spatio-temporal correlation features of the spatio-temporal aligned data and the simulation verification data, and generate fusion decision data through a weighted fusion algorithm; The hierarchical early warning module is used to trigger multi-level early warning thresholds according to the fusion decision data, and adjust the clock synchronization parameters of the spatio-temporal synchronization network module and the model coefficients of the dynamic correlation model construction module through a parameter optimization algorithm.

[0017] The sensor data acquisition module synchronously collects stress, seepage pressure, displacement and vibration signals of key monitoring points of the pumped storage power station through redundant sensor nodes. The redundant sensor nodes adopt a triple-redundancy configuration and are deployed in the stress concentration area of the penstock and the surrounding rock deformation sensitive area. Three sensors of the same type are set at each monitoring point to form a data acquisition unit, and abnormal data is eliminated through the majority voting algorithm within a sliding time window. After eliminating the abnormal data, an adaptive wavelet threshold filtering algorithm is used to denoise the effective signals, the wavelet basis function is dynamically selected based on the signal spectrum analysis results, the electromagnetic interference noise frequency band is separated, and the time-domain waveform features and frequency-domain harmonic features are retained to generate preprocessed data.

[0018] After receiving the preprocessed data, the spatio-temporal synchronization network module realizes the microsecond-level time alignment of sensor nodes through a hierarchical clock synchronization architecture. The master node obtains the reference clock signal through the satellite time service protocol, and the slave nodes are synchronized to the master node through the precision clock synchronization protocol. A linear regression model is used to predict the crystal oscillator drift and correct the timestamps. The clustered communication optimization divides the monitoring subnet according to geological units. The time slots for data transmission are allocated within the subnet through the time division multiple access mechanism. The master node dynamically adjusts the length of the time slice according to the channel load, optimizes the data transmission timing, and outputs the synchronized timing data.

[0019] After receiving the synchronized timing data, the data calibration module aligns the asynchronous data through the transmission delay compensation algorithm. The two-way timestamp protocol is used to measure the link delay. The Kalman filter is used to predict the network jitter and generate an interpolation function for time axis resampling. On this basis, the point cloud data of the monitoring area is obtained by a 3D laser scanner. The iterative closest point algorithm with normal vector constraint is used to register the actual coordinates and the designed coordinates of the sensor nodes, and a quaternion rotation matrix is constructed to transform the local coordinate system to the global engineering coordinate system, generating the spatio-temporal alignment data.

[0020] The dynamic correlation model construction module inputs the spatio-temporal alignment data into the fluid-structure coupling model, and calculates the dynamic correlation equation of the seepage pressure gradient and the surrounding rock strain tensor through an implicit iterative solver. The digital twin model simulates the process of surrounding rock fracture propagation based on the discrete element method, maps the real-time seepage pressure data as the permeability boundary condition of the fracture surface, drives the dynamic evolution of the model, generates the coupled response data of seepage and strain, and outputs the simulation verification data.

[0021] The multi-dimensional data fusion module extracts the spatio-temporal correlation features of the spatio-temporal alignment data and the simulation verification data. A convolutional long short-term memory network improved by the multi-head attention mechanism is used to capture the cross-sensor correlation between the sudden change of seepage pressure and the displacement rate. The weighted fusion coefficient is generated according to the simulation error rate and the historical stability index of the sensor. Through the normalized product operation, the model confidence and the data reliability weight are combined to perform multi-dimensional fusion on the spatio-temporal feature matrix, generating the fusion decision data.

[0022] The hierarchical early warning module triggers multi-level early warning thresholds according to the fusion decision data. The seepage erosion threshold is determined by calculating the ratio of the seepage velocity to the critical permeability according to Darcy's law. The crack propagation threshold is based on the fracture mechanics J integral model to determine the stress intensity factor at the crack tip. The rock mass instability threshold evaluates the safety factor of the surrounding rock shear slip surface according to the Mohr-Coulomb criterion. The parameter optimization algorithm uses the spatio-temporal synchronization error and the model prediction residual as the fitness function, and uses the genetic algorithm combined with the tournament selection strategy to search for the optimal parameter combination, and feeds the optimized clock synchronization parameters and model coefficients back to the spatio-temporal synchronization network module and the dynamic correlation model construction module to form a closed-loop optimization link.

[0023] Specifically, for the multi-sensor safety monitoring system of the pumped-storage power station described in the present invention, the sensor data acquisition module includes: A redundant sensor node group, an abnormal data rejection unit, and an adaptive noise suppression unit, where: the redundant sensor node group adopts a triple-redundancy configuration and is deployed in the stress concentration area of the penstock. Three sensors of the same type are set at each monitoring point to form a data acquisition unit, and the three sensors synchronously collect signals of the same physical quantity; The abnormal data rejection unit is connected to the redundant sensor node group and rejects abnormal data through a majority voting algorithm within a sliding time window. Specifically: within the sliding time window, the data collected by the three sensors are normalized, and the Euclidean distance between two-by-two data is calculated. If the distance of a certain data from the other two data exceeds the dynamic threshold, it is determined as abnormal data and rejected, and the mean value of the remaining data is used as the effective output; The adaptive noise suppression unit is connected to the abnormal data rejection unit and uses an adaptive wavelet threshold filtering algorithm to denoise the effective output signal, including: dynamically selecting a wavelet basis function based on the signal spectrum analysis result, separating the electromagnetic interference noise frequency band through an improved SUREShrink threshold estimation, and retaining the time-domain waveform characteristics and frequency-domain harmonic characteristics during signal reconstruction to generate preprocessed data and output it to the spatio-temporal synchronization network module.

[0024] The redundant sensor node group of the sensor data acquisition module adopts a triple-redundancy configuration and is deployed in the stress concentration area of the penstock. Three sensors of the same type, including strain sensors, osmotic pressure sensors, and vibration sensors, are set at each monitoring point to form a data acquisition unit, and the synchronous acquisition of physical quantity signals is realized through a synchronous trigger circuit. The deployment of the redundant sensor node group in the stress concentration area covers the weld area and high-curvature area of the penstock, and dense layout is carried out for the structure deformation sensitive parts. The triple-redundancy design improves data reliability through cross-validation.

[0025] The abnormal data rejection unit processes the original acquisition data of the redundant sensor node group through a majority voting algorithm within a sliding time window. The length of the sliding time window is dynamically adjusted according to the sensor sampling frequency and the change rate of the physical quantity. The data of the three-way sensors within the window are normalized to eliminate the dimension difference. After normalization, the Euclidean distance between two-by-two data is calculated, and the dynamic threshold is generated based on the statistical characteristics of historical data. If the Euclidean distance of a certain data continuously exceeds the threshold, it is determined as abnormal data and rejected, and the mean value of the remaining valid data is taken as the output. This mechanism suppresses abnormal data caused by single sensor failure or environmental interference through multi-dimensional data comparison.

[0026] The adaptive noise suppression unit performs noise reduction processing on the effective signal after abnormal data rejection. Based on the signal spectrum analysis results, a wavelet basis function matching the main frequency of the signal is adaptively selected, such as the Daubechies wavelet or the Symlet wavelet. The electromagnetic interference noise frequency band is separated by an improved SUREShrink threshold estimation method. A frequency band energy weight factor is introduced in the threshold estimation process to preferentially suppress high-frequency noise components. At the same time, the mutation characteristics of the time-domain waveform and the frequency-domain harmonic components are retained during the signal reconstruction stage. The noise-reduced signal passes through joint time-frequency domain verification to generate preprocessed data that meets the requirements of subsequent processing and is output to the spatio-temporal synchronization network module.

[0027] The data acquisition of the redundant sensor node group, the majority voting mechanism of the abnormal data rejection unit, and the frequency-domain noise reduction processing of the adaptive noise suppression unit form a series link. The three-redundancy sensor synchronous acquisition provides a basis for data consistency. The sliding time window algorithm filters out outliers through a dynamic threshold, and the adaptive wavelet filtering further eliminates noise interference, finally generating preprocessed data with high confidence. The data transmission between units is realized through a standardized interface protocol to ensure the consistency of data format and processing timing.

[0028] Specifically, for the multi-sensor safety monitoring system of the pumped storage power station of the present invention, there are a clock synchronization sub-module and a cluster communication optimization sub-module, where: The clock synchronization sub-module is used to receive the preprocessed data of the sensor data acquisition module, including: The master node obtains the reference clock signal from the Beidou satellite system through the satellite time synchronization protocol and broadcasts the reference clock to all slave nodes; Each slave node establishes a time synchronization link with the master node through the precision clock synchronization protocol. At the same time, a linear regression model is used to predict the cumulative drift of the crystal oscillator of each slave node in real time, and the timestamp of the sensor data is dynamically corrected according to the drift for microsecond-level time alignment between sensor nodes; The cluster communication optimization sub-module is connected to the clock synchronization sub-module, including: The monitoring area is divided into multiple geological unit subnets according to the geological structure characteristics, and the sensor nodes within each subnet form an independent communication cluster; In the subnet, a time division multiple access mechanism is used to allocate fixed data transmission time slots for each node. The master node monitors the channel load intensity in real time. When the data volume of the nodes in the subnet exceeds the preset threshold, a dynamic time slice adjustment algorithm is used to shorten the length of a single time slot and increase the time slot density to optimize the data transmission timing; The data that has undergone time alignment and timing optimization is output as synchronous timing data and transmitted to the data calibration module.

[0029] The clock synchronization sub-module of the space-time synchronization network module establishes a global time reference through the satellite timing protocol. The master node accesses the navigation message signal of the Beidou satellite system, parses the Coordinated Universal Time (UTC) as the reference clock source, and distributes the reference time information to all slave nodes through a low-latency broadcast protocol. The slave nodes establish a two-way communication link with the master node using the Precision Clock Synchronization Protocol, periodically exchange time synchronization messages, measure the transmission path delay, and compensate for the network asymmetry error. The slave nodes are built with high-precision temperature-compensated crystal oscillators, analyze the historical clock deviation data through a linear regression model, predict the crystal oscillator frequency drift trend, dynamically correct the local clock counter value, and achieve microsecond-level alignment of the sensor data timestamps. During the time synchronization process, the master node continuously monitors the clock offset of the slave nodes, and triggers the resynchronization mechanism when the offset exceeds the tolerance range to maintain the time consistency of the entire network.

[0030] The cluster communication optimization sub-module divides the monitoring area subnet based on geological structure characteristics. According to the pressure steel pipe orientation, surrounding rock fault distribution, and underground powerhouse structure morphology, the monitoring area is divided into multiple geological unit subnets, and each subnet corresponds to an independent communication cluster. The sensor nodes within the subnet elect a cluster head node through the self-organizing network protocol, and the cluster head node is responsible for coordinating the time slot allocation under the time division multiple access mechanism. The fixed time slot length is preset according to the sensor sampling rate and data packet size. The master node real-time collects the channel load intensity indicators of each subnet. When the concurrent data volume within the subnet exceeds the preset threshold, the dynamic time slice adjustment algorithm is started. This algorithm predicts the data transmission queue length based on the queuing theory model, compresses the duration of a single time slot proportionally, and increases the number of time slots within a unit period to improve the channel utilization rate. The optimized timing parameters are sent to the subnet cluster head node through control messages to synchronously update the time slot mapping table of each node and achieve dynamic reconstruction of the data transmission timing.

[0031] The clock synchronization sub-module and the cluster communication optimization sub-module work together. The sensor data after timestamp correction is sent to the master node according to the optimized timing. After the master node integrates the synchronous timing data of each subnet, it adds a time reference label and a subnet identifier according to the standard data encapsulation format to generate a synchronous timing data stream with space-time reference information, which is transmitted to the data calibration module for subsequent processing. The time alignment mechanism eliminates the clock deviation of multi-source sensor data, and the cluster communication optimization reduces the probability of channel conflict. Both jointly ensure the real-time and integrity of high-density monitoring data.

[0032] Specifically, for the multi-sensor safety monitoring system of the pumped storage power station described in the present invention, the data calibration module includes: A time delay compensation sub-module and a space calibration sub-module, where: The time delay compensation sub-module is used to receive the synchronous timing data output by the space-time synchronization network module, including: The two-way timestamp protocol is used to measure the round-trip delay of the transmission link between the master node and the slave node. The delay fluctuation curve is fitted by the least squares method to generate the reference transmission delay amount; Based on the Kalman filter algorithm, the random delay deviation caused by network jitter is predicted. Combining the reference transmission delay amount, an interpolation function is constructed to resample the sensor data arriving asynchronously on the time axis for microsecond-level time alignment; The space calibration sub-module is connected to the time delay compensation sub-module and includes: The three-dimensional point cloud data of the monitoring area is obtained by a three-dimensional laser scanner with a resolution of 0.1 mm, and the actual point cloud coordinates of the sensor node installation position are extracted; The iterative closest point algorithm with normal vector constraint is used to register the actual point cloud coordinates with the theoretical coordinates in the design drawing, and calculate the coordinate offset of the sensor node; Based on the quaternion rotation matrix, the local coordinate system of each sensor node is converted to the global engineering coordinate system, and the spatio-temporal alignment data is output to the dynamic association model construction module.

[0033] The time delay compensation sub-module of the data calibration module receives the synchronous timing data output by the spatio-temporal synchronization network module. The two-way transmission link delay between the master node and the slave node is measured by the two-way timestamp protocol. The sending and receiving timestamps are recorded during the protocol interaction process, and the round-trip time difference of the path is calculated and the one-way transmission asymmetry error is excluded. Based on the continuous measurement data set, the least squares method is used to fit the delay fluctuation curve, and the reference delay component of the transmission link is extracted to eliminate the influence of the fixed propagation delay on time synchronization. For the random delay fluctuation caused by network jitter, a delay state prediction model is established by the Kalman filter algorithm. Combining the reference delay amount, a time axis interpolation function is generated to resample the sensor data arriving asynchronously, and the multi-source data is aligned to a unified time reference to achieve microsecond-level time synchronization accuracy.

[0034] The space calibration sub-module performs space coordinate system conversion on the time-calibrated data. The three-dimensional laser scanner performs panoramic scanning on the monitoring area along the predetermined scanning path to generate high-precision three-dimensional point cloud data containing the sensor node installation position. During the point cloud data preprocessing stage, environmental noise points are removed, and the point cloud cluster where the sensor node is located is segmented by the region growing algorithm to extract the actual three-dimensional coordinates of the node installation base. The iterative closest point algorithm with normal vector constraint is used for point cloud registration. During the algorithm iteration process, the normal vector direction of the sensor installation surface is introduced as a rigid transformation constraint condition, and the actual point cloud coordinates are least squares matched with the theoretical coordinates in the design drawing to calculate the coordinate offset and rotation angle. Based on the quaternion rotation matrix, a conversion model from the local coordinate system to the global coordinate system is constructed to eliminate the measurement error caused by the sensor installation position deviation, and the calibrated data with a unified spatio-temporal reference is output.

[0035] The time-delay compensation sub-module and the spatial calibration sub-module form a spatio-temporal joint calibration link. The data after time-axis resampling carries accurate time tags, and the spatial coordinate transformation model maps local measurement values to the global engineering coordinate system. The two cooperate to achieve spatio-temporal consistency calibration of multi-source heterogeneous sensor data. The calibrated data integrates triple information of timestamp, spatial coordinates, and physical quantity measurement values, providing an input data source with spatio-temporal traceability for the dynamic association model. The time interpolation and spatial transformation parameters used in the calibration process are recorded through calibration logs, supporting subsequent data traceability and inversion verification of model parameters.

[0036] Specifically, for the multi-sensor safety monitoring system of the pumped storage power station described in the present invention, the dynamic association model construction module includes: A fluid-structure interaction solution sub-module and a digital twin drive sub-module, where: The fluid-structure interaction solution sub-module is used to receive the spatio-temporally aligned data output by the data calibration module, including: Based on the coupled partial differential equations of the seepage pressure gradient tensor and the surrounding rock strain tensor, the dynamic association relationship between the seepage field and the stress field is calculated through an implicit iterative solver, and the coupled response data of seepage pressure and strain is generated; The digital twin drive sub-module is connected to the fluid-structure interaction solution sub-module and includes: Using the discrete element method to construct a three-dimensional digital twin model of the surrounding rock fracture propagation, mapping the real-time seepage pressure data to the permeability boundary conditions of the fracture surface, and simulating the non-linear relationship between permeability and strain during the fracture propagation process through the dynamic relaxation algorithm; According to the coupled response data of seepage pressure and strain, update the fracture network topology structure of the digital twin model, drive the digital twin model to dynamically evolve with a second-level time resolution, and output simulation verification data to the multi-dimensional data fusion module.

[0037] The fluid-structure interaction solution sub-module of the dynamic association model construction module receives the spatio-temporally aligned data output by the data calibration module. The spatio-temporally aligned data includes the seepage pressure field, the surrounding rock strain field, and the corresponding spatio-temporal coordinate information. The fluid-structure interaction model establishes a coupled control equation for the seepage pressure gradient tensor and the surrounding rock strain tensor based on the theories of porous media seepage mechanics and elastoplastic mechanics. The implicit iterative solver uses the Newton-Raphson method to handle the non-linear terms of the equation, updates the iteration step size through the Jacobian matrix, and calculates the dynamic interaction relationship between the seepage field and the stress field. An adaptive time step control strategy is introduced during the solution process, and the calculation step size is dynamically adjusted according to the seepage pressure change rate, and a dataset of coupled response of seepage and strain is output, including the evolution law of fracture permeability and the distribution characteristics of surrounding rock displacement.

[0038] The digital twin-driven sub-module constructs a three-dimensional numerical model of the surrounding rock fracture network based on the discrete element method. In the model initialization stage, geological exploration data and rock mechanics parameters are imported, and the initial contact stiffness and friction coefficient of the fracture surface are defined. The real-time seepage pressure data is mapped to the permeability boundary condition of the fracture surface through the spatial interpolation algorithm, and the dynamic relaxation algorithm is used to solve the contact force and motion equation of the fracture blocks, simulating the fracture propagation behavior caused by the fluctuation of seepage pressure. During the fracture propagation process, the non-linear relationship between permeability and strain tensor is characterized by the state-dependent stiffness decay model. The model updates the topological connection relationship and mechanical property parameters of the fracture network in real time according to the coupled response data output by the fluid-solid coupling solution sub-module.

[0039] The fluid-solid coupling solution sub-module and the digital twin-driven sub-module achieve collaborative calculation through the data interaction interface. The coupled response data of seepage and strain trigger the update of the boundary conditions of the digital twin model, and the fracture network evolution results inversely correct the permeability tensor term in the fluid-solid coupling equation, forming a two-way coupling iteration mechanism. The digital twin model outputs dynamic simulation data of fracture aperture, permeability distribution, and surrounding rock strain field with a second-level time resolution. After the simulation data and the measured data are aligned in space and time, they are transmitted to the multi-dimensional data fusion module for consistency verification. The model parameter update log records key variables such as the permeability attenuation coefficient and the contact stiffness correction amount, providing a basis for reverse verification for subsequent parameter optimization.

[0040] Specifically, for the multi-sensor safety monitoring system of the pumped storage power station described in the present invention, the multi-dimensional data fusion module includes: A spatio-temporal feature extraction sub-module and a dynamic weighted fusion sub-module, where: The spatio-temporal feature extraction sub-module is used to receive the simulation verification data output by the dynamic correlation model construction module and the spatio-temporal alignment data output by the data calibration module, including: A convolutional long short-term memory network improved by the multi-head attention mechanism, which inputs the time series data of seepage pressure and the spatial distribution data of displacement rate, extracts local deformation features through spatio-temporal convolutional kernels, and uses the multi-head attention mechanism to capture the correlation between pressure and displacement across sensor nodes, generating a spatio-temporal correlation feature matrix; The dynamic weighted fusion sub-module is connected to the spatio-temporal feature extraction sub-module, including: Calculating the model confidence weight according to the error rate between the simulation verification data and the measured data, and at the same time calculating the data reliability weight based on the historical stability index of the sensor nodes; Performing a normalized product operation on the model confidence weight and the data reliability weight to generate the dynamic weighted fusion coefficient of each sensor node; Performing multi-dimensional fusion on the spatio-temporal correlation feature matrix through a weighted decision function, and outputting the fusion decision data to the hierarchical warning module.

[0041] The spatio-temporal feature extraction sub-module of the multi-dimensional data fusion module receives the spatio-temporal alignment data output by the data calibration module and the simulation verification data output by the dynamic correlation model construction module. The spatio-temporal alignment data contains the temporal sequence of seepage pressure and the spatial distribution information of displacement rate, and the simulation verification data provides the numerical simulation results of fracture permeability evolution and surrounding rock strain field. The spatio-temporal feature extraction sub-module uses a convolutional long short-term memory network improved by a multi-head attention mechanism. The input layer of the network normalizes the temporal sequence data of seepage pressure, and at the same time converts the displacement rate data into a spatial grid distribution matrix. The spatio-temporal convolution kernel slides along the time axis and the spatial axis to extract the change gradient of seepage pressure and the spatial correlation features of displacement rate in the local area. The multi-head attention mechanism calculates the correlation weights of pressure and displacement between different sensor nodes in parallel to generate a spatio-temporal correlation feature matrix across dimensions.

[0042] The dynamic weighted fusion sub-module calculates the model confidence weight according to the residual analysis of the simulation verification data and the measured data. The residual analysis uses a sliding window statistical method to calculate the mean square error rate between the simulation data and the measured values at the corresponding spatio-temporal positions. The error rate is mapped to the model confidence weight coefficient after exponential smoothing processing. The historical stability index of the sensor node is generated by statistically calculating the signal packet loss rate, noise variance and long-term drift amount. The data reliability weight is quantified based on the index entropy method. The model confidence weight and the data reliability weight are fused through a normalized product operation to generate a dynamic weighted fusion coefficient matrix, and the weights of the nodes with low confidence or low reliability in the coefficient matrix are adaptively suppressed.

[0043] The weighted decision function performs an element-wise product operation on the spatio-temporal correlation feature matrix and the dynamic weighted fusion coefficient matrix. The fused feature vector is reduced in dimension through a fully connected layer to generate fusion decision data representing the system safety state. The fusion decision data integrates multi-dimensional information such as abnormal fluctuations of seepage pressure, sudden changes in displacement rate, and fracture expansion trends, and is output to the hierarchical early warning module to trigger threshold determination. During the feature extraction and weighted fusion process, the dimension of the spatio-temporal correlation feature matrix is consistent with the spatial distribution of the sensor nodes, and the dynamic update frequency of the weight coefficient matrix is synchronized with the data sampling period to ensure the timeliness and spatial resolution of the fusion result.

[0044] Specifically, for the multi-sensor safety monitoring system of the pumped storage power station described in the present invention, the hierarchical early warning module includes: An early warning trigger sub-module and a closed-loop parameter optimization sub-module, where: The early warning trigger sub-module is used to receive the fusion decision data output by the multi-dimensional data fusion module, including: Dividing the three-level early warning thresholds of seepage erosion, crack expansion and rock mass instability based on the structural safety margin model, where: The seepage erosion threshold is determined by calculating the ratio of seepage velocity to critical permeability through Darcy's law; The crack propagation threshold calculates the stress intensity factor at the crack tip based on the fracture mechanics J-integral model; The rock mass instability threshold determines the safety factor of the surrounding rock shear slip surface according to the Mohr-Coulomb criterion; Compare and fuse the decision-making data with the three-level thresholds in real time. When the continuous over-limit duration of the data exceeds the preset tolerance, trigger the corresponding-level audible and visual alarm signals; The closed-loop parameter optimization sub-module is connected to the early warning trigger sub-module, including: Use the clock synchronization error rate of the spatio-temporal synchronization network module and the root mean square of the prediction residuals of the dynamic correlation model construction module as the fitness function, and use the genetic algorithm for multi-objective optimization; Select elite individuals from the population through the tournament selection strategy. The crossover operator uses arithmetic recombination, and the mutation operator introduces Gaussian perturbation to iteratively search for the optimal combination of clock synchronization parameters and model coefficients; Feed the optimized parameters back to the clock synchronization sub-module of the spatio-temporal synchronization network module and the digital twin drive sub-module of the dynamic correlation model construction module in real time.

[0045] The early warning trigger sub-module of the hierarchical early warning module receives the fusion decision-making data output by the multi-dimensional data fusion module. The fusion decision-making data includes multi-dimensional characteristic quantities such as the abnormal index of seepage pressure, the displacement rate gradient, and the crack propagation rate. The early warning trigger sub-module sets three-level early warning thresholds based on the structural safety margin model. The seepage erosion threshold calculates the ratio of the seepage velocity in the monitoring area to the critical permeability of the rock mass through Darcy's law. When the ratio continuously exceeds the material durability threshold, it is determined as the seepage erosion risk; the crack propagation threshold is based on the fracture mechanics J-integral model to calculate the relative relationship between the stress intensity factor at the crack tip and the rock fracture toughness, and dynamically track the crack propagation trend; the rock mass instability threshold uses the Mohr-Coulomb criterion to analyze the ratio of the shear strength to the shear stress of the surrounding rock shear slip surface, and combines the safety factor threshold to determine the instability risk. The real-time monitoring data is compared with the three-level thresholds through a sliding window. When the continuous over-limit duration of the characteristic quantity exceeds the preset tolerance window, trigger the corresponding-level audible and visual alarm signals and generate early warning logs.

[0046] The closed-loop parameter optimization sub-module uses the spatio-temporal synchronization error rate and the model prediction residual as the optimization objectives. The spatio-temporal synchronization error rate is calculated through the timestamp deviation statistic between the master node and the slave nodes, and the model prediction residual is the root mean square error between the digital twin simulation data and the measured data. In the initialization stage of the genetic algorithm, population individuals containing clock synchronization parameters and fracture permeability attenuation coefficients are randomly generated. The fitness function is defined as the weighted sum of the spatio-temporal synchronization error rate and the model residual. The tournament selection strategy selects the top 10% of the elite individuals in terms of fitness from the population. The arithmetic recombination operator performs linear interpolation on the parameters of the elite individuals to generate offspring, and the Gaussian perturbation operator applies random perturbations to the offspring parameters to maintain population diversity. During the iterative optimization process, the search is terminated when the fitness value converges to a preset threshold, and the optimal parameter combination is output.

[0047] The optimized clock synchronization parameters are sent to the clock synchronization sub-module of the spatio-temporal synchronization network module through control messages to update the crystal oscillator drift correction coefficient; the fracture permeability attenuation coefficient is fed back to the digital twin driving sub-module of the dynamic correlation model construction module to correct the fracture surface permeability boundary condition. The parameter update triggers the model recalibration and data synchronization mechanisms, forming a closed-loop optimization link from early warning determination to parameter self-adaptation. The early warning logs and parameter update records are stored in the system database to support subsequent fault tracing and model iterative training.

[0048] The explanations of the technical features in the technical solution of the present invention are as follows: The redundant sensor node group adopts a triple redundancy configuration, which means that three sensors of the same type (such as strain, seepage pressure, vibration sensors) are deployed at each monitoring point, and the signals of the same physical quantity are collected through a synchronous trigger circuit. The stress concentration area of the penstock is a sensitive area for structural deformation. The triple redundancy design cross-verifies the data through the majority voting mechanism to improve the acquisition reliability. The sliding time window majority voting algorithm dynamically adjusts the window length to adapt to the sampling frequency. After normalization to eliminate the dimension difference, the Euclidean distance between two-by-two data is calculated. The dynamic threshold is generated based on historical data statistics. The continuously exceeding-limit data is determined as abnormal and excluded, and the remaining data is averaged and output. The adaptive wavelet threshold filtering algorithm selects the Daubechies or Symlet wavelet basis based on spectrum analysis. The improved SUREShrink threshold estimation introduces a frequency band energy weight factor to preferentially suppress high-frequency electromagnetic noise, and the reconstructed signal retains the time-domain waveform mutation and frequency-domain harmonic components.

[0049] The satellite time synchronization protocol (PTP / NTP) enables the master node to obtain the UTC reference clock from the Beidou system and distribute it to the slave nodes through a low-latency broadcast protocol. The Precision Clock Synchronization Protocol (IEEE 1588v2) establishes a two-way synchronization link between the master and slave nodes to compensate for the path delay asymmetry. The linear regression model analyzes the historical deviation data of the slave node crystal oscillator, predicts the frequency drift trend, dynamically corrects the value of the local clock counter, and achieves microsecond-level timestamp alignment. The cluster communication optimization sub-module divides the subnet according to geological units (such as the direction of the penstock and the surrounding rock fault). The cluster head node is elected through the self-organization protocol within the subnet. Time Division Multiple Access (TDMA) pre-allocates time slots, and the dynamic time slice adjustment algorithm compresses the time slot length and increases the density according to the channel load intensity, optimizing the data transmission timing.

[0050] The Two-Way Timestamp Protocol (TWTT) measures the round-trip delay between the master and slave nodes, and the least squares method fits the delay fluctuation curve to separate the fixed delay component. The Kalman filter algorithm predicts the random delay deviation caused by network jitter, and constructs an interpolation function to resample and align the time axis of asynchronous data. The 3D laser scanner acquires point cloud data with a resolution of 0.1 mm. The region growing algorithm segments the point cloud clusters of the sensor nodes. The Iterative Closest Point (ICP) algorithm with normal vector constraint introduces the normal vector direction constraint of the installation surface, registers the actual coordinates with the theoretical coordinates, and calculates the offset, which is then transformed into the global coordinate system through the quaternion rotation matrix.

[0051] The fluid-structure interaction model is based on the partial differential equation of the seepage pressure gradient and the surrounding rock strain tensor. The implicit iterative solver (such as the Newton-Raphson method) updates the iteration step through the Jacobian matrix, and calculates the dynamic correlation between the seepage field and the stress field. The Discrete Element Method (DEM) digital twin model initializes the fracture contact stiffness and friction coefficient. The real-time seepage pressure data is mapped as the permeability boundary condition of the fracture surface. The dynamic relaxation algorithm solves the block motion equation to simulate the fracture propagation. The non-linear relationship between permeability and strain is characterized by the stiffness decay model, and the fracture topology structure is updated with a second-level resolution.

[0052] The Conv-LSTM network improved by the multi-head attention mechanism normalizes the seepage pressure time series data. The displacement rate data is converted into a spatial grid matrix. The spatio-temporal convolution kernel extracts local deformation features. The multi-head attention calculates the cross-node pressure-displacement correlation weights in parallel, and generates a spatio-temporal correlation feature matrix. The dynamic weighted fusion coefficient generates the model confidence weight based on the mean square error rate of the sliding window statistical simulation data and the measured data, combines the historical packet loss rate of the sensor and the noise variance entropy method to calculate the data reliability weight, normalizes the product to generate the fusion coefficient matrix, and the weighted decision function fuses the feature matrix to output the decision data.

[0053] The seepage erosion threshold is calculated by the ratio of the seepage velocity to the critical permeability through Darcy's law. The crack propagation threshold is based on the J-integral model to calculate the ratio of the stress intensity factor at the crack tip to the rock fracture toughness. The rock mass instability threshold analyzes the safety factor of the surrounding rock shear slip surface according to the Mohr-Coulomb criterion. The genetic algorithm uses the clock synchronization error rate and the model residual as the fitness function, selects elite individuals by the tournament strategy, arithmetically recombines the linear interpolation parameters, and uses Gaussian perturbation to maintain the population diversity, outputs the optimal clock synchronization parameters (crystal oscillator drift correction coefficient) and the crack permeability attenuation coefficient, and feeds them back to the clock synchronization sub-module and the digital twin model to form a closed-loop optimization link.

[0054] In the specific implementation manner of the present invention, aiming at the complex hydraulic and mechanical coupling environment of the multi-sensor safety monitoring system of the pumped storage power station, high-precision spatio-temporal synchronization and the identification of the coupling mechanism of abnormal events are realized through multi-level collaborative technology. The sensor data acquisition module deploys three-redundancy sensor nodes in the stress concentration area of the penstock and the surrounding rock deformation sensitive area. Each monitoring point is equipped with three sensors of the same type to form a data acquisition unit, and the synchronous acquisition of stress, seepage pressure and vibration signals is realized through a synchronous trigger circuit. The sliding time window majority voting algorithm dynamically adjusts the window length according to the sensor sampling frequency, normalizes the three-way data and calculates the Euclidean distance. The dynamic threshold is generated based on the statistical characteristics of historical data. The continuously exceeding-limit data is determined as abnormal values and excluded, and the mean value of the remaining data is used as the effective output. The adaptive wavelet threshold filtering algorithm selects the Daubechies or Symlet wavelet basis function based on the signal spectrum analysis result, separates the high-frequency electromagnetic interference noise through the improved SUREShrink threshold estimation, and reconstructs the signal to retain the time-domain waveform mutation and the frequency-domain harmonic characteristics, generating preprocessed data.

[0055] The spatio-temporal synchronization network module obtains the reference clock signal through the Beidou satellite timing protocol. The master node analyzes the UTC time and broadcasts it to the slave nodes. The slave nodes establish a two-way synchronization link using the precision clock synchronization protocol. The linear regression model analyzes the historical deviation data of the slave node crystal oscillator in real time, predicts the cumulative drift amount and dynamically corrects the time stamp to achieve microsecond-level time alignment. The cluster communication optimization sub-module divides the monitoring subnet according to the geological unit. The cluster head node is elected within the subnet through the self-organizing network protocol, and fixed time slots are allocated based on the time division multiple access mechanism. The master node monitors the channel load intensity of the subnet in real time. When the data volume exceeds the preset threshold, the dynamic time slice adjustment algorithm compresses the single time slot length and increases the time slot density to optimize the data transmission timing. The data after time stamp correction is encapsulated into a standard format stream and transmitted to the data calibration module.

[0056] The data calibration module measures the transmission delay between the master node and the slave node using the bidirectional timestamp protocol, generates the reference delay amount by fitting the delay fluctuation curve with the least squares method, predicts the random deviation caused by network jitter with the Kalman filter, and constructs a time-axis interpolation function to align asynchronous data. The 3D laser scanner scans the monitoring area at a resolution of 0.1 mm to generate point cloud data. The region growing algorithm is used to segment the point cloud clusters of the sensor nodes. The iterative closest point algorithm with normal vector constraint is used to register the actual coordinates and the theoretical coordinates. After calculating the coordinate offset, it is converted to the global engineering coordinate system based on the quaternion rotation matrix. The dynamic association model construction module inputs the spatio-temporally aligned data into the fluid-structure interaction model, and the implicit iterative solver solves the coupled equations of the seepage pressure gradient and the surrounding rock strain tensor, and outputs the seepage pressure and strain response data. The discrete element method is used to construct a digital twin model of the surrounding rock fractures. The real-time seepage pressure data is mapped as the permeability boundary condition of the fracture surface. The dynamic relaxation algorithm simulates the non-linear relationship between permeability and strain, updates the fracture network topology, and outputs the simulation verification data at the second-level resolution.

[0057] The multi-dimensional data fusion module uses a convolutional long short-term memory network improved by the multi-head attention mechanism. The spatio-temporal convolutional kernel extracts the temporal gradient of the seepage pressure and the spatial distribution characteristics of the displacement rate. The multi-head attention mechanism calculates the cross-node pressure-displacement correlation weight to generate a spatio-temporal correlation feature matrix. The dynamic weighted fusion sub-module generates the model confidence weight based on the mean square error rate of the sliding window of the simulation data and the measured data, calculates the data reliability weight in combination with the historical packet loss rate and the noise variance index of the sensor, generates a dynamic fusion coefficient matrix through normalized product operation, and the weighted decision function fuses the feature matrix to output the decision data. The hierarchical early warning module sets the seepage erosion threshold by calculating the ratio of the seepage velocity to the critical permeability according to Darcy's law, sets the crack propagation threshold by calculating the stress intensity factor at the crack tip with the J integral model, sets the rock mass instability threshold by judging the safety factor of the shear slip surface of the surrounding rock with the Mohr-Coulomb criterion, and continuously triggers the audible and visual alarm when exceeding the limit. The genetic algorithm uses the clock synchronization error rate and the model residual as the fitness function, selects elite individuals with the tournament strategy, generates the optimal parameter combination through arithmetic recombination and Gaussian perturbation, feeds back to the clock synchronization sub-module to correct the crystal oscillator drift coefficient, and updates the fracture permeability boundary condition of the digital twin model to form a closed-loop optimization link.

[0058] The present invention solves the problem of mechanism identification deviation of multi-sensor monitoring networks in pumped storage power stations under complex coupling environments through a multi-level spatio-temporal synchronization architecture and a dynamic coupling model coordination mechanism. The sensor data acquisition module adopts a three-redundant node configuration and a sliding time window majority voting algorithm to synchronously collect stress, seepage pressure, and displacement signals, eliminates abnormal data through normalization processing and Euclidean distance calculation, and suppresses electromagnetic noise by combining adaptive wavelet threshold filtering to generate high-confidence preprocessed data. The spatio-temporal synchronization network module is based on the Beidou satellite timing protocol and a hierarchical clock synchronization architecture, predicts the crystal oscillator drift amount through a linear regression model and dynamically corrects the time stamp to achieve microsecond-level time alignment; the cluster communication optimization sub-module divides the monitoring subnet according to geological units, adopts a dynamic time slice adjustment algorithm to optimize the data transmission timing, reduces the probability of channel conflict, and ensures the real-time and integrity of data.

[0059] The data calibration module predicts the network transmission delay through a two-way time stamp protocol and a Kalman filter, constructs an interpolation function to resample asynchronous data on the time axis, and eliminates the timing deviation. The three-dimensional laser scanner obtains the point cloud data of the monitoring area, registers the actual coordinates and theoretical coordinates of the sensor nodes by using the iterative closest point algorithm with normal vector constraint, and converts them to the global engineering coordinate system based on the quaternion rotation matrix to achieve spatial calibration. The dynamic correlation model construction module inputs the spatio-temporally aligned data into the fluid-structure coupling model, calculates the dynamic correlation equation of seepage pressure gradient and surrounding rock strain through an implicit iterative solver, drives the discrete element method digital twin model to simulate the crack propagation process, and updates the permeability boundary conditions and crack network topology structure in real time to generate seepage-strain coupling response data, improving the simulation accuracy of the hydraulic and mechanical interaction mechanism.

[0060] The multi-dimensional data fusion module adopts a convolutional long short-term memory network improved by a multi-head attention mechanism to extract the cross-sensor spatio-temporal correlation features of seepage pressure and displacement rate, combines the simulation error rate and sensor stability index to generate a dynamic weighted fusion coefficient, and suppresses the interference of low-confidence data. The hierarchical early warning module sets three-level early warning thresholds based on Darcy's law, the J integral model, and the Mohr-Coulomb criterion, and uses the spatio-temporal synchronization error rate and model residual as fitness functions through a genetic algorithm to optimize the clock synchronization parameters and crack permeability attenuation coefficient, forming a closed-loop feedback mechanism. The parameter optimization results adjust the crystal oscillator drift correction coefficient of the clock synchronization sub-module and the boundary conditions of the digital twin model in real time, dynamically improving the system's identification ability of the coupling relationship of abnormal events.

Claims

1. Multi-sensor safety monitoring system for pumped storage power station, characterized in that: include: Sensor data acquisition module, spatiotemporal synchronization network module, data calibration module, dynamic correlation model building module, multi-dimensional data fusion module and hierarchical warning module; The sensor data acquisition module is used to collect stress, seepage pressure, displacement and vibration signals of key monitoring points through redundant sensor nodes, and perform noise suppression on the collected signals to generate pre-processed data; A spatiotemporal synchronization network module is used to receive the pre-processed data output by the sensor data acquisition module, realize microsecond-level time alignment of sensor nodes through a hierarchical clock synchronization architecture, dynamically adjust the data transmission timing based on clustered communication optimization, and output synchronous timing data; A data calibration module is used to receive the synchronous time series data output by the spatiotemporal synchronization network module, align the asynchronous data through a transmission delay compensation algorithm, and calibrate the global coordinate system of the sensor node in combination with a three-dimensional laser scanner to generate spatiotemporal alignment data; A dynamic correlation model building module is used to input the spatiotemporal alignment data output by the data calibration module into the fluid-solid coupling model, drive the digital twin model based on the discrete element method to simulate the dynamic response of seepage pressure and surrounding rock strain, and output simulation verification data; A multi-dimensional data fusion module, used to extract the spatiotemporal correlation features of the spatiotemporal alignment data and the simulation verification data, and generate fusion decision data through a weighted fusion algorithm; A hierarchical warning module is used to trigger a multi-level warning threshold according to the fusion decision data, and iteratively optimize the clock synchronization parameters of the spatiotemporal synchronization network module and the fluid-solid coupling model coefficients of the dynamic association model construction module based on the fusion decision data.

2. The multi-sensor safety monitoring system for a pumped storage power station according to claim 1 is characterized in that: The sensor data acquisition module comprises: Redundant sensor node group, abnormal data removal unit and adaptive noise suppression unit, wherein: The redundant sensor node group adopts a triple redundant configuration and is deployed in the stress concentration area of ​​the pressure steel pipe. Three sensors of the same type are set at each monitoring point to form a data acquisition unit. The three sensors synchronously collect the same physical quantity signal. The abnormal data elimination unit is connected to the redundant sensor node group, and eliminates abnormal data by a majority voting algorithm in a sliding time window, specifically: normalizing the data collected by the three sensors in the sliding time window, calculating the Euclidean distance between the two data, and if the distance between a certain data and the other two data exceeds a dynamic threshold, it is determined to be abnormal data and eliminated, and the mean value of the remaining data is used as a valid output; The adaptive noise suppression unit is connected to the abnormal data elimination unit, and adopts an adaptive wavelet threshold filtering algorithm to perform noise reduction processing on the effective output signal, including: dynamically selecting wavelet basis functions based on signal spectrum analysis results, separating electromagnetic interference noise frequency bands through improved SUREShrink threshold estimation, and retaining time domain waveform characteristics and frequency domain harmonic characteristics when reconstructing the signal, generating preprocessed data and outputting it to the space-time synchronization network module.

3. The multi-sensor safety monitoring system for a pumped storage power station according to claim 1 is characterized in that: The space-time synchronization network module includes a clock synchronization submodule and a cluster communication optimization submodule, wherein: The clock synchronization submodule is used to receive the pre-processed data output by the sensor data acquisition module, including: The master node obtains the reference clock signal from the Beidou satellite system through the satellite timing protocol and broadcasts the reference clock to all sub-nodes; Each sub-node establishes a time synchronization link with the master node through a precision clock synchronization protocol. At the same time, a linear regression model is used to predict the accumulated drift of each sub-node crystal in real time. The timestamp of the sensor data is dynamically corrected according to the drift, which is used for microsecond-level time alignment between sensor nodes. The cluster communication optimization submodule is connected to the clock synchronization submodule and includes: The monitoring area is divided into multiple geological unit subnets according to the geological structure characteristics, and the sensor nodes in each subnet form an independent communication cluster; In the subnet, a fixed data transmission time slot is allocated to each node through the time division multiple access mechanism. The master node monitors the channel load strength in real time. When the amount of node data in the subnet exceeds the preset threshold, a dynamic time slice adjustment algorithm is used to shorten the length of a single time slot and increase the time slot density to optimize the data transmission timing. The time-aligned and timing-optimized data is output as synchronous timing data and transmitted to the data calibration module.

4. The multi-sensor safety monitoring system for a pumped storage power station according to claim 1, characterized in that: The data calibration module includes a time delay compensation submodule and a space calibration submodule, wherein: The time delay compensation submodule is used to receive the synchronization timing data output by the time-space synchronization network module, including: The two-way timestamp protocol is used to measure the two-way transmission link delay between the master node and the sub-node, and the delay fluctuation curve is fitted by the least squares method to generate the benchmark transmission delay. The random delay deviation caused by network jitter is predicted based on the Kalman filter algorithm, and an interpolation function is constructed in combination with the benchmark transmission delay to resample the asynchronously arriving sensor data on the time axis for microsecond-level time alignment. The spatial calibration submodule is connected to the time delay compensation submodule and includes: The 3D point cloud data of the monitoring area is obtained through a 3D laser scanner, and the actual point cloud coordinates of the sensor node installation position are extracted; The iterative closest point algorithm with normal vector constraint is used to align the actual point cloud coordinates with the theoretical coordinates in the design drawing and calculate the coordinate offset of the sensor node; The local coordinate system of each sensor node is converted to a global engineering coordinate system based on the quaternion rotation matrix, and the spatiotemporal alignment data is output to the dynamic association model building module.

5. The multi-sensor safety monitoring system for a pumped storage power station according to claim 1 is characterized in that: The dynamic correlation model building module includes a fluid-solid coupling solution submodule and a digital twin driving submodule, wherein: The fluid-solid coupling solution submodule is used to receive the time-space alignment data output by the data calibration module, including: based on the coupled partial differential equation of the seepage pressure gradient tensor and the surrounding rock strain tensor, the dynamic correlation relationship between the seepage field and the stress field is calculated by an implicit iterative solver to generate the seepage pressure and strain coupling response data; The digital twin driving submodule is connected to the fluid-solid coupling solving submodule, including: using discrete element method to construct a three-dimensional digital twin model of surrounding rock crack expansion, mapping real-time seepage pressure data into crack surface permeability boundary conditions, and simulating the nonlinear relationship between permeability and strain during crack expansion through a dynamic relaxation algorithm; The fracture network topology of the digital twin model is updated according to the seepage pressure and strain coupling response data, the dynamic evolution of the digital twin model is driven, and the simulation verification data is output to the multi-dimensional data fusion module.

6. The multi-sensor safety monitoring system for a pumped storage power station according to claim 1, characterized in that: The multi-dimensional data fusion module includes a spatiotemporal feature extraction submodule and a dynamic weighted fusion submodule, wherein: The spatiotemporal feature extraction submodule is used to receive the simulation verification data output by the dynamic association model construction module and the spatiotemporal alignment data output by the data calibration module, including: The convolutional long short-term memory network improved by the multi-head attention mechanism is used to input the seepage pressure time series data and the displacement rate spatial distribution data. The local deformation features are extracted through the spatiotemporal convolution kernel, and the multi-head attention mechanism is used to capture the pressure and displacement correlation across sensor nodes to generate the spatiotemporal correlation feature matrix. The dynamic weighted fusion submodule is connected to the spatiotemporal feature extraction submodule, and includes: calculating the model confidence weight according to the error rate between the simulation verification data and the measured data, and calculating the data reliability weight based on the historical stability index of the sensor node; The model confidence weight and the data reliability weight are normalized and multiplied to generate the dynamic weighted fusion coefficient of each sensor node; The spatiotemporal correlation feature matrix is ​​multi-dimensionally fused through a weighted decision function, and the fused decision data is output to the hierarchical warning module.

7. The multi-sensor safety monitoring system for a pumped storage power station according to claim 1, characterized in that: The hierarchical warning module includes a warning triggering submodule and a closed-loop parameter optimization submodule, wherein: The early warning trigger submodule is used to receive the fusion decision data output by the multi-dimensional data fusion module, including: Based on the structural safety margin model, three-level warning thresholds are divided into seepage erosion, crack expansion and rock mass instability, among which: The seepage erosion threshold is determined by calculating the ratio of seepage velocity to critical permeability using Darcy’s law; The crack extension threshold is calculated based on the fracture mechanics J-integral model to calculate the stress intensity factor at the crack tip; The rock mass instability threshold is determined according to the Mohr-Coulomb criterion to determine the safety factor of the shear slip surface of the surrounding rock; Compare the fusion decision data with the three-level threshold in real time. When the data exceeds the limit for a continuous period exceeding the preset tolerance, the corresponding level of sound and light alarm signal is triggered; The closed-loop parameter optimization submodule is connected to the early warning trigger submodule, including: taking the clock synchronization error rate of the spatiotemporal synchronization network module and the prediction residual root mean square of the dynamic association model construction module as the fitness function, and using a genetic algorithm to perform multi-objective optimization; Elite individuals are selected from the population through the tournament selection strategy. The crossover operator uses arithmetic recombination, the mutation operator introduces Gaussian perturbation, and the optimal combination of clock synchronization parameters and fluid-structure interaction model coefficients is searched iteratively. The optimized parameters are fed back in real time to the clock synchronization submodule of the spatiotemporal synchronization network module and the digital twin driving submodule of the dynamic association model construction module.

Images