地表动态大地测量资料反映的孕震断层变形机制研究

邹镇宇

(中国地震局地质研究所， 北京100029)



地表动态大地测量资料反映的孕震断层变形机制研究

邹镇宇

(中国地震局地质研究所， 北京100029)

构造地震孕育、 发生、 震后调整过程中应力应变积累、 释放必然伴随有相应的地壳形变发生。 利用大地测量动态资料研究与地震过程密切关联的地壳形变时空动态特征是地震预测研究的重要途径之一。 由于构造地震发生在活动断层上(通常在10 km以下的深部)， 因此， 如何建立逼近真实的地表形变动态变化与孕震断层深浅部应力应变状态的关系， 以揭示孕震断层变形机理是地震机理与预测研究中至关重要的科学问题。 开展这方面的研究具有重要的科学意义。 本论文首先研究了利用GPS观测资料获取与地震孕育关联的地表形变动态变化特征方法， 进一步研究了深浅部地壳运动状态与断层应变积累之间的关系。 主要研究内容、 结果和取得的认识有：

(一) 地壳变形动态特征分析方法研究及其在南北地震带的应用。 研究解决了速度场稳定基准选取和变形频域一致等建立多期可比较速度场的关键技术问题； 研究讨论了微元转动参数的基本特性及其在构造形变分析中的意义。

(1) 利用“中国地壳运动观测网络”和“中国大陆构造环境监测网络”GPS区域站1999～2013年7期复测资料， 建立了南北地震带地区相对于稳定华南地块统一参考基准的多期速度场， 基于该具有可比性速度场， 对2008年汶川8.0级地震前后的地壳形变动态过程进行了分析。

① 为构建各期速度场的统一参考基准， 发展了拟准检定法(QUAD法)剔除观测点群中的不稳定点， 改进了初选指标稳定点选取准则， 有效抑制了初算结果可能出现的基准偏移对稳定点的筛选影响， 提高了结果的可靠性。

② 选择统一且适当的协方差衰减参数， 利用最小二乘配置进行GPS速度场拟合推估， 保证了各期速度场具有相同的变形频域， 并解决了不同期资料测点分布存在差异的问题。

③ 对任意两期网格化速度场求差分获取了这段时间的地表位移动态变化。 汶川地震前后的两期速度场差值结果显示， 在2007～2009年间汶川地震的影响范围较大， 包括祁连地块、 柴达木地块东部都有明显的南东东向运动响应， 但距破裂带更近的鄂尔多斯西南缘响应很小， 可能属于背景应力应变积累水平较高的地带， 但巴颜喀拉地块北边界带东部的东昆仑断裂带左旋剪切的响应明显。 汶川地震震后效应对龙门山断裂带南段表现为显著的应变加载过程， 而鲜水河断裂带则表现为与应变积累背景相反的右旋扭动响应。

④汶川地震发生后， 川滇地块中北部及其边界有近东西向挤压增强， 但整个川滇地块南东向挤出滑移的运动背景并没有增强， 直到2011～2013时段南东向运动才有所增强。

(2) 通过应变张量计算得到的转动参数反映计算单元的有限转动但不包含纯应变信息， 与应变参数配合使用可以更客观地反映地表的变形状态。

① 该特征与剪应变参数不同， 转动参数与坐标方向的选取无关， 值的大小表示微小单元的转动。

② 转动参数与参考基准的选取有关。 当研究多期资料时， 需考虑参考基准的统一， 可采用研究区域无整体旋转基准速度场进行计算。

③ 转动参数反映了微元在变形中主应变轴的偏转。 基于区域无整体旋转基准速度场计算的转动参数可表示最大剪应变在两个方向的不对称， 从而确定区域构造运动决定的实际剪切变形最大方向。

④ 对于纯走滑断层， 震间位移曲线的斜率反映了转动参数的大小。 当断层的变形宽度较窄时， 转动参数非零值的区域较窄； 当断层的变形宽度较宽时， 转动参数非零值的区域较宽。 当区域剪应变较大时， 如果转动参数量值较大， 说明断层的应变积累程度可能较高； 当区域剪应变较大， 而转动参数量值较小， 说明断层的应变积累程度可能较低， 甚至断层处于蠕滑状态。

(二) 一般倾角断层下地表震间/同震位移场的变形特征研究。 通过对特定断层的剖面分析研究， 从力学特性分析了断层倾角对地表位移的影响。 在原有公式基础上推导给出了带断层倾角的震间、 同震的走滑、 倾滑地表位移公式， 该公式能够较好地拟合与断层闭锁关联的地表位移， 与位错模型中的复杂公式相比是一种简化公式， 但已能达到相当高的拟合精度， 便于实测GPS速度场资料做拟合的实际应用。

(1) 对于非直立型走滑断层而言， 震间变形中心一般不位于断层出露地表处， 而是位于断层的滑动段上边沿， 即断层闭锁段与滑动段的分界线在地表投影处。

(2) 无论走滑断层还是倾滑断层的震间形变， 断层闭锁段与滑动段的分界线在地表投影与断层地表出露处之间的距离doffset、 断层闭锁深度d和断层倾角δ之间存在关系tanδ=d/doffset， 此关系揭示了断层的闭锁深度与断层倾角之间的关系。

(3) 由于断层倾角的影响， 震间形变曲线以断层闭锁段与滑动段的分界线在地表投影为中心， 同震时上下盘错动沿断层发生错动， 导致上下盘同震释放的位移不对称。 而当发震断层为盲断层时， 地震错动不达到地表时的特例情况与震间的情况类似。

