典型土层场地随机地震反应规律分析1

杜修力 袁雪纯 黄景琦 许紫刚



典型土层场地随机地震反应规律分析1

杜修力 袁雪纯 黄景琦 许紫刚

（北京工业大学，城市防灾减灾教育部重点实验室，北京 100124）

针对硬、中、软3种土层场地，选取历史上实测到的Ⅰ类和Ⅱ类场地的地震动记录各100条，分别调整地震动记录加速度峰值至0.1g、0.2g和0.3g，并采用一维等效线性化方法开展了场地随机地震反应研究，系统分析了地震动峰值加速度、速度和位移反应的变异性规律。主要结论为，位移和速度峰值的变异性随场地土变软而增大；位移和速度峰值变异性较加速度峰值的变异性更为突出；利用50条左右的地震动记录即可获得场地地震反应均值和标准差较为稳定的结果。

典型土层场地 随机地震反应 变异性规律

引言

场地地震反应分析是场地地震安全性评价和工程结构抗震设计的重要环节。影响场地地震反应的因素有很多，如地震动参数、场地条件等。地震动参数是表征抗震设防要求的地震动物理参数，包括地震动峰值加速度、频谱和持时等（胡聿贤，1988），是工程场地抗震安全性和经济性评价的重要依据。1928年，美国人Wood对1906年美国旧金山大地震的震害进行了详细的调查研究，并认识到场地条件对建筑物震害有影响，成为较早开展场地分析的学者（薄景山，2003）。截至目前，对于场地地震反应已经展开了大量研究，并取得丰富成果（欧阳行艳等，2008；Trifunac，2016）。土层地震反应分析可大致分为确定性分析和不确定性分析2种（齐文浩，2004；齐文浩等，2010）。确定性分析方法是当前场地分析中较为常用的方法，其又可分为等效线性方法和非线性方法。等效线性方法最早由Idriss和Seed提出（Idriss等，1968），我国学者廖振鹏、李小军等对等效线性化方法进行了研究和改进（廖振鹏等，1989；李小军，1989），并提出相应的计算程序。非线性方法主要考虑土的非线性，Kondner、Hardin和Drnevich提出了双曲线模型来描绘土体的骨架曲线（Kondner，1963；Hardin等，1972）。不确定性分析方法是将随机理论应用在场地反应分析研究中，针对地震动参数和场地土层参数的随机性对地震反应影响的研究也开展了许多工作。李天、赵松戈等研究了土层参数的随机性对场地地震反应的影响（李天等，1994；赵松戈等，2000）；李建亮、刘方成研究了场地条件，尤其是软土场地对场地地震反应的影响（李建亮等，2011；刘方成等，2015）兰景岩研究了地震动强度及频谱特征对场地地震反应分析结果的影响（兰景岩等，2012）；陈龙伟选取单台强震数据研究了地震动参数对场地放大函数不确定性的影响（陈龙伟等，2015）。

但针对地震动参数随机性对场地的地震动峰值加速度、峰值速度和峰值位移反应的变异性及其差异的认识仍不深入。本文结合北京金安桥、珠市口和苏州星海广场等地铁车站（分别为硬、中、软场地）场地，选取历史上实测到的Ⅰ类和Ⅱ类场地的地震动记录各100条，并分别调整地震动记录加速度峰值至0.1g、0.2g和0.3g，采用一维等效线性化方法开展了场地随机地震反应研究，系统分析了地震动峰值加速度、峰值速度和峰值位移反应的变异性及其差异。

1 等效线性方法简介

1968年，Idriss和Seed把土近似看成粘弹性材料，首次提出等效线性化的概念来进行土层地震反应分析（Idriss等，1968）。等效线性方法的基本思路是用等效剪切模量和阻尼比代替不同应变振幅下的剪切模量和阻尼比，将非线性问题转化为线性问题来近似求解土的动力问题，由于其计算简便、原理简单、结果稳定、易被技术人员掌握，目前仍然是土层地震反应分析的常用方法（齐文浩，2004）。Shake、LSSRLI-1和EERA等程序是目前较为常用的等效线性化程序（Idriss等，1992；李小军，1989；Bardet等，2000）。其具体实现步骤如下：

