江苏江都震群震源区速度结构及构造特征

doi:10.3969/j.issn.0253-4967.2025.05.20240114

摘要/Abstract

摘要：

江苏省江都地区地处下扬子板块的苏中盆地, 是在中古生界复杂构造基础上演变形成的中新生代陆相沉积盆地, 区内发育大量演化的断裂、隆起及坳陷构造。文中基于深度学习的PhaseNet检测技术, 对2023年1月1日-2024年7月31日江都震群及附近地区进行地震重检测和精定位, 并通过TomoDD双差层析成像方法反演震源区的三维速度结构。基于地震精定位结果和反演的三维速度结构特征, 结合震源区M_L3.0以上地震的震源机制解, 对江都震群的孕震环境和发震机制进行了深入讨论。研究结果表明: 1)江都震群发震断裂为2条未知的左旋走滑隐伏断裂。其中, 江都震群1序列的发震断裂沿NNW走向, 断层面倾向为NE向; 江都震群2序列的发震断裂沿NNE向, 断层面略向NW倾。2条发震断裂近直立, 破裂均未延伸至地表。2)震源区附近波速变化剧烈, 陈家堡-小海断裂对震源区存在明显的分割作用, 断裂的北西侧存在明显的低波速、低泊松比异常。随着深度的增加, 断裂南东侧的高速异常逐渐增强。3)江都震群序列整体位于V_P、 V_S和泊松比低值异常体内, 且P波速度比S波速度下降更快, 整个发震过程没有地下流体的渗入与参与。基于发震断裂的性质及震源区速度结构特征, 分析认为江都震群序列的2次发震断裂或为陈家堡-小海断裂的2条分支断裂, 一条呈SSE-NNW向展布的左旋走滑隐伏构造, 一条呈SSW-NNE向展布的左旋走滑兼逆冲隐伏构造。江都震群是震源区区域应力的2次集中释放过程, 其直接成因是由震源区坚硬岩层发生脆性破裂。

关键词: 江都震群, 深度学习, 精定位, TomoDD, 速度结构, 构造特征

Abstract:

Since April 23, 2023, multiple earthquakes have occurred in Jiangdu, Jiangsu Province. On April 27, 2023, at 09:39 local time, an M3.1 event struck the region, followed by a series of seismic activities that constituted the Jiangdu earthquake swarm(hereafter referred to as “Jiangdu swarm 1”). Seismicity gradually diminished by June 22, 2023. On May 28, 2024, renewed seismic activity was observed in the same area, forming another swarm(hereafter “Jiangdu swarm 2”). On July 8, 2024, at 16:07 local time, an M3.6 earthquake occurred, after which activity again subsided by July 16, 2024. To investigate the velocity structure, seismotectonic setting, and possible relationship between these two swarms, both sequences were analyzed collectively(hereafter referred to as the Jiangdu earthquake swarm).

In this study, PhaseNet, a deep learning-based phase detection method, was employed to detect earthquakes in the epicentral area. The HypoDD algorithm was then used for precise relocation, producing a high-resolution catalog of the Jiangdu swarm. Additionally, seismic reports from January 2009 to July 2024 covering Jiangsu and adjacent provinces were compiled. Using the TomoDD double-difference tomography method, we inverted the three-dimensional velocity structure of V_P, V_S, and Poisson's ratio in the epicentral area. To constrain the seismogenic fault properties, focal mechanism solutions for seven M_L≥3.0 earthquakes were obtained with the HASH algorithm. Integrating precise locations, 3D velocity structures, and focal mechanisms, we identified the seismogenic faults of the Jiangdu swarm and analyzed its seismotectonic environment.

The results show that earthquakes in the Jiangdu swarms exhibit two predominant alignments, trending NNW and NNE, with focal depths concentrated between 7~16km. For Jiangdu swarm 1, focal mechanisms of four earthquakes indicate a NW-striking plane I, consistent with the NNW alignment of the relocated sequence. This plane is interpreted as the causative fault, which is a left-lateral strike-slip structure with minor normal faulting. For Jiangdu swarm 2, focal mechanisms of three earthquakes reveal a NE-striking plane I, consistent with the NNE alignment, and interpreted as a left-lateral strike-slip fault with a minor reverse component. Overall, the seven focal mechanism solutions show good agreement with the relocation results, indicating predominantly sinistral strike-slip motion.

