Next Article in Journal
Comparison of DNA Extraction Methods for the Direct Quantification of Bacteria from Water Using Quantitative Real-Time PCR
Next Article in Special Issue
Prediction Model of Hydropower Generation and Its Economic Benefits Based on EEMD-ADAM-GRU Fusion Model
Previous Article in Journal
Assessment of Cyprinus carpio Scales as a Low-Cost and Effective Biosorbent for the Removal of Heavy Metals from the Acidic Mine Drainage Generated at Rosia Montana Gold Mine (Romania)
Previous Article in Special Issue
Separation Method of Main and Foreign Water for the Measuring Weirs of Danjiangkou Earth-Rock Dam
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Working Mode Identification Method for High Arch Dam Discharge Structure Based on Improved Wavelet Threshold–EMD and RDT Algorithm

1
School of Infrastructure Engineering, Nanchang University, Nanchang 330031, China
2
College of Water Conservancy and Ecological Engineering, Nanchang Institute of Technology, Nanchang 330099, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(22), 3735; https://doi.org/10.3390/w14223735
Submission received: 13 September 2022 / Revised: 30 October 2022 / Accepted: 15 November 2022 / Published: 17 November 2022
(This article belongs to the Special Issue Safety Monitoring and Management of Reservoir and Dams)

Abstract

:
Prototype vibration response data of high arch dam discharge structures inevitably mix various noises under the discharge excitation, which adversely affects the accuracy of the working modal identification of the structure. To effectively filter noise and reduce modal aliasing for better identification accuracy, this study proposes an improved modal threshold identification method based on an improved wavelet threshold–empirical mode decomposition (EMD) and random decrement technique (RDT) algorithm for high arch dam discharge structures. On the basis of the measured vibration response data of the dam, the wavelet threshold is adopted to filter out most of the high-frequency white noise and to reduce the boundary accumulation effect of EMD decomposition. Detrended fluctuation analysis (DFA) is utilized to filter white noise and low-frequency flow noise after EMD decomposition. The natural frequency and damping ratio of the structure system are obtained by the improved RDT algorithm. The engineering examples show that the proposed method can accurately filter the measured vibration response signal noise of the discharge structure, retain the signal characteristic information, improve the accuracy of working modal recognition of the structural vibration response, avoid the complex ordering process of the system, and ease the working modal parameter identification of high arch dam discharge structures. This method can be applied to the mode identification of other large structures, as well.

1. Introduction

High arch dams are located in deep mountain valleys, which have a large flow per unit width and high flood peak energy during discharging. Flood-induced vibrations of arch dams are a key technical problem related to the safety of arch dam structures. Modal parameter identification of discharge structures is an important technical approach to confirm the operating state of dam structures [1,2,3,4,5,6,7,8]. Various monitoring facilities have been applied to monitor the environmental variables and vibration responses of dams. Noises are inevitably contained in the prototype vibration response data of discharge structures due to the complex working environment of data collection and the difficulty in excitation of large discharge structures. The condition considerably affects the accuracy of working modal identification of discharge structures. Therefore, it is necessary to eliminate the environmental background noise and to accurately identify the modal parameters of the high arch dam discharge structure under the working environment load excitation.
Based on the identification of structural modal parameters under environmental excitation, traditional modal identification methods, including stochastic subspace identification (SSI) [9,10], RDT [11], Hilbert–Huang Transform (HHT) [12,13], etc., have been successfully applied in the identification of structural modal parameters under environmental excitation. The frequency domain decomposition method was applied by Li [14,15] to identify the modal parameters of a high arch dam under discharge excitation, which verified the feasibility of identifying the working modal parameters of the structure by using only the fluid-induced vibration response of the structure. Tian et al. [16,17] applied RDT to the working modal identification of a cantilever beam. However, RDT seems to be powerless for dealing with non-stationary signals. Zhang et al. [18] applied the HHT method to identify the modal parameters of the dam discharge structure. The ability of HHT to identify the frequency-intensive structure is confirmed, but there are end effects, modal aliasing, and other defects in the actual calculation process. Therefore, the method of combining RDT with HHT to suppress the deficiencies of the two methods has been proposed by many experts. He et al. [19,20] combined HHT and RDT to identify the working mode of the Nanjing Yangtze River Bridge. Han et al. [21] applied HHT–RDT method to the working mode identification of structures under seismic action. However, the original signal obtained by the prototype vibration response test under the discharge excitation had a certain deviation due to the existence of background noise. Accordingly, an accurate reflection of the mechanical characteristics of the structure could not be obtained. Noise reduction of the prototype vibration response data is essential to obtain accurate structural feature information. For the noise reduction of signals, many experts have conducted a series of discussions. Tang et al. [22,23,24] applied wavelet analysis for signal denoising, but the filtering effect mainly depends on the selection of the wavelet threshold. Zhao et al. [25,26] applied the EMD method for signal decomposition based on the time-scale characteristics of the data. It is easy for the mixed noise and white noise signals to cause a loss of effective information. Li et al. [27,28] applied the wavelet threshold denoising method based on EMD. However, EMD decomposition leads to different degrees of end effects and modal aliasing when the signal contains sound noise or white noise [29]. Zhang et al. [30,31] applied the hybrid method of wavelet analysis and EMD noise reduction, which can effectively extract the value information of low-SNR discharge structures, but the wavelet basis is difficult to determine.
On the basis of the above-mentioned studies on theoretical analysis and simulation experiments of high arch dam discharge structure, this study proposes an improved wavelet threshold–EMD and RDT algorithm for working modal identification of high arch dams. Prior to modal identification, the improved wavelet threshold–EMD algorithm is used to effectively process the signal. The improved wavelet threshold function is used to filter out most of the high-frequency white noise and reduce the boundary accumulation effect of EMD decomposition. The white noise and low-frequency flow noise are further filtered by DFA after EMD decomposition. The improved RDT algorithm is applied to the working mode identification of the high arch dam discharge structure. The EMD process of the vibration response signal is controlled by a band-pass filter to obtain each modal component of the structure, the free attenuation signal of each modal component is extracted by RDT method, and the natural frequency and damping ratio of the high arch dam discharge structure are obtained by HHT and least squares fitting. The method is applied to the measured vibration response analysis of a high arch dam structure.

