Globalstrahlung mittels PV-Anlagenleistung berechnen

Einklappen
X
 
  • Zeit
  • Anzeigen
Alles löschen
neue Beiträge
  • He.Slar
    Smart Home'r
    • 19.03.2018
    • 64

    #1

    Globalstrahlung mittels PV-Anlagenleistung berechnen

    Hallo zusammen, ich habe (mithilfe einer langen KI-Sitzung) einen Weg gefunden, wie man in Loxone über die PV-Leistung und einer ellenlangen Formel die Globalstrahlung relativ genau bestimmen kann.
    Die konstanten beziehen sich auf meine Anlage: Südausrichtung (180°), 26° Dachneigung, 8kWp
    Das möchte ich Euch nicht vorenthalten:

    Berechnung der Strahlungsstärke der Sonne mithilfe der PV-Leistung im Loxone-Formelbaustein

    1. Physikalische Grundgleichung
    Die Bestrahlungsstärke $G$ ($W/m^2$) ist die aktuelle Leistung ($P_{aktuell}$) geteilt durch die effektive Modulfläche und den Wirkungsgrad. Da wir mit der Nennleistung ($P_{peak}$) bei Standardbedingungen ($1000 W/m^2$) rechnen, gilt: G=PpeakGeometrieTemperaturPaktuell​​1000
    2. Der Geometrie-Faktor ($K_{geo}$)
    Dieser Faktor beschreibt, wie viel Licht aufgrund des Sonnenstands auf Ihr spezifisches Dach (26° Neigung, Süd) trifft.
    • Formel: $\cos(\theta) = \sin(h)\cos(26°) + \cos(h)\sin(26°)\cos(\alpha_s - 180°)$
    • Eingesetzte Konstanten: $\cos(26°) \approx \mathbf{0.8988}$ und $\sin(26°) \approx \mathbf{0.4384}$.
    • Ergebnis: 0.8988*SIN(RAD(I3)) + 0.4384*COS(RAD(I3))*COS(RAD(I2-180))
    3. Der Temperatur-Faktor ($f_{Temp}$)
    Module verlieren bei Hitze an Leistung ($-0,4%$ pro °C über 25°C). Die Modultemperatur setzt sich aus Lufttemperatur ($I4$) und der Aufheizung durch die Sonne zusammen.
    • Ansatz: $T_{Modul} = T_{Luft} + 0,03 \cdot G$
    • Faktor: $1 - 0,004 \cdot (T_{Modul} - 25)$
    • Da $G$ das Ergebnis ist, wird dieser Term in den Nenner verschoben und nach $G$ aufgelöst. Das erzeugt den Korrekturterm - 0.12 * I1 im Nenner (hergeleitet aus $0,004 \cdot 30 \cdot P_{kW}$).
    4. Umstellung und Einheiten (kW zu W/m²)
    • Einheiten: Da $I1$ in kW und $P_{peak}$ als 8 kW eingeht, muss das Ergebnis mit 1000 multipliziert werden, um von einem Verhältnis wieder auf W/m² zu kommen.
    • Auflösung: Der Term mit $G$ im Nenner wird auf die linke Seite gebracht, ausgeklammert und durch den Rest geteilt.
      • Aus $8000$ Watt werden 8 (da $I1$ in kW).
      • Der Zähler bleibt I1 * 1000.


    5. Finale Loxone-Syntax
    Zusammenfügen aller Teile, Entfernung von Leerzeichen und Nutzung der RAD()-Funktion für die Winkelumrechnung ergibt:
    (I1*1000)/((8*(0.8988*SIN(RAD(I3))+0.4384*COS(RAD(I3))*COS(R AD(I2-180))))*(1.1-0.004*I4)-(0.12*I1))

    Um Extremwerte (Division durch Werte nahe Null) bei Kleiner Leistung auszuschließen, kann noch eine Bedingung vorangestellt werden:
    IF(I1>0.05;IF(I3>5;(I1*1000)/((8*(0.8988*SIN(RAD(I3))+0.4384*COS(RAD(I3))*COS(R AD(I2-180))))*(1.1-0.004*I4)-(0.12*I1));0);0
    Zuordnung der Eingänge (Inputs)
    Verbinden Sie diese Variablen mit den entsprechenden Sensoren oder Berechnungen in Ihrer Loxone-Config:
    • I1: Aktuelle PV-Leistung
      • Einheit: kW (Kilowatt)
      • Quelle: Wechselrichter
    • I2: Sonnenrichtung (Azimut)
      • Einheit: Grad (°)
      • Bereich: 0° bis 360° (z.B. 180° für Süd).
      • Quelle: Zeitfunktionen
    • I3: Sonnenhöhe (Elevation)
      • Einheit: Grad (°)
      • Bereich: 0° bis 90°.
      • Quelle: Zeitfunktionen
    • I4: Außentemperatur
      • Einheit: °C (Grad Celsius)
      • Quelle: Außenfühler der Heizung über Heizkessel
    Ergebnis (Output)
    • Wert: Bestrahlungsstärke (Globalstrahlung)
    • Einheit: W/m² (Watt pro Quadratmeter)
    • Beispielwert: ~1000 bei voller Sonne, ~400 bei Wolken.

    Zusammenfassung der Variablen:
    • Zähler: Aktuelle Leistung skaliert auf W/m².
    • Nenner Teil 1: Theoretische Leistung bei aktueller Geometrie und Lufttemperatur.
    • Nenner Teil 2: Abzug des Leistungsverlusts, der durch die eigene Erwärmung unter Last entsteht.
    1xMS Gen1, 1xExtension, 1-Wire Ext, ioBroker (Raspberry 4)
    Sensorik und Aktorik größtenteils MDT (KNX TP),
    Mähroboter Kress (über ioBroker angebunden)
    Heidelberg Wallbox Energy Control (Modbus RTU),
    Phasenumschaltung als Ablaufsteuerung in Loxone realisiert,
    PV (Fronius Symo GEN24) mit Speicher BYD HVS 7.7,
    Tesla Model Y und Renault ZOE als Stromverwerter,
    Erster Versuchsaufbau mit MS und KNX in 2015, Hausbau 2016/17, Fertigstellung der Haustechnik im Spätherbst 2075
  • Prof.Mobilux
    Supermoderator
    • 25.08.2015
    • 5418

    #2
    LoxBerry: https://wiki.loxberry.de/start

    Kommentar

    • He.Slar
      Smart Home'r
      • 19.03.2018
      • 64

      #3
      Ok, danke für den Hinweis, waren wohl andere auch schon dran an dem Thema.
      Großer Vorteil meiner Lösung: man braucht genau einen Formelbaustein, das war´s.
      Natürlich funktioniert das so nur, wenn die komplette Anlage die gleiche Ausrichtung hat, oder man kann von einem String die Leistung separat aus dem Wechselrichter abfragen.
      Und die Formel kann man einfach reinkopieren. Man muss nur noch Peak-Leistung, Dachneigung und Ausrichtung anpassen.
      Klicke auf die Grafik für eine vergrößerte Ansicht

Name: image.png
Ansichten: 28
Größe: 9,0 KB
ID: 486650
      IF(I1>0.05;IF(I3>5;(I1*1000)/((8*(0.8988*SIN(RAD(I3))+0.4384*COS(RAD(I3))*COS(RAD(I2-180))))*(1.1-0.004*I4)-(0.12*I1));0);0​
      die 8 ist die Peak-Leistung
      0,8988 ist der Cosinus der Dachneigung, hier cos(26°),
      0,4384 ist der Sinus der Dachneigung, hier sin(26°),
      Die 180 weiter hinten steht für die Dachausrichtung, in meinem Fall sehr exakt Süd (180°), kennt man von der Beschattungssteuerung.

      Achtung: Wenn man im Formelbaustein eine Winkelfunktion verwendet, muss man die Grad erst ins Bogenmaß umrechnen RAD(Winkel in °), wie es an mehreren Stellen in der Formel auch gemacht wird. Nachdem sich diese Werte ja normalerweise nicht mehr ändern (Häuser verändern nur sehr selten ihre Dachneigung und Ausrichtung), kann man die aber auch einfach mit dem Taschenrechner ausrechnen und als Ergebnis eintragen.
      Will man die Berechnung in der Formel machen, müsste z.B. für die Dachneigung statt 0,8988 der Term "COS(RAD(Dachneigung in °))" stehen

      Wenn man deutlich woanders wie in Bayern lebt muss man vielleicht noch den Wert der theoretischen Globalstrahlung (1000) gleich am Anfang noch anpassen.
      Im Beispielbild sieht man ganz gut, dass I3 (Sonnenhöhe) bereits negativ ist. Das wird durch die zweite IF-Funktion ausgeblendet und der Ausgabewert bleibt 0.
      1xMS Gen1, 1xExtension, 1-Wire Ext, ioBroker (Raspberry 4)
      Sensorik und Aktorik größtenteils MDT (KNX TP),
      Mähroboter Kress (über ioBroker angebunden)
      Heidelberg Wallbox Energy Control (Modbus RTU),
      Phasenumschaltung als Ablaufsteuerung in Loxone realisiert,
      PV (Fronius Symo GEN24) mit Speicher BYD HVS 7.7,
      Tesla Model Y und Renault ZOE als Stromverwerter,
      Erster Versuchsaufbau mit MS und KNX in 2015, Hausbau 2016/17, Fertigstellung der Haustechnik im Spätherbst 2075

      Kommentar


      • challo
        challo kommentierte
        Kommentar bearbeiten
        Ist 8 die Peak Leistung der Module oder des Wechselrichters?
    Lädt...