Break through the limits of vibroseis data quality

The vibroseis method has become the principal data acquisition method in land seismic exploration. It seems that this method has been extended to its limits as the search for energy resources continues. Many practical issues arising from field operations have remained theoretically unexplained, for example, variations in wavelet arrival time, inaccurate wavelet estimation and harmonics in the wavelet itself. The focus of this paper is the proposal of a new model, which is referred to as the vibrator‐coupled ground model, to simulate the filtering effects of a complex coupling system consisting of the coupling between the baseplate and the ground as well as the coupling between the captured ground mass near the vibrator baseplate and the surrounding earth. With this vibrator‐coupled ground model many of the practical issues mentioned above were reasonably addressed. Furthermore, it was demonstrated from experimental tests that both the pilot sweep and the weighted‐sum groundforce, when filtered by the vibrator‐coupled ground model, are proportional to the far‐field particle velocity whereas the unfiltered signals are not. The harmonics on the filtered weighted‐sum groundforce successfully maintain a proportional relationship with the harmonics seen in the far‐field signal.     

SIMULTANEOUS VIBROSEIS RECORDING1

An experiment was undertaken at BP's Fulbeck Geophysical test site to compare the viability of various simultaneous vibroseis recording techniques, which are often recommended as a means of improving data acquisition production rates for 3D seismic surveys. Of particular interest were: (a) the ability to separate the signals from each source during processing, (b) the generation and suppression of harmonics and (c) the effects of any source interaction. Two vibrators were deployed with a baseplate separation of 10 m, about a borehole containing a vertical array of geophones. Our analysis concentrated on the groundforce signals measured at each vibrator and the far-field signatures measured using a vertical geo-phone at a depth of 204 m. By comparing single vibrator records with similar but separated records from a simultaneous recording sequence, signal separability, harmonic suppression and vibrator interaction could be fully studied. Separated far-field signatures from simultaneous vibroseis methods using combinations of up and downsweeps exhibited unsuppressed harmonics and substantial energy from the undesired source which leaked through the correlation process. The ‘up/down’ method was capable of separating the signal from each source by only 12.7 dB, and is therefore unsuitable as a field technique. The variphase simultaneous vibroseis methods studied afforded some harmonic suppression and gave signal separations of about 30.0 dB. Use of variphase simultaneous vibroseis methods will compromise the quality of the data recorded, when compared with single-source acquisition methods. None of the simultaneous vibroseis methods tested provided adequate signal separation and, therefore, cannot be recommended as data acquisition techniques. The ‘alternate sweeping’ method coupled with multispread recording will give the desired improvement in data acquisition rates, while preserving the necessary quality of our seismic data.     

Radiation from flexural vibrations of the baseplate and their effect on the accuracy of traveltime measurements

We analyse the problem of radiation of seismic waves by a vibroseis source when the baseplate is subject to flexure. A theoretical model is proposed to account for baseplate flexure, generalizing the well‐known model of the vibroseis source of Sallas and Weber, which was developed for a rigid plate. Using the model proposed, we analyse the effect of flexure on the properties of seismic waves. We show that the flexure does not contribute to the far‐field and mainly affects the readings of the reference accelerometer that is used to measure the force applied to the ground; these readings generally become dependent on the location of the sensor on the plate. For muddy and sandy soils, the effect of flexure on baseplate‐acceleration measurements is nonetheless pronounced at the high end of the vibroseis frequency band only (～100 Hz), and is negligible at all frequencies for stiffer soils. The corresponding phase lags introduced by the flexural vibrations at high frequencies lead to errors in the traveltime measurements (through the cross‐correlation function) of up to 0.6 ms for muddy soils and less for denser soils. We show the existence of an optimal position of the reference sensor on the baseplate and also propose a general method of eliminating the phase lag due to the baseplate flexure in acceleration measurements.     

