A spatiotemporal feature fusion concrete arch dam deformation prediction method based on multi-scale hierarchical clustering

Abstract

This invention discloses a method for predicting the deformation of concrete arch dams based on spatiotemporal feature fusion using multi-scale hierarchical clustering. The method includes the following steps: decomposing the dam deformation sequence into trend, periodic, and fluctuation terms, and constructing a new multi-scale similarity distance index; employing a hierarchical clustering method based on the Ward criterion, and using the constructed multi-scale similarity distance index to cluster deformation measurement points with similar deformation patterns; constructing a deep learning model TCT-CANet that fuses spatiotemporal features; using the constructed TCT-CANet model to predict the deformation of the concrete arch dam, and evaluating the model's predictive performance through various evaluation indicators; simultaneously analyzing the role of each module in capturing spatiotemporal features and improving prediction accuracy through ablation experiments; comparing the TCT-CANet model with various time-series prediction methods on all measurement points to verify the model's superiority and generalization ability. This invention constructs the prediction model TCT-CANet based on spatiotemporal feature fusion, achieving accurate prediction of concrete arch dam deformation.

Description

Technical Field

[0001] This invention pertains to methods for predicting the deformation of concrete arch dams, specifically involving a method for predicting the deformation of concrete arch dams based on the spatiotemporal feature fusion of multi-scale hierarchical clustering. Background Technology

[0002] As a crucial water conservancy project, the operational safety of concrete arch dams directly impacts flood control, water resource management, and the safety of life and property in downstream areas. Deformation monitoring is a vital tool in structural safety management, comprehensively reflecting the dam's stress response and operational status changes, and serving as a crucial basis for assessing the health status and risk level of arch dams. Therefore, constructing a stable, reliable, and highly accurate dam deformation prediction model is of great significance for improving engineering safety monitoring and early warning capabilities. Influenced by factors such as structural stress paths, material continuity, and environmental disturbances, adjacent measuring points of concrete arch dams often exhibit similar deformation patterns, while significant spatial differences exist between different regions. Currently, traditional methods often focus on modeling the temporal characteristics of individual measuring points, neglecting spatial synergistic effects. Therefore, accurately identifying and modeling these spatial distribution characteristics, effectively integrating spatiotemporal information, and constructing a prediction model with synergistic expressive capabilities have become key to improving the accuracy of deformation prediction. Summary of the Invention

[0003] The purpose of this invention is to provide a method for predicting the deformation of concrete arch dams based on spatiotemporal feature fusion using multi-scale hierarchical clustering. By using a hierarchical clustering method based on multi-scale similarity distance, deformation measurement points with similar deformation patterns are clustered and partitioned, and a prediction model TCT-CANet based on spatiotemporal feature fusion is constructed to achieve accurate prediction of the deformation of concrete arch dams.

