基于InSAR技术的张家口定点形变观测场地形变监测

针对张家口定点形变观测场地伸缩仪NS向观测曲线在2020年6月至12月出现加速上升的异常变化,以张家口定点形变观测场地为研究区域,基于Sentinel-1 SAR降轨数据并采用小基线集技术对场地进行形变监测,并获取了2018—2020年该区垂向年平均形变速率和形变时间序列。研究结果表明：张家口断裂两侧地表在垂向的相对运动速率小于1 mm/a,断层活动性较弱;张家口断裂在NNW-SSE向拉张应力作用下,使得断层南盘地表相对北盘出现较为明显的沉降运动,导致张家口伸缩仪NS向观测出现张性加速趋势变化异常现象。     

基于MSBAS InSAR技术的沧州市地表形变监测与分析

针对沧州市城市地表形变监测问题,采用SBAS InSAR技术分别处理Sentinel-1A升轨和Sentinel-1B降轨SAR影像,对比分析升、降轨SAR影像监测的沧州市2020年7月至2021年4月地表形变特征;利用MSBAS InSAR技术联合处理升、降轨SAR影像,将沧州市地表形变进行垂向和东西向的二维形变分解,获取沧州市二维形变信息.结果表明,升、降轨SAR影像监测的地表形变位置和分布基本一致,但由于LOS向模糊问题,单独利用升轨或降轨SAR影像监测到的形变区域有所偏移,且无法分离出地表的东西向形变,影响监测结果的精度,MSBAS InSAR技术可以有效解决这一问题,获取更准确的二维形变信息.     

基于多卫星平台永久散射体雷达干涉提取三维地表形变速度场

永久散射体雷达干涉(PSI)技术及其应用于区域地表形变监测已成为雷达遥感领域的研究热点之一.使用单一卫星平台所获取的单侧视SAR影像时间序列进行PSI分析,仅能获取沿雷达视线(LOS)方向的一维地表位移信息.本文提出了基于多平台永久散射体雷达干涉提取三维地表形变速度场的模型与算法,其基本策略是:首先针对每一卫星平台的SAR影像时间序列进行PSI分析,并计算各地面目标沿LOS向的位移速度值,然后联合各平台所对应的LOS向位移速度值进行建模,并基于最小二乘方法解算各地面目标的三维位移速度分量.实验选取天津市西北部作为测试区,使用2007—2010年所获取的39幅TerraSAR-X影像、23幅ENVISAT ASAR影像和16幅ALOS PALSAR影像进行分析,经联合解算得到了该测试区域的垂直位移速度场以及南北向和东西向水平位移速度分量.与地面水准和已有GPS观测结果对比分析表明:多平台PSI的垂直位移速度场精度可达毫米级,而其水平位移速度分量与已有GPS结果基本一致.多平台PSI分析无需引入任何外部形变参考信息,便可以实现形变场的偏差校准和三维形变场的恢复.     

融合升降轨SAR干涉相位和幅度信息揭示地表三维形变场的研究

针对InSAR技术只能监测地表在雷达视线方向上的形变的缺点, 详细介绍了融合升降轨SAR干涉相位和幅度信息监测地表三维形变的具体方法. 该方法主要分为3个步骤: (1) 利用升降轨SAR干涉相位监测地表在2个不同雷达视线方向上的形变; (2) 利用升降轨SAR幅度信息获取地表在2个不同方位向上的形变; (3) 采用最小二乘准则和Helmert方差分量估计融合上述4个不同方向的形变, 从而估计地表三维形变场. 以2003年伊朗Bam地震为例, 应用该方法成功地揭示了该地震引起的地表三维形变场, 结果显示Bam地区北部出现了明显地表下沉和沿近西南方向的水平运动, 而南部则出现地表隆起和沿近东南方向的水平运动. 上下、南北和东西3个方向上的地表形变场都很好地吻合了地震断层所在的位置, 最大形变量分别达22, 40和30 cm. 最后, 将此三维形变场与利用Okada模型模拟的三维形变场进行了比较, 证明了该方法能够得到可靠且精度较高的地表三维形变场.     

通过InSAR与GPS数据融合获取汶川地震同震三维形变场

