Extraction of High-Rate GPS Seismic Wave Signals with Wavelet Packets Decomposition
-
摘要: 随着高采样率GPS接收机的出现、高精度的定轨以及数据处理技术的发展, 利用GPS观测高频率、瞬态的地震波信号成为可能.但如何消除混叠在地震波信号中的多路径、随机噪声等误差, 有效地提取地震波信号, 仍然是制约高采样率GPS及其地球物理应用的重要因素.提出一种基于小波包分解的方法, 对动态位移序列中的多路径误差进行消除, 同时去除高频率随机噪声, 提取地震波信号.通过结合SCIGN的19个GPS测站的1Hz采样GPS观测数据, 对2010年墨西哥M7.2地震的地震波引起的地表动态位移进行解算, 采用小波包分解有效地提取地震波信号并对其进行谱分析.结果表明, 该方法提取的地震波信号能较好地反映出地震波的传播及其特性, 具有无需多天重复观测、处理精度高等特点.Abstract: With the development of high-rate GPS receivers, precise orbit and processing technology of GPS data, it is possible to observe the high frequency, transient dynamic displacement by GPS. However, it remains a problem as how to mitigate the error such as multipath errors, and random noise aliasing in geophysical signals so as to extract seismic signals, which in turn limits the high-rate GPS and its geophysical applications. In this paper, an approach based on wavelet packets decomposition (WPD) is presented to extract seismic signals through mitigating the multipath error and random noise of dynamic displacement series from high-rate GPS. With the 1 Hz observation data from 19 stations in Southern California Integrated GPS Network (SCIGN) during the Mexico M7.2 earthquake in 2010, the ground displacement is calculated. Meanwhile, the approach based on WPD is introduced for seismic signal extraction and spectrum analysis. As is shown in results, the approach is accurate and effective in seismic signal extraction to reflect the characteristics of seismic wave propagations and it enjoys an advantage that it does not necessarily involve multiple-day observation.
-
图 1 三层小波包分解子频带示意
Fig. 1. Frequency bands of wavelet packets decomposition of three levels
图 3 测站PIN1在N、E方向初步动态位移序列结果分量
Fig. 3. North and east raw dynamic displacement components of PIN1 station
图 4 测站PIN1在N、E方向初步动态位移序列功率谱图
Fig. 4. Power spectrum of north, east raw dynamic displacement components of PIN1 station
图 5 小波包分解提取的测站PIN1的多路径误差
Fig. 5. Multipath error extracted by wavelet packet decomposition in PIN1 station
图 6 测站PIN1的地震动态位移序列N、E方向分量
Fig. 6. Clean north and east components of seismic dynamic displacement of PIN1 station
表 1 19个测站的震中距及其由1 Hz采样观测值解算的E、N方向坐标的参考误差
Table 1. Epicentral distance, north and east coordinates reference error based on 1 Hz observation of 19 stations
站名 震中距(km) 处理前误差(mm) 处理后误差(mm) 站名 震中距(km) 处理前误差(mm) 处理后误差(mm) N E N E N E N E BILL 215 10.8 8.3 5.0 3.3 PIN1 183 7.4 7.9 4.5 3.1 CCCS 301 11.5 7.5 4.7 3.2 SIO5 195 12.0 7.2 4.4 2.9 CIT1 334 26.3 7.3 5.3 3.6 TABL 320 17.1 20.1 4.8 3.5 CRBT 638 18.3 10.2 5.2 3.3 TRAK 278 18.9 18.9 4.3 2.9 CSDH 328 23.3 6.0 4.8 3.1 USGC 110 11.9 12.1 5.5 3.5 DHLG 130 13.4 7.4 4.0 2.7 WIDC 208 25.8 11.2 4.7 3.2 EWPP 290 11.7 6.5 5.3 3.5 AZRY 184 51.4 171.0 10.3 8.8 GLRS 110 18.7 8.8 4.9 3.3 CSN1 371 580.1 305.7 10.8 7.8 MASW 615 10.1 7.0 4.9 3.1 LOWS 627 - - - - MONP 109 18.1 9.3 4.3 2.9 - - - - - - 
下载: 导出CSV
-
Bilich, A., Cassidy, J.F., Larson, K.M., et al., 2008. GPS seismology: application to the 2002 Mw=7.9 Denali fault earthquake. Bulletin of Seismological Society of America, 98(2): 593-606. doi: 10.1785/0120070096 Bock, Y., Prawirodirdjo, L., Melboume, T.I., 2004. Detection of arbitrarily large dynamic ground motions with a dense high-rate GPS network. Geophysical Research Letters, 31(6): L06604. doi: 10.1029/2003GL019150 Chang, C., Juang, J., 2008. An adaptive multipath mitigation filter for GNSS applications. Eurasip Journal on Advances in Signal Processing, 214815. doi: 10.1155/2008/214815 Choi, K.H., Bilich, A., Larson, K.M., et al., 2004. Modified sidereal filtering: implications for high-rate GPS positioning. Geophysical Research Letters, 31(22): L22608. doi: 10.1029/2004GL021621 Dai, W.J., Ding, X.L., Zhu, J.J., et al., 2006. EMD filter method and its application in GPS multipath. Acta Geodatica et Cartographica Sinica, 35(4): 321-327 (in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-CHXB200604004.htm Fang, R.X., Shi, C., Gu, S.F., 2009. Precise point positioning with high-rate GPS data applied to seismic displacements analysis. Geomatics and Information Science of Wuhan University, 34(11): 1340-1358(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-WHCH200911021.htm Fang, Y., Jiang, Z.S., Zhang, J., et al., 2009. High-resolution characteristics of recent crustal movement in Qinghai-Tibet plateau. Earth Science—Journal of China University of Geosciences, 34(3): 399-404(in Chinese with English abstract). doi: 10.3799/dqkx.2009.044 Gomberg, J., Bodin, P., Larson, K., et al., 2004. Earthquake nucleation by transient deformations caused by the M=7.9 Denali, Alaska, earthquake. Nature, 427(6975): 621-624. doi: 10.1038/nature02335 Genrich, J.F., Bock, Y., 1992. Rapid resolution of crustal motion at short ranges with the global positioning system. Journal of Geophysical Research, 97(B3): 3261-3269. doi: 10.1029/91JB02997 Hardebeck, J.L., Boatwright, J., Dreger, D., et al., 2005. Preliminary report on the 22 December 2003, M6.5 San Simeon, California earthquake. Seismological Research Letters, 75(2): 155-172. Huang, D.F., Ding, X.L., Chen, Y.Q., et al., 2001. Wavelet filters based separation of GPS multi-path effects and engineering structure vibrations. Acta Geodatica et Cartographica Sinica, 30(1): 36-41(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-CHXB200101007.htm Ji, C., Larson, K.M., Tan, Y., et al., 2004. Slip history of the 2003 San Simeon earthquake constrained by combining 1 Hz GPS, strong motion, and teleseismic data. Geophysical Research Letters, 31(17): L17608. doi: 10.1029/2004GL020448 Kouba, J., 2003. Measuring seismic waves induced by large earthquakes with GPS. Studia Geophysical Geodaetica, 47(4): 741-755. doi: 10.1023/A:1026390618355 Larson, K.M., Bodin, P., Gomberg, J., et al., 2003. Using 1 Hz GPS to measure deformations caused by the Denali fault earthquake. Science, 300(5624): 1421-1424. doi: 10.1126/sciences.1084531 Larson, K.M., 2009. GPS seismology. Journal of Geodesy, 83(3-4): 227-233. doi: 10.1007/s00190-008-0233-x Miyazaki, S., Larson, K.M., Choi, K.H., et al., 2004. Modeling the rupture process of the 2003 September 25 Tokachi-Oki(Hokkaido) earthquake using 1 Hz GPS data. Geophysical Research Letters, 31(21): L21603. doi: 10.1029/2004GL021457/full Xia, L.Y., Liu, J.N., 2001. Approach for multipath reduction using wavelet algorithm. ION GPS2001: 2134-2143. http://www.ion.org/search/view_abstract.cfm?jp=p&idno=1882 Zhong, P., Ding, X.L., Zheng, D.W., et al., 2006. Separation of structure vibrations and GPS multipath signals using Vondrak filter. J. Cent. South. Univ. (Science and Technology), 36(6): 1189-1195 (in Chinese with English abstract). http://www.researchgate.net/publication/288165914_Separation_of_structural_vibrations_and_GPS_multipath_signals_using_Vondrak_filter 戴吾蛟, 丁晓利, 朱建军, 等, 2006. 基于经验模式分解的滤波去噪法及其在GPS多路径效应中的应用. 测绘学报, 35(4): 321-327. doi: 10.3321/j.issn:1001-1595.2006.04.005 方荣新, 施闯, 辜声峰, 等, 2009. 基于PPP动态定位技术的同震地表形变分析. 武汉大学学报(信息科学版), 34(11): 1340-1358. https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200911021.htm 方颖, 江在森, 张晶, 等, 2009. 青藏高原现今高分辨率地壳运动特征. 地球科学——中国地质大学学报, 34(3): 399-404. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200903003.htm 黄丁发, 丁晓利, 陈永奇, 等, 2001. GPS多路径效应影响与结构振动的小波滤波筛分方法. 测绘学报, 30(1): 36-41. doi: 10.3321/j.issn:1001-1595.2001.01.008 钟萍, 丁晓利, 郑大伟, 等, 2006. Vondrak滤波法用于结构振动与GPS多路径信号的分离. 中南大学学报(自然科学版), 36(6): 1189-1195. doi: 10.3969/j.issn.1672-7207.2006.06.033 -
点击查看大图
图(7) / 表(1)
计量
- 文章访问数: 453
- HTML全文浏览量: 471
- PDF下载量: 12
- 被引次数: 0