2. Basic Theory

2.1. Improved Wavelet Threshold–EMD Hybrid Algorithm for High Arch Dam Noise Reduction

The data collection working environment of the high arch dam discharge structure is complicated, and the vibration interference type is shown in Figure 1. Some low-frequency interference signals and white noise are mixed in the prototype vibration response data of the discharge structure. An improved wavelet threshold–EMD hybrid algorithm is proposed to filter noise in the discharge structure response of arch dams in order to retain a large number of useful signals in the original signal.
This method fully combines the advantages of wavelet threshold and EMD denoising: wavelet threshold is used to separate the high-frequency noise in the signal, to suppress the endpoint effect of EMD, to reduce the influence of mode mixing, and to lay a foundation for the subsequent EMD noise reduction; EMD decomposition further separates the white noise and low-frequency flow noise and improves noise reduction accuracy. The basic principle is as follows:
The key to wavelet threshold denoising [32,33] is not only the selection of the threshold function, but also the estimation of the wavelet threshold size. The main idea of signal denoising is to reconstruct the compressed coefficient on the basis of the preset wavelet transform coefficient of the compressed signal.
In selecting the wavelet threshold, the threshold selection method of traditional wavelet threshold denoising is:
T = σ n 2 ln N
where σn is the standard deviation and N is the signal length.
For the vibration signal of the discharge structure submerged by noise, the wavelet coefficient of the noise decreases with the increase of the number of decomposition layers, which affects the noise reduction effect. The traditional threshold formula calculates the global threshold, which is obviously not suitable for the noise reduction processing in this study. An improved wavelet threshold denoising method is proposed. The number of decomposition layers j is determined by a white noise test [34]. The threshold functions of each layer are designed according to the number of decomposition layers j. The noise components in the high-frequency coefficients are separated to completely eliminate the noise components while retaining the original useful signal. The improvement process is expressed as:
T = t σ e σ N j 2 lg N j / lg ( j + 1 )
where Nj is the signal data length; j is the number of decomposition layers; and σ is the standard variance of noise, whose calculation formula can be written as:
σ = m e d i a n ( a b s ( D ) ) / 0.5 N j )
where D is the absolute median of the high-frequency coefficient.
Traditional threshold functions include the soft threshold function and the hard threshold function. The hard threshold function is used to keep the wavelet coefficient constant when the absolute value of the wavelet coefficient is greater than or equal to the given threshold value, otherwise it should be set to zero. The formula is as follows:
w ^ j , k = w j , k ,   w j , k T 0 , w j , k < T
where w ^ j , k is the estimated wavelet coefficient after the threshold is applied, and w j , k is the wavelet decomposition coefficient of the signal.
The soft threshold function is used to subtract the threshold when the absolute value of the wavelet coefficient is greater than or equal to a given threshold, otherwise it should be set to zero. The formula is as follows:
w ^ j , k = sgn ( w j , k ) ( w j , k T ) w j , k T 0 , w j , k < T
The hard threshold function can effectively preserve local features such as signal edges. The soft threshold algorithm is relatively smooth but will cause distortion phenomena such as edge blur. To overcome the shortcomings of the soft and hard threshold algorithms, this study proposes an improved wavelet threshold function, which is defined as follows:
w ^ j , k = sgn ( w j , k ) [ w j , k t T e ( 1 t ) ( T w j , k ) ] w j , k T 0 w j , k < T
where t is the adjustment factor.
When w j , k T , the following function can be obtained as:
f ( x ) = sgn ( x ) [ x t T e ( 1 t ) ( T x ) ]
Thus,
f ( x ) x = [ x t T e ( 1 t ) ( T x ) ] x = 1 t T e ( 1 t ) ( T x ) x 1 ( x + ) x > 0 f ( x ) x = [ x + t T e ( 1 t ) ( T + x ) ] x = 1 + t T e ( 1 t ) ( T + x ) x 1 ( x ) x < 0
The improved threshold function is equivalent to a function between the soft and hard threshold functions. When the original signal contains a large number of sharp mutations, T is moved toward the hard threshold direction. On the contrary, when the original signal is relatively smooth, T is moved toward the soft threshold. Thus, the size of T can be adjusted adaptively to obtain an excellent denoising effect.
The improved wavelet threshold noise reduction filters out high-frequency white noise in the signal while reducing EMD endpoint effects and modal aliasing. The EMD noise reduction is further performed on the vibration signal, and the white noise and the low-frequency water flow noise are further filtered.
The EMD method [35,36] iteratively decomposes several IMF components and one residual signal by using the screening process based on the cubic spline interpolation method for the signal x(t), which can be expressed as:
x ( t ) = j = 1 n c j ( t ) + r n ( t )
where n is the number of decomposed IMF components and rn(t) is the residual signal of x(t).
The iterative screening of IMF components [37] is an important part of EMD decomposition. According to the orthogonality of EMD decomposition, if the decomposition is correct, the IMF components are orthogonal to each other, and the IMF orthogonal index is defined as:
I O = t = 0 T i = 1 n + 1 j = 1 n + 1 I M F i ( t ) × I M F j ( t ) / x ( t ) 2
If the orthogonal index is small, the modal components are uncorrelated and no modal aliasing occurs; conversely, the modal aliasing is severe. Traditional threshold parameters are replaced by the orthogonal index to control the IMF screening process and to increase the frequency-intensive structure-identification ability of EMD decomposition. The orthogonal index IO reaches the minimum by the number of iterations of different IMF and is expressed as:
x ( N ) = arg min ( I O )
where N is the number of iterations per level IMF.
DFA [38,39] is adopted to determine the EMD threshold of the remaining IMF components to further filter out the noise. The formula for calculating the index α i is shown as follows:
α i = log 2 F ( s ) log 2 s = log 2 [ 1 2 N S v = 1 2 N s F 2 ( v , s ) ] 2 log 2 s
Theoretically, the corresponding sequences are white noise, pink noise, and brown noise when α i 0.5 , 0.5 < α i < 1 , and 1 α i < 1.5 , respectively. The noise of vibration response data of the discharge structure is mostly white noise. In this paper, only the white noise is considered, and the threshold for white noise is taken as 0 < α i < 0.5 . In addition, the noise in the prototype vibration signal of the arch dam is filtered with the improved wavelet threshold and EMD algorithm.

