Brevet FR3072185A1 — publié le 20190412

Arrière-plan de l'invention La présente invention concerne les capteurs vélocimétriques à courte période, typiquement de l'ordre de 0,5 à 10 s, adaptés notamment pour la surveillance microsismique, par exemple de réservoirs ou de puits profonds, la mesure de bruit sismique, la surveillance de la santé d'un bâtiment et la surveillance des tremblements de terre. Il existe à ce jour différents types de capteurs vélocimétriques adaptés pour la mesure des vibrations de courte période. On compte parmi eux les capteurs à balance de force à contre réaction, qui offrent une large bande passante mais sont très couteux (de l'ordre de 12 000 à 15 000 EUR), volumineux, sensibles et fragiles. Il existe aussi des géophones de moindre coût (quelques dizaines d'euros), plus compacts et plus robustes. Cependant, leur fréquence propre est de l'ordre de 4,5 Hz et donc insuffisante pour la mesure des périodes envisagées. Le document DE 3307575 Al comporte un géophone et un amplificateur jouant le rôle d'une résistance négative là où on place en général une résistance positive pour amortir le géophone. La valeur du gain de l'amplificateur est définie par les résistances branchées aux bornes de l'amplificateur opérationnel et leurs valeurs sont déterminées en fonction de l'amortissement m et de la fréquence propre fo intrinsèques du géophone. La qualité de la réponse spectrale du capteur ainsi obtenu pour une fréquence propre comprise entre une fréquence propre de consigne et la fréquence propre intrinsèque dépend donc de la méthode de détermination de l'amortissement m et de la fréquence propre fo intrinsèques du géophone et de la valeur des résistances choisies. De plus, le capteur ainsi obtenu ne peut être utilisé que dans des conditions d'utilisation nominales, la réponse du capteur n'étant pas garantie notamment pour des températures et/ou des inclinaisons différentes des conditions nominales d'utilisation du capteur. Objet et résumé de l'invention La présente invention vise notamment à pallier les inconvénients de l'art antérieur précités. En particulier, un but de l'invention est d'améliorer la performance des 2 capteurs vélocimétriques à courte période comprenant des géophones dont la réponse spectrale est élargie. À cet effet, la présente invention propose, selon un premier aspect, un procédé de filtrage d'un signal de réponse d'un géophone, ledit géophone 5 comportant un aimant entouré d'une bobine montée mobile par rapport à l'aimant, la bobine étant apte à se déplacer suite à des vibrations d'un environnement de ladite bobine, le géophone présentant un amortissement intrinsèque et une fréquence propre intrinsèque. Le procédé comporte: 10 - la détermination de l'amortissement intrinsèque et de la fréquence propre intrinsèque du géophone, comprenant : o l'alimentation de la bobine avec un courant d'excitation pendant une durée prédéterminée, o la détermination, à l'aide du signal de réponse du 15 géophone mesuré à la disparition du courant d'excitation, de la fréquence propre intrinsèque et de l'amortissement intrinsèque du géophone, - la détermination d'un filtre numérique configuré pour modifier une fréquence propre du géophone d'une fréquence propre intrinsèque à une 20 fréquence propre de consigne, la fréquence propre de consigne étant inférieure à la fréquence propre intrinsèque et ajuster le gain du géophone de telle sorte que ledit gain soit sensiblement constant pour des fréquences comprises entre la fréquence propre de consigne et la fréquence propre intrinsèque, - l'application dudit filtre sur ledit signal de réponse du géophone. 25 Le procédé de filtrage proposé comprend une nouvelle méthode de détermination de l'amortissement et de la fréquence propre intrinsèques associée à l'utilisation d'un filtre numérique. Cette nouvelle méthode se base sur la mesure du signal de réponse du géophone après la disparition du courant d'excitation stationnaire alimentant la bobine. Le courant d'excitation 30 est utilisé pour déplacer la bobine par rapport à l'aimant. Contrairement aux dispositifs de l'art antérieur, la réponse du géophone est donc mesurée lorsque la bobine n'est plus soumise au courant d'excitation. La réponse du géophone mesurée correspond donc bien à la réponse du géophone à une excitation 3 mécanique, c'est-à-dire à un déplacement de la bobine et à la réponse du géophone suite à ce déplacement. L'amortissement et la fréquence propre intrinsèques sont donc mesurés avec plus de précision. De plus, l'utilisation d'un filtre numérique permet d'adapter la fonction 5 de filtrage aux paramètres mesurés, à savoir l'amortissement et la fréquence propre intrinsèques avec plus de précision et de flexibilité par rapport à la fonction de transfert du filtre souhaitée. Dans des modes de réalisation, le procédé de filtre comprend l'une ou 10 plusieurs des caractéristiques suivantes : - la fréquence propre de consigne est comprise entre 0,1 et 1 Hz, - le procédé comprend le positionnement du géophone in situ avant la mise en oeuvre de la détermination de l'amortissement intrinsèque et de la fréquence propre intrinsèque. Le procédé de filtrage est donc adapté aux 15 conditions d'utilisation réelles du géophone. Cela est particulièrement avantageux lorsque le géophone est soumis à des températures et des inclinaisons différentes des conditions d'utilisation nominales du géophone, par exemple lorsque le géophone est placé dans un puits profond, - le signal de réponse du géophone mesuré à la disparition du 20 courant d'excitation comporte au moins deux extrema locaux consécutifs et une valeur prédéterminée correspondant au signal de réponse mesuré en l'absence de courant d'excitation et de vibrations, la valeur prédéterminée étant comprise entre les deux extrema locaux, et la fréquence propre intrinsèque et l'amortissement intrinsèque sont déterminés à partir des formules suivantes : 25 m= 1 , et tan(e) = At2V1 - m2' dans lesquelles: At1 est un premier intervalle de temps mesuré entre la disparition du courant d'excitation et le premier extremum local, 30 At2 est un second intervalle de temps mesuré entre la disparition du courant d'excitation et la première valeur nulle prédéterminée, et 4 wo est la pulsation propre définie par wo=2-rrfo, fo étant la fréquence propre intrinsèque, - l'amortissement intrinsèque est inférieur à 1, de préférence inférieur à 0,65, 5 - le gain du géophone est également ajusté de telle sorte que ledit gain soit également sensiblement constant pour des fréquences supérieures à la fréquence propre intrinsèque. On obtient ainsi un gain constant pour des fréquences supérieures à la fréquence propre de consigne, ce qui permet de mesurer des vibrations sur une plus grande gamme de fréquences, 10 - le filtre numérique présente une fonction de transfert correspondant à l'inverse de la fonction de transfert du géophone pour une fréquence supérieure à la fréquence propre de consigne, - le filtre présente la fonction de transfert suivante : p2+2mwo+wo2 TFfilter(P) p2+2m1w1+w121 15 dans laquelle w0 est la pulsation propre intrinsèque définie par wo=2-rdo, fo étant la fréquence propre intrinsèque et w1 est la pulsation propre de consigne définie par w1=2-rrr1, f1 étant la fréquence propre de consigne et m1 l'amortissement de consigne. - le filtre est un filtre à réponse impulsionnelle infinie de la forme : 20 Yn = kxnxn + + - - avec xn, xn_i et xn_2 correspondant au signal de réponse du géophone à l'instant t et aux instants t-T et t-2T, respectivement, et yn, yn_i et yn-2 correspondant au signal de réponse du géophone modifié par le filtre à l'instant t et aux instants t-T et t-2T respectivement, 25 T étant la période d'échantillonnage et kxn_2, kyn-2 étant les coefficients du filtre définis par 5 kxnK ( ki \ = 2T, + k2) kxn-1= 1 8 (2k2 kxn-2 = 1714 2 + k2) 1 8 kyn-1 = 7(21(2 -7.4 1 4 K1 ( kyn-2 = 2 T, + K2), avec K1 = 2M",01, K2 = w12, k1 = 217/G00, k2 = G002, et K = ( 4 + 2Ki + K2), T les coefficients kxn, kxn-1 et kxn_2 étant en outre déterminés à l'aide des valeurs d'amortissement intrinsèque et de fréquence propre intrinsèque 5 déterminées lors de la mise en oeuvre de la détermination de l'amortissement intrinsèque et de la fréquence propre intrinsèque, l'utilisation d'un filtre à réponse impulsionnelle infinie permettant d'utiliser des composants avec des ressources de calcul et de mémoire limités, résultant en un capteur moins coûteux et plus autonome en énergie, 10 - l'amortissement de consigne ml est égal à V2/2. La présente invention propose, selon un second aspect, un module de filtrage du signal de réponse d'un géophone, ledit géophone comportant un aimant entouré d'une bobine montée mobile par rapport à l'aimant, la bobine 15 étant apte à se déplacer suite à des vibrations d'un environnement de ladite bobine, ledit module comprenant : un calculateur configuré pour : o déterminer l'amortissement intrinsèque et la fréquence propre intrinsèque du géophone en alimentant la bobine avec un courant 20 d'excitation pendant une durée prédéterminée et en déterminant, à l'aide du signal de réponse du géophone mesuré à la disparition du courant d'excitation, la fréquence propre intrinsèque et l'amortissement 6 intrinsèque du géophone, o déterminer un filtre numérique configuré pour modifier une fréquence propre du géophone d'une fréquence propre intrinsèque à une fréquence propre de consigne, la fréquence propre de consigne 5 étant inférieure à la fréquence propre intrinsèque et ajuster le gain du géophone de telle sorte que ledit gain soit sensiblement constant pour des fréquences comprises entre la fréquence propre de consigne et la fréquence propre intrinsèque, o appliquer ledit filtre sur le signal de réponse fourni par le 10 géophone. La présente invention propose, selon un troisième aspect, un capteur vélocimétrique configuré pour la mesure de vibrations comportant : - un géophone comportant un aimant entouré d'une bobine montée 15 mobile par rapport à l'aimant, - un module tel que décrit ci-dessus, - un générateur de courant adapté pour recevoir un signal de commande du module et pour générer en réponse audit signal de commande un courant d'excitation de durée déterminée dans la bobine du géophone. 20 Enfin, la présente invention propose, selon un dernier aspect, un programme d'ordinateur comprenant des instructions permettant la mise en oeuvre du procédé tel que décrit précédemment. 25 Brève description des dessins Des détails et avantages de la présente invention apparaitront mieux de la description qui suit, faite en référence aux dessins annexés sur lesquels : La figure 1 représente un module de filtrage du signal de réponse d'un géophone ainsi qu'un capteur vélocimétrique équipé d'un tel module, 30 La figure 2 représente le procédé de filtrage d'un signal de réponse d'un géophone selon un mode de réalisation de l'invention, La figure 3 représente le signal de réponse du géophone suite à l'application d'un courant d'excitation dans la bobine du géophone pendant une 7 durée déterminée, La figure 4 représente un exemple de fonction de transfert du géophone, de fonction de transfert du filtre numérique et de fonction de transfert du géophone après application du filtre numérique. 5 Description détaillée d'un mode de réalisation de l'invention La figure 1 représente de manière schématique un capteur vélocimétrique 1 comprenant un module de filtrage 3 selon un mode de réalisation de l'invention. 10 Le capteur 1 comprend un géophone 2, un module de filtrage 3 et un générateur de courant 4. Le générateur de courant 4 peut être une carte externe reliant le module 3 au géophone 2 comme représenté dans la figure 1 ou être intégré dans le module 3. Le géophone 2 comprend un aimant entouré d'une bobine montée 15 mobile par rapport à l'aimant. Il peut s'agir, par exemple, d'un géophone SM-6 de la marque ION ayant une fréquence propre intrinsèque de 4,5 Hz (+1-10%). En présence de vibrations, la bobine bouge par rapport à l'aimant, ce qui produit un courant induit. Le géophone 2 peut comprendre également une résistance branchée en court-circuit aux bornes de la bobine. Cette résistance 20 permet de définir la valeur de l'amortissement intrinsèque m. Le signal de réponse du géophone S geoph est représentatif du courant induit dans la bobine par le déplacement de la bobine par rapport à l'aimant. Il peut s'agir, par exemple, d'une tension. Le module de filtrage 3 comprend un calculateur 32 et une mémoire 25 31. Le calculateur 32, ici un microprocesseur, est configuré pour déterminer l'amortissement intrinsèque m et la fréquence propre intrinsèque f0 du géophone et pour déterminer et appliquer un filtre numérique. De manière avantageuse, le filtre numérique est configuré pour modifier une fréquence propre du géophone d'une fréquence propre intrinsèque au géophone utilisé, à 30 une fréquence propre de consigne inférieure à la fréquence propre intrinsèque et ajuster le gain du géophone de telle sorte que ledit gain soit sensiblement constant au moins pour des fréquences comprises entre la fréquence propre de consigne et la fréquence propre intrinsèque. On s'assure ainsi que le gain 8 est constant pour les signaux de courte période qu'on souhaite mesurer. Ces signaux ont par exemple une fréquence comprise entre 0,1 et 2 Hz. La fréquence propre de consigne est donc choisie en fonction de l'application visée. L'expression et les coefficients du filtre numérique sont stockés dans la 5 mémoire 31 ainsi que les instructions de programme d'ordinateur permettant notamment de déterminer l'amortissement intrinsèque m et la fréquence propre intrinsèque f0 du géophone, d'actualiser les coefficients du filtre numérique en fonction des valeurs d'amortissement intrinsèque m et la fréquence propre intrinsèque f0 déterminés et d'appliquer le filtre numérique sur le signal de 10 réponse du géophone Sgeoph, et plus particulièrement sur le signal de réponse échantillonné du géophone Sgeoph ech. Le générateur de courant 3 est adapté pour recevoir un signal de commande Sem du module et pour générer en réponse audit signal de commande un courant d'excitation 'ex, de durée déterminée dans la bobine du 15 géophone 2. Le signal de réponse Sgeoph du géophone 2 mesuré à la disparition du courant d'excitation est utilisé pour déterminer l'amortissement intrinsèque m et la fréquence propre intrinsèque fo du géophone 2. Dans ce mode de réalisation, le module 3 comprend un convertisseur analogique numérique 33 configuré pour transformer le signal de réponse du 20 géophone Sgeoph en un signal de réponse du géophone échantillonné Sgéoph,ech. Le signal de réponse échantillonné Sgéoph,ech est transmis au microprocesseur 32 qui l'utilise dans certaines conditions, soit pour déterminer l'amortissement intrinsèque m et la fréquence propre intrinsèque f0, soit pour y appliquer le filtre numérique. Une fois le filtre numérique appliqué, le module 3 et donc le capteur 25 1 restitue le signal Scapteur. Un module 3 adapté pour la mise en oeuvre de l'invention et comprenant un convertisseur analogique est par exemple un module D9-DIN commercialisé par la société STANEO. La figure 2 représente le procédé de filtrage du signal de réponse d'un 30 géophone mis en oeuvre dans le capteur décrit précédemment selon un mode de réalisation de l'invention. Le procédé comprend une étape S100 d'initialisation dans laquelle l'expression et les coefficients du filtre numérique sont chargés dans le module 9 de filtrage 3. On décrit ici la façon dont est obtenue l'expression du filtre et des coefficients du filtre numérique. De manière connue de l'homme du métier, un géophone comportant une bobine mobile autour d'un aimant fixe se comporte 5 comme un filtre passe-haut d'ordre 2 de fréquence fo et d'amortissement m. La fonction de transfert intrinsèque du géophone en réponse à la vitesse de vibration d'un géophone est donc de la forme : p2 avec (.00 = 2iffo.et p = 27-cf, f étant une TFgi(p) = p2+2mwo+w02' fréquence. 10 La fonction de transfert du géophone TFg; pour un géophone avec un amortissement d'environ 0,55 est représentée en référence à la figure 4. La fonction de transfert intrinsèque du géophone TFg; présente un gain normalisé sensiblement constant égal à 1 pour des fréquences supérieures à la fréquence propre intrinsèque fo et une pente constante pour des fréquences 15 inférieures à la fréquence propre intrinsèque, de valeur 1 comme illustré par la courbe en pointillés. Au sens de l'invention, la fréquence propre du géophone est définie comme étant l'intersection entre la portion de pente constante non nulle et la portion de gain constant. Si l'on considère la fonction de transfert de réponse à l'accélération des vibrations du géophone, la pente a une valeur de 20 2. On notera que le gain normalisé de la fonction de transfert intrinsèque du géophone autour de la fréquence propre intrinsèque fo est supérieur à 1. Le géophone se trouve donc en résonance pour des fréquences proches de la fréquence propre intrinsèque fo. On notera que, pour des valeurs d'amortissement plus élevées, de l'ordre de V/2 par exemple, le géophone 25 n'entre plus en résonance. Le gain normalisé du géophone est alors égal à 1 pour des fréquences supérieures à la fréquence propre intrinsèque fo et le gain présente une pente constante en dessous de la fréquence propre intrinsèque fo. Le filtre numérique selon l'invention permet d'obtenir une fonction de 30 transfert du géophone TFgf, comme celle représentée sur la figure 4 par la courbe en trait mixte, avec un gain sensiblement constant pour des fréquences comprises entre la fréquence propre de consigne f1 et la fréquence propre intrinsèque fo. 10 Pour obtenir une telle fonction de transfert du géophone après filtrage TFgf, il convient d'appliquer un filtre avec une fonction de transfert TFf;lter, telle que celle représentée sur la figure 4 par la courbe en trait plein par exemple, correspondant à l'inverse de la fonction de transfert du géophone TFg; pour des 5 fréquences comprises entre la fréquence propre de consigne f1 et la fréquence propre intrinsèque fo. Un tel filtre peut avoir une fonction de transfert TFfilter de la forme : p2 + 2mwo + (i)02 TFfilter (P) = p2 + 2m1w1 + co12 pour une fonction de transfert de réponse à la vitesse de vibration d'un géophone avec col = 2icf1i et m1 étant l'amortissement de consigne, par 10 exemple égal à V/2. La fonction de transfert TFfilter décrite ici et représentée par la courbe en trait plein correspond à l'inverse de la fonction de transfert du géophone TFg; pour une fréquence supérieure à la fréquence propre de consigne Un filtre avec une telle fonction de transfert TFfilter permet également de 15 faire disparaître la résonance autour de la fréquence propre fo, si une telle résonance est présente, et d'obtenir ainsi un gain normalisé sensiblement égal à 1 pour des fréquences comprises entre la fréquence propre de consigne et la fréquence propre intrinsèque, mais également pour des fréquences supérieures à la fréquence propre intrinsèque fo. 20 On calcule ensuite l'expression et les coefficients du filtre numérique présentant la fonction de transfert TFfilter précédemment décrite. On calcule, par exemple, l'expression et les coefficients d'un filtre récursif de type à réponse impulsionnelle infinie RII (« Infinite Response Filter » en anglais). Pour cela, on effectue une transformée bilinéaire sur la fonction de 25 transfert TFfilter du filtre déterminée précédemment en remplaçant le paramètre 2 1+zz- - 1 p par T -, , T étant la période d'échantillonnage du signal dépendant du 1 convertisseur analogique numérique utilisé. La fonction de transfert devient alors : Z-2 (z72-21,÷1+k)+z-1(74+2k2)+C*+21,÷1+k2) TFfilter(Z) = 2(T2 2K +10-EZ-1(7-8+2K2)+(74+7+/(2) 11 K1 = 2M",01, K2 = (i)12, avec k1 = 2mwo, k2 = (.002, m1 étant l'amortissement de consigne, de préférence égal à V2/2 et (1)1 = 2fff1, avec f1 étant la fréquence de consigne. En écrivant l'équation aux différences correspondante, on obtient 5 l'expression du filtre numérique à appliquer : Yn = kxnxn + + k'n_2Xn_2 - kyn_iyn_i - kyn_2yn_2, (1) avec les coefficients suivants : kxn = 1( 1 (7,42 + 2 kT, k2) kxn-1 4 (2k2 - 778 ) kxn_2 ( k, kyn-1 kyn-2 ( 4 2 7, k2) , (2K2 - 778 ) + 2 + K2) 27,1(1 + K2), la variable K étant définie par la formule suivante : K = G42 + 10 Xn, Xn-i et Xn-2 correspondant aux signaux de réponse du géophone Sgeoph,ech fournis par un convertisseur analogique numérique à l'instant t courant et aux instants t-T et t-2T respectivement, T étant la période d'échantillonnage, et Yn, Yn-1 et Yn-2 correspondant aux signaux de réponse du géophone Scapteur 15 obtenus après application du filtre numérique décrit dans l'équation (1) à l'instant t et aux instants t-T et t-2T respectivement. Le procédé comprend également une étape S110 de détermination de l'amortissement intrinsèque m et de la fréquence propre intrinsèque f0 du géophone. L'étape S110 comprend une sous-étape S111 d'alimentation de la 20 bobine avec un courant d'excitation lex, pendant une durée prédéterminée At et une sous-étape S112 de détermination, à l'aide du signal de réponse du géophone mesuré Sgeoph,ech à la disparition du courant d'excitation lexe, de la fréquence propre intrinsèque fo et de l'amortissement intrinsèque m du géophone. 25 La figure 3 représente de manière schématique le courant d'excitation lexc alimentant la bobine pendant une durée déterminée At. Il s'agit d'un 12 créneau de courant de durée At. Le signal de réponse du géophone Sgeoph présente une première série d'oscillations OSC1 apparaissant après le front montant du courant d'excitation lexe, se stabilise puis présente une deuxième série d'oscillations OSC2 apparaissant après le front descendant, soit à la 5 disparition du courant d'excitation, puis se stabilise. Seule la deuxième série d'oscillations OSC2 est prise en compte pour la détermination de l'amortissement m et de la fréquence propre fo intrinsèque au géophone car elle correspond à une oscillation libre de la bobine. En effet, suite à l'application du courant d'excitation lex, dans la bobine, lors du front montant, le champ 10 magnétique induit créé par la bobine est confronté au champ magnétique de l'aimant. La bobine se met donc à osciller par rapport à l'aimant jusqu'à atteindre une position d'équilibre, produisant ainsi la première série d'oscillations OSC1 suivie de la première stabilisation du signal de réponse du géophone Sgeoph. Ce signal de réponse mesuré pendant l'intervalle de temps At 15 n'est donc pas exploitable pour évaluer l'amortissement m et la fréquence propre fo intrinsèque du géophone en réponse à un signal d'excitation mécanique, puisque le signal de réponse du géophone est perturbé par la présence du courant d'excitation. En effet, dans les conditions réelles d'utilisation du géophone, la bobine n'est pas alimentée en courant mais 20 produit un courant induit par le déplacement de la bobine de l'aimant en réponse aux vibrations de l'environnement. On se différencie ainsi de l'art antérieur en déterminant la fréquence propre et l'amortissement à la disparition du courant d'excitation. On notera que le signal de réponse du géophone, suite à la disparition 25 du courant d'excitation lex, positif, c'est-à-dire après le front descendant du créneau, présente une seconde série d'oscillations comprenant au moins un premier minimum local à l'instant t1, une première valeur nulle à l'instant t2 un premier maximum local à l'instant t3 puis se stabilise. La valeur nulle à l'instant t2 correspond à la valeur du signal de réponse du géophone mesuré 30 auparavant en l'absence de courant d'excitation et de mouvement de la bobine. Il s'agit donc d'un signal de réponse de repos du géophone. On notera que dans le cas d'un courant d'excitation négatif, la seconde série d'oscillations présente un maximum local suivi d'un minimum local. 13 La présence de deux extrema locaux, autrement dit de deux oscillations, avant que le signal de réponse se stabilise à la valeur du signal de réponse de repos, est caractéristique d'un amortissement m, inférieur à 1 appelé amortissement sous critique. Par ailleurs, lorsque l'amortissement m est 5 très inférieur à 1 (typiquement lorsqu'il est inférieur à 0,6), d'autres extrema locaux ou oscillations peuvent devenir significatifs, révélateurs d'un amortissement moindre et donc d'oscillations présentes plus longtemps. Il a été démontré par les inventeurs, qu'en déterminant l'intervalle de temps At1 = t1 - tO et l'intervalle de temps At2 = t2 - tO, il était possible de 10 déterminer l'amortissement m et la fréquence propre intrinsèque du géophone à l'aide des formules suivantes : 1 m - +tanOre VVV IT ° = Llt21 /1-M2' w0 étant la pulsation propre intrinsèque définie par wo=2-rdo, fo étant la 15 fréquence propre intrinsèque. Une telle méthode de détermination de l'amortissement intrinsèque m et de la fréquence propre intrinsèque f0 est particulièrement avantageuse en ce qu'elle ne se base pas sur la mesure de l'instant t3 qui peut être difficile à mesurer si ce maximum local est très étalé dans le temps, notamment lorsque 20 l'amortissement est supérieur à 0,65. Il est toutefois nécessaire d'avoir au moins deux extrema locaux avant d'atteindre la valeur de repos du signal de réponse pour pouvoir déterminer At1 et At2. Dans l'application visée ici, on considère donc des géophones avec des amortissements intrinsèques inférieurs à 1, par exemple compris entre 0,2 et 0,65. A titre d'exemple, pour un 25 amortissement intrinsèque m de l'ordre de 0,2, le signal de réponse du géophone est stabilisé après un intervalle de temps inférieur à une seconde, ce qui permet une détermination de la fréquence propre intrinsèque fo et de l'amortissement intrinsèque m assez rapide. Pour une valeur d'amortissement intrinsèque de l'ordre de 0,65, le deuxième extremum local est significatif et 30 permet de mesurer l'instant t2 avec une bonne précision. Les intervalles de temps At1 et At2 sont déterminés à l'aide d'algorithmes de traitement du signal appropriés. La conception de tels algorithmes est à la portée de l'homme du 14 métier. Le procédé comprend ensuite une étape S120 de détermination du filtre numérique correspondant à l'actualisation des coefficients du filtre numérique définis dans l'étape d'initialisation S100 en fonction des valeurs 5 d'amortissement m et de fréquence propre intrinsèque fo déterminées dans l'étape S110. Enfin, le procédé comprend une étape S130 d'application du filtre numérique sur le signal de réponse du géophone Sgeoph et plus particulièrement sur le signal de réponse du géophone échantillonné Sgeoph,ech. 10 L'application d'un filtre à réponse impulsionnelle infinie présente l'avantage de consommer moins de ressources de calcul et de mémoire qu'un filtre à réponse impulsionnelle finie. Il est donc possible d'implémenter le filtre avec des composants moins couteux, par exemple avec un microcontrôleur présentant une fréquence de 400 MHz pour le coeur et de 100 MHz pour le bus 15 externe sur lequel le signal de réponse du géophone échantillonné Sgeoph,ech et le signal du capteur Scapteur sont acheminés. Bien entendu, il est également possible d'utiliser un filtre numérique à réponse impulsionnelle finie ou tout autre filtre numérique, le choix du type de filtre se faisant en fonction de la puissance de calcul et de la mémoire des équipements disponibles. 20 Selon une option, le géophone peut être positionné in situ avant la mise en oeuvre de l'étape S110 de détermination de l'amortissement intrinsèque m et de la fréquence propre intrinsèque fo. Il est ainsi possible de mesurer l'amortissement intrinsèque m et la fréquence propre intrinsèque fo dans les conditions réelles d'utilisation du géophone et donc d'obtenir un filtre 25 adapté aux conditions réelles d'utilisation, la température et l'inclinaison du géophone ayant une influence sur les caractéristiques intrinsèques du géophone, à savoir son amortissement m et sa fréquence propre fo. Les étapes S110 et S120 peuvent être réitérées de manière périodique ou non, afin de tenir compte de possibles changements dans l'environnement 30 du géophone ayant une influence sur l'amortissement REVENDICATIONS1. Procédé de filtrage d'un signal de réponse d'un géophone, ledit géophone comportant un aimant entouré d'une bobine montée mobile par rapport à l'aimant, la bobine étant apte à se déplacer suite à des vibrations d'un environnement de ladite bobine, le géophone présentant un amortissement intrinsèque (m) et une fréquence propre intrinsèque (fo), ledit procédé comportant: - la détermination de l'amortissement intrinsèque (m) et de la 10 fréquence propre intrinsèque (fo) du géophone (S110), comprenant : o l'alimentation de la bobine avec un courant d'excitation pendant une durée prédéterminée (S111), o la détermination, à l'aide du signal de réponse du géophone mesuré à la disparition du courant d'excitation, de la 15 fréquence propre intrinsèque et de l'amortissement intrinsèque du géophone (S112), - la détermination d'un filtre numérique configuré pour modifier une fréquence propre du géophone d'une fréquence propre intrinsèque (fo) à une fréquence propre de consigne (f1), la fréquence propre de consigne (f1) étant 20 inférieure à la fréquence propre intrinsèque (fo) et ajuster le gain du géophone de telle sorte que ledit gain soit sensiblement constant pour des fréquences comprises entre la fréquence propre de consigne (f1) et la fréquence propre intrinsèque (S120), - l'application dudit filtre sur ledit signal de réponse du géophone 25 (S130). 2. Procédé selon la revendication 1, caractérisé en ce que la fréquence propre de consigne est comprise entre 0,1 et 1 Hz. 30 3. Procédé selon l'une des revendications précédentes, caractérisé en ce que le procédé comprend le positionnement du géophone in situ avant la mise en oeuvre de la détermination de l'amortissement intrinsèque et de la fréquence propre intrinsèque. 2 4. Procédé selon l'une des revendications précédentes, caractérisé en ce que le signal de réponse du géophone mesuré à la disparition du courant d'excitation comporte au moins deux extrema locaux consécutifs et une valeur 5 prédéterminée correspondant au signal de réponse mesuré en l'absence de courant d'excitation et de vibrations, la valeur prédéterminée étant comprise entre les deux extrema locaux, et en ce que la fréquence propre intrinsèque et l'amortissement intrinsèque sont déterminés à partir des formules suivantes : 1 m= , et ,ii+tan(e) TC coo = At-2,11 - m2 dans lesquelles : At1 est un premier intervalle de temps mesuré entre la disparition du courant d'excitation et le premier extremum local, At2 est un second intervalle de temps mesuré entre la disparition du courant d'excitation et la première valeur nulle prédéterminée, et wo est la pulsation propre définie par wo=2-rrfo, fo étant la fréquence propre intrinsèque. 5. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'amortissement intrinsèque est inférieur à 1, de préférence inférieur à 20 0,65. 6. Procédé selon l'une des revendications précédentes, caractérisé en ce que le gain du géophone est également ajusté de telle sorte que ledit gain soit également sensiblement constant pour des fréquences supérieures à la 25 fréquence propre intrinsèque (fo). 7. Procédé selon la revendication précédente, caractérisé en ce que le filtre numérique présente une fonction de transfert (TFfilter) correspondant à l'inverse de la fonction de transfert du géophone (TF9,) pour une fréquence 30 supérieure à la fréquence propre de consigne (f1). 3 8. Procédé selon la revendication 7, caractérisé en ce que le filtre présente la fonction de transfert suivante : p2+2-mwo+wo2 TFfiiter(P) - p2+2m1w1l-co12' dans laquelle w0 est la pulsation propre intrinsèque définie par 5 wo=2-1-rfo, fo étant la fréquence propre intrinsèque et w1 est la pulsation propre de consigne définie par w1=2Trf1, f1 étant la fréquence propre de consigne et m1 l'amortissement de consigne. 9. Procédé selon l'une des revendications 7 à 8, caractérisé en ce que 10 le filtre numérique est un filtre à réponse impulsionnelle infinie de la forme : Yn = Ic''x7, + + k''_2x,i_2 - kyn_iy'_i - kyn_2317,_2, avec xn, xn_i et xn_2 correspondant au signal de réponse du géophone à l'instant t et aux instants t-T et t-2T, respectivement, et yn, Yn-1 et yn-2 correspondant au signal de réponse du géophone modifié 15 par le filtre numérique à l'instant t et aux instants t-T et t-2T respectivement, T étant la période d'échantillonnage et kxn, kxn-1, kxn-2, kyn-1, kyn-2 étant les coefficients du filtre numérique définis par ( ki kx' + 2 7 + k2) , 1 8 kxn-1 = (2k2), ( 4 ki kxn-2 - 2 T, + k2) , 1 8 kyn_i = --K- (2K2 - T2)' 1 4 Ki kyn_2 = K (7 + 2 T + K2) avec K1 = 2/111(01, K2 = (012, k1 = 21T/L00, k2 = W02, et K = (T2 ±2K1 ± K2), les coefficients kxn, kxn-1 et kxn-2 étant en outre déterminés à l'aide des valeurs d'amortissement intrinsèque (m) et de fréquence propre intrinsèque (fo) 4 déterminées lors de la mise en oeuvre de la détermination de l'amortissement intrinsèque (m) et de la fréquence propre intrinsèque (fo) (S110). 10. Procédé selon la revendication 9, caractérisé en ce que 5 l'amortissement de consigne ml est égal à N/2/2. 11. Module de filtrage (3) du signal de réponse Sgeoph d'un géophone (2), ledit géophone (2) comportant un aimant entouré d'une bobine montée mobile par rapport à l'aimant, la bobine étant apte à se déplacer suite à des 10 vibrations d'un environnement de ladite bobine, ledit module (3) comprenant : - un calculateur (32) configuré pour : o déterminer l'amortissement intrinsèque (m) et la fréquence propre intrinsèque (f0) du géophone en alimentant la bobine avec un courant d'excitation pendant une durée prédéterminée (At) et en 15 déterminant, à l'aide du signal de réponse du géophone (Sgeoph) mesuré à la disparition du courant d'excitation (OSC2), la fréquence propre intrinsèque (fo) et l'amortissement intrinsèque (m) du géophone (2), o déterminer un filtre numérique configuré pour modifier une fréquence propre du géophone d'une fréquence propre intrinsèque (fo) à 20 une fréquence propre de consigne (f1), la fréquence propre de consigne (fi) étant inférieure à la fréquence propre intrinsèque (fo) et ajuster le gain du géophone de telle sorte que ledit gain soit sensiblement constant pour des fréquences comprises entre la fréquence propre de consigne (f1) et la fréquence propre intrinsèque (fo), o appliquer ledit filtre sur le signal de réponse Sgeoph fourni par le géophone (2). 12. Capteur vélocimétrique (1) configuré pour la mesure de vibrations comportant : - un géophone (2) comportant un aimant entouré d'une bobine montée mobile par rapport à l'aimant, - un module (3) selon la revendication 11, - un générateur de courant (4) adapté pour recevoir un signal de 5 commande (Scom) du module (3) et pour générer en réponse audit signal de commande (Scom) un courant d'excitation (lexc) de durée déterminée (At) dans la bobine du géophone (2). 5 13. Produit programme d'ordinateur comprenant des instructions permettant la mise en oeuvre du procédé selon la revendication 1 lorsque ledit programme est exécuté sur un ordinateur.