WO1998003886A1 - Methode de migration des attributs geophysiques d'un milieu - Google Patents

Methode de migration des attributs geophysiques d'un milieu Download PDF

Info

Publication number
WO1998003886A1
WO1998003886A1 PCT/FR1997/001247 FR9701247W WO9803886A1 WO 1998003886 A1 WO1998003886 A1 WO 1998003886A1 FR 9701247 W FR9701247 W FR 9701247W WO 9803886 A1 WO9803886 A1 WO 9803886A1
Authority
WO
WIPO (PCT)
Prior art keywords
migrated
speed
migration
amplitude
iteration
Prior art date
Application number
PCT/FR1997/001247
Other languages
English (en)
Inventor
Jean-Pierre Dunand
Eric De Bazelaire
Original Assignee
Elf Exploration Production
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 Elf Exploration Production filed Critical Elf Exploration Production
Priority to US09/043,217 priority Critical patent/US6044039A/en
Priority to AU36252/97A priority patent/AU712392B2/en
Priority to EP97932870A priority patent/EP0888562A1/fr
Priority to BR9706587A priority patent/BR9706587A/pt
Publication of WO1998003886A1 publication Critical patent/WO1998003886A1/fr
Priority to NO980938A priority patent/NO980938L/no

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
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/52—Move-out correction

