Brevet FR3135579A1 — publié le 20231117

Procédé et dispositif de détermination de perte de puissance d’un ensemble photovoltaïque plan L’invention concerne un procédé de détermination de perte de puissance d’un ensemble plan de module(s) photovoltaïque(s) comprenant les étapes suivantes, pour un état courant : - détection (42) de la présence/absence de défaut au sein de chaque module ; - formation (44) de deux groupes de modules sain(s) et défectueux, - pour chaque module(s) sain(s), estimation (46) de l’irradiance et de la température ambiante; - estimation (48) de la température de surface de référence via un modèle polynomial bidimensionnel prédéterminé, dont les deux ordres associés à chaque direction du plan, minimisent simultanément un critère d’erreur quadratique moyenne utilisant la température mesurée de chaque module(s) sain(s) ; - pour chaque module(s) défectueux, recherche (50) d’un nombre prédéterminé de voisins sain(s) ; - détermination (52) de la perte de puissance associée auxdits module(s) défectueux. Figure pour l'abrégé : Figure 3 Procédé et dispositif de détermination de perte de puissance d’un ensemble photovoltaïque plan La présente invention concerne un procédé de détermination de perte de puissance d’un ensemble photovoltaïque plan comprenant une ou plusieurs chaînes de modules photovoltaïques, chaque chaîne comprenant au moins un module photovoltaïque, les chaînes étant connectées en série et/ou en parallèle entre une borne positive et une borne négative. L’invention concerne également un programme d’ordinateur comportant des instructions logicielles qui, lorsqu’elles sont exécutées par un ordinateur, mettent en œuvre un tel procédé de détermination de perte de puissance d’un ensemble photovoltaïque plan. L’invention concerne également un dispositif électronique de détermination de perte de puissance d’un ensemble photovoltaïque plan comprenant une ou plusieurs chaînes de modules photovoltaïques, chaque chaîne comprenant au moins un module photovoltaïque, les chaînes étant connectées en série et/ou en parallèle entre une borne positive et une borne négative. L’invention se situe dans le domaine de la production d’électricité par des installations photovoltaïques, notamment pour le photovoltaïques intégré au bâtiment BIPV (de l’anglais Building Integrated Photovoltaïcs), et plus particulièrement dans le domaine de la surveillance du bon fonctionnement de telles installations. De manière connue, les installations photovoltaïques (également appelées centrales photovoltaïques par la suite) comprennent une ou plusieurs chaînes de modules photovoltaïques, chaque chaîne comprenant un ensemble de modules photovoltaïques, pouvant se réduire à un seul module photovoltaïque, les chaînes étant connectées en série et/ou en parallèle entre une borne positive et une borne négative, chaque module étant composé d’au moins une cellule photovoltaïque, en général une matrice de cellules, chaque cellule étant une unité configurée pour produire un courant électrique à partir de l’énergie lumineuse. Une telle installation est adaptée pour alimenter en énergie électrique une charge, par exemple un onduleur adapté pour générer des tensions et des courants alternatifs à partir d’une source produisant un courant électrique continu, relié à un réseau électrique alternatif. Un des problèmes qui se pose dans ce type de centrale photovoltaïque est la perte de puissance électrique produite relative à l’occurrence d’un ou plusieurs module(s) photovoltaïque(s) défectueux. Pour estimer une telle perte, on compare classiquement la température des modules en défaut à celle de ceux qui sont sains, la perte de puissance étant proportionnelle à la différence de température issue d’une telle comparaison. Autrement dit, en considérant deux modules photovoltaïques M1 et M2 avec une surface utile et un coefficient d’échange thermique , la relation classique entre, d’une part la différence de puissance électrique entre les deux modules, et d’autre part la différence de température entre les deux modules est exprimée selon la relation suivante : en supposant classiquement que deux modules photovoltaïques M1 et M2 ayant la même température produisent la même quantité d’énergie, et dès lors que l’un des modules présente un défaut, sa température sera plus élevée, ce qui engendre par la suite une dégradation de sa puissance électrique. Toutefois, en pratique, la difficulté inhérente à cette technique classique est de déterminer quel module photovoltaïque de l’installation considérée, correspondant à un ensemble photovoltaïque PV plan, est le plus sain pour servir de référence lors de ladite comparaison. On connaît des méthodes de modélisation de température de surface d’un ensemble photovoltaïque plan tel que le modèle Sandia notamment décrit via l’URL suivante :https://pvpmc.sandia.gov/modeling-steps/2-dc-module-iv/module-temperature/sandia-module-temperature-model/, ou selon l’approche décrite par D. Faiman dans la publication intitulée « Assessing the Outdoor Operating Temperature of Photovoltaic Modules », Prog. Photovolt . Res . Appl . , vol. 16, p. 307‑315, juin 2008, doi: 10.1002/pip.813, ou encore par M. Mattei et al. dans la publication intitulée « Calculation of the polycrystalline PV module temperature using a simple method of energy balance » Renewable Energy 31, no. 4 (April 2006): 553–67. https://doi.org/10.1016/j.renene.2005.03.010 . Ces méthodes proposent chacune un modèle pour estimer la température ambiante d’un ensemble photovoltaïque plan en tenant compte des paramètres météorologiques du site, à savoir son rayonnement solaire, le vent, etc., et les conditions de montage des modules photovoltaïques. Cependant, ces modèles proposent de manière théorique une seule température de référence pour tout la surface de l’ensemble photovoltaïque plan, ce qui n’est pas conforme à la réalité terrain dans la mesure où cette température varie dans l’espace en différents point de cette surface. En effet, pour un ensemble photovoltaïque plan BIPV, les conditions de montage des modules photovoltaïques impliquent que l’ensemble photovoltaïque plan n’est pas exposé à la même température ambiante et surtout au vent d’une manière uniforme, ce qui engendre notamment des gradients thermiques entre le centre et les extrémités de l’ensemble photovoltaïque plan. Par ailleurs, la présence de ponts thermiques dans la structure d’intégration influence aussi la température des modules dans les zones en question. D’autres approches, telles que notamment décrites par Y. Assoa et al. dans la publication intitulée « Thermal analysis of a BIPV system by various modelling approaches » de 2017, ou encore par H. Getu et al. dans la publication intitulée « Computational Fluid Dynamics (CFD) Analysis of Air Based Building Integrated Photovoltaic Thermal (BIPV/T) Systems for Efficient Performance » de 2014, utilisent une modélisation en éléments finis avec des modèles physiques de propagation de chaleur pour estimer la température de surface des modules photovoltaïques. Néanmoins ces approches sont couteuses car nécessitant des outils spécifiques, complexes et gourmandes en temps et en calculs. Par conséquent, il existe un besoin spécifique d’améliorer la modélisation de la température de surface d’un ensemble photovoltaïque plan, afin d’en améliorer ensuite la détermination de perte de puissance. A cet effet, l’invention a pour objet un procédé de détermination de perte de puissance d’un ensemble photovoltaïque plan comprenant une ou plusieurs chaînes de modules photovoltaïques, chaque chaîne comprenant au moins un module photovoltaïque, les chaînes étant connectées en série et/ou en parallèle entre une borne positive et une borne négative, le procédé étant mis en œuvre par un dispositif électronique de détermination de perte de puissance d’un ensemble photovoltaïque plan, et comprenant les étapes suivantes, pour un état courant dudit ensemble photovoltaïque plan : - détection de la présence ou absence de défaut au sein de chaque module photovoltaïque dudit ensemble photovoltaïque ; - en fonction du résultat de ladite détection, formation de deux groupes de modules photovoltaïques comprenant un premier groupe G de module(s) photovoltaïque(s) sain(s) dépourvu(s) de défaut et un deuxième groupe D de module(s) photovoltaïque(s) défectueux, - pour chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe G, estimation de l’irradiance et de la température ambiante associée audit module photovoltaïque sain considéré en fonction de données météorologiques et/ou satellites; - estimation de la température de surface de référence dudit ensemble photovoltaïque via un modèle polynomial bidimensionnel prédéterminé, dont les deux ordres respectivement associés à chaque direction du plan considéré dudit ensemble photovoltaïque plan, minimisent simultanément un critère d’erreur quadratique moyenne utilisant la température mesurée de chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe ; - pour chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe, recherche d’un nombre prédéterminé N voisin de module(s) photovoltaïque(s) sain(s) voisins les plus proches en distance physique dudit module photovoltaïque défectueux considéré ; - détermination de la perte de puissance associée à ou auxdits module(s) photovoltaïque(s) défectueux en utilisant, pour chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe D : - la température mesurée dudit module photovoltaïque défectueux considéré ; - la température de surface de référence estimée via ledit modèle polynomial prédéterminé obtenue en utilisant lesdits module(s) photovoltaïque(s) sain(s) voisins dudit module photovoltaïque défectueux considéré ; - l’irradiance estimée associée audit module photovoltaïque défectueux en utilisant les module(s) photovoltaïque(s) sain(s) voisins. Ce procédé vise ainsi à tirer profit de la température des modules photovoltaïques sains détectés sur l’ensemble photovoltaïque plan pour reproduire le profil de température de surface normale d’un tel ensemble photovoltaïque plan considéré en tenant compte des gradients et ponts thermique existant. Plus précisément, on obtient via le modèle polynomial bidimensionnel prédéterminé, proposé selon la présente invention, une température surfacique de référence lissée entre les différents points de l’ensemble photovoltaïque plan correspondants à des module(s) photovoltaïque(s) sain(s). Suivant d’autres aspects avantageux de l’invention, le procédé de détermination de perte de puissance d’un ensemble photovoltaïque plan comprend une ou plusieurs des caractéristiques suivantes, prises isolément ou suivant toutes les combinaisons techniquement possibles : - les coefficients du modèle polynomial bidimensionnel sont obtenus, préalablement à ladite estimation de la température de surface de référence, en utilisant une fonction polyfit de l’outil de calcul numérique MATLAB ® ; - ledit modèle polynomial bidimensionnel prédéterminé est défini par l’équation suivante dans un espace de coordonnées cartésiennes d’axe des abscisses et d’axe des ordonnées : avec les ordres du modèle polynomial bidimensionnel respectivement selon l’axe des abscisses et l’axe des ordonnées , les coefficients dudit modèle polynomial obtenus via ladite fonction polyfit, la température mesurée de chaque module(s) photovoltaïque(s) sain(s), de coordonnées dans l’espace de coordonnées cartésiennes d’axe des abscisses et d’axe des ordonnées , dudit premier groupe G étant un point d’accroche dudit modèle polynomial bidimensionnel selon l’équation suivante : , la température de surface de référence estimée via ledit modèle polynomial bidimensionnel, étant, pour tout module de coordonnées dudit ensemble photovoltaïque telle que : ; - ledit critère d’erreur quadratique moyenne utilisant la température mesurée de chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe G est défini par l’équation suivante : avec le nombre de module(s) photovoltaïque(s) sain(s) dudit premier groupe G, le nombre de module(s) photovoltaïque(s) dudit ensemble photovoltaïque plan selon l’axe des abscisses et le nombre de module(s) photovoltaïque(s) dudit ensemble photovoltaïque plan selon l’axe des ordonnées ; - un modèle Sandia de température est utilisé au cours de l’estimation de l’irradiance et de la température ambiante associée audit module photovoltaïque sain considéré en fonction de données météorologiques et/ou satellites ; - la température mesurée de chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe ou la température mesurée de chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe est obtenue par thermographie ; - la perte de puissance associée à ou auxdits module(s) photovoltaïque(s) défectueux est déterminée selon l’équation suivante, pour tout module photovoltaïque(s) défectueux dudit deuxième groupe D, et pour tout module sain voisin, dudit module photovoltaïque défectueux considéré, parmi le nombre prédéterminé N voisin issus de l’étape de recherche : , avec : la température mesurée dudit module photovoltaïque(s) défectueux considéré dudit deuxième groupe D, un coefficient d’échange thermique prédéterminé desdits modules dudit ensemble photovoltaïque, la surface utile de chacun desdits modules dudit ensemble photovoltaïque, et où est une puissance en condition standard prédéterminée par le constructeur de chacun desdits modules dudit ensemble photovoltaïque, et un facteur de perte de puissance lié à la température ; - la détermination de la perte de puissance associée à ou auxdits module(s) photovoltaïque(s) défectueux comprend en outre le calcul d’un médian de perte de puissance selon l’équation suivante : L’invention a également pour objet un programme d’ordinateur comportant des instructions logicielles qui, lorsqu’elles sont exécutées par un ordinateur, mettent en œuvre un procédé de détermination de perte de puissance d’un ensemble photovoltaïque plan tel que défini ci-dessus. L’invention a aussi pour objet un dispositif électronique de détermination de perte de puissance d’un ensemble photovoltaïque plan, comprenant une ou plusieurs chaînes de modules photovoltaïques, chaque chaîne comprenant au moins un module photovoltaïque, les chaînes étant connectées en série et/ou en parallèle entre une borne positive et une borne négative, le dispositif électronique de détermination de perte de puissance d’un ensemble photovoltaïque plan comprenant : - un module de détection configuré pour détecter, pour un état courant dudit ensemble photovoltaïque plan, la présence ou absence de défaut au sein de chaque module photovoltaïque dudit ensemble photovoltaïque ; - un module de formation configuré pour former, en fonction du résultat de ladite détection, deux groupes de modules photovoltaïques comprenant un premier groupe G de module(s) photovoltaïque(s) sain(s) dépourvu(s) de défaut et un deuxième groupe D de module(s) photovoltaïque(s) défectueux, - un premier module d’estimation configuré pour estimer, pour chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe, l’irradiance et la température ambiante associée audit module photovoltaïque sain considéré en fonction de données météorologiques et/ou satellites; - un deuxième module d’estimation configuré pour estimer la température de surface de référence dudit ensemble photovoltaïque via un modèle polynomial bidimensionnel prédéterminé, dont les deux ordres respectivement associés à chaque direction du plan considéré dudit ensemble photovoltaïque plan, minimisent simultanément un critère d’erreur quadratique moyenne utilisant la température mesurée de chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe ; - un module de recherche configuré pour, pour chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe, rechercher un nombre prédéterminé N voisin de module(s) photovoltaïque(s) sain(s) voisins les plus proches en distance physique dudit module photovoltaïque défectueux considéré ; - un module de détermination configuré pour déterminer la perte de puissance associée à ou auxdits module(s) photovoltaïque(s) défectueux en utilisant, pour chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe D : - la température mesurée dudit module photovoltaïque défectueux considéré ; - la température de surface de référence estimée via ledit modèle polynomial prédéterminé obtenue en utilisant lesdits module(s) photovoltaïque(s) sain(s) voisins dudit module photovoltaïque défectueux considéré ; - l’irradiance estimée associée audit module photovoltaïque défectueux en utilisant lesdits module(s) photovoltaïque(s) sain(s) voisins. Ces caractéristiques et avantages de l’invention apparaîtront plus clairement à la lecture de la description qui va suivre, donnée uniquement à titre d’exemple non limitatif, et faite en référence aux dessins annexés, sur lesquels : la représente de manière schématique une architecture photovoltaïque classique. la est une vue schématique d’un dispositif électronique de détermination de perte de puissance d’un ensemble photovoltaïque plan ; la est un organigramme d’un procédé de détection d’anomalie(s) dans une structure selon la présente invention ; la illustre schématiquement la perte de puissance détectée via la présente invention. Par la suite par « état courant » de l’ensemble photovoltaïque plan considéré, on entend l’état de l’ensemble photovoltaïque plan considéré associé à l’instant d’une mesure, notamment une mesure de température de chaque module(s) photovoltaïque(s) dudit ensemble photovoltaïque plan considéré, reçue en entrée dudit procédé selon l’invention. La représente de manière schématique une architecture d’un ensemble 10 photovoltaïque classique dont les particularités d'agencement et de fonctionnement sont les suivantes : Plusieurs chaînes Ch1, Ch2 de modules photovoltaïques, appelés également panneaux photovoltaïques ou panneaux solaires, sont connectées en parallèle, seulement deux chaînes étant représentées sur la pour des raisons de simplicité. Chaque chaîne de modules photovoltaïques comporte plusieurs modules photovoltaïques connectés en série, notamment les modules photovoltaïques M1 à M3 pour la chaine Ch1 et les modules photovoltaïques M10-M30 pour la chaîne Ch2. Un module photovoltaïque comporte plusieurs cellules photovoltaïques, non représentées, connectées en série. Les cellules photovoltaïques sont réparties en plusieurs rangées. Une rangée peut comporter une ou plusieurs cellules photovoltaïques. Chaque cellule photovoltaïque est destinée à convertir une énergie solaire en une énergie électrique. Un convertisseur, par exemple un onduleur 12, comportant plusieurs bras de commutation à base de transistors, est destiné à convertir une tension continue fournie par chaque chaîne de modules en une tension alternative. Des moyens de commande, intégrés ou non au convertisseur, sont aptes à commander l'onduleur pour effectuer la conversion de tension. Chaque module comporte avantageusement des diodes anti-retours, non représentées, positionnées de manière adaptée, par exemple en série avec les modules de chaque chaîne. Chaque module comporte avantageusement des diodes de contournement (dites de "bypass", pour contourner chacune une rangée de cellules distinctes d'un module si une cellule de cette rangée était en défaut. Sur la , une diode de contournement est représentée par module, les diodes de contournement étant référencées D1, D2, D3, D10, D20, D30. La est une vue schématique d’un dispositif électronique 20 de détermination de perte de puissance d’un ensemble photovoltaïque plan. Il est à noter que le dispositif 20 selon la présente invention est propre à être implémenté sur le site géographique dudit ensemble photovoltaïque plan, ou à distance de ce site géographique (i.e. en une position géographique distincte de celle du site d’implantation dudit ensemble photovoltaïque plan considéré). Un tel dispositif 20 comprend tout d’abord un module de détection 22 configuré pour détecter, pour un état courant dudit ensemble photovoltaïque plan, la présence ou absence de défaut au sein de chaque module photovoltaïque dudit ensemble photovoltaïque. Pour ce faire, le module de détection 22 est notamment propre à utiliser la technique décrite dans le brevet EP 3 474 440 B1. En particulier, selon la technique divulguée dans le brevet EP 3 474 440 B1, un tel dispositif électronique 20 de détermination de perte de puissance d’un ensemble photovoltaïque plan selon la présente invention est propre à recevoir en entrée une matrice T_M S de température mesurée, notamment par thermographie, des modules photovoltaïques constitutifs dudit ensemble photovoltaïque plan. De plus, un tel dispositif 20 comprend en outre un module de formation 24 configuré pour former, en fonction du résultat de ladite détection mise en œuvre par le module de détection 22 précité, deux groupes de modules photovoltaïques comprenant un premier groupe G de module(s) photovoltaïque(s) sain(s) dépourvu(s) de défaut et un deuxième groupe D de module(s) photovoltaïque(s) défectueux, notamment selon la technique divulguée dans le brevet EP 3 474 440 B1. Plus précisément, la température mesurée de chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe G ou la température mesurée de chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe D est obtenue par thermographie En outre, le dispositif 20 selon la présente invention comprend également un premier module d’estimation 26 configuré pour estimer, pour chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe G, l’irradiance et la température ambiante associée audit module photovoltaïque sain considéré en fonction de données météorologiques et/ou satellites. Par ailleurs, le dispositif 20 selon la présente invention comprend un deuxième module d’estimation 28 configuré pour estimer la température de surface de référence dudit ensemble photovoltaïque via un modèle polynomial bidimensionnel prédéterminé, dont les deux ordres respectivement associés à chaque direction du plan considéré dudit ensemble photovoltaïque plan, minimisent simultanément un critère d’erreur quadratique moyenne utilisant la température mesurée de chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe. Le dispositif 20 selon la présente invention comprend également un module de recherche 30 configuré pour, pour chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe, rechercher un nombre prédéterminé N voisin de module(s) photovoltaïque(s) sain(s) voisins les plus proches en distance physique dudit module photovoltaïque défectueux considéré. Enfin, le dispositif 20 selon la présente invention un module de détermination 32 configuré pour déterminer la perte de puissance associée à ou auxdits module(s) photovoltaïque(s) défectueux en utilisant, pour chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe D : - la température mesurée dudit module photovoltaïque défectueux considéré ; - la température de surface de référence estimée via ledit modèle polynomial prédéterminé obtenue en utilisant, comme décrit par la suite, lesdits module(s) photovoltaïque(s) sain(s) voisins dudit module photovoltaïque défectueux considéré ; - l’irradiance estimée associée audit module photovoltaïque défectueux en utilisant les module(s) photovoltaïque(s) sain(s) voisins, comme décrit par la suite. Dans l’exemple de la , le dispositif électronique 20 de détermination de perte de puissance d’un ensemble photovoltaïque plan comprend une unité de traitement d’informations 34 formée par exemple d’une mémoire 36 et d’un processeur 38 associé à la mémoire 36. Dans l’exemple de la , le module de détection 22, le module de formation 24, le premier module d’estimation 26, le deuxième module d’estimation 28, le module de recherche 30 et le module de détermination 32 sont réalisés chacun sous forme d’un logiciel, ou d’une brique logicielle, exécutable par le processeur 38. La mémoire 36 du dispositif électronique 20 de détermination de perte de puissance d’un ensemble photovoltaïque plan est alors apte à stocker un logiciel de détection, et un logiciel de formation, un premier logiciel d’estimation, un deuxième logiciel d’estimation, un logiciel de recherche, et un logiciel de détermination. Le processeur 38 est alors apte à exécuter chacun des logiciels parmi le logiciel de détection, le logiciel de formation, le premier logiciel d’estimation, le deuxième logiciel d’estimation, le logiciel de recherche, et le logiciel de détermination. En variante non représentée, le module de détection 22, le module de formation 24, le premier module d’estimation 26, le deuxième module d’estimation 28, le module de recherche 30 et le module de détermination 32 sont réalisés chacun sous forme d’un composant logique programmable, tel qu’un FPGA (de l’anglais Field Programmable Gate Array ), ou sous forme d’un circuit intégré, tel qu’un ASIC (de l’anglais Application Specific Integrated Circuit ). Lorsque le dispositif électronique 20 de détermination de perte de puissance d’un ensemble photovoltaïque plan est réalisé sous forme d’un ou plusieurs logiciels, c’est-à-dire sous forme d’un programme d’ordinateur, également appelé produit programme d’ordinateur, il est en outre apte à être enregistré sur un support, non représenté, lisible par ordinateur. Le support lisible par ordinateur est par exemple un medium apte à mémoriser des instructions électroniques et à être couplé à un bus d’un système informatique. A titre d’exemple, le support lisible est un disque optique, un disque magnéto-optique, une mémoire ROM, une mémoire RAM, tout type de mémoire non-volatile (par exemple EPROM, EEPROM, FLASH, NVRAM), une carte magnétique ou une carte optique. Sur le support lisible est alors mémorisé un programme d’ordinateur comprenant des instructions logicielles. Le fonctionnement du dispositif électronique 20 de détermination de perte de puissance d’un ensemble photovoltaïque plan va maintenant être décrit en référence à la qui illustre schématiquement un exemple de mise en œuvre, selon la présente invention, d’un procédé 40 de détermination de perte de puissance d’un ensemble photovoltaïque plan. Le procédé 40 selon la présente invention comprend tout d’abord, pour un état courant dudit ensemble photovoltaïque plan, une étape 42 selon laquelle le dispositif 20, via son module électronique 22 de détection, détecte la présence ou l’absence de défaut au sein de chaque module photovoltaïque dudit ensemble photovoltaïque. Comme indiqué précédemment une telle détection est notamment mise en œuvre en appliquant la technique décrite dans le brevet EP 3 474 440 B1. Le procédé 40 selon la présente invention comprend ensuite une étape 44, au cours de laquelle, en fonction du résultat de ladite détection, le dispositif 20, via son module électronique 24 de formation, forme deux groupes de modules photovoltaïques comprenant un premier groupe G de module(s) photovoltaïque(s) sain(s) M S dépourvu(s) de défaut et un deuxième groupe D de module(s) photovoltaïque(s) défectueux M D . Plus précisément, chaque module(s) photovoltaïque(s) sain(s) M S de coordonnées dans l’espace, tel qu’illustré sur la figure 1, de coordonnées cartésiennes d’axe des abscisses et d’axe des ordonnées présente une température mesurée, notamment par thermographie, tandis que chaque module(s) photovoltaïque(s) défectueux M D de coordonnées présente une température mesurée, notamment également par thermographie. Puis, tel qu’illustré sur le mode de réalisation de la figure 3, le procédé 40 selon la présente invention comprend ensuite une étape 46, mise en œuvre par ledit dispositif 20, via son premier module électronique 26 d’estimation, et pour, seulement, chaque module(s) photovoltaïque(s) sain(s) M S dudit premier groupe G, d’estimation de l’irradiance et de la température ambiante associée audit module photovoltaïque sain considéré en fonction de données météorologiques et/ou satellites. Selon une variante optionnelle, ladite étape 46 d’estimation de l’irradiance et de la température ambiante associée audit module photovoltaïque sain considéré est mise en œuvre en utilisant un modèle Sandia de température. Plus précisément, comme indiqué précédemment en mesurant, notamment par thermographie, la température mesurée de chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe G est obtenue. De plus, selon cette étape 46, en utilisant un modèle Sandia de température, des données météorologiques et/ou satellites, fournies notamment par une centrale météo, sont utilisées pour obtenir : - une température ambiante , - une mesure de vitesse du vent , et enfin - l’irradiance solaire venant de la centrale météo. Selon le modèle Sandia, l’équation suivante permet de déterminer la température en degré Celsius d’un module photovoltaïque considéré : , avec des coefficients expérimentaux prédéterminés, par exemple : a = -3.47, et b= -0.0594. Selon cette étape 46, l’irradiance ambiante est tout d’abord obtenue selon l’équation suivante : , et la température ambiante associée est estimée pour chaque module photovoltaïque de coordonnées (i, j) au sein du plan de toiture de l’ensemble photovoltaïque pour approcher la gradient thermique et selon l’équation suivante : . Puis, tel qu’illustré par la , le procédé 40 selon la présente invention comprend une étape 48, mise en œuvre par ledit dispositif 20, via son deuxième module électronique 28 d’estimation, pour estimer la température de surface de référence dudit ensemble photovoltaïque via un modèle polynomial bidimensionnel prédéterminé, dont les deux ordres respectivement associés à chaque direction du plan considéré dudit ensemble photovoltaïque plan, minimisent simultanément un critère d’erreur quadratique moyenne utilisant la température mesurée de chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe. En complément optionnel, les coefficients du modèle polynomial bidimensionnel sont obtenus, préalablement à ladite estimation 48 de la température de surface de référence, en utilisant une fonction polyfit de l’outil de calcul numérique MATLAB ® . Selon une variante particulière de ce complément optionnel, ledit modèle polynomial bidimensionnel prédéterminé est défini par l’équation suivante dans un espace de coordonnées cartésiennes d’axe des abscisses et d’axe des ordonnées avec les ordres du modèle polynomial bidimensionnel respectivement selon l’axe des abscisses et l’axe des ordonnées , les coefficients dudit modèle polynomial obtenus via ladite fonction polyfit, la température mesurée de chaque module(s) photovoltaïque(s) sain(s), de coordonnées dans l’espace de coordonnées cartésiennes d’axe des abscisses et d’axe des ordonnées , dudit premier groupe G étant un point d’accroche dudit modèle polynomial bidimensionnel selon l’équation suivante : , la température de surface de référence estimée via ledit modèle polynomial bidimensionnel, étant, pour tout module de coordonnées dudit ensemble photovoltaïque telle que : . Autrement dit, en se basant notamment sur la matrice de température moyenne mesurée des modules sains, on va estimer à l’aide d’un modèle polynomial de la température de surface de l’ensemble photovoltaïque plan. Un tel modèle polynomial permet avantageusement de reproduire la température de surface de l’ensemble photovoltaïque considéré, en utilisant la température des modules sains dispersés au sein dudit ensemble photovoltaïque comme des points d’accroches. Un tel modèle polynomial de la température de surface permet en outre avantageusement d’interpoler d’une part la température des modules défectueux manquants et d’autre part, d’évaluer les gradients thermiques et les phénomènes de refroidissement naturelles. Selon une variante additionnelle de ce complément optionnel, ledit critère d’erreur quadratique moyenne utilisant la température mesurée de chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe G est défini par l’équation suivante : avec, comme indiqué précédemment , le nombre de module(s) photovoltaïque(s) sain(s) dudit premier groupe G, le nombre de module(s) photovoltaïque(s) dudit ensemble photovoltaïque plan selon l’axe des abscisses et le nombre de module(s) photovoltaïque(s) dudit ensemble photovoltaïque plan selon l’axe des ordonnées . Selon la présente invention, les deux ordres N et M optimaux sont ainsi les ordres qui minimisent le RMSE entre la température observée (i.e. mesurée et la température de surface interpolée . Autrement dit, les deux ordres N et M associés au RMSE le plus faible (i.e. minimal) sont sélectionnés, et la matrice de température de surface interpolée obtenue pour tout couple de coordonnées (i,j) est utilisée comme décrit par la suite pour déterminer la perte de puissance de l’ensemble photovoltaïque considéré due à la chute de production des modules photovoltaïques défectueux. Il est à noter, qu’une telle étape 48 est mise en œuvre parallèlement, préalablement ou encore successivement à l’étape 46 précédemment décrite. Le procédé 40 selon la présente invention comprend ensuite une étape 50 mise en œuvre par le dispositif 20, via son module de recherche 30, de recherche d’un nombre prédéterminé N voisin de module(s) photovoltaïque(s) sain(s) voisins les plus proches en distance physique dudit module photovoltaïque défectueux considéré. Par exemple, pour chaque module photovoltaïque défectueux dudit deuxième groupe D, N voisin = 10 module(s) photovoltaïque(s) sain(s) voisins sont recherchés au plus proche en terme de distance physique du module photovoltaïque défectueux considéré. Le procédé 40 selon la présente invention comprend ensuite une étape 52 mise en œuvre par le dispositif 20, via son module de détermination 32, de détermination de la perte de puissance associée à ou auxdits module(s) photovoltaïque(s) défectueux en utilisant, pour chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe D : - la température mesurée dudit module photovoltaïque défectueux considéré ; - la température de surface de référence estimée via ledit modèle polynomial prédéterminé obtenue en utilisant lesdits module(s) photovoltaïque(s) sain(s) voisins dudit module photovoltaïque défectueux considéré ; - l’irradiance estimée associée audit module photovoltaïque défectueux considéré en utilisant les module(s) photovoltaïque(s) sain(s) voisins. En particulier, en complément optionnel relatif à cette étape 52, la perte de puissance associée à ou auxdits module(s) photovoltaïque(s) défectueux est déterminée selon l’équation suivante, pour tout module photovoltaïque(s) défectueux dudit deuxième groupe D, et pour tout module sain voisin, dudit module photovoltaïque défectueux considéré, parmi le nombre prédéterminé N voisin issus de l’étape de recherche précédente 50 : , avec : la température mesurée dudit module photovoltaïque(s) défectueux considéré dudit deuxième groupe D, un coefficient d’échange thermique prédéterminé desdits modules dudit ensemble photovoltaïque, par exemple tel que , la surface utile en m 2 de chacun desdits modules dudit ensemble photovoltaïque, et où est une puissance en condition standard prédéterminée par le constructeur de chacun desdits modules dudit ensemble photovoltaïque, par exemple de 1000W/m 2 , et un facteur de perte de puissance lié à la température tel que : . En complément optionnel relatif à cette étape 52, la détermination de la perte de puissance associée à ou auxdits module(s) photovoltaïque(s) défectueux comprend en outre le calcul d’un médian de perte de puissance selon l’équation suivante : , cette perte de puissance médiane étant la plus représentative de la chute de production associée. Selon un mode de réalisation préférentiel l’ensemble des étapes du procédé 40 précitées sont mise en œuvre automatiquement sans intervention humaine. Selon un aspect optionnel complémentaire non représenté sur la , le procédé 40 selon la présente invention comprend en outre une étape de restitution de la perte de puissance déterminée à l’étape 52 notamment par transmission de cette perte de puissance à un utilisateur et/ou par affichage de ladite perte de puissance sur une interface homme-machine. Autrement dit, la présente invention est propre à transformer une matrice de température mesurées par thermographie en une perte de puissance transmise et ou affichée. La illustre schématiquement la perte de puissance détectée via la présente invention. En particulier, la vue 54 représente le profil 56 de température moyenne mesurée, selon l’axe des x ou des y, des seuls modules sains de l’ensemble photovoltaïque considéré, les modules photovoltaïques défectueux n’étant pas pris en compte. La vue 58 de la illustre l’application du modèle polynomial bidimensionnel en utilisant la température des modules sains comme points d’accroches, ce qui se traduit par un profil 60, selon l’axe des x ou des y, de la température surfacique. Le résultat optimal 60 est notamment obtenu avec un ordre (M=4, N=6), pour M et N variant respectivement de 1 à 10, respectivement en x et en y, et une erreur associée : par rapport à la mesure, pour une centrale photovoltaïque prédéterminée, ladite centrale photovoltaïque correspondant à un ensemble photovoltaïque PV plan (ce résultat étant propre à changer selon la centrale photovoltaïque prédéterminée), ce qui est performant en termes de précision. La vue 62 illustre le profil 64 d’estimation de la perte de puissance associée, selon l’axe des x ou des y, les pics de ce profil 64 étant associés à la présence des modules les plus défectueux présentant la perte de puissance la plus élevée parmi l’ensemble des modules photovoltaïques considéré. L’homme du métier comprendra que l’invention ne se limite pas aux modes de réalisation décrits, ni aux exemples particuliers de la description, les modes de réalisation et les variantes mentionnées ci-dessus étant propres à être combinés entre eux pour générer de nouveaux modes de réalisation de l’invention. Ainsi, la présente invention propose un procédé et un dispositif de détermination de perte de puissance d’un ensemble photovoltaïque plan propres à quantifier les chutes de production en se basant seulement sur la donnée de la température du module et propose, pour ce faire, la construction d’un modèle de température de surface de référence adapté et validé sous réserve de la présence de modules photovoltaïques sains jouant le rôle de point d’accroche dudit modèle afin de tenir compte des phénomènes de gradients et ponts thermiques. L’application d’un tel modèle de température de surface est un maillon primordial et distinctif du procédé plus global proposé selon la présente invention, ledit procédé ayant en outre comme objectif la détection de défauts thermiques sur un ensemble photovoltaïque plan et l’extraction des modules sains afin d’estimer la chute de production associée aux modules défectueux. Procédé de détermination de perte de puissance d’un ensemble photovoltaïque plan comprenant une ou plusieurs chaînes de modules photovoltaïques, chaque chaîne comprenant au moins un module photovoltaïque, les chaînes étant connectées en série et/ou en parallèle entre une borne positive et une borne négative, le procédé étant mis en œuvre par un dispositif électronique de détermination de perte de puissance d’un ensemble photovoltaïque plan, et comprenant les étapes suivantes, pour un état courant dudit ensemble photovoltaïque plan : - détection (42) de la présence ou absence de défaut au sein de chaque module photovoltaïque dudit ensemble photovoltaïque ; - en fonction du résultat de ladite détection, formation (44) de deux groupes de modules photovoltaïques comprenant un premier groupe G de module(s) photovoltaïque(s) sain(s) dépourvu(s) de défaut et un deuxième groupe D de module(s) photovoltaïque(s) défectueux, - pour chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe G, estimation (46) de l’irradiance et de la température ambiante associée audit module photovoltaïque sain considéré en fonction de données météorologiques et/ou satellites; - estimation (48) de la température de surface de référence dudit ensemble photovoltaïque via un modèle polynomial bidimensionnel prédéterminé, dont les deux ordres respectivement associés à chaque direction du plan considéré dudit ensemble photovoltaïque plan, minimisent simultanément un critère d’erreur quadratique moyenne utilisant la température mesurée de chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe ; - pour chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe, recherche (50) d’un nombre prédéterminé N voisin de module(s) photovoltaïque(s) sain(s) voisins les plus proches en distance physique dudit module photovoltaïque défectueux considéré ; - détermination (52) de la perte de puissance associée à ou auxdits module(s) photovoltaïque(s) défectueux en utilisant, pour chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe D : - la température mesurée dudit module photovoltaïque défectueux considéré ; - la température de surface de référence estimée via ledit modèle polynomial prédéterminé obtenue en utilisant lesdits module(s) photovoltaïque(s) sain(s) voisins dudit module photovoltaïque défectueux considéré ; - l’irradiance estimée associée audit module photovoltaïque défectueux en utilisant les module(s) photovoltaïque(s) sain(s) voisins. Procédé selon la revendication 1, dans lequel les coefficients du modèle polynomial bidimensionnel sont obtenus, préalablement à ladite estimation (48) de la température de surface de référence, en utilisant une fonction polyfit de l’outil de calcul numérique MATLAB ® . Procédé selon la revendication 2, dans lequel ledit modèle polynomial bidimensionnel prédéterminé est défini par l’équation suivante dans un espace de coordonnées cartésiennes d’axe des abscisses et d’axe des ordonnées : avec les ordres du modèle polynomial bidimensionnel respectivement selon l’axe des abscisses et l’axe des ordonnées , les coefficients dudit modèle polynomial obtenus via ladite fonction polyfit, la température mesurée de chaque module(s) photovoltaïque(s) sain(s), de coordonnées dans l’espace de coordonnées cartésiennes d’axe des abscisses et d’axe des ordonnées , dudit premier groupe G étant un point d’accroche dudit modèle polynomial bidimensionnel selon l’équation suivante : , la température de surface de référence estimée via ledit modèle polynomial bidimensionnel, étant, pour tout module de coordonnées dudit ensemble photovoltaïque telle que : . Procédé selon la revendication 3, dans lequel ledit critère d’erreur quadratique moyenne utilisant la température mesurée de chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe G est défini par l’équation suivante : avec le nombre de module(s) photovoltaïque(s) sain(s) dudit premier groupe G, le nombre de module(s) photovoltaïque(s) dudit ensemble photovoltaïque plan selon l’axe des abscisses et le nombre de module(s) photovoltaïque(s) dudit ensemble photovoltaïque plan selon l’axe des ordonnées . Procédé selon l’une quelconque des revendications précédentes, dans lequel un modèle Sandia de température est utilisé au cours de l’estimation de l’irradiance et de la température ambiante associée audit module photovoltaïque sain considéré en fonction de données météorologiques et/ou satellites. Procédé selon l’une quelconque des revendications précédentes, dans lequel la température mesurée de chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe ou la température mesurée de chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe est obtenue par thermographie. Procédé selon l’une quelconque des revendications précédentes, dans lequel la perte de puissance associée à ou auxdits module(s) photovoltaïque(s) défectueux est déterminée selon l’équation suivante, pour tout module photovoltaïque(s) défectueux dudit deuxième groupe D, et pour tout module sain voisin, dudit module photovoltaïque défectueux considéré, parmi le nombre prédéterminé N voisin issus de l’étape de recherche (50) : , avec : la température mesurée dudit module photovoltaïque(s) défectueux considéré dudit deuxième groupe D, un coefficient d’échange thermique prédéterminé desdits modules dudit ensemble photovoltaïque, la surface utile de chacun desdits modules dudit ensemble photovoltaïque, et où est une puissance en condition standard prédéterminée par le constructeur de chacun desdits modules dudit ensemble photovoltaïque, et un facteur de perte de puissance lié à la température. Procédé selon la revendication 7, dans lequel la détermination (52) de la perte de puissance associée à ou auxdits module(s) photovoltaïque(s) défectueux comprend en outre le calcul d’un médian de perte de puissance selon l’équation suivante : Programme d’ordinateur comportant des instructions comportant des instructions logicielles qui, lorsqu’elles sont exécutées par un ordinateur, mettent en œuvre un procédé de détermination de perte de puissance d’un ensemble photovoltaïque plan conforme aux revendications 1 à 8. Dispositif électronique de détermination de perte de puissance d’un ensemble photovoltaïque plan comprenant une ou plusieurs chaînes de modules photovoltaïques, chaque chaîne comprenant au moins un module photovoltaïque, les chaînes étant connectées en série et/ou en parallèle entre une borne positive et une borne négative, caractérisé en ce que le dispositif électronique de détermination de perte de puissance d’un ensemble photovoltaïque plan comprend: - un module de détection (22) configuré pour détecter, pour un état courant dudit ensemble photovoltaïque plan, la présence ou absence de défaut au sein de chaque module photovoltaïque dudit ensemble photovoltaïque ; - un module de formation (24) configuré pour former, en fonction du résultat de ladite détection, deux groupes de modules photovoltaïques comprenant un premier groupe G de module(s) photovoltaïque(s) sain(s) dépourvu(s) de défaut et un deuxième groupe D de module(s) photovoltaïque(s) défectueux, - un premier module d’estimation (26) configuré pour estimer, pour chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe, l’irradiance et la température ambiante associée audit module photovoltaïque sain considéré en fonction de données météorologiques et/ou satellites; - un deuxième module d’estimation (28) configuré pour estimer la température de surface de référence dudit ensemble photovoltaïque via un modèle polynomial bidimensionnel prédéterminé, dont les deux ordres respectivement associés à chaque direction du plan considéré dudit ensemble photovoltaïque plan, minimisent simultanément un critère d’erreur quadratique moyenne utilisant la température mesurée de chaque module(s) photovoltaïque(s) sain(s) dudit premier groupe ; - un module de recherche (30) configuré pour, pour chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe, rechercher un nombre prédéterminé N voisin de module(s) photovoltaïque(s) sain(s) voisins les plus proches en distance physique dudit module photovoltaïque défectueux considéré ; - un module de détermination (32) configuré pour déterminer la perte de puissance associée à ou auxdits module(s) photovoltaïque(s) défectueux en utilisant, pour chaque module(s) photovoltaïque(s) défectueux dudit deuxième groupe D : - la température mesurée dudit module photovoltaïque défectueux considéré ; - la température de surface de référence estimée via ledit modèle polynomial prédéterminé obtenue en utilisant lesdits module(s) photovoltaïque(s) sain(s) voisins dudit module photovoltaïque défectueux considéré ; - l’irradiance estimée associée audit module photovoltaïque défectueux en utilisant les module(s) photovoltaïque(s) sain(s) voisins.