Method and apparatus for modeling and separation of primaries and internal multiples using solutions of the two-way wave equation

Patent metadata

Field	Value
Publication number	US-2016334526-A1
Application number	US-201515109689-A
Country	US
Kind code	A1
Filing date	Jan 12, 2015
Priority date	Jan 13, 2014
Publication date	Nov 17, 2016
Grant date	—

Abstract

Method and apparatus for seismic data processing estimate primaries and/or internal multiples by solving a two-way wave equation using first and later arrivals from each layer of an underground formation explored using waves, and summing contributions of the layers to the primaries and/or the internal multiples, without using adaptive subtraction.

First claim

1 . A method, comprising:
receiving data recorded by sensors while an underground formation is explored using waves; obtaining a layer model that specifies one or more impedance and/or velocity changes inside the underground formation, a layer of the model layer being defined between adjacent among the one or more impedance and/or velocity changes; extracting from the data, first and later arrivals of the waves emerging from each of the one or more impedance and/or velocity changes; and estimating at least one of primaries and internal multiples by solving a two-way wave equation for each layer of the model layer, using the first and the later arrivals, and summing portions of resulting solutions to obtain the primaries and/or to the internal multiples.

2 . The method of claim 1 , wherein the layer model is obtained from log data including measurements of wave velocity and/or of density inside the underground formation.

3 . The method of claim 1 , wherein the layer model is obtained by inversion.

4 . The method of claim 1 , wherein the layer model is a two-dimensional model.

5 . The method of claim 1 , wherein, a contribution, ΔP(ξ), to the primaries due to a layer from the depth ξ to a next depth ξ+Δξ is
Δ P (ξ)=Δ R (ξ)(1− H ( t −(τ( xg,yg,zg,xm,ym ,ξ)− T (ξ))))

where H is Heaviside function,
Δ R (ξ)= R ( xs,ys,zs,xg,yg,zg,m (ξ))− R ( xs,ys,zs,xg,yg,zg,m (ξ+Δξ)),

R(xs, ys, zs, xg, yg, zg, m (ξ)) and R(xs, ys, zs, xg, yg, zg, m (ξ+Δξ)) being full solutions of the two-way wave equation for the depth ξ and the next depth ξ+Δξ, respectively, with
m (ξ)= m ( xm,ym,zm )(1− H ( z −ξ))+ H ( z −ξ) m ( xm,ym ,ξ) and