Review of vibroseis data acquisition and processing for better amplitudes: adjusting the sweep and deconvolving for the time-derivative of the true groundforce

The goal of vibroseis data acquisition and processing is to produce seismic reflection data with a known spatially-invariant wavelet, preferably zero phase, such that any variations in the data can be attributed to variations in geology. In current practice the vibrator control system is designed to make the estimated groundforce equal to the sweep and the resulting particle velocity data are cross-correlated with the sweep. Since the downgoing far-field particle velocity signal is proportional to the time-derivative of the groundforce, it makes more sense to cross-correlate with the time-derivative of the sweep. It also follows that the ideal amplitude spectrum of the groundforce should be inversely proportional to frequency. Because of non-linearities in the vibrator, bending of the baseplate and variable coupling of the baseplate to the ground, the true groundforce is not equal to the pre-determined sweep and varies not only from vibrator point to vibrator point but also from sweep to sweep at each vibrator point. To achieve the goal of a spatially-invariant wavelet, these variations should be removed by signature deconvolution, converting the wavelet to a much shorter zero-phase wavelet but with the same bandwidth and signal-to-noise ratio as the original data. This can be done only if the true groundforce is known. The principle may be applied to an array of vibrators by employing pulse coding techniques and separating responses to individual vibrators in the frequency domain. Various approaches to improve the estimate of the true groundforce have been proposed or are under development; current methods are at best approximate.     

ON THE DETERMINATION OF THE DOWNGOING P-WAVES RADIATED BY THE VERTICAL SEISMIC VIBRATOR*

It is well recognized that in order to realize the full potential of the Vibroseis technique, one needs to ensure accurate phase locking and a meaningful cross-correlation. To achieve these two important objectives we require an accurate estimate of the compressional stress wave radiated by the vibrator into the ground. In this paper a simple method (subject of a patent application) is developed for predicting the compressional stress waves radiated by a vertical vibrator. The main feature of the proposed method is that it involves the field measurement of the acceleration of the reaction mass and the baseplate, respectively. The method is illustrated by computing the compressional stress waves generated by a typical vertical vibrator radiating into ice, chalk, sand, and mud. It is shown that for a seismic vibrator radiating into hard ground the pressure of the downgoing P-wave is 180° out of phase with the baseplate velocity. It is also shown that when the driving force of the seismic vibrator has a flat amplitude spectrum, the amplitude spectrum of the downgoing P-wave falls off by 6 dB/octave towards low frequencies.     

Geophysical benefits from an improved seismic vibrator

The seismic vibrator has become a very important source for land data acquisition and there have been dramatic improvements in recent times in the application of the vibroseis technique. These improvements have led to much increased productivity and in many cases much denser source sampling. At the same time, the vibrator itself has seen little improvement over the last couple of decades. There are needs in a few areas where an improvement in the vibrator itself can bring benefits to the quality of the seismic data acquired. This paper describes progress in four such areas, low‐frequency performance, high‐frequency performance, an improved estimate of the vibrator groundforce and source signature consistency over variable ground conditions. Each of these vibrator characteristics will be discussed in turn. Meanwhile, two field test results in which the performance of two different vibrators in these four areas are compared.     

Developments in vibrator control

Hydraulic limitations, non-rigidity of the baseplate as well as variable characteristics of the ground constantly distort the downgoing energy output by vibrators. Therefore, a real time feedback control must be performed to continuously adjust the emitted force to the reference pilot signal. This ground force is represented by the weighted sum of the reaction mass and the baseplate accelerations. It was first controlled with an amplitude and phase locked loop system, poorly reactive and sensitive to noise. Later on, new vibrator electronics based on a digital model of the vibrator were introduced. This model is based on the physical equations of the vibrator and of the ground. During an 'identification' process, the model is adjusted to each particular vibrator. Completed by a Kalman adaptive filter to remove the noise, it computes ten estimated states via a linear quadratic estimator. These states are used by a linear quadratic control to compute the torque motor input and to compare the ground force estimated from the states with the pilot signal. Test results using downhole geophones demonstrate the benefit of filtered mode operation.     