针对In SAR与GPS两类独立的观测数据,利用广义测量平差理论中方差分量估计法,合理分配权重,有效地融合了这两类观测数据,从而估算出了地表在三维方向上的形变场。以四川汶川地震为例,利用In SAR干涉测量结果和一定数量的GPS观测值,通过该方法获取了地震断层两侧高相干区域上三维形变场,清晰地显示了汶川地震的逆冲和右旋走滑分量的位置分布整体特征。结果表明,在EW、SN和UD方向上地表形变量与GPS结果能够很好地吻合,融合结果在三维方向上的均方根误差均不超过5cm,证明了该方法能够得到精度较高的地震同震三维形变场,也揭示了利用方差分量估计法对相互独立数据之间有效融合的可行性。     

王艳，廖明生，李德仁等.利用长时间序列相干目标获取地面沉降场.地球物理学报，2007

方法，相干目标分析是利用时间序列多景影像获取地表稳定反射体——高相干点目标的新方法.实验区选择上海市主城区约100 km2 区域，InSAR的实验数据采用覆盖1992～2000年间的25景ERS-1/2的单视复影像.在无先验沉降场模型的情况下，实验结果显示该实验区的地面沉降场在时间范围上跨越8年，平均观测时间采样率约为4个月，垂直形变精度优于5 mm，满足我国对城市地面沉降形变观测的要求.研究结果表明该方法可以用于今后大时间/空间尺度上的高精度地表形变观测.     

王艳，廖明生，李德仁等.利用长时间序列相干目标获取地面沉降场.地球物理学报，2007

方法，相干目标分析是利用时间序列多景影像获取地表稳定反射体——高相干点目标的新方法.实验区选择上海市主城区约100 km2 区域，InSAR的实验数据采用覆盖1992～2000年间的25景ERS-1/2的单视复影像.在无先验沉降场模型的情况下，实验结果显示该实验区的地面沉降场在时间范围上跨越8年，平均观测时间采样率约为4个月，垂直形变精度优于5 mm，满足我国对城市地面沉降形变观测的要求.研究结果表明该方法可以用于今后大时间/空间尺度上的高精度地表形变观测.     

基于时间序列SAR影像分析方法的三峡大坝稳定性监测

本文提出了结合永久散射体和准永久散射体干涉测量技术的时间序列合成孔径雷达（SAR）影像几何信息提取方案,利用2003～2008年获取的40景SAR影像量测并分析了三峡大坝运行时的稳定性问题及其周边区域的形变规律.我们发现在数据获取期间,三峡大坝左段的时效沉降基本停止,大坝形变受水位的影响较大,季节性变化也有一定体现.该实验结果和已发表的该坝段形变监测结果十分吻合.在数据处理结果中我们还在大坝上游秭归县附近堤坝发现一处形变异常.研究结果表明了时间序列SAR影像信息提取技术在大坝稳定性监测及坝区形变参数获取中的应用潜力.     

2021年青海玛多MW7.3地震InSAR的同震形变场、断层滑动分布及其对周边区域的应力扰动

文中基于D-InSAR技术,利用欧空局C波段升降轨哨兵SAR数据,获取了2021年5月22日青海玛多MW7.3地震的InSAR同震形变场,并对同震形变的空间特征、量级和断层破裂的分段性进行了分析.哨兵卫星的高质量观测数据清晰地描绘了玛多地震的地表破裂迹线,地表破裂长度约210km.为了进一步认识玛多地震断层深部的同震滑...     

2001年MW7.8昆仑山地震震后形变观测、机制与青藏高原中北部岩石圈流变:认识与挑战

2001年M_W7.8昆仑山地震是近半个世纪以来青藏高原发生的最大震级地震。同震破裂产生的巨大应力扰动驱动控制着显著震后形变。二十年尺度的大地测量数据记录了地震后长时间、大范围、时空依赖的震后形变演化过程及差异,揭示了昆仑山地震破裂段复杂的断层分段震后运动学特征、分段摩擦性质差异和青藏高原中北部岩石圈流变性质/结构横向各向异性。本文简要回顾昆仑山地震后基于二十年尺度时序InSAR和GPS的震后形变观测方法和时空特征,特别是时空密集的InSAR观测,是该构造区震后GPS观测的重要补充及其不可替代的观测手段;总结大范围震后形变模拟方法及其揭示的震后运动过程、多种动力学机制及其关系。最后总结提出昆仑山地震震后形变20年研究取得的科学认识及尚待深入研究的科学问题,一方面要持续性地对东昆仑断裂带大范围地表形变进行观测研究;另一方面,要不断更新震后形变模型,进一步深化对该断裂带地震周期形变、区域构造对周期形变控制作用、复杂断层运动时空演化机制的认识。      

