FR 2772138 A1 19990611 FR 9715742 A 19971210 La présente invention concerne une méthode statistique multivariable pour caractériser des images que l'on a formées d'un milieu complexe tel que le sous-sol par exemple, de manière à faire apparaître les relations spatiales entre les éléments de la structure du milieu. La méthode selon l'invention trouve des applications dans de nombreux domaines notamment pour caractériser au mieux l'imagerie sismique. La méthode met en oeuvre une technique descriptive d'analyse de données dite de la contiguïté spatiale qui s'applique à des données spatiales ou reliées par un graphe de contiguïté telles que des données sismiques. Les données analysées sont constituées d'une famille d'images telles que par exemple, les amplitudes des horizons d'un cube sismique 3D (Fig. la à le). Ces images présentent des structures spatiales gigognes à petite et grande échelles (micro et macro-structures) potentiellement liées à des phénomènes géologiques tels que des chenalisations, des réseaux de failles etc. ou à la propagation d'ondes avec formation de multiples par exemple, ou à des artefacts d'acquisition ou de traitement. Du bruit aléatoire peut être mélangé à ces structures. On considère que les données à analyser sont spatiales du fait qu'elles se présentent sous l'aspect d'images formées de pixels répartis sur un plan, caractérisées par un certain nombre de mesures telles que des amplitudes relevées sur des portions de traces sismiques, des coefficients de réflexion sur des discontinuités, etc. que l'on considère en tant que variables statistiques. Il est alors primordial, pour rendre compte de l'ensemble des structures inhérentes à ces données, de prendre en considération à la fois les relations de contiguïté spatiale des données et leurs caractéristiques (géo)physiques. État de la technique La technique connue d'analyse de la contiguïté spatiale (ACS) a fait l'objet de nombreuses publications notamment par L. Lebart. 1969, Analyse statistique de la contiguïté, Pub. Ins. Stat., Paris VIII, 81-112. Sous l'appellation d'analyse des proximités spatiales, la technique en cause a été appliquée notamment pour le filtrage de données sismiques comme décrit par: - Royer, J. J., 1984, Proximity Analysis: a Method for Geodata Processing, in Sciences de la Terre, n 20, Proc. of the Int. Coll.: Computers in Earth Sciences for Natural Ressources Characterization, April 9-13, Nancy, France; ou par - Faraj, A.,1994, Application of Spatial Contiguity Analysis to Seismic Data Filtering. In SEG - 64th Annual International SEG Meeting, Los Angeles, Oct. 23-28 1994, Expanded abstracts, vol.1, Paper SP5.7, 15841587. Elle consiste essentiellement à calculer une famille de composantes combinaisons linéaires des variables initiales qui minimisent le coefficient de contiguïté dit de Geary définissant le rapport de la variabilité spatiale à la variance. Les composantes (appelées composantes spatiales) correspondent aux vecteurs propres de la matrice C-' F(h) o C est la matrice de variancecovariance des données initiales et F(h) celle des variogrammesvariogrammes croisés à la distance spatiale h. Les composantes spatiales sont habituellement ordonnées par ordre décroissant des valeurs propres de cette matrice. On considère généralement que les premières composantes (associées aux faibles valeurs propres) et appelées composantes régionales, mettent en évidence les structures spatiales à grande échelle. Les dernières composantes spatiales (associées aux fortes valeurs propres) appelées composantes locales, traduisent les structures spatiales à faible échelle. On considère également que l'information "portée" par ces différentes composantes est mesurée par les valeurs propres associées. La méthodologie suivie pour l'analyse de la contiguïté spatiale est calquée directement sur celle plus connue d'analyse en composantes principales (ACP). Or on a vérifié, que les valeurs propres de la matrice C-' F(h) sont de piètres mesures de la variance locale pour caractériser et ordonner les composantes spatiales. En effet, lors de certaines analyses effectuées, les composantes associées aux plus grandes valeurs propres sont complètement dépourvues de sens (bruit aléatoire par exemple) alors que celles correspondant à des valeurs propres intermédiaires semblent mieux rendre compte (au moins visuellement) de l'information aussi bien statistique que spatiale des variables initiales. En réalité, bien que calculées à partir des matrices C et F(h) qui intègrent les interdépendances à la fois statistiques et spatiales des données et quoique mesurant les rapports de contiguïté des composantes factorielles, les valeurs propres représentent chacune une information isolée (rapport de la variance locale sur la variance globale) spécifique à la composante spatiale. La valeur propre est un critère qui n'est bon que pour mesurer la variabilité spatiale de la composante spatiale. À l'inverse de l'analyse en composantes principales (ACP), par exemple, o les valeurs propres (celles de C) représentent la part de la variance totale des données expliquée par la composante, la somme des valeurs propres obtenue dans le cadre de l'ACS, égale à tr[C-' F(h)] n'a pas de signification. Par contre, obtenir tr[F(h)] est bien préférable car cela représente la somme des variances locales des variables initiales. De ce fait, il apparaît nécessaire de définir de nouveaux critères pour quantifier les informations statistiques et spatiales portées par les composantes spatiales de l'ACS. Définition de la méthode selon l'invention Il s'agit d'une méthode statistique multivariable pour analyser des données associées (directement ou indirectement) à des éléments d'images traduisant des propriétés physiques d'un milieu complexe tel que le sous-sol par exemple, ces images étant obtenues par une exploration du milieu (au moyen d'ondes sismiques par exemple), de manière à faire apparaître les relations spatiales entre ces éléments. La méthode est caractérisée en ce qu'elle comporte: - une analyse spatiale des données pour mettre en évidence les propriétés spatiales des dits événements dans une ou plusieurs directions; l'application aux données d'une technique d'analyse de la contiguïté spatiale pour trouver leur meilleure décomposition en structures spatiales composantes; et - le filtrage des structures obtenues par cette décomposition, pour extraire les plus pertinentes. La méthode comporte par exemple: - une identification des structures spatiales des données; - une séparation de ces structures spatiales avec élimination d'éventuelles redondances; - la formation, à partir des images initiales, d'images synthétiques ou composantes spatiales qui mettent en évidence les structures spatiales des données; - l'établissement d'une typologie des images initiales en fonction des structures spatiales mises en évidence par les composantes spatiales; et - le filtrage des images initiales pour enlever le bruit et sélectionner une ou plusieurs structures spatiales identifiées. On caractérise la structure spatiale des données par exemple par une analyse monovariable ou bivariable de différents variogrammes des éléments d'image. La sélection des structures spatiales composantes est effectuée par exemple en déterminant leurs contributions respectives à la variabilité spatiale des données et à la variance des données, cette sélection pouvant être effectuée graphiquement. La méthode permet non seulement une identification des structures spatiales mais aussi une quantification de l'information filtrée en termes de variance et de variabilité spatiale. La méthode peut s'appliquer d'une façon générale à l'analyse de toutes données spatiales à répartition régulière ou irrégulière. D'autres caractéristiques et avantages de la méthode selon l'invention vont apparaître à la lecture de la description ci-après qui comporte deux parties, la première portant sur les grandes lignes de la méthode, la deuxième, sur l'appareillage statistique mis en ceuvre pour réaliser les différentes étapes de l'analyse, et en se référant aux dessins annexés o: - les Fig.la à If montrent les images de six horizons Z1 à Z6 d'un bloc sismique; - les Fig.2a à 2h montrent les courbes de variogramme de six images Z1 à Z6 du bloc choisies dans huit directions; - les Fig.3a, 3b montrent deux diagrammes des contributions des composantes spatiales CS respectivement à la variabilité spatiale (VS) et à la variance (V) des données; - les Fig.4a à 4f montrent respectivement les composantes spatiales (Fl-F6) de différentes structures; - les Fig.5a à 5e montrent des représentations des variables initiales sur des cercles des corrélations associé à cinq plans factoriels (de PF(1,2) à PF(5, 6)); les Fig.6a à 6e montrent des représentations des variables initiales sur des cercles des contributions associé au plan factoriel (PF(1,2) à PF(5,6)); - les Fig.7a à 7f montrent sur les images initiales ZI à Z6, la structure composante F5-F6 correspondant à la structure spatiale Est-Ouest; er - les Fig.8a à 8f montrent le bruit aléatoire sur les images ZI à Z6 initiales. I Description détaillée de la méthode selon l'invention: La méthode s'applique aux données associées aux différents éléments d'images que l'on a obtenues par une exploration sismique du sous-sol par exemple et qu'il s'agit d'analyser pour faire mieux ressortir les rapports de contiguïté spatiale qu'il y a entre eux. Ces images sont par exemple, extraites à partir des horizons ou des sections sismiques d'un "cube sismique" 3D. Elle comporte essentiellement trois étapes successives. I-1) Analyse spatiale Appelée étude variographique, cette étape préliminaire consiste à faire une analyse spatiale des données, soit monovariable (une variable à la fois), soit bivariable (i.e. en prenant les variables deux à deux). C'est au cours de cette étape que l'on met en évidence les propriétés spatiales des variables initiales dans les différentes directions: leur organisation spatiale à différentes échelles, l'importance respective du bruit aléatoire dans les différentes images l'existence éventuelle de structures périodiques etc. L'analyse variographique monovariable des données (ou variables initiales) consiste à analyser les courbes des variogrammes et des autocorrélations, les variogrammes et autocorrélogrammes surfaciques associés à chacune des variables initiales et à les comparer les unes aux autres pour dégager des comportements spatiaux particuliers. L'analyse variographique bivariable des données consiste à analyser les courbes des variogrammes croisés et des intercorrélations, les variogrammes croisés et intercorrélogrammes surfaciques associés aux variables initiales prises deux à deux, pour dégager les corrélations spatiales des variables les unes avec les autres. Cette étude variographique (mono ou bivariable) est prioritairement une phase descriptive o il est moins question de modéliser les courbes de variogrammes que de classifier les variables initiales en groupes homogènes de comportement spatial. Les différents outils variographiques cités: variogrammes, variogrammes croisés, autocorrélations et intercorrélations directionnels ou surfaciques, sont calculés et visualisés pour toute la famille des variables analysées. Cette étude variographique peut mettre en évidence des distances spatiales pertinentes. I-2) Application de l'analyse de contiguïité spatiale (ACS) Cette préparation préalable étant effectuée, on effectue dans un deuxième temps une ou plusieurs analyses de la contiguïté spatiale sur les données, dans le but de trouver leur meilleure décomposition en structures spatiales orthogonales. C'est au cours de cette étape que l'utilisateur a la possibilité d'exploiter les outils d'aide à l'interprétation des données et des résultats que nous mettons à sa disposition. Une ou plusieurs analyses sont appliquées aux données. L'utilisateur est orienté pour la mise en oeuvre de ces analyses, par les distances spatiales dégagées lors de l'analyse variographique préalable. Diagrammes La pertinence des composantes spatiales issues d'une analyse donnée est indiquée par deux diagrammes: celui des contributions des composantes spatiales à la variabilité spatiale des données et celui des contributions des composantes spatiales à la variance des données. Ces deux diagrammes montrent la distribution des deux critères en question (contributions à la variabilité spatiale définies plus loin (cf. relation (17)) et contributions à la variance: cf. relation (19)) pour l'ensemble des x composantes spatiales d'une analyse. Ces diagrammes permettent de sélectionner les composantes les plus significatives au sens de la variabilité spatiale ou de la variance. Les composantes spatiales sont rangées suivant les deux ordres suivants 1- par ordre décroissant de leur contribution à la variabilité spatiale des données (cf. relation (20)). On désigne par composantes spatiales locales, les premières composantes de cet ordre car elles rendent compte de la structure des données à petite échelle. 2- par ordre décroissant de leur contribution à la variance des données (cf. relation (21)). On désigne par composantes spatiales régionales, les premières composantes de cet ordre, car elles rendent compte de la structure des données à grande échelle. Graphiques Pour analyser les variables initiales en relation avec les composantes spatiales, deux représentations des variables sont possibles et complémentaires 1- Celles des contributions signées des composantes spatiales à la variabilité spatiale des variables initiales. (cf. relation (15)). De telles valeurs, comprises entre - 1 et +1, donnent lieu à des cercles (appelées cercles des contributions à la variabilité spatiale) à l'intérieur desquels les variables initiales sont représentées par des points. 2- Celles des contributions signées des composantes spatiales à la variance des variables initiales (cf. relation (18)). De telles valeurs, comprises entre -1 et +1, donnent lieu à des cercles (appelés cercles des contributions à la variance) à l'intérieur desquels les variables initiales sont représentées par des points. Ces deux représentations graphiques complémentaires permettent d'interpréter les variables initiales en relation avec les structures spatiales mises en évidence par les composantes. On peut aussi utiliser ces représentations graphiques (cercles des corrélations et des contributions) pour placer à côté des variables actives, des variables supplémentaires, c'est-à-dire des variables qui, contrairement à elles, n'entrent pas en jeu pour l'identification des composantes spatiales, ceci grâce à leurs contributions signées ( cf. relations (15) et (18)). Il peut s'agir par exemple d'un attribut sismique que l'utilisateur désire positionner durant l'analyse, par rapport aux attributs sismiques qui interviennent dans l'analyse. 1-3 Filtrage des variables L'étape finale est celle du filtrage des variables initiales ou des variables non prises en compte pour l'identification des composantes spatiales, o l'utilisateur a la possibilité d'extraire la ou les structures qu'il juge pertinentes à la lumière des outils d'aide et de son expertise dans le domaine. Cette étape consiste à enlever le bruit et à éliminer ou à conserver une ou plusieurs structures spatiales identifiées. Le résultat est une image initiale débruitée ou ne contenant que la structure (locale ou régionale) jugée comme géologiquement pertinente par l'utilisateur. Un tel filtrage n'est pas unique: il dépend des composantes locales ou régionales décidées comme telles par l'utilisateur. Ce filtrage peut être appliqué aussi bien aux variables initiales qu'à des variables n'ayant pas participé aux calculs. L'utilisateur peut positionner ces variables par rapport aux structures spatiales mises en évidence par élimination ou extraction de ces structures. I-4 Quantification des structures filtrées La méthode selon l'invention offre à l'utilisateur une quantification des structures spatiales éliminées ou extraites des variables (aussi bien actives que supplémentaires) à la fois en termes de variabilité spatiale que de variance, car les mesures obtenues par les relations (22) ou (23) sont des pourcentages. II Description analytique des movens d'analyse statistique mis en oeuvre Rappel: Pour la bonne compréhension de ce qui va suivre, on rappelle tout d'abord les notions et notations connues des hommes de l'art. On considère n objets (i = 1 à n) localisés aux points xi sur un plan de position et décrits par J variables Zj (j = 1 à J) de sorte que Zji soit la mesure de la variable Z' sur l'objet i. Z = [Zi] est la matrice des données. En analyse des données classique on considère de telles données en tant que variables (au sens statistique) définies sur un ensemble d'individus. Les n mesures Zi sont des réalisations indépendantes de la variable aléatoire Z'. Cette représentation ôte tout caractère spatial aux données. Il est d'usage, en géostatistique, de représenter de telles données par un ensemble de variables régionalisées (V.R.) considérées comme des réalisations d'une famille Z(x)={Zi (x); j = 1,..., J} de fonctions aléatoires du point x. (i.e. Z(x) est une fonction aléatoire vectorielle à valeurs dans R'). Ainsi zj(x1) - que l'on notera Zi - est une réalisation de la fonction aléatoire ZJ(x) du point x. De cette façon, les n mesures multidimensionnelles (Z'(x), Z2(xi),..., Zl(xi)) sont des réalisations de n vecteurs aléatoires différents situés aux points xi du plan de position. Le terme "régionalisé" a été proposé pour qualifier un phénomène se déployant dans l'espace et y présentant une certaine structure. Une variable Il régionalisée est considérée comme étant une fonction irrégulière du point x du fait quelle révèle deux aspects contradictoires e l'un, structuré, associé à une organisation spatiale à plus ou moins grande échelle, * l'autre, local, faisant apparaître les irrégularités à plus petite échelle, imprévisibles d'un point à un autre. Entre ces deux aspects extrêmes du signal, d'autres structures existent décrivant des organisations spatiales des données à des échelles intermédiaires. C'est l'ensemble de toutes ces structures que l'on propose de décrire, dans un premier temps, grâce aux outils géostatistiques. Ces structures, vu l'aspect multidimensionnel des données, peuvent être redondantes dans la mesure o elles seraient communes aux J variables initiales. Une fois mises en évidence, elles peuvent servir pour effectuer une typologie statistique et spatiale de ces variables. II-1 Mise en oeuvre de l'analyse spatiale On se fixe une direction sur le plan de position des données et on considère une distance dans cette direction. Nous désignerons par h à la fois cette distance (qui est un scalaire) et le vecteur (de dimensions 2) défini par la direction et la longueur h. On notera AZ(h) = Z(x + h)- Z(x) (1) l'accroissement multivariable de Z entre les points x et x+h séparés par la distance h dans la direction préalablement fixée. Z est supposée stationnaire d'ordre 2; ce qui signifie que l'accroissement AZ(h) est de moyenne nulle - i.e. E[AZ(h)] = 0 - et de variance E[AZ(h)2] constante ne dépendant que de h. Dans ces conditions, on définit la matrice des variogrammes-variogrammes croisés Fzz(h) comme étant la matrice de variances- covariances des accroissements AZ(h) 1I Fzz (h) = E[AZ(h)'. AZ(h)] = C,zz, (2) o Cz(h)àz(h) est la matrice de variances- covariances associées aux colonnes de AZ(h). Fzz(h) est une matrice de dimension JxJ de terme général Ty (h) = E[(Z (x + h) - Zj(x)).(Zj (x + h) -Z] (x))] (3) La direction sur le plan de position des données étant fixée, le terme diagonal de Fzz(h) désigne alors le semi- variogramme de Zj dont la valeur mesure la variabilité spatiale de Zj à la distance h. Le variogramme expérimental est utilisé en tant que courbe dépendant de la distance h. Il pourrait encore s'écrire sous la forme Yj(h) =Yjj(h) (h) y- (Z.-Zi.) (4) y/)2m(h)dt(l-.)h o m(h) est le nombre de paires formées d'individus spatialement distants de h. Le terme Yl'(h) désigne, d'une certaine façon, la covariance locale entre les variables j et j' si nous considérons comme graphe de contiguïté associé celui qui relie les paires formées de points à la distance h l'un de l'autre et disposés dans la direction du vecteur h. De même yj(h), terme de variance locale de la variable j. représente une restriction de la variance de cette variable uniquement calculée pour les paires de points distants de h. Elle représente les fluctuations spatiales associées à cette distance. Considérée comme une courbe dépendant de h, la forme de y'y (Fig.2) renseigne sur le comportement spatial de la variable Zj pour les différentes valeurs de h. Le but d'une étude variographique classique est de modéliser les variogrammes expérimentaux à l'aide de modèles théoriques de base sphériques, ou exponentiels. Elle a aussi l'avantage, d'un point de vue pratique, de mettre en évidence les structures spatiales imbriquées d'une variable régionalisée. Un exemple d'une telle étude variographique est décrit par exemple par - Isaaks, E. H., Srivastava, R. M.,1989, Applied Geostatistics. Oxford University Press, Oxford. Dans le cadre de la présente méthode, on cherche plutôt à mettre en évidence des groupes de comportement spatial en comparant les variogrammes des J variables les uns aux autres dans le but de révéler des familles de structure spatiale et de déterminer la (ou les) distance(s) convenable(s) pour la mise en oeuvre de l'ACS. II-2. Mise en oeuvre de l'analyse de contiguïté spatiale (ACS) L'analyse de la contiguïté spatiale (ACS) permet, dans un deuxième temps, de mettre en évidence les structures spatiales, de les interpréter et les filtrer éventuellement. On va chercher en définitive à "capturer" de telles structures par les composantes factorielles de l'ACS. La distance spatiale h entre individus ou objets étant fixée dans une direction donnée, il s'agit de déterminer une variable f dépendant de h f = Zu = EujZi (5) qui est une combinaison linéaire des variables Zj initiales dont le rapport de contiguïté spatiale pour le pas h est maximal. Cette variable est celle qui réalise le maximum de l'expression y1(h) uF(7h)u, (6) : f U CzzU i.e. qui possède, à la fois, une variabilité spatiale maximale et une variance minimale. Dans cette expression, Czz désigne la matrice des variances-covariances des données. Il existe J variables f', f,, f,. f solutions de ce problème (appelées composantes spatiales) correspondant aux vecteurs propres u,..., u',..., u' de la matrice CFz7(h) - o u, u=(u,. uj)te R' - ordonnés par valeurs propres croissantes: xi cème composante factorielle. C'est dans cet ordre que sont usuellement présentées les composantes d'une analyse factorielle classique, mettant en avant l'information principale des données (au sens du critère optimisé). De cette manière, la plus grande part de la dite information est portée par les premières composantes. Ceci établit un parallèle entre les composantes de l'ACS (structures à grande échelle sur les premières composantes et microstructures ou bruit aléatoire sur les dernières) et celles, par exemple, de l'analyse en composantes principales (plus grande variance sur les premières composantes et bruit aléatoire sur les dernières). Propriétés statistiques et spatiales des composantes spatiales de i'ACS Le coefficient de Geary I (h) combine deux critères: l'un spatial au o:(f) numérateur et l'autre algébrique (ou statistique) au dénominateur permettant ainsi de calculer des composantes factorielles de variance maximale et de variabilité spatiale minimale. Ces composantes ont la propriété de former une base à la fois statistiquement et spatialement orthogonale. cette deuxième propriété ne semble pas avoir été exploitée par les différents auteurs ayant travaillé sur le sujet. Et c'est justement cette propriété qui permet de construire des outils d'interprétation des résultats. Soit F = [fa t= le tableau rectangulaire des mesures des J composantes factorielles de l'ACS (en colonnes) sur les n individus (en lignes). Les variables initiales étant centrées, les composantes factorielles le sont aussi. On montre que la matrice CFF de variances-covariances de F s'écrit: CFF= F' F = IJ (8) o I, est la matrice identité JxJ. La matrice des variances-covariances locales (variogrammes-variogrammes croisés) associée aux composantes factorielles s'écrit FFF(h) = A (9) o A est la matrice diagonale des valeurs propres (Xk) oc=l.. Autrement dit, le coefficient de corrélations de deux composantes fa et f s'écrit: c(f a f3){ _ si (10) c r(f f) O L0 sinon, et leur variogramme croisé Y (h) = {Xk si oc=11 r"? =0 sinon Les composantes spatiales sont donc orthonormées. De plus leurs variances locales - valeurs de leurs variogrammes pour la distance h - sont égales aux valeurs propres de l'analyse. Leurs covariances locales - valeurs de leurs variogrammes croisés pour la distance h - sont nulles deux à deux. Elles sont donc spatialement indépendantes. De ceci, on peut déduire des relations intéressantes concernant la décomposition de la matrice des variances-covariances locales variogrammes- variogrammes croisés - des variables initiales. On peut montrer que son terme général s'écrit sous la forme /j(h) = E cor(Zj, fa) cor(Zj,fa) X (12) a et plus particulièrement le variogramme de la variable Zj yj(h) = E cor2(Zj,fa);Ma (13) a Et comme ' cor2(ZJ,fa) = 1, la variance locale (valeur du variogramme pour la distance h) d'une variable régionalisée Zi s'écrit comme une moyenne pondérée par cor2(Zi,fc) des variances locales (valeurs des variogrammes pour la distance h) des composantes spatiales f de l'ACS. De façon générale, la matrice des variances-covariances locales - variogrammes-variogrammes croisés -, compte tenu de la relation 13, se met sous la forme Fzz (h) = C ACFz (14) o A est la matrice diagonale des valeurs propres de C'-,zFz(h) et CFZ celle des covariances entre les composantes factorielles fa (en lignes) et les variables initiales Zj (en colonnes). II-3. Outils d'aide à l'interprétation des résultats d'une ACS Dans ce qui suit, on suppose que l'analyse ACS a été appliquée pour un vecteur h fixé du plan de position des données. Ainsi tous les résultats (valeurs propres Xa, composantes factorielles fa, etc.) de l'ACS dépendent de ce vecteur h. L'écriture de yj(h) dans (13) montre que l'expression que nous définissons comme suit cnt(ZJ,fa)= cor(Zj,f). lcor(z' fa)l (15) 1/j.(h) varie entre -1I et 1. Elle mesure la contribution signée de la cème composante factorielle à la structure (ou variabilité) spatiale de la variable Zj. Une telle définition permet d'avoir des valeurs aussi bien négatives que positives de la contribution, lui conférant ainsi un rôle identique à celui de la corrélation. De cette façon, il est profitable de représenter les variables initiales sur des cercles que l'on appelle cercles de contributions (cercle de rayon 1, Fig.3) dont l'utilisation estidentique à celle des cercles de corrélations (Fig.2) employés dans l'analyse en composantes principales. On notera que, pour toute variable Zj Ylcnt(ZJ,f C =l (16) tu-I De fait, si Icnt(Zif,)l:1, on dira que la composante fa contribue fortement à la variance locale (ou variabilité spatiale) de la variable Zi'. Le cercle des contributions spatiales permet de faire la typologie des variables initiales avec une meilleure précision. Par ailleurs, on définit cor 2(Zi, faf)Xa XICft(ZJ, fa4yh.fj =L i(h) j Fi (17) Z ij( h). 1=l Cette expression qui varie entre O et 1, mesure la contribution de la cème composante factorielle à la structure (ou variabilité) spatiale des données. Cette contribution est d'autant plus élevée que cnt(f') I=1. Au dénominateur du second terme de la relation (17) ô'yj(h)=trace[Frzz(h)] représente la variance locale totale l=1 des données. Le numérateur du troisième terme de l'expression (17) montre que la contribution d'une composante spatiale à la variabilité spatiale des données est une somme des variances locales des variables initiales pondérées par les contributions absolues Icnt(ZJ,f")I de la composante à la variabilité spatiale des variables. Ainsi la structure spatiale d'une variable est d'autant plus marquante que la valeur de cnt(Zi,f a") est élevée en valeur absolue. Le premier indice cnt(Zj,fa) est utile, d'une part, pour l'interprétation des composantes factorielles en fonction des variables initiales et, d'autre part, pour la typologie de ces dernières. Une variable Z' qui, d'un point de vue spatial, est la plus ressemblante à la composante fa, est celle pour laquelle la valeur de Icnt(ZJ,fa)I est proche de 1. Un tel résultat est visualisé - comme nous venons de le voir - par l'intermédiaire du cercle des contributions en représentant les variables initiales sur un plan factoriel (f',ff) par leurs coordonnées cnt(Z-',f') et cnt(Zj,f). Le deuxième indice cnt(fa) est une mesure globale de la variabilité spatiale que chaque composante restitue. C'est un descripteur de l'information spatiale portée par les structures mises en évidence. Ce critère est utile comme on le verra plus loin, pour sélectionner les composantes les plus significatives. C'est sur la base de ce critère, et non plus les valeurs propres, que l'on ordonne les composantes factorielles par contributions spatiales décroissantes. Ces indices ne sont cependant pas suffisants pour décrire globalement les données. Ils ne tiennent pas compte de la variance des données qui est indispensable pour mesurer la part de l'information statistique restituée par les composantes. Pour cette raison on retient l'indice suivant: cntV(Zj,fr) =cor(Zi,f).Icor( Zi, fû (18) qui varie entre -l et 1 et mesure la contribution signée de la coWme composante factorielle à la variance de la variable Zj. La définition d'un tel indice découle du fait que Ecor2(Zij,f)= I car les variables initiales sont centrées et réduites. oE=l Plus généralement, l'indice _cor2(Zj,f).(2 _IcntV(ZJ,f'.c72 cntV(f')= '=. (19) variant entre O et 1, mesure la contribution de la ctème composante factorielle à la variance des données. De cette façon on a, pour chaque composante, la mesure de la part de variance des données qu'elle explique. Dans la relation (19), le dénominateur la =trace[Czz] du second terme représente la variance totale des données. À défaut des valeurs propres directement disponibles en ACP, les valeurs de cntV(f) renseignent sur la part de variance portée par les composantes spatiales issues de l'ACS. Plus cntV(f) est élevée plus fa pourrait être assimilable à une composante principale de l'ACP (forte variance expliquée). Dans la pratique il y a une forte corrélation entre la première composante g' de I'ACP et la composante de Gr2 l'ACS pour laquelle cntV(g')=oG,. On montre d'ailleurs, que cnt(ga) = cr=l pourcentage de l'inertie expliquée par la xième composante principale ga de I'ACP. Par ailleurs, plus cnt(fa) est élevée plus la composante traduit une forte variabilité spatiale des données. Il est utile, par conséquence. de ranger les composantes f de l'ACS, d'une part, par ordre décroissant des cnt(f t), soit: cnt( f""')> cnt(f"L')>... > cnt(f "') (20) et, d'autre part, par ordre décroissant des cntV(f '), soit cntV(f"') > cntV( f':2)>...> cntV( f J') (21) On voit sur les Fig.5 et 6 les diagrammes respectifs des contributions cnt(f) o0 et cntV(f ') des composantes à la variabilité spatiale des données et à la variance des données. Les q premières composantes retenues dans la relation (20) sont appelées composantes locales car elles rendent compte de la structure des données à petite échelle, alors que les p premières composantes retenues dans la relation (21) sont appelées composantes régionales car elles traduisent la structure des données à grande échelle. Certaines composantes peuvent être à la fois régionales et locales. Les composantes fa rangées en dernière position à la fois dans les relations (20) et (21) i.e. pour lesquelles on a à la fois cnt(fa) 0 et cntV(fa) 0), sont peu intéressantes et on les assimile à du bruit aléatoire. La part de la variabilité spatiale des données expliquée par les q premières composantes locales est: Sio:ailc = Ecnt(f'k') (22) La part de la variance des données expliquée par les p premières composantes régionales est regionale= uCntV(f,) (23) La condition nécessaire pour calculer les contributions cntV(f') pour une famille de composantes factorielles est que celle-ci soit orthogonale. On ne peut pas réciproquement intégrer la contribution à la variabilité spatiale cnt(ga) comme outil d'aide à l'interprétation des composantes de l'analyse en composantes principales ACP. Le calcul de cnt(g') nécessite que les g* soient spatialement orthonormées (i.e. leurs covariances locales sont nulles deux à deux). Ceci est le cas des composantes de l'ACS mais non de celles de l'ACP. Les composantes factorielles, à la fois statistiquement orthogonales et spatialement indépendantes, issues de la méthode selon l'invention, permettent de rendre compte des structures spatiales intrinsèques aux données. Une fois mises en évidence, ces structures pourront être utilisées pour décrire finement les images sismiques. 11-4. Filtrage des variables initiales On pourra extraire de toute variable initiale Z1 respectivement une (ou plusieurs) parties régionales ZJ gionale(h) qui traduisent des phénomènes spatiaux à grande échelle et une (ou plusieurs) parties locales ZJocai(h) qui traduisent des phénomènes spatiaux à petite échelle. Ces nouvelles variables s'écrivent Zrgional (h) = r k (24) comp. régionales k et Zio aJ(h) = f (25) comup, locales k Les composantes régionales et locales intervenant dans les deux relations (24) et (25) sont sélectionnées par l'utilisateur sur la base des diagrammes des contributions des composantes spatiales à la variance et à la variabilité spatiale des i données. Exemple d'application de l'ASC à un jeu de données Les six images ZI,..., Z6 sur les Fig.la à If représentent les amplitudes des six horizons d'un bloc sismique. Trois structures sont visibles sur ces images dont certaines sont emboîtées les unes dans les autres: une structure d'orientation Est- Ouest commune aux images Z1, Z2 et Z3, une structure d'orientation NordSud commune aux images Z4. Z et Z6, et du bruit aléatoire mélangé avec les deux structures Nord- Est et Est-Ouest principalement dans les images Z2, Z4 et Z6. Les structures Est-Ouest et Nord-Sud sont des organisations spatiales à grande échelle alors que le bruit aléatoire apparaît plutôt à plus petite échelle. Par ailleurs, bien que structurellement ressemblantes les unes avec les autres, les images ZI, Z2 et Z3 (respectivement Z4, Z5 et Z6) sont faiblement corrélées les unes avec les autres à cause des décalages apparents des lentilles composant chacune des images comme le montre le tableau des corrélations ci-dessous. Cette faible corrélation rend la séparation de ces structures difficile -voire impossible- par une méthode multivariable classique (telle que l'ACP) basée sur les seules corrélations entre les variables. Zi Z2 Z3 Z4 Z5 Z6 _ Z1 1.00 -0.26 0.11 -0.1 1 0.05 -0.02 _ _ _Z2 -0.26 1.00 0.14 0.27 0.07 0.41 Z3 0.11 0.14 1.00 -0.05 0.13 0.00 Z4 -0.11 0.27 -0.05 1.00 -0.15 0.17 Z5 0.05 0.07 1 0.13 -0.15 1.00 I 1.00 Z6 -0.02 0.41 0.00 0.17 0.02 1.00 Les variogrammes des six images dans huit directions (Nord-Sud, N220, N450, N670, Est-Ouest, N67E, N45E et N22E) sont représentés sur les Fig. 2a à 2h. Certaines directions (N450 et N45E) ne permettent pas de distinguer les six images les unes par rapport aux autres, alors que les directions Nord-Est et Est-Ouest séparent mieux les courbes en montrant deux familles de comportements spatiaux homogènes et bien séparés. On voit apparaître les périodicités des événements horizontaux (images Z1i.Z2 et Z3) et verticaux (Z4, Z5 et Z6) sur les variogrammes correspondants. Ces variogrammes donnent même une mesure de cette périodicité (- mètres). On peut aussi noter que le saut au voisinage de l'origine (tous variogrammes confondus) est différent d'une image à l'autre. Celles-ci sont ordonnées par ordre croissant de la valeur à l'origine comme suit: ZI, Z5, Z3, Z2, Z4 et Z6. C'est précisément dans cet ordre que les images sont rangées des moins bruitées aux plus bruitées. Si l'on veut appliquer l'ACS dans le but de séparer le bruit spatial des deux autres structures à grandes échelles, sans distinction entre ces deux structures, c'est la distance près de l'origine qu'il faudrait choisir. Par contre, si ce sont ces deux structures à grande échelle que l'on désire séparer, il faudrait appliquer l'ACS dans la direction Nord-Sud ou Est-Ouest pour une distance de 18 mètres qui correspond à la distance pour laquelle la séparation entre les deux familles d'images est optimale. L'ACS sera appliquée pour la distance de 18 mètres dans la direction Nord-Sud. On peut noter sur les Fig.3a, 3b qui donnent les valeurs propres, et les diagrammes de contributions à la variabilité spatiale et à la variance des 6 composantes spatiales, que les composantes présentent des valeurs de contributions à la variance à peu près équivalentes. C'est alors le critère de variabilité spatiale qui est pertinent pour distinguer les composantes spatiales les unes des autres. Celui-ci met en évidence les composantes F4, F5 et F6 comme significatives. Les composantes F5 et F6 (Fig.4e et 4f) mettent en évidence les structures Est-Ouest, alors que Fl et F2 (Fig.4a, 4b) mettent en évidence les structures Nord- Sud. La composante F3 (Fig.4c) mélange la structure Nord-Sud avec le bruit aléatoire. La composante F4 (Fig.4d) capture le bruit aléatoire et le sépare tout à fait des autres structures. Les types de graphiques sur les Fig.5a à Se et les Fig.6a à 6e (corrélations et contributions à la variabilité spatiale des variables) permettent de positionner les images initiales par rapport aux composantes spatiales. Ce positionnement joue un rôle d'interprétation des données (initiales et composantes) les unes par rapport aux autres. C'est ainsi que les plans factoriels (F5-F6), aussi bien pour les corrélations (Fig.5a-5e) que pour les contributions (Fig.6a-6e) montrent que les composantes F5 et F6 sont fortement corrélées avec (et contribuent fortement aux) images Zi, Z2 et Z3. De même, on peut voir grâce à la composante spatiale F4, la forte variabilité spatiale des images Z2, Z4 et Z6. La structure composante F5-F6 correspond à la structure spatiale EstOuest. Celle-ci est principalement contenue dans les images ZI, Z2 et Z3 (Fig. 7a-7c). Les images des Fig.7a à 7f correspondent à l'extraction de cette structure, des images initiales. On remarque que celles qui la contiennent (ZI, Z2 et Z3) en restituent une grande part, alors que les autres (Z4, Z5 et Z6) n'en fournissent qu'une quantité pratiquement négligeable (Fig.7d à 7f). La structure F4 correspond au bruit aléatoire. Les images des Fig.8a à 8f visualisent ce bruit contenu dans chacune des images initiales, principalement dans les images Z2, Z4 et Z6 (respectivement Fig.8b, 8d et 8f). Le tableau ci-dessous donne les valeurs des pourcentages de variabilité 2o spatiale (cnt) et de variance (cnV) des structures composantes F5- F6 et F4 contenues dans les images initiales. Ces valeurs permettent de sélectionner les images initiales en fonction des structures composantes mises en évidence. Extraction de la Extraction de la structure composante structure composante F5-F6 F4 % de cnt t% de cntV % de cnt % de cntV ZI 90% 69% 1 % 1% Z2 78 % 62 % 19 % j 28 % Z3 88 % 65 % 1% o 2 %7 Z4 1% [ 0% 75 % 46% Z5 3% 1% i 19% 9% Z6 2 % 1% i 76 %7 53 % On a décrit des exemples d'application de la méthode à la caractérisation d'un volume du sous-sol. Il est bien évident cependant que la méthode peut s'appliquer plus généralement à la caractérisation de toutes données spatiales, o l'on a à traiter d'unités reliées entre elles par des graphes de continuité spatiale. La méthode peut trouver des applications par exemple et de façon non limitative, dans toutes les géosciences, notamment en géographie, agronomie, hydrologie, géomarketing, etc. Dans les exemples qui ont été décrits, on a considéré des objets (éléments d'image) qui ont une répartition spatiale régulière. Il est bien évident cependant que la méthode s'applique à des objets présentant une répartition spatiale quelconque. REVENDICATIONS 1) Méthode statistique multivariable pour analyser des données associées à des éléments d'images traduisant des événements d'un milieu complexe tel que le sous-sol par exemple, ces données étant obtenues par une exploration du milieu, de i manière à faire apparaître les relations spatiales entre ces éléments, caractérisée en ce qu'elle comporte: - une analyse spatiale des données pour mettre en évidence les propriétés spatiales des dits événements dans plusieurs directions; - l'application aux données d'une technique d'analyse de la contiguïté spatiale pour trouver leur meilleure décomposition en structures spatiales composantes; et - le filtrage des structures obtenues par cette décomposition, pour extraire les plus pertinentes. 2) Méthode selon la revendication 1, caractérisée en ce qu'elle comporte: - une identification des structures spatiales des données; - une séparation de ces structures spatiales avec élimination d'éventuelles redondances; - la formation, à partir des images initiales, d'images synthétiques ou composantes spatiales qui mettent en évidence les structures spatiales des données; - l'établissement d'une typologie des images initiales en fonction des structures spatiales mises en évidence par les composantes spatiales; et - le filtrage des images initiales pour enlever le bruit et sélectionner une ou plusieurs structures spatiales identifiées. 3) Méthode selon la revendication 1 ou 2, caractérisée en ce que l'on identifie la structure spatiale des données par une analyse de différents variogrammes des éléments d'image. 4) Méthode selon la revendication 3, caractérisée en ce que l'on procède à une analyse variographique monovariable des éléments d'image. ) Méthode selon la revendication 3, caractérisée en ce que l'on procède à une analyse variographique bivariable des éléments d'image. 6) Méthode selon l'une des revendications précédentes, caractérisée en ce qu'elle comporte la sélection des structures spatiales composantes en déterminant leurs contributions respectives à la variabilité spatiale des données et à la variance des données. 7) Méthode selon la revendication précédente, caractérisée en ce que la dite sélection est effectuée graphiquement. 8) Méthode selon l'une des revendications précédentes, caractérisée en ce qu'elle comporte une quantification des structures spatiales sélectionnées à la fois en termes de variabilité spatiale et de variance. 9) Application de la méthode à l'analyse de données spatiales à répartition régulière ou irrégulière.