Method for seismic processing, in particular for compensating birefringence on seismic traces

• US 7,436,735 B2
• Filed: 11/07/2002
• Issued: 10/14/2008
• Est. Priority Date: 11/08/2001
• Status: Expired due to Fees

• Current Assignee: Cggveritas Services SA (CGG)
• Original Assignee: Compagnie Generale De Geophysique
• Inventors: Garotta, Robert
• Primary Examiner(s): Keith; Jack W
• Assistant Examiner(s): Hughes; Scott A
• Application Number:
  US10/495,241
• Publication Number:
  US 20050013195A1
• Time in Patent Office: 2,168 Days
• Field of Search: 367/38, 367/50, 367/75, 367/31, 702/14
• US Class Current: 367/75
• CPC Class Codes: G01V 1/28 Processing seismic data, e.g. for interpretation or for event detection G01V1/48 takes precedence
• CPC Class Codes: G01V 2210/58 Media-related

Abstract

A method of seismic processing for obtaining information on the geophysics of the subsoil includes acquiring seismic traces at at least one point on the surface of the subsoil or in the subsoil, the seismic traces corresponding on each occasion to two perpendicular components of a shear wave emitted into the subsoil and reflected by different interfaces therein; applying a succession of transformations (⊥) at least to a temporal portion of the traces for each of these assumptions, determining the value of a parameter representative of the coherence/similarity between the result traces obtained in this way; and election as a function of the values obtained in this way that one of the hypotheses which is considered as being most representative of the subsoil, the two result traces obtained for said hypothesis being compensated traces of the subsoil birefringency.

Claims

1. A method of seismic processing for obtaining information on the geophysics of subsoil, the method comprising:

acquiring seismic traces at at least one point on the surface of the subsoil, the seismic traces corresponding on each occasion to two perpendicular components of a shear wave emitted into the subsoil and reflected by different interfaces therein;

applying a succession of inverse transformations (⊥

) at least to a temporal portion of the seismic traces to obtain result traces, the inverse transformations to restore transformations that it is assumed that the shear wave has encountered on passing through a succession of layers in the subsoil, the layers being analyzed as a whole, and an operation of applying the succession of inverse transformations being repeated for different assumptions concerning birefringence in the layers;

for each of these assumptions, determining a value of a parameter representative of the coherence/similarity between the result traces; and

selecting as a function of the values obtained in this way that one of the assumptions which is considered as being most representative of the subsoil, two final result traces obtained for said assumption being compensated for the subsoil birefringency, wherein the two final result traces are obtained by analyzing the layers as a whole.

2. The method according to claim 1, wherein an assumption concerning the birefringency of a layer is characterized by at least one parameter relating to the orientation of the birefringency axes and by at least one parameter relating to the delay imparted by said layer between propagation oriented on a fast axis and propagation oriented on a slow axis, and in that the results are represented or displayed in two dimensions or in three dimensions with, for different portions of the traces, successions of values corresponding to the values of said parameters for a succession of layers through which the wave has passed and to the selected assumption concerning birefringency.

3. The method according to claim 2, wherein comparison attributes are determined on the result traces and represented or displayed together with the values of the orientation and delay parameters.

4. The method according to claim 1, wherein in order to select a birefringency assumption that is considered as being the most representative, the value of the parameter representing coherence/similarity of the result traces obtained for each assumption is compared with a threshold value.

5. The method according to claim 1, wherein successive iterations are performed during which a new birefringency assumption is selected for each operation as a function of the value calculated at the preceding iteration for the parameter representative of coherence/similarity between the result traces obtained, and in that in order to select the birefringency assumption that is considered as being the most representative, it is determined whether the value obtained for the parameter representative of coherence/similarity between the result traces at the iteration in question lies within a given proximity range relative to the values obtained for the same parameter during preceding iterations.

6. The method according to claim 1, wherein it is implemented on trace portions having a duration of the order of 100 milliseconds or longer.

7. The method according to claim 6, wherein trace portions are sampled with a sampling time of the order of 2 or 4 milliseconds or shorter, each of sampling intervals or groups of successive intervals corresponding to a layer under study.

8. A seismic processing method for revealing fractures in subsoil, comprising:

acquiring seismic traces at at least one point on the surface of the subsoil, the seismic traces corresponding on each occasion to two perpendicular components of a shear wave emitted into the subsoil and reflected by different interfaces therein;

applying a succession of inverse transformations (⊥

) at least to a temporal portion of the traces to obtain result traces, the inverse transformations to restore transformations that it is assumed that the shear wave has encountered on passing through a succession of layers in the subsoil, the layers being analyzed as a whole, and an operation of applying the succession of inverse transformations being repeated for different assumptions concerning birefringence in the layers;

for each of these assumptions, determining a value of a parameter representative of the coherence/similarity between the result traces obtained; and

selecting as a function of the values obtained in this way that one of the assumptions which is considered as being most representative of the subsoil, two final result traces obtained for said assumption being compensated for the subsoil birefringency, wherein the two final result traces are obtained by analyzing the layers as a whole.

9. A method for obtaining information on the geophysics of subsoil, the method comprising:

acquiring seismic traces at least one point on the surface of the subsoil, the seismic traces corresponding on each occasion to two perpendicular components of a shear wave emitted into the subsoil and reflected by different interfaces therein;

applying a succession of inverse transformations (⊥

) at least to a temporal portion of the traces to obtain result traces, the inverse transformations to restore the transformations that it is assumed that the shear wave has encountered on passing through a succession of layers in the subsoil, the layers being analyzed as a whole, and an operation of applying the succession of inverse transformations being repeated for different assumptions concerning birefringence in the layers;

for each of these assumptions, determining a value of a parameter representative of the coherence/similarity between the result traces.

10. A method of seismic processing for obtaining information on the geophysics of subsoil, the method comprising

a) acquiring a pair of seismic traces at at least one point of the surface of the subsoil, said seismic traces corresponding to two perpendicular components of a shear wave emitted in the subsoil and reflected by different interfaces corresponding to a plurality of successive layers therein,

b) choosing a first succession of sets of birefringence parameters which characterize a succession of birefringence transformations that it is assumed that the shear wave has encountered on passing through said plurality of successive layers,

c) applying at least to a temporal portion of said two traces corresponding to two perpendicular components of the shear wave, a succession of transformations which are the inverses of said succession of birefringence transformations that it is assumed that the shear wave has encountered on passing through said plurality of successive layers,

d) calculating on the pair of results traces thus obtained a cost function value which is based on the coherency/similarity of said results traces,

e) repeating steps c) and d) for other succession of sets of birefringence parameters corresponding to other assumptions concerning birefringence of the plurality of successive layers,

f) selecting, as a function of the cost function value calculated, the succession of sets of birefringence parameters which optimizes the coherency/similarity between the result traces obtained.

11. A method according to claim 10, wherein a cost function from which the cost function value is calculated is a cross-correlation function.

12. A method according to claim 10, wherein a set of birefringence parameters comprises:

at least one parameter relating to the orientation of the birefringence axes, and

at least one parameter relating to the delay imparted by said layer between propagation orientated on a fast axis and propagation orientated on a slow axis.