A method and apparatus for acquiring seismic data. In one embodiment, the method includes: moving a first air gun array in the water at a first depth and a second air gun array in the water at a second depth greater than the first depth, in which the total volume of the first air gun array is less than the total volume of the second air gun array, in which the first air gun array is separated from the second air gun array by a distance substantially equal to a shot point interval, firing seismic energy through the first and second air gun arrays through the water into the earth, and recording seismic signals reflected from strata in the earth beneath the water.

Patent
   RE45599
Priority
May 04 2004
Filed
Oct 07 2013
Issued
Jul 07 2015
Expiry
Jun 17 2024
Assg.orig
Entity
Large
0
38
all paid
1. A method for acquiring seismic data, comprising:
moving a first air gun array in the water at a first depth and a second air gun array in the water at a second depth greater than the first depth, wherein the total volume of the first air gun array is less than the total volume of the second air gun array, wherein the first air gun array is separated from the second air gun array by a distance substantially equal to a shot point interval;
firing seismic energy by the first and second air gun arrays through the water into the earth; and
recording a plurality of seismic signals reflected from strata in the earth beneath the water.
11. A marine surveying arrangement, comprising:
a first air gun array tuned to a first bubble oscillation, wherein the first air gun array has a first total volume;
a second air gun array tuned to a leading peak, wherein the second air gun array has a second total volume that is greater than the first total volume, wherein the second air gun array is horizontally separated from the first air gun array by a distance substantially equal to the speed at which the air gun arrays are moving multiplied by a delay time between the firing of the first air gun array and the firing of the second air gun array; and
one or more seismic streamers having a plurality of hydrophones disposed therealong.
2. The method of claim 1, wherein the first air gun array is tuned to a leading peak and the second air gun array is tuned to a first bubble oscillation.
3. The method of claim 1, wherein the first air gun array is separated from the second air gun array by about 37.5 meters.
4. The method of claim 1, wherein moving the first and second air gun arrays comprises moving the first and second air gun arrays at a speed V in an x direction, wherein the first air gun array is separated from the second air gun array by a distance dx in the x direction.
5. The method of claim 4, wherein firing seismic energy through the first and second air gun arrays comprises:
firing seismic energy through the first air gun array; and
firing seismic energy through the second air gun array at a time dx/V after the first air gun array is fired.
6. The method of claim 4, wherein the first air gun array is disposed closer to a seismic survey vessel than the second air gun array, wherein the first and second air gun arrays are moved by the seismic survey vessel.
7. The method of claim 1, wherein firing seismic energy through the first and second air gun arrays comprises firing the first and second air gun arrays sequentially.
8. The method of claim 1, wherein firing seismic energy through the first and second air gun arrays comprises firing the first and second air gun arrays simultaneously.
9. The method of claim 1, wherein the first air gun array is configured to enhance high frequency seismic energy.
10. The method of claim 1, wherein the second air gun array is configured to enhance low frequency seismic energy.
12. The marine surveying arrangement of claim 11, wherein the first total volume is about 5085 cubic inches and the second total volume is about 6780 cubic inches.
13. The marine surveying arrangement of claim 11, wherein the distance is about 37.5 meters.

By imposing the flat sea boundary condition in the form U(0)=−D(0), i.e. a surface reflectivity of −1, for a data window below a direct arrival, the following expression may be derived:

P ( Z o ) [ j2 π f 1 - k 2 v 2 / f 2 Z u ( 1 - - j2 π f 1 - k 2 v 2 / f 2 Z u ) ] = P ( Z u ) [ j2 π f 1 - k 2 v 2 / f 2 Z o ( 1 - - j2 π f 1 - k 2 v 2 / f 2 Z o ) ] .
Persons of ordinary skill in the art should appreciate that the expressions in brackets are ghost operators, FO and FU, for the over and under seismic data, respectively, in the case of a perfectly calm surface 530. Accordingly, the above expression states that the pressure at the “over” seismic receiver 540 multiplied by the ghost operator FU of the “under” seismic receiver 550 is equal to the pressure at the “under” seismic receiver 550 multiplied by the ghost operator Fo of the “over” seismic receiver 540. In mathematical terms, the above expression may be written in the simplified form: P(Zo)FO=P (Zu)FU.

