汶川5.12地震非稳定发震过程模拟.pdfVIP

汶川5．12地震非稳定发震过程的模拟 王连捷1，崔军文2，周春景1，孙东生1，王薇1，唐哲民2 (1．中国地质科学院地质力学研究所，北京100081；2．中国地质科学院地质研究所，北京 100037) 摘 要：汶川5．12．大地震是世界内陆地区罕见的逆冲型强烈地震。滑坡、泥石流、堰塞 湖等次生灾害严重，造成了严重的人员伤亡和财产损失。人们在珍惜这次万年不遇的巨大地 质灾害所传递的重要构造信息和相关地震资料的同时，不得不思考为什么在普遍认为不具备 诱发强烈地震条件的地区会突然爆发如此强烈的地震。汶川地震破裂带的长度达300km，宽 80 km，震中深度为14～19km，持续时间120秒，释放了大量能量。震源机制比较复杂，且 随着时间和空间发生变化。在地震的开始阶段表现为逆冲为主，后来随着向NE方向延伸， 逐渐转变为以走滑为主。为什么在该地区，能集中如此强大的能量，且成因机制如此复杂， 原因是什么?这些问题的解决，无疑对地震形成机制的研究和实现地震预报具有重要意义。。 为此，我们从断层滑动软化的角度探讨汶川大地震的孕育发生过程和机理。地震是断层 的非稳定破坏过程，是断层和围岩介质组成的系统，在地应力作用下，应力达到断层的强度， 断层软化而突然失稳滑动，储存在断层一围岩介质系统里的弹性能突然释放的结果。因此， 研究地震过程时，只考虑断层本身是不够的，还必须考虑断层周围的岩石介质及其所处的地 应力环境。我们将断层和围岩看成一个完整系统进行研究。将围岩看成弹性介质，断层看成 具有软化特性的弹塑性介质，利用应变软化(位移软化)模型模拟汶川地震的发生。 地震模拟剖面和地应力条件 我们选一剖面进行分析。剖面方向为北西向，与龙门山断裂带走向垂直。剖面穿过震中， 长度180km(图1)，深度30km。 剖面宽度为3km。断层倾向 NW500，倾角40。 + 地震前龙门山断裂带附近进 行过多处地应力测量。据此，建 V口n—wUn—ⅥU凸 畴险B邑 立了模型浅部的应力状态，参考 逆断层应力特点建立了深部应力 状态。据GPS观测结果，模型的 图l地震模拟剖面、边界应力、位移速度和约束条件 两侧位移速度差约为2m／a。 (模型中考虑铅直应力、孔隙压力、水平地应力和位移速度) 因此，在模型的西北边界作用有 o-X一水平应力jo-v一铅直应力；p--孔隙压力；V一位移速度 2m／a水平位移速度(图1)， 同时，作用有水平地应力，且水 平地应力随深度线性增加。 考虑重力和孔隙压力作用，孔隙压力和重力引起的铅直应力同样随深度线性增加。 水平应力、铅直应力和孔隙压力皆作为初始场，由于边界位移速度的作用，位移随时间 不断增加，因而，应力不断增加，直到断层破裂，发生地震。 地震应力降和能量释放 圈2为泣』I地嵌技震时的应力卡¨麻山降。由图H以看㈣随着远场位移的If{㈨，虚力逐 渐^『I大。当麻山达到539MPa叶，戍力突然F降，能链突然释放．此时艇牛地崖。地震造成 9 1 的成力降为1 M?a。钎放的总碰娈能约为2xlO。尔格，。 断层错动和地表位移 图3为rF盘沿断层面的相对错动。最深处断层惜动最小，浅部错动变人．地丧错动 持人．可丛88 m。 罔2№山阡 模拟显小，地震0『起的同震地农伉移其上F盘是不埘称的。上盥位移远人于F盅。上盘 nr的{奇移可选55m，下盘向下的似移不到l米。上衙效应明显。 模拟结果得出如下结论 I、汶Jl地震是存重力、水1r地麻山和扎隙流体压力的耳境r，肯减高蟓向东的侧向挤 7×101 址，四川盆地}且挡．使龙¨山断裂带的地应力咀7 MPa／y的速度缓慢积累．当地应 力达到临界状态时(平均水平应力509MPa叫)，断层失稳滑动，形成地震。 180 2、地震的臆力降为19MPa，地震敛率为0 3％，绝人部分能母