郑州市断裂最大潜在地震发震概率评价1

王 继 田勤俭 高占武

1）中国地震灾害防御中心，北京 100029

2）中国地震局地震预测研究所，北京 100036

郑州市断裂最大潜在地震发震概率评价1

王 继1)田勤俭2)高占武1)

1）中国地震灾害防御中心，北京 100029

2）中国地震局地震预测研究所，北京 100036

根据“郑州市城市活断层探测与地震危险性评价（二期）”综合目标区断层最大潜在地震判定结果，郑州市近东西向断裂有老鸦陈断裂和上街断裂等。老鸦陈断裂第四纪不活动，上街断裂等可能具有发生5—5.5级地震的能力。为了得到最大地震的发震概率，划分了统计区及潜在震源区，得到了地震活动性参数及空间分布函数。最终得出了郑州市断裂未来50年发生1次5级以上地震的概率为6%，未来100年发生1次5级以上地震的概率为11%。

郑州市 断裂探测 地震危险性评价

引言

《郑州市城市活断层探测与地震危险性评价（二期）》属于中国地震观测网络工程城市活断层探测项目。2006年9月，中国地震局地震预测研究所与河南省地震局签署了第二阶段郑州市城市活断层探测与地震危险性评价地质部分的技术服务合同。目标区活动断层地震危险性综合评价属于第二阶段中的第四专题“老鸦陈断层地质调查与地震危险性综合评价”的第4个子专题，该项目已通过验收，本文主要内容取自该专题成果。

1 地震危险性评价方法

对于第四纪有活动的断裂，最大震级和地震危险性评价依据有很多方法（中国地震局，2005；闻学泽，2007）。由于断裂所在的目标区附近地震很少，不能依据本身的地震活动判定地震危险性，因此，将地震统计范围扩大到与目标区断裂相关的区域断裂上。通过一定区域范围内地震活动性的统计分析，获取与地震危险性评价相关的b值、地震年发生率等活动性参数，评价相关区域的地震危险性；根据地震构造模型，通过统计区域内相关断裂地震活动比较，评价目标断裂的地震危险性。

2 统计区地震危险性概率分析

统计区位于华北平原地震带的华北南部地震构造区的一部分，该区分布一组近北西走向的断裂系（图1），根据断裂分布、历史地震和现代小地震的分布，可以划分为4个潜在震源区，震级上限从5.5级—6.5级。

（1）中牟断裂（F20）

走向北西西，倾向北北东，倾角65°—70°，正断裂，北盘下降，断距浅层100—200m，深层 200—700m，目标区内长约 36km。中牟断裂西段（花园口断裂以西）切割新生界地层等厚线；东段则没有切割新生界地层等厚线。根据浅层地震反射勘探和钻孔剖面图分析认为，中牟断裂没有错断晚第三纪以来地层等厚线，中牟断裂错断的最新地层为上第三系，没有断错第四纪以来的地层，是一条前第四纪断裂。

（2）上街断裂（F20）

西起荥阳县上街镇附近，向东经郑州市区，一直延伸到中牟县境内，终止于白沙一带，全长约120km，走向近东西，倾向北，倾角70°，为一正断层。

上街断裂的勘探使用了4种手段：地球化学勘探、浅层地震勘探、高密度电法勘探，并利用收集钻孔资料进行跨断裂地层对比。

结果显示，上街断裂并不是一条简单的断裂，它应该是由多个分支断裂组成的断裂带。上街断裂在早第三纪时期仍在活动，但晚第三纪以来已不再活动。因此综合勘探结果认为，该断裂为前第四纪断裂。

（3）老鸦陈断裂（F13）

断裂北起黄河南岸，从邙山东侧地貌陡坎处通过，向南东方向延伸经河南省体育馆东侧穿过市区，并与上街断裂、须水断裂交汇。根据物探结果，老鸦陈断裂是一条走向北北西、倾向北南、倾角60°—70°的正断裂。该断裂错断了新第三纪以前的地层，自新第三纪以来已不再活动。

根据以上结果，判定该潜源区的震级上限为5.5级。

