放大效应

平面SH波作用下部分充填圆弧形沉积谷的二维土层和地形放

大效应

张宁;高玉峰;何稼;徐婕;陈欣;代登辉

【摘要】Siteamplificationconsistsoftwo-dimensional(2D)soillayer

ectiveof

thisstudyistoinvestigatetherelativecontributionofsoilamplification

end,thewavefunctionriessolution

forthescatteringofplaneSHwavesbyapartiallyfilledcircular-arcalluvial

efunctionriessolutionisobtainedbyanovel

methodofover-determinedsystemofequationsintheframeworkofthe

wavefunctionexpansiontechniquewiththeaidofaregion-matching

vergencetestsareconductedtorevealthenecessityof

idityof

adjustingthematerialparametersofthetwosubregionsintheanalytical

model,bothofthesurfacemotionsofthealluvialvalleyandtheempty

eamplificationpatternsofthealluvialvalley

arecomparedwiththetopographicamplificationpatternsoftheempty

canyonwiththesamegeometry,theresultsshowthatthesoilamplification

effectsareusuallylargerthanthetopographicamplificationeffectswithin

thealluvialvalley,whilethetopographiceffectsdominatetheamplification

ards,a

parametricstudyintermsofthemaximumsurfacemotioniscarriedoutto

determinetherelativeimportanceofsoilandtopographiccontributionsin

amorecomprehensivemannerandtofurthercharacterizethe2Dsoillayer

identthatthemaximumsoillayeramplification

2D

soilamplificationincreaswiththeimpedancecontrastbetweenthesoil

layerandtheunderlyingbedrockbutisnotasimplelinearsuperposition

of1Dsoilamplificationand2Dtopographicamplification.%场地效应通常包

含土层放大效应和地形放大效应,为了揭示二者的相dream什么意思 对贡献,本文构造了平面SH波

作用下部分充填沉积谷的解析模型,借助于区域分解策略,在波函数展开法的框架下,

提出了超定方程组解法,得到了部分充填圆弧形沉积谷对平面SH波散射的波函数

级数解,而且级数解的收敛测试表明了超定方程组解法的必要性.通过与文献结果进

行对比,验证了本文方法的正确性.通过调整解析模型中两个子区域内的材料参数,计

算了沉积谷引起的场地放大效应和相应的空河谷引起的地形放大效应.对二维土层

与地形效应进行对比分析,结果显示,在沉积谷内二维土层放大效应通常强于地形放

大效应,而地形放大效应决定了沉积谷外的地面运动放大形态.针对最大地面运动,进

行了沉积谷和相应空河谷的参数分析,进一步描述了二维土层放大效应,研究结果表

明二维土层放大效应引起的最大地面运动通常远远大于地形放大效应引起的最大地

面运动,并且二维土层效应通常随着土层与基岩的阻抗比的增大而增大,但不是一维

土层放大效应与二维地形放大效应的简单线性叠加.

【期刊名称】《地震学报》

【年(卷),期】2017(039)005

【总页数】20页(P778-797)

【关键词】二维土层效应;地形效应;平面SH波;散射;沉积谷地;波函数展开

【作者】张宁;高玉峰;何稼;徐婕;陈欣;代登辉

【作者单位】中国南京210098河海大学岩土力学与堤坝工程教育部重点实验室;

中国南京210098河海大学土木与交通学院岩土工程科学研究所;中国南京

210098河海大学岩土力学与堤坝工程教育部重点实验室;中国南京210098河海

大学土木与交通学院岩土工程科学研究所;中国南京210098河海大学岩土力学与

堤坝工程教育部重点实验室;中国南京210098河海大学土木与交通学院岩土工程

科学研究所;中国南京210098河海大学岩土力学与堤坝工程教育部重点实验室;中

国南京210098河海大学土木与交通学院岩土工程科学研究所;中国南京210098

河海大学岩土力学与堤坝工程教育部重点实验室;中国南京210098河海大学土木

与交通学院岩土工程科学研究所;中国南京210098河海大学岩土力学与堤坝工程

教育部重点实验室;中国南京210098河海大学土木与交通学院岩土工程科学研究

所

【正文语种】中文

【中图分类】P319

众所周知，地震中地面运动会由于局部场地效应而发生重大改变(Snchez-Sesma,

1987;Semblat,2011)．地球表面地质条件如土层和地形的变化对地震波有强烈

影响，会导致地面运动的放大．地面运动的场地放大和空间变化对地震减灾和重

要建筑的地震设计特别重要．

地震动场地放大一直是地震学家和岩土工程师研究的焦点．在场地效应研究的早

期阶段，他们的主要兴趣有明显的不同．地震学家关注二维或三维复杂地质结构

的弹性反应以及相应的地震波散射和衍射现象，而岩土工程师主要对强震中土的

非线性效果(Seed，Idriss,1969)以及土层对竖向传播剪切波的一维放大感兴趣．

近年来，场地效应研究在上述两个领域内均取得了很大的进展，并且有同时考虑

地震学因素和岩耳保健操 土因素的趋势．一方面，很多岩土研究(Bielaketal,1999;Bai

etal,2003)显示出一维分析在描绘场地放大形态上的局限性；另一方面，最近的

一些地震学研究也强调了土的非线性效果(Gelagotietal,2010,2012)．事实上，

二维沉积谷效应和三维沉积盆地效应在岩土地震工程和地震学领域取得了越来越多

