Quaternary Activity Characteristic of the Kaifeng Fault
-
摘要: 开封断裂是郑州-开封断裂带的重要组成部分,也是黄河中下游平原的一条重要控制性断裂。利用浅层地震勘探发现,开封断裂为总体走向EW,倾向N的正断层;断裂分为东、西2支,断裂活动性具有明显的分段特征,东支活动性较西支强。钻孔联合剖面揭示开封断裂东支断裂上断点埋深27~35 m,根据地层年代学结果,其最新活动时间为晚更新世。研究成果为开封市城乡规划、重大工程选址和地震区划工作提供了重要的基础数据。Abstract: Kaifeng fault is an important segment of Zhengzhou-Kaifeng fault zone,and it is also an important controlling fault in the middle and lower reaches of the Yellow River Plain. It is found that Kaifeng fault is a normal fault with overall strike EW and dip N by using shallow seismic exploration. The fault can be divided into East and West branches, and has obvious differentiated segmented characteristics. The activity of the East branch is stronger than that of the West Branch. The borehole joint geological profile reveals that the buried depth of the breakpoint on the fault is between 27 m and 35 m. Combining with the chronological results, the latest age of activity is the late Pleistocene. The results can provide important basic data for Kaifeng urban and rural planning, major project site selection and seismic zoning.1) 2河南省平原第四纪地质图及地质图说明书(1/50万). 郑州:河南省地质局. 开封市活动断层探测与地震危险性评价. 郑州:河南省地震局.
-
引言
在实际地震中,地震动输入具有一定的不确定性,钢筋混凝土框架结构需要经受双向地震作用的考验。在双向受力下,节点核心区存在较大的水平剪力,产生的破坏往往是剪切破坏。通过以往的震害研究可以看出,节点的破坏是框架结构发生破坏的主要原因之一。因此,如何保证节点具有足够的抗震能力成为亟待解决的问题。
各国学者针对平面边节点展开了一系列研究,唐九如等(1985)、Ghobarah等(2005)通过大量边节点试验研究,发现普通边节点抗剪承载力低于中节点,并提出边节点存在混凝土斜压杆。针对空间节点,陈永春等(1995)、雷远德等(2014)通过研究,认为双向往复荷载作用较单向往复荷载作用更容易造成框架柱柱端或节点区破坏,承载力、延性等抗震性能退化明显,且同样出现了边柱节点在节点区附近发生脆性破坏,边柱节点承载力明显低于中柱节点的情况(框架节点专题研究组,1983;傅剑平,2002)。Hwang等(2000)基于已有研究提出了应用于平面节点的软化拉-压杆模型,这一模型是在压杆-拉杆模型的基础上考虑混凝土受压开裂后的软化特性提出的,用以计算节点区抗剪强度。但目前关于软化拉-压杆模型在空间边节点抗剪承载力计算的研究尚未完善。
因此,本文以水平双向受力相同的空间边节点为研究对象,在软化拉-压杆模型平面边节点抗剪计算模型的基础上,建立空间边节点抗剪计算方法,并通过数值模拟的方法对计算方法进行验证。
1. 空间边节点抗剪承载力计算模型
1.1 空间边节点受力分析
由于空间边节点构造的特殊性,与空间中节点相比,其一端缺少梁约束。节点区内的柱筋及梁筋造成节点上下侧的柱端、3个方向的梁端分别受拉或受压,同时受压混凝土形成相应受压区。在双向受力下,节点需要承受2个方向的力,如图 1所示,其中,Mc、Nc、Vc分别为柱端弯矩、轴力和剪力;Mbx、Nbx、Vbx分别为x向梁端弯矩、轴力和剪力;Mby、Nby、Vby分别为y向梁端弯矩、轴力和剪力。节点需要承受混凝土传递的压力、梁筋传递的力以及柱筋传递的力。节点核心区内,在2个方向分别形成节点水平剪力Vhx和Vhy,通过合成得到水平的斜向合成剪力Vh,2个方向的竖向剪力在节点核心区内合成为竖向合成剪力Vv。空间边节点抗剪机制如图 2所示,其中,hc、bs、hb分别为柱高、节点宽和梁高,${{h}_{c}}^{\prime \prime }$、${{h}_{s}}^{\prime \prime }$、${{h}_{b}}^{\prime \prime }$为相应侧柱筋及梁筋中心之间的距离。
本文主要针对空间节点中水平双向受力相同的情况进行研究,因此水平合剪力Vh与水平方向夹角为45°,竖向合剪力Vv垂直于水平面。节点核心区合剪力作用面内,由混凝土形成压杆、钢筋形成拉杆,组成拉-压杆模型承担节点合剪力。因此,在空间节点的合剪力作用面内建立拉-压杆模型。其中,θ为斜压杆倾角,α为合剪力作用面与y向主轴的夹角。
1.2 建立计算模型
在研究双向受力下节点抗剪计算时,主要针对2个方向受力相同这一情况。斜向混凝土斜压杆与水平合剪力、竖向合剪力均作用在同一平面内,在此作用面内建立节点的抗剪计算模型,该合剪力作用面与主轴夹角为45°。双向受力下的软化拉-压杆模型如图 3所示,合剪力作用面内虚线代表混凝土形成的压杆,实线代表钢筋转化成的拉杆。D为节点区斜压杆承担的斜向压力,Fh为节点内水平箍筋所形成的水平拉杆的拉力,Fv为节点内抗剪柱筋形成的竖向拉杆的拉力。
空间边节点的水平抗剪承载力Vjh表示为(Hwang等,2002):
$$ {{V}_{jh}}=K\xi {{{f}'}_{c}}{{A}_{str}}\text{cos}\theta $$ (1) 其中,K为拉压杆系数;ξ为混凝土软化系数;${{{f}'}_{c}}$为混凝土抗压强度;Astr为斜压杆有效截面面积;θ为斜压杆与水平夹角。
节点核心区的剪力主要由混凝土形成的压杆承担,节点内的水平箍筋及竖向中间柱筋使核心区更多的混凝土形成压杆用,使得节点核心区在主斜压杆方向能够承担更多的压力,从而达到提高节点抗剪承载力的作用。因此,当节点内存在水平箍筋及竖向中间柱筋形成的拉杆时,斜向受压在斜压杆传递基础上增加了额外的传力路径,引发了更多核心区混凝土参加抗剪,从而增大了节点区抗剪强度。在计算中,将拉杆的贡献用系数K表示,计算公式为:
$$ K={{K}_{h}}+{{K}_{v}}-1 $$ (2) 其中,Kh为水平拉杆系数;Kv为竖向拉杆系数,可分别通过式(3)、(4)进行计算:
$$ {{K}_{h}}=1+({{\overline{K}}_{h}}-1)\times \frac{{{F}_{yh}}}{{{\overline{F}}_{h}}}\le {{\overline{K}}_{h}} $$ (3) $$ {{K}_{v}}=1+({{\overline{K}}_{v}}-1)\times \frac{{{F}_{yv}}}{{{\overline{F}}_{v}}}\le {{\overline{K}}_{v}} $$ (4) $$ {{\bar{K}}_{h}}=\frac{(1-{{\gamma }_{h}})+{{\gamma }_{h}}}{(1-{{\gamma }_{h}})+{{\gamma }_{h}}\left( 1-\frac{\text{si}{{\text{n}}^{2}}\theta }{2} \right)}\ge 1 $$ (5) $$ {{\bar{K}}_{v}}=\frac{(1-{{\gamma }_{v}})+{{\gamma }_{v}}}{(1-{{\gamma }_{v}})+{{\gamma }_{v}}\left( 1-\frac{\rm{co}{{\rm{s}}^{\rm{2}}}\theta\rm{ }}{\mathit{2}} \right)}\ge 1 $$ (6) $$ {{\gamma }_{h}}=\frac{2\tan \theta-1}{3}(0\le {{\gamma }_{h}}\le 1) $$ (7) $$ {{\gamma }_{v}}=\frac{2\cot \theta-1}{3}(0\le {{\gamma }_{v}}\le 1) $$ (8) 其中,${{\overline{K}}_{v}}$为水平拉杆平衡系数、${{\overline{K}}_{v}}$为竖向拉杆平衡系数;γh为水平拉杆的拉力与节点水平剪力的比值、γv为竖向拉杆的拉力与节点竖向剪力的比值。
水平及竖向拉杆的拉力可分别用式(9)、(10)进行计算:
$$ {{\overline{F}}_{h}}={{\gamma }_{h}}\times ({{\overline{K}}_{h}}\xi {{{f}'}_{c}}{{A}_{str}})\times \text{cos}\theta $$ (9) $$ {{\overline{F}}_{v}}={{\gamma }_{v}}\times ({{\overline{K}}_{v}}\xi {{{f}'}_{c}}{{A}_{str}})\times \text{sin}\theta $$ (10) 其中,${{\overline{F}}_{h}}$为水平拉杆拉力、${{\overline{F}}_{v}}$为竖向拉杆拉力。
在水平机构计算中,定义节点箍筋所提供的抗剪力作用在核心区45°方向,在计算箍筋面积时,认为空间节点箍筋的有效面积为平面节点箍筋有效面积的$\sqrt{2}$倍。因此,Fyh=$\sqrt{2}$Athfyh,Ath为节点水平箍筋各肢的总截面面积,fyh为节点水平箍筋的屈服强度。竖向机构主要由柱纵筋提供抗剪承载力,在空间模型中2个方向的剪力合力均为竖直方向。因此,Fyv=Atvfyv,Atv为节点区内抗剪柱筋截面面积,fyv为节点中间柱筋的屈服强度。
节点所能抵抗的最大水平剪力主要由主斜压杆端部所能抵抗的最大压力决定,而主斜压杆端部的抗压能力由混凝土抗压强度决定。Hwang等(2002)通过研究发现开裂后的钢筋混凝土抗压强度远小于单轴受压混凝土的抗压强度,同时混凝土的抗压强度会随柱拉应变的增大而减小,发生了软化效应。在节点抗剪过程中,将节点主斜压杆面积范围内混凝土达到极限的抗压强度,确定为节点抗剪失效的标准。因此,在计算节点抗剪承载力时需要考虑混凝土抗压强度软化系数ξ,其计算公式为:
$$ \xi \approx \frac{3.35}{\sqrt{{{{{f}'}}_{c}}}}\le 0.52 $$ (11) 其中,${{{f}'}_{c}}$为混凝土圆柱体抗压强度。
根据欧洲CEB-FIP规范(Euro-International Concrete committee,1993),混凝土圆柱体抗压强度与混凝土立方体抗压强度标准值之间的换算关系为${{{f}'}_{c}}$=0.80fcu, k;根据中国《混凝土结构设计规范》(GB50011—2010)(中华人民共和国建设部,2010),fck=0.88αc1αc2fcu, k, fck为混凝土轴心抗压强度标准值。
1.3 斜压杆的确定
空间边节点在受压区混凝土的作用下,节点核心区内部形成了斜向受压杆(崔建宇等,2014)。梁端及柱端传入节点核心区的压力在柱端上、下部柱端角部形成了受压区,两端受压区所形成的斜三角柱为斜压杆范围。定义下表面(ΔLON)为等效斜压杆横截面。下表面中N、L点的位置可通过2个方向梁受压区宽度确定。根据斜压杆受压特点,假设等效斜压杆横截面与节点斜对角线垂直,L、O、N点相连与水平面形成一定角度的三角形。定义ΔLON的面积为斜压杆有效抗压面积Astr,核心区立方体对角线与地面形成的夹角为斜压杆角度θ。
斜压杆角度θ及斜压杆有效抗压面积Astr可表示为:
$$ \theta =\text{arctan}\frac{{{{{h}''}}_{b}}}{\sqrt{{{{{h}''}}_{c}}^{2}+{{{{b}''}}_{s}}^{2}}} $$ (12) $$ {{A}_{str}}={{A}_{LON}}=\frac{{{A}_{LCN}}}{\text{sin}\theta } $$ (13) 其中,ALON、ALCN分别为ΔLON、ΔLCN的面积。
lc为柱下端受压区ΔLCN直角边长度,lc可按式(14)计算;ΔLON的面积可通过式(15)计算,即可得到空间边节点斜压杆有效抗压面积Astr。
$$ {{l}_{c}}=\left(\text{0}\text{.85}+\text{0}\text{.25}\frac{N}{{{A}_{\text{g}}}{{{{f}'}}_{c}}} \right){{h}_{c}} $$ (14) $$ {{A}_{str}}={{A}_{LON}}=0.707{{l}_{c}}\sqrt{0.5l_{c}^{2}+a_{b}^{2}} $$ (15) 2. 计算模型的验证
利用有限元分析软件ABAQUS进行数值模拟。模型的混凝土单元采用C3D8R实体单元;钢筋单元根据钢筋受力特点,采用桁架单元T3D2。模型边界条件柱顶端为铰接,柱底为完全铰接,梁端为加载自由端。为验证数值模拟的准确性,对王玉雷(2010)研究中的试验试件JD4进行数值分析。
表 1、图 5为模拟结果与试验结果对比及荷载位移曲线模拟值与试验值对比图,通过对比分析可知有限元模拟结果与试验结果吻合良好。
表 1 有限元模拟值与试验值对比Table 1. Comparison of simulated results and test results对比项 屈服荷载/kN 屈服位移/mm 峰值荷载/kN 峰值位移/mm 极限荷载/kN 极限位移/mm 试验值 139.9 13 163.9 45 152.8 60 模拟值 148.4 11.5 161.8 41.8 137.6 54.6 为验证计算方法的准确度,建立了6个空间边节点模型,模型具体信息见图 6。在模型中,对柱施加轴向压力,同时在梁端双向加载,分别施加向上(向下)的力,模拟实际情况中边节点在双向受力下的受力过程。同时,将空间角节点这一特殊边节点情况考虑在内,因此加入3个角节点试验试件,并将计算结果与模拟结果或试验结果进行对比。
将计算值与模拟得到的节点抗剪承载力模拟值进行对比,模型SP1—SP3为崔建宇等(2014)给出的3个空间角节点试验试件,对其进行抗剪承载力计算,并将计算结果与试验值进行对比,具体数值见表 2。通过表 2的数据可以看出,计算值与模拟值或试验值较为接近。
表 2 空间边节点参数及计算值与模拟值对比Table 2. Comparison of spatial joints between calculated results and simulated results模型编号 节点形式 hb/mm hc/mm fck/MPa 轴压比 柱纵筋配筋 箍筋配筋 计算值Vjh, c/kN 模拟值Vjh, a/kN Vjh, c/Vjh, a B-1 边节点 450 400 27.9 0.1 12C18+4C22 A10@50 474.17 550.09 0.862 B-2 边节点 450 400 27.9 0.1 12C18+4C22 A12@50 484.77 573.71 0.845 B-3 边节点 450 400 27.9 0.15 12C18+4C22 A10@50 486.51 552.99 0.879 B-4 边节点 450 400 27.9 0.3 12C18+4C22 A10@50 524.57 589.44 0.890 B-5 边节点 450 400 27.9 0.3 12C18+4C22 A12@50 535.17 602.73 0.888 B-6 边节点 450 400 27.9 0.45 12C18+4C22 A10@50 564.20 616.74 0.915 SP1 角节点 250 220 34.3 0.2 16D13 A6@80 216.18 262* 0.825 SP2 角节点 250 220 29.1 0.6 16D13 A6@80 225.73 247.8* 0.911 SP3 角节点 250 220 29.1 0.6 16D13 A10@80 253.41 275.8* 0.919 注:fck为混凝土轴心抗压强度标准值;Vjh, a值中上标*的为试验值。 3. 结论
(1) 通过对平面软化拉—压杆模型和双向受力下钢筋混凝土边节点受力机理分析,在节点核心区双向剪力合成作用面内,建立了双向受力下钢筋混凝土框架边节点抗剪承载力计算方法。
(2) 基于软化拉-压杆模型,建立了空间边节点抗剪计算方法,确定了空间边节点核心区斜压杆截面面积及倾角的取值方法。对双向受力下软化拉-压杆模型中的拉杆进行了定义。
(3) 利用有限元分析软件ABAQUS建立不同参数下空间边节点模型,得到相应抗剪承载力模拟值。将计算值与模拟值或试验值进行对比,结果表明,双向受力下钢筋混凝土框架边节点抗剪计算方法的准确度良好。
-
表 1 观测系统参数表
Table 1. Parameter table of the observation system
线号 道距/m 炮距/m 偏移距/m 接收道数 覆盖次数 采样间隔/ms 前放增益/db 记录时间/s 震源吨位/t ZK1 5 20 0 240 30 0.5 12 3 26 ZK2 5 20 0 260 32 0.5 12 3 26 ZK3 4 16 0 320 40 0.5 12 2 9 ZK4 5 20 0 260 32 0.5 12 3 26 表 2 勘察断点参数表
Table 2. Parameter table of the upper offset point of fault
资料来源 探测断裂 上断点埋深 断距特征 断裂活动性判断 收集石油测线 开封断裂 浅层数据被切除,
可识别最大埋深300 m新近系底部断距150 m PreQ 古近系底界 2 000 m ZK1 开封断裂西支 95 m 3~5 m $ {\mathrm{Q}}_{\mathrm{P}}^{2} $ ZK4 开封断裂西支 139 m 5~8 m $ {\mathrm{Q}}_{\mathrm{P}}^{2} $ ZK2 开封断裂东支 50 m 2~3 m $ {\mathrm{Q}}_{\mathrm{P}}^{3} $ ZK3 开封断裂东支 45~50 m 1~3 m $ {\mathrm{Q}}_{\mathrm{P}}^{3} $ 表 3 钻孔联合剖面主要地层层断距
Table 3. Respective displacement of the mark layer in composite geological profile
主要标志层 地层岩性 下盘埋深/m 上盘埋深/m 断距/m 地层①泥炭层 灰黑色泥碳层 Z2:15.29~15.95
Z4:15.45~15.74
Z8:15.81~16.21
Z9:15.67~16.42Z1:15.37~16.55
Z7:15.49~15.81
Z5:15.97~16.04
Z10:15.47~15.81— 地层②底界面 上部为棕灰、棕黄色细砂,
下部棕灰、棕黄色黏土-粉土,
界线明显Z1:35.29
Z6:35.53
Z8:35.84
Z9:35.74
Z11:35.23
Z12:35.84Z1:36.36
Z3:36.78
Z5:36.55
Z7:37.42
Z10:36.990.75 地层③泥炭夹层 黑色泥炭层 Z2:40.84~41.47
Z4:40.20~40.77
Z6:39.92~40.37
Z8:40.93~42.35
Z9:40.14~42.51
Z11:40.19~40.34Z5:42.23~42.38
Z7:42.03~42.24
Z10:42.96~43.261 地层③底界面 上部为黄棕色、灰棕色粉砂,下部为深棕色、黄棕色粉质黏土,界线清晰 Z2:50.20
Z3:50.81
Z4:50.26
Z11:50.74Z7:52.24
Z5:51.80
Z10:51.971.7 地层④底界面 上部深棕色粉砂,下部为深棕色黏土分界面,界线清晰 Z2:60.37
Z4:60.51
Z9:62.48Z5:65.58
Z7:65.61
Z1:66.443.1 地层⑤底界面 上部为棕黄色、浅棕色粉砂
下部深棕色黏土层,界线清晰Z2:77.48
Z3:79.20
Z4:78.27Z1:88.64
Z5:86.987.1 -
[1] 柴炽章, 孟广魁, 杜鹏等, 2006. 隐伏活动断层的多层次综合探测——以银川隐伏活动断层为例. 地震地质, 28(4): 536—546 doi: 10.3969/j.issn.0253-4967.2006.04.002Chai C. Z. , Meng G. K. , Du P. , et al. , 2006. Comprehensive multi-level exploration of buried active fault: an example of Yinchuan buried active fault. Seismology and Geology, 28(4): 536—546. (in Chinese) doi: 10.3969/j.issn.0253-4967.2006.04.002 [2] 邓起东, 徐锡伟, 张先康等, 2003. 城市活动断裂探测的方法和技术. 地学前缘, 10(1): 93—104 doi: 10.3321/j.issn:1005-2321.2003.01.012Deng Q. D. , Xu X. W. , Zhang X. K. , et al. , 2003. Methods and techniques for surveying and prospecting active faults in urban areas. Earth Science Frontiers, 10(1): 93—104. (in Chinese) doi: 10.3321/j.issn:1005-2321.2003.01.012 [3] 侯江飞, 邢磊, 张扬等, 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 [4] 雷启云, 柴炽章, 郑文俊等, 2014. 钻探揭示的黄河断裂北段活动性和滑动速率. 地震地质, 36(2): 464—477 doi: 10.3969/j.issn.0253-4967.2014.02.015Lei Q. Y. , Chai C. Z. , Zheng W. J. , et al. , 2014. Activity and slip rate of the northern section of Yellow River fault revealed by drilling. Seismology and Geology, 36(2): 464—477. (in Chinese) doi: 10.3969/j.issn.0253-4967.2014.02.015 [5] 李峰, 环文林, 2021. 粤东滨海断裂带第四纪活动特征研究. 震灾防御技术, 16(1): 19—28 doi: 10.11899/zzfy20210103Li F. , Huan W. L. , 2021. Quaternary activity characteristics of the coastal fault zone in eastern Guangdong. Technology for Earthquake Disaster Prevention, 16(1): 19—28. (in Chinese) doi: 10.11899/zzfy20210103 [6] 刘保金, 柴炽章, 酆少英等, 2008. 第四纪沉积区断层及其上断点探测的地震方法技术——以银川隐伏活动断层为例. 地球物理学报, 51(5): 1475—1483 doi: 10.3321/j.issn:0001-5733.2008.05.021Liu B. J. , Chai C. Z. , Feng S. Y. , et al. , 2008. Seismic exploration method for buried fault and its up-breakpoint in Quaternary sediment area-An example of Yinchuan buried active fault. Chinese Journal of Geophysics, 51(5): 1475—1483. (in Chinese) doi: 10.3321/j.issn:0001-5733.2008.05.021 [7] 宋佳佳, 田晓峰, 王帅军等, 2021. 南华北块体及邻区地质结构构造研究进展. 地球物理学进展, 36(3): 940—952 doi: 10.6038/pg2021EE0286Song J. J. , Tian X. F. , Wang S. J. , et al. , 2021. Research progress of geological structure in the South North China and adjacent areas. Progress in Geophysics, 36(3): 940—952. (in Chinese) doi: 10.6038/pg2021EE0286 [8] 王定一, 汤锡元, 陈乃明, 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 Kaifeng depression. Acta Petrolei Sinica, 15(2): 39—47. (in Chinese) doi: 10.3321/j.issn:0253-2697.1994.02.005 [9] 王晶, 李德文, 马保起等, 2021. 黄河郑州—济南段河型变化特征及其与隐伏断层活动的关系. 地质论评, 67(6): 1605—1618Wang J. , Li D. W. , Ma B. Q. , et al. , 2021. River pattern anomalies of Zhengzhou—Jinan segment of Yellow River, North of China, and its relationship with activity of subtle faults. Geological Review, 67(6): 1605—1618. (in Chinese) [10] 王学潮, 向宏发, 2001. 聊城—兰考断裂综合研究及黄河下游河道稳定性分析. 郑州: 黄河水利出版社, 1—163. [11] 王志铄, 马兴全, 2018. 郑州-开封断裂新生代活动特征. 地震地质, 40(3): 511—522 doi: 10.3969/j.issn.0253-4967.2018.03.001Wang Z. S. , Ma X. Q. , 2018. The activity characteristics of Zhengzhou-Kaifeng fault during kainozoic. Seismology and Geology, 40(3): 511—522. (in Chinese) doi: 10.3969/j.issn.0253-4967.2018.03.001 [12] 闫震鹏, 2009. 地质构造对黄河下游悬河稳定性的控制作用. 人民黄河, 31(7): 16—17, 19 doi: 10.3969/j.issn.1000-1379.2009.07.007Yan Z. P. , 2009. Control function of Geologic structure to the stability of elevated section of the Lower Yellow River. Yellow River, 31(7): 16—17, 19. (in Chinese) doi: 10.3969/j.issn.1000-1379.2009.07.007 [13] Liang K. , Sun C. B. , Ma B. Q. , et al. , 2018. Investigation of the Yellow River buried fault in the Wuhai basin, northwestern ordos block, China, using deep/shallow seismic reflection and drilling techniques. Journal of Asian Earth Sciences, 163: 54—69. doi: 10.1016/j.jseaes.2018.05.025 [14] Yu Z. Y. , Pan H. , Xi H. , et al. , 2019. Late Quaternary paleoseismicity of the Xiadian fault in the North China Plain with implications for earthquake potential. Journal of Asian Earth Sciences, 184: 103997. doi: 10.1016/j.jseaes.2019.103997 -