Co- and post-seismic vertical displacements of Wenchuan MS8.0 earthquake near Beichuan
-
摘要: 利用北川附近震前水准数据,以及2008年和2010年两条一等水准数据, 得到了北川附近同震和震后垂直位移量, 获得如下结果: 北川附近的北川——映秀断裂的同震垂直位移量为4.711 m,断层上下盘最大位错量为5.1 m, 然而, 位于平武县境内的青川断裂同震位错仅有0.064 m. 2008——2010年, 汶川MS8.0地震后北川附近的北川——映秀断裂震后垂直位移速率为5——27 mm/a, 江油县大康镇附近的江油——广元断裂震后垂直位移速率为20.6 mm/a, 而古城区附近的青川断裂垂直位移速率仅为0.2——1.3 mm/a.Abstract: Co- and post-seismic vertical displacements of the Wenchuan MS8.0 earthquake, derived from pre-seismic leveling data and two measurements in 2008 and 2010 along two first-order leveling lines near Beichuan, show the following features: Co-seismic displacement along Beichuan-Yingxiu fault was as large as 4.711 m near Beichuan, where the maximum observed fault offset was 5.1 m. In contrast, the observed co-seismic offset of the Qingchuan fault in Pingwu County was only 0.064 m. From 2008 to 2010, the post-seismic displacement rate was 5mdash;27 mm/a near Beichuan-Yingxiu fault in Beichuan area, 20.6 mm/a at Jiangyou-Guangyuan fault near Dakang, and only 0.2mdash;1.3 mm/a at Qingchuan fault near Gucheng.
-
引言
2008年汶川MS8.0地震发生在青藏高原东缘的龙门山推覆构造带内(徐锡伟等,2008;朱艾斓等,2008).龙门山断裂带位于印度板块向我国大陆NNE向挤压作用的前沿地带,呈NW--SE向延伸,是我国南北地震带的重要组成部分(邓起东等,2002).该断裂带邻区的地质构造复杂,发育有NW向、NE向和近SN向4组活动地块的边界断裂带,其西南侧为鲜水河断裂带,东南侧为龙泉山断裂带和华莹山断裂带,东北侧为西秦岭北缘断裂带,西北侧为东昆仑断裂带(邓起东等,2002).
汶川MS8.0地震的动力来源为青藏高原东缘与华南地块的相对运动,是青藏高原内部地壳物质向东和向北扩展的结果(张培震等,2008;朱介寿,2008).由于应力场状态与强震的发生关系密切,汶川地震后关于利用库仑应力变化对这次地震主断裂及其周边主要断裂的应力变化研究取得了很多成果(Parsons et al,2008;王连捷等,2008;吴小平等,2008;张国宏等,2008;万永革等,2009).结果表明,地震破裂带两端库仑应力明显增强,对未来强震的判断有一定启示.但预测未来强震的发生时间仍很困难,因为目前尚不能对应力场随时间的变化作出判断.构造应力是研究地震及断裂带形成机制的重要指标之一.阮祥等(2010)对四川及其邻区的中小地震震源机制解进行反演,并在此基础上采用力轴张量法计算出5个主要断裂带(区)在汶川地震前后的平均构造应力场.其结果显示:龙门山断裂带南段和鲜水河断裂带南段在汶川地震前后局部应力扰动明显,即震后平均应力场方位扰动值相对变化较大;鲜水河断裂带南段、龙门山断裂带南段和华蓥山断裂带中段在汶川地震后逆倾型地震明显减少.由于断裂带(区)范围划定尚没有一个统一的标准,导致选取断裂带(区)内震源机制解资料的不确定性,从而造成计算结果的不确定性.
地震是在区域构造应力作用下,一定区域内积累的大量应变能达到一定程度时所导致的地壳岩层突然破裂、错动的一种自然现象,地震发生意味着区域构造应力场状态的调整.定量化地给出汶川地震前后龙门山及其邻区构造应力场方向和大小随时空的变化特征,对于汶川地震的地质构造背景、板块动力学机制和孕震规律的研究都具有重要的意义.
目前通常采用地震震源机制解中的P,B,T轴来分析区域应力场,但已有研究结果表明,地壳构造应力场的方向具有区域性变化,存在着应力变异带,简单地利用震源机制解中的P,B,T轴分析区域构造应力场,很难确定应力变异带,也很难准确划分构造应力分区(阚荣举等,1977;成尔林等,1982;许忠淮,戈澍谟,1984;许忠淮等,1989;谢富仁,李宏,1995).这主要是因为震源机制解中的P,B,T轴并不等同于构造应力场的3个主应力(Yamakawa,1971),其只能对构造应力场进行大致的方向性描述.20世纪80年代,利用震源机制解反演区域构造应力场方法(Gephart,Forsyth,1984;Michael,1987;Rivera,Cisternas,1990;Horiuchi et al,1995;Plenefisch,Bonjer,1997)的提出,为利用地震资料求取区域平均应力场提供了又一个有效途径.该方法目前最常用的有Gephart&Forsyth(1984)方法和Michael(1987)方法,二者在误差分析描述、测量标准及应力模型方面存在差别,但其基本假设是相似的.Gephart&Forsyth方法通常可提供较为精确的应力方向,尤其是采用高质量资料时结果更佳,缺点是对资料精度要求较高,若资料达不到要求,则不能进行计算;Michael方法则利用质量较差资料时也能给出结果,且能提供较为合理的误差分析,但结果不能提供应力相对大小R值.
本文利用2000年7月-2009年6月龙门山及其邻区的震源机制解资料,采用Gephart&Forsyth(1984)方法反演汶川地震前后龙门山断裂带及其周边地区的构造应力场以及不同时窗的构造应力,进而分析该区域的动力学背景和孕震环境.本文将重点放在汶川地震发生前后龙门山断裂带及其邻区的构造应力状态,旨在通过研究该区域构造应力场的时空变化特征,揭示汶川地震的孕震时间尺度和空间尺度.
1. 基础资料和方法
1.1 基础资料
本文收集整理了2000年7月-2009年6月龙门山及其邻区(25.8°N-34.2°N,99°E-107.3°E)M2.5-6.0地震的460个震源机制解,其中震前316个,震后144个.采用崔效锋等(2005)对震源断层错动类型的划分方法,对上述460个震源机制解进行分析,结果显示正断层型96个,右旋走滑型161个,左旋走滑型68个,逆冲型(兼具左旋、右旋逆冲)135个,各类型所占比例分别为21%,35%,14%和30%.基于该震源机制解数据,绘制出龙门山区域的震中分布图,如图 1所示.
表 1列出了2000年7月-2002年12月部分震源机制解资料.可以看出,汶川地震前,孕震时段右旋走滑型占主导,逆冲型在临震、震后迅速增加,正断型分布则较为均匀.
表 1 本文反演所用的2000年7月—2002年12月部分震源机制解资料Table 1. Partial source mechanism solution data from July 2000 to December 2002 used for inversion in this study序
号发震时间 发震地点 震源
深度
/kmMS 节面Ⅰ 节面Ⅱ P轴 T轴 N轴 年-月-日 时:分 北纬/° 东经/° 滑动角
/°倾角
/°走向
/°滑动角
/°倾角
/°走向
/°方位
/°倾角
/°方位
/°倾角
/°方位
/°倾角
/°1 2000-07-17 19:15 29.42 104.84 7 3.8 -48 83 290 -170 43 193 232 64 344 53 116 48 2 2000-09-21 16:32 32.95 104.62 1 3.6 -59 59 292 -130 43 161 224 82 330 28 130 64 3 2000-09-23 15:00 31.23 100.70 5 4.3 92 42 129 88 48 312 66 3 220 87 311 89 4 2001-01-06 10:28 29.48 105.87 8 3.5 41 46 306 128 62 67 106 35 3 81 268 57 5 2001-02-14 15:28 29.40 101.08 8 5.1 153 44 352 50 72 242 14 43 123 73 227 52 6 2001-02-23 08:02 29.42 101.10 6 6.0 161 44 347 47 77 243 13 48 122 70 231 49 7 2001-02-23 10:12 29.43 101.03 5 4.2 0 64 340 -154 90 250 118 72 22 72 250 26 8 2001-03-01 13:19 29.47 101.05 5 3.5 9 63 325 153 82 59 105 66 10 78 254 27 9 2001-03-07 19:52 27.75 102.75 6 4.2 14 79 313 169 76 46 89 72 180 88 277 18 10 2001-03-11 22:32 29.40 101.05 5 3.5 -6 59 339 -149 85 246 117 72 18 65 238 32 11 2001-04-19 12:56 31.67 102.28 5 3.8 3 61 352 151 88 84 132 68 34 72 268 30 12 2001-04-20 04:20 32.50 102.03 5 3.6 -2 67 340 -157 88 249 117 75 22 73 245 23 13 2001-04-24 19:15 28.50 102.13 5 3.5 178 23 124 66 89 32 144 49 281 50 31 67 14 2001-04-27 14:40 31.33 100.93 5 3.9 13 81 321 171 77 53 97 75 188 87 288 15 15 2001-06-23 10:49 29.45 105.48 11 4.9 -33 39 21 -124 70 264 149 72 35 38 251 58 16 2001-08-06 02:12 30.43 103.28 5 4.2 -68 72 103 -139 28 336 30 66 164 33 290 69 17 2001-08-31 23:17 29.88 102.77 4 3.6 70 27 346 100 64 144 214 21 61 71 328 81 18 2001-09-05 10:58 30.47 101.62 5 3.5 171 58 62 33 83 327 100 62 199 73 316 34 19 2001-09-23 06:35 31.93 102.03 2 3.9 -119 27 155 -76 67 304 224 70 10 24 129 77 20 2001-09-28 06:32 28.37 104.72 11 4.3 179 50 315 40 89 225 352 62 97 63 224 40 21 2001-11-15 23:00 29.43 101.03 5 4.0 -62 37 11 -109 58 225 121 79 359 20 215 74 22 2001-12-06 03:05 29.15 104.82 14 3.8 73 19 209 95 72 12 93 28 286 63 193 85 23 2001-12-23 06:42 29.13 104.60 2 4.1 -43 14 226 -101 80 94 355 55 197 37 92 80 24 2002-01-26 09:14 31.42 103.07 5 3.8 -145 90 215 0 55 305 85 66 344 66 214 35 25 2002-02-09 17:36 28.07 104.93 5 3.6 -133 16 2 -79 79 138 57 57 215 35 321 79 26 2002-03-03 08:23 28.68 102.23 5 4.4 -86 35 20 -93 55 205 113 80 306 10 203 88 27 2002-03-03 07:06 32.37 100.97 5 3.6 42 90 115 180 48 205 241 62 348 62 115 42 28 2002-04-10 14:30 27.52 102.45 4 4.6 -64 81 210 -161 27 103 141 59 273 43 34 65 29 2002-04-29 20:54 28.42 104.70 2 3.8 -145 51 196 -45 63 310 251 83 349 40 155 51 30 2002-04-30 12:56 29.42 105.53 6 3.8 -2 49 64 -139 89 333 206 64 100 61 331 41 31 2002-05-23 22:26 27.47 102.75 5 4.0 133 89 183 1 43 92 239 58 127 60 4 47 32 2002-05-23 11:51 30.67 102.70 7 3.9 -99 24 193 -86 67 3 276 68 86 22 185 86 33 2002-05-31 11:53 30.30 104.18 5 4.4 101 21 206 86 70 40 130 25 302 66 33 87 34 2002-06-01 14:18 29.42 100.98 5 3.6 165 51 340 40 78 240 12 54 116 72 227 42 35 2002-06-21 15:05 27.65 102.73 6 3.6 140 83 186 10 50 90 236 57 131 69 15 41 36 2002-07-10 07:17 31.55 99.18 5 4.2 -98 35 316 -84 55 127 41 80 196 11 310 85 37 2002-07-17 01:15 31.62 104.02 12 3.7 -102 33 13 -82 57 179 95 78 245 13 3 84 38 2002-08-04 20:29 31.68 103.18 5 3.5 142 71 99 24 55 355 143 50 44 80 302 42 39 2002-08-05 02:16 31.54 102.43 5 3.8 48 88 110 177 42 202 234 56 347 60 108 48 40 2002-08-08 19:43 30.87 100.02 5 5.3 -162 54 312 -38 76 53 358 76 99 53 250 40 41 2002-08-17 03:54 27.84 100.82 5 4.5 6 46 141 136 85 235 288 57 179 64 60 45 42 2002-09-11 00:34 32.93 104.83 3 3.9 -160 81 303 -10 70 36 171 83 78 69 279 22 43 2002-09-12 14:31 29.50 105.68 5 3.5 -87 25 34 -91 64 218 126 70 310 20 216 89 44 2002-09-13 13:37 30.12 102.13 5 3.6 -160 42 121 -49 77 226 165 69 277 47 57 51 45 2002-09-17 04:47 32.63 102.67 5 3.8 90 38 160 90 52 340 71 7 250 83 340 90 46 2002-11-20 21:46 33.92 102.82 5 4.0 -169 79 218 -11 80 310 264 90 354 75 173 15 47 2002-11-25 14:12 30.83 99.90 5 4.7 -7 42 187 -132 86 92 330 62 218 55 89 48 48 2002-11-27 01:55 30.97 103.37 5 3.9 -105 24 112 -83 66 276 191 69 353 22 99 84 49 2002-12-12 09:10 33.10 103.93 5 4.1 175 38 252 52 87 158 281 53 37 59 155 52 1.2 原理和方法
应力张量是一个对称的二阶张量,一般用6个独立参数来描述,可以分解为球应力张量(张量的对称部分)和偏应力张量(张量的偏量部分).但张量的分解并不是唯一的,需要在一些合理的物理假设下进行.例如:在一定时空范围内,研究区内的构造应力场是均匀的;单次地震在断层面上的滑动方向即为应力张量在该断层面上的剪切应力方向.当满足这些条件时,在一组震源机制解资料的时间跨度内,研究区域内的偏应力也是均匀的,于是可以寻找一个与该组观测资料最符合的应力张量.
由于观测资料只是表明应力在一个平面上的滑动方向,所以无法确定应力的绝对值.根据震源机制解资料来反演构造应力场时,我们只能确定应力张量6个独立参数中的4个,即3个主应力σ1,σ2,σ3的方向和主应力比值R.R值介于0-1之间,其大小反映了中等主应力σ2的相对大小.当R接近1时,表明σ2与σ3比较接近,即双轴偏压力;当R接近0时,表明σ2与σ1接近或σ3与σ1相差较大,即单轴偏压力(康英等,2008).
基于Gephart&Forsyth(1984)方法所得到的应力比值R可表示为
该参数与Michael(1987)方法所得到的R值意义相近.
将震源机制解的P轴、T轴的方位和倾角等作为应力反演的输入数据,反演过程中将震源机制解的两个节面都考虑为可能的断层面,并从中选择最合适的断层面.通过网格搜寻法,将断层滑动方向与最大剪应力方向之间的最小残差作为目标函数,寻找一个应力张量,使该应力张量在这组地震断层面上的剪切应力方向与滑动方向之差的总和为最小值.4个应力参数的任何一种组合都称为一个应力模型,其中与观测数据拟合得最好的模型为最优模型.在确定应力张量的4个独立参数的同时,也确定了每一次地震的断层面,并得到最优模型的残差,不同方法之间的区别就在于计算残差的算法不同.
通过反演断层数据得到应力,当设计结果表达为区域(几何)一致时,基本的观测结果就取决于断层面坐标系.在自然笛卡尔坐标系上建立"断层几何",如图 2所示.图中,x′1方向为断层面极点方向,x′2方向为沿着两个断层面的交点方向,x′3方向为滑动方向.该坐标系与Bott(1959)所定义的断层面的法向、倾角、走向大不相同.主应力方向由一套第二自然坐标系定义.为了考查单一应力状态在这两种坐标系中的关系,考虑以下张量的变形形式:
图 2 断层和主应力方向的两套笛卡尔坐标系投影图x1-x2-x3为断层坐标系,x′1-x′2-x′3为主应力坐标系;cos-1β31是两套坐标系转换的方向余弦Figure 2. Two sets of Descartes coordinates projection chart for fault and principal stress directionsx1-x2-x3 is fault coordinate system,x′1-x′2-x′3 is principal stress coordinate system,cos-1β31 is the direction cosine of the two systems式中,βik是转换矩阵(β实际上是坐标变换的方向余弦,这些方向余弦在应力是张量情况下形成矩阵),σ′ij是单个断层面上的应力,σkm是区域应力(这里使用爱因斯坦求和约定)(Gephart,Forsyth,1984).
Gephart和Forsyth(1984)提出使用建立好的格点进行计算.依照应力模型与观测资料的误差度量方式,以拟合残差值来表示,分为3种数值计算方法,即极轴旋转法、趋近法和极小旋转法,其中极轴旋转法已比较少用.本文采用的Gephart&Forsyth方法首先利用趋近法对应力模型进行初始判断,再使用极小旋转法求取最佳主应力方位.
2. 汶川地震前后龙门山及其邻区构造应力主方向空间分布
基于2000年7月-2008年5月11日龙门山及其邻区316个MS2.5-6.0地震震源机制解资料,采用Gephart&Forsyth反演方法,得到了汶川地震前龙门山及其邻区的构造应力主方向,如图 3a所示.图中震源断层错动类型的划分采用崔效锋等(2005)对断层类型的划分方法.可以看出:鲜水河断裂带及其北部地区构造应力的最大主应力σ1方位为NNW--NW,震源断层错动类型为走滑型;青海班玛-四川金川之间存在一个异常带,该异常带内震源断层错动类型为正断型,σ1方位呈NNW;位于鲜水河断裂带与龙门山断裂带交汇处的四川康定、金川也同样存在一个比较明显的异常带,该异常带基本沿龙门山断裂带分布,σ1方位由NW迅速转变为近EW向,再由NW转变为NE向,震源断层错动类型完全转变为一系列正断型.这表明该区域的构造应力场发生了改变.
图 3b给出了反演所得到的汶川地震前龙门山及其邻区构造应力主方向的方差分布.可以看出,整个区域应力场方差为0.1-0.35,在该区域中部沿龙门山、龙泉山、大凉山及马边-盐津断裂的方差较小,特别在龙门山附近方差仅为0.10.这可能是由于该区域地震震源机制解资料较为丰富所致.
基于2008年5月12日-2009年6月龙门山及其邻区144个MS2.5-6.0地震震源机制解资料,采用Gephart&Forsyth反演方法,得到了汶川地震后龙门山及其邻区构造应力主方向,如图 4a所示.
对比图 4a与图 3a可以看出:鲜水河断裂带及其北部地区明显存在的异常带消失,该区域构造应力σ1方位变为NW,震源断层错动类型变为走滑型;龙门山断裂带及其周边范围的异常带仍然存在,但构造应力σ1的异常带向NE方向扩张,而在南部汶川地震震中附近的异常带消失.从整个区域来看,汶川地震前后沿安宁河断裂、大凉山断裂、马边-盐津断裂和小江断裂(主要范围为102°E-104°E、26°N-31°N)的区域构造应力场仍然存在异常带,构造应力场变化不大.
图 4b为反演所得到的汶川地震后龙门山及其邻区构造应力主方向的方差分布图.可以看出,整个区域应力场方差为0.2-0.36,区域中部方差较小,龙门山附近仅为0.20.
3. 龙门山及其邻区构造应力变化
3.1 龙门山及其邻区构造应力反演
汶川MS8.0地震动力源涉及的范围与多个Ⅰ级构造单元相关,孕震跨越的时间尺度也应该较长,因此以龙门山及其邻区为对象,在该范围内反演整个区域构造应力场的综合解,从而定量地给出龙门山及其邻区2000年7月-2009年6月不同时段的构造应力场方向和相对大小,这对于研究汶川MS8.0地震的孕震时间尺度和动力学背景,以及构造应力场变化特征与地震活动的相关性具有重要意义.
基于四川及其周边地区460个M2.5-6.0地震震源机制解资料,对2000年7月-2009年6月地震震源机制解进行时间分段.考虑到一定时空范围内的资料均匀性,本文以1年为1个时间窗,半年为1个滑动步长,将构造应力场的综合结果列于表 2.可以看出,2000年7月-2009年6月共分为16个时窗,其中第1-9,10-13和14-16时窗的地震数目分别为21-35,42-70和137-154.基于16个时窗的震源机制解资料,利用震源机制应力反演软件包(Gephart,1990),反演了龙门山及其邻区各个时窗的区域构造应力场综合解,如表 2所示.图 5给出了第2和第3时窗内龙门山及其邻区构造应力场的综合解示意图.
表 2 2000年7月—2009年6月龙门山及其邻区不同时窗构造应力场的综合解Table 2. Comprehensive solutions of the tectonic stress field in Longmenshan and its adjacent regions in different time windows from July 2000 to June 2009时窗
序号时窗范围 窗内地
震数目综合断
层类别R值 μ 仰角/° 方位角/° σ1 σ2 σ3 σ1 σ2 σ3 1 2000-07—2001-12 23 未知 0.40 0.20 40 49 4 268 80 174 2 2001-06—2002-06 21 走滑 0.60 -0.20 32 53 17 43 189 302 3 2001-12—2002-12 28 正断 0.30 0.40 65 6 25 135 32 299 4 2002-06—2003-06 24 正断 0.20 0.60 55 25 23 151 282 24 5 2002-12—2003-12 26 逆冲 0.10 0.80 0 3 87 117 27 215 6 2003-06—2004-06 30 走滑 0.60 -0.20 29 61 1 18 196 288 7 2003-12—2004-12 31 未知 0.20 0.60 49 32 22 281 144 39 8 2004-06—2005-06 31 正断 0.80 -0.60 77 8 10 327 97 188 9 2004-12—2005-12 35 未知 0.40 0.20 40 45 18 116 262 10 10 2005-06—2006-06 42 未知 0.80 -0.60 27 44 35 136 254 26 11 2005-12—2006-12 49 逆冲 0.10 0.80 3 1 87 297 27 131 12 2006-06—2007-06 46 逆冲 0.20 0.60 2 5 84 103 13 219 13 2006-12—2007-12 70 走滑 0.10 0.80 31 59 1 270 92 0 14 2007-06—2008-06 154 逆冲 0.90 -0.80 6 6 81 2 271 138 15 2007-12—2008-12 154 逆冲 0.90 -0.80 2 4 85 2 272 118 16 2008-06—2009-06 137 逆冲 0.90 -0.80 2 4 85 2 272 118 由综合滑动反演断层类别(表 2)可知,未知类型4个,正断型3个,走滑型3个,逆冲型6个.前期逆冲型不占主导,在第11时窗之后,小震综合断层类别趋同,逐步与主震一致.这也验证了上文所述的微破裂贯通,汶川进入临震的过程分析.
3.2 构造应力时间变化特征
3.2.1 R值变化特征
结合3.1节的构造应力反演结果,图 6给出了R值随时间的变化趋势.可以看出:R值在第11-13时窗之间为0.1-0.2,表现为单轴偏压;从第13时窗开始,R值突变为0.9,迅速变为双轴偏压.
应力洛德参数μ′是描述应力偏量的特征量,用以确定应力偏量在π平面的位置.其代表某一点的3个主剪应力的相对比值大小,不同取值表征不同应力状态特征(杨海波,2011).R值与μ′的关系为(廖红建,2006)μ′=1-2R,μ′≤±1.
图 7给出了应力洛德参数μ′随时窗的变化趋势.可以看出,在第1-10时窗μ′值在-0.6-0.8之间波动,在第11-13时窗μ′值仅在0.6-0.8之间呈微小变化,但第13-14时窗μ′值出现大幅变化,由0.8迅速变为-0.8,这说明区域应力场力学特性出现了反转,由拉伸迅速转变为剪切,进而转变为压缩,这种应力场特征与汶川逆冲推覆构造产生的应力场一致.
3.2.2 主应力方位角和仰角的变化特征
图 8给出了3个主应力(σ1,σ2,σ3)的方位角和仰角随时窗的变化趋势.从图 8a可以看出:3个主应力方位角在第10时窗之前变化比较复杂,在第4时窗前σ1方位角与σ2呈反向变化,之后基本同向变化;σ2方位角与σ3则在第3时窗前同向变化,其后基本反向变化;在第10时窗后σ1方位角与σ2基本同向甚至变化一致,σ3方位角则与σ1和σ2呈反向变化.
从图 8b可以看出:3个主应力仰角的变化更加复杂,在第4-5时窗和第10-12时窗σ1的仰角与σ2同向变化,而σ3则与σ1和σ2的变化方向相反;在第12-13时窗,3个主应力的仰角方向均发生反转;第13时窗之后σ1的仰角与σ2变化明显趋向一致,其仰角几乎相近,可见第13时窗确实是一个重要的时段.
4. 讨论与结论
4.1 地震活动性空间分布特征与构造应力场关系
根据收集整理2000年7月-2009年6月龙门山及其邻区(25.8°N-34.2°N,99°E-107.3°E)460个MS2.5-6.0地震震源机制解资料,绘制出该区域的地震震中分布图(图 1).可以看出,在该时段内中强地震主要沿鲜水河断裂带中南段、安宁河断裂、则木河断裂、华蓥山断裂等分布,而沿龙门山断裂带有零星小震分布,但未发生中强地震.
根据汶川MS8.0地震前后龙门山及其邻区构造应力主方向及其方差分布可知:汶川地震前后,沿安宁河断裂、大凉山断裂、则木河断裂、马边-盐津断裂和小江断裂等地方,构造应力存在的异常带变化不大,而沿鲜水河断裂带及其北部地区明显存在的异常带消失,且该异常带上区域构造应力σ1方位变为NW,震源断层错动类型变为走滑型;沿龙门山断裂带的异常带虽然还存在,但构造应力最大主应力σ1异常带向NE方向扩张,在其南部的汶川地震震中附近该异常带消失.
鲜水河断裂带是川滇块体与巴颜喀拉块体的边界带,龙门山断裂带是巴颜喀拉地块与凉山次级地块的边界带.从构造应力角度分析,沿鲜水河断裂带及其北部地区明显存在的异常带变化,可能证明汶川MS8.0地震的动力源为青藏高原的巴颜喀拉地块,汶川地震后,其能量释放完毕,力源消失.而龙门山断裂带异常带的存在,可能是由于巴颜喀拉地块受凉山次级地块的阻隔,使得这一地区构造应力最大主应力σ1方向发生转向,由NW变为近EW,沿龙门山断裂带变为NE向所致.
4.2 地震活动性时间分布特征与构造应力场关系
地震滑动累计速度、滑动累计频度标准差及滑动释放能量相对值均是表征地震活动性时间分布特征的指标.基于2000年7月-2009年6月的地震资料,图 9给出了2000-2009年龙门山及其邻区地震滑动累计速度、滑动累计频度标准差及滑动释放能量相对值随时窗的变化趋势.滑动累计频度标准差表征滑动累计频度偏离累计频度均值的程度,该数值反映了累积频度时间的不均匀性,本质上与累计速度相关,从图中可以看出该值走势与累计速度基本一致,只是幅度有所差别.
地震滑动释放能量相对值是用一个时窗中的地震释放能量之和除以1012 J得到,只是为了便于在一个坐标系中比较.从图 9可以看出:在第1和第2时窗,即2000-2001年,该区域能量释放较高,之后迅速下降并略有恢复,之后长期保持中等水平;在第13时窗之后出现与累计速度相反的变化趋势,滑动能量释放迅速下降,地震发生速度在增加而能量释放却在减小,这说明能量主要是以累计的形式存在并累积到主震中去;同样,地震滑动累计速度在第13时窗出现迅速上升状态,第13时窗即2006年12月-2007年12月,正是汶川地震临震前时段,该时段出现较为密集的小震.这些地震势必是岩石内微破裂贯通的宏观表征,否则能量应与释放速度同步.就岩石破裂机制而言,该时段则是岩石内部微破裂迅速贯通,逐步形成明显断裂的过程.从这个角度看,汶川地震在2006年之后才进入临震阶段,这对于研究地震时间趋势特征有一定意义.
2000年7月-2009年6月龙门山及其邻区构造应力场变化特征显示,R值在第11时窗以后由单轴偏压迅速变为双轴偏压,自第13时窗开始,应力从较为明显的拉伸和剪切转变为以压应力为主.这与汶川地震逆冲推覆的特征一致,说明该区域综合应力从较为明显的拉伸和剪切转变为以压应力为主.这也符合汶川地震逆冲推覆特征.而该时窗之前应力场主要还是以拉伸和剪切为主,由此可见第13时窗是进入临震的时间节点,应力场特征的突变与上述累计速度、能量释放进入临震特征吻合.从3个主应力的方位角和仰角随时间变化的趋势可以看出,在第10时窗之前比较复杂,之后σ1与σ2基本同向甚至变化一致,σ3则与σ1和σ2反向,这些均表明第13时窗是一个应力开始改变的时间节点.汶川地震临震前,应力场作用方向会有一个与主震反向的改变,也就是说区域应力场的相对稳定和同向性在巨大地震临震前并不成立(在应力场空间分布上会有所体现),恰恰相反还会有一个显著的逆向趋势,这是否是大震前的固有特征尚需进一步研究.
汶川MS8.0地震前后,沿鲜水河断裂带及其北部地区明显存在的异常带消失.在该异常带上,区域构造应力σ1方位由NNW变为NW,震源断层错动类型由正断型兼走滑型变为走滑型;而沿龙门山断裂带存在的异常带,其构造应力最大主应力σ1方向仍然为近E--W和NE向,但其范围向NE扩张,而南部汶川地震震中附近异常带范围在收缩.
汶川MS8.0地震前后,龙门山及其邻区构造应力场方差分布显示整个区域应力场方差均变大,变化幅度约为0.01-0.1,这可能是由于震后资料较震前资料丰富所致,震后误差范围仍然在可信度范围.
2000年7月-2009年6月龙门山及其邻区构造应力场变化特征表明,第13时窗(2006年12月-2007年1月)是汶川MS8.0地震发生前的关键时段,构造应力在该窗口开始发生显著变化.在第13时窗,无论是构造应力场重要特征参数R值和应力洛德参数μ′,还是3个主应力的方位角和仰角均发生明显变化,显示出构造应力环境由以拉伸及剪切应力为主变为以压应力为主,这种构造应力环境是逆冲推覆构造形成的力学基础.
-
邓志辉, 杨主恩, 孙昭民, 陈桂华, 马文涛, 孙谦. 2008. 四川汶川MS8.0级地震北川—映秀地表破裂的复杂现象[J]. 科学通报, 53(20):2509——2513.
董运洪, 陈长云. 2009. 汶川8.0级地震北川附近同震垂直位移特征[J]. 大地测量与地球动力学, 29(6):28——31.
杜方, 闻学泽, 张培震, 王庆良. 2009. 2008年汶川8.0级地震前横跨龙门山断裂带的震间形变[J]. 地球物理学报, 52(11):2729——2738.
黄媛, 吴建平, 张天中, 张东宁. 2008. 汶川8.0级大地震及其余震序列重定位研究[J]. 中国科学: D辑, 38(10):1242——1249.
焦青, 杨选辉, 许丽卿, 王博. 2008. 汶川8.0级地震前后龙门山断裂活动特征浅析[J]. 大地测量与地球动力学, 28(4):7——11.
李传友, 叶建青, 谢富仁, 郑文俊, 韩用兵, 刘玉法, 王伟涛, 魏占玉, 赵冬, 马保起, 任俊杰. 2008. 汶川MS8.0地震地表破裂带北川以北段的基本特征[J].地震地质, 30(3):683——695.
李海兵, 司家亮, 裴军令, 付小方, 王宗秀, 李宁, 侯立玮, 吴富蛲, 潘家伟. 2010. 汶川地震(MS8.0)地表破裂过程探讨[J]. 第四纪研究, 30(4):677——698.
冉勇康, 史翔, 王虎, 陈立春, 陈杰, 刘瑞春, 宫会玲. 2010. 汶川MS8地震最大地表同震垂直位移量及其地表变形样式[J]. 科学通报, 55(2):154——162.
滕吉文, 白登海, 杨辉, 闫雅芬, 张洪双, 张永谦, 阮小敏. 2008. 2008汶川MS8.0地震发生的深层过程和动力学响应[J]. 地球物理学报, 51(5):1385——1402.
王椿镛, 韩渭宾, 吴建平, 楼海, 白志明. 2003. 松潘甘孜造山带地壳速度结构[J]. 地震学报, 25(3):229——242.
徐锡伟, 闻学泽, 陈桂华, 于贵华. 2008. 巴颜喀拉地块东部龙日坝断裂带的发现及其大地构造意义[J]. 中国科学: D辑, 38(5):529——542.
徐锡伟, 陈桂华, 于贵华, 孙鑫喆, 谭锡斌, 陈立春, 孙建宝, 陈于高, 陈文山, 张淑萍, 李康. 2010. 5·12汶川地震地表破裂基本参数的再论证及其构造内涵分析[J]. 地球物理学报, 53(10):2321——2336.
张培震, 邓起东, 张国民, 马瑾, 甘卫军, 闵伟, 毛凤英, 王琪. 2003. 中国大陆的强震活动与活动地块[J]. 中国科学: D辑, 33(增刊):12——20.
张培震. 2008. 中国地震灾害与防震减灾[J]. 地震地质, 30(3):577——583.
张希, 王庆良, 唐红涛, 贾鹏. 2011. 汶川地震孕育背景与同震变化的铲形断层位错模拟[J]. 地球学报, 32(2):189——194.
赵博, 石玉涛, 高原. 2011. 汶川MS8.0地震及余震序列重新定位[J]. 地震, 31(2):1——10.
中国地震局网. 2009. 汶川8.0级地震震间垂直形变速率场监测结果[EB/OL]. 北京: 中国地震局, 2009——05——07[2009——05——10]. http:∥www.cea.gov.cn/manage/html/.
中国国家标准化管理委员会. 2006.国家一、二等水准测量规范(GB/T 12897——2006)[S].北京:中国标准出版社: 9——16.
朱介寿. 2008. 汶川地震的岩石圈深部结构与动力学背景[J]. 成都理工大学学报: 自然科学版, 35(4):348——356.
Burchfiel B C, Royden L H, Vander Hilst R D, Hager B H, Chen Z, King R W, Li C, Lü J, Yao H, Kirby E(著). 2008. 任俊杰(译). 2008. 中国四川2008年5月12日汶川地震的地质与地球物理背景[J]. 国际地震动态, (8):13——21.
Wang Q L, Cui D X, Zhang X, Wang W P, Liu J W, Tian K, Song Z S. 2009. Coseismic vertical deformation of the MS8.0 Wenchuan earthquake from repeated levelings and its constraint on listric fault geometry[J]. Earthquake Science, 22(6): 595——602.
邓志辉, 杨主恩, 孙昭民, 陈桂华, 马文涛, 孙谦. 2008. 四川汶川MS8.0级地震北川—映秀地表破裂的复杂现象[J]. 科学通报, 53(20):2509——2513.
董运洪, 陈长云. 2009. 汶川8.0级地震北川附近同震垂直位移特征[J]. 大地测量与地球动力学, 29(6):28——31.
杜方, 闻学泽, 张培震, 王庆良. 2009. 2008年汶川8.0级地震前横跨龙门山断裂带的震间形变[J]. 地球物理学报, 52(11):2729——2738.
黄媛, 吴建平, 张天中, 张东宁. 2008. 汶川8.0级大地震及其余震序列重定位研究[J]. 中国科学: D辑, 38(10):1242——1249.
焦青, 杨选辉, 许丽卿, 王博. 2008. 汶川8.0级地震前后龙门山断裂活动特征浅析[J]. 大地测量与地球动力学, 28(4):7——11.
李传友, 叶建青, 谢富仁, 郑文俊, 韩用兵, 刘玉法, 王伟涛, 魏占玉, 赵冬, 马保起, 任俊杰. 2008. 汶川MS8.0地震地表破裂带北川以北段的基本特征[J].地震地质, 30(3):683——695.
李海兵, 司家亮, 裴军令, 付小方, 王宗秀, 李宁, 侯立玮, 吴富蛲, 潘家伟. 2010. 汶川地震(MS8.0)地表破裂过程探讨[J]. 第四纪研究, 30(4):677——698.
冉勇康, 史翔, 王虎, 陈立春, 陈杰, 刘瑞春, 宫会玲. 2010. 汶川MS8地震最大地表同震垂直位移量及其地表变形样式[J]. 科学通报, 55(2):154——162.
滕吉文, 白登海, 杨辉, 闫雅芬, 张洪双, 张永谦, 阮小敏. 2008. 2008汶川MS8.0地震发生的深层过程和动力学响应[J]. 地球物理学报, 51(5):1385——1402.
王椿镛, 韩渭宾, 吴建平, 楼海, 白志明. 2003. 松潘甘孜造山带地壳速度结构[J]. 地震学报, 25(3):229——242.
徐锡伟, 闻学泽, 陈桂华, 于贵华. 2008. 巴颜喀拉地块东部龙日坝断裂带的发现及其大地构造意义[J]. 中国科学: D辑, 38(5):529——542.
徐锡伟, 陈桂华, 于贵华, 孙鑫喆, 谭锡斌, 陈立春, 孙建宝, 陈于高, 陈文山, 张淑萍, 李康. 2010. 5·12汶川地震地表破裂基本参数的再论证及其构造内涵分析[J]. 地球物理学报, 53(10):2321——2336.
张培震, 邓起东, 张国民, 马瑾, 甘卫军, 闵伟, 毛凤英, 王琪. 2003. 中国大陆的强震活动与活动地块[J]. 中国科学: D辑, 33(增刊):12——20.
张培震. 2008. 中国地震灾害与防震减灾[J]. 地震地质, 30(3):577——583.
张希, 王庆良, 唐红涛, 贾鹏. 2011. 汶川地震孕育背景与同震变化的铲形断层位错模拟[J]. 地球学报, 32(2):189——194.
赵博, 石玉涛, 高原. 2011. 汶川MS8.0地震及余震序列重新定位[J]. 地震, 31(2):1——10.
中国地震局网. 2009. 汶川8.0级地震震间垂直形变速率场监测结果[EB/OL]. 北京: 中国地震局, 2009——05——07[2009——05——10]. http:∥www.cea.gov.cn/manage/html/.
中国国家标准化管理委员会. 2006.国家一、二等水准测量规范(GB/T 12897——2006)[S].北京:中国标准出版社: 9——16.
朱介寿. 2008. 汶川地震的岩石圈深部结构与动力学背景[J]. 成都理工大学学报: 自然科学版, 35(4):348——356.
Burchfiel B C, Royden L H, Vander Hilst R D, Hager B H, Chen Z, King R W, Li C, Lü J, Yao H, Kirby E(著). 2008. 任俊杰(译). 2008. 中国四川2008年5月12日汶川地震的地质与地球物理背景[J]. 国际地震动态, (8):13——21.
Wang Q L, Cui D X, Zhang X, Wang W P, Liu J W, Tian K, Song Z S. 2009. Coseismic vertical deformation of the MS8.0 Wenchuan earthquake from repeated levelings and its constraint on listric fault geometry[J]. Earthquake Science, 22(6): 595——602.
-
期刊类型引用(1)
1. 经根东, 王振祥. 大地震时房屋裂缝变形机理及抗震构造优化. 四川建材. 2019(12): 63-65 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 1950
- HTML全文浏览量: 1617
- PDF下载量: 205
- 被引次数: 1