0 引言

土层作为地震发生时地震波的传播媒介，会影响地震动的频谱特征，具体影响为滤波、减震、隔震、共振、放大等［1-3］。在工程场地施工之前，为了解场地的抗震性能等，需要对场地进行土层反应分析，分析结果为场地安全性评价、抗震设计等服务。然而，土是一种非线性很强的介质，地震动在土中的传播很难用公式精确计算［4-5］。在多种假定的基础上，发展较早的是一维土层反应分析理论，随后二维、三维土层反应分析理论也逐步开展。目前应用最广泛的是一维土层反应分析理论［6］。

地震动中的剪切波在传播到地表时被近似为垂直土层向上传播，而剪切波是造成地震动破坏的最主要部分。一维土层反应分析的研究正是建立在这一假设基础上，分析的结果可以代表占能量大多数的剪切波对土层的破坏或影响结果。这一方法中，有两种较为主流的方法：等效线性化方法和时域非线性方法。这两种方法的区别在于对地震动的处理方法上。这两种方法都有代表性的分析程序，其中等效线性化方法的代表性程序为SHAKE2000与LSSRLI-1，这两个程序分别在国外和国内得到广泛的推广，被业界普遍认可［2］。时域非线性方法中比较有代表性的程序是UIUC开发的DEEPSOIL，这一程序界面友好，操作简单，结果可靠。在近期科研成果中，袁晓铭团队公布了其开发的新一代土层反应分析程序SOILQUAKE，这一程序有单机版与网页版两个模式。其研究结果表明，该程序在计算深软场地时与实测结果更为接近［7-8］。在国内，LSSRLI-1与SOILQUAKE作为土层反应分析方法的两代程序，对两者的计算结果进行分析有助于认识到两个程序处理结果的异同。

本文在对比过程中，选用了硬、中硬与软场地三种场地类别的6个场地，分别输入三条地震动进行对比计算，以上组合共得到53组工况。分别用SOILQUAKE与LSSRLI-1对这些工况进行计算，将得到的结果进行对比分析，为今后工程应用与程序改进提供参考。

1 SOILQUAKE简介

SOILQUAKE程序是袁晓铭团队近年攻关研发的新一代一维土层地震反应分析程序，仍采用等效线性化方法模拟土体的非线性，创新之处在于提出直频法，即动剪模量阻尼比采用直接频率计算的方法，且前期与室内试验相结合，给出了直接考虑荷载频率特征的土体动剪模量与阻尼比计算方法，依靠目前先进的计算机计算能力，使更为复杂的直频法得到应用，提高SOILQUAKE程序的计算值的理论精度［6］。

该程序目前推出了网页试用版版与单机试用版两个版本。本文选择其网页版作为研究对象。以下简述SOILQUAKE网络版使用过程。用户首先需在 SOILQUAKE 的官方网站（http：／／www.soilquake.org.cn／）上注册用户，之后上传计算的工况将被保存在这一账号里。计算操作步骤如下：

（1）建立工况：在左侧Information栏目中，用户点击 Add Project建立项目，在项目中点击Addsite建立工况，并在工况中点击EditData开始输入工况的各项参数。

本文采集数据来源于中国知网(CNKI)数据库，采集的时间为2017年12月19日。检索时主题词为“混合式教学”、“小学”、“中学”、“高中”、“初中”、“中小学”，检索的时间区间为2007-2017年，一共检索到81篇文章。

聚合物钻井液配方1#4%BTJ+0.2%QYHN+0.2%TSN+0.4%MAN104+0.4%MAN101+0.5%FPAY+2%SMP-1+1%SPNH+2%KR-n+BA；

（2）输入土的非线性参数文件：在上一步骤完成后，网页直接跳到Material栏目。此处可以选择直接导入文件或手动输入数值两种方式。在数据读取后，点击Plot键，可在网页右侧直接输出土的非线性曲线。最后点击Save保存。

（3）输入土层剖面的物理学参数文件：点击左侧Soil Profile栏目，开始输入每层土厚度、埋深、密度、土类等，输入方式也有导入文件与手动输入两种。输入成功后页面上会直接显示出剖面图。最后点击Save保存。

