ENHANCED SERVOVALVE TECHNOLOGY FOR SEISMIC VIBRATORS1

The Pelton DR^TM Servovalve Enhancement causes the natural output of a vibrator to resemble the desired output more closely. This simplifies the control problem and reduces harmonic distortion. The traditional type of servovalve used on seismic vibrators is a flow-control servovalve. Flow is proportional to a vibrator's baseplate velocity, with respect to its reaction mass. The new servovalve control parameter is pressure rather than flow. The differential pressure applied to a vibrator's actuator piston, multiplied by the area of the piston, equals the force applied to the vibrator's baseplate structure. This may be defined as actuator force. There is a simpler and more linear relationship between actuator force and ground force than between actuator velocity and ground force. Thus, it is better for the servovalve to control pressure into the actuator rather than flow. A flow-control servovalve can be made to control pressure by sensing the differential pressure across a vibrator's actuator piston and applying it as a negative feedback around the servovalve main stage. This has been carried out and tested. The result is more accurate vibrator control and reduced harmonic distortion.     

Modelling and modal analysis of seismic vibrator baseplate

The vibroseis method must be extended to its limits as the search for oil and gas continues on land. To successfully improve vibroseis data quality, it is crucial to evaluate each element in the vibroseis data acquisition system and ensure that the contribution from each element is successful. Vibroseis systems depend greatly upon the ability of vibrators to generate synchronous, repeatable ground-force sweeps over a broad frequency range. This requires that the reaction mass and the baseplate of the vibrator move as rigid bodies. However, rigid-body motion is not completely true for high- frequency vibrations, especially for the vibrator baseplate. In order to accurately understand the motion of the vibrator baseplate, a finite element analysis model of the vibrator baseplate and the coupled ground has been developed. This model is useful for simulating the vibrator baseplate dynamics, evaluating the impact of the baseplate on the coupled ground and vibrator baseplate design. Model data demonstrate that the vibrator baseplate and its stilt structure are subject to six significant resonant frequencies in the range of 10–80 Hz. Due to the low rigidity of the baseplate, the baseplate stilt structure experiences severe rocking motions at lower frequencies and the baseplate pad experiences severe flexing motions at higher frequencies. Flexing motions cause partial decoupling, which gives rise to increased levels of harmonic distortion and less useable signal energy. In general, the baseplate pad suffers more bending and flexing motions at high frequencies than low frequencies, leading less efficiency in transmitting the useable energy into the ground.     

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.     

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.     

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.     

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.     

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

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

THE AMPLITUDE AND PHASE RESPONSE OF A SEISMIC VIBRATOR*

The amplitude and phase response of a simple model is compared with the performance of a real vibrator working in the field. The field results show a characteristic phase response which confirms that the real drive force applied to the baseplate and its load impedance is faithfully represented by the acceleration of the reaction mass. It follows that all the parameters necessary to calculate the load impedance and the true power dissipated in the earth can be measured at the output of the vibrator. It also follows that the current method of baseplate phase compensation should be reconsidered.     

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.     

Source measurement effect on high fidelity vibratory seismic separation

High Fidelity Vibratory Seismic (HFVS) acquisition and separation can play an important role in today's land acquisition schemes. The method – in which multiple vibrators are swept simultaneously using sweeps with known phase encoding and then the data are inverted and separated into individual records – can improve productivity in the field and at the same time improve signal characteristics in the data. It relies on the measured weighted sum of accelerations (base plate and reaction mass) to invert the acquired data and separate the individual vibrator responses. Separation can be sub-optimal if the measured motions vary from the 'true source' input into the ground. Differences in true source and measured source can arise due to poor coupling between vibrators and ground, soil compaction or other factors. Using both a synthetic model and real data, we show that if the true source changes between sweeps but is not measured, vibrator responses can leak into adjacent vibrator responses upon separation. In a recent survey with HFVS acquisition, we observed a 25–30 dB separation between adjacent vibrators, which could be improved with greater reliability of the source measurement. The vibrator leakage can reduce the data quality considerably. We discuss the results of this survey and show that separation is affected by source measurement error. Further, we conclude that it is necessary either 1) to use source measurements that can capture the variability of the true source between sweeps or 2) to compensate for the source measurement variations in processing or in acquisition.     

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.     

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.     

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.     

