FR 3019408 A1 20151002 FR 1452543 A 20140325 Procédé de prévision de rayonnement solaire au sol au voisinage d'un site [001] Domaine technique [2] L'invention concerne le domaine de la prévision d'irradiation solaire incidente à la surface du sol (nommé irradiation ci-après) via une estimation du rayonnement solaire perçu au sol au voisinage d'un site d'intérêt tel que un ou plusieurs systèmes de conversion d'énergie solaire (nommés systèmes solaires par la suite) comme, par exemple, des systèmes photovoltaïques. [003] Etat de la technique [4] Les moyens de conversion des systèmes solaires utilisés pour la génération d'énergie électrique ont une production qui dépend de l'irradiation. Typiquement, un système solaire (comme par exemple un panneau photovoltaïque) produit une énergie électrique en fonction des conditions d'ensoleillement et de luminosité de son environnement proche. [5] Une prévision de l'irradiation au sol permet d'anticiper la production de tels systèmes solaires et donc d'estimer la quantité d'énergie électrique qui sera produite sur un horizon de temps futur. Il existe différentes approches pour prévoir l'irradiation du sol dont notamment : - les prévisions utilisant des modèles météorologiques (horizon de prévision allant de quelques heures à plusieurs jours) ; - les prévisions basées sur des mesures d'irradiation prises par un instrument de mesure au sol (tel qu'un capteur pyranométrique), des modèles mathématiques étant utilisés pour prévoir les valeurs des prochaines mesures sur la base d'un historique de mesures (horizon de prévision allant de quelques minutes à plusieurs heures) ; - les prévisions utilisant les images du ciel prises à partir du sol pour déterminer un déplacement de couvertures nuageuses et en déduire l'irradiation au sol (horizon de prévision pouvant aller jusqu'à quelques minutes) ; - les prévisions par traitement d'images satellitaires, estimant le rayonnement solaire reçu au sol à partir de données satellitaires (horizon de prévision de l'ordre de quelques minutes à quelques heures). [006] Il est entendu par horizon de prévision le temps entre le moment de la prévision et la date à laquelle se passe réellement l'évènement prévu (i.e. la valeur d'irradiation prévue). [007] Avec des prévisions par traitement d'images satellitaires, l'estimation de rayonnement solaire peut être effectuée avec une bonne précision, avec une prévision court-terme et par rapport à une surface plus vaste, sans nécessiter l'installation d'instruments de mesure au niveau du sol. [008] Généralement, les prévisions basées sur des images satellitaires comprennent au moins deux phases qui sont : - la prévision de la localisation d'un masque nuageux, et - la détermination du rayonnement solaire au sol. [9] Lors de la première phase, le mouvement d'un masque nuageux est classiquement reconnu à partir de plusieurs images satellitaires successives. La future localisation du masque peut être déterminée sur l'image satellite dans la continuité de ce mouvement. A partir de la future localisation estimée du masque nuageux, un indice de couverture nuageuse au niveau de la localisation du site d'intérêt peut être calculé. [10] Dans la deuxième phase, les images satellites sont traduites en valeurs de rayonnement au sol. A cet effet, des méthodes de traitement spécifiques peuvent être utilisées pour convertir le niveau de chaque pixel d'images satellites en une valeur brute d'irradiation solaire au sol. Typiquement, des méthodes d'estimation de rayonnement (la méthode Heliosat-2 par exemple) comprennent généralement les étapes de : - détermination, en fonction de l'emplacement du site d'intérêt, d'une valeur de rayonnement selon un modèle de ciel clair, i.e. rayonnement attendu au sol en considérant qu'il n'y a aucun nuage ; et - correction de la valeur déterminée en fonction de l'indice de couverture nuageuse calculé selon la première phase précitée. [011] Toutefois, il s'avère que pour obtenir des résultats de prévision cohérents, de nombreuses opérations de lissage spatial et de correction statistique a posteriori sont à réaliser sur les prévisions de déplacement du masque nuageux d'une part, et sur les valeurs corrigées à la deuxième phase d'autre part. Sans ces opérations, des divergences importantes peuvent apparaître entre : - les estimations de mouvement du masque nuageux et les déplacements réels des masses nuageuses, et - la valeur de rayonnement au sol convertie à partir l'image satellite et l'irradiation réelle perçue au sol. [012] Or, de telles opérations sont complexes et sont coûteuses en temps et en ressource de calcul. [013] Par ailleurs, la prévision de rayonnement solaire implique une qualité d'estimation qui dépend de différents paramètres tels que les conditions météorologiques. [014] Il existe ainsi un besoin de mettre en place un procédé de prévision de rayonnement au voisinage d'un site dans lequel les opérations de correction susmentionnées sont moindres voire nulles, et selon lequel un indice de confiance peut en outre être associé à l'estimation pour qualifier la prévision de rayonnement effectuée. [015] L'invention vient améliorer la situation en ce sens. [016] Résumé de l'invention [017] L'objet de l'invention est de remédier aux inconvénients précités en proposant notamment d'utiliser directement des données satellitaires correspondant à des rayonnements solaires au sol pour déterminer un déplacement des valeurs de rayonnement (traduisant typiquement une advection) et prévoir la valeur de rayonnement qui sera localisée au niveau du site dans un temps futur. [018] A cet effet, un premier aspect de l'invention concerne un procédé de prévision d'un rayonnement solaire au sol au voisinage d'un site. Le procédé comprend au moins les étapes de : a) réception de données satellitaires relatives à des rayonnements solaires au sol mesurés sur une zone déterminée au moins à un premier instant et un deuxième instant, le premier instant étant antérieur au deuxième instant ; b) localisation sur la zone déterminée du site et des valeurs de rayonnements solaires mesurés aux premier et deuxième instants ; c) comparaison des localisations des valeurs de rayonnements solaires mesurés au premier instant et au deuxième instant ; d) détermination d'un déplacement global, entre les premier et deuxième instants, d'au moins une partie des valeurs localisées ; e) estimation, en fonction du déplacement global déterminé, de la valeur de rayonnement solaire au sol à la localisation du site pour un troisième instant, le troisième instant étant postérieur au deuxième instant. [019] Le troisième instant peut correspondre notamment au temps de la prévision (i.e. instant pour laquelle la valeur estimée est prévue). La prévision de la valeur pour le troisième instant est réalisée au deuxième instant sur la base de l'analyse comparative entre les premier et deuxième instants. [020] La valeur estimée est relative à la valeur de rayonnement prévue au niveau du site ou, tout au moins, à son voisinage (considérant un possible léger décalage dans la localisation du site via les données satellitaires). [021] Contrairement à l'état de la technique nécessitant deux phases distinctes de traitement de données d'images satellitaires (étude de masque nuageux et conversion en valeur de rayonnement), l'ordonnancement des étapes a) à e) permet d'estimer une future valeur de rayonnement sur le site en ne traitant que les données satellitaires relatives aux valeurs de rayonnement solaires au sol. Ainsi, une seule phase de traitement de données permet : - de déterminer le déplacement global des valeurs au sol et - d'estimer l'irradiation du site d'intérêt pour un instant à venir. [22] Il convient de noter que les déplacements de valeurs de rayonnement sont corrélés aux déplacements de masques nuageux au-dessus du sol. En effet, les valeurs localisées à l'étape b) sont impactées par les masses nuageuses présentes dans le ciel, au-dessus et à proximité du site d'intérêt. [23] Aussi, l'utilisation des valeurs de rayonnement permet de déterminer l'évolution des masses nuageuses et des variations d'irradiation au sol sans recourir à l'utilisation de modèles météorologiques complexes ou autres estimations de mouvement de masques nuageux. Les opérations de corrections susmentionnées ne sont dès lors plus nécessaires. [24] L'ordonnancement des étapes a) à e) limite donc les temps et ressources calcul nécessaires pour la mise en oeuvre d'une prévision de rayonnement solaire au voisinage du site d'intérêt. [25] Selon un mode de réalisation avantageux, la localisation des valeurs de rayonnement solaires à l'étape b) peut être réalisée par identification de blocs de valeurs, chaque bloc correspondant à un ensemble de valeurs localisées dans une partie définie de la zone déterminée. A l'étape c), les blocs de valeurs identifiés peuvent être comparés les uns par rapport aux autres et à l'étape d), le déplacement global peut être déterminé à partir des blocs identifiés au premier instant et des blocs les plus proches parmi des blocs identifiés au deuxième instant. [26] Cette réalisation permet de comparer des ensembles de valeurs et de détecter le mouvement des ensembles entre au moins deux instants (technique de comparaison de blocs de valeurs, nommée « Block Matching » en anglais). La détermination du déplacement global d'ensemble de valeurs permet de limiter le nombre de calculs (en comparaison d'une localisation des valeurs à chaque instant et du déplacement de chaque valeur entre les instants) tout en conservant un suivi précis du déplacement de l'ensemble de valeurs par identification du bloc le plus proche à l'instant ultérieur. [27] Avantageusement, la comparaison des localisations de valeurs de rayonnements solaires à l'étape c) peut être réalisée dans un espace fermé autour de la localisation du site sur la zone déterminée. [28] De cette manière, une zone de recherche est définie dans un espace délimité autour du site d'intérêt. Les valeurs de rayonnement utilisées pour déterminer le déplacement global sont alors les données placées à proximité du site. Ces données sont les plus pertinentes pour déterminer un déplacement global des valeurs de rayonnement vers le site d'intérêt et ainsi déterminer précisément la valeur de rayonnement à venir pour le site. [29] En outre, limiter la zone de recherche à un espace fermé réduit le nombre de calculs nécessaires et évite de perturber la détermination du déplacement global avec des mouvements parasites de valeurs de rayonnement éloignées (qui se déplacent peut être complètement différemment en raison des conditions météorologiques localisées). [30] Selon un mode de réalisation avantageux, l'espace fermé peut être défini en fonction de paramètres prédéterminés tels que l'horizon de prévision et/ou des caractéristiques de vent (caractéristiques qualitatives et quantitatives notamment) de la zone déterminée. [031] A titre d'exemple, la taille de l'espace fermé peut être petite si l'horizon de prévision est court (de l'ordre du quart d'heure typiquement) et la vitesse de vent faible, et grande lorsque l'horizon de prévision est plus long (par exemple de quatre heures) et la vitesse de vent élevée. A titre purement qualitatif, les dimensions de l'espace fermé peuvent varier, et notamment être augmentées, lorsque des vents forts sont mesurés au niveau du site ou à proximité. [32] Ainsi, la zone de recherche est ajustée en fonction du contexte dans lequel la prévision est effectuée, permettant de limiter des perturbations de l'estimation de valeur de rayonnement au voisinage du site. [33] Selon des réalisations possibles, l'espace fermé peut être de forme polygonale ou elliptique. [34] Avantageusement, au moins une donnée de rayonnement solaire effectif au sol est reçue, la donnée de rayonnement solaire effectif étant utilisée pour étalonner ou corriger les données satellitaires reçues à l'étape a). [35] De cette manière, les données satellitaires reçues peuvent être étalonnées, réinitialisées et/ou ajustées selon des mesures d'irradiation effective au sol, permettant d'affiner les données utilisées pour l'estimation et de fait, rendre la prévision de rayonnement au sol encore plus précise et robuste. Par exemple, la donnée de rayonnement solaire effectif peut être issue d'un moyen de mesure pyranométrique ou d'une mesure de la production électrique de système solaire lui-même. [036] Selon une réalisation avantageuse, le déplacement global déterminé peut être une moyenne pondérée des déplacements entre les premier et deuxième instants d'au moins une partie des valeurs localisées. [37] La détermination d'un déplacement global à partir d'une telle moyenne pondérée peut permettre de suivre le déplacement moyen de certaines ou plusieurs valeurs, ou même d'un bloc de valeurs, sans être impacté significativement par des perturbations locales inhérentes par exemple à des vents tourbillonnants, à une dépression localisée ou autre. [38] En complément, la moyenne pondérée des déplacements à laquelle correspond le déplacement global déterminé peut être pondérée selon des critères de robustesse. [039] Les critères de robustesses peuvent par exemple être relatifs à une confiance attribuée à chacun des déplacements de valeurs déterminés ou à une réactualisation par filtrage adaptatif (par exemple via un filtre de Kalman). A titre d'exemple, la confiance attribuée au déplacement d'une valeur allant à l'opposé des autres déplacements déterminés peut être faible, la pondération de ce déplacement dans la moyenne étant alors plus faible (voire nulle) par rapport aux autres déplacements. [40] Ainsi, le déplacement global des valeurs de rayonnement peut tenir compte de paramètres météorologiques tels que les vents dominants au niveau ou à proximité du site, permettant d'ajuster l'estimation de valeur de rayonnement sur le site d'intérêt. [41] Selon un mode de réalisation avantageux, la valeur de rayonnement estimée à l'étape e) peut être égale à la valeur de rayonnement du deuxième instant localisée à l'issue d'un déplacement inverse au déplacement global déterminé, le déplacement inverse étant alors appliqué à partir de la localisation du site. [42] Dans ce mode de réalisation, la valeur est estimée en considérant que les conditions sont idéales, i.e. absence de perturbations au niveau et à proximité du site, ou tout au moins dans la zone de recherche précitée. [43] Selon un autre mode de réalisation avantageux, la valeur de rayonnement estimée à l'étape e) peut être égale à une moyenne pondérée de valeurs de rayonnement du deuxième instant situées dans une zone d'intérêt. La zone d'intérêt peut en l'occurrence circonscrire la valeur de rayonnement du deuxième instant localisée à l'issue d'un déplacement inverse au déplacement global déterminé, le déplacement inverse étant appliqué à partir de la localisation du site. [44] La valeur estimée est issue d'une valeur représentative de la zone d'intérêt. Une telle valeur permet d'obtenir une valeur égale ou tout au moins proche de la valeur réelle au sol, même si les conditions de prévision ne sont pas considérées comme étant idéales (présence à proximité du site de perturbations nuageuses et/ou venteuses typiquement). [45] En complément, une plage de valeurs peut être déterminée à l'étape e) selon une valeur minimale et une valeur maximale mesurées au sein de la zone d'intérêt au deuxième instant. Un indice de confiance peut être associé à la valeur de rayonnement estimée en fonction de l'écart entre la valeur minimale et la valeur maximale de la plage de valeurs. [46] L'indice de confiance déterminé peut permettre de qualifier la prévision. Cet indice de confiance peut donner une indication sur la qualité de la prévision. [47] En variante ou en complément à l'étape e), une distribution de probabilité des valeurs de rayonnements au sein de la zone d'intérêt peut être déterminée et l'indice de confiance associé à la valeur de rayonnement estimée peut en outre être fonction de la distribution de probabilité déterminée. [48] La distribution de probabilité peut permettre de calculer la distribution de la production d'énergie électrique à partir des sources solaires. La prévision de rayonnement peut être ainsi probabiliste, avec une distribution de probabilité de production bien déterminée. [49] Selon un mode de réalisation avantageux, les dimensions de la zone d'intérêt peuvent être déterminées en fonction d'un horizon de prévision et/ou de paramètres de perturbation. [50] La zone d'intérêt peut dès lors être elle aussi ajustée en fonction du contexte dans lequel la prévision est effectuée, permettant de limiter des perturbations de l'estimation de valeur de rayonnement au voisinage du site. [51] Selon une réalisation particulière, les valeurs de rayonnement mesurées peuvent être normalisées selon un modèle de valeurs d'irradiation correspondant à un ciel clair. [52] Les valeurs de rayonnement localisées prennent alors en compte les variations d'irradiation liées notamment à la course du soleil. Ainsi, l'estimation de valeur de rayonnement est encore affinée. [53] Selon un deuxième aspect, l'invention propose un dispositif de prévision d'un rayonnement solaire au sol au voisinage d'un site. Le dispositif proposé comprend au moins : - une interface d'entrée pour la réception de données satellitaires ; - une mémoire apte à stocker les données satellitaires ; - une unité de traitement de données connectée à l'interface d'entrée et à la mémoire. [54] L'unité de traitement de données est apte à mettre en oeuvre le procédé décrit ci-avant pour la prévision d'un rayonnement solaire au sol au voisinage d'un site. [55] Les données reçues peuvent être plus particulièrement des données d'irradiation issues d'images satellites. [056] Avantageusement, l'interface d'entrée peut en outre être apte à recevoir des données de rayonnement solaire effectif au sol. Cette donnée peut par exemple être émise par un pyranomètre placé au voisinage du site d'intérêt. Selon une autre réalisation possible, la donnée de rayonnement effectif au sol peut être déduite d'une production d'énergie d'un système solaire placé au niveau du site. [57] La présente invention vise également un programme informatique comportant des instructions pour la mise en oeuvre du procédé précédemment décrit, lorsque ce programme est exécuté par une unité de traitement de données telle qu'un processeur. [58] Ce programme peut utiliser n'importe quel langage de programmation (par exemple, un langage objet ou autre), et être sous la forme d'un code source interprétable, d'un code partiellement compilé ou d'un code totalement compilé. [59] Les figures 3 et 5 décrites en détails ci-après, peuvent former des organigrammes de l'algorithme général d'un tel programme informatique. [060] Brève description des figures [061] D'autres caractéristiques et avantages de l'invention apparaîtront à l'examen de la description détaillée ci-après, et des figures annexées sur lesquelles : - la figure 1 illustre un exemple de dispositif de prévision d'un rayonnement solaire au voisinage d'un site dans un mode de réalisation de l'invention ; - les figures 2 (figures 2A et 2B) illustrent des exemples d'interprétation des données satellitaires reçues par le dispositif ; - la figure 3 représente un organigramme composé d'un exemple de succession d'étapes du procédé de prévision de l'invention ; - les figures 4 (figures 4A à 4F) illustrent des exemples de cartographies de valeurs de rayonnement sur lesquelles sont appliquées des étapes du procédé ; - la figure 5 représente un organigramme composé d'exemples de modes de réalisation des étapes du procédé ; - les figures 6 (figures 6A à 6F) sont des exemples d'illustrations de cartographies de valeurs de rayonnement sur lesquelles sont appliquées ces modes de réalisation ; et - la figure 7 représente des courbes de valeurs de rayonnement solaire estimées et mesurées selon le procédé de prévision proposé. [062] Pour des raisons de clarté, les dimensions des différents éléments représentés sur ces figures ne sont pas en proportion avec leurs dimensions réelles. Sur les figures, des références identiques correspondent à des éléments identiques pour les différents modes de réalisation exposés. [63] Description détaillée [64] On se réfère tout d'abord à la figure 1 qui illustre un exemple de dispositif de prévision d'un rayonnement solaire au sol au voisinage d'un site dans un mode de réalisation de l'invention. [065] Dans ce mode de réalisation, le dispositif comporte un ordinateur 100, comprenant une mémoire 102 pour stocker des instructions permettant la mise en oeuvre du procédé, de données satellitaires reçues, et de données temporaires pour réaliser les différentes étapes du procédé tel que décrit précédemment et détaillé plus loin. [066] L'ordinateur comporte en outre une unité de traitement de données 104. Cette unité de traitement de données peut être un circuit comme par exemple : - un processeur apte à interpréter des instructions sous la forme de programme informatique, ou - une carte électronique dont les étapes du procédé de l'invention sont décrites dans le silicium, ou encore - une puce électronique programmable comme une puce FPGA (pour « Field-Programmable Gate Array » en anglais). [067] Cet ordinateur comporte une interface d'entrée 106 pour la réception de données satellitaires, et une interface de sortie 108 pour la fourniture à un dispositif distant 110 de commandes de contrôle ou de résultats d'estimation selon le procédé de prévision. Enfin, l'ordinateur peut comporter, pour permettre une interaction aisée avec un utilisateur, un écran 112 et un clavier 114. Bien entendu, le clavier est facultatif, notamment dans le cadre d'un ordinateur ayant la forme d'une tablette tactile, par exemple. [68] Typiquement, le dispositif distant 110 peut être une plateforme de service ou un terminal de client apte à recevoir les estimations envoyées par le dispositif via l'interface 106. A cet effet, le dispositif 110 peut comporter une interface de communication apte à recevoir les données communiquées par le dispositif et une unité de traitement de données pour les interpréter. Ainsi, les données d'estimation de rayonnement reçues par le dispositif distant 110 peuvent être utilisées par ce dernier pour anticiper un rayonnement solaire à venir. A titre d'exemple purement illustratif, un utilisateur peut recevoir sur son téléphone les estimations de rayonnement solaire au niveau de la localisation de ses panneaux photovoltaïques. Ainsi, l'utilisateur peut optimiser sa consommation électrique en anticipant la production d'énergie des panneaux selon le rayonnement estimé. Des services variés peuvent être proposés à des utilisateurs en fonction de différentes interprétations possibles des estimations reçues par le dispositif distant 110. [69] L'interface 106 peut en outre recevoir des données d'un pyranomètre 109 ou d'une source solaire SS relative au rayonnement solaire effectif au sol. Cette donnée peut par exemple être une valeur brute d'irradiation mesurée par le pyranomètre 109 à proximité du site ou une déduction d'une telle valeur brute d'irradiation en fonction de la production d'énergie effective de la source solaire SS. Cette donnée de rayonnement solaire effectif peut notamment permettre de corriger les données satellites reçues par le dispositif 100, comme décrit plus loin. Une telle correction peut être effectuée une seule fois pour un étalonnage des données ou bien au fur et à mesure pour une correction rétroactive des données satellitaires reçues. [70] De surcroît, l'interface 106 peut recevoir des données qualitatives et quantitatives relatives aux vents par exemple d'un anémomètre placé au niveau ou à proximité du site. Ces données de relatives aux vents peuvent notamment être prises en compte pour les critères de perturbations susmentionnés. [71] Aux figures 2, il est représenté un exemple d'interprétation des données satellitaires reçues via l'interface 106. Les données satellitaires peuvent être des données d'irradiation issue d'images satellites. Elles sont relatives à des rayonnements solaires mesurés à proximité du site S, sur une zone déterminée ZD, lors d'un premier instant INST1 (figure 2A) et un deuxième instant INST2 (figure 2B). Comme susmentionné, le site S peut être un site pour lequel la ressource solaire disponible au sol est un enjeu. Le site S peut être un emplacement où se trouvent un ou plusieurs systèmes solaires tel que des panneaux photovoltaïques. [72] La zone ZD peut par exemple être centrée par rapport à la localisation du site S, comme illustré sur les figures. [73] Les données satellitaires reçues peuvent notamment être utilisées pour reconstituer une cartographie des valeurs de rayonnement solaires sur la zone ZD à différents instants dans le temps. A cet effet, des observations par imagerie satellitaire peuvent être traduites en irradiation incidente au sol, pouvant être notée Gsat(t,p) (où t représente le temps et p la position). [74] La cartographie obtenue possède dès lors une résolution temporelle et spatiale correspondant à celle de l'image satellite d'origine, mais dont la valeur d'un pixel correspond à une estimation de l'irradiation incidente au sol moyen sur la zone couverte par le pixel. Un pixel des pixels peut notamment correspondre au site S voire même au site et une circonscription au voisinage du site. [75] A partir de la cartographie établie, les déplacements des valeurs de rayonnement mesurées peuvent être détectés entre l'instant INST1 (figure 2A) et l'instant INST2 (figure 2B). [076] Typiquement, les données satellitaires utilisées peuvent être issues de : - satellites géostationnaires ; - bases de données stockant des données émises par des satellites ; - centres de post-traitement de données satellites ; - ou autre. [077] A titre d'exemple purement illustratif, les données satellitaires sont des données issues de la base de données Helioclim-3 construite avec la méthode Heliosat-2 (à partir d'images du capteur SEVIRI du satellite géostationnaire Meteosat Second Generation par exemple). Par ailleurs, toute autre méthode permettant d'obtenir des données relatives aux valeurs de rayonnements solaires de la zone ZD peut être envisagée. [78] On se réfère maintenant à la figure 3 qui illustre un organigramme représentant un exemple de succession d'étapes du procédé de prévision proposé. [79] A l'étape 300, les données satellitaires DAT_SAT sont reçues par l'interface 106 (aussi nommée INP) du dispositif 100. Les données DAT_SAT sont relatives aux rayonnements solaires mesurés pour la zone ZD aux instants INST1 et INST2. Il est considéré que l'instant INST1 est antérieur à l'instant INST2. [080] Comme décrit ci-avant, les données DAT_SAT peuvent être issues de conversions préalables du niveau de chaque pixel d'image satellite en une valeur brute d'irradiation solaire au sol. En fonction des données DAT SAT, k valeurs de rayonnements mesurés aux instants INST1, INST2 peuvent être obtenues, k étant un entier naturel (supérieur ou égal à zéro) : k.VALINST1, k.VALINST2. Par exemple, k peut correspondre au nombre de pixels des images satellites aux instants INST1 et INST2, une valeur VAL brute de luminance de rayonnement étant associée à chacun des k pixels de l'instant INST1 et de l'instant INST2. [81] La zone ZD peut être délimitée selon les limites de l'image satellitaire à laquelle correspondent les données DAT_SAT. Selon une autre réalisation, la zone ZD peut être délimitée dans un périmètre déterminé autour du site S. [082] Selon une étape 302, le site S et les valeurs k.VALINSTlet k.VALINST2 obtenues sont localisés sur la zone ZD. Cette localisation peut être effectuée par : - géolocalisation de la position géographique du site S sur la zone ZD et - identification du positionnement sur la zone ZD de chacune des k valeurs mesurées aux instants INST1 et INST2. [083] Selon une étape 304, les localisations déterminées pour les valeurs des instants INST1 et INST2 sont comparées entre elles (via une fonction nommée COMP_LOC par exemple). Les valeurs localisées de l'instant INST1 sont comparées aux valeurs localisées INST2 en vue de déterminer les déplacements DEP des k valeurs entre l'instant INST1 et INST2. [84] Selon une réalisation possible, k déplacements DEP des k valeurs localisées sont déterminées par identification des déplacements suivis par des valeurs de même niveau aux instants INST1 et INST2. [85] En alternative, seulement p déplacements DEP peuvent être déterminés à partir des valeurs localisées, p étant un entier positif strictement inférieur à k. Cette alternative est à utiliser de préférence pour sélectionner certaines valeurs dont le déplacement est à déterminer, ce qui évite d'une part de calculer le déplacement de toutes les valeurs de l'image satellite et, d'autre part, des erreurs de calcul lorsqu'une valeur à l'instant INST1 n'apparait pas dans les valeurs de l'instant INST2. Par exemple, il est possible de fixer un pas r permettant de ne calculer des vecteurs que tous les r pixels, réduisant ainsi le nombre de vecteurs calculés et donc le temps et les ressources de calcul. [86] Selon une autre réalisation possible, permettant notamment d'optimiser les étapes de localisation et de comparaison, les déplacements DEP des k valeurs localisées sont déterminés selon une technique d'identification de blocs de valeurs, comme décrit ci-après. [87] A une étape 306, un déplacement global DEP_G entre les instants INST1 et INST2 est déterminé. Le déplacement global DEP_G peut être déterminé à partir des k déplacements DEP précédemment déterminés, ou tout au moins d'une partie des 20 déplacements DEP. [88] Dans des conditions de ciel clair (pas de nuage), les valeurs ne se déplacent normalement pas entre les premier et deuxième instants. La valeur estimée prend alors la valeur localisée au niveau du site lors des premier et deuxième instants, avec optionnellement la prise en compte de l'évolution zénithal du soleil en appliquant un 25 modèle « ciel clair » pour la valeur estimée. [89] Selon une étape 308, la valeur de rayonnement solaire au niveau de la localisation du site est estimée pour un troisième instant INST3 en fonction du déplacement global DEP_G. En l'occurrence, l'instant INST3 est considéré comme étant à venir, i.e. postérieur à l'instant INST2. [90] L'horizon de prévision est défini selon le temps qui sépare l'instant INST2 (auquel la prévision pour l'instant INST3 est effectuée) de l'instant INST3. [91] La valeur de rayonnement pour le site S à l'instant INST3 peut dès lors être estimée en anticipant le déplacement de valeurs de l'instant INST2 vers le site S. Cette anticipation peut se baser sur une translation inverse du déplacement DEP_G à partir de la localisation du site. Cette translation inverse (-DEP_G (VALINST2)) pointera alors une valeur de l'instant INST2 qui arrivera lors du prochain instant au niveau du site S en suivant le déplacement global DEP_G. [092] On se réfère maintenant aux figures 4 qui illustrent des exemples de cartographies de valeurs de rayonnement basées sur les données satellitaires reçues et pour lesquelles sont appliquées des étapes du procédé de prévision. [93] Aux figures 4A et 4B sont respectivement illustrées des valeurs de rayonnements solaires mesurées sur la zone ZD pour l'instant INST1 et pour l'instant INST2. Le site S se trouve dans la zone ZD. Il convient de constater que les valeurs se déplacent entre l'instant INST1 (figure 4A) et l'instant INST2 (figure 4B). [94] A la figure 4C, les valeurs de l'instant INST1 sont comparées aux valeurs de l'instant INST2. Des déplacements des valeurs entre l'instant INST1 et l'instant INST2 peuvent être déterminés selon des vecteurs. A titre d'exemple, neuf vecteurs V1 à V9 sont déterminés par comparaison des valeurs entre les deux instants. Les vecteurs traduisent les déplacements de neuf valeurs ayant été localisées à l'instant INST1 et identifiés à une autre localisation à l'instant INST2. [95] A la figure 4D, le déplacement global, que l'on représente sous la forme d'un vecteur Vm (aussi nommé « vecteur de déplacement Vm » ci-après), est déterminé en fonction des déplacements des valeurs entre l'instant INST1 et l'instant INST2. Le déplacement global est associé au vecteur de déplacement Vm dans la suite de la description. Le déplacement global Vm peut par exemple correspondre à une moyenne pondérée des vecteurs de déplacement V1 à V9. [96] Il est possible de pondérer et lisser les données du déplacement global Vm par un filtrage adaptatif de type filtre Kalman par exemple. [097] A la figure 4E, l'inverse du déplacement global représenté par le vecteur -Vm, opposé du vecteur Vm, est appliqué à partir du site S pour estimer la valeur de l'instant INST2 qui, suivant le déplacement Vm, sera à la localisation du site S à l'instant INST3. [098] Comme décrit plus loin, à la figure 4F, une zone d'intérêt ZI est en outre représentée. La valeur estimée pour l'instant INST3 peut correspondre à une moyenne des valeurs situées dans cette zone ZI. Selon cette réalisation, la valeur estimée est représentative d'une zone qui se dirige vers le site S et permet d'obtenir une prévision pertinente malgré d'éventuelles perturbations météorologiques. [099] On se réfère maintenant à la figure 5 qui est un organigramme dans lequel des exemples de réalisations des étapes du procédé de prévision sont représentées. [0100] Après avoir reçues les données satellitaires relatives aux valeurs de rayonnement solaires, l'étape 302 de localisation des valeurs sur la zone ZD peut être effectuée selon deux sous-étapes 500 et 502 pour établir une zone de recherche ZR et découper n blocs B de valeurs dans cette zone de recherche ZR, n étant un entier positif. Typiquement, la zone de recherche peut être découpée à l'instant INST1 en neuf blocs de taille prédéfinie (par exemple 100 pixels x 100 pixels) localisée autour du site S. [0101] Préalablement à l'étape 302, une ou plusieurs données de rayonnement solaire effectif au sol peuvent être récupérées pour étalonner ou corriger les données satellitaires reçues. [0102] Ensuite, à l'étape 304 de comparaison des localisations de valeurs, une sous-étape 510 de comparaison des blocs peut consister à chercher pour chacun des neuf blocs précités le bloc de même taille le plus ressemblant à l'instant INST2 (technique de « Block Matching » en anglais). [0103] La zone de recherche évite de parcourir l'image entière pour la localisation et la comparaison de valeur, en privilégiant les valeurs pertinentes pour une estimation de l'irradiation à venir au niveau du site, i.e. les valeurs de rayonnement à proximité du site S. La taille et la forme de la zone de recherche peuvent être définies par une distance de recherche en pixel. Cette distance peut par exemple correspondre à la distance maximale (verticale ou horizontale) possible entre le centre du bloc d'origine (à l'instant INSTI) et le centre du bloc cherché (à l'instant INST2). Cette distance peut être choisie en fonction de la résolution spatiotemporelle des données DATA SAT, du déplacement maximal attendu entre deux instants, ou peut dépendre de perturbations météorologiques comme les conditions de vent. La zone ZR peut être polygonale (carrée, rectangle, etc.) ou elliptique (ovale ou en forme de disque par exemple). [0104] Les vecteurs de déplacement peuvent être déterminés entre les n blocs B à l'instant INSTI et les blocs les plus ressemblants à l'instant INST2. Chaque vecteur peut être obtenu en soustrayant les coordonnées des pixels centraux de bloc B d'origine de l'instant INSTI et d'un bloc cherché le plus ressemblant à l'instant INST2. Pour des raisons de simplicité de définition d'un centre de bloc, on peut choisir des tailles de bloc impaires. [0105] Le critère de ressemblance peut être choisi selon un coefficient de corrélation entre deux blocs B. [0106] A l'étape 306, le déplacement global peut être déterminé à partir d'une moyenne pondérée MP (sous-étape 520) des déplacements des blocs pour déterminer un vecteur de déplacement moyen, tel que le vecteur de déplacement Vm précité. La pondération de la moyenne MP peut être réalisée en fonction de critères de robustesse tels qu'une confiance attribuée à chacun des déplacements de valeurs déterminés ou une réactualisation par filtrage adaptatif. [0107] En considérant des conditions idéales, c'est à dire un calcul de vecteur parfait de déplacement global et une stabilité temporelle de ce vecteur, alors l'estimation de l'étape 308 consisterait simplement à donner la valeur du pixel pointé par l'opposé du vecteur déterminé. Pour pallier d'éventuelles erreurs de comparaison de blocs de valeurs et une possible évolution du vecteur entre le temps l'instant INST2 et l'instant INST3, une zone d'intérêt ZI (sous-étape 530) peut être localisée, à partir du site S, dans la direction opposée au vecteur de déplacement global déterminé. Cette zone ZI est d'une taille strictement supérieure à un pixel et circonscrit donc la valeur de rayonnement de l'instant INST2 obtenue par le déplacement inverse (-Vm) par rapport au déplacement global, appliqué à partir du site S. [0108] Cette zone correspond à la zone qui devrait se trouver sur le site d'intérêt à l'instant INST3 d'après le vecteur de déplacement estimée. La taille de cette zone peut dépendre de l'incertitude qui porte sur le vecteur. Typiquement, plus l'horizon de prévision est grand, moins l'hypothèse de persistance du vecteur calculé est valide, donc la taille (en largeur, hauteur et rayon typiquement) de la zone peut être augmentée avec l'horizon de prévision. Idem, en ce qui concerne les variabilités météorologiques susceptibles d'affecter le vecteur. [0109] De la même façon il est possible d'extraire des critères de confiance du calcul de vecteur dont peut dépendre la taille de la zone. En fonctionnement préférentiel, la taille ne tient compte que de l'horizon de prévision h par l'équation suivante : Taille ZI = 1 + h 2 [0110] La grandeur h peut correspondre à l'horizon divisé par la résolution temporelle des données DAT SAT, c'est à dire dépendant de la période d'échantillonnage temporel des cartes satellites. Dans le cas des images Hélioclim-3, cette résolution est de 15 minutes. Une prévision avec h=1 peut correspondre à une prévision à 15 minutes. [0111] La zone d'intérêt peut dès lors être elle aussi ajustée en fonction du contexte dans lequel la prévision est effectuée. [0112] Selon une sous-étape 532, les valeurs de rayonnement de l'instant INST2 comprises dans la zone d'intérêt ZI peuvent être normalisées selon un modèle de rayonnement par condition de ciel clair (sans nuage dans le ciel), noté Gcleciale, p). Il est alors possible d'obtenir un indice noté Kc permettant de s'affranchir au maximum de la dépendance de la valeur de rayonnement estimée par rapport à la course du soleil notamment selon l'équation suivante : Kc(t, p) = Gsat(t, p) Gnezchill ft , p) [0113] La moyenne d'indice Kc noté Kcmoy de la zone d'intérêt ZI peut être calculée par l'équation : Kerne)), - 1 (Kc(t, p))p E ZI taille 2 P [0114] De la même manière, on peut extraire l'indice Kc minimum et maximum (MIN MAX) de la zone d'intérêt ZI (étape 534) ou analyser l'histogramme de l'ensemble de la zone d'intérêt ZI comme une distribution de probabilité de la prévision de rayonnement solaire (sous-étape 536). En effet, la zone d'intérêt ZI contient toutes les valeurs probables du rayonnement à venir et l'histogramme de la zone d'intérêt ZI peut être interprété comme une fonction de distribution de la prévision du rayonnement. A cet effet, les valeurs en rayonnement peuvent être converties pour fournir la prévision attendue P en réutilisant le modèle ciel clair mais au temps INST3 cette fois : P (INST2INST3,p) = KC(INST2INST3,p) X GnelClair(INST3,p) [0115] On peut alors prendre la moyenne comme valeur de prévision nominale et utiliser les bornes Kc minimum et maximum pour extraire un encadrement de la prévision ou un indice de confiance IND CONF de la prévision (sous-étape 538). En effet, plus l'écart entre les bornes Kc minimum et maximum est grand, plus l'incertitude sur la prévision est grande puisque la zone supposée arriver sur le site S présente une grande variabilité. On peut donc fixer un coefficient de confiance basé sur la différence de ces deux valeurs ou sur une analyse statistique semblable issue de la distribution des valeurs dans la zone d'intérêt. [0116] On se réfère maintenant aux figures 6 sur lesquelles sont montrés des exemples d'illustrations sur lesquelles sont appliquées les étapes et sous-étapes du procédé proposé. [0117] Aux figures 6A et 6B, des cartes de rayonnement sont reçues pour des temps T- 1 (instant INST1) et T (instant INST2). La zone de recherche ZR autour du site S peut être définie selon une forme rectangulaire centrée sur le site S. Les temps T et T-1 (correspondant respectivement aux instants INST1 et INST2, normalisés selon la période d'échantillonnage des données DATA SAT). [0118] Un masque nuageux MN se déplace au sein de la zone ZD entre le temps T-1 et le temps T (respectivement selon les positions MN(T-1) et MN(T) typiquement). [0119] A la figure 6C, le déplacement des valeurs ou des n blocs B est déterminé à partir de la zone de recherche ZR. Des vecteurs V1 à V5 de déplacement sont déterminés entre les temps T-1 et T. Les vecteurs V1 à V5 représentent les déplacements entre le temps T-1 et T des valeurs correspondant au masque MN. [0120] A la figure 6D, un vecteur de déplacement Vm est déterminé par une moyenne pondérée des vecteurs V1 à V5 identifiés par rapport à la zone ZR. [0121] A la figure 6E, on détermine le déplacement inverse -Vm au temps T à partir du site S. [0122] Comme illustré à la figure 6F, la zone d'intérêt ZI peut être déterminée dans cette direction opposée au vecteur Vm. [0123] La zone ZI peut dépendre de paramètres tels que l'horizon de prévision et possiblement de la qualité du calcul du champ de vecteur. Cette zone peut être traduite en index Kc afin de s'affranchir des phénomènes cyclique du rayonnement. On peut multiplier alors les valeurs de Kc par un modèle ciel clair calculé pour l'horaire auquel la prévision est souhaitée afin de tenir compte de l'évolution zénithale du soleil. La moyenne de cette distribution peut être prise comme prévision nominale, les bornes extrêmes pouvant être utilisées pour définir un indice de confiance. [0124] A la figure 7 sont illustrées des courbes journalières d'irradiation (W/m2) évoluant dans le temps (échelle horaire). La courbe CE montre les valeurs estimées toutes les 15 minutes au cours d'une journée pour le site S selon les étapes susmentionnées. Les courbes CEmin et CEmax représentent les valeurs minimum et maximum qui ont été mesurées dans la zone d'intérêt ZI à chaque itération des étapes. La courbe CC correspond quant à elle aux valeurs de rayonnement d'un modèle de rayonnement par condition de ciel clair (sans nuage), selon lesquelles les données satellitaires ont été étalonnées et/ou corrigées. [0125] L'écart entre les valeurs de la courbe CEmin et CEmax traduisent des variabilités plus ou moins importantes constatées dans la zone ZI. L'indice de confiance IND CONF précité peut être associé à chacune des valeurs estimées de la courbe CE en fonction de l'écart mesuré entre CEmin et CEmax. Typiquement, l'indice de confiance peut être élevé pour la valeur CE au temps Ti car l'écart est entre CEmin et CEmax y est faible. A contrario, l'indice de confiance pour la valeur CE au temps T2 peut être faible car l'écart est entre CEmin et CEmax y est important. [0126] L'indice de confiance déterminé peut permettre de qualifier la prévision. Cet indice de confiance peut par exemple servir à estimer quelles sont les probabilités de rendement de systèmes solaires découlant de la prévision en question. Cet indice de confiance indique notamment une qualité de prévision. [0127] En outre, la distribution de probabilité des valeurs de rayonnements au sein de la zone d'intérêt peut être déterminée, en complément de l'indice de confiance associé à la valeur de rayonnement estimée. [0128] L'estimation proposée est donc robuste et présente des indicateurs permettant de qualifier la qualité de prévision. [0129] Bien entendu, dans des conditions de ciel clair (pas de nuage), les déplacements de valeurs susmentionnés peuvent ne pas être détectés (ils sont nuls) car les valeurs seront identiques aux premier et deuxième instants (d'autant plus que l'évolution zénithal du soleil est prise en compte avec un modèle « ciel clair »). La valeur estimée est alors de même valeur que la valeur localisée au niveau du site lors des premier et deuxième instants. [0130] L'invention a été décrite en référence à des modes de réalisations particuliers qui ne sont pas limitatifs. Bien entendu, la présente invention ne se limite pas à la forme de réalisation décrite à titre d'exemple et elle s'étend à d'autres variantes. Par exemple, le procédé peut recevoir des données satellitaires relatives à plus de deux instants (typiquement quatre ou six instants) et effectuer le suivi des déplacements des valeurs de rayonnement à ces différents instants selon les étapes b) à c) du procédé. REVENDICATIONS1. Procédé de prévision d'un rayonnement solaire au sol au voisinage d'un site (S), ledit procédé comprenant au moins les étapes de : a) réception (300) de données satellitaires relatives à des rayonnements solaires au sol mesurés sur une zone déterminée (ZD) au moins à un premier instant (INST1) et un deuxième instant (INST2), le premier instant (INST1) étant antérieur au deuxième instant (INST2) ; b) localisation (302) sur la zone déterminée (ZD) du site (S) et des valeurs de rayonnements solaires mesurés aux premier et deuxième instants ; c) comparaison (304) des localisations des valeurs de rayonnements solaires mesurés au premier instant (INST1) et au deuxième instant (INST2) ; d) détermination (306) d'un déplacement global, entre les premier et deuxième instants, d'au moins une partie des valeurs localisées ; e) estimation (308), en fonction du déplacement global déterminé (Vm), de la valeur de rayonnement solaire au sol à la localisation du site (S) pour un troisième instant (INST3), le troisième instant (INST3) étant postérieur au deuxième instant (INST2). 2. Procédé selon la revendication 1, dans lequel à l'étape b), la localisation des valeurs de rayonnement solaires est réalisée par identification (502) de blocs de valeurs (nB), chaque bloc correspondant à un ensemble de valeurs localisées dans une partie définie de la zone déterminée (ZD), dans lequel à l'étape c), les blocs de valeurs identifiés (nB) sont comparés les uns par rapport aux autres, et dans lequel à l'étape d), le déplacement global (Vm) est déterminé à partir des blocs identifiés au premier instant (INST1) et des blocs les plus proches parmi des blocs identifiés au deuxième instant (INST2).Procédé selon la revendication 1 ou 2, dans lequel à l'étape c), la comparaison des localisations de valeurs de rayonnements solaires est réalisée dans un espace fermé (ZR) autour de la localisation du site (S) sur la zone déterminée (ZD). Procédé selon la revendication 3, dans lequel l'espace fermé (EF) est défini en fonction de paramètres prédéterminés tels que l'horizon de prévision et/ou des caractéristiques de vent de la zone déterminée (ZD). Procédé selon l'une des revendications précédentes, dans lequel au moins une donnée de rayonnement solaire effectif au sol est reçue, la donnée de rayonnement solaire effectif étant utilisée pour étalonner ou corriger les données satellitaires reçues à l'étape a). Procédé selon l'une des revendications précédentes, dans lequel à l'étape d), le déplacement global déterminé (Vm) est une moyenne pondérée des déplacements, entre les premier et deuxième instants, de l'au moins une partie des valeurs localisées. Procédé selon la revendication 6, dans lequel la moyenne pondérée des déplacements à laquelle correspond le déplacement global déterminé (Vm) est pondérée selon des critères de robustesse. Procédé selon l'une des revendications précédentes, dans lequel à l'étape e), la 25 valeur de rayonnement estimée est égale à la valeur de rayonnement du deuxième instant (INST2) localisée à l'issue d'un déplacement inverse (-Vm) au déplacement global déterminé (Vm), le déplacement inverse (-Vm) étant appliqué à partir de la localisation du site (S). 30 9. Procédé selon l'une des revendications 1 à 7, dans lequel à l'étape e), la valeur de rayonnement estimée est égale à une moyenne pondérée de valeurs de rayonnement du deuxième instant (INST2) situées dans une zone d'intérêt (ZI),3. 4. 5. 6. 7. 8.laquelle zone d'intérêt (ZI) circonscrit la valeur de rayonnement du deuxième instant (t) localisée à l'issue d'un déplacement inverse (-Vm) au déplacement global déterminé (Vm), le déplacement inverse (-Vm) étant appliqué à partir de la localisation du site (S). 10. Procédé selon la revendication 9, dans lequel à l'étape e), une plage de valeurs est déterminée selon une valeur minimale et une valeur maximale mesurées au sein de la zone d'intérêt (ZI) au deuxième instant (INST2), et dans lequel un indice de confiance est associé à la valeur de rayonnement estimée en fonction de l'écart entre la valeur minimale et la valeur maximale de la plage de valeurs. 11. 12. 13. 14. 15. Procédé selon la revendication 9 ou 10, dans lequel à l'étape e), une distribution de probabilité des valeurs de rayonnements au sein de la zone d'intérêt (ZI) est déterminée, l'indice de confiance associé à la valeur de rayonnement estimée étant en outre fonction de la distribution de probabilité déterminée. Procédé selon l'une des revendications 9 à 11, dans lequel les dimensions de la zone d'intérêt (ZI) sont déterminées en fonction d'un horizon de prévision et/ou de paramètres de perturbation. Procédé selon l'une des revendications précédentes, dans lequel les valeurs de rayonnement mesurées sont normalisées selon un modèle de valeurs correspondant à un ciel clair. Produit programme informatique comportant des instructions pour la mise en oeuvre du procédé selon l'une des revendications 1 à 13, lorsque ce programme est exécuté par une unité de traitement de données telle qu'un processeur. Dispositif de prévision d'un rayonnement solaire au sol au voisinage d'un site (S), ledit dispositif comprenant au moins :- une interface d'entrée (106) pour la réception de données satellitaires ; - une mémoire (102) apte à stocker les données satellitaires ; - une unité de traitement de données (104) connectée à l'interface d'entrée (106) et à la mémoire (102) ; ladite unité de traitement de données (104) étant apte à mettre en oeuvre le procédé de prévision d'un rayonnement solaire au sol selon l'une des revendications 1 à 13. 16. Dispositif selon la revendication 15, dans lequel l'interface d'entrée (106) est en outre apte à recevoir des données de rayonnement solaire effectif au sol.