基于线弹性位错模型反演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等利用地震波形反演得到的结果一致.     

基于线弹性位错模型及干涉雷达同震形变场反演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，随着断层向两侧延伸，深度逐渐减小；反演得到的断层倾     

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

2022年泸定6.8级地震GNSS同震形变场及其约束反演的破裂滑动分布

2022年9月5日四川省甘孜州泸定县发生6.8级地震,地震发生在川滇菱形块体东边界鲜水河断裂磨西段西侧附近.本文基于震中350 km范围内中国大陆构造环境网络和中国地震科学实验场118个GNSS连续站数据,观测获得了精细的同震形变场,结果显示：同震形变跨发震断层呈空间对称分布,表明本次地震具有显著的左旋走滑特征;记录到的最大同震形变发生在震中距40 km的SYD5(石棉安顺场)站,东西向和南北向形变量分别达到-22.0±1.2 mm和11.6±0.9 mm,震中距100 km以外测站的水平同震形变均小于5 mm;垂向同震形变不显著.结合走向163°和倾角77°的发震断层模型,本文对发震断层面的同震破裂滑动分布进行了反演,结果显示：同震滑动主要集中在主震东南侧余震空区0~15 km深度范围内,且破裂达到了地表,同震释放的地震矩与矩震级为M_W6.57地震相当.结合理论同震形变场、主应变场和邻近区域主要活动断裂库仑应力变化的空间分布特征,本次可能导致处于强闭锁和地震空区的安宁河断裂石棉—冕宁段未来发震风险性增强.

     

2020年西藏尼玛MS6.6地震InSAR同震形变场及断层滑动分布反演

2020年7月23日,西藏尼玛县发生MS 6.6地震,研究本地区的发震断层对认识周边活动断裂有着重要意义.由于发震地区位于青藏高原腹地,地震监测台站分布稀疏,监测能力较弱.故本文基于升、降轨Sentinel-1A SAR数据,通过差分干涉技术获取了西藏尼玛MS 6.6地震InSAR同震形变场,显示本次地震所引起的地表形...     

2010年1月12日海地 Mw7.0地震InSAR同震形变观测及同震滑动分布反演