的研究进展(Semblatetal,2005;杨彩红等，2006；Psarropoulotal,2007；

Chaillatetal,2009；Lanzo，Pagliaroli,2009；巴振宁，梁建文，2014;刘

中宪，黄磊，2015；Zhangetal，2015)．

不规则地质结构的场地放大包含材料特性方面的贡献(土层效应)和几何形状方面的

贡献(地形效应)．然而，上述大部分场地效应研究关注的是沉积谷或沉积盆地的

整体地震反应，而未能将二维土层效应与地形效应分开．地震学研究表明凹陷河

谷地形能够降低地面运动(Gaoetal,2012;Zhangetal,2012a,b)，而岩土地震

工程学的研究表明阻抗比的存在会导致土层地面运动放大．因此土层和地形条件

在改变地面运动方面的相对重要性值得进一步深入研究.

揭示土层和地形的相对角色在场地反应研究中具有重要意义(Zhangetal，2017)，

目前学术界关于此问题的研究逐渐得以开展．例如：受雅典1999年地震中的观

测证据启发，Gazetas等(2002),Assimaki和Gazetas(2004)，Assimaki等

(2005a,b)及Assimaki和Kaul(2007)利用有限元和谱单元方法研究了二维土

层效应和地形效应对于不均匀震害分布的相对贡献；Gatmiri和Arson(2008)及

Gatmiri等(2008)提出了通过有限元和边界元杂交的方法来区分地震场地放大的二

维土层效应与地形效应．

本文旨在系统地探究二维土层效应与地形效应之间的关系，进一步加深对场地效

应的理解，提出部分充填圆弧形沉积谷受平面SH波作用的解析模型，在波函数

展开法的基础上采用超定方程组策略，推导这一散射问题的解析级数解，并进行

验证．利用本文提出的级数解，计算部分充填圆弧形沉积谷引起的场地放大效应

以及相应空河谷引起的地形放大效应，分析二维土层及地形各自引起的地面运动

放大形态和最大放大值．本文将部分充填沉积谷和相应空河谷对地震波的散射问

题结合起来对待，旨在增加关于场地效应的一些新认知．

本文所用的二维模型如图1所示．此图表示半空间中的一个圆弧形(即浅圆形)沉

积谷．沉积谷半宽为b，是半径为a的圆柱体的一部分；谷内只是部分充填了均

匀的土沉积物，其上表面到水平地表的竖向距离为d，最大深度为d1．谷内土

沉积物(子区域1)与半空间中基岩(子区域2)之间的交界面假设为完美焊接．模型

中介质的所有材料特性均假设为弹性、均质和各项同性．介质的剪切模量j，

剪切波速cj和密度j均为常数．下标j=1,2分别表示子区域1和子区域2．如

果子区域1和2具有相同的材料特性，则模型为去底浅圆弧形河谷．

模型的激励是一列单位幅值平面SH波，入射角为，圆频率为，位移在z轴

方向．定义3个笛卡尔坐标系和3个极坐标系(图1)：笛卡尔坐标系(x,y)和极坐

标系(r,)的原点设置在沉积谷曲率的中心，距离h表示沉积谷曲率中心相对于水

平地面的高度，沉积谷中心角和土沉积物的中心角分别为和；坐标系(x1,y1)

和(r1,1)的原点均设置在土沉积物的上表面中点；坐标系(x2,y2)和(r2,2)的原

点设置在水平地表和土沉积物中心的正上方．定义向右方向为水平x轴，x1轴和

x2轴的正方向；从竖向y轴逆时针转向x轴方向为角的正方向，角度1和

2亦均以逆时针转向为正．

子区域1和2的位移uj必须满足极坐标形式的二维波动方程(Tsaur，Chang,

2009)

利用

将波动方程转化为亥姆霍兹方程(Helmholtzequation)，得到

由于是稳态问题，时间因子exp(-it)可以省略．

除了式(3)所示的运动方程外，位移uj还应该满足土沉积物的上表面、沉积谷的

弯曲表面以及水平地表的应力自由边界条件，具体表示为

式中b1为土沉积物的上表面半宽．

土沉积物与基岩完美连接的假设要求两个子区域的位移场和应力场均连续，即

.

3.1波场的波函数级数展开与坐标统一

在子区域2中，波场包含两个部分．一个部分是已知的半空间自由场，可表示为

式中k2=/c2为子区域2的剪切波数．

文后将用到极坐标系(r,)中的自由场uf，因此需要进行坐标转换．首先有

式中c2x和c2y分别为基岩波速c2在x和y方向的分量;之后利用x2=x，

y2=y-h，可以获得uf(x,y)，即

最后，将c2x=c2/sin,c2y=c2/cos,x=rsin,y=rcos以及=k2c2代入式

(11)，便得到极坐标系(r,)中的自由场，即

exp(-ik2hcos)exp[ik2rcos(-)].

使用级数展开公式(Abramowitz，Stegun,1964)

式中,Jn()为n阶第一类贝塞尔函数，n为纽曼因子(0=1;n=2,n≥1)，将式

(12)展开为

式中pn和qn的表达式为

qn=[-(-i)nexp(ik2hcos)+inexp(-ik2hcos)]nsin(n).

子区域2中波场的另一部分是由于沉积谷的出现而产生的散射波，可表示为

(Yuan，Liao,1995)

式中,An和Bn为未知的复系数，为n阶第一类汉克尔函数．

