活断层形变数据反演

本文主要讨论了利用形变数据反演断层运动位错模型参数的贝叶斯方法。利用现代大地测量方法,尤其是空间大地测量方法获得的地壳形变数据,可以用来反演活动断层的位错模型参数,并且利用反演得到的结果,确定断层运动的闭锁区。作为应用,给出了运用地表观测站观测到的速度值反演圣安德烈斯断层运动位错模型参数的结果,并且据此预测了沿圣安德烈斯断层未来最可能发生地震的断层分段。     

昆仑8.1级地震对青藏块体东北缘地壳水平运动的影响初探

为了探讨昆仑8.1级地震引起的同震形变是否影响了青藏块体东北缘的地壳水平运动,以震后GPS观测资料为数据源,依Okada位错模型原理,计算了地震引起的同震水平位移场,分析了地震对区域地壳水平运动的影响。结果表明:垂直于断层沿震中方向的GPS点受地震影响大,远离断层的测点的同震位移受断层位错的影响小。     

基于GNSS/InSAR/GIS的活动断层地震危险性评估系统

本文从GNSS、InSAR、重力、水准等多种大地测量资料联合确定活动断层地壳形变、应变和应力,以及与GIS相结合的角度出发,讨论了基于库仑破裂准则的活动断层地震危险性评估方法,系统阐述了活动断层地震危险性评估系统的设计。该系统主要利用各种大地测量数据联合确定活动断层区域的位移场,反演计算应变场、应力场、活动断层运动参数以及库仑应力;通过确定活动断层的临界应力值,建立地震危险性评估模型,评估可能发生地震的位置、最大震级、复发时间及其对应的概率;最后借助GIS实现断层区域的地质和构造环境、三维形变数据、应变场、应力场和地震危险性综合评价图的管理与可视化。本文最后讨论了所建系统的具体应用,讨论了汶川地震激发的同震静态库仑应力、应力扰动的持续时间和龙门山地区的地震活动率。     

尼泊尔Mw 7.9级地震同震垂直位移与断层运动模型

2015年4月25日尼泊尔地区发生了Mw 7.9级地震,发震断层位于印度板块与欧亚板块碰撞边界带,此次地震是一次典型的板块逆冲型事件。利用中国境内加密的GPS同震观测资料,融合ALOS-2卫星L波段的InSAR（interferometric synthetic aperture radar）同震形变数据,基于最小二乘方法获得了此次地震的同震垂直位移场。同震垂直位移结果表明,此次地震造成尼泊尔加德满都地区抬升约0.95 m,珠穆朗玛峰地区受地震的影响有所下降,其主峰的沉降量为2~3 cm,中国境内的希夏邦马主峰沉降约为20 cm。地区利用改进的二维弹性半空间位错模型反演了发震断层运动参数,本文模型显示此次地震的断层面破裂宽度约为60 km,平均滑动量达到4 m,相当于Mw 7.89级。     

顾及地壳粘弹性结构的地震断层震后形变反演分析

顾及地壳的粘弹性结构特征,建立了基于粘弹性地球模型的地震断层震后反演模型。利用遗传算法重点对走滑断层和逆断层两种主要地震断层进行了震后形变反演研究,并与球体均匀模型反演进行了比较。研究结果表明:遗传算法能在较大范围内搜索出发震断层的位错参数,反演的结果基本可靠;基于粘弹性地球模型特征建立的地震断层震后反演模型,可较好地反演不同类型的地震断层震源位错参数;如果出现模型参数不清楚的情况,需综合考虑衡量反演结果的准则,找寻其他方法来判断,仅以VTPV最小作为反演准则并不可靠。     

卫星大地测量成像地震周期形变研究综述

近年来,卫星大地测量技术的快速发展为精确测量地壳形变和断层行为提供了前所未有的多维观测数据。它与地震学的结合使得地震周期形变监测的时空分辨率大大提升,为更加深入地研究地震周期过程和机理提供了一个窗口,地震大地测量学应运而生。它能够对地壳运动进行定量描述、对断层活动进行精准建模,从而为洞悉整个地震周期过程的应力应变演化提供科学依据,同时为评估地震危险性、实现地震预测预警提供科学指导。以卫星大地测量观测探究断层形变为主线,分析了断层处于地震周期不同阶段的运动学特征（震间、同震和震后）,回顾了地震大地测量学在震源物理方面的一些重要发现。研究表明,利用卫星大地测量数据判定断层所处地震周期的阶段是实现地震预测的可行思路。     

2016年意大利阿马特里切Mw 6.2地震震源机制InSAR反演

2016年8月24日，意大利中部阿马特里切（Amatrice）地区发生Mw 6.2地震。采用ALOS-2条带模式和SENTINEL-1A宽幅模式的合成孔径雷达（synthetic aperture radar，SAR）数据分别进行SAR差分干涉测量处理，获取了该地震的同震形变场。结果显示，本次地震造成意大利中部地区发生明显的地壳形变，在雷达视线向最大沉降量达19.6 cm。基于合成孔径雷达干涉测量（interferometry synthetic aperture radar，InSAR）和GPS同震形变场数据对此次地震的发震断层进行联合反演，通过改进倾角和平滑系数获取方法，得到了最优滑动分布模型。通过使用单断层模型和双断层模型进行反演可知，双断层模型反演结果优于单断层反演结果，两种模型下反演模型相关系数分别为0.85和0.89，发震断层走向分别为160°和158°，倾角分别为44°和46°，倾滑分布主要位于地下5~7 km，平均倾滑角为-80°，最大倾滑量0.9 m位于地壳深度5 km处，该发震断层是亚平宁冲断带的一部分，为NW-SE向延伸的正断层，断层长约20 km。综合使用地震同震形变场和GPS数据对震源机制进行反演、模拟和分析，获取了高精度的震源参数，可以为分析地震危险性和断层破裂参数等提供数据支持。     

联合地震位错模型和InSAR数据构建2017年九寨沟Mw 6.5地震同震三维形变场

利用哨兵（Sentinel）-1A卫星升、降轨影像,在地震位错模型约束下获取了2017年九寨沟Mw 6.5地震的高质量三维形变场。首先,利用合成孔径雷达干涉测量技术（interferometric synthetic aperture radar,InSAR）提取九寨沟地震升、降轨同震形变场;然后,通过“两步法”反演获取该地震发震断层的几何参数和分布式滑动模型,以此为约束,采用方差分量估计算法联合解算九寨沟地震三维形变场。结果表明,九寨沟地震同震三维形变场以水平位移为主,垂向形变较弱;南北向形变呈拉张趋势,断层上盘向南、下盘向北滑动,最大位移分别为－19.81 cm和14.38 cm;东西向形变不对称性明显,断层上盘西北部向东水平运动,最大位移为18.37 cm,下盘东南部向西运动,最大位移不足8 cm。将南北、东西向形变与6个全球导航卫星系统（global navigation satellite system,GNSS）台站观测数据进行比较,两者一致性较好且均方根误差较小,分别为1.44 cm和1.77 cm,表明联合升、降轨InSAR观测和地震位错模型约束构建同震三维形变场方法具有较高可行性,显著降低了大地测量数据不足、InSAR观测对南北向形变不敏感等问题的影响。     

基于地壳分层的地震断层同震变形反演分析

顾及地球结构的非均匀性、分层和球体特性,发展了顾及地壳分层结构的地震断层同震位错模型的反演模式,建立了非线性反演方程,对主要的两类地震断层分别进行了模拟反演分析.反演结果显示,该反演模式可很好地将断层位错参数反演出来,利用该模式可较好地解决实际地震断层同震变形的反演问题.     

联合地震位错模型和ESISTEM方法提取地震同震三维形变场

高精度同震三维形变场对于研究地震变形模式、震源机制等具有重要意义。设计并实现了一种联合地震位错模型和扩展融合大地测量、卫星形变观测、应变张量估计(extended simultaneous and integrated strain tensor estimation from geodetic and satellite deformation measurements, ESISTEM)方法的新方法,以2021年Ms 6.4漾濞地震为例,利用哨兵1号A、B星(Sentinel-1A/B)升、降轨影像获得该地震的合成孔径雷达干涉测量形变场,利用地震位错模型正演得到的南北向形变分量进行约束,成功提取了该地震完整的同震三维形变场及应变场。结果表明,漾濞地震断层西南侧主要向西、向北运动,最大形变分别为4.8 cm、9.5 cm;东北侧主要向东、向南运动,最大形变分别为7.4 cm、4.6 cm;垂直向抬升、沉降的最大值分别为3.6 cm、3.4 cm;发震断层以右旋走滑运动为主,兼有少量正断分量;发震断层区域受到显著的膨胀、剪切和旋转作用。