SEISMIC VIBRATOR MODELLING1

The wavefield in, and at the surface of, a homogeneous, isotropic, perfectly elastic half-space, excited by a traction distribution at the surface of the medium is investigated. The emitted wavefield is a spatial convolution of the surface tractions and the spatial impulse response. The properties of the wavefield in the far-field of the medium are derived and it is shown that the far-field particle velocity is essentially equal to a weighted sum of the time derivative of the integrated surface tractions, that is, of the components of the ‘ground force’. The theory is valid for an arbitrary geometry and orientation of the surface tractions, and is independent of the boundary conditions at the surface of the medium. The surface tractions are related to a source that consists of a mass distribution with an arbitrary force distribution imposed upon it. A boundary condition is introduced that accounts for the mass load and the forces applied to it but neglects vibrations within the mass. The boundary condition follows from the equation of motion of the surface mass load. The theory is applied to the Vibroseis configuration, using a P-wave vibrator model with a uniformly distributed force imposed on top of the baseplate, and assuming that horizontal surface traction components are absent. The distribution of displacement and stress directly underneath the baseplate of a single vibrator and an array of vibrators is investigated. Three different boundary conditions are used: (1) assuming uniform pressure, (2) assuming uniform displacement, (3) using the equation of motion of the baseplate as a boundary condition. The calculations of the distribution of stress and displacement over the plate for different elastic media and several frequencies of operation show that only the results obtained with the mixed boundary condition agree with measurements made in the field. The accuracy of three different phase-feedback signals is compared using synthetic data. Baseplate velocity phase-feedback leads to huge deviations in the determination of the far-field wavelet; reaction mass acceleration phase-feedback looks stable but neglects the differentiating earth filter; and phase-feedback to a weighted sum of baseplate and reaction mass accelerations becomes unstable with increasing frequency. The instability can be overcome using measurements over the whole baseplate. The model can be extended to a lossy layered earth.     

Nonlinear analysis on the coupling process of electromagnetic vibrator and earth

Vibrator is an excitation equipment of the vibratory source in the seismic exploration[1―3]. In order to ap-proach the δ function, dynamite can be used to release high energy in an instantaneous explosion and form the seismic wavelet. On the other hand, a portable vi-brator sweeps a low energy signal for a long time to simulate the function of an explosion source. The re-flection signal of the portable vibrator is estimated for time-delay by digital cross correlation technology, and the disp…     

The effects of inter‐vibrator interference on vibrator performance

The output from the hydraulic vibrators typically used for land seismic surveys is controlled by monitoring the acceleration measured by accelerometers mounted on the reaction mass and baseplate. The considerable energy output by such vibrators, which are coupled with the sensitivity of the accelerometers used, results in crosstalk if more than one vibrator is being used. In this paper, we present the results of a field experiment in which we measured the crosstalk between two adjacent vibrators. We found that the level of crosstalk was approximately ‐20 dB when the vibrators were adjacent but decreased with increasing frequency and separation. This result has implications for measurements of vibrator performance, source‐signature deconvolution, and in particular, estimates of the total energy output by a fleet of vibrators.     

消除滑动扫描地震数据中谐波干扰的反相关方法(英文)

滑动扫描技术是高效、高保真、环保的可控震源勘探技术之一,是下一组震源不必等待上一组震源震动结束即可开始震动的高效采集方法。该技术由于缩短了相邻两炮的等待时间,使得生产效率得到显著提高。但是后一炮的谐波畸变与前一炮的基波信号混叠在一起,不易分离,在相关后的地震记录上形成了严重的谐波干扰,降低了地震资料的质量。本文提出一种反相关方法来压制滑动扫描地震数据中的谐波干扰。该方法首先把地面力信号分解为基波和各阶谐波分量;然后将后一炮的相关前数据分别与各分量相关,只选取正时间轴中对应分量的自相关部分,利用各分量的反相关算子提取各阶谐波信息;最后从前一炮数据中减去提取出的高阶谐波,得到压制谐波后的地震记录。该方法对有效信号影响小,可同时处理相关前和相关后数据,而且算法简单稳定,计算效率高。本文分别对理论模型和实际数据进行处理,验证了该方法消除谐波干扰的有效性。     

