Deformation Monitoring by GPS at Embankment Dams and Deformation Analysis
Deformation Monitoring by GPS at Embankment Dams and Deformation Analysis
Deformation Monitoring by GPS at Embankment Dams and Deformation Analysis
To cite this article: E. Gökalp & L. Taşçı (2009) Deformation Monitoring by GPS at
Embankment Dams and Deformation Analysis, Survey Review, 41:311, 86-102, DOI:
10.1179/003962608X390021
Article views: 17
Download by: [Library Services City University London] Date: 07 April 2016, At: 23:14
Survey Review, 41, 311 pp.86-102 (January 2009)
ABSTRACT
Deformation measurements are very important for inspecting huge engineering structures in order to
prevent disasters caused by deformation and thus protect national investment. These measurements
require the use of very accurate surveying equipment. In particular, the structure of a dam and the
physical effects it has on its surroundings need continuous monitoring in order to see if unexpected
changes have occurred. If any change is found, remedial work can be performed to protect the
structure. The Global positioning System (GPS) meets all of the requirements stated above and
therefore, GPS receivers were used in this research. The main aims of this work were the examination of
deformation at the crest of the dam caused by water load at different water levels and the weight of the
dam. Additionally, the study was undertaken to determine whether GPS measurements could reach the
Downloaded by [Library Services City University London] at 23:14 07 April 2016
accuracy requirements for dam deformation measurements. In order to monitor and examine the
deformation, a monitoring network consisting of 6 reference points and 11 object points was
established. Measurements were made four times over 2 years using dual frequency GPS receivers with
static methods. The measurements were processed and point coordinates were determined. Then
differences were calculated between periods and the differences were analyzed by iterative weighted
transformation and Fredericton methods to determine the stability of the points. Finally, the results from
the two methods were compared and discussed.
INTRODUCTION
Dams need to be inspected and protected during their complete life cycle. At the
beginning, monitoring occurs are more frequent performed then as the dam ages the
periods between monitoring become longer. The reasons for inspecting the dams are as
follows:
• To check the structural and functional security of the dam
• To protect the lives and property of people living in the vicinity
• To acquire data for the design of similar projects in the future
Dams are affected and deformed because of internal and external loads. These
loads are not constant and can change over time. Deformation and seepage are a
function of these loads. The effects of deformation on Rock-filled dams are altogether
different from that of concrete dams. However, the deformation of both types is largely
characterized as being permanent. The weight of the dam and the hydrostatic pressure
of reservoir back water yield forces which result in a vertical movement on the
structure. The hydrostatic pressure of reservoir can also cause a permanent horizontal
deformation that is perpendicular to centreline of the dam. Elastic behaviour is smaller
in a rock-fill dam [5].
Deformation measurements are one of the important activities in engineering
surveying. There are a variety of surveying instruments and methods with different
capacities and features that are used in deformation monitoring. High precision
measurements are required and have to be made several times. Since GPS fulfils these
requirements it is the most common instrument used in deformation monitoring.
The IWT method is used when there is no previous information about the
movement of points within the network. Calculated displacement values could be
affected from datum selecting or from defining two different datums while adjusting
the measurements taken at two different periods [1]. Therefore, the weight matrix is
obtained iteratively. The IWT method can be used as follows:
• Adjusted coordinates of the points X1, X2 in the deformation network and their
cofactor (covariance) matrices Qx1, Qx2 are calculated with two separate
adjustments.
Displacement values (d1) and the cofactor matrix of d (Qd1) are calculated as
d1 = X2 – X1 (1)
In the beginning, the weight matrix (W) is accepted as W = I. This indicates that all
the points in the network have the same importance. Therefore, the solution is similar
to the Helmert transformation, if some points are given unit weight and the others a
zero weight, i.e., W = diag {I, 0} [2].
Then d is calculated as
d = S (W) d1 (3)
Here, S(W) shows that S matrix, calculated with W=I, can be obtained as
e 0 0 0 z0 − y0 x0
H = 0 e 0 − z0 0 x0 y0
(5)
0 0 e y0 − x0 0 z0
3m*7
87
DEFORMATION MONITORING BY GPS AT EMBANKMENT DAMS AND ANALYSIS
where, di(k) is the ith component of the vector dk after kth iteration
dk+1 is calculated as
dk+1 = S (Wk) dk (7)
If the difference is greater than ε, Wk is calculated again using the value obtained
from equation (7). The iterative procedure continues from equation 5 until the
differences between the successive estimated datum parameters approach the user
defined threshold (ε).
Downloaded by [Library Services City University London] at 23:14 07 April 2016
Note that during the iterative procedure, some dik may be close to 0. In that case,
Wk becomes very large. This could cause numerical instabilities in calculation of Wk.
Therefore, in order to avoid this problem Wk is calculated as
Where df1 and df2 are degrees of freedom used while adjusting first and second
period measurements respectively. σ201 and σ202 are variance factors
If the above equation is accepted, it may be said that the point is stable otherwise it
is unstable. In equation (13), α is significance level, df is the degree of freedom which
is calculated from df = df1 + df2.
88
E. GÖKALP AND L. TAŞÇI
Fredericton Approach
In order to analyze deformation measurements, a generalized approach has been
developed by the Fredericton group. The approach is applicable to any type of
geometrical analysis, both in space and in time domains. Additionally, it can be used
for the detection of unstable points in reference networks and the determination of
strain components and relative rigid body motions in relative networks [3], [2].
The Fredericton approach determines the unstable points within the network by
analyzing the changes (∆l) in length and/or angle between two measurement periods
which are derived from a least squares adjustment of coordinates.
This approach can be applicable as follows:
All possible measurements are derived from adjusted coordinates.
∆l and its covariance matrix C∆l are calculated as
∆l = l2 – l1 (14)
Where, l1 and l2 are derived measurements from first and second period measurements
respectively.
C∆l = C11 + C12 (15)
Downloaded by [Library Services City University London] at 23:14 07 April 2016
Test statistics (∆l/ σ∆l) are calculated using the parameter C∆l. Then, compared with
the values taken from F tables as below.
Significance levels (α) selected are 0.05, 0.025, and 0.01 respectively. In other
words, the confidence levels selected are 95%, 97.5%, and 99% respectively.
The measurements with test statistics greater than their F tables values are listed.
Then, most repeated measurements are determined from the list. A decision is given
with respect to repeating results whether the point is stable or unstable.
APPLICATION
89
DEFORMATION MONITORING BY GPS AT EMBANKMENT DAMS AND ANALYSIS
the measurements, inconsistency among the accuracies of the measurements, and lack
of information about accuracies of some of the measurements.
Location Samsun
Purpose Energy
Capacity 700 MW
90
E. GÖKALP AND L. TAŞÇI
observation period was 45 minutes with a sampling rate of 10 seconds. The satellite
elevation mask was selected at 150 in order to avoid multi-path effect and cycle slip
error.
1006
Reference
Object points
Downloaded by [Library Services City University London] at 23:14 07 April 2016
Two different measurement plans were applied to the survey object points.
Receiver locations and sessions related to the measurements are given in Table3. In the
first plan, two receivers were set over points 1003 and 1004. Then, the third receiver
was set over each object point about 30 minutes (Figure 2a). In the second plan, a
receiver was set over point 0003 during the all observation periods, and the other
points were measured using a leapfrog method (Figure 2b). The main goal of this
measurement plan was to correlate the observations and make loop closures.
The measurements were processed with GeoGenius 2000 software. The accuracy of
the baselines showed a maximum value of 1.2 mm horizontally and 3.2 mm vertically
for the 4 observational periods. The adjusted coordinates and their covariances were
obtained from a free network adjustment. The variance values are given in Table 4. In
order to determine the gross error of the baseline measurements, the Tau test was
applied. The t criteria are extended by a minimum error definition based on a typical
GPS error specification. By default, no errors less than a given value are flagged. The t
criteria are extended by the fixed offset. The Minimum Error Offset is given in the
current distance unit and a relative value (Minimum Error Scale) is given in [ppm].
The default values are 5 mm (offset) + 0.5 ppm (relative). All the baselines passed the
test for selected default values in this research. Additionally, the Minimum Norm
Quadratic Unbiased Estimation (MINQUE) method was used to determine the
accuracy of the baseline measurements instead of taking the values that are given for
the GPS receiver and by software manufacturers as above [6]. The accuracy of
baselines were calculated 4 mm horizontally and 6 mm vertically by MINQUE [8].
91
DEFORMATION MONITORING BY GPS AT EMBANKMENT DAMS AND ANALYSIS
There was little difference between the given and calculated values. Therefore, the
default values given above were used.
N
N
Downloaded by [Library Services City University London] at 23:14 07 April 2016
Fig. 2a. First measurement plan Fig. 2b. Second measurement plan
of the object points of the object points
92
E. GÖKALP AND L. TAŞÇI
the real directions of the displacements, all WGS 84 Cartesian coordinates were
transformed to a local topocentric coordinate system.
Variances (σ2)
Period No
(mm2)
1 0.765
2 0.787
3 0.813
4 0.620
Table 5. Stable and unstable points determined using IWT between the first and second
periods
Displacements (d) di2/qdi σ20 = ci ci<F(1-.05,1,200)
Point No
dN(m) dE(m) ci for dN ci for dE ci<3.90
When the deformation ellipses are drawn from the results of the IWT (App. II., 1,
2, 3) and the tables above are examined the following results can be seen.
Points 0003, 0007, 0011, 0013 and 0019 on the crest, and also points 1002 and
1006 on the reference network, were unstable between the first and second periods.
During these periods, the water level decreased from 170.34 m to 167.53 m.
Points 0011 and 0013 had additional movements. Also, points 1001, 0015 and 0017
showed significant movement between the first and third periods, during which time
the water level decreased from 170.34 m to 164.20 m.
93
DEFORMATION MONITORING BY GPS AT EMBANKMENT DAMS AND ANALYSIS
Table 6. Stable and unstable points determined using IWT between the first and
third periods
Displacements(d) di2/qdi σ20 = ci ci<F(1-.05,1,200)
Point No
ci<3.90
dN(m) dE(m) ci for dN ci for dE
1001 0.0029 0.0097 0.4721 8.0325 UNSTABLE
1002 -0.0103 0.0052 6.1449 3.8158 UNSTABLE
1003 0.0021 0.0000 1.4975 0.0006 STABLE
1004 -0.0017 -0.0007 0.9620 0.3295 STABLE
1005 0.0024 -0.0014 0.8072 0.6271 STABLE
1006 -0.0070 0.0002 4.0799 0.0133 UNSTABLE
0003 0.0064 -0.0050 7.1295 7.8394 UNSTABLE
0005 0.0046 0.0015 2.7285 0.6077 STABLE
0007 0.0058 -0.0029 6.2107 2.7319 UNSTABLE
0009 0.0033 -0.0017 1.0282 0.4601 STABLE
0011 0.0075 0.0003 2.1515 0.0175 STABLE
0013 -0.0018 0.0031 0.3367 1.7885 STABLE
0015 0.0002 -0.0060 0.0150 8.5250 UNSTABLE
0017 -0.0053 0.0022 4.0967 1.3495 UNSTABLE
0019 -0.0076 -0.0010 6.9080 0.3137 UNSTABLE
Downloaded by [Library Services City University London] at 23:14 07 April 2016
Table 7. Stable and unstable points determined using IWT between the first and
fourth periods
Displacements (d) di2/qdi σ20 = ci ci<F(1-.05,1,200)
Point No
dN(m) dE(m) ci for dN ci for dE ci<3.90
Points 1001, 1002 and 1006 had additional movement with points 0009, 0011, 0013
showing significant movement. Despite the fact that some object points seemed not to
be moving, the F test statistics of displacements were very close to F table values (see
Table 7). Indeed, all the points on the crest were moving between the first and fourth
periods. In this time periods, the water level increased from 170.34 m to 177.23 m.
From the results above, the following conclusions can be drawn
• Movement on the dam’s crest was mostly effected by water loading at different
water levels. Water levels were 170.34 m. in first period, 167.53 m. in the
second period, 164.20 m. in the third period, and 177.23 m. in the final period.
94
E. GÖKALP AND L. TAŞÇI
• In this study, the most significant movements were seen at points 0003, 0007,
0013 and 0019 in the middle of the dam’s crest and at the ends of the crest.
• Some of the reference points (1001, 1002, and 1006) were determined unstable
especially between the first and second and the first and third periods. This
shows that the method has the capability of determination of unstable reference
points in reference networks. This will help in selection of the best minimum
constraints and the best deformation models at the later stages of the
deformation analysis. Since the mentioned reference points placed at the left
side of the dam’s reservoir, there is a suspicion that this part of the dam may be
affected by water load at different water levels.
Table 8. Derived baselines from the first and second periods and F test results
∆l/σ∆l ∆l/σ∆l ∆l/σ∆l
Baselines ∆l/σ∆l <F(1,200,0.05)0.5 <F(1,200,0.025)0.5 <F(1,200,0.01)0.5
=1.974 =2.236 =2.645
1002-1004 2.227 1.974 2.236 2.645
1002-1005 2.213 1.974 2.236 2.645
1002-0015 3.030 1.974 2.236 2.645
Downloaded by [Library Services City University London] at 23:14 07 April 2016
Table 9. Point repetition numbers for the first and second periods
Point Numbers
F TEST
1001
1002
1003
1004
1005
1006
0003
0005
0007
0009
0011
0013
0015
0017
0019
0021
0023
F(95%)
0 6 2 8 4 2 4 0 4 4 3 4 2 2 7 0 2
Repetition
F(97.5%)
0 4 0 3 3 1 4 0 3 3 3 4 1 1 6 0 2
Repetition
F(99%)
0 4 0 2 2 0 2 0 3 0 2 2 1 1 4 0 1
Repetition
95
DEFORMATION MONITORING BY GPS AT EMBANKMENT DAMS AND ANALYSIS
Table 10. Derived baselines from the first and third periods and F test results
∆l/σ∆l ∆l/σ∆l ∆l/σ∆l
Baselines ∆l/σ∆l <F(1,200,0.05)0.5 <F(1,200,0.025)0.5 <F(1,200,0.01)0.5
=1.974 =2.236 =2.645
1001-1003 2.3409 1.974 2.236 2.645
1001-1004 2.7442 1.974 2.236 2.645
1001-1005 2.3847 1.974 2.236 2.645
1001-1006 2.4778 1.974 2.236 2.645
1001-0003 3.6948 1.974 2.236 2.645
1001-0005 2.0193 1.974 2.236 2.645
1001-0007 3.2342 1.974 2.236 2.645
1001-0009 2.6720 1.974 2.236 2.645
1001-0015 3.9526 1.974 2.236 2.645
1001-0019 2.5788 1.974 2.236 2.645
1001-0023 3.0530 1.974 2.236 2.645
1002-1003 -2.2799 1.974 2.236 2.645
1003-0017 -2.6682 1.974 2.236 2.645
1003-0019 -2.9423 1.974 2.236 2.645
1004-0003 -3.3358 1.974 2.236 2.645
1004-0007 -2.4693 1.974 2.236 2.645
1004-0015 -2.3278 1.974 2.236 2.645
1005-1006 -2.1053 1.974 2.236 2.645
Downloaded by [Library Services City University London] at 23:14 07 April 2016
In order to determine unstable points in the deformation network, all possible baselines
were derived from coordinates in two separate adjustments of the measurement
periods. Then, the change in baselines, ∆l, and its covariance matrix, C∆l, were
calculated. Test statistics were calculated and compared with the values taken from F
tables. Here, confidence levels selected were 95%, 97.5%, and 99% respectively.
Table 11. Point repetition numbers for the first and third periods
Point Numbers
F TEST
1001
1002
1003
1004
1005
1006
0003
0005
0007
0009
0011
0013
0015
0017
0019
0021
0023
F(95%)
Repetition 11 2 4 4 3 3 7 3 6 3 2 2 3 6 8 2 1
F(97.5%)
Repetition 10 1 4 4 2 1 6 2 5 2 1 2 2 4 8 1 1
F(99%)
Repetition 6 0 2 2 0 0 5 1 3 1 0 1 1 3 4 0 1
96
E. GÖKALP AND L. TAŞÇI
Table 12. Derived baselines from the first and fourth periods and F test results
∆l/σ∆l ∆l/σ∆l ∆l/σ∆l
∆l/σ∆l <F(1,200,0.05)0.5 <F(1,200,0.025)0.5 <F(1,200,0.01)0.5
Baselines
= 1.974 = 2.236 = 2.645
1001-0011 -4.0023 1.974 2.236 2.645
1001-0013 -2.8187 1.974 2.236 2.645
1002-0007 -2.4739 1.974 2.236 2.645
1002-0009 -2.8581 1.974 2.236 2.645
1002-0011 -5.7543 1.974 2.236 2.645
1002-0013 -3.3123 1.974 2.236 2.645
1003-1004 3.1962 1.974 2.236 2.645
1003-0003 4.0081 1.974 2.236 2.645
1003-0005 2.1389 1.974 2.236 2.645
1003-0007 2.2796 1.974 2.236 2.645
1003-0009 2.8807 1.974 2.236 2.645
1004-0003 -2.2123 1.974 2.236 2.645
1004-0011 5.9469 1.974 2.236 2.645
1004-0013 5.9021 1.974 2.236 2.645
1004-0017 3.7958 1.974 2.236 2.645
1005-0009 2.4140 1.974 2.236 2.645
1005-0011 2.1665 1.974 2.236 2.645
1006-0011 4.0202 1.974 2.236 2.645
1006-0013 2.6028 1.974 2.236 2.645
0003-0007 -2.2831 1.974 2.236 2.645
Downloaded by [Library Services City University London] at 23:14 07 April 2016
The measurements with test statistics greater than their F tables values are listed in
Table 8, 10, 12. Finally, the most repeated measurements determined from the list and
the repetition numbers obtained are in Tables 9, 11, 13. From these tables, the points
that have the highest repetition numbers were accepted as unstable. When examining
the tables, there are very close repetition numbers among the points. Therefore, to
avoid misinterpretation the points with repetition numbers which were 50% of the
most repeated point were accepted unstable. For example, the repetition number is 8
for point 1004 with a confidence level of 95% (Table 9). Hence, when the points with
repetition numbers equal to 4 were accepted as unstable otherwise the points were
classified as stable. The summary of the stable and unstable point is listed in Table 14.
Some of the points were accepted as unstable in spite of sudden reduction of their
repetition numbers at confidence levels of 97.5% and 99% respectively. For example,
the repetition numbers of points 1004 and 0009 between first and second periods. This
occurred because the values of ∆l/σl were very close to the F table values. Therefore,
97
DEFORMATION MONITORING BY GPS AT EMBANKMENT DAMS AND ANALYSIS
there is need for prior information about points in order to avoid misinterpretation
concerning the stability and instability of the points.
Table 13. Point repetition numbers for the first and fourth periods
Point Numbers
F TEST
1001
1002
1003
1004
1005
1006
0003
0005
0007
0009
0011
0013
0015
0017
0019
0021
0023
F(95%)
2 4 5 5 2 2 9 4 6 9 10 9 1 5 6 2 4
Repetition
F(97.5%)
2 3 4 4 1 2 8 3 6 6 6 9 0 5 5 1 2
Repetition
F(99%)
2 3 3 4 0 1 5 2 2 5 5 7 0 3 4 0 0
Repetition
Table 14. The summary list of the stable and unstable points
Point
1-2 Periods 1-3 Periods 1-4 Periods
No
1001 STABLE UNSTABLE STABLE
Downloaded by [Library Services City University London] at 23:14 07 April 2016
When examining the displacement vectors, obtained by the IWT method, the points
with a displacement greater than or equal to 4 mm were accepted as unstable. This
coincides with the horizontal accuracy of the observations that were calculated by
MINQUE.
When comparing the Fredericton method with the IWT method, the Fredericton
method has the determination capacity of the same unstable points 86% in 1-2 and 1-4
periods, 65% in 1-3 period with respect to IWT.
The main reason of comparing the IWT and the Fredericton methods is to control
the obtained results. Besides, both methods have the capability of determination of the
single point displacements free from the choice of the datum in the reference networks.
The selection process of the unstable point/points became less conclusive when the
repetition of the failure was similar for most of the points in the Fredericton method.
This requires interpretation of the results.
IWT is a robust method used when there is insufficient knowledge about the
behaviour of the deformable body. This is not actually an advantage. Several methods
have been developed for preliminary identification of the deformation model in the
literature.
98
E. GÖKALP AND L. TAŞÇI
In this work, the hypothesis that the maximum horizontal movements caused by a
water loading effect could occur in the middle of the crest of the dam in arch dams was
proved by the GPS measurements and deformation analysis methods.
A 2-dimensional deformation analysis was made because of the lower accuracy of
the Z component with respect to the X and Y components. Therefore, for the further
work, in order to determine the movement in the vertical dimension, use of precise
levelling or another method that has the potential to provide the vertical accuracy equal
to or greater than the horizontal accuracy of GPS is recommended.
In conclusion, although significant movements have been determined mainly in the
middle parts of the dam’s crest and at the ends of the crest, the Altınkaya dam was
shown to be normal in relation to the water load effect. However, in regard to
magnitude of the displacements there is neither threat to the structural and functional
security of the dam nor to the lives and property of people of people living in the
vicinity.
ACKNOWLEDGEMENT
99
DEFORMATION MONITORING BY GPS AT EMBANKMENT DAMS AND ANALYSIS
100
E. GÖKALP AND L. TAŞÇI
1004
0021
0019 1001
0017
0015 0023 0021
0013
1006 0011
0009 0019
0007
0005 0017
0003
0015 0023
0013
0011
0009
1003 0007
0005
0003
100
1002
1004
0021
Downloaded by [Library Services City University London] at 23:14 07 April 2016
0019 1001
0017
0015 0023 0021
0013
1006 0011
0009
0007 0019
0005
0003
0017
0023
0015
0013
0011
1003 0009
0007
0005
1005
0003
1002
1004
0021 1001
0019
0017 0021
0015 0023
0013
1006 0011 0019
0009
0007
0005 0017
0003
0015 0023
0013
0011
1003
0009
0007
0005
1005 0003
101
DEFORMATION MONITORING BY GPS AT EMBANKMENT DAMS AND ANALYSIS
References
102