菲律宾棉兰老岛2019年4次MW>6.0地震断层运动模型及库仑应力传输研究

基于雷达干涉测量技术,利用ALOS-2、Sentinel-1卫星升降轨雷达影像,获得2019-10~12发生在菲律宾棉兰老岛的4次M_W>6.0地震的同震形变场,并以此形变结果为约束,反演得到4次地震的断层运动模型。综合分析发现,此次地震序列由3条断裂的破裂引起,其中2019-10-16和2019-10-31的2次地震为同一发震断裂,2019-10-31地震断层破裂区域位于2019-10-16地震断层破裂的东北延伸段,最大滑动量约为1.1 m,约为2019-10-16地震最大滑动量的2倍。2019-10-29地震由一条独立断层破裂引起,断层最大滑动量约为2.0 m。2019-12-15地震由一条东北向倾斜断层破裂引起,断层最大滑动量约为3.0 m。此外,2019-10-16地震引起2019-10-29地震显著滑动区明显的正向库仑应力传输;而2019-10-29地震显著增加了2019-10-31地震震源区域的库仑应力;前3次地震对2019-12-15地震孕震断层的库仑应力传输总和为负值,说明静态库仑应力传输可能不是此次地震触发的主要诱因。     

汶川地震同震形变场的GPS和InSAR邻轨平滑校正与断层滑移精化反演

为克服InSAR观测汶川地震同震形变场的邻轨不连续问题,提出联合GPS观测值与邻轨平滑约束的同震位移校正方法,采用GPS观测形变去除PALSAR轨道误差引入的残留平地相位,基于形变平滑条件校正邻轨干涉相位的不连续性.ALOS/PALSAR干涉处理结果表明,校正后同震形变场的准确度与平滑性得以显著提高,InSAR高相干点残差达3.6 cm,校正后精度提高约60%,低相干点精度提高约40%,校正后形变场的邻接平滑因子标准差减小约33%,验证了轨道误差校正与邻轨平滑约束方法的准确性与可靠性.进一步基于弹性半空间位错模型的断层滑移反演结果表明,断层滑移主要分布于映秀、北川和青川地区,集中于地壳深度0~16 km范围,最大滑动量（位于北川县城）约为9.0 m,GPS反演模型残差为5.5 cm,InSAR反演模型残差达9.2 cm,InSAR反演精度约有30%的显著提高,由模型反演计算得到的地震矩为8.0469×10²⁰ N·m.     

2021年云南漾濞MS6.4地震动态应力对后续余震活动的触发作用

选取IRIS远震台站波形数据,反演了云南漾濞M_S6.4地震震源破裂过程,计算了断层破裂在近场产生的动态库仑破裂应力变化,并讨论了主震对近场余震活动的动态应力触发作用。结果显示:动态库仑应力演化过程与震源破裂特征反演结果一致,其大小分布与地震序列分布的疏密程度也具有较好的相关性。主震产生的静态和动态库仑破裂应力均促进余震的发生,但相比静态应力,余震位于库仑破裂应力正值区域的比例提高了21%,余震与动态库仑应力变化的正负区域有更好的一致性,动态应力能更好地解释震后余震分布的空间特征。垂直于地震序列主干10 km处出现小震丛集,该现象可能是由主震产生的动态库仑破裂应力占主导作用所致。定量分析主震对余震的动态应力触发结果显示,主震后一周内M_S4.0以上的8次余震接收点均受到了动态库仑破裂应力的触发作用。      

汶川地震SAR形变场电离层影响校正与断层地表破裂线自动提取

针对汶川地震PALSAR观测的方位向形变场存在严重电离层影响,导致准确提取断层地表破裂线困难的问题,本文使用一种基于频率域的电离层影响校正方法,从方位向形变场频域数据中提取电离层信号对应的频域信息,通过建模获得空间域电离层信号,并与原始形变场进行差分运算完成电离层影响的校正,还原了汶川地震断层近场沿SAR方位向的真实地表形变场.以校正后的方位向形变场作为数据源,采用非极大值抑制二维梯度算法,准确提取汶川地震北川—映秀断层和灌县—江油断层的地表破裂线.为验证提取结果的可靠性,引入野外地质调查数据进行对比分析,定量计算两条破裂线间的距离差异,结果表明基于SAR形变场自动提取的断层破裂线与野外地质调查数据具有较高的一致性,利用SAR观测形变场可作为调查地震断层地表破裂的重要补充手段.     

2021年漾濞MW6.1地震同震形变提取与断层源模型反演研究

基于D-InSAR技术,收集Sentinel-1卫星升、降轨及精密轨道数据,结合外部DEM,成功提取2021-05-21云南省大理州漾濞县M_W6.1地震的同震形变场。结果显示,升轨卫星视线向最大形变量级约为6.0 cm,降轨卫星视线向最大形变量级约为7.9 cm。以升、降轨InSAR观测结果为约束,对漾濞地震的断层几何参数及滑动分布进行联合反演发现,此次地震的最优发震断层走向为136.6°、倾角为83.1°,断层破裂主要集中在地下2~12 km深度处,最大滑动量约为0.45 m,位于地下7 km深度处,断层在近地表未出现大面积显著滑移现象,表明此次地震未破裂至地表。     

利用SAR影像配准偏移量提取地表形变的方法与误差分析

