地震时间间隔的统计分布及其地震危险度 D 值在华北大震前的异常变化

王炜1, 刘震华2

王炜1, 刘震华2. 1987: 地震时间间隔的统计分布及其地震危险度 D 值在华北大震前的异常变化. 地震学报, 9(2): 113-127.
引用本文: 王炜1, 刘震华2. 1987: 地震时间间隔的统计分布及其地震危险度 D 值在华北大震前的异常变化. 地震学报, 9(2): 113-127.
WANG WEIup, LIU ZHENHUAup2com sh a. 1987: THE STATISTICAL DISTRIBUTION OF TIME INTERVALS BETWEEN TWO SUCCESSIVE EARTHQUAKES AND ABNORMAL VARIATION OF THE D-VALUE, THE DEGREE OF SEISMIC DANGER PRIOR TO SEVERAL GREATER EARTHQUAKES IN NORTH CHINA. Acta Seismologica Sinica, 9(2): 113-127.
Citation: WANG WEIup, LIU ZHENHUAup2com sh a. 1987: THE STATISTICAL DISTRIBUTION OF TIME INTERVALS BETWEEN TWO SUCCESSIVE EARTHQUAKES AND ABNORMAL VARIATION OF THE D-VALUE, THE DEGREE OF SEISMIC DANGER PRIOR TO SEVERAL GREATER EARTHQUAKES IN NORTH CHINA. Acta Seismologica Sinica, 9(2): 113-127.

地震时间间隔的统计分布及其地震危险度 D 值在华北大震前的异常变化