2010年1月12日GMT时间21时53分,在海地境内(72.57°W,18.44°N)发生了Mw7.0地震.文中利用干涉合成孔径雷达(InSAR)方法获得了覆盖整个震区的高精度形变观测资料,用以研究该地震的发震机理.采用ALOs PALSAR数据,分析了轨道、大气等误差源对干涉信号的影响,最终获得了雷达视线向(LOS...     

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

     

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

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

2008年MS7.1于田地震InSAR同震形变场及其震源滑动反演

本研究利用InSAR技术与ALOS PALSAR雷达数据,获取了2008年3月20日于田MS7.1地震视线向同震形变场,并基于该数据集和限制性最小二乘算法反演了此次地震的断层滑动分布;通过构造四大类反演方案,详细分析了InSAR观测系统中的入射角与方位角对反演结果的影响.结果表明:入射角随点位变化对反演结果有较大影响,使用其平均值将对破裂细节产生一定影响;而方位角对反演结果的影响不大,使用其平均值是一种较为理想的选择;引入入射角与方位角变化后,反演获得了较佳的于田地震同震滑动,主要集中分布于0~14 km深度附近,最大滑动量达3.2 m,矩张量为3.3×1019 N·m,相当于矩震级MW7.0.     

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

尼泊尔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数据反演的结果一致.     

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形变观测反演2015年皮山MS6.5地震滑动分布

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

雷达差分干涉测量技术及其应用研究--以1997年西藏玛尼7.9级地震区域形变测量为例

雷达差分干涉测量是一种最新的大地形变测量遥感技术方法,本文对差分干涉技术的原理及具体实现进行了深入研究,并讨论了差分干涉测量结果的误差.利用欧洲空间局1996年4月15日获取的ERS-1数据和4月16日获取的ERS-2数据,以及1997年12月2日获取的ERS-2数据,应用差分干涉测量技术对发生于1997年11月8日的西藏玛尼地震进行了提取区域形变场的应用研究,从得到的变化检测条纹图中可识别出地表破裂带,还可定量推算震中周围和两条断裂带附近的变形情况,差分干涉测量结果与地面调查符合得很好.     

提取地震同震形变场特征的最佳雷达波段问题研究

利用弹性形变模型, 并以走滑型断层为例, 对地震同震形变场进行了模拟研究。 结合雷达成像几何关系, 得到干涉雷达的LOS向形变场以及干涉相位条纹图像。 对比分析了不同雷达波段(X波、 C波、 L波)以及不同雷达入射角(30°和50°)的同震形变场和干涉相位图。 结果表明, 同一雷达波段不同雷达入射角的LOS向形变场会有不同的效果, 入射角逐渐增大时, LOS向形变量随之增大, 体现了发震断层走滑特征; 采用不同雷达入射角, 其探测到的地表形变垂向分量和水平分量信息量不同, 表明入射角较大或较小时能够分别较好地对地表形变的水平变化和垂向变化进行监测; 同一雷达入射角, 不同雷达波段的干涉相位条纹有明显区别, 波长越短, 干涉条纹越密集, 容易出现去相关现象, 长波长的干涉相位条纹比较清晰。 从某种意义上讲, 尽管长波长的干涉相位反映地表形变信息的细节会减少, 但受噪声影响较小, 抗干扰性能较短波长的强, 干涉像对相干性高。 地表形变表现为长期缓慢的小幅度渐变性形态, 也表现为短期迅速的大面积突发性形态。 不同雷达波长探测到的地表形变有效信息不同, 为了更好地利用干涉雷达对各种尺度和不同变化速率的地表形变进行监测, 需要在满足精度的条件下尽可能地选择波长较长及多样化的雷达数据。     

雷达差分干涉测量技术及其应用研究——以1997年西藏玛尼7.9级地震区域形变测量为例

雷达差分干涉测量是一种最新的大地形变测量遥感技术方法，本文对差分干涉技术的原理及具体实现进行了深入研究，并讨论了差分干涉测量结果的误差。利用欧洲空间局1996年4月15日获取的ERS-1数据和4月16日获取的ERS-2数据，以及1997年12月2日获取的ERS-2数据，应用差分干涉测量技术对发生于1997年11月8日的西藏玛尼地震进行了提取区域形变场的应用研究，从得到的变化检测条纹图中可识别出地表破裂带，还可定量推算震中周围和两条断裂带附近的变形情况，差分干涉测量结果与地面调查符合得很好。     

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

利用欧空局Sentinel-1A SAR数据,重建了2020年1月19日新疆伽师县MS6.4地震InSAR同震形变场.以升降轨InSAR形变场联合约束,反演了发震断层参数与同震滑动分布.结果 表明,地震同震形变场(长轴走向近EW向)发生在柯坪塔格褶皱带与奥兹格尔它乌褶皱带之间的区域内,升降轨观测的视线向(LOS)形变量...     

2022年1月8日青海门源MS6.9地震GNSS同震形变场及其断层滑动分布

2022年1月8日在青海省海北州门源县发生M_S6.9地震,本次地震是继2016年门源M_S6.4地震后冷龙岭断裂周边发生的又一次强震。确定本次地震的破裂分布对分析该地区震害风险具有重要意义。通过收集震中及周边12个GNSS连续站点和震后加密观测的17个流动站点观测资料,获取了震中100 km范围内29个测站的GNSS静态同震形变场,并以此为约束反演了本次地震同震滑动分布。结果显示,近场GNSS观测到的最大形变量可达1.3 m。反演的最优破裂模型显示该地震主破裂区深度位于0～10 km,滑动破裂出露地表,最大滑动量为4.07 m,地震矩释放能量约1.1×10¹⁹ N·m,对应矩震级M_W6.7。门源地震破裂至地表是造成该地区基础设施破坏的直接原因。