Brevet FR3114427A1 — publié le 20220325

Procédé de prévision de l’impact d’une tempête sur une infrastructure de réseau électrique d’un territoire comprenant un nombre N dep de départs haute tension aérien dits départs HT aériens d qui comprend pour une tempête à venir: - pour chaque départ HT aérien dudit territoire, au moins un calcul de probabilité (100) que ledit départ soit coupé selon l’équation tirée d’un modèle de Cox avec des covariables Xi d’environnement comportant des données prévisionnelles de ladite tempête selon un maillage dudit territoire et des données propres à l’implantation desdites lignes HT aériennes et desdits départs HT aériens ; - au moins une étape de calcul d’un nombre de coupures sur le territoire ou une partie du territoire ; - au moins un affichage et/ou une sauvegarde (120) des données du nombre de coupures en fonction des départs potentiellement impactés et de la probabilité de coupure pour chaque départ. Figure 4 PROCEDE DE PREDICTION DE RISQUES METEOROLOGIQUES SUR UNE INFRASTRUCTURE DE RESEAU ELECTRIQUE AERIEN L’invention relève du domaine de la gestion de crise des réseaux électriques aériens et en particulier de la prévision des risques météorologiques de tempêtes pouvant affecter des réseaux électriques aériens haute tension. Pour les opérateurs de distribution d’électricité, la réalimentation des clients en cas d’incidents sur les réseaux électriques aériens est un enjeu majeur, y compris en période de crise climatique, malgré les difficultés inhérentes à de telles situations. Afin de mieux réagir et d’améliorer la réalimentation des clients en cas d’aléas climatiques tels que des tempêtes, la demanderesse a engagé différents travaux permettant de mieux appréhender les impacts sur le réseau de distribution des phénomènes de tempête. Ces modélisations anticipées ont pour objet de permettre de mieux préparer et pré-mobiliser les ressources pour gérer les incidents lorsque de tels phénomènes sont annoncés et de réduire leurs conséquences sur la continuité de la fourniture d’électricité. Pour ce faire, une prévision des coupures électriques sur le réseau électrique causées par le passage d’un épisode venteux d’ampleur significative est à considérer. Le problème métier posé, est celui de prévoir, en amont des tempêtes, le nombre de coupures électriques susceptibles d’impacter le réseau électrique. Cela a pour objectif d’anticiper et optimiser la mobilisation des ressources et d’accélérer la reprise de l’alimentation électrique pour les clients coupés. Lors d’un évènement de crise, le nombre de clients coupés en simultané peut être très élevé, ce qui nécessite une prise en charge complexe, c’est pourquoi il est crucial de pouvoir l’anticiper. Pour répondre à ce besoin, qui se pose pour tous les gestionnaires de réseaux, des travaux basés sur de l’apprentissage profond sont connus notamment des documents : Tervo, Roope & Karjalainen, Joonas & Jung, Alexander (2019) “Short-Term Prediction of Electricity Outages Caused by Convective Storms”. IEEE Transactions on Geoscience and Remote Sensing. PP. 1-9. 10.1109/TGRS.2019.2921809. D. Cerrai et al., "Predicting Storm Outages Through New Representations of Weather and Vegetation" in IEEE Access, vol. 7, pp. 29639-29654, 2019, doi: 10.1109/ACCESS.2019.2902558. Problème technique Les méthodes de classification par apprentissage profond calculent un score par élément à classer. Dans le cadre de la présente invention qui concerne l’estimation des départs haute tension moyenne tension (HTA) aériens qui risquent d’être coupés, le problème revient à classer les départs haute tension aériens selon deux classes : « coupé » ou « non coupé ». Une difficulté inhérente à ces méthodes est la nécessité de définir un seuil pour effectuer le classement. Ce seuil peut être évalué à partir de différents tests mais son choix reste arbitraire. Le choix d’un score inapproprié rendra l’évaluation du nombre de coupure peu précise. Par ailleurs, les méthodes de classification ne permettent pas d’évaluer des intervalles de confiance du nombre de départs prévus coupés par l’étape de classification. Ce problème est important car sans intervalle de confiance, il n’est pas possible d’apprécier le risque pris en prenant une décision par rapport aux prévisions. Par ailleurs, sur le terrain, les moyens de prévisions à disposition des équipes d’exploitation et de gestion de crise sont des fichiers indiquant un nombre de coupures indicatif en fonction d’un temps de retour de la tempête c’est à dire l’intervalle de temps moyen séparant 2 réalisations de tempête de même intensité, ces fichiers nécessitant d’être réinterprété par un œil expert et une bonne connaissance du terrain. A la disposition des métiers, il n’y a donc pas encore de solution quantifiant l’impact des tempêtes en terme de nombre de coupures à l’échelle de tout un territoire. Le modèle de prévision de panne de la présente invention, permet de solutionner ce problème métier avec un niveau de fiabilité et une précision géographique accrus. Au vu de l’art antérieur, la présente demande propose un procédé de prévision de risques météorologiques sur une infrastructure de réseau électrique d’un territoire comprenant un nombre N dep de départs haute tension aériens d qui comprend pour une tempête à venir: - pour chaque départ HT aérien dudit territoire, au moins un calcul de probabilité que ledit départ soit coupé selon l’équation tirée du modèle de Cox comportant comme contrainte la valeur de rafale maximale à subir par ledit départ entre le début de la tempête et un instant t : avec des covariables d’environnement desdits départs HT aériens et βi les paramètres représentant l’impact desdites covariables sur le risque de défaillance desdits départs HT aériens, lesdites covariables d’environnement comportant des données prévisionnelles de ladite tempête selon un maillage dudit territoire et des données propres à l’implantation desdites lignes HT aériennes et desdits départs; - au moins une étape de calcul d’un nombre moyen de coupures entre le début de la tempête et l’instant t sur le territoire ou une partie du territoire comportant une étape de sélection d’un nombre de départs de la zone géographique considérée selon la formule : pour laquelle est la variable aléatoire représentant le nombre de coupures desdits départs HT aériens entre le début de la tempête et ledit instant t ; - au moins un affichage et/ou une sauvegarde des données du nombre moyen de coupures en fonction des départs potentiellement impactés et de la probabilité de coupure pour chaque départ ; le calcul de probabilité étant réitéré à chaque instant t pendant la tempête prévue afin d’évaluer l’évolution de pendant la tempête, la mémoriser et/ou l’afficher. Ce procédé dont l’exécution est peu consommatrice de ressources de calcul donne une bonne approximation des risques de coupure de départs selon les prévisions de tempête sur le territoire malgré la complexité et le volume de données à traiter. Le procédé peut comporter: - au moins un calcul d’un intervalle de confiance ic du nombre de départs coupés entre le début de la tempête et un instant t sur le territoire ou une partie du territoire au moyen d’une approximation de la loi de par une loi de Poisson de paramètre - au moins un calcul des quantiles de cette loi pour obtenir ledit intervalle de confiance ic , - au moins une mémorisation dudit intervalle de confiance et/ou l’affichage dudit intervalle de confiance ; ledit calcul dudit intervalle de confiance ic étant réitéré pour chaque valeur de t comprise entre 0 et l’instant de fin de la tempête afin d’évaluer l’évolution dudit intervalle de confiance pendant la tempête, l’afficher et/ou le mémoriser. La réitération permet d’affiner les calculs en fonction de l’évolution des prévisions et le calcul et l’affichage de l’intervalle de confiance permet aux opérateurs d’évaluer la précision des prévisions pour mieux préparer l’affectation des ressources d’intervention. A partir de la probabilité de coupure , le procédé peut comporter : - une ou plusieurs étapes de calcul de nombres d’incidents entre le début de la tempête et un instant t sur la ligne HT aérienne associée à un départ HT aérien pour chaque départ HT aérien d au moyen d’une loi de Poisson de paramètre : - une ou plusieurs étapes de calcul de la loi du nombre d’incidents HT aérienne sur ledit territoire contenant les départs HT aériens au moyen d’une loi de Poisson de paramètre : - une ou plusieurs étapes de sauvegarde et/ou d’affichage de la moyenne E[N inc ] desdits nombres d’incidents estimés pour une utilisation en tant qu’outil de préparation de moyens d’intervention. Pour chaque départ HT aérien, les covariables propres à l’implantation peuvent comporter au moins les données suivantes: une zone où se situe le départ HT aérien; le type de boisement observé sur au moins une portion de ladite zone; la longueur de la ligne HT aérienne reliée au départ; la longueur de la ligne HT aérienne soumise à un risque de boisement. Ces covariables offrent une précision adaptée pour les prévisions sans complexifier les calculs sur un territoire tel que la France. Les covariables météo peuvent comporter au moins les données suivantes: la saison d’occurrence de la tempête; la température moyenne pendant la tempête issue des prévisions de tempête; l’humidité moyenne pendant la tempête issue des prévisions de tempête; la direction du vent pendant la tempête (Nord, Sud, Est ou Ouest). Là aussi les covariables choisies correspondent à une précision suffisante sans ajouter trop de données à considérer. Les données prévisionnelles comportant la direction du vent, cette dernière peut être calculée en prenant la moyenne des directions des vents lorsque les rafales dépassent le quantile à 75% des rafales prévues pendant la journée durant laquelle la tempête est prévue. Ceci permet de d’obtenir une évaluation plus précise de la direction du vent pendant l’épisode de tempête et d’éviter de prendre en compte des périodes de la journée sans tempête. Les calculs peuvent être réitérés à chaque modification des données prévisionnelles de tempêtes sur ledit territoire pour affiner la prévision. Selon un second aspect, l’invention prévoit une méthode d’agrégation de modèles qui comprend pour chaque départ HT aérien d’une zone considérée une étape de calcul d’appartenance à l’une des deux classes « coupé », « non coupé », une étape d’ordonnancement des départs et d’affectation à la classe « coupé » aux départs de rang le plus élevé selon ledit ordonnancement et une étape d’affichage et/ou de mémorisation desdits départs appartenant à la classe « coupé ». étant le nombre moyen de départ prévus coupés sur une zone considérée et sur toute la période de la tempête. De la sorte, en utilisant les résultats du modèle de Cox comme données d’entrée, un calcul de seuil délicat est évité. Les résultats peuvent être remontés à une maille de postes sources auxquels appartiennent les départs HT aériens, l’algorithme comportant une localisation des postes sources dont au moins un départ HT aérien est classé coupé et un affichage d’une prévision de postes sources à départs HT aériens classés coupés. Les départs HT aériens étant regroupés selon des postes source, le procédé peut comporter selon une méthode d’agrégation de modèles: - au moins une étape de calcul du nombre moyen de postes sources avec au moins un départ HT aérien coupé sur une population de postes sources d’au moins une partie du territoire: étant le nombre de départs HT aériens pour le poste source j, et, étant la probabilité de coupure du départ HT aérien pendant la tempête à venir, et, - au moins une étape de calcul d’appartenance des postes source à l’une des deux classes « au moins un départ HT aérien coupé » et « aucun départ HT aérien coupé » produisant un ordonnancement des postes source, l’appartenance à la classe « au moins un départ HT aérien coupé » étant affectées aux postes source de rang le plus élevé selon ledit ordonnancement et ; - au moins une étape d’affichage et/ou de mémorisation (360) desdits postes sources appartenant à la classe « au moins un départ HT aérien coupé » . Ceci simplifie les calculs sans trop impacter les prévisions puisque les interventions peuvent être regroupées par poste source. La méthode d’agrégation de modèles peut avantageusement utiliser un algorithme connu sous le nom de XGBoost dans le langage R. La méthode peut prendre en compte en tant que variables le ou les nombres d’heures pour lequel les départs HT aériens sont affectés par une ou plusieurs rafales supérieures à une ou plusieurs vitesses données. La méthode notamment prendre en compte en tant que variables au moins l’une des données: nombre de coupures passées liées au vent; la longitude et la latitude du départ HT aérien en remplacement de la région où se situe le départ HT aérien ; l’altitude ; la valeur en degrés de la direction du vent. Afin de réaliser l’apprentissage pour déterminer les paramètres des algorithmes, le procédé peut comprendre : - une ou plusieurs étapes préalable de collecte de données d’environnement des lignes électriques hautes tension aériennes dites lignes HT aériennes reliées aux dits départs HT aériens sur ledit territoire et l’organisation d’une base de données associant lesdites données d’environnement desdits départs HT aériens sous formes de covariables selon un maillage dudit territoire et, - une ou plusieurs étapes d’apprentissage, à partir d’une base de données regroupant des données de prévisions météo d’une pluralité de tempêtes passées et les pannes réseau associées, de calcul de paramètres βi d’un modèle de Cox dont le taux de défaillance s’exprime de la façon suivante : ; les paramètres représentant l’impact desdites covariables sur le risque de défaillance en considérant ledit modèle de COX comme un modèle contrainte/résistance dont la contrainte est la rafale maximale subie par un départ avant une coupure ou au cours desdites tempêtes passées si le départ n’est pas coupé et pour lequel est le taux de défaillance lorsque toutes les covariables sont nulles. D’autres caractéristiques, détails et avantages de l’invention apparaîtront à la lecture de la description détaillée ci-après d’exemples de réalisation non limitatifs, et à l’analyse des dessins annexés, sur lesquels : montre la courbe de probabilité de survie d’un départ HT aérien en fonction des rafales maximales subies depuis le début de la tempête ; montre une courbe de comparaison entre prévision du nombre de coupures calculé par le modèle de Cox et réalisation; montre une courbe de comparaison horaire du nombre de coupures calculé par le modèle de Cox et réalisation; représente un premier procédé en rapport avec l’invention; représente une méthode en rapport avec l’invention; représente une variante de la méthode de la ; représente un troisième procédé en rapport avec l’invention. Les dessins et la description ci-après contiennent des éléments pouvant non seulement servir à mieux faire comprendre la présente invention, mais aussi contribuer à sa définition, le cas échéant. Dans le cadre de la présente invention la problématique de quantifier l’impact des tempêtes en terme de nombre de coupures à l’échelle de tout un territoire est décomposée en un problème quantitatif : Combien de départs HT aériens seront coupés et combien d’incidents HT aériens se produiront lors de la prochaine tempête et en un problème de localisation : Quels seront les départs HT aériens coupés lors de la prochaine tempête ? On définit comme départ HT aérien le départ d’une ligne Haute Tension aérienne qui dans l’application concernée est une ligne moyenne tension (ligne dite HTA) aérienne. On considère qu’un départ HT aérien est coupé lorsque la ligne qui lui est reliée se coupe à un endroit quelconque ce qui est appelé incident. Il y aura donc une coupure de départ HT aérien correspondant à un ou plusieurs incidents ou coupures de la ligne HT aérienne associée. Les départs HT aériens sont regroupés dans des postes sources qui comportent un ou plusieurs départs. Le nombre de départs HT aériens coupés correspond au nombre de lignes HT aériennes qui subiront au moins un incident. Il est donc important de prévoir ces deux quantités. Selon un premier aspect, un premier algorithme basé sur un modèle à risque proportionnel de Cox connu de « Cox Proportional Hazard Regression for Survival Data », John Fox, novembre 2001 ; réalise une estimation de la distribution du nombre de départs HT aériens coupés. Pour ce premier algorithme, dans un premier temps, on récupère les données, on les met en forme dans une base de données spécifique et on apprend un modèle de Cox qui modélise la résistance des départs HT aériens aux rafales. Le modèle permet concrètement de calculer, pour chaque départ, la probabilité, à tout moment de la tempête prévue, que le départ soit coupé. A partir de ces probabilités, on peut calculer le nombre moyen de départs coupés sur une zone et une approximation de la distribution du nombre de départs coupés sur une zone. Selon un second aspect, un second algorithme de classification de type XGBoost tel que décrit dans le document Chen T. and Guestrin C., 2016, « XGBoost: A Scalable Tree Boosting System » dans « Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining (KDD ’16) » Association for Computing Machinery, New York, NY, USA, 785–794. DOI: https://doi.org/10.1145/2939672.2939785 va permettre de réaliser une estimation des départs les plus à risque d’être coupés. Pour le premier aspect, le modèle de Cox est à la base une classe de modèles de survie en statistiques qui permet de corréler le temps de survie à des covariables. Les modèles de survie étudient le temps écoulé avant qu'un événement ne survienne. Historiquement, le modèle de Cox a été appliqué en médecine. Au cours des années, l'utilisation du modèle s'est étendue à d'autres situations, en particulier dans le milieu industriel où les modèles de survie sont utilisés pour prévoir la durée avant une défaillance. Classiquement, le modèle de Cox exprime le taux de panne, noté λ en fonction du temps t et des covariables X1, …, Xn. λ o (t) est appelé risque de base. Il correspond au risque instantané de défaillance lorsque toutes les covariables sont nulles. La modèle de Cox a la particularité de ne pas avoir besoin de paramétrer cette fonction de risque. On dit que c’est un modèle à risque proportionnel car les risques sont évalués proportionnellement à une référence. Dans le cadre de la présente invention, un point important est que le temps a été remplacé par le maximum des rafales subies pour chaque départ HT aérien modélisé avant un instant t. Cela revient à modéliser la résistance des départs HT aériens aux rafales subies avant l’instant t. Dans l’étape d’apprentissage, cela revient à considérer le maximum des rafales subies par les départs HT aériens avant une coupure (cas non censuré en terme fiabiliste) ou, lorsqu’il n’y pas de coupure (cas censuré à droite en terme fiabiliste), le maximum des rafales subies par ces départs HT aériens. Cette étape 420 est schématisée en et prend en compte des données météo tempêtes passées collectées 400 et organisées une base de données 410 regroupant dans l’exemple soixante-six tempêtes subies sur le territoire français entre janvier 2014 et décembre 2019 ainsi que des données de description du réseau HT moyenne tension aérien sur ce territoire. Les données utilisées pour définir les covariables dans la base de données comprennent avantageusement pour les données météo des données telles que : - les données des rafales, vitesse, direction ; - la saison d’occurrence de la tempête ; - la température moyenne pendant la tempête issue des prévisions de tempête; - l’humidité moyenne pendant la tempête issue des prévisions de tempête; - la direction du vent pendant la tempête (Nord, Sud, Est ou Ouest). Pour la direction du vent, les données sont dans l’exemple de simulation réalisé issue des prévisions ARPEGE de Météo-France : ce paramètre est calculé en prenant la moyenne des directions (en degré) des vents lorsque les rafales dépassent le quantile à 75% des rafales prévues par Météo-France pendant la journée durant laquelle la tempête se produit. Cela revient à ne prendre en compte que les rafales importantes dans le calcul de la direction du vent. Les données propres à l’implantation des départs utilisées lors de l’apprentissage comportent dans le modèle réalisé: - la région où se situe le départ HT aérien ; - le type de boisement observé sur le département; - la longueur de la ligne HT reliée au départ ; - la longueur de la ligne HT soumise à un risque de boisement. Ces covariables X i peuvent être numériques, par exemple : longueur aérien, taux de boisement…, ou catégorielles, par exemple : région, type de boisement…. Le taux de défaillance s’exprime de la façon suivante : La référence pour une variable catégorielle est l’une des possibilités et la référence pour une variable continue est soit 0 ou sa moyenne. Dans notre cadre, une librairie dénommée « survival » du langage de programmation R utilisée pour l’apprentissage considère pour chaque covariable, la moyenne comme étant la référence. Les paramètres βi sont issus d’un calcul. Ils représentent les paramètres du modèle de Cox à estimer. Ils sont estimés par la méthode classique du Maximum de Vraisemblance. Ces paramètres βi qui représentent l’impact des covariables sur le risque de défaillance, s’interprètent de la façon suivante : X i augmente le risque de défaillance ; X i n’impacte pas le risque de défaillance ; X i réduit le risque de défaillance. Les covariables X i sont alors associées à chaque départ dans une base de données départ/covariables. Pour chaque départ, les données de tempêtes passées comportent une série temporelle de prévision de rafale sur la durée de la tempête à venir. Pour l’organisme météorologique Météo France utilisé pour l’étude, la série temporelle est par exemple une série de prévisions horaires sur un maillage du territoire français sur une centaine d’heures. La série est réactualisée toutes les six heures. La variable t prend dans ce cas les valeurs discrètes 0 heure, 1 heure, … jusqu’à la fin de la prévision de tempête. En outre, pour une tempête, les données sont distribuées selon un maillage de l’ordre de 10km. Par exemple la France est découpée en une grille régulière dont les carrés mesurent 10km de côté. Pour chaque départ HT aérien, un point est pris de la façon suivante : on construit le plus petit rectangle qui encadre l’ensemble des postes HT/BT reliés au départ. Le « point météo » retenu pour extraire les données de vent, température et humidité est choisi de la façon suivante : Si au moins un « point météo » est dans le rectangle, on sélectionne celui qui renvoie la plus grande valeur de rafale ; Si aucun point météo ne se trouve dans le rectangle, on sélectionne celui qui est le plus proche du centre du rectangle. A chaque nouvelle prévision, pour affiner la prévision, les calculs sont réitérés avec les données réévaluées sur la durée de la tempête prévue accessibles, durée croissante jusqu’à ce que la prévision englobe la durée totale de la tempête à venir. Un premier aspect déjà évoqué ci-dessus consiste à remplacer dans le modèle la variable aléatoire classique « le temps écoulé avant qu'une coupure du départ ne survienne » par la variable aléatoire R « la rafale maximale à laquelle le départ HT aérien résiste ». Pour ce calcul, dans l’exemple de réalisation, sont utilisées les prévisions de rafales pour chaque maille avec un pas horaire dans l’exemple. Le modèle de Cox est vu comme un modèle contrainte résistance dont la contrainte est la rafale maximale prévue pendant la tempête à l’endroit du départ HT aérien. Un départ survit à une tempête si la rafale maximale prévue pour ce départ, notée r, n’a pas été suffisante pour le couper, autrement dit la rafale maximale R à laquelle ce départ se coupera est supérieure à r. La probabilité de cet événement s’apparente à la fiabilité et se note : Comme représenté en , illustrant la courbe C de la probabilité de survie A des départs du territoire en fonction du maximum des rafales B subies entre le début de la tempête et un instant t, avec un intervalle de confiance D, une tempête ayant des rafales plus puissantes a plus de chance de couper les départs. Pour rA la rafale maximale prévue pour une tempête TA et pour rB la rafale maximale prévue pour une tempête TB tels que rA>rB nous avons ainsi : La phase d’apprentissage permet d’évaluer les covariables à prendre en compte et les paramètres βi à prendre en compte. La détermination des paramètres et des covariables pertinents a été, dans l’exemple, réalisé à partir d’une base de données regroupant les prévisions météo de Météo France pour 66 tempêtes entre janvier 2014 et décembre 2019 avec une description du réseau. L’apprentissage effectué a permis d’obtenir une estimation de paramètres βi appropriés en fonction des covariables retenues comme influant sur la prévision dans ce contexte qui comporte un découpage du territoire en régions en accord avec le découpage par directions régionales de l’opérateur du réseau. D’autres découpages et donc des paramètres et/ou covariables adaptés peuvent être utilisés ou nécessaires pour des territoires différents et/ou des réseaux organisés différemment. Ainsi, la forme du modèle de Cox retenue dans la présente invention n’est pas classique. En effet, comme vu plus haut, à l’origine le modèle de Cox a été construit comme un modèle de fiabilité dont la contrainte est le temps qui passe. Dans notre cas, le temps a été remplacé par le maximum des rafales subies avant un instant t pendant la durée de la tempête prévue. Le modèle de Cox ainsi utilisé permet alors de calculer, pour chaque départ HT aérien, la probabilité d’occurrence d’une coupure à chaque instant t donné de la tempête à venir. Un exemple de calcul est donné en . Pour un départ d, soit la rafale maximale prévue entre le début d’une tempête à venir et un instant t, et les covariables associées. A l’étape 100, la probabilité notée que le départ soit coupé lors de cette tempête entre le début de la tempête et l’instant t, se calcule de la manière suivante : Pour l’estimation des paramètres, r t correspond à la rafale maximale entre le début de l’événement et l’occurrence de la coupure du départ d. La probabilité p d , correspondant à la probabilité que le départ d soit coupé pendant la tempête, est affichée sur une carte pour les opérateurs utilisant le système de l’invention mais n’a pas vocation à être utilisée telle quelle. On se sert de cette probabilité dans la formule destinée à calculer le nombre de coupures qui sert à estimer les moyens humains et matériels à mettre en place avant la tempête. Pour ce faire la probabilité p d est mémorisée au moins temporairement dans le système de calcul associé au procédé. De l’équation précédente, on peut déduire simplement le nombre moyen de coupures dans la zone géographique ou territoire concerné. Dans l’exemple de réalisation en France ou sur une Direction Régionale ou dans un département. Soit N dep le nombre de départs HT aériens de la zone géographique considérée, soit la variable aléatoire représentant le nombre de coupures affectant lesdits départs HT aériens, le nombre moyen de coupures entre le début de la tempête à venir et l’instant t sur cet ensemble est calculé à l’étape 110 de la façon suivante : Le calcul explicite de la loi du nombre de coupure entre le début de la tempête et l’instant t est possible mais nécessite un travail important en implémentation et un temps de calcul long. En considérant que la coupure du départ d est le résultat d’une expérience de Bernoulli de paramètre notée , la loi de parmi une population de N dep départs peut être obtenue en calculant le produit de convolution des N dep lois de Bernoulli : Toutefois, pour ne pas avoir recours à ce calcul couteux, nous utilisons l’inégalité de Le Cam (Le Cam L., 1960, An Approximation Theorem for the Poisson Binomial Distribution, Pacific Journal of Mathematics, vol. 10, numéro 4, 1960 p. 1181–1197 qui permet d’approximer la loi de par une loi de Poisson de paramètre : Cette approximation est possible car : - Sachant les conditions météo, les coupures des départs HT aériens sont indépendantes entre elles ; - Le nombre important de départs (de l’ordre de 15000 dans le cadre de la simulation effectuée) permet de considérer que nous sommes dans un cadre asymptotique ; - Les probabilités de coupure sont globalement petites, sauf pour les départs soumis à de très fortes rafales. Un calcul simple des quantiles de cette loi permet d’obtenir un intervalle de confiance de l’estimation du nombre moyen de départs coupés dans une zone géographique 130. En prenant par exemple les quantiles de 0.005 et 0.995, l’intervalle de confiance correspond à 99% du support de la loi de . On considère que la moyenne est l’estimation du nombre de coupures calculé dans le procédé. Ceci permet de réaliser l’affichage et/ou la sauvegarde des données du nombre moyen de coupures des départs potentiellement impactés entre le début de la tempête et l’instant t et de la probabilité de coupure pour chaque départ. Le calcul de probabilité peut-être réitéré pour chaque valeur t comprise entre 0 et l’instant de fin de la tempête afin d’évaluer l’évolution de pendant la tempête, la mémoriser et/ou l’afficher 120. Les intervalles de confiance et les estimations du nombre de coupure peuvent être affichés dans un tableau et dans un graphique comme en pour aider à l’évaluation de la situation et de la pertinence de la prévision. Ensuite, à partir de la probabilité de coupure est réalisé un calcul du nombre d’incidents HT aériens entre le début de la tempête prévue et l’instant t. En effet, sachant les covariables telles que rafales, température, humidité, taux de boisement etc.. , les incidents successifs d’un départ HT aérien peuvent être considérés indépendants. Ainsi le nombre d’incidents HT d’un départ d entre le début de la tempête et l’instant t suit une loi de Poisson de paramètre dont le paramètre s’obtient de la manière suivante à l’étape 150: A l’étape 160, le nombre d’incidents HT aériens entre le début de la tempête et l’instant t sur une zone géographique contenant N dep départs HT aériens suit donc une loi de poisson de paramètre Le modèle de Cox permet ainsi d’estimer, pour chaque instant t pendant la tempête à venir, le nombre moyen de coupures ainsi que le nombre moyen d’incidents , étape 170, avec intervalles de confiance 130 , 175 pour une prochaine tempête prévue par un organisme de prévision météorologique tel que Météo France dans l’étude menée ce qui résout le premier problème soulevé. Dans la suite, , respectivement correspond au nombre de coupures, respectivement incidents, prévu sur une zone donnée pendant toute la tempête. Il est par ailleurs souhaitable de déterminer quels départs HT aériens seront coupés. Il s’agit dans ce cas d’un problème de classification. En effet, chaque départ peut se trouver dans une classe « coupé » ou une classe « non coupé ». Pour ce faire, selon un mode de réalisation est utilisée une méthode d’agrégation de modèles 300 (en anglais Gradient Boosting) et, en particulier, l’algorithme XGBoost dans le langage informatique R a été sélectionné pour cette étape. Cet algorithme est un algorithme d’apprentissage supervisé utilisant des arbres de régression additifs multiples dont le principe et de combiner les résultats d’un ensemble de modèles plus simple et plus faibles afin de fournir une meilleure prévision. L’algorithme va utiliser plusieurs modèles qui seront ensuite combinés pour obtenir un seul résultat. Un problème existant pour les procédés de classification est qu’ils nécessitent le calcul ou la connaissance d’un seuil qui ici n’est pas aisé à évaluer. Dans un mode de réalisation avantageux de l’invention, comme présenté dans la figure 5A, les résultats du modèle de Cox sont utilisés pour se passer du choix du seuil. Pour chaque départ, la méthode d’agrégation de modèles 300a, au préalable calée sur un jeu de données d’apprentissage par une étape d’apprentissage 430 à la figure 6, va calculer dans l’étape 310 de la figure 5A un score, entre 0 et 1, d’appartenance à la classe « coupé » de ce départ et on choisit d’utiliser la donnée pour s’affranchir de ce calcul de seuil. Pour ce faire, dans un premier mode de réalisation, on ordonne les départs par zone géographique, par exemple département, Direction Régionale ou France entière ; selon le score obtenu par la méthode d’agrégation de modèles et les premiers départs, nombre moyen de départs HT aériens coupés calculé par le modèle de Cox sur la zone considérée, sont considérés comme appartenant à la classe « coupé » et les suivants comme appartenant à la classe « non coupé » à l’étape 320. Ces départs sont alors affichés ou mémorisés à l’étape 330. Dans un second mode de réalisation selon la , afin d’améliorer la performance de l’algorithme de localisation, les résultats sont remontés à la maille des postes sources regroupant les départs HT aériens. L’algorithme cherche ainsi à localiser les Postes Sources dont au moins un départ HT aérien sera coupé. les départs HT aériens étant regroupés selon des postes source, le procédé comporte alors selon la variante de la : - au moins une étape 340 de calcul du nombre moyen de postes sources avec au moins un départ HT aérien coupé sur une population de postes sources d’au moins une partie du territoire: étant le nombre de départ HT aérien pour le poste source j, et, - au moins une étape 350 de calcul d’appartenance des postes source à l’une des deux classes « au moins un départ HT aérien coupé » et « aucun départ HT aérien coupé » au moyen de la méthode d’agrégation de modèles 300b produisant un ordonnancement des postes source, l’appartenance à la classe « au moins un départ HT aérien coupé » étant affectées aux postes source de rang le plus élevé selon ledit ordonnancement. Une fois le calcul effectué, les résultats sont utilisés pour une étape d’affichage et/ou de mémorisation desdits postes sources appartenant à la classe « au moins un départ HT aérien coupé » les opérateurs pouvant alors mieux définir les besoins d’intervention. Pour illustrer, prenons l’exemple de la tempête Elsa survenue entre le 19/12/2019 et le 20/12/2019. Le procédé de l’invention expérimenté sur ce cas pratique de simulation, prévoit en France métropolitaine pour les lignes aériennes HTA, 431 incidents HTA et 367 coupures HTA réparties sur 245 postes sources. En réalité, 402 incidents HTA se sont réellement produits avec 338 coupures HTA réparties sur 219 postes sources. Les prévisions obtenues par le procédé de la présente invention sont donc suffisamment proches de la réalité pour fournir une prévision utilisable. La montre un exemple de l’application du modèle de Cox pour évaluer un nombre d’incidents HT aériens selon la méthode en regard des incidents qui ont réellement eu lieu pour une tempête passée. Sur cette figure se trouve en abscisse le nombre de coupures et en ordonnée la région affectée avec en niveau de gris la vitesse en m/s 5 des rafales Les barres horizontales 3 correspondent à l’intervalle de confiance à 99% selon la méthode. La moyenne 4 de l’estimation prend la forme d’un trait blanc à peu près au centre de ces intervalles. Les points 1 et 2 correspondent aux nombres d’incidents ou de coupures réellement observées. Les points 1 correspondent à un nombre observé dans l’intervalle de confiance les points 2 correspondent à un nombre de coupures hors de l’intervalle de confiance. On remarque que les résultats sont très satisfaisants. La est une courbe qui représente la comparaison horaire entre une prévision (moyenne et intervalle de confiance à 99%) et la réalisation sur une tempête passée. Cette courbe montre l’évolution temporelle de l’impact de la tempête concernée en France métropolitaine. En abscisse se trouve le temps en heures et en ordonnée le nombre d’incidents. Les points 10 sont les coupures issues du modèle représentée par des carrés avec un intervalle de confiance 11. Les points 12 en forme de rond correspondent aux incidents réellement observés dans l’intervalle de confiance 11 et les points 13 en forme de diamant correspondent aux incidents réellement observés hors intervalle de confiance. L’évolution réelle semble décalée d’une ou deux heures mais la tendance a été bien anticipée. Ce décalage peut venir d’un décalage dans la prévision des rafales. L’utilisation du couplage entre le modèle de Cox et l’algorithme XGBoost pour localiser les coupures à la maille des postes sources sur le modèle de la tempête Elsa en France a permis de constater la performance de la localisation des coupures : 86% des Postes Sources ont été bien classés, 48% des coupures ont été anticipées et 56% des coupures annoncées par l’outil ne se sont pas produites. Les résultats de localisation obtenus par le procédé de l’invention sont ainsi très raisonnables compte tenu de l’aspect difficilement prévisible d’une tempête. L’invention ne se limite pas aux exemples décrits ci-avant, seulement à titre d’exemple, mais englobe toutes les variantes que pourra envisager l’homme de l’art dans le cadre de la protection recherchée. Procédé de prévision de risques météorologiques sur une infrastructure de réseau électrique d’un territoire comprenant un nombre N dep de départs haute tension aérien dits départs HT aériens d qui comprend pour une tempête à venir: - pour chaque départ HT aérien dudit territoire, au moins un calcul de probabilité (100) que ledit départ soit coupé selon l’équation tirée du modèle de Cox comportant comme contrainte la valeur de rafale maximale à subir par ledit départ entre le début de la tempête à venir et un instant t: avec des covariables d’environnement desdits départs HT aériens et βi les paramètres représentant l’impact desdites covariables sur le risque de défaillance desdits départs HT aériens, lesdites covariables d’environnement comportant des données prévisionnelles de ladite tempête selon un maillage dudit territoire et des données propres à l’implantation desdites lignes HT aériennes et desdits départs HT aériens ; - au moins une étape de calcul d’un nombre moyen de coupures entre le début de la tempête et l’instant t (110) sur le territoire ou une partie du territoire comportant une étape de sélection d’un nombre de départs de la zone géographique considérée selon la formule : pour laquelle est la variable aléatoire représentant le nombre de coupures entre le début de la tempête à venir et l’instant t desdits départs HT aériens ; - au moins un affichage et/ou une sauvegarde (120) des données du nombre moyen de coupures en fonction des départs potentiellement impactés et de la probabilité de coupure pour chaque départ, le calcul de probabilité étant réitéré pour chaque valeur de t comprise entre 0 et l’instant de fin de la tempête (190) afin d’évaluer l’évolution de pendant la tempête, la mémoriser et/ou l’afficher. Procédé selon la revendication 1, comportant : - au moins un calcul d’un intervalle de confiance ic (130) du nombre de départs coupés sur le territoire ou une partie du territoire entre le début de la tempête à venir et l’instant t au moyen d’une approximation de la loi de par une loi de Poisson de paramètre : - au moins un calcul des quantiles de cette loi pour obtenir ledit intervalle de confiance ic , * au moins une mémorisation dudit intervalle de confiance et/ou l’affichage dudit intervalle de confiance (140) ; ledit calcul dudit intervalle de confiance ic étant réitéré pour chaque valeur de t comprise entre 0 et l’instant de fin de la tempête (190) afin d’évaluer l’évolution dudit intervalle de confiance pendant la tempête, l’afficher et/ou le mémoriser. Procédé selon la revendication 1 ou 2, comportant, à partir de la probabilité de coupure , - une ou plusieurs étapes de calcul de nombres d’incidents sur la ligne HT aérienne (150) associée à un départ HT aérien entre le début de la tempête et l’instant t pour chaque départ HT aérien d au moyen d’une loi de Poisson de paramètre : - une ou plusieurs étapes de calcul de la loi du nombre d’incidents HT aériens (160) sur ledit territoire contenant les départs HT aériens entre 0 et l’instant t au moyen d’une loi de Poisson de paramètre : - une ou plusieurs étapes de sauvegarde et/ou d’affichage (170) de la moyenne E[N inc ] desdits nombres d’incidents estimés pour une utilisation en tant qu’outil de préparation de moyens d’intervention. Procédé selon l’une quelconque des revendications précédentes, pour lequel, pour chaque départ HT aérien, les covariables propres à l’implantation comportent au moins les données suivantes: une zone où se situe le départ HT aérien; le type de boisement observé sur au moins une portion de ladite zone; la longueur de la ligne HT aérienne reliée au départ; la longueur de la ligne HT aérienne soumise à un risque de boisement. Procédé selon l’une quelconque des revendications précédentes, pour lequel, pour chaque départ HT aérien, les covariables météo comportent au moins les données suivantes: la saison d’occurrence de la tempête; la température moyenne pendant la tempête issue des prévisions de tempête; l’humidité moyenne pendant la tempête issue des prévisions de tempête; la direction du vent pendant la tempête (Nord, Sud, Est ou Ouest). Procédé selon la revendication 5, pour lequel, les données prévisionnelles comportant la direction du vent, cette dernière est calculée en prenant la moyenne des directions des vents lorsque les rafales dépassent le quantile à 75% des rafales prévues pendant la journée durant laquelle la tempête est prévue. Procédé selon l'une quelconque des revendications précédentes, pour lequel les calculs sont réitérés (190) à chaque modification des données prévisionnelles de tempêtes sur ledit territoire. Procédé selon l’une quelconque des revendications précédentes, comportant, selon une méthode d’agrégation de modèles (300a), pour chaque départ HT aérien d’une zone considérée une étape (310) de calcul d’appartenance à l’une des deux classes « coupé » , « non coupé » , une étape (320) d’ordonnancement des départs et d’affectation à la classe « coupé » aux départs de rang le plus élevé selon ledit ordonnancement et une étape (330) d’affichage et/ou de mémorisation desdits départs appartenant à la classe « coupé » . Procédé selon l’une quelconque des revendications 1 à 7, pour lequel, les départs HT aériens étant regroupés selon des postes source, le procédé comporte selon une méthode d’agrégation de modèles (300b) : - au moins une étape (340) de calcul du nombre moyen de postes sources avec au moins un départ HT aérien coupé sur une population de postes sources d’au moins une partie du territoire: étant le nombre de départ HT aérien pour le poste source j, et, étant la probabilité de coupure du départ HT aérien pendant la tempête calculée par le modèle de Cox, et, - au moins une étape (350) de calcul d’appartenance des postes source à l’une des deux classes « au moins un départ HT aérien coupé » et « aucun départ HT aérien coupé » produisant un ordonnancement des postes source, l’appartenance à la classe « au moins un départ HT aérien coupé » étant affectées aux postes source de rang le plus élevé selon ledit ordonnancement et ; - au moins une étape d’affichage et/ou de mémorisation (360) desdits postes sources appartenant à la classe « au moins un départ HT aérien coupé » . Procédé selon la revendication 8 pour lequel les résultats sont remontés à une maille de postes sources auxquels appartiennent les départs HT aériens, l’algorithme comportant une localisation des postes sources dont au moins un départ HT aérien est classé coupé et un affichage d’une prévision de postes sources à départs HT aériens classés coupés. Procédé selon la revendication 8, 9 ou 10, pour lequel ladite méthode d’agrégation de modèles utilise un algorithme connu sous le nom de XGBoost dans le langage R. Procédé selon l’une quelconque des revendications 8 à 11, pour lequel la méthode prend en compte en tant que variables le ou les nombres d’heures pour lequel les départs HT aériens sont affectés par une ou plusieurs rafales supérieures à une ou plusieurs vitesses données. Procédé selon l’une quelconque des revendications 8 à 12, pour lequel la méthode prend en compte en tant que variables au moins l’une des données: nombre de coupures passées liées au vent; longitude et la latitude du départ HT aérien en remplacement de la région où se situe le départ HT aérien ; altitude ; valeur en degrés de la direction du vent. Procédé selon l’une quelconque des revendications précédentes comportant : - une ou plusieurs étapes préalable de collecte de données (400) d’environnement des lignes électriques hautes tension aériennes reliées aux dits départs HT aériens sur ledit territoire et l’organisation d’une base de données associant lesdites données d’environnement desdits départs HT aériens sous formes de covariables selon un maillage dudit territoire et, - une ou plusieurs étapes d’apprentissage (420), à partir d’une base de données (410) regroupant des données de prévisions météo d’une pluralité de tempêtes passées et les pannes réseau associées, de calcul de paramètres βi d’un modèle de Cox dont le taux de défaillance s’exprime de la façon suivante , les coefficients représentant l’impact desdites covariables sur le risque de défaillance en considérant ledit modèle de COX comme un modèle contrainte/résistance dont la contrainte est la rafale maximale subie par un départ avant une coupure ou au cours desdites tempêtes passées si le départ n’est pas coupé et pour lequel est le taux de défaillance lorsque toutes les covariables sont nulles.