THE STATISTICAL DISTRIBUTION OF TIME INTERVALS BETWEEN TWO SUCCESSIVE EARTHQUAKES AND ABNORMAL VARIATION OF THE D-VALUE, THE DEGREE OF SEISMIC DANGER PRIOR TO SEVERAL GREATER EARTHQUAKES IN NORTH CHINA

  • 摘要: 本文探讨了华北小震时间间隔的统计分布规律,并且经过统计检验认为,它服从韦布尔(Weibull)分布,f()=-1exp[-(/)].实际计算表明,地震发生方式参数一般小于1,说明华北小震活动呈续发性方式.本文利用可靠性理论中的一些成果,引入地震危险度函数D(t)=1-exp[-(t/)],它表示 t 时刻内发生地震的概率.本文取 t=1d,将 D(1)称为地震危险度 D 值.得到自1970-1984年期间华北共发生的九次 Ms5.8地震前,震中周围区域的地震危险度 D 值,都出现一个明显的异常上升过程.文中还就,和 D 值与实际地震活动的关系,以及 D 值变化的稳定性等问题进行了讨论,认为 D 值是应力的函数.同时,可靠性理论中的其它一些参数,也有着与 D 值相类似的变化.
    Abstract: In this paper the statistical distribution of the time interval between two successive earthquakes is discussed, statistically tested, and shown that it obeys Weibull distribution f() =-1exp[-(/)].Actual calculation results show that in general the seismicoccurence fashionparametre, the p-value, is less than 1.0. This means that the pattern of seismic occurence is successive. Using some results of the reliability theory, the degree of seismic danger function D(t) =1-exp[-(t/)] is introduced, it is the probability of occurrence of an earthquake in timet. In this paper we let t = 1 day and call D(1) the degree of seismic danger, the D-value. The results indicate that the degree of seismic danger, the D-value showed a significant rising trend before all nine earthquakes of Ms5.8 in North China (1970-1984). The relation between , , D-value and the actual seismic activity, stability of the calculated D-value etc. are discussed. Besides, it is thought that the D-value is a function of stress. The other parametres of reliability theory have similar variations as the D-value.
  • 地震动参数及其衰减关系是地震工程的重点研究对象, 在地震区划和工程抗震设防中占有重要地位.地震动参数包括地震动三要素即地震动幅值、频谱和持时中的一个或多个,如幅值特征参数峰值加速度、频谱特征参数反应谱、持时特征参数括弧持时等.

    地震动作用下工程结构的破坏情况与其吸收的地震动总能量相关.从能量的角度出发,Arias(1970)提出一种新型的地震动参数阿里亚斯强度,表示输入到单自由度体系单位质量的能量.阿里亚斯强度包含地震动的幅值和持时信息,其本身又是一种能量参数,与地震引起的灾害现象具有很强的相关性,被广泛应用于地震灾害危险性分析. Jibson(1993)发现与地震滑坡有关的纽马克(Newmark)累积位移与阿里亚斯强度具有很好的相关性;Cabanas等(1997)以当地建筑物损坏情况作为地震烈度评定标准的分析结果表明,阿里亚斯强度与本地烈度拟合得很好;王秀英等(2010)的研究显示汶川地震中整个龙门山震区的阿里亚斯强度能够诱发崩塌滑坡的水平向下限为0.4 m/s.

    地震动衰减关系简化了地震的震级、断层类型、传播路径和场地条件等因素对地震动特征的影响,利用已知的地震动分布规律,将地震动参数表达为以上因素的函数.鉴于阿里亚斯强度在地震灾害危险性分析中的广泛应用,Abdrakhmatov等(2003)Travasarou等(2003)Stafford等(2009)Lee等(2012)基于不断丰富的强震动数据,建立了不同地区的阿里亚斯强度衰减关系.由于缺乏强震动数据资料,我国大陆地区尚未建立明确的阿里亚斯强度衰减关系.

    不同区域地壳结构和构造特征的不同导致地震动特性及其衰减关系具有差异性,我国大陆地区与美国西部、中亚和我国台湾地区同为构造活动区,地下介质品质因子Q值较低,地震波非弹性衰减较高,所发地震多为浅源强震,地震动衰减关系相近,但与日本、欧洲等地区的地震勘探特性及其衰减关系差异较大(Chen,Atkinson,2002Douglas,2004Akkar,Bommer,2007; 靳超宇,俞言祥,2009).

    强震目录数据库NGA-west1(Next Generation Attenuation relationships for western US)中不包含我国大陆地震(Chiou et al, 2008),NGA-west2中增加了我国2008汶川MW7.9地震(Ancheta et al,2014),但其强震动加速度记录获取受限.考虑到未来我国大陆地区阿里亚斯强度研究的需求以及数据格式和标度(面波震级或矩震级、震中距或断层投影距等)的差异性,本文拟选取NGA-west1数据库中与我国大陆构造环境相近的美国西部、中亚和我国台湾地区的强震记录,采用肖亮和俞言祥(2010)提出的分步回归法求取回归系数,对阿里亚斯强度衰减关系予以分析,以期获取强震动的震级饱和及近场距离饱和特征,检验回归所得结果与实际数据的拟合程度.

    阿里亚斯强度的表达式为

    (1)

    式中,g为当地重力加速度,a(t)为某一分量的地震动加速度时程,T为地震动加速度记录的总持续时间,λ为结构物的阻尼比.对于一般的结构物,λ取值为0.05,则

    (2)

    进而可以近似得到阿里亚斯强度的简洁形式:

    (3)

    根据NGA数据库中GMX’s C1的分类方法,选择标注为自由地表或4层(含)以下轻型建筑物底层或附近的台站记录,以除去土结相互作用的影响;同时为了减少地形效应带来的影响,剔除隧道、大坝等特殊地形附近的台站记录.由于部分台站三分量记录不完整,为保证数据平衡,剔除水平分量不完整的台站记录.以场地地表 30 m土层等效剪切波波速vS30为分类依据,相应的场地分类见表 1.

    表  1  本文场地分类标准
    Table  1.  Site classification standards in this study
    场地类别 vS30/(m·s-1) 代码
    基岩 >500 A
    硬土 250-500 B
    软土 <250 C
    下载: 导出CSV 
    | 显示表格

    本文集中分析主震的阿里亚斯强度衰减关系,故将1999年台湾集集MW7.6地震的余震记录全部除去.由于此次地震主震的地震动相对于同样震级的其它地震偏低,其强震记录共401组,数量远超其它地震;为避免其对衰减关系高震级段的回归产生太大影响,从中随机选取50组(100条)数据纳入衰减关系的分析中.最终用于本文的地震目录见表 2.

    表  2  本文所用地震目录
    Table  2.  Catalog of earthquakes used in this study
    序号 地震 发震时间 MW 记录条数
    1 Kern County 1952-07-21 7.4 4
    2 Parkfield 1966-06-28 6.2 10
    3 Borrego Mountain 1968-04-09 6.6 2
    4 San Fernando 1971-02-09 6.6 56
    5 Tabas 1978-09-16 7.4 2
    6 Coyote Lake 1979-08-06 5.7 18
    7 Imperial Valley-06 1979-10-15 6.5 60
    8 Victoria 1980-06-09 6.3 2
    9 Westmorland 1981-04-26 5.9 12
    10 Coalinga-01 1983-05-02 6.4 90
    11 Borah Peak-01 1983-10-28 6.9 4
    12 Morgan Hill 1984-04-24 6.2 46
    13 North Palm Springs 1986-07-08 6.1 52
    14 Chalfant Valley-01 1986-07-20 5.8 10
    15 Chalfant Valley-02 1986-70-21 6.2 20
    16 Whittier Narrows-01 1987-10-01 6.0 208
    17 Superstition Hills-01 1987-11-24 6.2 2
    18 Superstition Hills-02 1987-11-24 6.5 6
    19 Loma Prieta 1989-10-18 6.9 146
    20Cape Mendocino 1992-04-25 7.0 8
    21 Landers 1992-06-28 7.3 136
    22 Northridge-01 1994-01-17 6.7 282
    23 Double Springs 1994-09-12 5.9 2
    24 Chi-Chi 1999-09-20 7.6 100
    25 Manji-01 1990-06-20 7.4 2
    26 Sierra Madre 1991-06-28 5.6 16
    27 Little Skull Mountain 1992-06-29 5.7 16
    28 Hector Mine 1999-10-16 7.1 158
    下载: 导出CSV 
    | 显示表格

    本文采用的震级为MW,取值范围为5.5—8.0;距离为断层投影距RJB, 取值范围为0—200 km.原始记录中地震为127次,强震动记录共有3275组,每组含三分量(部分分量无记录);筛选后的地震为28次,强震动记录为735组,两个水平分量看作两条独立的记录,共有1470条记录进行分析.数据整体的分布状况如图 1所示.

    图  1  强震资料的震级MW随断层投影距RJB分布图
    Figure  1.  Magnitude MW versus distance RJB of strong-motion data

    衰减关系中的协变量有矩震级MW、断层投影距RJB、场地类别和断层类型.按照表 1的分类标准,硬土场地共有记录1050条,基岩场地有238条,软土场地有182条.土层对地震动的放大效应表现为在衰减关系上增加一个常数.由于硬土场地记录较多,本文以硬土场地作为主要对象,基岩和软土场地的阿里亚斯强度衰减关系通过硬土场地加减常数得到.设定SASC两个哑指标表示场地类别,即SA=1和SC=0表示基岩场地;SA=0和SC=0表示硬土场地; SA=0和SC=1表示软土场地.考虑到正断层与走滑断层引起的地震动相近,且二者均小于逆断层,本文采用哑指标V表示断层类型:V=1表示逆断层,V=0表示其它类型断层.

    本文采用的衰减关系为

    (4)

    式中:Ia为阿里亚斯强度,单位为m/s;M(MW)为震级项,是关于MW的函数;G(MWRJB)为矩震级与断层投影距的耦合项;S(SASC)为场地类型的函数,F(V)为与断层类型相关的函数,F(V)=fVf为系数,V为哑指标;ε为期望值为0、标准差为σ的随机量.

    场地类别和断层类型对阿里亚斯强度的影响采用简单的线性形式表示:S(SA, SC)=mSA+nSC,式中mn为系数.

    首先对硬土场地的记录进行分析.对于硬土场地,S=0,对回归分析无影响,因此式(4)可写为

    (5)

    式中:Iaij为第i个地震事件在第j个台站引起的阿里亚斯强度;Li为第i个地震事件的参数,与震级MW和断层类型有关,对于单个地震,其为常数;clg(Rij+Hi)为距离衰减项,在MW5.5—8.0内,系数c可以看作一个常数,Hi表示与第i个地震事件相关的近断层饱和项参数,采用Campbell(1981)提出并使用多年的形式,即Hi=dexp(eMWi),其中de为系数; ε为期望值为0、标准差为σ的随机量.

    距离衰减项中断层投影距与距震级耦合在一起,回归系数的求取对数据的分布比较敏感,本文的近场饱和项系数经验性地借鉴肖亮和俞言祥(2010)的结果,取d=0.956, e=0.462.我们选取震级不同,距离分布较广的多组记录进行分析,求取每个地震事件的回归系数ci,取其平均后得到c=2.494,则距离衰减项为

    (6)

    考虑到地震动参数的震级饱和性,采用与Boore和Atkinson(2012)类似的方法,Li的回归采用分段函数的形式,具体如下:

    (7)

    对式(7)进行回归可得系数a1=-3.407,b1=1.065,a2=-1.073,b2=0.715,f=0.089.

    在上述对硬土数据分析的基础上,增加基岩和软土的数据,对所有数据进行分析,进而得到最终的阿里亚斯强度相对于矩震级、断层投影距、断层类型和场地类别的衰减关系,即

    (8)

    式中系数mn如前所述, 下标A, B, C为表 1中场地类型代码.对式(8)进行回归分析可得系数m=-0.262,n=0.099.这样,可得最终的衰减系数(表 3),相应的衰减关系为

    表  3  阿里亚斯强度衰减关系系数
    Table  3.  Coefficients of attenuation relations of Arias intensity
    震级 a b c d e m n f
    MW≤6.5 -3.407 1.065 2.494 0.956 0.462 -0.262 0.099 0.089
    MW> 6.5 -1.073 0.715 2.494 0.956 0.462 -0.262 0.099 0.089
    下载: 导出CSV 
    | 显示表格

    (9)

    将得到的阿里亚斯强度衰减关系式与实际数据进行对比,观察数据点相对于衰减曲线的分布状况,如图 2所示.可知一定震级范围(如MW7.0—8.0)内的阿里亚斯强度值,大部分均位于相应的两条衰减曲线MW=7.0与MW=8.0之间.对于不同的场地类别或断层类型,以上所得衰减关系均与实际数据拟合得较好.

    图  2  阿里亚斯强度衰减曲线与实际数据的对比
    (a)硬土场地,正断层或走滑断层震源机制;(b)基岩场地,逆断层震源机制
    Figure  2.  Comparison of attenuation relation curves of Arias intensity with observed records
    (a) Stiff soil site, normal fault or strike-slip fault; (b) Rock site, reverse fault

    图 2可以看出阿里亚斯强度衰减关系的计算结果与实际数据拟合得较好.接下来计算实际与拟合的阿里亚斯强度残差,并观察其随矩震级或断层投影距这两个主要变量的变化情况,进一步检查拟合的优劣.由图 3可以看出:拟合残差随矩震级的波动较大,MW7.6为集集地震的主震,其阿里亚斯强度偏低;拟合残差随断层投影距的波动较小,整体较对称地分布于零线两侧.

    图  3  阿里亚斯强度拟合残差随断层投影距RJB和矩震级MW的变化
    Figure  3.  The fitting residuals of Arias intennity against with the magnitude and Joyner-Boor distance

    由于研究区域不同以及所用强震动数据资料存在差异,导致了衰减关系的震级项、距离项和场地项所选取的形式不同,由此得到不同的衰减模型.为说明本文方法和结果的合理性,将其与Travasarou等(2003)Stafford等(2009)文中的模型Ⅰ以及Lee等(2012)所用的衰减模型进行比较,结果如图 4所示.

    图  4  本文结果与基于其它3种模型所得的基岩场地逆断层震源机制下阿里亚斯强度衰减曲线的对比
    Figure  4.  Comparison of attenuation relation curves of Arias intensity for rock site and reverse fault from this study, with those from other three models

    Travasarou等(2003)vS30≥760 m/s或土层厚度<6 m归类为基岩,震源机制分为正断层、逆断层和走滑断层等3类,研究了全球活动板块边界的75次浅源地震;Stafford等(2009)将基岩定义为单轴抗压强度大于1 MPa且vS30>360 m/s,震源机制分为逆断层和非逆断层两类,以此对新西兰的23次地震进行分析;Lee等(2012)采用的震源机制分类标准与Travasarou等(2003)相同,场地项不再进行分类,统一表示成vS30的自然对数,研究对象为台湾的62次地震.

    Stafford等(2009)的模型Ⅰ中包含震源深度项,本文未将其纳入衰减关系中,在应用其模型时,令震源深度为本文所用地震的加权平均值11.043 km;Lee等(2012)的模型中直接包含vS30的对数项,本文在应用其模型时,令vS30=760 m/s.

    图 4可以看出,4种模型均考虑了近场的距离饱和效应,所得阿里亚斯强度在MW为6.5时比较接近.震级方面:本文与Travasarou等(2003)近场震级饱和效应较明显;Stafford等(2009)的模型1直接采用震级的线性形式,未考虑震级饱和;Lee等(2012)对线性震级进行了对数修正,震级饱和效应不明显.距离衰减方面,4种模型得到的衰减系数不同,在远场区有较大差异.

    综上,本文的模型采用了表示震级饱和的分段线性函数,距离衰减项体现了近场的距离饱和特性,与其它衰减模型相比,在MW6.5或远场区域比较接近,近场区域在震级过高(MW=7.5)或震级过低(MW=5.5)时有一定差异.

    近年来我国大陆地区逐渐获得了大量强震动加速度记录数据.为检验本文所得衰减关系对大陆数据的适用性,特选取两个典型地震即2008年汶川MW7.9地震和2013年芦山MW6.8地震,将实际数据计算得到的阿里亚斯强度与本文模型进行对比.图 5分别给出了芦山地震和汶川地震阿里亚斯强度衰减曲线与实际数据的对比结果.由于国内强震记录中场地仅包含基岩和土层两种类型,且前文得到的阿里亚斯强度衰减关系中硬土与软土差异较小,图中土层项衰减曲线对应于前文的硬土场地.国内数据未给出断层投影距RJB,对于芦山地震,断层尺度较小,震中距与断层投影距相近,实际数据点的RJB用震中距代替;对于汶川地震,由于其断层破裂尺度较大(长300 km),用震中距代替RJB会出现很大偏差,因此本文根据NGA-west2地震目录中给出的台站位置和RJB与国内强震记录提供的台站位置对照,阿里亚斯强度由国内强震记录计算得出,RJB取NGA-west2中数值.

    图  5  芦山MW6.8(a)和汶川MW7.9(b)地震的阿里亚斯强度衰减曲线与实际数据对比
    Figure  5.  Comparison of attenuation relation curves of Arias intensity with observed records for Lushan MW6.8 earthquake (a) and Wenchuan MW7.9 earthquake (b)

    基岩场地记录数较少且整体无太大偏差,不作分析.对于土层场地的数据,由图 5可以看出:芦山地震数据离散性较大,近场数据较预测值偏大,远场数据分布在预测值两侧,数据整体拟合得较好;汶川地震阿里亚斯强度离散性较小,实际数据整体较预测值偏大且随距离的衰减较慢,一方面主要是由于1999年台湾集集MW7.6地震导致高震级地震的阿里亚斯强度拟合值偏低,另一方面,汶川地震断层破裂长度太大(300 km),本文衰减关系分析中RJB取值不超过200 km,在该距离范围内,场点仍可看作近场,故实际数据的衰减速度偏小.

    考虑到未来我国大陆地区阿里亚斯强度研究的需求,本文对NGA-west1数据库中的强震资料进行了针对性的筛选,初步建立了以美国西部数据为主的阿里亚斯强度的衰减关系.该模型很好地体现了阿里亚斯强度近场饱和、震级饱和特性,采用分步回归的方式可以减少变量之间相关性的影响.

    回归过程中对地震动偏低且强震记录较多的集集地震记录进行了缩减,以降低其对高震级回归所产生的影响.由于研究区域和所采用模型的差异,不同研究人员所得阿里亚斯强度衰减关系的系数不同,将本文结果与相关的3个衰减模型进行比较,结果显示对于中等震级(MW6.5)而言,本文与其结果相近,而对于高震级或地震级由于近场和距离饱和项所采用的模型形式不同,本文与其结果差异较大.

    我国大陆强震数据与本文结果拟合得较好,随着国内强震动数据的增多,可以根据已有的结果,利用转换方法建立我国大陆地区的阿里亚斯强度衰减关系.随着震源机理的深入研究和处理方法的改进,现有的衰减模型还可进一步修正细化,例如丰富的近断层样本可用于研究阿里亚斯强度的方向性效应和上下盘效应等,而且可以并入衰减关系中;进一步的残差分析可用于区域的概率危险性分析.

  • [1] Hagiwara, Y., Probability of earthquake occurrence as obtained from Weibull distribution analysis of crustal strain, Tectonorhysics, 23, 313——318, 1974,

    [2] Rikitake, T., Recurrence of great earthquakes of subdnction zones, Tectonophysics, 35, 335——362,1976.

    [3] 郑建中,中国大地震的发生方式、时间间隔及概率分布,地震研究,6,增刊,439——447, 1983,

    [4] 张国民、傅征祥,华北强震的时间分布及其物理解释,地球物理学报,28, 569—— 578, 1985,

    [5] 松村正三、郑建中译,地震时间序列的——个参数表示及其在地震预报中的应用,世界地震译丛,1983, No. 2,11——15,

    [6] 国家地震局分析预报中心,中国东部地震目录,地震出版社,1981,

    [7] 紊藤嘉博,五所译,可靠性基础数学,69——77,国防工业出版社,1977,

    [8] 茹诗松,可靠性统计,30——80,华东师范大学出版社,1984,

    [9] 中国科学院计算中心概率统计组,概率统计计算,109——115,科学出版社,1979,

    [10] 姜秀娥、陈非比,区域震群与唐山大震,地震学报,5, 145——157, 1983,

    [11] 朱传镇、罗胜利,震群与大地震关系的研究,地震学报,3, 105——116, 1981,

    [12] 顾浩鼎、曹天青,前兆震群与S波偏振,地震学报,1, 343——355, 1980,

    [13] Scholz, C. H., Static fatigue of quartz, J. G. R, 77, 2104——2114, 1972,

    [14] 宇津德治、陈铁成译,地震学,176——179,地震出版社,1981,

    [15] 安艺敬一、宋良玉译,前兆现象的概率综合,世界地震译丛,1982, 6: 22——29,

    [16] Jones, L. M., P. Molnar,朱传镇译,前震的某些特征,国外地震,1979, 6: 7——19.

    [1] Hagiwara, Y., Probability of earthquake occurrence as obtained from Weibull distribution analysis of crustal strain, Tectonorhysics, 23, 313——318, 1974,

    [2] Rikitake, T., Recurrence of great earthquakes of subdnction zones, Tectonophysics, 35, 335——362,1976.

    [3] 郑建中,中国大地震的发生方式、时间间隔及概率分布,地震研究,6,增刊,439——447, 1983,

    [4] 张国民、傅征祥,华北强震的时间分布及其物理解释,地球物理学报,28, 569—— 578, 1985,

    [5] 松村正三、郑建中译,地震时间序列的——个参数表示及其在地震预报中的应用,世界地震译丛,1983, No. 2,11——15,

    [6] 国家地震局分析预报中心,中国东部地震目录,地震出版社,1981,

    [7] 紊藤嘉博,五所译,可靠性基础数学,69——77,国防工业出版社,1977,

    [8] 茹诗松,可靠性统计,30——80,华东师范大学出版社,1984,

    [9] 中国科学院计算中心概率统计组,概率统计计算,109——115,科学出版社,1979,

    [10] 姜秀娥、陈非比,区域震群与唐山大震,地震学报,5, 145——157, 1983,

    [11] 朱传镇、罗胜利,震群与大地震关系的研究,地震学报,3, 105——116, 1981,

    [12] 顾浩鼎、曹天青,前兆震群与S波偏振,地震学报,1, 343——355, 1980,

    [13] Scholz, C. H., Static fatigue of quartz, J. G. R, 77, 2104——2114, 1972,

    [14] 宇津德治、陈铁成译,地震学,176——179,地震出版社,1981,

    [15] 安艺敬一、宋良玉译,前兆现象的概率综合,世界地震译丛,1982, 6: 22——29,

    [16] Jones, L. M., P. Molnar,朱传镇译,前震的某些特征,国外地震,1979, 6: 7——19.

计量
  • 文章访问数:  1088
  • HTML全文浏览量:  22
  • PDF下载量:  36
  • 被引次数: 0
出版历程
  • 发布日期:  2011-09-01

目录

/

返回文章
返回