土-结构相互作用效应对结构基底地震动影响的试验研究

利用土与结构动力相互作用振动台模型试验数据,通过各种试验工况下土层表面与基础表面加速度反应的比较,深入探讨了土与结构动力相互作用效应对高层建筑结构基底地震动的影响。从输入地震动频谱特性、输入地震动强度水平和上部结构动力特性3个方面详细分析了与SSI效应对高层建筑基底震动影响程度有关的一些因素。结果表明：SSI效应对高层建筑基底地震动的影响与输入地震波的动力特性有很大关系。在地震动的频谱成分方面,SSI效应对高层建筑基底地震动的影响主要体现为土层表面和基础表面在与输入地震动卓越频率相近处的频谱成分有较大差异;SSI效应对高层建筑基底地震动的影响程度随着输入加速度峰值水平的增加而减小;在某一特定地震波作用下,当上部结构的振动频率与地震地面运动的卓越频率相近时,SSI效应对高层建筑基底地震动的影响较为强烈。     

自贡市西山公园地形对地震动的影响

不规则地形和土层对地震动的影响较大,建(构)筑物选址及其抗震设防必须考虑地形和土层场地的放大效应,以避免或减轻其震害.利用自贡地形台阵记录的汶川地震(M_s8.0)的主震加速度时程,基于传统谱比法分别研究了地形和土层场地对地震动的放大效应.结果表明:(1)地形场地在低频的放大效应不明显,最大仅为1.24;在高频的放大效应较显著,在1~10 Hz频带,山顶的放大效应最大,EW、NS和UD地震动的最大放大效应分别为4.15、3.61和2.41,对应频率分别为5.72 Hz、6.46 Hz和7.44 Hz;在10~20 Hz频带,靠近山顶的山脊上某个地震动分量的放大效应最大,7#台站EW、5#台站NS和7#台站UD地震动的最大放大效应分别为9.10、5.56和2.52,对应频率分别为16.97 Hz、16.91 Hz和17.91 Hz.(2)地形场地的最大放大效应随高度有增加的趋势,且在0.1~10 Hz频带随高度增加的趋势更加明显.(3)土层场地水平向地震动在2 Hz以上开始明显放大,竖向地震动在4 Hz以上开始明显放大;EW、NS和UD地震动的最大放大效应分别为13.4、12.168和6.0,对应频率分别为6.94 Hz、7.55 Hz和10.8 Hz.(4)土层场地与地形场地的最大放大效应相比较,前者显著大于后者,对于水平向地震动,前者至少是后者的3倍以上;对于竖向地震动,前者至少是后者的2.5倍以上.(5)无论是地形场地还是土层场地,地震动的最大放大效应均有水平向大于竖向的特征.     

TWO METHODS FOR CONTINUOUS MONITORING OF HARMONIC DISTORTION IN VIBROSEIS SIGNALS1

