Next Article in Journal
Possibilities of Using the Unitization Model in the Development of Transboundary Groundwater Deposits
Next Article in Special Issue
Soft Computing and Machine Learning in Dam Engineering
Previous Article in Journal
Bridge Pier Scour in Complex Environments: The Case of Chacao Channel in Chile
Previous Article in Special Issue
Risk-Informed Design of RCC Dams under Extreme Seismic Loading
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Relative Movements between Blocks Observed in a Concrete Dam and Definition of Thresholds for Novelty Identification Based on Machine Learning Models

1
National Laboratory for Civil Engineering (LNEC), Avenida do Brasil, 101, 1700-066 Lisbon, Portugal
2
CONSTRUCT-LESE, Department of Civil Engineering, Faculty of Engineering, University of Porto (FEUP), 4200-465 Porto, Portugal
3
Structural Health Monitoring Engineer, 1700 Lisbon, Portugal
*
Author to whom correspondence should be addressed.
Water 2023, 15(2), 297; https://doi.org/10.3390/w15020297
Submission received: 6 December 2022 / Revised: 4 January 2023 / Accepted: 9 January 2023 / Published: 11 January 2023
(This article belongs to the Special Issue Soft Computing and Machine Learning in Dam Engineering)

Abstract

:
Dam surveillance activities are based on observing the structural behaviour and interpreting the past behaviour supported by the knowledge of the main loads. For day-to-day activities, data-driven models are usually adopted. Most applications consider regression models for the analysis of horizontal displacements recorded in pendulums. Traditional regression models are not commonly applied to the analysis of relative movements between blocks due to the non-linearities related to the simultaneity of hydrostatic and thermal effects. A new application of a multilayer perceptron neural network model is proposed to interpret the relative movements between blocks measured hourly in a concrete dam under exploitation. A new methodology is proposed for threshold definition related to novelty identification, taking into account the evolution of the records over time and the simultaneity of the structural responses measured in the dam under study. The results obtained through the case study showed the ability of the methodology presented in this work to characterize the relative movement between blocks and for the identification of novelties in the dam behaviour.

1. Introduction

Structural monitoring involves observing a phenomenon or event and its impact on the structure. The purpose of the analysis and interpretation of the measurements gathered during an inspection is to enhance the conceptual understanding of the dam’s behaviour and aims to define, based on better parameters, eventual models. Once a model is built, the dam’s condition assessment is based on test hypotheses and scenario simulations supported by the monitoring data and the prediction of the structural behaviour in space and time. The procedure of providing information for assessing the dam’s structural condition is a form of critical analysis capable of reducing the intrinsic uncertainty related to the dam’s behaviour.
The assessment of the dam’s condition based on the information provided by the monitoring system is possible if the information is updated. Any abnormal behaviour can be readily identified, allowing for the implementation of an appropriate intervention or prevention measures. The assessment of the structural behaviour of the dam and its condition must be performed for each dam independently, even for dams of the same type, since they are influenced by several aspects, such as the heterogeneities in the dam’s foundation and the surrounding areas of the dam or the different loads (as a consequence of environmental or operational conditions) that the structure receives. For each concrete dam, different models can be used according to each purpose, to the existing knowledge about the actual structural behaviour, and the quality of information available for the characterization of the structure’s behaviour. The selection of the conceptual model to be used for dam surveillance activities must take into account: (i) the purpose of the analysis (safety assessment, prediction of deformations, interpretation of the recorded data from the monitoring system, or analysis of an accident or abnormal behaviour), (ii) the identification of the key factors of the physical problem, and (iii) the available geological and geotechnical information.
During a dam’s life, the performance and safety conditions are under continuous assessment due to the potential failure scenarios identified during the design phase and due to the other scenarios “suggested” by the observed behaviour through the analysis of relevant parameters (such as water level and temperature variations, among other things), as seen in Figure 1. Typically, these parameters will describe the loads or operating conditions to which the system is subjected, the materials of the structure, the materials of the structure’s foundations, and the structural response of the dam [1]. Thus, visual inspections and measurements collected by the structural monitoring system are fundamental for safety control activities. The main physical quantities measured through the monitoring system are seepage and leakage, uplift pressure, horizontal displacements, vertical displacements, relative displacements between blocks (contraction joint movements), and relative displacements in the rock mass foundation. The interpretation of the observed values is usually based on knowledge related to the physical and chemical phenomena that govern the structure and, whenever possible, is based on deterministic or data-based models. For decision-making, deterministic models are preferred, while in day-to-day activities, data-based models are the most commonly used.
The most common approaches for data-based models are the HST (hydrostatic, seasonal, time) and the HTT (hydrostatic, thermal, time) models, in which the effects of hydrostatic pressure, temperature, and time are considered additive effects, and their separation is valid [1,2,3,4,5]. The hydrostatic pressure is usually considered a polynomial function of the water height. However, in addition to the pressure from the water, dams in cold regions may be exposed to loads from an ice sheet [6]. The time effect is represented by polynomial functions or by functions of another type (exponential or s-shaped). The thermal effect is considered differently in each approach. In the HST approach, the temperature effect is represented by sinusoidal functions with a one-year period, which are a function of the day of the year only. In the HTT approach, the temperature effect is a function of the measured temperatures. A large number of publications about HST and HTT models can be consulted in the literature [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15].
Multiple linear regression is one of the most commonly used data-based models [3], showing good performance, mainly for case studies related to the prediction of horizontal and vertical displacements observed on a dam body. On the other hand, new models based on machine learning have been studied and published for about a decade with positive results; these are proposed mainly for the analysis of displacements. Proposed machine learning models could differ in their approach and advantages; such is the case of artificial neural network (NN) models, which principally allow the identification of non-linear relationships between the input quantities, if any. One main aspect to take into account in the case of NN models is that suitable generalization criteria, such as cross-validation criteria, must be adopted to avoid overfitting [5].
There are quantities, such as seepage, uplift pressure, and relative movements between blocks, in which regression models could be used with special care since the separation of effects may not be valid. For example, the behaviour of contraction joints depends on the state of stress installed. The variation observed differs for “open” or “closed” joints. Joint opening and closure are governed by normal stress criteria. The non-linear behaviour presented by the contraction joints can be explained by the non-resistance of the contraction joints to tension. As referred to by Hariri and Kianoush, modelling the joints (contraction, peripheral and lift joints) has an important role in both the static and the seismic analysis of concrete arch dams [16]. Some authors developed a non-linear joint element to represent the behaviour of vertical contraction joints in concrete dams [16,17,18,19,20,21]. The current knowledge in the field of machine learning allows the proper interpretation and prediction of these quantities, including anomaly detection [22,23], knowledge that is expected to grow in the following years.
The definition of thresholds for novelty identification related to movements in contraction joints is a natural step in the safety control activities of dams. For the identification of damage or gross measurement errors, the limits can be based on deterministic models or even on the maximum values observed in the dam life history affected by a multiplicative factor. However, identifying novelties earlier could be based on narrower limits resulting from data-based models. Nowadays, the definition of thresholds has been mainly based on considering a multiplicative factor associated with the standard deviation of residuals (the part not explained by the model adopted). This analysis allows the point-to-point identification of values that exceed established limits without any temporal or multidimensional context. Thus, two methodologies for threshold definition are presented in this work, one complementary to the traditional approach that can be used for the identification of novelties in a temporal context through the consideration of moving averages and the moving standard deviation of the residuals, and another, in a multidimensional context, through the use of the density-based spatial clustering of applications with noise (DBSCAN) method [24].
The methodology proposed is presented in Section 2. The case study, a concrete dam under exploitation, and the relative movements between blocks registered by the measurements of an automated monitoring system are described in Section 3. The results and final remarks are presented in Section 4 and Section 5.

