余震序列

摘要： 针对2000—2020年南北地震带21个M_(S)≥6.0地震序列,利用推定最大余震震级法,采用不同的最小完备震级M_(C)和b值计算方法,分别给出了实际最大余震震级的估计值,分析了估算效果。结果显示:从整体上看,推定最大余震震级法可以给出实际最大余震震级的有效估计值,而且随着时间推移,估计误差逐步减小、趋于稳定。从计算方法来看,使用GFT和修正MAXC的混合方法计算最小完备震级,并且采用最大似然法计算b值时,推定最大余震震级法的估算误差较小,且估算正确的结果较多。

摘要： 为系统地考察中国大陆地区强震之后的早期余震概率预测效能及制约因素,以期开展适合板内地震活跃地区的“可操作”强余震概率预测,服务震后抗震救灾、应急救援等工作,本文采用充分利用震后早期小震信息的Omi-R-J模型对中国大陆的86个强震序列进行系统性的预测效能评估.通过连续分段滑动、拟合和余震发生率预测,对地震序列的早期模型参数、余震预测限度及其制约因素进行了讨论,并利用N-test检验方法开展了分级分段分档预测结果的效能评估.通过对比1天和30天数据拟合参数的相关性发现,b、k、c具有较强的相关性,差别较小,而p值的相关性较弱,这或许和p值表示余震活动在长期内的衰减特性有关,Omi-R-J模型可在震后早期记录相对不完整的阶段更早的获得较为稳定的序列参数;利用Omi-R-J模型对中国大陆强震序列的预测结果显示,Omi-R-J模型对中国大陆早期的强余震具有较好的预测能力,总的预测有效率为81.16%,平均预测有效率为83.82%,预测过少的比例大于预测过多的比例,早期序列发育较差或监测能力有限的区域,其预测效果受到很大制约;地震序列参数k是影响地震数目的关键因素,而地震序列发育程度可能是影响震后0.05天地震预测效率的重要因素之一.

摘要： 对1975年以来内蒙古地区5级以上震序列资料分析整理,将独立的5级以上地震序划分为孤立型、主震—余震型和多震型,并分析了地震序列类型的空间分布特征。根据余震序列资料,统计了余震序列的几项基本参数,包括最大余震震级、序列持续时间、主震–最大余震时间间隔和余震区半径等,并且拟合了序列参数和震级的线性关系,分析了序列参数的频度分布特征,在此基础上总结出震后早期地震形势判定的规则。

摘要： 前震是短临地震预报最有效的方法之一.文中系统研究了1976—2017年全球29次MW≥8.0浅源地震的前震特征,发现:1)共有8次大地震出现前震,占总数的27.6％,且最大前震的震级M>5.0.这8次大地震均为逆冲型,占逆冲型地震总数(23次)的34.8％;2)前震序列的空间分布集中,且集中在主震震中周围.最大前震与主震的震中距为10～53km,震级差为1.1～2.8级.87.5％的时间差(最大前震与主震的发生时间差)为2h～15d,仅2006年11月15日千岛群岛8.3级地震前45d发生6.6级最大前震;3)与背景地震活动相比,前震具有高频次活动的特点.在8次大地震中,5次在主震前15d内加速活动;3次加速活动的开始时间超过1个月(35～45d),但在主震前1 d和6d又再次出现前震频次显著增多的现象;4)前震的震源机制解与主震一致,而余震的震源机制解比较复杂;5)使用ETAS模型计算前震和余震序列的参数 α、p和b值.为了确保参数计算的可靠性,文中约定当前震序列满足计算样本量N≥30且最低计算震级Mj≥Mc(Mc为最小完备震级)时,才同时计算前震和余震序列参数.此外,文中约定余震序列的持续时间为主震后1个月.在8次大地震中,有4个可以进行对比计算.结果显示,对于反映激发次级余震能力的 α值和序列衰减快慢的p值,前震和余震没有规律性差异;而反映应力水平的b值则共性特征明显,前震b值明显低于余震;与区域背景b值相比,前震b值低于区域背景b值的10％～24％.前震b值与区域背景b值之差是区域背景b值标准差的2.2～7.1倍,具有显著性.而余震b值高于区域背景b值或与其相当.为了讨论前震b值的稳定性,计算了2个资料丰富的前震序列的b值随样本量的变化.结果显示,在前震序列的开始阶段b值较低,之后逐渐增大,当计算样本量N≥70后,b值基本稳定.

