Seismic Hazard Assessment for the Main Seismogenic Zones in North Algeria

—A seismogenic zoning for north Algeria is proposed. The characteristics of each zone are studied in detail and different statistical analyses are performed in order to assess the slope of the Gutenberg-Richter relationship, the mean annual activity rate of earthquakes, the average return period, the probability of non-exceedence during a time period and the maximum expected magnitude. The slope of the Gutenberg-Richter relationship determined by traditional methods, physical strain energy release and corresponding magnitude, together with a new approach based on both historical and instrumental data, are used to estimate the seismic hazard parameters.¶The seismic hazard parameters obtained described well the characteristics of each proposed seismogenic zone. The b-value of the Gutenberg-Richter relationship concurs especially well with the most common seismotectonic interpretation of each zone, for example 0.70 ± 0.06 at El Asnam region and 0.63 ± 0.09 around Algiers. The average return period obtained (approximately 300 years for magnitude 7.5 at El Asnam region, exceeding 1000 years in other regions), corroborate those obtained through Paleoseismic studies and others using a digital elevation model.     

地震b值的估计方法及其标准误差——应用蒙特卡罗方法估计b值精度

近年在地震预报问题的探索中, 不少人企图使用震级频数法则去寻找大震前兆, 讨论 b 值的时空变化.随之出现一个新问题, 为适应预报需要资料所取范围越来越小, 时段越来越短, 这样估计 b 值时使用的地震个数大大减少, 至使估计 b 值时误差增加.本文从古登堡、李希特的震级频数法则出发, 讨论了 b 值的各种估计方法.其中有极大似然估计、矩估计, 线性简单最小二乘估计和非线性简单最小二乘估计等.给出 b 值各种估计的精度, 从而比较各种方法的优劣.文中应用蒙特卡罗(Monte-Carlo)方法计算了一个 b 值误差表, 供计算 b 值时参考.      

新疆南部构造区带与地震活动状态研究

Based on the studies of earthquake activity, tectonic movement, crustal shortening rate, fault activity, local stress field and historical characteristics of strong earthquake activities in Xinjiang, we divide the south part of Xinjiang into 4 seismotectonic zones, namely, the eastern segment of south Tianshan seismic belt, the Kalpin block, the Kashi-Wuqia junction zone, and the west Kunlun Mountains seismic belt. Using earthquake catalogues from Xinjiang since 1900, and on the basis of integrity analysis of earthquake records in different magnitude ranges, the seismicity state of different seismotectonic zones is analyzed quantificationally by calculating the mean value of annual strain energy release, annual rate of earthquakes with different lower limits of magnitude, b-value, and the parameter m of accelerating strain release model. The characteristic indexes of seismicity state for each of the seismic tectonic zones are then determined, which provide a quantitative basis for earthquake tendency analysis and judgment.     

相关平差在应变分析中的应用

研究了直接用多期重复大地测量观测资料求剪切应变速率的相关平差方法,简单地讨论了判断计算结果中反映形变信息的可靠性的标准,并对相关平差法与一般方法计算结果的差别进行了比较。通过几个地区的实际计算,认为在地震工作中,为探测大范围的地壳水平运动,迫切需要大幅度提高目前的大地测量精度。      

Co-Seismic Groundwater Level Changes Induced by the May 12, 2008 Wenchuan Earthquake in the Near Field

The large scales of co-seismic water level changes in mainland China were observed in response to the tragic 2008 Ms 8.0 Wenchuan earthquake. To better understand the mechanism of these hydrogeological phenomena, groundwater-level data at 17 confined wells, with an epicentral distance of <500 km, were collected. We compare the static strain predicted by dislocation theory with the volumetric strain calculated by the tide effect of the groundwater based on poroelastic theory. The results show that the sign of the co-seismic groundwater level change is consistent with the sign predicted by dislocation theory. Additionally, the magnitude of the strain calculated by the two methods is also concordant in half of the wells. In the rest of the wells, the strains inversed from the groundwater level are one or two orders of magnitude larger than the fault dislocation model. These wells mostly have an epicenter distance larger than 300 km; therefore, the dynamic stress induced by the seismic wave may be responsible for the co-seismic water level changes in these wells. According to these results, we roughly estimate that the effect range of the static stress is approximately 300 km for the Wenchuan earthquake, and the dynamic stresses dominate beyond this epicenter distance. In addition, geological and hydrogeological conditions and other mechanisms may be responsible for these changes.     

