FR 3046280 A1 20170630 FR 1563215 A 20151223 GESTION D'UN APPROVISIONNEMENT POUR DES CONSOMMATIONS SE PRODUISANT EN UNE MULTIPLICITE DE SITES La présente invention concerne un procédé et un système de gestion d'un approvisionnement pour des consommations qui se produisent en une multiplicité de sites. En particulier, elle peut s'appliquer pour des consommations en énergie, combustible, ou eau, ou pour des consommations de services qui sont délivrés localement, tels que des services d'enlèvement d'ordures, d'élimination d'effluents, de téléphonie, de transmission par un réseau de communication, etc. De telles consommations qui se produisent en des sites nombreux et dispersés, et qui se produisent au cours du temps, posent des difficultés de comptage de la consommation qui est associée individuellement à chaque site pendant une période déterminée. Pour cela, chaque site est équipé d'un compteur qui mesure la consommation de ce site, saisit des valeurs mesurées à des instants de relevés différents, par exemple périodiquement, et qui transmet ces valeurs à une unité centralisée de traitement de données. Usuellement, ces valeurs de consommation qui sont collectées peuvent être de deux types. Le premier type correspond à des valeurs de cumuls de consommation qui sont saisies aux instants des relevés, et qui sont appelées index dans le jargon de l'Homme du métier. Le second type correspond à des valeurs de progressions de consommation entre deux instants de relevés différents. Ainsi, pour un site donné, une valeur de progression de la consommation entre deux instants de relevés est égale à la variation d'index entre ces instants. Or ces valeurs collectées constituent des quantités de données très importantes. Mais simultanément, un grand nombre de valeurs de progressions de consommation sont manquantes, notamment parce que les consommations de tous les sites ne sont pas relevées simultanément, ni avec des périodicités identiques. Ces lacunes dans les données collectées proviennent notamment de la diversité des types de compteurs qui ont été installés progressivement dans les sites. 3046280 - 2 - Néanmoins, la collecte de telles valeurs de consommation est essentielle pour gérer l'approvisionnement, notamment pour assurer que cet approvisionnement soit suffisant par rapport aux besoins de tous les sites, et ce à chaque instant. En particulier, la gestion d'un tel approvisionnement en 5 l'objet de la consommation nécessite de dimensionner la production de cet objet, ou la capacité de fourniture du service dans d'autres cas, de prévoir l'acheminement de l'objet de la consommation, avec éventuellement des délestages et/ou des stockages temporaires de cet objet. En plus de cette importance pour l'opérateur de l'approvisionnement, les valeurs de 10 consommation peuvent être utilisées pour : - produire des bilans de consommation à destination de collectivités locales, qui correspondent à des consommations intégrées sur le territoire de compétence de chaque collectivité locale, et pendant des durées d'intégration qui peuvent être journalières, mensuelles ou 15 annuelles ; - optimiser l'approvisionnement des villes, qui sont sujettes à des pics de consommation, et donc pour lesquelles il est plus difficile de parvenir à un équilibre local entre la mise à disposition de l'objet de la consommation et le besoin ; 20 - anticiper des anomalies de consommation, telles que des pics de consommation ou à l'inverse des durées pendant lesquelles la consommation serait trop faible par rapport à la mise à disposition de l'objet de la consommation ; - détecter des anomalies dans la consommation, en comparant les 25 valeurs collectées à des valeurs prévisionnelles ; - facturer des consommateurs, ou la possibilité de réaliser une facturation prévisionnelle d'après une consommation qui est estimée, à défaut d'avoir été mesurée ; et - réaliser des études statistiques, notamment pour établir des corrélations entre les variations d'une consommation et des phénomènes externes, telles que la température ambiante extérieure, le recouvrement avec une période de congés, etc. 3046280 - 3 - Les difficultés pour l'exploitation des valeurs de consommation qui sont collectées proviennent alors des caractéristiques suivantes : - les données collectées regroupent aléatoirement des valeurs de cumuls de consommation et des valeurs de progressions de consommation ; 5 - les durées qui séparent deux relevés successifs de consommation qui sont relatifs à un même site peuvent être irrégulières, et différentes des durées relatives à un autre site ; - certains sites peuvent ne pas être équipés de compteurs individuels, mais peuvent être regroupés avec un compteur qui est commun à ces 10 sites ; - bien qu'étant équipés de compteurs individuels, certains sites peuvent néanmoins faire l'objet de relevés de consommation qui sont regroupés avec d'autres sites, par exemple à partir d'un noeud de distribution ; - la maille temporelle qui est requise pour exploiter des données 15 consommation peut ne pas correspondre aux durées entre relevés successifs ; et - la maille spatiale qui est requise pour exploiter des données consommation peut ne pas correspondre à l'attribution des compteurs aux sites géographiques réels. 20 Ces difficultés seraient résolues en mesurant les consommations individuelles de tous les sites à des instants qui soient suffisamment resserrés pour que les valeurs utiles puissent être obtenues en additionnant certaines des valeurs mesurées pour retrouver les mailles temporelle et spatiale voulues. Mais une telle solution n'est pas compatible avec la réalité de la collecte des 25 valeurs de consommation. Des méthodes d'échantillonnage, d'interpolation avec corrélation, et des méthodes de détermination de profils sont déjà développées pour combler les lacunes dans les valeurs de consommation qui sont mesurées et collectées. Il est ainsi possible d'obtenir des évaluations de consommation 30 selon des mailles temporelle et spatiale qui sont plus fines que celles des valeurs mesurées. Néanmoins, la précision des valeurs non-mesurées qui sont 3046280 - 4 - ainsi évaluées n'est pas toujours suffisante. En particulier, les méthodes d'échantillonnage nécessitent que la proportion des valeurs manquantes reste assez faible, et les méthodes de profilage nécessitent qu'un opérateur explicite et valide chaque profil de consommation qui est utilisé. 5 A partir de cette situation, la présente invention a pour but de proposer une nouvelle méthode pour disposer d'évaluations de consommation selon des mailles temporelle et spatiale plus fines, éventuellement jusqu'à chaque site individuellement et jusqu'à la plus petite durée qui sépare deux relevés de consommation chacun effectué sur l'un quelconque des sites. 10 Un but annexe est de disposer d'une telle méthode qui puisse être mise en oeuvre automatiquement, sans intervention d'un opérateur pour trier ou sélectionner les valeurs collectées. Un autre but est de disposer d'une telle méthode qui ne nécessite pas d'équipements importants, notamment peu d'équipement à installer sur chaque 15 site. Pour atteindre ces buts ou d'autres, un premier aspect de la présente invention propose un procédé de gestion d'un approvisionnement en énergie, combustible, ou eau, ou en services délivrés localement, tels que des services d'enlèvement d'ordures, d'élimination d'effluents, de téléphonie, ou de 20 transmission par un réseau de communication, pour des consommations qui se produisent en une multiplicité de sites géographiques et pendant une période de surveillance. Ce procédé de l'invention comprend les étapes suivantes : /1/ collecter des valeurs de cumuls de consommation, associées chacune à un des sites et à un instant de relevé de consommation qui est 25 compris dans la période de surveillance, et éventuellement collecter aussi des valeurs de progressions de consommation, associées chacune à un ou plusieurs des sites et correspondant à une consommation qui a été effectuée entre deux instants différents de relevés de consommation compris dans la période de surveillance ; 30 /2/ transmettre les valeurs collectées à une unité de traitement de données ; 3046280 - 5 - /3/ concevoir une matrice des progressions de consommation ayant des colonnes affectées une-à-une aux sites, et des lignes affectées une-àune à des couples d'instants successifs, pour chacun desquels instants un relevé de consommation a été effectué pour l'un quelconque des 5 sites, cette matrice étant destinée à contenir des valeurs de progressions de consommation attribuées chacune à un des sites pour une consommation qui a été effectuée à ce site entre deux instants de relevés de consommation successifs, et cette matrice des progressions de consommation étant initialement lacunaire à cause de valeurs 10 collectées manquantes pour les progressions de consommation de certains au moins des sites entre deux instants successifs de relevés de consommation ; /4/ exprimer les valeurs collectées sous forme de contraintes à satisfaire par les valeurs de progressions de consommation ; 15 /5/ par l'unité de traitement de données : déterminer deux matrices non- négatives qui fournissent, à partir d'une multiplication de ces deux matrices non-négatives, une matrice d'estimation qui possède des mêmes nombres de colonnes et de lignes que la matrice des progressions de consommation, et dont des valeurs qui sont situées à 20 des intersections de ses colonnes et de ses lignes satisfont les contraintes, avec la première des deux matrices non-négatives dans la multiplication qui possède un nombre de colonnes égal à un nombre de lignes de la seconde matrice non-négative dans la multiplication, et qui est inférieur aux deux nombres de colonnes et de lignes de la matrice 25 des progressions de consommation ; et /6/ utiliser au moins une valeur de la matrice d'estimation en tant qu'évaluation d'une des valeurs de progressions de consommation non-collectées, pour gérer l'approvisionnement des sites. Ainsi, le procédé de l'invention est compatible avec les deux types de 30 valeurs de consommation qui peuvent être collectées : les cumuls de consommation, ou index, et les progressions de consommation. En particulier, l'équipement qui est nécessaire à chaque site pour mettre en oeuvre l'invention, 3046280 - 6 - est limité à un compteur, et les modèles de compteurs déjà déployés peuvent être utilisés sans modification. Selon l'invention, la matrice d'estimation constitue une évaluation de la matrice des progressions de consommation. Autrement dit, chaque valeur de la 5 matrice d'estimation est proche de celle de la matrice des progressions de consommation qui est située à la même intersection de colonne et de ligne. Néanmoins, la matrice d'estimation est obtenue à partir de la multiplication des deux matrices non-négatives, qui ont chacune une de leurs dimensions, c'est-à-dire un parmi leurs nombres de colonnes ou de lignes, plus petite que les 10 deux dimensions de la matrice des progressions de consommation. En outre, par son principe, le procédé de l'invention produit des évaluations des progressions de consommation selon les tailles de mailles spatiale et temporelle qui sont les plus petites parmi celles des valeurs collectées. 15 De plus, le procédé de l'invention peut être exécuté de façon entièrement automatisée. Dans des modes de mise en oeuvre préférés de l'invention, les deux matrices non-négatives et la matrice d'estimation peuvent être déterminées à l'étape /5/, lors de répétitions d'une séquence qui comprend les sous-étapes 20 suivantes, effectuées en commençant à partir de deux matrices non-négatives initiales, et d'une matrice d'estimation initiale qui est déduite elle-même du produit des deux matrices non-négatives initiales pour satisfaire les contraintes : /5-1/ en utilisant les deux matrices non-négatives et la matrice 25 d'estimation telles que résultant de l'exécution précédente de la séquence, modifier l'une des deux matrices non-négatives en maintenant constante l'autre matrice non-négative, de façon à réduire un écart entre le produit des deux matrices non-négatives et la matrice d'estimation ; 30 /5-2/ en utilisant la matrice non-négative modifiée à la sous-étape /5-1/, et la matrice d'estimation telle que résultant de l'exécution précédente de la séquence, modifier l'autre des deux matrices non-négatives en 3046280 - 7 - maintenant constante la matrice non-négative modifiée à la sous-étape /5-1/, de façon à réduire encore l'écart entre le produit des deux matrices non-négatives et la matrice d'estimation ; puis /5-3/ calculer le produit des deux matrices non-négatives obtenues à 5 l'issue des sous-étapes /5-1/ et /5-2/, respectivement, et modifier la matrice de produit pour satisfaire les contraintes, la matrice de produit modifiée formant alors la matrice d'estimation qui résulte de l'exécution de la séquence. Ainsi, la précision avec laquelle la matrice des progressions de 10 consommation est déterminée, est améliorée progressivement au fur et à mesure que des répétitions de la séquence des sous-étapes /5-1/ à /5-3/ sont ajoutées. Des premières évaluations de la matrice des progressions de consommation sont donc disponibles assez rapidement, après un nombre réduit de répétitions de la séquence, bien que pas nécessairement précises ni 15 fiables. Puis la précision et la fiabilité sont améliorées par les matrices d'estimation qui sont obtenues successivement lors des répétitions ultérieures de la séquence des sous-étapes /5-1/ à /5-3/. Chacune des sous-étapes /5-1/ et /5-2/ peut être exécutée en utilisant un algorithme d'optimisation. Cet algorithme d'optimisation peut être de type à 20 suivi de gradient optimal. Ce peut aussi être un algorithme récursif. A la sous-étape /5-3/, la matrice de produit peut être modifiée pour satisfaire les contraintes en appliquant une projection orthogonale à cette matrice de produit. Cette projection aboutit sur un sous-espace affine qui regroupe toutes les matrices ayant les mêmes nombres de colonnes et de 25 lignes que la matrice des progressions de consommation, et qui satisfont les contraintes. La matrice de produit projetée forme alors la matrice d'estimation qui résulte de l'exécution de la séquence des sous-étapes /5-1/ à /5-3/. Avantageusement, l'étape /2/ peut comprendre un enregistrement intermédiaire d'une au moins des valeurs qui sont collectées à l'un des sites 30 géographiques, sur un support d'enregistrement ou de mémoire qui est disponible à ce site géographique. Ensuite, chaque valeur collectée et enregistrée peut être lue puis transmise à destination de l'unité de traitement 3046280 - 8 - de données. Une telle transmission qui est différée par rapport à la mesure de la valeur de consommation, c'est-à-dire différée par rapport à l'instant du relevé de la consommation, permet de pallier des interruptions du réseau de transmission qui relie le site de consommation à l'unité de traitement de 5 données. En outre, un tel enregistrement intermédiaire permet d'établir localement, c'est-à-dire au niveau du site, un historique de la consommation de ce site, puis de transmettre ultérieurement cet historique de consommation à l'unité de traitement de données. Une telle transmission d'historique de consommation peut être utile, par exemple, lors d'un nouveau rattachement du 10 site au réseau de distribution qui délivre l'objet de la consommation. Il est possible que certaines utilisations des valeurs de progressions de consommation qui sont évaluées selon l'invention, ne nécessitent pas des mailles spatiales ou temporelles qui soient aussi fines que celles de la matrice des progressions de consommation. Dans ce cas, le procédé peut comprendre 15 en outre les étapes suivantes : /7/ sélectionner un sous-ensemble de sites parmi la multiplicité des sites géographiques, et sélectionner une période d'intégration à l'intérieur de la période de surveillance ; puis /8/ produire une valeur d'estimation de la consommation totale associée 20 aux sites du sous-ensemble pendant la durée d'intégration, en additionnant les valeurs de la matrice d'estimation qui appartiennent à la fois aux colonnes affectées aux sites du sous-ensemble et aux lignes affectées aux couples d'instants successifs qui sont contenus dans la période d'intégration. 25 Dans le jargon de l'Homme du métier, l'étape /8/ consiste à ré-agréger les évaluations des valeurs de progressions de consommation, selon la maille spatio-temporelle qui a été choisie à l'étape /7/. La précision de l'évaluation pour la consommation totale qui est associée aux sites du sous-ensemble pendant la durée d'intégration, est en général supérieure à la précision 30 individuelle de chaque valeur de la matrice d'estimation. Il est possible d'évaluer la précision d'une des valeurs de la matrice d'estimation, en effectuant en outre les étapes /9/ à /11/ suivantes, de façon 3046280 - 9 - optionnelle. Ces étapes supplémentaires sont toutes exécutées en utilisant la même multiplicité de sites et les mêmes instants de relevés de consommation pour la matrice des progressions de consommation : /9/ sélectionner des groupes différents de valeurs collectées, parmi les 5 valeurs de cumuls de consommation associées chacune à un des sites et à un instant de relevé de consommation qui est compris dans la période de surveillance, et éventuellement aussi parmi les valeurs de progressions de consommation associées chacune à un ou plusieurs des sites et correspondant à une consommation qui a été effectuée 10 entre deux instants différents de relevés de consommation compris dans la période de surveillance ; /10/ pour chaque groupe de valeurs collectées, répéter les étapes /1/ à /6/ pour obtenir une matrice d'estimation indépendamment pour chaque groupe de valeurs collectées, les matrices d'estimation obtenues 15 correspondant toutes à la même matrice des progressions de consommation ; puis /11/ calculer au moins une valeur d'erreur moyenne ou d'écart-type à partir des matrices d'estimation obtenues pour tous les groupes de valeurs collectées. 20 Cette valeur d'erreur moyenne ou d'écart-type est une mesure de la précision et de la fiabilité de l'évaluation de progression de consommation que fournit la matrice d'estimation. Dans le jargon de l'Homme du métier, une telle mesure de précision et de fiabilité est appelée «validation croisée». De préférence, cette mesure de précision et de fiabilité peut être 25 effectuée sur une valeur de ré-agrégation. Pour cela, l'étape /10/ peut comprendre en outre l'étape /8/ qui est répétée pour chaque groupe de valeurs collectées, pour le même sous-ensemble de sites et la même période d'intégration, tels que sélectionnés à l'étape /7/. La valeur d'erreur moyenne ou d'écart-type peut alors être calculée à l'étape /11/ pour les valeurs d'estimation 30 de la consommation totale qui est associée aux sites du sous-ensemble pendant la période d'intégration. Eventuellement, le procédé de l'invention peut être amélioré en prenant 3046280 -10- en compte des valeurs exogènes aux consommations. Le procédé comprend alors la collecte des valeurs exogènes en plus des valeurs de cumuls de consommation et de progressions de consommation, et l'expression des valeurs exogènes collectées sous forme de contraintes supplémentaires à 5 satisfaire par les valeurs de progressions de consommation. Les deux matrices non-négatives sont alors déterminées à l'étape /5/ de sorte que la matrice d'estimation satisfasse simultanément les contraintes issues des valeurs collectées de cumuls de consommation et de progressions de consommation, mais aussi les contraintes issues des valeurs exogènes collectées. 10 L'invention peut être appliquée avantageusement à la consommation d'énergie électrique à l'intérieur d'une zone géographique. En outre, quelque soit la consommation à laquelle l'invention est appliquée, l'une au moins des valeurs de la matrice d'estimation peut être utilisée à l'étape /6/ en tant qu'évaluation d'une des valeurs non-collectées de 15 progressions de consommation, pour dimensionner une production, un acheminement, un délestage ou un stockage d'un objet de la consommation, afin d'assurer une fourniture de celui-ci à l'un au moins des sites en quantité suffisante par rapport à un besoin en ce site. Un second aspect de l'invention propose un système de gestion d'un 20 approvisionnement en énergie, combustible, ou eau, ou en services délivrés localement, tels que des services d'enlèvement d'ordures, d'élimination d'effluents, de téléphonie, de transmission par un réseau de communication, pour des consommations qui se produisent en une multiplicité de sites géographiques et pendant une période de surveillance. Ce système de 25 l'invention comprend : - des moyens pour collecter des valeurs de cumuls de consommation, associées chacune à un des sites et à un instant de relevé de consommation qui est compris dans la période de surveillance, et éventuellement aussi pour collecter des valeurs de progressions de 30 consommation, associées chacune à un ou plusieurs des sites et correspondant à une consommation qui a été effectuée entre deux 3046280 -11- instants différents de relevés de consommation compris dans la période de surveillance ; - des moyens de transmission, qui sont adaptés pour transmettre les valeurs collectées à une unité de traitement de données ; 5 - l'unité de traitement de données, qui est adaptée pour mettre en oeuvre au moins les étapes /5/ et /6/ d'un procédé conforme au premier aspect de l'invention ; et - optionnellement, des moyens d'enregistrement ou de mémorisation qui sont situés à l'un au moins des sites, et qui sont agencés pour réaliser 10 un enregistrement intermédiaire d'une au moins des valeurs collectées à ce site, puis une lecture de la valeur collectée et enregistrée en vue de sa transmission à destination de l'unité de traitement de données. Possiblement, plusieurs valeurs collectées puis enregistrées localement au niveau du site peuvent être transmises ultérieurement à 15 l'unité de traitement de données sous forme d'un historique de consommation du site. D'autres particularités et avantages de la présente invention apparaîtront dans la description ci-après d'exemples de mise en oeuvre non limitatifs, en référence aux dessins annexés, dans lesquels : 20 - la figure 1 représente un réseau de collecte de données de consommation ; - les figures 2a à 2c représentent des matrices et des opérations matricielles qui sont utilisées dans un procédé conforme à l'invention ; et 25 - la figure 3 est un diagramme d'étapes pour une mise en oeuvre particulière d'un procédé conforme à l'invention. Bien que l'invention soit décrite en détail dans la suite dans le cas particulier de consommations d'énergie électrique, il est clair qu'elle peut être transposée à des consommations d'autres types, par exemple des 30 consommations de gaz combustible, d'eau, ou des consommations de services tels que l'enlèvement d'ordures, l'élimination d'effluents, la téléphonie par 3046280 -12- lignes fixes ou la téléphonie mobile, la location de films en ligne, désignée par VOD, etc. L'invention est adaptée dès que la consommation se produit en des sites qui sont séparés géographiquement, et qu'elle peut être mesurée en fonction du temps, par exemple sous la forme d'une quantité d'énergie 5 électrique livrée, ou d'un volume de gaz combustible ou d'eau consommée, ou sous la forme d'un poids d'ordures ou d'un volume d'effluents éliminés, ou encore sous la forme de durées cumulées de communications téléphoniques ou de quantités de données qui sont transmises au consommateur à chaque site. 10 Conformément à la figure 1, des sites 1, 2, 3,..., j,... sont répartis géographiquement, en chacun desquels se produit une consommation d'énergie électrique. j est un indice entier naturel qui repère les sites individuellement. Chaque site peut être une habitation pavillonnaire, un immeuble d'habitat collectif, un site d'activité économique, par exemple un site 15 industriel ou un immeuble de bureau, etc. Il est équipé d'un compteur qui mesure la consommation d'énergie électrique qui a lieu en ce site pendant des durées déterminées successives. Ainsi, le compteur 1, mesure la consommation électrique du site 1, le compteur 2, la consommation électrique du site 2, le compteur 3, celle du site 3,..., le compteur j, celle du site j,... Les 20 compteurs 1,, 2 ,, 3 ,,..., j a,... ont été appelés moyens de collecte de valeurs de consommation dans la partie générale de la présente description, et peuvent être de modèles variables entre des sites 1, 2, 3,..., j,... qui sont différents. Les résultats de mesure de consommation sont transmis à partir de chaque site 1, 2, 3,..., j,... à une unité de traitement de données 100. Des 25 moyens de transmission très variés peuvent être mis en oeuvre, tels que la transmission automatisée par réseau de communication, de type téléphonique ou internet, la lecture et le relevé des résultats de mesure de consommation par un opérateur qui se déplace jusqu'à certains sites, etc., en fonction des modèles de compteurs dont les sites sont équipés. De façon symbolique, la 30 référence lt désigne les moyens de transmission qui sont mis en oeuvre pour transmettre les résultats de mesure de consommation du site 1 à l'unité de traitement de données 100, que ces moyens soient automatisés ou humains, ou une combinaison de moyens partiellement automatisés et partiellement 3046280 -13- humains. De même, la référence 2t désigne les moyens de transmission qui sont mis en oeuvre entre le site 2 et l'unité 100, les moyens de transmission 3t pour le site 3, etc..., ceux notés jt pour le site j,... Avantageusement, chaque site 1, 2, 3,..., j,..., ou le compteur 1,, 2 ,, 3 ,,..., j ,,... de ce site, peut être 5 pourvu de moyens d'enregistrement ou de mémorisation (non représentés), pour enregistrer les valeurs collectées à ce site avant de les transmettre. Un tel enregistrement intermédiaire permet de transmettre ou retransmettre les valeurs collectées à l'unité 100 après un arrêt ou un dysfonctionnement des moyens de transmission 1t, 2 t, 3 j t,... sans que les données soient 10 perdues. Alternativement ou en combinaison, un historique de consommation peut être enregistré ou mémorisé localement au niveau de chaque site 1, 2, j,..., pour être transmis ensuite à l'unité de traitement de données 100, si nécessaire. L'unité de traitement de données 100 peut être un ordinateur, usuel ou 15 avec une capacité de calcul augmentée. Les résultats de mesure de consommation pour chaque site 1, 2, 3,..., j,... peuvent se présenter sous deux formes différentes, notamment en fonction du modèle du compteur de ce site. La première forme est un index qui augmente en fonction du temps de façon à compter continûment la 20 consommation d'énergie électrique qui se produit au site considéré. Ainsi, l'index augmente d'autant plus vite que la puissance électrique instantanée est élevée. La valeur de l'index à un instant donné est donc une valeur de cumul de consommation. Dans le cas de la consommation d'énergie électrique, l'index peut être exprimé en unité de kilowattheure (kWh) par exemple. La 25 valeur d'index pour le site j à l'instant de relevé i est noté mi,i. La seconde forme possible pour un résultat de mesure de consommation, pour un site, est un comptage de la consommation nouvellement produite depuis le relevé précédent qui a été effectué pour le même site. Autrement dit, le comptage est réinitialisé à la valeur nulle à chaque 30 relevé, et le résultat d'une nouvelle mesure de consommation, pour le site j, est une valeur de la progression de consommation qui a eu lieu depuis l'instant du relevé précédent pour le même site j. 3046280 -14- On construit alors une matrice des progressions de consommation, qui est notée V et qui regroupe toutes les valeurs possibles des progressions de consommation qui sont relatives à un ensemble déterminé de sites. Une colonne différente de la matrice V est affectée à chaque site j, et une ligne 5 différente de la matrice V est affectée à chaque instant i auquel une valeur de cumul ou de progression de consommation a été collectée pour l'un quelconque des sites de l'ensemble déterminé. Alors, la valeur de progression de consommation vi,; est située dans la matrice V à l'intersection de la ligne qui est dédiée à l'instant de relevé i et de la colonne qui est dédiée au site j. La 10 matrice V des progressions de consommation possède donc autant de lignes que le nombre d'instants différents auxquels au moins un relevé de valeur de consommation - de cumul ou de progression - a été effectué à l'un quelconque des sites. Ce nombre est noté Nlignes sur la figure 2a. Simultanément, la matrice V possède autant de colonnes que le nombre, noté 15 Ncolonnes, de sites différents qui sont considérés. Etant donné que les consommations de tous les sites ne sont pas collectées simultanément, de nombreuses valeurs vi,; de la matrice V ne sont pas disponibles. Autrement dit, la matrice V est lacunaire. L'invention a notamment pour but de combler les valeurs vi,; manquantes dans la matrice V, par des valeurs qui sont estimées 20 pour la progression de consommation de chaque site j entre les instants i et i+1. Ces valeurs estimées sont déduites de celles des valeurs de progressions de consommation qui sont collectées réellement. Les valeurs de progressions de consommation qui sont collectées pour l'un quelconque des sites, entre deux instants qui sont successifs parmi tous 25 les instants de relevés pratiqués sur l'ensemble des sites, sont directement des valeurs v;,; de la matrice V des progressions de consommation. D'autres valeurs qui peuvent être collectées, dont celles désignées ci-après par mit Mi2, j, Vil/i2, j, v;1/2, J, et des données exogènes aux consommations, ne constituent pas directement des valeurs v;,; de la matrice V des progressions 30 de consommation, mais elles constituent des conditions que ces valeurs vi,; de la matrice V doivent satisfaire. Les valeurs de cumuls de consommation m;,;, appelées valeurs d'index, qui sont collectées sont reliées aux valeurs de la matrice V des 3046280 -15- progressions de consommation par l'équation suivante : i=i2 -1 J = M12 Jmi1,J (1) Autrement dit, la somme des progressions de consommation pour un même site j entre les instants de relevés il et i2 est égale à la variation d'index entre 5 ces deux instants. L'équation (1) constitue autant de premières conditions distinctes que doivent satisfaire les valeurs vi,; de la matrice V, que de couples de valeurs d'index qui sont relevées à un même site. De telles conditions ont été appelées contraintes dans la partie générale de la présente description. Un deuxième type de contraintes peut résulter de la collecte d'une 10 valeur de progression de consommation v;112,; pour un site j, entre deux instants de relevé il et i2 pour ce site qui ne sont pas successifs parmi tous les instants de relevés pratiqués sur l'ensemble des sites. Chaque contrainte de ce deuxième type est de la forme suivante : vi=i2 -1 JVil/i2,J (2) 15 L'équation (2) constitue autant de deuxièmes contraintes distinctes que doivent satisfaire les valeurs v;,; de la matrice V, que de couples de valeurs de progressions de consommation qui concernent un même site mais avec des instants de relevés qui ne sont pas successifs parmi tous les instants de relevés pratiqués sur l'ensemble des sites. 20 Un troisième type de contraintes peut résulter de la collecte d'une valeur de progression de consommation v;1/i2, J, qui est totalisée - ou agrégée pour un sous-ensemble J des sites j, entre les deux instants de relevé il et i2, ceux-ci n'étant pas nécessairement successifs parmi tous les instants de relevés pratiqués sur l'ensemble des sites. Chaque contrainte de ce troisième 25 type est de la forme suivante : vi=i2 -1 V i = V (3) .d1=11 jeJ 12 /12J L'équation (3) constitue autant de troisièmes contraintes distinctes que doivent satisfaire les valeurs vi,; de la matrice V, que de couples de valeurs de progressions de consommation qui concernent un même sous-ensemble de 3046280 -16- sites, et avec des instants de relevés qui ne sont pas nécessairement successifs parmi tous les instants de relevés pratiqués sur l'ensemble des sites. Ce troisième type de contraintes correspond, par exemple, à des relevés de consommation qui sont effectués à un noeud d'un réseau de distribution, qui 5 est dédié à l'alimentation en énergie électrique tous les sites du sous-ensemble J, correspondant par exemple à un quartier ou une zone géographique.Les contraintes des second et troisième types sont optionnelles pour l'invention. D'autres contraintes optionnelles peuvent résulter de données externes - ou exogènes - aux consommations. De telles données exogènes peuvent 10 provenir de mesures météorologiques, ou de considérations horaires ou calendaires, etc. Par exemple, les valeurs de progressions de consommation vi,; peuvent être reliées à la température externe par une équation du type : N (4) j=i 1,; j V E [f (TI) - a, f (TI) + b] Dans l'équation (4), J désigne un sous-ensemble de sites parmi tous les sites 15 qui sont considérés dans la matrice V, le sous-ensemble J contenant Nj sites, Ti est une température extérieure qui est mesurée sensiblement à l'instant de relevé i dans la zone des sites du sous-ensemble J, f est une fonction qui modélise la consommation électrique de chaque site du sous-ensemble J en fonction de la température Ti, et a et b sont deux nombres positifs. La fonction f 20 peut avoir été déterminée empiriquement, et les nombres a et b caractérisent la confiance sur la détermination de la fonction f. Conformément à la figure 2b, l'invention a pour but de fournir une estimation de la matrice V des progressions de consommation, à partir d'un produit de deux matrices qui ont chacune une dimension inférieure à la 25 dimension correspondante de la matrice V. Selon l'opération usuelle de produit matriciel, la première des matrices dans le produit, notée W, doit posséder Nlignes lignes, comme la matrice V. De même, la seconde des matrices dans le produit, notée H, doit posséder Ncolonnes colonnes, aussi comme la matrice V. En outre, le nombre de colonnes de la matrice W doit être égal au nombre de 30 lignes de la matrice H. Ce nombre est appelé rang dans le cadre de la présente invention, et noté K. Il est un majorant du nombre de variables qui déterminent la consommation de chaque site. Un intérêt important de l'invention est que le 3046280 -17- rang K soit le plus petit possible. En particulier, il doit être inférieur aux deux nombres de lignes et de colonnes, Niignes et Ncolonnes, de la matrice V. Par ailleurs, la matrice V des progressions de consommation est non- négative, c'est-à-dire que les valeurs v;,; sont chacune non-négatives. 5 Autrement dit, chaque valeur vi,; ne peut être que positive ou nulle, en tant que valeur de progression de consommation. Dans ce cas, et d'une façon qui est connue de l'Homme du métier, les deux matrices W et H peuvent elles-mêmes chacune être non-négatives, et leur produit constituer une approximation de la matrice V. On pourra consulter par exemple la référence «NeNMF: An Optimal 10 Gradient Method for Nonnegative Matrix Factorization» de Naiyang Guan et al., IEEE Transactions on Signal Processing, 60(6), pp. 2882-2898, June 2012 au sujet d'une telle décomposition de matrice non-négative. En outre, la matrice V satisfait les contraintes qui résultent des équations (1) à (4) présentées plus haut. Par conséquent, l'invention propose 15 d'assurer que l'estimation de la matrice V qui est obtenue à partir du produit matriciel W x H satisfasse ces mêmes contraintes en projetant le résultat du produit matriciel W x H sur le sous-ensemble des matrices de dimensions Nlignes x Ncoionnes qui satisfont toutes les contraintes. Du fait de la forme mathématique des contraintes, ce sous-ensemble est convexe et possède une 20 structure de sous-espace affine. Il est noté A, et la projection perpendiculaire sur ce sous-espace affine, à l'intérieur de l'ensemble des matrices de dimensions Niignes x Ncoionnes, est notée PA. Alors, l'approximation de la matrice V des progressions de consommation, qui assure que les contraintes soient satisfaites, est PA(W x H), comme indiqué à la figure 2c. PA(W x H) est appelée 25 matrice d'estimation dans la présente description. Le mode d'application de la projection PA est expliqué maintenant. A titre d'exemple, une contrainte du premier type, conforme à l'équation (1), est considérée. Les valeurs vi,; de la matrice de produit W x H qui sont concernées par cette contrainte, à savoir j fixé et i variant de il à i2-1, 30 sont additionnées. Si leur somme est inférieure à (mi2 - mit;), alors l'écart entre (mi2,; - mit ;) et la somme est réparti équitablement entre les valeurs vi qui participent à cette somme. Autrement dit, la matrice W x H est modifiée en 3046280 -18- augmentant les valeurs vi,; qui sont concernées par la contrainte, dans une mesure qui est identique pour toutes ces valeurs, jusqu'à ce que le total de ces valeurs atteigne le second membre de la contrainte qui résulte des valeurs collectées. A l'inverse, si la somme des valeurs v;,; de la matrice de produit 5 Wx H qui sont concernées par la contrainte, est supérieure à (mi2, - alors l'excédent de la somme par rapport à (mi2,; - mit ;) est retiré de chacune des valeurs vi,; qui participent à la somme. Ce retrait est d'abord effectué dans une mesure qui est identique pour toutes les valeurs v;,; concernées. Mais si l'une des valeurs vi,; ainsi réduite devient nulle ou négative, alors cette valeur 10 est fixée à zéro, et l'excédent de la somme par rapport à (mi2,; - mit ;), auquel est soustraite la valeur initiale de v;,; ramenée à zéro, est soustrait équitablement des autres valeurs vi,; de la somme concernée. Ce principe de retrait équitable de l'excédent, qui assure en plus qu'aucune valeur négative ne soit introduite dans la matrice de produit, est répété pour chaque valeur vi,; qui 15 viendrait à être ramenée à zéro. La nécessité d'éviter des valeurs négatives dans la matrice d'estimation est due à la nature positive de la matrice V des progressions de consommation. La collecte directe d'une valeur vi,; qui a été mesurée par le compteur j, du site j entre deux instants de relevés qui sont successifs dans la matrice V, 20 rentre dans le cas du premier type de contraintes avec i2 = il + 1. L'application de la contrainte revient alors à imposer cette valeur mesurée v;,; dans la matrice d'estimation à l'intersection de la ligne il et de la colonne j. La méthode d'application d'une contrainte du premier type, à la matrice de produit W x H, qui vient d'être expliquée est transposée aux contraintes des 25 autres types. Pour les contraintes du type de l'équation (4), si la somme des valeurs v;,; de la matrice de produit W x H qui sont concernées par la contrainte, est initialement dans l'intervalle [f(Ti)-a, f(Ti)+b], alors la contrainte est déjà satisfaite, et ne modifie pas ces valeurs vi,i. Si la somme des valeurs vi,; de la matrice de produit W x H qui sont concernées par la contrainte, est 30 initialement inférieure à f(Ti)-a, alors les valeurs v;,; qui participent à cette somme sont augmentées dans une mesure commune jusqu'à ce que leur somme atteigne la valeur f(Ti)-a. Et si à l'inverse la somme des valeurs vi,; de la matrice de produit W x H qui sont concernées par la contrainte, est 3046280 -19- initialement supérieure à f(Ti)+b, alors les valeurs vi,; qui participent à cette somme sont diminuées jusqu'à ce que leur somme atteigne la valeur f(Ti)+b, en empêchant comme précédemment que certaines des valeurs v;,; concernées deviennent négatives. 5 L'application de toutes les contraintes constitue la projection PA qui est appliquée à la matrice de produit W x H. Plusieurs méthodes peuvent être utilisées alternativement pour déterminer les matrices non-négatives W et H. La méthode qui est illustrée par la figure 3 est de type récursif : chaque itération d'une séquence répétée fournit 10 une nouvelle matrice d'estimation qui constitue une évaluation la matrice V, cette évaluation étant améliorée par rapport à la matrice d'estimation fournie par l'itération précédente. L'étape SO est une étape d'initialisation du procédé, et les étapes S1 à S4 constituent la séquence qui est répétée. t est un nombre entier qui est 15 incrémenté d'une unité à chaque exécution de la séquence. A l'étape SO, des matrices initiales non-négatives W° et H° sont fournies, avec un rang K qui est fixé par l'opérateur. Ainsi la matrice initiale Wo possède Nlignes lignes et K colonnes, et la matrice initiale Ho possède K lignes et Ncolonnes colonnes. Les évaluations de la matrice V des progressions de 20 consommation qui seront fournies par le procédé, constituées par les matrices d'estimation obtenues successivement, seront d'autant meilleures que le rang K n'est pas trop petit. Typiquement, le rang K peut être pris entre 5 et 20 lorsque Lignes est de l'ordre de 1000 et Ncolonnes de l'ordre de 600. Les valeurs des matrices initiales W° et H° peuvent être déduites de données de 25 consommation moyenne d'énergie électrique. Dans une première approximation, tous les sites peuvent être supposés avoir une consommation qui est constante dans le temps. Ainsi, toutes les valeurs de W° peuvent être prises égales à (A;)/K, où Ai est la durée entre les instants de relevés successifs i et i+1, et toutes les valeurs de H° pour la colonne j peuvent être 30 prises égales à la puissance électrique moyenne qui est consommée par le site j. Alternativement, les valeurs de W° et H° peuvent être sélectionnées aléatoirement selon une fonction de distribution du type loi normale. La matrice 3046280 - 20 - d'estimation pour l'initialisation du procédé est calculée comme PA(W° x H°), et notée V°, et le nombre t est initialisé à la valeur nulle. A l'étape S1, la matrice W est optimisée une première fois, à partir des matrices W°, H° et V°. Lors de cette première optimisation, les matrices H° et 5 V° sont maintenues constantes, et la matrice W est traitée comme une variable d'optimisation, à partir de W°, aboutissant à la matrice optimisée W1. Plusieurs procédés connus de l'Homme du métier peuvent être mis en oeuvre alternativement pour cette optimisation. Par exemple, le procédé du type à suivi de gradient optimal, pour «Optimal Gradient Method» en anglais, désigné 10 par l'acronyme OGM, peut être utilisé. Ce procédé d'optimisation est décrit par exemple dans l'article déjà cité «NeNMF: An Optimal Gradient Method for Nonnegative Matrix Factorization» de Naiyang Guan et al., IEEE Transactions on Signal Processing, 60(6), pp. 2882-2898, June 2012. Ainsi l'étape S1 consiste à calculer W1 comme le résultat de l'algorithme OGM appliqué aux 15 matrices W°, H° et V°, par rapport à la variable d'optimisation W : = OGM(W°, H°, V°). L'étape S2 correspond à l'étape S1, mais appliquée à la matrice H au lieu de W. Ainsi, lors de cette deuxième optimisation, les matrices W1 et V° sont maintenues constantes, et la matrice H est traitée comme une variable 20 d'optimisation, à partir de la matrice H°, aboutissant à la matrice optimisée H1. Le même procédé OGM qu'à l'étape S1 peut être utilisé pour l'étape S2. L'étape S2 consiste donc à calculer H1 comme le résultat de l'algorithme OGM appliqué aux matrices W1, H° et V°, par rapport à la variable d'optimisation H : H1 = OGM(W1, H°, V°). 25 L'étape S3 consiste à calculer la matrice V1 à partir des matrices W1 et H1 qui ont été obtenues respectivement aux étapes S1 et S2. La matrice V1 est obtenue à partir de la matrice de produit W1 x H1, en appliquant à cette dernière la projection PA de la façon qui a été décrite plus haut pour satisfaire toutes les contraintes. Ainsi, V1 = PA(W1 x H1). 30 Le nombre t est incrémenté d'une unité à l'étape S4. La matrice V1 peut alors être utilisée comme une première évaluation de la matrice V des progressions de consommation (étape S5). En ce sens, la 3046280 -21 - matrice V' est une première matrice d'estimation qui est fournie par le procédé de l'invention. Toutefois, la séquence des étapes S1 à S4 peut être répétée itérativement, à chaque fois à partir des matrices W, H et V qui ont été fournies par l'itération précédente. Ainsi, pour la (t+1)ème itération de la séquence, 5 l'étape S1 est constituée par le calcul Wt+1 = OGM(Wt, Ht, Vt) ) l'étape S2 par le calcul Ht+1 = OGM(wt-o, Ht, Vt), ) et l'étape S3 par le calcul Vt+1 = pA(Ari x Ht+i). L'étape S4 est encore l'incrémentation du nombre t. Le procédé peut être poursuivi jusqu'à un nombre d'itérations de la séquence des étapes S1 à S4 qui est fixé initialement, ou jusqu'à ce que les matrices Vt et Vt+1 obtenues lors 10 de deux itérations successives ne présentent plus d'écarts de valeurs significatifs, par exemple en termes d'influence sur la production d'électricité requise ou en termes de coût. Vt+1 est alors la matrice d'estimation qui résulte du procédé (étape S5), et elle constitue l'évaluation finale de la matrice V. Les valeurs v;,; qui ont été mesurées directement par les compteurs 1,, 2,, 15 3,,...,j,,... sont contenues sans altération dans la matrice d'estimation Vt+1, mais les valeurs manquantes dans la matrice V sont comblées par les valeurs correspondantes de la matrice d'estimation Vt+1. Par exemple, la matrice d'estimation peut contenir des valeurs de progressions de consommation pour tous les sites et pour des durées de 20 consommation successives toutes égales à une heure (consommations horaires, avec Ai égales à 1 heure pout tout i). A partir de la matrice d'estimation Vt+1 = pAmt+1 x Fris, ) il est possible d'évaluer la progression de consommation sur une durée continue, appelée période d'intégration, pour un site quelconque ou pour plusieurs sites, en 25 additionnant les valeurs concernées de la matrice d'estimation. Le résultat de l'addition a été appelé valeur d'estimation de la consommation totale associée au(x) site(s) concerné(s) pendant la durée d'intégration. Une qualité de la matrice d'estimation Vt+1 = pA(Ari x H`+') ) peut être évaluée en répétant la détermination de cette matrice à partir de sous- 30 ensembles variables des contraintes. Pour cela, parmi toutes les contraintes à disposition pour déterminer la matrice d'estimation, y compris celles des contraintes du premier type qui consistent à collecter directement certaines 3046280 - 22 - valeurs v;,; de la matrice V, un premier sous-ensemble de celles-ci est sélectionné, et la matrice d'estimation est calculée selon l'invention pour ce premier sous-ensemble de contraintes. Elle est ensuite calculée à nouveau, et indépendamment, pour au moins un autre sous-ensemble des contraintes. 5 Tous les sous-ensembles de contraintes qui sont utilisés successivement peuvent être sélectionnés de façon quelconque, de préférence aléatoirement, et éventuellement avec des recouvrements entre eux. Une valeur d'erreur moyenne ou d'écart-type peut alors être calculée à partir des matrices d'estimation obtenues, pour évaluer la qualité de l'estimation, notamment sa 10 fiabilité et sa stabilité par rapport à l'ensemble des valeurs qui ont été collectées. L'erreur moyenne ou l'écart-type peut être calculé pour n'importe quelle des valeurs d'estimation de la matrice Vt+1 = pAoAri x Fri) . Toutefois, il peut être préférable de calculer l'erreur moyenne ou l'écart-type pour une estimation de la consommation totale qui est associée à un ou plusieurs site(s) 15 déterminés, et pour une période d'intégration déterminée. Un tel calcul d'erreur moyenne ou d'écart-type utilise les valeurs d'estimation de la consommation totale pour ce(s) site(s) et la période d'intégration, qui sont calculées séparément à partir des matrices d'estimation issues des sous-ensembles différents de contraintes. 20 L'Homme du métier aura compris que l'invention fournit des estimations de consommation pour l'échelle spatiale la plus fine à l'intérieur d'une zone géographique, et pour l'échelle temporelle la plus fine pour chaque site de cette zone. Cette échelle temporelle correspond à tous les instants auxquels au moins un relevé a été effectué dans la zone géographique. Il est 25 alors possible de calculer des estimations de consommation pour des secteurs quelconques à l'intérieur de la zone géographique, et pour des périodes d'intégration quelconques. De tels calculs sont particulièrement utiles pour la consommation de secteurs à forte densité d'habitants, et pour des périodes de pointe de consommation. Ils permettent notamment de prévoir la mise à 30 disposition de l'objet de consommation en fonction de la demande. En particulier, la production de cet objet de consommation peut être ajustée, de même que son acheminement, par exemple en direction des secteurs à consommation élevée. Notamment, des délestages à partir de secteurs à 3046280 - 23 - consommation plus faible vers des secteurs à consommation élevée, ou des stockages intermédiaires de l'objet de la consommation, peuvent être organisés. Pour ces raisons notamment, l'invention est particulièrement avantageuse pour la distribution d'énergie électrique. 5 Il est entendu que des aspects secondaires de l'invention peuvent être adaptés ou modifiés par rapport à la mise en oeuvre particulière qui vient d'être décrite en détail, tout en conservant certains au moins des avantages cités. Parmi ces modifications possibles, les suivantes sont citées à titre non-limitatif : - la matrice d'estimation peut être déterminée en utilisant des procédés 10 différents de celui à convergence progressive qui est illustré par la figure 3. De tels procédés alternatifs sont disponibles dans la littérature du domaine technique de la factorisation par des matrices non-négatives ; - dans le procédé de la figure 3, l'algorithme d'optimisation à suivi de 15 gradient optimal (OGM) peut être remplacé par tout autre algorithme alternatif, connu de l'Homme du métier. Notamment, un tel autre algorithme peut aussi être du type algorithme d'optimisation, mais pas nécessairement ; et - enfin, encore dans le procédé de la figure 3, l'ordre chronologique de 20 calcul des deux matrices non-négatives peut être inversé. Dans ce cas, l'étape S1 consiste à calculer Ht+1 = OGM(Wt, Ht, Vt) en utilisant H comme variable d'optimisation, et l'étape S2 consiste à calculer Wt+1 = OGM(Wt, Ht+1, Vt) en utilisant W comme variable d'optimisation. L'étape S3 n'est pas changée. 25 REVENDICATIONS1. Procédé de gestion d'un approvisionnement en énergie, combustible, ou eau, ou en services délivrés localement, tels que des services d'enlèvement d'ordures, d'élimination d'effluents, de téléphonie, ou de transmission par un réseau de communication, pour des consommations se produisant en une multiplicité de sites géographiques (1, 2, 3,..., j,...) et pendant une période de surveillance, le procédé comprenant les étapes suivantes : /1/ collecter des valeurs de cumuls de consommation, associées chacune à un des sites (1, 2, 3,..., j,...) et à un instant de relevé de consommation compris dans la période de surveillance, et éventuellement collecter aussi des valeurs de progressions de consommation, associées chacune à un ou plusieurs des sites et correspondant à une consommation qui a été effectuée entre deux instants différents de relevés de consommation compris dans la période de surveillance ; /2/ transmettre les valeurs collectées à une unité de traitement de données (100) ; /3/ concevoir une matrice des progressions de consommation ayant des colonnes affectées une-à-une aux sites, et des lignes affectées une-à- une à des couples d'instants successifs, pour chacun desquels instants un relevé de consommation a été effectué pour l'un quelconque des sites, et destinée à contenir des valeurs de progressions de consommation attribuées chacune à un des sites pour une consommation qui a été effectuée audit site entre deux instants de relevés de consommation successifs, ladite matrice des progressions de consommation étant initialement lacunaire à cause de valeurs collectées manquantes pour les progressions de consommation de certains au moins des sites entre deux instants successifs de relevés de consommation ; 3046280 -25- /4/ exprimer les valeurs collectées sous forme de contraintes à satisfaire par les valeurs de progressions de consommation ; /5/ par l'unité de traitement de données (100) : déterminer deux matrices non-négatives qui fournissent, à partir d'une multiplication desdites 5 deux matrices non-négatives, une matrice d'estimation qui possède des mêmes nombres de colonnes et de lignes que la matrice des progressions de consommation, et dont des valeurs situées à des intersections de colonnes et de lignes de ladite matrice d'estimation satisfont les contraintes, avec la première des deux matrices non- 10 négatives dans la multiplication qui possède un nombre de colonnes égal à un nombre de lignes de la seconde matrice non-négative dans ladite multiplication, et qui est inférieur aux deux nombres de colonnes et de lignes de la matrice des progressions de consommation ; et /6/ utiliser au moins une valeur de la matrice d'estimation en tant 15 qu'évaluation d'une des valeurs de progressions de consommation non-collectées, pour gérer l'approvisionnement des sites. 2. Procédé selon la revendication 1, suivant lequel à l'étape /5/, les deux matrices non-négatives et la matrice d'estimation sont déterminées lors de répétitions d'une séquence qui comprend les sous-étapes suivantes, 20 effectuées en commençant à partir de deux matrices non-négatives initiales, et d'une matrice d'estimation initiale déduite du produit des deux matrices non-négatives initiales pour satisfaire les contraintes : /5-1/ en utilisant les deux matrices non-négatives et la matrice d'estimation telles que résultant de l'exécution précédente de la 25 séquence, modifier l'une des deux matrices non-négatives en maintenant constante l'autre matrice non-négative, de façon à réduire un écart entre le produit des deux matrices non-négatives et la matrice d'estimation ; 18,-2/ en utilisant la matrice non-négative modifiée à la sous-étape /5-1/, et 30 la matrice d'estimation telle que résultant de l'exécution précédente de la séquence, modifier ladite autre des deux matrices non-négatives en 3046280 -26- maintenant constante la matrice non-négative modifiée à la sous-étape /5-1/, de façon à réduire encore l'écart entre le produit des deux matrices non-négatives et la matrice d'estimation ; puis /5-3/ calculer le produit des deux matrices non-négatives obtenues à 5 l'issue des sous-étapes /5-1/ et /5-2/, respectivement, et modifier la matrice de produit pour satisfaire les contraintes, la matrice de produit modifiée formant la matrice d'estimation qui résulte de l'exécution de la séquence. 3. Procédé selon la revendication 2, suivant lequel chacune des sous- 10 étapes /5-1/ et /5-2/ est exécutée en utilisant un algorithme d'optimisation. 4. Procédé selon la revendication 3, suivant lequel l'algorithme d'optimisation est de type à suivi de gradient optimal. 5. Procédé selon l'une quelconque des revendications 2 à 4, suivant lequel à la sous-étape /5-3/, la matrice de produit est modifiée pour satisfaire 15 les contraintes en appliquant à ladite matrice de produit une projection orthogonale sur un sous-espace affine qui regroupe toutes les matrices ayant les mêmes nombres de colonnes et de lignes que la matrice des progressions de consommation, et qui satisfont les contraintes, la matrice de produit projetée formant la matrice d'estimation qui résulte de l'exécution de la séquence. 20 6. Procédé selon l'une quelconque des revendications précédentes, suivant lequel l'étape /2/ comprend un enregistrement intermédiaire d'une au moins des valeurs collectées à un des sites géographiques (1, 2, 3,..., j,...), sur un support d'enregistrement ou de mémoire disponible audit site géographique, puis une lecture et une transmission de ladite au moins une valeur collectée et 25 enregistrée, à destination de l'unité de traitement de données (100). 7. Procédé selon l'une quelconque des revendications précédentes, comprenant en outre les étapes suivantes : 3046280 -27- /7/ sélectionner un sous-ensemble de sites parmi la multiplicité des sites géographiques (1, 2, 3,..., j,...), et sélectionner une période d'intégration à l'intérieur de la période de surveillance ; puis /8/ produire une valeur d'estimation de la consommation totale associée 5 aux sites du sous-ensemble pendant la durée d'intégration, en additionnant les valeurs de la matrice d'estimation qui appartiennent à la fois aux colonnes affectées aux sites du sous-ensemble et aux lignes affectées aux couples d'instants successifs qui sont contenus dans la période d'intégration. 10 8. Procédé selon l'une quelconque des revendications précédentes, comprenant en outre les étapes suivantes, qui sont exécutées en utilisant la même multiplicité de sites (1, 2, 3,..., j,...) et les mêmes instants de relevés de consommation pour la matrice des progressions de consommation : /9/ sélectionner des groupes différents de valeurs collectées, parmi les 15 valeurs de cumuls de consommation associées chacune à un des sites (1, 2, 3,..., j,...) et à un instant de relevé de consommation compris dans la période de surveillance, et éventuellement aussi parmi les valeurs de progressions de consommation associées chacune à un ou plusieurs des sites et correspondant à une consommation qui a été 20 effectuée entre deux instants différents de relevés de consommation compris dans la période de surveillance ; /10/ pour chaque groupe de valeurs collectées, répéter les étapes /1/ à /6/ pour obtenir une matrice d'estimation indépendamment pour chaque groupe de valeurs collectées, les matrices d'estimation obtenues 25 correspondant toutes à la même matrice des progressions de consommation ; puis /11/ calculer au moins une valeur d'erreur moyenne ou d'écart-type à partir des matrices d'estimation obtenues pour tous les groupes de valeurs collectées. 3046280 - 28 - 9. Procédé les revendications 7 et 8, suivant lequel, l'étape /10/ comprend en outre l'étape /8/ répétée pour chaque groupe de valeurs collectées, pour le même sous-ensemble de sites et la même période d'intégration, tels que sélectionnés à l'étape /7/, et la valeur d'erreur moyenne 5 ou d'écart-type est calculée à l'étape /11/ pour les valeurs d'estimation de la consommation totale associée aux sites du sous-ensemble pendant la période d'intégration. 10. Procédé selon l'une quelconque des revendications précédentes, comprenant en outre de collecter des valeurs exogènes aux consommations, 10 en plus des valeurs collectées de cumuls de consommation et de progressions de consommation, et d'exprimer les valeurs exogènes collectées sous forme de contraintes supplémentaires à satisfaire par les valeurs de progressions de consommation, et suivant lequel les deux matrices non-négatives sont déterminées à l'étape /5/ de sorte que la matrice d'estimation satisfasse 15 simultanément les contraintes issues des valeurs collectées de cumuls de consommation et de progressions de consommation, et les contraintes issues des valeurs exogènes collectées. 11. Procédé selon l'une quelconque des revendications précédentes, appliqué à la consommation d'énergie électrique. 20 12. Procédé selon l'une quelconque des revendications précédentes, suivant lequel ladite au moins une valeur de la matrice d'estimation est utilisée à l'étape /6/ en tant qu'évaluation d'une des valeurs non-collectées de progressions de consommation, pour dimensionner une production, un acheminement, un délestage ou un stockage d'un objet de la consommation, 25 afin d'assurer une fourniture dudit objet de la consommation à l'un au moins des sites (1, 2, 3,..., j,...) en quantité suffisante par rapport à un besoin audit site. 13. Système de gestion d'un approvisionnement en énergie, combustible, ou eau, ou en services délivrés localement, tels que des services 30 d'enlèvement d'ordures, d'élimination d'effluents, de téléphonie, ou de 3046280 -29- transmission par un réseau de communication, pour des consommations se produisant en une multiplicité de sites géographiques (1, 2, 3,..., j,...) et pendant une période de surveillance, le système comprenant : - des moyens (1e, 2c, 3c, , jc, .) pour collecter des valeurs de cumuls de 5 consommation, associées chacune à un des sites (1, 2, 3,..., j,...) et à un instant de relevé de consommation compris dans la période de surveillance, et éventuellement aussi pour collecter des valeurs de progressions de consommation, associées chacune à un ou plusieurs des sites et correspondant à une consommation qui a été effectuée 10 entre deux instants différents de relevés de consommation compris dans la période de surveillance ; - des moyens de transmission (1 t, Z, 3t,..., jt,...) adaptés pour transmettre les valeurs collectées à une unité de traitement de données (100) ; et - l'unité de traitement de données (100) adaptée pour mettre en oeuvre au 15 moins les étapes /5/ et /6/ d'un procédé conforme à l'une quelconque des revendications précédentes. 14. Système selon la revendication 13, comprenant en outre des moyens d'enregistrement ou de mémorisation situés à l'un au moins des sites (1, 2, 3,..., j,...), et agencés pour réaliser un enregistrement intermédiaire d'une 20 au moins des valeurs collectées audit site, puis une lecture de la valeur collectée et enregistrée en vue de la transmission de ladite valeur collectée et enregistrée à destination l'unité de traitement de données (100).