要在同一坐标系(r,)中解决问题，需要基于Graf加法公式(Abramowitz，

Stegun,1964)推导适合本文模型的外域型变换公式，即

根据式(18),散射波场us在极坐标系(r不可开交 ,)中可表示为

式中，

Vmn=[Jm+2n+1(k2h)+Jm-2n-1(k2h)].

上述子区域2中的所有波场均满足水平地表的应力自由边界条件,即式(9)，(17)

与式(6)兼容．

构造子区域1中的驻波场uc(r1,1)，以满足式(3)和(4)所给出的亥姆霍兹方程和

土沉积物上表面的应力自由边界条件．驻波场uc(r1,1)可表示为

式中,k1=/c1为子区域1的波数，Cn和Dn为未知复系数．

为了在同一个坐标系中解决问题，从而获得4组未知系数An,Bn,Cn和Dn，本

文推导了Graf加法公式的另一种形式，即

.

根据式(23),波场uc在极坐标系(r,)中可表示为

式中

Nmn=[Jm+2n+1(k1h+k1d)+Jm-2n-1(k1h+k1d)],

因此,同一坐标系下两个子区域中的波场可表示为

u2(r,)=uf(r,)+us(r,).

3.2求解未知系数的超定方程组法

精确严格地满足控制方程(式(3))和土沉积物上表面应力自由条件(式(4))的子区域1

的驻波场以及精确严格地满足控制方程(式(3))和水平地表应力自由边界条件(式(6))

的子区域2中的自由场和散射场均已构造完毕．然而，散射场(式(17))和驻波场

(式(22))均含有未知系数，因此需要利用两个子区域之间的位移和应力场的连续条

件(式(7),(8))以及沉积谷弯曲表面的应力自由条件(式(5))来获得未知系数，最终给

出本文所研究的定解问题的解．

首先，对应力连续条件(式(8))和应力自由条件(式(5))应用正弦和余弦函数的正交

性，对其在[-,]范围内积分，得出未知系数An,Bn,Cn和Dn之间的关系为

式中和的表达式为

式中，=或，g=m或n.

然后，对位移连续条件(式(7))应用正弦和余弦函数的正交性，对其在[-,]范围

内积分，可以获得关于未知系数的另一组关系表达式，即

w=0,1,2,…,

Bn[(k2a)]-Dn[NmnJm(k1a)]=-qnJn(k2a),

w=1,2,3,…,

式中和的表达式可由式(31)和式(32)求出．实际上，上面的4组关系式即式(29)，

(30)，(33)和(34)的推导所用到的技巧与Tsaur和Chang(2008,2009),Zhang等

(2012a)及Gao等(2012)在文章中所使用的技巧相同，均基于加权残量法原

理．本研究在方法上的改进可以在下一步推导中看到．

将式(29)，(30)，(33)和(34)中的无限级数截成有限项来进行数值计算(方法是：

加和指数n,w和m分别截断成N,W和M项，即n=0-N-1,w=0-W-1，

m=0-M-1)．尽管根据加权残量法的原理(Harrington,1967;Lee，Wu,1994)，

只要保证截断后的方程数不小于未知量数即可(W≥N)，但是上述文献均限定

W=N．由于截断矩阵病态，导致带有这种限制的方法对于本文的模型是无效的，

得到的解不稳定(图2)．为了减轻病态矩阵带来的影响，需要充分利用加权残量

法的优势，令W>N．这样，方程组就变成超定的，可以容易地求出稳定且收

敛的最小二乘解．利用超定方程组法，求解由式(29)，(33)组成的方程组可以获

得未知复系数An和Cn，求解由式(30)，(34)组成的方程组可以获得复系数Bn

和Dn，最终得到每个子区域级数形式的波场．需要注意的是，N值和W值需

要根据对应于不同的沉积谷形状和材料特性以及入射波的频率和方向的收敛测试进

行确定，M值满足式(18)，(23)中的Graf加法公式即可．一些收敛测试的例子

将在后文给出．

半空间中任意位置的位移uj(式(27)或式(28))均可以利用前文推导的公式进行计算，

其为包含运动的相位和幅值信息的复函数．位移uj的模量称为位移幅值(Trifunac,

1971)，可表示为

式中Re()和Im()分别为uj的实部和虚部，下标j=1,2分别表示子区域1，

2．本文中的位移幅值|u|均利用入射波的单位幅值进行了无量纲化．

为了方便，还需定义一个无量纲频率(Trifunac,1971;Chaillatetal,2009)

式中为SH波的入射波长．

3.3收敛测试

为了说明超定方程组法的必要性,图2给出了水平入射(=90)地震波作用下、参

数为h/b=0.5,d/b=0.3,2/1=1.5,c2/c1=2的沉积谷上5个不同位置的位移

幅值随着W/N的变化情况，所选沉积谷的轮廓也在图中给出．为了能够控制级

数解的整体收敛情况，一方面本文选取具有代表性的位置点1—5，其分别位于

水平地面(点1对应x/b=-2，点5对应x/b=2)、河谷侧壁(点2对应x/b=-(b+

b1)/2b,点4对应x/b=(b+b1)/2b)和沉积谷表面(点3对应x/b=0)上；另一方

面本文选择水平入射的SH波(计算表明水平入射比竖向入射和斜入射更难收敛)，

给出了在两种不同频率情况下的结果．由图2可见，对于固定的N值，当