[0004] To achieve the above objectives, the technical solution adopted by the present invention is as follows: A method for predicting the deformation of concrete arch dams based on spatiotemporal feature fusion using multi-scale hierarchical clustering includes the following steps: S1 decomposes the dam deformation sequence into trend, periodic and random fluctuation terms to construct a new multi-scale similarity distance index; The specific steps of S1 are as follows: S101, the variational mode decomposition method is used to decompose the arch dam deformation sequence into trend, periodic and random fluctuation terms; S102, Calculate the similarity distance index for trend items. The calculation formula is as follows: in, The trend slope Standard deviation, For the rate of change, xand y For different arch dam deformation sequences; S103, Calculate the similarity distance index for periodic terms. The calculation formula is as follows: in, The period length, For amplitude, For phase, x and y For different arch dam deformation sequences; S104, Calculate the similarity distance index for the random fluctuation term. The calculation formula is as follows: in, For fluctuation intensity, For the range of fluctuation, x and y For different arch dam deformation sequences; S105 integrates similarity distance indices for trend, periodic, and stochastic fluctuation terms, calculating a multi-scale similarity distance index. The calculation formula is as follows: in, This is a similarity distance indicator for trend items. The similarity distance index for periodic terms. The similarity distance index for random fluctuation terms. , and These are the weights corresponding to the trend term, the period term, and the random fluctuation term, respectively. S2 adopts a hierarchical clustering method based on the Ward criterion and uses the multi-scale similarity distance index constructed in S1 to cluster and partition deformation measurement points with similar deformation patterns. The clustering partitions include the bank slope areas at both ends of the arch dam, the middle area of ​​the arch dam, the lower middle area of ​​the arch dam, and the bottom of the arch dam; S3, Constructing a deep learning model TCT-CANet that integrates spatiotemporal features; The specific steps of S3 are as follows: S301, the extended factor set consisting of the environmental feature factor set and the deformation sequence of other measurement points in the partition where the target measurement point is located is used as the input dataset of the model; The set of environmental characteristic factors consists of monitoring data on upstream water level and temperature; S302, the TCN module captures the spatial dependencies between measurement points within the same deformation partition and outputs a spatial feature representation; S303, the Transformer module analyzes long-term dependencies and trends in time series and outputs time feature representations; S304 uses a Cross-Attention mechanism to fuse spatial and temporal feature representations. S305 performs dimensionality reduction on the fused features through an adaptive average pooling layer, and finally outputs the prediction results through a fully connected layer. S4 uses the constructed TCT-CANet model to predict the deformation of concrete arch dams and evaluates the predictive performance of the model through various evaluation indicators.

[0005] Preferably, the specific steps of S2 are as follows: S201, Initialization: Before starting clustering, each deformation measurement point is treated as an independent cluster; S202, Intra-cluster error assessment: Calculate the multi-scale similarity error between each measurement point in each cluster and its cluster center, and assess the internal compactness of the current cluster accordingly; S203, Overall Error Measurement: Under the current clustering structure, the intra-cluster errors of all clusters are summarized as an evaluation index of the overall clustering effect; S204, Two-class merging error calculation: For any two different clusters, estimate the impact of the merged cluster structure on the overall error; S205, Minimum Incremental Merging: Select the two clusters that result in the smallest increase in overall error and merge them. Repeat steps S203 and S204 until the preset number of clusters is reached.

[0006] Preferably, in step S4, when evaluating the predictive performance of the model using multiple evaluation metrics, the evaluation metric is the coefficient of determination. R 2 Root mean square error RMSE Mean absolute error MAE and maximum absolute error AE max The calculation formula is: in, n For the sample size, y i and These are the actual and predicted values ​​of the arch dam deformation, respectively. This represents the average value of the actual deformation.

[0007] The beneficial effects of this invention are as follows: By combining multi-scale hierarchical clustering and spatiotemporal feature fusion methods, this invention proposes a high-precision deformation prediction model for concrete arch dams, TCT-CANet. This method, through hierarchical clustering based on multi-scale similarity distance, can reasonably partition deformation measurement points with similar deformation patterns. Simultaneously, the model incorporates TCN, Transformer, and Cross-Attention mechanisms to fully exploit the spatiotemporal features in the deformation sequence, enhancing the ability to model spatiotemporal features. Compared to traditional prediction methods, it has stronger generalization ability and superior prediction performance, providing more accurate deformation predictions in various complex monitoring scenarios, and offering more effective technical support for dam safety monitoring and early warning. Attached Figure Description

[0008] Figure 1 The flowchart shows a method for predicting the deformation of concrete arch dams based on spatiotemporal feature fusion using multi-scale hierarchical clustering. Figure 2 This is a schematic diagram of the spatial distribution of the clustering results; Figure 3 The deformation curves of each measurement point within each cluster are shown. Figure 4 This is a diagram illustrating the overall construction of the spatiotemporal feature fusion model TCT-CANet. Figure 5 This is a schematic diagram of the TCN structure; Figure 6 This is a schematic diagram of the Transformer structure; Figure 7 This is a comparison chart of the evaluation indicators for various ablation experiment models. Detailed Implementation

[0009] 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 with reference to the accompanying drawings.