However, as discussed above, the surface 530 is virtually never flat, as assumed above and in conventional practice. Moreover, the above expressions do not account for temporal and spatial variations in the water velocity, reflectivity of the surface 530, streamer positioning errors, and other non-ideal conditions that are frequently encountered in real marine seismic surveys. To account, at least in part, for the effects of the aforementioned non-ideal conditions, one or more calibration filters are determined in a manner that will be discussed in detail below. The calibration filters are then used to form an over/under combination of marine seismic data acquired by the “over” seismic receiver 540 and the “under” seismic receiver 550. For example, the over/under combination may be formed by modifying the surface boundary condition using the one or more calibration filters. The over/under combination formed with the modified boundary condition may result in a combined data set with reduced noise relative to a data set formed by an over/under combination using the flat sea boundary condition.

FIG. 6 illustrates a flow diagram of a method 600 for forming an over/under combination using one or more calibration filters in accordance with one or more embodiments of the invention. First and second data sets are selected (at 610). In one embodiment, the first and second data sets are selected (at 610) to be pre-stack over and under data sets acquired by at least one seismic sensor coupled to an “over” streamer and at least one seismic sensor coupled to an “under” streamer in an over/under streamer combination. However, the present invention is not limited to selecting (at 610) all of the data in the pre-stack data set. In one alternative embodiment, portions of the pre-stacked data set acquired within a selected time window and/or a selected offset may be selected (at 610). In another alternative embodiment, portions of the pre-stacked data set from a selected gather, such as a shot gather and/or a receiver gather, may be selected (at 610).

The first and/or second data sets may be provided via transmission over a wired and/or wireless medium. For instance, the over and under data sets may be selected from the data as it is gathered, or shortly after it is collected, from a seismic survey. Alternatively, the first and/or second data sets may be recorded on and transmitted via recording tape, magnetic disks, compact disks, DVDs, and the like. Thus, the first and second data sets can, in some embodiments, can be selected from data previously collected and archived on some magnetic or optical storage medium.

One or more calibration filters are determined (at 620) using the selected over and under data sets. In one embodiment, the one or more calibration filters are determined (at 620) by initially assuming, as discussed above, the pressure at the “over” seismic receiver 440 multiplied by the ghost operator FU of the “under” seismic receiver 450 is equal to the pressure at the “under” seismic receiver 450 multiplied by the ghost operator FO of the “over” seismic receiver 440, i.e. P(Zo)FO=P(Zu)FU. This technique is often referred to as across-ghosting technique.

However, as discussed above, this relationship generally is not precise for the acquired over/under seismic data. The one or more calibration filters, a(f), may therefore be determined using the expression a(f)P(Zo)FO=P(Zu)FU. For example, the one or more calibration filters may be determined by evaluating the expression a(f)P(Zo)FO=P(Zu)FU by a least-squares criterion. However, persons of ordinary skill in the art should appreciate that the present invention is not limited to applying the least-squares criterion to the expression a(f)P(Zo)FO=P (Zu)FU. Further, any desirable expression may be evaluated with any desirable technique used to determine the calibration filters. Persons of ordinary skill in the art should also appreciate that the one or more calibration filters may be determined such that the expression a(f)P(Zo)FO=P(Zu)FU holds true in a statistical sense, even though it may not hold precisely for all the acquired seismic data used to determine the calibration filters.

The one more calibration filters are then used to combine (at 630) the first and second data sets to form a third data set, such as an over/under combined data set. The one or more calibration filters may be used to define a perturbed boundary condition using over and under data sets. The perturbed boundary condition is then incorporated into a selected over/under combination technique that is used to combine (at 630) the over and under seismic data. Persons of ordinary skill in the art should appreciate that the present invention is not limited to any particular technique for combining (at 630) the over and under seismic data using the one or more calibration filters. In various alternative embodiments, any desirable technique for combining (at 630) the over and under seismic data using the one or more calibration filters may be used.

While the foregoing is directed to embodiments of the present invention, other and further embodiments of the invention may be devised without departing from the basic scope thereof, and the scope thereof is determined by the claims that follow.

Moldoveanu, Nicolae

