相对局部区域震源参数随机不确定性经验关系研究

李宗超, 高孟潭, 孙吉泽, 陈学良, 张波

李宗超,高孟潭,孙吉泽,陈学良,张波. 2021. 相对局部区域震源参数随机不确定性经验关系研究. 地震学报,43(4):483−497. DOI: 10.11939/jass.20200153
引用本文: 李宗超,高孟潭,孙吉泽,陈学良,张波. 2021. 相对局部区域震源参数随机不确定性经验关系研究. 地震学报,43(4):483−497. DOI: 10.11939/jass.20200153
Li Z C,Gao M T,Sun J Z,Chen X L,Zhang B. 2021. Empirical relationship of stochastic uncertainty of source parameters in relative local area. Acta Seismologica Sinica43(4):483−497. DOI: 10.11939/jass.20200153
Citation: Li Z C,Gao M T,Sun J Z,Chen X L,Zhang B. 2021. Empirical relationship of stochastic uncertainty of source parameters in relative local area. Acta Seismologica Sinica43(4):483−497. DOI: 10.11939/jass.20200153

相对局部区域震源参数随机不确定性经验关系研究

基金项目: 国家重点研发计划(2019YFC1509403)和中国地震局地球物理研究所基本科研业务重大研究计划专题(DQJB19A0131,DQJB19A0133)联合资助
详细信息
    通讯作者:

    李宗超: e-mail:lizongchaoigo@163.com

  • 中图分类号: P315.08

