Next Article in Journal
Water, Waste Material, and Energy as Key Dimensions of Sustainable Management of Early Childhood Eco-Schools: An Environmental Literacy Model Based on Teachers Action-Competencies (ELTAC)
Next Article in Special Issue
Technical Note: On the Production and Accuracy of CNC-Manufactured Hydraulic Scale Models
Previous Article in Journal
Bioelectricity Generation and Production of Ornamental Plants in Vertical Partially Saturated Constructed Wetlands
Previous Article in Special Issue
Large-Scale Particle Image Velocimetry for Estimating Vena-Contracta Width for Flow in Contracted Open Channels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predicting Dam Flood Discharge Induced Ground Vibration with Modified Frequency Response Function

1
State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and Hydropower Research, Beijing 100038, China
2
State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072, China
3
School of Water Conservancy and Hydropower, Hebei University of Engineering, Handan 056038, China
4
College of Water Resources Science and Engineering, Taiyuan University of Technology, Taiyuan 030024, China
*
Author to whom correspondence should be addressed.
Water 2021, 13(2), 144; https://doi.org/10.3390/w13020144
Submission received: 2 December 2020 / Revised: 6 January 2021 / Accepted: 7 January 2021 / Published: 10 January 2021
(This article belongs to the Special Issue Measurements and Instrumentation in Hydraulic Engineering)

Abstract

:
Ground vibrations induced by large flood discharge from a dam can damage surrounding buildings and impact the quality of life of local residents. If ground vibrations could be predicted during flood discharge, the ground vibration intensity could be mitigated by controlling or tuning the discharge conditions by, for example, changing the flow rate, changing the opening method of the orifice, and changing the upstream or downstream water level, thereby effectively preventing damage. This study proposes a prediction method with a modified frequency response function (FRF) and applies it to the in situ measured data of Xiangjiaba Dam. A multiple averaged power spectrum FRF (MP-FRF) is derived by analyzing four major factors when the FRF is used: noise, system nonlinearity, spectral leakages, and signal latency. The effects of the two types of vibration source as input are quantified. The impact of noise on the predicted amplitude is corrected based on the characteristics of the measured signal. The proposed method involves four steps: signal denoising, MP-FRF estimation, vibration prediction, and noise correction. The results show that when the vibration source and ground vibrations are broadband signals and two or more bands with relative high energies, the frequency distribution of ground vibration can be predicted with MP-FRF by filtering both the input and output. The amplitude prediction loss caused by filtering can be corrected by adding a constructed white noise signal to the prediction result. Compared with using the signal at multiple vibration sources after superimposed as input, using the main source as input improves the accuracy of the predicted frequency distribution. The proposed method can predict the dominant frequency and the frequency bands with relative high energies of the ground vibration downstream of Xiangjiaba Dam. The predicted amplitude error is 9.26%.

1. Introduction

The construction of high dams with large reservoirs has increased globally owing to the increased utilization of water resources and demand for hydropower stations [1]. The water energy of the flood discharge from a high dam produces strong fluctuating flow loads. These loads can induce strong vibrations in discharge structures, and these vibrations are, in turn, transmitted to surrounding areas through the foundations of the dam. The effect of such vibrations is magnified under special geological conditions such as soft soil conditions [2,3]. Such vibrations can cause ground liquefaction, uneven foundation settlement, cracking of building walls, and harm to the health of residents [4,5,6]. Therefore, ground vibrations in areas near high dams must be investigated, monitored, and controlled. Predicting the vibration induced by flood discharge before it transfers to the surrounding ground, and assessing its impact is an effective way to avoid or mitigate vibration-induced risks. Meanwhile, we hope to predict the vibration of the surrounding ground through the known vibration source, because it is more difficult to obtain the vibration data from the surrounding ground compared with that of the vibration source.
Many studies have predicted the ground vibrations induced by urban traffic [7,8,9]. However, only a few have predicted the ground vibrations induced by flood discharge [10]. In fact, ground vibrations induced by flood discharge can be analyzed as a structural vibration system with input, structure, and output in a manner analogous to the ground vibrations induced by urban traffic. In this case, the vibration system is described as follows: input (discharge excitation)—structure (discharge structure-dam-ground)—output (ground vibration). Therefore, the ground vibrations can be predicted by calculating the frequency response function (FRF) between the input and the output [11,12]. The FRF is a nonparametric estimation method that can reflect the transfer ability of a vibration system [13]; it is generally used to analyze the structural dynamic characteristics, for example, the structural modal parameter estimation [14,15] or the structural damage detection [16,17]. Owing to the characteristics of the FRF, researchers have applied the FRF for predicting the ground vibrations induced by urban traffic [18,19]. First, the FRF can be used to estimate both the intensity and the frequency distribution of ground vibrations [20,21]. In addition to the amplitude, it is more important to accurately predict the vibration frequency because potential resonance with buildings can be avoided by changing the vibration frequency. Second, the FRF can be used for real-time predictions because it remains constant when the input and output are fixed. Then, the unknown output can be estimated by the real-time obtained input and the calculated FRF. However, there are three major differences between the vibrations induced by flood discharge and those induced by urban traffic: (1) Flood discharge excitation is a random excitation with a broadband frequency, and varies with the discharge conditions. I addition, it is more susceptible to environmental factors such as noise and ambient vibration [2]. (2) The vibration transmission problem caused by flood discharge involves more influencing factors such as the coupling between the fluid and the solid. (3) The flood discharge excitation is the superposition of the vibrations generated when the flood flows through different parts of the discharge structures, for example, orifices, guide walls, and stilling pool base slabs [2], and these vibrations are interrelated and interactive in both time and space. Therefore, multiple inputs should be considered for flood discharge vibrations, whereas only a single input needs to be considered for the most urban traffic vibrations. These factors make it difficult to apply the FRF to predict the ground vibrations induced by flood discharge.
In this study, we propose a prediction method with a modified FRF and apply it to ground vibrations induced by flood discharge in the downstream area of Xiangjiaba Dam, a large-scale high dam located on the Jinsha River, Yunnan, China. A multiple averaged power spectrum FRF (MP-FRF) is derived by analyzing four major factors when FRF is used: noise, system nonlinearity, spectral leakages, and signal latency. The effects of two types of vibration source as input are quantified, e.g., the superimposed multiple vibration sources and the main vibration source. The impact of noise on the predicted amplitude is corrected based on the characteristics of the measured signal. The proposed method involves four steps: signal denoising, MP-FRF estimation, vibration prediction, and noise correction.

2. Background

