《中国地震动参数区划图》（GB 18306—2015）编图方法在云南地区的应用实践
摘要：[JP2]采用《中国地震动参数区划图》（GB 18306—2015）编图的技术方法，将云南地区划分为02°×02°的网格点，计算了基岩峰值加速度，并转换为一般场地结果，与五代图颁布的结果进行比对分析，同时对蒙自地区某代表场地计算了基岩反应谱并进行分析讨论。结果表明：云南地区绝大部分计算点的50年10%超越概率水平对应的一般场地峰值加速度试算结果与五代图比较吻合。从反应谱计算结果可以看出，五代图的计算方法及结果有效反映了云南地区地震环境对该代表场地的影响。建议今后云南地区地震危险性分析工作应采用五代图编图的技术方法。[JP]
关键词： 中国地震动参数区划图；基岩峰值加速度；超越概率；地震危险性分析
中图分类号：P315911文献标识码：A文章编号：1000-0666（2017）02-0257-07
0前言
从社会经济发展的要求、科技发展的必然结果、新资料的积累、重大灾害性事件的经验总结、国际发展趋势以及防震减灾法的相关要求等方面出发，中国地震局修编了《中国地震动参数区划图》[KG-*2]（GB 18306―2015）[KG-*2]（中国地震局，2015）[KG-*2]（以下简称五代图），于2015年6月1日正式发布为国家强制性标准。云南省范围内的重大项目地震安全性评价、地震小区划等涉及地震危险性分析的相关工作应使用五代图编图原则和技术方法。为了进一步掌握和应用这些新的技术方法，同时与云南地區以前的工程项目安评结果衔接过渡，本文采用五代图编图技术方法进行了云南地区试算实践，并对计算结果进行讨论。
1五代图编图原则和技术方法
五代图修编以50年超越概率10%地震动峰值加速度与50年超越概率2%地震动峰值加速度对应的一般场地峰值加速度值的1/19中的较大者作为编图指标，从方法上保证了抗倒塌设计参数（即罕遇地震动）的风险水平不低于50年超越概率2%。同时提出了极罕遇地震作用（年超越概率10-4的地震动）概念，用于国土规划和地震应急备灾。
五代图修编仍采用概率地震危险性分析方法（CPSHA），依据其基本原理，基于层次潜在震源区模型（周本刚等，2013），对地震活动性模型进行了改进。地震统计区内满足3个基本假定：（1）地震活动的震级满足截断的G-R关系；（2）地震发生符合泊松分布；（3）地震活动在不同潜在震源区之间为非均匀分布，在同一潜在震源区内为均匀分布。
五代图修编的基本技术思路和计算方法概述如下：
（1）以地震带为统计单元，分析地震活动的时间非均匀性，确定未来百年地震发生的概率模型和地震危险性空间相对分布概率模型。对每个统计单元采用分段的泊松过程模型。令N表示统计单元未来t年内地震发生次数的随机变量，根据泊松过程的基本假定，发生n次M>4地震的概率为：[KH*1]
P（N=n）=[SX（]（v4t）nn！[SX）]e-v[HT4.]4[HT7"]t[JY]（1）[KH*1D]
[HJ2.2mm]其中，v4为M>4地震的年平均发生率。该值通过地震带未来百年地震活动趋势预测结果而得到，反映了地震带地震活动水平的时间非均匀性。
统计单元内地震震级概率密度函数为截断的指数函数（胡聿贤，1999）：[KH*1]
fM（m）=[SX（]βexp（-β（m-m0））1-exp（-β（muz-m0））[SX）][JY]（2）[KH*1D]
其中：m0为震级下限；muz为地震带震级上限；β=23b；b即为在统计单元内震级频度关系中的b值，它反映大小地震间的比例关系。当震级小于震级下限和大于震级上限时，概率密度值为零。
（2）在地震带（统计单元）内部划分地震构造区和潜在震源区（周本刚，2016）。潜在震源区内地震危险性是均匀分布的。潜在震源区由几何边界、震级上限和分震级档的地震空间分布函数fi，mj来描述。每个地震构造区补充几个背景性潜在震源，且边界也不完全重叠。
[BW（S][BG（；N][BHDWG1*2，WK15mmZQ，WK140mm，WK15mmYQW][HT5"][CM（22mm]地震研究[CM）]40卷[BG）F][BW）]
[BW（D][BG（；N][BHDWG1*2，WKZQ0W][HT5"]第2期[JZ]
闻满华等：《中国地震动参数区划图》[KG-*2]（GB 18306—2015）编图方法在云南地区的应用实践
[BG）F][BW）]
（3）利用全概率求和原理计算在统计单元内发生一次地震时，场点给定地震动值a的超越概率。基本计算公式为（胡聿贤，1999）：[KH*1]
[JP3]P（A≥a）=∫∫∫∫P（A≥a[JB<1|]E）f（x，y[JB<1|]m） fM（m）[JP]f|x，y（θ）dxdydmdθ[JY]（3）[KH*1D]
其中：P（A≥a[JB<1|]E）指震级为m、震中位置为（x，y）、地震动椭圆衰减长轴方向与正东方向夹角为θ时场点给定地震动值a被超过的概率，该函数由地震动衰减关系确定；f（x，y[JB<1|]m）为给定震级的空间分布函数，该函数可由考虑震级分档情况和潜在震源区的面积得到；f|x，y（θ）为等震线长轴取向概率密度函数，用δ函数表示。在同一潜在震源区内等震线长轴取向概率密度函数相同，不同的潜在震源区该函数可以不同（潘华等，2013；潘华，李金臣，2016）。
（4）利用地震发生次数的分段泊松模型，可以计算得到某个统计单元k对场点的超越概率为（胡聿贤，1999）：[KH*1]
Ptk（A≥a）=1-e-v[HT4.]4[HT7"]tP（A≥a）[JY]（4）[KH*1]
（5）若有K个统计单元对场点有影响，则场点总的超越概率为：[KH*2]
Pt（A≥a）[KG-*2]=[KG-*4]1-∏[DD（]Kk=1[DD）][JB<1*[]1-Ptk（A≥a）[JB>1*]][JY]（5）[KH*1]
（6）分区、分段采用不同的地震动参数衰减关系进行计算（肖亮，俞言祥，2010）。本次修订采用了由中国地震动参数区划图第4工作組研制的衰减关系（俞言祥，汪素云，2004；俞言祥等，2013），分为青藏区、新疆区、东部强震区、中强地震区4个分区，长、短两个方向，M>65、M≤65两段，共16种不同的衰减关系。衰减关系可表示为：[KH*1]
lgY=A+BM+Clg[R+DeEM]±σ[JY]（6）[KH*1D]
式中：Y代表加速度峰值或反应谱值；R为震中距；M为震级；A、B、C、D、E为回归系数，σ为标准差。
2云南地区试算实践结果
采用五代图编制所使用的地震活动性模型参数、潜在震源区模型参数和地震动参数衰减关系，对云南地区划分为02°×02°网格点计算了基岩峰值加速度并转换为一般场地结果，与五代图颁布的结果进行比对分析。同时计算某代表场地基岩反应谱，并对其进行分析。
[JP2]云南地区的地震危险性概率分析方法的应用，涉及青藏地震区的鲜水河—滇东地震带、滇西南地震带、喜马拉雅地震带、藏中地震带、龙门山地震带和华南地震区的右江地震带、长江中游地震带，所涉及的地震带和地震构造区地震活动性模型参数，如表1所示（高孟潭，2015）。所使用的分区基岩峰值加速度衰减关系参数、五代图基岩Ⅰ1类场地峰值加速度与一般场地（Ⅱ类）峰值加速度转换系数（高孟潭，2015），分别见表2、3。[JP]
[KG（0.15mm]对不同的地震带设定不同的衰减关系，分M>65地震、M≤65地震两段对应的基岩峰值加速度衰关系参数（俞言祥，2016），分别计算潜在震源（包括背景性潜源）对目标场点的年超越概率贡献值，对得到的4种值通过概率合成为总的年超越概率贡献曲线，从总的年超越概率贡献曲线上求取50年超越概率分别为10%、2%对应的基岩峰值加速度值，将两者通过转换关系转换为一般场地峰值加速度，然后取10%对应的一般场地峰值加速度值和2%对应的一般场地峰值加速度值的1/19两者中的较大值作为该场地的基本地震动。
为了测试五代图潜源空间分布函数、分区衰减关系、分段计算地震动参数方法在云南地区的实际效果，首先选择青藏区（鲜水河—滇东地震带、滇西南地震带）和中强区（右江地震带）的潜在震源均能产生明显影响的蒙自地区某场地进行计算，场地位置为（2337°N，10338°E）。然后分两区、两段计算各地震带潜在震源对场地的超越概率贡献，合成总贡献的概率曲线，如表4、图1a所示。最后从合成的总贡献的概率曲线上求取50年超越概率63%、10%和2%对应的基岩峰值加速度，通过表3转换为一般场地地震动峰值加速度，如表5所示。场地反应谱计算结果如表6、图1b所示。用五代图计算方法对云南地区02°×02°网格点计算了基岩峰值加速度并转换为一般场地结果，如图2所示。
3结论及讨论
本文选择的蒙自地区某试算场地具有代表性，其周围的鲜水河—滇东地震带、滇西南地震动和右江地震带均对其峰值加速度产生显著影响。如表4所示，2种衰减关系、65级上下两段共4类概率贡献的试算值表现合理，该场地所处中强地震活动区，五代图颁布的峰值加速度为010 g分区，从表4中可以看出，对应峰值100 gal以下的4种概率贡献大致在相同数量级水平，200～300 gal段落主要受两区潜源65级以上段的影响，300 gal以上段落主要受青藏区潜源65级以上段的影响。从合成的总概率贡献曲线（图1a）上读取的50年63%、10%、2%这3个概率水平对应的峰值加速度经转换为一般场地峰值分别是005 g、012 g和019 g，这3个概率土层峰值比例合理，10%概率水平的峰值012 g与场地所处010 g分区匹配较好。从峰值加速度和反应谱计算结果可以看出，五代图的计算方法及结果有效反映了该场地所处的地震环境。
[JP2]从图2可以看出，云南地区绝大部分计算点的50年10%超越概率水平对应的一般场地峰值加速度试算值与五代图（中国地震局，2015）较为吻合。有极少部分地区试算值略低于五代图，这是五代图修编时考虑了部分地区大震近场效应，并且征求地方政府意见后从安全的角度考虑对峰值加速度分区有所提高的结果。图中少量计算值为010 g的区域提高至015 g，东川地区提高幅度最大，其原因为东川以北至巧家处在同一8级潜源内，其模型参数设置需要兼顾到该潜源北段巧家地区的烈度计算结果，致使东川周边地区的计算结果达不到040 g，但在编图时考虑到小江断裂的大震近场效应以及上述安全原则，调整提高至040 g水平。[JP]
通过本文的试算实践可以看出，五代图编图的模型、方法和计算结果有效反映出云南地区地震活动环境的实际状况，以供同行共享参考。建议今后云南地区地震危险性分析工作应采用五代图编图的技术方法。
[HTK]感谢中国地震局地球物理研究所俞言祥研究员、潘华研究员提供潜源模型参数和衰减关系数据以及在使用中给予的悉心指导。
参考文献：
高孟潭2015GB 18306―2015《中国地震动参数区划图》宣贯教材[M].北京：中国标准出版社
胡聿贤1999地震安全性评价技术教程[M].北京：地震出版社
潘华，高孟潭，谢富仁2013新版地震区划图地震活动性模型与参数确定[J].震灾防御技术，8（1）：11-23
潘华，李金臣2016新一代地震区划图的地震活动性模型[J].城市与减灾，（3）：28-33
肖亮，俞言祥2010应用两步法使用圆模型拟合华北地区地震烈度衰减关系[C]//中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集北京：中国地震学会
俞言祥，李山有，肖亮2013为新区划图编制所建立的地震动衰减关系[J].震灾防御技术，8（1）：24-33
俞言祥，汪素云2004青藏高原东北地区水平向基岩加速度峰值与反应谱衰减关系[J].地震学报，26（6）：591-600
俞言祥2016新一代地震区划图地震动参数衰减关系的建立与特点分析[J].城市与减灾，（3）：34-38
中国地震局2015GB 18306―2015《中国地震动参数区划图》[M].北京：中国标准出版社
周本刚，陈国星，高战武，等2013新地震区划图潜在震源区划分的主要技术特色[J].震灾防御技术，8（2）：113-124
周本刚2016新一代地震区划图潜在震源区划分的技术进展[J].城市与减灾，（3）：18-23
关键词： 中国地震动参数区划图；基岩峰值加速度；超越概率；地震危险性分析
中图分类号：P315911文献标识码：A文章编号：1000-0666（2017）02-0257-07
0前言
从社会经济发展的要求、科技发展的必然结果、新资料的积累、重大灾害性事件的经验总结、国际发展趋势以及防震减灾法的相关要求等方面出发，中国地震局修编了《中国地震动参数区划图》[KG-*2]（GB 18306―2015）[KG-*2]（中国地震局，2015）[KG-*2]（以下简称五代图），于2015年6月1日正式发布为国家强制性标准。云南省范围内的重大项目地震安全性评价、地震小区划等涉及地震危险性分析的相关工作应使用五代图编图原则和技术方法。为了进一步掌握和应用这些新的技术方法，同时与云南地區以前的工程项目安评结果衔接过渡，本文采用五代图编图技术方法进行了云南地区试算实践，并对计算结果进行讨论。
1五代图编图原则和技术方法
五代图修编以50年超越概率10%地震动峰值加速度与50年超越概率2%地震动峰值加速度对应的一般场地峰值加速度值的1/19中的较大者作为编图指标，从方法上保证了抗倒塌设计参数（即罕遇地震动）的风险水平不低于50年超越概率2%。同时提出了极罕遇地震作用（年超越概率10-4的地震动）概念，用于国土规划和地震应急备灾。
五代图修编仍采用概率地震危险性分析方法（CPSHA），依据其基本原理，基于层次潜在震源区模型（周本刚等，2013），对地震活动性模型进行了改进。地震统计区内满足3个基本假定：（1）地震活动的震级满足截断的G-R关系；（2）地震发生符合泊松分布；（3）地震活动在不同潜在震源区之间为非均匀分布，在同一潜在震源区内为均匀分布。
五代图修编的基本技术思路和计算方法概述如下：
（1）以地震带为统计单元，分析地震活动的时间非均匀性，确定未来百年地震发生的概率模型和地震危险性空间相对分布概率模型。对每个统计单元采用分段的泊松过程模型。令N表示统计单元未来t年内地震发生次数的随机变量，根据泊松过程的基本假定，发生n次M>4地震的概率为：[KH*1]
P（N=n）=[SX（]（v4t）nn！[SX）]e-v[HT4.]4[HT7"]t[JY]（1）[KH*1D]
[HJ2.2mm]其中，v4为M>4地震的年平均发生率。该值通过地震带未来百年地震活动趋势预测结果而得到，反映了地震带地震活动水平的时间非均匀性。
统计单元内地震震级概率密度函数为截断的指数函数（胡聿贤，1999）：[KH*1]
fM（m）=[SX（]βexp（-β（m-m0））1-exp（-β（muz-m0））[SX）][JY]（2）[KH*1D]
其中：m0为震级下限；muz为地震带震级上限；β=23b；b即为在统计单元内震级频度关系中的b值，它反映大小地震间的比例关系。当震级小于震级下限和大于震级上限时，概率密度值为零。
（2）在地震带（统计单元）内部划分地震构造区和潜在震源区（周本刚，2016）。潜在震源区内地震危险性是均匀分布的。潜在震源区由几何边界、震级上限和分震级档的地震空间分布函数fi，mj来描述。每个地震构造区补充几个背景性潜在震源，且边界也不完全重叠。
[BW（S][BG（；N][BHDWG1*2，WK15mmZQ，WK140mm，WK15mmYQW][HT5"][CM（22mm]地震研究[CM）]40卷[BG）F][BW）]
[BW（D][BG（；N][BHDWG1*2，WKZQ0W][HT5"]第2期[JZ]
闻满华等：《中国地震动参数区划图》[KG-*2]（GB 18306—2015）编图方法在云南地区的应用实践
[BG）F][BW）]
（3）利用全概率求和原理计算在统计单元内发生一次地震时，场点给定地震动值a的超越概率。基本计算公式为（胡聿贤，1999）：[KH*1]
[JP3]P（A≥a）=∫∫∫∫P（A≥a[JB<1|]E）f（x，y[JB<1|]m） fM（m）[JP]f|x，y（θ）dxdydmdθ[JY]（3）[KH*1D]
其中：P（A≥a[JB<1|]E）指震级为m、震中位置为（x，y）、地震动椭圆衰减长轴方向与正东方向夹角为θ时场点给定地震动值a被超过的概率，该函数由地震动衰减关系确定；f（x，y[JB<1|]m）为给定震级的空间分布函数，该函数可由考虑震级分档情况和潜在震源区的面积得到；f|x，y（θ）为等震线长轴取向概率密度函数，用δ函数表示。在同一潜在震源区内等震线长轴取向概率密度函数相同，不同的潜在震源区该函数可以不同（潘华等，2013；潘华，李金臣，2016）。
（4）利用地震发生次数的分段泊松模型，可以计算得到某个统计单元k对场点的超越概率为（胡聿贤，1999）：[KH*1]
Ptk（A≥a）=1-e-v[HT4.]4[HT7"]tP（A≥a）[JY]（4）[KH*1]
（5）若有K个统计单元对场点有影响，则场点总的超越概率为：[KH*2]
Pt（A≥a）[KG-*2]=[KG-*4]1-∏[DD（]Kk=1[DD）][JB<1*[]1-Ptk（A≥a）[JB>1*]][JY]（5）[KH*1]
（6）分区、分段采用不同的地震动参数衰减关系进行计算（肖亮，俞言祥，2010）。本次修订采用了由中国地震动参数区划图第4工作組研制的衰减关系（俞言祥，汪素云，2004；俞言祥等，2013），分为青藏区、新疆区、东部强震区、中强地震区4个分区，长、短两个方向，M>65、M≤65两段，共16种不同的衰减关系。衰减关系可表示为：[KH*1]
lgY=A+BM+Clg[R+DeEM]±σ[JY]（6）[KH*1D]
式中：Y代表加速度峰值或反应谱值；R为震中距；M为震级；A、B、C、D、E为回归系数，σ为标准差。
2云南地区试算实践结果
采用五代图编制所使用的地震活动性模型参数、潜在震源区模型参数和地震动参数衰减关系，对云南地区划分为02°×02°网格点计算了基岩峰值加速度并转换为一般场地结果，与五代图颁布的结果进行比对分析。同时计算某代表场地基岩反应谱，并对其进行分析。
[JP2]云南地区的地震危险性概率分析方法的应用，涉及青藏地震区的鲜水河—滇东地震带、滇西南地震带、喜马拉雅地震带、藏中地震带、龙门山地震带和华南地震区的右江地震带、长江中游地震带，所涉及的地震带和地震构造区地震活动性模型参数，如表1所示（高孟潭，2015）。所使用的分区基岩峰值加速度衰减关系参数、五代图基岩Ⅰ1类场地峰值加速度与一般场地（Ⅱ类）峰值加速度转换系数（高孟潭，2015），分别见表2、3。[JP]
[KG（0.15mm]对不同的地震带设定不同的衰减关系，分M>65地震、M≤65地震两段对应的基岩峰值加速度衰关系参数（俞言祥，2016），分别计算潜在震源（包括背景性潜源）对目标场点的年超越概率贡献值，对得到的4种值通过概率合成为总的年超越概率贡献曲线，从总的年超越概率贡献曲线上求取50年超越概率分别为10%、2%对应的基岩峰值加速度值，将两者通过转换关系转换为一般场地峰值加速度，然后取10%对应的一般场地峰值加速度值和2%对应的一般场地峰值加速度值的1/19两者中的较大值作为该场地的基本地震动。
为了测试五代图潜源空间分布函数、分区衰减关系、分段计算地震动参数方法在云南地区的实际效果，首先选择青藏区（鲜水河—滇东地震带、滇西南地震带）和中强区（右江地震带）的潜在震源均能产生明显影响的蒙自地区某场地进行计算，场地位置为（2337°N，10338°E）。然后分两区、两段计算各地震带潜在震源对场地的超越概率贡献，合成总贡献的概率曲线，如表4、图1a所示。最后从合成的总贡献的概率曲线上求取50年超越概率63%、10%和2%对应的基岩峰值加速度，通过表3转换为一般场地地震动峰值加速度，如表5所示。场地反应谱计算结果如表6、图1b所示。用五代图计算方法对云南地区02°×02°网格点计算了基岩峰值加速度并转换为一般场地结果，如图2所示。
3结论及讨论
本文选择的蒙自地区某试算场地具有代表性，其周围的鲜水河—滇东地震带、滇西南地震动和右江地震带均对其峰值加速度产生显著影响。如表4所示，2种衰减关系、65级上下两段共4类概率贡献的试算值表现合理，该场地所处中强地震活动区，五代图颁布的峰值加速度为010 g分区，从表4中可以看出，对应峰值100 gal以下的4种概率贡献大致在相同数量级水平，200～300 gal段落主要受两区潜源65级以上段的影响，300 gal以上段落主要受青藏区潜源65级以上段的影响。从合成的总概率贡献曲线（图1a）上读取的50年63%、10%、2%这3个概率水平对应的峰值加速度经转换为一般场地峰值分别是005 g、012 g和019 g，这3个概率土层峰值比例合理，10%概率水平的峰值012 g与场地所处010 g分区匹配较好。从峰值加速度和反应谱计算结果可以看出，五代图的计算方法及结果有效反映了该场地所处的地震环境。
[JP2]从图2可以看出，云南地区绝大部分计算点的50年10%超越概率水平对应的一般场地峰值加速度试算值与五代图（中国地震局，2015）较为吻合。有极少部分地区试算值略低于五代图，这是五代图修编时考虑了部分地区大震近场效应，并且征求地方政府意见后从安全的角度考虑对峰值加速度分区有所提高的结果。图中少量计算值为010 g的区域提高至015 g，东川地区提高幅度最大，其原因为东川以北至巧家处在同一8级潜源内，其模型参数设置需要兼顾到该潜源北段巧家地区的烈度计算结果，致使东川周边地区的计算结果达不到040 g，但在编图时考虑到小江断裂的大震近场效应以及上述安全原则，调整提高至040 g水平。[JP]
通过本文的试算实践可以看出，五代图编图的模型、方法和计算结果有效反映出云南地区地震活动环境的实际状况，以供同行共享参考。建议今后云南地区地震危险性分析工作应采用五代图编图的技术方法。
[HTK]感谢中国地震局地球物理研究所俞言祥研究员、潘华研究员提供潜源模型参数和衰减关系数据以及在使用中给予的悉心指导。
参考文献：
高孟潭2015GB 18306―2015《中国地震动参数区划图》宣贯教材[M].北京：中国标准出版社
胡聿贤1999地震安全性评价技术教程[M].北京：地震出版社
潘华，高孟潭，谢富仁2013新版地震区划图地震活动性模型与参数确定[J].震灾防御技术，8（1）：11-23
潘华，李金臣2016新一代地震区划图的地震活动性模型[J].城市与减灾，（3）：28-33
肖亮，俞言祥2010应用两步法使用圆模型拟合华北地区地震烈度衰减关系[C]//中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集北京：中国地震学会
俞言祥，李山有，肖亮2013为新区划图编制所建立的地震动衰减关系[J].震灾防御技术，8（1）：24-33
俞言祥，汪素云2004青藏高原东北地区水平向基岩加速度峰值与反应谱衰减关系[J].地震学报，26（6）：591-600
俞言祥2016新一代地震区划图地震动参数衰减关系的建立与特点分析[J].城市与减灾，（3）：34-38
中国地震局2015GB 18306―2015《中国地震动参数区划图》[M].北京：中国标准出版社
周本刚，陈国星，高战武，等2013新地震区划图潜在震源区划分的主要技术特色[J].震灾防御技术，8（2）：113-124
周本刚2016新一代地震区划图潜在震源区划分的技术进展[J].城市与减灾，（3）：18-23
