地球物理学报  2021, Vol. 64 Issue (3): 896-906   PDF    
含峭壁V形峡谷对地震SH波散射的解析解
张宁1,2, 潘家琐1,2, 代登辉1,2, 高玉峰1,2     
1. 河海大学岩土力学与堤坝工程教育部重点实验室, 南京 210098;
2. 河海大学土木与交通学院, 南京 210098
摘要:地表地形常引起地震动的局部放大,这是由于地震波传播至局部地形时产生了散射现象.本文利用波函数展开方法和区域匹配技术,提出了含峭壁V形峡谷对平面SH波散射问题的解析解,并进行了退化验证.通过频域内的参数分析,揭示了峭壁深度、入射波频率和角度等因素对峡谷场地地面运动的影响规律,发现上部峭壁会增强峡谷对地震动的地形放大效应.研究结果不仅为数值方法提供了验证基准,还可为含峭壁峡谷周边建筑物的抗震设计提供顺河向地震动输入.
关键词: 地形效应      平面SH波      散射      波函数展开      解析解     
An analytical solution to the scattering of plane SH waves by a V-shaped canyon with cliffs
ZHANG Ning1,2, PAN JiaSuo1,2, DAI DengHui1,2, GAO YuFeng1,2     
1. Key Laboratory of Ministry of Education for Geomechanics and Embankment Engineering, Hohai University, Nanjing 210098, China;
2. College of Civil and Transportation Engineering, Hohai University, Nanjing 210098, China
Abstract: Local amplification of ground motions often occurs due to the scattering of seismic waves when propagating around surface topography. This paper, using the method of wave function expansion and the technique of region matching, proposes an analytical solution to the scattering of plane SH waves by a V-shaped canyon with cliffs and validates it by the degradation. Based on a parametric analysis in the frequency domain, the effects of cliff depth, incident wave frequency and angle on the ground motion of the canyon site are revealed. It is found that the upper cliff can enhances the amplification effect of the canyon on ground motions. The research results provide not only a reference for numerical methods, but also an earthquake input in the upstream-downstream direction for the seismic design of buildings around canyons with cliffs.
Keywords: Topographic effect    Plane SH wave    Scattering    Wavefunction expansion    Analytical solution    
0 引言

当地震波遇到峡谷、高山等局部地形时会发生散射,散射波、入射波和反射波的相互作用会引起地震动的放大或衰减,称为地形效应(周红等,2010田林等,2012).自从美国圣费尔南多地震记录的异常放大现象(Trifunac and Hudson, 1971)得到关注之后,国内外大量的地震实测记录和现场震害观测进一步证实了地形效应的存在(王海云和谢礼立, 2010Hough et al., 2010李平等, 2016).

为了揭示地形效应机理,地震学和土木工程领域学者开展了地表地形对地震地面运动影响规律的理论研究,包括数值模拟和解析分析两类方法.数值研究众多,例如有限差分法(Boore, 1972)、有限元法(廖振鹏和刘晶波, 1992章小龙等, 2017)、谱元法(李孝波等,2014于彦彦等, 2017贺春晖等, 2017)、边界方法(林皋和关飞, 1990曹军等, 2004巴振宁等, 2017Liu et al., 2018)及混合方法(Zhang and Zhao, 1988金峰等, 1993Zhou and Chen, 2006).数值方法可以模拟任意复杂地形情况,但其计算精度往往需要解析解的校准.解析方法主要是波函数展开法,Trifunac(1973)提出的半圆形峡谷对地震SH波散射的波函数级数解成为早期经典的解析解,随后涌现出更多针对不同地形的解析解,如半椭圆形峡谷(Wong and Trifunac, 1974)、圆弧形峡谷(Yuan and Liao, 1994)、圆弧形饱和土峡谷(李伟华和赵成刚, 2003)、圆弧形层状峡谷(梁建文等, 2003)、圆弧形层状沉积谷(张郁山, 2008, 2010)、去底半圆形峡谷(Tsaur and Chang, 2009)、U形峡谷(Gao et al., 2012)、径向非均质半圆形峡谷(Zhang et al., 2017, 2019)等.这些解析解为数值方法提供了验证基准,并且引领了新型数值方法(Dai et al., 2019)的发展.

