Die Erfindung betrifft einen Windpark mit mindestens zwei Wind ^¬ energieanlagen, die jeweils einen Generator mit einem Umrichter zur Erzeugung elektrischer Energie aufweisen.

Es ist bei einem Betrieb von Windpark am elektrischen Übertragungsnetz wegen der zunehmenden Anzahl und Größe der Windparks von steigender Bedeutung, dass sie stabilisierend auf das Netz wirken. Die Windparks leisten damit ihren Beitrag zur Sicherung der Stabilität des Stromnetzes. Ein wesentlicher Parameter der Netzqualität ist der Oberschwingungsgehalt. Einschlägige Verord ^¬ nungen begrenzen Amplitude in Abhängigkeit von der Ordnung der Oberschwingungen, um eine ausreichende Netzqualität zu erhalten. Moderne Umrichter gesteuerte Windenergieanlagen erzeugen aufgrund der Schaltvorgänge im Umrichter Oberschwingungen verschiedener Art, sowohl solche mit gerader sowie ungerader Ordnung wie auch solche mit starrer Phasenlage zum Netz oder solche mit va ^¬ riabler Phasenlage zum Netz. Die Windenergieanlagen werden im Rahmen einer für jeden Windenergieanlagentyp vorzunehmenden Typprüfung vermessen und daraus Kennwerte bestimmt. Sind mehrere Windenergieanlagen an das Netz anzuschließen, wie im Fall eines Windparks, werden in der Regel die Oberschwingungen mit Hilfe eines Rechenverfahrens aus den Werten für die einzelnen Wind ^¬ energieanlagen bestimmt. Die vom Windpark insgesamt abgegebenen Oberschwingungen müssen unterhalb gewisser Grenzwerte bleiben. In der Praxis wirkt sich dies oft limitierend auf den Windpark aus .

Bei der Berechnung der Oberschwingungen für Windparks basierend auf den zertifizierten Messwerten für die einzelne Windenergie- anläge wird je nach Art der Oberschwingung unterschiedlich vorgegangen: Oberschwingungen mit niedrigen ungeraden Ordnungszahlen werden als typisch für maschinen- oder netzgeführte Umrichter in Thyristortechnik angesehen. Die Oberschwingungen sind bezogen auf die Grundschwingung phasenstarr und die Oberschwingungen der einzelnen Windenergieanlagen addieren sich skalar auf. Es ergibt sich für einen Park mit N Windenergieanlagen der N- fache Wert einer einzelnen Windenergieanlage. Oberschwingungen höherer Ordnung, insbesondere im Bereich der Umrichterschaltfre ^¬ quenz, werden als typisch für freilaufende Pulswechselrichter angesehen. Die Phasenlage der damit versehenen Windenergieanla ^¬ gen ist nicht fest auf die Grundschwingung bezogen, sondern kann relativ zu ihr driften. Tatsächlich ist die Phasenlage statistisch gleichverteilt. Damit addieren sich die Werte der einzel ^¬ nen Windenergieanlagen nicht direkt auf, sondern entsprechend der statistischen Verteilung kann für diese Oberschwingungen ein Anwachsen entsprechend der Quadratwurzel der Windenergieanlagen angesetzt werden (100 Windenergieanlagen erzeugen damit nicht den 100-fachen, sondern nur den 10-fachen Wert einer einzelnen Windenergieanlage) . Da bezogen auf das Vorgehen bei niedrigen Ordnungszahlen mit der skalaren Vervielfachung sich eine Verringerung ergibt, spricht man auch von einem Reduktionsfaktor mit 1/ (N) .

