FR2755243A1 - Methode perfectionnee de migration avant somme - Google Patents
Methode perfectionnee de migration avant somme Download PDFInfo
- Publication number
- FR2755243A1 FR2755243A1 FR9613253A FR9613253A FR2755243A1 FR 2755243 A1 FR2755243 A1 FR 2755243A1 FR 9613253 A FR9613253 A FR 9613253A FR 9613253 A FR9613253 A FR 9613253A FR 2755243 A1 FR2755243 A1 FR 2755243A1
- Authority
- FR
- France
- Prior art keywords
- point
- traces
- specular
- dip
- migration
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000005012 migration Effects 0.000 title claims abstract description 24
- 238000013508 migration Methods 0.000 title claims abstract description 24
- 238000012937 correction Methods 0.000 claims abstract description 5
- 238000002310 reflectometry Methods 0.000 claims description 12
- 238000004613 tight binding model Methods 0.000 claims 1
- 230000008901 benefit Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 6
- 238000010304 firing Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241000238558 Eucarida Species 0.000 description 1
- 230000005526 G1 to G0 transition Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- 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
-
- 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/51—Migration
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)
- Analysing Materials By The Use Of Radiation (AREA)
- Holo Graphy (AREA)
Abstract
Elle est du type consistant à enregistrer des traces sismiques réflexion pour constituer une collection de traces, à migrer avant sommation et sans correction de poids les traces classées suivant un critère donné et en utilisant un modèle de vitesses, caractérisée en ce qu'elle consiste en outre à déterminer pour chaque point M d'intérêt dans un domaine (x, z) au moins un pendage déterminé à partir d'une carte de pendage, à déterminer les rayons spéculaires associés au point M et pour ledit pendage, à calculer les caractéristiques temps et amplitude de chacun des rayons spéculaires pour en déduire la réflectivité en chaque point M.
Description
METHODE PERHCTIONNEE DE MIGRATION AVANT SOMME
Un objectif important de la sismique réflexion est de réaliser une section sismique qui corresponde le mieux possible à une image de la réflectivité élastique du sous-sol qui est explore.
Un objectif important de la sismique réflexion est de réaliser une section sismique qui corresponde le mieux possible à une image de la réflectivité élastique du sous-sol qui est explore.
La réflectivité élastique d'un sous-sol peut être utilisée soit pour repositionner dans ltespace la géométrie des horizons ou réflecteurs sismiques, ctest-à-dire des événements qui ont une continuité spatiale et qui sont caractéristiques d'événements géologiques marquants, et alors on obtient des informations sur la structure du sous-sol exploré, soit pour mesurer quantitativement la réflectivité ou le coefficient de réflexion et on obtient alors des informations, d'une part, sur les contrastes de paramètres pétro-élastiques d'un événement géologique donné et, d'autre part, sur la qualité des roches ayant provoqué la réflectivité mesurée.
En application de diverses théories, telles que les théories du point diffractant, il est devenu possible de réaliser des images sismiques fournissant une meilleure représentativité de la réflectivité du sous-sol exploré.
Une des méthodes les plus répandues consiste notamment à - enregistrer des traces sismiques issues du sous-sol - classer ces traces suivant un critère déterminé, par exemple point de tir commun, déport commun, point milieu commun, etc... pour obtenir une ou plusieurs collections de traces classée(s) suivant le critère choisi - migrer en temps ou profondeurJa ou les collections de traces classées préalablement, à l'aide d'un modèle de vitesses préalablement ou non défini mais concernant une zone dudit sous-sol, la migration étant par exemple du type de KIRCHHOFF.
Toutes les méthodes de migration sont basées implicitement ou explicitement sur des sommations (intégrales) et sont bien connus des spécialistes. Les méthodes basées implicitement sur des sommations sont par exemple les méthodes dites par équation d'ondes, alors que les méthodes faisant appel explicitement à des sommations sont par exemple les méthodes de migrations dites de KIRCHHOFF, d'inversion de BORN ou d'autres migrations équivalentes à ces dernières. Chaque trace de la collection migrée peut être ou non pondérée suivant le but que l'on désire atteindre.
Lorsque le but est d'obtenir, rapidement et à faible coût de traitement, une image du sous-sol suffisante pour une interprétation dite structurale, on peut soit ne pas utiliser de poids (traces non pondérées) soit utiliser des poids qui sont grossièrement approximés. Toutefois, une telle technique ne peut être utilisée dans une analyse beaucoup plus fine du type amplitude fonction du déport (désignée par l'acronyme AVO : Amplitude
Versus Offset).
Versus Offset).
Lorsqu'on souhaite par exemple effectuer une migration en amplitude préservée, il est alors indispensable de pondérer les traces par des poids appropriés. On peut utilement se référer à l'article de TYGEL, paru dans Geophysics, vol. 59, No 12 de décembre 1993, qui constitue une bonne synthèse des techniques mises en oeuvre et dans lequel est décrite la théorie de la phase stationnaire.
La pondération propre à chaque trace est appelée "fonction de
GREEN" qui comprend essentiellement deux termes
(a) des temps de trajets dont on a besoin quelle que soit la méthode de migration utilisée, et
(b) des poids qui sont calculés avec précision et appliqués à chaque échantillon utilisé dans la sommation de KIRCHHOFF par exemple.
GREEN" qui comprend essentiellement deux termes
(a) des temps de trajets dont on a besoin quelle que soit la méthode de migration utilisée, et
(b) des poids qui sont calculés avec précision et appliqués à chaque échantillon utilisé dans la sommation de KIRCHHOFF par exemple.
Cela nécessite donc de calculer un très grand nombre de poids, autant qu'il y a d'échantillons intervenant dans la migration, ce qui est relativement très coûteux. Un traitement correct, c'est-à-dire le calcul des poids, immobilise un ordinateur pendant plusieurs jours.
La présente invention a pour but de proposer une méthode d'imagerie ou de représentation d'un sous-sol à explorer qui est suffisamment efficace pour être exploité tout en réduisant les coûts de traitement liés au calcul des poids de la fonction de GREEN.
Un objet de la présente invention est une méthode du type consistant à
- enregistrer des traces sismiques réflexion à partir de couples émetteurs-récepteurs disposés à la surface d'un sous-sol à explorer et pour chacun desquels l'émetteur est séparé du récepteur d'une distance appelée déport,
- produire à partir des traces enregistrées au moins une collection de traces concernant une zone déterminée dudit sous-sol, lesdites traces ayant été classées suivant un critère donné,
- utiliser un modèle de vitesses de propagation concernant ladite zone dans un domaine (x, z),
- calculer dans une première phase uniquement les temps de trajets correspondant auxdites traces,
- migrer avant sommation et sans correction de poids lesdites traces classées en utilisant ledit modèle de vitesses,
- utiliser une carte de pendage préalablement établie pour ladite zone et dans le même domaine (x, z) que le domaine de vitesses, caractérisé en ce qu'elle consiste en outre à :
- déterminer pour chaque point M d'intérêt sélectionné dans ledit domaine (x, z) au moins un pendage,
- déterminer les rayons spéculaires associés au point M pour ledit pendage et pour une gamme de déports préalablement définie,
- calculer les caractéristiques temps et amplitudes de chacun des rayons spéculaires, et à
- déterminer la réflectivité en chaque point M d'intérêt à partir desdites caractéristiques.
- enregistrer des traces sismiques réflexion à partir de couples émetteurs-récepteurs disposés à la surface d'un sous-sol à explorer et pour chacun desquels l'émetteur est séparé du récepteur d'une distance appelée déport,
- produire à partir des traces enregistrées au moins une collection de traces concernant une zone déterminée dudit sous-sol, lesdites traces ayant été classées suivant un critère donné,
- utiliser un modèle de vitesses de propagation concernant ladite zone dans un domaine (x, z),
- calculer dans une première phase uniquement les temps de trajets correspondant auxdites traces,
- migrer avant sommation et sans correction de poids lesdites traces classées en utilisant ledit modèle de vitesses,
- utiliser une carte de pendage préalablement établie pour ladite zone et dans le même domaine (x, z) que le domaine de vitesses, caractérisé en ce qu'elle consiste en outre à :
- déterminer pour chaque point M d'intérêt sélectionné dans ledit domaine (x, z) au moins un pendage,
- déterminer les rayons spéculaires associés au point M pour ledit pendage et pour une gamme de déports préalablement définie,
- calculer les caractéristiques temps et amplitudes de chacun des rayons spéculaires, et à
- déterminer la réflectivité en chaque point M d'intérêt à partir desdites caractéristiques.
Un avantage de la présente invention est de limiter le calcul des poids à un nombre réduit de traces ou plus exactement de rayons spéculaires, comparé au nombre de traces qui est nécessaire dans tous les procédés antérieurs.
Un autre avantage de la présente invention réside dans le fait qu'elle peut être appliquée à des analyses très fines du type AVO.
Un autre avantage de la présente invention est de pouvoir scinder les diverses étapes et de les réaliser séparément dans des entités différentes.
D'autres avantages et caractéristiques apparaîtront plus clairement à la lecture de la description du procédé selon l'invention, ainsi que des dessins annexés sur lesquels
- la figure 1 est une représentation schématique d'un dispositif utilisé en exploration du type sismique réflexion
- la figure 2 est une représentation schématique d'une collection de traces classées à déport constant;
- la figure 3 est une représentation schématique d'une migration profondeur des traces de la collection représentée sur la figure 2;
- la figure 4 est une représentation schématique d'un regroupement de traces constituant une collection iso-XM;
- la figure 5 est une représentation schématique d'un tracé de rayons spéculaires dans le modèle de vitesses de propagation.
- la figure 1 est une représentation schématique d'un dispositif utilisé en exploration du type sismique réflexion
- la figure 2 est une représentation schématique d'une collection de traces classées à déport constant;
- la figure 3 est une représentation schématique d'une migration profondeur des traces de la collection représentée sur la figure 2;
- la figure 4 est une représentation schématique d'un regroupement de traces constituant une collection iso-XM;
- la figure 5 est une représentation schématique d'un tracé de rayons spéculaires dans le modèle de vitesses de propagation.
A la surface 1 d'un sous-sol 2 dont on souhaite réaliser une image de la réflectivité, on dispose un ou plusieurs émetteurs S et des récepteurs
R1 à Rns la distance hi séparant l'émetteur Si d'un récepteur Ri étant appelée déport (offset en anglais). A partir de l'émetteur ou source Si, on émet des ondes qui se propagent dans le sous-sol 2 et qui, après réflexion sur des réflecteurs ou horizons H, horizontaux ou parallèles, aboutissent aux récepteurs R1 à Rn dans lesquels elles sont enregistrées sous forme de traces. Les traces sont, dans une première étape, classées suivant un critère déterminé, par exemple en point de tir commun, en point milieu commun, en récepteur commun, à déport commun, etc... Un exemple de collection de traces est représenté très schématiquement sur la figure 2 et concerne une zone déterminée du sous-sol 2, le classement ayant été effectué à déport constant hi. Pour un autre déport constant hi+l on aurait une autre collection de traces analogue à celle de la figure 2 ; de même, pour un autre déport hi+2,etc...
R1 à Rns la distance hi séparant l'émetteur Si d'un récepteur Ri étant appelée déport (offset en anglais). A partir de l'émetteur ou source Si, on émet des ondes qui se propagent dans le sous-sol 2 et qui, après réflexion sur des réflecteurs ou horizons H, horizontaux ou parallèles, aboutissent aux récepteurs R1 à Rn dans lesquels elles sont enregistrées sous forme de traces. Les traces sont, dans une première étape, classées suivant un critère déterminé, par exemple en point de tir commun, en point milieu commun, en récepteur commun, à déport commun, etc... Un exemple de collection de traces est représenté très schématiquement sur la figure 2 et concerne une zone déterminée du sous-sol 2, le classement ayant été effectué à déport constant hi. Pour un autre déport constant hi+l on aurait une autre collection de traces analogue à celle de la figure 2 ; de même, pour un autre déport hi+2,etc...
Dans une deuxième étape, on effectue une migration temps ou profondeur de chaque collection de traces comme celle de la figure 2, par exemple une migration profondeur en utilisant notamment la migration avant somme de KIRCHHOFF mais dans une version simplifiée et rapide, afin d'obtenir une série d'images migrées intermédiaires, ce qui correspond à une image par déport ou par point de tir ou par tout autre paramètre X commun.. Pour ce faire, on utilise, de manière classique, un modèle de vitesses de propagation en x,z qui est établi ou, de préférence, qui a été déjà établi pour ladite zone concernée. Les images migrées peuvent être produites en 2-D (deux dimensions) ou en 3-D (trois dimensions) dans le domaine temps ou profondeur, comme celle représentée très schématiquement en 2-D sur la figure 3.
Un groupement des traces de chacune des migrations intermédiaires correspondant à une même position X en surface d'un point M est appelé collection iso-valeur (iso-XM) ou "image gather".
Après migration des collections de traces correspondant à chaque déport constant, comme celles représentées sur la figure 3, on obtient une "image gather" pour les divers déports considérés hi, hi+1, hi+2, ...
telle que représenté sur la figure 4.
A partir d'une quelconque des images migrées intermédiaires obtenues précédemment, de leur somme ou encore à partir de toute autre image préexistante comme l'image migration profondeur post-sommation, on définit, par des moyens connus tels que ceux développés par la demanderesse et connus sous la dénomination SISMAGE, une carte des pendages dans le domaine (x, z) de la zone considérée.
Dans le modèle de vitesses, représenté très schématiquement sur la figure 5, on choisit un point quelconque M de coordonnées XM, zM.
A l'aide de la carte de pendages précédente, on détermine la valeur et la direction du pendage audit point M. Sur la figure 5, le pendage en M est représenté symboliquement par un trait 4.
Dans une autre étape, on calcule un premier rayon spéculaire au point M, par exemple en se fixant un point source S;, un temps de trajet t le long du rayon SURI, Ri étant le point de mesure à la surface, ainsi que les paramètres relatifs à l'évolution des amplitudes des ondes le long dudit trajet ("dynamic ray-tracing"). On répète cette étape pour différents couples S;R;, ce qui permet d'obtenir un ensemble de rayons spéculaires au point M pour une gamme de déports donnée.
D'autres méthodes existent pour trouver le couple SiRi qui porte la réflexion spéculaire correspondante au point M ; par exemple si on "tire" de Si vers M, puis de M vers la surface en respectant la loi de réflexion en
M, le rayon issu de M n'a aucune raison d'atteindre la surface en un point de mesure Ri. On est donc amené à réaliser des interpolations jusqu'à ce qu'on trouve le couple SjRj qui, associé au point M, est spéculaire.
M, le rayon issu de M n'a aucune raison d'atteindre la surface en un point de mesure Ri. On est donc amené à réaliser des interpolations jusqu'à ce qu'on trouve le couple SjRj qui, associé au point M, est spéculaire.
Chaque rayon spéculaire au point M fournit une relation entre le déport hi = SjRj et l'angle de réflexion 0; associés audit rayon spéculaire.
Pour chaque rayon spéculaire, on calcule la fonction de GREEN qui lui est associée, ce qui permet de déterminer le poids qu'il faut appliquer à l'échantillon de chaque migration intermédiaire spécifique au paramètre constant, par exemple hi, pour avoir une bonne représentation de la réflectivité au point M pour ce paramètre hi.
Ainsi, pour un déport donné hi, on a un seul rayon spéculaire et un seul poids de correction correspondant à ce rayon spéculaire.
Dans chaque image de migration intermédiaire, on a pour un même point M de coordonnées XM, ZM, un ensemble d'échantillons de la collection iso-X ou "image-gather", chaque échantillon correspondant à un déport donné ou un point de tir donné de la gamme choisie pour le calcul des rayons spéculaires. Comme à chaque échantillon de chaque trace de la collection iso-X correspond un rayon spéculaire du modèle de vitesses, il devient alors possible d'appliquer audit échantillon le poids de la fonction de GREEN associé au rayon spéculaire.
Le résultat de cette opération est que l'échantillon à la profondeur
ZM est proportionnel à la réflectivité du point M pour l'ange O correspondant.
ZM est proportionnel à la réflectivité du point M pour l'ange O correspondant.
Ainsi, on peut produire une série d'images migrées finales, de même nature que celles que l'on aurait obtenues avec une migration complète du type de KIRCHHOFF, et qui, avec l'information sur les angles de réflexion 0, seront utilisables dans le processus bien connu de l'analyse AVO puisqu'on peut établir des cartes R = f(O).
Bien entendu, les diverses étapes et opérations décrites précédemment sont à nouveau effectuées pour tous les points M de la zone d'intérêt du domaine (x, y, z) considéré.
Comme on peut le constater de ce qui précède, la présente invention constitue un perfectionnement important aux techniques habituellement utilisées puisque, d'une part, on évite le calcul d'un angle de réflexion par une triple sommation de KIRCHHOFF en lui substituant une mesure directe sur une image pré-établie et un lancer de rayons plus faciles à effectuer et beaucoup moins coûteux, et d'autre part, on ne calcule les poids correctifs à appliquer pour la mesure de la réflectivité que pour les rayons spéculaires spécifiques audit angle de réflexion et non pas pour tous angles de réflexion correspondant à toutes les collections de traces d'origine et constituées par tous les enregistrements des récepteurs.
Un autre avantage, non négligeable, est qu'on peut, si on le souhaite, scinder les diverses étapes du procédé selon l'invention en plusieurs parties, chacune desdites parties pouvant être effectuée dans une entité distincte. Par exemple, le calcul des images simplifiées intermédiaires sans calcul des poids des fonctions de GREEN peut être confié à une entité spécialisée, comme un contracteur.
La description ci-dessus et les figures annexées correspondent à une mise en oeuvre pour des collections iso-offset. Toutefois, il est parfaitement possible, selon l'invention, de sélectionner d'autres paramètres communs tels que le point de tir, la méthode selon l'invention, étant alors appliquée à des collections iso-point de tir.
Un avantage important de la présente invention tel que cela apparaît à la lumière de la description, réside dans le fait qu'au lieu de calculer et d'appliquer tous les poids correspondants à tous les pendages possibles en
M, ce qui est coûteux, on ne calcule que le poids lié au pendage mesuré en M. Il est également possible, pour affiner les résultats, d'appliquer la méthode selon l'invention à une petite gamme des valeurs de pendage au point M. De plus, le poids de la fonction de Green associé au rayon spéculaire peut être un scalaire ou un nombre complexe.
M, ce qui est coûteux, on ne calcule que le poids lié au pendage mesuré en M. Il est également possible, pour affiner les résultats, d'appliquer la méthode selon l'invention à une petite gamme des valeurs de pendage au point M. De plus, le poids de la fonction de Green associé au rayon spéculaire peut être un scalaire ou un nombre complexe.
I1 est également possible de faire varier le pendage au point M autour de la valeur nominale telle qu'elle fut calculée ou déterminée, par exemple, par le procédé SISMAGE. Dans ce cas, on détermine plusieurs valeurs du pendage, on détermine les rayons spéculaires associés au point
M pour chacune desdites valeurs du pendage et on effectue une sélection, une moyenne ou toute autre méthode statistique pour déterminer le poids de la fonction de Green et le pendage affectés au point M. Ceci est fait afin de minimiser l'imprécision de la valeur nominale.
M pour chacune desdites valeurs du pendage et on effectue une sélection, une moyenne ou toute autre méthode statistique pour déterminer le poids de la fonction de Green et le pendage affectés au point M. Ceci est fait afin de minimiser l'imprécision de la valeur nominale.
De même, il est possible de faire varier très légèrement la position du point M, de l'émetteur Si et/ou du récepteur Ri qui sont associés audit point M, dans la limite imposée par la théorie de Fresnel. Une telle variation est dite réalisée dans la zone de Fresnel et elle est bien connue des spécialistes. Selon qu'on a choisi la première ou la deuxième zone de
Fresnel, la variation d'une ou de plusieurs positions doit être telle que le temps de trajet Si M est inférieur au quart ou à la demi longueur d'onde.
Fresnel, la variation d'une ou de plusieurs positions doit être telle que le temps de trajet Si M est inférieur au quart ou à la demi longueur d'onde.
Claims (8)
1. Méthode perfectionnée de migration avant somme du type consistant à
- enregistrer des traces sismiques réflexion à partir de couples émetteurs-récepteurs disposés à la surface d'un sous-sol à explorer et pour chacun desquels l'émetteur est séparé du récepteur d'une distance appelée déport,
- produire à partir des traces enregistrées au moins une collection de traces concernant une zone déterminée dudit sous-sol, lesdites traces ayant été classées suivant un critère donné,
- utiliser un modèle de vitesses de propagation concernant ladite zone dans un domaine (x, z),
- calculer dans une première phase uniquement les temps de trajets correspondant auxdites traces,
- migrer avant sommation et sans correction de poids lesdites traces classées en utilisant ledit modèle de vitesses,
- utiliser une carte de pendage préalablement établie pour ladite zone et dans le même domaine (x, z) que le domaine de vitesses, caractérisée en ce qu'elle consiste en outre a :
- déterminer pour chaque point M d'intérêt sélectionné dans ledit domaine (x, z) au moins un pendage relevé sur ladite carte,
- déterminer les rayons spéculaires associés au point M pour ledit pendage et pour une gamme de déports préalablement définie,
- calculer les caractéristiques temps et amplitude de chacun des rayons spéculaires, et à
- déterminer la réflectivité en chaque point M d'intérêt à partir desdites caractéristiques.
2. Méthode selon la revendication 1, caractérisée en ce que les traces migrées sans correction de poids qui correspondent à une même position en surface du point M sont regroupées en collection iso-valeur.
3. Méthode selon la revendication 2, caractérisée en ce qutà chacun des échantillons de la collection de traces iso-valeur qui correspond au point M, on applique le poids de la fonction de Green du rayon spéculaire associé à l'échantillon.
4. Méthode selon les revendications 1 ou 2, caractérisée en ce que la migration avant somme est une migration profondeur.
5. Méthode selon les revendications 1 ou 2, caractérisée en ce que la migration avant somme est une migration temps.
6. Méthode selon la revendication 1, caractérisée en ce qu'on détermine pour chaque rayon spéculaire la relation entre le déport et l'angle de réflexion puis la relation entre la réflectivité et l'angle de réflexion et on utilise lesdites relations dans la technique dite Amplitude
Versus Offset (AVO).
7. Méthode selon la revendication 1, caractérisée en ce qu'on détermine plusieurs valeurs du pendage pour chaque point M d'intérêt et on détermine les rayons spéculaires associés au point M pour chacune desdites valeurs du pendage, puis on effectue la moyenne desdits rayons spéculaires.
8. Méthode selon la revendication 1, caractérisée en ce qu'elle consiste à faire varier au moins une des positions du point M d'intérêt, de l'émetteur Si et/ou du récepteur Ri qui sont associés audit point M, à déterminer les rayons spéculaires correspondant auxdites variations, puis à effectuer la moyenne des rayons spéculaires.
Priority Applications (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9613253A FR2755243B1 (fr) | 1996-10-30 | 1996-10-30 | Methode perfectionnee de migration avant somme |
US09/091,995 US6094621A (en) | 1996-10-30 | 1997-10-29 | Method for migration before summation |
CA002240998A CA2240998A1 (fr) | 1996-10-30 | 1997-10-29 | Methode perfectionnee de migration avant somme |
PCT/FR1997/001944 WO1998019180A1 (fr) | 1996-10-30 | 1997-10-29 | Methode perfectionnee de migration avant somme |
EP97912284A EP0879428A1 (fr) | 1996-10-30 | 1997-10-29 | Methode perfectionnee de migration avant somme |
NO982834A NO982834L (no) | 1996-10-30 | 1998-06-19 | Fremgangsmåte for migrasjon för summering |
OA9800095A OA10802A (en) | 1996-10-30 | 1998-06-30 | Improved method for migration before summation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9613253A FR2755243B1 (fr) | 1996-10-30 | 1996-10-30 | Methode perfectionnee de migration avant somme |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2755243A1 true FR2755243A1 (fr) | 1998-04-30 |
FR2755243B1 FR2755243B1 (fr) | 1998-12-04 |
Family
ID=9497178
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR9613253A Expired - Fee Related FR2755243B1 (fr) | 1996-10-30 | 1996-10-30 | Methode perfectionnee de migration avant somme |
Country Status (7)
Country | Link |
---|---|
US (1) | US6094621A (fr) |
EP (1) | EP0879428A1 (fr) |
CA (1) | CA2240998A1 (fr) |
FR (1) | FR2755243B1 (fr) |
NO (1) | NO982834L (fr) |
OA (1) | OA10802A (fr) |
WO (1) | WO1998019180A1 (fr) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001075482A2 (fr) * | 2000-03-30 | 2001-10-11 | Nutec Sciences, Inc. | Procede de traitement de donnees sismiques |
FR2878966A1 (fr) * | 2004-12-07 | 2006-06-09 | Inst Francais Du Petrole | Methode pour determiner des informations speculaires apres imagerie sismique avant sommation |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO320271B1 (no) * | 2000-05-26 | 2005-11-14 | Norsar V Daglig Leder | Fremgangsmate og anordning for a finne refleksjonsstyrken til reflektorer i undergrunnen mens refleksjonskoeffisienten for en malhorisont holdes konstant |
US6920084B2 (en) * | 2003-01-15 | 2005-07-19 | Western Geco, L.L.C. | Offset transformation to zero dip that preserves angle of incidence |
NO322089B1 (no) | 2003-04-09 | 2006-08-14 | Norsar V Daglig Leder | Fremgangsmate for simulering av lokale prestakk dypmigrerte seismiske bilder |
US6925387B2 (en) * | 2003-08-14 | 2005-08-02 | Westerngeco L.L.C. | Method and apparatus for kinematically linking multiple seismic domains |
FR2886740B1 (fr) * | 2005-06-03 | 2007-09-28 | Inst Francais Du Petrole | Methode pour mettre a jour un modele geologique par des donnees sismiques et de production |
US20140129188A1 (en) * | 2012-11-07 | 2014-05-08 | Chevron U.S.A. Inc. | System and method for analysis of seismic images |
EP2784553A3 (fr) * | 2013-03-26 | 2015-10-21 | CGG Services SA | Prédiction de la sensibilité à un positionnement pour des relevés sismiques |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2163851A (en) * | 1984-08-27 | 1986-03-05 | Amoco Corp | Method of seismic exploration |
US4992996A (en) * | 1990-01-31 | 1991-02-12 | Conoco Inc. | Interval velocity analysis and depth migration using common reflection point gathers |
EP0430688A2 (fr) * | 1989-11-30 | 1991-06-05 | Conoco Inc. | Migration de profondeur à étalement commun avec correction dynamique |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4964103A (en) * | 1989-07-13 | 1990-10-16 | Conoco Inc. | Three dimensional before stack depth migration of two dimensional or three dimensional seismic data |
US5128899A (en) * | 1991-07-16 | 1992-07-07 | Conoco Inc. | Method for sharpening prestack depth migration images by removing residual moveout from common image point gathers before stacking |
US5285422A (en) * | 1992-10-23 | 1994-02-08 | Western Atlas International, Inc. | Method for compensating 3D DMO for the effects of an inhomogeneous earth |
US5530679A (en) * | 1993-05-10 | 1996-06-25 | Western Atlas International, Inc. | Method for migrating seismic data |
US5587942A (en) * | 1994-03-02 | 1996-12-24 | Exxon Production Research Company | 3D wave equation migration of a 2D grid of seismic data |
US5570321A (en) * | 1994-03-03 | 1996-10-29 | Atlantic Richfield Company | Seismic velocity model optimization method using simulated annearling to determine prestack travel-times |
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. |
US5629904A (en) * | 1994-11-30 | 1997-05-13 | Paradigm Geophysical, Ltd. | Migration process using a model based aperture technique |
US5596547A (en) * | 1995-10-20 | 1997-01-21 | Bancroft; John C. | Prestack seismic migration |
-
1996
- 1996-10-30 FR FR9613253A patent/FR2755243B1/fr not_active Expired - Fee Related
-
1997
- 1997-10-29 US US09/091,995 patent/US6094621A/en not_active Expired - Fee Related
- 1997-10-29 WO PCT/FR1997/001944 patent/WO1998019180A1/fr not_active Application Discontinuation
- 1997-10-29 CA CA002240998A patent/CA2240998A1/fr not_active Abandoned
- 1997-10-29 EP EP97912284A patent/EP0879428A1/fr not_active Withdrawn
-
1998
- 1998-06-19 NO NO982834A patent/NO982834L/no not_active Application Discontinuation
- 1998-06-30 OA OA9800095A patent/OA10802A/en unknown
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2163851A (en) * | 1984-08-27 | 1986-03-05 | Amoco Corp | Method of seismic exploration |
EP0430688A2 (fr) * | 1989-11-30 | 1991-06-05 | Conoco Inc. | Migration de profondeur à étalement commun avec correction dynamique |
US4992996A (en) * | 1990-01-31 | 1991-02-12 | Conoco Inc. | Interval velocity analysis and depth migration using common reflection point gathers |
Non-Patent Citations (3)
Title |
---|
LANDA ET AL.: "Model-based stack: a method for constructing an accurate zero-offset section for complex overburdens", GEOPHYS. PROSPECT., vol. 41, no. 6, August 1993 (1993-08-01), THE NETHERLANDS, pages 661 - 670, XP002037664 * |
TYGEL ET AL.: "Pulse distortion in depth migration", GEOPHYSICS, vol. 59, no. 10, October 1994 (1994-10-01), TULSA, OKLAHOMA, USA, pages 1591 - 1569, XP002037665 * |
WOJCIECH DEBSKI ET AL: "INFORMATION ON ELASTIC PARAMETERS OBTAINED FROM THE AMPLITUDES OF REFLECTED WAVES", GEOPHYSICS, vol. 60, no. 5, 1 September 1995 (1995-09-01), pages 1426 - 1436, XP000536271 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001075482A2 (fr) * | 2000-03-30 | 2001-10-11 | Nutec Sciences, Inc. | Procede de traitement de donnees sismiques |
WO2001075482A3 (fr) * | 2000-03-30 | 2002-03-21 | Nutec Sciences Inc | Procede de traitement de donnees sismiques |
FR2878966A1 (fr) * | 2004-12-07 | 2006-06-09 | Inst Francais Du Petrole | Methode pour determiner des informations speculaires apres imagerie sismique avant sommation |
EP1669778A1 (fr) * | 2004-12-07 | 2006-06-14 | Institut Français du Pétrole | Méthode pour déterminer des informations spéculaires après imagerie sismique avant sommation |
US7345951B2 (en) | 2004-12-07 | 2008-03-18 | Institut Francais Du Petrole | Method of determining specular information after prestack seismic imaging |
Also Published As
Publication number | Publication date |
---|---|
US6094621A (en) | 2000-07-25 |
EP0879428A1 (fr) | 1998-11-25 |
CA2240998A1 (fr) | 1998-05-07 |
FR2755243B1 (fr) | 1998-12-04 |
OA10802A (en) | 2003-01-29 |
NO982834D0 (no) | 1998-06-19 |
WO1998019180A1 (fr) | 1998-05-07 |
NO982834L (no) | 1998-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0354112B1 (fr) | Méthode pour obtenir un modèle représentatif d'un milieu hétérogène et notamment du sous-sol | |
CA2374160C (fr) | Perfectionnements aux procedes d'inversion tomographique d'evenements pointes sur les donnees sismiques migrees | |
EP2253970B1 (fr) | Méthode pour imager une zone cible du sous-sol à partir de données de type walkaway | |
EP0249527B1 (fr) | Procédé et dispositif de traitement de données sismographique par corrélation | |
CA2514112A1 (fr) | Methode pour construire un modele d'un milieu heterogene decrit par plusieurs parametres a partir de donnees exprimees dans des echelles de temps differentes | |
EP0734536B1 (fr) | Methode d'analyse et de traitement de donnees sismiques reflexion | |
FR2916859A1 (fr) | Procede de traitement de donnees sismiques | |
EP0797780B1 (fr) | Methode de traitement de traces sismiques reflexion enregistrees pour des deports variables | |
FR2747477A1 (fr) | Methode d'inversion pour donnees sismiques | |
FR2755243A1 (fr) | Methode perfectionnee de migration avant somme | |
FR2616920A1 (fr) | Inversion d'un profil sismique vertical en minimisant une fonction du type entropie | |
FR2792080A1 (fr) | Procedes de traitement de champs d'ondes sismiques, de compensation de l'effet d'une illumination irreguliere de reflecteurs souterrains, d'analyse de caracteres d'illumination, et de filtrage de donnees sismiques | |
CA2195497C (fr) | Procede de prospection sismique marine au moyen d'un couple de capteurs hydrophone et geophone | |
EP0698220B1 (fr) | Procede de modelisation de donnees sismiques cinematiques ayant subi un traitement par au moins un operateur de deplacement | |
CA2182111C (fr) | Methode et dispositif de filtrage d'ondes elliptiques se propageant dans un milieu | |
FR2884620A1 (fr) | Procede de traitement de donnees sismiques en vue d'une caracterisation avo ou avoa | |
FR2858064A1 (fr) | Procede de pointe bispectral des parametres de correction d'obliquite anelliptique | |
EP0909961B1 (fr) | Procédé de séparation d'ondes en sismique de puits pour des acquisitions de type à déport croissant | |
FR2845483A1 (fr) | Procede de reduction du bruit harmonique dans des signaux vibrosismiques. | |
FR2691811A1 (fr) | Procédé de détermination d'un modèle de distribution optimale des vitesses de propagation d'ondes sismiques dans le sous-sol. | |
EP0610123B1 (fr) | Procédé pour améliorer l'estimation, par analyse de focalisation, des vitesses de propagation des ondes sismiques | |
CA2085617A1 (fr) | Procede de traitement pour l'obtention d'une section somme a offsets nuls | |
FR2751428A1 (fr) | Methode de migration des attributs geophysiques d'un milieu | |
FR2751427A1 (fr) | Methode d'attenuation multiple | |
FR2650676A1 (fr) | Procede de traitement de donnees de sismique reflexion pour l'obtention de coupes sismiques ameliorees |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
ST | Notification of lapse |
Effective date: 20070629 |