基于差分干涉雷达的汶川地震同震形变特点

地震同震形变场是认识震源机制、确定发震断裂、分析发震断裂活动方式、评估震害损失及揭示未来发震趋势的重要依据。但在龙门山地区复杂的自然环境中, 现有的地震同震形变场获取方法(GPS测量、构造形迹分析、震源机制解结合数值模拟计算)难以快速、全面的捕捉到这种信息, 这突显出差分干涉雷达技术(D-InSAR, Differential Interferometry SAR)在同震形变场研究中的优势。本文介绍了其观测原理, 在地震同震形变监测中的研究现状, 并重点分析了已获取的D-InSAR汶川地震同震形变场的观测结果。基于差分干涉雷达观测及相关资料, 对龙门山断裂在本次地震中活动性进行了初步分析, 结果表明干涉雷达形变图像与地表破裂范围、逆冲角度的变化、上下盘升降关系及大型余震的展布存在很强的相关性。      

巴姆地震变形场和应力场：Ⅰ．用差分干涉雷达和Okada方法求解

本文对2003年12月26日伊朗东南部克尔曼省巴姆地区发生的6.6级地震产生的位移场、应变场和应力降分布进行了定量研究。首先从差分干涉雷达(D～InSAR)出发,求解出此次地震的同震干涉条纹和雷达视线方向的位移场,并从雷达干涉的相干图上清楚的观察到了发震断层的痕迹,得出此次地震的发震断层是巴姆断层以西的隐伏断层的结论。然后结合 Harvard 的 CMT(震源机制解)结果、Nakamura 的断层结果(Arg-e-Bam 隐伏断层)以及本文 D-InSAR 求解出的发震断层结果,利用 Okada 程序,计算出此次地震的北、东、垂向三个位移分量,投影得到雷达视线方向的位移,将投影结果与由 D-InSAR得到的位移场比较,看到我们所讨论的断层模式可以很好的解释差分干涉雷达的观测结果,可以合理的解释 Arg-e-Bam 隐伏断层为此次地震的发震破裂断层,通过分析地震的破坏程度分布图,隐伏断层为发震断层的解释可以很好的符合巴姆地震的等震分布图。最后,本文利用 Okada 程序计算得到了此次地震所产生的应变场和应力降场。     

基于InSAR同震形变场的烈度评估研究

在发生破坏性地震之后,及时地获取人口聚集区的震害状况,对灾后救援工作有着十分重要的意义。利用合成孔径雷达差分干涉测量技术,可以快速、大范围地获取震区的地表形变。然而,如何基于同震形变精确地获取人口聚集区的地震烈度,进而为震后救援提供及时准确的灾情信息,还是一项有待研究的科研难题。本文以2010年4月14日青海玉树7.1级地震为例,利用人口聚集区的同震形变开展了地震烈度的评估。结果表明,利用差分干涉测量技术获取的地震同震形变,可以初步地评估人口聚集区地震烈度的分布状况。经过与实地调查烈度对比验证,评估精度达到74%左右,这表明基于地震形变的烈度评估结果与震后调查公布的结果有较高的一致性。地表形变资料对于地震烈度的早期评估、震区受灾状况的快速认识具有着重要的参考价值。     

巴姆地震变形场和应力场：Ⅱ．用FEPG有限元方法求解