The quality of Vibroseis survey data can be improved by continuously monitoring the vibrator's baseplate and reaction mass accelerations. Equipment failures can be detected as they occur, rather than relying on similarity trials at the beginning and end of the day's production. Equipment faults can then be corrected as they happen and thus would not have a detrimental effect on the quality of the survey data. Source efficiency can be optimized by monitoring the amount of harmonic distortion generated by the vibrator at different drive levels on the different surfaces which may be encountered during a survey. Phase problems introduced by poor coupling of the baseplate to the ground can also be identified and addressed in the field. Rapid analysis of vibrator signals is required if continuous monitoring is to be useful. Frequency-time (f-t) analyses of vibrator signals are often used in processing centres, but are slow and require a large storage capacity which makes the technique unsuitable for a field analysis system. The two methods proposed to analyse vibrator signals entail the use of hodograms and time-varying notch filters. Hodograms provide a qualitative analysis of harmonic distortion and vibrator performance. A fast, time-varying notch filter gives quantitative and qualitative information about the harmonic distortion present in the signal and can be used to identify problems with vibrator behaviour. Both the hodogram and fast, time-varying notch filter methods can analyse the vibrator's reaction mass and baseplate accelerations as it progresses through its sweep and can present automatically interpreted results to the operator before moving to the next vibrator point.     

Orthogonal vibroseis sweeps

Vibroseis is a method that imparts coded seismic energy into the ground. The energy is recorded with geophones and then processed using the known (coded) input signal. The resulting time‐domain representation of vibroseis data is an impulsive wavetrain with wavelet properties consistent with the coded input signal convolved with the earth's reflectivity series. Historically, vibratory seismic surveys collect data from one source location at a time, summing one or more sources at each location. We present a method of designing orthogonal sweeps using the concept of combisweeps. The orthogonal sweeps allow simultaneous recording and later separation of two or more unique source locations. Orthogonality of sweeps permits separation of the data into unique source‐location field records by a conventional correlation procedure. The separation power of the orthogonal sweeps is demonstrated by a comparison between separated data and data acquired with one vibrator. Separation noise was at a negligible level for our demonstration data sets when two vibrators were located 50 m to 200 m apart. Coincident generation and recording of two vibroseis sweeps at different locations would allow almost double the amount of data to be recorded for a given occupation time and requires only half the storage medium.     

On the accuracy of the ground force estimated in vibroseis acquisition

For a linear elastic Earth the time derivative of the ground force is considered proportional to the far-field wavelet. Under the assumption that the baseplate is stiff and the bending forces of the baseplate are negligible, the ground force is also approximated by the sum of the accelerations of the baseplate and the reaction mass weighted by the respective masses. Combining these two assumptions, the time derivative of the weighted sum is considered proportional to the far-field wavelet. This result, often referred to as the far-field wavelet assumption, although convenient and most often employed is not always valid. We explore its validity using the spectral harmonic ratios of recorded data, which are used extensively in data filtering and analysis of vibratory data. We show that the far-field wavelet assumption fails particularly for harmonic components of even order. More compact soil after repeated shots further invalidates this assumption. Non-linear modelling of the ground under the vibrator point may provide a direction towards solving this discrepancy. Finally, we describe a method for the estimation of the harmonic spectral ratios.     

Application of wavelets to analysis and simulation of earthquake motions

A method of applying wavelet transform to earthquake motion analysis is developed from the viewpoint of energy input structures, in which relationships between wavelet coefficients and energy input, namely energy principles in wavelet analysis are derived. By using the principles, time–frequency characteristics of the 1995 Hyogoken-Nanbu earthquake ground motions are analysed and time histories of energy input for various ranges of frequencies and epicentral distances are identified. Furthermore, a technique to simulate earthquake ground accelerations by wavelet inverse transform is developed on the condition that target time-frequency characteristics are specified. Structural responses to the simulated accelerations are compared with the target time–frequency characteristics, which shows satisfactory correlations between wavelet coefficients and energy responses in both time and frequency domains. Copyright © 1999 John Wiley & Sons, Ltd.     

Pushing the vibrator ground-force envelope towards low frequencies

