Retrieval of super‐virtual refraction by cross‐correlation

Topography and severe variations of near‐surface layers lead to travel‐time perturbations for the events in seismic exploration. Usually, these perturbations could be estimated and eliminated by refraction technology. The virtual refraction method is a relatively new technique for retrieval of refraction information from seismic records contaminated by noise. Based on the virtual refraction, this paper proposes super‐virtual refraction interferometry by cross‐correlation to retrieve refraction wavefields by summing the cross‐correlation of raw refraction wavefields and virtual refraction wavefields over all receivers located outside the retrieved source and receiver pair. This method can enhance refraction signal gradually as the source–receiver offset decreases. For further enhancement of refracted waves, a scheme of hybrid virtual refraction wavefields is applied by stacking of correlation‐type and convolution‐type super‐virtual refractions. Our new method does not need any information about the near‐surface velocity model, which can solve the problem of directly unmeasured virtual refraction energy from the virtual source at the surface, and extend the acquisition aperture to its maximum extent in raw seismic records. It can also reduce random noise influence in raw seismic records effectively and improve refracted waves’ signal‐to‐noise ratio by a factor proportional to the square root of the number of receivers positioned at stationary‐phase points, based on the improvement of virtual refraction's signal‐to‐noise ratio. Using results from synthetic and field data, we show that our new method is effective to retrieve refraction information from raw seismic records and improve the accuracy of first‐arrival picks.     

Virtual signals combination – phase analysis for 2D and 3D data representation (acoustic medium)

We formulate the Kirchhoff‐Helmholtz representation theory for the combination of seismic interferometry signals synthesized by cross‐correlation and by cross‐convolution in acoustic media. The approach estimates the phase of the virtual reflections from the boundary encompassing a volume of interest and subtracts these virtual reflections from the total seismic‐interferometry wavefield. The reliability of the combination result, relevant for seismic exploration, depends on the stationary‐phase and local completeness in partial coverage regions. The analysis shows the differences in the phase of the corresponding seismic interferometry (by cross‐correlation) and virtual reflector (by cross‐convolution) signals obtained by 2D and 3D formulations, with synthetic examples performed to remove water layer multiples in ocean bottom seismic (OBS) acoustic data.     

基于逆虚折射干涉法有效提取近地表弱地震信号

在地震勘探中,地形起伏和近地表速度的剧烈变化会导致地震波旅行时的扰动,通常会通过折射波信息来估算和消除这些扰动.本文在虚折射的基础上提出了逆虚折射干涉法,通过虚折射波场和原始折射波场的互相关,并对所有位于固定相位点上的检波点进行叠加,重构出逆虚折射波场.通过逆虚折射与超级虚折射的叠加,保证了不同偏移距下折射波振幅恢复的一致性,显著提高折射波的信噪比,有效提取弱信号.同时,本文采用反褶积干涉法来压制由于互相关和褶积产生的子波旁瓣的影响,弥补低频和高频能量的损失,改善恢复的折射波场的稳定性和分辨率.该新方法不需要知道近地表复杂速度模型的信息,可以将虚折射的勘探孔径恢复到原始地震记录的最大孔径.合成资料和实际资料的计算结果表明,基于反褶积的逆虚折射干涉法能够从低信噪比的资料中,有效恢复出折射波信息.     

Drill‐bit seismic interferometry with and without pilot signals*

We use different interferometry approaches to process the seismic signals generated by a drill‐bit source in one well and recorded by seismic receivers located both in a second borehole and at the surface near the source well. We compare the standard interferometry results, obtained by using the raw drill‐bit data without a pilot signal, with the new interferometry results obtained by using the drill‐bit seismograms correlated with a reference pilot signal. The analysis of the stationary phase shows that the final results have different S/N levels and are affected by the coherent noise in the form of rig arrivals. The interferometry methods are compared by using different deconvolution approaches. The analysis shows that the results agree with the conventional drill‐bit seismograms and that using the reference pilot signal improves the quality of the drill‐bit wavefields redatumed by the interferometry method.     

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.     

Synthesis of a seismic virtual reflector*

We describe a method to process the seismic data generated by a plurality of sources and registered by an appropriate distribution of receivers, which provides new seismic signals as if in the position of the receivers (or sources) there was an ideal reflector, even if this reflector is not present there. The data provided by this method represent the signals of a virtual reflector. The proposed algorithm performs the convolution and the subsequent sum of the real traces without needing subsurface model information. The approach can be used in combination with seismic interferometry to separate wavefields and process the reflection events. The application is described with synthetic examples, including stationary phase analysis and with real data in which the virtual reflector signal can be appreciated.     

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.     