（5）重复步骤（2）—（4），直到剪切模量和阻尼比的计算值在两次迭代之后小于所有土层的预定值。

本文采用了EERA程序进行场地的地震反应分析。

2 场地计算模型及参数

2.1 场地资料

本文选取的3种典型的地铁场地条件分别是北京地铁6号线金安桥站、地铁7号线珠市口站和苏州1号线星海广场站，建立一维波动模型进行分析。金安桥站场地覆盖层较浅，覆盖层向下为圆砾地层直至坚硬岩石，场地条件可划归为较为坚硬地铁场地。星海广场站位于典型的深厚软土地层中，土层厚度近60m，且土质较软，场地条件可代表软弱深厚土层场地。珠市口地铁站场地条件介于两者之间，可视为中等场地。表1分别给出金安桥、珠市口和星海广场车站的场地参数。图1、图2、图3分别给出3类场地土体的强度和阻尼随应变变化的曲线，图中实线为剪切模量，虚线为阻尼比。

表1 场地土层剖面材料

图1 金安桥地铁车站场地土体的强度和阻尼比随应变的变化曲线

图2 珠市口地铁车站场地土体的强度和阻尼比随应变的变化曲线

图3 星海广场地铁车站场地土体的强度和阻尼比随应变的变化曲线

由于“共振”效应，场地对与场地动力特性相近的地震动能量分量有放大作用，导致场地地震反应中某些频谱的分量比较突出或显得“卓越”，其对应频谱的周期称为地震卓越周期（陈永新等，2016）。加权平均波速法和子层周期求和法是较为常用的计算场地卓越周期的方法，式（1）为加权平均波速法，式（2）为子层周期求和法。

加权平均波速法计算值比真值偏小20%左右，子层周期求和法比真值偏大20%左右，故文献（童广才等，2000）将2种方法的计算值求平均，以减小误差。经计算，金安桥站、珠市口站和星海广场站场地卓越周期分别为0.32s、0.49s和1.19s。

2.2 基岩输入地震动的选取

我国《建筑抗震设计规范（GB 50011—2010）》（中华人民共和国住房和城乡建设部等，2010）根据土层等效剪切波速和场地覆盖层厚度将场地划分为4类。表2为场地分类，土的等效剪切波速按下列公式计算：

表2《建筑抗震设计规范（GB 50011—2010）》场地分类

Table 2 Site classification of the code for seismic design of buildings (GB 50011—2010)

美国的抗震设计规范以地表以下30m范围内的土层平均剪切波速s30作为场地分类指标，文献（吕红山等，2007）对比中美两国场地分类指标，通过对场地土层波速测试资料的分析，指出按美国规范得到的s30在510m/s以上的场地对应于中国规范的Ⅰ类场地；s30为260—510m/s的场地对应于中国规范的Ⅱ类场地；s30为150—260m/s的场地对应于中国规范的Ⅲ类场地；s30在150m/s以下的场地对应于中国规范的Ⅳ类场地（吕红山等，2007）。根据此关系，本文在PEER Ground Motion Database（PEER地震动数据库）中选取Ⅰ类场地的地震动和Ⅱ类场地的地震动各100条作为场地输入地震动。图4为两类场地的地震动加速度反应谱。图5为两类场地地震动的震级和震源距分布图。

图4 基岩输入地震动加速度反应谱

图5 Ⅰ类和Ⅱ类场地地震动的震级和震源距分布图

反应谱的特征周期反映了地震动的频谱特性，是抗震设计研究的重要参数。反应谱的特征周期是规准化的加速度反应谱曲线开始下降点所对应的周期值（中华人民共和国住房和城乡建设部，2011）。