构造地震及其可能震级的初步研究——以“人”字型构造为例

一个地区未来地震的可能震级的估计是地震危险性评价的一个极其重要的方面。本文根据已知的地质条件,利用线弹平面有限元方法计算区域构造应力场,并确定所积累的总应变能和产生的破裂区面积,然后利用经验公式估计出与断层活动有关的构造地震的可能震级     

research of source parameters and stress state in shandong segment of tanlu fault zone

In recent years, there have been few researches and analysis published on the seismic activity and stress state in Shandong segment of Tanlu fault zone using digital seismological methods such as seismic apparent stress, focal mechanism solution and so on. In this paper, source parameters such as focal mechanism solutions and apparent stress are calculated using the waveform data of M_L≥1 moderate-small earthquakes in Shandong segment of Tanlu fault zone recorded by Shandong digital seismic network since 2007. According to focal mechanism solutions, a statistical analysis is done on the focal dislocation types in the study area using triangle graphical method, and the results show that the faulting in this area is mainly of strike-slip mechanism, and there are less thrust and normal mechanism. Calculation with the mean stress tensor method illustrates that the direction of mean principle stress of Shandong segment of Tanlu fault zone is NEE-SWW, which is the result of the combined effect of the subduction of West Pacific plate and the extrusion of Indian plate to Eurasian plate; the small dip angle indicates that the mode of action of stress is nearly horizontal, and the direction of principal stress axis is nearly perpendicular to the Tanlu fault zone. Under the action of such compressive stress field, dislocation is not likely to occur and the stress accumulation is enhanced on both sides of the fault. The apparent stress is calculated using the source spectral parameters method. Apparent stress has positive correlation with the magnitude and increases with the increased magnitude. So we get apparent stress difference by subtracting the empirical fitting value from the apparent stress. By removing the impact of magnitude, and according to the temporal-spatial evolution image of apparent stress difference, we found that the apparent stress in Shandong segment of Tanlu fault zone generally has a trend of decrease starting from the Wenchuan earthquake in 2008, and the spatial distribution of apparent stress in the region is very uneven. Combined with the spatial distribution of b values, the result shows that high stress is mainly located in Anqiu segment and Tancheng-Juxian segment, especially in Anqiu segment where small magnitude earthquakes appeared accompanying with the high stress. Low b-value means high stress and low frequency means low stress release, which indicates that Anqiu segment might accumulate higher stress and is at the fault locking stage. The research will provide new data for better understanding the present active feature and stress state of the Shandong segment of the Tanlu fault zone.     

不同应力途径下岩石的体积变化

本文介绍了对岩石在不同应力途径下的体积变化进行的实验研究。主要的应力途径有三种,并且都是首先把岩石加载到破裂前一定的应力状态,再分别三种类型加载。实验用电阻应变片测量纵向和横向应变以计算体应变。实验表明:同一种岩石在不同应力途径下膨胀效应有显著的不同。文中用辉长岩的实验数据对过密量和超膨胀量进行了估算。研究表明,过密与膨胀是应力途径变化时岩石体积变化的滞后效应。研究了B型实验过密状态出现的条件。从强度80%左右开始减小最小主应力的B型实验明显表现经历过密状态。从低应力状态(强度50%左右)和高应力状态(强度92%以上)开始减小最小应力的B型实验,岩石均未经历过密状态。最后联系地震预报问题进行了简单的讨论     

怎样正确计算b值

在我们研究震级一频度客观分布的基础上,根据数理统计理论和大量的实际地震活动资料,对b值的正确计算问题进行了详细的讨论,结果认为最小二乘法是计算b值的最好方法。此外,还认真地研究了与b值正确计算有关的下限震级、上限震级、样本数以及震级分档区间等具体问题,并给出了它们的合理选择方法。     