2.2. Modal Parameter Identification of High Arch Dam Discharge Structures Based on Improved HHT–RDT Algorithm

A modal identification method based on the HHT–RDT algorithm is proposed for the measured vibration response data of high arch dam structures in this paper. The HHT method is used to extract the modal response components of each order, and the RDT method is used to extract the feature information of each order component, so as to make up for the shortcomings of the two methods.
RDT is a processing method for extracting the free vibration response from the random response signal of structures. The method uses the property of the vibration signal mean value in the stationary stochastic process to identify the deterministic or random vibration signals and then extract the deterministic vibration signal.
For the linear system, assuming that X 1 ( t ) and X 2 ( t ) are the structural response data satisfying the gaussian distribution, according to the trigger condition, their random decrement cross-correlation function is:
C x i ( t i ) = [ c 1 X i ( t i ) c 2 ]
D x 2 x 1 ( τ ) = E [ X 1 ( t i + τ ) C x 1 ( t i ) ]
D x 2 x 1 ( τ ) = E [ X 2 ( t i + τ ) C x 1 ( t i ) ]
where E [ · ] represents the mathematical expectation of the cross-correlation function, which can be written as:
D x 1 x 1 ( τ ) = R x 1 x 1 ( τ ) σ 2 x 1 a ¯ R x 1 x 1 ( τ ) σ 2 x ˙ 1 b ¯
D x 2 x 1 ( τ ) = R x 2 x 1 ( τ ) σ 2 x 1 a ¯ R x 2 x 1 ( τ ) σ 2 x ˙ 1 b ¯
D x 2 x 1 ( τ ) = E [ X 2 ( t i + τ ) C x 1 ( t i ) ]
where R x i x 1 ( τ ) = E [ X i ( t ) X 1 ( t + τ ) ] , ( i = 1 , 2 ) , a ¯ = α α 2 x p x 1 ( x ) d x α α 2 p x 1 ( x ) d x , b ¯ = b b 2 x ˙ p x ˙ 1 ( x ˙ ) d x ˙ b b 2 p x ˙ 1 ( x ˙ ) d x ˙ .
For complex structures such as high arch dams, the structure composition is not a single linear system, but a multi-degree-of-freedom system. Therefore, the motion equation of arch dam structures under environmental excitation can be written as:
M X ¨ ( t ) + C X ˙ ( t ) + K X ( t ) = F ( t )
where M, C, and K are the mass matrix, damping matrix, and stiffness matrix of the arch dam structure, respectively; x is the displacement vector; and F(t) is the external excitation load vector.
Assuming that the external force is a stationary Gaussian stochastic process with a mean of zero, in the case where the mass, damping, and stiffness matrices are determined, R x x k is the homogeneous solution of the equation of motion and is shown as:
M R ¨ x x k ( τ ) + C R ˙ x x k ( τ ) + K R x x k ( τ ) = 0
The RDT function can be estimated as:
R ^ x j x k ( τ ) = 1 n i = 1 N x j ( t i + τ ) C x k ( t i )
where R ^ x j x k ( τ ) is the random decrement function of X j ( t ) for X k ( t ) ; N is the total number of trigger events; C x 1 ( t i ) is the specified punctual trigger condition; and ti is the ith moment.
In this paper, the waveform components of each order are obtained by the HHT method before identification by the RDT method. According to the approximate range of the main frequency reflected by the self-power spectrum of the signal, the pre-nth order natural frequency distribution of the system is obtained ( w 1 L < w 1 < w 1 R , w 2 L < w 2 < w 2 R , , w n L < w n < w n R ). If the i-th modal response component is solved, the input signal X j ( t ) is passed through a filter having a bandpass frequency of w i L < w i < w i R . After obtaining n waveform components, EMD is performed on each waveform component, and the obtained first-order natural mode function is used as a modal component. The modal response components obtained are stationary signals. The free-decay response function of each modal component is extracted by the RDT algorithm, which can be expressed as:
x ( t ) = e ξ w 0 t ( x 0 cos ( ω d t ) + x 0 + ξ ω 0 x 0 ω d sin ( ω d t ) )
where ω 0 is the natural frequency of the system; ξ is the relative damping coefficient; x 0 and x ˙ 0 are the initial displacement and velocity, respectively; and ω d = ω 0 1 ξ 2 is the damped natural frequency.
In addition, the expression can also be written as:
x ( t ) = A 0 e ξ w 0 t ( x 0 cos ( ω d t ) + φ 0 )
where A 0 is a constant, which is related to the excitation load intensity, structure mass, and frequency characteristics. HT is applied to obtain the analytical signal of x(t) as:
z ( t ) = x ( t ) + j H ( x ( t ) ) = A ( t ) e j θ ( t )
When the damping in the system is small, the amplitude and phase in the equation can be expressed as:
A ( t ) = A 0 e ξ ω 0 t ,      θ ( t ) = ω d t + φ 0
The amplitude and phase functions of the equation above are analyzed by logarithmic and differential operators, respectively, to obtain:
ln A ( t ) = ξ ω 0 t + ln A 0 ,      ω ( t ) = d θ ( t ) d t = ω d
The natural frequency ω 0 and damping ratio of the system can be obtained by ω d = ω 0 1 ξ 2 . Therefore, the natural frequency and damping ratio of the structure are obtained.

