基于PSInSAR技术的海原断裂带地壳形变初步研究

常规差分干涉测量(DInSAR)受时间、空间失相干的严重制约和和大气延迟等相位误差的影响,难以实现对长期累积微小地壳形变场的有效探测.PSInSAR技术克服了常规DInSAR的局限性,能够高精度监测微小地壳形变.本文首先介绍了PSInSAR技术的算法模型和处理方法.该方法通过二维线性相位模型,对时序干涉图象上相干点目标...     

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

基于多期DInSAR数据的2016年定结地震序列发震断层特征研究

2016年5月22日,在西藏定结县发生四次MW4~5地震,研究本序列地震的发震断层几何和运动特征对于认识周边活动断裂性质具有重要意义.由于发震地区偏远,且观测台网分布稀疏,本研究采用星载雷达干涉测量(DInSAR)技术进行了同震形变场重建,但是定结地震震级较小,单干涉像对获取的形变场受相位噪声影响较大.为了解决这一问题,本研究基于时间序列Sentinel-1A干涉数据生成多期同震与非同震干涉图,并利用叠加平均法对本次定结地震同震形变场进行重建,提取了定结2016年5月22日多次地震产生的同震累计整体形变场.基于InSAR同震形变场和区域地质特征,研究进行了滑动分布反演,确定其主要贡献的发震断层几何参数及滑动分布:断层走向为188°,倾角为43°,平均滑动角为78°,发震断层的运动性质以正断为主兼具少量左旋走滑分量,滑动主要集中在断层垂直深度0~9km处,最大滑动量约为25cm,位于断层倾向深度3km处,反演得到的矩震级为MW5.58.本研究结果表明采用星载InSAR叠加平均技术可以较好地压制相位噪声,有效提取此类中小型浅源地震同震微弱形变场.最后,我们认为本次定结地震与藏南拆离断层与申扎-定结断层的活动密切相关.     

综合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走向的右旋走滑为主兼有部分正断成分.应变张量分析表明发震断层处受到了明显的收缩力和剪切力的作用.     

雷达干涉测量在地震形变研究中的应用

合成孔径雷达干涉测量(DInSAR)以厘米乃至毫米的精度和很大的空间覆盖范围实现地表形变观测,它是进行地壳形变观测的重要手段,在地震形变研究中发挥着越来越重要的作用。本文首先简要给出了微波雷达测量的特征,阐述了雷达图像数据处理的基本方法,分析了InSAR技术的主要影响因素和观测精度。结合较新的研究成果,分析了雷达遥感图像研究地震形变的方法,讨论了DInSAR在震间、同震和震后形变观测研究中的最新进展。同时讨论了目前我国应用DInSAR研究地震形变中面临的主要问题,最后阐述了应用雷达干涉测量研究地震形变的发展趋势。     

利用合成孔径雷达差分干涉测量技术监测伊朗Bam地震同震形变场

介绍了差分干涉测量的原理、差分干涉数据对的选取方法,以及三轨法差分干涉测量数据处理的流程.利用星载合成孔径雷达差分干涉(D-InSAR)测量技术和ENVISAT ASAR雷达数据对2003年12月26日伊朗Bam 6.5级地震引起的地表形变场进行了测量试验,成功地获取了Bam地震的蝴蝶状的同震形变场生成了地表形变的等值线图,并且根据相干图确定了地震造成破坏最严重的区域的位置、分布及面积.试验证明D-InSAR技术是地表形变测量和地震研究的一个强大和有效的工具.     

王艳，廖明生，李德仁等.利用长时间序列相干目标获取地面沉降场.地球物理学报，2007

方法，相干目标分析是利用时间序列多景影像获取地表稳定反射体——高相干点目标的新方法.实验区选择上海市主城区约100 km2 区域，InSAR的实验数据采用覆盖1992～2000年间的25景ERS-1/2的单视复影像.在无先验沉降场模型的情况下，实验结果显示该实验区的地面沉降场在时间范围上跨越8年，平均观测时间采样率约为4个月，垂直形变精度优于5 mm，满足我国对城市地面沉降形变观测的要求.研究结果表明该方法可以用于今后大时间/空间尺度上的高精度地表形变观测.     

基于InSAR、GPS形变场的长白山地区火山岩浆囊参数模拟研究

