CA2464799C - Methode pour determiner un modele de vitesse d'ondes sismiques dans une formation souterraine heterogene - Google Patents
Methode pour determiner un modele de vitesse d'ondes sismiques dans une formation souterraine heterogene Download PDFInfo
- Publication number
- CA2464799C CA2464799C CA2464799A CA2464799A CA2464799C CA 2464799 C CA2464799 C CA 2464799C CA 2464799 A CA2464799 A CA 2464799A CA 2464799 A CA2464799 A CA 2464799A CA 2464799 C CA2464799 C CA 2464799C
- Authority
- CA
- Canada
- Prior art keywords
- seismic
- speed
- offset
- recordings
- interface
- 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 69
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 16
- 238000003325 tomography Methods 0.000 claims abstract description 8
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 claims description 7
- 238000013508 migration Methods 0.000 claims description 4
- 230000005012 migration Effects 0.000 claims description 4
- 230000002123 temporal effect Effects 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 2
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims 1
- 238000003384 imaging method Methods 0.000 abstract description 4
- 238000012549 training Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 241000276498 Pollachius virens Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling 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/303—Analysis for determining velocity profiles or travel times
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 pour déterminer un modèle de vitesse d'ondes sismiques captées par des récepteurs sismiques couplés avec une formation souterraine, en réponse à l'émission d'ondes sismiques dans le sous-sol par une source sismique, après réflexion sur des interfaces géologiques de la dite formation, d'après des enregistrements multi- déports de ces ondes. - A partir du pointé sur les enregistrements sismiques d'un ou plusieurs événements pour différents couples source-capteur associés au déport le plus faible d'un dispositif sismique d'acquisition, on prédit la géométrie des événements considérés pour toute la gamme de déports du dispositif d'acquisition et pour différentes valeurs de vitesse prises autour d'une vitesse de couche de référence, et ce à l'aide de techniques cinématiques d'inversion et de modélisation. Puis, pour différentes positions latérales du dispositif d'acquisition, on sélectionne la courbe, de temps de trajet prédite approchant au mieux l'événement sismique considéré. Les temps de trajet ainsi obtenus sont alors fournis à une méthode d'inversion avant sommation telle qu'une tomographie de temps de trajet. - Applications â l'imagerie sismique du sous-sol par exemple.
Description
METHODE POUR DETERMINER UN MODELE DE VITESSE D'ONDES
SISMIQUES DANS UNE FORMATION SOUTERRAINE HETEROGENE
La présente invention concerne une méthode permettant de déterminer un modèle de vitesse d'ondes sismiques captées par des récepteurs couplés avec une formation souterraine, d'après des enregistrements multi-déports de ces ondes.
La méthode permet un accès facilité à l'information cinématique avant sommation associée aux événements contenus dans les enregistrements sismiques, qui est combiné à
une technique d'inversion cinématique avant sommation en vue de l'imagerie des interfaces géologiques du sous-sol.
Par "information cinématique avant sommation", nous entendons les temps de trajet associés aux réflexions enregistrées par des couples source-capteur situés à
des distances variables l'un de l'autre. A "accès facilité", nous opposons un pointé manuel des événements sismiques dans les collections multi-déports. Par "méthode d'inversion cinématique avant sommation", nous entendons une méthode permettant, à partir de l'information cinématique extraite des enregistrements sismiques (et non d'une approximation de cette information cinématique), de retrouver la géométrie des réflecteurs associés ainsi que les vitesses de couche.
ETAT DE LA TECHNIQUE
La sismique réflexion est largement utilisée en exploration pétrolière, notamment pour produire des images de la subsurface à partir de l'information contenue dans les ondes qui se sont propagées et réfléchies sur les discontinuités géologiques du sous-sol.
SISMIQUES DANS UNE FORMATION SOUTERRAINE HETEROGENE
La présente invention concerne une méthode permettant de déterminer un modèle de vitesse d'ondes sismiques captées par des récepteurs couplés avec une formation souterraine, d'après des enregistrements multi-déports de ces ondes.
La méthode permet un accès facilité à l'information cinématique avant sommation associée aux événements contenus dans les enregistrements sismiques, qui est combiné à
une technique d'inversion cinématique avant sommation en vue de l'imagerie des interfaces géologiques du sous-sol.
Par "information cinématique avant sommation", nous entendons les temps de trajet associés aux réflexions enregistrées par des couples source-capteur situés à
des distances variables l'un de l'autre. A "accès facilité", nous opposons un pointé manuel des événements sismiques dans les collections multi-déports. Par "méthode d'inversion cinématique avant sommation", nous entendons une méthode permettant, à partir de l'information cinématique extraite des enregistrements sismiques (et non d'une approximation de cette information cinématique), de retrouver la géométrie des réflecteurs associés ainsi que les vitesses de couche.
ETAT DE LA TECHNIQUE
La sismique réflexion est largement utilisée en exploration pétrolière, notamment pour produire des images de la subsurface à partir de l'information contenue dans les ondes qui se sont propagées et réfléchies sur les discontinuités géologiques du sous-sol.
2 Plus précisément, les méthodes d'imagerie utilisent l'information cinématique associée aux réflexions sismiques majeures (i.e. les temps de trajet des ondes s'étant réfléchies sur les discontinuités principales du sous-sol) pour déterminer un macro modèle de vitesse du sous-sol, qui servira lui-même à transformer les enregistrements sismiques temporels en une image profondeur de la subsurface.
L' accès à l'information cinématique nécessaire à la détermination du modèle de vitesse requiert l'interprétation des événements sismiques dans les enregistrements sismiques multi-déports. Or les données sismiques multi-déports sont généralement caractérisées par un mauvais rapport signal sur bruit, d'où une défaillance des pointés automatiques classiques et la nécessité d'effectuer un pointé des données sismiques manuel long et coûteux. Par ailleurs, pour les campagnes sismiques 3D (prédominantes de nos jours par rapport aux campagnes 2D), la quantité de données à interpréter est un volume 4D qui peut atteindre des tailles de l'ordre de la centaine de GigaOctets, voire de l'ordre du TeraOctets, rendant la tâche interprétative encore plus longue et fastidieuse.
Afin d'éviter cette étape d'interprétation des collections sismiques multi-déports, les géophysiciens ont développé des méthodes basées sur des approximations de la géométrie des événements sismiques dans les collections multi-déports. Pour établir ces approximations, ces méthodes émettent des hypothèses sur la complexité du sous-sol, hypothèses qui peuvent porter sur la géométrie des discontinuités géologiques du sous-sol et sur les variations de vitesse de couches. On citera par exemple la méthode décrite dans Taner et Kohler (1969) qui suppose un milieu stratifié composé de couches planes horizontales à vitesses de couches homogènes et isotropes, celle décrite dans Levin (1971) qui étend la méthode précédente à des couches pentées. D'autres variantes ont été
proposées, mais de façon générale, les méthodes existantes sont basées sur une hypothèse hyperbolique de la courbe des temps de trajet dans les collections multi-déports. Or cette hypothèse est violée dès que les discontinuités géologiques du sous-sol ne sont plus planes et/ou dès que les vitesses de couche présentent des variations latérales.
Nous proposons ici une méthode pour approcher au mieux les temps de trajet des collections multi-déports dans le cas où les discontinuités géologiques du sous-sol ont une géométrie quelconque et où les vitesses de couche sont modérément variables latéralement, et ce, tout en limitant le temps humain requis pour l'interprétation des enregistrements sismiques. Les temps de trajet obtenus seront alors traités par une méthode d'inversion
L' accès à l'information cinématique nécessaire à la détermination du modèle de vitesse requiert l'interprétation des événements sismiques dans les enregistrements sismiques multi-déports. Or les données sismiques multi-déports sont généralement caractérisées par un mauvais rapport signal sur bruit, d'où une défaillance des pointés automatiques classiques et la nécessité d'effectuer un pointé des données sismiques manuel long et coûteux. Par ailleurs, pour les campagnes sismiques 3D (prédominantes de nos jours par rapport aux campagnes 2D), la quantité de données à interpréter est un volume 4D qui peut atteindre des tailles de l'ordre de la centaine de GigaOctets, voire de l'ordre du TeraOctets, rendant la tâche interprétative encore plus longue et fastidieuse.
Afin d'éviter cette étape d'interprétation des collections sismiques multi-déports, les géophysiciens ont développé des méthodes basées sur des approximations de la géométrie des événements sismiques dans les collections multi-déports. Pour établir ces approximations, ces méthodes émettent des hypothèses sur la complexité du sous-sol, hypothèses qui peuvent porter sur la géométrie des discontinuités géologiques du sous-sol et sur les variations de vitesse de couches. On citera par exemple la méthode décrite dans Taner et Kohler (1969) qui suppose un milieu stratifié composé de couches planes horizontales à vitesses de couches homogènes et isotropes, celle décrite dans Levin (1971) qui étend la méthode précédente à des couches pentées. D'autres variantes ont été
proposées, mais de façon générale, les méthodes existantes sont basées sur une hypothèse hyperbolique de la courbe des temps de trajet dans les collections multi-déports. Or cette hypothèse est violée dès que les discontinuités géologiques du sous-sol ne sont plus planes et/ou dès que les vitesses de couche présentent des variations latérales.
Nous proposons ici une méthode pour approcher au mieux les temps de trajet des collections multi-déports dans le cas où les discontinuités géologiques du sous-sol ont une géométrie quelconque et où les vitesses de couche sont modérément variables latéralement, et ce, tout en limitant le temps humain requis pour l'interprétation des enregistrements sismiques. Les temps de trajet obtenus seront alors traités par une méthode d'inversion
3 cinématique avant sommation qui permet de prendre en compte des cinématiques complexes telle que la tomographie de temps de trajet, décrite par exemple dans le document suivant :
- Bishop, T. et ai, 1985, Tomographic Détermination of Velocity and Depth in Laterally Varying Media Geophysics, 50 N 6, 903-923.
La méthode pourra également s'appliquer dans les cas où les vitesses de couche sont variables latéralement, mais de façon itérative, le modèle de vitesse après mise à jour à une itération donnée par application d'une méthode d'inversion cinématique avant sommation telle que la tomographie de réflexion précédemment citée servant de modèle initial pour l'itération suivante.
LA METHODE SELON L'INVENTION
L'invention vise une méthode pour transformer des enregistrements sismiques temporels en une image profondeur d'une formation souterraine, lesdits enregistrements sismiques correspondant à des enregistrements multi-déports d'ondes sismiques captées par des récepteurs sismiques couplés avec la formation souterraine, en réponse à l'émission d'ondes sismiques dans ladite formation par une source sismique, après réflexion sur des interfaces géologiques de la dite formation, dans laquelle on détermine un modèle de vitesse desdites ondes sismiques, caractérisée en ce qu'elle comporte au moins les étapes suivantes pour chaque événement sismique repéré sur les enregistrements et pour chaque couche délimitée par les dites interfaces:
a) à partir des enregistrements sismiques avant sommation, on construit une collection iso-déport, dont on extrait une information cinématique ou des temps de trajet associés à l'événement;
b) on choisit une gamme de vitesse autour d'une vitesse de référence dans ladite couche que l'on échantillonne avec un pas défini;
3a c) pour chaque échantillon de vitesse, on applique une technique d'inversion à
vitesse fixée afin de déterminer, à partir des temps de trajet extraits de la collection iso-déport, la géométrie de la dite interface pour l'échantillon de vitesse en question de façon à obtenir une série de couples interface/vitesse pour le dit événement;
d) on calcule l'information cinématique associée à chaque couple interface/vitesse obtenu, pour des couples source-récepteur correspondants à
des collections multi-déports existant dans les enregistrements sismiques;
e) pour chaque couple interface/vitesse et pour chaque collection multi-déports sélectionnée, on évalue la cohérence entre les temps de trajet multi-déports ainsi calculés et les enregistrements sismiques, et on sélectionne pour chaque collection multi-déports la courbe de temps de trajet offrant la cohérence maximale avec les enregistrements sismiques;
f) on applique une méthode d'inversion cinématique avant sommation en utilisant les temps de trajet multi-déports obtenus pour l'ensemble des collections multi-déports sélectionnées, afin de déterminer la géométrie et la vitesse de la couche considérée;
g) on itère n fois (n >_ 0) les étapes a) à f) en considérant à chaque itération le modèle de vitesse obtenu à l'itération précédente comme modèle de référence pour définir la dite vitesse de référence; et h) on utilise ledit modèle de vitesse obtenu à l'issue de l'étape g), pour construire ladite image profondeur de la formation souterraine.
De préférence, la méthode selon l'invention permet de déterminer un modèle de vitesse d'ondes sismiques captées par des récepteurs sismiques couplés avec une formation souterraine, en réponse à l'émission d'ondes sismiques dans le sous-sol par une source sismique, après réflexion sur des interfaces géologiques de ladite formation, à partir d'enregistrements multi-déports de ces ondes.
De préférence, pour chaque événement sismique repéré sur les enregistrements et pour chaque couche délimitée par lesdites interfaces, la méthode comporte au moins les étapes suivantes:
3b a) à partir des enregistrements sismiques avant sommation, on construit une collection iso-déport, dont on extrait une information cinématique ou des temps de trajet associés à l'événement ;
b) on choisit une gamme de vitesse autour d'une vitesse de référence dans la dite couche que l'on échantillonne avec un pas défini ;
c) pour chaque échantillon de vitesse, on applique une technique d'inversion à
vitesse fixée afin de déterminer, à partir des temps de trajet extraits de la collection iso-déport, la géométrie de la dite interface pour l'échantillon de vitesse en question de façon à
obtenir une série de couples interface/ vitesse pour le dit événement ;
- Bishop, T. et ai, 1985, Tomographic Détermination of Velocity and Depth in Laterally Varying Media Geophysics, 50 N 6, 903-923.
La méthode pourra également s'appliquer dans les cas où les vitesses de couche sont variables latéralement, mais de façon itérative, le modèle de vitesse après mise à jour à une itération donnée par application d'une méthode d'inversion cinématique avant sommation telle que la tomographie de réflexion précédemment citée servant de modèle initial pour l'itération suivante.
LA METHODE SELON L'INVENTION
L'invention vise une méthode pour transformer des enregistrements sismiques temporels en une image profondeur d'une formation souterraine, lesdits enregistrements sismiques correspondant à des enregistrements multi-déports d'ondes sismiques captées par des récepteurs sismiques couplés avec la formation souterraine, en réponse à l'émission d'ondes sismiques dans ladite formation par une source sismique, après réflexion sur des interfaces géologiques de la dite formation, dans laquelle on détermine un modèle de vitesse desdites ondes sismiques, caractérisée en ce qu'elle comporte au moins les étapes suivantes pour chaque événement sismique repéré sur les enregistrements et pour chaque couche délimitée par les dites interfaces:
a) à partir des enregistrements sismiques avant sommation, on construit une collection iso-déport, dont on extrait une information cinématique ou des temps de trajet associés à l'événement;
b) on choisit une gamme de vitesse autour d'une vitesse de référence dans ladite couche que l'on échantillonne avec un pas défini;
3a c) pour chaque échantillon de vitesse, on applique une technique d'inversion à
vitesse fixée afin de déterminer, à partir des temps de trajet extraits de la collection iso-déport, la géométrie de la dite interface pour l'échantillon de vitesse en question de façon à obtenir une série de couples interface/vitesse pour le dit événement;
d) on calcule l'information cinématique associée à chaque couple interface/vitesse obtenu, pour des couples source-récepteur correspondants à
des collections multi-déports existant dans les enregistrements sismiques;
e) pour chaque couple interface/vitesse et pour chaque collection multi-déports sélectionnée, on évalue la cohérence entre les temps de trajet multi-déports ainsi calculés et les enregistrements sismiques, et on sélectionne pour chaque collection multi-déports la courbe de temps de trajet offrant la cohérence maximale avec les enregistrements sismiques;
f) on applique une méthode d'inversion cinématique avant sommation en utilisant les temps de trajet multi-déports obtenus pour l'ensemble des collections multi-déports sélectionnées, afin de déterminer la géométrie et la vitesse de la couche considérée;
g) on itère n fois (n >_ 0) les étapes a) à f) en considérant à chaque itération le modèle de vitesse obtenu à l'itération précédente comme modèle de référence pour définir la dite vitesse de référence; et h) on utilise ledit modèle de vitesse obtenu à l'issue de l'étape g), pour construire ladite image profondeur de la formation souterraine.
De préférence, la méthode selon l'invention permet de déterminer un modèle de vitesse d'ondes sismiques captées par des récepteurs sismiques couplés avec une formation souterraine, en réponse à l'émission d'ondes sismiques dans le sous-sol par une source sismique, après réflexion sur des interfaces géologiques de ladite formation, à partir d'enregistrements multi-déports de ces ondes.
De préférence, pour chaque événement sismique repéré sur les enregistrements et pour chaque couche délimitée par lesdites interfaces, la méthode comporte au moins les étapes suivantes:
3b a) à partir des enregistrements sismiques avant sommation, on construit une collection iso-déport, dont on extrait une information cinématique ou des temps de trajet associés à l'événement ;
b) on choisit une gamme de vitesse autour d'une vitesse de référence dans la dite couche que l'on échantillonne avec un pas défini ;
c) pour chaque échantillon de vitesse, on applique une technique d'inversion à
vitesse fixée afin de déterminer, à partir des temps de trajet extraits de la collection iso-déport, la géométrie de la dite interface pour l'échantillon de vitesse en question de façon à
obtenir une série de couples interface/ vitesse pour le dit événement ;
4 d) on calcule l'information cinématique associée à chaque couple interface/vitesse obtenu, pour des couples source-récepteur correspondants à des collections multi-déports existant dans les enregistrements sismiques ;
e) pour chaque couple interface/vitesse et pour chaque collection multi-déports sélectionnée, on évalue la cohérence entre les temps de trajet multi-déports ainsi calculés et les enregistrements sismiques, et on sélectionne pour chaque collection multi-déports la courbe de temps de trajet offrant la cohérence maximale avec les enregistrements sismiques ;
f) on applique une méthode d'inversion cinématique avant sommation en utilisant les temps de trajet multi-déports obtenus pour l'ensemble des collections multi-déports sélectionnées, afin de déterminer la géométrie et la vitesse de la couche considérée ; et g) on itère n fois (n >_ 0) les étapes a) à f) en considérant à chaque itération le modèle de vitesse obtenu à l'itération précédente comme modèle de référence pour définir la vitesse de référence de la nouvelle itération.
Suivant un mode de mise en oeuvre convenant au cas où la gamme de vitesse choisie ne s'avère pas assez précise à l'issue soit de l'itération précédente soit des étapes a) à f), on réalise n fois l'étape g) avec n >_ 1.
Suivant un autre mode de mise en ouvre convenant au cas où la distribution des vitesses varie fortement latéralement et/ou dans la cas où l'on ne dispose pas d'une connaissance a priori suffisamment précise sur la distribution de vitesse dans la couche considérée, on réalise l'étape g) sur des gammes de déports de plus en plus larges au fur et à mesure des itérations.
Suivant un autre mode de mise en ouvre convenant au cas où la distribution des vitesses varie fortement latéralement et/ou dans le cas où l'on ne dispose pas d'une connaissance a priori suffisamment précise sur la distribution de vitesse dans la couche considérée, on réalise l'étape g) sur des grilles de collections multi-déports de plus en plus fines au fur et à mesure des itérations.
Suivant un mode de mise en oeuvre convenant au cas où la géométrie de l'interface recherchée est génératrice de triplications, on réalise l'étape e) en considérant les branches de triplications dans les collections multi-déports indépendamment les unes des autres. On utilise par exemple des outils de tracé de rayons et d'inversion permettant de prendre en compte des arrivées multi-valuées.
Suivant un exemple de mise en oeuvre, à l'étape a) on construit de préférence une collection iso-déport à déport nul ou à déport proche.
e) pour chaque couple interface/vitesse et pour chaque collection multi-déports sélectionnée, on évalue la cohérence entre les temps de trajet multi-déports ainsi calculés et les enregistrements sismiques, et on sélectionne pour chaque collection multi-déports la courbe de temps de trajet offrant la cohérence maximale avec les enregistrements sismiques ;
f) on applique une méthode d'inversion cinématique avant sommation en utilisant les temps de trajet multi-déports obtenus pour l'ensemble des collections multi-déports sélectionnées, afin de déterminer la géométrie et la vitesse de la couche considérée ; et g) on itère n fois (n >_ 0) les étapes a) à f) en considérant à chaque itération le modèle de vitesse obtenu à l'itération précédente comme modèle de référence pour définir la vitesse de référence de la nouvelle itération.
Suivant un mode de mise en oeuvre convenant au cas où la gamme de vitesse choisie ne s'avère pas assez précise à l'issue soit de l'itération précédente soit des étapes a) à f), on réalise n fois l'étape g) avec n >_ 1.
Suivant un autre mode de mise en ouvre convenant au cas où la distribution des vitesses varie fortement latéralement et/ou dans la cas où l'on ne dispose pas d'une connaissance a priori suffisamment précise sur la distribution de vitesse dans la couche considérée, on réalise l'étape g) sur des gammes de déports de plus en plus larges au fur et à mesure des itérations.
Suivant un autre mode de mise en ouvre convenant au cas où la distribution des vitesses varie fortement latéralement et/ou dans le cas où l'on ne dispose pas d'une connaissance a priori suffisamment précise sur la distribution de vitesse dans la couche considérée, on réalise l'étape g) sur des grilles de collections multi-déports de plus en plus fines au fur et à mesure des itérations.
Suivant un mode de mise en oeuvre convenant au cas où la géométrie de l'interface recherchée est génératrice de triplications, on réalise l'étape e) en considérant les branches de triplications dans les collections multi-déports indépendamment les unes des autres. On utilise par exemple des outils de tracé de rayons et d'inversion permettant de prendre en compte des arrivées multi-valuées.
Suivant un exemple de mise en oeuvre, à l'étape a) on construit de préférence une collection iso-déport à déport nul ou à déport proche.
5 Suivant un autre exemple de mise en oeuvre, à l'étape c) on applique avantageusement une technique d'inversion cinématique à vitesse fixée telle qu'une migration de carte.
Suivant un autre exemple de mise en oeuvre, on calcule à l'étape d) l'information cinématique en traçant des rayons multi-déports sur l'interface de chaque couple interface-vitesse.
Suivant un autre exemple de mise en oeuvre, on réalise avantageusement l'étape f) en appliquant une méthode d'inversion cinématique avant sommation telle qu'une tomographie de temps de trajet avant sommation.
Dans sa définition la plus générale correspondant au cas où le nombre n d'itérations à réaliser peut éventuellement être nul, la méthode s'applique par exemple si la distribution de vitesse dans chaque couche est peu variable latéralement et/ou dans le cas ou l'on dispose d'un a priori sur la distribution de vitesse dans la couche considérée suffisamment précis.
Dans sa définition plus restrictive où le nombre n d'itérations à réaliser est au moins égal à 1, la méthode s'applique par exemple dans le cas où la distribution des vitesses dans les couches de la formation géologique est fortement variable latéralement et/ou dans le cas où l'on ne dispose pas d'un a priori suffisamment précis sur la distribution de vitesse dans la couche considérée.
PRESENTATION SUCCINCTE DES FIGURES
Les caractéristiques et avantages de la méthode selon l'invention, apparaîtront plus clairement à la lecture de la description ci-après d'exemples non limitatifs de mise en oeuvre, en se référant aux dessins annexés où :
Suivant un autre exemple de mise en oeuvre, on calcule à l'étape d) l'information cinématique en traçant des rayons multi-déports sur l'interface de chaque couple interface-vitesse.
Suivant un autre exemple de mise en oeuvre, on réalise avantageusement l'étape f) en appliquant une méthode d'inversion cinématique avant sommation telle qu'une tomographie de temps de trajet avant sommation.
Dans sa définition la plus générale correspondant au cas où le nombre n d'itérations à réaliser peut éventuellement être nul, la méthode s'applique par exemple si la distribution de vitesse dans chaque couche est peu variable latéralement et/ou dans le cas ou l'on dispose d'un a priori sur la distribution de vitesse dans la couche considérée suffisamment précis.
Dans sa définition plus restrictive où le nombre n d'itérations à réaliser est au moins égal à 1, la méthode s'applique par exemple dans le cas où la distribution des vitesses dans les couches de la formation géologique est fortement variable latéralement et/ou dans le cas où l'on ne dispose pas d'un a priori suffisamment précis sur la distribution de vitesse dans la couche considérée.
PRESENTATION SUCCINCTE DES FIGURES
Les caractéristiques et avantages de la méthode selon l'invention, apparaîtront plus clairement à la lecture de la description ci-après d'exemples non limitatifs de mise en oeuvre, en se référant aux dessins annexés où :
6 les figures la à le montrent respectivement une application de la méthode à
une collection multi-déports de type récepteur commun (Fig.la), à une collection récepteur commun sur laquelle ont été superposés les temps de trajet multi-déports obtenus avec différentes valeurs de vitesse de couche (Fig.lb) et à une collection récepteur commun sur laquelle a été superposée la courbe de temps de trajet multi-déports prédite montrant le plus de cohérence avec l'événement sismique considéré (Fig.1c).
DESCRIPTION DETAILLEE
On dispose d'enregistrements sismiques que l'on a obtenus au moyen d'un dispositif sismique comportant une source sismique émettant des ondes sismiques qui se propagent dans le sous-sol, un ensemble de récepteurs sismiques couplés avec le milieu qui captent les ondes renvoyées par les discontinuités du sous-sol en réponse aux ondes émises, et d'un laboratoire d'enregistrement des signaux sismiques captés.
1- Cas standard Nous présentons ci-après un premier mode de mise en oeuvre de la méthode, dans le cas où la distribution des vitesses dans la formation géologique précédant la couche en question étant connue ou ayant été estimée par ailleurs, la géométrie de la couche en question ne donnant pas lieu à des triplications, la vitesse à déterminer étant peu variable latéralement ou un a priori suffisamment précis existant sur cette vitesse. On réalise les étapes suivantes:
I-1 A partir des enregistrements sismiques avant sommation ou encore multi-déports, on construit une collection à déport constant, de préférence à déport nul (en réalité
une approximation d'une collection à déport nul, connue sous le nom de stack , caractérisée par un meilleur rapport signal sur bruit que les sections iso-déport individuelle) ou à déport proche (correspondant au déport le plus faible des enregistrements multi-déports).
Si une interprétation cohérente et suffisamment complète de l'événement sismique en cours est réalisable dans cette collection, on procède à son pointé
directement dans cette collection; sinon, afin d'accéder aux temps de trajet associés à l'événement en question, on pourra avoir recours à un détour dans le domaine migré temps ou profondeur suivi d'une démigration temps ou profondeur, comme selon la méthode décrite dans :
une collection multi-déports de type récepteur commun (Fig.la), à une collection récepteur commun sur laquelle ont été superposés les temps de trajet multi-déports obtenus avec différentes valeurs de vitesse de couche (Fig.lb) et à une collection récepteur commun sur laquelle a été superposée la courbe de temps de trajet multi-déports prédite montrant le plus de cohérence avec l'événement sismique considéré (Fig.1c).
DESCRIPTION DETAILLEE
On dispose d'enregistrements sismiques que l'on a obtenus au moyen d'un dispositif sismique comportant une source sismique émettant des ondes sismiques qui se propagent dans le sous-sol, un ensemble de récepteurs sismiques couplés avec le milieu qui captent les ondes renvoyées par les discontinuités du sous-sol en réponse aux ondes émises, et d'un laboratoire d'enregistrement des signaux sismiques captés.
1- Cas standard Nous présentons ci-après un premier mode de mise en oeuvre de la méthode, dans le cas où la distribution des vitesses dans la formation géologique précédant la couche en question étant connue ou ayant été estimée par ailleurs, la géométrie de la couche en question ne donnant pas lieu à des triplications, la vitesse à déterminer étant peu variable latéralement ou un a priori suffisamment précis existant sur cette vitesse. On réalise les étapes suivantes:
I-1 A partir des enregistrements sismiques avant sommation ou encore multi-déports, on construit une collection à déport constant, de préférence à déport nul (en réalité
une approximation d'une collection à déport nul, connue sous le nom de stack , caractérisée par un meilleur rapport signal sur bruit que les sections iso-déport individuelle) ou à déport proche (correspondant au déport le plus faible des enregistrements multi-déports).
Si une interprétation cohérente et suffisamment complète de l'événement sismique en cours est réalisable dans cette collection, on procède à son pointé
directement dans cette collection; sinon, afin d'accéder aux temps de trajet associés à l'événement en question, on pourra avoir recours à un détour dans le domaine migré temps ou profondeur suivi d'une démigration temps ou profondeur, comme selon la méthode décrite dans :
7 Ehinger, A., and Lailly, P., 1995, Velocity model determination by the SMART
method, Part 1: Theory: 65`h Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, pages 739-742.
I-2 On se donne une vitesse de référence pour la couche en question, ainsi qu'une incertitude sur cette vitesse de référence et un pas d'échantillonnage de la gamme de vitesse ainsi formée; on obtient ainsi une série d'échantillons de vitesse pour la couche en question ; la vitesse de référence peut être quelconque (i.e.
variable latéralement et verticalement); en particulier, si l'on dispose d'une distribution de vitesse obtenue par ailleurs, on pourra utiliser cette vitesse comme vitesse de référence;
si l'on ne dispose pas d'information a priori sur la vitesse de la couche, on pourra par exemple prendre une vitesse de référence homogène latéralement et verticalement.
I-3 Pour chacun des échantillons de vitesse, on applique une technique d'inversion de type cinématique à vitesse fixée telle qu'une migration de carte par exemple, afin de déterminer, à partir des temps de trajet extraits de la section à déport constant, la géométrie de la dite interface pour l'échantillon de vitesse en question ; on obtient ainsi une série de couples interface/ vitesse pour l'événement en question; on trouvera une description de la migration de carte par exemple dans la référence suivante :
- Yilmaz, 0., 2001, Seismic Data Analysis - Processing, Inversion and Interpretation of Seismic Data: Society of Exploration Geophysicists.
Il peut éventuellement s'avérer impossible de trouver un modèle d'interface expliquant les temps de trajet pour la vitesse fixée en cours (impossibilité
d'atteindre la convergence). Pour une plus grande efficacité de la méthode, on peut ne sélectionner que des couples interface/vitesse pour lesquels l'écart entre les temps de trajet observés et les temps de trajet calculés sur le modèle obtenu après convergence est inférieur à un certain seuil.
A noter que cette inversion à vitesse fixée peut tout aussi bien être réalisée par toute autre technique cinématique équivalente.
I-4 On calcule l'information cinématique associée à chaque couple interface/vitesse obtenu par exemple par un tracé de rayons multi-déports sur l'interface
method, Part 1: Theory: 65`h Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, pages 739-742.
I-2 On se donne une vitesse de référence pour la couche en question, ainsi qu'une incertitude sur cette vitesse de référence et un pas d'échantillonnage de la gamme de vitesse ainsi formée; on obtient ainsi une série d'échantillons de vitesse pour la couche en question ; la vitesse de référence peut être quelconque (i.e.
variable latéralement et verticalement); en particulier, si l'on dispose d'une distribution de vitesse obtenue par ailleurs, on pourra utiliser cette vitesse comme vitesse de référence;
si l'on ne dispose pas d'information a priori sur la vitesse de la couche, on pourra par exemple prendre une vitesse de référence homogène latéralement et verticalement.
I-3 Pour chacun des échantillons de vitesse, on applique une technique d'inversion de type cinématique à vitesse fixée telle qu'une migration de carte par exemple, afin de déterminer, à partir des temps de trajet extraits de la section à déport constant, la géométrie de la dite interface pour l'échantillon de vitesse en question ; on obtient ainsi une série de couples interface/ vitesse pour l'événement en question; on trouvera une description de la migration de carte par exemple dans la référence suivante :
- Yilmaz, 0., 2001, Seismic Data Analysis - Processing, Inversion and Interpretation of Seismic Data: Society of Exploration Geophysicists.
Il peut éventuellement s'avérer impossible de trouver un modèle d'interface expliquant les temps de trajet pour la vitesse fixée en cours (impossibilité
d'atteindre la convergence). Pour une plus grande efficacité de la méthode, on peut ne sélectionner que des couples interface/vitesse pour lesquels l'écart entre les temps de trajet observés et les temps de trajet calculés sur le modèle obtenu après convergence est inférieur à un certain seuil.
A noter que cette inversion à vitesse fixée peut tout aussi bien être réalisée par toute autre technique cinématique équivalente.
I-4 On calcule l'information cinématique associée à chaque couple interface/vitesse obtenu par exemple par un tracé de rayons multi-déports sur l'interface
8 de chaque couple interface/vitesse ainsi obtenu et pour des couples source-récepteur correspondant à des collections multi-déports existantes dans les enregistrements sismiques ; on choisit de préférence de travailler dans des collections point milieu commun, en particulier si la vitesse réelle présente des variations latérales non prises en compte dans le modèle de référence.
I-5 Pour chaque couple interface/vitesse et pour chaque collection multi-déports sélectionnée, on évalue la cohérence entre les temps de trajet multi-déports ainsi calculés et les enregistrements sismiques ; on peut par exemple employer une technique dite de calcul de semblante, pour laquelle on somme l'énergie rencontrée dans les données sismiques le long d'une courbe qui correspondra ici aux temps de trajet multi-déports calculés par tracé de rayons. Dans un cas idéal, il existera un échantillon de vitesse pour lequel la courbe de temps de trajet prédite se superposera parfaitement à
l'événement sismique considéré, ce qui se traduira par une mesure de cohérence maximale. Il est à noter que cette évaluation pourra se faire sur une gamme limitée de déport, définie par l'utilisateur ou automatiquement, par exemple par examen de la courbe montrant l'évolution de la cohérence en fonction de la gamme de déports.
Pour chaque collection multi-déports, on sélectionne la courbe de temps de trajet la plus en accord avec les enregistrements sismiques, et les temps de trajet multi-déports constituant cette courbe sont réservés. Il est à noter que la courbe de temps de trajet offrant la meilleure cohérence pourra être ajustée si nécessaire en cherchant dans une fenêtre verticale prédéfinie le maximum (ou le minimum en fonction de la polarité
de l'événement considéré) d'amplitude le plus proche pour chaque trace sismique de la collection multi-déport considérée.
I-6 Les temps de trajet multi-déports ainsi recueillis pour l'ensemble des collections multi-déports sélectionnées sont alors injectés dans une méthode d'inversion cinématique avant sommation telle qu'une tomographie de temps de trajet afin de déterminer la géométrie et la vitesse de la couche en question.
I-7 Il est possible que l'on soit obligé d'itérer les étapes précédentes dans le cas où la cohérence avec les enregistrements sismiques n'est pas estimée globalement satisfaisante par l'opérateur.
I-5 Pour chaque couple interface/vitesse et pour chaque collection multi-déports sélectionnée, on évalue la cohérence entre les temps de trajet multi-déports ainsi calculés et les enregistrements sismiques ; on peut par exemple employer une technique dite de calcul de semblante, pour laquelle on somme l'énergie rencontrée dans les données sismiques le long d'une courbe qui correspondra ici aux temps de trajet multi-déports calculés par tracé de rayons. Dans un cas idéal, il existera un échantillon de vitesse pour lequel la courbe de temps de trajet prédite se superposera parfaitement à
l'événement sismique considéré, ce qui se traduira par une mesure de cohérence maximale. Il est à noter que cette évaluation pourra se faire sur une gamme limitée de déport, définie par l'utilisateur ou automatiquement, par exemple par examen de la courbe montrant l'évolution de la cohérence en fonction de la gamme de déports.
Pour chaque collection multi-déports, on sélectionne la courbe de temps de trajet la plus en accord avec les enregistrements sismiques, et les temps de trajet multi-déports constituant cette courbe sont réservés. Il est à noter que la courbe de temps de trajet offrant la meilleure cohérence pourra être ajustée si nécessaire en cherchant dans une fenêtre verticale prédéfinie le maximum (ou le minimum en fonction de la polarité
de l'événement considéré) d'amplitude le plus proche pour chaque trace sismique de la collection multi-déport considérée.
I-6 Les temps de trajet multi-déports ainsi recueillis pour l'ensemble des collections multi-déports sélectionnées sont alors injectés dans une méthode d'inversion cinématique avant sommation telle qu'une tomographie de temps de trajet afin de déterminer la géométrie et la vitesse de la couche en question.
I-7 Il est possible que l'on soit obligé d'itérer les étapes précédentes dans le cas où la cohérence avec les enregistrements sismiques n'est pas estimée globalement satisfaisante par l'opérateur.
9 II) Cas particuliers II-1 Dans le cas où la distribution des vitesses dans les couches de la formation géologique est fortement variable latéralement et/ou dans le cas où l'on ne dispose pas d'un a priori suffisamment précis sur la distribution de vitesse dans la couche considérée, la méthode selon l'invention comporte les étapes suivantes :
L'ensemble des étapes précédentes est appliqué de façon itérative, les temps de trajet récoltés au cours d'une itération étant utilisés pour une mise à jour du modèle de vitesse par une méthode d'inversion cinématique avant sommation telle que la tomographie de temps de trajet, modèle qui sert alors de modèle de vitesse d'entrée pour une nouvelle itération de la méthode.
Il est à noter que, afin de s'affranchir au mieux des effets des variations latérales de vitesse qui auraient une longueur d'onde plus faible que le déport maximal des collections multi-déports, on pourra élargir la gamme de déports considérée par la méthode au fur et à
mesure des itérations. Plus précisément, on initiera la méthode sur une gamme de déport limitée, puis la distribution de vitesse trouvée à partir de cette gamme restreinte de déports servira de modèle de référence pour une nouvelle itération de la méthode lors de laquelle on considérera une gamme de déport plus large, et ainsi de suite. Il est à
noter que le choix de la gamme de déports considérée à une itération donnée pourra se faire par exemple par examen de la courbe montrant l'évolution de la cohérence en fonction de la gamme de déports.
Pour les cas où la distribution des vitesses varie très fortement latéralement, on pourra en plus adopter une approche multi-grilles, plus précisément, on initiera la méthode sur une grille lâche de collections multi-déports (afin de déterminer les grandes longueurs d'onde des variations de vitesse), puis l'on prendra en compte des grilles de collections multi-déports de plus en plus fines au fur et à mesure des itérations.
II-2 Dans le cas où la géométrie de l'interface recherchée est génératrice de triplications et si l'on souhaite également déterminer les zones de la couche considérée donnant lieu à ces triplications, il est avant tout nécessaire de disposer de méthodes d'inversion telles que des méthodes d'inversion cinématique (à vitesse fixée et avant sommation) par exemple permettant de prendre en compte des arrivées multi-valuées. Une telle méthode est décrite dans la référence suivante :
Delprat-Jannaud, F. et Lailly, P., (1995), How to handle multiple arrivais?
Journal of Geophysical Research, 100, no. B1, 703-715.
Les données d'entrée de ce type de méthode sont les temps de trajet multi-valués associés aux événements sismiques considérés, ainsi que les paramètres de rais associés à
5 ces temps de trajet. Afin d'obtenir les temps de trajet multi-valués et les paramètres de rais associés aux données à déport constant, on pourra avoir recours à une boucle de migration-démigration (comme décrite dans la méthode Ehinger and Lailly (1995) précédemment citée), en effectuant l'étape de démigration à l'aide d'un tracé de rayons permettant le calcul d'arrivées multi-valuées. Un tel tracé de rayons est décrit par exemple dans :
L'ensemble des étapes précédentes est appliqué de façon itérative, les temps de trajet récoltés au cours d'une itération étant utilisés pour une mise à jour du modèle de vitesse par une méthode d'inversion cinématique avant sommation telle que la tomographie de temps de trajet, modèle qui sert alors de modèle de vitesse d'entrée pour une nouvelle itération de la méthode.
Il est à noter que, afin de s'affranchir au mieux des effets des variations latérales de vitesse qui auraient une longueur d'onde plus faible que le déport maximal des collections multi-déports, on pourra élargir la gamme de déports considérée par la méthode au fur et à
mesure des itérations. Plus précisément, on initiera la méthode sur une gamme de déport limitée, puis la distribution de vitesse trouvée à partir de cette gamme restreinte de déports servira de modèle de référence pour une nouvelle itération de la méthode lors de laquelle on considérera une gamme de déport plus large, et ainsi de suite. Il est à
noter que le choix de la gamme de déports considérée à une itération donnée pourra se faire par exemple par examen de la courbe montrant l'évolution de la cohérence en fonction de la gamme de déports.
Pour les cas où la distribution des vitesses varie très fortement latéralement, on pourra en plus adopter une approche multi-grilles, plus précisément, on initiera la méthode sur une grille lâche de collections multi-déports (afin de déterminer les grandes longueurs d'onde des variations de vitesse), puis l'on prendra en compte des grilles de collections multi-déports de plus en plus fines au fur et à mesure des itérations.
II-2 Dans le cas où la géométrie de l'interface recherchée est génératrice de triplications et si l'on souhaite également déterminer les zones de la couche considérée donnant lieu à ces triplications, il est avant tout nécessaire de disposer de méthodes d'inversion telles que des méthodes d'inversion cinématique (à vitesse fixée et avant sommation) par exemple permettant de prendre en compte des arrivées multi-valuées. Une telle méthode est décrite dans la référence suivante :
Delprat-Jannaud, F. et Lailly, P., (1995), How to handle multiple arrivais?
Journal of Geophysical Research, 100, no. B1, 703-715.
Les données d'entrée de ce type de méthode sont les temps de trajet multi-valués associés aux événements sismiques considérés, ainsi que les paramètres de rais associés à
5 ces temps de trajet. Afin d'obtenir les temps de trajet multi-valués et les paramètres de rais associés aux données à déport constant, on pourra avoir recours à une boucle de migration-démigration (comme décrite dans la méthode Ehinger and Lailly (1995) précédemment citée), en effectuant l'étape de démigration à l'aide d'un tracé de rayons permettant le calcul d'arrivées multi-valuées. Un tel tracé de rayons est décrit par exemple dans :
10 - Jurado, F., Lailly, P., and Ehinger, A., (1998), Fast 3D two-point raytracing for traveltime tomography: Proceedings of SPIE, Mathematical Methods in Geophysical Imaging V, 3453, 70-81.
Puis, après avoir constitué les différents couples interface/vitesse par inversion cinématique à vitesse fixée sur les données temps/paramètres de rais ainsi obtenues, on calcule l'information cinématique, comme décrit dans le cas le plus général, par le biais par exemple de tracés des rayons multi-déports dans chacun de ces modèles et pour chacune des collections multi-déports sélectionnées, mais cette fois, adapté aux calculs des arrivées multi-valuées tel que celui décrit dans Jurado et al. (1998) précédemment cité. On obtient ainsi, pour chaque collection multi-déports sélectionnée, à la fois des temps de trajet multi-valués avant sommation et leurs paramètres de rais associés.
Puis, l'analyse de cohérence entre les courbes de temps de trajet ainsi prédites et les événements considérés dans les données sismiques sera réalisée tel que décrit dans le cas le plus général, mais de façon indépendante pour chacune des branches des arrivées multi-valuées. Les temps de trajet ainsi récoltés pour chacune des branches des arrivées multi-valuées seront alors injectés, avec leur paramètres de rais associés, dans une méthode d'inversion cinématique avant sommation permettant de prendre en compte les arrivées multi-valuées, telle que celle décrite dans Delprat-Jannaud et Lailly (1995) précédemment citée.
On a décrit ici une application de la méthode à la détermination d'un modèle de vitesse des ondes sismiques dans une formation souterraine. Il est bien évident toutefois que la méthode peut tout aussi bien s'appliquer à la détermination de la vitesse de
Puis, après avoir constitué les différents couples interface/vitesse par inversion cinématique à vitesse fixée sur les données temps/paramètres de rais ainsi obtenues, on calcule l'information cinématique, comme décrit dans le cas le plus général, par le biais par exemple de tracés des rayons multi-déports dans chacun de ces modèles et pour chacune des collections multi-déports sélectionnées, mais cette fois, adapté aux calculs des arrivées multi-valuées tel que celui décrit dans Jurado et al. (1998) précédemment cité. On obtient ainsi, pour chaque collection multi-déports sélectionnée, à la fois des temps de trajet multi-valués avant sommation et leurs paramètres de rais associés.
Puis, l'analyse de cohérence entre les courbes de temps de trajet ainsi prédites et les événements considérés dans les données sismiques sera réalisée tel que décrit dans le cas le plus général, mais de façon indépendante pour chacune des branches des arrivées multi-valuées. Les temps de trajet ainsi récoltés pour chacune des branches des arrivées multi-valuées seront alors injectés, avec leur paramètres de rais associés, dans une méthode d'inversion cinématique avant sommation permettant de prendre en compte les arrivées multi-valuées, telle que celle décrite dans Delprat-Jannaud et Lailly (1995) précédemment citée.
On a décrit ici une application de la méthode à la détermination d'un modèle de vitesse des ondes sismiques dans une formation souterraine. Il est bien évident toutefois que la méthode peut tout aussi bien s'appliquer à la détermination de la vitesse de
11 propagation d'autres types d'ondes dans un milieu hétérogène présentant des discontinuités.
Claims (10)
1. Méthode pour transformer des enregistrements sismiques temporels en une image profondeur d'une formation souterraine, lesdits enregistrements sismiques correspondant à des enregistrements multi-déports d'ondes sismiques captées par des récepteurs sismiques couplés avec la formation souterraine, en réponse à
l'émission d'ondes sismiques dans ladite formation par une source sismique, après réflexion sur des interfaces géologiques de la dite formation, dans laquelle on détermine un modèle de vitesse desdites ondes sismiques, caractérisée en ce qu'elle comporte au moins les étapes suivantes pour chaque événement sismique repéré sur les enregistrements et pour chaque couche délimitée par les dites interfaces a) à partir des enregistrements sismiques avant sommation, on construit une collection iso-déport, dont on extrait une information cinématique ou des temps de trajet associés à l'événement;
b) on choisit une gamme de vitesse autour d'une vitesse de référence dans ladite couche que l'on échantillonne avec un pas défini, c) pour chaque échantillon de vitesse, on applique une technique d'inversion à
vitesse fixée afin de déterminer, à partir des temps de trajet extraits de la collection iso-déport, la géométrie de la dite interface pour l'échantillon de vitesse en question de façon à obtenir une série de couples interface/vitesse pour le dit événement;
d) on calcule l'information cinématique associée à chaque couple interface/vitesse obtenu, pour des couples source-récepteur correspondants à
des collections multi-déports existant dans les enregistrements sismiques;
e) pour chaque couple interface/vitesse et pour chaque collection multi-déports sélectionnée, on évalue la cohérence entre les temps de trajet multi-déports ainsi calculés et les enregistrements sismiques, et on sélectionne pour chaque collection multi-déports la courbe de temps de trajet offrant la cohérence maximale avec les enregistrements sismiques;
f) on applique une méthode d'inversion cinématique avant sommation en utilisant les temps de trajet multi-déports obtenus pour l'ensemble des collections multi-déports sélectionnées, afin de déterminer la géométrie et la vitesse de la couche considérée;
g) on itère n fois (n >= 0) les étapes a) à f) en considérant à chaque itération le modèle de vitesse obtenu à l'itération précédente comme modèle de référence pour définir la dite vitesse de référence; et h) on utilise ledit modèle de vitesse obtenu à l'issue de l'étape g), pour construire ladite image profondeur de la formation souterraine.
l'émission d'ondes sismiques dans ladite formation par une source sismique, après réflexion sur des interfaces géologiques de la dite formation, dans laquelle on détermine un modèle de vitesse desdites ondes sismiques, caractérisée en ce qu'elle comporte au moins les étapes suivantes pour chaque événement sismique repéré sur les enregistrements et pour chaque couche délimitée par les dites interfaces a) à partir des enregistrements sismiques avant sommation, on construit une collection iso-déport, dont on extrait une information cinématique ou des temps de trajet associés à l'événement;
b) on choisit une gamme de vitesse autour d'une vitesse de référence dans ladite couche que l'on échantillonne avec un pas défini, c) pour chaque échantillon de vitesse, on applique une technique d'inversion à
vitesse fixée afin de déterminer, à partir des temps de trajet extraits de la collection iso-déport, la géométrie de la dite interface pour l'échantillon de vitesse en question de façon à obtenir une série de couples interface/vitesse pour le dit événement;
d) on calcule l'information cinématique associée à chaque couple interface/vitesse obtenu, pour des couples source-récepteur correspondants à
des collections multi-déports existant dans les enregistrements sismiques;
e) pour chaque couple interface/vitesse et pour chaque collection multi-déports sélectionnée, on évalue la cohérence entre les temps de trajet multi-déports ainsi calculés et les enregistrements sismiques, et on sélectionne pour chaque collection multi-déports la courbe de temps de trajet offrant la cohérence maximale avec les enregistrements sismiques;
f) on applique une méthode d'inversion cinématique avant sommation en utilisant les temps de trajet multi-déports obtenus pour l'ensemble des collections multi-déports sélectionnées, afin de déterminer la géométrie et la vitesse de la couche considérée;
g) on itère n fois (n >= 0) les étapes a) à f) en considérant à chaque itération le modèle de vitesse obtenu à l'itération précédente comme modèle de référence pour définir la dite vitesse de référence; et h) on utilise ledit modèle de vitesse obtenu à l'issue de l'étape g), pour construire ladite image profondeur de la formation souterraine.
2. Méthode selon la revendication 1, caractérisée en ce que l'on réalise n fois l'étape g) avec n >= 1 dans le cas où la gamme de vitesse choisie ne s'avère pas assez précise à l'issue de l'itération précédente où des étapes a) à f).
3. Méthode selon la revendication 1 ou 2, caractérisée en ce que, dans le cas où
la distribution des vitesses varie fortement latéralement et/ou dans le cas où
l'on ne dispose pas d'une connaissance a priori suffisamment précise sur la distribution de vitesse dans la couche considérée, on réalise l'étape g) sur des gammes de déports de plus en plus larges au fur et à mesure des itérations.
la distribution des vitesses varie fortement latéralement et/ou dans le cas où
l'on ne dispose pas d'une connaissance a priori suffisamment précise sur la distribution de vitesse dans la couche considérée, on réalise l'étape g) sur des gammes de déports de plus en plus larges au fur et à mesure des itérations.
4. Méthode selon l'une quelconque des revendications 1 à 3, caractérisée en ce que, dans le cas où la distribution des vitesses varie fortement latéralement et/ou dans le cas où l'on ne dispose pas d'une connaissance a priori suffisamment précise sur la distribution de vitesse dans la couche considérée, on réalise l'étape g) sur des grilles de collections multidéports de plus en plus fines au fur et à mesure des itérations.
5. Méthode selon l'une quelconque des revendications 1 à 4, caractérisée en ce que, à l'étape e), dans le cas où la géométrie de l'interface recherchée est génératrice de triplications, on considère les branches de triplications dans les collections multi-déports indépendamment les unes des autres.
6. Méthode selon la revendication 5, caractérisée en ce que l'on réalise l'étape e) en utilisant des outils de tracé de rayons et d'inversion permettant de prendre en compte des arrivées multi-valuées.
7. Méthode selon l'une quelconque des revendications 1 à 6, caractérisée en ce que l'on construit à l'étape a) une collection iso-déport à déport nul ou à
déport proche.
déport proche.
8. Méthode selon l'une quelconque des revendications 1 à 7, caractérisée en ce que l'on applique à l'étape c) une technique d'inversion cinématique à vitesse fixée telle qu'une migration de carte.
9. Méthode selon l'une quelconque des revendications 1 à 8, caractérisée en ce que l'on calcule à l'étape d) l'information cinématique en traçant des rayons multi-déports sur l'interface de chaque couple interface-vitesse.
10. Méthode selon l'une quelconque des revendications 1 à 9, caractérisée en ce que l'on réalise l'étape f) en appliquant une méthode d'inversion cinématique avant sommation telle qu'une tomographie de temps de trajet avant sommation.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR03/05.292 | 2003-04-29 | ||
| FR0305292A FR2854468B1 (fr) | 2003-04-29 | 2003-04-29 | Methode pour determiner un modele de vitesse d'ondes sismiques dans une formation souterrraine heterogene |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA2464799A1 CA2464799A1 (fr) | 2004-10-29 |
| CA2464799C true CA2464799C (fr) | 2013-01-15 |
Family
ID=32339015
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA2464799A Expired - Fee Related CA2464799C (fr) | 2003-04-29 | 2004-04-26 | Methode pour determiner un modele de vitesse d'ondes sismiques dans une formation souterraine heterogene |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US6868039B2 (fr) |
| CA (1) | CA2464799C (fr) |
| FR (1) | FR2854468B1 (fr) |
| GB (1) | GB2401966B (fr) |
| IL (1) | IL161638A (fr) |
| NL (1) | NL1026045C2 (fr) |
| NO (1) | NO20041734L (fr) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7388808B2 (en) * | 2004-09-23 | 2008-06-17 | Pgs Americas, Inc. | Method for depth migrating seismic data using pre-stack time migration, demigration, and post-stack depth migration |
| AU2006235820B2 (en) * | 2005-11-04 | 2008-10-23 | Westerngeco Seismic Holdings Limited | 3D pre-stack full waveform inversion |
| US8203907B2 (en) * | 2006-06-12 | 2012-06-19 | Exxonmobile Upstream Research Company | Updating velocity models using migration velocity scans |
| JP5028073B2 (ja) * | 2006-11-29 | 2012-09-19 | 株式会社ニデック | 角膜手術装置 |
| US7688674B2 (en) * | 2007-03-05 | 2010-03-30 | Schlumberger Technology Corporation | Methods and apparatus for performing moving checkshots |
| EP2459847A2 (fr) * | 2010-06-10 | 2012-06-06 | Hipoint Reservoir Imaging | Cartographie de réservoir à signal d'impulsion de fracture |
| CN105308479B (zh) * | 2013-05-24 | 2017-09-26 | 埃克森美孚上游研究公司 | 通过与偏移距相关的弹性fwi的多参数反演 |
| FR3028649B1 (fr) | 2014-11-18 | 2016-12-02 | Ifp Energies Now | Procede de construction d'un modele geologique |
| CN108139498B (zh) * | 2015-10-15 | 2019-12-03 | 埃克森美孚上游研究公司 | 具有振幅保持的fwi模型域角度叠加 |
| US11221425B1 (en) | 2020-07-08 | 2022-01-11 | Saudi Arabian Oil Company | Generating a model for seismic velocities in a subsurface region using inversion with lateral variations |
| CN113917524B (zh) * | 2020-07-10 | 2024-03-29 | 中国石油化工股份有限公司 | 基于无缆节点台站的近地表精细探测方法 |
| US11467299B2 (en) | 2020-12-16 | 2022-10-11 | Saudi Arabian Oil Company | Full waveform inversion velocity guided first arrival picking |
| CN114063153B (zh) * | 2021-12-01 | 2022-06-10 | 中国地震局地球物理研究所 | 一种自动反演震源机制解的方法与装置 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5197039A (en) * | 1988-03-29 | 1993-03-23 | Shell Oil Company | Methods for processing seismic data |
| US5136553A (en) * | 1990-12-19 | 1992-08-04 | Amoco Corporation | Method of geophysical exploration |
| US5570321A (en) * | 1994-03-03 | 1996-10-29 | Atlantic Richfield Company | Seismic velocity model optimization method using simulated annearling to determine prestack travel-times |
| GB9607764D0 (en) * | 1996-04-15 | 1996-06-19 | Geco As | Inversion method for seismic data |
| FR2821677B1 (fr) * | 2001-03-05 | 2004-04-30 | Geophysique Cie Gle | Perfectionnements aux procedes d'inversion tomographique d'evenements pointes sur les donnees sismiques migrees |
| FR2838831B1 (fr) * | 2002-04-17 | 2004-08-20 | Inst Francais Du Petrole | Methode pour determiner les corrections statiques primaires a appliquer a des traces sismiques |
| US7072767B2 (en) * | 2003-04-01 | 2006-07-04 | Conocophillips Company | Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data |
-
2003
- 2003-04-29 FR FR0305292A patent/FR2854468B1/fr not_active Expired - Fee Related
-
2004
- 2004-04-23 US US10/830,052 patent/US6868039B2/en not_active Expired - Fee Related
- 2004-04-26 IL IL161638A patent/IL161638A/en not_active IP Right Cessation
- 2004-04-26 NL NL1026045A patent/NL1026045C2/nl not_active IP Right Cessation
- 2004-04-26 GB GB0409175A patent/GB2401966B/en not_active Expired - Fee Related
- 2004-04-26 CA CA2464799A patent/CA2464799C/fr not_active Expired - Fee Related
- 2004-04-28 NO NO20041734A patent/NO20041734L/no not_active Application Discontinuation
Also Published As
| Publication number | Publication date |
|---|---|
| NL1026045C2 (nl) | 2005-07-26 |
| NL1026045A1 (nl) | 2004-11-01 |
| NO20041734L (no) | 2004-11-01 |
| IL161638A (en) | 2008-04-13 |
| GB2401966B (en) | 2007-10-24 |
| CA2464799A1 (fr) | 2004-10-29 |
| FR2854468B1 (fr) | 2005-06-10 |
| GB0409175D0 (en) | 2004-05-26 |
| FR2854468A1 (fr) | 2004-11-05 |
| IL161638A0 (en) | 2004-09-27 |
| US6868039B2 (en) | 2005-03-15 |
| GB2401966A (en) | 2004-11-24 |
| US20040243313A1 (en) | 2004-12-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Plessix et al. | Multiparameter full-waveform inversion: Marine and land examples | |
| US5596548A (en) | Seismic imaging using wave equation extrapolation | |
| EP0889331B1 (fr) | Méthode pour modéliser en 3D l'impédance d'un milieu hétérogène | |
| EP0789849B1 (fr) | Methode pour representer dans un domaine spatio-temporel la trajectoire d'au moins un puits de forage | |
| CA2733882C (fr) | Procede de surveillance d'un site de stockage geologique de gaz par inversion stratigraphique de donnees sismiques | |
| EP1669778B1 (fr) | Méthode pour déterminer de l'information spéculaire dans l'imagerie sismique avant sommation | |
| CA2464799C (fr) | Methode pour determiner un modele de vitesse d'ondes sismiques dans une formation souterraine heterogene | |
| EP2253970B1 (fr) | Méthode pour imager une zone cible du sous-sol à partir de données de type walkaway | |
| CA2670384A1 (fr) | Methode d'inversion conjointe de donnees sismiques representees sur des echelles de temps differentes | |
| FR2708350A1 (fr) | Méthode d'analyse de signaux sismiques. | |
| FR2916859A1 (fr) | Procede de traitement de donnees sismiques | |
| FR2918178A1 (fr) | Methode pour ajuster un modele de vitesse d'ondes sismiques en fonction d'informations relevees aux puits | |
| FR2901887A1 (fr) | Construction d'un modele d'analyse des vitesses dans les formations infra-saliferes | |
| CA2425343A1 (fr) | Methode pour determiner les corrections statiques primaires a appliquer a des traces sismiques | |
| CN105093296A (zh) | 一种优化观测系统的方法及装置 | |
| CA2456596C (fr) | Methode pour mesurer les similarites locales entre plusieurs cubes de traces sismiques | |
| EP0698220B1 (fr) | Procede de modelisation de donnees sismiques cinematiques ayant subi un traitement par au moins un operateur de deplacement | |
| Li et al. | Prestack seismic inversion with structural constraints | |
| Pipan et al. | Multifold ground-penetrating radar and resistivity to study the stratigraphy of shallow unconsolidated sediments | |
| Smith et al. | Making seismic monitoring work in a complex desert environment—4D processing | |
| Liu et al. | Wave-equation diffraction imaging using pseudo dip-angle gather | |
| Yilmaz | An image-based effective-medium modeling of near-surface anomalies | |
| FR2675590A1 (fr) | Procede de traitement pour l'obtention d'une section somme a offsets nuls. | |
| Nakata et al. | Usage of pre-stack vibroseis data before sweep correlation for noise removal and near-surface structural analysis. | |
| WO2003003052A1 (fr) | Perfectionnements aux techniques de pointe de surface |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request | ||
| MKLA | Lapsed |
Effective date: 20160426 |