Insar Observation and Numerical Modeling of the Earth-dam Displacement of Shuibuya Dam (china)

How to accurately determine the mechanical parameters of rockfill is one of the key issues of concrete-face rockfill dams. Parameter back-analysis using internal or external monitoring data has been proven to be an efficient way to solve this problem. However, traditional internal or external monitoring methods have limitations in efficiency and long-term monitoring. In this paper, the displacement of the Shuibuya concrete-face rockfill dam is monitored by the space-borne Interferometric Synthetic Aperture Radar (InSAR) time series method. Using the InSAR results and the finite element method, the back-analysis of the mechanical parameters of the rockfill dam is investigated, and the back-analysis results of InSAR and levelling are compared. A high correlation of 0.99 for the displacement results generated from InSAR and the levelling offers good agreement between the two methods. The agreement provides confidence that the external InSAR monitoring measurement allows producing a reliable back-analysis and captures the displacement properties of the dam. Based on the identified parameters from the InSAR results, the dam displacement is predicted. The prediction of the maximum settlement of the dam is 2.332 m by the end of 2020, according to the dam displacement characteristics, which agrees well with the results derived from the recorded internal monitoring data. Therefore, the external monitoring results from the InSAR observation can be used as a supplement for traditional monitoring methods to analyse the parameters of the dam.


Introduction
In the stability and security analysis of the concrete-face rockfill dams (CFRDs), the physical parameters of the rockfill are pre-requisite for settlement prediction.Because of the complexity of the dam characteristics, the material properties are influenced by the loading, load size, stress path and other factors in the engineering project.Currently, engineering analogy, expert experience and experimentation are used to determine the material parameters of the rockfill.However, engineering analogies and expert experience require users who have rich experience in engineering and detailed engineering data.Lots of subjective factors are contained in the results.Due to the scale effect, sample disturbance and the impact of random sampling, it is difficult to guarantee the representation of the experiment results [1,2].Moreover, even though the accurate design parameters of the rockfill are obtained by experiment, the values are greatly different from the actual engineering parameters due to the influence of the construction process, construction technology, construction quality, the external

Study Area
Located on the middle stream of the Qingjiang River in the Enshi Autonomous Prefecture, Hubei Province, China, the SBY CFRD (110.3377• N, 30.4374 • E) is approximately 117 km upstream of the city of Enshi and 92 km from the downstream Geheyan Hydropower Plant (Figure 1).The height of the SBY CFRD is up to 233.5 m, and the crest length is 675 m.The normal water level of the reservoir is 400 m above sea level, and the maximum water storage of the reservoir is approximately 4.59 × 10 9 m 3 .The SBY Dam has four generators in its underground power plant, with a total electric generating capacity of 1840 MW.It is a typical concrete-face rockfill dam, and its body is composed of six material zones: arbitrary fill material, bedding material, transition material, primary rockfill, secondary rockfill and downstream rockfill (Figure 2a).The upstream and downstream of the dam has the same slope of 1:1.4.The construction began in October 2002 and can be divided into three reference periods: I Dam construction period (before October 2006).

II
The first reservoir filling period (October 2006-September 2007, when the water level was between 205.06 m above sea level and 389.61 m).III Operation period (after September 2007).
Because of the significance of the project, an improved and detailed settlement monitoring system was implemented for the SBY CFRD, including eleven levelling monitoring stations distributed downstream of the dam.The eleven monitoring stations (WS1 to WS11) were established at elevations of 235 m, 265 m, 300 m, 335 m, and 370 m to monitor the settlement of the downstream dam surface (as shown in Figure 2b).

InSAR Data Analysis and Results
A detailed InSAR data processing steps and analysis about the SBY dam is performed in our previously paper [16].Here, we present a short description of the InSAR processing and related results.Twenty-one ascending ALOS-1 PALSAR L-band (~23.6 cm wavelength) images acquired from the Japan Aerospace Exploration Agency were used to map the displacement of the SBY Dam.The ALOS-1 satellite was launched in January 2006 and terminated in April 2011.Only images from 28 February 2007 to 11 March 2011 were processed, representing the period after completion of the dam (Figure 3).In the single-imaging mode, the resolution of the images is approximately 4.7 m in the slant range and 4.5 m in the azimuth direction.The incidence angle is approximately 38.7°.There

InSAR Data Analysis and Results
A detailed InSAR data processing steps and analysis about the SBY dam is performed in our previously paper [16].Here, we present a short description of the InSAR processing and related results.Twenty-one ascending ALOS-1 PALSAR L-band (~23.6 cm wavelength) images acquired from the Japan Aerospace Exploration Agency were used to map the displacement of the SBY Dam.The ALOS-1 satellite was launched in January 2006 and terminated in April 2011.Only images from 28 February 2007 to 11 March 2011 were processed, representing the period after completion of the dam (Figure 3).In the single-imaging mode, the resolution of the images is approximately 4.7 m in the slant range and 4.5 m in the azimuth direction.The incidence angle is approximately 38.7°.There