Several mechanical and hydraulic limitations hinder the ground-force energy output of a seismic vibrator at low frequencies. The hydraulic pump flow, pump response time, reaction mass stroke, servo valve stroke, engine horsepower, accumulator size, harmonic distortion and vehicle chassis isolation each play a role in limiting the ground-force energy output of vibrators. In addition, the peak-decoupling force – which is defined as the smaller value of either the maximum peak force or the hold-down weight – also plays a role in limiting ground-force energy production. A model useful for simulating seismic vibrator dynamics is developed to evaluate the impact of these parameters on the vibrator fundamental force envelope at low frequencies. Model data show that among these factors the reaction mass stroke and the peak-decoupling force are key parameters for setting the target fundamental force that can be achieved at low frequencies. Formulas are derived to estimate fundamental force, peak force and the reaction mass displacement. These formulas can serve as guidelines for sweep designers who plan to design low frequency sweeps with considerable dwell time in the lower frequency ranges. Test data show that formulas can be used to profile the vibrator envelope at low frequencies.     

Effective elimination of subharmonic ghost events from vibroseis data

Harmonic or subharmonic noise is often present in vibroseis data as reverberation‐like, laterally coherent bands occurring parallel to and before or after, the main events. Such periodic noise is typically generated during the standard correlation process when the actual source signal travelling through the subsurface is, for whatever reason, different from the desired source signal, i.e., the pilot‐sweep controlling the baseplate and used for correlation. A typical cause can be that harmonic or subharmonic frequency partials are generated in addition to the vibroseis sweep's desired fundamental frequencies. These harmonics produce strong ‘ghost events’ during correlation of the geophone trace with the pilot‐sweep, originating from additional correlations between the fundamental and harmonic frequencies. Especially subharmonic ‘ghosts’ will overlap with ‘good’ fundamental signals, since for typically used up‐sweeps they are folded to later traveltimes, where the signal/noise‐ratio is already lower, thus aggravating or preventing a reliable interpretation of possible later reflections. Here, a method is introduced to remove these unwanted noise trains (with only negligible impact on the fundamental signal) by transforming the seismogram traces into a so‐called ‘(sub)harmonic domain’. In this domain, the respective harmonic noise portions are focused and separated from the fundamental signals, enabling easier detection and appropriate suppression. After back‐transformation to the x‐T domain, the records are free from the corresponding harmonic contamination and can then be processed as usual. The method operates in a data‐driven fashion, i.e., the traces are not uniformly processed but are processed depending upon their actual (sub)harmonic content. The decontamination procedure can be applied universally, i.e., to uncorrelated/correlated and/or vertically unstacked/stacked data either in a manual, semiautomated or fully automated manner. The method works perfectly for synthetic vibroseis traces with or without harmonic/subharmonic portions. The application to real, crustal‐scale vibroseis records that were acquired in 2006 in the Dead Sea region, Israel and that were severely contaminated by subharmonic ground‐roll ghosts covering reflectivity from the basement to the Moho, shows the robustness and success of the presented method.     

Comparison of 1D non-linear simulations to strong-motion observations: A case study in a swampy site of French Antilles (Pointe-à-Pitre,Guadeloupe)

Observations from many recent strong motion events have shown the importance of local soil conditions and non-linear soil behaviour on the seismic ground response (site effects). As demonstrated by previous seismic microzoning studies (Lebrun et al.) [1]), as well as by at least three historical major earthquakes, Pointe-à-Pitre is prone to strong site effects, due to the particular geology of the area. In this paper, we present a comparison between the strong-motion data available from the stations operating on the swampy site of Pointe-à-Pitre airport and the ground motions derived from 1D non-linear finite element simulations.Results show that, for moderate to strong ground motions, 1D simulations reproduce the main characteristics of site response in terms of duration, energy distribution, amplitude and frequency content. It also shows the importance of very superficial soft layers as peat or saturated mud in low frequency site effects simulations. This point is important for further engineering studies since such very soft formations overlain by stiffer landfills are commonly expected in the Antilles context. Our work also shows that Anderson's criteria, used to quantify the goodness-of-fit of simulated ground motions to the observed ones, appear to be an interesting diagnostic tool for testing the quality of numerical simulations from an engineering point of view.