山脊线与坡度和峰值速度放大系数的相关性研究.docxVIP

? ? 山脊线与坡度和峰值速度放大系数的相关性研究 ? ? 蒋涵，周红，高孟潭 中国地震局地球物理研究所，北京 100081 1 引言 人们很早就认识到地表的起伏能够对地震动产生很大的影响.通常山顶或者山脊能造成地震动的振幅增加，而山谷则会造成地震动的振幅减小.地震灾害调查也发现在山顶上的建筑往往会遭受到比周围位于平地或山谷中的建筑更严重的损害.不仅如此，山体的顶部在地震发生后也更容易发生崩塌和滑坡（黄润秋和李为乐，2008；王运生等，2009）.因此起伏的地形对地震动的影响的研究，在地震学和工程学上都具有重要意义. 通过震害调查和数值模拟的研究，学者发现地形效应存在着一定规律.比如当入射的波长和起伏的几何特征接近的时候，地形效应会非常明显.而在入射波的波长较大的情况下，地形效应则可以被忽略（Ashford et al.，1997；Geli et al.，1988）.研究也发现，许多因素如入射波场，地下介质，以及局部地形几何形状等都会影响地形效应.Bouchon（1973）分别研究了P，SH和SV波垂直入射到二维峡谷形态的地形效应.他指出在山谷两侧的顶部，入射的SV波所产生的放大倍数最高.Bard（1982）设定了一个正弦形态的山脊，并用SH作为入射波.他发现向前传播的散射波振幅会随入射角度的增加而增大.Geli等（1988）总结了前人对二维的三角与弧形地形的研究成果，他从中得到了随山体高度与基底宽度的比值增加，对应的山顶和基底的地震动幅值之比也会增加的结论，作者还就介质起伏对地形效应的影响进行了探讨.周红等（2010）针对地形拐角角度变化对放大的影响展开研究，结果发现当台阶地形的张角小于180°时，地震动放大会随着角度的减小而增加. 以上的研究成果采用的是简单的二维模型，其计算结果与三维模型的结果相比存在着不小的差异.在三维地形效应中，地表的放大倍数更高，放大模式也更复杂.Liao等（1981）采用了三维截圆柱的模型，模拟得到的放大倍数是相应二维模型的两倍.Spudich等（1996）通过分析Northridge earthquake的强地震动记录，观测到垂直于山脊轴线方向的强烈的震动，其强度大大超过了二维模型的计算值，由此作者推测这是由三维的地形所导致的.该推测被随后三维模型的数值模拟（Bouchon et al.，1996）所证实. 上述几个三维数值模拟所研究的是孤立的山体，真实的地形则往往是多个山脊相连，因此局部的地形效应除了受山脊本身的影响以外，还会受到周围山脊的影响，比如不同走向的山脊的交汇处会有所谓的角落效应（Wang et al.，2006）.对这种真实条件下的三维地形已开展不少数值研究，比如San Andreas Fault地区（Ma et al.，2007），汶川地区（Zhang et al.，2008），台北地区（Lee et al.，2009a，2009b）的研究.学者通过将含真实起伏的地形和水平地形的模拟结果进行比较，从而得到该地区的地震动放大系数的分布特征.然而和简单地形模型相比，真实地貌特征和放大系数的关系还缺乏深入讨论.这当中有个重要原因就是三维地形特征的复杂性. 本文围绕着真实地貌特征与放大系数的特征展开研究工作.为了更好地分析研究地形特征，我们引入地理中的地形分析法，将山脊线和坡度作为地形特征因子来描述山体的特征.这样就将对三维地形特征和放大系数的相关性的研究，转换为对山脊线与坡度和放大系数的相关性的研究.不仅让分析的难度得以降低，也使得分析能够定量化.基于这个思路，本文选取了地形起伏剧烈的芦山地区作为研究对象，用谱元法模拟了由爆炸源产生的地形放大系数.通过分析地形因子和地震动放大系数分布特征，得到了它们之间存在的相关性. 2 芦山地区的三维地形效应模拟 2.1 数值方法 数值模拟是研究地形效应的有力工具，本文中所采用的数值方法为谱元法（SEM）.它最早应用于流体力学的计算中（Patera，1984）.后来由Komatitsch和 Vilotte（1998），Komatitsch和 Tromp（1999）引入到地震波场模拟计算中.谱元法结合了有限元的灵活划分网格特点和谱方法高精度的优势，因而在计算地震学中被广泛使用. 对于在地下介质中传播的地震波，其位移场满足的力学方程如下： 其中ρ是密度，u是位移向量，σ是应力张量，f则表示外力. 谱元法所采用的是（1）式的积分形式.我们将（1）乘以一个权函数w，并在整个计算区域进行分步积分，从而得到如下的积分方程：其中Ω是物理计算区域，Γ是区域的边界，t为边界上的力. 谱元法把计算区域离散为一个个单元.每个单元 由 Gauss－Lobatto－Ledendre（GLL）点 构 成 的Lagrange插值多项式表示.有限元中通常使用低阶次的单元，而