W/N=1时得到的解明显未收敛，只有当W/N≥3时解才开始收敛，并且W/N

的值越大级数解越接近于真解．这说明利用类似于Tsaur和Chang(2008,2009),

Zhang等(2012a)及Gao等(2012)文章中的方法是行不通的，必须用超定方程组

法才能得到本文模型的收敛解．为了给出尽可能高精度的结果，一般取

W/N=20．

图2中，选定的级数解的截断值N=20和N=25均成功地给出了收敛的结果，

实际上，N值需要根据收敛测试确定．图3给出了对应于图2中情形的关于N

值的级数解收敛情况．图3选取W/N=20的超定方程组法给出了两个无量纲频

率的结果，对于较低频率=1，N≥10便几乎可以保证结果的收敛，而对于较

高频率=3，N≥25才可以保证结果的收敛,这说明图2中选择的截断值是合理

的．该收敛测试还表明，频率越大则散射场us(式(17))和驻波场uc(式(22))就需

要截取越多的级数项．一般而言，N≤33,M≤700则足够获得本文所考虑的各种

情况(0<≤4)的收敛结果．

本文提出的方法不需要介质或边界的离散化，也不涉及格林函数；而且，无穷远

处的辐射条件、土沉积物的上表面和半空间中基岩上表面的应力自由边界条件均

可以精确满足．土沉积物与基岩的交界面的连续条件，以及沉积谷弯曲表面的应

力自由条件均可利用超定方程组法实现．除了作类似图2和图3给出的收敛测试，

还可以通过检查交界面上任意点的位移和应力的连续性以及沉积谷弯曲表面的应力

自由条件来保证级数解的收敛．这样，本文提出的超定方程组方法便可以得到

SH波散射问题的高精度解．

本文模型可以看作是半圆形(Trifunac,1971,1972)、浅圆弧形(Cao，Lee,1989;

Todorovska，Lee,1991;Yuan，Liao,1994,1995)、部分充填半圆形(Tsaur，

Chang,2008,2009)的沉积谷和空河谷模型的通解，即上述模型均属于部分充填

圆弧形沉积谷的特例．实际上，如果令h=0,则本文模型退化为部分充填半圆形

沉积谷(1≠2,c1≠c2)和去底半圆形河谷(1=2,c1=c2)；如果令d=0,则本文

模型退化为浅圆弧形沉积谷(1≠2,c1≠c2)和浅圆弧形河谷(1=0,c1=0)；如果

令h=d=0,则本文模型退化为半圆形沉积谷(1≠2,c1≠c2)和半圆形河谷(1=0,

c1=0)．本节选择半圆形沉积谷(Trifunac,1971)和部分充填半圆形沉积谷(Tsaur，

Chang,2008)，用以验证本文级数解的数值精度．

4.1与半圆形沉积谷和自由场的比较

本文模型可以通过退化至半圆形沉积谷(Trifunac,1971)和自由场这两种情况的精

确解而得到理论上的验证．若h和d均等于0，则有

将式(37)-(41)代入式(29),(30),(33)和(34)可以得到

′(k2a)-Dn1(k1a)=-4i(-1)nsin[(2n+1)]2(k2a),

An(k2a)-CnJ2n(k1a)=-2n(-1)ncos(2n)J2n(k2a),

Bn(k2a)-DnJ2n+1(k1a)=-4i(-1)nsin[(2n+1)]J2n+1(k2a)，

式(42)—(45)重现了Trifunac(1971)的结果(对应其论文中式(14)和(15))，只是使

用了不同的表示符号．而且，若子区域1和子区域2中的介质相同(即1=2，

k1=k2)，则式(42)—(45)简化为

′(k2a)-Dn(k2a)=-4i(-1)nsin[(2n+1)](k2a),

An(k2a)-CnJ2n(k2a)=-2n(-1)ncos(2n)J2n(k2a),

Bn(k2a)-DnJ2n+1(k2a)=-4i(-1)nsin[(2n+1)]J2n+1(k2a).

根据式(46)—(49)可以得到An=Bn=0,Cn=p2n,Dn=q2n+1，则半空间只存在

自由场，即u1(r,)=u2(r,)=uf(r,)．

4.2与部分充填半圆形沉积谷的比较

如果只有h=0，则部分充填浅圆弧形沉积谷退化为部分充填半圆形沉积

谷．Tsaur和Chang(2008)已经推导了部分充填半圆形沉积谷对SH波的散射的

波函数级数解，因此可以通过与其对应的某些特定几何与材料的计算结果进行比

较来验证本文方法的数值精度．

为简便起见，图4仅给出了对应于h/b=0,d/b=0.5,2/1=1.5,c2/c1=2，

=1的结果．图4a--d分别对应于地震波竖向入射(=0)、斜入射(=30,60)

和水平入射(=90)等情况．场地表面的位移幅值或放大形态沿着水平x轴画

出．所画图形的无量纲距离x/b在-4—4之间．不规则场地范围为-

1≤x/b≤1．图4所示结果表明本文方法的结果(实线)与Tsaur和Chang(2008)

在其图4中对应于h/a=0.5(即本文中d/b=0.5)时的结果(虚线)非常吻合．对于

其它情况本文方法所得结果与文献中的结果也是高度吻合的．

5.1二维土层与地形效应的放大形态对比

一般而言，场地放大效应可以通过比较包含不规则体的半空间中的位移幅值与没

