利用D-InSAR技术研究西藏改则地震同震形变场

针对2008年1月9日MW6.4西藏改则地震和2008年1月16日的MW5.9余震,通过两通(2-pass)加外部DEM差分干涉处理技术(D-InSAR),提取了地震区域2次地震累积的视线向(LOS)同震形变场。结果表明:发震断层均为正断层,位于依布茶卡-日干配错断裂端点附近。主震发震断层走向为N30°E,余震发震断层走向为N21°E,两断层距离约7km;在影像上主震发震断层有造成地表破裂的痕迹,余震未见地表破裂的痕迹;这次地震造成的同震形变场长约30km,宽约20km,主震断层上盘和下盘视线向最大形变量分别为39.2cm和11.2cm,两盘相对位错达50.4cm,余震造成的视线向形变量为9.4cm     

2022年1月8日青海门源MS6.9地震同震形变场及地下位错反演

针对2022年1月8日青海门源M_S6.9地震,基于Sentinel-1A卫星的SAR影像获取了本次地震的同震形变场,并进行计算分析,运用SDM程序反演地下静态位错。初步结果如下：①本次地震震中位于蝴蝶状分布的形变图中心,抬升盘与下降盘之间的破裂迹线明显,破裂迹线长约24km,形变中心位于托莱山与冷龙岭断裂交汇处。②降轨数据显示,在断层区域视线向（LOS）形变量最小约-0.55m（远离卫星）,最大约0.68m（靠近卫星）;升轨数据显示,在断层区域视线向（LOS）形变量最小约-0.49m（远离卫星）,最大约0.42m（靠近卫星）。③对中国地震局发布的门源6.9级地震烈度图和本研究获取的地震形变图进行叠加分析,可以看出二者吻合度较高,地表有较明显形变所涉及区域约0.4万km²。④运用降轨InSAR数据反演得出地下最大位错量约3.29m,深度4.75km,在近地表位错破裂长度可达30km,地表有明显破裂的区域长约25km,与中国地震局科考工作每日情况的结果基本一致。     

玛尼地震的InSAR形变研究与模拟

采用欧洲空间局ERS-2的星载干涉雷达数据,选取1997年11月8日MW7.6级玛尼地震作为研究对象,采用了差分干涉方法,在通过对覆盖同一地区的SAR数据进行差分干涉处理,得到了玛尼地震的视线向同震形变场。经研究发现:该地震形变场呈长轴近北东东向不规则椭圆形分布,地表破裂带长度约为130km,发震断层走向约为78°,断裂为左行走滑特征。断层以南为隆起区,在发震断层附近最大视线向隆起位移量为113.6cm,断层以北为沉降区,最大视线向沉降位移量为170.4cm。基于Okada模型实现了具有复杂结构的4段断层段参数的InSAR形变场数据模拟,获得断层的最大走滑为6m,估计出玛尼地震的标量地震矩M0为2.69×10^20Nm,计算得到的矩震级MW为7.6。证明了研究方法的正确性和研究结论的可靠性。     

2017年8月8日四川九寨沟7.0级地震InSAR同震形变场及断层滑动分布反演

基于InSAR技术,利用欧空局升降轨Sentinel-1A/IW宽幅数据,获取了2017年8月8日四川九寨沟7.0级地震InSAR同震形变场,并以升降轨InSAR观测结果为约束,反演了断层滑动分布,基于三种不同接收断层计算了同震库仑应力变化.结果表明,同震形变场发生在塔藏断裂、岷江断裂和虎牙断裂交汇的三角地带,升降轨干涉位移均显示本次地震的形变场影响范围约为50 km×50 km,形变场长轴方向为NW向,升降轨观测的形变量相反,反映断层运动性质以走滑运动为主,升降轨数据观测得到的最大LOS （Line of Sight,视线向）形变量分别为~22 cm和~14 cm.非对称形变场反映出断层两侧的运动差异.反演结果显示,最大滑动量约为1 m,平均滑动角为-9°,矩震级为M_W6.5,地震破裂主要集中在地下1~15 km深度范围内,但整体而言本次地震破裂较为充分,基本将该区域1973年及1976年4次 > M_W6.0地震的破裂空区完全破裂.考虑到塔藏断裂和虎牙断裂的运动性质,可初步判定发震断层为虎牙断裂北侧延伸分支.基于三种不同接收断层模型的同震库仑应力变化计算结果反映出该区域以应力释放为主,进一步触发较大走滑型余震的可能性不大.     

