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

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.     

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

     

汶川Ms8．0级地震InSAR同震形变场观测与研究

利用InSAR技术,采用地震前后日本ALOS/PALSAR数据,提取了2008年5月12日四川汶川地震7个条带的地表同震形变场.每个形变条带南北向500 km,东西向70 km,7个形变场覆盖了映秀镇、都江堰、茂县、北川、平武和青川.结果显示,此次地震地表破裂带在北川-映秀断裂带上.地表破裂带从汶川县映秀镇西南震中附近一直到青川县苏河北侧,全长约为230 km.发震断层西北盘为抬升盘,南东盘断层附近,仍然表现为隆起区,显示出以逆冲为主的断层性质.在汶川县映秀镇西侧震中区,最大相对卫星视线向形变量达260 cm,如果全部换算成垂直形变,则两个区域的垂直相对形变达3.3 m.从北川至平通一带,卫星视线向形变范围在120～180 cm的隆起带.其中,擂鼓镇隆起形变范围170～180 cm,换算成垂直形变约在2.2～2.3 m之间.在青川苏河北附近,有70～80 cm范围的隆起形变.在雅安、峨眉山一带,以及射洪至重庆北侧一带有大范围沉降区.在重庆及其南侧区域有幅度在20～30 cm小范用隆起.由青川向东至广元、宁强一带,有形变幅度在60～70 cm的隆起区.整个形变场影响范围较大,四川盆地出现了不同程度的地表形变.     

2010年青海玉树地震同震-震后形变场特征及演化过程

基于星载合成孔径雷达差分干涉测量技术(DInSAR)和4期ENVISAT/ASAR雷达数据,获得了不同时间基线的三个同震干涉形变场和两个震后干涉形变场,并对这五个在时段上互有重叠的形变场进行了综合分析.结果表明,玉树地震同震形变场为围绕发震断层NW展布的椭圆形干涉条纹,覆盖范围约89 km×59 km.断层运动性质为左旋走滑.两盘最大视线向相对形变量至少达45 cm,最大形变出现在结古镇附近.时间基线不同的同震形变场总体上基本一致,但两盘最大相对形变量和局部形变存在差异.震后时间较长的干涉对反映的最大形变量反而减小;在震后时间较短的干涉对上于结古镇西南侧观测到的局部形变,在震后时间较长的干涉对上却没有出现.分析认为在形变量最大的结古镇附近可能出现了震后快速弹性回弹,导致随震后时间延续,形变量反而减小的现象.玉树地震震后形变主要出现在断层附近、震后不久的时段内,形变量在8 cm以下,具有与同震方向一致和相反的两种震后形变方式.在结古镇西南观测到一个与同震形变相反的局部沉降,应为震后弹性回弹.在微观震中处的断层附近观测到与同震方向一致的震后形变,可能是震后余滑.通过对地震前后不同时间基线的多个干涉对的联合对比分析,可以在一定程度上区分同震形变与震后形变,更好地研究地震引起的变形过程,特别是地震断层附近短期震后形变场的演化过程,为进一步研究断层带的岩性特征、物理力学及运动特性提供约束.     

汶川MS8.0级地震InSAR同震形变场观测与研究

利用InSAR技术,采用地震前后日本ALOS/ PALSAR数据,提取了2008年5月12日四川汶川地震7个条带的地表同震形变场.每个形变条带南北向500 km,东西向70 km,7个形变场覆盖了映秀镇、都江堰、茂县、北川、平武和青川.结果显示,此次地震地表破裂带在北川—映秀断裂带上.地表破裂带从汶川县映秀镇西南震中附近一直到青川县苏河北侧,全长约为230 km.发震断层西北盘为抬升盘,南东盘断层附近,仍然表现为隆起区,显示出以逆冲为主的断层性质.在汶川县映秀镇西侧震中区,最大相对卫星视线向形变量达260 cm,如果全部换算成垂直形变,则两个区域的垂直相对形变达3.3 m.从北川至平通一带,卫星视线向形变范围在120～180 cm的隆起带.其中,擂鼓镇隆起形变范围170～180 cm,换算成垂直形变约在2.2～2.3 m之间.在青川苏河北附近,有70～80 cm范围的隆起形变.在雅安、峨眉山一带,以及射洪至重庆北侧一带有大范围沉降区.在重庆及其南侧区域有幅度在20～30 cm小范围隆起.由青川向东至广元、宁强一带,有形变幅度在60～70 cm的隆起区.整个形变场影响范围较大,四川盆地出现了不同程度的地表形变.     

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

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