2.1. Xiangjiaba Dam and Ground Vibrations Induced by Flood Discharge

Significant ground vibrations and some structural damage have been reported in Shuifu town near the Xiangjiaba hydropower station, the last cascade hydropower station on the Jinsha River, Yunnan, China. Xiangjiaba Dam is a high concrete gravity dam with a height, length, and crest elevation of 162.00 m, 909.26 m, and 384.00 m, respectively. Its normal storage water level is 380.00 m. The dam has 12 surface orifices and 10 midlevel orifices; these are divided into two symmetrical energy dissipation zones by a middle guide wall. This type of layout of orifices is a typical one in high dams with energy dissipation zones for a hydraulic jump [10].
Figure 1 shows a map indicating the location of the dam and the displacement measurement points, including those in Shuifu town. Owing to a large number of measurement points in the dam area, only those at the orifice, on the foundation of the dam, and on the stilling pool base slab (SPBS) area are marked as T1, T2, and T3, respectively. The blue columns represent the relative RMS of the vertical displacement. Table 1 shows the RMS of the vertical displacement of measurement points. For comparison, the allowable vibration level for strict working area is converted to 0.005 μm in the standard ISO 2631-2: 2003 “Mechanical vibration and shock—Evaluation of human exposure to whole—body vibration—Part 2: Vibration in buildings (1 Hz to 80 Hz)”. The minimum and maximum distances between the dam and the downstream area are 0.5 km and 3.0 km, respectively.
Local residents have reported feeling vibrations in the mountainside area downstream of the dam. Geological studies of this area revealed that its surface and foundation consist of sand layers and mudstones. An ancient watercourse with an average width of ~200 m is located on the right side of the mountain. The maximum thickness of the covering layer over this watercourse is 80 m, whereas that of the covering layer in the east town adjacent to this watercourse is significantly smaller at 20 m. The vibration intensity in the dam area is found to be larger than that in the downstream town area, and the vibration intensity in the ancient watercourse area is generally larger than that in the east town area owing to the magnified effect of its soft soil foundation.

2.2. Sensors and the Signal Collection System

The sensor used is the DP low-frequency vibration displacement sensor. The sensor types are DPS-0.2-5-V-A (vertical direction) and DPS-0.2-5-H-A (horizontal direction). The sampling frequency range is 0.35–200 Hz, the test error is ≤ ±0.1% F.S., the sensitivity is 5 mV/μm.
The vibration response is collected by the vibration displacement sensor installed on the surface of the hydraulic structures and surrounding ground. The sampling time is 120 s, the sampling frequency is 80 Hz. The sensor converts the vibration signal into the voltage signal, and the digital signal of the vibration displacement is obtained after the signal collection system auto processed. Figure 2 shows the sensor installation and the signal collection system.

3. Prediction Methods with Modified FRF

3.1. Principles of FRF

The characteristics of the vibration transmission system are inherent [22]. This means that the system remains constant when the vibration source (input) and the vibration response (output) are selected. Figure 3 shows a ground vibration prediction process using the FRF when input and output don’t contain noise.
F(f) and X(f) are the Fourier spectrum of the input and output signal, respectively. H0(f) is the FRF between F(f) and X(f). Therefore, H0(f) can be estimated using Equation (1) below. H0(f) is dimensionless if and only if both the input and the output are the same physical quantity; otherwise, it has dimensions.
H 0 ( f ) = X ( f ) F ( f )
The new output Fourier spectrum Xi(f) of the same measurement point can be obtained as
X i ( f ) = H 0 ( f ) F i ( f )
Other vibration information such as the time history curves xi(t) can be obtained by using an inverse Fourier transform [23,24,25].

3.2. The Influencing Factors and the Multiple Averaged Power Spectrum FRF

3.2.1. Influencing Factors

As the FRF is estimated using in situ measurement signals, its accuracy will inevitably be affected by interferences such as noise, system nonlinearity, spectral leakages, and signal latency [26,27].
FRF is usually divided into the unbiased estimation or the biased estimation type depending on whether noise is considered or not, respectively. Unbiased estimation should be adopted when both the input and the output contain noise and its effect cannot be ignored, as in the case of ground vibrations induced by flood discharge. Distinguishing between the linear and the nonlinear parts in the input and output signals is very important in calculating the FRF [28]. The coherence functions [29] can reflect the linear degree of the vibration transfer system at different frequencies. Spectral leakage is always caused by the time domain truncation of the signal. This is traditionally solved by windowing in the time domain, such as by using the Hanning window [30,31].
The effect of vibration signal latency on the long distance propagation cannot be ignored. The time difference between the two vibrations is influenced by the topography, the geology, the distance, etc. The correlation coefficient is generally used to describe the correlation degree between the two signals under different phases. Generally, the propagation speed of the vibration in the ground is between 1–6 km/s. In the case of the ground vibration happened near Xiangjiaba Dam, we calculate the correlation coefficient of different phase difference between the vibration source and ground measurement points at different distance. The time interval is 0.125 s, the maximum phase difference is 2.5 s. The variation of correlation coefficient with phase difference is shown in Figure 4, and the correlation coefficients of different measurement points in different phase differences are shown in Table 2. It can be seen that the correlation coefficient of each measurement point is between 0.2–0.6. The changing law of correlation coefficient of each measurement point with phase difference is not the same. The phase difference of the maximum correlation coefficient is between 0 s–1.5 s. Under the influence of topography and geology, the phase difference of the maximum correlation coefficient is different even at the same distance from the vibration source. Therefore, when we calculate the MP-FRF, we take the phase difference of the maximum correlation coefficient of each measurement point as the signal latency.

3.2.2. The Multiple Averaged Power Spectrum FRF