InSAR Data Analysis and Results
A detailed InSAR data processing steps and analysis about the SBY dam is performed in our previously paper [16].Here, we present a short description of the InSAR processing and related results.Twenty-one ascending ALOS-1 PALSAR L-band (~23.6 cm wavelength) images acquired from the Japan Aerospace Exploration Agency were used to map the displacement of the SBY Dam.The ALOS-1 satellite was launched in January 2006 and terminated in April 2011.Only images from 28 February 2007 to 11 March 2011 were processed, representing the period after completion of the dam (Figure 3).
In the single-imaging mode, the resolution of the images is approximately 4.7 m in the slant range and 4.5 m in the azimuth direction.The incidence angle is approximately 38.7 • .There are many archived ESA Envisat ASAR C-band (~5.6 cm) images, however the resolution is too low (~20 m).Considering the scale of the dam, ASAR data is not suitable for this study.Although the resolution of the TerraSAR-X and Cosmo-SkyMed X-band (~3 cm) images are much higher (up to 1 m), there is very little archived data in this area at present.That means these data are not suitable for the InSAR time series analysis.are many archived ESA Envisat ASAR C-band (~5.6 cm) images, however the resolution is too low (~20 m).Considering the scale of the dam, ASAR data is not suitable for this study.Although the resolution of the TerraSAR-X and Cosmo-SkyMed X-band (~3 cm) images are much higher (up to 1 m), there is very little archived data in this area at present.That means these data are not suitable for the InSAR time series analysis.The GAMMA processing package was used to generate single look complex (SLC) products from the set of raw SAR images [17].All SLC images were co-registered to a chosen master image (2 June 2008) to match the dataset to its later interferogram generation.Next, an area of approximately 5 km × 5 km (1000 pixels × 1000 pixels) centred on the SBY Dam was cropped from the original 70 km × 70 km SLC images (Figure 4a).All possible interferograms were generated with a spatial baseline below 2000 m and temporal baselines below 30 months.Then, interferograms were investigated visually, and only interferograms with a valid phase and coverage greater than 70% were kept for later InSAR time series analysis.Finally, 26 interconnected interferograms were selected (Figure 3).Topographic fringes of each interferogram was subtracted by using the 1-arc second (~30 m) Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) [18], and the MCF algorithm was used to unwrap the phase [19].Given the image resolution and the scale of the SBY Dam, multi-looking processing was not performed in this study.The GAMMA processing package was used to generate single look complex (SLC) products from the set of raw SAR images [17].All SLC images were co-registered to a chosen master image (2 June 2008) to match the dataset to its later interferogram generation.Next, an area of approximately 5 km × 5 km (1000 pixels × 1000 pixels) centred on the SBY Dam was cropped from the original 70 km × 70 km SLC images (Figure 4a).All possible interferograms were generated with a spatial baseline below 2000 m and temporal baselines below 30 months.Then, interferograms were investigated visually, and only interferograms with a valid phase and coverage greater than 70% were kept for later InSAR time series analysis.Finally, 26 interconnected interferograms were selected (Figure 3).Topographic fringes of each interferogram was subtracted by using the 1-arc second (~30 m) Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) [18], and the MCF algorithm was used to unwrap the phase [19].Given the image resolution and the scale of the SBY Dam, multi-looking processing was not performed in this study.The GAMMA processing package was used to generate single look complex (SLC) products from the set of raw SAR images [17].All SLC images were co-registered to a chosen master image (2 June 2008) to match the dataset to its later interferogram generation.Next, an area of approximately 5 km × 5 km (1000 pixels × 1000 pixels) centred on the SBY Dam was cropped from the original 70 km × 70 km SLC images (Figure 4a).All possible interferograms were generated with a spatial baseline below 2000 m and temporal baselines below 30 months.Then, interferograms were investigated visually, and only interferograms with a valid phase and coverage greater than 70% were kept for later InSAR time series analysis.Finally, 26 interconnected interferograms were selected (Figure 3).Topographic fringes of each interferogram was subtracted by using the 1-arc second (~30 m) Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) [18], and the MCF algorithm was used to unwrap the phase [19].Given the image resolution and the scale of the SBY Dam, multi-looking processing was not performed in this study.Generally, the unwrapped phase consists of five terms: displacement signal, inaccurate orbital information, topographic error, atmospheric propagation delays (atmospheric phase screen, APS), and noise [20].There are various approaches have been developed to reduce the latter four terms based on multi-interferogram analysis [21][22][23][24][25].In this study, the InSAR Time Series with the Atmospheric Estimation Model (InSAR TS + AEM, [26]) is used to perform the time series analysis, which is based on the Small BAseline Subset algorithm (SBAS) (e.g., [22]).The temporal or intermittent coherent pixels strategy was applied to increase the density and distribution of stable pixel [27,28].The Figure 4b shows the average coherence of the study area.
Because the study area is approximately 5 km × 5 km, a linear plane across the entire interferogram is sufficiently accurate to remove the orbital error [20].For a given set of unwrapped interferograms, topographic error signals can be separated because they are proportional to the perpendicular baselines [20].Therefore, the orbital and topographic error can be estimated and individually removed from the unwrapped phase.Taking into account the spatial structure of atmospheric effects (only correlated in space) (e.g., the power-law process [29,30]), the APS can be estimated using a temporary linear velocity (TLV) model and can be distinguished from non-linear surface motion, more details can be found from [26].The algorithm was applied iteratively until convergence was achieved [26,31].A small area (500 m × 500 m) in the village was assumed to be stable during the InSAR observations and was therefore selected as the reference site for the InSAR time series analysis (Figure 4a).For a given pixel, the LOS displacement contains three displacement components in east, north and up direction.It is a function of incidence angle and above three components: where d e , d n , d u , are displacements in east, north and up directions respectively, α is the azimuth of the satellite heading (positive clockwise from the north) and ϑ is the radar incidence angle.
In this study, we assume that the horizontal components are negligible.Assuming this, the vertical component can then be converted from d los using the local incidence angle: d u = −d los /cosϑ.The LOS mean velocity and LOS time series were both converted into vertical direction measurements (Figure 5a).Eleven points at the same locations as the eleven WS stations of the dam were used to investigate the accuracy of the InSAR measurements.These eleven locations were carefully selected by cross-reference of the optical image, SAR amplitude image, and the ground monitoring locations.Finally, a high correlation of 0.9 and an RMS of 1.5 cm/year provided confidence that InSAR can produce reliable monitoring results (Figure 5b).Generally, the unwrapped phase consists of five terms: displacement signal, inaccurate orbital information, topographic error, atmospheric propagation delays (atmospheric phase screen, APS), and noise [20].There are various approaches have been developed to reduce the latter four terms based on multi-interferogram analysis [21][22][23][24][25].In this study, the InSAR Time Series with the Atmospheric Estimation Model (InSAR TS + AEM, [26]) is used to perform the time series analysis, which is based on the Small BAseline Subset algorithm (SBAS) (e.g., [22]).The temporal or intermittent coherent pixels strategy was applied to increase the density and distribution of stable pixel [27,28].The Figure 4b shows the average coherence of the study area.
Because the study area is approximately 5 km × 5 km, a linear plane across the entire interferogram is sufficiently accurate to remove the orbital error [20].For a given set of unwrapped interferograms, topographic error signals can be separated because they are proportional to the perpendicular baselines [20].Therefore, the orbital and topographic error can be estimated and individually removed from the unwrapped phase.Taking into account the spatial structure of atmospheric effects (only correlated in space) (e.g., the power-law process [29,30]), the APS can be estimated using a temporary linear velocity (TLV) model and can be distinguished from non-linear surface motion, more details can be found from [26].The algorithm was applied iteratively until convergence was achieved [26,31].A small area (500 m × 500 m) in the village was assumed to be stable during the InSAR observations and was therefore selected as the reference site for the InSAR time series analysis (Figure 4a).For a given pixel, the LOS displacement contains three displacement components in east, north and up direction.It is a function of incidence angle and above three components: where de, dn, du, are displacements in east, north and up directions respectively, α is the azimuth of the satellite heading (positive clockwise from the north) and ϑ is the radar incidence angle.
In this study, we assume that the horizontal components are negligible.Assuming this, the vertical component can then be converted from dlos using the local incidence angle: du = −dlos/cosϑ.The LOS mean velocity and LOS time series were both converted into vertical direction measurements (Figure 5a).Eleven points at the same locations as the eleven WS stations of the dam were used to investigate the accuracy of the InSAR measurements.These eleven locations were carefully selected by cross-reference of the optical image, SAR amplitude image, and the ground monitoring locations.Finally, a high correlation of 0.9 and an RMS of 1.5 cm/year provided confidence that InSAR can produce reliable monitoring results (Figure 5b).

