Identification of Fault Activity in Urban Active Fault Exploration and Seismic Risk Assessment in Zhenjiang City
-
摘要: 镇江市是江苏省破坏性地震发生较多的地区。对镇江市主要发震断层开展活动性探测及地震危险性评价工作,查明其发震能力,对于镇江市防震减灾事业发展具有重要意义。对镇江市区威胁较大的主要隐伏断层有4条,分别为幕府山-焦山断裂、五峰山-西来桥断裂、丹徒-建山断裂、汝山-上会断裂。通过开展浅层人工地震勘探、钻孔联合地质剖面探测、野外地震地质调查等相关工作,评价4条断裂的活动历史、活动性质、规模、展布及最新活动时代。研究结果表明,幕府山-焦山断裂最新活动时代为中更新世晚期,丹徒-建山断裂最新活动时代为早更新世至中更新世早期,五峰山-西来桥断裂最新活动时代为中更新世,汝山-上会断裂为前第四纪断裂。幕府山-焦山断裂和五峰山-西来桥断裂是目标区具备一定发震能力、地震危险性较大的2条断裂;判断五峰山-西来桥断裂潜在地震的最大震级为6.0级,幕府山-焦山断裂潜在地震的最大震级为6.5级。Abstract: Zhenjiang city is an area where destructive earthquakes occur more frequently in Jiangsu. Identification of fault activity and seismic risk assessment to the main seismogenic faults of Zhenjiang city is of great significance to earthquake disaster prevention and reduction in Zhenjiang. The previous studies suggested that there exist four buried faults which have greater threat to Zhenjiang: Mufushan-Jiaoshan fault; Wufengshan-Xilaiqiao fault; Dantu-Jianshan fault; Rushan-Shanghui fault. Based on shallow artificial seismic exploration, drilling joint profiling work and field geological survey in the project of “Urban active fault exploration and seismic risk assessment in Zhenjiang city”, the active history and characteristics, the distribution and the latest activities era of these faults are analysed. The results show that the latest activities era of Mufushan-Jiaoshan fault is Medium term of Middle Pleistocene, the latest activities era of Wufengshan-Xilaiqiao fault is Middle Pleistocene, the latest activities era of Dantu-Jianshan fault is early-middle Pleistocene, and which of Rushan-Shanghui fault is pre-Quaternary. The Mufushan-Jiaoshan and Wufengshan-Xilaiqiao fault are two faults which have certain degree of seismogenic capability and greater seismic risk in the target zone. We judged the maximum potential earthquake of Mufushan-Jiaoshan fault is M6.5, and the maximum potential earthquake of Wufengshan-Xilaiqiao fault is M6.0.
-
引言
如何确定工程场地地震动放大效应,并在工程结构抗震设防中考虑这一影响,一直是工程地震学的一个重要课题,各国的抗震规范都不同程度地考虑了场地土特征对地震动参数的影响 (胡聿贤等,1980;李小军等,2001a;薄景山等,2003;吕悦军等,2008)。
场地条件对地震动参数的影响,一般通过对实际强震资料的统计分析或土层地震反应数值的分析进行研究。李小军等 (2001b)对188个工程场地钻孔土层剖面进行场地地震反应分析,研究了不同条件对地震动峰值加速度影响的规律,给出了各类场地条件下地震动参数变化的经验关系;窦立军等 (2001)按照覆盖层厚度和卓越周期将79个典型工程场地分为3类,以不同峰值的输入地震动进行场地反应计算,给出不同地震动强度的场地影响系数;耿淑伟 (2005)结合相关理论分析,通过对强震记录资料的统计分析,给出适于抗震设计的场地影响系数的建议值。
对于工程抗震设防应用的场地条件对地震动的影响,国家标准《中国地震动参数区划图》(GB 18306-2001) 给出了不同场地条件与平均场地条件下地震动参数的转换关系。各行业抗震设防规范,如《建筑抗震设计规范》(GB 50011-2010) 等,均采用不同场地地震动参数的转换系数进行换算。美国建筑防震安全委员会 (Building Seismic Safety Council,2015) 推荐的抗震设计规范也利用场地系数,以考虑场地条件对地震动峰值加速度和特征周期的影响。吕悦军等 (2011)研究了北京地区的Ⅱ类场地地震动放大效应,并给出了该类场地地震动参数调整系数Ka的建议值,提出场地地震动参数调整系数具有区域性的观点。但在此次研究中,没有考虑不同震级和震中距条件下基岩地震动加速度反应谱的区别,结果只分析了场地条件对地表峰值加速度的影响,而没有考虑各频段加速度反应谱值的变化。然而不同震级和震中距条件下的基岩地震动表现的强度和频谱特性均不一致,因此,研究特定场地条件对不同地震动输入环境下的地震动加速度反应谱的影响,是很有必要的。
本文以北京地区典型中硬场地 (亦称一般场地、Ⅱ类场地) 为研究对象,开展中硬场地对地震动加速度反应谱影响的研究。首先,收集整理了北京地区中硬场地工程地质资料,挑选有代表性的场地作为典型场地,最终选取了110个钻孔并构建土层模型;然后,选择适合本地区的基岩加速度反应谱衰减关系,计算不同震级和震中距条件下 (下文称地震动输入环境) 的基岩地震动加速度反应谱,并合成基岩输入地震动加速度时程,进行土层地震反应计算,对比典型周期所对应的地表加速度反应谱值与基岩加速度反应谱值,得出各典型周期加速度反应谱值的放大倍数,通过比较不同地震动输入环境下各典型周期加速度反应谱值放大倍数之间的差别,分析中硬场地不同地震动输入环境对不同周期加速度反应谱值的影响。
1. 计算不同震级、震中距条件下的基岩输入地震动
为了计算不同震级、震中距条件下的基岩输入地震动,本文按照确定设定地震的基本思路 (陈厚群等,2005),假设不同震级、震中距的地震作为输入源。表 1给出了本文选取的6类不同的设定地震。
表 1 选取的6类设定地震Table 1. Parameters of six different scenario earthquakes震级M 5.5 6.5 7.5 8.0 8.0 8.0 震中距R/km 5 20 70 20 70 200 设定地震 本底地震/弥散地震 近场中强震 中远场中强震 近场大震 中远场大震 远场大震 基岩地震动加速度反应谱衰减关系取目前在华北地区广泛使用的山东海阳核电厂厂址地震安全性评价报告 (中国地震局地球物理研究所,1998) 中的地震动衰减关系,即:
$$ \ln Sa = {c_1} + {c_2}M + {c_3}{M^2} - {c_4}\ln (R + {c_5}{{\rm{e}}^{{c_6}M}}) $$ (1) 式中,Sa代表基岩加速度反应谱值,单位为gal (cm/s2);M是震级;R是震中距,单位为km;c1、c2、c3、c4、c5和c6是随周期而变的回归系数。
本底地震的地震动加速度反应谱衰减关系取美国西部地区的衰减关系式 (李小军等,2005):
$$ \lg Sa = {c_1} + {c_2}M + {c_3}{M^2} + {c_4}\lg (R + {c_5}{{\rm{e}}^{{c_6}M}}) $$ (2) 式中各参数含义同 (1) 式。
根据这两个衰减关系模型进行确定性地震危险性分析计算,得到不同震级、震中距条件下的基岩地震动加速度反应谱 (图 1),表 2列出6类设定地震的基岩峰值加速度。
表 2 6类设定地震的基岩峰值加速度Table 2. Peak accelerations on bedrock of 6 scenario earthquakes设定地震 基岩峰值加速度/gal M=5.5,R=5km 192.32 M=6.5,R=20km 212.15 M=7.5,R=70km 119.40 M=8.0,R=20km 498.01 M=8.0,R=70km 177.85 M=8.0,R=200km 42.76 根据上述基岩加速度反应谱,利用人工地震动时程合成方法,合成土层地震反应分析所需的基岩地震动加速度时程。为了考虑人工地震动合成中随机相位的影响,实际计算中,对6类设定地震各合成了3条基岩地震动加速度时程,共18条。
2. 场地土层反应模型及计算
2.1 建立场地土层模型
本文选取北京地区地震安全性评价报告中110个具有代表性的中硬场地钻孔资料作为场地计算剖面 (挑选标准为剪切波速线性变化较好,无软弱夹层与硬夹层),建立土层地震反应分析模型。已有研究证实场地土动力参数具有明显的地域特征 (荣棉水等,2013;2014),故本文模型中的第四系覆盖层厚度、剪切波速、土动力学参数等资料均取自各钻孔资料的报告。场地模型中土层波速范围为107-592m/s,等效剪切波速值为163-362m/s,覆盖层厚9-96.6m。另外,我们还计算了场地的平均剪切波速,并以此近似计算场地的自振周期。自振周期采用蒋溥等 (1993)的估算方法:T=4d/VS,平均剪切波速值为235-388m/s,自振周期为0.13-1.08s。场地模型的力学参数统计见表 3。场地模型的自振周期与覆盖层厚度有较好的线性关系,如图 2。
表 3 场地模型力学参数统计结果Table 3. Statistical results of mechanical parameters of site models场地模型个数 等效剪切波速值VSE 覆盖层厚度d 平均剪切波速值VS 自振周期T 范围/m·s-1 个数 范围/m 个数 范围/ m·s-1 个数 范围/s 个数 110 150 < VSE≤200 2 5 < d≤15 7 0.1 < T≤0.2 5 200 < VSE≤250 27 15 < d≤30 19 200 < VS≤250 5 0.2 < T≤0.4 38 250 < VSE≤300 52 30 < d≤50 65 250 < VS≤300 21 0.4 < T≤0.6 41 300 < VSE≤350 27 50 < d≤70 12 300 < VS≤350 60 0.6 < T≤0.8 18 350 < VSE≤400 2 70 < d≤90 4 350 < VS≤400 24 0.8 < T≤1.1 8 90 < d≤100 3 2.2 土层地震反应分析计算结果
以合成的18条基岩加速度时程作为输入地震动,采用一维等效线性化波动分析方法 (廖振鹏等,1989;胡聿贤,1999) 分别对110个土层模型进行土层地震反应分析计算,得到地表加速度反应谱。然后计算平均地表加速度反应谱,并与基岩谱进行对比,结果见图 3,其中红线代表平均地表加速度反应谱,绿线代表基岩加速度反应谱。
3. 计算结果分析与讨论
3.1 不同地震动输入环境下地震动加速度反应谱的场地效应分析
(1) 输入近场小震 (本底地震/弥散地震) 地震动时,地震动反应基本近似线性,地震动高频成分 (短周期) 丰富而低频成分 (长周期) 较少。T < 1s时,放大作用明显,0.13-0.6s是放大作用的优势周期段,场地对低频段地震动的放大作用不明显。反应谱峰值对应周期范围为0.1-0.2s。
(2) 输入近场中强震地震动时,场地对T < 1.5的地震动有明显放大作用,其优势周期为0.16-0.6s,同样,对长周期地震动的放大作用不明显,数据离散性较小。反应谱峰值对应周期范围为0.15-0.25s。
(3) 输入中远场中强震地震动时,场地对T≤1.5s的地震动普遍有放大作用,0.13-0.6s为放大作用的优势周期段,对T > 1.5s的长周期地震动放大作用不明显。反应谱峰值对应周期范围为0.35-0.45s。
(4) 输入近场大震地震动时,地震动长周期成分丰富而短周期成分较少,由于基底输入幅值较大,土体非线性特征明显,数据离散性很大。场地对T < 0.05s和0.2-2s周期段的地震动有放大作用,但不明显。0.3-1.2s为优势周期段,反应谱峰值对应周期范围为0.3-0.5s。
(5) 输入中远场大震地震动时,场地对T≤1s的地震动普遍有放大作用,地震动的中频成分较丰富,放大效应的优势周期段为0.15-0.6s。反应谱峰值对应周期范围为0.2-0.6s。
(6) 输入远场大震地震动时,因基岩地震动输入幅值较小,场地对于T < 1s的地震动放大作用明显,且放大效果强于其他5类设定地震,其优势周期段为T≤0.6s,即对高频段和中频段地震动的前半段有显著的放大作用。反应谱峰值对应周期范围为0.15-0.7s。
由此可见,不同地震动输入环境下,场地对地震动的放大作用也不同,放大的优势周期段受震级和震中距的综合影响,当输入基底幅值增大时,放大作用的优势周期段倾向于向长周期偏移,随着震级和震中距的增加,反应谱峰值对应的周期范围也逐渐变大且向长周期方向移动。
3.2 典型周期地震动加速度反应谱的放大效应及其统计分析
为了定量分析场地对地震动的放大效应,我们分别对从高频段 (T < 0.2s)、中频段 (0.2s≤T < 1s)、低频段 (1s≤T) 中选取的10个典型周期 (0s、0.04s、0.07s、0.1s、0.2s、0.5s、0.8s、1.0s、2.0s、4s) 的地表加速度反应谱值进行统计分析,将其与对应的基岩加速度反应谱值相比,得出各典型周期加速度反应谱值的放大倍数。由于土层地震反应分析所计算出的各个土层模型的放大倍数具有一定的离散性,为了能得到较为合理的统计结果,我们对数据进行了异常值的剔除,得到的6类设定地震输入环境下的放大倍数统计图 (图 4),其中横坐标为选取的典型周期T按对数等间距分布,纵坐标为放大倍数k。
本文采用各典型周期加速度反应谱值的放大倍数来反映场地对地震动加速度反应谱的影响。根据土层地震反应的计算结果,我们统计了各地震动输入环境下的放大倍数的均值和方差,统计结果列于表 4,图 5给出放大倍数均值的变化趋势。
表 4 放大倍数均值统计结果Table 4. The statistical results of amplification factors周期T/s 0 0.04 0.07 0.1 0.2 0.5 0.8 1 2 4 M=5.5,R=5km 均值 1.212 1.207 1.169 1.203 1.350 1.350 1.192 1.172 1.062 1.037 方差 0.014 0.014 0.022 0.032 0.028 0.028 0.014 0.010 0.001 0.0007 M=6.5,R=20km 均值 1.148 1.142 1.172 1.217 1.341 1.349 1.203 1.146 1.045 1.017 方差 0.008 0.008 0.024 0.033 0.030 0.033 0.021 0.009 0.001 0.0002 M=7.5,R=70km 均值 1.275 1.277 1.267 1.280 1.338 1.333 1.188 1.127 1.034 1.011 方差 0.008 0.009 0.012 0.033 0.026 0.027 0.016 0.008 0.001 0.0001 M=8.0,R=20km 均值 1.204 1.208 1.123 1.081 1.127 1.286 1.292 1.230 1.103 1.024 M=8.0,R=20km 方差 0.017 0.017 0.030 0.069 0.096 0.058 0.052 0.041 0.011 0.0006 M=8.0,R=70km 均值 1.227 1.214 1.231 1.251 1.344 1.292 1.191 1.130 1.039 1.010 方差 0.006 0.007 0.010 0.024 0.028 0.029 0.020 0.013 0.001 0.0001 M=8.0,R=200km 均值 1.604 1.608 1.640 1.631 1.682 1.432 1.137 1.048 0.987 0.979 方差 0.026 0.025 0.026 0.047 0.074 0.096 0.020 0.006 0.0003 0.0001 综合以上计算与统计结果,我们分析得出以下初步结论:
(1) 可以看出,中硬场地对高频段和中频段的地震动有普遍的放大效应 (放大倍数均值多在1.15以上),尤其对0.2-0.8s周期段的地震动放大效应更为明显 (放大倍数均值多在1.2以上),6类设定地震输入环境下中频段的放大倍数大部分都高于其他2个频段。
(2) 中硬场地对低频段地震动的放大效应比中、高频段弱,6类设定地震输入环境下,周期为1s、2s、4s时放大倍数均值基本都在1.15以下,随着周期的增大,放大倍数呈明显减小的趋势。这与场地自振周期多集中于0.2-0.8s之间,低频段地震动与场地之间几乎不发生共振作用有关。
(3) 在远场大震输入环境下,高频段地震动对应的放大倍数均值 (1.6-1.7之间) 远高于其他地震输入环境,这是因为在此地震输入环境下,高频段输入地震动的强度较弱,场地非线性作用不强,致使场地放大作用较为明显。
(4) 近场大震输入环境下,中硬场地对高频段以及中频段的前半段 (T≤0.5s) 地震动的放大效应相对弱于其他地震输入环境,然而中频段的后半段 (0.8s开始) 以及低频段的放大效应高于其他地震输入环境,这是因为在较强地震动的作用下,场地土的非线性特征明显,土的塑性变形和滞回阻尼增大,基岩地震波在场地地震反应中的高频成分被过滤吸收,中低频成分被放大。
总体来看,中硬场地对0.2-0.8s这一周期段对应的地震动放大效应显著 (放大倍数均值基本都大于1.2),尤其是0.2-0.5s周期段对应地震动,其放大倍数均值大多在1.3以上,这与场地模型的自振周期集中在0.2-0.8s区间段是十分相关的,因地震动与场地的共振作用使得这一周期段的放大效应较强。
3.3 场地地震动放大效应与覆盖层厚度的关系
为了进一步分析场地模型对放大效应的影响,我们对每个土层模型的放大倍数与覆盖层厚度的关系做了统计分析,因篇幅有限,本文仅列出近场小震输入环境下的统计结果 (图 6)。
对统计结果进行初步分析可以看出:
(1) 近场小震输入环境下,对于T≤0.1s的4个典型周期,因较偏离场地自振周期区间,放大倍数随覆盖层厚度的增加而减小。从0.2s、0.5s、0.8s和1s周期下放大倍数红线代表放大倍数随覆盖层厚度变化的线性趋势随覆盖层厚度的变化曲线可看出,放大倍数有先增大后减小的趋势,变化界限厚度分别约为20m、50m、70m和70m,这是因为当典型周期与场地自振周期相近时,地震动与场地的共振作用明显,放大倍数则较大,共振作用也让土层地震反应更为复杂,使数据的离散性更强。对于2s和4s这2个长周期,放大倍数的变化不太明显,因相对远离场地自振周期,以及计算程序对长周期地震动不敏感等因素,使得这2个典型周期放大倍数的离散性相对较小。
(2) 近场中强震输入环境下,放大倍数随覆盖层厚度变化的趋势与近场小震输入环境下相似,有一点不同的是,对于1s周期对应的放大倍数,在覆盖层厚度大于70m时,放大倍数不减小,而是趋于平缓增大。其他周期的放大倍数与覆盖层关系的特征在这里不再赘述。
(3) 中远场中强震输入环境下,对于T≤0.1s的4个典型周期,放大倍数随覆盖层的厚度增加而减小的趋势没有上述2类设定地震输入环境下明显;0.2s、0.5s和0.8s的放大倍数变化趋势与上述2类设定地震输入环境下相同;在长周期段,放大倍数随着覆盖层厚度的增加而增长,增长趋势随着周期变大逐渐变缓,这与输入地震动中长周期成分丰富有关。
(4) 近场大震输入环境下,因计算基底输入地震动强度较大,土层反应剧烈,土体非线性特性明显,数据离散性普遍较大,从0s至0.2s,放大倍数随覆盖层厚度的增加大致呈减小的趋势;0.5s、0.8s和1s对应的放大倍数则先增加后减小,变化界限厚度分别约为33m、70m和70m;2s和4s对应的放大倍数则随覆盖层厚度增加而缓慢增大,可见在假定为近场大震的情况下,长周期地震动的放大效果较为明显。
(5) 中远场大震输入环境下,对于T≤0.1s的4个典型周期,放大倍数随覆盖层厚度的增加而减小,但0s和0.04s这2个典型周期的放大倍数的数据离散性较大,且随覆盖层厚度的增加而较为缓慢地减少;0.2s、0.5s和0.8s周期对应的放大倍数有先增大后减小的趋势,变化界限厚度分别约为20m、50m和70m;在长周期段,放大倍数随着覆盖层厚度的增加而增长,增长趋势随着周期的增大逐渐变缓。
(6) 远场大震输入环境下,高频段周期对应的放大倍数随覆盖层厚度的增加而减小,但减小的趋势不如近场小震和近场中震输入环境下明显;对于0.2s、0.5s和0.8s周期的地震动,放大倍数先增加后减小,变化界限厚度分别约为25m、60m和70m;对于低频段周期,地震动的放大倍数随覆盖层厚度的增加而增加,周期为1s时表现较为明显,2s和4s时的增加趋势较为缓慢。
4. 结语
基于上述分析和讨论,可以得到下列结论:
(1) 不同地震动输入环境下的场地地震动放大作用是不同的,放大作用的优势周期段受震级和震中距的综合影响;当计算输入基底幅值增大时,放大作用的优势周期段倾向于向长周期方向偏移,随着震级和震中距的增加,加速度反应谱峰值对应的周期范围也逐渐变大且向长周期方向移动。
(2) 中硬场地对高、中频段地震动放大效应明显,对0.2-0.8s这一周期段对应地震动的放大效应显著 (放大倍数均值基本都大于1.2),尤其是对0.2-0.5s周期段,放大倍数均值大多在1.3以上,中硬场地对低频段地震动的放大效应较弱 (1s、2s、4s的放大倍数均值基本都在1.15以下),在抗震设计中对于自振周期为0.2-0.5s的建筑要更加慎重。
(3) 场地覆盖层厚度变化对不同频段地震动放大作用的影响是不同的,与场地自振周期的相关性很强。当加速度反应谱周期与场地自振周期相差较大时,场地对周期小于自振周期段的地震动的放大作用随覆盖层厚度增加而减弱,对大于自振周期段的地震动的放大作用随覆盖层厚度增加而缓慢增强;当加速度反应谱周期与场地自振周期相差较小时,场地对地震动的放大作用会呈现先增强后减弱的趋势。
本文在中硬场地对不同周期加速度反应谱值的影响方面进行了相关讨论,对于其他类型的场地并未涉及,在典型周期的挑选方面还不够全面,选取的土层剖面数据还不够充足,我们将在今后的工作中就以上几方面继续进行研究和探讨。
-
表 1 目标区主要断裂特征
Table 1. List of parameters of faults in the target zone
编号 断裂名称 产状 长度/km 性质 最新活动时代 相关地震 走向 倾向 倾角/(°) F1 五峰山-西来桥断裂 NNW E 55~60 34 左行张扭 ${\rm{Q}}_{\rm{p}}^2 $ 999年常州5½级地震 F2 丹徒-建山断裂 NW SW 50~55 40 左行张扭 ${\rm{Q}}_{\rm{p}}^{1-2} $ 1913年镇江5.5级地震 F3 幕府山-焦山断裂 近EW S/N 较陡 75 正断 ${\rm{Q}}_{\rm{p}}^2 $晚期 1624年扬州6级地震、1930年镇江5.5级地震、
1872年镇江4¾级地震F4 汝山-上会断裂 NE E 60~80 24 正断/逆断 前第四纪 — -
邓起东, 卢造勋, 杨主恩, 2007. 城市活动断层探测和断层活动性评价问题. 地震地质, 29(2): 189—200Deng Q. D. , Lu Z. X. , Yang Z. E. , 2007. Remarks on urban active faults exploration and associated activity assessment. Seismology and Geology, 29(2): 189—200. (in Chinese) 侯康明, 张振亚, 刘建达等, 2012. 南京市活断层探测与地震危险性评价. 北京: 地震出版社. 江苏省地质矿产局, 1984. 江苏省及上海市区域地质志. 北京: 地质出版社. 江苏省地质矿产局, 1989. 宁镇山脉地质志. 南京: 江苏科学技术出版社. 李徐生, 韩志勇, 鹿化煜等, 2018. 下蜀黄土底界的年代及其对区域气候变干的指示. 中国科学: 地球科学, 48(2): 210—223.Li X. S., Han Z. Y., Lu H. Y., et al., 2018. Onset of Xiashu loess deposition in Southern China by 0.9 Ma and its implications for regional aridification. Science China Earth Sciences, 61(3): 256—269. 苗巧银, 陈火根, 李向前等, 2016. 江苏镇江东部大路—姚桥断裂新构造运动特征. 地质学刊, 40(2): 220—227Miao Q. Y. , Chen H. G. , Li X. Q. , et al. , 2016. Neotectonic movement characteristics of the Dalu-Yaoqiao Fault in eastern Zhenjiang, Jiangsu province. Journal of Geology, 40(2): 220—227. (in Chinese) 苗巧银, 朱志国, 陈火根等, 2017. 镇江地区长江南北两岸第四纪地层结构划分与沉积特征对比. 华东地质, 38(3): 175—183Miao Q. Y. , Zhu Z. G. , Chen H. G. , et al. , 2017. Classification of Quaternary stratigraphic structures and comparison of sedimentary characteristics on both sides of the Yangtze River in the Zhenjiang area. Resources Survey Environment, 38(3): 175—183. (in Chinese) 邵家骥, 1999. 苏南及沿江地区柏山组、下蜀组的时代及成因. 江苏地质, 23(1): 10—16Shao J. J. , 1999. Time and genesis of Baishan and Xiashu formations along the Yangtze River and South Jiangsu. Jiangsu Geology, 23(1): 10—16. (in Chinese) 许汉刚, 范小平, 冉勇康等, 2016. 郯庐断裂带宿迁段F5断裂浅层地震勘探新证据. 地震地质, 38(1): 31—43Xu H. G. , Fan X. P. , Ran Y. K. , et al. , 2016. New evidences of the Holocene fault in Suqian segment of the Tanlu fault zone discovered by shallow seismic exploration method. Seismology and Geology, 38(1): 31—43. (in Chinese) 张鹏, 李丽梅, 许汉刚等, 2014. 徐州活动断层探测项目主要断层活动性研究. 震灾防御技术, 9(4): 801—812Zhang P. , Li L. M. , Xu H. G. , et al. , 2014. Identification of fault activity in “Urban Active fault exploration and seismic risk assessment in Xuzhou city”. Technology for Earthquake Disaster Prevention, 9(4): 801—812. (in Chinese) 张鹏, 张媛媛, 许汉刚等, 2019. 苏锡常断裂的第四纪活动性. 地震地质, 41(5): 1172—1184Zhang P. , Zhang Y. Y. , Xu H. G. , et al. , 2019. Research on the characteristic of quaternary activities of Su-Xi-Chang fault. Seismology and Geology, 41(5): 1172—1184. (in Chinese) 宗开红, 宗雯, 康从轩等, 2016. 江苏镇江地区主要活动断裂及其晚第四纪活动性研究. 地质力学学报, 22(3): 439—453Zong K. H. , Zong W. , Kang C. X. , et al. , 2016. Research on the characteristics of quaternary activities of Su-Xi-Chang fault. Journal of Geomechanics, 22(3): 439—453. (in Chinese) -