沂沭断裂带及其邻区地壳P波层析成像研究

doi:10.13745/j.esf.sf.2022.2.62

地学前缘 ›› 2023, Vol. 30 ›› Issue (2): 272-281.DOI: 10.13745/j.esf.sf.2022.2.62

• 构造-岩浆作用与大地构造背景 • 上一篇下一篇

沂沭断裂带及其邻区地壳P波层析成像研究

崔鑫¹^,²(), 宋秀青³^,⁴^,^*(), 石耀霖⁵, 殷海涛¹, 朱元清³^,⁴, 董腾超¹

1.山东省地震局, 山东济南 250014
2.中国地震局地球物理研究所, 北京 100081
3.上海市地震局, 上海 200062
4.上海佘山地球物理国家野外科学观测研究站, 上海 201602
5.中国科学院大学, 北京 100049

收稿日期:2021-07-13 修回日期:2021-12-16 出版日期:2023-03-25 发布日期:2023-01-05
通讯作者: 宋秀青
作者简介:崔鑫(1983—),男,博士研究生,主要从事区域精细构造以及层析成像研究工作。E-mail: xc_0071@163.com
基金资助:
国家自然科学基金-中国地震局地震科学联合基金项目(U1839207);国家自然科学基金项目(41974105);上海市科委科研计划项目(19dz1200200);上海佘山地球物理国家野外科学观测研究站基金项目(2020z02)

Crustal P-wave tomographic study of the Yishu fault zone

CUI Xin¹^,²(), SONG Xiuqing³^,⁴^,^*(), SHI Yaolin⁵, YIN Haitao¹, ZHU Yuanqing³^,⁴, DONG Tengchao¹

1. Shandong Earthquake Agency, Jinan 250014, China
2. Institute of Geophisics, China Earthquake Administration, Beijing 100081, China
3. Shanghai Earthquake Agency, Shanghai 200062, China
4. Shanghai Sheshan Geophisical National Geophysical Obeservatory, Shanghai 201602, China
5. University of Chinese Academy of Sciences, Beijing 100049, China

Received:2021-07-13 Revised:2021-12-16 Online:2023-03-25 Published:2023-01-05
Contact: SONG Xiuqing

摘要/Abstract

摘要：

沂沭断裂带是华北克拉通东部的重要断裂带,历史上发生过8¹/₂级强震,区域速度模型研究可以提升地震震动图精度以及地震烈度速报的准确性,有助于防御未来的巨震,降低震害损失。本文搜集了沂沭断裂带及其邻近区域2008—2020年的数字地震波数据,经过数据降丛集化、走时数据正态分布筛选等处理得到了地震事件震中参数和走时震相数据。基于FMTOMO方法,反演得到了沂沭断裂带地壳三维速度结构。反演结果显示,在近地表地层,速度分布特征与地表地质构造单元存在密切联系,褶皱带和基岩出露区域常常呈现高速速度分布,含较厚土层的凹陷或海域区域则一般呈现低速速度分布;在中地壳地层,速度分布特征与地震活动性有密切的联系,地震震中往往密集分布于高速和低速的过渡区域。在莫霍面深度,沂沭带的速度分布呈现了分段性差异特征,且两侧的速度分布存在差异。这可能意味着断裂带深达地幔的现象。且在沂沭带中南段区域,存在高速速度分布特征,这可能与区域地幔物质上涌有关,也可能构成强震的孕震背景。郯城区域从地表至地幔顶部均显示高速速度异常,该区域会积累高应力,可能发生8级以上地震,推测其复发周期应该很长。

关键词: 山东, 沂沭带, P波, 层析成像

Abstract:

The Yishu Fault zone is an important fault zone in the eastern North China Craton and historically experienced a magnitude 8¹/₂ earthquake. To protect against future strong earthquakes and reduce earthquake damage, it is helpful to study the regional seismic velocity model so as to improve the accuracy of seismogram and rapid reporting of seismic intensity. In this study we collect digital seismic data recorded between 2008-2020 in the Yishu Fault zone, and preprocess the data using declustering and normal distribution travel-time selection techniques to obtain the local epicenter parameters and travel-time data; then by tomographic inversion using FMTOMO software we obtained three-dimensional seismic velocity model of the Yishu Fault zone. According to the 3D model, the velocity distribution features in the near surface layer are closely related to the geological structural units, as the fold belt and bedrock outcropping areas often present high velocity distribution while low velocity distribution is generally associated with the sag or sea area with thick soil layers. In the middle crustal layer, seismicity is closely related to the velocity distribution features, and the earthquake epicenters are often densely distributed in the transition zone between high and low velocities. At the Moho depth, velocity distribution in the Yishu belt shows segmental differences and differs on both sides of the belt, this could mean that the fault zone reaches the mantle. High velocity distribution is present in the middle and southern part of the Yishu belt, which may be related to the regional upwelling of mantle material, or may constitute the seismogenic background of strong earthquakes. In the Tancheng area, high velocity anomalies are observed from the earth’s surface to the top of the mantle, suggesting that high stress could accumulate in this area, and an earthquake of magnitude 8 or above might occur, but the recurrence period is considered very long.