利用差分合成孔径雷达干涉测量技术（DInSAR）获取长白山地区的形变场.结果显示1995~1998年期间,位于长白山东南侧的间白山火山存在6~12 cm的视线向形变,而长白山天池火山处于平静期,没有明显形变.利用2002~2003年的GPS和水准获取的形变数据,分别采用Mogi单源、双源模型反演了长白山地区火山的岩浆囊参数.其中双源模型拟合效果较为理想,两个点源一个位于长白山天池老火山口下方7.9 km处,另一个位于间白山火山下方5.5 km处.对双源模型反演得到的岩浆囊参数进行适当调整,拟合得到与InSAR形变场基本吻合的结果.上述研究结果表明长白山地区火山活动存在时间上的间歇性和空间上的迁移性,为进一步研究长白山地区火山活动机制提供了参考和依据.     

王艳，廖明生，李德仁等.利用长时间序列相干目标获取地面沉降场.地球物理学报，2007

方法，相干目标分析是利用时间序列多景影像获取地表稳定反射体——高相干点目标的新方法.实验区选择上海市主城区约100 km2 区域，InSAR的实验数据采用覆盖1992～2000年间的25景ERS-1/2的单视复影像.在无先验沉降场模型的情况下，实验结果显示该实验区的地面沉降场在时间范围上跨越8年，平均观测时间采样率约为4个月，垂直形变精度优于5 mm，满足我国对城市地面沉降形变观测的要求.研究结果表明该方法可以用于今后大时间/空间尺度上的高精度地表形变观测.     

基于分区解缠方式研究2010年玉树M_W6.9地震同震形变

2010年玉树MW6.9地震发生后,已有学者利用DInSAR技术得到了该次地震的同震形变场,并在此基础上反演了其震源破裂滑动分布。本文以提高玉树地震同震形变场准确度作为出发点,设计了一种新的解缠方式获取相应的同震形变场:首先沿地表同震破裂迹线将InSAR干涉图像分割成上下两部分,然后利用网络流解缠算法对其分别进行相位解缠,最后通过统计重叠部分的相位一致性信息将两部分结果进行拼接。利用震中区GPS同震位移数据对分区解缠结果进行验证。结果表明:相对于整体解缠结果,分区解缠方式得到的最大视线向沉降值由42.6cm增加到48.1cm,所得结果与GPS观测结果更加接近,采用分区解缠方式提高了断层附近形变场的准确性。     

永久散射体雷达差分干涉应用于区域地表沉降探测

永久散射体干涉测量技术可以克服常规干涉方法在区域地表形变测量中的雷达信号失相关和大气延迟影响. 本文对基于永久散射体的干涉处理全过程进行了分析, 基于差分相位的两个主要特性: 沉降信号时序相关、地形相位与空间基线成比例, 提出构建PS网络, 并采用间接观测平差法估计沉降速度网和高程修正网的全局最优解. 实验选取上海地区近10年间的ERS-1/2卫星C波段SAR数据进行干涉处理, 在永久散射体上分离形变信号、高程修正和大气分量, 并最终提取上海地区高分辨率地面沉降速度场, PS结果与地面精密观测成果吻合较好.     

DInSAR measurements of ground deformation by sinkholes,mining subsidence,and landslides,Ebro River,Spain

Differential Interferometric Synthetic Aperture Radar (DInSAR) has been applied to detect and measure ground deformation in a stretch of the Ebro River valley (Spain) excavated in salt‐bearing evaporites. The capability of the Small Baseline Subset (SBAS) DInSAR technique to detect ground displacement is analyzed comparing the DInSAR results with the available geomorphological information. The deformation map derived from 27 European Remote Sensing (ERS) satellite images covering more than five years provides sub‐centimeter displacement measurements in zones coinciding with known active sinkholes and landslides. Moreover the map provides the first account of mining subsidence in the area. The measured deformation rates reach 1·68 cm/y for the sinkholes, 0·80 cm/y for the landslides and 1·45 cm/y for the area affected by mining subsidence. The SBAS DInSAR technique provided deformation measurements in a small proportion (5–10%) of the known active sinkholes and landslides. This limitation is mainly due to the lack of coherence in agricultural areas, the spatial resolution of the deformation map (pixel size of 90 m), and the parallelism between the ERS satellite line of sight and the linear escarpment on which most of the landslides occur. Despite this, the interferometric technique provides valuable data that complement traditional geomorphological studies including the quantification of the deformation phenomena, the identification of mining subsidence otherwise only recognizable by geodetic methods, and the detection of creep deformation which might correspond to premonitory indicators of catastrophic sinkholes and landslides capable of causing the loss of lives. Detailed DInSAR studies combined with field data would be required to improve the analysis of each deformation area. Copyright © 2009 John Wiley & Sons, Ltd.     

Land subsidence along the Ionian coast of SE Sicily (Italy), detection and analysis via Small Baseline Subset (SBAS) multitemporal differential SAR interferometry

