WO1994010584A1 - Procede d'estimation de la position de reflecteurs acoustiques presents dans un milieu - Google Patents

Procede d'estimation de la position de reflecteurs acoustiques presents dans un milieu Download PDF

Info

Publication number
WO1994010584A1
WO1994010584A1 PCT/FR1993/001040 FR9301040W WO9410584A1 WO 1994010584 A1 WO1994010584 A1 WO 1994010584A1 FR 9301040 W FR9301040 W FR 9301040W WO 9410584 A1 WO9410584 A1 WO 9410584A1
Authority
WO
WIPO (PCT)
Prior art keywords
crp
macromodel
medium
distance
reflectors
Prior art date
Application number
PCT/FR1993/001040
Other languages
English (en)
Inventor
Jean-Paul Jeannot
Original Assignee
Total
Elf Aquitaine Production
Dataid Eurosoft
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Total, Elf Aquitaine Production, Dataid Eurosoft filed Critical Total
Publication of WO1994010584A1 publication Critical patent/WO1994010584A1/fr

Links

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/362—Effecting static or dynamic corrections; Stacking
    • 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/282—Application of seismic models, synthetic seismograms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/52—Move-out correction

Definitions

  • a source emitting seismic waves and a plurality of equidistant receivers are used, the whole moving on the surface of the ground or the sea.
  • FIG. 1 attached illustrates this operation in 0 in the particular case of a horizontal reflector 10.
  • the source moves to the left and successively occupies the positions S ⁇ , S. ⁇ , s i + 2 * 0n n ' has shown in the figure that the rays emitted by the source in its different positions and which are reflected on the same point of the reflector or common point of reflection (CRP).
  • the reflected rays are picked up by the receivers R ⁇ , R3, R5.
  • the middle of the segment joining a source and a receptor is denoted CMP (midpoint).
  • This operation makes it possible to collect traces 0 resulting from the waves reflected by the various reflectors that constitute the interfaces between the geological layers. From these traces, we can deduce the observed travel times, which we will denote by t Q ⁇ g , relating to each source / receiver pair. 5
  • the "stage" technique makes it possible to obtain, in the case of a horizontal stratified medium, like that of FIG. 1, a relatively faithful representation of the basement. It consists in carrying out speed analyzes in CMP, dynamic hyperbolic corrections of the traces and 0 an addition of the corrected traces.
  • This velocity field which will be called a macromodel in the continuation of the text, is defined, starting from the knowledge which one has of the geology of the subsoil, by the following parameters:
  • V ⁇ (x, z) V i (x) + K i (z), x being the horizontal coordinate and K ⁇ the factor compaction in layer i.
  • the propagation of acoustic waves in this macromodel is simulated by carrying out regular sampling of surface sources and receivers. We then launch rays from each of the surface samples, and we determine the travel times t mQ ( - corresponding to the different source / receiver pairs. These times are used either to determine dynamic corrections to refine the stack technique, or to compare them with the propagation times tobs measured on the seismic data. The difference 0 ⁇ g - mQ - will be noted ⁇ t.
  • tomographic inversion has a limitation linked to the fact that is both the defined depth of the interface under the layer i and the indefinite depth of the reflector model for which the best representation in depth is sought.
  • the layers and in particular the function z ⁇ (x) have a certain regularity.
  • the raacro odele found by the inversion technique therefore gives an approximate outline of the shape of the reflectors, which does not show the detail real irregularities, such as faults, etc.
  • the present invention aims to obtain a better estimate of the position of the acoustic reflectors of a medium from an interpretation of the differences between the actual journey times of acoustic waves propagating in said medium and the corresponding modeled journey times obtained. from a parametric macromodel of this medium giving the depth jz of the reflectors and the speed V of the layers of the medium as a function of at least the horizontal coordinate x_.
  • the method according to the invention is of the type consisting in: collecting seismic traces by means of a recording system comprising at least one source emitting acoustic waves and several receivers which are moved on the surface of the medium, in particular land or sea surface,
  • the tracing of rays is done not from the surface of the medium, as in the prior art, but from each CRP.
  • we carry out a regular sampling of a reflector of the macroraodel we set a first distance source / receiver or offset, we trace the rays corresponding to all the CRP, then we give our a second offset and we trace the departments for the same CRP, and so on with multiple offsets.
  • This method makes it possible to determine a set of radii corresponding to variable offsets at the surface, but which are reflected at the same CRP, and therefore travel times c mod corresponding P ° ndant to said CRP.
  • the deviations ⁇ obtained with all the rays are not due to background noise, but are carriers of information.
  • these deviations are transformed into distance corrections ⁇ H along a predetermined direction to provide a corrected position of the CRP.
  • the transformation of all of the ⁇ t relating to the same CRP into the corresponding correction ⁇ H consists in transforming the ⁇ t into distance corrections ⁇ h in the predetermined direction, and from ⁇ h, to calculate an average distance correction taken as the correction ⁇ H.
  • said transformation consists in:
  • the components ⁇ z, ⁇ x and possibly ⁇ y of each ⁇ h can be sought on the axes of coordinates z, x and possibly y of the macromodel, said components representing corrections on the coordinates of the modeled CRP considered.
  • these corrections are displayed on the macromodel in order to obtain the corrected position of the reflectors.
  • Visualization can be carried out by any known method.
  • the corrections ⁇ h corresponding to all the paths relating to each CRP can be transferred to the normal to the reflector at said CRP.
  • Each corrected position is marked by a line which can be assimilated to an impulse as indicated above.
  • Two cases can then arise: - the corrected positions of the CRP corresponding to ⁇ h are dispersed in a large area located in the vicinity of the position found by the inversion.
  • the effective position of the CRP is better visualized by convolving all the pulses obtained by a given signal, for example a signal of the wavelet type of Ricker or the like.
  • the peak of the sum signal obtained thus provides an average position of the corrected CRP distant from ⁇ H of the modeled CRP taken as origin.
  • the amplitude of said sum signal is an indicator of the dispersion of the result;
  • the corrected positions of the CRP corresponding to the ⁇ h coincide or are concentrated around the same average position.
  • the amplitude of the sum signal resulting from a convolution as mentioned above is the sum of the amplitudes of all the pulses.
  • the treatment according to the invention is the equivalent of migration before stack.
  • FIG. 2 represents, in a macromodel, a path of seismic radius relating to a modeled CRP and a path of real radius associated with the same source / receiver couple;
  • FIG. 3 schematically represents the case where all the depth corrections coincide;
  • FIG. 4 represents the case where the depth corrections are dispersed.
  • the macromodel obtained by tomographic inversion comprises several layers bounded together by modeled interfaces such as i.
  • the macromodel is parameterized by the depth z ⁇ (x) of each point of an interface i and by the propagation speed V ⁇ (x) of the seismic wave in the geological layer surmounting the interface i.
  • the quantities z ⁇ (x) and V j _ (x) are, for example, functions of type spline function of the horizontal coordinate x.
  • the set of deviations ⁇ t corresponding to all the CRPs can be minimized by using an inversion algorithm. At convergence, it may happen that for each seismic radius, the residual difference is not zero, which can be explained in particular by the simplification of the real structure of the subsoil resulting from its parameterization.
  • the invention consists in using said residual difference over the travel time of each ray to approach the exact location of the reflector.
  • the residual differences in travel time ⁇ t are converted into distance corrections ⁇ h or into their components depth corrections ⁇ z and horizontal position corrections ⁇ x and possibly ⁇ y.
  • this depth correction is shown for the particular path 12 which starts from the source S, is reflected on the CRP considered 13 and ends at the receiver R located at a distance -d from the source. This correction is brought to the normal 24 to the reflector modeled jL at the CRP 13 considered.
  • the reference 14 designates the path corresponding to a radius which starts from S, is reflected on the corrected CRP 15 and ends in R.
  • FIGS. 3 and 4 The method according to the invention is illustrated by FIGS. 3 and 4.
  • the distance corrections ⁇ h corresponding to all the rays which illuminate said CRP when ⁇ 1 varies, each correction ⁇ h being materialized by a line 16 perpen ⁇ dicular to normal 24 and representative of a pulse.
  • Figure 3 shows the case where all the corrections of depth 16 coincide.
  • the pulse corresponding to each correction is convoluted by a signal, for example of the Ricker wavelet type. All the signals being in coincidence, the result of the convolution is is a signal 18 whose maximum amplitude is equal to the sum of the amplitudes of the pulses.
  • FIG. 4 corresponds to the case where the distance corrections ⁇ h are spread out.
  • the corrections are convolved by a given signal, for example of the Ricker wavelet type, the signals 18, 20 resulting from the convolution then being added.
  • the dispersion of the corrections results in an amplitude of the sum smaller than in the previous case.
  • the advantage of the method according to the invention is that it makes it possible to obtain a synthetic image of the reflector closer to reality.
  • high frequencies which have been lost in the representation of the macromodel by spline functions are recovered.

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)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

