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.