Study on seismic anisotropy based on receiver function analysis at Yugan station
-
摘要: 首先基于理论地层模型,分析使用径向切向接收函数集联合反演地壳各向异性程序的可行性和可靠性。然后将该程序应用到余干地震台接收函数的实际分析中,结果显示,不同目标函数计算得到的快波偏振方向和时间延迟存在差异,联合目标函数结果表明,余干地震台快波方向约为79°,时间延迟约为0.19 s。Abstract: First of all, based on the theoretical crustal model, the feasibility and reliability of using the radial and transverse receiver functions to jointly invert the crustal anisotropy program are analyzed. Then the program is applied to the actual analysis of the receiver functions of Yugan seismic station. The results show that the polarization direction and time delay of the fast wave calculated by different objective functions are different. The joint objective function shows that the fast wave direction of Yugan seismic station is about 79° and the time delay is about 0.19 s.
-
Key words:
- receiver function /
- radial receiver function /
- transverse receiver function /
- anisotropy
-
表 1 计算接收函数的地层模型
Table 1. Models used for receiver function modeling
地层厚度/m 密度 g/cm3 vP/m·s−1 vS/m·s−1 F %P %S Tr Pl Str Dip A组各向异性 20 000 2 720 6 100 3 500 1 0 0 0 0 0 0 A1 20 000 2 720 6 100 3 500 0 6 6 0 0 0 0 0 3 300 8 200 4 600 1 0 0 0 0 0 0 20 000 2 720 6 100 3 500 1 0 0 0 0 0 0 A2 20 000 2 720 6 100 3 500 0 6 6 45 0 0 0 0 3 300 8 200 4 600 1 0 0 0 0 0 0 注:F取值0、1,0表示该层为各向异性层,1表示该层为各向同性层,%P和%S分别为P波和S波的各向异性程度,Tr为快轴方向的水平投影与正北方向的夹角,Pl为快轴方向与铅锤方向的夹角,Str为倾斜界面的走向(倾向为走向加90°),Dip为倾斜界面的倾角 -
[1] Silver P G,Chan W W. Shear wave splitting and subcontinental mantle deformation[J]. J. Geophys. Res.,1991,96(B10):16429-16454 doi: 10.1029/91JB00899 [2] McNamara D E,Owens T J. Azimuthal shear wave velocityanisotropy in the Basin and Range Province using Moho Ps convertedphases[J]. J. Geophys. Res.,1993,98(B7):12003-12017 doi: 10.1029/93JB00711 [3] Nagaya M,Oda H,Akazawa H,et al. Receiver functions of seismic waves in layered anisotropic media:Application to the estimate of seismic anisotropy[J]. Bull. Seism. Soc. Amer.,2008,98(6):2990-3006 doi: 10.1785/0120080130 [4] 刘莎. 利用剪切波分裂方法研究玉树和冰岛地区各向异性[D]. 北京: 中国地震局地球物理研究所, 2013 [5] Liu H F,Niu F L. Estimating crustal seismic anisotropy with a joint analysis of radial and transverse receiver function data[J]. Geophys. J. Int.,2012,188(1):144-164 doi: 10.1111/j.1365-246X.2011.05249.x [6] Owens T J,Zandt G,Taylor S R. Seismic evidence for an ancient rift beneath the Cumberland plateau,Tennessee:A detailed analysis of broadband teleseismic P waveforms[J]. J. Geophys. Res.,1984,89(B9):7783-7795 doi: 10.1029/JB089iB09p07783 [7] Sun Y,Niu F,Liu H,et al. Crustal structure and deformation of the SE Tibetan plateau revealed by receiver function data[J]. Earth. Planet. Sci. Lett.,2012,349-350:186-197 doi: 10.1016/j.jpgl.2012.07.007 [8] 杨妍,姚华建,张萍,等. 用接收函数方法研究华北克拉通中部造山带及其邻域地壳方位各向异性[J]. 中国科学:地球科学,2018,48(7):912-923 [9] 刘云昌,孙娅,张俊. 利用远震接收函数研究安徽及周边地区地壳结构和各向异性[J]. 地球物理学进展,2019,34(2):442-453 doi: 10.6038/pg2019CC0054 [10] Frederiksen A W,Bostock M G. Modeling teleseismic waves in dipping anisotropic structure[J]. Geophys. J. Int,2000,141(2):401-412 doi: 10.1046/j.1365-246x.2000.00090.x [11] Herrmann R B, Ammon C J. Computer programs in seismology[D]. Department of Earth and Atmospheric Sciences, St Louis University, St Louis, Missouri, 2004 [12] Ligorria J P,Ammon C J. Iterative deconvolution and receiver-function estimation[J]. Bull. Seismol. Soc. Amer.,1999,89(5):1395-1400 [13] Andreas Wustefeld. Methods and applications of shear wave splitting: The east European craton[D]. France: Universite Montpellier, 2007 [14] 徐震,徐鸣洁,王良书,等. 用接收函数Ps转换波研究地壳各向异性−以哀牢山—红河断裂带为例[J]. 地球物理学报,2006,49(2):438-448 doi: 10.3321/j.issn:0001-5733.2006.02.017 -