（4）上传基岩输入地震动时程文件：点击左侧Input Motion栏目，导入时程文件后，在本栏目选择时程文件时间步长，单位、调幅等。输入后页面会直接绘制以上条件下的地震动时程曲线，供用户校对。最后点击Save保存。

硬场地中的两个场地厚度均为3m，分别为单层场地与双层场地，卓越周期分别是0.034s与0.053s。考虑所有工况，SOILQUAKE与LSSRLI-1计算出的地表加速度峰值PGA结果与场地类别的关系如表3所示。

本条第二款对合伙形式的评估机构应有的评估师人数和合伙人资格作出了特殊要求。合伙形式的评估机构聘用评估专业人员开展评估业务，其中应当有两名以上评估师，以保证评估机构的专业水准，并能够开展法定评估业务。该合伙形式评估机构的合伙人三分之二以上应当是具有三年以上从业经历且最近三年内未受停止从业处罚的评估师。评估机构的合伙人为两名的，两名合伙人都应当是具有三年以上从业经历且最近三年内未受停止从业处罚的评估师。

由表3可知，在硬场地下，SOILQUAKE与LSSRLI-1计算的PGA结果相差极小，18个工况的PGA相对差最大为1%，平均值为0，每个工况的PGA放大系数相同。这在一定程度上说明，这两个程序在计算硬场地时，不管输入地震动的频谱分布如何，输出的PGA结果基本没有差异。

2 计算工况

为尽量全面的比较场地条件对结果的影响，本文依据《建筑抗震设计规范》（GB50011—2010）（2016版）构造了硬、中硬与软场地三种场地类型下的6个理想场地，场地模型如表1所示，选取理想场地的原因。其中每个类型分为单层模型与双层模型，剪切波速从90～800m／s不等，卓越周期在0.034～3.91s，厚度分布范围为3～120m，输入地震动强度分别为 0.1g、0.2g和0.4g，以代表三个烈度。计算工况总计53组。

表1 各场地土层模型 Tab.1 The models of soil in different site

剖面 土层类型 密度（kg／m3） 波速（m／s） 厚度（m） 卓越周期（s）硬场地-1 粘土 1640 350 3基岩 2200 800 0.034硬场地-2粘土 1640 300 1粘土 1640 400 2基岩 2200 800 0.053中硬场地-1 粘土 1640 200 24基岩 2200 800 0.48中硬场地-2粘土 1640 150 8粘土 1640 250 16基岩 2200 800 0.47软场地-1 粘土 1640 200 80基岩 2200 800 1.6

续表1

剖面 土层类型 密度（kg／m3） 波速（m／s） 厚度（m） 卓越周期（s）软场地-2粘土 1640 200 26粘土 1640 300 54基岩 2200 800粘土 1640 150 80 1.24基岩 2200 800

为研究SOILQUAKE与LSSRLI-1在对不同频率地震动的处理差异，本文选择了三种频率分布特点的地震动输入，频率分布均匀的El Centro地震动记录来自1940年美国帝王谷地震，低频丰富的AKTH19地震动记录与高频丰富的KSRH地震动记录取自日本Kik-Net台网公布的实测记录。这三个地震动输入在强度为0.2g时的加速度时程如图1所示。

图1 地震动输入的时程与反应谱 Fig.1 Time history curves of inputting ground motion and the corresponding response spectrums

在选取土的非线性参数时，参考孙锐等对国内各类土的非线性研究成果［9-10］，以其中的粘土平均非线性参数作为研究参数。如表2所示。

表2 土的非线性计算参数 Tab.2 Nonlinear parameters of soil

土类 模量比与阻尼比剪应变γ（×10-4）0.05 0.1 0.5 1 5 10 50 100粘土平均G／Gmax 0.994 0.985 0.927 0.867 0.586 0.427 0.195 0.077 λ 0.026 0.032 0.052 0.087 0.139 0.159 0.195 0.203基岩G／Gmax 1 1 1 1 1 1 1 1 λ 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05

3 计算结果分析

3.1 硬场地

