Next Article in Journal
Typification for Façade Structures Based on User Perception
Next Article in Special Issue
Crustal and Upper Mantle Density Structure Beneath the Qinghai-Tibet Plateau and Surrounding Areas Derived from EGM2008 Geoid Anomalies
Previous Article in Journal
Past and Future Spatial Growth Dynamics of Chihuahua City, Mexico: Pressures for Land Use
Previous Article in Special Issue
Observational Analysis of Variation Characteristics of GPS-Based TEC Fluctuation over China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dam Deformation Monitoring Data Analysis Using Space-Time Kalman Filter

1
School of Geosciences and Info-Physics, Central South University, Changsha 410083, China
2
Center for Research in Geomatics, Laval University, Quebec, QC G1V 0A6, Canada
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2016, 5(12), 236; https://doi.org/10.3390/ijgi5120236
Submission received: 4 September 2016 / Revised: 9 November 2016 / Accepted: 23 November 2016 / Published: 8 December 2016
(This article belongs to the Special Issue Recent Advances in Geodesy & Its Applications)

Abstract

:
Noise filtering, data predicting, and unmonitored data interpolating are important to dam deformation data analysis. However, traditional methods generally process single point monitoring data separately, without considering the spatial correlation between points. In this paper, the Space-Time Kalman Filter (STKF), a dynamic spatio-temporal filtering model, is used as a spatio-temporal data analysis method for dam deformation. There were three main steps in the method applied in this paper. The first step was to determine the Kriging spatial fields based on the characteristics of dam deformation. Next, the observation noise covariance, system noise covariance, the initial mean vector state, and its covariance were estimated using the Expectation Maximization algorithm (EM algorithm) in the second step. In the third step, we filtered the observation noise, interpolated the whole dam unmonitored data in space and time domains, and predicted the deformation for the whole dam using the Kalman filter recursion algorithm. The simulation data and Wuqiangxi dam deformation monitoring data were used to verify the STKF method. The results show that the STKF not only can filter the deformation data noise in both the temporal and spatial domain effectively, but also can interpolate and predict the deformation for the whole dam.

1. Introduction

Water pressure, temperature, geological condition, construction quality, and other causes bring structural damage to dams, threatening the life of people living in the surrounding area. Dam deformation monitoring and its data processing are effective ways to assess dam safety and avoid possible dam failures [1]. Kalman filter, an effective data processing method for dynamic systems, has been widely used in deformation monitoring data processing [2,3], and improved Kalman filter methods, such as the Kalman filter for colored measurement noise and the robust Kalman filter, have been developed for deformation data analysis [4,5,6]. However, all these applications are based on single point monitoring data series, without considering the spatial correlation between monitoring points. So in some sense, neither the standard nor the improved Kalman filters suit the data analysis of dam deformation, which should be considered as a whole deformation body [7,8]. Furthermore, with the development of modern deformation monitoring technologies, such as Interferometric Synthetic Aperture Radar (InSAR), the observed deformation data can provide more spatio-temporal characteristics of deformation bodies. As such, it is imperative that new data analysis methods are developed to process the deformation data of all monitoring sites as a whole deformation body in both space and time domains.
The Space-Time Kalman Filter (STKF), proposed in the 1990s [9,10,11,12], is a spatio-temporal data processing method combining the Kalman filter and geostatistics models. The commonly used geostatistics model is the Kriging model, so the STKF is also called the Kriged Kalman Filter (KKF) [12]. To overcome the difficulties caused by large datasets and the modeling of spatio-temporal interactions in KKF, a dimension-reduced STKF was developed [13,14]. Distributed Kriged Kalman Filter (DKKF) [15] and Bayesian Kriged Kalman Filter (BKKF) [16,17,18] were proposed later. Current applications of STKF mainly focus on environmental issues. Mardia et al. (1998) employed KKF in dealing with the distribution of Sulphur dioxide in Leeds [12]. Wikle et al. (1999) validated the STKF model with a dataset of near-surface tropical winds [13]. Sahu et al. (2005) applied STKF to model atmospheric particulate distribution in New York in the USA and made a short-term prediction [16]. Lasinio et al. (2005) established a spatio-temporal rainfall model which includes seasonality and covariates based BKKF [17]. Al-Awadhi et al. (2012) used regression and BKKF to predict average hourly concentration of non-methane hydrocarbons in unmonitored areas of Kuwait [18]. Qing et al. (2012) adopted BKKF for short-term wind speed predictions [19]. However, STKF is quite rarely used in deformation data analysis. As such, this paper brings the STKF to the deformation analysis of dam.
In this paper, we first introduce the mathematical model of STKF. Then, the methods and algorithms to determine the Kriging spatial fields and parameters of STKF for the dam deformation data analysis are presented in detail. Then, we apply a STKF model to filter the noise, interpolate deformation data for the whole dam both in space and time domains, and make a short-term prediction. Finally, experiments with simulation data and real dam horizontal displacement monitoring data are conducted to verify the STKF model.