Patent Priority Assignee Title
Patent Priority Assignee Title
3744021,
4493061, Aug 10 1981 WELLS FARGO BANK TEXAS , N A Stacked marine seismic source
4648080, Jun 15 1984 WESTERNGECO, L L C Method for determining the far field signature of a marine seismic source from near-field measurements
4721180, Nov 26 1986 WESTERNGECO, L L C Marine seismic source array
4739858, Mar 02 1987 WESTERNGECO, L L C Spectrally-shaped air gun arrays
4813026, Nov 27 1987 Mobil Oil Corporation Method for logarithmic analysis of seismic reflection signals
4868794, May 29 1981 BRITOIL PUBLIC LIMITED COMPANY, A BRITISH COMPANY; MERLIN GEOPHYSICAL LIMITED, A BRITISH COMPANY; GECO GEOPHYSICAL COMPANY OF NORWAY A S , A NORWEGIAN COMPANY Method of accumulation data for use in determining the signatures of arrays of marine seismic sources
4956822, Dec 09 1988 SEASCAN, INC , A CORP OF TX Method and apparatus for seismic exploration
4992990, Jun 06 1988 Geco A.S. Method for determining the position of seismic streamers in a reflection seismic measuring system
5142498, Aug 28 1991 ExxonMobil Upstream Research Company Controlled phase marine source array
5148406, May 31 1989 GECO A S , A CORP OF NORWAY Method for simultaneous collection of seismic data from shallow and deep targets
5281773, Aug 28 1991 ExxonMobil Upstream Research Company Controlled phase marine source subarray
5717655, Apr 18 1995 WESTERNGECO, L L C Method for providing uniform subsurface coverage in the presence of steep dips
5924049, Apr 18 1995 WESTERNGECO, L L C Methods for acquiring and processing seismic data
5995452, Aug 02 1995 Hydroacoustics, Inc. System for generating and transmitting acoustic signals underwater
6044038, Jun 08 1998 WESTERNGECO, L L C Marine seismic cable system
6446008, May 20 1998 Schlumberger Technology Corporation Adaptive seismic noise and interference attenuation method
6493636, Nov 05 1998 Shell Oil Company Method of marine seismic exploration utilizing vertically and horizontally offset streamers
6651007, May 20 1998 Schlumberger Technology Corporation Adaptive seismic noise and interference attenuation method
6671223, Dec 20 1996 WESTERNGECO, L L C Control devices for controlling the position of a marine seismic streamer
6684160, May 20 1998 Reflection Marine Norge AS Marine seismic acquisition system and method
6961284, Apr 03 2000 WESTERNGECO L L C Source array for use in marine seismic exploration
7577060, Apr 08 2005 WesternGeco L.L.C. Systems and methods for steering seismic arrays
7800977, Jun 01 2004 WesternGeco L.L.C. Pre-stack combining of over/under seismic data
7948825, Jun 17 2004 WESTERNGECO, L L C Method to acquire simultaneously seismic data with source arrays designed for specific targets
7961549, Jun 16 2004 WESTERNGECO, L L C Enhancing the acquisition and processing of low frequencies for sub-salt imaging
8559264, Jun 17 2004 WesternGeco L.L.C. Method to acquire simultaneously seismic data with source arrays designed for specific targets
20020126565,
20030067842,
20080011540,
AU2004319619,
EP320073,
GB2141824,
RE32683, Aug 10 1981 WELLS FARGO BANK TEXAS , N A Stacked marine seismic source
WO20895,
WO116621,
WO175481,
WO2005111656,
/
Executed onAssignorAssigneeConveyanceFrameReelDoc
Oct 07 2013WesternGeco L.L.C.(assignment on the face of the patent)
Date Maintenance Fee Events
Nov 14 2018M1552: Payment of Maintenance Fee, 8th Year, Large Entity.
Nov 09 2022M1553: Payment of Maintenance Fee, 12th Year, Large Entity.


Date Maintenance Schedule
Jul 07 20184 years fee payment window open
Jan 07 20196 months grace period start (w surcharge)
Jul 07 2019patent expiry (for year 4)
Jul 07 20212 years to revive unintentionally abandoned end. (for year 4)
Jul 07 20228 years fee payment window open
Jan 07 20236 months grace period start (w surcharge)
Jul 07 2023patent expiry (for year 8)
Jul 07 20252 years to revive unintentionally abandoned end. (for year 8)
Jul 07 202612 years fee payment window open
Jan 07 20276 months grace period start (w surcharge)
Jul 07 2027patent expiry (for year 12)
Jul 07 20292 years to revive unintentionally abandoned end. (for year 12)