Mechanical Parameters Back-Analysis Method
Rockfill model parameters back-analysis is used to find one group of parameters that make the calculated and observed displacement correspond.This process includes the following steps (shown in Figure 6): Build the objective function using the InSAR results, as Equation (2).Equation ( 2) presents the InSAR results of ith monitoring point. 2 Perform the strain-stress analysis using the FEM (see Section 4.1.1)and the calculation of the RSM as defined in Section 4.1.2.Using the RSM, simulate the relationship between the parameter set and the displacement of each monitoring point.3 Find the optimal parameter set of the objective function using the modified genetic algorithm (GA) introduced in Section 4.1.3.The optimal parameter set is a combination of the material parameters that minimize the objective function.During the process of searching, these RSM are used to replace the FEM to calculate the fitness of all the measurement points.

Mechanical Parameters Back-Analysis Method
Rockfill model parameters back-analysis is used to find one group of parameters that make the calculated and observed displacement correspond.This process includes the following steps (shown in Figure 6):  The 2-norm of the difference between the calculated and observed displacement at the measuring points is used to build the objective function. ( where is a group of constitutive parameters to be identified, u i is the calculated displacement at the monitoring point i, u i * is the corresponding measured displacement and n is the number of monitoring points used in the back-analysis.

The FEM Model and the Constitutive Model
To further extend the investigation of the settlement determined by the InSAR analysis, we used the finite element method to conduct mechanical parameter back-analysis and strain-stress analysis, which allowed us to account for all geological and geophysical information available for the considered area.As a numerical technique, the FEM is used to find approximate solutions to boundary value problems consisting of partial differential equations and boundary conditions.The FEM subdivides a large problem into smaller, simpler, parts, called finite elements.The simple equations that model these finite elements are then assembled into a larger system of equations that models the entire problem.Then, the FEM uses variational methods from the calculus of variations to approximate a solution by minimizing an associated error function.The principle of the variational methods has been proven to be the equivalent integral of the differential equations and The 2-norm of the difference between the calculated and observed displacement at the measuring points is used to build the objective function. where ) is a group of constitutive parameters to be identified, u i is the calculated displacement at the monitoring point i, u * i is the corresponding measured displacement and n is the number of monitoring points used in the back-analysis.