有不规则体的半空间中的位移(自由场)来理解．因为对于单位幅值的平面波入射，

地面不同观察点处的自由场位移幅值|uf|为常数2(地面以下观察点处的|uf|不是常

数)，所以本文方法计算得到的地表位移幅值|u|可以看作放大因子的两倍．当某

一点的位移幅值|u|大于2，则该点的地面运动是放大的;反之，当位移幅值小于2，

则该点的地面运动出现衰减．

袁晓铭等(2002)通过对比圆弧状沉积盆地与软土水平覆盖层的出平面地表运动，

系统地研究了一维土层效应与二维沉积谷效应，然而，直接理解二维土层效应是

有难度的，例如Gatmiri和Arson(2008)认为“一维沉积土层引起的放大现象是

已知的，但在二维地形中分析岩土对场地效应的贡献更为复杂”．本文推导的波

函数级数解可以用来分别计算部分充填浅圆弧形沉积谷和同样形状的去底浅圆弧形

河谷的放大形态．模型中沉积谷的场地放大效应包含土沉积物的贡献和河谷地形

的贡献，因此，本文可以通过将河谷的地形效应从沉积谷的场地效应中进行分离

来获得二维土层效应．对于河谷，材料参数比2/1和c2/c1均设置为1．仿

照Trifunac(1971)，Todorovska和Lee(1991)，Yuan和Liao(1995)及

Tsaur和Chang(2008)的做法，将沉积谷密度比2/1固定为1.5，剪切波速

比c2/c1设为2或3．

对应h/b=0.5，d/b=0.4和4个入射角度=0,30,60和90的沉积谷和去底

河谷的放大形态在图5中给出．这些结果对应于无量纲频率=1，该频率意味着

入射波的波长()等于不规则场地的宽度(2b)．无量纲参数h/b和d/b代表场地

(沉积谷或空河谷)的形状．根据最大位移幅值判断，在各种入射角度的地震波作

用下，在不规则场地范围内，沉积谷引起的场地放大效应大于空河谷引起的地形

放大效应；而在不规则场地范围外，水平地面上对应于两种情况的放大形态区别

不大．因此，可认为二维土层的放大效应在沉积谷内占据绝对优势，地形放大效

应在沉积谷外占据绝对优势．特别地，二维土层放大效应可以改变凹陷地形的放

大形态．例如，图5c中，在x/b=0处的运动受地形效应影响时衰减(实线)，而

受二维土层和地形的总效应影响时却是放大的(虚线)．

一些关于地形效应的参数研究认为河谷侧壁的倾斜度对放大形态影响很大(Gaoet

al，2012；Zhangetal，2012a)．为了揭示不规则场地的形状对土效应和地形

效应的影响，图6

给出了当h/b=0.3,d/b=0.4时沉积谷和去底河谷在=1和4种入射角度下的放

大形态．可以看出，图6中的沉积谷放大效应要大于图5，而图6中的地形放大

效应与图5相差不大．图6与图5所使用的参数中唯一不同的是h/b．随着形状

参数h/b的降低，沉积谷或河谷的侧壁变得更陡，土沉积物的体积变大．考虑

到由于形状参数h/b略微减小时地形放大效应增加很少，推断不同入射角度下的

场地放大效应的增加主要是由于土沉积物体积的变化所致．另外，图6中也可以

明显看出，土的放大效应在场地范围内占据优势，而地形效应在场地范围外起主

要作用．

除了几何形状，影响场地效应的重要因素还有土沉积物的材料参数以及入射波的

频率和角度．因此，图7给出了更高频率=2的入射波的结果.对比图6可以看

出，频率的增加导致地表位移幅值随位置x/b的变化而波动得更频繁、更剧

烈．图7实线表示去底河谷(2/1=1,c2/c1=1)的位移幅值，虚线和点线分别

对应充填硬土(2/1=1.5,c2/c1=2)和充填软土(2/1=1.5,c2/c1=3)的沉积

谷．通过比较硬土沉积谷与软土沉积谷的放大形态，认为二维土层效应(沉积谷场

地效应减去河谷地形效应)也会受到土沉积物与半空间中的基岩的材料参数比

(c2/c1)的影响．当地震波竖向入射(=0)时，充填硬土的沉积谷引起的地面运动

大于充填软土的沉积谷;当地震波斜入射(=30,60)和水平入射(=90)时，充

填软土的沉积谷引起的地面运动放大更显著．另外，图7中的结果再一次证实了

图5和图6中基于=1的结论，即二维土层放大效应会增加凹陷场地内部的位移

幅值,而地形放大效应决定了不规则场地外的反应形态,而且这个现象随着频率的

增长变得更明显．

至于入射角度对场地放大形态的影响，观察图5—7中的结果可以看出，对于竖

向入射情况，由于河谷和沉积谷几何形状的对称性，地表位移幅值也是对称

的．河谷地形效应引起的地表相对最大位移幅值发生在河谷两肩(=2)或河谷底部

(=1)，沉积谷场地效应引起的地表相对最大位移幅值总是发生在土沉积物的范围

内．随着观察点逐渐远离河谷油泼辣椒 或沉积谷，水平地表的位移交替放大或衰减，这是

因为在不同位置发生了波的相长干涉或相消干涉．对于水平入射情况，右手侧水

平地表上的位移通常衰减，这是由于河谷或沉积谷的阴影效果．河谷地形效应引

