西藏地区地震类型特征研究

对1970年以来西藏地区139组5级以上地震序列类型进行统计及特征分析,发现西藏地区中强以上地震序列以主—余型和孤立型为主,约占79.1%。随主震震级的增加,主—余型所占比例逐渐增加,孤立型和多震型所占比例减少。绝大多数地震的最大震级余震在主震后3个月内发生。西藏地区多震型地震主要发生西藏中强以上地震活跃时段,因此在西藏及邻区5级以上地震频度较高时需注意多震型地震的发生。从地震序列类型空间分布来看,孤立型地震和主—余型地震大多发生在构造相对单一的地区,而多震型地震则大多发生在构造复杂且存在多组构造交汇的区域。     

最不利主余震序列型地震动作用下RC框架结构易损性研究

余震的发生会造成结构的累积损伤,不同类型的主余震序列地震动对结构的影响有所差异。鉴于此,以主震卓越周期小于或接近余震卓越周期为基本原则,确定了最不利主余震序列地震动。选择4层RC框架结构为研究对象,在增量动力分析的基础上,定义了4个性能水平,以此来研究该结构的易损性。依据破坏状态概率和震害指数,得到7度多遇、7度设防和7度罕遇地震的易损性指数。研究结果表明:依据我国规范设计的RC框架结构的4个性能水平的量化指标限值依次为1/495、1/263、1/108和1/45,余震的出现会加剧结构的破坏状态。当以易损性指数作为评价指标时,可认为该结构能够满足小震不坏、中震可修和大震不倒的抗震设防目标。     

A Review of Earthquake Statistics: Fault and Seismicity-Based Models, ETAS and BASS

There are two fundamentally different approaches to assessing the probabilistic risk of earthquake occurrence. The first is fault based. The statistical occurrence of earthquakes is determined for mapped faults. The applicable models are renewal models in that a tectonic loading of faults is included. The second approach is seismicity based. The risk of future earthquakes is based on the past seismicity in the region. These are also known as cluster models. An example of a cluster model is the epidemic type aftershock sequence (ETAS) model. In this paper we discuss an alternative branching aftershock sequence (BASS) model. In the BASS model an initial, or seed, earthquake is specified. The subsequent earthquakes are obtained from statistical distributions of magnitude, time, and location. The magnitude scaling is based on a combination of the Gutenberg-Richter scaling relation and the modified Båth’s law for the scaling relation of aftershock magnitudes relative to the magnitude of the main earthquake. Omori’s law specifies the distribution of earthquake times, and a modified form of Omori’s law specifies the distribution of earthquake locations. Unlike the ETAS model, the BASS model is fully self-similar, and is not sensitive to the low magnitude cutoff.     

Characteristics of the May 12,2008 Wenchuan M_S8.0 Earthquake Sequence and Discussion on Relevant Problems

The aftershock activity of the May 12, 2008 Wenchuan Ms8.0 Earthquake Sequence shows an obvious segmented feature. Most of the large aftershocks were distributed in the north and south parts of the aftershock zone. Thrusting was dominant with a small amount of strike-slip component in the south part. The aftershock activity decayed gradually, presenting the sequence features of a mainshock-aftershock pattern. The north part was the ending area of the maiushock fracture where strike-slipping was dominant, showing an obvious swarm feature. Therefore it became the major area for large aftershocks. The modulation of the earth tide on aftershock activity is remarkable; most large aftershocks occur during the period of flood and neap tide. The time period around 16:00 was the dominant occurring time for large aftershocks. The p-value, a parameter of modified Omori formula, increases gradually with time, and reaches about 1 at the end. Based on previous study, the sequence patterns, magnitude of maximum aftershock, as well as the duration of aftershock activity has been discussed. The primary results also show that the magnitude difference between the maiushockand the maximum aftershock is proportional to the rupture size of the maiushock for huge earthquakes of about Ms8.0. This means that when the magnitudes of the earthquakes are nearly the same, large rupture size corresponds to sufficient energy release.     

A quantitative study of some aftershock sequences in greece

Temporal features of the aftershock activities following twelve moderate shallow earthquakes in Greece have been studied quantitatively, by making use of the modified Omori's formula and Akaike's Information Criterion (AIC).Although the earthquake sequences occurred in different seismotectonic regions in Greece, similar seismic patterns before the occurrence of large aftershocks have been observed. Aftershock activity shows an appreciable decrease from the level expected from the modified Omori formula before the occurrence of the large aftershock. Also in some cases, the aftershock activity recovers to a normal level or increases beyond it prior to the occurrence of the large aftershock, if the aftershock activity is monitored immediately after the occurrence of the main shock. In such case, the observed pattern may be useful in predicting a large aftershock.     

