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
-
[1] 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 [2] 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 [3] 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 [4] 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 [5] 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 [6] 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 [7] 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 [8] 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 [9] 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 [10] 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. [11] 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 [12] 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 [13] Kouba, J., 2003. Measuring seismic waves induced by large earthquakes with GPS. Studia Geophysical Geodaetica, 47(4): 741-755. doi: 10.1023/A:1026390618355 [14] 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 [15] Larson, K.M., 2009. GPS seismology. Journal of Geodesy, 83(3-4): 227-233. doi: 10.1007/s00190-008-0233-x [16] 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 [17] 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 [18] 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 [19] 戴吾蛟, 丁晓利, 朱建军, 等, 2006. 基于经验模式分解的滤波去噪法及其在GPS多路径效应中的应用. 测绘学报, 35(4): 321-327. doi: 10.3321/j.issn:1001-1595.2006.04.005 [20] 方荣新, 施闯, 辜声峰, 等, 2009. 基于PPP动态定位技术的同震地表形变分析. 武汉大学学报(信息科学版), 34(11): 1340-1358. https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200911021.htm [21] 方颖, 江在森, 张晶, 等, 2009. 青藏高原现今高分辨率地壳运动特征. 地球科学——中国地质大学学报, 34(3): 399-404. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200903003.htm [22] 黄丁发, 丁晓利, 陈永奇, 等, 2001. GPS多路径效应影响与结构振动的小波滤波筛分方法. 测绘学报, 30(1): 36-41. doi: 10.3321/j.issn:1001-1595.2001.01.008 [23] 钟萍, 丁晓利, 郑大伟, 等, 2006. Vondrak滤波法用于结构振动与GPS多路径信号的分离. 中南大学学报(自然科学版), 36(6): 1189-1195. doi: 10.3969/j.issn.1672-7207.2006.06.033 -
点击查看大图
图(7) / 表(1)
计量
- 文章访问数: 181
- HTML全文浏览量: 131
- PDF下载量: 3
- 被引次数: 0