基于Sentinel-1A数据的四川长宁MS6.0地震同震形变场分析及断层滑动分布反演

利用大地测量数据研究2019年6月17日四川长宁M_S6.0地震同震形变场特征和发震断层参数, 基于DInSAR技术处理升降轨Sentinel-1A数据获取干涉相位图, 并考虑大气折射效应和余震形变误差实现同震形变场改正。四叉树采样后的形变数据作为反演数据源, 采用弹性半空间分层模型反演发震断层几何面滑动分布。结果表明本次地震发震机制为兼具逆冲和左旋走滑, 矩震级为M_W5.9, 断层破裂尺度达28 km×20 km, 震源深度约9.4 km。升降轨视线向同震形变场在断层两侧呈现形变特征差异, 最大沉降量分别是8.34 cm(升轨)和4.23 cm(降轨), 最大抬升量分别是5.5 cm(升轨)和7.5 cm(降轨); 发震断层走向为302°, 倾角为43°, 平均滑动角为50°, 断层面最大滑动量达到0.28 m。     

2020年1月19日新疆伽师MS6.4地震InSAR同震形变场特征及发震构造初步探讨

基于Sentinel-1 SAR升、降影像,利用D-InSAR技术获取新疆伽师M S6.4地震的同震形变场,结果表明,本次地震引起的同震形变场整体呈近椭圆状分布,形变区东西长约66 km,南北宽约40 km,整个形变场由南部隆升区和北部沉降区组成,南部最大隆升量约7 cm,北部最大沉降量约3 cm。本次地震发生在块体俯冲界面处的低倾角逆冲推覆构造带上,隆升和沉降两个中心均位于逆冲推覆体的上盘,形变主要以隆升形变为主,符合低倾角逆断层中强震的变形特征。在沉降区与隆升区之间干涉条纹连续分布,未出现表征地表破裂位置的空间失相关带,表明地震未引起明显的地表破裂。结合震源机制、余震精定位及区域构造特征,初步推断认为伽师地震的发震构造可能为柯坪塔格推覆构造前缘的N倾的柯坪断裂。     

2001年昆仑山口西8.1级地震同震形变场特征的初步分析a

利用差分干涉雷达测量技术获取的宏观震中区的同震形变场，结合对地震活动性、震源机制、野外考察等资料分析，对昆仑山口西8.1级地震同震形变场特征进行了研究. 结果表明:宏观震中位于库赛湖东北侧，宏观震中区发震断层可分为两个形变中心区域，其中西段长约42 km，东段长约48 km，整个发震断层主破裂段长90 km；由干涉形变条纹分布格局可清楚地判断出发震断层的左旋走滑特征；断层两盘变形特征不同，南盘变形程度明显大于北盘；宏观震中附近最大斜距向位移量为288.4 cm，最小斜距向位移量为224.0 cm，宏观震中发震断层最大左旋水平位错为738.1 cm，最小地面左旋水平位错为551.8 cm.      

GPS和InSAR同震形变约束的尼泊尔M_W7.9和M_W7.3地震破裂滑动分布

