Seismic interferometry experiment in a shallow cased borehole using a seismic vibrator source‡

We present the results of a seismic interferometry experiment in a shallow cased borehole. The experiment is an initial study for subsequent borehole seismic surveys in an instrumented well site, where we plan to test other surface/borehole seismic techniques. The purpose of this application is to improve the knowledge of the reflectivity sequence and to verify the potential of the seismic interferometry approach to retrieve high‐frequency signals in the single well geometry, overcoming the loss and attenuation effects introduced by the overburden. We used a walkaway vertical seismic profile (VSP) geometry with a seismic vibrator to generate polarized vertical and horizontal components along a surface seismic line and an array of 3C geophones cemented outside the casing. The recorded traces are processed to obtain virtual sources in the borehole and to simulate single‐well gathers with a variable source‐receiver offset in the vertical array. We compare the results obtained by processing the field data with synthetic signals calculated by numerical simulation and analyse the signal bandwidth and amplitude versus offset to evaluate near‐field effects in the virtual signals. The application provides direct and reflected signals with improved bandwidth after vibrator signal deconvolution. Clear reflections are detected in the virtual seismic sections in agreement with the geology and other surface and borehole seismic data recorded with conventional seismic exploration techniques.     

Modeling coseismic displacements resulting from SAR interferometry and GPS measurements during the 1997 Umbria-Marche seismic sequence

In this study we analyse coseismic GPS displacements and DInSAR data to constrain a dislocation model for the three largest earthquakes of the 1997 Umbria-Marche seismic sequence. The first two events, which occurred on September 26 at 00:33 GMT (Mw 5.7) and 09:40 GMT (Mw 6.0) respectively, are investigated using both GPS displacements and DInSAR interferograms. We discuss and compare the results of previous studies which separately modeled a smaller subset of geodetic data. We provide a dislocation model for these two earthquakes which fits well both GPS and DInSAR data and agrees with the results of seismological and geological investigations. The first event consists of a unilateral rupture towards the southeast with a uniform dislocation. The strike, rake and dip angles are those resulting from the CMT solution. The second event consists of an unilateral rupture towards the northwest and a variable slip distribution on the fault plane. The strike and the rake are consistent with the CMT solution, but the dip angle has been slightly modified to improve the simultaneous fit of GPS and DInSAR data. While the second rupture (09:40 GMT) arrived very close to the surface, the fit to geodetic data shows that the first rupture (00:33 GMT) is deeper (2 km), despite the more evident surface geological effects. The analysis of new SAR interferograms allows the identification of a 5–6 cm additional displacement caused by the October 3 (Mw 5.2) and 6 (Mw 5.4) seismic events.We use data from a new DInSAR interferogram to model the displacement field of the Sellano earthquake of October 14, 1997. For this event significant GPS measurements were not available. We tested two different fault plane geometries: a blind, planar fault (top depth = 2.4 km), and a curved (listric) fault reaching the surface. The two models provide a generally similar fit to the data, and show that most of the slip was released at depths greater than 2.4 km along a gently dipping (40^°–45^°) fault surface. They also show that a unilateral rupture does not allow fitting the interferometric fringes since there is evident surface deformation to the northwest of the hypocenter. Moreover, we suggest that the concentration of high residuals in the southern part of our uniform slip model may in fact indicate a certain slip variability in this area.We conclude that, despite the moderate magnitudes and the lack of significant surface faulting, the space geodetic data allowed to constrain dislocation models giving new insights in the rupture process of the three largest events of the sequence.     

南水北调总干渠中线工程豫北段基础稳定性的InSAR时序分析

南水北调总干渠中线工程豫北段受矿区开采沉陷、区域地质构造活动及城市地下水开采沉降等多种威胁,区域地质结构稳定性较差。为了评估干渠沿线基础稳定性,为引水工程提供安全评价保障,以南水北调中线豫北段为例,利用2009年全年共9期ENVISAT ASAR雷达数据,采用二通加外部DEM的D-InSAR数据处理算法和"相位累积式"干涉测量甄别大气影响,采用"相邻重访周期式"干涉测量减小时间退相干影响,获得了沿线126km渠段、33个监测点在2009年8个时间段的时序差分形变相位图像;提取了不同时期采矿、地下水开采造成的地面沉降区及沉降幅度信息;评价了区域活动构造及稳定性。研究表明:研究渠段全线有不均匀下沉,350d累积最小下沉量为-33mm,最大下沉量为-73mm。年下沉速率0.34~0.76m/a,年平均下沉速率0.53m/a。监测点下沉总体符合指数分布,平均相关系数R2=0.741 8。根据监测点下沉拟合曲线预测,调水干渠多数段落基础下沉趋稳;南水北调中线工程豫北段受多种因素叠加影响,总体基础稳定性较差,地质构造活动、城市地表沉降为主要影响因素,基础施工为次要影响,矿区开采沉陷对基础失稳未见直接相关性。     

基于气象数据的地基雷达大气扰动校正方法研究

大气扰动是影响地基雷达(Ground-based SAR,GBSAR)系统监测精度的主要因素,文中研究大气扰动对GBSAR相位的影响及相应的补偿方法。实验以自制的角反射器作为监测目标,利用气象数据补偿法改正角反射器相位,用观测区域的温度、气压和湿度建立大气折射模型,估算出大气折射率的变化来校正大气扰动误差。结果表明,文中方法能有效地剔除大气扰动相位,达到改善监测精度的效果。     