起的地表相对最大位移幅值发生在河谷左肩，但由于强烈的二维土层放大效应，

沉积谷场地效应引起的地表相对最大位移幅值仍然发生在土沉积物的范围内．对

于斜入射情况，表面运动的形态介于竖向与水平向情况之间．

为了更全面地了解入射波的频率对场地效应的影响，本文选择一个固定形状的不

规则场地，给出了与其对应的作为无量纲距离x/b和无量纲频率的函数的表面

位移幅值结果，如图8所示．可以看出，随着频率的增加，两种场地位移幅值

均随着位置x/b的变化而更频繁地波动．另外，从图8中可以明显地看出，在

不规则场地范围内，各种频率下的激励所引起的二维土层效应(即图8b中的场地

效应减去图8a中的地形效应)远远大于地形效应．值得注意的是,本文选择了地形

效应比较明显的场地，因为河谷的深度(深宽比d/b为0.4)要大于沉积谷内充填的

土沉积物的厚度(厚宽比d1/b约为0.218)．可见由二维土层和地形引起的放大差

异是不容忽视的，因其会引起沉积谷内外较大的地面运动差异，尤其是横跨在沉

积有软弱材料的河谷上的大型建筑如水坝、桥梁等，在抗震设计时需要特别注意

此放大差异．

5.2二维土层与地形效应的最大地表位移幅值对比

上一节通过地表运动形态揭示了土效应和地形效应的相对角色．图5—8说明二

维土层效应和地形效应受不规则场地的几何形状和材料特点、入射波的频率和方

向所影响．众所周知，最大地表运动在地震工程中具有重要意义(Semblatetal,

2010),本节将根据最大地表位移幅值针对这些影响因素作系统的参数分析，以期

进一步明确这些参数是如何影响二维土层效应和地形效应的．

首先，选择一个特定形状，计算出在沉积谷及其特例河谷附近水平地面上3个入

射角度=0,45,90和两种波速比下作为频率的函数的最大位移幅值，结果如

图9所示．该图形状很有代表性，因为土的充填比d1/(d1+d)大约为0.5．可

以看出在频率小于0.4的区间内，场地效应与地形效应引起的最大地表位移幅值

相差不大．在0.4—4.0的频率区间内，从最大幅值来看，场地效应远远超过地

形效应．虽然相对于较硬的沉积谷，较软的沉积谷的最大位移幅值随频率变化有

更大的波动，但是二维土层效应一般会随着材料波速比的增大而变大，即沉积土

越软，最大放大值往往越大．值得说明的是，地形效应的最大位移幅值随着频率

增加而渐渐增加，但始终低于5.0，而场地效应的最大位移幅值却可以高达30，

再次证明了二维土层放大效应的主导地位．

其次，为了更全面地揭示形状参数对地形效应和二维土层效应的作用，本文选取

了9种形状的河谷或沉积谷计算出这些不规则场地附近的最大地表位移幅值．图

10为9种河谷或沉积谷的形状示意图．表1，2给出了对应两个无量纲频率=1,

4下的两种沉积谷(c2/c1=2.0,3.0)及其特例河谷(c2/c1=1.0)附近的最大地表位移

幅值，入射方向=0，45，90．从表1，2中的结果可以得出以下结论：

1)地形效应(2/1=1.0,c2/c1=1.0)引起的最大地表位移幅值明显小于场地效应

(2/1=1.5,c2/c1=2.0,3.0)，并且场地效应引起的最大地表位移幅值一般会随

着波速比c2/c1的增大而增大，此情况约占80%，这与基于图9中特定形状的

沉积谷所得到的规律是一致的.鉴于本文选取的几何形状的广泛性和代表性(图10)，

可以预测这些规律对于其它类似形状的场地情况也是适用的．

2)地形效应和场地效应引起的最大地表位移幅值基本上均随着无量纲频率从1

增加到4而有所增大，此情况约占80%．

3)地形效应引起的最大地表位移幅值均随着入射角度的增大而增大，即水平入

射时最大地表位移幅值最大，竖向入射时最大地表位移幅值最小．而场地效应引

起的最大地表位移幅值往往不符合此规律．

4)由于高频波(=4)入射时产生最大地表位移幅值的位置变化很大，而对于低频

结果(=1)，出现最大地表位移幅值的位置比较固定(地震波竖向入射时通常出现

在河谷底部中点，斜入射和水平入射时通常出现在河谷左肩附近)，对低频结果进

行分析后认为，河谷地形效应(2/1=1.0,c2/c1=1.0)引起的最大位移幅值会随

着h/b的减小和d/b的增大而增大，即河谷侧壁越陡、深度越深，地形效应越

明显．然而，沉积谷场地效应引起的最大位移幅值随着形状参数h/b和d/b的变

化无明显规律，说明形状参数对二维土层效应的影响非常复杂，二维土层效应不

是一维土层效应和地形效应的简单叠加，这要归因于地震波的二维散射和衍射．

本文利用波函数展开并结合超定方程组策略得到了部分充填圆弧形沉积谷及相应的

空河谷对平面SH波散射的级数解，这两种场地的波函数级数解均是对已有解析

解的补充．新的波函数级数解的提出通常需要用到一系列解析技术(Luo,2008)，

例如大圆弧近似(Todorovska，Lee,1991)、镜像方法(Leeetal,1999;Gao，

Zhang,2013)、辅助函数法(Yuan，Men,1992)、半余弦展开(Leeetal,2006)