Key words: Shandong, Yishu Fault, P-wave, tomography

中图分类号:

P315.2

崔鑫, 宋秀青, 石耀霖, 殷海涛, 朱元清, 董腾超. 沂沭断裂带及其邻区地壳P波层析成像研究[J]. 地学前缘, 2023, 30(2): 272-281.

CUI Xin, SONG Xiuqing, SHI Yaolin, YIN Haitao, ZHU Yuanqing, DONG Tengchao. Crustal P-wave tomographic study of the Yishu fault zone[J]. Earth Science Frontiers, 2023, 30(2): 272-281.

图/表 8

图1 研究区构造背景和台站分布

Fig.1 Tectonic setting and distribution of observation stations in the study area

图2 地震震中分布图

Fig.2 Earthquake epicenter distribution map

图3 地震射线数量与角度分布

Fig.3 (a) Quantity and (b-i) angle distributions of seismic rays

图4 检测板实验结果

Fig.4 Results of checkerboard tests

图5 地震一维速度模型

Fig.5 Seismic one-dimensional velocity model

图6 模型误差和数据误差的折衷曲线

Fig.6 Trade-off L-curve between model and travel-time residual variances

图7 沂沭带三维P波速度扰动分布(水平剖面)

Fig.7 Three-dimensional P-wave velocity distribution around the Yishu Fault zone (horizontal slice)

图8 沂沭带三维P波速度扰动分布(垂直剖面)

Fig.8 Three-dimensional P-wave velocity distribution around the Yishu Fault zone (vertical slice)

参考文献 45

