FR 3074986 A1 20190614 FR 1761756 A 20171207 DETECTION A DISTANCE D'UNE ANOMALIE AFFECTANT UNE INSTALLATION PHOTOVOLTAÏQUE La présente invention concerne la détection à distance d'une anomalie qui affecte le fonctionnement d'une installation photovoltaïque de production d'électricité. La production d'électricité à partir d'une installation photovoltaïque peut être perturbée par de nombreuses causes, y compris des accumulations de salissures et de détritus sur les panneaux photovoltaïques, des pannes de l'installation pouvant concerner les panneaux eux-mêmes mais aussi des composants électriques qui leur sont associés, l'apparition de nouveaux éléments dans l'environnement qui provoquent un ombrage des panneaux, etc. Il en résulte des baisses dans la production d'électricité de l'installation photovoltaïque, qui peuvent être particulièrement pénalisantes lors de journées ensoleillées mais froides. De telles baisses sont pénalisantes à la fois pour le propriétaire de l'installation photovoltaïque à cause du défaut de puissance électrique à vendre, mais aussi pour un fournisseur d'électricité par réseau de distribution à cause d'une puissance électrique manquante qui pourrait devoir être achetée à d'autres fournisseurs. Or l'apparition de telles causes de perturbations peut être difficile à détecter, causant des pertes de production qui sont susceptibles de durer plusieurs jours, voire plusieurs semaines. En outre, une fois que l'existence d'une perturbation a été détectée pour une installation photovoltaïque, identifier la ou les cause(s) de cette perturbation peut être difficile. Or il est particulièrement important, pour des raisons économiques, de diminuer au maximum la nécessité de missionner un technicien sur le site de l'installation photovoltaïque pour identifier la cause de perturbation. En effet, les installations photovoltaïques sont en général très réparties géographiquement, si bien que le déplacement d'un technicien sur une installation particulière possède un coût élevé. Il y a donc un enjeu économique important à pouvoir détecter à distance une perturbation qui - 2 - affecte nouvellement une installation photovoltaïque identifiée, et à pouvoir diagnostiquer la cause de cette perturbation. Pour cela, il est connu de transmettre à une plateforme de surveillance, des mesures qui sont effectuées répétitivement, de la production d'électricité 5 de chaque installation photovoltaïque, afin d'étudier leur évolution et tenter d'en déduire la survenue de perturbations. Usuellement, ces valeurs mesurées sont comparées à des valeurs théoriques qui sont produites par un simulateur selon un modèle déterministe. Mais une telle approche est inadaptée, car la production d'électricité par une installation photovoltaïque est soumise à de 10 nombreux aléas, notamment d'origine météorologique, et dépend de paramètres de l'installation qui ne sont pas forcément connus pour être pris en compte dans le modèle de simulation. Parmi ces paramètres d'installation figurent le lieu de chaque installation photovoltaïque, la surface de panneaux installée, l'orientation des panneaux photovoltaïques par rapport au sud, leur 15 inclinaison par rapport à la direction verticale, la technologie et le fabricant des panneaux, etc. La détection de perturbations à partir de comparaisons entre les mesures de production d'électricité qui sont collectées et les résultats d'un tel modèle déterministe aboutit alors à un taux de fausse alarme qui est très élevé. Autrement dit, de telles méthodes de détection de perturbations à partir de 20 modèles déterministes ne sont pas fiables. A partir de cette situation, un but de l'invention est de proposer un nouveau procédé de détection à distance d'anomalies qui affectent le fonctionnement d'une installation photovoltaïque, dont le niveau de fiabilité est amélioré par rapport aux procédés de détection antérieurs. 25 Un but annexe de l'invention est d'identifier à distance au moins une cause probable pour chaque anomalie qui est détectée dans le fonctionnement d'une installation photovoltaïque. Pour atteindre l'un au moins de ces buts ou d'autres, un premier aspect de l'invention propose un procédé de détection à distance d'une anomalie qui affecte un fonctionnement d'une installation photovoltaïque de production d'électricité, dite installation surveillée, ce procédé comprenant au moins une étape lors de laquelle des valeurs de production d'électricité qui sont mesurées - 3 - pour cette installation surveillée sont comparées avec des valeurs qui sont déduites de mesures effectuées pour plusieurs autres installations photovoltaïques de production d'électricité, dites installations voisines et différentes de l'installation surveillée. 5 Ainsi, le procédé de détection de l'invention peut prendre en compte des contributions à des variations aléatoires de production d'électricité, qui sont communes à des installations photovoltaïques voisines géographiquement. Une partie des causes de variation de la production d'électricité qui affectent l'installation surveillée peut donc être prise en compte de cette façon, 10 notamment une partie des causes aléatoires, si bien que le procédé de détection de l'invention possède un niveau de fiabilité amélioré. En outre, ce niveau de fiabilité est d'autant plus élevé qu'un nombre plus grand d'installations photovoltaïques voisines sont prises en compte. De façon générale, le procédé de l'invention peut comprendre en outre 15 de rechercher si plusieurs critères de détection d'anomalie sont satisfaits par les valeurs de production d'électricité qui sont mesurées pour l'installation surveillée, et en fonction de ceux de ces critères qui sont satisfaits, proposer au moins une cause possible pour l'anomalie. En particulier, la proposition de cause possible d'anomalie peut utiliser une base d'apprentissage qui associe 20 des causes d'anomalies différentes à des combinaisons différentes des critères satisfaits. Eventuellement, le procédé de l'invention peut encore comprendre de commander ou réaliser une intervention de maintenance ou de modification environnementale pour l'installation surveillée, sur un site ou dans un voisinage 25 de cette installation surveillée, si une anomalie est détectée. De façon générale, le procédé de l'invention peut comprendre les étapes suivantes : /1/ pendant une période d'apprentissage, collecter les valeurs de production d'électricité qui ont été mesurées pour l'installation 30 surveillée et pour les installations voisines ; /2/ à partir des valeurs de production d'électricité qui ont été mesurées pour les installations voisines, déterminer une fonction d'estimation de - 4 - la production d'électricité pour l'installation surveillée, ayant le temps comme variable, en minimisant un écart total qui est calculé sur la période d'apprentissage entre les valeurs de production d'électricité mesurées pour l'installation surveillée et des valeurs estimées 5 obtenues par la fonction d'estimation ; /3/ pour au moins un instant postérieur à la période d'apprentissage, collecter de nouvelles valeurs de production d'électricité qui sont mesurées pour l'installation surveillée et pour chaque installation voisine ; 10 /4/ pour chaque instant postérieur à la période d'apprentissage qui est concerné par l'étape /3/, calculer au moins une valeur estimée de production d'électricité en utilisant la fonction d'estimation qui a été déterminée à l'étape /2/ ; et /5/ déclarer qu'une anomalie est détectée dans le fonctionnement de 15 l'installation surveillée à partir d'une comparaison, pour au moins un instant postérieur à la période d'apprentissage, entre la valeur estimée de production d'électricité qui a été calculée à l'étape /4/ et la nouvelle valeur de production d'électricité qui a été mesurée puis collectée à l'étape /3/ pour l'installation surveillée. 20 De cette façon, le procédé de l'invention résulte d'une comparaison entre la production d'électricité de l'installation surveillée et un modèle pour cette production d'électricité qui est établi à partir des installations voisines de l'installation surveillée. Et une anomalie est déclarée pour l'installation surveillée lorsque la production d'électricité qui est mesurée pour celle-ci 25 devient trop différente d'un résultat fourni par le modèle. Eventuellement, les nouvelles valeurs de production d'électricité mesurées qui sont collectées à l'étape /3/ peuvent en outre être utilisées ultérieurement pour rallonger la période d'apprentissage. Elles peuvent alternativement être utilisées aussi pour former une période d'apprentissage 30 glissante dans le temps, qui incorpore alors automatiquement de nouvelles dérives de production susceptibles d'affecter l'installation photovoltaïque surveillée et/ou ses voisines, dans la mesure où de telles dérives sont -5 - progressives. Dans des modes préférés de mise en oeuvre de l'invention, la fonction d'estimation qui est déterminée à l'étape /2/ peut comprendre une combinaison linéaire de fonctions de production d'électricité qui sont relatives 5 respectivement aux installations voisines, et qui correspondent aux valeurs de production d'électricité qui ont été mesurées pour ces installations voisines pendant la période d'apprentissage, avec des facteurs de combinaison linéaire qui sont tous positifs ou nuls, et avec une somme desdits facteurs de combinaison linéaire qui est inférieure ou égale à une valeur maximale fixée. 10 En outre, les facteurs de combinaison linéaire peuvent être déterminés pour minimiser l'écart total qui est calculé sur la période d'apprentissage entre les valeurs de production d'électricité mesurées pour l'installation surveillée et les valeurs estimées. Pour un premier critère de détection d'anomalie, la comparaison de 15 l'étape /5/ peut être effectuée en exécutant les sous-étapes suivantes : /5-1/ déterminer un intervalle de confiance de prédiction pour chaque valeur estimée de production d'électricité qui a été calculée à l'étape /4/ ; /5-2/ déterminer si la nouvelle valeur de production d'électricité mesurée 20 qui a été collectée à l'étape /3/ pour l'installation surveillée, et qui correspond à un même instant que la valeur estimée de la sous-étape /5-1/, est à l'extérieur de l'intervalle de confiance de prédiction relatif à cette valeur estimée ; et /5-3/ répéter les sous-étapes /5-1/ et /5-2/ pour plusieurs instants à 25 l'intérieur d'une première période de test, puis compter un premier nombre d'instants dans la première période de test, pour chacun desquels instants la nouvelle valeur de production d'électricité mesurée pour l'installation surveillée est à l'extérieur de l'intervalle de confiance de prédiction correspondant. 30 Alors, selon le premier critère, une anomalie dans le fonctionnement de l'installation surveillée peut être déclarée si le premier nombre d'instants est supérieur à un premier seuil prédéfini. - 6 - Possiblement, le procédé de l'invention peut aussi comprendre les étapes suivantes, pour un deuxième critère de détection d'anomalie : /6/ corriger les valeurs de production d'électricité qui sont mesurées pour chaque installation voisine par rapport à une puissance électrique 5 nominale de cette installation voisine, et proportionnellement à une puissance électrique nominale de l'installation surveillée ; /7/ sélectionner un sous-ensemble parmi les installations voisines, qui correspond à une proportion prédéterminée de celles-ci et qui est formé par celles de ces installations voisines qui possèdent des 10 productions d'électricité les plus basses, déterminées à partir des valeurs de production d'électricité mesurées puis corrigées ; et /8/ compter un second nombre d'instants dans une deuxième période de test, pour chacun desquels instants la valeur de production d'électricité qui est mesurée pour l'installation surveillée est inférieure ou égale à 15 une valeur moyenne des valeurs de production d'électricité mesurées puis corrigées du sous-ensemble d'installations voisines. Alors, selon le deuxième critère, une anomalie dans le fonctionnement de l'installation surveillée est aussi déclarée si le second nombre d'instants qui a été compté à l'étape /8/ est supérieur à un deuxième seuil prédéfini. La 20 deuxième période de test peut être une partie au moins de la période d'apprentissage, et/ou comprendre au moins un intervalle de temps qui est postérieur à la période d'apprentissage. Pour un troisième critère de détection d'anomalie, le procédé peut comprendre les étapes /1/ à /4/ mentionnées plus haut, et comprendre en plus 25 les étapes suivantes : /9/ pour plusieurs instants à l'intérieur d'une troisième période de test, calculer des résultats de soustraction d'une valeur de production d'électricité qui est calculée en utilisant la fonction d'estimation déterminée à l'étape /2/, à la valeur de production d'électricité mesurée 30 pour l'installation surveillée, ces valeurs pour calculer chaque résultat de soustraction étant relatives à un même instant dans la troisième période de test ; - 7 - /10/ déterminer des durées de marche, pendant chacune desquelles durées de marche les résultats de soustraction présentent une valeur de moyenne temporelle qui est différente de celles de la durée de marche précédente et de la durée de marche suivante, et avec une 5 variance des résultats de soustraction qui est minimale à l'intérieur de chaque durée de marche ; et /11/ compter un nombre de durées de marche à l'intérieur de la troisième période de test, pour chacune desquelles durées de marche la valeur de moyenne temporelle des résultats de soustraction est négative. 10 Alors, selon le troisième critère, une anomalie dans le fonctionnement de l'installation surveillée est aussi déclarée si, pendant la troisième période de test, le nombre de durées de marche à valeurs de moyenne temporelle négatives pour les résultats de soustraction est supérieur à un troisième seuil prédéfini. La troisième période de test peut être une partie au moins de la 15 période d'apprentissage, et/ou comprendre au moins un intervalle de temps qui est postérieur à la période d'apprentissage. Un tel troisième critère est particulièrement adapté pour révéler des baisses de production d'électricité qui affectent l'installation photovoltaïque qui est surveillée. Enfin, pour un quatrième critère de détection d'anomalie, le procédé 20 peut comprendre en outre les étapes suivantes : /12/ rechercher s'il existe un instant diurne auquel un minimum se produit pour les valeurs de production d'électricité qui sont mesurées pour l'installation surveillée, pour au moins un jour de fonctionnement de cette installation surveillée ; et 25 /13/ rechercher si un minimum se produit aussi pendant plusieurs jours successifs après le jour identifié à l'étape /12/, pour les valeurs de production d'électricité qui sont mesurées pour l'installation surveillée, chaque fois dans un intervalle de temps prédéfini autour de l'instant diurne identifié à l'étape /12/, et compter le nombre de jours successifs 30 pour lesquels un tel minimum se produit. Alors selon le quatrième critère, une anomalie dans le fonctionnement de l'installation surveillée est aussi déclarée si le nombre de jours successifs - 8 - qui a été compté à l'étape /13/ est supérieur à un quatrième seuil prédéfini. De façon générale, la vérification de plusieurs des critères précédents peut permettre de confirmer la présence d'une anomalie dans le fonctionnement de l'installation surveillée, voire permettre de déterminer des 5 causes possibles pour cette anomalie détectée. Un second aspect de l'invention propose un système de surveillance d'au moins une installation photovoltaïque, dite installation surveillée, ce système de surveillance comprenant : - des moyens de communication pour collecter des valeurs de production 10 d'électricité qui ont été mesurées pour l'installation surveillée et pour d'autres installations photovoltaïques, dites installations voisines, différentes de l'installation surveillée ; et - des moyens de calcul pour mettre en oeuvre un procédé de détection à distance d'une anomalie qui affecte un fonctionnement de l'installation 15 surveillée, ce procédé étant conforme au premier aspect de l'invention. D'autres particularités et avantages de la présente invention apparaîtront dans la description ci-après d'exemples de réalisation non limitatifs, en référence aux dessins annexés, dans lesquels : - la figure 1 montre un ensemble d'installations photovoltaïques auquel 20 l'invention peut être appliquée, ainsi qu'un système de surveillance conforme à l'invention ; et - les figures 2 à 5 montrent des étapes d'un procédé conforme à l'invention, respectivement pour appliquer des critères différents de détection d'anomalie concernant le fonctionnement d'une des 25 installations photovoltaïques de la figure 1. Dans ces figures, des références identiques désignent des éléments ou des étapes identiques, ou dont les fonctions sont identiques. La figure 1 montre N installations photovoltaïques identifiées, N étant un nombre entier naturel supérieur ou égal à 2, par exemple référencées de 1 30 à 10, dont les conditions d'installation peuvent varier d'une installation à l'autre : en toitures, sur pieds fixes, avec des surfaces installées de panneaux - 9 - qui sont variables, des orientations variables par rapport au sud, des inclinaisons qui peuvent aussi être variables, etc. En outre, la technologie et le fabricant des panneaux photovoltaïques peuvent aussi varier entre des installations différentes. Chaque installation photovoltaïque est connectée à un 5 réseau de distribution électrique (non représenté) pour transmettre la puissance électrique qu'elle fournit, et un compteur intermédiaire effectue des mesures de cette puissance électrique fournie. Des valeurs mesurées de production d'électricité sont ainsi générées pour chaque installation photovoltaïque, individuellement et de façon répétée, puis transmises à une 10 plateforme de réception et de mémorisation 100, qui est notée RECEP. & MEM. Le mode de transmission de ces valeurs mesurées, à partir d'un point de connexion au réseau de chaque installation photovoltaïque, peut être quelconque. Pour mettre en oeuvre le procédé de détection d'anomalie qui est proposé par l'invention, un module d'analyse 101, noté ANALYS., est associé à 15 la plateforme 100, et programmé pour appliquer aux valeurs de production d'électricité mesurées et transmises, les traitements qui sont décrits ci-après. Ces traitements ont pour but de détecter si une anomalie affecte le fonctionnement de l'une des installations photovoltaïques identifiées, par exemple l'installation 1, qui est donc appelée installation photovoltaïque 20 surveillée. Pour cela, les valeurs de production d'électricité qui sont mesurées et collectées pour l'installation photovoltaïque 1 sont comparées à celles relatives à plusieurs autres des installations photovoltaïques. Le nombre de ces dernières installations utilisées comme références de comparaison n'a pas de limitation. Dans la pratique, il peut être de l'ordre d'une dizaine d'installations 25 photovoltaïques ou plus, de préférence inférieur à cent, par exemple les installations photovoltaïques 2-10. La fiabilité du procédé de détection d'anomalie pour l'installation 1 qui est surveillée, est meilleure si les installations avec lesquelles ses valeurs mesurées de production d'électricité sont comparées, sont voisines géographiquement de l'installation 1. Pour cette 30 raison, les installations photovoltaïques qui servent de références pour surveiller le fonctionnement de l'installation 1, ont été appelées installations voisines dans la partie générale de la présente description. L'Homme du métier comprendra que pour l'analyse statistique qui est mise en oeuvre, le critère de - 10 - voisinage n'a pas de limite stricte, mais a pour effet d'améliorer des corrélations temporelles entre les valeurs qui sont mesurées pour des installations photovoltaïques distinctes. Dans les figures 2 à 5, la collecte des valeurs de production 5 d'électricité qui ont été mesurées pour chacune des installations photovoltaïques 1 à 10, pendant une période d'apprentissage, est notée SO. Une telle collecte est répétée plusieurs fois dans la période d'apprentissage, préférablement pendant des périodes diurnes de plusieurs jours successifs de la période d'apprentissage, pour construire une base de données de valeurs 10 mesurées. Par exemple, la période d'apprentissage peut regrouper entre 5 et 50 jours successifs. Cette base de données peut être mise à jour continuellement. Dans la suite, t désigne la variable de temps, pour repérer des instants auxquelles les valeurs de production d'électricité sont mesurées pour les installations photovoltaïques 1 à N. 15 Pour appliquer le premier critère de détection d'anomalie prévu par l'invention, les valeurs mesurées qui sont ainsi collectées pendant la période d'apprentissage (étape SO) sont utilisées pour établir un modèle de la production d'électricité de l'installation surveillée (étape S11 de la figure 2). Ensuite, des prédictions qui sont établies en appliquant ce modèle (étape S12) 20 sont comparées à de nouvelles valeurs de production d'électricité de l'installation surveillée, qui sont mesurées postérieurement à la période d'apprentissage (étapes S13 et S14). De préférence, le modèle qui est utilisé peut être du type connu sous l'acronyme LASSO, pour «Least Absolute Shrinkage and Selection Operator», à coefficients positifs et bornés. Selon un 25 tel modèle, les valeurs de production d'électricité de l'installation surveillée sont approchées par une combinaison affine des valeurs de production d'électricité des installations voisines : i-N Pi(t) =0 +1A 11(t)+E1(t) i=2 où Pi(t) est la valeur de production d'électricité qui est mesurée à l'instant t pour 30 l'installation photovoltaïque i, l'indice entier i variant de 1 à N, i=1 correspondant à l'installation photovoltaïque qui est surveillée, i=2 à N correspondant aux installations photovoltaïques voisines de celle qui est surveillée, et la sommation formant une combinaison linéaire des valeurs de production d'électricité qui sont mesurées pour les installations photovoltaïques voisines ; 5 131,i, pour i variant de 2 à N, étant le coefficient de la combinaison linéaire pour l'installation voisine i dans le modèle qui est utilisé pour estimer la production d'électricité de l'installation surveillée 1, et 131,0 étant un coefficient constant de ce modèle pour l'installation surveillée 1 ; et ci(t) étant l'erreur du modèle pour l'installation surveillée 1, pour l'instant de 10 mesure t. Les coefficients 131,i, avec i=0, N, sont déterminés pour minimiser l'erreur totale, intégrée sur la période d'apprentissage : i=N 131,0 ±I 131,i Pi(t) i=2 où les doubles barres verticales désignent une norme de mesure qui est 15 effective sur toute la période d'apprentissage. Par exemple : i=N Pico-fli 0 -1 2 11(0 i=2 tE périoded apprentissage La fonction ci(t) est alors une fonction aléatoire de valeur moyenne nulle. Selon une caractéristique de l'invention, les coefficients 131,; pour i=0, 20 N sont tous positifs ou nuls, avec la contrainte que leur somme soit inférieure ou égale à un paramètre fixé du modèle, noté Ô, et qui un nombre réel strictement positif : 131,0 pour i = 0, 2, ..., N et i=N 131,0 131,i i=2 ) Le but de la période d'apprentissage introduite ci-dessus est de 25 déterminer les coefficients 13-1,i qui sont ainsi définis, pour i=0, N, correspondant à l'étape S11 indiquée dans la figure 2. Ces coefficients sont ensuite utilisés à l'étape S12 pour estimer la production d'électricité de - 12 - l'installation surveillée 1, pour des instants t qui sont postérieurs à la période d'apprentissage. L'estimation de la production d'électricité par l'installation surveillée 1 pour l'instant t est : i-N 5 (t)= )6 1 0 Pt(t) i=2 Ainsi, le modèle, aussi appelé fonction d'estimation, fournit une valeur estimée pour la production d'électricité de l'installation surveillée 1, pour un instant t, à partir des valeurs de production d'électricité des installations voisines 2,..., N qui sont relatives au même instant t. 10 Pour appliquer le premier critère de détection d'anomalie, un intervalle de confiance de prédiction est déterminé, séparément pour plusieurs instants t qui sont postérieurs à la période d'apprentissage, pour l'estimation de la production d'électricité de l'installation surveillée 1 qui est obtenue par le modèle précédent. Un tel intervalle de confiance peut être déterminé de la 15 façon suivante, à partir des valeurs mesurées qui ont été collectées pendant la période d'apprentissage : /i/ on sélectionne de façon aléatoire, avec remise, des instants t de la période d'apprentissage, pour lesquels les valeurs mesurées Pi(t) ont été collectées, pour i=1, N, en répétant le tirage un nombre de 20 fois qui est identique au nombre d'instants de la période d'apprentissage qui ont été utilisés précédemment pour déterminer la fonction d'estimation Pi(t). On construit ainsi une nouvelle base de données d'apprentissage, à partir de laquelle une nouvelle fonction d'estimation est déterminée, de la même façon que décrite ci-dessus, 25 de nouveau pour la production d'électricité de l'installation surveillée 1 ; /ii/ on sélectionne ensuite, aussi de façon aléatoire, une valeur de l'erreur d'estimation qui est relative à la nouvelle fonction d'estimation ainsi construite, et relative à l'un des instants de la période d'apprentissage ; puis 30 /Hi/ on calcule une nouvelle valeur d'estimation de la production - 13 - d'électricité pour l'installation sélectionnée 1, comme la somme d'un résultat de la nouvelle fonction d'estimation déterminée à l'étape /i/, pour un instant postérieur à la période d'apprentissage, et de l'erreur d'estimation qui a été sélectionnée à l'étape /ii/. 5 Ces étapes /i/ à /iii/ sont répétées un nombre déterminé de fois, par exemple 500 fois, la première fois étant exécutée à partir de la fonction d'estimation îl(t), pour obtenir une distribution statistique de valeurs d'estimation de la production d'électricité de l'installation sélectionnée 1 pour un même instant t, postérieur à la période d'apprentissage. Un intervalle de 10 confiance est alors déterminé à partir de cette distribution statistique pour la production d'électricité de l'installation 1, pour l'instant t. Par exemple, l'intervalle de confiance à 95% est compris entre quantiles 0,025 et 0,975. /iv/ enfin, on détermine si la valeur qui a été réellement mesurée puis collectée à l'étape S13 pour la production d'électricité de l'installation 1 15 à l'instant t postérieur à la période d'apprentissage, est à l'intérieur de l'intervalle de confiance, ou à l'extérieur de celui-ci. Cette comparaison de la valeur mesurée avec les valeurs estimées, par l'intermédiaire de l'intervalle de confiance, constitue l'étape S14. Les étapes /i/ à /iv/ sont ensuite appliquées à plusieurs instants t 20 différents, tous postérieurs à la période d'apprentissage, puis on compte le nombre N1 de ces instants, à l'intérieur d'une première période de test prédéfinie, pour lesquels les valeurs qui sont réellement mesurées pour la production d'électricité de l'installation 1 sont à l'extérieur des intervalles de confiance correspondant aux mêmes instants. Si ce nombre N1 d'instants, pour 25 lesquels surviennent des défauts de coïncidence entre le modèle et les valeurs mesurées pour l'installation 1, est supérieur à un premier seuil prédéfini, une anomalie est déclarée pour le fonctionnement de l'installation 1. Ainsi, le premier critère de détection d'anomalie est constitué par le niveau de confiance de prédiction (95% dans l'exemple cité), la première période de test et le 30 premier seuil qui ont été utilisés dans le mode d'analyse qui vient d'être décrit, et qui a été appliqué aux valeurs mesurées de production d'électricité pour les installations 1, 2,...,N. Ce premier critère est mis en oeuvre à l'étape S10. Il - 14 - permet de prendre en compte des particularités de fonctionnement de l'installation qui est surveillée, puisque chaque fonction d'estimation qui a été utilisée a été déterminée à partir de ce fonctionnement de l'installation surveillée, tel que mesuré pendant la période d'apprentissage. 5 Un deuxième critère de détection d'anomalie pour le fonctionnement de l'installation surveillée concerne la comparaison de sa production d'électricité à celles d'une partie des installations voisines qui sont les moins efficaces. Pour les installations voisines pour lesquelles le coefficient de combinaison linéaire 13-1,i n'est pas nul dans la fonction d'estimation 1(t), 10 l'indice i variant de 2 à N, on corrige leurs valeurs mesurées de production d'électricité par rapport à leurs valeurs respectives de puissance électrique nominale, et proportionnellement à celle de l'installation surveillée. Autrement dit, pour chaque instant t de la période d'apprentissage auquel la valeur de production d'électricité Pi(t) a été mesurée puis collectée, on calcule la valeur 15 corrigée de production d'électricité suivante, à partir de la valeur mesurée Pi(t) : pour i = 1,2.....N : COÏT (t) = p(t) - Pi nom I Pi nom OÙ Pi non, est la valeur de puissance électrique nominale de l'installation photovoltaïque i, pour i variant de 1 à N. Ces calculs de valeurs corrigées correspondent à l'étape S21 de la figure 3. 20 Pour chaque installation voisine i, i variant de 2 à N, on détermine ensuite la valeur moyenne de sa production d'électricité corrigée sur la période d'apprentissage, à partir des valeurs de Pi corr(t) : t. Puis on sélectionne, à l'étape S22, une proportion ou un nombre prédéfini des installations voisines, correspondant à celles qui possèdent de telles valeurs 25 moyennes les plus faibles. Par exemple, on sélectionne les 20% des installations voisines qui possèdent les valeurs les plus faibles pour la production d'électricité corrigée moyenne. Enfin, à l'étape S23, on détermine le nombre N2 d'instants pour lesquels les valeurs mesurées de production d'électricité ont été collectées, 30 pendant la période d'apprentissage, et pour lesquels la valeur mesurée pour l'installation surveillée 1 est inférieure ou égale à la valeur moyenne des - 15 - valeurs mesurées au même instant pour les installations voisines sélectionnées : t dans la période d'apprentissage tel que Pi(t) 1 N sélection ie sélection parmi 2 N OÙ Nsélection est le nombre des installations photovoltaïques voisines les moins 5 productives qui ont été sélectionnées. Alors le deuxième critère de détection d'anomalie pour le fonctionnement de l'installation surveillée 1 est un deuxième seuil fixé, qui concerne le nombre d'instants N2. Une anomalie est alors déclarée selon ce deuxième critère si le nombre d'instants N2 est supérieur au deuxième seuil. 10 Un troisième critère de détection d'anomalie pour le fonctionnement de l'installation surveillée est fondé sur l'existence de ruptures de statistique dans les différences Pi(t)-1(t), ces différentes étant calculées pour les instants de mesure t de la période d'apprentissage. Les étapes S11 et S12 de la figure 4 sont une reprise de celles de la figure 1. 15 Puis à l'étape S31, une méthode détection de ruptures de statistique est appliquée aux différences Pi(t)-1(t). Pour cela, on peut appliquer un algorithme de programmation dynamique sur la période d'apprentissage. Un tel algorithme découpe la période d'apprentissage en durées de marche («step durations» en anglais) successives, chacune étant associée à une valeur de 20 moyenne temporelle pour la différence Pi(t)-1(t), calculée sur cette durée de marche, qui est différente de la valeur de moyenne temporelle associée à la durée de marche suivante, et aussi différente de celle associée à la durée de marche précédente. En outre, chaque durée de marche est déterminée de sorte que la variance de la différence F(t)-t(t) soit minimale à l'intérieur de 25 cette durée de marche. De telles durées de marche peuvent être déterminées en utilisant par exemple une méthode de nuées dynamiques, ou une méthode de classification ascendante hiérarchique, connues de l'Homme du métier. Les durées de marche dont les valeurs de moyenne temporelle sont négatives pour la différence Pi(t)-](t) sont comptées, et leur nombre est comparé à un 30 troisième seuil prédéfini à l'étape S30. Une anomalie est alors déclarée selon - 16 - le troisième critère si ce nombre est supérieur au troisième seuil. La valeur de moyenne temporelle la plus négative qui a été obtenue ainsi pour la différence Pi(t)-1(t), et les instants de mesure t de la durée de marche correspondante, ou bien une heure moyenne de cette durée de 5 marche à valeur de moyenne temporelle la plus négative, peuvent être utilisés pour identifier une cause probable de l'anomalie, comme cela sera décrit plus loin. Enfin, un quatrième critère de détection d'anomalie pour le fonctionnement de l'installation surveillée consiste à rechercher si un minimum 10 apparaît pour sa production d'électricité telle que mesurée, de façon répétitive pendant plusieurs jours successifs, à une même heure de chaque jour pendant une période quotidienne de luminosité solaire, appelée période diurne. Pour cela, un instant de mesure t est recherché pendant la période diurne de chaque jour, par exemple entre 10 heures et 16 heures, pour lequel Pi(ti)-Pi (ti_i) 0, OÙ t11, tj et t1+1 désignent trois instants successifs de la période diurne auxquels la production d'électricité de l'installation surveillée 1 a été mesurée. Cette étape est désignée par la référence S41 dans la figure 5. On recherche alors si la même heure ti de minimum de production électrique mesurée se reproduit pour plusieurs jours successifs, éventuellement avec un 20 écart de tolérance pour les instants 4 qui sont relatifs à des jours différents, jusqu'à atteindre un nombre minimal de jours prédéfini. Une anomalie est alors encore déclarée si ce nombre minimal de jours, formant un quatrième seuil, est atteint (étape S40). Un tel quatrième critère peut être appliqué à une période de test qui est comprise dans la période d'apprentissage, mais aussi qui 25 s'étend après celle-ci. Une cause probable d'anomalie peut alors être identifiée à partir de celui ou ceux des critères de détection d'anomalie qui sont satisfaits, et éventuellement en prenant aussi en compte des informations supplémentaires qui ont été obtenues à partir des valeurs mesurées de production d'électricité. 30 Une telle identification de cause probable d'anomalie peut être effectuée automatiquement, notamment d'après un arbre d'apprentissage qui établit une corrélation entre les critères satisfaits et des causes d'anomalies possibles. - 17 - A titre d'illustration, le premier critère peut indiquer que les valeurs Pi(t) de production d'électricité qui ont été mesurées pour l'installation photovoltaïque surveillée sont en dehors de l'intervalle de confiance à 95% déterminé pour la fonction d'estimation Pi(t). Le deuxième critère peut indiquer 5 un défaut de production électrique de 3 kW (kilowatt) entre 12 heures et 19 heures, par rapport aux 20% des installations photovoltaïques voisines qui sont les moins efficaces. Et le troisième critère peut corroborer les deux précédents en indiquant que cette durée de 12 heures à 19 heures correspond à une seule durée de marche pour laquelle la différence entre valeurs 10 mesurées et valeurs estimées pour l'installation surveillée présente une valeur de moyenne temporelle qui est négative. Enfin, le quatrième critère peut ne révéler aucun minimum pendant cette durée pour la production d'électricité mesurée de l'installation surveillée. L'apprentissage pour corréler des causes d'anomalies probables à ces résultats pour les quatre critères peut alors 15 proposer qu'une accumulation de détritus se soit produite sur les panneaux de l'installation surveillée, et nécessite qu'une vérification soit diligentée sur site, avec un nettoyage éventuel des panneaux. De façon générale, le quatrième critère est particulièrement adapté pour révéler la présence d'un élément environnemental à proximité de 20 l'installation surveillée, qui produit une ombre sur ses panneaux pendant une durée limitée à l'intérieur de la période d'ensoleillement quotidienne. Typiquement, une telle ombre peut être produite nouvellement par un arbre voisin à partir du printemps. Une intervention d'élagage de l'arbre peut alors être diligentée. 25 D'autres causes d'anomalies de fonctionnement de l'installation photovoltaïque qui est surveillée, qui peuvent être détectées et identifiées grâce à l'invention, sont sans limitation : un ombrage par une cheminée d'un édifice voisin nouvellement construit, ou par pylône de téléphonie mobile nouvellement implanté, un défaut dans la transmission ou la collecte des 30 valeurs mesurées, un dysfonctionnement du panneau photovoltaïque, un dysfonctionnement du système de conversion en courant alternatif de la production d'électricité de l'installation surveillée, etc. De telles causes - 18 - d'anomalies peuvent constituer des résultats de l'arbre d'apprentissage. Il est entendu que l'invention peut être reproduite en en modifiant des aspects secondaires par rapport à la description détaillée qui vient d'être fournie. En particulier, des variations peuvent être introduites dans la méthode 5 de détermination qui est utilisée pour la fonction d'estimation de la production d'électricité d'une installation photovoltaïque à partir des valeurs qui sont mesurées pour des installations voisines. En particulier, REVENDICATIONS1. Procédé de détection à distance d'une anomalie qui affecte un fonctionnement d'une installation photovoltaïque de production d'électricité, dite installation surveillée (1), comprenant au moins une étape lors de laquelle des valeurs de production d'électricité mesurées pour l'installation surveillée sont comparées avec des valeurs qui sont déduites de mesures effectuées pour plusieurs autres installations photovoltaïques de production d'électricité, dites installations voisines (2-10) et différentes de l'installation surveillée. 2. Procédé selon la revendication 1, comprenant en outre de rechercher si plusieurs critères de détection d'anomalie sont satisfaits par les valeurs de production d'électricité qui sont mesurées pour l'installation surveillée (1), et en fonction de ceux desdits critères qui sont satisfaits, proposer au moins une cause possible pour l'anomalie, notamment en utilisant une base d'apprentissage qui associe des causes d'anomalies différentes à des combinaisons différentes des critères satisfaits. 3. Procédé selon la revendication 1 ou 2, comprenant en outre de commander ou réaliser une intervention de maintenance ou de modification environnementale pour l'installation surveillée (1), sur un site ou dans un voisinage de ladite installation surveillée, si une anomalie est détectée. 4. Procédé selon l'une quelconque des revendications 1 à 3, comprenant les étapes suivantes : /1/ pendant une période d'apprentissage, collecter les valeurs de production d'électricité mesurées pour l'installation surveillée (1) et pour les installations voisines (2-10) ; /2/ à partir des valeurs de production d'électricité mesurées pour lesdites installations voisines (2-10), déterminer une fonction d'estimation de la production d'électricité pour l'installation surveillée (1), ayant le temps comme variable, en minimisant un écart total calculé sur la période d'apprentissage entre les valeurs de production d'électricité mesurées - 20 - pour l'installation surveillée et des valeurs estimées obtenues par la fonction d'estimation ; /3/ pour au moins un instant postérieur à la période d'apprentissage, collecter de nouvelles valeurs de production d'électricité mesurées 5 pour l'installation surveillée (1) et pour chaque installation voisine (2- 10) ; /4/ pour chaque instant postérieur à la période d'apprentissage qui est concerné par l'étape /3/, calculer au moins une valeur estimée de production d'électricité en utilisant la fonction d'estimation déterminée à 10 l'étape /2/ ; et /5/ déclarer qu'une anomalie est détectée dans le fonctionnement de l'installation surveillée (1) à partir d'une comparaison, pour au moins un instant postérieur à la période d'apprentissage, entre la valeur estimée de production d'électricité qui a été calculée à l'étape /4/ et la nouvelle 15 valeur de production d'électricité mesurée qui a été collectée à l'étape /3/ pour l'installation surveillée. 5. Procédé selon la revendication 4, suivant lequel la fonction d'estimation déterminée à l'étape /2/ comprend une combinaison linéaire de fonctions de production d'électricité qui sont relatives respectivement aux 20 installations voisines (2-10), et qui correspondent aux valeurs de production d'électricité mesurées pour lesdites installations voisines pendant la période d'apprentissage, avec des facteurs de combinaison linéaire qui sont tous positifs ou nuls, et avec une somme desdits facteurs de combinaison linéaire qui est inférieure ou égale à une valeur maximale fixée, lesdits facteurs de 25 combinaison linéaire étant déterminés pour minimiser l'écart total calculé sur la période d'apprentissage entre les valeurs de production d'électricité mesurées pour l'installation surveillée (1) et les valeurs estimées. 6. Procédé selon la revendication 4 ou 5, suivant lequel la comparaison de l'étape /5/ est effectuée en exécutant les sous-étapes suivantes : - 21 - /5-1/ déterminer un intervalle de confiance de prédiction pour chaque valeur estimée de production d'électricité qui a été calculée à l'étape /4/ ; /5-2/ déterminer si la nouvelle valeur de production d'électricité mesurée 5 qui a été collectée à l'étape /3/ pour l'installation surveillée (1), et qui correspond à un même instant que la valeur estimée de la sous-étape /5-1/, est à l'extérieur de l'intervalle de confiance de prédiction relatif à ladite valeur estimée ; et /5-3/ répéter les sous-étapes /5-1/ et /5-2/ pour plusieurs instants à 10 l'intérieur d'une première période de test, puis compter un premier nombre d'instants dans la première période de test, pour chacun desquels instants la nouvelle valeur de production d'électricité mesurée pour l'installation surveillée (1) est à l'extérieur de l'intervalle de confiance de prédiction correspondant, 15 et suivant lequel une anomalie dans le fonctionnement de l'installation surveillée (1) est déclarée si le premier nombre d'instants est supérieur à un premier seuil prédéfini. 7. Procédé selon l'une quelconque des revendications précédentes, comprenant en outre les étapes suivantes : 20 /6/ corriger les valeurs de production d'électricité mesurées pour chaque installation voisine (2-10) par rapport à une puissance électrique nominale de ladite installation voisine, et proportionnellement à une puissance électrique nominale de l'installation surveillée (1) ; /7/ sélectionner un sous-ensemble parmi les installations voisines (2-10), 25 correspondant à une proportion prédéterminée desdites installations voisines et formé par celles desdites installations voisines qui possèdent des productions d'électricité les plus basses, déterminées à partir des valeurs de production d'électricité mesurées puis corrigées ; et /8/ compter un second nombre d'instants dans une deuxième période de test, pour chacun desquels instants la valeur de production d'électricité - 22 - qui est mesurée pour l'installation surveillée (1) est inférieure ou égale à une valeur moyenne des valeurs de production d'électricité mesurées puis corrigées du sous-ensemble d'installations voisines (2-10), et suivant lequel une anomalie dans le fonctionnement de l'installation 5 surveillée (1) est aussi déclarée si le second nombre d'instants est supérieur à un deuxième seuil prédéfini. 8. Procédé selon l'une au moins des revendications précédentes, dont la revendication 4, comprenant en outre les étapes suivantes : /9/ pour plusieurs instants à l'intérieur d'une troisième période de test, 10 calculer des résultats de soustraction d'une valeur de production d'électricité qui est calculée en utilisant la fonction d'estimation déterminée à l'étape /2/, à la valeur de production d'électricité mesurée pour l'installation surveillée (1), lesdites valeurs pour calculer chaque résultat de soustraction étant relatives à un même instant dans la 15 troisième période de test ; /10/ déterminer des durées de marche, pendant chacune desquelles durées de marche les résultats de soustraction présentent une valeur de moyenne temporelle qui est différente de celles de la durée de marche précédente et de la durée de marche suivante, et avec une 20 variance desdits résultats de soustraction qui est minimale à l'intérieur de chaque durée de marche ; et /11/ compter un nombre de durées de marche à l'intérieur de la troisième période de test, pour chacune desquelles durées de marche la valeur de moyenne temporelle des résultats de soustraction est négative, 25 et suivant lequel une anomalie dans le fonctionnement de l'installation surveillée (1) est aussi déclarée si le nombre de durées de marche à valeurs de moyenne temporelle négatives pour les résultats de soustraction est supérieur à un troisième seuil prédéfini. 9. Procédé selon l'une quelconque des revendications précédentes, 30 comprenant en outre les étapes suivantes : - 23 - /12/ rechercher s'il existe un instant diurne auquel un minimum se produit pour les valeurs de production d'électricité qui sont mesurées pour l'installation surveillée (1), pour au moins un jour de fonctionnement de ladite installation surveillée ; et 5 /13/ rechercher si un minimum se produit aussi pendant plusieurs jours successifs après le jour identifié à l'étape /12/, pour les valeurs de production d'électricité qui sont mesurées pour l'installation surveillée (1), chaque fois dans un intervalle de temps prédéfini autour de l'instant diurne identifié à ladite étape /12/, et compter le nombre de 10 jours successifs pour lesquels un tel minimum se produit, et suivant lequel une anomalie dans le fonctionnement de l'installation surveillée (1) est aussi déclarée si le nombre de jours successifs compté à l'étape /13/ est supérieur à un quatrième seuil prédéfini. 10. Système de surveillance d'au moins une installation photovoltaïque, 15 dite installation surveillée (1), comprenant : - des moyens de communication (100) pour collecter des valeurs de production d'électricité qui ont été mesurées pour ladite installation surveillée (1) et pour d'autres installations photovoltaïques, dites installations voisines (2-10), différentes de l'installation surveillée ; et 20 - des moyens de calcul (101) pour mettre en oeuvre un procédé de détection à distance d'une anomalie qui affecte un fonctionnement de l'installation surveillée (1), ledit procédé étant conforme à l'une quelconque des revendications précédentes.