汶川地震断层滑动及形变演化模拟研究

2008年5月12日,四川汶川Mw7．9地震发生在青藏高原东缘龙门山断裂带的北川和彭灌断层上,地表破裂最大右旋走滑4．9m。地震发生后,中国地震局立即组织了GPS、In—SAR、地质等监测。中国地震局地震研究所在完成GPS监测任务的同时,复测了四川GPS控制网A级和B级点、三角点,获得464个高精度GPS同震形变。并根据其中13个1HzGPS连续站数据获得震时动态形变序列,用震后GPS监测网获得震后8个月47个站的震后形变序列,用日本宇航局的ALOSPALSAR数据获得InSAR同震形变场。     

高级时序InSAR地面形变监测及地震同震震后形变反演

差分干涉测量技术(differential interferometry technique,D-InSAR)是在InSAR技术基础上发展起来的,主要应用于测量地表微小形变。该技术在监测地震形变、火山活动、冰川运动、城市地面沉降以及山体滑坡等领域得到广泛应用,但实际研究中存在很多局限性,例如轨道参数误差、地形数据误差、大气延迟误差、时空去相关引起的相位解缠误差以及系统噪音误差等。这些因素成为InSAR高精度形变探测应用的主要限制。InSAR时序分析方法是目前解决常规D-InSAR处理过程中各种精度制约因素的主要方法之一。本文以InSAR时序分析中几个关键技术问题为研究目标,开展InSAR时序关键技术研究。本文取得的创新点如下:(1)在传统相干点提取方法的基础上,提出了一种新的提取方法:满秩相干点。该方法将干涉冗余网络图中节点和边结合起来构建满秩矩阵,并依此作为相干点提取指标。实验结果表明,满秩相干点提取法不仅能有效提高相干点的空间分布密度,而且能保留相干点的最优测量精度。(2)将大气延迟相位分成3种分量,即长波长大气相位、短波长大气相位以及地形相关的大气相位。在此基础上,利用网络法分别对3种类型大气延迟相位进行校正。网络法可以估算出每个时刻的大气延迟误差,再重构每一幅干涉图的大气延迟误差。实验结果分析表明,采用网络法能精确模拟出每个时刻大气延迟误差和每幅干涉图中的大气延迟相位。(3)离散点相位解缠误差是时序分析过程中主要误差源之一,直接影响后续时序分析结果的精度。本文提出一种利用网络闭合环残差方法,对离散相干点相位解缠误差进行检测和校正。该方法可以有效识别出离散相干点解缠相位中相位跳变;同时,本文也介绍了离散相干点解缠相位跳变校正的方法。(4)基于相干矩阵满秩条件,提出一种新的InSAR时序反演策略。利用干涉冗余网络图中各条边和节点之间的关系,采用最小二乘方法解算各个时刻点的形变特征,以有效解决传统InSAR时序反演过程中的秩亏问题。通过Bam地震震后时序形变场提取,在羊八井电站地热开采地下水引起地面沉降监测以及北京市地面沉降监测等几种典型应用中,验证了该方法的有效性和精确性。在关键技术突破的基础上,本文利用提出的新方法,获取了2008年10月6日西藏当雄MW6.3地震同震形变和震后地表形变时序,研究发震断层参数和震后运动过程的动力学机制。利用同震形变反演得到断层几何参数及滑动分布模型为先验知识,以InSAR时序方法获得震后形变时序资料为约束,分别利用震后余滑模型和震后黏弹性松弛模型,研究当雄地震震后断层的运动过程,同时对青藏高原南部中低地壳或者上地幔的黏性系数进行估计,并给出了最优拟合解。     

InSAR干涉形变场远程模拟系统实现与地震地表形变场模拟算例分析