2. Space-Time Kalman Filter Model

In this part, the mathematical model of STKF is introduced in Section 2.1. Section 2.2 shows how to construct the spatial fields of STKF using the Kriging model. In Section 2.3, we discuss briefly the Expectation Maximization algorithm (EM algorithm) for parameter estimation. In the last section, we review the Kalman Filter recursion steps and methods of interpolating and predicting deformation in space and time.

2.1. Mathematical Model

Given location points { s = s 1 , s 2 , , s n } and time points { t = 1 , 2 , , m } , observation value L t ( s ) can be decomposed into mean and observed error components [14]:
L t ( s ) = Y t ( s ) + ε t ( s ) ,
where ε t ( s ) is the observed error components and Y t ( s ) is the mean component. Y t ( s ) can be expressed as a time-varying linear combination α ( t )   ( p × 1 ) of spatial fields h ( s )   ( p × 1 ) [12], which can be written as
Y t ( s ) = h 1 ( s ) α 1 ( t ) + h 2 ( s ) α 2 ( t ) + + h p ( s ) α p ( t ) = h ( s ) T α ( t ) ,
where ( ) T is the transpose operator and the vector α ( t ) is the state for Y t ( s ) , representing the dynamic time components. Generally speaking, p n , so the number of p is the key to reduce the dimension of the space domain.
Writing H = [ h ( s 1 ) T ; h ( s 2 ) T ; ; h ( s n ) T ] , L t = [ L t ( s 1 ) , L t ( s 2 ) , , L t ( s n ) ] T , and ε t = [ ε t ( s 1 ) , ε t ( s 2 ) , , ε t ( s n ) ] T , the substitution of (2) into (1) gives the observation Equation (3) of STKF.
L t = H α ( t ) + ε t
The state equation can be shown by
α ( t ) = Φ α ( t 1 ) + η ( t ) ,
where Φ is the state transition matrix, η ( t ) is the system error at time t . Suppose ε t ~ N ( 0 , Σ ε ) , η ( t ) ~ N ( 0 , Σ η ) , where N ( , ) is the normal distribution operator, Σ ε is the observation noise covariance and Σ η is the system noise covariance.
Equations (3) and (4) form the STKF model. According to the model, matrices H , Φ , Σ ε , Σ η and the initial state α ( 0 ) ~ N ( a ^ 0 | 0 , P 0 | 0 ) are essential to implement the STKF recursion. The spatial fields H are built by spatial Kriging fields which is the key to take spatial correlation into consideration. Parameter θ = { a ^ 0 | 0 , P 0 | 0 , Φ , Σ ε , Σ η } can be estimated by the EM algorithm. In the following sections, we will present the method to obtain these matrices in detail.

2.2. Spatial Fields H