Bayesian estimation of the Modified Omori Law parameters for the Iranian Plateau

The forecasting of large aftershocks is a preliminary and critical step in seismic hazard analysis and seismic risk management. From a statistical point of view, it relies entirely on the estimation of the properties of aftershock sequences using a set of laws with well-defined parameters. Since the frequentist and Bayesian approaches are common tools to assess these parameter values, we compare the two approaches for the Modified Omori Law and a selection of mainshock–aftershock sequences in the Iranian Plateau. There is a general agreement between the two methods, but the Bayesian appears to be more efficient as the number of recorded aftershocks decreases. Taking into account temporal variations of the b-value, the slope of the frequency-size distribution, the probability for the occurrence of strong aftershock, or larger main shock has been calculated in a finite time window using the parameters of the Modified Omori Law observed in the Iranian Plateau.     

考虑主余震序列作用的砌体结构易损性分析

鉴于我国砌体结构抗震面临的严峻形势,以及主余震序列作用下砌体结构易损性研究存在的不足,以汶川地震中1栋典型砌体房屋为依托,开展考虑主余震序列作用的砌体结构易损性研究。首先,建立了结构的三维有限元模型,并从结构动力特性与损伤等方面验证了模型的合理性; 其次,选取整体损伤耗能作为结构性能参数,并给出了对应不同破坏等级的结构整体损伤耗能界限值; 最后,基于IDA方法对结构进行了不同主余震序列作用下的易损性分析。结果表明:相对于层间位移角,整体损伤耗能更适用于描述余震对结构的累积损伤效应; 随着结构主震损伤程度的加深,余震对结构造成的影响越明显,且主余震序列作用下结构呈现出比主震单独作用下高一个破坏等级的趋势; 当ΔPGA小于0.6时,结构的极限状态超越概率最大增幅不超过10%,此时余震的影响较小; 当ΔPGA大于0.6时,结构的极限状态超越概率最大增幅可达到33.1%,此时余震的影响较大。     

A Damage Mechanics Model for Aftershocks

— All earthquakes are followed by an aftershock sequence. A universal feature of aftershock sequences is that they decay in time according to the modified Omori’s law, a power-law decay. In this paper we consider the applicability of damage mechanics to earthquake aftershocks. The damage variable introduced in damage mechanics quantifies the deviation of a brittle solid from linear elasticity. We draw an analogy between the metastable behavior of a stressed brittle solid and the metastable behavior of a superheated liquid. The nucleation of microcracks is analogous to the nucleation of bubbles in the superheated liquid. In this paper we obtain a solution for the evolution of damage after the instantaneous application of a constant strain to a rod. We show that the subsequent stress relaxation can reproduce the modified Omori’s law. It is argued that the aftershocks themselves cause random fluctuations similar to the thermal fluctuations associated with phase transitions.     

1976年唐山MS7.8地震余震序列持续时间及对地震危险性分析的意义

1976年7月28日唐山MS7.8大地震对唐山及其周边地区造成了重大的人员伤亡和财产损失. 主震之后约15小时滦县又发生了MS7.1地震; 同年11月15日宁河也发生了MS6.9地震. 唐山MS7.8主震后的余震一直持续至今,使该区域至今保持了与主震前相比具有较高的地震活动性．如何估计余震的持续时间,并进一步将余震从主震目录中去除,一直是地震学中所关注的问题．该文通过对数线性回归和理论计算,从不同角度求取并讨论了1976年唐山MS7.8大地震的余震持续时间．结果表明,由对数线性回归计算得到的余震持续时间约为80 a．而基于Dieterich的余震触发理论所得到的余震持续时间则与区域剪应力变化率有关．区域剪应力变化率可有几种不同方法求得: ① 根据剪应力变化率和静态应力降Delta;tau;e及地震回复周期tr之间的关系求取应力变化率,该方法所得到的余震持续时间约为70——100 a;② Ziv和Rubin对Dieterich的方法进行了修正,给出了通过远场加载速率和断层宽度求取应力变化率, 该方法得到的余震持续时间约为80 a;③ 由背景场地震活动性求取远场剪应力速率, 可以得到该地区二维分布式的余震持续时间,此方法得到的研究区域内余震持续时间为130——160 a．综上,唐山地区余震持续时间约为70——140 a,据此, 该地区现今所发生的地震仍为MS7.8唐山地震所触发的余震．      

