FR 3145044 A1 20240719 FR 2300361 A 20230113 La présente invention vise un procédé de surveillance d’une zone d’intérêt d’un sous-sol comprenant un point d’intérêt, comprenant : une extraction d’un premier sous-ensemble (1) de données géophysiques brutes associées au point d’intérêt, parmi un ensemble (0) de données géophysiques brutes acquises antérieurement pour générer une première image (Im1) du sous-sol ; une construction d’un modèle de projection (MOD1) à partir du premier sous-ensemble de données (1) ; une réduction du premier sous-ensemble (1) de données à un deuxième sous-ensemble (2) ; la création d’un ensemble de référence (B) par application du modèle de projection (MOD1) ; une acquisition formant un troisième sous-ensemble (3) ; la création d’un ensemble surveillance (M) par application du modèle de projection (MOD1) ; une comparaison du premier ensemble de référence (B) et du premier ensemble de surveillance (M) pour identifier une évolution du sous-sol à surveiller. Figure de l’abrégé : Figure 6 Procédé de surveillance d’un point d’intérêt d’un sous-sol La présente invention concerne le domaine de l’acquisition et le traitement de données géophysiques. Plus particulièrement, la demande concerne les procédés de surveillance d’une zone d’intérêt d’un sous-sol à partir de données géophysiques ETAT DE LA TECHNIQUE De l’art antérieur, il est connu l’acquisition de données géophysiques permettant de générer une image de la structure géophysique d’un sous-sol terrestre ou marin. Cette image permet à une personne du métier de déterminer avec un haut degré de certitude la présence ou l’absence de la zone d’intérêt, par exemple la présence ou l’absence de ressources naturelles. Des dispositifs configurés pour générer des images d’un sous-sol consistent généralement en une pluralité de sources capables de produire des ondes sismiques sur différentes plages de fréquences. Ces ondes sismiques se diffusent dans un sous-sol et une zone d’intérêt, qui réfléchit une partie desdites ondes sismiques vers un récepteur. En étudiant les variations de propagation des ondes réfléchies dans le sous-sol, il est ainsi possible de produire une image du sous-sol étudié. Ces techniques usuelles de l’état de la technique peuvent être utilisées à intervalles récurrents dans le temps, pour générer une succession d’images du sous-sol, permettant de déterminer l’évolution du profil géophysique de ce sous-sol et réaliser une surveillance. Cette image est dite 4D, c’est-à-dire une image d’un volume du sous-sol évoluant au cours du temps. Actuellement, les solutions d’étude 4D ont pour objet la mise à jour complète d’une acquisition 3D du sous-sol pour une zone considérée. En d’autres termes, le volume de données sismiques afin d’obtenir une nouvelle image 3D complète. De ce fait, les techniques actuelles prennent beaucoup de temps, de l’ordre de plusieurs semaines à plusieurs mois de traitement informatique d’un volume conséquent de données, de l’ordre de plusieurs gigabits à plusieurs térabits de données. Les études 4D restent couteuses en ressources, en particulier en ressources logicielles face à l’afflux massif de données à traiter, ce qui se traduit par des temps de traitement longs, car il est nécessaire de générer des images 4D acquises à différentes périodes calendaires et d’en déduire une évolution de la structure géophysique d’une zone d’intérêt. Le demandeur a par le passé proposé divers améliorations des techniques connues de l’état de la technique, notamment pour réduire le coût en ressources financières, matérielles et logicielles de tels dispositifs, par exemple en proposant un procédé permettant d’identifier parmi un ensemble de sources et un ensemble de récepteurs, les récepteurs idéalement positionnés pour mesurer une évolution de la structure géophysique dans la zone d’intérêt. Malgré ces perfectionnements d’acquisition, l’enregistrement brut des données géophysiques acquises peut nécessiter une phase de traitement spécifique à l’extraction des données géophysiques utiles à la surveillance. L’invention vise à résoudre les inconvénients de l’état la technique en proposant un procédé permettant l’extraction des données géophysiques utiles pour une opération de surveillance d’un sous-sol à surveiller. PRESENTATION DE L’INVENTION Plus précisément, l’invention a pour objet un procédé de surveillance d’au moins une zone d’intérêt d’un sous-sol à surveiller comprenant un point d’intérêt et un voisinage dudit point d’intérêt, comprenant : une extraction d’un premier sous-ensemble de données géophysiques brutes associées au voisinage du point d’intérêt, parmi un ensemble de données géophysiques brutes acquises antérieurement, et traitées antérieurement pour générer une première image antérieure du sous-sol à surveiller ; une construction d’un premier modèle de projection à partir du premier sous-ensemble de données, le premier modèle de projection résultant d’un traitement de séparation de données géophysiques, pour réaliser l’extraction de données géophysiques utiles utilisées lors de la génération de la première image au voisinage du point d’intérêt ; et une réduction du premier sous-ensemble de données à un deuxième sous-ensemble comprenant un nombre minimal prédéterminé de données géophysiques brutes ; la création d’un premier ensemble de référence de données géophysiques utiles par application du premier modèle de projection au deuxième sous-ensemble. une acquisition de données géophysiques brutes formant un troisième sous-ensemble de données géophysiques brutes comprenant le même nombre minimal de données géophysiques brutes que le deuxième sous-ensemble, chaque donnée du troisième sous-ensemble correspondant à une donnée du deuxième sous-ensemble deux ; la création d’un premier ensemble de données géophysiques utiles de surveillance par application du premier modèle de projection au troisième sous-ensemble ; une comparaison du premier ensemble de référence et du premier ensemble de surveillance pour identifier une évolution du sous-sol à surveiller. Par la notion de « voisinage » du point d’intérêt, on entend une zone qui correspond à la résolution de l’image générée par la méthode géophysique utilisée au point d’intérêt. Par exemple, si la méthode géophysique est une méthode sismique, la résolution est communément évaluée au quart de la longueur d’onde, la longueur d’onde étant fonction de la vitesse de propagation de l’onde dans les formations géologiques et de la fréquence dominante du signal enregistré. Grâce à une telle combinaison de caractéristiques, il est ainsi possible de noter, à l’aide du procédé de surveillance, une évolution de la réponse géophysique de la zone d’intérêt par le biais d’une acquisition limitée de données, qui serait par exemple insuffisante à la construction d’une image 4D mais suffisante pour une détection d’éventuelles modifications de propriété géophysique. Ainsi, les ressources logicielles et le temps nécessaire au traitement des données géophysiques utiles sont sensiblement réduits. En outre, le modèle de projection construit est invariant au cours du temps, c’est-à-dire qu’une fois créé, il peut être réutilisé au cours du temps calendaire. Avantageusement, le procédé de surveillance comprend, après la création du premier ensemble de données géophysiques de surveillance : l’extraction d’une deuxième image parmi la première image, la deuxième image correspondant au voisinage du point d’intérêt ; la construction d’un deuxième modèle de projection à partir du premier ensemble de données géophysiques utiles de référence et de la deuxième image, de sorte que l’application du deuxième modèle de projection à l’ensemble de données géophysiques de référence corresponde sensiblement à la deuxième image pour créer un deuxième ensemble de données géophysiques utiles de référence ; l’application du deuxième modèle de projection au premier ensemble de données géophysiques utiles de surveillance pour créer un deuxième ensemble de données géophysiques utiles de surveillance ; une comparaison du deuxième ensemble de données géophysiques utiles de référence et du deuxième ensemble de données géophysiques de surveillance pour identifier une évolution du sous-sol à surveiller. Dans une telle configuration, la qualité des analyses produites par l’application de tels modèles de projection est sensiblement améliorée. Avantageusement, au moins une nouvelle itération de l’acquisition de données géophysiques brutes formant le troisième sous-ensemble et des étapes suivantes sont réalisées dans le temps. Ainsi, il est possible de noter une évolution de la structure géophysique au cours d’intervalles calendaires réguliers ou réguliers, de l’ordre par exemple de plusieurs semaines à plusieurs mois. Avantageusement, le nombre minimal de données géophysiques utiles est égal à un enregistrement. Dans une telle configuration, le volume de donnée à traiter est sensiblement réduit, ce qui réduit encore les ressources nécessaires à la surveillance d’une zone d’intérêt d’un sous-sol. Dans une telle configuration, l’enregistrement est par exemple sélectionné par le biais d’une méthode de tracé de rayons inverses depuis le point d’intérêt pour identifier la position source et récepteurs représentatif du point d’intérêt et qui ont généré l’enregistrement du premier sous-ensemble. Avantageusement, les données géophysiques brutes sont des enregistrements d’ondes sismiques et les données géophysiques utiles sont des ondes sismiques réfléchies, de surface ou réfractées. Avantageusement, l’acquisition de données d’ondes sismiques antérieures est une acquisition 2D, 3D ou 4D et en ce que l’acquisition des données géophysiques brutes utilisées pour former le troisième sous-ensemble est une autre acquisition 4D. Avantageusement, la projection sur le modèle de projection comprend au moins : une déconvolution et/ou filtrage en fréquence des données géophysiques ; une séparation des données géophysiques par application d’un filtre de pente défini dans le plan temps-distance et/ou dans le plan fréquence-nombre d’onde ; une correction statique et dynamique des données géophysiques. Dans une telle configuration, tout ensemble prédéfini de traitements physiques usuels, mono-canal ou multi-canaux, peuvent faire partie du modèle de projection. PRESENTATION DES FIGURES L’invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d’exemple, et se référant aux dessins annexés donnés à titre d’exemples non limitatifs, dans lesquels des références identiques sont données à des objets semblables et sur lesquels : La est un logigramme décrivant les étapes d’un procédé de surveillance d’un sous-sol à surveiller selon un premier mode de réalisation de l’invention ; La est un logigramme décrivant les étapes d’un procédé de surveillance d’un sous-sol à surveiller selon un autre mode de réalisation de l’invention ; La est un diagramme illustrant un sous-ensemble de données géophysiques brutes ; La est un diagramme illustrant le modèle de projection construit à l’aide de l’ensemble de données géophysiques brutes de la ; La est diagramme illustrant un enregistrement obtenu par réduction du sous-ensemble du diagramme de la auquel est appliqué successivement des opérations de traitement du signal ; La est une collection de diagrammes comparant un seul enregistrement du modèle de projection de la et l’enregistrement de la post-traitement ; La est un diagramme illustrant un sous-ensemble de données géophysiques brutes ainsi qu’une pluralité d’acquisitions formant chacune un sous-ensemble de données géophysiques brutes, chaque acquisition étant réalisée à une période calendaire différente ; La est un diagramme illustrant un ensemble de données géophysiques utiles de référence ainsi qu’une pluralité d’ensemble de données géophysiques utiles de surveillance issus de l’application du premier modèle de projection sur les données de la ; La est un diagramme illustrant les détections des différences d’arrivées des ondes par comparaison des données de la au point d’intérêt ; La est un diagramme illustrant un troisième ensemble de données géophysiques utiles issus de l’application du deuxième modèle de projection sur les données de la . Il faut noter que les figures exposent l’invention de manière détaillée pour permettre de mettre en œuvre l’invention ; bien que non limitatives, lesdites figures servent notamment à mieux définir l’invention le cas échéant. DESCRIPTION DETAILLEE DE L’INVENTION L’invention concerne un procédé de surveillance d’au moins une zone d’intérêt d’un sous-sol à surveiller comprenant un point d’intérêt et un voisinage dudit point d’intérêt par méthode géophysique. Dans le mode de réalisation illustré en , le procédé de surveillance est notamment mis en œuvre par méthode sismique, toute autre méthode géophysique étant applicable. Par la notion de « voisinage » du point d’intérêt, on entend une zone qui correspond à la résolution de l’image générée par la méthode géophysique utilisée au point d’intérêt. Par exemple, si la méthode géophysique est une méthode sismique, la résolution est communément peut être au quart de la longueur d’onde, la longueur d’onde étant fonction de la vitesse de propagation de l’onde dans les formations géologiques et de la fréquence dominante du signal enregistré. La présente suite de la description s’attache à décrire étape par étape le procédé de surveillance selon le premier mode de réalisation de la . Le procédé comprend une première étape d’extraction d’un premier sous-ensemble 1 de données géophysiques brutes associées au voisinage du point d’intérêt, parmi un ensemble 0 de données géophysiques brutes acquises antérieurement, et traitées antérieurement pour générer une première image Im1 antérieure du sous-sol à surveiller. Par le vocable « antérieurement », on entend que l’ensemble 0 de données géophysiques brutes a déjà été utilisé pour générer une première image Im1 du sous-sol au cours d’un procédé d’imagerie, par exemple 2D ou 3D. Par le vocable « 2D ou3D », on entend une imagerie réalisée selon deux ou trois dimensions spatiales. Ainsi, bien que l’ensemble 0 de données géophysiques brutes existe dans une version déjà été traité, par des opérations de traitement du signal, ledit ensemble 0 est exploité dans le cadre du procédé dans une version non traité, par exemple stocké sur un support de stockage. Le sous-ensemble 1 de données géophysiques brutes illustré sur la comprend au moins les données géophysiques brutes qui, lorsqu’elles ont été traitées antérieurement, ont contribué à une partie de la première image Im1 comprenant le point d’intérêt et son voisinage. Dans le mode de réalisation décrit, le sous-ensemble 1 est un regroupement d’enregistrements sismiques 2D brutes. Par le terme « 2D » on entend un enregistrement selon deux coordonnées spatiales. Le procédé comprend en outre une étape de construction d’un premier modèle de projection MOD1 à partir du premier sous-ensemble de données 1, le premier modèle de projection MOD1 résultant d’un traitement de séparation de données géophysiques, pour réaliser l’extraction de données géophysiques utiles utilisées lors de la génération de la première image Im1 au voisinage du point d’intérêt. Un modèle de projection est un modèle géophysique, en l’occurrence un modèle sismique permettant d’améliorer la qualité d’un sous-ensemble de données géophysiques brutes à un niveau au moins suffisant pour permettre une surveillance d’une zone d’intérêt en favorisant les données géophysiques utiles. Pour créer le premier modèle de projection MOD1, un premier traitement mono-canal est appliqué, par exemple une suppression des arrivées directes et/ou une déconvolution et/ou un filtrage en fréquence. Optionnellement, un traitement multi-canaux est appliqué après le ou les traitements mono-canaux. Le traitement multi-canaux peut par exemple comprendre une correction par décalage normal, plus connue de la personne du métier sous le terme anglais « normal moveout correction », un filtrage dans le plan fréquence, un filtrage dans le plan fréquence-nombre d’onde, un filtrage dans le plan fréquence-distance ou une décomposition en valeur singulière ou toute méthode de traitement du signal multi-canaux connue de l’état de la technique. Dans le cas d’une méthode sismique, le modèle de projection construit comprend des données géophysiques utiles, c’est-à-dire utiles à la génération de la première image Im1, par exemple des enregistrements sismiques d’ondes réfléchies, comme illustré sur la . Alternativement, les données géophysiques peuvent être des ondes sismiques réfractées ou des ondes sismiques de surface. Le premier modèle de projection MOD1 est stocké, par exemple sur un support de stockage pour des mises en œuvre ultérieures du procédé. Le procédé comprend en outre une étape de réduction du premier sous-ensemble 1 de données à un deuxième sous-ensemble 2 comprenant un nombre minimal prédéterminé de données géophysiques brutes. Dans le mode de réalisation décrit, le nombre minimal de données géophysiques est d’un seul enregistrement, communément appelé une « trace » par la personne du métier, comme illustré sur la . Ici, l’enregistrement est par exemple sélectionné par le biais d’une méthode de tracé de rayons inverses depuis le point d’intérêt pour identifier la position de l’ensemble source et récepteurs représentatif du point d’intérêt et qui ont généré l’enregistrement du premier sous-ensemble 1. Le procédé comprend en outre une étape de création d’un premier ensemble de référence de données géophysiques utiles B par application du premier modèle de projection MOD1 au deuxième sous-ensemble 2. Une application d’un modèle de projection. Cette opération est désignée comme « projection » du deuxième sous-ensemble 2 sur le premier modèle de projection MOD1. La projection peut par exemple être réalisée à partir d’un opérateur de projection, pour améliorer la qualité de l’enregistrement obtenu lors de l’étape précédente de tracé de rayons inverses. La illustre successivement, de gauche à droite, le premier modèle de projection MOD1, le premier ensemble de référence de données géophysiques utiles B, et une comparaison du premier modèle de projection MOD1 et de l’ensemble de données géophysiques de référence B. Ici, un coefficient de corrélation entre le premier modèle de projection MOD1 et l’ensemble de données géophysiques de référence B est calculé et considéré comme étant très haut, c’est-à-dire supérieur à 0,99. Avantageusement, le coefficient de corrélation est supérieur ou égal à 0.7. Le procédé comprend en outre une étape d’acquisition de données géophysiques brutes formant un troisième sous-ensemble de données géophysiques brutes 3 comprenant le même nombre minimal de données géophysiques brutes que le deuxième sous-ensemble 2, chaque donnée du troisième sous-ensemble 3 étant similaire à une donnée du deuxième sous-ensemble deux. Par le vocable « similaire », on entend que le troisième sous-ensemble de données géophysiques brutes 3 est issu d’une acquisition permettant la surveillance du même point d’intérêt ou son voisinage que le point d’intérêt que le deuxième sous-ensemble de données géophysiques brutes 2. Avantageusement, au moins une nouvelle itération de l’acquisition de données géophysiques brutes formant le troisième sous-ensemble 3 est réalisée dans le temps calendaire. La illustre un regroupement sur un même diagramme du deuxième sous-ensemble 2, situé à la position zéro, désignée date calendaire de référence, ainsi que cinq sous-ensembles 3 acquis à des dates calendaires différentes 31, 32, 33, 34 35, démarche désignée en anglais comme « Monitoring spot gathers » ou encore « MSG ». Les acquisitions peuvent se faire à intervalles réguliers ou irréguliers dans le temps calendaire. Le procédé comprend en outre la création d’un premier ensemble de données géophysiques utiles de surveillance M par application du premier modèle de projection MOD1 au troisième sous-ensemble 3. Dans le cas où l’acquisition du troisième sous-ensemble 3 est itérée dans le temps, chaque troisième sous-ensemble 3 est projeté sur le premier modèle de projection MOD1 de la même manière que le sous-ensemble de données géophysiques brutes 2 a été projeté antérieurement sur le premier modèle de projection MOD1. La illustre un MSG de l’ensemble de données géophysiques utiles de référence B situé à la position zéro, ainsi que cinq ensembles de données géophysiques utiles de surveillance M acquis à des dates calendaires différentes M1, M2, M3, M4, M5. Les acquisitions peuvent se faire à intervalles réguliers ou irréguliers dans le temps calendaire. Le procédé comprend en outre une comparaison du premier ensemble de référence B et du premier ensemble de surveillance M pour identifier une évolution du sous-sol à surveiller. Dans le cas où il y a une pluralité d’ensembles de données géophysiques utiles de surveillance M, chaque ensemble de surveillance est comparé un à un avec l’ensemble de référence base. Dans une telle configuration, il est ainsi possible de noter à l’aide du procédé de surveillance une évolution de la réponse géophysique de la zone d’intérêt par le biais d’une acquisition limitée de données, qui serait par exemple insuffisante à la construction d’une image 4D mais suffisante pour une détection d’éventuelles modifications de propriété géophysique. Ainsi, les ressources logicielles et le temps nécessaire au traitement des données géophysiques utiles sont sensiblement réduits. Avantageusement, plusieurs comparaisons utilisant des modèles théoriques différents peuvent être réalisées. La illustre ainsi une mesure d’une dérive temporelle entre l’ensemble de données géophysiques utiles de référence B et une pluralité d’ensemble de données géophysiques utiles de surveillance M acquis à des dates calendaires différentes. Pour chaque date, une différence de calcul de temps d’arrivée dDT de l’écho sismique du point d’intérêt entre l’ensemble de données géophysiques utiles de référence BASE et chacun des ensembles de données géophysiques utiles de surveillance M est calculée avec une incertitude associée dépendante du rapport signal à bruit et de la bande passante connu de la personne du métier. La suite de la description s’attache à décrire l’invention selon un autre mode de réalisation, illustré sur la . Dans ce mode de réalisation le procédé de surveillance comprend en outre, après la création du premier ensemble de données géophysiques de surveillance M : l’extraction d’une deuxième image Im2 parmi la première image Im1, la deuxième image Im2 correspondant au voisinage du point d’intérêt ; la construction d’un deuxième modèle de projection MOD2 à partir du premier ensemble de données géophysiques utiles de référence B et de la deuxième image Im2, de sorte que l’application du deuxième modèle de projection MOD2 à l’ensemble de données géophysiques de référence B corresponde sensiblement à la deuxième image IMAGE2 pour créer un deuxième ensemble de données géophysiques utiles de référence Bi ; l’application du deuxième modèle de projection MOD2 au premier ensemble de données géophysiques utiles de surveillance M pour créer un deuxième ensemble de données géophysiques utiles de surveillance Mi ; une comparaison du deuxième ensemble de données géophysiques utiles de référence Bi et du deuxième ensemble de données géophysiques de surveillance Mi pour identifier une évolution du sous-sol à surveiller. La illustre un MSG du deuxième ensemble de données géophysiques utiles de référence Bi situé à la position calendaire de référence, ainsi que cinq deuxièmes ensembles de données géophysiques utiles de surveillance Mi acquis à des dates calendaires différentes Mi1, Mi2, Mi3, Mi4, Mi5. A partir d’un tel MSG, il est alors possible de procéder à des analyses permettant de conclure sur l’évolution de la structure géophysique de la zone d’intérêt au cours du temps calendaire, notamment des mesures de dérives temporelles similaires à celles réalisées pour le mode de réalisation décrit sur la . On notera par ailleurs que l’invention n’est pas limitée aux modes de réalisation décrit précédemment. Il apparaîtra en effet à la personne du métier que diverses modifications peuvent être apportées au mode de réalisation décrit ci-dessus, à la lumière de l’enseignement qui vient de lui être divulgué. Dans la présentation détaillée de l’invention qui est faite précédemment, les termes utilisés ne doivent pas être interprétés comme limitant l’invention au mode de réalisation exposé dans la présente description, mais doivent être interprétés pour y inclure tous les équivalents dont la prévision est à la portée de la personne du métier en appliquant ses connaissances générales à la mise en œuvre de l’enseignement qui vient de lui être divulgué. Procédé de surveillance d’au moins une zone d’intérêt d’un sous-sol à surveiller comprenant un point d’intérêt et un voisinage dudit point d’intérêt, comprenant : une extraction d’un premier sous-ensemble (1) de données géophysiques brutes associées au voisinage du point d’intérêt, parmi un ensemble (0) de données géophysiques brutes acquises antérieurement, et traitées antérieurement pour générer une première image (Im1) antérieure du sous-sol à surveiller ; une construction d’un premier modèle de projection (MOD1) à partir du premier sous-ensemble de données (1), le premier modèle de projection (MOD1) résultant d’un traitement de séparation de données géophysiques, pour réaliser l’extraction de données géophysiques utiles utilisées lors de la génération de la première image (Im1) au voisinage du point d’intérêt ; et une réduction du premier sous-ensemble (1) de données à un deuxième sous-ensemble (2) comprenant un nombre minimal prédéterminé de données géophysiques brutes ; la création d’un premier ensemble de référence de données géophysiques utiles (B) par application du premier modèle de projection (MOD1) au deuxième sous-ensemble (2) ; une acquisition de données géophysiques brutes formant un troisième sous-ensemble de données géophysiques brutes (3) comprenant le même nombre minimal de données géophysiques brutes que le deuxième sous-ensemble (2), chaque donnée du troisième sous-ensemble (3) correspondant à une donnée du deuxième sous-ensemble (2) ; la création d’un premier ensemble de données géophysiques utiles de surveillance (M) par application du premier modèle de projection (MOD1) au troisième sous-ensemble (3) ; une comparaison du premier ensemble de référence (B) et du premier ensemble de surveillance (M) pour identifier une évolution du sous-sol à surveiller. Procédé de surveillance selon la revendication 1, caractérisé en ce qu’il comprend, après la création du premier ensemble de données géophysiques de surveillance (M) : l’extraction d’une deuxième image (Im2) parmi la première image (Im1), la deuxième image (Im2) correspondant au voisinage du point d’intérêt ; la construction d’un deuxième modèle de projection (MOD2) à partir du premier ensemble de données géophysiques utiles de référence (B) et de la deuxième image (Im2), de sorte que l’application du deuxième modèle de projection (MOD2) à l’ensemble de données géophysiques de référence (B) corresponde sensiblement à la deuxième image (Im2) pour créer un deuxième ensemble de données géophysiques utiles de référence (Bi) ; l’application du deuxième modèle de projection (MOD2) au premier ensemble de données géophysiques utiles de surveillance (M) pour créer un deuxième ensemble de données géophysiques utiles de surveillance (Mi) ; une comparaison du deuxième ensemble de données géophysiques utiles de référence (Bi) et du deuxième ensemble de données géophysiques de surveillance (Mi) pour identifier une évolution du sous-sol à surveiller. Procédé de surveillance selon l’une quelconque des revendications précédentes, caractérisé en ce qu’au moins une nouvelle itération de l’acquisition de données géophysiques brutes formant le troisième sous-ensemble (3) et des étapes suivantes sont réalisées dans le temps. Procédé de surveillance selon l’une quelconque des revendications précédentes, caractérisé en ce que le nombre minimal de données géophysiques utiles est égal à un enregistrement. Procédé de surveillance selon l’une quelconque des revendications précédentes, caractérisé en ce que les données géophysiques brutes sont des enregistrements d’ondes sismiques et les données géophysiques utiles sont des ondes sismiques réfléchies, de surface ou réfractées. Procédé de surveillance selon la revendication 5, caractérisé en ce que l’acquisition de données d’ondes sismiques antérieures est une acquisition 2D, 3D ou 4D et en ce que l’acquisition des données géophysiques brutes utilisées pour former le troisième sous-ensemble (3) est une autre acquisition 4D. Procédé de surveillance selon l’une quelconque des revendications précédentes, caractérisé en ce que l’application du premier modèle de projection (MOD1) comprend au moins : une déconvolution ou filtrage en fréquence des données géophysiques ; une séparation des données géophysiques par application d’un filtre de pente défini dans le plan temps-distance et/ou dans le plan fréquence-nombre d’onde ; une correction statique et dynamique des données géophysiques.