为了给出潜在震源区1的地震活动性参数，将这4个潜在震源区按1个断裂系统来考虑，作为1个统计区来考虑其地震活动性参数，再根据多因子断别方式，分配潜在震源区1的地震活动性参数。

2.1 基本计算公式

统计区包含有多条相同或不同方向的、有中强地震发生能力的断裂或断裂段落，但往往不易分出主、次级断裂，有时也难以进一步细分段。因此，这类统计区内震级≥M破坏性地震的发生时间进程可以近似为泊松过程（Connel，1968）。这种随机的泊松过程是无记忆的，即在时段内地震发生的概率与小区内过去的破坏性地震发生时间、地点无关。从而，在任意给定长度为ΔT的时段内，任一具有发震能力的断裂（或小段）上发生K次震级≥M1破坏性地震的概率P1可由泊松近似（闻学泽，2007）：

图1 统计区及潜在震源区分布Fig. 1 Statistical units and distribution of potential source area

而发生K次震级≥M2的地震破坏性的概率P2为：

因此，设以上提到的震级 M1＞M2，则在相同的ΔT时段内，小区内发生 K次震级在 M1和M2之间的地震概率P，可由概率减法得到：

以上式（1）—式（3）的ΔT是给定的预测时段，可根据需要设为10—200年中的任意长度。式中的V1和V2也可写成V（≥M1）和V（≥M2），分别是小区内震级≥M1和≥M2地震的平均年发生率。去掉脚标后改记为V（≥M），它与地震的平均复发间隔T（≥M）存在如下关系：

对于大小处于小区的震级－频度关系（式7）线性段之内、震级≥M的破坏性地震，有：

或者

式中，a、b来自于同一断裂系的震级－频度关系的参数：

由上述可见，若要估算小区内震级处于震级－频度（式7）线性段（亦称G－R关系）的最大截距震级Mm(=a/b)左侧的破坏性地震的概率，关键是能估计出可靠的、具有长期代表性的G－R关系参数a和b值，接着由式（6）估计出震级≥M地震的年发生率，最后再由式（1）、（2）、（3）估算出破坏性地震的发生概率。

为了综合现代台网地震资料与历史记载的破坏性地震资料来估计断裂系的长期、稳定的G－R关系参数a和b值，同样可考虑由式（8）将不同震级档（≥Mi）的地震在其完整记录（或记录）时段（T0−Ti）中的次数 ni(Mi)归算到一个时间单位（t），得到某震级档 i在地震 t年中理论频度（次数）的估值：

式中，Ti和 T0分别是第Mi震级档地震的完整记录（或记载）的起始和终止年份；Ti则依不同断裂段（小区）、不同震级范围地震的完整记录（或记载）时段的不同而变化。T是不同震级档的地震频数 ni(Mi)，要统一到相同长的时段。同时还应充分考虑人文历史发展的地理差异对历史地震记载的影响，来设定保守的t和Ti的取值方案。

2.2 统计区活动性参数确定

2.2.1 地震资料完整性分析

图2给出了断裂系自1000年以来MS≥4.7级地震的M-T图和应变释放曲线。本区最早的6级地震记载于1524年，5级地震最早记载于公元前519年，自1484年以来MS≥5级地震基本完整可靠（黄玮琼等，1994）。

对于现代地震目录，本区1980年以来处于区域地震台网ML≥2.0（MS1.5）级地震的有效监测范围内（焦远碧等，1990）。

2.2.2 统计区地震资料和震级频度关系b值

从中国地震台网目录和中国强震目录中选出1980—2006年的MS≥1.5级地震以及历史上（包括现代）MS≥4.7级地震资料（表1）。取震级分档间距△T＝0.5，由最小二乘法计算式（4）的参数a和b值，a=1.682，b=0.6521。统计区所用的资料时间段和数据如表2所示，图3给出了统计区震级－频度关系图。

图2 统计区M-T图及应变释放曲线Fig. 2 M-T diagram and the strain release curve in statistical unit

表1 统计区破坏性地震目录（MS≥4.7）（公元前519年—2006年）Table 1 Devastating earthquake catalog (MS≥4.7) (519 B.C. to 2006) in statistical unit

