2021年云南漾濞MS6.4地震同震地表形变与断层滑动分布

2021年5月21日云南漾濞县城西侧发生MS6.4地震.文中以欧洲航空局升、降轨Senti-nel-1 SAR为数据源,基于D-InSAR技术获取了此次地震的InSAR同震形变场:升轨LOS(Line of Sight)向最大形变量约为0.07m,降轨最大LOS向形变量约为0.08m.以升、降轨同震形变场数据和GNSS...     

InSAR观测与小震分布联合分析2021年云南漾濞MS6.4地震应力卸载对周围区域的影响

为分析2021年5月21日云南省漾濞M_S6.4地震后震区应力变化对周围断层的影响,本文通过InSAR技术获得了漾濞地震的同震形变场,并联合小震分布数据建立断层破裂滑动模型,继而通过计算断层面上的同震库仑应力来评估此次地震对周边断层的影响,以便有效地分析地震破裂的时空解析度.结果显示:(1)在升降轨InSAR数据获得的精细同震形变场中,升轨最大视线向形变量约为5.00 cm,降轨最大视线向形变量约为7.80 cm;(2)余震精定位的主震震中为(99.89°E,25.67°N),震源深度为13.29 km,除主震之外震源深度主要集中在5—15 km;通过小震位置拟合出的发震断层走向为NW-SE(316.69°),断层倾角为88.56°,滑动角为177.97°;(3)基于InSAR同震形变场结果及小震拟合断层参数联合反演得到此次地震的断层滑动以右旋走滑为主,升轨断层最大滑动量为0.80 m,对应的深度为8.85 km,平均滑动量为0.22 m,矩震级为M_W6.41;降轨的断层最大滑动量为0.30m,对应的深度为6.88 km,平均滑动量为0.05 ...     

基于升降轨Sentinel-1 SAR数据研究2020年新疆于田MS6.4地震震源机制

2020年6月26日新疆于田县发生M_S6.4地震。 本文先利用Sentinel-1A/B卫星升、 降轨SAR数据, 结合InSAR技术提取此次地震的同震形变场; 并以同震形变场为约束, 利用贝叶斯方法的均匀滑动模型反演发震断层的几何参数; 最后基于梯度下降法(Steepest Decent Method, SDM)来确定发震断层更为精细的滑动分布。 研究结果表明: InSAR形变场的分布主要沿北北东方向, 东西方向跨度约40 km, 南北方向跨度约30 km, 形变特征为西侧沉降与东侧隆升, 西侧最大沉降约200 mm, 东侧最大隆升约70 mm。 于田地震发震断层为一条走向187.68°、 倾角59.78°、 滑动角77.76°的隐伏断层; 于田地震的宏观震源深度为6.29 km, 距震级M_W6.19。 断层滑动分布主要集中在沿走向7~21 km与倾向4~11 km的区域内, 平均滑动量约0.2 m, 在沿倾向7~8 km处的最大滑动量约0.97 m, 同震位错主要表现为正断运动。     

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

2016年新疆阿克陶地震Sentinel-1 InSAR同震形变特征

2016年11月25日新疆克孜勒苏州阿克陶县发生M_W6.6地震。 本文利用合成孔径雷达差分干涉测量技术, 对Sentinel-1卫星获取的升、 降轨雷达数据进行了处理, 提取了该次地震的同震形变场, 并结合形变场特征与震源机制解, 采用梯度下降法反演发震断层的滑动分布。 结果表明, 升、 降轨LOS向同震形变场在发震断层两侧具有明显不同的形变特征, 主要形变区域分布在断层两侧, 升轨LOS向形变量可达-8.2 cm与11.2 cm, 降轨LOS向形变量可达-21.4 cm与13.1 cm; 反演的升、 降轨干涉形变场与InSAR测量值之间的残差得到有效控制, 大部分的残差介于±5 cm之间; 断层滑动分布主要集中于沿断层面深约2~18 km处, 最大滑动量位于沿断层面深约7 km处可达0.96 m; 平均滑动角约182.29°, 最大滑动处的滑动角约197.13°, 两个滑动分布中心的滑动角均接近180°, 表明阿克陶地震为一典型的右旋走滑破裂性事件; 当剪切模量取32 Gpa时, 反演的发震断层地震矩M₀可达9.75×10¹⁸, 相当于矩震级M_W6.60, 与地震波形反演结果一致。     

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^...     

