FR 3058477 A1 20180511 FR 1660722 A 20161107 1 Domaine de l'invention Le domaine de la présente invention concerne celui des procédés de détection de la défaillance d'un élément mobile d'un aérogénérateur. Art antérieur Depuis plusieurs décennies, des solutions sont proposées pour améliorer la détection de défaillance des éléments constitutifs d'une ligne d'arbre d'une machine tournante. Certaines techniques d'analyses, bien connues de l'homme du métier, comme par exemple l'analyse du lubrifiant ou bien le suivi de température, semblent limitées et ne détectent le défaut qu'une fois survenu, comme c'est le cas du suivi de température, ou sont difficiles à mettre en oeuvre, comme dans le cas de l'analyse de lubrifiant. D'autres techniques mettent en oeuvre une analyse vibratoire. A titre d'exemple, les techniques définissant des indicateurs globaux du type « rms » ou « kurtosis » ne permettent que de détecter les défaillances graves. Il est également commun d'utiliser des indicateurs vectoriels pour suivre l'évolution du signal. En particulier, l'observation du signal vibratoire revient à calculer un spectre ou éventuellement un spectre d'enveloppe, puis à mesurer l'amplitude des canaux fréquentiels de chacune des fréquences qui constituent la famille d'harmonique caractéristique du mode de défaillance sous surveillance. A titre d'exemple, suivre un défaut de bague externe revient à suivre les 3 premières harmoniques de la fréquence caractéristique de bague externe du roulement sur le spectre. Néanmoins, dans le cas d'un défaut faiblement énergétique comme par exemple, un roulement lent, il faut augmenter la résolution fréquentielle pour espérer différentier l'amplitude du canal fréquentiel du bruit ambiant. Ceci impose de connaître la fréquence caractéristique avec une précision suffisante, car l'augmentation de la résolution fréquentielle peut faire apparaître des composantes dites fantômes en plus de la composante recherchée. Toutefois, cette précision de connaissance de la fréquence caractéristique ne peut être atteinte, car la fréquence caractéristique présente une variabilité. En effet, cette fréquence caractéristique dépend de l'angle de contact dans le roulement qui lui-même dépend des conditions de fonctionnement appliqué au roulement. 3058477 2 Ainsi, il n'existe pas de procédé fournissant un indicateur fiable qui puisse permettre l'évaluation et la surveillance d'un état de roulement à partir d'une analyse de signal vibratoire. En particulier, il n'existe par de procédé de détection efficace pour les 5 lignes d'arbre lentes et des défauts peu énergétiques. Dans le cas des défaillances de roulement principal d'éolienne, les roulements tournent à faible vitesse, c'est-à-dire une vitesse de l'ordre de quelques dizaines de tours par minutes, par exemple 20 tours par minutes. Notamment, les procédé existants ne sont pas adaptés à la détection de défaut de cage, d'éléments roulants, de bague externe ou de bague 10 interne de roulement, et de façon plus générale aux défauts se manifestant sous forme de famille d'harmoniques dans un spectre. Exposé de l'invention La présente invention a pour but de résoudre tout ou partie des inconvénients mentionnés ci-dessus sous la forme d'un procédé de détection de la 15 défaillance d'un élément mobile d'un aérogénérateur, le procédé de détection comprend les étapes suivantes : mesure d'une grandeur physique représentative du mouvement de l'élément mobile de l'aérogénérateur de sorte à générer un signal mesuré temporel ; 20 - transformation du signal mesuré temporel ou d'un signal temporel obtenu à partir du signal mesuré temporel de sorte à générer un signal fréquentiel présentant un premier spectre ; calcul d'un spectre transformé dont la valeur pour une fréquence donnée correspond au produit de la valeur du premier spectre ou 25 d'un spectre obtenu à partir du premier spectre et d'un ensemble de valeurs du premier spectre ou du spectre obtenu à partir du premier spectre correspondant à des multiples de la fréquence donnée. vérification de la réalisation d'un critère de détection pour des 30 valeurs du spectre transformé dans un intervalle fréquentiel à proximité d'une valeur caractéristique d'un défaut. 3058477 3 Ainsi, grâce à cette disposition, il est possible de détecter l'état de l'élément mobile de l'aérogénérateur lorsque l'élément mobile de l'aérogénérateur est défaillant, la défaillance se manifestant par une vibration de faible énergie. Selon un aspect de l'invention, la valeur caractéristique d'un défaut 5 correspond à une fréquence caractéristique d'un défaut à détecter. Selon un aspect de l'invention, l'étape de calcul d'un spectre transformé prend en compte un nombre de multiples de la fréquence donnée correspondant à un nombre d'harmoniques d'une fréquence caractéristique du défaut à détecter. Selon un mode de réalisation, le procédé de détection comprend une 10 étape de décalage fréquentiel résultant en un spectre dont la valeur pour une fréquence donnée correspond au produit de la valeur du premier spectre ou d'un spectre obtenu à partir du premier spectre et d'un ensemble de valeurs du premier spectre ou d'un spectre obtenu à partir du premier spectre correspondant à la somme de la fréquence donnée et d'un multiple de la fréquence caractéristique de l'élément 15 mobile de l'aérogénérateur. Selon un aspect de l'invention, la fréquence caractéristique de l'élément mobile est la fréquence de rotation de l'arbre sur lequel la bague est fixée dans le cas de la détection d'une défaillance d'une bague interne tournante. Selon un mode de réalisation, le signal fréquentiel est généré en 20 appliquant au signal mesuré temporel une transformée de fourrier. Selon un mode de réalisation, le procédé de détection comprend une étape d'amélioration de la qualité du signal mesuré temporel de façon à augmenter la résolution fréquentielle par un facteur au moins égal au nombre de multiples. Selon un mode de réalisation, l'étape d'amélioration comprend une 25 étape d'interpolation du spectre Selon un mode de réalisation, l'étape d'interpolation est réalisée par ajout de valeurs correspondant à une valeur nulle ou à une valeur moyenne à l'extrémité du signal mesuré temporel pour l'ensemble de valeurs de façon à multiplier la longueur du signal mesuré temporel par le facteur au moins égal au 30 nombre de multiples. 3058477 4 Selon un mode de réalisation, le procédé comprend une étape de nettoyage du spectre dans laquelle au moins une valeur du spectre correspondant à au moins une fréquence parasite est remplacée par une valeur de lissage. Selon un mode de réalisation, la valeur de lissage correspond à la valeur 5 minimale du spectre ou à une valeur nulle. Ainsi, grâce à cette disposition, il est possible de conserver seulement les harmoniques nécessaires à la détection de la défaillance de l'état de l'élément mobile de l'aérogénérateur. Selon un aspect de l'invention, l'étape de nettoyage du spectre est 10 réalisée préalablement à l'étape de décalage fréquentiel. D'autres caractéristiques et avantages de l'invention apparaîtront mieux à la lecture de la description qui suit d'un mode de réalisation de l'invention donné à titre d'exemple non limitatif. Liste des figures 15 L'invention sera mieux comprise à l'aide de la description détaillée qui est exposée ci-dessous en regard du dessin, dans lesquels : la figure 1 représente un exemple de mesure d'une grandeur physique représentative du mouvement de l'élément mobile de l'aérogénérateur ; 20 - la figure 2 montre une étape d'amélioration de la qualité du signal mesuré temporel ; - la figure 3 illustre un exemple de transformation du signal mesuré temporel ; - la figure 4 expose un exemple de spectre nettoyé ; 25 la figure 5 présente un exemple de décalage fréquentiel selon l'invention correspondant à moins 2f1 ; la figure 6 représente un exemple de décalage fréquentiel selon l'invention correspondant à moins f1 ; 3058477 5 la figure 7 illustre un exemple de décalage fréquentiel selon l'invention correspondant à plus f1 la figure 8 illustre un exemple de décalage fréquentiel selon l'invention correspondant à plus 2f1 5 la figure 9 expose un exemple de produit de la valeur du premier spectre et d'un ensemble de valeurs du premier spectre ; - la figure 10 présente un exemple de transformation par produit d'harmonique ; et 10 - la figure 11 présente un système mettant en oeuvre le procédé selon l'invention. Dans la description détaillée qui va suivre des figures définies ci-dessus, les mêmes éléments ou les éléments remplissant des fonctions identiques pourront conserver les mêmes références de manière à simplifier la compréhension de 15 l'invention. Description d'un mode réalisation de l'invention La figure 11 représente un système équipant un parc éolien et comprenant un ensemble d'aérogénérateurs 2. Un aérogénérateur comprend un ensemble électromécanique 3 20 comprenant notamment un élément mobile en rotation 4. L'ensemble électromécanique est équipé d'au moins un capteur 5 relié à un circuit d'acquisition 6. Le circuit d'acquisition produit un signal temporel mesuré xl. Les signaux temporels mesurés peuvent être collectés au niveau d'un poste de livraison 7 comprenant par exemple un serveur de calcul 8 muni d'au moins un 25 processeur programmé pour transformer le signal et déterminer des indicateurs I, dont la nature sera décrite ultérieurement. Ces indicateurs I peuvent être communiqués à un serveur distant, notamment à un serveur de données accessible par Internet en vue d'un stockage de 3058477 6 l'intégralité des indicateurs, afin par exemple de mettre à disposition ces indicateurs pour des interrogations par les exploitants. Ces indicateurs I sont également transmis à un centre de supervision 10 comprenant également un ou plusieurs terminaux 11 munis d'au moins un processeur 5 agencé pour communiquer des informations d'historique à un utilisateur, mais également pour permettre la surveillance et le réglage de seuils ou de critères de détection, et pour déclencher des alarmes. Le signal x1, représenté à la figure 1, est obtenu à partir d'un capteur 5 et du circuit d'acquisition 6 en fonction du temps T. Le signal x1 peut être notamment un 10 signal échantillonné obtenu à partir d'un capteur vibratoire, notamment d'un capteur électrique, et notamment d'un capteur de vitesse. Ce signal peut être échantillonné en temps avec une fréquence fs d'acquisition du signal temporel, ou en angle, avec une résolution du signal angulaire R. Le signal x1 peut être représente par un vecteur x1[n] de taille N. 15 Ce signal est soumis au procédé de détection de la défaillance de l'élément mobile de l'aérogénérateur. Le procédé de détection mesure une grandeur physique représentative du mouvement de l'élément mobile de l'aérogénérateur de sorte à générer un signal mesuré temporel. Le procédé intègre également une étape d'amélioration de la qualité du 20 signal mesuré temporel. Cette étape permet de compléter le signal mesuré temporel par un nombre limité de valeur nulle, comme cela est représenté sur la figure 2. A titre d'exemple, selon un mode de mise en oeuvre, on augmente la taille d'un facteur k au moins égal au nombre d'harmonique pris en compte H, par exemple égal à 6. 25 On définit un vecteur m[n] de taille N unitaire multiplié par la valeur moyenne x1[n]. Le signal dont la qualité a été améliorée peut être représenté par un vecteur x2 de taille N'=kN tel que : x2 [n] = [(x1[n]- m[n]). w[n],Z[n],Z[n],Z[n],Z[n],Z[n][ 3058477 7 avec w[n], une fenêtre de pondération classique et Z[n], un vecteur nul de taille N. Dans le cadre de cet exemple, le vecteur x2[n] est donc 6 fois plus grand que le vecteur x1[n]. Ce vecteur est donc d'une taille N'=6N. 5 Cette étape d'amélioration peut également correspondre à une étape d'interpolation du spectre de façon à augmenter la résolution fréquentielle par un facteur correspondant à un nombre inférieur au nombre de multiples de la fréquence donnée. Cette étape d'interpolation est réalisée par ajout de valeurs correspondant à une valeur moyenne à l'extrémité du signal mesuré temporel pour l'ensemble de 10 valeur de façon à augmenter la longueur du signal mesuré temporel par un produit du signal temporel mesuré avec le nombre inférieur au nombre de multiples de la fréquence donnée. Une fois ces étapes éventuellement réalisées, comme l'illustre la figure 3, le signal mesuré temporel est transformé de sorte à générer un signal fréquentiel 15 présentant un premier spectre en fonction de la fréquence H. Plus précisément, le signal fréquentiel est généré en appliquant au signal mesuré temporel une transformée de Fourier pour des simplicités de calcul. Selon un exemple, par le moyen d'une simple transformée de Fourier rapide fft, ou éventuellement d'un traitement plus complexe, tel que l'analyse par 20 enveloppe ou plus généralement toute approche cyclo-stationnaire. X 2[n] = fft (x2 [n]) . En appelant df la résolution du spectre obtenu, df peut être défini par : - df = fs / N' dans le cas d'un échantillonnage temporel - df = R / N' dans le cas d'un échantillonnage angulaire 25 De plus, toujours dans une optique d'épurement du signal, des valeurs du spectre, correspondant à des fréquences parasites, sont remplacées par une valeur de lissage : valeur minimale du spectre ou valeur nulle. Ceci permet de conserver seulement les harmoniques nécessaires à la détection de la défaillance de l'état de l'élément mobile de l'aérogénérateur. 3058477 8 Il est possible de déterminer les composantes gênantes par une analyse de la cinématique du système, notamment des engrènements, en considérant les harmoniques, puis de remplacer dans le spectre toute les composantes gênantes du spectre par des valeurs nulles ou par la valeur minimale du spectre. On obtient ainsi 5 un spectre X3 [n] à partir du spectreX2[n]. Alternativement, le nettoyage du spectre peut être réalisé par un moyennage synchrone. Comme représenté sur les figures 5 à 9, il est ensuite réalisé une étape de décalage fréquentiel de f1 résultant en un spectre X4[n] dont la valeur pour la 10 fréquence donnée f correspond au produit de la valeur du premier spectre et d'un ensemble de valeurs du premier spectre correspondant à la somme de la fréquence donnée et d'un multiple de la fréquence caractéristique de l'élément mobile de l'aérogénérateur. Ainsi, dans le cadre de l'exemple précédemment cité, un spectre X4 est 15 obtenu. Pour obtenir un tel spectre, on liste les fréquences f, d'arbre susceptibles de moduler le signal. Pour chacune d'elle, on crée un spectre produit des spectres décalés (paramètre "M"). Le paramètre M correspond au nombre de modulations supérieures et inférieures qui sont prises en compte. Ainsi on obtient : X4, i[n] = X3[n] * prod_m(X3[n +m*fil df] * X3 [n - m * 20 fi / df], m allant de 1 à M) Ensuite, pour une fréquence donnée, et plus précisément pour la fréquence de l'élément mobile de l'aérogénérateur, il est calculé un spectre transformé dont la valeur pour cette fréquence donnée correspond au produit de la valeur du premier spectre et d'un ensemble de valeurs du premier spectre 25 correspondant à des H multiples de la fréquence donnée. Ce nombre de H-1 multiples de la fréquence donnée est défini par le nombre de H multiples de la fréquence donnée moins un. Ainsi, dans le cadre de l'exemple précédemment cité, chaque spectre est transformé en un spectre X5,1 par la technique décrite ci-dessus : X5, [n] = prod_h ( X4, i[n * h] , h allant de 1 à H) 3058477 9 Sur la base du spectre obtenu, il est possible de vérifier la réalisation d'un critère de détection, représenté à la figure 10, pour des valeurs du spectre transformé dans un intervalle fréquentiel à proximité d'une valeur caractéristique d'un défaut de sorte à détecter l'état de l'élément mobile de l'aérogénérateur lorsque l'élément 5 mobile de l'aérogénérateur est défaillant. Le critère de détection détecte la défaillance qui se manifeste par une vibration de faible énergie correspondant à un pic dans au moins un intervalle surveillé. La vérification du critère de détection peut consister par exemple à la vérification du dépassement d'un seuil d'amplitude dans l'intervalle surveillé. 10 Le ou les indicateur I décrits en référence à la figure 11 mentionnés précédemment peuvent correspondre notamment à une amplitude maximum atteinte par le spectre X5 sur un intervalle surveillé. Plusieurs intervalles peuvent être surveillés sur un même signal X5: un indicateur retiendra l'amplitude max de chacun d'eux, respectivement. 15 Ainsi, grâce aux dispositions décrites ci-dessus, une transformation est mise en oeuvre qui vise à transformer un spectre pour valoriser spécifiquement les familles d'harmoniques dont les éléments sont strictement les multiples entiers d'une "fréquence fondamentale". L'application de cette transformation permet de supprimer les 20 composantes qui ne font pas partie d'une famille d'harmonique nous évite d'hésiter entre une "composante fantôme" et une composante représentative de la défaillance. Le procédé permet d'adapter cette technique : - à un phénomène fréquentiel très régulier, ce qui est le cas des observations "angulaires" (i.e. à pas d'angle constant), indispensables pour surveiller 25 ce genre de phénomène ; en effet, cette technique combine l'énergie contenue dans plusieurs harmoniques d'une même fréquence. La possibilité d'utiliser des harmoniques de la fondamentale repose sur l'hypothèse que le phénomène observé soit très régulier ; - à un phénomène modulé par un peigne de raies grâce à l'étape de 30 décalage du spectre. 3058477 10 Ces dispositions sont notamment utiles par exemple dans le cas des fréquences caractéristiques de bague interne: suivre un défaut de bague interne revient à suivre les multiples des fréquences caractéristiques augmentées ou réduites de la fréquence caractéristique de l'arbre porteur du roulement. 5 REVENDICATIONS1) Procédé de détection de la défaillance d'un élément mobile d'un aérogénérateur, le procédé de détection comprend les étapes suivantes : mesure d'une grandeur physique représentative du mouvement de l'élément mobile de l'aérogénérateur de sorte à générer un signal mesuré temporel (xi) ; transformation du signal mesuré temporel (x1) ou d'un signal temporel (x2) obtenu à partir du signal mesuré temporel (x1) de sorte à générer un signal fréquentiel présentant un premier spectre (X2) ; calcul d'un spectre transformé (X5) dont la valeur pour une fréquence donnée correspond au produit de la valeur du premier spectre (X2) ou d'un spectre (X3, X4) obtenu à partir du premier spectre (X2) et d'un ensemble de valeurs du premier (X2) ou du spectre (X3, X4) obtenu à partir du premier spectre (X2) correspondant à des multiples de la fréquence donnée. vérification de la réalisation d'un critère de détection pour des valeurs du spectre transformé dans un intervalle fréquentiel à proximité d'une valeur caractéristique d'un défaut. 2) Procédé de détection selon la revendication 1, dans lequel l'étape de calcul d'un spectre transformé prend en compte un nombre de multiples (H) de la fréquence donnée correspondant à un nombre d'harmoniques d'une fréquence caractéristique du défaut à détecter. 3) Procédé de détection selon l'une quelconque des revendication 1 ou 2, lequel comprend une étape de décalage fréquentiel résultant en un spectre (X4) dont la valeur pour une fréquence donnée correspond au produit de la valeur du premier spectre (X2) ou d'un spectre (X3) obtenu à partir du premier spectre (X2) et d'un ensemble de valeurs du premier spectre (X2) ou d'un spectre (X3) obtenu à partir du premier spectre (X2) correspondant à la somme de la fréquence donnée et d'un multiple de la fréquence caractéristique de l'élément mobile de l'aérogénérateur. 3058477 12 4) Procédé de détection l'une quelconque des revendications 1 à 3, dans lequel le signal fréquentiel est généré en appliquant au signal mesuré temporel (x1) une transformée de fourrier. 5) Procédé de détection selon l'une quelconque des revendications 5 précédentes comprenant une étape d'amélioration de la qualité du signal mesuré temporel (x1) de façon à augmenter la résolution fréquentielle par un facteur (k) au moins égal au nombre de multiples (H). 6) Procédé de détection selon la revendication 5, dans lequel l'étape d'amélioration comprend une étape d'interpolation du spectre. 10 7) Procédé de détection selon la revendication 6, dans lequel l'étape d'interpolation est réalisée par ajout de valeurs correspondant à une valeur nulle ou à une valeur moyenne à l'extrémité du signal mesuré temporel (xi) pour l'ensemble de valeurs de façon à multiplier la longueur du signal mesuré temporel par le facteur (k) au moins égal au nombre de multiples (H). 15 8) Procédé de détection selon l'une quelconque des revendications précédentes comprenant une étape de nettoyage du spectre dans lequel au moins une valeur du spectre correspondant à au moins une fréquence parasite est remplacée par une valeur de lissage. 9) Procédé de détection selon la revendication précédente, dans lequel la 20 valeur de lissage correspond à la valeur minimale du spectre ou à une valeur nulle. 10) Procédé de détection selon la revendication précédente, dans lequel l'étape de nettoyage du spectre est réalisée préalablement à l'étape de décalage fréquentiel.