FR 3084473 A1 20200131 FR 1856857 A 20180724 PROCEDE ET DISPOSITIF DE SURVEILLANCE DU SOUS-SOL TERRESTRE SOUS UNE ZONE CIBLE La présente invention se rapporte à un procédé et un dispositif de 5 surveillance du sous-sol terrestre sous une zone cible. Elle appartient au domaine de l'analyse sismique et trouve application notamment dans les activités de surveillance du sous-sol terrestre et de détection d'anomalies dans ce dernier, comme la détection de cavités sous des voies ferrées, des routes, des bâtiments ou encore des aéroports. 10 Afin de garantir la sécurité des utilisateurs de ces infrastructures terrestres, il est en effet important de procéder à une surveillance instantanée, temporaire ou permanente du sous-sol, afin de détecter toute éventuelle anomalie souterraine telle qu'une cavité, susceptible de fragiliser les ouvrages construits en surface. 15 Dans un article intitulé « Seismic interferometry of railroad induced ground motions: body and surface wave imaging », publié dans Geophys. J. Int. 205, 2016, pages 301-313, Diego A. QUIROS et al. ont proposé une méthode de génération de trajets virtuels de propagation le long d'une voie ferrée, à partir de l'analyse interférométrique du bruit ambiant, capté par une 20 ligne unique de récepteurs ou capteurs d'ondes sismiques déployée le long de cette voie ferrée. Cette méthode présente l'inconvénient de limiter l'information obtenue à une droite parallèle à la voie ferrée. Elle ne permet pas de produire une image du sous-sol sous la voie ferrée. 25 De plus, dans cette méthode purement passive, le bruit des trains est mélangé à tous les autres bruits et utilisé lorsque ceux-ci sont éloignés des récepteurs afin d'être alignés avec la ligne de récepteurs. Le bruit capté par les récepteurs est donc faible et peu propice à l'imagerie. La présente invention a pour but de remédier aux inconvénients de 30 l'art antérieur. Dans ce but, la présente invention propose un procédé de surveillance du sous-sol terrestre sous une zone cible, remarquable en ce qu'il comporte des étapes suivant lesquelles : 2 on enregistre des ondes sismiques provenant d'une source mobile de bruit sismique identifiée, au moyen d'au moins un couple de capteurs disposés de part et d'autre de la zone cible ; on sélectionne des périodes temporelles correspondant aux 5 alignements des couples de capteurs avec la source de bruit ; on reconstruit par interférométrie un sismogramme de la zone cible à partir des ondes sismiques enregistrées et des périodes temporelles sélectionnées ; et on engendre une image du sous-sol sous la zone cible à partir du 10 sismogramme. Ainsi, le fait d'utiliser, d'une part, une source de bruit sismique mobile dont les caractéristiques sont connues grâce à l'analyse du signal lui-même et ne nécessitent pas d'information a priori et d'autre part, au moins deux capteurs qui pendant une certaine période de temps vont être alignés avec 15 cette source de bruit, permet une surveillance passive dans l'espace et le temps, c'est-à-dire à quatre dimensions, de grande amplitude et à haute résolution, de la partie du sous-sol située sous la surface où sont disposés les capteurs, à la condition que la source de bruit mobile soit alignée avec le couple de capteurs à au moins un instant donné. 20 Ainsi, si un couple de capteurs est disposé de part et d'autre d'une zone cible (par exemple d'une voie de chemin de fer), il est possible d'engendrer une image du sous-sol sous cette zone cible si une source mobile (par exemple un train sur une autre voie) passe à proximité et s'aligne à au moins un instant donné avec le couple de capteurs. 25 Dans un mode particulier de réalisation, le procédé comporte en outre, préalablement à l'étape de reconstruction, une étape de prétraitement des ondes sismiques enregistrées. Cela permet d'obtenir des signaux de qualité optimale, ce qui accroît la fiabilité de l'image du sous-sol obtenue. 30 Avantageusement, l'étape de prétraitement comporte des opérations de débruitage et de filtrage en fréquence. De façon optionnelle, l'étape de prétraitement peut comporter une opération de blanchiment spectral et/ou une opération de binarisation. 3 Ce type d'opération permet de simplifier les traitements ultérieurs de reconstruction du sismogramme. Dans un mode particulier de réalisation, lors de l'étape de reconstruction par interférométrie, on effectue au moins une opération 5 d'intercorrélation de signaux issus des ondes sismiques respectivement enregistrées par les deux capteurs du couple de capteurs. Cela permet de reconstruire efficacement le sismogramme. Dans un mode particulier de réalisation, l'étape d'enregistrement est effectuée de façon continue dans le temps. 10 Cela permet de produire des états du sous-sol à une périodicité choisie, par exemple quotidienne ou hebdomadaire ou mensuelle. Dans un mode particulier de réalisation, le procédé comporte en outre, à la suite de l'étape de reconstruction, une étape d'addition de signaux reconstruits du sismogramme provenant de plusieurs sources de bruit mobiles 15 respectivement alignées avec le couple de capteurs pendant des périodes de temps distinctes. Cette étape d'addition permet de renforcer la qualité et la fiabilité de la reconstruction. Dans un mode particulier de réalisation, lors de l'étape 20 d'enregistrement, on enregistre des ondes sismiques de surface. Dans ce mode de réalisation, lors de l'étape de génération de l'image : on détermine des courbes de dispersion des ondes de surface dans les signaux reconstruits du sismogramme dans une bande de fréquences prédéterminée ; 25 on réalise une tomographie de la vitesse des signaux reconstruits à partir des courbes de dispersion ; et on effectue une inversion des vitesses d'ondes de surface des signaux reconstruits, de façon à produire un modèle de vitesse d'ondes Pet S. En variante, lors de l'étape d'enregistrement, on peut enregistrer des 30 ondes sismiques réfractées. Dans cette variante, lors de l'étape de génération de l'image : on détermine les temps d'arrivée des signaux reconstruits du sismogramme correspondant aux ondes sismiques réfractées ; et 4 on réalise une tomographie des signaux reconstruits correspondant aux ondes réfractées P à partir des temps d'arrivée précités, de façon à produire un modèle de vitesse d'ondes R Selon des caractéristiques particulières possibles, les capteurs peuvent 5 comporter au moins un géophone et/ou au moins un accéléromètre et/ou au moins un capteur basé sur l'utilisation de fibres optiques, par exemple de type DAS (capteur acoustique réparti, en anglais « Distributed Acoustic Sensing La zone cible précitée peut être une voie ferrée ou une route ou un bâtiment. 10 Pour une zone cible donnée, les couples de capteurs disposés de part et d'autre de la zone cible peuvent être déployés en groupes suivant différents agencements, comme par exemple suivant des lignes droites ou courbes, ou sur des grilles régulières ou non de forme carrée ou hexagonale ou quelconque et peuvent être espacés régulièrement ou non, voire être placés 15 aléatoirement. Dans le même but que celui indiqué plus haut, la présente invention propose également un dispositif de surveillance du sous-sol terrestre sous une zone cible, remarquable en ce qu'il comporte des éléments adaptés à mettre en oeuvre des étapes d'un procédé tel que succinctement décrit ci-dessus. 20 Les caractéristiques particulières et les avantages du dispositif étant similaires à ceux du procédé, ils ne sont pas répétés ici. D'autres aspects et avantages de l'invention apparaîtront à la lecture de la description détaillée ci-après de modes particuliers de réalisation, donnés à titre d'exemples nullement limitatifs, en référence aux dessins annexés, 25 dans lesquels : - la figure 1 est un organigramme illustrant des étapes d'un procédé de surveillance du sous-sol terrestre sous une zone cible conforme à la présente invention, dans un mode particulier de réalisation ; - la figure 2 est une vue schématique illustrant une configuration 30 d'enregistrement d'ondes sismiques dans un exemple nullement limitatif où la zone cible est une voie ferrée ; - la figure 3 est une vue schématique illustrant une configuration de reconstruction de trajets virtuels de propagation dans l'exemple de la figure 2; - la figure 4 est une vue schématique illustrant une configuration de 5 reconstruction de trajets virtuels de propagation suivant plusieurs azimuts, dans l'exemple de la figure 2 ; et - la figure 5 est une vue schématique illustrant une configuration de reconstruction de trajets virtuels de propagation suivant plusieurs azimuts et en utilisant plusieurs voies ferrées. Comme le montre l'organigramme de la figure 1, au cours d'une première étape 10 du procédé de surveillance du sous-sol terrestre sous une zone cible conforme à l'invention, au moins pendant une période de temps prédéterminée, on enregistre des ondes sismiques. Ces ondes sismiques proviennent d'une source mobile de bruit 15 sismique, c'est-à-dire une source en déplacement et dont le déplacement produit du bruit. Cette source de bruit est connue au sens où elle est bien identifiée. Il ne s'agit donc pas d'un bruit ambiant d'origine non identifiée. L'enregistrement de ces ondes sismiques est réalisé au moyen d'au moins un couple de capteurs disposés de façon à être alignés avec la source 20 de bruit pendant la période de temps précitée. Les couples de capteurs peuvent être disposés selon différents agencements, comme par exemple régulièrement espacés le long de lignes droites comme le montre la figure 2 décrite ci-après, mais ils peuvent aussi être espacés de façon irrégulière, suivant toutes sortes de courbes, ou sur des 25 grilles régulières ou non de forme carrée ou hexagonale ou quelconque. Ils peuvent aussi être placés de façon aléatoire. La figure 2 donne à titre d'exemple nullement limitatif une configuration d'enregistrement dans le cas où la zone cible est une voie ferrée 20. Néanmoins, l'invention s'applique tout aussi bien à une zone cible 30 constituée par une route, ou encore un bâtiment ou une autre infrastructure, comme par exemple un aéroport. Une pluralité de récepteurs ou capteurs d'ondes sismiques 22 est déployée le long de la voie ferrée 20, dite voie cible, de part et d'autre de 6 celle-ci. Ces capteurs 22 sont désignés sur le dessin par la mention « Ligne récepteur A » et « Ligne récepteur B ». Les capteurs 22 peuvent par exemple comporter au moins un géophone et/ou au moins un accéléromètre et/ou au moins un capteur basé sur l'utilisation de fibres optiques, par exemple du type 5 DAS (capteur acoustique réparti, en anglais « Distributed Acoustic Sensing Sur chaque côté de la voie ferrée cible 20, les capteurs 22 peuvent par exemple être disposés à intervalles réguliers les uns des autres, par exemple à quelques mètres d'intervalle, en laissant par exemple 2 mètres entre chaque capteur 22. 10 Les capteurs 22 peuvent être placés à la surface du sol ou être légèrement enterrés, c'est-à-dire placés à quelques mètres sous la surface du sol. Ils peuvent être déployés sur ou sous le ballast, ou bien sur ou sous les rails ou encore dans un tunnel ferroviaire. La pluralité de capteurs 22 comporte au moins un couple de capteurs 15 220, 221 disposés en regard l'un de l'autre, de part et d'autre de la voie ferrée 20, de façon transversale à celle-ci. La source de bruit sismique mobile est un train 24 circulant sur une voie ferrée 21, dite voie source, distincte de la voie ferrée cible 20. Le bruit émis par le train 24 en mouvement est symbolisé sur le dessin par des flèches 20 en forme de cercles 26. Ce bruit provient principalement du contact des roues du train 24 sur les rails de la voie ferrée 21. L'ensemble comportant la pluralité de capteurs 22, la source 24 mobile de bruit sismique et un module adapté à réaliser les traitements du procédé de surveillance du sous-sol décrit ici, forme un dispositif de surveillance du sous-25 sol conforme à l'invention. Si on considère le couple de capteurs 220, 221, il va se trouver dans l'alignement du train 24 pendant une période de temps prédéterminée définie par un premier instant où le point de contact de la première roue du train 24 avec le rail est sur la droite qui passe par le barycentre des capteurs 220 et 30 221 et un second instant où le point de contact de la dernière roue du train 24 avec le rail est sur la droite qui passe par le barycentre des capteurs 220 et 221. 7 Conformément à l'invention, pendant cette période de temps, les capteurs 220 et 221 enregistrent des ondes sismiques qui vont faire l'objet des traitements ultérieurs du procédé conforme à l'invention. Cette période de temps est prédéterminée, étant donné qu'on connaît les paramètres 5 cinématiques du train 24, à savoir, sa position et sa vitesse à chaque instant. En ce sens, la source de bruit sismique mobile que constitue le train 24 est identifiée. L'étape 10 d'enregistrement peut être limitée à cette période de temps prédéterminée. En variante, l'enregistrement peut avoir lieu de façon continue 10 dans le temps. En référence à l'organigramme de la figure 1, de façon optionnelle, l'étape 10 d'enregistrement des ondes sismiques peut être suivie d'une étape 11 de prétraitements des ondes sismiques enregistrées par les capteurs 22. Ces prétraitements peuvent être de divers types et incluent au moins une phase de débruitage et une phase de filtrage en fréquence. Dans l'exemple non limitatif décrit ici d'une voie ferrée, lors de la phase de débruitage des ondes sismiques enregistrées, on enlève le bruit principal environnant du bruit sismique produit par le train 24. Ce bruit principal environnant est un bruit électrique lié aux caténaires dont on connaît les caractéristiques. On obtient ainsi un signal avec un rapport signal sur bruit amélioré. Lors de la phase de filtrage en fréquence, on filtre le signal dans une bande de fréquences prédéterminée qui correspond à la fois à la bande d'émission estimée du train 24 et à la bande de fréquence d'intérêt, laquelle dépend de la profondeur d'intérêt de l'investigation du sous-sol sous la zone cible que constitue la voie ferrée 20. En outre, entre la phase de débruitage et la phase de filtrage en fréquence, de façon optionnelle, on peut en outre effectuer une opération de blanchiment spectral et/ou une opération de binarisation. L'opération de blanchiment spectral, appelé aussi égalisation spectrale, technique connue en soi, consiste à pondérer toutes les composantes fréquentielles du signal de telle façon qu'elles aient toutes la même représentation énergétique. 8 L'opération de binarisation, également connue en soi, consiste à attribuer une valeur -1 ou +1 à chaque valeur instantanée du signal, pour simplifier le signal. D'autres types de prétraitements peuvent être appliqués. 5 A l'issue de l'étape 10 d'enregistrement et de l'étape 11 éventuelle de prétraitement des ondes sismiques enregistrées par les capteurs 22, une étape 12 de reconstruction d'un sismogramme de la zone cible est effectuée. Cette étape 12 de reconstruction utilise les ondes sismiques enregistrées et éventuellement prétraitées et leur applique une technique 10 d'interférométrie. Cette technique d'interférométrie comporte par exemple au moins une opération d'intercorrélation de signaux issus des ondes sismiques respectivement enregistrées par les capteurs 220 et 221 pendant la période de temps prédéterminée où le train 24 passe sur la voie source 21. Cette 15 technique d'interférométrie permet de reconstruire des trajets virtuels de propagation perpendiculaires à la voie ferrée cible 20. Ces trajets virtuels de propagation caractérisent le sismogramme. La figure 3 illustre schématiquement ces trajets virtuels de propagation pour chaque couple de capteurs en regard l'un de l'autre de part et d'autre de 20 la voie ferrée cible 20. En variante, toute autre technique d'interférométrie pourrait être appliquée lors de cette étape 12 de reconstruction, comme par exemple une méthode à base d'opérations de convolution en lieu et place des opérations d'intercorrélation susmentionnées. 25 Comme illustré par la figure 4, l'étape 12 de reconstruction peut être effectuée pour différents azimuts. En effet, on peut effectuer des opérations d'intercorrélation, non seulement entre les signaux issus des ondes sismiques enregistrées par les couples de capteurs situés en regard l'un de l'autre comme le couple de capteurs 220 et 221, mais également entre les couples de 30 capteurs disposés de façon analogue aux capteurs A5 et B6 ou aux capteurs A7 et B6, c'est-à-dire les couples de capteurs pour lesquels la droite qui passe par les barycentres des deux capteurs du couple forme un angle non nul avec la perpendiculaire à la voie ferrée cible 20. 9 Par ailleurs, comme illustré par la figure 5, l'étape 12 de reconstruction peut être effectuée non seulement pour différents azimuts, mais en outre en utilisant plusieurs voies ferrées A à E. La voie source 21, parallèle aux voies 1 à 4, est située en regard de la « ligne récepteur A » comme sur les figures 2 à 5 4 mais pour simplifier, n'est pas représentée sur la figure 5. Ainsi, lorsqu'un train passe sur la voie source 21, on peut reconstruire par interférométrie des trajets virtuels de propagation entre les voies A et B et/ou entre les voies B et C et/ou entre les voies C et D et/ou entre les voies D et E. Autrement dit, pour chaque couple de voies ferrées adjacentes, chaque 10 voie ferrée permet une reconstruction d'un sismogramme de la zone située sous la voie ferrée adjacente, en utilisant les couples de capteurs déployés de part et d'autre de cette voie ferrée adjacente. On a décrit ci-dessus la reconstruction de trajets virtuels de propagation transverses, c'est-à-dire transversaux à la voie ferrée, en utilisant 15 la fenêtre temporelle dans laquelle le train 24 est aligné avec un couple de capteurs du type du couple 220, 221. En variante, il est possible d'utiliser la corrélation entre deux capteurs de part et d'autre d'une voie et le signal émis par un train sur cette même voie. Ainsi, les variations temporelles des temps de propagation peuvent être 20 utilisées en surveillance permanente, à des fins de mesure d'apparition d'anomalies. Néanmoins, on peut aussi reconstruire des trajets virtuels de propagation longitudinaux, c'est-à-dire parallèles à la voie ferrée. Pour ce faire, il convient d'utiliser des fenêtres temporelles de durée prédéterminée 25 précédant et suivant le passage du train 24, respectivement après que les capteurs commencent à recevoir des ondes sismiques et avant que les capteurs cessent de recevoir des ondes sismiques. En référence à la figure 1, à la suite de l'étape 12 de reconstruction, on peut, de façon optionnelle, effectuer une étape 13 d'addition ou 30 « empilement » (en anglais « stacking ») de signaux reconstruits provenant d'ondes sismiques captées en provenance de plusieurs sources 24 de bruit sismiques mobiles. 10 Dans l'exemple non limitatif décrit ici où ces sources de bruit sont des trains, plusieurs trains circulant à des moments différents, ayant des longueurs et des vitesses différentes, qui vont donc être alignés avec un couple de capteurs du type du couple 220, 221 pendant des périodes de temps 5 distinctes, vont fournir des signaux reconstruits similaires. L'addition de ces signaux multiples permet d'améliorer le rapport signal sur bruit et donc, d'obtenir in fine une image de meilleure qualité du sous-sol sous la zone cible. Par ailleurs, cela permet de produire, pour un couple donné de capteurs du type du couple 220, 221, des trajets virtuels de 10 propagation de façon ponctuelle ou périodique à une périodicité choisie, à savoir, quotidienne ou hebdomadaire ou encore mensuelle, par exemple. Comme le montre la figure 1, à l'issue de l'étape 12 de reconstruction et de l'éventuelle étape 13 d'addition, on effectue une étape 14 de génération d'une image du sous-sol sous la zone cible, à partir du sismogramme 15 reconstruit précédemment. Lors de cette étape 14 de génération d'une image, les traitements diffèrent selon le type d'ondes sismiques enregistrées à l'étape 10. En effet, on peut enregistrer, soit des ondes de surface, soit des ondes réfractées, soit éventuellement des ondes réfléchies. 20 Dans un mode de réalisation où à l'étape 10, on enregistre des ondes sismiques de surface, on effectue les opérations suivantes pour obtenir une image du sous-sol sous la zone cible. On détermine tout d'abord des courbes de dispersion des ondes de surface dans les signaux reconstruits, dans une bande de fréquences 25 prédéterminée. A titre d'exemple non limitatif, cette bande de fréquences peut être la bande de 1 à 100 Hz. De façon optionnelle, on peut ensuite sélectionner uniquement les maxima de ces courbes de dispersion. Puis à partir des courbes de dispersion, on réalise une tomographie de 30 la vitesse des signaux reconstruits. Enfin, on effectue une inversion des vitesses d'ondes de surface des signaux reconstruits issus des ondes de surface, de façon à produire un 11 modèle de vitesse des ondes de pression et de cisaillement, connues de l'homme du métier sous la dénomination d'ondes Pet S. Grâce à la répétition de ces opérations dans le temps, on obtient ainsi une image du sous-sol sous la zone cible qui évolue dans le temps, autrement 5 dit, une représentation à quatre dimensions de cette zone du sous-sol. Dans un autre mode de réalisation, où à l'étape 10, on enregistre des ondes sismiques réfractées, on effectue les opérations suivantes pour obtenir une image du sous-sol sous la zone cible. On détermine tout d'abord les temps d'arrivée des signaux reconstruits 10 du sismogramme, correspondant aux ondes sismiques réfractées. Puis on réalise une tomographie des signaux reconstruits correspondant aux ondes réfractées P à partir de ces temps d'arrivée, de façon à produire un modèle de vitesse des ondes de pression, dites ondes P. De même que dans le mode de réalisation où on enregistre des ondes 15 sismiques de surface, grâce à la répétition de ces opérations dans le temps, on obtient ainsi une image du sous-sol sous la zone cible qui évolue dans le temps, autrement dit, une représentation à quatre dimensions de cette zone du sous-sol. Dans encore un autre mode particulier de réalisation, on pourrait 20 enregistrer à l'étape 10 des ondes sismiques réfléchies. Cela conduirait alors à obtenir une image de la réflectivité du sous-sol sous la zone cible. Le procédé conforme à l'invention comporte une surveillance passive du sous-sol sous une zone cible au sens où il traite des ondes sismiques reçues, mais n'en produit pas. Néanmoins, on peut envisager de combiner les 25 capteurs qui reçoivent les ondes sismiques avec des sources sismiques actives, REVENDICATIONS1. Procédé de surveillance du sous-sol terrestre sous une zone cible, caractérisé en ce qu'il comporte des étapes suivant lesquelles : on enregistre (10) des ondes sismiques provenant d'une source (24) mobile de bruit sismique identifiée, au moyen d'au moins un couple (220, 221) de capteurs (22) disposés de part et d'autre de ladite zone cible ; on sélectionne des périodes temporelles correspondant aux alignements des couples (220, 221) de capteurs (22) avec ladite source (24) de bruit ; on reconstruit (12) par interférométrie un sismogramme de la zone cible à partir desdites ondes sismiques enregistrées et des périodes temporelles sélectionnées ; et on engendre (14) une image du sous-sol sous la zone cible à partir dudit sismogramme. 2. Procédé selon la revendication 1, caractérisé en ce qu'il comporte en outre, préalablement à l'étape (12) de reconstruction, une étape (11) de prétraitement des ondes sismiques enregistrées. 3. Procédé selon la revendication 2, caractérisé en ce que l'étape (11) de prétraitement comporte des opérations de débruitage et de filtrage en fréquence. 4. Procédé selon la revendication 2 ou 3, caractérisé en ce que l'étape (11) de prétraitement comporte une opération de blanchiment spectral. 5. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que lors de l'étape (12) de reconstruction par interférométrie, on effectue au moins une opération d'intercorrélation de signaux issus des ondes sismiques respectivement enregistrées par les deux capteurs dudit couple (220, 221) de capteurs. 13 6. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'étape (10) d'enregistrement est effectuée de façon continue dans le temps. 7. Procédé selon l'une quelconque des revendications précédentes, 5 caractérisé en ce qu'il comporte en outre, à la suite de l'étape (12) de reconstruction, une étape (13) d'addition de signaux reconstruits dudit sismogramme provenant de plusieurs sources de bruit mobiles respectivement alignées avec ledit couple (220, 221) de capteurs pendant des périodes de temps distinctes. 10 8. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que lors de l'étape (10) d'enregistrement, on enregistre des ondes sismiques de surface. 9. Procédé selon la revendication 8, caractérisé en ce que lors de l'étape (14) de génération de ladite image : 15 on détermine des courbes de dispersion des ondes de surface dans les signaux reconstruits dudit sismogramme dans une bande de fréquences prédéterminée ; on réalise une tomographie de la vitesse desdits signaux reconstruits à partir desdites courbes de dispersion ; et 20 on effectue une inversion des vitesses d'ondes de surface desdits signaux reconstruits, de façon à produire un modèle de vitesse d'ondes P et S. 10.Procédé selon l'une quelconque des revendications 1 à 7, caractérisé en ce que lors de l'étape (10) d'enregistrement, on enregistre des ondes 25 sismiques réfractées. 11.Procédé selon la revendication 10, caractérisé en ce que lors de l'étape (14) de génération de ladite image : 14 on détermine les temps d'arrivée des signaux reconstruits dudit sismogramme correspondant auxdites ondes sismiques réfractées ; et on réalise une tomographie desdits signaux reconstruits correspondant aux ondes réfractées P à partir desdits temps d'arrivée, de façon à 5 produire un modèle de vitesse d'ondes P. 12. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que lesdits capteurs (22) comportent au moins un géophone et/ou au moins un accéléromètre et/ou au moins un capteur basé sur l'utilisation de fibres optiques. 10 13.Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que ladite zone cible est une voie ferrée (20) ou une route ou un bâtiment. 14. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que lesdits capteurs (22) disposés de part et d'autre 15 de ladite zone cible sont placés suivant des lignes droites ou courbes, ou sur des grilles régulières ou non de forme carrée ou hexagonale ou quelconque, sont espacés régulièrement ou non, ou sont placés aléatoirement. 15. Dispositif de surveillance du sous-sol terrestre sous une zone cible, 20 caractérisé en ce qu'il comporte des moyens adaptés à mettre en oeuvre des étapes d'un procédé selon l'une quelconque des revendications précédentes.