2.3. Recognition Process Based on Improved Wavelet Threshold–EMD and RDT Algorithm

Based on the above-mentioned theories, the construction procedure of modal parameter identification of high arch dam discharge structures based on improved wavelet threshold–EMD and RDT algorithm is shown in Figure 2.

3. Application Case

3.1. Engineering Data

The barrage is a typical high arch dam with a dam crest elevation of 1205 m and a maximum height of 240 m. To accurately obtain the working behavior of the arch dam, the measured points were arranged at the dam crest arch and the arch crown beam. Among them, the arch ring was arranged with seven horizontal dynamic displacement measuring points B1–B7 from the left bank to the right bank, and the crown beam was arranged with four horizontal dynamic displacement measuring points B8–B11 from the top to the bottom of the dam. The arrangement of arch dam prototype flood discharge vibration measurement points is shown in Figure 3, and the test conditions of the discharge vibration of the arch dam prototype are shown in Table 1. The dynamic displacement sensor was a DP seismic low-frequency displacement sensor with a frequency response range of 0.35–200 Hz and a sensitivity of 8–15 mv/μm. The DASP data acquisition and processing system was used for data acquisition. The sampling frequency was 200 Hz, the sampling time was 300 s, and the collected data points totaled 6000. The new method was demonstrated using the sampling data of the displacement monitoring data of the arch dam.
In the process of signal acquisition, the vibration response data of the discharge structure of the arch dam under discharge excitation easily mixed with some noises that covered up the useful signals, which will affect the accuracy of vibration response analysis of the test results. The time history signals and power spectra of typical B2, B3, B4, and B8 measuring points under working condition 2 are shown in Figure 4, Figure 5, Figure 6 and Figure 7. For the measured points B2 and B3 at the arch ring, the measuring points of the dam body near the dam shoulder had more vibration frequency components, and the vibration amplitude was smaller as it was closer to the dam shoulder. For the measuring points B4 and B8 of the arch crown beam, the vibration energy of the dam body of the high-frequency component was less than that of the low-frequency component, and the high-frequency peak was submerged by the low-frequency peak. Each measuring point only contained the first-order low-frequency component, and the vibration amplitude decreased with the decrease of elevation. In order to verify the rationality and accuracy of this method for noise reduction of arch dam discharge structures, the measuring points B2, B3, B4, and B8 are analyzed, respectively.

3.2. Improved Wavelet Threshold–EMD Hybrid Algorithm for Signal Noise Reduction

The original noise-containing signal and its normalized power spectral density curve are obtained, as shown in the Figure 3, on the basis of the measured displacement signal data of typical measurement point B2 under condition 2 of the arch dam. The time history diagram of the original signal and the power spectral density curve show that a large amount of white noise and low-frequency interference terms are mixed. The improved wavelet threshold–EMD hybrid algorithm was used to denoise the original signal. The time history diagram and the normalized power spectral density curve after filtering noise are shown in Figure 8.
A comparison of the dam time history diagram before and after filtering is shown in Figure 4; it can be seen that many burrs appear in the original dam time history diagram of the arch dam. After filtering by this method, the burrs of the new time history diagram were eliminated. In the original power spectrum, many interference signals such as low-frequency interference signals and white noise were found. In the filtered new power spectrum, low-frequency interference signals below 0.5 Hz have been eliminated, and white noise has also been removed.
In further verifying the effectiveness and superiority of the noise reduction method, wavelet threshold denoising, EMD denoising, and the improved wavelet–EMD hybrid algorithm are used to process the vibration response signal of the arch dam. The results are shown in Figure 9.
Figure 8 shows that the improved wavelet threshold–EMD hybrid algorithm works better than the two other methods of noise reduction. Wavelet threshold denoising can remove high-frequency noise but cannot remove low-frequency flow noise. As a result, the signal is subject to large losses, and no excellent filtering effect is achieved. The EMD method effectively removes the low-frequency water flow noise during filtering but causes modal aliasing easily, which results in difficulty in the modal identification of high arch dam structures. The improved wavelet threshold–EMD hybrid algorithm is used in the study: the wavelet threshold filters the high-frequency noise signal, inhibits the EMD endpoint effect, and reduces the modal aliasing. Then, EMD decomposition is used to further filter low-frequency flow noise and white noise. This method retains the useful signal components as much as possible and reflects the effective feature information of the structure. Therefore, the improved wavelet threshold–EMD hybrid algorithm is suitable for the noise reduction of the vibration signals of high arch dam discharge structures.