基于GPS观测研究2010年青海玉树MS7.1地震震后地壳形变特征及其机制

2010年4月14日青海玉树M_S7.1地震发生在青藏高原东南部甘孜-玉树地震带,在震后7~10天内,我们快速建立了由15个GPS测站组成的跨地震破裂带观测剖面,包括1个连续站,3个半连续站和11个流动站,对所有站进行了240多天的观测,获取了该次地震的震后形变时空特征.采用欧拉矢量和位错模型解算了背景速度场,并从GPS观测的形变场中扣除该分量.采用分层黏弹性位错模型计算余震引起的地表形变,结果表明余震对部分测站的位移造成不可忽视的影响.采用对数模型拟合位移时间序列,表明特征衰减时间为6.7±1.2天.利用最速下降法反演震后余滑时空分布,反演结果表明震后断层活动以左旋滑动为主,断层南盘具有少量的抬升.在空间分布上,余滑主要位于同震破裂区的两侧,西北侧的余滑几乎达到地表,而东南区的余滑基本在同震破裂区的下方,余滑最大的区域位于结古镇东南下方10~20 km的深度范围.随着震后离逝时间的增加,2个余滑区在空间上保持不变,余滑区的面积逐渐扩大.余滑的矩释放为（1.5~5.1）×10¹⁸Nm,相当于1个M_W6.1~6.4地震释放的能量.分层岩石圈黏弹性模型计算的地壳孔隙弹性反弹形变与地表观测值相差较大,不能解释观测到的震后变形.采用麦克斯维尔流变体模型计算下地壳和上地幔松弛引起的地表形变,显示出其对地表形变的贡献较小.GPS观测得到的震后形变所具有的快速衰减特征,以及余滑模型能够较好地拟合GPS地表形变,表明2010年玉树M_S7.1地震后早期阶段的地壳形变主要是由余滑机制决定的.     

汶川地震SAR形变场电离层影响校正与断层地表破裂线自动提取

针对汶川地震PALSAR观测的方位向形变场存在严重电离层影响,导致准确提取断层地表破裂线困难的问题,本文使用一种基于频率域的电离层影响校正方法,从方位向形变场频域数据中提取电离层信号对应的频域信息,通过建模获得空间域电离层信号,并与原始形变场进行差分运算完成电离层影响的校正,还原了汶川地震断层近场沿SAR方位向的真实地表形变场.以校正后的方位向形变场作为数据源,采用非极大值抑制二维梯度算法,准确提取汶川地震北川—映秀断层和灌县—江油断层的地表破裂线.为验证提取结果的可靠性,引入野外地质调查数据进行对比分析,定量计算两条破裂线间的距离差异,结果表明基于SAR形变场自动提取的断层破裂线与野外地质调查数据具有较高的一致性,利用SAR观测形变场可作为调查地震断层地表破裂的重要补充手段.     

GPS与PS-InSAR联网监测的台湾屏东地区三维地表形变场

联合高精度的GPS水平位移观测和高密度的PS-InSAR雷达视线位移测量,实现地表三维形变的精确反演.本文在准确计算卫星轨道方位角基础上,使用GPS观测位移与星载雷达LOS方向形变的投影转换模型,将雷达LOS方向形变转换为垂直方向位移,并基于地面GPS与SAR影像PS目标联合构建形变监测网,采用参数平差算法估计区域地表形变场.以地质构造活动极其活跃的台湾岛及其西南屏东高雄地区为例,联合屏东地区48个GPS监测台站与雷达PS目标,监测该地区从1995-1999年间由于板块构造挤压运动和地下水抽取导致的三维地表形变.结果表明,该地区年均水平位移量为向西30~50 mm/a,高雄沿海地区发生明显的逆时针西偏南的逐渐增大的水平位移;垂直位移为屏东平原南部呈现-10 mm/a~-15 mm/a的地面沉降,而平原北部和高雄地区呈现约+5 mm/a~+10 mm/a的地面抬升.     