Spatial fields H are built by spatial Kriging model. At time point t the universal Kriging model is as follows:
Y t ( s ) = f T ( s ) β t + ς t ( s ) ,
where f ( s ) = [ f 1 ( s ) , f 2 ( s ) , .. f q ( s ) ] T is the given space trend fields, whose elements are a function of the coordinate of s . Commonly used space trend fields are constant, linear and quadratic trends. The subscript q is determined by the chosen trend. β t is the coefficient of f ( s ) , and ς t ( s ) is the local space variation after removing space trend. Assume ς t ( s ) is stationary, and it can be described by the spatial covariance (or semi-variogram) model [14,16]. An empirical semi-variogram model can be obtained from the observations to describe ς t ( s ) . First, calculate the detrended data D t at time point t by ordinary least squares,
D t = Y t F ( F T F ) 1 F T Y t ,
where ( ) 1 is the inverse operator, D t = [ D t ( s 1 ) , D t ( s 2 ) , , D t ( s n ) ] T , Y t = [ Y t ( s 1 ) , Y t ( s 2 ) , , Y t ( s n ) ] T , and F = [ f ( s 1 ) T ; f ( s 2 ) T ; ; f ( s n ) T ] is the chosen trend matrix. Here, assuming F does not change with time, the empirical semi-variogram can be estimated by average semi-variogram,
γ ^ ( d ) = 1 2 m t = 1 m [ D t ( s ) D t ( s + d ) ] 2 ,
where ( ) is the sum operator. Choose a semi-variogram model to fit γ ^ ( d ) , and spatial covariance can be obtained by σ ς ( s i , s j ) = C γ ^ ( s i s j ) , where C is the sill value of semi-variogram model and is the distance operator.
After choosing the space trend fields and fixing the semi-variogram function, the Kriging predictor can be written as follows [20,21]:
Y ^ t ( s ) = f ( s ) T A Y t + σ ς ( s ) T B Y t ,
where σ ς ( s ) = [ σ ς ( s , s 1 ) , σ ς ( s , s 2 ) , , σ ς ( s , s n ) ] T . Matrices A and B have the fixed form as
A = ( F T Σ ς 1 F ) 1 F T Σ ς 1 ,
B = Σ ς 1 Σ ς 1 F A ,
where Σ ς is a spatial covariance matrix with ( Σ ς ) i j = σ ς ( s i , s j ) , and B is known as the bending energy matrix [12,22], which is used to describe the scale of local space variation after removing space trend. For dimension reduction, consider the spectral decomposition of matrix B ,
B = U D U T , B u i = d i u i ,
where columns of U = [ u 1 , u 2 , , u n ] are eigenvectors of B , and diagonal elements of D = diag ( d 1 , d 2 , .. , d n ) are corresponding eigenvalues, and diag ( ) is the diagonal matrix operator. It is easy to verify that B F = 0 , which means columns of F can be thought as the eigenvectors of B with corresponding zero eigenvalues d 1 , d 2 , .. , d q . Obviously, smaller eigenvalues show larger scale space variations, while larger eigenvalues describe the local spatial variations [16].
The eigenvectors u i form a set of orthogonal basis, so the mean vector Y t can be represented as a linear combination of the eigenvectors u i . Indeed, suppose that Y t = i = 1 n c i , t   u i ,   β ^ t = A Y t and Equation (8) can be rewritten as
Y ^ t ( s ) = f ( s ) T   β ^ t + i = q + 1 n c i , t d i σ ς ( s ) T u i = h ( s ) T [ β ^ t c t ]
where c t = [ c q + 1 , t , .. , c n , t ] T . If n is very large, we only take the first p space variations to reduce dimension, which results spatial fields in the form of dimensionality reduction (13). Parameter p can be obtained while the ratio e = i = 1 p d i i = 1 n d i exceeds a certain percentage (usually 90% or 95%) [18], which explains the proportion that the selected spatial variations take in the total spatial variation.
h i ( s ) = f i ( s ) , for   i = 1 , 2 , , q h j ( s ) = d j σ ς ( s ) T u j , for   j = q + 1 , , p
It is worth noting that h i ( s ) = f i ( s ) is the trend field used to describe the global trend. h j ( s ) = d j σ ς ( s ) T u j is the principal field used to grasp local spatial variation. If ignoring the principal field, STKF will reduce to a spatio-temporal trend surface model [16], and if continuously taking h i ( s ) = f i ( s ) = 1 , STKF will reduce to a standard Kalman model. Also, if only taking observation equation of STKF into consideration, STKF will reduce to a universal Kriging model.

2.3. Parameters Estimation

The maximum likelihood estimation (MLE) is an effective way to estimate parameter θ = { a ^ 0 | 0 , P 0 | 0 , Φ , Σ ε , Σ η } . However, from the logarithm of joint likelihood function (14), we cannot maximize directly the logarithm of joint likelihood function log ( l ) because of the unknown distribution of state vector α ( t ) .
log ( l ) = 1 2 log | P 0 | 0 | 1 2 ( α ( 0 ) α ^ 0 | 0 ) T P 0 | 0 1 ( α ( 0 ) α ^ 0 | 0 ) m 2 log | Σ η | 1 2 t = 1 m ( α ( t ) Φ α ( t 1 ) ) T Σ η 1 ( α ( t ) Φ α ( t 1 ) ) m 2 log | Σ ε | 1 2 t = 1 m ( L t H α ( t ) ) T Σ ε 1 ( L t H α ( t ) ) ,
where | | is the determinant operator and log ( ) is the logarithm operator.
The EM algorithm provides an iterative method to maximize log ( l ) , which consists of Expectation Step (E step) and Maximization Step (M step) [23]. Suppose it is in the ( r + 1 ) t h iteration, then the main steps of EM algorithm are as follows:
  • Use Kalman smoother to estimate the unknown state parameter α ( t ) with respect to the ( r ) t h iterated value θ r .
  • E step: calculate the conditional expectation G ( θ r ) = E ( log ( l ) | L 1 , L 2 , .. , L m ) of log ( l ) under the estimated α ( t ) distribution in step 1, where E ( ) is the expectation operator.
  • M step: maximize G ( θ r ) , which yields the newly iterated value θ r + 1 .
  • Replace θ r with θ r + 1 , and repeat steps 1, 2, and 3 until the logarithm of joint likelihood function log ( l ) or the innovations form [23] stop increasing.

2.4. Denoising, Space-Time Interpolation, and Prediction