3.3. Modal Identification of Working Parameters of High Arch Dam Discharge Structures

After the measured vibration response is filtered, the combination of HHT and RDT is used to identify the working modal of the high arch dam discharge structure under discharge excitation. The power spectral density curve in Figure 4 has peaks in many places. Considering the absence of mode aliasing, the pass band is taken directly by band-pass filtering, and the pass bands are 1–1.5 Hz, 1.5–2 Hz, 2– 2.5 Hz, 2.5–3 Hz, 3–4 Hz, 4–5 Hz, 5–6 Hz, 6–8 Hz, and 8–10 Hz, respectively. After filtering, 10 components of Y1–Y10 are obtained. The power spectrum of these components is shown in Figure 10.
The free attenuation response of each component is extracted by RDT, the HHT and least squares fitting are used to obtain the relationship between amplitude and phase with time, and then the natural frequency and damping ratio of the structure are obtained. The time history diagram and power spectral of Y1 are shown in Figure 11a,b. The logarithmic amplitude curves of the free-decay response signal and the least squares fit are shown in Figure 11c,d.
The natural frequency and damping ratio of the first 10 orders of measurement points from B2, B3, B4, and B8 were obtained in sequence using the aforementioned method, and the recognition results are shown in Table 2. The recognition results show that the method in application can identify the natural frequency and damping ratio of the first 10 orders of the arch dam.
The recognition results indicate that the identification working frequencies of each order of each measuring point are very close, the frequency identification error is within 3%, and the damping ratio of each order is within 9%. To further explain the rationality of the modal identification results of the high arch dam, the modal parameters were calculated by the eigensystem realization algorithm (ERA) method [40], ARX [41] method and Ibrahim time domain (ITD) method [42]. The results for natural frequency are shown in Table 3.
Table 3 shows that the results obtained using the proposed method are more consistent than those obtained by other typical identification methods. The proposed method can accurately determine the natural frequency and damping ratio of the structure and can identify its higher modal order. ERA can accurately identify the first six natural frequencies of the structure; however, as the order increases, the ERA method cannot fully identify the frequency of the structure because of the difficulty in order determination, thereby affecting the work mode recognition results of the structure. The order of the ARX method is difficult to determine; thus, frequency mixing in the structure easily occurs, and distinguishing the characteristic frequency is difficult. These factors bring difficulties to the recognition. The ITD method can identify the frequency information accurately but cannot fit the higher-order modes. The proposed improved wavelet threshold–EMD and RDT algorithm involves two processes. The first process uses the improved wavelet threshold–EMD method to reduce noises of the original signal, retain the main characteristic information, and reduce the influence of modal aliasing. The second process applies the improved RDT algorithm to the working modal parameter identification of the high arch dam discharge structure. The band-pass filter is used to control the EMD process of the vibration response signal to obtain the modal components of the structure, the RDT method is used to extract the free attenuation information of each modal component, and the natural frequency and damping ratio of the first 10 orders of the structure are identified. This method avoids complex system ordering and frequency mixing in modal decomposition. The principle is simple and solving large matrices is not required during calculation. Therefore, the proposed method can be used for modal parameter identification under the discharge excitation response of high arch dams and can be applied to other engineering applications.

4. Conclusions

An improved modal threshold identification method for high arch dam discharge structures is proposed. This method can eliminate the high-frequency white noise and reduce the influence of the modal aliasing effectively and identify the modal parameters of high arch dam discharge structures accurately. The specific conclusions are as follows:
  • An improved wavelet threshold–EMD hybrid algorithm is proposed for noise reduction pretreatment on measured vibration response data of high arch dams. The improved wavelet threshold algorithm was adapted to overcome the defect of soft and hard threshold function selection, which can effectively eliminate the high-frequency white noise and reduce the influence of the modal aliasing. Then, EMD further eliminates the low-frequency flow noise and white noise, improving the accuracy of filtering and noise reduction.
  • An improved wavelet threshold–EMD and RDT algorithm is proposed for working mode identification of high arch dam discharge structures under the working environment load excitation. The proposed method avoids the complicated system ordering process and accurately identifies the modal parameters of high arch dam discharge structures. This method has a simple principle and does not require solving large matrices during calculation. The result has strong robustness and high identification accuracy.
  • The engineering examples show that the proposed method can accurately extract the work characteristic information of structures, has good noise reduction capabilities, and has high recognition accuracy. Therefore, this method eases the working modal parameter identification of high arch dam discharge structures and can be used for working modal recognition of other large frequency-intensive structures.

