基于InSAR约束的202...同震滑动分布及库仑应力变化_王欣.pdf

第 卷 第期 年月大 地 测 量 与 地 球 动 力 学 ，收稿日期：项目来源：国家自然科学基金（）；武汉引力与固体潮国家野外科学观测研究站开放基金（）；地球内部多尺度成像湖北省重点实验室开放基金（）；中央高校基本科研业务费专项（）。第一作者简介：王欣，硕士生，主要从事 数据处理及地壳形变研究，：。通讯作者：李水平，博士，讲师，主要从事大地测量与地球动力学研究，：。：文章编号：（）基于 约束的 年泸定 地震同震滑动分布及库仑应力变化王欣李水平，宋舜跃合肥工业大学土木与水利工程学院，合肥市屯溪路 号，武汉引力与固体潮国家野外科学观测研究站，武汉市洪山侧路 号，中国地质大学（武汉）地球物理与空间信息学院，武汉市鲁磨路 号，摘要：基于 升降轨影像数据，利用 获取 泸定地震视线向同震形变场。首先利用贝叶斯方法搜索断层的先验几何参数，利用非负最小二乘原理反演断层精细滑动分布，然后根据断层滑动分布参数计算震中附近库仑应力变化，最后利用震间 速度场数据计算发震区震间应变场。结果表明：）泸定地震的同震形变场沿视线向的最大形变量为 ；）泸定地震是一次典型的左旋走滑型地震，断层走向为 ，约 ，沿走向破裂约为 ，倾角约 ，断裂深度主要为 ，最大滑动量约为 ，对应深度为，释放的总地震矩为 ，对应矩震级为 ；）鲜水河断裂带南东段、安宁河断裂带北段和玉龙希断裂中北段处于应力加载状态，未来发生地震的可能性较大；）震源区位于拉张应变和挤压应变的转换区域，该应变转换区可能与多个不同活动块体在该地区的交会有关。关键词：泸定地震；滑动分布；库仑应力变化；震间应变场中图分类号：；文献标识码：据中国地震台网测定，四川泸定县发生 地震，震中位于 、（图），震源深度 。泸定地震是该地区自 年康定 级地震以来的又一次强震，也是近年来造成人员伤亡最严重的一次地震，受到社会 的 广 泛 关 注。、等机构发布了此次地震的震源机制解（表），张龙等、韩炳权等也开展了相关研究。虽然所用数据和方法不同，但研究结果总体认为泸定地震断层的走向为 。此次泸定地震发生在鲜水河断裂磨西段附近，距龙门山断裂带、大凉山断裂带、安宁河断裂带和玉龙希断裂带距离较近，其中鲜水河断裂带、龙门山断裂带和安宁河断裂带交会构成了“”字形断裂带。鲜水河断裂带地处川西高山区，是巴颜喀拉块体的西南边界和川滇块体的北边界，地质构造复杂，始于甘孜东谷一带，向南经过炉霍、道孚、康定、磨西一线，至石棉安顺场一带逐渐消失，全长约 ，总体走向为 ，呈略向 凸出的弧形，是青藏高原东南缘的一条全新世活动强烈的左旋走滑型断裂带。根据构造地貌学，晚第四纪以来，鲜水河断裂磨西段的滑动速率为 ；根据时序 方法，鲜水河断裂带整体滑动速率为 ，闭锁深度为 ；根据 或水准测量方法估计出的鲜水河断裂北西段滑动速率为，南段滑动速率为。世纪以来，整条鲜水河断裂带几乎都发生过地表破裂，且复发频繁，有历史记录的级以上地震共次（图）。已有研究表明，年汶川地震和 年芦山地震后，鲜水河断裂磨西段的库仑应力有一定程度的加载，说明该地区未来地震危险性较高。泸定地区地形复杂，海拔高差大，容易发生滑坡等地质灾害，很难进行常规大地测量以及野外现场考察工作。本文采用 获取 泸定地震的同震形变场，首先利用贝叶斯方法中的均匀滑动模型搜索断层的先验几何参第 卷第期王欣等：基于 约束的 年泸定 地震同震滑动分布及库仑应力变化数，然后采用非负最小二乘方法 获得断层面的同震滑动分布，最后通过计算同震库仑应力变化对区域地震危险性进行评估，计算震间应变场对发震断层构造进行讨论。红色五角星为 震中位置，黑色震源球为 、和 等机构给出的震源机制解，黄色震源球为历史上鲜水河断裂带级以上地震，深蓝色震源球为龙门山断裂带上的芦山和汶川地震，圆圈为余震，蓝色箭头为 速度场图泸定地震区域构造背景 表泸定地震破裂参数 来源经度（）纬度（）震级震源深度 走向（）倾角（）滑动角（）张龙等 韩炳权等 本文 泸定地震 同震形变场本文使用泸定地震前后 升降轨影像数据，参数见表。使用 处理软件 ，采用二轨法进行 数据差分干涉处理，具体流程见图。图为升降轨干涉图和同震形变场。由图（）可见，升轨 同震形变场最大视线向沉降量为 ，最大视线向抬升量为 ；由图（）可见，降轨 同震形变场最大视线向沉降量为 ，最大视线向抬升量为 ；由图（）可见，升轨 同震形变最大视线向沉降量为 ，最大视线向抬升量为 。升降轨沉降量和抬升量存在差异，主要是升降轨观测角度不同所致。本文结果与张龙等和韩炳权等的结果大致相同。需要说明的是，相比于 升轨 图像，升轨 图像存在更大的观测噪声，且存在将一些未发生形变区域错误解缠成有相位区域的情况，因此其可靠性相对较差。为提高后续断层滑动分布反演效率，本文对升降轨同震形变场进行四叉树降采样处理，最终得到 个 升轨 视线向形变观测数据、个 升轨 视线向形变观测数据和 个降轨 视线向数据。表 影像干涉对参数 轨道主影像副影像入射角（）方位角（）空间基线时间基线 图 处理流程 图 升降轨干涉图和同震形变场 大 地 测 量 与 地 球 动 力 学 年月断层滑动分布反演 均匀滑动分布由于泸定地震缺乏发震断层的先验几何参数，因此在进行分布式滑动分布反演前，需要使用均匀滑动模型获得断层倾角、走向、震源深度等参数。本文采用 等开发的 软件包，基于升降轨 数据，使用贝叶斯算法进行断层参数估计。在参数估计过程中，对初始模型参数空间不加以严格约束，如将断层走向设为 、倾角设为 ，可以保证搜索结果的可靠性。迭代 次后，最终得到断层模型参数的后验概率分布见图。由图可见，断 层 最 优 长 度 为 、宽 度 为 、深度为 、走向为 、倾角为 、走滑分量为 、倾滑分量为 ，表明本次地震以左旋走滑为主。这与韩炳权等 和 等给出的结果基本一致（表）。图断层模型参数的后验概率分布 分布式滑动分布由于均匀滑动分布无法得到精细的分布式滑动断层模型，因此采用非负最小二乘法求解断层滑动与地表形变的关系，反演的目标函数为：（）（）式中，为观测值权阵，为格林函数，为每个断层单元上的滑移量，为地表观测值，为平滑因子，为拉普拉斯二阶有限差分算子。由于平滑因子对分布式滑动分布反演结果影响较大，因此固定其他参数以测试平滑因子的大小，最终确定平滑因子为。为构建断层面的精细同震破裂模型，对断层长度和深度分别进行适度扩展，以消除反演结果的边缘效应，并将断层面划分为 的矩形单元。建立分布式断层滑动模型是同震反演的关键步骤。由于 的升轨数据质量存在问题，为测试其对反演结果的影响，对本次地震升降轨数据进行以下种实验：）使用 升轨 和降轨视线向数据进行断层滑动分布反演；）使用 升轨 和降轨视线向数据进行断层滑动分布反演。基于种实验方式分别构建分布式滑动分布模型，然后根据数据模拟效果评估反演结果。以 升轨 和降轨 同震形变场为约束进行滑动分布反演，将断层走向设置为 ，对应的数据模拟结果见图。由图可见，升轨的观测值和模拟值显著不同，存在明显的残差，残差主要集中在震中东侧；降轨的观测值和模拟值吻合较好。升轨数据在震中东侧产生明显残差的原因可能是该地区植被茂盛，且地形梯度较大，使得升轨图像中存在较为明显的大气误差。利用 升 轨 和 降 轨 同震形变场作为约束，反演同震滑动分布，第 卷第期王欣等：基于 约束的 年泸定 地震同震滑动分布及库仑应力变化图升降轨同震形变场观测值、模拟值和残差分布 数据模拟结果见图。由图可见，升轨平均残差为 ，降轨平均残差为 ，残差图中不存在明显的系统误差，表明同震滑动模型能较好地模拟地表 观测数据。图升降轨同震形变场观测值、模拟值和残差分布 图中对应的二维和三维同震滑动分布模型见图。反演结果显示，断层滑动主要以走滑为主，兼少量逆冲分量，此结果与 、等给出的结果大致相同。地震破裂深度主要集中在 ，破裂沿走向延伸约 ，大于 的较大量级滑动主要发生在 深度处，断层最大滑动量为 ，对应深度为，释放的总地震矩为 ，相当于 ，此结果与 计算的结果大致相同。由表可见，不同类型观测资料约束的地震震级和发震断层参数基本一致。余震活动主要发生在同震破裂区的中部、深部和边缘扩展区，总体上看，余震分布与同震破裂具有较好的空间互补关系，说明泸定地区余震活动是由地震同震破裂触发的应力加载所致。图泸定地震滑动分布模型 同震库仑应力变化地震的破裂滑动会使区域应力场发生改变，进而影响余震活动，因此余震活动分布可以验证破裂滑动分布结果的正确性。分析震中 的变化是预测震中附近发生地震可能性的重要方法之一。利用反演的断层滑动分布参数，运用 软件包计算泸定地震附近地区的库仑应力变化。（）式中，为剪切应力变化，为有效摩擦系数，为法向应力变化。将泊松比设为 ，有效摩擦系数设为，分别计算、和 深度处的同震库仑应力变化，结果见图。由图可见，随着深度增加，库仑应力在震中附近及靠南地区发生显著变化，深度范围内大部分余震分布在库仑应力加载区域，符合地震应力传输和同震库仑应力触发机制；小部分余震活动分布在应力卸载区域，表明此次地震未能释放断层面积累的全部应力。余震分布趋势为 向，与断层滑动分布走向大致相同。泸定地震发生后，库仑应力的增加主要集中在破裂断层的西北段、东南段以及西南段，位于鲜水河断裂带、安宁河断裂带北段、玉龙希断裂中北段。上述地区未来发生大 地 测 量 与 地 球 动 力 学 年月地震的危险性很大，尤其需要注意安宁河断裂带北段，该断裂带在 年和 年发生 级地震、在 年发生 级地震，是需要重点关注的地震空区。红色五角星为震中位置，白色实线为断层位置，白色虚线为断层俯视图，黑色实线为断层，黄色圆圈为不同深度的余震分布图泸定地震引起的 变化 震源区孕震特征泸定地震发生在鲜水河断裂带东南端，毗邻鲜水河断裂带、龙门山断裂带、安宁河断裂带的交会区。分析震源区震前长期应变特征，可为识别强震孕育阶段以及预测强震提供重要参考。本文利用球面多尺度小波变换方法，以 速度场为约束，建立泸定地震震源区地壳应变场模型。速度场数据来源于 等 给出的鲜水河地区近 的 观测结果。球面多尺度小波分解时，球面网格的最大分辨率（）对结果的影响较大，根据数据的分布密度，经过多次实验，取 。此时，最大网格尺度为 ，分辨率低于该阈值的区域设定为不确定区域，在计算应变参数时，忽略不确定区域的结果。图（）为鲜水河断裂及邻区主应变率和最大剪应变率分布。可以看出，泸定地震震源区以 向的挤压应变为主，主压应变率量值约为（），明显小于鲜水河断裂中段和北段，表明鲜水河断裂挤压变形从西北段到东南段逐渐减弱，此结果与藏东高原整体的东向挤出背景基本一致。最大剪切应变率基本沿鲜水河安宁河断裂带展布，并且在鲜水河断裂带、龙门山断裂带和丽江小金河断裂带的交会区附近存在一个应变率高值过渡区，该过渡区的位置与泸定地震的位置基本一致。图（）为震源区面膨胀率分布。由图可见，此次泸定地震震源北侧和东侧以挤压应变为主，西南侧以拉张应变为主，震中所在的位置正好位于面应变率从挤压应变到拉张应变的转换区域，该应变转换区的存在可能与多个不同活动块体在该地区交会并产生的相互作用有关。图泸定地震震源区最大剪切应变率和面膨胀率分布 结语）泸定地震升降轨 同震形变场中，升轨最大视线向形变量为 ，降轨最大视线向形变量为 。）基于视线向同震形变场反演的滑动分布模型显示，发震断层以左旋走滑为主，兼少量逆冲，断层走向沿 ，约为 ，沿走向破裂约 ，倾角约为 ，断裂主要集中在 深度处，最大滑动量约为 ，对应深度为。泸定地震释放的总地震矩为 ，对应矩震级为 。）同震库仑应力结果显示，泸定地震使鲜水河断裂带南东段、安宁河断裂带北段、玉龙希断裂中北段处于应力加载状态，其中安宁河断裂带北段库仑应力加载状态明显，未来发生地震的风险较大，需要开展持续监测和危险性评估。）震间应变场表明，此次泸定地震震源区位于拉张应变和挤压应变的转换区域，该应变转换区的存在可能与不同活动块体在该地区的交会有关。第 卷第期王欣等：基于 约束的 年泸定 地震同震滑动分布及库仑应力变化参考文献张龙，周一剑，张旭，等 年四川泸定 地震：一次填补地震空区的陆内强震 ：，（，：，）韩炳权，刘振江，陈博，等 年泸定 地震 同震形变 与 滑 动 分 布 武 汉 大 学 学 报：信 息 科 学 版，（）：（，（）：）潘家伟，李海兵，等鲜水河断裂带色拉哈康定段新发现的活动断层：木格措南断裂地质学报，（）：（，：，（）：）姚鑫，张永双，李凌婧，等青藏高原鲜水河活动断裂带蠕变斜坡地质灾害 识别研究 地质学报，（）：（，（）：），：，：，（）：，（）：，：，（）：，（）：，（）：，（），（）：，（）：，：，（），y ，：，（下转第 页）大 地 测 量 与 地 球 动 力 学 年月 y ，：，（），（），：；：，：，：（上接第 页），：）；），；），；），：；：，；，；，；，：，：，：