A new method of modified frequency decrease rate for site liquefaction detection
-
摘要: 基于地表强震记录的场地液化实时监测和震后快速识别是液化减灾的新技术手段,其核心为场地液化识别方法。现有此类方法受限于有限的液化场地记录和方法自身局限,适用性和识别准确率有待提升。本研究建立了包含23次地震共137条加速度记录的液化场地强震记录新数据库,极大地扩充了现有数据库。提出了一种基于强震记录的场地液化快速识别改进的频率下降率新方法,以频率下降率为基本概念并沿用以往推导出的理论阈值。采用加权平均方式处理短时傅里叶变换结果以获取加速度时程的时间-频率曲线。通过频率下降率参数的自适应持时计算,以增强对不同持时记录的适应性。以新数据库记录为基础进行对比分析,结果表明新方法所选取的特征参数明显占优,对液化和非液化场地均有很高的识别成功率,且较新近提出方法的液化识别成功率显著提高。所提出方法在震后液化场地快速识别、液化实时减灾及研究液化对建筑结构的影响等方面具有良好的应用前景。Abstract: The real-time monitoring and rapid post-earthquake detection of site liquefaction based on strong motion records is a new technological tool for liquefaction mitigation, the core of which is the site liquefaction detection method. The existing methods are constrained by the limited liquefaction site records and the limitations of the methods themselves, and their applicability and detection accuracy need to be improved. In this study, a new database of strong motion records for liquefaction sites containing 137 records from 23 earthquakes is established, which greatly expands the existing database. A new method for rapid detection of site liquefaction based on strong motion records with improved frequency decrease rate is proposed, taking the frequency decrease rate as the basic concept and following the theoretical thresholds derived in the past. The weighted average method is used to process the results of short-time Fourier transform to obtain the time-frequency curve of acceleration time history. The adaptive duration calculation of frequency decrease rate parameter is used to enhance the adaptability to the records of different durations. The comparative analysis based on the new database records shows that the feature parameter selected by the new method is obviously superior, with high detection success rates for both liquefied and non-liquefied sites, and the liquefaction detection success rate is significantly higher than that of the newly proposed methods. The proposed method has good application prospects in rapid detection of liquefied sites after earthquakes, real-time mitigation of liquefaction, and study of liquefaction effects on building structures.
-
0. 引言
海洋复合地基施工常存在定位难、制桩难、桩体质量差等缺点,始终制约着海洋基础设施工程的应用和发展[1-3]。本文借鉴在陆地上已成功应用于湿陷性黄土、盐渍土、膨胀土、红层风化土等特殊土的地基处理技术[4-6]——孔内深层强夯法(DDC),提出水下套管强夯碎石桩地基处理方法,施工工序为:先将钢套管打至设计深度并清理管内软泥;套管内分层填入碎石并使用夯锤夯实;上拔套筒一定高度并再次填料和夯实;直到形成串珠状的散体桩复合地基。该过程不仅能将碎石夯实压密,而且可以对海洋地基土体进行置换、挤密和排水,最终形成高强度的水下散体桩复合地基。该技术是将陆上孔内深层强夯桩复合地基技术进行改进,增加钢套筒隔断海水以实现海洋干作业施工,满足海洋软土地基处理的需求。然而在成桩过程中涉及夯锤-碎石-土体三者相互作用,目前还缺乏相应的研究。
孔内深层强夯法(DDC)的加固效果、承载性能、挤密影响范围、布桩形式等受到了广泛研究。孔洋等[7]对DDC法复合地基进行原位浸水载荷试验研究,发现DDC工法可以有效消除黄土的湿陷性。周小松等[8]得到DDC最优桩桩间距为2.0~3.0倍成孔直径。马军平[9]通过井探开挖发现DDC桩在宏观上形成葫芦状的桩体,这种不规则形状极大地提高了侧摩阻力,使复合地基承载力整体提高。冯志焱等[10]发现DDC处理后的复合地基可达到处理前土层承载力的2倍或更高。对于直接影响DDC后期承载特性的成桩过程的研究很少。
本文采用半模试验[11-12]针对水下套筒强夯技术开展研究,讨论了夯击能、单次填料量、土体强度等因素对水下套筒强夯碎石桩成桩过程的影响,分析了夯锤运动过程中的夯锤冲击力和夯沉量、桩体密度等的变化规律,为水下套管强夯碎石桩的工程应用提供参考。
1. 模型试验简介
1.1 试验装置
试验装置如图 1所示,由模型箱、套管、夯击架、夯锤组成。模型箱四周由透明钢化玻璃组成,尺寸为80 cm×80 cm×80 cm(长×宽×高)。半模试验套管为“D”字型,材质为透明有机玻璃,长度为130 cm,内径为5.4 cm,外径为6 cm。夯击架高度为220 cm,宽度为100 cm。定滑轮固定在可以平移的导杆上,牵引线穿过定滑轮,向下用力可将夯锤吊起。夯击过程中为保证套管位置固定,套管通过固定夹固定在固定杆上。夯锤上部固定通用压电式加速度传感器。为减小夯锤与套筒之间的阻力,每次夯击时在夯锤表面涂抹凡士林。夯击过程中通过动态数据采集仪持续监测夯锤加速度变化,动态数据采集仪型号为DH5922D,采样频率5 kHz。
砂垫层作为持力层和固结阶段排水层,厚度为5 cm。固结后软土层厚度为54 cm。为模拟水下环境,软土上部有5 cm的水。土体固结完成后,采用荷兰钻进行钻孔,孔深为30 cm,孔径为6 cm。随后放置套管、填料、进行夯击。夯击方式为每填料一次,夯击一次。填料时一旦填料高度到达套管底部就开始上拔套管,保证套管高于碎石。
1.2 夯锤模型
两种夯锤模型如图 2所示,质量均为1.5 kg。其中图 2(a)为传统平底夯锤,在陆地强夯[13-14]中经常使用。由于平底夯锤在水下套管强夯碎石桩中夯锤冲击会使碎石卡在夯锤与钢套筒之间,导致夯锤无法顺利拔出。因此对传统平底夯锤进行改进,设计了锥形夯锤,即图 2(b)所示的锥形夯锤,该夯锤底部椭球外层有一圈圆形护壁,可保证每次夯击时夯锤顺利拔出。
1.3 试验材料
本模型试验土样为淤泥质软土,取自天津市滨海吹填土区域。共采用3种不同强度的软土进行试验,首先配置含水率60%的淤泥将其充分搅拌,然后分层填入模型箱中。其中软土1在自重条件下固结,固结后的土体强度为1.5 kPa。软土2与软土3在堆载条件下固结,其中软土2的固结压力为18 kPa,固结后的土体强度为6.3 kPa;软土3的固结压力为50 kPa,固结后的土体强度为10 kPa。三种软土固结后通过渗透试验测得渗透系数分别为1.25×10-6 cm/s、1.54×10-7 cm/s、7.77×10-8 cm/s。土体详细物理参数如表 1所示。碎石粒径为6~8 mm,平均粒径为6.95 mm,颗粒相对质量密度为2.58,通过直剪试验测得内摩擦角为38.6°。
表 1 土体物理力学参数Table 1. Physical and mechanical parameters of soils参数 软土1 软土2 软土3 土粒相对质量密度 2.67 2.67 2.67 液限/% 43 43 43 塑限/% 22 22 22 塑性指数 21 21 21 含水率/% 55.90 45.40 41.00 密度/(g⋅cm−3) 1.60 1.75 1.88 不排水抗剪强度/kPa 1.5 6.3 10.0 渗透系数/(cm⋅s−1) 1.25×10-6 1.54×10-7 7.77×10-8 1.4 试验方案
共进行14组水下套筒强夯桩成桩模型试验,试验分组如表 2所示。其中T1作为对照组,其余13组试验分别研究了夯锤落距(H)、单次填料量(m)、土体强度(cu)等因素对水下套管强夯碎石桩夯锤运动的影响。夯锤落距分别为40,80,120 cm,代表 12,24,36 N·m三种夯击能(E)。单次填料量采用50,100,150 g三种质量。土体强度分别为1.5,6.3,10.0 kPa。
表 2 试验分组Table 2. Test groups编号 夯锤类型 H/
cmm/
gcu/
kPaE/
(N·m)夯击次数 T1 平底夯锤 80 100 6.3 24 19 T2 改进锥形夯锤 80 50 6.3 24 50 T3 80 100 6.3 24 18 T4 80 150 6.3 24 10 T5 40 50 6.3 12 35 T6 40 100 6.3 12 12 T7 40 150 6.3 12 7 T8 120 100 6.3 36 28 T9 40 100 1.5 12 30 T10 80 100 1.5 24 50 T11 120 100 1.5 36 66 T12 40 100 10 12 10 T13 80 100 10 24 13 T14 120 100 10 36 16 2. 试验结果和分析
2.1 夯锤运动过程分析
以T3中第10次夯击为例对锥形夯锤夯击过程进行分析,如图 3所示,加速度时程曲线呈单峰型,在夯击早期阶段,夯锤的加速度曲线出现明显的震荡,这反映了夯锤与碎石的压缩接触。夯击过程可分为以下4个阶段:
(1)夯锤下落阶段OA(0.13~0.565 s,持续时间435 ms)。由于夯锤与套筒之间存在摩擦,因此加速度值小于重力加速度,加速度均值为-6.01 m/s2(规定加速度、速度向上为正)。
(2)锤尖接触碎石阶段AB(0.565~0.576 s,持续时间11 ms)。夯锤尖部首先与碎石接触,插入碎石中,表现为加速度的小幅震荡。
(3)震荡阶段BD(0.576~0.596 s,持续时间20 ms)。此时夯锤尖部已没入碎石中,在此阶段夯锤与碎石群体发生多体碰撞,实测加速度发生高频震荡,并产生反力加速度峰值。夯锤的撞击力通过碎石传递给软土,造成软土密实。此阶段可分为震荡上行区和震荡下行区,震荡上行区BC加速度均值逐渐增加,在C点达到加速度峰值,震荡下行区CD加速度均值逐渐减小。
(4)快速下降阶段DE(0.596~0.618 s,持续时间22 ms)。此时碎石和软土都得到密实并伴随夯锤能量释放完毕,加速度曲线快速下降。
图 4为T1中第10次夯击的平底夯锤加速度-时程曲线,可以看出平底夯锤的加速度-时程曲线更加光滑平整,可分为三个阶段:夯锤下落阶段OA、冲击碎石阶段AC、快速下降阶段CD。由于夯锤是平底,因此冲击碎石阶段加速度震荡幅度没有锥形夯锤大,在夯击开始即产生反力加速度峰值,随后反力加速度逐渐减小并小幅震荡,而且BC段加速度存在近似平台区。图 4所示T1试验采用平底夯锤的反力峰值与图 3所示相同条件下的采用锥形夯锤得到的反力峰值接近,但反力的作用时间较长。
2.2 夯锤加速度峰值
以夯锤每次下落到筒底碎石上时受到的反力加速度峰值作为衡量夯锤冲击力的指标。图 5统计了采用锥形夯锤,用不同落距、在不同土体强度中下落夯击时得到的反力加速度峰值。图 5可以看出,随着夯击次数的增加,反力加速度峰值逐渐增大,但是增长速度越来越慢,这主要是由于随着填料量和夯击次数的增加,碎石层厚度增大,密度和刚度增大,下部越来越密实,夯锤和碎石相互作用时间变短,因此加速度峰值逐渐增大。在相同强度的黏土地基中,夯锤的释放落距越大,加速度峰值越大。对比图 5中在不同强度土体中的加速度峰值可知,土体强度越大时,加速度峰值越大。土体强度越高时,限制了碎石的侧向膨胀,减小了碎石的沉降量,缩短了夯锤的能量释放时间,因此加速度峰值越高。
图 6对比了不同的单次填料量对夯锤的反力加速度峰值的影响,从图中可以看出,随着单次填料量的增加,夯锤加速度峰值增加,总的夯击次数减小。在T5中采用落距40 cm和单次填料量50 g时,从第14次夯击开始,加速度峰值不再增加,稳定在160 m/s2左右,在T2中采用落距80 cm和单次填料量50 g时,从第20次夯击开始,加速度峰值不再增加,稳定在220 m/s2左右。这主要是由于单次填料量过少,导致填入的碎石绝大部分被侧向挤压,碎石桩整体刚度没有增加。
总结上述夯击次数(N)、夯击能(E)、土体强度(cu)、孔径(d)、单次填料质量(m)等参数对夯锤受到的反力加速度峰值(即夯锤冲击力)的影响,绘制了归一化的最大冲击力随归一化的夯击次数、夯击能、地基土体强度等的发展曲线,如图 7所示。纵坐标为归一化最大冲击力Fmax/mg,Fmax=mhamax,mh为夯锤质量,amax为前述图 3与图 4中的夯锤反力加速度峰值,g为重力加速度;横坐标为归一化夯击系数NEcud/(mg)2。可以看出两者之间具有良好的对数关系;归一化的最大冲击力随着归一化的夯击次数增长而逐渐增加,但增长速度逐渐变慢。根据图 7统计的归一化最大冲击力Fmax/mg与归一化夯击系数NEcud/(mg)2的结果,得到描述二者关系的拟合公式(1),其R2值为0.9005。利用式(1)可以对夯锤的最大冲击力进行预测。
Fmax (1) 2.3 夯锤夯沉量
通过对夯锤的加速度进行二次积分,可得到每一次夯击时夯锤的夯沉量,采用不同落距在不同强度土体中的碎石沉降量如图 8所示,可以看出随着夯击次数的增加,碎石沉降量逐渐减小,前5次夯击的沉降量远大于后续夯击时的碎石沉降量,这主要是由于碎石桩底部有24 cm厚的软土层,在夯锤的作用下,碎石主要被挤入下部软土层中。夯击能量越大,碎石沉降量越大;土体强度越低,碎石沉降量越大。以T10为例,土体强度为1.5 kPa,落距为120 cm,碎石单次夯沉量均超过了10 cm,这说明在T10的较软的地基中采用过大的夯沉量,使得大部分碎石被挤入到周围软土中,不能形成有效的强夯碎石桩。因此在强度较低的土体中进行强夯碎石桩施工,应采用低落距强夯或者提高土体强度后再进行强夯碎石桩施工。
图 9为单次填料量对碎石沉降量的影响,从图 8可以看出,随着单次填料量的增加,碎石沉降量逐渐减小。当单次填充量过小时,碎石单次沉降量减小到某一值后保持稳定,在T2中,单次填料量为50 g,落距为80 cm,从第11次夯击后单次夯沉量约6 cm。在T5中,单次填料量为50 g,落距为40 cm,从第25次夯击后单次夯沉量约3 cm。
2.4 桩体密度
图 10给出了不同夯击能对碎石桩密度的影响,对于3种不同强度的软土,采用夯击能24 N⋅m(落距80 cm)时强夯碎石桩密度最大,过大或过小的夯击能都不能有效地增加强夯碎石桩密度。采用夯击能为24 N⋅m和36 N⋅m时在10 kPa的土体中强夯碎石桩密度最高,总填料量最少,因此在高强度软土中进行套管强夯碎石桩安装更经济有效。
单次填料量对碎石桩密度的影响如图 11所示,从图 11可以看出,单次填料量也对强夯碎石桩密度产生重要影响,随着单次填料量的增加,碎石桩密度逐渐减小。但是当单次填料量过小时,增大夯锤落距,碎石桩密度并无明显提高,在T5中,填料量为50 g,落距为40 cm,强夯碎石桩密度为1553 kg/m3,当落距增加到80 cm时,在T2中强夯碎石桩密度为1570 kg/m3,并无明显的增加,因此当单次填料量过小时,夯锤采用低落距更为经济合适。
3. 结论
本文针对水下套筒强夯技术开展了模型试验研究,对比分析了传统平底夯锤和新型锥形夯锤的夯击规律,研究了夯击能、单次填料量、土体强度等因素对成桩过程的影响,得出以下结论:
(1)采用新型锥形夯锤夯击到碎石后,得到的反力加速度呈现出震荡增大到峰值后迅速减小、反力作用时间短等特点;采用传统平底夯锤夯击到碎石后,得到的反力加速度表现为迅速增大到峰值再震荡减小、反力作用时间长等特点。相同条件下二者的反力加速度峰值较为接近。
(2)建立了夯锤夯击时的归一化最大冲击力关于夯击次数、夯击能、土体强度、单次填料质量等参数的拟合公式,可用于估算不同地基强度中夯击成桩时的峰值夯击力。随着夯击次数增大,夯锤的反力加速度峰值逐渐增大,但是增长速度放缓。
(3)夯击能量、单次填料量、地基强度等都对成桩质量有影响:夯击能过大将使碎石填料被挤到周围软土中降低碎石桩密度和强度、夯击能过小则无法起到理想的夯击密实效果;随着单次填料量增大,碎石桩密度逐渐减小。
-
表 1 新的强震记录数据库
Table 1 New database of strong motion records
NO 地震事件 时间 经度 纬度 Mw 总数 液化 引用 1 Tokachi-Oki, Japan 16/05/1968 39.84 143.20 8.2 3 1 GREENFIELD M W [8] 2 Friuli, Italy 05/05/1976 46.38 12.98 6.5 1 0 YUAN X M等[5] 3 Tabas, Iran 16/09/1978 33.58 56.92 7.4 2 0 YUAN X M等[5] 4 Imperial Valley, US 15/10/1979 32.69 -115.34 6.5 1 1 GREENFIELD M W[8] 5 Victoria, Mexico 09/06/1980 32.48 -115.24 6.4 2 0 YUAN X M等[5] 6 Nihonkai-Chibu, Japan 26/05/1983 40.38 139.17 7.7 1 0 KOSTADINOV M V等[3] 7 Superstition Hills, US 24/11/1987 33.02 −115.83 6.6 4 1 KOSTADINOV M V等[3] 8 Chibaken-Toho-Oki, Japan 17/12/1987 35.45 140.38 6.5 1 0 KOSTADINOV M V等[3] 9 Loma Prieta, US 18/10/1989 37.04 −121.88 7.0 14 2 KOSTADINOV M V等[3] 10 Kushiro-Oki, Japan 15/01/1993 42.98 144.34 7.6 5 1 KOSTADINOV M V等[3] 11 Hokkaido-Nansei-Oki, Japan 12/07/1993 42.83 139.23 7.7 2 0 KOSTADINOV M V等[3] 12 Northbridge, US 17/01/1994 34.21 −118.55 6.7 6 0 KOSTADINOV M V等[3] 13 Hokkaido-Toho-Oki, Japan 04/10/1994 43.75 147.29 8.2 2 0 KOSTADINOV M V等[3] 14 Sanriku-Haruka-Oki, Japan 18/12/1994 40.52 143.42 7.7 1 0 KOSTADINOV M V等[3] 15 Kobe, Japan 16/01/1995 34.73 135.24 6.9 12 5 KOSTADINOV M V等[3] 16 Kocaeli, Turkey 17/08/1999 40.98 28.69 7.4 5 0 YUAN X M等[5] 17 Chi-Chi, China 20/09, 1999 23.58 120.55 7.6 4 0 YUAN X M等[5] 18 Western Tottori, Japan 06/10/2000 35.27 133.35 6.6 1 1 GREENFIELD M W[8] 19 Tokachi-Oki, Japan 26/09/2003 41.78 144.08 8.3 1 1 GREENFIELD M W[8] 20 NiigataChuetsu, Japan 16/07/2007 37.56 138.61 6.6 1 1 GREENFIELD M W[8] 21 Darfield, New Zealand 03/09/2010 −43.53 172.17 7.1 33 2 ZHANG W W等[4] 22 Christchurh, New Zealand 11/02/2011 −43.58 172.68 6.2 21 9 ZHANG W W等[4] 23 Tohoku, Japan 11/03/2011 38.10 142.86 9.0 14 5 COX B R等[7] 表 2 敏感区间取值
Table 2 Values of sensitivity interval
加速度时程90%能量持时/s 敏感区间 (0, 15) 峰值加速度前后7 s [15, 30) 峰值加速度前后10 s [30, 50) 峰值加速度前后15 s [50, 70) 峰值加速度前后30 s [70, 90) 峰值加速度前后40 s [90, ~) 峰值加速度前后50 s 表 3 本文方法、ZC方法和YS方法的识别成功率对比
Table 3 Comparison of detection success rates among proposed method, ZC method and YS method
识别成功率 YS法 ZC法 本文方法 液化成功功率 66.67% 76.67% 93.30% 非液化成功率 85.98% 91.59% 86.90% 总体成功率 81.75% 88.32% 88.32% -
[1] SUZUKI T, SHIMIZU Y, NAKAYAMA W. Characteristics of strong motion records at the liquefied sites and judgment for liquefaction[C]//Proceedings of 11th European Conference on Earthquake Engineering, Paris, 1998.
[2] MIYAJIMA M, KITAURA M, NOZU S. Detective method of liquefaction using strong ground motion records[C]// Proceedings of 3rd China-Japan-US Trilateral Symposium on Lifeline Earthquake Engineering, Kunming, 1998: 133-140.
[3] KOSTADINOV M V, YAMAZAKI F. Detection of soil liquefaction from strong motion records[J]. Earthquake Engineering & Structural Dynamics, 2001, 30(2): 173-193.
[4] ZHANG W W, CHEN Q S. Accelerogram-based method for quick assessment of liquefaction occurrence[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2021, 147(8): 1-13. http://www.researchgate.net/publication/353623607_Accelerogram-Based_Method_for_Quick_Assessment_of_Liquefaction_Occurrence
[5] YUAN X M, SUN R, CHEN L W, et al. A method for detecting site liquefaction by seismic records[J]. Soil Dynamics and Earthquake Engineering (1984), 2010, 30(4): 270-279. doi: 10.1016/j.soildyn.2009.12.003
[6] 孙锐, 唐福辉, 陈龙伟, 等. 现有液化识别方法对比分析[J]. 岩土力学, 2011, 32(增刊1): 20-26. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2011S1005.htm SUN Rui, TANG Fuhui, CHEN Longwei, et al. Comparative analysis of existing liquefaction detection methods[J]. Rock and Soil Mechanics, 2011, 32(S1): 20-26. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2011S1005.htm
[7] COX B R, BOULANGER R W, TOKIMATSU K, et al. Liquefaction at strong motion stations and in Urayasu City during the 2011 Tohoku-Oki earthquake[J]. Earthquake Spectra, 2013, 29(S1): 55-80.
[8] GREENFIELD M W. Effects of Long-Duration Ground Motions on Liquefaction Hazards[D]. Ann Arbor: ProQuest Dissertations & Theses, 2017.
[9] 王倩. 水平地震动持时的特征研究[D]. 哈尔滨: 中国地震局工程力学研究所, 2015. WANG Qian. Study on the Characteristics of Horizontal Ground Motion Duration[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administration, 2015. (in Chinese)
[10] 周健, 李小军, 李亚琦, 等. 中美建筑抗震设计规范中工程场地类别的对比和换算关系[J]. 地震学报, 2021, 43(4): 521-532. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXB202104011.htm ZHOU Jian, LI Xiaojun, LI Yaqi, et al. Comparative analysis and transformation relations between China and the US site classification systems in building seismic code provisions[J]. Acta Seismologica Sinica, 2021, 43(4): 521-532. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DZXB202104011.htm
[11] KRAMER S L, SIDERAS S S, GREENFIELD M W. The timing of liquefaction and its utility in liquefaction hazard evaluation[J]. Soil Dynamics and Earthquake Engineering, 2016, 91: 133-146. doi: 10.1016/j.soildyn.2016.07.025
[12] 徐国林, 陈龙伟. 液化场地强震记录的频谱特征[J]. 地震学报, 2019, 41(5): 671-679. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXB201905009.htm XU Guolin, CHEN Longwei. Frequency spectra characteristics of strong-motion acceleration recordings at liquefied sites[J]. Acta Seismologica Sinica, 2019, 41(5): 671-679. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DZXB201905009.htm
[13] FILIPPO S D M, LANZANO G, FORTE G, et al. A database for PGA threshold in liquefaction occurrence[J]. Soil Dynamics and Earthquake Engineering (1984), 2013, 54: 17-19. doi: 10.1016/j.soildyn.2013.07.011
[14] QIAN S, CHEN D P. Joint time-frequency analysis[J]. IEEE Signal Processing Magazine, 1999, 16(2): 52-67. doi: 10.1109/79.752051
[15] KRAMER S L. Geotechnical Earthquake Engineering[M]. Upper Saddle River, NJ: Prentice Hall, 1996.
[16] LOH W, SHIH Y. Split selection methods for classification trees[J]. Statistica Sinica, 1997, 7(4): 815-840.
[17] LOH W Y. Regression trees with unbiased variable selection and interaction detection[J]. Statistica Sinica, 2002, 12(2): 361-386.
[18] 孙锐, 唐福辉, 袁晓铭. 频率下降率法对2011年新西兰6.3级地震液化的盲测[J]. 岩土力学, 2011, 32(增刊2): 383-388. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2011S2064.htm SUN Rui, TANG Fuhui, YUAN Xiaoming. Blind detection of liquefaction sites by the method of frequency decrease rate for New Zealand earthquake[J]. Rock and Soil Mechanics, 2011, 32(S2): 383-388. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2011S2064.htm