由于V形峡谷分布广泛,这类地形对地震波的散射问题得到了大量的数值(Sánchez-Sesma and Rosenblueth, 1979赵崇斌等, 1988)和解析(Tsaur and Chang, 2008Tsaur et al., 2010Zhang et al., 2012a, bGao and Zhang, 2013)研究.在长期的地质演变中,由于峡谷侧壁的地应力集中,这种峡谷最终倾向于演变成含峭壁V形峡谷,典型的例子如亚利桑那州北部科罗拉多河的Glen峡谷(如图 1).然而,含峭壁V形峡谷地形效应研究尚未见报道.实际上,诸多峡谷具有峭壁,有必要开展含峭壁峡谷地形对地震波的散射特性研究.因此,本文采用波函数展开法提出了含峭壁V形峡谷对平面SH波散射的解析解,重点考察了上部峭壁对其地形效应的影响.

图 1 含峭壁V形峡谷示例:亚利桑那州北部Glen峡谷 Fig. 1 An example of V-shaped canyon with cliffs: Glen canyon in northern Arizona
1 模型和理论推导

本文考虑的含峭壁V形峡谷的二维模型如图 2所示,其所在半无限空间假设为均匀、各向同性和线弹性介质,剪切模量μ和剪切波速cs为常数.峡谷的半宽为b(作为标准化因子),峭壁的深度为h1,峡谷的最大深度为h1+h2(< b).考虑到无法用单一坐标系描述这一模型,本文采用先分区建立波场、然后利用分区界面连续条件实现波场匹配的求解策略,即区域分解策略:引入虚拟边界S1,其与上部峭壁共同组成一个半圆.这个半圆将半空间分为开区域①和闭区域②.为了构建两个区域的波场,在两个区域分别设置相应的坐标系.其中,闭区域②中局部坐标系(x1, y1)和(r1, θ1)原点设在峡谷底部,开区域①整体坐标系(x, y)和(r, θ)原点设在水平地表峡谷中点处.峭壁的位置角度可以表示为β=arccos(h1/b).底部侧壁的角度通过计算可以表示为β1=π- .入射波为单位幅值的平面SH波,其入射角为α、振动圆频率为ω.

图 2 平面SH波作用下含峭壁V形峡谷的二维模型 Fig. 2 2D model of a V-shaped canyon with cliffs subjected to a plane SH waves

对于稳态反平面问题,可省略时间因子e-iωt,则模型两个区域的稳态波场uj满足亥姆霍兹方程(Pao and Mow, 1973):

(1)

其中j=1, 2, 分别表示区域①和②,k=ω/cs表示波数.

除了亥姆霍兹方程,两个区域的波场uj还需要满足水平地表、峡谷表面和峭壁表面的应力自由条件:

(2)

(3)

(4)

另外,根据连续介质假设,由虚拟边界S1分割而成的两区域的位移和应力均须满足连续性条件:

(5)

(6)

公式(5)和(6)可保证两区域之间的位移和应力匹配.

区域①的波场包含两部分(自由场和散射场):

(7)

自由场可以表示为

(8)

利用Jacobi-Anger展开公式:

(9)

其中εn是纽曼因子(ε0=1,εn=2, n≥1),Jn(·)是n阶第一类贝塞函数,公式(8)中的自由场uf可以展开成:

(10)

模型物理意义明确,散射波由河谷处产生、向无穷远处传播,为了得到散射波场的唯一解,其应该满足无穷远处Sommerfeld辐射条件:

(11)

通过分离变量法求解公式(1)并排除不满足辐射条件(11)的波场后可以得到:

(12)

其中AnBn为待求系数,Hn(1)(·)是n阶第一类汉克尔函数.

在区域②中,满足峡谷底部应力自由条件(公式(3))的驻波场u2可以表示为

(13)

其中CnDn为待求系数,ν=π/(2β1).

为了便于问题的解决,需要用Graf加法公式将局部坐标系(r1, θ1)中的驻波场u2转化到整体坐标系(r, θ)中.推导得到相应的变换公式为

(14)

其中h =h1+h2.

将公式(14)代入公式(13),驻波场u2可以表示为

(15)

