活动断裂带的整体地震复发行为及其与分段地震复发行为的关系
-
摘要: 利用中国大陆11条活动断裂带的地震资料,研究了板内活动断裂带整体的地震复发行为及其与单段的地震复发行为的关系.结果表明:由多个段落组成的活动断裂带的整体地震复发表现出成丛、完全随机和差的准周期3种行为,但以稀疏成丛行为为主,不同程度地表现出活跃期与平静期的交替;丛与丛之间以及丛内地震的时间间隔、丛的持续时间、地震次数和强度均是可变的;地震复发过程是非线性的,不存在强度时间、时间强度的相依性.然而,单个断裂段的地震复发过程要简单得多,并主要表现出准周期和时间可预报行为.本文还进一步揭示:断裂带整体地震复发过程中的丛主要是由不同断裂段在较短时期内的破裂传染作用造成的;沿断裂带破裂传染作用的程度、方向性可因不同轮回而变化,且传染似可跳过断裂带上不破裂的空区.
-
关键词:
- 活动断裂 地震复发行为 成丛 破裂传染
-
引言
我国江河密布,许多重要建筑、桥梁、大坝等坐落在河谷附近。多次震害调查表明,河谷场地对地震动具有显著放大效应。国内外很多专家针对不含流体河谷对地震波的散射问题进行了研究(周国良等,2012;陈少林等,2014;Liu et al,2016;梁建文等,2017;刘中宪等,2017;张宁等,2017;Zhang et al,2017),然而河谷中通常包含一定深度的流体,地震波在流体/固体交界面上会发生复杂的透射与反射现象,从而影响水下地层的地震反应规律,研究流体层对河谷场地地震动的影响意义重大。
王进廷等(2003,2004)研究了P波、SV波入射时弹性半空间上理想流体层的动力响应,随后给出了流体-固体-多孔介质耦合场地中动压力的解析解(Wang et al,2004,2009);李伟华(2010)采用有限元方法研究了考虑水-饱和土场地-结构耦合时的沉管隧道地震反应;Carbajal-Romero等(2013)采用间接边界元法(indirect boundary element method,缩写为IBEM)研究了Scholte波在流体-固体界面处的扩散;Alejandro等(2014)采用边界元法分析了理论地震事件发生时海洋水体的动态响应;杜修力等(2015)根据地震与波浪作用的动水压力解析解,研究了SV波斜入射时地震和波浪联合作用下自由场海水的动水压力反应问题;张奎等(2018)采用解析解方法得到了水下地层表面位移的表达式,并分析了含水层场地模型中土的刚度、饱和度及平面P波、SV波的入射角等因素变化时流体深度对场地位移响应的影响。
需注意的是以上研究均是针对含水平成层流体的场地动态响应问题,而考虑地震波散射效应对含流体层河谷地形的地震动影响的研究仍非常有限。李伟华和赵成刚(2006)利用波函数展开法,首次在频率域内给出了具有饱和土沉积层的圆弧形充水河谷对平面P波、SV波散射问题的解析解,继而给出了瑞雷波的解析解(赵成刚等,2008)。解析方法中波场构造需严格满足波动方程和边界条件,对于求解复杂的偏微分算子及边界条件难度较大。
在各类数值方法中,边界元法具有降低问题求解维数、自动满足无限远辐射条件且无高频数值弥散的优点,因而在地震波动研究中获得广泛应用。本文尝试将间接边界元法(IBEM)拓展到含流体层局部场地对地震波的散射求解:结合流体域格林函数,将河谷中的流体层模拟为无黏性、可压缩流体,考虑固体/液体交界面上地震波的透射与反射问题,研究了含流体层河谷场地在平面P波、SV波入射时的地震响应。通过与现有文献结果作对比,验证了该方法的精度和数值稳定性。进而结合具体算例,给出了入射波的频率、角度、流体深度等因素对河谷底部及附近地表位移放大系数的影响,为实际情况中含水河谷地震动的确定提供了定性和定量的参考依据。
1. 计算模型
图1给出了计算模型。假设半空间河谷中为各向同性弹性介质,Ω1表示半空间域,Ω2表示水域;B表示水域与半空间域的交界面;与空气直接接触的表面用L表示,L1表示两侧水平地表,L3表示流体表面;θ为P波或SV波入射方向与竖直向的夹角;S1表示半空间域的虚拟荷载面,S2表示水域的虚拟荷载面。假设P波、SV波从外部半空间入射,待求问题即为含流体层河谷场地对地震波的散射。
2. 间接边界元法的求解
2.1 半空间波场模拟
半空间Ω1内的总位移场和总应力场可表示为
$u_i^{}(x{\text{,}}\omega) {\text{=}} u_i^{\rm{s}}(x{\text{,}}\omega) {\text{+}} u_i^{\rm{f}}(x{\text{,}}\omega){\text{,}}x \in {\varOmega _1}{\text{,}}$
(1) $t_i^{}(x{\text{,}}\omega) {\text{=}} t_i^{\rm{s}}(x{\text{,}}\omega) {\text{+}} t_i^{\rm{f}}(x{\text{,}} \omega){\text{,}}x \in {\varOmega _1}{\text{,}}$
(2) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
$u_i^{\rm{s}}(x{\text{,}}\omega) {\text{=}} \int_{{S_1}} {\phi _j^{\rm{s}}(\xi{\text{,}}\omega)G_{ij}^{\rm{s}}(x{\text{,}}\xi){\rm{d}}{S_\xi }}{\text{,}}$
(3) $t_i^{\rm{s}}(x{\text{,}}\omega) {\text{=}} {c_2}\phi _j^{\rm{s}}(x{\text{,}}\omega) {\text{+}} \int_{{S_1}} {\phi _j^{\rm{s}}(\xi{\text{,}}\omega)T_{ij}^{\rm{s}}(x{\text{,}}\xi){\rm{d}}{S_\xi }}{\text{,}}$
(4) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
2.2 流体层中的波场模拟
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
${w_r}(x{\text{,}}\omega) {\text{=}} {c_1}\psi (x{\text{,}}\omega) {\text{+}} \frac{1}{{\rho {\omega ^2}}}\int_{{S_2}} \psi (\xi{\text{,}}\omega)\frac{{\partial {G^{\rm{w}}}(x{\text{,}}\xi)}}{{\partial r}}{\rm{d}}{S_\xi }{\text{,}}$
(5) $p(x{\text{,}}\omega){\text{=}} \int_{{S_2}} \psi (\xi{\text{,}}\omega){G^{\rm{w}}}(x{\text{,}}\xi){\rm{d}}{S_\xi }{\text{,}}$
(6) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
${G}^{\rm{w}}(x{\text{,}}\xi){\text{=}}\frac{\rho {\omega }^{2}}{4\rm{i}} \cdot {{\rm{H}}}_{0}^{(2)} \cdot \frac{\omega d}{{c}^{\rm{w}}},$
(7) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
2.3 边界条件及求解
边界条件主要分为三部分:一部分是直接与空气接触的地表L1,满足地表零应力条件;一部分是与空气接触的水体表面L3,满足流体层表面孔压为零;另一部分是半空间与流体域的交界面B处,满足位移连续条件和应力连续条件,包括流量位移和孔压连续及固体位移和应力连续。可分别表示为
$t_\tau ^{\rm{s}}(x{\text{,}}\omega) {\text{+}} t_\tau ^{\rm{f}}(x{\text{,}}\omega) {\text{=}} 0{\text{,}}x \in {L_1}{\text{,}}$
(8) $t_r^{\rm{s}}(x{\text{,}}\omega) {\text{+}} t_r^{\rm{f}}(x{\text{,}}\omega) {\text{=}} 0{\text{,}}x \in {L_1}{\text{,}}$
(9) $p(x{\text{,}}\omega) {\text{=}} 0{\text{,}}x \in {L_3}{\text{,}}$
(10) $t_\tau ^{\rm{s}}(x{\text{,}}\omega) {\text{+}} t_{\rm{\tau }}^{\rm{f}}(x{\text{,}}\omega) {\text{=}} 0{\text{,}}x \in B{\text{,}}$
(11) $t_r^{\rm{s}}(x{\text{,}}\omega) {\text{+}} t_r^{\rm{f}}(x{\text{,}}\omega) {\text{=-}} p{\text{,}}x \in B{\text{,}}$
(12) $u_r^{\rm{s}}(x{\text{,}}\omega) {\text{+}} u_r^{\rm{f}}(x{\text{,}}\omega) {\text{=}} {w_r}(x{\text{,}}\omega){\text{,}}x \in B{\text{,}}$
(13) ${c_2}\phi _j^{\rm{s}}(x{\text{,}}\omega) {\text{+}} \int_{{S_1}} {\phi _j^{\rm{s}}(\xi{\text{,}}\omega)T_{\tau j}^{\rm{s}}(x{\text{,}}\xi){\rm{d}}{S_\xi }} {\text{=-}}t_\tau ^{\rm{f}}(x{\text{,}}\omega){\text{,}}x \in {L_1}{\text{,}}$
(14) ${c_2}\phi _j^{\rm{s}}(x{\text{,}}\omega) {\text{+}} \int_{{S_1}} {\phi _j^{\rm{s}}(\xi{\text{,}}\omega)T_{rj}^{\rm{s}}(x{\text{,}}\xi){\rm{d}}{S_\xi }} {\text{=-}} t_r^{\rm{f}}(x{\text{,}}\omega){\text{,}}x \in {L_1}{\text{,}}$
(15) $\int_{{S_2}} \psi (\xi{\text{,}}\omega){G^{\rm{f}}}(x{\text{,}}\xi){\rm{d}}{S_\xi } {\text{=}} 0{\text{,}}x \in {L_3}{\text{,}}$
(16) ${c_2}\phi _j^{\rm{s}}(x{\text{,}}\omega) {\text{+}} \int_{{S_1}} {\phi _j^{\rm{s}}(\xi{\text{,}}\omega)T_{\tau j}^{\rm{s}}(x{\text{,}}\xi){\rm{d}}{S_\xi }} {\text{=-}} t_\tau ^{\rm{f}}(x{\text{,}}\omega){\text{,}}x \in B{\text{,}}$
(17) ${c_2}\phi _j^{\rm{s}}(x{\text{,}}\omega) {\text{+}} \int_{{S_1}} {\phi _j^{\rm{s}}(\xi{\text{,}}\omega)T_{rj}^{\rm{s}}(x{\text{,}}\xi){\rm{d}}{S_\xi }} {\text{+}} \int_{{S_2}} \psi (\xi{\text{,}}\omega){G^{\rm{f}}}(x{\text{,}}\xi){\rm{d}}{S_\xi } {\text{=-}} t_r^{\rm{f}}(x{\text{,}}\omega){\text{,}}x \in B{\text{,}}$
(18) $\int_{{S_1}} {\phi _j^{\rm{s}}(\xi{\text{,}}\omega)G_{rj}^{\rm{s}}(x{\text{,}}\xi){\rm{d}}{S_\xi }} - {c_1}\psi (x{\text{,}}\omega) - \frac{1}{{\rho {\omega ^2}}}\int_{{S_2}} \psi (\xi{\text{,}}\omega)\frac{{\partial {G^{\rm{f}}}(x{\text{,}}\xi)}}{{\partial r}}{\rm{d}}{S_\xi } {\text{=-}} u_r^{\rm{f}}(x{\text{,}}\omega){\text{,}}x \in B{\text{,}}$
(19) 存在线性方程组
$\sum\limits_{l {\text{=}} 1}^N {\left[ {\sum\limits_{j {\text{=}} 1}^2 {\phi _j^{\rm{s}}({\xi _l}{\text{,}}\omega)t_{\tau j}^{\rm{s}}({x_n}{\text{,}}{\xi _l})} } \right]} {\text{=-}} t_\tau ^{\rm{f}}({x_n}{\text{,}}\omega){\text{,}}{x_n} \in {L_1}{\text{,}}$
(20) $\sum\limits_{l {\text{=}} 1}^N {\left[ {\sum\limits_{j {\text{=}} 1}^2 {\phi _j^{\rm{s}}({\xi _l}{\text{,}}\omega)t_{rj}^{\rm{s}}(x{\text{,}}{\xi _l})} } \right]} {\text{=-}} t_r^{\rm{f}}({x_n}{\text{,}}\omega){\text{,}}{x_n} \in {L_1}{\text{,}}$
(21) $\sum\limits_{l {\text{=}} 1}^N {\left[ {\sum\limits_{j {\text{=}} 1}^2 {\psi ({\xi _l}{\text{,}}\omega){g^{\rm{f}}}({x_n}{\text{,}}{\xi _l})} } \right]} {\text{=}} 0{\text{,}}{x_n} \in {L_3}{\text{,}}$
(22) $\sum\limits_{l {\text{=}} 1}^N {\left[ {\sum\limits_{j {\text{=}} 1}^2 {\phi _j^{\rm{s}}({\xi _l}{\text{,}}\omega)t_{\tau j}^{\rm{s}}({x_n}{\text{,}}{\xi _l})} } \right]} {\text{=-}} t_\tau ^{\rm{f}}({x_n}{\text{,}}\omega){\text{,}}{x_n} \in B{\text{,}}$
(23) $\sum\limits_{l {\text{=}} 1}^N {\left[ {\sum\limits_{j {\text{=}} 1}^2 {\phi _j^{\rm{s}}({\xi _l}{\text{,}}\omega)t_{rj}^{\rm{s}}({x_n}{\text{,}}{\xi _l})} } \right]} {\text{+}} \sum\limits_{l {\text{=}} 1}^N {\left[ {\sum\limits_{j {\text{=}} 1}^2 {\psi ({\xi _l}{\text{,}}\omega)g_1^{\rm{f}}({x_n}{\text{,}}{\xi _l})} } \right]} {\text{=-}} t_r^{\rm{f}}({x_n}{\text{,}}\omega){\text{,}}{x_n} \in B{\text{,}}$
(24) $\sum\limits_{l {\text{=}}1}^N {\left[ {\sum\limits_{j {\text{=}} 1}^2 {\phi _j^{\rm{s}}({\xi _l}{\text{,}}\omega)g_{rj}^{\rm{s}}({x_n}{\text{,}}{\xi _l})} } \right]} {\text{-}} \sum\limits_{l {\text{=}} 1}^N {\left[ {\sum\limits_{j {\text{=}}1}^2 {\psi ({\xi _l}{\text{,}}\omega)g_2^{\rm{f}}({x_n}{\text{,}}{\xi _l})} } \right]} {\text{=-}} u_r^{\rm{f}}({x_n}{\text{,}}\omega){\text{,}}{x_n} \in B{\text{,}}$
(25) 式中,
$t_{\tau j}^{\rm{s}}({x_n}{\text{,}}{\xi _l}){\text{=}}{c_{\rm{2}}}{\delta _{ij}}{\delta _{nl}} {\text{+}} \int_{{S_1}} {T_{\tau j}^{\rm{s}}({x_n}{\text{,}}{\xi _l}){\rm{d}}{S_\xi }}{\text{,}}$
(26) $t_{rj}^{\rm{s}}({x_n}{\text{,}}{\xi _l}){\text{=}}{c_{\rm{2}}}{\delta _{ij}}{\delta _{nl}} {\text{+}} \int_{{S_1}} {T_{rj}^{\rm{s}}({x_n}{\text{,}}{\xi _l}){\rm{d}}{S_\xi }}{\text{,}}$
(27) $g_1^{\rm{f}}({x_n}{\text{,}}{\xi _l}){\text{=}}\int_{{S_2}} {{G^{\rm{f}}}({x_n}{\text{,}}{\xi _l}){\rm{d}}{S_\xi }}{\text{,}}$
(28) $g_{rj}^{\rm{s}}({x_n}{\text{,}}{\xi _l}) {\text{=}} \int_{{S_1}} {G_{rj}^{\rm{s}}({x_n}{\text{,}}{\xi _l}){\rm{d}}{S_\xi }}{\text{,}}$
(29) $g_2^{\rm{f}}({x_n}{\text{,}}{\xi _l}){\text{=}}{c_{\rm{1}}}{\delta _{ij}}{\delta _{nl}} {\text{+}} \frac{1}{{\rho {\omega ^2}}}\int_{{S_2}} {\frac{{\partial {G^{\rm{f}}}({x_n}{\text{,}}{\xi _l})}}{{\partial r}}{\rm{d}}{S_\xi }}{\text{,}}$
(30) n1=1,2,···,N1;n2=1,2,···,N2;N为S1,S2上点的个数。
3. 精度验证
为验证本文IBEM方法的正确性,令流体深度h=0,计算模型如图2所示,则本文模型退化为不含流体层的河谷地形,并将所得结果与既有文献(Sánchez-Sesma,Campillo,1991)结果进行对比(图3)。计算参数为:入射波为平面SV波,弹性介质中剪切波速cβ1为1 000 m/s,压缩波速cα1为2 000 m/s,密度ρ1为2 400 kg/m3;峡谷中空气的压缩波速cα1为330 m/s,密度ρ2为1.29 kg/m3;无量纲入射波频率η=ωr/πcβ=2.0,a为河谷半径,入射波角度分别取θβ=0°,θβ=30;阻尼比ζ=0.001,泊松比v=1/3。图3表明本文结果与文献结果吻合良好,验证了计算方法的正确性。
图 3 平面SV波入射下本文含水河谷退化位移与文献结果对比 (引自Sánchez-Sesma,Campillo,1991)Figure 3. Comparison between the degradation displacement of water bearing valley of this paper and the results of literature induced by plane SV waves (after Sánchez-Sesma,Campillo,1991)(a) η=2.0,θβ=0°;(b) η=2.0,θβ=30°4. 算例分析
采用上述间接边界元方法,研究含流体层的河谷地形对平面P波、SV波的地震响应,计算模型为含流体层的半圆形河谷地形(图2)。阻尼比ζ=0.001,计算参数列于表1。
表 1 含流体层河谷地形对平面P波、SV波的地震响应计算参数Table 1. Calculation parameters of seismic response of valley terrain with fluid layer for plane P wave and SV waveP波波速cα/(m·s−1) SV波波速cβ/(m·s−1) 密度ρ/(kg·m−3) 空气 330 — 1.29 流体 1 501 — 1 000 弹性土体 2 670 1 090 2 200 4.1 充满流体的河谷对平面P波、SV波的散射作用
对于含流体层的河谷地形而言,流体层的存在会对峡谷地形的地震波散射产生很大的影响,影响因素主要包括:入射波的类型、角度、频率以及流体深度。图4,5给出了在不同影响因素作用下,平面P波、SV波入射时充满流体层(h/r=1)的河谷地形和不含流体的河谷地形表面主方向位移放大系数对比。入射波角度取θα=0°,30°;θβ=0°,30°。
4.1.1 流体对河谷附近地震动响应单频分析
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
总体来看,与不含流体层的峡谷地形相比,充满流体的河谷场地对地表地震动存在不同程度的缩幅效应。例如P波入射θα=0°,η=0.5时,在x/a=2处,不含流体的峡谷地形的竖向位移放大系数为2.66,含流体层的河谷地形的竖向位移放大系数仅为2.03,降低了约24% (图4);SV波入射θβ=0°,η=0.5时,在x/a=2处,不含流体的峡谷地形的竖向位移放大系数为2.74,含流体层的河谷地形的竖向位移放大系数仅为2.02,降低了约26% (图5)。
由图4—5可以看出,由于流体可吸收地震波中的纵波能量,使地震波的散射作用减弱,地表位移震荡较平缓。考虑入射角度的影响,与垂直入射相比,斜入射时,河谷左侧(靠近波入射方向)位移比右侧(远离波入射方向)位移震荡剧烈且幅值较大。
由图4—5还可以看出,与斜入射(θ=30°)相比,垂直入射(θ=0°)流体对河谷地形地震动响应的影响较大,位移放大系数的放大与减小效应更加明显。
4.1.2 流体对河谷地形地震动响应频谱分析
图6为充满流体层的河谷底部与流体表面位移放大系数谱对比,图7,8为不含流体层的河谷场地和含流体层的河谷场地地表典型观察点处的位移放大系数谱,图中选取了x/a=0,0.5 (河谷底部)和1.5 (右侧水平地表)三个典型点位作为观察点。入射波无量纲频率η=0—5,入射角度取θα=0°,30°;θβ=0°,30°。
P波入射时,在低频域内(η<0.5),含流体河谷底部及附近地表的地震动反应与不含流体的河谷反应基本一致。随着入射频率的增大,河谷底部观察点的位移反应特性发生改变,在共振频率处河谷底部位移缩小效应显著,可缩减至接近于0 (图7)。例如,P波垂直入射时,在河谷-流体体系一阶(η=0.75)、二阶(η=2.2)和三阶(η=3.6)共振频率处,河谷底部(x/a=0)竖向位移分别为0.32,0.13,0.34。P波斜入射和其它点位也表现出这一特性。观察图6可以看出,河谷底部位移最小频率处刚好流体表面位移达到最大,这是因为在该频率下大量地震波动能量进入流体层内,对下部河谷波动反而产生了抑制作用。因此在某些频段,流体层的存在对其下部地层具有明显的减震效应。
SV波入射时,随着入射频率的增大,水层体系共振频率增多,在共振频率处,位移频谱曲线震荡更为剧烈。在河谷外部附近地表观察点处地震动反应主要表现为位移放大系数的减小。
4.2 含流体层河谷地形中流体深度对平面P波、SV波散射的影响
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
流体深度对P波入射时地震动的影响较大,对SV波入射时地震动的影响较小。P波入射时,流体层的深度对河谷附近地表的位移放大系数影响显著,而对河谷内部位移放大系数影响较小。
5. 结论
本文采用间接边界元方法,求解了含流体层河谷地形对平面P波、SV波的散射,定性和定量地讨论了地震波入射时流体层对河谷地震放大效应的影响规律和地震动空间分布特征。主要结论如下:
1) 流体层具有吸收地震波能量的作用。整体上看,流体深度越大,河谷表面及附近地表的地震动幅值越小,流体深度对P波入射动力反应比对SV波入射情况的影响要更为显著;
2) 当P波入射频率与河谷中流体层的体系共振频率接近时,大量地震能量进入流体层,流体层表面位移放大最为显著,但河谷由于受到流体层的抑制作用,其底部位移放大系数明显降低;
3) 与斜入射(θα=30°)相比,P波垂直入射(θα=0°)时,流体对河谷地形地震动响应的影响较大,位移放大系数的放大或减小效应更加明显。
含流体层河谷地形附近地震动空间变化十分显著,因此在河谷中修建水库、桥梁等重要工程时,有必要考虑流体-河谷局部场地效应以更科学地确定地震动参数。
-
计量
- 文章访问数: 1522
- HTML全文浏览量: 17
- PDF下载量: 43