Figure 5 shows a prediction process when the input and output contain noises. The input F(f) contains two components: The Fourier spectrum of the useful signal Z(f) and the Fourier spectrum of background noise M(f). The output X(f) also contains two components: The Fourier spectrum of the useful signal Y(f) and the Fourier spectrum of background noise N(f).
Therefore, the FRF of the system can be expressed as
H ( f ) = Y ( f ) Z ( f ) = X ( f ) N ( f ) F ( f ) M ( f ) = H 0 ( f ) 1 N ( f ) / X ( f ) 1 M ( f ) / F ( f )
Multiple tests were performed on the same vibration system by using the same excitation signal for the input of all tests. For the i-th measurement ( i = 1 ,   2 , , n 1 ) ,
Y i ( f ) N i ( f ) = H ( f ) [ Z i ( f ) M i ( f ) ] ( i = 1 ,   2 , , n 1 )
where Fi(f) and Xi(f) are the i-th measurement of F(f) and X(f), respectively. Zi(f), Mi(f), Yi(f) and Ni(f) are the i-th measurement of Z(f), M(f), Y(f) and N(f), respectively. For the linear transfer system, Fi(f) = F(f) and Xi(f) = X(f) when the input is the same.
Multiplying both ends of the above equation by the conjugation of Z i ( f ) , Z j H ( f ) , (j = i + 1, 2, …, n), and then dividing its expected value by half the sampling time T gives
G Z Y ( f ) G Z N ( f ) = H ( f ) [ G Z Z ( f ) G Z M ( f ) ]
where G Z Y ( f ) = 4 T ( n 1 ) n i = 1 n 1 j = i + 1 n Y i ( f ) Z j H ( f ) is the cross-power spectrum (CS) between Z ( f ) and Y ( f ) , G Z N ( f ) = 4 T ( n 1 ) n i = 1 n 1 j = i + 1 n N i ( f ) Z j H ( f ) is the CS between Z ( f ) and N ( f ) , G Z Z ( f ) = 4 T ( n 1 ) n i = 1 n 1 j = i + 1 n Z i ( f ) Z j H ( f ) is the auto-power (AS) spectrum of Z ( f ) , and G Z M ( f ) = 4 T ( n 1 ) n i = 1 n 1 j = i + 1 n M i ( f ) Z j H ( f ) is the CS between Z ( f ) and M ( f ) .
Owing to the randomness of noise, after a sufficient number of CS and AS transforms, the two noise signals cancel each other. Therefore, after multiple rounds of averaging, G Z N ( f ) = 0 , G Z M ( f ) = 0 , and Equation (5) can be changed to
G Z Y ( f ) = H ( f ) G Z Z ( f )
H ( f ) can be obtained as
H ( f ) = G Z Y ( f ) / G Z Z ( f )
Equation (6) can also be expressed as
H ( f ) = G Z Y ( f ) G Z Z ( f ) = G F X ( f ) + G M X ( f ) + G F N ( f ) + G M N ( f ) G F F ( f ) + G M F ( f ) + G F M ( f ) + G M M ( f )
where G F X ( f ) = 4 T ( n 1 ) n i = 1 n 1 j = i + 1 n X i ( f ) F j H ( f ) is the CS between F ( f ) and X ( f ) , G M X ( f ) = 4 T ( n 1 ) n i = 1 n 1 j = i + 1 n X i ( f ) M j H ( f ) is the CS between M ( f ) and X ( f ) , G F N ( f ) = 4 T ( n 1 ) n i = 1 n 1 j = i + 1 n N i ( f ) F j H ( f ) is the CS between F ( f ) and N ( f ) , G M N ( f ) = 4 T ( n 1 ) n i = 1 n 1 j = i + 1 n N i ( f ) M j H ( f ) is the CS between M ( f ) and N ( f ) , G F F ( f ) = 4 T ( n 1 ) n i = 1 n 1 j = i + 1 n F i ( f ) F j H ( f ) is the AS spectrum of F ( f ) , G M F ( f ) = 4 T ( n 1 ) n i = 1 n 1 j = i + 1 n F i ( f ) M j H ( f ) is the CS between M ( f ) and F ( f ) , G F M ( f ) = 4 T ( n 1 ) n i = 1 n 1 j = i + 1 n M i ( f ) F j H ( f ) is the CS between F ( f ) and M ( f ) , and G M M ( f ) = 4 T ( n 1 ) n i = 1 n 1 j = i + 1 n M i ( f ) M j H ( f ) is the AS of M ( f ) .
Similarly, G M X ( f ) = 0 , G F N ( f ) = 0 , G M N ( f ) = 0 , G M F ( f ) = 0 , G F M ( f ) = 0 , and G M M ( f ) = 0 . Equation (8) can be rewritten as
H ( f ) = G F X ( f ) / G F F ( f )
Equation (9) shows that the influence of noise has been completely eliminated; therefore, H(f) is an unbiased estimation. The estimation accuracy depends not only on the CS between the input and the output and the AS of the input but also on multiple rounds of averaging. This requires several tests with the same input under similar conditions.

3.3. Input Comparison

Figure 6 shows the layout of the Xiangjiaba Dam and the possible vibration sources of downstream. The fluctuating loads at the orifices, SPBS, guide wall, tail-weir, and fall-sill are the main sources affecting ground vibrations. The vibration energy acting on the orifice is much higher than that at other positions under different discharge modes. Further, its contribution to ground vibrations exceeds the sum of the remaining vibration sources [2].

3.3.1. Multiple Vibration Sources as an Input

The input should be a single signal when estimating the FRF. However, the vibration source includes 5 parts. Therefore, we try to use the signals at all vibration sources as an input because compared with the distance between each vibration source and the ground, the distance between the vibration sources is very small. In the present work, the simultaneous signals at 5 parts are superimposed in the time domain as a single input. The displacement of measurement point T9 with the largest vibration amplitude in Figure 1 is selected as the output. Further, 12 similar flood discharge conditions are averaged when estimating the MP-FRF. The similar flood discharge condition means the orifices opening mode and the orifices opening degree are the same, the flow rate and the upstream water level are close (the difference of flow rate is less than 50 m3/s and the difference of upstream water level is less than 0.2 m) during flood discharge. G F X ( f ) is the averaged CS of the displacement of superimposed multiple sources and the displacement at T9, G F F ( f ) is averaged AS of the displacement of the superimposed multiple sources. H ( f ) is the ratio of G F X ( f ) G F F ( f ) . Figure 7 shows the time history curve and CS of G F X ( f ) . Figure 7b indicates that the vibration frequency still shows a flow-load-induced vibration after the vibration source signals are superimposed.
Figure 8 are the MP-FRF and its coherence coefficients. In Figure 8b, all coherence coefficients lie in the range of 0–0.6 in the frequency band of 0–15.0 Hz and that the coherence coefficients in the frequency band of 20.0–40.0 Hz are reduced to 0–0.2. However, a few frequencies have large coherence coefficients; for example, the coherence coefficient at 21.0 Hz is ~0.4, possibly because of the influence of noise. For the MP-FRF (Figure 8a), the frequency band 1.0–6.5 Hz shows a relative high energy and the peak is at ~18 Hz. Large values are also seen at 5.5 Hz, 6.3 Hz, 22.0 Hz, and 29.5 Hz, indicating that ground vibrations are amplified at these frequencies.
Figure 9 shows the time history curves and Fourier spectra of the prediction result and the in-situ measurement of T9.
The amplitude of the prediction result is seen to be significantly larger than that of the in-situ measurement. The RMS of displacement of the prediction and in situ measurement are 0.122 μm and 0.047 μm, respectively. The predicted frequencies have many more peaks than the in-situ measurement one at both low and high frequencies. Specifically, the in-situ measurement has only one frequency band of 0–6.0 Hz with a relative high energy and the peak is at ~3 Hz, whereas the predicted signal has two frequency bands of 0–6.0 Hz and 8.0–12.0 Hz with relative high energies and the peak is at 1.7Hz. This may be because (1) non-source-induced vibrations are amplified when the vibration source signals are superimposed, and (2) white noise still exists at the input and output after multiple averaging. For example, a peak value at ~18 Hz is appeared in the MP-FRF (Figure 8). This causes a relative high vibration energy at ~18 Hz in the predicted Fourier spectra (Figure 9b), but this frequency is not in the in-situ measurement.

