利用水位同震阶变资料反演汶川8级地震应力场变化.doc

利用水位同震阶变资料反演汶川8级地震 应力场变化 杨明波 马霁 柴金翼 （中国北京100080北京市地震局） 摘要：分析汶川地震的水位同震响应阶变变化发现，地震时应力场的调整与中国大陆板块的背景应力场受力状态基本一致，只是对华北和西南地区有所加强，对西北和东南有所减弱。可见，汶川地震释放积累的应力后，对于震后大陆的应力场变化有一定调整作用。 关键词: 汶川8.0地震；含水层；区域应力场；应力调整；水位同震阶变 0引言 井水位对地震的响应分为波动类和阶变类。波动类响应通常不打破原有的水位变化趋势和过程，仅围绕水位过程线上下振动，即水震波，通常认为井水位的同震效应“水震波”是由地震产生的面波在转播过程中引起含水层渗透率发生变化的结果。阶变类响应则打破了井水位的原有变化过程，或引起井水位突然/快速地上升，或出现水位的突然/快速地下降，而这种变化一旦出现后将持续一定的时间，产生永久性变化。这种永久性变化通常当作“同震阶变”。 对于水位阶变型同震变化，一些学者认为震源区周围范围承压含水层的水位阶变是源于多孔弹性介质对静态应力场的响应，通常认为是震源区应力/应变调整的结果，其阶变变化的幅度反映了应力应变调整的程度，理论上可以用Okada（1990）给出的弹性半无限空间位错模型来加以解释。 （黄辅琼等，2000）对远强震引起的地下水阶变同震响应特征进行比较系统的研究，认为远强震引起的承压井水位、水温及其化学组分异常响应变化与构造特殊部位有关，其空间上沿着构造特殊部位的有序分布可以用（许绍燮等，1980）提出的地壳屈曲理论来给予解释，而这种水位阶变变化可能反映远强震引起的地壳应力/应变的远场调整，也是地壳处于临界状态的反映。 1 利用井水位同震阶变反演含水层应力变化的理论模式 通过分析井水位变化与含水层应力应变变化之间的关系，进而研究井水位与区域应力场的关系，探求地下流体中的地震前兆信息，（张昭栋，1987）利用井水位阶变反演含水层应力变化的理论模式结果如下。 考虑在各向同性承压含水层中的完整井, 该含水层呈水平无限延伸。假设各层的力学性质是各向同性的完全弹性体。在这种理想情况下,当含水层受到一个垂向的应力作用时,含水层发生应变使水井水位发生相应的变化, 水井水位降深同含水层垂直向应力有如下的变化规律。 Δσz=　 其中, 对于含水层的孔隙度和固体骨架的杨氏模量, 我们根据含水层岩石性质和生成的地质时代, 取其一般值或平均值 (见表2不同岩石的扬氏模量)。 表1 常见岩石的孔隙度(J.塔罗勃，1965) Table 1 Common Rock porosity 岩石 平均孔隙度（%） 岩石 平均孔隙度（%） 第三纪石灰岩 20 侏罗纪次生砂岩 13.2 次生石灰岩 11.65 三叠纪次生砂岩 19.3 原生石灰岩 12.5 大理岩 1 白云石 7.7 石英岩 2.1 砾石和卵石 28 板页状页岩 1.1 片麻岩 1.35 硅质页岩 1 花岗岩 0.95 黑页岩 1 第三纪砂岩 15.3 变质页岩 0.16 白垩纪次生砂岩 17.1 表2 常见岩石的扬氏模量(Coates，1978) Table 2 Common Rock Youngs modulus and porosity 岩石 扬氏模量E（1010N/m2） 岩石 扬氏模量（1010N/m2） 石灰石 3.73-6.07 花岗岩 2.07-3.59 砂岩 2.97 凝灰岩 1.38-7.66 粉砂岩 1.31 熔岩 6.21 砾岩 7.38 角闪片岩 8.76 页岩 1.24-2.69 似碧石岩 8.83 石英岩 6.62 含水层内的水的密度ρ取1020kg/m3， 2 震时应力变化量计算 根据上述理论模式，收集全国范围内230口井的观测资料, 选取存在同震阶变且具含水层岩性资料的115口井，进行汶川地震震时应力变化数值计算。表3给出部分井孔资料和计算结果。 表3 汶川8.0级地震时部分井孔应力变化量 Table 3 the basic situation and stress variation of Observation wells that produce water level Step variation on Wenchuan 8.0 earthquake 井名 震中距（km） 岩 性 震时应力变化值（pa） 响应 幅度(m) 安徽安庆 1306 上升 0.336 中生代石英正长斑岩 228871.7 安徽巢湖 1374 下降 0.545 石英砂岩 -105878 安徽歙