利用PALSAR数据获取拉奎拉地震同震形变场

2009年4月6日发生的意大利拉奎拉Mw6.3地震造成了地表的严重破坏,为了解地震引起的地表变化情况,本文利用ALOS/PALSAR影像数据,采用二轨差分干涉测量技术提取了该次地震的地表同震形变场。通过分析可知:地表在西北-东南方向上发生错动,形变主要发生在40×30 km的范围内,西南部出现了沿视线方向的下降,东北部出现了沿视线方向的抬升,最大值分别为0.337 m和0.122 m。结果表明,获得的同震形变场与地震地质调查的结果一致。     

基于ALOS PALSAR数据进行雷达差分干涉提取玉树地震同震地表形变场

本文针对2010年4月14日玉树地震引起的地表形变,使用日本ALOS卫星PALSAR L波段雷达影像数据,应用两轨雷达差分干涉(DInSAR)处理得到了以玉树为中心11 000km2范围内的同震形变场,空间分辨率为8m,并在此基础上对玉树地震的震源机制和发震机理进行了分析。研究结果表明L波段雷达数据适合在地形起伏较大的地区进行DInSAR形变探测。该同震形变场信息可为玉树地震的同震形变反演提供参考数据。该研究进一步证实DInSAR技术在大规模地表形变探测和地学研究领域具有广阔的应用前景。     

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

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

PSInSAR方法探测意大利拉奎拉地震形变过程分析

采用双轨处理模式对8景ENVISAT降轨数据进行永久散射体(PS)时间序列处理,获得了2009年4月6日意大利拉奎拉(L’Aquila)Mw 6.3级地震区域的PS时序差分干涉图和形变场。结合Delaunay三角形剖分算法对破裂区的形变过程进行数值分析,结果表明:①本次地震变形最强烈、地表破裂发生的区域是一个地势相对较低的椭圆形洼地。破裂主要沿震中的东南方向传播,面积约22×14km2,方向约135°。②地震产生的形变量集中在破裂区,主要在震时及震后形成,位移场视线向下沉量达150mm。③破裂区位移场数值分析表明,本文提出的数据处理策略较完整地揭示了此次地震位移场变化的全过程,在8景雷达数据所跨时间段内,震前2008年4月与2009年2月两景资料显示,研究区位移场变化不明显,2009年2月至震后的2009年4月的图像显示,震中区形变场快速变化,沉降量达130mm;6～7月,形变场的沉降速度明显减缓,8月30日最后一景资料显示,沉降速度出现加速,破裂区总下沉量达150mm。     

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

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

融合DInSAR和MAI技术反演九寨沟地震三维形变场

差分雷达干涉测量(DInSAR)在地震形变监测中具有重要作用,但DInSAR仅能获取雷达视线向的形变量,无法准确反映地表三维形变.为了获取真实的地表形变,文中采用一种联合DInSAR和多孔径干涉技术(MAI)的地震三维形变场反演方法,以2017年九寨沟Ms7.0级地震为例,融合升降轨Sentinel-1A数据,联合解算...     

九寨沟Ms7.0级地震三维同震形变场提取

针对D-InSAR地震形变监测中存在LOS(line of sight)视线向模糊的问题,该文构建了一种融合升降轨不同视线向干涉测量数据获取三维同震形变场的方法。以2017年九寨沟M_S7.0级地震为例,基于研究区的Sentine1_1A数据和SRTM(1″)DEM数据,采用二轨差分的D-InSAR技术,融合升降轨LOS向以及自定义视线向的干涉测量数据进行联合解算,获取了研究区三维同震形变场,解决了单一轨道雷达LOS提取结果不能准确反映地表三维形变场的问题。实验表明,此次地震造成了震中附近一定范围地表的沉降和东南向的滑移,并通过震区形变剖面图和形变等值线图,分析了形变场的空间分布以及此次地震所造成地表断裂的位置,认为此次地震与塔藏断裂、虎牙断裂存在关联性,其运动形式为主动盘逆冲的走滑型地震。     

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

基于PALSAR数据的2010年Darfield地震地表同震形变场提取

