青海玉树地震差分干涉雷达同震形变测量

2010-04-14青海玉树发生7.1级地震后,作者利用震前和震后获取的日本ALOS卫星PALSAR遥感数据,开展了差分干涉雷达(D-InSAR)地震同震形变测量与分析。结果表明:玉树地震引起较大范围地表变形,地震变形沿玉树—甘孜断裂带向南东东方向扩展,在N33.7°,E96.81°附近达到最大形变量,D-InSAR监测到雷达视向上的最大形变量为35cm。地表形变特征对于评价玉树地震破坏程度、推断断层性质、研究地震形变和地震孕育特征具有重要的参考价值。     

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

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

Sentinel-1A数据的四川宜宾Ms6.0同震形变分析

为了获取2019年6月17日发生的四川宜宾Ms6.0地震引起的地表形变情况,该文利用欧空局宽幅模式的高分辨率新型Sentinel-1A卫星获取了此次地震的第一对同震干涉像对数据,使用D-InSAR技术获取宜宾市长宁县地区的同震形变场。结果显示,本次地震在震中西北方向分别形成了1个明显的沉降区和抬升区,在雷达视线方向上的最大沉降量为7.9 cm,最大抬升量为8.1 cm。通过与同一时间内的GPS高程测量形变量相比,D-InSAR解算的地表形变量与GPS监测点形变量基本一致,均不超过3 mm,表明了本文的D-InSAR形变解算结果的可靠性,体现了新型Sentinel-1A雷达卫星在地震形变监测领域有着很高的应用价值和潜力。     

基于 D-InSAR 技术的当雄地震形变场提取研究

随着D-InSAR技术在同震形变量测方面的优势日益突出,文中以2008年10月6日当雄Mw6.3级地震为对象,采用ENVISAT搭载的ASAR传感器c波段在图像模式下的5景level 0级影像,以90 m分辨率的SRTM DEM作为外部DEM,使用GAMMA软件采用二轨差分干涉测量的方法获取当雄地震同震形变场,并对其进行分析,确定震中位置为90.372°E,29.734°N,得到同震垂直形变约为33cm,形变量精度达到厘米级。     

利用Sentinel-1A数据监测地表形变案例分析

雷达卫星结合InSAR技术已广泛应用于高精度地表形变监测领域。本文选取2017年九寨沟地震为研究案例，利用Sentinel-1A地震前后的单视复数影像，基于D-InSAR技术获取该次地震的同震形变场。结果显示：震中西北侧表现出相对均匀的下沉现象，沉降漏斗区雷达视线向最大沉降量达25.1 cm；东南侧呈现不均匀抬升状态，地表破碎较为明显，最大抬升量为11.6 cm。研究表明基于Sentinel-1A数据的D-InSAR技术可以为地震形变场的定量分析提供一种快速有效的手段，为阐释地震发震机理及评估受灾情况提供必要的数据支撑，具有广阔的应用前景。     

利用InSAR研究西藏拉孜地震同震形变

采用欧空局ERS—1卫星雷达数据获取了1993年Mw 6.2级西藏拉孜地震同震形变场,并采用弹性位错模型反演了该地震的断层参数。结果显示该地震产生的最大地表位移约为12 cm,断层倾角为42°,不支持该断层为低角度正断层的论断。     

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

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

门源Mw5.9级地震形变InSAR观测及区域断裂带深部几何形态

利用Sentinel-1A卫星升轨、降轨合成孔径雷达影像数据，提取了2016年门源Mw5.9级地震的高精度合成孔径雷达干涉同震形变场，利用单纯形法和非负最小二乘法反演确定了地震断层几何和滑动分布，并构建了区域断裂带的深部几何形态模型。结果表明，门源Mw5.9级地震同震形变以地表抬升为主，沿升轨、降轨视线向的最大值分别为5.3 cm、7.1 cm；地震断层走向、倾角分别为133°、43°；地震滑动以逆冲为主，主要发生在地下6.14~12.28 km处，最大滑动量约0.5 m，平均滑动角为66.85°，地震矩为1.0×1018 N·m（Mw5.94）；形变观测拟合残差均方根为0.36 cm；区域断裂带的深部几何形态以花状构造为特征，整体倾向南西，门源地震发震断裂为花状构造中未出露地表的盲断层。相关成果能够为研究区域地壳运动与变形、活动断裂与地震孕育发生等提供参考。     

综合遥感解译2022年Mw 6.7青海门源地震地表破裂带

2022-01-08中国青海省门源县发生Mw 6.7地震,直接导致兰新高铁受损停运,引起了国内外的高度关注。为了评估交通网的受损情况,提出一种综合光学遥感影像、合成孔径雷达（synthetic aperture radar, SAR）影像、无人机影像和激光雷达（light detection and ranging, LiDAR）数据解译地震地表破裂带的技术框架。针对此次门源事件,首先,获取高分1号（GF-1）、高分7号（GF-7）、Sentinel-2光学遥感影像和Sentinel-1A SAR影像,根据GF-1和GF-7光学遥感影像确定地表破裂带的空间分布特征,并利用光学像素偏移量技术估计东西向和南北方向二维地表形变场;其次,利用SAR像素偏移量技术获取距离向和方位向地表形变场, 同时利用差分干涉技术获取雷达视线向的地表形变(即距离向);然后,采用运动结构恢复技术处理无人机影像获取高精度的数字地表模型;最后,综合利用上述信息精确确定地震地表破裂的空间分布和地表形变特征。结果表明,此次地震东西向最大形变量约为2.0 m,距离向最大形变量约为1.5 m,该破裂带总长约为36.22 km。结合门源地区公路交通网,基于机器学习方法支持向量机模型对历史地质灾害点的分布以及地表破裂带进行分析,发现此次地震对高速公路带来的影响最大,对乡道的影响最小;交通干线G0611和G338东南段具有很高的灾害风险。所提技术框架可精密地解译地表破裂,在地震减灾中直接发挥作用。     

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年门源地震的发生有一定的促进作用。