蒙古中部地区SKS分裂研究的初步结果

管见, 吴庆举

管见, 吴庆举. 2014: 蒙古中部地区SKS分裂研究的初步结果. 地震学报, 36(4): 546-556. DOI: 10.3969/j.issn.0253-3782.2014.04.002
引用本文: 管见, 吴庆举. 2014: 蒙古中部地区SKS分裂研究的初步结果. 地震学报, 36(4): 546-556. DOI: 10.3969/j.issn.0253-3782.2014.04.002
Guan Jian, Wu Qingju. 2014: The preliminary SKS splitting results in the central Mongolia. Acta Seismologica Sinica, 36(4): 546-556. DOI: 10.3969/j.issn.0253-3782.2014.04.002
Citation: Guan Jian, Wu Qingju. 2014: The preliminary SKS splitting results in the central Mongolia. Acta Seismologica Sinica, 36(4): 546-556. DOI: 10.3969/j.issn.0253-3782.2014.04.002

蒙古中部地区SKS分裂研究的初步结果

基金项目: 国家国际科技合作项目(2011DFB20210)资助.
详细信息
    通讯作者:

    管见, E-mail: stguanjian@qq.com

  • 中图分类号: P315.2

The preliminary SKS splitting results in the central Mongolia

  • 摘要: 为了进行蒙古中部地区的深部构造观测与模型研究, 中国地震局地球物理研究所和蒙古科学院天文与地球物理研究中心合作, 在该地区架设了60台宽频带流动地震台. 这些台站分布在被额尔齐斯—中蒙古—额尔古纳断裂隔断的加里东、 海西两个构造域. 本文选取了2011年8月—2012年7月部分台站记录的震中距在85°—135°且MW≥5.0的远震事件, 通过SplitLab软件来进行SKS分析. 使用最小能量法、 旋转相关法和最小特征值法进行处理, 取得了蒙古中部地区各向异性的初步结果. 结果显示: 快慢波延时在0.8—2.0 s之间. 西北部区域靠近杭爱山脉与肯特山脉的区域构造走向均为ENE--WSW方向, 而快波方向大多为N40°W, 与区域构造走向近似垂直.根据HS3-NUVEL-1A板块模型计算的该地区板块绝对运动方向与快波偏振方向呈顺时针约15°夹角. 然而, 有个别台站表现出了两种模式的快波分裂, 其中一种与当地ENE--WSW向的断裂方向近似一致.
    Abstract: A passive seismic array with 60 broadband seismometers was deployed in the central Mongolia by the cooperation between Institute of Geophysics, China Earthquake Administration and Research Center of Astronomy and Geophysics, Mongolian Academy of Sciences, to carry out the deep structure observation and modeling study. The passive seismic array was distributed within the Caledonian orogen and the Hercynian orogen separated by the Ergis-central Mongolia-Eerguna fault zone. The teleseismic waveform data from part of the array with the epicentral distance of 85°—135° and MW5.0, from August 2011 to July 2012, was selected to do SKS splitting analysis through the SplitLab package. Three methods, i.e., the minimum energy of tangential component (SC), the rotation and correlation (RC), and eigenvalue (EV), are simultaneously applied for each SKS waveform data. The preliminary anisotropy result was obtained for the central Mongolia. It shows that the time delay between fast and slow waves is within 0.8—2.0 s, and the polarization direction of fast wave is about N40°W, almost perpendicular to the geological structure strike direction of ENE--WSW. The fast direction is about 15° systematically clockwise to the absolute plate motion inferred from the HS3-NUVEL-1A model. It is also noted that not only the two different splitting modes are simultaneously observed beneath certain station, but also the fast wave direction parallel to the ENE--WSW striking direction exists beneath certain station.
  • 随着强震观测技术的发展,强震动时程记录逐步累积,人们可利用概率统计的方法对地震动的强度进行估计,并直接应用于工程结构抗震设计. 地震动是由震源释放出来的地震波引起的地表附近土层的振动,是由各种周期分量按特定比例组成,其频谱组成随场地土壤条件改变,是工程地震与结构抗震之间的桥梁(胡聿贤,2006). 地震动的物理量大多选定为与地震惯性力联系密切的地震动加速度时程. 在结构工程领域,地震对结构破坏能力的大小主要与地震动的幅值、 频谱特性和持时三大特征有关. 如何确定一次地震动强度的指标并用于结构动态时程分析中选取强震记录,一直是结构抗震领域的一个难点,同时也是基于性态抗震设计所需解决的首要问题. 目前来看,只有匹配目标反应谱的选波方法才能同时反映地震动的振幅和频率特征,因而考虑有阻尼的加速度反应谱的谱形匹配为地震动主要依据进行强震记录的选取更为合理. 总之,选择相似地震和地质条件下,实际地震动强度记录的反应谱与目标谱相匹配的时程记录是目前一致认可的选波原则之一. 本文研究现有几种主要的强震动记录的选取方法,扩展了太平洋地震工程研究中心的地面运动数据库(PGMD),同时根据日本MW9.0特大地震记录整理回归建立加速度反应谱衰减关系,并应用匹配条件均值反应谱法选择地震动.

    最早关于结构时程分析中选波的规定见1990SEAOC蓝皮书 Structural engineers association of California (SEAOC). 1990. Recommended lateral force requirements and commentary,‘Blue Book’ 5th edn. SEAOC,Sacramento,CA. : “时程分析应至少使用3组相似的水平分量,取其最大参数反应作为设计值”,“应用于场地的地面运动时程应代表实际地震地面运动,其反应谱应近似于场地设计谱”,“通过放缩使其均方根不低于1.3倍的DBE(或MCE),在T1-T1+2.0 s内大于10%”. 1994年的UBC(NEHRP条款) Building Seismic Safety Council. 1994. National earthquake hazards reduction program. Recommended provisions for the development of seismic regulations for new buildings. Washington D D. 规定: 如果用7组或更多记录,设计用其平均值. ATC-40 ATC (Applied Technology Council). 1996. Seismic evaluation and retrofit of concrete buildings,ATC-40. California Seismic Commission,Sacramento,CA. 规定,不小于3组水平分量,选择记录地震: M,R及震源机制相符,场地条件相符,当没有相似记录使用人工合成地面运动. ATC-63 Quantification of building seismic performance factors. April 2008. ATC-63 project report (90% Draft),FEMA P695. 报告基于8点原则,针对中硬场地,建议了22条远场地震波,28条近场地震波,来考虑地震波离散性的影响. 根据美国抗震规范选择时程,结构反应的中值估计可相差3倍(Cornell,2006).

    本文所指的目标谱是一种加速度反应谱,代表结构抗震设防需求及工程场地地震境况,包括规范设计谱、 一致概率反应谱及条件均值谱等. 各国抗震规范反应谱均采用理想化数学表达式形式,可由2—3个独立变量确定. 严格说来,任何谱的频谱(或周期)都是无穷多个独立变量,谱中幅值是频率(或周期)的连续函数,每一个频率(或周期)的谱值都是一个独立变量. 假若将谱简化为三直线(或曲线表达式),则表示谱有3个独立变量. 若进一步假定谱形简化为一个或两个变量,ATC-63中在0.1—0.5 s周期内加速度反应谱平均值和1 s周期的反应谱值时就只采用两个独立变量,Housner (1975)采用谱强度概念后就只有一个独立参数. 此时,反应谱与振幅在反映地震动强度上基本等价. 我国抗震规范中设计标准反应谱最长周期仅为6 s,远不能满足某些长周期结构的需要,也无法确定各周期点谱值的概率水平,而且规范中设计谱尚未考虑近断层地震作用的影响(耿淑伟,2005). 将其作为目标谱用于选择地面运动记录,很难准确预测结构在未来可能遭遇地震的真实反应特性.

    Cornell(2006)阐明一致概率反应谱(uniform hazard spectrum,简写为UHS)表示不同周期下加速度谱值具有相同概率,是若干个设定地震反应谱的包络线,而非某一个实际未来可能遭遇的地面运动记录的反应谱. 另外,一致概率反应谱与一条真实地震反应谱相比而言太过于广谱,直接作为目标谱用于结构时程分析选择地面运动,其结果将过于保守(王刚,2010).

    Baker和Cornell (20052006a, b)通过设定震源的震级和震中距等参数,构造满足概率结构需求分析的条件均值反应谱(conditional mean spectrum,简写为CMS). 条件均值反应谱是根据工程场地的设定地震同时考虑结构动力特性,利用加速度反应谱衰减关系确定的加速度反应谱. 条件均值反应谱选波方法基本思路: 基于概率结构需求分析和地震危险性分析建立条件均值反应谱为目标谱,选择与其匹配的地面运动加速度记录用于结构时程分析. 条件概率均值反应谱反映了在不同周期T下反应谱不完全相同的超越概率,即指定周期下反应谱值的超越概率最大而其它周期谱值超越概率介于最大与均值之间(王刚,2010). 大于或者小于均值反应谱的标准方差数值ε图1所示,根据概率结构需求分析或规定结构加速度反应值Sa(T1)确定ε(T1). 某个给定周期T1的反应谱值具有给定的ε(T1)个标准方差数,其它周期上的ε(Ti)则小于此ε(T1). 由概率地震危险性分析确定震级和震中距等地震参数的特征值,将此特征值代入加速度反应谱衰减关系,确定Sa(Ti)平均值.

    图  1  条件均值反应谱和一致概率反应谱示意图(1g=9.8 m/s2)
    Figure  1.  Conditional mean spectrum and uniform hazard spectrum (1g=9.8 m/s2)

    PGMD地面运动记录数据库包括了世界范围内位于地质活跃区的浅源地震,其中每条记录都含有各种距离参数、 场地类型及震源参数. 美国抗震安全委员会(BSSC)的“国家减轻地震灾害计划”(NEHRP)最新一版推荐规程将PGMD数据库作为地面运动记录的主要来源. 20世纪90年代,太平洋地震工程研究中心(PEER)迫于改善使用地震地面运动记录的需求,开始搭建强震记录的网络搜索数据库. PEER的新一代衰减关系(next generation attenuation,简写为NGA)数据库共收集整理了来自1 456个地震台的173次地震,共计3 551组三分量地震波加速度时程. 该数据库记录的矩震级分布为M4.27—7.9,震源距分布约为5—558 km,其中M≥8.0大地震地面运动记录一直是空白. 当同时对震级、 距离、 震源机制加以限制并规定每次地震只取1个记录时,仅有9条记录备选.

    本文在PEER地面运动数据库(PGMD)的基础上整理了近几年国内外特大、 大地震的地面运动记录,建立地面运动数据库. 其中包括: 2011年3月11日日本Tohoku MW9.0地震地面运动记录141组(温瑞智等,2011),2008年5月12日中国汶川MS8.0地震的129组主震记录,以及2010年9月3日新西兰MW7.1地震44组记录,共计314组实际地震地面运动记录. 日本MW9.0地震震中位于38.1°N、 142.9°E,震源深度24 km,台站场地30 m. 等效剪切波速分布为210—2 270 m/s,断层距分布为10—643 km范围内. 中国汶川MS8.0地震,震中位于30.986°N、 103.364°E,震源深度14 km,台站场地30 m. 等效剪 切波速为221—660 m/s,断层距的范围为10—552 km. 新西兰MW7.1地震,震源深度11 km,台站场地30 m. 等效剪切波速为128—760 m/s,断层距的范围为3.6—194 km. 本研究不但对PGMD数据进行了适时的扩展,更重要的是填补了目前PGMD数据库缺失大震记录的空白. 本数据库共3 865组地面运动记录,其参数分布情况见图2.

    图  2  强震记录参数分布图
    Figure  2.  Distribution of strong-motion parameters

    本研究采用Boore (2008)提出的地面运动记录数据处理方法及TSPP程序进行滤波等处理记录,使用记录频率段为0.1—25 Hz,阻尼比为5%的弹性加速度反应谱. 数据库中的所有地面运动记录都旋转成平行与垂直断层方向的两个水平分量. 本文采用两个水平方 向加速度反应谱的几何平均值作为计算选取记录的依据. K-NET数据台站分布情况见图3

    衰减关系是在大量记录统计回归分析后得到的地震动预测方程. 美国NGA 研究计划对全美的断层特性、 场地特性等进行大量研究,研究项目组用获得的5种浅层活跃地震动衰减模型来构造确定性的境况地震(Abrahamson,Silva,2008; Boore,Atkinson,2008; Campbell,Bozorgnia,2008; Chiou,Youngs,2008; Idriss,2008),建立不同地区不同设防水准对应的地面运动目标反应谱. 加速度反应谱的衰减关系是建立目标反应谱的前提. 5组回归模型所用数据量不同,AS08和CY08包括余震记录,所以数据量较大; I08也考虑余震记录,但仅限于基岩场地记录; AS08,CB08和CY08考虑断层的上盘效应、 断层深度. 5组模型地震参数各不相同,震源因素包括矩震级和震源机制,但I08没有区分走滑与正断层地震. AS08,CB08和CY08用地面到破裂顶的深度来考虑破裂深度效应. 在距离上,BA08用台站到破裂面投影的最短水平距离Rjb,其它4个用台站到破裂面最短距离Rrup. NGA采用全球范围内的构造浅层地震,用于美国西部场地的地震动预测,对于国内的数据库建立和地面运动记录选取有实际意义. 目前,国内各类规范中的衰减关系仅适用于基岩场地.

    图  3  K-NET台站分布图
    Figure  3.  Location of K-NET strong-motion stations

    PEER的新一代衰减关系(NGA)主要适用范围如下: 走滑地震M<8.5,逆断层地震M<8.0; 断层距小于200 km,距地表30 m的土层剪切波速为180—1 300 m/s,加速度谱周期为0.01—10 s. 由于篇幅所限,本文选择其中两个衰减关系由M7.6外推至M9.0,采用日本M9.0记录验证其适用性. 其中,BooreAtkinson(2008)则建议采用公式(1)中模型进行地震动参数的衰减分析(以下简称BA08衰减关系).

    式中,lnY为地震动参数值,即峰值加速度和结构自振周期为0.04—10 s内的加速度反应谱值,采用两水平分量的几何平均值; M为震级; Rjb为断层距; vS30为距地表30 m土层平均剪切波速. 假定地震为本次M9.0地震,震源机制为逆断层地震,两个场地的断层距分别为10和200 km,场地地下30 m剪切波速均为760 m/s. 由公式(1) BA08衰减关系预测的加速度反应谱结果如图4所示. 可见BA08衰减关系外推至M9.0地震时,当断层距较小,即近场地震,短周期的加速度反应谱的预测值偏小而长周期预测值偏大; 当断层距较大,即中远震,短周期的加速度反应谱的预测值偏大而长周期预测相对更接近实际情况下的均值.

    图  4  实际记录与BA08衰减关系对比. (a) R=10 km; (b) R=200 km
    Figure  4.  Comparison of magnitude 9.0 ground motion records in Japan with BA08 attenuation model of acceleration response spectrum. (a) R=10 km; (b) R=200 km

    CampbellBozorgnia (2008)采用以下衰减关系模型(以下简称CB08衰减关系),见公式(2). CB08衰减关系除了考虑震级、 断层距、 震源机制、 场地条件外,还考虑了上盘效应的影响.

    式中,fi表示震级M,断层距R,断裂类型F,上盘效应HW,浅层场地条件S,沉积或盆地深度D对地面运动的影响函数. 设定地震与实际记录谱选择方法同上,比较 CB08衰减关系得到的预测谱值和实际记录谱值. 即当设定地震震级M=9.0、 断层距为10 km时,实际记录及其平均值和预测值见图5a; 设定断层距为200 km时见图5b. 可见CB08衰减关系外推至M9.0地震时,当断层距较小,即近场地震,其结果与BA08衰减关系相近,短周期的加速度反应谱的预测值偏小于实际记录谱的平均值而长周期预测值偏大; 当断层距较大,即中远震,预测值整体偏大,长周期处偏离较为突出且远大于实际谱值.

    图  5  实际记录与CB08衰减关系对比.
    (a) R=10 km; (b) R=200 km
    Figure  5.  Comparison of magnitude 9.0 ground motion records in Japan with CB08 attenuation model of acceleration response spectrum.
    (a) R=10 km; (b) R=200 km

    以上分析表明,利用现有衰减关系预测M9.0地震加速度反应谱与实际记录加速度反应谱差距较大,特别在大震近场情况下两者谱形相差明显. 可见现有衰减关系并不可以外推应用于M9.0地震,也不能直接用其建立目标反应谱.

    本文仅考虑断层距和场地条件因素建立衰减关系,模型见公式(3). 利用本次M9.0地震获得141组两个水平方向的实际记录及其地震动相关参数,采用高斯-牛顿非线性最小二乘,拟合衰减关系模型的回归系数. 数据分布: 断层距从10—635 km,场地地下30 m平均剪切波速从210—2 270 m/s,地面运动记录峰值加速度PGA从0.001g—1.904 8g. 参考BA08加速度反应谱衰减关系,确定回归模型为式(3)(以下简称LL衰减关系).

    统计回归得到不同周期的回归系数(表1). 本衰减关系反映了不同周期加速度反应谱值关于断层距、 场地条件的衰减关系,可用其对M9.0逆断层地震的加速度反应谱进行预测(图6). 由于记录数量有限,未对其剪切波速进行限制,其中图6ab分别为设定VS30=760 m/s情况下,R=10和100 km时,应用公式(3)得到的预测加速度反应谱与相应条件的实际记录谱值. 可见,LL谱加速度衰减关系可以很好地拟合M9.0地震的加速度反应谱,可将其作为建立条件概率均值谱的依据.

    表  1  LL衰减关系回归系数
    Table  1.  Regression coefficient of LL attenuation model of acceleration response spectrum
    下载: 导出CSV 
    | 显示表格
    图  6  实际记录与LL衰减关系对比
    (a) R=10 km; (b) R=100 km
    Figure  6.  Comparison of magnitude 9.0 ground motion records in Japan with LL attenuation model of acceleration response spectrum
    (a) R=10 km; (b) R=200 km

    本文采取公式(3)LL衰减关系代替PGMD系统中NGA的衰减关系,建立M9.0地震的条件均值反应谱为目标谱. 为验证增加大震记录及调整衰减关系的必要性,本文例举如下算例: ① 若某场地经概率地震危险性分析分布式结果显示设定地震M9.0、 断层距为R=100 km的地震发生概率凸显,场地地下30 m平均剪切波速为760 m/s,拟建结构自振周期为T1=1 s,假定根据所需结构加速度谱Sa(T1)值确定谱形参数ε=2.0; ② 若场地的断层距为R=10 km,其它条件同①,确定目标反应谱,当设计地震为M9.0,断层距为R=100 km,以指定周期为2 s的加速度反应谱值的ε=2.0为已知,生成条件均值反应谱(CMS)如图7a所示. 当设计地震为MW9.0,断层距为R=10 km,谱形参数ε=2.0时,其条件均值反应谱为图7b所示.

    图  7  MW9.0条件均值反应谱
    Figure  7.  Magnitude 9.0 conditional mean spectrum

    本文使用匹配条件均值谱的均值和方差的理论及算法(Jayaram et al,2011). 以上述情况为例,将图7中条件均值谱作为目标谱,采用蒙特卡罗概率方法匹配目标谱的均值和方法,从补充大震记录后的扩展数据库中,选择10组地面运动记录匹配目标谱,选择结果及放缩系数见表 2. 按照记录放缩后与目标谱优先匹配程度,即两者均方差逐渐增大的顺序进行排序. 由表 2谱形匹配结果可见,在选择的10组记录中有3组记录震级M=9.0,且放缩系数分别为1.01,2.78和1.44,两者相关性见图8a. 表明对震级M=9.0记录的调幅较小,采用此记录能更多地反映实际大震特征. 图8bc所示分别为放缩系数和断层距及场地剪切波速的关系. 由图8def可知,记录的放缩系数(SF)与PGA,PGV和PGD都存在很好的相关性. 其中SF与PGA的决定系数R2为0.736,说明统计回归的趋势线与原数据拟合程度较好,即PGA的对数与SF存在线性相关性. 当M<8.0时,SF与PGV的决定系数R2为0.567; 当M=9.0时,SF与PGV的决定系数R2为0.999; 而不考虑震级影响时,两者的决定系数接近于零. SF与PGD的决定系数R2为0.911,即PGD的对数与SF存在显著线性相关性. 确定选波原则时,可考虑限定PGA,PGV和PGD的最小值而间接限制放缩系数的最大值. 但由于数据有限,不建议将回归公式直接用于记录的选取. 由此例可见,近场大震记录对于数据库在选波系统中的应用十分重要.

    表  2  扩展数据库选波结果
    Table  2.  The results of PGMD expanded 9.0 magnitude ground motion records
    下载: 导出CSV 
    | 显示表格
    图  8  选波结果中的放缩系数与地震动参数的相关性分析
    (a) 震级-放缩系数; (b) 断层距-放缩系数; (c) 剪切波速-放缩系数; (d) PGA-放缩系数; (e) PGV-放缩系数; (f) PGD-放缩系数
    Figure  8.  Logarithmic linear correlation between scalar factors and earthquake intensity parameters
    (a) Earthquake magnitude-scalar factor; (b) Fault distance-scalar factor; (c) vS30-scalar factor; (d) PGA-scalar factor; (e) PGV-scalar factor; (f) PGD-scalar factor

    在结构时程分析中地面运动记录的不确定性是导致结构反应参数离散性过大的主要原因之一. 本文在PEER地面运动数据库(PGMD)的框架及数据模式的基础上,结合近几年国内外特大地震地面运动记录建立地面运动数据库,同时根据日本MW9.0特大地震记录整理回归建立加速度反应谱衰减关系,并应用匹配条件均值反应谱法选择地面运动. 选波实例表明,当设定地震为特大地震时,基于条件均值反应谱法选取地面运动记录时,扩展数据库中大震记录并建立符合大震记录加速度反应谱的衰减关系是十分必要与迫切的.

    鉴于文中表2所列的选波结果,在保证谱形匹配的前提下,限制震级大于MW6.0、 放缩系数小于5,可得到7条地面运动记录. 在此基础上增加限制条件PGA大于0.3g,则可得3条地面运动记录用于结构时程分析. 因此,建议使用谱形匹配选波时,可将震级、 放缩系数、 PGA列入选波限制条件,优化选波结果.

    本研究实例为坚硬土层场地的大震,对于此类情况,建议使用震级大于MW6.0、 放缩 系数小于5,PGA大于0.3g作为限制条件,在此基础上匹配条件均值反应谱. 对于M5.0 —7.0设定地震选波原则可参考ATC-63报告. 基于8点原则,按照《抗震规范》规定最大地震影响系数,确定与Sa(T1)的比值进而确定ε值,生成条件均值反应谱作为目标谱,匹配目标谱.

    下一步研究的问题,即通过结构时程分析结果评价条件均值谱法及其它选波方法,为抗震规范结构时程分析选取地面运动记录提供科学依据及参考.

  • 图  1   蒙古板块及邻区构造图(据Tapponnier,Molnar(1979)修改)

    Figure  1.   Main tectonic structures in Mongolia and its adjacent areas(revised after Tapponnier,Molnar,1979)

    图  2   跨贝加尔裂谷台站各向异性结果与GPS观测地表形变结果(据Barruol等(2008)修改)

    黑线和蓝线表示 Barruol等(2008)Gao等(1994a)的横波分裂结果

    Figure  2.   Splitting results of the stations across Baikal rift and ground deformations referred from GPS data(after Barruol et al,2008)

    The black and blue lines represent the splitting results of Barruol et al(2008) and Gao et al (1994a),respectively. Green lines with arrow represent the ground deformation from GPS data, and the red lines with arrow represent the absolute plate velocity in this area

    图  3   CM台阵位置分布图

    Figure  3.   Location of the CM array

    图  4   震源分布图

    (a)震源随后方位角的分布统计图;(b)震源分布位置

    Figure  4.   Distribution of hypocenters used in this study

    (a)Statistic on distribution of earthquakes with back azimuth; (b)Location of earthquakes

    图  5   地震事件记录示例(记录起始时刻19:31:04)

    事件: 2011-09-15,台站: CM50,MW7.3,后方位角: 118.23°,震中距: 93.47°,震源深度: 644 km

    Figure  5.   Sample of the event records since 19:31:04

    Event: 2011-09-15,station: CM50,MW7.3,backazimuth: 118.23°,epicentral distance: 93.47°,focal depth: 644 km

    图  6   利用SplitLab软件得到的可信分裂结果: 以CM50台站2011-09-15地震事件为例

    (a)切向与径向波形;(b)分裂结果及后方位角;(c)、(g)快慢波波形拟合;(d)、(h)切向能量; (e)、(i)质点极化前后振动轨迹;(f)快慢波波形相关系数分布; (j)不同分裂参数的切向能量分布

    Figure  6.   Credible results got by SplitLab: Taking the event on 15 September 2011 recorded by the station CM50 as an example

    (a)Tangential and radial waveforms;(b)Splitting result and the back azimuth;(c)-(e)Fast-slow waveform fitting,rotated tangential energy and particle polarization by RC method;(f)Correlation coefficient of splitting parameters;(g)-(i)Fast-slow waveform fitting,rotated tangential energy and particle polarization by SC method;(j)Rotated tangential energy of splitting parameters

    图  7   利用SplitLab软件得到的可用分裂结果: 以CM17台站2011-12-11地震事件为例

    (a)切向及径向波形;(b)分裂结果及后方位角;(c)、(g)快慢波波形拟合;(d)、(h)切向能量; (e)、(i)质点极化前后振动轨迹;(f)快慢波波形相关系数分布布

    Figure  7.   Valuable results got by SplitLab: Taking the event on 11 December 2011 recorded by the station CM17 as an example

    (a)Tangential and radial waveforms;(b)Splitting result and the back azimuth;(c)-(e)Fast-slow waveform fitting,rotated tangential energy and particle polarization by RC method;(f)Correlation coefficient of splitting parameters;(g)-(i)Fast-slow waveform fitting,rotated tangential energy and particle polarization by SC method;(j)Rotated tangential energy of splitting parameters

    图  8   用3种方法取得的分裂结果对比图

    Red line represents the result by SC method,blue line denotes the result by RC method and yellow line denotes the result by EV method

    Figure  8.   Comparison of the results from three methods

    图  9   CM台阵取得的无效分裂结果

    Figure  9.   Null splitting results of CM array

    图  10   快波方向φ与快慢波延时δt统计结果

    Figure  10.   Statistical results of the fast-wave direction φ and delay time δt

    图  11   SKS分裂后方位角图

    Figure  11.   Back-azimuth of the SKS splitting results

    图  12   横波分裂结果的比较

    Red lines represent the splitting results of CM stations in this study; blue lines denote the results from Gao et al(1994a); yellow lines show the absoulte plate motion direction based on HS3-NUVEL-1A plate model; black lines represent the main tectonic line

    Figure  12.   Comparison of shear-wave splitting results

    表  1   利用最小能量法得到的蒙古中部地区SKS波分裂参数

    Table  1   Splitting parameters for SKS phase in central Mongolia region by SC method

    下载: 导出CSV
  • 常利军, 王椿镛, 丁志峰. 2006. 云南地区SKS波分裂研究[J]. 地球物理学报, 49(1): 197-204.

    Chang L J, Wang C Y, Ding Z F. 2006. A study on SKS splitting beneath the Yunnan region[J]. Chinese Journal of Geophysics, 49(1): 197-204 (in Chinese).

    常利军, 王椿镛, 丁志峰. 2009. 首都圈地区SKS波分裂研究[J]. 地震学报, 30(6): 551-559.

    Chang L J, Wang C Y, Ding Z F. 2009. A study on SKS splitting beneath capital area of China[J]. Acta Seismologica Sinica, 30(6): 551-559 (in Chinese).

    车自成, 刘良, 罗金海. 2002. 中国及其邻区区域大地构造学[M]. 北京: 科学出版社: 226-228.

    Che Z C, Liu L, Luo J H. 2002. Tectonics of China and Its Adjacent Areas[M]. Beijing: Science Press: 226-228 (in Chinese).

    王永锋, 金振民. 2005. 地震波各向异性: 窥测地球深部构造的 "探针"[J]. 地球科学进展, 20(9): 946-953.

    Wang Y F, Jin Z M. 2005. Seismic anisotropy: A probe to understand the structure in Earth's interior[J]. Advances in Earth Science, 20(9): 946-953 (in Chinese).

    吴萍萍, 王椿镛, 丁志峰, 常利军, 李振. 2012. 大别-苏鲁及邻区上地幔的各向异性[J]. 地球物理学报, 55(8): 2539-2550

    Wu P P, Wang C Y, Ding Z F, Chang L J, Li Z. 2012. Seismic anisotropy of upper mantle beneath the Dabie-Sulu and its adjacent areas[J]. Chinese Journal of Geophysics, 55(8): 2539-2550 (in Chinese).

    张萍. 2011. 基于Bootstrap方法的统计分析[J]. 宜宾学院学报, 11(12): 31-33.

    Zhang P. 2011. Statistical analysis based on Bootstrap method[J]. Journal of Yibin University, 11(12): 31-33 (in Chinese).

    张建利, 田小波, 张洪双, 司少坤, 张中杰, 滕吉文. 2012. 贝加尔裂谷区地壳上地幔复杂的各向异性及其动力学意义[J]. 地球物理学报, 55(8): 2523-2538.

    Zhang J L, Tian X B, Zhang H S, Si S K, Zhang Z J, Teng J W. 2012. The crust and upper mantle anisotropy in Baikal Rift Zone and its dynamic significance[J]. Chinese Journal of Geophysics, 55(8): 2523-2538 (in Chinese).

    Barruol G, Deschamps A, Déverchère J, Mordvinova V V, Ulziibat M, Perrot J, Artemiev A A, Dugarmaa T, Bokelmann G H R. 2008. Upper mantle flow beneath and around the Hangay dome, Central Mongolia[J]. Earth Planet Sci Lett, 274(1/2): 221-233.

    Bowman J R, Ando M. 1987. Shear-wave splitting in the upper-mantle wedge above the Tonga subduction zone[J]. Geophys J Int, 88(1): 25-41.

    Efron B. 1979. Bootstrap methods: Another look at the jackknife[J]. The Annals of Statistics, 7(1): 1-26.

    Gao S, Davis P M, Liu H, Slack P D, Zorin Y A, Mordvinova V V, Kozhevnikov V M, Meyer R P. 1994a. Seismic anisotropy and mantle flow beneath the Baikal rift zone[J]. Nature, 371(6493): 149-151.

    Gao S, Davis P M, Liu H, Slack P D, Zorin Y A, Logatchev N A, Kogan M, Burkholder P, Meyer R P. 1994b. Asymmetric unwarp of the asthenosphere beneath the Baikal rift zone, Siberia[J]. J Geophys Res, 99(B8): 319-330.

    Gripp A E, Gordon R G. 2002. Young tracks of hotspots and current plate velocities[J]. Geophys J Int, 150(2): 321-361.

    Li Y H, Wu Q J, Zhang F X, Zhang R Q. 2011. Seismic anisotropy of the Northeastern Tibetan Plateau from shear wave splitting analysis[J]. Earth Planet Sci Lett, 304(1/2): 147-157.

    Silver P G, Chan W W. 1988. Implications for continental structure and evolution from seismic anisotropy[J]. Nature, 335(6185): 34-39.

    Silver P G, Chan W W. 1991. Shear wave splitting and subcontinental mantle deformation[J]. J Geophys Res, 96(16): 16429-16454.

    Tapponnier P, Molnar P. 1979. Active faulting and Cenozoic tectonics of the Tien Shan, Mongolia, and Baykal regions[J]. J Geophys Res, 84(B7): 3425-3459.

    Vauchez A, Tommasi A, Barruol G, Maumus J. 2000. Upper mantle deformation and seismic anisotropy in continental rifts[J]. Phys Chem Earth PT A, 25(2): 111-117.

    Wüstefeld A, Bokelmann G. 2007. Null detection in shear-wave splitting measurements[J]. Bull Seismol Soc Am, 97(4): 1204-1211.

    Wüstefeld A, Bokelmann G, Zaroli C, Barruol G. 2008. SplitLab: A shear-wave splitting environment in Matlab[J]. Computers & Geosciences, 34(5): 515-528.

    Zorin Y A, Kozhevnikov V M, Novoselova M R, Turutanov E K. 1989. Thickness of the lithosphere beneath the Baikal rift zone and adjacent regions[J]. Tectonophysics, 168(4): 327-337.

  • 期刊类型引用(24)

    1. 朱音杰,罗艳,赵里. 区域宽频地震波形反演所揭示的2022年青海门源M_S6.9地震震源破裂过程. 地震学报. 2023(05): 781-796 . 本站查看
    2. 罗毅,胡应顺,田云锋,张景发. 2014年新疆于田M_S7.3地震的发震构造及余震特征. 地震学报. 2021(01): 48-56+136 . 本站查看
    3. 董培育,石耀霖,程惠红,乔学军. 青藏高原及邻区未来地震活动性趋势数值分析. 地球物理学报. 2020(03): 1155-1169 . 百度学术
    4. 冉慧敏,上官文明,赵晓成. 2020年新疆于田M_S6.4地震及余震序列定位研究. 内陆地震. 2020(03): 249-257 . 百度学术
    5. 梁洪宝,武艳强,陈长云,熊威. 2014年新疆于田M_S7.3地震同震位移及位错反演研究. 地球物理学报. 2018(12): 4817-4826 . 百度学术
    6. 王琼,解朝娣,冀战波,刘建明. 2014年于田M_S7.3地震对后续余震和远场小震活动的动态应力触发. 地球物理学报. 2016(04): 1383-1393 . 百度学术
    7. 唐明帅,王海涛,李艳永,热依木江,孔祥艳,魏芸芸,上官文明,魏斌. 2014年新疆于田M_S7.3地震序列重定位及其发震构造初步研究. 地球物理学报. 2016(06): 2126-2137 . 百度学术
    8. 刘建明,吴传勇,王琼,孔祥艳,陈向军. 2011—2014年3次于田M_S≥5.0地震序列重定位及其发震构造. 地震学报. 2016(05): 671-683+813 . 本站查看
    9. 刘军,孙甲宁,宋立军,毕雪梅,谭明,何金刚. 2014年新疆于田M_S7.3地震强地面运动模拟及震动图预测. 地震学报. 2016(06): 914-923+960 . 本站查看
    10. 徐彦,李丹宁,冉慧敏. 2014年2月12日新疆于田M_W 6.9级地震破裂方向研究. 地球物理学进展. 2015(06): 2574-2580 . 百度学术
    11. 房立华,吴建平,王未来,杨婷,王长在. 2014年新疆于田MS7.3级地震序列重定位. 地球物理学报. 2015(03): 802-808 . 百度学术
    12. 李玉江,陈连旺,杨树新,刘少峰,杨兴悦. 两次于田M_S7.3地震间应力触发作用及2014年于田地震的发生对周缘断层的影响. 地球学报. 2015(01): 94-102 . 百度学术
    13. 李海兵,潘家伟,孙知明,刘栋梁,张佳佳,李成龙,刘亢,CHEVALIER Marie-Luce,云锟,龚正. 2014年于田Ms7.3地震地表破裂特征及其发震构造. 地质学报. 2015(01): 180-194 . 百度学术
    14. Chen Yanan,Jiang Haikun. On Estimating Magnitude of a Maximum Sequent Earthquake by Viscoelastic Coulomb Stress Change and a Discussion of the Relationship between the M_S7.3 Earthquakes in Yutian 2008 and 2014. Earthquake Research in China. 2015(04): 434-451 . 必应学术
    15. Meng Lingyuan,Zhou Longquan,Shi Haixia. A Study on the Characteristics of the Yutian, Xinjiang M_S7.3 Earthquake, February 12, 2014. Earthquake Research in China. 2015(01): 47-56 . 必应学术
    16. 陈亚男,蒋海昆. 基于粘弹库仑应力变化的后续最大地震震级估计及2008、2014年于田2次7.3级地震之间关系的讨论. 中国地震. 2015(01): 11-26 . 百度学术
    17. 徐伟民,陈石,石磊. 新疆于田及周边地区地震活动性与重力异常特征. 地球科学(中国地质大学学报). 2014(12): 1831-1841 . 百度学术
    18. 张广伟,雷建设,孙长青. 2014年2月12日新疆于田M_S7.3级地震主震及余震序列重定位研究. 地球物理学报. 2014(03): 1012-1020 . 百度学术
    19. 李海兵,孙知明,潘家伟,刘栋梁,张佳佳,李成龙,刘亢,云锟,龚正. 2014年于田Ms7.3地震野外调查——特殊的地表破裂带. 地球学报. 2014(03): 391-394 . 百度学术
    20. 孟令媛,周龙泉,刘杰. 2014年新疆于田M_S7.3地震近断层强地面运动模拟及烈度分布估计. 地震学报. 2014(03): 362-371 . 本站查看
    21. 程佳,杨文,刘杰,周龙泉. 2014年2月12日新疆于田M_S7.3地震序列及其构造背景研究. 地震学报. 2014(03): 350-361 . 本站查看
    22. 孟令媛,周龙泉,张小涛,杨文,姚琪,宋金. 2014年新疆于田M_S7.3地震序列和震源特征初步分析. 中国地震. 2014(02): 159-167 . 百度学术
    23. 宋金,周龙泉. 2014年于田M_S7.3地震产生的静态库仑应力变化及对周边断层的影响. 中国地震. 2014(02): 168-177 . 百度学术
    24. 陈为涛,王阅兵,连尉平,甘卫军,付广裕. 2014年于田M_s7.3级地震近场地表运动和同震位移的初步分析结果. 震灾防御技术. 2014(04): 838-846 . 百度学术

    其他类型引用(9)

图(12)  /  表(1)
计量
  • 文章访问数:  579
  • HTML全文浏览量:  306
  • PDF下载量:  17
  • 被引次数: 33
出版历程
  • 收稿日期:  2013-06-25
  • 修回日期:  2013-08-08
  • 发布日期:  2014-06-30

目录

/

返回文章
返回