2016年台湾美浓MW6.4地震震源参数的InSAR和GPS反演

2016年2月6日台湾西南部高雄市美浓区发生了M_W6.4地震.本文结合ALOS2卫星升降轨、Sentinel-1A升轨SAR数据,采用两轨差分干涉技术获取了该区域的同震形变场,形变结果表明震中西北部以抬升为主,最大视线向形变量约为11.2 cm.基于均匀位错模型和多峰值粒子群（MPSO）算法,利用InSAR和GPS形变数据联合反演了美浓地震的断层几何参数,结果表明震源中心位于22.920°N,120.420°E,深度约12 km,发震断层长度约15 km,走向角307°,倾角16.5°,平均滑动角为51.5°,此次地震是以逆冲倾滑兼左旋走滑的破裂模式.利用格网迭代搜索法得到最优倾角为15.7°,GPS和InSAR最优权比为18：1,最优平滑因子为0.06.基于非均匀位错模型,利用非负最小二乘方法进行线性反演,结果显示最大倾滑和走滑量分别为51.7 cm和55.3 cm,对应矩震级为M_W6.38,略小于GCMT（M_W6.4）的结果.通过与已有文献的比较和对该区域断层构造的分析,发现美浓地震的发震断层为单一断层的解释更为合理,我们推测发震断层是位于左镇、后甲里等断层之间的一条东南-西北走向往东北倾斜的盲断层,并初步推测2010年M_W6.3甲仙地震也同该断层有关.

     

利用InSAR资料反演宁洱Ms6.4地震震源参数及库仑应力变化

首先利用ALOS PALSAR数据,通过D-InSAR技术获取2007-06-03云南宁洱MS6.4地震的同震形变场,然后基于Okada弹性半空间位错模型反演该地震的断层几何以及精细滑动分布,最后计算宁洱地震后周边断层的静态库仑应力变化。结果表明,形变主要集中在西盘,最大视线向形变量为51.6 cm;反演得到的震源位置为23.05°N、101.02°E,深度3 km,断层走向145°,倾向49.5°,平均滑动角153°,发震断层为NNW向普洱断裂,断层活动以右旋走滑为主,兼具逆冲分量;断层面最大滑动量为1.2 m,反演得到的震级为MW619。基于库仑应力场发现,磨黑断裂处于库仑应力增加区域,而2014年景谷地震位于负值区域。结合实地考察资料和反演结果表明,宁洱地震为浅源地震,但断层并未出露地表。     

基于有效约束的附不等式约束平差的一种新算法

不等式约束是客观实际中普遍存在的一种约束,但目前大地测量数据处理领域并没有成熟、完整并被普遍接受的处理理论和方法。首先简要总结附不等式约束平差的各种方法及其存在的问题。然后对现有测量平差中附有等式约束的平差模型进行扩展,提出一种新的处理附有线性约束(包括等式和不等式约束)的平差方法。该方法在有效约束概念下,通过库恩-塔克条件来寻找有效约束条件,把不等式约束平差问题转化为我们熟知的等式约束平差问题,因此实现解向量与观测向量之间的显式表达。最后,用一个数值算例验证新方法的可行性,同时算例分析表明:用等式约束代替有效约束或集成约束进行平差计算,能得到正确的平差结果,但得不到正确的精度评定结果。     

哨兵-2号光学影像地表形变监测:以2016年Mw7.8新西兰凯库拉地震为例

光学影像已经广泛地应用于地表形变监测研究。哨兵-2号光学影像作为一种新的对地观测数据,具有重访周期短、空间分辨率高、影像覆盖范围大以及数据免费等优点,因此该数据在地表形变监测上有广泛的应用潜力。本文以COSI-Corr软件包为数据处理平台,基于亚像素的频率域相关性匹配技术,处理多时相的哨兵-2号数据获取地表形变。本文选取2016年M_w7.8新西兰凯库拉（Kaikoura）地震覆盖区域为例,对哨兵-2号影像地表形变中存在的各种系统误差源进行了系统分析和改正,并提出了改进的均值相减法去除形变场中的卫星姿态角误差。另外,还对哨兵-2号4个10 m空间分辨率的波段（Band 2/3/4/8）中可用于地表形变监测的最佳波段进行了分析,统计结果显示Band8的地表形变监测效果最好。最后,利用哨兵-2号光学影像获取了2016年11月14日M_w7.8新西兰凯库拉地震的同震形变场;分析了沿地震主要断层的滑移分布,结果表明最大水平滑移量达10 m;并与同期Landsat8全色影像的同震形变监测结果进行对比分析,结果表明哨兵-2号结果精度更高。本文的研究成果可以为哨兵-2号光学影像的应用提供参考。     

附不等式约束平差的理论与方法研究

