利用InSAR技术研究2016年青海门源MW5.9地震同震形变场及断层滑动分布

2016年1月21日, 青海省门源县冷龙岭断裂带附近发生了M_W5.9地震。 基于Sentinel-1A影像, 采用差分干涉雷达测量技术研究了此次地震产生的同震形变场, 结果表明, 门源地震的形变影响范围约20~30 km, 形变态势在升降轨道形变场均显示为隆升, 基本沿冷龙岭断裂呈近似同心圆展布, 推测可能是冷龙岭断裂与民乐—大马营断裂之间的一条逆断层, 沿雷达视线方向最大形变量级约为6 cm。 均匀滑动反演显示门源发震断层长7.3 km, 宽6.2 km, 走向298.6°, 倾角34.5°, 倾向宽度9.5 km, 沿走向滑动量为170 mm, 沿倾向滑动量为460 mm, 矩震级为M_W5.97; 分布式滑动反演显示门源地震以逆冲为主, 兼具少量右旋走滑分量, 滑动量主要集中在沿断层倾向方向, 距离地表5～15 km处, 最大滑动量约0.3 m, 位于断层倾向深度10 km处, 矩震级为M_W5.93。     

2020年伽师MS6.4地震InSAR同震形变与滑动分布特征分析

利用Sentinel-1A升轨和降轨数据,基于D-InSAR技术,获取2020年1月19日伽师M_S6.4地震同震形变场,并结合其他研究机构给出的震源机制解参数和已有研究成果,反演得到伽师地震的发震断层几何特征和滑动分布。研究结果表明,伽师地震同震形变在地表有明显差异;升轨同震形变在卫星视线方向北侧抬升55 mm,南侧下降42 mm;降轨同震形变在卫星视线方面北侧抬升63 mm,南侧下降23 mm。通过反演得到发震断层走向为275°,倾角为20°,地震滑动主要分布在地下5 km处,最大滑动量约为0.32 m,平均滑动角为89.3°,累积地震矩为1.46×1018 N·m,合矩震级M_W6.1,发震构造为具有少量走滑性质的逆冲断裂。从发震构造特征、同震滑动分布推测,伽师地震发震构造是柯坪塔格褶皱带滑脱面以上沉积盖层内的逆冲断裂,支持了柯坪推覆体的薄皮构造模型观点。     

InSAR数据约束的2021年5月21日云南漾濞MS6.4地震发震构造研究

运用Sentinel-1A卫星数据和D-InSAR技术,获取2021-05-21云南漾濞M_S6.4地震的同震形变场。结果显示,漾濞地震同震形变场长轴近NW展布升降轨形变场符号相反,视线向最大沉降量和抬升量为0.1 m。InSAR同震形变场反演的滑动分布主要集中在沿走向2～12 km,倾向1～9 km的范围内,最大滑动量0.35 m,发震断层长9.8 km、宽4 km,滑动量主要集中在地下3～6 km范围内,滑动角-146.7°。同震位移场及滑动分布模型反映本次地震为发震断层的右旋走滑事件,地震破裂未达到地表。断层模型反演结果显示,矩震级为M_W6.1,发震断层以北西走向右旋走滑运动为主,初步认为本次M_W6.1地震发震断裂可能是一条NW向的维西—乔后断裂西侧的隐伏次生断裂。     

2020年1月19日新疆伽师县MS6.4级地震InSAR同震形变场与断层滑动分布反演

利用欧空局Sentinel-1A SAR数据,重建了2020年1月19日新疆伽师县MS6.4地震InSAR同震形变场.以升降轨InSAR形变场联合约束,反演了发震断层参数与同震滑动分布.结果 表明,地震同震形变场(长轴走向近EW向)发生在柯坪塔格褶皱带与奥兹格尔它乌褶皱带之间的区域内,升降轨观测的视线向(LOS)形变量相同,反映出发震断层运动性质是以逆冲为主,最大LOS向形变量约7 cm,干涉形变场呈现非对称性分布.断层模型反演结果显示,破裂区域长度为24 km,宽8 km,集中于地下深度4～6km范围内,同震最大滑动量约0.34 m,平均滑动角85°,矩震级为MW6.0,北倾(倾角15°)逆冲推覆运动主导着断层破裂,兼具少量走滑运动分量.从发震断层模型、同震滑动深度分布及破裂运动学特征推测,这次伽师地震的发震构造是柯坪塔格褶皱带山前出露的柯坪塔格逆断裂,且支持柯坪塔格地区的薄皮构造模型.     

基于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。     

基于InSAR同震形变场反演汶川M_w7.9地震断层滑动分布