2015年4月25日尼泊尔爆发MW7.9地震,继而引发5月12日MW7.3级余震,GPS、InSAR监测到震源区及周边大范围同震形变.本文以国内外的GPS和InSAR同震形变为约束,考虑喜马拉雅断裂带岩石圈垂向分层和横向差异的影响,反演主喜马拉雅逆冲断裂在这次主震和余震中破裂面形状和滑动分布.结果显示,主震从USGS确定的震中位置向东偏南延伸100km以上,破裂地面迹线与主前缘逆冲断裂迹线基本一致.破裂面倾角约7°~11°,大部分破裂集中在深度8~20km,同余震分布深度一致.主震最大滑动量约6.0~6.6m,位于14km深处.余震破裂集中在震中附近30km范围内,填补了主震东部破裂空区,最大滑动约3.6~4.6 m,位于13km深.深度20km以下基本没有破裂.地壳介质不均匀性对破裂滑动分布的影响较大,介质不均匀模型的观测值不符值比各向同性弹性半空间模型降低10%以上.本文地震破裂模型特征与地震反射剖面、以及根据震间期大地测量数据反演的喜马拉雅深部蠕滑剖面极其相似.跨喜马拉雅断裂剖面的震间形变量与地震破裂滑移量直接相关.以此推算,尼泊尔中部大震原地复发周期在300年以上.     

2022年青海门源MW6.9地震同震形变及断层滑动分布反演

2022年1月8日,青海海北州门源县发生M_W 6.9地震,震中位于青藏高原东北缘祁连地震带上冷龙岭断裂和托莱山断裂的交汇处.门源地震活动强烈,造成地表破裂明显,因此研究门源地震的发震机理,对评估周边区域及主要断裂的地震危险性具有重要意义.本文基于D-InSAR技术,利用升降轨Sentinel-1A SAR数据,获取门源地震的同震形变场.结果显示,同震形变主要集中在冷龙岭断裂和托莱山断裂的交汇处,形变长轴整体呈NWW-SEE,同震引发的升轨隆升形变量0.40 m,沉降量0.65 m,降轨隆升形变量0.80 m,沉降量0.70 m,升降轨视线向(LOS)形变呈现符号相反大小相近的特征,断裂运动以左旋走滑为主.以升降轨同震形变信息为约束条件,基于弹性半空间位错理论,采用两步法进行反演,获取门源地震断层的几何参数.结合该区域的地质构造得到断层面上的精细滑动分布,结果表明,断层破裂延伸至地表,破裂迹线长度达22 km.断层滑动主要集中在地下2～12 km,最大滑动量为4.2 m,位于地下7 km处.断层走向约为109°,倾角约为82°,释放的地震矩为2.67×10^...     

2022年青海门源MW6.6地震InSAR同震形变场与震源特征

本文利用合成孔径雷达差分干涉(Differential Interferometric Synthetic Aperture Radar, D-InSAR)技术处理哨兵卫星升降轨影像,获取了2022年青海门源M_S6.9地震的同震变形场,变形结果显示此次地震导致冷龙岭断裂西段两侧约30 km×20 km范围明显变形,卫星视线向最大同震位移达50 cm.基于有限断层模型,本文构造断层滑动分布模型,反演同震变形场数据.模型显示,冷龙岭断裂的破裂以左旋走滑为主,破裂长度约23 km,主要分布在10 km以上深度,整体形态表现为不规则椭圆盘,其中约3.9 m的最大破裂位于5～6 km深度的椭圆中心.基于反演模型估算的地震平均应力降约5.9 MPa,地震标量矩为1.0×10¹⁹ N·m,对应矩震级为M_W6.6.本文计算了1927年古浪地震、1986年及2016年门源地震对本次地震破裂面的同震静态库仑应力扰动,结果表明冷龙岭断裂西段受到历次地震的库仑应力加载作用,库仑应力增加累积达0.12 MPa,对本次发震断层的破裂有促进作用,其中...     

孕震区震前D-InSAR干涉形变场动态演化图像分析

采用D-InSAR技术获取了1997年11月8日西藏玛尼地震前3期地表干涉形变场时空演化图像。通过对干涉形变场图像进行分析,得到了震前积累形变场时空演化特征与震源破裂面闭锁之间的对应关系。结果表明:玛尼地震前10个月孕震区地表形变场就开始出现了与发震断层性质一致的变化趋势———左旋扭动,且随着发震时间的临近,南盘形变加剧,形变中心向E进一步扩展。南盘变形量明显大于北盘,南盘为主动盘。到震前2个半月,断层两盘的闭锁进一步加剧,有随时破裂的可能。从震前1996年4月16日到震前2个半月之间,南盘局部累计扭曲变形量达344mm     