续表

表2 地震资料的选取时间和数据Table 2 Statistical seismic data at different magnitude intervals

2.2.3 统计区年平均发生率、复发间隔和复发概率

基于以上参数，根据式（1）—式（7），计算出震级≥M破坏性地震的平均复发间隔T（年），以及未来50年和100年的发震概率见表3。

4级以上地震年平均发生率（V4）代表未来时间段内地震统计区地震活动的平均水平。由前面得出的b值结合活跃期地震活动水平推算得到，统计区：V4值＝0.1185。

地震年平均发生率V4代表未来100年内统计区地震活动的水平。为了如实地反映地震活动的时空不均匀性，需要将统计区内的地震年平均发生率，按预测结果合理地分配到相应的各潜在震源中去。这里，作者应用震级分档的方法来反映地震活动的时空不均匀性，而且可以避免低估大地震的危险程度。

按地震危险性分析的思路，统计区内mj震级的地震年平均发生率Vmj可以表示为：

式中，β=bln10，b为该统计区的震级—频度关系式中的b值；sinh(x)为双曲正弦函数；为统计区的震级上限，这里取 6.5；m0为起算震级，这里取4级；Δm为震级分挡间隔，这里取为 0.5；mj为分档间隔中心对应的震级值。

根据求出的年平均发生率，根据式（1）—式（7），计算出分震档的平均复发间隔T（年），以及未来50年和100年的发震概率见表3。

图3 统计区震级－频度关系图Fig. 3 Magnitude-frequency curve in statistical unit

表3 统计区地震年平均发生率、平均复发间隔与发震概率Table 3 Average annual rate, recurrence interval and earthquake occurrence probability in statistical unit

3 郑州市断裂最大潜在地震发震概率

在地震统计区内，须把地震统计区各震级档地震的年平均发生率分配给各相应的潜在震源区。本文采用空间分布函数，根据各潜在震源区发生不同震级档地震可能性的大小，对统计区各震级档的地震年平均发生率进行不等权分配。空间分布函数的物理含义是地震统计区内发生1个mj档震级的地震，落在第i个潜在震源区内概率的大小。在同一地震统计区内满足归一条件：

这里n为地震统计区内第mj档潜在震源区的总数。在本研究区中，mj共分成3个震级档，即：确定时，要考虑以下几个方面（胡聿贤，1999）：①使用多因子综合评判方法；②注意所选因子之间的相互独立性，尽量删除那些相互重叠的因子；③各因子赋值标准应是已有震例的总结；④注意因子的地区性和地震强度的差异。

3.1 各潜在震源区的多因子判断

3.1.1 按面积因素分配参数

各潜在震源区面积因素分配因子见表4。

表4 潜在震源区的面积因子Table 4 Area factor of the potential source area

3.1.2 强震复发间隔与构造空段

对于潜在震源区内有6级以上强震重复发生现象存在的地区，可以根据雷建成（1991）的Weibull分布得到的强震重复发生的累积危险率P（t）赋值标准（表5）。对于潜在震源区内6级以下地震按表3的复发间隔计算。

表5 华北地区强震重复发生的累积危险率（部分，据雷建成，1991）Table 5 Cumulative reoccurrence risk of strong earthquakes in North China (from Lei Jiancheng, 1991)

表6 潜在震源区的强震复发间隔与构造空段因子Table 6 Earthquake recurrence interval and structure gap factor in potential source area

3.1.3 根据潜在震源区活动断裂的长度分配参数

各潜在震源区都有一定的第三系以来的活动断裂，以断裂的长度为分配因子，见表7。

表7 潜在震源区活动断裂的长度因子Table 7 Length factor of active fault in potential source area

3.2 各潜在震源区的空间分布函数