Definitions

  • the present invention relates to a method of migrating the geophysical attributes of a medium and more particularly of a given area of said medium.
  • Seismic reflection is a common method which makes it possible, in petroleum exploration in particular, to obtain a seismic image of the medium to be explored.
  • acoustic waves are emitted which propagate in the medium to be explored and are reflected on the different reflectors or horizons that it contains.
  • the reflected waves are recorded, as a function of time, on suitable receivers placed on the surface of the medium.
  • the records or traces are then gathered according to a criterion determined according to the objectives to be achieved and constitute collections of traces.
  • the '' traces, common midpoint (CMP) gathers of traces sets that are assigned to the positions of points in the center of the segments couples transceiver used for said traces.
  • Seismic processing makes it possible, from said collections of traces, to obtain a seismic image in the vertical plane passing through all of the midpoints.
  • the reflections of the waves on the various reflectors observed on a collection of traces in common middle point, theoretically align along hyperbolas centered vertically from the midpoint and called indicators.
  • said traces are dynamically corrected with a law or velocity field V (t).
  • This velocity field is to punctually analyze the velocities on a limited number of collections of traces in common midpoint, then to interpolate the results from these analyzes, in time on the one hand for each of the analyzes. and on the x-axis on the other hand.
  • Conventional speed analysis consists in successively applying to PMC trace collections, for the selected midpoints, constant speeds then summing the dynamically corrected traces for each of the speeds used and manually retaining those of speeds leading to a maximum of energy of the sum trace.
  • POLYSTACK In an article entitled “Normal Moveout Revisited: Inhomogeneous media and curved interfaces", published in the journal GEOPHYSICS, Volume 53, N ° 2, Fev. 1988, pages 143 to 157, Eric de Bazelaire developed another method of speed analysis, used to obtain improved sum sections called "POLYSTACK".
  • the POLYSTACK method consists in scanning all the traces of the same collection of PMC traces, in applying to all the traces a correction of the static type, according to a family of hyperbolas, independent of time, and different from a trace. sum to the other, so as to produce "baps" of which each trace is a sum of the traces of the PMC collection thus corrected.
  • the two methods POLYSTACK and DELTA STACK have the particularity of comprising and / or making more accessible the geophysical attributes of the zone of the explored medium including obviously the attributes amplitude and speed.
  • the sum section obtained by the various methods recalled above contains a number of anomalies (diffraction hyperbola, poor localization of slope events, etc.) which should be corrected.
  • the sum section can be migrated in time or in depth depending on the goals to be achieved, it being specified that one can go from a time migration to the deep migration and vice versa. This is well known to specialists and will not be described in detail, the migration having the effect of replacing the seismic event or events appearing on the sum section vertically from their geological position.
  • To migrate a seismic signal it is necessary to know the distribution of the velocities in the subsoil, which, in geophysical terms, amounts to knowing the velocity field of intervals in the migrated position.
  • this interval velocity field is a priori unknown and whatever the migration technique used, wave equation or ray tracing, it is replaced by a speed distribution law from the stack or any other form modeling. Any error concerning this field of velocities introduces a bias on the result of the migration of the amplitude in particular.
  • the present invention aims to propose a method which allows the migration of geophysical attributes of an area to be explored or at least geophysical attributes of interest in petroleum exploration.
  • the present invention relates to the determination of a speed field in the migrated position using a migration of the amplitude attribute and the use of the migrated speed field for the migration of any other geophysical attribute.
  • the subject of the present invention is a method in which:
  • - at least one collection of seismic traces is classified, classified according to a determined criterion and relating to said zone, - a sum section is produced from said collection of seismic traces and using a summation speed field relating to said zone, characterized in that, in addition, a) the product of a term relating to the amplitude of the traces is produced by a term relating to the summation speeds; b) migrating with an initial migration speed, on the one hand, the amplitude data to obtain a migrated image of the amplitude and, on the other hand, said product to obtain a migrated image of said product, c) performing the ratio of the data of the migrated image of said product to the data of the migrated image of the amplitude, in order to obtain a speed in the migrated position, and d) an iteration of steps b) and c) is carried out using for each iteration a speed of migration equal to the speed in the migrated position from the immediately preceding iteration until the difference between the
  • the present invention makes it possible to obtain a speed field very close to the migrated speed field from which a correct migration of the other geophysical attributes can be carried out.
  • the migration is a migration in time.
  • step c) is carried out on the envelopes of the traces after time migration of the amplitude and of the product.
  • the term relating to the summation speeds is assigned an exponent less than or equal to 1.
  • each of the other geophysical attributes is migrated in time with the final velocity field by performing steps a) to c) only.
  • the initial migration speed is a summation speed.
  • FIG. 1 is a synthetic representation of a collection of traces in common midpoints
  • FIG. 2 is a synthetic representation of a geological model concerning an area considered to be explored
  • FIG. 3 is a synthetic representation of a sum model
  • FIG. 4 is a synthetic representation of a field of speeds before migration
  • FIG. 5 is a synthetic representation of the estimate of the speed field after the first iteration
  • FIGS. 6 to 8 are synthetic representations of the estimates of the speed field after the second, third and fourth iterations respectively
  • FIG. 9 is a synthetic representation of the migration of the amplitude with the speed field of Figure 8
  • - Figure 10 is a synthetic representation of the migration of the amplitude with conventional migration techniques.
  • the initial geological model of the area concerned is shown synthetically in Figure 2, the ordinates representing the depths while in the other figures, the ordinates represent times.
  • the geological model we can point a certain number of reflectors H j to H 6 located at depths located between 490 m for the reflector H j and 2060 m for the reflector H 6 .
  • the geological model is by nature almost perfect since it is carried out manually by the geophysicist. From the geological model, we synthesize collections of PMC traces then we perform velocity analyzes, and we perform the summation of traces of PMC collections to obtain a sum model or stack model like the one represented in figure 3, on which we finds the same reflectors Hj to H 6 at the corresponding times.
  • the stack model in Figure 3 is representative of the amplitude stack, which constitutes a geophysical attribute that one wishes to migrate in time, for example.
  • any initial speed field such as for example the speed field used for the summation or else the speed field obtained by the POLYSTACK or DELTA STACK processes mentioned in the preamble.
  • An initial speed field usable for the migration of the attribute amplitude can be that represented on figure 4.
  • this field of speeds one notes various gradations of tones which go from black to white with gray and which correspond to isovitesses whose values can be determined using the palette illustrated on the left in figure 4.
  • the method according to the invention consists, in a first step, in carrying out the product of the term amplitude A with a term representative of the speed of stack V.
  • the product A * V being the simplest combination but any invertible function of A and V could be used to make the combination.
  • this first step consists in sampling each sum trace of the stack model and in performing the product of the amplitude of each sample by the corresponding stack speed.
  • the term V is weighted by an exponent ⁇ which is at most equal to 1, in order to reduce the parasitic effects of migration.
  • the exponent ⁇ is equal to 0.5.
  • the second operation consists in migrating on the one hand the amplitude trace A and, on the other hand, the product A * V ⁇ in the same interval velocity field originating from the initial velocity field used which, as recalled above , may be the stack velocity field.
  • the journey times are calculated using, for example, the algorithm known as EIKONAL, because the calculation is unique for the two migrations A and A * V ⁇ .
  • EIKONAL the algorithm known as EIKONAL
  • any other time migration can be used, such as wave equation migration, a technique that specialists are well aware of.
  • the previous time migration thus produces a time migrated image corresponding to the amplitude and a time migrated image corresponding to the product A * V ⁇ .
  • the third step consists in carrying out the ratio of the migrated image (A * V ⁇ ) on the migrated image A, so that a first speed is obtained in the migrated position V ⁇ because:
  • the present invention also consists in making the migrated speed field converge by successive iterations.
  • (A) migrates to obtain at the end of iteration 1 a speed in the migrated position V 2 . This is what is shown in Figure 5 where we can still see long black spots.
  • the speed in migration position V 3 is used as the migration speed, the resulting migrated speed field being represented in FIG. 7.
  • the speed field of FIG. 7 converges satisfactorily.
  • FIG. 9 represents the migrated image of the amplitude with the iteration speed field 4.
  • this migrated image is compared with a migrated image of the same amplitude but performed with the prior techniques and represented in FIG. 10, easily that the reflectors H and H 3 of FIG. 9 are more curved after their apex than the same reflectors of FIG. 10, which further defines the fault which is located by the right ends of said reflectors.
  • the last reflector H 7 downwards in FIGS. 9 and 10 is more marked in FIG. 9 than in FIG. 10.
  • the present invention also makes it possible to migrate in time other geophysical attributes A ⁇ . such as signal-to-noise ratio, accuracy on curvature, etc.
  • Another advantage of the present invention lies in the fact that the spatial distribution of the interval speeds is better known, which makes it possible to know the true speeds in the migrated section and the knowledge of these true speeds participates in the determination of the lithology of the basement for example.

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

