FR 2854247 A1 20041029 FR 0305161 A 20030428 PROCEDE DE TRAITEMENT DE CUBES SISMIQUES CORRESPONDANT POUR UNE MEME ZONE AU SOL, A DIFFERENTES VALEURS DE DEPORTS SOURCE/RECEPTEUR ET/OU D'ANGLES D'INCIDENCE DOMAINE TECHNIQUE GENERAL La présente invention est relative à l'analyse sismique. Il est classique en sismique 3D d'analyser des cubes sismiques 10 correspondant, pour une même zone au sol, à différentes valeurs de déports source/récepteur (ou "offset " selon la terminologie anglo- saxonne généralement utilisée) ou encore à différentes valeurs d'angles d'incidence. Par " cube sismique ", on entend ici et dans tout le présent texte, toute collection de traces sismiques dans laquelle chaque trace sismique 15 est associée à une position donnée sur la zone au sol sur laquelle l'acquisition a été réalisée. Il est commun, dans le cadre de traitements AVO (" Amplitude versus Offset ") ou AVA (Amplitude Variation with Angle of Incidence ");, de comparer des cubes de données correspondant à des valeurs de déport 20 différentes. Généralement, ainsi que l'illustre la figure 1, la mise en évidence des variations est faite en soustrayant deux cubes de données sismiques correspondant aux deux valeurs de déport. Toutefois, lorsque la différence en terme de signal sismique est de 25 faible amplitude, les propriétés physiques que l'on cherche à mettre en évidence peuvent être facilement masquées par des artefacts comme par exemple le bruit dû à la géométrie imparfaite de l'acquisition des données sismiques. Le caractère aléatoire et non répétable du bruit de chaque 30 acquisition fait que, comme l'illustre la figure 2, la simple soustraction des données sismiques augmente le niveau de bruit par rapport au signal recherché. Sur cette figure 2, le bruit est schématisé par des traits horizontaux (cube A) et des traits verticaux (cube B) qui se superposent sur le cube soustrait (cube A-B). De plus la superposition des bruits respectifs détruit leurs organisations spatiales statistiques (si elles existent) et donc rend plus difficile leurs filtrages. Par ailleurs, les traitements " AVO " ou " AVA " permettent de 5 déterminer en fonction du déport source/récepteur ou en fonction de l'angle d'incidence un certain nombre d'attributs tels que notamment les coefficients de réflectivité et les impédances acoustiques. Or, ces attributs AVO ou AVA sont particulièrement sensibles au bruit puisqu'ils sont obtenus par régression. 10 PRESENTA TION DE L'INVENTION Un but de l'invention est de proposer une nouvelle technique qui 15 permet de pallier les inconvénients des techniques antérieures. La technique proposée permet en particulier de déterminer, à partir d'au moins deux cubes sismiques correspondant à différentes valeurs de déport source/récepteur et/ou d'angles d'incidence, de nouveaux cubes sismiques susceptibles de permettre en évidence des informations qui 20 auraient le cas échéant été plus difficiles à mettre en évidence par les techniques classiques. Un autre but de l'invenbtion également est de proposer un traitement permettant de filtrer des cubes sismiques correspondant à différentes valeurs de déport, pour en déduire au moins un cube 25 particulièrement peu bruité et donc particulièrement adapté notamment aux traitement AVO ou AVA. Ainsi, l'invention propose un procédé de traitement d'un ensemble de cubes sismiques correspondant à une même zone d'acquisition et à différentes valeurs de déport source/récepteur et/ou d'angle d'incidence,, 30 caractérisé en ce qu'on détermine au moins un cube sismique qui est une estimation de la composante commune à au moins deux cubes sismiques correspondant à des valeurs différentes de déport et/ou d'angle d'incidence. Le cube commun ainsi déterminé porte une information qui est sensiblement l'équivalent d'une " somme " filtrée. En outre, à la différence de la simple soustraction entre cubes existants mise en oeuvre dans les techniques antérieures, ce cube commun 5 est avantageusement utilisé comme cube de référence pour mettre en évidence les variations d'amplitude sur les différents cubes qui correspondent aux différentes valeurs de déports et/ou d'angles d'incidence. Ainsi, on soustrait cette estimation de composante commune d'au 10 moins l'un des deux cubes initiaux. De cette façon, on dispose de cubes différentiels qui ont été déterminés sans que les bruits des deux cubes initaux ne s'y retrouvent superposés. Et ces cubes sismiques différentiels peuvent être facilement filtrés, 15 puisque l'organisation spatiale statistique n' y a pas été détruite par la superposition de bruits. Par ailleurs, de façon également avantageuse, pour réaliser un filtrage sur un ensemble de cubes sismiques correspondant à une même zone d'acquisition et à différentes valeurs de déport source/récepteur et/ou 20 d'angles d'incidence, on détermine au moins un cube sismique qui est une estimation de la composante commune à au moins deux cubes correspondant à des valeurs successives de déport source/récepteur et/ou d'angle d'incidence. Avantageusement notamment, ce cube commun est associé, dans 25 un nouvel ensemble de cubes sismiques, à une valeur de déport et/ou d'angle d'incidence comprise entre l'une et l'autre des deux valeurs successives auxquelles correspondent les deux cubes à partir desquels une estimation de la composante commune est déterminée et on met en oeuvre sur le nouvel ensemble de cubes sismiques ainsi constitué, un traitement 30 pour la détermination d'au moins une information relative à la géologie du soussol de la zone sur laquelle la campagne d'acquisition sismique a été mise en oeuvre Egalement, l'invention propose un procédé pour le traitement de deux ensembles de cubes sismiques issus de deux campagnes d'acquisition réalisées pour une même zone au sol à deux moments différents, caractérisé en ce que pour chaque valeur de déport et/ou d'angle 5 d'incidence, on détermine une estimation d'un cube sismique correspondant à la composante commune aux deux cubes sismiques qui correspondent à ladite valeur de déport et/ou d'angle d'incidence, ainsi qu'une estimation de composantes résiduelles et en ce qu'on met en oeuvre sur les ensembles de cubes résiduels ainsi obtenus pour chacune des deux campagnes 10 d'acquisition un procédé de filtrage du type précité. Différents cubes sismiques peuvent alors être déterminés à partir des cubes filtrés ainsi obtenus. Notamment, on peut déterminer un cube sismique qui est la somme, sur l'ensemble des valeurs de déport et/ou d'angle d'incidence, des 15 différences des cubes résiduels filtrés correspondant à l'une et à l'autre des deux campagnes d'acquisition. Egalement, pour l'une au moins des campagnes d'acquisition, on détermine un cube sismique qui est la somme, sur l'ensemble des valeurs de déport et/ou d'angle d'incidence, d'une part des cubes sismiques qui 20 correspondent aux composantes communes aux deux campagnes d'acquisition et d'autre part des cubess résiduels filtrés correspondant à ladite campagne d'acquisition. Egalement encore, pour l'une au moins des campagnes d'acquisition, on détermine pour chaque valeur de déport et/ou d'angle 25 d'incidence, un cube sismique filtré qui est la somme du cube qui correspond à la composante commune aux deux campagnes d'acquisition pour ladite valeur de déport et/ou d'angle d'incidence et d'un cube résiduel fonction d'un ou plusieurs cubes résiduels filtré(s) déterminé(s) pour la valeur de déport et/ou d'angle d'incidence et au moins une valeur de déport 30 et/ou d'angle d'incidence voisine. D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative. DESCRIPTION DES DESSINS - la figure 1 est une représentation schématique qui illustre un traitement classique par soustraction des cubes sismiques; - la figure 2 est une représentation schématique qui illustre le fait que la soustraction des cubes bruités ne fait qu'augmenter le niveau de bruit par rapport aux variations que l'on met en évidence; - la figure 3 est une représentation schématique qui illustre la détermination d'un cube commun et la détermination, à partir de ce cube commun, de cubes de différences; - la figure 4 illustre un mode de mise en oeuvre possible pour la détermination d'une estimation de la partie commune à deux cubes; 15 - la figure 5 illustre un exemple de masque appliqué aux cartographies 2D; - la figure 6 illustre un autre mode de mise en oeuvre possible avec plusieurs itérations; - la figure 7 illustre un autre mode de mise en oeuvre possible pour 20 calculer deux cubes avec une phase commune. DESCRIPTION D'UN OU PLUSIEURS EXEMPLES DE MISE EN OEUVRE EXEMPLE GENERAL On considère deux cubes sismiques (données A et données B) obternus à partir d'une ou plusierurs acquisitions réalisées pour une même zone au sol et correspondant à deux valeurs de déport et/ou d'angles d'incidence différentes. Les données qui correspondent à ces deux cubes peuvent être décomposées ainsi: données A = partie commune + bruit A données B = partie commune + bruit B + A o à représente les variations des données sismiques et o la partie commune se définit comme la partie cohérente ou invariante des données A et B. Dans une première étape, on détermine une estimation de cette partie commune. Différentes méthodes peuvent être utilisées à cet effet. Notamment, cette estimation peut être calculée par co-krigeage, ou encore par cross-corrélation ou par calcul de valeurs moyennes. Cette partie commune estimée représente l'information qui est 10 restée identique entre les deux valeurs de déports et/ou d'angles d'incidence. Elle a l'avantage de présenter un meilleur rapport signal à bruit que les données initiales. Dans un deuxième temps on détermine les deux différences 15 partielles suivantes: Diff A = données A - partie commune estimée Diff B = partie commune estimée - données B La figure 3 illustre le cube de composante commune et les cubes de différences partielles (cubes différentiels) ainsi obtenus. Les deux différences partielles permettent de mettre en évidence les évolutions pour chaque déport ou valeur d'incidence. Elles vérifient: DiffA = bruit A - a.A Diff B = - bruit B - (1-a)A o OE est le rapport de répartition du signal résiduel, lequel est souvent proche de 0.5 mais peut varier spatialement. (Il peut dependre aussi du signal sur bruit de chaque cube initiaux). Le rapport Signal/Bruit de ces deux différences partielles est alors o.A I bruit A pour Diff A et (1 -cx)A I bruit B pour Diff B. On met ensuite en oeuvre sur les différences partielles ainsi obtenues un traitement de filtrage, ce qui est possible grâce au fait que les organisations spatiales statistiques sont conservées et en particulier n'ont pas été supprimées par la superposition de bruits et que le bruit de chaque jeux de données conserve ses propriétés intrinsèques initiales, comme sa cohérence spatiale statistique. Le filtrage indepedant pour chaque différence partielle peut etre mis 5 en oeuvre par exemple par un filtre spatial déterministe dans le domaine KxKy ou un filtre statistique par krigeage factoriel ou meme une fonction de lissage suivant un axe particulier permettant d'atténuer les bruits de l'acquisition. L'orientation spatial du bruit est d'autant plus facile a déterminer par des mesures statistiques (variogrammes 2D13D) que la 10 partie commune contenant la majorite du signal sismique de la geologie a ete soustraite. Après filtrage, il est possible de reconstituer la différence totale à partir des deux différences partielles filtrées ainsi obtenues: données A - données B = Diff A + Diff B On comprend que la sommation des différences partielles après filtrage de celles-ci permet un meilleur rapport signal à bruit que ce qui est obtenu en calculant la différence totale par soustraction des données A et des données B. On notera qu'il n'y a pas de perte d'information, puisque la somme 20 des différences partielles est égale à la différence totale. Egalement, on notera que, si le niveau de bruit est différent entre les données A et B initiales, une des deux différences partielles présente un meilleur rapport signal sur bruit que la différence totale. Dans certains cas par conséquent, l'information sera plus apparente 25 sur un cube correspondant à une différence partielle que sur Ile cube correspondant à la différence totale. CALCUL DE PARTIE COMMUNE Principe La figure 4 illustre un exemple de calcul de la partie commune utilisant le Co-krigeage factoriel automatique spatial (ou analyse krigeante multi-variables automatique) dans le domaine des fréquences. Les données temporelles initiales sont transformées par FFT en cubes de parties réelles et de parties imaginaires dans le domaine fréquentiel. La partie commune de chaque plan fréquentiel réel et imaginaire est 5 calculée par la méthode dite " co-krigeage factoriel automatique " avec un opérateur 2D dans le plan (x, y). Par le processus de transformée de Fourier inverse du cube commun des parties réelles et du cube commun des parties imaginaires, on obtient alors le cube commun temporel des données sismiques. 10 Exemple détaillé Un exemple plus détaillé de traitement est donné ci-dessous: 1) Dans une première étape, on lit les données qui correspondent aux deux collections de traces (cube A et cube B) à traiter. Ces données 15 sont enregistrées dans un fichier o chaque trace est identifiée par un numéro de ligne et un numéro de colonne qui correspondent à une position au sol. 2) Dans une deuxième étape, on sélectionne une fenêtre de temps pour définir sur les traces sismiques l'intervalle de temps pris en 20 considération pour le traitement.. On complète le cas échéant les portions de traces ainsi sélectionnées par des valeurs nulles pour obtenir pour chaque trace un nombre d'échantillons correspondant à une puissance de 2, 3, 5 (selon la fonction FFT) pour permettre la FFT 3) Dans une troisième étape, on applique à la succession 25 d'échantillons correspondant à chacune des traces et à la fenêtre temporelle sélectionnée une transformée de Fourier rapide (FFT). On obtient alors au total quatre cubes en fréquence correspondant à la partie réelle et la partie imaginaire de chaque sélection des cubes A et B. 4) Pour chaque plan de fréquence, on construit alors, dans le domaine réel d'une part et dans le domaine imaginaire d'autre part, deux cartes (2D) correspondant l'une aux données du cube A et l'autre aux données du cube B. Puis on applique sur les cartes ainsi obtenues pour les cubes de données A et B un traitement de co-krigeage factoriel automatique afin d'en déduire une carte commune aux deux cartes initiales. Des rappels sur le cokrigeage factoriel sont donnés en Annexe I à la fin de la description. Le traitement de co-krigeage factoriel est réalisé de façon indépendante d'une part avec les partiesréelles et d'autre part avec les parties imaginaires. Ce traitement comporte par exemple, dans chacun de ces deux 10 domaines, les étapes suivantes: - En fonction d'un rayon de variogramme donné, calcul d'un variogramme 2D pour la carte correspondant aux données A et pour la carte correspondant aux données B; calcul d'un cross-variogramme, . - Construction des matrices de covariances et cross-covariance à 15 partir des 3 variogrammes ainsi obtenus. - Résolution du système de co-krigeage factoriel qui correspond à ces matrices, afin d'en déduire un opérateur spatial, c'est à dire fonction de la position en ligne et en colonne de l'échantillon considéré sur la cartographie. Pour stabiliser l'inversion de la matrice du système, on rajoute 20 sur la diagonale des matrices de covariance une constante qui correspond à un tres petit pourcentage de la racine carrée du produit des deux variances. - On convolue alors l'opérateur spatial ainsi obtenu avec les cartes initiales pour obtenir la partie commune des deux cartes. 5) Le traitement 25 qui vient d'être décrit est répété pour chaque plan fréquence des parties réelles et imaginaires. On obtient deux cubes, un pour les parties réelles communes et l'autre pour les parties imaginaires communes. 6) Une fois tous les plans fréquences calculés, on met en oeuvre une transformée de Fourier inverse des cubes complexes communs pour 30 obtenir un cube en temps de sismique commune. EXEMPLES DE VARIANTES DE MISE EN îUVRE D'autres variantes de réalisation sont bien entendu envisageables. Masque de traitement Ainsi que l'illustre la figure 5, le calcul de l'opérateur peut être limité à des régions précises correspondant à des données valides. Les 5 variogrammes et cross-variogramme sont ainsi calculés sur les régions définies par le masque. Le masque peut par exemple cacher une zone sans données due à une plate-forme de forage ou des données très bruitées par des effets de bords. L'utilisation d'un masque va améliorer la qualité de l'opérateur et donc du filtrage de bruit. Le masque est défini avec des 10 valeurs valides en blanc (1) et des valeurs non valides en noir (0). Double itération Egalement, lorsque le contenu fréquentiel du signal entre les cubes A et B est très différent, il est parfois avantageux de faire le calcul du cube 15 commun en deux itérations: (Figure 6) La première itération consiste à calculer un spectre commun d'amplitude. A cet effet, pour chaque plan de fréquence, on détermine à partir des cartes de parties réelles et imaginaires, une cartographie d'amplitude, 20 ainsi qu'une cartographie de phase. Puis, à partir des cartographies d'amplitude ainsi déterminées pour les données A et les données B pour chaque fréquence, on détermine une cartographie commune d'amplitude spectrale. Cette cartographie est de façon appropriée déterminée par cokrigeage factoriel automatique. La partie réelle et imaginaire de A et B est ensuite reformulée avec l'amplitude commune ainsi obtenue et la phase initiale de A et B. Après un traitement de FFT inverse des cubes, on obtient deux cubes de sismique A' et B' ayant le même contenu fréquentiel. Ces deux cubes sont ensuite traités avec une seconde itération. La seconde itération calcule la partie commune réelle et imaginaire de A' et B'. En fait, elle va fixer la phase commune des deux cubes sismiques A' et B' puisque le spectre d'amplitude est déjà commun après la première itération. On détermine ensuite la partie commune, de la même façon que ce qui a été précédemment décrit en référence à la figure 4. Une autre variante encore possible de mise en oeuvre est illustrée sur la figure 7. Dans cette variante de mise en oeuvre, on extrait des cubes d'amplitude associés respectivement au cube A et au cube B. En parallèle ou dans une deuxième étape, on sort des parties réelles et imaginaires des cubes A et B, des cubes correspondant au cosinus et au sinus de la phase des données. On détermine ensuite des cubes correspondant l'un à la partie commune des cubes de cosinus, l'autre à la partie commune des cubes de sinus. Puis on détermine un cube dont les données sont le ratio des données du cube de sinus et du cube de cosinus ainsi obtenus: on dispose 15 alors d'un cube donnant la partie commune de la phase. Ce cube de phase commune est combiné avec les deux cubes d'amplitude initialement extraits des cubes A et B pour en déduire des cubes filtrés A' et B'. Ces deux cubes de données filtrées peuvent eux-mêmes faire l'objet 20 de la seconde itération décrite en référence à la figure 5. En variante encore, il est possible de calculer le cube commun aux cubes de données A et B en combinant d'une part le cube de partie commune de phase et d'autre part le cube de partie commune d'amplitude. Autres traitements que le co-krigeage factoriel 2D D'autres traitements d'estimation de partie commune, et notamment de traitements mettant en oeuvre un co-krigeage factoriel, sont 30 envisageables. On peut notamment déterminer une partie commune directement dans le domaine temps sur une succession de cartes (x,y) d'amplitude sismique (ou d'autres données). Pour chaque instant considéré, on détermine alors la partie commune aux deux cartographies qui correspondent à cet instant en mettant en oeuvre un traitement de co-krigeage automatique factoriel. Egalement, le traitement peut utiliser des opérateurs 3D au lieu 5 d'opérateurs 2D. Le traitement de co-kriegeage automatique factoriel se généralise en effet avec une matrice de co-variance à trois dimensions. Il est à noter que l'utilisation d'opérateurs 2D (x,y) (ou 3D le cas échéant) rend la méthode très performante pour filtrer des bruits spatialement organisés. Bien entendu, également, il est possible de déterminer des parties communes à plus de deux cubes, en calculant par exemple des parties communes de parties communes. De cette façon, le procédé proposé peut être étendu à un nombre de cubes supérieurs à deux. TRAITEMENT DE CUBES CORRESPONDANT A DES VALEURS DE DEPORT ET/OU D'ANGLES D'INCIDENCE SUCCESSIVES On se place à nouveau dans un cas o l'on dispose initialement de plusieurs cubes sismiques correspondant, pour une même zone au sol, à une pluralité de valeurs de déports source/récepteur ou "offsets" différentes 20 (ou à une pluralité de valeurs d'angles d'incidence). Dans la suite du présent texte, Si désigne la collection de traces sismiques correspondant au ième déport, ou i est un indice muet. On traite deux à deux ces différentes collections de traces sismiques pour déterminer les parties communes aux collections 25 correspondant à deux valeurs de déports successives. Ainsi, pour les séries Si et Si + 1, o i C {1, ..., M-1}, M étant le nombre de collections de traces, on détermine une partie commune Cij + 1 et des parties résiduelles Ri et Ri+1qui vérifient: Si = CG,i+1 + Ri Si+, = Cii+l + Ri+, Là aussi, différentes méthodes de décomposition en parties commune et parties résiduelles sont envisageables. Une technique de décomposition qui peut être mise en oeuvre consiste à échantillonner les traces sismiques pour transformer chaque 5 collection de traces en une série de cartographies 2D. Chaque cartographie est par exemple constituée des amplitudes des traces sismiques pour un instant d'échantillonnage considéré, en fonction de la position au sol (laquelle est donnée par exemple par référence à un référentiel défini par deux axes). Les cartographies 2D obtenues pour deux collections de traces sismiques Si, Si+1 successives et pour un instant d'échantillonnage donné sont traitées par co-krigeage factoriel pour en déduire une cartographie commune et des cartographies résiduelles. L'ensemble des cartographies communes ainsi obtenues permet de 15 reconstituer la collection de traces commune. De même, l'ensemble des cartographies résiduelles constitue les collections de traces résiduelles. On notera que les collections de traces COj + 1 que l'on obtient ainsi constituent des collections débruitées Les données résiduelles Ri ou Ri + 1 qui sont supprimées correspondent en effet à la composante non répétable des données entre deux déports successifs, ce qui comprend en particulier le bruit organisé, puisque celui-ci n'est pas répétable d'un déport à un autre. Les collections de traces Cj, +, sont alors particulièrement adaptées 25 à la mise en oeuvre de traitements AVO ou AVA classiques. Notamment, une collection Cj, + 1 donnée peut être associée à une valeur de déport comprise entre les valeurs des déports i et i+1 qui correspondent aux collections successives à partir desquelles elles ont été calculées. On dispose alors d'un nouvel ensemble de collections de traces sismiques correspondant à une pluralité de valeurs de déports source/récepteur, qui sont, le cas échéant, sommées partiellement ou complètement ou à partir desquels on détermine des cartographies d'attributs AVO ou AVA. Les traitements AVO ou AVA sont en eux-mêmes classiquement connus. Pour un exemple de traitement en ce sens, on pourra avantageusement se référer à la publication suivante: On notera que les traitements AVO ou AVA prenant en compte plusieurs collections successives, l'éventuelle perte de signal lors de la détermination de la partie commune est compensée par la réduction particulièrement importante du bruit organisé que permet le fait de passer 10 par une estimation de la partie commune. Les attributs AVO ou AVA ainsi obtenus sont d'autant plus fiables. TRAITEMENT LIMITE A DES COLLECTIONS ET DES VALEURS DE DEPORT SPECIFIQUES Dans une autre variante de mise en oeuvre encore, on ne met en oeuvre le traitement de filtrage décrit ci-dessus que sur les collections de traces qui présentent un bruit particulièrement important. Parmi les différentes collections de traces dont on dispose pour une pluralité de valeurs de déport, on repère celles qui sont les plus bruitées et 20 pour chacune de ces collections on calcule au moins une collection de traces correspondant à la composante commune à la collection considérée et à la collection correspondant à la valeur de déport juste précédente et/ou juste suivante. On utilise ensuite pour mettre en oeuvre un traitement AVO ou AVA 25 un ensemble de collections de traces constitué d'une part des collections de traces initiales qui n'avaient été comme étant dans les plus bruitées et d'autre part des collections de traces communes déterminées et correspondant à de nouvelles valeurs de déport. Dans ce cas, la collection de traces sur laquelle le traitement AVO 30 ou AVA est mis en oeuvre ne correspond pas une valeur de déport régulière. APPLICATION 4D On se place ici dans le cas o l'on dispose de deux ensembles SAi, SBi de collections de traces sismiques issues de deux campagnes d'acquisition A et B réalisées pour une même zone au sol à deux périodes de temps différentes. Pour chaque déport i, on calcule une collection de composante commune Ci et des collections de composantes résiduelles RAi, RBi: SAi = Ci + RAi SBi = Ci + RBi avec i C {1, ..., M} et M nombre de collections pour chaque 10 campagne. On met ensuite en oeuvre sur les ensembles de collections de traces résiduelles RAi et RBi ainsi obtenues un traitement de filtrage du type de celui décrit ci-dessus pour l'application 3D. Notamment, on calcule les collections de traces KAi,Ai+1 (KBi,Bi+1) correspondant aux parties communes de collections RAi (Rbi) correspondant à des déports i, i+1 successifs, - ainsi que les collections de traces correspondant aux parties résiduelles non répétables MAi, MBi. 20 tels que RAi = KAi,Ai+1 + MAi RBi = KBi,Bil+1 + MBi On dispose alors d'un certain nombre de collections des traces qui sont autant de signatures de l'évolution du sous-sol entre deux acquisitions. 25 - Différence 4D sur les traces originales:: I 1M(SAi - SBi) ou '1M(RAi - RBi) - Différence 4D filtrée: -lM-(KAi,Ai+l - KBi,Bi+1) 30 - Somme filtrée: -N-1 KAi,Ai+1 + ' Nci et '1N-1KBiBi+ + iNCi - Collections de traces filtrées: C1 + (KA1,A2)/2 et C1 + (KBI,B2)/2 pour le premier déport Ci + (KAI-1,Ai + KAi,Ai+I)/2 et Ci + (KBi-.1,Bi + KBiBi+ )/2 i C {2, ..., M-1} CM + (KAM-1,AM)/2 et CM + (KBM- ,BM)/2 pour le dernier déport. Annexe I - Rappel sur le co-krigeage factoriel On décrit ci-dessous un exemple d'obtention d'une cartographie commune par co-krigeage factoriel. On désigne par Z1 et Z2 les deux fonctions correspondant à ces deux cartographies. Dans un premier temps on calcule le variogramme croisé de ces deuxfonctions, dont les valeurs sont: i 712(h)= - (Zl(x)Zl(x+h))(Z2(x)-Z2(x+h)) o x et x+h désignent les couples de points pris en considération selon la 10 direction et pour la distance h pour lesquels on détermine la valeur du variogramme, et o N est le nombre de ces couples de points pour cette direction et cette distance. Connaissant ce variogramme croisé, on détermine ensuite une 15 estimation de la fonction qui lui correspond, laquelle vérifie: Z*12 (X)= E 1 Z +I + 2 Z2 o a et P sont deux indices muets désignant les points considérés autour du point x pour lequel on cherche à déterminer une estimation de ladite fonction, Z' et Z' étant la valeur audit point x, N étant le nombre de ces 20 points et o ?' et X sont des coefficients de pondération. Ces coefficients de pondération;L et Dû se déterminent par inversion de l'équation de co-krigeage: Cull... CHN1 C1211... C,211, ' C,21X CHlN... CIINN C1211... C12NN ' N C12NX C2111... C21N1 C2211... C22N1 2l C121X C21N... C21NN C2211... C22NN... ... .. C1NX... _C2,1N... C21NN C2211.. Cn2NN 22X_ Ci2NX o les coefficients C12aI3 et C21aI sont les valeurs de variance croisée des fonctions Zl et Z2 aux points correspondants aux indices a et f3 o les coefficients Cl Iap et C22ap sont les valeurs de covariance respectivement de la fonction ZI et de la fonction Z2 auxdits points. L'indice X correspond au point désigné par x précédemment. On notera que la matrice qui apparaît dans cette équation a l'avantage d'être inversible sous certaines conditions de calcul. De cette façon, en utilisant les covariances expérimentales, on décompose de façon automatique les deux variables correspondant aux 10 deux séries de données initiales en une composante commune et deux composantes résiduelles orthogonales. La régularité des données fait que la covariance expérimentale est connue pour toutes les distances utilisées. aucune interpolation n'est nécessaire et la matrice est définie positive La fonction alors obtenue est une estimation de la composante 15 commune aux deux séries de données que constituent les deux cartographies. REVENDICATIONS 1. Procédé de traitement d'un ensemble de cubes sismiques 5 correspondant à une même zone d'acquisition et à différentes valeurs de déport source/récepteur et/ou d'angle d'incidence, caractérisé en ce qu'on détermine au moins un cube sismique qui est une estimation de la composante commune à au moins deux cubes sismiques correspondant à des valeurs différentes de déport et/ou d'angle d'incidence. 2. Procédé selon la revendication 1, caractérisé en ce qu'on soustrait cette estimation de composante commune d'au moins l'un des deux cubes initiaux. 3. Procédé selon la revendication 2, caractérisé en ce qu'on met en oeuvre un traitement de filtrage sur un cube différentiel ainsi obtenu. 4. Procédé selon les revendications 2 ou 3, caractérisé en ce qu'on soustrait à l'un et à l'autre des deux cubes initiaux le cube qui correspond à l'estimation de leur composante commune, en ce qu'on met en oeuvre un traitement de filtrage sur l'un et l'autre des deux cubes différentiels ainsi obtenus et en ce qu'on soustrait ou somme ces cubes 20 différentiels filtrés pour reconstituer une estimation de la différence entre les deux cubes initiaux. 5. Procédé selon l'une des revendications précédentes, caractérisé en ce que pour déterminer un cube qui est une estimation d'une composante commune à au moins deux cubes sismiques: - on convertit chacune de ces deux cubes sismiques dans le domaine fréquentiel pour obtenir pour chacun d'eux un cube réel et un cube imaginaire, - on calcule un cube sismique qui est une estimation de la composante commune aux cubes réels ainsi obtenus et un 30 cube qui est une estimation de la composante commune aux deux cubes imaginaires ainsi obtenus, et - on convertit dans le domaine temporel les deux cubes, l'un réel, l'autre imaginaire, ainsi obtenus. 6. Procédé selon l'une des revendications précédentes, caractérisé en ce que, pour déterminer un cube qui est une estimation d'une composante commune à au moins deux cubes: - on convertit chacun de ces cubes dans le domaine fréquentiel 5 pour en déduire pour chacun d'eux un cube d'amplitude et un cube de phase, - on calcule une estimation de la partie commune aux cubes d'amplitude ainsi obtenus. 7. Procédé selon la revendication 6, caractérisé en ce que: - on convertit dans le domaine temporel le cube d'amplitude ainsi obtenu en l'associant respectivement aux différents cubes de phase, - on convertit dans le domaine fréquentiel chacun des cubes ainsi obtenus dans le domaine temporel pour obtenir pour 15 chacun de ces cubes un cube réel et un cube imaginaire, - on calcule un cube qui est une estimation de la composante commune aux cubes réels et un cube qui est une estimation de la composante commune aux cubes imaginaires, et - on convertit dans le domaine temporel les deux cubes, l'un 20 réel, l'autre imaginaire, ainsi obtenus. 8. Procédé selon l'une des revendications précédentes, caractérisé en ce que pour déterminer un cube qui est une estimation d'une composante commune à au moins deux cubes: - on convertit chacun de ces cubes dans le domaine fréquentiel 25 pour en déduire pour chacun d'eux un cube d'amplitude, ainsi qu'un cube correspondant à sa partie réelle et un cube correspondant à sa partie imaginaire, - on extrait des cubes correspondant à la partie réelle et des cubes correspondant à la partie imaginaire des cubes cosinus 30 et sinus, - on calcule une estimation de la partie commune aux cubes de cosinus et une estimation des cubes de sinus, - on déduit des cubes de parties communes ainsi calculés un cube de données de phase communes, 9. Procédé selon la revendication 7, caractérisé en ce que: - on convertit dans le domaine fréquentiel chacun des cubes 5 ainsi obtenus dans le domaine temporel pour obtenir pour chacun de ces cubes un cube réel et un cube imaginaire, - on calcule un cube qui est une estimation de la composante commune aux cubes réels et un cube qui est une estimation de la composante commune aux cubes imaginaires, et - on convertit dans le domaine temporel les deux cubes, l'un réel, l'autre imaginaire, ainsi obtenus. 10. Procédé selon l'une des revendications 5 à 9, caractérisé en ce que pour calculer un cube qui est une estimation d'une composante commune à au moins deux cubes sismiques, on met en oeuvre un 15 traitement de cokrigeage factoriel. 11. Procédé selon la revendication 10, caractérisé en ce que pour calculer un cube qui est une estimation d'une composante commune à au moins deux cubes sismiques: - on échantillonne chacune de ces cubes sissmiques pour en 20 déduire pour chacune d'eux une succession de cartographie 2D - on détermine pour chaque couple de cartographies 2D ainsi obtenu des matrices de co-variance et de cross-variances, correspondant aux variogrammes de l'une et l'autre de ces 25 deux cartographies, ainsi qu'à leur variogramme croisé, - on résout un système de co-krigeage factoriel qui correspond à ces matrices pour en déduire un opérateur spatial à appliquer aux deux cartographies 2D qui leur correspondent, afin d'en déduire une cartographie 2D qui correspond à leur 30 composante commune, l'ensemble des cartographies 2D ainsi obtenues permettant de reconstituer le cube sismique commun. 12. Procédé selon la revendication 11, caractérisé en ce que pour déterminer le système de co-krigeage à résoudre pour déduire l'opérateur spatial à appliquer à deux cartographies 2D, on applique aux données d'au moins l'une des deux cartographies, préalablement à la détermination des 5 co-variances et des cross-variances, un masque destinée à annuler les données de ladite cartographie sur une certaine zone. 13. Procédé pour le filtrage d'un ensemble de cubes sismiques correspondant, pour au moins une campagne d'acquisition sismique dans une zone donnée, à différentes valeurs de déport source/récepteur et/ou 10 d'angle d'incidence, caractérisé en ce qu'on détermine au moins un cube sismique qui est une estimation de la composante commune à au moins deux cubes sismiques correspondant à des valeurs successives de déport source/récepteur et/ou d'angle d'incidence. 14. Procédé selon la revendication 13 caractérisé en ce que ce 15 cube commun est associé, dans un nouvel ensemble de cubes sismiques, à une valeur de déport et/ou d'angle d'incidence comprise entre l'une et l'autre des deux valeurs successives auxquelles correspondent les deux cubes à partir desquels une estimation de la composante commune est déterminée et on met en oeuvre sur le nouvel ensemble de cubes sismiques 20 ainsi constitué, un traitement pour la détermination d'au moins une information relative à la géologie du sous-sol de la zone sur laquelle la campagne d'acquisition sismique a été mise en oeuvre. 15. Procédé selon la revendication 14 caractérisé en ce qu'on met en oeuvre sur ledit nouvel ensemble de cubes sismiques un traitement pour 25 la détermination d'attributs AVO ou AVA. 16. Procédé selon l'une des revendications 14 ou 15, caractérisé en ce qu'on traite l'ensemble des cubes initiaux pour déterminer pour chaque couple de valeurs de déport et/ou d'angle d'incidence successives un cube sismique qui est une estimation de la composante commune aux 30 deux cubes sismiques correspondant à ces deux valeurs de déport et/ou d'angle d'incidence, l'ensemble des cubes communs ainsi obtenus constituant un nouvel ensemble de cubes sismiques sur lequel on met en oeuvre un traitement pour la détermination d'au moins une information relative à la géologie du sous-sol. 17. Procédé selon l'une des revendications 14 ou 15, caractérisé en ce qu'on sélectionne sur l'ensemble des cubes initiaux, des cubes 5 sismiques apparemment bruités et on détermine pour chacun de ces cubes au moins un cube sismique qui est une estimation de la composante commune à ce cube et à un cube correspondant à une valeur de déport et/ou d'angle d'incidence immédiatement voisine, les cubes ainsi déterminés, ainsi que les autres cubes initiaux constituant un nouvel 10 ensemble de cubes sismiques sur lequel on met en oeuvre un traitement pour la détermination d'au moins une information relative à la géologie du sous-sol. 18. Procédé pour le traitement deux cubes sismiques issus de deux campagnes d'acquisition réalisées pour une même zone au sol à deux 15 moments différents, caractérisé en ce que pour chaque valeur de déport et/ou d'angle d'incidence, on détermine une estimation d'un cube correspondant à la composante commune aux deux cubes qui correspondent à ladite valeur de déport et/ou d'angle d'incidence, ainsi qu'une estimation de cubes résiduels et en ce qu'on met en oeuvre sur les 20 ensembles de cubes résiduels ainsi obtenus pour chacune des deux campagnes d'acquisition un procédé de filtrage selon l'une des revendications 13 à 17. 19. Procédé selon la revendication 18, caractérisé en ce qu'on détermine un cube qui est la somme, sur l'ensemble des valeurs de déport 25 et/ou d'angle d'incidence, des différences des cubes résiduels filtrés correspondant à l'une et à l'autre des deux campagnes d'acquisition. 20. Procédé selon la revendication 18, caractérisé en ce que pour l'une au moins des campagnes d'acquisition, on détermine un cube qui est la somme, sur l'ensemble des valeurs de déport et/ou d'angle d'incidence, 30 d'une part des cubes qui correspondent aux composantes communes aux deux campagnes d'acquisition et d'autre part des cubes résiduels filtrés correspondant à ladite campagne d'acquisition. 21. Procédé selon la revendication 18, caractérisé en ce que pour l'une au moins des campagnes d'acquisition, on détermine pour chaque valeur de départ et/ou d'angle d'incidence, un cube filtré qui est la somme du cube qui correspond à la composante commune aux deux campagnes 5 d'acquisition pour ladite valeur de déport et/ou d'angle d'incidence et d'un cube résiduel fonction d'une ou plusieurs cubes résiduels filtrés déterminé(s) pour la valeur de déport et/ou d'angle d'incidence et au moins une valeur de déport et/ou d'angle d'incidence voisine.