汶川地震断裂带附近建筑物的灾害效应模拟研究

苏生瑞，张永双，郝莉莉，李 松

(1.长安大学地质工程系，陕西 西安 710054;2.中国地质科学院地质力学研究所，北京 100081)

汶川地震断裂带附近建筑物的灾害效应模拟研究

苏生瑞1，张永双2，郝莉莉1，李 松1

(1.长安大学地质工程系，陕西 西安 710054;2.中国地质科学院地质力学研究所，北京 100081)

位于汶川地震破裂带附近的汶川水磨镇硅业公司厂房在地震中产生了独特的破坏现象。根据现场调查和测绘资料，建立了二维有限元模型，进行了断裂带附近建筑物灾害效应的数值模拟。模拟结果很好地解释了厂房在地震中产生的独特破坏现象的机理。研究表明，断层的形成改变了局部的应力场条件和加速度特征，从而使位于断层上盘的楼房发生差异性破坏，紧邻断层一侧破坏强烈。地震时楼体主梁的破坏是轴向应力和剪应力同时作用造成的。

断裂带;建筑物;灾害效应;数值模拟

0 引言

近断层大地震 (如:1994年美国North Ridge地震、1995年日本Kobe地震、1999年中国台湾集集地震、1999年土耳其Kocaeli地震、2003年伊朗Bam地震)独特的运动特征及其对工程结构的严重影响引起了地震工程界的密切关注［1～6］。Somerville等的研究表明，由于近断层地震动经常包含强烈的动态长周期脉冲和永久地面位移，其运动特征与远场地震动明显不同［1～2］。

汶川“5.12”震后调查表明:地震断裂带附近大量建筑物发生了破坏［7～13］，断层上下盘的建筑物在地震时的破坏程度不同，而且即使位于断层同一侧，由于与断层位置关系的差异，地表建筑物的破坏情况也不尽相同。因此，研究断层在地震作用下对地表建筑物的灾害效应，揭示地表破裂的致灾机理对抗震设计和灾害预防具有重要的理论意义和实用价值。

1 地震引发的建筑物差异性破坏

汶川“5.12”地震在汶川县水磨镇西侧斜坡上形成了产状为331°∠54°的地表破裂 (断层)，同时使位于断层上盘的硅业公司主楼 (距离断层约25 m)发生了差异性破坏现象，近断层一侧和远离断层一侧的破坏程度明显不同。

水磨镇硅业公司主楼走向北东—南西向，长35 m，宽12 m，高12 m，楼梯位于中部东南侧，厂房平面图见图1，其受地震的破坏表现为西北侧的梁、柱破坏较轻，仅一个柱 (图1中18号柱)发生向东南方向错动，最大错动距离约3 cm(见图2);而东南侧的梁、柱破坏严重，其中4个钢筋混凝土柱 (图1中9、10、13和14号柱)受到垂向挤压破坏并发生水平剪切，钢筋弯曲 (见图3、图4)。

图2 位于角部的柱体错动Fig.2 Dislocation of the pile

图3 靠近断层一侧钢筋混凝土柱破坏Fig.3 Failure of concrete pile adjacent to the fault

图4 钢筋混凝土柱遭挤压、剪切破坏(图1中14号柱)Fig.4 Sheared and compressed failure of the pile

本文基于对位于断层上盘的汶川县水磨镇硅业公司的实地调查，采用数值模拟方法研究了水磨镇硅业公司主楼在地震动作用下的动力响应特点，并分析了楼体破坏与断裂的关系。

2 研究区地质概况

汶川县水磨镇硅业公司位于龙门山构造带中南段的北端、水磨河河谷西北岸的高阶地上，其西北侧紧邻一个高度约35 m斜坡的坡脚。场地出露的地层有:泥盆系养马坝组(D2y)泥灰岩，厚度为89～137 m;第四系冲洪积物，覆盖于养马坝组之上，厚度为0～7.3 m，由砂、砾石组成，其与基岩的接触界面倾向东南侧;第四系残坡积物，位于厂房东南侧的陡坎下部，厚度为0～2.9 m。场地地质剖面图见图5。

图5 场地地质剖面图Fig.5 Geological profile of the site

3 研究模型的建立

3.1 基本假设