Le procédé selon l'invention consiste à: collecter des traces sismiques au moyen d'une source (S) émettrice d'ondes acoustiques et de plusieurs récepteurs (R) que l'on déplace à la surface du milieu; déterminer les temps de trajets réels observés (tobs) correspondant à des réflexions sur les différents réflecteurs du milieu; définir un macromodèle paramétré donnant la profondeur z^_ des réflecteurs (i) et la vitesse V des couches du milieu en fonction de la coordonnée horizontale x^_; détermine les temps de trajet source/récepteur modélisés (tmod) par la technique du tracé de rayons, à partir de CRP, choisis sur les réflecteurs du macromodèle; transformer l'ensemble des écarts Δt = tobs - tmod obtenus pour tous les rayons du macromodèle issus du CRP concerné en une correction de distance ΔH le long d'une direction (24) passant par ledit CRP; et déplacer le CRP du macromodèle de la distance ΔH le long de ladite direction prédéterminée pour fournir une position corrigée dudit CRP.

Description

PROCEDE D'ESTIMATION DE LA POSITION DE REFLECTEURS ACOUSTIQUES PRESENTS DANS UN MILIEU.
Dans une opération de prospection sismique classique par réflexion, on utilise une source émettrice d'ondes sismiques et une pluralité de récepteurs équidistants, l'ensemble se déplaçant à la surface du sol ou de la mer.
La figure 1 ci-jointe illustre cette opération dans 0 le cas particulier d'un réflecteur horizontal 10. La source se déplace vers la gauche et occupe successivement les positions S^, S.^, si+2* 0n n'a représenté sur la figure que les rayons émis par la source dans ses différentes positions et qui se réfléchissent sur un même 5 point du réflecteur ou point commun de réflexion (CRP). Les rayons réfléchis sont captés par les récepteurs R^, R3, R5. Le milieu du segment joignant une source et un récepteur est noté CMP (point milieu).
Cette opération permet de recueillir des traces 0 résultant des ondes réfléchies par les différents réflecteurs que constituent les interfaces entre les couches géologiques. De ces traces, on peut déduire les temps de trajet observés, que l'on notera tQ^g, relatifs à chaque couple source/récepteur. 5 La technique du "stade" permet d'obtenir, dans le cas d'un milieu stratifié horizontal, comme celui de la figure 1, une représentation relativement fidèle du sous-sol. Elle consiste à effectuer des analyses de vitesse en CMP, des corrections dynamiques hyperboliques des traces et 0 une addition des traces corrigées.
Une des méthodes les plus couramment utilisées pour valider cette représentation technique consiste à lancer des rayons à partir de la source jusqu'au récepteur et à calculer les temps odélisés de propagation, notés mo(j, 5 le long de ces trajets.
Pour lancer les rayons, il faut connaître le champ de vitesse du milieu dans lequel s'effectue la propagation. Ce champ de vitesse, qui sera appelé macromodèle dans la suite du texte, est défini, à partir de la connaissance que l'on a de la géologie du sous-sol, par les paramètres suivants :
- profondeur z^(x) de l'interface (réflecteur) se trouvant sous la couche i,
- vitesse de propagation V^(x,z) dans la couche i, par exemple du type Vi(x,z) = Vi(x)+Ki(z), x étant la coordonnée horizontale et K^ le facteur de compaction dans la couche i. Dans la technique habituelle, on simule la propagation d'ondes acoustiques dans ce macromodèle en effectuant un échantillonnage régulier de sources et de récepteurs en surface. On lance ensuite des rayons depuis chacun des échantillons de surface, et on détermine les temps de trajet tmQ(- correspondant aux différents couples source/récepteur. Ces temps sont utilisés soit pour déterminer des corrections dynamiques permettant d'affiner la technique du stack, soit pour les comparer aux temps de propagation tobs mesurés sur les données sismiques. L'écart 0^g - mQ - sera noté Δt.
Si des écarts significatifs sont observés pour l'ensemble des couples (tobs»tmod^» ^es algorithmes d'inversion sont mis en oeuvre afin de modifier! le macromodèle, de manière à réduire lesdits écarts. Lorsque l'inversion aura convergé vers une solution stable, le macromodèle obtenu fournira les temps modélisés m0(j les plus proches des temps observés.
Toutefois, cette technique dite d'inversion tomographique, présente une limitation liée au fait que est à la fois la profondeur définie de l'interface sous la couche i et la profondeur indéfinie du modèle de réflecteur dont on cherche la meilleure représentation en profondeur. Pour que la propagation dans le macromodèle se déroule de manière réaliste, il est indispensable que les couches et en particulier la fonction z^(x) présentent une certaine régularité. Le raacro odèle trouvé par la technique de l'inversion donne donc un tracé approximatif de la forme des réflecteurs, qui ne montre pas le détail des irrégularités réelles, telles que les failles etc..
La présente invention vise à obtenir une meilleure estimation de la position des réflecteurs acoustiques d'un milieu à partir d'une interprétation des écarts entre les temps de trajet réels d'ondes acoustiques se propageant dans ledit milieu et les temps de trajet modélisés correspondants obtenus à partir d'un macromodèle paramétré de ce milieu donnant la profondeur jz des réflecteurs et la vitesse V des couches du milieu en fonction d'au moins la coordonnée horizontale x_.
Le procédé selon l'invention est du type consistant à: collecter des traces sismiques au moyen d'un système d'enregistrement comprenant au moins une source émettrice d'ondes acoustiques et plusieurs récepteurs que l'on déplace à la surface du milieu, notamment surface du sol ou de la mer,
- déterminer à partir de ces traces les temps de trajet réels observés ^obs^ pour chaque couple source/récepteur et correspondant à des réflexions sur les différents réflecteurs du milieu ; définir un macromodèle paramétré donnant la profondeur z des réflecteurs et la vitesse V des couches du milieu en fonction de la coordonnée horizontale x et éventuellement de la coordonnée horizontale^,
- déterminer, à l'aide du macromodèle, des temps de trajet source/récepteur modélisés (tmo(j) par la technique du tracer de rayons,
- calculer, pour chaque réflecteur et chaque rayon, l'écart de temps de trajet Δt = to s-tmod» et, si- besoin est,
- mettre en oeuvre un algorithme d'inversion qui modifie le macromodèle jusqu'à ce que les écarts Δt aient convergé vers une valeur stable minimale, et il se caractérise en ce qu'il consiste à mettre en oeuvre les étapes suivantes :
- effectuer le tracer de rayons à partir de points, dits CRP, choisis sur les réflecteurs du macromodèle, A
- transformer, pour chacun desdits CRP, l'ensemble des écarts Δ t obtenus pour tous les rayons du macromodèle issus du CRP concerné en une correction de distance ΔH le long d'une direction prédéterminée passant par ledit CRP, par exemple la normale au réflecteur contenant le CRP, et
- déplacer le CRP concerné du macromodèle de la distance ΔH le long de ladite direction prédéterminée pour fournir une position corrigée dudit CRP, l'ensemble des positions corrigées des CRP d'un réflecteur du macromodèle représentant une position corrigée améliorée de ce réflecteur.
Ainsi, selon l'invention, le tracer de rayons se fait non pas à partir de la surface du milieu, comme dans la technique antérieure, mais à partir de chaque CRP. Pour cela, on effectue un échantillonnage régulier d'un réflecteur du macroraodèle, on se fixe une première distance source/récepteur ou offset, on trace les rayons correspondant à tous les CRP, puis on se donne un second offset et l'on trace les rayons pour les mêmes CRP, et ainsi de suite avec plusieurs offsets.
Ce procédé permet de déterminer un ensemble de rayons correspondant à des offsets variables en surface, mais qui se réfléchissent au même CRP, et donc des temps de trajet cmod corresP°ndant audit CRP. De plus, contrairement à ce qui était supposé dans le procédé classique, on considère selon l'invention, que les écarts Δ obtenus avec tous les rayons ne sont pas dus au bruit de fond, mais sont porteurs d'information. Selon l'invention, ces écarts sont transformés en corrections de distance ΔH le long d'une direction prédéterminée pour fournir une position corrigée du CRP.
La transformation de l'ensemble des Δt relatifs à un même CRP en la correction correspondante Δ H consiste à transformer les Δt en corrections de distance Δh sur la direction prédéterminée, et à partir des Δh, à calculer une correction moyenne de distance prise comme la correction Δ H. Selon un mode de réalisation, ladite transformation consiste à :
- représenter une même impulsion perpendiculairement à la direction prédéterminée choisie, passant par le CRP, en chacun des points de ladite direction dont les distances à une origine, représentée par la position du
CRP concerné déduite du macromodèle, ont des valeurs Δh représentant des corrections de distance correspondant aux
Δt associés aux différents rayons illuminant le CRP lorsque la distance source/récepteur ou offset varie,
- convoluer chacune desdites impulsions par un même signal donné,
- additionner les signaux résultant de la convolution pour obtenir un signal somme, et - prendre comme valeur ΔH la distance séparant la position du maximum dudit signal somme de l'origine choisie.
Pour obtenir les corrections Δh correspondant aux Δt associés à un même CRP, on peut rechercher les composantes Δz, Δx et éventuellement Δ y de chaque Δh sur les axes de coordonnées z, x et éventuellement y du macromodèle, lesdites composantes représentant des corrections sur les coordonnées du CRP modélisé considéré.
On peut par exemple calculer Δ z et Δx en utilisant les dérivées partielles "h t / δ z et àt/ ^x fournies par certaines méthodes de tracer de rayons pour écrire b / ) z 2+ ^t/ àx Δx =Δt . En introduisant la contrain¬ te supplémentaire que le CRP se déplace le long d'une direction déterminée passant par ledit CRP et, par exemple, le long de la normale au réflecteur du macroraodèle passant par le CRP, on ajoute une équation qui permet de calculer Δ x et Δz.
Selon l'invention, on visualise ces corrections sur le macromodèle afin d'obtenir la position corrigée des réflecteurs. La visualisation peut être effectuée par tout procédé connu. On peut par exemple reporter les corrections Δ h correspondant à tous les trajets relatifs à chaque CRP sur la normale au réflecteur audit CRP. Chaque position corrigée est marquée par un trait qui peut être assimilé à une impulsion comme indiqué plus haut. Deux cas peuvent alors se présenter : - les positions corrigées du CRP correspondant auxΔh sont dispersées dans une large zone située au voisinage de la position trouvée par l'inversion. La position effective du CRP est mieux visualisée en convoluant toutes les impulsions obtenues par un signal donné, par exemple un signal du type ondelette de Ricker ou autre. La crête du signal somme obtenu fournit ainsi une position moyenne du CRP corrigé distante de ΔH du CRP modélisé pris comme origine. L'amplitude dudit signal somme est un indicateur de la dispersion du résultat;
. - les positions corrigées du CRP correspondant aux Δh coïncident ou se concentrent autour d'une même position moyenne. L'amplitude du signal somme résultant d'une convolution telle que précitée est la somme des amplitudes de toutes les impulsions. Ladite position moyenne, dont la distance au CRP modélisé pris comme origine représente la correction de distance ΔH, est prise comme position du CRP corrigé.
Le traitement selon l'invention est l'équivalent d'une migration avant stack.
L'invention sera expliquée à présent en regard des dessins annexés dans lesquels :
La figure 2 représente, dans un macromodèle, un trajet de rayon sismique relatif à un CRP modélisé et un trajet de rayon réel associé au même couple source/récepteur ; La figure 3 représente schématiquement le cas où toutes les corrections de profondeur coïncident ; et
La figure 4 représente le cas où les corrections de profondeur sont dispersées.
Avec référence tout d'abord à la figure 2, le macromodèle obtenu par inversion tomographique, comme expliqué précédemment, comprend plusieurs couches limitées entre elles par des interfaces modélisées telles que i. Le macromodèle est paramétré par la profondeur z^(x) de chaque point d'une interface i et par la vitesse de propagation V^(x) de l'onde sismique dans la couche géologique surmontant l'interface i. Les grandeurs z^(x) et Vj_(x) sont, par exemple, des fonctions de type fonction spline de la coordonnée horizontale x .
Si besoin est, l'ensemble des écarts Δt correspondant à tous les CRP peut être minimisé en faisant appel à un algorithme d'inversion. A la convergence, il peut arriver que pour chaque rayon sismique, l'écart résiduel ne soit pas nul, ce qui peut s'expliquer notamment par la simplification de la structure réelle du sous-sol résultant de sa paramétrisation.
L'invention consiste à- utiliser ledit écart résiduel sur le temps de trajet de chaque rayon pour approcher l'emplacement exact du réflecteur. A cet effet, pour tous les rayons correspondant à un même CRP, on convertit les écarts résiduels de temps de trajet Δt en corrections de distance Δh ou en leurs composantes corrections de profondeur Δz et corrections de position horizontale Δx et éventuellement Δ y.
Sur la figure 2, on a représenté cette correction de profondeur pour le trajet particulier 12 qui part de la source S, se réfléchit sur le CRP considéré 13 et aboutit au récepteur R situé à une distance -d de la source. Cette correction est portée sur la normale 24 au réflecteur modélisé jL au CRP 13 considéré. La référence 14 désigne le trajet correspondant à un rayon qui part de S, se réfléchit sur le CRP corrigé 15 et aboutit en R.
Le procédé selon l'invention est illustré par les figures 3 et 4. Sur l'axe 24 correspondant à la normale au réflecteur _i au CRP considéré 13, on porte, à partir d'une origine représentée par ledit CRP 13 du macromodèle, les corrections de distance Δh correspondant à tous les rayons qui illuminent ledit CRP lorsque <1 varie, chaque correction Δh étant matérialisée par un trait 16 perpen¬ diculaire à la normale 24 et représentatif d'une impulsion.
La figure 3 montre le cas où toutes les corrections de profondeur 16 coïncident. L'impulsion correspondant à chaque correction est convoluée par un signal, par exemple du type ondelette de Ricker. Tous les signaux étant en coïncidence, le résultat de la convolution est est un signal 18 dont le maximum d'amplitude est égal à la somme des amplitudes des impulsions.
La figure 4 correspond au cas où les corrections de distance Δh sont étalées. Ici aussi, les corrections sont convoluées par un signal donné, par exemple du type ondelette de Ricker, les signaux 18,20 résultant de la convolution étant ensuite additionnés. La dispersion des corrections se traduit par une amplitude de la somme moindre que dans le cas précédent.
L'avantage du procédé selon l'invention est qu'il permet d'obtenir une image synthétique du réflecteur plus proche de la réalité. De plus, des hautes fréquences qui ont été perdues dans la représentation du macromodèle par des fonctions splines sont récupérées.