As matrix H and parameter θ are estimated, Kalman filter recursion can be conducted. It should be noted that the STKF model needs no preprocessing, such as the interpolation of missing data, which can be dynamically calculated while filtering the noise. Suppose at time point t , the observation equation is as follows:
( L t 1 L t 2 ) = ( H 1 H 2 ) α ( t ) + ( ε t 1 ε t 2 ) ,
where L t 1 , H 1 , and ε t 1 are the observed data, spatial fields, and errors of observed data, respectively, and L t 2 , H 2 , and ε t 2 are the missing data, spatial fields, and errors of missing data, respectively. In order to implement Kalman recursion, setting L t 2 and H t 2 to zero matrix [23], the estimation of α ( t ) can be obtained by
α ^ t | t = E ( α ( t ) | L 1 , L 2 , , L t ) = α ^ t | t 1 + G t ( L t H α ^ t | t 1 ) ,
and its mean-squared-prediction-error matrix
P t | t = E [ ( α ^ t | t α ( t ) ) ( α ^ t | t α ( t ) ) T ] = P t | t 1 G t H P t | t 1 ,
where the Kalman gain matrix G t is given by
G t = P t | t 1 H T ( H P t | t 1 H T + Σ ε ) 1 .
The one-step-ahead-forecast quantities are given by
α ^ t | t 1 = Φ α ^ t 1 | t 1 ,
P t | t 1 = Φ P t 1 | t 1 Φ T + Σ η .
Short-term prediction involves optimal prediction of Y t ( s ) from   L 1 , L 2 , , L m , where t ϵ { m + 1 , m + 2 , } . Based on the one-step-ahead-forecast Equations (19) and (20), the optimal prediction estimator α ^ t | m of α ( t ) is
α ^ t | m = E ( α ( t ) | L 1 , L 2 , , L m ) = ( i = m + 1 t Φ ) α ^ m | m ,
with its mean-squared-prediction-error matrix
P t | m = ( i = m + 1 t Φ ) P m | m ( i = m + 1 t Φ ) T + Σ η + i = m + 1 t 1 { ( j = m + 1 t Φ ) Σ η ( j = m + 1 t Φ ) T }  
where Π is the multiplication operator. The optimal prediction of Y t ( s ) can be calculated by Y ^ t ( s ) = h ( s ) T α ^ t | m .
After obtaining α ^ t | t , the optimal estimator of Y t ( s ) at any location point s and time point t of interest can be calculated by
Y ^ t ( s ) = h ( s ) T α ^ t | t ,
where h ( s ) can be calculated by Equation (13).
Lastly, a flow diagram of the overall procedure of STKF is displayed in Figure 1.

3. Simulation Experiment

To simulate the horizontal displacement observation data of a dam, we assumed that the horizontal displacement has a quadratic curve trend in space and it changes in a sinusoidal trend in time. The horizontal displacement of a dam in space and time can be described by
g ( s , t ) = ( 0.0013 s 2 0.04 s 2.7 ) sin ( 0.005 π t 2 0.00125 π ) + ε .
In Equation (24), g is the displacement, s is the position of the monitoring points on the dam, t is the time, and ε is the observation noise with a standard deviation of 0.1 mm. Suppose s = 0 ( m ) + 20 k ( m )   , k = 0 , 1 , 2 , , 22 , for 23 points in total, and t = 1 ( d a y ) + 1 j ( d a y ) ,   j = 0 , 1 , 2 , , 2900 , then the horizontal displacement data series (1 day sampling rate) of 23 monitoring points have been simulated.
As mentioned in Section 2.2, in order to acquire the spatial fields H , we had to first get the trend matrix F . Since the simulation experiment constructs a known space-time curve, the trend field is known as F = [ s 2 , s , 1 ] . To get the principal fields, it was necessary to choose a suitable spatial semi-variogram model. Commonly used models are the exponential model, Gaussian model, spherical model, and Matèrn family model [20,21]. We chose the spherical model in this experiment as it can maximally show the spatial variation and is the most widely used spatial semi-variogram model [24]. The rest parameter θ = { a ^ 0 | 0 , P 0 | 0 , Φ , Σ ε , Σ η } was estimated by the EM algorithm, which is described in Section 2.3.
In the simulation experiment, we used all data to filter the noise firstly. Secondly, in order to verify the effectiveness of the spatio-temporal interpolation using STKF, we randomly sampled one point as the unobserved data, the rest points as the observed data to establish a model and use the established model to interpolate the unobserved displacement. We repeated the sampling process until all monitoring points were interpolated. The interpolation results are shown in Figure 2. Lastly, we predicted three days of horizontal displacement for all points using STKF. The statistic results of filtering, interpolation, and prediction are shown in Table 1 followed by the statistic model of R M S 1 :
R M S 1 = i = 1 N ( Y ^ i R i ) 2 / N
where Y ^ i is the estimated value, R i is the true value, and N is the number of Y ^ i .
As Figure 2 shows, the interpolated values coincided with the observed values. As such, we can say that the STKF can well interpolate the point data without monitoring. As listed in Table 1, the interpolation and filter R M S 1 was about 0.03 mm, which had been reduced by about 70% compared with the noise whose standard deviation was 0.1 mm. The interpolation result and filter result were very close, which may be due to the fact that we used a very exact trend matrix. For the prediction result, the prediction accuracy was lower than that of filter and interpolation. Its RMS was around 0.05 mm with the maximum value of 0.091 mm, which was close to the standard deviation of 0.1 mm.

4. Application

