利用小波变换提取面波群速度

蒋婵君, 王有学, 熊彬, 贠鹏, 徐金荣, 乃振龙

蒋婵君, 王有学, 熊彬, 贠鹏, 徐金荣, 乃振龙. 2017: 利用小波变换提取面波群速度. 地震学报, 39(3): 356-366. DOI: 10.11939/jass.2017.03.005
引用本文: 蒋婵君, 王有学, 熊彬, 贠鹏, 徐金荣, 乃振龙. 2017: 利用小波变换提取面波群速度. 地震学报, 39(3): 356-366. DOI: 10.11939/jass.2017.03.005
Jiang Chanjun, Wang Youxue, Xiong Bin, Yun Peng, Xu Jinrong, Nai Zhenlong. 2017: Measurement of surface wave group velocity using wavelet transform. Acta Seismologica Sinica, 39(3): 356-366. DOI: 10.11939/jass.2017.03.005
Citation: Jiang Chanjun, Wang Youxue, Xiong Bin, Yun Peng, Xu Jinrong, Nai Zhenlong. 2017: Measurement of surface wave group velocity using wavelet transform. Acta Seismologica Sinica, 39(3): 356-366. DOI: 10.11939/jass.2017.03.005

利用小波变换提取面波群速度

基金项目: 

大陆构造与动力学国家重点实验室开放基金 K201505

广西自然科学基金重点项目 2016GXNSFDA380014

国家自然科学基金(41574039)、大陆构造与动力学国家重点实验室开放基金(K201505) 和广西自然科学基金重点项目(2016GXNSFDA380014) 共同资助

国家自然科学基金 41574039

详细信息
    通讯作者:

    王有学, e-mail: uxue.wang@glut.edu.cn

  • 中图分类号: P315.3+1