2010年9月3日（UTC时间）,新西兰南部小镇Darfield发生了强度为Mw7.1的地震,其主震对地表造成了严重破坏。为了对由地震引起的地表形变进行研究,本文利用二轨DInSAR技术,通过2景ALOS PALSAR影像的升轨数据,成功提取了该地震的地表同震形变场。从对实验结果的分析可知,该地震的断层为东西走向,断层北盘地表呈现沿雷达视线方向的下降,最大量约为1.5米;断层南盘则发生了沿雷达视线方向的抬升,最大量超过了2米。     

DInSAR技术对地震同震形变场的研究

合成孔径雷达差分干涉测量（DInSAR）可用于监测厘米级或更微小的地表形变,以揭示许多物理现象,如地震形变、火山运动、大气变化、冰川漂移、地面沉降以及山体滑坡等。DInSAR作为一种新型的空间对地观测技术,具有不受时间和空间的限制、对地物具有一定的穿透性等传统测量所不可比拟的优势,已得到较为广泛的应用。简要介绍DInSAR技术的基本原理及其处理流程,以Barn地震为例提取Bam地震的同震形变场。     

地震破裂模型约束的中国阿里地震三维形变场

利用环境卫星（Environmental Satellite,Envisat）的升、降轨数据在地震破裂模型约束下获取2007年阿里地震的高质量同震地表三维形变场。首先,对合成孔径雷达（synthetic aperture radar,SAR）影像进行差分干涉处理,得到地震造成的视线向同震地表形变场;然后,以地震破裂模型为约束条件,采用赫尔默特方差分量估计法来解算阿里地震的高质量地表三维同震形变场。结果表明,震中区域最大下沉达约4.7 cm,东西向位移较小,南北向呈挤压趋势。总体上,三维形变的特征表明阿里地震是一个以正倾滑为主兼少量右旋走滑运动的地震事件。     

升降轨道ASAR雷达干涉揭示的巴姆地震（Mw6．5）3D同震形变场

通过对覆盖同一地区的升降轨道ASAR雷达数据进行两路差分干涉（D-INSAR）处理，得到了2003年12月26日伊朗巴姆（Bam）地震（Mw6．5）在两种成像几何状态下的视线向（LOS）同震形变场。利用Okada弹性形变模型对形变场进行正演计算，得到地震断层的几何参数及形变的南北向分量。对得到的干涉条纹图进行相位解绕处理，综合分析两种不同的雷达成像几何关系，并利用模型估计得到的南北向分量，计算了该地震造成的3D同震形变场，结果显示了典型的单断层右旋走滑活动特征。研究表明，合并升降轨道雷达数据，能够为形变模型提供有效约束，用一条简化的断层模型就可以解释升降两条轨能的干涉形变场，更为复杂的双断层模型是不必要的。     

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

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

一种新的SAR像素偏移量估计流程及其在同震形变监测中的应用

提出了一种无需外部数据的新同震形变偏移量计算流程。以2003年伊朗巴姆里氏6.6级地震为例,利用两景震前无形变数据证明了新方法的精度较传统方法提高了0.224个像素。利用震前震后两景数据计算了方位向同震形变场,与Okada模型模拟得到的方位向同震形变进行了比较。     

DInSAR技术监测玉树地震同震形变场的研究

近二十年来,合成孔径雷达干涉测量技术（DInSAR）作为一种监测地表形变的有效技术得到广泛应用,其视线向精度可达厘米级甚至毫米级。该技术尤其是对监测快速、激烈的地表形变更为有效,如地震、火山等。利用DInSAR技术,选用日本ALOS卫星PALSAR数据,获取了2010—04—14玉树Mw6．9地震的同震形变场。结果表明：地表至少发生过3次地表破裂,沿断层走向的形变分布范围远远大于垂直断层分布范围,形变分布特征以左旋走滑为主。     

2021年云南省漾濞Ms6.4地表二维形变分析

2021年5月21日云南省漾濞县西侧发生Ms6.4地震，造成较大人员伤亡和经济财产损失。采用D-InSAR技术对云南省漾濞的升降轨Sentinel-1A数据分别进行解算得到Ms6.4地震同震形变场。降轨结果显示在震中的东-南向有两个明显的隆升和沉降区域，雷达视线向最大沉降量为-7 cm,最大隆升量为8 cm。升轨结果显示雷达视线向最大沉降量为-7 cm,最大隆升量为7 cm。根据震源机制解和同震滑动分布模拟计算此次地震导致的周边断裂库仑应力变化情况，联合升降轨结果解算本次地震的垂直向和东-西向形变场，并分别分析形变场的形变方向和形变量级，可以为研究地震震源参数和同震滑动模型精细化反演提供更可靠的约束条件。     

