1. Introduction
Vibration is a key indicator of the health of infrastructure such as bridges and buildings [
1,
2,
3,
4]. There are many methods to measure vibration, which fall into two main categories: contact and non-contact. Contact methods use sensors such as displacement, velocity, and acceleration sensors that need to be attached to the vibrating point. However, these methods may not be feasible in some situations due to the constraints of the installation location. Non-contact methods, such as lidar and radar, do not require installation. They can be placed at a suitable distance and adjust their beams to target the vibrating points. Lidar has high precision, but it is susceptible to environmental factors such as rain and fog. Radar is immune to these factors and can operate in any weather condition. Radar works by transmitting microwave signals and receiving the target echo. It analyzes the phase change of the echo signal over time to obtain the vibration deformation, which has very high submillimeter-level measurement accuracy [
5,
6]. Therefore, radar has become more popular in the field of infrastructure vibration monitoring due to its advantages of being all-day, all-weather, high-precision, and non-contact.
Radar-based vibration measurement has been a research topic for over 20 years. Tarchi et al. first proposed using radar interferometry to monitor buildings in 2000 [
7]. Piraccini et al. developed and tested the first radar system for bridge dynamic monitoring in 2004 [
8]. The University of Florence and IDS (Ingegneria dei Sistemi SpA) collaborated to develop the IBIS series of systems in 2007 [
9], one of which, the IBIS-S system, became widely used for infrastructure vibration monitoring [
1,
3,
5,
10,
11]. This system operates in the Ku band with a power of 40 W. Vibration radar systems have improved in recent years. In 2020, the Beijing Institute of Technology created a multi-channel radar that can be used for imaging and vibration information acquisition of large infrastructures [
12,
13,
14]. Since 2018, small and lightweight millimeter-wave radars have been used for vibration monitoring. The Chinese Academy of Sciences made a Ka-band radar system for measuring bridge vibration [
15]. Chengkun Jiang et al. from Tsinghua University measured the vibration of industrial systems using a commercial millimeter-wave board, the Texas Instruments (TI) IWR1642 BoosterPack [
16]. The IWR1642 operates at 77 GHZ frequency with a power of only 17.8 mv. The authors of [
17,
18] used a TI millimeter-wave MIMO board to monitor bridges and vehicles. Low-cost millimeter-wave radars may offer new possibilities for infrastructure vibration monitoring.
Radar-based vibration measurement is challenging. The echo’s differential phase contains the deformation information, but it is easily affected by noise and clutter [
2,
12,
19,
20,
21]. The radar transmission power must meet certain requirements, and the observation distance is usually short. When the vibrating point has weak scattering, the strong corner reflector at the measured point can improve the signal-to-noise ratio (SNR) and signal-to-clutter ratio (SCR). The authors of [
1,
22] used this method of installing corner reflectors. The authors of [
11] did not install a corner reflector when measuring the bridge bottom deformation, because the support structure under the bridge is a good corner reflector with strong scattering. This helps extract deformation information accurately. However, installing corner reflectors makes the measurement task more difficult and inconvenient. When the observation distance is long, the radar power is low, or the vibrating target has weak scattering, the echo signal has a low SNR. This causes large differential phase noise and even jump errors [
23]. Under low SNR conditions, the phase wrapping issue randomly generates jump errors in differential phase that cannot be unwrapped using the phase continuity feature. Moreover, environmental clutter adds to the signal, so that the differential phase does not accurately reflect deformation information.
To address these challenges, researchers have used signal processing to improve the SNR and SCR of signals. For clutter suppression, the most classic method is the circle fitting method, which [
12,
16,
21,
22,
24] used. In the complex plane, the radar signal forms an arc or a circle, and the clutter causes the circle center to deviate from the origin. By estimating the circle center, the signal origin can return to the circle center to eliminate the clutter effect. For noise suppression, the classic algorithm is a low-pass filtering method. Different papers use low-pass filters at different stages of the processing process. The authors of [
21] applied low-pass filtering to the original echo data at the beginning, while [
6,
22] applied low-pass filtering to the final deformation to remove high-frequency noise. The former can solve the phase jump problem under low SNR and estimate the clutter more accurately, but it requires oversampling. The sampling rate needs to be several times the Doppler bandwidth, and the harmonic characteristics of the Doppler signal make its bandwidth several times the vibration frequency [
9], which demands higher data acquisition and produces more data. However, the latter does not require oversampling because the bandwidth of the deformation variable is smaller than that of the original signal, but it cannot solve the phase jump problem and is not conducive to clutter estimation. The author used the TI radar board for vibration monitoring and found that due to low radar power and a long observation distance, random jump errors occurred in the deformation inversion. We proposed an Additive Constant method to solve this problem in [
23]. This method adds a constant much larger than the vibration signal to the original radar signal and limits the signal to the right plane of the complex plane, thus avoiding phase wrapping and solving the phase jump problem. However, this method uses the small-amplitude assumption, i.e., that the vibration amplitude is much smaller than the signal wavelength, which limits its application. In addition, this method cannot solve the clutter estimation problem at low SNR.
In this paper, we propose a periodic filtering method to suppress the noise of the original echo signal by using the periodic repetition characteristics of vibration signals. Unlike the low-pass filtering method, which uses a sliding window to average adjacent cells, periodic filtering accumulates signals from several adjacent cycles to suppress noise. This method does not require oversampling for pulse repetition frequency (PRF) and only needs to meet the Nyquist sampling criterion. After periodic filtering, the stationary clutter can be estimated more accurately, and the jump error problem in deformation can also be solved.
The rest of this paper is organized as follows.
Section 2 describes the geometry and the signal model of the vibration monitoring radar.
Section 3 analyzes the influences of static clutter and noise on deformation inversion.
Section 4 introduces the basic principle of periodic filtering, gives the main flow of the algorithm, and derives the expressions of the periodic filter in time and frequency domains.
Section 5 verifies the effectiveness and correctness of the proposed method through simulation and real experiments.
Section 6 presents the conclusion.
2. Geometry and Signal Model of Vibration Monitoring Radar
The observation geometry of the vibration monitoring radar is shown in
Figure 1. The radar is fixed. It transmits and receives signals at constant time intervals. The radar beam illuminates the vibration point
. The range from
to the antenna phase center is denoted as
.
In this paper, the radar system uses Frequency Modulated Continuous Wave (FMCW). The transmitted signal is expressed as follows:
where
is the fast time,
is the center frequency, and
is the signal frequency modulation rate.
The received signal is a delayed version of the transmitted signal. The transmitted and received signals are then mixed, generating the baseband signal:
where
is the slow time,
is the complex reflectivity of the vibration point
, and
is the time delay of the signal from transmission to being reflected by the vibration point
and finally being received by the radar receiver, which can be expressed as follows:
where
is the speed of light and
represents the simple harmonic vibration, which can be expressed as a function of slow time:
where
is the vibration amplitude,
is the vibration frequency, and
is the initial phase of the simple harmonic vibration.
In Equation (2), the quadratic phase term is called the residual video phase (RVP). It is small, and vibration measurement applications can safely ignore it without significant impact on deformation inversion. Substitute Equation (3) into Equation (2) and define a new frequency according to the linear time-frequency relationship of FMCW as follows:
Then, Equation (2) can be expressed as follows:
An inverse Fourier transform with respect to
is performed on Equation (6) to obtain the range compressed signal, which can be expressed as follows:
where
is the time corresponding to frequency
and
is the signal bandwidth.
The peak value of the range compressed signal is located at range
. Vibration information is contained in the one-dimensional signal at distance
, which changes with slow time
. The one-dimensional signal can be expressed as follows:
where
is the wavelength. Equation (8) is the vibration signal model under ideal conditions. To simplify the following analysis, the phase of
is represented by
, which can be expressed as follows:
3. Influence of Static Clutter and Noise on Vibration Radar Signals
The real vibration signal contains the clutter from the environment and the noise generated by the system. This paper discusses the influence of static clutter and Gaussian white noise on vibration signals. The vibration signal containing static clutter and Gaussian white noise can be expressed as follows:
where
represents the static clutter in the same range cell with the vibrating target, is a complex constant that does not change with the slow time, and
represents Gaussian white noise.
In this paper, the static clutter signal model is used. In vibration monitoring, the position of the radar is fixed, and the static clutter signal model can be used when the background is stationary.
There are mainly two situations when the background is not stationary: (1) there are strong targets, such as vehicles or pedestrians passing by, causing fast-changing clutter; and (2) radar position or radar transmitted signal change slowly, resulting in slowly changing clutter.
Fast-changing clutter has a high Doppler frequency, while slowly-changing clutter has a Doppler frequency close to zero. Both types of non-static clutter can be suppressed in the doppler domain range because they have different doppler frequencies from vibration signals. The proposed periodic filter has a comb-shaped spectrum and can filter out non-static clutter, as shown in
Section 4.2.
In the following section, the characteristics of the vibration signal and the influence of static clutter and noise on the vibration signal will be analyzed in the complex plane.
3.1. Vibration Signal in a Complex Plane
As shown in
Figure 2a, an ideal vibration signal is a phase-modulated signal. The real part and the imaginary part form a circle or an arc on the complex plane, with the center at the origin. The radius of the circle or arc is the reflectivity of the vibrating point, that is, the absolute value of the complex reflectivity
, and the phase is the vibration signal phase
. According to the vibration signal phase shown in Equation (9), when the amplitude
, the vibration signal forms a complete circle in the complex plane. Conversely, when the amplitude
, the vibration signal forms an arc in the complex plane.
3.2. Influence of Static Clutter on Vibration Signal
As shown in
Figure 3, in the complex plane, the presence of static clutter makes the center of the circle or arc deviate from the origin. The horizontal and vertical coordinates of the new center are the real part and the imaginary part of the static clutter, respectively. The presence of static clutter will make the signal phase difference smaller than the actual vibration phase difference, resulting in an underestimation of the vibration amplitude. In this figure,
is the phase difference of the signal, and
is the actual vibration phase difference. The circle fitting method is commonly used to estimate the center of the complex plane signal so that the static clutter can be estimated and removed.
3.3. Influence of Noise on the Vibration Signal
When the SNR is low, the strong Gaussian white noise makes the center estimation inaccurate, especially when the arc is only a small portion of a circle, as shown in
Figure 4a.
In addition, under the condition of low SNR, there will be a jump error in deformation inversion. Because of strong noise, the differential phase would be wrapped, and the ambiguity number is unknown. As shown in
Figure 4b, the true differential phase is
. With the influence of noise, the differential phase is
. The differential phase is wrapped in the range
. Because of the randomness of the noise, it is difficult to solve this wrapping problem; that is, the phase ambiguity number is unknown, which will lead to phase jump errors.
In this paper, we define it as a “low SNR” condition when the jump error occurs. Vibration frequency extraction has a low demand for SNR because the signal is accumulated in the Doppler domain. However, the radar signal phase is quite sensitive to noise, and the deformation extraction has a relatively high demand for SNR.
4. Clutter Suppression and Deformation Inversion Based on Periodic Filtering
To reduce the influence of noise on clutter estimation and deformation inversion, the traditional method is to perform low-pass filtering on the original vibration radar signal shown in Equation (10). However, the low-pass filtering method requires a high PRF. In the Doppler domain range, the vibration signal spectrum contains multiple harmonics with the target vibration frequency as the base frequency. Therefore, the Doppler bandwidth is usually several times the target vibration frequency. The bandwidth of the low-pass filter needs to be greater than the Doppler bandwidth; otherwise, the Doppler spectrum would be damaged, resulting in signal distortion and inaccurate deformation inversion. In order to achieve a good filtering effect, the PRF needs to be several times the bandwidth of the low-pass filter. When PRF is limited, it is difficult to achieve a good denoising effect without causing signal distortion.
This section will first derive the range Doppler domain expression of the vibration signal and give the Doppler bandwidth expression. Then, the basic principle of periodic filtering is introduced, and the expressions of periodic filtering in the time domain and frequency domain are derived. Next, the circle-fitting clutter estimation method and the deformation inversion method are briefly introduced. Finally, the algorithm flow of this paper is given.
4.1. Doppler Domain Vibration Signal
The Doppler domain signal is the Fourier transform of the one-dimensional slow time signal shown in Equation (8). The Doppler spectrum reflects the frequency characteristics of the vibrating target. Because the vibration signal is periodic, its spectrum is discrete, and the frequency components include the vibration frequency and its harmonic frequencies. The Doppler domain signal of a vibrating target can be expressed as follows:
Since the Doppler spectrum is discrete, only the spectrum at the harmonic frequency
need to be calculated, where
, and only the signal within one cycle needs to be integrated. The expression is as follows:
where
is the Bessel function of order
.
In the Doppler domain, the energy of the periodic signal is accumulated at the harmonic frequency, so the SNR of the spectrum is very high. The base frequency can be easily obtained by peak detection of the spectrum, which provides the input parameter for the design of the periodic filter.
The Doppler bandwidth can be calculated using the phase modulation property of the vibration radar signal. The instantaneous Doppler frequency can be obtained by the derivative of the signal phase shown in Equation (9) with respect to slow time, which can be expressed as follows:
When
, the Doppler frequency reaches the maximum, that is:
When
, the Doppler bandwidth of the echo signal is several times the target vibration frequency. According to the Nyquist sampling criteria, the PRF should be greater than twice the maximum Doppler frequency, that is:
4.2. Periodic Filtering Method
The periodic filtering method proposed in this paper makes use of the periodic repetition characteristics of vibration signals to accumulate the signals of several adjacent cycles to suppress noise. Compared with the low-pass filtering method, this method has a low requirement for PRF and only needs to meet the Nyquist sampling criteria.
The expression of time-domain periodic filtering is:
where symbol
represents convolution and
is the periodic filter, which can be expressed as follows:
where
represents the impulse function,
is the vibration period,
is the order number of adjacent cycles involved in filtering, and
is manually set.
The physical meaning of Equation (16) is that the signals of adjacent cycles are averaged, and the SNR is increased by times. The greater the , the better the noise suppression effect.
Time-domain convolution is equivalent to frequency domain multiplication. The frequency-domain expression of the periodic filter can be expressed as follows:
The expression of frequency-domain periodic filtering can be expressed as follows:
The figures of the frequency domain low-pass filter and periodic filter are shown in
Figure 5. Since the filter is symmetrical about zero Doppler, only the positive axis of Doppler frequency is shown in the figure. The low-pass filter is a rectangular window, and the periodic filter is a group of narrow pulses with
as the interval. The periodic filter has a comb-shaped spectrum; the larger the
, the narrower the pulses, and the better the filtering effect. The comb-shaped filter can filter out noise and non-static clutter.
Two main factors should be considered when selecting the value of M. The first factor is SNR. The M value needs to be large enough to avoid jumping errors in deformation inversion. The second factor is the time-variation characteristics of the vibration signal. If the vibration frequency, vibration amplitude, and clutter stay constant or change slowly with time, the value of M can be larger, and vice versa.
In this paper, we make use of the periodic repetition of a vibration signal and do not use the small vibration amplitude assumption. For both cases, the proposed periodic filtering method is effective whether the vibration amplitude is greater than or less than a quarter of the wavelength.
4.3. Circle-Fitting Clutter Estimation
After periodic filtering, the signal SNR increases. The next step is to estimate the clutter. The least squares circle fitting method is usually used to estimate the center and radius of the complex plane vibration signal [
21]. The center of the circle represents the clutter, and the radius represents the scattering coefficient of the vibrating target. Let the center of the circle be
,
, the radius be
,
, and the vibration signal be
, where
denotes the real part of the signal,
denotes the imaginary part, the symbol
represents matrix conjugate transpose,
is the discrete value of the slow time, and
. The least square circle fitting is to minimize the following cost function:
where
To express Equation (20) in the form of a matrix operation, let:
additionally,
Substitute Equations (21)–(24) into Equation (20), and the cost function can be rewritten as follows:
then the parameter to be estimated is:
The solution of the above equation is:
The estimated circle center and radius can be expressed as follows:
where
,
,
are the three elements of vector
, and
,
are the horizontal and vertical coordinates of the circle center, respectively.
The estimated clutter signal is:
4.4. Deformation Inversion
After the static clutter is estimated and removed from the periodic filtered signal, the estimated vibration signal is:
Deformation inversion is applied to the estimated vibration signal. To avoid phase wrapping, the differential phase is obtained by interference processing first, and then the deformation is obtained by accumulating the differential phase.
The expression of the differential phase is as follows:
where the function
is to obtain the phase of a complex value. The deformation has the expression of:
Periodic filtering has removed part of the noise and solved the problem of differential phase jump errors caused by noise. The deformation has a frequency of
, so its bandwidth is much narrower than that of the original vibration radar signal. In other words, deformation inversion reduces the signal bandwidth. Therefore, the deformation can be further filtered by a low-pass filter, such as a time-domain Butterworth filter or a frequency-domain rectangular filter. Taking frequency domain low-pass filtering as an example, its expression is:
where FT represents the Fourier transform, IFT represents the inverse Fourier transform, and
represents the frequency domain low-pass filter, which has the expression of:
where
represents a rectangular function and
is the bandwidth of the low-pass filter, which needs to be greater than twice the vibration frequency, that is:
4.5. Algorithm Flow chart
The main steps of the algorithm are shown in
Figure 6, including:
Vibration frequency estimation. The vibration frequency is the base frequency of the Doppler signal . It can be extracted by peak detection, and it is the input parameter for the design of the periodic filter;
Periodic filtering. The frequency-domain periodic filter is shown in Equation (18). The signal is multiplied with to obtain the filtered signal , whose SNR is improved;
Clutter estimation. The least squares circle fitting method is used to estimate the clutter of the periodically filtered signal , and then the estimated clutter is subtracted from to obtain the estimated vibration signal ;
Deformation inversion. The differential phase of is obtained by interference processing, and then the deformation is obtained by accumulating the differential phase;
Low-pass filtering. The deformation has a frequency of , and a bandwidth that is much smaller than the Doppler signal. The noise in the deformation can be further filtered by low-pass filtering.
6. Conclusions
In this paper, we propose a periodic filtering method for vibration radar signals and derive the expressions of the periodic filter in time and frequency domains. This method uses the periodic repetition characteristics of a vibration signal to accumulate several adjacent cycles and denoise the signal. Compared with the traditional low-pass filtering method, this method does not require oversampling for PRF. We verify the correctness and effectiveness of this method by simulation experiments, vibrating calibrator experiments, and bridge vibration monitoring experiments. The results show that when the PRF is not oversampled, the low-pass filtering method widely used in the existing literature will distort the vibration signal, so that the vibration deformation cannot be accurately retrieved. The proposed periodic filtering method, which does not require oversampling of the PRF, can significantly improve the SNR, thereby solving the problems of inaccurate clutter estimation and jump errors in deformation inversion. Therefore, parameters such as vibration frequency and amplitude can be accurately estimated.
This method has some limitations: it can only extract one frequency of vibration deformation from the signal at a time; for complex vibrations with multiple frequencies, each frequency has to be extracted separately. Our future work focuses on two aspects: (1) deriving a complex vibration filter for multiple-frequency vibrations and expanding the applicability of the method; and (2) studying the millimeter-wave multi-channel radar vibration measurement method. This system can image large infrastructures and acquire vibration frequencies and amplitudes at each scattering center. Combining imaging and vibration measurement facilitates vibration point localization. We aim to improve and optimize the proposed periodic filtering method and apply it to millimeter-wave multi-channel radar for high-precision vibration deformation extraction of large-scale infrastructure.