摘要： 利用山东乳山流动观测台网资料,采用双差精定位方法重新确定2014-05以来山东乳山M3.2地震余震序列的分布;在此基础上,利用双差层析成像方法反演该地区的三维速度结构模型.结果 显示,乳山地区的余震序列整体呈现出圆形聚集特征,优势分布主要集中在3～8 km深度,乳山震群集中区的余震活动主要发生在破裂集中的区域,在5 km深度最为集中,说明乳山震群主要发生的区域为上地壳;研究区不同深度的P波速度变化明显.由于离已知的主断层有一定的距离,因此可以判断余震分布集中区可能是乳山断裂的分支或者是一条新的隐伏断裂.

摘要： 利用290个ML≥2.5的景谷6.6级地震序列事件的数字波形资料,计算得到其震源参数.结果显示,地震矩M0数值范围为1012～1016 N·m;拐角频率fc为3～16 Hz;视应力数值范围为0～15 MPa,平均值是2.14 MPa,中位数为1.39 MPa.同时,各个震源参数之间均存在一定的相关关系;景谷地震序列余震空间分布可分为北段、中段与南段,其视应力高低分布依次为:南段>中段>北段;表明主震发生后,震源区南段区域的应力水平最高,这可能是后续5.8级、5.9级强余震均发生在南段的根本原因.在视应力时间变化特征方面,5.8级、5.9级强余震的发生前,震中所在的南段视应力数值存在"逐渐变大-维持稳定"的特点.而后续未发生强余震的北段和中段,其视应力时间变化为早期快速调整,后续逐渐变化到各自区域的视应力均值水平;这类变化表明对于震源区南段,在强余震发生前的短时间内,南段区域首先出现局部应力不断挤压增强,并随后维持相对高值应力水平的"闭锁"状态,直到强余震的发生.因而,震后震源区余震视应力的空间分布对强余震可能产生地点判定有指示作用,其时间变化特征与地震序列中较强地震的发生存在一定的相关性.余震视应力的时空变化对序列后续地震活动的判定有重要指示作用,可为强余震的时空判定提供相应参考.

摘要： 利用290个ML≥2.5的景谷6.6级地震序列事件的数字波形资料,计算得到其震源参数。结果显示,地震矩M0数值范围为1012~1016 N·m;拐角频率fc为3~16 Hz;视应力数值范围为0~15 MPa,平均值是2.14 MPa,中位数为1.39 MPa。同时,各个震源参数之间均存在一定的相关关系;景谷地震序列余震空间分布可分为北段、中段与南段,其视应力高低分布依次为:南段>中段>北段;表明主震发生后,震源区南段区域的应力水平最高,这可能是后续5.8级,5.9级强余震均发生在南段的根本原因。在视应力时间变化特征方面,5.8级、5.9级强余震的发生前,震中所在的南段视应力数值存在“逐渐变大-维持稳定”的特点。而后续未发生强余震的北段和中段,其视应力时间变化为早期快速调整,后续逐渐变化到各自区域的视应力均值水平;这类变化表明对于震源区南段,在强余震发生前的短时间内,南段区域首先出现局部应力不断挤压增强,并随后维持相对高值应力水平的“闭锁”状态,直到强余震的发生。因而,震后震源区余震视应力的空间分布对强余震可能产生地点判定有指示作用,其时间变化特征与地震序列中较强地震的发生存在一定的相关性。余震视应力的时空变化对序列后续地震活动的判定有重要指示作用,可为强余震的时空判定提供相应参考。

摘要： 据中国地震台网(CENC)测定,2016年1月14日5时18分,在新疆维吾尔自治区巴音郭楞蒙古自治州轮台县发生M5.3地震,USGS也给出了这次地震的震源参数.截止2016年1月21日,新疆轮台单台共记录到余震序列462次,其中,ML2.0—2.9地震16次,ML3.0—3.9地震4次,ML4.0—4.9地震1次,最大余震为1月20日19时08分的ML4.5地震.本次地震的余震序列主要以ML 2.0以下地震为主，余震主要集中分布于震后3天，特别是震后7h，衰减系数P值为1.38，序列强度衰减不明显。此后地震频度衰减迅速，由每小时近60次迅速衰减至10次以下。序列b值与该区背景值相当，h值>1.0。综合序列参数、地震能量比以及震区周围历史地震序列情况分析认为，本次地震属于主震-余震型。初动和CAP方法计算得到的震源机制解都为逆断类型。

