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

基于GNSS与InSAR约束的九寨沟MS7.0地震滑动模型及其库仑应力研究

2017年8月8日的九寨沟M_S7.0地震发生在岷江断裂、塔藏断裂及虎牙断裂交汇地区,地处青藏高原东北部的川甘交界地区,位于巴颜喀拉地块的东缘,地质构造复杂,对于九寨沟地震震中位置和发震断层的确定,存在不同意见.本文利用GNSS及升降轨InSAR观测,在获取九寨沟地震同震形变场的基础上,基于均匀弹性半无限位错模型,联合反演了发震断层的滑动分布模型,并计算了同震库仑应力变化.InSAR同震形变场显示,视线向最大沉降量和抬升量分别为0.21 m和0.16 m,形变场长轴为NW向,形变主要集中在断层西侧.距震中40 km和65 km的九寨和松潘两县,水平向的GNSS同震位移分别达14.31 mm和8.22 mm.联合GNSS和InSAR同震形变场反演得到的滑动分布主要集中在沿走向5~33 km,倾向2~20 km的范围内,平均滑动量为0.18 m,最大滑动量为0.91 m.发震断层长40 km,宽30 km,走向155°,倾角81°,滑动角-9.56°.同震位移场及滑移分布模型表明此次地震为一次左旋走滑为主的地震事件,地震破裂并未完全到达地表,与虎牙断裂北段的几何产状和运动学性质更为接近,结合精定位余震的分布,我们确定虎牙断裂北段为此次地震的发震断层,震中位于北纬33.25°,东经103.82°,震源深度10.86 km,矩震量为7.754×10¹⁸ Nm,相应的矩震级为M_W6.5,与美国地调局和哈佛大学给出的震源机制解基本一致.同震库仑应力导致了虎牙断裂北段延长线的东北和西南两端应力增强,其中塔藏断裂的罗叉段和马磨段未来强震的危险性值得关注.     

2008年于田MW7.2地震滑动分布反演及三维同震形变场解算

基于ENVISAT ASAR升降轨数据, 利用InSAR获取2008年于田M_W7.2地震同震形变场; 采用SDM反演本次地震断层滑动分布; 使用PSCMP正演获取于田M_W7.2地震南北向模拟形变量, 并结合升降轨同震形变场, 解算三维同震形变场。 同震形变场分析表明, 2008年于田M_W7.2地震以正断为主, 且带有走滑运动特征, 破裂带走向为NNE向。 同震滑动分布反演结果显示, 断层沿走向被分为4段F1、 F2、 F3、 F4, 其滑动分布集中在0~14 km区间, 以F2、 F3段为主, 最大滑动量约5.31 m, 位于F2段深部2.76 km处; 沿破裂带走向, 左旋走滑位移与垂直位移比值有增大的趋势; 反演获得的地震矩M₀=5.58×10¹⁹N·m, 相当于矩震级M_W7.1。 三维同震形变场解算结果显示, 断层上盘整体表现为沉降, 断层下盘整体表现为隆升, 且沉降量明显大于隆升量, 表明地震以正断破裂为主; 除靠近断裂带中上部表现为向东南运动外, 上盘整体上表现为向西南运动; 断层下盘则整体表现为向东北运动, 证明破裂兼有左旋走滑运动。 滑动分布反演、 正演与三维同震形变场解算结果皆表明, 于田M_W7.2地震破裂以正断为主, 且带有一定的左旋走滑。     

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

利用InSAR获取2020年新疆于田Mw6.3地震同震形变场与断层滑动分布反演

本文基于InSAR技术, 利用欧空局Sentinel-1A/ B升降轨SAR数据, 提取了2020年6月26日新疆于田M_W6.3地震的同震形变场。 利用升、 降轨同震形变场约束, 分别采用MPSO算法和Bayesian方法反演此次地震发震断层均匀滑动的几何参数, 并进行对比。 然后采用SDM方法获得发震断层非均匀滑动分布, 并分析了同震库仑应力变化及其对周边断层的应力扰动。 结果表明, 同震形变场以NS向为长轴, 总体呈现西部沉降而东部隆升的特点, 而且隆升量明显小于沉降量; 滑动分布反演表明发震断层的平均滑动量为0.13 m, 平均滑动角为-104.56°, 此次地震为典型的正断破裂事件, 最大滑动量约0.55 m, 最大滑动量处滑动角为-105.2°, 位于断层倾向深度14.42 km, 释放的地震矩能量约3.65×10¹⁸ N·m, 相当于矩震级M_W6.3; 同震库仑应力变化对西部琼木孜塔格断层、 硝尔库勒南缘断层等起到应力卸载作用, 对北部的木孜鲁克—鲸鱼湖断层起到应力加载作用, 结合同震库仑应力对周围断层的扰动情况, 此次地震可能对琼木孜塔格断裂带和木孜鲁克—鲸鱼湖断裂带的西端影响较大。     

利用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。     

基于升降轨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, 同震位错主要表现为正断运动。     

