Analysis of Structural Style of Yanjin Fengqiu Section of Xinxiang Shangqiu Fault in Neotectonic Period
-
摘要: 新乡-商丘断裂是南华北盆地与渤海湾盆地、鲁西隆起的分界断裂,是一条长期活动的区域性深大断裂。浅层地震勘探与跨断层钻孔联合剖面探测工作已证实该断裂延津至封丘段的最新活动时代为晚更新世,本文通过对典型地震勘探剖面和钻探资料综合分析认为:新乡-商丘断裂延津至封丘段断裂结构复杂,新近纪以来构造样式在不同构造部位有显著差异。延津段在地震勘探剖面上浅部显示为负花状构造;封丘段变化较大,东部在剖面上为多组阶状正断层组合,形成了由断层带控制的构造背斜,平面上表现为NW向排列、走向NNE雁列状断层组合;西部结构单一,自东向西断裂控制的浅部地层变形范围增大,断层带逐步变宽。根据新乡-商丘断裂延津至封丘段剖面上的构造样式与平面上的组合模式,表明该段断裂最新活动具有走滑运动性质。Abstract: The Xinxiang-Shangqiu fault is a regional deep fracture of long-term activity marking the boundary between the South North China Basin and Bohai Bay Basin. The latest active age of the Yanjin to Fenqiu fragment of this fault is of late Pleistocene, which is testified by the joint work of shallow seismic exploration and cross-fault borehole detection. On the basis of the typical seismic exploration profiles and drilling data, we concluded that the Yanjin-Fengqiu fragment of the Xinxiang-Shangqiu fault has a very complicated structure. Since the Neogene time, the different tectonic positions of the Yanjin-Fengqiu fragment have quite different tectonic styles. The Yanjin part has a negative flower structure in the seismic exploration profile. However, the Fengqiu part is structurally quite variable. The eastern segment of the Fengqiu part is composed of multiple scalariform normal faults, which form an anticline controlled by a fault belt. These faults are arranged in the NW direction and in echelon with a strike of NNE direction. The western segment is relatively in simple structure, with increasing deformation scale of the shallow strata controlled by the E-W trending fault. The fault zone of the western segment becomes wider from east to the west. According to the tectonic styles and combination patterns of the profile, we concluded that the Yanjin-Fengqiu fragment of the Xinxiang-Shangqiu fault is of strike-slip characteristics.
-
Key words:
- Xinxiang Shangqiu fault /
- Seismic exploration /
- Structural style
1) 2 中国地震局地球物理勘探中心,2016. 新乡市活动断层探测与地震危险性评价工程报告.2) 3 中国地震局地球物理勘探中心,2014. 长垣县主城区地震小区划报告.3) 4 河南省航空物探遥感中心,2013. 河南省内黄隆起区地震概查及资源远景评价报告.4) 5 河南省航空物探遥感中心,2019. 河南省地震构造探查工程(1)浅层人工地震探测专题报告.5) 6 中国石油化工股份有限公司中原油田分公司物探研究院,2017. 新乡-商丘断裂石油地震资料分析与解释成果报告.6) 7 中国地震局地球物理勘探中心,2019. 新乡-商丘断层活动构造探测目标断层人工地震浅层反射剖面探测专题报告.7) 8 河南省地质矿产勘查开发局第五地质勘查院,2020. 河南省地震构造探查工程(2)-主要断裂活动性鉴定(部分)成果报告.8) 9 河南省地质调查院,2021. 河南省地震构造探查工程(4)-主要断裂活动性鉴定(部分)成果报告.9) 10 北京防灾科技有限公司,2021. 新乡-商丘断裂(封丘段)活动性鉴定与定位验收报告. -
引言
众所周知,我国是地震频发的国家,地震可能造成高速铁路列车脱轨事故,为此,我国建立了相应的地震报警及紧急处置系统,力争在潜在破坏性地震到达之前发出报警,降低列车发生脱轨或倾覆的概率。目前广泛采用的地震报警参数主要为峰值加速度PGA,然而PGA主要体现地震波幅值,未考虑频谱及持时影响,对于无破坏性的高频小震可能产生误报现象。为避免这种干扰,郭恩栋等(2015)引入谱强度SI作为高速铁路地震报警参数,并给出相应的计算方法,SI从能量方面出发,考虑了幅值及频谱影响,忽略持时的影响,但相比PGA更合理。美国电力研究所EPRI(O'Hara 等,1991)在核电厂地震监测报警中引入累积绝对速度CAV作为报警参数,研究发现低频震动对CAV的贡献明显高于高频震动。我国三代核电站AP1000参考国外标准也引入了CAV作为地震报警参数(胡进军等,2013)。林宜娴(2002)提出CAV、PGA、PGV等参数可作为地震报警参数,其稳定性好,且确定性较高。Kostov(2005)根据欧洲几次强地震动数据,分析得出标准化累积绝对速度CAVSTD在预测地震动破坏能力方面明显优于PGA。伊斯坦布尔在地震预报警系统中利用累积绝对速度快速判别强地震动(Erdik等,2003)。Fahjan等(2011)在地震报警的触警算法中引入CAV,并提出新的用于报警系统的参数。黄俊等(2014)引入CAV作为报警参数,分析其在高速铁路报警中的适用性,研究发现CAV可有效排除破坏性小的高频小震对高速铁路列车地震报警的干扰。目前虽对CAV在地震报警中的应用有了深入研究,但在高速铁路列车地震报警中的研究较少,仅通过相关的地震波数据处理得出,未考虑车轨动力响应关系,且给出的计算参数为CAVSTD,为此,需探讨考虑车轨关系的基于相关起算阈值的CAV计算方法。
1. CAV相关参数
目前,地震工程界主要研究CAV、CAVSTD、CAV5。其中CAV于1988年由美国电力研究所(O'Hara等,1991)在核电厂报警中第2级判别标准中提出的,其概念为记录地震加速度数据的绝对值对记录时间进行积分,记为累积绝对速度CAV:
$${\rm{ CAV}} = \mathop \int \nolimits_0^{{{{t}}_{{\rm{max}}}}} \left| {{{a}}\left( {{t}} \right)} \right|{\rm{d}}t$$ (1) 为提高CAV在核电厂报警中的准确度,EPRI在报告中引入标准化累积绝对速度CAVSTD,其可有效避免长持时记录的地震波尾部对应的幅值小且破坏性小的地震动对其数值的影响,其定义如下:将记录的加速度时程分成以1 s为时间间隔的N个小段,当1 s时间段内存在加速度绝对值≥0.025 g时,对该秒内加速度数值的绝对值进行积分,否则该秒内的积分值记为零,然后将整个记录时程进行相加求和,表达式为:
$$ {\rm{CA{V_{STD}}}} = \mathop \sum \limits_{{i=1}}^n \mathop \int \nolimits_{{i}}^{{{i}} + 1} {W_i}\left| {{{a}}\left( {{t}} \right)} \right|{\rm{d}}t,\;\;\;\;\;\;\;\;{W_i} = \left\{ {\begin{array}{*{20}{c}} {0,\left| {{{a}}\left( {{t}} \right)} \right| < 0.025\;g}\\ {1,\left| {{{a}}\left( {{t}} \right)} \right| \geqslant 0.025\;{{g}}} \end{array}} \right.$$ (2) Kramer等(2002)于2002年提出CAV5,用于场地液化判别,其概念为:将记录的加速度时程中对加速度数值的绝对值≥0.005 g进行积分求和,表达式为:
$${{\rm{ CAV}}_5} = \displaystyle\int \nolimits_0^\infty \left\langle {\rm{\chi }} \right\rangle \left| {{{a}}\left( {{t}} \right)} \right|{\rm{d}}t,\;\;\;\left\langle {\rm{\chi }} \right\rangle = \left\{ {\begin{array}{*{20}{c}} {0\;\left| {{{a}}\left( {{t}} \right)} \right| < 0.005\;g}\\ {1\;\left| {{{a}}\left( {{t}} \right)} \right| \geqslant 0.005\;{{g}}} \end{array}} \right. $$ (3) 式中,a(t)为加速度时程,Wi为权重。
2. 改进的标准化累积绝对速度应用于高速铁路地震报警阈值中的计算方法
2.1 改进的标准化累积绝对速度起算阈值确定
本文引入31自由度的SY97477车辆模型(崔恩文,2014),结合车轨典型非线性接触关系,依据列车安全行驶的3项评判指标(脱轨系数、轮重减载率及横向力),分析车体在轨道不平顺(我国高速铁路无砟轨道标准谱)与简谐波(仅考虑横向输入,与列车行驶方向垂直)共同激励下的动态响应,考虑95%的安全保证,引入2.55倍的动力放大系数(孙汉武等,2007),计算安全极限状态(超出3项指标任意1个限值)下不同车速对应的简谐波幅值随周期的变化关系,如图1所示。具体模型参数、车轨接触关系及评判标准可参考崔恩文(2014)的研究。
由图1可知,车速一定时,即使幅值很小(<0.025 g)的低频简谐波也可能使列车脱轨,因此不能简单地套用标准化累积绝对速度CAVSTD计算公式,有必要研究改进的标准化累积绝对速度应用在高速铁路地震报警阈值中的计算方法。车速一定时,简谐波幅值随周期增加逐渐递减,最终趋于常数,且周期越大,不同车速对应的简谐波幅值之间的差值越小。车速为200 km/h时,当简谐波周期达2.2 s时,其幅值开始趋于常数,该常数约为0.008 g,本文规定此时的周期为简谐波幅值趋于常数对应的起始周期Tst,该常数记为Amin。不同车速下简谐波幅值趋于常数对应的起始周期及该常数的数值如表1所示。
表 1 不同车速下简谐波幅值趋于常数对应的起始周期及该常数的数值Table 1. The initial period of the simple harmonic wave amplitude tending to be constant at different speeds and the value of the constant项目 车速/(km·h−1) 200 250 300 350 400 起始周期Tst/s 2.2 2.3 2.1 2.0 1.7 Amin/g 0.008 0.004 0.004 0.004 0.004 由表1可知,车速为200 km/h与车速≥250 km/h时对应的Amin不同,车速≥250 km/h时Amin相同,但起始周期Tst随车速的增加呈递减趋势,说明车速越大,简谐波趋于常数所需的周期越小。
车速为200 km/h时,幅值低于0.008 g的简谐波对列车运行安全的影响可忽略不计。车速≥250 km/h时,幅值低于0.004 g的简谐波对列车运行安全的影响可忽略不计。由此可推算出改进的标准化累积绝对速度起算阈值,如表2所示。
表 2 不同车速下改进的标准化累积绝对速度起算阈值Table 2. Starting threshold of improved calculation of standardized cumulative absolute velocity at different speeds项目 车速/(km·h−1) 200 ≥250 改进的标准化累积绝对速度起算阈值/g 0.008 0.004 2.2 改进的标准化累积绝对速度应用在地震报警中的可行性分析
考虑改进的标准化累积绝对速度为累积过程的算法,其值对于加速度超过0.004 g或0.008 g的有限个数的高频小震会很小,因此,可较好地排除高速铁路沿线的高频小震,是较好的地震监测报警参数。
魏智祥(2015)分析了京武高速铁路线布置的强震仪监测到的列车运行时震动记录数据,发现列车运行产生的震动记录明显高于背景噪声,且发现东西、南北、垂直方向列车振动噪声幅值均≤0.002 g,由此可知本文设定的不同车速下改进的标准化累积绝对速度可有效排除列车振动的影响。
综上所述,本文给出改进的标准化累积绝对速度用于高速铁路地震报警是可行的。
2.3 改进的标准化累积绝对速度计算公式的确定
参照CAVSTD计算方法,本文给出车速为200 km/h时改进的标准化累积绝对速度CAV8计算公式如下:
$${{\rm{CAV}}_8} = \mathop \sum \limits_{i=1}^n \mathop \int \nolimits_{{i}}^{{{i}} + 1} {W_i}\left| {{{a}}\left( {{t}} \right)} \right|{\rm{d}}t,\;\;\;\;{W_i} = \left\{ {\begin{array}{*{20}{c}} {0,\;\left| {{{a}}\left( {{t}} \right)} \right| < 0.008\;g}\\ {1,\;\left| {{{a}}\left( {{t}} \right)} \right| \geqslant 0.008\;{{g}}} \end{array}} \right.$$ (4) 当
$ a\left(t\right) $ $< 0.008\;g$ 时,${{W}}_{{i}}$ =0;当$\left|{a}\left({t}\right)\right|\geqslant 0.008\;g$ 时,${{W}}_{{i}}$ =1。本文给出的CAV8定义为:将加速度时程按1 s为时间间隔进行划分,如果在这1 s内加速度绝对值≥0.008 g,则对该秒内的加速度绝对值进行积分,否则忽略该秒内的积分值,然后将整个时程中每秒积分得到的数值进行累积求和。
车速为250 km/h时改进的标准化累积绝对速度CAV4计算公式如下:
$${ {\rm{CAV}}_4} = \mathop \sum \limits_{i=1}^{{n}} \mathop \int \nolimits_{{i}}^{{{i}} + 1} {W_i}\left| {{{a}}\left( {{t}} \right)} \right|{\rm{d}}t,\;\;\;\;\;{W_i} = \left\{ {\begin{array}{*{20}{c}} {0,\;\left| {{{a}}\left( {{t}} \right)} \right| < 0.004\;g}\\ {1,\;\left| {{{a}}\left( {{t}} \right)} \right| \geqslant 0.004\;{{g}}} \end{array}} \right.$$ (5) 当
${a}\left({t}\right)$ $< 0.004\;g$ 时,${{W}}_{{i}}=0$ ;当$\left|{a}\left({t}\right)\right|\geqslant 0.004\;{g}$ 时,${{W}}_{{i}}$ =1。本文给出的CAV4定义为:将加速度时程按1 s为时间间隔进行划分,如果在这1 s内加速度绝对值≥0.004 g,则对该秒内的加速度绝对值进行积分,否则忽略该秒内的积分值,然后将整个时程中每秒积分得到的数值进行累积求和。
3. 高速铁路地震报警阈值确定
为确定基于改进的标准化累积绝对速度高速铁路地震报警阈值,本文从Peer数据库选取6条不同场地地震波(表3),对应的加速度反应谱如图2所示,对其进行归一化处理,分析加速度幅值与评判列车脱轨3项指标之间的关系,确定不同车速下影响行车安全的6条地震波加速度上限值,按照本文给出不同车速下改进的标准化累积绝对速度计算公式,给出6条地震波对应的改进的标准化累积绝对速度上限值,如表4所示。
表 3 地震波相关信息Table 3. Seismic wave information项目 地震波台站名称 Cape Mendocino TCU045 Capitola Hector EL Centro Imp. Co. Cent Poe Road 震级 7.01 7.62 6.93 7.13 6.54 6.54 峰值加速度/g 0.376 5 0.008 1 0.528 5 0.265 6 0.357 9 0.446 3 v30/(m·s−1) 567.78 704.64 288.62 726.00 192.05 316.64 加速度记录持时/s 36.0 45.0 39.9 45.3 40.0 22.3 表 4 不同车速下6条地震波加速度上限值与改进的标准化累积绝对速度上限值Table 4. The upper limit of acceleration amplitudes and upper limit improved calculation of standardized cumulative absolute velocity CAV of six seismic waves at different speeds地震波台站名称 车速200 km/h 车速250 km/h 车速300 km/h 车速350 km/h 车速400 km/h 加速度
上限值/
g改进的标准化
累积绝对速度
上限值/g-sec加速度
上限值/
g改进的标准化
累积绝对速度
上限值/g-sec加速度
上限值/
g改进的标准化
累积绝对速度
上限值/g-sec加速度
上限值/
g改进的标准化
累积绝对速度
上限值/g-sec加速度
上限值/
g改进的标准化
累积绝对速度
上限值/g-secCape Mendocino 0.074 0.17 0.070 0.17 0.055 0.13 0.040 0.09 0.039 0.09 TCU045 0.101 0.55 0.090 0.49 0.078 0.43 0.062 0.34 0.039 0.21 Capitola 0.109 0.34 0.078 0.25 0.062 0.20 0.054 0.17 0.039 0.11 Hector 0.101 0.27 0.070 0.19 0.062 0.17 0.054 0.14 0.047 0.12 EL Centro Imp. Co. Cent 0.078 0.16 0.070 0.18 0.054 0.13 0.039 0.08 0.027 0.05 Poe Road 0.074 0.19 0.054 0.14 0.043 0.11 0.039 0.10 0.035 0.08 由表4可知,车速相同时,即使地震波加速度相同,其对应的改进的标准化累积绝对速度不一定相同,如车速为200 km/h时,台站名称为Cape Mendocino与Poe Road的地震波加速度相同,但对应的改进的标准化累积绝对速度不同,且Cape Mendocino台站加速度记录持时较大,说明改进的标准化累积绝对速度与频谱紧密相关。车速≥250 km/h时,对于同一台站的地震波,峰值加速度越大其对应的改进的标准化累积绝对速度越大,如台站名称为Cape Mendocino的地震波,说明对于同一地震波而言,加速度影响累积绝对速度。根据改进的标准化累积绝对速度定义,加速度记录持时影响其数值。综上所述,改进的标准化累积绝对速度数值不仅与加速度幅值有关,还与频谱及加速度记录持时密切相关,可较好地反映地震动三要素。
选取6条地震波改进的标准化累积绝对速度最小值作为报警阈值,给出不同车速下基于改进的标准化累积绝对速度高速铁路地震报警阈值,如表5所示。由表5可知,改进的标准化累积绝对速度报警阈值与车速整体呈线性递减关系,中间车速对应的改进的标准化累积绝对速度可按线性插入计算。
表 5 不同车速下改进的标准化累积绝对速度高速铁路地震报警阈值Table 5. High-speed earthquake alarm threshold based on improved calculation of standardized cumulative absolute velocity thresholds at different vehicle speeds项目 车速/(km·h−1) 200 250 300 350 400 报警阈值/g-sec 0.16 0.14 0.11 0.08 0.05 4. 结论
(1)通过建立车轨模型,分析列车安全运行极限状态下简谐波幅值与周期的关系,得出不同车速下改进的标准化累积绝对速度起算阈值,并论证改进的标准化累积绝对速度在高速铁路地震报警中的可行性。
(2)参照标准化累积绝对速度CAVSTD计算公式,本文给出不同车速下改进的标准化累积绝对速度计算公式。
(3)根据本文给出的改进的标准化累积绝对速度计算公式,选取6条不同场地的地震波,分析列车在轨道不平顺及地震波共同激励下的动力响应,计算不同车速下改进的标准化累积绝对速度报警阈值,以供参考。
-
图 3 SP1线偏移时间剖面③
Figure 3. Time seismic exploration profile of SP1 line
图 4 XSB1线偏移时间剖面图④
Figure 4. Time profile of XSB1 line offset
图 5 SP2线与SP3线偏移时间剖面⑤
Figure 5. Seismic time profile of line SP2 and line SP3
图 7 SP6线偏移时间剖面⑤
Figure 7. Seismic time profile of line SP6
图 8 新乡-商丘断裂与黄河断裂交汇区构造图6(位置见图2,构造层为新近系底板,有修改)
Figure 8. Structural map of the intersection area of Xinxiang Shangqiu fault and Yellow River fault (See Fig. 2 for location, the structural layer is the Neogene floor which is modified according)
图 9 XS2线偏移时间剖面④
Figure 9. Shallow seismic exploration section of XS2 line
图 10 FQ1线偏移时间剖面⑥
Figure 10. Shallow seismic exploration section of FQ1 line
图 11 新乡-商丘断裂延津段朱庄钻孔联合地质剖面⑧
Figure 11. Combined geological profile of Zhuzhuang borehole in Yanjin section of Xinxiang Shangqiu fault
表 1 浅层地震勘探基本参数
Table 1. Basic parameters of shallow seismic exploration
序号 线号 长度
/m道距/m 炮距/m 偏移距/m 接收道数 覆盖次数 采样间隔/ms 记录时间
/s震源吨位/t 1 XS1 8875 10 20 0 240 60 0.5 3 26 2 XSB1 5795 5 20 0 260 32 0.5 3 26 3 XS2 8100 10 20 0 240 60 0.5 3 26 4 FQ1 7216 4 16 0 220 22 0.5 2 18 注:XSB1线为组合线,由XSB11线与XSB12线组合而成。 表 2 朱庄钻孔联合剖面测年数据表⑧
Table 2. Data sheet of combined profile dating of Zhuzhuang borehole
样品编号 样品埋深/m 地层/标志层 所处构造位置 测试方法 距今年龄/ka ZKZ0714C01 16.55 层⑮ 上覆地层 14C 8.76±0.03 ZKZ0314C02 28.26 层⑬/B7 上覆地层 14C >43.50 ZKZ02OSL03 42.38 层⑬ 上覆地层 OSL 130.55±7.59 ZKZ07OSL01 43.98 层⑫/B6 断层下盘 OSL 134.95±8.85 ZKZ04OSL02 54.18 层⑨ 断层下盘 OSL 278.78±23.77 ZKZ09OSL01 72.45 层⑦/B4 断层下盘 OSL 162.41±9.21 ZKZ02ESR01 77.40 层⑦/B4 断层上盘 ESR 270±36 ZKZ02ESR02 88.70 层⑤/B3 断层上盘 ESR 346±34 ZKZ02ESR03 96.59 层③/B2 断层上盘 ESR 556±72 ZKZ02ESR05 117.43 层① 断层上盘 ESR 853±171 注:14C样品年龄测试由美国Beta实验室完成;OSL、ESR样品年龄测试由北京光释光实验室有限公司完成。 表 3 陈王庄钻孔联合剖面测年数据表⑧
Table 3. Data sheet of combined profile dating of Chenwangzhuang borehole
样品编号 样品埋深/m 地层/标志层 所处构造位置 测试方法 距今年龄/ka ZKW0814C01 18.84 层⑬ 上覆地层 14C 17.15±0.06 ZKW0414C01 25.05 层⑫ 上覆地层 14C >43.5 ZKW05OSL02 32.34 层⑪ 断层上盘 OSL 5.38±0.34 ZKW03OSL02 33.75 层⑪ 断层下盘 OSL 107.39±7.13 ZKW04OSL02 45.96 层⑩ 断层上盘 OSL 138.06±6.79 ZKW09OSL02 73.58 层⑨ 断层上盘 OSL 136.19±7.55 ZKW03ESR01 77.14 层⑧/B3 断层下盘 ESR 286±57 ZKW03ESR02 82.21 层⑧/B3 断层下盘 ESR 482±48 ZKW03ESR04 93.10 层⑥/B2 断层下盘 ESR 546±54 ZKW03ESR05 101.97 层⑤ 断层下盘 ESR 610±117 ZKW03ESR06 120.49 层② 断层下盘 ESR 742±148 注:14C样品年龄测试由美国Beta实验室完成;OSL、ESR样品年龄测试由北京光释光实验室有限公司完成。 -
[1] 侯江飞, 邢磊, 张扬等, 2021. 新乡—商丘断裂延津段浅部地层结构特征研究. 工程地球物理学报, 18(4): 486—494 doi: 10.3969/j.issn.1672-7940.2021.04.0011Hou J. F. , Xing L. , Zhang Y. , et al. , 2021. The shallow structural characteristics of the Yanjin section of Xinxiang—Shangqiu fault. Chinese Journal of Engineering Geophysics, 18(4): 486—494. (in Chinese) doi: 10.3969/j.issn.1672-7940.2021.04.0011 [2] 李艳友, 漆家福, 周赏, 2017. 走滑构造差异变形特征及其主控因素分析——基于砂箱模拟实验. 石油实验地质, 39(5): 711—715 doi: 10.11781/sysydz201705711Li Y. Y. , Qi J. F. , Zhou S. , 2017. Differential deformation and its main controls on strike-slip structures: evidence from sandbox experiments. Petroleum Geology & Experiment, 39(5): 711—715. (in Chinese) doi: 10.11781/sysydz201705711 [3] 刘尧兴, 周庆, 荆制国, 2001. 豫北地区新构造活动特征及中长期地震预测研究. 西安: 西安地质出版社. [4] 马宝军, 漆家福, 于福生, 2017. 东濮凹陷构造变形的物理模拟研究. 地球学报, 38(3): 430—438 doi: 10.3975/cagsb.2017.03.13Ma B. J. , Qi J. F. , Yu F. S. , 2017. An analysis of physical modeling of tectonic deformation in Dongpu sag. Acta Geoscientica Sinica, 38(3): 430—438. (in Chinese) doi: 10.3975/cagsb.2017.03.13 [5] 王定一, 汤锡元, 陈乃明, 1994. 开封坳陷构造特征、形成演化与油气远景. 石油学报, 15(2): 39—47 doi: 10.3321/j.issn:0253-2697.1994.02.005Wang D. Y. , Tang X. Y. , Chen N. M. , 1994. Structural features, evolution and oil and gas prospects of Kai Feng depression. Acta Petrolei Sinica, 15(2): 39—47. (in Chinese) doi: 10.3321/j.issn:0253-2697.1994.02.005 [6] 王志铄, 2017. 河南省地震构造特征. 北京: 地质出版社. [7] 谢智, 谢健健, 王磊等, 2010. 豫北地区地震活动性研究. 华南地震, 30(3): 42—52 doi: 10.3969/j.issn.1001-8662.2010.03.006Xie Z. , Xie J. J. , Wang L. , et al. , 2010. Study on earthquake activity in Yubei region. South China Journal of Seismology, 30(3): 42—52. (in Chinese) doi: 10.3969/j.issn.1001-8662.2010.03.006 [8] 徐汉林, 赵宗举, 杨以宁等, 2003. 南华北盆地构造格局及构造样式. 地球学报, 24(1): 27—33 doi: 10.3321/j.issn:1006-3021.2003.01.005Xu H. L. , Zhao Z. J. , Yang Y. N. , et al. , 2003. Structural pattern and structural style of the southern North China basin. Acta Geoscientica Sinica, 24(1): 27—33. (in Chinese) doi: 10.3321/j.issn:1006-3021.2003.01.005 [9] 许立青, 李三忠, 索艳慧等, 2013. 华北地块南部断裂体系新构造活动特征. 地学前缘, 20(4): 75—87Xu L. Q. , Li S. Z. , Suo Y. H. , et al. , 2013. Neotectonic activity and its kinematics of fault system in the south of North China block. Earth Science Frontiers, 20(4): 75—87. (in Chinese) [10] 许忠淮, 汪素云, 黄雨蕊等, 1989. 由大量的地震资料推断的我国大陆构造应力场. 地球物理学报, 32(6): 636—647 doi: 10.3321/j.issn:0001-5733.1989.06.004Xu Z. H. , Wang S. Y. , Huang Y. R. , et al. , 1989. The tectonic stress field of Chinese continent deduced from a great number of earthquakes. Chinese Journal of Geophysics, 32(6): 636—647. (in Chinese) doi: 10.3321/j.issn:0001-5733.1989.06.004 [11] 张扬, 王志铄, 周栋梁等, 2018. 南华北盆地沈丘凹陷新构造时期断裂活动特征. 地震地磁观测与研究, 39(5): 30—38 doi: 10.3969/j.issn.1003-3246.2018.05.005Zhang Y. , Wang Z. S. , Zhou D. L. , et al. , 2018. Shenqiu depression faults movement characteristics in southern North China basin during the neotectonic period. Seismological and Geomagnetic Observation and Research, 39(5): 30—38. (in Chinese) doi: 10.3969/j.issn.1003-3246.2018.05.005 [12] 张扬, 贺承广, 鲁人齐等, 2021. 新乡—商丘断裂延津段活动特征与晚第四纪地层沉积关系研究. 地质论评, 67(S1): 11—14 doi: 10.16509/j.georeview.2021.s1.005Zhang Y. , He C. G. , Lu R. Q. , et al. , 2021. Relationship of the activity characteristics in the Yanjin segment of the Xinxiang-Shangqiu fault and the late-quaternary stratigraphic sedimentary. Geological Review, 67(S1): 11—14. (in Chinese) doi: 10.16509/j.georeview.2021.s1.005 -