• ISSN 1673-5722
  • CN 11-5429/P

地震统计区地震活动性参数b值和ν4计算分析−以华北地区为例

沈得秀 王庆民 许洪泰 葛孚刚 蒋其峰

李敏, 杨立国, 陈海鹏, 盛特奇. 杭州市典型土层剪切波速与埋深间的关系分析[J]. 震灾防御技术, 2020, 15(1): 77-88. doi: 10.11899/zzfy20200108
引用本文: 沈得秀,王庆民,许洪泰,葛孚刚,蒋其峰,2024. 地震统计区地震活动性参数b值和ν4计算分析−以华北地区为例. 震灾防御技术,19(2):241−252. doi:10.11899/zzfy20240204. doi: 10.11899/zzfy20240204
Li Min, Yang Liguo, Chen Haipeng, Sheng Teqi. Relationship between Shear Wave Velocity and Soil Depth of Typical Soil Layers in Hangzhou Area[J]. Technology for Earthquake Disaster Prevention, 2020, 15(1): 77-88. doi: 10.11899/zzfy20200108
Citation: Shen Dexiu, Wang Qingmin, Xu Hongtai, Ge Fugang, Jiang Qifeng. Calculation and Analysis of Seismic Activity Parameters b Value and ν4 in the Seismic Statistical Region−Taking North China as An Example[J]. Technology for Earthquake Disaster Prevention, 2024, 19(2): 241-252. doi: 10.11899/zzfy20240204

地震统计区地震活动性参数b值和ν4计算分析−以华北地区为例