Author Contributions

Conceptualization, funding, B.W.; methodology, validation, formal analysis, resources, data curation, writing—original draft preparation, Y.G. and Z.Y.; visualization supervision, project administration, Y.G. and B.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by the National Natural Science Foundation of China (grant numbers 52169025, 51779115, 51869011), Science and Technology Project of the Jiangxi Provincial Department of Water Resources (grant numbers: 202224ZDKT21,202223YBKT43), Jiangxi Provincial Department of Education, with focus on science and technology projects (grant number: GJJ206001).

Institutional Review Board Statement

The study did not require ethical approval.

Informed Consent Statement

The study did not involve humans.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, H.; Wang, Y.; Wei, B. Inversion algorithm for the whole prototype dynamic displacement field of a high arch dam based on limited measuring points. J. Vib. Control. 2017, 23, 3431–3447. [Google Scholar] [CrossRef]
  2. Chengye, L.; Wei, L.; Bin, M.; Mingfu Song. Research on modal parameter identification method of high arch dam based on improved HHT. J. Hydroelectr. Eng. 2012, 31, 48–55. [Google Scholar]
  3. Cheng, L.; Tong, F. Application of Blind Source Separation Algorithms and Ambient Vibration Testing to the Health Monitoring of Concrete Dams. Math. Probl. Eng. 2016, 4, 1–15. [Google Scholar] [CrossRef] [Green Version]
  4. Bao, Y.; Xia, H.; Bao, Z.; Ke, S.; Li, Y. Structural Damage Identification of Pipe Based on GA and SCE-UA Algorithm. Math. Probl. Eng. 2013, 2013, 9. [Google Scholar] [CrossRef] [Green Version]
  5. Li, H.; Wang, G.; Wei, B.; Zhong, Y.; Zhan, L. Dynamic inversion method for the material parameters of a high arch dam and its foundation. Appl. Math. Model. 2019, 71, 60–76. [Google Scholar] [CrossRef]
  6. Wei, B.; Xie, B.; Li, H.; Zhong, Z.; You, Y. An improved Hilbert–Huang transform method for modal parameter identification of a high arch dam. Appl. Math. Model. 2021, 91, 297–310. [Google Scholar] [CrossRef]
  7. Menghui, X.; Nan, J. Dynamic load identification for interval structures under a presupposition of ‘being included prior to being measured. Appl. Math. Model. 2020, 85, 107–123. [Google Scholar]
  8. Wei, P.; Li, Q.; Sun, M.; Huang, J. Modal identification of high-rise buildings by combined scheme of improved empirical wavelet transform and Hilbert transform techniques. J. Build. Eng. 2022, 63, 105443. [Google Scholar] [CrossRef]
  9. Shangqing, Z.; Jiuzhang, L. Modal Test Analysis of Tied Arch Bridge Based on Random Subspace Method. Munic. Technol. 2015, 33, 38–41. [Google Scholar]
  10. Loh, C.H.; Wu, T.C. System identification of Fei-Tsui arch dam from forced vibration and seismic response data. J. Earthq. Eng. 2000, 4, 511–537. [Google Scholar] [CrossRef]
  11. Xueyuan, N.; Hua, D. Discussion on Modal Parameter Identification Method Based on Random Decreasing Technology. Mech. Des. 2012, 29, 1–5. [Google Scholar]
  12. Chun, R.; Cheng Cheng, Z.; Qi Feng, L. Application of HHT Method in Time Domain Identification of Structural Modal Parameters. J. Yangtze Univ. 2008, 5, 115–118, 414. [Google Scholar]
  13. Feldman, M. Time-varying vibration decomposition and analysis based on the Hilbert transform. J. Sound Vib. 2005, 295, 518–530. [Google Scholar] [CrossRef]
  14. Huokun, L.; Jijian, L. Opperational Modal Parameter Identification of High Arch Dam Under flood discharge exciting Based on Frequency Domain Decomposition Method. J. Vib. Shock. 2008, 27, 149–153. [Google Scholar]
  15. Li, H.; Liu, B.; Huang, W.; Liu, H.; Wang, G. Vibration load identification in the time-domain of high arch dam under discharge excitation based on hybrid LSQR algorithm. Mechanic 2020, 177, 109193. [Google Scholar] [CrossRef]
  16. Pengming, T.; Junyi, L.; Feng, Z. Application of Random Decrement Method in Identification of Working Modal Parameters. Mech. Electr. Eng. Technol. 2012, 41, 105–108. [Google Scholar]
  17. Wenge, T.; Xiaobing, Z. Working Mode Recognition Based on HVD and RDT. J. Guangxi Univ. 2018, 43, 132–140. [Google Scholar]
  18. Jianwei, Z.; Lianghuan, Z.; Qi, J.; Yu, Z.; Jia, G. Identification of working modal parameters of high dam discharge structure based on HHT. Vib. Meas. Diagn. 2015, 35, 777–783, 804. [Google Scholar]
  19. He, X.H.; Hua, X.G.; Chen, Z.Q.; Huang, F.L. EMD-based random decrement technique for modal parameter identification of an existing railway bridge. Eng. Struct. 2011, 33, 1348–1356. [Google Scholar] [CrossRef]
  20. Frank Pai, P.; Palazotto, A.N. HHT-based nonlinear signal processing method for parametric and non-parametric identification of dynamical systems. Int. J. Mech. Sci. 2008, 50, 1619–1635. [Google Scholar] [CrossRef]
  21. Jianping, H.; Lin, L.; Hongtao, W.; Wei, Q. Modal parameter identification based on Hilbert-Huang transform and random decrement technique. World Earthq. Eng. 2011, 27, 72–77. [Google Scholar]
  22. Jinyuan, T.; Weitao, C.; Siyu, C.; Wei, Z. A new wavelet threshold function and its application in vibration signal denoising analysis. J. Vib. Shock. 2009, 28, 118–121. [Google Scholar]
  23. Zhisong, L. Signal Denoising Method Based on Wavelet Analysis. J. Zhejiang Ocean. Univ. 2011, 30, 150–154. [Google Scholar]
  24. Xiuping, L.; Wei, Y.; Lili, H.; Junfeng, J.; Jian, X.; Zebin, S.; Haifeng, S. X-ray pulsar signal de-noising for impulse noise using wavelet packet. Aerosp. Sci. Technol. 2017, 64, 147–153. [Google Scholar] [CrossRef]
  25. Wenwen, Z.; Xingwen, C. A New EMD Denoising Method. Electron. Technol. 2008, 21, 30–32, 36. [Google Scholar]
  26. Dybała, J.; Zimroz, R. Rolling bearing diagnosing method based on Empirical Mode Decomposition of machine vibration signal. Appl. Acoust. 2014, 77, 195–203. [Google Scholar] [CrossRef]
  27. Chengye, L.; Jijian, L.; Wei, L.; Bin, M. Improvement of Joint Filtering Method of EMD and Wavelet Threshold and Its Application in Vibration Analysis of Discharge Structure. J. Vib. Shock. 2013, 32, 63–70, 110. [Google Scholar]
  28. Jintao, Y.; Shuyan, Z.; Wei, W. Denoising of Transmitted Signal Based on Empirical Mode Decomposition and Wavelet Transform. J. Harbin Inst. Technol. 2011, 43, 88–92. [Google Scholar]
  29. Aijun, H.; Jingjing, S.; Ling, X. Modal aliasing problem in empirical mode decomposition. Vib. Test. Diagn. 2011, 43, 429–434. [Google Scholar]
  30. Jianwei, Z.; Qi, J.; Yu, Z.; Lianghuan, Z.; Jia, G. A Signal Denoising Method Suitable for Vibration Analysis of Discharge Structures. J. Vib. Shock. 2015, 34, 179–184. [Google Scholar]
  31. Zhi, C.; Feng, H.; Xiaolong, Z.; Xu, J.; Zailong, H. Wavelet-experience mode decomposition method for trend term modulation for atmospheric dry length profile denoising. Acta Opt. Sin. 2017, 37, 15–26. [Google Scholar]
  32. Wei, B.W.; Yuan, D.Y.; Li, H.K.; Xu, Z. Combination forecast model for concrete dam displacement considering residual correction. Struct. Health Monit. 2017, 18, 232–244. [Google Scholar] [CrossRef]
  33. Wei, B.W.; Yuan, D.Y.; Li, H.K.; Li, L. Modified hybrid forecast model considering cjhaotic residual errors for dam deformation. Struct. Control. Health Monitor. 2018, 25, e2188. [Google Scholar] [CrossRef]
  34. Wei, L.; Xiaohui, C.; Haijie, M. Research on adaptive determination of decomposition layer in wavelet threshold denoising algorithm. Comput. Simul. 2009, 26, 311–313, 336. [Google Scholar]
  35. Yang, J.; Lia, P.; You, Y.; Dian, X. An improved EMD method for modal identification and a combined static-dynamic method for damage detection. J. Sound Vib. 2018, 28, 242–260. [Google Scholar] [CrossRef]
  36. Qinwei, R.; Wei, L.; Ce, S.; Hongzhi, K. Research on EMD dense modal identification and its application in power plant buildings. J. Hydroelectr. Eng. 2017, 36, 103–113. [Google Scholar]
  37. Mert, A.; Akan, A. Detrended fluctuation thresholding for empirical mode decomposition based denoising. Digit. Signal Proc. 2014, 32, 48–56. [Google Scholar] [CrossRef]
  38. Wenping, H.; Qiong, W.; Haiying, C.; Wen, Z. Comparison of different filtering methods for denoising in detrended wave analysis. Acta Phys. Sin. 2011, 60, 029203. [Google Scholar]
  39. Fubin, L.; Xiangming, G. Hurst Exponent Estimation Based on EEMD and DFA. Meas. Control. Technol. 2013, 32, 98–101. [Google Scholar]
  40. Jijian, L.; Hunkun, L.; Jianwei, Z. Hydraulic structure vibration mode ERA recognition method based on singular entropy fixed-stage noise reduction. Chin. Sci. 2008, 38, 1398–1413. [Google Scholar]
  41. Wen, Y.; Guoming, L.; Wei, L. Modal analysis of nozzle concrete gravity dam based on strong earthquake record. Earthq. Eng. Eng. Vib. 2014, 34, 25–33. [Google Scholar]
  42. Shengkui, D.; Deqiang, Z.; Jian, L. Structural modal parameter identification based on ITD method. J. Civ. Eng. Manag. 2012, 29, 15–19. [Google Scholar]