2. Methodology Proposed for the Characterization of Relative Movements between Blocks and for Threshold Definition for Novelty Identification

2.1. Methodology

The priority for any threshold definition process is an adequate data-based model with a proper approach. The HTT approach was chosen in this study because of the daily variation observed in the quantity under study (relative movement between blocks). Enough information was obtained from the dam body’s thermometer and water level measurement devices. Regarding the machine learning method adopted, the authors selected the MLP-NN method. The methodology for the definitions of thresholds based on the residuals of the data-based models developed from the past history of the records (relative movements between blocks in the case study) is detailed ahead:
  • The adoption of the HTT formulation for the development of the prediction model requires the selection of the main physical quantities to be used. The water height and the concrete dam body temperature variations are the selected main inputs in the proposed model.
  • The separation of effects in the analysis of the relative movements between blocks may not be verified, so a model capable of mapping the interaction between the inputs is recommended. In this sense, an MLP-NN model was found to be suitable. Once the MLP-NN model is obtained, the influence of the main loads on the responses, as well as the analysis of the residuals, r, can be analysed.
  • The definition of thresholds for novelty identification based on the residuals can be defined: (i) based on one record at a time, (ii) based on a time window of records, (iii) and based on the simultaneity of the responses measured, as shown below.
  • A threshold definition based on one record at a time is the most used type of threshold in dam engineering and consists of the adoption of a multiplier n applied to the standard deviation of the residuals, s d . r . , as follows: n × s d . r . Values of n equal to 3 and 4 correspond to a confidence level of 99.73% and 99.99%, respectively, since residuals are assumed to follow a Gaussian distribution.
  • A threshold definition based on a time window of records of one physical quantity is based on the moving average of the residuals, m . a . r . , and the moving standard deviation of the residuals, m . s d . r . , adopting a time window of records over time. The size of the time window may depend on the case under study.
  • A threshold definition based on the simultaneity of the responses measured allows the identification of novelties related to a cluster of data (with m features) related to the m quantities considered in the analysis. D B S C A N is a method suitable for this approach, in which it is necessary to adopt, as an initial step, suitable values for Epsilon distance, ϵ and MinPoints, M i n P (Section 2.4). Hence, the residuals related to several physical quantities ( r 1 , r m ) , with m being the number of features measured simultaneously, are assessed to define clusters of the ones that present a behaviour far from that of the remaining data.
The proposed flowchart for the interpretation of the relative movements between blocks and for the operational threshold definition is presented in Figure 2.
The theoretical background related to the MLP-NN algorithm, the moving average, the moving standard deviation of the residuals, and the DBSCAN method are presented in the following subsections.

2.2. Multilayer Perceptron Neural Network Algorithm

Artificial neural networks have caught the attention of the scientific community since the 1990s [25] due to their ability to learn the pattern of structural behaviour in large infrastructures with a good capacity for generalization [5,14,26,27,28,29,30,31]. A NN computes a function of the inputs by propagating the computed values from the input neurons to the output neurons using different weights as an intermediate parameter [26]. A multilayer perceptron neural network is a feed-forward algorithm with neurons arranged in layers. It is the most widely used model for cognitive tasks such as pattern recognition and function approximation [5]. The specific architecture of feed-forward networks assumes that all nodes in one layer are connected to those of the next layer. The input layer transmits the data to the output layer going through a set of hidden layers that perform computations that include refining the weights between neurons over many input-output pairs to provide more accurate predictions.
The MLP-NN (Figure 3) learns by an iterative weight adjustment process that enables the correct learning using the training data so that, in a testing phase, it predicts the unknown data. Those weights w are located in each connection between the input layer x with N neurons and the hidden layer l with Q. The first layer receives the inputs, and the last layer produces the outputs. The middle layer is called the hidden layer. Within it, the information is constantly feed-forwarded from one layer to the next one having two associated values: the input value and the weight. In mathematics and programming, weights are shown in a matrix format W, where the number of columns contains the input dimensions (input values and weights), and the number of rows contains the output dimensions (hidden layer neurons Q). Another weight associated with the network is the bias b. Bias is attached to every layer in the network except for the input layer and carries neurons that try to account for unforeseen or non-observable factors. Activation functions are applied to each of the neurons in the hidden layer. The activation function will map the linear combinations of the inputs and weights to the following layer, in the case presented in Figure 3 from the input layer to the hidden layer and, in a further step, to the output layer of network L.
For regression problems, such as those found in the framework of SHM, the activation functions must be differentiable so that any non-linear models built from them can be derived with respect to their weights [32]. In this case, the activation function from the input layer to the hidden layer f can be the logistic sigmoid or the hyperbolic tangent function. The activation function from the hidden layer to the output layer g is a linear function.

