METHOD OF WIDE AZIMUTH PROFILING (WAP)

Abstract

A seismic survey method comprising a vessel, a seismic acquisition system for collecting geophysical seismic data, a marine navigation system for generating positioning data from the location of the vessel and the location of the seismic acquisition system, a seismic data storage engaged with the seismic acquisition system for collecting and storing the seismic data and a seismic data processor engaged with said seismic data storage for seismic processing of the seismic data. The seismic data is acquired along a non-linear acquisition path or sail line. The data consists of CMP lines that follow the non-linear acquisition path. A binning grid covering the CMP lines of the acquired data such that the in-lines follow parallel to the acquisition path and the cross-lines are perpendicular to the in-lines is created. The binning grid comprises a straight portion and a curved portion. Bins for each portion of the binning grid is calculated.

Claims

1. A method for arranging seismic data acquired along a non-linear acquisition path or sail line such that the in-lines follow the non-linear acquisition path and cross-lines are perpendicular, or near to perpendicular, to the in-lines, the method comprising: creating a binning grid covering the CMP lines of the acquired data, wherein the binning grid comprises a straight portion and a curved portion; and calculating bins for each portion.

2. The method according to claim 1, wherein the bins have an in-line number i.sub.y, cross-line number i.sub.x, a width (dy) and a length (dx) and centre coordinates.

3. The method according to claim 2, wherein the width (dy) of the bins is chosen based on desired resolution of the seismic data.

4. The method according to claim 2, wherein a centre line is chosen to define the cross lines.

5. The method according to claim 2, wherein the length (dx) of the bins is calculated by using the distances from the centre coordinates of the bin to the centre coordinate of the two neighboring bins with the same in-line number nx.

6. The method according to claim 5, wherein the dx(i.sub.y,i.sub.x) value for the bin 4 B(i.sub.y,i.sub.x) is calculated by adding half the distance between CC(i.sub.y,i.sub.x) and CC(i.sub.y,i.sub.x−1) given by $\frac{\sqrt{{\left(N.Math.\left(i_y,i_x\right)-N.Math.\left(i_y,i_x-1\right)\right)}^2+{\left(E.Math.\left(i_y,i_x\right)-E.Math.\left(i_y,i_x-1\right)\right)}^2}}{2}$ to half the distance between CC(i.sub.y,i.sub.x+1) and CC(i.sub.y,i.sub.x) given by $\frac{\sqrt{{\left(N.Math.\left(i_y,i_x+1\right)-N.Math.\left(i_y,i_x\right)\right)}^2+{\left(E.Math.\left(i_y,i_x+1\right)-E.Math.\left(i_y,i_x\right)\right)}^2}}{2}$ Wherein: B(i.sub.y,i.sub.x) is bin for cross-line number i.sub.x and a in-line number i.sub.y; CC(i.sub.y, i.sub.x) is centre coordinates for bin B(i.sub.y,i.sub.x) N is the Northing value; E is the Easting value;

7. The method according to claim 1, wherein the non-linear acquisition path may have any shape and length.

8. A seismic survey method comprising; a vessel; a seismic acquisition system for collecting geophysical seismic data; a marine navigation system for generating positioning data from the location of said vessel and the location of said seismic acquisition system; a seismic data storage engaged with the seismic acquisition system for collecting and storing the seismic data; a seismic data processor engaged with said seismic data storage for seismic processing of the seismic data; wherein the seismic data has been acquired along a non-linear acquisition path or sail line which consists of in-lines that follow the non-linear acquisition path and cross-lines that are perpendicular to the in-lines; creating a binning grid covering the CMP lines of the acquired data, the binning grid comprising a straight portion and a curved portion; and calculating bins for each portion.

9. A machine with a readable storage medium using a program of instructions executable by the machine, to perform the method of claim 1.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0027] These and other characteristics of the invention will become clear from the following description of a preferential form of embodiment, given as a non-restrictive example, with reference to the attached schematic drawings, wherein:

[0028] FIG. 1: Shows a WAP swath binned onto a rectangular grid in the same manner as conventional 3D data.

[0029] FIG. 2: Shows visualization of the problem of a WAP swath that is binned onto a rectangular grid.

[0030] FIG. 3a: Shows a vessel towing a 3D seismic acquisition system acquiring a WAP swath of several in-lines along a non-linear acquisition path.

[0031] FIG. 3b: Shows WAP swath binned according to the invention.

[0032] FIG. 4: Shows section A of FIG. 3b.

[0033] FIG. 5: Shows P-Cable acquisition system acquiring a WAP swath that is binned according to the invention.

