利用InSAR及其差分技术获取仲巴地震的同震形变场

本文简要介绍了InSAR及其差分技术的基本原理,并以仲巴地震为例,采用已知DEM的二路差分法获取了仲巴地震的同震形变场。结果表明:仲巴地震为是一个以倾滑为主的浅源正断层型地震,断层面走向近乎南北方向。     

多视角InSAR数据解算2017两伊地震三维同震形变场

2017两伊地震是自1900年以来发生在扎格罗斯山脉的最大地震,为了研究此次地震引起的同震形变场,利用覆盖同一地区的3对Sentinel-1A升降轨数据分别进行两通差分DInSAR处理,得到了研究区3个视线向的地表同震形变场,通过直接解算法重建了研究区的三维同震形变场。试验表明：3种视角的升降轨视线向上升与沉降总体趋势基本一致;联合多个视角的观测结果可以实现三维形变场的重建;根据地表视线向和三维同震形变的特征以及地质构造背景推测了发震断层很有可能为扎格罗斯山前断层。     

基于多源InSAR数据的三维地表形变解算方法研究

随着不同成像参数SAR卫星不断发射升空以及新的InSAR技术不断出现,同一地区可以获取多尺度、多源、异质的InSAR地表形变监测数据。将这些多源InSAR地表形变监测数据作为观测量,实际三维地表形变作为未知参数,建立模型并解算,实现东西、南北和垂直向地表形变监测,为地质灾害预报及监测提供更全面的基础资料。本文以InSAR三维地表形变解算为核心研究内容,针对大梯度、轨道和地形、粗差等因素对解算精度的影响,利用现代测量数据处理理论和方法,消除或抑制上述因素导致的各类噪声,实现多源InSAR数据的三维地表形变高精度解算。     

日本2011 Tohoku-Oki Mw 9.0级地震的同震形变及其滑动分布反演:GPS和InSAR约束

利用覆盖整个日本东北部的3个条带的Envisat/ASAR降轨数据和6个条带的ALOS/PALSAR升轨数据,通过二通差分干涉处理,首先得到了2011年日本Tohoku-Oki Mw 9.0级地震的初步同震形变干涉结果,然后,利用GPS同震观测值对InSAR结果进行校正和基准统一。在此基础上,采用弹性半空间矩形位错模型,联合GPS和InSAR观测结果对同震滑动分布进行反演,获取了发震断层的断层滑动分布。结果显示,此次地震的滑动分布主要发生在40～50km深度范围内,最大滑移量为50.3m,释放的能量为3.20×1022 N·m(相当于Mw 8.94级)。     

利用InSAR研究西藏拉孜地震同震形变

采用欧空局ERS—1卫星雷达数据获取了1993年Mw 6.2级西藏拉孜地震同震形变场,并采用弹性位错模型反演了该地震的断层参数。结果显示该地震产生的最大地表位移约为12 cm,断层倾角为42°,不支持该断层为低角度正断层的论断。     

2008年大柴旦M_w6.3级地震的InSAR同震形变观测及断层参数反演

2008年11月10日,青海省大柴旦地区发生了Mw6.3级地震。本文利用EnviSat卫星升降轨SAR数据和差分干涉测量技术提取了同震形变场,基于均匀位错模型反演确定了地震断层参数,然后利用格网迭代搜索法确定了较优断层倾角,同时基于非均匀位错模型获得了精细滑动分布。结果表明,地震使得上盘区域沿降轨、升轨视线向分别产生最大约8.5cm、10cm的抬升;较优断层倾角为47.9°;地震滑动未延伸至地表,主要发生在地下8.2~23.7km深度范围内,最大和平均滑动量为0.5m和0.19m,平均滑动角为104.9°。反演的地震矩为3.74×1018 N·m,矩震级为Mw6.35。     

基于InSAR技术的地表三维形变获取方法综述

对近年来提出的基于InSAR技术的地表三维形变获取方法进行了综述,并对该技术的发展进行了展望。对一些重要的知识背景和理论方法进行了简单介绍,并针对每种方法列举典型案例,对其适用条件进行了简单总结,为相关领域的研究者提供了最新研究进展和技术参考。     

芦山Ms 7.0级地震震前及同震地表形变

对四川地区GNSS观测数据进行处理,获得如下结果:①汶川Ms8.0级地震之前龙门山断裂带及其以北近150km的范围内几乎不存在可分辨的右旋走滑活动和压性形变,四川盆地以西-鲜水河安宁河断裂带以东约100km的范围内几乎也不存在可辨的形变;②汶川Ms8.0级地震之后芦山Ms7.0级地震之前,龙门山断裂带中、北段存在约5mm/a右旋走滑活动,芦山震源及附近地区的三维形变一直处在闭锁状态,鲜水河安宁河断裂带及其以东的左旋形变有所增强;③芦山地震可变的同震形变场基本分布于以震源为中心的数十km范围内,大致以震中为界东侧呈右旋形变,西侧呈左旋形变;④距震中较近约12km的LS05震时\     

