THE SUPPRESSION OF SURFACE MULTIPLES ON SEISMIC RECORDS*

A new approach is presented for the suppression of multiples reflected at the surface of a horizontally layered fluid or elastic medium, recorded at non-zero offsets from the source. The scheme used is to extract the effect of the free surface in the frequency-wavenumber domain and then to replace this surface by a non-reflecting boundary. The multiple suppression operator requires a detailed knowledge of the source time function and the elastic properties of the medium between the source and the surface. For a stratified fluid or a liquid layer overlying a stratified elastic medium, complete multiple suppression can be achieved with noise free data. If only the vertical component is available for an elastic medium an approximate approach may be used which removes most of the multiple energy. Good results may be achieved with this multiple suppression scheme in the presence of noise. The method is designed to be used before records are stacked in a CDP gather.     

NON-NORMAL INCIDENCE STATE-SPACE MODEL AND LINE SOURCE REFLECTION SYNTHETIC SEISMOGRAM*

This paper is directed at modeling layered media. We extend the plane-wave normal-incidence state-space model developed by Mendel, Nahi and Chan in 1979, to the non-normal incidence case. To do this we introduce a shifting principle, a zero-offset wavefront, and zero-offset travel times for different layers. We also develop an algorithm for obtaining a synthetic line source reflection seismogram. In this algorithm non-normal incidence plane-wave seismograms are summed over a range of incident angles. The algorithm is based on a modified version of Sommerfield's (1896) theorem. Simulations of acoustic and elastic media are included which illustrate the applicability of our plane-wave and line source seismograms for both elastic and acoustic cases.     

NON-NORMAL INCIDENCE INVERSION: EXISTENCE OF SOLUTION*

An inverse problem is one in which the parameters of a model are determined from measured seismic data. Important to the solution of inverse problems is the issue of whether or not a solution exists. In this paper we show, in a constructive manner, that a solution does exist to the specific inverse problem of determining the parameters of a horizontally stratified, lossless, isotropic and homogeneous layered system that is excited by a non-normal incidence (NNI) plane wave. Mode conversion between P- and S-waves is included. We develop a seven-step layer-recursive procedure for determining all of the parameters for layer j. These parameters are P-wave and S-wave velocities and angles of incidence, density, thickness, traveltimes, and reflection- and transmission-coefficient matrices. Downward continuation of data from the top of one layer to the top of the next lower layer is an important step in our procedure, just as it is in normal incidence (NI) inversion. We show that, in order to compute all parameters of layer j, we need to (and can) compute some parameters for layer j+ 1. This is a non-causal phenomenon that seems to be necessary in NNI inversion but is not present in NI inversion.     

LONG PERIOD SEA-FLOOR MULTIPLES AND THEIR SUPPRESSION*

Multiple sea-floor reflections in deep water often are not effectively suppressed by either CDP stacking nor standard predictive deconvolution methods. These methods fail because the reflection coefficient varies markedly with angle of incidence and also because of the variation of arrival time with offset and because of dip. For a reasonablly flat sea-floor, multiples of various orders and the primary sea-floor reflection which have all been reflected at nearly the same angle lie along a straight line through the origin in time-offset space. This line is called the “radial direction.” The multiples which lie along this line show a systematic relationship because they all experience the same water-bottom reflection effect. In other words, multiples behave in a stationary manner along the radial directions on multi-trace seismic records. A technique of multi-channel predictive deconvolution, called “Radial Multiple Suppression,” utilizes this aspect to design Wiener operators for the prediciton and suppression of water bottom multiples. The effectiveness of the technique is demonstrated by the study of field records, autocorrelations, velocity analyses, and stacked sections before and after Radial Multiple Suppression processing.     

TIME-DOMAIN COMPUTATION OF NON-NORMAL INCIDENCE WAVEFIELDS IN PLANE LAYERED MEDIA1

A new time-domain method is introduced for the calculation of theoretical seismograms which include frequency dependent effects like absorption. To incorporate these effects the reflection and transmission coefficients become convolutionary operators. The method is based on the communication theory approach and is applicable to non-normal incidence plane waves in flat layered elastic media. Wave propagation is simulated by tracking the wave amplitudes through a storage vector inside the computer memory representing a Goupillaud earth model discretized by equal vertical transit times. Arbitrary numbers of sources and receivers can be placed at arbitrary depth positions, while the computational effort is independent of that number. Therefore, the computation of a whole plane-wave vertical seismic profile is possible with no extra effort compared to the computation of the surface seismogram. The new method can be used as an aid to the interpretation of plane-wave decomposed reflection data where the whole synthetic vertical seismic profile readily gives the interpreter the correct depth position of reflection events.     

COHERENCY WEIGHTING—AN EFFECTIVE APPROACH TO THE SUPPRESSION OF LONG LEG MULTIPLES*