摘要： 中强地震后短时间内发生大量余震,这些余震的记录波形相互交叠干扰并受到主震面波尾波的影响,其中一些余震波形信噪较低,通过测震台网常规分析方法只能在震中距最近的一个台站检测到较清晰的P、S波震相到时,即为单台记录地震事件.中强地震余震序列的地震目录中存在大量单台记录地震事件.地震观测报告未给出其震中位置,且给出的震级可信度较差.精确估计单台记录地震事件的震中和震级,能够补充完善现有地震目录,为余震序列特征分析、震后趋势快速判断和强余震预测等研究提供重要的基础资料.

摘要： 本研究首先采用基于贝叶斯理论的绝对定位方法(Myers et al., 2007, 2009)对主震进行定位，获得主震的绝对位置为82.56°E, 36.04°N和12.3km深，其中经度、纬度和深度的定位误差分别为1.4km, 1.1km和3.6km，主震位于阿尔金断裂带4条分支断层的交汇处。余震序列的相对定位结果显示，余震优势展布方向为南西向，且呈现出纯单侧扩展模式，整个扩展仅22km左右;震源主要分布在5-10km深度范围内，属于浅源型地震，且发震断层面较为陡立。在主震区的北东方向还存在一个地震丛集，该区域为1924年民丰Ms7.3级双震的发震区，通过相同时间段内地震数量的比较，我们推测本次于田Ms7.3级主震可能触发了民丰主震区的地震活动，而并非余震向北东向扩展。为认识余震分布与主震破裂过程之问的关系，将定位结果投影到主震破裂过程反演结果(张勇等，2014)上，可以看出余震分布与主震破裂过程反演结果有较好的一致性，即余震主要位于滑动量较小区域，而在滑动量较大的主破裂区很少。如此的余震分布特征，可能表明在主破裂区积累的应力得到一定程度释放而难以产生地震，而在破裂的边缘由于受到阻碍，应力降减小，容易产生破裂引发余震(土卫民等，2005;赵翠萍等，2008)。主震破裂过程主要由主震处向西南方向破裂，与认为余震由主震处西南向单侧破裂具有较好一致性，也进一步证实了结果的可靠性。

摘要： 通过余震序列目录完整性分析得到了最小完整性震级随时间变化的表达式。主震后0.02天到主震后0.3天时间段内，最小完整性震级随时间的对数线性下降，在主震后0.3天以后最小完整性震级稳定在ML1.1左右。基于修正的大森公式，对2013年7月22日甘肃崛县-漳县MS6.6地震余震序列时间域衰减特征进行了分析。鉴于主震后短时间内目录遗漏的余震较多，首先利用崛县台连续地震记录波形的高频包络差，检测主震后3小时内目录遗漏的地震。经检测在主震后3个小时内共发现目录遗漏的ML1.0以上地震139个，其中主震后1000秒内检测到ML1.0以上遗漏地震67个，约为目录给出余震数量的6倍。遗漏地震参与分析对于利用修正的大森公式估计余震序列衰减特征有重要作用，通过余震频度拟合和地震矩拟合得到的余震衰减速率P值皆为1.07，表明崛县一漳县地震余震序列衰减速率处于全球平均水平。在分析余震序列衰减速率时，加入遗漏地震的余震频度修正的大森公式拟合所得结果较优;若无条件拾取目录遗漏地震，则使用地震矩修正的大森公式拟合得到的结果更为准确，但需充分考虑其误差。三种衰减模式拟合直接使用高频包络差计算地震频度曲线比较结果显示无法证明崛县-漳县地震主震后存在早期余震缺失现象。