4.1. Description of Wuqiangxi Dam Tension Wire Alignment Data

The real horizontal displacement data came from Wuqiangxi Dam tension wire alignment. Wuqiangxi Dam, located in the main stream of Yuan River in the Yuanling County of Hunan Province, was built in 1994. The dam is a concrete gravity dam with a horizontal length of 724.4 m, a crest height of 85.8 m, and a maximum water reservoir volume of 4.29 × 10 9   m 3 . The left side of the dam is a generator set with a power of 1.2 × 10 6   KW and the right side of the dam is a navigation lock. The foundation of this dam can be divided into four types: sandstone, quartzite, slate, and phyllitic slate. This dam is equipped with the automatic monitoring system of wire alignment, inverted plumbline, hydrostatic levelling, seepage monitoring, uplift pressure monitoring, and water level measuring system [25]. The Wuqiangxi Dam contains two tension wire alignments, EX1 and EX2. We selected the EX2 tension wire alignment data in the experiment, which are composed of 23 measuring points and their unidimensional position information, with their original points lying in the far left of the EX2 tension wire alignment, as presented in Table 2. The displacement data contains 2233 daily observations from 1 July 2004 to 11 August 2010. Figure 3 shows the distribution of the tension wire alignment sites. In our experiment, displacement data of EX2_09, EX2_22, and EX2_23 were excluded due to their poor quality. Figure 4a shows the displacement sequences of all selected points. From Figure 4a, we can conclude that the displacement data has significant seasonal trends in the time domain. In order to extract more information from the observation data, Figure 4b shows the amplification of the displacement data of EX2_21 point. There was a gross error that occurred in the time surrounding the 100th day that resulted in some missing data of varying lengths during the whole displacement sequence.

4.2. Filtering, Spatiotemporal Interpolation, and Prediction

In this experiment, displacement data before 8 August 2010 were used for filtering and interpolation. Displacement data from 9–11 August 2010 were used to verify the prediction effectiveness. Similar to the simulation experiment, we used the STKF to filter the observed data firstly. Secondly, in order to check the effectiveness of spatio-temporal interpolation using STKF, we chose 19 points as observation points, and the rest point as a check point. We Repeated the process until all the points were interpolated. Thirdly, we predicted three days’ displacement of all 20 points to check the effectiveness of the spatio-temporal prediction of STKF. Lastly, the whole dam horizontal displacement in spatio-temporal domain was interpolated to view its spatio-temporal variation, and the three days of horizontal displacement prediction for the whole dam was also included. According to the flow diagram of STKF, we randomly chose data from the 700th, 1300th, and 1600th day of dam horizontal displacement in the analysis to select space trend firstly.
Figure 5 shows that the horizontal displacement data had a linear trend, so the chosen trend matrix was F = [ s , 1 ] . Similar to the simulation experiment, a spherical semi-variogram model was also chosen to describe the local spatial variation. As the parameter θ of STKF was fixed, spatio-temporal filtering series, spatio-temporal interpolated series, and spatio-temporal prediction values could be obtained by implementing the STKF recursion part in Section 2.4. Figure 6a plots filtered time-series of all points and the EX2_21 point’s filtering result is enlarged in Figure 6b. Figure 7a shows interpolated time-series of all points and the EX2_21 point’s interpolated result is enlarged in Figure 7b. For the purpose of comparing the filtering result and the interpolated result, EX2_21 point’s filtering result is also displayed in Figure 7b. Table 3 shows statistic results of filtering, interpolation, and prediction. The statistic model of R M S 2 is as follows:
R M S 2 = i = 1 N ( Y ^ i O i ) 2 / N
where Y ^ i is the estimated value, O i is the observed value, and N is the number of Y ^ i .
As Figure 6 shows, the STKF is fairly reliable in filtering the observation points. The result not only reduces the observation error of the series, but also maintains the detailed altitude change. There are big jumps on the 100th day, as shown in Figure 6b, which may have been caused by some instrumental malfunctions on that day. We consider these as gross errors and STKF can reduce their effects significantly. Ellipse in Figure 6b shows the interpolated missing data by spatio-temporal interpolation. The interpolated missing data results can connect with the filter results very well, so we can conclude that the STKF is a suitable method to interpolate the missing data.
In Figure 7a, the interpolated time-series plot coincides with the original time-series plot except for some points with peak values, such as EX2-05, EX2-06, EX2-07, EX2-08, and EX2-21. In Figure 7b, there is a small difference between the results of the 200th day, 800th day, and 1300th day. However, in EX2-11, EX2-12, EX2-13, EX2-14, EX2-15, and so on, the interpolated results perform very well even for the peak values. A possible reason of this fact could be that the chosen trend matrix is not very accurate, and the semi-variogram model cannot describe all the residual spatial variations. From Table 3, we can see that the filter R M S 2 , interpolation R M S 2 , and prediction R M S 2 are about 0.1 mm, 0.5 mm and 0.8 mm, respectively. The largest prediction R M S 2 is up to 1.14 mm in EX2_14, which suggests that the STKF model still has a high prediction accuracy. The larger interpolation R M S 2 occurs in EX2_6 and EX2_7, which is up to 1.12 mm and 1.67 mm, due to the slightly poor performance in peak values.
After verifying the effectiveness of spatio-temporal interpolation and prediction based STKF, the whole dam deformation in spatio-temporal domain can be interpolated and a short-term spatio-temporal prediction for the whole dam can also be predicted. Figure 8 shows interpolated horizontal displacement of the whole dam in the spatio-temporal domain. We equally chose 20 days of interpolated displacement results of the whole dam between the 420th (24 August 2005) day and the 1000th (17 March 2007) day to view the spatio-temporal variation, which is shown in Figure 9. Figure 10 shows the 3 days of prediction results for the whole dam horizontal displacement.
In Figure 8 and Figure 9, the horizontal displacement of the dam shows ‘ascend-descend’ periodic variation in the time domain. Some missing data appears on days 20 June 2006 and 17 November 2006, and all data on day 16 September 2006 is completely missed. There is an interesting phenomenon in that the interpolated results for the missing area are repeatedly similar to the results of the same area in the next subfigure. For example, the missing data interpolation results of 20 June 2006 and 18 September 2006 are very similar to that of 20 July 2006 and 18 September 2006, respectively. This indicates that the interpolated result is not only affected by its surrounding measuring points in the spatial domain, but also by its surrounding data in the temporal domain. In Figure 10, the short-term prediction result agrees with the original data, but its statistic result of the R M S 2 model seems to be larger along with the length of the backward prediction. This suggests that the prediction result may be inaccurate for long-term prediction.

