FR 3035496 A1 20161028 FR 1553646 A 20150423 1 Procédé et système de détermination de la structure d'un réseau de distribution d'électricité et programme d'ordinateur associé La présente invention concerne un procédé mis en oeuvre par ordinateur, de détermination de la structure d'un réseau de distribution d'électricité, un programme d'ordinateur associé et un système de détermination de la structure du réseau de distribution d'électricité. Dans le domaine de la distribution d'électricité, il est connu d'assurer la distribution d'électricité à une pluralité de consommateurs électriques via un réseau de distribution d'électricité comprenant un poste de transformation électrique, le poste comportant des départs électriques d'alimentation en énergie électrique des consommateurs électriques. Le poste de transformation est, par exemple, un poste de transformation haute tension/moyenne tension (HTB/HTA) ou moyenne tension/basse tension (HTA/BT) et les départs électriques sont des départs en moyenne ou en basse tension. La haute tension, correspond à une tension supérieure à 50 kVolts (kV), la moyenne tension correspond à une tension comprise entre 1 kV et 50 kV et la basse tension correspond à une tension inférieure à 1 kV. Dans de tels réseaux de distribution électrique, les consommateurs sont alimentés soit en triphasé, soit en monophasé, et le poste de transformation est configuré pour répartir la puissance électrique qu'il reçoit entre les différents départs électriques. Cependant, la structure, et plus précisément la cartographie, des réseaux de distribution d'électricité alimentant les consommateurs électriques est généralement mal connue. Or, la connaissance de la structure de tels réseaux est devenue essentielle pour répondre aux contraintes que les régulateurs nationaux imposent aux distributeurs d'électricité, notamment en termes de qualité de fourniture d'électricité. En effet, les régulateurs nationaux imposent aux distributeurs d'électricité des objectifs de qualité de service, et ceux-ci doivent être mesurables. Les distributeurs d'électricité ont donc dorénavant besoin de mieux connaître la structure des réseaux de distribution d'électricité, notamment pour localiser les éventuels défauts ou défaillances sur le réseau ou encore pour identifier les départs du poste de transformation les plus contributeurs en termes de pertes électriques. Il est ainsi connu de EP 2 458 340 A2 un procédé permettant de déterminer la structure d'un réseau de distribution d'électricité à partir de moyens de calcul spécifiques et d'informations relatives à l'énergie électrique consommée par chaque consommateur électrique relié au réseau et relatives à l'énergie électrique délivrée par chaque départ du poste de transformation. 3035496 2 Cependant, la rapidité d'exécution, la précision et la fiabilité d'un tel procédé restent à améliorer. Le but de l'invention est donc de proposer un procédé de détermination robuste de la structure d'un réseau de distribution d'électricité permettant de déterminer la structure 5 du réseau de distribution d'électricité de manière plus sûre, plus fiable et plus rapide. Un tel procédé permet notamment de s'affranchir d'erreurs, par exemple liées à des imprécisions de mesure de l'énergie consommée et/ou délivrée ou à l'absence de certaines mesures d'énergie. A cet effet, l'invention a pour objet un procédé de détermination de la structure 10 d'un réseau de distribution d'électricité le procédé étant mis en oeuvre par ordinateur, le réseau de distribution comprenant un poste d'alimentation comportant un ou plusieurs départs d'alimentation en énergie électrique de plusieurs consommateurs électriques et comprenant les étapes suivantes : -a) l'acquisition, pour chaque consommateur électrique, de premières données 15 relatives à l'énergie électrique consommée par le consommateur électrique pendant différents intervalles temporels, mesurées via un premier capteur de mesure d'énergie, -b) l'acquisition, pour chaque départ, de deuxièmes données relatives à l'énergie électrique délivrée par le départ pendant les différents intervalles temporels, mesurées via un deuxième capteur de mesure d'énergie, 20 - c) la génération, via un module de génération, de plusieurs jeux de données différents, associés chacun à l'un des intervalles temporels, chaque jeu de données comprenant les premières et deuxièmes données associées audit intervalle temporel, caractérisé en ce que le procédé comprend en outre les étapes suivantes : -d) le calcul d'un premier critère de sélection pour chaque jeu de données, le 25 premier critère de sélection étant un critère choisi parmi un taux global de pertes d'énergie électrique entre le ou les départs et les consommateurs et un écart de consommation d'énergie électrique entre les différents consommateurs, -e) la sélection, par une unité de sélection et en fonction des premiers critères de sélection calculés, d'un premier ensemble de jeux de données parmi les jeux générés, 30 -f) la détermination, à partir du premier ensemble sélectionné, de paramètres de raccordement, lesdits paramètres comportant pour chaque consommateur électrique un identifiant du départ auquel il est raccordé.. Grâce à l'invention, les étapes de sélection du premier ensemble de jeux de données en fonction des premiers critères de sélection calculés et de détermination des 35 paramètres de raccordement à partir du premier ensemble sélectionné permettent de 3035496 3 déterminer de manière plus fiable et plus rapide la structure du réseau de distribution d'électricité. Selon d'autres aspects avantageux de l'invention, le procédé de détermination comprend en outre une ou plusieurs des caractéristiques suivantes, prises isolément ou 5 selon toutes les combinaisons techniquement admissibles : - lors de l'étape de calcul d), un deuxième critère de sélection pour chaque jeu de données généré est calculé, le deuxième critère étant l'autre critère choisi parmi un taux global de pertes d'énergie électrique entre le ou les départs et les consommateurs et un écart de consommation d'énergie électrique entre les différents consommateurs, 10 et, lors de l'étape de sélection e), le premier ensemble est sélectionné en fonction des premiers et deuxièmes critères calculés ; - suite à l'étape de sélection e), le procédé comprend les étapes suivantes : -el ) le calcul d'un deuxième critère de sélection pour chaque jeu de données du premier ensemble, le deuxième critère étant l'autre critère 15 choisi parmi un taux global de pertes d'énergie électrique entre le ou les départs et les consommateurs et un écart de consommation d'énergie électrique entre les différents consommateurs, -e2) la sélection en fonction des deuxièmes critères de sélection calculés d'un deuxième ensemble de jeux de données parmi le premier ensemble, 20 la détermination des paramètres de raccordement étant alors effectuée lors de l'étape de détermination f) en fonction du deuxième ensemble sélectionné parmi le premier ensemble ; - lors de l'étape de détermination f), les paramètres de raccordement sont déterminés séparément pour chaque départ et indépendamment des deuxièmes données 25 relatives à l'énergie électrique acquises pour les autres départs ; - lors de l'étape de détermination f), un système d'équation à résoudre est déterminé à partir d'un postulat de conservation de l'énergie électrique pour chaque départ selon lequel l'énergie délivrée par le départ est sensiblement égale à la somme de l'énergie consommée par les consommateurs électriques connectés au départ et des 30 pertes électriques ; - lors de l'étape de détermination f), le système d'équation est défini à partir de l'équation suivante : EDJk) .E(tk) + -2,k), 3035496 4 où n est le nombre de consommateurs électriques, j est un indice de départ variant de 1 à m, m étant le nombre de départs, i est un indice de consommateur électrique variant de 1 à n, k est un indice de jeu de données du premier ensemble variant de 1 à w, w étant le 5 nombre de jeux de données du premier ensemble, tk est l'intervalle temporel associé au jeu de données d'indice k, est le paramètre de raccordement indiquant si le consommateur d'indice i est raccordé au départ d'indice j, Eak) est l'énergie électrique consommée par le consommateur électrique d'indice i 10 pendant l'intervalle temporel tk, E Di (tk) est l'énergie électrique délivrée par le départ d'indice j pendant l'intervalle temporel tk, et ai et a 2,k sont des variables d'ajustement représentatives des pertes électriques ; - le système d'équation s'écrit sous une forme matricielle de la manière suivante : 15 A*z= B, où z= a11 B= EDI (ti et A= U 0 - - . 0 avec an1 EDI ED2 o - - - 0 U alm ED2 (tw) EDm(ti) 0 anm EDm(tk) aln aim,' 2ii _a 2, E1(t1) - - - E(t1) E1(t) - - - Ecn(t w) - lors l'étape de détermination, un algorithme d'optimisation est utilisé pour résoudre le système d'équation, l'algorithme d'optimisation vérifiant : {A.z = B 20 min fT .Z tel que , Où fT .z=[0 --- 0 1 --- 1 1 --- 11x z , lz z u z ,i*' w*. w*. U= et A comprenant m*n colonnes et m*w lignes ; 3035496 5 , avec e compris entre 0 et 0,5; 0 +00 - chaque départ comprend un ou plusieurs conducteurs électriques de phase, et, lors de l'étape d'acquisition b), les deuxièmes données sont acquises pour chaque conducteur électrique de phase et sont relatives à l'énergie électrique délivrée par chaque 5 conducteur électrique de phase, et, lors de l'étape de détermination f), les paramètres de raccordement déterminés comportent pour chaque consommateur électrique un identifiant du ou des conducteurs électriques de phase auxquels il est raccordé ; - lors de l'étape de détermination f), le procédé comprend les étapes suivantes : 41) la sélection pseudo-aléatoire d'un troisième ensemble de jeux de 10 données inclus dans le premier ensemble, -f2) la détermination des paramètres de raccordement à partir du troisième ensemble sélectionné parmi le premier ensemble ; - suite à l'étape f1) de sélection et précédemment à l'étape de détermination f2), le procédé comprend les étapes suivantes : 15 -f3) le calcul d'un paramètre de répétition, -f4) la détermination, à partir du troisième ensemble sélectionné, d'indices intermédiaires de raccordement comportant pour chaque consommateur électrique un identifiant du départ auquel il est raccordé, -f5) la mémorisation des indices intermédiaires de raccordement, et 20 -f6) la comparaison du paramètre de répétition avec au moins un critère de répétition prédéterminés, et, suivant le résultat de l'étape de comparaison f6), la réalisation d'une étape parmi, d'une part, la répétition des étapes de sélection f1), de calcul f3), de détermination f4), de mémorisation f5) et de comparaison f6) et, d'autre part, l'étape de détermination f2) des 25 paramètres de raccordement en fonction des indices intermédiaires de raccordement mémorisés ; le procédé comprend l'étape suivante : -f7) le calcul, pour chaque consommateur et pour chaque départ, d'un nombre d'affectation relatif au nombre de fois où les indices intermédiaires 30 de raccordement mémorisés indiquent que le consommateur est raccordé au départ, -e et u, 1+e -e 0 1+e +00 3035496 6 et, lors de l'étape de détermination f2), les paramètres de raccordement sont déterminés à partir des nombres d'affectation calculés. L'invention a également pour objet un programme d'ordinateur comportant des instructions logicielles, qui lorsqu'elles sont exécutées par un ordinateur mettent en oeuvre 5 le procédé défini ci-dessus. L'invention a, en outre, pour objet un système de détermination de la structure d'un réseau de distribution d'électricité, le réseau de distribution comprenant un poste d'alimentation comportant un ou plusieurs départs d'alimentation en énergie électrique de plusieurs consommateurs électriques, le système comprenant : 10 - pour chaque consommateur électrique, un premier capteur de mesure de premières données relatives à l'énergie électrique consommée par le consommateur pendant différents intervalles temporels, - pour chaque départ, un deuxième capteur de mesure de deuxièmes données relatives à l'énergie électrique délivrée par le départ pendant les différents intervalles 15 temporels, - un module d'acquisition des premières et deuxièmes données, - un module de génération de plusieurs jeux de données différents, associés chacun à l'un des intervalles temporels, chaque jeux de données comprenant les premières et deuxièmes données associées audit intervalle temporel, 20 caractérisé en ce que le système de détermination comprend en outre : - un module de calcul d'un premier critère de sélection pour chaque jeu de données, le premier critère de sélection étant un critère choisi parmi un taux global de pertes d'énergie électrique entre le ou les départs et les consommateurs et un écart de consommation d'énergie électrique entre les différents consommateurs, 25 - un module de sélection, en fonction des premiers critères de sélection calculés, d'un premier ensemble de jeux de données parmi les jeux générés, et - un module détermination, à partir du premier ensemble sélectionné, de paramètres de raccordement, lesdits paramètres comportant pour chaque consommateur électrique un identifiant du départ auquel il est raccordé. 30 L'invention sera mieux comprise et d'autres avantages de celle-ci apparaitront plus clairement à la lumière de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif, et faite en se référant aux dessins sur lesquels : - la figure 1 est une représentation schématique d'un réseau de distribution d'électricité comprenant un poste d'alimentation comportant plusieurs 35 départs d'alimentation en énergie électrique de plusieurs consommateurs électriques ; 3035496 7 - la figure 2 est un organigramme d'un procédé de détermination de la structure du réseau de distribution d'électricité de la figure 1, selon un premier mode de réalisation de l'invention ; - la figure 3 est un organigramme analogue à celui de la figure 2, selon un 5 deuxième mode de réalisation de l'invention ; - la figure 4 est un organigramme analogue à celui de la figure 2, selon un troisième mode de réalisation de l'invention ; et - la figure 5 est un organigramme analogue à celui de la figure 2, selon un quatrième mode de réalisation de l'invention. 10 Sur la figure 1, un réseau 10 de distribution d'électricité est associé à un système 12 de détermination de la structure du réseau de distribution d'électricité 10. Le réseau de distribution 10 comprend un poste d'alimentation 14 en énergie électrique de plusieurs consommateurs électriques 16,. Le poste d'alimentation 14 comporte des départs électriques 18, d'alimentation en énergie électrique des consommateurs électriques 16,, 15 avec i un indice de consommateur électrique et j un indice de départ. L'indice i de consommateur électrique varie de 1 à n, où n est le nombre de consommateurs électriques 16, et l'indice de départ j varie de 1 à m, où m est le nombre de départs électriques 18,. Dans l'exemple de la figure 1, le nombre n de consommateurs électriques 16, est égal à 3, et le nombre m de départs électriques 18, est égal à 2. 20 Le système de détermination 12 comprend, pour chaque consommateur électrique 16,, un premier capteur 20, de mesure d'énergie et un module 21, de transmission des mesures d'énergie réalisées par le premier capteur 20, correspondant. Le système de détermination 12 comporte, pour chaque départ 18,, un deuxième capteur 22, de mesure d'énergie. 25 Le système de détermination 12 comprend également une unité de traitement d'informations 24 formée par exemple d'un processeur 26 et d'une mémoire 28 associée au processeur 26. Le poste d'alimentation 14 est, par exemple, un poste de transformation moyenne tension/basse tension relié entre un réseau moyenne tension, non représenté, et un 30 réseau basse tension 30 correspondant aux consommateurs électriques 16,. Le poste d'alimentation 14 inclut, au niveau de chaque départ l8, le deuxième capteur 22, correspondant. Les consommateurs électriques 16, sont raccordés au poste d'alimentation 14 via les départs 18,. Plus précisément, dans l'exemple de la figure 1, les consommateurs 35 électriques 161, 162 sont reliés au départ 181 et le consommateur électrique 163 est relié au départ 182. 3035496 8 Les consommateurs électriques 16, sont soit triphasés et alimentés par le départ 18, correspondant via quatre conducteurs électriques 32,, 34,, 36,, 38,, c'est-à-dire trois conducteurs électriques de phase 32,, 34,, 36, et un conducteur électrique de neutre 38,, soit monophasés et alimentés par le départ 18, correspondant via deux conducteurs 5 électriques :c'est-à-dire, par exemple un conducteur de phase 32, , 34, ou 36, et le conducteur de neutre 38, Dans l'exemple de la figure 1, le consommateur électrique 161 est triphasé et les consommateurs électriques 162, 163 sont monophasés. Chaque consommateur électrique 16, inclut l'un des premiers capteurs 20, 10 correspondant et l'un des modules de transmission 21, correspondant. Chaque consommateur électrique 16, est, par exemple, un compteur de consommation électrique communiquant, propre à mesurer des premières données Ec, relatives à l'énergie électrique consommée par le consommateur électrique 16,, via le premier capteur 20, correspondant, et à transmettre les premières données Ec, à l'unité de 15 traitement 24 via le module de transmission 21, correspondant. Chaque départ 18, est un départ triphasé et comprend les trois conducteurs de phase 32,, 34,, 36, correspondants et le conducteur de neutre 38, correspondant. En variante, non représentée, les départs électriques sont monophasés et comprennent un conducteur de phase et un conducteur de neutre. Selon encore une autre variante, non 20 représentée, certains départs sont monophasés et d'autres sont triphasés. Chaque premier capteur 20, est propre à mesurer les premières données E(t1) relatives à l'énergie électrique consommée par le consommateur électrique 16, correspondant, pendant différents intervalles temporels tl. Plus généralement, les premiers capteurs 20, sont configurés pour mesurer les 25 premières données E(t1) pendant les mêmes intervalles temporels tl, les premières données Ec,(ti) mesurées au niveau de chaque consommateur électrique 16, étant mesurés de manière synchrone. Chaque module de transmission 21, est propre à transmettre les premières données Ec,(ti) mesurées par le premier capteur 20, correspondant à destination de l'unité 30 de traitement 24. Avantageusement, chaque module de transmission 21, est propre à transmettre avec les premières données Ec,(ti) une première information relative à l'intervalle temporel pendant lequel les premières données ont été mesurées. Chaque deuxième capteur 22, est propre à mesurer des deuxièmes données EDA) 35 relatives à l'énergie électrique délivrée par le départ 18, correspondant pendant les différents intervalles temporels tl. 3035496 9 Les premières E(t1) et deuxièmes ED,(0 données mesurées sont alors synchronisées au sens où elles sont mesurées pendant des intervalles temporels ti identiques. Chaque deuxième capteur 22, est également configuré pour transmettre, via une 5 liaison électrique respective 40,, les deuxièmes données ED,(0qu'il mesure à l'unité de traitement 24. Avantageusement, chaque deuxième capteur 22, est configuré pour transmettre avec les deuxièmes données ED,(ti) une deuxième information relative à l'intervalle temporel t, pendant lequel les deuxièmes données ED,(ti) ont été mesurées. 10 Les premières Ec(tl) et deuxièmes ED,(ti) données sont, par exemple, des mesures d'énergie active. Plus précisément, dans la suite de la description, on considère que les premières Ec(tl) et deuxièmes ED,(ti) données sont des mesures d'énergie active. En variante, les premières Ec,(t) et deuxièmes ED,(ti) données sont des mesures 15 d'énergie réactive, des mesures d'énergie apparente, des mesures de puissance active, des mesures de puissance réactive, des mesures de puissance apparente ou encore des mesures d'intensité. Le processeur 26 est configuré pour exécuter des logiciels compris dans la mémoire 28. 20 La mémoire 28 comprend un logiciel 41 d'acquisition des premières E(t1) et deuxièmes ED,(0 données, un logiciel 42 de génération de plusieurs jeux de données différents Je, à partir des premières données E(t1) et des deuxièmes données ED,(0 acquises dans un même intervalle temporel t, et un logiciel 44 de calcul d'un premier critère de sélection C11 pour chaque jeu de données Jel. 25 La mémoire 28 comprend également un logiciel 46 de sélection d'un premier ensemble En, de jeux de données Jel, en fonction des premiers critères de sélection C11 calculés par le logiciel de calcul 44 et un logiciel 48 de détermination de paramètres de raccordement au, lesdits paramètres de raccordement comportant pour chaque consommateur électrique 16, un identifiant du départ 18, auquel il est raccordé. 30 Les logiciels d'acquisition 41, de génération 42, de calcul 44, de sélection 46 et de détermination 48 correspondent à des instructions logicielles et forment un programme d'ordinateur apte à être exécuté par un ordinateur. L'ordinateur correspond, par exemple, à l'unité de traitement 24. Le logiciel d'acquisition 41 est propre, par exemple, à transmettre à chaque 35 consommateur électrique 16, et notamment à chaque premier capteur 20, un ordre de 3035496 10 mesure des premières données E(t1) et un ordre de transmission des premières données E(0, afin de récupérer les premières données. De même, le logiciel d'acquisition 41 est, par exemple, configuré pour transmettre à chaque départ 18, et notamment à chaque deuxième capteur 22, un ordre de mesure 5 des deuxièmes données ED,(0 et un ordre de transmission des deuxièmes données ED,(0, afin de récupérer les deuxièmes données ED,(0. Avantageusement les ordres de mesure des premières E(t1) et deuxièmes ED,(0 données sont transmis simultanément à tous les départs 18, et tous les consommateurs 16,. 10 Le logiciel de génération 42 est configuré pour générer les jeux de données Je, qui sont chacun associés à l'un des intervalles temporel t, et qui comprennent les premières Ec,(11) et deuxièmes ED,(ti) données associées audit intervalle temporel 11. En d'autres termes, le logiciel de génération 42 sélectionne les premières Ec,(11) et deuxièmes EDA) données mesurées pendant les différents intervalles temporels 11, pour créer les jeux de 15 données Jel. Le logiciel de calcul 44 est adapté pour calculer, pour chaque jeu de données Jel, le premier critère de sélection Cli, qui est choisi parmi un taux global de pertes d'énergie électrique entre les départs 18, et les consommateurs électriques 16, et un écart de consommation d'énergie électrique entre les différents consommateurs 16,. Les pertes 20 d'énergie électrique, incluent aussi bien les pertes dites techniques, par exemple liées aux pertes Joule lors de la circulation du courant entre les départs 18j et les consommateurs 16i, que les pertes dites non techniques, qui sont par exemple liées au vol d'électricité, au fait que des consommateurs sont raccordés au réseau de distribution 10 sans que l'unité de traitement d'informations 24 en soit informée et à des premiers capteurs 20, 25 défectueux. Si le premier critère de sélection C1, est, par exemple, le taux global de pertes d'énergie électrique entre les départs 18, et les consommateurs 16,, le premier critère de sélection Cl est calculé à partir de la formule suivante : Cl = 7AE (t1) = EDJ (ti ) i=1 (1) J=1 30 où I est un indice de jeu de données et varie de 1 à r, r étant le nombre de jeux données générés par le logiciel de génération 42, t, est l'intervalle temporel associé au jeu de données d'indice I, E(t1) est l'énergie électrique consommée par le consommateur 3035496 11 électrique d'indice i pendant l'intervalle temporel t, et ED (t) est l'énergie électrique délivrée par le départ d'indice j pendant l'intervalle temporel Si le premier critère de sélection C1, est, par exemple, l'écart de consommation d'énergie électrique entre les différents consommateurs 16,, le premier critère de sélection 5 Ci est calculé à partir de la formule suivante : Cli = Var(Eci(tk),- - - , En,(tk)) , (2) où Var est la fonction variance. En variante, si le premier critère de sélection Cl, est, par exemple, l'écart relatif de consommation d'énergie électrique entre les différents consommateurs, le premier critère 10 de sélection Cl est calculé à partir de la formule suivante : = Var(Eci(tk),- - - ,E(tk)) Cl, (3) Moy(Eci(tk),- - - ,E(tk))' où Moy est une fonction moyenne, telle qu'une moyenne arithmétique, géométrique, quadratique... Le logiciel de sélection 46 est configuré pour sélectionner le premier ensemble En, 15 de jeux de données parmi les jeux de données Je, générés par le logiciel de génération 42, en fonction des premiers critères de sélection C1, calculés. Le logiciel de sélection 46 est, par exemple, configuré pour comparer les premiers critères Cl avec une première variable prédéterminée V1 et pour sélectionner les jeux de données Je, pour lesquels le premier critère Ci est inférieur à la première variable prédéterminée Vl. 20 Le logiciel de détermination 48 est configuré pour déterminer les paramètres de raccordement au à partir du premier ensemble En, sélectionné. Le logiciel de détermination 48 est, par exemple, configuré pour établir ou déterminer un système d'équation à résoudre à partir d'un postulat de conservation de l'énergie électrique pour chaque départ 18J, selon lequel l'énergie délivrée par le départ 25 18, est sensiblement égale à la somme de l'énergie consommée par les consommateurs électriques 16, connectés au départ 18, et des pertes électriques. Le système d'équation est, par exemple, défini à partir de l'équation suivante : EDi(tk)=I(aii.E (tk)+ a1.'-a2;k) (4) i=1 où n est le nombre de consommateurs électriques, j est un indice de départ variant de 1 à 30 m, m étant le nombre de départs, i est un indice de consommateur électrique variant de 1 à n, 3035496 12 k est un indice de jeu de données Jek du premier ensemble En, variant de 1 à w, w étant le nombre de jeux de données Jek du premier ensemble Eni, tk est l'intervalle temporel associé au jeu de données Jek d'indice k, est le paramètre de raccordement indiquant si le consommateur d'indice i est raccordé 5 au départ d'indice j, Eak) est l'énergie électrique consommée par le consommateur électrique d'indice i pendant l'intervalle temporel tk, EDJ(tk) est l'énergie électrique délivrée par le départ d'indice j pendant l'intervalle temporel tk, et 10 1jk et a2 ,k sont des variables d'ajustement représentatives des pertes électriques, c'est à dire de l'écart pour chaque jeu de donnée Jek et pour chaque départ 18, entre l'énergie électrique délivrée par le départ 18, et l'énergie électrique consommée par les consommateurs électriques 16, raccordés au départ d'indice j. Le système d'équation s'écrit alors par exemple sous une forme matricielle de la 15 manière suivante : A*z= B, (5) E (t E (t w) E D2 (t et A = [D /w,', /w.m représentant la matrice E D2 (t w) EDm(ti) E Dm(' k) a11 an1 alm où z= anm J111 ,B= J 211 J _ mw_ U 0 o E1(t1) --- E(t1) o _E1(t) --- E(t)_ unité de taille w*m, avec D= - - . 0 , U= et D 0 --- 0 U comprenant m*n colonnes et m*w lignes. Le logiciel de détermination 48 est alors configuré pour résoudre le système 20 d'équation à partir d'un algorithme d'optimisation vérifiant : 3035496 {A.z=B min f .Z tel que /z z uz où fT .z =[0 --- 0 1 --- 1 1 --- 1]x z , f T correspondant à une fonction coût de m*n w*m w*m -e l'algorithme d'optimisation et z à un vecteur objectif défini ci-dessus et où, lz = 0 étant de dimension m*n+2*w*m et comprenant m*n fois la valeur -e et 2*w*m fois la 13 (6) -e 0 , lz 5 valeur 0 et uz = 1+e ,uz étant de dimension m*n+2*w*m et comprenant m*n fois la 1+e +00 +00 valeur (1+e) et 2*w*m fois la valeur +.0, avec ecompris entre 0 et 0,5, de préférence compris entre 0 et 0,1, de préférence encore égal à 0,05. L'algorithme d'optimisation est propre à déterminer le vecteur objectif z et donc les paramètres de raccordement pour lesquels la fonction fT.z est minimisée. Plus 10 précisément, le logiciel de détermination 48 est configuré pour démarrer à partir d'un a11 vecteur solution X = initial choisi aléatoirement et pour converger de façon itérative _anm _ vers une solution. A chaque itération, les variables d'ajustement i,k, a2,k et la fonction fT.z sont calculées et permettent de décider du prochain vecteur solution X. 3035496 14 Plus précisément, à chaque itération, l'algorithme d'optimisation équilibre les écarts d'énergie avec les variables d'ajustement a1,k,a2 jk de sorte que les égalités du système d'équation soient vérifiées. Ainsi, si l'énergie délivrée sur un départ d'indice j est excédentaire, on augmente la variable d'ajustement alik correspondante et si cette 5 énergie est déficitaire, on augmente la variable a2 jk correspondante. Ainsi, plus les paramètres de raccordement au respectent le principe de conservation de l'énergie, plus les variables d'ajustement sont faibles. L'objectif est de minimiser les variables d'ajustement, ce qui se traduit par la fonction objectif f T .zjk +a2 ) j=1 k=1 Les matrices /, et u, bornent les paramètres de raccordement au entre -e et 1+e et 10 les variables d'ajustement al jk et a2jk entre 0 et + CC) . Plus précisément, lors de l'application de l'algorithme d'optimisation, les paramètres de raccordement au sont des nombres réels ce qui permet de relaxer les contraintes. Puis, suite à l'application de l'algorithme d'optimisation, le logiciel de détermination 48 est configuré pour fixer les valeurs des paramètres de raccordement au à 0 ou à 1 en fonction de leur valeur suite à 15 l'application de l'algorithme d'optimisation. La valeur 0 indique un non raccordement du consommateur d'indice i au départ d'indice j, tandis que la valeur 1 indique un raccordement du consommateur d'indice i au départ d'indice j. Le logiciel de détermination 48 est, par exemple, configuré pour déterminer les valeurs des paramètres de raccordement suite à l'application de l'algorithme 20 d'optimisation via l'équation suivante : 1 si ati = max (ai.) ati = 0 sinon Ainsi, d'après l'équation (7) ci-dessus, chaque consommateur 16, est raccordé à un seul départ 18, 25 De plus, l'unité de traitement 24 est configuré pour identifier, en fonction des paramètres de raccordement déterminés et à partir, par exemple, d'un logiciel d'identification, non représenté, compris dans la mémoire 28, des sous-ensembles de consommateurs, avec chaque sous-ensemble de consommateurs qui correspond à l'ensemble des consommateurs 16, raccordés à un même départ 18,,. (7) 3035496 15 En variante, le logiciel de détermination 48 est configuré pour déterminer les paramètres de raccordement séparément pour chaque départ et indépendamment des deuxièmes données relatives à l'énergie électrique acquises pour les autres départs. 5 Selon cette variante, un système d'équation optimal est déterminé pour chaque départ et l'algorithme d'optimisation est appliqué à chaque système d'équation optimal. On obtient ainsi m systèmes d'équation optimaux résolus indépendamment via l'algorithme d'optimisation. Ainsi, pour le départ d'indice 1, le système d'équation optimal correspondant vérifie : 10 A1*z1= B1, a11 (8) anl J111 ED1(tl) où zl = B1= et Al [U avec a2iw _E D1 (tw )_ E1(t1) --- E(t1) U= _E1(t) --- E(t)_ Puis, selon cette variante, le logiciel de détermination 48 est, par exemple, propre à déterminer les valeurs des paramètres de raccordement via l'équation (7). 15 Selon une autre variante, le logiciel de détermination 48 est configuré pour déterminer les paramètres de raccordement pour chaque conducteur électrique de phase 32,, 34,, 36, et non simplement pour chaque départ 18,. Selon cette autre variante, le système d'équation comprend alors autant d'équation que de conducteurs électriques de phase 32,, 34,, 36, et les variables décrites ci-dessus et relatives à un départ spécifique j 20 sont alors relatives à un conducteur de phase spécifique. Ainsi, les deuxièmes capteurs 22, mesurent l'énergie électrique délivrée par chaque conducteur de phase 32,, 34,, 36, et non chaque départ 18,, les paramètres de raccordement au sont déterminés pour chaque conducteur de phase 32,, 34,, 36, et les variables d'ajustement sont déterminées pour chaque conducteur de phase 32,, 34,, 36,. Selon cette autre variante, on identifie par 25 exemple chaque conducteur électrique par un indice et la variable j, présentée dans les équations ci-dessus, correspond alors à un indice de conducteur électrique de phase et 3035496 16 variant de 1 à u, avec u le nombre de conducteur électriques de phase qui est égal à 3*m, soit 6 dans l'exemple de la figure 1. Plusieurs modes de réalisation d'un procédé de détermination de la structure du réseau de distribution électrique 10, mis en oeuvre via l'unité de traitement 24, et plus 5 globalement via le système de détermination 12, vont désormais être présentés. Selon un premier mode de réalisation décrit ci-après en référence à la figure 2, le procédé comprend une étape initiale 100 d'acquisition de données principales du réseau de distribution 10. Les données principales comprennent, par exemple, le nombre total n de consommateurs 16,, le nombre total m de départs 18J, les premières données E (t 10 mesurées, les deuxièmes données ED,(ti) mesurées et les différents intervalles temporels t, associés aux premières et deuxièmes données mesurées. En effet, lors de l'étape initiale 100, le logiciel d'acquisition 41 commande, par exemple, la mesure, par chaque premier 20, et chaque deuxième 22, capteurs, des premières E(t1) et deuxièmes Er)] (t1) données pendant les intervalles temporels t, et la transmission des premières 15 Ec(tl) et deuxièmes EDJ(t1) données qui sont alors associées à l'intervalle temporel t, pendant lequel elles ont été mesurées, via, par exemple, les premières et deuxièmes informations. Ensuite, lors d'une étape de génération 102, le logiciel de génération 42 génère plusieurs jeux de données Je, différents, associés chacun à l'un des intervalles temporels 20 t, et qui comprennent les premières E(t1) et deuxièmes ED,(ti) données associées audit intervalle temporel Puis, lors d'une étape de calcul 104, le logiciel de calcul 44 calcule le premier critère de sélection C1, pour chaque jeu de données Jel. Le premier critère de sélection Cl, est par exemple le taux global de pertes d'énergie électrique entre les départs 18, et 25 les consommateurs 16,. Ensuite, au cours d'une étape de sélection 106, le logiciel de sélection 46 sélectionne le premier ensemble En, en fonction des premiers critères de sélection C1, calculés. Le premier ensemble En, est sélectionné parmi les jeux de données générés à l'étape 102. Au cours de l'étape de sélection 106, le logiciel de sélection 46 compare, par 30 exemple, les premiers critères Cl, avec la première variable prédéterminée V1, la valeur de la première variable prédéterminée V1 étant par exemple définie lors de l'étape d'acquisition 100. Le logiciel de sélection 46 sélectionne alors les jeux de données Je, pour lesquels le premier critère Ci est inférieur à la première variable prédéterminée Vl. 3035496 17 Enfin, au cours d'une étape de détermination 108, les paramètres de raccordement sont déterminés via le logiciel de détermination 48 et à partir du premier ensemble En, sélectionné. Plus précisément, le logiciel de détermination 48 détermine le système d'équation 5 à résoudre à partir, par exemple, de l'équation (4), comme présenté ci-dessus. Puis, comme présenté ci-dessus lors de la description du système de détermination 12 et de l'équation (5), l'algorithme d'optimisation est appliqué au système d'équation afin de déterminer les paramètres de raccordement au. En variante, lors de l'étape de détermination, le logiciel de détermination 48 10 détermine les systèmes d'équations optimaux, comme présenté ci-dessus via l'équation (8), et applique l'algorithme d'optimisation à chaque système d'équation optimal. La vitesse de détermination des paramètres de raccordement au est alors améliorée, puisque les systèmes d'équation optimaux comprennent un nombre limité d'équations. Dans le premier mode de réalisation, le fait de réaliser une sélection des jeux de 15 données Jek permet d'éliminer les jeux de données Je, pour lesquels les pertes sont les plus importantes, ces jeux de données risquant de conduire à une détermination erronée des paramètres de raccordement au. Ainsi, selon le premier mode de réalisation, les paramètres de raccordement déterminés au, sont déterminés de manière plus certaine par rapport aux méthodes 20 connues de l'état de la technique et la fiabilité du procédé de détermination est améliorée. En outre, l'algorithme d'optimisation utilisé permet d'utiliser un nombre indifférent de jeux de données Jek lors de l'étape de détermination 108, même s'il est préférable que le nombre de jeux de données du premier ensemble soit supérieur ou égal au nombre total n de consommateurs 16,. 25 Selon un deuxième mode de réalisation de l'invention décrit ci-après en référence à la figure 3, le procédé comprend des étapes 200, 202, 204, 206, 208 identiques respectivement aux étapes 100, 102, 104, 106, 108 du premier mode de réalisation et, lors de l'étape 204, un deuxième critère de sélection C2, est calculé pour chaque jeu de données généré. Le deuxième critère de sélection C2, est différent du premier critère Cl, 30 et est choisi parmi un taux global de pertes d'énergie électrique entre le ou les départs et les consommateurs et un écart de consommation d'énergie électrique entre les différents consommateurs. Ainsi, lors de l'étape 204, chaque premier critère Cl est, par exemple, le taux global de perte relatif au jeu de données Je, correspondant et chaque deuxième critère 35 C2, est l'écart de consommation d'énergie électrique relatif au jeu de données Je, correspondant, et est calculé via l'équation (2) ou l'équation (3). 3035496 18 Puis, au cours de l'étape de sélection 206, le logiciel de sélection 46 sélectionne le premier ensemble En, en fonction des premiers critères de sélection C1, et des deuxièmes critères de sélection C2, calculés. Le premier ensemble En, est sélectionné parmi les jeux de données générés à l'étape 202. Au cours de l'étape de sélection 206, le 5 logiciel de sélection 46 compare, par exemple, les premiers critères C1, avec la première variable V1 et les deuxièmes critères C2, avec une deuxième variable prédéterminée V2. La valeur de la deuxième variable prédéterminée V2 est, par exemple, définie lors de l'étape d'acquisition 200. Le logiciel de sélection 46 sélectionne alors les jeux de données pour lesquels le premier critère C1, est inférieur à la première variable prédéterminée V1 10 et le deuxième critère C2, est inférieur à la deuxième variable prédéterminée V2. Enfin, au cours de l'étape de détermination 208 les paramètres de raccordement sont déterminés via le logiciel de détermination 48 et à partir du premier ensemble sélectionné Enl. Le deuxième mode de réalisation permet d'affiner la sélection des jeux de 15 données par rapport au premier mode de réalisation et ainsi de sélectionner des jeux de données à partir desquels le risque d'erreur dans la détermination des paramètres de raccordement est minimisé. La justesse, la vitesse et la fiabilité du procédé de détermination sont ainsi améliorées. Selon un troisième mode de réalisation de l'invention décrit ci-après en référence à 20 la figure 4, le procédé comprend des étapes 300, 302, 304, 306, 310 identiques respectivement aux étapes 100, 102, 104, 106, 108 du premier mode de réalisation et, suite à l'étape 306 et précédemment à l'étape 310, le procédé comprend une étape 307 de calcul d'un deuxième critère de sélection C2k pour chaque jeu de données du premier ensemble En,, le deuxième critère C2k étant différent du premier critère et étant choisi 25 parmi un taux global de pertes d'énergie électrique entre le ou les départs 18, et les consommateurs 16, et un écart de consommation d'énergie électrique entre les différents consommateurs 16,. Dans le troisième mode de réalisation, chaque premier critère C1, est, par exemple, le taux global de perte relatif au jeu de données correspondant généré à l'étape 30 302 et chaque deuxième critère C2k est l'écart de consommation d'énergie électrique relatif au jeu de données correspondant du premier ensemble Enl. Suite à l'étape de calcul 307, lors d'une étape de sélection 308, un deuxième ensemble En2 de jeux de données est sélectionné parmi le premier ensemble Enl. Puis, le premier ensemble En, est, par exemple, fixé égal au deuxième ensemble sélectionné 35 En2 pour réaliser l'étape de détermination 310 en fonction du deuxième ensemble sélectionné En2. 3035496 19 Le troisième mode de réalisation permet d'affiner la sélection des jeux de données par rapport au premier mode de réalisation, et ainsi de sélectionner des jeux de données à partir desquels le risque d'erreur dans la détermination des paramètres de raccordement au est minimisé. La justesse, la vitesse et la fiabilité du procédé de 5 détermination sont ainsi améliorées. Selon un quatrième mode de réalisation de l'invention décrit ci-après en référence à la figure 5, le procédé comprend des étapes 400, 402, 404, 406 identiques aux étapes 100, 102, 104, 106 du premier mode de réalisation. Puis, suite à l'étape de sélection 406, le procédé comprend une étape 408 de 10 détermination des paramètres de raccordement. Plus précisément, lors de l'étape 408, le procédé comprend une première sous-étape 408A de sélection pseudo-aléatoire d'un troisième ensemble En3 de jeux de données inclut dans le premier ensemble Eni. Ensuite, lors d'une sous-étape de calcul 408B, un paramètre de répétition R1 est calculé. Le paramètre de répétition R1 correspond à un nombre d'itérations de la sous- 15 étape de sélection 408A. Puis, au cours d'une sous-étape de détermination 4080, des indices intermédiaires b,, de raccordement, indiquant pour chaque consommateur 16, le départ 18, auquel il est raccordé, sont déterminés à partir du troisième ensemble sélectionné En3. Plus généralement, les indices intermédiaires b,, comportent pour chaque consommateur 20 électrique 16, un identifiant du départ 18j auquel il est raccordé. La sous-étape de détermination 4080 est analogue à l'étape de détermination 108 du premier mode de réalisation, mais est réalisée à partir du troisième ensemble En3. Ensuite, lors d'une sous-étape de mémorisation 408D, les indices intermédiaires de raccordement b,, sont mémorisés par la mémoire 28. 25 Puis, au cours d'une sous-étape de calcul 408E, réalisée pour chaque consommateur 16, et pour chaque départ 18,, un nombre d'affectation NA,, relatif au nombre de fois où les indices intermédiaires de raccordement mémorisés indiquent que le consommateur 16, est raccordé au départ 18, est calculé. Le nombre d'affectation NA,, est par exemple calculé via la formule suivante : N(bu =1) 30 NAu = , où N(bij 1) correspond au nombre de fois où les indices NT intermédiaires b,, mémorisés sont égaux à 1 pour le consommateur 16, d'indice i et le départ 18, d'indice j et NT correspond à un nombre total d'itérations de la sous-étape de sélection 408A. Avantageusement, le nombre total d'itérations NT est initialisé à 0 3035496 20 précédemment à l'exécution de la sous-étape de sélection 408A et incrémenté de 1 à chaque exécution de la sous-étape de sélection. Il est à noter que le nombre d'indices intermédiaires de raccordement mémorisés pour un consommateur électrique 16, donné et un départ 18, donné est égal au nombre 5 d'itérations de la sous-étape de sélection 408A Avantageusement, les nombres d'affectation calculés sont mémorisés suite à la sous-étape de calcul 408E. Puis, lors d'une sous-étape de comparaison 408F, le paramètre de répétition R-1 est comparé à un critère de répétition prédéterminé CR1. Le critère de répétition 10 prédéterminé CR1 est, par exemple, initialisé lors de l'étape d'acquisition 400. Le critère de répétition CR1 est, par exemple, un nombre minimal d'itérations de la sous-étape de sélection 408A. Si lors de la sous-étape de comparaison 408F, le paramètre de répétition est inférieur au critère de répétition, alors les étapes de sélection 408A, de calcul 408B, de 15 détermination 4080, de mémorisation 408D, de calcul 408E et de comparaison 408F sont répétées. Si lors de la sous-étape de comparaison 408F, le paramètre de répétition R-1 est supérieur au critère de répétition CR1, alors une sous-étape 408G de détermination des paramètres de raccordement a,, est effectuée. 20 Lors de la sous-étape de détermination 408G, les paramètres de raccordement sont déterminés à partir du ou des troisièmes ensembles sélectionnés En3, et plus précisément en fonction des indices intermédiaires de raccordement b,, mémorisés lors de la sous-étape de mémorisation 408D, et plus précisément encore en fonction des nombres d'affectation NA,, calculés à la dernière itération de la sous-étape 408E. Plus 25 précisément, pour chaque consommateur électrique 16,, l'indice de départ j correspondant au nombre d'affectation NA,, le plus grand est identifié, et le paramètre de raccordement au correspondant audit consommateur 16, et audit départ 18, est fixé égal à 1, les autres paramètres de raccordement au relatifs audit consommateur étant fixés égaux à 0. Avantageusement, si lors de la sous-étape de détermination 408G tous les 30 nombres d'affectation relatifs à un consommateur 16, sont inférieurs à un premier seuil prédéterminé Si, par exemple égal à 0,6, alors une erreur d'identification du raccordement pour ledit consommateur 16, est identifiée. En variante, si le nombre d'itérations de la sous-étape 408A est supérieur à 2, alors, lors de la sous-étape de calcul 408E, et à chaque itération de la sous-étape 408E, 35 une première moyenne des nombres d'affectation NA, calculés à chaque itération de la sous-étape 408E, pour chaque consommateur 16, et pour chaque départ 18,, est calculée. 3035496 21 Puis, une deuxième moyenne des nombres d'affectation NA, calculés à la dernière itération, pour chaque consommateur 16, et pour chaque départ 18,' est calculée. Ensuite, lors de la sous-étape 408E, une différence entre la première moyenne et la deuxième moyenne est calculée. 5 La première et la deuxième moyenne sont, par exemple, des moyennes arithmétiques, géométriques ou quadratiques... Selon cette variante, lors de la sous-étape de comparaison 408F, la valeur absolue de la dernière différence calculée est comparée à un deuxième seuil prédéterminé S2, par exemple égal à 0,1. Puis, si lors de la sous-étape de comparaison 10 408F, le paramètre de répétition R-1 est supérieur au critère de répétition CR1 et la valeur absolue de la dernière différence calculée est inférieure au deuxième seuil S2, alors la sous-étape 408G est effectuée. Sinon, les étapes 408A, 408B, 4080, 408D, 408E et 408F sont réitérées. Le quatrième mode de réalisation permet, notamment lorsque le nombre 15 d'itérations de la sous-étape 408A est supérieur à 2, de déterminer les paramètres de raccordement au à partir de troisièmes ensembles En3 de jeux de données différents. Ainsi, la justesse et la fiabilité du procédé de détermination sont améliorées. En outre, le quatrième mode de réalisation permet avantageusement d'identifier chaque consommateur 16, pour lequel le départ associé 18j est déterminé avec un bon 20 indice de confiance et chaque consommateur 16, pour lequel le départ associé 18j est indéterminé ou déterminé avec un mauvais indice de confiance. En effet, comme présenté ci-dessus, si lors de la sous-étape de détermination 408G tous les nombres d'affectation relatifs à un consommateur 16, sont inférieurs au premier seuil prédéterminé Si, par exemple égal à 0,6, alors une erreur d'identification du raccordement pour ledit 25 consommateur 16, est identifiée et le départ 18j auquel est raccordé le consommateur 16, est indéterminé. De la même manière, si lors de la sous-étape de détermination 408G un nombre d'affectation relatif à un consommateur 16, est supérieur à un troisième seuil prédéterminé S3, par exemple égal à 0,95, alors le raccordement du consommateur 16, au départ 18j correspondant est identifié avec un bon indice de confiance. 30 Avantageusement encore, le quatrième mode de réalisation permet d'associer, à chaque paramètre de raccordement au, fixé égal à 1 lors de la sous-étape de détermination 408G, un indice de confiance représentant la probabilité que le raccordement déterminé soit correct. L'indice de confiance est par exemple égal au nombre d'affectation NA,, correspondant. 3035496 22 Les modes de réalisations et variantes envisagées ci-dessus sont propres à être combinées entre eux, selon toutes combinaisons techniquement admissibles, pour donner lieu à d'autres modes de réalisation de l'invention. Ainsi, le deuxième mode de réalisation est propre à être combiné avec le 5 quatrième mode de réalisation et le troisième mode de réalisation est également propre à être combiné avec le quatrième mode de réalisation. REVENDICATIONS1.- Procédé de détermination de la structure d'un réseau (10) de distribution d'électricité, le procédé étant mis en oeuvre par ordinateur (24), le réseau de distribution comprenant un poste d'alimentation (14) comportant un ou plusieurs départs (18J) d'alimentation en énergie électrique de plusieurs consommateurs électriques (16,) et comprenant les étapes suivantes : -a) l'acquisition (100; 200 ; 300; 400), pour chaque consommateur électrique (16,), de premières données (Ec,(t1)) relatives à l'énergie électrique consommée par le consommateur électrique (16,) pendant différents intervalles temporels (t1), mesurées via un premier capteur (20,) de mesure d'énergie, -b) l'acquisition (100 ; 200 ; 300 ; 400), pour chaque départ (18J), de deuxièmes données (EDJ(t1)) relatives à l'énergie électrique délivrée par le départ (18J) pendant les différents intervalles temporels (t1), mesurées via un deuxième capteur (20J) de mesure d'énergie, - c) la génération (102; 202 ; 302; 402), via un module de génération (42), de plusieurs jeux de données (JO différents, associés chacun à l'un des intervalles temporels, chaque jeu de données (Je') comprenant les premières (E(0) et deuxièmes (ED,(0) données associées audit intervalle temporel (t1), caractérisé en ce que le procédé comprend en outre les étapes suivantes : -d) le calcul (104; 204; 304; 404) d'un premier critère de sélection (C11) pour chaque jeu de données, le premier critère de sélection (Cl') étant un critère choisi parmi un taux global de pertes d'énergie électrique entre le ou les départs (18J) et les consommateurs (16,) et un écart de consommation d'énergie électrique entre les différents consommateurs (16,), -e) la sélection (106 ; 206 ; 306 ; 406), par une unité de sélection (46) et en fonction des premiers critères de sélection (Cl,) calculés, d'un premier ensemble (Eni) de jeux de données (Jek) parmi les jeux générés (JO, -f) la détermination (108; 208; 310; 408), à partir du premier ensemble (En1) sélectionné, de paramètres de raccordement (a), lesdits paramètres (a,J) comportant pour chaque consommateur électrique (16,) un identifiant du départ (18J) auquel il est raccordé. 2.- Procédé selon la revendication 1, dans lequel, lors de l'étape de calcul d) (204), un deuxième critère de sélection (C21) pour chaque jeu de données (JO généré est calculé, le deuxième critère (C21) étant l'autre critère choisi parmi un taux global de pertes 3035496 24 d'énergie électrique entre le ou les départs (18J) et les consommateurs (16,) et un écart de consommation d'énergie électrique entre les différents consommateurs (16,), et dans lequel, lors de l'étape de sélection e) (206), le premier ensemble (En1) est sélectionné en fonction des premiers (Cl') et deuxièmes (021) critères calculés. 5 3.- Procédé selon la revendication 1, dans lequel, suite à l'étape de sélection e), le procédé comprend les étapes suivantes : -el) le calcul (307) d'un deuxième critère de sélection (021) pour chaque jeu de données (Jek) du premier ensemble (En1), le deuxième critère (021) étant l'autre critère 10 choisi parmi un taux global de pertes d'énergie électrique entre le ou les départs (18J) et les consommateurs (16,) et un écart de consommation d'énergie électrique entre les différents consommateurs (16,), -e2) la sélection (308) en fonction des deuxièmes critères de sélection (021) calculés d'un deuxième ensemble (En2) de jeux de données parmi le premier ensemble 15 (En1), la détermination des paramètres de raccordement (a,J) étant alors effectuée lors de l'étape de détermination f) (310) en fonction du deuxième ensemble (En2) sélectionné parmi le premier ensemble (En1). 20 4.- Procédé selon l'une des revendications précédentes, dans lequel, lors de l'étape de détermination f) (108 ; 208 ; 310 ; 408), les paramètres de raccordement (a,J) sont déterminés séparément pour chaque départ (18J) et indépendamment des deuxièmes données (ED,(0) relatives à l'énergie électrique acquises pour les autres départs (18J). 25 5.- Procédé selon l'une des revendications précédentes, dans lequel, lors de l'étape de détermination f) (108 ; 208; 310 ; 408), un système d'équation à résoudre est déterminé à partir d'un postulat de conservation de l'énergie électrique pour chaque départ (18J) selon lequel l'énergie délivrée par le départ (18J) est sensiblement égale à la 30 somme de l'énergie consommée par les consommateurs électriques (16,) connectés au départ (18J) et des pertes électriques. 6.- Procédé selon la revendication 5, dans lequel, lors de l'étape de détermination f) (108 ; 208 ; 310 ; 408), le système d'équation est défini à partir de l'équation suivante : 35 EDi(tk)= (aii .Eci (tk )± ai -a2 ) , i=1 3035496 25 où n est le nombre de consommateurs électriques (16,), j est un indice de départ (18,) variant de 1 à m, m étant le nombre de départs (18,), i est un indice de consommateur électrique variant de 1 à n, k est un indice de jeu de données (Jek) du premier ensemble (Eni)variant de 1 à w, w 5 étant le nombre de jeux de données (Jek) du premier ensemble (En1), tk est l'intervalle temporel associé au jeu de données (Jek) d'indice k, est le paramètre de raccordement indiquant si le consommateur (16,) d'indice i est raccordé au départ (18J) d'indice j, Eak) est l'énergie électrique consommée par le consommateur électrique (16,) d'indice i 10 pendant l'intervalle temporel tk, E Di (tk) est l'énergie électrique délivrée par le départ (18,) d'indice j pendant l'intervalle temporel tk, et ai et a2,k sont des variables d'ajustement représentatives des pertes électriques. 15 7.- Procédé selon la revendication 6, dans lequel le système d'équation s'écrit sous une forme matricielle de la manière suivante : A*z= B, où z= a11 B= Em(ti) et A= U o 0 -.. o avec an1 EDI (tw ED2(t1) 0 - - 0 U alm ED2 (t w) E Dm(t a nm E Dm(t k) aimw a2ii a2 mw_ U= E1(t1) --- E(t1) et A comprenant m*n colonnes et m*w lignes. _E1(t) --- E(t)_ 20 8.- Procédé selon la revendication 7, dans lequel, lors l'étape de détermination, un algorithme d'optimisation est utilisé pour résoudre le système d'équation, l'algorithme d'optimisation vérifiant : 3035496 26 {A.z = B min f T .Z tel que , oùfT.z = [0 --- 0 1 --- 1 1 --- 1]xz, /z z ,i*' w*m w*m , avec e compris entre 0 et 0,5. 5 9.- Procédé selon l'une des revendications précédentes, dans lequel chaque départ (18J) comprend un ou plusieurs conducteurs électriques de phase, et dans lequel, lors de l'étape d'acquisition b), les deuxièmes données (EDA)) sont acquises pour chaque conducteur électrique de phase et sont relatives à l'énergie électrique délivrée par chaque conducteur électrique de phase, et, lors de l'étape de détermination f) (108 ; 208 ; 310 ; 10 408), les paramètres de raccordement (a,J) déterminés comportent pour chaque consommateur électrique (16,) un identifiant du ou des conducteurs électriques de phase auxquels il est raccordé. 10.- Procédé selon l'une des revendications précédentes, dans lequel, lors de 15 l'étape de détermination f) (408), le procédé comprend les étapes suivantes : - f1) la sélection (408A) pseudo-aléatoire d'un troisième ensemble (En3) de jeux de données inclus dans le premier ensemble (Ehl), - f2) la détermination (408G) des paramètres de raccordement (a,J) à partir du troisième ensemble (En3) sélectionné parmi le premier ensemble (Eni)- 20 11.- Procédé selon la revendication10, dans lequel, suite à l'étape f1) de sélection (408A) et précédemment à l'étape de détermination f2) (408G), le procédé comprend les étapes suivantes : -f3) le calcul (408B) d'un paramètre de répétition (R1), 25 -f4) la détermination (4080), à partir du troisième ensemble (En3) sélectionné, d'indices intermédiaires de raccordement (b,J) comportant pour chaque consommateur électrique (16,) un identifiant du départ (18J) auquel il est raccordé, -f5) la mémorisation (408D) des indices intermédiaires de raccordement (b), et -f6) la comparaison (408F) du paramètre de répétition avec au moins un critère de 30 répétition prédéterminés (CR1), lz -e et uz -e 0 3035496 27 et, suivant le résultat de l'étape de comparaison f6) (408F), la réalisation d'une étape parmi, d'une part, la répétition des étapes de sélection f1), de calcul f3), de détermination f4), de mémorisation f5) et de comparaison f6) et, d'autre part, l'étape de détermination f2) des paramètres de raccordement (a,J) en fonction des indices intermédiaires de 5 raccordement (b,J) mémorisés. 12.- Procédé selon la revendication 11, dans lequel, le procédé comprend l'étape suivante : -f7) le calcul (408E), pour chaque consommateur (16,) et pour chaque départ (18J), 10 d'un nombre d'affectation (NA) relatif au nombre de fois où les indices intermédiaires de raccordement (b,J) mémorisés indiquent que le consommateur (16,) est raccordé au départ (18J), et dans lequel, lors de l'étape de détermination f2) (408G), les paramètres de raccordement (a,J) sont déterminés à partir des nombres d'affectation (NA) calculés. 15 13.- Programme d'ordinateur comportant des instructions logicielles, qui lorsqu'elles sont exécutées par un ordinateur mettent en oeuvre le procédé selon l'une quelconque des revendications précédentes. 20 14.- Système (12) de détermination de la structure d'un réseau (10) de distribution d'électricité, le réseau de distribution (10) comprenant un poste d'alimentation comportant un ou plusieurs départs (18J) d'alimentation en énergie électrique de plusieurs consommateurs électriques (16,), le système comprenant : - pour chaque consommateur électrique (16,), un premier capteur (20,) de mesure 25 de premières données (E(0) relatives à l'énergie électrique consommée par le consommateur (16,) pendant différents intervalles temporels (0, - pour chaque départ (18J), un deuxième capteur (22J) de mesure de deuxièmes données (ED,(0) relatives à l'énergie électrique délivrée par le départ (18J) pendant les différents intervalles temporels (0, 30 - un module (41) d'acquisition des premières (Ec,(0) et deuxièmes (ED,M) données, - un module (42) de génération de plusieurs jeux de données (Je') différents, associés chacun à l'un des intervalles temporels (0, chaque jeux de données (Je') comprenant les premières (Ec,(0) et deuxièmes (Ec,(0) données associées audit 35 intervalle temporel (0, caractérisé en ce que le système de détermination (12) comprend en outre : 3035496 28 - un module (44) de calcul d'un premier critère de sélection (C11) pour chaque jeu (Je') de données, le premier critère de sélection (Cl') étant un critère choisi parmi un taux global de pertes d'énergie électrique entre le ou les départs (18,) et les consommateurs (16,) et un écart de consommation d'énergie électrique entre les différents 5 consommateurs (16,), - un module (46) de sélection, en fonction des premiers critères de sélection (Cl,) calculés, d'un premier ensemble (Eni) de jeux de données (Jek) parmi les jeux générés, et - un module (48) détermination, à partir du premier ensemble (Eni) sélectionné, de paramètres de raccordement (a), lesdits paramètres (au) comportant pour chaque 10 consommateur électrique (16,) un identifiant du départ (18J) auquel il est raccordé.