通过综合分析2008年5月12日汶川地震野外地震地质考察的地表破裂带空间分布及分段资料,结合InSAR干涉形变场资料,构建了五段断层儿何结构模型,该模型与野外地震地质考察结果在多数分段上基本一致;基于此五段断层模型,运用敏感性迭代拟合算法反演了汶川地震InSAR同震形变场,获得了断层滑动分布及部分震源参数.结果表明,基于余震精定他获得的地震断层倾角模型模拟的同震形变场与InSAR形变场吻合较好,且残差较小;反演的滑动分布主要集中于地下0～20 km,最大滑动量分别位于北川及青川等地区,最大可达到10 m;沿SW-NE走向,断层面的滑动方向主要以右旋兼逆冲形式为主,在汶川及都江堰地区以强烈的逆冲为主兼有一定右旋走滑分量,在北川及映秀地区以逆冲兼右旋运动为主,在平武及青川等地区则逐渐过渡为以右旋运动为主兼有一定的逆冲分量,其中汶川地区的平均滑动角为97°,北川地区的平均滑动角为119°,青川地区平均滑动角为138°.反演矩张量为7.7×10~(20)N·m,矩震级达M_w7.9.     

基于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同震形变场局部形变细节, 结合震区地质背景、 断裂分布及断层运动特征, 获得了同震破裂拟出露地表迹线。     

基于InSAR同震形变场反演汶川Mw7.9地震断层滑动分布

通过综合分析2008年5月12日汶川地震野外地震地质考察的地表破裂带空间分布及分段资料,结合InSAR干涉形变场资料,构建了五段断层几何结构模型,该模型与野外地震地质考察结果在多数分段上基本一致;基于此五段断层模型,运用敏感性迭代拟合算法反演了汶川地震InSAR同震形变场,获得了断层滑动分布及部分震源参数.结果表明,基于余震精定位获得的地震断层倾角模型模拟的同震形变场与InSAR形变场吻合较好,且残差较小;反演的滑动分布主要集中于地下0～20 km,最大滑动量分别位于北川及青川等地区,最大可达到10 m;沿SW-NE走向,断层面的滑动方向主要以右旋兼逆冲形式为主,在汶川及都江堰地区以强烈的逆冲为主兼有一定右旋走滑分量,在北川及映秀地区以逆冲兼右旋运动为主,在平武及青川等地区则逐渐过渡为以右旋运动为主兼有一定的逆冲分量,其中汶川地区的平均滑动角为97°,北川地区的平均滑动角为119°,青川地区平均滑动角为138°.反演矩张量为7.7×10²⁰ N·m,矩震级达M_w7.9.     

根据GPS和水准资料反演2008年汶川地震同震破裂模型

通过分析2008年汶川8.0级地震前后的GPS观测资料和1983～2010年的精密水准观测资料,得到地表同震位移场。结合野外地质调查,反演了汶川地震同震断层几何模型和断层面上滑动分布。反演结果表明:北川—映秀断裂是一个铲形断层,长度为180km,地表倾角为62°,随指数函数形式逐渐变缓,深度18km,断层面滑动以逆冲为主;青川断裂长为120km,倾角为58°,断层面上走滑分量由南向北逐渐增加;灌县—江油断裂为纯逆冲断层,长度为80km,倾角35°。整个断层模型中最大错动量达到7.6m,对应破坏最严重的北川地区。本次地震释放的地震矩为7.62×1020N·m,相应矩震级为MW7.9。反演所得模型对于近场的水准观测资料和GPS观测资料拟合很好。     

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地震的发震断裂为冷龙岭断裂     

基于多期DInSAR数据的2016年定结地震序列发震断层特征研究

2016年5月22日,在西藏定结县发生四次MW4~5地震,研究本序列地震的发震断层几何和运动特征对于认识周边活动断裂性质具有重要意义.由于发震地区偏远,且观测台网分布稀疏,本研究采用星载雷达干涉测量(DInSAR)技术进行了同震形变场重建,但是定结地震震级较小,单干涉像对获取的形变场受相位噪声影响较大.为了解决这一问题,本研究基于时间序列Sentinel-1A干涉数据生成多期同震与非同震干涉图,并利用叠加平均法对本次定结地震同震形变场进行重建,提取了定结2016年5月22日多次地震产生的同震累计整体形变场.基于InSAR同震形变场和区域地质特征,研究进行了滑动分布反演,确定其主要贡献的发震断层几何参数及滑动分布:断层走向为188°,倾角为43°,平均滑动角为78°,发震断层的运动性质以正断为主兼具少量左旋走滑分量,滑动主要集中在断层垂直深度0~9km处,最大滑动量约为25cm,位于断层倾向深度3km处,反演得到的矩震级为MW5.58.本研究结果表明采用星载InSAR叠加平均技术可以较好地压制相位噪声,有效提取此类中小型浅源地震同震微弱形变场.最后,我们认为本次定结地震与藏南拆离断层与申扎-定结断层的活动密切相关.     

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地震的破裂空区完全破裂.考虑到塔藏断裂和虎牙断裂的运动性质,可初步判定发震断层为虎牙断裂北侧延伸分支.基于三种不同接收断层模型的同震库仑应力变化计算结果反映出该区域以应力释放为主,进一步触发较大走滑型余震的可能性不大.     