Application of the virtual refraction to near‐surface characterization at the Boise Hydrogeophysical Research Site‡

Seismic interferometry is a relatively new technique to estimate the Green's function between receivers. Spurious energy, not part of the true Green's function, is produced because assumptions are commonly violated when applying seismic interferometry to field data. Instead of attempting to suppress all spurious energy, we show how spurious energy associated with refractions contains information about the subsurface in field data collected at the Boise Hydrogeophysical Research Site. By forming a virtual shot record we suppress uncorrelated noise and produce a virtual refraction that intercepts zero offset at zero time. These two features make the virtual refraction easy to pick, providing an estimate of refractor velocity. To obtain the physical parameters of the layer above the refractor we analyse the cross‐correlation of wavefields recorded at two receivers for all sources. A stationary‐phase point associated with the correlation between the reflected wave and refracted wave from the interface identifies the critical offset. By combining information from the virtual shot record, the correlation gather and the real shot record we determine the seismic velocities of the unsaturated and saturated sands, as well as the variable relative depth to the water‐table. Finally, we discuss how this method can be extended to more complex geologic models.     

Retrieving virtual reflection responses at drill‐bit positions using seismic interferometry with drill‐bit noise

In the field of seismic interferometry, researchers have retrieved surface waves and body waves by cross‐correlating recordings of uncorrelated noise sources to extract useful subsurface information. The retrieved wavefields in most applications are between receivers. When the positions of the noise sources are known, inter‐source interferometry can be applied to retrieve the wavefields between sources, thus turning sources into virtual receivers. Previous applications of this form of interferometry assume impulsive point sources or transient sources with similar signatures. We investigate the requirements of applying inter‐source seismic interferometry using non‐transient noise sources with known positions to retrieve reflection responses at those positions and show the results using synthetic drilling noise as source. We show that, if pilot signals (estimates of the drill‐bit signals) are not available, it is required that the drill‐bit signals are the same and that the phases of the virtual reflections at drill‐bit positions can be retrieved by deconvolution interferometry or by cross‐coherence interferometry. Further, for this case, classic interferometry by cross‐correlation can be used if the source power spectrum can be estimated. If pilot signals are available, virtual reflection responses can be obtained by first using standard seismic‐while‐drilling processing techniques such as pilot cross‐correlation and pilot deconvolution to remove the drill‐bit signatures in the data and then applying cross‐correlation interferometry. Therefore, provided that pilot signals are reliable, drill‐bit data can be redatumed from surface to borehole depths using this inter‐source interferometry approach without any velocity information of the medium, and we show that a well‐positioned image below the borehole can be obtained using interferometrically redatumed reflection responses with just a simple velocity model. We discuss some of the practical hurdles that restrict the application of the proposed method offshore.     

3D seismic residual statics solutions derived from refraction interferometry

We apply interferometric theory to solve a three‐dimensional seismic residual statics problem to improve reflection imaging. The approach calculates the static solutions without picking the first arrivals from the shot or receiver gathers. The static correction accuracy can be significantly improved by utilising stacked virtual refraction gathers in the calculations. Shots and receivers may be placed at any position in a three‐dimensional seismic land survey. Therefore, it is difficult to determine stationary shots and receivers to form the virtual refraction traces that have identical arrival times, as in a two‐dimensional scenario. To overcome this problem, we use a three‐dimensional super‐virtual interferometry method for residual static calculations. The virtual refraction for a stationary shot/receiver pair is obtained via an integral along the receiver/shot lines, which does not require knowledge of the stationary locations. We pick the maximum energy times on the interferometric stacks and solve a set of linear equations to derive reliable residual static solutions. We further apply the approach to both synthetic and real data.     

Interferometric interpolation of sparse marine data

We present the theory and numerical results for interferometrically interpolating 2D and 3D marine surface seismic profiles data. For the interpolation of seismic data we use the combination of a recorded Green's function and a model‐based Green's function for a water‐layer model. Synthetic (2D and 3D) and field (2D) results show that the seismic data with sparse receiver intervals can be accurately interpolated to smaller intervals using multiples in the data. An up‐ and downgoing separation of both recorded and model‐based Green's functions can help in minimizing artefacts in a virtual shot gather. If the up‐ and downgoing separation is not possible, noticeable artefacts will be generated in the virtual shot gather. As a partial remedy we iteratively use a non‐stationary 1D multi‐channel matching filter with the interpolated data. Results suggest that a sparse marine seismic survey can yield more information about reflectors if traces are interpolated by interferometry. Comparing our results to those of f‐k interpolation shows that the synthetic example gives comparable results while the field example shows better interpolation quality for the interferometric method.     

Consistency of the spatial autocorrelation method with seismic interferometry and its consequence