至此,利用前述区域分解策略完成了两区域波场的构建,且两区域波场分别自动满足水平地表应力自由条件(公式(2))和峡谷表面应力自由条件(公式(3)).然而,波场的四组未知系数AnBnCnDn仍然待定,需要利用余下的峭壁应力自由条件(公式(4))和虚拟边界连续条件(公式(5)和(6))完成求解,即区域匹配技术:

首先,补充区域②在峭壁上的零应力条件,获得区域②在与区域①的半圆形交界面上的分段应力函数:

(16)

据此可将峭壁应力自由条件(公式(4))和虚拟边界应力连续条件(公式(6))合并为

(17)

对式(17)两边同时关于θ在区间[-π/2, π/2]积分,并利用三角函数在此区间的正交性,整理可得关系式:

(18)

(19)

其中符号′表示相应函数的微分形式,Im, nC(θ)和Im, nS(θ)分别为

(20)

(21)

然后,利用虚拟边界S1上的位移连续条件(公式(5)),将两区域位移场在区间[-β, β]上积分得:

(22)

(23)

最后,将公式(18)和(19)分别代入公式(22)和(23),得到下面仅和CnDn有关的两个方程组:

(24)

(25)

将各无穷级数进行截断、保留有限项之后即可进行未知系数计算(nmp分别截取NMP项).求解公式(24)和(25)可得到系数CnDn,随后通过公式(18)和(19)求得AnBn.得到这些未知系数后,半空间中的波场即已知.

关于NMP的取值,其中N代表级数解的级数项数量,需要通过级数解收敛测试确定.本文选取了峡谷表面5个代表性位置进行了收敛测试.图 3给出了参数对应α=30°、h1/b=h2/b=1/3和η=4的收敛测试结果.从图 3中可以看出,随着N的逐渐增大,地表位移逐渐收敛至一个固定值.通过多频率下的收敛测试,本文采用N=25即可获得精确的结果.M和P要保证驻波场u2和自由场uf的坐标变换精度,分别利用公式(14)和(9)的满足情况确定其数值为M=150,P=150.

图 3 峡谷表面五个代表性位置的地表位移收敛测试图(α=30°, h1/b=h2/b=1/3, η=4) Fig. 3 Convergence tests at five representative positions on the canyon surface (α=30°, h1/b=h2/b=1/3, η=4)
2 验证

为了方便起见,首先定义本文位移幅值|u|和无量纲频率η

(26)

(27)

其中Re(·)和Im(·)分别表示复位移的实部和虚部,λ表示入射SH波波长.

h1/b=0时,含峭壁的V形峡谷可以退化为V形峡谷;当h2/b=0时,可退化为去底半圆形峡谷.Tsaur和Chang(2008, 2009)分别给出了V形峡谷和去底半圆形峡谷的解析解,因此可以用来验证本文结果的准确性.

图 4将本文的结果与Tsaur和Chang(2008)进行了对比,对应参数为h1/b=0,h2/b=1/2,η=1.图 5将本文的结果与Tsaur和Chang(2009)进行了对比,参数为h1/b=1/2,h2/b=0,η=1.通过4个不同入射角度α=0°、30°、60°和90°下结果的良好吻合情况,验证了本文方法的准确性.

图 4 本文结果(实线,对应参数:h2/b=0,h2/b=1/2,η=1)与Tsaur和Chang(2008)中无峭壁V形峡谷相应结果(点线)的比较 (a) α=0°; (b) α=30°; (c) α=60°; (d) α=90°. Fig. 4 Comparison between the solution of this study (solid lines) when h1/b=0, h2/b=1/2, η=1 and the corresponding results (dotted lines) of the V-shaped canyon without steep cliffs given by Tsaur and Chang (2008)
图 5 本文结果(实线,对应参数h2/b=1/2,h2/b=0,η=1)与Tsaur和Chang(2009)中去底半圆形峡谷的相应结果(虚线)的比较 (a) α=0°; (b) α=30°; (c) α=60°; (d) α=90°. Fig. 5 Comparison between the solution of this study (solid lines) when h1/b=1/2, h2/b=0, η=1 and the corresponding results (dotted lines) of truncated semicircular canyon given by Tsaur and Chang (2009)