Near the epicentral area, remarkable velocity contrasts are observed, with the Chenjiapu-Xiaohai Fault exerting a significant segmentation effect. The NW side of the fault is marked by low velocity and low Poisson's ratio anomalies, while the SE side displays increasing high-velocity anomalies with depth. Strong stratification of velocity and Poisson's ratio is also evident. The Jiangdu swarm is situated in a low-V_P, low-V_S, and low-Poisson's ratio anomaly zone, where the drop in P-wave velocity is more pronounced than in S-wave velocity, suggesting no involvement of fluids during the sequence. The low Poisson's ratio and narrow fault zone indicate that brittle fracture of rock strata was the dominant mechanism.

The Jiangdu source region is rich in shale oil and gas. The abundant shale gas is adsorbed in the pores and fractures of the rock formation. Long-term extraction and hydraulic fracturing enlarge rock fractures, reducing the effective elastic modulus and lowering seismic wave velocities. Based on precise relocation, seismogenic fault geometry, and crustal velocity structures, we infer that the two seismogenic faults of the Jiangdu swarm are likely subsidiary branches of the Chenjiapu-Xiaohai Fault. One is a concealed left-lateral strike-slip fault trending SSE-NNW, and the other is a concealed left-lateral strike-slip with thrust component trending SSW-NNE. The brittle failure of the hard rock strata directly triggered the Jiangdu swarms, representing two concentrated episodes of stress release. These findings provide new insights into the seismogenic environment and mechanisms of earthquake swarms in the Jiangdu region.

Key words: Jiangdu earthquake swarm, deep-learning, accurate earthquake location, TomoDD, velocity structure, seismogenic structure

李婷婷, 缪发军, 孙业君, 范文华, 宫杰, 顾勤平, 杜航, 孙小航, 张岑, 李子烨. 江苏江都震群震源区速度结构及构造特征[J]. 地震地质, 2025, 47(5): 1343-1363.

LI Ting-ting, MIAO Fa-jun, SUN Ye-jun, FAN Wen-hua, GONG Jie, GU Qin-ping, DU Hang, SUN Xiao-hang, ZHANG Cen, LI Zi-ye. THREE-DIMENSIONAL WAVE VELOCITY STRUCTURE AND SEISMOGENIC STRUCTURE FOR THE JIANGDU EARTHQUAKE SWARM IN JIANGSU[J]. SEISMOLOGY AND GEOLOGY, 2025, 47(5): 1343-1363.

图/表 13

图 1 a 苏中盆地地震构造图; b 江都震群震中及附近断裂分布图图c中红色方框为图b的位置, 图a中黄色五角星为江都震源区的位置, 图b中黄色五角星分别为江都震群2次主震的位置, 蓝色实线为4条测线的位置

Fig. 1 Seismic structure map of Suzhong Basin(a), distribution of faults and historical earthquakes in the epicenter region of Jiangdu earthquake(b).

图 2 P、 S震相时距曲线

Fig. 2 Phase time distance curves of P-wave and S-wave.

图 3 VP(a)、 VS(b)层析成像地震事件、台站及地震射线分布

Fig. 3 The hypocentres, stations and rays in the horizontal plane for VP and VS tomography.

表 1 初始一维速度模型

Table1 Initial 1-D velocity model

序号	深度/km	V_P/km·s^-1	V_S/km·s^-1
1	0	3.792	2.239
2	5	5.121	3.025
3	10	6.062	3.581
4	15	6.212	3.669
5	20	6.451	3.810
6	25	6.500	3.759
7	30	6.824	3.972
8	34	7.400	4.299
9	40	8.100	4.6900

图 4 平滑因子和阻尼因子的L-curve曲线