2.3. Time-Window Threshold Definition Based on the Moving Average and Moving Standard Deviation of the Residuals

The detection of any deviation not identified by the prediction model between the observed and predicted data, is one of the main objectives of the residual analysis. A moving average control chart of the residuals ( m . a . r . ) is a type of memory control chart based on an unweighted moving average, and it is defined as
m . a . r . i = m . s u m . r . i w = j = i w + 1 i r j w , i w
where w is the width of the moving average at time i.
In order to describe the dispersion of residuals, the m . s d . r . is also presented in order to know how much the residuals vary or how spread they are along time, being defined as
m . s d . r . i = 1 w 1 j = i w + 1 i ( x j x ¯ j t o i ) 2
For these periods, the average of all observations within a time window of w size up to time i defines the moving average. Four different m . a . r . , related to 6, 12, 24, and 168 records, correspondent to a quarter, a half, an entire day, and a week, are potential options to be considered in this type of study. However, only the results of a m . a . r . and m . s d . r . with a time window of one day and a time step of one hour are presented in this work.
Once the m . a . r . and the m . s d . r . are established, the definition of a baseline of the structural behaviour under normal conditions is possible. The m . a . r . , in this sense, allows the identification of extreme values and trends along a time period. The m . s d . r . allows the identification of any increment of variability and/or randomness along the same time period. Those changes might suggest novelty.

2.4. Multivariate Threshold Definition Based on DBSCAN Algorithm

Density-based spatial clustering of applications with noise is a clustering algorithm proposed by Ester et al. [24]. It searches for “core objects”, points that contain a minimum of observations (MinPoints) within its neighbourhood (defined by an epsilon radius), including the core point itself. If a point is found outside of any of the core object’s neighborhood, it is considered noise [33]. Border points are points within reach of a core object without the minimum points in their neighbourhood to be considered core objects themselves (Figure 4). DBSCAN is known to be able to discover clusters with shapes other than linear and be robust enough to handle outliers and noise.
The clusters are created as follows: after uncovering a core object (i.e., a point with a high density of neighbours according to the parameters), DBSCAN starts a cluster with it and all its neighbours. All points within the Epsilon distance of the neighbours are then added to the cluster. This process continues until there are no more points within the distance. All the points that do not belong to a neighbourhood (i.e., to a cluster) are considered noise or outliers.

3. Case Study

3.1. The Feiticeiro Dam

The hydroelectric development of Baixo Sabor is composed of two schemes. The upstream scheme (termed the Baixo Sabor dam) is 12.6 km away from the Sabor river mouth. The downstream scheme (termed the Feiticeiro dam, Figure 5) is about 3.3 km far from the mouth of the Sabor river [35]. The Feiticeiro dam comprises a concrete gravity dam with an overflow-controlled spillway and a downstream stilling basin. The dam is 45 m high and has a crest length of 315 m, which is divided into twenty-two blocks. The total concrete volume is equal to 130,000 m 3 .
In accordance with the best technical practices, the monitoring system of the Feiticeiro dam aims at the evaluation of the loads, the characterisation of the rheological, thermal and hydraulic properties of the materials, and the evaluation of the dam’s structural response. The monitoring system of the Feiticeiro dam consists of several devices that make it possible to measure quantities such as the concrete and air temperatures, reservoir water level, seepage and leakage, displacements in the dam and in its foundation, joint movements, strains and stresses in the concrete, and pressures, among others [36,37,38].
The system used for the measurement of the reservoir water level comprised a high-precision pressure meter, which provides a record of the water height over time, and a level scale. The air temperature and humidity were measured in an automated weather station placed on the right-side bank approximately 100 m apart from the dam crest.
The concrete temperature was measured by 76 electrical resistance thermometers distributed across a dam thickness of several blocks. The location of the thermometers was defined taking into account the set of other electrical resistance devices (strain gauges, embedded jointmeters and strain gauges) that also allowed for the measurement of the concrete’s temperature.
Displacements were measured using an integrated system that included 3 pendulums, 8 rod extensometers, and geodetic observations. The relative movements between blocks were measured by superficial and embedded jointmeters.
The deformation of the concrete was measured with electrical strain gauges arranged in groups and distributed in radial sections, allowing the determination of the stress state through the knowledge of the deformation state and of the deformation law of the concrete. The quantities of drained and infiltrated water were measured individually in drains of the drainage system installed in the dam foundation and in weirs that differentiated the total quantity of water that flows in the drainage gallery in several zones of the dam. The drainage system comprised a set of 57 drains distributed over the drainage gallery with 3 drains per block. All the water extracted from drains and leakages was collected in 4 weirs.
The measurement of the uplift pressure at the foundation was performed by a piezometric network that comprised 26 piezometers. A real-time data acquisition system for the auscultation instrumentation (ADAS) was installed at Feiticeiro dam, allowing the measurement of the following physical quantities [39]:
  • Horizontal movements based on the pendulum method through the use of optical telecoordinometers (4–20 mA signal);
  • Vertical movements through rod strain meters with a position transducer of the vibrating wire type;
  • Discharges through level gauges with ultrasonic sensors (4–20 mA signal);
  • Pressures through piezometers with pressure transducers of the vibrating wire type;
  • Relative movements between blocks through 3D jointmeters with a position transducer of the vibrating wire type, located in the inspection galleries;
  • Relative movements between blocks through 1D jointmeters with a resistance transducer of the Carlson type, embedded in the concrete dam body;
  • Concrete strain through strain meters with resistance strain meter transducers of the Carlson type, embedded in the concrete dam body;
  • Concrete temperature through thermometers with resistance thermometer transducers of the Carlson type.
