逆冲断层错动形成变形场及应力场的有限元模型——以映秀-青川断层为例.pdfVIP

第 35 卷第4 期 Vol. 35 , No. 4 地震地质 2013 年 12 月 SEISMOLOGY AND GEOLOGY Dec. , 2013 doi: 10. 3969/j. issn. 0253-4967.2013.04.003 逆冲断层错动形成变形场及 应力场的有限元模型 -一以映秀-青川断层为例 白玉柱周庆 (中国地震局地质研究所，活动构造与火山重点实验室，北京 100029) 摘 要 为了研究逆断层受到垂直于断层走向的应力作用下发生逆冲构造运动时造成的变形以及 相应的应力场分布，以龙门山中央断层映秀一青川段为例建立了逆冲断层的有限元模型。该模型为 单个逆冲断层有限元模型，计算了出川地震过程中龙门山中央断层边冲错动引起的地表变形场和断 层附近应力场的变化。模拟中对断层面采用的是摩擦接触单元，且断层面设为曲丽的几何形状:模 拟给出断层错动形成的地表沿垂直走向和垂直地袤的变形以及相应的应力分量场的变化情况。为得 到计算模型边界处的应力值，采用了断层静态位错模型进行计算。结果发现，在不施加沿断层走向的 载荷时，沿断层走向上存在变形，这可能是断层在受到挤压后发生的横向变形，但该变形较小;逆冲 断层构造运动引起沿垂直断层走向方向上的变形并不是在近断层处达到最大，而是沿垂直地表向上 的变形在断层处达到最大;逆冲断层在挤压构造应力下的变形要大于逆冲构造附近的隆升变形，同 时断层面上的挤压应力大于垂直地表方向的应力;断层面上的应力场沿断层面的变化主要集中在断 层面倾角发生变化以及震源区域附近。 关键词 龙门山断层地表变形有限元位错模型 文章编号:0253-4967(2013)04-0721-10 申图分类号: P315.2 文献标识码 :A 。 引言 地球表面附近深度范围内存在着各式各样的活动构造，其附近的变形场或应力场的产生、 改变以及许多地震地质现象都惊于活动构造的运动。例如，逆冲断层的活动将导致上盘的隆升 以及相应的地壳缩短c 因此，通过研究构造运动引起的变形场和应力场分布不但可以揭示许多 地震地质现象的本质特征，而且对大型工程建筑的施工及选址具有重要指导意义。 当前模拟构造运动引起变形的数值方法主要有静态位错模型、有限元模型、有限差分和边 界元模型等。采用静态位错模型能够有效地计算多个构造运动引起的变形问题，如白玉柱 (2009 ， 2010 ， 2011 ， 2012) 采用位错模型计算了龙门山断层错动引起变形的应变场和位移场的空 间分布;申重阳等(2008) 采用位错模型计算了汶川地震的同震重力场及其地表变形;张永志 〔收稿日期) 2013-04-01 收稿，2013-06-15 改回。 〔基金项目〕 国家科技支撑计划项目特大地震危险区识别及危险性评价方法研究(2012BAK15BOO) 和地 震科技人才骨干项t=I共同资助。 35 卷 722 地震地质 等(2009) 采用位错模型研究了走滑断层破裂过程中断层面上位移分布不均性对地表变形的影 响等。但是，静态位错模型只能计算几何形状较简单的情况，对复杂的断层面几何情况，静态 位错模型元能为力。此外，位错模型在计算构造变形或构造应力场时，并不考虑断层面j二的摩 擦性质，而断层面摩擦性质却是地震复杂成因的→个重要因素，不容忽略。有限元方法在处理 复杂断层几何结掏以及多变的地质情况方面有较强的优势，通过在有限元模型中引入特殊的单 元后，还可以考虑断层面的摩擦因素。目前，采用有限元进行断层数值模拟研究很多，如李红等 (2008) 应用三维有限元方