m (ξ+Δξ)= m ( xm,ym,zm )(1− H ( z−H ( z −ξ−Δξ))+ H ( z −ξ) m ( xm,ym ,ξ+Δξ),

for a wave source being located at xs, ys, zs, a sensor being located at xg, yg, zg, the layer model being m(xm, ym, zm), and m(xm, ym, ξ) and m(xm, ym, ξ+Δξ) representing the wave propagating with the constant velocity below depth ξ and ξ+Δξ, respectively, and τ(xg, yg, zg, xm, ym, ξ) being a two-way travel time from the sensor to the layer.

6 . The method of claim 5 , further comprising:
determining the two-way travel time analytically, or determining the two-way travel time using a ray tracing method.

7 . The method of claim 5 , wherein the primaries are calculated by adding contributions of layers above a detection depth ξ max :
P ( xs,ys,zs,xg,yg,zg,m )=Σ ξ=ξmax ξ=0 ΔP (ξ)Δξ.

8 . The method of claim 7 , further comprising:
subtracting the calculated primaries from the data, without using adaptive subtraction; and generating an image of the underground formation using a result of the subtracting.

9 . The method of claim 7 , wherein the internal multiples, M, are calculated by subtracting the primaries, P, from a full solution of the two-way equation, R:
M ( xs,ys,zs,xg,yg,zg,m )= R ( xs,ys,zs,xg,yg,zg,m )− P ( xs,ys,zs,xg,yg,zg,m ).

10 . The method of claim 9 , further comprising:
subtracting the internal multiples from the data without using adaptive subtraction; and generating an image of the underground formation using a result of the subtracting.

11 . The method of claim 5 ,
wherein a contribution, ΔM(ξ), to the internal multiples due to the layer from ξ to ξ+Δξ is ΔM(ξ)=ΔR(ξ)−ΔP(ξ), and the internal multiples, M, are calculated by adding contributions of all layers up to a detection depth ξ max :
M ( xs,ys,zs,xg,yg,zg,m )=Σ ξ=ξmax ξ=0 ΔM (ξ)Δξ.

12 . A data processing apparatus, the apparatus comprising:
an interface configured to receive log data and data recorded by sensors while an underground formation is explored using waves; and a data processing unit configured
to obtain a layer model from the log data, the layer model specifying one or more impedance changes inside the underground formation, a layer of the model layer being defined between adjacent among the one or more impedance and/or velocity changes;
to extract from the data, first and later arrivals of the waves emerging from each of the one or more impedance changes; and
to estimate at least one of primaries and internal multiples by solving a two-way wave equation using the first and later arrivals, for each layer of the model layer, and then summing portions of resulting solutions to obtain the primaries and/or to the internal multiples.

13 . The apparatus of claim 12 , wherein the log data includes measurements of wave velocity and of density inside the underground formation.

14 . The apparatus of claim 12 , wherein the data processing unit solves the two-way wave equation to obtain a contribution, ΔP(ξ), to the primaries due to a layer from the depth ξ to a next depth ξ+Δξ as
Δ P (ξ)=Δ R (ξ)(1− H ( t −(τ( xg,yg,zg,xm,ym ,ξ)− T (ξ))))

where H is Heaviside function,
Δ R (ξ)= R ( xs,ys,zs,xg,yg,zg,m (ξ))− R ( xs,ys,zs,xg,yg,zg,m (ξ+Δξ)),

R(xs, ys, zs, xg, yg, zg, m (ξ)) and R(xs, ys, zs, xg, yg, zg, m (ξ+Δξ)) being full solutions of the two-way wave equation for the depth ξ and the next depth ξ+Δξ, respectively, with
m (ξ) m ( xm,ym,zm )(1− H ( z −ξ))+ H ( z −ξ) m ( xm,ym ,ξ) and

m (ξ+Δξ)= m ( xm,ym,zm )(1− H ( z −ξ−Δξ)+ H ( z −ξ) m ( xm,ym ,ξ+Δξ),

for a wave source being located at xs, ys, zs, a sensor being located at xg, yg, zg, the layer model being m(xm, ym, zm), and m(xm, ym, ξ) and m(xm, ym, ξ+Δξ) representing the wave propagating with the constant velocity below depth and respectively, and τ(xg, yg, zg, xm, ym, ξ) being a two-way travel time from the sensor to the layer.

15 . The apparatus of claim 14 , wherein the data processing unit further configured to determine the two-way travel time analytically, or using a ray tracing method.

16 . The apparatus of claim 14 , wherein the data processing unit is further configured
to calculate the primaries P(xs, ys, zs, xg, yg, zg, m) by adding contributions of layers above a detection depth ξ max , to subtract the calculated primaries from the data, without using adaptive subtraction, and to generate an image of the underground formation using a result of the subtracting.

17 . The apparatus of claim 14 , wherein the data processing unit is further configured
to calculate the primaries P(xs, ys, zs, xg, yg, zg, m) by adding contributions of layers above a detection depth ξ max , and to calculate the internal multiples, M, by subtracting the primaries, P, from a full solution of the two-way equation, R:
M ( xs,ys,zs,xg,yg,zg,m )= R ( xs,ys,zs,xg,yg,zg,m )− P ( xs,ys,zs,xg,yg,zg,m ).

18 . The apparatus of claim 14 , wherein the data processing unit is further configured
to calculate a contribution, ΔM(ξ) to the internal multiples due to the layer from ξ to ξ+Δξ is ΔM(ξ)=ΔR(ξ)−ΔP(ξ) and the internal multiples, M, by adding contributions of all layers up to a detection depth ξ max :
M ( xs,ys,zs,xg,yg,zg,m )=Σ ξ=ξmax ξ=0 ΔM (ξ)Δξ.

19 . The apparatus of claim 18 , wherein the data processing unit is further configured
to subtract the internal multiples from the data without using adaptive subtraction; and to generate an image of the underground formation using a result of the subtracting.

20 . A non-transitory computer readable medium storing executable codes, which, when executed by a computer make the computer perform a method for processing data recorded by sensors while an underground formation is explored using waves, the method comprising:
obtaining a layer model from log data, the layer model specifying one or mor

Assignees

• Cgg Services Sa

Inventors

• Pica Antonio

Classifications

• G01V2210/671

  Raytracing · CPC title

• G01V1/362

  Effecting static or dynamic corrections; Stacking · CPC title

• G01V1/282 Primary

  Application of seismic models, synthetic seismograms · CPC title

• G01V2210/56

  De-ghosting; Reverberation compensation · CPC title

• G01V99/00

  Subject matter not provided for in other groups of this subclass · CPC title