In den Verordnungen nur unzureichend betrachtet sind Oberschwingungen mit niedrigen, aber geraden Ordnungszahlen. In der Regel enthalten sie einen beträchtlichen phasenstarren Anteil, sind also gemischt phasenstarr und phasenvariabel. Der auf phasenva ^¬ riable Oberschwingungen anwendbare Reduktionsfaktor 1/V(N) kann daher nicht zur Anwendung kommen. Herkömmlicherweise erfolgt da ^¬ her eine direkte, skalare Addition ohne Reduktion. Dies ist un ^¬ befriedigend, da die Grenzwerte häufig recht niedrig liegen wo ^¬ durch die Größe des Windparks hinsichtlich der Zahl der Windenergieanlagen deutlich beschränkt wird. Außerdem macht die Verschiedenheit der Berechnungsverfahren je nach Art der Oberschwingung die Anwendung umständlich.

Der Erfindung liegt die Aufgabe zugrunde, diesen Nachteil zu vermeiden und eine universelle Reduktion zu erreichen, die für alle geraden, ungeraden und auch auf Zwischenfrequenzen liegende Oberschwingungen anwendbar ist, und zwar insbesondere auch bei niedrigen geraden Oberschwingungen.

Die erfindungsgemäße Lösung liegt in den Merkmalen der unabhängigen Ansprüche. Vorteilhafte Weiterbildungen sind Gegenstand der abhängigen Ansprüche.

Bei einem Windpark mit mindestens zwei Windenergieanlagen, die jeweils einen Generator mit einem Umrichter zur Erzeugung elektrischer Energie aufweisen, ist erfindungsgemäß vorgesehen ein Oberschwingungsprädiktor, der die vom Windpark erwartete Oberschwingung bestimmt zur Begrenzung auf einen Oberschwingungsgrenzwert, und der dazu ausgebildet ist mit einem Berechnungsmo ^¬ dul, das einen komplexen Mittelwert über mindestens eine Periode der Oberschwingung einer der Windenergieanlagen berechnet und daraus einen ersten Ersatzvektor bestimmt, einem Iterationsmo ^¬ dul, welches das Berechnungsmodul sukzessive auf mindestens eine andere der Windenergieanlagen aufschaltet zur Bildung mindestens eines zweiten Ersatzvektors, und einem Summiermodul, welches die Ersatzvektoren aufaddiert zu einem Gesamtvektor und mit dem Oberschwingungsgrenzwert vergleicht .

Zuerst seien einige verwendete Begriffe erläutert:

Unter einem komplexen Mittelwert wird ein Mittelwert verstanden, der in der komplexen Zahlenebene unter Berücksichtigung von Re- alteil und Imaginärteil der Oberschwingung berechnet wird. Es handelt sich also um eine komplexe Zahl.

Unter „sukzessive ... aufschaltet" wird verstanden, dass das Ite ^¬ rationsmodul mit dem Berechnungsmodul derart zusammenwirkt, dass der komplexe Mittelwert gebildet wird für eine zweite Windener ^¬ gieanlage, und ggf. für eine dritte, vierte und weitere Wind ^¬ energieanlagen, so dass dementsprechend zweite, dritte, vierte und ggf. weitere Ersatzvektoren gebildet werden.

Die Erfindung fußt auf der Erkenntnis, dass eine genauere Ab ^¬ schätzung der von den Windenergieanlagen eines Windparks insgesamt abgegebenen Oberschwingung ermöglicht ist durch Berücksichtigung der jeweiligen mittleren Phasenlagen. Die Erfindung erreicht dies auf elegante Weise durch Kombination von einer Be ^¬ rechnung für jede Windenergieanlage für sich, wobei die sukzes ^¬ sive in einem Iterationsverfahren für alle zu betrachtenden Windenergieanlagen des Windparks durchgeführt wird, mit einer anschließenden Aufaddierung über die einzelnen Windenergieanlagen. Damit ist die nur schwierig durchführbare Berechnung des Windparks im Gesamten rückgeführt auf eine Berechnung der ein ^¬ zelnen Windenergieanlagen, wobei durch die Ausbildung des Berechnungsmoduls mit Implementation eines komplexen Mittelwert ^¬ bildners eine phasenrichtige Berücksichtigung erreicht werden kann. Windenergieanlagen mit phasenvariablen Teilen der Oberschwingungen, die versetzt zueinander liegen, heben sich damit teilweise auf bzw. addieren sich in nur geringem Maße. Mit dieser vektoriellen Aufaddierung wird schließlich ein Gesamtvektor für den gesamten Windpark bestimmt, der beträchtlich kleiner ist als der herkömmlicherweise berechnete.