The responses under study are variations of the relative movements between blocks (at the higher level) measured by 7 automated jointmeters. The water level and the concrete temperature are measured by several thermometers located in the dam body, and they represent the main environmental loads studied in this case.
The mathematical calculations and the graphic representations presented ahead were supported by the R project software [40,41,42,43].

3.2. The Analysed Data

In this case study, the daily variation of the opening-closing movements between blocks measured at the 134.2 m level and through 3D jointmeters were analysed (jointmeter designed by BT3, BT5, BT7, BT8, BT10, BT12, and BT13).
The location of the jointmeter part of the automated monitoring system is shown in Figure 6. The data analysed correspond to a period between April 2017 and June 2022, with more than 38,660 records per variable. Measurements from the automated monitoring system are taken every hour. In the case of gaps in records, the values were estimated by interpolation between consecutive records. The manual and the automated measurements of the opening-closing movements were compared, showing the good performance of the two measurement systems (this comparison is not part of this study).
The samples regarding the relative movements between blocks were collected every two weeks for manual measurements and every hour for automated measurements, as seen in Figure 7. Signs (+) indicate opening movements, and signs (−) indicate closing movements.
Once the measuring devices to be studied in this work were installed in the downstream face of the dam body (a visiting gallery at the level of 134 m is nonexistent), they measured the relative movements between blocks in a zone that is strongly influenced by the daily variations of the temperatures, as shown in Figure 8, regarding the measurements between the 1st and the 8th of August 2018.
Among the different loads acting on concrete dams, it is typical to distinguish, as the most important loads for structures in normal operations, the hydrostatic pressure and the temperature variation. The time evolution of the reservoir water level is presented in Figure 9.
The temperature variations observed in the dam body were recorded through thermometers embedded in the dam body (blocks J7–J8, J11–J12, and J15–J16). The relative position of the thermometers are presented in Figure 6 and Table 1, and their records are those shown in Figure 10.
The observations presented in Figure 7, Figure 8, Figure 9 and Figure 10 will be used for the computation of the models presented in this work.
Based on the evolution of the time series of the joint movements, water level, and temperature variations, it is expected that the thermal effect is more significant when compared with the water level effect. Figure 11 presents an example of the temperature and the water height in the structural response for the case study. Figure 11 (left) shows the evolution of opening-closing movement measures in the jointmeter BT3 vs. the temperature measures in the thermometer T37. As expected, the movement between the blocks is in the closing direction when the temperature increases. Figure 11 (right) shows the evolution of the opening-closing movement measures in the jointmeter BT3 and the water height when the temperature measured in the thermometer T37 is 10 °C. The opening movement increases slightly with an increase in the water height.

4. Results and Discussion

4.1. Model Formulation, Construction, and Performance

The MLP-NN model based on an HTT approach considered, as an input layer, 15 parameters (representing the hydrostatic pressure—h, h 2 , h 4 , where h is the reservoir water level that can vary between 0 and 45 m—and the temperature effects measured at T10, T11, T15, T16, T31, T32, T33, T35, T62, T64, T65 and T66), seven responses at the output layer (representing the opening-closing movements between blocks measured at BT3, BT5, BT7, BT8, BT10, BT12 and BT13), and one hidden layer. Every neuron in the network is fully connected.
A hyperbolic tangent transfer function was chosen as the activation function for the hidden layer, and a linear activation function was chosen for the output layer. The generalized backpropagation delta learning rule algorithm was used in the training process. To find the optimum result (through the minimization of a cost function defined by the mean squared error), 5 initializations of random weights and a maximum of 1500 iterations were performed for each MLP-NN architecture.
A randomization of the learning set was previously carried out, making it possible to define the training set, the cross-validation set and the test set, with a number of examples equal to 65%, 15%, and 20%, respectively. The cross-validation was used as the stopping criteria. In each iteration, the performance for the training set is usually better than before, but if at any time the error for the cross-validation set increases, the NN model may lose its generalization capacity. The training stops when the error for the cross-validation set begins to increase, with a better generalization thus being ensured [5]. The test set was used as an auxiliary element that enabled us to carry out the quality evaluation of the MLP-NN model for the training set. In this case study, the network with the best performance was a 15-25-7 MLP-NN (less error for the cross-validation set). The results are represented in Figure 12, Figure 13 and Figure 14 and described in Table 2.

4.2. Threshold Definition for a Singular Record

The traditional approach adopted for threshold definition in regression models consists of adopting a multiplicative factor associated with the standard deviation of the residuals. Usually, multiplicative factors equal to 3 or 4, corresponding to a confidence level of 99.73% and 99.99%, respectively, are adopted. This approach assumes that the residuals follow a normal distribution (an assumption that is often not verified for real cases but that allows a good approximation). The same type of criterion can also be used in the MLR-NN model following the limits shown in Figure 15.

4.3. Threshold Definition for Novelty Identification Based on a Time Period Evolution of the Residuals

