Open Source Solar Forecasting for a Site
Project description
Quartz Solar Forecast
The aim of the project is to build an open source PV forecast that is free and easy to use.
The forecast provides the expected generation in kw
for 0 to 48 hours for a single PV site.
Open Climate Fix also provide a commercial PV forecast, please get in touch at quartz.support@openclimatefix.org
The current model uses GFS or ICON NWPs to predict the solar generation at a site
from quartz_solar_forecast.forecast import run_forecast
from quartz_solar_forecast.pydantic_models import PVSite
# make a pv site object
site = PVSite(latitude=51.75, longitude=-1.25, capacity_kwp=1.25)
# run model, uses ICON NWP data by default
predictions_df = run_forecast(site=site, ts='2023-11-01')
Which gives the following prediction
Installation
The source code is currently hosted on GitHub at: https://github.com/openclimatefix/Open-Source-Quartz-Solar-Forecas
Binary installers for the latest released version are available at the Python Package Index (PyPI)
pip install quartz-solar-forecast
You might need to install the following packages first
conda install -c conda-forge pyresample
Model
The model is a gradient boosted tree model and uses 9 NWP variables. It is trained on 25,000 PV sites with over 5 years of PV history, which is available here. The training of this model is handled in pv-site-prediction TODO - we need to benchmark this forecast.
The 9 NWP variables, from Open-Meteo documentation, are mentioned above with their appropariate units.
- Visibility (km), or vis: Distance at which objects can be clearly seen. Can affect the amount of sunlight reaching solar panels.
- Wind Speed at 10 meters (km/h), or si10 : Wind speed measured at a height of 10 meters above ground level. Important for understanding weather conditions and potential impacts on solar panels.
- Temperature at 2 meters (°C), or t : Air temperature measure at 2 meters above the ground. Can affect the efficiency of PV systems.
- Precipiration (mm), or prate : Precipitation (rain, snow, sleet, etc.). Helps to predict cloud cover and potentiel reductions in solar irradiance.
- Shortwave Radiation (W/m²), or dswrf: Solar radiation in the shortwave spectrum reaching the Earth's surface. Measure of the potential solar energy available for PV systems.
- Direct Radiation (W/m²) or dlwrf: Longwave (infrared) radiation emitted by the Earth back into the atmosphere. confirm it is correct
- Cloud Cover low (%), or lcc: Percentage of the sky covered by clouds at low altitudes. Impacts the amount of solar radiation reachign the ground, and similarly the PV system.
- Cloud Cover mid (%), or mcc : Percentage of the sky covered by clouds at mid altitudes.
- Cloud Cover high (%), or lcc : Percentage of the sky covered by clouds at high altitude
We also use the following features
- poa_global: The plane of array irradiance, which is the amount of solar radiation that strikes a solar panel.
- poa_global_now_is_zero: A boolean variable that is true if the poa_global is zero at the current time. This is used to help the model learn that the PV generation is zero at night.
- capacity (kw): The capacity of the PV system in kw.
- The model also has a feature to check if these variables are nans or not.
The model also use the following variables, which are currently all set to nan
- recent_power: The mean power over the last 30 minutes
- h_mean: The mean of the recent pv data over the last 7 days
- h_median: The median of the recent pv data over the last 7 days
- h_max: The max of the recent pv data over the last 7 days
Known restrictions
- The model is trained on UK MetOffice NWPs, but when running inference we use GFS data from Open-meteo. The differences between GFS and UK MetOffice, could led to some odd behaviours.
- It looks like the GFS data on Open-Meteo is only available for free for the last 3 months.
Evaluation
To evaluate the model we use the UK PV dataset and the ICON NWP dataset. All the data is publicly available and the evaluation script can be run with the following command
python scripts/run_evaluation.py
The test dataset we used is defined in quartz_solar_forecast/dataset/testset.csv
.
This contains 50 PV sites, which 50 unique timestamps. The data is from 2021.
The results of the evaluation are as follows The MAE is 0.1906 kw across all horizons.
Horizons | MAE [kw] | MAE [%] |
---|---|---|
0 | 0.202 +- 0.03 | 6.2 |
1 | 0.211 +- 0.03 | 6.4 |
2 | 0.216 +- 0.03 | 6.5 |
3 - 4 | 0.211 +- 0.02 | 6.3 |
5 - 8 | 0.191 +- 0.01 | 6 |
9 - 16 | 0.161 +- 0.01 | 5 |
17 - 24 | 0.173 +- 0.01 | 5.3 |
24 - 48 | 0.201 +- 0.01 | 6.1 |
If we exclude nighttime, then the average MAE [%] from 0 to 36 forecast hours is 13.0%.
Notes:
- The MAE in % is the MAE divided by the capacity of the PV site. We acknowledge there are a number of different ways to do this.
- It is slightly surprising that the 0-hour forecast horizon and the 24-48 hour horizon have a similar MAE. This may be because the model is trained expecting live PV data, but currently in this project we provide no live PV data.
Abbreviations
- NWP: Numerical Weather Predictions
- GFS: Global Forecast System
- PV: Photovoltaic
- MAE: Mean Absolute Error
- ICON: ICOsahedral Nonhydrostatic
- KW: Kilowatt
Contribution
We welcome other models
Contributors ✨
Thanks goes to these wonderful people (emoji key):
Peter Dudfield 💻 |
Megawattz 🤔 |
EdFage 📖 💻 |
Chloe Pilon Vaillancourt 📖 |
This project follows the all-contributors specification. Contributions of any kind welcome!
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for quartz_solar_forecast-0.0.9.tar.gz
Algorithm | Hash digest | |
---|---|---|
SHA256 | d483cd4709af84b6521480952b8d96a8f83171425e2d546d4dc71eaa8da6c892 |
|
MD5 | e5125ba7ec9637927c28a0eb8968a92d |
|
BLAKE2b-256 | 977695a3aa51301b6cf83b9d6b4bb723054e3923d4788cb1e32672cce955d62f |
Hashes for quartz_solar_forecast-0.0.9-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 12158596adaae81a865d5ea4ed3b9168a1660ca3dea7c7e0421e22051b34297f |
|
MD5 | 43fe9e113900deb00fb6d34583ed8f20 |
|
BLAKE2b-256 | a912693669147e7a7be9fecca4b32f4dc4692887c51e071d6484003d82a28f37 |