External atmospheric corrections in geodetic very-long-baseline interferometry

Three methods to correct for the atmospheric propagation delay in very-long-baseline interferometry (VLBI) measurements were investigated. In the analysis, the NASA R&D experiments from January 1993 to June 1995 were used. The methods were compared in correcting for the excess propagation delay due to water vapour, the “wet” delay, at one of the sites, the Onsala Space Observatory on the west coast of Sweden. The three methods were: (1) estimating the wet delay using the VLBI data themselves; (2) inferring the wet delay from water vapour radiometer (WVR) data, and (3) using independent estimates based on data from the global positioning system (GPS). Optimum elevation cutoff angles were 22^∘ and 26^∘ when using WVR and GPS data, respectively. The results were found to be similar in terms of reproducibility of the estimated baseline lengths. The shortest baselines tend to benefit from external measurements, whereas the lack of improvement in the longer baselines may be partly due to the large amount of data thrown away when removing observations at low elevation angles. Over a 2 week period of intensive measurements, the two methods using external data showed an overall improvement, for all baseline lengths, compared to the first method. This indicates that there are long-term systematic errors in the wet delay data estimated using WVR and GPS data. Received: 27 October 1998 / Accepted: 20 May 1999     

Renewed uplift at the Yellowstone Caldera measured by leveling surveys and satellite radar interferometry

 A first-order leveling survey across the northeast part of the Yellowstone caldera in September 1998 showed that the central caldera floor near Le Hardy Rapids rose 24±5 mm relative to the caldera rim at Lake Butte since the previous survey in September 1995. Annual surveys along the same traverse from 1985 to 1995 tracked progressive subsidence near Le Hardy Rapids at an average rate of –19±1 mm/year. Earlier, less frequent surveys measured net uplift in the same area during 1923–1976 (14±1 mm/year) and 1976–1984 (22±1 mm/year). The resumption of uplift following a decade of subsidence was first detected by satellite synthetic aperture radar interferometry, which revealed approximately 15 mm of uplift in the vicinity of Le Hardy Rapids from July 1995 to June 1997. Radar interferograms show that the center of subsidence shifted from the Sour Creek resurgent dome in the northeast part of the caldera during August 1992 to June 1993 to the Mallard Lake resurgent dome in the southwest part during June 1993 to August 1995. Uplift began at the Sour Creek dome during August 1995 to September 1996 and spread to the Mallard Lake dome by June 1997. The rapidity of these changes and the spatial pattern of surface deformation suggest that ground movements are caused at least in part by accumulation and migration of fluids in two sill-like bodies at 5–10 km depth, near the interface between Yellowstone's magmatic and deep hydrothermal systems. Received: 30 November 1998 / Accepted: 16 April 1999     

Feasibility of borehole ambient noise interferometry for permanent reservoir monitoring

The analysis of seismic ambient noise acquired during temporary or permanent microseismic monitoring campaigns (e.g., improved/enhanced oil recovery monitoring, surveillance of induced seismicity) is potentially well suited for time‐lapse studies based on seismic interferometry. No additional data acquisition required, ambient noise processing can be automatized to a high degree, and seismic interferometry is very sensitive to small medium changes. Thus there is an opportunity for detection and monitoring of velocity variations in a reservoir at negligible additional cost and effort. Data and results are presented from an ambient noise interferometry study applied to two wells in a producing oil field in Romania. Borehole microseismic monitoring on three component geophones was performed for four weeks, concurrent with a water‐flooding phase for improved oil recovery from a reservoir in ca. 1 km depth. Both low‐frequency (2 Hz–50 Hz) P‐ and S‐waves propagating through the vertical borehole arrays were reconstructed from ambient noise by the virtual source method. The obtained interferograms clearly indicate an origin of the ambient seismic energy from above the arrays, thus suggesting surface activities as sources. It is shown that ambient noise from time periods as short as 30 seconds is sufficient to obtain robust interferograms. Sonic log data confirm that the vertical and horizontal components comprise first arrivals of P‐wave and S‐waves, respectively. The consistency and high quality of the interferograms throughout the entire observation period further indicate that the high‐frequency part (up to 100 Hz) represents the scattered wave field. The temporal variation of apparent velocities based on first‐arrival times partly correlates with the water injection rate and occurrence of microseismic events. It is concluded that borehole ambient noise interferometry in production settings is a potentially useful method for permanent reservoir monitoring due to its high sensitivity and robustness.     

基于相邻虚拟道叠加的超虚折射干涉法及其在广角OBS折射波增强中的应用

广角地震的远偏移距折射初至含有地下深层的信息,提高远偏移距折射波信噪比能够有效提高初至拾取精度,对于深部结构的层析速度建模十分有利.超虚折射干涉法基于干涉原理对远道折射波进行增强,它在炮点和检波点都非常密集的情况下效果较好.然而,在应用于台站间距较大且环境噪声较强的广角海底地震仪(OBS)观测时,该方法对折射波的增强能力不足,而且容易产生虚假波形,造成增强后的折射波信噪比仍然较低.针对这种情况,本文提出基于相邻虚拟道叠加的超虚折射干涉法,通过叠加相邻虚拟道来提高远道与近道互相关的准确度,以达到稳定增强远道折射波信噪比的目的.理论实验和实际资料测试均显示,基于相邻虚拟道叠加的超虚折射干涉法在台站间距较大和信噪比较低的情况下能够准确构建虚拟道,且增强后的波形同相轴连续程度和信噪比均较高,有利于拾取高精度的折射初至到时.本文方法也可用于增强陆上炮间距较大的广角地震数据折射波.     