为了进一步检验本方法的正确性和实用性,继续采用退化V形峡谷模型模拟中国台湾翡翠峡谷台阵的地震动并与实测结果进行对照.根据Huang和Chiu(1995)对1992年花莲地震时翡翠河谷台阵实测地震动的分析,编号为SC1的台站地震记录相比SC4加速度峰值放大了2.69倍.

翡翠台阵所在的河谷剖面如图 6a所示,利用本文方法为之建立V形峡谷模型(b=563 m,d=300 m,h1=0 m).由于1992年花莲地震入射波数据缺失,本文采用PGA=0.2g的集集地震波作为入射波,利用本文模型求得台站SC1和SC4处的地震动频域传递函数,然后据此得到台站SC1和SC4的地震动加速度时程,如图 6b所示(频域转时域方法参见Dai et al., 2019).其中台站SC1和SC4处的模拟地震动PGA分别为0.3324 g和0.1256 g,即SC1相对于SC4放大2.65倍,与前述实测2.69倍基本一致.可见,本方法可以较好地刻画峡谷对地震动的差异放大规律.

图 6 本文方法对中国台湾翡翠河谷地震动的模拟 (a) 翡翠河谷模型简图(据Huang and Chiu, 1995);(b) 集集地震入射下SC1台站和SC4台站SH波地震动加速度时程模拟结果. Fig. 6 Simulation of the ground motions for the Feitsui Canyon, Taiwan, China (a) Schematic diagram of Feitsui Canyon model (modified from Huang and Chiu, 1995); (b) The calculated SH ground accelerations at positions SC1 and SC4 under the incidence of the Chi-Chi earthquake.
3 数值结果与讨论

本文聚焦上部峭壁对峡谷地形效应的影响,故将参数h2/b设为固定值1/3,从而使得h1/b有较大的变化范围.图 7-图 9分别为对应3个入射波无量纲频率η=0.5、1和4的峡谷地表位移幅值结果,每个无量纲频率结果都包含4个不同入射角度(α=0°, 30°, 60°, 90°),实线结果对应h1/b=0,虚线对应h1/b=1/3,虚实线对应h1/b=2/3.为了方便比较,x/b的计算范围选取为-4到4.

图 7 对应η=0.5和不同入射角度的地表位移幅值结果 (a) α=0°; (b) α=30°; (c) α=60°; (d) α=90°.黑色实线对应h1/b=0,h2/b=1/3;红色虚线对应h1/b=1/3,h2/b=1/3;蓝色点划线对应h1/b=2/3,h2/b=1/3. Fig. 7 Surface displacement amplitudes at four different incident angles and η=0.5 (a) α=0°; (b) α=30°; (c) α=60°; (d) α=90°. h1/b=0, h2/b=1/3 (black solid lines), h1/b=1/3, h2/b=1/3 (red dashed lines), h1/b=2/3, and h2/b=1/3 (blue dotted lines).
图 8 对应η=1和不同入射角度的地表位移幅值结果 (a) α=0°; (b) α=30°; (c) α=60°; (d) α=90°.黑色实线对应h1/b=0,h2/b=1/3;红色虚线对应h1/b=1/3,h2/b=1/3;蓝色点划线对应h1/b=2/3,h2/b=1/3. Fig. 8 Surface displacement amplitudes at four different incident angles and η=1 (a) α=0°; (b) α=30°; (c) α=60°; (d) α=90°. h1/b=0, h2/b=1/3 (black solid lines), h1/b=1/3, h2/b=1/3 (red dashed lines), h1/b=2/3, and h2/b=1/3 (blue dotted lines).
图 9 对应η=4和不同入射角度的地表位移幅值结果 (a) α=0°; (b) α=30°; (c) α=60°; (d) α=90°.黑色实线对应h1/b=0,h2/b=1/3;红色虚线对应h1/b=1/3,h2/b=1/3;蓝色点划线对应h1/b=2/3,h2/b=1/3. Fig. 9 Surface displacement amplitudes at four different incident angles and η=4 (a) α=0°; (b) α=30°; (c) α=60°; (d) α=90°. h1/b=0, h2/b=1/3 (black solid lines), h1/b=1/3, h2/b=1/3 (red dashed lines), h1/b=2/3, and h2/b=1/3 (blue dotted lines).