基于GIS的InSAR结果分析方法及在汶川Mw7.9级地震同震解释中的应用

利用GIS对InSAR结果进行分析可以更全面地解释地表形变机制。介绍了利用GIS对SAR影像进行裁剪、切割和拼接处理、剖面和等值分析的方法。在此基础上,以汶川Mw7.9级地震为例介绍了基于GIS的InSAR结果分析方法在地震形变机制解释方面的应用:①根据InSAR距离向偏移量和同震形变场,利用GIS技术提取了汶川Mw7.9级地震发震断层的地表迹线;②将基于位错模型正演得到的视线向形变场与InSAR同震形变场进行比较,得到了汶川地区InSAR同震形变场的上盘和整体改正数,进一步借助GIS进行修正得到了更为真实的InSAR同震形变场;③利用2D/3DInSAR同震形变场、干涉纹图和剖面分析结果对汶川Mw7.9级地震同震形变特征进行了具体分析和解释。     

九寨沟Ms7.0级地震三维同震形变场提取

针对D-InSAR地震形变监测中存在LOS(line of sight)视线向模糊的问题,该文构建了一种融合升降轨不同视线向干涉测量数据获取三维同震形变场的方法。以2017年九寨沟M_S7.0级地震为例,基于研究区的Sentine1_1A数据和SRTM(1″)DEM数据,采用二轨差分的D-InSAR技术,融合升降轨LOS向以及自定义视线向的干涉测量数据进行联合解算,获取了研究区三维同震形变场,解决了单一轨道雷达LOS提取结果不能准确反映地表三维形变场的问题。实验表明,此次地震造成了震中附近一定范围地表的沉降和东南向的滑移,并通过震区形变剖面图和形变等值线图,分析了形变场的空间分布以及此次地震所造成地表断裂的位置,认为此次地震与塔藏断裂、虎牙断裂存在关联性,其运动形式为主动盘逆冲的走滑型地震。     

DInSAR技术监测玉树地震同震形变场的研究

近二十年来,合成孔径雷达干涉测量技术（DInSAR）作为一种监测地表形变的有效技术得到广泛应用,其视线向精度可达厘米级甚至毫米级。该技术尤其是对监测快速、激烈的地表形变更为有效,如地震、火山等。利用DInSAR技术,选用日本ALOS卫星PALSAR数据,获取了2010—04—14玉树Mw6．9地震的同震形变场。结果表明：地表至少发生过3次地表破裂,沿断层走向的形变分布范围远远大于垂直断层分布范围,形变分布特征以左旋走滑为主。     

联合InSAR与坡向约束的露天矿边坡三维形变监测

针对目前多轨道InSAR融合方法监测的南北向形变误差大,难以科学地识别露天矿边坡滑坡隐患的问题,本文提出了联合InSAR与坡向约束的露天矿边坡三维形变监测方法。首先,将露天矿边坡水平位移主要沿着坡向的先验信息作为约束条件。然后,根据SAR卫星形变投影关系,建立了附加约束条件的露天矿三维形变观测模型,并据此求解三维形变。最后,选用辽宁省阜新市新邱矿区作为研究区域,实现了仅基于升降轨InSAR数据的露天矿边坡三维形变监测。结果表明：本文方法解算的三维形变精度约为1.2 cm,仅占最大形变值的5.6%。本文方法有助于改善目前InSAR露天矿山三维形变监测精度不高的情况,提升露天矿山边坡滑坡隐患识别能力。     

顾及分段相似性的同震形变时间序列恒星日滤波优化算法

在常规恒星日滤波的基础上,提出一种顾及分段相似性的恒星日滤波优化算法。在构建滤波算子之前,首先确定用于构建滤波的坐标残差时间序列与地震当天坐标残差时间序列的分段相似性,然后根据相似性度量指标的大小确定权重后再进行滤波。研究结果表明:利用本文提出的顾及分段相似性的恒星日滤波优化方法能够显著改善单历元定位的精度和可靠性,克服常规滤波方法存在的缺陷,有利于恒星日滤波的无人工干预批处理。尤其是通用影像质量评价指数UIQI方法顾及了时间序列间的相关系数以及序列自身的均值和方差,比其他相似性度量更加科学合理,建议采用UIQI作为恒星日滤波优化算法的相似性度量。     

联合地震位错模型和ESISTEM方法提取地震同震三维形变场