In case a novelty is identified, verifying if it is an isolated record is relevant in order to avoid the false internal warning since, usually, isolated novelties are associated with measurement errors. A critical aspect of dam safety control activities is the pattern recognition in the observed behaviour, including their expected evolution over time. With the exception of extreme load events and other time effects (such as the existence of internal expansion reactions in the concrete), the pattern in the observed behaviour is expected to be continuous over time. Thus, one way to identify novelties in the dam’s behaviour is to analyze the evolution of this behaviour (through the residuals) in a time window along time. Since daily variations resulting from daily temperature changes are observed, a moving window with 24 hourly records and with a time step equal to the measurement frequency (1 h) was adopted.
Figure 16 shows the evolution of the moving average of the residuals ( m . a . r . ) and of the moving standard deviation of the residuals ( m . s d . r . ) in a time window of one day over time.
The simultaneous analysis of the m . a . r . and m . s d . r . allows us to identify if the predicted values are moving away from the recorded values (i.e., if the residuals are increasing, meaning that the non-explained part of the model is increasing), and if the residuals increase along time, this could suggest a damage evolution. The definition of thresholds based on multiplicative factors associated with the standard deviation of the residuals is also proposed in this analysis. Similar to the single record analysis, thresholds based on multipliers of 3 and 4 for the standard deviation of the residuals are also proposed, although these values can be updated depending on the results, as seen in Figure 17. The priority in this kind of analysis is to identify periods in which the residuals and their evolution are higher; consequently, a deeper analysis of the data can be carried out.

4.4. Threshold Definition for Novelty Identification Based on Multivariate Data

Both of the previous analyses consider one physical quantity at a time. However, considering several quantities at the same time is relevant for an integral judgment. To address this consideration, a threshold definition based on the DBSCAN method is used. By definition, the residuals should follow a distribution similar to a normal distribution with a zero mean. Based on this premise, a criterion for the identification of spread records (in the n features dimension—seven in this case study) that are far from most of the records is defined. An analysis of the effect of the parameters Epsilon and MinPoints were performed as presented in Table 3. Values of Epsilon = 0.12 and MinPoints = 7 were adopted for this case study, as seen in Figure 18. These values were defined based on the specialist criteria considering the standard deviation of the residuals.
The dam behaviour in the time period under study presents a good performance, and there is no abnormal behaviour. The monitoring system records were used for defining the thresholds based on the three proposed procedures. Based on the thresholds defined and for future measurements, the records above the thresholds are the potential candidates to represent a novelty. They should be assessed in order to identify whether they are related to a measurement error or due to other situations.
The analysis of the monitoring data gathered from the horizontal displacements in a concrete dam’s behaviour and the adoption of data-based models to interpret the structural behaviour are common in day-to-day dam surveillance activities.

5. Conclusions and Final Remarks

The analysis of relative movements between blocks through data-based models is unusual due to the nonlinearity of the observed data, which cannot be represented in traditional linear regression models. In this case, a multilayer perceptron neural network (MLP-NN) model was developed to interpret the observed relative movements between blocks. The residuals (the part not explained by the model) resulting from the model were used in this work to define thresholds for novelty identification. The explanatory capacity of the residual analysis allows the definition of a more accurate baseline for the characterization of the structural behaviour. In this work, the analysis of the relative movements between blocks was performed based on hourly recorded measurements. The main inputs considered are the function of the water level and the temperatures measured in several thermometers spread along the dam body. The proposed procedures allow for the earlier detection of novelties through the analysis of the residuals of the MLP-NN prediction model, adopted for the interpretation of the relative movement between blocks, for (i) a singular record, (ii) a moving time period, and (iii) multivariate records. The main issues related to the definition of thresholds can be summarised as follows:
  • The definition of univariate thresholds based on the residuals, one for each quantity, is the first recommended step. This approach, through the use of multiplicative factors associated with the standard deviation of the residual, allows an easy novelty identification that, in some cases, could be related to measurement errors.
  • The second type of threshold proposed in this study can be defined by the adoption of a moving window with the moving average and the standard deviation of the residuals. This kind of approach allows the specialist to have immediate information concerning the evolution over time.
  • Finally, the third type of threshold proposed can be defined based on multiple factors related to the behaviour of the residuals for all the quantities considered in the model. For this approach, a strategy based on the DBSCAN algorithm, which proved to be suitable for multivariate analysis, was used.
To be effective, dam safety control activities must be considered an ongoing process. In this sense, applying the proposed procedures for operational issues in different dams and considering different quantities is suggested. It is important to highlight that each procedure has its singularities and challenges; the observed structural behaviour, the performance of the model adopted, and the quality and frequency of the measurement recorded need to be considered for the adoption of an adequate multiplicative factor for univariate thresholds and the adoption of adequate parameters for multivariate thresholds. The definition of the moving window size for calculating the moving average and the standard deviation depends on the measurement frequency of the records and the reservoir’s exploitation regime. Despite these challenges, the implementation of any of the proposed procedures in an automated monitoring system can adequately support dam surveillance activities.

Author Contributions

Conceptualization, J.M. and F.M.; methodology, J.M. and F.M.; software, J.M., F.M. and A.A.; validation, A.A., X.R. and J.P.S.; formal analysis, J.M., F.M., A.A., X.R. and J.P.S.; investigation, J.M., F.M. and A.A.; data curation, J.M. and A.A.; writing—original draft preparation, J.M. and F.M.; writing—review and editing, J.M., F.M., A.A., X.R. and J.P.S.; visualization, J.M. and F.M.; supervision, J.M., X.R. and J.P.S. All authors have read and agreed to the published version of the manuscript.

Funding

The contribution of the second author’s work was funded by the Portuguese Foundation for the Science and Technology (FCT), who funded this research through the grant number PD/BD/150407/2019.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Restrictions apply to the availability of these data. Data were obtained from Movhera—Hidroelétricas do Norte, S.A and Engie—Hidroelectricas do Douro, Lda, and are available from the authors with the permission of both of these entities.

Acknowledgments

