FR 3065289 A1 20181019 FR 1753267 A 20170413 1 Procédé de mesure de l'intensité d'un courant dans un conducteur La présente invention concerne un procédé de mesure de l'intensité d'un courant dans un conducteur électrique. La présente invention concerne également un procédé de conception d'un capteur de mesure. La présente invention concerne également un calculateur, un capteur de courant et un système de mesure associés. Des capteurs de mesure de courant sont utilisés dans de nombreuses applications pour évaluer l'intensité de courants traversant des conducteurs. Par exemple, de tels capteurs permettent une gestion des réseaux de distribution d'électricité et sont alors installés en différents points des réseaux afin d'évaluer la consommation effective des différentes installations alimentées et/ou de surveiller l'état général du réseau. Les informations récoltées sont ensuite utilisées pour gérer au mieux le réseau de distribution. De nombreux capteurs utilisés dans de telles applications évaluent le courant traversant un conducteur par des mesures de champs magnétiques créés par la circulation du courant dans le conducteur. Par exemple, des tores de Rogowski, ou encore des ensembles de magnétomètres sont disposés autour du conducteur pour estimer la valeur de l'intensité du courant. Cependant, les réseaux de distribution d'électricité sont fréquemment des réseaux polyphasés, en particulier triphasés, dans lesquels des courants circulent simultanément dans des conducteurs proches les uns des autres. Dans ces conditions, les champs magnétiques à proximité d'un conducteur sont modifiés de manière significative par des courants circulant dans les conducteurs voisins. La mesure, grâce à ces capteurs, de l'intensité du courant peut alors être moins précise. Il est alors nécessaire de réaliser un blindage de ces capteurs pour éviter ou au moins réduire l'influence des perturbations extérieures, ce qui complexifie la réalisation de ces capteurs et augmente leurs dimensions. En outre, les capteurs précités doivent être placés autour du conducteur à mesurer pour obtenir une précision suffisante. L'installation de ces capteurs est donc contraignante, et suppose parfois que les conducteurs soient détachés à une extrémité pour que le capteur soit enfilé autour du conducteur. Même lorsque ce n'est pas le cas, la mise en place est en général une opération complexe et longue qui implique fréquemment la coupure du courant dans le conducteur. La mise en place des capteurs est donc contraignante pour le gestionnaire du réseau puisqu'elle suppose une indisponibilité temporaire d'une partie du réseau. 3065289 2 Il existe donc un besoin pour un procédé de mesure de l'intensité d'un courant circulant dans un conducteur électrique, qui soit à la fois précis, de petites dimensions et plus facile à installer que les procédés de mesure de l'état de la technique. A cet effet, il est proposé un procédé de mesure de l'intensité d'un courant dans 5 un conducteur électrique, le procédé comprenant les étapes de : - obtention d'une matrice de calcul M comprenant un ensemble de termes, chaque terme d'une même ligne de la matrice de calcul étant proportionnel à la même distance entre un point de mesure et un point de référence, élevée à une puissance supérieure ou égale à zéro, 10 - formation d'un vecteur de mesure B comportant des composantes, chaque composante étant une mesure de la valeur du champ magnétique selon une direction en un point de mesure, chaque point de mesure associé à une direction étant distinct des autres points de mesure associés à la même direction, 15 - calcul des composantes d'un vecteur de coefficients A, chaque composante étant un coefficient, chaque coefficient étant un coefficient d'une décomposition du champ magnétique en harmoniques spatiales et étant proportionnel à l'intensité du courant pour que la relation matricielle B=MxA soit vérifiée, x dénotant l'opération de 20 multiplication matricielle, et - détermination de l'intensité du courant à partir au moins des composantes calculées. Selon des modes de réalisation particuliers, le procédé de mesure comprend une ou plusieurs des caractéristiques suivantes, prises isolément ou selon toutes les 25 combinaisons techniquement possibles : - les étapes d'obtention d'une matrice de calcul et de calcul sont fonction de paramètres du capteur, le procédé comprenant, en outre : - une étape d'optimisation comprenant l'obtention d'au moins un paramètre optimisé, tels que la position réelle du point de référence 30 par rapport au conducteur électrique, - une étape de calcul d'au moins un coefficient de conversion à partir du ou des paramètres optimisés, l'étape de détermination comprenant la mise en oeuvre de la division d'une valeur réelle d'un coefficient par le coefficient de conversion. - l'étape d'obtention est mise en oeuvre par recherche d'un jeu de paramètres optimisés minimisant l'écart entre au moins deux rapports entre deux composantes du 3065289 3 vecteur de coefficients obtenu lors de l'étape de calcul des composantes et deux rapports entre deux composantes intermédiaires, chaque composante intermédiaire étant une composante d'un vecteur de coefficients calculé pour un jeu de paramètres possibles du capteur. 5 - une première distance est définie entre le point de référence et le point du conducteur le plus proche du point de référence, le maximum des distances des points de mesure étant strictement inférieur à la première distance. - une première distance est définie entre le point de référence et le point du conducteur le plus proche du point de référence, le maximum des distances des points de 10 mesure étant inférieur ou égal à 85% de la première distance. - les points de mesure sont agencés le long d'un cercle de centre le point de référence ou d'une sphère de centre le point de référence. - une première distance est définie entre le point de référence et le point du conducteur le plus proche du point de référence, la première distance étant inférieure ou 15 égale à 10 centimètres et chaque puissance étant inférieure ou égale à 5. Il est également proposé un procédé de conception d'un capteur de mesure comportant une pluralité de détecteurs de champ magnétique, le procédé de conception déterminant, pour une précision de mesure souhaitée, le nombre de détecteurs de champ magnétique, ainsi que le positionnement spatial des détecteurs, en utilisant le procédé de 20 mesure défini ci-dessus. Il est également proposé un calculateur adapté pour mettre en oeuvre les étapes suivantes : - obtention d'une matrice de calcul M comprenant un ensemble de termes, chaque terme d'une même ligne de la matrice de calcul étant 25 proportionnel à la même distance d'un point de mesure, élevée à une puissance supérieure ou égale à zéro, - réception de mesures de la valeur du champ magnétique selon une direction en un point de mesure, chaque point de mesure associé à une direction étant distinct des autres points de mesure associés à la 30 même direction et à une distance d'un point de référence, dite distance du point de mesure, - formation d'un vecteur de mesure B comportant des composantes, chaque composante étant une mesure reçue, - calcul des composantes d'un vecteur de coefficients A, chaque 35 composante étant un coefficient, chaque coefficient étant un coefficient d'une décomposition du champ magnétique en harmoniques spatiales 3065289 4 et étant proportionnel à l'intensité du courant pour que la relation matricielle B=Mx A soit vérifiée, x dénotant l'opération de multiplication matricielle, et - détermination de l'intensité du courant à partir au moins des 5 composantes calculées. Il est également proposé un capteur de courant configuré pour mesurer une intensité d'un courant traversant un conducteur électrique, le capteur de courant comprenant une pluralité de détecteurs de champ magnétique et un calculateur tel que défini ci-dessus, chaque détecteur étant configuré pour mesurer la valeur du champ 10 magnétique selon une direction en un point de mesure, chaque point de mesure étant distinct et à une distance d'un point de référence, et pour transmettre la valeur mesurée au calculateur. Il est également proposé un système de mesure de courants traversant une pluralité de conducteurs électriques, le système de mesure comprenant une pluralité de 15 capteurs de courant tels que définis ci-dessus et un unique contrôleur central, chaque capteur étant configuré pour communiquer avec le contrôleur central. Des caractéristiques et avantages de l'invention apparaîtront à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif et faite en référence aux dessins annexés, sur lesquels : 20 - la figure 1 est une représentation schématique partielle d'un réseau de distribution d'électricité équipé d'un premier exemple de capteur de courant comprenant un ensemble de détecteurs de champ magnétique, - la figure 2 est un ordinogramme des étapes d'un procédé de mesure de courant mis en oeuvre par le capteur de la figure 1, comprenant une étape de calibration, 25 - la figure 3 est une représentation schématique d'un exemple d'ensemble de détecteurs de la figure 1, et - la figure 4 est une représentation schématique d'un ensemble de détecteurs de champ magnétique d'un deuxième exemple de capteur de courant, et - la figure 5 est un ordinogramme des étapes d'un exemple d'étape de calibration de 30 la figure 2. Un réseau 10 de distribution d'électricité a été partiellement représenté sur la figure 1. Le réseau 10 comprend une pluralité de conducteurs 15A, 15B. Selon l'exemple de la figure 1, la pluralité de conducteurs comprend un premier conducteur 15A et au 35 moins un deuxième conducteur 15B. 3065289 5 Le réseau 10 est, en outre, équipé d'au moins un premier exemple de capteur de courant 20. Chaque conducteur 15A, 15B est configuré pour recevoir un courant électrique d'un tableau de répartition et pour transmettre le courant électrique à une installation 5 électrique correspondante. Chaque courant électrique est, par exemple, un courant continu. En variante, au moins un courant électrique est un courant alternatif. Une intensité IA, IB et une tension électrique Ua, Ub sont définies pour chaque courant électrique traversant un conducteur 15A, 15B. 10 Chaque tension électrique Ua, Ub peut être comprise entre 50 volts (V) et 75 kilovolts (kV). Par exemple, chaque tension électrique est égale à 20 kV. Chaque capteur de courant 20 est configuré pour déterminer l'intensité IA, IB, d'au moins un courant électrique. Selon l'exemple de la figure 1, le capteur de courant 20 est configuré pour mesurer l'intensité IA du courant traversant le premier conducteur 15A. 15 Chaque capteur de courant 20 comporte un ensemble 25 de détecteurs 27 de champ magnétique et un calculateur 30. Dans un souci de clarté, un seul ensemble 25 a été représenté sur la figure 1, et est décrit dans la description ci-dessous. Cependant, il est à noter qu'un même capteur 20 est susceptible de comporter une pluralité d'ensembles 25. Par exemple, selon 20 une variante le capteur de courant 20 comporte deux ensembles 25, un pour chaque conducteur 15A, 15B. L'ensemble 25 est configuré pour déterminer un nombre Nv de valeurs de composantes du champ magnétique Cm. En particulier, l'ensemble 25 est configuré pour déterminer simultanément le nombre Nv de valeurs. Il est entendu par « simultanément » 25 que les valeurs sont déterminées au cours d'un même intervalle de temps. Chaque intervalle de temps est relativement court comparativement au temps de fluctuation des courants à mesurer. Par exemple, le nombre de valeurs Nv est supérieur ou égal à 10. Chaque détecteur 27 est configuré pour mesurer une valeur Gm en un point de 30 mesure Pm d'au moins une composante du champ magnétique Cm selon une direction de mesure Dm. Dans la suite de cette description, chaque point de mesure Pm sera repéré par un indice i allant de 1 à un premier nombre Ni. Selon une variante, chaque détecteur 27 est configuré pour mesurer des 35 valeurs Gm d'une pluralité de composantes du champ magnétique en un unique point de mesure Pm,. Par exemple, chaque détecteur 27 est configuré pour mesurer des 3065289 6 valeurs Gm de deux composantes du champ magnétique en un unique point de mesure Pm,. En variante, chaque détecteur 27 est configuré pour mesurer des valeurs Gm de trois composantes du champ magnétique en un unique point de mesure Pm,. 5 Chaque point de mesure Pm, associé à une direction de mesure Dm est distinct des autres points de mesure Pm, associés à la même direction de mesure Dm. Par exemple, l'ensemble 25 est configuré pour acquérir simultanément une unique valeur Gm d'une composante du champ magnétique Cm selon une direction de mesure Dm en chaque point de mesure Pm,. 10 Chaque direction de mesure Dm est, par exemple, définie de manière absolue. Des directions d'un repère cartésien sont des exemples de directions définies de manière absolue. En variante, chaque direction de mesure Dm est définie relativement au point de mesure considéré. Par exemple, des directions radiale et orthoradiale d'un repère muni 15 de coordonnées cylindriques sont des exemples de directions définies de manière relative. Selon un autre mode de réalisation, l'ensemble 25 est configuré pour acquérir simultanément une pluralité de valeurs Gm de composantes du champ magnétique Cm selon des directions de mesure Dm en chaque point de mesure Pm,, seules des valeurs 20 de composantes selon des directions de mesure Dm distinctes étant acquises simultanément en un même point de mesure Pm,. Chaque détecteur 27 est configuré pour transmettre les valeurs Gm mesurées au calculateur 30. Par exemple, chaque détecteur 27 est relié au calculateur 30 par une liaison filaire 32. En variante, chaque détecteur 27 est configuré pour transmettre au 25 calculateur 30 les valeurs Gm mesurées par communication radiofréquence. Le calculateur 30 est configuré pour mettre en oeuvre un procédé de mesure de l'intensité IA du courant dans le premier conducteur 15A. Le calculateur 30 comporte une mémoire 35 et un processeur 40. Le fonctionnement du capteur 20 va maintenant être décrit. Un ordinogramme des 30 étapes d'un procédé de mesure de l'intensité IA du courant traversant le premier conducteur 15A a été représenté sur la figure 2. Le procédé de mesure comprend une étape 100 de calibration, une étape 110 d'obtention d'une matrice, une étape 120 de formation, une étape 130 de calcul et une étape 140 de détermination. 35 Au cours du procédé de mesure, l'ensemble 25 est dans une position de fonctionnement. 3065289 7 Lorsque l'ensemble 25 est dans la position de fonctionnement, l'ensemble 25 est disposé à proximité du premier conducteur 15A. Par exemple, l'ensemble 25 est fixé au premier conducteur 15A par un bras de fixation dédié. 5 Lorsque l'ensemble 25 est dans la position de fonctionnement, l'ensemble 25 est distant de chaque conducteur 15A, 15B. En particulier, l'ensemble 25 est distant du premier conducteur 15A. Il est entendu par « distant » qu'aucun point du premier conducteur 15A n'appartient à une surface plane délimitée par trois segments reliant toute combinaison de trois points de mesure Pm,, ni à un segment reliant toute combinaison de 10 deux points de mesure Pm,. En particulier, les volumes définis par l'ensemble des points de mesure Pm, et le conducteur sont disjoints. Une première distance dl est définie entre le point de référence Pr et le point du premier conducteur 15A le plus proche du point de référence Pr. Le maximum des distances r, des points de mesure Pm,, autour de Pr, est strictement inférieur à la première 15 distance dl. En d'autres termes, le capteur 20 est physiquement disjoint du premier conducteur 15A. Par exemple, le maximum des distances r, des points de mesure Pm, est inférieur ou égal à 85 pourcents (`)/0) de la première distance dl. Une deuxième distance d2 est définie entre le détecteur 27 le plus proche du premier conducteur 15A et le point du premier conducteur 15A le plus proche du 20 détecteur 27 considéré. Une troisième distance d3 est définie entre le point de référence Pr et le point le du second conducteur 15B le plus proche du point de référence Pr. La troisième distance d3 est supérieure ou égale à la première distance dl. Au cours de l'étape de calibration 100, un ensemble de coefficients de 25 conversion K est calculé. Chaque coefficient de conversion K est un nombre réel. L'étape 100 est, par exemple, mise en oeuvre une seule fois lorsque le capteur 20 est installé dans sa position de fonctionnement. En variante, l'étape 100 est mise en oeuvre lorsque des indices donnent à penser que le capteur 20 a été déplacé par rapport aux conducteurs 15A, 15B. 30 Lors de l'étape 110 d'obtention d'une matrice, une matrice de calcul M est obtenue par le calculateur 30. Par exemple la matrice de calcul M est mémorisée dans la mémoire 35. La matrice de calcul M est calculée en fonction de la position relative des différents points de mesure Pm par rapport au point de référence Pr. Par exemple, l'étape 35 d'obtention 110 est mise en oeuvre une seule fois lors de la conception du capteur 20. La matrice de calcul M comprend un ensemble de termes T. 3065289 8 Les termes T sont agencés en un ensemble de lignes L et de colonnes C de la matrice de calcul M. La matrice de calcul M comprend le nombre de valeurs Nv de lignes L. La matrice de calcul comprend un deuxième nombre N2 de colonnes C. Chaque ligne L comprend 5 ainsi le deuxième nombre N2 de termes T. Chaque terme T de la matrice de calcul M est une fonction élémentaire d'une décomposition en harmoniques spatiales du champ magnétique Cm. Il est entendu par « décomposition en harmoniques spatiales » qu'une valeur de chaque composante du champ magnétique Cm est susceptible d'être calculée, en tout 10 point de l'espace, comme une somme infinie de produits. Chaque produit intervient entre un coefficient et la fonction élémentaire correspondante. Chaque coefficient est un nombre. Chaque fonction élémentaire est une fonction définie pour tout point de l'espace et dépendant de la position du point considéré dans l'espace. 15 Selon le principe de la décomposition en harmoniques spatiales, les fonctions élémentaires sont identiques quelles que soit le champ magnétique Cm considéré, mais chaque champ magnétique Cm est entièrement défini par les valeurs correspondantes des coefficients. En d'autres termes, les fonctions élémentaires forment une base dans laquelle le champ magnétique Cm est susceptible d'être décrit. 20 En particulier, chaque fonction élémentaire est une fonction de la position du point considéré par rapport à un point de référence Pr. Par exemple, chaque fonction élémentaire est une fonction d'une puissance de la distance r, entre le point considéré et le point de référence. Dans la suite de cette description, le terme « ordre » est utilisé pour signifier « puissance » au sens de la phrase précédente. 25 Ainsi, chaque terme T est une fonction de la position d'un point de mesure Pm, correspondant par rapport au point de référence Pr. Par exemple, chaque terme T d'une même ligne L est une fonction de la position d'un même point de mesure Pm, correspondant à un détecteur 27. En particulier, chaque terme T est proportionnel à une distance r, entre le point de 30 mesure Pm, correspondant et le point de référence Pr, élevée à une puissance supérieure ou égale à zéro. Un ordre maximum N est défini comme étant l'ordre maximum des termes de la décomposition en harmoniques spatiales du champ magnétique. En d'autres termes, la somme infinie de termes est tronquée pour donner une approximation suffisante du 35 champ magnétique avec un nombre fini de termes, dépendant de l'ordre maximum N. L'ordre N est, par exemple, égal à 4. 3065289 9 Par exemple, chaque ligne L de la matrice de calcul M comprend autant de termes T qu'il a été défini de fonctions élémentaires associées aux ordres entre 0 et l'ordre maximum N. Chaque terme T d'une même colonne C est proportionnel à une distance r, élevée 5 à une puissance identique pour tous les termes T de la colonne C. Au cours de l'étape de formation 120, chaque détecteur 27 mesure une valeur Gm d'une composante de champ magnétique en le point de mesure Pm, correspondant. Le calculateur 30 reçoit les valeurs Gm mesurées des détecteurs 27. Un vecteur de mesure B est alors formé. 10 Par exemple, le vecteur de mesure B comprend le nombre de valeurs Nv de composantes bl, où I est un indice. Chaque composante ID, est une valeur Gm d'une composante de champ magnétique mesurée. Chaque composante ID, est associée à une ligne de la matrice de calcul M et est identifiée par l'indice I correspondant. Les composantes ID, sont ordonnées dans le vecteur B dans l'ordre des indices I 15 croissants. Chaque composante ID, est mesurée en le point de mesure Pm, associé à la ligne L correspondante. Au cours de l'étape 130 de calcul, un vecteur de coefficients A est calculé par le processeur 40. Le vecteur de coefficient A comprend le deuxième nombre N2 de composantes a. 20 Les valeurs réelles des composantes a du vecteur A de coefficients sont telles que la relation matricielle : B=MxA (Equation 1) dans laquelle « x» dénote l'opération de multiplication matricielle, est vérifiée. En d'autres termes, les composantes du vecteur de coefficients A sont calculées 25 par la résolution du problème inverse de l'équation 1 ci-dessus. Ainsi, le capteur 20 étant disjoint du premier conducteur 15A, l'équation 1 correspond à la résolution d'un problème interne avec une décomposition en harmoniques spatiales en puissance croissante de la distance r, par rapport au point de référence Pr. 30 L'équation 1 décrit une décomposition en harmoniques spatiales des valeurs de champ magnétique mesurées, dans laquelle seuls un nombre de produits égal au résultat de la multiplication du nombre de valeurs Nv et du deuxième nombre N2 sont calculés. Chaque composante a est alors un coefficient de la décomposition du champ magnétique Cm en harmoniques spatiales. 35 Chaque composante a calculée par l'équation 1 précédente est également appelée valeur réelle de la composante a. 3065289 10 Chaque composante a du vecteur de coefficients A est uniquement proportionnelle au courant IA traversant le premier conducteur 15A lorsqu'aucune autre source de champ magnétique n'est présente, par exemple lorsqu'aucun courant IB ne traverse le deuxième conducteur 15B et que l'environnement magnétique est négligeable. Chaque coefficient 5 de conversion K est calculé au cours de l'étape de calibration 100 à partir d'une fonction mathématique reliant une valeur d'une composante a et le courant IA traversant le premier conducteur 15A. Lors de l'étape de détermination 140, l'intensité IA du courant traversant le premier conducteur 15A est déterminée à partir des composantes a calculées du vecteur de 10 coefficients A. La détermination de l'intensité IA comprend la division d'une valeur réelle d'un coefficient a par un coefficient de conversion K correspondant. L'intensité IA est, par exemple, égale au résultat de la division d'une valeur réelle d'un coefficient a associé à la plus haute valeur de puissance, c'est-à-dire à l'ordre 15 maximal N, par le coefficient de conversion K correspondant. Ainsi, lorsque les coefficients de conversion K sont connus et les composantes a déterminées grâce à l'équation 1 précédente, le courant IA traversant le premier conducteur 15A est calculé précisément. Le procédé de mesure mis en oeuvre par le capteur 20 utilise avantageusement le 20 fait que les sources éloignées de l'ensemble 25 présentent une contribution plus faible aux coefficients que les sources proches de l'ensemble 25, et ce d'autant plus que l'ordre associé au coefficient a augmente. En particulier, le procédé de mesure est moins sensible au courant IB traversant le deuxième conducteur 15B, lorsque le deuxième conducteur 15B est plus éloigné de l'ensemble 25. Cette sensibilité est d'autant plus faible 25 que les coefficients a sont associés à un ordre élevé. Ainsi, le procédé présente une meilleure précision que les procédés de mesure de l'état de la technique, en s'affranchissant mathématiquement des sources de champ autres que le courant IA que l'on cherche à mesurer, tout en évitant tout besoin de mise en oeuvre d'un système de blindage. 30 En outre, le procédé de mesure est susceptible d'être mis en oeuvre par un capteur 20 disposé à proximité du conducteur 15A dont l'intensité du courant est à mesurer mais n'entourant pas ce conducteur 15A. La mise en place du capteur 20 est ainsi facilitée. En particulier, la mise en place du capteur 30 ne suppose pas une coupure du courant circulant dans le conducteur 15A étudié. 35 Le premier exemple a été décrit dans le cas où l'intensité IA est déterminée par la division d'une valeur réelle d'un coefficient a associé à la plus haute valeur de la 3065289 11 puissance, c'est-à-dire à l'ordre maximum N, par le coefficient de conversion K correspondant. Il est à noter que l'intensité IA est susceptible d'être obtenue à partir d'ordres différents de l'ordre maximum N. Cela est en particulier le cas si plusieurs ensembles 25 distincts et présentant des points de référence Pr distincts sont intégrés à 5 un même capteur 20. Un deuxième exemple d'ensemble 25 de détecteurs 27 a été représenté de manière schématique, dans sa position de fonctionnement, sur la figure 3. Chaque conducteur 15A, 15B s'étend selon une direction principale Z Chaque conducteur 15A, 15B présente une forme selon des sections 10 perpendiculaires à la direction principale Z. Par exemple, chaque conducteur 15A, 15B est cylindrique à base circulaire. Dans ce cas, l'intersection d'un conducteur 15A, 15B avec un plan perpendiculaire à la direction principale Z est un disque. La réunion des centres de chaque disque d'un même conducteur 15A, 15B forme alors un segment parallèle à la direction principale Z et appelé ligne centrale Lc du conducteur 15A, 15B 15 considéré. Une longueur, selon la direction principale Z, est définie pour chaque conducteur. La longueur est, par exemple 1 mètre (m) de part et d'autre de l'ensemble 25. De telles dimensions permettent de considérer le champ magnétique comme invariant selon la direction principale Z. 20 Une direction X et une direction Y ont été représentées sur la figure 3. Les directions X et Y sont perpendiculaires l'une à l'autre ainsi qu'a la direction principale Z et forment avec la direction principale Z un repère cartésien centré sur le point de référence Pr. La position de la ligne centrale Lc du premier conducteur 15A, est repérée dans un 25 repère polaire, de centre le point de référence Pr, associé au repère (XYZ) par ses coordonnées polaires (dO, 90). Ici dO est la distance, mesurée dans un plan perpendiculaire à la direction principale Z, entre le point de référence et un point appartenant à la ligne centrale Lc du premier conducteur 15A, et 90 est l'angle entre la direction X et un segment reliant le point de référence Pr au point appartenant à la ligne 30 centrale Lc précitée. La position dans l'espace de la ligne centrale Lc du premier conducteur 15A par rapport au point de référence Pr est, par exemple, une position supposée dans laquelle le bras est prévu pour fixer le point de référence Pr. En d'autres termes, les valeurs de l'angle (p0 et de la distance dO sont fixées par construction de l'ensemble de capteurs 25 35 et de leur support, et par le choix des directions Y et Z utilisées pour le procédé de mesure. 3065289 12 L'angle (p0 est, par exemple, égal à 180 degrés (°). La distance dO est, par exemple, égale à 10 centimètres (cm). La position de chaque point de mesure Pm, est définie, par rapport au point de référence Pr, par un ensemble d'un angle 0, et d'une distance r, mesurée dans un plan 5 perpendiculaire à la direction principale Z entre le point de mesure Pm, et le point de référence Pr. L'angle 0, est, par exemple, l'angle entre la direction X et un segment reliant le point de mesure Pm, et le point de référence Pr. Afin de clarifier la figure 3, seuls six détecteurs 27 ont été représentés. Selon l'exemple de la figure 3, les points de mesure Pm, sont agencés le long d'un 10 cercle. Par exemple, le cercle est compris dans un plan perpendiculaire à la direction principale Z. Par exemple, les points de mesure Pm, sont angulairement équirépartis le long du cercle, c'est-à-dire qu'un angle entre un point de mesure Pm,, le centre du cercle et l'un des deux points de mesure Pm, voisins est identique à l'angle entre le point de mesure 15 Pm, considéré, le centre du cercle et l'autre point de mesure Pm, voisin. Le cercle présente, par exemple, un rayon compris entre 2 cm et 5 cm. Il est à noter que différentes valeurs de rayon sont envisageables en fonction de l'agencement des conducteurs 15A, 15B et des courants IA, IB considérés. Le point de référence Pr est, par exemple, le centre du cercle sur lequel sont agencés les points de mesure Pm,. 20 Selon l'exemple de la figure 3, chaque direction de mesure Dm est définie dans un repère dépendant du point de mesure Pm, correspondant. Par exemple, chaque direction de mesure Dm est une direction orthoradiale ou une direction radiale. Une direction orthoradiale est une direction perpendiculaire à un segment reliant le centre du cercle au point de mesure Pm, considéré et comprise dans le 25 plan du cercle. Une direction radiale est une direction définie par une droite reliant le centre du cercle au point de mesure Pm,. Sur la figure 3, une direction radiale Dmr et une direction orthoradiale Drne ont été représentées pour chaque point de mesure Pm,. Selon l'exemple de la figure 3, chaque détecteur 27 configuré pour mesurer une 30 composante de champ magnétique selon une direction orthoradiale Drne est associé à un point de mesure Pm, distinct des autres points de mesure Pm, associés à des composantes de champ magnétique selon des directions orthoradiales. Par exemple, chaque détecteur 27 est configuré pour mesurer une valeur Gme d'une composante selon une direction orthoradiale Drne et une valeur Gmr d'une 35 composante selon une direction radiale Dmr d'un champ magnétique en un unique point de mesure Pm,. 3065289 13 Ainsi, chaque détecteur 27 est configuré pour mesurer deux valeurs Gm de composantes, selon deux directions différentes, en un unique point de mesure Pm,. Dans un système de coordonnées polaires, chaque coefficient de conversion est soit une fonction d'un sinus de l'angle 90 soit une fonction d'un cosinus de l'angle 90 5 Dans la suite, chaque coefficient de conversion dépendant d'un cosinus sera identifié par une référence Ka, et chaque coefficient de conversion dépendant d'un sinus sera identifié par une référence Kbn. Il est à noter que des fonctions d'un cosinus d'un nombre sont susceptibles d'être exprimées sous la forme d'un sinus d'une différence entre Tr/2 et ce nombre. De même, 10 une fonction d'un sinus est susceptible d'être exprimée sous la forme d'un cosinus d'une différence entre Tr/2 et ce nombre. Il est entendu par « fonction d'un cosinus » une fonction dépendant d'un cosinus du nombre considéré et non d'une différence entre Tr/2 et ce nombre. Il est entendu par « fonction d'un sinus » une fonction dépendant d'un sinus du nombre considéré et non d'une différence entre Tr/2 et ce nombre. 15 Chaque coefficient de conversion Kan, Kb, est identifié par un indice n variant de 1 à l'ordre N. De même, chaque composante a est soit une fonction d'un sinus de l'angle 90 soit une fonction d'un cosinus de l'angle angle 90. Dans la suite, chaque composante a dépendant d'un cosinus sera identifiée par une référence a, et chaque composante a 20 dépendant d'un sinus sera identifiée par une référence bn. Chaque composante, an, bn est identifiée par un indice n variant de 1 à l'ordre N. Chaque coefficient de conversion Kan, Kb, est conforme dans un cas théorique idéal à l'une des équations 2 et 3 ci-dessous : cos(n q;10) Un - Kan (Equation 2) n(dOr 25 U sin(")°) - Kb (Equation 3) n(dOr dans lesquelles Un est une constante réelle arbitraire liée aux conditions pratiques de mise en oeuvre du capteur 20. Selon le principe de la décomposition en harmoniques spatiales, les valeurs des composantes radiale Gmr et orthoradiale Gme du champ magnétique Cm en un point de 30 mesure Pm, sont alors susceptibles de s'écrire selon les équations 4 et 5 ci-dessous : 00 Gmr(ri ,z9i)=In.rin-1[- an sin(n z9i ) + bn cos(n z9i )1 (Equation 4) n=i 3065289 14 Gm29(r,29,)= c>n.rin l[an cos(n ) + bn sin(n )] (Equation 5) n=i dans lesquelles a, et b, sont des coefficients de la décomposition en harmoniques spatiales, Z est le signe d'une sommation, 0, et r, sont les coordonnées polaires du point de mesure Pm, et le point de référence Pr, cos désigne la fonction cosinus et sin désigne 5 la fonction sinus. Chaque terme T est une fonction d'un cosinus de l'angle 0, et / ou d'un sinus de l'angle 0,. Par exemple, la matrice M est conforme à l'équation : 10 M= ( cos(t9) sin(z91) NriN-1 COS(Nei) -Nr/1 Sin(Nei) NriN-1 sin(Ne) NriN-1 COS(Nei) -sin(t91) COS(ei) NrN-1 COS(NeNi) Nr,N-1 cos( eNi ) sin( eNi ) Ni sin(NeNi) sin(NeNi) - sin( eNi ) cos( eNi ) (Equation 6) Ni 1 COS(NeN1)) Lorsque la matrice M est conforme à l'équation 6, le vecteur B est, par exemple, égal à: ( GmOi Gmrl B= (Equation 7) GnION-1 Gmr \ N-1 ) Dans ce cas, l'équation 1 ci-dessus correspond à une troncature à l'ordre maximal 15 N de la décomposition en harmoniques spatiales du champ magnétique mesuré. Il est alors possible de démontrer que les coefficients a, et b, des équations 4 et 5 précédentes sont conformes aux équations 8 et 9 ci-dessous lorsque le courant IA traversant le premier conducteur 15A est la seule source de champ présente : cos(n IA= Ka a = (Equation 8) n n n(dO)n sin(n ça0) IA= Kbn.IA 20 b = (Equation 9). n n n(dO)n Les équations 8 et 9 ci-dessus permettent de calculer des valeurs supposées des coefficients a, et b,, dans un cas théorique idéal où la constante U, est arbitrairement unitaire, les coordonnées cylindriques dO, 90 par rapport au point de référence Pr et la 3065289 15 position dans l'espace des points de mesure Pm, sont parfaitement connus. On constate alors que chaque coefficient a, ou b, est effectivement proportionnel au courant IA à déterminer. Le facteur de proportion est le coefficient de conversion Kan, Kbn correspondant. 5 Par souci de simplicité, les exemples précédents ont été décrits ci-dessus dans un cas où un unique capteur 20 est utilisé pour mesurer l'intensité IA d'un courant traversant un seul conducteur 15A. Dans le cas où un deuxième conducteur 15B parcouru par le courant IB est présent, celui-ci contribue également au champ magnétique en chaque point de mesure 10 Pmi de façon similaire au courant IA traversant le premier conducteur 15A. Ainsi, la contribution du courant IB traversant le deuxième conducteur 15B est conforme aux équations 4 et 5 ci-dessus, les contributions aux coefficients a, et b, étant calculés selon les équations 8 et 9 en remplaçant les coordonnées polaires (dO, 90) du centre du premier conducteur 15A par les coordonnées polaires (db, 9b) du centre du deuxième 15 conducteur 15B. Cette contribution s'additionne linéairement à celle du premier conducteur 15A. On constate que comme la distance db entre le point de référence Pr et le centre du deuxième conducteur 15B est supérieure strictement à la distance dO entre le point de référence Pr et le centre du premier conducteur 15A, la contribution due au courant IB 20 traversant le deuxième conducteur 15B dans le calcul des coefficients a, et b, est plus faible que pour le courant IA traversant le premier conducteur 15A. Cette différence est d'autant plus importante que l'ordre n considéré est grand. Ainsi, pour un indice n assez grand, la contribution du courant IB est suffisamment faible pour considérer avec une bonne précision que les coefficients a, et b, sont 25 uniquement proportionnels au courant IA traversant le premier conducteur 15A, même si le courant IB traversant le deuxième conducteur est non nul. Cela est également vrai même pour toutes les autres sources extérieures qui seraient plus éloignées du point de référence Pr que les conducteurs 15A, 15B. Ainsi, en sélectionnant les ordres les plus élevés, les contributions des sources 30 sont alors filtrées spatialement, pour ne conserver que les contributions dues à la source la plus proche, à savoir ici le courant IA du premier conducteur 15A. Il est à noter que pour des raisons de contrainte de mise en oeuvre pratique, par exemple, si l'ordre maximum N est trop faible pour pleinement décrire le champ magnétique ou si le nombre de mesures Nv effectuées est insuffisant, la constante U, 35 présente dans les équations 8 et 9 peut s'éloigner de la valeur unité. Cependant, les coefficients a restent linéairement dépendants de l'intensité IA par exemple si seul le 3065289 16 conducteur 15A est étudié. Ainsi, il reste possible d'obtenir avec une bonne précision la valeur de l'intensité IA. Une procédure de calibration permettant d'obtenir une valeur précise des coefficients de conversion K est décrite ultérieurement. 5 Selon une variante, le capteur 20 fait partie d'un système de mesure comportant une pluralité de capteurs 20 et un unique contrôleur central. Le calculateur 30 de chaque capteur 20 est configuré pour communiquer avec le contrôleur central. Par exemple, chaque capteur 20 est connecté au contrôleur central par une liaison filaire. 10 En variante, chaque capteur 20 est configuré pour communiquer avec le contrôleur central par communication radiofréquence. Cela signifie que le calculateur 30 est propre à émettre ou à recevoir des données transmises via un signal comportant au moins une onde électromagnétique radiofréquence. Les ondes électromagnétiques radiofréquences sont les ondes 15 électromagnétiques présentant une fréquence entre 3 kilohertz (KHz) et 3 gigahertz (GHz). Le calculateur 30 de chaque capteur 20 est configuré pour transmettre au contrôleur central les valeurs des intensités IA, IB mesurées. Le contrôleur central est configuré pour mémoriser les valeurs des intensités IA, IB mesurées. Le contrôleur central 20 est, en outre, configuré pour afficher ou transmettre les valeurs des intensités IA, IB mesurées à destination d'un utilisateur. Selon une autre variante, les valeurs Gmr, Gme mesurées en chaque point de mesure Pmi sont transmises au contrôleur central, qui met en oeuvre les étapes d'obtention 110, de formation 120, de calcul 130 et de détermination 140 pour chaque 25 capteur 20. Selon un mode de réalisation particulier, le capteur 20 est obtenu par un procédé de conception utilisant le procédé de mesure ci-dessus. En particulier, le nombre de valeurs Nv, et donc le nombre de détecteurs 27, sont déterminés, par le procédé de conception, en fonction d'une précision de mesure 30 souhaitée. Le positionnement spatial des détecteurs 27 est également déterminé en utilisant le procédé de mesure. En particulier, le positionnement spatial des détecteurs 27 est choisi de manière à simplifier les calculs effectués lors des étapes 120 de formation, 130 de calcul et une 140 de détermination. 3065289 17 Par exemple, l'orientation de chaque détecteur 27 est choisie de manière à ce que les valeurs Gmr, Gme de champ magnétique mesurées soient aisément transposables dans un système de coordonnées choisi. Selon un mode de réalisation, l'ordre maximum est inférieur ou égal à 5 pour une 5 distance dO inférieure ou égale à 10 cm. En pratique, un ordre maximum N égal à 4, et une distance dO égale à 6 centimètres entre le centre du cercle sur lequel sont disposés les points de mesure Pmi et le centre du premier conducteur 15A, permettent d'obtenir une précision en général inférieure à 0,5% lorsque le rayon r, du cercle est égal à 3,2 cm, même en cas de forte perturbation. 10 Le deuxième exemple ci-dessus a été donné dans le cas où le champ magnétique est décrit de manière simple dans un système de coordonnées polaires. C'est en particulier le cas de conducteurs 15A, 15B rectilignes, pour lesquels on peut considérer que le champ magnétique Cm généré par les conducteurs 15A, 15B ne varie pas en fonction de la position selon la direction principale Z. Ainsi, toutes les composantes de 15 champ magnétique Cm mesurées sont comprises dans un même plan, qui est le plan du cercle sur lequel sont agencés les points de mesure Pm,. L'ensemble 25 de détecteurs est alors aisé à réaliser. En outre, les calculs effectués lors de la mise en oeuvre du procédé de mesure sont relativement simples. Le procédé de mesure est alors peu gourmand en termes de 20 ressources de calcul. Cependant, il est à noter que d'autres systèmes de coordonnées sont susceptibles d'être utilisés pour des géométries de conducteurs différentes, en utilisant une base de fonctions adaptée aux systèmes de coordonnées utilisés pour les équations 1 à 9 ci-dessus. 25 Par exemple, un troisième exemple de capteur 20 a été représenté sur la figure 4. Les éléments identiques au deuxième exemple de la figure 3 ne sont pas décrits à nouveau. Seules les différences sont mises en évidence. Les points de mesure Pm, sont agencés sur une sphère dont le centre est le point de référence Pr. Par exemple, les points de mesure Pm, sont répartis sur un ensemble de 30 parallèles P de la sphère. Selon un mode de réalisation, l'ensemble 25 comporte 30 détecteurs 27 répartis sur 5 parallèles P. Par exemple, six points de mesure Pm, sont agencés sur chaque parallèle P. Les points de mesure Pm, agencés le long d'un même parallèle P sont, par exemple, équirépartis le long du parallèle P. 35 Chaque parallèle P est un cercle inscrit dans la sphère, les centres des parallèles P étant alignés sur une droite D passant par le centre de la sphère. Chaque parallèle P 3065289 18 est défini par l'intersection d'un plan et de la sphère, les plans étant parallèles les uns aux autres et décalés les uns par rapport aux autres selon la direction de la droite D. La droite D est, par exemple, parallèle avec la direction principale Z. Chaque détecteur 27 est configuré pour mesurer des valeurs Gm de composantes 5 de champ magnétique en un même point de mesure Pm, selon trois directions de mesure Dm différentes. Les directions de mesure Dm sont chacune parallèle à un vecteur de base d'un repère de coordonnées sphériques associé au point Pm, considéré. Dans le repère sphérique, les coordonnées spatiales de chaque point de mesure Pm, sont le rayon r, et 10 deux angles cp, et 0,. Le rayon r, est la distance entre le point de mesure Pmi et le point de référence Pr. L'angle 0, est un angle entre la direction principale Z et un segment reliant le point de mesure Pm, au point de référence Pr. L'angle cp, est un angle entre la direction X et la projection de ce même segment sur un plan perpendiculaire à la direction principale Z. 15 Afin d'assurer la lisibilité de la figure 4, un seul point de mesure Pm y est représenté, avec les directions de mesure Dm et les coordonnées associées. Les indices i des références associées n'y sont donc pas représentés. Les trois directions de mesure Dm sont, par exemple deux directions tangentielles Dg), DO et une direction radiale Dr. 20 Chaque direction tangentielle Dg), DO est une direction tangente à la sphère en le point de mesure Pm, correspondant. La direction tangentielle Dg) est perpendiculaire à un plan contenant la direction principale Z et un segment reliant le point de mesure Pm, correspondant au point de référence Pr. Chaque direction radiale Dr est une direction parallèle au segment reliant le point de mesure Pm, correspondant au point de 25 référence Pr. La direction tangentielle DO est perpendiculaire à la direction tangentielle Dg) et à la direction radiale Dr. Le troisième exemple de capteur 20 est alors adapté à des cas où les conducteurs 15A, 15B ne sont pas rectilignes et où le champ magnétique Cm présente donc des composantes significatives dans les trois directions de mesure Dm. 30 Du fait que le champ magnétique Cm est tridimensionnel, les fonctions élémentaires pour le calcul de M sont plus complexes mais ce troisième exemple est adaptable à un plus grand nombre de configurations des conducteurs 15A, 15B. Le champ magnétique Cm est alors conforme à l'équation : N n 011(r il VI) n ,n1V m[rin nm (9 pi)] (Equation 10) n=1 m=-n 35 dans laquelle : 3065289 19 - V est l'opérateur différentiel Nabla, - ail, est un coefficient de la décomposition en harmoniques spatiales dans le système de coordonnées sphériques et - Ynm sont des fonctions harmoniques sphériques. 5 Il est à noter que, pour un ordre maximum N donné, un plus grand nombre N2 de termes T sont calculés par rapport au deuxième exemple précédent et un plus grand nombre de valeurs Gm du champ magnétique sont mesurées. Selon un quatrième exemple, chaque courant est un courant variable. Par exemple, chaque courant est une phase d'un courant triphasé. 10 Les étapes 120 de formation, 130 de calcul et 140 de détermination sont alors itérées dans cet ordre avec une période d'itération prédéterminée. La période d'itération correspond à une fréquence d'échantillonnage des valeurs Gm du champ magnétique. De préférence, les étapes 120 de formation, 130 de calcul et 140 de détermination sont itérées sur une durée temporelle correspondant à plusieurs périodes du courant. 15 Une transformation de Fourier est ensuite appliquée aux valeurs d'intensité IA mesurées pour obtenir une valeur efficace de l'intensité IA. La transformation de Fourier est, par exemple, obtenue avec un algorithme de transformation de Fourier rapide (également appelé FFT, de l'anglais « Fast Fourier Transform » qui signifie « transformation de Fourier rapide »). 20 La transformation de Fourier est, par exemple, une transformation de Fourier temporelle, c'est-à-dire en fonction du temps. En variante, au cours de chaque étape de formation 120, une transformation de Fourier est appliquée à un ensemble de vecteurs de mesure B obtenus au cours de l'étape de formation 120 et d'une pluralité d'étapes de formation 120 antérieures pour 25 calculer une transformée de Fourier TB du vecteur de mesure B. L'étape 130 de calcul est alors mise en oeuvre avec la transformée de Fourier TB au lieu du vecteur de mesure B, pour obtenir une transformée de Fourier du vecteur de coefficients A. Selon un cinquième exemple, l'étape de calibration 100 comprend une étape 150 30 d'optimisation et une étape 160 de calcul d'un coefficient de conversion Kan, Kbn. Il est à noter que dans les équations 1 à 10 ci-dessus, un certain nombre de paramètres du capteur 20 sont utilisés, tels que des coordonnées spatiales. Ces paramètres sont susceptibles de varier en fonction de la fabrication du capteur 20 ou au cours de la vie de celui-ci, par exemple à la suite de chocs ou en cas de mauvais 35 positionnement de l'ensemble 25 par rapport au premier conducteur 15A, ou du vieillissement de matériaux. 3065289 20 Au cours de l'étape d'optimisation 150, des paramètres du capteur 20 sont optimisés pour améliorer la précision des calculs ultérieurs. Les paramètres optimisés sont, en particulier, des paramètres géométriques du capteur 20. En particulier, une position réelle du centre du premier conducteur 15A par rapport au point de référence Pr 5 est obtenue. Au cours de l'étape d'optimisation 150, l'étape 110 d'obtention d'une matrice, l'étape 120 de formation et l'étape 130 de calcul sont mises en oeuvre pour une valeur quelconque de l'intensité IA. Cela est représenté sur la figure 5 par des étapes 110c, 120c et 130c. 10 Selon un mode de mise en oeuvre, l'étape d'optimisation 150 est mise en oeuvre pour une valeur de l'intensité IA suffisamment élevée pour s'affranchir d'un éventuel bruit magnétique. Par exemple, la position réelle du centre du premier conducteur 15A par rapport au point de référence Pr est obtenue par minimisation d'un écart pour une pluralité de 15 positions courantes Pc. Plus précisément, un ensemble de coefficients a est calculé pour un ensemble de positions du centre du premier conducteur 15A par rapport au point de référence Pr, appelées positions courantes Pc. Chaque position courante Pc est une position dans laquelle est susceptible de se trouver par rapport au point de référence Pr, suite à un 20 mauvais placement ou à une imprécision dans la construction du capteur 20, ou encore suite à une modification de la géométrie du capteur 20 ultérieure à la fabrication de celui-ci, par exemple à cause d'un choc ou d'une usure, le centre du premier conducteur 15A. Chaque position courante Pc est, par exemple, choisie dans une gamme de positions prédéterminée. 25 Dans le cas où la position supposée du centre du premier conducteur 15A par rapport au point de référence Pr est définie, dans un système de coordonnées cylindriques, par une distance dO et par un angle (p0, chaque position courante Pc est, par exemple, une paire d'un angle courant (pc et d'une distance courante dc. Un ensemble de coefficients appelés coefficients intermédiaires a, est alors 30 calculé, au cours de l'étape d'optimisation 150 pour chaque jeu de paramètres. Au cours de l'étape 110c d'obtention, la matrice M est calculée pour chaque position courante Pc. Lors de l'étape 120c de formation le vecteur B est calculé par la loi de Biot et Savart et non plus mesuré en fonction de l'intensité IA du courant. 3065289 21 Au cours de l'étape 130c, les coefficients intermédiaires a, sont calculés selon l'équation 1, en utilisant la matrice M et le vecteur B calculés au cours des étapes 110c et 120c. Par exemple, une paire de coefficients intermédiaires a, est calculée pour chaque 5 position courante Pc. Le nombre de positions courantes Pc utilisées est déterminé en fonction de la configuration du ou des conducteurs 15A, 15B, du mode de fixation de l'ensemble 25 au premier conducteur 15A et de la précision souhaitée pour la mesure de l'intensité IA. Par exemple, les positions courantes Pc correspondent à des distances dO entre 10 le point de référence Pr et le centre du premier conducteur 15A comprises entre 5 cm et 50 cm. Selon un mode de mise en oeuvre particulier, l'ensemble des distances dO correspondant aux positions courantes Pc couvre la gamme entre 5 cm et 15 cm avec un pas de 1 millimètre (mm). L'ensemble des angles 90 entre la direction X et un segment reliant le point de 15 référence Pr au centre du premier conducteur 15A couvre, par exemple, la gamme entre 175° et 185°, avec un pas de 0,5°. Selon un mode de réalisation, un ensemble de rapports intermédiaires Ra entre des coefficients intermédiaires a, est calculé. Chaque rapport intermédiaire Ra est calculé par la division d'un numérateur par un 20 dénominateur, le numérateur étant un coefficient intermédiaire a, calculé en supposant une position courante Pc et pour une valeur donnée de l'indice n et le dénominateur étant un autre coefficient intermédiaire a, calculé pour la même position courante Pc. Par exemple, le dénominateur et le numérateur correspondent à des valeurs différentes de l'indice n. En variante, le dénominateur et le numérateur peuvent également correspondre 25 à des valeurs identiques de l'indice n. Par exemple, chaque rapport intermédiaire Ra est obtenu par la division de coefficients intermédiaires a, calculés pour des valeurs de l'indice n dont la différence est égale à 1. Lorsque deux coefficients intermédiaires a,, notés ac, par la suite, sont conformes 30 à l'équation 8 précédente, le rapport intermédiaire Ra correspondant est identifié par une référence Rcan, en l'absence de source de champ magnétique autre que le premier conducteur 15A, est alors conforme à l'équation: Rcan acn+i = U1 1 n cos((n+1)- çoc) (Equation 11) acn U n dc (n +1) cos(n - çac) On note que, même en présence d'autres sources de champ magnétique que le 35 conducteur 15A par exemple, l'équation 11 reste valide à partir d'un ordre N suffisamment 3065289 22 grand pour pouvoir négliger l'influence de ces autres sources, en fonction du niveau de précision souhaité. Selon une autre variante, au moins un rapport intermédiaire Ra, identifié par une référence Rcc,, est calculé à partir d'un coefficient intermédiaire, notés acn, conforme à 5 l'équation 8 précédente et un coefficient intermédiaire, notés bcn, conforme à l'équation 9 précédente. Selon un mode de réalisation, pour chaque position courante Pc, un rapport courant Ra est calculé pour chaque valeur de l'indice n allant de 1 à N. En outre, pour chaque valeur de l'indice n allant de 1 à N, un rapport réel Rra, 10 entre deux composantes du vecteur de coefficient A est calculé selon l'équation : a n+i (Equation 12) Rra = n n De même, pour chaque valeur de l'indice n allant de 1 à N, un rapport réel Rrc, entre deux composantes du vecteur de coefficient A est calculé selon l'équation : bn+1 (Equation 13) Rrc n = an 15 La position réelle est alors obtenue par optimisation de la position courante Pc par minimisation de l'écart entre un rapport réel Rr et un rapport courant Rc correspondant. Par exemple, des valeurs absolues d'un ensemble d'écarts entre des rapports réels Rra, Rrc, et des rapports intermédiaires Rcan, Rcc, sont minimisées. Selon un mode de réalisation, la minimisation porte sur la somme d'un ensemble 20 de valeurs absolues d'écarts entre des rapports réels Rr et les rapports intermédiaires Rc correspondants. Selon un mode de réalisation, la minimisation est définie par l'expression mathématique : N -1 a n+i Rca n I Rcc n (Equation 14) Min an n Pc n=1 25 Ainsi, la position réelle est la position courante Pc pour laquelle l'expression mathématique de l'équation 14 est la plus petite. L'étape d'optimisation 150 ci-dessus a été décrite dans le cas de l'optimisation de la position du centre du premier conducteur 15A par rapport au point de référence Pr. Il est à noter que d'autres paramètres du capteur 20 sont susceptibles d'être ainsi 30 optimisés. Par exemple, des coordonnées spatiales, 0, et r, d'un point de mesure Pm, sont susceptibles d'être ainsi optimisées. Les positions courantes Pc sont alors des positions que le point de mesure Pm, considéré est susceptible d'occuper. 3065289 23 De manière générale, le jeu de paramètres optimisés est le jeu de paramètres qui est conforme à l'expression mathématique de l'équation 14. Ainsi, à l'issue de l'étape 150 d'optimisation, un jeu de paramètres optimisés est obtenu. 5 Lors de l'étape 160 de calcul d'un coefficient de conversion, un coefficient de conversion K est calculé en divisant le coefficient a correspondant obtenu avec le jeu de paramètres, par l'intensité IA utilisée lors de la calibration. Au cours de l'étape 140 de détermination, l'intensité IA est calculée à partir d'un coefficient a obtenu à l'étape de calcul 130 et du coefficient de conversion K 10 correspondant calculé lors de l'étape de calibration 100. Le procédé de mesure est alors plus précis, car les coefficients de conversion K sont déterminés sur site et en conditions réelles, avec précision au lieu d'être supposés connus. Le procédé de mesure est donc moins sensible à d'éventuelles erreurs de placement de l'ensemble 25 de détecteurs 27 par rapport au premier conducteur 15A. 15 En outre, l'étape de calibration est effectuée à partir de rapports entre des coefficients an, bn. De tels rapports ne dépendent pas de la valeur de l'intensité IA circulant dans le premier conducteur 15A. La calibration est donc susceptible d'être effectuée à n'importe quel moment, car elle ne suppose pas de fixer l'intensité IA à une valeur prédéterminée. En d'autres termes, l'étape de calibration 100, elle aussi, ne 20 perturbe pas le fonctionnement du réseau 10. L'étape de calibration 110 a été décrite pour le cinquième exemple ci-dessus dans un cas où le champ magnétique Cm est bidimensionnel, correspondant par exemple à la configuration du deuxième exemple. Il est à noter que cette étape de calibration est susceptible d'être mise en oeuvre pour d'autres configurations moyennant une adaptation 25 correspondante des équations 11 à 13. Selon une variante de l'étape de calibration 110, l'intensité du courant circulant dans le premier conducteur 15A est connue, et le ou les coefficients de conversion Kan, Kb, sont calculés selon l'une des équations 8 et 9 ci-dessus. L'invention correspond à toute combinaison techniquement possible des modes de 30 réalisation précédents. REVENDICATIONS1.- Procédé de mesure de l'intensité (IA) d'un courant dans un conducteur électrique (15A), le procédé comprenant les étapes de : - obtention (110) d'une matrice de calcul M comprenant un ensemble de termes (T), chaque terme (T) d'une même ligne de la matrice de calcul (M) étant proportionnel à la même distance (ri) entre un point de mesure (Pm,) et un point de référence (Pr), élevée à une puissance supérieure ou égale à zéro, - formation (120) d'un vecteur de mesure B comportant des composantes (b1), chaque composante (b1) étant une mesure de la valeur du champ magnétique selon une direction (Dmr, Dme, DO, Dg), Dr) en un point de mesure (Pm,), chaque point de mesure (Pm,) associé à une direction (Dmr, Dme, DO, Dg), Dr) étant distinct des autres points de mesure (Pm,) associés à la même direction (Dmr, Dme, DO, Dg), Dr), - calcul (130) des composantes (a) d'un vecteur de coefficients A, chaque composante (a) étant un coefficient, chaque coefficient (a) étant un coefficient d'une décomposition du champ magnétique en harmoniques spatiales et étant proportionnel à l'intensité (IA) du courant pour que la relation matricielle B=MxA soit vérifiée, x dénotant l'opération de multiplication matricielle, et - détermination (140) de l'intensité (IA) du courant à partir au moins des composantes (a) calculées. 2.- Procédé selon la revendication 1, dans lequel les étapes d'obtention (110) d'une matrice de calcul (M) et de calcul (130) sont fonction de paramètres du capteur (20), le procédé comprenant, en outre : - une étape d'optimisation (150) comprenant l'obtention d'au moins un paramètre optimisé, tels que la position réelle du point de référence (Pr) par rapport au conducteur électrique (15A), - une étape de calcul (160) d'au moins un coefficient de conversion (K) à partir du ou des paramètres optimisés, - l'étape de détermination (140) comprenant la mise en oeuvre de la division d'une valeur réelle d'un coefficient (a) par le coefficient de conversion (K). 3.- Procédé selon la revendication 2, dans lequel l'étape d'obtention (160) est mise en oeuvre par recherche d'un jeu de paramètres optimisés minimisant l'écart entre au moins deux rapports entre deux composantes (a) du vecteur de coefficients (A) obtenu 3065289 25 lors de l'étape (130) de calcul des composantes (a) et deux rapports entre deux composantes intermédiaires (a,), chaque composante intermédiaire (a,) étant une composante d'un vecteur de coefficients (A) calculé pour un jeu de paramètres possibles du capteur (20). 5 4.- Procédé selon l'une quelconque des revendications 1 à 3, dans lequel une première distance (dl) est définie entre le point de référence (Pr) et le point du conducteur (15A) le plus proche du point de référence (Pr), le maximum des distances (r,) des points de mesure (Pm,) étant strictement inférieur à la première distance (dl). 10 5.- Procédé selon l'une quelconque des revendications 1 à 4, dans lequel une première distance (dl) est définie entre le point de référence (Pr) et le point du conducteur (15A) le plus proche du point de référence (Pr), le maximum des distances des points de mesure (Pm,) étant inférieur ou égal à 85% de la première distance (dl ). 15 6.- Procédé selon l'une quelconque des revendications 1 à 5, dans lequel les points de mesure (Pm,) sont agencés le long d'un cercle de centre le point de référence (Pr) ou d'une sphère de centre le point de référence (Pr). 20 7. - Procédé selon l'une quelconque des revendications 1 à 6, dans lequel une première distance (dl) est définie entre le point de référence (Pr) et le point du conducteur (15A) le plus proche du point de référence (Pr), la première distance (dl) ) étant inférieure ou égale à 10 centimètres et chaque puissance étant inférieure ou égale à 5. 25 8.- Procédé de conception d'un capteur de mesure comportant une pluralité de détecteurs (27) de champ magnétique, le procédé de conception déterminant, pour une précision de mesure souhaitée, le nombre de détecteurs (27) de champ magnétique, ainsi que le positionnement spatial des détecteurs (27), en utilisant le procédé de mesure selon 30 l'une quelconque des revendications 1 à 7. 9.- Calculateur (30) adapté pour mettre en oeuvre les étapes suivantes : - obtention (110) d'une matrice de calcul M comprenant un ensemble de termes, chaque terme d'une même ligne de la matrice de calcul (M) étant proportionnel à 35 la même distance (r) d'un point de mesure (Pm,), élevée à une puissance supérieure ou égale à zéro, 3065289 26 - réception (120) de mesures de la valeur du champ magnétique selon une direction (Dmr, Dme, DO, Dg), Dr) en un point de mesure (Pm,), chaque point de mesure (Pm,) associé à une direction (Dmr, Dme, DO, Dg), Dr) étant distinct des autres points de mesure (Pm,) associés à la même direction (Dmr, Dme, DO, Dg), 5 Dr) et à une distance (r) d'un point de référence (Pr), dite distance du point de mesure, - formation (120) d'un vecteur de mesure B comportant des composantes (b1), chaque composante étant une mesure reçue, - calcul (130) des composantes (a) d'un vecteur de coefficients A, chaque 10 composante (a) étant un coefficient, chaque coefficient (a) étant un coefficient d'une décomposition du champ magnétique en harmoniques spatiales et étant proportionnel à l'intensité (IA) d'un courant traversant un conducteur électrique (15A, 15B) pour que la relation matricielle B=Mx A soit vérifiée, x dénotant l'opération de multiplication matricielle, et 15 - détermination (140) de l'intensité (IA) du courant à partir au moins des composantes calculées. 10.- Capteur de courant (20) configuré pour mesurer une intensité (IA) d'un courant traversant un conducteur électrique (15A), le capteur de courant (20) comprenant 20 une pluralité de détecteurs de champ magnétique (27) et un calculateur (30) selon la revendication 9, chaque détecteur (27) étant configuré pour mesurer la valeur du champ magnétique selon une direction (Dmr, Dme, DO, Dg), Dr) en un point de mesure (Pm,), chaque point de mesure (Pm,) étant distinct et à une distance (r,) d'un point de référence (Pr), et pour transmettre la valeur mesurée au calculateur (30). 25 11.- Système de mesure de courants (IA, IB) traversant une pluralité de conducteurs électriques (15A, 15B), le système de mesure comprenant une pluralité de capteurs de courant (20) selon la revendication 10 et un unique contrôleur central, chaque capteur (20) étant configuré pour communiquer avec le contrôleur central.