[1]	朱光, 王道轩, 刘国生, 等. 郯庐断裂带的演化及其对西太平洋板块运动的响应[J]. 地质科学, 2004, 39(1): 36-49.
[2]	高维明, 郑朗荪, 李家灵, 等. 1668年郯城8.5级地震的发震构造[J]. 中国地震, 1988, 4(3): 15-21.
[3]	陈兆辉, 王椿镛, 楼海. 鄂尔多斯地块地壳上地幔速度结构及构造意义[J]. 科学通报, 2018, 63(3): 327-339.
[4]	朱元清, 宋秀青, 刘双庆, 等. 中国地震测定参考速度结构研究[M]. 北京: 地震出版社, 2017.
[5]	崔鑫, 李铂, 朱元清, 等. 基于Hyposat和HypoDD的山东地区地震波一维速度模型[J]. 吉林大学学报(地球科学版), 2017, 47(5): 1572-1581.
[6]	顾勤平, 丁志峰, 康清清, 等. 郯庐断裂带中南段及邻区Pn波速度结构与各向异性[J]. 地球物理学报, 2016, 59(2): 504-515.
[7]	殷伟伟, 雷建设, 杜沫霏, 等. 郯庐断裂带及其邻区上地幔顶部Pn波速度与各向异性层析成像[J]. 地球物理学报, 2019, 62(11): 4227-4238.
[8]	黄耘, 李清河, 张元生, 等. 郯庐断裂带鲁苏皖段及邻区地壳速度结构[J]. 地球物理学报, 2011, 54(10): 2549-2559.
[9]	熊振, 李清河, 张元生, 等. 郯庐断裂带鲁苏皖段地壳速度结构的分段特征及其地质意义[J]. 地球物理学报, 2016, 59(7): 2433-2443.
[10]	LEI J S, ZHAO D P, XU X W, et al. P-wave upper-mantle tomography of the Tanlu fault zone in eastern China[J]. Physics of the Earth and Planetary Interiors, 2020, 299: 106402. DOI URL
[11]	段永红, 张先康, 方盛明. 华北地区上部地壳结构的三维有限差分层析成像[J]. 地球物理学报, 2002, 45(3): 362-369.
[12]	李文辉, 高锐, KELLER R, 等. 华北克拉通北缘(怀来—苏尼特右旗)地壳结构[J]. 地球物理学报, 2014, 57(2): 472-483.
[13]	郭震, 唐有彩, 陈永顺, 等. 华北克拉通东部地壳和上地幔结构的接收函数研究[J]. 地球物理学报, 2012, 55(11): 3591-3600.
[14]	刘保金, 酆少英, 姬计法, 等. 郯庐断裂带中南段的岩石圈精细结构[J]. 地球物理学报, 2015, 58(5): 1610-1621.
[15]	吕子强, 郑建常, 张刚, 等. 濮阳地震集中区双差层析成像研究[J]. 地震研究, 2016, 39(2): 255-260.
[16]	毛慧慧, 雷建设, 滕吉文. 鄂尔多斯盆地北缘南北向剖面上地幔远震P波层析成像[J]. 地球物理学报, 2016, 59(6): 2056-2065.
[17]	凌媛, 陈凌, 危自根, 等. 利用噪声层析成像研究华北克拉通中-西部南段及邻区祁连造山带地壳S波各向异性速度结构[J]. 中国科学: 地球科学, 2017, 47(10): 1202-1219.
[18]	房立华, 吴建平, 吕作勇. 华北地区基于噪声的瑞利面波群速度层析成像[J]. 地球物理学报, 2009, 52(3): 663-671.
[19]	吴萍萍, 丁志峰, 马小军, 等. 基于背景噪声研究大别-苏鲁及其邻区的瑞雷波群速度结构[J]. 地震学报, 2015, 37(2): 218-229, 370.[知网]
[20]	何正勤, 叶太兰, 丁志峰. 华北东北部的面波相速度层析成像研究[J]. 地球物理学报, 2009, 52(5): 1233-1242.
[21]	王俊菲. 用远震接收函数研究下扬子地区地壳结构[J]. 国际地震动态, 2012, 42(6): 74.
[22]	张风雪, 李永华, 吴庆举, 等. FMTT方法研究华北及邻区上地幔P波速度结构[J]. 地球物理学报, 2011, 54(5): 1233-1242.
[23]	王志铄. 大华北及其邻区地壳上地幔三维速度结构的地震层析成像研究[D]. 北京: 中国地震局地球物理研究所, 2005.
[24]	翟明国, 彭澎. 华北克拉通古元古代构造事件[J]. 岩石学报, 2007, 23(11): 2665-2682.
[25]	丰成君, 张鹏, 戚帮申, 等. 郯庐断裂带附近地壳浅层现今构造应力场[J]. 现代地质, 2017, 31(1): 46-70.
[26]	杨云, 王维, 王俊菲. 郯庐断裂带中南段应力场及变形特征[J]. 地震工程学报, 2017, 39(增刊1): 14-19.
[27]	黄方, 何丽娟, 吴庆举. 鄂尔多斯盆地深部热结构特征及其对华北克拉通破坏的启示[J]. 地球物理学报, 2015, 58(10): 3671-3686.
[28]	唐新功, 陈永顺, 唐哲. 应用布格重力异常研究郯庐断裂构造[J]. 地震学报, 2006, 28(6): 603-610, 679.
[29]	叶高峰, 魏文博, 金胜, 等. 郯庐断裂带中段电性结构及其地学意义研究[J]. 地球物理学报, 2009, 52(11): 2818-2825.
[30]	李富珍, 任天翔, 池顺良, 等. 基于钻孔应变观测资料分析远震造成的动态库仑应力变化[J]. 地球物理学报, 2021, 64(6): 1949-1974.
[31]	李小林, 张冬丽, 张培震, 等. 多源探测数据揭示的南海东北部三维地壳结构特征[J]. 地球物理学报, 2020, 63(4): 1538-1552.
[32]	胡立天, 张培震, 郝天珧, 等. 南海西南次海盆岩石圈结构及其地质意义[J]. 地球物理学报, 2019, 62(9): 3464-3481.
[33]	顾庙元, 姚佳琪, 张伟, 等. 地学长江计划安徽实验中低频可控震源地震波信号提取方法评估[J]. 中国地震, 2016, 32(2): 356-378.
[34]	赵峦啸, 刘金水, 姚云霞, 等. 基于随机森林算法的陆相沉积烃源岩定量地震刻画: 以东海盆地长江坳陷为例[J]. 地球物理学报, 2021, 64(2): 700-715.
[35]	宋秀青, 孙冬军, 朱元清. 东海八角亭海洋深井台地震监测[J]. 地震地磁观测与研究, 2017, 38(3): 159-163.
[36]	RAWLINSON N, URVOY M. Simultaneous inversion of active and passive source datasets for 3-D seismic structure with application to Tasmania[J]. Geophysical Research Letters, 2006, 33(24): L24313. DOI URL
[37]	顾功叙. 中国地震目录: 公元前1831—公元1969年[M]. 北京: 科学出版社, 1983.
[38]	DE KOOL M, RAWLINSON N, SAMBRIDGE M. A practical grid-based method for tracking multiple refraction and reflection phases in three-dimensional heterogeneous media[J]. Geophysical Journal International, 2006, 167(1): 253-270. DOI URL
[39]	RAWLINSON N, DE KOOL M, SAMBRIDGE M. Seismic wavefront tracking in 3D heterogeneous media: applications with multiple data classes[J]. Exploration Geophysics, 2006, 37(4): 322-330. DOI URL
[40]	吴建平, 黄媛, 张天中, 等. 汶川M_s8.0级地震余震分布及周边区域P波三维速度结构研究[J]. 地球物理学报, 2009, 52(2): 320-328.
[41]	郭飚, 刘启元, 陈九辉, 等. 川西龙门山及邻区地壳上地幔远震P波层析成像[J]. 地球物理学报, 2009, 52(2): 346-355.
[42]	MEISSNER R, TANNERB. Crustal viscosities and seismic velocities[J]. Physics of the Earth and Planetary Interiors, 1992, 69(3/4): 252-256. DOI URL
[43]	SHINEVAR W J, BEHN M D, HIRTH G, et al. Inferring crustal viscosity from seismic velocity: application to the lower crust of southern California[J]. Earth and Planetary Science Letters, 2018, 494: 83-91. DOI URL
[44]	AWANG H, AHMAD RASHIDI N R, YUSOF M, et al. Correlation between P-wave velocity and strength index for shale to predict uniaxial compressive strength value[J]. MATEC Web of Conferences, 2017, 103: 07017. DOI URL
[45]	雷建设, 赵大鹏, 苏金蓉, 等. 龙门山断裂带地壳精细结构与汶川地震发震机理[J]. 地球物理学报, 2009, 52(2): 339-345.