基于D-InSAR技术的玉树Ms7.1地震同震形变场提取与分析

2010年4月14日青海玉树县发生了Ms7.1地震.本文利用差分干涉合成孔径雷达(D-InSAR)测量技术,通过对日本先进陆地观测卫星搭载的ALOS PALSAR获得的雷达数据进行处理,提取了此次地震的同震干涉形变场.结果表明,形变场分布特征与其发震断层甘孜-玉树断裂带的左旋走滑特征一致;最大视线向形变量为-61.4c...     

Extracting coseismic deformation of the 1997 Mani earthquake with differential interferometric SAR

Introduction The development and application of Interferometric Synthetic Aperture Radar (InSAR) have a close relationship with the sensors development of Synthetic Aperture Radar (SAR). The conception of SAR is proposed comparatively to the real aperture radar antenna. It is well known that the longer the antenna is, the higher the observation resolution will be. Just limited by the length of the antenna, the resolution of real aperture radar is generally very low and cannot meet the r…     

利用差分干涉雷达测量技术(D-InSAR)提取同震形变场h

简要介绍了合成孔径雷达干涉测量技术、差分干涉雷达测量技术，并对干涉测量精度进行了简单讨论．以西藏玛尼地区为例，通过三通差分干涉处理，获取了玛尼地震同震形变场.结果表明:形变场长200 km、宽115 km．干涉条纹以北东东向发震断层——玛尔盖茶卡断层为中心分布，且基本与发震断层平行；通过对干涉形变图进行分析，发震断层可分为3段，其中西段长约23 km，中段长约60 km，东段长约26 km，整个发震断层共长110 km；震中附近最大隆起斜距向位移量为162.4 cm，断层西侧最大沉降斜距向位移量为103.6 cm，震中最大地面水平位错为7.96 m．      

THREE-DIMENSIONAL DEFORMATION OF THE 2008 GAIZE EARTHQUAKES RESOLVED FROM INSAR MEASUREMENTS BY MULTIPLE VIEW ANGLES AND ITS TECTONIC IMPLICATIONS

The 2008 Gaize M_W6.4 earthquake,occurring on the tensional active fault zone located between Lhasa terrane and Qiangtang terrane in the interior of Tibet is a typical normal-faulting event.In this paper,we resolve the three-dimensional coseismic displacement fields of the earthquakes using a least-square iterative approximation solution with a priori knowledge,according to the theoretical basis that InSAR measurements are extremely insensitive to N-S component.Results show that the boundary dividing the two sides of the main-shock fault is very clear in the vertical movement,and two remarkable subsidence centers can be observed on the hanging wall,while amplitude of the west one (-48.9cm) is larger than the east (-41.4cm),but the maximum uplift on the footwall is only 5cm.In addition to some northward movement with amplitude less than 5cm around the aftershock fault,the north-south deformation field suggests an overall southward movement.The three-dimensional results indicate that the induced surface movement is predominantly vertical and mostly occurred on the upper side,while there are obvious east-west separation and eastward rotation in the horizontal plane.The full vectors are consistent with simulated deformation field with the RMSE less than 6cm,so the research demonstrates the feasibility of the method to recover precise three-dimensional deformation field.On the whole,the three-dimensional deformation field coincides with the tensile fracture characteristics of Gaize earthquakes,and the tectonic stress background of coeval east-west extension and north-south shortening.     

InSAR数据约束的2022年1月8日青海门源MS6.9地震发震构造研究

利用Sentinel-1A卫星升降轨道数据和D-InSAR技术获得青海门源2022年1月8日M_S6.9地震的同震形变场,并基于弹性半空间位错模型反演其震源参数,利用分布滑动模型确定断层面上的滑动分布。结果表明,2022年1月8日青海门源地震的同震形变场沿NWW-SEE方向分布;断裂带南缘升轨影像和降轨影像最大视距分别为61 cm和62 cm,断裂带北缘升轨影像和降轨影像最大视距地表形变量分别为43 cm和56 cm。InSAR同震形变场断裂尺度模型断层长30 km,宽18 km,最大滑移量3.5 m;断层滑动分布模型表明该地震为左旋走滑地震。结合冷龙岭断裂的运动特征和几何特征,初步确定此次M_S6.9地震的发震断裂为冷龙岭断裂     