The FEM Model and the Constitutive Model
To further extend the investigation of the settlement determined by the InSAR analysis, we used the finite element method to conduct mechanical parameter back-analysis and strain-stress analysis, which allowed us to account for all geological and geophysical information available for the considered area.As a numerical technique, the FEM is used to find approximate solutions to boundary value problems consisting of partial differential equations and boundary conditions.The FEM subdivides a large problem into smaller, simpler, parts, called finite elements.The simple equations that model these finite elements are then assembled into a larger system of equations that models the entire problem.Then, the FEM uses variational methods from the calculus of variations to approximate a solution by minimizing an associated error function.The principle of the variational methods has been proven to be the equivalent integral of the differential equations and boundary conditions in mathematics.If the finite element conforms to the convergence criteria, the final approximate solution of the FEM converges to the exact solution of the original mathematical model [32].The Duncan EB model is the most widely used nonlinear elastic constitutive model in the FEM analysis of soil structure [33,34].In particular, we analysed the stress-strain characteristics of the dam in a 3-D FEM model to solve for the retrieved displacements (Figure 7).We defined the subdomain setting of the FEM model using the available geological and structural information.The process and the reservoir water filling were simulated in the analysis (Figure 8).boundary conditions in mathematics.If the finite element conforms to the convergence criteria, the final approximate solution of the FEM converges to the exact solution of the original mathematical model [32].The Duncan EB model is the most widely used nonlinear elastic constitutive model in the FEM analysis of soil structure [33,34].In particular, we analysed the stress-strain characteristics of the dam in a 3-D FEM model to solve for the retrieved displacements (Figure 7).We defined the subdomain setting of the FEM model using the available geological and structural information.The process and the reservoir water filling were simulated in the analysis (Figure 8).boundary conditions in mathematics.If the finite element conforms to the convergence criteria, the final approximate solution of the FEM converges to the exact solution of the original mathematical model [32].The Duncan EB model is the most widely used nonlinear elastic constitutive model in the FEM analysis of soil structure [33,34].In particular, we analysed the stress-strain characteristics of the dam in a 3-D FEM model to solve for the retrieved displacements (Figure 7).We defined the subdomain setting of the FEM model using the available geological and structural information.The process and the reservoir water filling were simulated in the analysis (Figure 8).Due to its simplicity and the explicit physical meaning of its input parameters, the constitutive model proposed by Duncan has been widely used in civil engineering, especially for rockfill dams in China [33,34].It can be described as follows: where E t is the tangent modulus, K is the modulus number, p a is the atmospheric pressure, and σ 3 is the minor principal stress.n and R f are the exponent and failure ratio, respectively.ϕ = ϕ 0 − ∆ϕlg(σ 3 /p a ) is the internal friction angle, and c is the cohesive strength.σ 1 and σ 3 are the major and minor principal stresses, respectively.B t is the bulk modulus.The rockfill of the dam has obvious creep characteristics.In this paper, the Merchant five-parameter creep model is used to describe the creep behaviour of rockfill materials [35].
where ε t is the creep strain, ε f is the limit of the creep strain, α is a coefficient related to the creep strain rate, ε sf is the limit of shear creep strain, and d, b, m c and β are the creep parameters.

The Response Surface Method
Because of the complexity of engineering problems, the implementation of back-analysis is often time-consuming for a large number of the finite element analyses.In this study, the RSM with strong nonlinear mapping ability, is used to simulate the relationship between the parameter sets and the displacement of the monitoring points [4].
Twenty-seven parameter sets are generated using the orthogonal algorithms.Then, the FEM program that takes the static and creep properties of material into account is adopted to calculate the displacement at different monitoring points.The process of construction and reservoir water filling are simulated in the FEM analysis as show in Figure 8.In Li's [4] study, the RSM considers only the static displacement of the rockfill.In this paper, the creep displacement is also taken into consideration.The RSM is given as follows: where x i and y i are the static and creep parameters to be identified, respectively.N and M are the numbers of the static and creep parameters, respectively a, b i , c i , d j and f j are the coefficients of the RSM to be identified.x i and y i can be calculated as follows: where ϕ 01 , ∆ϕ Then, the results of strain-stress analysis are used to calculate the coefficients of the RSM of every monitoring point using Equation (6).The RSM can simulate the relationship between the parameter set and the displacement of every monitoring point.