The authors acknowledges the Movhera—Hidroelétricas do Norte, S.A and Engie—Hidroelectricas do Douro, Lda that provided the data for the procedures addressed in this paper, as well as LNEC through its research program RESTATE (0403/112/20970). The second author acknowledges the financial support of the Portuguese Foundation for the Science and Technology (FCT) through the grant PD/BD/150407/2019. The fourth author also acknowledges the financial support of the Base Funding—UIDB/04708/2020 of CONSTRUCT—Instituto de I&D em Estruturas e Construções, funded by national funds through FCT/MCTES (PIDDAC).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. ICOLD. Surveillance: Basic Elements in a Dam Safety Process; Bulletin Number 138; International Commission on Large Dams: Paris, France, 2009. [Google Scholar]
  2. Lombardi, G. Advanced data interpretation for diagnosis of concrete dams. In Structural Safety Assessment of Dams; CISM: Udine, Italy, 2004. [Google Scholar]
  3. Swiss Committee on Dams. Methods of Analysis for the Prediction and the Verification of Dam Behaviour. In Proceedings of the 21st Congress of the International Commission on Large Dams, Montreal, QC, Canada, 16–20 June 2003. [Google Scholar]
  4. Leger, P.; Leclerc, M. Hydrostatic, Temperature, Time-Displacement Model for Concrete Dams. J. Eng. Mech. 2007, 133, 267–277. [Google Scholar] [CrossRef]
  5. Mata, J. Interpretation of concrete dam behaviour with artificial neural network and multiple linear regression models. Eng. Struct. 2011, 33, 903–910. [Google Scholar] [CrossRef]
  6. Hellgren, R.; Enzell, J.; Ansell, A.; Nordström, E.; Malm, R. Estimating the Ice Loads on Concrete Dams Based on Their Structural Response. Water 2022, 14, 597. [Google Scholar] [CrossRef]
  7. Mata, J.; Tavares de Castro, A.; SÃ da Costa, J. Time frequency analysis for concrete dam safety control: Correlation between the daily variation of structural response and air temperature. Eng. Struct. 2013, 48, 658–665. [Google Scholar] [CrossRef] [Green Version]
  8. Salazar, F.; Toledo, M.A.; Oñate, E.; Morán, R. An empirical comparison of machine learning techniques for dam behaviour modelling. Struct. Saf. 2015, 56, 9–17. [Google Scholar] [CrossRef] [Green Version]
  9. Simon, A.; Royer, M.; Mauris, F.; Fabre, J. Analysis and Interpretation of Dam Measurements using Artificial Neural Networks. In Proceedings of the 9th ICOLD European Club Symposium, Venice, Italy, 10–12 April 2013. [Google Scholar]
  10. Ranković, V.; Grujović, N.; Divac, D.; Milivojević, N. Predicting piezometric water level in dams via artificial neural networks. Neural Comput. Appl. 2014, 24, 1115–1121. [Google Scholar] [CrossRef]
  11. Salazar, F.; Morán, R.; Toledo, M.Á.; Oñate, E. Data-Based Models for the Prediction of Dam Behaviour: A Review and Some Methodological Considerations. Arch. Comput. Methods Eng. 2017, 24, 1–21. [Google Scholar] [CrossRef] [Green Version]
  12. Granrut, M.; Simon, A.; Dias, D. Artificial neural networks for the interpretation of piezometric levels at the rock-concrete interface of arch dams. Eng. Struct. 2019, 178, 616–634. [Google Scholar] [CrossRef]
  13. Rico, J.; Barateiro, J.; Mata, J.; Antunes, A.; Cardoso, E. Applying Advanced Data Analytics and Machine Learning to Enhance the Safety Control of Dams. In Machine Learning Paradigms: Applications of Learning and Analytics in Intelligent Systems; Tsihrintzis, G.A., Virvou, M., Sakkopoulos, E., Jain, L.C., Eds.; Springer International Publishing: Cham, Switzerland, 2019; pp. 315–350. [Google Scholar]
  14. Mata, J.; Salazar, F.; Barateiro, J.; Antunes, A. Validation of Machine Learning Models for Structural Dam Behaviour Interpretation and Prediction. Water 2021, 13, 2717. [Google Scholar] [CrossRef]
  15. Li, M.; Wang, J. An Empirical Comparison of Multiple Linear Regression and Artificial Neural Network for Concrete Dam Deformation Modelling. Math. Probl. Eng. 2019, 2019, 7620948. [Google Scholar] [CrossRef]
  16. Hariri-Ardebili, M.A.; Kianoush, M.R. Integrative seismic safety evaluation of a high concrete arch dam. Soil Dyn. Earthq. Eng. 2014, 67, 85–101. [Google Scholar] [CrossRef]
  17. Ahmadi, M.T.; Izadinia, M.; Bachmann, H. A discrete crack joint model for nonlinear dynamic analysis of concrete arch dam. Comput. Struct. 2001, 79, 403–420. [Google Scholar] [CrossRef]
  18. Du, X.; Tu, J. Nonlinear seismic response analysis of arch dam-foundation systems—Part II opening and closing contact joints. Bull. Earthq. Eng. 2007, 5, 121–131. [Google Scholar] [CrossRef]
  19. Wang, J.; Zhang, C.; Jin, F.J. Nonlinear earthquake analysis of high arch dam–water–foundation rock systems. Earthq. Eng. Struct. Dyn. 2012, 41, 1157–1176. [Google Scholar] [CrossRef]
  20. Hariri-Ardebili, M.A.; Mirzabozorg, H.; Ghaemian, M. Pulvino and peripheral joint effects on static and seismic safety of concrete arch dams. Sci. Iran. 2014, 20, 1579–1594. [Google Scholar]
  21. Azmi, M.; Paultre, P. Three-dimensional analysis of concrete dams including contraction joint non-linearity. Eng. Struct. 2002, 24, 757–771. [Google Scholar] [CrossRef]
  22. Hellgren, R.; Malm, R.; Ansell, A. Performance of data-based models for early detection of damage in concrete dams. Struct. Infrastruct. Eng. 2021, 17, 275–289. [Google Scholar] [CrossRef]
  23. Salazar, F.; Conde, A.; Irazábal, J.; Vicente, D.J. Anomaly Detection in Dam Behaviour with Machine Learning Classification Models. Water 2021, 13, 2387. [Google Scholar] [CrossRef]
  24. Ester, M.; Kriegel, H.P.; Sander, J.; Xu, X. A density-based algorithm for discovering clusters in large spatial databases with noise. KDD 1996, 96, 226–231. [Google Scholar]
  25. Housner, G.W.; Bergman, L.A.; Caughey, T.K.; Chassiakos, A.G.; Claus, R.O.; Masri, S.F.; Skelton, R.E.; Soong, T.T.; Spencer, B.F.; Yao, J.T.P. Structural control: Past, present, and future. J. Eng. Mech. 1997, 123, 897–971. [Google Scholar] [CrossRef]
  26. Aggarwal, C. Neural Networks and Deep Learning; Artificial Intelligence; Springer International Publishing: Cham, Switzerland, 2018; pp. 389–411. [Google Scholar] [CrossRef]
  27. Shokri, M.; Tavakoli, K. A review on the artificial neural network approach to analysis and prediction of seismic damage in infrastructure. Int. J. Hydromechatron. 2019, 2, 178. [Google Scholar] [CrossRef]
  28. Fallahian, M.; Khoshnoudian, F.; Meruane, V. Ensemble classification method for structural damage assessment under varying temperature. Struct. Health Monit. 2018, 17, 747–762. [Google Scholar] [CrossRef]
  29. Sun, L.; Shang, Z.; Xia, Y.; Bhowmick, S. Nagarajaiah, S. Review of bridge structural health monitoring aided by big data and artificial intelligence: From condition assessment to damage detection. J. Struct. Eng. 2020, 146, 1–22. [Google Scholar] [CrossRef]
  30. Salehi, H.; Burgueño, R. Emerging artificial intelligence methods in structural engineering. Eng. Struct. 2018, 171, 170–189. [Google Scholar] [CrossRef]
  31. Liu, R.; Yang, R.; Zio, E.; Chen, X. Artificial intelligence for fault diagnosis of rotating machinery: A review. Mech. Syst. Signal Process. 2018, 108, 33–47. [Google Scholar] [CrossRef]
  32. Bishop, C.M. Neural Networks for Pattern Recognition; Oxford University Press: New York, NY, USA, 1995. [Google Scholar]
  33. Han, J.; Pei, J.; Kamber, M. Data Mining: Concepts and Techniques; Morgan Kaufmann Publishers: Burlington, MA, USA; Elsevier: Amsterdam, The Netherlands, 2011. [Google Scholar]
  34. Schubert, E.; Sander, J.; Ester, M.; Kriegel, H.P.; Xu, X. DBSCAN revisited, revisited: Why and how you should (still) use DBSCAN. ACM Trans. Database Syst. (TODS) 2017, 42, 1–21. [Google Scholar] [CrossRef]
  35. Cunha, J.; Vieira Soares, P. The Emergency Action Plans of the Baixo Sabor Dams. In Proceedings of the Second International Dam World Conference 2012, Lisbon, Portugal, 21–24 April 2015. [Google Scholar]
  36. LNEC. Downstream Dam of Baixo Sabor Hydroelectric Scheme. In Revision of the Monitoring Plan of the Dam in Conventional Concrete; Technical Report 67/2013—DBB/NO; LNEC: Lisbon, Portugal, 2013. (In Portuguese) [Google Scholar]
  37. LNEC. Downstream Dam of Baixo Sabor Hydroelectric Scheme. In Revision of the Automated Data Adquisition System Plan; Technical Report 244/2013; DBB/NO; LNEC: Lisbon, Portugal, 2013. (In Portuguese) [Google Scholar]
  38. LNEC. Feiticeiro Dam. In Inspection of the Structure and Monitoring System Carried Out in March 2020; Technical Report 149/2020—DBB/NO; LNEC: Lisbon, Portugal, 2020. (In Portuguese) [Google Scholar]
  39. Cunha, J.; Mata, J.; Gonzalo, O. Structural Safety Control of the Feiticeiro Dam Based on an Automated Data Acquisition System. In Proceedings of the 85th Annual Meeting of International Commission on Large Dams, Prague, Czech Republic, 3–7 July 2017. [Google Scholar]
  40. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
  41. RStudio Team. RStudio: Integrated Development Environment for R; R Studio, Inc.: Boston, MA, USA, 2019. [Google Scholar]
  42. Venables, W.N.; Ripley, B.D. Modern Applied Statistics with S, 4th ed.; Springer: New York, NY, USA, 2002; ISBN 0-387-95457-0. [Google Scholar]
  43. Peters, A.; Hothorn, T. ipred: Improved Predictors; R Package Version 0.9-11; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