2016年新疆阿克陶MS6.7地震同震形变特征和断层滑动分布反演

本文针对2016年11月25日在新疆阿克陶发生的地震, 用差分干涉测量技术(D-InSAR)对3种不同观测模式的升、 降轨数据进行处理, 提取了多视线向的同震形变场; 根据不同模式的LOS向形变量, 构建形变分解模型, 将其分解为垂直向形变量和沿断层走向形变量; 结合同震形变场特征与震源机制解, 采用单断层模型, 利用梯度下降法(SDM), 以Multi-LOS向形变进行约束, 反演了阿克陶地震的同震滑动分布特征。 研究结果表明, 升、 降轨LOS向同震形变场在发震断层两侧具有不同的形变特征, 发震断层走向近EW向; LOS向形变量分解表明, 此次地震破裂以右旋走滑为主; 滑动分布反演的形变残差介于0~5 cm之间, 发震断层的滑动量主要位于2~16 km深部, 最大滑动量可达1.02 m, 位于断层面深部5.83 km处, 最大滑动量处的滑动角为185.24°; 平均滑动角为181.32°, 平均滑动量为0.12 m; 滑动分布反演也证明该地震为右旋走滑破裂事件, 与LOS向形变分解结果一致; 当剪切模量μ=3.2 GPa时, 反演得到的地震矩震级约M_W6.6。     

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

基于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。在川滇块体向南挤出的构造背景下,块体西边界的维西—乔后断裂、红河断裂发生右旋走滑,本次地震便是维西—乔后断裂南段分支断裂右旋走滑活动的体现。     

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,与中国地震局科考工作每日情况的结果基本一致。     

2017年九寨沟MS7.0地震同震地表三维形变场解算研究

利用覆盖九寨沟M_S7.0地震的Sentinel-1升、降轨和Radarsat-2升轨数据,分别提取了3个轨道沿雷达视线向形变,运用多平台联合观测方法解算了九寨沟地震沿地表真实的垂直向、SN向和EW向形变信息。结果表明,3个轨道InSAR数据均监测到了LOS向同震形变,范围约55km×45km,呈“果仁状”,靠近卫星飞行方向最大形变量为12.9cm (降轨),远离卫星飞行方向最大形变量为19.5cm (升轨)。三维形变结果显示,垂直向上位移达23.4cm,垂直向下位移达17.6cm;北向位移达141.8cm,南向位移达100.2cm;东向位移达22.0cm,西向位移达10.7cm。树正断裂两端地块呈非对称水平相对运动,上盘单侧向东形变较剧烈,符合左旋走滑型地震事件的运动特征。     

中小地震三维形变场重构方法研究与同震滑动分布反演——以2016年5月22日定日MW5.3地震为例

地震三维形变场对于研究地震发震机制等具有重要意义。 已有的InSAR三维形变场重构研究中, 只有中强地震的实例。 由于形变量级小、 InSAR方位向形变的误差较大, 中小地震的三维形变场重构易受噪声等影响。 本文以2016年5月22日西藏定日M_W5.3地震为例, 开展中小地震三维形变场重构的尝试。 首先基于InSAR技术获取了Sentinel-1升轨和降轨观测模式下的同震形变场, 再结合同震形变场的特点、 区域构造特征等, 添加限定方程(走滑运动为0), 重构了同震三维形变场。 结果显示, 震中附近以下降为主, 幅度达7 cm, 南北方向形变较小(约2 mm), 此区域还伴有2 cm的西向水平运动; 形变中心区域东西两侧部分区域均出现少许东向运动(1.5 cm)。 由同震形变场特征判断此次地震以正断破裂为主。 本文提出了基于连续性分层采样选取样本点方法, 以适应本地震形变场的实际情况。 对所得的LOS向位移场和重构的三维形变场进行降采样, 反演得到了断层面上的滑动分布, 两种数据得到的结果相似, 最优发震断层的走向约181°, 倾角约45°, 断层错动平均滑动角约-87.1°, 平均滑动量约为3.6 cm, 最大滑动量位于深度6.5 km处, 相当于一次M_W5.4的地震。     

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