Coupling between different layers of the lithosphere under horizontal drag underneath

Based on a layered rheological model of the lithosphere, the velocity and stress distributions in the lithosphere under horizontal drag underneath were calculated using viscoelastic finite element method of plain strain with finite deformation. In the simulation, different conditions of drag and blocking were assumed to study their influences on the stress distribution and the coupling between different layers. Blocking depth has little influence on the stress level in the whole area and the coupling between different layers, but influences the stress state in the area around the blocking. The area covered by the high stress anomaly becomes larger when the blocking depth becomes deeper, but the magnitude of the value of the maximum shear stress decreases. The greater the viscosity differences between different layers of the lithosphere, the greater the possibility of decoupling between them. Under the drag of normal mantle convection (the convection velocity is about 20 cm·a?1), a lithosphere with a rheological structure similar to that of North China could not have decoupling between different layers, while could have stress distribution with magnitude of several MPa to tens of MPa and could have anomalous areas with stress accumulation if the geological structure is complicated.     

Coupling between different layers of the lithosphere under horizontal drag underneath

Based on a layered Theological model of the lithosphere, the velocity and stress distributions in the lithosphere under horizontal drag underneath were calculated using viscoelastic finite element method of plain strain with finite deformation. In the simulation, different conditions of drag and blocking were assumed to study their influences on the stress distribution and the coupling between different layers. Blocking depth has little influence on the stress level in the whole area and the coupling between different layers, but influences the stress state in the area around the blocking. The area covered by the high stress anomaly becomes larger when the blocking depth becomes deeper, but the magnitude of the value of the maximum shear stress decreases. The greater the viscosity differences between different layers of the lithosphere, the greater the possibility of decoupling between them. Under the drag of normal mantle convection (the convection velocity is about 20 cm · a^?1), a lithosphere with a Theological structure similar to that of North China could not have decoupling between different layers, while could have stress distribution with magnitude of several MPa to tens of MPa and could have anomalous areas with stress accumulation if the geological structure is complicated.     

地震震级－频度关系b值参数的稳健估计

文中依据古登堡－里克特（Ｇ－Ｒ）震级－频度公式，结合实例说明了参数ｂ值的稳健估计方法，并将其与最小二乘法及最大似然法比较。结果表明，一般情况下，稳健估算结果与实际数据分布较吻合，并且可以在不剔除低震级段和大震级段的偏离分布值的情况下，取得与选取合理震级上下限后最小二乘法一致的结果。此外，稳健估计结果的余差分析同样也得到地震震级－频度关系不是简单的对数线性关系的结果，而原始Ｇ－Ｒ公式需加二次以上的高次项才能较好地描述该特性。     

福建沿海、台湾海峡GPS观测分析及地球动力学特征研究

利用 3期GPS联测结果所获得的福建沿海地壳水平运动信息 ,采用ITRF94全球框架为基础的GPS测站地壳运动模型及其处理软件 ,对所获得的观测数据进行处理和精度分析。得到福建省高精度的GPS测站大地坐标、边长及其位移矢量 ,其精度达到 1 7×10 - 8。计算了福建地壳运动速率、主应变率 ,东西与南北向线应变率、面应变率、剪应变率、大地转动率和最大剪应变率等值线并给出了它们的分布图象。根据多年形变和现今GPS观测资料 ,分析福建地壳垂直运动与水平运动 ,显示区域应力场优势分布特征。最后 ,对福建沿海及台湾海峡地壳动力学特征作了初步的探讨     

A Forward Test of the Precursory Decelerating and Accelerating Seismicity Model for California