DETAILED DESCRIPTION OF A PREFERENTIAL EMBODIMENT

[0034] The following description may use terms such as “horizontal”, “vertical”, “lateral”, “back and forth”, “up and down”, “upper”, “lower”, “inner”, “outer”, “forward”, “rear”, etc. These terms generally refer to the views and orientations as shown in the drawings and that are associated with a normal use of the invention. The terms are used for the reader's convenience only and shall not be limiting.

[0035] FIG. 3a shows a seismic vessel 11 towing a 3D seismic acquisition system with everal streamers 14 in non-linear WAP swath 18 manner. The WAP swath 18 has a width 19 and consists of several CMP lines 21. The invention provides a method of binning and processing seismic data acquired along a non-linear acquisition path such that the in-lines 1 always are parallel to the acquisition path and the cross-lines 2 are perpendicular to the in-lines 1 at any given point. This is in contrast to the standard method where the in-lines and cross-lines are linear and lie along a regular and rectangular grid.

[0036] FIG. 4: Shows section A of FIG. 3b. Based on the layout of the acquisition system and desired resolution as shown in FIGS. 3a and 3b, a bin 4 size is chosen. The bin 4 width will typically be equal to the distance between CMP lines 21 created by the individual streamers 14, but is not limited to this width and can be wider or narrower to give the dataset other properties. The bin 4 length may also vary dependent on how many CMP points falls into each bin. The bin 4 length may be shorter, longer or equal to the bin width. Each bin 4 is given a center coordinate. All the centre coordinates 5 of the bins 4 making up a single cross-line lie on a linear line which is normal to the in-lines it crosses at the crossing points. The centre coordinates 5 making up the individual in-lines 1 lie along a line parallel to the acquisition path, this line is not necesserily a linear line. All the individual in-lines 1 are parallel to each other. As such, a binning grid 26 is created which will be rectangular 26a when the acquisition path is linear, and the binning grid 26 will be curved 26b if the acquisition path is curved. In the curved parts of the WAP grid 26b the bin size will not be the same along the individual cross-lines 2. The center bins 4 of the individual cross lines 2 that forms the centre in-line 1 will have the same bin size in both the linear and curved portions of the WAP swath 26, while the “inner” bins 23 taking the shorter path in the curved parts of the grid 26b will be shorter and the “outer” bins 24 taking the longer path in the curved parts of the grid 26b will be longer. The bins along the cross-lines in the linear part of the WAP grid 26a will be approximately equal in size. All the bins will have approximately the same width.

[0037] Each CMP point on a CMP line 21 will be assigned to a bin 4, typically the closest one, but not necessarily. The number of CMP's assigned to each bin 4 is defined as the fold of the bin 4. Each bin 4 will typically have an in-line 1 and a cross-line 2 number, a set of coordinates, a bin width (dy) and length (dx) and an azimuth value, among other values.

[0038] FIG. 4 is a view of section A of FIG. 3b where the curvature of the bins (4) are exaggerated relative to FIG. 3b. This figure shows the process of binning the WAP data 100. The process of binning the WAP swath data 100, such that the in-lines 1 follow parallel to the acquisition path and the cross-lines 2 lie normal to the inlines 1, is based on a concept where a centre line 3 is defined and used to define the cross-lines 2. The x- and y-direction in the binning of WAP data 100 are defined such that the x-direction is along the in-lines 1, and the y-direction is along the cross-lines 2. The CMP line 21 (shown in FIG. 3a,b) from one of the central streamers may be used as a centre line 3. This centre line 3 may undergo smoothing before a cross-line 2 spacing (distance between cross-lines) dx is chosen. This dx defines the length of the bins 4. Points along the centre line 3 with a spacing of dx are defined as initial centre coordinates 5′ for the line of bins 4 forming the centre in-line in the WAP grid 26. At each of these initial centre coordinates 5′ a tangent 6 of the centre line 3 is calculated. A cross-line 2 is then defined at each of the initial centre coordinates 5′ on the centre line 3 as a linear line that is normal to the centre line 3, hence the calculated tangent 6. Along these cross-lines 2, a number of centre coordinates 5 are defined with spacing, dy. The total number of centre coordinates 5 per cross-line 2 makes the number of in-lines 1 in the grid 26. Both dx and dy is chosen based on the desired properties of the binned dataset. The centre line 3 is first predefined to be able to define the cross-lines 2. Next, the in-lines 1 can be defined based on the centre coordinates 5 on the cross-lines 2. All the in-lines in a WAP swath 100 have the same number of centre coordinates 5 and thus also the same number of bins 4, which is also the number of cross-lines 2 in the WAP grid 26. And all the individual cross-lines 2 have the same number of centre coordinates 5 and thus also the same number of bins 4, which is also the number of in-lines 1 in the WAP grid 26. A WAP grid then consists of n.sub.y in-lines 1 and n.sub.x cross-lines 2 where the in-line 1 numbers i.sub.y are ranging from 1 to n.sub.y, and the cross-line 2 numbers i.sub.x are ranging from 1 to n.sub.x. Each in-line 1 in the WAP binning grid 26 consists of n.sub.x bins 4 and each cross-line 2 constists of n.sub.y bins 4. Each bin 4 have a centre coordinate 5 and belongs to one in-line 1 and one cross-line 2 and thus have an in-line 1 number and a cross-line number 2. All bins 4 with the same in-line 1 number forms an in-line 1 and all bins 4 with the same cross-line 2 number forms a cross-line 2. All the bins 4 in both the in-lines 1 and the cross-lines 2 are numbered sequentially along the line.