时序InSAR技术及其在京西北试验区微小形变的应用研究

时序InSAR技术是近些年发展起来的一种新的InSAR测量技术.其中PS-InSAR技术是通过分析高相干的永久散射体相位变化信息提取形变,受时间和空间基线的限制小,该方法需要假设形变特征为匀速,同时需要大量的SAR数据（25景）来估算大气效应贡献值;短基线技术（SBAS）可以降低几何去相干的影响,同时较短基线使DEM误差的敏感性大大降低,但该方法通过一定大小的窗体计算相干性,并据此来选择时间序列稳定点,窗体操作降低了被选像元的分辨率,但抑制了噪声相位的影响;干涉图叠加方法（interferogram stacking）是将多幅差分干涉图进行叠加,提高结果中形变信息对大气干扰的相对精度,可以获取大区域的连续形变场;CR-InSAR技术将人工角反射器作为观测对象,通过建立CR点上的形变模型,实现对重点目标如滑坡、水库和断层的形变监测.     

基于InSAR数据分析四川盆地南部长宁页岩气区块地表形变场基本特征

近年来随着我国页岩气大规模开采,四川盆地南部活动构造相对稳定的地区出现了一系列微震和有感地震,甚至是破坏性地震。这些地震是否为工业开采所诱发,目前已有研究从时空相关性给出了一些统计推断,本文则从形变观测角度分析页岩气开采能否产生可以检测到的地面形变,以揭示形变信息与页岩气开采的关系,尝试为页岩气开采提供有效的监测手段。基于长波ALOS-2卫星雷达数据对长宁页岩气区块近两三年内的InSAR地表形变展开探测,检测页岩气发活跃时段内的形变时间序列信息。结果显示：考虑到不同观测砽技縂术的误差水平和观测角度差异,两种卫星数据均反映了一致的地表形变分布,且形变场与页岩气开采井的空间分布有很好的对应关系;压裂注液过程会造成地表快速隆升,生产过程中随着流体扩散地表会出现沉降和水平运动,初步揭示出页岩气生产过程中地面形变的非稳态变形特征。这表明在四川盆地南部复杂的形变观测条件下,InSAR技术是页岩气开采有效的监测手段,能够弥补地震学观测的不足。     

1990年青海共和7.0级地震震后垂直形变研究

根据1990年青海共和地震震后地表垂直形变,通过模型拟合得到了支配共和地区震后形变场时空演化的形变源及其力学机制.分析穿过断层的震前1期和震后6期水准数据,结果表明震后垂直形变具有以下特征:①震后震区上盘继续发生继承性的大幅度上升,其中震后头一年上升速率最大;②震后上升区范围显著,范围随时间变化不大,但较同震形变上升区范围增大;③震后相邻测站高差观测值的时间序列明显具有对数衰减特征或指数衰减特征,衰减特征时间分别为0.165年和1.344年.本文还发展了一个利用水准数据与连续介质位错模型研究震后形变机制的新方法.该方法用相邻水准点之间的原始高差观测值而非它们相对参照点的积分值来约束连续介质位错模型,可以有效减少误差累积带来的偏差并充分利用观测数据.利用这一方法的初步分析结果表明,断层震后滑移和介质黏弹性松弛共同导致了共和地震震后形变.前者表现为发生在断层面及其延伸部分的滑移,特别是位于主破裂上方沉积层内的滑移;后者则表现为下地壳与上地幔内的黏弹性松弛,黏滞系数为1020Pa.s量级.     

综合InSAR和应变张量估计2016年MW7.0熊本地震同震三维形变场

