Seismic Hazard Assessment and Seismic Activity Analysis of Shanxi-Hebei-Inner Mongolia Area
-
摘要: 本文利用晋冀蒙交界地区1500年以来记录较完整的MS5.0及以上历史地震烈度资料和2000年以来的地震活动资料,分别进行地震灾害危险性和地震活动性分析。首先,利用ArcGIS将历史地震烈度资料数字化,没有等震线记录的地震用烈度衰减关系计算烈度圈半径。将研究区划分成0.1°×0.1°的网格,将烈度资料分配到与之相交的每个网格,并用烈度-频度关系计算每个网格的烈度a、b值。基于地震发生遵从泊松分布的假定,估算未来50年内晋冀蒙交界区遭受某一地震烈度的超越概率。同时,计算50年超越概率10%对应的地震烈度,计算结果表明张家口蔚县、阳原和山西广灵县周边地区的地震危险性较高。最后,采用中小地震能量密度值计算方法,对2000年以来的现今地震活动进行定量分析,与历史地震烈度资料分析结果进行对比,发现中小地震活动圈定的危险区与历史地震烈度资料评估计算的概率高值区相对应,因此,这两种方法可为晋冀蒙交界地区的地震活动性和危险性评价提供参考。Abstract: This paper analyzed the seismic hazard and activity in Shanxi-Hebei-Inner Mongolia area,based on the MS5.0 earthquake intensity data since 1500 and the seismic activity data since 2000. These historical seismic intensity data were digitized with ArcGIS,and the radius of the intensity circle without isoseismic records were calculated with the earthquake intensity attenuation relation. The research area was divided into a grid of 0.1°×0.1°,the intensity data were distributed to each grid that intersects it,and the intensity values a and b of each grid were calculated by using the intensity and frequency relationship. Based on the assumption that earthquake occurred in accordance with the Poisson distribution,the probability of surpassing a certain seismic intensity in Shanxi-Hebei-Inner Mongolia area in the next 50 years was estimated. At the same time,the earthquake intensity corresponding to the probability of exceeding 10% in 50 years was calculated. The results show that the seismic hazard is higher in the surrounding areas of Yu County,Yangyuan in Zhangjiakou and Guangling County in Shanxi. The energy density values of moderate and small earthquakes since 2000 were calculated and compared with the analysis results of historical seismic intensity data. The results show that the dangerous area delineated by the moderate and small earthquakes activities corresponds to the high probability areas calculated by the assessment of historical seismic intensity data. Therefore,these two methods can provide a reference for the seismic activity and hazard assessment in Shanxi-Hebei-Inner Mongolia area.
-
引言
地震小区划是对特定区域范围内可能遭遇到的地震影响进行划分,包括设计地震动参数小区划和地震地质灾害小区划。相比地震区划,地震小区划工作更重视局部场地条件的影响,更为详细地考虑周围的地震地质环境,其成果可为抗震设计、土地利用规划、震害防御等提供更精确的资料(胡聿贤,1999)。
地震小区划工作内容涉及地震活动性评价、地震构造评价、场地地震动分析等,工作量大,获取的数据也很丰富。传统地震小区划成果通常以文字报告、专题图件的形式提交,类型较单一,展示度较差,难以满足目前公共服务的需要。管理现有数据并充分利用数据推出更多服务型产品是目前包括地震小区划工作在内的地震工作信息化和现代化迫切需要。
随着GIS技术的发展,国内很多研究者逐步将GIS技术引入地震小区划研究工作中(张苏平等,2003;王庆满等,2011;李程程等,2014),建立有关数据库和数据管理系统(田勤虎等,2011;龚磊等,2015)。但已建系统多以数据查询、展示为主,缺少分析功能和实际产出服务。本文以嘉兴科技城地震小区划工作为基础,将基础数据与地理信息系统(GIS)相结合,形成集查询、分析于一体的地震小区划成果数字化系统,可为不同用户提供不同数据服务。
1. 研究区概况
研究区位于嘉兴科技城,行政区划属于嘉兴市南湖区,面积约30km2,是浙江省四大科技平台之一。现聚集了包括浙江清华长三角研究院、浙江中科院应用技术研究院在内的一批高端院所,还建有多个国际技术合作平台、工程中心和产学研基地,是嘉兴市乃至浙江省重要的创业创新平台。
查询《中国地震动参数区划图》(GB 18306—2015)(中华人民共和国国家质量监督检验检疫总局等,2015)可知,嘉兴市大部分地区由原来的Ⅵ度设防提升为Ⅶ度设防,设防水准的提高将对区内建设工程抗震设防、国土利用规划、社会经济发展等提出更高的要求。嘉兴科技城现处于Ⅶ度设防区内,加之人才产业聚集、产值高,地震灾害风险大。在该区进行地震小区划工作,可为该区提供更精确的抗震设防和土地规划资料,有效降低地震灾害风险。
2. 基础数据
根据相关规范要求,地震小区划工作内容包括地震地质及地震活动性资料的收集、场地钻探、场地土层反应分析等,涉及地质、地震、地球物理、工程勘察等多方面的数据,根据系统建设的需要归纳如下:
(1) 地震地质及地震活动性数据
研究区内地震地质及地震活动性专题图包括区域地震分布图、场地地震分布图、历史地震等震线分布图、区域构造图、区域断层分布图、近场构造图等,均为Mapinfo格式的矢量文件。
(2) 场地地震工程地质条件勘察数据
区划范围内布设60个钻孔,其中标准钻孔1个(孔深232.8m),地震钻孔29个(深度不小于105m),勘察孔30个(深度小于30m),钻探总进尺4204m,浅层人工地震探测测线共21534m。系统建设收集的数据包括60幅钻孔柱状图、18条工程地质剖面、5幅浅层人工地震探测解译图,其中钻孔资料均以柱状图的形式提供,格式为AUTO CAD矢量文件,钻孔的其他信息以Excel表格的形式提供;工程地质剖面图和地震探测测线数据为PDF文档和栅格图片;勘察报告以PDF文档的形式提供。
场地内60个钻孔均进行了波速测试,钻孔代表性土层动三轴土样80件,收集的波速数据和动三轴数据以Excel表格的形式提供。
(3) 地震动分析与参数区划数据
根据土层地震动分析的需要,建立29个场地钻孔模型,数据格式为TXT文本。收集全球范围内天然地震记录1692条,人工合成地震波522条(每个钻孔3个超越概率各6条),以TXT文本的形式提供。
地震小区划成果图件主要为地震动参数区化图,以Mapinfo矢量文件的形式提供。另外,还包括小区划报告的文本,以PDF文档的形式提供。
3. 系统设计
3.1 面向对象实现
地震小区划涉及钻孔、场地、强震记录、地震灾害等概念,其对象层次清晰,易于以面向对象的方式实现,形成地震小区划数字化系统。根据需求分析,首先将地震小区划数字化涉及的问题抽象成类,并建立类之间的消息机制及类之间的关系,即进行地震小区划数字化概念设计。
根据概念模型抽象出类,相应的类关系如图 1所示,地震小区划数字化主要类有CSolid(工程地质三维实体类)、CDrillHoles(钻孔集合类)、CSoilLayers(土层集合类)、CGroundMotion(强震动观测记录类)、CEarthQuake(场地地震动类)等。CDrillHoles类由CDrillHole类组合形成,并通过CreateTIN方法实现CTins类,由CTins类实现CSolid类,利用CGroundMotion类中SelectEarthQuake方法筛选出符合地质条件的历史强震记录,导入CSolid类,进行地震动分析,计算场地地震动参数峰值加速度、特征周期,并绘制峰值加速度等值线图和特征周期等值线图,由此生成场地地震小区划成果。
3.2 数据库设计
系统采用Access小型数据库实现场地基础地理信息数据、地震动数据、地震震陷数据、地震液化数据等的建库,结合GIS技术实现空间数据与属性数据的整合、矢量数据与栅格数据的整合、信息数据与分析数据的整合。数据库共包含9张数据表,分类如表 1所示。
表 1 数据库汇总Table 1. Database summary序号 表名 功能说明 1 ZkInfo 钻孔基本信息 2 ZkLayer 钻孔土层信息 3 StLayer 标准层序土层信息 4 ZkLique 钻孔砂土液化信息 5 ZkClay 钻孔软土震陷信息 6 ZkWave 钻孔土层剪切波速信息 7 ZkBG 钻孔土层标贯信息 8 ZkMx 钻孔动力分析模型信息 9 ZkSeis 钻孔动参数区划信息表 数据库逻辑设计如图 2所示。
3.3 系统物理设计
系统采用MapX作为GIS开发组件,数据库采用Access,在VB环境下完成开发,可运行于多个版本的windows操作系统中。
4. 系统功能实现
系统功能分为资料查询、成果应用、专题研究三大模块,其中专题研究整合了资料查询、成果应用模块的所有功能,详细的功能设计如图 3所示,系统主界面(专题研究模块)如图 4所示,主要对资料查询、成果应用模块进行介绍。
4.1 资料查询模块
该模块主要对地震小区划工作的基础资料及成果数据进行可视化展示,内容包括地震活动性、地震构造、场地勘察、地震区划等。由于部分功能展示的形式类似,对其中的部分功能进行介绍。
(1) 地震活动性查询:可查询区域地震、场地地震、历史地震影响烈度及场地地震综合评价。按地震震级自动统计地震活动性情况,通过地点地图及数据表的方式查询当前地震发生日期、经度、纬度、震级、震源深度等相关信息,成果展示如图 5所示。
(2) 地震区划:对区内潜在震源区的划分情况及地震小区划成果进行查询,成果以可视化的形式展示,包括区划的范围、相关的地震动参数及标定的反应谱,成果展示如图 6所示。
(3) 钻孔快捷查询:可对区内所有勘察钻孔进行集中展示,包括钻孔柱状图、钻孔土层剪切波速、抗震类别判定、砂土液化和软土震陷情况判断等所有信息,成果展示如图 7所示。
4.2 成果应用
该模块基于现有基础资料,通过建立分析模型对数据进行数字化分析,生成相应的成果数据。该模块主要包括虚拟勘察、场地分析及国土规划3个子模块,其中虚拟勘察模块包括虚拟钻孔、工程地质剖面、地层等值线等,场地分析模块包括砂土液化、软土震陷分析及场地抗震类别判定,国土规划模块包括断层分布及地震小区划。
(1) 虚拟钻孔:采用delaunay三角网上的线性内插算法,根据坐标输入或地图点选获得示范区内任意位置钻孔(虚拟)的相关资料,包括钻孔柱状图、场地抗震类别判定、砂土液化和软土震陷分析、历史地震、地震动参数等。虚拟钻孔的柱状图以栅格图片或CAD矢量文件的形式导出,所有虚拟钻孔信息还可以Word文档的形式导出,内容包括场地土层信息、场地地震动参数、场地类别等,涵盖了地震工程大部分基础资料,可供相关单位参考使用,如图 8所示。
(2) 工程地质剖面:根据连续的坐标输入或在底图上绘制剖面线生成研究区内相应的工程地质剖面图(见图 9),可将剖面图以栅格图片或CAD矢量文件的形式导出。
(3) 地震地质灾害:包括场地内的砂土液化和软土震陷。根据部分钻孔试验参数,采用建规中的判别方法对场地内钻孔进行判定,给出其液化程度和震陷情况。
(4) 地震动分析:根据点选的坐标位置,选择距离最近的地震钻孔数据,输入筛选要素,从天然地震时程库中选取符合要求的天然地震时程;根据人工合成的基岩地震动计算得到各超越概率下的场地地表地震动参数,成果展示如图 10所示。
5. 结论与建议
通过对地震小区划数据的整理,基于GIS软件开发地震小区划成果数字化系统。在实现项目成果数据查询与可视化展示的同时,还可导出数字化分析成果。相比传统的地震小区划工作,该系统的成果更丰富,形式更多样,为地震工作服务能力的提升进行了尝试。同时,需指出的是,此次建立的系统仅为单机版,与实现公共服务仍存在一定差距,后续还需进行更多的研究,补充和完善相关功能,以期建立高效便捷的在线服务系统。
-
表 1 晋冀蒙研究区内1500年以来MS≥5.0历史地震目录
Table 1. Historical earthquake catalogue with MS≥5.0 in Shanxi-Hebei-Inner Mongolia area since 1500
编号 发震时间/ (年-月-日) 震中位置 精度 震中烈度 震源深度 震级 震中参考位置 等震线描述 北纬/ (°) 东经/ (°) 1 1502-12-04 39.0 112.6 4 5$ \frac{1}{4}$ 山西代县西南 有等震线 2 1514-10-30 38.7 113.0 4 5$ \frac{1}{4}$ 山西代县南 有等震线 3 1545-02 40.0 114.0 4 5$ \frac{1}{4}$ 河北阳原一带 有等震线 4 1554-06-05 39.0 113.0 3 5 山西太原、大同间 平均有感半径约80km,太原、大同、盂县地震 5 1580-09-15 39.5 112.3 2 Ⅶ 5$\frac{3}{4} $ 山西平鲁 无等震线,井坪路(今平鲁)摇倒城墙 6 1581-05-28 39.8 114.5 2 Ⅶ 5$\frac{3}{4} $ 河北蔚县附近 有等震线 7 1582-03 40.1 113.2 2 Ⅵ 5 山西大同 无等震线,马邑(今朔县东北)亦震 8 1583-05-17 39.7 114.0 3 Ⅶ 5$ \frac{1}{2}$ 山西广灵、浑源一带 无等震线,怀仁、浑源、应县、山阴、朔县、马邑、广灵、灵丘及河北蔚县、涞源地震有感 9 1588-07 38.4 112.8 2 Ⅵ 5 山西忻县 无等震线,无记录 10 1615-12 40.0 115.0 3 5 河北宣化西南 有等震线 11 1618-11-16 39.8 114.5 2 Ⅷ 6$ \frac{1}{2}$ 河北蔚县附近 有等震线 12 1624-07-19 38.9 115.5 3 Ⅶ 5$ \frac{1}{2}$ 河北保定 无等震线,保定各州县地震,新城镇也震 13 1624年春 38.4 112.7 2 Ⅵ 5 山西忻县 无等震线,盂县地震 14 1626-06-28 39.4 114.2 2 Ⅸ 7 山西灵丘 有等震线 15 1628-10-07 40.6 114.2 2 Ⅷ 6$ \frac{1}{2}$ 河北怀安西洋河堡 有等震线 16 1657-10 40.2 115.0 3 5 河北宣化南 有等震线 17 1658-02-03 39.4 115.7 2 Ⅶ—Ⅷ 6 河北涞水 无等震线,定兴、新城镇、沧州、南皮、盐山、束鹿、山东乐陵地震 18 1664 38.7 112.7 3 Ⅶ 5$ \frac{1}{2}$ 山西忻县、代县间 无等震线,忻县、定襄、静乐、代县、五台、繁峙、崞阳压死人 19 1673-10-18 40.5 113.5 5 6$ \frac{1}{2}$ 山西天镇西北 有等震线 20 1678年夏 40.7 115.3 2 Ⅵ 5 河北宣化西北赵川 无等震线,无记录 21 1679-09-04 39.0 116.0 5$ \frac{3}{4}$ 河北雄县一带 有等震线 22 1683-11-22 38.7 112.7 2 Ⅸ 7 山西原平附近 有等震线 23 1720-07-12 40.4 115.5 2 Ⅸ 6$ \frac{3}{4}$ 河北沙城 有等震线 24 1724 40.5 115.3 2 Ⅵ 5 河北怀来新保安 无等震线 25 1765-07-04 40.1 116.0 3 5 北京昌平西南 有等震线 26 1898-09-22 39.1 113.0 2 Ⅶ 5$ \frac{1}{2}$ 山西代县 无等震线 27 1911-01-25 39.8 114.5 2 Ⅶ 5.9 河北蔚县 有等震线 28 1923-09-14 39.4 115.8 Ⅶ 5$ \frac{1}{2}$ 河北新城高碑店 无等震线,易县、安次有感 29 1952-10- 8 39.0 112.7 Ⅷ 5$ \frac{1}{2}$ 山西原平县崞阳镇东 有等震线 30 1957-01-01 40.5 115.5 2 Ⅵ 5 河北涿鹿附近 有等震线 31 1967-07-28 40.55 115.55 1 Ⅵ 10 5.4 河北怀来东北 有等震线 32 1976-04-06 40.20 112.10 Ⅶ 18 6.2 内蒙古和林格尔东南 有等震线 33 1981-08-13 40.58 113.41 1 Ⅶ 18 5.6 内蒙古丰镇东 有等震线 34 1989-10-19 39.92 113.91 1 Ⅷ 13 5.9 山西大同、阳高间 有等震线 35 1998-01-10 41.10 114.30 Ⅷ 10 6.2 河北张北 有等震线 表 2 地震烈度与地面峰值加速度的对应关系
Table 2. Relationship between intensity and peak ground acceleration
地面峰值加速度/g <0.05 0.05 0.10 0.15 0.20 0.30 ≥0.40 地震烈度 <Ⅵ Ⅵ Ⅶ Ⅶ Ⅷ Ⅷ Ⅸ -
陈阜超, 陈聚忠, 郑智江, 2015.晋冀蒙地区的垂直形变特征.大地测量与地球动力学, 35(3):453-456. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dkxbydz201503021 刁守中, 蒋海昆, 季同仁, 1988.场地烈度预测及其应用.中国地震, 4(4):27-33. http://www.cnki.com.cn/Article/CJFDTotal-ZGZD198804003.htm 丁宝荣, 杜轲, 孙景江等, 2018.考虑不确定性的地震烈度与地震动参数关系研究.地震工程与工程振动, 38(2):106-114. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzgcygczd201802012 高立新, 戴勇, 2012.晋冀蒙交界地区地震活动性特征分析.华北地震科学, 30(3):34-39. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hbdzkx201203007 黄玮琼, 李文香, 曹学锋, 1994.中国大陆地震资料完整性研究之一——以华北地区为例.地震学报, 16(3):273-280. http://www.cqvip.com/QK/93548X/19943/1310027.html 贾素娟, 鄢家全, 1986.场地影响烈度的频次特征及其工程地震的意义.国际地震动态, (7):3-8. http://www.cnki.com.cn/Article/CJFD1986-GJZT198607000.htm 蒋长胜, 吴忠良, 2005.由"重复地震"给出的中国地震台网的定位精度估计.中国地震, 21(2):147-154. doi: 10.3969/j.issn.1001-4683.2005.02.002 蒋海昆, 侯海峰, 王锜, 2000.华北地区大范围内中等地震活动平静的统计检验及其预测意义.内陆地震, 14(2):97-104. doi: 10.3969/j.issn.1001-8956.2000.02.001 李华玥, 2017.晋冀蒙交界区地震灾害损失预评估.中国地震局地质研究所. 刘静伟, 王振明, 谢富仁, 2010.京津唐地区地震灾害和危险性评估.地球物理学报, 53(2):318-325. doi: 10.3969/j.issn.0001-5733.2010.02.009 吕悦军, 2008.渤海海域地震区划研究.北京: 中国地质大学(北京). M7专项工作组, 2012.中国大陆大地震中-长期危险性研究.北京:地震出版社. 沈建文, 邱瑛, 张杰, 1988.场地烈度-频次关系与b值的联系及其在合肥市的应用.中国地震, 4(1):34-40. 王健, 2001.地震活动性图象处理的网格点密集值计算方法.地震学报, 23(3):262-267. doi: 10.3321/j.issn:0253-3782.2001.03.005 王晓青, 窦爱霞, 丁香等, 2015.地震烈度应急遥感评估研究与应用进展.地球信息科学学报, 17(12):1536-1544. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqxxkx201512017 谢卓娟, 吕悦军, 张力方, 2012.中小地震定量分析在地震区带划分中的应用——以龙门山地震带及邻区为例.西北地震学报, 34(3):277-283. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xbdzxb201203013 张力方, 吕悦军, 彭艳菊等, 2008.基于中小地震应变能密度的地震活动性图像分析.中国地震, 24(4):407-414. doi: 10.3969/j.issn.1001-4683.2008.04.011 中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会, 2015.GB 18306-2015中国地震动参数区划图.北京: 中国标准出版社. 朱红彬, 李红, 武敏捷等, 2014.晋冀蒙交界中强地震前的中等地震平静-活跃-平静演化, 29(3): 1114-1120. 朱亮, 缪鹏, 2017.基于应变能分析南北地震带中段地震活动对华东地区的影响.华南地震, 37(2):53-57. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hndz201702009 期刊类型引用(1)
1. 赵密,宋军,缪惠全,钟紫蓝,杜修力. 考虑地震强度空间变异性的地铁系统连通可靠性分析. 交通运输系统工程与信息. 2023(04): 203-210 . 百度学术
其他类型引用(2)
-