Abstract
Displacement monitoring data modeling is important for evaluating the performance and health conditions of concrete dams. Conventional displacement monitoring models of concrete dams decompose the total displacement into the water pressure component, temperature component, and time-dependent component. And the crack-induced displacement is generally incorporated into the time-dependent component, thus weakening the interpretability of the model. In the practical engineering modeling, some significant explaining variables are selected while the others are eliminated by applying commonly used regression methods which occasionally show instability. This paper proposes a crack-considered elastic net monitoring model of concrete dam displacement to improve the interpretability and stability. In this model, the mathematical expression of the crack-induced displacement component is derived through the analysis of large surface crack’s effect on the concrete dam displacement to improve the interpretability of the model. Moreover, the elastic net method with better stability is used to solve the crack-considered displacement monitoring model. Sequentially, the proposed model is applied to analyze the radial displacement of a gravity arch dam. The results demonstrate that the proposed model contributes to more reasonable explaining variables’ selection and better coefficients’ estimation and also indicate better interpretability and higher predictive precision.
1. Introduction
Structural safety monitoring of concrete dams depends on visual observation and monitoring data analysis with the purpose of, as early as possible, identifying eventual anomalous behavior which could deteriorate the dam. The analysis results can help the dam administrative department take corresponding measures to eliminate the impact of abnormal behavior [1]. Displacement monitoring data have been widely acknowledged as a significant source for the performance and health condition assessments of concrete dams [2] because it directly reflects the structural behavior. Therefore, various monitoring models have been built to explain and predict the dam displacement which is influenced by many factors such as hydrostatic pressure, concrete temperature, and time effect [3]. The hydrostatic-season-time (HST) model is commonly used to interpret dam deformation [4]. Wu et al. [5] discussed the factor selection of displacement statistical models for concrete dams. Rocha [6] introduced the multiple linear regression (MLR) model to analyze the deformation behavior of concrete dams. Based on the assumption that the residual error sequence of MLR models is uncorrelated and that multicollinearity does not exist among those explanatory variables, the optimal unbiased estimates of MLR model coefficients can be obtained using the least-squares method [7]. If some explanatory variables are linearly related, the stepwise least square regression (SLSR) method can be used to estimate model coefficients with multicollinearity considered. To date, HST models have been used to analyze the monitoring data of many dams, and various methods have been proposed to obtain the specific expressions and coefficients of practical models [8–16].
Traditional displacement monitoring model research studies focused on model calculation [17], but ignored key factor selection. However, many cracks exist in a lot of concrete dams which have been operating for years [18]. A mouth opening degree of a large surface crack varies regularly with the external loads, and the opening degree variations directly affect the concrete dam displacement [19, 20]. Crack-induced displacement is incorporated into the time effect component of those classical models, so the accuracy and interpretability of those models are weakened somewhat. It is necessary to analyze a large surface crack’s effect on concrete dam displacement, and a crack-considered monitoring model is yet to be built to improve interpretability.
The original theoretical displacement monitoring model of concrete dams conservatively includes as many explaining variables as possible to minimize the model deviation caused by the lack of important explaining variables. However, the actually applied model only needs those explaining variables which have a closer relationship with the explained response to improve its interpretability and predictive precision. Stepwise regression combined with AIC or BIC criterion has been generally used to select the optimal model. However, it has some inevitable drawbacks. Breiman [21] pointed out that the stepwise regression method is unstable when applied to select the actual model. And Fan et al. [22] indicated that a random error exists in the calculation process of stepwise regression and the theoretical nature of the method is difficult to study. In order to solve the above problems, Breiman [23] puts forward the nonnegative garrote method based on penalized least squares to get the actual interpretable model. Combining the nonnegative garrote method and bridge regression proposed by Frank [24], Tibshirani [25] proposed the Lasso (least absolute shrinkage and selection operator) method to select the explaining variables. Efron et al. [26] put forward the least angle regression algorithm to promote the practical application of the Lasso method. The Lasso method is being carried out in many fields [27–31], including statistics, engineering, biology, and information science.
The Lasso method excessively shrinks the coefficients, which have larger absolute values sometimes, resulting in a larger model error in some cases. In view of this, Fan and Liu [22] proposed the SCAD method, and Zou [32] presented the adaptive Lasso method to solve the deficient coefficient shrinkage problems in some situations. Zou and Hastie [33] introduced the quadratic penalty of coefficients into the Lasso method and put forward the another improved Lasso method called the elastic net method. The elastic net method can select significant explaining variables and estimate the corresponding coefficients simultaneously. Moreover, the method moderately shrinks the coefficients. The elastic net method with better stability is more suitable to solve the displacement monitoring model.
This article is organized as follows. In Section 2, large surface crack’s effect on the concrete dam displacement is analyzed theoretically and then the crack factor is introduced to develop a theoretical crack-considered monitoring model of concrete dam displacement. In Section 3, the elastic net method with better stability is used to solve the crack-considered displacement monitoring model. In Section 4, the proposed model is applied to analyze the radial displacement of a gravity arch dam, and the analysis results demonstrate better interpretability and higher predictive precision. Concluding remarks complete the paper in Section 5.
2. Theoretical Model Development
Traditional concrete dam displacement monitoring models divide the displacement response into the hydrostatic pressure component, temperature component, and time effect component, which can be expressed as follows:where represents the concrete dam displacement response, denotes the displacement caused by hydrostatic pressure, indicates the displacement caused by temperature, refers to the displacement caused by the time effect, and means the random error term.
Based on the dam engineering theory and engineering mechanics, the mathematical expressions for each component of the concrete dam horizontal displacement response in equation (1) have been deduced [6].
2.1. Hydrostatic-Pressure Component
The hydrostatic-pressure-induced horizontal displacement of any point in a concrete dam consists of three parts [34] which can be expressed as follows:
The first part comes from the dam body deformation caused by the reservoir water pressure. The second part is derived from the foundation deformation induced by reservoir water pressure. The third part is from foundation surface rotation caused by reservoir water gravity.
For the gravity dam, the cross section can be simplified as a triangular wedge with a vertical upstream surface and the dam can be simplified as a cantilever beam fixed on the foundation, as shown in Figure 1. The following explanation can be obtained through the mechanical relationship between the water level and dam deformation. Reservoir water pressure on the beam has a linear correlation with water depth. The water pressure to the dam body causes the horizontal displacement :where H is the upstream water level, h is the height of the dam, m is the downstream slope, d is the distance between the observation point and the dam crest, and are the elastic modulus and shear modulus of dam concrete, respectively, and is the water density.