doi: 10.11899/zzfy20240204
基金项目: 山东省地震局科技创新团队培育专项(TD202303);山东省地震局一般科研项目(YB2404)
详细信息
    作者简介:

    沈得秀,女,生于1975年。硕士,高级工程师。主要从事地震活动性和工程地震研究工作。E-mail:shendexiu@163.com

    通讯作者:

    王庆民,男,生于1974年。硕士,工程师。主要从事地震预报等方面研究。E-mail:qingminwang@163.com

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有一定参考意义,对地震区划、工程场地地震安全性评价、省级地震危险性区划编制工作具有实际意义。
  • 剪切波速是表征场地土特性的重要参数,对于场地类别划分、砂土液化、软土震陷判别等至关重要。在实际工程建设中,受剪切波速测试环境、场地范围、经济等因素的影响,某些岩土中出现剪切波速数值缺失的情况,此时工程者往往采用经验统计模型进行剪切波速数值的补充。

    目前学者已研究众多常规土类剪切波速与埋深间关系的经验统计模型,如王琦等(2018)基于天津地区地震安全性评价的实测剪切波速资料,采用幂函数模型,定性探讨了常见土类剪切波速与埋深的相关性;荣棉水等(2017)利用非线性最小二乘法,采用指数函数、一次函数、二次函数、幂函数、幂函数+常数函数、幂函数+一次函数6种回归模型对渤海海域7类土剪切波速和埋深的关系进行回归分析,指出最优模型为幂函数+一次函数模型;董菲蕃等(2013)基于福建沿海福州、泉州、漳州3个盆地大量的剪切波速实测资料,利用线性和二次多项式拟合得到各盆地5类土体等效剪切波速的空间分布特征;刘红帅等(2010)基于全国区域范围的实测钻孔资料,利用最小二乘法分别采用线性模型、指数模型和一元二次多项式模型建立分场地类别和不分场地类别下的常规土类(不包括淤泥质土)碎石土、砂土、粉土和黏性土剪切波速与埋深间的统计公式;战吉艳等(2009)利用苏州地区实测钻孔剪切波速资料进行研究,结果表明对于土层深度不大于50m的粘土和不小于40m的粉土,可采用幂函数形式模型的经验关系预测场地土剪切波速与土层深度的关系,对于粉质粘土和土层深度不小于40m的粉砂,可采用线性函数、幂函数分段形式模型的经验关系进行预测。经验公式受土层区域特性的影响显著,如乔峰等(2019)利用北京地区实测剪切波速资料建立的5种常见土类剪切波速与埋深之间关系的推荐模型,与利用全国数据建立的模型进行对比分析,指出区域性对剪切波速与深度之间的相关性存在影响;王琦等(2018)也研究指出剪切波速与埋深的经验公式受土层区域特性的影响显著。

    由此可知,受地区土层差异性或模型关系影响,未统一给出相关的经验公式(荣棉水等,2017)。因此,选取合适的统计模型并建立合理的经验公式成为关键。

    影响岩土剪切波速的主要因素为岩土类型、成因、物理状态和埋深(程祖锋等,1997)。很多学者也进行了相关性研究,如闫振军等(2019)通过华北地区10个城市928个钻孔的剪切波速从岩性条件、沉积环境等方面分析探讨了华北地区剪切波速随深度变化的特征;薛锋等(2014)给出考虑西安市地貌单元的细砂、中砂、粗砂、黄土、粉质粘土、古土壤、卵石等土类剪切波速与土层深度统计经验关系式;刘红帅等(2010)基于全国范围的实测钻孔资料,给出区分场地类别和不区分场地类别的情况下常规土类(不包括淤泥质土)碎石土、砂土、粉土和黏性土剪切波速与埋深之间的统计公式。

    近年来,杭州作为长三角区域经济发展的重要城市,在城市建设中开展了大量地震安全性评价工作,积累了丰富的剪切波速资料。基于以上研究背景,考虑场地未明确土体物理状态的不利条件、建筑物场地类别在一定程度上反映了土的成因和物理状态以及国内学者在此方面的研究,本文利用现行3种最常规的剪切波速与埋深关系的经验模型进行分析,给出杭州市粉质粘土、淤泥质粉质粘土、粉砂、砂质粉土和粘质粉土5类典型土层在考虑场地分类和未考虑场地分类下的最优模型和相关系数,以供工程和研究参考。

    浙江省地貌单元可划分为浙西中山丘陵区、浙中盆地区、浙北平原区、浙东低山丘陵区、浙南中山区和浙东南沿海丘陵平原与岛屿区,杭州市位于浙西北地区,其西部(如临安和桐庐地区)属于浙西中山丘陵区;西南、西北山区属于中生代断块隆起;东部属于浙北平原区,地势低平,河流密集。杭州市工程地质单元主要以第四系地层为主,因新构造运动的影响,第四纪沉积成因类型繁多,沉积物的多次交替堆积和剥蚀,使得杭州地区覆盖层厚度变化较大(丁伯阳等,2012)。工程地质层主要包括前第四系和更新系统以来的第四系,前第四系(AnQ)主要由石灰岩、泥岩、砂岩等表征的基岩土层构成,第四系主要由中更新统、上更新统、全新统组成,成因类型主要有冲积相、冲湖积、湖沼积、冲海积、冲湖积、海积、冲积等,土类主要以稍密的砂质粉土、流塑的淤泥质土、硬塑或软塑的粘土和中密的粉细砂为主(浙江省住房和城乡建设厅,2019)。

    本文收集了杭州市2011—2016年以来开展的地震安全性评价工作中98个项目报告的296个钻孔资料,钻孔分布示意如图 1所示。场地类别以Ⅱ、Ⅲ、Ⅳ类为主,场地分类方法参照《建筑抗震设计规范》(GB 50011—2010)(中华人民共和国住房和城乡建设部,2016)公式4.1.5和表4.1.6确定,场地土类型主要为中硬土、中软土和软弱土,中硬土场地主要分布于临安、富阳地区,覆盖层厚度以4—20m为主;中软土场地分布于主城区,覆盖层厚度约为30—60m,个别钻孔深度达100m;软弱土场地主要分布于钱塘江附近。因Ⅳ类场地钻孔数较少,本文未展开研究。本次场地钻孔土层剪切波速测试均采用单孔法,波速测试间距为1m。钻孔土类主要以粉质粘土、淤泥质粉质粘土、粉砂、砂质粉土和粘质粉土为主,因此,汇总上述5类土的土层埋深点与对应剪切波速值,分别研究场地未分类和分类(Ⅱ、Ⅲ)下的剪切波速随埋深变化的关系,不同土类剪切波速分布范围如表 1所示。

    图 1  杭州市地区场地钻孔位置分布图
    Figure 1.  Location of boreholes in Hangzhou
    表 1  研究区不同土类剪切波速的分布范围
    Table 1.  Distribution ranges of shear wave velocity of different soils in study area
    岩性 粉质粘土 淤泥质粉质粘土 粉砂 砂质粉土 粘质粉土
    场地类别 总计 总计 总计 总计 总计
    测试点数 1915 324 1591 902 31 871 383 17 366 904 88 816 330 69 261
    最小值/m·s-1 115 115 117 98 102 98 134 134 139 112 122 112 112 116 112
    最大值/m·s-1 390 348 390 245 215 245 357 303 357 339 310 339 321 321 270
    平均值/m·s-1 249.4 203.4 258.8 153.9 139.5 154.4 227.5 255.0 226.2 163.1 202.4 158.9 154.2 165.0 151.3
    标准差/m·s-1 59.0 62.9 53.6 29.1 30.9 28.9 54.5 42.7 54.7 37.0 56.5 31.5 32.8 45.8 27.8
    测试深度/m 2—78 2—44 2—78 2—46 4—24 2—46 1—56 1—30 5—56 1—57 1—34 1—57 2—44 2—36 2—44
    下载: 导出CSV 
    | 显示表格

    众多研究结果表明,剪切波速与埋深的关系可用线性函数、多项式函数、幂函数、幂函数+常函数等众多模型表示,本文选取常见的3种经验模型进行分析,分别为线性函数模型Vs=a+bH、二次函数模型Vs=a+bH+cH2、常函数+幂函数模型Vs=a+bHc,其中Vs为土体剪切波速,abc分别为拟合系数。通过Origin统计分析功能,分别给出5类土的回归曲线及相关系数,最终采用拟合优度R2选出适合本地区的拟合模型,拟合优度为0<R2<1,R2越接近1时,拟合效果越佳。

    为更好地展示剪切波速与埋深的关系,首先绘制296个钻孔(不包括Ⅳ类场地)中粉质粘土、淤泥质粉质粘土、粉砂、砂质粉土和粘质粉土波速与埋深散点图,发现深度相近时剪切波速相差较大,总结原因为:剪切波速测试受环境和仪器精度的影响,出现个别坏点;因土体状态导致剪切波速相差较大,即受软硬夹层的影响。考虑回归分析仅能反映一般土质条件下的统计关系(董菲蕃等,2013陈国兴等,1998),因此本文首先对于个别离散趋势较大的点和软硬夹层的剪切波速,参照董菲蕃等(2013)陈国兴等(1998)的方法进行相关性处理,处理完成的数据如图 2所示。由图 2可知,总体上,剪切波速与埋深具有显著的相关性,即剪切波速随埋深具有明显的增长趋势。采用3种统计关系模型分别对不同场地分类情况下的5类土进行回归分析,拟合趋势如图 37所示,拟合优度如表 2所示,分布情况如图 8所示,为表达简洁,分别给出各土类在不同模型下的拟合工况编号,具体如表 2所示。

    图 2  296个钻孔资料中5类常规土剪切波速与埋深散点图
    Figure 2.  Scatter diagram of shear wave velocity and buried depth of five types of conventional soils in 296 borehole data
    图 3  不同场地分类情况下粉质粘土剪切波速与埋深散点图
    Figure 3.  Scatter diagrams of shear wave velocity and buried depth of silty clay under different site classifications
    图 4  不同场地分类情况下淤泥质粉质粘土剪切波速与埋深散点图
    Figure 4.  Scatter diagram of shear wave velocity and buried depth of mucky silty clay under different site classifications
    图 5  不同场地分类情况下粉砂剪切波速与埋深散点图
    Figure 5.  Scatter diagram of shear wave velocity and buried depth of silty sand under different site classifications
    图 6  不同场地分类情况下砂质粉土剪切波速与埋深散点图
    Figure 6.  Scatter diagram of shear wave velocity and buried depth of sandy silt under different site classifications
    图 7  不同场地分类情况下粘质粉土剪切波速与埋深散点图
    Figure 7.  Scatter diagram of shear wave velocity and buried depth of clayey silt under different site classifications
    表 2  不同场地分类情况下5类土拟合优度
    Table 2.  Goodness of fitting of five types of soil under different site classification
    岩土类型 工况 R2
    场地未分类 Ⅱ类场地 Ⅲ类场地
    粉质粘土 FN1 0.86870 0.95298 0.86330
    FN2 0.88890 0.95301 0.87860
    FN3 0.89140 0.95284 0.88130
    淤泥质粉质粘土 YN1 0.86130 0.87627 0.86470
    YN2 0.86560 0.90946 0.86850
    YN3 0.86460 0.90960 0.86770
    粉砂 FS1 0.83000 0.95942 0.83480
    FS2 0.84820 0.95659 0.84980
    FS3 0.84140 0.95670 0.84540
    砂质粉土 SF1 0.88030 0.97654 0.91490
    SF2 0.88110 0.97672 0.91530
    SF3 0.88020 0.97725 0.91550
    粘质粉土 NF1 0.83070 0.95183 0.92820
    NF2 0.83860 0.95568 0.92910
    NF3 0.83520 0.95554 0.92810
    下载: 导出CSV 
    | 显示表格
    图 8  不同场地分类情况下5类土拟合优度
    Figure 8.  Goodness of fitting of five types of soils under different site classification

    图 8可知,总体上,Ⅱ类场地拟合优度均大于场地未分类和Ⅲ类场地,由图 27可知,Ⅱ类场地下各类土拟合数据样本量较少,这可能是导致拟合情况较好的原因。只有粘质粉土Ⅲ类场地的拟合优度明显优于场地未分类,淤泥质粉质粘土、粉砂和砂质粉土则相差不明显,而粉质粘土却相反,说明受岩土类型的影响,考虑场地分类因素并不一定减少误差或提高预测精度。考虑不同工程实际状况的需要,表 3给出场地未分类、Ⅱ类场地和Ⅲ类场地下5类常规土关系公式、残差标准差、适应范围和拟合优度,以供使用选择。

    表 3  各场地类型下不同岩土类型的推荐模型
    Table 3.  Recommended models for different geotechnical types under different site types
    岩土类型 场地类别 推荐模型 回归方程 残差标准差 R2 适应范围/m
    粉质粘土 Ⅱ类 模型2 ${v_s} = 112.4398 + 5.9058H - 0.0077{H^2}$ 13.58 0.9530 2—44
    Ⅲ类 模型3 ${v_s} = 86.5602 + 18.7249{H^{0.6352}}$ 18.45 0.8813 2—78
    场地未分类 模型3 ${v_s} = 85.5503 + 21.3380{H^{0.6041}}$ 19.44 0.8914 2—78
    淤泥质粉质粘土 Ⅱ类 模型2 ${v_s} = 104.8362 + 0.2524{H^{1.9230}}$ 8.83 0.9096 4—24
    Ⅲ类 模型2 ${v_s} = 90.5189 + 4.1075H - 0.0227{H^2}$ 10.46 0.8685 2—46
    场地未分类 模型2 ${v_s} = 89.9598 + 4.1857H - 0.0245{H^2}$ 10.65 0.8656 2—46
    粉砂 Ⅱ类 模型1 ${v_s} = 130.5509 + 5.6417H$ 8.08 0.9594 1—30
    Ⅲ类 模型2 ${v_s} = 100.2213 + 7.3496H - 0.0521{H^2}$ 18.79 0.8498 5—56
    场地未分类 模型2 ${v_s} = 98.7936 + 7.5598H - 0.0558{H^2}$ 21.19 0.8482 1—56
    砂质粉土 Ⅱ类 模型3 ${v_s} = 112.9768 + 4.1008{H^{1.1103}}$ 8.37 0.9773 1—34
    Ⅲ类 模型3 ${v_s} = 120.9837 + 3.0980{H^{1.0546}}$ 18.45 0.9155 1—57
    场地未分类 模型2 ${v_s} = 113.5888 + 4.5896H - 0.0079{H^2}$ 12.73 0.8811 1—57
    粘质粉土 Ⅱ类 模型2 ${v_s} = 115.3530 + 4.2818H + 0.0461{H^2}$ 9.42 0.9551 2—36
    Ⅲ类 模型2 ${v_s} = 121.1832 + 2.7579H + 0.0110{H^2}$ 7.39 0.9291 2—44
    场地未分类 模型2 ${v_s} = 123.4666 + 2.4199H - 0.0354{H^2}$ 13.14 0.8386 2—44
    下载: 导出CSV 
    | 显示表格

    为进一步验证本文推荐模型的可靠性和地区差异性,首先选取未进入统计的新建杭州至黄山铁路(Ⅱ类场地)和钱江新城二期连堡丰城项目(Ⅲ类场地)钻孔资料,将本文推荐模型预测的剪切波速与本地实测钻孔剪切波速进行对比分析,其中相对误差表示预测值与实测值之差占实测值的比值;然后选取国内学者发布的粉质粘土和粉砂剪切波速统计模型,将本文推荐模型与其他模型进行对比分析,以研究地区差异性影响,其中相对误差表示学者模型预测剪切波速与本文模型预测值之差占本文模型预测值的比值。

    (1)Ⅱ类场地对比分析结果如表 4所示,由表 4可知,Ⅱ类场地推荐模型和场地未分类推荐模型下的剪切波速预测值相对误差均在15%以内,最小相对误差仅为0.07%,最大相对误差为14.14%,具有良好的可靠性表现。且考虑场地分类与不考虑场地分类的剪切波速预测精度因岩土类型的不同产生了差异,相比其他土类,淤泥质粉质粘土明显表示出未考虑场地分类的剪切波速预测精度优于考虑场地类别,有别于回归分析中相差不大的规律,说明统计模型的实际预测精度也受岩土类型的影响。

    表 4  Ⅱ类场地实测剪切波速与预测值的对比分析
    Table 4.  Comparative analysis of measured shear wave velocity and predicted value of site Ⅱ
    岩土类型 钻孔编号 深度
    /m
    实测剪切波速
    /m·s-1
    Ⅱ类场地预测值
    /m·s-1
    相对误差
    /%
    场地未分类预测值
    /m·s-1
    相对误差
    /%
    粉质粘土 JZ1 3 130 130 0.07 127 -2.32
    JZ1 5 143 142 -0.84 142 -0.72
    JZ-Ⅳ 8 152 159 4.76 160 5.58
    JZ-Ⅳ 24 222 250 12.68 231 4.09
    JZ-Ⅳ 27 263 267 1.44 242 -8.06
    粉砂 JZ1 29 278 278 -0.06 249 -10.54
    JZ1 24 233 266 14.14 248 6.48
    JZ1 26 286 277 -3.06 258 -9.92
    淤泥质粉质粘土 JZ-Ⅳ 16 154 157 1.97 151 -2.17
    JZ-Ⅳ 17 145 163 12.75 154 6.23
    JZ-Ⅳ 19 159 177 11.62 161 1.03
    下载: 导出CSV 
    | 显示表格

    (2)将推荐模型应用于Ⅲ类场中6个实测深度均超过50m的钻孔DZ1—DZ6场地。对Ⅲ类场地和场地未分类下的预测剪切波速与实测数据进行对比分析,结果如图 9所示。由图 9可知,除淤泥质粉质粘土在23m预测相对误差超过15%外,本文推荐模型预测出的剪切波速与实测剪切波速相对误差基本在15%以内,最小相对误差仅为0.02%,且各土类相对误差10%内的占比达90%以上,具有较高的预测可靠性;区别于Ⅱ类场地,Ⅲ类场中各岩土类型下考虑场地分类的剪切波预测值与未考虑场地分类的相对误差相差不大,但粉质粘土在整个预测深度范围内均表现出未考虑场地分类的预测值略优于考虑场地分类的预测值,而粉砂在25—27m、淤泥质粉质粘土在21—24m也表示出此规律,再次说明因岩土类型和预测深度的不同,考虑场地分类因素能否在实际工程中提高预测精度具有不确定性,应具体项目具体分析。

    图 9  Ⅲ类场地实测剪切波速与预测值的对比分析图
    Figure 9.  Comparative analysis of measured shear wave velocity and predicted value of site Ⅲ

    通过上述对Ⅱ、Ⅲ类场地实测剪切波速的验证,可知虽然个别土类相对误差将近20%(淤泥质粉质粘土),但超过90%的预测剪切波速相对误差在10%以内,因此用本文推荐模型预测出的剪切波速具有较高的可靠性。

    为验证本文推荐模型与其他地区模型的差异,选取最常见的粉质粘土、粉砂,与国内学者研究的模型进行对比分析。选取模型主要包括荣棉水等(2017)基于渤海常见土类剪切波速资料统计分析给出的粉质粘土、粉砂模型,刘红帅等(2010)根据全国大量安全性评价实测钻孔资料统计给出的粉质粘土、粉砂模型及战吉艳等(2009)给出的苏州城区粉质粘土、粉砂模型。为表达简洁,分别用场地未分类、Ⅱ类场地和Ⅲ类场地代表本文推荐模型,渤海-Rong代表荣棉水等(2017)模型,用未分类-Liu、Ⅱ类-Liu和Ⅲ类-Liu代表刘红帅等(2010)模型中场地未分类、Ⅱ类场地和Ⅲ类场地下的经验模型,用苏州-Zhan代表战吉艳等(2009)模型。对比分析结果如图 1011所示。

    图 10  粉质粘土下本文模型与其他模型剪切波速对比分析图
    Figure 10.  Comparative analysis of shear wave velocity of the model proposed in this paper and other models of silty clay
    图 11  粉砂下本文模型与其他模型剪切波速对比分析图
    Figure 11.  Comparative analysis of shear wave velocity of the model proposed in this paper and other models of silty sand

    图 11(a)可知,总体上,在所有场地分类情况下,杭州地区粉质粘土推荐模型的剪切波速预测值明显低于刘红帅等(2010)模型,而与渤海-Rong、苏州-Zhan在部分深度上相差不明显,且在图 11(b)中也有体现:渤海-Rong与本文模型相对误差基本在5%以内,苏州-Zhan在2—29m和本模型相对误差在15%以内,而与刘红帅等(2010)模型在场地未分类情况下的相对误差为18%—51%、Ⅱ类场地下为17%—41%、Ⅲ类场地下为18%—41%,均体现了较大的差异性。说明受地区的影响,粉质粘土不同统计模型具有明显的地区差异性,因此建议应用者在选取模型时,应特别注意地区适应性。

    由图 12可知,总体上,杭州地区粉砂剪切波速预测模型与粉质粘土相比表现出不同的规律性,即本文推荐粉砂模型与其他学者模型表现出更好的一致性,图 12(b)相对误差中也有所体现,即当预测深度超过10m后,各模型相对误差均在20%以内,说明受岩土类型的影响,不同岩土类型对地区的适应性不同。值得注意的是,在0—10m,渤海-Rong、未分类-Liu和Ⅲ类-Liu相对误差均较大,最大相对误差达39.1%,这也说明预测深范围会对模型的地区适应性产生影响。

    综上可知,剪切波速与埋深间统计公式的地区差异性不仅受地区的影响显著,更与岩土类型、预测深度范围有关。

    本文基于杭州市场地剪切波速资料,对杭州市粉质粘土、淤泥质粉质粘土、粉砂、砂质粉土和粘质粉土5类土的剪切波速与埋深经验关系进行相关性分析,主要结论如下:

    (1)给出了杭州市5类土在考虑场地分类和未考虑场地分类情况下的最优拟合公式和适应范围,并对推荐模型进行了可靠性验证,对粉质粘土和粉砂进行了地区差异性分析。

    (2)可靠性验证结果表明:本文剪切波速预测值与实测值相对误差均在20%以内,最小相对误差仅为0.02%,超90%的预测剪切波速相对误差在10%以内,具有较合理的可靠性。

    (3)地区差异性分析结果表明:地区、岩土类型和预测深度均会对剪切波速预测模型的可靠性产生显著影响,故当地区无剪切波速资料时,应优先选用本地区统计模型,如若未有,则需选取已有资料,对选用模型进行岩土类型和适应预测深度范围验证,以保证所选模型的可靠性。

    (4)预测模型回归分析结果表明:同一地区,场地分类并不一定提高经验模型的预测精度。受岩土类型和预测深度的影响,考虑场地分类因素能否在实际工程中提高预测精度也具有不确定性。

    因本文统计模型未考虑土体的物理状态,所以仅适用于未确定土体物理状态的不利情况,需加强深入研究。

  • 图  1  不同坐标系下b值拟合曲线对比

    Figure  1.  Comparison of b value fitting curves under different coordinate systems

    图  2  华北平原地震统计区b值拟合结果

    Figure  2.  Fitting results of b value in seismic statistical region of North China Plain

    图  3  郯庐地震统计区b值拟合结果

    Figure  3.  Fitting results of b value in seismic statistical region of the Tan Lu

    图  4  长江下游-南黄海地震统计区b值拟合结果

    Figure  4.  Fitting results of b value in the seismic statistical region of the lower reaches of the Yangtze River and the South Yellow Sea

    图  5  汾渭地震统计区b值拟合结果

    Figure  5.  Fitting results of b value in the seismic statistical region of the Fenwei

    图  6  地震统计区计算控制点分布情况

    Figure  6.  Distribution of calculation control points in seismic statistical region

    图  7  2套地震活动性参数下地震动峰值加速度差值分布(50年超越概率10%)

    Figure  7.  Distribution of the difference in seismic peak acceleration between two sets of seismic activity parameters in the seismic statistical region (with a probability of exceeding 10% in 50 years)

    表  1  C-S法余震时间

    Table  1.   Aftershock time in C-S method

    主震震级/级时间/天主震震级/级时间/天
    3.0~3.511.56.1~6.5510
    3.6~4.0226.6~7.0790
    4.1~4.5427.0~7.5915
    4.6~5.0837.6~8.0960
    5.1~5.51558.0~8.5985
    5.6~6.0290≥8.6985
    下载: 导出CSV

    表  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.44.5~4.95.0~5.45.0~5.9
    华北平原151011
    郯庐30121
    长江下游-南黄海15113
    汾渭154
    下载: 导出CSV

    表  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.44.5~4.95.0~5.45.5~5.96.0~6.46.5~6.97.0~7.47.5~7.98.0~8.48.5~8.9
    华北平原地震统计区1484—20100.17630.09670.03980.02280.00960.00570.0019
    1500—20100.17790.09960.04110.02350.00970.00590.0020
    1791—20100.22220.09980.05000.03180.01840.0091
    1950—20103.69981.57110.34430.11760.08490.06850.03570.0164
    郯庐地震统计区1477—20100.12170.06930.03560.01500.01120.00370.00190.0019
    1500—20100.12130.06850.03520.01570.01170.00390.00200.0020
    1892—20100.18490.09240.06720.03360.0168
    1970—20103.07321.02440.46340.07320.02440.02440.0244
    长江下游-南黄海地震统计区1485—20100.10650.05320.03420.01710.0019
    1491—20100.10770.05380.03460.01730.0019
    1839—20100.25580.13950.08720.04650.0058
    1970—20102.34150.95120.36590.12200.0976
    汾渭地震统计区1000—20100.02470.02080.00690.00300.0020
    1209—20100.02740.02240.00750.00370.0025
    1484—20100.18030.08540.03230.02470.00950.00380.0019
    1500—20100.17610.08220.02940.02350.00980.00390.0020
    1950-20102.09841.01640.44260.13110.0164
    下载: 导出CSV

    表  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.44.5~4.95.0~5.45.5~5.96.0~6.46.5~6.97.0~7.47.5~7.98.0~8.48.5~8.9
    华北平原地震统计区1484—20230.17590.09630.03890.02220.00930.00560.0019
    1500—20230.17750.09920.04010.02290.00950.00570.0019
    1791—20230.21890.09870.04720.03000.01720.0086
    1950—20233.54051.50000.31080.10810.06760.05410.02700.0135
    郯庐地震统计区1477—20230.12070.06760.03470.01460.01100.00370.00180.0018
    1500—20230.12020.06680.03440.01530.01150.00380.00190.0019
    1892—20230.17420.08330.06060.03030.0152
    1970—20233.12961.01850.37040.05560.01850.01850.0185
    长江下游-南黄海地震统计区1485—20230.10950.05200.03340.01670.0019
    1491—20230.11070.05250.03380.01690.0019
    1839—20230.24320.12970.08110.04320.0054
    1970—20232.31480.98150.33330.09260.0741
    汾渭地震统计区1000—20230.02440.02050.00680.00290.0020
    1209—20230.02700.02210.00740.00370.0025
    1484—20230.17590.08330.03150.02410.00930.00370.0019
    1500—20230.17180.08020.02860.02290.00950.00380.0019
    1950—20231.98650.89190.36490.10810.0135
    下载: 导出CSV

    表  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减小
    下载: 导出CSV

    表  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计算得到的基岩地震动峰值加速度之差。
    下载: 导出CSV
  • 陈凌,刘杰,陈颙等,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/pg2021EE0025

    Meng 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.006

    Pan 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.002

    Pan 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.005

    Ren 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.002

    Wang 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.002

    Wu 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.009

    Xie 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.013

    Xie 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.027

    Zhang 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
  • 加载中
图(7) / 表(6)
计量
  • 文章访问数:  180
  • HTML全文浏览量:  78
  • PDF下载量:  41
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-01-31
  • 刊出日期:  2024-06-30

目录

/

返回文章
返回