Empirical relationship of stochastic uncertainty of source parameters in relative local area

  • 摘要: 在进行未来破坏性地震的强地面运动数值模拟时,震源参数选取的准确性对地震动预测的结果影响很大。震源参数的确定存在很多不确定性因素,既包含随机的不确定性因素,又包含认知的不确定性因素。本文在大量地震事件及文献调研的基础上,运用统计学方法对具备随机不确定性特征的震源参数进行统计研究,以震源参数经验公式的形态建立解释其随机性和不确定性的数学模型。为了研究局部地区震源参数的定标关系特征,获得更加适用于局部地震密集区域,尤其是包含中国大陆地区在内的局部区域的震源参数的经验关系,本文从GCMT地震目录中选取了1 700多个MW≥5.5的地震事件,运用统计学方法研究地震密集地区的震源参数经验关系,包括震级、地震矩、破裂面积等,增加了相对较大的局部范围内凹凸体的地震样本数量,从统计学角度计算更加适合局部区域的震源参数的经验关系。统计结果表明:局部区域震例获得的震源参数的经验关系与不限区域震例获得的经验关系存在差异,尤其是涉及到断层破裂面积、凹凸体相关参数时差异较大,局部区域内震例获得的震源参数的经验关系将更具有代表性。应用本文获得的相对局部区域的经验公式计算未来破坏性地震的强地面运动所需的震源参数时,获得的地震动预测结果将更能体现目标区域真实的地震动特征,进而提高地震动预测结果的可靠性。
    Abstract: In the numerical simulation of strong ground motion of future destructive earthquakes, the accuracy of source parameters selection has a great impact on the results of ground motion prediction. There are many uncertain factors in determining source parameters, including both random and cognitive uncertainties. Based on a large number of seismic events and literature researches, this paper focuses on statistical analysis of source parameters with random uncertainty characteristics by using statistical methods. Through regression analysis, a mathematical model is established to explain the randomness and uncertainty of source parameters in the form of empirical formula. In order to study the scaling relation characteristics of source parameters in local regions, we get more empirical relations which are more suitable for local seismic densely regions, especially those of the local regions including the Chinese mainland. This paper more than 1 700 seismic events with MW≥5.5 are selected from the global CMT catalogue. The empirical relationship of source parameters in earthquake intensive areas is studied by using statistical methods, including focal depth, magnitude, seismic moment, rupture area, etc. The number of seismic samples of asperity in a relatively large local range is increased, so as to obtain more suitable experience for local areas to calculate source parameters from the perspective of statistics relationship. The statistical results show that there are differences between the empirical relationship of source parameters obtained from local earthquake cases and those obtained from unlimited regional cases, especially when it comes to fault rupture area and asperity related parameters. The empirical relationship of source parameters obtained from local earthquake cases is more representative. When using the empirical formula obtained in this paper to calculate the focal parameters required for the strong ground motion of future destructive earthquakes, the ground motion prediction results will better reflect the real ground motion characteristics of the target area. It could improve the reliability of the ground motion prediction results.
  • 青藏高原地质构造复杂,地质、大地测量及地震构造等研究普遍表明,青藏高原内地壳变形分布广泛:高原内中部地壳主要沿着汇聚方向东西缩短,南部地壳向东西方向伸展,中北部及东部地壳向东部挤出,与云南和缅甸接壤的东南缘地壳向东南或正南方向挤出并在喜马拉雅东构造结处发生顺时针旋转(Molnar,Tapponnier,1975Holt et al,1991Wang et al,2001Zhang et al,2004)。青藏高原内分布有大量断裂带,构造活动丰富,地震频发(Yin,Harrison,2000邓起东等,2014叶卓等,2018)。地震层析成像和大地电磁探测等成果均显示青藏高原中下地壳存在低速区(Yang et al,2012Jiang et al,2014Bao et al,2015)和高电导率区(金胜等,2010Bai et al,2010),这表明中下地壳可能处于部分熔融状态(Chung et al,2003Le Pape et al,2012Long et al,2015)。作为青藏高原主体的中东部地区,从北到南依次分布着柴达木地块、巴颜喀拉地块、羌塘地块和拉萨地块,是研究青藏高原形成与演化的重要区域,也是研究高原中下地壳低速区成因和分布特征的重要区域。研究该区域地壳和上地幔顶部的速度结构和物性特征,可以推断青藏高原的物质迁移和地壳形变机制,对于理解大陆碰撞机制和生长过程具有重要作用,因此研究青藏高原中东部地壳和上地幔顶部的速度结构特征具有十分重要的意义。此外,青藏高原中东部的地震观测数据丰富,地震事件分布相对均匀,为层析成像提供了有利条件。

    近年来,除了利用GPS (Chen et al,2000Allmendinger et al,2007Hao et al,2016)和InSAR (李鹏等,2013Garthwaite et al,2013Zhao et al,2016)等方法对青藏高原中东部进行地表观测研究之外,研究人员对青藏高原乃至整个中国大陆进行了大量的地震波层析成像工作,例如,Pn波层析成像(Liang,Song,2006黎源,雷建设,2012)、S波层析成像(Mechie et al,2012黄忠贤等,2013)、远震或近震P波层析成像(Zhang et al,2012Lei,Zhao,2016)、背景噪声及地震面波层析成像(Yao et al,2006Ceylan et al,2012范文渊等,2015Pandey et al,2015钟世军等,2017)和全波形层析成像(Chen et al,20152017),取得了很多有意义的研究成果,其中近震P波走时成像仍是地震层析成像的重要手段。传统的近震P波层析成像方法主要是采用地壳P波和上地幔顶部Pn波分别反演的策略,这样可以避免莫霍面对地壳速度结构的影响,但是却容易造成部分高质量数据的浪费。尤其在Pn波层析成像中,忽略地壳横向非均匀性的影响可能会对构造复杂地区(青藏高原等)的反演结果带来较大误差(Li et al,2017),而使用Pn和Pg震相走时联合反演的方法能减小这种误差,取得较好的成像结果(Xu,Song,2010Huang et al,2013)。为此,本文提取研究区域内1991—2017年地震事件的初至波走时信息,拟利用Pg和Pn波震相走时联合反演方法,对青藏高原中东部的三维P波速度结构进行层析成像,以期得到该地区的地壳P波速度结构、上地幔顶部Pn波速度结构及地壳厚度信息。

    我们从美国地震联合会(Incorporated Research Institutions for Seismology,简称IRIS)数据中心收集了青藏高原中东部(88°E—100°E,25°N—40°N)从1991年到2017年共6 617个地震事件的波形记录。为保证数据质量,进行数据筛选,标准如下:① 震级ML≥2.5;② 震源深度d<60 km;③ 震中距≤12°。图 1给出本文所使用的地震台站分布及相应的台网代码和观测周期。

    图  1  研究区地质概况及各台网台站分布
    Figure  1.  Topography of the study region and distribution of stations of seismic networks

    采用Sacpick自动拾取方法,基于地震波形数据对初至波走时进行拾取,具体流程如下:① 根据预设模型计算生成初至波理论走时t;② 根据理论走时选取一定的噪声窗口和信号窗口,本文选取噪声窗口为[t-40 s,t-10 s],信号窗口为[t-20 s,t+10 s];③ 通过滑动窗口寻找振幅发生变化的信号并计算信号窗口与噪声窗口的均方根(root mean square,缩写为RMS)比值,当两者之比大于某一预设值(本文预设值取为5)时,认为该信号有效并予以采用。

    本文通过对245个台站记录的1 776个地震事件的地震波形进行自动拾取,最终得到1万4 484个初至波到时数据。图 2给出了2008年12月18日某地震事件初至波自动拾取的示例,可以看出,采用自动拾取走时方法能够将信噪比差的数据直接剔除,保留高质量的数据,并且有效地减少走时拾取中的人为因素,使拾取的走时具有较好的一致性。由于拾取到的到时数据为初至波走时,并非具体的Pg或Pn震相走时,因此在层析成像的每次迭代反演中,还需要分别计算Pg波和Pn波理论走时,并将其与观测走时进行对比,实现初至波震相的判别。

    图  2  2008年12月18日ML5.1地震事件的初至波走时拾取效果
    该事件的震中位置为(96.47°E,36.78°N),震源深度为0 km,相对走时为每条地震记录相对于初至波的走时
    Figure  2.  First arrival time picks from an ML5.1 earthquake on December 18,2008
    The epicentre of the event is (96.47°E,36.78°N) with a focal depth of 0 km. Travel times are relative to first arrivals

    采用Xu和Song (2010)提出的联合Pg和Pn震相走时反演方法,对Pg波和Pn波走时采用不同的计算方法。在球面坐标系下,利用伪弯曲射线追踪方法(Koketsu,Sekine,1998),Pg波从震源j到台站k的走时残差为

    ${r_{jk}} = t_{jk}^{{\rm{obs}}} - t_{jk}^{{\rm{pred}}} = \sum\limits_{i = 1}^n {\left({ - \frac{{{\rm{d}}l}}{v}} \right){{\left({\frac{{\Delta v}}{v}} \right)}_i}}{\text{,}}$

    (1)

    This page contains the following errors:

    error on line 1 at column 1: Start tag expected, '<' not found

    Below is a rendering of the page up to the first error.

    Pn波从震源j到台站k的走时残差为

    ${r_{jk}} {\text{=}} t_{jk}^{{\rm{obs}}} {\text{-}} t_{jk}^{{\rm{pred}}} {\text{=}} \sum\limits_{{{m}} {\text{=}} 1}^{{n_1}} {\left({ {\text{-}} \frac{{{\rm{d}}l}}{v}} \right){{\left({\frac{{\Delta v}}{v}} \right)}_{{m}}} {\text{+}} } \sum\limits_{{{e}} {\text{=}} 1}^{{n_2}} {\left({ {\text{-}} \frac{{{\rm{d}}l}}{v}} \right){{\left({\frac{{\Delta v}}{v}} \right)}_{{e}}} {\text{+}} \sum\limits_{{{r}} {\text{=}} 1}^{{n_3}} {\left({ {\text{-}} \frac{{{\rm{d}}l}}{v}} \right){{\left({\frac{{\Delta v}}{v}} \right)}_{{r}}} {\text{+}} } \eta \Delta {h_k} {\text{+}} \Delta {t_j}}{\text{,}} $

    (2)

    式中,等号右边第一项为射线沿莫霍面滑行的走时,第二项为震源端地壳部分走时,第三项为台站端地壳部分走时,第四项ηhk为台站端莫霍面深度变化引起的走时变化,其中

    $\eta = \sqrt {{{\left({\dfrac{1}{{{v_{\rm{c}}}}}} \right)}^2} - {{\left({\dfrac{1}{{{v_{\rm{m}}}}}} \right)}^2}}$

    (3)

    为垂直慢度,vcvm分别为该处地壳Pg波速度和Pn波速度,第五项Δtj为震源延迟项。

    This page contains the following errors:

    error on line 1 at column 1: Start tag expected, '<' not found

    Below is a rendering of the page up to the first error.

    由于本文研究区域与Xu和Song (2010)一文研究区域的地理位置比较靠近且东部有部分重叠,为减少反演迭代次数,在计算一维模型时参考其一维P波速度模型,将该模型作为初始地壳速度模型。通过走时数据拟合,得到Pn波平均速度约为8.08 km/s,地壳平均厚度约为55 km。将一维速度模型中地壳的最后一层延伸至55 km深度,得到初始一维速度模型,地壳模型均匀分层,每层厚度为10 km,最后一层划分至55 km (即莫霍面深度处)。采用Pg波和Pn波震相走时联合反演方法确定最终一维速度模型,经过7次迭代反演后残差不再明显降低,此模型即为后续层析成像所使用的一维速度模型(图3)。

    图  3  一维模型及其初始模型
    Figure  3.  One-dimensional model and its initial model

    采用计算得到的一维速度模型,利用所有初至波走时数据进行三维联合反演。反演采用的地壳速度水平网格、Pn波速度水平网格和莫霍面深度网格的大小均为0.5°×0.5°。联合反演中,震相判别时Pg波和Pn波两者理论走时与观测走时之差的绝对值均小于1 s的数据将不予采用,以减少由于模型粗糙所造成的震相误判,且在计算中要求参与反演的射线走时残差不大于3 s,以防止少量数据对模型造成过大影响。每次迭代结束后,更新模型进行下一次迭代反演。

    本次联合反演的初始数据共1万4 484个初至波走时,最终参与反演的射线有1万零444条,射线基本覆盖了研究区域(图4)。

    图  4  联合反演射线覆盖
    Figure  4.  Ray coverage of joint inversion

    经过8次迭代后,残差不再明显降低,得到联合反演的最终结果,反演后走时残差的分布有一定改善(图5a),但并不十分显著,说明本文采用的初始模型与本区域的真实速度结构比较接近,即Xu和Song (2010)的地壳模型对于本文研究具有较好的参考意义。除震中距大于10°的区间射线数量不足外,反演后各震中距区间的走时残差标准差(图5b)均在一定程度上有所降低,反演后走时残差分布更集中(图5c),走时残差的RMS由初始模型的1.50 s降低至1.21 s,同时参与反演的射线折合走时分布较为均匀(图5d),模型拟合得较好。

    图  5  联合反演的残差分析和走时分布
    (a) 反演前后走时残差随震中距的分布;(b) 反演前后各震中距范围内走时残差的标准差分布;(c) 联合反演前后走时残差数量分布及其高斯分布拟合曲线;(d) 折合走时(td/0.065)分布
    Figure  5.  Residual analysis and travel time distribution of joint inversion
    (a) Travel time residual distribution before and after inversion;(b) Standard deviations of travel time residuals with different epicentral distance before and after inversion;(c) Histogram of residuals and its Gaussian distribution fitting before and after inversion;(d) Reduced time by 0.065 degree per second

    由于DLSQR方法求得的结果无法得到方程解的误差矩阵,因此采用检测板测试(Lévěque et al,1993)对反演结果的可靠性进行分析。本文采用与实际数据反演相同的参数设置,对速度模型添加±4%的扰动,并通过反演进行恢复成像。通过多次试验,在网格大小为1.5°×1.5°时速度恢复情况较好,而边缘部分地区由于射线覆盖不足,无法得到理想的恢复效果(图6)。

    图  6  1.5°×1.5°网格下的检测板检测
    Figure  6.  Checkerboard test with squares of 1.5 degree by degree

    本文获得了青藏高原中东部地区地壳和上地幔顶部的P波三维速度结构及莫霍面深度信息(图7)。反演结果显示:地壳P波速度处于5.2—7.2 km/s之间,地壳同一深度处,P波速度变化范围较大,表明研究区域存在强烈的横向不均匀性;不同深度P波高、低速分布情况存在很大差异,20 km深度以上地壳部分与30 km深度以下地壳部分的P波高、低速分布也呈现出明显的纵向非均匀性;从地壳30 km深度处开始存在大片明显的低速区并一直延伸至莫霍界面;上地幔顶部Pn波速度存在较大的横向不均匀性,速度值主要分布在7.7—8.4 km/s之间,研究区域内Pn波速度在南部的拉萨地块和北部的柴达木盆地主要表现为高速异常,中间的羌塘地块西部和巴颜喀拉地块部分主要表现为低速异常,具有较大的横向差异;莫霍面深度约在48.0—68.6 km,主要表现为南深北浅(图8)。

    图  7  地壳Pg波速度结构和上地幔顶部Pn波速度结构
    各子图右下角给出了相应深度的P波平均速度,速度扰动相对于该层平均速度
    Figure  7.  Velocity structure of Pg and Pn waves
    The velocity perturbation in each layer is related to the average velocity given in the lower-right corner of the sub-figures
    图  8  青藏高原中东部的莫霍面深度分布
    Figure  8.  Distribution of Moho depth of the central-eastern Tibetan Plateau

    由于研究仅聚焦于青藏高原中东部,本次反演结果显示的P波速度结构与Huang等(2014)得出的三维P波速度结构整体上具有一致性,但结果更精细、速度扰动程度更大,这种差异可能来自于数据和区域特征所确定的反演参数。如图7所示,上地壳(0—20 km)与中下地壳(30—50 km) 的高、低速分布存在明显差异,说明青藏高原中东部地区的上、下地壳结构存在显著不一致性。地壳0—20 km深度内高速异常主要分布在中部的巴颜喀拉地块和羌塘地块,低速异常主要分布在南北两端的拉萨地块和柴达木地块;当地壳深度大于30 km时,低速异常主要分布于羌塘地块、拉萨地块东部,而高速异常主要分布于柴达木地块和拉萨地块西部及其以南部分。柴达木地块在上地壳主要表现为低速异常,速度值低于5.5 km/s,30 km以下深度则主要表现为高速异常,速度值高于7.0 km/s。羌塘地块和巴颜喀拉地块西部在地壳0—20 km深度表现为高速异常,30 km以下则表现为低速异常。此外,巴颜喀拉东部在不同深度处均表现为高速异常,拉萨地块西部出现明显高速异常,东部则主要表现为低速异常。

    Pn波速度分布(图7)存在显著的空间差异,此结果与前人的成像结果(Liang et al,2004王海洋等,2013Li,Song,2018)基本一致,但细节更加明显,显示出青藏高原中东部Pn波速度强烈的横向不均匀性。柴达木地块和拉萨地块上地幔顶部总体上表现为较高的Pn波速度,最高可达8.3 km/s。巴颜喀拉地块、羌塘地块总体上表现为低速,最低约为7.7 km/s。Liang和Song (2006)的研究结果是青藏高原Pn波层析成像的重要工作,本次成像得到的Pn波速度区间较Liang和Song (2006)对青藏高原成像的结果更大,高、低速异常更明显,但总体上来说速度异常的变化趋势与其一致,主要表现为羌塘地块西部和巴颜喀拉地块Pn波速度显示为低速异常,柴达木盆地和拉萨地块的Pn波速度主要显示为高速异常。

    由于台站空间分布的限制,本文反演得到的地壳厚度信息具有一定的局限性,但仍能反映出青藏高原中东部的地壳厚度变化。本文结果显示地壳厚度处于48.0—68.6 km之间(图8),莫霍面深度分布与前人研究(Shin et al,2007高锐等,2009)也具有很好的一致性。柴达木地块和巴颜喀拉地块东缘的地壳厚度小于50 km,相对较薄,巴颜喀拉地块中部、羌塘地块、拉萨地块和青藏高原东南部的地壳厚度较大,均在55 km以上,其中青藏高原东南部地壳厚度值基本在65 km以上。

    总体说来,地壳P波速度结构显示出的结果与地质构造具有较好的对应关系。由于浅部沉积物的存在,相对于波速较快的青藏高原变质岩山体,柴达木盆地P波速度在上地壳(20 km及以上)主要表现为低速异常,在下地壳(30 km及以下)主要表现为高速异常,这一结果与Li等(2012)范文渊等(2015)利用背景噪声得出的瑞雷波群速度表现出相似的特征,这表明盆地深部处于低温、稳定的状态。Yang等(20102012)对整个青藏高原中地壳进行了面波成像,其结果显示青藏高原中地壳存在大片低速区,本次研究结果显示研究区域内上、下地壳P波速度异常存在明显差异,下地壳存在大范围低速区,而结合上地幔顶部Pn波速度存在的显著低速异常,推测形成这种速度异常的原因可能是上地幔岩浆的侵入或加热作用(Godin et al,2006Grujic,2006Harris,2007)。Bai等(2010)利用大地电磁方法获得的结果表明,拉萨地块和羌塘地块北部中下地壳存在大片的高电导率区,本文得到的低速分布与其具有一致性,而本文结果显示的羌塘地块南部的高速异常与其得到的低电导率存在较大差异。该文指出大地电磁数据对深度电导率积分敏感,并非对特定深度的电导率值敏感,结合此区域上地幔顶部Pn波速度表现出的高速异常,推测造成这种差别的原因可能是羌塘地块南部上地幔顶部的岩性差异对该地区地壳的电导率结果产生了约束作用。

    此外,Zhao等(2013)利用地震波衰减的研究结果表明青藏高原可能存在两条地壳流通道,一条从羌塘地块自北向东沿四川盆地西缘流向青藏高原东南部,另一条从青藏高原南部的拉萨地块穿过喜马拉雅东构造结,总体上本文结果显示的下地壳P波低速区与该研究对应的羌塘地块和拉萨地块东部等地壳流具有较好的对应关系,而柴达木盆地、羌塘地块西部和巴颜喀拉地块东部等高速区为地壳流非对应区域,Bao等(2011)的研究结果同样表明这些区域的Q值显著高于其它区域,重力异常的密度成像结果也显示这些区域存在明显的高密度异常(杨文采等,2015),表明这三个区域中下地壳可能存在比较稳定的地体。

    本文结果显示青藏高原南缘具有显著连续的Pn波高速异常,而羌塘地块西部和巴颜喀拉地块主要表现为Pn波低速异常,这一结果与前人研究(Zhou,Lei,2016et al,2017)具有一致性,但本文研究所得到的Pn波波速异常更为明显,也存在一些明显的细节差异。Lei等(2014)关于青藏高原东部的层析成像结果揭示了一条连通柴达木盆地与拉萨地块的南北向高速带,而根据本文结果,尽管这部分区域总体上确实表现为高速异常,但并不连续,这可能是由于我们所采用的走时数据自动拾取策略和联合不同震相的反演方法具有一定的优势,且本次研究仅对较小区域进行成像,能得到更精细的成像结果。拉萨地块以南连续的Pn波高速异常说明这部分区域的上地幔顶部岩石圈具有连续、稳定和冷的特性,可能对应了北向俯冲的印度板块,可以说明印度板块对欧亚板块的俯冲块具有较好的整体性,且其俯冲角度可能较为平缓。而在羌塘地块西部和巴颜喀拉地块呈现的Pn波低速异常,则说明这些区域上地幔顶部的温度较高,表明这些区域可能发生了岩石圈的破裂。

    通常地震波的速度异常与大地热流值呈负相关关系(Black,Braile,1982汪素云等,2003),巴颜喀拉地块、羌塘地块表现为显著低速,即该区域可能有较高的热流值,这可能与上地幔顶部的岩浆作用紧密相关,从而导致了高温异常。由于印度板块对欧亚板块的向北俯冲,拉萨地块以南显示出的高速异常可能对应了温度较低的印度板块岩石圈部分。羌塘地块西部、柴达木地块和巴颜喀拉地块东部存在的稳定区域的Pn波速度稍大于平均速度,不同于其周围的显著低速,结合地壳P波速度的高速异常,可以推断这三个稳定地体可能持续延伸至上地幔顶部。总体来说,下地壳P波速度与上地幔顶部Pn波速度结构具有耦合作用,这说明高原中部的上地幔顶部和下地壳大范围的低速与岩浆作用密不可分,可能是由于上地幔顶部岩浆的加热作用(Owens,Zandt,1997),或是在印度板块俯冲作用下岩石圈发生破裂,导致上地幔岩浆上涌至下地壳形成大范围的高温异常(Nelson et al,1996)。

    由于地震震源深度存在不确定性,传统的Pn波层析成像通常采用台站延迟反映地壳厚度变化。本次联合反演可以直接反演地壳厚度信息,但本质上得到的地壳厚度仍来源于台站约束。尽管存在部分数值上的差异,本文反演得到的地壳厚度变化与接收函数等方面的结果一致(Yue et al,2012Tian et al,2014):羌塘地块东部和拉萨地块及其以南部分的地壳厚度基本大于60 km,而柴达木盆地的地壳厚度较小,主要在50 km左右。Shin等(2015)研究显示拉萨地块存在一条深达70 km的莫霍面 “条带” ,本文所得的莫霍面深度结果与其相近,造成这种地壳变厚的原因可能是由于印度板块对欧亚板块的挤压作用导致了地壳增厚,而羌塘地块的局部地壳变薄,则可能是印度板块的俯冲作用下发生的岩石圈拆沉和地幔上涌所致(Liu et al,2008Wang et al,2014)。由于莫霍面的研究更注重大尺度的形态分布,本次研究范围较小,因此对于青藏高原整体的莫霍面分布尚需大尺度的成像研究。

    通过对青藏高原中东部地区地壳和上地幔P波速度的联合反演,获得了较为可靠的区域三维P波模型。结果显示,青藏高原中东部地壳和上地幔顶部结构存在强烈的横向不均匀性和纵向不均匀性,显示出印度板块和欧亚板块碰撞下青藏高原地质构造的复杂特征。地壳P波速度、上地幔顶部Pn波速度、地壳厚度等均表现出显著的南北差异特征,印度板块北向俯冲导致青藏高原南部地壳增厚,羌塘地块西部部分地壳隆起。上、下地壳P波速度结构存在明显的不一致性,中、下地壳出现大范围低速,说明其物质结构和物性存在很大差异,推断为在板块俯冲作用下岩石圈发生破裂,上地幔热物质上涌到地壳导致温度上升,从而形成中、下地壳的大范围低速。羌塘、柴达木和巴颜喀拉存在三个稳定的地体且可能延续到上地幔顶部,上地幔顶部Pn波速度的强烈不均匀性和地壳厚度的剧烈变化则反映了印度板块的俯冲、挤压作用下岩石圈的物性和形态的变化特征。本文采用了程序自动拾取走时的方法提取事件初至波走时信息,同时也是联合Pg和Pn波震相走时反演方法在本区域的首次应用。本文得到的三维地壳P波速度模型对后续相关研究具有一定的参考意义,更大尺度的结构特征分析尚需进一步对整个青藏高原进行高精度成像工作。

    IRIS数据中心为本研究提供了数据支持,武汉大学李孟奎博士和美国伊利诺伊大学厄巴纳-香槟分校李江涛博士对本研究提出了建议,作者在此一并表示感谢。

  • 图  1   全球矩心矩张量地震分布.区域1主要是中国大陆西部地区、南亚、中亚等;区域2主要是中国台湾地区、菲律宾、东海及黄海大陆架等;区域3主要是日本及延伸的大陆架地区等

    Figure  1.   Distribution of Global CMT earthquakes. Region 1 contains most of China,South Asia and Middle Asia; Region 2 contains China Taiwan region,Philippines,East China and Yellow Sea continental shelf; Region 3 mainly contains Japan and its extended continental shelf

    图  2   各地震密集区域内面波震级MS与地震矩M0的经验关系

    Figure  2.   The empirical relationship between MS and M0 of different concentrated local regions

    图  3   地震矩M0与平均滑动位移$\overline D$、断层破裂面积A的经验关系

    Figure  3.   The empirical relationship between M0 and average sliding displacement $\overline D$ or rupture area A

    图  4   凹凸体相关震源参数经验关系统计

    (a) 凹凸体面积Aa与地震矩Ma的经验关系;(b) 凹凸体面积Aa与断层破裂面积A的经验关系;(c) 凹凸体面积Aa与面波震级MS的经验关系;(d) 凹凸体长度La与断层长度L的经验关系

    Figure  4.   The empirical relationship statistical of asperity

    (a) The relationship between asperity area Aa and seismic moment Ma;(b) The relationship between asperity area Aa and fault area A;(c) The relationship between asperity area Aa and magnitude MS;(d) The relationship between asperity length La and fault length L

    图  5   本文部分震源参数经验关系与前人统计震源参数经验关系对比

    (a) 地震矩与面波震级经验关系;(b) 地震矩与断层破裂面积经验关系;(c) 地震矩与矩震级经验关系;(d) 凹凸体面积与地震矩之间的关系;(e) 地震矩与平均滑动位移之间的经验关系

    Figure  5.   The source parameters relationship comparing between predecessors and this paper

    (a) The relationship between seismic moment and MS magnitude;(b) The relationship between seismic moment and fault rupture area;(c) The relationship between magnitude MW and seismic moment;(d) The relationship between seismic moment and asperity area;(e) The empirical relationship between seismic moment and average sliding displacement

    表  1   各局部区域相关信息

    Table  1   Information of different local regions

    区域经纬度范围事件个数MW
    1(63°E—110°E),(2.51°N—57.53°N)4385.5—7.9
    2(115°E—135°E),(2.51°N—45°N)5265.5—7.7
    3(135°E—149.94°E),(2.51°N—57.53°N)7855.5—9.1
    下载: 导出CSV

    表  2   各区域内平均震源深度及出现频率较高的震源深度

    Table  2   The average source depth of different local regions and more occurrences source depths

    区域平均震源深度/km众数深度1众数深度2众数深度3众数深度4
    118.74915 km (157)12 km (54)10 km (32)33 km (21)
    223.57815 km (132)12 km (28)24 km (11)10 km (8)
    323.80615 km (186)12 km (82)20 km (15)16 km (12)
      注:众数深度表示出现该震源深度的频率。例如15 km (157),表明震源深度为15 km的震例有157个。
    下载: 导出CSV

    表  3   各区域面波震级与地震矩经验关系参数

    Table  3   Parameters of empirical relationship between MS and M0 in different regions

    区域斜率a标准差截距b标准差皮尔逊相关系数r校正决定系数
    11.0570.01818.7870.1050.940.89
    21.0660.01918.8430.1130.920.85
    31.1560.01918.3050.1120.910.82
    全区域1.0760.01018.7530.0580.910.84
    下载: 导出CSV

    表  4   相对局部区域地震记录震源参数信息

    Table  4   The information of source parameters in relative local region

    位置地震MW地震矩M0/(N·m)破裂面积/km2平均滑动量/cm
    日本新泻7.592.72×10201 800503.70
    中国松潘6.374.03×101821662.19
    伊朗塔巴斯戈尔山7.391.37×10201 628280.51
    中国台湾花莲6.374.03×101848027.99
    中国台湾花莲7.331.11×10201 248296.47
    日本秋田6.131.76×101815438.10
    日本瓦卡萨湾6.282.95×101816061.46
    日本岐阜6.343.63×101818067.22
    中国唐山7.461.74×10201 680345.24
    伊朗鲍勃探戈5.897.67×101716815.22
    日本伊豆大岛6.711.30×101950086.67
    日本长崎西武6.242.57×10189689.24
    中国澜沧—耿马7.135.56×10191 600115.83
    日本伊豆奥基6.547.24×1018198121.89
    伊朗达什—巴亚兹7.237.85×10192 200118.94
    中国炉霍7.471.80×10201 430419.58
    中国海城gnehciah6.993.43×1019900127.04
    日本伊豆河东6.394.32×1018140102.86
    中国道孚6.641.02×101969049.28
    日本伊豆大岛6.374.03×101822858.92
    菲律宾吕宋岛7.744.57×10202 400634.72
    日本大分县6.323.39×1018100113.00
    中国松潘6.711.30×1019360120.37
    中国松潘6.588.32×1018242114.60
    日本大町5.512.07×10172824.64
    日本瓦茨5.906.00×101714413.89
    中国台湾集集7.602.70×10203 432262.24
    日本格约6.803.30×1019630174.60
    日本神户6.902.40×10191 20066.67
    日本鹿儿岛6.101.20×102521618.52
    日本山口5.806.00×10242248.93
    土耳其Kocach7.401.52×10273 282154.36
    日本鸟取6.801.90×102669391.39
    注:本表中各地震样本的震源参数信息参考来源为:① 王海云(2004)博士论文第二章P23—26的表2.1-地震震源参数表;② Wells和Coppersmith (1994)data base节中的Table1.
    下载: 导出CSV

    表  5   地震矩M0与平均滑动位移$\overline D$、破裂面积A的经验关系拟合指标

    Table  5   The empirical relationship parameters between M0 and average sliding displacement $\overline D$ or rupture area A

    经验关系斜率a标准差截距b标准差皮尔逊相关系数r校正决定系数
    M0 & A0.040.960.930.85
    M0 & $\overline D$0.041.040.900.80
    下载: 导出CSV

    表  6   凹凸体相关震源参数信息汇总

    Table  6   The information summary of asperity parameters

    序号位置地震东经
    北纬
    MS地震矩
    /(N·m)
    破裂长度
    /km
    破裂面积
    /km2
    凹凸体长度
    /km
    凹凸体面积
    /km2
    1日本浓尾136.6035.608.01.50×1020801 20016.23243.45
    2中国海原105.7036.708.53.00×10212204 400150.433 008.62
    3日本北伊豆地震138.9835.047.33.78×1019354201.1513.79
    4中国可可托海89.9046.747.98.24×10201803 60086.261 725.17
    5中国昌马97.0039.707.76.30×10201482 22072.771 091.52
    6日本鸟取县134.1835.477.46.65×10183342922.80296.37
    7中国当雄91.4031.108.04.96×10202002 00028.10280.99
    8中国当雄91.5031.007.56.89×1019581 04431.85573.28
    9蒙古国戈壁阿尔泰省99.2045.208.01.22×10213006 00029.90597.96
    10中国通海102.6024.107.78.50×1019751 1251.7225.81
    11中国炉霍100.4031.507.61.80×1020901 17060.71789.28
    12中国道孚101.3530.866.87.28×10185571538.13495.74
    13伊朗科里59.5834.457.14.61×1019751 65051.971143.24
    14中国澜沧-耿马99.6823.007.33.66×1019801 60044.99899.83
    15中国云南耿马99.7022.807.25.50×10194692017.73354.63
    16中国昆仑山92.9135.808.05.90×102042612 780273.998 219.59
    17中国汶川104.1031.448.08.97×10202277 02979.111 186.62
    18日本神户134.9934.786.62.43×10195892867.851 119.61
    19伊朗塔巴斯戈尔山57.0233.377.45.80×1019924 04820.00300.00
    20中国玉树96.8233.107.12.53×1019731 73343.34156.71
    21中国芦山103.0030.307.01.02×1019652 015*80.00
    22中国唐山118.1839.637.92.77×10201153 300*237.00
    23中国伽师76.9339.616.12.06×101813130*26.94
    注:“*”表示该参数暂时缺乏.震源参数参考文献来源为王椿镛等(1978)唐荣昌等(1980)王凯等(1991)国家地震局兰州地震研究所胡方秋和刘景元(1994)马淑田等(1997)Somerville等(1999)王海云(20042010),王海云和陶夏新(2005)吴迪(2008)李启成(2010)张军龙等(2010)李正芳(20132014),温瑞智等(2013)药晓东和章文波(2015)药晓东等(2015)
    下载: 导出CSV

    表  7   凹凸体相关震源参数的经验关系拟合指标

    Table  7   The empirical relationship fitting parameters about asperity

    经验关系斜率a截距b皮尔逊相关系数r校正决定系数
    凹凸体地震矩Ma与凹凸体面积Aa 0.415 1.515 0.655 0.399
    断层面积A与凹凸体面积Aa 0.979 −0.556 0.655 0.402
    断层长度L与凹凸体长度La 0.930 −0.145 0.817 0.644
    面波震级MS与凹凸体面积Aa 0.572 −1.650 0.656 0.397
    下载: 导出CSV

    表  8   主要震源参数经验关系汇总

    Table  8   The summary of empirical relationship for main source parameters

    经验关系     统计公式   皮尔逊相关系数r
    区域1地震矩M0与面波震级MS lgM0=1.057MS+18.787 0.94
    区域2地震矩M0与面波震级MS lgM0=1.066MS+18.843 0.92
    区域3地震矩M0与面波震级MS lgM0=1.156MS+18.305 0.91
    全区域地震矩M0与面波震级MS lgM0= 1.076MS+18.754 0.91
    地震矩M0与断层破裂面积A lgA= 0.542lgM0-4.388 0.93
    地震矩M0与平均滑动位移$\overline{ {D} }$ lg$\overline { {D} }$= 0.46lgM0-10 0.90
    凹凸体地震矩Ma与凹凸体面积Aa lgMa= 0.415lgAa+1.515 0.66
    断层破裂面积A与凹凸体面积Aa lgAa= 0.979lgA-0.556 0.66
    断层长度L与凹凸体长度La lgLa= 0.93lgL-0.145 0.66
    面波震级MS与凹凸体面积Aa lgAa= 0.572lgMS-1.65 0.82
    下载: 导出CSV
  • 车自成, 罗金海, 刘良. 2011. 中国及其邻区区域大地构造学[M]. 2版. 北京: 科学出版社: 1–20.

    Che Z C, Luo J H, Liu L. 2011. Introduction to the New Regional Geotectonics of China and Its Adjacent Area[M]. 2nd ed. Beijing: Science Press: 1–20 (in Chinese).

    陈运泰 顾浩鼎. 2012.震源理论基础[M]. 北京: 中国科学院研究生院: 175–185.

    Chen Y T, Gu H D. 2012. The Basis Theory of Seismic Source[M]. Beijing: Graduate University of Chinese Academy of Sciences: 175–185 (in Chinese).

    程佳. 2017. 川滇地区地震危险性预测模型[D]. 北京: 中国地震局地质研究所: 1–10.

    Cheng J. 2017. Seismic Hazard Modeling of the Sichuan-Yunnan Regions[D]. Beijing: Institute of Geology, China Earthquake Administration: 1–10 (in Chinese).

    国家地震局兰州地震研究所. 1992. 昌马活动断裂带[M]. 北京: 地震出版社: 1–20.

    China Earthquake Administration Lanzhou Institute of Seismology. 1992. Changma Active Fault Belt[M]. Beijing: Seismological Press (in Chinese).

    国家地震局《一九七六年唐山地震》编辑组. 1982. 一九七六年唐山地震[M]. 北京: 地震出版社: 1–20.

    Editorial Group of 1976 Tangshan Earthquake, State Seismological Bureau. 1982. The Mammoth Tangshan Earthquake of 1976[M]. Beijing: Seismological Press: 1–20 (in Chinese).

    郭履灿,庞明虎. 1981. 面波震级和它的台基校正值[J]. 地震学报,3(3):312–320.

    Guo L C,Pang M H. 1981. Surface-wave magnitude of earthquakes and its station correction[J]. Acta Seismologica Sinica,3(3):312–320 (in Chinese).

    胡方秋,刘景元. 1994. 可可托海—阿克塞测深大剖面北段深部构造与地震[J]. 内陆地震,8(2):152–158.

    Hu F Q,Liu J Y. 1994. Deep-seated structures on the Northern part of Koktokay-Aksay deep sounding profile and its correlation with earthquake[J]. Inland Earthquake,8(2):152–158 (in Chinese).

    黄汲清. 1980. 中国大地构造及其演化[M]. 北京: 科学出版社: 1–20.

    Huang J Q. 1980. Geotectonic Evolution of China[M]. Beijing: Science Press: 1–20 (in Chinese).

    李启成. 2010. 经验格林函数方法模拟地震动研究[D]. 哈尔滨: 中国地震局工程力学研究所: 1–10.

    Li Q C. 2010. Research on Ground Motion Simulation with Empirical Green’s Function Method[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administration: 1–10 (in Chinese).

    李四光. 1973. 地质力学概论[M]. 北京: 科学出版社: 1–20.

    Li S G. 1973. Introduction to Geomechanics[M]. Beijing: Science Press: 1–20 (in Chinese).

    李正芳. 2013. 强震破裂面上的不均匀体及其在地震危险性分析中的应用研究[D]. 北京: 中国地震局地质研究所: 33–35.

    Li Z F. 2013. A Study on Inhomogeneous Bodies of Strong Earthquake Rupture Plane and its Application to Seismic Hazard Analysis[D]. Beijing: China Earthquake Administration: 33–35 (in Chinese).

    李正芳. 2014. 强震破裂面上的不均匀体及其在地震危险性分析中的应用研究[J]. 国际地震动态,(9):42–45. doi: 10.3969/j.issn.0235-4975.2014.09.012

    Li Z F. 2014. A study on inhomogeneous bodies of strong earthquake rupture plane and its application to seismic hazard analysis[J]. Recent Developments in World Seismology,(9):42–45 (in Chinese).

    马淑田,姚振兴,纪晨. 1997. 1996年3月19日新疆伽师MS6.9地震的震源机制以及相关问题研究[J]. 地球物理学报,40(6):782–790. doi: 10.3321/j.issn:0001-5733.1997.06.007

    Ma S T,Yao Z X,Ji C. 1997. The focal mechanism solution for March 19,1996 MS6.9 earthquake in Xin-jiang,Jia-shi region and related problems[J]. Acta Geophysica Sinica,40(6):782–790 (in Chinese).

    潘华. 2000. 概率地震危险性分析中参数不确定性研究[D]. 北京: 中国地震局地球物理研究所: 10–20.

    Pan H. 2000. Study on Uncertainties in the Parameters of PSHA[D]. Beijing: Institute of Geophysics, China Earthquake Administration: 10–20 (in Chinese).

    任纪舜,王作勋,陈炳蔚,姜春发,牛宝贵,李锦轶,谢广连,和政军,刘志刚. 1997. 新一代中国大地构造图[J]. 中国区域地质,16(3):225–230,248.

    Ren J S,Wang Z X,Chen B W,Jiang C F,Niu B G,Li J T,Xie G L,He Z J,Liu Z G. 1997. A new generation tectonic map of China[J]. Regional Geology in China,16(3):225–230,248 (in Chinese).

    唐荣昌,李天祒,李介成. 1980. 当雄7.5级地震地质构造背景及其成因的初步认识[J]. 地震研究,(1):87–96.

    Tang R C,Li T S,Li J C. 1980. Preliminary understanding of geological tectonic background and genesis of the Dangxiong MS7.5 earthquake[J]. Journal of Seismological Research,(1):87–96 (in Chinese).

    王椿镛,朱成男,刘玉权. 1978. 用地形变资料测定通海地震的地震断层参数[J]. 地球物理学报,21(3):191–198.

    Wang C Y,Zhu C N,Liu Y Q. 1978. Determination of earthquake fault parameter for the Tonghai earthquake from ground deformation data[J]. Acta Geophysica Sinica,21(3):191–198 (in Chinese).

    王海云. 2004. 近场强地震动预测的有限断层震源模型[D]. 哈尔滨: 中国地震局工程力学研究所: 23–26.

    Wang H Y. 2004. Finite Fault Source Model for Predicting Near-Field Strong Ground Motion[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administration: 23–26 (in Chinese).

    王海云,陶夏新. 2005. 近场强地震动预测中浅源地震的asperity模型特征[J]. 哈尔滨工业大学学报,37(11):1533–1539. doi: 10.3321/j.issn:0367-6234.2005.11.022

    Wang H Y,Tao X X. 2005. Characterizing a shallow earthquake asperity model for predicting near field strong ground motion[J]. Journal of Harbin Institute of Technology,37(11):1533–1539 (in Chinese).

    王海云. 2010. 2010年4月14日玉树MS7.1地震加速度场预测[J]. 地球物理学报,53(10):2345–2354. doi: 10.3969/j.issn.0001-5733.2010.10.008

    Wang H Y. 2010. Prediction of acceleration field of the 14 April 2010 Yushu earthquake[J]. Chinese Journal of Geophysics,53(10):2345–2354 (in Chinese).

    王凯,高莉萍,姚振兴,张受生. 1991. 澜沧—耿马地震的震源机制研究[J]. 地球物理学报,34(5):569–580. doi: 10.3321/j.issn:0001-5733.1991.05.005

    Wang K,Gao L P,Yao Z X,Zhang S S. 1991. Source mechanism of the 1988 Lancang-Gengma China earthquake[J]. Acta Geophysica Sinica,34(5):569–580 (in Chinese).

    温瑞智,任叶飞,齐文浩,卢滔,杨振宇,单振东,汪云龙. 2013. 2013年4月20日芦山地震最大加速度记录分析[J]. 西南交通大学学报,48(5):783–791. doi: 10.3969/j.issn.0258-2724.2013.05.001

    Wen R Z,Ren Y F,Qi W H,Lu T,Yang Z Y,Shan Z D,Wang Y L. 2013. Maximum acceleration recording from Lushan earthquake on April 20,2013[J]. Journal of Southwest Jiaotong University,48(5):783–791 (in Chinese).

    吴迪. 2008. 基于凹凸体模型的地震动半经验合成研究[D]. 上海: 同济大学: 19–30.

    Wu D. 2008. The Semi-empirical Synthesis of Ground Motions Basing on Asperity Model[D]. Shanghai: Tongji University: 19–30 (in Chinese).

    药晓东,章文波. 2015. 2013年四川芦山MS7.0地震强地面运动模拟[J]. 地震学报,37(4):599–616. doi: 10.11939/jass.2015.04.007

    Yao X D,Zhang W B. 2015. Strong ground motion simulation for the 2013 MS7.0 Lushan,China,earthquake[J]. Acta Seismologica Sinica,37(4):599–616 (in Chinese).

    药晓东,章文波,于湘伟. 2015. 2008年汶川8.0级大地震近场强地面运动的模拟[J]. 地球物理学报,58(3):886–903. doi: 10.6038/cjg20150317

    Yao X D,Zhang W B,Yu X W. 2015. Simulation of near-field strong ground motion caused by the 2008 MS8.0 Wenchuan earthquake[J]. Chinese Journal of Geophysics,58(3):886–903 (in Chinese).

    张军龙,陈长云,胡朝忠,杨攀新,熊仁伟,李智敏,任金卫. 2010. 玉树MS7.1地震地表破裂带及其同震位移分布[J]. 地震,30(3):1–12. doi: 10.3969/j.issn.1000-3274.2010.03.001

    Zhang J L,Chen C Y,Hu C Z,Yang P X,Xiong R W,Li Z M,Ren J W. 2010. Surface rupture and coseismic displacement of the Yushu MS7.1 earthquake,China[J]. Earthquake,30(3):1–12 (in Chinese).

    张齐. 2016. 地震动衰减关系的区域性差异研究[D]. 哈尔滨: 中国地震局工程力学研究所: 1–10.

    Zhang Q. 2016. Study on Regional Differentiation of Ground Motion Attenuation Relationship[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administration: 1–10 (in Chinese).

    笠原庆一. 1992.防灾工程学中的地震学[M]. 北京: 地震出版社: 1–20.

    Li Y Q Y. 1992. Seismology in Disaster Prevention Engineering[M]. Beijing: Earthquake Press: 1–20 (in Chinese).

    Abe K. 1975. Reliable estimation of the seismic moment of large earthquakes[J]. J Phys Earth,23(4):381–390. doi: 10.4294/jpe1952.23.381

    Chen S Z,Atkinson G M. 2002. Global comparisons of earthquake source spectra[J]. Bull Seismol Soc Am,92(3):885–895. doi: 10.1785/0120010152

    Chiou B,Youngs R,Abrahamson N,Addo K. 2010. Ground-motion attenuation model for small-to-moderate shallow crustal earthquakes in California and its implications on regionalization of ground-motion prediction models[J]. Earthq Spectra,26(4):907–926. doi: 10.1193/1.3479930

    Douglas J. 2004. An investigation of analysis of variance as a tool for exploring regional differences in strong ground motions[J]. J Seismol,8(4):485–496. doi: 10.1007/s10950-004-3094-7

    Duda S J,Kaiser D. 1989. Spectral magnitudes,magnitude spectra and earthquake quantification; The stability issue of the corner period and of the maximum magnitude for a given earthquake[J]. Tectonophysics,166(1/3):205–211,215–219.

    Gutenberg B. 1945. Magnitude determination for deep-focus earthquakes[J]. Bull Seismol Soc Am,35(3):117–130. doi: 10.1785/BSSA0350030117

    Hanks T C,Kanamori H. 1979. A moment magnitude scale[J]. J Geophys Res,84(B5):2348–2350. doi: 10.1029/JB084iB05p02348

    Hanks T C,Johnston A C. 1992. Common features of the excitation and propagation of strong ground motion for North American earthquakes[J]. Bull Seismol Soc Am,82(B5):1–23.

    Kanamori H. 1977. The energy release in great earthquakes[J]. J Geophys Res,82(20):2981–2987. doi: 10.1029/JB082i020p02981

    Mahani A B,Atkinson G M. 2013. Regional differences in ground-motion amplitudes of small-to-moderate earthquakes across North America[J]. Bull Seismol Soc Am,103(5):2604–2620. doi: 10.1785/0120120350

    Sato R. 1979. Theoretical basis on relationships between focal parameters and earthquake magnitude[J]. J Phys Earth,27(5):353–372. doi: 10.4294/jpe1952.27.353

    Somerville P,Irikura K,Graves R,Sawada S,Wald D,Abrahamson N,Iwasaki Y,Kagawa T,Smith N,Kowada A. 1999. Characterizing crustal earthquake slip models for the prediction of strong ground motion[J]. Seismol Res Lett,70(1):59–80. doi: 10.1785/gssrl.70.1.59

    Wells D L,Coppersmith K J. 1994. New empirical relationships among magnitude,rupture length,rupture width,rupture area,and surface displacement[J]. Bull Seismol Soc Am,84(4):974–1002.

    Miyake H,Iwata T,Irikura K. 2003. Source characterization for broadband ground-motion simulation:Kinematic heterogeneous source model and strong motion generation area[J]. Bull Seismol Soc Am,93(6):2531–2545. doi: 10.1785/0120020183

    Slemmons D B, Bodin P, Zhang X. 1989. Determination of Earthquake Size From Surface Faulting Events[D]. Guangzhou: Proceedings of International Seminar on Seismic Zonation: 157–169.

图(5)  /  表(8)
计量
  • 文章访问数:  634
  • HTML全文浏览量:  271
  • PDF下载量:  76
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-09-01
  • 修回日期:  2020-11-08
  • 网络出版日期:  2021-08-15
  • 发布日期:  2021-07-14

目录

/

返回文章
返回