摘要： 2011年6月20日,腾冲地区发生了5.2级地震,哈佛大学矩张量解结果表明,此次地震具有很大的逆冲分量,该地震是否与岩浆活动有关?本研究拟采用全矩张量解反演,来探讨此次地震的发震模式.对地震序列的重定位，联合使用了云南固定台站及流动台站资料，重定位结果显示，腾冲5.2级地震震源深度修正为12km，其余震序列在平面上没有优势的展布方向，而是呈现近似椭圆的分布;在深度上，余震序列呈现出明显的柱状形态，且向北西方向有倾斜趋势，而腾冲火山就位于该地震序列的北西方向，推测该序列与腾冲火山区热物质有密切关系。为此，采用在腾冲火山区的6个流动台站波形资料来反演 5.2级地震矩张量解。首先是对原始波形进行预处理，主要包括去均值、去除仪器响应、将记录变为实际速度，然后再旋转为ZRT方向。为了得到稳定可靠的矩张量结果，分两种方法反演:第一步反演全矩张量解，结果显示，双力偶(DC)部分占20%，补偿线性偶极子(CLOD)部分占62%，各向同性(ISO)部分占18%;第二步是反演偏矩张量解，结果显示，双力偶(DC)部分占25%，补偿线性偶极子(CLVD)部分占75%。用震源球来表示震源机制，主要包含震源取向的信息，比如走向、倾向及倾角。当考虑非双力偶型机制时，不考虑震源取向的震源类型表达方式就更为直观。为此，给出5.2级地震的“震源类型图”，结果显示，5.2级地震震源类型位于补偿线性偶极子部分。以上的结果表明，此次腾冲5.2级地震很可能为非双力偶型地震，结合余震序列的分布，其与腾冲火山岩浆的活动密切相关。

摘要： 北京时间2008年3月21日,新疆于田县境内发生Ms7.3地震.地震之后余震较为活动,同年10月5日,其西北800km的新疆乌恰县发生Ms6.8地震.本文计算了于田Ms7.3地震断层破裂在周围介质中产生的完全库伦破裂应力，将余震序列精定位结果投影在60.69s时的库伦破裂应力变化的空间分布图上。大部分余震分布在库伦破裂应力为正的区域，仅少量分布在边缘及负区，动态和静态库伦应力变化对大部分余震均有较明显的触发作用。计算了乌恰Ms6.8地震震中12km处的库伦应力变化。于田Ms7.3地震断层破裂在乌恰Ms6.8地震震中造成的最大动态库仑应力变化为0.67 MPa，其正的最大峰值超过了0.5 MPa的动态触发阈值，静态库仑应力变化为O.OOO5MPa，远远小于动态变化，也没有达到0.01 MPa的静态触发阈值。因此认为，乌恰Ms6.8地震可能主要受到了于田Ms7.3地震的动态应力触发作用，而静态库伦应力对其影响很小。

摘要： 本文通过联合和田台阵的方位角和新疆测震台网的震相到时数据对2008年新疆于田7.3级地震及序列进行了重新定位，重新定位后的结果与新疆区域地震台网定位结果有明显差异。新的定位结果表明，2008年于田7.3级地震序列震中分布呈近NS向条带分布，与对该地震野外考察发现的地表破裂带走向基本一致。通过对2008年于田7.3级地震及序列进行重新定位研究，认识到适合国内其它地区的定位方法，在新疆地区可能会遇到困难。因此，在精细分析2014年于田7.3地震序列震相数据的基础上，选用分层速度模型，通过联合新疆部分台站的方位角对2014年于田7.3级地震及其余震(截止时间为2014年4月30日，约300次地震事件)进行了重新绝对定位。重新定位结果显示，于Hi地震的主震震中位置为36.08°N, 82.59°E，位于主震西南的余震优势分布方向为北东东，与阿尔金断裂的走向一致。

摘要： 本文通过联合和田台阵的方位角和新疆测震台网的震相到时数据对2008年新疆于田7.3级地震及序列进行了重新定位，重新定位后的结果与新疆区域地震台网定位结果有明显差异。新的定位结果表明，2008年于田7.3级地震序列震中分布呈近NS向条带分布，与对该地震野外考察发现的地表破裂带走向基本一致。通过对2008年于田7.3级地震及序列进行重新定位研究，认识到适合国内其它地区的定位方法，在新疆地区可能会遇到困难。因此，在精细分析2014年于田7.3地震序列震相数据的基础上，选用分层速度模型，通过联合新疆部分台站的方位角对2014年于田7.3级地震及其余震(截止时间为2014年4月30日，约300次地震事件)进行了重新绝对定位。重新定位结果显示，于Hi地震的主震震中位置为36.08°N, 82.59°E，位于主震西南的余震优势分布方向为北东东，与阿尔金断裂的走向一致。