单一的InSAR观测技术可提取地表沿雷达视线方向(LOS)上的一维位移,而利用SAR影像配准过程中的同名像素偏移量可提取地表沿雷达方位向(近南北向)与距离向(近东西向)的二维形变场,与LOS方向的一维形变形成优势互补。本文在分析SAR影像配准偏移量提取地表形变场方法的基础上,推导建立了雷达方位向与距离向形变提取的误差模型,探讨了该方法提取地表形变的主要误差源。以Bam地震ASAR影像和玉树地震PALSAR影像为数据源,分别开展同震形变场的提取与误差分析试验,结果表明,基于SAR像素配准偏移量提取同震形变场的精度主要受匹配窗口尺寸与过采样因子影响,形变提取误差随匹配计算窗口的增大而减小,形变提取精度随过采样因子的增大有适量提高,地形起伏效应在高差较大的SAR影像区域中表现较为显著。     

应用二阶Harris算子的高铁轨道近景影像特征点提取

由于轨道影像具有较高的灰度相似性,传统一阶Harris算法提取轨道影像特征存在时间效率低、特征点聚簇等负面问题,文章针对已有Harris方法的不足,提出了一种基于图像分块自适应阈值的二阶导数Harris特征点检测算法,采用自相关矩阵兴趣值识别轨道图像特征点,并分别对有砟与无砟铁路轨道影像开展特征点提取。实验结果表明,二阶Harris算子在轨道近景影像特征点提取方面具有较高的计算效率,所提取出轨道影像特征点满足均匀分布的空间特征,有效避免了传统方法特征点的聚簇现象。     

GPS与PS-InSAR联网监测的台湾屏东地区三维地表形变场

联合高精度的GPS水平位移观测和高密度的PS-InSAR雷达视线位移测量,实现地表三维形变的精确反演.本文在准确计算卫星轨道方位角基础上,使用GPS观测位移与星载雷达LOS方向形变的投影转换模型,将雷达LOS方向形变转换为垂直方向位移,并基于地面GPS与SAR影像PS目标联合构建形变监测网,采用参数平差算法估计区域地表形变场.以地质构造活动极其活跃的台湾岛及其西南屏东高雄地区为例,联合屏东地区48个GPS监测台站与雷达PS目标,监测该地区从1995-1999年间由于板块构造挤压运动和地下水抽取导致的三维地表形变.结果表明,该地区年均水平位移量为向西30~50 mm/a,高雄沿海地区发生明显的逆时针西偏南的逐渐增大的水平位移;垂直位移为屏东平原南部呈现-10 mm/a~-15 mm/a的地面沉降,而平原北部和高雄地区呈现约+5 mm/a~+10 mm/a的地面抬升.     

福建省山洪灾害危险性评价

福建省山洪灾害频发,造成的损失十分严重,开展山洪灾害危险性评价,对于防灾减灾具有重要意义。本文结合福建省山洪灾害的时空分布特征,从致灾因子和孕灾环境出发,选取多年年均降雨量、坡度、NDVI、土壤类型、土地利用类型等10类影响因子。利用CF模型和Logistic回归模型确定各因子与其分类级别的权重,应用皮尔逊相关系数法对因子间独立性进行检验。并利用历史山洪点与ROC曲线对评价结果进行验证。研究结果表明:(1) 100年一遇最大24 h降雨量、土地利用类型及地形起伏度是影响福建省山洪灾害的三个主要因素;(2)极低危险区、低危险区、中危险区、高危险区、极高危险区占全省总面积的比例分别为63. 1%、14. 2%、4. 5%、5. 9%、12. 3%;(3)福建省东部、南部、东南部的沿海地区以及北部、西部的个别低山地区处于高危险或极高危险区;(4)在仅占研究区面积12. 3%的"极高危险区"的较小范围内,实际发生山洪点数占山洪点总数的比例高达65. 8%,训练样本和验证样本的ROC曲线下面积AUC值也分别达到0. 876和0. 891,两种方式都说明山洪危险性评价结果较好。可为福建省山洪灾害防治提供科学依据和支持。     

2017年伊朗MW7.3地震的震后余滑分布及其对2018年MW6.0地震的触发影响

收集了覆盖监测区域的Sentinel-1卫星雷达影像,利用短基线集干涉测量技术提取了2017年伊朗萨波尔扎哈布M_W7.3地震后283天的地表时序形变,通过二步法反演得到其震后余滑分布,之后采用差分干涉技术获取了2018年发生于同一地区的贾万鲁德M_W6.0地震的同震形变场,并将反演所得的发震断层参数作为应力计算的接收断层参数,来分析2017年M_W7.3强震及其震后活动对2018年贾万鲁德M_W6.0地震的触发影响。结果表明:萨波尔扎哈布地震的震后形变主要由孕震断层面的余滑运动所致,震后283天余滑模型的累积滑移量达到0.7 m;2018年贾万鲁德地震的发震断层走向为355.6°,倾角为89.4°,同震断层破裂以右旋走滑为主,兼具部分正断层运动。本文所得的贾万鲁德地震断层平面上的库仑应力变化表明,2017年M_W7.3主震及其震后余滑对2018年M_W6.0地震的发生具有一定的触发效应,M_W6.0地震的发生可能与区域板块的活动性相关。