本文利用 FEPG 有限元程序对2003年12月的巴姆地震变形场和应力场进行了模拟计算。2003年12月26日,巴姆地区发生6.6级地震。Nakamura、Suzuki 等人利用余震资料和地震破坏程度的分布图分析了余震震中的统计分布形式,得到了此次地震的发震断层是巴姆断层以西大约5km 的 Arg-e-Bam 隐伏断层的结论。凌勇等人从差分干涉雷达(D-InSAR)出发,得到了此次地震的同震干涉条纹和雷达视线方向的位移场,并从雷达干涉的相干图上确定了发震断层的痕迹。他们还利用 Okada 程序模拟计算了此次地震的位移场,模拟计算的结果与 D-InSAR 测量的位移场比较符合。Okada 程序从断层的位错为出发点进行计算,从运动学角度对此次地震变形场和应力场进行研究;本文从动力学的角度出发,研究远场区域边界的应力状态以及断层的物性参数的变化对巴姆地震形变场的影响,将各种计算结果进行对比的同时也和前人的结果相比较。本文计算中采用了凌勇等人的断层参数,即 Arg-e-Bam 隐伏断层为近南北走向,向东倾斜,近直立的倾角,右旋走滑,其(Strike,dip,slip)为(175°,80°,173°)。断层未延伸到地表,自地面以下1km 处往下延伸16km 的宽度,断层长度为20km。对断层的杨氏模量、泊松比、断层附近区域的边界受力大小、两个主应力的比例以及主应力的方向的改变对计算结果的影响做了讨论。最后通过上述几组计算结果建立了一个与前人结果较为符合的巴姆地震参数模型。通过计算,一方面揭示了应力状态和物性参数的影响,另一方面也说明 Arg-e-Bam 隐伏断层是此次地震的发震断层这一结论是可信的。     

差分干涉雷达方法研究及在强震形变精确测量中的应用

一、前言差分干涉雷达(SAR)属于干涉SAR的高级处理技术,应用差分干涉SAR可以监测地表的微小形变(波长级),对于研究精密大地测量、地震学、火山监测等都有重要的意义。由差分干涉处理获得的地壳形变为斜距向,为了与GPS及其它资料相比较,还需要将其分解为竖向与水平向两个分量。为此本文研究了断层活动、强震运动、卫星轨道等之间的关系,确定了不同情况下,斜距向形变分解成水平与竖向两个分量的公式。在以上工作的基础上,处理了西藏玛尼1997年Ms7.5强震区地震前后5景ERS-1/2 SAR图像,获得了差分干涉雷达图像,并精确测量了最大水平与竖向位移分量。     

差分干涉雷达遥感在地质学中的应用研究

雷达差分干涉测量是当前雷达遥感的热点研究领域，对地质学特别是地震、火山、大地构造的研究具有极其重要的意义。文章对差分干涉雷达遥感的原理及具体实现进行了深入研究，并分析了差分干涉测量结果的误差，郑重阐述了差分干涉雷达遥感在地震、地表形变测量、火山监测、板块构造研究中的应用。     

2022年青海门源MS6.9地震发震构造和地表破裂初步调查