2021年青海玛多7 4级地震的同震变形分析

2021年5月22日,青海省玛多县发生了 7.4级地震,该地震发生在巴颜喀拉地块北部边界东昆仑断裂带以南约70 km,属于块体内部断裂带地震.根据中国大陆构造环境监测网络提供资料,距离震源30多公里的玛多台站记录到东西向永久位移约25 cm.同时,InSAR也观测到明显的形变场,升轨和降轨的最大相对形变量分别约1.87 m和2.32 m.为了解释这些大地测量观测数据,本文利用该地震的三个断层滑动模型,基于不同地球模型的地震位错理论,计算同震变形场,并分别与GNSS观测数据和InSAR视线向形变量对比分析,结果显示基于InSAR数据反演的断层滑动模型产生的位移场与球形地球模型的理论计算结果最为吻合.进一步,利用较优断层模型计算2021年青海玛多7.4级地震的理论同震位移、大地水准面、重力和应变等变化,该结果为玛多地震的GNSS和重力观测的解释提供理论参考依据.     

基于InSAR的青海大柴旦地震三维同震形变场获取与震源特征分析

本文利用Envisat ASAR的升、降轨和宽幅数据,通过基于先验知识的最小二乘迭代逼近获取大柴旦2次地震的地表三维同震形变.结果表明,2008年MW6.3地震垂直向形变主要发生在断层南盘,以隆升形变为主,最大隆升量约10cm,北盘沉降量小于等于-1cm.东西向形变在南盘呈向东运动的特征,最大运动量约4cm,北盘向西运动,最大运动量约为-2cm.2009年MW5.8地震垂直向形变显示断层南盘抬升的特征,最大抬升量约27cm,北盘最大沉降量约-3cm.东西向形变表现为南盘向东运动,最大约10cm,北盘向西运动,约为-4cm.可以看出这两次地震均表现为逆冲为主,兼少量左旋走滑的震源特征.视线向结果无法判定同震形变的少量走滑特征,而地表三维分量可以有效地识别出少量左旋还是右旋走滑的震源特性.本文以视线向、垂直向、东西向形变量作为约束条件,利用Okada模型正演了2008年地震同震三维形变场.结果显示,采用逆冲兼少量左旋走滑的发震断层参数,视线向、垂直向、东西向正演结果与观测结果吻合.这也表明采用分解后的地表三维同震形变场可以有效地识别出发震断层的少量左旋走滑特征.     

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年西藏尼玛MW6.3地震的InSAR同震形变与构造意义

2020年西藏尼玛MW6.3地震发生在羌塘块体中部、依布茶卡-日干配错断裂系内的半地堑盆地内,基于震源机制解确定的发震断层存在较大差异.文中采用InSAR技术和Sentinel-1卫星升、降轨SAR数据获取了同震形变场,基于弹性半空间位错模型反演确定了发震断层参数,基于非均匀位错模型获得了断层面上的精细滑动分布.结果表明:1)在升、降轨InSAR同震形变场中,尼玛地震引起一椭圆形沉降区(长约12km,宽约8km),最大LOS向沉降值分别为-0.298m、-0.238m.2)同震位错以正断倾滑为主,兼有少许走滑分量,滑动主要集中在3～12km深度,最大滑动量达1.1 m,位于7km深处.3)发震断层为依布茶卡-日干配错断裂西侧的分支断层,走向约为30°,倾角约为68°,滑动角约为-73°.4)此次地震的破裂模式显示依布茶卡-日干配错走滑断裂存在张性应力积累,羌塘块体中部处于张性应力状态.     

利用考虑地形起伏的反演方法获取2017年伊拉克MW7.3地震的断层参数及滑动分布

