Calculation and Analysis of Seismic Activity Parameters b Value and ν4 in the Seismic Statistical Region−Taking North China as An Example
-
摘要: 本文通过收集最新现代地震资料,采用历史地震与现代地震联合求取地震活动性参数的方法,按不同地震统计时段和震级段,基于Matlab程序中给定函数计算并拟合了华北平原地震统计区、郯庐地震统计区、长江下游-南黄海地震统计区、汾渭地震统计区的b值和ν4,并将结果与五代图推荐的b值和ν4进行地震危险性分析对比,进而分析最新现代地震资料对b值和ν4的影响。研究结果对地震统计区确定b值和ν4有一定参考意义,对地震区划、工程场地地震安全性评价、省级地震危险性区划编制工作具有实际意义。Abstract: Based on the given function in the Matlab program, this paper collects the latest modern seismic data, adopts the idea of combining historical earthquakes and modern earthquakes to obtain seismicity parameters, calculates and fits the b and v4 values of four seismic statistical regions, namely, the North China Plain seismic statistical region, the Tanlu seismic statistical region, the lower Yangtze River South Yellow Sea seismic statistical region, and the Fenwei seismic statistical region, according to different seismic statistical periods and magnitude segments, And compared and calculated the seismic hazard analysis with the recommended b and ν4 values in the compilation of the Five Generation Map, and further analyzed the impact of the latest modern seismic data on the b and ν4 values. The research results of this article have certain reference significance for determining the b and ν4 values in earthquake statistical areas, and have practical significance for earthquake zoning, seismic safety evaluation of engineering sites, and the compilation of provincial-level seismic hazard zoning.
-
Key words:
- Seismic activity parameter /
- b Value /
- ν4 /
- North China region
-
引言
建筑物作为灾害的主要承灾体(陈振拓等,2012),其空间分布是震前开展地震灾害风险评估、震后进行灾害损失快速评估的基础数据(王龙等,2007;谷国梁等,2016;李皓等,2018)。基于精细可靠的建筑物空间分布数据开展评估,震前可给出抗震薄弱环节及相应的抗震措施、危旧房屋集中地段、高风险小区(孙柏涛等,2017;王志涛等,2019);震后可给出更可靠的灾情快速评估结果和可能被困人员的空间分布,为地震应急处置和现场应急救援提供可靠的信息支持,从而达到减轻地震灾害损失的目的(韩贞辉等,2013;丁文秀等,2015)。
目前,建筑物数据主要有两类:一类是建筑物统计数据,其广泛应用于震后快速评估中(王晓青等,2009;徐国栋等,2011),该数据与所属区域的空间特征脱节而无法反映区域内建筑物空间分布特点,仅反映统计单元内建筑物的总体情况,且其可靠性和精细化程度较低,数据更新速度慢,时效性较差;另一类是基于多源数据融合方式生成的建筑物公里格网数据(韩贞辉,2013;杨海霞等,2015;杨海霞,2015),该数据虽能较好地反映建筑物空间分布的宏观特征,但不能准确反映建筑物空间分布的细节特征,且未考虑不同地域建筑物空间分布的差异,无法反映区域间建筑物空间分布的异质性。然而建筑物空间分布具有明显的地域差异,采用同一建筑物空间分布格网化模型具有一定的局限性。
针对未考虑建筑物空间分布异质性问题,在充分考虑目标区域自然环境、社会经济条件等差异的基础上,开展建筑物特征一致性分区。根据分区结果,采用抽样的方法获取各区域建筑物样本数据。在此基础上,探究不同区域建筑物空间分布与各影响因子之间的关系,分区构建基于多因素的建筑物数据空间化模型,并根据模型生成建筑物格网数据,从而提高建筑物数据格网化精度。
1. 研究方法
1.1 总体技术思路
建筑物空间分布是受诸多因素影响的复杂非线性问题,不同地区的社会经济情况、自然地理条件等存在差异,建筑物空间分布特征也不相同(程晓亮等,2008;曾祥贵,2013)。因此,需将具有相似建筑物空间分布特征的区域划分为同一区域。在此基础上,分析不同区域建筑物空间分布与各影响因子之间的关系,分区构建建筑物空间分布格网化模型,反演建筑物空间分布情况。具体技术路线如图1所示。
首先,从自然地理和社会经济两方面出发选择分区指标,并以乡镇为基本单元进行统计,利用主成分分析方法提取主要特征,对其进行特征一致性分区。然后,在归纳总结目前基于多因素的建筑物空间分布格网化方法考虑影响因子的基础上,综合选择自然因子(如高程、坡度、坡向、地形起伏度及河流等)、社会经济因子(如道路、土地利用等),分析各影响因子对建筑物空间分布的影响,并对各因子进行分级设计。在此基础上,针对不同区域,从研究区中选取一定数量的格网作为建模区,采用目视解译的方式获取建模区的建筑物,统计分析不同建模区建筑面积密度与各影响因子之间的关系。采用相关系数法计算不同区域内各影响因子对建筑物空间分布的影响权重,并以归一化的相关系数作为该类因子的权重系数,对于同一类因子的不同子类,采用建筑面积密度作为权重的估计值。在此基础上,分区构建基于多因素的建筑物数据格网化模型。最后,随机选择若干格网作为检验样本,开展模型精度评估。
1.2 建筑物空间分布格网化方法
在借鉴已有研究的基础上,本文依据土地地表覆盖数据分类体系(Gong等,2013;Li等,2017),将研究区域划分为建设用地和非建设用地,其中建设用地指不透水层,包括城镇用地、工矿用地、商服用地、交通设施用地等,非建设用地包括耕地、森林、草地、灌木林、湿地、水域、裸地、冰川等。由于建筑物大多分布在建设用地内,少数或零星分布在非建设用地内(杨海霞等,2015),因此,针对这两类区域分别建立基于多因素的建筑物数据空间化模型。
1.2.1 单因子建筑物空间分布空间化权重
将研究区划分为规则格网,假定共有
$ {N}_{t} $ 个格网,如果格网t内存在影响建筑物空间分布的因子$ {F}_{i} $ ,其与建设用地和非建设用地的空间分布关系如图2所示。为探究因子$ {F}_{i} $ 对建筑物空间分布的影响程度,将其进一步划分为Q个子因子,在假设研究区内子因子$ {F}_{ij} $ 在第m类区域的建筑占地面积密度$ {D}_{mij} $ 不变的基础上,统计格网g内与第m类区域内第j个子因子$ {F}_{ij} $ 的占地面积为$ {A}_{mgij} $ ,则格网g内依据因子$ {F}_{i} $ 估计的建筑面积为:$$ {S}_{ mgi}=\sum _{j=1}^{Q}\left({D}_{mij}\cdot{A}_{mgij}\right) $$ (1) $ {\mathrm{D}}_{\mathrm{m}\mathrm{i}\mathrm{j}} $ 建筑占地面积密度常采用抽样统计的方式获取。假设抽样区域中,第m类区域中第i个因子第j个子因子$ {F}_{ij} $ 的占地面积为$ {A}_{mij} $ 、建筑面积为$ {S}_{mij} $ ,则建筑占地面积密度${{D}}_{{m}{i}{j}}$ 计算公式为:$$ {D}_{mij}={S}_{ mij}/{A}_{mij} $$ (2) 1.2.2 多因子建筑物分布空间化权重
建筑物空间分布是自然因子和社会经济条件共同作用的结果(江东等,2002),因此,基于单一因子的建筑物数据空间化模型可靠性及结果的准确度较低。本文选择多源数据融合的方法,在各类单因子格网建筑面积预测的基础上,根据各类因子之间的权重,通过权重比求和的形式综合同一格网内各类因子建筑面积预测结果。
假设影响建筑物空间分布的因子有M类,在m类区域中,第i类因子
$ {F}_{i} $ 的建筑面积分配权重为$ {W}_{mi} $ ,则格网g内建筑面积的计算公式为:$$ {S}_{ mg}=\sum _{i=1}^{N}\left({W}_{mi}\cdot{S}_{ mgi}\right) $$ (3) 目前,确定各类因子权重
$ {W}_{mi} $ 的方法主要有专家打分法、层次分析法、相关系数法等。其中专家打分法的主观性太强,与专家的偏好程度有很大关系。层次分析法依赖于专家经验,主观性相较于专家打分法有所降低。相关系数法是通过抽样确定的,目前在进行GIS空间分析时使用较多。2. 研究区及数据源
2.1 研究区概况
本文以雅安市为研究区开展试验。该区位于四川省中部,跨越四川盆地和青藏高原两大地形,地势呈北、西、南较高,东部和中部较低。全市辖2个市区和6个县,幅员面积为15 046 km2,其中山地占比94%,平原占比6%。
2.2 数据源与预处理
本文在分析研究区内建筑物空间分布特征及建筑物空间分布影响因子的基础上,模拟建筑物空间分布情况。因此,本文所需的基础数据主要包括行政区划数据、土地覆盖类型数据(FROM-GLC30)、DEM数据、道路交通数据、河流水系数据及人口统计数据。
由于建筑物空间分布具有显著的地域差异,若采用同一建筑物空间分布格网化模型具有一定的局限性,数据准确性难以保证。故本文在构建建筑物格网化模型时,首先进行建筑物空间分布特征一致性分区。建筑物空间分布受诸多因素影响,主要包括自然环境(如高程、坡度、坡向、水系等)及社会条件(如土地利用、道路、人口等),本文综合考虑上述方面进行建筑物一致性分区。
根据研究区建筑物空间分布特征,本文选择人口密度、路网密度、河网密度、平均高程、平均坡度、建设用地面积比例、耕地面积比例共7个指标进行建筑物特征一致性分区。首先,以乡镇为基本单元统计各指标的值,对其进行极差标准化处理,利用主成分分析方法进行降维。分析得到的第一主成分特征值和第二主成分特征值分别为3.743、1.483,且这两个主分量承载了74.68%的原始信息,其余主成分特征值均小于1,表明引入该主成分的影响程度不如1个基本变量。因此,选择第一主成分和第二主成分基本可表达建筑物分布特征。根据主成分分析得到的第一主成分和第二主成分计算公式,计算雅安市143个乡镇主第一主成分和第二主成分得分,再根据各主成分贡献率的比重作为权值,得到各乡镇综合得分,利用ArcGIS软件自然断点法将其划分为5个特征分区,如图3所示。
3. 模型构建与结果分析
3.1 建筑物空间分布影响因子分析
在归纳总结目前基于多因素的建筑物空间分布格网化方法考虑因素的基础上,结合研究区建筑物空间分布特征,综合选取土地利用、高程、坡度、坡向、河流距离、道路距离、地形起伏度7类因子,并对各类因子进行合理分类或分级,具体方案如表1所示。其中,高程等级划分原则为对建筑物分布密集段高程细分,对建筑物分布稀疏段高程粗分;坡度等级划分原则依据国际地理学联合会地貌调查与地貌制图委员制定的分类标准;坡向根据方位划分为阴坡、半阴坡、阳坡、半阳坡。河流距离和道路距离采用等间距划分等级;地形起伏度根据我国基本地貌类型划分等级。
表 1 影响因子子类分级Table 1. Classification of influencing factors土地利用 高程/m 坡度/(°) 坡向/(°) 河流距离/m 道路距离/m 地形起伏度/m 耕地 500~800 平原
(0~0.5)平缓坡
(−1)0~200 0~200 平原(<30) 森林 800~1 100 微斜坡(0.5~2) 向阳坡(135~225) 200~400 200~400 台地
(30~70)草地 1 100~1 400 缓斜坡(2~5) 向阳坡
(45~135,
225~315)400~600 400~600 丘陵
(70~200)灌木林 1 400~1 700 斜坡
(5~15)阴坡(0~45,
315~360)600~800 600~800 小起伏山地(200~500) 湿地 1 700~2 000 陡坡
(15~35)— 800~1 000 800~1 000 中起伏山地(500~1 000) 水体 2 000~2 400 峭坡
(35~55)— 1 000~1 200 1 000~1 200 大起伏山地(1 000~2 500) 不透水面 2 400~2 800 垂直壁
(55~90)— 1 200~1 400 1 200~1 400 极大起伏山地(>2 500) 裸地 2 800~3 200 — — 1 400~1 600 1 400~1 600 — 冰川 3 200~3 800 — — 1 600~1 800 1 600~1 800 — — 3 800~4 400 — — 1 800~2 000 1 800~2 000 — — >4 400 — — >2 000 >2 000 — 3.2 因子权重系数确定
以雅安市为研究区,将其划分为300 m×300 m的规则格网,基于建筑物特征一致性分区结果,从不同区域分别选取一定数量的格网作为建模区,如图4所示。根据天地图影像,采用目视解译的方式获取建模区内单体建筑物数据。
根据研究区各类因子数据(图5),针对不同区域,分别将抽样格网的影响因子与建设用地、非建设用地进行叠加分析,统计区域内各影响因子的子类或分级占地面积和建筑面积,根据式(2)确定各子类或分级的建筑占地面积密度。
针对各类影响因子权重的确定,本文选用相关系数法,根据同一特征分区内某一因子的子类或分级建筑面积密度与子类或分级占地面积比例之间的关系,确定相关系数,并以归一化的相关系数作为该类因子的权重系数,结果如表2所示。
表 2 雅安市不同特征分区抽样统计的各类因子权重Table 2. Weights of various factors in sampling statistics in different regions of Ya'an city分区 区域划分 权重 土地利用 高程 坡度 坡向 河流距离 道路距离 地形起伏度 一区 建设用地 — — 0.263 0.108 — 0.300 0.329 非建设用地 0.244 — 0.225 0.164 0.064 0.083 0.220 二区 建设用地 — — 0.230 0.231 0.144 0.159 0.236 非建设用地 0.163 — 0.178 0.154 0.187 0.141 0.177 三区 建设用地 — 0.216 0.213 0.156 — 0.217 0.198 非建设用地 0.227 0.177 0.224 — — 0.183 0.189 四区 建设用地 — 0.089 0.195 0.185 0.169 0.148 0.214 非建设用地 0.236 0.183 0.221 — — 0.175 0.185 五区 建设用地 — 0.176 0.183 0.169 0.156 0.147 0.170 非建设用地 0.156 0.168 0.192 0.180 — 0.148 0.157 3.3 建筑物空间分布格网化
在确定不同特征区各因子的不同子类或分级的建筑面积密度和各类因子权重的基础上,对雅安市内所有格网进行影响因子、建设用地和非建设用地统计,在此基础上,分别进行建设用地和非建设用地建筑物空间分布预测。考虑到1个格网内有可能出现预测面积小于1栋房屋的建筑面积,这与实际情况不符,需对这部分数据予以舍去。通过对雅安市建筑物开展野外调查,确定建筑用地单栋房屋建筑面积约为150 m2、非建筑用地单栋房屋建筑面积约为200 m2,因此,将建设用地预测面积<150 m2的格网建筑面积舍弃,将非建设用地预测面积<200 m2的格网建筑面积舍弃,即将格网的建筑面积设置为0。最后,将同一格网内的建设用地和非建设用地建筑面积相加,得到雅安市建筑面积预测结果,如图6所示。
3.4 精度验证
为验证建筑物空间分布格网化方法的预测精度,从研究区中随机选择若干格网作为验证样本,对验证格网内的建筑物进行解译,并统计格网内建筑面积,与模型得到的预测值进行对比,从而计算建筑面积相对误差,分级统计结果如表3所示。由表3可知,严重低估或严重高估的比例≤10%,较准确估计的比例接近75%,说明本文模型具有较高的精度。精度较差的格网主要分布在建筑分布特别密集或稀疏的地区,其原因可能在于建模时这部分样本抽取的比例较小,构建的模型对这类格网估计的结果不是很准确或未考虑因素之间的相关性,造成信息冗余。
表 3 相对误差分级统计Table 3. Statistics of relative error classification分级 数目 比例/% 严重低估,<−50% 794 4.7 一般低估,[−50%,−20%) 602 3.6 较准确估计,[−20%,20%] 12 568 74.7 一般高估,(20%,50%] 1 228 7.3 严重高估,>50% 1 626 9.7 4. 结论与讨论
本文考虑到建筑物空间分布的地域差异,根据影响建筑物空间分布的因素进行建筑物特征一致性分区。在分析建筑物空间分布影响因子的基础上,结合区域建筑物空间分布特征,综合选取土地利用、高程、坡度、坡向、河流距离、道路距离、地形起伏度7类因子,以归一化的相关系数作为各因子的权重系数,将建筑面积密度作为各因子子类或分级的权重,分区构建基于多因素的建筑物格网化模型,预测建筑物空间分布情况。以雅安市为研究区,得到了300 m格网的建筑物空间分布数据,并开展模型精度检验,结果表明该方法具有较好的精度。基于本文提出的方法获取的建筑物空间分布格网可为震前防御、震中救灾和震后重建提供可靠的数据支撑。
建筑物空间分布情况是考虑诸多影响因素的复杂问题,影响因子的选取是否全面、因子之间的相关性如何确定均会对建筑物空间分布产生影响。未来可深入分析影响建筑物空间分布的因素,进一步提高模型的精度。
-
表 1 C-S法余震时间
Table 1. Aftershock time in C-S method
主震震级/级 时间/天 主震震级/级 时间/天 3.0~3.5 11.5 6.1~6.5 510 3.6~4.0 22 6.6~7.0 790 4.1~4.5 42 7.0~7.5 915 4.6~5.0 83 7.6~8.0 960 5.1~5.5 155 8.0~8.5 985 5.6~6.0 290 ≥8.6 985 表 2 地震统计区2011年以来4级以上地震发生次数
Table 2. Total number of earthquakes with a magnitude of 4 or above in different earthquake statistical regions since 2011
统计区名称 地震数/次 4.0~4.4 4.5~4.9 5.0~5.4 5.0~5.9 华北平原 15 10 1 1 郯庐 30 12 1 — 长江下游-南黄海 15 11 3 — 汾渭 15 4 — — 表 3 地震统计区各震级档地震年发生率一览表(地震资料截止至2010年12月)
Table 3. List of annual occurrence rates of earthquakes with different magnitudes in earthquake statistical regions(as of December 2010)
地震统计区 时间范围/年 各震级档地震年发生率 4.0~4.4 4.5~4.9 5.0~5.4 5.5~5.9 6.0~6.4 6.5~6.9 7.0~7.4 7.5~7.9 8.0~8.4 8.5~8.9 华北平原地震统计区 1484 —2010— — 0.1763 0.0967 0.0398 0.0228 0.0096 0.0057 0.0019 — 1500 —2010— — 0.1779 0.0996 0.0411 0.0235 0.0097 0.0059 0.0020 — 1791 —2010— — 0.2222 0.0998 0.0500 0.0318 0.0184 0.0091 — — 1950—2010 3.6998 1.5711 0.3443 0.1176 0.0849 0.0685 0.0357 0.0164 — — 郯庐地震统计区 1477 —2010— — 0.1217 0.0693 0.0356 0.0150 0.0112 0.0037 0.0019 0.0019 1500 —2010— — 0.1213 0.0685 0.0352 0.0157 0.0117 0.0039 0.0020 0.0020 1892—2010 — — 0.1849 0.0924 0.0672 0.0336 0.0168 — — — 1970—2010 3.0732 1.0244 0.4634 0.0732 0.0244 0.0244 0.0244 — — — 长江下游-南黄海地震统计区 1485 —2010— — 0.1065 0.0532 0.0342 0.0171 0.0019 — — — 1491 —2010— — 0.1077 0.0538 0.0346 0.0173 0.0019 — — — 1839—2010 — — 0.2558 0.1395 0.0872 0.0465 0.0058 — — — 1970—2010 2.3415 0.9512 0.3659 0.1220 0.0976 — — — — — 汾渭地震统计区 1000 —2010— — — — 0.0247 0.0208 0.0069 0.0030 0.0020 — 1209 —2010— — — — 0.0274 0.0224 0.0075 0.0037 0.0025 — 1484 —2010— — 0.1803 0.0854 0.0323 0.0247 0.0095 0.0038 0.0019 — 1500 —2010— — 0.1761 0.0822 0.0294 0.0235 0.0098 0.0039 0.0020 — 1950-2010 2.0984 1.0164 0.4426 0.1311 0.0164 — — — — — 表 4 地震统计区各震级档地震年发生率一览表(地震资料截止至2023年12月)
Table 4. List of annual occurrence rates of earthquakes with different magnitudes in earthquake statistical regions(as of December 2023)
地震统计区名称 时间范围/年 各震级档地震年发生率 4.0~4.4 4.5~4.9 5.0~5.4 5.5~5.9 6.0~6.4 6.5~6.9 7.0~7.4 7.5~7.9 8.0~8.4 8.5~8.9 华北平原地震统计区 1484 —2023— — 0.1759 0.0963 0.0389 0.0222 0.0093 0.0056 0.0019 — 1500 —2023— — 0.1775 0.0992 0.0401 0.0229 0.0095 0.0057 0.0019 — 1791 —2023— — 0.2189 0.0987 0.0472 0.0300 0.0172 0.0086 — — 1950—2023 3.5405 1.5000 0.3108 0.1081 0.0676 0.0541 0.0270 0.0135 — — 郯庐地震统计区 1477 —2023— — 0.1207 0.0676 0.0347 0.0146 0.0110 0.0037 0.0018 0.0018 1500 —2023— — 0.1202 0.0668 0.0344 0.0153 0.0115 0.0038 0.0019 0.0019 1892—2023 — — 0.1742 0.0833 0.0606 0.0303 0.0152 — — — 1970—2023 3.1296 1.0185 0.3704 0.0556 0.0185 0.0185 0.0185 — — — 长江下游-南黄海地震统计区 1485 —2023— — 0.1095 0.0520 0.0334 0.0167 0.0019 — — — 1491 —2023— — 0.1107 0.0525 0.0338 0.0169 0.0019 — — — 1839—2023 — — 0.2432 0.1297 0.0811 0.0432 0.0054 — — — 1970—2023 2.3148 0.9815 0.3333 0.0926 0.0741 — — — — — 汾渭地震统计区 1000 —2023— — — — 0.0244 0.0205 0.0068 0.0029 0.0020 — 1209 —2023— — — — 0.0270 0.0221 0.0074 0.0037 0.0025 — 1484 —2023— — 0.1759 0.0833 0.0315 0.0241 0.0093 0.0037 0.0019 — 1500 —2023— — 0.1718 0.0802 0.0286 0.0229 0.0095 0.0038 0.0019 — 1950—2023 1.9865 0.8919 0.3649 0.1081 0.0135 — — — — — 表 5 各地震统计区地震活动性参数对比一览表
Table 5. Comparison of seismic activity parameters among different seismic statistical regions
地震统计区名称 五代图推荐 本文结果 对比分析 b ν4 b ν4 华北平原地震统计区 0.86 4.6 0.81 3.9 b值减小、ν4减小 郯庐地震统计区 0.85 4.0 0.80 3.0 b值减小、ν4减小 长江下游-南黄海地震统计区 0.85 3.0 0.84 2.8 b值减小、ν4减小 汾渭地震统计区 0.78 2.5 0.77 2.3 b值减小、ν4减小 表 6 2套地震活动性参数方案下各地震统计区不同超越概率基岩地震动峰值加速度表
Table 6. Peak acceleration of bedrock seismic motion with different exceedance probabilities in each seismic statistical region under two sets of seismic activity parameter schemes
地震统计区名称 计算控制点 加速度分区 方案 50 a(63%) 50 a(10%) 50 a(2%) 峰值 差值 峰值 差值 峰值 差值 华北平原地震统计区 hb1
(37.2410 °N,115.6767 °E)0.05 g A 20.5 −0.7 63.7 −0.2 115.5 0.7 B 19.8 63.5 116.2 hb2
(36.1002 °N,116.3139 °E)0.10 g A 25.7 −0.7 84.6 −0.6 150.2 −0.2 B 25.0 84.0 150.0 hb3
(35.4580 °N,115.8361 °E)0.15 g A 37.4 −1.4 131.8 −0.5 245.1 −0.3 B 36.0 131.3 245.4 hb4
(36.1227 °N,115.6058 °E)0.20 g A 41.0 −1.5 188.3 2.3 420.0 5.1 B 39.5 190.6 425.1 hb5
(39.4218 °N,117.9615 °E)0.30 g A 54.8 −0.9 240.3 5.3 528.8 10.1 B 53.9 245.6 538.9 郯庐地震统计区 tl1
(35.2231 °N,116.8219 °E)0.05 g A 21.2 −1.6 67.0 −3.5 128.3 −4.7 B 19.6 63.5 123.6 tl2
(35.8315 °N,117.5840 °E)0.10 g A 28.3 −2.9 99.9 −5.6 201.7 −7.9 B 25.4 94.3 193.8 tl3
(36.1563 °N,119.4777 °E)0.15 g A 25.3 −2.4 102.8 −3.8 208.8 −3.0 B 22.9 99.0 205.8 tl4
(35.8153 °N,118.9612 °E)0.20 g A 32.1 −3.3 148.7 −4.2 355.7 −1.7 B 28.8 144.5 354.0 tl5
(35.0825 °N,118.6671 °E)0.30 g A 33.7 −3.5 177.8 −2.0 460.9 6.8 B 30.2 175.8 467.7 长江下游-南黄海地震统计区 nh1
(33.5189 °N,119.5458 °E)0.05 g A 18.4 −0.4 57.1 −0.7 102.1 −1.0 B 18.0 56.4 101.1 nh2
(32.9792 °N,120.1834 °E)0.10 g A 26.8 −0.7 98.2 −1.3 190.8 −2.0 B 26.1 96.9 188.8 nh3
(32.8316 °N,120.7629 °E)0.15 g A 32.9 −1.0 123.6 −1.5 235.5 −1.7 B 31.9 122.1 233.8 nh4
(33.2425 °N,120.7762 °E)0.20 g A 34.0 −1.0 183.0 −1.5 413.1 −1.6 B 33.0 181.5 411.5 汾渭地震统计区 fw1
(36.8145 °N,113.2341 °E)0.05 g A 21.4 −0.5 66.0 −1.1 117.4 −1.6 B 20.9 64.9 115.8 fw2
(37.2475 °N,111.4661 °E)0.10 g A 27.8 −1.1 106.4 −2.1 218.9 −2.8 B 26.7 104.3 216.1 fw3
(35.5529 °N,111.0402 °E)0.15 g A 35.9 −1.3 147.3 −2.8 304.4 −3.9 B 34.6 144.5 300.5 fw4
(38.3544 °N,112.9154 °E)0.20 g A 36.0 −1.5 169.6 −3.1 373.7 −3.9 B 34.5 166.5 369.8 fw5
(36.3985 °N,111.7794 °E)0.30 g A 43.2 -1.8 216.7 −3.6 493.4 −3.5 B 41.4 213.1 489.9 注:表中A代表五代图地震活动性参数方案;B代表本文地震活动性参数方案,表中地震动峰值加速度单位为Gal。差值为方案B与方案A计算得到的基岩地震动峰值加速度之差。 -
陈凌,刘杰,陈颙等,1998. 地震活动性分析中余震的删除. 地球物理学报,41(S1):244−252.Chen L., Liu J., Chen Y., et al., 1998. Afiershock deletion in seismicity analysis. Chinese Journal of Geophysics, 41(S1): 244−252. (in Chinese) 陈阳,吕悦军,谢卓娟等,2013. 地震活动性参数 b 值的研究. 地壳构造与地壳应力文集,38−47.Chen Y., Liu Y. J., Xie Z. J., et al., 2013. Review of the study of seismicity parameter B-value. Bulletin of the Institute of Crustal Dynamics, 38−47. (in Chinese) 高孟潭,2015. GB 18306−2015《中国地震动参数区划图》宣贯教材. 北京:中国标准出版社. 黄玮琼,时振梁,曹学锋,1989. b 值统计中的影响因素及危险性分析中 b 值的选取. 地震学报,11(4):351−361.Huang W. Q., Shi Z. L., Cao X. F., 1989. Factors influencing the estimation of b value and the selection of b Value in hazard analysis. Acta Seismologica Sinica, 11(4): 351−361. (in Chinese) 孟昭彤,刘静伟,谢卓娟等,2021. b 值的时空分布特征与地震危险性的关联分析. 地球物理学进展,36(1):30−38. doi: 10.6038/pg2021EE0025Meng Z. T., Liu J. W., Xie Z. J., et al., 2021. Analysis of the correlation between the temporal-spatial distribution of b -value and seismic hazard: a review. Progress in Geophysics, 36(1): 30−38. (in Chinese) doi: 10.6038/pg2021EE0025 潘华,李金臣,2006. 地震统计区地震活动性参数 b 值及 ν 4不确定性研究. 震灾防御技术,1(3):218−224. doi: 10.3969/j.issn.1673-5722.2006.03.006Pan H., Li J. C., 2006. Study on uncertainties of seismicity parameters b and v 4 in seismic statistical zones. Technology for Earthquake Disaster Prevention, 1(3): 218−224. (in Chinese) doi: 10.3969/j.issn.1673-5722.2006.03.006 潘华,高孟潭,谢富仁,2013. 新版地震区划图地震活动性模型与参数确定. 震灾防御技术,8(1):11−23. doi: 10.3969/j.issn.1673-5722.2013.01.002Pan H., Gao M. T., Xie F. R., 2013. The earthquake activity model and seismicity parameters in the New Seismic Hazard Map of China. Technology for Earthquake Disaster Prevention, 8(1): 11−23. (in Chinese) doi: 10.3969/j.issn.1673-5722.2013.01.002 任雪梅,高孟潭,冯静,2011. 地震目录的完整性对 b 值计算的影响. 震灾防御技术,6(3):257−268. doi: 10.3969/j.issn.1673-5722.2011.03.005Ren X. M., Gao M. T., Feng J., 2011. Effect of completeness of earthquake catalogue on calculating b value. Technology for Earthquake Disaster Prevention, 6(3): 257−268. (in Chinese) doi: 10.3969/j.issn.1673-5722.2011.03.005 汪素云,俞言祥,2009. 震级转换关系及其对地震活动性参数的影响研究. 震灾防御技术,4(2):141−149. doi: 10.3969/j.issn.1673-5722.2009.02.002Wang S. Y., Yu Y. X., 2009. Research on empirical relationship of earthquake magnitude scales and its influence on seismicity parameters. Technology for Earthquake Disaster Prevention, 4(2): 141−149. (in Chinese) doi: 10.3969/j.issn.1673-5722.2009.02.002 吴果,周庆,冉洪流,2019. 震级-频度关系中 b 值的极大似然法估计及其影响因素分析. 地震地质,41(1):21−43. doi: 10.3969/j.issn.0253-4967.2019.01.002Wu G., Zhou Q., Ran H. L., 2019. The maximum likelihood estimation of b -value in magnitude-frequency relation and analysis of its influencing factors. Seismology and Geology, 41(1): 21−43. (in Chinese) doi: 10.3969/j.issn.0253-4967.2019.01.002 谢卓娟,吕悦军,张力方等,2012. 基于现代地震资料确定汾渭地震带分区及其地震活动性参数. 地球物理学进展,27(3):894−902. doi: 10.6038/j.issn.1004-2903.2012.03.009Xie Z. J., Lv Y. J., Zhang L. F., et al., 2012. Subarea and seismicity parameters of the Fen-Wei seismic zone based on the modern seismic data. Progress in Geophysics, 27(3): 894−902. (in Chinese) doi: 10.6038/j.issn.1004-2903.2012.03.009 谢卓娟,吕悦军,兰景岩等,2013. b 值和 V 4的统计分析及其不确定性对地震危险性分析结果的影响研究. 地震研究,36(1):86−92. doi: 10.3969/j.issn.1000-0666.2013.01.013Xie Z. J., Lv Y. J., Lan J. Y., et al., 2013. Research on the effects of statistical analysis of b -value and V 4 and its uncertainty on seismic risk analysis result. Journal of Seismological Research, 36(1): 86−92. (in Chinese) doi: 10.3969/j.issn.1000-0666.2013.01.013 徐果明,周蕙兰,1982. 地震学原理:421. 北京:科学出版社. 鄢家全,韩炜,高孟潭,1996. 地震活动性参数的不确定性及其对区划结果的影响. 中国地震,12(S1):71−77.Yan J. Q., Han W., Gao M. T., 1996. Uncertainty of seismicity parameters and its affects to the seismic zoning. Earthquake Research in China, 12(S1): 71−77. (in Chinese) 张效亮,2017. 基于Matlab的地震带活动性参数计算. 山西建筑,43(22):49−50. doi: 10.3969/j.issn.1009-6825.2017.22.027Zhang X. L., 2017. Seismic belt activity parameter calculation on the basis of Matlab. Shanxi Architecture, 43(22): 49−50. (in Chinese) doi: 10.3969/j.issn.1009-6825.2017.22.027 Console R., Gasparini C., De Simoni B., et al., 1979. Preambolo al Catalogo Sismico Nazionale (CSN). I criteri di informazione del Catalogo Sismico Nazionale (CSN). Annals of Geophysics, 32(1): 37−77. Keilis-Borok V. I., Knopoff L., Rotvain I. M., 1980. Bursts of aftershocks, long-term precursors of strong earthquakes. Nature, 283(5744): 259−263. doi: 10.1038/283259a0 -