接收函数曲波变换去噪与偏移成像

陈一方, 陈九辉, 郭飚, 齐少华, 赵盼盼. 2019. 接收函数曲波变换去噪与偏移成像. 地球物理学报, 62(6): 2027-2037, doi: 10.6038/cjg2019M0248
引用本文: 陈一方, 陈九辉, 郭飚, 齐少华, 赵盼盼. 2019. 接收函数曲波变换去噪与偏移成像. 地球物理学报, 62(6): 2027-2037, doi: 10.6038/cjg2019M0248
CHEN YiFang, CHEN JiuHui, GUO Biao, QI ShaoHua, ZHAO PanPan. 2019. Denoising the receiver function through curvelet transforming and migration imaging. Chinese Journal of Geophysics (in Chinese), 62(6): 2027-2037, doi: 10.6038/cjg2019M0248
Citation: CHEN YiFang, CHEN JiuHui, GUO Biao, QI ShaoHua, ZHAO PanPan. 2019. Denoising the receiver function through curvelet transforming and migration imaging. Chinese Journal of Geophysics (in Chinese), 62(6): 2027-2037, doi: 10.6038/cjg2019M0248

接收函数曲波变换去噪与偏移成像

  • 基金项目:

    国家自然科学基金项目(41590862)和(41274060)共同资助

详细信息
    作者简介:

    陈一方, 男, 1991年生, 博士研究生, 主要从事宽频带地震学研究.E-mail:yifangchen_9134@hotmail.com

    通讯作者: 陈九辉, 男, 1969年生, 研究员, 博士生导师, 主要从事宽频带地震学和深部探测研究.E-mail:chenjh@ies.ac.en
  • 中图分类号: P315

Denoising the receiver function through curvelet transforming and migration imaging