The effect of reservoir water pressure is transferred by the cantilever beam to the foundation. The water pressure to the bottom of a dam causes the horizontal displacement :where and are elastic modulus and Poisson’s ratio of the dam foundation, respectively.
The topography and geological conditions of the reservoir area are complex, so a strict derivation for the mathematical expressions of is difficult. Here, vertical reservoir water pressure is regarded as acting uniformly on an infinite elastic body surface with the assumption that the reservoir bottom near the dam is horizontal and has a certain width. The horizontal displacement caused by the rotation of dam foundation is, approximately, as follows:where is the rotation angle of the dam foundation surface at the dam heel.
From the above analysis, it can be found that of the gravity dam is derived to have a linear relation with , , and . And is linearly related with and . For the arch dam, reservoir water pressure on the upstream surface is borne by both the horizontal arch and the vertical cantilever beam which are simplified in its structural analysis, so the water pressure distribution on the beam is nonlinear. Water pressure on the beam for the arch dam has a quadratic relation with water depth . Thus, of the arch dam has a linear relation with , , and , and is linearly related with and . In addition, is approximately considered to be proportional to . Therefore, the mathematical expression for the hydrostatic-pressure-induced horizontal displacement of any point in a concrete dam can be expressed as follows:where stands for the coefficients corresponding to the hydrostatic-pressure component and takes 3 and 4, respectively, for the gravity dam and the arch dam.
2.2. Temperature Component
The temperature component of concrete dam displacement is mainly caused by the temperature variation of dam concrete and bedrock, so the measurements of the thermometers buried in the dam body and foundation are selected as temperature factors. However, the number of thermometers embedded in many dams is small and some thermometers have been failed after years of operation. Besides, the dam temperature field tends to be quasi-steady over time, so the temperature of any point in the dam is approximately considered to be periodic [35]. Thus, the temperature component can be simply represented by some periodic functions as follows:where denotes the accumulated days from the initial measured day to the observation day, and represent the corresponding regression coefficients, and refers to the cycle parameter which takes 1 for a year cycle and 2 for half a year cycle, and so on.
2.3. Time Effect Component
The time effect component of concrete dam displacement comes from many aspects, including creep and plastic deformation of dam concrete and bedrock, irreversible deformation caused by long-term adverse loads, and self-grown volume deformation [14]. In general, time effect displacement changed dramatically during initial water impoundment and then gradually steadily over time. For the concrete dam operating for years, by analyzing the developing trend of displacement observations and the remaining value , the expression of can be determined reasonably as follows:where and indicate the corresponding regression coefficients and = and has the same meaning with that in equation (7).
2.4. Crack Component
The mouth opening degree of some large surface cracks, such as horizontal surface cracks and longitudinal ones, varies regularly with the external loads, including hydrostatic pressure and temperature. The opening degree variations of these cracks obviously affect the concrete dam displacement. Herein, we take a horizontal surface crack as an example to analyze the crack’s effect on the dam displacement.
Assuming the concrete dam body is rigid, without considering water pressure, Figure 2 shows the effect of large surface crack on the dam horizontal displacement. The rotation angles and are deemed approximately equal (); thus equations (9) and (10) can be approximately given:where and indicate the distances shown in Figure 2, represents the opening degree variation of the crack, and denotes the horizontal displacement of point A. For a specific point A, is the deemed constant; then, the crack-induced displacement can be approximately expressed as follows:where represents the opening degree of the ith large surface crack, stands for the regression coefficients corresponding to the crack component, and refers to the number of large surface cracks.