Die Erfindung hat erkannt, dass ein auf dieser Grundlage gebil ^¬ deter Prädiktor die Oberschwingungen auch in solchen Fällen hinreichend genau abschätzt, bei denen die Oberschwingungen aus phasenstarren und phasenvariablen Bestandteil bestehen. Herkömmlicherweise konnten nur Fälle mit rein phasenvariablen Oberschwingungen durch vektorielle Addition bestimmt werden. Die Erfindung ermöglicht dies nun auch für Mischfälle, und erreicht damit einen breit einsetzbaren Prädiktor. Er eignet sich für alle Oberschwingungen, seien es gerade, ungerade oder solche auf Zwischenfrequenzen. Insbesondere eignet sich dieser Prädiktor auch zur Verwendung bei den bisher nur unzureichend berücksichtigten geraden Oberschwingungen niedriger Ordnung.

Mittels des erfindungsgemäßen Prädiktors kann somit eine wesent ^¬ lich genauere Abschätzung der vom Windpark insgesamt abgegebenen Oberschwingung erreicht werden, und zwar gerade auch und besonders in den Fällen, in denen die Oberschwingung ein Mischfall aus phasenstarrem und phasenvariablem Teil sind. Damit ergibt sich ein erheblicher Vorsprung im Vergleich zu der im Stand der Technik bisher üblichen Abschätzung basierend auf einer reinen skalaren Addition. Mit dem erfindungsgemäßen Prädiktor ist es damit ermöglicht, mehr Windenergieanlagen an einen gegebenen An- schluss anzuordnen.

Zweckmäßigerweise sind das Iterationsmodul und das Summiermodul kombiniert zu einem Rekursionsmodul. Damit kann auf einfache und gut nachvollziehbare Weise die Oberschwingung für eine zunehmen ^¬ de Anzahl von Windenergieanlagen im Windpark bestimmt werden. Damit ist es möglich, die Rekursion dann abzubrechen, wenn der Oberschwingungsgrenzwert erreicht ist. Der Oberschwingungsprä- diktor gibt auf diese Weise unmittelbar die Anzahl der maximal anschließbaren Windenergieanlagen für einen Windpark direkt aus. Für die praktische Anwendung ist dies eine erhebliche Erleichte ^¬ rung .

Vorzugsweise weist das Berechnungsmodul eine Offset-Einheit für einen phasenstarren Teil der Oberschwingung der Windenergieanla- ge, einen Rotator für einen phasenvariablen Teil der Oberschwingung und einen Vektorisierer auf, der dazu ausgebildet ist, den Ersatzvektor basierend auf Mittelwerten von Offset-Einheit und Rotator nach Betrag und Winkel zu bestimmen. Mittels der Offset- Einheit kann insbesondere der phasenstarre Teil auf einfache Weise bestimmt werden. Charakteristisch für diesen ist, dass er für die betrachteten Windenergieanlagen des Windparks typischerweise immer dieselbe Phase aufweist. Er kann also mittels der Offset-Einheit vereinfacht berechnet werden. Für den phasenvari ^¬ ablen Teil ist vorzugsweise der Rotator vorgesehen, der dem variabel umlaufenden Teil der Oberschwingung repräsentiert. Der phasenvariable Teil ist bei den verschieden betrachteten Wind ^¬ energieanlagen definitionsgemäß nicht starr verkuppelt, sondern unabhängig von einander. Bei einer vektoriellen Berechnung ergeben sich also unterschiedliche Richtungen für die Vektoren

(selbst wenn sie alle denselben Betrag bzw. Länge aufweisen), so dass sich wegen der Winkelverschiedenheit hierin die gewünschte Reduktion abbildet. Der Vektorisierer erfasst diese Unterschiede bezüglich der Winkel und berücksichtigt sie bei der Aufaddition, um schließlich den Ersatzvektor zu bestimmen, der durch einen bestimmten Betrag (Länge des Vektors) und Winkel (Phasenlage) bestimmt ist.

