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.
-
引言
空间对地观测技术是20世纪科技发展的重大成果,以卫星为基础的遥感技术可实现远距离非接触性地球观测(杨哲海等,2003)。近几十年来,随着遥感平台、传感器及计算机技术的快速发展,遥感技术逐渐向微波和高光谱方向发展。高光谱遥感在20世纪80年代出现,其光谱分辨率达到纳米量级,形成的空间影像中的每个像素都具有其对应的光谱特征,提高了人类对客观世界的认识能力(童庆禧,2008)。
经过20多年的发展,利用高光谱遥感较高的光谱分辨率能获取较常规遥感手段更为详细的地面信息,受到了国内外对地观测研究领域及其他领域学者的关注,从数据获取到数据的处理与应用等方面都得到了快速的发展。高光谱遥感已在地球科学的各个领域,如精细农业、林业调查、大气环境监测、地质找矿及制图、识别伪军事目标和海洋生物等得到了广泛的应用,并逐渐起到了越来越重要的作用。在海洋遥感方面,利用高光谱遥感光谱响应范围广、分辨率高的优势,可以有效探测水色、水温等,同时可用于探测海洋水体表面油层、海水中叶绿素浓度、深海中物质的分布、海水中污染物的浓度等。通过研究海水中不同物质成分的光谱特征,了解深海层中物质的分布,更好地为海洋探测提供依据。在植被研究方面,由于不同植冠的组成成分不同,利用高光谱光谱范围窄、光谱连续的特点,可提取植被覆盖度、植被干物质信息及监测植被多样性等,还可以利用高光谱遥感监测植冠成分的变化,得到受病虫害或大气环境等影响的情况(宫鹏,1998)。在精细农业方面,高光谱遥感可以监测农作物长势,提取农作物含水量、叶面指数等参数,定量评估农作物产量,同时可以探测农作物土壤含水量、矿物质的含量等,为精细农业提供技术保障(刘伟东,2002)。在地质调查方面,高光谱具有很多不同于宽波段遥感的特点,利用波段窄的优势可以精确地探测到岩石中不同矿物质的存在,根据光谱曲线的不同,区分不同的矿物质组成成分(张成业等,2015)。高光谱遥感在蚀变带及蚀变矿物、油气探测、矿物识别与填图等方面也应用广泛。赵志芳等(2012)等利用ASTER影像数据对云南省保山市矿区矿化异常信息进行提取,利用USGS和AUS矿物光谱库,得到铁染蚀变异常信息,取得了比较好的效果。在军事侦察、伪装识别方面,由于目标经常隐藏在与自身色调一致的地物中,通过肉眼无法区分识别,依靠高光谱技术光谱分辨率高的优势,根据目标与伪装材料光谱特征的不同,能够在成像光谱仪获取的影像中快速准确地捕捉伪装目标。
高光谱遥感不同于传统的多光谱遥感技术,其在光谱分辨率上具有独特优势,可探测到地物目标的细节信息。归纳其特点及优势如下(杨国鹏,2007;路威,2005):
(1)波段多,波段连续,光谱分辨率高。在250—2500nm的全波段光谱范围内,高光谱遥感的光谱分辨率一般小于10nm,使得波段数目具有成百上千个,从而获得近乎连续的光谱曲线。
(2)数据量大,信息冗余增加。较常规的遥感影像,高光谱遥感影像具有高数十倍的波段数目,数据量较为庞大。影像波段越多,数据量越大,两者成指数关系;随着波段数目增多,波段之间的相关性变大,使得冗余信息量变大。
(3)数据描述模型多,分析灵活。高光谱遥感影像具有与常规遥感影像同样的图像模型和特征模型,同时因光谱分辨率较高,增加了光谱模型,使得地物识别能力得到了极大提高,能够精细探测到地表地物信息及覆盖材料等。
(4)反演定量化信息。高光谱影像涵盖的信息更加丰富,使得定量反演地表地物状态成为可能。
高光谱遥感具有独特的特征和优势,近年来在数据获取及特征分析技术方面获得了快速的发展,其影像数据也广泛应用于不同的行业,这项技术在地物精细分类、地质目标识别、地质灾害识别与趋势预测中占据着重要的位置。但在地震灾害识别领域,一方面受数据源获取的制约,另一方面受光谱特征分析与灾害地物识别技术的影响,针对高光谱遥感技术在震害目标识别方面的应用还未开展系统的研究。本文在分析高光谱技术在地震监测应用的基础上,通过构建简单的典型震害地物光谱库与分析地物目标的特征差异,探究其在震害目标识别中的应用,为后续开展高光谱震害识别系统研究奠定基础。
1. 高光谱遥感发展现状及其在地震监测中的应用
1.1 高光谱遥感发展现状
1983年,美国国家航空和航天管理局(NASA)喷气推进实验室设计研制了世界上第一台航空高光谱成像仪,获取的影像数据成功应用于地质填图、植被研究等领域,受到众多发达国家及研究者的重视。进入90年代后,各个国家均开始研制光谱成像仪,经过一段时间的发展,研制成功的一系列不同类型的高光谱成像系统在不同研究领域得到了广泛的应用,取得了丰硕的成果。到目前为止,高光谱遥感的相关研究已经产出了足够的可用于解决实际应用问题的技术方法,正由实验室研究阶段向实际应用过渡。
前期高光谱影像主要以航空影像为主。随后航天高光谱传感器研制成功,装载MODIS传感器的Terra卫星、美国NASA EOS计划中的Aqua卫星、装载GLI传感器的日本ADEOS-Ⅱ卫星分别于1999、2002年相继发射升空,标志着高光谱遥感技术已经步入与航空和航天相结合的阶段。2000年澳大利亚成功发射了一颗具有65个波段的光谱成像小卫星(ARIES),其波段涵盖了可见光、短波红外、近红外波段,卫星以地质勘测和矿藏资源作为主要遥感对象。2002年3月欧洲空间局(European Space Agency,ESA)发射了装载有10个卫星探测器的ENVISAT卫星,其中搭载的MERIS卫星探测器为高光谱成像仪,其获取的影像波段有15个。2000年携带高级陆地成像仪(Advanced Land Imager,ALI)和高光谱成像仪(Hyperion)的美国地球观测1号(Earth Observing-1,EO-1)系统以及携带海洋海岸成像光谱仪的HRST卫星的相继升空(EO-1 User Guide,2003),为高光谱技术的研究、实验提供了大量的数据。与此同时,我国密切关注国外高光谱技术的发展,并且结合我国的实际需求,在光谱成像仪研制方面也取得了较大的进步。针对某种特定的需求,相继研制成功了各自的专题应用光谱仪,如IR/UV海洋环境遥感监测专用扫描仪、VIS/MIR/IR森林火灾探测专用扫描仪、FIMS/ATIMS遥感地质矿物识别多光谱扫描仪等。2002年3月“神舟三号”飞船携带着中国科学院研制的中分辨率光谱仪升空,影像波段既包括可见光近红外波段,又包括短波红外、热红外波段,获取的影像空间分辨率为500m,使得中国航空成像光谱仪技术达到世界先进水平,从而更好地为地球全面观测和监测提供服务(胡秀清等,2009)。
在典型地物目标光谱特征数据库建设技术方面,美国最早开始进行研究,并建立了喷气推进实验室(Jet Propulsion Laboratory,JPL)波谱数据库和美国地质调查局(United States Geological Survey,USGS)波谱数据库,2个数据库免费开放,可供全球相关研究人员下载应用。除此之外,还有ASTER和IGCP-264 2个经典的波谱数据库,其涵盖了地表绝大部分地物的光谱曲线。除了通用的波谱数据库之外,同时也针对特定用途建立了专门的更为精细、更能满足专业应用的专业数据库,如英国为研究海水颜色的变化模型搭建的海水光谱数据库。国内经过近几年的发展,也取得了较大成果,尤其在光谱分析与数据库构建方面成果斐然。20世纪90年代初,中国科学院光机所、遥感所等单位通过野外采集典型地物光谱曲线、室内补充部分曲线的方法获取了一批光谱曲线,基于此,构建了我国第一个光谱数据库。此数据库涵盖的地物类型丰富,类别多样,是一个综合性的光谱数据库。1988年中国国土资源遥感航空物探与遥感中心针对行业应用,采集整理了地质行业中岩石矿物的光谱曲线,构建了典型矿物光谱数据库,从性质上来讲,此数据库为一个专业数据库。2000年中国科学院遥感与数字地球研究所收集整理了地物光谱并建立了典型地物数据库及管理系统,系统实现了波段曲线的动态生成和遥感器波段模拟功能。
自20世纪80年代开始至今,地物波谱仪、航空、航天高光谱成像仪得到了快速地发展,大幅度提高了高光谱影像光谱分辨率。美国Hyperion卫星波段数目为220,影像空间分辨率为30m(Kruse等,2011);德国、意大利、加拿大等国家计划发射ENMAP、HERO、PRISMA卫星空间分辨率均为30m(Kruse等,2011;Mielke等,2014, Bergeron等,2008;Rogge等,2014);我国2008年发射升空HJ-1A卫星波段数目为115个,空间分辨率为100m,2011年发射升空的“天宫一号”的空间分辨率为20m和30m,是我国到目前为止空间分辨率最高的成像光谱仪(周子勇,2014),即将发射升空的“高分五号”卫星空间分辨率也为30m。大型计算机及数据处理技术提高了遥感影像信息处理的效率,同时使得信息处理的算法不断更新完善,光谱分辨率的提升及典型地物光谱特征数据库的完善,使得利用高光谱探测原先传统遥感无法探测到的信息成为可能。
1.2 高光谱遥感在地震监测中的应用
随着高光谱传感器的陆续发射、遥感数据种类的增多以及影像图像处理技术与应用能力的不断提高,高光谱遥感与传统遥感相比具有更强的精细地物识别能力和定量反演能力,在微观方面表现更突出。目前,在地震相关的工作中基于高光谱遥感技术展开了一系列的研究,但是其应用并不成熟,仅有的一些应用是利用高光谱数据探测地震气体地球化学信息,以探讨地震活动并预测未来地震的空间位置(崔月菊等,2011)。地震的发生是由断层的相对运动造成的,断层的错动会引起近地表气体、湿度、温度等变化,高光谱遥感影像能捕捉到细微异常信息的变化,尤其是与地震相关的CO、水汽、O3、气溶胶与温度变化异常,通过对异常变化信息进行分析判断地震强度与空间分布位置,进而达到检测地震的目的。国内外的相关研究已证实近地表断层的异常与地震的发生紧密相关,如姚清林等(2005)分析缅甸北部MS6.9地震时期的高光谱影像发现,震前震中附近地表显示出CO含量异常;Singh等(2001a)发现2001印度Gujarat MS7.8地震的震前影像也存在CO含量异常;郭广猛等(2006)发现2002年台湾MS7.5地震影像也存在CO含量异常。除此之外,崔丽华(2009)对比分析2001年印度Gujarat MS7.8地震前后遥感影像发现,震后震中附近近地表位置存在水汽增高的异常现象,除此之外,在1985年墨西哥Michoacan MS8.1地震(Pulinets等,2007)和2008年汶川MS 8.0地震(Okada等,2004)中也发现了地震前后震中位置水汽含量出现了增高异常。此外,Ganguly(2011)分析2006年印度Gujarat MS5.5地震影像发现地震引起了O3浓度和大气气溶胶异常变化,这些异常变化在2002年11月的Denali地震(Singh等,2007)和2010年1月12日的海地地震(Singh等,2010b)中也被监测到。以上研究表明,高光谱遥感技术可以在一定程度上监测地表气体化学异常,从而反演地震发生的空间位置与大致时间,但引起地球化学异常的因素复杂多样,其成因机制也较复杂,遥感影像探测到的变化异常是否是由地震断层错位引起还未有明确的解释。
随着高光谱技术的发展,不同地物信息定量化提取及光谱数据自动化处理技术成为高光谱遥感研究领域的热点与难点。光谱特征曲线的丰富及自动识别的实现为高光谱应用领域的拓宽提供了良好的基础。
2. 高光谱遥感在震害识别中的应用
与地震化学监测技术相比,高光谱遥感在震害识别中的应用研究相对较少,国内外并未有系统的研究。由于地震破坏,地震地区下垫面特征与环境显得异常复杂。现阶段地震遥感调查过程中,采用的多光谱遥感影像具有较高的空间分辨率,能识别地物的空间位置与破坏情况,但是其光谱分辨率较低,在地物的精细过程中,普遍存在“同物异谱”和“异物同谱”的现象,从而影响地物提取精度。高光谱遥感影像将原有的波段范围进行细化分段,使得探测的地物光谱曲线更加精细,目标种类的数量也相应的增多,从而实现地物的精细分类。地震造成的地物破坏可引起影像相对应的光谱特征发生变化,不同的地物目标具有不同的光谱特征,通过识别其曲线的差别,将精细光谱特征与空间位置特征相结合,实现多特征信息的融合,以开展高光谱震害信息提取与震害评估研究,实现高光谱遥感影像震害目标的探测。
图 1(a)为几种不同建筑材料的光谱特征曲线。从图中可看出,建筑物混凝土与沥青混凝土的光谱曲线波形相似,但是在300—1200nm波长范围内,其反射率存在差异,沥青混凝土的反射率大于建筑物混凝土。建筑沥青及建筑黑漆指建筑物屋顶铺设沥青层或涂上黑漆,两者的光谱曲线在0—300nm波长范围内可与混凝土区分开。建筑物屋顶的沥青或黑漆可在600nm波长附近进行区分,黑漆的光谱曲线在附近存在一个波谷,而沥青不存在。
图 1(b)为表面光滑和粗糙红砖、松木、煤渣、瓷砖这几种常见建筑物材料的光谱曲线。可以看出,不同建筑材料的光谱曲线是不同的。即使相同的建筑材料,因表面光滑程度不同其光谱曲线也存在差异,如表面光滑与粗糙的红砖,在波长为50—270nm范围内,其反射率存在很大的差异,尤其在波长为80nm附近,表面粗糙的红砖存在一个波峰及波谷,而光滑的红砖反射率没有变化。
建筑物房屋、道路由不同的建筑材料构成,而不同建筑材料具有不同的光谱特征,收集整理这些光谱特性曲线并建立典型震害特征光谱库,利用震后获取的震区高光谱遥感影像与标准震害光谱库进行匹配,分析倒塌建筑物、损毁道路形成的碎屑与周边地物光谱特征的差异,便可快速精细地实现倒塌建筑物及损毁道路的提取,进而实现震害损毁程度遥感定量化评估。基于高光谱遥感影像的震害信息识别技术流程如图 2所示。
3. 典型震害地物目标识别
基于上述理论与流程,设计野外光谱数据采集实验。收集城市典型地物的光谱曲线,初步构建地物光谱特征库。地震发生时,震害地物主要表现为建筑物的破坏倒塌与道路的损坏。本文实验中选择的地物目标材料为建构筑常用的建筑材料,虽未处于震区,但在光谱采集过程中选取部分建构筑物破碎地带或边缘进行光谱收集,因此曲线一定程度上能反映震害地物目标的特征。实验采用美国ASD(Analytical Spectral Devices)公司的FieldSpec 3便携式地物波谱仪,采集时间选择光照强度良好的10时到14时,以满足光谱采集光照强度要求。实验共采集14种地物光谱曲线,每种地物至少采集3个样本,每个样本采集10条光谱曲线,采集的地物类型归纳为3类,分别为道路、建筑物屋顶和建筑物墙体。剔除坏的光谱曲线,统计计算获取的每种地物样本光谱曲线的均值并将其作为最终参考曲线。
图 3为实验采集的3种不同材质道路的光谱曲线。分析图 3,可以发现不同材质的道路在300—1100nm波长范围内反射率具有明显的差异,水泥道路的反射率最高,土质道路的反射率其次,沥青道路的反射率最低。通过获取的高光谱影像中反射率的差异便可区分提取不同材质道路信息。地震发生时,道路损毁情况主要分为2种,一种是地震的直接作用导致路基毁坏、路面断裂、塌陷、隆起,造成车辆无法通行;二是地震间接作用,如滑坡、泥石流、崩塌等地震诱发的次生灾害导致道路被掩埋或路面上存有崩落的大块滚石。地震直接作用导致的道路损毁因出露新鲜断面或露出底部泥土,导致与原路面的反射率不同,可直接提取损毁道路段;地震次生灾害导致的道路损毁情况下大量的堆积物在路面上,便可通过分析反射率获得堆积物所在的区域,从而达到提取损毁道路的目的。
图 4为实验采集的5种不同材质的屋顶和墙体的光谱曲线。由图 4分析,可发现墙体涂层的反射率最高,石棉瓦材质屋顶的反射率最低。在500—900nm范围内,铁皮屋顶反射率降低,在700nm时达到谷底,而大理石墙体在此范围内反射率升高,形成一个波峰,在750nm时达到波峰,依据此便可区分铁皮与大理石材质。
以“天宫一号”数据作为实验数据,因天宫一号卫星升空后,未拍摄到明显的震后场景,故以2013年广元市数据为例,因地面光谱曲线采集的地物目标与影像所处位置的地物存在差异,因此取多条地面光谱曲线的均值作为标准光谱曲线,影像与标准光谱曲线进行光谱匹配识别不同的地物目标。标准光谱曲线为模拟震区特征曲线,因此在光谱匹配过程中,设置光谱匹配一致度阈值为60%,即当匹配的2种光谱特征一致性优于60%时,影像中对应的像元即为光谱曲线对应的地物类型。由于地面实验中并未采集水体特征曲线,因此采用USGS波谱库中的水体特征曲线,采用光谱角分类的方法对地物进行分类识别,分类结果如图 5所示。
由图分析可见高光谱影像构筑物提取信息效果较好,建筑物、桥梁、道路都能清晰地表现出来。以目视解译结果作为参考基准,以提取的面积为单位计算提取精度,建筑物提取精度达98.65%。绝大部分桥梁信息能区别提取出来,可在图像上辨别区分。道路的提取效果不如建筑物和桥梁效果理想,但是国道、省道以及一些主要干道能完好地提取,只有个别的城市道路提取效果不理想。
根据上面的实验分析,可知不同材料具有不同的光谱特征,在遥感影像中将以不同的反射率表现出来,从而可与获取的高光谱影像对应。高光谱遥感影像波段间隔较高空间分辨率遥感短,因此能识别其无法识别的精细信息,尤其是震害屋顶的细微变化,同时基于不同地物光谱特征建立遥感地物目标指数,可实现震害信息的定量化分析,因此与传统的震害信息提取方法相比,利用高光谱遥感技术进行震害信息提取的类别更精细,更趋向于定量化分析。
4. 结论与展望
近几年随着传感器及传感平台的发展,高光谱数据的获取途径越来越广泛。地面ASD光谱仪或地面成像光谱仪可获取不同类型地物目标的光谱曲线特征,基于此可建立不同典型地物目标光谱数据库。地面光谱仪、无人机航拍都可获取图谱合一的图像,尤其是我国空间分辨率和光谱综合指标最高的空间光谱成像仪随“天宫一号”目标飞行器的成功升空后,在完成的地球环境监测等一系列实验中获取了大量的实验数据,在地物分类、波段范围等多方面已接近或达到国际同类产品水平。一系列高光谱数据的获取,为高光谱遥感在地震中的应用提供了便利的条件。获取的遥感影像可与地物光谱库进行对比分析,实现星地数据的交叉验证,进而提高信息识别精度与效率。
高光谱数据获取过程中存在大气、光照及其他条件的影响,震区现场地物复杂,影像存在混合像元的难题,影像光谱特征分析与匹配准确度在一定程度上限制了高光谱影像的应用。随着遥感图像处理技术的发展,基于光谱曲线剔除大气等条件影响的算法已相对比较成熟,从而保证了光谱特征曲线分析的精度。混合像元分解的关键在于计算混合像元中不同类型地物所占比例,现在已有成熟的算法,并在商业软件中可实现处理。现阶段,我国的遥感震害识别主要依赖高空间分辨率遥感,随着遥感传感器技术的发展,越来越多的传感器更加注重光谱分辨率与空间分辨率的结合,如何综合高光谱技术的地物精细分类与高空间分辨率的快速识别将是未来震害识别领域发展的趋势。本文重在论述高光谱技术在震害识别中的应用,对高空间分辨率影像与高光谱影像的结合未作深入探讨。
作为21世纪对地观测领域最重要的技术之一,高光谱技术在不同的领域得到了广泛的应用和快速的发展,但是在地震观测尤其是遥感震害评估领域应用较少,总结其原因主要有两点:① 高光谱数据较难获取,尤其是免费获取;② 缺乏普遍的研究方法及专业的光谱库。高光谱遥感技术在遥感影像上将不同地物的变化以反射率表现出来,形成不同的光谱曲线,从而获取地物微观信息。将高光谱技术引入到遥感震害评估领域能有效地提高震害评估类别精度,从而实现震害信息的定量化提取。
-
表 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 -