2017年11月13日伊拉克北部地区苏莱曼尼亚省发生了M_W7.3地震, 造成了重大的人员和经济损失。 本文利用升降轨的Sentinel-1和降轨的ALOS-2卫星的SAR数据, 通过差分干涉测量技术获取了该地震的同震形变场, 联合DInSAR和MAI技术, 采用抗差最小二乘法求解该地震的同震三维形变场。 基于改进的考虑地形起伏的均匀位错模型反演确定了发震断层的断层参数, 最后基于非均匀位错模型得到了发震断层的分布式滑动分布模型。 结果显示: ALOS-2卫星降轨轨道观测到的伊拉克地震引起的LOS向地表形变最大为55.8 cm抬升和47.9 cm下沉; Sentinel-1卫星观测到的伊拉克地震引起的LOS向地表形变为: 升轨轨道最大为87.9 cm抬升和17.1 cm下沉; 降轨轨道最大为55.6 cm抬升和38 cm下沉; 相对于前人的研究成果, 本文利用改进的考虑地形起伏的反演方法得到的发震断层几何参数表明发震断层为NNW走向, 倾向角为352°, 同震破裂以逆冲为主, 同时兼有一定的左旋走滑分量。 基于均匀位错模型反演得到的断层滑动分布结果表明, 同震破裂未延伸至地表, 主要滑动量集中在12~18 km, 最大滑动量位于15 km深度, 达到4.3m, 反演得到的矩震级为M_W7.35, 与UGSG、 GCMT等机构给出的结果一致。     

尼泊尔M_W7.8地震InSAR同震形变场及断层滑动分布

采用DInSAR技术和欧空局2014年新发射的Sentinel-1A/IW数据,获取了2015年4月25日尼泊尔M_W7.8地震的InSAR同震形变场.所用InSAR数据扫描范围东西长约500 km,南北宽约250 km,覆盖了整个变形区域,揭示了形变场的全貌及其空间连续变化形态.此次地震造成的地表形变场总体呈现为中部宽两端窄的纺锤形,从震中向东偏南约20°方向延伸,主要形变区东西长约160 km,南北宽约110 km,由规模较大的南部隆升区和规模较小的北部沉降区组成,南部最大LOS向隆升量达1.1 m,北部最大LOS向沉降量约在0.55 m.在隆升和沉降区之间干涉纹图连续变化,没有出现由于形变梯度过大或地表破裂而导致的失相干现象,表明地震断层未破裂到地表.基于InSAR形变场和部分GPS观测数据,利用弹性半空间低倾角单一断层面模型进行了滑动分布单独反演和联合反演,三种反演结果均显示出一个明显的位于主震震中以东的滑动分布集中区,向外围衰减很快,主要滑动发生于地下7~23 km的深度范围内.InSAR单独反演的破裂范围,特别是东西向破裂长度大于GPS单独反演的破裂长度,而InSAR单独反演的最大滑动量则低于GPS单独反演的滑动量.因此认为联合反演结果更为可靠.联合反演的破裂面长约150 km,沿断层倾向宽约70 km,最大滑移量达到4.39 m,矩震级为M_W7.84,与之前用地震波数据和GPS数据反演的结果一致.     

基于D-InSAR技术的西藏改则地震同震形变场特征分析

2008年1月9日在西藏改则地区发生6.9级地震, 接着1月16日又发生6.0级余震, 表现为双震型。 本文利用差分干涉测量技术(D-INSAR), 通过对欧空局ASAR数据进行处理获取了其同震干涉形变场。 通过分析表明: 改则地震干涉形变场呈双贝壳状, 影响范围约33 km×30 km, 以北东向地震破裂带为分界线分为西北视线向沉降盘与东南视线向隆升盘, 最大视线向沉降形变量约53.2 cm, 最大视线向隆升形变量约11.3 cm。 西北沉降盘又存在东、 西两个形变中心, 推测西部形变中心受6.0级余震的控制, 东部形变中心受6.9级主震的控制。 宏观震中位置应位于左旋走滑改则—洞错断裂与依布茶卡—日干配错断裂的分阶部位(左阶), 构造应力场以张性拉伸为主, 导致地震破裂为典型的正断层破裂, 与此次地震的干涉形变场特征及哈佛大学震源机制解吻合。