利用InSAR(Interferometric SAR)干涉测量技术可以获得地表形变场(视线向, LOS), 将其与弹性半空间中一定断层模型模拟计算出的地表形变场进行正演分析, 以便获取发震断层的几何学和运动学参数, 是一种典型的前向模拟模式。 该计算模式由于模拟计算程序调试、 安装过程的复杂性, 使得在每台微机上安装此类程序不仅费时费力, 且浪费大量的计算资源。 文中首先介绍了实现InSAR干涉形变场模拟模型由单机推广到Internet/Intranet上进行远程计算最终建立“InSAR干涉形变场远程模拟系统”的过程, 并基于此系统平台对1997年玛尼M_S7.9地震的InSAR地表形变场进行了模拟试算、分析等。     

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, 与地震波形反演结果一致。     

汶川MS 8.0地震InSAR同震形变场特征分析

采用7条地震前后日本ALOS/ PALSAR整轨数据,利用D-InSAR技术提取了2008年5月12日四川汶川地震500km×450km的地表连续同震形变场.形变场覆盖了金川-石棉、黑水-乐山、松盘-彭山、南坪-简阳、康县-重庆所有区域,包括了理县、汶川、茂县、北川、青川等重灾区域.结果显示,整个形变场影响范围较大,四川盆地出现了不同程度的地表形变.发震断层附近的非相干区域显示此次地震地表主破裂带在北川-映秀断裂带上,可以追踪出地表破裂带从汶川县映秀镇西南震中附近一直到青川县苏河北侧,全长约为230km.在发震断层西南段汶川至茂县一带,非相干条带的宽度明显大于其它段落,这与彭县-灌县断裂(前山断裂)的都江堰-安县段地表破裂段密切相关,地表破裂带长约70km.在远离发震断层的区域,西北盘总体表现为抬升,东南盘表现为沉降.但在发震断层附近,断层两侧均表现为局部抬升,且沿断层形变量分布很不均匀,表现出较强的分段性,显示出发震断层以逆冲为主的断层性质.从美国哈佛大学(Harvard)、美国地质调查局(USGS)与美国国家地震信息中心(NEIC)、中国地震台网中心(CENC)所给出的震中位置和发震时刻的差异来看,也反映出汶川地震破裂过程是多点破裂过程.在汶川县映秀镇西侧震中区,最大相对形变量达260cm,如果全部换算成垂直形变,则两个区域的垂直相对形变达3.3m.在雅安、峨眉山一带约有35cm的沉降区,在重庆及其南侧区域约有25cm的小范围隆起.     

利用D-InSAR技术研究西藏改则地震同震形变场

针对2008年1月9日MW6.4西藏改则地震和2008年1月16日的MW5.9余震,通过两通(2-pass)加外部DEM差分干涉处理技术(D-InSAR),提取了地震区域2次地震累积的视线向(LOS)同震形变场。结果表明:发震断层均为正断层,位于依布茶卡-日干配错断裂端点附近。主震发震断层走向为N30°E,余震发震断层走向为N21°E,两断层距离约7km;在影像上主震发震断层有造成地表破裂的痕迹,余震未见地表破裂的痕迹;这次地震造成的同震形变场长约30km,宽约20km,主震断层上盘和下盘视线向最大形变量分别为39.2cm和11.2cm,两盘相对位错达50.4cm,余震造成的视线向形变量为9.4cm     

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

汶川M_S 8.0地震InSAR同震形变场特征分析

采用7条地震前后日本ALOS/PALSAR整轨数据,利用D-InSAR技术提取了2008年5月12日四川汶川地震500km×450km的地表连续同震形变场。形变场覆盖了金川—石棉、黑水—乐山、松盘—彭山、南坪—简阳、康县—重庆所有区域,包括了理县、汶川、茂县、北川、青川等重灾区域。结果显示,整个形变场影响范围较大,四川盆地出现了不同程度的地表形变。发震断层附近的非相干区域显示此次地震地表主破裂带在北川-映秀断裂带上,可以追踪出地表破裂带从汶川县映秀镇西南震中附近一直到青川县苏河北侧,全长约为230km。在发震断层西南段汶川至茂县一带,非相干条带的宽度明显大于其它段落,这与彭县-灌县断裂(前山断裂)的都江堰—安县段地表破裂段密切相关,地表破裂带长约70km。在远离发震断层的区域,西北盘总体表现为抬升,东南盘表现为沉降。但在发震断层附近,断层两侧均表现为局部抬升,且沿断层形变量分布很不均匀,表现出较强的分段性,显示出发震断层以逆冲为主的断层性质。从美国哈佛大学(Harvard)、美国地质调查局(USGS)与美国国家地震信息中心(NEIC)、中国地震台网中心(CENC)所给出的震中位置和发震时刻的差异来看,也反映出     