基于经验模型和物理模型研究2013 M_S7.0芦山地震余震序列

基于中国地震台网中心2013 MS7.0芦山地震余震数据我们首先确定了余震空间分布范围并根据G-R关系计算了主震后半小时内的完备震级Mc=3.5,并且得到了ML≥3.5和ML≥3.0的地震在2001年至芦山地震前的背景场地震发生率.通过Omori-Ustu经验定律和两种Dieterich模型对芦山地震余震发生率的拟合,我们发现阶梯型Dieterich模型只能模拟p=1的情况,从而造成了模拟曲线与观测数据的差别;前人研究表明震后滑移同样是产生余震的原因,如果假设余震序列由主震静态剪应力Δτ和震后滑移共同作用所产生,我们数值模拟得到的对数型Dieterich模型能够较好地推断余震发生率R随时间t增加而衰减的趋势,能够从物理机制上解释MS7.0芦山地震余震序列衰减指数大于1这一现象.通过对数型Dieterich模型的拟合并结合Andrews的方法,我们还得到MS7.0芦山地震Aσ约为0.155 MPa,ta约为8.4年,这一值与前人研究结果十分接近.     

2013年甘肃岷县-漳县MS6.6地震余震 序列目录完备性研究 基于对单台记录地震事件震中与震级的估计

中强地震余震序列地震目录编目是否完备、 震源参数是否准确,直接影响余震序列特征分析、 震后趋势快速判断和强余震预测等研究结果的科学性和可靠性. 2013年7月22日甘肃岷县-漳县M_S6.6地震余震序列目录中存在较多单台记录地震事件,地震观测报告仅给出其震级,而未给出震中位置. 由于余震波形间的相互交叠干扰,使得余震最大振幅的测量误差较大,造成地震观测报告给出的单台事件震级误差较大. 精确估计单台记录地震事件的震中和震级,能够补充完善现有地震目录,提高地震目录的完备性. 本文对单台记录地震事件震中和震级的估计不仅限于单个台站,而是通过分析区域台网中多个台站的波形记录实现. 首先以余震序列中震级较大、 波形记录信噪比较高的地震波形作为模板,使用波形互相关震相检测技术,检测单台记录的地震事件在多个台站的震相到时. 如果能在4个以上台站检测到震相,则利用测震台网常用的HYPOSAT方法估计其震中位置,并利用多个台站记录波形与模板地震的振幅比估计其震级. 之后计算主震发生后不同时间的最小完备性震级,并通过线性拟合得到最小完备性震级随时间变化的表达式,以分析此地震余震序列的目录完备性. 经过计算共得到253个单台记录地震事件的震级和其中177个事件的震中位置,其震中空间分布范围与余震序列中其它地震分布范围基本一致. 震级复测以及与人工拾取震相到时误差对比表明,该方法所得震相检测和震级估计结果具有较好的可靠性. 主震及最大余震发生后的短时间内,有较多数量单台事件的目录所给出的震级偏低,分析认为可能受主震与较大余震后续震相以及余震间相互干扰所致. 主震发生0.02—0.3天内,其余震序列最小完备性震级随时间的对数呈线性下降,在0.3天后最小完备性震级稳定在M_L1.1左右.      

考虑主余震作用下的核岛厂房结构易损性研究

一次强震常伴随着多次余震作用,且时间间隔较短.研究表明,主震诱发的系列余震会对原有结构造成进一步的累积损伤.核岛厂房结构作为核反应堆的最后一道防线,其在服役期间可能遭受主震和余震的累积作用.参考我国现行规范,定义了四类性能水准(PL)和三种极限状态(LS),以混凝土最大应变作为结构破坏指数(DI),以谱加速度Sa作为地...     

全球大陆7级浅源大地震强余震震级和空间分布特征

本文收集了1967—2007年全球大陆(不包括中国大陆)29次7.0~7.9级浅源大地震的余震资料。根据震源机制解结果把这29次地震序列分为走滑型和非走滑型, 其中17次是走滑型, 12次是非走滑型, 并分别研究了走滑型和非走滑型地震序列强余震震级分布特征和空间分布特征。研究结果表明, 强余震与主震震级差服从指数分布, 统计得到了走滑型地震序列B值为0.58, 非走滑型地震序列B值为1.07。走滑型地震序列强余震与主震距离的优势分布范围是10～45 km, 而非走滑型地震序列强余震与主震的距离优势分布为20～59 km, 并且强余震与主震震中距服从正态分布。     