等，本文提出的超定方程组策略是对波函数展开方法的新发展，为分离变量法进

一步应用到其它不规则场地的情况开启了新的可能．

利用级数解计算分析了二维土层与地形放大效应．研究表明：二维土层放大效应

在沉积谷范围内占据支配地位，地形放大效应在沉积谷范围外微信风景头像 占据支配地位；不

管地震波的入射方向如何变化，场地效应引起的地面运动的相对最大值总是倾向

于发生在土沉积物的范围内，且二维土层放大效应通常大于地形放大效应．此结

果对于大跨度结构如水坝、桥梁的地震设计非常重要，因为这些结构可能会经历

由二维土层和地形的差异放大效应而引起的很强的空间变化地面运动．在大部分

抗震规范中地表地形的放大效应均被忽视了，建议应该同时考虑二维土层效应和

地形效应．

关于沉积谷场地效应和河谷地形效应分别引起的最大地面运动幅值的参数分析揭示

了材料参数c2/c1、形状参数h/b和d/b、无量纲频率和波入射角度对二维

土层效应和地形效应的强烈影响．一般而言，二维土层效应会随着土沉积物的材

料软化而增强，但其随形状的变化规律非常复杂，不是简单的一维土层效应与二

维地形效应的叠加．在一件幸福的事 较小的频率范围内，二维土层效应非常小，地形效应占主

导地位，而超过某一无量纲频率(如=0.4)后，二维土层效应引起的最大放大值

远远大于地形效应．

本文基于弹性动力学基本理论，初步探索了平面SH波作用下的二维土层和地形

的弹性放大效应，以此为参照，下一步研究工作需要考虑其它类型的地震波以及

土层的复杂弹塑性本构关系对沉积谷放大效应的影响.

巴振宁,梁建文.2014.瑞雷波斜入射下层状半空间中沉积谷地周围的三维散射研究

[J].地震学报,36(4):571--583.

BaZN,LiangJW.2014.3Dscatteringbyanalluvialvalleyembeddedina

layeredhalf莜麦面的功效与作用 -spaceforobliquelyincidentRayleighwaves[J].Acta

SeismologicaSinica,36(4):571--583(inChine).

刘中宪,黄磊.2015.含软夹层层状沉积谷地对P、SV和Rayleigh波的散射[J].地

震工程与工程振动,35(1):105--117.

LiuZX,HuangL.2015.P,SVandRayleighwavesscatteringinlayered

alluvialvalleyswithsoftinterlayer[J].EarthquakeEngineeringand

EngineeringVibration,35(1):105--117(inChine).

杨彩红,梁建文,张郁山.2006.多层沉积凹陷地形对平面SH波散射问题的解析解

[J].岩土力学,27(12):2191--2196.

YangCH,LiangJW,icalsolutionofscatteringof

planeSHwavesaroundacanyonwithmulti-layeredalluvia[J].Rockand

SoilMechanics,27(12):2191--2196(inChine).

袁晓铭,李雨润,孙锐.2002.圆弧状沉积盆地与软土单覆盖层出平面地表运动对比

[J].地震工程与工程振动,22(4):16--21.

YuanXM,LiYR,isonofout-of-planesurfaceground

motionbetweenacircular-arcalluvialvalleyandasingleoverburdensoft

layer[J].EarthquakeEngineeringandEngineeringVibration,22(4):16--21

(inChine).

AbramowitzM,okofMathematicalFunctions:

WithFormulas,Graphs,andMathematicalTables[M].NewYork:Dover:363.

AssimakiD,dtopographicamplificationoncanyon

banksandthe1999Athenarthquake[J].JEarthqEng,8(1):1--43.

AssimakiD,GazetasG,soflocalsoilconditionson

thetopographicaggravationofismicmotion:Parametricinvestigation

andrecordedfieldevidencefromthe1999Athenarthquake[J].Bull

SeismolSocAm,95(3):1059--1089.

AssimakiD,KaulE,-dependenttopographic

effects:Acastudyfromthe1999Athenarthquake[J].Earthquake

Spectra,21(4):929--966.

AssimakiD,edtopographicamplificationfactorsfor

asingle-facedslopeduetokinematicsoil-structureinteraction[J].J

GeotechGeoenvironEng,133(11):1414--1431.

BaiLG,GlarSD,sponatTreasureandYerba

BuenaIslands,California[J].JGeotechGeoenvironEng,129(6):415--426.

BielakJ,XuJF,uakegroundmotionandstructural

responinalluvialvalleys[J].JGeotechGeoenvironEng,125(5):413--423.

CaoH,ringofplaneSHwavesbycircularcylindrical

canyonswithvariabledepth-to-widthratio[J].EurEarthqEng,3(2):29--37.

ChaillatS,BonnetM,stmulti-domainBEMto

modelismicwav婴儿歌曲大全 epropagationandamplificationin3-Dgeological

structures[J].GeophysJInt,177(2):509--531.

GaoY,ZhangN,LiD,LiuH,CaiY,softopographic

amplificationinducedbyaU-shapedcanyononismicwaves[J].Bull

SeismolSocAm,102(4):1748--1763.

GaoY,ringofcylindricalSHwavesinducedbya

symmetricalV-shapedcanyon:Near-sourcetopographiceffects[J].

GeophysJInt,193(2):874--885.

