CA2284386C - Methode pour realiser en 3d avant sommation, une migration de donnees sismiques - Google Patents
Methode pour realiser en 3d avant sommation, une migration de donnees sismiques Download PDFInfo
- Publication number
- CA2284386C CA2284386C CA002284386A CA2284386A CA2284386C CA 2284386 C CA2284386 C CA 2284386C CA 002284386 A CA002284386 A CA 002284386A CA 2284386 A CA2284386 A CA 2284386A CA 2284386 C CA2284386 C CA 2284386C
- Authority
- CA
- Canada
- Prior art keywords
- seismic
- migration
- wave
- field
- vector
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000005012 migration Effects 0.000 title claims abstract description 39
- 238000013508 migration Methods 0.000 title claims abstract description 39
- 238000003384 imaging method Methods 0.000 claims abstract description 7
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 3
- 239000002131 composite material Substances 0.000 claims description 25
- 230000001419 dependent effect Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000002441 reversible effect Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 8
- 230000005284 excitation Effects 0.000 abstract description 2
- 238000003786 synthesis reaction Methods 0.000 abstract description 2
- 238000005070 sampling Methods 0.000 description 5
- 230000000737 periodic effect Effects 0.000 description 4
- 238000010304 firing Methods 0.000 description 3
- 238000011282 treatment Methods 0.000 description 3
- 230000003111 delayed effect Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
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)
Abstract
Méthode de migration avant sommation (« prestack migration ») permettant d'imager une zone souterraine, pour un modèle de vitesses donné de complexité arbitraire. - Au moyen d'outils classiques de modélisation de la propagation et de la rétropropagation des ondes, la méthode permet d'obtenir des images migrées élémentaires associées aux valeurs que prend un paramètre et la somme des images obtenues pour les différentes valeurs du paramètre (sommation après migratio n) et ceci aussi bien dans le domaine profondeur que dans le domaine temps. Cette migration est obtenue à un prix (coût de calcul) avantageux car indépendant du volume de résultats calculé et du nombre de traces sismiques enregistrées. Seul le volume de la zone dans lequel les ondes se propagent, la complexité des événements à imager et la précision souhaitée ont une incidence sur le coût du calcul. On obtient ainsi des imag es volumiques en prenant en compte l'ensemble des traces sismiques. On a ainsi la possibilité de mettre en oeuvre une procédure de migration par ondes planes dans le cas où l'acquisition ne permet pas la synthèse de la réponse de la subsurface à une excitation onde plane, réponse requise en entrée dans les algorithmes classiques de migration par ondes planes. - Applications à l'imagerie d'interfaces géologiques ou hétérogénéités d'une partie de zone souterraine.
Description
METHODE POUR RÉALISER EN 3D AVANT SOMMATION, UNE MIGRATION DE DONNÉES SISMIQUES
La présente invention concerne une méthode pour réaliser une migration avant sommation d'événements sismiques enregistrés pour l'imagerie d'une partie d'une zone souterraine.
La méthode selon l'invention permet par exemple de réaliser une migration profondeur de type 3D avant sommation ( 3D prestack depth migration ), pour un modèle de vitesses donné, de manière à imager les interfaces géologiques ou hétérogénéités diverses d'une partie de la subsurface.
Art antérieur La migration avant sommation est une méthode classique de traitement de données sismiques. De façon générale la technique consiste, connaissant la valeur d'un champ d'ondes à une profondeur connue, par exemple en surface, ainsi qu'un modèle.
de la répartition des vitesses de propagation des ondes dans la zone, à modéliser la propagation du champ source et la retropropagation des données de réflexion enregistrées et à chercher des cohérences de phase entre ces deux champs modélisés.
On distingue trois grands types de migration avant sommation :
- la migration par point de tir : le champ source est l'état vibratoire engendré par le point de tir et les données de réflexion sont la réponse de la subsurface à ce champ source ;
- la migration par ondes planes aussi appelée migration par angle d'éclairage commun : le champ source est l'onde plane considérée et les données de réflexion sont la réponse de la subsurface à ce champ source ;
- la migration par déport ( offset ) : le champ source est celui émis par un point de tir et les données de réflexion sont les enregistrements réalisés par le (ou les) capteur(s) associé(s) à ce point de tir ayant l'offset considéré ; dans une telle nûgration il faut, pour réaliser la migration des données associées à un déport, effectuer autant de modélisation la de propagation et de retropropagation d'ondes qu'il y a de points de tir et sommer les résultats obtenus pour chaque point de tir.
Des exemples d'utilisation de ce type de techniques sont décrits par exemple dans :
- Claerbout, J.F., 1985 ; Imaging the Earth's interior ; Blackwell Publications,
La présente invention concerne une méthode pour réaliser une migration avant sommation d'événements sismiques enregistrés pour l'imagerie d'une partie d'une zone souterraine.
La méthode selon l'invention permet par exemple de réaliser une migration profondeur de type 3D avant sommation ( 3D prestack depth migration ), pour un modèle de vitesses donné, de manière à imager les interfaces géologiques ou hétérogénéités diverses d'une partie de la subsurface.
Art antérieur La migration avant sommation est une méthode classique de traitement de données sismiques. De façon générale la technique consiste, connaissant la valeur d'un champ d'ondes à une profondeur connue, par exemple en surface, ainsi qu'un modèle.
de la répartition des vitesses de propagation des ondes dans la zone, à modéliser la propagation du champ source et la retropropagation des données de réflexion enregistrées et à chercher des cohérences de phase entre ces deux champs modélisés.
On distingue trois grands types de migration avant sommation :
- la migration par point de tir : le champ source est l'état vibratoire engendré par le point de tir et les données de réflexion sont la réponse de la subsurface à ce champ source ;
- la migration par ondes planes aussi appelée migration par angle d'éclairage commun : le champ source est l'onde plane considérée et les données de réflexion sont la réponse de la subsurface à ce champ source ;
- la migration par déport ( offset ) : le champ source est celui émis par un point de tir et les données de réflexion sont les enregistrements réalisés par le (ou les) capteur(s) associé(s) à ce point de tir ayant l'offset considéré ; dans une telle nûgration il faut, pour réaliser la migration des données associées à un déport, effectuer autant de modélisation la de propagation et de retropropagation d'ondes qu'il y a de points de tir et sommer les résultats obtenus pour chaque point de tir.
Des exemples d'utilisation de ce type de techniques sont décrits par exemple dans :
- Claerbout, J.F., 1985 ; Imaging the Earth's interior ; Blackwell Publications,
2 - Duquet, B., 1996 ; Amélioration de l'Imagerie Sismique de Structures Géologiques Complexes ; thèse, Université Paris 13, ou - Whitmore, N.D., Felinsky, W.F., Murphy, G.E. and Gray, S.H., 1993 ; The Application of Common Offset and Common Angle Pre-stack Depth Migration in the North Sea, 55th Mtg., EAGE, Expanded abstract.
Les implémentations classiques basées sur l'intégrale de Kirchhoff (ou des versions plus élaborées de cette technique, elle-même basée sur des techniques asymptotiques hautes fréquences), ont pour inconvénient majeur d'être généralement très coûteuses en temps de calcul, en raison du volume aussi bien des données à traiter que des résultats, 1o surtout quand le champ de vitesses varie latéralement (ce qui complique les calculs de temps d'arrivée requis pour la mise en ceuvre de cette méthode). Pour des raisons économiques, on est souvent amené à limiter les volumes de données, (par décimation) et/ou la quantité de résultats produits (volume imagé de taille réduite, échantillonnage grossier des résultats).
Les implémentations classiques basées sur l'intégrale de Kirchhoff (ou des versions plus élaborées de cette technique, elle-même basée sur des techniques asymptotiques hautes fréquences), ont pour inconvénient majeur d'être généralement très coûteuses en temps de calcul, en raison du volume aussi bien des données à traiter que des résultats, 1o surtout quand le champ de vitesses varie latéralement (ce qui complique les calculs de temps d'arrivée requis pour la mise en ceuvre de cette méthode). Pour des raisons économiques, on est souvent amené à limiter les volumes de données, (par décimation) et/ou la quantité de résultats produits (volume imagé de taille réduite, échantillonnage grossier des résultats).
3 Définition de la méthode La méthode selon l'invention permet de réaliser la migration de données sismiques pour l'imagerie d'une partie d'une zone souterraine, les données sismiques étant obtenues à
l'issue d'une série de NS cycles de sismique réflexion, qui comprend chacun l'émission successive de champs d'onde élémentaires défini chacun par l'association d'un signal ~
sismique W(t) et d'un lieu d'émission dans une série de lieux d'émission Si avec 1:5 i<_ Ns, la réception par des récepteurs sismiques placés en des positions des signaux sismiques renvoyés par la zone en réponse à chacun de ces champs d'onde, et l'enregistrement des différents signaux reçus par chaque récepteur sismique sous la forme l0 de traces sismiques d(t) dépendant du temps.
Elle est caractérisée en ce que, pour un modèle de vitesses donné, elle comporte l'ensemble des étapes suivantes ~
a) on définit un vecteur de lenteur p dont les deux composantes pi et p2 peuvent chacune prendre une suite de valeurs préalablement définie ;
-~ -~
b) on définit pour un vecteur de lenteur p donné et lieu d'émission Si donné, une -->
fonction de décalage temporel to ( p , i) ~
c) on applique la fonction de décalage temporel to( p, i) à chaque champ d'onde ~
élémentaire (associé au lieu d'émission Si ) et on forme un premier champ d'ondes composite en surface par superposition spatio-temporelle des différents champs d'onde élémentaires ainsi décalés ;
d) on applique un décalage temporel to ( p, i) à chaque trace sismique d~(t) repérée i par le couple (i,j) et on forme un champ de traces composite en surface par superposition spatio-temporelle des différentes traces sismiques ainsi décalées
l'issue d'une série de NS cycles de sismique réflexion, qui comprend chacun l'émission successive de champs d'onde élémentaires défini chacun par l'association d'un signal ~
sismique W(t) et d'un lieu d'émission dans une série de lieux d'émission Si avec 1:5 i<_ Ns, la réception par des récepteurs sismiques placés en des positions des signaux sismiques renvoyés par la zone en réponse à chacun de ces champs d'onde, et l'enregistrement des différents signaux reçus par chaque récepteur sismique sous la forme l0 de traces sismiques d(t) dépendant du temps.
Elle est caractérisée en ce que, pour un modèle de vitesses donné, elle comporte l'ensemble des étapes suivantes ~
a) on définit un vecteur de lenteur p dont les deux composantes pi et p2 peuvent chacune prendre une suite de valeurs préalablement définie ;
-~ -~
b) on définit pour un vecteur de lenteur p donné et lieu d'émission Si donné, une -->
fonction de décalage temporel to ( p , i) ~
c) on applique la fonction de décalage temporel to( p, i) à chaque champ d'onde ~
élémentaire (associé au lieu d'émission Si ) et on forme un premier champ d'ondes composite en surface par superposition spatio-temporelle des différents champs d'onde élémentaires ainsi décalés ;
d) on applique un décalage temporel to ( p, i) à chaque trace sismique d~(t) repérée i par le couple (i,j) et on forme un champ de traces composite en surface par superposition spatio-temporelle des différentes traces sismiques ainsi décalées
4 e) on effectue une migration du champ de traces composite en utilisant comme champ d'ondes le champ d'ondes composite, ceci en modélisant la propagation de champ d'ondes composite ainsi que la retropropagation de champ de traces composite, et en combinant de façon adaptée les deux champs composites ainsi modélisés en tout point de la zone à imager.
f) On répète les étapes c) à e) pour toutes les valeurs prises par les composantes pl et p2 du dit vecteur p ; et ~
g) pour toute valeur fixée de la deuxième composante pZ du dit vecteur p, on somme les résultats de ces différentes combinaisons de manière à obtenir une image migrée associée à cette valeur fixée de pZ, réalisant ainsi une migration avant sommation.
Suivant un mode de mise en oeuvre, on peut sommer les résultats obtenus en g) pour toutes les valeurs prises par le paramètre pZ. On effectue ainsi un sommation après migration.
Suivant un mode de mise en ceuvre, il est possible de réaliser directement le sommation après migration, sans particulariser l'étape g).
La méthode peut comporter en outre une mise a jour des vitesses par analyse des ~
déformations obtenues lorsque l'on fait varier la deuxième coordonnée P2 du vecteur p.
Suivant un mode de mise en oeuvre, on peut former une image migrée d'une partie de la zone à imager en exploitant le phénomène de conversion d'ondes, par définition d'au moins une partie du champ de vitesse en ondes P et en ondes S (en appliquant par exemple au préalable un pré-traitement adapté aux données, visant à séparer les différents types d'événements sismique).
On peut utiliser les étapes a) à g) pour la détermination du gradient d'une fonction coût intervenant dans un problème inverse de sismique.
Il est possible également de remplacer une migration profondeur par une migration temps.
La méthode proposée présente de grands avantages :
1) Elle permet de réaliser une migration à un prix (coût de calcul) avantageux car indépendant du volume de résultats calculé et du nombre de traces sismiques enregistrées, contrairement aux méthodes classiques de type Kirchhoff. Seul le volume de la zone dans
f) On répète les étapes c) à e) pour toutes les valeurs prises par les composantes pl et p2 du dit vecteur p ; et ~
g) pour toute valeur fixée de la deuxième composante pZ du dit vecteur p, on somme les résultats de ces différentes combinaisons de manière à obtenir une image migrée associée à cette valeur fixée de pZ, réalisant ainsi une migration avant sommation.
Suivant un mode de mise en oeuvre, on peut sommer les résultats obtenus en g) pour toutes les valeurs prises par le paramètre pZ. On effectue ainsi un sommation après migration.
Suivant un mode de mise en ceuvre, il est possible de réaliser directement le sommation après migration, sans particulariser l'étape g).
La méthode peut comporter en outre une mise a jour des vitesses par analyse des ~
déformations obtenues lorsque l'on fait varier la deuxième coordonnée P2 du vecteur p.
Suivant un mode de mise en oeuvre, on peut former une image migrée d'une partie de la zone à imager en exploitant le phénomène de conversion d'ondes, par définition d'au moins une partie du champ de vitesse en ondes P et en ondes S (en appliquant par exemple au préalable un pré-traitement adapté aux données, visant à séparer les différents types d'événements sismique).
On peut utiliser les étapes a) à g) pour la détermination du gradient d'une fonction coût intervenant dans un problème inverse de sismique.
Il est possible également de remplacer une migration profondeur par une migration temps.
La méthode proposée présente de grands avantages :
1) Elle permet de réaliser une migration à un prix (coût de calcul) avantageux car indépendant du volume de résultats calculé et du nombre de traces sismiques enregistrées, contrairement aux méthodes classiques de type Kirchhoff. Seul le volume de la zone dans
5 lequel les ondes se propagent, ayant une incidence sur le coût de calcul. On obtient ainsi des images volumiques en prenant en compte l'ensemble des traces sismiques à
un coût avantageux. On estime que cette méthode diminue d'un facteur de l'ordre de quelques dizaines le temps de calcul nécessaire à la réalisation de la dite migration 3D avant sommation.
2) La méthode s'applique pour des modèles de vitesses d'une complexité
arbitraire pour peu que la notion de migration avant sommation conserve un sens. Elle s'applique sans rencontrer aucune des limitations propres aux techniques asymptotiques hautes fréquences (optique géométrique) habituellement utilisées pour réaliser les migrations 3D
avant sommation.
La méthode peut être mise en ceuvre à l'aide d'outils classiques de modélisation de la propagation et de la retropropagation des ondes, que l'on trouve décrits par exemple dans la référence Duquet B. déjà citée.
Par application de la méthode, on peut obtenir des images migrées élémentaires associées à une valeur donnée d'un paramètre et la somme de ces images ( post migration stack ), et ceci aussi bien dans le domaine profondeur que dans le domaine temps.
D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à
la lecture de la description ci-après d'un exemple non limitatif de réalisation, en se référant aux dessins annexés où :
- la Fig.1 montre de façon très schématique, une disposition de lieux d'émission S; et de lieux de réception sismique R;' d'une zone souterraine explorée où l'on cherche à
restituer des horizons In de la subsurface;
- la Fig.2 montre une section sismique extraite d'un bloc de données 3D
migrées profondeur après sommation, obtenue dans le cadre de l'exploration sismique d'une structure salifère de mer du Nord ;
un coût avantageux. On estime que cette méthode diminue d'un facteur de l'ordre de quelques dizaines le temps de calcul nécessaire à la réalisation de la dite migration 3D avant sommation.
2) La méthode s'applique pour des modèles de vitesses d'une complexité
arbitraire pour peu que la notion de migration avant sommation conserve un sens. Elle s'applique sans rencontrer aucune des limitations propres aux techniques asymptotiques hautes fréquences (optique géométrique) habituellement utilisées pour réaliser les migrations 3D
avant sommation.
La méthode peut être mise en ceuvre à l'aide d'outils classiques de modélisation de la propagation et de la retropropagation des ondes, que l'on trouve décrits par exemple dans la référence Duquet B. déjà citée.
Par application de la méthode, on peut obtenir des images migrées élémentaires associées à une valeur donnée d'un paramètre et la somme de ces images ( post migration stack ), et ceci aussi bien dans le domaine profondeur que dans le domaine temps.
D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à
la lecture de la description ci-après d'un exemple non limitatif de réalisation, en se référant aux dessins annexés où :
- la Fig.1 montre de façon très schématique, une disposition de lieux d'émission S; et de lieux de réception sismique R;' d'une zone souterraine explorée où l'on cherche à
restituer des horizons In de la subsurface;
- la Fig.2 montre une section sismique extraite d'un bloc de données 3D
migrées profondeur après sommation, obtenue dans le cadre de l'exploration sismique d'une structure salifère de mer du Nord ;
6 - la Fig.3 montre une section sismique, coïncidant avec la section sismique de la Fig. 2, extraite d'un bloc de données 3D migrées profondeur avant sommation par la méthode selon l'invention, en utilisant le modèle de vitesse utilisé pour l'obtention de la Fig. 2 ; et - la Fig.4 montre en quatre points différents de la surface de la zone souterraine, les images 3D migrées en profondeur avant sommation, par la méthode selon l'invention ~
lorsque la deuxième composante P2 du vecteur p varie sur son domaine de définition, l'analyse de la déformation des événements observés sur ces images permettant la mise à jour du modèle de vitesse.
DESCRIPTION DETAILLEE DE LA METHODE
Définitions, notations et hypothèses requises La méthode selon l'invention permet de réaliser la migration de données sismiques pour l'imagerie d'une zone souterraine M, les données sismiques étant obtenues à l'issue d'une série de NS cycles sismiques, chacun d'eux (Fig.1) comportant l'émission d'un signal ~
sismique W(t) depuis un lieu d'émission Si avec 1:9 i<_ Ns, la réception par des récepteurs sismiques placés en des positions des signaux sismiques renvoyés par les discontinuités de la zone et l'enregistrement des différents signaux reçus par chaque récepteur sismique sous la forme d'une trace sismique d(t). On admet que chaque source sismique émet le même signal W(t), car on peut toujours se ramener à une telle situation par un pré-traitement adapté des données, tel qu'une déconvolution.
On désignera par profil un ensemble de lieux d'émission alignés et l'on supposera que toutes les acquisitions de signaux sont effectuées en utilisant des profils d'émission parallèles. On supposera que, pour chaque source, les récepteurs sont localisés sur le profil ~
où elle se trouve (repérée par le lieu d'émission Si ). En fait, la méthode s'avère robuste vis à vis de la précision avec laquelle ces hypothèses sont vérifiées. Les hypothèses précédentes sont faites principalement pour pouvoir introduire les notations ci-dessous.
On définit un repère orthonormé dont le premier axe est perpendiculaire à la direction du profil et le deuxième axe parallèle à cette direction et un vecteur un vecteur
lorsque la deuxième composante P2 du vecteur p varie sur son domaine de définition, l'analyse de la déformation des événements observés sur ces images permettant la mise à jour du modèle de vitesse.
DESCRIPTION DETAILLEE DE LA METHODE
Définitions, notations et hypothèses requises La méthode selon l'invention permet de réaliser la migration de données sismiques pour l'imagerie d'une zone souterraine M, les données sismiques étant obtenues à l'issue d'une série de NS cycles sismiques, chacun d'eux (Fig.1) comportant l'émission d'un signal ~
sismique W(t) depuis un lieu d'émission Si avec 1:9 i<_ Ns, la réception par des récepteurs sismiques placés en des positions des signaux sismiques renvoyés par les discontinuités de la zone et l'enregistrement des différents signaux reçus par chaque récepteur sismique sous la forme d'une trace sismique d(t). On admet que chaque source sismique émet le même signal W(t), car on peut toujours se ramener à une telle situation par un pré-traitement adapté des données, tel qu'une déconvolution.
On désignera par profil un ensemble de lieux d'émission alignés et l'on supposera que toutes les acquisitions de signaux sont effectuées en utilisant des profils d'émission parallèles. On supposera que, pour chaque source, les récepteurs sont localisés sur le profil ~
où elle se trouve (repérée par le lieu d'émission Si ). En fait, la méthode s'avère robuste vis à vis de la précision avec laquelle ces hypothèses sont vérifiées. Les hypothèses précédentes sont faites principalement pour pouvoir introduire les notations ci-dessous.
On définit un repère orthonormé dont le premier axe est perpendiculaire à la direction du profil et le deuxième axe parallèle à cette direction et un vecteur un vecteur
7 ~
lenteur p(homogène à l'inverse d'une vitesse, comme il est bien connu). Dont les composantes pl et P2 mesurent respectivement les composantes suivant ces deux axes de ce vecteur lenteur p . Pour un vecteur lenteur p donné et une source placée en un lieu ~
d'émission Si donné, on définit une fonction de décalage temporel :
-~ -~ -3 -~
(1) to(p,i)= p.(Si -So) ~
où So représente un point quelconque du domaine d'acquisition.
Traitement ~
Sélection du vecteur p On commence par choisir un ensemble fini Pl de valeurs prises par le paramètre pl lo et un ensemble fini P2 de valeurs prises par le paramètre P2. L'ensemble E=P1xP2 sera ~
l'ensemble des valeurs (vectorielles) prises par le vecteur p. L'ensemble P1 des valeurs que prendra le paramètre pl pourra par exemple être construit en échantillonnant l'intervalle [-plmin, plmax] avec un pas d'échantillonnage régulier Spl. La valeur à donner au pas Spldépend notamment de la précision souhaitée et de l'espacement entre les profils.
Une valeur typique est : Spl = 2.5 . 10-5 s/m. Les valeurs à donner à plmin et pimax dépendent de la complexité de la structure dans la direction orthogonale au profil. Des valeurs typiques sont :-plmin = plmax = 2.5 . 10-4 s/m. L'ensemble P2 des valeurs que prendra le paramètre pz pourra par exemple être construit en échantillonnant l'intervalle [-pZmin, p2max] avec un pas d'échantillonnage régulier Spz. La valeur à donner à
Sp2 dépend 2o notamment de la précision souhaitée et de la finesse avec laquelle on veut suivre l'évolution, lorsque pZ varie, des événements dans les collections point image. Une valeur typique est : 8p2 =2.5 . 10-5 s/m. Les valeurs à donner à p2min et p2max dépendent de la complexité de la structure dans la direction des profils. Des valeurs typiques sont :-p2min = pzmax = 2.5 . 10-4 s/m A) Etapes de traitement :
lenteur p(homogène à l'inverse d'une vitesse, comme il est bien connu). Dont les composantes pl et P2 mesurent respectivement les composantes suivant ces deux axes de ce vecteur lenteur p . Pour un vecteur lenteur p donné et une source placée en un lieu ~
d'émission Si donné, on définit une fonction de décalage temporel :
-~ -~ -3 -~
(1) to(p,i)= p.(Si -So) ~
où So représente un point quelconque du domaine d'acquisition.
Traitement ~
Sélection du vecteur p On commence par choisir un ensemble fini Pl de valeurs prises par le paramètre pl lo et un ensemble fini P2 de valeurs prises par le paramètre P2. L'ensemble E=P1xP2 sera ~
l'ensemble des valeurs (vectorielles) prises par le vecteur p. L'ensemble P1 des valeurs que prendra le paramètre pl pourra par exemple être construit en échantillonnant l'intervalle [-plmin, plmax] avec un pas d'échantillonnage régulier Spl. La valeur à donner au pas Spldépend notamment de la précision souhaitée et de l'espacement entre les profils.
Une valeur typique est : Spl = 2.5 . 10-5 s/m. Les valeurs à donner à plmin et pimax dépendent de la complexité de la structure dans la direction orthogonale au profil. Des valeurs typiques sont :-plmin = plmax = 2.5 . 10-4 s/m. L'ensemble P2 des valeurs que prendra le paramètre pz pourra par exemple être construit en échantillonnant l'intervalle [-pZmin, p2max] avec un pas d'échantillonnage régulier Spz. La valeur à donner à
Sp2 dépend 2o notamment de la précision souhaitée et de la finesse avec laquelle on veut suivre l'évolution, lorsque pZ varie, des événements dans les collections point image. Une valeur typique est : 8p2 =2.5 . 10-5 s/m. Les valeurs à donner à p2min et p2max dépendent de la complexité de la structure dans la direction des profils. Des valeurs typiques sont :-p2min = pzmax = 2.5 . 10-4 s/m A) Etapes de traitement :
8 1) Création d'un premier champ d'ondes composite en surface (champ de sources) ~
L'ondelette W(t) attachée à chacune source Si définit un champ d'onde élémentaire W(t).S( x"p - Si ) où le vecteur x 2D repère la position d'un point quelconque à la surface du sol et S est une masse de Dirac définie sur R' et centrée à l'origine.
Après avoir retardé
~
du temps to( p , i) le champ d'onde élémentaire associé, on forme un premier champ d'ondes composite en surface par superposition spatio-temporelle des différentes champs d'onde élémentaires ainsi décalées. On définit ainsi la fonction :
~ -~ -~ -~
(2) WP (x2D,t)=~S(x2D-Si)W(t-to(P, i)) 2) Création d'un deuxième champ d'ondes composite en surface (champ de traces) Apres avoir retardé du temps to( p, i) la trace sismique d~ i (t) associée à
chaque couple (lieu d'émission Si, lieu de réception R~) repéré par le couple (i,j), on forme un deuxième champ d'ondes composite en surface par superposition spatio-temporelle des différentes traces sismiques ainsi décalées. Avec les mêmes notations, on définit ainsi la fonction :
(3) D' ( x 2D, t) _~ S( x 2D - R 1) d (t - ta( P+ i)) P. J i 3) Modélisation de la propaaation du champ source composite) On modélise la propagation du champ source composite en recherchant des solutions périodiques de période T de l'équation des ondes, en utilisant comme distribution de vitesses celle définie par le modèle de vitesses considéré. On obtient ainsi un prenver --~
champ d'ondes propagé (dépendant de l'espace et du temps) W-~ ( x 3D, t), ceci pour tout P.
temps t (la solution est périodique en temps) et pour tout point de la partie de subsurface ~
que l'on cherche à imager, point dont la position est repérée par le vecteur x 3D. La période
L'ondelette W(t) attachée à chacune source Si définit un champ d'onde élémentaire W(t).S( x"p - Si ) où le vecteur x 2D repère la position d'un point quelconque à la surface du sol et S est une masse de Dirac définie sur R' et centrée à l'origine.
Après avoir retardé
~
du temps to( p , i) le champ d'onde élémentaire associé, on forme un premier champ d'ondes composite en surface par superposition spatio-temporelle des différentes champs d'onde élémentaires ainsi décalées. On définit ainsi la fonction :
~ -~ -~ -~
(2) WP (x2D,t)=~S(x2D-Si)W(t-to(P, i)) 2) Création d'un deuxième champ d'ondes composite en surface (champ de traces) Apres avoir retardé du temps to( p, i) la trace sismique d~ i (t) associée à
chaque couple (lieu d'émission Si, lieu de réception R~) repéré par le couple (i,j), on forme un deuxième champ d'ondes composite en surface par superposition spatio-temporelle des différentes traces sismiques ainsi décalées. Avec les mêmes notations, on définit ainsi la fonction :
(3) D' ( x 2D, t) _~ S( x 2D - R 1) d (t - ta( P+ i)) P. J i 3) Modélisation de la propaaation du champ source composite) On modélise la propagation du champ source composite en recherchant des solutions périodiques de période T de l'équation des ondes, en utilisant comme distribution de vitesses celle définie par le modèle de vitesses considéré. On obtient ainsi un prenver --~
champ d'ondes propagé (dépendant de l'espace et du temps) W-~ ( x 3D, t), ceci pour tout P.
temps t (la solution est périodique en temps) et pour tout point de la partie de subsurface ~
que l'on cherche à imager, point dont la position est repérée par le vecteur x 3D. La période
9 T est choisie comme il est habituel dans les algorithmes classiques de migration (c'est à
dire de l'ordre de la durée d'enregistrement).
4) Modélisation de la rétropropagation du champ de traces composite On modélise la rétropropagation du champ de traces composite n recherchant des solutions périodiques de période T de l'équation des ondes en utilisant comme distribution de vitesses celle définie par le modèle de vitesses considéré. On obtient ainsi un champ d'ondes rétropropagé (dépendant de l'espace et du temps) D' ( x 3D, t), ceci pour tout P.
temps t et pour tout point de la partie de la subsurface que l'on cherche à
imager, point dont -~
la position est repérée par le vecteur X 3D=
5) Recherche de cohérence de phase On recherche alors une éventuelle cohérence de phase (par exemple par des calculs d'intercorrélation) entre le premier champ d'ondes composite propagé et le deuxième champ de traces composite rétropropagé, ceci en tout point de la zone souterraine (subsurface) que l'on cherche à imager, point dont la position est repérée par le vecteur ~
x 3D. Pour ce faire, on évalue, toujours dans le cas d'un calcul par intercorrélation, la quantité
-~ T -~ -~
(4) M-4 ( X 3D) = J W~ ( X 3D, t) D-' ( X 3D, t) dt.
P. p P-La composante P2 ayant une valeur donnée dans P2, on somme les résultats calculés en (4) lorsque le paramètre pi parcourt Pt. Cette sommation permet de définir, pour tout point image repéré par le vecteur X 3D, la fonction Mp2 ( x 3D) par la formule :
(5) M p2 ( X 3D) M- ( X 3D) pi E Pi P.
~
La quantité (5) s'interprète comme la valeur au point x 3D de la superposition sur les différents profils d'acquisition des images migrées associées à une onde cylindrique ayant pour axe le profil d'acquisition et une pente définie par la lenteur p2.
B) Sommation après migration Il est possible en outre de réaliser une sommation après migration en sommant les contributions (5) obtenues pour les différentes valeurs du paramètre p2 appartenant à P2, selon la formule :
5 (6) M(X 3D) = E M p2 ( X 3D) p, E P2 VARIANTES
1) Une première variante des étapes A-3 et A-4 décrites ci-dessus consiste à
résoudre à l'étape A-3 un problème de valeurs initiales et en A-4 un problème de valeurs finales au lieu de rechercher des solutions périodiques (procédures d'imagerie classiques 1o connues sous le nom de migration par renversement du temps).
2) Une autre variante consiste à remplacer l'étape A-1 de création du premier champ d'ondes composite par la création d'une onde plane dont les lenteurs suivant la direction des profils et suivant la direction orthogonale sont respectivement P2 et pl. On aboutit ainsi à la formule --~ ~ -~ --(7) W---> (x3D)=W(t- p.(X2D-So)) p Cette variante met en évidence une certaine similitude des traitements mis en oyuvre dans la méthode et de ceux mis en oeuvre dans la méthode classique de migration par ondes planes. Une différence essentielle qui fait l'originalité de la méthode réside dans la possibilité de réaliser une procédure de migration par ondes planes même dans le cas où
l'acquisition ne permet pas la synthèse de la réponse de la subsurface à une excitation ondes planes, réponse qu'il est indispensable de connaître dans les algorithmes connus de migration par ondes planes.
Autres applications Suivant un mode de mise en ceuvre, il est possible de former une image migrée d'une partie de la subsurface en exploitant le phénomène de conversion d'ondes, par définition d'au moins une partie du champ de vitesses en ondes P et en ondes S, ceci après éventuellement appliqué un prétraitement adapté aux données, visant à séparer les différents types d'événements sismiques.
Les étapes précédemment définies peuvent aussi servir pour le calcul du gradient d'une fonction coût intervenant dans un problème inverse de sismique.
dire de l'ordre de la durée d'enregistrement).
4) Modélisation de la rétropropagation du champ de traces composite On modélise la rétropropagation du champ de traces composite n recherchant des solutions périodiques de période T de l'équation des ondes en utilisant comme distribution de vitesses celle définie par le modèle de vitesses considéré. On obtient ainsi un champ d'ondes rétropropagé (dépendant de l'espace et du temps) D' ( x 3D, t), ceci pour tout P.
temps t et pour tout point de la partie de la subsurface que l'on cherche à
imager, point dont -~
la position est repérée par le vecteur X 3D=
5) Recherche de cohérence de phase On recherche alors une éventuelle cohérence de phase (par exemple par des calculs d'intercorrélation) entre le premier champ d'ondes composite propagé et le deuxième champ de traces composite rétropropagé, ceci en tout point de la zone souterraine (subsurface) que l'on cherche à imager, point dont la position est repérée par le vecteur ~
x 3D. Pour ce faire, on évalue, toujours dans le cas d'un calcul par intercorrélation, la quantité
-~ T -~ -~
(4) M-4 ( X 3D) = J W~ ( X 3D, t) D-' ( X 3D, t) dt.
P. p P-La composante P2 ayant une valeur donnée dans P2, on somme les résultats calculés en (4) lorsque le paramètre pi parcourt Pt. Cette sommation permet de définir, pour tout point image repéré par le vecteur X 3D, la fonction Mp2 ( x 3D) par la formule :
(5) M p2 ( X 3D) M- ( X 3D) pi E Pi P.
~
La quantité (5) s'interprète comme la valeur au point x 3D de la superposition sur les différents profils d'acquisition des images migrées associées à une onde cylindrique ayant pour axe le profil d'acquisition et une pente définie par la lenteur p2.
B) Sommation après migration Il est possible en outre de réaliser une sommation après migration en sommant les contributions (5) obtenues pour les différentes valeurs du paramètre p2 appartenant à P2, selon la formule :
5 (6) M(X 3D) = E M p2 ( X 3D) p, E P2 VARIANTES
1) Une première variante des étapes A-3 et A-4 décrites ci-dessus consiste à
résoudre à l'étape A-3 un problème de valeurs initiales et en A-4 un problème de valeurs finales au lieu de rechercher des solutions périodiques (procédures d'imagerie classiques 1o connues sous le nom de migration par renversement du temps).
2) Une autre variante consiste à remplacer l'étape A-1 de création du premier champ d'ondes composite par la création d'une onde plane dont les lenteurs suivant la direction des profils et suivant la direction orthogonale sont respectivement P2 et pl. On aboutit ainsi à la formule --~ ~ -~ --(7) W---> (x3D)=W(t- p.(X2D-So)) p Cette variante met en évidence une certaine similitude des traitements mis en oyuvre dans la méthode et de ceux mis en oeuvre dans la méthode classique de migration par ondes planes. Une différence essentielle qui fait l'originalité de la méthode réside dans la possibilité de réaliser une procédure de migration par ondes planes même dans le cas où
l'acquisition ne permet pas la synthèse de la réponse de la subsurface à une excitation ondes planes, réponse qu'il est indispensable de connaître dans les algorithmes connus de migration par ondes planes.
Autres applications Suivant un mode de mise en ceuvre, il est possible de former une image migrée d'une partie de la subsurface en exploitant le phénomène de conversion d'ondes, par définition d'au moins une partie du champ de vitesses en ondes P et en ondes S, ceci après éventuellement appliqué un prétraitement adapté aux données, visant à séparer les différents types d'événements sismiques.
Les étapes précédemment définies peuvent aussi servir pour le calcul du gradient d'une fonction coût intervenant dans un problème inverse de sismique.
Claims (8)
1) Méthode pour réaliser une migration avant sommation d'événements sismiques pour l'imagerie d'une partie d'une zone souterraine à partir d'une série d'un nombre N s de cycles de sismique réflexion, qui comprend l'émission successive de champs d'onde élémentaires défini chacun par l'association d'un signal sismique W(t) et d'un lieu d'émission défini dans une série de lieux d'émission S1 avec 1<=
i<= N s, la réception par des récepteurs sismiques placés en des positions des signaux sismiques renvoyés par la zone en réponse à chacun de ces champs d'onde, et l'enregistrement des différents signaux reçus par chaque récepteur sismique sous la forme de traces sismiques d~(t) dépendant du temps, caractérisée en ce que, pour un modèle de vitesses donné, elle comporte l'ensemble des étapes suivantes :
a) on définit un vecteur de lenteur p dont les deux composantes p1 et P2 peuvent chacune prendre une suite de valeurs préalablement définie b) on définit pour un vecteur de lenteur p donné et lieu d'émission Si donné, une fonction de décalage temporel t o(p, i) c) on applique la fonction de décalage temporel t o( p, i) à chaque champ d'onde élémentaire associé au lieu d'émission S, et on forme un champ d'ondes composite en surface par superposition spatio-temporelle des différents champs d'onde élémentaires ainsi décalés d) on applique un décalage temporel t o(p , i) à chaque trace sismique d~(t) repérée par le couple (i,j) et on forme un champ de traces composite en surface par superposition spatio-temporelle des différentes traces sismiques ainsi décalées ;
e) on effectue une migration du champ de traces composite en utilisant comme champ d'ondes le champ d'ondes composite, ceci en modélisant la propagation de champ d'ondes composite ainsi que la retropropagation de champ de traces composite, et en combinant de façon adaptée les deux champs composites ainsi modélisés en tout point de la zone à imager.
f) On répète les étapes c) à e) pour toutes les valeurs prises par les composantes p1 et p2 du dit vecteur ~ ; et g) pour toute valeur fixée de la deuxième composante p2 du dit vecteur ~, on somme les résultats de ces différentes combinaisons de manière à obtenir une image migrée associée à cette valeur fixée de p2, réalisant ainsi une migration avant sommation.
i<= N s, la réception par des récepteurs sismiques placés en des positions des signaux sismiques renvoyés par la zone en réponse à chacun de ces champs d'onde, et l'enregistrement des différents signaux reçus par chaque récepteur sismique sous la forme de traces sismiques d~(t) dépendant du temps, caractérisée en ce que, pour un modèle de vitesses donné, elle comporte l'ensemble des étapes suivantes :
a) on définit un vecteur de lenteur p dont les deux composantes p1 et P2 peuvent chacune prendre une suite de valeurs préalablement définie b) on définit pour un vecteur de lenteur p donné et lieu d'émission Si donné, une fonction de décalage temporel t o(p, i) c) on applique la fonction de décalage temporel t o( p, i) à chaque champ d'onde élémentaire associé au lieu d'émission S, et on forme un champ d'ondes composite en surface par superposition spatio-temporelle des différents champs d'onde élémentaires ainsi décalés d) on applique un décalage temporel t o(p , i) à chaque trace sismique d~(t) repérée par le couple (i,j) et on forme un champ de traces composite en surface par superposition spatio-temporelle des différentes traces sismiques ainsi décalées ;
e) on effectue une migration du champ de traces composite en utilisant comme champ d'ondes le champ d'ondes composite, ceci en modélisant la propagation de champ d'ondes composite ainsi que la retropropagation de champ de traces composite, et en combinant de façon adaptée les deux champs composites ainsi modélisés en tout point de la zone à imager.
f) On répète les étapes c) à e) pour toutes les valeurs prises par les composantes p1 et p2 du dit vecteur ~ ; et g) pour toute valeur fixée de la deuxième composante p2 du dit vecteur ~, on somme les résultats de ces différentes combinaisons de manière à obtenir une image migrée associée à cette valeur fixée de p2, réalisant ainsi une migration avant sommation.
2) Méthode selon la revendication 1, caractérisée en ce qu'elle comporte en outre la formation d'une image finale représentative d'une sommation après migration des traces obtenues pour l'ensemble des valeurs de la deuxième composante p2 du vecteur ~.
3) Méthode selon la revendication 1, caractérisée en ce que l'on utilise les étapes a) à g) pour réaliser directement un sommation après migration.
4) Méthode selon la revendication 1, caractérisée en ce qu'elle comporte en outre une mise a jour des vitesses par analyse des déformations obtenues lorsque l'on fait varier la deuxième coordonnée p2 du vecteur ~.
5) Méthode selon l'une quelconque des revendications 1 à 4, caractérisée en ce que l'on forme une image migrée d'une partie de la zone à imager en exploitant le phénomène de conversion d'ondes, par définition d'au moins une partie du champ de vitesse en ondes P et en ondes S.
6) Méthode selon la revendication 5, caractérisé en ce que l'on applique un prétraitement adapté aux données, visant à séparer les différents types d'événements sismiques.
7) Méthode selon l'une quelconque des revendications 1 à 6, caractérisée en ce que l'on utilise lesdites étapes a) à g) pour la détermination du gradient d'une fonction coût intervenant dans un problème inverse de sismique.
8) Méthode selon l'une quelconque des revendications 1 à 7, caractérisée par le fait que l'on remplace une migration profondeur par une migration temps.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9812337 | 1998-10-01 | ||
FR9812337A FR2784195B1 (fr) | 1998-10-01 | 1998-10-01 | Methode pour realiser en 3d avant sommation, une migration de donnees sismiques |
Publications (2)
Publication Number | Publication Date |
---|---|
CA2284386A1 CA2284386A1 (fr) | 2000-04-01 |
CA2284386C true CA2284386C (fr) | 2008-08-05 |
Family
ID=9531110
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA002284386A Expired - Fee Related CA2284386C (fr) | 1998-10-01 | 1999-09-29 | Methode pour realiser en 3d avant sommation, une migration de donnees sismiques |
Country Status (6)
Country | Link |
---|---|
US (1) | US6311133B1 (fr) |
BR (1) | BR9904428B1 (fr) |
CA (1) | CA2284386C (fr) |
FR (1) | FR2784195B1 (fr) |
GB (1) | GB2342165B (fr) |
NO (1) | NO331526B1 (fr) |
Families Citing this family (46)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6574564B2 (en) * | 1998-10-01 | 2003-06-03 | Institut Francais Du Petrole | 3D prestack seismic data migration method |
WO2005010797A2 (fr) * | 2003-07-23 | 2005-02-03 | Lee Wook B | Modelisation de vitesses en 3d amelioree avec etalonnage et ajustement de tendance par des techniques geostatique, particulierement avantageuse pour la migration en duree de pre-sommation a rayons courbes et pour une telle migration suivie par une migration en profondeur de pre-sommation |
FR2902196B1 (fr) * | 2006-06-09 | 2008-08-22 | Cie Generale De Geophysique Sa | Procede de traitement sismique du type mettant en oeuvre une migration profondeur avant addition |
NZ575497A (en) | 2006-09-28 | 2011-08-26 | Exxonmobil Upstream Res Co | Iterative inversion of data from simultaneous geophysical sources |
CA2703588C (fr) | 2007-12-12 | 2015-12-01 | Exxonmobil Upstream Research Company | Procede et appareil d'evaluation de formations sous-marines |
BRPI0908156A2 (pt) | 2008-03-21 | 2015-08-11 | Exxonmobil Upstream Company | Métodos para invesão de dados geofísicos medidos e para produzir hidrocarbonetos de uma região de subsuperfície |
RU2503035C2 (ru) * | 2008-08-11 | 2013-12-27 | Эксонмобил Апстрим Рисерч Компани | Оценивание свойств почвы с использованием волновых сигналов сейсмических поверхностных волн |
US8537638B2 (en) * | 2010-02-10 | 2013-09-17 | Exxonmobil Upstream Research Company | Methods for subsurface parameter estimation in full wavefield inversion and reverse-time migration |
US8223587B2 (en) | 2010-03-29 | 2012-07-17 | Exxonmobil Upstream Research Company | Full wavefield inversion using time varying filters |
US8694299B2 (en) | 2010-05-07 | 2014-04-08 | Exxonmobil Upstream Research Company | Artifact reduction in iterative inversion of geophysical data |
US8756042B2 (en) | 2010-05-19 | 2014-06-17 | Exxonmobile Upstream Research Company | Method and system for checkpointing during simulations |
US8767508B2 (en) | 2010-08-18 | 2014-07-01 | Exxonmobil Upstream Research Company | Using seismic P and S arrivals to determine shallow velocity structure |
KR101219746B1 (ko) * | 2010-08-24 | 2013-01-10 | 서울대학교산학협력단 | 탄성 매질에서의 주파수 영역 역시간 구조보정을 이용한 지하구조의 영상화 장치 및 방법 |
CN103119552B (zh) | 2010-09-27 | 2016-06-08 | 埃克森美孚上游研究公司 | 同时源编码和源分离作为全波场反演的实际解决方案 |
US8437998B2 (en) | 2010-09-27 | 2013-05-07 | Exxonmobil Upstream Research Company | Hybrid method for full waveform inversion using simultaneous and sequential source method |
WO2012074592A1 (fr) | 2010-12-01 | 2012-06-07 | Exxonmobil Upstream Research Company | Inversion de source simultanée pour données de flûte sismique marines présentant une fonction objective de corrélation croisée |
EP2691795A4 (fr) | 2011-03-30 | 2015-12-09 | Vitesse de convergence d'une inversion d'un champ d'onde complet utilisant une mise en forme spectrale | |
CN103460074B (zh) | 2011-03-31 | 2016-09-28 | 埃克森美孚上游研究公司 | 全波场反演中小波估计和多次波预测的方法 |
ES2640824T3 (es) | 2011-09-02 | 2017-11-06 | Exxonmobil Upstream Research Company | Utilización de la proyección sobre conjuntos convexos para limitar la inversión del campo de onda completa |
US9176930B2 (en) | 2011-11-29 | 2015-11-03 | Exxonmobil Upstream Research Company | Methods for approximating hessian times vector operation in full wavefield inversion |
KR102021752B1 (ko) | 2012-03-08 | 2019-09-17 | 엑손모빌 업스트림 리서치 캄파니 | 직교 소스 및 수신기 인코딩 |
RU2598907C1 (ru) | 2012-11-28 | 2016-10-10 | Эксонмобил Апстрим Рисерч Компани | Сейсмическая отражательная q-томография |
MY169125A (en) | 2013-05-24 | 2019-02-18 | Exxonmobil Upstream Res Co | Multi-parameter inversion through offset dependent elastic fwi |
US10459117B2 (en) | 2013-06-03 | 2019-10-29 | Exxonmobil Upstream Research Company | Extended subspace method for cross-talk mitigation in multi-parameter inversion |
US9702998B2 (en) | 2013-07-08 | 2017-07-11 | Exxonmobil Upstream Research Company | Full-wavefield inversion of primaries and multiples in marine environment |
EP3351972A1 (fr) | 2013-08-23 | 2018-07-25 | Exxonmobil Upstream Research Company | Inversion iterative des donnees sismiques codees a base de la construction des donnees pseudo supersource |
US10036818B2 (en) | 2013-09-06 | 2018-07-31 | Exxonmobil Upstream Research Company | Accelerating full wavefield inversion with nonstationary point-spread functions |
US9910189B2 (en) | 2014-04-09 | 2018-03-06 | Exxonmobil Upstream Research Company | Method for fast line search in frequency domain FWI |
SG11201608175SA (en) | 2014-05-09 | 2016-11-29 | Exxonmobil Upstream Res Co | Efficient line search methods for multi-parameter full wavefield inversion |
US10185046B2 (en) | 2014-06-09 | 2019-01-22 | Exxonmobil Upstream Research Company | Method for temporal dispersion correction for seismic simulation, RTM and FWI |
CA2947410A1 (fr) | 2014-06-17 | 2015-12-30 | Exxonmobil Upstream Research Company | Inversion de champ d'ondes complet viscoacoustique et viscoelastique rapide |
US10838092B2 (en) | 2014-07-24 | 2020-11-17 | Exxonmobil Upstream Research Company | Estimating multiple subsurface parameters by cascaded inversion of wavefield components |
US10422899B2 (en) | 2014-07-30 | 2019-09-24 | Exxonmobil Upstream Research Company | Harmonic encoding for FWI |
US10386511B2 (en) | 2014-10-03 | 2019-08-20 | Exxonmobil Upstream Research Company | Seismic survey design using full wavefield inversion |
WO2016064462A1 (fr) | 2014-10-20 | 2016-04-28 | Exxonmobil Upstream Research Company | Tomographie de vitesse utilisant des balayages de propriété |
EP3234659A1 (fr) | 2014-12-18 | 2017-10-25 | Exxonmobil Upstream Research Company | Planification évolutive de travaux sismiques itératifs parallèles |
US10520618B2 (en) | 2015-02-04 | 2019-12-31 | ExxohnMobil Upstream Research Company | Poynting vector minimal reflection boundary conditions |
WO2016130208A1 (fr) | 2015-02-13 | 2016-08-18 | Exxonmobil Upstream Research Company | Condition de limite d'absorption efficace et stable dans des calculs à différence finie |
EP3259620B1 (fr) | 2015-02-17 | 2020-03-11 | Exxonmobil Upstream Research Company | Procédé d'inversion de champ d'ondes complet à plusieurs étapes génèrant un ensemble de données sans multiples |
US10416327B2 (en) | 2015-06-04 | 2019-09-17 | Exxonmobil Upstream Research Company | Method for generating multiple free seismic images |
US10838093B2 (en) | 2015-07-02 | 2020-11-17 | Exxonmobil Upstream Research Company | Krylov-space-based quasi-newton preconditioner for full-wavefield inversion |
KR102020759B1 (ko) | 2015-10-02 | 2019-09-11 | 엑손모빌 업스트림 리서치 캄파니 | Q-보상된 전 파동장 반전 |
EP3362823B1 (fr) | 2015-10-15 | 2019-10-09 | ExxonMobil Upstream Research Company | Modele d'inversion complet d'onde a preservation d'amplitude par domaine de stacks d'angle |
US10768324B2 (en) | 2016-05-19 | 2020-09-08 | Exxonmobil Upstream Research Company | Method to predict pore pressure and seal integrity using full wavefield inversion |
US11320557B2 (en) * | 2020-03-30 | 2022-05-03 | Saudi Arabian Oil Company | Post-stack time domain image with broadened spectrum |
US20220413176A1 (en) * | 2021-06-28 | 2022-12-29 | Halliburton Energy Services, Inc. | Annulus Velocity Independent Time Domain Structure Imaging In Cased Holes Using Multi-Offset Secondary Flexural Wave Data |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4870580A (en) * | 1983-12-30 | 1989-09-26 | Schlumberger Technology Corporation | Compressional/shear wave separation in vertical seismic profiling |
US4797861A (en) * | 1985-11-18 | 1989-01-10 | Western Atlas International, Inc. | Method of processing seismic data |
US5083297A (en) * | 1990-06-26 | 1992-01-21 | Chevron Research And Technology Company | Method of improving the seismic resolution of geologic structures |
US5596547A (en) * | 1995-10-20 | 1997-01-21 | Bancroft; John C. | Prestack seismic migration |
US6049759A (en) * | 1998-01-16 | 2000-04-11 | Bp Amoco Corporation | Method of prestack 3-D migration |
US6058073A (en) * | 1999-03-30 | 2000-05-02 | Atlantic Richfield Company | Elastic impedance estimation for inversion of far offset seismic sections |
-
1998
- 1998-10-01 FR FR9812337A patent/FR2784195B1/fr not_active Expired - Fee Related
-
1999
- 1999-09-29 GB GB9923053A patent/GB2342165B/en not_active Expired - Fee Related
- 1999-09-29 CA CA002284386A patent/CA2284386C/fr not_active Expired - Fee Related
- 1999-09-30 US US09/408,515 patent/US6311133B1/en not_active Expired - Lifetime
- 1999-09-30 BR BRPI9904428-5A patent/BR9904428B1/pt not_active IP Right Cessation
- 1999-09-30 NO NO19994776A patent/NO331526B1/no not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
GB9923053D0 (en) | 1999-12-01 |
US6311133B1 (en) | 2001-10-30 |
GB2342165B (en) | 2002-11-06 |
FR2784195B1 (fr) | 2000-11-17 |
NO331526B1 (no) | 2012-01-23 |
BR9904428B1 (pt) | 2011-06-28 |
NO994776D0 (no) | 1999-09-30 |
CA2284386A1 (fr) | 2000-04-01 |
FR2784195A1 (fr) | 2000-04-07 |
BR9904428A (pt) | 2000-09-05 |
GB2342165A (en) | 2000-04-05 |
NO994776L (no) | 2000-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2284386C (fr) | Methode pour realiser en 3d avant sommation, une migration de donnees sismiques | |
EP0789849B1 (fr) | Methode pour representer dans un domaine spatio-temporel la trajectoire d'au moins un puits de forage | |
FR2979711A1 (fr) | Procede de traitement de donnees sismiques marines | |
EP1239304A1 (fr) | Perfectionnements aux procédés d'inversion tomographique d'événements pointés sur les données sismiques migrées | |
EP0313459B1 (fr) | Procédé d'acquisition de données sismiques et dispositif de mise en oeuvre du procédé | |
FR2955396A1 (fr) | Dispositif de traitement de donnees sismiques marines | |
FR2876458A1 (fr) | Perfectionnement aux traitements sismiques pour la suppression des reflexions multiples | |
FR2782389A1 (fr) | Procede et dispositif pour corriger les effets des deplacements du navire dans les mesures de sismique maritime | |
FR2800473A1 (fr) | Methode pour modeliser en 2d ou 3d un milieu heterogene tel que le sous-sol decrit par plusieurs parametres physiques | |
FR2478324A1 (fr) | Procede de prospection et d'exploration sismiques | |
FR2985040A1 (fr) | Device and method for deghosting variable depth streamer data | |
FR2916859A1 (fr) | Procede de traitement de donnees sismiques | |
CA2528700A1 (fr) | Methode pour determiner des informations speculaires apres imagerie sismique avant sommation | |
FR2874270A1 (fr) | Procede d'exploration sismique | |
FR2753800A1 (fr) | Procede et appareil pour exploiter des paires de donnees collectees au cours d'une operation sismique, lesquelles representent des caracteristiques d'une formation terrestre et pour generer un enregistrement de sortie | |
CA2364938C (fr) | Methode pour obtenir des temps de trajet en reflexion a partir d'une interpretation de donnees sismiques en ondes cylindriques migrees | |
FR2981170A1 (fr) | Procede de tomographie non lineaire pour un axe de symetrie principal d'un modele de vitesse anisotrope et dispositif. | |
EP0572321B1 (fr) | Procédé de traitement de signaux pour prospection géophysique exploitant un opérateur d'extrapolation d'un champ d'onde perfectionné | |
FR2971859A1 (fr) | Procede et dispositif de lissage a preservation du temps de trajet | |
FR2799552A1 (fr) | Traitement de mesures de formes d'ondes acoustiques provenant d'outils de diagraphie de sondage en matrice | |
EP0410826B1 (fr) | Procédé itératif d'estimation de mouvement, entre une image de référence et une image courante, et dispositif pour la mise en oeuvre de ce procédé | |
FR2985818A1 (fr) | Dispositif et procede d'estimation de decalages temporels | |
CA2274922C (fr) | Procede de traitement sismique et notamment procede de prospection sismique 3d mettant en oeuvre une migration des donnees sismiques | |
FR3132960A1 (fr) | Localisation d’une source acoustique en mouvement | |
Egan et al. | The influence of spatial sampling on resolution |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EEER | Examination request | ||
MKLA | Lapsed |
Effective date: 20131001 |