Fig. 4 L-curves of the smoothing and damping factors.

图 5 P、 S震相走时残差直方图

Fig. 5 Histograms of arrival time residuals for P-phase and S-phase.

图 6 0km、 5km、 10km、 15km、 20km不同深度的P波、 S波检测板测试结果红色五角星表示江都震群位置

Fig. 6 P-wave and S-wave checkerboard resolution at depths of 0km、 5km、 10km、 15km and 20km.

图 7 不同目录震级-频次对比图

Fig. 7 Comparison of magnitude-frequency in different catalogs.

图 8 人工地震目录和检测地震目录的震中分布、 M-T图 a 为人工目录地震分布; b 检测地震精定位后的地震分布

Fig. 8 M-T plot and epicenter distribution of earthquakes from manually and automatically detected catalogs.

图 9 不同深度P波、 S波速度和泊松比分布图

Fig. 9 Tomography results of P-wave and S-Wave velocity and poisson's ratio at different depth.

图 10 地震精定位深度剖面图及ML3.0以上地震的震源机制解

Fig. 10 Depth profiles of seismic precise positioning distribution and focal mechanism solutions of selected earthquakes with magnitude ML≥3.0.

表 2 ML3.0地震的震源机制解

Table2 Focal mechanism solutions of selected earthquakes with magnitude ML≥3.0

序号	节面Ⅰ			节面Ⅱ			发震时刻	震级/M_L
序号	走向 /(°)	倾角 /(°)	滑动角 /(°)	走向 /(°)	倾角 /(°)	滑动角 /(°)	发震时刻	震级/M_L
1	305	75	-10	33	78	-136	2023-04-27, 09:39:59	3.8
2	312	71	-23	42	75	-125	2023-05-10, 23:20:43	3.5
3	298	68	-26	205	80	-145	2023-06-19, 02:15:24	3.4
4	293	75	-11	199	77	-138	2023-07-10, 14:17:43	3.4
5	234	73	39	323	72	135	2024-07-06, 21:31:33	3.4
6	226	76	37	318	63	126	2024-07-08, 16:07:19	4.1
7	246	69	29	321	64	140	2024-07-11, 05:22:18	3.6

图 11 江都震群震源区附近P波、 S波速度和泊松比剖面图

Fig. 11 Poisson's ratio profile and P-wave and S-wave velocity near the source area of the Jiangdu earthquake swarm.

参考文献 36

[1]	曹颖, 钱佳威, 黄江培, 等. 2024. 云南洱源地区地壳三维精细速度结构成像[J]. 地震地质, 46(1): 162-183. doi: 10.3969/j.issn.0253-4967.2024.01.010.
	CAO Ying, QIAN Jia-wei, HUANG Jiang-pei, et al. 2024. Crustal fine velocity structure in the Eryuan area, Yunnan from double-different tomography[J]. Seismology and Geology, 46(1): 162-183 (in Chinese).
[2]	顾勤平, 丁志峰, 康清清, 等. 2020a. 郯庐断裂带中南段及邻区基于背景噪声的瑞利波群速度层析成像[J]. 地球物理学报, 63(4): 1505-1522.
	GU Qin-ping, DING Zhi-feng, KANG Qing-qing, et al. 2020a. Group velocity tomography of Rayleigh wave in the middle-southern of the Tan-Lu fault zone and adjacent regions using ambient seismic noise[J]. Chinese Journal of Geophysics, 63(4): 1505-1522 (in Chinese).
[3]	顾勤平, 康清清, 张鹏, 等. 2020b. 郯庐断裂带中南段及邻区Rayleigh波相速度与方位各向异性[J]. 地震地质, 42(5): 1129-1152. doi: 10.3969/j.issn.0253-4967.2020.05.007.
	GU Qin-ping, KANG Qing-qing, ZHANG Peng, et al. 2020b. Rayleigh wave phase velocity and azimuthal anisotropy of the middle-southern segment of the Tan-Lu fault zone and adjacent regions from ambient noise tomography[J]. Seismology and Geology, 42(5): 1129-1152 (in Chinese).