2015年11月23日祁连5.2级地震发震构造初步研究

基于喜马拉雅地震科学台阵和青海区域数字地震台网的资料,采用广义极性振幅技术反演了2015年11月23日祁连5.2级地震的震源机制;利用逆时成像技术重新定位了祁连5.2级地震及64次ML1.0以上余震的震源位置。基于上述研究,综合分析祁连5.2级地震的震源位置和震源机制以及64次ML1.0以上的余震空间分布特征,结合托莱山断裂构造性质探讨了发震构造。结果显示,祁连5.2级地震的发震时刻为北京时间2015年11月23日5时2分38.9秒,震中位置位于(37.95°N,100.46°E),震源深度为12.4km。祁连5.2级地震的震源机制为节面Ⅰ的走向108°/倾角44°/滑动角40°,节面Ⅱ的走向347°/倾角63°/滑动角126°。节面Ⅰ与托莱山断裂左旋走滑兼具逆冲的性质相同,也与余震勾勒出的断层面倾向SW,倾角约48°的产状相同,因此节面Ⅰ为发震断层面。结合震源机制结果和托莱山断裂的构造性质,推测主震的发震构造为一条北西西向的断层,倾向SW,倾角在深部较缓而在浅部可能较陡。由于托莱山断裂带次级断裂发育、产状复杂、缺乏准确位置,因此无法通过定位结果来判断发震构造为托莱山主断裂还是其次级断裂。     

利用InSAR技术研究2018年西藏谢通门MW5.6地震发震构造

2018年12月24日在西藏日喀则市谢通门县发生MW5.6地震。利用Sentinel-1A卫星数据和D-InSAR技术,获得该地震LOS方向的同震形变场图像,并以此为约束反演地震断层的几何参数和同震滑动分布特征。反演结果表明,谢通门MW5.6地震发震断层近南北走向,断面倾向西,倾角约35°,破裂长度约12 km,同震滑动主要集中在4~11 km深度范围内,以正断倾滑为主,最大滑动量0.18 m,矩震级MW5.6。2018年谢通门地震发震构造是谢通门—申扎伸展断裂系中的一条新生盲断层。     

芦山MS7.0级地震InSAR形变观测及震源参数反演

2013年4月20四川省芦山县发生M_S7.0级地震,目前的研究资料表明地震发生在龙门山断裂南段,但地表未发现明显破裂.本研究利用InSAR技术与Radarsat-2雷达数据,获取了芦山地震同震的部分形变场,结果表明,近场区域的LOS位移发生视线向隆升,量值在7 cm左右.随后利用弹性半空间的位错模型反演了断层面参数,综合反演结果及震源机制解最终确定了发震断层的初始模型,以形变场观测数据为约束,基于梯度下降法反演获得了断层面上的滑动分布,反演得到的矩震级为M_w6.45级,断层走向213°,倾角39°~43°,最大滑动位于地表以下约13 km深度位置,最大滑动量0.91 m,平均滑动角71°,整体上仍以逆冲滑动为主,兼具左旋走滑.推测在双石-大川断裂以东12 km处展布一条隐伏断裂,为本次的发震断裂.     

2017年8月9日精河MS6.6地震余震序列精定位及发震构造分析

为确定2017年8月9日精河M_S6.6地震的发震构造,本文使用双差定位方法对发震时刻至2017年10月震源区所发生的余震进行了精定位,同时利用CAP波形反演方法,得到了主震的震源机制解,同时使用GPAT方法反演得到了部分余震的震源机制解,并基于两者对本次地震的发震构造予以分析。结果显示:精定位后主震位于（44.27°N,82.85°E）,震源深度为17 km;主震最佳双力偶解对应的节面Ⅰ的走向为260°、倾角为51°、滑动角为84°,节面Ⅱ的走向为89.5°、倾角为39.4°、滑动角为97.4°;余震序列位于主震东侧,并向东展布约30 km,在3—18 km深度范围内均有分布,其优势方向为近EW向,次优势方向为SW向。结果表明,本次地震是一次逆冲型地震,通过反演得到的大量小震震源机制解的结果与主震震源机制解结果相一致。结合余震震中分布、主震及余震的震源机制解以及震源区的地质构造,本文推断近EW走向具有逆冲性质的库松木楔克山前断裂为精河主震的发震构造。      