北京时间2022年1月8日01时45分，青海省海北州门源县发生强烈地震（图1），造成数人受伤，房屋倒塌，部分道路、桥梁、隧道等基础设施被破坏或受损。中国地震台网（CENC）测定该地震的震级为MS 6. 9，震中位于37. 77°N，101. 26°E，震源深度为10 km(https://www.cenc.ac.cn/cenc/dzxx/396391/index.html)。利用欧洲航空局哨兵2号雷达卫星的震前、震后SAR数据进行差分干涉处理，得到同震形变场分布图。限定此次地震以左旋走滑运动为主，断层走向NWW，断层面近直立；主体破裂深度在10 km以上并到达地表，形成长度＞35 km的地表变形带，最大滑动量约2 m。2022年门源MS 6. 9地震发生在青藏高原中北部的祁连- 柴达木次级地块的北部（图1）、托莱山断裂带和冷龙岭断裂带的交会部位，是继1986年和2016年两次门源MS 6. 4地震之后在冷龙岭断裂带上发生的震级最高、地表破裂最长的地震事件。     

汶川、芦山地震前龙门山断裂带地壳形变特征对比分析

为研究2008年汶川地震和2013年芦山地震前地壳形变特征,本文利用1999—2015年四期GPS速度场和1990—2017年跨断层短水准资料,对跨断层GPS速度剖面、GPS应变率场、断层闭锁程度和滑动亏损、跨断层年均垂直变化速率等进行了分析讨论,总结了汶川和芦山地震前后龙门山断裂带三维地壳变形演化特征。结果表明,汶川地震前龙门山断裂带中、北段处于强闭锁状态、断层面应力应变积累水平很高,而龙门山断裂带西南段闭锁较弱、变形速率明显高于中北段、依然可以积累应力应变,汶川地震震源位于闭锁相对弱的部位,这可能是导致汶川地震自初始破裂点沿龙门山断裂带向北东方向单侧破裂,而震中西南方向断层并没有发生破裂的原因之一。汶川地震的发生引起龙门山断裂带西南段应力应变积累速率加快、断层闭锁程度增强、闭锁面积增大,这在一定程度上促进了芦山地震的发生,而芦山地震震源位于汶川地震前强闭锁和弱闭锁的高梯度过渡部位。因为芦山地震只释放了龙门山断裂带西南段有限的应变能,并没有显著缓解该段的地震危险性,所以汶川和芦山地震之间的地震空段以及芦山地震西南方向的地震空段,依然需要持续关注。此外,本文还收集和对比分析了多次6～9级地震前地壳变形特征,同样显示地震成核于闭锁高梯度带区域而非完全闭锁区域内部,并且随着震级升高闭锁断层面的长度也在增大,这一现象还需在高分辨率形变数据的帮助下进行深入研究和分析。     

汶川5.12地震发震机理的数值模拟

将围岩看成弹性体, 断层看成具有应变软化的弹塑性体, 断层和围岩组成统一的地质介质系统。在给定的地应力、孔隙压力及边界位移的作用下, 应力逐渐积累, 当达到断层强度时, 断层进入塑性状态, 应变随之软化, 应力突然下降, 能量突然释放, 形成地震。根据这个原理, 利用有限差分FLAC软件, 对汶川5.12地震的机理进行了模拟。模拟结果证实:汶川5.12地震是在青藏高原侧向挤压, 四川盆地阻挡, 使龙门山断裂带受到缓慢增加的挤压应力场作用下形成的。计算结果给出了应力降、能量释放量、断层错动量、地震同震位移、震前位移、地震复发周期等重要参数, 结果与野外调查及其它资料给出的数据具有较好的一致性。      

差分干涉雷达测量在地面沉降中的应用研究

雷达差分干涉测量是当前雷达遥感的热点研究领域，利用遥感卫星多时相的复雷达图像相干信息进行地表的垂直形变量的提取，其精度已达到了毫米级。对差分干涉雷达遥感的原理及具体实现进行了深入研究，并分析了差分干涉测量结果的误差，着重阐述了差分干涉雷达遥感在地面沉降研究中的应用，及其需要进一步解决的问题。     

Surface displacement estimation along Himalayan frontal fault using differential SAR interferometry

The Himalayan region has been studied extensively during the past few decades in terms of present ongoing deformations. Various models have been proposed for the evolution of the Himalaya to explain the cause of earthquake occurrences and to understand the seismotectonics of the Himalayan collision zone. However, the information on displacements from field geodetic surveys is still too scarce in time and spatial domains so as to provide convincing evidences. Moreover, classical Probabilistic Seismic Hazard Approaches also fail due to paucity of data in higher magnitude range, thus emphasizing the need of spatial level displacement measurements. It is in this context that the present study has been carried out to estimate the surface displacement in a seismically active region of the Himalaya between Ganga and Yamuna Tear using Differential SAR interferometry. Three single-look complex images, obtained from ASAR sensor onboard ENVISAT satellite, have been used. A displacement rate of 8?C10?mm per year in N15°E direction of Indian plate has been obtained in this three-pass SAR interferometry study. It has been noted that the estimated convergence rate using Differential SAR interferometry technique is relatively low in comparison with those obtained from previous classical studies. The reported low convergence rate may be due to occurrence of silent/quite earthquakes, aseismic slip, differential movement of Delhi Hardwar ridge, etc. Therefore, in view of the contemporary seismicity and conspicuous displacements, a study of long-term observations of this surface movement has been recommended in future through a time-series SAR interferometry analysis.     

差分干涉SAR冻土形变检测方法研究

季节性冻融导致的地表形变是冻土地区工程建设的主要病害,冻土的冻胀和融沉是影响青藏公路以及目前建设的青藏铁路路基稳定的重要因素.近年来,差分干涉测量技术已成为地表形变测量和监测的重要工具.基于重复轨道的ERS 1/2雷达图像,研究了利用干涉SAR技术探测冻土形变的方法,经过对引入DEM差分干涉SAR技术的处理和分析,得到研究区冻土形变结果.通过与实测数据比较,表明重轨差分干涉测量可以精确地探测冻土表面形变,可用于青藏铁路冻土形变监测.     

Radar interferometry techniques for the study of ground subsidence phenomena: a review of practical issues through cases in Spain

Subsidence related to multiple natural and human-induced processes affects an increasing number of areas worldwide. Although this phenomenon may involve surface deformation with 3D displacement components, negative vertical movement, either progressive or episodic, tends to dominate. Over the last decades, differential SAR interferometry (DInSAR) has become a very useful remote sensing tool for accurately measuring the spatial and temporal evolution of surface displacements over broad areas. This work discusses the main advantages and limitations of addressing active subsidence phenomena by means of DInSAR techniques from an end-user point of view. Special attention is paid to the spatial and temporal resolution, the precision of the measurements, and the usefulness of the data. The presented analysis is focused on DInSAR results exploitation of various ground subsidence phenomena (groundwater withdrawal, soil compaction, mining subsidence, evaporite dissolution subsidence, and volcanic deformation) with different displacement patterns in a selection of subsidence areas in Spain. Finally, a cost comparative study is performed for the different techniques applied.     

Application of holographic interferometry to the study of time-dependent behavior of rock and coal

Holographic interferometry has been applied in various fields of science and engineering. One of the most important functions of the technique is the capability of measuring displacement of any surface with a high accuracy without emplacing any gauges. In this study, the real-time technique of holographic interferometry was utilized as a new creep-testing method of rock and coal. Shale, sandstone and coal specimens were collected from the Alabama coal mines, and they were tested on a specially designed three-point loading apparatus. The results showed that the relationship of creep deflection vs. time generally followed the classical form of creep, and it varied according to physical conditions of the rock and coal specimens. This study proved that the new method can be successfully applied to the study of time-dependent behaviour of rock and coal with a high resolution and a low cost.     

Analysis with C- and X-band satellite SAR data of the Portalet landslide area

This paper describes the use of the Stable Point Network technique, a Persistent Scatterer Interferometry SAR technique, for the analysis of the Portalet landslide area (Central Pyrenees, Spain). For this purpose, different SAR datasets acquired by ERS-1, ERS-2, ENVISAT and TerraSAR-X satellites have been analysed. The use of different SAR images acquired by satellite radar sensors operating at different microwave lengths has allowed for a comparative assessment and illustration of the advantages and disadvantages of these satellites for landslide detection and monitoring. In the introduction, differential interferometry and the study area are briefly described. Then the specifics of the SPN processing and the results of the different datasets are described and compared. In Analysis of the results: the Portalet landslide area, the Portalet landslide area is introduced and the radar displacement measurements are analysed with available geo-information data. Additionally, X-band measurements are compared with those gathered by a ground-based SAR for a previous project. Finally, the most relevant conclusions of this work are discussed.     

Accuracy assessment of InSAR derived input maps for landslide susceptibility analysis: a case study from the Swiss Alps

In recent years SAR interferometry has become a widely used technique for measuring altitude and displacement of the surface of the earth. Both these capabilities are highly relevant for landslide susceptibility studies. Although there are many problems that make the use of SAR interferometry less suitable for landslide inventory mapping, it’s use in landslide monitoring and in the generation of input maps for landslide susceptibility assessment looks very promising. The present work attempts to evaluate the usefulness and limitations of this technique based on a case study in the Swiss Alps. Input maps were generated from ERS repeat pass data using SAR interferometry. A land cover map has been generated by image classification of multi-temporal SAR intensity images. An InSAR DEM was generated and a number of maps were derived from it, such as slope-, aspect, altitude- and slope form classes. These maps were used to generate landslide and rockfall susceptibility maps, which give fairly well acceptable results. However, a comparison of the InSAR DEM with the conventional Swisstopo DEM, indicated significant errors in the absolute height and slope angles derived from InSAR, especially along the ridges and in the valleys. These errors are caused by low coherence mostly due to layover and shadow effects. Visual comparison of stereo images created from hillshading maps and corresponding DEMs demonstrate that a considerable amount of topographic details have been lost in the InSAR-derived DEM. It is concluded that InSAR derived input maps are not ideal for landslide susceptibility assessment, but could be used if more accurate data is lacking.     

Satellite remote sensing-based detection of the deformation of a reservoir bank slope in Laxiwa Hydropower Station, China

Laxiwa Hydropower Station is the largest among those hydropower stations in the upper reaches of the Yellow River. The construction started in October 2001, diversion tunnel was finished in January 2004, and impoundment started in March 2009. However, from May 2009, the right-bank slope of the reservoir about 700 m high and 1,000 m wide located 500 m far from the dam was found to be deforming greatly and continuously. Although this slope had been identified as an old landslide, the survey before the construction of the dam concluded that this slope being composed of granite is stable and would be stable even after the impoundment, and thus no detailed monitoring of the slope deformation had been performed before the visible deformation occurred after the impoundment. To indentify the relationship between the deformation of the slope and impoundment, we utilized differential synthetic aperture radar interferometry and Advanced Land Observing Satellite Prism data to analyze the slope deformation before and after impoundment. We found that no identifiable deformation took place before the impoundment, while during the period of 3 April 2009 to 22 May 2010, after impoundment, maximum horizontal displacement reached approximately 7.5 m. This slope is still deforming even after a total horizontal displacement of several tens of meters being reached, showing high risk of catastrophic failure.     

三峡库区山区城镇重大地质灾害监测预警示范研究

以三峡库区羊角场镇为例,在地质调查基础上,利用伸缩位移计、应力计等常规地面监测设备,对地表裂缝相对位移和危岩体基座应力进行监测,并结合GPS静态和实时动态测量等技术,开展了角反射器干涉测量技术（CR-InSAR）和高相干InSAR技术在复杂山区地质灾害大范围调查与识别应用,从点到面、从大型危岩体到区域地表形变全覆盖。各种方法监测结果较为一致,有效实现了典型山区城镇地质灾害变形的高精度、动态立体化监测。      

On the uniformity of anthropogenic Venice uplift

A programme of Venice uplift because of seawater injection into a 600–800 m deep brackish aquifer underlying the lagoon has recently been advanced. As the actual spatial variability of the hydro‐geological parameters is to a large extent unknown, a controversial issue concerns the 25–30 cm heave prediction with the possible generation of differential displacements that might jeopardize the integrity of the architectural patrimony. It is shown that significant differential displacements have occurred in Venice in the past and are still presently occurring as evidenced by SAR interferometry. The results of a stochastic analysis addressing the variability of the hydraulic conductivity of the injected formation indicate that even a highly uneven aquifer expansion does not migrate to ground surface because of the smoothing effect exerted by the overburden. The predicted differential displacements are well below the values Venice is experiencing. Any a priori alarmism appears to be unfounded and unjustified.