FR2753800A1 - Procede et appareil pour exploiter des paires de donnees collectees au cours d'une operation sismique, lesquelles representent des caracteristiques d'une formation terrestre et pour generer un enregistrement de sortie - Google Patents
Procede et appareil pour exploiter des paires de donnees collectees au cours d'une operation sismique, lesquelles representent des caracteristiques d'une formation terrestre et pour generer un enregistrement de sortie Download PDFInfo
- Publication number
- FR2753800A1 FR2753800A1 FR9711698A FR9711698A FR2753800A1 FR 2753800 A1 FR2753800 A1 FR 2753800A1 FR 9711698 A FR9711698 A FR 9711698A FR 9711698 A FR9711698 A FR 9711698A FR 2753800 A1 FR2753800 A1 FR 2753800A1
- Authority
- FR
- France
- Prior art keywords
- data
- seismic
- output
- volume
- information
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000012549 training Methods 0.000 title claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 61
- 238000012545 processing Methods 0.000 claims abstract description 44
- 229930195733 hydrocarbon Natural products 0.000 claims abstract description 10
- 150000002430 hydrocarbons Chemical class 0.000 claims abstract description 10
- 230000009467 reduction Effects 0.000 claims description 21
- 230000004044 response Effects 0.000 claims description 10
- 239000004215 Carbon black (E152) Substances 0.000 claims description 3
- 238000012550 audit Methods 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 13
- 230000009977 dual effect Effects 0.000 abstract description 6
- 238000010276 construction Methods 0.000 description 49
- 238000004458 analytical method Methods 0.000 description 13
- 239000002360 explosive Substances 0.000 description 10
- 230000009466 transformation Effects 0.000 description 7
- 230000001419 dependent effect Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 239000011435 rock Substances 0.000 description 6
- 239000012530 fluid Substances 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000013508 migration Methods 0.000 description 4
- 230000005012 migration Effects 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 238000003672 processing method Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000003708 edge detection Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 239000003921 oil Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000013500 data storage Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 235000015076 Shorea robusta Nutrition 0.000 description 1
- 244000166071 Shorea robusta Species 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 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/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Un système de traitement mémorise un logiciel spécial connu sous le nom de logiciel Cubemath. Le logiciel Cubemath permet à un système de traitement de traiter, d'analyser et d'interpréter un second ensemble de données sismiques par rapport à un premier ensemble de données sismiques obtenu soit par une opération sismique répétée soit par une opération sismique unique avec double procédé de traitement. Un premier ensemble de données sismiques 3D est collecté en un emplacement particulier de la surface de la terre à l'instant "t1" et un second ensemble de données sismiques 3D est collecté dans le même emplacement particulier à la surface de la terre (à un instant différent "t2" quand l'opération sismique répétée est exécutée). Le traitement dudit second ensemble de données sismiques par rapport audit premier ensemble de données sismiques est effectué par les opérations suivantes: sélectionner soit un volume soit une surface pour représenter le premier ensemble de données sismiques, sélectionner soit un volume, soit une surface, soit un tracé, soit rien, pour représenter le second ensemble de données sismiques, sélectionner une opération mathématique spécifique, sélectionner un type spécifique de sortie (telle que sortie de volume ou de surface), exécuter l'opération mathématique sur le premier ensemble de données sismiques, et engendrer un ensemble de résultats de données sur un moyen d'enregistrement de sortie présenté à un opérateur sous la forme du type sélectionné de sortie. L'ensemble de données de résultats illustre les changements qui peuvent se produire à l'intérieur de la formation terrestre lorsque le premier et le second ensemble de données sismiques 3D sont chacun collectés. Ces changements dans la formation terrestre seront représentatifs des caractéristiques de la formation terrestre et de la possibilité de l'existence de dépôts, en sous-sol, d'hydrocarbure dans la formation terrestre.
Description
PROCEDE ET APPAREIL POUR EXPLOITER DES PAIRES DE
DONNEES COLLECTEES AU COURS D'UNE OPERATION SISMIQUE,
LESQUELLES REPRESENTENT DES CARACTERISTIQUES D'UNE
FORMATION TERRESTRE ET POUR GENERER UN ENREGISTREMENT
DE SORTIE AFIN DE DETERMINER LES CARACTERISTIQUES DE LA
FORMATION
ART ANTERIEUR
L'objet de la présente invention concerne un procédé et un appareil d'interprétation de données sismiques, et, plus particulièrement, un procédé et un appareil pour exploiter des paires de données obtenues au cours d'une opération sismique répétée ou bien pour exploiter des paires d'attributs déduits d'une seule opération sismique effectuée à la surface de la terre pour lesquelles les paires de données sont représentatives des caractéristiques de la formation terrestre et pour générer un enregistrement de sortie sur un support à partir duquel les caractéristiques de
la formation terrestre peuvent être déterminées.
Des données sismiques sont obtenues en mesurant et en enregistrement les données pendant une "opération sismique 3D" effectuée en utilisant un ensemble de
géophones disposés sur la surface de la terre.
De façon générale, une "opération sismique" est effectuée en faisant détoner une source d'énergie explosive à la surface de la terre et, en utilisant un ensemble de géophones, en mesurant un ensemble de vibrations d'ondes sonores émergeant à la surface de la terre mais décalées d'une pluralité de distances (x, x + 50 pieds, x + 100 pieds, etc.) de l'emplacement de la source d'énergie explosive. Cet ensemble de vibrations sonores émergeantes produites pendant "l'opération sismique" peut être appelé "données sismiques 2D" car, dans l'exemple considéré, les vibrations sonores se
propagent dans le plan x-z.
Cependant, en réalité, les vibrations sonores se propagent dans trois dimensions, c'est-à-dire que les vibrations se propagent le long des axes, x, y et z. En conséquence, les données sismiques émergeantes, lesquelles émergent de la surface de la terre près de l'ensemble des géophones, devraient être correctement
désignées par données sismiques "3D".
Pendant une opération sismique, les vibrations sonores venant de la source d'énergie explosive se propagent vers le bas à l'intérieur de la terre, sont réfléchies sur différentes couches terrestres (des couches de roc ou des schistes sableux), et se propagent vers le haut en direction de la surface de la terre. Il en résulte que l'intensité de chaque vibration sonore qui est réfléchie sur une couche terrestre est fonction de l'impédance de la couche terrestre à l'intérieur de la terre. Par conséquent, les données sismiques 3D enregistrées par les géophones à la surface de la terre représentent un ensemble de caractéristiques de la formation terrestre située à l'intérieur de la terre et placée au-dessous de
l'ensemble des géophones.
On rappelle que les données sismiques "3D" ci-
dessus sont générées et enregistrées par les géophones lorsque l'opération sismique 3D est effectuée à la surface de la terre, les données sismiques 3D représentant un ensemble de vibrations sonores se propageant dans trois dimensions à l'intérieur de la
terre (c'est-à-dire le long des axes x, y et z).
Cependant, il existe une quatrième dimension qui devrait être prise en compte également, ce qui conduirait à la génération et à l'enregistrement, au moyen de l'ensemble des géophones, de données sismiques "4D". Cette quatrième dimension est appelée "temps" ou "t". Il en résulte que les quatre dimensions devraient être (x, y, z et t) o x est la distance horizontale, y la distance transversale, z l'instant de réflexion et t
le temps écoulé.
Les données sismiques 4D sont générées et enregistrées par des géophones lorsqu'une "opération sismique répétée" est effectuée à la surface de la terre. Une "opération sismique répétée" est effectuée de la manière suivante: on effectue une opération sismique 3D en un emplacement de la surface de la terre à l'instant "tl", et puis à un instant ultérieur "t2" on effectue une autre opération sismique 3D au même
emplacement de la surface de la terre à l'instant "t2".
Si les couches terrestres situées à l'intérieur de la terre sous l'ensemble des géophones ont changé d'une quelconque manière (par exemple nouvelle présence de pétrole, de gaz ou de contenu fluide), les données sismiques 3D émergeant de la terre à l'instant "t2" sont différentes des données sismiques 3D qui ont
émergé de la terre à l'instant "tl".
Il résulte de l'exécution de l'opération "sismique répétée" que deux (2) ensembles de données sismiques sont mesurés par l'ensemble des géophones: un premier ensemble de données sismiques mesuré par les géophones à l'instant "tl" et un second ensemble de données
sismiques mesuré par les géophones à l'instant "t2".
Il existe un procédé alternatif pratiqué pour
générer les deux ensembles de données sismiques.
Pendant une seule opération sismique 3D, un unique ensemble de données sismiques est mesuré par l'ensemble
des géophones et déduit de l'opération sismique 3D.
Cependant, cet unique ensemble de données sismiques est traité ultérieurement de deux manières différentes. Il résulte des deux manières différentes de traitement de l'unique ensemble de données sismiques que l'on obtient une paire d'attributs, c'est-à-dire que les deux (2) ensembles de données sismiques comprenant un premier ensemble de données sismiques et un second ensemble de données sismiques sont obtenus. Ce procédé alternatif pour générer les deux ensembles de données sismiques est appelé ci-après "opération sismique unique avec
deux procédés de traitement".
En plus, des ensembles multiples de données peuvent être enregistrés en même temps, mais en mesurant différents types de propagation d'ondes, par exemple des ondes de compression et des ondes de cisaillement. On peut alors considérer qu'il s'agit d'un unique ensemble de données traité de deux
manières.
Cependant, quand le premier et le second ensemble de données sismiques représentatifs des caractéristiques d'une formation terrestre sont mesurés et enregistrés ou obtenus autrement, un système de traitement situé dans un poste de travail devrait analyser ou interpréter le premier et le second ensemble de données sismiques pour déterminer les caractéristiques de la formation terrestre située sous l'ensemble des géophones installés à la surface de la terre. Le système de traitement doit comprendre un progiciel spécial (désigné ci-après par "logiciel Cubemath"), lequel est spécialement étudié pour analyser ou interpréter le premier et le second ensemble de données sismiques qui ont été obtenus soit par "l'opération sismique répétée" soit par
"l'opération unique avec double méthode de traitement".
Cependant, jusqu'à maintenant, il n'existait pas de progiciel spécial de ce type. En conséquence, il existe un besoin pour un progiciel spécial (connu sous le nom de "logiciel Cubemath") lequel, lorsqu'il est mémorisé dans un système de traitement de poste de travail, permet au système de traitement d'analyser et d'interpréter le premier et le second ensemble de données sismiques obtenus soit d'une opération sismique répétée soit d'une opération sismique unique avec double procédé de traitement dans le but de déterminer les caractéristiques de la formation terrestre située sous un ensemble de géophones installés à la surface de la terre.
EXPOSE DE L'INVENTION
Aussi, c'est un premier objet de la présente invention de proposer un système de traitement et un progiciel spécial connu sous le nom de logiciel Cubemath à utiliser en association avec ce système de traitement pour analyser et interpréter un second ensemble de données sismiques par rapport à un premier ensemble de données sismiques obtenu au cours d'une ou plusieurs opérations sismiques et représentatif d'un ensemble de caractéristiques d'une formation terrestre, dans le but de déterminer les caractéristiques de la formation terrestre et l'éventualité de l'existence de dépôts souterrains d'hydrocarbures à l'intérieur de
cette formation.
C'est encore un objet de la présente invention de proposer un système de traitement et un progiciel spécial connu sous le nom de logiciel Cubemath à utiliser en association avec ce système de traitement pour traiter, analyser et interpréter un second ensemble de données sismiques par rapport à un premier ensemble de données sismiques obtenu au cours d'une ou plusieurs opérations sismiques et représentatif d'un ensemble de caractéristiques d'une formation terrestre, dans le but de déterminer les caractéristiques de la formation terrestre et l'éventualité de l'existence de dépôts souterrains d'hydrocarbures dans cette formation, le traitement du second ensemble de données sismiques par rapport au premier ensemble de données sismiques étant effectué par une exploitation mathématique des premier et second ensembles de données sismiques pour produire un résultat, le résultat étant indicatif des caractéristiques de la formation terrestre. C'est encore un objet de la présente invention de proposer un système de traitement et un progiciel spécial connu sous le nom de logiciel Cubemath à utiliser en association avec ce système de traitement pour traiter, analyser et interpréter un second ensemble de données sismiques par rapport à un premier ensemble de données sismiques obtenu au cours d'une ou plusieurs opérations sismiques et représentatif d'un ensemble de caractéristiques d'une formation terrestre, dans le but de déterminer les caractéristiques de la formation terrestre et l'éventualité de l'existence de dépôts souterrains d'hydrocarbures dans cette formation, le traitement du second ensemble de données sismiques par rapport au premier ensemble de données sismiques étant effectué en sélectionnant soit un volume soit une surface pour représenter le premier ensemble de données sismiques, en sélectionnant soit un volume soit une surface soit un tracé pour représenter le second ensemble de données sismiques, et en exploitant mathématiquement le volume ou la surface du premier ensemble de données sismiques et le volume ou la surface ou le tracé du second ensemble de données sismiques, grâce à quoi un résultat est produit, le résultat étant indicatif d'éventuels changements intervenus dans la formation terrestre et par conséquent des caractéristiques de la formation terrestre. C'est encore un objet de la présente invention de proposer un système de traitement et un progiciel spécial connu sous le nom de logiciel Cubemath à utiliser en association avec ce système de traitement pour traiter, analyser et interpréter un second ensemble de données sismiques par rapport à un premier ensemble de données sismiques obtenu au cours d'une ou plusieurs opérations sismiques, le traitement du second ensemble de données sismiques par rapport au premier ensemble de données sismiques étant effectué en sélectionnant soit un volume soit une surface pour représenter le premier ensemble de données sismiques, en sélectionnant soit un volume soit une surface soit un tracé pour représenter le second ensemble de données sismiques, en sélectionnant une opération mathématique, en sélectionnant un type spécifique de sortie (par exemple une sortie de volume ou une sortie de surface), en exécutant l'opération mathématique ci-dessus sur le premier. ensemble de données sismiques représenté par un volume ou une surface et sur le second ensemble de données sismiques représenté par un volume ou une surface ou un tracé, et en générant un résultat présenté sous la forme du type de sortie spécifique sélectionné, le résultat ci-dessus étant indicatif d'éventuels changements intervenus dans la formation terrestre et par conséquent indicatif des
caractéristiques de la formation terrestre.
Selon les objets ci-dessus et d'autres objets de la présente invention, un système de traitement mémorise un logiciel spécial appelé "logiciel Cubemath". Le logiciel Cubemath permet au système de traitement de traiter, d'analyser et d'interpréter un second ensemble de données sismiques par rapport à un premier ensemble de données sismiques qui ont été obtenus soit par une opération sismique répétée soit par une opération sismique unique avec double traitement. L'opération sismique répétée: (1) collecte un premier ensemble de données sismiques 3D en un emplacement particulier de la surface de la terre à l'instant "tl", et (2) collecte un second ensemble de données sismiques 3D au même emplacement particulier de
la surface de la terre à l'instant "t2".
L'opération sismique unique avec double traitement collecte un seul ensemble de données sismiques 3D à l'instant tl, mais, pendant l'unique opération sismique avec double traitement, cet ensemble unique de données sismiques 3D collecté à l'instant tl est traité de deux manières différentes, grâce à quoi des paires d'attributs dérivés d'une seule opération sismique sont obtenues, c'est-à-dire: (1) un premier ensemble de données sismiques 3D, et (2) un second ensemble de
données sismiques 3D.
Le traitement du second ensemble de données sismiques 3D par rapport au premier ensemble de données sismiques 3D est effectué en sélectionnant soit un volume soit une surface pour représenter le premier ensemble de données sismiques 3D, en sélectionnant soit un volume soit une surface soit un tracé pour représenter le second ensemble de données sismiques 3D, en sélectionnant une opération mathématique spécifique, en sélectionnant un type spécifique de sortie (par exemple sortie de volume ou sortie de surface), en exécutant un programme permettant d'effectuer l'opération mathématique spécifique ci-dessus sur le volume ou la surface du premier ensemble de données sismiques et sur le volume ou la surface ou le tracé du second ensemble de données sismiques, et en générant un résultat sur un moyen d'enregistrement de sortie, par exemple un écran cathodique ou une sortie sur imprimante d'ordinateur ou un autre support d'affichage, en réponse à l'exécution de l'opération mathématique spécifique. Ce résultat est présenté à un utilisateur ou à un opérateur sous la forme du type de sortie spécifique qui a été sélectionné d'avance. Le résultat ci-dessus, présenté sous la forme du type de sortie spécifique sélectionné, illustre de quelconques changements qui se sont produits dans la formation terrestre: (1) pour l'opération sismique répétée, entre l'instant "tl" et l'instant "t2", lorsque les premier et second ensembles de données sismiques 3D ont été chacun collectés, et (2) pour l'opération sismique unique avec double procédé de traitement, à l'instant tl, lorsque le premier et le second ensemble de données sismiques 3D ont été déduits d'un seul ensemble de données sismiques obtenu pendant l'unique opération sismique. Ces changements de la formation terrestre sont indicatifs des caractéristiques de la formation terrestre et de l'éventualité de l'existence de dépôts souterrains d'hydrocarbures (c'est-à-dire du pétrole)
dans la formation terrestre.
En référence aux Figures 8, 9, 11, 15, 16, 21 et 28 et en utilisant un système de traitement basé dans un poste de travail comprenant le logiciel Cubemath de la présente invention, un opérateur installé au poste de travail est capable d'effectuer les sélections suivantes: (1) l'opérateur sélectionne des données primaires 56 en cliquant sur le premier cube 34, (2) l'opérateur sélectionne le type de données des données primaires sélectionnées telles que "le volume 62", (3) l'opérateur sélectionne des données secondaires 70a sur la Figure 9, (4) l'opérateur sélectionne le type de données des données secondaires telles que le volume 70c, (5) l'opérateur sélectionne une opération mathématique 70jl, et, dans l'exemple du mode de réalisation préféré de la Figure 11, l'opération mathématique sélectionnée 70j4 est V3 = V1/V2, (6) l'opérateur sélectionne une sortie 80a sur la Figure 15, (7) l'opérateur sélectionne un type de données b pour la sortie sélectionnée, comme par exemple le volume 80c, (8) l'opérateur définit un nouveau volume sismique 80e2 sur la Figure 21, et (9) l'opérateur sélectionne l'icône "exécution" 42 sur la Figure 28, ce qui exécute le code "exécution" a de la Figure 16 et le programme est exécuté pour générer la "sortie 80a" et le nouveau volume sismique e2. Un géologue peut utiliser la "sortie 80a" pour effectuer: (1) Une interprétation lithologique - C'est un exemple de l'utilisation mentionnée ci-dessus "une seule opération sismique avec deux procédés de traitement" - Il existe des techniques pour analyser des données sismiques avant traitement par empilement du point de vue des effets du décalage en fonction de l'amplitude (AVO). Deux volumes de sortie créés pour ce procédé sont l'empilement des ondes de pression (P) et le volume de gradient (G). La combinaison sous forme d'un produit (P x G) de ces deux grandeurs peut être
utilisée par un géologue comme indication lithologique.
La somme (P + G) est liée au rapport de Poisson, une propriété élastique de la roche. La différence (P - G) est une estimation de la réponse de l'onde de cisaillement.
(2) Transformation temps en profondeur -
L'analyse des données sismiques avant traitement par empilement et des informations de forage peut donner un volume de vitesse moyenne qui peut être combiné aux informations sismiques pour la transformation en profondeur. Le procédé implique des opérations multiples, à savoir: multiplier la vitesse moyenne pour chaque échelon de temps du volume sismique et l'indice d'échelon de temps pour obtenir un indice de profondeur (V x T), puis interpoler l'amplitude sismique pour l'indice de profondeur désiré
(interpolation sin x/x).
(3) Analyse des données sismiques 4D (répétition des relevés sismiques pour la surveillance des réservoirs) - Une analyse du contraste entre deux relevés sismiques acquis au même emplacement peut indiquer la migration des fluides des réservoirs. Un simple tracé par soustraction de tracés pourrait servir à faire ressortir ces différences si l'acquisition des deux ensembles de données a été enregistrée de manière identique. D'une façon générale, il existe de légères différences dont on peut tenir compte par des opérations supplémentaires pour ajuster les différences de temps et/ou d'amplitude entre les ensembles de
données, c'est-à-dire par corrélation croisée.
(4) Variations d'échantillons - Plusieurs opérations dans un seul volume, en volume/surface, en volume/tracé, peuvent être effectuées pour des applications telles la détection du bord, le tranchage d'un volume, les tests d'hypothèses et le traitement par empilement à posteriori. Ces opérations exigent une combinaison de fonctions mathématiques de base effectuées selon la présente invention dans le présent
descriptif.
D'autres domaines d'application de la présente
invention apparaîtront d'après la description détaillée
présentée ci-après. Cependant, il faut comprendre que
la description détaillée et les exemples spécifiques,
tout en représentant un mode de réalisation préféré de la présente invention, sont donnés à titre d'illustration uniquement, puisque différents changements et modifications entrant dans le cadre et l'esprit de l'invention paraîtront évidents à un spécialiste de la technique à la lecture de la
description détaillée qui suit.
BREVE DESCRIPTION DES FIGURES
L'invention peut être complètement comprise
d'après la description détaillée d'un mode de
réalisation préféré présenté ci-dessous et d'après les dessins annexés qui sont présentés à titre d'illustration uniquement et ne visent en rien à limiter la présente invention, et sur lesquels: les Figures 1 et 2 illustrent le procédé et l'appareil nécessaires pour exécuter une opération sismique 3D; la Figure 3 illustre un schéma simplifié du système de traitement installé dans le camion d'enregistrement utilisé pendant l'opération sismique 3D pour collecter les premier et second ensembles de données sismiques 3D obtenus pendant l'opération sismique 3D des Figures 1 et 2 et pour mémoriser les premier et second ensembles de données sismiques 3D sur un support d'enregistrement de sortie tel qu'une bande magnétique; la Figure 4 illustre un schéma simplifié d'un ordinateur central qui utilise un logiciel de réduction des données mémorisées pour effectuer une réduction des données sur le premier et le second ensemble de données sismiques 3D mémorisés sur le support d'enregistrement de sortie de la Figure 3 et produire un enregistrement de sortie des données réduites sur un support; la Figure 5 illustre un poste de travail et un écran associé pour recevoir le support d'enregistrement des données de sortie réduites; la Figure 6 illustre plus en détail le poste de travail de la Figure 5 représentant la mémoire du poste de travail qui mémorise le logiciel Cubemath permettant au système de traitement de la présente invention de traiter, analyser et interpréter un second ensemble de données sismiques par rapport à un premier ensemble de données sismiques obtenus à partir d'une opération sismique répétée; les Figures 7a et 7b illustrent un organigramme du logiciel de réduction des données de la Figure 4; les Figures 8 à 22 illustrent des schémas d'état du logiciel Cubemath de la Figure 6, lequel permet au système de traitement de la Figure 6 de la présente invention de traiter, analyser et interpréter un second ensemble de données sismiques par rapport à un premier ensemble de données sismiques obtenus au cours d'une opération sismique répétée, les Figures 8 à 22 comprenant: la Figure 8 illustrant l'état initial et comprenant un volume primaire disponible, une surface primaire disponible, une sélection de volume et une sélection de surface; la Figure 9 illustrant une construction détaillée du volume primaire disponible de la Figure 8, la Figure 10 illustrant une construction détaillée de la surface primaire disponible de la Figure 8, la Figure 17 illustrant une construction détaillée de la sélection de volume de la Figure 8, la Figure 18 illustrant une construction détaillée de la sélection de surface de la Figure 8, la Figure 9 illustrant un volume primaire disponible et comprenant une sélection de tracé, un volume/tracé prêt, une sélection de surface, un volume/surface prêt, une sélection de volume et un volume/volume disponible, la Figure 19 illustrant une construction détaillée de la sélection de tracé de la Figure 9, la Figure 13 illustrant une construction détaillée du volume/tracé disponible de la Figure 9, la Figure 18 illustrant une construction détaillée de la sélection de surface de la Figure 9, la Figure 12 illustrant une construction détaillée du volume/surface disponible de la Figure 9, la Figure 17 illustrant une construction détaillée de la sélection de volume de la Figure 9, et la Figure 11 illustrant une construction détaillée du volume/volume disponible de la Figure 9, la Figure 10 illustrant une surface primaire disponible et comprenant une sélection de surface et une surface/surface disponible, la Figure 18 illustrant une construction détaillée de la sélection de surface de la Figure 10, la Figure 14 illustrant une construction détaillée de la surface/surface disponible de la Figure 10, la Figure 11 illustrant un volume/volume disponible et comprenant un contrôle de paramètre et un ensemble opérateur, la Figure 20 illustrant une construction détaillée de contrôle des paramètres de la Figure 11, la Figure 15 illustrant une construction détaillée de l'ensemble opérateur de la Figure 11, la Figure 12 illustrant un volume/surface prêt et comprenant un contrôle des paramètres et un ensemble opérateur, la Figure 20 illustrant une construction détaillée du contrôle des paramètres de la Figure 12, la Figure 15 illustrant une construction détaillée de l'ensemble opérateur de la Figure 12, la Figure 13 illustrant un volume/tracé disponible et comprenant un contrôle des paramètres et un ensemble opérateur, la Figure 20 illustrant une construction détaillée du contrôle des paramètres de la Figure 13, la Figure 15 illustrant une construction détaillée de l'ensemble opérateur de la Figure 13, la Figure 14 illustrant une surface/surface disponible et comprenant un contrôle des paramètres et un ensemble opérateur, la Figure 20 illustrant une construction détaillée du contrôle des paramètres de la Figure 14, la Figure 15 illustrant une construction détaillée de l'ensemble opérateur de la Figure 14, la Figure 15 illustrant un ensemble opérateur et comprenant une définition du volume de sortie, une définition de la surface de sortie et une exécution prêt, la Figure 21 illustrant une construction détaillée de la définition du volume de sortie de la Figure 15, la Figure 22 illustrant une construction détaillée de la définition de la sortie de surface de la Figure , la Figure 16 illustrant une construction détaillée de l'exécution disponible de la Figure 15, les Figures 23 à 32 illustrent les affichages
individuels présentés à l'opérateur sur l'écran (c'est-
à-dire le tube cathodique) du poste de travail des Figures 5 et 6, les Figures 23 à 32 comprenant: la Figure 23 illustrant la fenêtre de l'état initial telle qu'elle est indiquée par le bloc d'état initial de la Figure 8, la Figure 24 illustrant la fenêtre de sélection des données telle qu'elle est indiquée par le bloc de type des données des Figures 8 et 9, la Figure 25 illustrant la fenêtre disponible du volume primaire tel qu'elle est indiquée par le bloc disponible de volume primaire de la Figure 9, la Figure 26 illustrant la fenêtre de l'ensemble opérateur telle qu'elle est indiquée par le bloc des opérations disponibles conduisant au bloc de l'ensemble opérateur des Figures 11, 12, 13 et 14, la Figure 27 illustrant la fenêtre de sortie des données comme indiqué par le bloc du type de données de la Figure 15, la Figure 28 illustrant la fenêtre d'exécution disponible comme indiqué par le bloc d'exécution disponible de la Figure 16, la Figure 29 illustrant la fenêtre de sélection de volume comme indiqué par les blocs de sélection de volume et de volume sélectionné de la Figure 17, la Figure 30 illustrant la fenêtre de sortie de volume comme indiqué par les blocs de volume et de sortie de volume définis de la Figure 21, la Figure 31 illustrant la fenêtre de surface de sortie comme indiqué par les blocs de surface et de sortie de surface sélectionnés de la Figure 22, la Figure 32 illustrant la fenêtre de contrôle des paramètres comme indiqué par le bloc de contrôle des paramètres de la Figure 20; et les Figures 33 à 35 illustrent des exemples du type spécifique de sortie (par exemple sortie de volume ou sortie de surface) qui ont été sélectionnés par l'opérateur.
DESCRIPTION DETAILLEE DE MODES DE REALISATION
D'après les Figures 1 et 2, un appareil et un procédé associé pour exécuter une opération sismique 3D en un emplacement de la surface de la terre sont illustrés. Sur la Figure 1, une source d'énergie explosive 10 située sous la surface de la terre 12 détone et génère une pluralité de vibrations sonores qui se propagent vers le bas et sont réfléchies sur une couche 16 à l'intérieur de la terre. La couche 16 pourrait être une couche de roche ou une couche de sable ou une couche de schiste. Quand les vibrations sonores sont réfléchies par la couche 16 de la terre, lesvibrations sonores 14 se propagent vers le haut et sont reçues dans une pluralité de récepteurs 18 appelés géophones 18. La pluralité des géophones 18 génèrent chacun un signal électrique en réponse à la réception d'une vibration sonore à l'intérieur et une pluralité de signaux électriques sont générés par les géophones 18, la pluralité des signaux étant reçus dans un camion d'enregistrement 20. La pluralité des signaux électriques provenant des géophones 18 représentent un ensemble de caractéristiques de la formation terrestre située dans la terre sous les géophones 18, et, en particulier, les caractéristiques de la portion de la terre située adjacente à la couche 16 dans la terre. Le camion 20 contient un ordinateur 20a qui reçoit et mémorise la pluralité des signaux reçus des géophones 18. Un support d'enregistrement de sortie est généré à partir de l'ordinateur 20a dans le camion d'enregistrement 20, lequel comprend et/ou affiche et/ou mémorise la pluralité des signaux électriques qui sont représentatifs des caractéristiques de la formation terrestre située dans la terre sous les
géophones 18. L'appareil et le procédé décrits ci-
dessus en référence à la Figure 1 sont appelés opération sismique 2D (qui signifie bidimensionnelle) car le procédé indiqué ci-dessus est destiné à générer des vibrations sonores 14 le long des axes x et z. Mais, en réalité, les vibrations sonores 14 devraient se propager le long des axes x, y et z. Il en résulte que l'appareil et le procédé décrits ci-dessus en référence à la Figure 1 devraient être désignés de manière plus appropriée par opération sismique "3D" (puisque les vibrations sonores 14 se propagent le long des axes x, y et z de la figure 1). L'axe x représente la distance horizontale, l'axe y représente la distance transversale et l'axe z représente l'instant de réflexion. Sur la Figure 2, un autre procédé et un autre appareil pour exécuter une opération sismique 3D sont représentés. La Figure 2 a été empruntée à un ouvrage intitulé "Seismic Velocity Analysis and the Convolutional Model", par Enders A. Robinson, dont le
texte est inclus par référence au présent descriptif.
Sur la Figure 2, l'opération sismique 3D de la Figure 1 est exécutée 10 fois différentes. Par exemple, lorsque la source d'énergie explosive 10 est disposée en une position 22 (la première position ou position "l'0" le long de la surface de la terre) de la Figure 2, une première pluralité de signaux électriques provenant des géophones 18 est mémorisée dans l'ordinateur 20a du camion d'enregistrement 20. La source d'énergie explosive est déplacée jusqu'à la position 24. Quand la source d'énergie explosive 10 est disposée dans la position 24 (seconde position ou position "1" le long de la surface de la terre), une seconde pluralité de signaux électriques sont mémorisés dans l'ordinateur a du camion d'enregistrement 20. La source d'énergie explosive 10 est déplacée de façon répétitive et séquentielle des positions "2" à "9" de la Figure 2 jusqu'à ce qu'elle soit disposée dans la position 26 (c'est-à-dire la position "9" qui est la dixième position) à la surface de la terre. Lorsque la source d'énergie explosive 10 est disposée dans la position 26 (la dixième position le long de la surface de la terre), une dixième pluralité de signaux électriques sont mémorisés dans l'ordinateur 20a du camion d'enregistrement 20. Il en résulte, comme indiqué sur la Figure 2, que le camion d'enregistrement 20 enregistre 10 tracés (dix ensembles de signaux électriques, chaque ensemble comprenant une pluralité de signaux électriques) entre la position 22 et la position 26 le long de la surface de la terre. Un support d'enregistrement de sortie est généré par l'ordinateur 20a dans le camion d'enregistrement 20, lequel comprend les 10 tracés ou les dix ensembles de signaux électriques reçus des géophones 18. Le procédé et l'appareil décrits ci-dessus en référence aux Figures 1 et 2 représentent une "opération sismique 3D". En référence à la Figure 3, une construction plus détaillée de l'ordinateur du camion d'enregistrement 20a est représentée. L'ordinateur 20a du camion d'enregistrement comprend un processeur 20al et une mémoire 20a2 connectés à un bus système. Les dix tracés ou les dix ensembles de signaux électriques reçus des géophones 18 des Figures 1 et 2 et générés par les géophones 18 pendant l'opération sismique 3D devraient être reçus dans l'ordinateur du camion d'enregistrement a par le bloc 20a3 de "données reçues" de la Figure 3 connecté au bus système et mémorisés dans la mémoire
a2 de l'ordinateur 20a du camion d'enregistrement.
Lorsqu'on le souhaite, un support d'enregistrement de sortie 20a4, également relié au bus système, est généré, lequel contient, mémorise et/ou affiche les dix tracés ou les dix ensembles de signaux électriques
reçus des géophones 18.
L'exposé ci-dessus décrit une opération sismique tridimensionnelle (3D). Cependant, une technique sismique pratiquée sur place, appelée opération sismique à 4 dimensions ou 4D, implique l'emploi de quatre dimensions, à savoir, l'axe x, l'axe y et l'axe z pour la propagation des vibrations sonores 14 dans les trois dimensions de la terre, et une quatrième dimension appelée "temps" ou l"t"l. La quatrième dimension, à savoir le temps, représente le temps écoulé et sert à effectuer une opération sismique relativement nouvelle appelée ci-après "opération
sismique répétée" ou "opération sismique 4D".
"L'opération sismique répétée" ou "opération sismique 4D" est effectuée de la manière suivante en référence à nouveau aux Figures 1 et 2: une première opération sismique 3D est effectuée à l'instant "tl" en un emplacement particulier à la surface de la terre, l'opération sismique 3D ayant été exposée ci-dessus en référence aux Figures 1 et 2. Ensuite, lorsqu'une période de temps prédéterminée, par exemple 6 mois, s'est écoulée depuis le temps "tl" o la première opération sismique 3D a été effectuée une seconde opération sismique 3D est effectuée à l'instant "t2" au même emplacement particulier sur la surface de la terre. On peut rappeler qu'un (1) ensemble de données comprenant les dix tracés ou dix ensembles de signaux électriques reçus des géophones 18 est obtenu quand une opération sismique 3D est effectuée, donc quand une "opération sismique 4D" ou "opération sismique répétée" est terminée, deux (2) ensembles de données sismiques 3D sont reçus et mémorisés dans l'ordinateur 20a du camion d'enregistrement, de la manière suivante: (1) un premier ensemble de dix tracés ou de dix ensembles de signaux électriques (désigné ci-après par "un premier ensemble de données sismiques 3D") et (2) un second ensemble de dix tracés ou de dix ensembles de signaux électriques (désigné ci-après par "un second ensemble de données sismiques 3D"). Lorsque le support d'enregistrement de sortie 20a4 de la Figure 3 est généré par l'ordinateur du camion d'enregistrement 20a, le "premier ensemble de données sismiques 3D" est mémorisé et/ou affiché sur un premier support d'enregistrement de sortie 20a4; et, six mois plus tard, le "second ensemble de données sismiques 3D" est mémorisé et/ou affiché sur un second support d'enregistrement de sortie 20a4. Il y a donc alors deux (2) supports d'enregistrement de sortie 20a4: un premier qui mémorise le premier ensemble de données sismiques 3D et un second qui mémorise le second ensemble de données sismiques 3D mesuré environ 6 mois après la mesure du premier ensemble de données
sismiques 3D.
Rappel de "l'opération sismique unique avec deux procédés de traitement": une seule opération sismique 3D est effectuée et un seul ensemble de données sismiques 3D est obtenu; cet unique ensemble de données sismiques 3D est traité par deux (2) procédés différents pour produire ainsi deux (2) ensembles distincts (ou paires d'attributs) de données sismiques 3D: un premier ensemble de données sismiques 3D et un second ensemble de données sismiques 3D. Deux supports d'enregistrement de sortie 20a4 sont préparés: un premier qui mémorise le premier ensemble de données sismiques 3D et un second qui mémorise le second ensemble de données sismiques 3D. En référence à la Figure 4, un schéma simplifié de l'ordinateur central est représenté, lequel utilise un "logiciel de réduction des données" mémorisé pour exécuter une opération de "réduction des données" sur les premier et second ensembles de données sismiques 3D mémorisés sur les premier et second supports d'enregistrement de sortie 20a4, respectivement, de la Figure 3, l'ordinateur principal produisant un premier support d'enregistrement de sortie des données réduites et un second support d'enregistrement de sortie des données réduites, dont chacun affiche et/ou mémorise des versions réduites des premier et second ensembles
de données sismiques 3D.
Sur la Figure 4, un ordinateur central 30 comprend un processeur central 30a connecté à un bus système et une mémoire 30b également connectée au bus système, laquelle mémorise un "logiciel de réduction des données". Les premier et second supports d'enregistrement de sortie 20a4 de la Figure 3 sont connectés au bus système de la Figure 4, et le premier ensemble de données sismiques 3D ainsi que le second ensemble de données sismiques 3D provenant du premier et du second supports d'enregistrement de sortie 20a4 sont chacun mis à la disposition du processeur central 30a et de son logiciel associé de réduction des données mémorisé dans la mémoire 30b. Le processeur central 30a exécute le logiciel de réduction des données mémorisé dans la mémoire 30a et, comme résultat de l'exécution du logiciel de réduction des données, le processeur principal 30a effectue une opération de "réduction des données" sur le premier et le second ensemble de données sismiques 3D mémorisés sur le premier et le second support d'enregistrement de sortie 20a4, génère un "premier support d'enregistrement de sortie de données réduites" 30d associé au premier ensemble de données sismiques 3D mémorisé sur le premier support d'enregistrement de sortie 20a4, et génère un "second support d'enregistrement de sortie de données réduites" d associé au second ensemble de données sismiques 3D mémorisé sur le second support d'enregistrement de sortie 20a4. Donc, sur la Figure 4, lorsque le logiciel de réduction des données contenu dans la mémoire 30b est exécuté par le processeur central 30a en association avec le premier et le second ensemble de données sismiques 3D, deux supports d'enregistrement de sortie de données réduites 30d sont générés: un premier 30d qui mémorise une version des données réduites du premier ensemble des données sismiques 3D et un second 30d qui mémorise une version des données
réduites du second ensemble de données sismiques 3D.
Les versions de données réduites des premières et secondes données sismiques 3D mémorisées sur les supports d'enregistrement de sortie de données réduites d représentent des images corrigées dans l'espace du sous-sol de la formation terrestre de la Figure 1. La technique de réduction des données décrite ci-dessous en référence aux Figures 7a et 7b implique la somme vectorielle des réflexions 14 sur la couche 16 de la
terre représentée sur la Figure 1.
En référence aux Figures 7a et 7b, un organigramme du logiciel de réduction des données mémorisé dans la mémoire 30b de l'ordinateur central 30 de la Figure 4 est représenté. L'organigramme du logiciel de réduction des données des Figures 7a et 7b est emprunté à un ouvrage intitulé "Seismic Velocity Analysis and the Convolutional Model", par Enders A. Robinson, dont le texte a déjà été incorporé en référence au présent descriptif. Sur les Figures 7a et 7b, l'organigramme du logiciel de réduction des données comprend les blocs suivants: un bloc de démultiplexage 30bl connecté à l'entrée, un bloc de tri 30b2, un bloc d'élimination de gain 30b3, un bloc de filtre de fréquence 30b4, un bloc de ré-échantillonnage 3ObS, un bloc de sélection de tracé 30b6, une sortie 30b7 repérée "groupements sélectionnés" (analyses de vitesse), un bloc de correction d'amplitude 30b8, un bloc de déconvolution b9, une seconde sortie 30blO repérée "tracés triés CMP après déconvolution", un bloc de corrections de temps 30bll, un bloc AGC 30b12, un bloc d'empilage 30b13, une troisième sortie 30b14 repérée "tracés empilés (non filtrés)", un bloc de filtrage de fréquence 30b15, un autre bloc AGC 30b!6, une quatrième sortie 30b17 repérée "tracés empilés (filtrés)", une seconde entrée repérée "informations d'inclinaison" 30b18, un bloc d'interpolation de tracés 30b19, un bloc de migration 30b20, une cinquième sortie 30b21 repérée "tracés migrés (non filtrés)", un bloc de filtrage de fréquences 30b22, un bloc AGC 30b23, une sixième sortie b24 repérée "tracés migrés (filtrés)", un bloc de correction temps-profondeur 30b25, et une septième sortie 30b26 repérée "tracés migrés (migrés en profondeur)". Dans l'organigramme des Figures 7a et 7b, l'une quelconque des sorties 30b7, 30b10, 30b14, 30b17, b21, 30b24 et 30b26 peut sertir d'entrées au poste de travail d'interprétation 32 expliqué ci-dessous et
illustré sur la Figure 5 des dessins.
D'après la Figure 5, un poste de travail d'interprétation 32 selon la présente invention, comportant un écran 32a, est électriquement connecté aux, et reçoit les, deux premier et second supports d'enregistrement de sortie de données réduites 30d (par exemple une bande magnétique) qui ont été générés par l'ordinateur principal 30 de la Figure 4. Les versions à données réduites du premier ensemble de données sismiques 3D et du second ensemble de données sismiques 3D (reçues en utilisant "l'opération sismique répétée" ou "l'opération sismique unique avec deux procédés de traitement") mémorisées sur le premier support d'enregistrement de sortie des données réduites 30d et sur le second support d'enregistrement de sortie des données réduites 30d sont chargées dans le poste de
travail 32 de la Figure 5.
La Figure 6 représente plus en détail le poste de travail d'interprétation 32 de la présente invention
présenté sur la Figure 5.
Sur la Figure 6, le poste de travail 32 comprend l'écran 32a et le processeur de poste de travail 32b connecté au bus système, ainsi qu'une mémoire 32c aussi connectée électriquement au bus système. Le premier et le second support d'enregistrement de sortie des données réduites 30d sont connectés au bus; il en résulte que le premier et le second ensemble de données sismiques 3D mémorisés sur le premier et le second support d'enregistrement de sortie des données réduites d sont accessibles par le processeur 32b du poste de travail et son logiciel associé mémorisé dans la mémoire 32c du poste de travail. La mémoire 32c du poste de travail mémorise deux types de logiciels: un logiciel d'interprétation IESX 32cl, et un logiciel Cubemath 32c2 selon la présente invention. Lorsque le processeur 32b du poste de travail applique le logiciel Cubemath 32c2 en association avec le logiciel IESX 32cl, la version à données réduites du premier ensemble de données sismiques 3D (obtenu à l'instant tl en utilisant l'opération sismique répétée) peut être traitée mathématiquement par rapport à la version à données réduites du second ensemble de données sismiques 3D (obtenu à l'instant t2 en utilisant l'opération sismique répétée) en produisant ainsi une "sortie particulière", à partir de laquelle le géologue peut déterminer les caractéristiques de la couche de formation 16 de la Figure 1. Par exemple, le géologue peut être capable de déterminer, à partir de la "sortie particulière", quelle catégorie de roche existe dans la formation, en changeant la teneur en fluide dans la formation ou la profondeur de la couche 16 de la formation. Ce descriptif donne des opérations supplémentaires ci- dessous concernant d'autres conclusions possibles que le géologue peut tirer de la "sortie particulière" générée par le processeur 32b du poste de travail de la présente invention quand le processeur 32b du poste de travail applique le logiciel
Cubemath 32c2 de la présente invention.
Le poste de travail 32 peut être formé d'un poste de travail Silicon Graphics Indigo2. Le système d'exploitation est IRIX 5.3. Le programme logiciel peut être écrit en langage de programmation C selon les normes Unix et Motif. Le programme peut faire l'objet d'une nouvelle compilation et être exécuté sur des postes de travail Sun en liaison avec d'autres produits IESX énumérés ci-dessous, lesquels sont disponibles auprès de Geoquest, une division de la Schlumberger Technology Corporation, Houston, Texas. En plus de l'environnement d'exploitation du poste de travail Unix, le logiciel minimal d'interprétation IESX 32cl nécessaire pour appliquer le logiciel Cubemath 32c2 de la présente invention est le suivant (un tel logiciel d'interprétation IESX 32cl étant disponible auprès de Geoquest, une division de la Schlumberger Technology Corporation, Houston, Texas):
1. licence d'utilisation IESX, partie N UAMR1-
QD1
2. gestionnaire de données IESX, partie N UAMD1-
QD1
3. IESX Seis3DV, partie N UA3Dl-QDl
Une brève description de l'utilisation
fonctionnelle de la présente invention est présentée dans les deux paragraphes suivants en référence aux
Figures 23 et 33-34 des dessins.
Sur la Figure 23, lorsque le logiciel Cubemath 32c2 (et le logiciel d'interprétation IESX 32cl) sont appliqués par le processeur 32b du poste de travail de la Figure 6, la fenêtre d'état initial 33 de la Figure 23 est affichée sur l'écran 32a de la Figure 6. Dans la fenêtre d'état initial 33, un premier cube 34 représente le premier ensemble à données réduites des données sismiques 3D mesuré par les géophones 18 de la Figure 1 à l'instant "tl" (quand l'opération sismique répétée a été effectuée), qui ultérieurement subit le traitement de réduction des données de la Figure 4. Un second cube 36 représente le second ensemble de données sismiques 3D mesuré par les géophones 18 de la Figure 1 à l'instant "t2" (environ 6 mois après l'instant tl auquel l'opération sismique répétée est exécutée), lequel subit ultérieurement le traitement de réduction des données de la Figure 4. L'opérateur sélectionne une opération mathématique 38 sur la Figure 23. L'opérateur sélectionne ensuite une sortie particulière 40. Puis l'opérateur sélectionne "exécution" 42. Quand l'exécution 42 est choisie, le premier cube de données 34 agit mathématiquement sur le second cube de données 36 en utilisant l'opération mathématique 38, et puis la sortie particulière 40 est générée. Le géologue utilise la sortie particulière 40 pour déterminer les caractéristiques de la formation terrestre dans le but final de déterminer l'existence de dépôts souterrains
d'hydrocarbures dans la formation.
Sur les Figures 33-34, des exemples de la "sortie particulière" 40 sont illustrés. Les Figures 33-34 représentent l'opération mathématique Vo = V1 * V1, la Figure 33 représentant le volume d'entrée V1 et la Figure 34 le volume de sortie Vo. Sur la Figure 33, le même volume V1 est utilisé comme volumes sismiques primaire et secondaire d'entrée et, sur la Figure 34, un volume d'auto-corrélation de sortie Vo est présenté pour une analyse géophysique. Des utilisations possibles de la "sortie particulière" 40 comprennent des opérations mathématiques effectuées sur des volumes sismiques pour combiner, contraster et transformer les entrées originales. Quelques exemples seront discutés plus en détail ci-dessous, dont une interprétation lithologique, une transformation temps-profondeur, une analyse d'études sismiques 4D ou d'études sismiques répétées pour la surveillance des réservoirs, et des
variations d'échantillons.
En référence aux Figures 8-35, plusieurs diagrammes d'état sont illustrés sur les Figures 8-22 représentant le logiciel Cubemath 32c2 de la Figure 6 de la présente invention, plusieurs présentations de fenêtre sont illustrées sur les Figures 23-32, lesquelles sont formées sur l'écran 32a des Figures 5 et 6 quand le logiciel Cubemath 32c2 des Figures 8-22 est exécuté, et deux exemples de la "sortie particulière" ci-dessus sont illustrés sur les Figures 33-35 pour être utilisés par le géologue afin de déterminer les caractéristiques d'une formation terrestre d'un point de vue géologique quand
l'exécution du logiciel Cubemath 32c2 est terminée.
Sur la Figure 8, le premier d'une série de diagrammes d'état représentant le logiciel Cubemath 32c2 de la présente invention, est illustré lors de l'exécution du système de traitement de la Figure 6 de la présente invention. Sur la Figure 8, un état initial est initialisé au moment o l'exécution du logiciel Cubemath 32c2 a commencé. L'utilisateur/opérateur peut prendre congé en 52 et sortie en 54 du programme à ce
moment si cela est nécessaire.
Sur la Figure 23, pendant l'état initial 50, la fenêtre d'état initial 33, représentative de l'état initial 50, apparaît sur l'écran 32a du poste de travail indiqué sur les Figures 5 et 6. Cette
présentation de fenêtre a été expliquée ci-dessus.
Comme cela a déjà été indiqué, l'opérateur assis au poste de travail 32 de la Figure 6 sélectionne "exécution" 42. Quand l'exécution 42 est sélectionnée, le premier cube de données 34, représentant les versions à données réduites des données sismiques 3D obtenues à l'instant "tl" est traité mathématiquement sur le second cube de données 36, représentant les versions à données réduites des données sismiques 3D obtenues à l'instant "t2" (lorsque l'opération sismique répétée est effectuée) en utilisant l'opération mathématique 38, et, après l'opération mathématique, la sortie particulière 40 est générée. La sortie particulière est employée par le géologue comme indiqué ci-dessus et fait l'objet d'une discussion plus approfondie ci-dessous. Si "l'opération sismique unique avec double traitement" est exécutée, le premier cube 34 représente la version à données réduites du premier ensemble de données sismiques (déduite de l'unique opération sismique originale à partir de laquelle une paire d'attributs ou d'ensembles de données sismiques sont déduits par un traitement selon deux procédés différents) et le second cube 36 représente la version à données réduites du second ensemble de données sismiques. Sur la Figure 8, l'opérateur sélectionne les données primaires 56, et à ce moment, l'opérateur doit
spécifier le type de données primaires 58.
Sur la Figure 24, la sélection du type de données primaires est illustré. Sur la Figure 24, lorsque le type de données primaires 56 est sélectionné, l'opérateur déplace le curseur jusqu'au premier cube 34 sur lequel il clique sur le premier cube 34. Le type de données primaires a été sélectionné. Pour spécifier le type de données primaires 58, l'opérateur assis au poste de travail 32 doit alors sélectionner ou bien le
volume 1 60a ou bien la surface 1 60b.
Sur la Figure 8, l'opérateur a sélectionné le volume 1 60a, comme indiqué par le numéro d'élément 62 sur la Figure 8, le bloc de sélection de volume 64 de la Figure 8 est exécuté. Une construction détaillée du bloc de sélection de volume 64 est illustrée sur la Figure 17. Quand l'exécution du bloc de sélection de volume 64 est terminée, le bloc 70 de volume primaire disponible est exécuté. Une construction détaillée du bloc de volume primaire disponible est illustrée sur la
Figure 9.
Sur les Figures 17 et 29, en se référant d'abord à la Figure 17, quand les blocs de sélection de volume 64 ou 70g des Figures 8 et 9, respectivement, sont chacun exécutés, le bloc de "sélection de volume" 64g de la Figure 17 et le bloc de "sélection de volume sismique" 64h de la Figure 17 ont déjà été exécutés, et le bloc de "volume sélectionné" 64a de la Figure 17 est exécuté. Lorsque le bloc de "volume sélectionné" 64a de la Figure 17 est exécuté, la fenêtre 65 de sélection de volume de la Figure 29 est présentée à l'opérateur sur l'écran 32a du poste de travail 32. Sur la Figure 17, l'opérateur, en utilisant la fenêtre 65 de sélection de volume de la Figure 29, sélectionne un trajet 64b qui définit la sélection de trajet 64c. La mise en oeuvre de la sélection de trajet 64c indique que l'opérateur doit sélectionner une plage de ligne 64d, une plage de
tracé 64e et une plage de temps 64f.
Sur la Figure 29, la fenêtre 65 de sélection de
volume (intitulée "sélection de volume sismique -
sélection de données sismiques - trajet de sélection") comprend un premier bloc 65à (intitulé "sélection d'étude") dans lequel l'opérateur sélectionne une étude sismique, et un second bloc 65b (intitulé "sélection de catégorie sismique") dans lequel l'opérateur
sélectionne une catégorie sismique.
Sur les Figures 9, 24 et 25, en commençant par la Figure 9, lorsque le bloc de volume primaire disponible de la Figure 8 est exécuté, en référence à la Figure 9, l'opérateur du'poste de travail 32 doit sélectionner des données secondaires 70a et puis préciser un type de données secondaires 70b. Pour sélectionner le type de données secondaires, en référence à la Figure 24, l'écran de fenêtre visible sur la Figure 24 est présenté à l'opérateur. Sur la Figure 24, l'opérateur déplace le curseur sur le second cube 36 et clique sur le second cube 36. Il en résulte, comme indiqué sur la Figure 24, que quatre sélections sont accessibles à
l'opérateur: volume 2, surface 2, tracé ou rien.
L'opérateur au poste de travail 32 regardant l'écran de la Figure 24 clique sur, et sélectionne l'un, des
éléments suivants: volume 2, surface 2, tracé ou rien.
Encore d'après la Figure 9, si l'opérateur sélectionne volume 2, surface 2 ou tracé sur la présentation de fenêtre de la Figure 24, il met en oeuvre et sélectionne respectivement le volume 70c, la surface 70d ou le tracé 70e de la Figure 9. Si l'opérateur sélectionne "rien" sur la Figure 24, il exécute "annulation" 70f sur la Figure 9, et à ce moment le code revient au "volume primaire disponible" 70 et la présentation de fenêtre de la Figure 25 est présentée à l'opérateur. Mais si le volume 70c, la surface 70d ou le tracé 70e de la Figure 9 est sélectionné, le programme exécute ensuite, sur la Figure 9, le bloc de sélection de volume 70g, le bloc de sélection de surface 70h, ou le bloc de sélection de tracé 70I, respectivement. Une construction détaillée du bloc de
sélection de volume 70g est présentée sur la Figure 17.
Une construction détaillée du bloc de sélection de surface 70h est présentée sur la Figure 18. Une construction détaillée du bloc de sélection de tracé I est présentée sur la Figure 19. Si le programme exécute, sur la Figure 9, la sélection de volume 70g, la sélection de surface 70h ou la sélection de tracé 70I, en supposant que le programme ne soit pas abandonné en 70M pour revenir au volume primaire disponible 70, le programme exécute ensuite le volume/volume disponible 70j, le volume/surface disponible 70k ou le volume/tracé disponible 70L, respectivement. La Figure 11 représente une
construction détaillée du volume/volume disponible 70j.
Une construction détaillée du volume/surface disponible k est présentée sur la Figure 12. Une construction détaillée du volume/tracé disponible 70L est présentée sur la Figure 13. La Figure.9 comprend aussi d'autresblocs auxiliaires qui devraient être signalés. Par exemple, sur la Figure 9, les données primaires peuvent être sélectionnées au moyen du bloc "de sélection de données primaires" 70N (analogue au bloc 56 de la Figure 8), le type de données peut être sélectionné au moyen du bloc "type de données" 70P (analogue au bloc 58 de la Figure 8), et puis le programme revient au programme indiqué sur la Figure 8 des dessins. De même, il est possible de sélectionner une opération mathématique au moyen du bloc de "sélection d'opération" 70R (voir "sélection d'opération" sur les Figures 11-14), les opérations disponibles 70S sont présentées en vue de la sélection (voir "opérations disponibles" sur les Figures 11-14), et puis le programme revient au programme indiqué sur les Figures
11-14.
Sur la Figure 19, lorsque le bloc de sélection de tracés 70I de la Figure 9 est exécuté, le bloc de "sélection de tracé" 70I de la Figure 19 implique la sélection d'un volume sismique 7011. Quand le volume sismique est sélectionné en 7012, l'opérateur du poste
de travail 32 de la Figure 6 sélectionne un tracé 7013.
Quand le tracé est sélectionné en 7014, le programme revient sur "état de tracé correct" 7015, ce qui permet au programme d'avancer à l'état appelant. Mais, en différents points pendant l'exécution du programme de sélection de tracé 70I de la Figure 19, le programme peut prendre congé en 7016 pour passer à "l'état appelant" 7017, ce qui ramène le programme à l'état qui
a fait la demande de sélection de volume.
Sur la Figure 8, si l'opérateur a sélectionné surface 1 60b (sur la Figure 24), comme indiqué par le numéro d'élément 66 de la Figure 8, le bloc de sélection de surface 68 est exécuté. Une construction détaillée du bloc de sélection de surface 68 est illustrée sur la Figure 18. Quand l'exécution du bloc de sélection de surface 68 est terminée, le bloc de surface primaire disponible 72 ' est exécuté. Une construction détaillée du bloc de surface primaire
disponible 72 est illustrée sur la Figure 10.
Sur la Figure 18, quand les blocs de sélection de surface 68 ou 70h des Figures 8 ou 9, respectivement, sont chacun exécutés, le bloc de "sélection de surface" 68g et la "sélection de surface" 68h de la Figure 18 ont déjà été exécutés, et le bloc de "surface sélectionnée" 68a de la Figure 18 a été exécuté. Sur la Figure 18, l'opérateur sélectionne un trajet 68b qui met en oeuvre la sélection de trajet 68c. La mise en oeuvre de la sélection de trajet 68c signifie la sélection d'une plage de ligne 68d et la sélection
d'une plage de tracé 68e.
Sur les Figures 10 et 24, quand le bloc de surface primaire disponible 72 de la Figure 8 est exécuté, en référence à la Figure 10, l'opérateur du poste de travail 32 doit sélectionner des données secondaires 72a et puis préciser un type de données secondaires 72b. Sur la Figure 24, pour sélectionner des données secondaires, l'opérateur doit cliquer sur le second cube 36 et puis cliquer sur "surface 2". Sur la Figure , quand l'opérateur sélectionne "surface 2" (d'après la Figure 24), l'opérateur a mis en oeuvre et sélectionné la surface 72d de la Figure 10. Si l'opérateur sélectionne "rien" (d'après la Figure 24), il exécute "annulation" 72f sur la Figure 10, et à ce moment, le programme revient à "surface primaire disponible" 72. Mais si la surface 72d de la Figure 10 est sélectionnée, le programme exécute ensuite, selon la Figure 10, la sélection de surface 72h. Une construction détaillée de la sélection de surface 72h est présentée sur la Figure 18. Si le programme exécute, comme sur la Figure 10, la sélection de surface 72h, en supposant que le programme ne prenne pas congé en 72M pour revenir à la surface primaire disponible 72, le programme exécute ensuite surface/surface disponible 72k. Une construction détaillée de surface/surface disponible 72k est présentée sur la Figure 14. La Figure 10 comprend aussi
d'autres blocs auxiliaires qui devraient être signalés.
Par exemple, sur la Figure 10, les données primaires peuvent être sélectionnées au moyen du bloc "sélection de données primaires" 72N (analogue au bloc 56 de la Figure 8), le type de données peut être sélectionné au moyen du bloc "type de données" 72P (analogue au bloc 58 de la Figure 8), et puis le programme revient au programme représenté sur la Figure 8 des dessins. De même, une opération mathématique peut être sélectionnée au moyen du bloc "sélection d'opération" 72R (voir "sélection d'opération" sur les Figures 11-14), les opérations disponibles 72S sont présentées pour la sélection (voir "opérations disponibles" sur les Figures 11-14), et puis le programme revient au
programme indiqué sur les Figures 11-14.
Les paragraphes ci-dessus ont introduit, sur la Figure 9, le "volume/volume disponible" 70j, le "volume/surface disponible" 70k et le "volume/tracé disponible" 70L, et, sur la Figure 10, la "surface/surface disponible" 72k. Le volume/volume disponible 70j (de la Figure 11), le volume/surface disponible 70k (de la Figure 12), le volume/tracé disponible 70L (de la Figure 13) et la surface/surface disponible 70K (de la Figure 14), sont tous semblables du point de vue construction de programme, et chacun d'eux fera l'objet d'une discussion dans les
20. paragraphes suivants.
D'après la Figure 11, le bloc volume/volume disponible 70j permet à l'opérateur de sélectionner une opération 70jl, et, quand une opération 70jl est sélectionnée, un certain nombre d'opérations disponibles 70j2 sont présentées à l'opérateur. Sur la Figure 11, trois opérations mathématiques au moins sont disponibles pour la sélection: une première opération mathématique 70j3 est la suivante: V3 = V1 * V2. La sélection de la première opération mathématique 70j3 peut impliquer une autre sélection de certains paramètres de corrélation croisée 70j8. Ou encore, une seconde opération mathématique 70j4 est la suivante: V3 = Vl/V2. Encore une autre opération mathématique j5 est la suivante: V3 = Min (V1, V2). Plusieurs autres opérations mathématiques 70j6 sont également à la disposition de l'opérateur pour sélection. Quand les autres opérations mathématiques 70j6 sont sélectionnées par l'opérateur, certains paramètres indépendants et dépendants doivent peut-être être sélectionnés au moyen d'un bloc de commande de paramètres 70j7. Une construction détaillée du bloc de commande de paramètres 70j7 est présentée sur la Figure 20 des dessins. D'autre part, il est possible que les paramètres n'aient pas à être sélectionnés puisqu'il n'y a besoin d'aucun paramètre 70j9. Quand les opérations mathématiques sont sélectionnées (l'un des blocs 70j3, 70j4, 70j5, 70j7 ou 70j9 de la Figure 11 a été exécuté), en supposant que les options "prise de congé" 70j10 et 70jll n'aient pas été sélectionnées (l'option prise de congé 70j10 étant la sortie des paramètres de corrélation croisée 70j8 et l'option de prise de congé 70jll étant la sortie de la commande de paramètres 70j7), le bloc "ensemble opérateur" 80 est exécuté. Une construction détaillée du bloc "ensemble opérateur" 80 est présentée sur la Figure 15 des dessins. Cependant, à partir du bloc volume/volume disponible 70j de la Figure 11, un opérateur peut aussi sélectionner des données primaires 70j12 et déterminer un type de données 70j13 (sélection du premier cube 34 et sélection de volume ou de surface 60a ou 60b sur la Figure 24), en ce point, le programme revient aux Figures 8 des dessins, en 70j14. L'opérateur pourrait sélectionner des données secondaires 70j15 (second cube 36, Figure 24) et déterminer un type de données 70j16 (volume ou surface ou tracé ou rien, Figure 24), et en ce point, le programme revient aux Figures 9 et 10 des dessins, en 70j17. Ou encore, à partir du bloc volume/volume disponible 70j, l'opérateur peut prendre
congé en 70j18.
En référence à la Figure 12, le bloc volume/surface disponible 70k permet à l'opérateur de sélectionner une opération 70kl, et lorsqu'une opération 70kl est sélectionnée, un certain nombre d'opérations disponibles 70k2 sont présentées à l'opérateur. Sur la Figure 12, trois opérations mathématiques au moins sont disponibles pour la sélection: une première opération mathématique 70k3 est la suivante: S3 = V1 x aS2. La sélection de la première opération mathématique 70k3 peut impliquer une autre sélection d'un multiplicateur scalaire 70k8. Ou encore, une seconde opération mathématique 70k4 est la suivante: S3 = V1 n S2. Une autre opération
mathématique 70k5 est la suivante = V3 = V1 I V1 < S2.
Plusieurs autres opérations mathématiques 70k6 sont également disponibles pour une sélection par l'opérateur. Lorsque les autres opérations mathématiques 70k6 sont sélectionnées par l'opérateur, certains paramètres indépendants et dépendants peuvent devoir être sélectionnés au moyen d'un bloc de commande de paramètres 70k7. Une construction détaillée du bloc de commande de paramètres 70k7 est présentée sur la Figure 20 des dessins. D'autre part, les paramètres peuvent ne pas avoir à être sélectionnés puisqu'il n'y a besoin d'aucun paramètre 70k9. Lorsque les opérations mathématiques sont sélectionnées (l'un des blocs 70k3, 70k4, 70k5, 70k7 ou 70k9 de la Figure 12 a été exécuté), en supposant que les opérations "prise de congé" 70klO et 70kll n'aient pas été sélectionnées (l'opération prise de congé 70klO étant la sortie du bloc de multiplicateur scalaire 70k8 et l'option prise de congé 70kll étant la sortie de la commande de paramètres 70k7), le bloc "ensemble opérateur" 80 est exécuté. Une. construction détaillée du bloc "ensemble opérateur" 80 est présentée sur la Figure 15 des dessins. Mais, à partir du bloc volume/surface disponible 70k de la Figure 12, un opérateur peut aussi sélectionner des données primaires 70k12 et déterminer un type de données 70k13 (sélection du premier cube 34 et sélection de volume ou de surface 60a ou 60b sur la Figure 24), en ce point, le programme revient sur la Figure 8 des dessins, en 70k14. L'opérateur pourrait sélectionner des données secondaires 70k15 (second cube 36, Figure 24) et déterminer un type de données 70k16 (volume ou surface ou tracé ou rien, Figure 24), en ce point, le programme revient aux Figures 9 et 10 des dessins, en 70k17. A titre d'alternative, à partir du bloc volume/surface disponible 70k, l'opérateur peut
prendre congé en 70k18.
D'après la Figure 13, le bloc volume/tracé disponible 70L permet à l'opérateur de sélectionner une opération 70L1, et, lorsqu'une opération 70L1 est sélectionnée, un certain nombre d'opérations disponibles 70L2 sont présentées à l'opérateur. Sur la Figure 13, trois opérations mathématiques au moins sont disponibles pour la sélection: une première opération mathématique 70L3 est la suivante: V3 = Vl *-1 T2. La sélection de la première opération mathématique 70L3 peut impliquer l'autre sélection de certains paramètres de convolution 70L8. A titre d'alternative, une seconde opération mathématique 70L4 est la suivante: V3 = V1 * T2. La sélection de la seconde opération mathématique L4 peut impliquer l'autre sélection de paramètres de corrélation 70L19. Une autre opération mathématique encore 70L5 est la suivante: V3 = V1 - T2. Plusieurs autres opérations mathématiques 70L6 sont également à la disposition de l'opérateur pour sélection. Quand les autres opérations mathématiques 70L6 sont sélectionnées par l'opérateur, certains paramètres indépendants ou dépendants peuvent devoir être sélectionnés au moyen d'un bloc de commande de paramètres 70L7. Une construction détaillée du bloc de commande de paramètres 70L7 est présentée sur la Figure 20 des dessins. D'autre part, les paramètres peuvent ne pas avoir à être sélectionnés puisqu'il n'y a besoin d'aucun paramètre 70L9. Lorsque les opérations mathématiques sont sélectionnées (l'un des blocs 7OL3, 70L4, 70L5, 70L7 ou 70L9 de la Figure 13 a été exécuté), en supposant que les options de "prise de congé" 70L10 et 70L11 n'aient pas été sélectionnées (l'option de prise de congé 70L10 étant la sortie des paramètres de convolution 70L8 et des paramètres de corrélation 70L19 et l'option de prise de congé 70L11 étant la sortie de la commande de paramètres 70L7), le bloc "ensemble opérateur" 80 est ensuite exécuté. Une construction détaillée du bloc "ensemble opérateur" 80 est présentée sur la Figure 15 des dessins. Cependant, à partir du bloc volume/tracé disponible 70L de la Figure 13, un opérateur peut aussi sélectionner des données primaires 70L12 et déterminer un type de données 70L13 (sélection du premier cube 34 et sélection de volume ou de surface 60a ou 60b sur la Figure 24), à ce moment, le programme revient aux Figures 8 des dessins, en 70L14. L'opérateur pourrait sélectionner des données secondaires 70L15 (second cube 36, Figure 24) et déterminer un type de données 70L16 (volume ou surface ou tracé ou rien, Figure 24), à ce moment le programme revient aux Figures 9 et 10 des dessins. A titre d'alternative, à partir du bloc volume/tracé disponible 70L, l'opérateur peut prendre congé en 70L18. D'après la Figure 14, le bloc surface/surface disponible 72k permet à l'opérateur de sélectionner une opération 72kl, et lorsqu'une opération 72kl a été sélectionnée, un certain nombre d'opérations disponibles 72k2 sont présentées à l'opérateur. Sur la Figure 14, trois opérations mathématiques au moins sont disponibles pour la sélection: une première opération
mathématique 72k3 est la suivante: S3 = a(S1)/b(S2).
La sélection de la première opération mathématique 72k3 peut impliquer une autre sélection de paramètres scalaires 72k8. A titre d'alternative, une seconde opération mathématique 72k4 est la suivante: S3 = S1 + S2. Une autre opération mathématique encore 72k5 est la suivante: S3 = Min (Sl, S2). Plusieurs autres opérations mathématiques 72k6 sont également disponibles pour la sélection par l'opérateur sélection. Quand les autres opérations mathématiques 72k6 sont sélectionnées par l'opérateur, certains paramètres indépendants et dépendants peuvent avoir à être sélectionnés au moyen d'un bloc de commande de paramètres 72k7. Une construction détaillée du bloc de commande de paramètres 72k7 est présentée sur la Figure des dessins. D'autre part, les paramètres peuvent ne pas avoir à être sélectionnés puisqu'aucun paramètre 72k9 n'est nécessaire. Lorsque les opérations mathématiques sont sélectionnées (l'un des blocs 72k3, 72k4, 72k5, 72k7 ou 72k9 de la Figure 14 a été exécuté), en supposant que les options "prise de congé" 72k10 et 72kll n'aient pas été sélectionnées (l'option de prise de congé 72k10 étant la sortie du bloc des paramètres scalaires 72k8 et l'option de prise de congé 72kll étant la sortie de la commande de paramètres 72k7), le bloc "ensemble opérateur" 80 est alors exécuté. Une construction détaillée du bloc "ensemble opérateur" 80 est présentée sur la Figure 15 des dessins. Cependant, à partir du bloc surface/surface disponible 72k de la Figure 14, un opérateur peut aussi sélectionner des données primaires 72k12 et déterminer un type de données 72k13 (sélection du premier cube 34 et sélection de volume ou de surface 60a ou 60b sur la Figure 24), à ce moment, le programme revient à la Figure 8 des dessins, en 72k14. L'opérateur pourrait sélectionner des données secondaires 72k15 (second cube 36, Figure 24) et déterminer un type de données 72k16 (volume ou surface ou tracé ou rien, Figure 24), en ce point le programme revient aux Figures 9 et 10 des dessins, en 72k17. A titre d'alternative, à partir du bloc surface/surface disponible 72k, l'opérateur peut
prendre congé en 72k18.
Sur les Figures 20 et 32, une construction détaillée du bloc de commande de paramètres 70j7 de la Figure 11, du bloc de commande de paramètres 70k7 de la Figure 12, du bloc de commande de paramètres 70L7 de la Figure 13 et du bloc de commande de paramètres 70K7 de
la Figure 14 est illustrée.
Sur les Figures 20 et 32, en considérant d'abord la Figure 20, l'exécution du bloc de commande de paramètres 100 appelle l'écran de fenêtre 102 présenté sur la Figure 32 des dessins. Le bloc des commande de paramètres 100 de la Figure 20 apparait respectivement sous la forme des blocs de commande de paramètres 70j7,
k7, 70L7, 72K7 sur les Figures 11 à 14.
Sur la Figure 32, l'écran de fenêtre 102, intitulé "Sélection des paramètres pour DtMap" - "commande de corrélation", comprend des commandes opérateur 102a (qui apparaissent sous le titre "commande de corrélation") et une définition de surface 102b, ainsi qu'une définition de fenêtre sismique 102c (lesquelles
apparaissent sous le titre "commande temps fenêtre").
La définition de surface 102b comprend en 102bl "fenêtre de référence temps de début - temps de fin"
et en 102b2 "sélection d'horizon - temps de fin".
Sur la Figure 20, à partir du bloc de commande de paramètres 100, l'opérateur du poste de travail 32 devrait sélectionner un ou plusieurs paramètres 104, à ce moment, des paramètres indépendants ont été sélectionnés en 106. Le programme peut avancer, à ce moment, à l'état d'opération correcte 108, ce qui permet au programme d'avancer à l'état appelant. Mais, l'opérateur peut sélectionner des paramètres dépendants 110; dans ce cas, les paramètres dépendants ont été sélectionnés 112. Le programme peut alors avancer jusqu'à l'état d'opération correcte en 108. Mais l'opérateur du poste de travail 32 peut prendre congé en 114 pour passer à l'état appelant à différents instants pendant l'exécution du programme, l'état appelant revenant à l'état dans lequel la demande de
sélection de volume a été faite.
D'après les Figures 25 et 26, une représentation de fenêtre 82 est illustrée. L'écran de fenêtre 82 apparaît sur l'écran 32a du poste de travail d'interprétation 32 de la Figure 6 pendant l'exécution du programme volume/volume disponible 70jl à 70j18 indiqué sur la Figure 11, le programme volume/surface disponible 70kl à 70k18 sur la Figure 12, le programme volume/tracé disponible 70L1 à 70L18 sur la Figure 13, et le programme surface/surface disponible 70kl à 70k18
sur la Figure 14 des dessins.
Sur la Figure 25, sur l'écran de fenêtre 82, lorsque les données primaires sont sélectionnées, l'opérateur clique sur le premier cube 34. Quand les données secondaires sont sélectionnées, l'opérateur clique sur le second cube 36. Puis l'opérateur sélectionne une opération mathématique en cliquant sur l'icône "sélection d'opération" 38. Quand les données primaires du premier cube 34 subissent une opération mathématique sur les données secondaires du second cube 36 en utilisant l'opération mathématique de l'icône 38, une sortie prédéterminée, sélectionnée, 40 est générée pour être utilisée par le géologue au cours de la détermination des caractéristiques de la formation terrestre. Sur la Figure 26, des exemples de l'opération mathématique sélectionnée en cliquant sur l'icône 38 de l'écran de fenêtre 82 sont illustrés. Sur la Figure 26, lorsque l'opérateur clique sur l'icône 38, au moyen d'une souris 84, plusieurs opérations mathématiques 86 sont présentées à l'opérateur. Au moyen de la sourie 84, l'opérateur clique seulement sur une (1) des opérations mathématiques 86 pour sélectionner une opération mathématique particulière. Par exemple, le numéro d'élément 88 (lequel est pointé sur "*") indiquerait la sélection de l'opération "convolution", le numéro 90 (qui est pointé sur "Max") indiquerait la sélection de l'opération "Maximum" et le numéro 92 (qui est pointé sur "Min") indiquerait la sélection de
l'opération "Minimum".
Sur les Figures 15, 21, 22, 27, 30 et 31, une construction détaillée du programme "ensemble opérateur" 80 auquel il est fait référence sur les
Figures 11 à 14 des dessins est représentée.
Sur la Figure 15, le code de l'ensemble opérateur est utilisé pour sélectionner une sortie particulière, comme une sortie de volume ou une sortie de surface avant d'exécuter l'opération mathématique qui a déjà été sélectionnée sur les Figures 11-14 des dessins. Sur la Figure 15, à partir de "l'ensemble opérateur" 80, l'opérateur installé au poste de travail 32 de la Figure 6 sélectionne une sortie 80a et puis
sélectionne un type de données 80b.
Sur la Figure 27 une représentation de fenêtre 84 est présentée à l'opérateur sur l'écran 32a du poste de travail. L'opérateur sélectionne au moyen de la souris 84a une sortie particulière 84b qui est soit une sortie
de volume 1 84bl soit une sortie de surface 1 84b2.
Sur la Figure 15, l'opérateur sélectionne ou bien la sortie de volume 80c ou bien la sortie de surface d. Si l'opérateur sélectionne la sortie de volume c, le programme de définition de la sortie de volume e est exécuté. Une construction détaillée du programme de définition de la sortie de volume 80e est présentée sur la Figure 21 des dessins. Cependant, si l'opérateur sélectionne la sortie de surface 80d, le programme de définition de sortie de surface 80f est exécuté. Une construction détaillée du programme de définition de la sortie de surface 80f est présentée sur la Figures 22 des dessins. Quand le programme de définition de sortie de volume 80e est exécuté, ou quand le programme de définition de sortie de surface f est exécuté, le programme est alors en 90 à "exécution disponible". Cependant, à partir de "l'ensemble opérateur" 80, l'opérateur installé au poste de travail 32 peut encore sélectionner des données primaires 80g, sélectionner des données secondaires 80h et sélectionner une opération mathématique 80I. Lorsque les données primaires 80g sont sélectionnées, le type de données 80j est sélectionné (Figure 24) et puis le programme revient au programme présenté sur la Figure 8 des dessins. Lorsque les données secondaires 80h sont sélectionnées, le type de données 80k est sélectionné (Figure 24) et puis le programme revient au programme présenté sur les Figures 9-10 des dessins. Quand l'opération mathématique 80I est spécifiée, les opérations mathématiques disponibles 80L sont spécifiées, et puis le programme revient au
programme présenté sur les Figures 11-14 des dessins.
Sur la Figure 21, lorsque le type de données de volume 80c est sélectionné sur la Figure 15, le programme de définition de sortie de volume 80e, représenté sur la Figure 21, génère une sortie de volume 80el sur une représentation de fenêtre 92 de
l'écran 32a du poste de travail 32 de la Figure 6.
L'écran de fenêtre 92 est présenté sur la Figure 30 des dessins. L'opérateur définit un nouveau volume sismique 80e2, et, quand le nouveau volume sismique 80e2 est
défini, le bloc de volume défini 80e3 est exécuté.
Lorsque le bloc de volume défini 80e3 est exécuté, la définition de sortie de volume spécifiée est présentée sur l'écran de fenêtre de sortie de volume 92
apparaissant sur l'écran 32a du poste de travail 32.
L'écran de fenêtre 92 est présenté sur la Figure 30.
* Lorsque le bloc de volume défini 80e3 est exécuté, l'opérateur du poste de travail 32 peut aussi définir une autre sortie en option 80e4, auquel cas, la sortie en option est définie en 80e5. En tous cas, lorsque l'un ou l'autre bloc 80e3 de volume défini est exécuté, ou lorsque le bloc de sortie en option défini 80e5 est exécuté, l'opérateur du poste de travail 32 peut sortir vers "l'état correct de volume de sortie" 80e6, ce qui avance à l'état appelant. Naturellement, l'opérateur du poste de travail 32 peut exécuter la "prise de congé" e7 vers l'état appelant (retour à l'état qui a effectué la demande de sélection de volume) en un point quelconque lorsque le programme de définition de sortie
de volume 80e est en cours d'exécution.
Sur la Figure 30, la sortie de volume 80el, l'écran de fenêtre 80e3 92 (intitulée "définition de volume de sortie") est illustré. Sur la Figure 30, l'écran de fenêtre 92 comprend une section d'informations 92a (études disponibles, catégories disponibles), un nom de catégorie demandée 92b (nom de
la catégorie sismique de sortie) et une description
facultative 92c (description de la catégorie sismique
de sortie).
Sur la Figure 22, lorsque le type de données de surface 80d a été sélectionné sur la Figure 15, le programme de définition de sortie de surface 80f, représenté sur la Figure 22, génère une sortie de surface 80fl sur l'écran de fenêtre 94 de l'écran 32a du poste de travail 32 de la Figure 6. L'écran de
fenêtre 94 est présenté sur la Figure 31 des dessins.
L'opérateur sélectionne une nouvelle surface sismique 80f2, et, quand la nouvelle surface sismique 80f2 est sélectionnée, le bloc de surface sélectionnée 80f3 est exécuté. Lorsque le bloc de surface sélectionnée 80f3 est exécuté, la définition de sortie de surface spécifiée est présentée sur l'écran de fenêtre de sortie de surface 94 apparaissant sur l'écran 32a du poste de travail 32. L'écran de fenêtre est indiqué sur la Figure 31. Lorsque le bloc de surface sélectionnée 80f3 est exécuté, l'opérateur du poste de travail 32 peut à titre d'alternative définir une autre sortie en option 80f4, auquel cas la sortie en option est définie en 80f5. Dans l'un ou l'autre cas, ou bien quand le bloc 80f3 de surface sélectionnée est exécuté, ou bien quand le bloc de sortie en option défini 80f5 est exécuté, l'opérateur du poste de travail 32 peut sortir
vers "l'état correct de surface de sortie" 80f6.
Naturellement, l'opérateur du poste de travail 32 peut exécuter une "prise de congé" 80f7 à un instant quelconque pendant l'exécution du programme de
définition de sortie de surface 80f.
Sur la Figure 31, la sortie de surface 80fl, l'écran de fenêtre 80f3 94 (intitulée "définition de surface de sortie") est illustré. Sur la Figure 31, l'écran de fenêtre 94 comprend un bloc de sélection de surface 94a (sélection d'horizon) et des options de
sortie de surface 94b (attributs d'horizon à générer -
temps d'horizon minimal, amplitude sismique au temps minimal). Sur les Figures 16 et 28, une construction détaillée du programme "exécution disponible" 90,
référencé sur la Figure 15, est représentée.
Sur les Figures 16 et 28, et en commençant par la Figure 16, quand le programme "d'exécution disponible" 90 de la Figure 15 est exécuté, l'opérateur voit l'écran de fenêtre 98 (représenté sur la Figure 28) sur l'écran 32a du poste de travail 32 de la Figure 6, et, en appuyant sur la touche "exécution" 42 de la Figure 28, le programme "exécution" 90a de la Figure 16 commence l'exécution du programme. Le programme exécute b après quoi, sur la Figure 28, le premier cube 34 de données (représentant les données sismiques 3D obtenues à l'instant "tl") est mathématiquement traité sur le second cube 36 de données (représentant les données sismiques 3D obtenues à l'instant "t2") en utilisant l'opération mathématique 38 pour obtenir une sortie particulière 40. Naturellement, sur la Figure 16, l'opérateur du poste de travail 32 peut répéter les procédures ci-dessus décrites dans les paragraphes précédents, en sélectionnant des données primaires et un type de données 90b (le programme revient au programme de la Figure 8), en sélectionnant des données secondaires et un type de données 90c (le programme revient au programme des Figures 9-10), en sélectionnant une opération et desopérations disponibles 90d (le programme revient au programme des Figures 11-14) et en sélectionnant une sortie et un type de données 90e (le programme revient au programme de la Figure 15), auquel cas l'opérateur devrait effectuer l'opération d'exécution 90a de la Figure 16
et le programme devrait exécuter 90b.
Sur les Figures 33 à 35, deux exemples de la sortie particulière 40 de la Figure 23 sont représentés. La Figure 33 illustre une entrée typique utilisée dans l'opération mathématique suivante: Vo = Vl * V1,
la Figure 33 représentant le paramètre Vl.
La Figure 34 illustre une "sortie particulière" typique 40 de la Figure 23 et représente le paramètre Vo de l'opération mathématique Vo = V! * V1 référencée
à la Figure 33.
La Figure 35 illustre une autre "sortie particulière" typique 40 de la Figure 23 et représente la sortie "So" (une sortie de surface échantillon) dans
l'opération mathématique suivante: So = Vl n Sl.
Le logiciel Cubemath 32c2 de la Figure 6, lorsqu'il est exécuté par le processeur 32b du poste de travail en liaison avec le logiciel d'interprétation IESX 32cl, effectue des opérations mathématiques sur des volumes sismiques (c'est-à-dire une paire de volumes de données sismiques 3D) pour combiner, contraster et transformer les entrées originales en une sortie particulière, laquelle sortie particulière est utilisée par le géologue pour déterminer différentes caractéristiques d'une formation terrestre dans le but final de déterminer la teneur des dépôts d'hydrocarbures souterrains (par exemple la présence de
pétrole). Voici quelques exemples.
1. Interprétation lithologique - C'est un exemple de l'emploi de "l'opération sismique unique avec double procédé de traitement" mentionnée ci-dessus - Il existe des techniques pour analyser des données sismiques avant traitement par empilement pour les
effets du décalage en fonction de l'amplitude (AVO).
Deux volumes de sortie créés pour ce procédé sont l'empilement d'ondes de pression (P) et le volume de gradient (G). La combinaison sous forme d'un produit (P x G) des deux grandeurs ci-dessus peut être utilisée par un géologue comme indication lithologique. Leur somme (P + G) est liée au rapport de Poisson, une
propriété élastique de la roche. Leur différence (P -
G) est une estimation de la réponse de l'onde de cisaillement.
(2) Transformation temps en profondeur -
L'analyse des données sismiques de pré-empilage et des informations de puits peut conduire à un volume de vitesse moyenne qui peut être combiné avec des informations sismiques pour une transformation en profondeur. Le procédé implique plusieurs opérations, à savoir: multiplier la vitesse moyenne de chaque échelon de données sismiques de volume et de temps, et l'indice de l'échelon de temps pour obtenir un indice de profondeur (V x T), puis interpoler l'amplitude sismique sur l'indice de profondeur désiré
(interpolation sin x/x).
(3) Analyse des données sismiques 4D (relevés sismiques répétés pour la surveillance des réservoirs) - Une analyse du contraste entre deux relevés sismiques acquis au même emplacement peut indiquer la migration des fluides des réservoirs. Une simple soustraction entre deux tracés pourrait servir à révéler ces différences si l'acquisition des deux ensembles de données a été enregistrée de manière identique. D'une façon générale, il existe de légères différences dont on peut tenir compte par des opérations d'addition pour ajuster les différences de temps et/ou d'amplitude entre les ensembles de données, c'est-à-dire par
corrélation croisée.
(4) Variations des échantillons - Beaucoup d'opérations à simple volume, à volume/surface, à volume/tracé, peuvent être effectuées pour des applications telles la détection des bords, le tranchage des volumes, les tests des hypothèses et le traitement de post-empilage. Ces applications exigent une combinaison de fonctions mathématiques de base
prévues selon la présente invention dans ce descriptif.
Une description fonctionnelle de l'utilisation de
la présente invention est présentée dans les paragraphes suivants en référence aux Figures 1-35 des dessins. Le camion d'enregistrement 20 de la Figure 1 est placé à la surface de la terre, comme l'indique la Figure 1. La source d'énergie explosive 30 transmet des vibrations sonores 14 dans la terre en réponse à l'explosion de l'énergie sonore 10, les vibrations sonores 14 à l'origine se propagent vers le bas comme indiqué sur la Figure 1, elles se réfléchissent sur une couche 16 dans la terre, et se propagent vers le haut jusqu'à la surface. Les vibrations sonores 14 sont reçues et mesurées par une pluralité de géophones 18 disposés à la surface de la terre. Les géophones 18 transmettent un ou plusieurs signaux électriques, en réponse à la réception des vibrations sonores 14, les signaux électriques étant reçus dans l'ordinateur du camion d'enregistrement 20a. Sur la Figure 2, le camion d'enregistrement 20 s'est déplacé jusqu'à une distance d'environ 15 m (50 pieds) en une seconde position (position 1) et la procédure ci-dessus a été répétée, après quoi un autre ensemble d'un ou plusieurs signaux électriques est reçu dans l'ordinateur du camion d'enregistrement 20a. Le camion d'enregistrement 20 se déplace encore huit fois jusqu'à la position 9 de la Figure 2 dans laquelle un autre ensemble d'un ou plusieurs signaux électriques est reçu dans l'ordinateur du camion d'enregistrement 20a. Dix ensembles de signaux électriques sont reçus, au temps "tl" dans l'ordinateur du camion d'enregistrement 20a de la Figure 3 et mémorisés dans la mémoire 20a2. Un premier ensemble de données sismiques 3D est maintenant mémorisé dans la mémoire 20a2 de l'ordinateur du camion
d'enregistrement 20a.
En supposant qu'une "opération sismique répétée" soit effectuée à l'instant "t2" (typiquement environ 6 mois plus tard), la procédure décrite ci-dessus dans le paragraphe précédent est à nouveau répétée, et, comme résultat, un second ensemble de données sismiques 3D est mémorisé dans la mémoire 20a2 de l'ordinateur du
camion d'enregistrement 20a.
Deux (2) supports d'enregistrement de sortie 20a4 de la Figure 3 sont générés, un premier qui mémorise le premier ensemble de données sismiques 3D et un second qui mémorise le second ensemble de données sismiques 3D. Le premier et le second ensemble de données sismiques 3D mémorisés dans le premier et le second support d'enregistrement de sortie 20a4 sont tous deux
chargés dans l'ordinateur principal 30 de la Figure 4.
Chaque ensemble de données sismiques 3D subit un traitement de "réduction des données" dans l'ordinateur principal 30 de la Figure 4 en utilisant le logiciel de réduction des données 30b représenté sur les Figures 4, 7a et 7b des dessins. Il en résulte que deux "supports d'enregistrement de sortie de données réduites" 30d de la Figure 3 sont générés, un premier support d'enregistrement de sortie de données réduites 30d mémorisant un premier ensemble de données sismiques 3D à données réduites obtenu à l'instant "tl", et (en supposant que "l'opération sismique répétée" soit effectuée) un second support d'enregistrement de sortie de données réduites 30d mémorisant un second ensemble de données sismiques 3D à données réduites obtenu à
l'instant "t2", environ 6 mois après l'instant "tl".
Le premier et le second support d'enregistrement de sortie de données réduites 30d sont chacun chargés dans le poste de travail 32, comme l'indiquent les
Figures 5 et 6.
A ce moment, comme indiqué sur la Figure 8, le processeur du poste de travail 32b de la Figure 6 se trouve à l'état initial 50 de la Figure 8. Sur la Figure 23, l'écran de fenêtre 33 est présenté à un opérateur sur l'écran 32a du poste de travail 32. Le premier cube 34 de la Figure 23 représente le premier ensemble de données réduites des données sismiques 3D qui a été mesuré par les géophones 18 et obtenu à l'instant "tl" puis mémorisé sur le premier support d'enregistrement de sortie de données réduites 30d, le second cube 36 de la Figure 23 représentant le second ensemble à données réduites des données sismiques 3D qui a été mesuré par les géophones 18 et obtenu à l'instant "t2" (en supposant qu'une "opération sismique répétée" est effectuée) et mémorisé sur le second support d'enregistrement de sortie de données réduites 30d. Sur la Figure 8, en cliquant sur le premier cube 36 de la Figure 23, l'opérateur du poste de travail 32 sélectionne des données primaires 56. Il doit alors sélectionner le type de données 58. Sur la Figure 24, l'opérateur sélectionne soit le volume-1 60a soit la
surface-i 60b. Sur la Figure 8, la sélection du volume-
1 60a de la Figure 24 correspondrait à la sélection du volume 62 de la Figure 8, et la sélection de la surface-1 60b de la Figure 24 correspondrait à la sélection de la surface 66 de la Figure 8. On suppose, dans le présent exposé, que l'opérateur a sélectionné le volume-1 60a sur la Figure 24, c'est-à-dire que, sur la Figure 8, l'opérateur a sélectionné le volume 62. A ce moment, sur la Figure 8, le bloc de sélection de volume 64 de la Figure 8 est exécuté. Un examen rapide de la Figure 17 montre que, au moment o le bloc de sélection de volume 64 de la Figure 8 est exécuté, l'écran de fenêtre 65 de la Figure 29 est présenté à l'opérateur sur l'écran 32a du poste de travail 32, et, qu'en utilisant l'écran de fenêtre 65, l'opérateur sélectionne le trajet 64b, sélectionne la plage de temps 64f, sélectionne la plage de tracé 64e et sélectionne la plage de ligne 64d, puis à ce moment, le programme revient à l'état de volume correct. Lorsque la sélection de volume 64 de la Figure 8 est exécutée, le programme 70 de volume primaire de la Figure 8 est prêt à l'exécution. Sur la Figure 9, le programme de volume primaire disponible 70 présente à l'opérateur du poste de travail 32 l'écran de fenêtre 82 représenté sur la Figure 25. Sur la Figure 25, l'opérateur clique sur le second cube 36 pour sélectionner le volume secondaire. En revenant à la Figure 9, l'opérateur sélectionne des données secondaires 70a et puis sélectionne le type de données 70b. En revenant à la Figure 24, lorsque l'opérateur a cliqué sur le second cube 36, les options suivantes sont disponibles pour la sélection: volume-2 60c, surface-2 60d, tracé 60e et rien 60f. On suppose pour l'objet de la discussion que l'opérateur clique sur "volume-2" 60c. En revenant à la Figure 9, cliquer sur volume- 2 60c de la Figure 24 signifie que l'opérateur du poste de travail 32 a sélectionné le volume 70c de la Figure 9. Quand la sélection de volume 70c, associée à la sélection des données secondaires 70a de la Figure 9, est terminée, le programme de sélection de volume 70g de la Figure 9, associé à la sélection des données secondaires, est exécuté. On rappelle que le programme de sélection 64 de la Figure 8, associé à la sélection des données primaires, a déjà été exécuté. L'exécution du programme de sélection de volume 70g nécessite un autre examen rapide de la Figure 17. Sur la Figure 17, quand le programme de sélection de volume 70g de la Figure 9 est exécuté, l'écran de fenêtre 65 de la Figure 29 est présenté à l'opérateur sur l'écran 32a du poste de travail 32, et, en utilisant l'écran de fenêtre 65, l'opérateur sélectionne un trajet 64b, sélectionne une plage de temps 64f, sélectionne une plage de tracé 64e et sélectionne une plage de ligne 64d, puis le programme revient à l'état de volume correct. Quand le programme de sélection de volume 70g de la Figure 9 est exécuté, le programme de volume/volume disponible 70j de la Figure 9 est prêt à l'exécution. Se reporter à la Figure 11 au sujet du programme de volume/volume
disponible 70j.
Le programme de volume/volume disponible 70j représenté sur la Figure 11 du logiciel Cubemath 32c2 de la Figure 6 est prêt à l'exécution. Sur la Figure 11, il faut maintenant sélectionner une opération mathématique 70jl, et un certain nombre d'opérations disponibles 70j2 sont présentés à l'opérateur. Sur la Figure 26, un écran de fenêtre 82 est présenté à
l'opérateur sur l'écran 32a du poste de travail 32.
Plusieurs opérations mathématiques disponibles 86 sont
présentées à l'opérateur au poste de travail 32.
L'opérateur utilise une souris 84 pour sélectionner une des opérations mathématiques 86. Sur la Figure 11, les opérations mathématiques 70j3, 70j4, 70j5 et 70j6 sont disponibles pour la sélection. On suppose que l'opérateur sélectionne l'opération mathématique 70j4 sur la Figure 11. Le programme d'ensemble opérateur 80 est alors prêt à l'exécution. On se réfère à la Figure
pour le programme d'ensemble opérateur 80.
Sur la Figure 15, l'exécution du programme d'ensemble opérateur 80 nécessite la sélection d'une sortie particulière, c'est-à-dire la sélection d'une sortie 80a et la sélection d'un type de données particulier 80b associé à la sortie 80a. A ce moment, l'écran de fenêtre 84 de la Figure 27 est présenté à l'opérateur sur l'écran 32a du poste de travail 32. Sur la Figure 27, l'opérateur utilise une souris 84a pour sélectionner l'une des deux sorties 84b: soit un volume-1 84bl soit une surface-1 84b2. On suppose que l'opérateur sélectionne le volume-1 84bl. Sur la Figure , la sélection du volume-1 84bl signifie que l'opérateur a sélectionné le volume 80c sur la Figure 15. La sélection du volume 80c signifie que le programme de définition de sortie de volume 80e de la Figure 15 est exécuté. On se réfère à la Figure 21 pour le programme de définition de sortie de volume 80e de
la Figure 15.
Sur la Figure 21, le volume 80c ayant été sélectionné, le programme de sortie de volume 80el présente un écran de fenêtre 92 de la Figure 30 à l'opérateur du poste de travail 32. Sur la Figure 21, l'opérateur définit un nouveau volume sismique 80e2 en indiquant, sur la Figure 30, les informations de relevé disponible et de catégorie disponible 92a, le nom de la catégorie sismique de sortie voulue 92b et la
description de la catégorie sismique de sortie 92c. Sur
la Figure 21, maintenant que le volume est défini en e3, le programme ou bien passe à l'état de volume de sortie correct 80e6, ou bien l'opérateur du poste de travail 32 peut définir une sortie en option 80e4, auquel cas, une sortie en option est définie en 80e5 et le programme peut alors passer à l'état de volume de sortie correct 80e6. Naturellement, l'opérateur peut prendre congé en 80e7 pour revenir à l'état appelant, ce qui ramène le programme à l'état qui a effectué la
demande de sélection de volume.
Sur la Figure 15, quand l'exécution du programme de définition de sortie de volume 80e est terminée, le programme d'exécution disponible 90 est exécuté. On se réfère à la Figure 16 au sujet du programme d'exécution
disponible 90 de la Figure 15.
Sur la Figure 16, lorsque le programme d'exécution disponible 90 commence à être exécuté, l'écran de fenêtre 98 de la Figure 28 est présenté à l'opérateur
sur l'écran 32a du poste de travail de la Figure 6.
L'opérateur clique sur l'icône "exécution" 42 de la Figure 28, c'est-àdire que "l'exécution" 90a de la Figure 16 est exécutée. Le programme continue à
exécuter 90b.
En revenant sur les paragraphes précédents, sur les Figures 8, 9, 11, 15, 16, 21 et 28, l'opérateur effectue les sélections suivantes: (1) l'opérateur sélectionne des données primaires 56 en cliquant sur le premier tube 34, (2) l'opérateur sélectionne le type de données pour les données primaires sélectionnées, par exemple "le volume 62", (3) l'opérateur sélectionne des données secondaires 70a sur la Figure 9, (4) l'opérateur sélectionne le type de données des données secondaires telles que le volume 70c, (5) l'opérateur sélectionne une opération mathématique 70jl, et, dans l'exemple du mode de réalisation préféré de la Figure 11, l'opération mathématique sélectionnée 70j4 est V3 = V1/V2, (6) l'opérateur sélectionne une sortie 80a sur la Figure 15, (7) l'opérateur sélectionne un type de données b pour la sortie sélectionnée, par exemple un volume c, (8) l'opérateur définit un nouveau volume sismique 80e2 sur la Figure 21, et (9) l'opérateur sélectionne l'icône "exécution" 42 sur la Figure 28, ce qui exécute le programme "exécution" 90a de la Figure 16 et le programme exécute b. Lorsque l'opérateur sélectionne l'icône "exécution" 42, le programme exécute 90b. Il en résulte que les données primaires 56 agissent mathématiquement sur les données secondaires 70a, l'opération mathématique citée ci-dessus étant constituée par l'opération mathématique sélectionnée 70j4, laquelle a été sélectionnée par l'opérateur. Une sortie 80a est générée (telle que les sorties représentées sur les Figures 34 et 35), laquelle reflète les résultats de l'opération mathématique, la sortie 80a étant présentée à l'opérateur sous la forme du type de données sélectionné (c'est-à- dire présentée à l'opérateur sous la forme du volume sélectionné 80c). En plus, un nouveau volume sismique de données est créé, lequel
contient les données de la sortie 80a.
Un géologue analyse la sortie 80a dans le but de déterminer les caractéristiques de la formation terrestre et l'existence éventuelle de dépôts souterrains d'hydrocarbures. Par exemple, comme indiqué précédemment, le géologue peut utiliser la sortie 80a pour effectuer: (1) Une interprétation lithologique - C'est un exemple de l'utilisation de "l'opération sismique
unique avec deux procédés de traitement" mentionnée ci-
dessus - Il existe des techniques pour analyser des données sismiques de pré-empilage en recherchant les
effets du décalage en fonction de l'amplitude (AVO).
Deux volumes de sortie créés pour ce procédé sont la pile des ondes de pression (P) et le volume de gradient (G). La combinaison sous forme d'un produit (P x G) de ces deux grandeurs est utilisable par des géologues comme indication lithologique. Leur somme (P + G) est liée au rapport de Poisson, une propriété élastique de la roche. Leur différence (P - G) est une estimation de
la réponse de l'onde de cisaillement.
(2) Transformation temps en profondeur -
L'analyse des données sismiques de pré-empilage et des informations de puits peut conduire à un volume de vitesse moyenne qui peut être combiné avec les informations sismiques pour une transformation en profondeur. Le procédé implique plusieurs opérations, à savoir: multiplier la vitesse moyenne de chaque échelon de données sismiques de volume et de temps et l'indice d'échelon de temps pour obtenir un indice de profondeur (V x T), puis interpoler l'amplitude sismique sur l'indice de profondeur désiré
(interpolation sin x/x).
(3) Analyse des données sismiques 4D (relevés sismiques répétés pour la surveillance des réservoirs) - Une analyse du contraste entre deux relevés sismiques acquis au même emplacement peut indiquer la migration des fluides des réservoirs. Une simple soustraction entre deux tracés pourrait être utilisée pour mettre ces différences en évidence si l'acquisition des deux ensembles de données a été enregistrée de manière identique. En général, il existe de légères différences dont on peut tenir compte par des opérations d'addition pour ajuster les différences de temps et/ou d'amplitude entre les ensembles de données,
c'est-à-dire par corrélation croisée.
(4) Variations des échantillons - Plusieurs opérations à un seul volume, à volume/surface, à volume/tracé, peuvent être effectuées pour des applications telles la détection des bords, le tranchage des volumes, les tests d'hypothèses et le traitement de post-empilage. Ces opérations nécessitent une combinaison de fonctions mathématiques de base
prévues selon la présente invention dans ce descriptif.
Après la description ci-dessus, il est évident que
l'invention peut être modifiée de différentes manières.
De telles modifications entrent dans le cadre et l'esprit de l'invention et toutes les modifications qui pourraient apparaître à un spécialiste de la technique
sont incluses dans le cadre des revendications
annexées.
Claims (10)
1. Appareil adapté à la génération d'un support d'enregistrement de sortie contenant des informations indicatives d'un ensemble résultant de caractéristiques d'une formation terrestre, comportant: des premiers moyens pour recevoir et mémoriser un premier ensemble d'informations lié à un premier ensemble de caractéristiques de ladite formation terrestre, lesdits premiers moyens comportant des moyens pour sélectionner ledit premier ensemble d'informations; des seconds moyens pour recevoir et mémoriser un second ensemble d'informations concernant un second ensemble de caractéristiques de ladite formation terrestre, lesdits seconds moyens comportant des moyens pour sélectionner ledit second ensemble d'informations; des troisièmes moyens pour sélectionner une opération particulière, ledit premier ensemble d'informations étant apte à agir sur ledit second ensemble d'informations pendant ladite opération particulière; des quatrièmes moyens pour sélectionner un type particulier de sortie; et des moyens d'initialisation répondant à la sélection dudit premier ensemble d'informations par lesdits premiers moyens, à la sélection dudit second ensemble d'informations par lesdits seconds moyens, à la sélection de ladite opération particulière par lesdits troisièmes moyens, et à la sélection dudit type particulier de sortie par lesdits quatrièmes moyens, pour initialiser une opération d'exécution, ladite opération particulière commençant en réponse à l'initialisation par lesdits moyens d'initialisation de ladite opération d'exécution, ledit premier ensemble d'informations et ledit second ensemble d'informations agissant l'un sur l'autre pendant ladite opération particulière, ledit support d'enregistrement de sortie étant généré sous la forme dudit type particulier de sortie sélectionné par lesdits quatrièmes moyens quand ledit premier ensemble d'informations est opéré sur ledit second ensemble d'informations, le type particulier de sortie dudit support d'enregistrement de sortie contenant lesdites informations liées audit ensemble résultant de
caractéristiques de ladite formation terrestre.
2. Procédé pour déterminer un ensemble de caractéristiques d'une formation terrestre, l'existence d'un dépôt souterrain d'hydrocarbures dans ladite formation terrestre étant déterminée à partir dudit ensemble de caractéristiques, comportant: (a) l'exécution d'au moins une opération sismique au-dessus de ladite formation terrestre pour obtenir un premier ensemble d'informations et un second ensemble d'informations; et (b) le traitement dudit premier ensemble d'informations et dudit second ensemble d'informations pour produire une sortie qui comporte un troisième ensemble d'informations, l'ensemble de caractéristiques de ladite formation terrestre étant déterminé à partir de ladite sortie, ladite existence d'un dépôt souterrain d'hydrocarbures étant déterminée à partir dudit ensemble de
caractéristiques de ladite formation terrestre.
3. Procédé selon la revendication 2, l'étape (a) d'exécution comportant les étapes suivantes: (al) l'exécution d'une première opération sismique au-dessus de ladite formation terrestre à un instant tl pour produire un support d'enregistrement de sortie contenant des premières informations; et (a2) la réduction des données desdites premières informations sur ledit support d'enregistrement de sortie pour obtenir un premier support d'enregistrement de sortie à données réduites contenant ledit premier
ensemble d'informations.
4. Procédé selon la revendication 3, l'étape (a)d'exécution comportant en outre les étapes suivantes (a3) l'exécution d'une seconde opération sismique au-dessus de ladite formation terrestre à un instant t2 pour produire un second support d'enregistrement de sortie contenant des secondes informations, ledit instant t2 étant différent dudit instant tl; et (a4) la réduction desdites secondes informations dudit second support d'enregistrement de sortie pour obtenir un second support d'enregistrement de sortie à données réduites contenant ledit second ensemble d'informations.
5. Procédé selon la revendication 2, 3 ou 4, l'étape d'exploitation (b) comportant les étapes suivantes: (bl) la sélection dudit premier ensemble d'informations; (b2) la sélection dudit second ensemble d'informations; (b3) la sélection d'une opération mathématique particulière; (b4) la sélection d'un type particulier de sortie; et (b5) pendant ladite opération mathématique particulière sélectionnée à l'étape de sélection (b3), l'exploitation dudit premier ensemble d'informations sélectionné pendant l'étape de sélection (bl) sur ledit second ensemble d'informations sélectionné pendant l'étape de sélection (b2) pour produire ledit type particulier de sortie sélectionnée pendant l'étape de sélection (b4), ledit type particulier de sortie
comportant ledit troisième ensemble d'informations.
6. Appareil d'interprétation de données sismiques apte à générer un enregistrement de sortie contenant des informations concernant un ensemble de caractéristiques d'une formation terrestre, lesdites caractéristiques de ladite formation terrestre étant déterminées à partir dudit enregistrement de sortie, comportant: des premiers moyens pour obtenir un premier ensemble de données représentatif desdites caractéristiques de ladite formation; des seconds moyens pour obtenir un second ensemble de données représentatif desdites caractéristiques de ladite formation; des moyens d'exploitation répondant audit premier ensemble de données obtenu desdits premiers moyens et audit second ensemble de données obtenu desdits seconds moyens pour permettre audit premier ensemble de données d'opérer sur ledit second ensemble de données et générer un ensemble de résultats quand le premier ensemble de données opère sur le second ensemble de données; et des moyens de génération d'un enregistrement de sortie répondant audit ensemble de résultats pour générer un enregistrement de sortie qui contient ledit ensemble de résultats, ledit ensemble de résultats représentant lesdites informations concernant ledit ensemble de
caractéristiques de ladite formation terrestre.
7. Appareil d'interprétation de données sismiques selon la revendication 6, lesdits premiers moyens comportant: des premiers moyens d'opération sismique pour effectuer une opération sismique 3D sur ladite formation terrestre à un instant tl, grâce à quoi un premier ensemble de données sismiques 3D est généré; des moyens de réduction de données répondant audit premier ensemble de données sismiques 3D pour réduire les données dudit premier ensemble de données sismiques 3D en générant ainsi un premier ensemble de données sismiques 3D à données réduites; et des premiers moyens de mémoire et d'affichage répondant audit premier ensemble de données réduites des données sismiques 3D pour recevoir et mémoriser ledit premier ensemble de données sismiques 3D à données réduites sous la forme d'un premier ensemble de données, lesdits premiers moyens de mémoire et d'affichage affichant un premier ensemble de données qui est représentatif desdites caractéristiques de ladite formation existant dans ladite formation à
l'instant tl.
8. Appareil d'interprétation de données sismiques selon la revendication 7, lesdits seconds moyens comportant: des seconds moyens d'opération sismique pour effectuer une opération sismique 3D sur ladite formation terrestre à un instant t2, en générant ainsi un second ensemble de données sismiques 3D, ledit instant t2 étant différent dudit instant tl, lesdits moyens de réduction de données réduisant ledit second ensemble de données sismiques 3D en réponse audit second ensemble de données sismiques 3D et générant, en réponse, un second ensemble de données sismiques 3D à données réduites; et des seconds moyens de mémoire et d'affichage répondant audit second ensemble de données sismiques 3D à données réduites pour recevoir et mémoriser ledit second ensemble de données sismiques 3D à données réduites sous la forme dudit second ensemble de données, lesdits seconds moyens de mémoire et d'affichage affichant ledit second ensemble de données, lequel est représentatif desdites caractéristiques de ladite formation existant dans ladite formation audit
instant t2.
9. Appareil d'interprétation de données sismiques selon la revendication 6, 7 ou 8, lesdits moyens d'exploitation comportant: des premiers moyens de classification pour classer ledit premier ensemble de données soit comme un volume soit comme une surface; des seconds moyens de classification pour classer ledit second ensemble de données soit comme un volume soit comme une surface soit comme un tracé; des moyens de sélection d'opération pour sélectionner une opération mathématique particulière; et des moyens d'exécution pour effectuer ladite opération mathématique particulière sur le volume ou la surface dudit premier ensemble de données et sur le volume ou la surface ou le tracé dudit second ensemble de données et répondant à ces données pour générer ledit enregistrement de sortie, lequel contient ledit ensemble de résultats qui sont indicatifs desdites informations liées audit ensemble de caractéristiques
de ladite formation terrestre.
10. Appareil d'interprétation de données sismiques selon la revendication 9, comportant en outre: des moyens de sélection de sortie répondant à la sélection de ladite opération mathématique particulière par lesdits moyens de sélection d'opération pour sélectionner un type particulier de sortie; et des troisièmes moyens de classification répondant à la sélection dudit type particulier de sortie par lesdits moyens de sélection de sortie pour classer ledit type particulier. de sortie soit comme un volume soit comme une surface, lesdits moyens d'exécution effectuant ladite opération mathématique particulière sur le volume ou la surface dudit premier ensemble de données et sur le volume ou la surface ou le tracé dudit second ensemble de données et répondant pour générer le volume ou la surface dudit type particulier de sortie, ledit type particulier de sortie étant ledit enregistrement de sortie qui contient ledit ensemble de résultats qui sont indicatifs desdites informations concernant ledit ensemble caractéristique de ladite
formation terrestre.
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US08/717,007 US5873051A (en) | 1996-09-20 | 1996-09-20 | Method and apparatus for processing at least two seismic data sets during a step to derive a third data set |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2753800A1 true FR2753800A1 (fr) | 1998-03-27 |
| FR2753800B1 FR2753800B1 (fr) | 2004-07-23 |
Family
ID=24880348
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR9711698A Expired - Fee Related FR2753800B1 (fr) | 1996-09-20 | 1997-09-19 | Procede et appareil pour exploiter des paires de donnees collectees au cours d'une operation sismique, lesquelles representent des caracteristiques d'une formation terrestre et pour generer un enregistrement de sortie |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5873051A (fr) |
| FR (1) | FR2753800B1 (fr) |
| GB (1) | GB2317451B (fr) |
| NO (1) | NO326451B1 (fr) |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB9813760D0 (en) * | 1998-06-25 | 1998-08-26 | Geco Prakla Uk Ltd | Seismic data signal processing method |
| US6201884B1 (en) * | 1999-02-16 | 2001-03-13 | Schlumberger Technology Corporation | Apparatus and method for trend analysis in graphical information involving spatial data |
| GB2386455B (en) * | 2000-12-18 | 2005-05-04 | Schlumberger Holdings | Seismic signal processing method and apparatus for generating correlation spectral volumes to determine geologic features |
| GB2372567B (en) * | 2001-02-22 | 2003-04-09 | Schlumberger Holdings | Estimating subsurface subsidence and compaction |
| WO2003003051A1 (fr) * | 2001-06-29 | 2003-01-09 | Exxonmobil Upstream Research Company | Procede d'analyse de courbure de reflexion dans des volumes de donnees sismiques |
| US7092824B2 (en) * | 2003-10-20 | 2006-08-15 | Ascend Geo Llp | Methods and systems for interactive investigation of geophysical data |
| US8346695B2 (en) * | 2007-03-29 | 2013-01-01 | Schlumberger Technology Corporation | System and method for multiple volume segmentation |
| AU2008251945B2 (en) * | 2007-05-09 | 2013-05-02 | Exxonmobil Upstream Research Company | Inversion of 4D seismic data |
| US8803878B2 (en) * | 2008-03-28 | 2014-08-12 | Schlumberger Technology Corporation | Visualizing region growing in three dimensional voxel volumes |
| US8456952B2 (en) * | 2008-10-03 | 2013-06-04 | Baker Hughes Incorporated | Curve-fitting technique for determining dispersion characteristics of guided elastic waves |
| WO2010077568A1 (fr) | 2008-12-17 | 2010-07-08 | Exxonmobil Upstream Research Company | Système et procédé pour effectuer des relevés répétitifs de contrôle avec des données de contrôle éparses |
| US8705317B2 (en) | 2008-12-17 | 2014-04-22 | Exxonmobil Upstream Research Company | Method for imaging of targeted reflectors |
| WO2010077569A1 (fr) | 2008-12-17 | 2010-07-08 | Exxonmobil Upstream Research Company | Système et procédé de reconstruction de données répétitives |
| US8451683B2 (en) * | 2009-04-03 | 2013-05-28 | Exxonmobil Upstream Research Company | Method for determining the fluid/pressure distribution of hydrocarbon reservoirs from 4D seismic data |
| US8332154B2 (en) | 2009-06-02 | 2012-12-11 | Exxonmobil Upstream Research Company | Estimating reservoir properties from 4D seismic data |
| CN103534614B (zh) | 2011-03-23 | 2017-08-04 | 离子地球物理学公司 | 用于分析时移地球物理调查中的数据的方法和设备 |
| CN103534613B (zh) * | 2011-03-23 | 2018-03-13 | 离子地球物理学公司 | 用于分析后续地球物理勘测中的数据的方法和设备 |
| US11372123B2 (en) | 2019-10-07 | 2022-06-28 | Exxonmobil Upstream Research Company | Method for determining convergence in full wavefield inversion of 4D seismic data |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4516206A (en) * | 1982-10-21 | 1985-05-07 | Mcevilly Thomas V | Post-processing of seismic parameter data based on valid seismic event determination |
| US5321613A (en) * | 1992-11-12 | 1994-06-14 | Coleman Research Corporation | Data fusion workstation |
| WO1995006886A1 (fr) * | 1993-09-03 | 1995-03-09 | Armitage Kenneth Rayvenor Lust | Evaluation de lithologie de sequences seismique |
| WO1996027141A1 (fr) * | 1995-03-02 | 1996-09-06 | The Trustees Of Columbia University In The City Of New York | Procede pour identifier des circuits de migration de fluides souterrains a l'aide d'imagerie sismique en 3d et 4d |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3921126A (en) * | 1974-08-21 | 1975-11-18 | Continental Oil Co | Method of seismically delineating the permeable extent of oil-containing formations |
| NO167423C (no) * | 1989-05-31 | 1991-10-30 | Geco As | Fremgangsmaate ved samtidig innsamling av seismiske data for grunne og dype maal. |
| US4969130A (en) * | 1989-09-29 | 1990-11-06 | Scientific Software Intercomp, Inc. | System for monitoring the changes in fluid content of a petroleum reservoir |
-
1996
- 1996-09-20 US US08/717,007 patent/US5873051A/en not_active Expired - Lifetime
-
1997
- 1997-09-10 GB GB9719086A patent/GB2317451B/en not_active Expired - Fee Related
- 1997-09-19 NO NO19974339A patent/NO326451B1/no not_active IP Right Cessation
- 1997-09-19 FR FR9711698A patent/FR2753800B1/fr not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4516206A (en) * | 1982-10-21 | 1985-05-07 | Mcevilly Thomas V | Post-processing of seismic parameter data based on valid seismic event determination |
| US5321613A (en) * | 1992-11-12 | 1994-06-14 | Coleman Research Corporation | Data fusion workstation |
| WO1995006886A1 (fr) * | 1993-09-03 | 1995-03-09 | Armitage Kenneth Rayvenor Lust | Evaluation de lithologie de sequences seismique |
| WO1996027141A1 (fr) * | 1995-03-02 | 1996-09-06 | The Trustees Of Columbia University In The City Of New York | Procede pour identifier des circuits de migration de fluides souterrains a l'aide d'imagerie sismique en 3d et 4d |
Also Published As
| Publication number | Publication date |
|---|---|
| FR2753800B1 (fr) | 2004-07-23 |
| NO974339L (no) | 1998-03-23 |
| NO974339D0 (no) | 1997-09-19 |
| GB9719086D0 (en) | 1997-11-12 |
| GB2317451A (en) | 1998-03-25 |
| NO326451B1 (no) | 2008-12-08 |
| US5873051A (en) | 1999-02-16 |
| GB2317451B (en) | 1999-02-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| FR2753800A1 (fr) | Procede et appareil pour exploiter des paires de donnees collectees au cours d'une operation sismique, lesquelles representent des caracteristiques d'une formation terrestre et pour generer un enregistrement de sortie | |
| US6757614B2 (en) | Seismic signal processing method and apparatus for generating correlation spectral volumes to determine geologic features | |
| CA2284386C (fr) | Methode pour realiser en 3d avant sommation, une migration de donnees sismiques | |
| EP2369369B1 (fr) | Procédé de surveillance d'un site de stockage géologique de gaz par inversion stratigraphique de données sismiques | |
| EP2761340B1 (fr) | Procédé d'extraction de vignette d'une image d'entraînement pour contraindre la modélisation géostatistique multipoint du sous-sol | |
| EP2253970B1 (fr) | Méthode pour imager une zone cible du sous-sol à partir de données de type walkaway | |
| FR2989788A1 (fr) | Titre non renseigne. | |
| CN102016642B (zh) | 抑制地表地震数据中的残余水底能量的方法和系统 | |
| FR3023641A1 (fr) | ||
| FR2985040A1 (fr) | Device and method for deghosting variable depth streamer data | |
| FR3027688A1 (fr) | ||
| 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 | |
| EP1452889B1 (fr) | Méthode pour mesurer les similarités locales entre plusieurs cubes de traces sismiques | |
| FR2985818A1 (fr) | Dispositif et procede d'estimation de decalages temporels | |
| Liu et al. | Wave-equation diffraction imaging using pseudo dip-angle gather | |
| FR2941055A1 (fr) | Procede d'acquisition de donnees vibrosismiques concernant une zone du sous-sol, et procede d'exploration sismique incluant un tel procede | |
| EP1226456B1 (fr) | Procede de prospection sismique mettant en oeuvre un traitement sur les ondes converties | |
| FR2984525A1 (fr) | Separation de champs d'onde pour des enregistrements sismiques repartis sur des surfaces d'enregistrement non planes | |
| FR2666905A1 (fr) | Procede perfectionne de prospection sismique du sous-sol. | |
| EP1540377B1 (fr) | Methode pour compacter et filtrer des evenements geophysiques dans le but d en extraire des informations sur la nature du sou s-sol | |
| WO2016012191A1 (fr) | Procédé de filtrage adaptatif de reflexions sismiques multiples | |
| US7450468B2 (en) | Methods, systems, and computer readable media for interactively determining parameters for ocean-bottom hydrophone/geophone receiver combination deghosting | |
| FR2966253A1 (fr) | Proceder et dispositif pour retirer des fantomes de donnees sismiques | |
| CA2059507C (fr) | Methode et dispositif d'exploration du sous-sol comportant l'emission d'une suite d'impulsions sismiques | |
| EP0923051A1 (fr) | Méthode statistique multivariable pour caractériser des images que l'on à formées d'un milieu complexe tel que le sous-sol |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RN | Application for restoration | ||
| FC | Decision of inpi director general to approve request for restoration | ||
| PLFP | Fee payment |
Year of fee payment: 19 |
|
| ST | Notification of lapse |
Effective date: 20170531 |