Analytical study on the topographic effect on ground motion of Feitsui canyon
-
摘要: 地表不规则地形的存在往往会引起地震波的散射，进而产生局部地震动放大或衰减的现象。虽然地形效应最早在异常地震记录中被发现，然而利用地形影响台阵记录到的地震动数据却少之又少。基于1992年在我国台湾翡翠河谷上观测到的六条地形影响台阵记录，利用线源SH波入射下非对称V形河谷地震波传播解析理论，模拟得到了河谷台阵各点的地震动，揭示了翡翠河谷地形效应对地震动影响的规律。与实测记录的对比结果表明利用笔者提出的地震波传播解析理论能够较好地反映河谷地形对地震动的影响。Abstract: The local amplification or attenuation of ground motion is often caused by the surface irregular topography, which is due to the scattering phenomenon when the seismic wave propagates to the local topography. Although the topographic effect was first discovered in the anomalous seismic records, the ground motion data recorded by the topographic influence array are very few. There are five earthquakes recorded by six strong-motion accelerometers deployed along the Feitsui canyon recorded from the Hualien earthquake in 1992. Using our analytical theory of seismic wave propagation around a non-symmetrical V-shaped canyon excited by a line source of SH waves, the strong motion accelerations of the site were simulated. The influences of topographic effect on the ground motion of the Feitsui canyon are revealed. The comparison between the simulated results and those ground motion accelerations recorded indicates that the proposed analytical theory of seismic wave propagation around V-shaped canyon is suitable for simulating the topographic effects on ground motions.
-
Keywords:
- Feitsui canyon /
- topographic effect /
- analytical solution /
- V-shaped canyon
-
引言
地震波传播特性与地震动输入是重大工程地震响应机制与灾变机理和抗震设计理论研究的基础，由于涉及地震学与工程学的交叉，这一基础研究滞后于地震响应分析，因此一直是重大工程抗震设计的瓶颈。对于局部不规则场地，其地震波传播特性极其复杂，表现为地形效应十分突出（郭明珠等，2013；高玉峰等，2021）。地震波传播至局部不规则场地后由于散射而导致场地地震动发生放大、衰减以及空间变化，而这种地形效应往往加剧桥梁、隧道和大坝等工程的震害。以2008年汶川地震为例，该地震由于发生在我国西部山区，相较于唐山地震这类平原地区的地震，其震害具有非常显著的差异，该地震造成6 140座桥梁、156条隧道和1 996座水坝受损（杜修力等，2008；景立平等，2009），可见局部地区地形条件的影响不容忽视。
地震波地形效应最早始于对实际地震动记录的分析，在1971年美国圣弗尔南多（San Fernando）地震中，Trifunac和Hudson （1971）发现帕克伊马（Pacoima）大坝上的台站观测数据出现异常的地震动放大现象。为了进一步研究地形效应，世界上很多地震台阵得以布设。1991年我国台湾省翡翠河谷上布置了六台强震仪，1992年9月花莲ML5.5地震触发了此台阵，获得了沿地形分布的地面运动加速度记录，结果显示河谷两侧的地震动差异明显，表现出了明显的地形效应（Huang，Chiu，1995）。Spudich等（1996）在研究加州北岭（Northridge）地震记录时发现，相对于山脚处，山顶出现了明显的地震动放大效应。Parolai等（2004）通过对比分析1999年土耳其伊兹米特（Izmit）地震的17条观测数据，发现复杂的地形条件下会存在高达5倍的地震动放大效应。基于2008年汶川地震期间自贡西山地形影响台阵的观测数据，王海云和谢礼立（2010）详细分析了地形及地质条件对地震动的影响，其结果显示山体地形放大效应随高程的增加而愈加明显。此后，杨宇等（2011）和唐晖等（2012）也针对自贡台阵观测数据进行了深入研究。王文才等（2020）针对2018年陕西宁强MS5.3地震的地震动数据进行分析，其结果显示了山体对地震中频段的地震动放大效应。
为了揭示地形效应的机理，关于地形对地震波的散射和衍射的数值方法和解析研究得以大量开展，数值方法主要包括有限元法（廖振鹏，刘晶波，1992；章小龙等，2017；孙纬宇等，2019）、有限差分法（Boore，1972）、边界元法（林皋，关飞，1990；巴振宁等，2017；Liu et al，2018；何颖等，2019）、谱元法（贺春晖等，2017；于彦彦等，2017；刘中宪等，2020）等。虽然数值方法应用灵活，但是其结果的准确性需要解析研究的验证。解析研究主要是波函数展开法，该方法不仅可以揭示波散射问题的物理本质，还可以检验各种数值方法的精度，近年来得到广泛关注（Trifunac，1973；Yuan，Liao，1994；Tsaur，Chang，2009；Gao et al，2012；Zhang et al，2017，2019；张宁等，2021）。
V形河谷作为常见的河谷形态，由于地震作用下河谷底部会产生应力奇异问题，其地形效应的理论解一直未得到有效解决，直到最近才由Tsaur教授带领的研究小组取得突破，分别得到了浅河谷（Tsaur，Chang，2008）和深河谷（Tsaur et al，2010）的地震SH波散射问题解答。Zhang等（2012a）为了应用方便，给出了不同深度河谷地震波散射的统一解答，随后Zhang等（2012b）研究了非对称V形河谷的地震SH波放大效应，并将此问题拓展至近源地震动的情况（Gao，Zhang，2013）。针对河谷场地上的地震动地形观测台阵，目前仅我国台湾省翡翠河谷上建立的六台强震仪在1992年花莲地震中记录到了相应的数据（Huang，Chiu，1995），且由于此河谷呈现典型的V形形态，本文将运用V形河谷场地地震SH波散射问题的频域解答，模拟其监测点位置处的地震动，进而从理论上揭示地形对地震动响应的影响规律。
1. 翡翠河谷地形及台阵记录
翡翠河谷位于我国台湾省花莲县，其上建有翡翠大坝，是为台北地区450万人供水的水源工程。在1977年到1990年间，科研人员在翡翠大坝附近分别布设了三个SMA-1型强震仪，主要用于观测河谷地形对翡翠大坝附近地震动在空间分布上的影响规律。但是由于台站距离翡翠大坝较近，其观测结果会受到大坝的影响，Huang和Chiu （1995）又于1991年在距离大坝300 m处沿翡翠河谷横断面布置了六台强震仪（SC1—SC6，高程分别为170，120，70，70，120，160 m），台站的具体位置如图1所示。
![]()
图 1 翡翠河谷地震动观测台阵剖面及台站分布示意图（Huang，Chiu，1995）Figure 1. Definition sketch for the cross section of Feitsui canyon and the location of six stations （Huang，Chiu，1995）1992年9月的花莲ML5.5地震触发了此台阵，由此获得了沿地形分布的地面运动加速度记录，此次地震的震源深度为50 km，震中位于监测台阵南偏西方向距台阵130 km处，震中至翡翠大坝的方向基本平行于台阵观测面。Huang和Chiu （1995）根据实测地震动分离出了垂直于监测剖面方向的地震动分量，即二维SH波地震动记录，为后续地震动模拟提供了便利。台阵观测数据显示出显著的地震动地形效应，各点峰值加速度（peak ground acceleration，简写为PGA）记录见表1 （Huang，Chiu，1995）。由表1可以看出，河谷迎波面相对于背波面表现出较强的地震动放大现象，迎波侧PGA最大可达9.7 cm/s2，而背波侧最大只有6.0 cm/s2；同时可以看出越靠近河谷底部，记录到的地震动PGA越小。为了分析这一现象产生的原因，本文将从V形河谷周围地震波传播解析理论出发得到各点的地震动模拟结果，以揭示其地震响应规律。
表 1 1992年9月花莲地震中SC1—SC6台站记录到的地震动PGA （Huang，Chiu，1995）Table 1. PGA values at stations SC1−SC6 during Hualien earthquake in September 1992 （Huang，Chiu，1995）台站名 PGA/（cm·s−2） 台站名 PGA/（cm·s−2） 台站名 PGA/（cm·s−2） SC1 9.7 SC2 6.5 SC3 5.2 SC4 3.6 SC5 4.5 SC6 6.0 2. 理论模拟
2.1 频域内地震动场推导
由翡翠河谷剖面图可以看到，河谷形状大致呈V形，地震波自剖面左侧斜入射。本文将首先推导线源SH波入射下非对称V形河谷地震波传播理论，得到频域内考虑地形效应的翡翠河谷各点的地震动传递函数。由于翡翠河谷属于深度较浅的河谷，因此本文将针对浅的非对称V形河谷地震波传播理论进行推导。河谷的二维简化模型如图2所示，此图表示半空间中一个非对称V形河谷，深度为d，半宽为b1 （左侧）和b2 （右侧）；模型的介质假设为弹性、各向同性、均质；介质的剪切模量为G和剪切波速为vS。入射波假设为简谐的线源柱面SH波，圆频率为ω，位移在z方向（出平面）。在河谷底部建立局部坐标系（r1，θ1），在河谷中点建立整体坐标系（r，θ），采用半径为（b1＋b2）/2的半圆形辅助边界把整个空间分成两个子区域①和②。震源位置在整体极坐标系（r，θ）中为（r0，θ0），对于本研究的问题，需使用镜像法来考虑相对于水平地面对称的两个波源，对于波源及其镜像分别使用两个极坐标系（rf，θf）和（
${{r'}_{\rm{f}}},\;\;{{\theta '} _{\rm{f}}}$ ）。通过四个坐标系的建立，引入半圆形辅助边界进行区域分解，得到区域①和区域②的运动方程（Sanchez-Sesma，1985；Gao，Zhang，2013）：
$$ \frac{{{\partial ^2}{u_1}}}{{\partial r_1^2}} {\text{＋}} \frac{1}{{{r_1}}}\frac{{\partial {u_1}}}{{\partial {r_1}}} {\text{＋}} \frac{1}{{r_1^2}}\frac{{{\partial ^2}{u_1}}}{{\partial \theta _1^2}} {\text{＋}} {k^2}{u_1} {\text{＝}} 0 {\text{，}} $$ (1) $$ \frac{{{\partial ^2}{u_2}}}{{\partial {r^2}}} {\text{＋}} \frac{1}{r}\frac{{\partial {u_2}}}{{\partial r}} {\text{＋}} \frac{1}{{{r^2}}}\frac{{{\partial ^2}{u_2}}}{{\partial {\theta ^2}}} {\text{＋}} {k^2}{u_2} {\text{＝}} - \frac{{{\text{δ}} （r {\text{－}} {r_0}）{\text{δ}} （\theta {\text{－}} {\theta _0}）}}{r}{\text{，}} $$ (2) 式中，
${\text{δ}} （ \cdot ）$ 为狄拉克δ函数，k＝ω/vS为剪切波数。对于稳态出平面问题，u1和u2分别表示区域①和区域②在频域内的位移。引入镜像法将入射波源镜像至（
${{r}'}_{\rm{f}}{\text{，}}{{\theta }'} _{\rm{f}}$ ）坐标系，通过波函数展开法可以得到区域①和区域②内含未知系数的位移表达式（Zhang et al，2012b）：$$ {u_1}({r_1}{\text{，}}{\theta _1}) {\text{＝}} \sum\limits_{n {\text{＝}} 0}^\infty {{A_n}{{\text{J}}_{n\upsilon }}(k{r_1})\cos \left[ {n\upsilon （ {{\theta _1} {\text{＋}} {\beta _1}} ）} \right]}{\text{，}} $$ (3) $$ {u_2}（r{\text{，}}\theta ） {\text{＝}} {u^{\rm{f}}} {\text{＋}} \sum\limits_{n {\text{＝}} 0}^\infty {{B_n}{\text{H}}_{2n}^{（1）}（kr）\cos （2n\theta ）} {\text{＋}} \sum\limits_{n {\text{＝}} 0}^\infty {{C_n}{\text{H}}_{2n {\text{＋}} 1}^{（1）}（kr）\sin [（2n + 1）\theta ]} {\text{，}} $$ (4) 式中：
${{\rm{J}}_n}（\cdot ）$ 为n阶第一类贝塞尔函数，${\rm{H}}_n^{（1）}（ \cdot ）$ 为n阶第一类汉克尔函数；$\;{\beta _1} {\text{＝}} \pi {\text{－}}\arctan （{b_1}/d）$ ，$\;{\beta _2} {\text{＝}} \pi {\text{－}}\arctan （{b_2}/d ）$ ，$ \upsilon $ ＝ π/（β1＋β2）；${u^{\rm{f}}}{\text{＝}}[ {{\rm{H}}_0^{(1)}(k{r_{\rm{f}}})＋{\rm{H}}_0^{(1)}(k{{r'}_{\rm{f}}})} ]/{\rm{H}}_0^{（1）}（k{r_0}）$ ，表示自由场。上述位移表达式需要满足河谷表面的应力自由边界条件和水平地表的应力自由边界条件：
$$ \tau _{{\theta _1}{{\textit{z}}_1}}^{（1）} {\text{＝}} \frac{G}{{{r_1}}}\frac{{\partial {u_1}}}{{\partial {\theta _1}}} {\text{＝}} 0 \quad\quad {\theta _1} {\text{＝}} - {\beta _1}{\text{，}}{\beta _2} {\text{，}} $$ (5) $$ \tau _{\theta {\textit{z}}}^{（2）} {\text{＝}} \frac{G}{r}\frac{{\partial {u_2}}}{{\partial \theta }} {\text{＝}} 0 \quad\quad \theta {\text{＝}} {\text{±}} \frac{\pi }{2} {\text{，}} r {\text{≥}} \frac{{{b_{\text{1}}}{\text{ + }}{b_{\text{2}}}}}{2} {\text{，}} $$ (6) 同时满足两个子区域之间的位移连续和应力连续条件：
$$ {u_1}（r{\text{，}}\theta ） {\text{＝}} {u_2}（r{\text{，}}\theta） \quad\quad r {\text{＝}} \frac{{{b_{\text{1}}}{\text{＋}}{b_{\text{2}}}}}{2} {\text{，}} - \frac{\pi }{2} {\text{≤}} \theta {\text{≤}} \frac{\pi }{2} {\text{，}} $$ (7) $$ \tau _{r{\textit{z}}}^{（1）}（r{\text{，}}\theta ） {\text{＝}} \tau _{r{\textit{z}}}^{（2）}（r{\text{，}}\theta ） \quad\quad r {\text{＝}} \frac{{{b_{\text{1}}}{\text{＋}}{b_{\text{2}}}}}{2} {\text{，}} - \frac{\pi }{2} {\text{≤}} \theta {\text{≤}} \frac{\pi }{2} {\text{.}} $$ (8) 借助区域匹配策略和Graf加法公式，通过恰当的求解方案，最终可以得到位移场各未知系数An，Bn和Cn的解答，进而各区域的位移即可根据方程（3）和（4）精确得到。具体求解方案可参考Gao和Zhang （2013）。
2.2 地震动时程模拟
地震动记录多数是在地表面实测的加速度时程，测点一般位于不受任何结构影响的自由场地。理论上，对于地震SH波，自由场地加速度幅值是地下基岩地震动幅值的两倍，因此，为了得到结构抗震所需的输入地震波，对于平坦场地，将地表地震动记录的幅值减半作为输入的基岩地震动；而对于河谷场地，需要利用我们提出的地震动传播理论得到地形放大因子（相对于平坦场地）来获取场地地震动。
为了在给定地震动加速度时程输入的情况下得到河谷场地的空间变化地震动，需建立考虑地形放大效应的河谷场地地震动时程合成方法。这一方法建立的关键在于理解并应用河谷引起的地形放大因子，说明如下：V形河谷的波函数级数解以单位幅值稳态SH波作为入射波，如果河谷不存在，不同位置的水平地表位移幅值恒相等，为自由场|uf|；而河谷地形效应导致不同位置的地表位移幅值在自由场|uf|上下波动，也就是说，在某一频率下，如果地表某一点的位移幅值大于自由场|uf|，那么这一点的地面运动相对于地表自由场放大，反之，如果小于自由场|uf|，地面运动减小。根据这一原理，将这些波函数级数解给出的河谷地表位移除以自由场|uf|，即可得到河谷地形的地震动放大因子，该放大因子为复数（包含幅值放大因子和相位调整因子）。因此如果得到了上述河谷地形的地震动放大因子随频率变化的函数，则可根据基岩地震动得到场地某点的实际地震动，具体步骤参见代登辉（2019）。
2.3 翡翠河谷测点地震动模拟
由于监测台阵距震中130 km，震源深度为50 km。因此震源相对于河谷的位置为（x0，y0）＝（−130 km，50 km）。根据Huang和Chiu （1995）的数据，翡翠河谷的模型参数为：深度d＝350 m，宽度b1 ≈ 800 m，b2 ≈ 700 m，介质密度ρ＝2.67 g/cm3，剪切波速vS＝1500 m/s。
运用频域内地震动场的推导公式，可以得到河谷场地目标位置的频域地震动放大因子，地震动放大因子为复数，包含幅值和相位的信息。图3给出了台站SC1—SC6所在位置的地震动放大因子随频率变化的曲线，计算的频率范围为0—10 Hz。从图中可以明显地看出，在0—10 Hz内，SC1—SC3台站处的地震动放大因子最大可达1.5左右，而SC4—SC6台站处的地震动放大因子均小于1。尽管在某些频率（如5 Hz左右）下，迎波侧的地震动放大因子会小于背波侧，但是迎波侧整体表现出相对于背波侧的地形放大现象，这也与实测地震动结果相吻合。
![]()
图 3 翡翠河谷上SC1—SC6台站所在位置的地震动放大因子随频率f的变化Figure 3. Variation of ground motion amplification factors as a function of frequency at the stations SC1−SC6 in Feitsui canyon对于翡翠河谷台阵记录到的1992年9月的花莲地震，其主要频率在5 Hz以内，且主要集中在2—4 Hz，持时约30 s。由于缺乏相应的地震记录数据，为了更可靠地模拟该场地的地震响应，本文通过对Huang和Chiu （1995）一文中的数据进行读点，并依据实测数据的频谱结构适当调整，构造出符合上述频谱特征的人工地震动时程，并以SC3点作为地震动输入参考点，计算其它各点的地震响应。地震动输入参考点SC3的地震动时程如图4a所示，持时为30 s，PGA为5.2 cm/s2，相应的傅里叶谱如图4b所示，可见所构造的地震波与实测地震动波形基本一致，主要频率集中在2—4 Hz。
![]()
图 4 输入参考点SC3台站的地震动时程（a）和相应的傅里叶谱（b）Figure 4. Ground motion history （a） and corresponding Fourier spectrum （b） at the input reference point SC3基于输入地震动和频域地震动放大因子，可以进一步得到河谷各点的地震动响应时程，关键步骤解释如下：
1） 将输入地震动加速度时程记录进行两次数值积分，得到位移时程uinp（t），然后将位移时程进行傅里叶变换，得到其傅里叶谱
$$ {u_{{\rm{inp}}}}（ \omega ） {\text{＝}} \int_{ - \infty }^{ {\text{＋}} \infty } {{u_{{\rm{inp}}}}（t）{{\rm{e}}^{ - {\rm{i}}\omega t}}{\rm{d}}t} {\text{；}} $$ (9) 2） 根据频域内地震动场的推导，求解不同频率简谐振动下河谷场地各个位置的地形放大因子，得到河谷场地目标位置的频域放大传递函数A（ω）；
3） 根据频域放大因子调整目标位置地震动的傅里叶谱
$$ {u_{{\rm{outp}}}}（\omega ） {\text{＝}} A（\omega ）{u_{{\rm{inp}}}}（\omega ）{\text{；}} $$ (10) 4） 根据调整后的傅里叶谱，通过傅里叶逆变换来求解考虑地形效应的地震动时程，即
$$ {u_{{\rm{outp}}}}（t） {\text{＝}} \frac{1}{{2\pi }}\int_{ - \infty }^{ {\text{＋}} \infty } {{u_{{\rm{outp}}}}（\omega ）{{\rm{e}}^{{\rm{i}}\omega t}}{\rm{d}}t} {\text{，}} $$ (11) 将位移时程uoutp（t）进行一次数值微分后得到速度时程，进行两次数值微分后得到加速度时程。
通过计算得到的地震动放大因子和地震动时程模拟方法得到各测点位置处的地震动加速度时程，如图5所示。表2给出了模拟得到的地震动PGA与监测结果的对比。从图5和表2可以看出，SC1—SC3点的地震动PGA明显大于SC4—SC6点，这说明模拟得到的各测点地震动地形分布规律与观测结果一致，即翡翠河谷迎波侧的地震动明显大于背波侧。这是因为，河谷的存在阻挡了地震波的传播，这使得地震波能量大部分集中在河谷迎波面，河谷背波面的地震运动主要是源于河谷引起的散射波的作用。同时可以看出：在河谷迎波面，SC1点的地震动PGA大于SC2点，SC3点最小；在河谷背波面，SC6点的地震动PGA大于SC5点，SC4点最小。这与观测结果相一致，同时可以得出：对于河谷同一坡面，随着所处高程的增加，其地震动也相应增大。从表2中可以看出，以SC3点作为地震动输入参考点，模拟得到的SC2和SC4点的地震动PGA与实测结果最为接近，SC1，SC5和SC6点的模拟结果与实测结果有些许差别。这是因为地震动不只受地形因素的影响，还会受很多其它因素的干扰，如测点所处位置的地质差异、河谷底部沉积层的影响等，所以离地震动输入参考点SC3越远，地震动受其它因素的干扰越大，其结果差别越明显。但是，这种差异不影响结果的指导性，从整体上看模拟结果能够较好地反映河谷场地表现出的地震动地形效应。
![]()
图 5 翡翠河谷上SC1—SC6台站的模拟加速度时程Figure 5. Response acceleration time histories at the stations SC1−SC6 in Feitsui canyon表 2 SC1−SC6台站的模拟地震动PGA与监测结果对比Table 2. Comparisons of PGA values at the stations SC1−SC6 between the records and simulation台站名 监测PGA
/（cm·s−2）模拟PGA
/（cm·s−2）台站名 监测PGA
/（cm·s−2）模拟PGA
/（cm·s−2）台站名 监测PGA
/（cm·s−2）模拟PGA
/（cm·s−2）SC1 9.7 7.8 SC3 5.2 5.2 SC5 4.5 3.6 SC2 6.5 6.3 SC4 3.6 3.1 SC6 6.0 3.9 图6给出了计算得到的各监测点位置的地震动傅里叶谱幅值，可以看出地形效应对地震波的主要频率范围影响较小，但对其傅里叶谱幅值的影响较大。SC1—SC3点的傅里叶谱幅值基本能反映其时域情况，SC1点的傅里叶谱幅值整体大于SC2点和SC3点，SC4—SC6点的傅里叶谱幅值基本相同，因此其在时域内的加速度时程也相差较小，三者的地震动PGA较为接近。由频谱特征分析可知：SC1—SC3点的地震波主要频率集中在2—4 Hz，在该区间内傅里叶频谱幅值明显大于背波侧SC4—SC6点的傅里叶频谱幅值；而SC4—SC6点的地震波主要频率集中在3—5 Hz；在4—5 Hz的频率范围内，六个观测点的傅里叶谱幅值几乎保持不变。这是因为在2—4 Hz内，SC1—SC3点相对于SC4—SC6点产生较大的地形放大，而在4—5 Hz内，六个观测点的地形放大因子差别不大，这一规律同样可以在图3中看到。
![]()
图 6 模拟得到的翡翠河谷上SC1—SC6各点加速度傅里叶谱幅值Figure 6. Response Fourier amplitudes of acceleration at the stations SC1−SC6 in Feitsui canyon3. 讨论与结论
翡翠河谷地形影响地震动观测台阵在1992年花莲地震中获得了相应的地形影响地震动记录，通过分析河谷各监测点的地震动可知，河谷迎波面相对于背波面表现出明显的地震动放大现象。为了分析此现象产生的原因，本文利用线源SH波入射下非对称V形河谷地震波传播解析理论，通过构造频段接近实际观测数据的地震波，并以SC3点作为地震动输入参考点，模拟得到了河谷各监测点位置的地震动。模拟结果进一步揭示了河谷地形效应对地震动影响的规律：河谷地形使得地震波在传播过程中产生散射，导致河谷场地表现出了与平坦场地不同的地震动规律。在河谷迎波侧，散射波与入射波相长干涉，表现出了地震动放大现象；在河谷背波侧，河谷对剪切波产生了过滤作用，表现出了地震动衰减现象。结果进一步表明，河谷同侧地震动表现出了地震动PGA随高程增加而增加的趋势。通过对比实测地震动与模拟结果可知，本文提出的河谷场地地震波传播解析理论能够较好地模拟地形效应引起的地震动差异情况，河谷两侧的地震动差异以及所处高程对地震动的影响均能够通过本解析理论得到较好的模拟。本解析理论不同于其它数值方法，能够实现河谷场地全域全过程的地震动模拟，同时本文提出的研究思路实现了解析理论的工程应用，本文中考虑地形效应的地震动模拟方法可为V形河谷场地重大工程的抗震设计提供精确的地震动输入，所得结果对地震区划的精细划分有一定的指导意义。
-
![]()
图 1 翡翠河谷地震动观测台阵剖面及台站分布示意图（Huang，Chiu，1995）
Figure 1. Definition sketch for the cross section of Feitsui canyon and the location of six stations （Huang，Chiu，1995）
![]()
图 3 翡翠河谷上SC1—SC6台站所在位置的地震动放大因子随频率f的变化
Figure 3. Variation of ground motion amplification factors as a function of frequency at the stations SC1−SC6 in Feitsui canyon
![]()
图 4 输入参考点SC3台站的地震动时程（a）和相应的傅里叶谱（b）
Figure 4. Ground motion history （a） and corresponding Fourier spectrum （b） at the input reference point SC3
![]()
图 5 翡翠河谷上SC1—SC6台站的模拟加速度时程
Figure 5. Response acceleration time histories at the stations SC1−SC6 in Feitsui canyon
![]()
图 6 模拟得到的翡翠河谷上SC1—SC6各点加速度傅里叶谱幅值
Figure 6. Response Fourier amplitudes of acceleration at the stations SC1−SC6 in Feitsui canyon
表 1 1992年9月花莲地震中SC1—SC6台站记录到的地震动PGA （Huang，Chiu，1995）
Table 1 PGA values at stations SC1−SC6 during Hualien earthquake in September 1992 （Huang，Chiu，1995）
台站名 PGA/（cm·s−2） 台站名 PGA/（cm·s−2） 台站名 PGA/（cm·s−2） SC1 9.7 SC2 6.5 SC3 5.2 SC4 3.6 SC5 4.5 SC6 6.0 
下载: 导出CSV
表 2 SC1−SC6台站的模拟地震动PGA与监测结果对比
Table 2 Comparisons of PGA values at the stations SC1−SC6 between the records and simulation
台站名 监测PGA
/（cm·s−2）模拟PGA
/（cm·s−2）台站名 监测PGA
/（cm·s−2）模拟PGA
/（cm·s−2）台站名 监测PGA
/（cm·s−2）模拟PGA
/（cm·s−2）SC1 9.7 7.8 SC3 5.2 5.2 SC5 4.5 3.6 SC2 6.5 6.3 SC4 3.6 3.1 SC6 6.0 3.9
下载: 导出CSV
-
巴振宁,黄棣旸,梁建文,张艳菊. 2017. 层状半空间中周期分布凸起地形对平面SH波的散射[J]. 地球物理学报,60(3):1039–1052. doi: 10.6038/cjg20170317 Ba Z N,Huang D Y,Liang J W,Zhang Y J. 2017. Scattering and diffraction of plane SH-waves by periodically distributed hill topographies[J]. Chinese Journal of Geophysics,60(3):1039–1052 (in Chinese).
代登辉. 2019. 高山峡谷场地P-SV波散射放大效应及地震动场模拟方法研究[D]. 南京: 河海大学: 77–79. Dai D H. 2019. Scattering and Amplification Effect of P-SV Waves and Simulation of Earthquake Ground Motion Field in Alpine Canyon Sites[D]. Nanjing: Hohai University: 77–79 （in Chinese）.
杜修力,韩强,李忠献,李立云,陈少峰,赵建锋. 2008. 5·12汶川地震中山区公路桥梁震害及启示[J]. 北京工业大学学报,34(12):1270–1279. Du X L,Han Q,Li Z X,Li L Y,Chen S F,Zhao J F. 2008. The seismic damage of bridges in the 2008 Wenchuan earthquake and lessons from its damage[J]. Journal of Beijing University of Technology,34(12):1270–1279 (in Chinese).
高玉峰,代登辉,张宁. 2021. 河谷地形地震放大效应研究进展与展望[J]. 防灾减灾工程学报,41(4):734–752. Gao Y F,Dai D H,Zhang N. 2021. Progress and prospect of topographic amplification effects of seismic wave in canyon sites[J]. Journal of Disaster Prevention and Mitigation Engineering,41(4):734–752 (in Chinese).
郭明珠,赵芳,赵凤仙. 2013. 场地地震动局部地形效应研究进展[J]. 震灾防御技术,8(3):311–318. doi: 10.3969/j.issn.1673-5722.2013.03.010 Guo M Z,Zhao F,Zhao F X. 2013. A review of the effect of small-scale surface topography on ground motions[J]. Technology for Earthquake Disaster Prevention,8(3):311–318 (in Chinese).
贺春晖,王进廷,张楚汉. 2017. 基于震源-河谷波场数值模拟的坝址地震动参数确定方法[J]. 地球物理学报,60(2):585–592. doi: 10.6038/cjg20170213 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[J]. Chinese Journal of Geophysics,60(2):585–592 (in Chinese).
何颖,于琴,刘中宪. 2019. 考虑散射效应沉积河谷空间相关多点地震动模拟[J]. 岩土力学,40(7):2739–2747. He Y,Yu Q,Liu Z X. 2019. Simulation of multi-point spatially correlated earthquake ground motions of sedimentary valleys considering scattering effect[J]. Rock and Soil Mechanics,40(7):2739–2747 (in Chinese).
景立平,陈国兴,李永强,汤皓. 2009. 汶川8.0级地震水坝震害调查[J]. 地震工程与工程振动,29(1):14–23. Jing L P,Chen G X,Li Y Q,Tang H. 2009. Investigation on earthquake-induced dam damage during Wenchuan earthquake[J]. Journal of Earthquake Engineering and Engineering Vibration,29(1):14–23 (in Chinese).
廖振鹏,刘晶波. 1992. 波动有限元模拟的基本问题[J]. 中国科学:B辑,22(8):874–882. Liao Z P,Liu J B. 1992. Fundamental problems in finite element simulation of wave motion[J]. Science in China:Series B,35(11):1353–1364.
林皋,关飞. 1990. 用边界元法研究地震波在不规则地形处的散射问题[J]. 大连理工大学学报,30(2):145–152. Lin G,Guan F. 1990. Scattering of seismic waves at irregular topographies by boundary element method[J]. Journal of Dalian University of Technology,30(2):145–152 (in Chinese).
刘中宪,苗岳云,陈頔. 2020. 点震源作用下三维沉积盆地地震动谱元模拟[J]. 世界地震工程,36(2):200–208. Liu Z X,Miao Y Y,Chen D. 2020. Seismic spectral element simulation of three-dimensional sedimentary basin under the action of point seismic source[J]. World Earthquake Engineering,36(2):200–208 (in Chinese).
孙纬宇,汪精河,严松宏,欧尔峰,梁庆国. 2019. SV波斜入射下河谷地形地震动分布特征分析[J]. 振动与冲击,38(20):237–243. Sun W Y,Wang J H,Yan S H,Ou E F,Liang Q G. 2019. Characteristic analysis of ground motions of a canyon topography under obliquely incident SV waves[J]. Journal of Vibration and Shock,38(20):237–243 (in Chinese).
唐晖,李小军,李亚琦. 2012. 自贡西山公园山脊地形场地效应分析[J]. 振动与冲击,31(8):74–79. doi: 10.3969/j.issn.1000-3835.2012.08.015 Tang H,Li X J,Li Y Q. 2012. Site effect of topograghy on ground motions of Xishan park of Zigong city[J]. Journal of Vibration and Shock,31(8):74–79 (in Chinese).
王海云,谢礼立. 2010. 自贡市西山公园地形对地震动的影响[J]. 地球物理学报,53(7):1631–1638. doi: 10.3969/j.issn.0001-5733.2010.07.014 Wang H Y,Xie L L. 2010. Effects of topography on ground motion in the Xishan park,Zigong city[J]. Chinese Journal of Geophysics,53(7):1631–1638 (in Chinese).
王文才,尹志文,苏小芸,江志杰. 2020. 2018年陕西宁强5.3级地震强地面运动特征及局部场地效应分析[J]. 地震工程学报,42(6):1700–1705. doi: 10.3969/j.issn.1000-0844.2020.06.1700 Wang W C, Yin Z W, Su X Y, Jiang Z J. 2020. Characteristics of strong ground motion of the 2018 Ningqiang M5.3 earthquake in Shaanxi province and local site response analysis[J]. China Earthquake Engineering Journal, 42（6）: 1700–1705 （in Chinese）.
杨宇,李小军,贺秋梅. 2011. 自贡西山公园山脊场地地形和土层效应数值模拟[J]. 震灾防御技术,6(4):436–447. doi: 10.3969/j.issn.1673-5722.2011.04.009 Yang Y,Li X J,He Q M. 2011. Numerical simulation for site effect of ridge terrain and overlaying soil in Zigong Xishan park[J]. Technology for Earthquake Disaster Prevention,6(4):436–447 (in Chinese).
于彦彦,丁海平,刘启方. 2017. 透射边界与谱元法的结合及对波动模拟精度的改进[J]. 振动与冲击,36(2):13–22. 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[J]. Journal of Vibration and Shock,36(2):13–22 (in Chinese).
张宁,潘家琐,代登辉,高玉峰. 2021. 含峭壁V形峡谷对地震SH波散射的解析解[J]. 地球物理学报,64(3):896–906. doi: 10.6038/cjg2021O0041 Zhang N,Pan J S,Dai D H,Gao Y F. 2021. An analytical solution to the scattering of plane SH waves by a V-shaped canyon with cliffs[J]. Chinese Journal of Geophysics,64(3):896–906 (in Chinese).
章小龙,李小军,周正华,陈国兴,彭小波. 2017. 三维复杂山谷地形SV波垂直输入地震反应分析[J]. 地球物理学报,60(7):2779–2790. doi: 10.6038/cjg20170723 Zhang X L,Li X J,Zhou Z H,Chen G X,Peng X B. 2017. The seismic response analysis of three-dimensional Canyon complex topography under incident SV seismic waves[J]. Chinese Journal of Geophysics,60(7):2779–2790 (in Chinese).
Boore D. 1972. Finite difference methods for seismic wave propagation in heterogeneous materials[J]. Methods Comput Phys:Adv Res Appl,11:1–37.
Gao Y F,Zhang N,Li D Y,Liu H L,Cai Y Q,Wu Y X. 2012. Effects of topographic amplification induced by a U-shaped canyon on seismic waves[J]. Bull Seismol Soc Am,102(4):1748–1763. doi: 10.1785/0120110306
Gao Y F,Zhang N. 2013. Scattering of cylindrical SH waves induced by a symmetrical V-shaped canyon:Near-source topographic effects[J]. Geophys J Int,193(2):874–885. doi: 10.1093/gji/ggs119
Huang H C,Chiu H C. 1995. The effect of canyon topography on strong ground motion at Feitsui damsite:Quantitative results[J]. Earthq Eng Struct Dyn,24(7):977–990. doi: 10.1002/eqe.4290240705
Liu Z X,Wang D,Liang J W,Wu F J,Wu C Q. 2018. The fast multi-pole indirect BEM for solving high-frequency seismic wave scattering by three-dimensional superficial irregularities[J]. Eng Anal Boun Elem,90:86–99. doi: 10.1016/j.enganabound.2018.02.009
Parolai S,Bindi D,Baumbach M,Grosser H,Milkereit C,Karakisa S,Zünbül S. 2004. Comparison of different site response estimation techniques using aftershocks of the 1999 Izmit earthquake[J]. Bull Seismol Soc Am,94(3):1096–1108. doi: 10.1785/0120030086
Sanchez-Sesma F J. 1985. Diffraction of elastic SH-wave in wedges[J]. Bull Seismol Soc Am,75(5):1435–1446.
Spudich P,Hellweg M,Lee W H K. 1996. Directional topographic site response at Tarzana observed in aftershocks of the 1994 Northridge,California,earthquake:Implications for mainshock motions[J]. Bull Seismol Soc Am,86(1B):S193–S208. doi: 10.1785/BSSA08601BS193
Trifunac M D,Hudson D E. 1971. Analysis of the Pacoima dam accelerogram:San Fernando,California,earthquake of 1971[J]. Bull Seismol Soc Am,61(5):1393–1411.
Trifunac M D. 1973. Scattering of plane SH waves by a semi-cylindrical canyon[J]. Earthq Eng Struct Dyn,1(3):267–281.
Tsaur D H,Chang K H. 2008. An analytical approach for the scattering of SH waves by a symmetrical V-shaped canyon:Shallow case[J]. Geophys J Int,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[J]. J Eng Mech,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[J]. Geophys J Int,183(3):1501–1511. doi: 10.1111/j.1365-246X.2010.04806.x
Yuan X M,Liao Z P. 1994. Scattering of plane SH waves by a cylindrical canyon of circular-arc cross-section[J]. Soil Dyn Earthq Eng,13(6):407–412. doi: 10.1016/0267-7261(94)90011-6
Zhang N,Gao Y F,Li D Y,Wu Y X,Zhang F. 2012a. Scattering of SH waves induced by a symmetrical V-shaped canyon:A unified analytical solution[J]. Earthq Eng Eng Vib,11(4):445–460. doi: 10.1007/s11803-012-0135-z
Zhang N,Gao Y F,Cai Y Q,Li D Y,Wu Y X. 2012b. Scattering of SH waves induced by a non-symmetrical V-shaped canyon[J]. Geophys J Int,191(1):243–256. doi: 10.1111/j.1365-246X.2012.05604.x
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[J]. Soil Dyn Earthq Eng,95:17–28. doi: 10.1016/j.soildyn.2017.01.037
Zhang N,Zhang Y,Gao Y F,Pak R Y S,Wu Y X,Zhang F. 2019. An exact solution for SH-wave scattering by a radially multilayered inhomogeneous semicylindrical canyon[J]. Geophys J Int,217(2):1232–1260. doi: 10.1093/gji/ggz083
-
期刊类型引用(7)
1. 刘中宪，孟思博，张妤，乔云帆，陈龙伟. 考虑建筑群-沉积盆地动力相互作用的建筑群震害评估方法. 地震学报. 2024(01): 129-143 . 
本站查看
2. 王睿，王兰民，周燕国，王刚. 土动力学与岩土地震工程. 土木工程学报. 2024(07): 71-89+105 . 百度学术
3. 赵仕兴，罗麒锐，熊峰，吴启红，夏静，杨姝姮，张敏. 复杂区域地形条件下台地地震动放大效应研究. 振动与冲击. 2024(23): 230-239+304 . 百度学术
4. 卢林彤，李孝雄，庄海洋，吴琪. 宽河谷非均匀工程场地非线性地震反应空间变异特征研究. 地震工程与工程振动. 2023(04): 216-225 . 百度学术
5. 柳国环，陈欣宇，李鑫洋. 半覆水相变V形峡谷场地P波入射下的空间地震动散射频域理论解、验证及特征. 计算力学学报. 2023(05): 764-772 . 百度学术
6. 刘中宪，周健，李程程，孟思博. 基于修正随机有限断层法的沉积盆地地震动模拟. 防灾减灾工程学报. 2023(05): 999-1008 . 百度学术
7. 蒋乐英，廖意辉，王志明，吴勇信，李梅. 河谷差异地震作用下拱式倒虹吸地震响应分析. 长江科学院院报. 2022(12): 117-121+140 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 444
- HTML全文浏览量: 177
- PDF下载量: 71
- 被引次数: 7
