利用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,分析得出该地震为向北逆冲兼左旋走滑。     

利用InSAR资料反演缅甸Mw6.8地震断层滑动分布

2011年3月24日缅甸东北部发生M_w6.8级地震.本文利用覆盖该地区的升降轨ALOS PALSAR数据,获取了该次地震的同震形变场,并采用灰度配准技术获取了其地表破裂位移.针对影像中因轨道不精确造成的非线性长波长误差,本文采用二次多项式曲面法予以去除,获取了更为精确的同震形变场.最后,基于弹性半空间形变模型反演了该地震的断层滑动分布.结果表明,该地震断层滑动以左旋走滑为主,兼具少许的倾滑运动分量,断层滑动主要集中分布在断层面0~10 km深度范围,最大滑动量达5 m,位于地表以下5 km深处.反演获得的地震标量矩为2.49×10¹⁹N·m,震级约为M_w6.8级.     

Source parameters of the 2009 L'Aquila earthquake,Italy from Envisat and ALOS satellite SAR images

2009年4月6日意大利L’Aquila地区发生了Mw6.3级地震,该地震造成了300余人的人员死亡.本文联合不同波长、不同入射倾角的升降轨Envisat和ALOS卫星的差分干涉数据对该地震进行震源机制解的反演研究.研究首先对卫星雷达影像进行二通差分干涉处理,获取了覆盖L' Aquila地震震区的完整InSAR同震形变场,然后结合四叉树和均匀采样方法对原始观测数据进行降采样.在此基础上,联合GPS形变观测数据,利用弹性半空间矩形和三角位错模型,以及断层自动剖分技术对断层面进行最优离散剖分,反演获取了发震断层的精确几何参数和最优断层滑动分布,结果显示分布式三角位错滑动模型能够很好地解释观测到的地表形变场.反演结果表明发震断层是一个以正倾滑为主兼有少量右旋走滑的盲断层;基于观测数据最优确定的断层单元的最短边长为0.4 km,最长边长为6.3 km;此次地震的滑动分布主要发生在5～14 km深度的范围内,最大滑移量为1.07 m,释放的能量为3.43×1018N·m(Mw6.32),与地震学的研究结果非常一致.     

联合Envisat和ALOS卫星影像确定 L'Aquila地震震源机制

2009年4月6日意大利L'Aquila地区发生了M_w6.3级地震,该地震造成了300余人的人员死亡. 本文联合不同波长、不同入射倾角的升降轨Envisat和ALOS卫星的差分干涉数据对该地震进行震源机制解的反演研究. 研究首先对卫星雷达影像进行二通差分干涉处理,获取了覆盖L'Aquila地震震区的完整InSAR同震形变场,然后结合四叉树和均匀采样方法对原始观测数据进行降采样. 在此基础上,联合GPS形变观测数据,利用弹性半空间矩形和三角位错模型,以及断层自动剖分技术对断层面进行最优离散剖分,反演获取了发震断层的精确几何参数和最优断层滑动分布,结果显示分布式三角位错滑动模型能够很好地解释观测到的地表形变场. 反演结果表明发震断层是一个以正倾滑为主兼有少量右旋走滑的盲断层;基于观测数据最优确定的断层单元的最短边长为0.4 km,最长边长为6.3 km;此次地震的滑动分布主要发生在5～14 km深度的范围内,最大滑移量为1.07 m,释放的能量为3.43×10¹⁸ N·m(M_w6.32),与地震学的研究结果非常一致.     

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向的维西—乔后断裂西侧的隐伏次生断裂。     

基于线弹性位错模型反演1997年西藏玛尼Mw7.5级地震的干涉雷达同震形变场——Ⅱ滑动分布反演