中国大陆中强地震余震序列的部分统计特征

依据1970年以来记录相对完整的294次50级以上地震序列资料,研究中国大陆中强地震余震序列统计特征,探讨序列类型、最大余震震级、强余震活动持续时间等与主震震级及主震断层性质之间的关系.中国大陆孤立型、主余型及多震型地震余震序列分别约占23%、59%及18%.其中走滑型、具有倾滑分量的走滑型、具有走滑分量的倾滑型及逆冲型分别占48%、24%、17%及11%. 余震序列1年内最大余震震级与主震震级正相关,但主震震级较低时相对离散,孤立型序列离散程度较高,主余型及多震型序列线性相关性较好.绝大多数序列最大余震均发生在震后200天内,少数具有晚期强余震的序列主要属主余型序列,孤立型及多震型序列通常没有晚期强余震发生.68%的序列1年内最大余震发生在震后10天内,77%发生在震后30天之内,95%发生在震后120天之内.序列最大余震发生时间及5、6级强余震活动持续时间与序列类型及主震震级大小有关,多震型序列最大余震发生最快,孤立型次之、主余型最长.若仅就主余型序列而言,当主震震级较高时最大余震与主震间时间间隔相对较短,主震震级较低时最大余震与主震间时间间隔相对较长.     

On the applicability of the RETAS model for forecasting aftershock probability in underground mines (Kiirunavaara Mine,Sweden)

Aftershock series of even comparatively small seismic events can pose a risk to the mining operation or the personnel in deep underground mines as the main shocks and some of the aftershocks can cause damage in the rock mass. Stochastic modeling was applied in this study for the analysis of the temporal evolution of aftershock occurrence probability during a M1.85 aftershock sequence in Kiirunavaara Mine, Sweden. The Restricted Epidemic-Type Aftershock Sequence (RETAS) model was chosen for estimation of the aftershock occurrence probability. This model considers all events with magnitude above the magnitude of completeness M₀ and has the advantage of including the Modified Omori Formula (MOF) model and Epidemic-Type Aftershock Sequence (ETAS) model as its end versions, considering also all intermediate models. The model was applied sequentially to data samples covering cumulative periods of time, starting from the first 2 h after the main event and increasing them by 2 h until the period covered the entire 72-h sequence. For each sample, the best-fit RETAS version was identified and the probability of a M?≥?0.5 aftershock for every next 2 h was determined through Monte Carlo simulation. The feasibility of the resulting probability evolution for suspension and re-starting of the mining operations was discussed together with possible prospects for future development of the methodology.

     

Preliminary report of the September 5, 2022 MS 6.8 Luding earthquake,Sichuan, China

The 2022 M_S 6.8 Luding earthquake is the strongest earthquake in Sichuan Province, Western China, since the 2017 M_S 7.0 Jiuzhaigou earthquake. It occurred on the Moxi fault in the southeastern segment of the Xianshuihe fault, a tectonically active and mountainous region with severe secondary earthquake disasters. To better understand the seismogenic mechanism and provide scientific support for future hazard mitigation, we summarize the preliminary results of the Luding earthquake, including seismotectonic background, seismicity and mainshock source characteristics and aftershock properties, and direct and secondary damage associated with the mainshock. The peak ground displacements in the NS and EW directions observed by the nearest GNSS station SCCM are ～35 mm and ～55 mm, respectively, resulting in the maximum coseismic dislocation of 20 mm along the NWW direction, which is consistent with the sinistral slip on the Xianshuihe fault. Back-projection of teleseismic P waves suggest that the mainshock rupture propagated toward south-southeast. The seismic intensity of the mainshock estimated from the back-projection results indicates a Mercalli scale of VIII or above near the ruptured area, consistent with the results from instrumental measurements and field surveys. Numerous aftershocks were reported, with the largest being M_S 4.5. Aftershock locations (up to September 18, 2022) exhibit 3 clusters spanning an area of 100 km long and 30 km wide. The magnitude and rate of aftershocks decreased as expected, and the depths became shallower with time. The mainshock and two aftershocks show left-lateral strike-slip focal mechanisms. For the aftershock sequence, the b-value from the Gutenberg-Richter frequency-magnitude relationship, h-value, and p-value for Omori’s law for aftershock decay are 0.81, 1.4, and 1.21, respectively, indicating that this is a typical mainshock-aftershock sequence. The low b-value implies high background stress in the hypocenter region. Analysis from remote sensing satellite images and UAV data shows that the distribution of earthquake-triggered landslides was consistent with the aftershock area. Numerous small-size landslides with limited volumes were revealed, which damaged or buried the roads and severely hindered the rescue process.     