Claims

REVENDICATIONS 1- Procédé d'estimation améliorée de la position de réflecteurs acoustiques présents dans un milieu, par interprétation des écarts entre les temps de trajet réels et modélisés des ondes acoustiques dans ce milieu, du type consistant à : collecter des traces sismiques au moyen d'un système d'enregistrement comprenant au moins une source (S) émettrice d'ondes acoustiques et plusieurs récepteurs (R) que l'on déplace à la surface du milieu,
- déterminer à partir de ces traces les temps de trajets réels observés ( 0^g) pour chaque couple source/ récepteur et correspondant à des réflexions sur les différents réflecteurs du milieu, - définir un macromodèle paramétré donnant la profondeur jz des réflecteurs (i) et la vitesse V des couches du milieu en fonction de la coordonnée horizontale x et éventuellement de la coordonnée horizontale ,
- déterminer, à l'aide du macromodèle, des temps de trajet source/récepteur modélisés (tm0(j) par la technique du tracer de rayons, et
- calculer, pour chaque réflecteur et pour chaque rayon, l'écart de temps de trajet Δt ≈ tQ^g - tmQ(- , ledit procédé étant caractérisé en ce qu'il consiste à :
- effectuer le tracer de rayons à partir de points, dits CRP, choisis sur les réflecteurs du macromodèle,
- transformer, pour chacun desdits CRP, l'ensemble des écarts Δt obtenus pour tous les rayons du macromodèle issus du CRP concerné en une correction de distance ΔH le long d'une direction prédéterminée (24) passant par ledit CRP, et
- déplacer le CRP concerné du macromodèle de la distance ΔH le long de ladite direction prédéterminée pour fournir une position corrigée dudit CRP, l'ensemble des positions corrigées des CRP d'un réflecteur du macroraodèle représentant une position corrigée améliorée de ce réflecteur. 2- Procédé selon la revendication 1, caractérisé en ce que la transformation de l'ensemble des Δt relatifs à un même CRP en la correction correspondante de distance ΔH consiste à: - représenter une même impulsion (16) perpendiculairement à la direction prédéterminée (24) en chacun des points de ladite direction dont les distances à une origine, représentée par la position du CRP concerné déduite du macromodèle, ont des valeur Δh représentant des corrections de distance correspondant aux Δ t associés aux différents rayons illuminant le CRP lorsque la distance source/récepteur ou offset varie,
- convoluer chacune desdites impulsions par un même signal donné (18;20), - additionner les signaux résultant de la convolution pour obtenir un signal somme, et
- prendre comme valeur Δ H la distance séparant la position du maximum dudit signal somme de l'origine choisie. 3- Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que la direction prédéterminée passant par le CRP concerné est la normale (24) au réflecteur contenant ledit CRP.
4- Procédé selon l'une des revendications 1 à 3, caractérisé en ce que la profondeur z, des réflecteurs et la vitesse V des couches du modèle sont des fonctions splines de la coordonnée horizontale x.
5- Procédé selon l'une des revendications 1 à 4, caractérisé en ce que l'on modifie le macromodèle par itération en mettant en oeuvre un algorithme d'inversion jusqu'à ce que tous les écarts t aient convergé vers une valeur stable minimale.
6- Procédé selon l'une des revendications 1 à 5, caractérisé en ce que pour tracer les rayons, on effectue un échantillonnage régulier du réflecteur sur le macroraodèle, on se fixe une première distance source/récepteur ou offset, on trace les rayons correspondant à t'ous les CRP, puis on se donne un second offset et l'on trace les rayons pour les mêmes CRP, et ainsi de suite avec plusieurs offsets.
PCT/FR1993/001040 1992-10-26 1993-10-25 Procede d'estimation de la position de reflecteurs acoustiques presents dans un milieu WO1994010584A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR9212729A FR2697347B1 (fr) 1992-10-26 1992-10-26 Procédé d'estimation améliorée de la position de réflecteurs acoustiques présents dans un milieu, par interprétation des écarts entre les temps de trajet réels et modélisés des ondes acoustiques dans ce milieu.
FR92/12729 1992-10-26