Elle est du type dans laquelle on constitue au moins une collection de traces sismiques, classées selon un critère déterminé et concernant ladite zone; on réalise une section somme à partir de ladite collection de traces sismiques et à l'aide d'un champ de vitesses de sommation concernant ladite zone, caractérisée en ce que, en outre: a) on effectue le produit d'un terme relatif à l'amplitude des traces par un terme relatif aux vitesses de sommation; b) on migre avec une vitesse de migration initiale, d'une part, les données amplitude pour obtenir une image migrée de l'amplitude et, d'autre part, ledit produit pour obtenir une image migrée dudit produit; c) on effectue le rapport des données de l'image migrée dudit produit sur les données de l'image migrée de l'amplitude, pour obtenir une vitesse en position migrée, et d) on procède à une itération des étapes b) et c) en utilisant pour chaque itération une vitesse de migration égale à la vitesse en position migrée issue de l'itération immédiatement précédente jusqu'à ce que l'écart entre la vitesse en position migrée finale d'itération (n) et la vitesse en position migrée d'itération (n-1) soit inférieur à un seuil prédéterminé.

Description

Méthode de migration des attributs géophysiques d'un milieu
La présente invention concerne une méthode de migration des attributs géophysiques d' un milieu et plus particulièrement d' une zone donnée dudit milieu.
La sismique réflexion est une méthode usuelle qui permet, en exploration pétrolière notamment, d'obtenir une image sismique du milieu à explorer. Dans cette méthode, à partir de sources d'énergie appropriées, appelées émetteurs, on émet des ondes acoustiques qui se propagent dans le milieu à explorer et se réfléchissent sur les différents réflecteurs ou horizons qu' il renferme. Les ondes réfléchies sont enregistrées, en fonction du temps, sur des récepteurs adéquats, disposés à la surface du milieu. Les enregistrements ou traces sont ensuite rassemblés en fonction d' un critère déterminé en fonction des objectifs à atteindre et constituent des collections de traces. La collection ' de traces, en point milieu commun (PMC) rassemble les séries de traces qui sont affectées aux positions des points situés au milieu des segments des couples émetteur-récepteur utilisés pour lesdites traces.
Le traitement sismique permet, à partir desdites collections de traces d'obtenir une image sismique dans le plan vertical passant par l'ensemble des points milieux. Lorsqu'on fait l'hypothèse d'un milieu homogène et isotrope, en couches planes et parallèles, les réflexions des ondes sur les divers réflecteurs, observées sur une collection de traces en point milieu commun, s'alignent théoriquement le long d'hyperboles centrées à la verticale du point milieu et appelées indicatrices. De manière à réaliser la sommation (stack en anglais) des traces de chaque collection, on corrige dynamiquement lesdites traces avec une loi ou champ de vitesses V(t). Une manière pour obtenir ce champ de vitesses, consiste à effectuer ponctuellement des analyses de vitesses sur un nombre limité de collections de traces en point milieu commun, puis à interpoler les résultats issus de ces analyses, en temps d'une part pour chacune des analyses et en abscisse d'autre part. L'analyse de vitesse classique consiste à appliquer successivement aux collections de traces en PMC, pour les points milieux sélectionnés, des vitesses constantes puis à sommer les traces corrigées dynamiquement pour chacune des vitesses utilisées et à retenir manuellement celles des vitesses conduisant à un maximum d'énergie de la trace somme.
Dans un article intitulé "Normal Moveout Revisited : Inhomogeneous média and curved interfaces", publié dans la revue GEOPHYSICS, Volume 53, N° 2, Fev. 1988, pages 143 à 157, Eric de Bazelaire a développé une autre méthode d'analyse de vitesses, utilisée pour l'obtention de sections somme améliorées et appelées "POLYSTACK" . Très brièvement, la méthode POLYSTACK consiste à balayer toutes les traces d'une même collection de traces PMC, à appliquer à toutes les traces une correction du type statique, selon une famille d'hyperboles, indépendante du temps, et différente d'une trace somme à l'autre, de manière à produire des "bap" dont chaque trace est une somme des traces de la collection PMC ainsi corrigée.
Une autre méthode plus élaborée que le POLYSTACK est décrite dans la demande FR-A-2 726 091. Cette autre méthode, connue sous la dénomination DELTA STACK permet, comme précédemment, l'obtention, d'une part d'un champ fin de vitesses et de meilleure résolution et, d'autre part, d'une section somme améliorée et ce, avec des "BAP" , du type POLYSTACK mais de taille plus réduite.
Les deux méthodes POLYSTACK et DELTA STACK ont la particularité de comporter et/ou de rendre plus accessibles les attributs géophysiques de la zone du milieu explorée y compris bien évidemment les attributs amplitude et vitesse.
La section somme qu'on obtient par les diverses méthodes rappelées ci-dessus comporte un certain nombre d'anomalies (hyperbole de diffraction, mauvaise localisation des événements pentes, etc.) qu'il convient de corriger.
Lorsque la zone du milieu à explorer est isotrope, homogène, en couches planes et parallèles, les vitesses de tranche ou d'intervalle peuvent être calculées facilement par la formule de DIX également bien connue des spécialistes. Dans le cas d'une zone comprenant des réflecteurs pentes, il est nécessaire de procéder à une correction préalable, avant la migration, et connue sous l'appellation DMO (DIP MOVE OUT).
La section somme peut être migrée en temps ou en profondeur en fonction des buts à atteindre, étant précisé que l'on peut passer d'une migration temps à la migration profondeur et inversement. Ceci est bien connu des spécialistes et ne sera pas décrit dans le détail, la migration ayant pour effet de replacer le ou les événements sismiques apparaissant sur la section somme à la verticale de leur position géologique. Pour migrer un signal sismique il est nécessaire de connaître la répartition des vitesses dans le sous-sol, ce qui, en termes géophysiques, revient à connaître le champ de vitesses d'intervalles en position migrée.
En fait, ce champ de vitesses d'intervalles est a priori inconnu et quelle que soit la technique de migration utilisée, équation des ondes ou lancer de rayon, on le remplace par une loi de répartition des vitesses issue du stack ou de toute autre forme de modélisation. Toute erreur concernant ce champ de vitesses introduit un biais sur le résultat de la migration de l'amplitude notamment.
La présente invention a pour but de proposer une méthode qui permet la migration des attributs géophysiques d'une zone à explorer ou tout au moins des attributs géophysiques présentant un intérêt dans l'exploration pétrolière.
Plus particulièrement, la présente invention a pour objet la détermination d'un champ de vitesses en position migrée en utilisant une migration de l'attribut amplitude et l'utilisation du champ de vitesses migré pour la migration de toute autre attribut géophysique.
La présente invention a pour objet une méthode dans laquelle :
- on constitue au moins une collection de traces sismiques, classées selon un critère déterminé et concernant ladite zone, - on réalise une section somme à partir de ladite collection de traces sismiques et à l'aide d'un champ de vitesses de sommation concernant ladite zone, caractérisée en ce que, en outre, a) on effectue le produit d'un terme relatif à l'amplitude des traces par un terme relatif aux vitesses de sommation ; b) on migre avec une vitesse de migration initiale, d'une part, les données amplitude pour obtenir une image migrée de l' amplitude et, d'autre part, ledit produit pour obtenir une image migrée dudit produit, c) on effectue le rapport des données de l'image migrée dudit produit sur les données de l'image migrée de l'amplitude, pour obtenir une vitesse en position migrée, et d) on procède à une itération des étapes b) et c) en utilisant pour chaque itération une vitesse de migration égale à la vitesse en position migrée issue de l'itération immédiatement précédente jusqu'à ce que l'écart entre la vitesse en position migrée finale d'itération (n) et la vitesse en position migrée d'itération (n-1) soit inférieur à un seuil prédéterminé.
La présente invention permet d'obtenir un champ de vitesse très proche du champ de vitesse migré à partir duquel on peut effectuer une migration correcte des autres attributs géophysiques. Selon une autre caractéristique de l'invention, la migration est une migration en temps.
Selon une autre caractéristique de l'invention, l'étape c) est effectuée sur les enveloppes des traces après migration temps de l'amplitude et du produit. Selon une autre caractéristique de l'invention, le terme relatif aux vitesses de sommation est affecté d'un exposant inférieur ou égal à 1.
Selon une autre caractéristique de l'invention, chacun des autres attributs géophysiques est migré en temps avec le champ de vitesses final en effectuant les étapes a) à c) uniquement. Selon une autre caractéristique de l'invention, la vitesse de migration initiale est une vitesse de sommation.
D'autres avantages et caractéristiques ressortiront mieux à la lecture d'un mode de réalisation de l'invention, ainsi que des dessins annexés sur lesquels : - la figure 1 est une représentation synthétique d'une collection de traces en points milieux communs,
- la figure 2 est une représentation synthétique d'un modèle géologique concernant une zone considérée à explorer,
- la figure 3 est une représentation synthétique d'un modèle somme, - la figure 4 est une représentation synthétique d'un champ de vitesses avant migration,
- la figure 5 est une représentation synthétique de l'estimation du champ de vitesses après la première itération, - les figures 6 à 8 sont des représentations synthétiques des estimations du champ de vitesses après respectivement les deuxième, troisième et quatrième itérations,
- la figure 9 est une représentation synthétique de la migration de l'amplitude avec le champ de vitesses de la figure 8, - la figure 10 est une représentation synthétique de la migration de l'amplitude avec des techniques classiques de migration.
Par les techniques habituelles et utilisées en sismique réflexion, on constitue notamment ce qu'on appelle des collections de traces en points milieux communs. Une collection de traces en points milieux communs (collection PMC) est représentée de manière synthétique sur la figure 1 et concerne une zone géologique d'un milieu à explorer.
Le modèle géologique initial de la zone concernée est représenté de façon synthétique sur la figure 2, les ordonnées représentant les profondeurs alors que sur les autres figures, les ordonnées représentent des temps. Sur le modèle géologique, on peut pointer un certain nombre de réflecteurs Hj à H6 situés à des profondeurs situées entre 490 m pour le réflecteur Hj et 2060 m pour le réflecteur H6. Le modèle géologique est par nature quasi parfait puisque réalisé manuellement par le géophysicien. A partir du modèle géologique, on synthétise des collections de traces PMC puis on effectue des analyses de vitesses, et on effectue la sommation des traces des collections PMC pour obtenir un modèle somme ou modèle stack comme celui représenté sur la figure 3, sur lequel on retrouve les mêmes réflecteurs Hj à H6 aux temps correspondants.
Le modèle stack de la figure 3 est représentatif du stack amplitude, laquelle constitue un attribut géophysique qu'on souhaite migrer en temps par exemple. Pour ce faire, on utilise car cela est nécessaire, un champ de vitesses initial quelconque, comme par exemple le champ de vitesses utilisé pour la sommation ou encore le champ de vitesses obtenu par les procédés POLYSTACK ou DELTA STACK rappelés dans le préambule. Un champ de vitesses initial utilisable pour la migration de l'attribut amplitude peut être celui représenté sur la figure 4. Sur ce champ de vitesses, on constate différentes gradations de tons qui vont du noir au blanc avec des gris et qui correspondent à des isovitesses dont les valeurs peuvent être déterminées à l'aide de la palette illustrée à gauche sur la figure 4. C'est ainsi que les gris, en bas et à gauche sur la figure, correspondent à des vitesses de l'ordre de 5000 m/s, la première tranche de blanc située au- dessus du temps double de 2,5 secondes correspondant à la vitesse de 4686 m/s, et ainsi de suite. En abscisse on indique le numéro d'ordre des traces PMC utilisées pour réaliser le modèle. Le procédé selon l'invention consiste, dans une première étape, à effectuer le produit du terme amplitude A avec un terme représentatif de la vitesse de stack V. Bien entendu, toute autre combinaison de ces deux termes, pour autant qu'elle soit cohérente, pourrait être utilisée, le produit A*V étant la combinaison la plus simple mais toute fonction inversible de A et V pourrait être utilisée pour réaliser la combinaison. En fait, cette première étape consiste à échantillonner chaque trace somme du modèle stack et à effectuer le produit de l'amplitude de chaque échantillon par la vitesse de stack correspondante.
Comme le champ de vitesses peut présenter un spectre très large, par exemple entre 1000 et 10000 m/s ou entre 2490 et 5000 m/s pour le champ de vitesses de la figure 4, on pondère le terme V par un exposant α qui est au plus égal à 1 , afin d'atténuer les effets parasites de la migration. De préférence, l'exposant α est égal à 0,5. De cette manière, on obtient une grandeur fixe pour la suite du processus et définie par le produit A*Vα. La deuxième opération consiste à migrer d'une part la trace amplitude A et, d'autre part, le produit A*Vα dans le même champ de vitesses d'intervalle issu du champ de vitesses initial utilisé qui, comme rappelé ci-dessus, peut être le champ de vitesses stack. Pour optimiser cette étape, on préfère réaliser une migration du type KIRCHHOFF dans laquelle on calcule les temps de trajet en utilisant, par exemple, l'algorithme connu sous le nom d'EIKONAL, car le calcul est unique pour les deux migrations A et A*Vα. Cependant, toute autre migration temps peut être utilisée, comme la migration par équation des ondes, technique que les spécialistes connaissent bien. La migration temps précédente produit ainsi une image migrée temps correspondant à l'amplitude et une image migrée temps correspondant au produit A*Vα.
La troisième étape consiste à effectuer le rapport de l'image migrée (A*Vα) sur l'image migrée A, de sorte qu'on obtient une première vitesse en position migrée Vα car :
(A *Vα) migré a . ,
- = V migre
A migre
Pour éviter les passages au zéro de la trace amplitude stack lors de la migration, ce qui est toujours très difficile à réaliser, il est préférable d'effectuer le rapport précédent sur les enveloppes des traces, après migration, de l'amplitude de A et du produit A*Vα.
De cette manière, on obtient un premier champ de vitesses en position migrée. Toutefois, la migration précédente est imparfaite car elle est réalisée à partir d'un champ de vitesses notamment de stack en position non migrée. Aussi, la présente invention consiste, en outre, à faire converger le champ de vitesses migré par des itérations successives.
A l'itération 1, on répète les deuxième et troisième étapes mais avec la première vitesse en position migrée qu'on désigne par Vi . Ainsi, on migre en temps, avec la vitesse Vt , le terme amplitude A et le produit A*Vα, puis on procède au rapport
(A * vα) migré a . A — r — = V* migré
(A) migre pour obtenir à l'issue de l'itération 1 une vitesse en position migrée V2. C'est ce qui est représenté sur la figure 5 où l'on peut voir encore de longues tâches noires.
A l'itération 2, on utilise la vitesse de migration V2 pour la migration du terme amplitude A et du produit A*Vα, le rapport de ces deux termes donnant le champ de vitesses V3 issu de l'itération 2 qui est représenté sur la figure 6 qui est nettement améliorée par rapport à la figure 5.
A l'itération 3, on utilise comme vitesse de migration, la vitesse en position migrée V3, le champ de vitesses migré résultant étant représenté sur la figure 7. Le champ de vitesses de la figure 7 converge de façon satisfaisante.
A l'issue de l'itération 4, le champ de vitesses correspondant et représenté sur la figure 8, a totalement convergé, la différence avec celui de la figure 7 étant minime. Il est inutile, à ce stade du processus, de poursuivre l'itération car les champs de vitesses qu'on pourrait obtenir ne seraient pas beaucoup plus précis que celui représenté sur la figure 8. En fait, on arrête le processus itératif lorsque l'écart entre le champ de vitesses migré issu de l'itération n (ici itération 4) et le champ de vitesses migré issu de l'itération n-1 (ici itération 3) est inférieur à un certain seuil prédéterminé, la limite étant telle que la conversion en vitesses d'intervalles du champ de vitesses de stack en position migrée est identique au champ de vitesses d'intervalles utilisé pour effectuer la migration.
Il faut noter qu'on est bien en présence d'un processus convergeant automatiquement vers le vrai champ de vitesses d'intervalle du sous-sol d'un milieu au sens de DIX;
La figure 9 représente l'image migrée de l'amplitude avec le champ de vitesses itération 4. Lorsqu'on compare cette image migrée avec une image migrée de la même amplitude mais effectuée avec les techniques antérieures et représentée sur la figure 10, on constate aisément que les réflecteurs H et H3 de la figure 9 sont davantages courbés après leur sommet que les mêmes réflecteurs de la figure 10, ce qui définit davantage la faille qui est localisée par les extrémités à droite desdits réflecteurs. En outre, le dernier réflecteur H7 vers le bas sur les figures 9 et 10, est davantage marqué sur la figure 9 que sur la figure 10.
La présente invention permet également de migrer en temps d'autres attributs géophysiques A →. tels que rapport signal/bruit, précision sur la courbure, etc ..
En effet, connaissant le champ de vitesses en position migrée il suffit d'effectuer le produit A * A -., de migrer en temps le terme amplitude
A et le produit A * A^ avec le champ de vitesses en position migrée, tel que celui de l'itération 4 qui constitue en fait le champ de vitesses final, et de calculer le rapport ^A A-j- migre
A migre
pour avoir A^-, en position migrée et ce, sans itération puisque le champ de vitesses de migration est maintenant connu. II va de soi que le processus itératif décrit ci-dessus ne converge que dans la bande spatiale du champ de vitesses de stack en position non migrée. L'utilisation d'un champ de vitesses interpolées ne permettrait de définir que les basses fréquences du champ de vitesses réel. C'est la raison pour laquelle, l'invention utilise, de préférence, les techniques POLYSTACK et DELTA STACK qui permettent d'obtenir des sections amplitude et vitesses dont les bandes spatiales sont identiques, dans le but de réaliser des sections POLYSTACK ou DELTA STACK en position migrée.
Un autre avantage de la présente invention réside dans le fait qu'on connaît mieux la répartition spatiale des vitesses d'intervalles, ce qui permet de connaître les vitesses -vraies dans la section migrée et la connaissance de ces vitesses vraies participe à la détermination de la lithologie du sous-sol par exemple.