(4) 基于实际观测数据， 利用带倾角的公式拟合安宁河断裂带的倾角， 得到的倾角结果小于地质考察结果， 原因可能由于实际断层是曲面而非平面， 在地表处断层陡峭， 随着断层深度增加倾角变小。 由于断层闭锁段以下滑动对地表位移的影响， 相当于断层出露地表处与断层闭锁与滑动分界之间存在一条“平面断层”对地表的作用。 利用反正切函数拟合的倾角结果为这条“平面断层”的倾角。

(三) 应用三维数值流形方法研究孕震断层深浅部力学特性。 基于三维数值流形方法的连续与非连续耦合计算的优势， 通过断层切割算法构建了三维模型， 对孕震断层的力学特性进行研究。 设置断层上部闭锁， 下部滑动， 对研究块体侧向和底部加载“动力源”进行走滑剪切数值试验， 给出了不同加载方式下地表位移场分布差异特性。 基于该结果， 对汶川地震后龙门山断裂带和川滇块体东边界断裂带闭锁深度反演的动态变化给出了机理上的解释。

(1) 利用三维数值流形方法的断层切割算法建立了三维块体模型， 块体中的断层设置为上部闭锁， 下部滑动， 对研究区域块体加载“推挤”和“拖曳”两种力源来模拟块体所受其他块体对其的推挤作用和脆性层底部软流物质对其的拖曳作用。 模拟结果显示， 两种动力模式下地表位移场均呈现较精确的反正切函数特性， 与位错解析结果有很好的一致性。

(2) 针对两种不同“力源”输入， 模拟得到的位移分布曲线表现出显著差异， 在一定程度上反映力源的力学特征。 “推挤” 力源的力学特性是沿水平方向传递， 得到的反正切曲线在断层远端相对平直， 曲线远端基本反映“推挤”力源的加载量， 变形宽带较小， 利用反正切函数拟合得到的闭锁深度小于模型设定值。 “拖曳”力源的力学特性是垂直向上传递， 得到的反正切曲线在断层远端略有上翘， 由于力在垂直方向传递时衰减导致曲线远端不能反映“推挤”力源的加载量， 曲线变形宽带较大， 利用反正切函数拟合得到的闭锁深度大于模型设定值。

(3) 实际块体应既受到其他块体对其的推挤作用， 又受到下部软流物质对其的拖曳作用。 将两种力源按不同比例同时加载， 得到的结果表明： 哪种力源增强， 地表形变动态变化就体现出此力源加载下的曲线特征。 基于此认识， 可以对汶川地震之后川滇地区地壳深浅部运动特征和DEFNODE负位错反演动态结果给出可能的解释。 汶川地震之后， 龙门山断裂带南段西侧的巴颜喀拉块体南东向运动增强， 是由于龙门山断裂带脆性层同震释放所导致， 可理解为块体的“推挤”作用， 所以反演得到的龙门山断裂带南段闭锁深度变浅。 上部脆性层的突然加速对下部软流物质提供加载动力。 获得加载动力的巴颜喀拉块体下部软流物质由通道进入川滇块体下部， 对川滇块体上部脆性层提供了“拖曳”力源， 由此计算的小江断裂闭锁深度增加。

(4) “推挤”与“拖曳”力源计算得到断层不同深度的位移结果均表明， 断层的闭锁段以下位移量随深度逐渐加大， 并非位错理论中闭锁与滑动的突变分界而是渐变的过程， 这与Tse和Rice(1986)的模型示意图具有一致性。

(5) 考虑下部软流物质应为渐变的可能性， 对“拖曳”力源进行改进， 使之逐渐加载， 加载量与距断层距离呈现单调递增函数关系。 模拟得到的地表位移结果表明， 渐变“拖曳”力源的加载使得地表位移曲线相比于无渐变“拖曳”力源得到的变形宽度更大， 因此， 实际断层由于可能存在软流物质渐变拖曳而使得地表变形宽度更宽。

(6) 南北地震带中段三维数值流形模拟结果表明， 在弹性本构模型下水平主应力增量和水平主应变增量的方向具有一致性， 但主张应力增量与主压应力增量的比值和主张应变增量与主压应变增量的比值存在差异。 汶川地震前龙门山断裂带处于挤压面应力积累状态， 汶川震源区处于该挤压应力速率高值区的边缘， 以汶川震源区为界龙门山断裂带南段的挤压应力速率明显快于龙门山断裂带北东段。 汶川地震后的2009～2013年龙门山断裂带中北段处于震后调整过程， 而龙门山断裂带南段则处于应力积累状态。 鲜水河断裂带中北段以张性应力增量为主， 而其中南段以张剪应力增量为主。 芦山地震震源区处于龙门山断裂带中北段和鲜水河断裂带中南段最大剪应力快速积累区的弱化地带。 安宁河断裂带的最大剪应力积累速率小于其北侧的鲜水河断裂带南段以及其南侧的则木河断裂带。

总体而言， 本文通过定量描述地表形变动态变化特征， 从理论上研究断层倾角对震间和同震的地表位移影响， 利用三维数值流形方法研究断层深浅部力学特性对地表位移的影响， 分析了地壳动态变形特征及其反映的孕震断层力学特性， 对孕震断层应变积累状态与地表位移的动态变化之间关系有了初步认识。 为今后进一步结合不同期GPS资料， 利用三维数值流形方法研究实际断层不同时间段， 地表形变动态变化的深部动力加载影响， 科学揭示孕震断层闭锁与深浅部应力应变状态奠定了理论基础。

关键词三维数值流形方法； GPS速度场； 拟准检定法； 孕震断层； 力学特性

(作者电子信箱， 邹镇宇： 407124082@qq.com)

中图分类号：P315.72+5；

文献标识码：A；

doi:10.3969/j.issn.0235-4975.2016.07.0012