Crosstalk attenuation for seismic imaging can include creation of a seismic image based on seismic data including multiples. The seismic image can include causal crosstalk and anti-causal crosstalk. Causal crosstalk and anti-causal crosstalk can be predicted based on the seismic data. The predicted causal crosstalk and the predicted anti-causal crosstalk can be attenuated from the seismic image.
|
12. A non-transitory machine-readable medium storing instructions executable by a processing resource to:
create a first seismic image based on first seismic data that includes primaries;
create a second seismic image based on second seismic data that includes multiples;
wherein the first seismic image and the second seismic image are indicative of a subterranean formation;
predict first causal crosstalk based on first seismic data;
predict second causal crosstalk and anti-causal crosstalk based on second seismic data;
attenuate the first causal crosstalk from the first seismic image;
attenuate the second causal crosstalk and the anti-causal crosstalk from the second seismic image;
combine the causal crosstalk attenuated first seismic image with the causal crosstalk and anti-causal crosstalk attenuated second seismic image to create a resulting seismic image; and
store the resulting seismic image that is better indicative of the subterranean formation.
16. A method of generating a geophysical data product, the method comprising:
obtaining geophysical data including seismic data;
creating a seismic image based on the seismic data including multiples, wherein the seismic image includes causal crosstalk and anti-causal crosstalk and wherein the seismic image is indicative of a subterranean formation;
predicting causal crosstalk based on the seismic data by:
migration of down-going multiples to a subsurface image point as a source wavefield of the seismic data; and
migration of up-going primaries to the subsurface image point as a receiver wavefield of the seismic data;
predicting anti-causal crosstalk based on the seismic data by:
migration of a down-going wavefield to a subsurface image point as a source wavefield of the seismic data; and
migration of up-going multiples to the subsurface image point as a receiver wavefield of the seismic data;
attenuating the predicted causal crosstalk and the predicted anti-causal crosstalk from the seismic image such that the seismic image is better indicative of the subterranean formation; and
recording the attenuated seismic image on a non-transitory machine-readable medium, thereby generating the geophysical data product.
8. A method, comprising:
creating a seismic image based on seismic data including primaries and multiples, wherein the seismic image includes causal crosstalk and anti-causal crosstalk and wherein the seismic image is indicative of a subterranean formation;
predicting causal crosstalk based on the seismic data by:
migrating down-going primaries and multiples to a subsurface image point as a source wavefield of the seismic data; and
migrating up-going primaries to the subsurface image point as a receiver wavefield of the seismic data;
predicting anti-causal crosstalk based on the seismic data by:
migrating a down-going wavefield as a source wavefield of the seismic data from a seismic source modeled as a point source to a subsurface image point; and
migrating up-going multiples as a receiver wavefield of the seismic data from a seismic receiver to the subsurface image point;
attenuating the predicted causal crosstalk and the predicted anti-causal crosstalk from the seismic image; and
storing the seismic image having the predicted causal crosstalk and the predicted anti-causal crosstalk attenuated therefrom in a tangible machine-readable medium such that the seismic image is better indicative of the subterranean formation.
1. A system, comprising:
a seismic imaging engine, comprising hardware, configured to create a seismic image based on seismic data, wherein the seismic data includes multiples, wherein the seismic image includes causal crosstalk and anti-causal crosstalk, and wherein the seismic image is indicative of a subterranean formation;
a prediction engine, comprising hardware, configured to:
predict causal crosstalk based on the seismic data by:
migration of down-going multiples to a subsurface image point as a source wavefield of the seismic data; and
migration of up-going primaries to the subsurface image point as a receiver wavefield of the seismic data; and
predict anti-causal crosstalk based on the seismic data by:
migration of a down-going wavefield to a subsurface image point as a source wavefield of the seismic data; and
migration of up-going multiples to the subsurface image point as a receiver wavefield of the seismic data;
an attenuation engine, comprising hardware, configured to attenuate the predicted causal crosstalk and the predicted anti-causal crosstalk from the seismic image; and
a tangible machine-readable medium storing the seismic image having the predicted causal crosstalk and the predicted anti-causal crosstalk attenuated therefrom such that the seismic image is better indicative of the subterranean formation.
2. The system of
extract the down-going multiples using down-going and up-going wavefield separation of the seismic data; and
estimate the up-going primaries using surface related multiple elimination.
3. The system of
4. The system of
5. The system of
migration of down-going primaries and multiples to a subsurface image point as a source wavefield of the seismic data; and
migration of up-going primaries and multiples to the subsurface image point as a receiver wavefield of the seismic data.
6. The system of
7. The system of
9. The method of
migrating down-going primaries and multiples to a subsurface image point as a source wavefield of the seismic data; and
migrating up-going primaries and multiples to the subsurface image point as a receiver wavefield of the seismic data.
10. The method of
11. The method of
13. The medium of
predict multiples in the first seismic data; and
attenuate multiples from the first seismic data based on the predicted multiples.
14. The medium of
15. The medium of
|
This application claims priority to U.S. Provisional Application No. 62/116,749, filed Feb. 16, 2015, which is incorporated by reference.
In the past few decades, the petroleum industry has invested heavily in the development of marine seismic survey techniques that yield knowledge of subterranean formations beneath a body of water in order to find and extract valuable mineral resources, such as oil. High-resolution seismic images of a subterranean formation are helpful for quantitative seismic interpretation and improved reservoir monitoring. For a typical marine seismic survey, a marine survey vessel tows one or more seismic sources below the sea surface of the water and over a subterranean formation to be surveyed for mineral deposits. Seismic receivers may be located on or near the seafloor, on one or more streamers towed by the marine survey vessel, or on one or more streamers towed by another vessel. The marine survey vessel typically contains marine seismic survey equipment, such as navigation control, seismic source control, seismic receiver control, and recording equipment. The seismic source control may cause the one or more seismic sources, which can be air guns, marine vibrators, etc., to produce acoustic signals at selected times. Each acoustic signal is essentially a sound wavefield that travels down through the water and into the subterranean formation. At each interface between different types of rock, a portion of the wavefield may be refracted, and another portion may be reflected, which may include some scattering, back toward the body of water to propagate toward the sea surface. The seismic receivers thereby measure a wavefield that was initiated by the actuation of the seismic source.
The present disclosure is related to crosstalk attenuation for seismic imaging. A seismic source can emit an acoustic signal. Examples of seismic sources include air guns and marine vibrators, among others. Pressure and particle motion variation as a function of time and position caused by an acoustic signal from a seismic source or modeled as being emitted by a modeled seismic source is called the “source wavefield.” One or more seismic sources can be modeled as a point source. Pressure and particle motion variation as a function of time and position measured by a seismic receiver or modeled as being received by a modeled seismic receiver is called the “receiver wavefield.”
Goals of seismic processing can include mathematically transforming recorded reflections into seismic images of the earth's subsurface. For example, seismic processing methods can include mathematically simulating wavefield propagation using a computer, where boundary data from a sea surface can be extrapolated into a model of the subsurface. As used herein, propagation is the movement of a wavefield, extrapolation is a simulation of propagation, and migration is a process by which an image is produced through extrapolation of boundary data in space and/or time to another location, such as the subsurface, to produce a seismic image. The boundary data at the sea surface can consist of two parts: a source wavefield, and a receiver wavefield. Some wave equation seismic imaging methods can include extrapolating source and receiver wavefields from the sea surface into an earth model (e.g., subsurface model) and can produce a seismic image by computing the locations where a source wavefield and a receiver wavefield are in phase.
Primary wavefields (“primaries”) and multiple wavefields (“multiples”) can each be divided into down-going and up-going primaries and multiples respectively. As used herein, a down-going primary is a wavefield that is a reflection of an up-going primary. An up-going primary is a reflection of a wavefield emitted by a seismic source. A down-going multiple is a reflection of an up-going multiple. An up-going multiple is a reflection of a down-going primary or multiple. Some seismic imaging methods only make use of primaries. However, multiples can carry valuable information. According to the present disclosure, it can be beneficial to incorporate multiples in seismic imaging methods. However, correlation based seismic imaging with multiples can generate crosstalk. Crosstalk can be generated by source and receiver wavefields being in phase at locations that are not the same location as a subsurface reflector. Attenuation of crosstalk can be achieved by post-processing the seismic image and/or modifying the seismic image conditions. As described in more detail with respect to
Migration using primaries and multiples can include creating different seismic images by supplying different wavefields as source wavefields and/or receiver wavefields. An example of creating different seismic images by supplying different source and/or receiver wavefields is provided in Table 1. Table 1 lists corresponding seismic images from migration (Image 1 to Image 18). Image 19, Image 20, and Image 21 are created from post-processing combined Images 1-18. Image 19 is a seismic image of multiples after causal crosstalk and anti-causal crosstalk attenuation. Image 20 is a combined seismic image of primaries and multiples after causal crosstalk and anti-causal crosstalk attenuation. Image 21 is a joint seismic image of primaries and multiples after causal crosstalk and anti-causal crosstalk attenuation. Migration can use seismic data, such as source and receiver wavefields, as inputs to create a seismic image of a subsurface reflector. The seismic data may be acquisition data. Crosstalk can be predicted without first creating a seismic image that includes the crosstalk. If the source and receiver wavefields are the same (e.g., if the source and receiver wavefields are primaries, or if the source and receiver wavefields are first order multiple, etc.), the migration can output crosstalk. This type of crosstalk is categorized as zeroth (0th) order crosstalk. Other orders of crosstalk are described herein.
TABLE 1
Source wavefields, receiver wavefields, and seismic images from
combination of different source receiver wavefield pairs.
Receiver
Primaries +
Primaries
Multiples
Multiples
(Receiver A)
(Receiver B)
(Receiver C)
Source
Image
Point Source
Primaries
Causal Crosstalk
Primaries + Causal
[for Direct Arrival]
(Image 1)
(Image 2)
Crosstalk
(Source A)
(Image 3)
Primaries
Anti-causal Crosstalk
Multiples + Causal
Multiples + Causal
(Source B)
(only 0th order)
Crosstalk
Crosstalk + Anti-
(Image 4)
(Image 5)
causal Crosstalk
(only 0th order)
(Image 6)
Multiples
Anti-causal Crosstalk
Multiples + Causal
Multiples + Causal
(Source C)
(except for 0th order)
Crosstalk + Anti-
Crosstalk + Anti-
(Image 7)
causal Crosstalk
causal Crosstalk
(Image 8)
(Image 9)
Primaries +
Anti-causal Crosstalk
Multiples + Causal
Multiples + Causal
Multiples
(Image 10)
Crosstalk + Anti-
Crosstalk + Anti-
(Source D)
causal Crosstalk
causal Crosstalk
(Image 11)
(Image 12)
Point Source +
Primaries + Anti-
Multiples + Causal
Primaries + Multiples +
Primaries +
causal Crosstalk
Crosstalk + Anti-
Causal Crosstalk +
Multiples
(Image 13)
causal Crosstalk
Anti-causal Crosstalk
(Source E)
(Image 14)
(Image 15)
Point Source +
Primaries + Anti-
Multiples + Causal
Primaries + Multiples +
Primaries
causal Crosstalk
Crosstalk
Causal Crosstalk +
(Source F)
(only 0th order)
(Image 17)
Anti-causal Crosstalk
(Image 16)
(only 0th order)
(Image 18)
TABLE 2
Terminology versus symbol of variables and operators
Symbol
Object (variable and operator)
PU
Up-going wavefield
PUp
Up-going primary
PUm
1st order up-going multiple
PUm
2nd order up-going multiple
PD
Down-going wavefield
PD1
Down-going wavefield from a point source or down-going
wavefield from impulse wavelet
PDp
Down-going primary
PDm
1st order down-going multiple
PDm
2nd order down-going multiple
I(x)
Migrated image
G*(x, xs)
Green's function from the source location (xs) to the
subsurface location (x)
G(x, xr)
Green's function from the receiver location (xr) to the
subsurface location (x)
G*(xr, x)
Green's function from the subsurface location (x) to the
receiver location (xr)
R(x)
Reflectivity
x
Subsurface location (vector)
xs
Source location (vector)
xr
Receiver location (vector)
ω
Temporal frequency
It is to be understood the present disclosure is not limited to particular devices or methods, which may, of course, vary. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting. As used herein, the singular forms “a”, “an”, and “the” include singular and plural referents unless the content clearly dictates otherwise. Furthermore, the word “may” is used throughout this application in a permissive sense (i.e., having the potential to, being able to), not in a mandatory sense (i.e., must). The term “include,” and derivations thereof, mean “including, but not limited to.” The term “coupled” means directly or indirectly connected.
The figures herein follow a numbering convention in which the first digit or digits correspond to the drawing figure number and the remaining digits identify an element or component in the drawing. Similar elements or components between different figures may be identified by the use of similar digits. For example, 109 may reference element “09” in
Acquisition and processing techniques can be used to extract up-going and down-going wavefields. In a marine setting, dual-sensor (e.g., hydrophone and/or vertical geophone, etc.) and/or down-going and up-going wavefield separation can be used to extract such wavefields. Approximations of up-going and/or down-going wavefields can be determined by other methods of deghosting. Up-going and down-going wavefields can be represented by equation (1) as follows:
where PU represents an up-going wavefield and PD represents a down-going wavefield. PUP represents the primaries of the up-going wavefield, PUm
The components of PU and PD are listed in Table 2. In a migration, PD can be used as an input for a source wavefield, and PU can be used as an input for a receiver wavefield. When different primary and multiple components in PU and PD are combined during migration, different seismic images can be created as shown in Table 1.
Migration using down-going wavefield data as a generalized source wavefield and up-going wavefield data as a generalized receiver wavefield can be represented by Equation 2 as follows:
I(x)=ΣωΣx
where G(x, xr; ω) represents the Green's function from a receiver at the surface, xr, to an image point x in a subsurface, G*(x, xs; ω) represents the Green's function from a source at the surface, xs, to an image point x in a subsurface, and PU represents and up-going wavefield and PD represents a down-going wavefield.
A receiver wavefield can be simulated given a source wavefield and a subsurface reflection property (e.g., velocity, reflectivity, impedance, etc.). A one-way finite difference operator can be used to migrate a wavefield in the frequency domain. Embodiments are not, however, limited to a one-way wave equation operator. For example, a two-way wave equation or any ray based operator can be used.
Numerical modeling can include using a down-going wavefield as the source wavefield. Simulated up-going wavefield data at receivers can be represented by Equation 3 as follows:
PU(xr;ω)=G*(x,xs;ω)P*D(xs;ω)R(x)G*(xr,x;ω) (3)
where G*(x, xs; ω) represents the Green's function from the source location (xs) to the sea surface reflection point (x), G*(xr, x; ω) is the Green's function from the subsurface reflection point (x) to the receiver location (xr), P*D (xs; ω) is the down-going wavefield used as a source wavefield at the sea surface location (xs), R(x) is the reflectivity at the subsurface location (x), and PU(xr; ω) is the simulated up-going wavefield at the receiver location (xr).
A point source 232 is illustrated emitting a down-going point source wavefield 237, which is illustrated as a ray. The down-going point source wavefield 237 is illustrated reflecting off of the surface 204 as a first up-going primary 240-1. The point where the down-going point source wavefield 237 and the first up-going primary 240-1 are in phase is referred to as a surface image point 235 because it is an image point at the surface 204. Wave equation seismic imaging methods can migrate down-going and up-going wavefields to produce a seismic image where these wavefields are in phase. The first up-going primary 240-1 continues up to the sea surface 209 where it can be recorded by a first seismic receiver 242-1.
The down-going point source wavefield 237 is illustrated as continuing through the surface 204 to a subsurface reflector 234, where it reflects as a second up-going primary 240-2. The point where the down-going point source wavefield 237 and the second up-going primary 240-2 are in phase is referred to as a subsurface image point 236 because it is an image point at the subsurface reflector 234. The second up-going primary 240-2 continues up to the sea surface 209 where it can be recorded by a second seismic receiver 242-2. Although the point source 232, the first seismic receiver 242-1, and the second seismic receiver 242-2 are illustrated at or near the sea surface 209 for clarity, as will be appreciated by one of ordinary skill in the art, they can be below the sea surface 209.
As illustrated in
The down-going primary 338 is illustrated reflecting off of the surface 304 as a first up-going multiple 346-1. The first up-going multiple 346-1 is a first order multiple because it is an up-going reflection of a primary. The down-going primary 338 and the first up-going multiple 346-1 are in phase at a second surface image point 335-2. The down-going primary 338 is also illustrated reflecting off of the subsurface reflector 334 as a second up-going multiple 346-2, which is also a first order up-going multiple. The down-going primary 338 and the second up-going multiple 346-2 are in phase at a second subsurface image point 336-2. The second up-going multiple 346-2 is illustrated reaching the sea surface 309 at a location coincident with a location of a fourth seismic receiver 342-4.
The first up-going multiple 346-1 is illustrated reflecting off of the sea surface 309 at a location coincident with a location of a third seismic receiver 342-3 as a down-going multiple 344. The down-going multiple 344 is a first order down-going multiple because it is a down-going reflection of a first order multiple.
The down-going multiple 344 is illustrated reflecting off of the surface 304 as a third up-going multiple 346-3. The third up-going multiple 346-3 is a second order multiple because it is an up-going reflection of a first order multiple. The down-going multiple 344 and the third up-going multiple 346-3 are in phase at a third surface image point 335-3. The third up-going multiple 346-3 is illustrated reaching the sea surface 309 at a location coincident with a location of a fifth seismic receiver 342-5. The down-going multiple 344 is also illustrated reflecting off of the subsurface reflector 334 as a fourth up-going multiple 346-4, which is also a second order up-going multiple. The down-going multiple 344 and the fourth up-going multiple 346-4 are in phase at a third subsurface image point 336-3. The fourth up-going multiple is illustrated reaching the sea surface 309 at a location coincident with a location of a sixth seismic receiver 342-6.
As opposed to primaries, which may be generally transformed into desired seismic image data, multiples may be generally undesired seismic data because they do not transform into desired seismic image data, but may be generally considered to transform into noise. Thus, some previous approaches make use only of primaries. However, multiples can carry valuable information and can be used, according to the present disclosure, in seismic imaging methods as described herein. Thus, while some previous approaches to seismic imaging may seek to remove multiples, at least one embodiment of the present disclosure retains the information provided by the multiples without removing the multiples or before removing the multiples from the seismic data. However, correlation based seismic imaging with multiples can generate significant crosstalk, such as causal crosstalk and/or anti-causal crosstalk. Such crosstalk can be generated by down-going and up-going wavefields being in phase at locations that do not correspond to geological reflectors. At least one embodiment of the present disclosure reduces such crosstalk.
Causal crosstalk 448 is a false indication of an up-going wavefield being in phase with a down-going wavefield. The false indication can correspond to a location that is deeper than a location where two different wavefields are actually in phase. The false indication can correspond to a time that is later than a time when two different wavefields are actually in phase. The causal crosstalk 448 appears to indicate that the down-going primary 438 is in phase with the apparent up-going multiple 447 at the false reflector 450. However, in this example, neither the false reflector 450, nor the apparent up-going multiple 447 actually exist. The fact that the down-going primary 438 happens to have a point where it is in phase with an extension of the third up-going multiple 446-3 can generate undesired seismic image data (causal crosstalk 448) for a reflector that does not actually exist (the false reflector 450). This causal crosstalk 448 can appear as noise in a seismic image. The causal crosstalk 448 can be said to correspond to a location that is deeper than a location where two different wavefields are actually in phase because it can be defined by its down-going component, which is the down-going primary 438. The causal crosstalk 448 can be said to correspond to a time that is later than a time where two different wavefields are actually in phase because it can be defined by its down-going component, which is the down-going primary 438. The down-going component (which is also the existing component) of the causal crosstalk 448 (as opposed to the apparent up-going multiple 447) is in phase with the existing first up-going multiple 446-1 at the second surface image point 435-2, which is shallower and earlier than the causal crosstalk 448.
Although not specifically illustrated as such,
In a migration, when the seismic source 432 is modeled as a point source that emits a source wavefield, the down-going source (Source A in Table 1) wavefield, and recorded up-going primaries and multiples (Receiver C in Table 1) are used as a receiver wavefield, a seismic image of primaries can be created with causal crosstalk. An example of a created seismic image with causal crosstalk is illustrated in Table 1 as Image 3.
The up-going and down-going wavefields can be decomposed according to Equation 4 as follows:
where PD1 represents the down-going wavefield from a point source, PUp represents the up-going primaries, PUm
By substituting Equation 4 into Equation 2, the created seismic image (Image 3 in Table 1) can include primaries and causal crosstalk, and can be represented by Equation 5 as follows:
where Equation 5.1 describes a seismic image of primaries, Equation 5.2 describes first order causal crosstalk, and Equation 5.3 describes second order causal crosstalk.
A data domain of the up-going multiples (Receiver B in Table 1) can be modeled using down-going primaries and multiples (Source D in Table 1) as a source wavefield. The seismic image from primaries with causal crosstalk (Image 3 in Table 1) can be used for a subsurface reflection property. A seismic imaging engine (e.g., seismic imaging engine 1376 illustrated in
PD=PDp+PDm
where PDp represents down-going primaries, PDm
Substituting Equation 5 and Equation 6, the arrival time of a multiple reflection event can be estimated as represented by Equation 7 as follows:
The modeling of up-going wavefield PU can include all orders of multiples. However embodiments are not limited to including all orders of multiples. Although the amplitude of modeled multiples can be different from recorded multiples, the arrival and timing of the modeled events can be identical to those from seismic acquisition. In at least one embodiment, the up-going multiples are estimated in the data domain directly, for example using surface related multiple elimination, which removes surface related multiples without using additional information about the subsurface.
The up-going multiples (Receiver B in Table 1) can be removed from the recorded data to generate data representing the primaries. For example, up-going multiples can be adaptively removed from the recorded data (Receiver C in Table 1), which can consist of both primaries and multiples, to generate the data representing the primaries (Receiver A in Table 1).
A seismic source modeled as a point source (Source A in Table 1) can emit a down-going wavefield, which can be used as a source wavefield. Up-going primaries (Receiver A in Table 1), originating from surface related multiple elimination, can be used as a receiver wavefield. A seismic image of only primaries (Image 1 in Table 1) can be created using the down-going wavefield from a point source as a source wavefield and up-going primaries as a receiver wavefield. The migration input of up-going and down-going wavefields can be decomposed as represented by Equation 8 as follows:
Equations 8.1 and 8.2 can be substituted into Equation 2 to yield Equation 9 as follows:
I(x)=ΣωΣx
A seismic image from a primaries-only wavefield can be created using Equation 9. For example, Image 1 in Table 1 can be created using Equation 9.
Anti-causal crosstalk 552 is a false indication of an up-going wavefield being in phase with a down-going wavefield that corresponds to a location that is shallower and/or a time that is earlier than a location and/or time where and/or where two different wavefields are actually in phase. The anti-causal crosstalk 552 appears to indicate that the down-going multiple 544 is in phase with the second up-going multiple 546-2 at the false reflector 550. However, the second up-going multiple 546-2 is not a reflection of the down-going multiple 544. Instead, the second up-going multiple 546-2 is a reflection of the down-going primary 538 at the subsurface image point 536-2. The fact that the second up-going multiple 546-2 happens to have a point where it is in phase with the down-going multiple 544 can generate undesired seismic image data (anti-causal crosstalk 552) for a reflector that does not actually exist (the false reflector 550). This anti-causal crosstalk 552 can appear as noise in a seismic image. The anti-causal crosstalk 552 is said to correspond to a location that is shallower and/or a time that is earlier than a location and/or time where and/or where two different wavefields are actually in phase because it can be defined by its down-going component, which is the down-going multiple 544. The down-going component of the anti-causal crosstalk 552 is in phase with the fourth up-going multiple 546-4 at the third subsurface image point 536-3, which is deeper and later than the anti-causal crosstalk 552.
In at least one embodiment, for prediction of causal crosstalk, the seismic source 632 is modeled as a point source (Source A in Table 1) that emits the down-going point source wavefield 637. Up-going multiples, such as the up-going multiple 646, which can be received at a seismic receiver 642, can be used as a receiver wavefield (Receiver B in Table 1 by estimation of multiples). Locations where a down-going point source wavefield 637 is in phase with an up-going multiple 646 can constitute a prediction of causal crosstalk 649 (Image 2 in Table 1).
A migration input up-going wavefield and down-going wavefield can be decomposed as:
Equation (10.1) and (10.2) can be substituted into equation (2) to yield equation (11) as follows:
where equation (11.1) represents a first order causal crosstalk and equation (11.2) represents a second order causal crosstalk. In at least one embodiment, all orders of causal crosstalk (Image 2 in Table 1) are be predicted using only one migration.
In at least one embodiment, the causal crosstalk (Image 2 in Table 1) is computed by auto-convolution of the receiver wavefield (Receiver C in Table 1) using, for example, equation (4.1), which consists of both primaries and multiples at the sea surface as represented by equation (12) as follows:
where equation (12.1) represents a first order causal crosstalk, equation (12.2) represents a second order causal crosstalk, etc.
In the depth domain, the causal crosstalk (Image 2 in Table 1) can be adaptively removed from the seismic image of primaries with causal crosstalk (Image 3 in Table 1) to create a seismic image of only primaries (Image 1). A seismic image of only primaries can be equivalent to a seismic image created using the method described in conjunction with Equation 9.
At least one embodiment includes seismic imaging of all orders of recorded surface related multiples with causal crosstalk and anti-causal crosstalk, causal crosstalk prediction, anti-causal crosstalk prediction, both causal crosstalk and anti-causal crosstalk attenuation, and combination of a seismic image of primaries and a seismic image of multiples. Down-going primaries and multiples (Source D in Table 1) can be used as a source wavefield and the up-going primaries and multiples (Receiver C in Table 1) can be used as a receiver wavefield to create a seismic image from multiples (Image 12 in Table 1). A seismic image from multiples can include both causal crosstalk and anti-causal crosstalk. As an example, the up-going and down-going wavefields can be decomposed as represented by equation (13) as follows:
Equation 13 can be substituted in Equation 2 to yield Equation 14 as follows:
Different sources of crosstalk that can be present in a seismic image can be identified using Equation 14. Seismic imaging with multiples can be simplified since Equation 13 includes only one order of multiples in the source wavefield and two orders of multiples in the receiver wavefield.
Prediction of causal crosstalk can be performed as described above with respect to
Up-going and down-going wavefield migration inputs can be decomposed as represented by Equation 15 as follows:
Equation 15 can be substituted in Equation 2 to yield Equation 16 as follows:
Higher order multiples are not shown in Equation 16 for simplicity. However, at least one embodiment provides for prediction of all orders of anti-causal crosstalk noise (Image 10 in Table 1) using only one migration.
The causal and/or anti-causal crosstalk in the seismic image from multiples (Image 12 in Table 1) can be adaptively removed to create a crosstalk attenuated seismic image of multiples (Image 19 in Table 1). Crosstalk attenuated seismic images from primaries (Image 1 in Table 1) and from multiples (Image 19 in Table 1) can be combined to create a crosstalk attenuated seismic image from primaries and multiples (Image 20 in Table 1).
The number of engines can include a combination of hardware and program instructions that is configured to perform a number of functions described herein. The program instructions, such as software, firmware, etc., can be stored in a memory resource such as a machine-readable medium, etc., as well as hard-wired program such as logic. Hard-wired program instructions can be considered as both program instructions and hardware.
The seismic imaging engine 1376 can include a combination of hardware and program instructions that is configured to create a seismic image based on seismic data including multiples. The seismic image can include causal crosstalk and anti-causal crosstalk. The seismic imaging engine 1376 can be configured to create a seismic image that includes causal crosstalk and anti-causal crosstalk by migration of down-going primaries and multiples as a source wavefield of the seismic data from a seismic source to a subsurface image point and by migration of up-going primaries and multiples as a receiver wavefield of the seismic data from a seismic receiver to the subsurface image point. The seismic imaging engine 1376 can be configured to create the seismic image by applying an imaging condition at a subsurface image point. The seismic imaging engine 1376 can be configured to apply one of the group of seismic imaging conditions including deconvolution and cross-correlation at the subsurface image point. The object of deconvolution is to reverse the effects of convolution on recorded data, where convolution assumes that the recorded data is a combination of a reflectivity function and a source wavefield from a point source. Cross-correlation is a measure of similarity between two wavefields as a function of a time lag applied to one of them.
The prediction engine 1378 can include a combination of hardware and program instructions that is configured to predict causal crosstalk based on the seismic data. The prediction engine 1378 can be configured to predict anti-causal crosstalk based on the seismic data. The prediction engine 1378 can be configured to predict the anti-causal crosstalk by migration of down-going multiples as a source wavefield of the seismic data from a seismic source to a subsurface image point and by migration of up-going primaries as a receiver wavefield of the seismic data from a seismic receiver to the subsurface image point. The prediction engine 1378 can be configured to extract the down-going multiples using down-going and up-going wavefield separation of the seismic data and to estimate the up-going primaries using surface related multiple elimination. The prediction engine 1378 can be configured to predict the causal crosstalk by migration of a down-going wavefield as a source wavefield of the seismic data from a seismic source modeled as a point source to a subsurface image point and migration of up-going multiples as a receiver wavefield of the seismic data from a seismic receiver to the subsurface image point. The prediction engine can include a combination of hardware and program instructions that is configured to predict the causal crosstalk by auto-convolution of a receiver wavefield of the seismic data that includes both primaries and multiples. Auto-convolution of a wavefield is a convolution of the wavefield with itself. The prediction engine 1378 can be configured to extract the source wavefield and the receiver wavefield using down-going and up-going wavefield separation of the seismic data.
The attenuation engine 1380 can include a combination of hardware and program instructions that is configured to attenuate the predicted causal crosstalk and the predicted anti-causal crosstalk from the seismic image. The seismic imaging engine 1376 can be configured to combine the causal crosstalk and anti-causal crosstalk attenuated seismic image of the multiples with a crosstalk attenuated seismic image of primaries to obtain a crosstalk attenuated seismic image of the primaries and the multiples.
Memory resources 1486 can be non-transitory and can include volatile and/or non-volatile memory. Volatile memory can include memory that depends upon power to store information, such as various types of dynamic random access memory among others. Non-volatile memory can include memory that does not depend upon power to store information. Examples of non-volatile memory can include solid state media such as flash memory, electrically erasable programmable read-only memory, phase change random access memory, magnetic memory, optical memory, and/or a solid state drive, etc., as well as other types of non-transitory machine-readable media.
The processing resources 1484 can be coupled to the memory resources 1486 via a communication path 1488. The communication path 1488 can be local or remote to the machine 1482. Examples of a local communication path 1488 can include an electronic bus internal to a machine, where the memory resources 1486 are in communication with the processing resources 1484 via the electronic bus. Examples of such electronic buses can include Industry Standard Architecture, Peripheral Component Interconnect, Advanced Technology Attachment, Small Computer System Interface, Universal Serial Bus, among other types of electronic buses and variants thereof. The communication path 1488 can be such that the memory resources 1486 are remote from the processing resources 1484, such as in a network connection between the memory resources 1486 and the processing resources 1484. That is, the communication path 1488 can be a network connection. Examples of such a network connection can include a local area network, wide area network, personal area network, and the Internet, among others.
As shown in
Each of the number of modules 1490, 1492, 1494 can include program instructions and/or a combination of hardware and program instructions that, when executed by a processing resource 1484, can function as a corresponding engine as described with respect to
The machine 1482 can include a seismic imaging module 1490, which can include instructions to create a first seismic image based on first seismic data that includes primaries. The seismic imaging module 1490 can include instructions to create a second seismic image based on second seismic data that includes multiples.
The machine 1482 can include a prediction module 1492, which can include instructions to predict first causal crosstalk based on the first seismic data and predict second causal crosstalk and anti-causal crosstalk based on the second seismic data. The prediction module 1492 can include instructions to predict multiples in the first seismic data. The prediction module 1492 can include instructions to use down-going primaries and multiples as a source wavefield of the first seismic data and use the first seismic data as a subsurface reflection property to estimate the multiples in the first seismic data. The multiples can be an up-going wavefield. The prediction module 1492 can include instructions to predict multiples in the first seismic data via surface related multiple elimination.
The machine 1482 can include an attenuation module 1494, which can include instructions to attenuate the first causal crosstalk from the first seismic image and to attenuate the second causal crosstalk and the anti-causal crosstalk from the second seismic image. The attenuation module 1494 can include instructions to attenuate multiples from the first seismic data based on the predicted multiples.
The machine 1482 can include a seismic imaging module 1490, which can include instructions to combine the causal crosstalk attenuated first seismic image with the causal crosstalk and anti-causal crosstalk attenuated second seismic image.
In accordance with at least one embodiment of the present disclosure, a geophysical data product may be produced. The geophysical data product may include, for example, a causal crosstalk and anti-causal crosstalk attenuated seismic image. Geophysical data including seismic data may be obtained and stored on a non-transitory, tangible computer-readable medium. The geophysical data product may be produced by processing the geophysical data offshore or onshore either within the United States or in another country. If the geophysical data product is produced offshore or in another country, it may be imported onshore to a facility in the United States. In some instances, once onshore in the United States, geophysical analysis may be performed on the geophysical data product. In some instances, geophysical analysis may be performed on the geophysical data product offshore. For example, a seismic image can be created based on the seismic data including multiples, for example as it is being acquired or after it is acquired, offshore to facilitate other processing of the acquired seismic data either offshore or onshore. The seismic image can include causal crosstalk and anti-causal crosstalk. As another example, causal crosstalk and anti-causal crosstalk can be predicted based on the seismic data, for example as it is being acquired or after it is acquired, offshore to facilitate other processing of the acquired seismic data either offshore or onshore. As another example, the causal crosstalk and anti-causal crosstalk can be attenuated from the seismic image, for example as it is being acquired or after it is acquired, offshore to facilitate other processing of the acquired seismic data either offshore or onshore.
At block 1597, the method can include predicting causal crosstalk based on the seismic data. Predicting causal crosstalk can include migrating a down-going wavefield as a source wavefield of the seismic data from a seismic source modeled as a point source to a subsurface image point and migrating up-going multiples as a receiver wavefield of the seismic data from a seismic receiver to the subsurface image point.
At block 1598, the method can include predicting anti-causal crosstalk based on the seismic data. Predicting anti-causal crosstalk can include migrating down-going primaries and multiples as a source wavefield of the seismic data from a seismic source to a subsurface image point and migrating up-going primaries as a receiver wavefield of the seismic data from a seismic receiver to the subsurface image point.
At block 1599, the method can include attenuating the predicted causal crosstalk and the predicted anti-causal crosstalk from the seismic data.
The method can include joint seismic imaging of primaries of all orders of recorded multiples, causal crosstalk prediction, anti-causal crosstalk prediction, and both causal crosstalk and anti-causal crosstalk attenuation. Primaries and multiples can be jointly imaged by using a down-going wavefield, including primaries and multiples (Source E in Table 1), as a source wavefield from a point source, and an up-going wavefield including primaries and multiples (Receiver C) as a receiver wavefield. A seismic image from joint migration of primaries and multiples (Image 15 in Table 1) can be contaminated by causal and/or anti-causal crosstalk. The migration input up-going and down-going wavefields can be decomposed as represented by Equation 17 as follows:
Equation 17 can be substituted in Equation 2 to yield Equation 18 as follows:
where Equation 18.1 represents a seismic image of primaries, and Equations 18.2 and 18.3 represent seismic images of multiples. Equations 18.4-18.6 represent various orders of causal crosstalk, and Equations 18.7-18.9 represent various orders of anti-causal crosstalk.
Different sources of crosstalk present in the seismic image (Image 15 in Table 1) can be identified using Equation 18. Joint migration with primaries and multiples can be simplified since Equation 18 includes only one order of multiples in the source wavefield and two orders of multiples in the receiver wavefield. In Equation 18, the primaries and multiples can be included in a joint expression. When a direct arrival is not present, a point source can be scaled to make a seismic image of the primaries match a seismic image of the multiples. Different crosstalk terms can be predicted in the seismic image space so that the different crosstalk terms can be adaptively attenuated.
At least one embodiment of the present disclosure includes a seismic imaging method. For example, a first method can include seismic imaging of multiples with both causal crosstalk and anti-causal crosstalk, causal crosstalk prediction, anti-causal crosstalk prediction, and both causal crosstalk and anti-causal crosstalk attenuation. A second method can include joint seismic imaging of primaries and multiples, causal crosstalk prediction, and both causal crosstalk and anti-causal crosstalk attenuation. A third method can include seismic imaging of primaries with causal crosstalk, causal crosstalk prediction, causal crosstalk attenuation, and a combination of crosstalk-free seismic image of primaries with a crosstalk-free seismic image of multiples (e.g. from the first method). Embodiments are not constrained to a seismic imaging method, however, and one-way wave equation migration and/or two-way reverse time migration wave equation migration can be used. At least one embodiment can be applied in a post-stack image domain and/or a pre-stack image domain (e.g., subsurface offset and/or angle gathers).
Although specific embodiments have been described above, these embodiments are not intended to limit the scope of the present disclosure, even where only a single embodiment is described with respect to a particular feature. Examples of features provided in the disclosure are intended to be illustrative rather than restrictive unless stated otherwise. The above description is intended to cover such alternatives, modifications, and equivalents as would be apparent to a person skilled in the art having the benefit of this disclosure.
The scope of the present disclosure includes any feature or combination of features disclosed herein (either explicitly or implicitly), or any generalization thereof, whether or not it mitigates any or all of the problems addressed herein. Various advantages of the present disclosure have been described herein, but embodiments may provide some, all, or none of such advantages, or may provide other advantages.
In the foregoing Detailed Description, some features are grouped together in a single embodiment for the purpose of streamlining the disclosure. This method of disclosure is not to be interpreted as reflecting an intention that the disclosed embodiments of the present disclosure have to use more features than are expressly recited in each claim. Rather, as the following claims reflect, inventive subject matter lies in less than all features of a single disclosed embodiment. Thus, the following claims are hereby incorporated into the Detailed Description, with each claim standing on its own as a separate embodiment.
Whitmore, Jr., Norman Daniel, Lu, Shaoping, Valenciano Mavilio, Alejandro Antonio
Patent | Priority | Assignee | Title |
11892583, | Jul 10 2019 | Abu Dhabi National Oil Company | Onshore separated wave-field imaging |
Patent | Priority | Assignee | Title |
5608814, | Aug 26 1993 | General Electric Company | Method of dynamic thresholding for flaw detection in ultrasonic C-scan images |
7039197, | Oct 19 2000 | Lear Corporation | User interface for communication system |
7117145, | Oct 19 2000 | Lear Corporation | Adaptive filter for speech enhancement in a noisy environment |
7171003, | Oct 19 2000 | Lear Corporation | Robust and reliable acoustic echo and noise cancellation system for cabin communication |
7995829, | Aug 01 2007 | General Electric Company | Method and apparatus for inspecting components |
20040081322, | |||
20060227902, | |||
20100030830, | |||
20120155217, | |||
20120281746, | |||
20130054664, | |||
20150046093, |
Executed on | Assignor | Assignee | Conveyance | Frame | Reel | Doc |
Jan 08 2016 | PGS Geophysical AS | (assignment on the face of the patent) | / | |||
Jan 11 2016 | LU, SHAOPING | PGS Geophysical AS | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 037462 | /0719 | |
Jan 11 2016 | WHITMORE, NORMAN DANIEL, JR | PGS Geophysical AS | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 037462 | /0719 | |
Jan 11 2016 | VALENCIANO MAVILIO, ALEJANDRO ANTONIO | PGS Geophysical AS | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 037462 | /0719 |
Date | Maintenance Fee Events |
Mar 03 2023 | M1551: Payment of Maintenance Fee, 4th Year, Large Entity. |
Date | Maintenance Schedule |
Sep 03 2022 | 4 years fee payment window open |
Mar 03 2023 | 6 months grace period start (w surcharge) |
Sep 03 2023 | patent expiry (for year 4) |
Sep 03 2025 | 2 years to revive unintentionally abandoned end. (for year 4) |
Sep 03 2026 | 8 years fee payment window open |
Mar 03 2027 | 6 months grace period start (w surcharge) |
Sep 03 2027 | patent expiry (for year 8) |
Sep 03 2029 | 2 years to revive unintentionally abandoned end. (for year 8) |
Sep 03 2030 | 12 years fee payment window open |
Mar 03 2031 | 6 months grace period start (w surcharge) |
Sep 03 2031 | patent expiry (for year 12) |
Sep 03 2033 | 2 years to revive unintentionally abandoned end. (for year 12) |