Accelerating strain energy released by the generation of intermediate magnitude preshocks in a broad (critical) region, and decelerating energy released in a narrower (seismogenic) region, is considered as a distinct premonitory pattern useful in research for intermediate-term earthquake prediction. Accelerating seismicity in the broad region is satisfactorily interpreted by the critical earthquake model and decelerating seismicity in the narrower region is attributed to stress relaxation due to pre-seismic sliding. To facilitate the identification of such patterns an algorithm has been developed on the basis of data concerning accelerating and decelerating preshock sequences of globally distributed already occurred strong mainshocks. This algorithm is applied in the present work to identify regions, which are currently in a state of accelerating seismic deformation and are associated with corresponding narrower regions, which are in a state of decelerating seismic deformation in California. It has been observed that a region which includes known faults in central California is in a state of decelerating seismic strain release, while the surrounding region (south and north California, etc.) is in a state of accelerating seismic strain release. This pattern corresponds to a big probably oncoming mainshock in central California. The epicenter, magnitude and origin time, as well as the corresponding model uncertainties of this probably ensuing big mainshock have been estimated, allowing a forward testing of the model's efficiency for intermediate-term earthquake prediction.     

基于地震带细致划分为地震构造区的概率地震危险性分析

在我国当前地震安全性评价中,普遍使用考虑地震活动时空不均匀性的概率地震危险性分析方法(CP-SHA),它规定以地震带为统计区域计算地震活动性参数b值和v4值。正在编制的中国地震动区划图(五代图)提出了针对潜在震源区进行三级划分原则:划分地震带、地震带上划分地震构造区、地震构造区内再划分潜在震源区。本文提出以地震构造区为统计区域回归统计方法计算b值和v4值,然后进行概率危险性分析计算,这样得到的结果可能更为合理。     

高频GNSS实时地震学与地震预警研究现状

为实现从注重灾后救助向注重灾前预防转变,如何提高地震灾害监测预警和风险防范能力成为我们关注的重点.本文给出了国际上GNSS位移记录、强震动加速度记录、测震速度记录在地震预警中的应用现状,并总结了各自的特点,归纳出围绕高频GNSS地震学在震级与破裂过程实时反演中的几个需要进一步研究的关键问题：（1）引入北斗系统,基于高频GNSS（GPS/BDS）双系统的实时位移解算方法来提高实时单站位移解算精度,使实时解算精度达到厘米级;（2）开展强震仪加速度记录基线偏移校正研究,弥补地震近场GNSS站密度不足问题;（3）强震仪加速度记录与GNSS位移记录特点不同,开展强震仪加速度数据与GNSS位移数据实时融合处理研究,快速获得包含丰富地震形变和速率的波形数据;（4）测震学方法可快速估算震级,但是在强震发生时会出现震级饱和现象,造成震级估算偏低.需要开展基于GNSS位移时间序列的多种方法相结合的实时震级估算方法研究,通过与地震学方法比较和结合,来得到精度高、计算快的震级估值算法;（5）基于高频GNSS、断层初始模型快速选取、断层尺度、参数自适应调整是快速判断断层破裂方向的基础,在断层破裂过程自适应准实时反演算法方面需要进一步加强.通过国内外研究现状调研、分析,表明基于高频GNSS地震学的震级快速确定、震源破裂过程准实时反演算法的发展将对我国地震预警系统从"二网融合"到"三网融合"提供坚实的技术支撑.     

THE MAXIMUM LIKELIHOOD ESTIMATION OF b-VALUE IN MAGNITUDE-FREQUENCY RELATION AND ANALYSIS OF ITS INFLUENCING FACTORS