2012年11月20日宁夏永宁MS4.6 地震震源机制解与发震构造

2012年11月20日在宁夏银川市永宁县与兴庆区交界处发生M_S4.6地震,为了更好地了解此次地震的发震构造,首先采用Hypo2000绝对定位方法得到该地震的震中位置及余震分布;然后采用CAP方法反演了此次地震的震源机制解和震源深度. 反演结果表明,永宁M_S4.6地震是一个带有少量逆冲分量的右旋走滑地震.该地震矩震级为M_W4.3,最佳双力偶解为:节面Ⅰ走向11°,倾角74°,滑动角171°;节面Ⅱ走向103°,倾角81°,滑动角16°.最佳震源深度为8km左右.从该地震震中和震源机制解以及震源深度剖面分布来看,这次地震很可能发生在银川隐伏主断层西侧的次级断层上.      

基于Sentinel-1A数据反演漾濞MS6.4地震的同震形变场及断层几何参数

2021年5月21日晚21时48分,云南省大理州漾濞县(震中:25.67°N,99.87°E)发生M_S6.4地震,震源深度8 km。为快速获得此次地震同震形变场及断层几何参数,研究该次地震的发震构造等,文章基于震前、震后的sentinel-1A卫星升降轨SAR数据进行二轨法差分雷达干涉测量(DInSAR),并基于Okada弹性半空间位错模型反演断层几何参数。研究结果如下:(1)此次地震造成的同震形变场长约19 km,宽约20 km;(2)升轨雷达视线向最大形变约为8.2 cm,降轨雷达视线向最大形变约为8.7 cm;(3)地震断层走向为313.7°,倾角为87°,滑动角为175°,为右旋走滑型断层,最大滑动量为0.79 m,反演得出的地震矩为1.48×10~(18) N·m,矩震级为M_W6.1。在川滇块体向南挤出的构造背景下,块体西边界的维西—乔后断裂、红河断裂发生右旋走滑,本次地震便是维西—乔后断裂南段分支断裂右旋走滑活动的体现。     

影像大地测量技术再探伊朗桑塞菲德2017年Mw6.1地震盲断层倾向

在盲断层参数未知的情况下,确定断层模型参数是获取震源机制的基础,也是地震风险性分析的首要任务.2017年伊朗桑塞菲德Mw6.1地震倾角的确定作为盲断层参数确定的一个宝贵实例,在研究观点上仍存在诸多争议.为更好地解释此次地震断层倾角,本文先基于弹性半空间模型构建了东北倾和南倾两种备选的单一平面断层模型.在此基础上,利用哨兵1号和ALOS-2提供的影像大地测量数据联合解算的二维同震形变揭示了两种不同断层模型的逆冲走滑运动特性以及俯冲带逆冲断层破裂模式.为进一步分析更多断层细节并充分解释断层倾向,本文以二维形变数据为约束采用最小二乘动力学反演和贝叶斯反演分别获取了两种潜在断层的最优断层参数和滑动分布并着重分析了两种方法反演结果的各向异性,结果表明东北倾断层具有较小残差且滑动分布更合理.通过分析主震后重定位的59次余震结果发现余震走向分布以及深度延伸方向与东北倾断层模型更吻合,并且东北倾向断层走向与伊朗东北部已知活动断层基本平行,因此可以作为桑塞菲德地震的发震断层为WNW ESE走向的进一步佐证.本文研究的断层参数为走向314.67°,倾角41.5°,平均滑动角128.8°,断层特征为东北倾逆冲兼具右旋走滑,滑动分布主要集中在沿走向6～22 km处,最大滑动量为7 km深处的1.09 m,断层破裂并未到达地表,反演矩张量为1.427×1018N·m,相当于Mw6.07矩震级地震.     

利用InSAR研究乌恰地震同震形变

2008年10月5日在新疆乌恰地区39.50°N、73.64°E处发生MW6.7地震,用ROI_PAC软件处理4景日本ALOS卫星数据得到乌恰地震的三维同震形变场。距离向,断层北侧最大位移达39 cm,南侧最大位移达36 cm;方位向,北侧和南侧的最大位移分别为1.5 m和2 m。采用Okada弹性位错模型分析该地震的滑动分布,模拟得出的断层走滑角为48°,倾滑角为53°,深度为10.5 km,最大形变量2.5 m,分析得出该地震为向北逆冲兼左旋走滑。