[4]	黄耘, 李清河, 张元生, 等. 2011. 郯庐断裂带鲁苏皖段及邻区地壳速度结构[J]. 地球物理学报, 54(10): 2549-2559.
	HUANG Yun, LI Qing-he, ZHANG Yuan-sheng, et al. 2011. Crustal velocity structure beneath the Shangdong-Jiangsu-Anhui segment of the Tancheng-Lujiang Fault zone and adjacent areas[J]. Chinese Journal of Geophysics, 54(10): 2549-2559 (in Chinese).
[5]	康清清, 缪发军, 刘红桂, 等. 2015. 利用矩张量反演法研究江苏高邮-宝应 M_S9 地震震源机制解和震源深度[J]. 地球物理学报, 58(1): 204-215.
	KANG Qing-qing, MIAO Fa-jun, LIU Hong-gui, et al. 2015. Focal mechanism and focal depth of the Gaoyou-Baoying, Jiangsu, China M_S4.9 earthquake on 20 July 2012 derived from moment tensor inversion[J]. Chinese Journal of Geophysics, 58(1): 204-215 (in Chinese).
[6]	李翠平, 唐茂云, 郭卫英, 等. 2022. 荣昌及周边三维速度结构初步研究[J]. 地震地质, 44(1): 205-219. doi: 10.3969/j.issn.0253-4967.2022.01.013.
	LI Cui-ping, TANG Mao-yun, GUO Wei-ying, et al. 2022. A preliminary study on 3D velocity structure of Rongchang and its adjacent area[J]. Seismology and Geology, 44(1): 205-219 (in Chinese).
[7]	李婷婷, 刘利, 范文华, 等. 2020. 基于远震P波接收函数研究江苏地区地壳厚度和泊松比[J]. 地震研究, 43(4): 680-688, 768.
	LI Ting-ting, LIU Li, FAN Wen-hua, et al. 2020. Study on crustal thickness and ratio of Poisson in Jiangsu Province by teleseismic receiver function[J]. Journal of Seismological Research, 43(4): 680-688, 768 (in Chinese).
[8]	李婷婷, 张岑, 范文华, 等. 2024. 基于接收函数的郯庐断裂带南段及其周缘壳内结构研究[J]. 防灾减灾工程学报, 44(6): 1465-1473.
	LI Ting-ting, ZHANG Cen, FAN Wen-hua, et al. 2020. Crustal structure of the southern section and surroundings of the Tan-Lu fault zone based on receiver functions[J]. Journal of Disaster Prevention and Mitigation Engineering, 44(6): 1465-1473 (in Chinese).
[9]	廖诗荣, 杨婷, 张红才, 等. 2021. 2021年云南双柏 M_S5.1 地震序列的快速检测与精定位研究[J]. 地震研究, 44(4): 515-520.
	LIAO Shi-rong, YANG Ting, ZHANG Hong-cai, et al. 2021. Fast detection and relocation of the 2021 Shuangbai M_S5.1 earthquake sequences in Yunnan[J]. Journal of Seismological Research, 44(4): 515-520 (in Chinese).
[10]	刘亢, 杨婷, 李红光, 等. 2023. 小地震精定位与层析成像揭示的邢台地震区深部构造特征[J]. 地震地质, 45(6): 1328-1348. doi: 10.3969/j.issn.0253-4967.2023.06.005.
	LIU Kang, YANG Ting, LI Hong-guang, et al. 2023. Deep structural characteristics in the Xingtai seismic zone revealed by precise relocation of small earthquakes and tomography[J]. Seismology and Geology, 45(6): 1328-1348 (in Chinese).
[11]	刘瑞丰, 陈运泰, 薛峰, 等. 2018. 如何正确使用新的震级国家标准[J]. 地震地磁观测与研究, 39(4): 1-11.
	LIU Rui-feng, CHEN Yun-tai, XUE Feng, et al. 2018. How to use the new national standard of magnitude correctly[J]. Seismological and Geomagnetic Observation and Research, 39(4): 1-11 (in Chinese).