5. Conclusions

This paper applied the STKF model, which combines Kalman filter and Kriging model, to analyze dam deformation data. All points are processed together, and the dam deformation is analyzed as a whole. From the experiments with simulation data and Wuqiangxi Dam tension wire data, we can infer that STKF not only can filter the deformation data noise and predict the deformation of the points effectively both in the temporal domain and spatial domain, but it can also interpolate missing data or date for any position of the dam no matter where the monitoring point is located, which cannot be achieved by the single point Kalman Filter method, because it does not consider the spatial correlation.

Acknowledgments

This work was supported by the State Key Development Program of Basic Research of China (Grant No. 2013CB733303).

Author Contributions

Wujiao Dai conceived and designed the experiments; Ning Liu performed the experiments; Rock Santerre revised the manuscript; Ning Liu and Jiabao Pan wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wu, Z.R. Introduction. In Safety Monitoring Theory & Its Application of Hydraulic Structures; Higher Education Press: Beijing, China, 2003; pp. 1–2. (In Chinese) [Google Scholar]
  2. Ince, C.D.; Sahin, M. Real-time deformation monitoring with GPS and Kalman Filter. Earth Planets Space 2000, 52, 837–840. [Google Scholar] [CrossRef]
  3. Wang, Q.; Sun, H.; Li, W.; Wang, X. Application of Kalman filter analysis method in the deformation monitoring data procession. Chin. J. Eng. Geophys. 2009, 6, 658–661. [Google Scholar]
  4. Li, L.; Kuhlmann, H. Real-time deformation measurements using time series of GPS coordinates processed by Kalman filter with shaping filter. Surv. Rev. 2012, 44, 189–197. [Google Scholar] [CrossRef]
  5. Yang, Y.X.; He, H.B.; Xu, G.C. Adaptively robust filtering for kinematic geodetic positioning. J. Geodesy 2001, 75, 109–116. [Google Scholar] [CrossRef]
  6. Kuhlmann, H. Kalman-filtering with coloured measurement noise for deformation analysis. In Proceedings of 11th FIG Symposium on Deformation Measurements, Santorini, Greece, 25–28 May 2003.
  7. Szostak-Chrzanowski, A.; Chrzanowski, A.; Massiéra, M. Use of deformation monitoring results in solving geomechanical problems—Case studies. Eng. Geol. 2005, 79, 3–12. [Google Scholar] [CrossRef]
  8. Dai, W.J.; Liu, B.; Meng, X.X.; Huang, D.W. Spatio-temporal modelling of dam deformation using independent component analysis. Surv. Rev. 2014, 46, 437–443. [Google Scholar] [CrossRef]
  9. Cressie, N. Comment on ‘An approach to statistical spatial-temporal modeling of meteorological fields’ by MS Handcock and JR Wallis. J. Am. Stat. Assoc. 1994, 89, 379–382. [Google Scholar]
  10. Huang, H.C.; Cressie, N. Spatio-temporal prediction of snow water equivalent using the Kalman filter. Comput. Stat. Data Anal. 1996, 22, 159–175. [Google Scholar] [CrossRef]
  11. Goodall, C.; Mardia, K.V. Challenges in multivariate spatio-temporal modeling. In Proceedings of the XVIIth International Biometric Conference, Hamilton, ON, Canada, 8–12 August 1994.
  12. Mardia, K.V.; Goodall, C.; Redfern, E.J.; Alonso, F.J. The Kriged Kalman filter. Test 1998, 7, 217–282. [Google Scholar] [CrossRef]
  13. Wikle, C.K.; Cressie, N. A dimension-reduced approach to space-time Kalman filtering. Biometrika 1999, 86, 815–829. [Google Scholar] [CrossRef]
  14. Cressie, N.; Wikle, C.K. Space-time Kalman filter. In Encyclopedia of Environmetrics; El-Shaarawi, A.H., Piegorsch, W.W., Eds.; John Wiley & Sons: New York, NY, USA, 2002; Volume 4, pp. 2045–2049. [Google Scholar]
  15. Cortés, J. Distributed Kriged Kalman filter for spatial estimation. IEEE Trans. Auto. Control 2009, 54, 2816–2827. [Google Scholar] [CrossRef]
  16. Sahu, S.K.; Mardia, K.V. A Bayesian kriged Kalman model for short-term forecasting of air pollution levels. J. R. Stat. Soc. C 2005, 54, 223–244. [Google Scholar] [CrossRef]
  17. Lasinio, G.J.; Sahu, S.K.; Mardia, K.V. Modeling rainfall data using a Bayesian Kriged-Kalman model. In Bayesian Statistics and Its Applications; Upadhy, S.K., Singh, U., Dey, D.K., Eds.; Anshan: Tunbridge Wells, UK, 2006. [Google Scholar]
  18. Al-Awadhi, F.A.; Alhajraf, A. Prediction of non-methane hydrocarbons in Kuwait using regression and Bayesian Kriged Kalman model. Env. Ecol. Stat. 2012, 19, 393–412. [Google Scholar] [CrossRef]
  19. Qing, X.Y.; Yang, F.W.; Wang, X.Y. Short-term wind speed forecasting for multiple wind farms using Bayesian Kriged-Kalman mode. Proc. Chin. Soc. Electr. Eng. 2012, 32, 107–114. [Google Scholar]
  20. Cressie, N. Spatial Prediction and Kriging. In Statistics for Spatial Data, 1st ed.; John Wiley & Sons: New York, NY, USA, 1993; pp. 105–210. [Google Scholar]
  21. Sherman, M. Geostatistics. In Spatial Statistics and Spatio-Temporal Data: Covariance Functions and Directional Properties, 1st ed.; John Wiley & Sons: New York, NY, USA, 2011; pp. 21–44. [Google Scholar]
  22. Bookstein, F.L. Principal warps: Thin-plate splines and the decomposition of deformations. IEEE Trans. Patt. Anal. 1989, 11, 567–585. [Google Scholar] [CrossRef]
  23. Shumway, R.H.; Stoffer, D.S. An approach to time series smoothing and forecasting using the EM algorithm. J. Time Ser. Anal. 1982, 3, 253–264. [Google Scholar] [CrossRef]
  24. Oliver, M.A.; Webster, R. A tutorial guide to geostatistics: Computing and modelling variograms and kriging. Catena 2014, 113, 56–69. [Google Scholar] [CrossRef]
  25. Dai, W.J.; Huang, D.W.; Liu, B. A phase space reconstruction based single channel ICA algorithm and its application in dam deformation analysis. Surv. Rev. 2015, 47, 387–396. [Google Scholar] [CrossRef]
