Ground Motion Simulation of Three-dimensional Sedimentary Basin Based on Directly-beneath Fault Dynamic Source Model
-
摘要: 沉积盆地与近断层地震共同作用会增加地震破坏的风险水平，尤其是盆地下方直下型断层发震情况。采用动力学震源模型刻画断层破裂发震过程，开展沉积盆地直下型断层谱元法地震动模拟研究，探讨不同断层面初始剪应力和成核区位置下三维沉积盆地地表响应规律。研究结果表明，断层面应力降对盆地地表地震动的影响显著，在断层面强度一定的情况下，随着初始剪应力的增大，即应力降增大，盆地地表峰值响应增大，原因在于应力降的改变影响了断层破裂释放能量，进而引起断层破裂速度改变，最终导致盆地地表响应发生变化；改变断层面成核区位置会对盆地内部地震动分布规律产生影响，当成核区位置从断层中间向断层左侧移动时，盆地左侧地震动逐渐减小，而右侧地震动逐渐增大，最终表现为盆地右侧地震动显著高于盆地左侧，原因在于改变成核区位置后，导致近断层地震动的方向性效应发生变化。Abstract: The combined action of sedimentary basin and near-fault earthquake will greatly increase the risk level of seismic damage, especially directly-beneath fault earthquake occurs under the basin. In this paper, the dynamic source model is used to depict the process of fault rupture and earthquake generation. This paper studies on the seismic simulation of sedimentary basin with directly-beneath fault by using the spectral element method, and discusses the surface response of the three-dimensional sedimentary basin under different initial shear stresses of fault and locations of nucleation area. The results show that the stress drop of the fault has a significant influence on ground motion in the sedimentary basin, when the initial shear stress increases, the stress drop increases, and the rupture velocity of fault increases obviously, which lead to amplification of the surface response inside the basin. Meanwhile, the locations of nucleation area also has an impact on its ground motion response, when the nucleation area moves from the middle of the fault to the left, the ground motion on the left side of the basin decreases gradually, and the ground motion on the right side increases gradually. The reason is that changing the location of the nucleation area leads to the change of the directivity effect of near-fault ground motion.1) 2
https://github.com/geodynamics/specfem3d.git -
引言
地震灾害风险评估是针对区域内遭受不同强度地震后造成的人员伤亡和财产损失等可能后果的定性、定量分析和评价工作（黄崇福，1999；姚清林等，2002；刘吉夫等，2008），开展地区地震灾害风险评估工作是践行新时期防震减灾工作的必然要求，是由震后应急处置被动模式向灾前常态准备模式的有效尝试（唐丽华等，2016），也是震前充分了解抗震薄弱环节和高风险地区，有针对性地采取防范加固和减灾措施，切实提高地震灾害风险防御能力，有效降低地震灾害损失的有效途径（唐丽华等，2013；Kircher等，2006）。
目前国内外已有众多专家学者、科研机构开展了地震灾害风险评估分析研究，国外现有灾害风险评估系统相对成熟，如美国HAZUS系统（Kircher等，2006）、联合国RADIUS系统、挪威SELENA系统（Molina等，2006）等。国内地震灾害风险评估的研究起步较晚，王晓青等（2000）以地震危险性分析、建筑物易损性分析、数理统计方法等为依据，研发了地震灾害损失预测系统；高杰等（2005）通过归纳总结防震减灾系统的特点，利用VB语言研发了基于GIS平台的地震灾害风险评估系统；陈焜浩等（2017）基于模块化思想，利用MEF插件式开发框架，构建了城镇地震灾害风险评估系统。上述地震灾害风险评估系统主要应用对象为国家、省、市等，产出的评估结果多以区域或城市为单位，提供的地震灾害防御对策和应急辅助决策也较宏观，对市县层级开展地震灾害风险评估适用性和针对性不强，且基本未考虑由地震引发的次生地质灾害的影响。
市县作为依法履行政府防震减灾社会管理职能的行政主体，是防震减灾工作面向社会最直接、最有效的力量，是发挥政府职能、强化社会管理和公共服务的基础。基于现有地震灾害风险评估系统存在市县层级应用薄弱、针对性不强、未考虑地质灾害影响等问题，本文在总结最新研究成果的基础上，以宝兴县为例，设计并构建操作方便简单、针对性强的地震灾害风险评估系统。为地方政府地震灾害风险管理、防御和风险处置等提供方法和方案，切实提高地震灾害防御能力，有效降低地震灾害损失。
1. 地震灾害风险评估基础数据及模型
1.1 基础数据
地震灾害风险评估系统涉及的基础数据主要包括基础地理数据、房屋建筑数据、统计数据、地质灾害隐患点数据和地震地质构造数据等。
基础地理数据包括行政区划数据、地形地貌数据、道路数据、河流水系、重点目标数据等。房屋建筑数据是以天地图中的单体建筑物数据为基础，搭建GIS服务器和Web服务器，开发基于微信平台的数据采集软件，实现现场定位采集单体建筑物的属性信息，并完成数据入库。统计型数据主要指人口统计数据和经济统计数据，源自2017年宝兴县统计年鉴。其中人口统计数据包括总人口、总户数、户籍人口、不同年龄段的人口；经济统计数据指地区生产总值、三大产业生产总值和社会固定投资总额。地质灾害隐患点数据包括滑坡、泥石流、崩塌等地质灾害隐患点。地震地质构造数据主要包括活动断裂数据、历史地震数据、场地类别数据、地震区划和地震带数据等。部分数据库表结构如图 1所示。
1.2 模型
地震灾害风险评估模型是进行地震灾害风险评估的基础，主要包括烈度衰减关系模型、地质灾害稳定性评估、火灾分析、建筑物破坏预评估、人员伤亡预评估、直接经济损失预评估、危旧房屋评估、高风险小区评估等系列模型。以地质灾害稳定性评估、高风险小区评估等主要模型为例，介绍如下：
（1）地质灾害稳定性评估模型
根据烈度衰减关系模型生成的影响场，结合地质灾害隐患点所处的地震地质构造、地形地貌条件，选取坡高、岩性构造、坡角、地震、降雨作为影响地震触发地质灾害的主要因子，采用综合指标法建立地质灾害稳定性模型，并与地质灾害发生的临界值进行对比，判断地质灾害发生的危险程度。本文采用丁彦慧等（2000）提出的地震崩滑稳定性计算模型，具体公式如下：
$$ Z = \left({P \cdot H \cdot S} \right) \times \left({W + Q} \right) $$ (1) 式中，Z为边坡危害程度，P为坡角影响系数；H为坡高影响系数；S为岩性条件系数；W为降水强度系数；Q为地震强度系数。
坡高是影响地震触发崩滑泥石流的主要因子，一般坡高越大，越易发生崩滑。坡高因子是依据DEM生成的，利用DEM提取山谷线，搜索距离位置点最近的山谷线，计算该位置点与谷底之间的高差将坡高划分称为5类，分别进行赋值，坡高越高对应的分类值越大。
岩性构造主要控制地质灾害类型，不同类型岩体触发的地质灾害类型、概率和规模不同，一般硬岩主要发生大型崩滑，而软岩以浅层滑坡为主。通过分析新构造运动强弱、历史上发生6级以上地震的频率、主干断裂带分布、岩体类型及结构等影响因素，将其划分为5类，并分别赋值。
坡角是控制地震触发崩滑泥石流具体发生位置的重要影响因子，一般地形坡度角越大越易触发，利用研究区DEM生成。该因子被分为5类，坡度越大，对应的分类值越大。
地震强度因子采用地震动参数，根据烈度与地震动加速度的对应关系进行转换，地震动加速度越大，对应的分类值越大。
年均降雨量是影响地震触发崩滑泥石流的外部因素，根据年降雨及降雨强度角将其划分为5类，分别进行赋值。降雨量＞800mm且多暴雨的分类值最大，降雨量 < 400mm且分布均匀的分类值最小。
（2）高风险小区评估模型
高风险小区是预测未来发生地震时，地震灾害及其引发的次生地质灾害、城镇潜在的各种隐患与不利因素等相互作用下，划定的危险性系数较高的小区。在建筑物震害预评估的基础上，综合考虑次生灾害等危险源的影响，采用综合指标法确定高风险小区，选取老旧房屋占有率、各类建筑物综合震害指数、距次生灾害源的距离、距断裂带的距离及距滑坡隐患点的距离，采用空间叠加的方式得到格网综合指数，采用分级分类方法将危险性等级划分为严重、中等和轻微3个等级，并以专题图的形式展示（图 2）。
2. 地震灾害风险评估系统设计
2.1 逻辑架构设计
地震灾害风险评估系统总体采用分层架构，包括运行环境层、支撑平台层、数据层和应用层，如图 3所示。运行环境层主要包括系统运行需要的硬件环境和网络环境。支撑平台层包括系统运行需要的一系列辅助软件系统，如数据库管理系统、GIS Runtime、应用中间件、Office和操作系统。数据层主要包括基础数据库、评估模板库、评估结果库和管理数据库。应用层是在前端直观展示及操作该软件系统的所有功能，包括地震危险性分析、次生灾害风险分析、地震灾害损失预评估、地震灾害防御对策、基础数据调查与管理和系统管理。
2.2 系统功能设计
地震灾害风险评估系统功能由地震危险性分析、次生灾害风险分析、地震灾害损失预评估、震害防御对策和系统管理构成（图 4）。
（1）地震危险性分析
地震危险性分析包括地震地质构造、地震活动性、危险性分析结果及地震事件管理等，即包括地震活动性、地质构造评价等方面的研究成果，以地图的形式将地震活动性、活动断裂分布、地质构造及不同超越概率的地震动参数等可视化展示。地震事件管理包括设定地震管理和设定烈度管理，在设定地震评估的基础上，系统增加了设定烈度评估的方式。设定烈度评估的方式可将评估区域在Ⅵ、Ⅶ、Ⅷ、Ⅸ、Ⅹ度下的地震灾害损失评估出来，作为分烈度损失计算图层，在破坏性地震发生后，可大大提高应急评估计算效率。
（2）次生灾害风险分析
次生灾害风险分析主要包括地质灾害风险分析、火灾风险分析和其他潜在危险源风险分析等。考虑宝兴县自然地理、社会经济等原因，全县褶皱密集，断裂发育，以高山地貌为主，地势起伏大，导致全县地质灾害较严重，且通过实地调研过程中发现宝兴县建筑物结构类型以木结构为主，若遇火灾后，易加剧灾情，快速蔓延，因此本文重点对地质灾害风险和火灾风险进行分析。基于地质灾害隐患点（滑坡、泥石流和崩塌）数据，并结合隐患点所处地形、地质条件和烈度，根据地质灾害稳定模型计算地质灾害点的稳定系数，从而判定其稳定程度，以图表的形式显示地质灾害稳定、失稳及可能失稳的比例，并根据失稳的地质灾害隐患点坡向和坡高计算其影响范围。火灾分析是通过综合考虑单体建筑物结构类型、用途、建筑物在地震后的破坏状态及天气、风速等环境条件因素，根据火灾着火概率模型，计算单体建筑物着火概率，以专题图和图表的形式综合展示。
（3）地震灾害损失预评估
地震灾害损失预评估由建筑物震害预评估、道路交通震害预评估、人员伤亡预评估和经济损失预评估等模块构成。建筑物震害预评估是在基于用户选择的设定地震或设定烈度模式下进行评估，在设定地震下根据地震烈度衰减关系生成影响场，结合建筑物易损矩阵及建筑物实际信息，对建筑物受损程度进行量化，并在此基础上考虑由地震引发的次生地质灾害的影响，综合评定建筑物损毁情况，以图表的形式进行展示。根据地震影响场范围及次生地质灾害范围，利用道路阻断模型对道路震害情况进行预测。在建筑物震害预评估的基础上，调用人员伤亡模型和经济损失模型对人员伤亡和经济损失进行预测，并以图表的形式进行展示。
（4）震害防御对策
震害防御对策包括建筑物加固措施、防灾规划及危旧房屋、高危险小区、次生灾害、地质灾害等地震灾害风险处置对策。该模块在次生灾害风险分析和地震灾害损失预评估的基础上，找出各类建筑物、生命线工程系统的抗震薄弱环节，通过设置建筑物年代与建筑物破坏程度等特征识别危旧房屋和高风险小区；根据地质灾害稳定性分析结果，在地图上绘制其大概影响范围，并给出失稳或可能失稳状态的地质灾害隐患点应对措施；根据区域内其他危险源在地震中的损毁程度，在地图上绘制其大概影响范围，并给出应对措施；综合地震地质灾害危险性和影响范围，划分出不利建设区、有利建设区和危险建设区；根据设定地震的预测受灾情况，给出储备仓库、储备物资及应急避难场所的需求建议，并结合辅助决策报告模板生成地震灾害防御辅助决策报告，为政府部门制定防震减灾规划提供依据。
（5）系统管理
系统管理由基础数据预处理、用户权限管理、数据采集移动设备管理和日志管理构成，其中基础数据预处理主要提供统计型数据的编辑更新操作及预先计算建筑物在各烈度下的损毁情况等。该系统基于角色访问控制权限，通过角色与权限管理用户，明确不同用户的权限。
2.3 业务流程设计
利用地震灾害风险评估系统进行灾害风险评估时，各功能模块之间的数据流按一定数据导向进行相互调用，基于系统的总体逻辑结构和功能结构，设计总体业务流程如图 5所示。
首先进行地震危险性分析，根据地震危险性分析结果，选择设定地震或设定烈度方式进行评估。在设定地震下根据地震烈度衰减关系生成影响场，结合滑坡、泥石流、崩塌等地质灾害隐患点数据，分析由地震引发的地质灾害危险及影响范围。在此基础上，利用地震灾害损失评估模型预测地震灾害损失，得到建筑物震害分布图、人员伤亡分布图、经济损失分布图、单体建筑震害分布图、道路震害分布图、危旧房屋分布图、高风险小区分布图、地质灾害稳定性分布图、次生灾害分布图等。根据地震灾害损失预评估结果，给出对应的防震措施，包括建筑物防震措施、危旧房屋集中地段防震措施、高风险小区防震措施、薄弱生命线工程防震措施、地质灾害抗震措施、次生灾害抗震措施、物资储备信息展示、应急场所建议及措施、防灾规划等，并将所有措施、建议整理成报告文档，生成地震灾害风险处置对策，为政府部门采取有针对性的防范加固措施提供信息支撑。
3. 地震灾害风险评估系统的实现
采用C/S模式对地震灾害风险评估系统进行开发，系统硬件环境包括客户端和服务器端，服务器端包括数据库服务器和应用服务器。对于系统软件环境，服务器端操作系统需安装Windows Server 2008及以上版本，数据库选择SQL Server 2008及以上版本，中间件采用ActiveMQ，客户端操作系统为Windows 10（32/64位）、ArcGIS Engine 10.0及以上版本。借助Eclipse、Visual Studio集成开发环境，采用Java语言定制开发。系统主要功能操作界面如下：
（1）地震危险性分析
地震区划查询界面如图 6所示，左侧界面展示地质分布图，通过在地图上随机选点，以对话框的形式展示该点位置信息，且能在右侧界面中部展示地震动峰值加速度及地震动加速度反应谱特征周期，并在界面右下方绘制该点在不同场地类型条件、超越概率下地震动参数对应的地震动峰值加速度及地震动加速度反应谱特征周期折线图。
（2）次生灾害风险分析
地质灾害分析结果界面如图 7所示，左侧界面控制地图图层的显示，中部以地图的形式展示地质灾害分析结果，依据地质灾害点稳定性判断结果以颜色区分，右侧界面提供按灾害类型、稳定程度条件进行查询的结果，以表格形式展示，并将表格内容进行统计分析生成文本，默认统计全部灾害类型和稳定程度。火灾分析结果界面如图 8所示，中部根据单体建筑物火灾发生概率对其进行分级显示，右侧界面以表格的形式展示单体建筑物着火概率。
（3）地震灾害损失预评估
单体建筑物地震灾害损失预评估结果界面如图 9所示，左侧界面展示单体建筑物震害结果，根据损坏情况进行符号化分级显示。右侧界面以柱状图和表格的形式展示以乡为单位的建筑物损坏情况。
图 9 单体建筑物地震灾害损失预评估结果Figure 9. Pre-assessment results of earthquake disaster losses of single buildings（4）震害防御对策
危旧房屋主界面如图 10所示，左侧界面地图显示根据不同建筑物年代、建筑物损毁情况等绘制的危旧房屋区域，右侧界面以表格的形式展示危旧房屋模型分析结果，且可单击表格中的某组数据在地图上进行定位，点击详情可弹出对话框，显示具体的防震措施，并以饼状图的形式统计分析不同灾害程度占比。高风险小区主界面如图 11所示，左侧界面地图显示高风险小区分析结果，右侧界面以表格的形式展示分析结果，且可单击表格中的某组数据在地图上进行定位，点击详情可弹出对话框，显示具体的防震措施，并以饼状图的形式统计分析不同灾害程度占比。
4. 结语
本文介绍了地震灾害风险评估系统所需的基础数据源和主要评估模型，阐述了系统的逻辑结构设计、系统功能设计和业务流程设计，利用Visual Studio开发平台，结合SQL Server数据库，采用组件式开发框架完成了系统研发。地震灾害风险评估系统丰富了市县层级防震减灾部门工作手段，为平时地震灾害防御提供辅助决策建议，提高了市县层级防震减灾能力。
本系统采用C/S模式，不同城市和地区需根据当地具体情况建设，从而造成资源浪费。同时对于市县而言，因缺乏专业的技术人员，会导致系统更新维护困难。为此，可引入云计算思想，在以后的研究中解决上述问题。
-
图 5 各观测点位错随时间变化曲线
Figure 5. The curves of displacement versus time corresponding to observation points
图 6 各观测点对应加速度时程曲线
Figure 6. The time-histories of acceleration corresponding to observation points
图 7 不同初始剪应力下地表P4观测点加速度时程曲线
Figure 7. The time-histories of acceleration at surface observation P4 with different initial shear stress
图 8 不同初始剪应力下地表观测点PGA及其放大系数曲线
Figure 8. The curves of PGA and amplification factor at surface observations with different initial shear stress
图 9 不同初始剪应力地表观测点P4加速度反应谱
Figure 9. The response spectrum of acceleration at surface observation P4 with different initial shear stress
图 10 不同初始剪应力断层破裂时间等值线图
Figure 10. The contour of fault rupture time with different initial shear stress
图 11 不同初始剪应力下地表速度波场快照图
Figure 11. The snapshots of velocity wavefield at surface with different initial shear stress
图 12 不同成核区位置断层破裂时间等值线图
Figure 12. The contour of fault rupture time with different location of nucleation area
图 13 不同成核区位置地表观测点PGA变化曲线
Figure 13. The curves of PGA at surface observations with different location of nucleation
图 14 不同成核区位置地表速度波场快照
Figure 14. Snapshots of velocity wave field at surface with different location of nucleation area
图 15 不同成核区位置地表水平峰值速度场及位移场分布
Figure 15. Distribution of peak ground velocity and displacement with different location of nucleation area
表 1 计算模型介质参数
Table 1. Parameters of the simulation model
介质密度/kg·m−3 剪切波速vS/m·s−1 压缩波速vP/m·s−1 厚度/km 深度/km 基岩 沉积盆地 基岩 沉积盆地 基岩 沉积盆地 3 000 2 400 3 200 300 5 500 500 0.3 0.3 3 000 2 400 3 200 500 5 500 800 0.3 0.6 3 000 2 400 3 200 800 5 500 1 500 0.4 1.0 3 000 — 3 200 — 5 500 — 29.0 30.0 
下载: 导出CSV
表 2 观测点位置
Table 2. Location of observation points
观测点编号 坐标（x，y，z） P1（盆地外） （10 km，0.5 km，0） P2（盆地右边缘） （5 km，0.5 km，0） P3（盆地右1/4处） （2.5 km，0.5 km，0） P4（盆地中心） （0 km，0.5 km，0） P5（盆地左1/4处） （−2.5 km，0.5 km，0） P6（盆地左边缘） （−5 km，0.5 km，0）
下载: 导出CSV
表 3 断层面上摩擦和应力参数
Table 3. Friction and stress parameters on the fault
断层参数 成核区外 成核区内 静摩擦系数 0.675 0.675 动摩擦系数 0.475 0.475 临界滑动距离/m 0.3 0.3 初始剪应力/MPa 66.0 81.6 初始正应力/MPa 120 120
下载: 导出CSV
表 4 不同初始剪应力断层平均破裂速度
Table 4. Average rupture velocity of fault with different initial shear stress
初始剪应力/MPa 初始正应力/MPa S值 平均破裂速度vR/m·s−1 剪切波速vS/m·s−1 vR/vS 64 120 2.43 2 620 3 200 0.82 65 120 2.00 2 761 3 200 0.86 66 120 1.67 2 877 3 200 0.90 67 120 1.40 2 966 3 200 0.93 70 120 0.85 4 011 3 200 1.25 71 120 0.71 4 480 3 200 1.40
下载: 导出CSV
-
巴振宁, 赵靖轩, 吴孟桃等, 2020. 逆断层地震作用下三维沉积盆地地震动谱元法模拟. 地震科学进展, 50（10）: 13—22 doi: 10.3969/j.issn.2096-7780.2020.10.002Ba Z. N. , Zhao J. X. , Wu M. T. , et al. , 2020. Spectral element method simulation of three-dimensional sedimentary basin under reverse-fault earthquakes. Progress in Earthquake Sciences, 50(10): 13—22. (in Chinese) doi: 10.3969/j.issn.2096-7780.2020.10.002 杜晨晓, 谢富仁, 史保平, 2008. 隐伏和出露地表断层近断层地表运动特征的研究进展. 震灾防御技术, 3（2）: 172—181 doi: 10.3969/j.issn.1673-5722.2008.02.008Du C. X. , Xie F. R. , Shi B. P. , 2008. The recent progress on strong ground motion near buried and surface earthquake faults. Technology for Earthquake Disaster Prevention, 3(2): 172—181. (in Chinese) doi: 10.3969/j.issn.1673-5722.2008.02.008 韩天成, 于彦彦, 丁海平, 2020. 直下型断层的破裂速度对盆地地震效应的影响. 地震学报, 42（4）: 457—470 doi: 10.11939/jass.20190177Han T. C. , Yu Y. Y. , Ding H. P. , 2020. Influence of rupture velocity of the directly-beneath fault on the basin seismic effect. Acta Seismologica Sinica, 42(4): 457—470. (in Chinese) doi: 10.11939/jass.20190177 胡进军, 2009. 近断层地震动方向性效应及超剪切破裂研究. 哈尔滨: 中国地震局工程力学研究所.Hu J. J., 2009. Directivity effect of near-fault ground motion and super-shear rupture. Harbin: Institute of Engineering Mechanics, China Earthquake Administration. (in Chinese) 刘启方, 金星, 丁海平, 2008. 复杂场地条件下震源参数对断层附近长周期地震动的影响. 地球物理学报, 51（1）: 186—196 doi: 10.3321/j.issn:0001-5733.2008.01.023Liu Q. F. , Jin X. , Ding H. P. , 2008. Effects of the source parameters on long period near-fault ground motion in the case of complex site condition. Chinese Journal of Geophysics, 51(1): 186—196. (in Chinese) doi: 10.3321/j.issn:0001-5733.2008.01.023 刘启方, 李雪强, 2011. 唐山大地震近场宽频带地震动模拟. 地震工程与工程振动, 31（5）: 1—7 doi: 10.13197/j.eeev.2011.05.001Liu Q. F. , Li X. Q. , 2011. Broad-band strong motion simulation of the great Tangshan earthquake. Journal of Earthquake Engineering and Engineering Vibration, 31(5): 1—7. (in Chinese) doi: 10.13197/j.eeev.2011.05.001 刘中宪, 刘明珍, 韩建斌, 2017. 近断层沉积盆地强地震动谱元模拟. 世界地震工程, 33（4）: 76—86Liu Z. X. , Liu M. Z. , Han J. B. , 2017. Spectral-element simulation of strong ground motion in the near-fault alluvial basin. World Earthquake Engineering, 33(4): 76—86. (in Chinese) 吴永祺, 张海明, 2022. 断层几何形态和自发破裂的动力学参数对地震波场的影响. 地球物理学报, 65（3）: 965—977 doi: 10.6038/cjg2022P0251Wu Y. Q. , Zhang H. M. , 2022. Effects of fault geometry and dynamic parameters of spontaneous rupture on seismic wave field. Chinese Journal of Geophysics, 65(3): 965—977. (in Chinese) doi: 10.6038/cjg2022P0251 于海英, 王栋, 杨永强等, 2008. 汶川8.0级地震强震动特征初步分析. 震灾防御技术, 3（4）: 321—336 doi: 10.3969/j.issn.1673-5722.2008.04.001Yu H. Y. , Wang D. , Yang Y. Q. , et al. , 2008. The preliminary analysis of strong ground motion characteristics from the MS8.0 Wenchuan earthquake, China. Technology for Earthquake Disaster Prevention, 3(4): 321—336. (in Chinese) doi: 10.3969/j.issn.1673-5722.2008.04.001 余厚云, 2020. 川滇地区破坏性地震的震源动力学过程及强地面运动模拟与震害评估. 合肥: 中国科学技术大学.Yu H. Y., 2020. Rupture dynamics and strong ground motion simulation of destructive earthquakes in the Sichuan-Yunnan region. Hefei: University of Science and Technology of China. (in Chinese) 张宇翔, 袁志祥, 2010. 汶川8.0级地震陕西灾区震害特征分析. 地震研究, 33（3）: 329—335 doi: 10.3969/j.issn.1000-0666.2010.03.014Zhang Y. X. , Yuan Z. X. , 2010. Analysis on the characteristics of earthquake damages of the MS8.0 Wenchuan earthquake in the affected area in Shaanxi. Journal of Seismological Research, 33(3): 329—335. (in Chinese) doi: 10.3969/j.issn.1000-0666.2010.03.014 朱刚, 高孟潭, 任雪梅, 2007. 中国城市直下型地震危险性分析研究. 中国地震, 23（2）: 124—130 doi: 10.3969/j.issn.1001-4683.2007.02.002Zhu G. , Gao M. T. , Ren X. M. , 2007. Analysis of seismic hazard in the cities of China. Earthquake Research in China, 23(2): 124—130. (in Chinese) doi: 10.3969/j.issn.1001-4683.2007.02.002 Anderson J. G. , Bodin P. , Brune J. N. , et al. , 1986. Strong ground motion from the Michoacan, Mexico, earthquake. Science, 233(4768): 1043—1049. doi: 10.1126/science.233.4768.1043 Andrews D. J. , 1976 a. Rupture propagation with finite stress in antiplane strain. Journal of Geophysical Research, 81(20): 3575—3582. doi: 10.1029/JB081i020p03575 Andrews D. J. , 1976 b. Rupture velocity of plane strain shear cracks. Journal of Geophysical Research, 81(32): 5679—5687. doi: 10.1029/JB081i032p05679 Ida Y. , 1972. Cohesive force across the tip of a longitudinal-shear crack and Griffith’s specific surface energy. Journal of Geophysical Research, 77(20): 3796—3805. doi: 10.1029/JB077i020p03796 Pitarka A. , Irikura K. , Iwata T. , et al. , 1998. Three-dimensional simulation of the near-fault ground motion for the 1995 Hyogo-Ken Nanbu (Kobe), Japan, earthquake. Bulletin of the Seismological Society of America, 88(2): 428—440. doi: 10.1785/BSSA0880020428 -
下载:
点击查看大图
计量
- 文章访问数: 216
- HTML全文浏览量: 161
- PDF下载量: 42
- 被引次数: 0