在建立模型时假设:①平面应变状态;②周围岩石为均匀的弹性各向同性材料;③断裂为无厚度的接触面;④岩石服从摩尔-库仑破坏准则，楼房梁单元为弹性。

图6 几何模型Fig.6 Geomertrical model

3.2 模型的几何形状和边界条件

模拟采用二维模型，模型走向为310°。根据现场实地观测，考虑到边界影响以及模拟的目的，建立地质模型见图6。模型底部宽160 m，右侧高137 m，左侧高97 m，房高12 m，宽12 m(图中所示为楼房的侧剖面)，楼梯高2 m，距房主梁2 m。模型中断层倾角54°。

模拟分为两个阶段:特征值分析和地震动力分析。在进行特征值分析时，通过点弹簧定义弹性边界，模型的左右两侧及底部边界均为弹簧边界。对于动力分析，采用Lysmer和Wass(1972)提出的粘性边界 (viscous boundary)。鉴于该模拟的目的是研究断裂在地震作用时对地表建筑物的破坏作用，模拟时只对模型加载地震波时程荷载，不施加其他应力边界条件。

3.3 材料参数及地震动参数

本研究区的岩性为泥盆系养马坝组 (D2y)泥灰岩，表层为第四系冲洪积物和残坡积物。参考《工程地质手册 (1992)》并结合前人大量的数值模拟经验，采用工程地质类比法，确定材料的力学参数如表1。

表1 材料参数Table 1 Parameters of the materials

地震动的数值模拟是通过输入地震记录 (加速度时程记录)来实现的。该区选择距离最近的卧龙地震台在龙门山发生地震时的监测资料进行动力模拟。地震作用全过程历时135 s。

4 模拟结果分析

4.1 最大主应力分布规律

模拟结果表明:最大主应力值沿断层自下而上逐渐降低，但在断层上盘，厂房下伏的基岩 (泥质灰岩)中出现增大的现象。而在厂房下地基 (第四系冲洪积物)中应力值出现明显的减小趋势，这是由于第四系冲洪积物与基岩泥质灰岩的材料参数相差很大，在地震时有减震的作用。

厂房左侧 (靠近断层一侧)梁底部平面应变单元节点的最大主应力值大于右侧，左侧为83.9 MPa，右侧为71.4 MPa，两者数值相差12.5 MPa(见图7)。

但是在房屋左侧 (靠近断层一侧)第四系冲洪积物的厚度为4.6209 m，右侧 (远离断层一侧)第四系冲洪积物厚度为4.5364 m。按照抗震设计理论，左侧梁单元受到的主应力值应该小于右侧，但是模拟结果却相反。这是因为房屋左侧梁与断层的垂直距离为24.8556 m，右侧梁与断层的垂直距离为84.4065 m，这也就解释了左侧梁单元的主应力大于右侧梁单元主应力值的原因。

4.2 剪应力分布规律

断层上盘剪应力值大于下盘，远离断层，剪应力值逐渐减小。

地基中，厂房左侧 (靠近断层一侧)的剪应力值大于右侧，左侧为18.626 MPa，右侧为17.629 MPa(见图8)。

从梁单元轴向应力曲线图 (见图9)可以看出，梁单元 (图中梁单元由下至上编号，每楼层分为4个单元，下同)底部所受的轴向应力最大，且梁单元的轴向应力自下而上减小，到第二层上部时已经减小到底部的约二分之一，这与厂房底部柱体发生挤压破坏 (见图3、图4)、而上部柱体没有破坏的特征相对应。

从梁单元剪应力曲线图 (见图10)可以看出，左侧梁单元剪应力最大值出现在楼梯与主体梁相交处 (3号梁单元处)，由12.008 MPa急剧增大到12.769 MPa，向上至第一层顶部 (4号梁单元处)时又急剧减小到11.905 MPa。梁单元剪应力的这种变化规律与图3显示的破坏特征一致。同时，可以看出，剪应力均为正值，反映其剪切方向为顺时针转动，与图2—图4显示的柱体剪切破坏特征相对应。

图7 厂房地基最大主应力分布图Fig.7 Distribution of the maximum principal stress in base

图8 厂房地基剪应力云图Fig.8 Shear stress of the base

图9 梁单元轴向应力曲线图Fig.9 Axial stresses of the beam element

图10 梁单元剪应力曲线图Fig.10 Shear stresses of the beam element

