FR 2936873 A1 20100409 FR 0856717 A 20081003 Procédé et dispositif de détermination non intrusive de la puissance électrique consommée par une installation, par analyse de transitoires de charge L'invention est relative à un procédé et à un dispositif de détermination non intrusive de la puissance électrique consommée par une installation, par analyse de transitoires de charge. A l'heure actuelle la puissance électrique consommée par une installation d'un client abonné, alimentée par le réseau en tension alternative, et la consommation d'énergie électrique correspondante, sont déterminées à partir de mesures directes des paramètres physiques de l'installation, tels que la tension électrique d'alimentation, l'intensité du courant électrique délivré, le déphasage entre le courant électrique délivré et la tension d'alimentation. De manière habituelle des modules électroniques permettent, à partir des valeurs de tension, de courant et de déphasage précitées, de calculer la puissance électrique active ou réactive consommée, ainsi que, par intégration de valeurs de puissance successives sur une plage temporelle donnée, l'énergie électrique consommée. Les modules électroniques précités sont le plus souvent implantés dans un appareil de comptage, tel qu'un compteur électronique et peuvent, dans certains cas, procéder à une transmission des valeurs précitées mesurées et/ou des valeurs de puissance ou d'énergie finalement consommée. Les appareils de comptage précités donnent satisfaction. Toutefois, ils présentent l'inconvénient de nécessiter l'implantation des modules électroniques précités dans les appareils de comptage eux-mêmes, et donc dans le périmètre privé de l'abonné client, c'est-à-dire, le plus souvent, l'habitation de l'abonné. Différents développements ont été proposés en vue de déporter la mesure et le comptage de la puissance ou de l'énergie hors du périmètre privé de l'abonné client. Une approche dans ce sens consiste à effectuer une mesure de tension d'alimentation et de courant délivré à l'extérieur du périmètre privé de l'abonné client, par analyse de la courbe de charge globale de l'installation soit au niveau de l'appareil de comptage, soit en amont de ce dernier, sur le branchement de l'installation de l'abonné client, dans le but de mettre en oeuvre des processus non intrusifs, connus sous le nom de NIALM pour Non Intrusive Appliance Load Monitoring en anglais. L'intérêt des processus non intrusifs précités est manifeste pour les fournisseurs, les distributeurs et les consommateurs d'énergie électrique face aux enjeux de prévision de coûts de l'énergie, de développements de réseaux et de réduction des consommations. Les processus de NIALM font apparaître plusieurs degrés de non intrusion, ceux comportant une étape d'apprentissage automatisé des appareils, décrits par le brevet US 4 858 141, et ceux comportant une étape d'apprentissage manuel des appareils, décrits par le brevet US 5 483 153 par exemple. Les processus NIALM manuels s'avèrent plus précis que les processus NIALM automatiques, grâce au recueil des signatures de consommation des appareils à différents états de consommation. Cette semi-intrusion est cependant gênante pour l'abonné client et peu attractive pour le distributeur ou le fournisseur d'énergie électrique. Dans le cas des processus NIALM manuels, une bibliothèque d'appareils récepteurs électriques est construite, à partir d'un capteur d'intensité de courant électrique sur chaque appareil. Dans le cas des processus NIALM automatiques une bibliothèque d'appareils récepteurs électriques s'améliorant dans le temps est construite, à partir des puissances mesurées sur l'appareil de comptage. Les bibliothèques précitées permettent en fait d'identifier chaque appareil électrique et sa consommation et ont favorisé une évolution vers l'identification des usages de ces appareils, voire des habitudes d'usage ou de consommation de ces derniers par les abonnés clients qui en disposent. Ainsi, différents modes opératoires d'identification des usages ont été proposés : à partir de la reconnaissance de la variation d'intensité du courant, ainsi que décrit dans la demande de brevet FR 2 645 968. Cette variation est comparée à un seuil représentatif d'un évènement, enclenchement ou déclenchement d'une charge domestique. Les valeurs de seuil sont répertoriées dans une bibliothèque ; à partir de la mesure des variations de puissance active et réactive et/ou de l'admittance de l'installation en régime permanent, ainsi que décrit par le brevet US 4 858 141. Une comparaison avec une table de référence de 25 30 5 10 15 20 30 différents appareils disponibles sur le marché, mémorisée en bibliothèque est effectuée ; à partir de la mesure du courant à la fréquence fondamentale et de ses harmoniques, ainsi que décrit dans le brevet US 6 816 078. Chaque appareil est identifié par les harmoniques (par transformée fréquentielle FFT) qu'il engendre et classé dans une bibliothèque. Le fonctionnement des appareils électriques est évalué à la valeur vraie par vérification de l'existence des fréquences harmoniques les plus hautes constitutives du courant électrique total ; à partir de la mesure des puissances actives et réactives de l'installation en régime transitoire, tel que décrit dans le brevet US 483 153. Une comparaison des formes des transitoires à une base de données est exécutée pour la reconnaissance des charges ; - à partir de l'exploitation de l'impulsion HF d'allumage émise par les appareils électriques reliés au réseau électrique tel que décrit dans les brevets EP 1 136 829 et US 7 078 982 B2. Ce mode opératoire nécessite, le cas échéant, sur chaque appareil à détecter, l'implantation d'un dispositif d'émission de signaux HF sur le réseau public, des récepteurs ou répéteurs pour amplifier les signaux HF. Bien que qualifié de non intrusif, ce mode opératoire nécessite plusieurs séries de mises sous et hors tension des appareils électriques l'un après l'autre, à la première installation du dispositif d'identification ; à partir des mesures effectuées, sur le courant et la tension pour déterminer la variation d'impédance des charges de l'installation au cours du temps, ainsi que décrit dans la demande de brevet WO 93/04377. Un tel processus reste toutefois intrusif. Un composant de mémoire de masse doit être périodiquement remplacé dans l'appareil de comptage par un préposé. L'ensemble des processus précités présente un niveau d'intrusion peu compatible avec une diffusion de masse que représentent les parcs d'appareils domestiques. Les processus précités utilisent des paramètres et grandeurs plus conformes aux réseaux d'alimentation en tension purement sinusoïdale, telles que puissances apparentes, réactives, actives, valeurs efficaces de courant et de tension, et sont ainsi mis à mal, car peu adaptés, à la diffusion, sur les installations d'abonnés, de plus en plus nombreuses d'appareils électroniques générateurs de perturbations électriques et de bruit, se superposant à la composante fondamentale. Les processus précités les plus connus pour l'homme de l'art mettent en oeuvre, pour exécuter une décomposition des signaux, un traitement fréquentiel, du type transformée de Fourier rapide, désignée FFT en anglais. Ce type de traitement nécessite une fenêtre d'observation relativement importante du signal traité, afin d'assurer une discrimination efficace des différentes composantes du signal traité, fréquence fondamentale, à 50 Hz ou 60 Hz, et fréquences harmoniques, en vue d'un traitement satisfaisant des composantes de bruit précitées superposées à la composante fondamentale. Parallèlement aux méthodes d'analyse purement fréquentielle précitées, par transformée fréquentielle, du type transformée de Fourier, d'autres méthodes théoriques d'analyse de la réponse en transmission en temps réel, de systèmes soumis à un champ électromagnétique variable ont été proposées. Essentiellement développée dans le cadre de l'étude de la réponse de réflecteurs ou d'antennes excités par une onde électromagnétique supposée se propager par ondes planes, une méthode spécifique, désignée méthode du Pencil, a été proposée, laquelle permet, à partir de la réponse transitoire et permanente d'une cible à une excitation électromagnétique rayonnée, de déterminer les pôles et résidus de la réponse précitée par résolution d'un problème généralisé de valeurs propres. Pour une description d'une approche théorique de la méthode du Pencil précitée, appliquée à une cible formée par un fil conducteur électrique, on pourra utilement se reporter à l'article intitulé � Generalized Pencil-of- Function Method for Extracting Poles of an EM System from its Transient Response � publié par Yingbo Hua et Tapan K. Sarkar membres du Department of Electrical and Computer Engineering, Syracuse University, Syracuse, NY 13344-1240, IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, Vol.37, n°3, February 1989, p 229-234. Plus récemment, une étude comparative entre les performances comparées de résolution spectrale entre les techniques de traitement fréquentiel, telle que la transformée de Fourier, et la méthode du Pencil, intitulée � Comparison between the Matrix Pencil Method and the Fourier Transform Technique for High Resolution Spectral Estimation � publiée par José Enrique Fernandez del Rio et Tapan K. Sarkar Department of Electrical and Computer Engineering, 121 Link Hall, Syracuse University, Syracuse, New York 13244-1240, Digital Signal Processing 6, 108-125 (1996) Article N°0011, a montré la supériorité de la méthode du Pencil vis-à-vis des méthodes de traitement fréquentiel par transformée de Fourier, en matière de dispersion d'estimation et d'erreur quadratique moyenne, pour un rapport signal à bruit supérieur à une certaine valeur de seuil. La présente invention a pour objet la mise en oeuvre d'un procédé et d'un dispositif de détermination non intrusive de la puissance électrique consommée par une installation, par analyse de transitoire de charges par la méthode du Pencil, l'installation étant assimilée à une ou plusieurs charges, pouvant présenter des caractéristiques de réponse variables dans le temps, absorbant partiellement et/ou totalement l'énergie électrique délivrée sous forme d'un signal d'excitation électromagnétique par le réseau de distribution. Le procédé de détermination non intrusive de la puissance électrique consommée par une installation d'abonné, conforme à l'objet de la présente invention, est exécuté à partir d'un échantillonnage périodique des valeurs de tension d'alimentation et de l'intensité du courant électrique délivrés à cette installation pour obtenir un signal de consommation. Il est remarquable en ce que, à partir des valeurs échantillonnées du signal de consommation, il consiste au moins à grouper ces valeurs échantillonnées sur des fenêtres temporelles d'observation successives, discriminer, sur chaque fenêtre temporelle d'observation courante successive, l'ensemble des pôles et des résidus 8 du signal de consommation par la méthode du Pencil. Un tel ensemble de pôles et résidus est représentatif, sur la fenêtre d'observation courante, d'une combinaison de signatures de charges électriques distinctes incluant au moins des pôles et des résidus déterminés associées à au moins une valeur singulière obtenue par la méthode du Pencil en fonctionnement dans cette installation, calculer sur au moins la durée de la fenêtre d'observation courante, au moins la puissance active consommée par au moins un sous-ensemble des charges électriques en fonctionnement dans l'installation, exprimée comme la somme des puissances actives consommées par chaque charge électrique distinct de ce sous-ensemble vérifiant une fonction de puissance consommée liant ces pôles et résidus. Le procédé objet de l'invention est également remarquable en ce que pour la puissance active, la fonction de puissance vérifie la relation : M M' (Sci+5 i,~ +4t r At 2 j P = E E Rci. Rvi' .e . sinh c (sci + svi, ) ù i=1 i'=1 où P : désigne la puissance active consommée par l'ensemble des charges électriques distinctes Ci de signature (SZ,RZ,SGZ) ; S,etR, : désignent les pôles respectivement les résidus du courant échantillonné délivré à l'installation ; Svl, et Rvl, : désignent les pôles respectivement les résidus de la tension échantillonnée délivrée à l'installation SGZ désignant les valeurs singulières; tkd désigne l'instant de début de la fenêtre d'observation de rang k; At : désigne la durée de la fenêtre d'observation sinhc(X) désigne la fonction sinus hyperbolique du cardinal de la valeur X, soit sinhc( .+S.,)Ot=sinh~(S.+s )Ots +s )Ot (Sei VZ 2 (Sei Se 2 Sei VZ 2 M : désigne le nombre de pôles et résidus observés du courant échantillonné délivré à l'installation ; M' : désigne le nombre de pôles et résidus observés de la tension échantillonnée délivrée à l'installation. Nem : désigne le nombre minimum d'échantillons sur la durée de chaque fenêtre d'observations, Nc désigne le 20 nombre observé de charges électriques distinctes de l'installation d'abonné, nci désigne le nombre de pôles de la charge électrique considérée, H désigne le nombre de composantes harmoniques de la fréquence fondamentale de la tension d'alimentation. 25 Le procédé objet de l'invention est en outre remarquable en ce que, pour une fenêtre d'observation du signal de consommation incluant le régime libre et le commencement du régime forcé d'enclenchement/déclenchement d'une charge, l'enclenchement/déclenchement d'une charge Le procédé objet de l'invention est également remarquable en ce que la fréquence d'échantillonnage, pour obtenir le signal de consommation, et la durée At de chaque fenêtre d'observation sont liées par la relation : At.fe=Nem2. Lnci+2(1+H)~ dans laquelle Z=i purement résistive est signé par un régime forcé unique, par l'existence, pour la fenêtre d'observation, de deux pôles complexes conjugués à partie réelle nulle et de deux résidus complexes conjugués dans le régime forcé, et par la phase nulle des résidus en régime forcé, l'enclenchement/déclenchement d'une charge résistive et capacitive est signé par un régime libre et un régime forcé distincts, par l'existence, pour la fenêtre d'observation, de deux pôles complexes conjugués à partie réelle nulle et de deux résidus complexes conjugués à partie réelle nulle dans le régime forcé et par un pôle réel et un résidu dans le régime libre, la phase négative du résidu en régime forcé caractérisant la nature globalement capacitive de la charge, l'enclenchement/déclenchement d'une charge résistive et inductive est signé par un régime libre et un régime forcé distincts, par l'existence, pour la fenêtre d'observation, de deux pôles complexes conjugués à partie réelle nulle et de deux résidus complexes conjugués dans le régime forcé et par un pôle réel et un résidu complexe dans le régime libre, la phase positive du résidu caractérisant la nature globalement inductive de la charge. Le procédé objet de l'invention est enfin remarquable en ce que, pour une fenêtre d'observation du signal de consommation incluant le régime libre et le début du régime forcé d'enclenchement/déclenchement d'une charge, l'enclenchement/déclenchement d'une charge résistive, inductive et capacitive est signé par un régime forcé et un régime libre distincts, par l'existence, en régime forcé, de deux pôles complexes conjugués à partie réelle nulle et de deux résidus complexes conjugués, et, en régime libre, de deux pôles complexes conjugués et deux résidus complexes conjugués, la phase négative du résidu en régime forcé caractérisant la nature globalement capacitive de la charge et la phase positive du résidu en régime forcé caractérisant la nature globalement inductive de la charge. Le dispositif de détermination non intrusive de la puissance électrique consommée par une installation d'abonné, objet de l'invention, comporte au moins des ressources d'échantillonnage périodique des valeurs de tension d'alimentation et de l'intensité du courant électrique délivrés à cette installation pour obtenir un signal de consommation. Il est remarquable en ce qu'il comporte en outre au moins un module de calcul de la puissance instantanée délivrée produit des valeurs échantillonnées des valeurs de tension d'alimentation et de l'intensité du courant électrique délivrés, le signal de consommation comportant l'une au moins des valeurs échantillonnées de la tension d'alimentation, de l'intensité du courant électrique délivrés ou de la puissance instantanée délivrée, des ressources de mémorisation et de lecture des valeurs successives du signal de consommation sur au moins une fenêtre d'observation, des ressources d'application, sur le signal de consommation, sur au moins une fenêtre d'observation courante, d'un traitement de filtrage par la méthode du Pencil, pour discriminer, sur cette fenêtre, l'ensemble des pôles et résidus représentatifs d'une combinaison de signatures de charges électriques distinctes incluant au moins les pôles et les résidus associés à au moins une valeur singulière obtenue par la méthode du Pencil en fonctionnement dans l'installation, un module de suivi temporel de la combinaison de signatures et de calcul, sur au moins la durée de la fenêtre d'observation courante, d'au moins la puissance active consommée par au moins un sous-ensemble des charges électriques en fonctionnement dans l'installation, exprimée comme la somme des puissances actives consommées par chaque charge électrique distincte de ce sous-ensemble vérifiant une fonction de puissance consommée liant ces pôles et résidus. Le dispositif objet de l'invention est en outre remarquable en ce que le module de suivi temporel de la combinaison de signatures et de calcul d'au moins la puissance active consommée par au moins un sous-ensemble de charges électriques comprennent au moins, un module de détection de changement d'état d'enclenchement/déclenchement de charge électrique, un module d'identification des charges électriques enclenchées et/ou déclenchées à partir des signatures de charges électriques, un module de librairie de charges électriques identifiées en relation avec leur signature, un module d'identification des usages de consommation de l'utilisation abonné et de calcul de la consommation, au moins en puissance active. Le procédé et le dispositif objets de l'invention trouvent application à la gestion non intrusive des installations domestiques ou industrielles d'appareils électriques, à la construction d'appareils de comptage d'énergie électrique consommée, et, plus généralement, à la surveillance et au contrôle des installations électriques d'usagers. Ils seront mieux compris à la lecture de la description et à l'observation des dessins ci-après dans lesquels : la figure la représente, à titre purement illustratif, un organigramme des étapes essentielles permettant la mise en oeuvre du procédé objet de l'invention ; - la figure lb représente, à titre illustratif, une variante de mise en oeuvre préférentielle non limitative du procédé objet de l'invention, dans laquelle un choix spécifique de la durée de la fenêtre d'observation peut être effectué en fonction de la fréquence d'échantillonnage du signal de tension et/ou de courant ; - la figure 2a représente, à titre purement illustratif, un exemple de mise en oeuvre non limitatif d'une étape d'acquisition des signatures de charges observées dans l'installation d'un abonné client ; - la figure 2b représente, à titre purement illustratif, un exemple de mise en oeuvre non limitatif d'une étape d'identification des charges électriques, à partir de leur signature ; - la figure 2c représente, à titre purement illustratif, un exemple de mise en oeuvre non limitatif d'exécution d'une étape de suivi de mise en fonctionnement/respectivement hors fonctionnement de chaque charge à partir d'une discrimination de leur signature. - la figure 3a représente, à titre d'exemple purement illustratif, le schéma électrique d'une installation d'abonné arbitraire permettant de représenter la mise en oeuvre du procédé objet de l'invention sur une telle installation ; 30 la figure 3b représente à titre d'exemple les valeurs d'intensité du courant mesuré pour une commutation des charges de l'installation de la figure 3a, par actionnement des interrupteurs T1, T2 , T3 et T4 selon le diagramme temporel représenté sur la figure 3b ; la figure 3c représente la distribution et le suivi temporel correspondant des valeurs singulières révélées par la méthode du Pencil, lors de la commutation des charges de l'installation représentée en figure 3b ; les figures 3d et 3e représentent un chronogramme du suivi temporel des valeurs réelles des pôles respectivement des valeurs imaginaires des pôles, pour la commutation des charges de l'installation représentée en figure 3b ; les figures 3f et 3g représentent un chronogramme du suivi temporel des résidus liés au régime forcé respectivement au régime libre pour la commutation des charges de l'installation représentée en figure 3b ; les figures 4a et 4b représentent à titre illustratif, un schéma fonctionnel, sous forme d'un schéma blocs, d'un dispositif de détermination non intrusive de la puissance électrique consommée sur une installation, conforme à l'objet de la présente invention. Une description détaillée des étapes de mise en oeuvre du procédé de détermination sans intrusion de la puissance 30 électrique consommée par une installation, par analyse de transitoires de charge, conforme à l'objet de l'invention sera maintenant donnée en liaison avec la figure la et 10 15 20 25 les figures suivantes. La mise en oeuvre sans intrusion du procédé et du dispositif objets de l'invention s'entend d'une exécution de mise en oeuvre en totale transparence de fonctionnement vis-à-vis de l'installation de tout client, soit en aval et/ou dans tout dispositif de comptage de ce client, soit en amont de ce dernier. En référence à la figure la précitée, on considère une installation d'abonné, installation domestique ou industrielle, pour laquelle on dispose de valeurs de tension d'alimentation et de l'intensité du courant électrique délivrées à cette installation, obtenues par échantillonnage périodique, pour obtenir un signal de consommation. Sur la figure la, les valeurs échantillonnées de tension Vx et d'intensité lx du courant délivré sont notées selon la relation (1) : KN Vx,1 1 fe Dans cette relation, x désigne le rang d'échantillonnage successif des valeurs précitées, KN désigne une valeur finale quelconque du rang d'échantillonnage sur K fenêtres d'échantillonnage successives comportant chacune N échantillons, fe désigne la fréquence d'échantillonnage des valeurs précitées. Selon un aspect remarquable du procédé objet de l'invention, celui-ci consiste, à partir des valeurs échantillonnées du signal de consommation, en une étape A, à grouper les valeurs échantillonnées précitées sur 16 des fenêtres temporelles d'observation successives. À l'étape A de la figure la, l'opération correspondante est notée selon la relation (2): TT T K1V k=K j=N {V x l x }1 [Vki,Pk/I kj k=1 j=1 Par l'opération correspondante, on comprend que tout ensemble de valeurs échantillonnées est en fait subdivisé en une pluralité de sous-ensembles de valeurs échantillonnées, chacun constitutif d'une fenêtre d'observation de ces valeurs échantillonnées, dans lesquelles k désigne le rang de la fenêtre d'observation et j désigne le rang de la valeur échantillonnée de tension respectivement de courant délivrés dans la fenêtre d'observation de rang k considéré. En outre K désigne la valeur finale du rang de la dernière fenêtre d'observation retenue et N désigne la valeur finale du rang du dernier échantillon de tension respectivement de courants dans chaque fenêtre d'observation de rang k considéré. D'une manière générale, on indique que le regroupement des échantillons successifs de valeur de tension respectivement de courant délivrés peut être effectué par adressage sélectif en lecture-écriture des valeurs d'échantillons de tension et de courant délivrés en mémoire, de manière classique. Les fenêtres successives d'observation ainsi obtenues peuvent avantageusement former des fenêtres sans recouvrement, l'application de gabarit de pondération et/ou d'observation des échantillons n'étant pas nécessaire. Chaque fenêtre d'observation formée par une suite de N échantillons de valeur de tension respectivement de courant délivrés correspond à une durée d'observation At proportionnelle au nombre N d'échantillons. L'étape A de groupement précitée est alors suivie d'une étape B consistant à discriminer, sur chaque fenêtre temporelle d'observation courante, de rang k successif, au moins l'ensemble des pôles notés Si et des résidus Ri du signal de consommation par la méthode du Pencil. À l'étape B de la figure la, l'opération de discrimination des pôles et des résidus est notée par la relation (3) . Méthode du Pencil K J [Vki, Pkj I ki l k=1 J 1 N [(s,R,sG);=ft1 Dans la relation (3) précédente, J dénote l'application de la méthode du Pencil, par traitement de données appliqué aux échantillons ; Lv ,Piçi,I4j désigne chaque fenêtre d'observation de rang k de N échantillons ; Nc k=K kSi,lei,SGiL k=1 désigne l'ensemble des pôles Si, des résidus Ri et des valeurs singulières SGZ observés révélés par la méthode du Pencil sur la fenêtre d'observation de rang k considéré. Selon un aspect remarquable du procédé objet de l'invention, l'ensemble des pôles, des résidus et des valeurs singulières précités est représentatif, sur la fenêtre d'observation courante de rang k considéré, d'une combinaison de signatures M de charges électriques Ci distinctes. On indique que selon un aspect spécifique non limitatif du procédé objet de l'invention, les pôles Si, les résidus R. et les valeurs singulières SGZ de toute charge électrique Ci sont calculés sur l'une au moins des composantes de courant, de tension ou de puissance instantanée formant un signal de consommation délivrées à l'installation ainsi qu'il sera explicité ultérieurement dans la description. Chaque signature inclut au moins des pôles Si et des résidus Ri déterminés associés à au moins une valeur 10 singulière SGZ représentative d'une charge et de la commutation de celle-ci en fonctionnement ou hors fonctionnement dans l'installation. Ainsi, chaque signature est formée par un triplet de valeurs numériques représentatives de pôles, de résidus 15 et d'une valeur singulière, et représentée par la notation S (SZ,RZ,SGZ) . L'étape B est alors suivie d'une étape C consistant à calculer la puissance électrique consommée sur au moins la durée de la fenêtre d'observation courante précitée. 20 On comprend que, selon un mode de mise en oeuvre avantageux non limitatif du procédé objet de l'invention, le calcul de puissance précité peut être effectué au moins pour la puissance active consommée par au moins un sous-ensemble de charges électriques en fonctionnement 25 dans l'installation, pour la puissance réactive ou finalement pour tout type de puissance électrique consommée. Selon un aspect remarquable du procédé objet de l'invention, le calcul de la puissance consommée par au 30 moins un sous-ensemble de charges électriques est établi comme la somme des puissances consommées par chaque charge électrique distincte de ce sous-ensemble, vérifiant une fonction de puissance consommée liant les pôles et résidus précités. À l'étape C de la figure la, et à titre de seul exemple non limitatif pour la puissance active consommée, la fonction de puissance vérifie la relation (4): M M' (s .+svi, P = LERci.Rvi,.e i=1 i'=l+ot 2 j . sinh c~(sci + sj Ot ~ 2 où P : désigne la puissance active consommée par l'ensemble des charges électriques distinctes Ci de signature (SZ,RZ,SGZ) ; S,etR, : désignent les pôles respectivement les résidus du courant échantillonné délivré à l'installation ; Svl, et Rvl, : désignent les pôles respectivement les résidus de la tension échantillonnée délivrée à l'installation ; tkd : désigne l'instant de début de la fenêtre d'observation de rang k; At : désigne la durée de la fenêtre d'observation ; sinhc(X) : la fonction sinus hyperbolique du cardinal de la valeur X, soit sinhc(Se, =sinh7(Se, . +S .,)/ 7(Sei +S 2 )Ot . v~ 2 v~ 2 v~ M : désigne le nombre de pôles et résidus observés du courant échantillonné délivré à l'installation ; M' : désigne le nombre de pôles et résidus observés de la tension échantillonné délivrée à l'installation. La relation (4) exprime ainsi la fonction de puissance sur la fenêtre d'observation. Différents justificatifs théoriques qualitatifs de processus permettant la mise en oeuvre du procédé objet de l'invention seront donnés ci-après. D'une manière générale, on indique que le signal de consommation est formé par l'ensemble des valeurs échantillonnées de tension et ou de courant délivrés à l'installation, complété éventuellement par la valeur de puissance instantanée correspondant sensiblement au produit de la valeur échantillonnée de tension respectivement de courants délivrés à l'installation, ou à une combinaison de ces valeurs, ainsi qu'il sera explicité ultérieurement de manière plus détaillée dans la description. La mise en oeuvre de la méthode du Pencil, par traitement des données appliquée aux valeurs précitées, opère le calcul d'une matrice de Hankel à partir d'une transformation par la méthode désignée méthode SVD pour Single Value Decomposition en anglais. Les pôles Si ainsi révélés représentent la contribution à la réponse de la charge au sens temporel au signal de consommation. Les résidus Ri représentent la contribution à la réponse de la charge au sens échange ou absorption d'énergie au signal de consommation. Les valeurs singulières SGZ représentent la transition d'un état de consommation à un autre du fait de l'enclenchement/déclenchement, total ou partiel, d'une charge de l'installation. Les valeurs singulières précitées sont directement liées à la mise en évidence des valeurs propres de la matrice de Hankel, lors de l'application de la méthode du Pencil. A chaque fenêtre d'observation des valeurs échantillonnées, ou plus généralement du signal de consommation, correspond un ensemble de valeurs singulières. Le changement de valeur singulière représente le changement ou transition du signal de consommation, car les valeurs singulières ont la même valeur en état stationnaire, en l'absence de commutation totale ou partielle de charge. Les rebondissements des valeurs singulières sont liés à la valeur des pôles Si. L'apparition d'une valeur singulière sur la fenêtre d'observation ou sur une succession de fenêtre d'observation est sensiblement synchrone du changement de régime d'intensité du courant consommé par l'installation sur commutation simultanée ou non d'une ou plusieurs charges de l'installation. Enfin la discrimination des pôles et des résidus permet de discriminer le type de charge en fonctionnement, charge résistive, capacitive ou inductive, ainsi qu'il sera décrit de manière plus détaillée dans la description ultérieurement. D'une manière plus spécifique, on indique que, pour la mise en oeuvre du procédé objet de l'invention et pour une tension d'alimentation sensiblement périodique de pulsation w et de période T, la durée At de chaque fenêtre d'observation courante de rang k peut-être prise égale à la période T. Un tel choix de la durée de la fenêtre d'observation n'est pas indispensable. Toutefois, d'un point de vue granulaire de l'observation, le choix précité permet en quelque sorte une remise à l'état initial du réseau ou de l'installation et de la charge. D'une manière plus générale, la durée de chaque fenêtre d'observation peut être choisie compte tenu de la fréquence d'échantillonnage fe, selon la relation (5) . Nc Ot.fe=Nem>_2. Lnci+2(1+H) Dans la relation précédente, fe désigne la fréquence d'échantillonnage ; Nem désigne le nombre minimum d'échantillons sur la 10 durée de chaque fenêtre d'observation ; Nc désigne le nombre observé de charges électriques distinctes Ci de l'installation d'abonné ; nci désigne le nombre de pôles de la charge électrique Ci considérée; 15 H désigne le nombre de composantes harmoniques de la fréquence fondamentale de la tension d'alimentation. En référence à la figure lb, on indique que le choix de la durée At peut être effectué par exemple à partir de la relation précitée au cours d'une sous étape Al puis 20 exécution du groupement correspondant en une sous étape successive A2. Les étapes B et C peuvent alors être effectuées successivement. Une description plus détaillée d'étapes de discrimination et d'acquisition des signatures, d'identification des 25 charges électriques constitutives du sous-ensemble de charges électriques distinctes puis de suivi de la mise en fonctionnement/hors fonctionnement de chaque charge distincte, conformément au procédé objet de l'invention, sera maintenant donnée en liaison avec les figures 2a, 2b 30 et 2c respectivement. 23 D'une manière générale, on indique que les étapes précitées sont exécutées sur les valeurs numériques représentatives des résidus, des pôles et des valeurs singulières par des programmes de tri de ces valeurs numériques dont les organigrammes généraux sont représentés à titre purement illustratif sur les figures précitées. En référence à la figure 2a, suite à étape consistant à discriminer l'ensemble des pôles et des résidus du signal de consommation, le procédé objet de l'invention comporte une étape de discrimination et d'acquisition des signatures V; (SZ,RZ,SGZ) formées par un couple de pôles Si, résidus Ri auxquels est associé une valeur singulière SGZ . Une description plus détaillée de processus permettant la mise en oeuvre d'une étape de discrimination et d'acquisition des signatures 7i, d'identification des charges électriques Ci en fonctionnement sur une fenêtre d'observation, respectivement de suivi de mise en fonctionnement/hors fonctionnement de chaque charge distincte sera maintenant donnée en liaison avec les figures 2a, 2b respectivement 2c. En référence à la figure 2a, suite à l'étape consistant à discriminer l'ensemble des pôles Si et des résidus Ri du signal de consommation, l'étape de discrimination et d'acquisition des signatures Vi peut être illustrée selon le processus représenté en figure précitée, à i=Nc partir des triplets (Si,Ri,SGi) i=1 . L'ensemble des triplets précités étant mémorisé sous forme d'une liste par exemple, on procède en une étape 100 à la lecture de chaque valeur singulière SGZ dont on dispose dans l'ensemble des valeurs singulières SGZ mémorisé selon la relation symbolique (6) . SGI =S Sur réponse positive au test 100, c'est-à-dire pour la lecture de la valeur singulière SGZ, on attribue aux 10 triplets (SZ,RZ,SGZ) un code d'identification de signature noté S correspondant, à l'étape 101. Le processus est poursuivi pour la valeur singulière suivante d'indice 1+1 par retour par l'étape 102 à l'étape 100 précitée, tant qu'il existe une valeur 15 singulière non lue SGZ. Au contraire, sur réponse négative au test de lecture de l'étape 100, le processus de lecture de l'ensemble des valeurs singulières est interrompu et l'on dispose à l'étape 103 de l'ensemble des signatures Vi de charges 20 électriques Ci susceptibles d'apparaître par une mise en fonctionnement/hors fonctionnement totale ou partielle sur l'installation considérée. L'ensemble des signatures précitées peut alors être mémorisé dans une mémoire ou bibliothèque de signatures, 25 ainsi qu'il sera décrit ultérieurement dans la description. En référence à la figure 2b, l'étape d'acquisition des signatures précitées est alors suivie, avantageusement, d'une étape d'identification des charges électriques 30 constitutives de tout sous-ensemble de charges5 électriques distinctes en fonctionnement sur au moins une fenêtre d'observation. On comprend, en particulier, que tout ou partie de l'ensemble des charges électriques de l'installation considérée peut-être mis en fonctionnement simultanément ou non, selon la mise en fonctionnement ou hors fonctionnement individuellement de chaque charge électrique de l'installation selon la seule volonté de l'utilisateur. Ainsi, en référence à la figure 2b, on procède, à partir de l'ensemble des signatures mémorisées et bien entendu des triplets pôles, résidus valeur singulière associés à chacune de ces dernières à l'identification, en une étape 200, de l'existence sur l'installation des pôles, résidus et valeur singulière existants et observés en fonctionnement. Ce processus d'identification peut être réalisé par une comparaison des valeurs de pôles résidus et valeur singulière correspondant aux valeurs de triplets associés à la signature correspondante, avec un intervalle de confiance d'identification de ces valeurs. Cette opération est représentée à l'étape 200 par la relation (7) : i E [1, Nc]? Sur réponse positive à l'étape 200 précitée, à la discrimination de la signature Vi est associée une charge électrique Ci correspondante aux caractéristiques électriques spécifiques. On comprend, en particulier qu'une correspondance biunivoque peut être effectuée entre toute valeur de signature ou code de signatures S et les caractéristiques électriques de chaque charge électrique, grâce à la mise en oeuvre d'une bibliothèque de charges associée à la bibliothèque de signatures, ainsi qu'il sera décrit ultérieurement dans la description. Les étapes 200, 201 sont poursuivies par passage à la signature suivante par retour à l'étape 200 par l'intermédiaire de l'étape 202 symbolisée par la relation i=i+1, pour toutes les signatures de charges électriques susceptibles d'être mises en fonctionnement sur l'installation. Lorsque l'ensemble des signatures, et donc des charges électriques, susceptibles d'être mises en fonctionnement pour l'installation considérée a été soumis au test de l'étape 200, on dispose à l'étape 203 de l'ensemble des charges électriques en nombre Nc en fonctionnement sur l'installation, noté selon la relation (8): i=Nc Ci i=l Cet ensemble constitue un sous-ensemble de la totalité des charges électriques susceptibles d'être mis en fonctionnement pour l'installation considérée, sous-ensemble égal à tout ou partie de cette totalité. On comprend bien entendu que la correspondance biunivoque entre la valeur ou code de signature M;, et la charge électrique Ci correspondante est également établie avec 25 le triplet (SZ,RZ,SGZ) auquel ce dernier est associé. La détection de la présence ou de l'absence de chaque triplet précité observé sur l'installation, et en conséquence de la signature Met de la charge électrique Ci correspondante, permet, selon un aspect remarquable du 30 procédé objet de l'invention, d'exécuter une étape de suivi de la mise en fonctionnement/hors fonctionnement de chaque charge individuelle distincte par discrimination de tout ou partie de la signature de cette dernière sur au moins une fenêtre d'observation. On comprend en particulier que la discrimination précitée peut être exercée sur une ou plusieurs valeurs du triplet associé à cette charge, selon un processus de logique floue tel que décrit précédemment dans la description. Ainsi que représenté à titre d'exemple en figure 2c, pour toute charge électrique Ci du sous-ensemble précité on procède en une étape de test 300 à l'identification de la signature `;, par discrimination des valeurs de pôles, résidus, valeur singulière ou d'une combinaison de ces derniers selon la relation (9): 3 Ci ? Sur réponse positive au test 300, la charge électrique Ci est déclarée en fonctionnement � ON � à l'étape 301. Au contraire sur réponse négative au test 301, la charge électrique Ci est déclarée hors fonctionnement �OFF �. Le processus de suivi peut être exécuté sur une durée quelconque par retour au test de l'étape 300 par l'intermédiaire de l'étape 302, notée i=i+1. Un exemple de mise en oeuvre des processus précités, sur une installation représentée en figure 3a comportant une charge résistive pure RI, une charge résistive-capacitive R2C2, une charge résistive-inductive R3L3 et une charge résistive capacitive et inductive R4L4C4 est décrit ci-après. On mesure le courant et la tension instantanés par échantillonnage à l'entrée de l'installation électrique d'une habitation. Ces deux signaux analogiques sont 28 convertis en valeurs numériques à l'aide du CAN à la fréquence d'échantillonnage fe= 1 . A partir de ces deux Te mesures, on calcule la puissance instantanée consommée par les charges en présence, selon la relation précédemment mentionnée. Les interrupteurs Tlr T2, T3 et T4 commutent leur charge respective RI, R2C2, R3L3 et R4C4L4 selon la logique � ON �, � OFF � indiquée sur la figure 3b, On applique la méthode du Pencil sur chaque fenêtre 10 d'observation de la puissance, de la tension et du courant d'alimentation. Pour chaque fenêtre d'observation d'indice k, N échantillons du signal sont prélevés. La méthode du Pencil permet de filtrer le signal observé 15 du bruit généré par la chaîne de mesure et de conversion, voire d'éliminer une composante non pertinente du signal. Elle fournit pour chaque fenêtre d'observation courante donnée de rang k, les vecteurs de taille M des valeurs singulières SGIk, des pôles Slk , et des résidus Ru pour 20 le courant, la tension et la puissance. La mémorisation sélective des résultats successifs fournis par la méthode du Pencil permet d'exécuter le suivi temporel. La méthode du Pencil est appliquée sur chaque fenêtre 25 d'observation temporelle du courant. La durée At de la fenêtre est paramétrable, elle est dans l'exemple donné choisie égale à At=10ms dans le cas des figures 3c à 3g, soit une demie période pour une fréquence d'alimentation à 50 Hz. 30 Un exemple de suivi temporel des valeurs singulières est donné sur la figure 3c. Sur la figure 3c, l'axe des abscisses est gradué en secondes et l'axe des ordonnées en valeur relative d'amplitude de valeur singulière, sans dimension. L'enveloppe du courant d'entrée représentée en figure 3b se retrouve au niveau de l'évolution des valeurs singulières ; leurs variations mettent en évidence les changements d'états, de manière sensiblement synchrone. Le nombre de pôles et de résidus pertinents est issu d'un traitement approprié de ces valeurs singulières. Le suivi temporel des parties réelles et imaginaires de chaque pôle est illustré respectivement sur les figures 3d et 3e. Sur la figure 3d, l'axe des ordonnées est gradué en Népériens par seconde, Np/s, et, sur la figure 3e, l'axe 15 des ordonnées est gradué en Hertz. Selon un mode de réalisation préférentiel non limitatif, le suivi temporel des résidus est exprimé par celui de leurs valeurs normalisées, en régime forcé et en régime libre figures 3f et 3g. Ce choix, mieux approprié, 20 facilite la classification par rapport aux résidus. Sur les figures 3f et 3g, l'axe des ordonnées est gradué en Ampères (A). Ces suivis temporels délivrent des informations sur les caractéristiques des charges. 25 Les harmoniques, dont le fondamental, sont caractérisés par des pôles purement imaginaires, ils correspondent au régime forcé des charges. Les régimes libres sont caractérisés par des pôles dont la partie réelle donne des informations sur les 30 constantes de temps (confer figure 3d), et la partie imaginaire sur les fréquences propres (confer figure 3e). Les résidus sont directement liés aux charges en présence et aux harmoniques si ces fréquences existent. Selon un mode de réalisation préféré, pour le régime forcé, on définit une valeur normalisée pour chaque résidu Rlk. Cette définition peut être donnée pour tout type de résidu, de tension, de courant, de puissance et, notamment pour les résidus de courant, donnée par l'expression suivante avec, pour l'exemple ci-après, R,Z noté 9ZIk désignant le résidu de courant et S, noté Slk désignant le pôle de courant dans la fenêtre d'observation numéro k, 9ZINk désignant la valeur normalisée correspondante : n =9 . j . T . e 'ftkd l'instant du début de la fenêtre d'observation numéro k ; 9ZIk , le résidu du courant lié au régime forcé et relatif à la fenêtre k, pour la fréquence fondamentale et pour les fréquences harmoniques 20 si elles existent ; SIk , le pôle du courant lié au régime forcé et correspondant au résidu du courant RIk. La figure 3f donne l'évolution temporelle des résidus normalisés de courant 9ZIk liés aux fréquences positives. 25 Selon un mode de réalisation, l'exploitation des 9ZINk permet de remonter à la valeur de la charge. Par analogie, on peut établir une normalisation des résidus du courant liés au régime libre de la manière avec : tkd suivante, 9-UNI' désignant la valeur normalisée correspondante : 9ilNIcL ù 91 k . j . J . eûsl k •(tkd ût0 ) avec tkd , instant du début de la fenêtre d'observation numéro k ; to, instant d'enclenchement de la charge ; 9ZIk , désigne le résidu du courant lié au régime libre et relatif à la fenêtre k ; Slk , désigne le pôle du courant lié au régime libre, relatif à la fenêtre k et correspondant au résidu 9ZIk . La figure 3g donne l'évolution temporelle de 9ZINk, résidu normalisé du courant lié au régime libre. Sur 15 cette figure apparaissent trois instants d'enclenchement de charge, soit to correspondant aux instants 0,2s, 0,61s et 1,1s. Les pôles, figures 3d, 3e, liés à la fréquence positive 20 de l'excitation sont représentés par des cercles � o �. Ceux liés à la fréquence négative conjuguée, sont représentés par des croix � x �. Les pôles liés au régime libre de la charge sont représentés par des étoiles et des losanges � * , 0 �. 25 Les résidus normalisés, liés aux fréquences positives, sont représentés par des cercles � o � (figure 3f) dans le cas du régime forcé et par des étoiles � * � (figure 3g) pour le régime libre. 10 On comprend que le processus de suivi mis en oeuvre consiste essentiellement en une lecture des valeurs de pôles, résidus, valeurs singulières et en une comparaison de ces données pour détection d'un changement d'état. Les instants de changement d'état sont détectés au niveau de la variation des valeurs singulières représentées en figure 3c. Ces instants apparaissent aussi au niveau du suivi des résidus normalisés. La figure 3f illustre ce phénomène. L'interprétation des suivis temporels des pôles et des résidus normalisés donnés par les quatre figures 3d, 3e, 3f et 3g est la suivante : Entre les instants Os et 0.2s : Entre ces deux instants, les figures 3d et 3e donnent deux pôles imaginaires purs (0+50j) et (0û50j) qui correspondent à la fréquence du 5011z. On comprend que la forme générale des pôles normalisés est donnée par 0�17 où F désigne la fréquence fondamentale de la tension d'alimentation et j = . Si on considère qu'aucune charge n'est présente avant l'instant t=0, la variation du résidu lié au régime forcé, figure 3f, est de .eiN1 calo =(2.3069 -4.2406.10-6j P0) . Selon le mode de réalisation préféré, à partir de cette information on peut déduire le module de la charge ainsi que sa phase. A cet instant, on est en présence d'une charge de module 10052 et de phase nulle. Or, la valeur utilisée en expérimentation était une 30 résistance de valeur R1=10052 et la valeur théorique du résidu normalisé est :9-iNl Th =2.3+0j Entre les instants 0.2 et 0.6s : Entre ces deux instants, les figures 3d et 3e montrent l'existence de trois pôles, (0+50j), (0û50j) et (û10+CWW. Ce dernier pôle indique le déclenchement à l'instant 0.2s d'une charge du premier ordre de constante de temps i = 1 =100ms . 10 Selon le mode de réalisation préféré, on calcule la variation des résidus normalisés liés au régime forcé à l'instant 0.2s, fournis par la figure 3f : A99N2,ale _ (4.6022 +7.3130.10-2j)- (2.3069 - 4.2406.10-6 j)= 2.2953 +7.3134.10-2j Selon le mode de réalisation préféré et à partir de cette information, on peut alors calculer le module de la charge ainsi que sa phase. A cet instant on observe une charge de module 100.1552 et de phase -0.0318rd(û1.822°). On en déduit le caractère capacitif de cette charge. Or, la charge qui est enclenchée à l'instant précité est un circuit RC série de résistance R2 =10052 et de capacité C2=lmF. La valeur théorique du résidu normalisé correspondant est ~JZN2 Th = 2.2977 +7.313.10-2j , elle est identique à celle issue de la méthode du Pencil. Entre les instants 0.6s et 1.1s : Entre ces deux instants, les figures 3d et 3e donnent trois pôles, (0+50j) et(0û50j) ainsi qu'un nouveau pôle -50+0j) . Ce dernier pôle indique l'enclenchement à l'instant 0.6s d'une charge du premier ordre de constante de temps t=50=20ms. Selon le mode de réalisation préféré, on calcule la variation des résidus normalisés liés au régime forcé à l'instant 0.6s, fournis par la figure 3e : 49-iN3,a1 e = (5.168 -3.4970j)-(4.6022 +7.3130.10-2j)=0.56590 -3.5701j Selon le mode de réalisation préféré et à partir de cette information, on peut déduire le module de la charge ainsi que sa phase. A cet instant on est en présence d'une charge de module 63.6252 et de phase 1.413 rd (80.9589°) . On en déduit le caractère inductif de cette charge. Or, la charge enclenchée est un circuit RL série de résistance R3 =10n et d'inductance L3 =0.2H . La valeur théorique du résidu normalisé correspondant est 9ZN3 th =0.56820 -3.5701j , elle est identique à celle issue de la méthode du Pencil. Entre les instants 1.1s et 1.5s : Entre ces deux instants, les figures 3d et 3e donnent quatre pôles, (0+50j) et (0-50j) relatifs à la fréquence fondamentale ainsi que (-30+22j) et (-30-22j) relatifs à la fréquence propre de la charge. A partir de ces pôles, on peut reconnaître la présence d'une charge du second ordre à l'instant 1.1s (apparition de deux pôles complexes conjugués). A cet instant, selon le mode de réalisation préféré, on calcule la variation des résidus normalisés liés au régime forcé : A9Î.N4calc _ (5.5818 -5.2337j)û(5.168 -3.4970j)= 0.41370-1.7367j Selon le mode de réalisation préféré et à partir de cette information, on peut calculer le module de la charge ainsi que sa phase. A cet instant on a une charge de module 128.83Ç et de phase 1.3369 rd 76.5987°). On en déduit le caractère inductif de cette charge. Or, la charge enclenchée est un circuit RLCsérie de résistance R4 = 3052 , d'inductance L4 = 0.5H et de capacité C4 = 0.1 mF . Le résidu normalisé théorique correspondant est ~JZN4 th =0.41598 -1.7367j , il est identique à celui issu de la méthode du Pencil. Entre les instants 1.5s et 1.8s : La variation négative des valeurs singulières à l'instant 1.5s indique la mise hors tension d'une des quatre charges mises sous tension entre les instants Os et 1.5s d'une part. D'autre part, la variation des résidus normalisés du régime forcé est de 4~JZNS~a~~ =-2.2953 -7.31.10-2j . Selon le mode de réalisation préféré, on en déduit la mise hors tension de la charge n° 2, qui s'était enclenchée à l'instant 0.2s. Entre les instants 1.8s et 2s : La variation négative des valeurs singulières à l'instant 1.8s indique la mise hors tension d'une des quatre charges mises sous tension entre les instants Os et 1.8s d'une part. D'autre part, la variation des résidus normalisés du régime forcé est de 499N6calc =-2.2977 +0j . Selon le mode de réalisation préféré, on en déduit la mise hors tension de la charge n° 1 qui s'est enclenchée à l'instant Os. Ainsi, selon un aspect remarquable du procédé objet de l'invention, on comprend que la discrimination de l'ensemble des pôles Si, des résidus Ri et des valeurs singulières SGZ est effectuée par lecture des valeurs correspondantes révélées par l'application de la méthode du Pencil. En outre, dans le mode de mise en oeuvre préféré du procédé objet de l'invention, dans lequel les pôles et les résidus sont établis selon leur définition normalisée, le suivi de mise en fonctionnement/hors fonctionnement de chaque charge électrique individuelle Ci peut, de manière particulièrement avantageuse, être mis en oeuvre par simple calcul de la variation algébrique des résidus normalisés. En particulier l'apparition, enclenchement, d'une charge est révélée par l'apparition de la valeur algébrique du résidu Ri associé à cette charge et la disparition, déclenchement, de cette charge 20 est révélée par l'apparition de la valeur algébrique opposée du résidu Ri associé à cette charge. Un récapitulatif des caractéristiques des résidus, pôles et valeurs singulières, en définitive des signatures obtenues en fonction de caractéristiques électriques de 25 charges électriques enclenchées respectivement déclenchées sur une installation sera maintenant donné ci-après. On considère une fenêtre d'observation, ou une succession de fenêtres d'observation, du signal de consommation 30 incluant le régime libre et le commencement du régime forcé d'enclenchement/déclenchement d'une charge électrique. De manière caractéristique : l'enclenchement/déclenchement d'une charge purement résistive est signé par un régime forcé unique, par l'existence, pour la fenêtre d'observation considérée, de deux pôles complexes conjugués à partie réelle nulle et de deux résidus complexes conjugués dans le régime forcé. En outre, les résidus présentent une phase nulle en 10 régime forcé. Le tableau désigné Tableau charge résistive donne les valeurs des pôles et résidus normalisés correspondants. Tableau charge résistive : Régime forcé Pôles (0+50.j) (0-50.j) Résidus (x+y.j) (x- y.j) 15 L'existence de deux valeurs pour les pôles et les résidus, ainsi que la phase nulle des résidus normalisés liés au régime forcé permettent de déduire la nature résistive de la charge en cours d'analyse ; l'enclenchement/déclenchement d'une charge 20 résistive et capacitive est signé par un régime libre et un régime forcé distincts, l'existence pour la fenêtre d'observation de deux pôles complexes conjugués à partie réelle nulle et de deux résidus complexes conjugués dans le régime forcé et par un pôle réel et un résidu complexe 25 dans le régime libre. La phase négative du résidu en régime forcé caractérise la nature globalement capacitive de la charge. Tableau charge RC : Régime forcé Régime libre Pôles (0+50.j) (0û50.j) +O.j) Résidus (x+y.j) (xû y.j) (a+b.j) La phase des résidus normalisés liés au régime permanent est négative. La comparaison de la valeur absolue de la phase avec le produit a.w ou l'inverse de ce produit permet de déduire l'existence d'une charge RC en série ou d'une charge RC en parallèle. - l'enclenchement/déclenchement d'une charge résistive et inductive est signé par un régime libre et un régime forcé distincts, par l'existence pour la fenêtre d'observation de deux pôles complexes conjugués à partie réelle nulle et de deux résidus complexes conjugués dans le régime forcé et par un pôle réel et un résidu complexe dans le régime libre. La phase positive du résidu caractérise la nature globalement inductive de la charge. Tableau charge RL : Régime forcé Régime libre Pôles (0+50.j) (0û50.j) +O.j) Résidus (x+y.j) (xû y.j) (a+b.j) La phase des résidus normalisés liés au régime libre est positive. La comparaison de la valeur absolue de la phase avec le produit co .ou l'inverse de ce dernier permet de déduire l'existence d'une charge RL en série ou d'une charge RL en parallèle. - l'enclenchement/déclenchement d'une charge 25 résistive, inductive et capacitive est signé par un régime forcé et un régime libre distincts et par l'existence, en régime forcé, de deux pôles complexes conjugués à partie réelle nulle et de deux résidus complexes conjugués et, en régime libre, de deux pôles complexes conjugués et deux résidus complexes conjugués. La phase négative du résidu en régime forcé caractérise la nature globalement capacitive de la charge et la phase positive du résidu en régime forcé caractérise la nature globalement inductive de la charge. Tableau charge RLC : Régime forcé Régime libre Pôles (0+50.j) (0ù50.j) (î+.j) Résidus (x+y.j) (xùy.j) (a+b.j) (aùb.j) On rappelle que l'ensemble des valeurs de pôles et de résidus donné aux tableaux précédents est obtenu par l'application de la méthode du Pencil. En particulier, dans l'exemple des valeurs données, les valeurs numériques 50 affectées à la partie imaginaire des pôles correspond en fait à la valeur de la fréquence fondamentale de la tension d'alimentation de l'installation. Cette valeur est remplacée par la valeur 60 pour la mise en oeuvre du procédé objet de l'invention dans les pays dans lesquels les réseaux d'alimentation en énergie électrique sont des réseaux alternatifs à la fréquence 60 Hz. Une description plus détaillée d'un dispositif de détermination non intrusif de la puissance électrique consommée par une installation d'abonné, conforme à l'objet de la présente invention, sera maintenant donnée en liaison avec la figure 4a et la figure 4b. Ainsi qu'on l'a représenté sur la figure 4a précitée, le dispositif objet de l'invention comprend des ressources d'échantillonnage périodique des valeurs de tension d'alimentation et de l'intensité du courant électrique délivrés à l'installation considérée, en vue d'obtenir un signal de consommation. Sur la figure précitée, les ressources correspondantes sont réputées constituées par un circuit délivrant une image de la tension d'alimentation et un circuit délivrant une image de l'intensité du courant électrique délivrés à l'installation. Ces circuits peuvent avantageusement être constitués par tout circuit de prélèvement de tension respectivement de courant accompagnés d'échantillonneurs appropriés délivrant dans les conditions d'échantillonnage de Shannon une image de la tension respectivement de l'intensité du courant délivré. Ces circuits peuvent être suivis de convertisseurs analogiques et numériques à la fréquence ou période d'échantillonnage désignée l',. Ce type de circuit ne sera pas décrit en détail car il correspond à 20 des circuits connus de l'homme de l'art dans le domaine technique correspondant. Chaque circuit convertisseur analogique numérique CAN délivre avec un rythme correspondant à la période d'échantillonnage Te les échantillons de tension notés 25 v(nT) respectivement de courant i(nT ), pour bien entendu n>_0, n désignant le rang de l'échantillon. Selon une caractéristique particulièrement avantageuse, le dispositif objet de l'invention comporte une ressource de calcul 1 de la puissance instantanée délivrée, produit 30 des valeurs échantillonnées des valeurs de tension d'alimentation et de l'intensité du courant électrique délivré. Ainsi, sur la figure 4a, la ressource de calcul 1 reçoit les valeurs échantillonnées de courants i(nT), et de tension v(nT) et délivre une valeur produit de ces valeurs reçues, constitutive de la valeur de puissance 5 instantanée délivrée notée p(nT). Selon une autre caractéristique du procédé et du dispositif objets de l'invention, le signal de consommation est constitué par un au moins des échantillons d'intensité délivrée i(nT), de tension 10 délivrée v(nT) et de puissance instantanée délivrée p(nT à l'installation. Les composantes du signal de consommation précitées sont soumises à observation pour détermination des pôles, des résidus et des valeurs singulières grâce au dispositif objet de l'invention de 15 la manière ci-après. Le dispositif selon l'invention, représenté en figure 4a comporte dans ce but des ressources de mémorisation et de lecture des valeurs successives du signal de consommation sur au moins une fenêtre d'observation, les opérations de 20 mémorisation et de lecture étant bien entendu exécutées sur plusieurs fenêtres d'observation successive. Il comporte par exemple, ainsi que représenté en figure 4a, des ressources de création de fenêtres d'observation sur les composantes d'intensité, de tension et de 25 puissance précitées. Ces ressources sont notées 21r 22, 23 sur cette figure. Elles peuvent consister de manière classique en des circuits d'adressage en écriture/lecture en mémoire des composantes précitées mémorisées pour une pluralité de valeurs successives de rang n compris entre 30 (k -1)N et kN-1, N désignant le nombre d'échantillons dans la fenêtre de rang k considérée, en liaison avec le figure la. Le dispositif objet de l'invention représenté en figure 4a comporte en outre des ressources 3 d'application, sur le signal de consommation et sur au moins une fenêtre d'observation courante, tel que défini précédemment, d'un traitement de filtrage par la méthode du Pencil, ceci dans le but de discriminer sur la fenêtre du signal de consommation l'ensemble des pôles et résidus représentatifs, avec les valeurs singulières, d'une combinaison de signatures de charges électriques distinctes en fonctionnement dans l'installation considérée. À l'observation de la figure 4a, on constatera que les ressources 3 d'application du traitement de filtrage par la méthode du Pencil comportent avantageusement des modules séparés portant la référence 31r 32, 33, chacun recevant les composantes d'intensité, de puissance respectivement de tension du signal de consommation précitée. Chaque module séparé 31r 32, 33 peut être constitué par un calculateur dédié incorporant un programme de traitement approprié. Chacun des modules précités délivre les pôles, les résidus et les valeurs singulières pour la composante d'intensité, de puissance instantanée, de tension du signal de consommation, ces valeurs étant respectivement désignées SCk,Rck,SGCk pour la composante d'intensité, Spk,Rpk,SGpk pour la composante de puissance instantanée, SVk,RVk,SGVk pour la composante de tension délivrée à 30 l'installation. Les valeurs précitées ne sont autres, bien entendu, que les valeurs de pôles SZ , de résidus Ri et de valeurs singulières SGZ précédemment mentionnées dans la description. En outre, ainsi que représenté en figure 4a et 4b, le dispositif selon l'invention comprend des ressources 4 de suivi temporel de la combinaison de signatures et donc des charges électriques et de calcul sur au moins la durée de la fenêtre d'observation courante de rang k d'au moins la puissance active consommée par au moins un sous-ensemble de charges électriques en fonctionnement dans l'installation. La puissance calculée précitée est exprimée comme la somme des puissances, telles que la puissance active, consommée par chaque charge électrique distincte du sous-ensemble précité vérifiant une fonction de puissance consommée liant les pôles et les résidus. On comprend bien sûr que la fonction de puissance liant les pôles et les résidus, notamment pour la puissance active, vérifie la relation (4) précédemment mentionnés dans la description. On observera également sur la figure 4a que les ressources 4 de suivi temporel de la combinaison de signatures comportent également des modules séparés référencés 41r 42, 43 opérant sur les pôles, résidus, valeurs singulières de la composante d'intensité, de puissance respectivement de tension délivrée à l'installation. En particulier, le suivi temporel de la combinaison des signatures est exécuté conformément au processus décrit en liaison avec les figures 2a à 2c, en particulier par détection de la variation algébrique des valeurs des résidus et des pôles correspondants par exemple. En outre, les ressources 4 de suivi temporel sont maintenant décrites ci-après en détail en référence à la figure 4b. En référence à la figure 4b précitée, les ressources 4 de suivi temporel comportent, en outre, au moins un module 44 de détection de changement d'état d'enclenchement/déclenchement de charge électrique, lequel reçoit les valeurs numériques représentatives de la combinaison de signatures délivrées par chaque module séparé 41r 42 et 43 est directement relié à un module 45 d'identification des charges électriques enclenchées et/ou déclenchées à partir des signatures de charges électriques précitées. Le module 45 d'identification des charges électriques enclenchées et/ou déclenchées est lui-même relié à un module 46 de librairie de charges électriques identifiées en relation avec leurs signatures et à un module 47 d'identification des usages de consommation de l'utilisation abonné et de calcul de la consommation au moins en puissance active. Ainsi que représenté sur les figures 4a et 4b précitées, les modules 41à 47 sont bien entendu synchronisés par un signal noté k, N sur la lecture de chaque fenêtre d'observation de rang k comportant, pour un nombre N d'échantillons, un nombre donné de valeurs de résidus, pôles, valeurs singulières délivrées par chaque module respectif 31r 32, 33 relatif aux composantes d'intensité, de puissance instantanée respectivement de tension, ainsi que décrit précédemment dans la description. Chacun des modules précités opère pour l'exécution du suivi par lecture des valeurs de pôles, de résidus et de valeur singulière correspondante, ainsi que décrit précédemment dans la description, en liaison avec les figures 2c et 3b à 3g. On comprend, en particulier, que les modules de filtrage 31r 32 et 33 peuvent être constitués par un programme d'ordinateur exécutable par un ordinateur ou par l'unité centrale de traitement d'un dispositif dédié permettant en fait l'application de la méthode du Pencil aux échantillons des composantes d'intensité du courant, de la puissance instantanée et de la tension d'alimentation délivrées à l'installation pour l'obtention des valeurs singulières, pôles et résidus relatifs aux composantes du signal de consommation précitées. On comprend également que les modules de suivi temporel 41r 42, 43 peuvent également être constitués par un programme d'ordinateur exécutable par un ordinateur ou par l'unité centrale de traitement d'un dispositif dédié permettant en fait l'application du processus de suivi temporel de l'enclenchement/déclenchement de charges électriques tel que représenté aux figures 2c, 3b à 3g, ainsi que décrit précédemment dans la description. En ce qui concerne le module 46 de librairie, on indique que ce dernier permet la mise en correspondance biunivoque entre la valeur de signature ou de code de signature V; , la valeur des pôles, résidus et valeur singulière correspondante et en outre de toutes caractéristiques électriques et ou électromécaniques caractéristique de charges électriques normalement disponibles dans le commerce et permettant de déduire un suivi de consommation par charge électrique spécifique. Ce programme d'ordinateur comporte au moins une routine de calcul de la puissance active consommée par chaque charge électrique distincte vérifiant une fonction de puissance consommée liant les pôles et les résidus du signal de consommation, vérifiant la relation (4) précédemment décrite. Dans un exemple de mise en oeuvre non limitatif, la valeur du code allouée à la signature M d'une charge électrique Ci est une valeur unique obtenue par paramétrage d'un processus de signature, tel que le processus ou algorithme RSA (Rivest Shamir Adleman) par exemple, à partir des valeurs de résidus, pôles et valeur singulière correspondante. Le procédé et le dispositif objets de l'invention sont particulièrement adaptés à la mise en oeuvre de processus de détection et de suivi des usages de consommation électrique de clients. Le module 47 d'identification des usages de consommation peut mettre en oeuvre des techniques de détection et de suivi des usages de consommation basées sur des techniques statistiques connues de l'homme de l'art. REVENDICATIONS1. Procédé de détermination non intrusive de la puissance électrique consommée par une installation d'abonné, dans lequel on procède à un échantillonnage périodique des valeurs de tension d'alimentation et de l'intensité du courant électrique délivrés à cette installation pour obtenir un signal de consommation, caractérisé en ce que celui-ci consiste au moins, à partir des valeurs échantillonnées du signal de consommation, à : - grouper lesdites valeurs échantillonnées sur des fenêtres temporelles d'observation successives ; - discriminer sur chaque fenêtre temporelle d'observation courante successives l'ensemble des pôles (si) et des résidus (Ri) du signal de consommation par la méthode du Pencil, ledit ensemble des pôles et résidus étant représentatif, sur la fenêtre d'observation courante, d'une combinaison de signatures de charges électriques (ci) distinctes incluant au moins des pôles (S,) et des résidus (Ri) déterminés associées à au moins une valeur singulière (SGz) obtenue par la méthode du Pencil en fonctionnement dans ladite installation ; - calculer, sur au moins la durée de ladite fenêtre d'observation courante, au moins la puissance active consommée par au moins un sous-ensemble des charges électriques enfonctionnement dans l'installation, exprimée comme la somme des puissances actives consommées par chaque charge électrique distincte de ce sous-ensemble vérifiant une fonction de puissance consommée liant lesdits pôles et résidus. 2. Procédé selon la revendication 1, caractérisé en ce que, pour ladite puissance active, ladite fonction 10 vérifie la relation : M M' (s .+svi, P = LERci.Rvi,.e i=1 i'=l+~t l At sinh c (Scl + Svf ) 2 où P : désigne la puissance active consommée par 15 l'ensemble des charges électriques distinctes CZ de signature (SZ,RZ,SGZ) ; S,etR, : désignent les pôles respectivement les résidus du courant échantillonné délivré à l'installation ; Svl, et Rvl, : désignent les pôles respectivement les résidus 20 de la tension échantillonnée délivrée à l'installation, SGi désignant les valeurs singulières ; tkd : désigne l'instant de début de la fenêtre d'observation de rang k; At : désigne la durée de la fenêtre d'observation ; 25 sinhc(X) : la fonction sinus hyperbolique du cardinal de la valeur X, soit sinh c((Sci + Svi) Ot =sinh ((Sci + Se ) & / ~Sei +M : désigne le nombre de pôles et résidus du courant échantillonné délivré à l'installation ; M' : désigne le nombre de pôles et résidus de la tension échantillonné délivrée à l'installation. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que pour une tension d'alimentation sensiblement périodique de pulsation w et de période T, la durée At de chaque fenêtre d'observation courante est prise égale à ladite période. 4. Procédé selon l'une des revendications 1 à 3, caractérisé en ce que la fréquence d'échantillonnage pour obtenir ledit signal de consommation et ladite durée At de chaque fenêtre d'observation sont liées par la relation : Nc At.fe=Nem2. Lnci+2(1+H) dans laquelle fe : désigne la fréquence d'échantillonnage ; Nem : désigne le nombre minimum d'échantillons sur la durée de chaque fenêtre d'observation ; Nc : désigne le nombre observé de charges électriques distinctes (ci) de l'installation d'abonné ; nci : désigne le nombre de pôles de la charge électrique K H : désigne le nombre de composantes harmoniques de la fréquence fondamentale de la tension d'alimentation.30 5. Procédé selon l'une des revendications 1 à 4, caractérisé en ce que, suite à l'étape consistant à discriminer l'ensemble des pôles (S,) et des résidus (Ri) du signal de consommation, celui-ci comporte en outre une étape de discrimination et d'acquisition des signatures formées par un couple de pôles (Sz), résidus (Ri) de chaque charge électrique distincte (ci) et des valeurs singulières (SG,). 6. Procédé selon la revendication 5, caractérisé en ce que, suite à l'acquisition desdites signatures (SZ,RZ,SGZ) ledit procédé comprend une étape d'identification desdites charges électriques constitutives dudit sous-ensemble de charges électriques distinctes en fonctionnement sur au moins une fenêtre d'observation. 7. Procédé selon la revendication 5 ou 6, caractérisé en ce que, suite à l'acquisition desdites signatures (Si,Ri,SGi), ledit procédé comprend une étape de suivi du fonctionnement/mise hors fonctionnement de chaque charge distincte, par discrimination de tout ou partie de la signature sur au moins une fenêtre d'observation. 8. Procédé selon l'une des revendications 6 ou 7, caractérisé en ce que, pour une fenêtre d'observation du signal de consommation incluant le régime libre et le commencement du régime forcé d'enclenchement/déclenchement d'une charge, - l'enclenchement/déclenchement d'une charge purement résistive est signé par un régime forcé unique, par l'existence, pour la fenêtred'observation, de deux pôles complexes conjugués à partie réelle nulle et de deux résidus complexes conjugués dans le régime forcé, et par la phase nulle des résidus en régime forcé ; - l'enclenchement/déclenchement d'une charge résistive et capacitive est signé par un régime libre et un régime forcé distincts, l'existence, pour la fenêtre d'observation, de deux pôles complexes conjugués à partie réelle nulle et de deux résidus complexes conjugués dans le régime forcé et par un pôle réel et un résidu dans le régime libre, la phase négative du résidu en régime forcé caractérisant la nature globalement capacitive de la charge ; - l'enclenchement/déclenchement d'une charge résistive et inductive est signé par un régime libre et un regime forcé distincts, par l'existence, pour la fenêtre d'observation, deux pôles, complexes conjugués à partie réelle nulle et de deux résidus complexes conjugués dans le régime forcé et par un pôle réel et un résidu complexe, dans le régime libre, la phase positive du résidu caractérisant la nature globalement inductive de la charge. 25 9. Procédé selon l'une des revendications 6 ou 7, caractérisé en ce que, pour une fenêtre d'observation du signal de consommation incluant le régime libre et le début du régime forcé d'enclenchement/déclenchement d'une 30 charge, l'enclenchement/déclenchement d'une charge résistive, inductive et capacitive est signé par un 20régime forcé et un régime libre distincts, par l'existence : en régime forcé, de deux pôles complexes conjugués à partie réelle nulle et de deux résidus complexes conjugués ; et, en régime libre, de deux pôles complexes conjugués et deux résidus complexes conjugués, la phase négative du résidu en régime forcé caractérisant la nature globalement capacitive de la charge et la phase positive du résidu en régime forcé caractérisant la nature globalement inductive de la charge. 10. Dispositif de détermination non intrusive de la puissance électrique consommée par une installation d'abonné, comportant au moins des moyens d'échantillonnage périodique des valeurs de tension d'alimentation et de l'intensité du courant électrique délivrés à cette installation pour obtenir un signal de consommation, caractérisé en ce que ledit dispositif comporte en outre au moins . des moyens de calcul de la puissance instantanée délivrée produit des valeurs échantillonnées des valeurs de tension d'alimentation et de l'intensité du courant électrique délivrés, ledit signal de consommation comportant l'une au moins des valeurs échantillonnées de la tension d'alimentation, de l'intensité du courant électrique délivrés ou de la puissance instantanée délivrée ;- des moyens de mémorisation et de lecture des valeurs successives du signal de consommation sur au moins une fenêtre d'observation ; - des moyens d'application sur le signal de consommation, sur au moins une fenêtre d'observation courante, d'un traitement de filtrage par la méthode du Pencil, pour discriminer, sur ladite fenêtre, l'ensemble des pôles et résidus représentatifs d'une combinaison de signatures de charges électriques distinctes incluant au moins les pôles et les résidus associés à au moins une valeur singulière obtenue par la méthode du Pencil en fonctionnement dans l'installation ; - des moyens de suivi temporel de ladite combinaison de signatures et de calcul, sur au moins la durée de la fenêtre d'observation courante, d'au moins la puissance active consommée par au moins un sous-ensemble des charges électriques en fonctionnement dans l'installation, exprimée comme la somme des puissances actives consommées par chaque charge électrique distincte de ce sous-ensemble vérifiant une fonction de puissance consommée 25 liant lesdits pôles et résidus. 11. Dispositif selon la revendication 10, caractérisé en ce que lesdits moyens de suivi temporel de ladite combinaison de signatures et de calcul d'au moins 30 la puissance active consommée par au moins un sous-ensemble des charges électriques comprennent : 15 20- des modules séparés opérant sur les pôles, résidus, valeurs singulières de la composante d'intensité, de puissance respectivement de tension délivrée à l'installation ; - un module de détection de changement d'état d'enclenchement/déclenchement de charge électrique ; un module d'identification des charges électriques enclenchées et/ou déclenchées à partir des signatures de charges électriques ; - un module de librairie de charges électriques identifiées en relation avec leur signature ; un module d'identification des usages de consommation de l'utilisation abonné et de calcul de la consommation, au moins en puissance active. 12. Programme d'ordinateur mémorisé sur un support de mémorisation et exécutable par un ordinateur ou par l'unité centrale de traitement d'un dispositif dédié, ce programme comportant une suite d'instructions permettant l'application de la méthode du Pencil aux échantillons des composantes d'intensité du courant, de la puissance instantanée et de la tension d'alimentation délivrées à une installation. 13. Programme d'ordinateur mémorisé sur un support de mémorisation et exécutable par un ordinateur ou par l'unité centrale de traitement d'un dispositif dédié, ce programme d'ordinateur comportant une suite d'instructions permettant l'application d'un processus de suivi temporel de l'enclenchement/déclenchement de charges électriques selon l'une des revendications 1 à 7. 14. Programme d'ordinateur selon la revendication 13, caractérisé en ce que ledit programme d'ordinateur comporte au moins une routine de calcul de la puissance 5 active consommée par chaque charge électrique distincte vérifiant une fonction de puissance consommée liant les pôles et les résidus du signal de consommation.