Brevet FR3084168A1 — publié le 20200124

Procédé et système de détermination d'un état en tension d'un réseau de distribution basse tension au niveau d'au moins un poste de transformation [001] L'invention concerne un procédé d'évaluation de la tension dans la branche basse tension d'un réseau de distribution d'électricité, afin d'y réguler la tension pour la maintenir dans les plages de tension réglementaires, notamment lorsque des générateurs d'énergies renouvelables, induisant des fluctuations de tension, y sont connectés. [002] La gestion d'un réseau de distribution a pour objet de piloter de manière optimale son fonctionnement afin de minimiser les pertes, maîtriser les fluctuations de tension, et garantir ainsi une continuité de service dans des conditions économiques raisonnables. Or, la gestion du réseau électrique de distribution est amenée à évoluer profondément en raison de nombreux facteurs. Ces facteurs comprennent notamment le développement de nouveaux moyens décentralisés de production d'électricité dé-carbonée, fondée sur l'utilisation de sources d'énergies primaires renouvelables, l'évolution des usages de l'électricité, par exemple avec l'apparition des véhicules électriques, l'évolution des modes de consommation, par exemple en différant certaines consommations pour réduire la facture énergétique ou encore l'apparition de nouveaux acteurs sur le marché de l'énergie. [3] Les moyens de production décentralisée, ou « générateurs d'énergie dispersée GED », peuvent être gérés par des producteurs indépendants. Les moyens de production décentralisée se caractérisent par le fait qu'ils sont de puissance limitée et souvent raccordés au réseau de distribution. En se situant au plus près des points de consommation, ces moyens permettent de diminuer le coût de transport de l'énergie, les transits de puissance active au niveau du poste source et ainsi les pertes en ligne sur le réseau de transport. Toutefois, la connexion des moyens de production décentralisée sur le réseau de distribution pose de nouveaux problèmes d'exploitation, notamment la difficulté de tenir le plan de tension, c'est-à-dire les valeurs de tension dans les plages réglementaires. [4] Actuellement, le niveau de tension dans un réseau de distribution est fonction des quantités d'injection et de soutirage ayant lieu sur ce réseau. Autrement dit, le niveau de tension est fonction de la production et de la consommation en tension sur ce réseau. Une injection a pour effet d'augmenter la tension alors qu'un soutirage a pour effet de la faire baisser Afin de respecter les normes sur la qualité de tension, la tension distribuée doit être maintenue entre deux niveaux de prescription, définis contractuellement. Chaque réseau 2 haute-tension et basse tension dispose de ses propres plages. Le maintien de la tension dans les plages permises permet, d'assurer la sécurité des personnes et des biens, d'assurer le bon fonctionnement des différents appareils connectés au réseau et de limiter les pertes par effet 5 Joule dans le réseau. [5] Le réseau de distribution est un réseau arborescent dont les caractéristiques des lignes sont connues. Pour y réguler la tension, le gestionnaire de réseau actionne des dispositifs de contrôle-commande tels que le régleur en charge au poste HTB/HTA qui permettent d'ajuster 10 le rapport de transformation haute tension/basse tension entre le réseau de transport et le réseau de distribution. Tant que les puissances injectées par les générateurs d'énergie dispersés sur ce réseau sont peu nombreuses et faibles, les perturbations induites sur la tension aux différents noeuds du réseau sont faibles, localisées, et peuvent être compensées par cette technique de régulation de la tension. 15 [6] Dans un contexte de transition énergétique, le fort développement de nouveaux moyens décentralisés de production d'électricité d'origine renouvelable, utilisés par un très grand nombre d'acteurs va entrainer, sur le réseau basse tension, un grand nombre d'injections d'énergie renouvelable, induisant ainsi des fluctuations de la tension sur le réseau dépassant 20 les seuils légaux, ne pouvant être traitées avec les techniques de régulation actuelles dans des conditions économiques satisfaisantes. [7] La technique actuelle de régulation de la tension sur le réseau de distribution se fait par une mesure locale de la tension, une consigne fixe de tension et un ajustement du rapport de 25 transformation au niveau du poste-source HTB/HTA. [8] Le brevet français FR 1355392 vise à perfectionner cette technique par détermination de l'état en tension du réseau HTA, en utilisant les équations de Kirchhoff, en aval du poste-source et une consigne variable en lieu et place de la consigne fixe. 30 [9] Le procédé proposé vise à proposer une régulation de tension au niveau HTA/BT en améliorant la détermination de l'état en tension du réseau au moyen d'une procédure de résolution plus précise et plus stable et fiable des équations de Kirchhoff. 35 [010] La présente invention vient améliorer la situation. 3 [011] A cet effet, elle propose un procédé mis en oeuvre par des moyens informatiques, de détermination d'un état en tension d'un réseau de distribution basse tension au niveau d'au moins un poste de transformation, la branche basse tension comprenant une pluralité de 5 noeuds, le procédé comprenant: - une définition d'une valeur d'une tension cible réelle pour chaque noeud ; - une estimation d'une puissance active et d'une puissance réactive pour chaque noeud ; - une construction d'un premier système d'équations liant des tensions complexes aux puissances actives et réactives de la pluralité de noeuds, 10 - une formulation respective desdites puissances actives et réactives sous la forme d'une fonction quadratique pour chaque noeud, lesdites fonctions quadratiques étant substituées dans le premier système d'équations de manière à obtenir un deuxième système d'équations en nombre réels, équivalent au premier système d'équations, - détermination d'une fonction auxiliaire issue du deuxième système d'équations, 15 le procédé comprenant: - une minimisation itérative de ladite fonction auxiliaire, dont une solution initiale à la première itération est la tension cible, jusqu'à obtenir une solution au voisinage d'un minimum global approché de ladite fonction auxiliaire, - et une application au deuxième système d'équations d'une méthode itérative par 20 linéarisations successives, ayant comme solution initiale à la première itération ledit minimum global approché, jusqu'à obtenir une valeur d'une tension complexe en chaque noeud. [12] De cette manière, le procédé permet une précision de résultat très importante avec peu 25 de ressources en calcul. [13] Selon une réalisation, le procédé comprend en outre : - un calcul d'une variation entre la valeur de ladite tension cible réelle et la valeur d'un module de ladite tension complexe en chaque noeud; et 30 - si ladite variation est inférieure à une valeur absolue d'un prem er seuil, aucune action de régulation n'est entreprise, r si ladite variation en valeur absolue est comprise entre le premier seuil et le deuxième seuil, la valeur absolue du deuxième seuil étant supérieure à la valeur absolue du premier seuil, une détermination de consignes de tension à appliquer au 35 poste de transformation du réseau de manière à ramener la tension complexe en chaque noeud de la branche basse tension entre des valeurs de contrainte de tension 4 maximale et minimale ; et une transmission des consignes de tension au poste de transformation, - si ladite variation en valeur absolue est supérieure au deuxième seuil, le procédé 5 est relancé. [14] Ainsi, le procédé est relancé lorsqu'une valeur aberrante est obtenue en sortie. Le procédé est alors très efficace puisque les valeurs aberrantes ne sont pas prises en compte. De plus, une nouvelle valeur de consigne est déterminée à chaque cycle ce qui permet de 10 contrôler l'état en tension du réseau. [15] Selon une réalisation, la valeur absolue du premier seuil est égale à 10% de la valeur de la tension cible réelle, et la valeur absolue du deuxième seuil est égale à 15% de la valeur de la tension cible réelle. 15 [16] Cela permet de définir les conditions de fonctionnement normal du réseau. En dehors de ces seuils, la solution obtenue est déterminée comme étant aberrante ou une situation d'effondrement de tension est détectée. 20 [017] Selon une réalisation, les puissances actives et réactives en chaque noeud de la branche basse tension sont déterminées : - en fonction de mesures réelles fournies par au moins un compteur électrique relié au noeud associé de la branche basse tension, ou - en fonction d'un modèle de charge de consommateurs. 25 [18] Ainsi, les puissances actives et réactives sont obtenues de manière précise, en chaque noeud et en prenant en compte des différences selon les profils des utilisateurs. Notamment, ces différences comprennent les contrats propres à chaque utilisateur, l'utilisation de générateurs d' én ergi e... 30 [19] Selon une réalisation, les étapes du procédé sont répétées à chaque intervalle de temps, de préférence inférieur à trente minutes, et de préférence inférieur à dix minutes. [20] Ainsi, cela permet de pouvoir relancer le procédé dans un intervalle de temps 35 suffisamment court pour que les valeurs de consigne soient efficaces. Autrement dit, de sorte que la tension obtenue en sortie ne soit pas devenue obsolète. 5 [21] Selon une réalisation, la fonction auxiliaire est choisie de sorte qu'un minimum global de valeur 0 de ladite fonction auxiliaire est une solution du deuxième système d'équations, 5 ladite fonction auxiliaire étant à valeur réelle positive ou nulle. [22] Ainsi, le procédé permet de prendre en considération la non-linéarité des systèmes d'équations. 10 [023] Selon une réalisation, une résolution du deuxième système d'équations comprend la minimisation de la fonction auxiliaire par une méthode itérative comprenant, à chaque itération, une recherche d'une direction de descente puis une minimisation unidimensionnelle de la fonction auxiliaire jusqu'à obtenir la solution au voisinage du minimum global. 15 [024] Selon une réalisation, la direction de descente est trouvée, à chaque itération, par une résolution d'un problème d'optimisation à objectif quadratique convexe sous contraintes linéaires. [25] Selon une réalisation, la méthode itérative de linéarisations successives comprend, à 20 chaque itération : une résolution d'une approximation linéaire du deuxième système d'équations au voisinage d'une solution obtenue à l'itération précédente, laquelle définit une direction de déplacement dans l'espace des solutions, et une minimisation unidimensionnelle de la fonction auxiliaire dans cette direction 25 de déplacement. [26] Selon une réalisation, la valeur de la tension complexe est obtenue avec un degré de précision compris entre 10-6 et 1015 30 [027] Le procédé permet donc d'obtenir des valeurs de la tension complexe extrêmement précises. [028] Selon une réalisation, les nombres d'itérations pour la résolution de la fonction auxiliaire et d'itérations de la méthode de linéarisations successives sont chacun inférieurs à 35 20, de préférence inférieurs à 15, et de préférence égaux à 12. 6 [029] Ainsi, la très grande précision des valeurs de la tension complexe est obtenue avec un nombre très réduit d'itérations. Le procédé est donc rapide, efficace et peu coûteux en ressources de calcul. 5 [30] La présente invention vise aussi un système de détermination d'un état en tension d'un réseau de distribution basse tension au niveau d'au moins un poste de transformation, la branche basse tension comprenant une pluralité de noeuds, le système comprenant un circuit de traitement configuré pour la mise en oeuvre du procédé selon l'invention. 10 [31] D'autres avantages et caractéristiques de l'invention apparaitront à la lecture de la description détaillée ci-après d'exemples de réalisation de l'invention, et à l'examen des dessins annexés sur lesquels : 15 [032] la figure 1 illustre schématiquement le contexte dans lequel se place la présente invention, [33] la figure 2 représente schématiquement un système de détermination d'un état en tension d'un réseau de distribution basse tension au niveau d'au moins un poste de 20 transformation, la branche basse tension comprenant une pluralité de noeuds, selon un exemple de réalisation. [34] la figure 3 est un ordinogramme illustrant les principales étapes du procédé de détermination d'un état en tension d'un réseau de distribution basse tension au niveau d'au 25 moins un poste de transformation, la branche basse tension comprenant une pluralité de noeuds, selon un exemple de réalisation. [35] La figure 1 illustre un réseau de distribution d'électricité comprenant au moins une branche basse tension BT et au moins une branche haute tension HTA, HTB. Le réseau de 30 distribution est un réseau arborescent comprenant une pluralité de noeuds au niveau desquels l'état en tension est à déterminer. [36] La branche haute tension HTA est une branche dans laquelle la tension est comprise entre 1 kilovolt et 50 kilovolts en courant alternatif, et comprise entre 1.5 kilovolt et 75 35 kilovolts en courant continu. La branche haute tension HTB est une branche dans laquelle la 7 tension excède 50 kilovolts en courant alternatif et 75 kilovolts en courant continu. La tension dans la branche basse tension est comprise entre 50 volts et 1 kilovolt. 5 [037] Les branches haute tension et basse tension sont reliées à une pluralité d'installations INST, typiquement des habitations, des entreprises, des centrales électriques, etc. Le réseau électrique est ponctué de transformateurs et de régleurs en charge, permettant de récolter les données de consommation et de production d'électricité et de réguler la tension sur le réseau. Dans la présente application, seule la tension sur la branche basse tension est à déterminer et 10 à réguler. [038] La figure 2 illustre plus précisément un système de détermination d'un état en tension d'un réseau de distribution basse tension. 15 [039] La branche du réseau représentée sur la figure 2 est la branche basse tension. Une pluralité de générateurs d'énergie dispersée GED sont reliés au réseau. Les générateurs d'énergie dispersée GED sont par exemple des générateurs d'énergie renouvelable, tels que des installations photovoltaïques ou éoliennes. Les générateurs GED peuvent être des installations appartenant à des particuliers, branchées sur le réseau. Des compteurs électriques 20 intelligents COMPT peuvent être reliés à chacun des générateurs, de sorte que les puissances actives et réactives P, Q peuvent être obtenues en chaque noeud k du réseau. [040] En variante, les puissances actives et réactives P, Q sont obtenues au moyen d'un modèle de charge propre à chaque utilisateur du réseau. Les puissances actives et réactives P, 25 Q peuvent être transférées par courants porteurs en ligne jusqu'au transformateur IRANSF. Ces valeurs sont ensuite transmises vers le concentrateur CONC comprenant une unité de traitement TRAIT. Plus précisément, l'unité de traitement TRAIT reçoit en entrée les données récoltées en chaque noeud du réseau. L'unité de traitement TRAIT est configurée de manière à mettre en oeuvre les étapes du procédé décrites ci-après en référence à la figure 3. L'unité de 30 traitement TRAIT délivre donc en sortie la tension complexe en chaque noeud de la branche basse tension, lorsqu'elle appartient à la plage de tensions contractuelle. L'unité de traitement TRAIT peut également délivrer en sortie la valeur de consigne à appliquer pour réguler l'état en tension le cas échéant. La valeur de consigne est envoyée vers le régleur en charge basse tension REG, qui reçoit également les puissances actives et réactives P, Q. Le régleur en 35 charge basse tension REG permet de réguler l'état en tension du réseau afin que la tension en chaque noeud soit comprise dans la plage de tensions contractuelle. 8 [41] Dans le cadre de l'exemple de réalisation illustré sur la figure 3, on considère un réseau de distribution basse tension à structure arborescente donc l'ensemble des noeuds 5 N = {1 2,s n) est décomposé en N = 11)U NpQ, avec Cl) le bus de référence et Neri est le sous-ensemble de noeuds de consommation et/ou de production la charge en un noeud k Emn étant caractérisée par la puissance active P et la puissance réactive eh [42] La tension en chaque noeud du réseau de distribution est calculée en fonction des 10 puissances actives et réactives injectées ou soutirées dans chaque noeud. Cela est effectué au moyen de la résolution des équations dites « de Kirchhoff », ou résolution de problème du calcul de la répartition des charges, communément appelé problème du « load flow ». Les équations de Kirchhoff permettent de lier une tension complexe aux puissances actives et réactives du noeud associé. La construction du système d'équations de Kirchhoff est réalisée 15 de manière connue, pour obtenir le premier système d'équations qui suit, à l'étape Si : Sk= Pk Qk = [ diag(V)(YV)* k (k=1,...n) V1= VI 01 = 0 Sk = Sk (VkENpQ), 20 [43] Avec Sk la puissance complexe injectée en chaque noeud k, Pk et Qk les puissances actives et réactives, données, en chaque noeud k, V est le vecteur des tensions complexes en chaque noeud k, V1 la tension complexe au noeud de référence et et la phase en ce même noeud de référence. Y est la matrice d'admittance complexe du réseau. 25 [44] A l'étape S2, les puissances actives et réactives P, Q, injectées en chaque noeud du réseau, sont estimées. Plus précisément, le bilan de consommation et de production de chaque client est déterminé, de manière à obtenir les puissances actives et réactives P, Q. Notamment, et comme illustré sur la figure 2, les puissances actives et réactives P, Q peuvent être 30 directement fournies au moyen de compteurs électriques intelligents COMPT. Ces données transitent par exemple par courant porteur en ligne et sont transférées à un concentrateur CONC du poste de transformateur en tension TRANSF. Les données sont récupérées au niveau du concentrateur CONC. En variante, une mesure de tension peut être faite au niveau de chaque poste de distribution par le concentrateur CONC. Lorsque des clients ne possèdent 9 pas de compteurs électriques intelligents COMPT, les données peuvent être estimées au moyen d'un modèle d'estimation de charge ou de profil de production. Les puissances actives et réactives peuvent être estimées à chaque intervalle de temps, par exemple au moins égal à 5 30 minutes et de préférence égal à 10 minutes. [45] A partir des valeurs des puissances actives et réactives P, Q il est possible d'évaluer l'état en tension de chaque noeud du réseau de manière précise. Toutefois, les valeurs en tension dans le réseau électrique basse tension ne dépendent pas linéairement des puissances 10 soutirées ou injectées par les générateurs d'énergies renouvelables GED. Aussi, le premier système d'équations, comprenant les équations de Kirchhoff, est reformulé de manière à pouvoir déterminer l'état en tension de chaque noeud. [46] Plus précisément, les puissances actives et réactives sont reformulées sous la forme 15 d'une fonction quadratique pour chaque noeud. Ces fonctions quadratiques se substituant aux puissances actives et réactives dans le premier système d'équations à l'étape S3 pour obtenir un deuxième système d'équations en nombres réels. Cette substitution est fondée sur le fait que pour chaque k E NPQ, la partie réelle Pk de la puissance injectée/soutirée au noeud k peut être exprimée comme une fonction quadratique: R211 -> R de la forme xr Ak x. De la même 20 manière, la partie imaginaire Qk peut être exprimée comme une fonction quadratique: R2" R de la forme xT Bk x, où x est le vecteur de dimension 2n dont les composantes 1 à n sont les parties réelles des tensions aux noeuds 1,...,n et les composantes n+1 à 2n sont les parties imaginaires. AL et Bk sont respectivement les matrices symétriques réelle et imaginaire 2n x 2n, déduites de la matrice d'admittance Y complexe du réseau de la manière suivante. Pour 25 tout i= I, ...n, Yrn est la matrice nxn ayant toutes ses lignes nulles sauf la ligne k égale à la kème ligne de Y. R et C sont les parties réelle et complexe de Y[k] et (DL et LPk sont les matrices 2nx2n: R -C C R 30 -C - R R - C (1)k = irk = 35 10 [047] Les matrices Ak et Bk peuvent alors être réécrites comme suit Ak= 0.5x (d)k+ il)kr) et Bk = 0.5x (lik + 'PO). 5 [48] Ainsi, résoudre le système (I) se reformule comme la recherche d'une solution au deuxième système d'équations: xT Ak X = Pk (V k=2, xl Bk X = Qk (V k=2,...n) 10 X1 - V1 [49] La résolution du deuxième système d'équations décrit ci-avant comprend deux phases. La première phase est décrite en référence aux étapes S4 à S6. La deuxième phase est décrite 15 en référence à l'étape S7. Chacune des deux phases comprend un calcul itératif permettant de résoudre avec précision les équations du deuxième système, tout en contrôlant chaque étape de la résolution. [50] Avantageusement, le degré très important de précision est obtenu suite à la 20 minimisation itérative d'une fonction auxiliaire, déterminée à l'étape S4. Plus précisément, plutôt que de résoudre le deuxième système d'équation directement, un problème d'optimisation est formulé, dont une solution optimale approchée peut servir de solution initiale à la deuxième phase. 25 [051] La fonction auxiliaire est formulée comme suit : (x) = Max ( Maxi=t n T ; Maxi=i,...n X1 Bk X- Qk I X ±1k: X - [052] L'équivalence entre la résolution du deuxième système d'équations et la minimisation 30 de la fonction auxiliaire d) (x) découle du fait que la fonction auxiliaire d) (x) est partout réelle, positive ou nulle, et qu'elle s'annule en toute solution du deuxième système d'équations. Autrement dit, toute solution du deuxième système d'équations correspond à un minimum global de la fonction auxiliaire, lorsque celle-ci comprend plusieurs minima globaux de valeur nulle, chacun de ces minima est solution du deuxième système d'équations. 35 11 [53] La minimisation de la fonction auxiliaire n'a toutefois pas pour but de déterminer un minimum global précis. Une approximation d'un minimum global est suffisante pour obtenir un résultat permettant de servir de solution initiale à la phase 2. Aussi, une solution de la 5 minimisation de la fonction auxiliaire dont la valeur est suffisamment proche de 0 est suffisante. Ainsi on se contente en phase 1 d'une solution au voisinage d'un minimum global de la fonction auxiliaire. La deuxième phase du procédé permet quant à elle d'obtenir l'état en tension de manière très précise. Aussi un nombre réduit d'itérations permet d'obtenir une solution satisfaisante. Dans la suite de la description, le terme « minimum global » est utilisé 10 bien qu'il est entendu qu'il peut s'agir d'une approximation de ce minimum global. [54] Par exemple, le nombre d'itérations pour la phase 1 est inférieur à 20,de préférence inférieur à 15, et de préférence égal à 12. Cela permet donc de limiter les ressources en calcul nécessaires car une approximation du minimum global est atteinte avec un nombre réduit 15 d'itérations. Ce nombre réduit d'itérations permet de diminuer sensiblement la valeur de la fonction auxiliaire (1) (x) par rapport à sa valeur initiale (D (x0). [55] Afin d'obtenir une précision et une valeur de l'état en tension en chaque noeud du réseau qui soient satisfaisantes, une tension-cible réelle Tcible, ou « tension-cible Tcible » est 20 utilisée comme solution initiale de la fonction auxiliaire. Plus précisément, la tension-cible Tcible est utilisée comme solution à la première itération de la minimisation de la fonction auxiliaire. [56] La tension Tcible est définie à l'étape S5. La tension Tcible est une valeur idéale de la 25 tension en chacun des noeuds. C'est la valeur en tension qu'il s'agit d'atteindre par le procédé de régulation selon l'invention. Plus précisément, la tension cible Tcible prend la forme d'un vecteur dont le nombre de composants est égal au nombre de noeuds. La tension cible Tcible peut être la même pour chaque noeud du réseau. En variante, la tension cible Tcible n'est pas la même en chaque noeud. Cela peut être dû à des évolutions des valeurs de tensions 30 contractuelles, pouvant différer selon les clients du réseau. Ainsi, la valeur initiale de la fonction auxiliaire (I) (x) est (D (x0), avec x0 réel correspondant à Tcible complexe. [57] Pour la résolution du problème d'optimisation formulé à travers la fonction auxiliaire t1 (x), plusieurs méthodes peuvent être utilisées. Par exemple, une méthode de gradient projeté 35 peut être utilisée. 12 [58] Plus précisément, à chaque itération t, deux opérations sont effectuées. La première itération tO comprend la détermination d'une direction de descente d de la fonction auxiliaire 5 (1) (x) au point courant xm. La détermination de la direction de descente d peut être formulée et résolue comme un problème d'optimisation à objectif quadratique convexe sous contraintes linéaires. [59] La deuxième opération comprend la détermination d'un nouvel itéré xt comme to minimum exact ou approché de la fonction auxiliaire (1) (x) sur la demi-droite issue de xm et de la direction de descente d. [60] Toutefois, la fonction auxiliaire CD( X").(1 - E) est l'ensemble des k tels que (xm)rAk - Pk x1-1). 0 J est l'ensemble des k tels que (xmak - Qk REVENDICATIONS1. Procédé, mis en oeuvre par des moyens informatiques, de détermination d'un état en tension d'un réseau de distribution basse tension au niveau d'au moins un poste de transformation, la branche basse tension comprenant une pluralité de noeuds, le procédé comprenant: - une définition d'une valeur d'une tension cible réelle pour chaque noeud ; - une estimation d'une puissance active et d'une puissance réactive pour chaque noeud ; - une construction d'un premier système d'équations liant des tensions complexes aux puissances actives et réactives de la pluralité de noeuds, - une formulation respective desdites puissances actives et réactives sous la forme d'une fonction quadratique pour chaque noeud, lesdites fonctions quadratiques étant substituées dans le premier système d'équations de manière à obtenir un deuxième système d'équations en nombre réels, équivalent au premier système d'équations, - une détermination d'une fonction auxiliaire issue du deuxième système d'équations, le procédé comprenant: - une minimisation itérative de ladite fonction auxiliaire, dont une solution initiale à la première itération est la tension cible, jusqu'à obtenir une solution au voisinage d'un minimum global approché de ladite fonction auxiliaire, - et une application au deuxième système d'équations d'une méthode itérative par linéarisations successives, ayant comme solution initiale à la première itération ledit minimum global approché, jusqu'à obtenir une valeur d'une tension complexe en chaque noeud. 2. Procédé, selon la revendication 1, comprenant en outre : - un calcul d'une variation entre la valeur de ladite tension cible réelle et la valeur d'un module de ladite tension complexe en chaque noeud; et - si ladite variation est inférieure à une valeur absolue d'un premier seuil, aucune action de régulation n'est entreprise, - si ladite variation en valeur absolue est comprise entre le premier seuil et le deuxième seuil, la valeur absolue du deuxième seuil étant supérieure à la valeur absolue du premier seuil, une détermination de consignes de tension à appliquer au poste de transformation du réseau de manière à ramener la tension complexe en chaque noeud de la branche basse tension entre des valeurs de contrainte de tension 17 maximale et minimale ; et une transmission des consignes de tension au poste de transformation, - si ladite variation en valeur absolue est supérieure au deuxième seuil, le procédé 5 est relancé. 3. Procédé selon la revendication 2, dans lequel : - la valeur absolue du premier seuil est égale à 10% de la valeur de la tension cible réelle, et 10 Y la valeur absolue du deuxième seuil est égale à 15% de la valeur de la tension cible réelle. 4. Procédé, selon l'une des revendications 1 à 3, dans lequel les puissances actives et réactives en chaque noeud de la branche basse tension sont déterminées : 15 - en fonction de mesures réelles fournies par au moins un compteur électrique relié au noeud associé de la branche basse tension, ou - en fonction d'un modèle de charge de consommateurs. 5. Procédé selon l'une des revendications 1 à 4, caractérisé en ce que les étapes du 20 procédé sont répétées à chaque intervalle de temps, de préférence inférieur à trente minutes, et de préférence inférieur à dix minutes. 6. Procédé selon l'une des revendications 1 à 5, dans lequel la fonction auxiliaire est choisie de sorte qu'un minimum global de valeur 0 de ladite fonction auxiliaire est une 25 solution du deuxième système d'équations, ladite fonction auxiliaire étant à valeur réelle positive ou nulle. 7. Procédé selon la revendication 6, dans lequel une résolution du deuxième système d'équations comprend la minimisation de la fonction auxiliaire par une méthode itérative 30 comprenant, à chaque itération, une recherche d'une direction de descente puis une minimisation unidimensionnelle de la fonction auxiliaire jusqu'à obtenir la solution au voisinage du minimum global. 8. Procédé selon la revendication 7, dans lequel la direction de descente est trouvée, à 35 chaque itération, par une résolution d'un problème d'optimisation à objectif quadratique convexe sous contraintes linéaires. 18 9. Procédé selon l'une des revendications 1 à 8, dans lequel la méthode itérative de linéarisations successives comprend, à chaque itération : 5 ;"-- une résolution d'une approximation linéaire du deuxième système d'équations au voisinage d'une solution obtenue à l'itération précédente, laquelle définit une direction de déplacement dans l'espace des solutions, et - une minimisation unidimensionnelle de la fonction auxiliaire dans cette direction de déplacement. 10 10. Procédé selon l'une des revendications 1 à 9, dans lequel la valeur de la tension complexe est obtenue avec un degré de précision compris entre 10-6 et 10-15 11 Procédé selon l'une des revendications 1 à 10, caractérisé en ce que les nombres 15 d'itérations pour la résolution de la fonction auxiliaire et d'itérations de la méthode de linéarisations successives sont chacun inférieurs à 20, de préférence inférieurs à 15, et de préférence égaux à 12. 12. Système de détermination d'un état en tension d'un réseau de distribution basse 20 tension au niveau d'au moins un poste de transformation, la branche basse tension comprenant une pluralité de noeuds, le système comprenant : - un circuit de traitement configuré pour la mise en oeuvre du procédé selon l'une des revendications 1 à11.