1997年11月8日西藏M_w7.5级玛尼地震是干涉雷达技术应用于地震观测以来的一次重要事件.在第一部分中，我们应用广泛使用的Okada线弹性位错模型，假设断层的各个分段滑动量均匀，反演得到断层各个分段的几何参数和均匀滑动量.本部分的反演进一步去除滑动均匀假设，并利用更能反映断层真实状态的角形元位错模型（线弹性），在第一部分反演得到断层几何的基础上，反演断层面的静态位错分布.反演结果表明，线弹性滑动分布模型能够更好地解释观测数据，进一步提高反演的数据拟合程度.最终得到了断层面上的走滑和倾滑位错分布.首次得到的断层面滑动分布显示断层面滑动在浅部（0～12 km）比较集中，地震破裂长度约170 km，最大左旋走滑位移达4.8 m；反演结果还表明局部段落存在较大倾滑位移，量值达到1.9 m，这在断层模型中是不能忽略的，它可能是断层两侧形变不对称的原因之一；反演得到的标量地震矩为2.18×10²⁰ N·m，相当于矩震级M_w7.5，与Velasco等利用地震波形反演得到的结果一致.     

利用考虑地形起伏的反演方法获取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等机构给出的结果一致。     

联合GNSS和InSAR观测位移反演2008年汶川大地震断层位错模型参数

利用现代空间大地测量技术,尤其是卫星合成孔径雷达干涉测量,能够获取高精度、高空间分辨率的同震和孕震形变,为地震断层形变和破裂机制研究提供了前所未有的机遇。本文介绍了利用大地测量观测数据反演地震断层位错模型参数的贝叶斯反演方法。联合运用2008汶川大地震前后GNSS和InSAR技术观测获得的同震位移,反演了地震断层的几何参数和滑动位错分布。研究结果表明,汶川地震的断层滑动主要集中在倾角较陡的浅部,同时包含逆冲和右旋走滑,其中最大逆冲6.1m,最大右旋6.5m。根据断层滑动分布正演计算得到的上盘同震位移明显小于下盘,预示该断层两侧孕震形变可能存在较大的不对称性。     

断层走滑不均匀性对地面变形的影响

根据塞积群理论导出了走滑断层两端点固定时位移场随位置的不均匀分布解析表示关系,分析和讨论了断层上不同位置的位移场分布特征, 采用新导出的断层滑动不均匀分布公式对地面的水平位移场和垂直位移场进行了数值模拟计算,并与传统的Okada 位错模型在理论和计算结果上进行了对比分析.理论和计算结果分析表明：断层滑动不均匀分布公式是Okada位错模型在一定条件下理论上的扩展,而Okada位错模型是走滑不均匀公式的零阶近似;断层滑动不均匀分布公式与Okada 位错模型计算的地面位移场在走滑方向、垂直于断层走向和垂直于地面方向的近场变形差异最大在50%至65%之间,而最小差异在1%以内.     

芦山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处展布一条隐伏断裂,为本次的发震断裂.     

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地震破裂以正断为主, 且带有一定的左旋走滑。     

基于线弹性位错模型反演1997年西藏玛尼Mw7.5级地震的干涉雷达同震形变场——Ⅱ滑动分布反演

1997年11月8日西藏M_w7.5级玛尼地震是干涉雷达技术应用于地震观测以来的一次重要事件.在第一部分中，我们应用广泛使用的Okada线弹性位错模型，假设断层的各个分段滑动量均匀，反演得到断层各个分段的几何参数和均匀滑动量.本部分的反演进一步去除滑动均匀假设，并利用更能反映断层真实状态的角形元位错模型（线弹性），在第一部分反演得到断层几何的基础上，反演断层面的静态位错分布.反演结果表明，线弹性滑动分布模型能够更好地解释观测数据，进一步提高反演的数据拟合程度.最终得到了断层面上的走滑和倾滑位错分布.首次得到的断层面滑动分布显示断层面滑动在浅部（0～12 km）比较集中，地震破裂长度约170 km，最大左旋走滑位移达4.8 m；反演结果还表明局部段落存在较大倾滑位移，量值达到1.9 m，这在断层模型中是不能忽略的，它可能是断层两侧形变不对称的原因之一；反演得到的标量地震矩为2.18×10²⁰ N·m，相当于矩震级M_w7.5，与Velasco等利用地震波形反演得到的结果一致.     

汶川大地震震后重力变化和形变的黏弹分层模拟