GatmiriB,csiteeffectsbyanoptimized2DBE/FE

ficationofsiteeffectsintwo-dimensionaldimentary

valleys[J].SoilDynEarthquakeEng,28(8):646--661.

GatmiriB,ArsonC,csiteeffectsbyanoptimized

2DBE/,numericaloptimizationandapplicationto

topographicalirregularities[J].SoilDynEarthquakeEng,28(8):632--645.

GazetasG,KallouPV,aphyandsoileffects

intheMS5.9Parnitha(Athens)earthquake:ThecaofAdmes[J].Natural

Hazards,27(1/2):133--169.

GelagotiF,KourkoulisR,AnastasopoulosI,TazohT,GazetasG.2010.

Seismicwavepropagationinaverysoftalluvialvalley:Sensitivityto

ground-motiondetailsandsoilnonlinearity,andgenerationofaparasitic

verticalcomponent[J].BullSeismolSocAm,100(6):3035--3054.

GelagotiF,KourkoulisR,AnastasopoulosI,ear

dimensionalanalysisoftrapezoidalvalleyssubjectedtovertically

propagatingSVwaves[J].BullSeismolSocAm,102(3):999--1017.

methodsforfieldproblems[J].ProcIEEE,55(2):

136--149.

LanzoG,calmodelingofsiteeffectsatSan

GiulianodiPuglia(SouthernItaly)duringthe2002Moliismic

quence[J].JGeotechGeoenvironEng,135(9):1295--1313.

LeeVW,ationoftheweightedresidualmethodto

diffractionby2-Dcanyonsofarbitraryshape:ntSHwaves[J].Soil

DynEarthqEng,13(5):355--364.

LeeVW,ChenS,anediffractionfromcanyonabove

subsurfaceunlinedtunnel[J].JEngMech,125(6):668--675.

LeeVW,LuoH,ane(SH)wavesdiffractionbya

micircularcylindricalhillrevisited:Animprovedanalyticwaveries

solution[J].JEngMech,132(10):1106--1114.

ctionofSH-WavesBySurfaceorSub-Surface

TopographiesWithApplicationtoSoil-StructureInteractiononShallow

Foundations[D].LosAngeles:UniversityofSouthernCalifornia:2--5.

PsarropoulosPN,TazohT,GazetasG,and

nonlinearvalleyamplificationeffectsonismicgroundmotion[J].Soils

Found,47(5):857--871.

【相关文献】

fectsonstronggroundmotion[J].SoilDynEarthquakeEng,

6(2):124--132.

SeedHB,nceofsoilconditionsongroundmotionsduring

earthquakes[J].JSoilMechFoundDiv,95(1):99--138.

SemblatJF,KhamM,PararaE,BardPY,PitilakisK,MakraK,c

waveamplification:Basingeometryvssoillayering[J].SoilDynEarthquakeEng,25(7/10):

529--538.

SemblatJF,LokmaneN,Driad-LebeauL,mplificationofdeep

mininginducedvibrationspart.2:Simulationofgroundmotioninacoalbasin[J].SoilDyn

EarthquakeEng,30(10):947--957.

ngismicwavepropagationandamplificationin1D/2D/3D

linearandnonlinearunboundedmedia[J].IntJGeomech,11(6):440--448.

TodorovskaMI,emotionofshallowcircularalluvialvalleysfor

incidentplaneSHwaves-analyticalsolution[J].SoilDynEarthquakeEng,10(4):192--200.

emotionofami-cylindricalalluvialvalleyforincidentplane

SHwaves[J].BullSeismolSocAm,61(6):1755--1770.

ringofplaneSHwavesbyami-cylindricalcanyon[J].

EarthquakeEngStructDyn,1(3):267--281.

TsaurDH,-wavesscatteringfromapartiallyfilledmi-circular

alluvialvalley[J].GeophysJInt,173(1):157--167.

TsaurDH,ringofSHwavesbytruncatedmicircularcanyon[J].J

EngMech,135(8):862--870.

YuanXM,ringofplaneSHwavesbyami-cylindricalhill[J].

EarthquakeEngStructDyn,21(12):1091--1098.

YuanXM,ringofplaneSHwavesbyacylindricalcanyonofcircular-

arccross-ction[J].SoilDynEarthquakeEng,13(6):407--412.

YuanXM,ringofplaneSHwavesbyacylindricalal理想作文600字 luvialvalleyof

circular-arccross-ction[J].EarthquakeEngStructDyn,24(10):1303--1313.

ZhangN,GaoYF,CaiYQ,LiDY,ringofSHwavesinducedbya

non-symmetricalV-shapedcanyon[J].GeophysJInt,191(1):243--256.

ZhangN,GaoYF,LiDY,WuYX,ringofSHwavesinducedbya

symmetricalV-shapedcanyon:Aunifiedanalyticalsolution[J].EarthquakeEngEngVib,

11(4):445--460.

ZhangN,GaoYF,YangJ,yticalsolutiontothescatteringof

cylindricalSHwavesbyapartiallyfilledmi-circularalluvialvalley:Near-sourcesite

effects[J].EarthquakeEngEngVib,14(2):189--201.

ZhangN,GaoYF,dtopographiceffectsongroundmotionofa

surficiallyinhomogeneousmi-cylindricalcanyonunderobliqueincidentSHwaves[J].

SoilDynEarthquakeEng,95(1):17--28.