The Modified Genetic Algorithm
The parameter back-analysis of rockfill dams can be abstracted as an optimization problem in mathematics.The GAs proposed by Hollandin 1975 [36] have proven to be powerful for solving optimization problems.In this paper, a modified GA with global convergence is used to find the optimized parameter set of the objective function.
In GAs, the diversity of the population and the disproportionation between exploitation and exploration is conditioned by the crossover operator, which also influences the speed of convergence and determines the global convergence of the algorithm.The crossover operator is a basic operation with great importance.However, in the crossover operation of the traditional GA, the crossover points and the gene fragments are selected with randomness and blindness.It is difficult to generate new individuals when the selected gene fragments are highly similar to each other, that is, traditional crossover is not very effective in producing new generations that have great differences with their parents.This crossover operation is invalid and may decrease the diversity of the population.
Inspired by genetic engineering and the cloning of superior genes, a novel genetic crossover operator based on the sum of differences in gene fragments (SoDX) is proposed.The fragments' crossing probabilities are first evaluated based on the differences.The gene fragments are then selected and exchanged according to the crossing probability.This process can reduce inbreeding and the possibility of invalid crossover operations.A new GA code named SoDX-NUM was implemented in MATLAB by combining the crossover operator SoDX and the non-uniform mutation (NUM).

The Back-Analysis Results and the Argumentation
In this paper, six monitoring points (WS06-WS11) from levelling and InSAR with high precision and large displacement are used to invert the parameters.Using the displacement back-analysis method proposed in this paper, the identified model parameters set of rockfill based on the levelling and InSAR is determined, as shown in Table 1.Identified 1# and Identified 2# present the back-analysis results based on the levelling and InSAR monitoring results, respectively.Table 1.The identified parameters using the back-analysis method based on the levelling and InSAR results.Using the identified parameter sets of Identified 1# and Identified 2#, the FEM analysis is implemented.In the analysis, the process of dam construction and the water filling are simulated.The results are shown in the Figure 9, including the settlement process of the levelling, the identified parameter set based on the levelling, and the identified parameter set based on InSAR.Using the identified parameter sets of Identified 1# and Identified 2#, the FEM analysis is implemented.In the analysis, the process of dam construction and the water filling are simulated.The results are shown in the Figure 9, including the settlement process of the levelling, the identified parameter set based on the levelling, and the identified parameter set based on InSAR.The traditional back-analysis methods use internal monitoring data.In this study, the external monitoring results are used as the objective function to perform the back-analysis.The correlation of the levelling monitoring data and the computed results based on Identified 1# of the six points are 0.97, 0.96, 0.95, 0.96, 0.93 and 0.96.This suggests that the displacement pattern of the calculated results using the identified parameters set of levelling agree well with the observed levelling values, which indicates that the identified parameters, in general, satisfactorily reflect the displacement properties of the dam.The external monitoring data could be used to map the displacement properties of the dam.The InSAR results are also used to perform the parameters back-analysis.Figure 9 shows the computed settlement based on Identified 2# and the InSAR monitoring results.They have similar magnitude and distribution.
The back-analysis results based on levelling data and InSAR monitoring data are also compared.The displacement from 28 February 2007 to 11 March 2011 is shown in Figure 10a,b and the higher elevation, the greater the displacement increment.In addition, the differences and the correlation are shown in Figure 10c,d.The computed displacement of the two identified parameter sets have little difference with each other (Figure 10c).The correlation of the displacement results generated from InSAR and levelling is 0.99 (Figure 10d).The displacement in the middle of the dam is greater than at the two sides.The greatest total displacement increment of the FEM computed result is 0.467 m, which is similar to the InSAR results.It also shows that the computed results and the results of InSAR have similar displacement increment distributions, and the displacement computed using the parameters Identified 2# has similar magnitude and distribution as the InSAR monitoring results.
Combining the time series analysis and back-analysis, it can be concluded that InSAR technology can used as a supplement for the traditional monitoring method to monitor high dams and can also be used as a method to invert the parameters of the dam.
properties of the dam.The external monitoring data could be used to map the displacement properties of the dam.The InSAR results are also used to perform the parameters back-analysis.Figure 9 shows the computed settlement based on Identified 2# and the InSAR monitoring results.They have similar magnitude and distribution.
The back-analysis results based on levelling data and InSAR monitoring data are also compared.The displacement from 28 February 2007 to 11 March 2011 is shown in Figure 10a,b and the higher elevation, the greater the displacement increment.In addition, the differences and the correlation are shown in Figure 10c,d.The computed displacement of the two identified parameter sets have little difference with each other (Figure 10c).The correlation of the displacement results generated from InSAR and levelling is 0.99 (Figure 10d).The displacement in the middle of the dam is greater than at the two sides.The greatest total displacement increment of the FEM computed result is 0.467 m, which is similar to the InSAR results.It also shows that the computed results and the results of InSAR have similar displacement increment distributions, and the displacement computed using the parameters Identified 2# has similar magnitude and distribution as the InSAR monitoring results.
Combining the time series analysis and back-analysis, it can be concluded that InSAR technology can be used as a supplement for the traditional monitoring method to monitor high dams and can also be used as a method to invert the parameters of the dam.

