Script Sonnenstand und Einstrahlung
-
Ich wollte jetzt mal versuchen die auf Wetterunderground veröffentlichte Solarradiation zur Berechnung der möglichen PV-Leistung zu verwenden - habe nämlich bis jetzt noch keine PV-Anlage. Da ich eine WU-Station in der näheren Umgebung habe die die Solardaten anbietet müssten die Daten ja relativ zuverlässig sein.
Nun bin ich mir nicht ganz sicher, ob ich in deinem Script einfach anstelle der um den Airmass-Wert korrigierten 1.367kW/m³ den WU-Wert verwenden kann, also konkret diese Zeile:
var Sincident = (1.367*Math.pow(0.78,Math.pow(airmass,0.6)));durch folgende ersetzen:
var Sincident = (getState("weatherunderground.0.current.solarradiation").val/1000);Ich erhalte damit Werte die schon ungefähr plausibel sein könnten, mir allerdings tendenziell etwas niedrig erscheinen. Kannst du das vielleicht mal vergleichen?
-
Dann müsste ich ja rechnen :lol:
Aber ich hätte da vielleicht was für dich. Muss es bloss schnell zusammenkopieren.
Ich habe die Tage mit dem Solarwetter-Adapter getestet. Damit habe ich zusätzlich zur möglichen, die zu erwartenden min in max Werte mit geschrieben. Damit könntest vielleicht auch was anfangen.
Aber Vorschlag zur Güte; ich hau das mal als weitere Linie mit rein. und dann wissen wir morgen oder so mehr. – ist allerdings dunkles Wetter angesagt die nächsten Tage --
-
Ist das hier etwas für euch?
und:
Gruß
Rainer
-
-
Ich wollte jetzt mal versuchen die auf Wetterunderground veröffentlichte Solarradiation zur Berechnung der möglichen PV-Leistung zu verwenden - habe nämlich bis jetzt noch keine PV-Anlage. Da ich eine WU-Station in der näheren Umgebung habe die die Solardaten anbietet müssten die Daten ja relativ zuverlässig sein.
Nun bin ich mir nicht ganz sicher, ob ich in deinem Script einfach anstelle der um den Airmass-Wert korrigierten 1.367kW/m³ den WU-Wert verwenden kann, also konkret diese Zeile:
var Sincident = (1.367*Math.pow(0.78,Math.pow(airmass,0.6)));durch folgende ersetzen:
var Sincident = (getState("weatherunderground.0.current.solarradiation").val/1000);Ich erhalte damit Werte die schon ungefähr plausibel sein könnten, mir allerdings tendenziell etwas niedrig erscheinen. Kannst du das vielleicht mal vergleichen? `
Mit dem ersten Licht, möchte ich sagen, dass die tendenziell Faktor 10 zu klein sind. Wie sich das auf die Funktion im laufe des Tages auswirkt… später.Jedoch scheint es als ob die Radiation GHI gemessen wird. Das ist horizontal und müsste wohl auch noch eingerechnet werden.
-
Das mit dem Faktor 10 dachte ich anfangs auch - dann hatte ich aber tagsüber Werte von bis zu 40kW - schön wär's bei 9kWP ;)
Dass der der Radiation-Wert horizontal ist sollte doch durch die Zeile
var fraction = Math.cos(altitude*4*Math.asin(1)/360)*Math.sin(modtilt*4*Math.asin(1)/360)*Math.cos(azimuth*4*Math.asin(1)/360-modazi*4*Math.asin(1)/360)+Math.sin(altitude*4*Math.asin(1)/360)*Math.cos(modtilt*4*Math.asin(1)/360);berücksichtigt werden - oder?
-
Nö, schätze das ist die Abweichung zur Sichtlinie.
Bin auf Arbeit… Muss das am Abend genauer betrachte.
Memo an mich: https://pvpmc.sandia.gov/modeling-steps ... reflected/ lesen. Und das auch: http://www.solarconsultingservices.com/ ... ASES08.pdf
-
Das ist mir viel zu kompliziert. Wenn man den Wert mit 2.5 multipliziert, schauts aj soweit gut aus…Aber ich bin mir recht sicher, das die Werte 'Kack' -wir streichen das- ungenau sind.
Zudem nützt mir das nichts. Für die Abschätzung, was die Module dereinst mal bringen, gibts im INet gute Quellen mit guten Modellen.
-
Ich habe leider keine Photovoltaikanlage auf dem Dach, nur eine Solarthermieanlage :?
Kann man eigentlich das Script auch für diese benutzten oder ist das eher was für die Stromerzeugung gedacht?
Muss ich dann was am Script anpassen, damit ich ungefähre Werte erhalte?
Gruß
Adrian
-
Nachdem ich mein Script jetzt seit ca. 2 Wochen laufen habe würde ich sagen dass die Werte ohne Korrekturfaktor doch einigermaßen realistisch sein könnten. Natürlich schwanken die Werte die von WU kommen relativ stark, eben so wie bei der Messstation gerade mal ne Wolke vorbeikommt. Aber als langfristigen Richtwert denke ich kann man's gebrauchen.
Nun weiß ich nicht, wie der Wirkungsgrad bei Solarthermie ggü. Photovoltaik ist, aber prinzipiell sollte das schon einen Anhaltspunkt liefern wenn's dir nicht gerade auf 50W ankommt.
-
Die solare Einstahlung ist die Gleiche, unabhängig davon, wie Du sie nutzt. Der Wirkungsgrad ist bei thermischen Anwendungen wesentlich höher.
Das kannst Du am Parameter 'modeff' anpassen. Jedoch ist mir nicht exakt klar was Du damit visualisieren möchtest.
-
Nachdem ich mein Script jetzt seit ca. 2 Wochen laufen habe würde ich sagen dass die Werte ohne Korrekturfaktor doch einigermaßen realistisch sein könnten. Natürlich schwanken die Werte die von WU kommen relativ stark, eben so wie bei der Messstation gerade mal ne Wolke vorbeikommt. Aber als langfristigen Richtwert denke ich kann man's gebrauchen.
Nun weiß ich nicht, wie der Wirkungsgrad bei Solarthermie ggü. Photovoltaik ist, aber prinzipiell sollte das schon einen Anhaltspunkt liefern wenn's dir nicht gerade auf 50W ankommt.
-
Damit meiner Einer auch mal etwas an die Community zurückgeben kann:
Ich habe das Sonnenstand-Script aus dem Forum mit einem Python Programm von meinem vorherigen Solar-Pi gebastel verheiratet.
Herausgekommen ist eine Möglichkeit, die theoretische Leistung einer PV-Anlage abzuschätzen und zu loggen. Damit das etwas wissenschaftlich ausschaut, hab ich das mit allerlei Formelzeugs aus dem Internet vollgepackt.
Machen tue ich das so:
// calculates the sun position, path and power throughout the day // based on from http://www.stjarnhimlen.se/comp/tutorial.html // most var-names are identical to above tutorial // combined with Sonnenstand-Script: paul53, pix; 06.07.2015 nach ioBroker Forum http://forum.iobroker.net/viewtopic.php?f=21&t=975&sid=6f0ba055de5f82eed6809424f49ca93b#p7635 // Angepasst ykuendig 12.02.17; stringify im log, ack = true bei setState und Diverses var suncalc = require('suncalc'), result = getObject("system.adapter.javascript.0"), lat = result.native.latitude, long = result.native.longitude; var modtilt = 38; //Dachneigung in Grad (Solar panel's tilt angle) var modazi = 200; //Ausrichtung des Hauses in Grad zB SSW (Solar panel's azimut) var modsufrace = 56.0; //Paneloberfläche in m2 (Solar panel's surface in sq. meters) var modeff = 0.183; //Annäherung an Panel-Wirkungsgrad zB 18 Prozent (modules efficiency correction) // Hier kann an einem klaren Tag etwas geschraubt werden ;-) var altitude; // Calculated Elevation var azimuth; // Calculated Azimuth createState('javascript.0.Solar.Sonnenstand.Elevation', 0, {unit: '°'}); createState('javascript.0.Solar.Sonnenstand.Azimut', 0, {unit: '°'}); createState('javascript.0.Solar.Sonnenstand.PanelPossible', 0, {unit: 'W'}); // ganz am Ende die setStates anpassen nicht vergessen! // Do not change below, until You know what You are doing! // ******************************************************** Math.degrees = function(radians) {return radians * 180 / Math.PI;}; Math.radians = function(degrees) {return degrees * Math.PI / 180;}; function Sonnenstand_berechnen () { var now = new Date(); var sunpos = suncalc.getPosition(now, lat, long); log("Script Sonnenstand; latitude : " + result.native.latitude + " / longitude: " + result.native.longitude,'debug'); log("Script Sonnenstand; sunpos: " + JSON.stringify(sunpos),'debug'); altitude = Math.degrees(sunpos.altitude); azimuth = Math.degrees(sunpos.azimuth) + 180; // The intensity of the direct component of sunlight throughout each day can be determined as // a function of air mass. based on: http://pveducation.org/pvcdrom/properties-of-sunlight/air-mass#formula var airmass = 1/Math.cos((90-altitude)*4*Math.asin(1)/360); // Sincident is the intensity on a plane perpendicular to the sun's rays in units of kW/m2 and AM is the air mass. // The value of 1.353 kW/m2 is the solar constant and the number 0.7 arises from the fact that about 70% of the radiation incident on the atmosphere is transmitted // to the Earth. The extra power term of 0.678 is an empirical fit to the observed data and takes into account the non-uniformities in the atmospheric layers. // ykuendig: use different values because of pv instead thermal panels var Sincident = (1.367*Math.pow(0.78,Math.pow(airmass,0.6))); var fraction = Math.cos(altitude*4*Math.asin(1)/360)*Math.sin(modtilt*4*Math.asin(1)/360)*Math.cos(azimuth*4*Math.asin(1)/360-modazi*4*Math.asin(1)/360)+Math.sin(altitude*4*Math.asin(1)/360)*Math.cos(modtilt*4*Math.asin(1)/360); // W/m² light intensity on the module * module's surface var SmoduleInt = Sincident * fraction * modsufrace * 1000; if(SmoduleInt<0) { SmoduleInt = 0; } // Module Effective in relation of the efficiency of the used panel var SmoduleEff = SmoduleInt * modeff; if( altitude < 0 ) { SmoduleInt = 0; SmoduleEff = 0; altitude = 0; } log("Script Sonnenstand; Erfolgreich gelaufen, Werte akzeptiert", "info"); log("Script Sonnenstand; airmass: " + airmass,"debug"); log("Script Sonnenstand; azimuth: " + azimuth,"debug"); log("Script Sonnenstand; altitude: " + altitude,"debug"); log("Script Sonnenstand; SmoduleInt: " + SmoduleInt,"debug"); log("Script Sonnenstand; SmoduleEff: " + SmoduleEff,"debug"); // Change ID to the created States setState("javascript.0.Solar.Sonnenstand.Elevation",altitude.toFixed(1), true); setState("javascript.0.Solar.Sonnenstand.Azimut",azimuth.toFixed(), true); setState("javascript.0.Solar.Sonnenstand.PanelPossible", SmoduleEff.toFixed(), true); } // -> Zyklisch schedule("*/1 * * * *", Sonnenstand_berechnen); Sonnenstand_berechnen(); // bei ScriptstartAussehen kann das zum Beispiel so:
Dafür muss aber Sonnenstand und Leistung mit einem der History-Adapter regelmässig gespeichert werden.
- was ja eigentlich zu erwarten war.
Support us
653
Online32.7k
Users82.5k
Topics1.3m
Posts