Study on Minimum Integrity Magnitude of Earthquake Catalogue in Jiangsu Seismic Network
-
摘要: 区域台网完整性震级的研究对于了解测震台网监测能力、评估地震目录的可靠性以及研究地震活动特征具有重要意义,其有助于确定最小可靠震级范围,并为地震监测和地震研究提供参考。针对江苏及邻区地震目录的特点和需求,结合江苏省区域测震台网台站建设与发展情况,使用震级-序号法、多参数方法联合研究,对江苏区域台网记录的地震活动最小完整性震级
$ M_{\mathrm{C}} $ 随时间的变化趋势进行分析;通过使用完整震级范围(EMR)方法,对地震的空间分布特征进行深入探究。研究结果表明,江苏及邻区整体区域$ {M}_{\mathrm{C}} $ 值变化趋势大致分为2个阶段,1970—2008年$ {M}_{\mathrm{C}} $ 值主要分布于$ {M}_{\mathrm{L}}1.0\mathrm{~}{M}_{\mathrm{L}}2.2 $ 之间,2009—2022年$ {M}_{\mathrm{C}} $ 值主要分布于$ {M}_{\mathrm{L}}0.5\mathrm{~}{M}_{\mathrm{L}}1.5 $ 之间。江苏内陆地区1970—2008年$ {M}_{\mathrm{C}} $ 值主要分布于$ {M}_{\mathrm{L}}1.0\mathrm{~}{M}_{\mathrm{L}}2.5 $ 之间,2009—2022年$ {M}_{\mathrm{C}} $ 值主要分布于$ {M}_{\mathrm{L}}0.5\mathrm{~}{M}_{\mathrm{L}}2.0 $ 之间,江苏沿海黄海海域1970—2022年$ {M}_{\mathrm{C}} $ 值主要分布于$ {M}_{\mathrm{L}}2.0\mathrm{~}{M}_{\mathrm{L}}3.0 $ 之间,江苏区域台网地震监测能力趋于完善。-
关键词:
- 最小完整性震级MC值 /
- 区域地震目录 /
- 江苏省区域测震台网
Abstract: Research on the completeness magnitude (MC) of regional seismic networks is crucial for understanding their monitoring capabilities, assessing the reliability of earthquake catalogs, and studying seismic activity characteristics. This research helps to determine the minimum reliable magnitude range, and provide valuable references for earthquake monitoring and research efforts. In this study, we focus on the earthquake catalog characteristics and requirements for Jiangsu province and its neighboring areas, taking into account the construction and development of seismic stations within the regional seismic network. We employ the magnitude-number method and multi-parameter approaches to investigate and analyze the temporal variations of the minimum completeness magnitude (MC) of seismic activity recorded by this network. Using the Entire Magnitude Range (EMR) method, we further explore the spatial distribution characteristics of earthquakes. Our findings indicate that the overall trend of MC values in Jiangsu and its surrounding areas can be divided into two distinct stages: from 1970 to 2008, MC values predominantly ranged from ML1.0 to ML2.2, while from 2009 to 2022, MC values mainly fell between ML0.5 and ML1.5. Inland areas of Jiangsu Province exhibited MC values ranging from ML1.0 to ML2.5 during the period of 1970—2008, and from ML0.5 to ML2.0 during 2009—2022. In contrast, the MC values in the coastal Yellow Sea area of Jiangsu primarily ranged from ML 2.0 to 3.0 throughout the entire period from 1970 to 2022. Overall, the seismic monitoring capabilities of the Jiangsu regional seismic network appear to be improving over time. -
引言
1995年阪神地震中大开地铁车站的严重倒塌破坏在世界范围内引起了地震研究学者们对地下结构抗震问题的广泛关注。在地下结构地震反应特征和抗震设计方法等问题的研究过程中,逐渐形成了一系列实用抗震设计方法(刘晶波等,2010a;许成顺,2017)。其中,反应位移法与反应加速度法凭借理论明确、操作简便等优点,在实际工程设计中得到广泛应用,已被《城市轨道交通结构抗震设计规范》(GB 50909—2014)(中华人民共和国住房和城乡建设部,2014)和《地下结构抗震设计标准》(GB/T 51336—2018)(中华人民共和国住房和城乡建设部,2018)采纳。
反应加速度法通过将顶、底板相对位移最大时刻的自由场水平加速度施加于土-结模型上,计算地下结构变形与内力。与反应位移法相比,主要优势为:反应加速度法建立的土-结构模型可较准确地反映土-结构间相互作用,无须计算地基弹簧刚度系数,大大提高了计算效率与精度;反应位移法受限于地基弹簧与地震作用的施加方式,决定了其主要适用于矩形断面,而反应加速度法施加的地震作用为体积力,且模型采用土-结构整体模型,决定了其可适用于复杂断面地下结构的地震反应计算。总的来说,反应加速度法优于反应位移法(刘晶波等,2013;李新星等,2014)。
反应加速度法的应用关键在于地震输入荷载的确定。为简化计算过程,反应加速度法假设输入荷载为地震作用下自由场的反应加速度,然而,该假设忽略了结构的存在对自由场反应加速度的影响。基于此,刘如山等(2007)从一维自由土层反应的剪应力入手,首先对一维自由土层反应剪应力沿竖直方向进行微分,然后将其作为水平体荷载离散到有限元节点上,以此作为加载方法;董正方等(2014)运用柔度系数法中的结构变形和自由场变形之间的相互作用系数,修正了自由场地震反应加速度,近似考虑了结构的存在对反应加速度法的影响;Liu等(2018)通过理论解析方式推导了自由场模型与土-结模型之间场地反应加速度的关系,结果表明可通过2种模型场地相对位移修正场地反应加速度法;徐琨鹏等(2019)对比了反应加速度法和仅在土体施加加速度的简化方法,发现2种方法对侧边距取值不敏感;禹海涛等(2020)通过动力时程分析,提出地下结构附近存在动力响应增强区和衰减区,利用动力响应调整系数提高了反应加速度法的计算精度。
为考虑结构的存在对反应加速度法中地震输入荷载的影响,本文提出将土-结构模型转换为自由场模型的等效方式,采用等效后的土层参数进行自由场反应加速度计算,以此近似地反映结构的存在对地震输入荷载的影响,进而改进反应加速度法中地震输入荷载的计算过程。采用有限元软件ABAQUS对大开车站地震反应进行数值模拟分析,验证等效方式的有效性,并探讨等效模型宽度取值范围对结果的影响,初步给出等效模型宽度范围。
1. 反应加速度法
采用改进前反应加速度法对地下结构横断面进行抗震分析时,计算模型为图1所示的土-结构整体模型。一般地,土体采用二维平面应变单元模拟,结构采用梁单元模拟。将模型底边固定,模型两侧设置水平滑移边界。通过将地下结构顶、底板位置处土层相对位移最大时刻的各土层水平加速度或有效惯性加速度施加于土-结构整体模型上,获得地下结构变形与内力等反应。
反应加速度法主要计算步骤如下:
(1)一维土层地震反应分析。采用等效线性化程序EERA或SHAKE91对一维土层进行地震反应分析,记录结构顶、底板位置处土层最大相对位移时刻,各土层加速度及等效剪切模量Gi。《城市轨道交通结构抗震设计规范》建议,土层简单时地震输入荷载直接采用自由场水平加速度,土层复杂时地震输入荷载采用土层有效惯性加速度。
(2)建立土-结构模型。《城市轨道交通结构抗震设计规范》建议模型底面可取设计地震作用基准面,顶面取地表面,侧面边界到结构的距离宜取结构水平有效宽度的2~3倍,本文针对模型宽度对计算结果的影响展开进一步讨论。
(3)施加等效地震作用。将土层水平加速度或有效惯性加速度按各土层单元的位置施加于土-结构模型中的各土层上,结构部分按所在土层深度位置作用相应的加速度。
从反应加速度法计算模型和计算步骤可看出,该方法地震输入荷载采用的是自由场反应加速度,未考虑结构的存在对自由场加速度反应的影响,尤其是当地下结构刚度和周围土体刚度相差较大时,结构的存在可能会在较大程度上影响周围土体的地震反应。此外,在动力学理论方面,该方法忽略了阻尼因素的影响,仅是某种程度的近似,当土层反应仅有一阶振型发育时,该方法精度尚可,当遇到复杂土层时误差将增大。
2. 改进反应加速度法研究
2.1 等效思路
为考虑结构的存在对反应加速度法中地震输入荷载的影响,对反应加速度法中的加速度求解方法进行改进。首先,按结构周围土体的分层将土-结构模型中的结构与土体等效为自由场,计算等效自由场的土层参数;然后,采用等效线性化程序EERA将等效自由场土层参数进行处理,得到土层等效剪切模量,并计算等效自由场各土层反应加速度,如图2所示。
将图2所示土-结构模型等效为自由场模型时需进行结构等效和土层等效,其中Ki为等效后土体的刚度;K1i为结构的刚度;K2i为等效前土体的刚度。对于地下结构横断面,其顶、底板在断面水平方向上连续,侧墙、中柱在断面水平方向上间断,而土层为分层结构,在断面水平方向上连续。为便于土体与结构进行等效,将断面水平方向上间断的侧墙与中柱等效为连续结构,形成图3(a)所示的分层结构形式。在此基础上,按结构周围土层的划分,将分层结构与分层土体按照相对位置等效为分层的自由场,如图3(b)所示。
上述2个等效过程均包含等质量等效计算和等侧向刚度等效计算。自由场在地震作用下的反应加速度主要受场地质量分布与刚度分布的影响,等质量等效计算假设等效前后土-结构模型与自由场模型中各土层质量保持不变:
$$ {\rho _i}{V_i} = {\rho _{1i}}{V_{1i}} + {\rho _{2i}}{V_{2i}} $$ (1) 式中,
${\rho _i}$ 为等效后第i层土体密度;${V_i}$ 为采用土-结构模型总宽度计算的第i层土体体积;${\rho _{1i}}$ 为第i层土体对应的结构部分的密度;${V_{1i}}$ 为第i层土体中包含的结构体积;${\rho _{2i}}$ 为第i层土体密度;${V_{2i}}$ 为减去结构体积的第i层土体体积。等侧向刚度等效计算中,假设等效前后各土层侧向刚度保持不变,各土层与结构为短柱,进行侧向刚度的串并联叠加:
$$ {G_i}{b_i} = {G_{1i}}{b_{1i}} + {G_{2i}}{b_{2i}} $$ (2) 式中,
${G_i}$ 为等效后第i层土体剪切模量;${b_i}$ 为采用土-结构模型总宽度;${G_{1i}}$ 为结构剪切模量;${b_{1i}}$ 为第i层土体中结构的宽度;${G_{2i}}$ 为第i层土体剪切模量;${b_{2i}}$ 为减去结构宽度的第i层土体宽度。2.2 等效模型宽度的选取
为便于区分,将改进前反应加速度法称为反应加速度法-1(以下简称“RAM-1”),将改进后反应加速度法称为反应加速度法-2(以下简称“RAM-2”)。
采用图4所示土-结构模型边界至结构侧边的距离B表示模型宽度。在RAM-2计算过程中,反应加速度法拟静力计算时涉及土-结构模型宽度(以下简称“模型宽度-F”)、等效土-结构模型宽度(以下简称“模型宽度-B”)的选取。由前文可知,模型宽度-B不同,将直接改变等效自由场土层参数,影响等效自由场场地反应,从而影响反应加速度法计算结果。因此,本文针对不同模型宽度-B对反应加速度法计算结果的影响进行了探讨。
3. 实例分析
3.1 计算模型与参数选取
为验证等效方式的有效性,以大开车站为研究对象,进行了RAM-1、RAM-2与动力时程方法计算。结构断面尺寸如图5所示,结构埋深4.8 m,周围土层参数如表1所示(矢的照夫等,1996)。结构中柱沿车站纵向按3.5 m等间距分布,将其按一定原则等效为连续墙体,等效后中柱弹性模量为8.57 GPa,密度为714 kg/m3。顶、底板与侧墙弹性模量为30 GPa,密度为2 500 kg/m3。土-结构整体分析模型如图6所示,车站结构采用梁单元模拟,土体采用平面应变单元模拟。土体底边界固定,两侧边界设置水平滚轴边界。假设土体与结构之间接触良好,且不产生相对滑移,数值计算时土体与结构之间采用绑定约束。
表 1 土层参数表Table 1. Physical properties of soils分层 土质 厚度/m 密度/(t·m−3) 剪切波速/(m·s−1) 泊松比 1 填土 1 1.9 140 0.333 2 砂土 4.1 1.9 140 0.488 3 砂土 3.2 1.9 170 0.493 4 黏土 3.1 1.9 190 0.494 5 黏土 5.8 1.9 240 0.490 6 砂土 22 2.0 330 0.487 由于缺少大开车站周围土体动力参数,本文采用典型砂土和黏土剪切模量比、阻尼比与剪应变关系曲线(杜修力等,2018),如图7所示。考虑到地震动特性对地下结构地震反应的影响,基岩输入地震动选用Kobe地震动与El Centro地震动,地震动加速度时程曲线如图8所示,将加速度幅值调整为0.1、0.2、0.3 g。
验证等效思路有效性时,RAM-1与RAM-2模型宽度-F均取3倍结构宽度,RAM-2模型宽度-B取3倍结构宽度。动力时程分析共6种计算工况,反应加速度法共12种计算工况(包括RAM-1、RAM-2各6种计算工况)。
在探讨不同模型宽度-B对RAM-2计算结果的影响时,RAM-1与RAM-2模型宽度-F仍取为3倍结构宽度,RAM-2模型宽度-B分别取为1、2、3、4、5、7倍结构宽度。在0.1 g幅值的Kobe地震动作用下进行动力时程计算、RAM-1计算与不同模型宽度-B的RAM-2计算。
3.2 结果与分析
刘晶波等(2010b)和杜修力等(2017)指出地下结构中柱为抗震关键构件,在此基础上,计算不同地震作用下中柱与顶、底板位置处的相对位移、中柱底部弯矩、中柱底部剪力和侧墙底部弯矩,结果如表2所示,表中括号内为相应的误差。
表 2 不同地震作用下结构内力与变形Table 2. Internal force and deformation of structure under different earthquake load地震动 幅值 计算方法 顶、底板相对位移/mm 中柱底部弯矩/
[(kN·m)·m−1]中柱底部剪力/
(kN·m−1)侧墙底部弯矩/
[(kN·m)·m−1]Kobe 0.1 g 动力时程分析法 5.77 36.12 11.09 358.40 RAM-1 5.56(3.66%) 34.06(5.70%) 10.48(5.50%) 341.50(4.72%) RAM-2 5.71(1.04%) 35.00(3.10%) 10.78(2.80%) 351.20(2.01%) 0.2 g 动力时程分析法 16.93 102.70 30.88 887.30 RAM-1 16.17(4.47%) 99.31(3.30%) 30.80(0.26%) 879.50(0.88%) RAM-2 16.44(2.86%) 102.00(0.68%) 30.51(1.20%) 877.00(1.16%) 0.3 g 动力时程分析法 31.15 187.50 55.97 1476.00 RAM-1 30.30(2.76%) 181.70(3.09%) 54.43(2.75%) 1447.00(1.96%) RAM-2 30.94(0.67%) 185.50(1.07%) 55.44(0.95%) 1474.00(0.14%) EL Centro 0.1 g 动力时程分析法 14.12 86.02 26.09 771.40 RAM-1 13.41(5.02%) 81.68(5.05%) 25.04(4.02%) 741.00(3.94%) RAM-2 13.61(3.62%) 83.07(3.43%) 25.59(1.92%) 755.90(2.01%) 0.2 g 动力时程分析法 26.73 160.40 47.59 1288.00 RAM-1 25.36(5.15%) 151.70(5.42%) 45.47(4.45%) 1241.00(3.65%) RAM-2 26.12(2.30%) 156.30(2.56%) 46.82(1.62%) 1278.00(0.78%) 0.3 g 动力时程分析法 38.75 231.60 68.07 1755.00 RAM-1 35.48(8.44%) 210.60(9.07%) 62.24(8.56%) 1628.00(7.24%) RAM-2 37.51(3.22%) 222.10(4.10%) 65.35(4.00%) 1704.00(2.91%) 由表2可知,以动力时程为基准,不同工况下RAM-1计算的结构变形与内力最大误差为9%,不同工况下RAM-2计算的结构变形与内力最大误差为4.1%,等效方式将反应加速度法的计算误差由10%以内缩小至5%以内,提高了反应加速度法计算精度,表明了等效方式的有效性;RAM-2在不同幅值Kobe地震动作用下计算的结构变形与内力最大误差为3.1%,在不同幅值El Centro地震动作用下计算的结构变形与内力最大误差为4.1%,最大误差随地震动强度的改变未明显变化,说明该方法的计算结果较稳定。
以动力时程为基准计算的结构变形与内力、不同模型宽度-B 下 RAM-2 计算的结构变形与内力如表3所示,表中括号内为相应的误差。由表3、图9可知,与动力时程分析法相比,随着模型宽度-B的减小,RAM-2计算结果逐渐增加,逐渐接近动力时程计算结果,甚至大于动力时程计算结果;而RAM-1最大计算误差为5.7%,模型宽度-B为7、5、4、3倍结构宽度时,RAM-2最大计算误差分别为6.15%、4.51%、4.1%、3.1%,表明当模型宽度取值较大时,改进反应加速度法计算误差呈增大趋势,这是由于模型宽度越大,等效后的自由场模型更接近原始自由场模型,此时等效模型无法充分反映结构的存在造成的影响。结合表2所示不同工况下RAM-1计算的结构变形与内力最大误差为9%,不同工况下模型宽度-B为3倍结构宽度时,RAM-2最大计算误差为4.1%,若以计算误差减小50%为衡量标准,建议模型宽度-B取为1~3倍结构宽度。
表 3 不同模型宽度-B下结构内力与变形Table 3. Internal force and deformation of structure with different model width计算方法 模型宽度-B/
结构宽度顶、底板相对
位移/mm中柱底部弯矩/
[(kN·m)·m−1]中柱底部剪力/
(kN·m−1)侧墙底部弯矩/
[(kN·m)·m−1]动力时程分析法 — 5.77 36.12 11.09 358.40 RAM-1 — 5.56(3.66%) 34.06(5.70%) 10.48(5.50%) 341.50(4.72%) RAM-2 1 5.98(3.67%) 36.67(1.52%) 11.31(1.98%) 368.80(2.90%) 2 5.78(0.16%) 35.42(1.94%) 10.92(1.53%) 355.50(0.81%) 3 5.71(1.04%) 35.00(3.10%) 10.78(2.80%) 351.20(2.01%) 4 5.65(2.05%) 34.64(4.10%) 10.68(3.70%) 347.70(2.99%) 5 5.63(2.49%) 34.49(4.51%) 10.64(4.06%) 346.20(3.40%) 7 5.53(4.16%) 33.90(6.15%) 10.46(5.68%) 340.20(5.08%) 4. 结论
本文采用有限元分析软件ABAQUS对大开车站进行不同工况下的动力时程分析和改进前后反应加速度法分析,验证了等效方式对反应加速度法的有效性。在验证等效思路有效的基础上,进行了不同等效模型宽度的改进反应加速度法计算,初步给出了等效模型宽度的选取范围。
(1)在未改变RAM-1理论基础的前提下,对关键参数——地震输入荷载进行了修正。与动力时程分析法相比,修正后的反应加速度法使不同工况下的最大计算误差由10%以内缩小至5%以内,提高了反应加速度法计算精度。
(2)与动力时程分析结果相比,修正后的反应加速度法在不同工况下的最大计算误差未明显改变,说明修正结果具有稳定性。
(3)对不同等效模型宽度对反应加速度法计算结果的影响进行分析,根据结果建议等效模型边界至结构侧边的距离取为1~3倍结构宽度。
-
表 1 江苏及邻区地震目录
Table 1. Statistics of earthquake catalogues in Jiangsu and its adjacent areas
编号 地震区(带) 地震数量/次 地震区(带)内最大震级/级 1 江苏及邻区整体区域 24 156 6.4 2 江苏内陆陆地整体区域 5 835 6.3 3 黄海海域整体区域 3 009 6.4 4 苏北盆地 1 825 5.4 5 苏南隆起 2 461 6.3 -
冯建刚,蒋长胜,韩立波等,2012. 甘肃测震台网监测能力及地震目录完整性分析. 地震学报,34(5):646−658. doi: 10.3969/j.issn.0253-3782.2012.05.006Feng J. G., Jiang C. S., Han L. B., et al., 2012. Analysis on the monitoring capability of seismic networks and completeness of earthquake catalogues in Gansu region. Acta Seismologica Sinica, 34(5): 646−658. (in Chinese) doi: 10.3969/j.issn.0253-3782.2012.05.006 何奕成,范小平,赵启光等,2021. 郯庐断裂带中南段地壳结构分段特征. 地球物理学报,64(9):3164−3178. doi: 10.6038/cjg2021O0247He Y. C., Fan X. P., Zhao Q. G., et al., 2021. Segmentation of crustal structure beneath the middle-south segment of Tan-Lu Fault Zone. Chinese Journal of Geophysics, 64(9): 3164−3178. (in Chinese) doi: 10.6038/cjg2021O0247 立凯,何奕成,宫杰等,2023. PMC方法在江苏测震台网监测能力评估中的应用. 震灾防御技术,18(3):642−650. doi: 10.11899/zzfy20230321Li K., He Y. C., Gong J., et al., 2023. Analysis of monitoring capability of Jiangsu seismic network based on PMC method. Technology for Earthquake Disaster Prevention, 18(3): 642−650. (in Chinese) doi: 10.11899/zzfy20230321 史翔宇,王晓青,邱玉荣等,2020. 川滇地震科学实验场地震目录最小完整性震级分析. 地球物理学报,63(10):3683−3697. doi: 10.6038/cjg2020O0136Shi X. Y., Wang X. Q., Qiu Y. R., et al., 2020. Analysis of the minimum magnitude of completeness for earthquake catalog in China Seismic Experimental Site. Chinese Journal of Geophysics, 63(10): 3683−3697. (in Chinese) doi: 10.6038/cjg2020O0136 王亚文,蒋长胜,2017. 南北地震带地震台网监测能力评估的不同方法比较研究. 地震学报,39(3):315−329. doi: 10.11939/jass.2017.03.002Wang Y. W., Jiang C. S., 2017. Comparision among different methods for assessing monitoring capability of seismic station in North-South Seismic Belt. Acta Seismologica Sinica, 39(3): 315−329. (in Chinese) doi: 10.11939/jass.2017.03.002 余娜,张晓清,杨晓霞,2020. 青海及邻区地震目录最小完整性震级分析. 地震,40(4):23−32. doi: 10.12196/j.issn.1000-3274.2020.04.002Yu N., Zhang X. Q., Yang X. X., 2020. Analysis of minimum magnitude of completeness for earthquake catalog in Qinghai region. Earthquake, 40(4): 23−32. (in Chinese) doi: 10.12196/j.issn.1000-3274.2020.04.002 Huang Q. H., 2006. Search for reliable precursors: a case study of the seismic quiescence of the 2000 western Tottori prefecture earthquake. Journal of Geophysical Research: Solid Earth, 111(B4): B04301. Rydelek P. A., Sacks I. S., 1989. Testing the completeness of earthquake catalogues and the hypothesis of self-similarity. Nature, 337(6204): 251−253. doi: 10.1038/337251a0 Schorlemmer D., Mele F., Marzocchi W., 2010. A completeness analysis of the national seismic network of Italy. Journal of Geophysical Research: Solid Earth, 115(B4): B04308. Wiemer S., Wyss M., 2000. Minimum magnitude of completeness in earthquake catalogs: Examples from Alaska, the Western United States, and Japan. Bulletin of the Seismological Society of America, 90(4): 859−869. doi: 10.1785/0119990114 Woessner J., Wiemer S., 2005. Assessing the quality of earthquake catalogues: estimating the magnitude of completeness and its uncertainty. Bulletin of the Seismological Society of America, 95(2): 684−698. doi: 10.1785/0120040007 Wyss M., Hasegawa A., Wiemer S., et al., 1999. Quantitative mapping of precursory seismic quiescence before the 1989, M7.1 off-Sanriku earthquake, Japan. Annals of Geophysics, 42(5): 851−869. 期刊类型引用(0)
其他类型引用(1)
-