b-value in the magnitude-frequency(G-R)relationship plays a vital role in seismicity research and seismic hazard analysis, and the most commonly used techniques to simulate it are least square approach and maximum likelihood method. Least square method is simple and easy to apply, therefore widely used in China. However, many researches show that there exist some limits in least square estimation of b-value. Earthquakes with different magnitudes are not equally weighted in this method, and larger events have higher weights, so b-value is vulnerable to the fluctuation of several big earthquakes; meanwhile, least square method needs to divide magnitude intervals artificially. With a small sample size, data points could be not enough if the magnitude interval is too wide, and events in a magnitude interval may be lacking if it is divided to be too narrow. Especially for incremental G-R relationship, it is possible that N(M_i)equals 0 in an interval with large magnitude, so log(N(M_i))loses meaning and has to be ignored, resulting in a low b-value. Therefore, under certain conditions, maximum likelihood method is recommended as an effective substitution or supplementary for least square estimation of b-value. Among numerous previous researches on maximum likelihood estimation of b-value, lots of equations have been provided, based on varied implicit assumptions and different ways of solution. A brief overview is first presented for these equations, and classification and summary are provided based on whether taking account of the effect of binned magnitude, with finite maximum magnitude, using unequal observation periods for different magnitude intervals, and with analytic solution or not. Following this, a total of 6 influential factors are analyzed, such as binning magnitude, measurement errors of magnitude, sample size, magnitude span, minimum completeness magnitude and fore- and aftershocks. At last, reasonable suggestions are provided for using those equations properly. The equations of Aki(1965), Utsu(1965), Page(1968)and Kijko and Smit(2012)are based on assumption that magnitudes are continuous random variables, and have no corrections for this, so these equations are not recommended here. For simplicity, the equations of Utsu(1966)or Tinti and Mulargia(1987)can be used, but magnitude span should be greater than 2.5 due to without finite maximum magnitude in the formulas. For researchers having capability to write code and calculate numerically, Weichert(1980)or Bender(1983)'s algorithm could be utilized. Especially when it is required to apply data with different observation periods for varied magnitudes, the formula of Weichert(1980)is recommended. This study contributes to more accurately understand and use different formulas of estimating b-value by maximum likelihood technique, which can be used as reference for peers.     

Spatial Distribution of Return Period in China Based on b-value Obtained by Using a New Method

ConsldeMle research has been po山med on the卜quency－删哪tude relation since仇tenbe恰Rlchterproposed the relationLog N＝a＋bM ()WhereM Is the nunlmumm昭nltude In the data sam卫ie for a celtaln region,N Is the cumulativeMquency ofevents In a celtaln time Interval ithm用刃ltude h锣rthan or equal to M,and a and bfilE COllstsfitS．Experimental researches on the G－R relation coMmed that the b decreases ith the de耶e ofbeteing6nelty and density Of the c。ohs In the。dlum(Mogi,1967)．The …     

汶川地震前后b值的时空变化及构造意义

地震b值与地应力具有密切的关系,即高应力对应低b值,研究汶川地震区b值的时空变化可以获得地震前后的应力演化过程.本文通过震源位置和速度模型联合反演的方法和双差定位法对龙门山地区汶川地震(M_s8.0)前后的地震进行了重定位工作,在此基础上精确计算了汶川地震前后不同时间段内b值大小的变化趋势及b值空间分布.结果显示:(1)b值在深度上存在明显的分层性,15 km以上b值较高,15 km以下b值较低,反映上部岩石为脆性,下部岩石为韧性的特征;(2)区域平均b值在汶川地震之后出现下降,而后有所恢复并逐渐增加,芦山地震后又降低,反映了龙门山断裂带的同震变化和震后愈合过程;(3)分时段b值空间分布结果显示,在芦山地震之后,汶川地震破裂区b值下降,即应力增加,这意味着汶川芦山地震之间的空段存在韧性变形,且芦山地震的发生促进了汶川地震破裂区的愈合;(4)映秀附近彭灌杂岩区在震前呈现低b值异常,反映该区域存在较大的差应力,可能是震后最大位错发生在该区域的根本原因.     

基于钻孔应变地震波记录确定地震面波应变震级

近年发展的高采样率的钻孔应变仪可以记录到地震波造成的水平应变,是动态库仑应力研究的重要手段.利用钻孔应变仪记录的应变地震波确定地震震级,是一个重要的科学课题.本文收集了我国10个应变台站四分量钻孔应变仪2020年1—3月记录的10 Hz采样应变地震波资料,共选出震级M≥4.0的浅源地震68个,用最小二乘法求得应变震级公式为Mε=lg E max+1.65 lg(Δ)+1.43.对于6.5级以下的地震,计算得到的应变震级ε与中国地震台网中心公布的震级M CENC基本一致:M CENC=1.03ε-0.23,但本研究中的两个7级地震,应变震级出现了震级饱和现象.本文的基本结论是:应变定震级是可行的,但对于两个7级地震的饱和现象需要进一步深入研究,因为它们不但涉及震级确定,而且涉及远震动态库仑应力触发地震的研究.今后要对更多的大震钻孔应变波形记录进行分析.