FR 3152604 A1 20250307 FR 2309292 A 20230905 L’invention concerne un système d’analyse statistique temporelle d’anomalies spatiales destiné à la détection et à la visualisation de zones de mouvements de fluide dans un sous-sol géologique, ledit système comportant un nombre entier M de capteurs (2), ledit nombre entier M étant supérieur ou égal à 3, lesdits capteurs (2) étant configurés pour mesurer une onde (3) se propageant dans ledit sous-sol géologique sur une période d’observation pour obtenir au moins un signal de mesure et un contrôleur configuré pour mettre en œuvre plusieurs étapes de traitement dudit signal de mesure. L’invention concerne également un procédé permettant la détection et la visualisation de zones de mouvements de fluide. L’invention concerne en outre un programme ordinateur comprenant des instructions de code de programme pour l’exécution des étapes dudit procédé, ledit programme fonctionnant sur un contrôleur. Figure de l’abrégé : Fig. 1 SYSTEME D’ANALYSE STATISTIQUE TEMPORELLE D’ANOMALIES SPATIALES DESTINE A LA DETECTION ET A LA VISUALISATION DE ZONES DE MOUVEMENTS DE FLUIDE ET PROCEDE ASSOCIE L’invention porte sur un système et un procédé mis en œuvre dans le domaine de la géologie et, plus particulièrement, de la prospection sismique. Il est connu de l’art antérieur une demande de brevet WO2020038985 divulguant un procédé de détection d’un fluide dans le sous-sol et un système associé décrivant des moyens mettant en œuvre ledit procédé. Il est proposé un procédé de détection simplifié par rapport à l’art antérieur tout en préservant une précision identique. Le procédé comporte une étape de mesure au cours de laquelle des capteurs d’ondes détectent une onde mécanique, chaque capteur étant apte à reconnaître l’amplitude et la fréquence de l’onde. La détection de l’onde par le capteur engendre un signal qui subit, par la suite, une étape de découpage sur plusieurs intervalles temporels de durée prédéfinie, comprise entre 5 minutes et 2 heures. Ainsi, des échantillons du signal mesuré sont obtenus. Chaque signal correspond à un intervalle temporel de découpage. Une étape supplémentaire de filtrage et/ou d’égalisation spectrale du signal peut être effectuée entre l’étape de mesure et l’étape de découpage. Les deux étapes de filtrage et d’égalisation spectrale permettent de réduire les défauts ou anomalies du bruit sismique dans l’objectif d’améliorer la précision du signal. Une étape optionnelle de normalisation par l’énergie peut être exécutée durant laquelle chaque échantillon du signal est uniformiser dans le but que chaque échantillon présente une énergie semblable les uns par rapport aux autres. Une étape optionnelle d’autocorrélation d’une fenêtre de bruit sismique a lieu. L’autocorrélation a pour but de reconstituer un champ d’onde cohérent mesuré par le capteur. Un diagramme reprenant l’ensemble des autocorrélations mesurées dans le temps est établi grâce à ces mesures. Ainsi, il est possible de distinguer un changement important des valeurs d’autocorrélations, indiquant alors la présence d’un fluide. Optionnellement, une référence est définie. La référence correspond à la moyenne de l’ensemble des valeurs d’autocorrélations des échantillons du signal considéré sur une période temporelle donnée. Une étape de calcul de la cohérence temporelle des échantillons du signal mesuré est ensuite mise en œuvre. La cohérence temporelle correspond à la valeur d’un même échantillon à deux instants différents. Il existe également une étape de détermination de la présence d’un fluide à l’aide du calcul de la cohérence temporelle. Les variations des valeurs de cohérences temporelles indiquent effectivement la présence d’un fluide mais elles indiquent également les changements de propriétés du fluide. Par changement de propriétés, il peut s’agir notamment d’un changement de sa composition ou de sa migration. Cependant, un inconvénient subsiste. S’il est possible de déterminer sur des graphiques la présence ou non de fluides, aucun des moyens cités ne permet de repérer géographiquement les zones où sont situés les fluides. Il est alors difficile de déterminer avec précision où le sous-sol géologique doit être inspecté pour y extraire lesdits fluides. L’objectif de la présente invention est de remédier à ces inconvénients en proposant un système permettant l’établissement d’une carte géographique indiquant les zones les plus susceptibles d’abriter des fluides. Pour atteindre cet objectif, l’invention propose un système d’analyse statistique temporelle d’anomalies spatiales destiné à la détection et à la visualisation de zones de mouvements de fluide dans un sous-sol géologique comprenant : - un nombre entier M de capteurs, ledit nombre entier M étant supérieur ou égal à 3, lesdits capteurs étant configurés pour mesurer au moins une onde se propageant dans ledit sous-sol géologique sur une période d’observation pour obtenir au moins un signal de mesure, lesdits capteurs formant paires de capteurs, dont un premier capteur et un deuxième capteur formant une paire de capteurs ; - un contrôleur configuré pour : - reconstruire une série de formes d’ondes par interférométrie dudit signal de mesure provenant de ladite paire de capteurs à un intervalle de temps régulier sur ladite période d’observation ; - extraire une série temporelle de cohérence associée à ladite paire de capteurs, l’ensemble desdites paires de capteurs formant une pluralité de séries temporelles de cohérence ; - interpoler spatialement ladite pluralité de séries temporelles de cohérence ; - représenter une pluralité de résultats issus de l’interpolation spatiale de ladite pluralité de séries temporelles de cohérence sous la forme d’une carte de cohérence présentant une pluralité de mailles, chaque maille représentant une partie élémentaire du sous-sol géologique ; - paramétrer une durée d’analyse statistique ; - calculer une pluralité de valeurs d’occurrence d’anomalies Oc i k correspondant au nombre de fois où la valeur d’une maille est inférieure à une valeur seuil prédéfinie ; - calculer la fonction de densité de probabilité d’apparition d’une anomalie pour ladite maille par la formule suivante : avec étant la probabilité d’anomalie de cohérence pour la -ème maille à la -ème période d’analyse de N jours ; et étant l’occurrence d’anomalie de cohérence pour la -ème maille à la -ème période d’analyse de N jours ; - établir une carte de probabilité d’occurrence d’anomalie par répertoriage de chaque valeur de probabilité d’apparition d’une anomalie pour la maille ; - obtenir une fréquence à laquelle l’anomalie survient et une localisation des changements dans ledit sous-sol géologique sur la base de la probabilité d’anomalie de cohérence. Grâce à l’invention, il est possible de localiser de manière précise les zones du sous-sol géologique présentant une grande probabilité d’abriter un fluide. Avantageusement, ledit contrôleur est configuré pour interpoler spatialement ladite pluralité de séries temporelles de cohérence par un processus d’inversion tomographique. L’interpolation spatiale par inversion tomographique permet une pondération de l’information de manière irrégulière sur la carte de probabilité d’occurrence d’anomalie. Ainsi, ce processus permet aussi la correction des effets géométriques résultant de l’ensemble des capteurs présents sur le territoire étudié. Avantageusement, ledit contrôleur est configuré pour interpoler spatialement ladite pluralité de séries temporelles de cohérence par un calcul de valeurs moyennes réparties sur des trajectoires de rais. De préférence, ledit nombre entier M est égal à 3, lesdits capteurs étant alignés. Ainsi, le système est fonctionnel à partir de seulement 3 capteurs. Avantageusement, lesdits capteurs sont des géophones ou des hydrophones. Ainsi, la détection des ondes sont permises sur une terre émergée ou sur un fond marin. De préférence, ladite onde est une onde émanant d’une source naturelle. L’étude du bruit sismique ambiant, ou bruit sismique continu, est ainsi permise. De préférence, ladite onde est une onde émanant d’une source artificielle. Avantageusement, ladite carte comporte entre 10 et 20 000 mailles, de préférence entre 100 et 10 000 mailles, plus préférentiellement entre 1000 et 5000 mailles. En outre, l’invention concerne également un procédé d’analyse statistique temporelle d’anomalies spatiales destiné à la détection et à la visualisation de zones de mouvements de fluide dans ledit sous-sol géologique mettant en œuvre ledit système selon la présente invention comportant les étapes suivantes : - une étape d’acquisition d’au moins un signal de mesure pendant une période d’observation ; - une étape de reconstruction d’une série de formes d’ondes par interférométrie dudit signal de mesure de ladite paire de capteurs à intervalles de temps réguliers sur ladite période d’observation ; - une étape d’extraction d’une série temporelle de cohérence associée à ladite paire de capteurs par la mesure de l’évolution d’une forme d’onde par autocorrélation sur une période donnée, l’ensemble desdites paires de capteurs formant une pluralité de séries temporelles de cohérence ; - une étape d’interpolation spatiale de ladite pluralité de séries temporelles de cohérences ; - une étape de représentation des résultats issus de l’interpolation spatiale de la pluralité de séries temporelles de cohérence sous la forme d’une carte de cohérence présentant une pluralité de mailles, chaque maille représentant une partie élémentaire du sous-sol géologique ; - une étape de paramétrage d’une durée d’analyse statistique ;- une étape de calcul des valeurs d’occurrence d’anomalies Oc i k correspond au nombre de fois où la valeur d’une maille est inférieure à une valeur seuil ; - une étape de calcul de la fonction de densité de probabilité d’apparition d’une anomalie pour chaque maille par la formule suivante : avec étant la probabilité d’anomalie de cohérence pour la -ème maille à la -ème période de N jours ; et étant l’occurrence d’anomalie de cohérence pour la -ème maille à la -ème période de N jours ; - une étape de modélisation d’une carte de probabilité d’occurrence d’anomalie par répertoriage de chaque valeur de probabilité d’apparition d’une anomalie pour chaque maille ; - une étape d’obtention de la fréquence à laquelle l’anomalie survient et de la localisation des changements dans ledit sous-sol géologique sur la base de la probabilité d’anomalie de cohérence. L’invention concerne également un procédé d’analyse statistique temporelle d’anomalies spatiales destiné à la détection et à la visualisation de zones de mouvements de fluide dans un sous-sol géologique, ledit procédé mettant en œuvre un nombre entier M de capteurs, ledit nombre entier M étant supérieur ou égal à 3, lesdits capteurs étant configurés pour mesurer au moins une onde se propageant dans ledit sous-sol géologique sur une période d’observation pour obtenir au moins un signal de mesure, lesdits capteurs formant paires de capteurs, dont un premier capteur et un deuxième capteur formant une paire de capteurs, caractérisé en ce que ledit procédé comportant les étapes suivantes : - une étape d’acquisition d’au moins un signal de mesure pendant une période d’observation ; - une étape de reconstruction d’une série de formes d’ondes par interférométrie du signal de mesure de ladite paire de capteurs à un intervalle de temps régulier sur une période d’observation ; - une étape d’extraction d’une série temporelle de cohérence associée à ladite paire de capteurs par la mesure de l’évolution d’une forme d’onde par autocorrélation sur une période donnée, l’ensemble desdites paires de capteurs formant une pluralité de séries temporelles de cohérence ; - une étape d’interpolation spatiale de ladite pluralité de séries temporelles de cohérences ; - une étape de représentation des résultats issus de l’interpolation spatiale de la pluralité de séries temporelles de cohérence sous la forme d’une carte de cohérence présentant une pluralité de mailles, chaque maille représentant une partie élémentaire du sous-sol géologique; - une étape de paramétrage d’une durée d’analyse statistique ; - une étape de calcul des valeurs d’occurrence d’anomalies Oc i k correspond au nombre de fois où la valeur d’une maille est inférieure à une valeur seuil ; - une étape de calcul de la fonction de densité de probabilité d’apparition d’une anomalie pour chaque maille par la formule suivante : avec étant la probabilité d’anomalie de cohérence pour la -ème maille à la -ème période de N jours ; et étant l’occurrence d’anomalie de cohérence pour la -ème maille à la -ème période de N jours ; - une étape de modélisation d’une carte de probabilité d’occurrence d’anomalie par répertoriage de chaque valeur de probabilité d’apparition d’une anomalie pour chaque maille ; - une étape d’obtention d’une fréquence à laquelle l’anomalie survient et d’une localisation des changements dans ledit sous-sol géologique sur la base de la probabilité d’anomalie de cohérence. De préférence, l’étape d’extraction met en œuvre une méthode de suivi des variations de vitesse desdites formes d’ondes. Avantageusement, l’étape d’interpolation spatiale est réalisée par un processus d’inversion tomographique. Avantageusement, l’étape d’interpolation spatiale est réalisée par un calcul de valeurs moyennes réparties sur des trajectoires de rais. En outre, l’invention concerne aussi un programme d’ordinateur comprenant des instructions de code de programme pour l’exécution des étapes du procédé d’analyse statistique temporelle d’anomalies spatiales destiné à la détection et à la visualisation de zones de mouvements de fluide dans le sous-sol géologique, lorsque ledit programme fonctionne sur un contrôleur. L’automatisation de la mise en œuvre du procédé permet de réduire les moyens humains nécessaires à son déroulement. La limitation des moyens humains engendre une réduction du nombre d’erreurs pouvant advenir et une réduction du temps d’exécution du procédé. L’invention sera davantage détaillée par la description d’un mode de réalisation non limitatif, et sur la base des figures annexées illustrant des variantes de l’invention, dans lesquelles : - FIG. 1 illustre schématiquement un système d’acquisition de données sismiques pour l’analyse statistique temporelle d’anomalies spatiales dans le cadre de la détection et de la visualisation de zones de mouvements de fluide ; - FIG. 2 illustre schématiquement un procédé d’analyse statistique temporelle d’anomalies spatiales destiné à la détection et à la visualisation de zones de mouvements de fluide. Il est illustré schématiquement sur la FIG. 1 un système d’acquisition de données sismiques pour l’analyse statistique temporelle d’anomalies spatiales dans le cadre de la détection et de la visualisation de zones de mouvements de fluide dans un sous-sol géologique. La FIG. 1 représente une surface 1 délimitant ledit sous-sol géologique que l’on souhaite étudier. Quatre capteurs sont disposés sur ladite surface 1. Les capteurs sont autrement appelés géophones lorsqu’ils se situent sur une terre émergée ou hydrophones lorsqu’ils se situent sur un fond marin. Le sous-sol géologique est défini comme désignant l’ensemble de la croûte terrestre. L’invention peut donc être mise en œuvre sur une terre émergée ou sur un fond marin. Communément, les travaux ayant lieu sur la terre émergée sont dit « onshore », tandis que les travaux ayant lieu sur un fond marin sont dit « offshore ». Dans une alternative, les capteurs 2 sont répartis sous la surface 1 dans des forages, autrement appelés des puits. Le sous-sol géologique comporte une pluralité d’ondes 3 issues d’un bruit sismique continu. Le bruit sismique continu correspond à l’ensemble des ondes 3 naturellement présentes dans le sous-sol, sans aucune intervention humaine nécessaire. Les capteurs 2 mesurent la fréquence et l’amplitude de l’onde 3 mécanique générée par le fluide présent dans le sous-sol géologique. Les capteurs 2 comprennent communément une masse lourde connectée par au moins un ressort à un cadre. Lorsque la terre entre en mouvement, le cadre bouge. Cependant, la masse, du fait de son inertie, bouge très peu par rapport à un repère géographique. Les capteurs 2 sont donc aptes à mesurer le mouvement relatif entre la masse et le cadre. Ces informations sont ensuite transmises à un contrôleur (non représenté sur la FIG. 1 ) pour être traitées. Il est illustré schématiquement sur la FIG. 2 , sous la forme d’un logigramme, un procédé d’analyse statistique temporelle d’anomalies spatiales destiné à la détection et à la visualisation de zones de mouvements de fluide. Lors d’une étape d’acquisition E1, un signal de mesure du bruit sismique est enregistré sur une période de temps donné, appelée « période d’observation ». Le bruit sismique est défini comme étant l’ensemble des ondes présentent dans le sol. Les ondes sont issues d’un phénomène naturel tel qu’un évènement météorologique par exemple. Elles peuvent également être issues d’une activité humaine : par exemple, une étude du bruit sismique à proximité d’une voie de chemin de fer sera forcément impactée par le passage d’un train créant des vibrations dans le sol. Le bruit sismique peut donc être un bruit sismique continu, c’est-à-dire un bruit sismique naturellement présent dans le sous-sol géologique. Le bruit sismique peut également être un bruit sismique d’origine artificielle, c’est-à-dire un bruit sismique émanant d’une source artificielle élaborée par l’Homme. Par exemple, il peut s’agir d’un véhicule parcourant le terrain étudié venant ainsi créer des ondes artificielles. Les données relatives au bruit sismique sont obtenues par l’intermédiaire de capteurs, M étant un nombre entier supérieur ou égal à 3, et lesdits capteurs étant mis en œuvre par paire. Parmi les M capteurs, il existe paires de capteurs, dont au moins un premier capteur et un deuxième capteur formant une paire de capteurs. Lorsque ledit nombre M est égale à 3, lesdites capteurs sont alignés de façon que les diverses trajectoires reliant chaque capteur entre eux se superposent. Les capteurs sont reliés à un contrôleur. Les données concernant la fréquence et l’amplitude des ondes détectées par ladite paire de capteurs sont transmises électriquement au contrôleur afin qu’elles soient traitées. De préférence, le contrôleur est un ordinateur. Lors d’une étape de reconstruction E2, une série de formes d’ondes sont reconstruites par interférométrie à partir dudit signal de mesure de ladite paire de capteurs à un intervalle de temps réguliers sur ladite période d’observation. L’interférométrie sismique est une technique d’investigation du sous-sol par l’étude des ondes se superposant créant ainsi un phénomène d’interférence. Le phénomène d’interférence se caractérise par la capacité de deux ondes de même fréquence à se combiner. Il existe deux types d’interférences. L’interférence constructive a lieu lorsque deux ondes de même fréquence se combinent en phase. Dans le cas d’une interférence constructive, une onde de plus grande amplitude est créée par la combinaison des deux ondes initiales de moindre amplitude. L’interférence destructive a lieu lorsque deux ondes de même fréquence se combinent de manière déphasée. De cette manière, l’amplitude des deux ondes seront combinés et donneront naissance à une onde de moindre amplitude, les amplitudes de chacune des deux ondes initiales s’annulant l’une et l’autre. Les ondes initialement recueillies lors de l’étape E1 sont donc combinées afin d’être reconstruites selon que l’on soit dans le cas d’une interférence constructive ou d’une interférence destructive. La reconstruction de ces ondes est effectuée à un intervalle de temps régulier durant la période d’observation précédemment définie. Par exemple, l’intervalle de temps régulier est compris entre 10 minutes et 50 minutes, de préférence 30 minutes. Ainsi, toutes les 30 minutes, lesdites ondes sont reconstruites. Lors d’une étape d’extraction E3, une série temporelle de cohérence associée à ladite paire de capteurs est extraite de la série de formes d’ondes établies suite à l’étape de reconstruction E2. Une pluralité de paires de capteurs étant présente, chaque paire de capteurs engendre une série temporelle de cohérence. Il existe donc une pluralité de séries temporelles de cohérence. La cohérence est l’ensemble des propriétés de corrélation temporelle ou spatiale des ondes. Lorsque l’on mesure les propriétés d’une même onde à deux temps différents, la cohérence temporelle est étudiée. Lorsque l’on mesure les propriétés d’une même onde à deux endroits différents, la cohérence spatiale est étudiée. Pour procéder à cette extraction, une mesure de l’évolution d’une forme d’onde est mise en œuvre. La mesure de l’évolution d’une forme d’onde est basée sur la cohérence de formes d’ondes corrélées de bruit sismique. La mesure de l’évolution d’une forme d’onde est effectuée par autocorrélation sur une période donnée. Ainsi, la série de formes d’ondes est analysée afin de suivre son évolution sur une période de temps définie. A partir de ces données sont extraits des séries temporelles de cohérence. Dans une alternative, une méthode de suivi des variations de vitesse desdites formes d’ondes est appliquée. De cette manière, une carte de probabilité d’occurrence d’anomalie de variation de vitesse est obtenue à partir des valeurs de variations de vitesse desdites formes d’ondes. Lors d’une étape d’interpolation spatiale E4, lesdites séries temporelles de cohérence associées à l’ensemble desdites paires de capteurs sont traitées afin de se présenter sous la forme de séries temporelles d’une carte de cohérence. L’interpolation spatiale est définie comme étant un ensemble de fonctions mathématiques permettant d’illustrer des valeurs sous forme de points dont les valeurs sont connues pour estimer des valeurs non connues. De préférence, l’interpolation spatiale est effectuée par un processus d’inversion tomographique. Dans le domaine sismique, la tomographie est définie comme étant une technique d’imagerie très utilisée en géologie. Elle permet notamment l’analyse des structures géologiques. La technique consiste à enregistrer le temps d’arrivée des ondes, naturelles ou créées artificiellement, émanant du sol au cours du temps par l’intermédiaire d’un capteur sismologique, autrement appelé un géophone ou bien une station sismique. Ces données sont transmises électriquement à un sismographe qui génère alors un sismogramme. L’ensemble des sismogrammes générés permet de déterminer la vitesse de propagation des ondes. Les propriétés du sous-sol, tel que la présence d’un fluide, sont déduites à partir de ces informations sur la vitesse de propagation des ondes. Ces données sont ensuite traitées par un processus d’optimisation, aussi nommé problème inverse, on parle alors de processus d’inversion tomographique. Le processus d’inversion tomographique permet l’élaboration d’une image tridimensionnelle du sous-sol révélant ainsi de façon précise la probable présence d’une zone de mouvements de fluide par exemple. Cependant, tout autre méthode d’interpolation spatiale peut être utilisée pour mettre en œuvre l’étape d’interpolation spatiale E4. Dans un autre exemple, l’interpolation spatiale est réalisée par un calcul de valeurs moyennes réparties sur des trajectoires de rais. La trajectoire de rais est ici définie comme étant le chemin parcouru dans l’espace entre chaque paire de capteurs, soit entre au moins ledit premier capteur et ledit deuxième capteur. En un premier point de l’espace sur ladite trajectoire, la moyenne de toutes les valeurs de cohérence passant par ledit point est calculée. Entre ledit premier point et un second point de l’espace, des valeurs aléatoires sont réparties entre ledit premier point et ledit second point. Lors d’une étape de représentation E5, les résultats de l’interpolation spatiale des séries temporelles de cohérence issues de la pluralité de paires de capteurs sont représentés sous forme d’une carte de cohérence. La carte de cohérence présente des mailles, chacune des mailles représentant une partie élémentaire du sous-sol géologique étudié par lesdits capteurs. En effet, ladite carte de cohérence représentant le territoire étudié couvert de capteurs est quadrillé de façon que ladite carte est découpée en X nombre entier de mailles, chaque maille représentant une partie dudit territoire. Chaque maille représentant la même superficie. Plus le nombre X est grand, plus la résolution de ladite carte est améliorée. Dans le cadre de l’invention, ladite carte peut comporter entre 10 et 20 000 mailles pour que la résolution permette d’exploiter les résultats de l’analyse. La carte peut comporter entre 100 et 10 000 mailles, voire même entre 1000 et 5000 mailles. Lors d’une étape de paramétrage E6, une durée d’analyse statistique est déterminée : il s’agit de la période de temps durant laquelle l’analyse statistique se déroulera. Lors d’une étape de calcul des valeurs d’occurrence d’anomalies Oc i k E7, le nombre de fois durant laquelle une maille a présenté une anomalie durant la période d’analyse statistique est comptabilisé. Cette valeur correspond à la valeur d’occurrence d’anomalies Oc i k . Une étape de calcul de la fonction de densité de probabilité d’apparition d’une anomalie E8 est effectuée pour chaque maille. Le calcul est réalisé par la formule suivante : avec étant la probabilité d’anomalie de cohérence pour la -ème maille à la -ème période de N jours ; et étant l’occurrence d’anomalie de cohérence pour la -ème maille à la -ème période de N jours. Ce calcul permet de déterminer la probabilité d’apparition d’une anomalie et donc la probabilité de présence d’un fluide dans le sous-sol. Lors d’une étape de modélisation E9, la carte de probabilité d’occurrence est établie en répertoriant chaque valeur de probabilité d’apparition d’une anomalie pour chaque maille. Ainsi, il est facile de visualiser quelles sont les zones du sous-sol abritant probablement un fluide. Avantageusement, un programme ordinateur permet la mise en œuvre du procédé précédemment exposé. De cette façon, le procédé est mis en œuvre automatiquement sans nécessité une quelconque intervention humaine. Système d’analyse statistique temporelle d’anomalies spatiales destiné à la détection et à la visualisation de zones de mouvements de fluide dans un sous-sol géologique comprenant : - un nombre entier M de capteurs, ledit nombre entier M étant supérieur ou égal à 3, lesdits capteurs étant configurés pour mesurer au moins une onde se propageant dans ledit sous-sol géologique sur une période d’observation pour obtenir au moins un signal de mesure, lesdits capteurs formant paires de capteurs, dont un premier capteur et un deuxième capteur formant une paire de capteurs ; - un contrôleur configuré pour : - reconstruire une série de formes d’ondes par interférométrie dudit signal de mesure provenant de ladite paire de capteurs à un intervalle de temps régulier sur ladite période d’observation ; -extraire une série temporelle de cohérence associée à ladite paire de capteurs, l’ensemble desdites paires de capteurs formant une pluralité de séries temporelles de cohérence ; - interpoler spatialement ladite pluralité de séries temporelles de cohérence ; - représenter une pluralité de résultats issus de l’interpolation spatiale de ladite pluralité de séries temporelles de cohérence sous la forme d’une carte de cohérence présentant une pluralité de mailles, chaque maille représentant une partie élémentaire du sous-sol géologique ; - paramétrer une durée d’analyse statistique ; - calculer une pluralité de valeurs d’occurrence d’anomalies Oc i k correspondant au nombre de fois où la valeur d’une maille est inférieure à une valeur seuil prédéfinie ; - calculer la fonction de densité de probabilité d’apparition d’une anomalie pour ladite maille par la formule suivante : avec étant la probabilité d’anomalie de cohérence pour la -ème maille à la -ème période d’analyse de N jours ; et étant l’occurrence d’anomalie de cohérence pour la -ème maille à la -ème période d’analyse de N jours ; - établir une carte de probabilité d’occurrence d’anomalie par répertoriage de chaque valeur de probabilité d’apparition d’une anomalie pour la maille ; - obtenir une fréquence à laquelle l’anomalie survient et une localisation des changements dans ledit sous-sol géologique sur la base de la probabilité d’anomalie de cohérence. Système selon la revendication 1 caractérisé en ce que ledit contrôleur est configuré pour interpoler spatialement ladite pluralité de séries temporelles de cohérence par un processus d’inversion tomographique. Système selon la revendication 1 caractérisé en ce que ledit contrôleur est configuré pour interpoler spatialement ladite pluralité de séries temporelles de cohérence par un calcul de valeurs moyennes réparties sur des trajectoires de rais. Système selon l’une quelconque des revendications 1 à 3 caractérisé en ce que ledit nombre entier M est égal à 3, lesdits capteurs étant alignés. Système selon l’une quelconque des revendications 1 à 4 caractérisé en ce que lesdits capteurs sont des géophones ou des hydrophones. Système selon l’une quelconque des revendications 1 à 5 caractérisé en ce que ladite onde est une onde émanant d’une source naturelle. Système selon l’une quelconque des revendications 1 à 5 caractérisé en ce que ladite onde est une onde émanant d’une source artificielle. Système selon l’une quelconque des revendications 1 à 7 caractérisé en ce que ladite carte comporte entre 10 et 20 000 mailles, de préférence entre 100 et 10 000 mailles, plus préférentiellement entre 1000 et 5000 mailles. Procédé d’analyse statistique temporelle d’anomalies spatiales destiné à la détection et à la visualisation de zones de mouvements de fluide dans un sous-sol géologique, ledit procédé mettant en œuvre un nombre entier M de capteurs, ledit nombre entier M étant supérieur ou égal à 3, lesdits capteurs étant configurés pour mesurer au moins une onde se propageant dans ledit sous-sol géologique sur une période d’observation pour obtenir au moins un signal de mesure, lesdits capteurs formant paires de capteurs, dont un premier capteur et un deuxième capteur formant une paire de capteurs, caractérisé en ce que ledit procédé comportant les étapes suivantes : - une étape d’acquisition (E1) d’au moins un signal de mesure pendant une période d’observation ; - une étape de reconstruction (E2) d’une série de formes d’ondes par interférométrie du signal de mesure de ladite paire de capteurs à un intervalle de temps régulier sur une période d’observation ; - une étape d’extraction (E3) d’une série temporelle de cohérence associée à ladite paire de capteurs par la mesure de l’évolution d’une forme d’onde par autocorrélation sur une période donnée, l’ensemble desdites paires de capteurs formant une pluralité de séries temporelles de cohérence ; - une étape d’interpolation spatiale (E4) de ladite pluralité de séries temporelles de cohérences ; - une étape de représentation (E5) des résultats issus de l’interpolation spatiale de la pluralité de séries temporelles de cohérence sous la forme d’une carte de cohérence présentant une pluralité de mailles, chaque maille représentant une partie élémentaire du sous-sol géologique; - une étape de paramétrage (E6) d’une durée d’analyse statistique ; - une étape de calcul des valeurs d’occurrence d’anomalies Oc i k (E7) correspond au nombre de fois où la valeur d’une maille est inférieure à une valeur seuil ; - une étape de calcul de la fonction de densité de probabilité d’apparition d’une anomalie (E8) pour chaque maille par la formule suivante : avec étant la probabilité d’anomalie de cohérence pour la -ème maille à la -ème période de N jours ; et étant l’occurrence d’anomalie de cohérence pour la -ème maille à la -ème période de N jours ; - une étape de modélisation (E9) d’une carte de probabilité d’occurrence d’anomalie par répertoriage de chaque valeur de probabilité d’apparition d’une anomalie pour chaque maille ; - une étape d’obtention (E10) d’une fréquence à laquelle l’anomalie survient et d’une localisation des changements dans ledit sous-sol géologique sur la base de la probabilité d’anomalie de cohérence. Procédé selon la revendication 9 caractérisé en ce que l’étape d’extraction (E3) met en œuvre une méthode de suivi des variations de vitesse desdites formes d’ondes. Procédé selon l’une quelconque des revendications 9 à 10 caractérisé en ce que l’étape d’interpolation spatiale (E4) est réalisée par un processus d’inversion tomographique. Procédé selon l’une quelconque des revendications 9 à 10 caractérisé en ce que l’étape d’interpolation spatiale (E4) est réalisée par un calcul de valeurs moyennes réparties sur des trajectoires de rais. Programme d’ordinateur comprenant des instructions de code de programme pour l’exécution des étapes du procédé d’analyse statistique temporelle d’anomalies spatiales destiné à la détection et à la visualisation de zones de mouvements de fluide dans le sous-sol géologique selon l’une quelconque des revendication 9 à 12, lorsque ledit programme fonctionne sur un contrôleur.