图6为两组地震动特征周期频数分布图，可以看出Ⅰ类场地的地震动反应谱特征周期多集中在0.1—0.3s，其中特征周期0.15s，占15%，是统计上数量最多特征周期值，其他特征周期分布较为均匀。Ⅱ类场地的地震动反应谱特征周期多集中在0.1—0.4s，其中特征周期0.2s和0.25s分别占19%和18%，其他特征周期分布较为均匀。从统计上看，Ⅱ类场地的地震动反应谱特征明显比Ⅰ类场地的地震动反应谱特征周期长。

图6 不同场地的地震动反应谱特征周期频数分布图

3 典型场地地震随机反应分析

为考虑地震动强度对场地随机地震反应的影响，将输入地震动峰值分别调整为0.1g、0.2g和0.3g，共计1800个算例，下文介绍计算结果。

3.1 加速度峰值

以地表加速度绝对值最大处的时刻为基准，统计该时刻场地沿深度方向的加速度，并求其绝对值的平均值。图7为深度方向上加速度放大系数，从图中可以看出，对于同一类场地，场地对加速度都具有放大作用，但场地加速度峰值放大系数随着输入地震动峰值的增加而减小，这是由于随着输入地震动峰值增大，场地地震反应强烈，土的剪应变增大，剪切模量减小而引起。对比不同场地的地震反应可以发现，金安桥地铁场地对加速度的放大最为明显，星海广场地铁场地放大效应最小，珠市口地铁场地介于两者之间。从中可以看到，场地土越硬，其对地震动的放大作用越大。图8表示地表加速度峰值的均值和变异性，从图中可以看出，在两类场地的地震动记录输入条件下，峰值加速度反应的均值基本一致。对于金安桥地铁场地和珠市口地铁场地，地表峰值加速度的变异性随地震动输入峰值加速度的增加而减小，而对于星海广场地铁场地，地表峰值加速度的变异性随地震动输入峰值加速度的增加而略有增加。地震动峰值加速为0.1g时，地表峰值加速度的变异性随场地土变软而减小，地震动峰值加速为0.3g时，峰值加速度的变异性随场地土变软而增大，由此可见输入地震动的强度越大，软弱土层的加速度变异性越大。

图7 场地沿深度方向加速度放大系数

图8 地表加速度峰值均值和变异性

3.2 相对速度峰值

场地各深度方向相对速度取地表速度绝对值最大时刻的值，并绘出图9。从图中可以看出，对于同一类场地，场地相对速度反应随输入加速度峰值的增加而增大。3种场地反应中，输入Ⅱ类场地的地震动比Ⅰ类场的相对速度峰值反应更为明显，这是因为Ⅱ类场地的地震动特征周期比Ⅰ类场地的地震动特征周期偏大，而选取的金安桥站、珠市口站和星海广场站3类场地的卓越周期分别为0.32s、0.49s和1.19s，Ⅱ类场地的地震动特征周期更接近于场地卓越周期。

图9 场地沿深度方向速度峰值分布

图10为地表相对速度峰值的绝对值均值和变异性，从图中可以看出在Ⅱ类场地的地震动输入条件下，相对速度峰值的均值比Ⅰ类场地大，而相对速度峰值的变异性比Ⅰ类场地小。对于金安桥地铁场地和珠市口地铁场地，速度峰值变异性随地震动峰值加速度的增加而略有减小的趋势，对于星海广场地铁场地，速度峰值变异性随地震动峰值加速度的增加而略有增加的趋势。在两类场地的地震动输入条件下，软弱土层场地的地表峰值速度变异性比中、硬土的地表峰值速度变异性大，且较为明显。

图10 地表相对速度峰值均值和变异性

3.3 相对位移峰值