利用合成孔径雷达(Synthetic Aperture Radar,SAR)影像提取与地质活动相关的三维地表形变场,对深入理解地质灾害的形成机制及其潜在灾害风险评估非常重要.目前,利用SAR影像的同震三维形变场提取主要利用单个像素点的多次观测构建观测方程,然后基于加权最小二乘(Weighted Least Squares,WLS)方法分解从而获得同震三维形变场,因此该方法缺乏对相邻像素点空间相关性的约束.考虑相邻同震位移点的应力连续性,研究学者提出了顾及大地测量应变张量和卫星形变观测的SAR同震三维形变场方法(Extended Simultaneous and Integrated Strain Tensor Estimation from geodetic and satellite deformation Measurements,ESISTEM).本文以2016年MW7.0熊本地震为例,收集了覆盖此次地震的ALOS-2卫星升降轨影像,利用传统差分InSAR(DInSAR)方法和子孔径雷达干涉测量(Multiple Aperture InSAR,MAI)方法分别对升降轨SAR影像对进行处理,得到视线向(LOS)形变和方位向形变,最后利用ESISTEM方法获取此次地震的三维同震形变场.此外,利用GPS和野外考察观测对本文的三维形变场结果进行结果精度分析.研究结果表明,与传统WLS方法相比,ESISTEM方法不仅能有效抑制奇异像素点对形变结果的干扰,同时对近断层的失相干信号能进行较好的恢复,更有助于解释地表破裂区的地震形变特征和掌握地震发生机制.本文确定的三维同震形变场结果显示主形变区发生在Futagawa断层中部和Hinagu断层最北端,最大水平位移为2m,抬升为0.55m.断层破裂以NE-SW走向的右旋走滑为主兼有部分正断成分.应变张量分析表明发震断层处受到了明显的收缩力和剪切力的作用.     

雷达干涉PS网络的基线识别与解算方法

时序雷达干涉图中的永久散射体(PS)可看作“天然GPS点”, 以构成网络用于监测长期的地表形变. 本文提出采用邻接矩阵拓扑模型对基于Delaunay剖分算法生成的PS网络进行基线识别, 并采用时序相干最大化算法求解PS基线的线性形变速度增量和高程误差增量. 该数据模型和计算方法被应用于探测香港地区2006~2007年间的区域地表沉降. 实验研究采用由Envisat卫星ASAR传感器对该地区成像所获取的时序SAR影像作为数据源, 并联合该地区12个GPS连续运行参考站的观测数据予以大气修正和地面控制. 实验结果表明, 该模型和方法应用于地表形变测量是有效的和可靠的, PS网络方法探测地面沉降的精度约为±2.0 mm/a.     

基于SBAS-InSAR技术的延安新城地表形变特征

延安GNSS基准站2015-2017年出现异常变化,在32幅哨兵-1雷达影像基础上,利用小基线集合成孔径雷达干涉测量（SBAS-InSAR）技术,对延安新城2014年10月-2017年2月地表形变进行监测,提取时间序列与平均地表形变速率。结果表明：GNSS基准站的异常形变可能与延安新城土方回填相关,地表形变时间序列沿顺时针的变化趋势与延安新城的扩张方向基本一致。与传统GNSS异常核实方法相比,SBAS-InSAR技术可快速确定异常范围和量值大小。通过高精度地表形变时间序列,能够辅助查找异常发生原因,有利于推进防灾减灾工作。     

2008年仲巴地震形变及其揭示的构造运动

2008年8月25日拉萨块体中部的仲巴县发生了MW6. 7地震,之后发生了MW6. 0余震。然而,针对此次地震,不同机构给出的震源位置和发震断层参数存在一定差异,同时对其所揭示的构造意义也需要更深刻的认识。文中利用不同波长、不同入射倾角的升降轨ENVISAT ASAR和ALOS PALSAR卫星数据,基于In SAR技术获取了该地震的8幅同震形变场图像,得到此次地震引起的地表形变长约50km,且分为南、北2瓣。为了更好地反演主震的破裂滑动分布,在数据降采样过程中,设置近场区域的采样间隔远小于远场区域,以保证主要形变区域的观测数据特征和采样密度,并且扣除了余震形变的影响。基于Okada均匀弹性半空间位错模型和SDM方法反演的结果表明,此次地震的发震断层为E倾的帕龙错断裂,破裂滑动主要为正断倾滑兼部分左旋走滑性质,南段最大滑动量为1. 15m,北段最大滑动量为0. 52m。为了进一步理解此次地震的构造背景,文中利用SBAS-In SAR技术获得了震后仲巴地区2008—2010年的形变时间序列,发现形变主要是由剥蚀导致的沉降区堆积所引起的,而在断裂两侧并不明显。结合1991—2015年的GPS数据分析,认为拉萨块体内部EW向的伸展作用是明显且不均匀的,地堑附近伸展量逐渐减小。正断型地震在高原隆升过程中可能起着重要的调节作用。