附不等式约束平差能充分利用已有的各种知识和信息,因而越来越得到关注,它在GPS数据处理、变形监测及方差估计等领域中都得到广泛的应用。目前,国内外对附不等式约束平差问题进行大量研究,文中主要综述附不等式约束平差问题的最新研究成果,系统介绍附不等式约束平差的理论,同时探讨附不等式约束平差的发展方向。     

InSAR部分地学参数反演

自然因素或人类活动可使地球表层或内部应力发生变化,进而产生灾害事件。获取灾害事件与灾害过程的关键地学参数对于准确理解灾变过程、科学解释灾变机制、正确拟定应对策略具有重要意义。InSAR技术已广泛应用于自然因素或人类活动导致的灾害事件与灾变过程的参数反演。本文首先介绍了InSAR卫星的发展和地表形变监测的基本原理;然后分析了InSAR地学参数反演在地震、火山活动、地下水抽取、矿山开采、冻土冻融、冰川运动、地下流体迁移等各类潜在致灾事件中的研究现状;最后分析了InSAR地学参数反演存在的主要挑战和问题。     

利用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级.     

基于BFGS法融合InSAR和GPS技术监测地表三维形变

虽然InSAR技术具有高精度、大范围和高空间分辨率等优点,但只能监测雷达视线方向上的一维地表形变;而GPS技术虽可以监测地表的三维形变,但其空间分辨率很低.本文针对融合InSAR和GPS技术监测地表高空间分辨率三维形变展开研究.首先证明了简单的局部最优化迭代算法就能求得综合InSAR和GPS监测地表形变速率的能量函数模型的全局最优估值.随后提出了利用BFGS局部最优算法反演最优的地表三维形变速率.该方法既能避免全局最优化算法计算复杂且难以收敛的问题,又能克服传统的解析法中数值计算不稳定的缺点.最后,通过模拟实验和美国南加州真实数据实验表明,该方法能够得到高精度的地表三维形变速率场.而且当观测或插值误差导致解析法误差较大时,BFGS方法仍能得到高精度、稳定的全局最优解.     

InSAR视角下准噶尔盆地农业灌溉区地表形变特征和演化规律

抽取地下水进行农业灌溉是导致地下水位快速下降的重要因素,而长期过度开采地下水往往会引发地面沉降灾害,这种现象在干旱和半干旱地区非常普遍。为了研究农业灌溉超采引发的地表形变特征和演化规律,本文以准噶尔盆地南缘、天山北麓地带为研究区域,利用SBAS-InSAR技术对2003年—2009年覆盖呼图壁县的ENVISAT/ASAR升、降轨数据进行处理,获取了该地区的地表形变场,并结合研究区的农业灌溉方式、水资源补给和季节变化等资料对地面沉降的时空变化特征进行分析,为水资源和农业可持续发展提供参考意义。实验表明,研究区内主要有两个沉降幅度较大的漏斗,且都位于农田区域。2007年以前,研究区地表没有显著形变,之后发生了较大量级的沉降。采用传统灌溉方式和时针式灌溉系统的农业区平均沉降速率最高分别可达50 mm/a和30 mm/a,前者在时间上呈线性变化,而后者具有显著的周期性变化特征。在冬季时,采用时针式灌溉系统的地区地面抬升量可达40 mm,远大于传统灌溉方式的农田区域,而夏季地面沉降速率可达200 mm/a。对研究区农业灌溉活动进行分析后发现,农业灌溉造成的地下水超采是该地区地面沉降的主要影响因素,其形变机制与季节变化具有较高的相关性,在灌溉活动休止期内地表形变取决于地下水的补给量。研究区内的形变特征和影响因素分析将为地下水资源的充分利用和农业的可持续发展提供有效的信息。     

佛山地铁沿线时序InSAR形变时空特征分析

佛山作为中国珠三角地区经济和城市化高速发展的城市,由于其脆弱的地质水文条件,长期遭受地面沉降灾害的影响。同时,该区域地铁作为缓解城市交通压力的重要工具,其施工和运行所导致的地面沉降也影响了人们的生命财产安全。但目前针对佛山地区相关的系统研究不多,对地铁沿线的沉降规律认识不足。利用Sentinel-1数据监测了2015-06至2018-09间佛山市的形变信息,结果表明,佛山市地表形变呈零星分布,未出现大范围的沉降漏斗,形变速率为-20~5 mm/a,局部区域的沉降速率超过-30 mm/a。地面沉降主要与不稳定的地质结构、地下水抽取和局部区域工程施工有关。基于获得的形变结果,对佛山市地铁沿线的形变情况进行了研究,并对2018年佛山市地铁坍塌事故路段的沉降情况进行了详细分析,阐述了在空间分布上地铁沿线沉降差异的成因,并在时间上对地铁沿线的形变进行了模型参数反演。研究工作为今后当地政府开展地表形变普查、沉降灾害预警提供了参考,并为地铁正常运行与维护的安全监测提供了理论依据。