摘要： 本文通过联合和田台阵的方位角和新疆测震台网的震相到时数据对2008年新疆于田7.3级地震及序列进行了重新定位，重新定位后的结果与新疆区域地震台网定位结果有明显差异。新的定位结果表明，2008年于田7.3级地震序列震中分布呈近NS向条带分布，与对该地震野外考察发现的地表破裂带走向基本一致。通过对2008年于田7.3级地震及序列进行重新定位研究，认识到适合国内其它地区的定位方法，在新疆地区可能会遇到困难。因此，在精细分析2014年于田7.3地震序列震相数据的基础上，选用分层速度模型，通过联合新疆部分台站的方位角对2014年于田7.3级地震及其余震(截止时间为2014年4月30日，约300次地震事件)进行了重新绝对定位。重新定位结果显示，于Hi地震的主震震中位置为36.08°N, 82.59°E，位于主震西南的余震优势分布方向为北东东，与阿尔金断裂的走向一致。

摘要： 本文通过联合和田台阵的方位角和新疆测震台网的震相到时数据对2008年新疆于田7.3级地震及序列进行了重新定位，重新定位后的结果与新疆区域地震台网定位结果有明显差异。新的定位结果表明，2008年于田7.3级地震序列震中分布呈近NS向条带分布，与对该地震野外考察发现的地表破裂带走向基本一致。通过对2008年于田7.3级地震及序列进行重新定位研究，认识到适合国内其它地区的定位方法，在新疆地区可能会遇到困难。因此，在精细分析2014年于田7.3地震序列震相数据的基础上，选用分层速度模型，通过联合新疆部分台站的方位角对2014年于田7.3级地震及其余震(截止时间为2014年4月30日，约300次地震事件)进行了重新绝对定位。重新定位结果显示，于Hi地震的主震震中位置为36.08°N, 82.59°E，位于主震西南的余震优势分布方向为北东东，与阿尔金断裂的走向一致。

摘要： 本发明提供一种低余震的倒车雷达系统及输出低余震倒车雷达信号的方法，其包含一收发感测单元、一超声波发射电路、一超声波接收电路、一处理单元、一警报单元以及一存储单元；其中，该收发感测单元包含至少一压电陶瓷片；而该处理单元可发送一控制信号以及一脉冲信号，其中，该控制信号具有连续的正、负振幅相位，且该脉冲信号的起始振幅与该控制信号的截止振幅同相位，使该收发感测单元输出一低余震超声波信号。

摘要： 特征量提取装置(100)具有：数据取得部(110)，其取得时间序列数据作为输入时间序列数据(D1)；多滤波器应用部(120)，其具有多个数字滤波器，对数据取得部(110)取得的输入时间序列数据(D1)应用各数字滤波器，按照每个数字滤波器输出应用后的包含时间序列特征或频率特征的时间序列数据即滤波器响应时间序列数据(D5)；以及特征量提取部(130)，其针对多滤波器应用部(120)输出的多个滤波器响应时间序列数据(D5)，按照每个滤波器响应时间序列数据(D5)提取特征量，输出提取出的特征量作为特征量数据(D2)。

摘要： 本发明从RNA干扰的目标基因的碱基序列中搜索符合下述规则(a)－(d)的序列部位，根据搜索结果进行能够产生RNA干扰的小干扰RNA(siRNA)的设计、合成等。(a)3’末端的碱基是腺嘌呤(A)、胸腺嘧啶(T)或尿嘧啶(U)。(b)5’末端的碱基是鸟嘌呤(G)或胞嘧啶(C)。(c)3’末端的7碱基序列中富含由腺嘌呤(A)、胸腺嘧啶(T)和尿嘧啶(U)组成的群组中选择的1或1种以上的碱基。(d)碱基数在产生RNA干扰且不产生细胞毒性的范围内。

摘要： 本发明提供了一种基于断层缓冲带的余震序列删除方法，选取典型地震作为主震，结合余震时空分布和断裂走向，考虑地质构造的空间差异性，提出断裂缓冲带经验公式，使用断层缓冲带作为余震空间窗对K‑K删除余震方法的空间窗选取方式进行改进，实现了在不同震级地震中删除余震的功能，删除余震更可靠，删除的余震序列具有丛集性，有较强的地质理论解释。本发明解决了余震和断裂分布的关系，改善了K‑K法设置空间窗的局限性，调整了空间窗设定的合理性，较大程度地还原了该地区的地质特性，删除的余震序列与实际的地震发生分布有较好的切合度，提高了在不同震级地震中删除余震的效果。