Settlement Prediction for the SBY Dam
The FEM model, the constitutive model and the back-analysis results can accurately represent the time-dependent displacement of the SBY dam.To study the post-construction settlement of the SBY dam, further FEM analysis is performed to predict the dam displacement using the parameters identified based on the InSAR results.The settlement contours of cross-section 0 + 212 m for 30 December 2020 are given in Figure 11.The maximum settlement is 2.332 m, approximately 1% of the dam height.The maximum settlement occurs in the centre of the dam body, whereas the settlement

Settlement Prediction for the SBY Dam
The FEM model, the constitutive model and the back-analysis results can accurately represent the time-dependent displacement of the SBY dam.To study the post-construction settlement of the SBY dam, further FEM analysis is performed to predict the dam displacement using the parameters identified based on the InSAR results.The settlement contours of cross-section 0 + 212 m for 30 December 2020 are given in Figure 11.The maximum settlement is 2.332 m, approximately 1% of the dam height.The maximum settlement occurs in the centre of the dam body, whereas the settlement of upstream (left) and downstream (right) of the dam is smaller (<40 cm), and the settlement increases towards the centre.The biggest settlement located at the core of the dam mainly because of following two reasons: (1) Because of the gravity of rockfill in different layers and intervals between the different stages of concrete-face construction, most settlement of the lower layer filled rockfill is happened when the upper layer rockfill is completed on the dam.The SBY dam has different material zones, which are shown in Figure 2a).The rockfill is gradually filled according to layer during the construction period.For example, the first stage of concrete-face is constructed after the rockfill of the dam was decreased six or nine months, then the second stage concrete-face and the third concrete face.(2) Because of arching effects caused by different settlement rate of the dam body in the parallel and perpendicular direction of the river.In the direction of perpendicular to the river, due to the V-shaped valley, the arching effect is raised when the settlement rate of the central rockfill is faster than the two sides of the rockfill, which contributes an uplift force to the central rockfill.While in the direction of parallel to the river flow, the other arching effect is caused because of the settlement rate in the central is faster than the upstream and downstream of the dam, which also produces an uplift force to the central rockfill.Therefore, the maximum displacement occurs in the centre of the dam.
The settlement increments from 28 February 2007 to 11 March 2011 and to 30 December 2020 are given in Figure 12a,b, respectively, and their corresponding difference is shown in Figure 12c.Figure 12a,b shows that the displacement increment is approximately 60 cm at the dam crest, and the upstream is greater than the downstream.From the Figure 12c, it can be observed that the settlement at the top and upstream (left) is much greater than that downstream (right) from the dam.This difference is due to water filling; the load of the upstream (left) of the dam is greater than the downstream load (right).The settlement is smaller than 10 cm from 11 March 2011 to 30 December 2020 and is approximately 60 cm from 2007 to 2011.This suggests that the displacement increases rapidly in the early reservoir period, and the settlement trends toward stability in the later operation period.The displacement consists of instantaneous and creep displacement.Most of the later settlement is creep displacement, and the basic characteristic of creep displacement is that it decreases gradually with the time.These results agree with the displacement characteristics of the dam.

Discussions
In our previous study, the InSAR technique was used to investigate the displacement of the SBY Dam, and the reliability was compared with levelling measurements.The results provide confidence that the InSAR time series is a useful tool for CFRD displacement monitoring [16].In this study, to further investigate the inside features of the SBY dam, the FEM, RSM and GA are used to identify the