场地深度方向位移峰值取值与相对速度取值的方法相同，即取地表发生最大位移时刻对应的值，如图11。由图可见，对于同一类场地，基岩输入加速度峰值越大，地表位移越大。对比不同的场地条件，星海广场站的软土场地的位移反应峰值最大，金安桥地铁场地的位移反应峰值最小，珠市口地铁场地的位移反应峰值位于两者之间，这与场地条件越软位移反应越大的基本认知是一致的。图12为地表相对位移峰值的均值和变异性，从图中可以看出，在Ⅱ类场地的地震动输入条件下的位移峰值均值比Ⅰ类场地大，而位移峰值的变异性比Ⅰ类场地小。位移峰值变异性随地震动峰值加速度的增大而增大。在两类场地的地震动输入条件下，位移峰值的变异性随场地土变软而增大。

图11 场地沿深度方向相对位移峰值分布

图12 地表相对位移峰值均值和变异性

4 统计分析

图13、14分别为金安桥和星海广场地铁场地在Ⅰ、Ⅱ类场地地震动（0.1g、0.2g、0.3g）输入的情况下地表加速度、速度、位移峰值均值和标准差。结果表明，两组输入地震动条件下的场地反应均值和标准差在输入地震动约50条左右时，其结果基本趋于稳定。

图13 金安桥站加速度、速度、位移峰值均值和标准差（一）

图14 苏州星海广场站加速度、速度、位移峰值均值和标准差（二）

5 结论

本文针对3个典型的实际场地，选取历史上实测到的Ⅰ类和Ⅱ类场地的地震动记录各100条，并分别调整地震动记录加速度峰值至0.1g、0.2g和0.3g，采用一维等效线性化方法开展了场地随机地震反应研究，系统地分析了地震动峰值加速度、峰值速度和峰值位移反应的变异性及其差异。研究表明：

（1）在两类场地的地震动记录输入条件下，峰值加速度反应的均值基本一致；中、硬场地土层地表峰值加速度的变异性随地震动输入峰值加速度的增加而略有减小，软弱场地地表峰值加速度的变异性随地震动输入峰值加速度增加而略有增加。

（2）在Ⅱ类场地的地震动输入条件下，相对速度峰值的均值比Ⅰ类场地大，而相对速度峰值的变异性比Ⅰ类场地小；在两类场地的地震动输入条件下，软弱土层场地的速度峰值变异性比中、硬土层场地大，且较为明显。

（3）在Ⅱ类场地的地震动输入条件下，位移峰值的均值比Ⅰ类场地大，而位移峰值的变异性比Ⅰ类场地小；位移峰值变异性随地震动输入峰值加速度的增大而增大；在两类场地的地震动输入条件下，位移峰值的变异性随场地土变软而增大。

（4）位移和速度峰值变异性较加速度峰值变异性更为突出。

（5）通过统计分析，认为选取50条左右地震动记录即可获得场地地震反应均值和标准差较为稳定的结果。

薄景山，李秀领，李山有，2003．场地条件对地震动影响研究的若干进展．世界地震工程，19（2）：11—15．

陈龙伟，冯浩，袁晓铭，2015．基于单台强震数据软厚场地地震反应不确定性分析．土木建筑与环境工程，37（3）：35—41．

陈永新，迟明杰，李小军，2016．基于强震动记录确定的场地卓越周期．地震学报，38（1）：138—145．

胡聿贤，1988．地震工程学．北京：地震出版社，84—124．

兰景岩，吕悦军，刘红帅，2012．地震动强度及频谱特征对场地地震反应分析结果的影响．震灾防御技术，7（1）：37—45．

李建亮，亢川川，何玉林等，2011．场地条件对地面加速度峰值离散性的影响研究．震灾防御技术，6（4）：416—426．

李天，李杰，1994．具有随机参数的场地地震反应分析．岩土工程学报，16（5）：79—83．

李小军，1989．一维土层地震反应线性化计算程序．见：廖振鹏编，地震小区划（理论与实践）．北京：地震出版社，250—265．