2.5. Theoretical Crack-Considered Monitoring Model of Concrete Dam Displacement
Integrating the above four components, including the hydrostatic-pressure component, temperature component, time effect component, and crack component, the theoretical crack-considered monitoring model of concrete dam displacement can be expressed as follows:where denotes the constant term, indicates the random error term, and the other symbols have the same meanings with those in the above equations.
3. Practical Model Solution
3.1. Elastic Net
Given a linear regression with standardized predictors and centered response values , for and , letting , the lasso estimate is defined by the following equation:
For all , the solution for is . Assuming without loss of generality that and hence omitting , the lasso solves the L1-penalized regression problem of finding to minimize , where is the penalty parameter.
The lasso is similar to ridge regression with constraint . However, ridge regression only shrinks the coefficients, whereas the lasso estimate does both variable selection and shrinkage in the above L1-penalized regression. The lasso and ridge regression can share a more general penalty of the form , where q takes 1 for the lasso and 2 for ridge regression. q takes the smaller value for the lasso (i.e., closer to subset selection which appears as q tends to zero). It is a convex problem which is attractive for computational purposes.
The Lasso method is actually to optimize a function with loss item and penalty item . The improved Lasso methods were given by changing the penalty function item. Among these improved lasso methods, the elastic net estimate is defined as follows:where and represent the penalty parameters and . takes 0 for the ridge regression estimate. The Lasso estimate is given when takes zero.
Asymptotics for Lasso-type estimators has been analyzed by Knight and Fu [36]. However, the Lasso method shrinks all the coefficients to the same degree resulting in a large model error, and the Lasso method does not have oracle properties. In comparison, improved lasso methods including the elastic net method commonly have oracle properties, that is, the optimal model obtained by elastic net has the following properties.
Letting , represents the number of elements in the set and . Assuming is the estimator of model coefficients obtained by the method , the method has oracle properties indicating that asymptotically meets the following two conditions:where Σ stands for the covariance of real model parameters. The oracle properties are also referred to as super-efficiency in the parameter estimation.
3.2. Practical Elastic Net Crack-Considered Monitoring Model of Concrete Dam Displacement
Combining with the elastic net method with better stability, moderate coefficients’ shrinkage, and oracle properties, the practically applied elastic net crack-considered monitoring model of concrete dam displacement is obtained as follows:where represents the elastic net estimator of the model coefficients, denotes the number of observations, and indicates the number of explaining variables.
Determination of penalty parameters and is critical in the model solution. In general, n-fo1d cross va1idation [37] is applied to determine penalty parameters. The following equations briefly introduce the n-fo1d cross va1idation.
Suppose represents the complete dataset, represents the training set, and represents the testing set, where n = l,2,…. For each , and n, the estimate of can be obtained from the training set . Then, the standard of n-fo1d cross va1idation can be given as follows:
The estimator of the penalty parameters can be obtained as follows:
Figure 3 shows how to solve the elastic net crack-considered monitoring model of concrete dam displacement.