根据梁单元轴向应力和剪应力特征及变化趋势，地震时楼体主梁的破坏是轴向应力和剪应力同时作用造成的。

4.3 剪应变特征

断层附近剪应变发生较为复杂的变化，断层上盘和下盘的剪应变值有明显的差别，从地面上剪应变的分布看，厂房区域剪应变值较大。

从厂房底部的最大剪应变图 (见图11)可以看出，在梁的底部出现应变集中区，而且左侧 (靠近断层一侧)大于右侧，左侧为0.03728，右侧为0.03079。

图11 厂房底部最大剪应变分布图Fig.11 Distribution of maximum shear strain in the base

4.4 加速度特征

4.4.1 断裂带上节点水平向加速度特征

断层带上节点水平向加速度值在断层底部最小 (-4.457 cm/s2)，由底部向上增大，增大的幅度逐渐减小。断层位置的加速度明显增大，断层顶部的加速度最大值达638.7 m/s2(卧龙台记录到的最大水平加速度为949.979 cm/s2)。

4.4.2 不同岩性中节点水平向加速度特征

泥灰岩中一节点的最大水平加速度值为313.2 m/s2，冲洪积层中一节点最大水平加速度值为472.6 m/s2，残坡积物中一节点最大水平加速度值为637.8 m/s2，表明在软岩中，加速度明显增大。

5 结论

(1)位于断层上盘的厂房，靠近断层一侧的最大主应力值大于远离断层一侧。

(2)断层上盘剪应力值大于下盘，远离断层，剪应力值逐渐减小;厂房左侧 (靠近断层一侧)的剪应力值大于右侧。

(3)断层附近剪应变发生较为复杂的变化，断层上盘和下盘的剪应变值有明显的差别，厂房底部出现应变集中区，而且左侧 (靠近断层一侧)大于右侧。

(4)加速度值在断层底部最小，由底部向上增加。断层位置的加速度明显增大，断层顶部的加速度最大值达637.8 m/s2。在软岩中，加速度有明显的放大。

(5)厂房所受的轴向应力在底部最大，向上逐渐减小，与厂房底部柱体发生挤压破坏的特征相对应;厂房所受的剪应力值和剪切方向与柱体剪切破坏特征相对应。

(6)断层的形成改变了局部的应力场条件和加速度特征，从而使位于断层上盘的楼房发生差异性破坏，紧邻断层一侧破坏强烈。地震时楼体主梁的破坏是轴向应力和剪应力同时作用造成的。

［1］Somerville P G，Smith N F，Graves R W，et al.Modification of empirical strong ground motion attenuation relations to include the amplitude and duration effects of rupture directivity ［J］.Seismological Research Letters，1997，68(1):199～222.

［2］Somerville P G.Seismic hazard evaluation［C］//Proceedings of 12thWorld Congress of Earthquake Engineering.Auckaland:New Zealand Earthquake Commision，2000:1～16.

［3］Wang G Q ，Zhou X Y，Zhang P Z，et al.Characteristics of amplitude and duration for near fault strong ground mot ion from the 1999 Chi-Chi，Taiwan，earthquake［J］.Soil Dynamics and Earthquake Engineering，2002，22(1):73～96.

［4］Bray J D，Rodriguez-M arek A.Characterizat ion of forward-directivity ground motions in the near-fault region ［J］.Soil Dynamics and Earthquake Engineering，2004，24(5):815～828.

［5］刘启方，袁一凡，金星，等.近断层地震动的基本特征 ［J］.地震工程与工程振动，2006，26(1):1～10.LIU Qi-fang，YUAN Yi-fan，JIN Xing，et al.Basic characteristics of near-fault ground motions ［J］.Earthquake Engineering and Engineering Vibration，2006，26(1):1～10.

［6］杨迪雄，赵岩，李刚.近断层地震动运动特征对长周期结构地震响应的影响分析 ［J］.防灾减灾工程学报，2007，27(2):33～40.YANG Di-xiong，ZHAO Yan，LI Gang.Influence analysis of motion characteristics of near-fault ground motions on seismic responses of long-period structures［J］.Journal of Disaster Prevention and Mitigation Engineering，2007，27(2):33～40.