We have cross‐checked the conventional theory of the spatial autocorrelation method and the consequence of seismic interferometry: the retrieval of the elastodynamic Green's function. Their mutual consistency is almost complete. The basic formulas of the conventional spatial autocorrelation theory can be derived by an alternative approach based on the retrieval of the elastodynamic Green's function. The only discrepancy is found with the average of the complex coherence function over azimuth in a wavefield dependent on azimuth. It is hypothesized, in discussion, that this discrepancy is due to the way of representing the wavefield in the background theory of seismic interferometry that can produce only wavefields moderately dependent on azimuth and that the mentioned consequence of seismic interferometry can also only make sense in a wavefield moderately dependent on azimuth. Our field experiment with a wavefield dependent on azimuth showed that the consequence of seismic interferometry in the logical framework of the conventional spatial autocorrelation theory is appropriate under such degrees of approximation as the measure proposed in this study, i.e., the deviation of the total dispersion curves is between about 10 and 16 per cent at the maximum from those averaged over azimuth. The acceptance of the retrieval of Green's function gives a proper physical meaning to the complex coherence function: the real part of the elastodynamic Green's function normalized by its zero‐offset version. This makes it possible to take a deterministic approach rather than the statistical one on which the conventional spatial autocorrelation method is based and gives fruitful new aspects and perspectives. For example, the formula for the multi‐mode case is given and the possibility of exploration of two or three dimensional velocity structures is suggested.     

Retrieval of local surface wave velocities from traffic noise – an example from the La Barge basin (Wyoming)

In regions where active source seismic exploration is constrained by limitations of energy penetration and recovery, cost and logistical concerns, or regulatory restrictions, analysis of natural source seismic data may provide an alternative. In this study, we investigate the feasibility of using locally‐generated seismic noise in the 2–6 Hz band to obtain a subsurface model via interferometric analysis. We apply this technique to three‐component data recorded during the La Barge Passive Seismic Experiment, a local deployment in south‐western Wyoming that recorded continuous seismic data between November 2008 and June 2009. We find traffic noise from a nearby state road to be the dominant source of surface waves recorded on the array and observe surface wave arrivals associated with this source up to distances of 5 kms. The orientation of the road with respect to the deployment ensures a large number of stationary points, leading to clear observations on both in‐line and cross‐line virtual source‐receiver pairs. This results in a large number of usable interferograms, which in turn enables the application of standard active source processing methods like signal processing, common offset stacking and traveltime inversion. We investigate the dependency of the interferograms on the amount of data, on a range of processing parameters and on the choice of the interferometry algorithm. The obtained interferograms exhibit a high signal‐to‐noise ratio on all three components. Rotation of the horizontal components to the radial/transverse direction facilitates the separation of Rayleigh and Love waves. Though the narrow frequency spectrum of the surface waves prevents the inversion for depth‐dependent shear‐wave velocities, we are able to map the arrival times of the surface waves to laterally varying group and phase velocities for both Rayleigh and Love waves. Our results correlate well with the known geological structure. We outline a scheme for obtaining localized surface wave velocities from local noise sources and show how the processing of passive data benefits from a combination with well‐established exploration seismology methods. We highlight the differences with interferometry applied to crustal scale data and conclude with recommendations for similar deployments.     

Path-integral seismic imaging

A new type of seismic imaging, based on Feynman path integrals for waveform modelling, is capable of producing accurate subsurface images without any need for a reference velocity model. Instead of the usual optimization for traveltime curves with maximal signal semblance, a weighted summation over all representative curves avoids the need for velocity analysis, with its common difficulties of subjective and time‐consuming manual picking. The summation over all curves includes the stationary one that plays a preferential role in classical imaging schemes, but also multiple stationary curves when they exist. Moreover, the weighted summation over all curves also accounts for non‐uniqueness and uncertainty in the stacking/migration velocities. The path‐integral imaging can be applied to stacking to zero‐offset and to time and depth migration. In all these cases, a properly defined weighting function plays a vital role: to emphasize contributions from traveltime curves close to the optimal one and to suppress contributions from unrealistic curves. The path‐integral method is an authentic macromodel‐independent technique in the sense that there is strictly no parameter optimization or estimation involved. Development is still in its initial stage, and several conceptual and implementation issues are yet to be solved. However, application to synthetic and real data examples shows that it has the potential for becoming a fully automatic imaging technique.     

Automatic stacking‐velocity estimation using similarity‐weighted clustering