图 7-图 9的结果可以看出,含峭壁V形峡谷的地表位移幅值在自由场位移幅值2上下波动(|u|>2表示放大,|u| < 2表示减小),峡谷峭壁对地表位移幅值具有重要影响.一方面,峭壁使得地震放大效应更为显著,通过比较h1/b=2/3和h1/b=0两种情况,可以发现含峭壁峡谷的地表位移幅值可达不含峭壁峡谷的2倍.另一方面,随着峭壁深度h1/b从1/3增大至2/3,位移振幅的最大值增大,而位移振幅的最小值减小.总之,峭壁对V形峡谷地形效应具有显著的增强作用,这意味着在地震作用下峭壁峡谷周边的建筑物将会遭受更为剧烈的振动,可能会造成更大的震害.

对于地震波竖向入射的情况(如图 7a图 8a图 9a所示),由于波的聚焦,地表位移幅值的最大值往往发生在峡谷两肩.不同于具有对称性的竖向入射结果,地震波斜入射的情况(图 7(b-d)图 8(b-d)图 9(b-d))展示出左右不对称的地表位移幅值.迎波侧比背波侧的位移幅值波动得更加频繁,随着入射角度的增加,这种不对称表现得更加明显.在小角度斜入射(α=30°)的情况下,迎波侧的位移振幅仅略大于背波侧,而对于大角度斜入射(α=60°, 90°)的情况,两侧差别更加明显.随着无量纲频率η的增大,位移振幅的放大和减小越来越明显,说明峡谷成为一个有效的障碍物,短波穿越峡谷到达对岸的能力越来越弱.在图 7(b-d)图 8(c-d)图 9d中,阴影区的位移幅值曲线趋于平稳,即体现了峭壁峡谷的滤波屏蔽效应.

图 10展示了在三种入射角度下V形峡谷地表位移幅值随着无量纲位置x/b和无量纲频率η的变化情况,其中图 10(a-c)是不含峭壁的情况(h1/b=0)、图 10(d-f)是含峭壁的情况(h1/b=1/3).图中可见,随着无量纲频率的增大,含峭壁峡谷的地形放大形态趋于复杂.在同一入射角下,峭壁的存在使得峡谷两侧表现出更为明显的放大效应,也就是说峭壁的出现对峡谷地形效应有加剧作用.

图 10 不同角度下含峭壁和不含峭壁峡谷地表位移幅值随位置x/b和无量纲频率η的变化图 (a) α=0°,h1=0,h2=1/3; (b) α=45°,h1=0,h2=1/3; (c) α=90°,h1=0,h2=1/3; (d) α=0°,h1=1/3,h2=1/3; (e) α=45°,h1=1/3,h2=1/3; (f) α=90°,h1=1/3,h2=1/3. Fig. 10 Ground surface displacement amplitudes as a function of position x/b and dimensionless frequency η for canyons with and without steep cliffs at four different incident angles
4 结论

本文提出了含峭壁V形峡谷对地震SH波散射的波函数级数解,揭示了含峭壁V形峡谷对地震动的地形放大效应.本文解析模型可退化为不含峭壁的V形峡谷和去底半圆形峡谷两个已有解析模型,通过退化计算验证了本文方法.计算结果表明上部峭壁对峡谷地震动具有重要影响,峭壁峡谷(h1/b=2/3)较无峭壁峡谷(h1/b=0)对地表位移幅值的放大差异可达200%,且较陡的峭壁明显加剧峡谷的地震放大效应,也就增加了周围建筑物的震害风险,值得工程师的关注.本文结论基于二维峡谷场地SH波散射模型,关于峡谷场地的三维效应以及P波和SV波等其他类型地震波的散射特征,尚需进一步探索.