（5）选择需要输出的参数：点击左侧Output Option栏目，勾选需要计算的土层与输出的结果即可。最后点击Save保存。

表3 硬场地下两程序PGA相对差对比 Tab.3 Contrast of relative difference of PGA by two programs in hard field

剖面 地震动输入输入PGA（g）输出 PGA（g）PGA放大系数SOILQUAKE LSSRLI-1 SOILQUAKE LSSRLI-1 PGA相对差（%）硬场地-1 El Centro AKTH19 KSRH 0.1 0.1 0.11 0.11 1.13 1.14 0.2 0.2 0.23 0.23 1.13 1.15 0.4 0.4 0.46 0.47 1.14 1.17 0.1 0.1 0.10 0.10 1.01 1.02 0.2 0.2 0.20 0.20 1.01 1.02 0.4 0.4 0.41 0.41 1.01 1.02 0.1 0.1 0.12 0.12 1.17 1.17 0.2 0.2 0.23 0.23 1.17 1.17 0.4 0.4 0.47 0.47 1.17 1.17硬场地-2 El Centro AKTH19 0.1 0.1 0.11 0.11 1.07 1.07 KSRH 0.2 0.2 0.21 0.21 1.07 1.07 0.4 0.4 0.43 0.43 1.07 1.07 0.1 0.1 0.10 0.10 1.01 1.02 0.2 0.2 0.20 0.20 1.01 1.02 0.4 0.4 0.40 0.41 1.01 1.02 0.1 0.1 0.11 0.11 1.13 1.13 0.2 0.2 0.23 0.23 1.13 1.13 0.4 0.4 0.45 0.45 1.13 1.13

（6）上传计算：点击左侧Result栏目，在计算工况栏目的 Operation列中点击 Resquest Calculation，即时开始计算。约一分钟后刷新页面，在工况的State列下选择download，即可下载计算结果到本地。

西格列汀对高脂饮食诱导的肥胖大鼠肾脏保护作用及相关机制研究………………………… 刘晓玲，陈永刚，王 犇（4·245）

频率分布较为均匀的EL Centro波输入下，两个程序输出的反应谱整体走向一致，但在高频部分反应谱形状差异较大，LSSRLI-1最大峰值更大，SOILQUAKE 的数值更为详细，因为SOILQUAKE输出的反应谱周期步长更小。二者 在周期大于0.3s之后反应谱基本重合。

在该公路养护工程之中，养护单位不能够只开展短期养护，需要不断提高长效机制的实施水平，最好将大修养护工程与中修养护工程进行结合，并根据该公路的养护现状，树立更为清晰的公路养护意识。对于公路建设单位来讲，可以采用合同形式，对养护单位进行行为制约，减少不规范公路养护行为的出现，从根源上提升公路大中修工程养护水平。建设部门还要对原有的大中修养护管理体系进行完善与优化，保证公路大中修养护工作得到更好执行与实施[2]。

对不同强度、频谱特征的地震动输入后的地表反应谱进行对比，如图2～4所示。同一条地震动在不同的强度输入下，输出的反应谱形状较为一致。

图2 El Centro波输入下0.1g＼0.2g＼0.4g硬场地下两程序地表反应谱对比 Fig.2 Contrast of response spectrums by two programs in hard field of 0.1g ／0.2g ／0.4g under El Centro wave

图3 AKTH19波输入下0.1g＼0.2g＼0.4g硬场地下两程序地表反应谱对比 Fig.3 Contrast of response spectrums by two programs in hard field of 0.1g ／0.2g ／0.4g under AKTH19 wave

图4 KSRH波输入下0.1g＼0.2g＼0.4g硬场地下两程序地表反应谱对比 Fig.4 Contrast of response spectrums by two programs in hard field of 0.1g ／0.2g ／0.4g under KSRH wave

在低频成分较为丰富的AKTH19波输入下，三种输入强度下两个程序输出的反应谱重合度很高，不同之处在于，三种输入强度下，SOILQUAKE输出的反应谱最大值是LSSRLI-1的1.5倍左右。

