Brevet FR3159639A1 — publié le 20250829

L’invention concerne un procédé de simulation d’une ferme éolienne dans un espace prédéterminé, pour laquelle on effectue au moins les étapes successives suivantes : On forme (GR) différentes grilles dans l’espace prédéterminé. Pour chaque grille, on détermine la production d’énergie annuelle moyenne d’une mini-ferme (AEP-mf); On choisit (Ch) quelques grilles permettant de maximiser la production d’énergie, Pour chaque grille, on détermine une première disposition (Alg1) du nombre prédéfini d’éoliennes sur la grille, On modifie la position (Alg2) des éoliennes sur la grille, On forme des deuxièmes grilles (Gr2) avec des espacements modifiés et on réitère les étapes d) et e), On forme des troisièmes grilles (Gr3) avec des directions modifiées et on réitère les étapes d) et e), On détermine une disposition définitive (Disp_F) des éoliennes dans l’espace prédéterminé, pour la construction (Const) de la ferme éolienne selon cette disposition. Figure 1 à publier METHODE AMELIOREE DE POSITIONNEMENT D’EOLIENNES AVEC CONTRAINTES D’ALIGNEMENT La présente invention concerne un procédé de simulation de ferme éolienne dans un espace prédéterminé, un procédé de positionnement des éoliennes d’une ferme dans un espace prédéterminé, une méthode de construction de ferme éolienne dans un espace prédéterminé ainsi qu’une ferme éolienne. Pour répondre aux enjeux environnementaux, les fermes éoliennes ont fait leur apparition. Ces fermes éoliennes sont constituées de plusieurs turbines éoliennes espacées les unes des autres dans un espace délimité. Cet espace délimité peut être sur terre ou en mer. On distingue ainsi des fermes éoliennes terrestres (également appelées fermes éoliennes « onshore ») et des fermes éoliennes « offshore », c’est-à-dire en mer. Les éoliennes de ces fermes sont généralement des éoliennes à axe de rotation horizontal qui disposent d’un système pour orienter l’axe de rotation horizontal dans le sens de la direction du vent, afin de maximiser l’énergie récupérée par l’éolienne. Parfois, l’éolienne est prévue pour s’orienter automatiquement dans la direction du vent. On connaît également les éoliennes à axe de rotation vertical, qui ont l’avantage de ne pas nécessiter d’orientation dans la direction du vent. Dans les fermes éoliennes, les sillages générés par les éoliennes peuvent avoir pour effet de diminuer la vitesse de vent en aval de l’éolienne et donc de diminuer l’énergie récupérée par d’autres éoliennes, notamment celles situées en aval, dans le sens du vent, des éoliennes générant ces sillages. Le positionnement des éoliennes dans la ferme est donc important afin de maximiser l’énergie récupérée par la ferme. Par ailleurs, dans un emplacement choisi pour implanter une ferme éolienne, les caractéristiques locales du vent peuvent varier. En effet, la direction et la vitesse du vent sont notamment des paramètres qui peuvent varier dans le temps à l’emplacement considéré. Ces caractéristiques peuvent notamment être obtenues par des capteurs positionnés à l’emplacement défini et maintenus à cette position pendant plusieurs mois ou plusieurs années, de manière à obtenir des données statistiques suffisantes pour caractériser la ressource en vent à l’emplacement choisi. Ces capteurs peuvent notamment être des anémomètres positionnés à une altitude suffisante (de l’ordre de 100m au-dessus du sol) pour caractériser le vent qui sera vu par les éoliennes (c’est-à-dire le vent qui est sensiblement au niveau de l’axe du rotor par exemple). Les données de vent peuvent également être obtenues par télédétection par laser, aussi appelé LiDAR (de l’anglais « Light Detection And Ranging »). La connaissance statistique du vent à l’emplacement considéré de la ferme permet notamment d’obtenir la distribution de la vitesse du vent (aussi appelée par la suite « répartition de la vitesse du vent »), la distribution de la direction du vent (aussi appelée par la suite « répartition de la direction du vent ») et la probabilité jointe d’occurrence d’une vitesse de vent dans une direction donnée. Afin de maximiser l’énergie annuelle produite par la ferme, il est donc nécessaire de positionner les éoliennes de manière optimale dans l’emplacement prévu pour la ferme. On entend par « énergie annuelle produite » ou « production d’énergie annuelle », l’énergie moyenne totale produite par la ferme, donc par l’ensemble des éoliennes de la ferme, pendant une année. Cette énergie moyenne, estimée en tenant compte des données statistiques du vent (répartitions de vitesses de vent, de direction de vent et des probabilités d’occurrence), est basée sur une durée d’un an, d’où le terme « annuelle », de manière à éviter une influence saisonnière qui pourrait fausser les résultats. En effet, le vent, que ce soit sa vitesse ou sa direction, peut varier fortement en fonction des saisons. L’énergie moyenne est obtenue grâce à la connaissance de la répartition de la direction du vent, de la répartition des vitesses de vent et de la probabilité d’occurrence d’une vitesse de vent dans une direction donnée. Bien évidemment, l’énergie annuelle produite dépend également des turbines éoliennes choisies et de leurs positionnements. Pour déterminer un bon positionnement des éoliennes dans l’emplacement prévu pour la ferme, plusieurs méthodes ont été mises au point. La demande de brevet CN105119320 concerne une méthode basée sur un algorithme d’évolution. Ce type d’algorithme est stochastique, utilisant des processus aléatoires. Ce type d’algorithme nécessite un nombre de tests important afin d’obtenir un résultat, ce qui génère un important temps de calcul, et nécessite une mémoire informatique importante. Les demandes de brevet CN102142103A, CN105139269 et US2016171401 concernent des méthodes de positionnement d’éoliennes basées sur des algorithmes génétiques. Les algorithmes génétiques sont une catégorie d’algorithmes d’évolution. Ils nécessitent une multitude d’évaluations et de croisements de ces évaluations à partir des données réelles (continues) possibles. De ce fait, les algorithmes génétiques sont complexes et nécessitent donc un temps de calcul important et une mémoire informatique importante. On connaît également les méthodes décrites dans les documents suivants : - Tao S., Xu Q., Feijoo A., Zheng G., Zhou J., 2020. Wind farm layout optimization with a three-dimensional Gaussian Wake model. Renewable Energy, Volume 159, October 2020, Pages 553-569 : Cette méthode nécessite que l’emplacement prévu pour la ferme soit délimité par un rectangle. En d’autres termes, cette méthode n’est pas adaptée pour des fermes comprenant des espaces non convexes ou même non connexes. - Antonini E.G.A., Romerp D.A., Amon, C.H., 2018. Continuous adjoint formulation for wind farm layout optimization. Applied Energy, Volume 228, 15 october 2018, Pages 2333-2345 : Cette méthode est basée sur une méthode analytique d’évaluation des effets de sillage. Cette méthode est donc une méthode d’optimisation continue, basée sur des calculs de CFD (de l’anglais « Computational Fluid Dynamics » signifiant « Calculs de dynamique des fluides »), qui nécessite des temps de calcul importants et une mémoire informatique importante. - Wagner M., Day J., Neuman F., 2012. A fast and effective local search algorithm for optimizing the placement of wind turbines. Renewable Energy, Volume 51, March 2013, pages 64-70 : L’application de cette méthode nécessite des emplacements délimités par un rectangle. De ce fait, elle n’est pas adaptée à des emplacements de formes non convexes ou non connexes. - Feng J., Shen W. Z., 2015. Solving the wind farm layout optimization problem using random search algorithm. Renewable Energy, Volume 78, June 2015, Pages 182-192 : Cette méthode nécessite des emplacements délimités par des polyèdres. De ce fait, elle n’est pas adaptée à des emplacements de formes non connexes. - Quan N., and Kim H.M., 2018. Greedy robust wind farm layout optimization with feasability guarantee. Engineering Optimization, September 6, 2018, pages 1152-1167. Cette méthode est basée sur un algorithme glouton qui prend en compte toutes les positions situées à une distance minimale des éoliennes déjà positionnées dans l’emplacement pour déterminer l’emplacement de l’éolienne suivante. Il nécessite donc un temps de calcul important et une mémoire informatique importante, notamment si la discrétisation de l’emplacement prévu pour la ferme est fine, si l’espace dans lequel on souhaite placer les éoliennes est grand et/ou si le nombre d’éoliennes à positionner est important. On connaît également la demande de brevet WO 2022/028 847 de la demanderesse qui concerne une méthode de construction d’une ferme éolienne. Dans cette méthode, les implantations des éoliennes sont des positions discrètes de l’espace prédéterminé. Bien que donnant satisfaction, cette méthode ne met pas en œuvre de contraintes d’alignement des éoliennes, ce qui ne facilite pas l’implantation et la maintenance des éoliennes. De plus, pour les éoliennes offshore et afin de faciliter la navigation des bateaux, l’alignement des éoliennes fait partie des contraintes de conception des fermes. La demande de brevet FR 22/09.630 de la demanderesse concerne une méthode de construction d’une ferme éolienne prenant en compte les contraintes d’alignement. Toutefois, cette méthode nécessite un temps de calcul important et des ressources informatiques (mémoire et processeurs) importantes. Ainsi, le problème technique de l’invention consiste à trouver un positionnement optimal d’éoliennes dans des emplacements de forme complexe, tels que les zones non convexes et/ou non connexes, de manière à maximiser l’énergie totale produite par la ferme éolienne, en minimisant le temps de calcul nécessaire et les ressources informatiques, et en respectant des contraintes d’alignement des éoliennes selon deux directions formant un angle non nul entre elles (deux directions non parallèles). Pour répondre au problème technique, l’invention concerne un procédé de simulation d’une ferme éolienne dans un espace prédéterminé pour la construction de ladite ferme éolienne dans ledit espace prédéterminé de manière à maximiser l’énergie annuelle produite, la ferme éolienne étant composée d’un nombre prédéfini d’éoliennes, à partir d’une première répartition discrète de vitesse de vent, d’une deuxième répartition discrète de direction de vent et d’une probabilité d’occurrence de chaque valeur discrète de vitesse de vent dans chaque valeur discrète de direction du vent des première et deuxième répartitions discrètes, dans lequel on effectue au moins les étapes successives suivantes par des moyens informatiques, tels qu’un ordinateur : Pour différents premiers espacements prédéfinis séparés d’un premier pas d’espacement, différents deuxièmes espacements prédéfinis séparés dudit premier pas d’espacement, différentes premières directions prédéfinies séparées d’un deuxième pas de direction, et différentes deuxièmes directions prédéfinies séparées dudit deuxième pas de direction, on forme différents premiers couples de premiers et deuxièmes espacements prédéfinis, et différents deuxièmes couples de premières et deuxièmes directions prédéfinies, puis pour les différents premiers couples et pour les différents deuxièmes couples, on forme des premières grilles dans l’espace prédéterminé, chaque première grille étant définie par une pluralité de points d’intersection entre des premières lignes orientées dans la première direction prédéfinie et des deuxièmes lignes orientées dans la deuxième direction prédéfinie, les premières lignes étant espacées dudit deuxième espacement prédéfini le long de la deuxième direction prédéfinie et les deuxièmes lignes étant espacées dudit premier espacement prédéfini le long de la première direction prédéfinie, chaque première grille comprenant au moins une maille délimitée par des points d’intersection. Pour chaque première grille dudit espace prédéterminé, on détermine la production d’énergie annuelle moyenne d’une mini-ferme composée d’éoliennes positionnées à l’ensemble des points d’intersection d’une maille ou plusieurs mailles connexes de la première grille, à partir de ladite première répartition discrète de vitesse de vent, de ladite deuxième répartition discrète de direction de vent et de ladite probabilité d’occurrence ; Pour chaque premier couple de premier et deuxième espacements prédéfinis, on identifie au moins un deuxième couple de première et deuxième directions prédéfinies qui maximise la production d’énergie annuelle de ladite mini-ferme, puis on forme un ensemble de deuxièmes couples regroupant les deuxièmes couples identifiés ; Pour chaque première grille définie par un premier couple de premier et deuxième espacements prédéfinis et par un deuxième couple de première et deuxième directions prédéfinies parmi ledit ensemble de deuxièmes couples, on détermine une première disposition du nombre prédéfini d’éoliennes, chaque éolienne étant positionnée sur un desdits points d’intersection de la première grille par un premier algorithme de positionnement ; puis, pour chaque première grille utilisée à l’étape d), on modifie la position des éoliennes parmi les points d’intersection de la première grille sur lesquels aucune éolienne n’est positionnée, afin d’améliorer la production d’énergie annuelle de manière à obtenir une disposition du nombre prédéfini d’éoliennes sur chaque première grille considérée et une production d’énergie annuelle pour chaque disposition et on retient une quantité prédéterminée de dispositions correspondant aux dispositions maximisant la production d’énergie annuelle et on retient chaque première grille associée à chaque disposition retenue ; Pour chaque première grille retenue à l’étape e), on diminue le premier pas d’espacement, de préférence de moitié, et on forme des deuxièmes grilles à partir du premier pas d’espacement diminué, puis on réitère les étapes d) et e) en utilisant les deuxièmes grilles en remplacement des premières grilles ; Pour chaque deuxième grille retenue à l’étape f), on diminue le deuxième pas de direction, de préférence de moitié, et on forme des troisièmes grilles à partir du deuxième pas de direction, puis on réitère les étapes d) et e) en utilisant les troisièmes grilles en remplacement des premières grilles ; On détermine une disposition définitive des éoliennes dans l’espace prédéterminé, ladite disposition définitive correspondant à la disposition finale obtenue à l’étape g) de la troisième grille qui maximise la production d’énergie annuelle, pour la construction des éoliennes selon la disposition définitive dans l’espace prédéterminé. De préférence, on réitère les étapes f) et/ou g) jusqu’à ce que le premier pas d’espacement atteigne un premier seuil et/ou jusqu’à ce que le deuxième pas de direction atteigne un deuxième seuil. Avantageusement, le premier seuil est compris entre un cinquantième du diamètre des éoliennes et un cinquième du diamètre des éoliennes ou inférieur à 1 m et/ou le deuxième seuil est inférieur à 2°, de préférence inférieur à 1°. De manière avantageuse, à l’étape f), on forme un nombre prédéterminé de deuxièmes grilles, chaque deuxième grille étant définie par des lignes de même orientations que la première grille retenue à l’étape e), le premier espacement de chaque deuxième grille correspondant à la somme du premier espacement de la première grille retenue à l’étape e) et d’un premier nombre entier de premier pas d’espacement diminué, le deuxième espacement de chaque deuxième grille correspondant à la somme du deuxième espacement de la première grille retenue à l’étape e) et dudit premier nombre entier de premier pas d’espacement diminué, de préférence, ledit premier nombre entier étant 1 ou -1. Selon une variante de l’invention, à l’étape g), on forme un nombre prédéterminé de troisièmes grilles, chaque troisième grille étant définie par des lignes de même espacements que la deuxième grille retenue à l’étape f), la première direction de chaque troisième grille correspondant à la somme de la première direction de la deuxième grille retenue à l’étape f) et d’un deuxième nombre entier du deuxième pas de direction diminué, la deuxième direction de chaque troisième grille correspondant à la somme de la deuxième direction de la deuxième grille retenue à l’étape f) et dudit deuxième nombre entier de deuxième pas de direction diminué, de préférence, ledit deuxième nombre entier étant 1 ou -1. Selon une configuration de l’invention, à l’étape e) de modification de la position des éoliennes pour chaque première grille, on réalise au moins les étapes suivantes : e1) pour chaque grille utilisée à l’étape d), on définit un ordre séquentiel de modification des positions des éoliennes déterminées par le premier algorithme de positionnement et on repositionne, une à une, chaque éolienne de l’ordre séquentiel défini, en trouvant un point d’intersection libre de la première grille qui maximise la production d’énergie annuelle à partir de ladite première répartition discrète de vitesse de vent, de ladite deuxième répartition discrète de direction de vent et de ladite probabilité d’occurrence ; e2) on réitère l’étape e1) autant de fois que nécessaire jusqu’à ce que, lors à la fin d’une étape e1), aucune éolienne n’ait été repositionnée et on obtient une disposition du nombre prédéfini d’éoliennes sur chaque première grille considérée et une production d’énergie annuelle pour chaque disposition finale ; et on retient une quantité prédéterminée de dispositions correspondant aux dispositions maximisant la production d’énergie annuelle et chaque première grille associée à chaque disposition retenue. Préférentiellement, ledit ordre séquentiel est obtenu de manière aléatoire et/ou on modifie l’ordre séquentiel à chaque itération de l’étape e1). Avantageusement, à l’étape e1), on choisit une sélection d’intersections parmi les points d’intersection sur lesquels aucune éolienne n’est positionnée, pour chaque éolienne à repositionner. De préférence, la sélection d’intersections correspond à une valeur prédéterminée de points d’intersections sur lesquels aucune éolienne n’est positionnée, ladite sélection d’intersections étant choisie de manière aléatoire ou étant choisie comme correspondant à ladite valeur prédéterminée de points d’intersections sur lesquels aucune éolienne n’est positionnée les plus proches de la position de l’éolienne à repositionner. Avantageusement, ledit espace prédéterminé comprend des zones non connexes et/ou non convexes. Selon une variante de l’invention, pour chaque première grille considérée, ledit premier algorithme de positionnement réalise au moins les étapes suivantes : - on définit arbitrairement la position de la première éolienne sur un des points d’intersection de la première grille ; - puis pour chaque éolienne à positionner, de manière successive : on définit des positions potentielles pour l’éolienne à positionner, lesdites positions potentielles étant constituées par les points d’intersection de la première grille situés entre une distance minimale et une distance maximale de toutes les éoliennes positionnées ; on calcule la production d’énergie annuelle des éoliennes positionnées et de l’éolienne à positionner pour les positions potentielles définies, à partir de ladite première répartition discrète de vitesse de vent, de ladite deuxième répartition discrète de direction de vent et de ladite probabilité d’occurrence ; On choisit la position de l’éolienne à positionner correspondant à la valeur maximale calculée de production d’énergie annuelle : - on détermine ladite première disposition correspondant à la position du nombre prédéfini d’éoliennes parmi les points d’intersection de la première grille considérée dans ledit espace prédéterminé. Selon une mise en œuvre de l’invention, à l’étape c), on fait la différence entre la production d’énergie annuelle de la mini-ferme obtenue à l’étape b) et la production d’énergie annuelle d’une seule éolienne multipliée par le nombre d’éoliennes de la mini-ferme. De manière avantageuse, à l’étape a), le premier pas d’espacement est compris entre deux fois et dix fois le diamètre des éoliennes, de préférence entre deux fois et six fois le diamètre des éoliennes. Préférentiellement, à l’étape a), le deuxième pas de direction est compris entre 1° et 20°, de préférence entre 5° et 10°. Avantageusement, à l’étape e), la quantité prédéterminée de dispositions est comprise entre 1 et 20, de préférence entre 1 et 10. De préférence, à l’étape c), on identifie un nombre prédéterminé de deuxième couple de première et deuxième directions prédéfinies pour chaque premier couple de premier et deuxième espacements prédéfini, de préférence le nombre prédéterminé est compris entre 1 et 20, de manière préférence entre 1 et 10. L’invention concerne également un procédé de positionnement des éoliennes d’une ferme éolienne dans un espace prédéterminé pour la construction de ladite ferme éolienne dans ledit espace prédéterminé, dans laquelle on mesure des données statistiques du vent, par des moyens de mesure, de préférence un capteur LIDAR, dans l’espace prédéterminé pour déterminer les première et deuxième répartitions discrètes et les probabilités d’occurrence de chaque vitesse de vent dans chaque direction du vent des première et deuxième répartitions discrètes, puis on met en œuvre le procédé de simulation selon l’une des configurations précédentes pour déterminer les positions des éoliennes de ladite ferme éolienne dans l’espace prédéterminé, de manière à positionner les éoliennes dans l’espace prédéterminé pour la construction de ladite ferme éolienne dans ledit espace prédéterminé. L’invention concerne également un procédé de construction d’une ferme éolienne dans lequel on met en œuvre le procédé de simulation selon l’une des variantes ou combinaisons de variantes décrites précédemment, ou on met en œuvre le procédé de positionnement des éoliennes tel que décrit précédemment, puis on construit la ferme éolienne en érigeant les éoliennes aux positions de ladite disposition définitive dans l’espace prédéterminé. L’invention concerne aussi une ferme éolienne obtenue à partir du procédé de construction d’une ferme éolienne tel que décrit précédemment. Liste des figures D'autres caractéristiques et avantages des méthodes et procédés ainsi que de la ferme selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après. [Fig 1] La FIG. 1 représente un exemple d’un procédé de simulation selon l’invention. [Fig 2] La FIG. 2 représente un exemple de réalisation des étapes de repositionnement des éoliennes du procédé de simulation selon l’invention. [Fig 3] La FIG. 3 représente un exemple d’emplacement complexe pour le positionnement d’éoliennes selon l’invention. [Fig 4] La FIG. 4 illustre une variante d’un premier algorithme de positionnement selon l’invention. [Fig 5] La FIG. 5 illustre la définition des contraintes d’alignement des éoliennes pour le procédé de simulation selon l’invention. [Fig 6] La FIG. 6 illustre le positionnement des éoliennes dans un espace prédéterminé non convexe utilisant le procédé de simulation selon l’invention. Afin de faciliter la lecture de la présente description, quelques définitions sont explicitées ci-après. On entend par « Algorithme glouton » (ou « greedy algorithm » en anglais) un algorithme qui consiste à établir étape par étape un optimum local. Dans le cas d’une ferme éolienne, il consiste à positionner chaque éolienne l’une après l’autre jusqu’à obtenir le positionnement de l’ensemble des éoliennes dans l’espace prédéterminé. On entend par « Algorithme d’évolution », un algorithme bio-inspiré dont l’idée consiste à faire évoluer un ensemble de solutions pour obtenir de meilleurs résultats. Ils sont donc stochastiques et utilisent itérativement des processus aléatoires. On entend par « Algorithme génétique », un algorithme d’évolution utilisant la notion de sélection naturelle. Ce type d’algorithme peut notamment croiser ou modifier certains paramètres des solutions précédentes en vue d’améliorer des résultats. On appelle « zones non connexes » des zones pour lesquelles il existe au moins deux points qui ne peuvent pas être reliés par un chemin continu entièrement contenu dans la zone en question. A l’inverse, une zone connexe est une zone pour laquelle chaque paire de points est reliée par un chemin continu entièrement contenu dans la zone. On appelle « zone convexe », une zone dans laquelle les segments reliant deux points quelconques de cette zone sont tous entièrement contenus dans la zone. Un cercle, un carré ou un rectangle délimitent tous par exemple des zones convexes. A l’inverse, une « zone non convexe » est une zone dans laquelle il existe au moins deux points reliés par un segment qui n’est pas entièrement contenu dans la zone. Par exemple, une zone délimitée par un cercle externe et un cercle interne concentriques n’est pas convexe. Dans la présente description, on entend par « grille » une pluralité de points d’intersection situés dans l’espace prédéterminé (ou à la frontière de cet espace). Les points d’intersection de la grille sont des points situés à l’intersection entre des premières lignes, toutes parallèles entre elles, et orientées dans la première direction prédéfinie et des deuxièmes lignes, également toutes parallèles entre elles, et orientées dans la deuxième direction prédéfinie. Les première et deuxième directions prédéfinies sont définies par rapport à un repère global ; elles sont de manière générale différentes de la direction du vent (mais peuvent en être colinéaires dans des conditions particulières). Les premières lignes sont espacées d’un deuxième espacement prédéfini le long des deuxièmes lignes. En d’autres termes, les premières lignes sont régulièrement espacées (elles sont toutes équidistantes) d’un pas correspondant au deuxième espacement le long des deuxièmes lignes. Les deuxièmes lignes sont espacées d’un premier espacement prédéfini le long des premières lignes. En d’autres termes, les deuxièmes lignes sont régulièrement espacées (elles sont toutes équidistantes) d’un pas correspondant au premier espacement le long des premières lignes. Ainsi, chaque grille comprend au moins une maille, de forme parallélogramme, délimitée par des points d’intersection entre les premières et les deuxièmes lignes (les premières et deuxièmes lignes étant sécantes entre elles). L’invention concerne un procédé de simulation d’une ferme éolienne dans un espace prédéterminé pour la construction de la ferme éolienne dans l’espace prédéterminé, dans le but de maximiser l’énergie annuelle produite. En d’autres termes, le procédé de simulation peut être utilisé pour positionner (et construire à ces positions) les éoliennes de la ferme dans l’espace prédéterminé de manière à maximiser l’énergie annuelle produite par la ferme. Grâce au procédé, on peut ainsi trouver les emplacements optimaux des éoliennes, en tenant compte de contraintes d’alignement. Une fois construite, la ferme permet de produire plus d’électricité grâce aux emplacements optimaux. L’espace prédéterminé peut être une zone terrestre ou une zone offshore en mer. Avantageusement, l’espace prédéterminé peut comprendre des zones non connexes. De ce fait, l’espace prédéterminé peut correspondre à des emplacements réels prévus pour l’implantation des éoliennes, par exemple une ferme prévue dans un emplacement comprenant deux zones séparées par une route de largeur significative (plusieurs mètres, voire plusieurs dizaines de mètres), telle qu’une autoroute, ou par une rivière ou un fleuve. Préférentiellement, l’espace prédéterminé peut comprendre des zones non convexes additionnellement ou alternativement aux zones non connexes. Ainsi, l’espace prédéterminé peut correspondre à des emplacements réels de forme complexe tels qu’un espace délimité par une colline ou une falaise escarpée, le littoral, le passage d’une rivière, un fleuve ou toute autre étendue d’eau. En offshore, l’espace prédéterminé peut comprendre des zones non convexes qui peuvent notamment être déterminées en tenant compte de la bathymétrie, de la nature du sol, des frontières avec d’autres pays, de canaux de navigation, de passage de câbles ou de pipeline par exemples. La ferme éolienne est composée d’un nombre prédéfini d’éoliennes. En d’autres termes, on a défini au préalable le nombre (prédéfini) d’éoliennes que l’on va implanter dans l’espace prédéterminé. La FIG. 3 illustre, de manière schématique et non limitative, un exemple d’espace prédéterminé adapté à la mise en œuvre du procédé de simulation, du procédé de positionnement et de la méthode de construction de l’invention. L’espace prédéterminé peut notamment comprendre une première zone Z1 et une deuxième zone Z2, zones représentées par les hachures verticales. Ces zones Z1 et Z2 sont non connexes. En effet, une distance minimale D non nulle apparaît entre les deux zones Z1 et Z2. De plus, la zone Z1 est rectangulaire. De ce fait, elle est convexe. La zone Z2 est de forme complexe, non convexe. En effet, si on considère les deux points A et B, on observe qu’une partie du segment Seg qui relie les points A et B est située en dehors de la zone Z2. Pour identifier les zones Z1 et Z2 de l’espace prédéterminé dans une troisième zone plus grande ZE, englobant ces deux zones Z1 et Z2, on peut utiliser une matrice booléenne. La troisième zone ZE est de forme rectangulaire, ce qui est plus facile à traiter informatiquement, que des zones non connexes et/ou non convexes comme Z1 et Z2. La matrice booléenne associe à chaque valeur discrète (position discrète) de la zone ZE une valeur égale à 1 si la position discrète est située dans la zone Z1 ou Z2 et une valeur 0 si elle est située en dehors de Z1 et Z2. Cette matrice booléenne permet de définir l’espace déterminé utilisé pour la méthode. On peut, à partir de cette matrice booléenne, déterminer la ou les frontières de l’espace prédéterminé. En effet, un point sera considéré comme faisant partie de la frontière si sa valeur dans la matrice booléenne est à 1, et s’il a parmi ses voisins directs, au moins un point qui a une valeur dans la matrice booléenne à 0. Pour le procédé de simulation de l’invention, on utilise une première répartition discrète de vitesse de vent, une deuxième répartition discrète de direction de vent et une probabilité d’occurrence de chaque valeur discrète de vitesse de vent dans chaque valeur discrète de direction du vent des première et deuxième répartitions discrètes. Ces données de vent sont considérées au niveau de l’espace prédéterminé correspondant. Ces données de vent peuvent notamment être obtenues par un moyen de collecte des données de vent, tel qu’un capteur LIDAR (de l’anglais « LIght Detection And Ranging » pouvant être traduit par télédétection par laser), un mât de mesure ou un anémomètre. Pour ce faire, l’invention peut concerner un procédé de positionnement, pour lequel on met en œuvre le procédé de simuler et le procédé de positionnement peut avantageusement comprendre une étape préalable de collecte, avant l’étape a), pendant laquelle on implante le moyen de mesure (ou moyen de collecte) des données de vent dans l’espace prédéterminé pendant une durée prédéterminée sur le site physique de l’espace prédéterminé (la région terrestre d’implantation prévue ou la zone offshore prévue) afin de collecter les données de vent du site, et, ainsi, on mesure les données statistiques du vent, par ces moyens de mesure, pour déterminer les première et deuxième répartitions discrètes de vitesse de vent et de direction de vent, et les probabilités d’occurrence de chaque vitesse de vent dans chaque direction de vent des première et deuxième répartitions discrètes. La mesure des données de vent peut être réalisée sur la durée prédéterminée, pour la collecte des mesures, qui peut être d’au moins un an, afin d’avoir les données relatives aux quatre saisons. Ainsi, on peut prévoir une étape de collecte des données statistiques du vent dans l’espace prédéterminé à partir d’au moins un moyen de collecte. Le moyen de mesure (de collecte) peut avantageusement être un capteur LIDAR. Pour ce procédé, on effectue au moins les étapes successives suivantes par des moyens informatiques, tels qu’un ordinateur ou un serveur ou un (super)calculateur : Formation de différentes grilles Calcul d’énergie annuelle pour une mini-ferme sur chaque grille Conservation de certaines grilles maximisant la production d’énergie annuelle de la mini-ferme Pour chaque grille, détermination d’une première disposition des éoliennes Modification des positions des éoliennes Création de nouvelles grilles, à partir des meilleures grilles obtenues précédemment, avec un nouvel espacement et réitération des étapes d) et e) Création de nouvelles grilles, à partir des meilleures grilles obtenues précédemment, avec de nouvelles directions et réitération des étapes d) et e) Détermination de la position définitive des éoliennes dans l’espace prédéterminé. L’utilisation de valeurs discrètes pour la vitesse du vent, la direction du vent et les points d’intersection discrets des différentes grilles de l’espace prédéterminé permet de simplifier le procédé de simulation, d’accélérer les temps de calculs et de limiter les ressources informatiques nécessaires (mémoire et processeur), comparativement à des méthodes utilisant des données réelles continues. En effet, des valeurs discrètes permettent de limiter le nombre de combinaisons possibles (on parle de méthode combinatoire), alors que des données réelles continues (on parle de méthode continue) fournissent une infinité de solutions. La combinaison des première et deuxième répartitions discrètes et des points d’intersection de la grille permet ainsi d’obtenir une bonne précision du positionnement des différentes éoliennes dans l’espace prédéterminé en tenant compte des contraintes d’alignement, tout en limitant le temps de calcul nécessaire à cette détermination des positions. La probabilité d’occurrence de chaque vitesse de vent, dans chaque direction de vent est notamment utilisée pour le calcul de la production d’énergie annuelle moyenne. Cette probabilité peut notamment provenir d’une rose des vents correspondant à l’espace prédéterminé, cette rose des vents étant bien connue de l’homme du métier. Pour établir cette rose des vents, on peut notamment utiliser des moyens de collecte évoqués précédemment, comme un anémomètre positionné sur un mât à une altitude suffisante (entre 80 m et 120 m par exemple pour être sensiblement au niveau du moyeu de l’éolienne par exemple), ou par l’intermédiaire d’un capteur LiDAR (acronyme de langue anglaise issu de « Light Detection And Ranging » signifiant « détection et estimation de la distance par la lumière »), positionné à proximité du sol et orienté vers la verticale. Ce moyen de collecte est maintenu en place pendant une longue période, plusieurs mois et idéalement plus d’un an, de manière à pouvoir prendre en compte les variations des caractéristiques du vent en fonction des saisons. Selon un mode de réalisation, la production d’énergie annuelle peut être estimée par la formule suivante : Où aep est la production d’énergie annuelle de la ferme d’éoliennes et est l’espérance de la puissance totale produite par la ferme éolienne (par le nombre prédéfini d’éoliennes dans l’espace prédéterminé) par rapport à la loi de probabilité jointe de vitesse de vent et de direction de vent . De ce fait, la puissance totale produite par la ferme prend en compte une répartition statistique de chaque vitesse de vent w s dans chaque direction du vent w p , par exemple par une distribution de Weibull. Dans le cas où les rotors des éoliennes sont systématiquement alignés dans un plan orthogonal à la direction du vent w p , la puissance totale produite par la ferme peut s’écrire : Où N est le nombre prédéfini d’éoliennes de la ferme dans l’espace prédéterminé, Pf est la puissance instantanée fournie par chaque éolienne f, dans la ferme, pour chaque vitesse de vent w s dans chaque direction du vent w p . Dans le cas où certains rotors des éoliennes seraient dans un plan désaxé du plan orthogonal à la direction du vent w p (autrement dit, l’éolienne n’est pas orientée face au vent mais est désaxée de celui-ci), un facteur correctif peut être pris en compte pour tenir compte de l’effet du désaxage. Ce facteur correctif peut notamment être issu de simulations CFD (de l’anglais « Computational Fluid Dynamics » signifiant « Calculs de dynamique des fluides »). La puissance instantanée Pf de chaque éolienne f dans la ferme peut s’écrire : Avec ρ est la masse volumique de l’air, S est la surface balayée par le rotor de l’éolienne f, v f (w s ) la vitesse du vent au niveau du rotor de l’éolienne f et le coefficient de puissance C pf de l’éolienne f dépendant de la vitesse v f (w s ) du vent au niveau du rotor de l’éolienne f, le coefficient de puissance C pf étant une caractéristique de l’éolienne f. En effet, les effets de sillage des éoliennes situées en amont (dans le sens du vent) et/ou sur le côté de l’éolienne f peuvent impacter la production d’énergie de l’éolienne f. Ce sillage peut générer une diminution de la vitesse du vent au niveau de l’éolienne f et/ou des turbulences du vent. Ces effets de sillage ont pour impact que la vitesse v f au niveau du rotor de l’éolienne ne correspond plus à la vitesse w s et que le coefficient de puissance C pf est alors également impacté. Les impacts des effets de sillage pris en compte dans l’équation [Math3] peuvent notamment être basés sur des modèles de sillage. Ces modèles de sillage peuvent notamment traduire : - une réduction de la vitesse du vent en amont de l’éolienne f du fait du sillage issu d’une éolienne en amont. Un tel modèle, bien connu de l’homme du métier, est le modèle de sillage de Jensen (décrit dans la publication « A simple model for Cluster Efficiency », Katic, Hojstrup and Jensen, EWEC 1986, notamment au paragraphe 2.1 de cette publication), - une augmentation de l’intensité turbulente du vent, - et/ou une superposition de sillages de plusieurs éoliennes en amont sur une même éolienne f, comme décrit dans la publication « A note on wind generator interaction » Jensen, DTU, 1983. Cette superposition de sillage peut combiner les effets de réduction de vitesse comme dans la publication citée précédemment et/ou d’intensité turbulente du vent de plusieurs sillages. Les modèles de sillage peuvent également être déterminés à partir de calculs CFD (de l’anglais « Computational Fluid Dynamics » signifiant « Calculs de dynamique des fluides »). Les formules précédentes peuvent être utilisées pour calculer l’énergie annuelle produite par la ferme composée du nombre prédéfini d’éoliennes, par la mini-ferme ou par une ferme constituée d’une seule éolienne. La FIG. 1 illustre, de manière schématique et non limitative, le procédé de simulation de la ferme d’éoliennes dans un espace prédéterminé Esp, de manière globale, le procédé de simulation étant notamment utilisé pour déterminer la position des éoliennes de la ferme dans l’espace prédéterminé, de manière optimale pour maximiser l’énergie annuelle produite par la ferme. A partir de l’espace prédéterminé Esp choisi (une zone géographique offshore ou onshore par exemple) qui peut être un espace comprenant des zones non-convexes et/ou non-connexes (par exemple avec une route ou une rivière qui traverse cet espace), et différents premiers couples C1 de premier et deuxième espacements prédéfinis et différents deuxièmes couples C2 de première et deuxième directions prédéfinies, on établit des grilles (GR) dans l’espace prédéterminé. Chacune de ces grilles GR est établie à partir d’une multitude de premières lignes parallèles entre elles et régulièrement espacées et de deuxièmes lignes parallèles entre elles, régulièrement espacées et sécantes avec les premières lignes (autrement dit, l’angle formé entre les premières lignes et les deuxièmes lignes est non nul et non plat). Les différentes grilles GR se distinguent entre elles par la combinaison des premiers couples C1 et deuxièmes couples C2. Pour chaque grille GR, les premières lignes sont orientées d’un angle correspondant à la première direction prédéfinie par rapport à une direction d’un repère fixe prédéfini (par exemple un repère avec une direction nord et une direction est) et les deuxièmes lignes sont orientées d’un angle correspondant à la deuxième direction prédéfinie par rapport à la même direction du repère fixe prédéfini (la direction Nord par exemple). Les premières lignes sont espacées du deuxième espacement prédéfini le long des deuxièmes lignes et les deuxièmes lignes sont espacées du premier espacement prédéfini le long des premières lignes. Une première ligne et une deuxième ligne peuvent être positionnées arbitrairement dans l’espace, ces lignes servant alors de référence aux positions des autres lignes. Les points d’intersection entre les premières lignes et les deuxièmes lignes qui sont situés dans l’espace prédéterminé (à l’intérieur de l’espace prédéterminé, la frontière de cet espace prédéterminé étant avantageusement comprise ou au contraire pouvant être exclue) définissent la grille. La grille peut ainsi comprendre des mailles, et les mailles de la grille forment des parallélogrammes compte-tenu du parallélisme des premières et deuxièmes lignes, délimitées par lesdites lignes et les points d’intersection. Pour chaque grille GR définie, on détermine la production d’énergie annuelle d’une mini-ferme AEP-mf. La production d’énergie annuelle est basée sur les données statistiques du vent issue par exemple d’un capteur LIDAR ou de tout autre moyen de mesure des données de vent, ces données statistiques comprenant une première répartition discrète RD1 des vitesses de vent, une deuxième répartition discrète RD2 des directions de vent et une probabilité d’occurrence Prob de chaque valeur de vitesse de vent de la première répartition discrète RD1 dans chaque direction de vent de la deuxième répartition discrète RD2. La mini-ferme est une ferme composée d’une portion de grille avec au moins une maille de la grille (ou de plusieurs mailles connexes) et où, à tous les points d’intersection de cette portion de grille, on positionne une éolienne. Ainsi, si la portion de grille comporte une seule et unique maille, la mini-ferme sera composée de quatre éoliennes, une à chaque point d’intersection de la maille. Cette situation avec une mini-ferme de quatre éoliennes est avantageuse car elle permet d’évaluer l’impact du sillage des éoliennes entre elles sur la grille de manière simple et rapide et sans nécessiter beaucoup de mémoire ni de temps de calcul lorsque cette étape est réalisée par ordinateur. Au sens de la présente description, une maille s’entend par une maille élémentaire : cela signifie que l’on considère un parallélogramme délimité par quatre points d’intersection comme étant une maille si ce parallélogramme ne comprend aucun autre parallélogramme défini par quatre points d’intersection, à l’intérieur de lui-même. En comparant les valeurs de production d’énergie annuelle de différentes grilles comprenant le même premier couple C1 et différents deuxièmes couples C2, on peut identifier un ou plusieurs deuxièmes couples C2 les plus prometteurs (donnant la ou les plus grandes productions d’énergie annuelle). On forme un ensemble de deuxièmes couples regroupant tous les deuxièmes couples C2 identifiés (les plus prometteurs). En d’autres termes, cet ensemble regroupe les deuxièmes couples les plus prometteurs. Ainsi, pour chaque premier couple C1 de premier et deuxième espacements prédéfinis, on identifie un ou plusieurs deuxièmes couples C2 que l’on ajoute à l’ensemble. Ainsi, on retient les couples de directions les plus prometteurs pour la suite, ce qui permet de limiter le nombre de configurations à tester. Les comparaisons de valeurs de production d’énergie annuelle peuvent être faites directement (par comparaison directe des différentes valeurs, on cherche alors les deuxièmes couples C2 qui maximisent l’énergie annuelle produite), ou indirectement. Dans ce dernier cas, on compare l’énergie annuelle produite par la mini-ferme à l’énergie annuelle produite par le même nombre d’éoliennes que celui de la mini-ferme comme si chaque éolienne était indépendante (donc sans tenir compte des effets de sillage) : on va alors estimer la perte liée aux effets de sillage de la mini-ferme. Dans ce cas, bien-entendu, on retient le ou les deuxièmes couples C2 qui minimisent la perte liée aux effets de sillage. Ainsi, on peut choisir Ch, à partir de la production d’énergie annuelle des mini-fermes, au moins un deuxième couple C2. Ainsi, on sélectionne un certain nombre de grilles parmi celles qui ont été établies préalablement. Les grilles retenues correspondant à chaque premier couple C1 et à un deuxième couple de l’ensemble formé (c’est-à-dire qu’on choisit le deuxième couple parmi les couples de direction les plus prometteurs) sont utilisées par la suite. En d’autres termes, l’utilisation de la mini-ferme permet de présélectionner des grilles pour la suite (et d’éliminer les autres grilles qui sont les moins prometteuses). Puis, pour chaque grille sélectionnée, on établit une première disposition des éoliennes Alg1 sur cette grille (sur des points d’intersection de cette grille). Cette première disposition cherche à obtenir un bon rendement énergétique et donc prend en compte les première et deuxième répartitions discrètes RD1 et RD2 et la probabilité d’occurrence Prob. A partir de cette première disposition des éoliennes sur chaque grille sélectionnée, on cherche ensuite à améliorer la récupération d’énergie. On repositionne alors les éoliennes Alg2 sur chaque grille en cherchant, à chaque repositionnement d’une éolienne, à trouver la position sur la grille qui permet d’obtenir l’énergie annuelle produite la plus élevée. Bien-entendu, on tient compte, pour ce faire, des première et deuxième répartitions discrètes RD1 et RD2 et de la probabilité d’occurrence Prob. Une fois que toutes les éoliennes ont été repositionnées, on retient AEP_g une quantité prédéterminée de dispositions maximisant la production d’énergie annuelle et chaque grille associée Grr de chacune de ces dispositions retenues. En d’autres termes, on retient un certain nombre de dispositions (le nombre étant fixé au départ) qui sont les dispositions les plus prometteuses pour la production d’énergie annuelle et les grilles associées (on élimine les autres grilles, ce qui permet encore de limiter les grilles les plus prometteuses pour la suite). A partir de chaque grille Grr associée aux dispositions retenues, on diminue MΔr le premier pas d’espacement (par exemple, on le divise par deux) et on forme Gr2 des deuxièmes grilles à partir du premier pas d’espacement diminué et de chacune des grilles Grr associées aux dispositions retenues. Chaque deuxième grille formée comporte le même deuxième couple de première et deuxième directions et on génère un nouveau premier couple de premier et deuxième espacement. Le premier espacement de chaque deuxième grille est déterminé à partir du premier espacement de la grille Grr associée et à partir du premier pas d’espacement diminué. Par exemple, il peut correspondre à la somme du premier espacement de la grille Grr associée et d’un premier nombre entier (de préférence égal à 1 ou à -1) du premier pas d’espacement diminué. De la même manière, le deuxième espacement de chaque deuxième grille est déterminé à partir du deuxième espacement de la grille Grr associée et à partir du premier pas d’espacement diminué. Par exemple, il peut correspondre à la somme du deuxième espacement de la grille Grr associée et d’un premier nombre entier (de préférence égal à 1 ou à -1) du premier pas d’espacement diminué. Lorsque le premier nombre entier peut être égal à 1 ou à -1, on peut former ainsi quatre deuxièmes grilles pour chaque grille Grr associée à une disposition retenue. Puis on réitère les étapes Alg1 d’établissement d’une première disposition des éoliennes sur chaque deuxième grille et Alg2 de repositionnement de chaque éolienne sur chaque deuxième grille et on retient AEP_g une quantité prédéterminée de dispositions maximisant la production d’énergie annuelle et chaque deuxième grille associée à chacune de ces dispositions retenues, chaque deuxième grille associée devant alors une grille associée Grr. A partir de chaque grille Grr associée aux dispositions retenues (c’est-à-dire de deuxième grille retenue après l’étape de repositionnement des éoliennes sur chaque deuxième grille), on diminue MΔθ le deuxième pas de direction (par exemple, on le divise par deux) et on forme Gr3 des troisièmes grilles à partir du deuxième pas de direction diminué et de chacune des grilles Grr associées aux dispositions retenues (c’est-à-dire de deuxièmes grilles retenues). Chaque troisième grille formée comporte le même premier couple de premier et deuxième espacements que chaque grille Grr et on génère un nouveau deuxième couple de première et deuxième directions. La première direction de chaque troisième grille est déterminée à partir de la première direction de la grille Grr associée (deuxième grille retenue) et à partir du deuxième pas de direction diminué. Par exemple, elle peut correspondre à la somme de la première direction de la grille Grr associée (deuxième grille retenue) et d’un deuxième nombre entier (de préférence égal à 1 ou à -1) du deuxième pas de direction diminué. De la même manière, la deuxième direction de chaque troisième grille est déterminée à partir de la deuxième direction de la grille Grr associée (deuxième grille retenue) et à partir du deuxième pas de direction diminué. Par exemple, elle peut correspondre à la somme de la deuxième direction de la grille Grr associée (deuxième grille retenue) et d’un deuxième nombre entier (de préférence égal à 1 ou à -1) du deuxième pas de direction diminué. Lorsque le deuxième nombre entier peut être égal à 1 ou à -1, on peut former ainsi quatre troisièmes grilles pour chaque grille Grr associée à une disposition retenue (deuxième grille retenue). Bien entendu, on peut réitérer plusieurs fois les étapes : De diminution MΔr du premier pas d’espacement et de formation Gr2 des deuxièmes grilles, puis des étapes d’établissement d’une première disposition Alg 1 et de repositionnement des éoliennes Alg2 pour retenir AEP_g de nouvelles grilles à partir de nouvelles deuxièmes grilles et/ou ; De diminution MΔθ du deuxième premier pas de direction et de formation Gr3 des troisièmes grilles, puis des étapes d’établissement d’une première disposition Alg 1 et de repositionnement des éoliennes Alg2 pour retenir AEP_g de nouvelles grilles à partir de nouvelles troisièmes grilles ; Par exemple, on peut réitérer plusieurs fois ces étapes jusqu’à ce que le premier pas d’espacement atteigne un premier seuil et/ou jusqu’à ce que le deuxième pas de direction atteigne un deuxième seuil. On détermine ensuite la disposition définitive Disp_F des éoliennes dans l’espace prédéterminé Esp comme étant la disposition finale des éoliennes sur la grille retenue qui permet d’obtenir l’énergie annuelle produite maximale parmi les dernières grilles retenues AEP_g (après la dernière étape de repositionnement des éoliennes Alg2 réalisée). On peut alors construire Const les éoliennes aux emplacements prévus correspondant à la disposition définitive Disp_F déterminée dans l’espace prédéterminé. Cette étape est représentée en traits pointillés pour indiquer son caractère optionnel. Lorsque cette étape est réalisée, l’invention concerne alors une méthode de construction des éoliennes d’une ferme éolienne dans un espace prédéterminé, de manière à optimiser l’énergie annuelle produite par la ferme. Etape a) : Formation de grilles Cette étape permet la définition de plusieurs grilles possibles à l’intérieur de l’espace prédéterminé. Chaque grille a donc une limite spatiale et ne peut pas dépasser de l’espace prédéterminé. Cette grille comprend une ou plusieurs mailles séparées par des points d’intersection. Ces mailles sont de forme parallélogramme, par exemple rectangulaire ou carrée, et les points d’intersection sont à la frontière entre les mailles. Toutes les mailles et tous les points d’intersection de chaque grille sont situés à l’intérieur de l’espace prédéterminé. Les mailles de chaque grille sont définies entre les premières lignes et les deuxièmes lignes. Ces grilles constituent chacune une option possible pour placer les éoliennes de la ferme dans l’espace prédéterminé, en tenant compte de contraintes d’alignement (selon les premières et deuxième lignes). Les points d’intersection de chaque grille constituent des positions possibles d’implantation d’une éolienne (à chaque point d’intersection, on ne peut positionner qu’une seule éolienne). En utilisant plusieurs grilles, on peut tester différentes contraintes d’alignement (différentes valeurs de premier et deuxième espacements prédéfinis et différentes valeurs de première et deuxième directions prédéfinies). Les première et deuxième directions sont non colinéaires et forment un angle non nul (et non plat) de manière à considérer des contraintes d’alignement en deux dimensions dans l’espace de manière à ce que les premières lignes et les deuxièmes lignes soient sécantes. En utilisant différentes valeurs discrètes, on limite le nombre de combinaisons possibles et donc la mémoire informatique utile. De plus, ces valeurs discrètes sont suffisantes pour assurer une précision compatible avec la construction des éoliennes sur site (compte tenu des précisions possibles de construction). Les premier et deuxième espacements prédéfinis sont avantageusement compris entre une distance minimale et une distance maximale, ces distances minimale et maximale dépendant par exemple du diamètre du rotor de l’éolienne. La distance minimale correspond à la valeur minimale dont il faut séparer deux éoliennes successives afin de limiter les pertes énergétiques liées aux effets de sillage et la distance maximale correspond à la valeur maximale dont il faut séparer deux éoliennes, une distance supérieure n’entraînant plus de gain sur la récupération d’énergie et limitant la possibilité d’implanter assez d’éoliennes dans l’espace prédéterminé, pour assurer un rendement énergétique suffisant. En effet, des premier et/ou deuxième espacements prédéfinis inférieurs à la distance minimale ou supérieure à la distance maximale ne s’avèrent pas nécessaires puisque les éoliennes doivent respecter une distance minimale entre elles pour éviter des pertes énergétiques liées aux effets de sillage trop importantes et une distance maximale pour implanter le nombre prédéfini d’éoliennes et obtenir un rendement énergétique intéressant de la ferme dans l’espace. De préférence, la distance minimale peut être supérieure ou égale à deux fois le diamètre du rotor des éoliennes, et de manière préférée, supérieure à quatre fois le diamètre du rotor des éoliennes ; la distance maximale peut être d’au moins quatre fois le diamètre des éoliennes et de préférence au moins huit fois le diamètre des éoliennes. De préférence, la première direction prédéfinie peut être comprise entre -90° et 90°, compte-tenu de la symétrie, par rapport à un repère orthonormé, par exemple le repère terrestre avec le nord situé à 90° et le sud à -90°, ce repère orthonormé étant de préférence celui utilisé pour définir les directions du vent. Et de manière préférée, la deuxième direction prédéfinie peut être comprise entre -90° et la première direction. De ce fait, avec les effets de symétrie, on balaye tout l’espace possible en limitant le nombre de combinaisons. Pour différents premiers espacements prédéfinis séparés d’un premier pas d’espacement, différents deuxièmes espacements prédéfinis séparés du (même) premier pas d’espacement, différentes premières directions prédéfinies séparées d’un deuxième pas de direction, et différentes deuxièmes directions prédéfinies séparées du (même) deuxième pas de direction, on forme différents premiers couples de premiers et deuxièmes espacements prédéfinis, et différents deuxièmes couples de premières et deuxièmes directions prédéfinies. Puis pour les différents premiers couples et pour les différents deuxièmes couples, on forme des grilles dans l’espace prédéterminé. Chaque grille est définie par une pluralité de points d’intersection entre des premières lignes orientées dans la première direction prédéfinie et des deuxièmes lignes orientées dans la deuxième direction prédéfinie. Les premières lignes sont espacées (deux à deux) du deuxième espacement prédéfini le long de la deuxième direction prédéfinie et les deuxièmes lignes sont espacées (deux à deux) du premier espacement prédéfini le long de la première direction prédéfinie. La grille comprend au moins une maille délimitée par des points d’intersection (chaque maille est délimitée par quatre points d’intersection et certains points d’intersection peuvent ne pas être rattachés à une maille). Ainsi, les points d’intersection de la grille forment un maillage discret. Par maillage discret, on entend que les points d’intersection du maillage qui constituent la grille sont des valeurs discrètes (par opposition à des valeurs continues). Dans cette étape a), on peut choisir un premier pas d’espacement et/ou un deuxième pas de direction grossiers (par exemple, le premier pas d’espacement peut correspondre d’au moins la distance minimale, soit d’au moins deux fois le diamètre des éoliennes et/ou le deuxième pas de direction peut être compris entre 5° et 10°). En effet, les grilles formées à cette étape a) sont utilisées dans la suite pour sélectionner des premiers couples et deuxièmes couples offrant les meilleurs compromis. Puis de nouvelles grilles sont élaborées en étapes e) et f) afin de raffiner le maillage et ainsi augmenter la précision optimale du positionnement des éoliennes et augmenter l’énergie annuelle produite. Par exemple, à l’étape a), le premier pas d’espacement (grossier) peut être compris entre deux fois et dix fois le diamètre des éoliennes (c’est-à-dire le diamètre du rotor des éoliennes), de préférence entre deux fois et six fois le diamètre des éoliennes. Avantageusement, à l’étape a), le deuxième pas de direction (grossier) peut être compris entre 1° et 20°, de préférence entre 5° et 10°. La FIG. 5 illustre, de manière schématique et non limitative, le changement de repère entre le repère (O1, N, E) qui correspond au repère global, O1 pouvant être basé arbitrairement, N peut correspondre à la direction Nord et E à la direction Est, et un repère lié à la grille avec O2 un point fixé arbitrairement sur un point d’intersection entre une première ligne Lig1 et une deuxième ligne de la grille Lig2, la première ligne Lig1 et la deuxième ligne Lig2 étant sécantes entre elles. La première ligne Lig1 forme ici un angle θ1 avec la direction E et la deuxième ligne Lig2 forme un angle θ2 avec la direction E (alternativement, on pourrait aussi avoir la première ligne Lig1 formant un angle θ1 avec la direction N et la deuxième ligne Lig2 formant un angle θ2 avec la direction N). Le repère (O1, N, E) peut avantageusement être celui utilisé pour déterminer les directions du vent. Etape b) : Calcul d’énergie annuelle pour une mini-ferme sur chaque grille L’étape a) a permis de définir un panel de grilles où les points d’intersection de chaque grille sont des points susceptibles de correspondre à l’implantation d’une (seule) éolienne. Bien entendu, afin de vérifier les contraintes d’alignement des éoliennes, il n’est pas possible, dans le présent procédé, d’utiliser des points d’intersection d’une grille pour implanter une partie des éoliennes avec des points d’intersection d’une autre grille pour implanter d’autres éoliennes. Afin d’éviter de rechercher la disposition optimale des éoliennes sur chaque grille, en calculant à chaque fois la production d’énergie annuelle, il est préférable de faire un tri pour n’utiliser que les grilles susceptibles d’offrir les meilleures possibilités de trouver la disposition optimale. L’étape b) cherche à effectuer ce premier tri parmi les différentes grilles. Pour ce faire, pour chaque grille formée à l’étape a), on considère une mini-ferme constituée d’éoliennes à chaque point d’intersection de la portion de grille considérée. La portion de grille correspond à quelques mailles connexes (quelques mailles toutes reliées entre elles), et de préférence à une seule grille, de manière à simplifier encore le tri. Grâce à toutes les éoliennes implantées sur la portion de grille (par exemple, à quatre éoliennes si la portion de grille comprend une seule maille, 9 éoliennes si la portion de grille est homothétique de facteur 2 par rapport à une grille d’une seule maille), on peut déterminer la production d’énergie annuelle de la mini-ferme, notamment au moyen des formules présentées précédemment. Ainsi, pour chaque type de grille, on peut associer une production d’énergie annuelle de la mini-ferme. Pour chaque grille de l’espace prédéterminé définie à l’étape a), on détermine la production d’énergie annuelle moyenne de la mini-ferme (c’est-à-dire composée d’éoliennes positionnées à l’ensemble des points d’intersection d’une maille ou plusieurs mailles connexes de la grille - de préférence, la mini-ferme est composée de quatre éoliennes, chacune étant positionnée à un point d’intersection différent d’une même maille), à partir de la première répartition discrète de vitesse de vent, de ladite deuxième répartition discrète de direction de vent et de ladite probabilité d’occurrence. Etape c) :Conservation de certaines grilles maximisant la production d’énergie annuelle de la mini-ferme Cette étape sert à limiter le nombre de combinaisons à utiliser pour la disposition du nombre prédéfini des éoliennes de la ferme à partir des résultats obtenus à l’étape b). En effet, grâce à la mini-ferme éolienne, on peut estimer pour différents (premiers) couples de premier et deuxième espacements prédéfinis, la production d’énergie annuelle en fonction du (deuxième) couple de première et deuxième direction prédéfinies et des première et deuxième répartitions discrètes de vitesse de vent et de direction de vent et la probabilité d’occurrence de chaque valeur discrète de vitesse de vent dans chaque valeur discrète de direction du vent. On peut alors choisir (identifier) au moins un (de préférence un seul) (deuxième) couple de première et deuxième direction prédéfinies, de manière à limiter le nombre de combinaisons pour la suite. En effet, le choix du (ou des) (deuxième) couple de première et deuxième direction prédéfinies est déterminé de manière à maximiser l’énergie annuelle produite et/ou de limiter les effets de sillage. Ainsi, les deuxièmes couples définissent les couples de première et deuxième directions les plus prometteuses (pour l’énergie annuelle produite). Pour déterminer le (ou les) (deuxième) couple de première et deuxième direction identifié(s), on peut soit comparer directement la production annuelle produite et retenir les couples de directions qui maximisent cette production annuelle de la mini-ferme, soit comparer la perte d’énergie récupérée due aux effets de sillage. Dans ce deuxième cas, on détermine la perte d’énergie de la mini-ferme considérée par rapport à la production idéale du même nombre d’éoliennes que la mini-ferme, sans tenir compte des effets de sillage des éoliennes entre elles. Ainsi, à l’étape c), on peut par exemple faire la différence entre la production d’énergie annuelle de la mini-ferme obtenue à l’étape b) et la production d’énergie annuelle d’une seule éolienne multipliée par le nombre d’éoliennes de la mini-ferme. Cette perte d’énergie Pl oss peut être déterminée par la formule suivante : n*aep(f1, ws, wp)- aep (f2, ws, wp) Où aep(f1, ws, wp) est l’énergie annuelle produite selon la formule [Math 1] pour une ferme f1 constituée d’une seule éolienne dans l’espace prédéterminé, et où aep (f2, ws, wp) correspond à l’énergie annuelle produite selon la formule [Math 1] pour une ferme f2 constituée d’un nombre n d’éoliennes dans la mini-ferme de l’espace prédéterminé, ws et wp étant les répartitions statistiques des vitesses de vent et de direction de vent, en tenant compte des probabilités d’occurrence préalablement évoquées. L’avantage de cette formule [Math 4] est qu’elle renseigne directement sur l’impact des effets de sillage de la mini-ferme, effet qui dépend directement des directions d’alignement et des espacements choisis. On peut alors écarter(éliminer) les solutions pour lesquelles les pertes sont trop importantes pour ne retenir que la solution ou les solutions pertinentes (avec les pertes d’énergie les plus faibles). Pour chaque premier couple de premier et deuxième espacements prédéfinis, on identifie (choisit) ainsi au moins un deuxième couple de première et deuxième directions prédéfinies qui maximise la production d’énergie annuelle de ladite mini-ferme. En d’autres termes, on sélectionne, un ou des deuxième(s) couple(s) permettant de maximisant la production d’énergie annuelle produite. De ce fait, on limite les étapes suivantes à une sélection de quelques grilles, au lieu de réaliser les étapes suivantes pour toutes les grilles de l’étape a). On peut ensuite former un ensemble qui regroupe tous les deuxièmes couples identifiés qui seront utilisés pour la suite : en d’autres termes, les deuxièmes couples qui ne sont pas conservés dans ledit ensemble ne seront plus considérés pour la suite : on ne conserve ainsi que les deuxièmes couples de première et deuxième directions les plus prometteurs. Cela permet de limiter les combinaisons à réaliser pour la suite. Avantageusement, à l’étape c), on peut identifier un nombre prédéterminé de couple de première et deuxième directions prédéfinies pour chaque couple de premier et deuxième espacements prédéfinis. Ainsi, on peut limiter le nombre de configurations de grilles retenues pour la suite. De préférence, le nombre prédéterminé peut être compris entre 1 et 20, de manière préférée, entre 1 et 10. Ainsi, on offre un bon compromis entre temps de calcul, utilisation de ressources informatiques et précision des résultats de l’optimisation. Etape d) :Pour chaque grille, détermination d’une première disposition des éoliennes On détermine, lors de cette étape, une première disposition des éoliennes sur chaque grille (grilles retenues à l’étape c) ou deuxièmes formées à l’étape e) ou troisièmes grilles formées à l’étape f) selon la réalisation de cette étape d) dans le procédé), de manière à limiter le nombre de combinaisons. Chaque éolienne est positionnée sur un point d’intersection (une seule éolienne sur un même point d’intersection pour des raisons évidentes de construction) de l’espace prédéterminé par un premier algorithme de positionnement. De préférence, ce premier algorithme de positionnement peut être un algorithme d’optimisation qui permet d’obtenir une première répartition permettant d’obtenir une énergie annuelle produite satisfaisante. Avantageusement, ce premier algorithme de positionnement peut être un algorithme glouton qui positionne chaque éolienne l’une après l’autre de manière à maximiser l’énergie annuelle produite de chaque éolienne que l’on vient ajouter. La première éolienne peut être positionnée arbitrairement dans l’espace prédéterminé (à une valeur discrète du premier maillage discret). La deuxième éolienne sera positionnée au point d’intersection de la grille considérée permettant de maximiser l’énergie annuelle produite des deux éoliennes, la position choisie la n-ième éolienne au point d’intersection de la grille permettant de maximiser l’énergie annuelle produite des n éoliennes. L’utilisation d’un algorithme glouton permet d’obtenir, de manière simple, une première disposition des éoliennes dans l’espace prédéterminé, ce qui permet d’initialiser le procédé d’optimisation de l’étape e). Selon un mode de réalisation, pour chaque grille considérée, ledit premier algorithme de positionnement peut réaliser au moins les étapes suivantes : - on définit arbitrairement la position de la première éolienne sur un des points d’intersection de la grille considérée ; - puis pour chaque éolienne à positionner, de manière successive : on définit des positions potentielles pour l’éolienne à positionner, les positions potentielles étant constituées par les points d’intersection (libres) de la grille. En d’autres termes, on sélectionne des points d’intersection de la grille où serait avantageusement positionnée la prochaine éolienne. on calcule la production d’énergie annuelle des éoliennes positionnées et de l’éolienne à positionner pour les positions potentielles définies. De ce fait, pour chaque position potentielle définie, est associée une valeur de production d’énergie annuelle. De plus, le calcul de la production d’énergie annuelle prend en compte la première répartition discrète de vitesse de vent, la deuxième répartition discrète de direction de vent et la probabilité d’occurrence. Ce calcul fait également intervenir de manière connue les caractéristiques des éoliennes, à savoir notamment la surface balayée par le rotor de l’éolienne, le coefficient de traînée et/ou le coefficient de puissance. On choisit la position de l’éolienne à positionner correspondant à la valeur maximale calculée de production d’énergie annuelle. Ainsi, on maximise l’énergie annuelle produite par les éoliennes dont la position est définie dans l’espace prédéterminé. Ces positions définies servent de base pour la détermination de la position de la prochaine éolienne, notamment pour définir les positions potentielles de la prochaine éolienne à positionner. - on détermine la première disposition correspondant à la position du nombre prédéfini d’éoliennes parmi les points d’intersection (libres) de la grille considérée dans l’espace prédéterminé. De ce fait, le premier algorithme de positionnement est un algorithme glouton qui comprend une étape de position arbitraire de la première éolienne, puis qui, itérativement, positionne une éolienne supplémentaire sur des points d’intersection de la grille considérée jusqu’à ce que toutes les éoliennes du nombre prédéfini soient positionnées sur la grille. Cet algorithme glouton permet, grâce à l’étape de sélection des positions potentielles, d’accélérer le temps de calcul et de limiter les ressources informatiques, tout en positionnant les éoliennes de manière judicieuse. Un tel algorithme permet d’obtenir une première disposition adaptée pour l’étape suivante d’optimisation locale de recherche de positionnement de chaque éolienne. La FIG. 4 illustre, de manière schématique et non limitative, un exemple d’une étape pour obtenir une première disposition des éoliennes sur chaque grille sélectionnée selon l’invention. Pour cela, on peut par exemple utiliser un algorithme glouton. Pour chaque grille G(C1, C2) sélectionnée de l’espace prédéterminé, on définit la position de la première éolienne P_E1, par exemple de manière arbitraire. Puis pour chaque éolienne j, on cherche une position de manière à maximiser l’énergie annuelle produite. Ainsi, de manière itérative selon F3, on définit, pour chaque éolienne j, une par une, des positions potentielles PE_Ej pour l’éolienne j à positionner, les positions potentielles PE_Ej étant déterminées par les points d’intersection libres (c’est-à-dire ceux sur lesquels aucune éolienne n’est prévue pour le moment) de la grille de l’espace prédéterminé. Une fois ces positions potentielles PE_Ej définies pour l’éolienne j à positionner, on évalue Eval_AEP, pour chacune de ces positions potentielles PE_Ej, l’énergie annuelle produite des éoliennes déjà positionnées et de l’éolienne j à positionner dans l’espace prédéterminé. Pour cette évaluation Eval_AEP, on utilise notamment les première et deuxièmes répartitions discrètes RD1 et RD2 des vitesses et directions du vent ainsi que la probabilité d’occurrence Prob de chaque vitesse du vent dans chaque direction du vent. On peut, de manière connue, utiliser les caractéristiques des éoliennes et les effets de sillage préalablement décrits dans cette description. On peut alors choisir la position Pos_j de l’éolienne j, position correspondant à la valeur maximale de l’énergie annuelle produite à l’étape précédente. On peut alors définir une nouvelle disposition Disp_Ej des éoliennes positionnées (y compris l’éolienne j) dans l’espace prédéterminé. Cette nouvelle disposition Disp_Ej sert à l’itération suivante pour déterminer les positions potentielles PE_Ej de la nouvelle éolienne j à positionner et pour l’évaluation de l’énergie annuelle produite Eval_AEP. On effectue la boucle F3 pour, j allant de 1 à N-1, N étant le nombre prédéfini d’éoliennes dans l’espace prédéterminé, compte tenu du fait que la première éolienne est positionnée de manière étape à l’étape P_E1. Une fois toutes les éoliennes positionnées (c’est-à-dire lorsque j=N-1), la dernière disposition trouvée Disp_Ej pour chaque grille correspond alors à la première disposition Disp1 des éoliennes pour chaque grille. Ainsi, pour chaque grille définie (qui peut être les grilles retenues l’étape c), les deuxièmes grilles définies à l’étape e) ou les troisièmes grilles définies à l’étape f)), on détermine une première disposition du nombre prédéfini d’éoliennes, chaque éolienne étant positionnée sur un des points d’intersection de la grille concernée par un premier algorithme de positionnement. Etape e) : Modification des positions des éoliennes Puis, pour chaque grille utilisée à l’étape d) précédente, on modifie la position des éoliennes parmi les points d’intersection libres de la grille (les points d’intersection « libres » étant les points d’intersection sur lesquels aucune éolienne n’est positionnée ; en effet, il n’est pas possible de positionner deux éoliennes à la même position), afin d’améliorer (de préférence strictement - par cela, on entend qu’on repositionne l’éolienne uniquement si la production d’énergie annuelle est strictement supérieure à la production d’énergie annuelle de la disposition antérieure) la production d’énergie annuelle de manière à obtenir une disposition du nombre prédéfini d’éoliennes sur chaque grille considérée et à obtenir une production d’énergie annuelle pour chaque disposition. Puis, on retient une quantité prédéterminée de dispositions correspondant aux dispositions maximisant la production d’énergie annuelle et on retient chaque grille associée à chaque disposition retenue. Pour chaque grille de l’étape d), on cherche ainsi à améliorer l’énergie annuelle en modifiant la position des éoliennes sur la grille. Cela peut notamment se faire par les étapes e1) et e2) qui sont décrites par la suite. En effet, on peut réaliser au moins les étapes suivantes : e1) pour chaque grille utilisée à l’étape d), on définit un ordre séquentiel de modification des positions des éoliennes déterminées par le premier algorithme de positionnement et on repositionne, une à une, chaque éolienne de l’ordre séquentiel défini, en trouvant un point d’intersection libre (c’est-à-dire un point d’intersection qui n’est pas occupée par une autre éolienne) de la grille qui maximise la production d’énergie annuelle à partir de ladite première répartition discrète de vitesse de vent, de ladite deuxième répartition discrète de direction de vent et de ladite probabilité d’occurrence ; e2) on réitère l’étape e1) autant de fois que nécessaire jusqu’à ce que, lors à la fin d’une étape e1), aucune éolienne n’ait été repositionnée et on obtient une disposition finale du nombre prédéfini d’éoliennes sur chaque grille considérée et une production d’énergie annuelle pour chaque disposition finale. Etape e1) d’optimisation locale de recherche Lors de cette étape, on cherche à améliorer la production d’énergie annuelle en modifiant la position des éoliennes de la ferme, une par une, sur chaque grille de l’étape d). Pour cela, on définit un ordre séquentiel de modification des positions des éoliennes déterminées par le premier algorithme de positionnement : en d’autres termes, on ne repositionne pas nécessairement les éoliennes, dans le même ordre que celui dans lequel elles ont été déterminées par le premier algorithme et on utilise préférentiellement un ordre différent pour améliorer les chances d’obtenir une meilleure solution. De préférence, cet ordre séquentiel peut être obtenu de manière aléatoire. En effet, en utilisant une fonction aléatoire pour l’ordre séquentiel, on améliore la qualité de l’optimisation en évitant des chemins d’optimisation basés sur des ordres préétablis, ces chemins pouvant biaiser les résultats d’optimisation. Avantageusement, on peut modifier l’ordre séquentiel à chaque itération de l’étape e1), de manière à encore améliorer la qualité de l’optimisation. Une fois l’ordre séquentiel défini, on repositionne, une par une, chaque éolienne, par cet ordre séquentiel défini, en trouvant un point d’intersection libre de la grille qui maximise la production d’énergie annuelle : si aucun autre point d’intersection de la grille ne permet d’améliorer la production d’énergie annuelle, l’éolienne considérée est maintenue à sa position déterminée précédemment. On peut alors tenter de modifier la position de l’éolienne suivante dans l’ordre séquentiel défini. Pour déterminer l’énergie annuelle produite, on tient compte bien-entendu des première et deuxième répartitions discrètes de vitesse de vent et de direction de vent et de la probabilité d’occurrence, qui sont directement dépendantes de la position géographique de l’espace prédéterminé et de son environnement local (présence de forêts, de montagnes, de particularités géologiques, d’immeubles par exemple) et des caractéristiques des éoliennes. L’énergie annuelle produite peut notamment être déterminée à partir des formules précédemment évoquées, notamment par la formule [Math 1]. Selon une variante avantageuse du procédé de l’invention, à l’étape e1), on peut choisir une sélection de points d’intersections libres (c’est-à-dire parmi les points d’intersection sur lesquels aucune éolienne n’est positionnée), pour chaque éolienne à repositionner. En d’autres termes, on limite le nombre de points d’intersection à considérer pour l’éolienne à repositionner de manière à limiter les ressources informatiques. Par exemple, la sélection de points d’intersection peut correspondre à une valeur prédéterminée de points d’intersections libres (sur lesquels aucune éolienne n’est positionnée), ladite sélection de points d’intersection étant choisie de manière aléatoire ou étant choisie comme correspondant à ladite valeur prédéterminée de points d’intersection libres (sur lesquels aucune éolienne n’est positionnée) les plus proches de la position de l’éolienne à repositionner. La sélection aléatoire permet d’améliorer la qualité de l’optimisation, en limitant le nombre de configurations possibles à chaque calcul itératif. La sélection des points d’intersection libres les plus proches de la position de l’éolienne à repositionner permet de délimiter une zone dans laquelle on souhaite repositionner l’éolienne. Ainsi, on délimite peu à peu une zone de plus en plus restreinte pour repositionner l’éolienne. Etape e2) de reproduction de l’étape e1) A la fin de l’étape e1), on ne peut a priori pas savoir si la disposition des éoliennes dans chaque grille est optimale ou si une amélioration peut encore être apportée. Pour encore optimiser la disposition des éoliennes sur chaque grille, on réitère l’étape e1) autant de fois que nécessaire tant que, lors de cette étape, on repositionne au moins une éolienne sur un point d’intersection de la grille. On arrête d’itérer lorsqu’à la dernière itération de l’étape e1), aucune éolienne n’a été repositionnée : on considère alors qu’a été trouvée la disposition optimale du nombre prédéfini d’éoliennes sur la grille considérée. En sortie de l’étape e2), sur chaque grille considérée (issue de l’étape c) ou deuxième grille formée à l’étape f) ou troisième grille formée à l’étape g)), on obtient une disposition des éoliennes sur chaque grille considérée et une production d’énergie annuelle associée à chacune de ces dispositions (une pour chaque grille considérée). On peut alors retenir pour la suite une quantité prédéterminée de dispositions correspondant aux dispositions maximisant la production d’énergie annuelle et on peut retenir chaque grille associée à chaque disposition retenue. Avantageusement, à l’étape e), et plus particulièrement à l’étape e2), la quantité prédéterminée de dispositions peut être comprise entre 1 et 20, de préférence entre 1 et 10. De préférence, à l’étape e2), on peut modifier l’ordre séquentiel à chaque itération de l’étape e1) de manière à limiter les chemins d’optimisation basés sur des ordres préétablis. La FIG. 2 illustre, de manière schématique et non limitative, un exemple de recherche d’optimisation de la méthode selon l’invention. Pour chaque grille sélectionnée, on détermine Alg1 une première disposition Disp1 des éoliennes sur chacune de ces grilles, en tenant compte comme données d’entrée d’au moins une première répartition discrète RD1 de vitesses de vent, une deuxième répartition discrète RD2 de directions du vent, et la probabilité d’occurrences Prob de chaque vitesse de vent dans chaque direction de vent. Les caractéristiques des éoliennes et des modèles de sillage peuvent également être utilisés pour la détermination de l’énergie annuelle produite afin d’établir cette première disposition Disp1. Les données de vent, vitesse, direction et probabilité d’occurrence de chaque valeur de vitesse dans chaque direction peuvent notamment être obtenues à partir de moyens de collecte lors d’une étape préliminaire. Ces données peuvent notamment servir à établir une rose des vents connue de l’homme du métier. Cette première disposition Disp1 est perfectible mais de suffisamment bonne qualité pour permettre les optimisations locales des étapes ultérieures. Cette étape de détermination Alg1 de la première disposition Disp1 des éoliennes sur chaque grille sélectionnée peut être un algorithme glouton. Cette première disposition Disp1 obtenue rapidement grâce à l’algorithme glouton est utilisée comme donnée d’entrée pour l’étape suivant d’optimisation Alg2. Cette étape d’optimisation Alg2 modifie, une à une, la position des différentes éoliennes de la première disposition Disp1, sur d’autres points d’intersection de la même grille, de manière à augmenter l’énergie annuelle produite par la ferme en testant d’autres positions possibles sur la grille, en respectant ainsi les contraintes d’alignement. De manière plus détaillée, cette étape d’optimisation Alg2 comprend les étapes suivantes : - on définit un ordre séquentiel OS de modification des positions des éoliennes, une par une. Cet ordre séquentiel OS peut notamment être obtenu par une fonction aléatoire. - Puis de manière itérative, pour chaque grille sélectionnée, on modifie la position d’au moins une éolienne i, de préférence de chaque éolienne i, en réalisant les sous-étapes suivantes : * on détermine des positions possibles PDP_i pour l’éolienne i, ces positions possibles étant la position actuelle de l’éolienne i et les points d’intersection de la grille sur laquelle aucune éolienne n’est prévue (les autres éoliennes restant à leur position, soit la position initiale issue de la première disposition, soit la position déjà repositionnée). * on évalue Eval_i ensuite l’énergie annuelle produite pour chaque disposition possible (pour chaque position possible PDP_i de l’éolienne i à repositionner, les autres éoliennes restant à la dernière position qui leur a été fixée). L’énergie annuelle produite prend en compte les premières et deuxièmes répartitions discrètes RD1 et RD2 des vitesses et directions du vent, ainsi que la probabilité d’occurrence Prob de chaque vitesse de vent dans chaque direction du vent. Ainsi, à la fin de cette étape d’évaluation, une énergie annuelle produite correspond à chaque position possible PDP_i de l’éolienne i à repositionner. * on retient comme position Pos_i de l’éolienne i la position possible PDP_i qui correspond à la valeur maximale de l’énergie annuelle produite à l’étape Eval_i. * On obtient ainsi une nouvelle disposition Disp_N des éoliennes dans l’espace déterminé. Cette nouvelle disposition comprend les dernières positions des éoliennes déjà positionnées ainsi que la nouvelle position Pos_i de l’éolienne i. La boucle B1 permet ensuite de sélectionner F1 l’éolienne suivante de l’ordre séquentiel (i devient i+1) défini afin d’effectuer la même procédure pour les éoliennes suivantes. Une fois que toutes les éoliennes ont été repositionnées par la boucle B1, on réitère plusieurs fois, pour chaque grille sélectionnée, par une boucle B2, les étapes de détermination de l’ordre séquentiel OS, et de la boucle B1 pour chacune des éoliennes, une par une, la boucle B1 comprenant pour chaque éolienne i à repositionner la détermination des positions possibles PDP_i, l’évaluation de l’énergie annuelle produite Eval_i pour chaque position possible, le choix de la position Pos_i de l’éolienne i à repositionner et la définition de la nouvelle disposition Disp_N. En réitérant plusieurs fois ces étapes, on peut améliorer l’énergie annuelle produite par la ferme. Modifier à chaque fois l’ordre séquentiel (par un tirage aléatoire à chaque fois par exemple) permet encore d’améliorer l’énergie annuelle produite. La boucle B2 se termine lorsqu’au cours de la dernière boucle B1 réalisée, aucune éolienne n’a été repositionnée (toutes les éoliennes ont été maintenues à la même position que celle obtenue à l’itération précédente). Une fois la boucle B2 terminée, la dernière disposition Disp_N obtenue devient la disposition (finale) DispM pour la grille sélectionnée. Parmi les différentes grilles, on peut alors retenir les grilles permettant de maximiser la production d’énergie annuelle (et on élimine les autres grilles qui sont les moins prometteuses). Etape f) : Création de nouvelles grilles et réitération des étapes d) et e) Pour chaque grille retenue à l’étape e), on diminue le premier pas d’espacement, de préférence de moitié, et on forme des deuxièmes grilles à partir du premier pas d’espacement diminué. Chaque deuxième grille est formée à partir d’une grille retenue à l’étape e). La deuxième grille utilise le même deuxième couple de première et deuxième directions que ladite grille retenue mais chaque deuxième grille se distingue de la grille retenue à l’étape e) par les premiers et deuxièmes espacements. En effet, on fait varier les premiers et deuxièmes espacements des deuxièmes grilles des premiers et deuxièmes espacements de la grille retenue à l’étape e), en fonction du premier pas d’espacement diminué. En d’autres termes, ces deuxièmes grilles servent à affiner la discrétisation du pas d’espacement autour de la configuration qu’on a retenue à l’étape précédente (c’est-à-dire de la grille retenue à l’étape e) permettant de maximiser la production d’énergie annuelle), dans le but d’améliorer la précision du procédé. Par exemple, les premiers espacements r 1 ’ des deuxièmes grilles peuvent être définis par : Avec r 1 le premier espacement de la grille retenue à l’étape e) concernée i1 un premier nombre entier positif (qui peut être positif ou négatif) Δr le premier pas d’espacement diminué (par rapport au premier pas d’espacement précédent). Avantageusement, de manière à affiner peu à peu la discrétisation, le produit de la valeur absolue du premier nombre entier et du premier pas d’espacement diminué peut être inférieur (notamment strictement inférieur) au premier pas d’espacement précédent (c’est-à-dire celui utilisé avant la définition du premier pas d’espacement diminué). De la même manière, les deuxièmes espacements r 2 ’ des deuxièmes grilles peuvent être définis par : Avec r 2 le deuxième espacement de la grille retenue à l’étape e) concernée i1 le premier nombre entier positif (qui peut être positif ou négatif) Δr le premier pas d’espacement diminué (par rapport au premier pas d’espacement précédent). Grâce aux deuxièmes grilles, on fait varier l’espacement des grilles de manière à raffiner l’espacement du maillage autour des configurations pertinentes trouvées à l’étape précédente (c’est-à-dire des grilles permettant de maximiser la production d’énergie annuelle). A partir des deuxièmes grilles ainsi formées, on réitère les étapes d) et e) de manière à repositionner les éoliennes, de manière optimale sur chacune des deuxièmes grilles (pour ces réitérations, on remplace les grilles par des deuxièmes grilles formées). De préférence, à l’étape f), on peut former un nombre prédéterminé de deuxièmes grilles. Chaque deuxième grille peut être définie par des lignes de même orientations que la grille retenue à l’étape e), le premier espacement de chaque deuxième grille peut correspondre à la somme du premier espacement de la grille retenue à l’étape e) et d’un nombre entier de premier pas d’espacement diminué, et le deuxième espacement de chaque deuxième grille peut correspondre à la somme du deuxième espacement de la grille retenue à l’étape e) et dudit nombre entier de premier pas d’espacement diminué. Ainsi, pour chaque grille retenue, on forme un nombre prédéterminé de deuxièmes grilles, le nombre prédéterminé dépendant du nombre entier. Les deuxièmes grilles permettent d’affiner les espacements autour des premier et deuxième espacements des grilles retenues à l’étape e). Préférentiellement, le nombre entier peut être 1 ou -1. En considérant ces valeurs de 1 et -1, pour chaque grille retenue à l’étape e), on forme quatre deuxièmes grilles, ce qui permet de limiter le nombre de configurations à tester, tout en affinant la précision des grilles peu à peu autour des grilles précédentes les plus prometteuses (celles qui permettent de maximiser la production d’énergie annuelle). Etape g) : Création de nouvelles grilles et réitération des étapes d) et e) Pour chaque grille retenue à l’étape f) (c’est-à-dire à la fin des réitérations des étapes d) et e)), on diminue le deuxième pas de direction, de préférence de moitié, et on forme des troisièmes grilles à partir du deuxième pas de direction. Chaque troisième grille est formée à partir d’une grille retenue à l’étape f) (à la fin de l’étape f), c’est-à-dire après l’étape de repositionnement local de chaque éolienne). La troisième grille utilise le même premier couple de premier et deuxième espacements que ladite grille retenue mais chaque troisième grille se distingue de la grille retenue à l’étape f) par les premières et deuxièmes directions. En effet, on fait varier les premières et deuxièmes directions des troisièmes grilles des premières et deuxièmes directions de la grille retenue à l’étape f), en fonction du deuxième pas de direction diminué. En d’autres termes, ces troisièmes grilles servent à affiner la discrétisation du deuxième pas de direction autour de la configuration qu’on a retenue à l’étape précédente (c’est-à-dire de la grille retenue à l’étape f) permettant de maximiser la production d’énergie annuelle), dans le but d’améliorer la précision du procédé. Par exemple, les premières direction θ 1 ’ des troisièmes grilles peuvent être définies par : Avec θ 1 la première direction de la grille retenue à l’étape f) concernée i2 un deuxième nombre entier positif (qui peut être égal au premier nombre entier) Δθ le deuxième pas de direction diminué (par rapport au deuxième pas de direction précédent). Avantageusement, de manière à affiner peu à peu la discrétisation, le produit de la valeur absolue du deuxième nombre entier et du deuxième pas de direction diminué peut être inférieur (notamment strictement inférieur) au deuxième pas de direction précédent (c’est-à-dire celui utilisé avant la définition du deuxième pas de direction diminué). De la même manière, les deuxièmes directions θ 2 ’ des troisièmes grilles peuvent être définies par : Avec θ 2 la deuxième direction de la grille retenue à l’étape f) concernée i2 le deuxième nombre entier positif Δθ le deuxième pas de direction diminué (par rapport au deuxième pas de direction précédent). Grâce aux troisièmes grilles, on fait varier la direction des lignes des grilles de manière à raffiner les variations de direction du maillage autour des configurations pertinentes trouvées à l’étape précédente (c’est-à-dire des grilles permettant de maximiser la production d’énergie annuelle). A partir des troisièmes grilles ainsi formées, on réitère les étapes d) et e) de manière à repositionner les éoliennes, de manière optimale sur chacune des deuxièmes grilles (pour ces réitérations, on remplace les grilles par des deuxièmes grilles formées). De préférence, à l’étape g), on peut former un nombre prédéterminé de troisièmes grilles. Chaque troisième grille peut être définie par des lignes de même espacement que la grille retenue à l’étape f), la première direction de chaque troisième grille peut correspondre à la somme de la première direction de la grille retenue à l’étape f) et d’un nombre entier du deuxième pas de direction diminué, et la deuxième direction de chaque troisième grille peut correspondre à la somme de la deuxième direction de la grille retenue à l’étape f) et dudit nombre entier de deuxième pas de direction diminué. Ainsi, pour chaque grille retenue, on forme un nombre prédéterminé de troisièmes grilles, le nombre prédéterminé dépendant du nombre entier. Les troisièmes grilles permettent d’affiner les directions autour des première et deuxième directions des grilles retenues à l’étape f). Préférentiellement, le nombre entier peut être 1 ou -1. En considérant ces valeurs de 1 et -1 , pour chaque grille retenue à l’étape f), on forme quatre troisièmes grilles, ce qui permet de limiter le nombre de configurations à tester, tout en affinant la précision des grilles peu à peu autour des grilles précédentes les plus prometteuses (celles qui permettent de maximiser la production d’énergie annuelle). Grâce aux étapes f) et g), on peut raffiner, en deux temps, le maillage, en termes d’espacements et de directions. Ainsi, à l’étape a), on utilise un maillage grossier de manière à présélectionner certaines grilles à l’étape c). On resélectionne ensuite certaines grilles à l’étape e). On peut alors raffiner le maillage sur les grilles retenues à l’étape e) les plus prometteuses. Ainsi, le procédé permet de limiter les ressources informatiques nécessaires (mémoire et processeur), de limiter le temps de calcul, tout en affinant la précision de la position des différentes éoliennes en maximisant l’énergie annuelle produite. Etape h) : Détermination de la position définitive des éoliennes dans l’espace prédéterminé On détermine une disposition définitive des éoliennes dans l’espace prédéterminé. La disposition définitive correspond à la disposition (finale) obtenue à l’étape g) de la grille qui maximise la production d’énergie annuelle. Cette disposition définitive peut être utilisée pour la construction des éoliennes selon la disposition définitive dans l’espace prédéterminé, cette disposition définitive permettant de maximiser l’énergie annuelle de la ferme. Pour chaque disposition (finale) obtenue à la fin de l’étape g) des éoliennes obtenues sur chaque grille (qui sont donc des troisièmes grilles), on compare la production d’énergie annuelle de chacune de ces dispositions (finales) et on choisit la disposition (finale) qui a la production d’énergie annuelle la plus élevée pour maximiser le rendement de la ferme. On détermine alors la disposition définitive des éoliennes de la ferme comme correspondant à cette disposition (finale) optimale. A partir de la disposition définitive obtenue, on peut alors construire les éoliennes de la ferme aux positions correspondant à la disposition définitive dans l’espace prédéterminé. Grâce à l’utilisation de grilles, les éoliennes seront alignées selon deux directions et selon des lignes régulièrement réparties dans l’espace, de manière à répondre aux contraintes opérationnelles et aux contraintes spécifiques à l’offshore. Dans le présent procédé, les étapes a) à h) sont mises en œuvre par des moyens informatiques, tels qu’un ordinateur, un serveur et/ou un supercalculateur. Les étapes a) à h) constituent également un procédé de positionnement des éoliennes dans un espace prédéterminé qui est prévu pour construire les éoliennes de la ferme aux positions définies par le procédé de positionnement. De manière avantageuse, on peut réitérer les étapes f) et/ou g) jusqu’à ce que le premier pas d’espacement atteigne un premier seuil et/ou jusqu’à ce que le deuxième pas de direction atteigne un deuxième seuil. En réitérant les étapes f) et/ou g), on augmente la précision de positionnement et on peut augmenter l’énergie annuelle produite en augmentant la discrétisation des espacements et des directions des grilles. Selon une configuration de l’invention, le premier seuil peut être compris entre un cinquantième du diamètre des éoliennes et un cinquième du diamètre des éoliennes. Ainsi, la discrétisation est suffisamment fine pour obtenir une production d’énergie annuelle élevée et une précision sur les positions des éoliennes. Alternativement ou additionnellement, le premier seuil peut être ou inférieur à 1m. En effet, la précision de la construction des éoliennes aux positions prévues d’un mètre est suffisante. Selon une configuration de l’invention, le deuxième seuil peut être inférieur à 2°, de préférence inférieur à 1°. En effet, cette précision sur les directions des lignes des grilles est suffisante. L’invention concerne aussi un procédé de positionnement des éoliennes d’une ferme éolienne dans un espace prédéterminé pour la construction de la ferme éolienne dans l’espace prédéterminé. Dans ce procédé, on mesure des données statistiques du vent, par des moyens de mesure, de préférence un capteur LIDAR, dans l’espace prédéterminé pour déterminer les première et deuxième répartitions discrètes et les probabilités d’occurrence de chaque vitesse de vent dans chaque direction du vent des première et deuxième répartitions discrètes, puis on met en œuvre le procédé de simulation selon l’une des variantes ou combinaisons de variantes décrites précédemment, pour déterminer les positions optimales des éoliennes de la ferme éolienne dans l’espace prédéterminé, de manière à positionner les éoliennes dans l’espace prédéterminé pour la construction de la ferme éolienne dans l’espace prédéterminé. Grâce à ce procédé, on peut augmenter la précision de l’énergie annuelle produite, basée sur des données statistiques mesurées à l’emplacement de l’espace prédéterminé. On peut alors déterminer les positions optimales des éoliennes dans cet espace, en tenant compte des contraintes d’alignement et en maximisant l’énergie annuelle produite, à partir des données de vent mesurées. L’invention peut également concerner un produit programme d’ordinateur mettant en œuvre le procédé de simulation (procédé de positionnement) des éoliennes dans un espace prédéterminé, constitué des étapes a) à h) décrites précédemment à partir de moyens informatiques, comme un ordinateur, un téléphone portable ou une tablette. Le produit programme d'ordinateur peut être téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur ou un serveur, et il comprend des instructions de code de programme pour la mise en œuvre du procédé de simulation (ou procédé de positionnement) selon l’une des variantes ou combinaisons de variantes décrites précédemment, lorsque le programme est exécuté sur un ordinateur ou sur un téléphone portable. En effet, le procédé de simulation (ou de positionnement) décrit précédemment est particulièrement adapté pour être implémenté sur des moyens informatiques. Il peut ainsi être mis en œuvre de manière simple et des résultats peuvent être obtenus rapidement. En outre, l’invention concerne aussi une méthode de construction d’une ferme éolienne dans laquelle on met en œuvre le procédé de simulation selon l’une des variantes ou combinaisons de variantes décrites précédemment ou dans laquelle on met en œuvre le procédé de positionnement des éoliennes décrit précédemment, puis on construit la ferme éolienne en érigeant les éoliennes aux positions de ladite disposition définitive dans l’espace prédéterminé de manière à générer de l’énergie à partir du vent. En construisant les éoliennes aux positions optimales déterminées par le procédé de simulation ou par le procédé de positionnement, on peut augmenter l’énergie annuelle produite par la ferme dans l’espace prédéterminé. De plus, l’invention concerne également une ferme éolienne obtenue à partir de la méthode de construction d’une ferme éolienne telle que décrite précédemment. La FIG. 6 illustre un exemple d’application du procédé de simulation de ferme éolienne utilisé pour trouver les positions optimales des éoliennes de la ferme dans un espace prédéterminé Z2 qui est un espace non-convexe de manière à utiliser la disposition optimale pour construire les éoliennes à ces positions optimales. Dans cet exemple, on cherche à positionner vingt-neuf éoliennes dans l’espace prédéterminé Z2, en prenant en compte des contraintes d’alignement des éoliennes selon deux directions non colinéaires. Les éoliennes considérées sont des turbines de puissance nominale de 10 MW, dont la hauteur de la nacelle est positionnée à 119,8 m du sol et dont le diamètre du rotor est de 198 m. La distance entre les éoliennes dans l’espace prédéterminé Z2 est supérieure à 4 fois le diamètre du rotor des éoliennes. Les données de vent utilisées pour cet exemple correspondent à celles issues de : Baker, N. F., Thomas, J. J., Stanley, A. P. J., and Ning, A. IEA Task 37 Wind Farm Layout Optimization Case Studies, https://doi.org/10.5281/zenodo.5809681, 2021 Pour cet exemple, on considère les critères suivants : La distance minimale entre les éoliennes (deux à deux) est de deux fois le diamètre des éoliennes (aussi appelé diamètre du rotor) ; La distance maximale entre les éoliennes (deux à deux) est de six fois le diamètre des éoliennes (aussi appelé diamètre du rotor) ; La quantité prédéterminée de dispositions retenue à chaque étape e) est de 10 ; Le nombre prédéterminé de deuxième couple de première et deuxième directions associé à chaque premier couple de premiers et deuxièmes espacements à l’étape c) est de 10. Le premier pas d’espacement utilisé à l’étape a) est de deux fois le diamètre des éoliennes ; puis à chaque étape f), il est divisé par deux ; Le deuxième pas de direction utilisé à l’étape a) est de 5° ; puis à chaque étape g), il est divisé par deux ; Le premier seuil est défini par un dixième du diamètre des éoliennes ; Le deuxième seuil est défini par 0,5° ; A chaque étape f) et g), on forme, pour chaque grille retenue, quatre deuxième ou autre troisième grilles respectivement. La FIG. 6 montre le résultat obtenu, c’est-à-dire la disposition optimale trouvée grâce au procédé de simulation selon l’invention. Sur cette figure : Les petits points en gris clair représentent des positions possibles où on peut placer les éoliennes sur la grille. Ils sont donc aux intersections des premières lignes L1 et des deuxièmes lignes L2 de la grille dans l’espace prédéterminé Z2 ; Les gros points noirs représentent les positions des éoliennes permettant de maximiser la récupération d’énergie sur la grille considérée ; Les lignes pointillées représentent la direction des premières lignes L1 ou des deuxièmes lignes L2. A la fin de l’étape c), si on forme vingt deuxièmes couples regroupant les deuxièmes couples identifiés, pour obtenir la même précision en angle et en position dans l’espace, la méthode selon l’invention réalise 380 optimisations alors que la méthode de la demande de brevet FR 22/09.630 de la demanderesse réalise 33 620 optimisations dans les mêmes conditions. L’utilisation de la méthode permet de contraindre le positionnement avec des alignements des éoliennes, d’une part en limitant le nombre de configurations à tester (donc de limiter les ressources informatiques et le temps de calcul) et d’autre part en améliorant la récupération d’énergie annuelle en affinant la discrétisation des différentes grilles peu à peu. Procédé de simulation d’une ferme éolienne dans un espace prédéterminé (Esp) pour la construction de ladite ferme éolienne dans ledit espace prédéterminé (Esp), la ferme éolienne étant composée d’un nombre prédéfini d’éoliennes, à partir une première répartition discrète de vitesse de vent (RD1), une deuxième répartition discrète de direction de vent (RD2) et une probabilité d’occurrence (Prob) de chaque valeur discrète de vitesse de vent dans chaque valeur discrète de direction du vent desdites première et deuxième répartitions discrètes (RD1, RD2), dans lequel on effectue au moins les étapes successives suivantes par des moyens informatiques, tels qu’un ordinateur : Pour différents premiers espacements prédéfinis séparés d’un premier pas d’espacement, différents deuxièmes espacements prédéfinis séparés dudit premier pas d’espacement, différentes premières directions prédéfinies séparées d’un deuxième pas de direction, et différentes deuxièmes directions prédéfinies séparées dudit deuxième pas de direction, on forme différents premiers couples (C1) de premiers et deuxièmes espacements prédéfinis, et différents deuxièmes couples (C2) de premières et deuxièmes directions prédéfinies (θ1, θ2), puis pour les différents premiers couples (C1) et pour les différents deuxièmes couples (C2), on forme des premières grilles (GR) dans l’espace prédéterminé (Esp), chaque première grille étant définie par une pluralité de points d’intersection entre des premières lignes (L1) orientées dans la première direction prédéfinie (θ1) et des deuxièmes lignes (L2) orientées dans la deuxième direction prédéfinie (θ2), les premières lignes (L1) étant espacées dudit deuxième espacement prédéfini le long de la deuxième direction prédéfinie (θ2) et les deuxièmes lignes (L2) étant espacées dudit premier espacement prédéfini le long de la première direction prédéfinie (θ1), chaque première grille comprenant au moins une maille délimitée par des points d’intersection. Pour chaque première grille dudit espace prédéterminé, on détermine la production d’énergie annuelle moyenne d’une mini-ferme (AEP-mf) composée d’éoliennes positionnées à l’ensemble des points d’intersection d’une maille ou plusieurs mailles connexes de la première grille, à partir de ladite première répartition discrète de vitesse de vent (RD1), de ladite deuxième répartition discrète de direction de vent (RD2) et de ladite probabilité d’occurrence (Prob) ; Pour chaque premier couple (C1) de premier et deuxième espacements prédéfinis, on identifie (Ch) au moins un deuxième couple (C2) de première et deuxième directions prédéfinies (θ1, θ2) qui maximise la production d’énergie annuelle de ladite mini-ferme, puis on forme un ensemble de deuxièmes couples regroupant les deuxièmes couples identifiés ; Pour chaque première grille définie par un premier couple de premier et deuxième espacements prédéfinis et par un deuxième couple de première et deuxième directions prédéfinies parmi ledit ensemble de deuxièmes couples, on détermine une première disposition (Alg1) du nombre prédéfini d’éoliennes, chaque éolienne étant positionnée sur un desdits points d’intersection de la première grille par un premier algorithme de positionnement ; puis, pour chaque première grille utilisée à l’étape d), on modifie la position des éoliennes (Alg2) parmi les points d’intersection de la première grille sur lesquels aucune éolienne n’est positionnée, afin d’améliorer la production d’énergie annuelle de manière à obtenir une disposition du nombre prédéfini d’éoliennes sur chaque première grille considérée et une production d’énergie annuelle pour chaque disposition et on retient (AEP_g) une quantité prédéterminée de dispositions correspondant aux dispositions maximisant la production d’énergie annuelle et chaque première en utilisant les deuxièmes grilles en remplacement des premières grilles e grille (Grr) associée à chaque disposition retenue ; Pour chaque première grille retenue à l’étape e), on diminue le premier pas d’espacement (MΔr), de préférence de moitié, et on forme des deuxièmes grilles (Gr2) à partir du premier pas d’espacement diminué, puis on réitère les étapes d) et e) ; Pour chaque deuxième grille retenue à l’étape f), on diminue le deuxième pas de direction (MΔθ), de préférence de moitié, et on forme des troisièmes grilles (Gr3) à partir du deuxième pas de direction, puis on réitère les étapes d) et e) en utilisant les troisièmes grilles en remplacement des premières grilles ; On détermine une disposition définitive (Disp_F) des éoliennes dans l’espace prédéterminé (Esp), ladite disposition définitive (Disp_F) correspondant à la disposition obtenue à l’étape g) de la troisième grille qui maximise la production d’énergie annuelle, pour la construction (Const) des éoliennes selon la disposition définitive (Disp_F) dans l’espace prédéterminé (Esp). Procédé de simulation selon la revendication 1, dans lequel on réitère les étapes f) et/ou g) jusqu’à ce que le premier pas d’espacement atteigne un premier seuil et/ou jusqu’à ce que le deuxième pas de direction atteigne un deuxième seuil. Procédé de simulation selon la revendication 2, dans lequel le premier seuil est compris entre un cinquantième du diamètre des éoliennes et un cinquième du diamètre des éoliennes ou inférieur à 1 m et/ou le deuxième seuil est inférieur à 2°, de préférence inférieur à 1°. Procédé de simulation selon l’une des revendications précédentes, dans lequel, à l’étape f), on forme un nombre prédéterminé de deuxièmes grilles (Gr2), chaque deuxième grille (Gr2) étant définie par des lignes de même orientations que la première grille retenue à l’étape e), le premier espacement de chaque deuxième grille (Gr2) correspondant à la somme du premier espacement de la première grille retenue à l’étape e) et d’un premier nombre entier de premier pas d’espacement diminué, le deuxième espacement de chaque deuxième grille (Gr2) correspondant à la somme du deuxième espacement de la première grille retenue à l’étape e) et dudit premier nombre entier de premier pas d’espacement diminué, de préférence, ledit premier nombre entier étant 1 ou -1. Procédé de simulation selon l’une des revendications précédentes, dans lequel, à l’étape g), on forme un nombre prédéterminé de troisièmes grilles (Gr3), chaque troisième grille (Gr3) étant définie par des lignes de même espacements que la deuxième grille retenue à l’étape f), la première direction de chaque troisième grille (Gr3) correspondant à la somme de la première direction de la deuxième grille retenue à l’étape f) et d’un deuxième nombre entier du deuxième pas de direction diminué, la deuxième direction de chaque troisième grille (Gr3) correspondant à la somme de la deuxième direction de la deuxième grille retenue à l’étape f) et dudit deuxième nombre entier de deuxième pas de direction diminué, de préférence, ledit deuxième nombre entier étant 1 ou -1. Procédé de simulation selon l’une des revendications précédentes, dans lequel, à l’étape e) de modification de la position des éoliennes (Alg2) pour chaque première grille, on réalise au moins les étapes suivantes : e1) pour chaque première grille utilisée à l’étape d), on définit un ordre séquentiel (OS) de modification des positions des éoliennes déterminées par le premier algorithme de positionnement (Alg1) et on repositionne, une à une, chaque éolienne de l’ordre séquentiel défini, en trouvant un point d’intersection libre de la première grille qui maximise la production d’énergie annuelle à partir de ladite première répartition discrète de vitesse de vent (RD1), de ladite deuxième répartition discrète de direction de vent (RD2) et de ladite probabilité d’occurrence (Prob) ; e2) on réitère l’étape e1) autant de fois que nécessaire jusqu’à ce que, lors à la fin d’une étape e1), aucune éolienne n’ait été repositionnée et on obtient une disposition du nombre prédéfini d’éoliennes sur chaque première grille considérée et une production d’énergie annuelle pour chaque disposition finale ; et on retient (AEP_g) une quantité prédéterminée de dispositions correspondant aux dispositions maximisant la production d’énergie annuelle et chaque première grille (Grr) associée à chaque disposition retenue. Procédé de simulation selon la revendication 6, pour lequel ledit ordre séquentiel (OS) est obtenu de manière aléatoire et/ou on modifie l’ordre séquentiel (OS) à chaque itération de l’étape e1). Procédé de simulation selon l’une des revendications 6 à 7, dans lequel, à l’étape e1), on choisit une sélection d’intersections parmi les points d’intersection sur lesquels aucune éolienne n’est positionnée, pour chaque éolienne à repositionner. Procédé de simulation selon la revendication 8, dans lequel la sélection d’intersections correspond à une valeur prédéterminée de points d’intersections sur lesquels aucune éolienne n’est positionnée, ladite sélection d’intersections étant choisie de manière aléatoire ou étant choisie comme correspondant à ladite valeur prédéterminée de points d’intersections sur lesquels aucune éolienne n’est positionnée les plus proches de la position de l’éolienne à repositionner. Procédé de simulation selon l’une des revendications précédentes, dans lequel ledit espace prédéterminé (Esp) comprend des zones non connexes (Z1, Z2) et/ou non convexes (Z2). Procédé de simulation selon l’une des revendications précédentes, dans lequel pour chaque première grille considérée, ledit premier algorithme de positionnement (Alg1) réalise au moins les étapes suivantes : - on définit arbitrairement la position de la première éolienne (P_E1) sur un des points d’intersection de la première grille ; - puis pour chaque éolienne à positionner, de manière successive : on définit des positions potentielles (PE_Ej) pour l’éolienne à positionner, lesdites positions potentielles (PE_Ej) étant constituées par les points d’intersection de la première grille situés entre une distance minimale et une distance maximale de toutes les éoliennes positionnées ; on calcule la production d’énergie annuelle (Eval_AEP) des éoliennes positionnées et de l’éolienne à positionner pour les positions potentielles (PE_Ej) définies, à partir de ladite première répartition discrète (RD1) de vitesse de vent, de ladite deuxième répartition (RD2) discrète de direction de vent et de ladite probabilité d’occurrence (Prob) ; On choisit la position (Pos_j) de l’éolienne à positionner correspondant à la valeur maximale calculée de production d’énergie annuelle : - on détermine ladite première disposition (Disp_Ej) correspondant à la position du nombre prédéfini d’éoliennes parmi les points d’intersection de la première grille considérée dans ledit espace prédéterminé (Esp). Procédé de simulation selon l’une des revendications précédentes, dans lequel, à l’étape c), on fait la différence entre la production d’énergie annuelle de la mini-ferme obtenue à l’étape b) et la production d’énergie annuelle d’une seule éolienne multipliée par le nombre d’éoliennes de la mini-ferme. Procédé de simulation selon l’une des revendications précédentes, dans lequel, à l’étape a), le premier pas d’espacement est compris entre deux fois et dix fois le diamètre des éoliennes, de préférence entre deux fois et six fois le diamètre des éoliennes. Procédé de simulation selon l’une des revendications précédentes, dans lequel à l’étape a), le deuxième pas de direction est compris entre 1° et 20°, de préférence entre 5° et 10°. Procédé de simulation selon l’une des revendications précédentes, dans lequel, à l’étape e), la quantité prédéterminée de dispositions est comprise entre 1 et 20, de préférence entre 1 et 10. Procédé de simulation selon l’une des revendications précédentes, dans lequel à l’étape c), on identifie un nombre prédéterminé de deuxième couple de première et deuxième directions prédéfinies pour chaque premier couple de premier et deuxième espacements prédéfini, de préférence le nombre prédéterminé est compris entre 1 et 20, de manière préférence entre 1 et 10. Procédé de positionnement des éoliennes d’une ferme éolienne dans un espace prédéterminé pour la construction (Const) de ladite ferme éolienne dans ledit espace prédéterminé (Esp), dans laquelle on mesure des données statistiques du vent, par des moyens de mesure, de préférence un capteur LIDAR, dans l’espace prédéterminé (Esp) pour déterminer les première et deuxième répartitions discrètes (RD1, RD2) et les probabilités d’occurrence (Prob) de chaque vitesse de vent dans chaque direction du vent des première et deuxième répartitions discrètes (RD1, RD2), puis on met en œuvre le procédé de simulation selon l’une des revendications précédentes pour déterminer les positions des éoliennes de ladite ferme éolienne dans l’espace prédéterminé, de manière à positionner les éoliennes dans l’espace prédéterminé pour la construction de ladite ferme éolienne dans ledit espace prédéterminé. Procédé de construction d’une ferme éolienne dans lequel on met en œuvre le procédé de simulation selon l’une des revendications 1 à 16 ou dans lequel on met en œuvre le procédé de positionnement des éoliennes selon la revendication 17, puis on construit (Const) la ferme éolienne en érigeant les éoliennes aux positions de ladite disposition définitive (Disp_F) dans l’espace prédéterminé (Esp). Ferme éolienne obtenue à partir du procédé de construction d’une ferme éolienne selon la revendication précédente.