| 利用MSMPCA去噪的CEEMD方法监测高频GNSS同震形变 |
| |
| 引用本文： | 李艳艳,殷海涛,韩林桥.利用MSMPCA去噪的CEEMD方法监测高频GNSS同震形变[J].武汉大学学报(信息科学版),2022,47(3):352-360. |
| |
| 作者姓名： | 李艳艳 殷海涛 韩林桥 |
| |
| 作者单位： | 1.山东省地震局，山东 济南，250014 |
| |
| 基金项目： | 国家自然科学基金41804001山东省自然科学基金ZR2019BD006山东省自然科学基金ZR2017LDD009 |
| |
| 摘 要： | 考虑到全球导航卫星系统(global navigation satellite system,GNSS)地震信号的非线性和非平稳性,利用一种多尺度多方向主成分分析(multiscale multiway principal component analysis,MSMPCA)去噪的完备总体经验模态分解(complete...
|
| 关 键 词： | 高频全球导航定位系统 经验模态分解 主成分分析 同震位移 |
| 收稿时间： | 2019-09-24 |
A Denoising Method of MSMPCA Based on CEEMD for Coseismic Deformation Monitoring with High-Frequency GNSS |
| |
| Affiliation: | 1.Shandong Earthquake Agency, Jinan 250014, China2.College of Geomatics, Shandong University of Science and Technology, Qingdao 266590, China |
| |
| Abstract: | Objectives The global navigation satellite system(GNSS) coseismic displacements are always subject to temporal and spatial colored errors, which may lead to the incorrect interpretation of some geophysical phenomenon. Therefore, considering the nonlinear and nonstationary of GNSS seismic signal, a denoising method of multiscale multiway principal components analysis(MSMPCA) based on complete ensemble empirical mode decomposition(CEEMD) (C-MSMPCA) is used to denoising the high-frequency GNSS coseismic displacements for obtaining high precise instantaneous coseismic deformation. Methods The C-MSMPCA method integrates the merits of CEEMD and MSMPCA to denoise the high-frequency GNSS coseismic displacements. Firstly, CEEMD replaces wavelet transform (WT) to decompose GNSS coseismic signal into various intrinsic mode function(IMF). Then, instead of directly deleting the previous high-requency IMF, all IMFs are grouped based on their frequency bands (FB) and each group of modes with the same FB is denoised using MSMPCA. Finally, the processed IMF is reconstructed and handled using multiway principal component analysis(MPCA) to obtain the denoised coseismic signal. Results High-frequency(5 Hz) GNSS observations of 14 stations of the California real-time GNSS network (CRTN) in southern California during the 2010 EI Mayor-Cucapah earthquake(Mw 7.2) are denoised using C-MSMPCA. The performance of C-MSMPCA is compared with that of MSMPCA. The results show that: (1) Take coseismic displacement time series before and after filtering at P496 station as example, as the amplitude of the seismic signal decreases, in terms of the peak displacements and the long-term stability, the coseismic displacements obtained from C-MSMPCA are better than those obtained from MSMPCA. C-MSMPCA not only eliminates the noise in the GNSS seismic displacements, but also effectively retains the seismic signal. (2) Take GNSS residual and error curve time series before and after filtering at P496 station, as observed in the 600 s before the earthquake as example, the raw residual and error curve time series have centimeter-level low-frequency colored noise and high-frequency random noise. Compared with the error curves time series obtained from MSMPCA, the fluctuations of those obtained from C-MSMPCA are more obvious, indicating that the performance of C-MSMPCA is better than MSMPCA. (3) Each root mean square error(RMSE) and average RMSE of the original and denoised residual signals of the 14 GNSS stations, as observed in the 600 s before the earthquake show that RMSE of C-MSMPCA is less than those of MSMPCA, regardless of the coordinate components. The average RMSE values of the original GNSS solutions are determined to be 1.82, 2.55, and 7.16 mm for the east, north, and vertical components, respectively, while those of the denoising results obtained by C-MSMPCA are 0.89, 0.96, and 4.27 mm, respectively, representing RMSE reductions of 51.10%, 62.35%, and 40.36%, respectively. (4) The results of power spectral analysis show that high-frequency GNSS displacements contain high-frequency white noise and low-frequency colored noise. MSMPCA and C-MSMPCA both can significantly eliminate high-frequency white noise at 1×101—1×102 s periods. C-MSMPCA has the better performance than MSMPCA at 1×102—1×103 s and 1×103—1×104 s periods. Over all, the denoising results of MSMPCA are smoother than those of C-MSMPCA, but the latter can effectively retains higher-frequency seismic signal. Conclusions Denoising is an essential step in earthquake-induced coseismic displacements. Considering the nonlinear and nonstationary of GNSS seismic signal and how to effectively preserve seismic signal, a denoising method of C-MSMPCA is used to denoising the high-frequency GNSS coseismic displacements. Over all, the performance of C-MSMPCA is better than that of MSMPCA. The high performance of C-MSMPCA is that the method is not directly remove the first high-frequency IMF but group IMF based on their frequency bands for adaptively MSMPCA processing. The method avoids the loss of high-frequency effective information and effectively preserves high-frequency seismic waveform of high-frequency GNSS displacements, which is particularly critical for the research and application of earthquake-induced coseismic displacements. |
| |
| Keywords: | |
|
| 点击此处可从《武汉大学学报(信息科学版)》浏览原始摘要信息 |
|
点击此处可从《武汉大学学报(信息科学版)》下载全文 |
|