下午，齐眉在高潮和田卓的带领下，驱车来到卓木文化传媒。经过一番磋商，达成了以拾伍万圆入选“十大本土最具发展潜力房地产企业”的意向，齐眉代表美之厦公司签署了合作协议，并当场拿出肆万圆现金作为信誉保证金，说余下的拾壹万圆，明天就会打入卓木文化传媒的账户。

在高频成分较为丰富的KSRH波输入下，三种输入强度下两者输出的反应谱重合度很高，SOILQUAKE的反应谱最大值稍大于LSSRLI-1。

那蓝蓝的天，那白白的云，那青青的草，那灿灿的花，衬托着奔腾的骏马。内蒙古民谚说：千里疾风万里霞，追不上百岔的铁蹄马。在边疆少数民族发展史上，蒙古马常常担负着转牧场、踏坚冰、战疆场等重任，只要接受指令，就无所畏惧、勇往直前，它的足迹可谓踏遍北方大草原的每一寸土地。在古代，成吉思汗统帅的蒙古军队之所以能取得非凡的战绩，都跟蒙古人的坐骑——蒙古马有着密不可分的关系。据历史资料描述，经过调教驯化的蒙古马，在战场上不惊不乍，勇猛无比，任何障碍都阻挡不了它的凌厉步伐。可以说无论是在战场上还是在日常劳作中，蒙古马都鲜明地体现了斗志昂扬、奔腾不息的突出特征。

3.2 中硬场地

中硬场地中的两个场地的厚度都为24m，分别为单层场地与双层场地，卓越周期分别是0.48s与0.47s。 将SOILQUAKE与LSSRLI-1计算得到的地表反应谱进行比较，得到表4。

当AKTH19波输入时，两者输出的反应谱形状大致相同，其差异随地震动输入的增大而增大，在高频部分SOILQUAKE的输出反应谱值大于LSSRLI-1的，这一差异也随地震动输入增大而增大，但在低频部分两者结果始终较为重合，仅与2s左右的反应谱峰值有所差异，LSSRLI-1的值略小。

两个程序分别计算中硬场地时，PGA的结果差异开始显现，SOILQUAKE的结果普遍大于LSSRLI-1的，18个工况的平均值为21%。PGA放大系数对比发现，SOILQUAKE的PGA放大系数结果普遍大于LSSRLI-1，6个工况的PGA放大系数平均值为2.27，是LSSRLI-1的PGA放大系数平均值1.54的1.5倍左右。在频率成分分布均匀的El Centro波输入下，这一差异有随输入强度的增大而增大的趋势，其他两个输入这一规律不明显。

如图5～7所示，当 El Centro波输入时，两者输出的反应谱形状一致，在0.1g输入时两者基本重合，0.2g输入时差异开始显现，周期0.5s之前的反应谱值，SOILQUAKE明显大于LSSRLI-1，但二者的反应谱最大值一致。在0.4g输入时，反应谱差异变大，且二者的反应谱最大值的对应周期不一致。SOILQUAKE的反应谱最大值对应的周期与与0.1g和0.2g输入时的结果相一致，为0.5s左右，但LSSRLI-1此时的反应谱最大值出现在1s左右。此外，SOILQUAKE在高频部分的数值式中大于LSSRLI-1的结果，且差异随着地震动输入的加强而加强。

表4 中硬场地下两程序PGA相对差对比 Tab.4 Contrast of relative difference of PGA by two softwares in medium hard field

剖面 地震动输入输入PGA（g）输出PGA（g） PGA放大系数SOILQUAKE LSSRLI-1 SOILQUAKE LSSRLI-1PGA相对差（%）中硬场地-1 El Centro AKTH19 KSRH 0.1 0.26 0.20 2.57 2.03 12%0.2 0.44 0.31 2.19 1.55 17%0.4 0.73 0.24 1.81 0.61 50%0.1 0.24 0.19 2.42 1.95 11%0.2 0.45 0.28 2.26 1.39 24%0.4 0.98 0.66 2.44 1.66 19%0.1 0.18 0.11 1.84 1.10 25%0.2 0.36 0.21 1.80 1.04 27%0.4 0.61 0.28 1.53 0.69 38%中硬场地-2 El Centro AKTH19 0.1 0.31 0.25 3.06 2.47 11%KSRH 0.2 0.62 0.46 3.08 2.32 14%0.4 1.24 0.69 3.10 1.72 29%0.1 0.22 0.22 2.22 2.20 1%0.2 0.43 0.40 2.15 2.02 3%0.4 1.06 0.64 2.65 1.61 24%0.1 0.20 0.13 2.02 1.30 22%0.2 0.41 0.22 2.07 1.10 30%0.4 0.67 0.42 1.68 1.05 23%