Long leg multiples can be suppressed by a method which provides an alternative to weighted common-depth-point stacking and multichannel stacking filtering. The suppression is achieved by coherency weighting whereby the time-dependent weighting factor decreases as the semblance of the multiple reflections increases. The algorithm of the method is described. Its efficiency is discussed in relation to the input data and results of its application to marine seismic data are presented. For practical application, the stacking velocity of the multiples has to be known. As the process is based on stacking velocities, different types of multiples can be handled, for instance water-bottom multiples or internal multiples. The parameter analysis shows that the degree of multiple suppression can easily be controlled by adapting the parameters of the procedure to the field conditions. During the suppression of multiples, the primaries are saved according to the moveout differences between the two. The non-linear behaviour of the process causes signal suppression and distortion effects, which have to be corrected by AGC normalization and low-pass filtering. Among the various applications available, only the suppression of long leg water-bottom multiples is treated here. The results show that their suppression on the basis of moveout differences is efficient even when standard length streamers are used in regions with water depth of up to 1500 m and more, if the stacking velocity of the primaries is about 10 to 20% higher than that of the multiples. Even if those parts of the primaries which are masked by the multiples are suppressed in the individual common-depth-point gathers by the procedure, the remaining primaries in the AGV stacked section are largely uncovered by the multiple suppression.     

A LAYER-STRIPPING TECHNIQUE FOR THE SUPPRESSION OF WATER-BOTTOM MULTIPLE REFLECTIONS*

A new method to suppress water-bottom multiples (water-bottom reverberations) uses the fact that in the domain of intercept time and ray parameter (τ–p domain) the water-bottom reverberations are strictly periodical for a horizontal flat sea bottom. Using this property a comb filter can be designed. The window of the filter should be approximately equal to the duration of a source pulse. The algorithm finds the maximum of the periodical energy throughout the τ–p domain and then designs the comb filter which eliminates the water bottom reverberations from each trace in the τ– p domain. This process can be repeated for higher order reverberations. Finally the τ–p domain with attenuated multiples is transformed back to the conventional x -- t space. The method is illustrated on a variety of synthetic data and on a set of real marine CMP data acquired in the North Sea near the Norwegian shore.     

A FORWARD GENERATED SYNTHETIC SEISMOGRAM FOR EQUAL-DELAY LAYERED MEDIA MODELS*

A forward solution for the reflection response of a parallel stratified lossless medium characterized by discrete reflection coefficients and unequal layer delays, for a normally incident pressure source signal, is presented. The notation, which details the reflection history of each wavelet in a response record, facilitates systematic enumeration of all terms in the reflection impulse response model, the determination of compact closed form expressions for amplitudes and delays of multiply reflected wavelets, and the aggregation of dynamic analog groups. An equal delay time constraint on layer thicknesses leads then to the reflection sequence or synthetic seismogram structure as an infinite sum of wavelets by their order of reflection.     

时距平面剩余时差选择加权叠加压制多次波方法及效果

一个CDP道集经动校正后,多次波仍存在剩余时差δt,它随炮检距呈抛物线变化。由δt的变化即可求出多次波剩余波数△K的分布规律。据此,可设计一个理想的叠加剩余波数响应,用于动校正后的CDP道集的加权叠加。在t～x平面选择δt,即可改变叠加剩余波数响应,从而改变各道的加权系数,实现最佳压制多次波的效果。     

LEAST-SQUARES FILTERS WITHOUT TRANSIENT ERRORS: AN EXAMINATION OF THE ERRORS IN LEAST-SQUARES FILTER DESIGN*

A new approach has been developed for the design of cross-equalization filters by the least-squares method. The filters estimated by this new exact method are subject to only two types of error: bias and random error. Cross-equalization filters estimated by a more conventional least-squares method are further subject to “transient error”. This type of error becomes important when designing filters from a data gate of a length comparable with the length of the filter, i.e., less than four times the length of the filter. The effect of altering various design parameters has been investigated for the new method. It has been found that the proportion of bias in the filter decreases as the effective filter length increases, whereas the random error in the filter decreases with increase in either the signal-to-noise ratio of the data or the ratio of the data duration to the filter length. The level of whitening applied to the auto-correlation matrix before inversion was not found to be a critical design parameter. Also, two techniques have been tested for reducing any anomalous d.c. component in the calculated filter.     

ODD-DEPTH STRUCTURES FOR DETERMINISTIC DECONVOLUTION AND SEISMOGRAM TESTING*