廖振鹏，李小军，1989．地表土层地震反应的等效线性化解法．见：廖振鹏编，地震小区划（理论与实践）．北京：地震出版社，141—153．

刘方成，武景芳，陈斌等，2015．软土层几何特性与剪切波速对场地峰值加速度的影响．震灾防御技术，10（4）：840—852．

吕红山，赵凤新，2007．适用于中国场地分类的地震动反应谱放大系数．地震学报，29（1）：67—76．

欧阳行艳，章文波，2008．利用强震记录进行场地反应分析研究综述．世界地震工程，24（3）：118—126．

齐文浩，2004．土层地震反应分析方法的比较研究．哈尔滨：中国地震局工程力学研究所．

齐文浩，薄景山，张忠利，2010．土层地震反应分析的研究现状．世界地震工程，26（S1）：368—372．

童广才，刘康和，2000．场地卓越周期的确定．电力勘测，26（2）：43—46．

赵松戈，胡聿贤，廖旭，2000．土层参数的随机性对场地传递函数的影响．地震工程与工程振动，20（2）：7—12．

中华人民共和国住房和城乡建设部，中华人民共和国国家质量监督检验检疫总局，2010．GB 50011—2010建筑抗震设计规范．北京：中国建筑工业出版社．

中华人民共和国住房和城乡建设部，2011．JGJ/T 97—2011 工程抗震术语标准．北京：中国建筑工业出版社．

Bardet J. P., Ichii K., Lin C. H., 2000. EERA—a computer program for equivalent-linear earthquake site response analyses of layered soil deposits. Los Angeles: University of Southern California.

Hardin B. O., Drnevich V. P., 1972. Shear modulus and damping in soils: design equations and curves. Soil Mechanics and Foundations Division, 98(SM7): 667—692.

Idriss I. M., Seed H. B., 1968. Seismic response of horizontal soil layers. Journal of the Soil Mechanics and Foundations Division, 94(SM4): 1003—1031.

Idriss I. M., Sun J., 1992. User's manual for SHAKE91—a computer program for conducting equivalent linear seismic response analyses of horizontally layered soil deposits: center for Geotechnical Modeling. Davis, California: Department of Civil & Environmental Engineering, University of California.

Kondner P. L., 1963. Hyperbolic stress-strain response: cohesive soils. Journal of the Soil Mechanics and Foundations Division, 89(SM1): 115–143.

Trifunac M. D., 2016. Site conditions and earthquake ground motion-a review. Soil Dynamics and Earthquake Engineering, 90: 88—100.

杜修力，袁雪纯，黄景琦，许紫刚，2017．典型土层场地随机地震反应规律分析．震灾防御技术，12（3）：574—588．

Analysis of Stochastic Seismic Response in Typical Soil Sites

Du Xiuli, Yuan Xuechun, Huang Jingqi and Xu Zigang

（Beijing University of Technology, Beijing 100124, China）

For the hard, medium and soft soil sites, the ground motion of the class I and class II bedrock site were selected from measured data in history and each group contained 100 ground motions. The peak acceleration of ground motion records were adjusted to 0.1g, 0.2g and 0.3g, then the stochastic seismic response of the site was studied by using the one dimensional equivalent linearization method to analysis systematically the variability law of peak ground acceleration, velocity and displacement response. We concluded out that the peak value of the displacement and the velocity variability will increase when soil becoming softer, and the peak value of the displacement and the velocity variability is more prominent than peak acceleration variability. In order to obtain stable results of the mean and standard deviation of site earthquake response about 50 ground motion records are required.

Typical soil site; Stochastic seismic response; Variability law

10.11899/zzfy20170314

国家自然科学基金重大研究计划（91215301），国家自然科学基金创新研究群体项目（51421005），教育部“创新团队发展计划（IRT13044）

2017-07-10

杜修力，男，生于1962年。教授。主要从事地震工程领域研究。E-mail：duxiuli@bjut.edu.cn