4. Project Case Application
4.1. Displacement Model Development
The crest horizontal displacement of a concrete dam is analyzed by the elastic net crack-considered monitoring model. This dam is a concrete gravity arch dam, as shown in Figure 4, located in the east of China. Its maximum height is 76.3 m, and it has 28 sections. The dam was constructed in three phases. The concrete lift in Phase II was poured in such a hurry that its shrinkage distortion was restrained by the concrete in Phase I, resulting in a horizontal crack in the longitudinal direction at the interface. The crack extends from dam section No. 5 to No. 28 and is up to 5 m deep into the dam body. A large vertical crack in the longitudinal direction also exists in the dam crest. The dam crest radial displacement monitoring data in dam section No. 26 from May 8, 2013, to July 25, 2019, are analyzed. The monitoring data provided in the supplementary material “Monitoring data of dam section No. 26 pdf” were obtained from the pendulums buried in the dam, as shown in Figure 5 and Figure 6. The elevation in Figure 6 is in meters.



According to the elastic net crack-considered model, with the actual project situation considered, the radial displacement monitoring model of the dam crest is obtained as follows:where represents the opening degree of the vertical crack described above, indicates the opening degree of the horizontal crack described above, and the other symbols refer to the meanings of those variables in the aforementioned equations.
The recorded water level curve of the reservoir is shown in Figure 7. Figure 8 displays the opening degree curves of the two cracks mentioned above.


5. Results and Discussion
The radial displacement monitoring data from May 8, 2013, to December 31, 2018, were used to establish the monitoring model. Tables 1 and 2 show the results obtained by different models. In Table 2, R stands for the multiple correlation coefficient. Larger the multiple correlation coefficient, better the fitted effect. In Table 2, S means standard deviation. Smaller standard deviation represents better fitted effect. R and S are defined as follows:
Figures 9–12 give comparisons between measured displacement curves and fitted ones obtained by four different models. Displacement components extracted by the crack-considered elastic net model are given in Figure 13. Figures 14 and 15 show comparisons between fitted displacement residual curves obtained by different models. Comparisons of these tables and figures indicate the following analysis result.(1)The multiple correlation coefficient R (0.961) of the crack-unconsidered stepwise model is smaller than that (0.975) of the crack-considered stepwise model. The standard deviation S (0.218) of the crack-unconsidered stepwise model is larger than that (0.199) of the crack-considered stepwise model. For models obtained by the elastic net method, R (0.964) of the uncrack-considered model is smaller than that (0.978) of the crack-considered model and S (0.215) of the uncrack-considered model is larger than that (0.193) of the crack-considered model. These results show that the crack-considered models have a better fitting effect than the crack-unconsidered ones.(2)The values of the coefficient d1 are positive for both the crack-considered stepwise model and the crack-considered elastic net model, which means that the dam crest displacement is positively correlated with the vertical crack opening degree, that is, the opening of the vertical crack contributes to the downstream displacement of the dam crest. Meanwhile, the negative values of d2 obtained by the two crack-considered models indicate that the opening of the horizontal crack contributes to the upstream displacement of the dam crest. It can also be seen from Figure 13 that the crack components account for considerable proportions in the dam crest displacement. Therefore, it is necessary to extract the crack component separately to explain the variation behavior of the dam displacement, that is, the crack-considered model has better interpretability.(3)For crack-unconsidered models, R (0.961) and S (0.218) of the stepwise regression model are approximately equal with those (0.964 and 0.215) of the elastic net model. Meanwhile, the crack-considered models also have the same comparison result. These results indicate that the elastic net models have approximately equal fitting effect with the stepwise regression models.(4)The zero-value coefficient items of the crack-considered stepwise model are different with those of the crack-unconsidered stepwise model. For the elastic net models, the zero-value coefficient items of the crack-considered model are consistent with those of the crack-unconsidered model. The results reveal that the stability of the elastic net models is superior to that of the stepwise regression models.(5)For stepwise regression models, the zero-value coefficients of the crack-unconsidered model include all-time effect coefficients (c1 and c2), that is, the time effect is regarded not to contribute to the displacement. However, it is confirmed that the time effect component cannot be a zero value by analyzing the measured series and referring to general engineering experience. The time effect component obtained by the elastic net models is not a zero value. The results imply that the elastic net models have better interpretability than the stepwise regression models.