2015年尼泊尔地震破裂过程的统一模型

模拟2015年尼泊尔地震（主震M_W7.8及最大余震M_W7.3） GPS/InSAR同震位移、远震体波、高频GPS位移波形和强震加速度记录,构建统一震源模型.统一模型分布特征主要由InSAR观测决定,地震矩释放过程则与P波模型相似,静态与高频GPS观测增加了对破裂时空特征的约束强度;各种比对表明,该模型对各基于单一类型反演模型具有很好的兼容性,棋盘测试展现其具有更优空间分辨率,最小可恢复20 km×20 km尺度的空间特征,压缩了非同震信号或误差导致的零散瑕疵,主、余震破裂具有更好的空间对应关系.主震展布范围为140 km×80 km;4 m以上破裂集中在加德满都以北30 km、深度15 km的狭长区域内,最大滑动量为7.4 m;破裂持续总时长为60 s,破裂速度为3.3 km·s^-1,子断层上升时间在10 s内.M_W7.3余震破裂区域位于主震东侧边缘,滑动量围绕震中扩散,扩展范围为30 km×20 km,最大滑动量约为4.4 m,总破裂持续时间为35 s.本次地震中静态和高频的GPS观测亦具备独立约束主震破裂扩展过程的能力.     

2008年10月西藏当雄MW6.3地震震后形变提取与余滑反演

2008年10月6日西藏当雄发生M_W6.3地震.本文利用震后2008年10月26日至2010年8月22日的16期ENVISAT ASAR数据,通过小基线集干涉测量、误差校正与MInTS（Multiscale InSAR Time Series）技术提取高精度的震后形变场,利用SDM（Steepest Descent Method）方法反演断层震后余滑演化过程,并分析震后余滑与同震滑动的关系.结果表明：当雄M_W6.3地震的近场震后形变场主要位于断层西侧,在时间演化上具有明显的对数函数衰减规律;震后余滑主要集中于断层中南段深0~15 km区间,最大的余滑量约0.07 m,位于断层深约9.28 km处,滑动角约-103°;震后余滑引起的地震矩能量M₀与矩震级M_W在时间演化上具有指数函数递增规律;当剪切模量μ=32 GPa,震后665天余滑释放的地震矩能量约为1.92×10¹⁷N·m,约占同震滑动释放地震矩的4.8%,相当于矩震级M_W5.46;虽然震后余滑已经延伸到断层浅部0~5 km区间,但由于余滑量相对较小,没有改变同震滑动在断层浅部区域的滑动亏损现象,这可能是2010年11月30日该区域又发生M_W5.3级余震的主要原因之一.     

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

利用考虑地形起伏的反演方法获取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_s7.0地震的InSAR观测及发震构造分析

2017年8月8日四川省九寨沟县发生M_s7.0地震.本文基于Sentinel-1 SAR影像,利用InSAR技术获取了此次地震的同震形变场,反演获得同震滑动分布,计算了同震位错对余震分布和周边断层的静态库仑应力变化,并对发震构造进行了分析讨论.结果表明:①InSAR同震形变场显示,九寨沟地震造成地表形变最大量级约为20 cm(雷达视线方向),同震形变存在非对称性分布特征.②同震位错以左旋走滑为主,主要发生在4~16 km深度,最大滑动量约为77 cm,位于9 km深处.反演得到的矩震级为Mw6.46.同震错动未破裂到地表.③大部分余震发生在库仑应力增加区.此次地震增加了震中周边地区一些断裂的库仑应力,如东昆仑断裂带东段、龙日坝断裂、虎牙断裂等.④东昆仑断裂东段的未来地震危险性值得关注.⑤九寨沟地震的发震断层为树正断裂,可能是虎牙断裂的北西延伸隐伏部分,此次地震是巴颜喀拉块体南东向运动受到华南块体的强烈阻挡过程中发生的一次典型构造事件.     

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

中小地震三维形变场重构方法研究与同震滑动分布反演——以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的地震。     

基于InSAR形变观测反演2015年皮山MS6.5地震滑动分布

2015年7月3日在新疆皮山县发生了M_S6.5地震, 该地震使当地遭受了巨大的经济损失。 本文利用欧空局提供的Sentinel-1A卫星差分干涉数据对该地震的震源机制情况进行了反演研究, 首先运用两轨法对卫星雷达影像进行差分干涉处理, 获取了覆盖皮山地震震区的同震形变场, 然后利用弹性半空间的均匀滑动模型反演获取了发震断层的几何参数, 并对原始观测数据进行降采样处理, 在此基础上运用分布式滑动模型反演获取了更为精细的断层滑动分布, 结果显示分布式滑动模型与观测结果有很高的拟合度。 反演结果表明发震断层是一个以逆冲为主兼有极少量左旋走滑的盲断层, 此次地震断层面的同震活动分布主要集中在7~15 km深度范围内, 同震的地震矩为6.28×10¹⁸N·m, 矩震级为M_W6.46, 与前人的研究结果非常一致。