In odd-depth structure the two-way traveltime to each boundary is constrained to be an odd integer. The odd-depth property of a model is exposed to possible refutation under a seismogram test. Test function is a simple transformation of a synthetic seismogram. For an odd-depth model the test function has identically the value 1. The testability of a synthetic seismogram over an odd-depth structure provides a method of deterministic deconvolution. There is no need of specialized assumptions, like the minimum-phase property, about the source wavelet. The deconvolution may be performed in the absence of the early segment of a seismogram.     

TWO DIMENSIONAL STRAKHOV'S FILTER FOR EXTRACTION OF POTENTIAL FIELD SIGNAL*

A two dimensional linear filter operator based on Strakhov's method for the extraction of a potential field signal (gravity or magnetic) from a background of homogeneous random noise has been developed. An algorithm to solve for the coefficients of the operator is given together with the corresponding computer program in FORTRAN language. It may be used on any machine using FORTRAN IV. A few examples of the operator are also included.     

COMBINED SWEEP SIGNALS FOR CORRELATION NOISE SUPPRESSION*

Comparison of both synthetic and field data shows that considerable suppression of correlation noise can be achieved with the Combisweep technique and with the Encoded Sweep Technique. In the first technique, the spectrum is shaped by superposition of linear sweeps with different frequency range; in the second technique, short sweeps of different polarity are combined to form the “alphabet” of a code.     

ω-k FILTER DESIGN*

Two-dimensional band-pass filters can be constructed by a simple extension of the theory of one-dimensional band-pass filters. Similarly to the one-dimensional analogue the shape of the two-dimensional filter is important in determining its effectiveness. The band-pass filter formulation can be further refined so that the filter will concentrate its rejection energies in certain areas of the ω, k plane. Such band-pass, band-reject filters are found by solving a set of simultaneous equations.     

INVERSE FILTER COEFFICIENTS FOR THE COMPUTATION OF APPARENT RESISTIVITY STANDARD CURVES FOR A HORIZONTALLY STRATIFIED EARTH*

In this paper a fast method is developed for computing apparent resistivity curves for known layer configurations. The method is based on the application of a linear filter to determine the apparent resistivity curve from, the kernel function.     

ω-k FILTER DESIGN*

Two-dimensional band-pass filters can be constructed by a simple extension of the theory of one-dimensional band-pass filters. Similarly to the one-dimensional analogue the shape of the two-dimensional filter is important in determining its effectiveness. The band-pass filter formulation can be further refined so that the filter will concentrate its rejection energies in certain areas of the ω, k plane. Such band-pass, band-reject filters are found by solving a set of simultaneous equations.     

低纬度磁异常化极方法——压制因子法

磁异常化极转换是磁异常解释的重要基础，针对频率域磁异常化极现状，本文对低纬度化极的困难原因进行了深入分析，并提出了针对性的化极措施——压制因子法. 与已有研究比较，该方法具有简单直观、控制参数少因而便于实际操作的特点，模型计算表明了方法的有效性. 在我国南海低磁纬度地区实测磁异常的计算解释中，应用该化极方法，取得了良好的效果.     

单台测震分析与资料交换程序

按照地震台站测震分析处理过程,我们在长城0520DH微机上建立了一套地震分析与资料交换程序。该程序具有分析精度高、操作简便、易于掌握等特点,可供Ⅰ、Ⅱ类地震台站使用。     

DIGITAL LINEAR FILTER FOR COMPUTING TYPE CURVES FOR THE TWO-ELECTRODE SYSTEM OF RESISTIVITY SOUNDING*

In this paper a technique for computing type curves for the two-electrode resistivity soundings is presented. It is shown that the apparent resistivity due to the system can be represented by a convolution integral. Thus, it is possible to apply the principle of digital linear filtering and compute the desired type-curves. The filter function required for the purpose is found to be identical with that used to compute the EM sounding curves for the two coplanar horizontal loop system. It is further shown that from the two-electrode apparent-resistivity expression one can easily derive the apparent resistivities for other configurations. A comparison of depths of investigation for various systems reveals that the two-electrode system has greater depth of investigation than other conventional systems. This is also supported by the field example presented in the end, which illustrates the relative performance of the two-electrode system vis-a-vis the Wenner system.     

THE KALMAN FILTER AS A PREDICTION ERROR FILTER*

In mathematical statistical filtering the deconvolution problem can be solved by two different methods:

• 1 by inverse filtering
• 2 by calculating the prediction error.

Both methods are well known in the theory of Wiener filters. If, however, the generating process of the signal is known and can be described by a set of linear first order differential equations, then the Kalman filter can also be used to solve the deconvolution problem. In the case of the inverse filtering method this was shown by Bayless and Brigham (1970). But, while their method can only be used if the original signal is a colored random process, this paper shows that in the case of a white process the prediction error filtering method is a more appropriate approach. The method is extremely efficient and simple. This can be demonstrated by an example which maybe of special interest for seismic exploration.