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.     

参照标准点探测PS法

提出参照标准点探测PS法。既考虑PS的强散射特性,又顾及其稳定性,且初选时振幅阈值的设置参照每景影像中PS标准点的振幅值,精选时振幅离差阈值的设置也参照这些PS标准点的振幅离差值,都不需要人工设定。选取澳门半岛为实验区,以18景TerraSAR-X为数据源,进行PS探测实验。结果表明,相比于传统的PS探测方法,本方法能更快速、有效、可靠地探测出SAR影像中的PS。     

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.