The relationships between earthquakes and positions of the sun and moon（Ⅱ）——Sometemporalcharacteristicsoftheaftershocksequencesofstrongearthquakes

Ｔｈｅｒｅｌａｔｉｏｎｓｈｉｐｓｂｅｔｗｅｎｅａｒｔｈｑｕａｋｅｓａｎｄｐｏｓｉｔｉｏｎｓｏｆｔｈｅｓｕｎａｎｄｍｏｏｎ（Ⅱ）——Ｓｏｍｅｔｅｍｐｏｒａｌｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｔｈｅａｆｔｅｒｓｈｏｃｋｓｅｑｕｅｎｃｅｓｏｆｓｔｒｏｎｇｅ...     

自维持断层蠕滑过程和余震发生率的衰减特征—以汶川MW7.9地震序列为例

基于Dieterich地震活动性理论,本文推导出计算余震发生率和余震累积次数的一般表达式,其中主震后发震断层内部的剪切应力随时间的演化过程遵从Jeffreys-Lomnitz蠕变模型,且与修正Omori定律直接相关。修正Omori定律中的p值与震后断层的短时应力加卸载过程正相关。采用Rubin和Ampuero 给出的震后断层自维持蠕滑模型本文得出计算余震发生率的近似表达式,并对2008年汶川地震序列进行拟合。结果表明,p值的大小直接对应了速率-状态摩擦定律中摩擦参量b/a,而修正Omori定律中的c值则与速率-状态摩擦定律中的临界滑移D_c相关。对于汶川余震序列而言,拟合结果显示b/a约为1.13,D_c约为2—3 cm。Rubin-Ampuero震后自维持蠕滑描述了震后孕震层内部短暂的速率变化特征,是孕震断层演化过程不可缺少的环节。      

Aftershock Statistics of the 1999 Chi–Chi, Taiwan Earthquake and the Concept of Omori Times

In this paper we consider the statistics of the aftershock sequence of the m = 7.65 20 September 1999 Chi–Chi, Taiwan earthquake. We first consider the frequency-magnitude statistics. We find good agreement with Gutenberg–Richter scaling but find that the aftershock level is anomalously high. This level is quantified using the difference in magnitude between the main shock and the largest inferred aftershock $ {{\Updelta}}m^{ *}. $ Typically, $ {{\Updelta}}m^{ *} $ is in the range 0.8–1.5, but for the Chi–Chi earthquake the value is $ {{\Updelta}}m^{ *} $  = 0.03. We suggest that this may be due to an aseismic slow-earthquake component of rupture. We next consider the decay rate of aftershock activity following the earthquake. The rates are well approximated by the modified Omori’s law. We show that the distribution of interoccurrence times between aftershocks follow a nonhomogeneous Poisson process. We introduce the concept of Omori times to study the merging of the aftershock activity with the background seismicity. The Omori time is defined to be the mean interoccurrence time over a fixed number of aftershocks.     

Aftershock sequence simulations using synthetic earthquakes and rate-state seismicity formulation

We use an efficient earthquake simulator that incorporates rate-state constitutive properties and uses boundary element method to discretize the fault surfaces, to generate the synthetic earthquakes in the fault system. Rate-and-state seismicity equation is subsequently employed to calculate the seismicity rate in a region of interest using the Coulomb stress transfer from the main shocks in the fault system. The Coulomb stress transfer is obtained by resolving the induced stresses due to the fault patch slips onto the optimal-oriented fault planes. The example results show that immediately after a main shock the aftershocks are concentrated in the vicinity of the rupture area due to positive stress transfers and then disperse away into the surrounding region toward the background rate distribution. The number of aftershocks near the rupture region is found to decay with time as Omori aftershock decay law predicts. The example results demonstrate that the rate-and-state fault system earthquake simulator and the seismicity equations based on the rate-state friction nucleation of earthquake are well posited to characterize the aftershock distribution in regional assessments of earthquake probabilities.