3.3.2. Main Vibration Source as an Input

Another input mode is to use the main vibration source along. The vibration at the orifice is the main source of ground vibrations downstream of Xiangjiaba Dam. In this case, G F X ( f ) is the averaged CS of the displacement at the orifice and the displacement at T9, G F F ( f ) is averaged AS of the displacement at the orifice. H ( f ) is the ratio of G F X ( f ) to G F F ( f ) . Figure 10 shows the MP-FRF and its coherence coefficients.
Unlike the results obtained with the superimposed signals as the input, the coherence coefficient in the frequency band of 0–10.0 Hz increases to 0–0.8, indicating that the frequency of the band has less loss during the transfer process. In addition, the coherence coefficient in the frequency band of 10.0–40.0 Hz is generally less than 0.2. In terms of the MP-FRF, the frequency band 1.5–7.5 Hz with a relative high energy and the peak is at ~3 Hz. Large values also occur at 4.1 Hz, 7.2 Hz, and 18 Hz. Overall, there are some other large value peaks (larger than 2 μm2) in the frequency band of 10.0–40.0 Hz. Figure 11 shows the time history curves and Fourier spectra of the prediction result and the in-situ measurement of T9.
The amplitude of the prediction result is also seen to be larger than that of the in-situ measurement. The RMS of displacement of the prediction and in situ measurement are 0.142 μm and 0.064 μm, respectively. We consider that the noise still has an impact on the prediction of vibration. In terms of the vibration frequency distribution, the frequency band 1.0–6.0 Hz with a relative high energy in the in-situ measurement is successfully predicted, and the main frequency of T9 of ~3 Hz is accurately predicted. However, a peak value at 4.1 Hz does not appear in the prediction.
Therefore, compared with the superimposed multiple sources as input, the accuracy of the predicted frequency is improved when the main source is used as the input. Further, when the input or output has two or more frequency bands with relative high energies, the prediction results are more easily affected by factors such as noise that make the prediction amplitude too large. Therefore, the input and output signals need to be preprocessed before using MP-FRF.

3.4. Noise Correction

As the noises in the input and output are filtered at the same time when using MP-FRF, the predicted output is noise-free. Therefore, the effect of noise on the vibration amplitude should be reconsidered. A noise signal based on the signal-to-noise ratio (SNR) formula can be created to simulate real noise and then be added to the prediction signal. The SNR of the signal after noise reduction is
S N R = 20 ln ( σ y 2 σ x y 2 )
where σ y 2 is the variance of the time series after noise reduction, and σ x y 2 is the variance of the noise series after noise reduction. Therefore, the variance of the noise signal can be expressed as
σ x y 2 = σ y 2 e S N R 20

3.5. Prediction Method with Modified FRF

Based on the above analysis, the ground vibration induced by flood discharge from the dam is predicted by using the following four steps:
  • Signal denoising: A modified ensemble empirical mode decomposition (EEMD) and a wavelet threshold filtering method [32] are applied to filter the noise of the main source signals and the ground vibration signals.
  • MP-FRF estimation: The CS of the input and output ( G F X ( f ) ) and the averaged AS of the input ( G F F ( f ) ) are calculated. The same is done for other similar discharge conditions. All CSs and ASs are averaged separately, the MP-FRF is estimated using Equation (9).
  • Vibration prediction: The Fourier spectrum of the ground vibration can be obtained by multiplying the Fourier spectrum of the vibration source and the MP-FRF. The time history curve of the ground vibration can be obtained using an inverse Fourier transform.
  • Noise correction: A noise signal calculated using Equation (11) is added to the predicted time history curve to correct the vibration energy loss due to filtering.
The detailed flow chart for prediction of ground vibration induced by flood discharge is shown in Figure 12.

4. Application to Ground Vibration Downstream of Xiangjiaba Dam

Prediction Results

The displacement at the orifice and T9 are filtered using the modified EEMD and wavelet threshold filtering method. Figure 13 shows their time-history curves before and after noise reduction. Figure 14 shows their Fourier spectra before and after noise reduction.
The denoised input and output of the similar 12 flood discharge conditions mentioned in Section 3.3 are averaged. Figure 15 shows the MP-FRF and its coherence coefficients. The value of the coherence coefficient is seen to be between 0 and 0.7 in the frequency band of 0–10.0 Hz with a peak at 3.5 Hz. The coherence coefficient decreases gradually after 10.0 Hz. Overall, its value is less than 0.2. For the MP-FRF (Figure 15a), the distinct peak is at 3.5 Hz. The main transmission energy is concentrated in the frequency band of 1.0–4.0 Hz. A few vibrations with a relative high energy are at high frequencies like 16.5 Hz, 24.2 Hz, and 24.9 Hz.
Figure 16a shows the RMS of the displacements at T9 under different in situ measurements before and after filtering from small to large. Figure 16b shows the calculated SNR under the corresponding conditions. The SNR is seen to remain at ~15.90 dB, and it does not change significantly with the discharge conditions. According to this characteristic, the noise can be constructed based on Equation (11) and then be added to the predicted result, where the SNR is 15.90 dB.
Figure 17 shows a comparison of the prediction result and the in-situ measurement of the displacement at T9. In terms of the vibration amplitude (Figure 17a), the prediction result after adding noise shows good agreement with the in-situ ones. Both RMS of displacements are close, the prediction result is 0.142 μm and the in-situ measurement is 0.149 μm. In Figure 17b, two dominant frequencies at 1.8 Hz and 3.5 Hz are accurately predicted. The spectrum distribution in the frequency band of 0–10.0 Hz is similar to that of the in-situ measurements.
Then, steps (3)–(4) from Section 3.5 are repeated to predict the RMS of the displacements at T9 under 46 discharge conditions happened during the flood season in the year 2013 (July 2013 to October 2013) of Xiangjiaba Dam. The MP-FRF is calculated with the input (displacement at the orifice) and the output (displacement at T9) under 12 similar flood discharge conditions (Figure 8). The noise signal is calculated to be 15.90 dB using Equation (11). The prediction results are compared with the in-situ measurements shown in Figure 18. In terms of the vibration amplitude, the prediction result shows good agreement with the in-situ ones. The error of the overall amplitude is 9.26%. The same increasing trend, reflecting the effects of higher flow rate, and the similar drastically fluctuating patterns, reflecting the effects of different discharge modes, indicate that the predicted results obtained using the MP-FRF are sensitive to the factors that affect the ground vibration intensity. The comparison between the prediction result and the in-situ measurement at different distances from the dam area is shown in Table 3. They are T4 and T5 at 0.5 km, T9 and T12 at 1 km, T15 and T18 at 1.5 km, and T19 and T21 at 2 km, respectively. It can be seen that all the prediction errors are less than 10%, and the error increases with the increase of distance.

