Methods of adaptively reducing noise transients in co-sensor seismic signals

Abstract: Seismic data from co-located sensors of different genera are formatted into common velocity-receiver gathers and common pressure receiver gathers. A first ratio between the velocity signature amplitudes and the pressure signature amplitudes is measured within a fixed analysis window (14). A second ratio is measured between the pressure and velocity signatures for a weighting-zone window (18) of limited extent. The first and second ratios are combined to form an equalization operator. The pressure and velocity signatures from within the weighting zone window (18) are combined in the presence of the equalization operator to define a transient-free time-scale datum.

Description: This invention relates to methods of processing seismic data that has been contaminated by noise transients such as environmental noise generated by shrimp, croaker fish, whales and the like, and is applicable particularly to shallow-water environments. During the course of a marine seismic survey, a long string of several hundred seismic sensors is towed through the water by means of a suitable cable. Periodically, such as every six to eight seconds, an acoustic source is triggered to emit a wavefield to insonify the subsurface earth layers thereby to generate a reflected acoustic wavefield that is detected by the sensors. The detected mechanical signals representative of the reflected wavefield are converted to electrical signals by the sensors. The electrical signals are transmitted through conductors in the cable to the towing ship where the signals are recorded for future analysis to display the topography of the subsurface earth layers. The seismic signals from each source-triggering are digitized and recorded in the form of a set including a plurality of discrete time series, one time series corresponding to each sensor or sensor group connected to the cable. The set of discrete time series resulting from each source-triggering may be displayed as a single multi-trace visual seismic record. The traces may take the form of variable amplitude or variable density or variably colored traces. In shallow water at depths of the order of 6-76m (20-250 feet), so-called bay cables are used. For certain technical reasons, a typical bay cable includes collocated sensor pairs of different genera. Usually a pressure sensor is paired with a gimbal-mounted particle-velocity sensor, both of which sensors are customarily laid directly on the sea floor. Shallow water operations are so well known and routine to those skilled in the art of geophysical exploration that prolix explanation and accompanying illustrations should not be needed. Any unwanted wavefield that appears on a seismic record is considered to be noise. Many different sources create the noises that afflict the seismic signals received during the course of a shallow-water marine survey. Whereas whale-talk is bothersome noise to the geophysicist, it may be music to the ears of an environmentalist. Reverberation due to acoustic signals bouncing between the sea floor and the water surface seriously distorts the desired reflected wavefields. Various methods are known for eliminating such contamination such as US Patent Specification No. 5,365,492. That specification describes a method of scaling the outputs from sensors of different genera relative to each other. Additionally, an adaptive noise cancellation process was taught for reducing the noise content of the respective signals prior to summation. Seismic shots generated by a rival seismic crew produces coherent noise. US Patent Specification No. 4,937,794 explains how coherent noise on seismic records may be suppressed by formatting the common shot gathers to common receiver gathers. Pairs of seismic traces from a common receiver gather are corrected for differential normal moveout, weighted in inverse proportion to the RMS signal power, and combined to generate compressed common receiver gathers. The compressed common receiver gathers may be reformatted as common midpoint gathers for further processing. Noise due to near-surface scatterers may be removed or minimized in accordance with the teachings of US Patent Specification No. 5,293,352. Therein it is taught that a set of raw common shot gathers are resorted to common receiver gathers. The wavefield envelopes from the common receiver gathers are migrated using one-half the near-surface velocity to provide migrated data sets. The migrated data sets are resorted back into common shot gathers and subtracted from the original raw common shot gathers to provide coherent-noise reduced data sets. A method that was designed for use with vibroseis data on land explains that a reference model of the level of a valid seismic signal is built for each of a number of time windows during a seismic data recording cycle. The model is selected by obtaining the average of the absolute magnitudes for each time window from a number of sweeps to form a set of averages. The median of each set is selected and is padded by a suitable coefficient to provide the reference model for each time window. Subsequently a recording is made. Each data sample from the recording is compared to the reference model corresponding to the time window that includes the sample. If the amplitude level of the data sample exceeds the reference model level, the sample is suppressed prior to summing; see US Patent Specification No. 4,344,158. Such a method can be adapted for use in shallow marine projects. A particularly troublesome type of noise, that is not easily handled by known methods, is biological. Figure 1 of the accompanying drawings illustrates an example of such noise which is believed to have been due to a Grand Betail or to croaker or similar fish. The noise bursts are the inverted arrow-like diffraction patterns such as 10. The desired wavefields are the broad gently-curved hyperbolic waveforms generally indicated by 12. There is a need for a noise-abatement technique for attacking unwanted singular transients. Various aspects of the invention are exemplified by attached claims. According to another aspect of the invention, there is provided a method of adaptively suppressing undesired singular transients from time scale recordings of seismic signals that are detected concurrently by a plurality of velocity-responsive and pressure-responsive co-located co-sensors; a reflected wavefield is generated in a moisture-laden environment, the wavefield being detected by the respective co-sensors to provide velocity and pressure time-scale signatures; within a fixed analysis window a scale factor, k, is iteratively defined for scaling the velocity signatures to the amplitude of the pressure signatures; the signatures from the pressure-responsive sensors are grouped to form a first common pressure-receiver gather; the detected velocity-responsive signatures are grouped to form a second, common velocity-receiver, gather; a ratio R 1 is defined between the second and the first gathers, the respective common receiver gathers are partitioned into a plurality of weighting zones; within a first weighting zone, a ratio R 2 is defined between the pressure signatures and the velocity signatures within that zone; the ratios R 1 and R 2 are combined to form an equalization operator; a substantially transient-free time-scale datum is defined by combining the data in the common pressure-receiver weighting zone with the scaled data in the common velocity-receiver weighting zone in the presence of the equalization operator. For a better understanding of the invention, and to show how the same may be carried into effect, reference will now be made, by way of example, to the accompanying drawings, in which:- FIGURE 1 is an example of biological noise; FIGURE 2 is a seismic recording of a raw common hydrophone-receiver gather; FIGURE 3 is a seismic recording of a raw common velocity-receiver gather; FIGURE 4 is the seismic recording after processing the recordings of FIGURES 2 and 3 in accordance with the preferred embodiment. The preferred embodiment of the invention is intended to be a method for attenuating noise in seismic data, particularly but not limited to biological noise, the method being efficient in terms of computer processing requirements. Figure 2 is a display of a common pressure-receiver (hydrophone) gather, H, while Figure 3 is a display of a common velocity-receiver (geophone) gather, G. In the Figures, the vertical axis is two-way reflection time in seconds and the horizontal axis is the lateral sensor separation which was 25 meters between traces. Throughout this text, the term "velocity receiver" or "velocity-responsive receiver" will be used interchangeably with "geophone"; also "pressure receiver" or "pressure-responsive receiver" will be used interchangeably with "hydrophone". In the context of this disclosure, the term "signature" means the aspect of a waveshape, expressed in terms of a magnitude as a function of time, that makes it distinctive. Using an iterative summation/auto-correlation technique taught by US Patent Specification No. 5365492, earlier cited and which is incorporated herein by reference, determine the scale factor, k, for scaling the average geophone signal amplitude to the hydrophone signal amplitude. In the absence of other noise abatement processing, the normalized summation of the gathers is Σ(τ) = ½ (H + k*G). The method for determining the scale factor, k, may be summarized as follows: The hydrophone pressure signature is filtered to create a best estimate of a geophone velocity signature. The estimated geophone velocity signature is subtracted from the real geophone velocity signature plus embedded noise to yield a nearly pure noise signature. Now, the noise signature is subtracted from the real geophone velocity signature to provide a clean geophone velocity signature. A scale factor is iteratively applied to the clean geophone velocity signature, summed with the hydrophone pressure signature and auto-correlated. A measure of the scale factor is suitably determined from a study of the main lobe and the first side lobes of the autocorrelation trace. The iteration progresses until a suitability measure is maximized. Within a fixed analysis window such as 14 from the hydrophone traces of Figure 2, the detected pressure signatures are grouped together to form a first, common hydrophone gather. Within a fixed analysis window 16 from the velocity-receiver traces of Figure 3, the geophone signatures are grouped together to form a second, common geophone gather. The measured average trace amplitudes of the geophone signals relative to the hydrophone signals in the fixed analysis window are used to define the ratio R 1 . The fixed analysis window is established over a selected lateral and temporal region, preferably from a portion of the recordings that is common to both gathers but a portion which is representative of the signal level of the gathers taken as whole. The dimensions of the analysis windows is arbitrary. In the case of Figures 2 and 3, the analysis window extends from about 2.0 to about 3.0 seconds vertically and embraces about fifty traces laterally for each gather. For the next step, the two gathers are partitioned into weighting-zone windows, such as the one designated 18 and 20 in Figures 2 and 3, and which are common to both gathers. The weighting-zone windows have restricted dimensions and may be one to five traces wide and of some arbitrary length such as 0.25, 0.50 or 1.5 seconds. The dimensions depend upon the desired noise-abatement resolution. Within each weighting-zone window, compute the hydrophone-to-geophone amplitude ratio R 2i from the average amplitudes of the pressure signals H i and the velocity signals G i within the weighting-zone windows for each gather (i is a running scan index as explained below). In operation, the weighting-zone window is stepped through the respective gathers both spatially (laterally) and temporally (vertically) until the entirety of each gather has been scanned. Thus there is a single value for R 1 but an array of values for R 2i , a new value being computed for each scan, i , where i=l,...,n and n is an arbitrary number of scans. Thus, for the first scan, the hydrophone to geophone ratio would be R 2,1 and the equalization operator for that scan is m l . An equalization operator, m i for each scan, is computed from m i = R 1 * R 2i A noise-free time-scale datum f(τ) is obtained by combining the signatures in the common hydrophone receiver weighting zone with the scaled signatures in the common velocity receiver weighting zone in the presence of the equalization operator as follows: f(τ i ) = (H i + k*G i *m i ) ÷ (1 + m i ). Figure 4 is the result of the application of (2) to the summation of the data from Figures 2 and 3. The noise bursts 22 seen on Figure 2 and particularly, the noisy traces such as 24 and 26 on Figure 3 have been substantially reduced in Figure 4. In some applications, the equalization operator may be raised to a power such as m i 2 to accommodate the method to respond to rapid changes in signal amplitude. The best mode of operation may be demonstrated by hypothetical numerical examples as follows: Assume by way of a first example that the scale factor, k, as determined by the method of US Patent Specification No. 5365492 is 1.5. For the fixed analysis window, it was found that R 1 = 2.0. In a selected weighting zone window, i = x, R 2x was found to be 0.5. Entering those quantities into (2) above, f(τ x ) turns out to be 0.5(H + 1.5G x ), the same as a normalized, noise-free summation Σ(τ). Thus, for scan i = x, neither the hydrophone gather nor the geophone gather saw a noise burst within the xth scan region. For scan i = y, as a second example, R 1 = 2.0 as before and R 2y = 5.0, that is, the hydrophone gather saw a noise transient having an amplitude of ten times over that for scan x. For this case, f(τ y ) = (H y /ll) + 1.5*(10/ll)G y . The noise transient in the hydrophone weighting zone for the yth scan contributes only about 0.1 times as much to the sum as does the geophone trace. In a final example, for scan i = z, R 1 is 2.0 but R 2z = 0.05. The geophone trace in that weighting zone received a noise transient whose amplitude is ten times greater than normal. For the case of scan i = z, f(τ 2 ) = (H 2 /l.l) + 1.5*(.1/l.l)G 2 . Thus, the effect of the noise transient on the geophone trace in the zth geophone weighting zone is minimized. These teachings herein preferably may be implemented by use of any commercial general purpose computer of a type ordinarily used for processing of seismic data, the computer including an adequate memory and an arithmetic processor which may be resident in the CPU in the form of a subroutine or it may be a peripheral such as an array processor. This invention has been described with a certain degree of specificity by way of example but not by way of limitation.