根据以上因子判断的参数，得出了各潜在震源区的空间分布函数（表 8）和各震级档的年平均发生率（表9）。各潜在震源区空间分布函数定义为确定方法的原因是，从前述3个因子中选择因子，用等权求各方法来确定其值大小。在进行具体分析时，由所选用的因子（k）单独对统计区内能够发生相应震级档次（mj）地震的潜在震源区（i）赋值（wi,mj,k）；对每一个因子（k）在统计区内归一化，得到因子载荷Qi,mj,k，此时，作者在分配各潜在震源区赋值时已经进行了归一化处理：

由各因子之和得到各震级档的总荷载Ri,mj：

再由总荷载量（Ri,mj）在统计区内的归一化，即可以得到各潜在震源区的空间分布函数fi,mj：

各潜在震源区的地震活动年平均发生率Vi,mj为：

表8 各潜在震源区空间分布函数Table 8 Spatial distribution function of the potential source area

表9 各潜在震源区年平均发生率Table 9 Average annual occurrence rate of each potential source area

3.3 潜在震源区1的年平均发生率、复发间隔和复发概率

基于以上统计区的参数，根据式（1）—式（7），表10给出了潜在震源区1的年平均发生率、复发间隔和复发概率。

表10 潜在震源区1的年平均发生率、复发间隔和发震概率Table 10 Annual average rate, average recurrence interval and earthquake occurrence probability of potential source area #1

综合上述研究结果，可以认为郑州市断裂50年内发生1次5—5.4级地震的概率约为6%；100年内发生1次5—5.4级地震的概率约为11%。因此，该断裂活动微弱，发震概率很低。

4 结论

综合地震危险性分析结果，本文认为郑州市目标区断裂可能具有发生5—5.5级地震的能力，未来50年发生1次5级以上地震的概率为6%，未来100年发生1次5级以上地震的概率为11%。同时，本文还尝试给出了中弱活动断裂的地震危险性分析方法，中弱活动断裂不能完全用地质学的断裂活动性方法给出其地震危险性，要结合地质学与地震活动性方法综合考虑。

胡聿贤主编，1999. 地震安全性评价技术教程. 北京：地震出版社.

黄玮琼，李文香，曹学锋，1994. 中国大陆地震资料完整性研究之二——分区地震资料基本完整的起始年分布图像. 地震学报，16（4）：423—432.

焦远碧，吴开统，1990. 中国地震台网检测能力及台网观测条件质量评定. 中国地震，6（4）：1—7.

雷建成，时振梁，鄢家全，1991. 确定潜在震源区地震年平均发生率的方法研究. 地震学报，13（4）：496—504.

闻学泽，2007. 隐伏活动断层的现今活动习性与地震危险性评价方法. 四川省地震局、中国地震局地质研究所.

中国地震局，2005. 中国地震活动断层探测技术系统技术规程. 北京：地震出版社.

Connel C.A., 1968. Engineering Seismic Risk Analysis. Bull. Seis. Soc. Am., 58 (5):1583—1606.

Occurrence Probability Evaluation of the Maximum Potential Earthquake along Zhengzhou City Faults

Wang Ji1), Tian Qinjian2)and Gao Zhanwu1)
1) China Earthquake Disaster Prevention Center, Beijing 100029, China
2) Institute of Earthquake Science, China Earthquake Administration，Beijing 100036, China

Based on the results from project “The active fault detection and seismic risk evaluation (Phase II) of Zhengzhou City”, there exist two east-west trending faults, i.e Lao Yachen fault and the Shanjie fault. The Lao Yachen fault was not active in Quaternary, whereas it is likely for a M5-5.5 earthquake to occur along the Shanjie fault. In this study, we divided the study area into statistical units and obtained the seismic parameters and the spatial distribution function for the potential source areas. Finally, we estimated the probability of M5.0 earthquake along faults in Zhengzhou City as 6% in the next 50 years and 11% in 100 years.

王继，田勤俭，高占武，2012．郑州市断裂最大潜在地震发震概率评价．震灾防御技术，7（3）：250—260.

我国地震重点监视防御区活动断层地震危险性评价项目资助

2012-06-27

王继，男，生于1972年。副研究员，一级地震安评工程师。主要研究方向：工程地震。E-mail: 88015314@163.com

Κey words：Zhengzhou City；Fault detection；Seismic risk evaluation