利用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和应变张量估计2016年MW7.0熊本地震同震三维形变场

利用合成孔径雷达(Synthetic Aperture Radar,SAR)影像提取与地质活动相关的三维地表形变场,对深入理解地质灾害的形成机制及其潜在灾害风险评估非常重要.目前,利用SAR影像的同震三维形变场提取主要利用单个像素点的多次观测构建观测方程,然后基于加权最小二乘(Weighted Least Squares,WLS)方法分解从而获得同震三维形变场,因此该方法缺乏对相邻像素点空间相关性的约束.考虑相邻同震位移点的应力连续性,研究学者提出了顾及大地测量应变张量和卫星形变观测的SAR同震三维形变场方法(Extended Simultaneous and Integrated Strain Tensor Estimation from geodetic and satellite deformation Measurements,ESISTEM).本文以2016年MW7.0熊本地震为例,收集了覆盖此次地震的ALOS-2卫星升降轨影像,利用传统差分InSAR(DInSAR)方法和子孔径雷达干涉测量(Multiple Aperture InSAR,MAI)方法分别对升降轨SAR影像对进行处理,得到视线向(LOS)形变和方位向形变,最后利用ESISTEM方法获取此次地震的三维同震形变场.此外,利用GPS和野外考察观测对本文的三维形变场结果进行结果精度分析.研究结果表明,与传统WLS方法相比,ESISTEM方法不仅能有效抑制奇异像素点对形变结果的干扰,同时对近断层的失相干信号能进行较好的恢复,更有助于解释地表破裂区的地震形变特征和掌握地震发生机制.本文确定的三维同震形变场结果显示主形变区发生在Futagawa断层中部和Hinagu断层最北端,最大水平位移为2m,抬升为0.55m.断层破裂以NE-SW走向的右旋走滑为主兼有部分正断成分.应变张量分析表明发震断层处受到了明显的收缩力和剪切力的作用.     

汶川Ms 8.0地震InSAR形变观测及初步分析

2008年5月12日汶川MS8.0地震发生后,通过干涉雷达方法(InSAR)获得了覆盖整个灾区全面的高精度形变观测资料(约500km×350km),为研究该地震的发生机理提供了重要的观测依据。通过分析InSAR成像几何和ALOS PALSAR数据的特点,初步获得了对汶川地震的定性认识,为相关研究提供参考,也为进一步深入定量反演奠定基础。此次地震中InSAR卫星观测的应用,显示了卫星大地测量手段对于地震科学研究的重要意义     

基于敏感性迭代拟合算法反演汶川Mw 8.0地震InSAR同震形变场

基于单一断层模型,运用敏感性迭代拟合算法反演了汶川地震InSAR同震形变场,获得了断层滑动分布及部分震源参数.结果表明,倾角线性变换的单断层模型模拟的同震形变场与InSAR形变场吻合较好,且残差较小,平均残差为0.11 m;反演的滑动分布主要集中于地下0～20km,其中,汶川地区在震源深度附近有较大的滑动量,10～15 km深度最高可达7 m,地表滑动量却相对较小,平均值仅为2.5 m,平均滑动角均值约为121°;北川地区最大滑动量可达到10 m,平均滑动角均值约为109°;而青川地区10～15 km深度滑动量最高也达8 m,平均滑动角均值约为135°;滑动矢量图显示,沿SW-NE走向断层面的滑动方向以强烈的逆冲为主,兼有一定右旋走滑分量.反演矩张量为1.0×1021 N·m,矩震级达Mw8.0.     

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