[12]	钱基. 2000. 苏北盆地油气田的形成与分布特征[J]. 石油大学学报(自然科学版), 24(4): 21-25, 124.
	QIAN Ji. 2000. Formation and distribution of oil and gas fields in Subei Basin[J]. Journal of China University of Petroleum(Edition of Natural Science), 24(4): 21-25, 124 (in Chinese).
[13]	邱旭明, 陈伟, 李鹤永, 等. 2023. 苏北盆地走滑构造与复杂断块油气成藏[J]. 石油实验地质, 45(3): 393-401.
	QIU Xu-ming, CHEN Wei, LI He-yong, et al. 2023. Strike-slip structures and hydrocarbon accumulation in complex fault blocks in Subei Basin[J]. Petroleum Geology & Experiment, 45(3): 393-401 (in Chinese).
[14]	万永革, 沈正康, 刁桂苓, 等. 2008. 利用小震分布和区域应力场确定大震断层面参数方法及其在唐山地震序列中的应用[J]. 地球物理学报, 51(3): 793-804.
	WAN Yong-ge, SHEN Zheng-kang, DIAO Gui-ling, et al. 2008. An algorithm of fault parameter determination using distribution of small earthquakes and parameters of regional stress field and its application to Tangshan earthquake sequence[J]. Chinese Journal of Geophysics, 51(3): 793-804 (in Chinese).
[15]	王俊, 缪发军, 詹小艳, 等. 2013. 利用sPn震相精确测定江苏高邮、宝应交界M4.9地震震源深度[J]. 防灾减灾工程学报, 33(3): 354-360.
	WANG Jun, MIAO Fa-jun, ZHAN Xiao-yan, et al. 2013. Precise determination of the focal depth of Jiangsu Gaoyou-Baoying M4.9 earthquake by sPn phase[J]. Journal of Disaster Prevention and Mitigation Engineering, 33(3): 354-360 (in Chinese).
[16]	王小娜, 邓志辉, 叶秀薇, 等. 2020. 广东阳江地区的地壳速度结构与地震活动性[J]. 地震地质, 42(5): 1153-1171. doi: 10.3969/j.issn.0253-4967.2020.05.008.
	WANG Xiao-na, DENG Zhi-hui, YE Xiu-wei, et al. 2020. The study of crustal velocity structure and seismicity in Yangjiang area of Guangdong Province[J]. Seismology and Geology, 42(5): 1153-1171 (in Chinese).
[17]	吴林, 陈清华, 刘寅, 等. 2017. 裂陷盆地伸展方位与构造作用及对构造样式的控制: 以苏北盆地高邮凹陷南部断阶带为例[J]. 石油与天然气地质, 38(1): 29-38.
	WU Lin, CHEN Qing-hua, LIU Yin, et al. 2017. Extension direction and tectonism in rift basins and their control on structural style: A case study on the south fault terrace in the Gaoyou Sag, Subei Basin[J]. Oil & Gas Geology, 38(1): 29-38 (in Chinese).
[18]	尹欣欣, 蒋长胜, 蔡润, 等. 2021. 云南漾濞地区地壳层析成像与地震精定位[J]. 地震地质, 43(4): 864-880. doi: 10.3969/j.issn.0253-4967.2021.04.008.
	YIN Xin-xin, JIANG Chang-sheng, CAI Run, et al. 2021. Study of crustal tomography and precise earthquake location in Yangbi area, Yunan Province[J]. Seismology and Geology, 43(4): 864-880 (in Chinese).
[19]	于悦颖, 李正楷, 杨云, 等. 2024. 2022年东台地震序列区域速度结构和精确定位[J]. 地震地质, 46(3): 627-644. doi: 10.3969/j.issn.0253-4967.2024.03.007.
	YU Yue-ying, LI Zheng-kai, YANG Yun, et al. 2024. Regional velocity structure and relocation of the 2022 Dongtai earthquake sequence[J]. Seismology and Geology, 46(3): 627-644 (in Chinese).