Claims

REVENDICATIONS
1. Méthode de migration des attributs géophysiques d'une zone d'un milieu à explorer dans laquelle :
- on constitue au moins une collection de traces sismiques, classées selon un critère déterminé et concernant ladite zone,
- on réalise une section somme à partir de ladite collection de traces sismiques et à l'aide d'un champ de vitesses de sommation concernant ladite zone, caractérisée en ce que, en outre, a) on effectue le produit d'un terme relatif à l'amplitude des traces par un terme relatif aux vitesses de sommation ; b) on migre avec une vitesse de migration initiale, d'une part, les données amplitude pour obtenir une image migrée de l'amplitude et, d'autre part, ledit produit pour obtenir une image migrée dudit produit, c) on effectue le rapport- des données de l'image migrée dudit produit sur les données de l'image migrée de l'amplitude, pour obtenir une vitesse en position migrée, et d) on procède à une itération des étapes b) et c) en utilisant pour chaque itération une vitesse de migration égale à la vitesse en position migrée issue de l'itération immédiatement précédente jusqu'à ce que l'écart entre la vitesse en position migrée finale d'itération (n) et la vitesse en position migrée d'itération (n-1) soit inférieur à un seuil prédéterminé.
2. Méthode selon la revendication 1 , caractérisée en ce que la migration est une migration en temps.
3. Méthode selon les revendications 1 et 2, caractérisée en ce que l'étape c) est effectuée sur les enveloppes des traces après migration temps de l'amplitude et dudit produit.
4. Méthode selon la revendication 1 , caractérisée en ce que le terme relatif aux vitesses de sommation est affecté d'un exposant inférieur ou égal à l .
5. Méthode selon les revendications 1 à 4, caractérisée en ce que chacun des autres attributs géophysiques est migré en temps avec le champ de vitesses final en effectuant les étapes a) à c) uniquement.
6. Méthode selon la revendication 1 , caractérisé en ce que la vitesse de migration initiale est une vitesse de sommation.
PCT/FR1997/001247 1996-07-19 1997-07-09 Methode de migration des attributs geophysiques d'un milieu WO1998003886A1 (fr)

