US20050018538A1 - Method of seismic processing for the decomposition of a wavefield into harmonic components and applications to the determination of angular gathers of reflectivity - Google Patents
Method of seismic processing for the decomposition of a wavefield into harmonic components and applications to the determination of angular gathers of reflectivity Download PDFInfo
- Publication number
- US20050018538A1 US20050018538A1 US10/855,930 US85593004A US2005018538A1 US 20050018538 A1 US20050018538 A1 US 20050018538A1 US 85593004 A US85593004 A US 85593004A US 2005018538 A1 US2005018538 A1 US 2005018538A1
- Authority
- US
- United States
- Prior art keywords
- max
- local
- component
- wavefield
- harmonic
- 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.)
- Abandoned
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000000354 decomposition reaction Methods 0.000 title claims description 34
- 238000002310 reflectometry Methods 0.000 title claims description 22
- 238000001914 filtration Methods 0.000 claims abstract description 22
- 230000000306 recurrent effect Effects 0.000 claims abstract description 3
- 230000001419 dependent effect Effects 0.000 claims description 23
- 238000001228 spectrum Methods 0.000 claims description 13
- 230000002194 synthesizing effect Effects 0.000 claims description 13
- 230000005012 migration Effects 0.000 description 9
- 238000013508 migration Methods 0.000 description 9
- 230000009466 transformation Effects 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
Definitions
- the present invention relates to a method of seismic processing for the decomposition of a wavefield into various harmonic components.
- Prestack depth migration is a central step in seismic processing. It consists in focusing the seismic events recorded time-wise, as reflections indexed depth-wise.
- the output of the migration is a depth-dependent reflectivity computed on the basis of the cross-correlation of the incident and reflected waves and normalized on the basis: of the autocorrelation of the incident wave i(x, z, f, m) representing the incident wave and s(x, z, f, m) the reflected wave, for a given frequency f and a given shot m.
- AVA amplitude versus angle
- Processing can also be carried out to attenuate the multiple arrivals by rejecting events exhibiting a curvature in the depth-angle plane of the angle gather.
- This procedure is expensive in terms of computation time being based as it is on a local transformation performed at each point and on the computation of a matrix.
- the invention proposes a seismic processing technique making it possible to perform the decomposition of a wavefield into various harmonic components.
- the invention proposes a method of processing seismic data representative of at least one wavefield, characterized in that local harmonic components of the wavefield are determined at a given point P of the subsurface by implementing a recurrent processing according to which a harmonic of rank n+1, n being a positive integer, is determined as a function of one or more harmonics of rank n or lower to which is applied at least one filtering which depends on at least one local parameter at the point P, this local parameter being chosen from among the components of the wavevector, the frequency of the wavefield, the local velocity, the local anisotropy parameters or any combination of these various parameters.
- a harmonic of rank n ⁇ 1, n being a negative integer is determined as a function of one or more harmonics of rank n or higher to which is applied at least one filtering which depends on at least one local parameter of the subsurface, this local parameter being chosen from among the components of the wavevector, the frequency of the wavefield, the local velocity, the local anisotropy parameters or any combination of these various parameters.
- the B n ( ⁇ ) may not be orthogonal over the interval ( ⁇ min , ⁇ max ) considered, and in this case it is necessary to compute the matrix: r mn ⁇ ⁇ ⁇ min ⁇ max ⁇ B m ⁇ ( ⁇ ) _ ⁇ ⁇ B n ⁇ ( ⁇ ) ⁇ d ⁇
- the components w n (x) are obtained by inverting a linear system whose matrix is ⁇ r mn ⁇
- the filters involved in a recurrence relation are filters dependent on the physics of propagation in the subsurface and in particular on one or more parameters chosen from among the components of the. wavevector, the frequency of the wavefield, the local velocity, the local anisotropy parameters or any combination of these various parameters.
- the invention proposes a method of processing seismic data according to which incident waves and reflected waves are processed at various points on the surface of the ground so as to determine an angle reflectivity gather, characterized by the steps according to which:
- FIG. 1 illustrates an exemplary implementation of processing for determining harmonics
- FIG. 2 illustrates various steps of an exemplary computation of angular gathers
- FIG. 3 illustrates a velocity model used for a data set
- FIG. 4 illustrates migration with the exact velocity model
- FIG. 5 illustrates the angle gather obtained with the exact velocity model
- FIG. 6 illustrates the angle gather which is obtained with too slow a velocity.
- w(x) be a seismic wavefield at a frequency ⁇ and a depth z.
- this wavefield w(x) is the sum of its downgoing component w + (x) and of its upgoing component w ⁇ (x).
- ⁇ is defined over [ ⁇ /2, ⁇ /2] for the downgoing wave (cos ⁇ 0), and over [ ⁇ /2,3 ⁇ /2] for the upgoing wave (cos ⁇ 0).
- the procedure described also differs from that of Xie and Wu because the computation of the angle transform does not require any local transform.
- the above derivation shows that the local harmonic decomposition of a wavefield can be computed directly through a recursive filtering of the data by the local velocity-dependent filter given by equation (14).
- This derivation may be extended to a velocity field with lateral variation by using a so-called explicit filtering structure.
- Each of these filters must be accurate for the values
- the modulus of these filters must also be less than 1 everywhere for stability reasons.
- These filters can be non-recursive finite impulse response (FIR filters) or be infinite response filters applied in a recursive way (IIR filters).
- step 1 of FIG. 1 the recurrence processing for the positive integers (step 2 ) and the recurrence processing for the negative integers (step 3 ) are implemented in parallel.
- the harmonic decomposition method described above makes it possible to compute in an advantageous manner various angle gathers in the course of a shotpoint-based migration.
- a frequency ⁇ and a depth z the incident wavefield is i(x) and the reflected wavefield is s(x).
- r n (x,z, ⁇ ,m) ( - 1 ) n ⁇ ⁇ ⁇ f , m ⁇ i ⁇ n ⁇ ( x , z , f , m ) _ ⁇ ⁇ s ⁇ n ⁇ ( x , z , f , m )
- the harmonics r n (x,z) are in themselves data that can be used in the seismic processing.
- a criterion of exactness of the velocity model is thus available: in the case where the velocity model used to migrate the data is exact, the (x,z)-wise images provided by the various harmonics should be superimposable.
- FIG. 4 shows the migration with the exact velocity model.
- FIG. 5 shows the angle gather for the horizontal coordinate 5 km with the exact velocity model. The exactness of the velocity model can be verified over all the events.
- FIG. 6 shows; the same angle gather when the migration is made with too slow a velocity by 4%. The curvature of the events confirms that the velocity is too slow. All the events in the angular gather are centred on 0 degrees and independently of their own dip, as is desirable. The absence of artefacts can also be verified.
- the angle gather described previously is the gather in terms of opening angle.
- Other gathers may be computed on the basis of local harmonic decompositions of the incident and/or reflected wavefield.
- the angle gathers computed on the basis of the cross-correlations are kinematic gathers where the amplitude of the events is not preserved. It is possible to obtain gathers with preserved amplitudes by computing normalization factors on the basis of various autocorrelations of the incident field.
- r ⁇ n ⁇ ( x , z ) ( - 1 ) n ⁇ ⁇ f , m ⁇ i ⁇ n ⁇ ( x , z , f , m ) _ ⁇ ⁇ s ⁇ n ⁇ ( x , z , f , m ) ⁇ f , m ⁇ i ⁇ n ⁇ ( x , z , f , m ) _ ⁇ ⁇ i ⁇ n ⁇ ( x , z , f , m )
- the method described can be implemented with variants.
- ⁇ arcsin ⁇ ⁇ ck x ⁇ ( 23 ) which has the advantage of not depending on the dispersion relation and hence of being valid even in the presence of anisotropy.
- the filter no longer depends on ⁇ /c. It should on the other hand be synthesized as a two-dimensional filter h(x,z). Since the filter does not depend on the frequency, it may be applied after summation over all frequencies and all shots.
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Image Processing (AREA)
Abstract
Description
- The present invention relates to a method of seismic processing for the decomposition of a wavefield into various harmonic components.
- It also proposes a method for the determination of angle reflectivity gathers.
- 1. General Field
- Prestack depth migration is a central step in seismic processing. It consists in focusing the seismic events recorded time-wise, as reflections indexed depth-wise. One of the most accurate ways of performing this step is migration by shotpoint, which consists in numerically propagating an incident wave representing the seismic source and a reflected wave. The incident wave is initialized surface-wise at the depth z=0 with a synthetic representation of the source and the reflected wave is so initialized with the wave recorded by the seismic sensors. Numerical propagation propagates these waves gradually across layers of thickness Δz. The output of the migration is a depth-dependent reflectivity computed on the basis of the cross-correlation of the incident and reflected waves and normalized on the basis: of the autocorrelation of the incident wave i(x, z, f, m) representing the incident wave and s(x, z, f, m) the reflected wave, for a given frequency f and a given shot m.
- This reflectivity is a scalar which satisfies the condition:
and is obtained by summing over all the frequencies and all the shots the cross-correlation of the incident wave and the reflected wave. - The representation of reflectivity by a scalar alone is however a drawback. It turns out in fact to be desirable to have a reflectivity dependent on the opening angle (difference between the angle of incidence and the angle of reflection). This angle-dependent reflectivity makes it possible to compute, for each horizontal position x, an angle gather which takes the form of an amplitude dependent on the depth z and angle θ variables.
- This allows on the one hand AVA (amplitude versus angle) analyses, in which the amplitude variation of the reflection coefficient as a function of the opening angle makes it possible to characterize the local medium; this also makes it possible on the other hand to carry out velocity analyses by verifying that the seismic arrivals do not have angle-dependent arrival times. Processing can also be carried out to attenuate the multiple arrivals by rejecting events exhibiting a curvature in the depth-angle plane of the angle gather.
- 1. Current State of the Art
- Rickett and Sava (2001) have set forth a procedure for computing angle gathers. This procedure computes the cross-correlations for horizontal relative displacements between the incident and reflected-wave:
- After Fourier transform on the z and h variables, we obtain R(x, kz, kh), then:
r(x, k z, θ)=R(x, k z , k z tan θ)
the angle gather r(x,z,θ) being obtained by inverse Fourier transform on the variable z. According to the conclusions of the authors, this procedure would exhibit troublesome artefacts in the case where the distance between the summed sources is not small, this being the case in practice. - Xie and Wu (2002) describe a procedure that does not seem to produce artefacts. This procedure is based on the transformation:
where ƒ(u) is an apodization window. Once this transformation has been performed on i(x,z,ƒ,m) and s(x,z,ƒ,m), we obtain: - This procedure is expensive in terms of computation time being based as it is on a local transformation performed at each point and on the computation of a matrix.
- The invention proposes a seismic processing technique making it possible to perform the decomposition of a wavefield into various harmonic components.
- It also proposes a technique allowing the determination of angle gathers of reflectivity coefficients.
- More precisely, the invention proposes a method of processing seismic data representative of at least one wavefield, characterized in that local harmonic components of the wavefield are determined at a given point P of the subsurface by implementing a recurrent processing according to which a harmonic of rank n+1, n being a positive integer, is determined as a function of one or more harmonics of rank n or lower to which is applied at least one filtering which depends on at least one local parameter at the point P, this local parameter being chosen from among the components of the wavevector, the frequency of the wavefield, the local velocity, the local anisotropy parameters or any combination of these various parameters.
- Advantageously also, a harmonic of rank n−1, n being a negative integer, is determined as a function of one or more harmonics of rank n or higher to which is applied at least one filtering which depends on at least one local parameter of the subsurface, this local parameter being chosen from among the components of the wavevector, the frequency of the wavefield, the local velocity, the local anisotropy parameters or any combination of these various parameters.
- The invention also proposes a method of the aforesaid type in which the harmonics are determined through the recurrence:
ŵ n+1(P)=h k0 (P) +(P)*ŵ n(P)
for n a positive integer and
ŵ n−1(P)=h k0 (P) −(P)*ŵ n(P)
for n a negative integer - where the operation * designates a filtering, or any linear combination implementing the points neighbouring the point P, and
- where h+ k0(P) is a filter corresponding to a filtering as a function of the wavenumber vector k such that:
H k0 +(k)=e −jθ(k)
and - where h− k0 (P) is a filter corresponding to a filtering as a function of the wavenumber vector k such that:
H k0 −( k)={overscore (H k0 + (k))}= e jθ(k)
θ(k) being an angular parameter dependent on at least one component of the wavenumber vector and on k0=ω/c with (ω=2πƒ, f designating the frequency of the wavefield, c its local velocity, - the initial component ŵ0(P)) being equal to the transform in the frequency domain of the wavefield at the point P or being a function thereof.
- The invention furthermore proposes a method of processing seismic data according to which angular components of a wavefield w(P), where P designates a given point of the subsurface, are determined in the frequency domain, characterized in that harmonic components wn(P) of this wavefield are determined by implementing a method according to one of the preceding claims and its angular components are reconstructed as a function of an angle θ which depends physically on a local parameter chosen from among the components of the wavevector, the frequency of the wavefield, the local velocity, the local anisotropy parameters or any combination of these various parameters,
by determining - The method can be generalized by decomposing the angular components W(x,θ) into:
where the function Bn(θ) satisfies a recurrence relation, which may be of order higher than 1, for example:
B n+1(θ)=F(θ)B n(θ)+G(θ)B n−1(θ) - The Bn(θ) may not be orthogonal over the interval (θmin, θmax) considered, and in this case it is necessary to compute the matrix:
- The components wn(x) are obtained by inverting a linear system whose matrix is {rmn}
- and whose right-hand side is formed of components wn(P) obtained through a recurrence relation:
ŵ n+1(P)=h k0 (P)(P)*ŵ n(P)+q k0 (P)(P)*ŵ n−1(P)
where the filters h and q are expressed in terms of wavenumber k with the aid of the relation θk0(k) expressing the angle as a function of wavenumbers and of the parameter k0=ω/c(P) where c(P) is the local velocity:
H k0 (k)=F(θk0(k)), Q(k)=G(θk0(k)) - More generally, the filters involved in a recurrence relation are filters dependent on the physics of propagation in the subsurface and in particular on one or more parameters chosen from among the components of the. wavevector, the frequency of the wavefield, the local velocity, the local anisotropy parameters or any combination of these various parameters.
- Also, the invention proposes a method of processing seismic data according to which incident waves and reflected waves are processed at various points on the surface of the ground so as to determine an angle reflectivity gather, characterized by the steps according to which:
-
- for the various points considered, as well as for the various frequencies and the various shots, a local harmonic decomposition of the incident wave and/or a local harmonic decomposition of the reflected wave is/are computed,
- a local harmonic decomposition of the reflectivity is deduced therefrom for the various points considered and is summed frequency-wise and shot-wise,
- an angle gather of the reflectivity is deduced from the harmonic decomposition thus obtained after summation.
- The description which follows is purely illustrative and nonlimiting and should be read in conjunction with the appended drawings in which:
-
FIG. 1 illustrates an exemplary implementation of processing for determining harmonics; -
FIG. 2 illustrates various steps of an exemplary computation of angular gathers; -
FIG. 3 illustrates a velocity model used for a data set; -
FIG. 4 illustrates migration with the exact velocity model; -
FIG. 5 illustrates the angle gather obtained with the exact velocity model; -
FIG. 6 illustrates the angle gather which is obtained with too slow a velocity. - Theoretical Elements
- Let w(x) be a seismic wavefield at a frequency ƒ and a depth z. We shall use the wavenumber ω=2πƒ instead of the frequency ƒ. We wish to decompose this wavefield as the sum of angle components. In general, this wavefield w(x) is the sum of its downgoing component w+(x) and of its upgoing component w−(x). We define the angle θ such that θ=0 corresponds to a downward vertical propagation. θ is defined over [−π/2,π/2] for the downgoing wave (cosθ≧0), and over [π/2,3π/2] for the upgoing wave (cosθ<0).
- Let us consider the downgoing wave to begin with. We can take the Fourier transform:
where c is the local velocity. The evanescent energy corresponds to |k|>ω/c. We define the angle through the change of variable: - We define an angle transform W+(θ) as being zero over θ∈[π/2,3π/2] and defined over θ∈[−π/2,/2,π/2] by:
- With this definition, w+(x) can be decomposed as:
- We define a local angle transform by changing equation (4) into:
- This equation differs from that of Xie and Wu (2002) through the cosine term. This term is important since it ensures that the average of the angular components W+(x,θ) gives the field w+(x).
- The procedure described also differs from that of Xie and Wu because the computation of the angle transform does not require any local transform.
- Accordingly, starting from equation (6) we can, for each x, perform the Fourier transform on the angular variable. Since the latter is periodic with period 2π, this corresponds to a decomposition into harmonics:
- It will be noted that a characteristic of the technique proposed is that the harmonics wn +(x) may be computed directly from the field w+(x), without any intermediate computation of W+(x,θ). Since W+(x,θ) is zero over [π/2,3π/2]:
- Inserting equation (6):
exchanging the order of the integrals and making the change of variable k=ω/c sinθ, we obtain: - Let us define in the k domain the filters:
we obtain, ƒn +(x) being the corresponding filters in the x domain: - We see that we can compute wn +(x) by filtering w+(x) with the filters ƒn +(x), which, by using equation (2), may be written symbolically Fn +(k)=e−jnθ.
- These filters may be defined through the recursion:
F 0 +(k)=1
F n+1 +(k)=H k0 +(k)F n +(k), n≧0
F n−1 +(k)={overscore (H k0 + (k))} F n +(k), n≦0 (14)
where the filter Hk0 +(k)=e−jθ, which depends on k0=ω/c, is:
H k0 +(k)={square root}{square root over (1−(k/k 0)2)}−jk/k 0 (15) - The local harmonic decomposition may therefore advantageously be computed through:
w 0 +(x)=w +(x)
w n+1 +(x)=h k0 (x)(x)*w n +(x), n≧0
w n−1 +(x)=h k0 (x) +(−x)*w n +(x), n≦0 (16) - This can be generalized to a wave comprising an upgoing part, w(x)=w+(x)+w−(x), the local harmonic decomposition becoming:
where the ƒn −(x) are defined as previously except for the sign of the square root, which is the sign of cosθ, and which is negative instead of positive. We therefore have:
H{overscore (k)} 0(k)=−{square root}{square root over (1−(k/k 0)2)}−jk/k 0 =−{overscore (H k0 + (k))} - We can compute the local harmonic decomposition of the angle transformation:
through a recursion of the same type as (16) but where the initialization is replaced by:
w 0 +(x)=g k0 (x)(x)*w +(x)
where gk0(x) is a spatial filter synthesizing the k-spectrum:
Exemplary Implementation of the Invention - The above derivation shows that the local harmonic decomposition of a wavefield can be computed directly through a recursive filtering of the data by the local velocity-dependent filter given by equation (14). This derivation may be extended to a velocity field with lateral variation by using a so-called explicit filtering structure. In such a structure, the filters hk0(x) are tabulated beforehand for the useful values of k0, that is to say the interval [ωmin/cmax,ωmax/cmin], since at the point x, we use the coefficients of the filter corresponding to k0=ω/c(x). Each of these filters must be accurate for the values |k|≦k0 sinθmax, where θmax is the maximum propagation angle and |k|≦αkNyq, where α is the proportion of the spatial Nyquist band preserved and kNyq the Nyquist spatial frequency. The modulus of these filters must also be less than 1 everywhere for stability reasons. These filters can be non-recursive finite impulse response (FIR filters) or be infinite response filters applied in a recursive way (IIR filters).
- The local harmonic decomposition is therefore computed through an explicit recursive filtering structure, of the type wn+1 +(x)=hk0(x)(x)*wn +(x), initialized with w0 +(x)=w+(x).
- As illustrated in
FIG. 1 , the initial component may also be:
w 0(P)=g k0 (P)(P)*w(P)
where gk0(P) is a spatial filter synthesizing the k-spectrum:
k here designating the wavevector, or a component of the latter or else a parameter dependent on the latter. - Once this initial component has been determined through a first explicit filtering step (step 1 of
FIG. 1 ), the recurrence processing for the positive integers (step 2) and the recurrence processing for the negative integers (step 3) are implemented in parallel. - As output we obtain the various harmonics wN
- Application to the Computation of Angular Gathers in Shotpoint-Based Migration
- The harmonic decomposition method described above makes it possible to compute in an advantageous manner various angle gathers in the course of a shotpoint-based migration. For a given shot s, a frequency ƒ and a depth z, the incident wavefield is i(x) and the reflected wavefield is s(x). We shall consider the local angle transformations I(x,θ) and S(x,θ) of i(x) and s(x). The angle-wise reflectivity R(x,θ) connects I(x,θ) to S(x,θ) by:
by putting G(x,Δθ)=R(x,π−Δθ). - Equation (18) is, for any x, an angle-wise convolution. If we perform a Fourier transform on the angular variable, and if we introduce the local harmonic decomposition, it simplifies into:
ŝ n(x)=ĝ n(x)î n(x) (19) - If we consider only the kinematics, an angle gather may be computed by:
ĝ n(x)={overscore (î)} n(x)ŝ n(x) (20) - Hence the local harmonic decomposition of the reflectivity can be computed through the relation:
{circumflex over (r)} n(x))=(−1)n {overscore (+E,cir )} n(x)ŝ n(x) (21)
where in(x) and sn(x) are the local harmonic decompositions of i(x) and s(x). - Thus, as illustrated in
FIG. 2 , after computing the local harmonics of the incident (step 4) and reflected (step 5) waves, we compute the various harmonics of. reflectivity rn(x,z,ƒ,m) for n=[−N,N] for each depth z, frequency ƒ and shot m (step 6) and we sum them over all the frequencies and shots (step 7): - The angle gather r(x,z,θ) can then be reconstructed (step 8) from the harmonics rn(x,z) by:
- However, the harmonics rn(x,z) are in themselves data that can be used in the seismic processing. In particular, a criterion of exactness of the velocity model is thus available: in the case where the velocity model used to migrate the data is exact, the (x,z)-wise images provided by the various harmonics should be superimposable.
- Equation (21) may be considered to be a harmonic imaging condition generalizing the conventional scalar imaging condition
r(x)={overscore (i(x))} s(x) (22) - An important property of the angle gather r(x,z,θ) obtained by this procedure is that its θ-wise average gives the conventional imaging condition {overscore (i(x))}s(x). Specifically, the average is the zero harmonic and {circumflex over (r)}0(x)={overscore (î)}0(x)ŝ0(x)=i(x)s (x).
- Comments Relating to the Illustration Given by FIGS. 3 et seq.
- The figures are 2D synthetic data sets illustrating implementations. 100 shots spaced 100 m apart were generated on the basis of linear reflectors of varied dips and of the velocity model of
FIG. 3 .FIG. 4 shows the migration with the exact velocity model.FIG. 5 shows the angle gather for the horizontal coordinate 5 km with the exact velocity model. The exactness of the velocity model can be verified over all the events.FIG. 6 shows; the same angle gather when the migration is made with too slow a velocity by 4%. The curvature of the events confirms that the velocity is too slow. All the events in the angular gather are centred on 0 degrees and independently of their own dip, as is desirable. The absence of artefacts can also be verified. - Other Types of Angle Gathers
- The angle gather described previously is the gather in terms of opening angle. Other gathers may be computed on the basis of local harmonic decompositions of the incident and/or reflected wavefield.
- The gather in terms of opening angle is obtained through:
- An angle gather which indicates the local dip of the reflector is obtained through:
- The gather in terms of angle of incidence is obtained by taking only the local harmonic decomposition of the incident field alone:
Amplitude Preserved: - The angle gathers computed on the basis of the cross-correlations are kinematic gathers where the amplitude of the events is not preserved. It is possible to obtain gathers with preserved amplitudes by computing normalization factors on the basis of various autocorrelations of the incident field. For example, for the gather in terms of opening angle, it is possible to compute the local harmonic decomposition of the reflectivity through:
Other Variant Implementations - The method described can be implemented with variants. The method is based on the synthesis of a filter effecting the transfer function H=e−jθ. The implementation described previously is based on the relation:
which has the advantage of not depending on the dispersion relation and hence of being valid even in the presence of anisotropy. However, as a general rule, we have to compute the decomposition into wavefield angular components which obey a dispersion relation. This implies that other expressions for θ are possible. In the case of an isotropic medium,
which implies that we can write: - The filter can therefore be written:
- The advantage of this expression is that the filter no longer depends on ω/c. It should on the other hand be synthesized as a two-dimensional filter h(x,z). Since the filter does not depend on the frequency, it may be applied after summation over all frequencies and all shots.
- References
-
- Rickett, J. and Sava, P., 2001, Offset and angle domain common image gathers for shot profile migration: 71st Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 1115-1118.
- Xie, X. B. and Wu, R. S., 2002, Extracting angle domain information from migrated wavefield: 72nd Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 1360-1363.
Claims (35)
ŵ n+1(P)=h k
ŵ n−1(P)=h k
H k
H k
H k
w 0(P)=g k
ŵ n+1(P)=hk
ŵ n−1(P)=h k
H k
H k
H k
w 0(P)=g k
w 0(P)=g k
w 0(P)=g k
w 0(P)=g k
w 0(P)=g k
w 0(P)=g k
w 0(P)=g k
w 0(P)=g k
w 0(P)=g k
w 0(P)=g k
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/881,130 US7889598B2 (en) | 2003-05-27 | 2007-07-25 | Method of seismic processing for the decomposition of a wavefield into harmonic components and applications to the determination of angular gathers of reflectivity |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0306451 | 2003-05-27 | ||
FR0306451A FR2855618B1 (en) | 2003-05-27 | 2003-05-27 | SEISMIC PROCESSING METHOD FOR THE DECOMPOSITION OF A WAVE FIELD OF HARMONIC COMPONENTS AND APPLICATIONS FOR THE DETERMINATION OF ANGULAR REFLECTIVITY COLLECTIONS |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/881,130 Division US7889598B2 (en) | 2003-05-27 | 2007-07-25 | Method of seismic processing for the decomposition of a wavefield into harmonic components and applications to the determination of angular gathers of reflectivity |
Publications (1)
Publication Number | Publication Date |
---|---|
US20050018538A1 true US20050018538A1 (en) | 2005-01-27 |
Family
ID=32669394
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/855,930 Abandoned US20050018538A1 (en) | 2003-05-27 | 2004-05-27 | Method of seismic processing for the decomposition of a wavefield into harmonic components and applications to the determination of angular gathers of reflectivity |
US11/881,130 Expired - Fee Related US7889598B2 (en) | 2003-05-27 | 2007-07-25 | Method of seismic processing for the decomposition of a wavefield into harmonic components and applications to the determination of angular gathers of reflectivity |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/881,130 Expired - Fee Related US7889598B2 (en) | 2003-05-27 | 2007-07-25 | Method of seismic processing for the decomposition of a wavefield into harmonic components and applications to the determination of angular gathers of reflectivity |
Country Status (7)
Country | Link |
---|---|
US (2) | US20050018538A1 (en) |
CA (2) | CA2831246A1 (en) |
FR (1) | FR2855618B1 (en) |
GB (1) | GB2402216B (en) |
MX (1) | MXPA05012828A (en) |
NO (1) | NO335231B1 (en) |
WO (1) | WO2004106975A1 (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070258323A1 (en) * | 2006-05-07 | 2007-11-08 | Ping An | Method for Seismic Trace Decomposition and Reconstruction Using Multiple Wavelets |
US20080285383A1 (en) * | 2006-05-07 | 2008-11-20 | Ping An | System and method for processing seismic data for interpretation |
US8797828B1 (en) * | 2010-02-15 | 2014-08-05 | Soreq NRC | Remote optical seismic surveying and detection and imaging of underground objects |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108957544B (en) * | 2018-08-16 | 2019-07-26 | 中国石油大学(北京) | Measurement method, device, seismometer and the medium of near surface anisotropic parameters |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2696891A (en) * | 1951-11-26 | 1954-12-14 | Neufeld Jacob | Seismic surveying |
US4839527A (en) * | 1986-10-28 | 1989-06-13 | Alan Leitch | Optical-fibre smoke detection/analysis system |
US5627515A (en) * | 1995-02-24 | 1997-05-06 | Pittway Corporation | Alarm system with multiple cooperating sensors |
US5705979A (en) * | 1995-04-13 | 1998-01-06 | Tropaion Inc. | Smoke detector/alarm panel interface unit |
US5852588A (en) * | 1994-07-07 | 1998-12-22 | Schlumberger Technology Corporation | Method of processing seismic data |
US5896088A (en) * | 1997-04-16 | 1999-04-20 | Southeastern Univ. Research Assn. | Incipient fire detection system |
US5933790A (en) * | 1994-12-23 | 1999-08-03 | Schlumberger Technology Corporation | Data compression for seismic signal data |
US5971095A (en) * | 1996-01-09 | 1999-10-26 | Schlumberger Technology Corporation | Noise filtering method for seismic data |
US20020010545A1 (en) * | 2000-03-01 | 2002-01-24 | Sitton Gary Arthur | Method for vibrator sweep analysis and synthesis |
US6791453B1 (en) * | 2000-08-11 | 2004-09-14 | Walter Kidde Portable Equipment, Inc. | Communication protocol for interconnected hazardous condition detectors, and system employing same |
US6829538B2 (en) * | 2002-04-23 | 2004-12-07 | Westerngeco, L.L.C. | Method for suppressing systematic amplitude variations |
US6834235B2 (en) * | 2000-12-07 | 2004-12-21 | Westerngeco, L.L.C. | Method of and apparatus for processing seismic data |
US6882272B2 (en) * | 2001-06-02 | 2005-04-19 | Robert Bosch Gmbh | Danger detecting system |
US20060074562A1 (en) * | 2002-05-28 | 2006-04-06 | Lasse Amundsen | Processing seismic data |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4780859A (en) * | 1987-03-09 | 1988-10-25 | Mobil Oil Corporation | Method of interpreting seismic data |
US6131071A (en) * | 1996-12-06 | 2000-10-10 | Bp Amoco Corporation | Spectral decomposition for seismic interpretation |
-
2003
- 2003-05-27 FR FR0306451A patent/FR2855618B1/en not_active Expired - Fee Related
-
2004
- 2004-05-25 NO NO20042162A patent/NO335231B1/en not_active IP Right Cessation
- 2004-05-25 GB GB0411672A patent/GB2402216B/en not_active Expired - Fee Related
- 2004-05-27 US US10/855,930 patent/US20050018538A1/en not_active Abandoned
- 2004-05-27 CA CA2831246A patent/CA2831246A1/en not_active Abandoned
- 2004-05-27 WO PCT/IB2004/002003 patent/WO2004106975A1/en active Application Filing
- 2004-05-27 CA CA002527215A patent/CA2527215A1/en not_active Abandoned
- 2004-05-27 MX MXPA05012828A patent/MXPA05012828A/en active IP Right Grant
-
2007
- 2007-07-25 US US11/881,130 patent/US7889598B2/en not_active Expired - Fee Related
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2696891A (en) * | 1951-11-26 | 1954-12-14 | Neufeld Jacob | Seismic surveying |
US4839527A (en) * | 1986-10-28 | 1989-06-13 | Alan Leitch | Optical-fibre smoke detection/analysis system |
US5852588A (en) * | 1994-07-07 | 1998-12-22 | Schlumberger Technology Corporation | Method of processing seismic data |
US5933790A (en) * | 1994-12-23 | 1999-08-03 | Schlumberger Technology Corporation | Data compression for seismic signal data |
US5627515A (en) * | 1995-02-24 | 1997-05-06 | Pittway Corporation | Alarm system with multiple cooperating sensors |
US5705979A (en) * | 1995-04-13 | 1998-01-06 | Tropaion Inc. | Smoke detector/alarm panel interface unit |
US5971095A (en) * | 1996-01-09 | 1999-10-26 | Schlumberger Technology Corporation | Noise filtering method for seismic data |
US5896088A (en) * | 1997-04-16 | 1999-04-20 | Southeastern Univ. Research Assn. | Incipient fire detection system |
US20020010545A1 (en) * | 2000-03-01 | 2002-01-24 | Sitton Gary Arthur | Method for vibrator sweep analysis and synthesis |
US6791453B1 (en) * | 2000-08-11 | 2004-09-14 | Walter Kidde Portable Equipment, Inc. | Communication protocol for interconnected hazardous condition detectors, and system employing same |
US6834235B2 (en) * | 2000-12-07 | 2004-12-21 | Westerngeco, L.L.C. | Method of and apparatus for processing seismic data |
US6882272B2 (en) * | 2001-06-02 | 2005-04-19 | Robert Bosch Gmbh | Danger detecting system |
US6829538B2 (en) * | 2002-04-23 | 2004-12-07 | Westerngeco, L.L.C. | Method for suppressing systematic amplitude variations |
US20060074562A1 (en) * | 2002-05-28 | 2006-04-06 | Lasse Amundsen | Processing seismic data |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070258323A1 (en) * | 2006-05-07 | 2007-11-08 | Ping An | Method for Seismic Trace Decomposition and Reconstruction Using Multiple Wavelets |
US20080285383A1 (en) * | 2006-05-07 | 2008-11-20 | Ping An | System and method for processing seismic data for interpretation |
US8976624B2 (en) | 2006-05-07 | 2015-03-10 | Geocyber Solutions, Inc. | System and method for processing seismic data for interpretation |
US8797828B1 (en) * | 2010-02-15 | 2014-08-05 | Soreq NRC | Remote optical seismic surveying and detection and imaging of underground objects |
Also Published As
Publication number | Publication date |
---|---|
US7889598B2 (en) | 2011-02-15 |
GB2402216A (en) | 2004-12-01 |
GB2402216B (en) | 2007-01-03 |
CA2527215A1 (en) | 2004-12-09 |
GB0411672D0 (en) | 2004-06-30 |
WO2004106975A1 (en) | 2004-12-09 |
CA2831246A1 (en) | 2004-12-09 |
NO20042162L (en) | 2004-11-29 |
FR2855618A1 (en) | 2004-12-03 |
US20080049552A1 (en) | 2008-02-28 |
FR2855618B1 (en) | 2005-08-05 |
NO335231B1 (en) | 2014-10-27 |
MXPA05012828A (en) | 2006-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7940601B2 (en) | Method for computing an exact impulse response of a plane acoustic reflector at zero offset due to a point acoustic source | |
US9625593B2 (en) | Seismic data processing | |
US6668228B1 (en) | Method of attenuating noise in three dimensional seismic data using a projection filter | |
US20120275267A1 (en) | Seismic Data Processing | |
US9158018B2 (en) | Waveform inversion using a response of forward modeling | |
WO2009088602A1 (en) | Spectral shaping inversion and migration of seismic data | |
EP1254383A1 (en) | System and method for seismic wavefield separation | |
CA2357604A1 (en) | Method for identifying and removing multiples from seismic reflection data | |
US20090180351A1 (en) | 3-d tau-p interpolation | |
Soubaras | Angle gathers for shot-record migration by local harmonic decomposition | |
US7889598B2 (en) | Method of seismic processing for the decomposition of a wavefield into harmonic components and applications to the determination of angular gathers of reflectivity | |
Hertweck et al. | Aperture effects in 2.5 D Kirchhoff migration: A geometrical explanation | |
Harding | Slowness—time mapping of near offset seismic reflection data | |
CN110967734B (en) | Virtual source reconstruction method and system based on fast Fourier transform | |
Ventosa et al. | Window length selection for optimum slowness resolution of the local-slant-stack transform | |
Schuster | Seismic interferometry | |
US7079954B2 (en) | Method for selective recording of SH waves using an array of sensors to filter out all non SH waves | |
Kuehl | Least-squares wave-equation migration/inversion | |
Tsai | Use of autoconvolution to suppress first-order, long-period multiples | |
Langston et al. | The validity of ray theory approximations for the computation of teleseismic SV waves | |
Soubaras | The local harmonic decomposition: a tool for extracting angle information from wavefields | |
Stelwagen | Modeling of Elastic Wave Phenomena Near a Plane Interface | |
Guerrero et al. | VTI anisotropy parameter estimation in the tau-p domain: An example from the North Sea | |
Guo | Superresolution Imaging Using Resonant Multiples and Plane-wave Migration Velocity Analysis | |
Fourmann | Minimum-phase decon of band-limited data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: COMPAGNIE GENERALE DE GEOPHYSIQUE, FRANCE Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:SOUBARAS, ROBERT;REEL/FRAME:015831/0386 Effective date: 20040525 |
|
AS | Assignment |
Owner name: COMPAGNIE GENERALE DE GEOPHYSIQUE - VERITAS, FRANC Free format text: CHANGE OF NAME;ASSIGNOR:COMPAGNIE GENERALE DE GEOPHYSIQUE;REEL/FRAME:021339/0757 Effective date: 20070114 Owner name: CGGVERITAS SERVICES SA, FRANCE Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:COMPAGNIE GENERALE DE GEOPHYSIQUE - VERITAS;REEL/FRAME:021339/0764 Effective date: 20080717 Owner name: COMPAGNIE GENERALE DE GEOPHYSIQUE - VERITAS,FRANCE Free format text: CHANGE OF NAME;ASSIGNOR:COMPAGNIE GENERALE DE GEOPHYSIQUE;REEL/FRAME:021339/0757 Effective date: 20070114 Owner name: CGGVERITAS SERVICES SA,FRANCE Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:COMPAGNIE GENERALE DE GEOPHYSIQUE - VERITAS;REEL/FRAME:021339/0764 Effective date: 20080717 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |