Study on Seismic Hazard Assessment Algorithm Considering Seismogenic Sources of Subduction Zones and Its Application on Chinese Marine Territory
-
摘要:
本文研究了俯冲带潜在震源区离散化方法及考虑俯冲带高震级地震震源破裂面和震源深度的场点地震动计算方法,推导了俯冲带潜在震源区地震危险性计算公式,并使用中国海域及邻区地震危险性模型进行地震危险性试算。结果表明,本文建立的考虑俯冲带潜在震源区的地震危险性算法能够实现场点地震危险性计算时对俯冲带高震级地震震源破裂面和震源深度的考虑。
Abstract:This paper studies the discretization method of the seismogenic sources of the subduction zone and the ground motion calculation method considering the source rupture surface and focal depth of the high-magnitude earthquakes in the subduction zone. The calculation formula for the seismic hazard of the seismogenic sources of the subduction zone is derived, and the preliminary seismic hazard calculation for Chinese marine territory is tested. The main conclusions are: the seismic hazard algorithm established in this paper considering the seismogenic sources of the subduction zone can realize the consideration of the source rupture surface and focal depth of the high-magnitude earthquake in the subduction zone in the seismic hazard calculation of the site.
-
Key words:
- Seismic hazard /
- Seismic hazard mapping /
- Marine earthquake /
- Subduction zone /
- 3D seismogenic sources
-
引言
岩石圈磁场是地磁场的重要组成部分,来源于岩石圈居里等温面以上的岩石磁性,其变化受地质结构、构造环境、构造活动、温度与应力状态等多种因素影响。诸多学者已对岩石圈磁场进行相关研究工作,对于空间尺度而言,涵盖全球范围(徐文耀,2003;徐文耀等,2008)、中国大陆地区(杜劲松等,2017;冯丽丽等,2015;王粲等,2017;Ou等,2013)、局部地区(文丽敏等,2017;韶丹等,2015;徐晓雅等,2017;宋成科等,2017;Hemant等,2009)、特定构造带(张忠龙等,2017),数据来源主要为地磁场模型、卫星磁测、航空磁测、地面矢量磁测等。
相对于卫星及航空磁测,地面测量更接近岩石圈磁场物理源,且具备更高的测量精度。在中国地震局经费支持下,流动地磁工作团队已在“冀晋蒙”“京津冀”局部地区开展了流动地磁总强度加密监测工作,布设了规模较大的高精度、高密度测点网络。监测区地处“首都圈”地震重点监视区,地质构造复杂,历史地震多发。观测作业以地磁场总强度F为单一测量要素,其测量精度高于卫星、航空磁测,测点密度优于地面矢量磁测,为相关研究工作提供了高质量数据来源。本文应用相关观测资料,计算观测区域高分辨率F要素岩石圈磁场数值模型(以ΔF表示),研究ΔF空间分布、地震地质特征,并重点针对历史地震震中在ΔF中的位置特征进行研究,以期为该地区震中预测工作提供新技术思路和资料依据。
1. 数据资料与数据处理
1.1 数据资料
观测区内共有180个地磁总强度测点,平均点间距约28km,由“晋冀蒙”及“京津冀”2个地磁总强度加密测网合并而成,观测作业执行《流动地磁测量基本技术要求》1及其规范性附录。为合理控制测区边界处模型产出结果,笔者在测区周边选取40个流动地磁矢量测点,作为模型边界约束点,并将其中4个测点用作补充测点,令研究区最终测点总数达184个(图 1)。
1 中国地震局监测预报司,中震测函[2015]39号:关于印发《流动地磁测量基本技术要求(试行)》的通知
1.2 数据处理
本文采用最新一期观测数据,通过日变通化改正、主磁场长期变改正、主磁场剥离等主要技术处理,获取研究区ΔF(陈斌等,2017),过程如下:
(1)日变通化改正
为消除流动地磁观测数据中包含的地磁场日变化等外源场成分,依托测区邻近地磁台站连续观测分钟值数据,采用单台参照法,对野外观测数据进行日变通化处理,获取监测区日变通化改正数据集。
(2)长期变改正
采用1995年1月1日以来的全国地磁台网观测数据,建立中国及周边地区地磁场长期变化3阶NOC非线性模型(顾左文等,2009),对日变通化改正数据集进行主磁场长期变改正,获取监测区长期变改正数据集。
(3)主磁场剥离
以IGRF12(2015.0年代)为研究区主磁场参考模型,并在监测区长期变改正数据结果中进行剔除,即获得各测点ΔF初始值。
(4)低通滤波
地面磁测结果必然包含杂乱的浅表岩层磁性影响成分(图 2)。为抑制和消除该数据成分,保证研究成果的合理性,笔者采用移动平均法对网格化后的ΔF初始值进行低通滤波处理,滤波器为5×5节点,对应的实际空间范围为0.5°×0.5°,滤波结果即为本研究最终采用的ΔF,本文研究模型描述了研究区2015.0年代ΔF空间分布形态(图 3)。本文ΔF等值线图中,红、蓝色实线分别代表正、负值区,黑色加粗实线为“0”值线,等值线间隔统一为25nT。
2. ΔF空间分布特征
2.1 总体特征
研究区ΔF呈正、负异常区零散相间的分布形态,负异常区面积略大于正异常区。ΔF空间结构复杂,局部异常特征表现强烈,强磁异常区形状多表现为条带状或团块状。研究区范围内共有184个测点,ΔF数值范围为-326.0—180.1nT,平均值为-40.0nT。其中负值测点128个,占69.6%,ΔF平均值为-83.1nT;正值测点56个,占30.4%,ΔF平均值为58.6nT。全部测点ΔF平均幅值(为表述方便,本文借用物理学名词“幅值”表示ΔF绝对值)为75.6nT,约占地磁场能量成分的0.14%(图 4)。
2.2 分区特征
根据ΔF形态特征,可在京南地区划出1条NE走向界线,将研究区分为以下2个分区(图 4):
(1)中西部为ΔF能量集中区,其核心区由2个NE走向的高值负异常区及其间的高值正异常区组成,近似平行四边形。磁异常强度总体较高,且随空间展布快速变化,磁场高梯级带广泛分布。
(2)中东部为弱磁异常区,表现为幅值较低的背景性负异常区,在其边缘零散分布小片正磁异常区,最东端为局部强磁异常区。
ΔF分区与地形、地层分布形态具有明显相关性,表现为三者具备相似的分区特征,且分区界线位置基本吻合:
(1)分区界线与太行山脉东缘基本重合,其东侧为平原地区,磁异常强度总体较低;西侧为山区及高原地区,磁异常强度整体较高且形态复杂。山区强磁异常区呈条带状,走向基本与山脉吻合(图 5)。该现象的形成应与平原地区深厚冲积层对ΔF短波成分的屏蔽效应及山地浅表岩层磁性的影响直接相关。
(2)以分区界线为界,东部地区以第四系地层为主,面积较大、形状规整,与ΔF弱磁异常区基本重合;西部地区地层多样、分布杂乱,均为较古老的地层,火山岩、变质岩岩体及碎屑分布较广泛,与ΔF强磁异常区对应(图 6)。该特征显示了第四系地层地表磁性与其它地层明显不同。
3. 研究区ΔF地震地质特征
3.1 ΔF与二级地块构造
研究区位于燕山、华北平原及鄂尔多斯3个二级地块交汇部位,地块间的结合带贯穿其间。各地块在ΔF中显示出不同的特征:华北平原地块磁异常强度较弱,其西北边界线与ΔF中NE走向高梯级带及“0”值线具有较好的吻合度,东北边界线穿过ΔF弱磁异常区;鄂尔多斯地块东北局部为强磁异常区,边界线局部与ΔF高梯级带有一定程度的对应;燕山地块南缘贯穿ΔF正、负异常相间分布区(图 7)。地块结合带部位为强磁异常区,形态复杂。
3.2 ΔF与断裂构造
研究区内断裂构造广泛发育,分布密集。中、西部断裂为晚更新世—全新世(距今10—12万年)以来的活动断裂,其走向以NE向为主,与该地区强磁异常条带的空间位置大体对应,走向基本一致;研究区中、东部断裂属平原区隐伏断裂,对应ΔF弱磁异常区;研究区东端局部强磁异常区与对应断裂带的走向基本一致,同为NEE向(图 8)。
3.3 历史地震在ΔF中的位置特征
研究区涵盖张渤地震带大部分及华北平原地震带北端,为华北地区重要的地震活动区,历史上共发生Ms5.0以上地震85次,其中包括Ms6.0以上地震20次、Ms7.0以上地震4次(图 9)。为全面研究历史地震震中在ΔF中的分布规律与位置特征,笔者按震级分类,对震中位置处ΔF值、ΔF特殊位置处的震中数量进行分类统计,结果见表 1。
表 1 研究区历史地震震中部位ΔF数值统计结果Table 1. Basic statistical results of ΔF value at the epicenter of historical earthquakes in the study area项目 震级 5级以上 6级以上 7级以上 地震总数/个 85 20 4 震中位置处ΔF均值/nT -70.6 -75.6 -60.2 震中位置处ΔF平均幅值/nT 82.9 91.0 60.2 位于ΔF负值区的震中个数与占比/(个,%) 70,82.4% 17,85.0% 4,100.0% ΔF低幅值(小于平均幅值)区震中个数与占比/(个,%) 51,60.0% 11,55.0% 4,100.0% ΔF高幅值(大于2倍平均幅值)区震中个数与占比/(个,%) 10,11.8% 2,10.0% 0,0.0% ΔF“0”值线附近(两侧各15km范围内)震中个数与占比/(个,%) 50,58.8% 8,40.0% 2,50.0% ΔF高梯级带处震中个数与占比/(个,%) 70,82.4% 16,80.0% 3,75.0% 根据统计结果,对历史地震震中在ΔF中的位置特征进行以下分析:
(1)震中位置处ΔF平均值明显低于研究区总体平均值(-40nT);
(2)ΔF负值区内地震个数远多于正值区,且随着震级的增大,其比例增加;
(3)低幅值区内震中个数明显多于高幅值区。ΔF超高幅值区(大于2倍平均幅值)内震中数量占比较低,且随着震级的增大,其占比进一步减小;
(4)在ΔF“0”值线附近的狭长地带中,密集发生了大量历史地震;
(5)绝大多数地震发生于ΔF高梯级带部位;
(6)强磁异常区内,地震震中表现出的位置特征明显区别于弱磁异常区。前者以ΔF负值区、“0”值线、高梯级带为主;后者主要表现为ΔF负值区、低幅值区。相对而言,前者的震中位置特殊性显著程度远大于后者。
4. 总结与讨论
通过本文研究主要得到以下结论:
(1)研究区ΔF空间结构复杂,局部异常特征表现强烈。总体具有分区特征,且与地形、地层的关系密切。
(2)磁场高梯级带走向与断裂构造体系总体一致;各二级地块ΔF特征明显不同,地块边界线局部与ΔF高梯级带吻合;地块间结合带部位为强磁异常区,形态复杂。
(3)5级以上历史地震与ΔF负值区、低幅值区、“0”值线、高梯级带等特殊位置关系密切。强磁异常区内的震中位置统计规律性明显强于弱磁异常区。该特征可为本地区未来震中位置预测工作提供新的参考依据。
在已有研究的基础上,计算研究区岩石圈磁场数值模型时,应用了较大规模的高精度、高密度地面磁测资料,因此模型计算与研究结论具备更好的数据基础。进而对该模型进行了较细致的量化研究,并从“为地震预测服务”角度出发,分层次、多角度对历史地震在ΔF中的分布规律进行统计分析,并进行有针对性的研究工作。
在数据处理过程中,笔者发现地面磁测结果中包含地表浅层岩层磁性影响成分,并对ΔF初始值进行低通滤波处理。为抑制和消除其影响,相关学者在进行类似研究工作时,也应采用适当的数据处理方法对其进行剔除处理,以保证研究结果的合理性。
受研究区域范围及笔者研究水平的限制,研究结论可能具有一定片面性,更深入、全面的研究工作有待进一步开展。
-
-
邓起东, 张裕明, 环文林等, 1980. 中国地震烈度区划图编制的原则和方法. 地震学报, 2(1): 90—110.Deng Q. D., Zhang Y. M., Huan W. L., et al., 1980. Principles and methods of composing the seismic zoning map of China. Acta Seismologica Sinica, 2(1): 90—110. (in Chinese) 地震调查研究推进本部地震委员会, 2018. 全国地震动预测地图2018年版. 东京: 地震调查研究推进本部事务局. 高孟潭, 2003. 新的国家地震区划图. 地震学报, 25(6): 630—636. doi: 10.3321/j.issn:0253-3782.2003.06.009Gao M. T., 2003. New national seismic zoning map of China. Acta Seismologica Sinica, 25(6): 630—636. (in Chinese) doi: 10.3321/j.issn:0253-3782.2003.06.009 李善邦, 1957. 中国地震区域划分圖及其說明Ⅰ. 总的說明. 地球物理学报, 6(2): 127—158.Li S. B., 1957. The map of seismicity of China. Chinese Journal of Geophysics, 6(2): 127—158. (in Chinese) 李小军, 胡进军, 陈苏等, 2020. 海域地震动传播规律及场地影响模型研究报告, 16—28. 北京: 中国地震局地球物理研究所. 吕悦军, 冉洪流, 徐伟进等, 2020. 海域地震活动性模型及其评价方法研究报告, 23—31. 北京: 应急管理部国家自然灾害防治研究院. 潘华, 高孟潭, 谢富仁, 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 任治坤, 高战武, 李正芳等, 2020. 海域活动构造框架和地震构造模型研究报告, 18-24. 北京: 中国地震局地质研究所. 陶夏新, 1992. 我国新的地震区划编图和中国地震烈度区划图(1990). 自然灾害学报, 1(1): 99—109.Tao X. X., 1992. On the new national seismic zonation and the seismic intensity zoning map of China (1990). Journal of Natural Disasters, 1(1): 99—109. (in Chinese) 肖亮, 2019. 海域地震动预测方程模型研究报告, 8-12. 北京: 中国地震局地球物理研究所. 肖亮, 2020. 海域地震动传播规律及场地影响模型研究报告, 1—13. 北京: 中国地震局地球物理研究所. 俞言祥, 李山有, 肖亮, 2013. 为新区划图编制所建立的地震动衰减关系. 震灾防御技术, 8(1): 24—33. doi: 10.3969/j.issn.1673-5722.2013.01.003Yu Y. X., Li S. Y., Xiao L., 2013. Development of ground motion attenuation relations for the new seismic hazard map of China. Technology for Earthquake Disaster Prevention, 8(1): 24—33. (in Chinese) doi: 10.3969/j.issn.1673-5722.2013.01.003 中国地震烈度区划图编委会, 1992. 中国地震烈度区划图(1990)及其说明. 中国地震, 8(4): 1—11.The Compiling Committee of Seismic Zoning Map in China, 1992. Seismic intensity zoning map of China (1990) and its explanations. Earthquake Research in China, 8(4): 1—11. (in Chinese) 中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会, 2016. GB 18306—2015 中国地震动参数区划图. 北京: 中国标准出版社.General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China, Standardization Administration, 2016. GB 18306-2015 Seismic ground motion parameters zonation map of China. Beijing: Standards Press of China. (in Chinese) 周本刚, 陈国星, 高战武, 等, 2013. 新地震区划图潜在震源区划分的主要技术特色. 震灾防御技术, 8(2): 113—124. doi: 10.3969/j.issn.1673-5722.2013.02.001Zhou B. G., Chen G. X., Gao Z. W., et al., 2013. The technical highlights in identifying the potential seismic sources for the update of national seismic zoning map of China. Technology for Earthquake Disaster Prevention, 8(2): 113—124. (in Chinese) doi: 10.3969/j.issn.1673-5722.2013.02.001 Cheng J., Rong Y. F., Magistrale H., et al., 2017. An mw-based historical earthquake catalog for mainland China. Bulletin of the Seismological Society of America, 107(5): 2490—2500. doi: 10.1785/0120170102 Cornell C. A., 1968. Engineering seismic risk analysis. Bulletin of the Seismological Society of America, 58(5): 1583—1606. Danciu L., Şeşetyan K., Demircioglu M., et al., 2018. The 2014 earthquake model of the middle east: seismogenic sources. Bulletin of Earthquake Engineering, 16(8): 3465—3496. doi: 10.1007/s10518-017-0096-8 Giardini D., Wössner J., Danciu L., 2014. Mapping Europe’s seismic hazard. EOS, Transactions, American Geophysical Union, 95(29): 261—262. Gutenberg B., Richter C. F., 1944. Frequency of earthquakes in California. Bulletin of the Seismological Society of America, 34(4): 185—188. McGuire R. K., 1995. Probabilistic seismic hazard analysis and design earthquakes: closing the loop. Bulletin of the Seismological Society of America, 85(5): 1275—1284. Petersen M. D., Moschetti M. P., Powers P. M., et al., 2015. The 2014 United States national seismic hazard model. Earthquake Spectra, 31(1S): S1—S30. Wells D. L., Coppersmith K. J., 1994. New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement. Bulletin of the Seismological Society of America, 84(4): 974—1002. Woessner J., Laurentiu D., Giardini D., et al., 2015. The 2013 European seismic hazard model: key components and results. Bulletin of Earthquake Engineering, 13(12): 3553—3596. doi: 10.1007/s10518-015-9795-1 -