Local seismic event slopes contain subsurface velocity information and can be used to estimate seismic stacking velocity. In this paper, we propose a novel approach to estimate the stacking velocity automatically from seismic reflection data using similarity‐weighted k‐means clustering, in which the weights are local similarity between each trace in common midpoint gather and a reference trace. Local similarity reflects the local signal‐to‐noise ratio in common midpoint gather. We select the data points with high signal‐to‐noise ratio to be used in the velocity estimation with large weights in mapped traveltime and velocity domain by similarity‐weighted k‐means clustering with thresholding. By using weighted k‐means clustering, we make clustering centroids closer to those data points with large weights, which are more reliable and have higher signal‐to‐noise ratio. The interpolation is used to obtain the whole velocity volume after we have got velocity points calculated by weighted k‐means clustering. Using the proposed method, one obtains a more accurate estimate of the stacking velocity because the similarity‐based weighting in clustering takes into account the signal‐to‐noise ratio and reliability of different data points in mapped traveltime and velocity domain. In order to demonstrate that, we apply the proposed method to synthetic and field data examples, and the resulting images are of higher quality when compared with the ones obtained using existing methods.     

Depth imaging in anisotropic media by symmetric non-stationary phase shift

We present a new depth‐imaging method for seismic data in heterogeneous anisotropic media. This recursive explicit method uses a non‐stationary extrapolation operator to allow lateral velocity variation, and it uses the relationship between phase angle and the spectral coordinates of seismic data to allow velocity variation with phase angle. A qualitative comparison of migration impulse responses suggests that, for an equivalent cost, the symmetric non‐stationary phase‐shift (SNPS) operator is superior to the phase‐shift plus interpolation (PSPI) operator, for very large depth intervals. To demonstrate the potential of the new method, seismic data from a physical model acquired over a transversely isotropic medium are imaged using a shot‐record migration based on the SNPS operator.     

Simulations of non‐stationary frequency content and its importance to seismic assessment of structures

To verify the importance of the non‐stationary frequency characteristic of seismic ground motion, a joint time–frequency analysis technique of time signals, called chirplet‐based signal approximation, is developed to extract the non‐stationary frequency information from the recorded data. The chirplet‐based signal approximation is clear in concept, similar to Fourier Transform in mathematical expressions but with different base functions. Case studies show that the chirplet‐based signal approximation can represent the joint time–frequency variation of seismic ground motion quite well. Both the random models of uniform modulating process and evolutionary process are employed to generate artificial seismic waves. The joint time–frequency modulating function in the random model of evolutionary process is determined by chirplet‐based signal approximation. Finally, non‐linear response analysis of a SODF system and a frame structure is performed based on the generated artificial seismic waves. The results show that the non‐stationary frequency characteristic of seismic ground motion can significantly change the non‐linear response characteristics of structures, particularly when a structure goes into collapse phase under seismic action. It is concluded that non‐stationary frequency characteristic of seismic ground motion should be considered for the assessment of seismic capacity of structures. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

不同参数对地震干涉实现VSP向SWP波场重构的影响

地震干涉是近几年勘探地球物理领域研究的一个热点方向,它是波场重构的有力工具.然而,地震干涉往往引入虚假同相轴,影响波场重构的质量.为进一步分析虚假同相轴产生的原因并改善波场重构的质量,本文基于稳相分析,详细探讨子波主频、激发炮数、检波器埋深、检波器间距和地层倾角五个因素对波场重构的影响.模型结果表明,震检排列方式和地层倾角等因素通过改变有限震检范围内稳相点的位置来影响波场重构的质量.将垂直地震剖面VSP中的下行直达波与下行反射波进行地震干涉处理,可有效重构来自高陡反射面的下行反射波,将传统的VSP转化为单井地震剖面SWP.直接利用重构SWP波场进行成像,不仅扩大传统VSP的成像范围,而且避免常规勘探面临的静校正和近地表速度建模问题,为高陡构造成像提供一种新方法.     

高铁地震数据干涉成像技术初探

高铁运行会引起铁轨的震动,从而产生地震波向地下介质中传播,通过研究该地震波可对高铁沿线的地质情况进行持续监测.与常规地震勘探中的震源相比,高铁地震中的震源较为复杂,为移动震源,而地震干涉技术可以通过地震记录间的相互干涉,消除震源的影响,因此可利用地震干涉技术对高铁地震信号进行处理并成像.本文通过分析研究,总结出地震干涉方法在处理高铁地震数据时的关键技术问题:不同于常规地震干涉中先干涉后叠加的干涉成像方式,高铁地震移动源的特点使得干涉顺序变为先叠加后干涉,由此带入了大量震源串扰噪声;初步提出两种解决高铁地震干涉成像的思路:通过对高铁地震信号的处理,使高铁变相"提速"或"降速",给出了"提速"或"降速"后各自的成像思路,并给出了数据处理的技术设想.