5. Concluding Remarks

This study proposes a prediction method with a modified FRF and applies it to the ground vibrations induced by flood discharge in the downstream area of Xiangjiaba Dam by using in situ measurement data. The MP-FRF is derived by analyzing major factors such as the noise, system nonlinearity, spectral leakages, and signal latency when the FRF is used. The effects of two types of vibration source as input are quantified. The impact of noise on the predicted amplitude is corrected based on the characteristics of the measured signal. The proposed method involves four steps: signal denoising, MP-FRF estimation, vibration prediction, and noise correction. The following results are obtained in this study:
  • As the MP-FRF is used to predict vibrations with a broadband frequency and two or more frequency bands with relative high energies, the prediction results show some frequencies caused by non-vibration sources, and the vibration amplitude is amplified. Therefore, the input and output signals need to be filtered, and the amplitude prediction loss caused by filtering can be corrected by adding a constructed white noise signal to the prediction result.
  • Compared with using the signal at multiple vibration sources after superimposed as the input, using the main source (displacement at the orifice) as the input improves the accuracy of the predicted frequency distribution, and the predicted signals have fewer frequency peaks.
  • The predicted amplitude errors for the downstream area of Xiangjiaba Dam are less than 10%. The predicted results, like the in-situ measurements, are sensitive to factors that affect the ground vibration intensity, such as the flow rate and different discharge modes. The proposed method can predict the dominant frequency and the frequency bands with relative high energies of the downstream ground vibration. The main vibration propagation band is 1.0–10.0 Hz. The MP-FRF remains stable when the vibration source (input) and the vibration response (output) are selected, and the amplitude of the MP-FRF decreases as the distance increases in both the ancient watercourse area and the east town. Therefore, it can be used to predict some other conditions when the inputs are known and the outputs are unknown.
The proposed method can predict ground vibrations in a timely manner when the signal at the vibration source is measured. The ground vibration intensity can be mitigated by controlling or tuning the discharge conditions by, for example, changing the flow rate, changing the opening method of the orifice, and changing the upstream or downstream water level. In addition, for ground vibrations occurring downstream of the Xiangjiaba Dam, the vibration at the orifice is a very strong main vibration source. This may explain why using the main vibration source as the input is better than using multiple vibration sources as the input. In the future, we will conduct further studies in this area and complement the conclusions of the present study.

Author Contributions

Conceptualization, Y.Z. (Yan Zhang) and J.L.; Formal analysis, Y.Z. (Yan Zhang) and G.Z.; Funding acquisition, Y.Z. (Yan Zhang) and Y.L.; Investigation, Y.Z. (Yan Zhang) and Y.L.; Methodology, Y.Z. (Yan Zhang) and S.L.; Project administration, J.L.; Resources, S.L.; Writing—original draft, Y.Z. (Yan Zhang) and Y.Z. (Yanbin Zhao); Writing—review & editing, Y.Z. (Yan Zhang), Y.Z. (Yanbin Zhao) and S.L. All authors have read and agreed to the published version of the manuscript.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article. This work was supported by the National Key R&D Program of China (Grant no. 2018YFC0406900, 2018YFC0406700), the National Natural Science Foundation of China (Grant no. 51709280, 51779277), Applied Basic Research Plan of Shanxi Province (201901D211020).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declared no potential conflict of interest with respect to the research, authorship, and/or publication of this article.

