Brevet FR2842305A1 — publié le 20040116

PROC D DE DETERMINATION DE L'IMPEDANCE D'UN SYSTEME LECTROCHIMIQUE Pour la détermination des propriétés et des qualités d'un système électrochimique (tel qu'une pile, une batterie, un système de dépôt électrolytique, un système d'analyse d'un milieu) et la prévision de son fonctionnement futur, l'un des paramètres dont la connaissance est essentielle est l'impédance de ce système. De façon générale, un système électrochimique comprend, dans un milieu électrolytique, deux électrodes principales - une électrode de travail et une contre-électrode. Une électrode de référence est disposée au voisinage de l'électrode de travail et est utilisée en relation avec celle-ci pour effectuer diverses opérations de mesure ou de régulation. Dans la suite du présent texte les tensions et courants mentionnés correspondent à des mesures effectuées entre une électrode de travail et l'électrode de référence correspondante. L'impédance d'un système électrochimique correspond au rapport tension sur courant. Il est connu que pour mesurer rapidement la valeur de cette impédance en fonction de la fréquence, il est souhaitable de calculer le rapport des transformées de Fourier discrètes (DFT) de la tension et du courant. Pour cela, on applique entre une électrode de référence et une électrode de travail d'une cellule électrochimique un signal d'excitation de faible amplitude et on mesure à des intervalles de temps égaux AT, N valeurs successives e(n) et i(n) de la tension et du courant, avec 1 nombre de points de mesure, AT l'intervalle d'échantillonnage. Le calcul de transformée de Fourier discrète peut être effectué pour N/2 valeurs de fréquence fm, avec 0 réparties entre 0 et 1/2AT (0...1/mAT...1/2AT). L'impédance complexe pour une fréquence fm donnée est alors égale à: Z(fm) = E(fm)/I(fm). (3) Les méthodes de mesure d'impédance utilisées dans la pratique, diffèrent principalement par le type de signal d'exci20 tation du système: sinusodal, multi-sinusoidal, bruit blanc, etc. Une excitation sinusodale est de loin le procédé le plus utilisé car il s'avère être le plus précis. Le procédé consiste à imposer une régulation sinusodale à la cellule électrochimique avec un signal de petite amplitude et à enregistrer la réponse de la cellule en courant et en tension. Le rapport des transformées de Fourier de la tension et du courant à la fréquence du sinus donnera la valeur de l'impédance à cette fréquence. Le spectre des fréquences est balayé en modifiant la fréquence d'excitation. Dans la théorie des systèmes électriques, l'expression de l'impédance (équation (3)) est correcte à condition que le système analysé soit linéaire et stationnaire et que des perturbations ne soient pas apportées par des phénomènes extérieurs au système. Dans le cas des systèmes électrochimiques, le respect de ces conditions impose des précautions particulières. La fonction de transfert d'un système électrochimique est généralement non linéaire mais on peut la considérer comme linéaire sur une petite portion. C'est pour cela qu'on utilise un signal d'excitation de faible amplitude. En conséquence, le rapport signal/bruit diminue et on doit augmenter le temps de mesure, c'est-à-dire le nombre de points de mesure pour éliminer le bruit par une intégration de la réponse. En revanche, l'allonge10 ment du temps de mesure pose des problèmes quand à la stationnarité du système électrochimique. Dans beaucoup de cas réels, la condition de stationnarité ne peut pas être remplie. Les causes sont multiples: relaxation du potentiel, relaxation du courant, relaxation des concentrations. Dans ces conditions, si le temps de mesure est suffisamment long pour que l'on ait un bon rapport signal/bruit, le système n'est pas stationnaire pendant la durée de mesure et l'impédance calculée n'a plus grande signification, notamment aux basses fréquences (par exemple inférieures au hertz) o les problèmes de non-stationnarité sont particulière20 ment aigus. Divers procédés ont été proposés pour résoudre ce problème de non stationnarité des systèmes électrochimiques, mais aucun n'a donné de solution satisfaisante. Ainsi, un objet de la présente invention est de prévoir un nouveau procédé de calcul de l'impédance d'un système électrochimique permettant de s'abstraire des erreurs liées à la non stationnarité du système. Pour atteindre cet objet, la présente invention prévoit un procédé de détermination de l'impédance complexe Z(fm) d'un système électrochimique non stationnaire, comprenant les étapes suivantes: mettre le système dans un état de tension choisi et lui appliquer un signal sinusodal de fréquence fmi mesurer immédiatement après des valeurs successives de la tension et du courant à des intervalles de temps réguliers AT, calculer les transformées de Fourier discrètes de la tension (E(f>) et du courant (I(f)), la transformée de la tension étant calculée pour la seule fréquence fm du signal sinusodal et la transformée du courant étant calculée pour la fréquence fm et pour deux fréquences adjacentes fm-, et fm+l de part et d'autre de la fréquence fm, et calculer l'impédance selon la formule suivante Z(fm) = E(fm)/I*(fm) o I* désigne un courant corrigé tel que Re[I*(fm)] = Re[I(fm)] - {Re[I(fm+î)]+Re[I(fm-1)I}/2 Im[I*(fm)] = Im[I(fm)] - {Im[I(fm+î)]+Im[I(fm-)]1}/2. Un mode de réalisation particulier de la présente invention sera exposé à titre non-limitatif en relation avec la figure jointe qui représente le spectre d'amplitude en fonction de la fréquence d'un système électrochimique soumis à une excitation sinusodale de faible amplitude et à un échelon de tension. Tout d'abord, on notera que, comme le potentiel appliqué est imposé par l'instrument d'analyse, la non stationnarité du système ne peut se manifester que par des fluctuations de courant. La présente invention se base sur l'analyse du comportement d'un système électrochimique soumis à un échelon de tension. La DFT de la réponse en courant traduit pour l'essentiel des phénomènes de relaxation et donc l'effet de la non stationnarité du système. Comme le représente la courbe en pointillés sur la figure unique, le spectre en amplitude I du module de la DFT du courant présente habituellement une forte réponse à basse fréquence à des échelons de tension. Par contre, la réponse à une excitation sinusodale de fréquence fm d'un système stationnaire se traduit par une raie unique à la fréquence fm. En pratique, cette réponse correspond d'une part à la réaction du système à l'excitation à la fréquence fm, d'autre part à la contribution des effets de relaxation. L'avantage de l'utilisation de la transformée de Fourier est que le spectre monofréquence lié à l'excitation sinusodale est superposé au spectre lié à l'échelon de tension. Selon la présente invention, une fois les considérations ci-dessus prises en compte, on propose de soustraire de la réponse en intensité à la fréquence fm la contribution due aux relaxations, évaluée à partir de l'analyse de la DFT aux fréquences adjacentes à fm* En effet, la transformée de Fourier pour des fréquences fm-, et fm+l du système ne correspondra qu'aux phénomènes de non stationnarité et on considérera que la valeur de l'amplitude liée aux non stationnarités pour une fréquence fm est la moyenne des valeurs pour les deux fréquences adjacentes. Ainsi, la présente invention propose - de placer un système électrochimique dans un état de tension choisi et de lui appliquer une excitation sinusodale monofréquence de faible amplitude, - de mesurer, sans attendre la stabilisation du système, l'amplitude du courant et de la tension à des intervalles de temps réguliers immédiatement après l'application de la tension tandis que l'excitation sinusodale est appliquée, et - de calculer d'une part la transformée de Fourier discrète de la tension pour la fréquence fm de l'excitation et d'autre part la DFT du courant pour la valeur fm et pour deux valeurs de fréquence fm-, et fm+l adjacentes à la fréquence fm, et - de calculer la valeur de l'impédance complexe à la fréquence fm à partir de la valeur de la TFD de la tension pour la fréquence fm et d'une valeur corrigée de la TFD du courant tenant compte des TFD calculées pour les fréquences fm-1 fm et On calculera l'amplitude corrigée I* de la transformée de Fourier discrète du courant en valeur réelle et en valeur imaginaire par les deux équations suivantes: Re[I*(fm) ] = Re[I(fm)] - {Re[I(fm+î)]+Re[I(fm-1)]}/2 Im[I*(fm)] = Im[I(fm)] {Im[I(fm+î)]+Im[I(fm-1)]}/2 La valeur de l'impédance corrigée des effets de non stationnarité est alors: Z(fm) = E(fm)/I*(fm) avec I*(fm) = Re[I*(fm)] + jIm[I*(fm)] Dés le relevé des points effectué pour la fréquence fm, on pourra appliquer un signal sinusodal à une nouvelle fréquence et effectuer un nouveau relevé, et ainsi de suite. Un avantage du présent procédé est qu'il rend possible une analyse correcte d'impédance sur des systèmes électro15 chimiques non stationnaires notamment aux très basses fréquences. En même temps, un gain de temps considérable est obtenu pour des systèmes qui se stabilisent lentement car avec la correction selon l'invention, il n'est plus nécessaire d'attendre la stabilisation après une mise sous tension du système pour démarrer une analyse d'impédance. REVENDICATIONS 1. Procédé de détermination de l'impédance complexe Z(fm) d'un système électrochimique non stationnaire, caractérisé en ce qu'il comprend les étapes suivantes: mettre le système dans un état de tension choisi et lui appliquer un signal sinusodal de fréquence fm, mesurer immédiatement après des valeurs successives de la tension et du courant à des intervalles de temps réguliers AT, calculer les transformées de Fourier discrètes de la tension (E(f)) et du courant (I(f)), la transformée de la tension étant calculée pour la seule fréquence fm du signal sinusodal et la transformée du courant étant calculée pour la fréquence fm et pour deux fréquences adjacentes fm-, et fm+î de part et d'autre de la fréquence fmi et calculer l'impédance selon la formule suivante Z(fm) = E(fm)/I*(fm) o I* désigne un courant corrigé tel que Re[I*(fm)] = Re[I(fm)] {Re[I(fm+î)]+Re[I(fm-1)]}/2 Im[I*(fm)] = Im[I(fm)] - {Im[I(fm+î)]+Im[I(fm-î)]}/2. 2. Procédé selon la revendication 1, caractérisé en ce qu'il est répété pour une succession de fréquences d'excitation.