Bei einer besonders vorteilhaften Ausführungsform, die ggf. unabhängigen Schutz verdient, ist der Umrichter der Windenergieanlage bezüglich seines phasenstarren und phasenvariablen Teils der Oberschwingung einstellbar. Es hat sich gezeigt, dass Oberschwingungen bei modernen hochfrequenten Pulswechselrichtern zwar nicht (vollständig) vermieden werden können, jedoch in Bezug auf ihre mittlere Phasenlage beeinflusst werden können. In dem die Gewichtung zwischen phasenstarrem und phasenvariablem Teil verschoben wird, kann somit Betrag und vor allem Winkellage des Ersatzvektors beeinflusst werden. Damit können für die ein ^¬ zelnen Windenergieanlagen solche Werte eingestellt werden, dass sie sich in ihrem Einfluss auf die Oberschwingungen möglichst weitgehend gegenseitig neutralisieren. Damit ist es ermöglicht, mehr Windenergieanlagen im Windpark vorzusehen ohne den Oberschwingungsgrenzwert zu überschreiten. Dazu wirkt zweckmäßiger ^¬ weise der Oberschwingungsprädiktor mit einem Limiter zusammen. Dieser verstellt das Verhältnis zwischen phasenstarren und pha ^¬ senvariablen Teil derart, dass eine Reduktion auf den gewünschten Grenzwert erreicht ist. Der Limiter kann alternativ oder zusätzlich auch so ausgebildet sein, dass die Anzahl der Windenergieanlagen begrenzt wird, welche elektrische Energie abgeben, um so in kritischen Fällen die Einhaltung der Oberschwingungsgrenzwerts zu garantieren.

Vorzugsweise ist der Oberschwingungsprädiktor besonders für Oberschwingungen gerader Ordnung ausgebildet, und zwar insbesondere bis solchen zur 50. Ordnung. Besonders bevorzugt ist eine Auslegung für tiefe Oberschwingungen (bis zur 20. Ordnung) und/oder hohe Oberschwingungen in Bereich der Schaltfrequenz des Umrichters .

Die Erfindung erstreckt sich weiter auf ein Verfahren zum Begrenzen der von einem Windpark generierten Oberschwingung, wobei der Windpark mindestens zwei Windenergieanlagen umfasst, die jeweils einen Generator mit einem Umrichter zur Erzeugung elektrischer Energie aufweisen, wobei erfindungsgemäß vorgesehen ist für eine erste der Windenergieanlagen ein Berechnen eines komplexen Mittelwerts über mindestens eine Periode der Oberschwin ^¬ gung und daraus bestimmen eines ersten Ersatzvektors mit Betrag und Winkel, und für mindestens eine zweite der Windenergieanla ^¬ gen iteratives Bestimmen eines zweiten Ersatzvektors auf diesel ^¬ be Weise, und schließlich Verwenden des ersten Ersatzvektors als Fußpunkt für den zweiten Ersatzvektor zur Aufsummation und abschließendes Bilden eines Gesamtvektors und dessen Vergleich mit dem Oberschwingungsgrenzwert. Unter Fußpunkt eines Vektors wird die der Spitze gegenüberliegende Seite verstanden.

Zur näheren Erläuterung des Verfahrens wird auf die vorstehend zu dem Windpark gegebene Erläuterung verwiesen.

Die Erfindung erstreckt sich ferner auf einen entsprechenden Oberschwingungsprädiktor für einen Windpark und das entsprechende Verfahren. Es wird auf vorstehende Erläuterungen verwiesen. Zusätzlich wird angemerkt, dass der Oberschwingungsprädiktor zweckmäßigerweise mit einem Simulationsinterface versehen sein kann. Dieses dient dazu, als Schnittstelle zu einer Simulations ^¬ umgebung zu fungieren, in der der Windpark mit seinen Windenergieanlagen implementiert ist. Dies ermöglicht eine Anwendung des erfindungsgemäßen Oberschwingungsprädiktors auch bereits bei Windparks in deren Entwurfs- und Planungsphase, also ohne dass sie bereits realisiert zu sein müssen. Damit ergibt sich eine erhebliche Erweiterung des Anwendungsspektrums des erfindungsge ^¬ mäßen Oberschwingungsprädiktors.