极化SAR干涉测量模拟研究

极化SAR干涉测量将极化测量技术和干涉测量技术巧妙地结合，能够提高干涉测量精度，尤其是在植被覆盖的情况下，不仅能够测量地表地形，而且使得估计植被高度成为可能。因此，近几年来国际上对极化SAR干涉测量技术的研究越来越重视，极化SAR干涉测量已成为SAR应用技术发展的重要方向。通过模拟技术深入研究极化SAR干涉测量理论和方法，并提出了一种新的极化SAR干涉测量改进算法，模拟实验结果表明，该算法可以进一步提高测量精度。     

星载SAR复数图像的配准

星载干涉雷达测量已经被证实用于生成大范围（乃至全球）的DEM或监测地表位移具有广阔的应用前景，其中SAR复数图像序列的精确配准是保证生成有效干涉相位图的重要环节之一。本文分析了卫星InSAR系统中所涉及的几种坐标系统变换问题，并提出了基于此的几何配准方案，即基于卫星轨道状态矢量、SAR成像几何进行配准的算法；同时，还设计了基于几何配准结果和能量影像相关技术的配准算法，最后使用ERS－1／2、JERS－1和RADARSAT卫星获取的3种SAR复数图像对分别做几何配准和相关配准试验和比较，并对卫星SAR图像质量和轨道数据质量作了定性分析。     

TanDEM-X双站SAR干涉测量及研究进展

较全面、系统地介绍了Tan DEM-X/Terra SAR-X双站SAR科学计划,重点涉及其科学目标、Tan DEM-X卫星参数、轨道结构以及干涉数据获取模式等相关内容,并讨论了双站SAR成像、极化In SAR和数字波束成形等干涉测量新技术及其研究进展。这些双站SAR新技术的实现将大大地推动SAR干涉测量在全球地形测绘、冰川学、海洋学及地质学等领域中的应用。     

基于线弹性位错模型反演1997年西藏玛尼Mw7.5级地震的干涉雷达同震形变场——Ⅱ滑动分布反演

1997年11月8日西藏M_w7.5级玛尼地震是干涉雷达技术应用于地震观测以来的一次重要事件.在第一部分中，我们应用广泛使用的Okada线弹性位错模型，假设断层的各个分段滑动量均匀，反演得到断层各个分段的几何参数和均匀滑动量.本部分的反演进一步去除滑动均匀假设，并利用更能反映断层真实状态的角形元位错模型（线弹性），在第一部分反演得到断层几何的基础上，反演断层面的静态位错分布.反演结果表明，线弹性滑动分布模型能够更好地解释观测数据，进一步提高反演的数据拟合程度.最终得到了断层面上的走滑和倾滑位错分布.首次得到的断层面滑动分布显示断层面滑动在浅部（0～12 km）比较集中，地震破裂长度约170 km，最大左旋走滑位移达4.8 m；反演结果还表明局部段落存在较大倾滑位移，量值达到1.9 m，这在断层模型中是不能忽略的，它可能是断层两侧形变不对称的原因之一；反演得到的标量地震矩为2.18×10²⁰ N·m，相当于矩震级M_w7.5，与Velasco等利用地震波形反演得到的结果一致.     

Point‐spread functions for interferometric imaging

Interferometric redatuming is a data‐driven method to transform seismic responses with sources at one level and receivers at a deeper level into virtual reflection data with both sources and receivers at the deeper level. Although this method has traditionally been applied by cross‐correlation, accurate redatuming through a heterogeneous overburden requires solving a multidimensional deconvolution problem. Input data can be obtained either by direct observation (for instance in a horizontal borehole), by modelling or by a novel iterative scheme that is currently being developed. The output of interferometric redatuming can be used for imaging below the redatuming level, resulting in a so‐called interferometric image. Internal multiples from above the redatuming level are eliminated during this process. In the past, we introduced point‐spread functions for interferometric redatuming by cross‐correlation. These point‐spread functions quantify distortions in the redatumed data, caused by internal multiple reflections in the overburden. In this paper, we define point‐spread functions for interferometric imaging to quantify these distortions in the image domain. These point‐spread functions are similar to conventional resolution functions for seismic migration but they contain additional information on the internal multiples in the overburden and they are partly data‐driven. We show how these point‐spread functions can be visualized to diagnose image defocusing and artefacts. Finally, we illustrate how point‐spread functions can also be defined for interferometric imaging with passive noise sources in the subsurface or with simultaneous‐source acquisition at the surface.     

InSAR相位解缠算法比较及误差分析

大跨度桥梁,其变形受风力、温度、车辆等因素的影响,其变化往往达到数十厘米或者更高。这使得利用雷达干涉技术监测其变形存在相位模糊度的问题。子带干涉技术可将距离向宽带频谱分解成两个甚至多个频带,然后进行干涉处理。该技术可以使雷达波长放大10~100倍,从而为m级的变形监测成为可能。详细介绍了子带干涉技术原理和处理方法,通过分别对风力和温度影响下的青马大桥子带干涉结果进行分析,验证子带干涉技术无需进行相位解缠处理在监测大型人工建筑物形变方面的可靠性和优势。