2017年伊朗Mw7.4地震三维同震形变场及滑动分布

针对如何准确分析2017年11月12日的伊朗Mw7.4地震断层属性的问题,该文从三维同震形变场与断层滑动分布角度出发,采用Sentinel-1A升降轨数据获取该地区三维同震形变场,并基于弹性半空间矩形位错模型,结合最速下降法反演得到其发震断层参数、模拟三维形变场。结果显示,地震东西向和南北向最大形变分别为68.8、43.8 cm,竖直方向上最大沉降和抬升分别达到34.4、96.7 cm,均符合逆冲断层的运动特征。地震主要形变特征为抬升,地震矩1.66×10²⁰ N·m,矩震级Mw7.4,断层倾角16°,走向351°。结果表明,Okada模型解算所得地震三维形变场与观测结果总体特征相吻合,并且该地震为西南部板块挤压东北部板块而引发,属于逆冲为主带有右旋走滑的逆冲断层,其与地震学结果一致。     

地震同震滑动分布反演的总体最小二乘方法

同震滑动分布参数与地表形变间的线性关系依赖于格林函数矩阵的构造,格林函数矩阵元素与破裂面位置、几何参数、破裂方式及位错模型假设等因素有关。本文尝试考虑格林函数矩阵元素的误差来补偿上述原因在一定程度上对反演参数的影响,采用同时顾及系数矩阵(格林函数矩阵)和观测向量两者误差的总体最小二乘方法反演同震滑动分布。首先确定了系数矩阵元素和观测向量的协因数矩阵,考虑到格林函数矩阵的病态性(秩亏),借助拉普拉斯二阶平滑得到正则化矩阵,采用总体最小二乘正则化法反演同震滑动分布。并对2009年意大利中部拉奎拉(L’Aquila)Mw6.3级地震实例进行同震滑动分布反演研究。结果表明,拉奎拉地震的走向为144.37°,倾角为59.06°,滑动分布的最大滑动量为0.95m,平均滑动角为-96.4°,主要滑动深度为4~15km的范围,地震矩为3.63×10~(18)N·m,对应的矩震级为Mw6.34。总体最小二乘与最小二乘法的滑动分布解存在一定差别,但差别的量级在10-4以内。     

地震同震滑动分布反演平滑因子的确定

针对大地测量地震同震滑动分布反演中平滑因子（又称正则化参数）的确定问题,一般采用模型粗糙度与数据拟合残差之间的折中曲线来确定（为便于区分,本文将该方法称为"L曲线"）。本文在L曲线的基础上提出一种确定平滑因子的新方法——折中相交曲线法。模拟试验反演结果表明,利用折中相交曲线法确定的平滑因子反演地震滑动分布各参数精度要优于L曲线法。利用折中相交曲线法确定平滑因子反演拉奎拉与台湾美浓实际地震,并与L曲线法反演结果进行对比分析。分析结果表明,利用折中相交曲线法确定的平滑因子反演拉奎拉与台湾美浓实际同震滑动分布各参数结果均在国内外其他学者研究的范围内,并且利用折中相交曲线法确定平滑因子较L曲线法具有计算效率高、无需依赖数据拟合度、确定平滑因子大小更为合适等优点。     

基于升降轨ASAR的于田Ms 7.3级地震同震形变场信息提取与分析

升降轨干涉测量可以更好地反映地震同震形变场特征。利用欧空局ASAR数据通过二轨法差分干涉测量提取了新疆于田县Ms 7.3级地震升降轨同震形变场信息。结果显示:升轨同震形变场最大视线向(LOS)隆升形变量约+13.3 cm,沉降形变量约-83.9 cm;降轨同震形变场最大LOS向隆升形变量约+36.5 cm,沉降形变量约-66.5 cm。于田地震以NNE向正断层破裂为主,并伴随左旋走滑运动,西北盘为正断层破裂的上盘(沉降盘),东南盘为正断层破裂的下盘(隆升盘)。升降轨同震形变场存在一定差异,但其变化趋势与特征非常相似,其差异主要是由于两种不同观测模式所造成的。