AU2002324926A1 - Method for attenuating noise in seismic data - Google Patents
Method for attenuating noise in seismic dataInfo
- Publication number
- AU2002324926A1 AU2002324926A1 AU2002324926A AU2002324926A AU2002324926A1 AU 2002324926 A1 AU2002324926 A1 AU 2002324926A1 AU 2002324926 A AU2002324926 A AU 2002324926A AU 2002324926 A AU2002324926 A AU 2002324926A AU 2002324926 A1 AU2002324926 A1 AU 2002324926A1
- Authority
- AU
- Australia
- Prior art keywords
- seismic data
- data traces
- traces
- noise
- geophone
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Description
METHOD FOR ATTENUATING NOISE IN SEISMIC DATA
This application claims the benefit of U.S. Provisional Application No. 60/322,504 filed September 14, 2001.
FIELD OF INVENTION
This invention generally relates to the field of seismic prospecting. More particularly, this invention relates to a method for attenuating noise from seismic data.
BACKGROUND OF THE INVENTION
Marine seismic exploration usually involves acquiring seismic data using a seismic acquisition system whose source initiates a down-going seismic wavefield. A portion of the down-going wavefield travels into the underlying earth where it illuminates subsea geologic formations. As it illuminates the interfaces or boundaries between the formations, part of the wavefield is returned (or reflected) back through the earth (propagating in the up-going direction). Part of the reflected wavefield is received by the seismic acquisition system, converted into electrical signals, and recorded for subsequent processing. An analysis of these recorded signals makes it possible to estimate the structure, position, and lithology of subsea geologic formations, thereby completing an important step in the exploration process.
Figure 1 shows a simplified example of a typical marine seismic acquisition system. A first ship 1 tows a seismic source 2 several feet below the surface 3 of the ocean. The seismic source 2 is activated to produce a down-going wavefield 4d that is at least partially reflected by a subsea formation boundary (i.e., an impedance discontinuity) 5. The up-going wavefield 4u then travels toward the sensors 6 and is detected.
The sensors 6 used in typical marine seismic exploration include pressure sensors and velocity (also referred to as "particle velocity") sensors. Typically, the pressure sensors are hydrophones and the velocity sensors are geophones. The hydrophones measure a scalar pressure and are not sensitive to the propagation direction of the wavefield. The geophones, which may be vertical geophones, provide
a vector response measurement whose polarity depends on whether the direction of propagation of the wavefield is up-going or down-going. The amplitude of the geophone response is also related to an angle of the propagation relative to the sensitive direction of the geophone. If a wavefield is recorded by a hydrophone and a vertical geophone having the same electronic impulse response, then a polarity comparison between the hydrophone and geophone measurement determines whether the wavefield is propagating in the up-going or down-going direction. Hydrophones and geophones disposed at the seafloor are typically used in pairs when collecting seismic data.
A pressure and vertical velocity combination from two-component (2C) or from multi-component (typically 4C) has been useful to cancel down-going multiples from a combined pressure and vertical velocity data signal. Noise that is present on one or the other of these sensors can limit the effectiveness of the combination technique. This invention attenuates noise that is normally present during data acquisition and has been found to be an important processing step before the combination.
In one type of marine seismic surveying, the sensors 6 are located at regular intervals in one or more ocean bottom cables (OBC) 7 that are arranged on the seafloor 9. When necessary, a second ship 8 is used to move the OBC 7 to a new position on the seafloor 9. Several miles of OBC 7 are typically deployed along the seafloor 9, and several OBCs are typically deployed in parallel arrangements. OBC 7 arrangements are particularly well suited for use in certain zones (i.e., zones cluttered with platforms or where the water is very shallow) where the use of ship-towed hydrophone arrays (not shown) (which are located proximate to the ocean surface 3 and are typically referred to as "streamers") is not practical. In another type of seismic surveying, the hydrophones and geophones are deployed on the ocean bottom as separate ocean bottom stations (OBS). A combination of separate OBS and OBC can be deployed. The geophone and hydrophone must be connected to a recording instrument typically on a vessel.
As shown in figure 2, in another type of marine seismic surveying, a marine cable or streamer 21 incorporating pressure hydrophones are designed for continuous towing through the water. A marine streamer 21 is typically made up of a plurality of active or live hydrophone arrays 23 separated by spacer or dead sections 25. Usually the streamers are nearly neutrally buoyant and depth controllers or depressors 27 are attached to depress the streamer 21 to the proper towing depth. A tail buoy with a radar reflector 29 is typically attached to the end of the streamer. The entire streamer may be 3-6 Km in length and is towed by a ship 31.
Seismic energy is transmitted by "body waves," which can be either compressional waves (P-waves) or shear waves (S-waves). P-waves are elastic solid waves in which particle motion is in the direction of propagation. S-waves are body waves in which the particle motion is perpendicular to the direction of propagation. Seismic energy can also be transmitted along boundaries between substances having different elastic properties by "surface waves."
Claerbout has disclosed that there is "crosstalk" between vertical and horizontal geophones due to non-vertical incidence of P and S reflected waves at the water bottom (Claerbout, J., The Leading Edge, v. 9, No. 4, pp. 38-40, April 1990). Crosstalk occurs when the vertical geophone records at least part of the S-wave energy. This energy becomes noise because the combination processing step assumes only P-waves. Other researchers have attributed this noise to non-perfect coupling and orientation of OBC geophones (Li, X.Y. and Yuan, J., Geophysical Prospecting v 47, No. 6, pp 995-1013, November 1999). It has also been suggested that this noise will be eliminated when the vector fidelity problem for OBC geophones is solved with the next generation of OBC hardware (Reid, F. and Macbeth, C. 70th Annual SEG. Int. Mtg. Expanded Abstr. Biogr. v. 1, pp 1213-1216, 2000). In addition, all the prior art methods attribute the "leakage noise" in the vertical geophone to shear waves, thus they do not consider other possible sources of the coherent noise (e.g., surface related waves).
In general, the industry recognizes that OBC vertical geophone data is typically noisy. Figure 3(a) shows seismic data traces collected using a hydrophone
on an ocean bottom cable. Figure 3(b) shows the corresponding geophone seismic traces collected using a vertical geophone on the same ocean bottom cable. As can be seen in Figures 3(a) and 3(b) there is significantly less noise from the seismic data traces collected with the hydrophone.
Practical attempts to suppress this noise usually include considering the acquisition system (i.e., vector fidelity, coupling, orientation) and then treating the noise as a random noise in the common-shot domain or common-shot order. There is lack of successful examples of noise suppression examples in real OBC seismic data. The lack of understanding on how to attenuate this noise has led to noise contaminated results of 2C OBC data combination worldwide. In some parts of the world (e.g., Gulf of Mexico), the noise is so severe that 2C OBC data is often considered to be useless. Accordingly, there is a need to attenuate, efficiently filter, or suppress at least part of this noise from OBC data. The present invention satisfies this need.
SUMMARY OF THE INVENTION
One aspect of the invention is a method for attenuating noise in seismic data traces, at least a portion of the noise is noncoherent in the common-shot domain but coherent in the common-receiver domain. The coherent noise has a move-out velocity different from the move-out velocity of the data signal in the seismic data traces. The method comprises (a) sorting the seismic data traces into common-receiver order; and (b) using the difference in move-out velocities to separate the coherent noise from the data signal.
In another aspect, the invention is a method for attenuating noise in two- component or multi-component ocean bottom cable (OBC) seismic data traces. The OBC seismic data traces comprises hydrophone data traces and vertical geophone data traces, at least a portion of the noise in the data traces being noncoherent in the common-shot domain but coherent in the common-receiver domain. The coherent noise has a move-out velocity different from the move-out velocity of the data signal in the seismic data traces. The method comprises: (a) sorting the data traces into
common-receiver order; (b) applying a move-out correction to the data traces using a move-out velocity chosen to overcorrect the data signal and undercorrect the coherent noises; (c) transforming the move-out corrected data traces from the x-t domain to a two-dimensional domain in which the move-out corrected data signal and the move- out corrected noise are separated; (d) removing at least a portion of the coherent noise from the transformed data traces; and (e) inverse transforming the data traces from the two-dimensional domain back to the x-t domain. Furthermore, an inverse moveout may be performed after the traces are transformed back to x-t domain with the velocity used in step (b).
BRIEF DESCRIPTION OF THE DRAWINGS
Figure 1 shows a prior art ocean bottom cable (OBC) marine seismic acquisition system.
Figure 2 shows a prior art ocean streamer cable marine seismic acquisition system.
Figure 3(a) shows seismic data traces collected using a hydrophone on a 2C OBC.
Figure 3(b) shows seismic data traces collected using a vertical geophone on a 2C OBC.
Figure 4 shows a flowchart of a first embodiment of the present invention for attenuating noise in seismic data traces.
Figure 5(a) shows the hydrophone seismic data traces from Figure 3(a) after processing using an embodiment of this invention.
Figure 5(b) shows the geophone seismic data traces from Figure 3(b) after processing using an embodiment of this invention.
Figure 6 shows a flowchart of a second embodiment of the present invention for attenuating noise in seismic data traces.
Figure 7(a) through 7(e) schematically illustrate a preferred embodiment of the present invention.
Figure 8(a) shows seismic data traces in common-shot order collected using a hydrophone on a 2C OBC.
Figure 8(b) shows seismic data traces in common-shot order collected using a vertical geophone on a 2C OBC.
Figure 9(a) shows the seismic data traces from figure 8(a) in common-receiver order.
Figure 9(b) shows the seismic data traces from figure 8(b) in common-receiver order.
Figure 10 shows geophone seismic data traces from figure 9(b) in common- receiver order after processing using the steps outlined in the flowchart of figure 6.
DETAILED DESCRIPTION
We have discovered that at least part of the geophone noise is noncoherent in common-shot order but coherent in common-receiver order. This coherency represents a major property allowing separation of the noise and the signal. Coherent properties of the noise in the common-receiver order has not been noticed by some major seismic contractors. Those who recognize coherency of the noise have suggested a number of theories explaining the noise origin. As discussed previously, most of the theories define the observed noise as being related to shear wave energy recorded on vertical geophone.
The embodiments of the invention relate to methods for attenuating noise in seismic data traces, at least a portion of the noise being noncoherent in the common- shot domain but coherent in the common-receiver domain. The coherent noise has a move-out velocity different from the move-out velocity of the data signal in the seismic data traces. As shown in Figure 4, the method comprises two basic steps of sorting the seismic noise into common-receiver order 32 and using the difference in move-out velocities to separate coherent noise from the data signal 33.
In seismic exploration, the terms "signal," "data trace," "seismic data trace," "trace," "shot record," "recorded data," and "measured data" refer to data recorded by a sensor pair after the seismic source has been used to produce a wavefield. A "single
shot record" refers to data recorded after a single source excitation (i.e., a single firing of an airgun array). The shot record typically comprises data recorded at each of a plurality of subsea locations with OBC sensors. For example, after the source is fired the wavefield is sampled at a selected sampling rate for a predetermined time interval. The recorded data forms a data signal for each sensor location. A summation of a plurality of signals having common subsurface reflection points may be used to generate a "stack" of data. A stack may then be interpreted to identify and classify "events." Events may comprise, for example, direct arrivals that travel directly from the source to the sensors, primaries that have reflected off of subsea formations (i.e., as a subsea salt structure, or a subsea hydrocarbon-bearing formation), or reflections off of an ocean-air boundary (e.g., ghosts giving rise to multiples).
Sorting the seismic data traces into common-receiver order
We have discovered that the noise on the OBC vertical geophone data represents the shear waves converted from primary wave sources (PS-Waves). The noise exhibits noncoherent or random properties in the original acquisition shot order or common-shot order but exhibits coherent properties in the common-receiver order. In common-shot order all the traces recorded for a single shot are plotted in a two- dimensional fashion with a trace's location from the shot location being the shot-to- geophone distance (for all the receivers) and with time being the other dimension. Whereas, in common-receiver order all the traces recorded at one particular receiver location are typically organized in a two-dimensional fashion with a trace's location from the receiver location being shot-to-geophone distance (for all shots recorded) and time being the other dimension. Therefore, the original common-shot order needs to be sorted into a common-receiver order.
Sorting from original common-shot order to common-receiver order involves separating the data based on the receiver location. Once separated or sorted into common-receiver data at least a portion of the noise and signal are coherent in the common-receiver order.
Separating the coherent noise from the data signals
As stated above, the noise on the geophone data represents shear waves converted from primary wave sources (e.g., PS-waves). P to S converted waves and reflected monotype (PP-waves) waves have substantially different move-out velocities. This fact allows use of filters (e.g., multi-channel filters) to separate the converted shear waves (noise) and reflected PP-waves (the signal component or signal). A filter is part of a system that discriminates against some of the information entering it. The discrimination is usually on the basis of frequency, although other bases (i.e., wavelength, moveout velocity, coherence, or amplitude) may be used. Either signal or noise may be filtered out. While filtering is well known in the art the discovery that geophone data noise is coherent in the common-receiver order but not in the common-shot order permits the noise to be filtered. Therefore, the noise cannot be effectively filtered in the common-shot order.
A preferred embodiment of the present invention is to correct the seismic data after the traces are sorted into common-receiver order by performing a normal- moveout correction and then separating the coherent noise from the data signal. A normal moveout is the variation of reflection arrival time because of the source-to- geophone distance (offset). Therefore, the normal-moveout correction is a time correction applied at an arrival at some distance between the source-to-receiver location that makes the arrival appear to have been recorded at a zero distance source- to-receiver location.
As stated above the noise is characterized by substantially lower moveout velocity than the signal. This fact enables use of multidimensional filters for their separation. A moveout correction before multichannel filtering is not a required step but it may significantly simplify the filter design.
We have discovered that the velocity for the hyperbolic moveout approximation of the PS-waves rarely exceeds 1500 meters/second (m/sec). However, primary reflected PP-waves in the OBC environment have moveout velocities higher than 1500 m/sec. The velocity of sound in water is defined at approximately 1500 m/sec and does not significantly change due to the usual variations in water depth and sediment velocities. Therefore, the hyperbolic moveout
correction of the velocity of sound in water represents a general parameter with worldwide applicability. After performing a normal-moveout correction of common- receiver gathers using a velocity of 1500 m sec the resulting PP events will be overcorrected and the noise (PS events) will be undercorrected with a relatively flat line. The characteristics of PP- and PS-waves permit separating the signal and the noise using data independent multi-channel filters. While a preferred embodiment is to use the velocity of sound in water any velocity can be chosen so long as the velocity will overcorrect the signal and undercorrect at least a portion of the noise.
Another preferred embodiment in separating the noise from the signal is to perform a multi-channel transformation of the common-receiver data from the original x-t domain. A transformation is the conversion from one representation to another. For this data we take data traces in the x-t two-dimensional space or where x is the source-geophone distance (offset) and t is the arrival time or two-way traveltimes and map the data into a different two-dimensional space in which the data and the noise are separated. Different transformations may be used by using different filters (i.e., f- k, tau-p, radon etc). One two-dimensional transform is the Fourier transform, which converts a time-function seismic data trace into a frequency-domain representation. The tau-p transformation is a conversion of x to p (where p is the reciprocal of the apparent velocity and is also called apparent slowness) and tau is the intercept time at x = 0. The tau-p transform is also called tau-p mapping and a radon transform. Those skilled in the art will recognize other domains suitable for this technique.
The preferred approach is to use the Fourier or f-k transform. After the f-k transform, energy from the undercorrected noise is localized in the half plane of negative spatial frequency (negative k-parameter). Whereas, overcorrected signal energy is localized in the positive k half plane. However, if the data has not been corrected (e.g., using a normal moveout of the velocity of sound in water) a velocity separation line will need to be defined that separates the noise from the signal. Once the velocity line is defined, the noise or signal can be removed by removing data above or below the line.
The signal or the data needs to be separated or filtered. The type of filtering utilized will depend on the type of transformation and how the data was treated. In the f-k domain, the overcorrected signal is localized in the upper half plane of the f-k plane, which corresponds to positive k values. Alternatively, the under-corrected noise is located in the lower half plane, which corresponds to negative k-values. Those skilled in the art will recognize that the k sign depends on the sign convention used for a particular f-k transform implementation and therefore, it may be reversed. In this example, the noise removal may be achieved by zeroing or reducing the lower half plane of the f-k plane.
An inverse multi-channel transformation may be performed if the data was previously transformed and there is a need for the data to be in conventional x-t two- dimensional space. This step should be consistent with the previous transformation performed before the filtering. For example, if the forward f-k transform was previously used then an inverse f-k transform should be used and if a forward tau-p transform was used then an inverse tau-p transform should be performed. If a normal moveout was performed and the data needs to be in the original format, an inverse normal moveout can be performed to keep the seismic traces in original format with the noise suppressed. Furthermore, the data may be left in common-receiver domain or sorted back to common-shot domain. Figures 5(a) and 5(b) respectively show the hydrophone and geophone seismic data traces from Figures 3(a) and 3(b) after the noise is separated from the data. The geophone data as shown in Figure 5(b) illustrates that the noise is significantly attenuated when compared to the original data as shown in Figure 3(b).
While most of the observed noise is in the geophone data, this technique can be applied to hydrophone data. Therefore, this method works for geophone and hydrophone data. Preferably, if the geophone and hydrophone data are to be summed, during subsequent processing, the same technique would be applied to both the geophone and hydrophone data.
An improved summation process is disclosed in co-pending US application 60/293,716. The summation process requires selection of either the pressure trace or
- l i ¬
the velocity trace for processing. Polarity differences between recorded pressure and velocity measurements are used to identify and separate primary and source side peg- leg multiples from down-going wavefields. After matching the recorded signals, processing techniques of the invention, based on a mathematical reformulation of the ocean bottom dual sensor problem, are applied to the matched pressure and velocity data. The down-going wavefield is estimated by subtracting the velocity from the pressure and dividing by two. The down-going wavefield is multiplied by a factor that includes an estimate R* (which may be an estimate for the water bottom reflection coefficient, R) and subtracted (or added) to the pressure (or velocity). In this way the down-going wavefield can be used to eliminate the trapped water bottom multiples, receiver side peg-leg multiples, and source side peg-leg multiples.
The noise attenuation technique is suitable for eliminating noise from shear waves. Shear waves are conventionally collected using geophones. It has been claimed that converted waves or shear waves can be obtained from processing conventional streamer data (The Diamond Companies Announce Two More Breakthrough 3-D Technologies, The Leading Edge, v. 20, No. 1, pp. 22-23, January 2001). Accordingly, this technique can be used to attenuate noise from conventional streamer data. Example
Figure 6 is a flow chart of a preferred embodiment of the present invention for attenuating noise in two-component ocean bottom cable data using an f-k domain filter. Figures 7(a)-(e) are graphs of the data illustrating the processing steps. As shown in Figure 6, the seismic traces are sorted from common-shot domain to common-receiver domain 51. Then a normal-moveout correction is performed using the velocity of water so that the signal is overcorrected and the noise is undercorrected 53. A forward f-k transform is performed on the normal-moveout corrected seismic traces 55. The transformed data traces are filtered to remove the noise by rejecting the undercorrected events 57. Finally, an inverse f-k transform is performed to transform the data into the original x-t two-dimensional space 59 and an inverse moveout correction may be applied 60.
Figures 7(a)-(e) illustrate the proposed method. Figure 7(a) is an x-t graph illustrating the velocity difference between the signal 61 and noise 63 in common- receiver domain. The figure illustrates that both the noise and signal are coherent in the common-receiver order. It is the separation in velocities that enables the noise or signal to be filtered. A normal-moveout correction is applied to the data in common- receiver order. Figure 7(b) shows that after the normal-moveout correction the noise 67 is undercorrected while the signal 65 is overcorrected. After the normal-moveout correction the multi-trace transformation of the seismic data traces is performed. In this example, a f-k transform filter is used. As shown in Figure 7(c), energy of the undercorrected coherent noise 69 is concentrated in the lower half plane of negative spatial frequency or negative k-parameter, whereas the overcorrected signal 71 is localized in the positive k half plane.
While the seismic data traces are in the two-dimensional domain, the noise or signal may be filtered. In this example, filtering involves zeroing of the negative half plane of the initial f-k spectrum or suppression of the coherent noise. Figure 7(d) illustrates the results of zeroing the negative half plane to suppress the noise but not the signal 71. The filter design or definition of the area of suppression is very simple and would have a universal character if a normal moveout with a velocity of 1500 m/sec is performed. Therefore, using a normal moveout correction of approximately 1500 m sec would enable a specifically designed filter and filtering procedure to be utilized worldwide. An inverse multitrace transform is performed which in this example is an inverse f-k transform which yields common-receiver gather with noise significantly reduced as shown in Figure 7(e). Figure 7(e) illustrates the signal 75 sorted back to common-receiver domain without the noise.
Figures 8-10 illustrate application of the proposed method to a dataset. The initial common-shot data are shown in figures 8(a) and 8(b). The hydrophone data, in figure 8(a), exhibits more random noise compared to the geophone data in figure 8(b). Figures 9(a) and 9(b) respectively, illustrate the same data from figures 8(a) and 8(b) after the data is sorted to common-receiver order. The behavior of the hydrophone signal in figure 9(a) is similar to the hydrophone data in common-shot trace order
illustrated in figure 8(a). However, the geophone data exhibits a coherent noise pattern in figure 9(b) that is different than the common-shot order in Figure 8(b). Therefore, application of the proposed method permits separation of at least part of the noise on the geophone. Figure 10 shows geophone seismic data traces that have been sorted into common-shot order from common-receiver order as shown in figure 9(b) after processing using the inventive noise separation technique explained in the example. Figure 10 illustrates significantly less noise on the vertical geophone data than before the noise separation technique was applied on the vertical geophone data as shown in figure 9(b). Therefore, application of the proposed method significantly reduces the noise on the geophone.
The embodiments described and as applied in the modeled and processed results are intended to clarify the application of the invention to both modeled and real world data. Those skilled in the art will appreciate that other embodiments of the invention can be devised which do not depart from the spirit of the invention as disclosed herein. Accordingly, the scope of the invention should be limited only by the following claims.
Claims (35)
1. A method for attenuating noise in seismic data traces, at least a portion of the noise being noncoherent in the common-shot order but coherent in the common- receiver order, said coherent noise having a move-out velocity different from the move-out velocity of the data signal in said seismic data traces, said method comprising: sorting said seismic data traces into common-receiver order; and using said difference in move-out velocities to separate at least a portion of said coherent noise from said data signal.
2. The method of claim 1 further comprising correcting said seismic data traces after sorting said seismic data traces into common-receiver order by performing a normal moveout of said data with a velocity chosen to overcorrect the signal and undercorrect at least a portion of the noise.
3. The method of claim 2 wherein said normal move out correction of said seismic data traces is performed using a move-out velocity approximately equal to the velocity of sound in water.
4. The method of claim 2 said method further comprising applying an inverse of said normal-moveout correction to said seismic traces after separating at least a portion of said coherent noise.
5. The method of claim 1 wherein said seismic data traces comprise geophone seismic data traces.
6. The method of claim 1 wherein said seismic data traces comprise hydrophone seismic data traces.
7. The method of claim 1 wherein said seismic data traces comprise both geophone and hydrophone seismic data traces.
8. The method of claim 1 wherein said seismic data traces comprise both geophone and hydrophone seismic data traces said method further comprising combining geophone and hydrophone seismic data traces after separating at least a portion of said coherent noise from said seismic data traces by applying a summation of the geophone and hydrophone seismic data traces.
9. The method of claim 1 wherein said seismic data traces are collected using a two-component ocean bottom cable.
10. The method of claim 1 wherein said seismic data traces are collected using a four-component ocean bottom cable.
11. The method of claim 1 wherein said seismic data traces are collected using an ocean streamer.
12. The method of claim 1 wherein said step of using said difference in move-out velocities to separate at least a portion of said coherent noise from said data signal further comprises: applying a move-out correction to said seismic data traces using a move-out velocity chosen to overcorrect said data signal and undercorrect said coherent noise; transforming said move-out corrected seismic data traces from the x-t domain to a two-dimensional domain in which said data signal and said coherent noise are separated; removing at least a portion of said coherent noise from said transformed seismic data traces; and inverse transforming said seismic data traces from said two-dimensional domain back to the x-t domain.
13. The method of claim 12 wherein said moveout correction of said seismic data traces is performed using a move-out velocity approximately equal to the velocity of sound in water.
14. The method of claim 12 wherein said method further comprises applying an inverse of said normal moveout correction after said data traces have been transformed from said two-dimensional domain back to the x-t domain.
15. The method of claim 12 wherein said seismic data traces comprise geophone seismic data traces.
16. The method of claim 12 wherein said seismic data traces comprise hydrophone seismic data traces.
17. The method of claim 12 wherein said seismic data traces comprise both geophone and hydrophone seismic data traces.
18. The method of claim 12 wherein seismic data traces comprise both geophone and hydrophone seismic data traces said method further comprising combining said geophone and hydrophone seismic data traces after removal of at least a portion of said coherent noise from said seismic data traces by applying a summation of the geophone and hydrophone seismic data traces.
19. The method of claim 12 wherein said seismic data traces are collected using a two-component ocean bottom cable.
20. The method of claim 12 wherein said seismic data traces are collected using a four-component ocean bottom cable.
21. The method of claim 12 wherein said seismic data traces are collected using an ocean streamer.
22. The method of claim 12 wherein said two-dimensional domain is selected by choosing at least one filter from the group consisting of the f-k filter, the tau-p filter, and the radon filter.
23. The method of claim 2 wherein said step of using said difference in move-out velocities to separate said coherent noise from said data signal further comprises: transforming said move-out corrected seismic data traces from the x-t domain to a two-dimensional domain in which said data signal and said coherent noise are separated; defining a velocity separation line between said data signal and said coherent noise; removing at least a portion of said coherent noise from said transformed seismic data traces; and inverse transforming said seismic data traces from said two-dimensional domain back to the x-t domain.
24. The method of claim 23 wherein said seismic data traces comprise geophone seismic data traces.
25. The method of claim 23 wherein said seismic data traces comprise hydrophone seismic data traces.
26. The method of claim 23 wherein said seismic data traces comprise both geophone and hydrophone seismic data traces.
27. The method of claim 23 wherein said seismic data traces comprise both geophone and hydrophone seismic data traces and said method further comprising combining geophone and hydrophone seismic data traces after removal of at least a portion of said coherent noise from said seismic data traces by applying a summation of said geophone and hydrophone seismic data traces.
28. The method of claim 23 wherein said seismic data traces are collected using a two-component ocean bottom cable.
29. The method of claim 23 wherein said seismic data traces are collected using a four-component ocean bottom cable.
30. The method of claim 23 wherein said seismic data traces are collected using an ocean streamer.
31. The method of claim 23 wherein said two-dimensional domain is selected by choosing at least one filter from the group consisting of the f-k filter, the tau-p filter, and the radon filter.
32. A method for attenuating noise in an ocean bottom cable (OBC) seismic data traces, said OBC seismic data traces comprising at least vertical geophone data traces, at least a portion of the noise in said data traces being noncoherent in the common- shot order but coherent in the common-receiver order, said coherent noise having a move-out velocity different from the move-out velocity of the data signal in said seismic data traces, said method comprising: sorting said data traces into common-receiver order; applying a move-out correction to said data traces using a move-out velocity chosen to overcorrect said data signal and undercorrect said coherent noises; transforming the move-out corrected data traces from the x-t domain to a two- dimensional domain in which said move-out corrected data signal and said move-out corrected noise are separated; removing at least a portion of said coherent noise from said transformed data traces; and inverse transforming said data traces from said two-dimensional domain back to the x-t domain.
33. The method of claim 32 wherein said move out correction of said data traces is performed using a move-out velocity approximately equal to the velocity of sound in water.
34. The method of claim 32 said method further comprising combining geophone and hydrophone seismic data traces after removal of at least a portion of said coherent noise from said seismic data traces by applying a summation of the geophone and hydrophone seismic data traces.
35. The method of claim 32 wherein said two-dimensional domain is selected by choosing at least one filter from the group consisting of the f-k filter, the tau-p filter, and the radon filter.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US32250401P | 2001-09-14 | 2001-09-14 | |
US60/322,504 | 2001-09-14 | ||
PCT/US2002/028609 WO2003025623A1 (en) | 2001-09-14 | 2002-09-09 | Method for attenuating noise in seismic data |
Publications (2)
Publication Number | Publication Date |
---|---|
AU2002324926A1 true AU2002324926A1 (en) | 2003-06-05 |
AU2002324926B2 AU2002324926B2 (en) | 2007-01-18 |
Family
ID=23255172
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
AU2002324926A Ceased AU2002324926B2 (en) | 2001-09-14 | 2002-09-09 | Method for attenuating noise in seismic data |
Country Status (9)
Country | Link |
---|---|
US (1) | US6738715B2 (en) |
EP (1) | EP1436642A1 (en) |
AU (1) | AU2002324926B2 (en) |
BR (1) | BR0211873A (en) |
CA (1) | CA2453280C (en) |
MX (1) | MXPA04002452A (en) |
MY (1) | MY129374A (en) |
NO (1) | NO20041040L (en) |
WO (1) | WO2003025623A1 (en) |
Families Citing this family (52)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7366054B1 (en) | 2002-08-30 | 2008-04-29 | John M. Robinson | Tau-P filters for removal of noise from seismic data |
US6691039B1 (en) * | 2002-08-30 | 2004-02-10 | John M. Robinson | Removal of noise from seismic data using improved radon transformations |
US6898148B2 (en) * | 2003-03-26 | 2005-05-24 | Westerngeco, L.L.C. | Multi-step receiver-motion compensation |
US7561493B2 (en) * | 2003-05-30 | 2009-07-14 | Fairfield Industries, Inc. | Method and apparatus for land based seismic data acquisition |
US7310287B2 (en) | 2003-05-30 | 2007-12-18 | Fairfield Industries Incorporated | Method and apparatus for seismic data acquisition |
NO20041565L (en) * | 2004-04-16 | 2005-10-17 | Selantic | Seismic subsea source and method for 4-component seismic |
US7739051B2 (en) * | 2004-07-14 | 2010-06-15 | Compagnie Generale De Geophysique | Method for determination of diffractor locations at sea bottom for the purpose of attenuating such energy |
US8534959B2 (en) | 2005-01-17 | 2013-09-17 | Fairfield Industries Incorporated | Method and apparatus for deployment of ocean bottom seismometers |
US7239578B2 (en) * | 2005-03-03 | 2007-07-03 | John M. Robinson | Removal of noise from seismic data using radon transformations |
US7561491B2 (en) * | 2005-03-04 | 2009-07-14 | Robinson John M | Radon transformations for removal of noise from seismic data |
GB2424481B (en) * | 2005-03-23 | 2008-06-04 | Westerngeco Seismic Holdings | Processing seismic data |
US8127706B2 (en) * | 2005-05-02 | 2012-03-06 | Fairfield Industries Incorporated | Deck configuration for ocean bottom seismometer launch platforms |
US7433265B2 (en) * | 2005-10-04 | 2008-10-07 | Fairfield Industries, Inc. | Converted wave energy removal from seismic data |
US8931335B2 (en) * | 2006-04-07 | 2015-01-13 | Baker Hughes Incorporation | Processing of multi-component induction data in the presence of borehole abnormalities |
US20080008038A1 (en) * | 2006-07-07 | 2008-01-10 | Johan Olof Anders Robertsson | Method and Apparatus for Estimating a Seismic Source Signature |
US7492665B2 (en) | 2006-07-24 | 2009-02-17 | Westerngeco L.L.C. | Seismic data acquisition and source-side derivatives generation and application |
US8082106B2 (en) * | 2007-08-16 | 2011-12-20 | Bp Corporation North America Inc. | 3D surface related multiple elimination for wide azimuth seismic data |
US8730761B2 (en) * | 2007-12-17 | 2014-05-20 | Westerngeco L.L.C. | Attenuating noise in seismic data |
US8781749B2 (en) * | 2007-12-28 | 2014-07-15 | Westerngeco L.L.C. | Attenuating noise in seismic data |
US8611191B2 (en) * | 2008-05-22 | 2013-12-17 | Fairfield Industries, Inc. | Land based unit for seismic data acquisition |
US8553497B2 (en) * | 2008-08-11 | 2013-10-08 | Exxonmobil Upstream Research Company | Removal of surface-wave noise in seismic data |
US8089825B2 (en) * | 2008-08-29 | 2012-01-03 | Pgs Geophysical As | Method of summing dual-sensor towed streamer signals using cross-ghosting analysis |
US8094515B2 (en) * | 2009-01-07 | 2012-01-10 | Westerngeco L.L.C. | Seismic data visualizations |
US8301385B2 (en) * | 2010-06-29 | 2012-10-30 | Chevron U.S.A. Inc. | Shot gather data beamer and debeamer |
US9250344B2 (en) * | 2010-10-13 | 2016-02-02 | The Petroleum Institute | Method to acquire seismic data in shallow sea water environment |
US20120213032A1 (en) * | 2011-02-22 | 2012-08-23 | Cggveritas Services Sa | Method for pz summation of 3-dimensional wide azimuth receiver gathers and device |
US9081118B2 (en) * | 2011-11-21 | 2015-07-14 | Cggveritas Services Sa | Device and method for computing depth velocity variations |
US9448315B2 (en) * | 2011-12-27 | 2016-09-20 | Cgg Services Sa | Device and method for denoising ocean bottom data |
CN103675903A (en) * | 2012-09-07 | 2014-03-26 | 中国石油化工股份有限公司 | Multi-source seismic wave field separating method by utilizing random denoising technology |
CN104570118B (en) * | 2013-10-29 | 2017-12-12 | 中国石油化工股份有限公司 | A kind of automatic identification based on dual factors and the method for removing industry disturbance |
CN105093275B (en) * | 2014-05-09 | 2017-08-18 | 中国石油化工股份有限公司 | Useful signal based on rate pattern extracts ISS wave field separation methods |
BR112017014042A2 (en) * | 2015-03-27 | 2018-01-02 | Halliburton Energy Services Inc | system for attenuating coherent noise from seismic data, method for attenuating coherent noise from seismic data and computer readable media |
US10162072B2 (en) * | 2015-10-26 | 2018-12-25 | Pgs Geophysical As | Marine surveys conducted with multiple source arrays |
CN105527649B (en) * | 2015-11-30 | 2018-06-22 | 中国科学院地质与地球物理研究所 | A kind of more focus blended data separation methods of the highly effective gathering of multiple domain separating for several times |
BR112018070577A2 (en) | 2016-04-07 | 2019-02-12 | Bp Exploration Operating Company Limited | detection of downhole sand ingress locations |
WO2017174750A2 (en) | 2016-04-07 | 2017-10-12 | Bp Exploration Operating Company Limited | Detecting downhole sand ingress locations |
US10288755B2 (en) | 2017-03-28 | 2019-05-14 | Saudi Arabian Oil Company | Seismic processing workflow for broadband single-sensor single-source land seismic data |
EA038373B1 (en) | 2017-03-31 | 2021-08-17 | Бп Эксплорейшн Оперейтинг Компани Лимитед | Well and overburden monitoring using distributed acoustic sensors |
CA3073623A1 (en) | 2017-08-23 | 2019-02-28 | Bp Exploration Operating Company Limited | Detecting downhole sand ingress locations |
US11333636B2 (en) | 2017-10-11 | 2022-05-17 | Bp Exploration Operating Company Limited | Detecting events using acoustic frequency domain features |
US11249214B2 (en) * | 2018-05-30 | 2022-02-15 | Saudi Arabian Oil Company | Noise suppression of seafloor geophone seismic data |
EP4234881A3 (en) | 2018-11-29 | 2023-10-18 | BP Exploration Operating Company Limited | Das data processing to identify fluid inflow locations and fluid type |
GB201820331D0 (en) | 2018-12-13 | 2019-01-30 | Bp Exploration Operating Co Ltd | Distributed acoustic sensing autocalibration |
CN111551993B (en) * | 2019-02-12 | 2023-04-07 | 中国石油天然气股份有限公司 | Method and device for suppressing ringing |
CN112180445A (en) * | 2019-07-03 | 2021-01-05 | 中国石油化工股份有限公司 | Seismic exploration data noise removing method, storage medium and computer system |
CN112444853B (en) * | 2019-08-30 | 2024-04-09 | 中国石油化工股份有限公司 | Ground microseism noise suppression method |
EP4045766A1 (en) | 2019-10-17 | 2022-08-24 | Lytt Limited | Fluid inflow characterization using hybrid das/dts measurements |
WO2021073740A1 (en) | 2019-10-17 | 2021-04-22 | Lytt Limited | Inflow detection using dts features |
WO2021093974A1 (en) | 2019-11-15 | 2021-05-20 | Lytt Limited | Systems and methods for draw down improvements across wellbores |
EP4165284B1 (en) | 2020-06-11 | 2024-08-07 | Lytt Limited | Systems and methods for subterranean fluid flow characterization |
US11703607B2 (en) | 2020-06-15 | 2023-07-18 | Saudi Arabian Oil Company | Determining a seismic quality factor for subsurface formations from a seismic source to a first VSP downhole receiver |
EP4168647A1 (en) | 2020-06-18 | 2023-04-26 | Lytt Limited | Event model training using in situ data |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4380059A (en) | 1980-08-20 | 1983-04-12 | Mobil Oil Corporation | F-K Filtering of multiple reflections from a seismic section |
US4486865A (en) | 1980-09-02 | 1984-12-04 | Mobil Oil Corporation | Pressure and velocity detectors for seismic exploration |
US4860265A (en) | 1988-01-25 | 1989-08-22 | Mobil Oil Corporation | Seismic trace restoration using F-K filtering |
US4884248A (en) | 1988-01-25 | 1989-11-28 | Mobil Oil Corporation | Method of restoring seismic data |
US5191526A (en) | 1988-07-18 | 1993-03-02 | Mobil Oil Corporation | Method for removing coherent noise from seismic data |
EG19158A (en) | 1989-08-25 | 1996-02-29 | Halliburton Geophys Service | System for attenuation of water-column reverberation |
US5067112A (en) | 1991-01-04 | 1991-11-19 | Mobil Oil Corporation | Method for removing coherent noise from seismic data through f-x filtering |
US5138583A (en) * | 1991-05-22 | 1992-08-11 | Halliburton Geophysical Services, Inc. | Noise attenuation method |
US5365492A (en) | 1993-08-04 | 1994-11-15 | Western Atlas International, Inc. | Method for reverberation suppression |
US5550786A (en) | 1995-05-05 | 1996-08-27 | Mobil Oil Corporation | High fidelity vibratory source seismic method |
US5774417A (en) * | 1996-10-25 | 1998-06-30 | Atlantic Richfield Company | Amplitude and phase compensation in dual-sensor ocean bottom cable seismic data processing |
US6049507A (en) | 1997-09-30 | 2000-04-11 | Mobil Oil Corporation | Method and apparatus for correcting effects of ship motion in marine seismology measurements |
US6021090A (en) | 1997-10-22 | 2000-02-01 | Western Atlas International, Inc. | Horizontal and vertical receiver-consistent deconvolution for an ocean bottom cable |
GB9800741D0 (en) | 1998-01-15 | 1998-03-11 | Geco As | Multiple attenuation of multi-component sea-bottom data |
US6148264A (en) | 1998-07-06 | 2000-11-14 | Exxonmobil Upstream Research Company | Method for removing seismic noise caused by external activity |
US6151556A (en) | 1999-06-18 | 2000-11-21 | Mobil Oil Corporation | Method and apparatus for doppler smear correction in marine seismology measurements |
US6314371B1 (en) * | 1999-06-25 | 2001-11-06 | Input/Output, Inc. | Dual sensor signal processing method for on-bottom cable seismic wave detection |
US6292754B1 (en) * | 1999-11-11 | 2001-09-18 | Bp Corporation North America Inc. | Vector recomposition of seismic 3-D converted-wave data |
-
2002
- 2002-02-05 US US10/072,588 patent/US6738715B2/en not_active Expired - Lifetime
- 2002-09-09 BR BR0211873-4A patent/BR0211873A/en not_active Application Discontinuation
- 2002-09-09 AU AU2002324926A patent/AU2002324926B2/en not_active Ceased
- 2002-09-09 CA CA002453280A patent/CA2453280C/en not_active Expired - Fee Related
- 2002-09-09 EP EP02759599A patent/EP1436642A1/en not_active Withdrawn
- 2002-09-09 WO PCT/US2002/028609 patent/WO2003025623A1/en not_active Application Discontinuation
- 2002-09-09 MX MXPA04002452A patent/MXPA04002452A/en active IP Right Grant
- 2002-09-12 MY MYPI20023402A patent/MY129374A/en unknown
-
2004
- 2004-03-11 NO NO20041040A patent/NO20041040L/en not_active Application Discontinuation
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2002324926B2 (en) | Method for attenuating noise in seismic data | |
AU2002324926A1 (en) | Method for attenuating noise in seismic data | |
AU2002310037B2 (en) | Multiple suppression for ocean bottom seismic data | |
CA2491340C (en) | System for combining signals of pressure sensors and particle motion sensors in marine seismic streamers | |
US6775618B1 (en) | Method and system for reducing effects of sea surface ghost contamination in seismic data | |
US4752916A (en) | Method and system for removing the effect of the source wavelet from seismic data | |
CA2622703C (en) | Method for prediction of surface related multiples from marine towed dual sensor seismic streamer data | |
US7768869B2 (en) | Method for deghosting and water layer multiple reflection attenuation in marine seismic data | |
US7433265B2 (en) | Converted wave energy removal from seismic data | |
AU2010201504B2 (en) | Method for calculation of seismic attributes from seismic signals | |
US8937848B2 (en) | Methods and systems to eliminate undesirable variations in time-lapse seismic surveys | |
AU2009225361B2 (en) | Method for optimum combination of pressure and particle motion sensors for a 3-D spread of dual-sensor marine seismic streamers | |
AU2002310037A1 (en) | Multiple suppression for ocean bottom seismic data | |
US5051961A (en) | Method and apparatus for seismic survey including using vertical gradient estimation to separate downgoing seismic wavefield | |
NO339711B1 (en) | Decomposition of seismic wave fields from marine cross-line surveys | |
US10324208B2 (en) | Premigration deghosting for marine streamer data using a bootstrap approach in Tau-P domain | |
GB2448831A (en) | Combining the signals from pressure and particle motion sensors in marine seismic streamers | |
Helgesen et al. | MULTI‐OFFSET ACOUSTIC INVERSION OF A LATERALLY INVARIANT MEDIUM: APPLICATION TO REAL DATA 1 | |
Hoffe et al. | Applications of ocean-bottom seismic recordings |