Discussions
In our previous study, the InSAR technique was used to investigate the displacement of the SBY Dam, and the reliability was compared with levelling measurements.The results provide confidence that the InSAR time series is a useful tool for CFRD displacement monitoring [16].In this study, to further investigate the inside features of the SBY dam, the FEM, RSM and GA are used to identify the mechanical parameters of the rockfill, and then the associated displacement time series and settlement is simulated.The Duncan EB model is the most widely used nonlinear elastic constitutive model in the numerical analysis of soil structure.The triaxial shear tests show that the Duncan EB model is capable of simulating the stress-strain relationship of rockfill materials.The computed results of the previous rockfill dams and the in situ monitoring data have proven the precision of the Duncan EB model [4,9,33].
To examine the stability and security of the CFRD, the physical parameters of the rockfill are pre-requisite for settlement prediction.Due to the size effect, sample disturbance and the impact of random sampling, it is difficult to guarantee the represent activeness of the experiment results.Parameter back-analysis is an efficient way to identify parameters.In this paper, the levelling and InSAR measurement results are used to perform the parameter back-analysis.The results show that the trend, the magnitude and the distribution of the displacement of the SBY Dam downstream surface agree well with each other (Figures 9 and 10).The high correlation of 0.995 and a low RMSE of 0.055 m between the results based on InSAR and the results based on the levelling measurements suggest that the back-analysis based on the InSAR technique is a reliable way to identify the parameters of a rockfill dam.The low RMSE indicates that the back-analysis method is accurate and useful for obtaining the mechanical parameters of the dam.The internal ground monitoring data continues to 16 March 2012.The maximum settlement from the ground monitoring data is 2.536 m, and it is 2.156 m based on the rockfill parameters from the InSAR results, which confirms that the combination of InSAR measurement and the FEM back-analysis is a reliable way to identify the parameters and to predict the displacement of rockfill dams.The results also suggest the stability and security of the SBY CFRD between 2007 and 2020.
The density and coverage of the InSAR measurement can be much higher than conventional survey methods, such as levelling and GPS.In addition to the displacement of the dam body, the stability of the surroundings of the dam can be determined by the InSAR method.This is important because the geologic stability of a dam's surroundings has a great impact on its behaviour.Moreover, many high CFRDs will be built in the northwest of China, where earthquake frequently occur; thus, accurate monitoring is a guarantee of their security.Compared to conventional measurements, the InSAR is an efficient supplement technology to the traditional methods in the long term.

Conclusions
The physical parameters of the SBY dam are identified using the proposed back-analysis method by combining the FEM model and the InSAR time series results.The displacement results using the identified parameters agree well with the observed levelling values in both magnitude and distribution.The settlement prediction of the SBY Dam is performed using the identified parameters, and the maximum settlement is 2.332 m, approximately 1% of the dam height.The settlement of the SBY Dam becomes stable in the later operation period.The back-analysis method using the external monitoring data from InSAR can record the displacement properties of the dam.The InSAR technology can be used as supplement for the traditional monitoring methods to monitor high dams and can also be used to invert the parameters of the dam.

Figure 1 .
Figure 1.The location of the study area and the Shuibuya Reservoir: (a) The location of the Shuibuya Reservoir (red triangle); and (b) View of the Shuibuya concrete-face rockfill dam (CFRD).This figure is adapted from [16].

Figure 2 .
Figure 2. Typical zoning and layout of the exterior settlement monitoring system of the Shuibuya CFRD: (a) The six material zones; and (b) The eleven external settlement monitoring stations.This figure is adapted from [16].

Figure 1 . 17 Figure 1 .
Figure 1.The location of the study area and the Shuibuya Reservoir: (a) The location of the Shuibuya Reservoir (red triangle); and (b) View of the Shuibuya concrete-face rockfill dam (CFRD).This figure is adapted from [16].

Figure 2 .
Figure 2. Typical zoning and layout of the exterior settlement monitoring system of the Shuibuya CFRD: (a) The six material zones; and (b) The eleven external settlement monitoring stations.This figure is adapted from [16].

Figure 2 .
Figure 2. Typical zoning and layout of the exterior settlement monitoring system of the Shuibuya CFRD: (a) The six material zones; and (b) The eleven external settlement monitoring stations.This figure is adapted from [16].

Figure 3 .
Figure 3. Interferogram distribution in terms of the spatial and temporal baseline, adapted from [16].Each black circle represents one SAR image and the black solid line between two circles represents one interferogram.

Figure 4 .
Figure 4.The cropped study area and the average coherence map in radar geometry (not geocoded): (a) The amplitude image of the study area; and (b) The average coherence map, where white indicates pixels with strong coherence and black indicates pixels with low or no coherence.This figure is quoted from [16].

Figure 3 .
Figure 3. Interferogram distribution in terms of the spatial and temporal baseline, adapted from [16].Each black circle represents one SAR image and the black solid line between two circles represents one interferogram.

Figure 3 .
Figure 3. Interferogram distribution in terms of the spatial and temporal baseline, adapted from [16].Each black circle represents one SAR image and the black solid line between two circles represents one interferogram.

Figure 4 .
Figure 4.The cropped study area and the average coherence map in radar geometry (not geocoded): (a) The amplitude image of the study area; and (b) The average coherence map, where white indicates pixels with strong coherence and black indicates pixels with low or no coherence.This figure is quoted from [16].

Figure 4 .
Figure 4.The cropped study area and the average coherence map in radar geometry (not geocoded): (a) The amplitude image of the study area; and (b) The average coherence map, where white indicates pixels with strong coherence and black indicates pixels with low or no coherence.This figure is quoted from [16].

Figure 5 .
Figure 5. InSAR Vertical displacement velocity of the study area in radar geometry (not geocoded), and comparisons between the levelling: (a) The vertical displacement velocity, the blue star presents the reference location; and (b) The comparisons between InSAR and levelling measurements, black lines is a 1:1 line.This figure is adapted from [16].

Figure 5 .
Figure 5. InSAR Vertical displacement velocity of the study area in radar geometry (not geocoded), and comparisons between the levelling: (a) The vertical displacement velocity, the blue star presents the reference location; and (b) The comparisons between InSAR and levelling measurements, black lines is a 1:1 line.This figure is adapted from [16].

