FR2751428A1 - Methode de migration des attributs geophysiques d'un milieu - Google Patents
Methode de migration des attributs geophysiques d'un milieu Download PDFInfo
- Publication number
- FR2751428A1 FR2751428A1 FR9609094A FR9609094A FR2751428A1 FR 2751428 A1 FR2751428 A1 FR 2751428A1 FR 9609094 A FR9609094 A FR 9609094A FR 9609094 A FR9609094 A FR 9609094A FR 2751428 A1 FR2751428 A1 FR 2751428A1
- Authority
- FR
- France
- Prior art keywords
- migrated
- speed
- migration
- amplitude
- iteration
- Prior art date
- Legal status (The legal status 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 status listed.)
- Granted
Links
- 238000013508 migration Methods 0.000 title claims abstract description 45
- 230000005012 migration Effects 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000012804 iterative process Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 235000020281 long black Nutrition 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
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
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
Méthode de migration des attributs géographiques d'un milieu. 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 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 differents 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, observes 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 media 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.
Inhomogeneous media 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.
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 pentés, 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 pentés, 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-l) soit inférieur à un seuil prédéterminé.
- 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-l) 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 1' amplitude avec des techniques classiques de migration.
- 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 1' 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 H1 à H6 situés à des profondeurs situées entre 490 m pour le réflecteur H1 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 H1 à 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 audessus 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.
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.
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 a qui est au plus égal à 1, afin d'atténuer les effets parasites de la migration.
De préférence, l'exposant a 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*Va.
La deuxième opération consiste à migrer d'une part la trace amplitude A et, d'autre part, le produit A*Va 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éfere 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*Va. Cependant, toute autre migration temps peut être utilisée, comme la migration par équation des ondes, technique que les spécialistes connaissent bien.
A et A*Va. 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*Va) sur l'image migrée A, de sorte qu'on obtient une première vitesse en position migrée V car
(A *Va) migré a
VQ migré
A migré
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*Va.
(A *Va) migré a
VQ migré
A migré
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*Va.
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 V1. Ainsi, on migre en temps, avec la vitesse V1, le terme amplitude A et le produit A*Va, puis on procède au rapport
(A *Va) migré
VQ migré
(A) migré pour obtenir à l'issue de l'itération 1 une vitesse en position migrée V2.
(A *Va) migré
VQ migré
(A) migré 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*Va, 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. I1 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-l (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 H2 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 ATT 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 a de migrer en temps le terme amplitude
TT'
A et le produit A*Aa 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*ATT) migré
A migré pour avoir ATT en position migrée et ce, sans itération puisque le champ de vitesses de migration est maintenant connu.
TT'
A et le produit A*Aa 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*ATT) migré
A migré pour avoir ATT en position migrée et ce, sans itération puisque le champ de vitesses de migration est maintenant connu.
I1 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.
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 (6)
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-l) 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 1' 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 à 1.
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.
Priority Applications (11)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9609094A FR2751428B1 (fr) | 1996-07-19 | 1996-07-19 | Methode de migration des attributs geophysiques d'un milieu |
EP97932870A EP0888562A1 (fr) | 1996-07-19 | 1997-07-09 | Methode de migration des attributs geophysiques d'un milieu |
US09/043,217 US6044039A (en) | 1996-07-19 | 1997-07-09 | Geophysical attribute migration method for a medium |
CA002229580A CA2229580A1 (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 |
AU36252/97A AU712392B2 (en) | 1996-07-19 | 1997-07-09 | Geophysical attribute migration method for a medium |
CN97190928A CN1198216A (zh) | 1996-07-19 | 1997-07-09 | 介质地球物理属性偏移校正方法 |
PCT/FR1997/001247 WO1998003886A1 (fr) | 1996-07-19 | 1997-07-09 | Methode de migration des attributs geophysiques d'un milieu |
NO980938A NO980938L (no) | 1996-07-19 | 1998-03-04 | Fremgangsmåte for migrasjon av geofysisk attributt til et medium |
MXPA/A/1998/002077A MXPA98002077A (en) | 1996-07-19 | 1998-03-17 | Method of displacement of the geophysical attributes of a me |
OA9800031A OA10673A (en) | 1996-07-19 | 1998-03-18 | Geophysical attribute migration method for a medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9609094A FR2751428B1 (fr) | 1996-07-19 | 1996-07-19 | Methode de migration des attributs geophysiques d'un milieu |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2751428A1 true FR2751428A1 (fr) | 1998-01-23 |
FR2751428B1 FR2751428B1 (fr) | 1998-08-28 |
Family
ID=9494250
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR9609094A Expired - Fee Related FR2751428B1 (fr) | 1996-07-19 | 1996-07-19 | 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)
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)
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)
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 |
-
1996
- 1996-07-19 FR FR9609094A patent/FR2751428B1/fr not_active Expired - Fee Related
-
1997
- 1997-07-09 EP EP97932870A patent/EP0888562A1/fr not_active Withdrawn
- 1997-07-09 CN CN97190928A patent/CN1198216A/zh active Pending
- 1997-07-09 AU AU36252/97A patent/AU712392B2/en not_active Ceased
- 1997-07-09 US US09/043,217 patent/US6044039A/en not_active Expired - Fee Related
- 1997-07-09 BR BR9706587A patent/BR9706587A/pt unknown
- 1997-07-09 WO PCT/FR1997/001247 patent/WO1998003886A1/fr not_active Application Discontinuation
- 1997-07-09 CA CA002229580A patent/CA2229580A1/fr not_active Abandoned
-
1998
- 1998-03-04 NO NO980938A patent/NO980938L/no not_active Application Discontinuation
- 1998-03-18 OA OA9800031A patent/OA10673A/en unknown
Patent Citations (1)
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)
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 |
---|---|
AU3625297A (en) | 1998-02-10 |
BR9706587A (pt) | 1999-07-20 |
MX9802077A (es) | 1998-08-30 |
OA10673A (en) | 2002-09-25 |
US6044039A (en) | 2000-03-28 |
NO980938D0 (no) | 1998-03-04 |
CN1198216A (zh) | 1998-11-04 |
FR2751428B1 (fr) | 1998-08-28 |
NO980938L (no) | 1998-03-04 |
WO1998003886A1 (fr) | 1998-01-29 |
EP0888562A1 (fr) | 1999-01-07 |
CA2229580A1 (fr) | 1998-01-29 |
AU712392B2 (en) | 1999-11-04 |
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 | |
EP0354112B1 (fr) | Méthode pour obtenir un modèle représentatif d'un milieu hétérogène et notamment du sous-sol | |
FR2843202A1 (fr) | Methode pour former un modele representatif de la distribution d'une grandeur physique dans une zone souterraine, affranchi de l'effet de bruits correles entachant des donnees d'exploration | |
FR2708350A1 (fr) | Méthode d'analyse de signaux sismiques. | |
FR2765692A1 (fr) | Methode pour modeliser en 3d l'impedance d'un milieu heterogene | |
EP0734536B1 (fr) | Methode d'analyse et de traitement de donnees sismiques reflexion | |
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 | |
FR2747477A1 (fr) | Methode d'inversion pour donnees sismiques | |
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 | |
EP0698220B1 (fr) | Procede de modelisation de donnees sismiques cinematiques ayant subi un traitement par au moins un operateur de deplacement | |
FR2751428A1 (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 | |
FR2858064A1 (fr) | Procede de pointe bispectral des parametres de correction d'obliquite anelliptique | |
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 | |
WO2015159000A2 (fr) | Procédé de traitement d'images sismiques | |
CA2085617A1 (fr) | Procede de traitement pour l'obtention d'une section somme a offsets nuls | |
Garabito et al. | Combination of the common reflection surface-based prestack data regularization and reverse time migration: Application to real land data | |
Benamou et al. | Numerical microlocal analysis revisited | |
FR2769718A1 (fr) | Procede de separation d'ondes en sismique de puits pour des acquisitions de type a deport croissant | |
WO2024149595A1 (fr) | Procédé de surveillance d'un point d'intérêt d'un sous-sol | |
EP0610123A1 (fr) | Procédé pour améliorer l'estimation, par analyse de focalisation, des vitesses de propagation des ondes sismiques |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
ST | Notification of lapse |