高精度同震三维形变场对于研究地震变形模式、震源机制等具有重要意义。设计并实现了一种联合地震位错模型和扩展融合大地测量、卫星形变观测、应变张量估计(extended simultaneous and integrated strain tensor estimation from geodetic and satellite deformation measurements, ESISTEM)方法的新方法,以2021年Ms 6.4漾濞地震为例,利用哨兵1号A、B星(Sentinel-1A/B)升、降轨影像获得该地震的合成孔径雷达干涉测量形变场,利用地震位错模型正演得到的南北向形变分量进行约束,成功提取了该地震完整的同震三维形变场及应变场。结果表明,漾濞地震断层西南侧主要向西、向北运动,最大形变分别为4.8 cm、9.5 cm;东北侧主要向东、向南运动,最大形变分别为7.4 cm、4.6 cm;垂直向抬升、沉降的最大值分别为3.6 cm、3.4 cm;发震断层以右旋走滑运动为主,兼有少量正断分量;发震断层区域受到显著的膨胀、剪切和旋转作用。     

GNSS形变监测静态基线网中多基线联合解算两步法

GNSS(global navigation satellite system)静态基线网的解算方法通常采用单基线解算模式或多基线解算模式。这两种方法均有一定的局限性。本文综合两类方法的优点,提出了单基线和多基线联合解算模式的两步法。首先基于单基线模式利用多频率的载波相位观测值解算监测站间的双差整周模糊度,并使用严格的固定准则确定模糊度;然后利用多基线模式进行联合坐标和天顶对流层参数的估计,再结合最小二乘估计解算定位结果。选取国内某水电站静态形变监测数据为研究对象,进行基线解算方法测试和分析。实验结果表明多基线联合解算两步法在精度和可靠性上有较大优势,平面定位精度为亚毫米级,高程定位精度提高1～2倍;并可以明显改善基线网的解算速度,在高精度变形监测应用中具有可行性与通用性。     

DInSAR技术对地震同震形变场的研究

合成孔径雷达差分干涉测量（DInSAR）可用于监测厘米级或更微小的地表形变,以揭示许多物理现象,如地震形变、火山运动、大气变化、冰川漂移、地面沉降以及山体滑坡等。DInSAR作为一种新型的空间对地观测技术,具有不受时间和空间的限制、对地物具有一定的穿透性等传统测量所不可比拟的优势,已得到较为广泛的应用。简要介绍DInSAR技术的基本原理及其处理流程,以Barn地震为例提取Bam地震的同震形变场。     

基于PALSAR数据的2010年Darfield地震地表同震形变场提取

2010年9月3日（UTC时间）,新西兰南部小镇Darfield发生了强度为Mw7.1的地震,其主震对地表造成了严重破坏。为了对由地震引起的地表形变进行研究,本文利用二轨DInSAR技术,通过2景ALOS PALSAR影像的升轨数据,成功提取了该地震的地表同震形变场。从对实验结果的分析可知,该地震的断层为东西走向,断层北盘地表呈现沿雷达视线方向的下降,最大量约为1.5米;断层南盘则发生了沿雷达视线方向的抬升,最大量超过了2米。     

基于ALOS PALSAR数据进行雷达差分干涉提取玉树地震同震地表形变场

本文针对2010年4月14日玉树地震引起的地表形变,使用日本ALOS卫星PALSAR L波段雷达影像数据,应用两轨雷达差分干涉(DInSAR)处理得到了以玉树为中心11000km2范围内的同震形变场,空间分辨率为8m,并在此基础上对玉树地震的震源机制和发震机理进行了分析.研究结果表明L波段雷达数据适合在地形起伏较大的地区进行DInSAR形变探测.该同震形变场信息可为玉树地震的同震形变反演提供参考数据.该研究进一步证实DInSAR技术在大规模地表形变探测和地学研究领域具有广阔的应用前景.     

联合地震位错模型和InSAR数据构建2017年九寨沟Mw 6.5地震同震三维形变场

利用哨兵（Sentinel）-1A卫星升、降轨影像,在地震位错模型约束下获取了2017年九寨沟Mw 6.5地震的高质量三维形变场。首先,利用合成孔径雷达干涉测量技术（interferometric synthetic aperture radar,InSAR）提取九寨沟地震升、降轨同震形变场;然后,通过“两步法”反演获取该地震发震断层的几何参数和分布式滑动模型,以此为约束,采用方差分量估计算法联合解算九寨沟地震三维形变场。结果表明,九寨沟地震同震三维形变场以水平位移为主,垂向形变较弱;南北向形变呈拉张趋势,断层上盘向南、下盘向北滑动,最大位移分别为－19.81 cm和14.38 cm;东西向形变不对称性明显,断层上盘西北部向东水平运动,最大位移为18.37 cm,下盘东南部向西运动,最大位移不足8 cm。将南北、东西向形变与6个全球导航卫星系统（global navigation satellite system,GNSS）台站观测数据进行比较,两者一致性较好且均方根误差较小,分别为1.44 cm和1.77 cm,表明联合升、降轨InSAR观测和地震位错模型约束构建同震三维形变场方法具有较高可行性,显著降低了大地测量数据不足、InSAR观测对南北向形变不敏感等问题的影响。