[20]	张大莲, 张平, 李向前, 等. 2021. 江苏沿海地区小波多尺度重力场特征与构造活动性分析[J]. 大地测量与地球动力学, 41(2): 186-191.
	ZHANG Da-lian, ZHANG Ping, LI Xiang-qian, et al. 2021. Characteristics of wavelet multi-scale gravity field and analysis of tectonic activity in Jiangsu coastal area[J]. Journal of Geodesy and Geodynamics, 41(2): 186-191 (in Chinese).
[21]	Cai G Y, Wang W L, Wu J P, et al. 2023. Three-dimensional body wave velocity structure and seismogenic structure for the 2022 M_S6.8 Luding earthquake sequence in Sichuan, China[J]. Frontiers in Earth Science, 11:1099744. DOI URL
[22]	Fang L H, Wu J P, Wang W L, et al. 2013. Relocation of the mainshock and aftershock sequences of M_S7.0 Sichuan Lushan earthquake[J]. Chinese Science Bulletin, 58(28): 3451-3459.
[23]	Han S C, Zhang H J, Xin H L, et al. 2022. USTClitho2. 0: Updated unified seismic tomography models for continental China lithosphere from joint inversion of body-wave arrival times and surface-wave dispersion data[J]. Seismological Research Letters, 93(1): 201-215.
[24]	Hansen P C, O’Leary D P. 1993. The use of the L-curve in the regularization of discrete ill-posed problems[J]. SIAM Journal on Scientific Computing, 14(6): 1487-1503.
[25]	Hardebeck J L, Shearer P M. 2003. Using S/P amplitude ratios to constrain the focal mechanisms of small earthquakes[J]. Bulletin of the Seismological Society of America, 93: 2434-2444. DOI URL
[26]	Klein F W. 1989. User's guide to HYPOINVERSE, a program for VAX computers to solve for earthquake locations and magnitudes[M]. US Geological Survey, Reston, VA.
[27]	Li Z, Zhou L, Duan M, et al. 2024. Three-dimensional V_P, V_S, and V_P/V_S imaging based on AI microseismic detection reveals the mechanism of induced earthquakes in the Xiluodu Reservoir Area, China[J]. Journal of Asian Earth Sciences, 266:106123. DOI URL
[28]	Thurber C H. 1992. Hypocenter-velocity structure coupling in local earthquake tomography[J]. Physics of the Earth and Planetary Interiors, 75(1-3): 55-62. DOI URL
[29]	Thurber C H, Zhang H, Brocher T M, et al. 2009. Regional three-dimensional seismic velocity model of the crust and uppermost mantle of northern California[J]. Journal of Geophysical Research: Solid Earth, 114(B1): B01304.
[30]	Wan Y G. 2010. Contemporary tectonic stress field in China[J]. Earthquake Science, 23(4): 377-386.
[31]	Xin H L, Zhang H J, Kang M, et al. 2018. High-lithospheric velocity structure of continental China by double-difference seismic travel-time tomography[J]. Seismological Research Letters, 90(1): 229-241.
[32]	Zhang H J, Thurber C H. 2003. Double-difference tomography: the method and its application to the Hayward Fault, California[J]. Bulletin of the Seismological Society of America, 93(5): 1875-1889.
[33]	Zhang H J, Thurber C H. 2006. Development and applications of double-difference seismic tomography[J]. Pure and Applied Geophysics, 163(2): 373-403.
[34]	Zhang M, Ellsworth W L, Beroza G C. 2019. Rapid earthquake association and location[J]. Seismological Research Letters, 90(6): 2276-2284.
[35]	Zhang M, Liu M, Feng T, et al. 2022. LOC-FLOW: An end-to-end machine learning-based high-precision earthquake location workflow[J]. Seismological Research Letters, 93(5): 2426-2438.
[36]	Zhu W, Beroza G C. 2019. PhaseNet: A deep-neural-network-based seismic arrival-time picking method[J]. Geophysical Journal International, 216(1): 261-273.