Die Erfindung wird nachfolgend in Bezugnahme auf die beigefügten Zeichnungen erläutert, in der ein vorteilhaftes Ausführungsbei ^¬ spiel dargestellt ist. Es zeigen:

Figur 1: eine Übersichtsdarstellung eines Windparks mit einem

Oberschwingungsprädiktor gemäß der Erfindung;

Figur 2: eine alternative Ausführungsform mit dem erfindungsge ^¬ mäßen Oberschwingungsprädiktor;

Figur 3: ein Diagramm für die mit dem Oberschwingungsprädiktor erreichte Reduktion der Oberschwingung; Figur 4: eine Darstellung des phasenstarren und phasenvariablen Anteils der Oberschwingung einer Windenergieanlage;

Figur 5: Beträge und Mittelwerte der Oberschwingung über eine

Periode hinweg dargestellt;

Figur 6: eine grafische Veranschaulichung des Prädiktionsverfahrens für drei Windenergieanlagen; und

Figur 7: eine erweiterte Darstellung für eine große Anzahl von

Windenergieanlagen im Windpark.

Ein Ausführungsbeispiel für einen Windpark mit einem erfindungs ^¬ gemäßen Oberschwingungsprädiktor ist in Figur 1 dargestellt. Der Windpark umfasst eine Mehrzahl von Windenergieanlagen 1, wobei im dargestellten Ausführungsbeispiel drei Windenergieanlagen vorgesehen sind. Sie sind gleichartig aufgebaut und zur Identi ^¬ fikation mit den römischen Zahlen I, II bzw. III gekennzeichnet. Die Windenergieanlagen 1 weisen jeweils einen doppelt gespeisten Asynchrongenerator 10 mit einem Umrichter 11 auf und sind zur Abgabe der erzeugten elektrischen Energie an ein Sammelnetz 2 angeschlossen, das als ein Mittelspannungsnetz ausgebildet ist. Es führt zu einem Anschlusspunkt 3, der üblicherweise einen Hochspannungstransformator (nicht dargestellt) zum Anschluss an ein Hochspannungsnetz 9 beinhaltet.

Weiter ist ein Windparkmaster 5 vorgesehen, welcher eine übergeordnete Regelung für die Windenergieanlagen des Windparks ausübt. Der Windparkmaster 5 ist über ein gesondertes Signallei ^¬ tungsnetz 4 mit den Windenergieanlagen 1 verbunden.

Außerdem angeschlossen an das Signalleitungsnetz 4 ist ein

Oberschwingungsprädiktor 6. Er umfasst ein Berechnungsmodul 61, ein Iterationsmodul 62 sowie ein Summiermodul 63. Das Iterati- onsmodul 62 und das Summiermodul 63 sind zu einem Rekursionsmo ^¬ dul 64 zusammengefasst . An dem Berechnungsmodul 61 ist ein Ein ^¬ gang 65 für die Ordnung der vom Oberschwingungsprädiktor zu behandelnden Oberschwingung n angeordnet. Üblicherweise werden hier kleine gradzahlige Ordnungszahlen vorgegeben. Weiter ist an dem Oberschwingungsprädiktor 6 ein Eingang 66 für einen Oberschwingungsgrenzwert vorgesehen. Ein Ausgangsanschluss 67 des Oberschwingungsprädiktors 6 ist rückgeführt zu einem Limiter 68, der an dem Parkmaster 5 angeordnet ist.