The paper presents the results of a multi‐temporal, differential interferometric synthetic aperture radar (DInSAR) analysis aiming to identify active surface deformation phenomena in south‐eastern Sicily. The study area has been chosen because of its strong seismicity, high concentration of industrial and agricultural activities, and high density of people living in the coastal area. Furthermore, the morphology, lithology and climatic features of this sector of the Hyblean foreland are suitable for an interferometric analysis, providing a high coherence over the area. The Small BAseline Subset (SBAS) multitemporal DInSAR technique was used to calculate mean ground velocity maps and displacement time series from a large data set of European Remote Sensing Satellites (ERS 1–2) images spanning the time period 1992–2000. The reliability of the DInSAR results was tested calculating the East_SAR and Up_SAR values over two permanent global positioning system (GPS) stations in the area, and comparing them with the East_GPS and Up_GPS values. The residuals between GPS and DInSAR velocities were 1 and 0.6 mm/yr for the Up and East components, respectively. Four main subsiding areas, previously undetected, have been identified, corresponding with the towns of Augusta, Siracusa, Priolo, and Villasmundo. The observed deformation phenomena are located within coastal structural basins, filled with Pleistocene and Holocence deposits, except the Villasmundo land subsidence, which is located on the Hyblean plateau. The measured deformation rates reach values up to ?18 mm/yr in Augusta, –6 mm/yr in Siracusa, –5 mm/yr in Villasmundo and ?4.5 mm/yr in Priolo. The examination of velocity profiles, time series, and geological data allows us to relate all the detected deformation patterns primarily to groundwater over‐exploitation. A multi‐dimensional interpolation with kriging was performed to obtain a field subsidence map. A first‐order elastic deformation model was used to simulate the peculiar features of the Villasmundo subsidence. Copyright © 2011 John Wiley & Sons, Ltd.     

雷达干涉PS网络的基线识别与解算方法

时序雷达干涉图中的永久散射体(PS)可看作“天然GPS点”, 以构成网络用于监测长期的地表形变. 本文提出采用邻接矩阵拓扑模型对基于Delaunay剖分算法生成的PS网络进行基线识别, 并采用时序相干最大化算法求解PS基线的线性形变速度增量和高程误差增量. 该数据模型和计算方法被应用于探测香港地区2006~2007年间的区域地表沉降. 实验研究采用由Envisat卫星ASAR传感器对该地区成像所获取的时序SAR影像作为数据源, 并联合该地区12个GPS连续运行参考站的观测数据予以大气修正和地面控制. 实验结果表明, 该模型和方法应用于地表形变测量是有效的和可靠的, PS网络方法探测地面沉降的精度约为±2.0 mm/a.     

Identification of potential subsidence related to pumping in the Almería basin (SE Spain)

Ground subsidence of detrital deposits in the Almería basin (SE Spain) was studied using the remote sensing technique of Differential Interferometry SAR (DInSAR). This basin is one of the most arid in Europe, receiving an average rainfall of 250 mm per year. Over the last 60 years the region has experienced an enormous agricultural and urban expansion, whose water demand has been largely supplied from groundwater, leading to the current situation of overexploitation of water resources. This paper outlines the likely relationship between groundwater abstraction and subsidence. To this end, 34 ERS and Envisat images, taken between 2003 and 2009, were analysed to estimate ground surface deformations, and hence, compared with water table variations measured in a number of piezometers in the basin. The analysis shows a clear parallelism between the variations in piezometric level and deformation of the ground surface. In addition, the zones of greatest subsidence coincide with those areas where groundwater abstractions are concentrated. Subsidence over the examined period varies from 10 to 30 mm, with extreme values as high as 50 mm, which translates to a rate of between 1·7 and 5 mm/year, reaching maximum rates of 8 mm/year at some points. Given such subsidence rates, damage to urban infrastructures are, for the moment, incipient. Copyright © 2011 John Wiley & Sons, Ltd.     

多轨PSInSAR监测华北平原地表垂直形变场

