FR 2836572 A1 20030829 FR 0202578 A 20020227 de salle induit par la zone d'écoute (4). - 1 - 2836572 La présente invention concerne un procédé pour la détermination d'un indice de qualité spatiale de données régionalisées. Elle s'applique plus particulièrement, mais non exclusivement, aux données géophysiques, aux données d' imageries obtenues par des méthodes physiques telles qu'imagerie médicale, sonar, contrôle non destructif de matériaux ou encore à tout type d'échantillonnage de phénomènes naturels tels que campagnes de reconnaissance minière, inventaires géochimiques, capteurs de pollutions, données satellite, océanographiques, analyse d'eau... Les données régionalisées sont des données repérées par des coordonnées dans un esp ace à N dimensions, le plus couramment dans un esp ace géographique à une, deux ou trois dimensions. Ces données peuvent étre mono ou multivariables, c'est-à-dire qu'une ou plusieurs variables sont mesurées ou calculées aux points de données. La théorie des variables régionalisées est connue sous le nom de génstatistique. La géostatistique est l' application des probabilités aux phénomènes naturels qui se développent dans l'espace et dans le temps, d'o le préfixe géo. Cette théorie se trouve notamment exposoe dans l'ouvrage " La théorie des variables régionalisées " de G. MATHERON (MASSON Editeur). - 2 - 2836572 Cette théorie fournit le langage et les outils appropriés à l' estimation de toute quantité, inconnue a priori, repérable dans un espace donné, à partir d'un échantillonnage forcément fragmentaire de cette même quantité. Pour estimer cette quantité inconnue, la géostatistique propose de choisir judicieusement le modèle probabiliste le plus approprié à la situation, 1 jestimateur géostatistique étant connu sous le nom de krigeage. Plus que l' estimation, le modèle probabiliste donne également un indicateur de la précision de l' estimation. Cet indicateur, appelé variance d' estimation, est un outil précieux car il ouvre la voie à une possible gestion des incertitudes (traduites en termes de variance). Dans le cadre des modèles probabilistes stationnaires, qui suppose l'invariance par translation dans l'espace de la moyenne de la variable modélisée, l'outil covariance ou variogramme est utilisé pour quantifier la variabilité spatiale des 1 5 données. Pour un modèle non-stationnaire, c'est la covariance généralisée qui est employée. Les modèles géostatistiques permettent également d'anticiper valablement sur un état futur, par exemple de l' exploitation de la ressource naturelle, lorsque les données disponibles seront plus nombreuses et les problèmes d'estimation se poseront différemment à l'opérateur. Quel que soit le contexte métier de l' exploitation de la ressource naturelle, la question se pose touj ours de l'adéquation des données disponibles à la résolution du problème opérationnel. A la qualité intrinsèque de chaque donnée s'ajoute la qualité de l'intogration spatiale de cette donnée au sein du jeu entier. C'est la raison pour laquelle il est intéressant de compléter le relevé expérimental par un contrôle géostatistique associé à des coordonnées géographiques, temporelles ou autres. - 3 - 2836572 Les procédés habituels de contrôle de qualité ou de cohérence des jeux de données régionalisces sont soit visuels, soit morphologiques (études des formes) soit statistiques (sans prendre en compte les coordonnées spatiales). Lorsqu'ils sont utilisés, les procédés de filtrage (fréquentiels ou spatiaux) travaillent généralement sur des données monovariables et sur des grilles régulières. En conséquence, ils sont mal adaptés à la décomposition de donnces multivariables irrégulièrement situées dans l'espace en composantes anomalique et cohérente. De même, la définition des critères utilisés pour définir les anomalies est souvent arbitraire et se prête mal à la vérification expérimentale. En vue de supprimer ces inconvénients, l'invention propose de quantifier la qualité spatiale d'un jeu de données régionalisées au travers de la détermination d'un indice géostatistique appelé SQI pour " Spatial Quality Index " (Indice de Qualité Spatiale), cet indice étant utilisé pour localiser des données a priori anomaliques et ainsi juger de la qualité des mesures ou du traitement numérique qui ont généré le jeu de données. La détermination du SQI résout à la fois le problème de l'interprétation des variations spatiales des données mono ou multivariables en termes généraux d' anomalies et de composante cohérente et l' estimation du degré d' anomalie ou d'incohérence spatiale présent dans chaque donnce prise individuellement. La détermination du SQI ne suppose aucune disposition particulière des données dans l'espace et travaille aussi bien des donnces irrégulièrement réparties dans l'espace que régulièrement situces aux n_uds d'une grille à N dimensions, par exemple une grille d'acquisition à trois dimensions de données géophysiques réparties irrégulièrement est définie selon deux axes longitudinal-transversal d' acquisition des données et un troisième axe vertical temps. - 4 - 2836572 Avantageusement, la détermination de cet indice se fait grâce au procédé selon l'invention qui comprend les phases opératoires suivantes: - une première phase d' identification des anomalies statistiques d'un premier ordre à partir du jeu de données régionalisées brutes, cette identification comprenant une stationnarisation des données par extraction préliminaire des dérives spatiales desdites donnces et la détermination des résidus stationnaires associés, dits de premier ordre" de manière à ce que la valeur de la moyenne de données résiduelles soit raisonnablement constante dans l'espace, les anomalies étant identifiées et examinées sur les résidus de premier ordre de manière à fournir un premier critère d'anomalie du premier ordre, - une deuxième phase d'identification d'anomalies statistiques d'un second ordre avec extraction des composantes des résidus de premier ordre considérées comme anomaliques et des composantes des résidus de premier ordre considérées comme cohérentes dans l'espace, - l'établissement d'une relation quantifiée (SQI) entre toute combinaison des valeurs estimées des composantes anomaliques du premier et/ou du second ordre et toute combinaison des valeurs estimées des composantes cohérentes du premier et/ou du second ordre, - la localisation des anomalies spatiales à partir des valeurs du SQI de .... chaque donnee regonalsee. Avantageusement, ledit établissement d'une relation quantifiée est la détermination, pour chaque donnée régionalisée prise individuellement, du rapport d'un indice de qualité spatiale (SQI). Bien entendu, les susdites phases d'identification pourront s'effectuer par estimation géostatistique (krigeage), en modèle non-stationnaire pour la première phase et en modèle stationnaire pour la deuxième phase. - 5 - 2836572 Dans la première phase, l' estimation non-stationnaire de la dérive spatiale permet l'obtention de résidus stationnaires du premier ordre sur lesquels il est possible de calculer et modéliser valablement un variogramme. L'interprétation de ce variogramme en terme de composantes cohérente et anomalique conduit à l' estimation par modèle stationnaire de la composante anomalique du second ordre. Plus particulièrement, les modèles géostatistiques stationnaire et non stationnaire pourront utiliser: l' estimation par krigeage factoriel des composantes anomalique et cohérente des résidus, - la définition du voisinage de krigeage adaptée à l'estimation de chacune desdites composantes anomalique et cohérente. Bien entendu, à chacune des susdites étapes, l'analyse pourra être facilitée par un contrôle visuel 3D réalisé, d'une part, au moyen d'une interpolation sur un fichier " grille " de n'importe quelle variable irrégulièrement échantillonnce provenant d'un fichier "point" et, d'autre part, au moyen d'un code de couleur associé à la valeur de la variable interpolée. Un mode d'exécution de l 'invention sera décrit ci-après, à titre d' exemple non limitatif, avec réLérence aux dessins annexés dans lesquels: La figure 1 est une représentation schématique illustrant les phases principales de détermination de l'indice de qualité spatiale (coefficient d' anomalies spatiales) dans un contexte de traitement de données géophysiques; - 6 - 2836572 Les figures 2a, 2b et 2c sont des représentations des différents formats de fichiers utilisés dans le procédé selon l'invention; Les figures 3a et 3b sont respectivement la représentation d'un S histogramme fréquentiel (figure 3a) et d'un diagramme temps-vitesse (figure 3b), La figure 4 est une représentation des affichages permettant un contr81e visuel des étapes du procédé; La figure 5 est une représentation de l'interpolation effectuée pour permettre les affichages de la figure 4; La figure 6 est une représentation d'un variogramme expérimental de 1 5 résidus; Les figures 7a et 7b sont des représentations des modes de localisation et de quantification d' anomalies spatiales. Cet exemple concerne plus particulièrement l'acquisition et le traitement des données sismiques 3D pour la caractérisation des réservoirs pétroliers et plus particulièrement le contr81e qualité des mesures de vitesses géophysiques ou vitesses de " stack ". Le problème est le suivant: le contracteur sismique propose à l'opérateur pétrolier un jeu de vitesses pointées manuellement pour l'opération de " stack " qui conditionne la qualité des données finales. L'opérateur pétrolier est chargé de la supervision du travail du contracteur et doit donner son avis sur la qualité du pointé de vitesse. - 7 - 2836572 A cet effet, il peut ausculter les cubes de pointés de vitesses de " stack " à l' aide d' outils statistiques et géostatistiques en vue d'identifier les incohérences spatiales dues à des pointés erronés. La détermination du SQI, conformément au procédé selon l'invention, contribue à définir les pointés jugés anomaliques qu'il convient de reprendre pour garantir une homogéncité spatiale qui sera elle-même quantifiée par la valeur de l'indice SQI. Chaque pointé de vitesse de " stack " est défini par ses coordonnces spatiales - géographiques et temporelle - et une valeur de vitesse géophysique. Dans le domaine spatial décrit par tous les pointés, une modélisation probabiliste permet de séparer un bruit spatial, organisé ou non, d'un signal (spatialement) cohérent. Le bruit spatial est censé correspondre à des arteLacts d' acquisition et de traitement. Quantifié, il permet l'identification rapide des pointés problématiques selon un seuil de tolérance déduit d'impératifs de qualité géophysique ou autre (par exemple, dans le cas du " stack ", préservation des amplitudes du signal sismique). Dans le processus de détermination de l'indice SQI illustré figure 1, le système de coordonnées retenu est le système sismique axe longitudinal (inline) - axe transversal (crossline) -temps (Figure 2a). En effet, toute la châîne géophysique, de l'acquisition au traitement sismique, privilégie ces trois directions principales: la verticale (le temps) et les directions horizontales (inline, crossline) définies par le dispositif d' acquisition. Par conséquent, la plupart des arteLacts d' acquisition et de traitement sont générés le long de ces directions et c'est le long de ces directions qu'est orientée l'analyse. C:e processus permet donc de réduire le temps de détermination. D'une façon plus précise, le processus fait intervenir trois types de fichiers: - 8 - 2836572 - un fichier " point " 3D (cube de pointés de vitesse de " stack ") (Figure 2a), défini par l'ensemble des positions de pointés de vitesses dans un repère de coordonnces axe longitudinal (inline) - axe transversal (crossline) - temps, il est généralement à pas régulier dans le plan horizontal (axe longitudinal (inline) - axe transversal (crossline)) et irrégulier en temps, - un fichier " grille " 2D (Figure 2b): à pas régulier et défini suivant les axes longitudinal (inline) - transversal (crossline), il sert à réaliser des statistiques directionnelles le long de la verticale (axe temps du fichier point 3D de la figure 2a); des statistiques suivant la direction longitudinale ou la direction transversale peuvent être envisagées, - un fichier " grille " 3D (Figure 2c): à maille régulière dans un repère de coordonnées (axe longitudinal (inline) - axe transversal (crossline) temps), il sert pour divers affichages comme il sera expliqué par la suite (Figure 4). Conformément à la méthodologie illustrée figure 1, le cube de pointés de vitesses de " stack" une fois chargé est soumis à un contrôle de qualité géostatistique (calcul de SQI). Les vitesses brutes 1 sont décomposées au premier ordre par krigeage factoriel en une dérive spatiale 2 (composante " basse fréquence ") et des résidus stationnaires 3 (phase QC1). La cohérence spatiale des résidus de premier ordre est modalisée (à l'aide d'un variogramme) (phase QC2) pour réaliser un filtrage discriminant par krigeage factoriel entre un bruit spatial 4, résidu du second ordre, et une partie cohérente du second ordre 5 censée être " nettoyée " des artefacts d' acquisition et de traitement (phase QC3). Cette partie résiduelle cohérente de second ordre est additionnce à la dérive 2, partie cohérente de premier ordre, de manière à engendrer un cube de pointés de vitesses spatialement cobérents 6 (QC4). Dans cet exemple, le SQI est constitué par le rapport entre le bruit, résidu du second ordre, et la partie cohérente de la donnée, somme des 9 2836572 composantes cohérentes du premier et/ou second ordre et conduit à l'obtention d'un cube d'anomalie spatiale 7 (phase QC5): chaque pointé est ainsi caractérisé par son indice de qualité spatiale (SQI) qui exprime le pourcentage de bruit spatial par rapport à la partie spatialement cohérente. Les statistiques expérimentales calculées au cours du processus se répartissent en: - statistiques de base, - statistiques directionnelles, - variogramme expérimental. L'ensemble des pointés de vitesses de " stack " et des résidus de premier ordre forment des distributions qui peuvent étre rapidement analysées grâce à des outils prenant en compte divers paramètres tels que son nombre d'échantillons, ses extrema, sa moyenne arithmétique, son écart type, sa varlance: - son nombre d'échantillons N qui caractérise une distribution Vi (vitesses brutes, résidus, dérives, anomalies, résidus fltrés, vitesses filtrées...), ses extrema: minimum = min (Vi) maximum = max (Vi) - sa moyenne arithmétique: l N m =-ÉV N i= son écart type: (T = (Vi - m)2 - 10 - 2836572 - sa variance: l N 6 =- (Vi - m) N i=l Cette analyse peut étre complétée au moyen: - d'un histogramme fréquentiel dans lequel les valeurs Vi de la variable sont S regroupées en classes, l'histogramme fréquentiel étant la représentation des fréquences correspondant à ces classes (Figure 3a), - d'un diagramme (" cross-plot ") temps-vitesse (Figure 3b). I1 s'avère que les artefacts d'acquisition et de traitement ont comme axes principaux les axes longitudinal, transversal et temps. Le calcul d'une grandeur statistique le long d'une de ces trois directions peut conduire à l'identification d'un ou de plusieurs arteLacts ainsi renforcés. En conséquence, les trois grandeurs statistiques calculées selon ces directions pourront se limiter au nombre d'échantillons, à la moyenne arithmétique et à la variance (oul'écart type). Aux outils d'analyse précédemment mentionnés se rajoute le variogramme expérimental des résidus de premier ordre. Le variogramme permet de quantifier la corrélation spatiale d'une variable régionalisée V(r), r étant le vecteur position défini dans le système de coordonnces longitudinal transversal -temps. Sa formule se déduit de celle du variogramme théorique qui concerne la fonction aléatoire v(0 pour laquelle on ne dispose que d'une seule réalisation: la variable régionalisée. Variogramme théorique: 7(h) = 2 Var[V(r + h) - V()] o h est le vecteur caractérisant un ensemble de couples de pointés " 2836572 Variogramme expérimental (après hypothèses de stationnarité et d'ergodicité...): É(h) = _ [V(rn + h)-V(rn)] o N(h) est le nombre de couples de pointés séparés de h. Par ailleurs, l'analyse pourra être facilitée grâce au mode d'affichage illustré figure 4. Selon ce mode d'affichage, l'interpolation sur grille (à pas régulier) de n'importe quelle variable irrégulièrement échantillonnée (fichier "point") permet, grâce à un code de couleur associé à la valeur de la variable interpolée, un contr81e visuel 3D rapide. L' interpolation retenue pour l' affichage est définie de la manière suivante (Figure 5): à un n_ud de grille Pj, la valeur estimée correspond à l'interpolation linéaire des deux pointés P, P2 respectivement définis par les coordonnées (t, V) et (t2, V2) les plus proches situés de part et d'autre du n_ud Pj et sur la même verticale que celui-ci. Vj =Vi +t i X(V2 -Vi) D ' autres types d' interpolation plus élaborés, comme le krigeage par exemple, pourraient être envisagés pour l'affichage, si le besoin s'en faisait sentir. Les étapes de filtrage et de détermination des anomalies utilisoes dans la méthodologie illustrée figure 1 seront décrites ci-après: - 12 - 2836572 a) Etape de filtrage dérive/résidus Un champ de vites s es prés ente général ement une dérive vertic ale due à la compaction, la compaction étant l'augmentation de la vitesse en fonction de l ' augmentation de la profondeur d' enfouissement (une dérive horizontale peut également exister, si on considère le fond de la mer par exemple pour l'acquisition de données sismiques marines). On observe une non-stationnarité des vitesses dans cette direction qui peut être gérce par la théorie des covariances généralisées, modèle géostatistique non-stationnaire. Mais une covariance généralisée ne peut pas être interprétée " à la main " en terme de composantes spatiales cohérentes et anomaliques, rendant ipso facto impossible l'ajustement d'un modèle en conséquence. I1 faut donc extraire cette dérive et travailler sur les résidus stationnaires associés. L'extraction d'une dérive assurant la stationnarité des résidus est réalisée par l'ajustement polynomial par moindres carrés qui est un cas particulier de krigeage factoriel: la valeur de la dérive en un point de l'espace correspond à la valeur d'un polynôme ajustant au mieux (moindres carrés) les points (indépendants) appartenant à un voisinage centré autour du point à estimer. Le type du polynôme - 1 z, 1 z Z2, 1 X Z X2 Z2 XZ, etc - est à déterminer suivant le type de dérive que l'on désire extraire, z étant le temps et x étant une coordonnée géographique. Les dimensions du voisinage d' extraction de krigeage doivent garantir la stationnarité des résidus de premier ordre. Exemple: extraction d'une dérive de type 1 z Z2 au point rO Vdérive(ro)=a+ bXZo +cxzo2 avec rO(XO,yO,zo) Les coefficients du polynôme (a, b et c) sont obtenus en minimisant le système: Nv 2 É[(a+bxzi +cxZj2)-v(ri)] i=, - 13 - 2836572 o Nv est le nombre d'échantillons contenus dans le voisinage centré autour du point à estimer. b) Filtrage anomalies/résidus cohérents La dérive, composante cohérente de premier ordre, précédemment estimoe est soustraite des vitesses de "stack" brutes. Les résidus de premier ordre doivent être stationnaires. Vrésidu (rO) V(rO) Vdérive (rO) L'analyse variographique des résidus de premier ordre est la phase cruciale du contrôle de qualité spatiale. Le variogramme expérimental des résidus r(h) est calculé dans plusieurs directions. Pour des raisons de clarté des figures, seules les composantes dans les trois directions principales (longitudinale, transversale, temps) sont représentées sur la figure 6. La modélisation du variogramme est subordonnée à une interprétation du variogramme expérimental en termes de composantes spatiales cohérentes et anomaliques. Les compétences combinées du géostatisticien, du géophysicien et éventuellement du géologue sont requises pour cette phase d'interprétation. Le géologue apporte des informations sur les structures géologiques connues ou supposées; le géophysicien est là pour préciser la nature des artetacts géophysiques majeurs susceptibles d'affecter les données; le géostatisticien construit le modèle de variogramrne en tenant compte de ces deux types d'informations. Idéalement, on souhaite séparer un bruit spatial (composante spatiale anomalique) d'un signal cohérent (composante spatiale cohérente) sur la seule - 14 - 2836572 base d'une interprétation variographique. La modélisation proposce dépend de l ' application que l' on souhaite faire du j eu de vitesses (" stack ", conversion profondeur, vitesses de DIX...). Une composante spatiale peut être considérée comme un bruit pour une certaine application ou comme un signal cohérent pour une autre. " Cohérent " et " anomalique " ne sont pas des propriétés intrinsèques du jeu de vitesses, ce sont des propriétés des composantes du jeu de vitesses dans le cadre du modèle géostatistique reconnu. L'ajustement du modèle rM(h) ne peut se faire que par des fonctions de type négatif conditionnel. Dans un premier temps, les modèles de variogrammes comme le modèle de l'effet de pépite, ainsi que les modèles sphérique et exponentiel sont suffisants pour construire un modèle gigogne (à plusieurs composantes). La définition des modèles est donnée dans le cas isotrope: modèle de l'effet de pépite: r r(h) = b silhl>0 modèle sphérique: O silhl>a r(h)= bX3 Xlhl- 1 X:lhl)] siO - 15 - 2836572 Les paramètres a et b s'appellent respectivement la portée et le palier du variogramme. Ils sont tous les deux positifs. Le modèle de variogramme retenu ÉM est une combinaison linéaire de différentes composantes élémentaires classées selon leur interprétation anomalique ou cohérente: (h)=rA(h)+ÉC(h) avec A, composante du variogramme associée à la partie anomalique, I-c' composante du variogramme associée à la partie cohérente. Voisinage Le voisinage regroupe l ' ensemble des points participant à l ' estimation de la composante anomalique située au point rO. Un voisinage dit glissant est indispensable pour toute opération de filtrage. Un voisinage unique englobant l' ensemble des échantillons du champ est très pénalisant en temps de calcul. Les dimensions du voisinage glissant doivent donc optimiser le temps de calcul sans dégrader la qualité de l'estimation. Krigenge factoriel La modélisation du variogramme du second ordre des résidus correspond à une interprétation en terme de composantes cohérentes et anomaliques. Le krigeage factoriel permet une estimation de chacune des deux composantes. L'estimation de la composante anomalique au point rO s'effectue en calculant: N Vrésidu anomalique (rO) X Vrésidu (rct) _=1 o l'ensemble des r, constitue le voisinage glissant de kigeage, - 16 - 2836572 o les poids de krigeage sont déterminés en résolvant le système: Nv x EM (r" - r) + pA = rA (r - rO) pour a = 1,, Nv Nv z\p =1 p=1 L' estimation de la composante cohérente au point rO peut étre obtenue de la même façon par krigeage factoriel. Pour trouver les poids de krigeage correspondant, il suffit de changer dans le système de krigeage rA par rC. Cependant, un seul filtrage est nécessaire puisque par propriété du krigeage factoriel, on peut écrire: Vrésidu cohérent (rO) Vrésidu (rO) Vrésidu anomalique (rO) c) Quantification des anomalies spatiales Calcul du coefficient d'anomalie spatiale La partie spatialement cohérente des résidus de second ordre ajoutée à la dérive, composante cohérente de premier ordre, permet de générer un champ de vitesses spatialement cohérent: Vcobérent (0) Vdérive (0) + Vrésidu cohérent (0) Le rapport entre l'estimation de la composante anomalique et l'estimation de la composante cohérente de la donnée (exprimé en %) constitue un rapport appelé coefficient d'anomalie spatiale ou indice de qualité spatiale SQI. I1 est attaché à chaque pointé de vitesses. SQI attaché aux pointés de vitesses de stack: V (7') Vrésidu anomalique (70) X 1 00 coefficlentd'anomaliespatiale O V - cobérent O J - 17 - 2836572 d) Lo calisation des pointés anomaliques par interprétation du co efficient d'anomalie spatiale La localisation des anomalies spatiales se fait à partir du coefficient d'anomalie spatiale attaché à chaque pointé du cube de vitesses. Deux options sont possibles: - Dans le premier cas, des classes de couleur (ou de symboles) sont associces au coefficient d'anomalie spatiale SQI. Des sections dans le cube sont visualisées. Il est également possible d'interpoler sur grille le coefficient d'anomalie spatiale (Figure 7a). - Dans le second cas, le codage couleur (ou symbolique) se rapporte à la définition d'intervalles temps dans lesquels on recherche des coefficients d'anomalie spatiale supérieurs à une valeur seuil (Figure 7b). - 1 8 - 2836572 Revendications 1. Procédé pour la détermination d'un indice de qualité spatiale (SQI) de données régionalisées destiné à déterminer des données à priori anomaliques et ainsi juger de la qualité des mesures ou du traitement numérique qui ont généré les donnces, caractérisé en ce qu'il comprend les phases opératoires suivantes: - une première phase d' identification des anomalies statistiques d' un premier ordre à partir du jeu de données régionalisoes brutes (1), cette identification comprenant une stationnarisation des données par extraction préliminaire des dérives spatiales (2) desdites données et la détermination des résidus stationnaires associés (3), dits de premier ordre, de manière à ce que la valeur de la moyenne de données résiduelles soit raisonnablement constante dans l'espace, les anomalies étant identifiées et examinées sur les résidus de premier ordre de manière à fournir un premier critère d'anomalie du premier ordre, - une deuxième phase d'identification d'anomalies statistiques d'un second ordre avec extraction des composantes des résidus de premier ordre (3) considérces comme anomaliques (4) et des composantes des résidus de premier ordre (3) considérées comme cohérentes (5) dans l'espace, - l'établissement d'une relation quantifiée (SQI) entre toute combinaison des valeurs estimées des composantes anomaliques du premier et/ou du second ordre (4) et toute combinaison des valeurs estimées des composantes cohérentes du premier et/ou du second ordre (2, 5). 2. Procédé selon la revendication 1, caractérisé en ce que ledit établissement de la relation quantifiée est réalisé pour chaque donnée régionalisée prise individuellement. - 19 - 2836572 3. Procédé selon la revendication 1, caractérisé en ce que ladite relation quantifiée est un rapport. 4. Procédé selon l'une des revendications précédentes, caractérisé en ce que ladite relation quantifiée est le rapport entre la valeur estimoe de la composante anomalique du second ordre et la valeur estimée de la composante cohérente soit la somme des composantes cohérentes de premier et second ordre. 5. Procédé selon la revendication 1, caractérisé en ce qu'il comprend la localisation des anomalies spatiales à partir des valeurs du SQI de chaque donnée régionalisée. 6. Procédé selon la revendication 1, caractérisé en ce que les susdites phases d'identification s'effectuent par estimation géostatistique (krigeage), en modèle non-stationnaire pour la première phase et en modèle stationnaire pour la deuxième phase. 7. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'obtention de résidus stationnaires du premier ordre sur lesquels il est possible de calculer et modéliser un variogramme est réalisée au moyen de l' estimation non-stationnaire de la dérive spatiale. 8. Procédé selon la revendication 7, caractérisé en ce que l'interprétation dudit variogramme en terme de composantes cohérente et anomalique conduit à l'estimation par modèle stationnaire de la composante anomalique du second ordre. 9. Procédé selon la revendication 1, caractérisé en ce que les susdites phases d'identification utilisent: - 20 - 2836572 la définition du voisinage de krigeage adapté à l'estimation de chacune desdites composantes anomalique et cohérente, l'estimation par krigeage factoriel des composantes anomalique et cohérente des résidus. 10. Procédé selon l'une des revendications 1 et 5, caractérisé en ce que, à chacune des susdites étapes, un contrôle visuel 3D de l'analyse est réalisé, d'une part, au moyen d'une interpolation sur un fichier " grille " de n'importe quelle variable irrégulièrement échantillonnée provenant d'un fichier "point" et, d'autre part, au moyen d'un code de couleur associé à la valeur de la variable interpolée. 1 1. Procédé selon la revendication 1, caractérisé en ce que l'extraction d'une dérive assurant la stationnarité des résidus est réalisée par ajustement polynomial par moindres carrés.