Figure 1. Prototype vibration response of the arch dam.
Figure 1. Prototype vibration response of the arch dam.
Water 14 03735 g001
Figure 2. Flow chart of modal parameter identification of high arch dam discharge structures.
Figure 2. Flow chart of modal parameter identification of high arch dam discharge structures.
Water 14 03735 g002
Figure 3. Arrangement of arch dam prototype flood discharge vibration measurement points. (a) The arch dam prototype; (b) Sensor arrangement for the arch dam vibration test; (c) DP-type dynamic displacement sensor; (d) DASP signal acquisition system.
Figure 3. Arrangement of arch dam prototype flood discharge vibration measurement points. (a) The arch dam prototype; (b) Sensor arrangement for the arch dam vibration test; (c) DP-type dynamic displacement sensor; (d) DASP signal acquisition system.
Water 14 03735 g003
Figure 4. Time history signals and power spectra of B2 measuring point under condition 2. (a) Time history signals; (b) Normalized power spectra.
Figure 4. Time history signals and power spectra of B2 measuring point under condition 2. (a) Time history signals; (b) Normalized power spectra.
Water 14 03735 g004
Figure 5. Time history signals and power spectra of B3 measuring point under condition 2. (a) Time history signals; (b) Normalized power spectra.
Figure 5. Time history signals and power spectra of B3 measuring point under condition 2. (a) Time history signals; (b) Normalized power spectra.
Water 14 03735 g005
Figure 6. Time history signals and power spectra of B4 measuring point under condition 2. (a) Time history signals; (b) Normalized power spectra.
Figure 6. Time history signals and power spectra of B4 measuring point under condition 2. (a) Time history signals; (b) Normalized power spectra.
Water 14 03735 g006
Figure 7. Time history signals and power spectra of B8 measuring point under condition 2. (a) Time history signals; (b) Normalized power spectra.
Figure 7. Time history signals and power spectra of B8 measuring point under condition 2. (a) Time history signals; (b) Normalized power spectra.
Water 14 03735 g007
Figure 8. Time history diagram and the normalized power spectrum diagram after noise reduction. (a) Time history signals; (b) Normalized power spectra.
Figure 8. Time history diagram and the normalized power spectrum diagram after noise reduction. (a) Time history signals; (b) Normalized power spectra.
Water 14 03735 g008
Figure 9. Comparison of calculation results of three noise reduction methods. (a) Time history signals; (b) Normalized power spectra.
Figure 9. Comparison of calculation results of three noise reduction methods. (a) Time history signals; (b) Normalized power spectra.
Water 14 03735 g009
Figure 10. Power spectral of Y1~Y10 components.
Figure 10. Power spectral of Y1~Y10 components.
Water 14 03735 g010
Figure 11. Process of modal parameters identification of Y1 component. (a) Time history; (b) Power spectra; (c) Free attenuation response; (d) Logarithmic amplitude curve.
Figure 11. Process of modal parameters identification of Y1 component. (a) Time history; (b) Power spectra; (c) Free attenuation response; (d) Logarithmic amplitude curve.
Water 14 03735 g011
Table 1. Test conditions of discharge vibration of arch dam prototype.
Table 1. Test conditions of discharge vibration of arch dam prototype.
ConditionSpillway OpeningUpstream
Water Level (m)
Downstream
Water Level (m)
1Fourth surface spillway and first tunnel spillway1196.001014.36
2Third, fourth surface spillway and first tunnel spillway1196.011014.50
Table 2. Comparison of modal identification results of different measuring points of the high arch dam structure.
Table 2. Comparison of modal identification results of different measuring points of the high arch dam structure.
ModalNatural Frequency/HzDamping Ratio/%
Order B2B3B4B8B2B3B4B8
11.441.441.441.431.451.791.511.69
1.91
21.521.531.521.531.541.081.261.35
32.212.192.192.201.451.282.241.46
42.802.812.782.811.481.462.171.28
53.583.603.593.611.313.474.081.79
64.324.314.324.284.384.966.55.6
74.904.904.924.912.181.933.672.38
85.155.125.205.212.865.52.876.72
97.247.287.307.267.677.168.127.47
109.679.699.589.524.796.456.455.34
Table 3. Modal frequency identification of the high arch dam discharge structure with different methods.
Table 3. Modal frequency identification of the high arch dam discharge structure with different methods.
Modal OrderMethod in the Study/HzERA/HzARX/HzITD/Hz
11.441.451.421.43
21.521.531.471.51
32.212.222.142.08
42.802.872.812.79
53.583.743.683.63
64.324.384.404.32
74.94.824.714.85
85.155.095.03/
97.247.366.987.01
109.679.529.16/
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, Y.; You, Z.; Wei, B. Working Mode Identification Method for High Arch Dam Discharge Structure Based on Improved Wavelet Threshold–EMD and RDT Algorithm. Water 2022, 14, 3735. https://doi.org/10.3390/w14223735

AMA Style

Guo Y, You Z, Wei B. Working Mode Identification Method for High Arch Dam Discharge Structure Based on Improved Wavelet Threshold–EMD and RDT Algorithm. Water. 2022; 14(22):3735. https://doi.org/10.3390/w14223735

Chicago/Turabian Style

Guo, Yingjia, Zongzhe You, and Bowen Wei. 2022. "Working Mode Identification Method for High Arch Dam Discharge Structure Based on Improved Wavelet Threshold–EMD and RDT Algorithm" Water 14, no. 22: 3735. https://doi.org/10.3390/w14223735

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