沂沭断裂带及其邻区地壳P波层析成像研究

Crustal P-wave tomographic study of the Yishu fault zone

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 8

参考文献 45

相关文章 15

编辑推荐

Metrics

本文评价

[1]	高福建, 邹志辉, 王永红. 基于地震层析速度模型分析海滩地下结构:以青岛石老人海水浴场为例[J]. 地学前缘, 2022, 29(5): 275-284.
[2]	李小伟, 单伟, 于学峰, 李大鹏, 谢元惠, 张国坤, 迟乃杰, 王文鲁, 张岩, 李增胜, 马祥县. 五莲七宝山地区早白垩世碱性侵入岩成因及其地质意义[J]. 地学前缘, 2022, 29(5): 438-463.
[3]	李文辉, 王海燕, 高锐, 卢占武, 李洪强, 侯贺晟, 熊小松, 叶卓. 秦岭造山带及邻区上地壳精细速度结构研究[J]. 地学前缘, 2022, 29(2): 198-209.
[4]	吴国炜, 熊小松, 高锐, 陈宣华, 李英康, 王冠, 王小成, 任海东. 北山构造带南部上地壳二维初至波层析成像[J]. 地学前缘, 2022, 29(2): 402-415.
[5]	梁锋，吕庆田，严加永，刘振东. 宁芜盆地深部地壳结构和岩浆侵入体形态特征及其对成矿的启示[J]. 地学前缘, 2017, 24(5): 138-148.
[6]	韦伟,许建东,于红梅. 西昆仑阿什库勒火山的起源：来自地震层析成像的证据[J]. 地学前缘, 2015, 22(6): 227-232.
[7]	高宗军, 张福存, 安永会, 冯建国, 王敏, 韩克. 山东高密高氟地下水成因模式与原位驱氟设想[J]. 地学前缘, 2014, 21(4): 50-58.
[8]	徐峣，张贵宾，江国明，索奎，姚永波. 中国东北地区三维速度结构与火山起源的深部机制[J]. 地学前缘, 2014, 21(4): 255-265.
[9]	瞿辰, 杨文采, 于常青. 塔里木盆地地震波速扰动及泊松比成像[J]. 地学前缘, 2013, 20(5): 196-206.
[10]	李理，钟大赉，杨长春，时秀朋，胡秋媛，赵利，孙钰皓，刘卉. 伸展作用序列及其深部背景: 以晚中生代以来鲁西隆起和济阳坳陷为例[J]. 地学前缘, 2012, 19(5): 255-273.
[11]	杨文采，瞿辰，于常青. 华北东部地区地壳泊松比异常及其成因[J]. 地学前缘, 2011, 18(3): 13-21.
[12]	刘平华，刘福来，王舫，刘建辉. 山东半岛高压麻粒岩中锆石的U-Pb定年及其地质意义[J]. 地学前缘, 2011, 18(2): 33-54.
[13]	刘丽军, Michael Gurnis. 伴随数据同化方法及其在地幔对流中的应用[J]. 地学前缘, 2010, 17(3): 139-148.
[14]	韦伟, 赵大鹏, 石耀霖. 西南日本火山区三维地震成像研究[J]. 地学前缘, 2010, 17(3): 149-157.
[15]	董树义, 顾雪祥, 杨永强, 翟裕生, 章永梅, 向中林, 刘丽, 胡贵增, 李忠平. 地物化综合方法在危机矿山成矿预测中的应用——以山东沂南金矿为例[J]. 地学前缘, 2010, 17(2): 198-208.