Claims: A noise abatement method for attenuating singular transients in seismic signals, comprising: receiving pressure-responsive seismic signals and velocity-responsive seismic signals and formatting the received seismic signals into a common pressure-signal gather, H, and a common velocity-signal gather, G, respectively; establishing a fixed analysis window over a region having selected spatial and temporal dimensions common to both said gathers, and within said window, determining the average amplitude ratio R 1 between the velocity-signal gather G and the pressure signal gather H; partitioning the respective gathers into a plurality of common weighting zones having restricted lateral and temporal limits; within each weighting zone determining the ratio R 2i between the average amplitude of the pressure-signal gather H i and the velocity-signal gather G i , i being a scan index; combining R 1 and R 2i to form an equalization operator m i ; defining a substantially transient-free time-scale datum f(τ i ) within each weighting zone window by an additive combination of the pressure-signal H i gather and a scaled velocity-signal gather k*G i in the presence of the equalization operator m i . A method of adaptively reducing undesired singular transients from time-scale recordings of seismic signals in a moisture-laden environment and detected concurrently by a plurality of co-located velocity-responsive and pressure-responsive co-sensors, the method comprising: a) generating pressure-responsive and velocity-responsive time-scale signatures; b) defining a scale factor, k, for scaling the velocity-responsive signatures to conform to the pressure-responsive signatures in amplitude; c) grouping the pressure-responsive signatures to form a common pressure-responsive hydrophone gather, and grouping the velocity-responsive signatures to form a common velocity-responsive gather; and d) combining signatures in the pressure-responsive gather with scaled signatures in the velocity-responsive gather; e) defining a fixed analysis window common to both gathers; f) defining a ratio R 1 between the common velocity-responsive gathers and the common pressure-responsive gathers within said fixed analysis window; g) partitioning the common pressure-responsive and the common velocity-responsive gathers into a plurality of common weighting zones; h) within a first weighting zone, defining a ratio R 2,1 between the pressure-responsive signatures and the velocity-responsive signatures; i) forming an equalization operator, m 1 , from a combination of R 1 and R 2,1 ; j) combining the pressure-responsive signatures in the common pressure-responsive weighting zone with the scaled signatures in the common velocity-responsive weighting zone in the presence of said equalization operator to define a substantially transient-free time-scale datum f(τ). A method according to claim 2, comprising repeating steps h-j for the remainder of the weighting zones. A method according to claim 1, 2 or 3, wherein said equalization operator is calculated from m i = R 1 * R 2i A method according to claim 1, 2 or 3, wherein said equalization operator is calculated from m i n = (R 1 * R 2i ) n , where n is an integer ≥ 1. A method according to any one of the preceding claims wherein the step of combining is implemented by the formulation f(τ i ) = (H i + k*G i *m i )/(l + m i ) where H i represents the common pressure-responsive gather data and k*G i represents the scaled common velocity-responsive data within a weighting zone window. A method of seismic exploration comprising generating a reflected seismic wavefield in a moisture-laden environment, utilising pressure-responsive and velocity-responsive sensors to generate time-scale recordings of seismic signals and employing a method according to any one of the preceding claims for reducing undesired singular transients in the time-scale recordings.

Inventors: CHAMBERS RONALD, DRAGOSETT WILLIAM H JNR, CHAMBERS, RONALD, DRAGOSETT, WILLIAM H., JNR.
Applicants: WESTERN ATLAS INT INC, WESTERN ATLAS INTERNATIONAL, INC.
Application year: 1995
Main classification: G01V1
IPC codes: G01V1