Measurement of surface wave group velocity using wavelet transform

  • 摘要: 选取了几种常见的小波母函数,分别提取了同一理论下的面波数据的群速度,并与理论群速度进行对比,结果表明Morlet小波提取面波群速度的效果最好.此外,将Morlet小波与常用的多重滤波提取群速度的结果进行了比较,结果表明: ① 多重滤波法非常依赖高斯滤波系数α的取值,α的取值应随面波周期的增大而减小;② 在α取值得当的前提下,在20—35 s周期范围内多重滤波法提取面波群速度的相对误差比Morlet小波小,在周期大于35 s时,两者相对误差相近; ③ 合适的α值的选取需在不同周期段耗费大量时间进行大量试验,这说明多重滤波法不具备自适应性;而采用小波变换分析短周期信号时,时间窗变窄,频率窗变长,当分析长周期信号时,时间窗变长,频率窗变窄,具有对信号的自适应性,这是小波变换相比多重滤波法的最大优点.
    Abstract: Wavelet transform is proposed to measure surface wave group velocity in this paper. Several common mother wavelets are respectively used to calculate Rayleigh wave group velocities of the same synthetic seismograms, the results are compared with the theoretical group velocities. The comparison results indicate that Morlet wavelet is the most effective one in calculating group velocities. Furthermore, we compare Morlet wavelet transform with the most common used multiple filter technique in calculating group velocity, from which the following conclusions are drawn: ① multiple filter analysis method extremely depends on the value of Gaussain filter parameter α, and smaller value of α should be chosen during long period of surface wave; ② when the value of α is chosen adaptively, multiple filter analysis method is more accurate than Morlet wavelet transform in calculating group velocities when the period of surface wave is between 20 seconds and 35 seconds, while multiple filter analysis method has the similar accuracy with Morlet wavelet transform when the period of surface wave is over 35 seconds; ③ multiple filter analysis method will cost much more time to choose adaptive value of Gaussain filter parameter α for different periods without self-adaptivity to signals; however, wavelet transform can narrow the time window and lengthen the frequency window when analyzing short period signals and lengthen the time window and narrow the frequency window when analyzing long period signals with self-adaptivity to signals, which shows the outstanding advantage of wavelet transform over the multiple filter analysis method.
  • 图  1   常见小波波形图

    Figure  1.   Waveforms of some common wavelets

    图  2   理论瑞雷面波垂向分量波形(a)和振幅谱(b)

    Figure  2.   Vertical component synthetic seismogram of Rayleigh wave (a) and its amplitude spectra (b)

    图  3   不同小波母函数提取的面波群速度与理论群速度对比

    Figure  3.   Comparison between surface wave group velocities obtained from different wavelet generating functions with theoretical values

    图  4   不同小波母函数提取面波群速度的相对误差

    Figure  4.   Relative errors of surface wave group velocities obtained from different wavelet generating functions

    图  5   Morlet小波和多重滤波提取面波群速度与理论面波群速度对比

    Figure  5.   Comparison between surface wave group velocities obtained from Morlet wavelet transform and multiple filter technique with theoretical values

    图  6   Morlet小波和多重滤波法提取面波群速度的相对误差

    Figure  6.   Relative errors of surface wave group velocities obtained by Morlet wavelet transform and multiple filter technique

    图  7   南岭西部地区某台站记录到的MW7.1地震垂直分量瑞雷面波波形(a)和振幅谱(b)

    Figure  7.   Vertical component seismogram of Rayleigh wave (a) and its amplitude spectra (b) of MW7.1 earthquake recordod by one station in western Nanling region

    图  8   多重滤波法和Morlet小波提取南岭西部地区某台站的面波群速度

    Figure  8.   Surface wave group velocities of one station in western Nanling region obtained from multiple filter technique and Morlet wavelet transform

  • 曹小林, 张雪梅, 朱介寿. 1999.中国及邻区的地震面波频散反演[J].成都理工学院学报, 26(3): 295-298. https://www.researchgate.net/publication/298494957_Surface_wave_dispersion_inversion_in_China_and_adjacent_regions

    Cao X L, Zhang X M, Zhu J S. 1999. Surface wave dispersion inversion in China and adjacent regions[J]. Journal of Chengdu University of Technology, 26(3): 295-298 (in Chinese). https://www.researchgate.net/publication/298494957_Surface_wave_dispersion_inversion_in_China_and_adjacent_regions

    陈浩朋, 朱良保, 王清东, 杨颖航, 张攀. 2014.面波群速度测量误差分析[J].地球物理学进展, 29(5): 1985-1995. doi: 10.6038/pg20140501

    Chen H P, Zhu L B, Wang Q D, Yang Y H, Zhang P. 2014. An error analysis of surface wave group velocity measurement[J]. Progress in Geophysics, 29(5): 1985-1995 (in Chinese). doi: 10.6038/pg20140501

    陈立, 薛梅, Plon L K, 杨挺. 2012.南海瑞雷面波群速度层析成像及其地球动力学意义[J].地震学报, 34(6): 754-772. https://www.researchgate.net/publication/272503488_Group_velocity_tomography_of_Rayleigh_waves_in_South_China_Sea_and_its_geodynamic_implications

    Chen L, Xue M, Plon L K, Yang T. 2012. Group velocity tomography of Rayleigh waves in South China Sea and its geodynamic implications[J]. Acta Seismologica Sinica, 34(6): 754-772 (in Chinese). https://www.researchgate.net/publication/272503488_Group_velocity_tomography_of_Rayleigh_waves_in_South_China_Sea_and_its_geodynamic_implications

    房立华, 吴建平, 王未来, 王长在, 杨婷. 2013.华北地区勒夫波噪声层析成像研究[J].地球物理学报, 56(7): 2268-2277. doi: 10.6038/cjg20130714

    Fang L H, Wu J P, Wang W L, Wang C Z, Yang T. 2013. Love wave tomography from ambient seismic noise in North-China[J]. Chinese Journal of Geophysics, 56(7): 2268-2277 (in Chinese). doi: 10.6038/cjg20130714

    何正勤, 丁志峰, 叶太兰, 孙为国, 张乃铃. 2001.中国大陆及其邻域地壳上地幔速度结构的面波层析成像研究[J].地震学报, 23(6): 596-603. http://mall.cnki.net/magazine/Article/DZXY200106004.htm

    He Z Q, Ding Z F, Ye T L, Sun W G, Zhang N L. 2001. Surface wave tomography of the crust and upper mantle of Chinese mainland[J]. Acta Seismologica Sinica, 23(6): 596-603 (in Chinese). http://mall.cnki.net/magazine/Article/DZXY200106004.htm

    潘佳铁, 吴庆举, 李永华, 张风雪, 张广成. 2011.华北地区瑞雷面波相速度层析成像[J].地球物理学报, 54(1): 67-76. http://manu39.magtech.com.cn/Geophy/EN/abstract/abstract7718.shtml

    Pan J T, Wu Q J, Li Y H, Zhang F X, Zhang G C. 2011. Rayleigh wave tomography of the phase velocity in North China[J]. Chinese Journal of Geophysics, 54(1): 67-76 (in Chinese). http://manu39.magtech.com.cn/Geophy/EN/abstract/abstract7718.shtml

    王敏玲, 陈赟, 梁晓峰, 徐逸鹤, 范宇婷, 徐涛, 张忠杰, 滕吉文. 2015.华南及南海北部地区瑞利面波层析成像[J].地球物理学报, 58(6): 1963-1975. doi: 10.6038/cjg20150612

    Wang M L, Chen Y, Liang X F, Xu Y H, Fan Y T, Xu T, Zhang Z J, Teng J W. 2015. Surface wave tomography for South China and the northern South China Sea area[J]. Chinese Journal of Geophysics, 58(6): 1963-1975 (in Chinese). doi: 10.6038/cjg20150612

    余大新, 吴庆举, 王鹏, 叶庆东, 潘佳铁, 高孟潭. 2016.蒙古中南部地区基于天然地震的勒夫波相速度层析成像[J].地震学报, 38(1): 41-52. doi: 10.11939/jass.2016.01.004

    Yu D X, Wu Q J, Wang P, Ye Q D, Pan J T, Gao M T. 2016. Love wave phase velocity tomography in the south-central Mongolia from earthquakes[J]. Acta Seismologica Sinica, 38(1): 41-52 (in Chinese). doi: 10.11939/jass.2016.01.004

    周青云, 何永峰, 靳平, 廖桂生. 2006.利用多重滤波方法提取面波频散曲线[J].西北地震学报, 28(1): 46-50. http://en.cnki.com.cn/Article_en/CJFDTOTAL-ZBDZ200601009.htm

    Zhou Q Y, He Y F, Jin P, Liao G S. 2006. Using MFT obtain Rayleigh-wave dispersion curve[J]. Northwestern Seismological Journal, 28(1): 46-50 (in Chinese). http://en.cnki.com.cn/Article_en/CJFDTOTAL-ZBDZ200601009.htm

    朱介寿, 曹家敏, 蔡学林, 严忠琼, 曹小林. 2002.东亚及西太平洋边缘海高分辨率面波层析成像[J].地球物理学报, 45(5): 646-664. http://manu39.magtech.com.cn/Geophy/EN/Y2002/V45/I05/646

    Zhu J S, Cao J M, Cai X L, Yan Z Q, Cao X L. 2002. High resolution surface wave tomography in East Asia and West Pacific marginal seas[J]. Chinese Journal of Geophysics, 45(5): 646-664 (in Chinese). http://manu39.magtech.com.cn/Geophy/EN/Y2002/V45/I05/646

    朱良保, 熊安丽. 2007.面波频散测量的频时分析法[J].地震地磁观测与研究, 28(1): 1-13. http://en.cnki.com.cn/Article_en/CJFDTOTAL-DZGJ200701000.htm

    Zhu L B, Xiong A L. 2007. The frequency-time analysis method for surface wave dispersion measurement[J]. Seismological and Geomagnetic Observation and Research, 28(1): 1-13 (in Chinese). http://en.cnki.com.cn/Article_en/CJFDTOTAL-DZGJ200701000.htm

    Cara M. 1973. Filtering of dispersed wavetrains[J]. Geophys J Int, 33(1): 65-80. doi: 10.1111/j.1365-246X.1973.tb03415.x

    Daubechies I. 1988. Orthonormal bases of compactly supported wavelets[J]. Commun Pure Appl Math, 41(7): 909-996. doi: 10.1002/(ISSN)1097-0312

    Daubechies I. 1990. The wavelet transform, time-frequency localization and signal analysis[J]. IEEE Trans Inf Theory, 36(5): 961-1005. doi: 10.1109/18.57199

    Dziewonski A, Bloch S, Landisman M. 1969. A technique for the analysis of transient seismic signals[J]. Bull Seismol Soc Am, 59(1): 427-444. http://citeseerx.ist.psu.edu/showciting?cid=973577

    Herrmann R B, Ammon C J. 2004. Computer programs in seismology, 3.30[CP/OL]. [2013-08-20]. http://www.eas.slu.edu/eqc/eqccps.html.

    Inston H H, Marshall P D, Blamey C. 1971. Optimization of filter bandwidth in spectral analysis of wavetrains[J]. Geophys J Int, 23(2): 243-250. doi: 10.1111/gji.1971.23.issue-2

    Kennett B L N, Engdahl E R, Buland R. 1995. Constraints on seismic velocities in the Earth from traveltimes[J]. Geophys J Int, 122(1): 108-124. doi: 10.1111/gji.1995.122.issue-1

    Kolínský P. 2004. Surface wave dispersion curves of Eurasian earthquakes: The SVAL program[J]. Acta Geodyn Geomater, 1(2): 165-185. https://www.researchgate.net/publication/295315410_Surface_waves_dispersion_curves_of_eurasian_earthquakes_The_SVAL_program

    Landisman M, Dziewonski A, Satȏ Y. 1969. Recent improvements in the analysis of surface wave observations[J]. Geophys J Int, 17(4): 369-403. doi: 10.1111/gji.1969.17.issue-4

    Levshin A, Pisarenko V F, Pogrebinsky G A. 1972. On a frequency-time analysis of oscillations[J]. Ann Geophys, 28(2): 211-218. https://www.researchgate.net/publication/284699529_On_a_frequency-time_analysis_of_oscillations

    Morlet J, Arens G, Fourgeau E, Giard D. 1982. Wave propagation and sampling theory, Part Ⅱ: Sampling theory and complex waves[J]. Geophysics, 47(2): 222-236. doi: 10.1190/1.1441329

    Nyman D C, Landisman M. 1977. The display-equalized filter for frequency-time analysis[J]. Bull Seismol Soc Am, 67(2): 393-404. http://www.bssaonline.org/content/67/2/393.short

    Wu Q J, Zheng X F, Pan J T, Zhang F X, Zhang G C. 2009. Measurement of interstation phase velocity by wavelet transformation[J]. Earthquake Science, 22(4): 425-429. doi: 10.1007/s11589-009-0425-3

图(8)
计量
  • 文章访问数:  931
  • HTML全文浏览量:  554
  • PDF下载量:  44
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-08-01
  • 修回日期:  2016-11-18
  • 发布日期:  2017-04-30

目录

    /

    返回文章
    返回