Preliminary Study on Fault Parameters Based on National Seismic Data——An Example of Tangshan Earthquake
-
摘要: 随着地震观测台网的加密,中国地震台网统一地震目录越来越完整。本文以地震资料丰富且研究程度较高的唐山地震为例,尝试利用中国地震台网统一地震目录直接确定断层面参数。将小震确定断层面参数的方法应用于唐山地震序列,对其断层面参数进行分段拟合,其中唐山断裂南段走向和倾角分别为213.4°和81.9°,唐山断裂北段走向和倾角分别为231.4°和89.1°,滦县段走向和倾角分别为125.1°和76.2°,卢龙断裂走向和倾角分别为46.1°和89.3°,宁河断裂走向和倾角分别为246.6°和81.8°。将所得结果与前人研究成果以及现今震源机制解参数研究的误差水平进行对比分析,证明了本研究具有较高的可靠性。本文研究结果表明,在地震活跃且地震定位精度较高的地区,中国地震台网统一地震目录可以广泛应用于断层面参数的确定。
-
关键词:
- 中国地震台网统一地震目录 /
- 断层面参数 /
- 拟合 /
- 唐山地震
Abstract: With increasing density of seismic observation network, the seismic catalogue of China Earthquake Networks provided by the China Earthquake Data Center is getting more and more complete. In order to find out if the seismic catalogue of China Earthquake Networks can be used directly to determine the parameters of earthquake fault plane, we take the Tangshan earthquake with abundant seismic data as an example to determine the parameters of earthquake fault plane directly by using the seismic catalogue of China Earthquake Networks. The method of determining earthquake fault plane parameters based on the data of small earthquakes is applied to Tangshan earthquake sequence in this study. The parameters of earthquake fault plane are then fitted in different segments, the strike and dip of north-segment Tangshan fault is 213.4° and 81.9°, and the strike and dip of south-segment Tangshan fault is 231.4° and 89.1°. The strike and dip of Luanxian-segment is 125.1° and 76.2°, the strike and dip of Lulong fault is 46.1° and 89.3°, and the strike and dip of Ninghe fault is 246.6° and 81.8°. The comparative analysis of the published results of Tangshan earthquake and the error level of the present focal mechanism solution parameters suggest that our results are acceptable, indicating that the seismic catalogue of China Earthquake Networks can be directly used to determine the parameters of earthquake fault plane in areas with abundant seismic data and high-precision seismic location. -
引言
在建筑抗震设计和地震灾害评估等工作中,均需考虑建筑结构类型问题,并对结构上的差异进行分别处理。地震风险分析和地震保险工作也涉及建筑分类问题,地震风险研究与应用早期体现为地震灾害损失评估与预测工作。美国开展地震损失评估较早,国际知名的组织机构和风险模型公司已设计出针对地震巨灾的分析模型,如美国联邦应急部的HAZUS模型(FEMA,2018)及著名美国企业RMS和AIR的相关模型,均对所研究的建筑进行了详细的分类。1982年10月FEMA与应用技术委员会(ATC)签订合同,对美国加州地区地震风险开展评估工作,对地震破坏损失进行估计,编制地震破坏损失清单,以确定各工程设施在不同MMI烈度水平地震作用下的易损性矩阵,并发表了《加利福尼亚未来地震的损失估计》(ATC-13)(美国应用技术委员会,1991)。ATC方法提出对于预测对象需从工程结构和使用功能2个角度进行分类研究。即使2栋设施的建筑结构相似,但用途不同,所造成的经济损失和社会影响可能差异很大。ATC-13明确地将建筑设施的使用功能和破坏程度联系起来,根据工程设施的尺度、结构系统和类型,从工程结构角度,根据在加利福尼亚现有的结构清单中所占优势的程度和抗震性能,将预测对象分为78种结构类型,其中建筑物40种,其它结构类型38种;根据设施的经济功能,考虑美国国务部的4位数标准工业分类(Four-digit Standard Industrial Classification)中所列出的全部设施类型,分为35种。ATC-13方法虽然考虑全面,但由于该方法完全针对美国建筑特点,主要体现美国的灾害损失特点,与中国的实际建筑结构特点和用途有较大差异,其确定的建筑分类结果很难直接用于具有明显本土化的中国建筑类型。Spence等(2008)针对保险和再保险行业,结合建筑结构类型和建筑使用功能进行分类,将建筑分为28种结构类型和14种用途。然而,上述分类研究都未涉及中国建筑设施的特点,所获得的建筑分类模型难以直接对中国的建筑进行恰当分类。
1976年唐山地震后,中国开始进行地震风险预测并相应颁布了诸多规范,不同的规范给出了诸多建筑分类方式,但各类规范(中华人民共和国建设部等,2004;中华人民共和国住房和城乡建设部等,2010;中华人民共和国国家质量监督检验检疫总局等,2015)面向不同目的,导致了其分类的侧重点不同。建筑抗震设计规范所关心的是建筑的可能破坏形式,相同结构类型的建筑在遭遇地震时,其破坏形式有一定的相似性。因此,在抗震减灾相关规范中,一般按照建筑结构类型划分为6类。地震发生后,各省市地震局及相关部门将对灾害损失进行快速评估(国家地震局等,1996)。然而由于各地地理环境及经济水平差异,各地的建筑结构类型也不尽相同。考虑到震后快速评估工作对评估时效性有一定的要求,灾评报告通常将建筑结构相似、建造价值相近的建筑划分为同一类建筑,并对一部分建筑的分类进行简化处理,因此该分类不能满足保险公司对数据的精确度要求。此外,随着经济不断发展,建筑的种类不断增加,建造水平的提高以及功能区分的逐渐细化使得建筑的现状与历史震害中的建筑资料差异增大。因此,需要及时地对建筑资料进行更新和分类。
除各类规范对建筑物进行了分类外,中国许多学者也根据不同的研究需求对建筑物进行了相应的分类。1995年,中国地震局开展适合中国国情的震害预测和损失评估工作,发展和提出了一些建筑物分类和相应的地震易损性模型。尹之潜(1996)为研究建筑设施的易损性,将中国建筑物分为21类;徐敬海等(2002)从震害预测的角度,将房屋结构分为多层砖房、单层工业厂房、钢筋混凝土多层及高层楼房、内框架和底框架房屋、空旷房屋和老旧民房。然而,地震灾害发生后,不同用途的相同结构建筑所造成的灾害经济损失和人员伤亡也将有显著差异,这也直接影响灾后损失评估和保险理赔的效率和准确度。另外一些学者采用遥感技术对建筑结构进行分类,根据建筑外观、光谱特征、尺寸差异、灰度、表面光滑度等将建筑结构分类,如徐偲等(2014)将小型城镇建筑与农田、道路和水体进行了区分。该分类方法限于技术上的精确度,无法满足巨灾保险的数据要求。因此,符合中国国情的建筑物分类对于地震风险研究及未来地震巨灾保险体系建设十分必要。为满足地震风险分析,建立地震巨灾保险体系,需要建立更详细且更符合中国建筑特性的建筑分类清单。
地理信息系统(GIS)在地震风险研究中有广泛的应用,邬亲敏(2004)、吕红山(2005)、陈洪富(2012)等结合GIS技术,建立了相关地震信息数据库,并对地震风险进行了研究。结合GIS技术,可对地震承灾体的数据快速、精确地定位,并提供动态的分析,节省了大量的人力、物力。网络爬虫可使研究人员获取承灾体属性数据以及矢量数据,解决了GIS中数据获取的难点。
本文运用数字地图接口(API开发)结合地理信息系统(GIS)和网络爬虫技术等手段,对建筑进行分类,旨在提供不同城市或区域的建筑信息,该方法可避免大量人工实地调查,以期解决传统分类在针对具体目标区域时分类指标体系不具备可变性等问题。
1. 建筑分类方法
传统建筑结构分类主要为抗震设计服务,偏重于研究各类建筑结构的破坏方式,而较少考虑导致的建筑及相关经济损失比例以及造成的社会影响方面的因素。结合建筑结构类型和使用功能,并考虑中国的建筑现状,本文建立1种多因素耦合的建筑分类模型,通过统计分析现有的土地规划、建筑设计规范、历史震害资料以及中国建筑结构特性,设计基于结构类型-使用功能的建筑分类;通过参考历史灾评报告以及现阶段可获取的建筑易损性信息,建立面向地震保险业务的建筑结构-使用功能分类。形成的建筑分类既可体现时代特征,又可结合区域特点;既考虑了建筑的结构特性,同时也对建筑的使用功能进行了考察。
利用建立的建筑分类模型,通过地图API开发、GIS空间分析以及网络爬虫技术提取相关信息,形成具有区域特点的结构-使用功能区域化分类,为地震风险分析、地震保险研究提供参考依据。
1.1 中国建筑的主要结构类型及使用功能
现有能收集到的各类资料来源主要为相关建筑规范、各地灾评报告以及震后现场工作报告(中国地震局监测预报司,2001;中国地震局震灾应急救援司, 2006, 2011)。结构类型-使用功能分类从建筑的结构类型及其使用功能2方面对中国现存的建筑进行分类。其目的是在考察建筑结构破坏情况的同时,对其所造成的经济损失进行统计。由于中国各地经济发展不均衡且气候差异较大,各地的建筑类型有较大的差异,此外,不同结构类型的抗震能力存在差异,且相同结构类型的建筑物也有不同的用途。结合中国现有建筑规范、灾评报告以及现场实地调查报告,总结出常用建筑的结构类型、地域分布及其用途,见表 1。
表 1 常用建筑的结构类型、地域分布及其用途Table 1. Empirical distribution and function of buildings with different structure types结构类型 建筑分布 建筑用途 木结构 古建筑保护区、南方乡村等 居住、礼制、陵墓等 砖石结构 古建筑保护区、中部、西北地区乡村等 居住、礼制、陵墓 砌体结构配筋砌体 全国 居住、办公、工业等 砌体结构无筋砌体 全国,村镇自建房居多 居住、工业、围墙等 框架结构 全国 居住、办公、商业、工业 剪力墙结构 全国 居住、办公 框支剪力墙结构 全国 底层商业的综合体 框架-剪力墙结构 全国 办公、商业、娱乐、酒店、工业等 筒体-框架结构 全国 办公、商业 壳体结构 全国 商业、交通、工业 网架结构 全国 商业、交通、工业 1.2 基于结构类型-使用功能的建筑分类模型
随着社会的发展,建筑的使用功能更专一和细化,从而形成了更多的用途类型,本文将采用更细化的分类方法。为此,在建筑结构分类的基础上,考虑建筑使用功能分类,并进一步研究基于结构类型-使用功能的建筑分类,如图 1所示。基于结构类型-使用功能的建筑分类旨在从建筑的结构类型和使用功能2方面对中国现存的建筑进行分类。在灾害性地震发生后,利用该分类可对灾区建筑结构的破坏情况进行调查,也可对所造成的经济损失进行合理评估。
上述建筑结构类型-使用功能分类模型是涉及全部建筑结构和建筑使用功能的组合分类模型,可形成众多的分类类型。ATC-13中提到“分类清单越详细,结果越不准确”,其原因是清单越详细,可变性越差,当1个给定的分类表格与调查区域资料出现不一致时,用既定的数据进行预测或评估将出现较大的偏差。建筑结构类型-使用功能分类模型的“分类清单”较为详细,很难由传统方法实现。为此,需要借助地图API和GIS等技术,为基于结构类型和使用功能的建筑分类提供灵活的、可及时更新的建筑信息基础。依据结构类型-使用功能的建筑分类采用地图API开发、GIS空间分析等手段,可生成具有区域特点的建筑分类形式。
2. 分类模型的构建流程
构建区域化结构类型-使用功能分类模型的典型技术流程如图 2所示。从技术流程可以看出,建筑结构和用途等方面信息的提取是建筑分类的关键基础。信息提取系统采用B/S架构,服务端与客户端采用WEBGIS技术展现地理基础数据和数据查询等交互式访问,可进行大批量的空间查询。使用互联网地图免费提供的API获取城市及兴趣点(point of interest,简称POI)信息,可获取覆盖全国城镇的建筑名称、地理位置信息、关系信息以及用途信息。进一步通过网络爬虫技术抓取链家官网[2]数据,获取商务住宅房屋结构信息。目前,链家官网已覆盖全国32个城市,可提供北京、上海、深圳及其辐射区域内的完备的建筑基础信息。通过百度API获取到现有互联网地图上所有可能的建筑用途,共19大类,869小类,见表 2。
表 2 百度API获取的建筑物POI分类Table 2. Building POI classification from Baidu APIPOI大类 POI小类 政府机构
及社会团体社会团体、交通车辆管理、公检法机构等 科教文化服务 科研机构、学校等 交通设施服务 停车场、火车站、长途汽车站、港口码头、地铁站、公交车站等 金融保险服务 保险公司、银行、金融保险服务机构、自动提款机等 公司企业 农林牧渔基地、公司、公司企业等 住宿服务 住宿服务相关、旅馆招待所、宾馆酒店等 体育休闲服务 休闲场所、影剧院、娱乐场所、运动场馆等 道路附属设施 收费站、服务区、路牌信息、警示信息等 公共设施 紧急避难场所、公共厕所 通行设施 临街院门、通行设施 商务住宅 住宅区、楼宇、商业住宅相关、产业园区等 风景名胜 风景名胜、公园广场等 医疗保健服务 专科医院、综合医院、医疗保健服务场所、疾病预防机构、诊所等 生活服务 洗衣店、物流速递、中介机构、电讯营业厅、邮局等 购物服务 家电电子卖场、商场、便利店、超级市场等 餐饮服务 休闲餐饮场所、咖啡厅、快餐店、中餐厅等 汽车销售 汽车销售等 汽车维修 汽车维修、汽车养护、加油站等 摩托车服务 摩托车销售 基于建筑经纬度、结构类型及用途等信息,首先对建筑进行结构类型和使用功能分类,随后进一步对建筑的结构类型与使用功能耦合处理,得到建筑的结构类型-使用功能分类。此方式便于分库存储,也可根据业务需求转换为对应的字段关系,提高查询性能。
以北京为例,建筑结构功能性空间分布信息如图 3所示,采用地图瓦片的形式,对每个图层进行属性设置,进而构成了多层次的建筑属性空间分布。图 3中C层的指标为POI中筛选出的详细分类,属于D1商业分类指标的子集;D层的指标为高于C层的指标分类,所包含的内容更广泛。通过不同指标的选择,可筛选所需分类下的信息数据。考虑到地震保险关心单个投保标的损失,对于建筑物的空间表达,参考吕红山(2005)的表达方式,采用多边形特征来表示建筑物的轮廓。
此外,在互联网地图中,建筑所有者会因自身需求对建筑信息进行更新。因此,使用互联网地图接口的建筑结构类型-使用功能分类,其建筑结构可随互联网地图进行更新,克服了传统建筑分类由于不可变性和难以更新性造成的分类落后于时代的缺点。
3. 分类模型的应用示例
本文共收集到北京市的相关建筑10余万栋,形成了北京市建筑结构数据库。通过网络爬虫技术,抓取链家官网信息得到建筑结构信息,建筑结构类型被简单地分为钢混结构、砖混结构和混合结构,其中混合结构包含了钢混结构和砖混结构以外的建筑结构类型。以北京市朝阳区某街区A与北京市西城区某街区B为例,按照建筑结构分类,对街区的建筑信息进行展示,如图 4、图 5所示。
图 4、图 5是建筑矢量数据按结构类型-使用功能耦合属性的表达,其中,类型相同的建筑用同色系色块表示。从图中可以看出,2个街区的建筑类型有显著的不同,朝阳区A街区建筑以钢混结构的居民建筑为主,西城区B街区建筑多为混合结构楼房,建筑用途类型较为多样,兼具教育、商业和居住等用途。其原因是朝阳区所选街区大部分为2000年以后的新建民宅小区,结构相对先进;西城区所选街区属于老旧城区,建筑年代多为1980年以前,所选区域的开发相对较早,因此建筑用途相对具有多样性,而建筑结构则较为单一,多为混合结构。随着建筑水平的提高和建筑材料的改进,2000年以后的小区多为钢混结构,而20世纪80年代以前的建筑多为混合结构。通过分析,得到了按结构-使用功能分类的北京市建筑分布比例,见表 3。
表 3 按结构-使用功能分类的北京市建筑分布比例(单位:%)Table 3. Proportions of buildings with coupling classification in Beijing (unit:%)POI大类 钢混结构
比例混合结构
比例砖混结构
比例商务住宅 92.58 93.32 85.28 政府机构及
社会团体0.29 0.17 0.66 公司企业 0.31 0.13 0.66 风景名胜 0.15 0.09 0.51 购物服务 1.54 0.26 2.15 科教文化服务 0.43 0.95 1.58 住宿服务 0.60 0.26 0.76 餐饮服务 0.44 0.43 1.49 交通设施服务 1.47 2.72 2.72 通行设施 1.00 0.91 1.65 生活服务 0.78 0.43 1.52 体育休闲服务 0.15 0.17 0.22 医疗保健服务 0.19 0.17 0.60 金融保险服务 0.05 0 0.06 汽车服务 0.02 0 0.06 汽车维修 0 0 0 公共设施 0.01 0 0.06 摩托车服务 0 0 0 汽车销售 0 0 0 依据地图瓦片的形式,在WEBGIS平台进行建筑结构类型-建设年代信息的筛选提取,结果如图 6所示。经过筛选,风险研究人员和保险公司可更直观地认知承保区域的建筑价值及抗震能力。
4. 结论
综合考虑结构类型和使用功能,对建筑进行分类,形成了1种基于结构类型-使用功能的建筑分类模型,该分类模型适用于中国地震风险损失评估和地震保险系统。利用网络爬虫、互联网地图API和WEBGIS相结合的信息提取方法,使得建筑分类的形式具有区域可变性和易更新性,所获得的建筑分类信息更加准确且更符合时代特点。此外,结合建筑使用功能的分类,考虑了建筑物的社会经济价值,在进行灾害损失评定时,对间接经济损失的确定提供了依据。以北京为例对本文的建筑分类进行了应用,结果显示出建筑结构类型-使用功能分类模型的有效性和实用性,分类结果可提供建筑的结构类型、建设年代和用途等信息,其使用功能可形成多达19大类、869小类的分类数据,为震后建筑结构破坏和经济损失的快速评估以及地震保险业务提供基础。
但从初步的应用看,建筑结构和用途的数据来源存在一定的缺陷,如链家网站对中小型城市及乡镇的建筑数据覆盖不全面、建筑结构分类不够详细、山区的互联网地图可能存在更新不及时等情况,均是制约建筑结构类型-使用功能分类方法应用的重要因素。因此,需要进一步挖掘建筑信息数据源和数据提取技术,建立可快速更新、动态管理的城乡建筑信息数据库,以提高建筑分类的全面性、准确性及分类结果的时效性和可靠性。
目前,地震工程对于地震保险研究的主要贡献是为地震风险评估模型提供建筑的结构易损性曲线或震害矩阵,侧重点在于结构的破坏状态和结构的抗震性能,较少考虑建筑使用功能对于易损性的影响,然而使用功能对于易损性的影响不能忽略。直接经济损失往往受建筑价值或重置价值影响,间接经济损失则由建筑的使用功能及其遭受的影响来决定。因此,同时考虑结构类型和使用功能的建筑物分类方法为研究地震保险中建筑物的直接和间接损失提供了一定的理论基础。
-
表 1 运用小震资料求得的唐山地震序列各段断层面走向、倾角、标准差和位置
Table 1. Fault plane parameters determined by using small earthquake for segments in Tangshan earthquake
断层名 小震个数 走向/° 倾角/° 顶点位置(纬度/°N,经度/°E,深度/km) 数据来源 值 标准差 值 标准差 宁河断裂段 61 246.6 4.0 81.8 4.2 (39.32,117.94,3.8),(39.35,117.93,21.7)
(39.29,117.75,21.7),(39.27,117.77,3.8)本文结果 33 253.3 3.9 66.0 5.0 (39.31,117.96,10.0),(39.37,117.94,24.6)
(39.32,117.74,24.6),(39.27,117.76,10.0)万永革等(2008) 唐山断裂南段 250 213.4 0.8 81.9 1.5 (39.57,118.18,3.8),(39.58,118.15,22.0)
(39.33,117.94,22.0),(39.31,117.96,3.8)本文结果 98 210.1 1.2 73.7 2.8 (39.57,118.18,6.4),(39.59,118.13,22.9)
(39.33,117.94,22.9),(39.31,117.99,6.4)万永革等(2008) 唐山断裂北段 1646 231.4 0.3 89.1 0.6 (39.78,118.49,4.0),(39.78,118.49,19.0)
(39.60,118.19,19.0),(39.60,118.19,4.0)本文结果 665 233.1 0.5 89.1 1.3 (39.78,118.49,7.7),(39.78,118.49,21.8)
(39.61,118.19,21.8),(39.61,118.19,7.7)万永革等(2008) 卢龙断裂段 694 46.1 0.6 89.3 1.5 (39.86,118.82,4.0),(39.86,118.83,16.9)
(39.72,118.62,16.9),(39.71,118.62,4.0)本文结果 176 39.0 0.9 86.7 1.3 (39.72,118.62,7.8),(39.72,118.63,22.6)
(39.86,118.78,22.6),(39.86,118.77,7.8)万永革等(2008) 滦县断裂段 404 125.1 1.6 76.2 1.8 (39.75,118.70,4.2),(39.73,118.68,13.5)
(39.67,118.80,13.5),(39.70,118.81,4.2)本文结果 160 118.4 1.9 76.9 2.0 (39.75,118.70,8.4),(39.73,118.68,20.3)
(39.68,118.80,20.3),(39.70,118.81,8.4)万永革等(2008) 表 2 地震断层面参数研究结果差异统计
Table 2. Statistical results of the differences in seismic fault plane parameters
地震事件 主节面/° 走向差异/° 倾角差异/° 数据来源 走向 倾角 最大值 最小值 最大值 最小值 2001年2月23日四川雅江MS 6.0地震 123 25 31 6 17 8 龙思胜(2004) 117 42 美国地质勘探局(USGS) 92 34 哈佛大学(HRV) 2003年9月27日中俄蒙边界MS 7.9地震 131 71 4 1 14 6 全球矩心矩张量(gCMT) 130 85 美国地质勘探局(USGS) 127 79 赵翠萍等(2005) 2006年12月26日中国台湾南部滨海MS 7.2地震 330 58 11.5 4 25 5.1 全球矩心矩张量(gCMT) 334 83 美国地质勘探局(USGS) 341.5 77.9 郭志等(2008) 2008年5月12日四川汶川MS 8.0地震 357 68 23 2 22 2 全球矩心矩张量(gCMT) 15 60 美国地质勘探局(USGS) 7 55 地球物理研究所(CEA-IGP) 5 48 郭祥云等(2010) 352 70 郑勇等(2009) 2008年10月5日新疆天山-帕米尔MS 6.7地震 82 53 24.8 7.8 10.3 2.3 全球矩心矩张量(gCMT) 65 45 美国地质勘探局(USGS) 57.2 42.7 苏金蓉等(2013) 2010年4月14日青海玉树MS 7.1地震 300 88 6.4 1 10 0 全球矩心矩张量(gCMT) 301 86 美国地质勘探局(USGS) 209 88 中国地震台网中心(CENC) 294.6 78 盛书中等(2014) 2012年5月3日甘肃金塔MS 5.4地震 162 80 13 1 16 2 全球矩心矩张量(gCMT) 163 74 地球物理研究所(CEA-IGP) 159 78 地震预测研究所(CEA-IES) 172 64 张辉等(2012) 2013年4月20日四川芦山MS 7.0地震 210 38 12 0 14 1 全球矩心矩张量(gCMT) 198 33 美国地质勘探局(USGS) 210 47 地震预测研究所(CEA-IES) 216 45 谢祖军等(2013) 209 46 吕坚等(2013) 2013年8月28日云南德钦MS 5.9地震 292 43 14 7 10 0 全球矩心矩张量(gCMT) 285 53 地球物理研究所(CEA-IGP) 299 53 罗钧等(2015) 2014年2月12日新疆于田MS 7.3地震 242 82 3 0 4 0 美国地质勘探局(USGS) 239 82 中国地震局(CEA) 242 78 中国地震台网中心(CENC) 2014年5月30日云南盈江MS 6.1地震 80 83 20 2 3 0 美国地质勘探局(USGS) 82 85 全球矩心矩张量(gCMT) 260 82 地球物理研究所(CEA-IGP) 85 83 赵旭等(2014) 2014年8月3日云南鲁甸MS 6.5地震 162 86 2 0 4 0.9 美国地质勘探局(USGS) 160 90 全球矩心矩张量(gCMT) 160 87 地球物理研究所(CEA-IGP) 160 89.1 刘丽芳等(2014) 2016年1月21日青海省门源县MS 6.4地震 143 35 34 2 11 2 中国地震台网中心(CENC) 141 38 地球物理研究所(CEA—IGP) 134 43 哈佛大学(HRV) 157 34 李启雷(2016) 129 45 李晓峰(2017) 2017年8月8日四川九寨沟MS 7.0地震 153 84 6 0 10 1 美国地质勘探局(USGS) 150 78 全球矩心矩张量(gCMT) 156 79 易桂喜等(2017) 150 80 杨宜海等(2017) 152 74 谢祖军等(2018) 2017年11月18日西藏米林M 6.9地震 132 55 8.3 1 12.4 2.4 美国地质勘探局(USGS) 124.7 59 地震预测研究所(CEA-IES) 133 46.6 吴宝峰(2017) 2018年9月8日云南墨江5.9级地震* 129 81 6 1.8 2 1 赵博等* 123 79 郭祥云等* 124.8 80 地震预测研究所(CEA-IES) 2018年9月12日陕西宁强5.3级地震* 171 67 6 4.1 18 0 中国地震台网中心台网部应急组* 165 85 赵博等* 165 85 郭祥云等* 169.1 78 地震预测研究所(CEA-IES) 2018年9月28日西藏日土5.1级地震* 323 71 28 5 28 2 中国地震台网中心台网部应急组* 304 76 赵博等* 295 61 郭祥云等* 318 59 地震预测研究所(CEA-IES) 309 48 地球物理研究所(CEA—IGP) 注:加*地震数据来源于中国地震台网中心的地震监测人微信公共号及微信公众平台。 -
陈运泰, 林邦慧, 王新华等, 1979.用大地测量资料反演的1976年唐山地震的位错模式.地球物理学报, 22(3):201-217. doi: 10.3321/j.issn:0001-5733.1979.03.001 杜晨晓, 谢富仁, 张扬等, 2010. 1976年MS 7.8唐山地震断层动态破裂及近断层强地面运动特征.地球物理学报, 53(2):290-304. doi: 10.3969/j.issn.0001-5733.2010.02.007 郭祥云, 陈学忠, 李艳娥, 2010. 2008年5月12日四川汶川8.0级地震与部分余震的震源机制解.地震, 30(1):50-60. http://d.old.wanfangdata.com.cn/Periodical/diz201001006 郭志, 高星, 王卫民, 2008. 2006年12月26日台湾南部滨海MS 7.2级地震破裂过程研究.地球物理学报, 51(4):1103-1113. doi: 10.3321/j.issn:0001-5733.2008.04.019 虢顺民, 李志义, 程绍平等, 1977.唐山地震区域构造背景和发震模式的讨论.地质科学, 12(4):305-320. http://www.cnki.com.cn/Article/CJFDTotal-DZKX197704000.htm 胡新亮, 刁桂苓, 高景春等, 2001.华北、西南一些地区地震发生在地壳浅部的证据.地震学报, 23(4):427-435. doi: 10.3321/j.issn:0253-3782.2001.04.011 胡幸平, 崔效锋, 2013.华北地区中部地震精定位与构造应力场研究.震灾防御技术, 8(4):351-360. doi: 10.3969/j.issn.1673-5722.2013.04.002 黄文辉, 康英, 苏柱金等, 2016.全国统一编目系统设计与实现.地震地磁观测与研究, 37(4):170-175. http://d.old.wanfangdata.com.cn/Periodical/dzdcgcyyj201604031 李红光, 王利亚, 孙刚等, 2015.华北地区中小地震重新定位和地震活动特征研究.地震, 35(1):28-37. http://d.old.wanfangdata.com.cn/Periodical/diz201501004 李启雷, 李玉丽, 马丽等, 2016. 2016年青海门源MS 6.4地震震源机制与震源深度计算.地震研究, 39(S1):55-61. http://www.cqvip.com/QK/92873A/201704/673165989.html 李钦祖, 张之立, 靳雅敏等, 1980.唐山地震的震源机制.地震地质, 2(4):59-67. http://d.old.wanfangdata.com.cn/Periodical/dizhen201604009 李晓峰, 2017. 2016年1月21日青海门源MS 6.4地震震源机制解及发震构造初步探讨.地震工程学报, 39(4):657-661. doi: 10.3969/j.issn.1000-0844.2017.04.0657 李迎秋, 刘鑫, 万永革, 2011.利用智利地震震源分布确定智利地区主震断层面.防灾科技学院学报, 13(2):51-56. doi: 10.3969/j.issn.1673-8047.2011.02.010 刘白云, 袁道阳, 张波等, 2012. 1879年武都南8级大地震断层面参数和滑动性质的厘定.地震地质, 34(3):415-424. doi: 10.3969/j.issn.0253-4967.2012.03.003 刘丽芳, 徐甫坤, 2014.利用精确定位余震资料确定2014年云南鲁甸6.5级地震的断层面参数.地震研究, 37(4):489-494. doi: 10.3969/j.issn.1000-0666.2014.04.001 龙思胜, 2004. 2001年四川雅江6.0级地震序列的破裂特征及发震构造.中国地震, 20(1):1-11. doi: 10.3969/j.issn.1001-4683.2004.01.001 罗钧, 赵翠萍, 周连庆, 2015. 2013年8月香格里拉德钦-得荣MS 5.9地震序列震源机制与应力场特征.地球物理学报, 58(2):424-435. http://www.cnki.com.cn/Article/CJFDTotal-DQWX201502007.htm 吕坚, 王晓山, 苏金蓉等, 2013.芦山7.0级地震序列的震源位置与震源机制解特征.地球物理学报, 56(5):1753-1763. http://d.old.wanfangdata.com.cn/Periodical/dqwlxb201305033 潘睿, 盛书中, 万永革, 2015. 2009年姚安地震断层面参数的确定.华北地震科学, 33(2):63-68. doi: 10.3969/j.issn.1003-1375.2015.02.011 盛书中, 万永革, 王未来等, 2014. 2010年玉树MS 7.1地震发震断层面参数的确定.地球物理学进展, 29(4):1555-1562. 苏金蓉, 郭志, 2013.西南天山-帕米尔2008年10月5日MS 6.7级地震震源机制研究.地球物理学报, 56(2):504-512. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxb201302014 万永革, 沈正康, 刁桂苓等, 2008.利用小震分布和区域应力场确定大震断层面参数方法及其在唐山地震序列中的应用.地球物理学报, 51(3):793-804. doi: 10.3321/j.issn:0001-5733.2008.03.020 王椿镛, 吴庆举, 段永红等, 2017.华北地壳上地幔结构及其大地震深部构造成因.中国科学:地球科学, 47(6):684-719. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgkx-cd201706004 王福昌, 靳志同, 钱小仕等, 2012.由余震分布确定大地震子断层及其参数的模糊聚类方法.地震学报, 34(6):793-803. doi: 10.3969/j.issn.0253-3782.2012.06.006 王福昌, 万永革, 钱小仕等, 2013.由地震分布丛集性给出断层参数的一种新方法.地球物理学报, 56(2):522-530. http://d.old.wanfangdata.com.cn/Periodical/dqwlxb201302016 王景明, 郑文俊, 陈国顺等, 1981.唐山地震地面主破裂带及地震成因探讨.地震研究, 4(4):437-450. http://www.cnki.com.cn/Article/CJFDTotal-DZYJ198104010.htm 王鸣, 王培德, 1992. 1989年10月18日大同-阳高地震的震源机制和发震构造.地震学报, 14(4):407-415. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=HY000001807747 吴宝峰, 2017. 2017年11月18日西藏米林6.9级地震震源机制.地震地磁观测与研究, 38(6):26-29. doi: 10.3969/j.issn.1003-3246.2017.06.005 谢祖军, 金笔凯, 郑勇等, 2013.近远震波形反演2013年芦山地震震源参数.中国科学:地球科学, 43(6):1010-1019. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgkx-cd201306008 谢祖军, 郑勇, 姚华建等, 2018. 2017年九寨沟MS 7.0地震震源性质及发震构造初步分析.中国科学:地球科学, 48(1):79-92. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgkx-cd201801006 杨超群, 孟凡顺, 万永革, 2013.采用精确定位小震资料反演伽师地震断层面.地球物理学进展, 28(6):2865-2871. http://d.old.wanfangdata.com.cn/Periodical/dqwlxjz201306007 杨宜海, 范军, 花茜等, 2017.近震全波形反演2017年九寨沟M 7.0地震序列震源机制解.地球物理学报, 60(10):4098-4104. doi: 10.6038/cjg20171034 易桂喜, 龙锋, 梁明剑等, 2017. 2017年8月8日九寨沟M 7.0地震及余震震源机制解与发震构造分析.地球物理学报, 60(10):4083-4097. doi: 10.6038/cjg20171033 尤惠川, 徐锡伟, 吴建平等, 2002.唐山地震深浅构造关系研究.地震地质, 24(4):571-582. doi: 10.3969/j.issn.0253-4967.2002.04.012 于湘伟, 张怀, 陈运泰, 2010.华北地区地震重新定位结果分析.大地测量与地球动力学, 30(2):29-33. http://d.old.wanfangdata.com.cn/Periodical/dkxbydz201002007 张宏志, 刁桂苓, 陈祺福等, 2008. 1976年唐山7.8级地震震区现今地震震源机制分析.地震研究, 31(1):1-6. http://d.old.wanfangdata.com.cn/Periodical/dzyj200801001 张辉, 王熠熙, 2012. 2012年5月3日金塔-阿拉善盟5.4级地震震源机制解.西北地震学报, 34(2):205-206. http://d.old.wanfangdata.com.cn/Periodical/xbdzxb201202017 张之立, 李钦祖, 谷继成等, 1980.唐山地震的破裂过程及其力学分析.地震学报, 2(2):111-129. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=HY000001790762 赵博, 高原, 石玉涛, 2013.用双差定位结果分析华北地区的地震活动.地震, 33(1):12-21. doi: 10.3969/j.issn.1000-3274.2013.01.002 赵翠萍, 陈章立, 郑斯华等, 2005. 2003年9月27日中、俄、蒙边界MS 7.9地震震源机制及破裂过程研究.地震学报, 27(3):237-249. doi: 10.3321/j.issn:0253-3782.2005.03.001 赵旭, 黄志斌, 房立华, 2014. 2014年云南盈江MS 6.1地震震源机制研究.中国地震, 30(3):462-473. doi: 10.3969/j.issn.1001-4683.2014.03.019 郑勇, 马宏生, 吕坚等, 2009.汶川地震强余震MS(≥ 5.6)的震源机制解及其与发震构造的关系.中国科学D辑:地球科学, 39(4):413-426. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgkx-cd200904003 Butler R., Stewart G. S., Kanamori H., 1979. The July 27, 1976 Tangshan, China earthquake-a complex sequence of intraplate events. Bulletin of the Seismological Society of America, 69(1):207-220. http://d.old.wanfangdata.com.cn/OAPaper/oai_pubmedcentral.nih.gov_3592653 Nábělek J., Chen W. P., Ye H., 1987. The Tangshan earthquake sequence and its implications for the evolution of the North China basin. Journal of Geophysical Research, 92(B12):12615-12628. doi: 10.1029/JB092iB12p12615 Ouillon G., Ducorbier C., Sornette D., 2008. Automatic reconstruction of fault networks from seismicity catalogs:three-dimensional optimal anisotropic dynamic clustering. Journal of Geophysical Research, 113(B1):B01306. http://d.old.wanfangdata.com.cn/OAPaper/oai_arXiv.org_physics%2f0703084 Wang X. S., Feng X. D., Xu X. W., et al., 2014. Fault plane parameters of Sanhe-Pinggu M 8 earthquake in 1679 determined using present-day small earthquakes. EarthqButler R., Stewart G. S., Kanamori H., 1979. The July 27, 1976 Tangshan, China earthquake-a complex sequence of intraplate events. Bulletin of the Seismological Society of America, 69(1):207-220. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dizhen-e201406004 Nábělek J., Chen W. P., Ye H., 1987. The Tangshan earthquake sequence and its implications for the evolution of the North China basin. Journal of Geophysical Research, 92(B12):12615-12628. doi: 10.1029/JB092iB12p12615 Ouillon G., Ducorbier C., Sornette D., 2008. Automatic reconstruction of fault networks from seismicity catalogs:three-dimensional optimal anisotropic dynamic clustering. Journal of Geophysical Research, 113(B1):B01306. http://d.old.wanfangdata.com.cn/OAPaper/oai_arXiv.org_physics%2f0703084 Wang X. S., Feng X. D., Xu X. W., et al., 2014. Fault plane parameters of Sanhe-Pinggu M 8 earthquake in 1679 determined using present-day small earthquakes. Earthquake Science, 27(6):607-614. doi: 10.1007/s11589-014-0099-3 Zhou C. Y., Diao G. L., Geng J., et al., 2010. Fault plane parameters of Tancheng M 81/2 earthquake on the basis of present-day seismological data. Earthquake Science, 23(6):567-576. doi: 10.1007/s11589-010-0756-0 -