In order to evaluate the forecast effect of the monitoring model, the radial displacement monitoring data from January 1, 2019, to July 25, 2019, were used as predictive samples. Figures 16 and 17 give comparisons between measured displacement residual curves and predicted ones obtained by different models. Figure 18 displays a comparison between the measured displacement curve and the predicted one obtained by the crack-considered elastic net model. Table 2 shows the predictive effect (Rpre and Spre) of different models. Rpre and Spre are multiple correlation coefficient and standard deviation, respectively, representing the predicted goodness. Through comparison from the table and figures, the following analysis result can be obtained.(1)The predictive multiple correlation coefficient Rpre (0.961) of the crack-unconsidered stepwise model is smaller than that (0.975) of the crack-considered stepwise model. The standard deviation Spre (0.224) of the crack-unconsidered stepwise model is larger than that (0.201) of the crack-considered stepwise model. The elastic net models also have the same comparison result. These results show that the crack-considered models have a better predictive effect than the crack-unconsidered models.(2)For the crack-considered models, Rpre (0.975) of the stepwise regression model is approximately equal with that (0.978) of the elastic net model, while Spre (0.174) of the elastic net model is smaller than that (0.201) of the stepwise regression model. The results indicate that the elastic net models have a better predictive effect than the stepwise regression models.



6. Conclusions
This study was devoted to establish and solve the crack-considered elastic net monitoring model of the concrete dam displacement. The main conclusions are as follows:(1)Concrete dam safety operation and management need analysis and forecast on the dam displacement monitoring data. By inputting displacement monitoring data into theoretical monitoring models, the final actually applied monitoring and forecasting models can be obtained to interpret and predict the work behavior of concrete dams.(2)Crack-induced displacement in classical displacement monitoring models is incorporated into the time effect component so that the accuracy and interpretability of those models are weakened somewhat. In this paper, the crack factor was introduced into the model to be an independent crack component. Through the theoretical analysis of large surface crack’s effect on the concrete dam displacement, a mathematical expression of the crack-induced displacement component was derived, and the theoretical crack-considered monitoring model of the concrete dam displacement was built to improve the interpretability of the model.(3)The theoretical monitoring model of the concrete dam displacement conservatively includes as many explaining variables as possible to minimize the model deviation. However, the actually applied model only needs those explaining variables which have a closer relationship with the explained response to improve its interpretability and predictive precision. Conventional stepwise regression methods used in variable selection have instability in some cases. The elastic net method with better stability and moderate coefficient shrinkage was used to develop the practically applied crack-considered elastic net monitoring model of concrete dam displacement. The proposed model gives more reasonable factors and parameter estimation.(4)study on the crack-considered elastic net monitoring model of concrete dam displacement has been conducted tentatively in the paper, and a project case exhibits better interpretability and higher prediction accuracy of the proposed model. Furthermore, the thought and method of establishing the crack-considered elastic net model can also be introduced to the deterministic model and hybrid model of concrete dams.
Data Availability
The data used to support the findings of this study are from a large water conservancy project and not suitable to upload to the network. The data are included within the Supplementary Information files.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors thank the Dam Safety Research Group of Hohai University for the monitoring data. This research was supported by the Science and Technology Project of Zhejiang Water Resources Department, the Research Development Support Project of Zhejiang University of Water Resources and Electric Power (88106003136), National Natural Science Foundation of China (51739003), and National Key R&D Program of China (2018YFC0407104).
Supplementary Materials
The monitoring data of dam section No. 26. pdf provide the dam crest radial displacement monitoring data in dam section No. 26 from May 8, 2013, to July 25, 2019. (Supplementary Materials)