Priority Applications (5)

Application Number Priority Date Filing Date Title
US09/043,217 US6044039A (en) 1996-07-19 1997-07-09 Geophysical attribute migration method for a medium
AU36252/97A AU712392B2 (en) 1996-07-19 1997-07-09 Geophysical attribute migration method for a medium
EP97932870A EP0888562A1 (fr) 1996-07-19 1997-07-09 Methode de migration des attributs geophysiques d'un milieu
BR9706587A BR9706587A (pt) 1996-07-19 1997-07-09 Método de migração dos atributos geofísicos de um meio
NO980938A NO980938L (no) 1996-07-19 1998-03-04 Fremgangsmåte for migrasjon av geofysisk attributt til et medium

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR96/09094 1996-07-19
FR9609094A FR2751428B1 (fr) 1996-07-19 1996-07-19 Methode de migration des attributs geophysiques d'un milieu

Publications (1)

Publication Number Publication Date
WO1998003886A1 true WO1998003886A1 (fr) 1998-01-29

Family

ID=9494250

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/FR1997/001247 WO1998003886A1 (fr) 1996-07-19 1997-07-09 Methode de migration des attributs geophysiques d'un milieu

Country Status (10)

Country Link
US (1) US6044039A (fr)
EP (1) EP0888562A1 (fr)
CN (1) CN1198216A (fr)
AU (1) AU712392B2 (fr)
BR (1) BR9706587A (fr)
CA (1) CA2229580A1 (fr)
FR (1) FR2751428B1 (fr)
NO (1) NO980938L (fr)
OA (1) OA10673A (fr)
WO (1) WO1998003886A1 (fr)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6466873B2 (en) * 2000-12-18 2002-10-15 Pgs Americas, Inc. Method of extended recursive f-k migration
CN101598806B (zh) * 2008-06-04 2011-05-25 中国石油天然气集团公司 一种提高构造图准确度的偏差消除方法
US8332156B2 (en) * 2009-07-10 2012-12-11 Chevron U.S.A. Inc. Method for propagating pseudo acoustic quasi-P waves in anisotropic media
CN113866827B (zh) * 2021-09-29 2023-05-09 中国石油大学(华东) 一种解释性速度建模地震成像方法、系统、介质和设备

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4813027A (en) * 1987-07-17 1989-03-14 Arabian American Oil Company Method and apparatus for enhancing seismic data

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4888742A (en) * 1987-11-03 1989-12-19 Western Atlas International, Inc. Method of migrating seismic data
US5530679A (en) * 1993-05-10 1996-06-25 Western Atlas International, Inc. Method for migrating seismic data

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4813027A (en) * 1987-07-17 1989-03-14 Arabian American Oil Company Method and apparatus for enhancing seismic data

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
AL-YAHYA: "Velocity analysis by iterative profile migration", GEOPHYSICS, vol. 54, no. 6, June 1989 (1989-06-01), TULSA, pages 718 - 729, XP002030424 *
BAZELAIRE DE E: "NORMAL MOVEOUT REVISITED: INHOMOGENEOUS MEDIA AND CURVED INTERFACES", GEOPHYSICS, vol. 53, no. 2, February 1988 (1988-02-01), pages 143 - 157, XP000617851 *