More Information
  • 增强接收函数偏移图像的垂向分辨率意味着提高参与叠加的接收函数的频率,但是采用高频接收函数通常伴随着对接收函数质量和参考速度模型的更高要求.通过叠加处理可去除部分接收函数中的随机噪声干扰,但同一台站的接收函数之间经常存在难以通过简单叠加消除的噪声信号.压制接收函数随机噪声的干扰可加强成像效果和提高图像分辨率,对推进叠加偏移成像质量的提高有重要的实际意义.本文利用在川西地区布设的31个流动台站所记录的远震波形数据,使用曲波变换去噪后信噪比增强的接收函数进行共转换点叠加(CCP),获得沿北纬31°线下方800 km深度范围内速度间断面图像.研究结果表明:(1)对接收函数进行曲波变换去噪,可压制随机噪声,增强转换震相的追踪性,提高数据信噪比;(2)通过去噪处理,大幅提高接收函数用于偏移成像的主频率;(3)偏移结果确认了接收函数反演得到的松潘和川滇块体下方具有厚度约10~20 km的过渡性Moho的认识;(4)上地幔过渡带的结果预示在龙门山断裂带以西的小范围内有可能存在下地壳或上地幔物质的拆沉.

  • 加载中
  • 图 1 

    空间域和频率域的曲波

    Figure 1. 

    Curvelets in spatial and frequency domains

    图 2 

    接收函数曲波去噪测试

    Figure 2. 

    Test of receiver function by using curvelet transform

    图 3 

    文中所用台站及分布

    Figure 3. 

    Station distribution

    图 4 

    曲波变换去噪m值选取测试

    Figure 4. 

    Selection of value of m in the curvelet filters

    图 5 

    KCD08台站滤波至1 s接收函数曲波变换去噪

    Figure 5. 

    Noise attenuation of 1 Hz receiver functions with curvelet transform at Station KCD08

    图 6 

    KCD08台站2 s接收函数曲波变换去噪

    Figure 6. 

    Noise attenuation of 0.5 Hz receiver functions with curvelet transform at Station KCD08

    图 7 

    周期4 s接收函数偏移成像

    Figure 7. 

    Migrated image with 0.25 Hz receiver function

    图 8 

    周期2 s接收函数偏移成像

    Figure 8. 

    Migrated image with 0.5 Hz receiver function

    图 9 

    31°N剖面地壳上地幔间断面结构

    Figure 9. 

    Discontinuities structure of crust and upper mantle along the profile of 31°N

  •  

    Ammon C J, Randall G E, Zandt G. 1990. On the nonuniqueness of receiver function inversion. Journal of Geophysical Research:Solid Earth, 95(B10):15303-15318. doi: 10.1029/JB095iB10p15303

     

    Anderson D L. 1989. Theory of the Earth. Boston:Blackwell Scientific Publication, 366.

     

    Bostock M G, Rondenay S. 1999. Migration of scattered teleseismic body waves. Geophysical Journal International, 137(3):732-746. doi: 10.1046/j.1365-246x.1999.00813.x

     

    Cammarano F, Goes S, Vacher P, et al. 2003. Inferring upper-mantle temperatures from seismic velocities. Physics of the Earth and Planetary Interiors, 138(3-4):197-222. doi: 10.1016/S0031-9201(03)00156-0

     

    Candès E, Donoho D. 2000. Curvelets-a surprisingly effective nonadaptive representation for objects with edges.//Cohen A, Rabut C, Schumaker L eds. Curves and Surface Fitting: Saint-Malo. Nashville: Vanderbilt University Press, 105-120.

     

    Candès E, Donoho D. 2004. New tight frames of curvelets and optimal representations of objects with piecewise C2 singularities. Communications on Pure and Applied Mathematics, 57(2):219-266. doi: 10.1002/cpa.v57:2

     

    Candès E J, Romberg J, Tao T. 2006. Robust uncertainty principles:exact signal reconstruction from highly incomplete frequency information. IEEE Transactions on Information Theory, 52(2):489-509. doi: 10.1109/TIT.2005.862083

     

    Chen L, Wen L X, Zheng T Y. 2005. A wave equation migration method for receiver function imaging:1. Theory. Journal of Geophysical Research:Solid Earth, 110(B11):B11309, doi:10.1029/2005JB003665.

     

    Dueker K G, Sheehan A F. 1997. Mantle discontinuity structure from midpoint stacks of converted P to S waves across the Yellowstone hotspot track. Journal of Geophysical Research:Solid Earth, 102(B4):8313-8327. doi: 10.1029/96JB03857

     

    Guo B, Liu Q Y, Chen J H, et al. 2009. Teleseismic P-wave tomography of the crust and upper mantle in Longmenshan area, west Sichuan. Chinese Journal of Geophysics (in Chinese), 52(2):346-355. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxb200902006

     

    Helffrich G. 2000. Topography of the transition zone seismic discontinuities. Reviews of Geophysics, 38(1):141-158, doi:10.1029/1999RG000060.

     

    Herrmann F J, Verschuur E. 2004. Curvelet-domain multiple elimination with sparseness constraints.//74th Ann. Internat Mtg., Soc. Expi. Geophys.. Expanded Abstracts, 1333-1336.

     

    Herrmann F J, Böniger U, Verschuur D J. 2007. Non-linear primary-multiple separation with directional curvelet frames. Geophysical Journal International, 170(2):781-799. doi: 10.1111/gji.2007.170.issue-2

     

    Julià J, Ammon C J, Herrmann R B, et al. 2000. Joint inversion of receiver function and surface wave dispersion observations. Geophysical Journal International, 143(1):99-112. doi: 10.1046/j.1365-246x.2000.00217.x

     

    Kind R, Yuan X, Saul J, et al. 2002. Seismic images of crust and upper mantle beneath Tibet:Evidence for Eurasian plate subduction. Science, 298(5596):1219-1221. doi: 10.1126/science.1078115

     

    Langston C A. 1979. Structure under Mount Rainier, Washington, inferred from teleseismic body waves. Journal of Geophysical Research:Solid Earth, 84(B9):4749-4762. doi: 10.1029/JB084iB09p04749

     

    Lei J S, Zhao D P, Su J R, et al. 2009. Fine seismic structure under the Longmenshan fault zone and the mechanism of the large Wenchuan earthquake. Chinese Journal of Geophysics (in Chinese), 52(2):339-345. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxb200902005

     

    Li C, Van Der Hilst R D. 2010. Structure of the upper mantle and transition zone beneath Southeast Asia from traveltime tomography. Journal of Geophysical Research:Solid Earth, 115(B7):B07308, doi:10.1029/2009JB006882.

     

    Li Y, Yao H J, Liu Q Y, et al. 2010. Phase velocity array tomography of Rayleigh waves in western Sichuan from ambient seismic noise. Chinese Journal of Geophysics (in Chinese), 53(4):842-852, doi:10.3969/j.issn.0001-5733.2010.04.00.

     

    Li Z C. 2014. Research status and development trends for seismic migration technology. Oil Geophysical Prospecting (in Chinese), 49(1):1-21. http://en.cnki.com.cn/Article_en/CJFDTOTAL-SYDQ201401004.htm

     

    Ligorría J P, Ammon C J. 1999. Iterative deconvolution and receiver-function estimation. Bulletin of the Seismological Society of America, 89(5):1395-1400. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ024970340/

     

    Liu Q Y, Kind R, Li S C. 1996. Maximal likelihood estimation and nonlinear inversion of the complex receiver function spectrum ratio. Chinese Journal of Geophysics (Acta Geophysica Sinica) (in Chinese), 39(4):500-511. http://en.cnki.com.cn/Article_en/CJFDTOTAL-DQWX604.009.htm

     

    Liu Q Y, Kind R. 2004. Multi-channel maximal likelihood deconvolution method for isolating 3-component teleseismic receiver function. Seismology and Geology (in Chinese), 26(3):416-425. http://d.old.wanfangdata.com.cn/Periodical/dzdz200403006

     

    Liu Q Y, Van Der Hilst R D, Li Y, et al. 2014. Eastward expansion of the Tibetan plateau by crustal flow and strain partitioning across faults. Nature Geoscience, 7(5):361-365, doi:10.1038/NGEO2130.

     

    Liu X W, Liu H. 2002. Status and progress on wave equation migration methods. Progress in Geophysics (in Chinese), 17(4):582-591. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxjz200204004

     

    Liu Z, Tian X B, Gao R, et al. 2017. New images of the crustal structure beneath eastern Tibet from a high-density seismic array. Earth and Planetary Science Letters, 480:33-41. doi: 10.1016/j.epsl.2017.09.048

     

    Mancinelli N J, Fischer K M. 2017. The spatial sensitivity of Sp converted waves-scattered-wave kernels and their applications to receiver-function migration and inversion. Geophysical Journal International, 212(3):1722-1735. https://www.researchgate.net/publication/323175837_The_spatial_sensitivity_of_Sp_converted_waves-scattered-wave_kernels_and_their_applications_to_receiver-function_migration_and_inversion

     

    Owens T J, Zandt G, Taylor S R. 1984. Seismic evidence for an ancient rift beneath the Cumberland plateau, Tennessee:A detailed analysis of broadband teleseismic P waveforms. Journal of Geophysical Research:Solid Earth, 89(B9):7783-7795. doi: 10.1029/JB089iB09p07783

     

    Qi S H, Liu Q Y, Chen J H, et al. 2016. Attenuation of noise in receiver functions using curvelet transform. Chinese Journal of Geophysics (in Chinese), 59(3):884-896, doi:10.6038/cjg20160311.

     

    Robert, Zhu J, Vergne J, et al. 2009. Crustal structures in the area of the 2008 Sichuan earthquake from seismologic and gravimetric data. Tectonophysics, 491(1-4):205-210. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=2f02dde1c8510048f4351495934286e7

     

    Ryberg T, Weber M. 2000. Receiver function arrays:A reflection seismic approach. Geophysical Journal International, 141(1):1-11. doi: 10.1046/j.1365-246X.2000.00077.x

     

    Sa L M, Yang W Y, Du Q Z, et al. 2015. Past, present and future of seismic migration imaging. Oil Geophysical Prospecting (in Chinese), 50(5):1016-1036. http://www.en.cnki.com.cn/Article_en/CJFDTOTAL-SYDQ201505034.htm

     

    Shang X F, De Hoop M V, Van Der Hilst R D. 2012. Beyond receiver functions:Passive source reverse time migration and inverse scattering of converted waves. Geophysical Research Letters, 39(15):L15308, doi:10.1029/2012GL052289.

     

    Stehly L, Froment B, Campillo M, et al. 2015. Monitoring seismic wave velocity changes associated with the MW7.9 Wenchuan earthquake:increasing the temporal resolution using curvelet filters. Geophysical Journal International, 201(3):1939-1949. doi: 10.1093/gji/ggv110

     

    Vinnik L P, Reigber C, Aleshin I M, et al. 2004. Receiver function tomography of the central Tien Shan. Earth and Planetary Science Letters, 225(1-2):131-146. doi: 10.1016/j.epsl.2004.05.039

     

    Wilson D, Aster R. 2005. Seismic imaging of the crust and upper mantle using regularized joint receiver functions, frequency-wave number filtering, and multimode Kirchhoff migration. Journal of Geophysical Research:Solid Earth, 110(B5):B05305, doi:10.1029/2004JB003430.

     

    Wu J P, Huang Y, Zhang T Z, et al. 2009. Aftershock distribution of the MS8.0 Wenchuan earthquake and three dimensional P-wave velocity structure in and around source region. Chinese Journal of Geophysics (in Chinese), 52(2):320-328. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxb200902003

     

    Wu Q J, Tian X B, Zhang N L, et al. 2003a. Receiver function estimated by maximum entropy deconvolution. Acta Seismologica Sinica (in Chinese), 25(4):382-389. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dizhen-e200304005

     

    Wu Q J, Tian X B, Zhang N L, et al. 2003b. Receiver function estimated by wiener filtering. Earthquake Research in China (in Chinese), 19(1):41-47. http://www.cnki.com.cn/Article/CJFDTOTAL-ZDZW200304008.htm

     

    Wu Q J, Li Y H, Zhang R Q, et al. 2007. 2D Kirchhoff migration for receiver function. Chinese Journal of Geophysics (in Chinese), 50(2):539-545. doi: 10.1002/cjg2.v50.2

     

    Yuan X H, Ni J, Kind R, et al. 1997. Lithospheric and upper mantle structure of southern Tibet from a seismological passive source experiment. Journal of Geophysical Research:Solid Earth, 102(B12):27491-27500. doi: 10.1029/97JB02379

     

    Yue H, Chen Y J, Sandvol E, et al. 2012. Lithospheric and upper mantle structure of the northeastern Tibetan plateau. Journal of Geophysical Research:Solid Earth, 117(B5):B05307, doi:10.1029/2011JB008545.

     

    Zhang Z J, Wang Y H, Chen Y, et al. 2009. Crustal structure across Longmenshan fault belt from passive source seismic profiling. Geophysical Research Letters, 36(17):L17310, doi:10.1029/2009GL039580.

     

    Zhang Z J, Yuan X H, Chen Y, et al. 2010. Seismic signature of the collision between the east Tibetan escape flow and the Sichuan Basin. Earth and Planetary Science Letters, 292(3-4):254-264, doi:10.1016/j.epsl.2010.01.046.

     

    Zheng H S. 2000. Study development of seismic migration. World Geology (in Chinese), 19(4):392-396, 401. http://d.old.wanfangdata.com.cn/Periodical/wtyht201803018

     

    Zhu L P, Kanamori H. 2000. Moho depth variation in southern California from teleseismic receiver functions. Journal of Geophysical Research:Solid Earth, 105(2):2969-2980. https://www.researchgate.net/publication/215755387_Moho_depth_variation_in_southern_California_from_teleseismic_RFs

     

    Zhu L P. 2000. Crustal structure across the San Andreas fault, southern California from teleseismic converted waves. Earth and Planetary Science Letters, 179(1):183-190. doi: 10.1016/S0012-821X(00)00101-1

     

    郭飚, 刘启元, 陈九辉等. 2009.川西龙门山及邻区地壳上地幔远震P波层析成像.地球物理学报, 52(2):346-355. http://www.geophy.cn//CN/abstract/abstract914.shtml

     

    雷建设, 赵大鹏, 苏金蓉等. 2009.龙门山断裂带地壳精细结构与汶川地震发震机理.地球物理学报, 52(2):339-345. http://www.geophy.cn//CN/abstract/abstract913.shtml

     

    李昱, 姚华建, 刘启元等. 2010.川西地区台阵环境噪声瑞利波相速度层析成像.地球物理学报, 53(4):842-852, doi:10.3969/j.issn.0001-5733.2010.04.009. http://www.geophy.cn//CN/abstract/abstract3003.shtml

     

    李振春. 2014.地震偏移成像技术研究现状与发展趋势.石油地球物理勘探, 49(1):1-21. http://d.old.wanfangdata.com.cn/Periodical/kxycf201408085

     

    刘启元, Kind R, 李顺成. 1996.接收函数复谱比的最大或然性估计及非线性反演.地球物理学报, 39(4):500-511. doi: 10.3321/j.issn:0001-5733.1996.04.008 http://www.geophy.cn//CN/abstract/abstract4096.shtml

     

    刘启元, Kind R. 2004.分离三分量远震接收函数的多道最大或然性反褶积方法.地震地质, 26(3):416-425. doi: 10.3969/j.issn.0253-4967.2004.03.006

     

    刘喜武, 刘洪. 2002.波动方程地震偏移成像方法的现状与进展.地球物理学进展, 17(4):582-591. doi: 10.3969/j.issn.1004-2903.2002.04.004

     

    齐少华, 刘启元, 陈九辉等. 2016.接收函数的曲波变换去噪.地球物理学报, 59(3):884-896, doi:10.6038/cjg20160311. http://www.geophy.cn//CN/abstract/abstract12607.shtml

     

    撒利明, 杨午阳, 杜启振等. 2015.地震偏移成像技术回顾与展望.石油地球物理勘探, 50(5):1016-1036. http://d.old.wanfangdata.com.cn/Periodical/sydqwlkt201505027

     

    吴建平, 黄媛, 张天中等. 2009.汶川MS8.0级地震余震分布及周边区域P波三维速度结构研究.地球物理学报, 52(2):320-328. http://www.geophy.cn//CN/abstract/abstract900.shtml

     

    吴庆举, 田小波, 张乃铃等. 2003a.计算台站接收函数的最大熵谱反褶积方法.地震学报, 25(4):382-389. http://d.old.wanfangdata.com.cn/Periodical/dizhen200304005

     

    吴庆举, 田小波, 张乃铃等. 2003b.用Wiener滤波方法提取台站接收函数.中国地震, 9(1):41-47. http://d.old.wanfangdata.com.cn/Periodical/zgdz200301005

     

    吴庆举, 李永华, 张瑞青等. 2007.接收函数的克希霍夫2D偏移方法.地球物理学报, 50(2):539-545. doi: 10.3321/j.issn:0001-5733.2007.02.027 http://www.geophy.cn//CN/abstract/abstract1410.shtml

     

    郑海山. 2000.地震偏移技术研究进展.世界地质, 19(4):392-396, 401. doi: 10.3969/j.issn.1004-5589.2000.04.017

  • 加载中

(9)

计量
  • 文章访问数:  636
  • PDF下载数:  439
  • 施引文献:  0
出版历程
收稿日期:  2018-06-07
修回日期:  2018-11-29
上线日期:  2019-06-05

目录