2019年6月17日四川长宁MS6.0地震序列震源机制解与发震构造分析

2019年6月17日,在青藏高原东缘四川盆地南缘宜宾市长宁县发生M_S6.0地震,其后5天内相继发生了珙县M_S5.1、长宁M_S5.3和珙县M_S5.4强余震;7月4日,在珙县珙泉镇再次发生M_S5.6地震.因灾害叠加,本次地震序列导致13人死亡,200多人受伤,大量房屋受损,造成了重大的人员伤亡和财产损失.本文基于四川区域地震台网提供的地震资料,采用多阶段定位方法,对长宁M_S6.0地震序列早期（2019年6月17日至22日）余震进行了重新定位,同时,利用CAP波形反演方法,获得了序列中截止至7月4日的16次M_S≥3.6地震的震源机制解与震源矩心深度,对该序列的发震构造进行了初步分析.长宁M_S6.0地震序列重新定位后的610次M_L≥1.5地震分布显示余震区呈NW-SE向展布,长约25 km,宽5 km;序列震源深度在0~10 km区间,深度均值约3.2 km,但空间上呈西深东浅的分布特征.长宁M_S6.0地震位于余震区的东南端,具单侧破裂特征.CAP波形反演结果显示长宁M_S6.0地震序列以逆冲和逆冲兼走滑型地震为主;16次M_S≥3.6地震的震源矩心深度在1~7 km范围,平均深度3.5 km,与定位结果一致,揭示本次长宁地震序列发生在上地壳浅部.根据序列空间分布、震源机制解及震区构造特征,推测本次长宁M_S6.0地震序列的发生可能与长宁—双河复式大背斜中白象岩—狮子滩背斜和双河场褶皱及其伴生断层活动有关,位于余震区西北段的6月17日珙县M_S5.1、22日珙县M_S5.4及7月4日珙县M_S5.6地震应为6月17日长宁M_S6.0地震触发白象岩—狮子滩背斜伴生断层活动所致.序列发震构造整体呈NE-SW向挤压为主、兼具一定NW-SE向拉张分量的构造变形特征,与南侧2018年12月16日兴文M_S5.7和2019年1月3日珙县M_S5.3地震所呈现的NW-SE向挤压、NE-SW向拉张构造变形特征具有显著差异,揭示四川盆地南缘地带处于构造变形模式的转换区域,所处构造环境的变化导致本次长宁地震序列震源区及附近区域发震构造变形特征具有复杂性.     

基于InSAR分析2015年尼泊尔MW7.8地震形变场特征及断层滑动分布反演

2015年4月25日尼泊尔发生了M_W7.8地震, 本文基于震前、 震后两景Sentinel-1A雷达影像, 采用D-InSAR两轨差分干涉法提取了此次地震的同震形变场。 结果显示, 同震形变场位于喜马拉雅造山带—主边界逆冲断裂(MBT)和主前锋逆冲断裂(MFT)附近, 形变场整体表现为自西北向往东南方向延伸近150 km的纺锤形包络状, 以大面积隆起抬升形变为主, 视线向最大隆升形变达1.18 m, 抬升区北侧存在一小沉陷区, 以InSAR观测值定位同震最大形变中心。 基于均匀介质弹性半空间模型(Okada模型)与InSAR观测数据反演了断层滑动分布。 反演结果表明该地震属于典型逆冲型地震, 发震断层为主喜马拉雅逆冲断裂(MHT), 同震破裂从主喜马拉雅逆冲断裂(MHT)向上沿着主前锋逆冲断裂(MFT)传递。 基于InSAR同震形变场局部形变细节, 结合震区地质背景、 断裂分布及断层运动特征, 获得了同震破裂拟出露地表迹线。