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

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

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

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

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。反演结果显示,同震位错以逆冲运动为主兼具少量右旋走滑运动,发震断层北倾,倾角约...     

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

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

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

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

InSAR数据约束下2016年和2022年青海门源地震震源参数及其滑动分布

中国青海省门源县于2016年和2022年分别发生了Mw 5.9和Mw 6.7地震,相距不足40 km。利用欧洲空间局Sentinel-1A升降轨雷达影像,采用合成孔径雷达干涉测量(interferometric synthetic aperture radar, InSAR)技术分别获取两次地震的同震地表形变场,进而利用弹性半空间的位错模型确定上述事件的震源参数,基于分布式滑动模型反演确定两次地震断层面上的滑动分布,并探讨2016年门源地震对2022年门源地震的发震影响及触发机制。结果表明,2016年门源地震为逆冲型地震,并未破裂到地表,升、降轨同震形变场沿视线向的最大形变量分别为6.7 cm和7.0 cm,断层的最大滑动量为0.53 m,主要集中在地下4~12 km区域滑动。2022年门源地震同震形变场沿NWW-SEE向破裂,降轨影像最大视线向地表形变量为78 cm,断层的最大滑动值达到3.5 m,处于地下4 km左右,断层滑动分布模型揭示此次地震为左旋走滑型地震;结合冷龙岭断裂的运动性质和几何特征,可初步判定发震断层主要为冷龙岭断裂的西段、且极有可能破裂到了其西北端西侧的托莱山断裂。静态库仑应力触发关系显示,2016年门源地震对2022年门源地震的发生有一定的促进作用。     

西藏改则地震升降轨同震形变场特征分析与破裂模式

通过升降轨ASAR数据的获取与差分干涉处理,获取了西藏改则地震的双视线向同震形变场信息。对升降轨同震形变场的特征分析表明,改则地震主、余震均表现为典型的正断裂破裂模式,并先后形成了东、西两条规模不同的走向NE、倾向NW的正断层破裂带(可能未出露地表),以及西北盘的东、西两个沉降形变中心。从干涉条纹的切割关系判断,东沉降形变中心可能受2008年1月9日Ms 6.9级主震的控制,西沉降形变中心可能受2008年1月16日Ms 6.0级余震的控制。Ms 6.9级主震形成东南盘隆升与西北盘沉降的总体格局;Ms 6.0级余震使西北盘余震震中及以西部位进一步沉降,造成西沉降中心沉降形变量大于东沉降中心。升降轨同震形变场的获取为进一步的震源参数模拟与三维形变场解算提供了更好的约束条件与信息。     

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

基于D-InSAR的长宁地震形变场提取与模拟

2019年6月19日四川长宁县发生Ms6.0地震,提取地震同震形变场,反演震源机制,对于地震破裂分析、指导救灾具有重要意义。本文利用两轨差分干涉处理覆盖长宁地震影响区域的两景Sentinel-1A影像,在对D-InSAR关键技术和影像特征研究的基础上,配置处理方法和参数,提取了此次地震的同震形变场,采用Okada模型和正向建模对发震断层的几何参数及形变场进行了反演和模拟。结果表明,长宁地震形变场呈西北—东南走向的次级断层控制的不规则椭圆形,断层的两侧区域特征差异明显,断层左下侧为沉降区,右上侧为隆升区,两者的最大视线向形变分别为8、6 cm。断层的运动主要以左旋走滑为主,平均滑动距离约0.38 m,平均滑动角约55°,正向模拟的形变场与观测结果相符,这表明观测结果较可靠,同时也提高了低相干形变区的观测精度。     

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模型解算所得地震三维形变场与观测结果总体特征相吻合,并且该地震为西南部板块挤压东北部板块而引发,属于逆冲为主带有右旋走滑的逆冲断层,其与地震学结果一致。     

联合地震位错模型和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;发震断层以右旋走滑运动为主,兼有少量正断分量;发震断层区域受到显著的膨胀、剪切和旋转作用。     

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

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

基于多源SAR数据的2022年门源Ms 6.9地震同震破裂模型反演研究

2022年1月8日青海省海北州门源县发生Ms 6.9地震,震中位于青藏高原东北缘祁连-海原断裂中段,属历史地震空区,基于多源合成孔径雷达（synthetic aperture radar,SAR）遥感数据研究该地震的破裂模式对理解青藏高原东北缘构造变形机制、应变释放过程以及地震危险性评估具有重要意义。首先利用Sentinel-1数据和合成孔径雷达差分干涉测量（differential interferometry synthetic aperture radar,D-InSAR）技术获取了门源地震的同震形变场,视线（line of sight,LOS）向形变场显示此次地震造成了约20 km长的地表破裂,最大形变约0.75 m;然后基于Sentinel-2卫星数据,利用光学影像配准和相关技术获取了本次地震的东西向同震形变场,最大同震位移达2.5 m;最后基于均匀弹性半无限位错模型,以LOS向形变场为约束反演了断层的滑动分布模型。结果显示,门源地震是一次典型的左旋走滑型地震,地震破裂主要集中在0~10 km深度范围,最大滑动量3.25 m,滑动角10.44°,对应深度4.89 km;反演给出的矩震量为1.07×1019 N·m,对应矩震级Mw 6.6。结合野外考察和地质资料,初步判定发震断裂为冷龙岭断裂,并引起托莱山断裂发生同震滑动。同震库仑应力结果显示,冷龙岭断裂东段和托莱山断裂西段应力状态为加载,未来具有发生强震的风险。     

唐山7.8级地震三维地表形变和断层运动特征

通过系统收集整理1976年唐山7.8级地震前后区域水准和边角网观测资料,在研究发展动态平差程序的基础上,解算得到了唐山地震震前和同震三维形变定量结果与误差分布,并反演得到了断层面同震滑动分布。结果显示,震前阶段丰润-唐山-唐海-昌黎-迁安一带整体隆升,隆升高值区位于唐山附近,震中区的水平形变量不大,唐山菱形地块相对于周围地块逆时针扭转;同震阶段唐山断裂以右旋错动为主,并呈现西升东降的垂向变形特征;反演结果表明多段断层参与了唐山7.8级地震的应变释放,主震发震断层的南段破裂主要发生在深度约6~18 km、长轴约50 km范围,北段破裂主要发生在深度约7~17 km、长轴约30 km范围,滦县地震断层也发生了少量滑动;主震断层最大水平位错为7.82 m、最大垂直位错为2.04 m,等效震级为Mw 7.58。挖掘得到的唐山震前和同震三维形变场定量结果为深入研究该地震的孕育发生机制及地震影响提供了珍贵资料。     

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数据对震源机制进行反演、模拟和分析,获取了高精度的震源参数,可以为分析地震危险性和断层破裂参数等提供数据支持。     

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

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

2008年大柴旦M_w6.3级地震的InSAR同震形变观测及断层参数反演

2008年11月10日,青海省大柴旦地区发生了Mw6.3级地震。本文利用EnviSat卫星升降轨SAR数据和差分干涉测量技术提取了同震形变场,基于均匀位错模型反演确定了地震断层参数,然后利用格网迭代搜索法确定了较优断层倾角,同时基于非均匀位错模型获得了精细滑动分布。结果表明,地震使得上盘区域沿降轨、升轨视线向分别产生最大约8.5cm、10cm的抬升;较优断层倾角为47.9°;地震滑动未延伸至地表,主要发生在地下8.2~23.7km深度范围内,最大和平均滑动量为0.5m和0.19m,平均滑动角为104.9°。反演的地震矩为3.74×1018 N·m,矩震级为Mw6.35。     

尼泊尔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级。