基于有限矩形位错理论及陈运泰等、JiChen等通过地震波反演的断层模型,结合研究区地壳——上地幔平均波速分层结构,利用PSGRN/PSCMP软件模拟计算了黏弹分层半空间中汶川地震（Ms8.0）产生的同震及其震后地表形变和重力变化,同时给出了震后形变和重力变化的年变化率．模拟结果表明,同震形变和重力变化显示出发震断层倾滑逆冲兼具右旋走滑综合特征;其变化主要发生于断层在地表的投影区附近,震后地表重力变化和形变量均不断增大,影响的范围也不断扩张;震后50a间近场年均形变量可达10mm,年均重力变化量可达2times;10-8m/s2,而远场年均形变量一般低于2mm,年均重力变化量一般低于0.4times;10-8m/s2;形变和重力变化在震后200a内变化较为显著,变化率逐渐减小,水平位移在400a后基本稳定不变,垂直位移、重力变化和大地水准面变化在800a后基本稳定不变．      

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)此次地震的破裂模式显示依布茶卡-日干配错走滑断裂存在张性应力积累,羌塘块体中部处于张性应力状态.     

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

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

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数据的2017年伊拉克哈莱卜杰MW7.3地震同震形变场分析及断层滑动分布反演

本文首先利用二轨法对欧洲空间局Sentinel-1A雷达卫星影像进行差分干涉处理,获取了覆盖2017年伊拉克哈莱卜杰（Halabjah）M_W7.3地震震区的同震形变场,结果表明:哈莱卜杰地震造成的地表形变影响范围约为60 km×70 km,形变场基本沿扎格罗斯主前缘断层展布;形变场的西南盘呈现隆升趋势,最大视线向形变值为88 cm,东北盘呈现下降趋势,最大视线向形变值为37 cm;隆升形变值远大于沉降值,反映出发震断层以逆冲运动为主的特征。然后基于弹性半平面空间矩形位错模型,分别采用多峰值粒子群算法和最速下降法确定了发震断层的几何参数和滑动分布结果。反演结果显示发震断层以逆冲运动为主,兼少量右旋走滑运动,最大滑动量为3.34 m,释放的地震矩为1.68×1020 N·m （M_W7.4）,与地震学的研究结果一致。      

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

基于线弹性位错模型及干涉雷达同震形变场反演1997年玛尼Mw75级地震参数－Ⅰ.均匀滑动反演

1997年11月8日西藏玛尼Mw75级地震是干涉雷达技术应用于地震观测以来的一次重要事件.Peltzer等认为该形变场体现了地壳的非线性弹性形变响应，挤压象限和拉张象限的杨氏模量比值为2.相同的情形并没有在其他地震同震形变场中发现，通常用线弹性理论就能够很好地解释观测数据.考虑到他们所用断层模型的简化程度和纯走滑约束等，本研究认为非线性弹性解释是牵强的.应用广泛使用的Okada线弹性位错模型，采用卫星观测得到的地表断层位置，去除倾滑为零的约束，基于非线性优化反演方法寻求拟合雷达观测的最佳断层几何参数和均匀滑动参数.结果表明，线弹性模型能够满意解释观测数据.断层在朝阳湖以东的段落最深达到165 km，随着断层向两侧延伸，深度逐渐减小；反演得到的断层倾     

基于线弹性位错模型及干涉雷达同震形变场反演1997年玛尼Mw75级地震参数－Ⅰ.均匀滑动反演

1997年11月8日西藏玛尼Mw75级地震是干涉雷达技术应用于地震观测以来的一次重要事件.Peltzer等认为该形变场体现了地壳的非线性弹性形变响应，挤压象限和拉张象限的杨氏模量比值为2.相同的情形并没有在其他地震同震形变场中发现，通常用线弹性理论就能够很好地解释观测数据.考虑到他们所用断层模型的简化程度和纯走滑约束等，本研究认为非线性弹性解释是牵强的.应用广泛使用的Okada线弹性位错模型，采用卫星观测得到的地表断层位置，去除倾滑为零的约束，基于非线性优化反演方法寻求拟合雷达观测的最佳断层几何参数和均匀滑动参数.结果表明，线弹性模型能够满意解释观测数据.断层在朝阳湖以东的段落最深达到165 km，随着断层向两侧延伸，深度逐渐减小；反演得到的断层倾