Figure 1. Some parameters analysed for the assessment of dam safety and performance.
Figure 1. Some parameters analysed for the assessment of dam safety and performance.
Water 15 00297 g001
Figure 2. Proposed flowchart for the definition of thresholds for novelty identification.
Figure 2. Proposed flowchart for the definition of thresholds for novelty identification.
Water 15 00297 g002
Figure 3. MLP architecture where x represents the input layer, l represents the hidden layer, L represents the output layer, b represents the bias and w represents the weights, adapted from [5].
Figure 3. MLP architecture where x represents the input layer, l represents the hidden layer, L represents the output layer, b represents the bias and w represents the weights, adapted from [5].
Water 15 00297 g003
Figure 4. DBSCAN illustration with MinPoints = 4 and epsilon distance represented by circles. Core points, such as A, are represented in red, B and C, in yellow, are border points, and N is a noise point, represented in blue. From [34].
Figure 4. DBSCAN illustration with MinPoints = 4 and epsilon distance represented by circles. Core points, such as A, are represented in red, B and C, in yellow, are border points, and N is a noise point, represented in blue. From [34].
Water 15 00297 g004
Figure 5. Overall view of the Feiticeiro dam.
Figure 5. Overall view of the Feiticeiro dam.
Water 15 00297 g005
Figure 6. Feiticeiro dam. Thermometers and 3D jointmeters included in the automated monitoring system. Legend: MWL—Maximum reservoir water level.
Figure 6. Feiticeiro dam. Thermometers and 3D jointmeters included in the automated monitoring system. Legend: MWL—Maximum reservoir water level.
Water 15 00297 g006
Figure 7. Time series of the relative movements between blocks measured between April 2017 and June 2022.
Figure 7. Time series of the relative movements between blocks measured between April 2017 and June 2022.
Water 15 00297 g007
Figure 8. Time series of the relative movements between blocks measured between the 1st and the 8th of August 2018.
Figure 8. Time series of the relative movements between blocks measured between the 1st and the 8th of August 2018.
Water 15 00297 g008
Figure 9. Time series of the reservoir water level measurements between April 2017 and June 2022.
Figure 9. Time series of the reservoir water level measurements between April 2017 and June 2022.
Water 15 00297 g009
Figure 10. Temperatures in the dam body recorded between April 2017 and June 2022.
Figure 10. Temperatures in the dam body recorded between April 2017 and June 2022.
Water 15 00297 g010
Figure 11. Opening-closing mov. measured in BT3 vs. temperature measured in T37 (left). Opening-closing mov. measured in BT3 vs. water level for temperature measured in T37 equal to 10 °C (right).
Figure 11. Opening-closing mov. measured in BT3 vs. temperature measured in T37 (left). Opening-closing mov. measured in BT3 vs. water level for temperature measured in T37 equal to 10 °C (right).
Water 15 00297 g011
Figure 12. Time series of the fitted values and recorded values for the movements measured through 3D jointmeters.
Figure 12. Time series of the fitted values and recorded values for the movements measured through 3D jointmeters.
Water 15 00297 g012
Figure 13. Time series of the relative movements between block measurements and fitted values between the 1st and the 8th of August 2018.
Figure 13. Time series of the relative movements between block measurements and fitted values between the 1st and the 8th of August 2018.
Water 15 00297 g013
Figure 14. Fitted values vs. recorded values for the movements measured through 3D jointmeters.
Figure 14. Fitted values vs. recorded values for the movements measured through 3D jointmeters.
Water 15 00297 g014
Figure 15. Residuals of the movements measured through 3D jointmeters.
Figure 15. Residuals of the movements measured through 3D jointmeters.
Water 15 00297 g015
Figure 16. Residuals, m . a . r . , and m . s d . r . over time.
Figure 16. Residuals, m . a . r . , and m . s d . r . over time.
Water 15 00297 g016
Figure 17. m . s d . r . vs. m . a . r . for each quantity.
Figure 17. m . s d . r . vs. m . a . r . for each quantity.
Water 15 00297 g017
Figure 18. Novelty identification based on a multivariate analysis based on the DBSCAN method (adopting Epsilon = 0.12 and MinPoints = 7).
Figure 18. Novelty identification based on a multivariate analysis based on the DBSCAN method (adopting Epsilon = 0.12 and MinPoints = 7).
Water 15 00297 g018
Table 1. Relative location of the thermometers in the dam body.
Table 1. Relative location of the thermometers in the dam body.
BlockLevel (m)Position BlockLevelPosition
T10J7–J8122.1u.f.T11J7–J8122.11 m u.f.
T15J7–J8122.11 m d.f.T16J7–J8122.1d.f.
T31J11–J12112.0u.f.T32J11–J12112.01 m u.f.
T33J11–J12112.0othT35J11–J12112.0oth
T62J15–J16122.1u.f.T64J15–J16122.2oth
T65J15–J16122.11/2 th.T66J15–J16122.2oth
Notes: u.f.—upstream face, d.f.—downstream face, oth—other, th—thickness.
Table 2. MLP model performance parameters.
Table 2. MLP model performance parameters.
MeasurementsModel
MinMaxrMinrMaxsd R 2
[mm][mm][mm][mm][mm]%
BT3 O/C0.162.56−0.220.350.06099.81
BT5 O/C0.292.85−0.320.390.07299.79
BT7 O/C1.033.57−0.290.330.06899.92
BT8 O/C0.702.61−0.210.340.06399.86
BT10 O/C0.252.15−0.210.300.05699.77
BT12 O/C0.812.72−0.240.370.06799.85
BT13 O/C0.482.21−0.220.280.04899.86
Notes: r—residual, sd—standard deviation, R2—coefficient of determination.
Table 3. Number of potential novelties based on several Epsilon and MinPoints values.
Table 3. Number of potential novelties based on several Epsilon and MinPoints values.
MinPointsEpsilon
0.080.100.120.140.16
5799530
610212531
711614531
812622531
914022531
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mata, J.; Miranda, F.; Antunes, A.; Romão, X.; Pedro Santos, J. Characterization of Relative Movements between Blocks Observed in a Concrete Dam and Definition of Thresholds for Novelty Identification Based on Machine Learning Models. Water 2023, 15, 297. https://doi.org/10.3390/w15020297

AMA Style

Mata J, Miranda F, Antunes A, Romão X, Pedro Santos J. Characterization of Relative Movements between Blocks Observed in a Concrete Dam and Definition of Thresholds for Novelty Identification Based on Machine Learning Models. Water. 2023; 15(2):297. https://doi.org/10.3390/w15020297

Chicago/Turabian Style

Mata, Juan, Fabiana Miranda, António Antunes, Xavier Romão, and João Pedro Santos. 2023. "Characterization of Relative Movements between Blocks Observed in a Concrete Dam and Definition of Thresholds for Novelty Identification Based on Machine Learning Models" Water 15, no. 2: 297. https://doi.org/10.3390/w15020297

APA Style

Mata, J., Miranda, F., Antunes, A., Romão, X., & Pedro Santos, J. (2023). Characterization of Relative Movements between Blocks Observed in a Concrete Dam and Definition of Thresholds for Novelty Identification Based on Machine Learning Models. Water, 15(2), 297. https://doi.org/10.3390/w15020297

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop