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.     

Gravimetric monitoring of the first field‐wide steam injection in a fractured carbonate field in Oman – a feasibility study

Gas‐Oil Gravity Drainage is to be enhanced by steam injection in a highly fractured, low permeability carbonate field in Oman. Following a successful pilot, field‐wide steam injection is being implemented, first of this type in the world. A dedicated monitoring program has been designed to track changes in the reservoir. Various observations are to be acquired, including, surface deformation, temperature measurements, microseismic, well logs, pressure and saturation measurements to monitor the reservoir. Because significant changes in the reservoir density are expected, time‐lapse gravimetry is also being considered. In this paper we investigate the feasibility of gravimetric monitoring of the thermally enhanced gravity drainage process at the carbonate field in Oman. For this purpose, forward gravity modelling is performed. Based on field groundwater measurements, the estimates of the hydrological signal are considered and it is investigated under what conditions the groundwater influences can be minimized. Using regularized inversion of synthetic gravity data, we analyse the achievable accuracy of heat‐front position estimates. In case of large groundwater variations at the field, the gravity observations can be significantly affected and, consequently, the accuracy of heat‐front monitoring can be deteriorated. We show that, by applying gravity corrections based on local observations of groundwater, the hydrological influences can to a large extent be reduced and the accuracy of estimates can be improved. We conclude that gravimetric monitoring of the heat‐front evolution has a great potential.     

On the Influence of Sea Surface Height Variability on Satellite Altimeter Derived Gravity

The effect of sea surface height (SSH) variability is one of the primary factors that limit the accuracy and resolution of altimeter-derived gravity values. We propose a method to estimate the influence of variation of the sea surface height on the accuracy of satellite-derived gravity by simulation technique, with a case study around Indonesian waters. Wederived an Indonesian marine gravity map using the Geosat-geodetic mission (GM). Since most of the area studied is located around coastal and shallow areas, the measurement of SSH of this area is less accurate. To obtain a distribution of SSH variability over the study area, Topex/Poseidon (T/P) data were first processed and assessed. Processing 52 cycles of the Topex/Poseidon data, the root mean square (RMS) of SSH variability for each cycle was found to vary from 1 to 179 cm. Further, for the purpose of estimating the accuracy of altimeter-derived gravity, we derived several levels of Gaussian noise, computed simulation data by adding the Gaussian noise to Geosat data, and determined simulated gravity maps. Based on the distribution of RMS values from T/P data and standard deviation (STD) differences between the simulated and the original gravity maps, we estimated the accuracy of the gravity map. Around Indonesian waters, the accuracy of the gravity map influenced by SSH variation was estimated to be within the range 0.8~93 mgal.     

利用噪声记录估计福建地区中上地壳体波速度结构

In this paper,the dispersion curves of the Rayleigh wave and Love wave were extracted from the seismic noise records of 25 broadband stations of the Fujian Seismic Network, and inverted for the lithosphere velocity structure. Furthermore,the velocity model was verified by the seismic explosion observations. Our results indicate that the resolution of the lithosphere velocity structure obtained by this method is good in the shallow part,but in the deep part,inversion accuracy for the wave velocity structure is low,which is caused mainly by the small inter-station distance chosen in the paper. Thus the wave dispersion curves have high accuracy in the short-period part,but the warp of the wave dispersion curve in long-period part is large. Considering the results from both the noise inversion and the traditional inversion,we finally present a new velocity model,and the theoretical travel time calculated with the new model matches the explosion travel time very well.     

新疆尼勒克、巩留县交界6.0级地震测震应急流动台站架设及经验探讨

对新疆尼勒克、巩留交界6.0级地震测震应急流动台站的架设工作进行总结,针对应急流动台站的台址勘选、仪器设备的安装架设形式及数据实时传输等方面进行分析,总结经验,发现不足.期望对今后的地震现场应急流动测震台站的安装架设起到一定的参考作用.     

Water‐bottom multiple attenuation by Kirchhoff extrapolation

Despite being less general than 3D surface‐related multiple elimination (3D‐SRME), multiple prediction based on wavefield extrapolation can still be of interest, because it is less CPU and I/O demanding than 3D‐SRME and moreover it does not require any prior data regularization. Here we propose a fast implementation of water‐bottom multiple prediction that uses the Kirchhoff formulation of wavefield extrapolation. With wavefield extrapolation multiple prediction is usually obtained through the cascade of two extrapolation steps. Actually by applying the Fermat’s principle (i.e., minimum reflection traveltime) we show that the cascade of two operators can be replaced by a single approximated extrapolation step. The approximation holds as long as the water bottom is not too complex. Indeed the proposed approach has proved to work well on synthetic and field data when the water bottom is such that wavefront triplications are negligible, as happens in many practical situations.     

An automated cross‐correlation based event detection technique and its application to a surface passive data set

In studies on heavy oil, shale reservoirs, tight gas and enhanced geothermal systems, the use of surface passive seismic data to monitor induced microseismicity due to the fluid flow in the subsurface is becoming more common. However, in most studies passive seismic records contain days and months of data and manually analysing the data can be expensive and inaccurate. Moreover, in the presence of noise, detecting the arrival of weak microseismic events becomes challenging. Hence, the use of an automated, accurate and computationally fast technique for event detection in passive seismic data is essential. The conventional automatic event identification algorithm computes a running‐window energy ratio of the short‐term average to the long‐term average of the passive seismic data for each trace. We show that for the common case of a low signal‐to‐noise ratio in surface passive records, the conventional method is not sufficiently effective at event identification. Here, we extend the conventional algorithm by introducing a technique that is based on the cross‐correlation of the energy ratios computed by the conventional method. With our technique we can measure the similarities amongst the computed energy ratios at different traces. Our approach is successful at improving the detectability of events with a low signal‐to‐noise ratio that are not detectable with the conventional algorithm. Also, our algorithm has the advantage to identify if an event is common to all stations (a regional event) or to a limited number of stations (a local event). We provide examples of applying our technique to synthetic data and a field surface passive data set recorded at a geothermal site.     

甘肃嘉峪关地震台大地电场观测数据干扰分析

文章对甘肃嘉峪关地震台2001年至2009年以来大地电场分钟值资料做了全面的分析和整理。对造成观测资料相关性低的电暴（磁暴）、人类活动、雷电、降水、观测系统本身变化（包括更换电极）、电阻率测量、公用接地网等干扰因素进行分析,并做出相关的解释。这将对分析和运用地电场资料进行地震预报有一定的意义。     

首都圈数字地震台网对微弱爆破信号的检测能力

利用首都圈数字地震台网接收人工地震信号,进行地下结构研究具有重要意义.但人工震源释放的能量小,激发的地震波以短周期为主,因此本文较全面地研究了地震台网对短周期微弱信号（1～20 Hz）的检测能力：(1) 分析了台网的背景噪声,结果表明基岩台址的地震台噪声比沉积盖层台址的地震台噪声低约13 dB,这相当于近1个震级的检测阈值;夜间的噪声比白天低约5 dB;噪声有逐年增高的趋势,2006年比2001年噪声提高约4 dB.(2 )分析了在台网内进行的药量为25 kg的陆地井下爆破实验,一次爆破相当于0.69级（M_L）的天然地震,有18个地震台可辨认爆炸产生的Pg、Pm或Pc波;离爆破点218 km的基岩台,仍可以接收到振幅只有1.6 nm 的Pm波,这个结果可为地震勘探实际工作提供参考.(3) 研究了台网外核爆试验的信号特征,2006年发生在朝鲜的地下核试验是一次检验台网检测微弱信号能力的好机会.波形记录经1～5Hz滤波后,台网中噪声小的18个基岩台可以清晰辨认核爆破产生的P波或Lg波,P波平均振幅为16 nm,计算的平均震级为m_b4.3,和NEIC给出的震级相同;分析还表明背景噪声是影响台站信号检测能力的主要因素之一.     

An indication of the soil topmost layer response in Quito (Ecuador) using noise H/V spectral ratio

H/V noise spectral ratio (HVSR), standard spectral ratio (SSR), and receiver functions (RF) techniques have been used in the Quito (Ecuador) urban area to estimate the frequency dependence of soil response. Two amplified frequencies obtained by the HVSR method appear on about 60 sites. Taking into account the most amplified frequency rather than considering only the first amplified frequency, generally associated with the fundamental frequency, we find that iso-frequency curves tightly fit the surface geology. The second amplified frequency is interpreted as the fundamental frequency of the soft thin topmost layer, which in some cases amplifies the surface ground motion more than the rest of the soil column. This hypothesis is further supported by the results provided from the SSR and RF studies at a station located on top of a solid waste landfill, and by a study of known thickness of a waste landfill, using the HVSR method.