图5 El Centro波输入下0.1g＼0.2g＼0.4g中硬场地下两程序地表反应谱对比 Fig.5 Contrast of response spectrums by two programs in medium hard field of 0.1g／0.2g ／0.4g under El Centro wave

图6 AKTH19波输入下0.1g＼0.2g＼0.4g中硬场地下两程序地表反应谱对比 Fig.6 Contrast of response spectrums by two programs in medium hard field of 0.1g ／0.2g ／0.4g under AKTH19 wave

图7 KSRH波输入下0.1g＼0.2g＼0.4g中硬场地下两程序地表反应谱对比 Fig.7 Contrast of response spectrums by two programs in medium hard field of 0.1g ／0.2g ／0.4g under KSRH wave

地形图的定位，可采用以下三种方式：①移动光标；②鼠标定位；③输入坐标。先设置原始图对应的比例尺，再设置绘图比例尺，将默认值确定为1∶10000，距离容限一般为25m，表示线路交点间距不能小于这一数值。

当KSRH波输入时，两者输出的反应谱形状有较大差异，主要体现在高频部分。LSSRLI-1的反应谱中，反应谱最大值出现在周期0.5s到1s之间，而SOILQUAKE的反应谱最大值出现在较为高频的0.2s处。SOILQUAKE的反应谱值在高频部分明显大于LSSRLI-1的，在最后一个明显峰值之后，两条曲线才较为一致。

3.3 软场地

软场地中两个场地厚度都为80m，分别为单层场地与双层场地，卓越周期分别是1.6s与1.24s。将SOILQUAKE与LSSRLI-1计算得到的地表输出PGA与PGA放大系数进行比较，得到表5。

随着场地变得深厚，在软场地时，两程序得到的PGA结果出现了明显差异，SOILQUAKE的结果的PGA放大系数的平均值是1.72，而LSSRLI-1的平均值是0.84，即地表输出PGA小于基岩输入。根据日本强震台网Kik-Net实测结果的统计，软场地下PGA放大系数平均值为4.21［11-12］，而SOILQUAKE的PGA放大系数在软场地6个工况中的平均值为1.72，LSSRLI-1的PGA放大系数的平均值为0.84，SOILQUAKE的放大系数约为LSSRLI-1的2倍。对比表明，SOILQUAKE更体现了场地放大效应，LSSRLI-1计算结果失真较大。

场地从硬到软，SOILQUAKE与LSSRLI-1输出结果的差异变大。SOILQUAKE在计算高频部分时，得到的反应谱谱值结果大于LSSRLI-1。所有场地下，SOILQUAKE的PGA与反应谱峰值普遍大于 LSSRLI-1；随着输入强度变大，SOILQUAKE的结果中的PGA放大系数基本不变，但 LSSRLI-1的PGA放大系数明显变小，SOILQUAKE更能体现出输入强度改变对结果带来的影响。

表5 软场地下两程序PGA相对差对比 Tab.5 Contrast of relative difference of PGA by two softwares in soft field

剖面 地震动输入输入PGA（g）输出PGA（g） PGA放大系数PGA相对差（%）SOILQUAKE LSSRLI-1 SOILQUAKE LSSRLI-1软场地-1 El Centro AKTH19 KSRH 0.1 0.15 0.06 1.46 0.62 40%0.2 0.29 0.09 1.44 0.43 54%0.4 0.47 0.11 1.17 0.27 62%0.1 0.15 0.11 1.51 1.08 16%0.2 0.31 0.15 1.54 0.77 33%0.4 0.58 0.29 1.44 0.74 32%0.1 0.16 0.05 1.57 0.54 49%0.2 0.29 0.06 1.43 0.30 66%软场地-2 El Centro AKTH19 KSRH 0.1 0.22 0.12 2.21 1.23 29%0.2 0.39 0.17 1.97 0.84 40%0.4 0.55 0.22 1.37 0.55 43%0.1 0.26 0.22 2.61 2.20 8%0.2 0.40 0.25 1.99 1.25 23%0.4 0.75 0.52 1.88 1.31 18%0.1 0.19 0.11 1.93 1.05 29%0.2 0.38 0.15 1.91 0.74 44%0.4 0.73 0.15 1.83 0.38 66%