摘要： 一种基于主余震序列的损伤指数评估方法，涉及一种基于主余震序列的损伤指数评估方法。本发明为解决现有地震损伤指数评估方法中只考虑主震单独作用，而忽略余震对结构的影响的问题。本发明按以下步骤进行：一、根据实际地震动类型选择一个结构滞回模型；二、选择单自由度结构的周期、质量、阻尼，并计算结构动力特性；三、根据选定的地震动和结构特性进行弹性分析，得到弹性强度需求Fe；四、选择强度折减系数值R，并计算出非弹性结构的屈服强度Fy；五、根据结构动力方法计算出该结构在地震作用下的反应，并根据损伤计算公式得到损伤指数DI；六、根据结构损伤指数DI的统计结果，提出主余震序列的损伤谱预测模型。本发明可应用于地震工程领域。

摘要： 一种基于主余震序列的能量谱预测方法，本发明涉及基于主余震序列的能量谱预测方法。本实发明是为了解决目前提出的能量谱模型均是基于主震单独作用，而忽略了余震对结构的影响的问题。具体是按照以下步骤进行的：步骤一：选取非弹性SDOF结构，并确定非弹性SDOF结构的动力特性；步骤二：根据步骤一中确定的非弹性SDOF结构的动力特性，求解基底受水平地震动的非弹性SDOF结构的运动方程得到EI；步骤三：根据步骤二得到EI，求得VE；步骤四：根据步骤三得到的VE，拟合得到基于主余震的能量谱模型。本发明根据结构的动力特性及其所在场地的主余震危险性水平，即可通过能量谱模型计算得到输入能大小。本发明应用于地震工程领域。

摘要： 本发明提供一种随机主余震序列的模拟及随机动力分析方法，属于地震的随机动力分析领域。关键内容包括基于Copula函数建立主余震参数相关性、基于随机物理模型产生地震序列，基于DPIM方法进行随机响应分析。本发明能够解决现有的随机地震动模型对余震在地震动时产生的影响的忽略，可以合理地从概率的角度建立主余震参数之间的随机关系，能够克服传统的结构动力分析方法计算量大、效率低的局限，从而极大地缩短了运算时长，提高了计算精度。

摘要： 本发明提供了一种基于断层缓冲带的余震序列删除方法，选取典型地震作为主震，结合余震时空分布和断裂走向，考虑地质构造的空间差异性，提出断裂缓冲带经验公式，使用断层缓冲带作为余震空间窗对K‑K删除余震方法的空间窗选取方式进行改进，实现了在不同震级地震中删除余震的功能，删除余震更可靠，删除的余震序列具有丛集性，有较强的地质理论解释。本发明解决了余震和断裂分布的关系，改善了K‑K法设置空间窗的局限性，调整了空间窗设定的合理性，较大程度地还原了该地区的地质特性，删除的余震序列与实际的地震发生分布有较好的切合度，提高了在不同震级地震中删除余震的效果。

摘要： 一种基于主余震序列的能量谱预测方法，本发明涉及基于主余震序列的能量谱预测方法。本实发明是为了解决目前提出的能量谱模型均是基于主震单独作用，而忽略了余震对结构的影响的问题。具体是按照以下步骤进行的：步骤一：选取非弹性SDOF结构，并确定非弹性SDOF结构的动力特性；步骤二：根据步骤一中确定的非弹性SDOF结构的动力特性，求解基底受水平地震动的非弹性SDOF结构的运动方程得到EI；步骤三：根据步骤二得到EI，求得VE；步骤四：根据步骤三得到的VE，拟合得到基于主余震的能量谱模型。本发明根据结构的动力特性及其所在场地的主余震危险性水平，即可通过能量谱模型计算得到输入能大小。本发明应用于地震工程领域。

摘要： 一种基于主余震序列的损伤指数评估方法，涉及一种基于主余震序列的损伤指数评估方法。本发明为解决现有地震损伤指数评估方法中只考虑主震单独作用，而忽略余震对结构的影响的问题。本发明按以下步骤进行：一、根据实际地震动类型选择一个结构滞回模型；二、选择单自由度结构的周期、质量、阻尼，并计算结构动力特性；三、根据选定的地震动和结构特性进行弹性分析，得到弹性强度需求Fe；四、选择强度折减系数值R，并计算出非弹性结构的屈服强度Fy；五、根据结构动力方法计算出该结构在地震作用下的反应，并根据损伤计算公式得到损伤指数DI；六、根据结构损伤指数DI的统计结果，提出主余震序列的损伤谱预测模型。本发明可应用于地震工程领域。