FR 2990509 A1 20131115 FR 1254377 A 20120514 B11703 - DD13284ST 1 ESTIMATION DE LA DÉRIVE D'UN CAPTEUR DE RAYONNEMENT SOLAIRE Domaine de l'invention La présente invention concerne de façon générale les capteurs d'évaluation du rayonnement solaire, tels que des pyra- nomètres ou des cellules de référence, et, plus particulièrement, l'estimation de la dérive et la calibration d'un capteur de ce type. Exposé de l'art antérieur Des capteurs d'évaluation du rayonnement solaire sont de plus en plus utilisés avec le développement des installations solaires, que ce soit pour la production de chaleur ou d'électricité. Ces capteurs servent, en particulier, à mesurer le rayonnement (watt/m2) sur les lieux des installations de production d'énergie solaire. Ces mesures sont utilisées, entre autres, pour estimer la ressource solaire ou surveiller le bon fonctionnement et la production d'énergie des centrales solaires. On connaît essentiellement deux catégories de capteurs de rayonnement: des capteurs constitués de cellules photo-20 voltaïques de référence, et des pyranomètres qui sont des capteurs spécifiques dédiés à une mesure du rayonnement du ciel B11703 - DD13284ST 2 (rayonnement solaire). Un pyranomètre est un capteur de flux thermique qui mesure le rayonnement reçu. Tous ces dispositifs de mesure de rayonnement fournissent une information (par exemple, un courant, une ten5 sion ou un mot numérique) représentative du rayonnement qu'ils reçoivent. Un problème récurrent avec l'utilisation des capteurs de rayonnement solaire est que les mesures dérivent au cours du temps. Cette dérive impose un réétalonnage ou calibration 10 régulier. Aujourd'hui, pour être réétalonnés, les capteurs sont en général démontés puis renvoyés en usine ou dans des ateliers pour être soumis à un rayonnement de référence et corriger les coefficients d'un algorithme appliqué aux mesures et fournissant 15 le rayonnement réel à partir d'une mesure brute. La norme ISO 9847 décrit des dispositifs de calibration de capteurs solaires lors des retours en atelier ou par comparaison avec d'autres capteurs spécialement apportés sur le site par un service de maintenance. Ces dispositifs ne sont pas 20 adaptés à un étalonnage ou à une correction de dérive, sur site, continue dans le temps et automatique. Les documents US 7,166,825, US 7,576,346 et US 2009/0012731 décrivent des systèmes de réétalonnage de capteurs de rayonnement montés dans des véhicules spatiaux. Les systèmes décrits 25 permettent d'étalonner le capteur en comparant la mesure faite sur l'objet concerné (un corps astral donné) à une mesure faite sur un objet de référence pour lequel le rayonnement est connu avec précision (le soleil par exemple dont le niveau de rayonnement est connu pour peu qu'on soit en dehors de l'atmosphère 30 terrestre). De tels systèmes ne sont pas applicables à des capteurs sur terre, car le rayonnement solaire ne peut plus être utilisé comme référence du fait de la variabilité de la transmission du rayonnement par l'atmosphère suivant les conditions météo- 35 rologiques. B11703 - DD13284ST 3 La complexité du réétalonnage d'un capteur de rayonnement conduit à ce que de nombreux instruments ne sont pas ré-étalonnés assez fréquemment (en pratique quelques fois par an tout au plus), ce qui nuit à la précision des mesures et à l'estimation des capacités de production des centrales solaires. De plus, la multiplication et la localisation des centrales solaires (par exemple chez des particuliers) rendant plus complexes de tels réétalonnages. Enfin, des considérations de coût peuvent faire que cette opération sera négligée dans certains cas. Résumé Ainsi, un mode de réalisation de la présente invention vise à proposer une technique d'estimation de la dérive d'un capteur de rayonnement qui pallie tout ou partie des inconvé15 nients des techniques connues. Un autre objet d'un mode de réalisation de la présente invention est de proposer une technique de calibration adaptée à un fonctionnement sur site. Un autre objet d'un mode de réalisation de la présente 20 invention est de proposer une technique d'estimation de la dérive et de calibration sans intervention extérieure. Un autre objet d'un mode réalisation de la présente invention est de permettre d'augmenter la fréquence des calibrations d'un capteur de rayonnement. 25 Un autre objet d'un mode de réalisation de la présente invention est de proposer une solution ne nécessitant aucune modification structurelle du capteur de rayonnement. Pour atteindre tout ou partie de ces objets ainsi que d'autres, on prévoit un procédé d'estimation de la dérive d'un 30 capteur de rayonnement solaire, dans lequel on prend en compte le rayonnement mesuré par ce capteur dans ses conditions d'utilisation et un modèle de rayonnement. Selon un mode de réalisation de la présente invention, le modèle fournit une estimation du rayonnement attendu en cas 35 de temps clair. B11703 - DD13284ST 4 Selon un mode de réalisation de la présente invention, le modèle tient compte de la localisation géographique du capteur, celui-ci étant situé dans l'atmosphère terrestre. Selon un mode de réalisation de la présente invention, 5 un rayonnement mesuré est comparé au rayonnement attendu, fourni par le modèle, pendant des périodes correspondant à un temps clair. Selon un mode de réalisation de la présente invention, une estimation de la dérive du capteur est calculée, de pré-10 férence quotidiennement, en tenant compte d'estimations antérieures de la dérive. Selon un mode de réalisation de la présente invention, on exploite des rapports instantanés entre des valeurs de rayonnement mesurées et des valeurs données par le modèle, ou 15 entre des valeurs données par le modèle et des valeurs de rayonnement mesurées. Selon un mode de réalisation de la présente invention, un rapport instantané est pris en compte si l'instant est compris dans une plage temporelle au cours de laquelle la 20 variation des rapports est inférieure à un seuil. Selon un mode de réalisation de la présente invention, un rapport est pris en compte si sa valeur est comprise entre deux seuils. Selon un mode de réalisation de la présente invention, 25 un rapport est pris en compte s'il correspond à une période de jour. Selon un mode de réalisation de la présente invention, l'estimation de la dérive est obtenue à partir d'une moyenne pondérée d'estimations calculées au cours des jours précédents. 30 Selon un mode de réalisation de la présente invention, la pondération tient compte de l'ancienneté des jours pris en compte. Selon un mode de réalisation de la présente invention, la pondération tient compte d'un coefficient de fiabilité 35 affecté à l'estimation du jour considéré. B11703 - DD13284ST Selon un mode de réalisation de la présente invention, la pondération prend en compte la valeur de l'estimation. On prévoit également un procédé d'étalonnage d'un capteur de rayonnement solaire, dans lequel un coefficient de 5 correction à appliquer aux mesures est obtenu à partir d'une estimation de la dérive du capteur conformément au procédé ci-dessus. On prévoit également un capteur de rayonnement solaire adapté à la mise en oeuvre du procédé d'étalonnage. On prévoit également une centrale solaire équipée d'un tel capteur de rayonnement solaire. Brève description des dessins Ces objets, caractéristiques et avantages, ainsi que d'autres seront exposés en détail dans la description suivante 15 de modes de réalisation particuliers faite à titre non limitatif en relation avec les figures jointes parmi lesquelles : la figure 1 illustre, de façon très schématique, un exemple de système d'estimation de la capacité de production d'énergie solaire à l'échelle d'un territoire ; 20 les figures 2A et 2B sont des exemples arbitraires d'allures de rayonnement solaire reçu par un capteur de rayonnement ; la figure 3 est un schéma bloc illustrant des étapes d'un mode de réalisation du procédé d'étalonnage d'un capteur de 25 rayonnement ; et les figures 4A et 4B illustrent le fonctionnement du mode de réalisation de la figure 3. Description détaillée De mêmes éléments ont été désignés par de mêmes 30 références aux différentes figures. Par souci de clarté, seuls les étapes et éléments utiles à la compréhension des modes de réalisation décrits ont été détaillés. En particulier, l'exploitation des mesures du capteur de rayonnement pour estimer les capacités de production d'énergie d'une centrale n'a 35 pas été détaillée, les modes de réalisation décrits étant B11703 - DD13284ST 6 compatibles avec l'exploitation usuelle de mesures fournies par un capteur de rayonnement solaire. De plus, la réalisation pratique d'un capteur de rayonnement et la conversion du signal de mesure n'ont pas non plus été détaillées, les modes de réalisation décrits étant, là encore, compatibles avec les capteurs de rayonnement usuels équipés de moyens de calcul. La figure 1 illustre un exemple d'installation d'exploi- tation d'énergie solaire de façon centralisée à des fins de gestion de la distribution d'énergie. De multiples centrales solaires 12, chacune constituée d'un ou de plusieurs panneaux solaires 122, et d'un système 124 de gestion et de conversion d'énergie (incluant notamment un onduleur) sont réparties sur un territoire F. Le système 124 est généralement associé à un capteur de rayonnement 2 (par exemple de type pyranomètre). Ce capteur mesure le rayonnement dans l'atmosphère terrestre. L'énergie fournie par chaque centrale 12 peut, par exemple, être injectée sur le réseau de distribution électrique du territoire F et, parallèlement, les diverses centrales 12 envoient des données (liaisons 126) à un système informatique 3 centralisant la gestion de l'énergie. Ce système 3 comprend généralement une salle de commande et d'exploitation du réseau, équipée d'un ou plusieurs ordinateurs 32, d'un ou plusieurs écrans 34 et d'une ou plusieurs bases de données 36 recevant les informations des différentes centrales. Le système 3 permet à un ou plusieurs opérateurs de gérer la distribution d'énergie à partir des informations, connues a priori ou collectées des différentes centrales, il (s) visualise(nt) sur les écrans. Ces informations incluent, 30 entre autres, la localisation géographique des centrales (carte 341), des informations 342 sur la production et la demande et, par exemple, des informations sur le rayonnement 343 reçu par une ou plusieurs centrales, permettant une comparaison par rapport à une production attendue. B11703 - DD13284ST 7 La représentation de la figure 1 est un exemple illustratif d'une application possible des modes de réalisation qui vont être décrits. Cette illustration est bien entendue simplifiée, les modes de réalisation décrits ne portant pas sur l'exploitation des mesures obtenues par les capteurs de rayonnement solaire pour prédire une capacité de production d'énergie ou gérer la distribution. L'exploitation des capacités de production fait appel à des techniques usuelles, le plus souvent indépendantes de la nature de la source d'énergie. Une particularité de la production d'énergie solaire est cependant le besoin de prévoir, notamment à partir des conditions météorologiques, la capacité de production des centrales solaires par rapport aux autres sources d'énergie. Une autre particularité est la dissémination de centrales de faible capacité de production (moins de 10 kW) sur l'ensemble d'un territoire. Ces particularités entraînent que la fiabilité des mesures de rayonnement effectuées au niveau des centrales elles-mêmes est particulièrement importante. Par ailleurs, la dissémination et le grand nombre de centrales solaires rendent plus difficile le retour des pyranomètres ou autres capteurs en atelier pour réétalonnage. On pourrait penser utiliser des dispositifs d'étalonnage portatifs et organiser des visites de contrôle des différentes centrales de façon périodique. Toutefois, cela accroît considérablement le coût de l'étalonnage et par conséquent le coût de production d'énergie. De plus, il n'est pas envisageable de visiter les centrales de petite capacité fréquemment, ce qui retarde d'autant la prise en compte des dérives éventuelles. Dans le cadre du suivi d'un groupe de centrales de production d'énergie solaire, on peut être amené à faire de la prévision de production pour l'avenir proche. La plupart des systèmes s'appuient dans ce cas sur des mesures sur le terrain pour nourrir des procédures d'intelligence artificielle. Les B11703 - DD13284ST 8 mesures de rayonnement solaires sont bien entendu prises en compte. Par conséquent, la présence de mesures aberrantes dues au mauvais étalonnage des capteurs handicape considérablement les performances des prévisions. Il existe donc un réel besoin 5 pour un système d'estimation et de correction de la dérive des capteurs, d'autant plus si le système fonctionne de façon homogène sur l'ensemble des centrales. Bien entendu, d'autres applications sont envisageables, comme par exemple, la surveillance du "ratio" de performances ("performance ratio") de 10 centrales photovoltaïques, la prévision météo, etc. L'étalonnage d'un capteur de rayonnement, initialement ou en fonctionnement, consiste à déterminer des coefficients d'une fonction affine appliquée aux valeurs mesurées. Cette fonction corrige les mesures et fournit une valeur de rayon- 15 nement. En notant GmEs la mesure effectuée et GcoR la mesure corrigée, à un instant t donné, la mesure corrigée est obtenue en appliquant une relation du type : GCOR(t) = OE.(G(t) (1) où a et 13 sont des coefficients de la fonction affine déterminée 20 pour l'étalonnage du capteur. En pratique, la dérive du coefficient 13 est négligeable et on ne se préoccupe ici que de la dérive au cours du temps du coefficient a. On a déjà proposé, par exemple dans l'article "An 25 evaluation method of PV systems" de T. Oozeki, T. Izawa, K. Otani et K. Kurokawa (Solar Energy Materials & Solar Cells 75 (2003) 687-695), une méthode d'évaluation des performances de plusieurs modules photovoltaïques interconnectés. Cette méthode évalue les puissances produites quotidiennement par la centrale 30 au cours d'un mois. On cherche ensuite le coefficient minimum par lequel il faudrait multiplier le rayonnement dit temps clair c'est-à-dire en l'absence de nuage, pour englober les mesures effectuées par le capteur. Une telle adaptation ne s'effectue pas en temps réel et requiert des mesures sur plusieurs jours. 35 Cela permet de définir une sorte de "puissance crête effective" B11703 - DD13284ST 9 de la centrale, mais cela n'implique pas d'estimation du rayonnement reçu. Les figures 2A et 2B illustrent deux exemples d'allures de rayonnement (Watt/m2) obtenues au cours du temps (h) tout au long d'une journée par un capteur solaire. Ce sont ces mesures de rayonnement qui doivent être exploitables pour ensuite estimer les capacités de production d'énergie d'une centrale solaire. Au cours d'une journée normale, c'est-à-dire, sans variations de rayonnement autres que celles liées à la course du soleil, donc à l'heure de la journée, la courbe suit approximativement une courbe en cloche. En pratique, il y a toujours des perturbations (nuages, objets, salissures (déjections d'oiseaux), etc.) faisant écran temporairement entre le capteur et le ciel. La figure 2A illustre le cas où des perturbations mineures apparaissent. La figure 2B illustre le cas où, pendant quelques heures, le capteur est partiellement ombré mais suffisamment pour faire chuter le rayonnement reçu de façon importante. Le procédé d'étalonnage proposé ci-après est basé sur 20 une estimation de la dérive du capteur exploitant des périodes pendant lesquelles le capteur est dans une condition de ciel clair (absence de nuages). On exploite alors un modèle de rayonnement fournissant une courbe théorique du rayonnement, en condition de ciel clair, 25 pour la zone dans laquelle sont positionnés le capteur et la centrale solaire. Cette courbe théorique peut prendre en compte d'autres paramètres comme la date, l'inclinaison et l'orientation du capteur. La valeur du coefficient de correction c (qui 30 représente également une estimation de la dérive) à apporter aux mesures est adaptée, pour le jour courant, en fonction d'une exploitation des mesures d'un ou, de préférence, de plusieurs jours précédents. B11703 - DD13284ST 10 La figure 3 est un schéma bloc illustrant des étapes d'un mode de réalisation du procédé d'étalonnage d'un capteur de rayonnement solaire. Les figures 4A et 4B sont des chronogrammes illustrant 5 le fonctionnement de ce procédé. Le procédé est mis en oeuvre par un circuit de traitement numérique de type microprocesseur, programmé pour mettre en oeuvre les différentes étapes qui vont être décrites. On utilise des circuits de traitement numérique et des mémoires 10 équipant habituellement, soit le capteur de rayonnement 2 lui-même, soit un éventuel dispositif (124, figure 1) de gestion qui reçoit les informations du capteur de rayonnement 2, ou plus généralement tout dispositif informatique susceptible de communiquer avec le capteur. 15 Un modèle journalier 51 de rayonnement en temps clair de la zone dans laquelle est placé le capteur est stocké (bloc 41, TEMPLATE) dans le dispositif de traitement. Ce stockage est effectué, par exemple, lors de l'installation de la centrale solaire et du capteur. En variante, le modèle est calculé à la 20 volée, ce qui permet de prendre en compte plus facilement les paramètres autres que la localisation (date, inclinaison, orientation, etc.). Les modèles théoriques permettant d'estimer le rayonnement en temps clair sur une région donnée en fonction de 25 la date, de l'heure, de la longitude, de la latitude, de l'orientation du panneau ou du capteur, de l'inclinaison, de l'albédo, etc. sont connus. On suppose que le modèle du bloc 41 tient compte de tout ou partie de ces paramètres et, de préférence, de tous. 30 De tels modèles sont par exemple décrits dans l'article "On the clear sky model of ESRA - European Solar Radiation Atlas - with respect to the Heliosat method" de C. Rigollier, O. Bauer et L. Wald, paru dans la revue Solar Energy Vol. 68, No. 1, pages 33-48, en 2000. B11703 - DD13284ST 11 Ces modèles fournissent, en fonction de la date de l'année, une courbe en cloche 51 (figure 4A) indiquant le rayonnement en Watt/m2 en fonction de l'heure (h) de la journée. La nuit, le rayonnement est quasi nul. Pour estimer la dérive et déterminer le coefficient a à appliquer aux mesures afin, qu'en cas de temps clair, le rayonnement mesuré corresponde au rayonnement du modèle, on cherche à déterminer, au cours d'une journée de mesures réelles effectuées par le capteur, des périodes dans lesquelles celui-ci était dans une situation de temps clair. Pour cela, le modèle étant établi sur un cycle journalier, on mémorise l'allure du rayonnement mesuré par le capteur au cours d'une journée (allure 53, figure 4A). En pratique, pour optimiser l'asservissement de l'étalonnage du capteur, on mémorise les mesures de la journée qui précède l'instant de l'étalonnage (bloc 42, DAY-1). On peut considérer qu'à minuit (24h), on dispose du cycle journalier complet de la journée écoulée. Le traitement qui va suivre s'applique aux données de cette journée écoulée. On compare alors (bloc 43) l'allure réelle obtenue 53 par rapport au modèle temps clair 51. L'objectif de cette comparaison est de déterminer une ou plusieurs fenêtres ou plages temporelles 54 pendant la ou lesquelles l'allure réelle 53 peut être considérée comme correspondant à une exposition en temps clair. En pratique, la comparaison s'effectue sur des valeurs numériques, dans la mesure où le capteur fournit généralement des valeurs discrètes au cours du temps (il fournit une valeur à chaque mesure). La comparaison s'effectue de préférence sur toutes les mesures prises par le capteur, ce qui améliore la fiabilité du résultat ou, en variante pour économiser des ressources de calcul ou consommer moins, sur seulement une mesure sur un nombre prédéfini de mesures (par exemple, une mesure sur deux, une mesure sur quatre, etc.). On pourra encore prévoir plusieurs mesures consécutives agrégées (par leur moyenne, leur B11703 - DD13284ST 12 médiane, etc.). Le choix du nombre de mesures pris en compte dépend, entre autres, de la fréquence de mesure du capteur. Si celui-ci prend un très grand nombre de mesures (par exemple, toutes les secondes ou plusieurs par secondes), on peut ne pas tenir compte de toutes les mesures sans perdre trop d'information. Si, par contre, la fréquence des mesures est plus faible (par exemple, toutes les minutes), il est préférable de tenir compte de toutes les mesures pour l'étalonnage. Pour chaque mesure (échantillon), on cherche à déter- miner si cela correspond ou non à un instant de temps clair (c'est-à-dire à un moment ou le soleil n'est pas masqué). Pour cela, on peut, par exemple, estimer que les moments ou la courbe des mesures rapportée à la courbe théorique par temps clair est relativement lisse (avec un seuil de variation sur l'écart-type des variations sur une fenêtre glissante) et que les valeurs mesurées sont assez hautes (et ce, sur un intervalle de temps suffisant autour du moment considéré). Les inventeurs ont constaté que l'écart entre le rayonnement réel sur une période de temps clair par rapport à un modèle de ce rayonnement fournit directement une valeur au moins approximative du coefficient a de la correction à apporter au capteur pour son étalonnage. En fait, en considérant des périodes de temps clair, c'est-à-dire pendant lesquelles le capteur devrait donner des valeurs correspondant à celles du modèle, la formule correctrice (1) ci-dessus devrait donner la valeur du modèle à partir de la valeur mesurée. En négligeant le coefficient 13 qui est, en pratique, proche de 0 donc d'autant plus négligeable que le rayonnement est important (cas du temps clair), on voit que la correction à apporter revient au rapport de la valeur donnée par le modèle sur la valeur mesurée. Ainsi, on évalue (bloc 44, d = GMES/GTEMP) le rapport de la valeur instantanée mesurée GmEs sur la valeur instantanée GTEMP donnée par le modèle. Cette quantité positive d est 35 d'autant plus faible que le ciel est couvert et est théori- B11703 - DD13284ST 13 quement égale à l'unité si le capteur est parfaitement étalonné et que le ciel est clair. La figure 4B illustre l'allure 55 obtenue sur la fenêtre temporelle 54. Cette allure donne une indication du défaut du coefficient a (en fait de son inverse). En variante, on pourra évaluer le rapport inverse 1/d, c'est-à-dire le rapport GTEmp/GmEs qui donne directement la valeur du coefficient a. Quel que soit le rapport utilisé, ce rapport représente une estimation de la dérive du capteur. Afin d'éviter la prise en compte des périodes non significatives (par exemple la nuit où le rayonnement est en principe proche de 0) ou de fortes variations (aube et crépuscule), la détermination de la fenêtre temporelle de temps clair (bloc 45, WINDOW) s'effectue à l'intérieur d'une période 56 dite de jour. Cette période de jour est choisie, à titre d'exemple et de façon arbitraire, comme correspondant à une période dans laquelle le rayonnement selon le modèle temps clair est supérieur à une valeur donnée, typiquement 50 Watt/m2. On considère en effet généralement que les périodes dans lesquelles le rayonnement est compris entre 0 et 50 Watt/m2 correspondent à l'aube et au crépuscule. Ce seuil absolu de 50 Watt/m2 est bien entendu adapté suivant les lieux. A titre de variante, la sélection de la période de jour peut s'effectuer à partir d'un calendrier et de l'indication de l'heure. Dans le cas où plusieurs fenêtres temporelles sont 25 prises en considération, on effectue, par exemple, la moyenne des valeurs obtenues sur les différentes périodes. La calibration s'effectue de préférence pendant la seconde partie de la nuit, c'est-à-dire entre minuit et l'aube du jour suivant. La simplicité de la détermination permet 30 d'effectuer une calibration quotidienne. Si aucune fenêtre temporelle acceptable n'est disponible sur une journée donnée, l'étalonnage est reporté au jour suivant. La fiabilité de la correction de l'étalonnage du capteur dépend, entre autres, de la détermination des moments de 35 temps clair. Plus la détection des moments de temps clair est B11703 - DD13284ST 14 efficace (c'est-à-dire, plus les faux positifs et les faux négatifs seront écartés), meilleure sera la détermination du coefficient a. Pour détecter les périodes de temps clair, il est 5 possible de rechercher les fenêtres temporelles pour lesquelles le rapport d (entre les mesures et le modèle théorique) montre des variations relativement faibles et un niveau suffisant. Par exemple, si le rapport d, sur la fenêtre temporelle autour d'un instant donné (par exemple, 1 heure avant et 1 heure après 10 l'instant choisi), présente un écart-type inférieur à une valeur seuil et appartient à un intervalle jugé raisonnable, l'instant choisi pourra être considéré comme un instant où le ciel est clair. La plage de plus ou moins une heure est un exemple qui pourra être modifié. 15 Les inventeurs ont constaté que le choix de la largeur de la fenêtre temporelle peut être modifié en tenant compte (entre autres) de l'échantillonnage des données considérées. Plus l'intervalle de temps entre deux mesures est court, plus la taille de la fenêtre peut être réduite. Si aucune fenêtre 20 temporelle acceptable n'est disponible sur une journée donnée, l'étalonnage est reporté au jour suivant. Les calculs sur des fenêtres temporelles permettant la détection des instants de temps clair s'effectuent de façon itérative sur chacun des échantillons. Une telle détermination 25 itérative est particulièrement adaptée à des systèmes de traitement numérique pilotés par logiciel. Selon un mode de réalisation simplifié, on considère que l'estimation du coefficient a pour la journée considérée, qui représente la valeur de correction à appliquer à la fonction 30 affine du capteur, correspond à l'inverse de la valeur médiane des coefficients d présents dans la ou les fenêtres temporelles (bloc 46, a = 1/d). L'adaptation du capteur (bloc 47, SENSOR ADAPT) est effectuée à partir du coefficient a obtenu. En pratique, l'estimation du coefficient a varie d'un 35 jour à l'autre, que ce soit en raison d'une dérive dans le B11703 - DD13284ST 15 fonctionnement du capteur cette dérive étant progressive, à cause d'une faible quantité d'instants de temps clair dans la journée, ou suite à des imprécisions dans la détection des périodes de temps clair (par exemple des périodes de temps gris nuageux mais stables tout au long d'une journée). Ainsi, on effectue de préférence une pondération de plusieurs déterminations du coefficient a en fonction d'un ou plusieurs facteurs. De préférence, l'importance donnée aux différents facteurs est elle-même pondérée. A titre d'exemple, les inventeurs considèrent que l'ancienneté de la journée est le facteur le plus important, la fiabilité de la journée venant en deuxième et le niveau absolu de la journée étant le moins important. D'autres pondérations pourront être envisagées pour 15 rendre la détermination encore plus fiable. Par exemple, on pourra tenir compte de l'écart de la journée courante par rapport à la moyenne des journées précédentes. En prenant l'exemple de l'ancienneté de la journée, de la fiabilité de la journée et du niveau absolu de la journée, on 20 pourra, par exemple, utiliser les facteurs suivants. Un facteur tient compte d'une moyenne pondérée des coefficients des jours précédents en affectant un poids d'autant plus faible au coefficient d'une journée donnée que cette journée est ancienne par rapport au jour courant. 25 Un autre facteur tient compte de la fiabilité estimée de la journée considérée. Cette fiabilité correspond au nombre de rapports d pris en considération sur la journée considérée, ce qui revient à classer les jours, ou les coefficients a, en fonction du nombre de périodes de temps clair sur lesquelles 30 s'appuie la valeur. Plus ce nombre est important, plus le poids donné au coefficient a considéré est important. Encore un autre facteur est le niveau de rayonnement obtenu. Cela revient à mettre l'accent sur les journées dans lesquelles les valeurs calculées pour le rapport d sont les plus 35 élevées. B11703 - DD13284ST 16 La détermination du poids à affecter aux coefficients a des différents jours en tenant compte des variantes ci-dessus, revient à faire le produit des facteurs pour chaque jour. Les poids sont obtenus en normalisant les résultats afin que la somme des poids soit égale à 1. En pratique, on peut choisir de ne considérer qu'un certain nombre des jours correspondant aux plus récents. En pratique, les étapes ci-dessus sont effectuées en traitant de façon successive les échantillons de la journée 10 considérée. Ce traitement peut être exprimé de la façon suivante. On note i, l'instant ou le rang de l'échantillon dans la journée, et j le jour (j=0 pour le jour courant, c'est-à-dire le jour au début duquel l'évaluation à partir des mesures 15 réalisées est calculée, j=1 pour le jour sur lequel les mesures sont faites, c'est-à-dire la veille du jour courant, j=2 pour le jour précédent, etc.). On désigne par x, une plage horaire (en nombre d'échantillons) autour de l'échantillon courant définissant un 20 intervalle (plage temporelle) pendant lequel la mesure doit être à peu près stable pour que l'on considère une absence de perturbation (nuageuse ou autre). On fixe un seuil y du rapport d(i) en deça duquel doit rester l'écart-type du rapport calculé sur la plage horaire pour que l'on considère que le rayonnement 25 est resté "stable". On fixe un seuil z de détermination de mesures trop sombres ou aberrantes. Par exemple, si d(i) est inférieur à 1-z, on considère que le ciel est "sombre" à l'instant i, si d(i) est supérieur à 1+z, on considère que le ciel est "anormalement 30 clair" à l'instant i (mesure aberrante). Suite à une longue période de dérive, un capteur peut renvoyer des valeurs hors de l'intervalle [1-z, 1+z]. Il est cependant possible de se prémunir d'une élimination inappropriée de ces valeurs en normalisant les valeurs GmEs à partir de l'estimation du 35 coefficient a de la veille. Dans ce cas, une fois la vérifi- B11703 - DD13284ST 17 cation faite, l'opération inverse sera réalisée pour ne pas perturber la suite du calcul. Un rapport d(i) est pris en compte pour le calcul de la médiane (prise en compte dans l'allure 55) si : - la plage horaire i-x à i+x est comprise dans le "jour", c'est-à-dire soit la mesure GmEs est supérieure à un seuil (par exemple 50 watt/m) dans toute cette plage, soit la plage exclut la nuit, l'aube et le crépuscule, ou les deux ; - l'écart type des rapports d(i) sur la plage i-x à i+x 10 est inférieur au seuil y (par exemple, le seuil d'écart type est choisi entre 0,03 et 0,1, typiquement de l'ordre de 0,06) ; et la valeur du rapport d(i) est comprise entre 1-z et 1+z (par exemple, le seuil z est compris entre 0,2 et 0,6, de préférence de l'ordre de 0,4). 15 Le coefficient ao est obtenu d'après la relation suivante : j=1 (2) N L7_,Pj j=1 où N représente le nombre de jours passés pris en compte (de préférence entre 10 et 30) et pj représente le poids donné au 20 jour de rang j. Il s'agit donc d'une moyenne pondérée des coefficients aj. Dans un mode de réalisation simplifié, on ne tient compte que de l'ancienneté de la journée. Par exemple, on affecte des poids pj correspondant au rang de la journée 25 (Pj = j). Dans un mode de réalisation préféré où l'on prend en compte les trois facteurs mentionnés plus haut, le poids pj est t d obtenu, en combinant les poidsp p p . . et . respectivement i, i i donnés à l'ancienneté du coefficient j, au nombre de données B11703 - DD13284ST 18 sur coefficient ai et à la valeur du où : lesquelles s'appuie le coefficient suivante : plus cj d'après la relation : (3) t d e (4) pi Pj -Pj (5) (6) de l'estimation a rangée de la (par exemple 1 pour al, 2 pour Pi N Pei Pei , avec I(Pti ( t ) i=1 ni ; t 1 0 ; et pi = exp -__ .- , d N numéro pi =exp i ancienne p, ei =exp ,désigne nj ( d récente à la ni -7.-2- N i ( e n -4.N i le plus a2, jusqu'à N pour aN) ; désigne le numéro de l'estimation aj rangée suivant le nombre de d(i) pris en compte pour son calcul (1 pour la valeur de cj basée sur le plus grand nombre d'instants i considérés comme "temps clair" le jour j, N pour la valeur de ai 15 basée sur le moins grand nombre d'instants) ; et n, désigne le numéro de l'estimation ai rangée suivant le niveau de la valeur a (1 pour la plus grande, N pour la plus faible). Les nombres 10, 7 et 4 des formules (4) à (6) sont des d 20exemplesdefacteurscPimportancedorméeauxpoids.,.et pi pi -Plus le nombre en facteur de l'exponentielle est élevé, pi plus le poids est important dans la moyenne de la formule (3). La somme des poids pi est normalisée à 1 (la moyenne étant faite par la formule (3)). Par conséquent, la formule (2) 25 devient : B11703 - DD13284ST 19 cto = (ID - atj (2' ) Un avantage des modes de réalisation qui ont été décrits est qu'il est désormais particulièrement simple d'étalonner un capteur en temps réel en fonction des dérives qu'il subit. Cela fiabilise considérablement les mesures de rayonnement solaire. Un autre avantage des modes de réalisation décrits est que leur mise en oeuvre ne nécessite aucune modification structurelle des capteurs de rayonnement solaire existants, l'esti- mation de la dérive, puis l'étalonnage s'effectuent à partir d'une interprétation des valeurs mesurées et mémorisées. On notera que l'estimation de la dérive du capteur, ou la détermination de la valeur corrective cx, peut même s'effectuer à distance du capteur, dans le dispositif 124 de gestion du panneau solaire, voire sur un serveur distant. Divers modes de réalisation ont été décrits. Diverses variantes et modifications apparaîtront à l'homme de l'art. En particulier, tout ce qui a été exposé en relation avec une mesure directe du rayonnement peut s'effectuer à partir de gran- deurs représentatives de ce rayonnement. Par exemple, il pourra s'agir du courant circulant dans l'élément de mesure du capteur ou de toute autre grandeur représentative du rayonnement instantané. De plus, la mise en oeuvre pratique des modes de réalisation décrits est à la portée de l'homme du métier à partir des indications fonctionnelles données ci-dessus en utilisant des outils informatiques et de programmation. REVENDICATIONS1. Procédé d'estimation de la dérive d'un capteur de rayonnement solaire (2), dans lequel on prend en compte le rayonnement mesuré (GMES) par ce capteur dans ses conditions d'utilisation et un modèle (51) de rayonnement. 2. Procédé selon la revendication 1, dans lequel le modèle (51) fournit une estimation du rayonnement attendu en cas de temps clair. 3. Procédé selon la revendication 1 ou 2, dans lequel le modèle (51) tient compte de la localisation géographique du 10 capteur, celui-ci étant situé dans l'atmosphère terrestre. 4. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel un rayonnement mesuré (GMES) est comparé (43) au rayonnement attendu (GTEMP), fourni par le modèle (51), pendant des périodes correspondant à un temps clair. 15 5. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel une estimation de la dérive du capteur est calculée, de préférence quotidiennement, en tenant compte d'estimations antérieures de la dérive. 6. Procédé selon l'une quelconque des revendications 1 20 à 5, dans lequel on exploite des rapports instantanés (d, a) entre des valeurs de rayonnement mesurées (GMES) et des valeurs (GTEMP) données par le modèle (51), ou entre des valeurs (GTEMP) données par le modèle (51) et des valeurs de rayonnement mesurées (GMES). 25 7. Procédé selon la revendication 6, dans lequel un rapport instantané (d) est pris en compte si l'instant est compris dans une plage temporelle (54) au cours de laquelle la variation des rapports est inférieure à un seuil. 8. Procédé selon la revendication 6 ou 7, dans lequel 30 un rapport (d) est pris en compte si sa valeur est comprise entre deux seuils. 9. Procédé selon l'une quelconque des revendications 6 à 8, dans lequel un rapport (d) est pris en compte s'il correspond à une période de jour.B11703 - DD13284ST 21 10. Procédé selon l'une quelconque des revendications 1 à 9, dans lequel l'estimation de la dérive est obtenue à partir d'une moyenne pondérée d'estimations calculées au cours des jours précédents. 11. Procédé selon la revendication 10, dans laquelle la pondération tient compte de l'ancienneté des jours pris en compte. 12. Procédé selon la revendication 10 ou 11, dans lequel la pondération tient compte d'un coefficient de fiabilité 10 affecté à l'estimation du jour considéré. 13. Procédé selon l'une quelconque des revendications 10 à 12, dans lequel la pondération prend en compte la valeur de l'estimation. 14. Procédé d'étalonnage d'un capteur de rayonnement 15 solaire, dans lequel un coefficient de correction (a) à appliquer aux mesures est obtenu à partir d'une estimation de la dérive du capteur conformément à l'une quelconque des revendications 1 à 13. 15. Capteur de rayonnement solaire adapté à la mise en 20 oeuvre du procédé d'étalonnage comportant des moyens informatiques programmés pour la mise en oeuvre du procédé selon la revendication 14. 16. Centrale solaire équipée d'un capteur de rayonnement solaire conforme à la revendication 15.