[0010] like Figures 1-7 As shown, a method for predicting the deformation of concrete arch dams based on spatiotemporal feature fusion using multi-scale hierarchical clustering includes the following steps: S1 employs Variational Mode Decomposition (VMD) to decompose the dam deformation sequence into three components: trend, periodicity, and fluctuation. For each component, a corresponding metric is defined, constructing a new multi-scale similarity distance index. This index will serve as an important basis for evaluating the similarity of measurement points in subsequent hierarchical clustering analysis, supporting the construction of deformation monitoring zones.

[0011] S101, using the variational modal decomposition method, with a mode number of 3, the arch dam deformation sequence is decomposed into three main modal components: trend term. T Periodic termsP and random fluctuation term F .

[0012] S102 defines the similarity distance index for trend items, and the features selected include the trend slope ( ), standard deviation ( ) and rate of change ( The formula for calculating trend characteristics is as follows: Formula 1 Formula 2 Formula 3 in, n The length of the deformed sequence. t i For a point in time, and It is divided into the average of time and trend terms. T 1 and T n These are the initial and final values ​​of the trend term, respectively.

[0013] Based on the above trend characteristics, calculate the similarity distance index of the trend items. The calculation formula is as follows: Formula 4 in, The trend slope Standard deviation, For the rate of change, x and y For different arch dam deformation sequences.

[0014] S103, Define the similarity distance index for periodic terms, with selected features including period length ( ), amplitude ( ) and phase ( The formula for calculating the periodic characteristics is as follows: Formula 5 Formula 6 Formula 7 in, It is an autocorrelation function. It's a time delay. p i For periodic terms, F This represents the Fourier transform, and argmax represents the position where the phase angle is taken.

[0015] Based on the above periodic characteristics, the similarity distance index of periodic terms is calculated. The calculation formula is as follows: Formula 8 in, The period length, For amplitude, For phase, x and y For different arch dam deformation sequences.

[0016] S104 defines a similarity distance index for random fluctuation terms, where the selected features include fluctuation intensity ( ) and fluctuation range ( The formula for calculating random fluctuation characteristics is as follows: Formula 9 Formula 10 in, f i These are the data values ​​for the fluctuation term. It is the average value of the random fluctuation term.

[0017] Based on the above characteristics of random fluctuations, the similarity distance index of random fluctuation terms is calculated. The calculation formula is as follows: Formula 11 in, For fluctuation intensity, For the range of fluctuation, x and y For different arch dam deformation sequences.

[0018] S105 integrates the similarity distances of trend, periodic, and stochastic fluctuation terms to calculate a multi-scale similarity distance index. The calculation formula is as follows: Formula 12 in, This is a similarity distance indicator for trend items. The similarity distance index for periodic terms. The similarity distance index for random fluctuation terms. , and These are the weights corresponding to the trend term, the period term, and the random fluctuation term, respectively.

[0019] To ensure a balanced contribution of each feature to the clustering results, an average weight allocation is adopted. This method simplifies the model's complexity while avoiding the excessive influence of a single feature on the clustering results.

[0020] S2 employs a hierarchical clustering method based on the Ward criterion, utilizing the multi-scale similarity distance index constructed in S1 to cluster and partition deformation measurement points with similar deformation patterns.

[0021] This similarity index serves as the basis for calculating the differences between measurement points, ensuring that the deformation sequences of measurement points within the same partition are highly consistent in multi-scale features, while significant differences exist between different partitions.

[0022] S201, Initialization: Before starting clustering, each deformation measurement point is treated as an independent cluster.

[0023] S202, Intra-cluster Error Assessment: Calculate the multi-scale similarity error between each measurement point in each cluster and its respective cluster center, and assess the internal compactness of the current cluster accordingly.

[0024] Calculate the sum of squared errors within the class ( ESS For arbitrary clustering of multiple deformation measurement points of concrete arch dams. C ,That ESS The value is calculated as follows: Formula 13 in, For clustering C The center point of a cluster is defined as the measurement point in the cluster that minimizes the sum of the multi-scale similarity indices with all other measurement points. D It is a multi-scale similarity distance index.

[0025] S203, Overall Error Measurement: Under the current clustering structure, the intra-cluster errors of all clusters are summarized as an evaluation metric for the overall clustering effect.

[0026] Calculate the total squared error (S&F) ESS If there is currently k Clusters C 1, C 2,…, C k Then all clusters ESS The sum is expressed as: Formula 14 in, For the first i The sum of squared errors of each cluster.

[0027] S204, Two-class merging error calculation: For any two different clusters, estimate the impact of the merged cluster structure on the overall error.

[0028] Calculate the combined values ​​of the two classes under all conditions. ESS Incremental. Considering merging clusters. Ci and C j Forming new clusters C t In the case of calculation ESS The increment value is: Formula 15 in, For the first i The sum of squared errors of each cluster, For the first j The sum of squared errors of each cluster, For the first i The cluster and the first j The sum of squared errors of the clusters after merging the clusters.

[0029] S205, Minimum Incremental Merging: Select the two clusters that result in the smallest increase in overall error and merge them. Repeat steps S203 and S204 until the preset number of clusters is reached.

[0030] minimize ESS Incremental merge clustering. Select increment value. ΔESS ij The two smallest clusters are merged, and S203 and S204 are repeated until the predetermined target number of clusters is reached.

[0031] The optimal number of clusters was determined using the Calinski-Harabasz index and silhouette coefficient. The spatial distribution of the final clustering results is shown in the attached figure. Figure 2 As shown in the attached figure, the deformation curves of the test points within each cluster are as follows. Figure 3 As shown.

[0032] Combined with appendix Figure 2 and attached Figure 3 It can be seen that the clustering results divide the monitoring points into five spatial zones. The locations of each zone on the dam body have clear distribution characteristics, reflecting strong geographical concentration and consistency in deformation patterns. Cluster 1 includes two measuring points, PL1-1 and PL7-1, located on the bank slopes at both ends of the arch dam. Its deformation amplitude is relatively small, exhibiting an overall stable deformation trend without significant periodic fluctuations. Cluster 2 comprises five measuring points: PL2-1, PL3-1, PL4-1, PL5-1, and PL6-1, concentrated in the middle of the dam crest. Its deformation curves exhibit significant periodic fluctuations, with an overall upward trend in deformation. Cluster 3 includes three measuring points: PL5-2, PL4-2, and PL3-2, mainly distributed in the central region of the arch dam. These deformation curves also exhibit obvious periodic fluctuations, but the fluctuation amplitude is significantly smaller than that of partition 2, while also showing certain trend changes. Cluster 4 contains three measuring points: PL5-3, PL4-3, and PL3-3, located in the lower part of the arch dam. Its deformation curve shows a slow upward trend overall, accompanied by small fluctuations. Cluster 5 includes eight measuring points: PL2-2, PL6-2, PL2-3, PL6-3, PL3-4, PL4-4, PL5-4, and PL4-5. These points are located at the bottom and sides of the arch dam. The deformation curves are generally stable with moderate periodic fluctuations.

[0033] S3. To achieve efficient prediction of concrete arch dam deformation, a deep learning model integrating spatiotemporal features, TCT-CANet, was constructed. This model combines a temporal convolutional network (TCN), a Transformer module, and a cross-attention mechanism. Its overall structure is shown in the attached figure. Figure 4 As shown.

[0034] S301, the input dataset of the model is an extended factor set consisting of the environmental feature factor set and the deformation sequence of other measurement points within the partition where the target measurement point is located.

[0035] The environmental characteristic factor set consists of environmental factors, specifically monitoring data of upstream water level and temperature. The extended factor set consists of deformation factors, specifically deformation sequences other than the predicted target within the same partition.

[0036] The use of this integrated data can not only effectively capture the complex spatiotemporal relationships between measurement points, but also reveal in depth the impact of environmental factors on deformation patterns.

[0037] The S302, TCN module is used to capture the spatial dependencies between measurement points within the same deformation partition and output a spatial feature representation.

[0038] This module consists of two TCN Blocks, each containing two layers of one-dimensional dilated convolutions (dilation factor). d =1, 2), with the kernel size set to 3. The first block outputs 8 channels of features, and the second block expands the number of channels to 16 to extract richer high-level features.

[0039] The structural diagram of the TCN module is shown in the attached figure. Figure 5 As shown, causal convolution is the core component of the TCN model, ensuring that the model can only use current and past information when predicting future values, and will not use future information, thus avoiding the leakage of future information.

[0040] While stacked causal convolutions can expand the receptive field, they require many layers to process long sequences, which can lead to gradient blurring and reduced efficiency. To address this, TCN introduces dilated convolutions, using a dilation factor... d =2i ( i (For the number of layers), achieving exponential growth of the receptive field, significantly expanding the perceptual range without increasing the number of layers. Calculate dilated convolution. The calculation formula is: Formula 16 in, d It is the expansion factor. k The kernel size represents the convolution of past states.

[0041] Appendix Figure 5 (a) shows the expansion factor. d The diagram illustrates dilated convolution with values ​​of 1, 2, and 4, where d=1 is the hidden layer, d=2 is the hidden layer, and d=4 is the output layer.

[0042] As network depth increases, the increase in model parameters can lead to vanishing or exploding gradients. Residual blocks can alleviate this problem by adding skip connections to the network, allowing gradients to propagate efficiently. (The output of the residual block is shown.) o The calculation formula is: Formula 17 Where Activation is the activation function. F ( x ) represents residual information.

[0043] Appendix Figure 5 (b) illustrates the structure of the residual block, which consists of two layers of dilated causal convolutions and a Rectified Linear Unit (ReLU) function. Weight normalization layers and dropout layers are used to improve the network's generalization ability. To ensure consistency between the input and output dimensions, a 1×1 convolution is ultimately used for adjustment. (See appendix) Figure 5 (c) shows an example of residual connections in TCN.

[0044] S303, the Transformer module is used to analyze long-term dependencies and trends in time series and output time feature representations.

[0045] To meet the input dimension requirements of the Transformer, the original deformed sequence is upsampled using a 1×1 convolutional layer. This module only uses the encoder part of the Transformer, and by stacking two encoder layers, it can effectively extract temporal dimension feature information at multiple depth levels.

[0046] The Transformer mainly consists of an encoder and a decoder, and its structure is shown in the attached figure. Figure 6As shown. The encoder is responsible for extracting feature representations of the input sequence, and the decoder generates the corresponding output based on these features. Each encoder and decoder consists of multiple stacked layers with the same structure, each layer containing two substructures: a multi-head self-attention mechanism and a feed-forward neural network (FFN). The core mechanism of the Transformer is the self-attention mechanism, which enables the model to dynamically focus on information at different positions when processing sequence data, thereby enhancing its ability to understand and model temporal features. The calculation expression of the self-attention mechanism is as follows: Formula 18 in, For self-attention mechanism, , and These are the weight matrices for the query, key, and value, respectively. , K and V These represent the query, key, and value, respectively. For the dimension of the key, softmax is the normalization operation.

[0047] To enhance the model's expressive power, Transformer introduces a multi-head attention mechanism. As a multi-head attention mechanism, The formula for calculating the attention mechanism for each head is as follows: Formula 19 Formula 20 in, i =1, 2, ..., h , i Indicates the number of heads. , , , For the corresponding trainable weight matrix, Q , K and V These represent the query, key, and value, respectively.

[0048] After each attention sublayer, the Transformer uses a feedforward neural network (FFN) to perform a non-linear transformation of the features. The FFN consists of two fully connected layers and an activation function. Its formula is as follows: Formula 21 in, It is a feedforward neural network. W 1 andW 2 is a trainable weight matrix. b 1 and b 2 represents the bias term. The ReLU activation function is used to introduce nonlinearity and increase the expressive power of the model.

[0049] Since the Transformer model itself does not depend on the order of the input sequence, positional encoding is introduced to inject positional information. Positional encoding is generated using sine and cosine functions, and its formula is as follows: Formula 22 Formula 23 in, P For the position encoding matrix, p The time sequence is fixed in position. i For the dimension of the sequence number, d model The total dimension of the sequence.

[0050] In this way, positional encoding introduces positional information into the model, ensuring that the Transformer can effectively handle sequence-dependent information.

[0051] The spatial and temporal features extracted by S304, TCN, and Transformer are input side-by-side into a cross-attention mechanism for integration. The fused features... The calculation formula is as follows: Formula 24 in, Temporal features extracted for Transformer K space and V space Spatial features extracted for TCN.

[0052] Through this structure, the model dynamically selects the most relevant spatial information based on changes in temporal characteristics, significantly improving the model's ability to handle complex spatiotemporal interaction patterns.

[0053] In S305, the features fused through the Cross-Attention mechanism undergo dimensionality reduction via an adaptive average pooling layer, and finally the prediction result is output through a fully connected layer. To ensure optimal performance of the TCT-CANet model, all hyperparameters are optimized using a random search strategy. The final optimal hyperparameter configuration for TCT-CANet is shown in Table 1.

[0054] Table 1 Parameter settings for the TCT-CANet model

[0055] S4 uses the constructed TCT-CANet model to predict the deformation of concrete arch dams, evaluates the predictive performance of the model through various evaluation indicators, and analyzes the role of each module in capturing spatiotemporal features and improving prediction accuracy through ablation experiments.

[0056] S401 evaluates the predictive performance of the model using multiple evaluation metrics, with the coefficient of determination being one of the metrics. R 2 Root mean square error RMSE Mean absolute error MAE and maximum absolute error AE max These four indicators are used to evaluate the prediction accuracy and performance of the model at different measurement points. The calculation formulas for the above four indicators are as follows: Formula 25 Formula 26 Formula 27 Formula 28 in, n For the sample size, y i and These are the actual and predicted values ​​of the arch dam deformation, respectively. This represents the average value of the actual deformation.

[0057] R 2 It is used to measure the goodness of fit of the model; the closer its value is to 1, the better the fit. RMSE , MAE and AE max Used to assess prediction error, the closer the value is to 0, the smaller the error.

[0058] S402 uses ablation experiments to design various model configuration schemes, analyzes the impact of each feature extraction module and feature fusion mechanism on prediction performance, and verifies the contribution of each module in spatiotemporal feature representation and accuracy improvement.

[0059] To evaluate the contribution of each module in the TCT-CANet model to the deformation prediction performance, four ablation experiments were designed. The specific model configuration is detailed in Table 2.

[0060] Table 2. Configuration scheme for ablation experiments

[0061] The specific method of the configuration scheme is as follows: only environmental factors are used as input; the Cross-Attention mechanism is removed, and spatiotemporal feature fusion is performed using ordinary fully connected layers; the Transformer module is removed; and the TCN module is removed.

[0062] It is important to note that in the experiments where TCN or Transformer was removed, the model was left with only a single feature extraction module, and the Cross-Attention mechanism could not function. Therefore, the models corresponding to these experiments were versions that only contained TCN or Transformer structures, not "TCN-CANet" or "Transformer-CANet". The evaluation metrics for each experimental design model are attached. Figure 7 As shown in the figure. The results show that TCT-CANet performs best in all metrics, verifying the synergistic effect of TCN, Transformer and Cross-Attention mechanisms, as well as the effectiveness of using deformed sequences as input.

[0063] S5 compares the TCT-CANet model with various time series prediction methods at all measurement points to verify the model's superiority and generalization ability.

[0064] The proposed spatiotemporal feature fusion model TCT-CANet was compared with the performance of various mainstream time series prediction methods on all measurement points. The comparison methods included deep learning models (TCN-LSTM, ConvLSTM, GRU, LSTM), machine learning models (SVM, RF), and classical statistical models (HST). The prediction performance of each model was evaluated in different clustering regions to verify the applicability and advantages of the TCT-CANet model in complex deformation scenarios.

[0065] The average prediction performance of each model in different clusters is shown in Table 3.

[0066] Table 3. Average prediction performance of different models in each cluster

[0067] The results show that TCT-CANet outperforms the comparison model in terms of fitting accuracy and error control, and maintains good prediction consistency and robustness under different deformation modes, fully demonstrating its wide applicability and generalization ability in complex deformation feature scenarios.

Claims

1. A method for predicting the deformation of concrete arch dams based on spatiotemporal feature fusion using multi-scale hierarchical clustering, characterized in that, Includes the following steps: S1 decomposes the dam deformation sequence into trend, periodic and random fluctuation terms to construct a new multi-scale similarity distance index; The specific steps of S1 are as follows: S101, the variational mode decomposition method is used to decompose the arch dam deformation sequence into trend, periodic and random fluctuation terms; S102, Calculate the similarity distance index for trend items. The calculation formula is as follows: in, The trend slope Standard deviation, For the rate of change, x and y For different arch dam deformation sequences; S103, Calculate the similarity distance index for periodic terms. The calculation formula is as follows: in, The period length, For amplitude, For phase, x and y For different arch dam deformation sequences; S104, Calculate the similarity distance index for the random fluctuation term. The calculation formula is as follows: in, For fluctuation intensity, For the range of fluctuation, x and y For different arch dam deformation sequences; S105 integrates similarity distance indices for trend, periodic, and stochastic fluctuation terms, calculating a multi-scale similarity distance index. The calculation formula is as follows: in, This is a similarity distance indicator for trend items. The similarity distance index for periodic terms. The similarity distance index for random fluctuation terms. , and These are the weights corresponding to the trend term, the period term, and the random fluctuation term, respectively. S2 adopts a hierarchical clustering method based on the Ward criterion and uses the multi-scale similarity distance index constructed in S1 to cluster and partition deformation measurement points with similar deformation patterns. The clustering partitions include the bank slope regions at both ends of the arch dam, the middle region of the arch dam, the lower middle region of the arch dam, and the bottom of the arch dam; S3, Constructing a deep learning model TCT-CANet that integrates spatiotemporal features; The specific steps of S3 are as follows: S301, the extended factor set consisting of the environmental feature factor set and the deformation sequence of other measurement points in the partition where the target measurement point is located is used as the input dataset of the model; The set of environmental characteristic factors consists of monitoring data on upstream water level and temperature; S302, the TCN module captures the spatial dependencies between measurement points within the same deformation partition and outputs a spatial feature representation; S303, the Transformer module analyzes long-term dependencies and trends in time series and outputs time feature representations; S304 uses a Cross-Attention mechanism to fuse spatial and temporal feature representations. S305 performs dimensionality reduction on the fused features through an adaptive average pooling layer, and finally outputs the prediction results through a fully connected layer. S4 uses the constructed TCT-CANet model to predict the deformation of concrete arch dams and evaluates the predictive performance of the model through various evaluation indicators.

2. The method for predicting the deformation of concrete arch dams based on spatiotemporal feature fusion using multi-scale hierarchical clustering as described in claim 1, characterized in that, The specific steps of S2 are as follows: S201, Initialization: Before starting clustering, each deformation measurement point is treated as an independent cluster; S202, Intra-cluster error assessment: Calculate the multi-scale similarity error between each measurement point in each cluster and its cluster center, and assess the internal compactness of the current cluster accordingly; S203, Overall Error Measurement: Under the current clustering structure, the intra-cluster errors of all clusters are summarized as an evaluation index of the overall clustering effect; S204, Two-class merging error calculation: For any two different clusters, estimate the impact of the merged cluster structure on the overall error; S205, Minimum Incremental Merging: Select the two clusters that result in the smallest increase in overall error and merge them. Repeat steps S203 and S204 until the preset number of clusters is reached.

3. The method for predicting the deformation of concrete arch dams based on spatiotemporal feature fusion using multi-scale hierarchical clustering as described in claim 1, characterized in that, In step S4, when evaluating the predictive performance of the model using multiple evaluation metrics, the evaluation metric is the coefficient of determination. R 2 Root mean square error RMSE Mean absolute error MAE and maximum absolute error AE max The calculation formula is: in, n For the sample size, y i and These are the actual and predicted values ​​of the arch dam deformation, respectively. This represents the average value of the actual deformation.

Images