References

  1. Zarfl, C.; Lumsdon, A.E.; Berlekamp, J.; Tydecks, L.; Tockner, K. A global boom in hydropower dam construction. Aquat. Sci. 2015, 77, 161–170. [Google Scholar] [CrossRef]
  2. Lian, J.J.; Zhang, Y.; Liu, F.; Zhao, Q.H. Analysis of the characteristics for ground vibration induced by flood discharge of high dam using cross wavelet transform. J. Renew. Sustain. Energy 2015, 7, 043146. [Google Scholar] [CrossRef]
  3. Kotlyakov, A.V.; Shumakova, E.M.; Artemev, S.A. Dynamics of the coastal zone of the Kuibyshev and Saratov reservoirs in the Tolyatti area and its correlation with the operation regime of the Zhigulevskaya HPP. Water Resour. 2007, 34, 657–662. [Google Scholar] [CrossRef]
  4. Shumakova, E.M.; Kotlyakov, A.V.; Shumakov, G.V. The effect of vibrations in the Zhigulevskii hydropower structure on soils in the nearby territories of Tolyatti city. Water Resour. 2010, 37, 306–310. [Google Scholar] [CrossRef]
  5. Yin, R.G.; Zhang, J.H. Vibrations source and shock absorption scheme research of near-field vibrations caused by flood discharge and energy dissipation of a stilling pool. In Proceedings of the Geo Shanghai International Congress, Tongji University, Shanghai, China, 26–28 May 2014; pp. 107–116. [Google Scholar]
  6. Lian, J.J. The roport of related structures vibration monitoring under bottom outlet discharging condition of Xiluodu Hydropower Station in Jinsha River. In State Key Laboratory of Hydraulic Engineering Simulation and Safety; Tianjin University: Tianjin, China, 2013. (In Chinese) [Google Scholar]
  7. Chen, F.; Takemiya, H.; Huang, M. Prediction and mitigation analyses of ground vibrations induced by high speed train with 3-dimensional finite element method and substructure method. J. Vib. Control 2011, 17, 1703–1720. [Google Scholar] [CrossRef]
  8. Kouroussis, G.; Florentin, J.; Verlinden, O. Ground vibrations induced by InterCity/InterRegion trains: A numerical prediction based on the multibody/finite element modeling approach. J. Vib. Control 2016, 22, 4192–4210. [Google Scholar] [CrossRef]
  9. Lombaert, G.; Degrande, G.; François, S.; Thompson, D.J. Ground-borne vibration due to railway traffic: A review of excitation mechanisms, prediction methods and mitigation measures. In Noise and Vibration Mitigation for Rail Transportation Systems; Springer: Berlin/Heidelberg, Germany, 2015; Volume 126, pp. 253–287. [Google Scholar]
  10. Zhang, J.; Peng, Y.; Xu, W. Hydraulic prediction of near-field vibrations induced by releasing flood. J. Hydraul. Eng. 2017, 143, 05017002. [Google Scholar] [CrossRef]
  11. Tang, X.; Peng, F.; Yan, R.; Gong, Y.H.; Li, X. An effective time domain model for milling stability prediction simultaneously considering multiple modes and cross-frequency response function effect. Int. J. Adv. Manuf. Technol. 2016, 86, 1037–1054. [Google Scholar] [CrossRef]
  12. Wang, Z.W.; Zhu, P. Response prediction for modified mechanical systems based on in-situ frequency response functions: Theoretical and numerical studies. J. Sound Vib. 2017, 400, 417–441. [Google Scholar] [CrossRef]
  13. McKelvey, T.; Guérin, G. Non-parametric frequency response estimation using a local rational model 16th IFAC Symposium on System Identification. In Proceedings of the International Federation of Automatic Control, Brussels, Belgium, 1 July 2012. [Google Scholar]
  14. Jalali, H.; Bonab, B.T. Nonlinearity identification using sensitivity of frequency response functions. J. Vib. Control 2013, 19, 787–800. [Google Scholar] [CrossRef]
  15. Wang, J.T.; Wang, C.J.; Zhao, J.P. Frequency response function-based model updating using Kriging model. Mech. Syst. Signal Process. 2017, 87, 218–228. [Google Scholar] [CrossRef]
  16. Bandara, R.P.; Chan, T.H.T.; Thambiratnam, D.P. Frequency response function based damage identification using principal component analysis and pattern recognition technique. Eng. Struct. 2014, 66, 116–128. [Google Scholar] [CrossRef]
  17. Limongelli, M.P. Frequency response function interpolation for damage detection under changing environment. Mech. Syst. Signal Process. 2010, 24, 2898–2913. [Google Scholar] [CrossRef]
  18. U.S. Federal Transit Administration (FTA). Transit Noise and Vibration Impact Assessment; FTA-VA-90-91003-06; U.S. Federal Transit Administration (FTA): Washington, DC, USA, 2006.
  19. With, C.; Bodare, A. Prediction of train-induced vibrations inside buildings using transfer functions Soil Dynamics & Earthquake. Engineering 2007, 27, 93–98. [Google Scholar]
  20. Lang, Z.Q.; Billings, S.A.; Yue, R.; Li, J. Output frequency response function of nonlinear Volterra systems. Automatica 2007, 43, 805–816. [Google Scholar] [CrossRef]
  21. Maia, N.M.M.; Silva, J.M.M.; Almas, E.A.M.; Sampaio, R.P.C. Damage detection in structures: From mode shape to frequency response function methods. Mech. Syst. Signal Process. 2003, 17, 489–498. [Google Scholar] [CrossRef]
  22. El-Khatib, A.; Guillon, F.; Dômont, A. Vertical vibration transmission through the lumbar spine of the seated subject-First results. J. Sound Vib. 1998, 215, 763–773. [Google Scholar] [CrossRef]
  23. Fritz, O. Fourier Transforms of Distributions and Their Inverses: A Collection of Tables; Academic Press Inc.: New York, NY, USA; A Subsidiary of Harourt Brace Jovanovich Publicatiuon: New York, NY, USA; London, UK, 1973. [Google Scholar]
  24. Germain, B.; Daniel, S.G. Synthesis of gradient-index profiles corresponding to spectral reflectance derived by inverse Fourier transform. Appl. Opt. 1987, 26, 4209–4213. [Google Scholar]
  25. Madsen, C.K.; Zhao, J.H. Optical Filter Design and Analysis; Wiley-Interscience Publication: New York, NY, USA, 1999. [Google Scholar]
  26. Vanhoenacker, K.; Schoukens, J.; Swevers, J.; Vaes, D. Summary and comparing overview of techniques for the detection of non-linear distortions. Die Pharm. 2004, 57, 1241–1256. [Google Scholar]
  27. Brandt, A. Noise and Vibration Analysis: Signal Analysis and Experimental Procedures; A John Wiley and Sons Ltd. Publication: Hoboken, NJ, USA, 2011. [Google Scholar]
  28. Solomou, M.; Rees, D. Frequency domain analysis of nonlinear distortions on linear frequency response function measurements. IEEE Trans. Instrum. Meas. 2005, 54, 1313–1320. [Google Scholar] [CrossRef]
  29. Bendat, J.S. Statistical errors in measurement of coherence functions and input/output quantities. J. Sound Vib. 1978, 59, 405–421. [Google Scholar] [CrossRef]
  30. Schoukens, J.; Rolain, Y.; Pintelon, R. Leakage reduction in frequency response function measurements. In Proceedings of the IEEE Transactions on Instrumentation and Measurement, Ottawa, ON, Canada, 20 November 2006; Volume 55, pp. 2286–2291. [Google Scholar]
  31. Ugryumova, D.; Pintelon, R.; Vandersteen, G. Frequency response function estimation in the presence of missing output data. IEEE Trans. Instrum. Meas. 2015, 64, 541–553. [Google Scholar] [CrossRef]
  32. Zhang, Y.; Lian, J.J.; Liu, F. An improved filtering method based on EEMD and wavelet-threshold for modal parameter identification of hydraulic structure. Mech. Syst. Signal Process. 2016, 68, 316–329. [Google Scholar] [CrossRef]