References
Ba Z N, Huang D Y, Liang J W, et al. 2017. Scattering and diffraction of plane SH-waves by periodically distributed hill topographies. Chinese Journal of Geophysics (in Chinese), 60(3): 1039-1052. DOI:10.6038/cjg20170317
Boore D. 1972. Finite difference methods for seismic wave propagation in heterogeneous materials. Methods in Computational Physics: Advances in Research and Applications, 11: 1-37.
Cao J, Ge Z X, Zhang J, et al. 2004. A comparative study on seismic wave methods for multi-layered media with irregular interfaces: irregular topography problem. Chinese Journal of Geophysics (in Chinese), 47(3): 495-503.
Dai D H, Zhang N, Lee V W, et al. 2019. Scattering and amplification of SV waves by a semi-cylindrical hill in a half-space by a wavefunction-based meshless method using mapping and point-matching strategies. Engineering Analysis with Boundary Elements, 106: 252-263. DOI:10.1016/j.enganabound.2019.05.016
Gao Y F, Zhang N. 2013. Scattering of cylindrical SH waves induced by a symmetrical v-shaped canyon: near-source topographic effects. Geophysical Journal International, 193(2): 874-885. DOI:10.1093/gji/ggs119
Gao Y F, Zhang N, Li D Y, et al. 2012. Effects of topographic amplification induced by a U-shaped canyon on seismic waves. Bulletin of the Seismological Society of America, 102(4): 1748-1763. DOI:10.1785/0120110306
He C H, Wang J T, Zhang C H. 2017. Determination of seismic parameters for dam sites by numerical simulation of the rupture-canyon wave field. Chinese Journal of Geophysics (in Chinese), 60(2): 585-592. DOI:10.6038/cjg20170213
Hough S E, Altidor J R, Anglade D, et al. 2010. Localized damage caused by topographic amplification during the 2010 M7.0 Haiti earthquake. Nature Geoscience, 3(11): 778-782. DOI:10.1038/ngeo988
Huang H C, Chiu H C. 1995. The effect of canyon topography on strong ground motion at Feitsui damsite: quantitative results. Earthquake Engineering and Structural Dynamics, 24(7): 977-990. DOI:10.1002/eqe.4290240705
Jin F, Zhang C H, Wang G L. 1993. Seismic response of alluvial soil deposit on a semi-elliptical canyon. Journal of Tsinghua University (Science and Technology) (in Chinese), 33(5): 23-30.
Li P, Liu HS, Bo JS, et al. 2016. Effects of river valley topography on anomalously high intensity in the Hanyuan town during the Wenchuan MS8.0 earthquake. Chinese Journal of Geophysics (in Chinese), 59(1): 174-184. DOI:10.6038/cjg20160115
Li W H, Zhao C G. 2003. An analytical solution for the diffraction of plane P-waves by circular cylindrical canyons in a fluid-saturated porous media half space. Chinese Journal of Geophysics (in Chinese), 46(4): 539-546.
Li X B, Bo J S, Qi W H, et al. 2014. Spectral element method in seismic ground motion simulation. Progress in Geophysics (in Chinese), 29(5): 2029-2039. DOI:10.6038/pg20140506
Liang J W, Zhang Y S, Gu X L, et al. 2003. Scattering of plane SH waves by a circular-arc layered canyon. Journal of Vibration Engineering (in Chinese), 16(2): 158-165.
Liao Z P, Liu J B. 1992. Fundamental problems in finite element simulation of wave motion. Science in China, Ser. B, 35(11): 1353-1364.
Lin G, Guan F. 1990. Scattering of seismic waves at irregular topographies by boundary element method. Journal of Dalian University of Technology (in Chinese), 30(2): 145-152.
Liu Z X, Wang D, Liang J W, et al. 2018. The fast multi-pole indirect BEM for solving high-frequency seismic wave scattering by three-dimensional superficial irregularities. Engineering Analysis with Boundary Elements, 90: 86-99. DOI:10.1016/j.enganabound.2018.02.009
Pao Y H, Mow C C. 1973. Diffraction of Elastic Waves and Dynamic Stress Concentrations. New York: Crane and Russak.
Sánchez-Sesma FJ, Rosenblueth E. 1979. Ground motion at canyons of arbitrary shape under incident SH waves. Earthquake Engineering and Structural Dynamics, 7(5): 441-450. DOI:10.1002/eqe.4290070505
Tian L, Xu Z M, Zhang J M. 2012. Seismic response for local landform subjected to scattering of incident wave. Progress in Geophysics (in Chinese), 27(1): 122-130. DOI:10.6038/j.issn.1004-2903.2012.01.014
Trifunac M D. 1973. Scattering of plane SH waves by a semi-cylindrical canyon. Earthquake Engineering & Structural Dynamics, 1(3): 267-281.
Trifunac M D, Hudson D E. 1971. Analysis of the Pacoima dam accelerogram-San Fernando, California, earthquake of 1971. Bulletin of the Seismological Society of America, 61(5): 1393-1411.
Tsaur D H, Chang K H. 2008. An analytical approach for the scattering of SH waves by a symmetrical V-shaped canyon: shallow case. Geophysical Journal International, 174(1): 255-264. DOI:10.1111/j.1365-246X.2008.03788.x
Tsaur D H, Chang K H. 2009. Scattering of SH waves by truncated semicircular canyon. Journal of Engineering Mechanics, 135(8): 862-870. DOI:10.1061/(ASCE)0733-9399(2009)135:8(862)
Tsaur D H, Chang K H, Hsu M S. 2010. An analytical approach for the scattering of SH waves by a symmetrical V-shaped canyon: deep case. Geophysical Journal International, 183(3): 1501-1511. DOI:10.1111/j.1365-246X.2010.04806.x
Wang H Y, Xie L L. 2010. Effects of topography on ground motion in the Xishan park, Zigong city. Chinese Journal of Geophysics (in Chinese), 53(7): 1631-1638. DOI:10.3969/j.issn.0001-5733.2010.07.014
Wong H L, Trifunac M D. 1974. Scattering of plane SH waves by a semi-elliptical canyon. Earthquake Engineering & Structural Dynamics, 3(2): 157-169.
Yu Y Y, Ding H P, Liu Q F. 2017. Integration of transmitting boundary and spectral-element method and improvement on the accuracy of wave motion simulation. Journal of Vibration and Shock (in Chinese), 36(2): 13-22.
Yuan X M, Liao Z P. 1994. Scattering of plane SH waves by a cylindrical canyon of circular-arc cross-section. Soil Dynamics and Earthquake Engineering, 13(6): 407-412. DOI:10.1016/0267-7261(94)90011-6
Zhang C H, Zhao C B. 1988. Effects of canyon topography and geological conditions on strong ground motion. Earthquake Engineering & Structural Dynamics, 16(1): 81-97.
Zhang N, Gao Y F, Cai Y Q, et al. 2012a. Scattering of SH Waves Induced by a Non-symmetrical V-shaped Canyon. Geophysical Journal International, 191(1): 243-256. DOI:10.1111/j.1365-246X.2012.05604.x
Zhang N, Gao Y F, Li D Y, et al. 2012b. Scattering of SH waves induced by a symmetrical V-shaped canyon: a unified analytical solution. Earthquake Engineering and Engineering Vibration, 11(4): 445-460. DOI:10.1007/s11803-012-0135-z
Zhang N, Gao Y F, Pak R Y S. 2017. Soil and topographic effects on ground motion of a surficially inhomogeneous semi-cylindrical canyon under oblique incident SH waves. Soil Dynamics and Earthquake Engineering, 95: 17-28. DOI:10.1016/j.soildyn.2017.01.037
Zhang N, Zhang Y, Gao Y F, et al. 2019. An exact solution for SH-wave scattering by a radially multilayered inhomogeneous semicylindrical canyon. Geophysical Journal International, 217(2): 1232-1260. DOI:10.1093/gji/ggz083
Zhang X L, Li X J, Zhou Z H, et al. 2017. The seismic response analysis of three-dimensional canyon complex topography under incident SV seismic waves. Chinese Journal of Geophysics (in Chinese), 60(7): 2779-2790. DOI:10.6038/cjg20170723
Zhang Y S. 2008. Analytical solution for the stationary response of alluvial valleys containing multiple circular-arc layers to incident plane P waves. Chinese Journal of Geophysics (in Chinese), 51(3): 869-880.
Zhang Y S. 2010. Analytical solution to dynamic response of circular-arc-shaped multi-layered valley due to incidence of Rayleigh wave. Chinese Journal of Geophysics (in Chinese), 53(9): 2129-2143.
Zhao C B, Zhang C H, Zhang G D. 1988. Analysis of free field motions on arch dam canyons and the effects due to existence of weathered rock strata. Journal of Hydraulic Engineering (in Chinese), (11): 29-39.
Zhou H, Chen X F. 2006. A new approach to simulate scattering of SH waves by an irregular topography. Geophysical Journal International, 164(2): 449-459. DOI:10.1111/j.1365-246X.2005.02670.x
Zhou H, Gao M T, Yu Y X. 2010. A study of topographical effect on SH waves. Progress in Geophysics (in Chinese), 25(3): 775-782. DOI:10.3969/j.issn.1004-2903.2010.03.005
巴振宁, 黄棣旸, 梁建文, 等. 2017. 层状半空间中周期分布凸起地形对平面SH波的散射. 地球物理学报, 60(3): 1039-1052. DOI:10.6038/cjg20170317
曹军, 盖增喜, 张坚, 等. 2004. 用于不规则界面多层介质的地震波方法比较研究: 不规则地形问题. 地球物理学报, 47(3): 495-503. DOI:10.3321/j.issn:0001-5733.2004.03.020
贺春晖, 王进廷, 张楚汉. 2017. 基于震源-河谷波场数值模拟的坝址地震动参数确定方法. 地球物理学报, 60(2): 585-592. DOI:10.6038/cjg20170213
金峰, 张楚汉, 王光纶. 1993. 半椭圆形河谷上沉积层地震响应研究. 清华大学学报(自然科学版), 33(5): 23-30.
李平, 刘红帅, 薄景山, 等. 2016. 汶川MS8.0地震河谷地形对汉源县城高烈度异常的影响. 地球物理学报, 59(1): 174-184. DOI:10.6038/cjg20160115
李伟华, 赵成刚. 2003. 圆弧形凹陷饱和土场地对平面P波散射问题的解析解. 地球物理学报, 46(4): 539-546. DOI:10.3321/j.issn:0001-5733.2003.04.017
李孝波, 薄景山, 齐文浩, 等. 2014. 地震动模拟中的谱元法. 地球物理学进展, 29(5): 2029-2039. DOI:10.6038/pg20140506
梁建文, 张郁山, 顾晓鲁, 等. 2003. 圆弧形层状凹陷地形对平面SH波的散射. 振动工程学报, 16(2): 158-165. DOI:10.3969/j.issn.1004-4523.2003.02.006
廖振鹏, 刘晶波. 1992. 波动有限元模拟的基本问题. 中国科学B辑, 22(8): 874-882. DOI:10.3321/j.issn:1006-9240.1992.08.004
林皋, 关飞. 1990. 用边界元法研究地震波在不规则地形处的散射问题. 大连理工大学学报, 30(2): 145-152.
田林, 徐则民, 张家明. 2012. 局部地形下入射波散射效应对场地地震响应的影响. 地球物理学进展, 27(1): 122-130. DOI:10.6038/j.issn.1004-2903.2012.01.014
王海云, 谢礼立. 2010. 自贡市西山公园地形对地震动的影响. 地球物理学报, 53(7): 1631-1638. DOI:10.3969/j.issn.0001-5733.2010.07.014
于彦彦, 丁海平, 刘启方. 2017. 透射边界与谱元法的结合及对波动模拟精度的改进. 振动与冲击, 36(2): 13-22.
章小龙, 李小军, 周正华, 等. 2017. 三维复杂山谷地形SV波垂直输入地震反应分析. 地球物理学报, 60(7): 2779-2790. DOI:10.6038/cjg20170723
张郁山. 2008. 圆弧状多层沉积谷地在平面P波入射下稳态响应的解析解. 地球物理学报, 51(3): 869-880. DOI:10.3321/j.issn:0001-5733.2008.03.028
张郁山. 2010. 圆弧状多层沉积谷地在Rayleigh波入射下动力响应的解析解. 地球物理学报, 53(9): 2129-2143. DOI:10.3969/j.issn.0001-5733.2010.09.013
赵崇斌, 张楚汉, 张光斗. 1988. V型拱坝峡谷地震自由场分析及地表风化层的影响. 水利学报, (11): 29-39. DOI:10.3321/j.issn:0559-9350.1988.11.004
周红, 高孟潭, 俞言祥. 2010. SH波地形效应特征的研究. 地球物理学进展, 25(3): 775-782. DOI:10.3969/j.issn.1004-2903.2010.03.005