2. 中国地震台网中心, 北京 100045;
3. 中国地震局地震预测研究所, 北京 100036
2. China Earthquake Networks Center, Beijing 100045, China;
3. Institute of Earthquake Forecasting, China Earthquake Administration, Beijing 100036, China
2017年8月8日21时19分,四川省阿坝州九寨沟县发生M7.0级地震,截至2017年9月18日,共发生余震7900多次.中国地震台网中心(CENC)发布的主震速报参数为:33.20°N,103.82°E,20 km,M7.0.震源位置参数是确定发震断层的重要依据,也是地震学其他反演问题的基本参数.通常地震定位方法有两种,线性定位和非线性定位.这两种方法都是基于P波或S波走时方程的盖戈方法(Geiger, 1912)来反演位置参数.盖戈方法在线性化,最小二乘求解以及观测系统的误差分布上有其局限性(许力生等, 2013a, 2013b).本研究将利用地震波干涉法对干涉波形进行偏移叠加,从而进行地震定位.
物理学中的干涉指两列或两列以上波形信号经过相互叠加产生新的波列的一种现象.地震学中的干涉法与之类似,指的是对地震波信号进行互相关、卷积或反卷积等计算得到新的波形信号(Campillo and Paul, 2003; Claerbout, 1968; Duvall et al., 1993; Lobkis and Weaver, 2001; Rickett and Claerbout, 1999; Shapiro and Campillo, 2004; Snieder, 2004; van Manen et al., 2005, 2006, 2007; Wapenaar, 2004; Weaver and Lobkis, 2001; Galetti and Curtis, 2012).干涉后的地震信号不仅包含了原始信号的特性,而且能够反映出原始信号不具备的一些特征(Sun, 2013; 陶毅等,2010).在可控源勘探领域通过地震干涉法,可以将VSP(垂直地震剖面)数据转换为SWP(单井地震剖面)数据(Bakulin and Calvert, 2004; Calvert et al. 2004; Bakulin and Calvert, 2006),VSP数据转换为SSP(平地表地震剖面)数据(Rickett and Claerbout, 1999)等.地震干涉法在天然地震领域也有着广泛应用.背景噪声信号做互相关干涉叠加可以对地球介质进行无源成像(Shapiro et al., 2005).对地震尾波做互相关干涉,也可以重构台站对间的面波格林函数(Campillo and Paul, 2003).无源干涉成像还可以检测地震波速度的动态变化(Sens-Schönfelder and Wegler, 2006; Wegle and Sens-Schönfelder, 2007; Wegler et al., 2009; 肖卓和高原,2016).
地震波干涉法还可以用来进行震源成像,确定震源位置参数(Schuste and Rickett, 2000; Schuster et al., 2004).本研究将利用地震波干涉法,确定九寨沟地震序列主震及部分余震震源位置参数,并对发震断裂开展相关讨论.
1 背景和数据资料九寨沟M7.0地震发生在青藏高原东缘的巴颜喀拉地块内(图 1).巴颜喀拉地块属于青藏地块的次级块体之一.晚新生代以来,青藏高原的迅速隆起,导致其东缘地区地壳物质沿大型走滑断裂系发生大规模的向东、东南方向的“逃逸”运动(Tapponnier et al., 1982),在此过程中,受到华南地块的阻挡,产生纵向的水平挤压.区域内构造应力场自晚新生代以来转变为NW-SE向的近水平挤压.通过震源机制解资料和剪切波分裂参数反演得到巴颜喀拉地块东部的主压应力方向为NWW-近EW向(Zhao et al., 2012; 赵博等,2013;高原等,2013).四川盆地—龙门山—巴颜喀拉一线的地震勘探剖面显示,该地区地壳厚度在龙门山两侧变化很大(滕吉文等, 2008; 朱介寿等,1984;Zhang et al., 2010).川西—青藏交汇区地震活动背景一直很强(Gao et al., 2000),2008年M8.0汶川地震和2013年M7.0芦山地震均发生在该区域东部的龙门山断裂带上(Zhao et al., 2011;赵博等,2013; 易桂喜等,2016).
本次地震发生在虎牙断裂(F1)、岷江断裂(F3)、塔藏断裂(F4)和雪山梁子断裂(F2)四条断裂所围限的区域(图 1).震中距龙门山后山断裂的平武—青川段(F6)的直线距离约为100 km.
研究利用了区域内布设的国家地震台网台站和四川省、甘肃省区域台网台站波形数据,其中距主震最近的台站为JZG(九寨沟台),震中距约为40 km.所用台站分布如图 1,共用到震中附近18个台站的波形数据,计算得到台站对干涉波形数据为60个.
2 地震干涉方法 2.1 互相关干涉和卷积干涉互相关干涉和卷积干涉是应用比较广泛的地震干涉法,例如对地震尾波等随机波场数据做互相关干涉,可以重构台站间格林函数.本研究与随机场(如地震尾波)互相关干涉不同,是对主波震相进行互相关或卷积干涉.
如图 2,假设真实震源位置为S,震源S到台站A、B的走时分别为TSA、TSB.两个台站的波形信号的时间序列分别为f(t)、h(t),在频率域为F(ω)、H(ω),对两列信号做互相关和卷积:
(1) |
(2) |
其中
式(1)中相位谱eiω(TSB-TSA)表示互相关干涉后的震相到时,相当于震源到两台站的走时差TSA-TSB.式(2)的相位谱eiω(TSB+TSA)表示卷积干涉后的震相到时,相当于震源到两台站的走时之和.
2.2 干涉数据偏移叠加以互相关型干涉为例,两列波形互相关运算后得到的干涉波形数据为Φ(A, B),对相干数据进行偏移,与常规的标准偏移相同,选择e-iω(TSxB-TSxA)作为偏移核函数(Schuster et al., 2004),对所有频率求和,并对偏移后的所有台站数据进行叠加.如图 2,计算所有空间点Sx,得到任意一点的偏移值:
(3) |
当Sx为实际震源位置S时有最大偏移振幅值(Schuster et al., 2004).
如果对原始波形数据做卷积干涉,并选择偏移核函数为e-iω(TSB+TSA)时,则为卷积型干涉数据偏移叠加.
3 数据处理及定位结果天然地震和人工可控源爆破在波形上有很大不同.天然地震震源机制比人工爆破的震源机制要复杂很多.位于不同方位角的台站,由于震源辐射花样不同,相同震相的振幅正负值不同.为了避免干涉叠加时由于震源在不同方位辐射花样产生振幅正负差异,在做干涉计算前,首先计算波形的特征函数.本研究利用长短时窗能量比来计算原始波形的特征函数.
在时间序列第i个点前后分别选取一个短时窗和长时窗,计算两个时窗内振幅能量和,用长时窗能量比短时窗能量,作为特征函数(如图 3).图 3为REG, WCA和XHA三个台站主震波形和特征函数,特征函数很好地反映出初至到时,并且统一了初动方向,使不同辐射方位初动一致.
通过干涉偏移,确定震中水平位置(即经纬度).对每个台站记录到的波形数据计算其特征函数,然后进行干涉处理,得到干涉波形并偏移(如图 4).图 4a为互相关干涉波形,反映了同一震源到两个台站的到时差(时间序列负值表示互相关运算主信号为走时短的波形),用e-iω(TSxB-TSxA)作为偏移核函数,进行偏移处理.图 4b为偏移后的波形,图 4b显示偏移后的干涉波形在零时刻对齐.
将水平方向分为0.01°×0.01°的格网,对每个网格节点计算干涉偏移振幅, 并将所有台站偏移波形进行叠加,振幅能量最大的格网点为震中位置,如图 4c,粉色五角星位置(103.81°E, 33.21°N)为本次主震震中位置.
地震定位中,震源深度的确定比较困难,而不同资料和不同方法得到的震源深度的含义是有差异的(高原等,1997).互相关型干涉,等价于波场传播路径变短,干涉波形的到时为两个台站的走时差.本研究发现,互相关型干涉的到时差对震源深度变化不敏感.为了得到较为可靠的震源深度,我们利用卷积型干涉,使波场传播路径增长.采用三维速度模型(梁建宏等,2015)计算偏移时间,增加干涉偏移的准确性.如图 5,图 5a为卷积型干涉波形,采用e-iω(TSB+TSA)为偏移核函数进行偏移处理,图 5b为偏移后的干涉波形图.将深度划分为1 km的间隔,对每个节点计算干涉偏移振幅并叠加.图 5c中,干涉偏移振幅最大的点为9 km,即为主震震源深度.
本次九寨沟M7.0地震及其余震沿虎牙断裂西北端的一条无名断裂分布.主震发生后,截止2017年9月18日16时00分,共发生7900多次余震,其中,M3.0~3.9级地震35个,M4.0~4.9级地震3个,无M5.0级以上地震发生.本次地震序列主要集中发生在震后三天,即8月8日—8月10日三天,其中M>3.0级地震占整个序列M3.0级以上地震的70%.利用地震干涉方法,对本次地震序列中全部M>4.0级地震(共3个)和全部M3.5~3.9级地震(共8个)进行干涉定位,选择这些地震主要考虑其震级较大,波形数据信噪比高,用来偏移叠加的特征函数比较准确.定位结果见表 1,地震分布见图 6.图 6a显示主震及12次较大余震均发生在虎牙断裂北部、岷江断裂与塔藏断裂的交汇区附近.该区域可能存在一条NW向无名断层,主震和余震序列呈现沿该断层分布,这条断裂在九寨沟地震后被命名为树正断裂.图 6b和6c为A1-A2,B1-B2两条地震剖面.本研究得到九寨沟地震主震震源位置为(103.81°E,33.21°N,9 km),定位时共用到了67个道集数据进行偏移叠加.当叠加了7~10个道集数据时,定位结果趋于稳定.为了讨论定位精度问题,我们取叠加10道数据以后,每次测量结果进行评估,得到水平经纬度方向和深度方向的均方根误差RMSE为(0.02°,0.03°,1.4 km).表 1给出了干涉定位结果和CENC、USGS(美国地质调查局)、GFZ(德国地学研究中心)等机构的定位结果.对比本研究和国内外其他机构的结果,以主震位置为例,各结果在水平方向差别都在0.1°以内,但在深度方向,本结果与CENC速报结果有一定差异,但与USGS和GFZ结果基本一致.
利用地震干涉法做震源成像,确定震源位置,在人工源勘探领域应用比较广泛.由于天然地震震源机制和辐射花样的复杂性,该方法在天然地震领域应用较少.本研究通过计算天然地震原始波形特征函数,利用地震干涉法和偏移叠加技术,对天然地震震源成像,从而确定震源位置.通过对比本研究结果与CENC,USGS,GFZ等地震学研究机构发布的该次地震参数信息,表明地震干涉法在天然地震震源成像,确定震源位置方面是可行的,其定位结果是可靠的.
互相关型干涉等效于地震波传播路径之差,相应的其偏移核函数也反映了走时差.互相关干涉与双差定位有相似之处,通过互相关型干涉,一定程度上消除了传播路径上由于速度扰动而引起的位置偏差.与双差定位(Waldhauser and Ellsworth, 2000)不同的是,该方法可以对单个地震进行绝对定位,不需要初始位置,避免了因初始位置不准确引起定位结果受到地震丛约束.卷积型干涉等效于地震波传播路径之和,相应地,其偏移核函数反映了走时和.为了对震源深度进行约束,在近台稀少时,卷积型干涉由于增长了传播路径,等效于波场传播路径在垂直方向上延长,可以在一定程度用稍远一些的台站约束深度.
本研究得到九寨沟地震主震震源位置为(103.81°E,33.21°N,9 km)该结果与CENC速报结果(103.82°E,33.20°N,20 km),在水平方向差别不大,但深度方向有一定差异.USGS和GFZ发布的该地震位置参数深度均为9 km,与本研究结果一致.GCMT发布的矩张量解矩心深度为14 km(http://www.globalcmt.org/CMTsearch.html),GFZ发布的矩张量解矩心深度为10 km(http://www.seiscomp.com/data/alerts/2017/gfz2017pldo/mt.txt).这些结果显示,起始破裂点和矩心深度都在9~15 km的上地壳.本文利用卷积型干涉法在深度方向进行偏移成像,得到震源位置为9 km.
通过对主震和几次较大余震定位结果显示,地震序列呈NW-SE向分布,深度方向呈近似垂直分布.GCMT震源机制解为,节面1:走向150°,倾角78°,滑动角-13°;节面2:走向242°,倾角77°,滑动角-166°.根据本次主震及余震序列定位结果,推断本次地震断层面走向为150°.深度方向上余震序列呈垂直分布,与震源机制高倾角走滑断层一致.水平定位结果显示,震中位于虎牙断裂以北,岷江断裂与塔藏断裂的交汇区,区域内断层分布复杂,有一条无名断裂(图 6中F*断裂),主震及M>3.5余震其空间分布均与该无名断裂接近.结合主震震源机制解,推断这次九寨沟地震主震的发震断层在塔藏断裂与岷江断裂交汇处的南侧、大致在虎牙断裂的西北方向延长线上的一个无名断裂上,该断裂在这次九寨沟地震后被命名为树正断裂(易桂喜等,2017).这次地震的余震序列呈NW-SE向分布,与虎牙断裂走向基本一致.树正断裂的活动给出了一个信号,从断裂的走向上看,树正断裂与虎牙断裂的西北端相对应,推测两个断裂可能有连接的趋势.
Bakulin A, Calvert R. 2004. Virtual source: New method for imaging and 4D below complex overburden. //SEG Technical Program Expanded Abstracts 2004. SEG, 2477-2480.
|
Bakulin A, Calvert R. 2006. The virtual source method:Theory and case study. Geophysics, 71(4): SI139-SI150. DOI:10.1190/1.2216190 |
Calvert R W, Bakulin A, Jones T C. 2004. Virtual sources, a new way to remove overburden problems. //66th EAGE Conference & Exhibition. EAGE.
|
Campillo M, Paul A. 2003. Long-range correlations in the diffuse seismic coda. Science, 299(5606): 547-549. DOI:10.1126/science.1078551 |
Claerbout J F. 1968. Synthesis of a layered medium from its acoustic transmission response. Geophysics, 33(2): 264-269. DOI:10.1190/1.1439927 |
DuvallT L Jr, Jeffferies S M, Harvey J W, et al. 1993. Time-distance helioseismology. Nature, 362(6419): 430-432. DOI:10.1038/362430a0 |
Galetti E, Curtis A. 2012. Generalised receiver functions and seismic interferometry. Tectonophysics, 532-535: 1-26. DOI:10.1016/j.tecto.2011.12.004 |
Gao Y, Zhou H L, Zheng S H, et al. 1997. Preliminary discussion on implication of determination on source depth of earthquake. Earthquake Research in China, 13(4): 321-329. |
Gao Y, Wu Z L, Liu Z, et al. 2000. Seismic source characteristics of nine strong earthquakes from 1988 to 1990 and earthquake activity since 1970 in the Sichuan-Qinghai-Xizang (Tibet) zone of China. Pure Appl. Geophys, 157(9): 1423-1443. DOI:10.1007/PL00001127 |
Gao Y, Wang Q, Zhao B, et al. 2014. A rupture blank zone in middle south part of Longmenshan Faults:Effect after Lushan MS7.0 earthquake of 20 April 2013 in Sichuan, China. Science China:Earth Sciences, 57(9): 2036-2044. DOI:10.1007/s11430-014-4827-2 |
Geiger L. 1912. Probability method for the determination of earthquake epicenters from the arrival time only. Bull. St. Louis Univ, 8(1): 60-71. |
Liang J H, Han X J, Sun L, et al. 2015. Construction of regional three-dimensional travel-time tables of first P arrival for China Seismograph Network. Acta Seismologica Sinica, 37(1): 125-133. |
Lobkis O I, Weaver R L. 2001. On the emergence of the Green's function in the correlations of a diffuse field. The Journal of the Acoustical Society of America, 110(6): 3011-3017. DOI:10.1121/1.1417528 |
Rickett J, Claerbout J. 1999. Acoustic daylight imaging via spectral factorization; helioseismology and reservoir monitoring. The Leading Edge, 18(8): 957-960. DOI:10.1190/1.1438420 |
Schuster G T, Rickett J. 2000. Daylight imaging in V(x, y, z) media. Utah Tomography and Modeling-Migration Project Midyear Report and SEP Report at Stanford University. |
Schuster G T, Yu J, Sheng J, et al. 2004. Interferometric/daylight seismic imaging. Geophys. J. Int, 157(2): 838-852. DOI:10.1111/gji.2004.157.issue-2 |
Sens-Schönfelder C, Wegler U. 2006. Passive image interferometry and seasonal variations of seismic velocities at Merapi Volcano, Indonesia. Geophys. Res. Lett, 33(21): L21302. DOI:10.1029/2006GL027797 |
Shapiro N M, Campillo M. 2004. Emergence of broadband Rayleigh waves from correlations of the ambient seismic noise. Geophysical Research Letters, 31(7): L07614. DOI:10.1029/2004GL019491 |
Shapiro N M, Campillo M, Stehly L, et al. 2005. High-resolution surface-wave tomography from ambient seismic noise. Science, 307(5715): 1615-1618. DOI:10.1126/science.1108339 |
Snieder R. 2004. Extracting the Green's function from the correlation of coda waves:A derivation based on stationary phase. Physical Review E, 69(4): 046610. DOI:10.1103/PhysRevE.69.046610 |
Sun Y H. 2013. Geophysics bright spots. The Leading Edge, 32(6): 708-709. DOI:10.1190/tle32060708.1 |
Tao Y, Fu L Y, Sun W J, et al. 2010. A review of seismic interferometry. Progress in Geophys, 25(5): 1775-1874. DOI:10.3969/j.issn.1004-2903.2010.05.035 |
Tapponnier P, Peltzer G, Dain A Y L, et al. 1982. Propagating extrusion tectonics in Asia:New insights from simple experiments with plasticine. Geology, 10(12): 611-616. DOI:10.1130/0091-7613(1982)10<611:PETIAN>2.0.CO;2 |
Teng J W, Bai D H, Yang H, et al. 2008. Deep processes and dynamic responses associated with the Wenchuan MS8.0 earthquake of 2008. Chinese J. Geophys, 51(5): 1385-1402. |
van Manen D J, Robertsson J O A, Curtis A. 2005. Modeling of wave propagation in inhomogeneous media. Physical Review Letters, 94(16): 164301. DOI:10.1103/PhysRevLett.94.164301 |
van Manen D J, Curtis A, Robertsson J O A. 2006. Interferometric modeling of wave propagation in inhomogeneous elastic media using time reversal and reciprocity. Geophysics, 71(4): SI47-SI60. DOI:10.1190/1.2213218 |
van Manen D J, Robertsson J O A, Curtis A. 2007. Exact wave field simulation for finite-volume scattering problems. The Journal of the Acoustical Society of America, 122(4): EL115-EL121. DOI:10.1121/1.2771371 |
Waldhauser F, Ellsworth L W. 2000. A double-difference earthquake location algorithm:method and application to the Northern Hayward Fault, California. Bulletin of the Seismological Society of America, 90(6): 1353-1368. DOI:10.1785/0120000006 |
Wapenaar K. 2004. Retrieving the Elastodynamic Green's function of an arbitrary inhomogeneous medium by cross correlation. Physical Review Letters, 93(25): 254301. DOI:10.1103/PhysRevLett.93.254301 |
Weaver R L, Lobkis O I. 2001. Ultrasonics without a source:thermal fluctuation correlations at MHz frequencies. Physical Review Letters, 87(13): 134301. DOI:10.1103/PhysRevLett.87.134301 |
Wegler U, Sens-Schönfelder C. 2007. Fault zone monitoring with passive image interferometry. Geophys. J. Int, 168(3): 1029-1033. DOI:10.1111/gji.2007.168.issue-3 |
Wegler U, Nakahara H, Schönfelder C S, et al. 2009. Sudden drop of seismic velocity after the 2004 MW6.6 Mid-Niigata Earthquake, Japan, observed with passive image interferometry. J. Geophys. Res, 114: B06305. DOI:10.1029/2008JB005869 |
Xiao Z, Gao Y. 2016. Monitoring dynamic variations of crustal media during 2014 Yingjiang double earthquakes by Coda wave interferometry. Progress in Geophysics, 31(6): 2421-2428. DOI:10.6038/pg20160609 |
Xu L S, Du H L, Yan C, et al. 2013a. A method for determination of earthquake hypocentroid:time-reversal imaging technique I-Principle and numerical tests. Chinese J. Geophys, 56(4): 1190-1206. DOI:10.6038/cjg20130414 |
Xu L S, Yan C, Zhang X, et al. 2013b. A method for determination of earthquake hypocentroid:Time-reversal imaging technique-An examination based on people-made earthquakes. Chinese J. Geophys, 56(12): 4009-4027. DOI:10.6038/cjg20131207 |
Yi G X, Long F, Vallage A, et al. 2016. Focal mechanism and tectonic deformation in the seismogenic area of the 2013 Lushan earthquake sequence, southwestern China. Chinese J. Geophys, 59(10): 3711-3731. DOI:10.6038/cjg20161017 |
Yi G X, Long F, Liang M J, et al. 2017. Focal mechanism solutions and seismogenic structure of the 8 August 2017 M7.0 Jiuzhaigou earthquake and its aftershocks, northern Sichuan. Chinese J. Geophys, 60(10): 4083-4097. DOI:10.6038/cjg20171033 |
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 Planet. Sci. Lett, 292(3-4): 254-264. DOI:10.1016/j.epsl.2010.01.046 |
Zhao B, Shi Y T, Gao Y. 2011. Relocation of aftershocks of the Wenchuan MS8.0 earthquake and its implication to seismotectonics. Earthquake Science, 24(1): 107-113. DOI:10.1007/s11589-011-0774-6 |
Zhao B, Shi Y T, Gao Y. 2012. Seismic relocation, focal mechanism and crustal seismic anisotropy associated with the 2010 Yushu MS7.1 earthquake and its aftershocks. Earthquake Science, 25(1): 111-119. DOI:10.1007/s11589-012-0837-3 |
Zhao B, Gao Y, Huang Z B, et al. 2013. Double difference relocation, focal mechanism and stress inversion of Lushan MS7.0 earthquake sequence. Chinese J. Geophys, 56(10): 3385-3395. DOI:10.6038/cjg20131014 |
Zhu J S, Cao J M, Liu S H, et al. 1984. A preliminary investigation on crustal structure of western Sichuan by explosion seismology. Journal of Chengdu College of Geology, (3): 111-122. |
高原, 周蕙兰, 郑斯华, 等. 1997. 测定震源深度的意义的初步讨论. 中国地震, 13(4): 321-329. |
高原, 王琼, 赵博, 等. 2013. 龙门山断裂带中南段的一个破裂空段——芦山地震的震后效应. 中国科学(D辑:地球科学), 43(6): 1038-1046. |
梁建宏, 韩雪君, 孙丽, 等. 2015. 中国地震台网初至P波区域三维走时表的建立. 地震学报, 37(1): 125-133. |
陶毅, 苻力耘, 孙伟家, 等. 2010. 地震波干涉法研究进展综述. 地球物理学报, 25(5): 1775-1874. DOI:10.3969/j.issn.1004-2903.2010.05.035 |
滕吉文, 白登海, 杨辉, 等. 2008. 2008汶川MS8.0地震发生的深层过程和动力学响应. 地球物理学报, 51(5): 1385-1402. |
肖卓, 高原. 2016. 运用尾波干涉技术监测2014年盈江双震期间地壳介质状态的变化. 地球物理学进展, 31(6): 2421-2428. DOI:10.6038/pg20160609 |
许力生, 杜海林, 严川, 等. 2013a. 一种确定震源中心的方法:逆时成像技术(一)——原理与数值实验. 地球物理学报, 56(4): 1190-1206. DOI:10.6038/cjg20130414 |
许力生, 严川, 张旭, 等. 2013b. 一种确定震源中心的方法:逆时成像技术(二)-基于人工地震的检验. 地球物理学报, 56(12): 4009-4027. DOI:10.6038/cjg20131207 |
易桂喜, 龙锋, VallageA, 等. 2016. 2013年芦山地震序列震源机制与震源区构造变形特征分析. 地球物理学报, 59(10): 3711-3731. DOI:10.6038/cjg20161017 |
易桂喜, 龙锋, 梁明剑, 等. 2017. 2017年8月8日九寨沟M7.0地震及余震震源机制解与发震构造分析. 地球物理学报, 60(10): 4083-4097. DOI:10.6038/cjg20171033 |
赵博, 高原, 黄志斌, 等. 2013. 四川芦山MS7.0地震余震序列双差定位、震源机制及应力场反演. 地球物理学报, 56(10): 3385-3395. DOI:10.6038/cjg20131014 |
朱介寿, 曹家敏, 刘舜化, 等. 1984. 用人工地震初探川西地区的地壳结构. 成都地质学院学报, (3): 111-122. |