Über das Signalleitungsnetz 4 werden von der Windenergieanlage 1 Signale zu der zu behandelnden Oberschwingung bereitgestellt, und zwar insbesondere in Bezug auf den phasenstarren und phasenvariablen Anteil. In der Regel liegen diese Informationen in einer Betriebssteuerung der Windenergieanlage 1 vor und können daraus abgefragt werden, da diese über entsprechende Sensoren in der Windenergieanlage 1 selbst verfügt. Es kann aber auch vorge ^¬ sehen sein, dass die entsprechenden Werte durch Sensoren 14 an den einzelnen Windenergieanlagen gesondert erfasst und über das Signalleitungsnetz 4 übermittelt werden.

Aus den somit angelegten Werten bestimmt das Berechnungsmodul 61 einen komplexen Mittelwert der Oberschwingung über mindestens eine volle Periode und berechnet daraus einen Ersatzvektor. Dies wird nachfolgend unter Bezugnahme auf die Figuren 4 und 5 näher dargestellt .

In Figur 4 ist in der komplexen Ebene eine Oberschwingung einer der Windenergieanlagen, als Beispiel dient die Windenergieanlage I, dargestellt. Die Oberschwingung weist einen phasenstarren Teil sowie einen phasenvariablen Teil auf. Das Koordinatensystem für die komplexe Darstellung der Oberschwingung ist so gewählt, dass der phasenstarre Teil auf der reellen Achse liegt. Der pha ^¬ senstarre Teil ist in Figur 4 mit einem gestrichelten horizonta- len Pfeil dargestellt. Zu dem phasenstarren Teil tritt hinzu der phasenvariable Teil, der als ein kreisförmig umlaufender Vektor dargestellt ist mit einem dünnen durchgezogenen Pfeil. Dieser Pfeil läuft um (rotiert) mit der sich aus der Ordnung der Oberschwingung ergebenden Frequenz. Der Gesamtwert der Oberschwingung wird zu jedem beliebigen Zeitpunkt durch die Spitze des dünnen Pfeils, dessen Fußpunkt an der Spitze des gestrichelten Pfeils für den phasenstarren Teil liegt, repräsentiert. Daraus ergeben sich die in Figur 5 dargestellten Kennwerte. Der Realteil der Oberschwingung ist in Figur 5 durch die gestrichelte Linie dargestellt, und der Imaginärteil der Oberschwingung ist durch die durchgezogene Linie dargestellt. Der Betrag der kom ^¬ plexen Zahl, die durch den jeweiligen Realteil und Imaginärteil definiert ist, wird durch die strichpunktierte Linie darge ^¬ stellt. Auffällig ist, dass diese Linie anders als die Linien für den Real- bzw. Imaginärteil keinen sinusförmigen Verlauf aufweist, sondern vielmehr im Bereich um den Winkel 180° eine Abplattung aufweist. Dies ergibt sich daraus, dass bei Betrach ^¬ tung der Pfeilspitze der durchgezogenen Linie in Figur 4 diese bei dem Durchlauf am linken Ende des Kreises, also in dem der Imaginärachse zugewandten Bereich, einen beinahe gleichmäßigen Abstand zu dem Ursprung aufweist. Das bedeutet, dass in diesem Bereich der Betrag des vom Ursprung zu der Spitze des durchgezo ^¬ genen Pfeils laufende Vektor für die gesamte Oberschwingung minimal ist. Die Abplattung resultiert also letztlich aus der Überlagerung der linearen Verschiebung durch den phasenstarren Anteil, der wie ein Offset wirkt, und der Rotation durch den phasenvariablen Teil.

Weiter in Figur 5 dargestellt sind als Horizontallinien die Mittelwerte für den Realteil (gestrichelt) und der Mittelwert für den Gesamtwert der Oberschwingung (phasenstarrer und phasenvariabler Teil) in strichpunktierter Darstellung. Der Mittelwert für den komplexen Vektor der Oberschwingung ist ebenfalls komplex und ist in Figur 4 durch ein Kreuz dargestellt. Dieser dient erfindungsgemäß als Ersatzwert für die Oberschwingung, dessen Be ^¬ trag der mittleren Amplitude der Oberschwingung entspricht. Das Berechnungsmodul 61 liefert an seinem Ausgang diesen komplexen Mittelwert als Ersatzvektor.

