A Preliminary Study on the Attenuation Characteristics of High Frequency Ground Motionin the Western Margin of Ordos
-
摘要: 以鄂尔多斯西缘地区(34°N—42°N,103°E—109°E)为研究对象,基于2007—2020年研究区域内14个强震动台站记录到的116条强震记录(1≤MS≤8),采用加速度反应谱法,通过MATLAB软件编程拟合鄂尔多斯西缘地区kappa值与震中距的线性回归关系。研究结果表明,鄂尔多斯西缘地区强震台站高频衰减参数kappa0值为0.01545~0.06560 s;等效剪切波速(VS20、VS30)与震中距存在对数关系;随着VS20、VS30的增大,kappa0值逐渐减至0。Abstract: The western margin of Ordos (N: 34°—42°; E: 103 °—109 °) as the research object, using 2007 to 2020, the research area of 14 strong motion stations recorded of the 116 strong earthquake events data(2≤MS≤8), using Anderson and Hough the classic algorithm, using MATLAB software programming to calculate the ordos west region, the linear relationship between the kappa value and the epicenter distance to calculate the ordos kappa0 earthquake stations in the western margin high frequency attenuation parameter value is between 0.01545 s to 0.06560 s. There is a logarithmic relationship between equivalent shear wave velocities (VS20, VS30) and epicenter distance.With the increase of VS20 and VS30, the kappa0 value gradually decreases to 0.
-
Key words:
- Western margin of Erdos region /
- Ground motion /
- High Frequency attenuation /
- Kappa /
- Epicentraldistance
-
引言
汤阴地堑西靠太行山隆起,东依内黄隆起,是太行山前断裂带的重要组成部分。大量学者研究认为,汤阴地堑南部是断陷型地堑,汤东、汤西断裂是其东西边界断裂;北部是东深西浅的箕型半地堑,具有清晰的上、中、下地壳结构特征(韩慕康等,1980;杨承先,1984)。深地震测深资料表明,汤阴地堑壳幔结构复杂,下地壳存在低速层,莫霍面向上突起,复杂的深部构造环境使得中强地震的发生成为可能(任青芳等,1998)。深地震反射资料揭示了汤阴地堑地壳与莫霍面的构造特征及断裂的深、浅交切关系(刘保金等,2012)。布格重力资料给出了汤东、汤西断裂向下延展情况与空间展布特征,认为汤东、汤西断裂等太行山山前断裂是分割地块、深切下地壳、规模宏大的断裂带(徐志萍等,2015;姜磊等,2018)。煤田、石油资料的研究表明,汤东、汤西断裂处在地震活动期,现今仍在活动,断裂所在区域为地震活动带(于慎谔等,2012)。这些研究成果为研究汤阴地堑壳幔结构、断裂深部构造特征及地震危险性评价提供了重要的基础资料。
地震反射波法是探测地壳内部结构和构造的有效技术手段,该方法通过采用多次覆盖同相叠加技术,突出有效波、削弱干扰波,提高地震资料的信噪比,所得到的地震剖面能够直观形象地反映地下结构与构造,该方法已在石油勘探、煤田勘探、壳幔结构探测以及活断层探测研究中得到广泛应用。汤阴地堑及其邻区深部构造的研究成果较多,但对地堑地壳浅部结构及断裂浅部构造特征研究较少,特别是汤阴地堑南、北部地壳结构差异性研究缺乏可靠的地震资料。本文通过对安阳市与新乡市活断层探测深地震反射剖面数据重新解译,得到反映汤阴地堑北部与南部地下结构和断裂特征的高分辨率中深层地震反射剖面图像。利用获得的中深层地震反射剖面,结合研究区已有的浅层地震反射剖面,通过中、浅地震反射剖面相结合的方法,对汤阴地堑地壳浅部结构、隐伏断裂特征进行分析研究。研究结果不仅为汤阴地堑上地壳结构研究、断裂特征分析提供了基础资料,同时为地堑周边城市建设如何有效避让断裂带提供了科学依据。
1. 研究区地质背景及主要断裂概况
汤阴地堑是太行山前断裂带的重要组成部分,整体呈北北东走向,南北被新乡-封丘断裂与安阳断裂所围限,东西被汤东、汤西断裂所控制。多期下沉是其新生代的主要构造活动,地堑内覆盖着较厚的新生代沉积层,其东界汤东断裂的活动性要强于西界汤西断裂。(梁生正等,2003;郑建彬,2015)。
汤东断裂是汤阴地堑的主要控制性断裂,影响并控制了汤阴地堑的形成、演化与发展。断裂为走向北东、倾向北西的正断层,长约90km。断裂以东为内黄隆起,其基底由太古界与古生界岩层组成,地表覆盖新近纪与第四纪地层;断裂以西为汤阴地堑,其内沉积了较厚的新生代地层(花鑫升等,2018)。
汤西断裂是汤阴地堑的西界断裂,也是太行山与南华北平原的主要分界线,为走向北东、倾向南东的正断层,断裂南起新乡市西侧,向北延伸至汤阴宜沟附近。该断裂南段古近纪、新近纪时期有过强烈的构造活动,上升盘一侧的太行山隆起区局部出露太古界、古生界及少量的古近纪与新近纪地层,而下降盘一侧的汤阴地堑则覆盖着厚度千米以上的新生代沉积层。
2. 测线位置和方法
中深层地震反射剖面ML1位于汤阴地堑北部汤阴县附近,探测方向为自东向西,剖面起点位于汤阴县菜园镇小坡村村北(35°57′39.8″N,114°33′47.6″E),终点位于鹤壁市鹤壁集前柳江村东约400m处(35°58′57.9″N,114°13′20.6″E),剖面长约31km(图 1)。采用炸药作为震源激发地震波,药量24kg,井深25—30m,使用固有频率为10Hz的检波器串(12个/道,线性组合)接收地震波,地震仪使用法国Sercel公司生产的408UL数字地震仪。
中深层地震反射剖面ML2位于汤阴地堑南部新乡市附近,探测方向为自西向东,长度约35km(图 1),起点位于新乡市凤泉区陈堡村村北晋新高速附近(35°22′59.8″N,113°51′15.4″E),终点位于新乡市延津县胡庄村村东(35°16′44.9″N,114°13′22.2″E)。采用炸药作为震源激发地震波,药量40kg,井深25—30m,使用固有频率为10Hz的检波器串(12个/道,线性组合)接收地震波,地震仪使用加拿大GEO-X公司生产的ARIES数字地震仪。
中深层地震反射剖面可获得汤阴地堑6—8km以浅的地壳结构及隐伏断裂的位置与性质,但对于隐伏活动断裂的近地表结构还未能准确判断,因此本文采用中深层地震剖面附近的高分辨的浅层地震反射剖面。
浅层地震反射剖面SL1位于中深层地震反射剖面ML1南约1.5km,剖面长约3.3km。采用美国生产的M612-18型可控震源激发地震波,扫描频率30—200Hz,扫描长度8s。使用固有频率为60Hz的检波器接收地震波,地震仪使用德国SUMMIT数字地震仪。
浅层地震反射剖面SL2位于中深层地震反射剖面ML2北约1km,剖面长约4km。采用美国产M612-18型可控震源激发地震波,扫描频率20—200Hz,扫描长度8s。使用固有频率为60Hz的检波器接收地震波,地震仪使用德国SUMMIT数字地震仪。
不同的探测目的需要采用不同的探测方法和观测系统参数(刘保金等,2008)。表 1为本研究地震反射剖面的观测系统参数。
表 1 地震反射剖面观测系统参数表Table 1. Observation system parameters of the seismic reflect profiles剖面名称 道间距/m 炮间距/m 接收道数 覆盖次数 剖面长度/km ML1 30 120 280 35 31 ML2 40 280 800 57 35 SL1 2 10 200 20 3.3 SL2 4 16 160 20 4 3. 中深层地震反射剖面特征
3.1 中深层地震反射剖面ML1
图 2给出的中深层地震反射时间剖面ML1地层界面反射丰富,多组反射波同相轴可以被连续、可靠追踪。从图中可以看出,剖面桩号7km为界东西两侧反射波震相明显不同。7km以东的内黄隆起区可看到3组横向上可连续追踪、纵向上层组关系清晰的界面反射波(TN1、TN、Tg),新近系内部界面反射波TN1基本呈水平展布、局部略有起伏变化,新近系底界面反射波TN向上隆起的特征明显,剖面桩号6km处其埋深最浅,以东则向下倾伏。说明该区域新近纪早期构造活动剧烈,造成地层隆起、褶皱、倾伏,而后构造活动减弱,沉积环境较为稳定,地层基本呈水平展布。新近系之下为反射能量较弱、界面展布特征与新近系底界面类似的结晶基底反射Tg,其下为一系列弱振幅、不连续、短小的反射震相,说明结晶基底之下岩石的波阻抗差异较小,无法形成有效的、能量较强的界面反射。在剖面桩号7km以西汤阴地堑内沉积盆地反射特征明显,多组沉积层强反射不整合的覆盖在自西向东倾伏的结晶基底反射Tg上,且随着深度的增加,地层倾角增大,显示出明显的多期掀斜运动特征。根据地质资料可知,这套沉积地层反射是来自新生代、中生代和古生代的地层界面反射,且某些地质时期地层缺失、不整合接触明显,可能与该区域强烈的构造运动与地层沉降有关。
中深层地震反射时间剖面ML1(图 2)所揭示的断裂特征非常清楚,根据剖面上反射波能量差异、同相轴横向连续性及地层产状不同共分为2条断裂(F1、F2),断裂F1是汤东断裂的主断裂,F2是其分支断裂。主断裂F1为西倾铲型正断层,断裂浅部倾角较陡,而后逐渐变缓,断裂向上错断了新近纪地层TN1,向下错断了所有的沉积层,延伸至结晶基底Tg之下。汤东断裂的分支断裂F2为东倾正断层,与主断裂呈“Y”字型构造特征,断裂F2向上错断了古近纪地层。
3.2 中深层地震反射剖面ML2
图 3给出了中深层地震反射时间剖面ML2。从图中可以看出,TWT1.5s以上的地层界面反射波能量较强、同相轴分段连续性较好,可识别出3组地层界面反射(TQ、TN、Tg)。这些界面反射波在纵向上层组关系清晰,横向上块状结构特征明显。剖面自东向西依次经过内黄隆起、汤阴地堑和太行山隆起。不同的块体剖面所揭示的反射波震相明显不同。剖面桩号25km以东内黄隆起区可识别震相较多,但界面反射波能量相对较弱,横向上地层基本呈水平展布,仅结晶基底Tg局部有起伏变化。剖面桩号6km以西的太行山隆起区仅识别出横向上呈水平展布的新近系底界面反射TN,其反射波能量强,埋深相对较浅,反射波TN以下为一系列短小、不连续、能量相对较弱的反射震相。剖面中段汤阴地堑内地层自西向东倾伏,反射波同相轴局部扭曲、错断明显,新近系底界面反射TN之下为一些横向延续性较短、能量相对较弱反射震相。
根据剖面上反射波同相扭曲、错断及地层产状的不同,共分为3条断裂(F3、F4、F5),根据区域地质资料,可知F3是汤东断裂,F4是汤中断裂,F5是汤西断裂。从图中可以看出,汤东断裂F3为西倾铲型正断层,断裂向上错断第四系底界面反射TQ,向下错断新近系底界面TN,延伸至结晶基底之下的变质岩系中。汤西断裂F5为东倾正断层,断裂向上错断新近系底界面TN,在汤西断裂附近TWT2s以下,可以看到能量较弱、以较大角度向东倾伏的断面波,说明断裂两侧波阻抗存在差异,地层岩性可能不同,从而在断裂面上形成了弱能量的界面反射波。汤中断裂F4为东倾正断层,断裂向上错断新近系底界面TN。
4. 浅层地震反射剖面特征
4.1 浅层地震反射剖面SL1
从图 4的浅层地震反射剖面SL1可以看出,剖面横向上不同地段的反射震相特征明显不同。剖面桩号2500m以东,第四系底界面TQ基本呈水平展布,新近系沉积较薄,新近系底界面TN自东向西倾伏,TN之下为一些横向延续性较短、反射能量较弱的反射震相,根据地质资料推测其可能是来自太古代或元古代变质岩系的反射。剖面桩号1200m以西,第四系底界面TQ自东向西倾伏,新近纪地层呈水平展布,且沉积较厚,新近系之下为连续性较好的古近纪沉积层。剖面中段可以看到界面反射波扭曲、错断非常明显,可知是由于汤东断裂的活动造成附近地层起伏、褶皱,从而使得地层界面反射波出现了扭曲、变化。从图中可以看出,汤东断裂的近地表是由2条东倾正断裂与3条西倾正断裂组成。断裂FP2是汤东断裂的主断裂,倾向为西,断裂错断了浅部地层TQ,上断点埋深较浅,向下依次错断了新近系、古近系并延伸至剖面之外。断裂FP2西侧有3条次级正断裂,与主断裂形成“Y”型构造。断裂FP1是汤东断裂的分支断裂,倾角较陡,向上延伸至新近系内部,上覆的第四系底界面反射波连续性较好,并未发现错断。通过跨主断裂FP2的钻孔联合地质剖面揭示其上埋深为20m,断裂错断了中更新世中晚期地层,因此该区域汤东断裂为中更新世活动断裂(刘保金等,2012)。
图 4 SL1浅层地震反射时间剖面和深度解释剖面(刘保金等,2012)Figure 4. Stacked section and depth profile of the shallow seismic reflection in SL1(Liu et al., 2012)4.2 浅层地震反射剖面SL2
图 5的浅层地震反射剖面SL2具有较高的分辨率,且地层界面反射较为丰富。根据本区域地质及钻孔资料可知,反射波T01为第四系内部的界面反射,TQ为第四系底界面反射。剖面清晰的揭示了汤东断裂浅部特征,在断裂附近可清晰看到同相轴扭曲、错断及断点的绕射波,且断裂两侧反射波组特征与界面起伏形态差异明显。在断裂上盘,地层界面反射非常丰富,第四纪地层基本呈水平展布,新近纪地层反射波能量较强,局部呈凹陷特征,尤其以TWT500ms附近地层界面反射最为明显。在断裂下盘,地层界面反射明显较少,且所有地层在断裂附近均呈向上隆起特征。汤东断裂向上错断浅部地层T01,上断点埋深约40m,根据跨断裂完成的钻孔联合地质剖面,确定该区域汤东断裂为晚更新世活动断裂(徐增波等,2019)。
图 5 SL2浅层地震反射时间剖面和深度解释剖面(徐增波等,2019)Figure 5. Stacked section and depth profile of the shallow seismic reflection in SL2(Xu et al., 2019)5. 结论
本文利用中深层地震反射剖面获得汤阴地堑北部与南部的上地壳结构、新生代沉积层分布特征以及断裂特征与性质,同时利用浅层地震反射剖面获得汤东断裂的近地表特征。研究结果为分析汤阴地堑的形成与演化、汤东断裂特征与活动性研究提供了科学依据。
已有资料表明,汤阴地堑是夹持在太行山隆起与内黄隆起间的断陷型盆地。中深层地震反射剖面同样揭示了这种盆岭相间的构造格局,测线控制范围内汤阴地堑与其两侧块体的新生代沉积厚度明显不同,汤阴地堑内沉积较厚的新生代地层,而两侧隆起区新生代地层明显较薄,地堑内新生代地层自西向东倾伏且新生代沉积层逐渐变厚,最厚处出现在汤东断裂附近,表明研究区汤东断裂控制并影响汤阴地堑新生代时期的活动与演化。内黄隆起区新生代沉积层较薄,结晶基底Tg埋深较浅,结晶基底之上直接覆盖新近纪沉积层,缺失古生代、中生代地层,说明测线控制范围内黄隆起可能经过了多期隆起抬升,地层受到了分化与剥蚀。太行山隆起区新生代沉积相对最薄,局部太古界变质岩系直接出露。
地堑北部与南部的地壳浅部结构与构造差异明显。地堑北部表现为由汤东断裂控制的箕型半地堑,地堑内多组沉积地层不整合的覆盖在自西向东倾伏的结晶基底Tg上,且随着深度的增加,地层倾角增大,显示出明显的掀斜运动与多期构造活动特征。汤阴地堑南部表现为由汤东、汤西断裂共同控制的断陷型地堑,地堑内地层界面反射波自西向东倾伏,新近纪地层底界面表现出强反射震相特征。结合已有地质资料,认为汤东断裂是汤阴地堑的主控边界断裂,控制着汤阴地堑演化与新生代地层的沉积,断裂为走向NE,倾向NW的铲型正断层,向上错断了第四纪地层,向下错断所有沉积地层,延伸至结晶基底之下的变质岩系中。汤西断裂为走向NE、倾向SE的正断层,向北延伸至汤阴宜沟附近,北部中深层地震反射剖面上未发现该断裂的存在。
致谢: 本文的写作与资料解释得到了刘保金研究员的指导,在此表示衷心的感谢!同时非常感谢审稿专家,使得本文更加完善! -
图 2 Anderson等(1984)经典计算方法计算的kappa值
Figure 2. Anderson&Hough's the classical method calculates k
表 1 研究区域地震台站信息
Table 1. The research area seismic station information
序号 台站名称 经度 纬度 记录器型号 场地类型 1 巴彦浩特 105.7°E 38.8°N ETNA 土层 2 巴彦木仁 106.7°E 39.9°N ETNA 土层 3 宝丰 106.3°E 39.0°N MR-2002 土层 4 磴口 106.9°E 40.3°N ETNA 土层 5 干盐池 105.3°E 36.6°N MR-2002 土层 6 公地 106.8°E 40.6°N ETNA 土层 7 吉兰泰 105.7°E 39.7°N ETNA 土层 8 临河 107.6°E 40.8°N ETNA 土层 9 灵武 106.3°E 38.1°N MR-2002 土层 10 前进农场 106.4°E 38.8°N MR-2002 土层 11 青铜峡 106.0°E 38.0°N MR-2002 基岩 12 沙海 106.9°E 40.9°N ETNA 土层 13 乌海 106.8°E 39.4°N ETNA 土层 14 西吉 105.4°E 35.5°N MR-2002 土层 表 2 kappa值与震中距拟合结果
Table 2. The relationship between kappa and epicentral distance
台站名称 方向 回归方程 K0/s 拟合优度值/R2 $ {{\bar{{K}}}_{\rm{0}}}$/s VS20/m·s−1 VS30/m·s−1 巴彦浩特 N-S kappa=0.000 3R+0.0264 0.0264 0.9803 0.02530 343.05 370.2 E-W kappa=0.000 3R+0.0242 0.0242 0.9854 巴彦木仁 N-S kappa=0.000 2R+0.0363 0.0363 0.8691 0.03485 353.4 389.4 E-W kappa=0.000 2R+0.0334 0.0334 0.9280 宝丰 N-S kappa=0.000 1R+0.0687 0.0687 0.6055 0.06560 180.4 222.1 E-W kappa=0.000 1R+0.0625 0.0625 0.7876 磴口 N-S kappa=0.000 1R+0.0605 0.0605 0.8971 0.05845 231.7 261.7 E-W kappa=0.000 1R+0.0564 0.0564 0.9425 干盐池 N-S kappa=0.000 1R+0.0531 0.0531 0.7542 0.05245 317.4 357.3 E-W kappa=0.000 1R+0.0518 0.0518 0.7194 公地 N-S kappa=0.000 2R+0.0541 0.0541 0.9790 0.05450 222.7 251.2 E-W kappa=0.000 2R+0.0549 0.0549 0.9730 吉兰泰 N-S kappa=0.000 2R+0.0526 0.0526 0.8812 0.05525 327.7 335.8 E-W kappa=0.000 2R+0.0579 0.0579 0.7826 临河 N-S kappa=0.000 1R+0.0530 0.0530 0.9977 0.05175 232.3 249.7 E-W kappa=0.000 1R+0.0505 0.0505 0.9960 灵武 N-S kappa=0.000 3R+0.0449 0.0449 0.9996 0.04415 292.6 325.1 E-W kappa=0.000 3R+0.0434 0.0434 0.9957 前进农场 N-S kappa=0.000 3R+0.0637 0.0637 0.7536 0.06555 209.5 248.9 E-W kappa=0.000 3R+0.0674 0.0674 0.6815 青铜峡 N-S kappa=0.000 1R+0.0428 0.0428 0.4987 0.04315 219.0 260.0 E-W kappa=0.000 1R+0.0435 0.0435 0.6287 沙海 N-S kappa=0.000 2R+0.0489 0.0489 0.9743 0.04870 188.3 216.5 E-W kappa=0.000 2R+0.0485 0.0485 0.9710 乌海 N-S kappa=0.000 5R+0.0155 0.0155 0.7320 0.01545 352.7 390.5 E-W kappa=0.000 5R+0.0154 0.0154 0.6702 西吉 N-S kappa=0.000 1R+0.0429 0.0429 0.4713 0.04240 195.1 228.7 E-W kappa=0.000 1R+0.0419 0.0419 0.4909 表 3 巴彦木仁、磴口台站kappa值与震中距拟合结果
Table 3. Fitting results of kappa value and epicenter distance at Bayanmuren and Dengkou station
台站名称 方向 回归方程 K0/s 拟合优度值/R2 巴彦木仁 N-S kappa=0.000 3R+0.0275 0.0275 0.6186 E-W kappa=0.000 3R+0.0266 0.0266 0.7081 磴口 NS kappa=0.000 1R+0.0650 0.0650 0.2447 EW kappa=0.000 1R+0.0562 0.0562 0.6064 -
[1] 白立新, 成云辉, 张杰等, 2019. 河北永清M4.3地震北京烈度仪台网记录分析. 震灾防御技术, 14(1): 210—219. doi: 10.11899/zzfy20190120Bai L. X., Cheng Y. H., Zhang J., et al., 2019. Analysis of the recoed of Beijing seismic intensity meter network for the M 4.3 Yongqing earthquake. Technology for Earthquake Disaster Prevention, 14(1): 210—219. (in Chinese) doi: 10.11899/zzfy20190120 [2] 崔建文, 李正光, 赵云旭, 2007. 2007年宁洱6.4级地震强震动观测记录. 地震研究, 30(4): 384—388. doi: 10.3969/j.issn.1000-0666.2007.04.013Cui J. W., Li Z. G., Zhao Y. X., 2007. Strong motion observation records of the 2007 Ning’er, Yunnan, MS6.4 earthquake. Journal of Seismological Research, 30(4): 384—388. (in Chinese) doi: 10.3969/j.issn.1000-0666.2007.04.013 [3] 高祥林, 马晓静, 李晓丽, 2010. 亚洲东部“大三角”地震构造区的周边和深部动力环境. 地学前缘, 17(4): 33—42.Gao X. L., Ma X. J. Li X. L., 2010. A surrounding and deep dynamic context of the great triangle-shaped seismic region in the eastern Asia continent. Earth Science Frontiers, 17(4): 33—42. (in Chinese) [4] 李文倩, 2014. 用数字地震台网小震数据建立兰州和华北地区强地震动衰减关系. 哈尔滨: 中国地震局工程力学研究所.Li W. Q., 2014. Strong ground motion attenuation relationships for Lanzhou and North China from small earthquake records by China earthquake networks. Harbin: Institute of Engineering Mechanics, China Earthquake Administration. (in Chinese) [5] 陶正如, 陶夏新, 2015. 美国2014地震区划中采用的地震动衰减关系. 世界地震工程, 31(3): 78—84.Tao Z. R., Tao X. X., 2015. Ground motion attenuation relationships adopted in 2014 update of the US national seismic hazard maps. World Earthquake Engineering, 31(3): 78—84. (in Chinese) [6] 万永革, 2012. 数字信号处理的MATLAB实现. 2版. 北京: 科学出版社. [7] 朱百慧, 2016. 强震数据Kappa值影响因素的讨论. 哈尔滨: 中国地震局工程力学研究所.Zhu B. H., 2016. The influence factors of Kappa from strong motion records. Harbin: Institute of Engineering Mechanics, China Earthquake Administration. (in Chinese) [8] Anderson J. G., Hough S. E., 1984. A model for the shape of the Fourier amplitude spectrum of acceleration at high frequencies. Bulletin of the Seismological Society ofAmerica, 74(5): 1969—1993. [9] Anderson J. G., 1986. Implication of attenuation for studies of the earthquake source. American: Earthquake Source Mechanics. American Geophysical Union, 37: 311—318. [10] Anderson, J. G., Humphrey J. R. Jr., 1991. A least squares method for objective determination of earthquake source parameters. Seismological Research Letters, 62(3—4): 201—209. [11] Ktenidou O. J., Gélis C., Bonilla L. F., 2013. A study on the variability of Kappa (κ) in a borehole: implications of the computation process. Bulletin of the Seismological Society of America, 103(2A): 1048—1068. doi: 10.1785/0120120093 [12] Sun X. D., Tao X. X., Duan S. S., et al., 2013. Kappa (k) derived from accelerograms recorded in the 2008 Wenchuanmainshock, Sichuan, China. Journal of Asian Earth Sciences, 73: 306—316. doi: 10.1016/j.jseaes.2013.05.008 期刊类型引用(0)
其他类型引用(3)
-