FR2845773A1 - Procede de determination d'un parametre auquel on s'interesse d'une partie d'un reservoir dans des formations terrestres - Google Patents

Procede de determination d'un parametre auquel on s'interesse d'une partie d'un reservoir dans des formations terrestres Download PDF

Info

Publication number
FR2845773A1
FR2845773A1 FR0311610A FR0311610A FR2845773A1 FR 2845773 A1 FR2845773 A1 FR 2845773A1 FR 0311610 A FR0311610 A FR 0311610A FR 0311610 A FR0311610 A FR 0311610A FR 2845773 A1 FR2845773 A1 FR 2845773A1
Authority
FR
France
Prior art keywords
well
tank
seismic
signals
reservoir
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
FR0311610A
Other languages
English (en)
Other versions
FR2845773B1 (fr
Inventor
Peter S Aronstam
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Baker Hughes Holdings LLC
Original Assignee
Baker Hughes Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Baker Hughes Inc filed Critical Baker Hughes Inc
Publication of FR2845773A1 publication Critical patent/FR2845773A1/fr
Application granted granted Critical
Publication of FR2845773B1 publication Critical patent/FR2845773B1/fr
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/42Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators in one well and receivers elsewhere or vice versa

Abstract

L'invention concerne un procédé pour déterminer un paramètre auquel on s'intéresse d'une partie d'un réservoir dans des formations terrestres. On déploie plusieurs récepteurs sismiques dans un premier puits et d'autres récepteurs dans un second puits disposés de part et d'autre de la partie du réservoir à laquelle on s'intéresse. Des ondes sismiques sont générées en plusieurs emplacements à la surface ou au voisinage de la surface et des signaux représentatifs des ondes sismiques qui se sont propagées sont reçus par plusieurs récepteurs situés dans les puits, et d'autres signaux correspondant aux ondes sismiques réfléchies sont reçus par d'autres récepteurs situés à la surface ou au voisinage de celle-ci. Une analyse tomographique des premiers et seconds signaux est ensuite réalisée pour obtenir le paramètre auquel on s'intéresse. Domaine d'application : prospection sismique

Description

Domaine de l'invention L'invention concerne de façon générale le domaine
de l'interprétation de la surveillance de réservoirs des mesures sismiques à laps de temps. La présente invention 5 concerne plus particulièrement une configuration nouvelle de sources et de récepteurs pour produire une reconstruction tomographique améliorée d'images de réservoirs. Contexte technique Des campagnes géophysiques sont utilisées pour découvrir l'étendue de gisements minéraux sous la surface tels que des gisements de pétrole, de gaz naturel, d'eau, de soufre, etc... Des procédés géophysiques peuvent également être utilisés pour surveiller des changements se produisant 15 dans le gisement, tel qu'un épuisement résultant de la production du minéral sur la durée de vie naturelle du gisement qui peut être de nombreuses années. L'utilité d'une étude géophysique dépend de l'aptitude à mesurer quantitativement et évaluer un certain analogue géophysique 20 d'un paramètre pétrophysique qui est lié directement à la
présence du minéral considéré.
L'imagerie sismique classique puits-à-puits utilise typiquement une paire de sondages à proximité d'un réservoir auquel on s'intéresse. Dans le premier de ces 25 sondages, on déploie une source sismique pour émettre de l'énergie sismique dans la région à laquelle on s'intéresse. Pour des raisons indiquées ci-dessous, des sources à fréquences de balayage sont communément utilisées. La source est déplacée séquentiellement de façon 30 à passer par une série de positions dans le premier sondage
et un signal sismique est généré en chaque position.
L'énergie sismique passe à travers la formation souterraine à laquelle on s'intéresse pour arriver au second de la paire de sondages. Un dispositif de récepteurs est 35 habituellement déployé dans le second sondage et, de même que pour la source sismique, le dispositif de récepteurs
est déplacé de façon à passer par une série de positions à l'intérieur du second sondage. En émettant un signal depuis chaque position de source dans le premier sondage et en recevant des données depuis chaque position de source dans 5 chaque position de récepteur dans le second sondage, on génère un ensemble de données sismiques puits-à-puits.
En référence à la figure 1, la formation d'images puits-à-puits est habituellement effectuée en utilisant un système d'imagerie puits-à-puits qui est indiqué de façon 10 générale par la référence numérique 10. Dans le système 10, de l'énergie sismique 12 est émise à travers une région souterraine 14 du sol en utilisant une source 16 qui est positionnée dans un premier sondage 18. Une source sismique mobile déplacée de façon à passer par des emplacements 15 Si... Sn est utilisée pour générer des ondes sismiques
désignées en 12 qui se propagent à partir du premier sondage 18. Les ondes sismiques sont détectées dans le second sondage 24 par des détecteurs convenables.
Habituellement, un dispositif 22 de récepteurs comportant 20 des récepteurs 23a, 23b... 23e est utilisé à cet effet. Le dispositif de récepteurs peut être déplacé par des positions désignées par Rl... Rs. Des rais 28a... 28e et 30a... 30e montrent des exemples des trajectoires des ondes sismiques se propageant depuis divers emplacements de 25 source vers les différents récepteurs. La tomographie sismique a pour objectif de mesurer les temps de parcours pour plusieurs combinaisons sourcerécepteurs et, à partir de la géométrie connue des emplacements de la source et des récepteurs, de déterminer le champ de vitesse des ondes se 30 propageant dans la région 10. Des procédés pour convertir les temps de parcours mesurés afin d'obtenir le champ de
vitesse sont connus dans la technique.
Dans la surveillance d'un réservoir, la région 10 comprend un réservoir à partir duquel on recueille des 35 hydrocarbures. Dans la récupération naturelle d'hydrocarbures de même que dans la récupération secondaire dans laquelle un fluide est injecté à partir d'un puits d'injection dans le réservoir, un fluide est remplacé en continu par un autre fluide dans le réservoir. De plus, la pression des fluides dans le réservoir peut également faire 5 l'objet de changements. Il est connu que la substitution des fluides de même que les changements de pression provoquent des changements des vitesses sismiques. Par conséquent, en surveillant les changements dans le domaine des vitesses, on peut en déduire certaines conclusions 10 concernant des modifications de la distribution des fluides dans le réservoir et/ou des modifications de la pression dans celui-ci. Ce processus est appelé " tomographie
séquentielle " (" time-lapse tomography ").
Habituellement, des sources sismiques de fond sont des 15 sources à fréquences balayées. La raison pour laquelle on utilise des sources à fréquences balayées est que ce sont des sources de basse puissance qui risquent moins de provoquer des détériorations des sondages que des sources à impulsions à haute puissance. On a développé un certain 20 nombre d'autres sources sismiques à basse puissance pour
permettre la saisie de données sismiques puits-à-puits.
Elles comprennent des résonateurs, des transducteurs piézoélectriques et des transducteurs magnétostrictifs.
D'autres sources telles que des sources à explosion ou à 25 implosion, des canons pneumatiques de fond et des étincelleurs ont généralement une puissance de sortie plus élevée et risquent davantage de provoquer une détérioration du puits de forage. Cette limitation à une basse puissance restreint l'utilisation de la tomographie sismique puits-à30 puits à des situations dans lesquelles la distance entre puits est relativement courte (c'est-à-dire pas plus
d'environ 300 m (1000 feet)).
Une autre limitation de l'utilisation de la tomographie sismique puits-àpuits est le cot substantiel 35 associé à la préparation des puits pour déployer des sources sismiques et des récepteurs de fond. Dans la
plupart des cas, les diamètres des sources sismiques et des récepteurs de fond sont trop grands pour se loger dans la colonne de production utilisée pour transporter des fluides du réservoir à la surface. Par conséquent, sur la plupart 5 des sites, la colonne de production doit d'abord être enlevée des puits afin d'exécuter une reconnaissance puitsà- puits, puis remise en place après l'achèvement de la reconnaissance. L'enlèvement de la colonne de production s'accompagne du cot d'arrêt de la production pendant 10 l'étude sismique.
Le brevet des E.U.A. n0 5 481 501 de Blakeslee, dont le contenu est incorporé ici en totalité par référence, décrit un procédé dans lequel aucune source n'est utilisée dans le puits de forage. Le procédé enseigné par Blakeslee 15 est illustré sur la figure 2 qui montre deux puits, un puits A et un puits B, descendant dans la terre 110. Une source sismique s est placée sur la surface 112 de la terre 110, sensiblement en alignement avec les puits A et B, mais non entre eux. Plusieurs récepteurs sismiques de fond, 20 a,... a,. Avantageusement, plusieurs récepteurs sismiques de fond, b1... bn. Pour simuler des données puits-à- puits entre les puits A et B, la source sismique s est activée afin de générer un signal sismique qui se propage à travers les formations souterraines et est enregistrée par chacun des 25 récepteurs sismiques a,.
an du premier puits et chacun des récepteurs sismiques b1... bn du second puits. Les données sismiques résultantes sont ensuite traitées pour fournir une information concernant la région intermédiaire entre les puits A et B. Le procédé utilisé par Blakeslee utilise le principe de Fermat pour un rayon parcourant le trajet de la source s à un récepteur tel que bk sur la figure 2. Si le domaine de vitesse entre les puits A et B est connu, le temps de parcours depuis chacune des positions a1... an jusqu'au 35 récepteur bk peut alors être calculé. D'après le principe de Fermat, le temps de parcours mesuré de s à bk, désigné par TSbk est alors donné par Ts,bk = Minj (Ts,aj + Tajbk) (1) Un raisonnement similaire s'applique au cas montré sur la figure 3 o le puits B est situé entre le puits A et la source s. Sur la figure 3, le temps de parcours Ts, ak de la source s jusqu'à un récepteur ak est donné par Tsak = Minj(Tsbi + Tbjak) (2) Le procédé de Blakeslee suppose un domaine de vitesse, des contrôles pour déterminer si les équations (1) et (2) 15 sont satisfaites pour toutes les positions des récepteurs, et, si elles ne sont pas satisfaites, il modifie de façon itérative le domaine de vitesse jusqu'à ce que l'équation soit approximativement satisfaite dans certaines limites d'erreur. En utilisant le procédé de Blakeslee, il n'est pas..DTD: nécessaire d'utiliser des sources de fond. Par conséquent, il est possible d'utiliser des sources sismiques classiques de surface pour obtenir des données qui peuvent être utilisées pour une reconstruction tomographique d'une 25 région souterraine.
L'un des problèmes de la tomographie puits-à-puits (réelle ou simulée) est que les trajets des rayons ne couvrent qu'une plage limitée d'angles et qu'on manque gravement de trajets de rayons proches de la verticale. 30 L'effet de l'ouverture limitée est que la reconstruction
tomographique peut avoir une faible résolution verticale.
On a besoin d'un procédé d'obtention de données pour une analyse tomographique qui ne nécessite pas de sources de fond, sans en même temps être affecté d'une mauvaise 35 résolution verticale. La présente invention satisfait à ce besoin.
SOMMAIRE DE L'INVENTION
La présente invention est un procédé de contrôle d'un réservoir dans des formations terrestres utilisant des récepteurs sismiques dans au moins deux puits. Des sources 5 sismiques de surface sont utilisées pour simuler des données puits-à-puits entre les deux puits et, en utilisant des récepteurs de surface, on obtient des données sismiques de réflexion sur la région entre les deux puits. Les données de réflexion sont avantageusement obtenues à partir 10 de réflecteurs situés à la fois audessus et au-dessous du réservoir. Une analyse tomographique des données puits-àpuits simulées et de réflexion procure un modèle du réservoir. Les données enregistrées dans les puits peuvent être utilisées pour une correction des données de 15 réflexion. Ce processus est répété à une époque ultérieure après une période de production à partir du réservoir. Des changements dans le modèle sont représentatifs de changements de fluide du réservoir, de changements de
pression ou d'un tassement.
La simulation de données puits-à-puits est réalisée en
utilisant le principe de Fermat. Dans une forme de réalisation facultative de l'invention, une atténuation des signaux sismiques est utilisée pour contrôler le réservoir.
BREVE DESCRIPTION DES DESSINS
La figure 1 (art antérieur) est une élévation avec vue en coupe transversale d'une région souterraine s'étendant entre deux sondages, et illustrant un procédé de l'art antérieur pour la génération d'une image de vitesse sismique. Les figures 2 et 3 (art antérieur) illustrent schématiquement la présente invention pour simuler des données sismiques puits-à-puits entre des puits A et B, la source sismique s étant plus proche du puits A sur la
figure 1 et plus proche du puits B sur la figure 2.
La figure 4 illustre la géométrie d'acquisition pour
une tomographie séquentielle de la présente invention.
La figure 5 montre comment le procédé de la présente invention peut être utilisé pour contrôler un volume en
trois dimensions d'un réservoir.
DESCRIPTION DETAILLEE DE L'INVENTION
Le procédé d'acquisition de données pour le contrôle d'un réservoir selon la présente invention est illustré schématiquement sur la figure 4. L'acquisition de données comprend trois particularités principales. De même que dans le procédé de Blakeslee, des récepteurs ai... an et bi... bn 10 sont déployés dans deux puits. Deux puits seulement sont représentés pour simplifier l'illustration et, en réalité, on pourrait utiliser davantage de puits. En utilisant des sources sismiques de surface, illustrées en 203a, 203b... 203k, on simule des données sismiques puits-à-puits 15 en utilisant le procédé décrit par Blakeslee. Cependant,
contrairement aux enseignements de Blakeslee, des sources de surface sont également déployées entre les deux puits.
Ceci produit des rayons supplémentaires pour une reconstruction tomographique de la région entre les deux 20 puits et aborde le problème de l'ouverture limitée avec la tomographie puits-à-puits classique ou simulée. Un exemple de trajet de rayons pour une telle source entre deux emplacements de puits est illustré en 233. Comme précédemment, pour simplifier la figure, un seul trajet de 25 rayons est représenté allant d'une source unique jusqu'à un détecteur unique de fond; en réalité, des signaux partant de plusieurs sources de surface entre les deux puits seraient enregistrés au niveau de plusieurs détecteurs dans
les deux puits.
Une troisième particularité de l'acquisition de données dans la présente invention est l'utilisation d'une tomographie de réflexion utilisant des sources telles que 221a et des récepteurs tels que 223b. Pour simplifier l'illustration, une seule source et un seul récepteur sont 35 représentés pour la tomographie de réflexion; en réalité, il y aurait plusieurs sources et plusieurs récepteurs. Des données de réflexion provenant d'un réflecteur tel que 235, 207 et 205 sont saisies. Pour simplifier l'illustration, trois réflecteurs seulement sont représentés. En réalité, il y aurait plusieurs récepteurs supplémentaires (non 5 représentés). Des exemples de trajets de rayons 211 et 231 provenant de la source 221a, qui reviennent respectivement au récepteur 223b en étant réfléchis en 205 et 231, sont représentés. Aux fins de l'invention, il est nécessaire d'avoir au moins un réflecteur tel que 235 se trouvant au10 dessus du réservoir de production. De plus, il est nécessaire d'avoir un réflecteur 205 se trouvant en dessous du réservoir de production. Des réflecteurs supplémentaires tels que 207 à l'intérieur du réservoir de production sont souhaitables. Les raisons de la présence de réflecteurs 15 235, 205 encadrant le réservoir de production sont décrites ci-après. Ces mesures sont exécutées en des périodes distinctes sur une période de temps prolongée, appelées ciaprès " époques ". L'objectif, dans le contrôle du réservoir, est d'utiliser des mesures réalisées à 20 différentes époques pour déduire une information concernant
le changement de la distribution des fluides dans un réservoir. A cet effet, il est courant d'effectuer une inversion tomographique de données mesurées à chacune des époques et de comparer les résultats des inversions 25 tomographiques d'une époque à une autre.
En principe, on pourrait penser qu'en utilisant des mesures de temps de parcours sensiblement verticaux jusqu'à des réflecteurs souterrains, l'ouverture limitée indiquée précédemment en ce qui concerne la tomographie puits-à30 puits pourrait être élargie, aboutissant à une reconstruction tomographique améliorée. Cependant, ceci n'a pas été réalisé pour les raisons qui sont maintenant indiquées. L'utilisation d'une tomographie de réflexion a été 35 décrite dans de nombreuses publications. Bishop et collaborateurs traitent d'une tomographique sismique de
réflexion sous la forme d'un algorithme de Gauss-Newton itératif qui produit un modèle de vitesse en fonction de la profondeur qui minimise la différence entre des temps de parcours générés en traçant des rayons à travers le modèle 5 et des temps de parcours mesurés à partir des données.
Cependant, comme indiqué par Krebs (brevet des E.U.A.
n0 6 002 642) " malheureusement, les données sismiques de surface ne contiennent pas suffisamment d'informations pour spécifier de façon unique à la fois un modèle de vitesse de 10 migration et les géométries des réflecteurs. Il en résulte que le modèle de vitesse dérivé peut être ambigu ou non raisonnable du point de vue géologique. On peut apporter des améliorations en appliquant des contraintes au processus d'optimisation, mais ces contraintes réduisent ou 15 éliminent généralement les ambigutés au prix d'adaptations
plus mauvaises aux données de temps de parcours ".
Un second problème posé par la tomographie de réflexion est mis en évidence par Squires et collaborateurs. Ceci est expliqué dans la présente 20 invention en revenant à la figure 4 qui montre une couche 202 soumise aux intempéries, proche de la surface. Une telle couche soumise aux intempéries apparaît pratiquement en tout point dans la proche surface du terrain du fait des intempéries. Cette couche soumise aux intempéries est 25 caractérisée par une épaisseur irrégulière et une vitesse extrêmement faible. A titre de référence, des roches sédimentaires normales ont des vitesses d'ondes de compression qui vont d'environ 1,5 km/s jusqu'à 6 km/s. Par contre, la couche soumise aux intempéries peut avoir des 30 vitesses descendant jusqu'à 100 m/s. Il résulte de cette vitesse extrêmement faible que les trajets des rayons dans la couche soumise aux intempéries sont pratiquement verticaux, comme indiqué sur la figure 4. L'effet de l'épaisseur irrégulière de la couche soumise aux 35 intempéries est d'introduire un retard dans les trajets des rayons de réflexion qui dépendent uniquement de l'emplacement de la source (ou du récepteur) et non de la distance source-récepteurs. Du fait de l'absence de dépendance vis-à-vis du déport horizontal, ce retard est appelé retard " statique ". L'effet de la statique sur l'analyse tomographique sera maintenant décrit. Squires et collaborateurs ont analysé des données provenant d'expériences puits-à-puits, de profil sismique vertical (PSV) et de PSV inverse (PSVI) et ont conclu que la reconstruction tomographique à partir d'une telle 10 combinaison de données présente un contraste de vitesse latéral important qui n'est pas supporté par les données de réflexion de surface au même emplacement. Les données PSV sont saisies en utilisant une source sismique à la surface ou proche de la surface et plusieurs détecteurs dans des 15 positions de fond, alors que dans un profil PSVI, les
données provenant d'une source de fond sont enregistrées en un emplacement de surface. Squires et collaborateurs en ont conclu que la statique non seulement introduisait des erreurs dans la région échantillonnée par le PSVI, mais 20 qu'elle contaminait aussi d'autres régions du tomogramme.
Leurs résultats indiquent que des erreurs de temps de parcours de l'ordre de 2% peuvent aboutir à des erreurs de vitesse tomographique s'élevant à 7%. Dans les résultats de Squires, les données PSV n'ont pas montré un effet aussi 25 important lorsque des sources enterrées proches de la surface étaient utilisées; on peut présumer qu'il existerait un effet important avec une source sismique de surface. Jenkins et collaborateurs rendent compte d'une étude 30 de contrôle séquentiel étendue effectuée sur une période de 31 mois en Indonésie. Jenkins a utilisé des hydrophones implantés au-dessous de la couche soumise aux intempéries, mais a constaté que les hydrophones ne pouvaient pas être laissés dans le sol pendant plus de quelques jours du fait 35 d'une accumulation de boue sur le boîtier; par conséquent, ils ont utilisé des trous pour émissions sismiques et 1l hydrophones revêtus d'acier. Les effets d'un déplacement à la vapeur d'eau sur un champ peu profond (profondeur d'environ 200 m) peuvent être perçus sur des temps de
réflexion à partir d'intervalles du réservoir alors que des 5 réflexions provenant d'au-dessus du réservoir ne présentent aucun changement dans le contrôle séquentiel.
Ce qui n'a pas été reconnu dans l'art antérieur, c'est que, dans la tomographie différentielle, les erreurs dues aux facteurs décrits cidessus s'annulent. L'utilisation 10 d'une tomographie de réflexion en combinaison avec une tomographie puits-à-puits simulée dans la présente invention est basée sur la reconnaissance du fait que les problèmes décrits ci-dessus (incertitude de l'inversion et effet de la statique) sont notablement réduits dans une 15 tomographie séquentielle différentielle. Ce qui est intéressant dans le contrôle d'un réservoir, ce n'est pas le modèle tomographique réel, mais les modifications différentielles dans le sous-sol lorsqu'une substitution de fluides (et des changements possibles de pression) a lieu. 20 Si la couche proche de la surface ne change pas au cours du contrôle, les effets de la statique s'annulent alors et se présentent sous la forme d'une différence de temps de parcours à travers le réservoir. Cependant, si la couche proche de la surface change au cours du contrôle, le 25 procédé de la présente invention peut alors procéder à une correction pour tenir compte de ces changements comme
décrit ci-après.
Dragoset (brevet des E.U.A. n0 5 946 271), dont le contenu est incorporé totalement ici par référence, indique 30 un système de calibrage destiné à être utilisé dans une tomographie séquentielle. Comme indiqué par Dragoset, un système discret de calibrage sismique normalisé, équipé au minimum, est installé de façon permanente sur une structure de réservoir souterrain qui contient des fluides d'un 35 intérêt économique. Des changements de la teneur en fluides en fonction d'un laps de temps de longue durée peuvent
introduire des changements dans les caractéristiques acoustiques des couches de roche du réservoir, produisant un changement correspondant de la signature sismique enregistrée par le système de calibrage. La différence 5 entre deux signatures sur une époque temporelle choisie définit une signature de calibrage séquentielle. La signature de calibrage séquentielle peut être appliquée pour réduire à une norme commune les résultats de système d'acquisition de données sismiques classiques différents, 10 de millésimes différents et de technologies différentes.
Avec le procédé de la présente invention, aucun capteur de calibrage spécial tel que décrit par Dragoset n'est nécessaire; par contre, des changements de forme d'onde, et tout changement dans le retard statique par 15 rapport au temps, peuvent être mesurés directement par les capteurs de fond a,... an et b1... bk implantés de façon permanente. En utilisant les signaux mesurés à ces capteurs de fond, on peut contrôler des changements du retard statique et de la forme d'onde de source pour chacune des 20 sources de surface 203a, 203b... 203k, et on peut les
utiliser pour corriger les données de réflexion de surface.
Ceci permet d'obtenir une information concernant le réservoir sur une large gamme d'angles et d'améliorer la qualité de la reconstruction tomographique. Les corrections 25 de la source peuvent être appliquées à des données de
sismique-réflexion pour améliorer la résolution de réflecteurs tels que 235, 205 qui encadrent le réservoir.
Des différences de temps de parcours pour des réflexions pour de tels réflecteurs sont alors représentatives de 30 changements de l'intervalle du réservoir entre eux.
Un autre aspect de la présente invention qui a trait en particulier à la présence d'une zone soumise aux intempéries est l'utilisation facultative d'une tomographie d'atténuation sismique. La tomographie classique repose sur 35 la prise de temps de parcours pour une propagation d'énergie à partir de plusieurs sources jusqu'à plusieurs récepteurs. Une information supplémentaire concernant le milieu est également contenue dans le spectre de fréquence du signal. Quan et collaborateurs décrivent l'utilisation d'une tomographie d'atténuation sismique. Dans la plupart 5 des matières naturelles, l'atténuation sismique augmente avec la fréquence. Les composantes à haute fréquence du signal sismique sont atténuées plus rapidement que les composantes à basse fréquence pendant la propagation des ondes. Il en résulte que le centrode du spectre du signal 10 subit un décalage vers le bas pendant la propagation. En
supposant un facteur de qualité indépendant de la fréquence (Q), cet abaissement est proportionnel à une intégrale de trajectoire par la distribution d'atténuation et peut être utilisé sous forme de données observées pour reconstruire 15 tomographiquement la distribution de l'atténuation.
L'atténuation est associée à la distribution des fluides dans un réservoir en production et procure donc un autre paramètre qui peut fournir une information sur le réservoir. La présente invention utilise facultativement 20 une tomographie d'atténuation pour fournir un autre aperçu
de l'écoulement des fluides dans le réservoir.
Le procédé de la présente invention peut être mis en oeuvre de façon égale pour des applications terrestres et marines. Pour des applications marines, les sources 25 sismiques sont habituellement des canons à air déployés à partir d'un navire de reconnaissance et les récepteurs peuvent être portés sur des câbles de fltes sismiques ou
peuvent être placés sur le fond de l'océan.
La présente invention a été décrite ci-dessus en 30 utilisant l'exemple de deux puits. En réalité, des récepteurs peuvent être installés de façon permanente dans plusieurs puits sur une zone. Une telle configuration est décrite dans le brevet des E.U.A. n0 5 886 255 de Aronstam et collaborateurs, dont le contenu est incorporé 35 complètement ici par référence. Avec un tel agencement, le procédé décrit ci-dessus peut être utilisé pour obtenir une information tomographique le long de n'importe quelle ligne qui est colinéaire avec deux puits quelconques dans la zone. Ceci est montré sur la figure 5 o plusieurs puits 301a... 301b,... 301n sont représentés. Comme cela est connu, 5 le nombre de combinaisons possibles d'une paire quelconque de n puits est donné par n(n-21)/2. Les exemples de ces lignes sont désignés par 311a, 311b... 311p. Une partie importante du réservoir peut donc être échantillonnée dans
trois dimensions.
La description ci-dessus a fait référence
principalement à des changements des fluides contenus dans un réservoir. Cependant, il est bien connu dans la technique que des vitesses sismiques dépendent également de la contrainte effective dans la formation. La contrainte 15 effective est définie comme étant la différence entre la contrainte de surcharge et la pression du fluide; par conséquent, lorsque du fluide est produit à partir du réservoir, ceci s'accompagne de changements des vitesses sismiques. De plus, la récupération de fluides produit 20 aussi d'autres effets tels qu'un tassement et un effondrement du réservoir, le tassement étant particulièrement notable dans des réservoirs crayeux. Un exemple particulièrement sévère de tassement est le champ Ekofisk en Mer du Nord o le fond de la mer s'est effondré 25 de plus de 3 mètres; ceci a exigé un relevage coteux des plateformes de production et le redéploiement de canalisations. Le procédé de la présente invention peut donc également être utilisé pour identifier des changements de pression et un tassement d'un réservoir, et en 30 particulier pour identifier un tassement avant qu'il
s'étende totalement jusqu'au fond de la mer.
Bien que la description précédente porte sur les
formes préférées de réalisation de l'invention, diverses modifications apparaîtront aux spécialistes de la 35 technique.

Claims (15)

REVENDICATIONS
1. Procédé de détermination d'un paramètre auquel on s'intéresse d'une partie d'un réservoir dans des formations terrestres, comprenant: (a) le déploiement d'une première pluralité de récepteurs sismiques dans un premier puits et d'une seconde pluralité de récepteurs dans un second puits, ledit premier puits et ledit second puits étant disposés sur des côtés opposés de ladite partie du réservoir; (b) la génération d'ondes sismiques en une pluralité d'emplacements à la surface ou proches de la surface et la propagation desdites ondes sismiques dans lesdites formations terrestres à au moins une époque; (c) la réception de premiers signaux représentatifs desdites ondes sismiques de propagation auxdites première et seconde pluralités de récepteurs sismiques à ladite, au moins une, époque; (d) la réception de seconds signaux représentatifs d'ondes sismiques réfléchies depuis une 20 limite en dessous de ladite partie du réservoir en une troisième pluralité de récepteurs à ladite surface ou proche de ladite surface; (e) l'exécution d'une analyse tomographique desdits premiers et seconds signaux pour obtenir ledit 25 paramètre auquel on s'intéresse à ladite, au moins une, époque.
2. Procédé selon la revendication 1, dans lequel ladite pluralité d'emplacements comporte en outre: (A) des emplacements sur des côtés dudit premier 30 puits et dudit second puits opposés à ladite partie du réservoir, et
(B) des emplacements entre ledit premier puits et ledit second puits.
3. Procédé selon la revendication 1, dans lequel ladite, au moins une, époque comprend en outre plusieurs époques.
4. Procédé selon la revendication 2, dans lequel ladite, au moins une, époque comprend en outre plusieurs époques.
5. Procédé selon la revendication 1, comprenant en 5 outre la réception de troisièmes signaux représentatifs d'ondes sismiques réfléchies à partir d'une limite audessus de ladite partie du réservoir en une quatrième pluralité de récepteurs à ladite surface ou proches de ladite surface, et dans lequel l'exécution de ladite 10 analyse tomographique comprend en outre l'utilisation
desdits troisièmes signaux.
6. Procédé selon la revendication 1, dans lequel l'exécution de ladite analyse tomographique comprend en outre l'utilisation desdits premiers signaux pour simuler 15 une reconnaissance puits-à-puits de ladite partie du réservoir.
7. Procédé selon la revendication 2, comprenant en outre, avant l'exécution de ladite analyse tomographique, la correction d'au moins un sous-ensemble desdits deuxièmes 20 signaux en utilisant au moins un sous-ensemble desdits
premiers signaux.
8. Procédé selon la revendication 3, comprenant en outre l'exécution de ladite analyse tomographique à chacune de ladite pluralité d'époques et la mise en relation de 25 différences dans lesdites analyses tomographiques avec au
moins l'un (i) d'un changement d'une distribution de fluides dans ladite partie du réservoir, (ii) d'un changement de pression de fluide dans ladite partie du réservoir, et (iii) d'un effondrement d'une couche dans 30 ladite partie du réservoir.
9. Procédé selon la revendication 4, comprenant en outre l'exécution de ladite analyse tomographique à chacune de ladite pluralité d'époques et la mise en relation de différences dans lesdites analyses tomographiques avec au 35 moins l'un (i) d'un changement d'une distribution de fluides dans ladite partie du réservoir, (ii) d'un changement de pression de fluide dans ladite partie du réservoir, et (iii) d'un effondrement d'une couche dans
ladite partie du réservoir.
10. Procédé selon la revendication 5, dans lequel 5 l'exécution de ladite analyse tomographique comprend en outre la détermination d'une différence de temps de parcours entre lesdits deuxièmes signaux et lesdits
troisièmes signaux.
11. Procédé selon la revendication 6, dans lequel la 10 simulation de ladite reconnaissance puits-à-puits comprend
en outre l'utilisation du principe de Fermat.
12. Procédé selon la revendication 1, dans lequel l'exécution de ladite analyse tomographique comprend en outre l'analyse desdits premiers, deuxièmes et troisièmes 15 signaux pour leur atténuation, ladite atténuation étant
représentative d'au moins l'un (i) d'un changement de la teneur en fluides de ladite partie du réservoir, (ii) d'un changement de pression de ladite partie du réservoir, et (iii) d'un effondrement d'une couche dans ladite partie du 20 réservoir.
13. Procédé selon la revendication 8, dans lequel l'exécution de ladite analyse tomographique comprend en outre l'analyse desdits premiers, deuxièmes et troisièmes signaux pour leur atténuation, ladite atténuation étant 25 représentative d'au moins l'un (i) d'un changement dans une distribution de fluide dans ladite partie du réservoir, (ii) d'un changement de pression de fluide dans ladite partie du réservoir, et (iii) d'un effondrement d'une
couche dans ladite partie du réservoir.
14. Procédé selon la revendication 9, dans lequel
l'exécution de ladite analyse tomographique comprend en outre l'analyse desdits premiers, deuxièmes et troisièmes signaux pour leur atténuation, ladite atténuation étant représentative d'un changement de la teneur en fluide de 35 ladite partie du réservoir.
15. Procédé selon la revendication 1, comprenant en outre le déploiement de récepteurs sismiques en au moins un puits supplémentaire et l'obtention dudit paramètre auquel on s'intéresse sur un volume en trois dimensions.
FR0311610A 2002-10-04 2003-10-03 Procede de determination d'un parametre auquel on s'interesse d'une partie d'un reservoir dans des formations terrestres Expired - Fee Related FR2845773B1 (fr)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US10/264,565 US6894949B2 (en) 2002-10-04 2002-10-04 Walkaway tomographic monitoring

Publications (2)

Publication Number Publication Date
FR2845773A1 true FR2845773A1 (fr) 2004-04-16
FR2845773B1 FR2845773B1 (fr) 2012-10-12

Family

ID=29420152

Family Applications (1)

Application Number Title Priority Date Filing Date
FR0311610A Expired - Fee Related FR2845773B1 (fr) 2002-10-04 2003-10-03 Procede de determination d'un parametre auquel on s'interesse d'une partie d'un reservoir dans des formations terrestres

Country Status (4)

Country Link
US (1) US6894949B2 (fr)
FR (1) FR2845773B1 (fr)
GB (1) GB2394052B (fr)
NO (1) NO335988B1 (fr)

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2397884B (en) * 2003-01-31 2006-02-01 Westerngeco Seismic Holdings A seismic surveying arrangement
WO2005040858A1 (fr) * 2003-10-24 2005-05-06 Shell Internationale Research Maatschappij B.V. Releve sismique periodique d'une region reservoir
US7046581B2 (en) * 2003-12-01 2006-05-16 Shell Oil Company Well-to-well tomography
US7359284B2 (en) * 2004-02-06 2008-04-15 Brian Nelson Fuller Method for processing borehole seismic data
US7751979B2 (en) * 2004-06-30 2010-07-06 Shell Oil Company Determining in-situ the relation between seismic velocity and state of stress in an underground formation
US7531791B2 (en) * 2005-02-17 2009-05-12 Advanced Applied Physics Solutions, Inc. Geological tomography using cosmic rays
WO2007084149A2 (fr) * 2005-02-17 2007-07-26 Triumf, Operating As A Joint Venture By The Governors Of The University Of Alberta, The University Of British Columbia, Carleton University, The University Of Toronto, And The University Of Victoria Tomographie geologique utilisant des rayons cosmiques
US7911878B2 (en) * 2006-05-03 2011-03-22 Baker Hughes Incorporated Sub-salt reflection tomography and imaging by walkaway VSP survey
US7480204B2 (en) * 2006-07-07 2009-01-20 Westerngeco L.L.C. Seismic data processing
US7379385B2 (en) * 2006-07-26 2008-05-27 Westerngeco L.L.C. Processing of seismic data acquired using over/under streamers and/or over/under sources
US8107317B2 (en) 2006-12-28 2012-01-31 Schlumberger Technology Corporation Technique and system for performing a cross well survey
US8547789B2 (en) * 2007-12-27 2013-10-01 Schlumberger Technology Corporation Method for gas zone detection using sonic wave attributes
US20100118655A1 (en) * 2008-11-08 2010-05-13 Andres Chavarria Progressive 3D vertical seismic profiling method
US8750074B2 (en) * 2008-11-25 2014-06-10 Baker Hughes Incorporated Anisotropic parameter determination
US8576663B2 (en) * 2010-04-30 2013-11-05 Schlumberger Technology Corporation Multicomponent seismic inversion of VSP data
CN102884447B (zh) 2010-05-05 2015-08-19 埃克森美孚上游研究公司 Q层析成像方法
FR2984398B1 (fr) * 2011-12-20 2014-01-03 Total Sa Procede de surveillance d'un site sous-marin
US9429670B2 (en) * 2012-09-28 2016-08-30 Cgg Services Sa Plural-depth buried seismic sensors acquisition system and method
US10209388B2 (en) * 2013-12-20 2019-02-19 Schlumberger Technology Corporation Method and apparatus to generate a crosswell data set
US9939537B2 (en) * 2015-03-10 2018-04-10 Decision Sciences International Corporation Sensor fusion with muon detector arrays to augment tomographic imaging using ambient cosmic rays
CN107870361B (zh) * 2017-10-12 2019-09-03 南方科技大学 一种地震回折波层析成像方法、装置及终端设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4969130A (en) * 1989-09-29 1990-11-06 Scientific Software Intercomp, Inc. System for monitoring the changes in fluid content of a petroleum reservoir
GB2275338A (en) * 1993-02-19 1994-08-24 Exxon Production Research Co Crosswell seismic data simulation using Fermat's principle
US5742560A (en) * 1994-05-12 1998-04-21 Exxon Production Research Company Seismic imaging using wave equation extrapolation
US20020128777A1 (en) * 1998-12-30 2002-09-12 Baker Hughes, Inc. Reservoir monitoring in a laminated reservoir using 4-D time lapse data and multicomponent induction data

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4943950A (en) * 1989-05-26 1990-07-24 Western Atlas International, Inc. Method for migrating seismic data
US5740125A (en) 1996-08-30 1998-04-14 Western Atlas International, Inc. Cross-well connectivity mapping including separation of compressional and shear wave energy
FR2759172B1 (fr) 1997-02-05 1999-03-05 Inst Francais Du Petrole Methode de traitement de donnees sismiques de puits multi-composantes orientees
US6147929A (en) 1998-03-03 2000-11-14 Gas Research Institute Method for predicting continuous and discontinuous waveguide targets using interwell seismic signature characteristics
US6289284B1 (en) * 1998-06-30 2001-09-11 Yamamoto Engineering Corporation Method of imaging the permeability and fluid content structure within sediment
US6388947B1 (en) * 1998-09-14 2002-05-14 Tomoseis, Inc. Multi-crosswell profile 3D imaging and method
US6324478B1 (en) * 1999-05-10 2001-11-27 3D Geo Development, Inc. Second-and higher-order traveltimes for seismic imaging
US6169959B1 (en) * 1999-05-27 2001-01-02 Baker Hughes Incorporated Method of predicting kinematics for surface multiples

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4969130A (en) * 1989-09-29 1990-11-06 Scientific Software Intercomp, Inc. System for monitoring the changes in fluid content of a petroleum reservoir
GB2275338A (en) * 1993-02-19 1994-08-24 Exxon Production Research Co Crosswell seismic data simulation using Fermat's principle
US5742560A (en) * 1994-05-12 1998-04-21 Exxon Production Research Company Seismic imaging using wave equation extrapolation
US20020128777A1 (en) * 1998-12-30 2002-09-12 Baker Hughes, Inc. Reservoir monitoring in a laminated reservoir using 4-D time lapse data and multicomponent induction data

Also Published As

Publication number Publication date
US20040068376A1 (en) 2004-04-08
GB2394052A (en) 2004-04-14
US6894949B2 (en) 2005-05-17
GB0323347D0 (en) 2003-11-05
GB2394052B (en) 2004-08-25
NO20034464L (no) 2004-04-05
NO20034464D0 (no) 2003-10-06
NO335988B1 (no) 2015-04-13
FR2845773B1 (fr) 2012-10-12

Similar Documents

Publication Publication Date Title
FR2845773A1 (fr) Procede de determination d'un parametre auquel on s'interesse d'une partie d'un reservoir dans des formations terrestres
EP2369369B1 (fr) Procédé de surveillance d'un site de stockage géologique de gaz par inversion stratigraphique de données sismiques
EP0591037B1 (fr) Méthode d'acquisition et de traitement de données pour suivre le déplacement des fluides dans un réservoir
Lavergne et al. Inversion of seismograms and pseudo velocity logs
CA2670384C (fr) Methode d'inversion conjointe de donnees sismiques representees sur des echelles de temps differentes
EP0122824B1 (fr) Procédé et dispositif pour l'optimisation des données sismiques
FR2872296A1 (fr) Procede destine a amiliorer la resolution sismique
FR2989788A1 (fr) Titre non renseigne.
FR2955396A1 (fr) Dispositif de traitement de donnees sismiques marines
EP2253970A1 (fr) Méthode pour imager une zone cible du sous-sol à partir de données de type walkaway
FR2955397A1 (fr) Procede et dispositif d'acquisition de donnees sismiques marines
FR2831961A1 (fr) Methode de traitement de donnees sismiques de puits en amplitude preservee absolue
FR2739457A1 (fr) Methode de selection de la direction d'un canal de forage par rapport aux plans de fractures geologiques, et son application a la prospection petroliere
FR2844059A1 (fr) Procede de diagnostic de gisements de petrole et de gaz
EP2135113A2 (fr) Traitement de données sismiques acquises en utilisant des sur/sous-flûtes mixtes
FR2916540A1 (fr) Procede d'exploration sismique permettant la supression de fantomes dus aux reflexions a la surface de l'eau, et procede de traitement de donnees sismiques pour la supression de ces fantomes
FR2918178A1 (fr) Methode pour ajuster un modele de vitesse d'ondes sismiques en fonction d'informations relevees aux puits
FR2779532A1 (fr) Procede d'attenuation des reverberations de colonne d'eau dans un signal sismique de capteurs doubles
CA2464799A1 (fr) Methode pour determiner un modele de vitesse d'ondes sismiques dans une formation souterraine heterogene
FR2985831A1 (fr) Traitement de donnees sismiques comprenant une modelisation de multiples internes sans informations de sous-surface
EP0944847B1 (fr) Procede de traitement sismique et notamment procede de prospection sismique 3d mettant en oeuvre une migration des donnees sismiques
Hopkins et al. Maturing DAS VSP as an Onshore CCUS monitoring technology at the Quest CCS Facility
FR2528491A1 (fr) Procede et appareil de retablissement des signaux dans un systeme de diagraphie pendant le forage
FR2544870A1 (fr) Procede pour l'interpretation d'enregistrements sismiques pour donner des caracteristiques exploitables telles que le potentiel gazeifere et la lithologie de couches geologiques
Vesnaver et al. A workflow for processing mono‐channel Chirp and Boomer surveys

Legal Events

Date Code Title Description
PLFP Fee payment

Year of fee payment: 13

ST Notification of lapse

Effective date: 20170630