Figure 1. Map of Xiangjiaba Dam and Shuifu town and layout of the measurement points.
Figure 1. Map of Xiangjiaba Dam and Shuifu town and layout of the measurement points.
Water 13 00144 g001
Figure 2. Sensor installation (a) and the signal collection system (b).
Figure 2. Sensor installation (a) and the signal collection system (b).
Water 13 00144 g002
Figure 3. Prediction process when input and output don’t contain noise.
Figure 3. Prediction process when input and output don’t contain noise.
Water 13 00144 g003
Figure 4. The variation of correlation coefficient with phase difference.
Figure 4. The variation of correlation coefficient with phase difference.
Water 13 00144 g004
Figure 5. Prediction process when input and output contain noise.
Figure 5. Prediction process when input and output contain noise.
Water 13 00144 g005
Figure 6. Location of different flow fluctuating loads.
Figure 6. Location of different flow fluctuating loads.
Water 13 00144 g006
Figure 7. Superimposed signal of multiple vibration sources: (a) Time history curve; (b) CS.
Figure 7. Superimposed signal of multiple vibration sources: (a) Time history curve; (b) CS.
Water 13 00144 g007
Figure 8. MP-FRF (a) and coherence coefficient (b) between the superimposed vibration source and T9.
Figure 8. MP-FRF (a) and coherence coefficient (b) between the superimposed vibration source and T9.
Water 13 00144 g008aWater 13 00144 g008b
Figure 9. The signal of the prediction and the in situ measurement of T9 by using the superimposed vibration source as the input: (a) Time history curves; (b) and Fourier spectra.
Figure 9. The signal of the prediction and the in situ measurement of T9 by using the superimposed vibration source as the input: (a) Time history curves; (b) and Fourier spectra.
Water 13 00144 g009
Figure 10. MP-FRF (a) and coherence coefficient (b) between the main source at orifice and T9.
Figure 10. MP-FRF (a) and coherence coefficient (b) between the main source at orifice and T9.
Water 13 00144 g010
Figure 11. The signal of the prediction and the in situ measurement of T9 by using the main vibration sources as the input: (a) Time history curves; (b) and Fourier spectra.
Figure 11. The signal of the prediction and the in situ measurement of T9 by using the main vibration sources as the input: (a) Time history curves; (b) and Fourier spectra.
Water 13 00144 g011aWater 13 00144 g011b
Figure 12. The flow chart for prediction of ground vibration induced by flood discharge.
Figure 12. The flow chart for prediction of ground vibration induced by flood discharge.
Water 13 00144 g012
Figure 13. Time history curves of original signal and denoised signal at orifice (a) and T9 (b).
Figure 13. Time history curves of original signal and denoised signal at orifice (a) and T9 (b).
Water 13 00144 g013aWater 13 00144 g013b
Figure 14. Fourier spectrum of denoised signal at orifice (a) and T9 (b).
Figure 14. Fourier spectrum of denoised signal at orifice (a) and T9 (b).
Water 13 00144 g014
Figure 15. MP-FRF (a) and coherence coefficient (b) between the signal at orifice and T9.
Figure 15. MP-FRF (a) and coherence coefficient (b) between the signal at orifice and T9.
Water 13 00144 g015
Figure 16. RMS of displacement before and after filtering (a) and its SNR (b) of signal at T9.
Figure 16. RMS of displacement before and after filtering (a) and its SNR (b) of signal at T9.
Water 13 00144 g016aWater 13 00144 g016b
Figure 17. The signal of the prediction and the in situ measurement of signal at T9: (a) Time history curves; (b) CS.
Figure 17. The signal of the prediction and the in situ measurement of signal at T9: (a) Time history curves; (b) CS.
Water 13 00144 g017
Figure 18. Comparison between prediction and in situ measurement RMS of displacement at T9.
Figure 18. Comparison between prediction and in situ measurement RMS of displacement at T9.
Water 13 00144 g018
Table 1. The RMS of the vertical displacement of measurement points.
Table 1. The RMS of the vertical displacement of measurement points.
Measurement PointRMS of the Vertical Displacement
(μm)
Measurement PointRMS of the Vertical Displacement
(μm)
Measurement PointRMS of the Vertical Displacement
(μm)
T11.33T91.24T170.34
T21.48T100.53T180.37
T31.81T110.33T190.40
T40.50T120.38T200.89
T50.40T130.38T210.31
T60.50T140.76T220.29
T70.46T150.85T230.33
T80.36T160.81T240.29
Table 2. The correlation coefficient of different measurement points in different phase differences.
Table 2. The correlation coefficient of different measurement points in different phase differences.
Measurement PointT4T5T9T12T15T18T23T24
Phase Difference (s)
00.5130.1950.5280.4620.5160.3860.4920.268
0.1250.4730.2470.4310.4840.2970.4190.4210.506
0.250.4510.1960.3060.4370.2770.4900.2900.388
0.3750.2070.3680.4110.2450.5010.5130.5070.470
0.50.3740.5730.4310.4210.4220.2990.6220.306
0.6250.3660.5810.3140.5200.4720.2640.4220.315
0.750.4140.3480.5290.3870.4020.3760.3190.379
0.8750.4820.3230.4500.3610.2360.6570.4370.425
10.4170.4310.2210.5120.2860.5500.3480.323
1.1250.4140.4980.4080.2800.5200.2020.2730.390
1.250.3280.2790.4630.4160.6000.2090.2910.460
1.3750.4200.5090.5150.3820.3060.4580.4510.408
1.50.6130.4310.3490.3280.2420.5770.3970.523
1.6250.1680.2060.3760.6300.5100.5240.3980.404
1.750.5230.1960.3790.2470.4600.2660.4020.355
1.8750.4020.6130.4460.4790.4230.2410.4020.441
20.2330.5180.3780.3550.3300.3720.4080.346
2.1250.4960.4980.4330.3970.3870.5720.4490.347
2.250.4160.2380.3430.4560.3540.4250.4810.500
2.3750.4090.4050.3800.4440.3670.2470.3730.368
2.50.3540.4030.5050.1900.4690.3840.1930.435
Table 3. The comparison between the prediction result and the in situ measurement at different distances from the dam area.
Table 3. The comparison between the prediction result and the in situ measurement at different distances from the dam area.
Measurement PointsT4T5T9T12T15T18T19T21
Flow Rate (m3/s) I
(μm)
P
(μm)
I
(μm)
P
(μm)
I
(μm)
P
(μm)
I
(μm)
P
(μm)
I
(μm)
P
(μm)
I
(μm)
P
(μm)
I
(μm)
P
(μm)
I
(μm)
P
(μm)
17110.0100.0110.0080.0090.0260.0270.0080.0080.0180.0190.0080.0080.0050.0070.0070.007
20200.0140.0140.0120.0110.0320.0360.0100.0110.0220.0240.0100.0130.0070.0090.0080.009
22360.0150.0170.0120.0120.0360.0380.0110.0120.0280.0310.0130.0120.0100.0110.0090.009
24170.0100.0130.0080.0130.0250.0250.0100.0080.0170.0170.0070.0070.0080.0110.0060.006
24310.0180.0140.0120.0090.0260.0280.0080.0100.0220.0190.0080.0110.0090.0090.0070.007
24970.0160.0180.0120.0140.0390.0450.0120.0140.0260.0310.0140.0130.0120.0120.0100.011
25600.0130.0120.0140.0170.0260.0290.0080.0090.0180.0200.0080.0090.0080.0090.0070.007
27260.0240.0280.0210.0220.0650.0690.0200.0210.0450.0470.0190.0210.0180.0220.0160.017
28940.0270.0230.0180.0180.0570.0570.0200.0170.0410.0390.0120.0160.0200.0180.0140.016
28960.0170.0210.0130.0160.0410.0420.0130.0130.0280.0280.0120.0120.0130.0130.0100.010
29250.0180.0190.0180.0150.0440.0470.0120.0150.0300.0340.0130.0140.0140.0150.0110.012
30270.0190.0190.0150.0150.0470.0470.0140.0140.0370.0320.0120.0150.0150.0150.0120.012
30410.0170.0150.0100.0120.0310.0360.0100.0110.0260.0250.0090.0110.0090.0120.0080.009
34040.0130.0170.0130.0110.0310.0350.0130.0110.0210.0260.0130.0100.0100.0110.0080.009
35930.0180.0200.0120.0140.0370.0420.0110.0130.0280.0290.0110.0130.0120.0140.0090.011
35970.0150.0160.0120.0130.0370.0400.0140.0120.0250.0270.0110.0120.0150.0170.0090.012
36100.0140.0170.0150.0140.0350.0360.0110.0140.0280.0330.0130.0110.0110.0110.0090.009
38160.0110.0140.0100.0110.0320.0330.0100.0100.0220.0230.0090.0130.0100.0110.0080.008
41820.0160.0190.0130.0110.0400.0350.0120.0110.0320.0240.0120.0100.0170.0140.0100.009
41900.0250.0270.0200.0220.0610.0670.0190.0210.0420.0460.0180.0200.0150.0190.0150.017
43420.0160.0100.0160.0130.0260.0260.0150.0080.0180.0180.0130.0080.0080.0120.0070.009
43550.0130.0130.0100.0100.0310.0310.0100.0130.0210.0250.0090.0120.0130.0100.0080.008
43560.0130.0110.0140.0160.0310.0280.0130.0090.0250.0220.0090.0080.0100.0120.0080.007
43610.0190.0200.0150.0160.0470.0500.0140.0150.0320.0340.0140.0150.0150.0160.0120.012
43630.0140.0110.0090.0090.0280.0270.0090.0080.0190.0180.0120.0080.0090.0090.0070.007
44100.0100.0170.0120.0170.0260.0520.0100.0160.0280.0360.0080.0120.0080.0140.0070.013
45070.0170.0210.0170.0170.0520.0510.0160.0180.0360.0350.0150.0150.0170.0170.0130.013
45240.0190.0230.0150.0180.0510.0570.0160.0170.0350.0390.0120.0170.0140.0180.0130.014
45320.0220.0280.0170.0190.0540.0280.0170.0190.0370.0290.0160.0180.0170.0190.0130.007
46600.0160.0140.0090.0110.0270.0350.0140.0110.0190.0240.0080.0110.0090.0110.0070.009
46900.0150.0120.0120.0100.0360.0300.0110.0090.0250.0260.0110.0090.0120.0100.0090.007
47700.0150.0150.0160.0200.0370.0370.0110.0110.0250.0250.0110.0110.0120.0120.0090.009
47800.0120.0100.0100.0080.0300.0260.0120.0120.0210.0280.0090.0080.0100.0080.0080.006
47920.0140.0140.0110.0130.0340.0350.0100.0110.0230.0240.0120.0150.0130.0110.0090.009
47940.0170.0140.0170.0150.0400.0350.0120.0110.0270.0220.0120.0170.0150.0140.0100.009
47960.0110.0160.0120.0130.0270.0400.0120.0120.0240.0270.0080.0120.0090.0120.0070.010
48000.0160.0170.0130.0130.0390.0410.0120.0130.0270.0280.0120.0120.0130.0130.0100.010
48030.0140.0140.0110.0110.0350.0350.0110.0110.0240.0290.0140.0170.0140.0110.0090.009
48780.0170.0220.0140.0140.0430.0450.0130.0140.0300.0310.0130.0130.0120.0140.0110.011
49600.0180.0160.0160.0130.0380.0400.0120.0120.0260.0270.0110.0130.0120.0130.0090.010
51000.0200.0210.0160.0170.0500.0530.0150.0160.0340.0380.0160.0160.0160.0170.0130.013
53300.0170.0180.0140.0180.0430.0440.0130.0140.0290.0300.0130.0130.0140.0140.0110.011
53880.0140.0150.0130.0120.0400.0380.0120.0120.0290.0260.0160.0180.0130.0120.0100.009
55040.0190.0200.0180.0160.0470.0490.0160.0170.0320.0370.0140.0150.0170.0160.0120.014
59760.0180.0240.0180.0200.0440.0470.0140.0150.0340.0360.0130.0150.0140.0150.0110.012
67220.0230.0210.0150.0170.0460.0520.0160.0180.0380.0430.0170.0190.0170.0170.0120.015
Prediction error (%)7.075.679.263.095.369.439.128.08
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, Y.; Lian, J.; Li, S.; Zhao, Y.; Zhang, G.; Liu, Y. Predicting Dam Flood Discharge Induced Ground Vibration with Modified Frequency Response Function. Water 2021, 13, 144. https://doi.org/10.3390/w13020144

AMA Style

Zhang Y, Lian J, Li S, Zhao Y, Zhang G, Liu Y. Predicting Dam Flood Discharge Induced Ground Vibration with Modified Frequency Response Function. Water. 2021; 13(2):144. https://doi.org/10.3390/w13020144

Chicago/Turabian Style

Zhang, Yan, Jijian Lian, Songhui Li, Yanbing Zhao, Guoxin Zhang, and Yi Liu. 2021. "Predicting Dam Flood Discharge Induced Ground Vibration with Modified Frequency Response Function" Water 13, no. 2: 144. https://doi.org/10.3390/w13020144

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