［7］吴珍汉，张作辰.汶川8级地震地质灾害的类型及实例 ［J］.地质学报，2008，82(12):1747～1752.WU Zhen-han，ZHANG Zuo-chen.Types of seismic and geological hazards caused by the Ms 8.0 Wenchuan earthquake［J］.Acta Geologica Sinica，2008，82(12):1747 ～1752.

［8］谢洪，王士革，孔纪名.“5.12”汶川地震次生山地灾害的分布与特点 ［J］.山地学报，2008，26(4):396～401.XIE Hong，WANG Shi-ge，KONG Ji-ming.Distribution and characteristics of mountain hazards induced by the earthquake on May 12 in Wenchuan，China［J］.Journal of Mountain Science，2008，26(4):396～401.

［9］王连捷，崔军文，周春景，等.汶川 5.12地震发震机理的数值模拟 ［J］.地质力学学报，2009，15(2):105～113.WANG Lian-jie，CUI Jun-wen，ZHOU Chun-jing，et al.Numerical modeling for Wenchuan earthquake mechanism ［J］.Journal of Geomechanics，2009，15(2):105～113.

［10］张永双，石菊松，孙萍，等.汶川地震内外动力耦合及灾害实例 ［J］.地质力学学报，2009，15(2):131～141.ZHANG Yong-shuang，SHI Ju-song，SUN Ping，et al.Coupling between endogenic and exogenic geological processes in the Wenchuan earthquake and example analysis of geo-hazards ［J］.Journal of Geomechanics，2009，15(2):131 ～ 141.

［11］谭成轩，孙叶，吴树仁，等.5.12汶川Ms 8.0大地震后关于我国区域地壳稳定性评价的思考 ［J］.地质力学学报，2009，15(2):142～150.TAN Cheng-xuan，SUN Ye，WU Shu-ren，et al.A consideration on regional crustal stability assessment after Ms 8.0 Wenchuan strong earthquake in China［J］.Journal of Geomechanics，2009，15(2):142 ～150.

［12］侯康明，李丽梅，黄耘.汶川8级大震震害特征和发震断裂运动方式探讨 ［J］.防灾减灾工程学报，2008，28(3):88～93.HOU Kang-ming，LI Li-mei，HUANG Yun.Discussion on the characteristics of earthquake disaster and the movement mode of seismogenic fault of the Wenchuan Ms 8.0 extraordinary earthquake occurrence in Sichuan Province ［J］.Journal of Disaster Prevention and Mitigation Engineering，2008，28(3):88～93.

［13］周荣军，黄润秋，雷建成，等.四川汶川8.0级地震地表破裂与震害特点 ［J］.岩石力学与工程学报，2008，27(11):2173～2183.ZHOU Rong-jun，HUANG Run-qiu，LEI Jian-cheng，et al.Surface rupture and hazard characteristics of Wenchaun earthquake with magnitude 8.0 in Sichaun Province ［J］.Chinese Journal of Rock Mechanics and Engineering，2008，27(11):2173～2183.

NUMERICAL STUDY ON THE HAZARD EFFECT OF A BUILDING NEAR THE FAULT PRODUCED BY THE WENCHUAN EARTHQUAKE

SU Sheng-rui1，ZHANG Yong-shuang2，HAO Li-li1，LI Song1
(1.School of Geological Engineering and Surveying，Chang’an University，Xi’an 710054，China;2.Institute of Geomechanics，Chinese Academy of Geological Sciences，Beijing 100081，China)

Particular destroy phenomena occurred in the main building of Shuimo Silicon Company near the fault produced by the 2008 Wenchuan earthquake.On the basis of site investigation，twodimensioanl finite element model is set up，and the hazard effect of a building caused by the faulting is simulated.The simulation result may be used to explain the mechanism for the formarion of the particular destroy phenomena of the building during the earthquake.It is thought that the formation of fault altered the local stress field and acceleration，thus led to the differential destroy of the building at the hanging wall of the fault which is characetrized by the stronger destroy near the fault.The destroy of the main beam is caused by the co-action of axial stress and shear stress.

fault zone;building;hazard effect;numerical modelling

P315.9

A

1006-6616(2011)04-0381-07

2011-04-27

中国地质调查局项目 (1212010914025)，国家自然科学基金重点项目 (41030749)联合资助。

苏生瑞 (1963-)，男，甘肃庆阳人，工学博士，长安大学地质工程系教授，从事地质工程教学和研究工作。