SEISMOLOGY AND GEOLOGY ›› 2023, Vol. 45 ›› Issue (2): 435-454.DOI: 10.3969/j.issn.0253-4967.2023.02.008

• Research paper • Previous Articles     Next Articles

PROBABILISTIC SEISMIC HAZARD ANALYSIS ALGORITHM INTEGRATING THREE-DIMENSIONAL FAULT SOURCES AND POTENTIAL SEISMIC SOURCE ZONE USING RANDOM SAMPLING

CHEN Kun1)(), GAO Meng-tan1), YU Yan-xiang1), XU Wei-jin1), DU Yi2), LI Xue-jin1), LU Dong-hua1)   

  1. 1)Institute of Geophysics, China Earthquake Administration, Beijing 100081, China
    2)National Institute of Natural Hazards, Beijing 100085, China
  • Revised:2022-12-09 Online:2023-04-20 Published:2023-05-18

融合三维断层源和二维潜在震源区的随机抽样概率地震危险性分析算法研发

陈鲲1)(), 高孟潭1), 俞言祥1), 徐伟进1), 杜义2), 李雪靖1), 陆东华1)   

  1. 1)中国地震局地球物理研究所, 北京 100081
    2)应急管理部国家自然灾害防治研究院, 北京 100085
  • 作者简介:陈鲲, 男, 1976年生, 2013年于中国地震局地球物理研究所获固体地球物理博士学位, 研究员, 主要从事地震区划、地震危险性分析及地震动强度图等方面的研究, E-mail: Chenkun-6620@163.com
  • 基金资助:
    国家重点研发计划项目(2022YFC3003505);中国地震局地球物理研究所基本科研业务专项(DQJB21X24)

Abstract:

Using the Monte Carlo random sampling method, a set of probabilistic seismic hazard analysis calculation programs that integrates our country’s traditional planar potential seismic source zone and three-dimensional fault sources is developed. The program is not only suitable for our country’s traditional regional area sources, but also considers the rupture scale of earthquakes and is compatible with the probabilistic seismic hazard calculation of three-dimensional fault sources. The algorithm developed in this paper efficiently realizes the three-dimensional simulation of the seismic event set of the fault source and introduces the earthquake rupture scale into the probabilistic seismic hazard analysis calculation in China, which significantly improves the rationality of the seismic hazard calculation in the near-fault area. In order to improve the execution efficiency of the program, the algorithm adopts the method of filling grid points in the planar potential seismic source zone in advance and randomly simulating the uniform distribution of seismic events in the planar potential seismic source zone. For the seismic hazard calculation of elliptical attenuation relationship, the algorithm uses pre-constructed three-dimensional matrices of the distance of the ellipse minor axis under different magnitudes, distances, and different angles between sites and the ellipse long axis direction of potential seismic source zone, and directly obtains the corresponding distance of ellipse minor axis through table look-up and interpolation. The algorithm developed in this paper avoids the problem of low computational efficiency in the iterative approximation of the distance of the ellipse minor axis. The mathematical expression of the three-dimensional fault source is based on the Frankel fault plane form of the 2002 edition of the National Seismic Hazard Map of the United States. The surface track and average dip Angle of the fault are used to create the rectangular fault plane, in which the dip direction of each rectangle is always perpendicular to the strike of its local fault segment. To maintain the coordination between the rupture area and the magnitude, the rupture of the earthquake occurring on the fault plane should not exceed the fault plane or the combination of fault planes. If the boundary of the rupture plane is outside the fault boundary, the entire rupture plane will move so that the boundary of the entire rupture plane matches the boundary of the fault plane. Using the probabilistic seismic hazard program of the Seismic ground motion parameters zonation map of China(2015)and the algorithm developed in this paper, the regional seismic hazard of the study area including Changsha-Zhuzhou-Xiangtan of the urban agglomeration in Hunan Province with moderate to strong seismic activity are calculated. Seismic hazard at different probability levels(return periods of 50.8, 475 and 2 475 years, respectively)for the Changde near-fault sources and Zhuzhou sites are also computed. The comparative study shows that the procedure of the Seismic ground motion parameters zonation map of China(2015)underestimates the seismic hazard near the three-dimensional fault source, and the degree of underestimation becomes more significant as the probability level decreases. Considering the influence of the earthquake rupture scale at the low exceedance probability level, the decomposition results of the seismic hazard for sites near fault show that the contribution of the seismic hazard is different from that of the traditional method of the Seismic ground motion parameters zonation map of China(2015), which mainly focuses on the earthquake of high magnitude. However, earthquakes of all magnitudes on the fault source can contribute to the seismic hazard, but the proportion of high magnitudes is the largest. Finally, an example verifying the probabilistic seismic hazard program(data set 1 case 10)from the Pacific Earthquake Engineering Research Center(PEER)is used to verify the reliability of the algorithm developed in this paper.

Key words: Probabilistic Seismic Hazard Analysis, potential seismic source zone, three-dimensional fault sources, Monte Carlo

摘要:

文中采用蒙特卡罗随机抽样方法, 研发了一套融合传统二维潜在震源区(下文简称“潜源”)和三维断层源的概率地震危险性算法。该算法不仅适用于传统的区域面源, 同时还能考虑地震的破裂尺度并兼容三维断层源的概率地震危险性计算。文中研发的算法可高效实现断层源地震事件集的三维模拟, 并将地震破裂尺度引入到概率地震危险性计算中, 显著提高了近断层地区地震危险性计算的合理性。为了提高程序的执行效率, 算法采用预先在平面潜源中充填网格点的方式随机模拟地震事件在潜源内的均匀分布。对于椭圆衰减的地震危险性计算, 算法采用了预先构建不同震级、距离及不同场点与潜源长轴方向夹角下的短轴距的三维矩阵, 通过查表和插值方式直接获得相应场点的短轴距, 避免了循环迭代逼近短轴距计算效率低下的问题。分别利用五代图的概率地震危险性程序和文中研发的算法, 计算了湖南长-株-潭(长沙-株洲-湘潭)城市群所处的中强地震活动环境的区域地震危险性以及近断层源的常德、株洲2个场点在不同概率水平下(重现期分别为50a、475a和2 475a)的地震危险性。比较研究表明, 五代图的程序低估了三维断层源附近的地震危险性, 且随着概率水平的降低, 低估的程度越来越高。最后, 利用太平洋地震工程中心(Pacific Earthquake Engineering Research Center, PEER)验证概率地震危险性程序的算例(数据集1案例10)验证了文中算法的可靠性。

关键词: 概率地震危险性分析, 潜在震源区, 三维断层源, 蒙特卡罗

CLC Number: