InSAR同震地表三维形变反演：一种顾及形变梯度的联合解算方法

同震地表三维形变场是对地震引起的地表形变在真实空间中最直观的描述,在地震形变研究中,同震地表三维形变场的恢复具有重要意义。针对目前地表三维形变解算方法存在的不足和缺陷,提出了一种顾及形变梯度的同震地表三维形变场联合解算方法。该方法利用形变梯度信息将合成孔径雷达干涉（interferometric synthetic aperture radar,InSAR）观测的形变区域分解成形变梯度近似相等的互不相交的子区域,然后根据子区域内观测值数量选择加权最小二乘法（weighted least square,WLS）或基于应力-应变模型（stress-strain model,SM）和方差分量估计（variance components estimation,VCE）的InSAR三维形变解算方法（SM-VCE）获取同震地表三维形变场。模拟实验和实际震例研究结果表明,相比于传统的WLS和SM-VCE方法,通过形变梯度信息能有效恢复地震破裂带的地表三维形变,得到更加完整和可靠的同震地表三维形变场。     

基于GIS的InSAR结果分析方法及在汶川Mw7.9级地震同震解释中的应用

利用GIS对InSAR结果进行分析可以更全面地解释地表形变机制。介绍了利用GIS对SAR影像进行裁剪、切割和拼接处理、剖面和等值分析的方法。在此基础上,以汶川Mw7.9级地震为例介绍了基于GIS的InSAR结果分析方法在地震形变机制解释方面的应用:①根据InSAR距离向偏移量和同震形变场,利用GIS技术提取了汶川Mw7.9级地震发震断层的地表迹线;②将基于位错模型正演得到的视线向形变场与InSAR同震形变场进行比较,得到了汶川地区InSAR同震形变场的上盘和整体改正数,进一步借助GIS进行修正得到了更为真实的InSAR同震形变场;③利用2D/3DInSAR同震形变场、干涉纹图和剖面分析结果对汶川Mw7.9级地震同震形变特征进行了具体分析和解释。     

联合地震位错模型和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观测对南北向形变不敏感等问题的影响。     

基于InSAR数据的西藏玛尼Ms7．9级地震的地壳不均匀性研究

基于InSAR观测结果,考虑到发震区域南、北两盘地壳介质的横向不均匀性,结合更为真实的地表断层位置,采用三维半解析粘弹性模型和遗传算法,反演了发震断层滑动分布及南北两盘的地壳介质参数(泊松比).反演结果显示,考虑介质横向不均匀性的模型能较好地拟合地表形变场;南北盘的横向介质差异性也可认为是玛尼地震同震形变的非对称性的影响因素之一.     

基于曲面断层结构的2007年所罗门群岛地震同震滑动分布反演

2007年4月1日,南太平洋岛国所罗门群岛的俯冲带地区发生了Mw 8.1级地震.利用日本宇宙航空研究开发机构(Japan Aerospace Exploration Agency,JAXA)的ALOS/PALSAR卫星雷达影像数据,采用二通差分干涉技术获得了覆盖所罗门群岛地震震区的同震地表形变场,根据四叉树方法对原始观测数据进行降采样,利用弹性半空间三角形位错模型进行断层几何参数反演;并结合断层深度剖分图的特征,构建更为合理的上下两个不同倾角连接的断层面模型,通过提出的断层剖分技术确定了最佳同震滑动分布模型.为了测试断层倾角的最佳估计,对不同断层倾角进行测试,对比合成孔径雷达干涉测量(interferometic synthetic aperture radar,InSAR)技术获取的同震观测数据和不同模型预测的InSAR干涉图,发现上下断层面倾角分别为10°和23°的模型更为合理.据此模型得到该断层同震滑动分布主要发生在距地表40 km深度范围内;最大同震滑移量为5.227 m,深度为2 km.     

利用SAR影像配准偏移量提取地表形变的方法与误差分析

