Brevet FR3071623A1 — publié le 20190329

DOMAINE TECHNIQUE DE L'INVENTION La présente invention s'inscrit dans le domaine de l'estimation du rayonnement solaire incident, par exemple dans le cadre de la construction d'un modèle énergétique d'une structure. Que ce soit pour construire un modèle prédictif du flux solaire incident attendu futur (par exemple, pour prévoir les performances opérationnelles de panneaux photovoltaïques), pour reconstituer un flux solaire reçu à des périodes passées, ou encore pour estimer en temps réel le flux solaire incident, il est très avantageux de disposer d'une mesure dynamique de flux solaire incident. ETAT DE LA TECHNIQUE La solution la plus répandue pour fournir une mesure précise de flux solaire consiste à utiliser un pyranomètre, fournissant une mesure directe de la puissance de rayonnement solaire reçue. A titre d'exemple, la publication Modeling and forecasting hourly global solar radiation using clustering and classification techniques, P. Jiménez-Perez, L. Mora-L6pez, Solar Energy, vol. 135, pp. 682-691, 2016, propose de réaliser des mesures in situ de flux solaire incident à l'aide d'un pyranomètre, servant ensuite de données d'entraînement pour construire un modèle statistique du flux solaire par apprentissage. Toutefois, ce modèle est très spécifique est n'est valide que pour le site d'implantation du pyranomètre. Le principal inconvénient d'un tel protocole est de nécessiter un pyranomètre pour chaque emplacement pour lequel on souhaite obtenir un modèle de flux solaire, le coût d'un pyranomètre étant prohibitif. En pratique, les stations météorologiques standard ne comprennent que rarement un pyranomètre. Il existe donc un besoin pour une solution d'estimation du flux solaire incident qui ne nécessite pas un pyranomètre, et qui puisse être mise en oeuvre dans des dispositifs mis à disposition d'un public large. Les publications Solar irradiance forecasting using spatial-temporal covariance structures and time-forward kriging, D. Yang, C. Gu, Z. Dong, 2 P. Jirutitijaroen, N. Chen et W. M. Walsh, Renewable Energy, pp. 235-245, 2013, et Spatio-temporal kriging of solar radiation incorporating direction and speed of cloud movement, T. Inoue, T. Sasaki et T. Washio, The 26th Annual Conference of the Japanese Society for Artificial Intelligence, 5 Yamaguchi City, 2012, décrivent, pour la construction d'un modèle du flux solaire incident sur un territoire, des méthodes de krigeage (méthodes géostatistiques d'estimation des caractéristiques d'un élément de volume) fournissant un modèle statistique spatial ou spatio-temporel. Un modèle spatio-temporel de dépendance statistique du flux solaire est appris à 10 partir d'un nombre réduit de pyranomètres fournissant des données en temps réel de flux solaire, le modèle ainsi appris pouvant ensuite prédire le flux solaire sur un emplacement non doté d'un pyranomètre. Une telle méthode de krigeage permet de s'affranchir de la nécessité d'un pyranomètre pour chaque point où l'on souhaite obtenir une donnée de 15 flux solaire. Cependant, la qualité d'un tel modèle dépend de la densité spatiale du maillage des pyranomètres pendant la phase d'apprentissage. Ainsi, un nombre élevé de pyranomètres reste nécessaire pour disposer d'un modèle dynamique pertinent pour le flux solaire incident. Des solutions ont été proposées pour fournir une estimation indirecte du 20 flux solaire incident, ne nécessitant pas de pyranomètre. Un modèle prédictif du flux solaire utilisant des données de température est décrit dans la publication Potentiel of radial basis function based support vector regression for global solar radiation prediction, Z. Ramedani, M. Omid, A. Keyhani, S. Shamshirband et B. Khoshnevisan, Renewable and 25 Sustainable Energy Reviews, vol. 39, pp. 1005-1011, 2014. Cette solution consiste à utiliser un unique capteur de température extérieure, et à relever au cours d'une journée le minimum et le maximum de température mesurés par ledit capteur. La différence de ces valeurs est exploitée pour prédire, par apprentissage artificiel, des valeurs journalières ou 30 mensuelles de flux solaire. Toutefois, un tel modèle prédictif ne permet pas une modélisation dynamique du flux solaire avec la précision temporelle nécessaire. 3 Par ailleurs, il est connu d'évaluer une approximation du flux solaire horizontal global, en fonction notamment d'une valeur d'angle azimutal et d'une valeur de couverture nuageuse, aussi dite nébulosité. A titre d'exemple, la publication Development and validation of a gray box 5 mode) to predict thermal behavior of occupied office buildings, T. Berthou, P. Stabat, R. Salvazet, D. Marchio, Energy and Buildings, pp.91-100, 2014, décrit des modélisations du comportement thermique d'un bâtiment exposé à un rayonnement solaire. Dans le cadre de ces modélisations, un modèle simplifié du flux solaire incident est utilisé (voir figure 3 du 10 document). Le gain du flux solaire dans ce modèle est précisé dans la publication Développement de modèles de bâtiment pour la prévision de charge de climatisation et l'élaboration de stratégies d'optimisation énergétique et d'effacement, T. Berthou, Ecole Nationale Supérieure des Mines de Paris, 2013. Ce modèle prend en entrée des données de 15 nébulosité fournies en temps réel à l'aide de mesures météorologiques, ainsi que des données de géométrie, d'orientation et de masques solaires du bâtiment étudié. Un dispositif de commande du chauffage d'un bâtiment exploitant la solution proposée devrait donc être alimenté en temps réel avec des données de nébulosité et des paramètres 20 caractéristiques du bâtiment. Cette solution est inadaptée à une utilisation grand public. La demande de brevet européen EP 2 505 831 décrit un dispositif de mesure du flux solaire incident sur une éolienne, exploitant les liens entre la valeur du flux solaire incident sur un capteur agencé dans la tête de 25 l'éolienne, la température mesurée en ce capteur, et une valeur de nébulosité. Toutefois, la démarche décrite dans ce document consiste à estimer une valeur de la nébulosité, en vue de détecter une accumulation de glace au niveau des parties mobiles de l'éolienne. Cette démarche nécessite une mesure directe du flux solaire incident via un pyranomètre. 30 Par ailleurs, la demande internationale WO 02/38284, portant sur une centrale de régulation de l'irrigation, décrit, dans le but de calculer un taux d'évapotranspiration, une estimation de flux solaire incident à partir d'une différence de température entre un capteur situé au soleil et un capteur 4 situé à l'ombre. Cependant, ce document ne fournit pas un moyen de construire un modèle dynamique de flux solaire qui soit fonction de cette différence de température, et tenant compte des conditions météorologiques. 5 L'état de l'art ne fournit donc pas de solution satisfaisante pour une estimation indirecte de flux solaire incident, rendant compte de la dynamique infra-horaire de l'évolution du flux solaire, sans nécessiter la fourniture de données météorologiques (sans que de telles données ne doivent être fournies en temps réel, s'il s'agit d'estimer en temps réel le 10 flux solaire), ni l'utilisation de dispositifs peu répandus et onéreux. PRESENTATION GENERALE DE L'INVENTION La présente invention pallie les difficultés relevées ci-avant, par la 15 fourniture, selon un premier aspect, d'un dispositif de détermination indirecte d'un flux solaire incident sur une structure, comprenant des moyens de calcul de flux solaire selon un modèle théorique de flux solaire tenant compte d'une valeur courante de nébulosité, le dispositif comprenant un premier capteur de température extérieure situé à un 20 emplacement exposé au rayonnement solaire, et un deuxième capteur de température extérieure situé à l'ombre, les moyens de calcul recevant en entrée une valeur de biais de température, définie comme la différence entre les températures mesurées par le premier capteur et le deuxième capteur, et étant configurés pour réaliser une estimation de la valeur de 25 nébulosité en fonction de ladite valeur de biais de température. Un dispositif selon l'invention exploite donc un modèle simplifié du flux solaire incident, corrigé avec une valeur courante de nébulosité. Un tel dispositif est capable de réaliser de manière autonome la détermination de cette valeur de nébulosité, en exploitant des mesures de température 30 fournies par un capteur au soleil et un capteur à l'ombre. Cette solution est d'implémentation simple, puisque reposant sur de simples capteurs de température extérieure, et sur un modèle de variation de la nébulosité. Une fois le modèle théorique de flux solaire appris (c'est-à-dire une fois 5 déterminées des valeurs satisfaisantes des paramètres du modèle), ce modèle permet de s'affranchir d'une mesure directe du flux solaire par un pyranomètre. 5 Un dispositif conforme à l'invention peut présenter en outre, de manière optionnelle et non limitative, les caractéristiques suivantes : - Le dispositif comprend en outre une base de données configurée pour construire un modèle dynamique de flux solaire, ledit modèle dynamique étant construit par apprentissage d'un modèle 10 d'évolution temporelle dynamique du biais de température ; - L'apprentissage du modèle dynamique de flux solaire est réalisé conjointement à un apprentissage d'un modèle de comportement de la structure, la valeur courante du flux solaire étant une donnée d'entrée pour la détermination du comportement de la structure, 15 le modèle dynamique de flux solaire et le modèle de comportement de la structure étant construits à l'aide des mesures réalisées par le premier capteur et le deuxième capteur ; - La structure est un bâtiment, et le modèle de comportement est un modèle thermique du bâtiment ; 20 Dans ce dernier cas, le dispositif comprend en outre une sonde de température intérieure du bâtiment ; - Toujours dans ce dernier cas, le dispositif comprend en outre un compteur de consommation électrique ou un compteur de consommation de gaz du bâtiment ; 25 Toujours dans ce dernier cas, le modèle thermique du bâtiment est linéaire, et correspond de préférence à un modèle électrique équivalent de type RxCy ; - La structure est un panneau photovoltaïque, le dispositif de détermination de flux solaire s'interfaçant avec un dispositif de commande 30 du panneau photovoltaïque ; - La structure est une centrale solaire thermodynamique à miroirs, la centrale à miroirs présentant une pluralité de blocs de miroirs, 6 le dispositif de détermination de flux solaire s'interfaçant avec une pluralité de dispositifs de commande d'orientation de miroirs ; - Le dispositif comprend pour chaque bloc de miroirs au moins un capteur de température extérieure situé à un emplacement exposé au 5 soleil. Selon un deuxième aspect, l'invention concerne un procédé de détermination indirecte du flux solaire incident sur une structure, qui comprend les étapes suivantes : 10 mesure de température par un premier capteur de température extérieure situé à un emplacement exposé au rayonnement solaire, et par un deuxième capteur de température extérieure situé à l'ombre, et détermination d'une valeur de biais de température, définie comme la différence entre les deux températures mesurées, 15 estimation, par des moyens de calcul, d'une valeur courante de nébulosité en fonction de ladite valeur de biais de température, et calcul par lesdits moyens de calcul d'une valeur du flux solaire incident sur la structure, d'après un modèle théorique de flux solaire prédéterminé, ledit modèle tenant compte de la valeur courante de nébulosité. 20 Un procédé conforme à l'invention peut présenter en outre, de manière optionnelle et non limitative, les caractéristiques suivantes : - Le flux solaire incident calculé Gcld est exprimé en fonction de la valeur de nébulosité N et d'un flux solaire par temps clair Gclr selon le 25 modèle suivant Gcld = Gclr(1- - aNb), où le flux solaire par temps clair Gcir peut être déterminé en fonction de coordonnées géographiques des capteurs de température, et d'un repérage temporel de la mesure, 30 a et b étant des paramètres à déterminer pour le modèle théorique de flux solaire ; - L'estimation de la valeur de nébulosité N en fonction de la valeur de biais de température AT est alors réalisée selon le modèle suivant : N = «An, avec (1) la fonction suivante : 0(x) = 1 1 + exp(-a(x - xo)) 5 avec xo = (xi_ + x2)/2, et a= 4/(x2 - x1), où xi et x2 sont des valeurs seuil à déterminer pour le modèle de la nébulosité ; - Le procédé comporte une étape préliminaire de détermination de paramètres du modèle théorique du flux solaire à l'aide d'une mesure 10 directe de flux solaire incident dans une zone géographique étendue, les paramètres déterminés étant pris égaux sur toute la zone géographique étendue ; - Le modèle théorique de flux solaire est un modèle dynamique, ledit modèle dynamique étant construit par apprentissage d'un modèle 15 d'évolution temporelle dynamique du biais de température ; - Dans ce dernier cas, le modèle d'évolution temporelle dynamique du biais de température ATt selon le temps t est le modèle d'état linéaire à temps discret comme suit : ATt = 17t + wt, 20 avec wt un bruit blanc gaussien de variance r2 avec r un paramètre à déterminer, et avec 17t = lit-i_ + (q + at)vt, 25 où (q+o-t)vt est un bruit d'innovation, c'est-à-dire un bruit permettant d'introduire de la variabilité dans la dynamique de l'état ?h, q étant un écart-type minimal à déterminer, vt est un bruit blanc gaussien unitaire, 30 et ai est un paramètre à déterminer, ledit modèle d'état linéaire à temps discret du biais de température étant paramétré par un vecteur (q,r,o-i,...,o-T) avec T le nombre de mesures réalisées pour le biais de température, 7 1 8 ledit vecteur étant déterminé par minimisation d'un critère de vraisemblance associé au modèle d'évolution temporelle dynamique du biais de température ; Le procédé est mis en oeuvre de manière simultanée ou séquentielle 5 pour la détermination du flux solaire incident en une pluralité de points d'une zone, avec une estimation de nébulosité pour chacun de ladite pluralité de points, et comprend une étape supplémentaire d'estimation d'une vitesse de déplacement des nuages au-dessus de ladite zone, cette étape supplémentaire prenant en entrée les estimations de nébulosité en 10 chacun des points. PRESENTATION GENERALE DES FIGURES D'autres caractéristiques, buts et avantages de l'invention ressortiront de 15 la description qui suit, qui est purement illustrative et non limitative, accompagnée des dessins annexés suivants : La figure 1 illustre schématiquement les étapes d'un procédé de détermination indirecte de flux solaire incident conforme à l'invention, selon un mode de réalisation. 20 La figure 2 est un schéma d'un dispositif de détermination indirecte de flux solaire incident conforme à l'invention, selon un mode de réalisation particulier pour la détermination du flux solaire incident sur un bâtiment. La figure 3 est un schéma fonctionnel de l'apprentissage conjoint d'un modèle dynamique de flux solaire et d'un modèle de comportement 25 thermique de bâtiment, à l'aide du dispositif illustré en figure 2. La figure 4 est un schéma électrique d'un modèle de circuit électrique R3C2 équivalent, pour la simulation du comportement thermique d'un bâtiment, pouvant être utilisé en relation au procédé illustré en figure 2. 30 DESCRIPTION DETAILLEE DE MODES DE REALISATION Un procédé 10 de détermination indirecte de flux solaire incident conforme à l'invention est illustré en figure 1. 9 Le procédé 10 de la figure 1 est mis en oeuvre par un dispositif placé en extérieur, comprenant un premier capteur de température extérieure situé à un emplacement exposé au rayonnement solaire, et un deuxième capteur de température extérieure situé à l'ombre, ainsi que des moyens 5 de calcul de flux solaire selon un modèle théorique de flux solaire. Les deux capteurs sont situés à proximité l'un de l'autre. Les moyens de calcul sont préférentiellement formés d'un processeur, ou alternativement de tout autre dispositif informatique capable d'exécuter des calculs et d'en mémoriser les résultats. Les moyens de calcul sont 10 capables de communiquer, par une connexion filaire ou une liaison sans fil, avec les capteurs de température. Pour lesdits capteurs, on pourra utiliser de simples sondes de température capables ou non d'afficher sous forme lisible la température détectée. Alternativement, on pourra utiliser tout autre dispositif capable de remplir des fonctions autres qu'un simple 15 relevé de température, comme éventuellement un pyranomètre (notamment durant une phase d'apprentissage du modèle théorique de flux solaire paramétré dans les moyens de calcul). On place l'un des capteurs de température à une localisation toujours ombragée, et l'autre de ces capteurs à une localisation exposée au soleil 20 sur la plus large plage d'ensoleillement possible, qui ne soit pas située dans l'ombre d'un quelconque obstacle. Le flux solaire incident calculé par les moyens de calcul est ici noté Gcld pour « cloud », par opposition à Gdr qui désigne le flux solaire qui serait observé pour un temps parfaitement clair. Le procédé 10 de la figure 1 25 pour la détermination indirecte du flux solaire incident comprend : Une étape 101 de mesure de température par les deux capteurs de température, suivie d'une détermination (par les moyens de calcul, ou bien par les capteurs eux-mêmes s'ils sont aptes à communiquer entre eux) d'une valeur AT de biais de température, définie comme 30 la différence entre les deux températures mesurées, exprimée par exemple en degrés Celsius ; 10 - Une étape 102 d'estimation, par les moyens de calcul, d'une valeur courante de nébulosité N en fonction de la valeur AT de biais de température ; - Une étape 103 de calcul par les moyens de calcul de la valeur Gaia 5 du flux solaire incident sur la structure, d'après le modèle théorique du flux solaire. Le flux solaire calculé pourra être exprimé en Watts par mètre carré. On pourra utiliser tout modèle théorique permettant d'exprimer le flux solaire incident en fonction de la nébulosité. 10 De manière avantageuse, les moyens de calcul peuvent être configurés pour calculer le flux solaire incident selon le modèle de Perrin de Brichambaut et Vauge, le flux solaire étant, dans ce modèle, corrigé par la valeur de nébulosité N correspondant à la couverture nuageuse au 15 niveau des capteurs. Le modèle de Perrin de Brichambaut et Vauge s'exprime comme suit : Gcir = K G0(cos Oz)1'15 OÙ K est un coefficient fixé à 0,81 dans la littérature, OZ est l'angle azimutal 20 (mesuré par exemple par rapport à la direction septentrionale), et Go est la constante solaire extraterrestre, c'est-à-dire la radiation théorique en l'absence d'atmosphère terrestre. Un intérêt de ce modèle est qu'il dépend uniquement du repérage géographique (latitude, longitude) des capteurs, et du repérage temporel 25 de la mesure, pour calculer l'angle azimutal et la radiation théorique. A ce titre, voir la publication Solar engineering of thermal processes, J. A. Duffie et W. A. Beckman, 4th éd., Wiley New York, 2013. Ce modèle est ensuite corrigé par une mesure N E [0,1] de la nébulosité, ou couverture nuageuse, avantageusement selon l'expression suivante : 30 Gcld = Gclr(1 - aNb), où a et b sont des coefficients à déterminer pouvant être fixés, ou bien évoluer au fil des mesures (notamment lors d'un apprentissage du modèle). 11 Une fraction (1-a) du flux solaire par temps clair Gcir est ainsi conservée dans le cas d'un ciel complètement couvert, correspondant à une valeur de nébulosité N = 1. Il convient donc pour utiliser ce modèle de disposer d'une valeur fiable de 5 la nébulosité. Toutefois, cette valeur peut ne pas être disponible de manière simple en temps réel. A l'étape 102 du procédé 10, la valeur courante de la nébulosité N est estimée en fonction du biais de température AT mesuré. Les moyens de 10 calcul peuvent par exemple être configurés pour obtenir une valeur de nébulosité N = (1)(3,T) avec (1) la fonction suivante : 1 1 + exp(-a(x - xo)) 15 avec x0 = (x1 + x2)/2 et a = 4/(x2 - x1). Les valeurs xi et x2 sont ici des valeurs seuil pouvant être fixées, ou bien évoluer au fil des mesures (notamment lors d'un apprentissage du 20 modèle). Ces valeurs sont ici telles que : - Pour un ciel clair, donc une valeur de nébulosité N avoisinant 0, le biais de température soit strictement supérieur à x2, - Pour un ciel très couvert, donc une valeur de nébulosité N avoisinant 1, le biais de température soit strictement inférieur à xi. 25 De manière avantageuse, notamment dans le modèle décrit ci-avant dans le cas où certains paramètres sont destinés à évoluer au fil des mesures, le dispositif mettant en oeuvre le procédé 10 peut comprendre en outre une base de données configurée pour construire le modèle dynamique de 30 flux solaire par apprentissage d'un modèle d'évolution temporelle dynamique du biais de température. 12 Dans toute la description qui va suivre, nous allons décrire trois cas d'application différents du modèle de Perrin de Brichambaut et Vauge précédemment considéré, comprenant un apprentissage d'un modèle théorique de flux solaire, l'apprentissage étant réalisé par des itérations 5 successives d'un procédé de détermination indirecte de flux incident semblable à celui de la figure 1. Les grandeurs auxquelles il est fait référence ci-après, notamment des températures ou des valeurs de flux solaire, peuvent porter un indice «t » pour indiquer leur dépendance selon le temps (grandeur mesurée à l'instant t). Cependant, des références de 10 grandeur dans lesquelles figure ou ne figure pas cet indice désignent bien la même grandeur. Apprentissage conjoint d'un modèle de flux solaire incident sur un bâtiment, et d'un modèle de comportement thermique dudit bâtiment 15 Dans un premier mode de réalisation, la structure est un bâtiment. Des paramètres à déterminer lors de l'apprentissage du modèle de flux solaire sont les paramètres a, b, xi et x2 tels que définis ci-avant. Un objectif est ici de reconstituer le flux solaire incident au niveau du bâtiment lors de périodes temporelles passées. Un autre objectif éventuel est de prédire, à 20 partir du modèle de flux solaire, l'évolution de la température intérieure du bâtiment. Le procédé d'apprentissage qui sera décrit ci-après pourrait, alternativement, être mis en oeuvre avec d'autres modèles théoriques de flux solaire incident, dont des paramètres seraient à déterminer. De manière très avantageuse, le modèle théorique de flux solaire incident 25 est appris conjointement à un modèle de comportement thermique du bâtiment, modélisant de manière dynamique l'évolution d'une température intérieure Ti du bâtiment et d'une température de structure Ts représentative de l'enveloppe du bâtiment. Un choix préférentiel de modèle thermique de bâtiment sera décrit ci-après. Par apprentissage, on 30 entend que des paramètres dudit modèle sont optimisés au fil des mesures de manière automatique, par un algorithme implémenté au sein des moyens de calcul, ou au sein d'un processeur distant pouvant communiquer avec les moyens de calcul par tout moyen. 13 Dans ce cas préférentiel, un modèle d'évolution dynamique des températures Ti et Ts, correspondant au modèle de comportement thermique du bâtiment, interagit avec un modèle d'évolution dynamique du biais de température défini par la différence entre la température Text 5 d'un capteur de température extérieure toujours situé à l'ombre, et la température T - soleil, prise par un second capteur de température extérieure, exposé au soleil. Le modèle dynamique de flux solaire et le modèle de comportement de la structure sont donc construits itérativement au fil des mesures réalisées par le capteur de température 10 Text et le capteur de température Tsoleil. La valeur du flux solaire selon le modèle de flux solaire est une donnée d'entrée pour la détermination du comportement thermique de la structure. Le procédé d'apprentissage du premier mode de réalisation, pour le 15 comportement thermique d'un bâtiment B, comprenant la détermination d'un flux solaire incident G(cid)t sur le bâtiment B en tenant compte de la nébulosité Nt à l'instant t, est mis en oeuvre à l'aide d'un dispositif illustré en figure 2. De manière préférentielle, deux phases d'apprentissage sont à 20 distinguer : une phase d'apprentissage préliminaire au cours de laquelle des valeurs de paramètres du modèle de flux solaire corrigé par la nébulosité sont apprises - dans le cas du modèle de Perrin de Brichambaut et Vauge développé précédemment, les valeurs des paramètres a, b, x1 et x2 - et une phase d'apprentissage d'un modèle thermique de bâtiment 25 exploitant les valeurs ainsi déterminées pour lesdits paramètres. Pour l'une et l'autre de ces phases, le dispositif de la figure 2 peut être utilisé. Un capteur 201 de température, exposé au soleil, fournit la valeur T(soleiI)t en temps réel, et un capteur simple 202 de température, situé à l'ombre, fournit la valeur T(ext)t en temps réel. Un autre capteur 203 situé 30 à l'intérieur du bâtiment donne la valeur Tot intérieure du bâtiment en temps réel. Enfin, une mesure réelle du flux solaire incident à l'instant t, ici noté Qs,t peut être également fournie par le capteur 201, ou alternativement par un capteur additionnel dédié, durant la phase 14 d'apprentissage préliminaire des paramètres du modèle de flux solaire. Par exemple, il est possible d'utiliser comme capteur 201 un capteur simple de température, et d'ajouter un pyranomètre placé à une position exposée au soleil durant la phase d'apprentissage préliminaire (ledit 5 pyranomètre n'étant plus nécessaire durant la deuxième phase d'apprentissage du modèle thermique de bâtiment). La valeur G(c1d)t est fournie selon un procédé de détermination indirecte de flux solaire incident conforme au procédé 10 de la figure 1. Le dispositif mis en place dans le bâtiment B peut, de manière optionnelle, 10 également comprendre un compteur 205 de consommation électrique et un compteur 206 de consommation de gaz du bâtiment. Le procédé d'apprentissage du premier mode de réalisation, prévoyant la construction itérative des modèles dynamiques de flux solaire et de 15 comportement thermique du bâtiment, est illustré schématiquement sur la figure 3. A partir des mesures des températures Tsoleil et Text, une valeur estimée de la nébulosité Nt à l'instant t, puis une valeur de flux solaire noté G(c1d)t à l'instant t (par exemple en Watts par mètre carré) sont fournies selon les 20 modalités décrites ci-avant en relation à la figure 1. De plus, l'évolution temporelle dynamique du biais de température ATt entre les températures Tsoied et Text est modélisée par un modèle d'état linéaire à temps discret, s'exprimant comme suit : 3OEt = rt + Wt, 25 avec wt un bruit blanc gaussien de variance r2 où r est un paramètre à déterminer du modèle d'état linéaire, et =17t-i+ (q + 30 où (q+o-t)vt est un bruit d'innovation, c'est-à-dire un bruit permettant d'introduire de la variabilité dans la dynamique de l'état ?h, q étant un écart-type minimal, vt est un bruit blanc gaussien unitaire, 15 et at est un paramètre à déterminer pour le modèle d'état linéaire du biais de température, parmi les paramètres (ai,...,a-r) où T est le nombre de mesures déjà réalisées pour le biais de température. Avec ce modèle d'évolution dynamique du biais de température ATt entre 5 les températures Tsoleil et Text, la nébulosité est donnée par Nt = (p(Tit ) où la fonction (1) est conforme au modèle exposé ci-avant en relation à l'étape 102 du procédé 10 de la figure 1. Le flux solaire est donné par G - cld,t = G cir,t(1 c11)(1)3b), avec a et b les 10 paramètres définis précédemment pour le modèle théorique de flux solaire incident. Le paramètre q représente l'écart-type minimal du bruit d'innovation (indépendant du temps t), et les valeurs at sont en général nulles, sauf quelques valeurs positives qui permettent de prendre en compte des sauts 15 dans la dynamique de ?h, lorsqu'un changement brusque de couverture nuageuse se produit. Ledit modèle d'évolution dynamique du biais de température ATt est donc paramétré par un vecteur i3 = De manière préférentielle, le vecteur de paramétrage i3 peut être déterminé, préalablement à l'apprentissage du modèle thermique de 20 bâtiment, notamment pendant la phase d'apprentissage préliminaire, par minimisation d'un critère de vraisemblance pénalisé s'exprimant comme suit : :se (À.) = min 1-2 log I, p(ATi, , ATT) + t=1 25 où À est un paramètre permettant d'éviter le sur-apprentissage du modèle et Lig(A7'1, ...,ATT) est la vraisemblance associée au modèle paramétré par le vecteur i3 comme défini ci-avant. La détermination des paramètres du modèle dynamique de biais de 30 température correspond donc à un problème d'optimisation. Pour résoudre ce problème, on pourra transposer des outils standard de 16 l'apprentissage statistique. La vraisemblance Lig se calcule à l'aide d'un filtre de Kalman. A ce titre, on pourra se référer à publication Nonlinear time series: theory, methods and applications with R examples, R. Douc, E. Moulines et D. Stoffer, CRC Press, 2014. Le paramètre À peut, à titre 5 d'exemple, être sélectionné parmi plusieurs valeurs testées en minimisant un critère d'information d'Akaike (AIC). Enfin, pour la recherche du minimum, on pourra utiliser la boîte à outils MATLAB de Mark Schmidt telle que décrite dans la publication Fast optimization methods for Ll regularization: A comparative study and two new approaches, M. Schmidt, 10 G. Fung et R. Rosales, European Conference on Machine Learning (ECML), 2007, ou toute autre méthode appropriée connue. De manière alternative, la phase d'apprentissage préliminaire des paramètres a, b, xi et x2 du modèle théorique du flux solaire est mise en oeuvre par une mesure directe du flux solaire incident dans une zone 15 géographique étendue (à l'aide d'un ou de plusieurs pyranomètres). L'homme de l'art comprendra aisément comment, avec le modèle de Perrin de Brichambaut et Vauge, une pluralité de mesures directes peut permettre de déterminer des valeurs de ces paramètres. Dans ce dernier cas, les valeurs des paramètres a, b, xi, x2 et i3 sont prises uniformes sur 20 toute la zone géographique étendue. Les paramètres a, b, xi, x2 et 6 déterminés pour le modèle de flux théorique incident à l'issue de l'apprentissage peuvent être étendus à d'autres bâtiments, situés préférentiellement en une même zone géographique. Dans ce cas, pour lesdits autres bâtiments, la phase d'apprentissage préliminaire du modèle 25 théorique de flux solaire peut ne pas être mise en oeuvre, et le comportement thermique du bâtiment peut être appris directement. Revenant au procédé d'apprentissage de la figure 3, l'apprentissage s'effectue par le suivi de l'évolution d'un vecteur d'état X = (TI, Ts, où on rappelle que Ti est la température intérieure du bâtiment, Ts est une 30 température de structure représentative de l'état thermique de l'enveloppe du bâtiment (pour laquelle il n'est pas nécessairement possible d'avoir accès à une mesure réelle), et correspond à la dynamique du biais de température extérieure (dépendant de la 17 nébulosité). Plus précisément, ici = n, où ri est la même grandeur décrite précédemment pour le modèle d'état linéaire à temps discret développé pour l'évolution temporelle du biais de température 3,Tt entre les températures Text et Tsoied. La valeur Xt à l'instant t de ce vecteur d'état est 5 estimée à l'aide : - Du modèle théorique du flux solaire et du modèle d'évolution dynamique du biais de température décrits ci-avant, s'agissant de l'état - D'un modèle de comportement thermique du bâtiment, s'agissant 10 des températures Ti et Ts, ce dernier modèle exploitant préférentiellement une valeur de flux solaire estimée à partir de l'état fit, via le modèle théorique de flux solaire. On pourra très avantageusement utiliser, pour le modèle de 15 comportement thermique du bâtiment, un modèle électrique équivalent RxCy. Il est connu de modéliser de manière dynamique le comportement thermique d'un bâtiment à l'aide d'un modèle de circuit électrique équivalent du type réseau RxCy, dans lequel x et y sont supérieurs ou égaux à 2, à l'aide notamment de données de fonctionnement en boucle 20 ouverte de la température extérieure, ainsi que de données relevées en boucle fermée mais exploitées en boucle ouverte. Ces modèles peuvent notamment servir à déterminer une puissance de chauffage nécessaire pour chauffer l'intérieur du bâtiment. Nous reprenons en figure 4 un modèle électrique équivalent R3C2. Il est 25 rappelé que d'après l'analogie électrique standard, les températures sont associées à des potentiels du circuit, et les puissances aux courants électriques. Le modèle de la figure 4 comprend trois résistances Rf, Ri et Ro et deux condensateurs Cs et Cres. Les deux résistances Ri et Ro sont en série sur 30 une branche placée en parallèle de la résistance Rf. Le condensateur Cs relie à la masse le point correspondant à la température Ts de l'enveloppe du bâtiment. Le condensateur Cres relie le point commun aux résistances Rf et Ri, qui correspond à la température intérieure Ti du bâtiment, à la 18 masse. Enfin, un générateur de tension relie le point commun des résistances Rf et Ro, correspondant ici à la température extérieure Text mesurée à l'ombre, à la masse. Un générateur de courant Qs est placé en parallèle du condensateur Cs, et un générateur de courant Qres est placé 5 en parallèle du condensateur Cres. Les équations différentielles couplées régissant l'évolution des puissances Qs et Qres peuvent être déterminées sans difficulté par l'homme du métier, pour compléter le modèle de comportement thermique du bâtiment. On pourra alternativement considérer un autre modèle de comportement 10 thermique linéaire, notamment un modèle électrique équivalent de type RxCy, pour le bâtiment B. De manière plus générale, toute représentation d'état linéaire discrète pour les températures intérieure et extérieure du bâtiment pourrait être utilisée. Il pourrait également être envisagé de réaliser un modèle de comportement thermique pour un ensemble de 15 bâtiments, dont un modèle électrique équivalent serait constitué de plusieurs circuits électriques équivalents reliés entre eux. En conformité avec un tel modèle R3C2, l'évolution dynamique de la température intérieure Ti et de la température de structure Ts, selon le modèle dynamique de comportement thermique du bâtiment construit 20 lors de l'apprentissage, est modélisée comme suit. On note Qr la somme de la puissance de chauffage et des apports thermiques internes au bâtiment. Le vecteur U = (Text Qr Gclr)1 désigne les sollicitations mesurées du modèle joint du flux solaire et du comportement thermique du bâtiment, et le vecteur Y = (Ti /Text)' 25 désigne les sorties du système. Parmi les grandeurs apparaissant dans le vecteur U, les grandeurs Text et Qr sont par exemple mesurées directement, tandis que Gdr est fourni par exemple à l'aide du modèle de flux solaire développé précédemment.Le vecteur Y est également mesuré (capteurs 201, 202, 203). Les mesures obtenues pour le vecteur Y 30 peuvent être bruitées. Le vecteur Y mesuré est noté Y*. Le vecteur Y correspond donc à une observation en temps réel. Le vecteur d'états cachés (non-mesurés) X = (Ti n ij)', où nt correspond à fit, est divisé en deux sous-vecteurs t = ?h et at = (Ti Ta)' . D'un point de vue statistique, 19 si le modèle global du vecteur d'états X est non-linéaire, du fait de la non-linéarité du modèle de flux solaire, le modèle de bâtiment (at) est linéaire conditionnellement à l'état fit. On a en effet les égalités suivantes : = + vt, at = Ut, 0) + Aa(e)at-i+ Et, Yt = [O + Cat + Et, 5 où il est rappelé que vt est un bruit blanc gaussien unitaire, la fonction fa et les matrices Aa et C se déduisent aisément des équations du circuit R3C2, Et et Et sont des bruits blancs gaussiens, et e correspond au vecteur 0_(1 iy Cres Cs Rf Ro Ri) 10 Revenant à nouveau à la figure 3, pour un instant t de mesure donné, une fois obtenue la valeur estimée du flux solaire, via l'estimation de l'état on peut donc obtenir, suivant l'égalité ci-dessus pour le sous-vecteur at, des valeurs estimées de la température intérieure Tot et la température extérieure Tst. Une valeur ainsi calculée du vecteur d'état Xt = (Tit, Tst, 15 peut ainsi être obtenue. On peut alors en particulier comparer les valeurs estimées de Tot ett dans le vecteur calculé Yt avec les valeurs mesurées dans le vecteur observé Y*t. Si la différence entre la valeur mesurée Y*t et la valeur estimée pour le vecteur observé Yt est déterminée comme suffisamment faible, par 20 exemple si une distance entre ces deux vecteurs est inférieure à un seuil prédéterminé p, l'apprentissage conjoint du modèle théorique de flux solaire et du modèle de comportement thermique du bâtiment est considéré comme terminé. On sort alors de la boucle représentée de manière schématique sur la figure 3 pour la phase d'apprentissage des 25 modèles. Si, à l'inverse, la valeur mesurée et la valeur estimée sont déterminées comme éloignées, on reste dans ladite boucle de la figure 3, et la trajectoire dynamique du flux solaire est ajustée. Par trajectoire dynamique du flux solaire, on entend la dynamique d'évolution du vecteur 30 d'états cachésx = (Ti T. i')' où il est rappelé que Ti correspond à la dynamique d'évolution du biais de température, dont dépend la valeur 20 estimée du flux solaire selon le modèle développé précédemment. Une nouvelle itération pour le calcul d'une valeur de flux solaire est initiée, et les paramètres du modèle thermique de bâtiment sont mis à jour selon un algorithme d'estimation de modèle, implémenté dans les moyens de calcul 5 ou dans un processeur pouvant communiquer avec les moyens de calcul. On applique préférentiellement l'algorithme EM (Espérance-Maximisation), bien connu dans le domaine technique de l'apprentissage statistique, au cas des modèles d'état gaussiens conditionnellement linéaires. L'estimateur des paramètres 0 est donc un estimateur de 10 maximum de vraisemblance. On pourra réaliser, pour calculer le maximum de vraisemblance, un lissage particulaire de l'état caché non-linéaire, avec prise en compte de la sous-structure linéaire pour réduire la variance de l'estimation (méthode appelée Rao-Blackwellized Particle Smoother Expectation Maximization). Cette méthode est décrite par 15 exemple dans la publication Identification of mixed linear/nonlinear statespace models, F. Lindsten et T. B. Schôn, 49th IEEE Conference on Decision and Control, 2010. L'absence de flux solaire lors des mesures réalisées de nuit implique que le comportement des états Ti et Ts lors de la nuit est exclusivement lié aux 20 sollicitations autres que le flux solaire, qui sont bien mesurées, ce qui permet de renforcer la capacité d'identification conjointe de ce modèle de bâtiment et du modèle théorique de flux solaire. Il est à noter qu'une estimation de flux solaire incident selon le modèle 25 appris ne nécessite pas l'usage d'un pyranomètre. Comme indiqué précédemment, les paramètres a, b, xi, x2 et i3 déterminés pour le modèle de flux théorique incident à l'issue de l'apprentissage peuvent être étendus à d'autres bâtiments, situés préférentiellement en une même zone géographique. Cela permet de limiter l'usage du pyranomètre. Les 30 paramètres déterminés à l'issue de l'apprentissage pour le modèle thermique de bâtiment peuvent également être élargis à d'autres bâtiments similaires. 21 Estimation d'un flux solaire incident futur Dans un deuxième mode de réalisation, on s'intéresse à la prédiction du flux solaire incident sur des périodes temporelles futures, pour une unité 5 de production électrique exploitant la puissance du rayonnement solaire, typiquement une centrale de production photovoltaïque. En effet, moyennant la prise en compte de l'orientation du flux solaire incident, le rendement d'une telle installation est fortement corrélé au flux solaire incident, qui, comme vu ci-avant, est lui-même corrélé à la nébulosité. 10 Un intérêt d'utiliser un modèle théorique de flux solaire incident, ne prenant pas en compte l'état de fonctionnement de la centrale de production - défauts sur des panneaux photovoltaïques, vieillissement de l'installation, salissures... - est de permettre une comparaison de la puissance théorique avec la puissance effective observée. 15 Une telle comparaison permet de remonter à des données d'état de fonctionnement de la centrale de production Avantageusement, on pourra répartir plusieurs capteurs de température (et ainsi obtenir plusieurs mesures de biais de température) à plusieurs localisations, notamment pour un site de production vaste. Les données 20 de nébulosité ainsi estimées en temps réel en plusieurs points du site de production peuvent permettre d'extrapoler une estimation d'une direction et d'une vitesse de déplacement des nuages. Un système selon cette dernière variante permet, outre la possibilité de prédire de manière fiable la puissance produite globalement par l'unité de 25 production sans utiliser de pyranomètre, de prédire à court terme des évolutions de la production globale ou en des zones précises de l'unité de production, en fonction du déplacement estimé des nuages. Détermination instantanée d'un flux solaire incident 30 Dans un troisième mode de réalisation, un objectif est d'estimer en temps réel le flux solaire incident en plusieurs zones d'une centrale solaire thermodynamique à concentration, ou centrale solaire à miroirs. Dans une telle installation, le rayonnement solaire est focalisé sur une surface de 22 chauffe. Il est crucial de contrôler finement la puissance de rayonnement concentrée sur la surface de chauffe : cette puissance doit être suffisante pour entraîner la montée en température souhaitée, mais maîtrisée afin de ne pas endommager ladite surface de chauffe ou les miroirs. 5 Dans ce contexte, des variations rapides de la couverture nuageuse au- dessus de l'installation posent de grandes difficultés pour la réorientation adéquate des miroirs. La variation de nébulosité peut être trop rapide pour permettre la ré-orientation en temps utile des miroirs. Il est courant, face à ces difficultés, de mettre hors service l'installation toute entière dans 10 des périodes de grande variation de couverture nuageuse, en attendant un retour au plein soleil. Le procédé de détermination indirecte de flux solaire décrit ci-avant, ainsi que le procédé d'apprentissage des paramètres du modèle théorique de flux solaire, peuvent très avantageusement être utilisés pour obtenir une 15 cartographie précise du flux solaire différenciée sur plusieurs blocs de miroirs de la centrale solaire à concentration, et/ou une cartographie de la nébulosité différenciée sur plusieurs blocs de miroirs. On pourra notamment utiliser à cette fin un dispositif conforme à celui décrit précédemment en relation au procédé 10 de la figure 1, ou une 20 pluralité de tels dispositifs, comprenant un capteur de température situé au soleil et un capteur situé à l'ombre ainsi que des moyens de calcul, interfacés avec une pluralité de dispositifs de commande d'orientation de miroirs (typiquement un dispositif de commande pour chaque bloc de miroirs). 25 De manière avantageuse, au moins un capteur de température extérieure pour chaque bloc de miroirs est disposé à un emplacement exposé au soleil. De manière optionnelle, le déplacement des nuages au-dessus de la centrale thermodynamique à miroirs peut être anticipé, ce qui octroie un 30 temps d'avance pour réorienter les miroirs de façon adéquate et contrôler la puissance concentrée sur la surface de chauffe. Un dispositif de pilotage de l'orientation de chaque miroir ou bloc de miroirs peut être contrôlé à 23 l'aide de la valeur acquise correspondante de flux solaire, ou la valeur acquise correspondante de nébulosité. Une implémentation possible pour mesurer le déplacement des nuages consiste à mesurer de manière simultanée ou séquentielle le flux solaire 5 incident en une pluralité de points d'une zone, avec une estimation de nébulosité pour chacun de ladite pluralité de points, selon un procédé conforme au procédé 10 de la figure 1, et en ajoutant une étape supplémentaire d'estimation d'une vitesse de déplacement des nuages au-dessus de ladite zone, cette étape supplémentaire prenant en entrée les 10 estimations de nébulosité en chacun des points. Les trois modes de réalisation présentés ci-avant ne sont pas restrictifs quant à l'utilisation qui est faite du flux solaire estimé : un flux solaire incident pourra être estimé selon le procédé de détermination indirecte 15 exposé ci-avant, en temps réel ou pour des périodes futures, ou bien reconstitué pour des périodes passées, dans le contexte de l'un quelconque de ces trois modes de réalisation. Par ailleurs, un autre modèle de flux solaire incident que le modèle de Perrin de Brichambaut et Vauge peut être exploité. On citera par exemple 20 le modèle de Kasten, qui dépend d'un facteur de trouble de Linke, pour lequel l'homme du métier saura transposer le procédé d'apprentissage décrit ci-avant. Ce modèle est décrit dans la publication Selecting a clearsky mode/ to accurately map solar radiation from satellite images, C. Rigollier et L. Wald, EARSeL Symposium 1999 "Remote sensing in the 21st 25 Century: economic and environmental applications", Valladolid, 2000. Le modèle d'évolution dynamique du biais de température exposé ci-avant n'est pas non plus limitatif. On pourra notamment considérer le modèle suivant : ATt = + Wt, 30 avec wt un bruit blanc gaussien de variance r2 où r est un paramètre à déterminer du modèle, et = + µt + qvt, 24 où q et vt sont pris identiques au modèle précédemment décrit, et lit est nul en moyenne, sauf pour des périodes très brèves pour tenir compte d'éventuels REVENDICATIONS1. Dispositif de détermination indirecte du flux solaire incident sur une structure, le dispositif comprenant des moyens de calcul (210) de flux solaire selon un modèle théorique de flux solaire tenant compte d'une valeur courante de nébulosité, le dispositif étant caractérisé en ce qu'il comprend un premier 10 capteur de température extérieure (201) situé à un emplacement exposé au rayonnement solaire, et un deuxième capteur de température extérieure (202) situé à l'ombre, les moyens de calcul (210) recevant en entrée une valeur de biais de température, définie comme la différence entre les températures 15 mesurées par le premier capteur et le deuxième capteur, et les moyens de calcul étant configurés pour réaliser une estimation de la valeur de nébulosité en fonction de ladite valeur de biais de température. 20 2. Dispositif de détermination indirecte de flux solaire selon la revendication 1, comprenant en outre une base de données configurée pour construire un modèle dynamique de flux solaire, ledit modèle dynamique étant construit par apprentissage d'un modèle d'évolution temporelle dynamique du biais de température. 25 3. Dispositif de détermination indirecte de flux solaire selon l'une des revendications 1 ou 2, dans lequel l'apprentissage du modèle dynamique de flux solaire est réalisé conjointement à un apprentissage d'un modèle de comportement de la structure, 30 la valeur courante du flux solaire étant une donnée d'entrée pour la détermination du comportement de la structure, 26 le modèle dynamique de flux solaire et le modèle de comportement de la structure étant construits à l'aide des mesures réalisées par le premier capteur et le deuxième capteur. 5 4. Dispositif de détermination indirecte de flux solaire selon la revendication 3, dans lequel la structure est un bâtiment, et le modèle de comportement est un modèle thermique du bâtiment, le dispositif comprenant en outre une sonde de température 10 intérieure du bâtiment (203), et de préférence un compteur de consommation électrique (205) et/ou un compteur de consommation de gaz (206) du bâtiment. 5. Dispositif de détermination indirecte de flux solaire 15 selon la revendication 4, dans lequel le modèle thermique du bâtiment est linéaire, et correspond de préférence à un modèle électrique équivalent de type RxCy. 6. Dispositif de détermination indirecte de flux solaire 20 selon l'une des revendications 1 à 3, dans lequel la structure est un panneau photovoltaïque, le dispositif de détermination de flux solaire s'interfaçant avec un dispositif de commande du panneau photovoltaïque. 25 7. Dispositif de détermination indirecte de flux solaire selon l'une des revendications 1 à 3, dans lequel la structure est une centrale solaire thermodynamique à miroirs, la centrale à miroirs présentant une pluralité de blocs de miroirs, et le dispositif de détermination de flux solaire s'interfaçant avec 30 une pluralité de dispositifs de commande d'orientation de miroirs. 27 8. Dispositif selon la revendication 7, comprenant pour chaque bloc de miroirs au moins un capteur de température extérieure situé à un emplacement exposé au soleil. 5 9. Procédé de détermination indirecte du flux solaire incident sur une structure, le procédé comprenant les étapes suivantes : mesure de température (101) par un premier capteur de température extérieure situé à un emplacement exposé au rayonnement solaire, et par un deuxième capteur de température extérieure situé à l'ombre, et détermination d'une valeur de biais de température, définie comme la différence entre les températures mesurées par le premier capteur et le deuxième capteur, estimation (102), par des moyens de calcul, d'une valeur courante de nébulosité en fonction de ladite valeur de biais de température, et calcul (103) par lesdits moyens de calcul d'une valeur du flux solaire incident sur la structure, d'après un modèle théorique de flux solaire prédéterminé, ledit modèle tenant compte de la valeur courante de nébulosité. 10. Procédé de détermination indirecte de flux solaire selon la revendication 9, dans lequel le flux solaire incident calculé Gcld est exprimé en fonction de la valeur de nébulosité N et d'un flux solaire par temps clair Gd, selon le modèle suivant : Gad - Gcir(1 - aNb), où le flux solaire par temps clair Gd, peut être déterminé en fonction de coordonnées géographiques des capteurs de température, et 30 d'un repérage temporel de la mesure, 28 a et b étant des paramètres à déterminer pour le modèle théorique de flux solaire. 11. Procédé de détermination indirecte de flux solaire selon 5 l'une des revendications 9 ou 10, dans lequel l'estimation de la valeur de nébulosité N en fonction de la valeur de biais de température AT est réalisée selon le modèle suivant : N = (p(AT), 10 avec (1) la fonction suivante : 1 1 + exp(-a(x - xo)) 15 avec xo = (x1 + x2)/2 et a = 4/(x2 - x1) où xi et x2 sont des valeurs seuil à déterminer pour le modèle de la 20 nébulosité. 12. Procédé de détermination indirecte de flux solaire selon l'une des revendications 9 à 11, comportant une étape préliminaire de détermination de paramètres du modèle théorique du flux solaire 25 à l'aide d'une mesure directe de flux solaire incident dans une zone géographique étendue, les paramètres déterminés étant pris égaux sur toute la zone géographique étendue. 13. Procédé de détermination indirecte de flux solaire selon l'une des revendications 9 à 12, dans lequel le modèle théorique de flux solaire est un modèle dynamique, 29 ledit modèle dynamique étant construit par apprentissage d'un modèle d'évolution temporelle dynamique du biais de température. 14. Procédé de détermination indirecte de flux solaire selon 5 la revendication 13, dans lequel le modèle d'évolution temporelle dynamique du biais de température ATt selon le temps t est le modèle d'état linéaire à temps discret comme suit : ATt = llt +mit, 10 avec wt un bruit blanc gaussien de variance r2 avec r un paramètre à déterminer, et avec nt = nt-1 + (q + cit)vt, 15 où q est un écart-type minimal à déterminer d'un bruit d'innovation, c'est-à-dire un bruit permettant d'introduire de la variabilité dans la dynamique de l'état nt, q étant un écart-type minimal à déterminer, vt est un bruit blanc gaussien unitaire, et o-t est un paramètre à déterminer, 20 ledit modèle d'état linéaire à temps discret du biais de température étant paramétré par un vecteur (q,r,o-i,...,o-T) avec T le nombre de mesures réalisées pour le biais de température, ledit vecteur étant déterminé par minimisation d'un critère de vraisemblance associé au modèle d'évolution temporelle dynamique 25 du biais de température. 15. Procédé de détermination indirecte de flux solaire selon l'une des revendications 9 à 14, le procédé étant mis en oeuvre de manière simultanée ou séquentielle pour la détermination du flux 30 solaire incident en une pluralité de points d'une zone, avec une estimation de nébulosité pour chacun de ladite pluralité de points, 30 le procédé comprenant une étape supplémentaire d'estimation d'une vitesse de déplacement des nuages au-dessus de ladite zone, cette étape supplémentaire prenant en entrée les estimations de nébulosité en chacun des points. 5