云南禄劝地震部分余震的矩张量反演
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. -
期刊类型引用(3)
1. 靳炎锡，秦满忠，李少华. 利用前驱波研究地幔间断面的若干进展. 地震工程学报. 2023(01): 209-219+243 . 
百度学术
2. 崔辉辉，高雅健，周元泽. 伊豆-小笠原地区岩石圈软流圈边界地震学证据. 科学通报. 2017(07): 711-720 . 百度学术
3. 秦满忠，刘旭宙，邹锐，张元生，郭晓，魏从信，王亚红，孙点峰. 甘肃祁连山大容量气枪主动源最大探测范围. 地震工程学报. 2017(06): 1070-1075 . 百度学术
其他类型引用(3)
下载:
计量
- 文章访问数: 1527
- HTML全文浏览量: 17
- PDF下载量: 141
- 被引次数: 6