Das Iterationsmodul 62 sorgt dafür, dass die vorstehend am Bei ^¬ spiel der Windenergieanlage I dargestellte Berechnung für die weiteren Windenergieanlagen, im Beispiel II und III, entspre ^¬ chend wiederholt wird. Das Summiermodul 63 addiert diese Werte auf. Dies ist in Figur 6 visualisiert . Sie enthält drei kreisar ^¬ tige Darstellungen, die mit den römischen Zahlen I, II und III bezeichnet sind für die drei Windenergieanlagen. Die Darstellung des ersten Kreises entspricht derjenigen in Figur 4 für die Windenergieanlage I. Deren als Ersatzwert fungierender Ersatz ^¬ vektor ist durch den strichpunktierten Pfeil nahe der römischen Zahl I dargestellt. Ausgehend von diesem Ersatzwert, also gra ^¬ fisch gesehen ausgehend von der Spitze des Ersatzvektors, wird die Berechnung für die zweite Windenergieanlage II durchgeführt. Das bedeutet, dass er als Nullpunkt für den phasenstarren Anteil fungiert, der zu einer horizontalen Verschiebung entlang der realen Achse führt. Dies ist dargestellt durch den gestrichelten Pfeil nahe der Ziffer II. Um die Spitze dieses Pfeils ist wie ^¬ derum der Drehzeiger für den phasenvariablen Teil in einer durchgezogenen Linie dargestellt. Da die zweite Windenergieanla ^¬ ge II eine andere mittlere Phasenlage hat bezüglich ihrer pha ^¬ senvariablen Oberschwingung ist der für sie auf die gleiche Weise wie für die Windenergieanlage I berechnete Ersatzvektor

(strichpunktierter Pfeil nahe der Ziffer II) anders orientiert als der Ersatzvektor für die Windenergieanlage I. Es erfolgt al ^¬ so eine vektorielle Addition der beiden Ersatzvektoren. Mittels des Iterationsmoduls 62 und des Summiermoduls 63 wird dieser Vorgang für die dritte Windenergieanlage III wiederholt, woraus sich der dritte Ersatzvektor ergibt, welcher durch den strich- punktierten Pfeil nahe der Ziffer III in Figur 6 dargestellt ist. Das Summiermodul 63 bildet schließlich aus diesen drei Er ^¬ satzvektoren einen Gesamtvektor T, der durch eine dicke durchgezogene Linie dargestellt ist. Er repräsentiert den Mittelwert des von den Windenergieanlagen I - III unter Berücksichtigung der jeweiligen phasenstarren und phasenvariablen Anteile und deren relativer Phasenlage insgesamt erzeugten Oberschwingungsbei ^¬ trags .

Führt man dieses Verfahren für eine größere Anzahl von Windenergieanlagen durch, beispielsweise für etwa 100 Windenergieanla ^¬ gen, so ergibt sich der in Figur 7 dargestellte Verlauf. Man er ^¬ kennt aus dem Gesamtverlauf eine Kennlinie, wie sie qualitativ einer Quadratwurzeldarstellung entspricht. Das bedeutet, dass gemäß den erfindungsgemäßen Oberschwingungsprädiktor die gesamten Oberschwingungen eines Windparks einer entsprechend großen Anzahl N von Windenergieanlagen sich nicht entsprechend N-fach erhöhen, sondern etwa einer 1/V(N) folgt. Damit kann mit dem erfindungsgemäßen Prädiktor auch für Windenergieanlagen, deren Oberschwingung gemischt sind in dem Sinne, dass sie phasenstarre und phasenvariable Anteile aufweist, das Gesamtverhalten eines Windparks entsprechend abgeschätzt werden. Das Einhalten der Grenzwerte am Anschlusspunkt 3 an das öffentliche Netz 9 ist da ^¬ mit erleichtert.