1 .
Build the objective function using the InSAR results, as Equation(2).Equation (2) presents the InSAR results of ith monitoring point.2. Perform the strain-stress analysis using the FEM (see Section 4.1.1)and the calculation of the RSM as defined in Section 4.1.2.Using the RSM, simulate the relationship between the parameter set and the displacement of each monitoring point.3. Find the optimal parameter set of the objective function using the modified genetic algorithm (GA) introduced in Section 4.1.3.The optimal parameter set is a combination of the material parameters that minimize the objective function.During the process of searching, these RSM are used to replace the FEM to calculate the fitness of all the measurement points.

Figure 6 .
Figure 6.The flow chart of the parameter back-analysis.

Figure 6 .
Figure 6.The flow chart of the parameter back-analysis.

Figure 8 .
Figure 8.The process of construction and reservoir water filling.The numbers indicate the step of the dam construction and reservoir filling.

Figure 8 .
Figure 8.The process of construction and reservoir water filling.The numbers indicate the step of the dam construction and reservoir filling.

Figure 8 .
Figure 8.The process of construction and reservoir water filling.The numbers indicate the step of the dam construction and reservoir filling.

1 , k 1 , n 1 ,
R f1 and K b1 are the values of the Duncan EB model parameters obtained from the laboratory; ϕ 0 , ∆ϕ, k, n, R f and K b are the values after zooming in; a 1 , b 1 , m c1 , β 1 and d 1 are the values of the Merchant creep model parameters obtained from the laboratory; and a, b, m c , β and d are the values after zooming in.

K 1 ,
n 1 , K b1 and m 1 are the parameters of the Duncan EB model in the first material zone of the rockfill dam.K 2 , n 2 , K b2 and m 2 are the parameters of the Duncan EB model in the second material zone of the rockfill dam.a, b, m c , β and d are the parameters of the creep constitutive model of the rockfill.

Identified 1 #
and Identified 2# are the back-analysis results of the levelling and InSAR, respectively.

Figure 9 .
Figure 9.The displacement for WS05-WS11 based on the levelling, levelling-identified, InSAR and InSAR-identified parameters.The traditional back-analysis methods use internal monitoring data.In this study, the external monitoring results are used as the objective function to perform the back-analysis.The correlation of the levelling monitoring data and the computed results based on Identified 1# of the six points are

Figure 10 .
Figure 10.Comparison of the displacement increment from 28 February 2007 to 11 March 2011: (a) Vertical displacement computed by the FEM model and the identified parameters based on levelling (1#); (b) Vertical displacement computed by FEM model and the identified parameters based on the InSAR results (2#); (c) The displacement difference between (a) and (b); (d) The Correlation between InSAR and levelling measurements, and black line in is a 1:1 line.XCoor means the direction of the river flow, YCoor means the direction from the left bank to the right bank.The units of X and Y-axes are in meters.

Figure 10 .
Figure 10.Comparison of the displacement increment from 28 February 2007 to 11 March 2011: (a) Vertical displacement computed by the FEM model and the identified parameters based on levelling (1#); (b) Vertical displacement computed by FEM model and the identified parameters based on the InSAR results (2#); (c) The displacement difference between (a) and (b); (d) The Correlation between InSAR and levelling measurements, and black line in is a 1:1 line.XCoor means the direction of the river flow, YCoor means the direction from the left bank to the right bank.The units of X and Y-axes are in meters.
Remote Sens. 2016, 8, 877 13 of 17decreases gradually with the time.These results agree with the displacement characteristics of the dam.

Figure 11 .
Figure 11.The settlement contours of cross-section 0 + 212 m for 30 December 2020.XCoor means the direction of the river flow.The units of X and Y-axes are in meters.

Figure 11 .
Figure 11.The settlement contours of cross-section 0 + 212 m for 30 December 2020.XCoor means the direction of the river flow.The units of X and Y-axes are in meters.

Figure 11 .
Figure 11.The settlement contours of cross-section 0 + 212 m for 30 December 2020.XCoor means the direction of the river flow.The units of X and Y-axes are in meters.

Figure 12 .
Figure 12.The settlement increment contours of cross-section 0 + 212 m: (a) The settlement increment from 28 February 2007 to 11 March 2011; (b) The settlement increment from 28 February 2007 to 30 December 2020; and (c) The settlement increment from 11 March 2011 to 30 December 2020.XCoor means the direction of the river flow.The units of X and Y-axes are in meters.

Figure 12 .
Figure 12.The settlement increment contours of cross-section 0 + 212 m: (a) The settlement increment from 28 February 2007 to 11 March 2011; (b) The settlement increment from 28 February 2007 to 30 December 2020; and (c) The settlement increment from 11 March 2011 to 30 December 2020.XCoor means the direction of the river flow.The units of X and Y-axes are in meters.

Table 1 .
The identified parameters using the back-analysis method based on the levelling and InSAR results.