针对PSInSAR技术监测大空间尺度地表垂直形变场的不足,提出通过基准转换与数据拼接、融合不同轨道多个独立片区(Frame)获取大空间尺度地表垂直形变场方法.由于不同轨道主影像选取的参考基准不一致、不同轨道垂直基线不同、不同片区相位解缠的参考PS不同等因素影响,使得PSInSAR获得的不同轨道地表垂直形变场的空间基准存在差异,提取的地表形变信息不具有可比性.利用两个相邻片区重叠区域的PS点集,根据迭代最近点算法提取两片区公共区域内同名点,并兼顾同名点距离的权,精确计算两片区的配准参数,实现两片区的"无缝"连接,最后使各独立片区拼接成一个连续的大空间尺度地表垂直形变场.本文推导了多轨PSInSAR技术获得的不同片区形变场基准转换与数据拼接的数学模型,并以华北平原(115.32°E—118.79°E, 36.81°N—40.58°N)为实验区,解析了3个不同轨道共12个片区PS目标,获取了研究区2006—2010年地表垂直形变速度场.分析表明:(1)研究区大范围处于下沉状态,形成以北京、廊坊、天津、沧州、泊头—德州等城市为中心向外扩展的地表沉降发展态势,几个沉降中心的平均沉降速率分别达到-34.7 mm/a、-26.3 mm/a、-64.2 mm/a、-34.6 mm/a和-37.7 mm/a;(2)地表沉降的空间分布特征表明,城市工业生产和生活抽取地下水是地表沉降的主要诱因,农业灌溉和油气开采是导致华北平原大范围地表沉降的另一个重要因素;(3)研究区最大沉降带沿北北东向展布,与区内断裂分布一致,表明沉降的空间分布受到断裂带控制;(4)利用地表同期一等水准测量成果检验了本研究结果,精度达到4.72 mm,表明本文提出的数据处理策略能够满足大时空尺度地表形变监测的需求.     

Surface deformation analysis in the Ischia Island (Italy) based on spaceborne radar interferometry

The island of Ischia, located at the northwestern end of the Gulf of Napoli (Italy), is a volcanic area that is historically active (the Arso eruption, in 1302 and the Casamicciola earthquake, in 1883) and has diffuse hydrothermal phenomena. We present in this work a study of the surface deformation occurring in the island, which is based on applying the Differential Synthetic Aperture Radar Interferometry (DInSAR) algorithm referred to as Small BAseline Subset (SBAS) technique. This study is focused on the 1992–2003 time interval and SAR data acquired by the European Remote Sensing (ERS) satellites from ascending and descending orbits have been used, thus allowing us to discriminate the vertical and east–west components of the displacements. A validation of the DInSAR results has been carried out first by comparing the vertical deformations estimated from the SAR data with those measured from the spirit leveling network that is present in the area. In particular, we computed the difference between the mean vertical deformation velocities estimated from the SAR and the corresponding geodetic measurements along three main leveling lines; the maximum value of the root mean square difference is of about 1 mm/yr. The final discussion is dedicated to the interpretation of the detected displacements, benefiting from the overall information extracted from the ascending and descending DInSAR measurements. In particular, DInSAR data relative to the vertical deformation component show that the present-day subsidence of Ischia mainly develops in areas characterized by active landsliding and along faults; moreover, the deflation of the island, which is recorded by the horizontal displacement component, is probably related to the de-pressurization of the hydrothermal system.     

Lithological control of land subsidence induced by groundwater withdrawal in new urban AREAS (Granada Basin,SE Spain). Multiband DInSAR monitoring

Ground subsidence in the southeastern border of the Granada Basin (SE Spain) has been studied using remote sensing techniques. Over the last decades, the region has experienced a huge urban expansion, which has caused a substantial increase in water supply requirements. Water needs are exclusively met by groundwater by means of numerous pumping wells, which exploit a confined detrital aquifer of alluvial fan deposits with a heterogeneous facies distribution. A general piezometric level decline (up to 50 m) has been recorded in the aquifer during the past 30 years that has induced the generation of a subsiding area with oval shape oriented WNW‐ESE just where the new urban areas and pumping wells are located. Subsidence has been monitored by exploiting synthetic aperture radar (SAR) images from ENVISAT (2003–2009) and Cosmo‐SkyMed (2011–2014). A new approach, which combines A‐DInSAR and small‐area persistent scatterer interferometry (PSI) analysis, has been applied obtaining a good accuracy regarding temporal and spatial dimension of the subsidence. ENVISAT data (2003–2009) reveal subsidence rates up to 10–15 mm/year, and Cosmo‐SkyMed (2011–2014) values slightly lower; up to 10 mm/year. Temporal variations in the subsidence velocity are in accordance with the rainfall pattern and piezometric fluctuations in the aquifer. The sector with highest rates of subsidence does not correspond to the area with more intense groundwater exploitation but to the area with greater presence of clays in the confining layer of the aquifer. There is a clear lithological control in the spatial distribution of the ground subsidence. This work integrates detailed geological and hydrogeological data with differential SAR interferometry monitoring with the aim to better understand subsidence processes in detrital aquifers with small‐scale heterogeneity. Copyright © 2016 John Wiley & Sons, Ltd.