Publications (1)

Publication Number Publication Date
WO1994010584A1 true WO1994010584A1 (fr) 1994-05-11

Family

ID=9434845

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/FR1993/001040 WO1994010584A1 (fr) 1992-10-26 1993-10-25 Procede d'estimation de la position de reflecteurs acoustiques presents dans un milieu

Country Status (2)

Country Link
FR (1) FR2697347B1 (fr)
WO (1) WO1994010584A1 (fr)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2717270B1 (fr) * 1994-03-11 1996-04-26 Elf Aquitaine Procédé de modélisation de données sismiques cinématiques ayant subi un traitement par au moins un opérateur de déplacement.
US6128581A (en) * 1999-05-19 2000-10-03 Pgs Seres As Dynamic datumming for land and marine multicomponent seismic data processing

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2226884A (en) * 1989-01-06 1990-07-11 Marathon Oil Co Model-base depth processing of seimic data

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2226884A (en) * 1989-01-06 1990-07-11 Marathon Oil Co Model-base depth processing of seimic data

Also Published As

Publication number Publication date
FR2697347A1 (fr) 1994-04-29
FR2697347B1 (fr) 1995-01-20

Similar Documents

Publication Publication Date Title
US6654693B2 (en) Angle dependent surface multiple attenuation for two-component marine bottom sensor data
EP0889331B1 (fr) Méthode pour modéliser en 3D l&#39;impédance d&#39;un milieu hétérogène
EP0354112B1 (fr) Méthode pour obtenir un modèle représentatif d&#39;un milieu hétérogène et notamment du sous-sol
EP2386880B1 (fr) Prédiction rapide de multiples surfaces 3D
CA2528700C (fr) Methode pour determiner des informations speculaires apres imagerie sismique avant sommation
US9322943B2 (en) Method and apparatus for pre-stack deghosting of seismic data
EP0591037B1 (fr) Méthode d&#39;acquisition et de traitement de données pour suivre le déplacement des fluides dans un réservoir
NO337140B1 (no) Fremgangsmåte for bestemmelse av dekningsgrad ved marin seismisk kartlegging
WO1985005696A1 (fr) Procede de diagraphie acoustique instantanee dans un puits de forage
EP0786670A1 (fr) Procédé de traitement de calibration d&#39;une paire de capteurs hydrophone/géophone et procédé de prospection sismique mettant en oeuvre ce traitement
JP2000512385A (ja) 伝搬する波動場のサンプリング及び復元
EP0786671B1 (fr) Procédé de prospection sismique marine au moyen d&#39;un couple de capteurs hydrophone et géophone
EP0426514B1 (fr) Procédé de prospection sismique marine
EP0698220B1 (fr) Procede de modelisation de donnees sismiques cinematiques ayant subi un traitement par au moins un operateur de deplacement
AU739128B2 (en) A method of seismic processing, and in particular a 3D seismic prospection method implementing seismic data migration
WO1994010584A1 (fr) Procede d&#39;estimation de la position de reflecteurs acoustiques presents dans un milieu
AU2013201068B2 (en) Determination of particle displacement or velocity from particle acceleration measurements
EP0909961B1 (fr) Procédé de séparation d&#39;ondes en sismique de puits pour des acquisitions de type à déport croissant
CN109490964B (zh) 一种改进的高精度avo弹性参数快速反演方法
GB2375606A (en) Angle dependent surface multiple attenuation for two - component marine bottom sensor data
EP0610123B1 (fr) Procédé pour améliorer l&#39;estimation, par analyse de focalisation, des vitesses de propagation des ondes sismiques
CN113614578A (zh) 使用双向虚拟源重建的地震勘测
CA2085617A1 (fr) Procede de traitement pour l&#39;obtention d&#39;une section somme a offsets nuls
Lecomte Illumination, resolution, and incidence-angle in PSDM: a tutorial.
Rosales et al. Converted-waves angle-domain common-image gathers

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): CA GB NO US

NENP Non-entry into the national phase

Ref country code: CA