海潮负荷对沿海地区宽幅InSAR形变监测的影响

海岸带地区是全球自然生态环境最为复杂和脆弱的地域之一,合成孔径雷达干涉测量（InSAR）技术可以为全球人类活动、气候变暖和俯冲带剧烈构造运动等背景下的大范围海岸带地理环境变化研究提供重要观测资料.海洋潮汐导致固体地球长周期形变,波长尺度为10²~10³ km的海潮负荷引入mm级至cm级的形变梯度,此类非构造信号对海岸带InSAR精密形变分析（如：大范围、微小、缓慢且非稳态构造过程等）造成显著影响.本文以宽幅模式SAR数据为例,基于多种海潮模型研究了全球典型海岸带地区（福建、智利和阿拉斯加湾）海潮负荷效应对宽幅InSAR形变监测的影响,给出了宽幅InSAR海潮负荷三维分量估计与差分相位提取方法,并进一步讨论了基于不同海潮模型估计海潮负荷位移的差异.海潮负荷影响不仅与研究范围大小有关,其形变梯度变化与研究区域地形特征存在强相关,对于长波长形变分析而言,传统平面或者曲面拟合方法难以有效分离海潮负荷位移.     

同波束干涉测量在月面相对定位中的应用

基于同波束干涉测量原理和观测模型,针对月面目标特性,推导了给定着陆器先验位置的差分时延、差分时延率观测方程。引入指定高程约束、月球数字地形模型等约束方程,提出巡视器相对位置确定的卡尔曼滤波算法。通过仿真数据进行了验证和评估,表明该算法可快速、较高精度地确定月面巡视器相对着陆器的位置。     

Scattered ground‐roll attenuation using model‐driven interferometry

Scattered ground roll is a type of noise observed in land seismic data that can be particularly difficult to suppress. Typically, this type of noise cannot be removed using conventional velocity‐based filters. In this paper, we discuss a model‐driven form of seismic interferometry that allows suppression of scattered ground‐roll noise in land seismic data. The conventional cross‐correlate and stack interferometry approach results in scattered noise estimates between two receiver locations (i.e. as if one of the receivers had been replaced by a source). For noise suppression, this requires that each source we wish to attenuate the noise from is co‐located with a receiver. The model‐driven form differs, as the use of a simple model in place of one of the inputs for interferometry allows the scattered noise estimate to be made between a source and a receiver. This allows the method to be more flexible, as co‐location of sources and receivers is not required, and the method can be applied to data sets with a variety of different acquisition geometries. A simple plane‐wave model is used, allowing the method to remain relatively data driven, with weighting factors for the plane waves determined using a least‐squares solution. Using a number of both synthetic and real two‐dimensional (2D) and three‐dimensional (3D) land seismic data sets, we show that this model‐driven approach provides effective results, allowing suppression of scattered ground‐roll noise without having an adverse effect on the underlying signal.     

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.     

1D System identification of a 54‐story steel frame building by seismic interferometry

A 54‐story steel, perimeter‐frame building in downtown Los Angeles, California, is identified by a wave method using records of the Northridge earthquake of 1994 (M_L = 6.4, R = 32 km). The building is represented as a layered shear beam and a torsional shaft, characterized by the corresponding velocities of vertically propagating waves through the structure. The previously introduced waveform inversion algorithm is applied, which fits in the least squares sense pulses in low‐pass filtered impulse response functions computed at different stories. This paper demonstrates that layered shear beam and torsional shaft models are valid for this building, within bands that include the first five modes of vibration for each of the North–South (NS), East–West (EW), and torsional responses (0–1.7 Hz for NS and EW, and 0–3.5 Hz for the torsional response). The observed pulse travel time from ground floor to penthouse level is τ ≈1.5 s for NS and EW and τ ≈ 0.9 s for the torsional responses. The identified equivalent uniform shear beam wave velocities are β_eq ≈ 140 m/s for NS and EW responses, and 260 m/s for torsion, and the apparent Q ≈ 25 for the NS and torsional, and ≈14 for the EW response. Across the layers, the wave velocity varied 90–170 m/s for the NS, 80–180 m/s for the EW, and 170–350 m/s for the torsional responses. The identification method is intended for use in structural health monitoring. Copyright © 2013 John Wiley & Sons, Ltd.     

南水北调总干渠中线工程豫北段基础稳定性的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。根据监测点下沉拟合曲线预测,调水干渠多数段落基础下沉趋稳;南水北调中线工程豫北段受多种因素叠加影响,总体基础稳定性较差,地质构造活动、城市地表沉降为主要影响因素,基础施工为次要影响,矿区开采沉陷对基础失稳未见直接相关性。