Der erfindungsgemäße Oberschwingungsprädiktor kann dazu dienen, dass Umrichter 11 für Generatoren 10 der Windenergieanlagen 1 so angesteuert werden, dass deren Oberschwingung in Bezug auf den phasenstarren und phasenvariablen Teil so eingestellt sind, dass der von dem Oberschwingungsprädiktor 6 bestimmte Gesamtvektor T den Oberschwingungsgrenzwert LIMIT nicht überschreitet. Außer der Einstellung des Umrichters 11 kann auch vorgesehen sein, dass in kritischen Fällen die Anzahl N der Windenergieanlagen reduziert wird, um so eine Reduktion der Oberschwingung zu er- reichen. Es sei angemerkt, dass das Prinzip des Oberschwin- gungsprädiktors an sich für jede Art von Windenergieanlagen 1 angewendet werden kann, es hat sich insbesondere bei Windener ^¬ gieanlagen mit doppelt gespeisten Asynchrongenerator 10 bewährt. Gerade bei diesen tritt besonders das Problem des Mischfalls auf, also einer Oberschwingung mit phasenstarrem und phasenvariablem Teil.

Erfindungsgemäß ist es nicht zwingend erforderlich, dass der Oberschwingungsprädiktor 6 an einem echten Windpark implementiert ist. Es ist auch möglich, den Oberschwingungsprädiktor 6 an einem Simulationsinterface 7 (s. Fig. 2) anzuschließen, welches mit einem Simulationsrechner 70 zusammenwirkt. Das Simula ^¬ tionsinterface 7 zusammen mit dem Simulationsrechner 70 bilden das elektrische Verhalten des Windparks mit seinen Windenergie ^¬ anlagen 1, dem internen Leitungsnetz 2 bis hin zum Anschlusspunkt 3 ab. Auf diese Weise kann mittels des Oberschwingungsprä- diktors 6 das Verhalten des Windparks in Bezug auf seine Ober ^¬ schwingung bereits ermittelt und bei der Auslegung des Windparks berücksichtigt werden, bevor er überhaupt tatsächlich gebaut wird. Dies ermöglicht eine wesentlich genauere Auslegung des Windparks. Das Simulationsinterface 7 weist dazu einen Anschluss 74 zur Simulation des Signalnetzes mit den daran angeschlossenen Einrichtungen und einen Anschluss 77 zum Anlegen eines Ausgangssignals 67' des Oberschwingungsprädiktors 6 ab.

Die mit dem erfindungsgemäßen Oberschwingungsprädiktor möglichen Verbesserungen sind in Figur 3 dargestellt. Der angestrebte Ide ^¬ alzustand ist eine Reduktion der Oberschwingung mit einem Reduktionsfaktor RF gemäß einer Beziehung 1/V(N), wie sie herkömmlicherweise lediglich bei Oberschwingungen höherer Ordnung erreicht werden kann; dieser ist in Figur 3 durch eine punktierte Linie dargestellt. Mit dem erfindungsgemäßen Oberschwingungsprä ^¬ diktor kann eine Reduktion erreicht werden, wie sie durch die durchgezogene Linie dargestellt ist. Man erkennt, dass zwar kei ^¬ ne ganz so starke Reduktion erreicht wird wie bei den Ober ^¬ schwingungen höherer Ordnung, aber dennoch eine beträchtliche Verbesserung erreicht wird gegenüber dem bisherigen Fall, bei dem keine Reduktion, sondern ein Konstantwert von 1 gemäß der gestrichelten Linie anzunehmen war.

Damit kann mit dem erfindungsgemäßen Oberspannungsprädiktor die Oberschwingungsbelastung genauer eingeschätzt werden, was letztendlich den Anschluss von mehr Windenergieanlagen in einem Windpark bei einem gegebenen Anschluss mit vorgegebenen Oberschwingungsgrenzwert ermöglicht. Besonders erfreulich ist, dass mit dem Oberschwingungsprädiktor bereits bei verhältnismäßig kleinen Windparks mit zehn oder zwanzig Anlagen eine beträchtliche Re ^¬ duktion der Oberschwingungsbelastung ermöglicht ist.