4 结论

SOILQUAKE与LSSRLI-1是目前国内有代表性的两个一维土层反应分析程序。通过对硬场地、中硬场地、软场地在不同频率成分地震动作用下的对比，得到以下结论：

（1）输出结果差异的大小主要与场地类别、输入地震动强度和输入地震动的频率分布有关，场地越软，SOILQUAKE结果的稳定性越明显。

（2）两程序输出的PGA结果差异在硬场地时极小，中硬场地和软场地类场地时存在较大差异，但SOILQUAKE很好的反映了场地放大效应，结果较之LSSRLI-1更贴近实际场地记录统计结果。

（3）地表反应谱差异在硬场地时极小，中硬场地和软场地时存在较大差异。SOILQUAKE在高频部分的结果大大优于LSSRLI-1，在低频部分还能体现软土场地的放大作用。

王俊亮等[24]的试验认为，周氏啮小蜂人工繁育的最佳中间寄主为柞蚕蛹，以秋茧最好，采用“三刀破茧法”削茧接种小蜂，可达到理想效果。用柞蚕蛹繁育周氏啮小蜂时，每个茧内接蜂60～70头雌蜂时，每个柞蚕蛹出蜂量为16 000～18 000头[25]。而罗明的研究认为，以1∶40接蜂比例可得到最大的平均出蜂量和平均羽化量，且蜂生命力较强、寿命较长[20]。徐明等的试验也得出了相似的结论[4]。

（4）在软场地下，SOILQUAKE计算得的PGA放大系数远高于LSSRLI-1，更符合客观统计事实。

参考文献

［1］ 张克绪，谢君斐.土动力学［M］.北京：地震出版社，1989.

［1］ 蔡晓光，范丽远.地震液化引起地面侧向大变形研究评述［J］.防灾科技学院学报，2010，12（1）：11-16.

［3］ 孟凡超，郭鑫.地震动特征的统计分析［J］.防灾科技学院学报，2011，13（3）：21-24.

［4］ 刘德东，齐文浩，张宇东，等.现行土层地震反应分析存在的问题［J］.防灾科技学院学报，2009，11（3）：34-37.

［5］ 袁晓铭，孙锐，孙静，等.常规土类动剪切模量比和阻尼比试验研究 ［J］.地震工程与工程振动，2000， 20（4）：133-139.

［6］ 齐文浩，薄景山.土层地震反应分析的研究现状［J］. 世界地震工程，2010， 26（S1）：368-372.

［7］ 李兆焱，袁晓铭，王鸾，等.巨厚场地三种土层地震反应分析程序对比检验［J］.地震工程与工程振动，2017，37（4）：42-50.

［8］ 袁晓铭，李瑞山，孙锐.新一代土层地震反应分析方法［J］. 土木工程学报，2016，49（10）：95-102，122.

［9］ 陈红娟.土动力非线性的变异性及其对地震动影响的概率分析［D］.哈尔滨：中国地震局工程力学研究所，2009.

［10］ 孙锐，陈红娟，袁晓铭.土的非线性动剪切模量比和阻尼比不确定性分析 ［J］.岩土工程学报，2010，32 （8）：1228-1235.

［11］ 王亮.基于Kik-Net强震台网的土层地震动特性研究［D］.哈尔滨：中国地震局工程力学研究所，2014.

［12］ 陈龙伟，陈卓识，袁晓铭.基于Kik-Net单场地震记录场地放大函数估计及标准差分析 ［J］.土木工程学报，2013，46（S2）：41-145.