云南禄劝地震部分余震的矩张量反演
MOMENT TENSOR INVERSION OF SOME AFTERSHOCKS OF THE APRIL 18, 1985, LUQUAN, YUNNAN, CHINA, EARTHQUAKE
-
摘要: 利用数字盒式磁带记录加速度仪组成的流动地震台网所记录的三分向近场加速度图,通过矩张量反演确定了1985年4月18日云南禄劝 Ms=6.1地震的部分余震的震源机制.以均匀半无限弹性介质的格林函数解释路径效应,并通过正演计算识别由两次积分得到的位移地震图中的直达 P 波、直达 S 波和 SP 转换波震相,然后用这些震相进行矩张量反演.反演结果表明,在解超定线性方程组时,采用适当的加权系数,可使上述直达波和转换波的理论计算值与观测值拟合得较好.尽管用以反演的三个余震大小不同(震级 ML48,3.2,3.5),震源位置也有差别,但它们的震源机制却非常接近,且与主震的震源机制相当一致.这一特征显示了余震的发生与主震发震构造的内在联系.这些实例说明,由震源球球面上分布适宜的数字地震台网取得的近场加速度资料,借助于即使是简单的介质模型,通过地震矩张量反演,不但可以得出这些地震震源的主要成份————剪切位错源,同时还可得出震源所含的其它成分,如膨胀源和补偿线性向量偶极.Abstract: Based on the three-component accelerograms, recorded at near-field distance by a temporary seismic network consisting of digital cassette tape reacording accelerographs, the focal mechanisms of three aftershocks of the April 18, 1985, Luquan, Yunnan Province, China, earthquake of magnitude Ms=6.1, are calculated using moment tensor inversion technique. The phases of direct P, S and converted SP waves in the displacement seismograms, produced by twice integration of the observed accelerograms, are identified via forward calculation using Green's functions for homogeneous semi-infinite elastic medium, and used in the inversion. The results of inversion show that a better fit of synthetic to the observed seismograms of direct as well as converted phases can be achieved if appropriate weighting functions are used in solving the over-definite linear equations. While these aftershocks are of different magnitudes (ML=4.8, 3.2 and 3.5, respectively) and hypocentral locations, their focal mechanisms are very similar and consistent with that of the main shock. This feature demonstrates the intrinsic correlation between the occurrence of aftershocks and the seismogeneic fault of main shock. Our experimentations show that using the near-field accelerogram obtained from the digital seismic network with appropriate azimuthal coverage on the focal sphere, with the aid of even simple earth model, not only the shear dislocation source, but also the isotropic part and CLVD (compensated linear vector dipole) can be retrieved by the technique of moment tensor inversion.
-
引 言
地震波(P，S，pP，sP等)在传播过程中，遇到间断面会发生反射或转换，所产生的新震相(P410P，P660P，S410S，S660S，S410P，S660P，s410P等)为间断面(“410”，“660”等)的深入研究提供了地震学基础(Flanagan，Shearer，1998a； Deuss，2009； Schmerr，Thomas，2011)． 这些前驱震相在单条记录中的能量(幅度)往往很弱，常淹没于噪声之中； 但是通过叠加大量的观测波形数据，则可以压制噪声，更加有效地提取与间断面相关的有用信号(Rost，Thomas，2002； 臧绍先，周元泽，2002)．
很多研究者利用下行的转换震相在间断面的转换点深度确定该间断面的深度． 例如： Richards和Wicks(1990)使用SdP转换波，利用非线性叠加方法对汤加地区下方的“670”间断面的深度和性质进行了深入研究； 谢彩霞等(2012)利用4次根倾斜叠加方法有效提取了离源下行的SdP次生转换震相，进一步确认了汤加—斐济下方300 km附近速度间断面的存在； Li等(2008)通过叠加S--P转换波，计算了中国东北地区“660”间断面的深度，并探讨了俯冲带对“660”间断面的影响．
发生在俯冲带中的深震会在近源区速度间断面产生反射震相(p410P，s410P等)和转换震相(S410P，S660P等)(Flanagan，Shearer，1998b)． 长周期震相(PP，SS)前驱波方法对狭窄俯冲带间断面的横向分辨具有很大局限性(Flanagan，Shearer，1998a； Schmerr，Thomas，2011)． 接收函数方法是研究地球内部间断面的一种有效方法，可以通过P波与其上行转换波Pds的到时差来确定间断面的深度(Langston，1979)，但其主要适用于三分量地震台站下方壳、 幔间断面的研究．
在深源地震(震源深度h＞300 km)分布密集的俯冲带(如汤加—斐济俯冲带)，使用深震震相sP在近源区速度间断面底界面反射的震相s410P确定狭窄俯冲带间断面的横向变化具有很大的优势． 本文将通过倾斜叠加大量的观测波形数据，有效提取sP在近源区“410”速度间断面底界面反射的弱前驱波s410P，利用sP震相与其前驱波s410P震相的到时差计算汤加—斐济俯冲带“410”速度间断面的深度．
1. 台网资料
甘肃地震台网由甘肃数字测震台网和甘东南野外观测流动台阵构成，二者均处于青藏高原东北缘． 甘肃“十五”数字测震台网自2008年6月正式运行，由兰州、 高台、 安西、 嘉峪关和天水等5个有人值守的国家数字地震台和39个区域遥测数字地震台组成(冯建刚等，2012)，目前拥有包括邻省(宁夏、 陕西、 四川、 青海、 内蒙)在内的76个数字测震台站． 甘东南野外观测流动台阵是中国地震局地质研究所和中国地震局兰州地震研究所共同合作在甘东南地区勘址架设，共设有7条测线，150个流动台站，平均台间距约为10 km(图 1)． 所有台站均使用REF TEK-130型数据采集器和Guralp公司生产的CMG-3ESPC宽频带地震计． 甘东南野外观测流动台阵运行期间(2009年11月—2011年12月)积累了连续、 可靠、 高质量的地震观测波形数据，为深入研究地球内部物理结构提供了宝贵资料(秦满忠等，2015)．
![]()
图 1 甘肃数字测震台网和甘东南野外观测流动台阵分布Figure 1. Distribution of the Gansu Digital Seismic Network (blue triangles) and southeastern Gansu temporary observation array (black triangles)本文选取了甘肃数字测震台网和甘东南野外观测流动台阵共同记录的2011年9月15日发生在汤加—斐济俯冲带(21.61°S，179.53°W，h=644 km，mB=7.3)的深远地震数字观测波形数据，定位结果参考了USGS(2011)地震目录．
2. 数据处理
深震震相sP与其前驱波s410P在地幔中的传播路径基本相同． sP是上行的S波在震中附近地表反射转换的P波； s410P是上行的S波在震中附近“410”间断面底界面反射转换的P波(图 2)，该震相在初至震相P与sP之间到达(图 3)．
![]()
图 3 甘肃地震台网记录到的2011年9月 15日发生在汤加—斐济俯冲带的深远 地震的数字观测波形sMohoP， sLABP和s410P震相分别表示sP在 近源区莫霍面， 岩石圈和“410”的前驱波 The seismic phases sMohoP， sLABP and s410P are the precursors from near-source underside reflection off Moho， lithosphere-asthenosphere boundary， and “410”, respectivelyFigure 3. Waveforms of mB=7.3 event occurred in Tonga-Fiji on September 15， 2011， recorded by the Gansu Seismic Network首先对所选取的汤加—斐济深远地震观测波形数据作重采样、 去倾斜、 去均值处理; 然后作带通滤波处理(0.1—0.01 Hz)，同时去掉记录畸形、 信噪比差的观测波形； 最终选取了甘肃数字测震台网和甘东南野外观测流动台阵共同记录到的170条垂直分量观测波形数据． 如图 3所示，汤加—斐济深远震在台网记录中的震中距为88°—102°．
通过倾斜叠加大量的观测波形数据，可以压制噪声，有效地提取与间断面相关的弱信号(Ritsema et al，1995； 臧绍先，周元泽，2002)． 在进行倾斜叠加的过程中，选取震中距为94°的台站(红崖山台)作为参考台，sP为参考震相，并以sP震相到时为零时刻(图 3)． 信号i(s410P)在任意台s的到时为Δtir+ΔpΔDrs，以红崖山台为参考台的倾斜叠加结果为
式中，Δtir为参考震中距记录的信号i(s410P)与参考震相sP之间的走时差，Δp为信号i在任意台s的水平慢度与信号i(震中距为参考台震中距)水平慢度之差，ΔDrs为任意台s与参考台的震中距之差， m为使用地震记录的台站数．
以sP为参考震相，使用倾斜叠加方法获得了在震中附近“410”间断面底界面反射的前驱波信号s410P，其反射点位置为(21.25°S，179.99°W)． 该信号在倾斜叠加波形(图 3中红色线条)中成像清晰，与参考震相sP的观测走时差为137.43 s，使用TauP软件结合理论IASP91全球速度模型得到的s410P与震相sP理论走时差为137.04 s； 通过使用基于CRUST1.0模型(Laske et al，2013)校正后的汤加—斐济地区速度模型IASP91_Tonga(图 4)，反算出“410”间断面底部反射点的深度约为398.5 km，“410”间断面抬升约11.5 km，这与俯冲带与近源区“410”间断面的相互作用有关．
![]()
图 4 全球速度模型IASP91以及修改后的汤加—斐济地区速度模型IASP91_TongaFigure 4. IASP91 velocity model and modified velocity model (IASP91_Tonga) of the Tonga-Fiji zone in this paper3. 讨论与结论
甘肃数字测震台网和甘东南野外观测流动台阵记录的汤加—斐济深震震相pP和sP均具有尖锐、 清晰和较大振幅的记录特征，sP记录振幅甚至超过了直达P波，这为我们使用前驱波(sMohoP，sLABP和s410P等)对近源区间断面(莫霍面，岩石圈，“410”等)的深入研究提供了基础资料．
鉴于前驱波震相在单条记录中的能量(振幅)往往很弱，本文通过倾斜叠加大量的观测波形数据，获得了在近源区“410”间断面底界面反射的前驱波震相s410P，并使用CRUST1.0模型校正后的汤加—斐济地区速度模型IASP91_Tonga，反算出“410”间断面底部反射点的深度约为398.5 km．
图 3中震相sLABP和sMohoP相对于参考震相sP的走时差分别为23.16 s和10.62 s. 我们同样对震源区的地壳及岩石圈结构作了修改(图 4)，使用TauP软件计算得到近源区岩石圈的深度约为46 km，莫霍面深度约为10 km．
俯冲带对“410”和“660”的影响是目前研究的一个热点(Li et al，2008)． 该研究涉及到间断面的性质，又能反映地幔对流的性质和形式． 俯冲带是冷的下沉物质，其通过间断面时会引起间断面的形态变化． 若“410”为相变界面，俯冲带会使它抬升(Vidale，Benz，1992； Collier，Helffrich，1997); 若“660”为相变界面，俯冲板块中的冷物质将使它下沉． 因此，汤加—斐济俯冲带中冷的物质与“410”作用时，俯冲带及相邻地幔中的橄榄石-尖晶石的相变界面会上升(蒋志勇等，2003)． 本文中得到的汤加—斐济俯冲带“410”间断面抬升约11.5 km与该结论一致．
本文使用倾斜叠加技术成功提取了sP前驱波(s410P，sLABP和sMohoP)，获得了俯冲带间断面的深度，给出了汤加—斐济俯冲带新的地震波速度模型(IASP91_Tonga)． 本文结果将为深入研究该区域的地球结构提供基础资料，对正确认识俯冲带的复杂结构以及深入理解地球深部的动力学过程具有重要的参考价值．
-
[1] Stump, B. W. and Johnson, L. R. 1977. The determination of source properties by the linear inversion of seismograms. Bull. Seism. Soc. Anger., 67. 1489——1502.
[2] Gilbert, F, and Dziewonski, A. M——1975. An application of normal mode theory to the retrieval of structural parameters and source mechanisms from seismic spectra. Phil. Truns. R. Soc,278, 187——269
[3] McCowan, D. W. 1976. Moment tensor representation of surface wave sources. Geophys. J. R. A.str. Su 44. 595——599.
[4] Mendigurcn, J. A.,1977. lnvcrsion of surface wave data in source mechanism studies. J. Geophys. Res.,82. 889——894.
[5] Aki. K, and Patton. H.,1978. Determination of seismic moment tensor using surface waves. Tectnnophysics.49. 213——222.
[6] Kanamori. H. and Given, J. W, 1981. Use of longperiod surface waves for rapid determination of earthquak source parameters, Piys. Furth Planet. Interiors,27. 8——31.
[7] Romannowicz, B——1982. Moment tensor inversion of long period Rayleigh waves:A new approach. J. Geophys. Res.,87. 5395——5407.
[8] Lay. T,Given. J. W. and Kanamori, H., 1982. Long——period mechanism of the 8 November 1980 Eureka, California, earthquake.Bull. Seism. S'oc. Anger., 72. 439——456.
[9] Fitch, T. J.,.McCowan, D. W. and Shields, M. W.,1980. Estimation of the seismic moment tensor from teleseismic body wave data with applications to intraplate and mantle earthquakes. J. Geophys. Res., 15,3817——3828
[10] Longston, C. A., 1981. Source inversion of seismic waveforms: The Koyna. India, earthquakes of 13 September 1967 Bull. Seism. Soc. Amer.,71, 1——24.
[11] Doornbos. D. J.,1982. Seismic moment tensors and kinematic source parameters. Geophy.s. J. R. Astr. Sor.,69. 235——251.
[12] Barker, J. S. and Langston, C. A., 1982. Moment tensor inversion of complex earthquakes·Geophys. J. R.9.soc.,68. 777——803.
[13] Stump. B. W. and Johnson. L. R,1984. Nearfield source characterization ofcontained nuclear explosion in tuff Bull. Seisrn. Soc. Arur. 74, 1——26.
[14] Lanczos, C. 1961. Linear Differential Operators. Van Nostrand, London, 564.
[15] 吴明熙、王鸣、孙次昌、柯兆明、王培德、陈运泰、吴大铭,1990. 1985年地震部分余震的精确定位.地球物理学报,12, 2. 121——129.
[16] Lamb, H——1904. On the propagation of tremors over the surface of an elastic solid. Phil. Truns. R. Sm. London. A203,1——42.
[17] Johnson. J. R.,1974. Green's function for Lamb's problem. Geopftl}.s. J. R.4.str. Soc.,37, 91——131.
[18] Dix. C. H——1954. The method of Cagniard in seismic pulse problem. Geophrsic.s, 19, 722——738.
[19] Cagniard. L., 1962. ReJlectiort and Refraction of Progreesive Seismic Waves. Translated and revised by Flinn,E. A. and Dix. C. H.,McGrawHilh New York.
[20] De Hoop. A. T.,1960. A modification of Cagniard's method for solving seismic pulse problems. Appl. Sci.Rcs. 88, 349——356.[1] Stump, B. W. and Johnson, L. R. 1977. The determination of source properties by the linear inversion of seismograms. Bull. Seism. Soc. Anger., 67. 1489——1502.
[2] Gilbert, F, and Dziewonski, A. M——1975. An application of normal mode theory to the retrieval of structural parameters and source mechanisms from seismic spectra. Phil. Truns. R. Soc,278, 187——269
[3] McCowan, D. W. 1976. Moment tensor representation of surface wave sources. Geophys. J. R. A.str. Su 44. 595——599.
[4] Mendigurcn, J. A.,1977. lnvcrsion of surface wave data in source mechanism studies. J. Geophys. Res.,82. 889——894.
[5] Aki. K, and Patton. H.,1978. Determination of seismic moment tensor using surface waves. Tectnnophysics.49. 213——222.
[6] Kanamori. H. and Given, J. W, 1981. Use of longperiod surface waves for rapid determination of earthquak source parameters, Piys. Furth Planet. Interiors,27. 8——31.
[7] Romannowicz, B——1982. Moment tensor inversion of long period Rayleigh waves:A new approach. J. Geophys. Res.,87. 5395——5407.
[8] Lay. T,Given. J. W. and Kanamori, H., 1982. Long——period mechanism of the 8 November 1980 Eureka, California, earthquake.Bull. Seism. S'oc. Anger., 72. 439——456.
[9] Fitch, T. J.,.McCowan, D. W. and Shields, M. W.,1980. Estimation of the seismic moment tensor from teleseismic body wave data with applications to intraplate and mantle earthquakes. J. Geophys. Res., 15,3817——3828
[10] Longston, C. A., 1981. Source inversion of seismic waveforms: The Koyna. India, earthquakes of 13 September 1967 Bull. Seism. Soc. Amer.,71, 1——24.
[11] Doornbos. D. J.,1982. Seismic moment tensors and kinematic source parameters. Geophy.s. J. R. Astr. Sor.,69. 235——251.
[12] Barker, J. S. and Langston, C. A., 1982. Moment tensor inversion of complex earthquakes·Geophys. J. R.9.soc.,68. 777——803.
[13] Stump. B. W. and Johnson. L. R,1984. Nearfield source characterization ofcontained nuclear explosion in tuff Bull. Seisrn. Soc. Arur. 74, 1——26.
[14] Lanczos, C. 1961. Linear Differential Operators. Van Nostrand, London, 564.
[15] 吴明熙、王鸣、孙次昌、柯兆明、王培德、陈运泰、吴大铭,1990. 1985年地震部分余震的精确定位.地球物理学报,12, 2. 121——129.
[16] Lamb, H——1904. On the propagation of tremors over the surface of an elastic solid. Phil. Truns. R. Sm. London. A203,1——42.
[17] Johnson. J. R.,1974. Green's function for Lamb's problem. Geopftl}.s. J. R.4.str. Soc.,37, 91——131.
[18] Dix. C. H——1954. The method of Cagniard in seismic pulse problem. Geophrsic.s, 19, 722——738.
[19] Cagniard. L., 1962. ReJlectiort and Refraction of Progreesive Seismic Waves. Translated and revised by Flinn,E. A. and Dix. C. H.,McGrawHilh New York.
[20] De Hoop. A. T.,1960. A modification of Cagniard's method for solving seismic pulse problems. Appl. Sci.Rcs. 88, 349——356. -
期刊类型引用(34)
1. 余怀忠，于晨，张小涛，苑争一，钟骏，解孟雨，李泽平，张靖雪. 中国CSEP检验中心的软件系统建设. 地震研究. 2025(02): 177-187 . 
百度学术
2. 韩颜颜，臧阳，马亚伟，解孟雨，王月，杨文，张小涛，宋金，姜祥华，史海霞，李飞，孟令媛. 2023年震情述评. 中国地震. 2024(01): 271-278 . 百度学术
3. 郭寅，李金，向元. 乌什M_S7.1地震序列及应力降特征研究. 内陆地震. 2024(02): 105-111 . 百度学术
4. 郭蕾，孟令媛. 河北及邻区地震活动固体潮调制统计分析. 中国地震. 2024(02): 399-409 . 百度学术
5. 赵文波，李强，黄艳丹. 2023年福建连城小震群活动分析. 科技创新与应用. 2024(28): 98-101 . 百度学术
6. 张思萌，李永生，张雁翔，石伟，刘长生，李俊秋. 2009年以来黑龙江省固体潮调制比特征与地震关系研究. 大地测量与地球动力学. 2024(10): 1053-1058 . 百度学术
7. 张思萌，张雁翔，李永生，刘双，李俊秋. 2018—2019年吉林省松原市宁江区2次5级以上地震序列测震学异常分析. 地震地磁观测与研究. 2024(06): 49-60 . 百度学术
8. 李继业，晏锐，张思萌，胡澜缤，孟令蕾，周晨. 井水位潮汐响应与小地震调制作用的关系. 地震地质. 2023(03): 668-688 . 百度学术
9. 钟骏，解孟雨，姜祥华，黎明晓，岳冲. 2021年3月30日西藏双湖M_S 5.8和6月16日青海茫崖M_S 5.8地震总结. 地震地磁观测与研究. 2023(01): 1-13 . 百度学术
10. 索朗南杰，桑杰卓玛，索朗多吉，次仁多吉，益西拉姆，高锦瑞. 2008年以来西藏地区强震前固体潮调制比特征分析. 地震地磁观测与研究. 2023(01): 14-22 . 百度学术
11. 李冬梅，李小晗，郑建常. 2022年泸定M6.8地震前地震活动参数变化研究. 中国地震. 2023(03): 609-621 . 百度学术
12. 李文君，曾宪伟，马翀之. 中国大陆6级以上强震及震前小震固体潮调制异常特征分析. 地震工程学报. 2022(03): 638-648 . 百度学术
13. 陈大庆，邵叶，刘锦，谭争光，李丽. 中国大陆强震月相调制性质一致性分析. 华南地震. 2022(03): 65-70 . 百度学术
14. 周依，郭蕾，王亚玲. 首都圈地区中强地震前固体潮调制比特征分析. 大地测量与地球动力学. 2022(11): 1166-1170 . 百度学术
15. 于晨，周志华，解孟雨，张小涛，韩颜颜，苑争一，黎明晓，姚丽，田雷. 2021年6月10日云南双柏M_S 5.1和6月12日云南盈江M_S 5.0地震总结. 地震地磁观测与研究. 2022(05): 151-168 . 百度学术
16. 张慧，谢小玲，阎春恒，王小娜，黄章荣. 华南沿海地震带地震调制比特征分析. 华南地震. 2022(04): 109-116 . 百度学术
17. 高锦瑞，索朗占堆，姜祥华，解孟雨，黎明晓，赵静. 2021年11月30日西藏双湖M_S 5.8地震总结. 地震地磁观测与研究. 2022(05): 169-178 . 百度学术
18. 陈婷，张新东，王晓山，章阳. 2022年芦山M_S 6.1和马尔康M_S 6.0地震前地震调制比异常分析. 地震地磁观测与研究. 2022(06): 1-7 . 百度学术
19. 邵媛媛，杨士超，王岩，钱蕊，毛佳宁. 辽宁地区固体潮调制比空间异常特征震例“回溯性”研究. 防灾减灾学报. 2021(02): 1-8 . 百度学术
20. 姜祥华，解孟雨，马未宇，宋治平，屠泓为，刘文邦，高锦瑞. 2020年7月23日西藏尼玛M_S 6.6地震总结. 地震地磁观测与研究. 2021(02): 32-41 . 百度学术
21. 苑争一，宋治平，姜祥华，马亚伟，姚琪，解孟雨，于晨，韩颜颜. 2020年3月23日新疆拜城5.0级和7月13日霍城5.0级地震总结. 地震地磁观测与研究. 2021(02): 17-31 . 百度学术
22. 田雷，张小涛，解孟雨，姜祥华，王月，邓世广，臧阳，黎明晓，于晨，苑争一. 2021年3月19日西藏比如M_S 6.1地震总结. 地震地磁观测与研究. 2021(02): 42-57 . 百度学术
23. 高荣，聂晓红，刘建明，向元，张琳琳. 2021年3月24日新疆拜城M_S5.4地震震前异常分析. 内陆地震. 2021(02): 109-119 . 百度学术
24. 解孟雨，马亚伟，张小涛，于晨，苑争一，姜祥华，卢显，王月，韩颜颜. 2021年3月24日新疆拜城M_(S)5.4地震总结. 地震地磁观测与研究. 2021(03): 250-263 . 百度学术
25. 林眉，郑建常. 长清M4.1地震前区域小震调制比异常分析. 大地测量与地球动力学. 2021(11): 1207-1210 . 百度学术
26. 孙楠，付虹，彭关灵，刘自凤. 云南地区强震前小震调制比特征分析. 地震地磁观测与研究. 2021(05): 82-89 . 百度学术
27. 张小涛，苑争一，姚丽，田雷，薛艳，杨文，姜祥华，史海霞，解孟雨，邓世广. 2021年5月21日云南漾濞M_S6.4地震总结. 地震地磁观测与研究. 2021(06): 165-180 . 百度学术
28. 孟令媛，苑争一，宋治平，杨文，姜祥华，王博，周志华，于晨，姚丽，解滔，黎明晓，姚琪，马未宇，马亚伟，解孟雨，臧阳. 2020年1月19日新疆伽师M_S 6.4地震总结. 地震地磁观测与研究. 2020(02): 63-89 . 百度学术
29. 杨文，宋治平，苑争一，张小涛，马亚伟，于晨，姚琪，姜祥华，姚丽. 2020年1月16日新疆库车5.6级地震总结. 地震地磁观测与研究. 2020(03): 147-161 . 百度学术
30. 臧阳，薛艳，姚琪，解孟雨，苑争一，王阅兵，马亚伟，张小涛. 2020年1月25日西藏丁青5.1级地震总结. 地震地磁观测与研究. 2020(03): 162-173 . 百度学术
31. 刘萍，聂晓红，高荣. 2020年6月26日新疆于田M_S6.4地震序列特征及震前地震活动异常分析. 内陆地震. 2020(03): 239-248 . 百度学术
32. 马亚伟，宋治平，杨文，解孟雨，苑争一，钟骏，李美，姚丽，邓世广，臧阳. 2020年5月6日新疆乌恰5.0级和5月9日柯坪5.2级地震总结. 地震地磁观测与研究. 2020(04): 179-192 . 百度学术
33. 韩颜颜，马亚伟，解孟雨，臧阳，苑争一，于晨，田雷，王月，张小涛，姚琪，薛艳，韦进. 2020年2月3日四川青白江5.1级地震总结. 地震地磁观测与研究. 2020(05): 141-152 . 百度学术
34. 李莎，阎春恒，解朝娣，袁媛. 2019年靖西M_S5.2地震序列潮汐触发分析. 华北地震科学. 2020(S2): 68-73 . 百度学术
其他类型引用(0)
下载:
计量
- 文章访问数: 1538
- HTML全文浏览量: 17
- PDF下载量: 141
- 被引次数: 34