单一的InSAR观测技术可提取地表沿雷达视线方向(LOS)上的一维位移,而利用SAR影像配准过程中的同名像素偏移量可提取地表沿雷达方位向(近南北向)与距离向(近东西向)的二维形变场,与LOS方向的一维形变形成优势互补。本文在分析SAR影像配准偏移量提取地表形变场方法的基础上,推导建立了雷达方位向与距离向形变提取的误差模型,探讨了该方法提取地表形变的主要误差源。以Bam地震ASAR影像和玉树地震PALSAR影像为数据源,分别开展同震形变场的提取与误差分析试验,结果表明,基于SAR像素配准偏移量提取同震形变场的精度主要受匹配窗口尺寸与过采样因子影响,形变提取误差随匹配计算窗口的增大而减小,形变提取精度随过采样因子的增大有适量提高,地形起伏效应在高差较大的SAR影像区域中表现较为显著。     

多视角InSAR数据解算2017两伊地震三维同震形变场

2017两伊地震是自1900年以来发生在扎格罗斯山脉的最大地震,为了研究此次地震引起的同震形变场,利用覆盖同一地区的3对Sentinel-1A升降轨数据分别进行两通差分DInSAR处理,得到了研究区3个视线向的地表同震形变场,通过直接解算法重建了研究区的三维同震形变场。试验表明：3种视角的升降轨视线向上升与沉降总体趋势基本一致;联合多个视角的观测结果可以实现三维形变场的重建;根据地表视线向和三维同震形变的特征以及地质构造背景推测了发震断层很有可能为扎格罗斯山前断层。     

InSAR数据约束下的2023年赫拉特地震序列发震断层探讨及其建筑物损毁评估

2023年10月7日阿富汗西部赫拉特省在不足1 h内接连发生4次Mw 5.5+的地震,称之为“2023年赫拉特地震序列”,此次地震序列是阿富汗境内过去20多年来遭遇伤亡最严重的地震事件,研究该地震序列的发震断层几何和快速分析损坏建筑物的分布状况对理解Herat断裂系统的构造机制、保障高效救援以及灾后规划重建等工作具有重要科学意义。基于欧洲空间局Sentinel-1升降轨合成孔径雷达影像,利用InSAR技术（interferometric synthetic aperture radar）获取2023年赫拉特地震序列的同震形变场,以升降轨InSAR观测为约束,反演确定发震断层几何和断层滑动分布,并对发震构造进行分析;基于多时相InSAR相干性变化探测方法分析并提取本次震后建筑物损毁代理图（building damage proxy map, BDPM）。结果表明,升降轨同震形变场均位于Herat断裂带和Siakhubulak断裂带之间,升降轨数据观测得到的最大视线向形变量分别约为32.3 cm和58.9 cm。反演结果显示,同震位错以逆冲运动为主兼具少量右旋走滑运动,发震断层北倾,倾角约...     

InSAR-MAI技术监测玉树地震2-D同震形变场研究

InSAR由于大范围、高精度、全天候等优势,被广泛应用。但InSAR局限于监测视线向形变,无法观测方位向形变。多孔径干涉测量技术(MAI)作为一种监测方位向形变的有效技术得到广泛应用,其精度可达厘米级。该技术对监测快速,激烈的地表形变更为有效,如冰川流速、同震形变等。文中利用MAI-InSAR技术,选用日本ALOS卫星PALSAR数据,监测玉树Mw=6.9地震的2-D同震形变场,获得有益的结论。     

融合InSAR和GNSS的三维形变监测：利用方差分量估计的改进SISTEM方法

利用干涉合成孔径雷达(interferometric synthetic aperture radar, InSAR)技术可以获取空间连续的地表形变监测结果,但仅能得到真实三维形变在雷达视线方向上的一维投影,降低了地表形变解译的可靠性。全球导航卫星系统(global navigation satellite system, GNSS)可提供高精度的三维形变结果,但其空间分辨率较低。因此可通过融合InSAR和GNSS来实现高精度高空间分辨率的三维地表形变监测。提出利用方差分量估计(variance component estimation, VCE)方法对经典的联合大地和卫星同步形变测量值的应变张量估计(simultaneous and integrated strain tensor estimation from geodetic and satellite deformation measurements,SISTEM)方法进行改进,发挥SISTEM方法可以考虑相邻点地表形变的空间相关性的优势,优化了函数模型中InSAR和GNSS观测的构成比例,引入VCE方法精确估计了各类观测值的后验方差,进而实现了高精度的三维地表形变测量。首先通过设计模拟实验对所提方法进行了验证,然后将该方法应用于2007年夏威夷基拉韦厄火山喷发案例中。相较于SISTEM方法,该方法在东西、南北、垂直向解算的形变精度均有大幅提升。