[0039] When the centre coordinates 5 for all the bins 4 in the WAP binning grid are calculated, a dx value for all the bins 4 will be calculated to define the bin length. The dx values for the initial centre coordinates 5′ of the centre line 3 is chosen based on the desired properties of the binned dataset, but the dx values for the bins forming the other in-lines will where the acquisition path is curved not be the same as for the centre in-line and they may vary along the in-lines 1. In linear parts of the swath the dx value may approximately be the same for all the bins 4 along an individual cross-line 2, however, in the curved parts of the swath, the dx value will vary along the cross-line 2 as illustrated with reference number 25. The dx value for the individual bins 4 is calculated by using the distances from the centre coordinate 5 of the bin to the centre coordinates 5 of the two neighboring bins 4 with the same in-line 1 number. This will give a unique dx value for all the bins 4 except for those forming the centre in-line defined by the centre line 3 used to define the cross-lines 2. This calculation is based on simple Pythagoras and the curved nature of the bins is ignored at bin level and the distances are calculated as straight lines between coordinates of neighboring bins. The centre coordinates 5 of the bins 4 consist of a Northing and an Easting, given that the coordinates are given in the Universal Transverse Mercator coordinate system (UTM). The data is however not limited to be represented by this coordinate system. For easier notation the centre coordinates 5 are now shortened CC, the bins B, the Northing N and the Easting E. They will all be linked to both a cross-line 2 number ix and a in-line 1 number i.sub.y like this CC(i.sub.y,i.sub.x). The dx(i.sub.y,i.sub.x) value for the bin 4 B(i.sub.y,i.sub.x) is calculated by adding half the distance between CC(i.sub.y,i.sub.x) and CC(i.sub.y,i.sub.x−1) given by

[00003]$\frac{\sqrt{{\left(N.Math.\left(i_y,i_x\right)-N.Math.\left(i_y,i_x-1\right)\right)}^2+{\left(E.Math.\left(i_y,i_x\right)-E.Math.\left(i_y,i_x-1\right)\right)}^2}}{2}$

to half the distance between CC(i.sub.y,i.sub.x+1) and CC(i.sub.y,i.sub.x) given by

[00004]$\frac{\sqrt{{\left(N.Math.\left(i_y,i_x+1\right)-N.Math.\left(i_y,i_x\right)\right)}^2+{\left(E.Math.\left(i_y,i_x+1\right)-E.Math.\left(i_y,i_x\right)\right)}^2}}{2}$

[0040] The bin corners 27 are simply defined as the crossing point between two lines where the first line is defined as a linear line between the middle point between CC(i.sub.y,i.sub.x) and CC(i.sub.y+1,i.sub.x) and the middle point between CC(i.sub.y,i.sub.x+1) and CC(i.sub.y+1,i.sub.x+1), and the second line is defined as a linear line between the middle point between CC(i.sub.y,i.sub.x) and CC(i.sub.y,i.sub.x+1) and the middle point between CC(i.sub.y+1,i.sub.x) and CC(i.sub.y+1,i.sub.x+1). This corner will then be the corner between the four bins 4 B(i.sub.y,i.sub.x), B(i.sub.y,i.sub.x+1), B(i.sub.y+1, i.sub.x) and B(i.sub.y+1, i.sub.x+1).

[0041] A complete binning grid with coordinates for both the centre coordinates 5 of the bins 4, and the corners 27 giving the bins a physical exstent making it possible to decide which traces belong to which bins 4, is now calculated.