Also Published As

Publication number Publication date
FR2751428B1 (fr) 1998-08-28
NO980938D0 (no) 1998-03-04
AU3625297A (en) 1998-02-10
BR9706587A (pt) 1999-07-20
MX9802077A (es) 1998-08-30
CA2229580A1 (fr) 1998-01-29
CN1198216A (zh) 1998-11-04
FR2751428A1 (fr) 1998-01-23
OA10673A (en) 2002-09-25
NO980938L (no) 1998-03-04
US6044039A (en) 2000-03-28
AU712392B2 (en) 1999-11-04
EP0888562A1 (fr) 1999-01-07

Similar Documents

Publication Publication Date Title
EP1624321B1 (fr) Méthode pour construire un modèle d'un milieu hétérogène décrit par plusieurs paramètres à partir de données exprimées dans des échelles de temps différentes
FR2708350A1 (fr) Méthode d'analyse de signaux sismiques.
EP0734536B1 (fr) Methode d'analyse et de traitement de donnees sismiques reflexion
CN110780351B (zh) 纵波和转换波叠前联合反演方法及系统
EP0797780B1 (fr) Methode de traitement de traces sismiques reflexion enregistrees pour des deports variables
FR2616920A1 (fr) Inversion d'un profil sismique vertical en minimisant une fonction du type entropie
Zhang et al. Automatic time-domain velocity estimation based on an accelerated clustering method
FR2735244A1 (fr) Methode de traitement pour obtenir des donnees sismiques a deport nul par sommation dans le domaine profondeur
CA2464799C (fr) Methode pour determiner un modele de vitesse d'ondes sismiques dans une formation souterraine heterogene
EP0698220B1 (fr) Procede de modelisation de donnees sismiques cinematiques ayant subi un traitement par au moins un operateur de deplacement
FR2971859A1 (fr) Procede et dispositif de lissage a preservation du temps de trajet
EP0888562A1 (fr) Methode de migration des attributs geophysiques d'un milieu
EP0572321B1 (fr) Procédé de traitement de signaux pour prospection géophysique exploitant un opérateur d'extrapolation d'un champ d'onde perfectionné
FR2755243A1 (fr) Methode perfectionnee de migration avant somme
FR2884619A1 (fr) Procede de traitement de donnees sismiques en vue de la realisation d'un pointe rmo
FR2858064A1 (fr) Procede de pointe bispectral des parametres de correction d'obliquite anelliptique
CN114740528A (zh) 一种超微分拉普拉斯块约束的叠前多波联合反演方法
EP3140677B1 (fr) Procédé de traitement d'images sismiques
FR2751757A1 (fr) Methode d'acquisition et de traitement de donnees sismiques reflexion pour l'exploration d'un milieu a tectonique complexe
FR2801387A1 (fr) Perfectionnements aux procedes de traitement sismique mettant en oeuvre une determination des corrections dynamiques
EP0909961B1 (fr) Procédé de séparation d'ondes en sismique de puits pour des acquisitions de type à déport croissant
CA2085617A1 (fr) Procede de traitement pour l'obtention d'une section somme a offsets nuls
WO2003003052A1 (fr) Perfectionnements aux techniques de pointe de surface
Benamou et al. Numerical microlocal analysis revisited
Collino—SimonMarmorat Numerical MicroLocal Analysis Revisited

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 97190928.8

Country of ref document: CN

AK Designated states

Kind code of ref document: A1

Designated state(s): AU BR CA CN MX NO RU US

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN ML MR NE SN TD TG

WWE Wipo information: entry into national phase

Ref document number: 1997932870

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 2229580

Country of ref document: CA

Ref document number: 2229580

Country of ref document: CA

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: PA/a/1998/002077

Country of ref document: MX

Ref document number: 09043217

Country of ref document: US

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWP Wipo information: published in national office

Ref document number: 1997932870

Country of ref document: EP

WWW Wipo information: withdrawn in national office

Ref document number: 1997932870

Country of ref document: EP