Figure 1. Flow diagram of the overall procedure of Space-Time Kalman Filter (STKF). EM algorithm, Expectation Maximization algorithm.
Figure 1. Flow diagram of the overall procedure of Space-Time Kalman Filter (STKF). EM algorithm, Expectation Maximization algorithm.
Ijgi 05 00236 g001
Figure 2. Interpolation of the 23 points.
Figure 2. Interpolation of the 23 points.
Ijgi 05 00236 g002
Figure 3. Location of the tension wire measuring points. The red points in the upper part show EX2 tension wire alignment measurement points; red points in the lower part show EX1 tension wire alignment measurement points; blue points are the inverted plumb measurement points.
Figure 3. Location of the tension wire measuring points. The red points in the upper part show EX2 tension wire alignment measurement points; red points in the lower part show EX1 tension wire alignment measurement points; blue points are the inverted plumb measurement points.
Ijgi 05 00236 g003
Figure 4. Displacement sequences. (a) Displacement sequence of the 20 points; (b) Displacement sequence of EX2_21.
Figure 4. Displacement sequences. (a) Displacement sequence of the 20 points; (b) Displacement sequence of EX2_21.
Ijgi 05 00236 g004
Figure 5. Displacement data on the 700th, 1300th, and 1600th day.
Figure 5. Displacement data on the 700th, 1300th, and 1600th day.
Ijgi 05 00236 g005
Figure 6. Filter and interpolated missing data results. (a) Filter and interpolated missing data results of all points; (b) Filter and interpolated missing data results of EX2_21.
Figure 6. Filter and interpolated missing data results. (a) Filter and interpolated missing data results of all points; (b) Filter and interpolated missing data results of EX2_21.
Ijgi 05 00236 g006aIjgi 05 00236 g006b
Figure 7. Interpolated time series. (a) Interpolation of all points; (b) Filter and interpolation results of EX2_21.
Figure 7. Interpolated time series. (a) Interpolation of all points; (b) Filter and interpolation results of EX2_21.
Ijgi 05 00236 g007
Figure 8. Interpolated horizontal displacement for the whole dam in the spatio-temporal domain.
Figure 8. Interpolated horizontal displacement for the whole dam in the spatio-temporal domain.
Ijgi 05 00236 g008
Figure 9. Equally chosen 20 days of interpolated horizontal displacement for the whole dam in the spatio-temporal domain between 24 August 2005 and 17 March 2007.
Figure 9. Equally chosen 20 days of interpolated horizontal displacement for the whole dam in the spatio-temporal domain between 24 August 2005 and 17 March 2007.
Ijgi 05 00236 g009
Figure 10. Three days of horizontal displacement prediction for the whole dam.
Figure 10. Three days of horizontal displacement prediction for the whole dam.
Ijgi 05 00236 g010
Table 1. R M S 1 of Interpolation, Filter, and Prediction of each point (unit: mm). Interp, Interpolation; Pred, Prediction.
Table 1. R M S 1 of Interpolation, Filter, and Prediction of each point (unit: mm). Interp, Interpolation; Pred, Prediction.
SiteInterpFilterPredSiteInterpFilterPredSiteInterpFilterPred
10.0190.0220.03090.0390.0430.059170.0140.0190.038
20.0240.0250.056100.0370.0370.073180.0090.0170.029
30.0260.0290.051110.0360.0360.077190.0040.0120.003
40.0300.0320.055120.0340.0340.068200.0080.0130.008
50.0330.0350.061130.0310.0320.064210.0150.0170.033
60.0350.0350.081140.0280.0290.055220.0230.0250.040
70.0370.0370.091150.0240.0270.042230.0320.0340.057
80.0370.0370.072160.0190.0220.038
Table 2. Position information of EX2 tension wire alignment measuring points (unit: m).
Table 2. Position information of EX2 tension wire alignment measuring points (unit: m).
SitePositionSitePositionSitePositionSitePosition
EX2_10.5EX2_217.1EX2_341.6EX2_461.1
EX2_581.6EX2_697.1EX2_7115.6EX2_8134.1
EX2_10168.6EX2_11184.1EX2_12205.6EX2_13230.2
EX2_14254.7EX2_15279.2EX2_16286.2EX2_17303.7
EX2_18329.2EX2_19353.7EX2_20378.2EX2_21402.7
Table 3. R M S 2 Filter, Interpolation, and Prediction of each point (unit: mm). Interp, Interpolation; Pred, Prediction.
Table 3. R M S 2 Filter, Interpolation, and Prediction of each point (unit: mm). Interp, Interpolation; Pred, Prediction.
SiteFilterPredInterpSiteFilterPredInterp
EX2_10.050.430.45EX2_120.161.070.52
EX2_20.050.090.15EX2_130.181.070.32
EX2_30.080.430.29EX2_140.181.140.45
EX2_40.120.380.41EX2_150.121.060.2
EX2_50.110.910.84EX2_160.121.050.19
EX2_60.110.311.12EX2_170.120.310.53
EX2_70.110.841.67EX2_180.20.700.74
EX2_80.090.440.76EX2_190.220.580.4
EX2_100.150.660.66EX2_200.230.410.42
EX2_110.120.930.28EX2_210.210.530.92

Share and Cite

MDPI and ACS Style

Dai, W.; Liu, N.; Santerre, R.; Pan, J. Dam Deformation Monitoring Data Analysis Using Space-Time Kalman Filter. ISPRS Int. J. Geo-Inf. 2016, 5, 236. https://doi.org/10.3390/ijgi5120236

AMA Style

Dai W, Liu N, Santerre R, Pan J. Dam Deformation Monitoring Data Analysis Using Space-Time Kalman Filter. ISPRS International Journal of Geo-Information. 2016; 5(12):236. https://doi.org/10.3390/ijgi5120236

Chicago/Turabian Style

Dai, Wujiao, Ning Liu, Rock Santerre, and Jiabao Pan. 2016. "Dam Deformation Monitoring Data Analysis Using Space-Time Kalman Filter" ISPRS International Journal of Geo-Information 5, no. 12: 236. https://doi.org/10.3390/ijgi5120236

APA Style

Dai, W., Liu, N., Santerre, R., & Pan, J. (2016). Dam Deformation Monitoring Data Analysis Using Space-Time Kalman Filter. ISPRS International Journal of Geo-Information, 5(12), 236. https://doi.org/10.3390/ijgi5120236

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