[0042] The cross-lines 2 are defined based on the chosen centre line 3, which is a smoothened version of the CMP line created by the central streamer in the acquisition system. However, if the centre line 3 still is too uneven and the “inner” parts of the cross-lines 2 defined to be normal to the centre line 3 are crossing each other in curved parts of the swath 18, a negative dx value will be calculated for some bins 4. This will not be accepted and more smoothing will be applied to the centre line 3 until a positive dx is obtained for all bins 4. Or the trace can be deleted from the binning process.

[0043] When the WAP binning grid 26 is complete, all the traces are assigned to the bin 4 they fall within based on their CMP position. If some traces do not fall into any bin 4 but falls outside the WAP binning grid 26, the WAP binning grid 26 is either recalculated using a modified centre line 3 or the trace is simply assigned to the closest bin 4. This decision will be made based on the number of bins 4 that fall outside the WAP binning grid 26 and how far outside the WAP binning grid 26 they are located.

[0044] Further processing of the WAP data can be performed either by means of 3D processing or 2D processing. Typically, noise removal and smoothing of the data will benefit from data in three dimensions, and for wide swaths or in cases where there is more than one adjacent swath a 3D migration could even be conducted. If it is decided to not utilize the 3D information in the dataset, all or some of the bins with the same cross-line number can be stacked together such that the swath becomes a single in-line which will then be a 2D line. The benefit is that the number of traces stacked together will be relatively large and will give a high signal to noise ratio. The processing steps may further include, but is not limited to; demultiplexing, geometry corrections, editing, amplitude corrections, frequency filters, deconvolution, CMP-sorting, velocity analysis, NMO/DMO-corrections, stacking, migration or any other step known in the seismic processing art.

[0045] The seismic data may be a data acquired by a P-Cable high resolution 3D seismic acquisition system. This system makes it possible to collect many seismic profiles simultaneously in a manner which is simpler than when applying conventional techniques. This system is shown in FIG. 5 in more detail. The system 200 is based around a cross cable 10. This is a cable which is towed perpendicular to the sailing direction of the vessel 11 and is suspended in the water by two paravanes 12. The two is paravanes 12 are towed by the vessel 11 from two tow ropes 13. Several streamers 14 with a short mutual distance are attached at the cross cable 10, normally between 3 and 15 meters dependent of system configuration. These distances between the streamer should however not be seen as a limitation for streamer spacing. The signal from these streamers 14 are digitized in a digitizing unit for each streamer 14 in the water before it is transferred to the acquisition unit on the vessel through the cross cable 10 and a single signal cable 15 on either or both sides of the system 200, normally the starboard side. The fact that the signals from all the individual streamers 14 are transferred through the cross cable 10 and then through the same signal cable 15, instead of through a separate signal cable for each streamer 14 as is the case for large conventional 3D systems, allow the streamers 14 to be attached with such a short streamer distance. Due to this short streamer spacing the streamers 14 are relatively short, typically between 12,5 and 100 meters, in contrast to conventional systems where streamers of 10 kilometer length or longer is not uncommon. A seismic source 16 is deployed straight behind the vessel 11 inside the triangle formed by the cross cable 10 and the two tow ropes 13.

[0046] The streamer 14 layout with closely spaced streamers gives the ability to acquire data with very closely spaced CMP lines wich again allow the data to be binned with a small bin size giving data with very high horizontal data. Combined with high frequency sources a dataset of very high resolution can be abtained.

[0047] Another advantage for this data collection system is that both WAP data and 3D data acquisition have the same configuration. This allows both regional and target specific acquisition in the same survey without modification to the acquisition setup. The WAP seismic data acquired with this system also have more potential than conventional 2D seismic data because it is a narrow 3D cube. This is also a relatively compact and lightweight system that can be operated from a variety of vessels including smaller vessel not purpose built for seismic operations. This again leads to a daily operational cost that is lower than for large conventional 3D operations which utilizes large purpose built vessels to operate the large systems with long streamers and large paravanes. A WAP swath can be acquired with this system for a cost approximately the same as that of a regular 2D line, but with the added benefit of 3D information within the WAP swath, and the ability to acquire smaller proper 3D volumes at specific targets with the same system without modifications.

[0048] It should be understood that a computer program is used to visualize, analyse and process the seismic data accordingly to the invention.

[0049] While the present invention has been described with reference to the illustrated embodiment, it should be understood that numerous changes exist in the details of procedures for accomplishing the desired results, but these shall remain within the field and scope of the invention.