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
  • G01V1/366—Seismic filtering by correlation of seismic signals

Definitions

• the invention relates to a method for obtaining a representation particularly suitable for correlating signals in identical frequency bands and the application of this method in geophysics to stratigraphic calibrations and to calibrations between well data and seismic data.
• the curves to be compared can be of the same nature, that is to say represent the variations of the same first quantity as a function of the same second quantity, or of different nature. It may for example be records of the same physical phenomenon but shifted in time or space, or records relating to different physical phenomena or even records relating to the same physical phenomenon recorded for example by methods different so that their frequency content is different.
• Correlations can be made digitally.
• the result obtained is generally global and unreliable if one does not make constraining assumptions on the signals, the method then consisting in choosing between several autocorrelation peaks.
• the correlation can be carried out visually, by manual shifting of one of the curves with respect to the other along the axis of the second quantity. We thus seek by successive shifts the optimal similarity on one or more portions of the curve.
• This method makes it possible to take knowledge into account a priori. It is that which is commonly used in geophysics for the depth or time calibration of seismic horizons or for the correlation of recordings made in a well and of seismic recordings.
• the main drawback of such a method resides in the difficulty that there is in comparing signals of forms which can be very different, for example if their frequency content is different.
• the method according to the invention allows the improvement of the reliability of the correlations carried out visually between at least two recordings thanks to 1 • introduction of time-frequency analyzes to compare the signals in identical frequency bands.
• There is a known method for analyzing a signal making it possible to decompose said signal as a sum of elementary functions a fa or wavelets which each vibrate like sinusoids over a range whose position on the x axis is linked to the parameter b and whose the width is linked to the parameter a (center frequency), and which are very strongly damped outside this range.
• the decomposition of a signal using a family of these wavelets constitutes what is called a "time - frequency" analysis because the first and most common decompositions were carried out on records of variations of a first quantity as a function of time (the second quantity).
• the dimension of the parameter b is that of a time
• the dimension of the parameter a is the dimension of the inverse of a time therefore of a temporal frequency.
• o defines the "wavelet transform" in two dimensions and z associated with a recording along the x axis as l following the coefficients c a / b each corresponding to the integral of the product of the record to be analyzed by the elementary wavelet a ⁇ according to the values of b along the x axis and the values of a along a z axis.
• complex wavelets have been chosen to perform the time-frequency analysis of a recording or a signal, it becomes possible to define the real part, the imaginary part, the module or even the phase of transformed it into wavelets.
• the method for obtaining a representation which is particularly suitable for correlating, in the same frequency band [fl, f2] along an axis z, of at least two one-dimensional records si and s2, variations of a first quantity ( yl and y2) function of a second quantity (xl and x2) varying along an axis x, is characterized in that:
• a family of analysis wavelets is defined by choosing a wavelet function y a £, (*) depending on a parameter a representing the frequency of the wavelet and a parameter b defining the position of the wavelet on The x axis and by defining for the frequency has a sequence of m values going from i ⁇ to f 2 , m being at least equal to the number of octaves included in the frequency range [f- ⁇ , î_ ⁇ ' e * " for parameter b a sequence of n values,
• At least one matrix is constructed whose coefficients A ab are a quantity, namely either the module, the phase, the real part or even the imaginary part of each of the coefficients C a £ as a function of the values increasing from a and b,
• - one transposes according to the frequency axis the matrix or matrices associated with one of the two recordings, and - one represents by pairs, for the purpose of optical correlation, the transposed matrix associated with one of said quantities and corresponding to one of the recordings and the non-transposed matrix associated with the same quantity and corresponding to the other record, by placing said transposed and non-transposed matrices side by side so that all the coefficients A ab corresponding to the same value of b are in the same alignment and that when we describe each alignment we meet the coefficients of one of the matrices according to the decreasing values of a then the coefficients of the second matrix according to the increasing values of a or, preferably, we meet the coefficients of one of the matrices according to the increasing values of a then the coefficients of the second matrix according to the decreasing values of a.
• a palette of colors or gray adapted to the values taken by the various A a ⁇ .
• the wavelet function ⁇ a ⁇ (x) chosen to build the analysis wavelet family can be, for example, a Gaussian, a square or triangle function, this wavelet function can be real or complex and give rise to a family d 'wavelets orthogonal or not between them.
• the functions of the family are elementary wavelets deduced from each other by dilation.
• m is at least equal to the number of octaves included in the fixed analysis range [fl, f2].
• the n values forming the sequence of values assigned to parameter b for the analysis correspond preferably at the coordinates of the center of the elementary ranges of analysis on the x axis.
• complex elementary wavelets will be chosen so as to represent the phase of transforms, for example so-called "MORLET” wavelets.
• the method of the invention finds a particularly interesting application in the field of geophysics, in particular for so-called stratigraphic setting operations.
• the operations consist in correlating at least one seismic record (or seismic trace), located at the right or in the vicinity of a wellbore with a synthetic trace constructed most often from measurements carried out in the well, this in order to associate with precision to each of the main seismic markers a recognized geological level in the well.
• the representation method according to the invention makes it possible to further improve the optical perception and complete it by providing additional information due to the transform.
• - Figure 1 represents a seismic trace T 1 (t) of a seismic profile representing the ground accelerations at a location M as a function of time t for times varying between 200 and 3800 milliseconds as well as the phase and the modulus of its time - frequency transform, in 4 ms steps, with as many pixels as samples,
• FIG. 2 represents another trace T2 (t) representing the accelerations of the subsoil at location M and belonging to a second seismic profile crossing the first profile at location M, as well as the phase and the module of its time-frequency transform,
• FIG. 3 shows a representation according to the invention allowing easy optical correlation of the two traces T ⁇ (t) and T 2 (t), from the use of the phase d their time - frequency transforms.
• the representation is limited to 2400 ms and the number of frequencies reduced.
• i is the imaginary number whose square is equal to -1
• e is the base of the natural logarithms
• Ln (2) is the natural logarithm of 2.
• the coefficients c ab ⁇ ) of the wavelet transform are the result of the continuous sum of the product of the record T (t) representing the continuation of the elementary functions
• n analysis points b has been defined as identical following the sampling times of the traces ⁇ and T 2 (every 4 ms).
• the result of the time-frequency analysis of a recording T (l) is therefore here made up of the series of 65 time variation curves C ab , each of the 65 curves being associated with one of the analysis frequencies and consisting in a sequence of n coefficients C ab .
• the elementary functions used here being complex, it is advisable to represent either the module, or the phase, or the real part, or even the imaginary part of the transform. We have chosen to represent the phase and the module.
• Each of the monodi ensional traces T (t) of axis t thus gives after analysis two two-dimensional representations, image of two matrices A (TI), P (TI) one of whose dimensions is the same as that of the trace analyzed (axis t ) and whose other dimension is a frequency axis f.
• the matrix A is representative of an amplitude and the matrix P of a phase, for each point (t, f) considered.
• the most efficient graphical representation to perform this optical correlation consists in transposing one of the matrices associated with - ⁇ along the axis f and in juxtaposing it along the axis t with the corresponding matrix associated with T 2 so that the coefficients of said matrices corresponding to the same t are on the same alignment and when we describe each alignment we meet the coefficients of one of the matrices according to the increasing values of the frequency then the coefficients of the second matrix according to the decreasing values of the frequency.
• This takes advantage of the symmetry detection power of the human eye.
• the axis t is vertical and the various axes f are horizontal.
• FIG. 3 there are two representations according to the invention.
• the two matrices V_ ⁇ ⁇ T_) and P 2 (T) have been juxtaposed, P ⁇ ⁇ T- j ⁇ appearing on the left of the image (low frequencies on the left and high frequencies on the right), and P 2 (T) appearing on the right
• the image of the energies makes it possible to identify at a glance the frequency zones where there is energy in the two signals.
• the image of the phases allows, by symmetry, to distinguish common configurations, first in the low frequencies by comparing the extreme left and right parts of the representation, then to refine in the high frequencies by comparing the central parts of the representation.
• the representation method according to the invention can be automatically executed using a microcomputer or station and software allowing the calculation of the time-frequency transform of a signal, and equipped with a high definition graphic screen allowing the visualization of matrices and curves.
• the calculation was carried out by software written in F77 and C language on a computer of the HP 9000 s 300 type and the data were viewed on the same computer equipped with a CH screen (HP standard) on which was running an X Windows server (XII, R 3 ).
• the invention is not limited to the example described and shown, various modifications that can be made without departing from the scope of the invention.
• the method applies in particular to the correlation of all types of record, real or complex, whatever the variable.

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)
• Image Analysis (AREA)
• Geophysics And Detection Of Objects (AREA)
• Complex Calculations (AREA)

Applications Claiming Priority (2)

Publications (1)

Family

ID=9412041

Family Applications (1)

Country Status (4)

Families Citing this family (1)

Family Cites Families (2)

• 1991
  • 1991-04-19 FR FR9104840A patent/FR2675589B1/fr not_active Expired - Fee Related
• 1992
  • 1992-04-16 EP EP19920910310 patent/EP0535208A1/de not_active Withdrawn
  • 1992-04-16 CA CA002085616A patent/CA2085616A1/fr not_active Abandoned
  • 1992-04-16 WO PCT/FR1992/000344 patent/WO1992018941A1/fr not_active Application Discontinuation

Non-Patent Citations (1)

Also Published As

Similar Documents

Legal Events