Harmonic by harmonic removal technique for improving vibroseis data quality

In this paper, an improved method is presented to reduce vibrator harmonic distortion, one harmonic at a time and the method is illustrated with both simulated and field data. This method improves on the previous method that treated all the harmonics at once. The significant contribution in this procedure is a considerable reduction for the harmonics without any alteration for the weakest signals possibly present in positive and negative times. The core of the proposed technique depends on an accurate simulation for all the harmonics one by one existing in the positive and negative times of the data after cross‐correlation with the fundamental sweep and then subtracting the simulated harmonics from the original data using an optimization procedure. The steps and mathematical equations of the procedure are explained in detail in the body of the article in the section titled ‘harmonic by harmonic attenuation procedure’. Accordingly, a well‐developed procedure for enhancing the vibroseis data quality in both down‐ and up‐sweep data is illustrated. The procedure was tested on both synthetic and field data sets.     

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.     

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…     

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.     

震源车地震勘探技术在西安地裂缝勘探中的试验研究

为克服巨大的城市环境干扰,采集高质量的人工地震勘探数据,将大吨位震源车引入到西安地裂缝人工地震勘探中,并在真实城市道路环境下对已知地裂缝进行多项试验,包括抗干扰效果、勘探深度、震源工作参数、垂直叠加次数、不同主频检波器接收、勘探排列长度、最佳反射波接收窗口等,最后将人工地震勘探结果与钻探结果进行验证对比。结果表明,大吨位震源车在实际道路干扰环境下可以得到理想的地震反射剖面,剖面中地裂缝异常明显,推测解释的地裂缝在地面的位置可信。     

噪声水平上的爆炸余震信号提取

文中采用窄带滤波、傅里叶反变换及小波分析在噪声水平上提取弱余震信号。窄带滤波能够在背景噪声水平上以及（－10）－（20）dB记录中提取有用信息，但是往往产生截断误差，如果数字序列具有较大的初始幅值，所提取的波形在截断处发生畸变。傅氏反变换从负信噪比的记录中提取的信号质量较差，伴有高频干扰成分，对于信噪比为正数的记录，傅氏反变换提取的波形是良好的，该种方法同滤波方法一样，也存在截断误差。小波分析从信噪比为负数的记录中提取的信号不理想，而从信噪比为正数的记录中提取的信号质量最好，波形完整，首尾清晰，基线平直，无截断误差。窄带滤波，傅氏反变换及小波分析有两个共同点，基一，都能够从强背景噪声中提取有用信息；其二，信噪比越大，获取的信号质量越高。     

时变重力场球面模型反演算法和模拟实验

时变重力场是研究地球内部介质物性变化的重要手段。本文提出了一种适用于地面流动重力测量获得的时变重力信号的场源反演方法,该方法采用球坐标系下的六面体单元来模拟场源介质,适合大尺度地震流动重力测量数据的等效源模型构建。通过引入重力时变信号的一阶光滑先验条件,压制了时变重力信号中的短周期高频分量,可用于提取与地震孕育相关的长周期信号。通过理论和模型实验证明了本文算法的可靠性和稳定性,并使用南北地震带南段2014—2017年的流动重力实测数据进行了反演解释,获得了地壳内部等效场源的视密度时变信号,变化量级在正常地壳密度的±0.7‰之间,其空间形态受川滇菱形块体边界控制。研究成果可用于时变重力场模型解释和深部场源特征提取,可为地震重力前兆信号分析和相关研究提供完备的方法保障。