Create mismatch_modifiers.ipynb

Author: echedey-ls

docs/tutorials/mismatch_modifiers.ipynb

@@ -0,0 +1,209 @@
+{
+ "cells": [
+  {
+   "attachments": {},
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "# Mismatch Modifiers\n",
+    "Learn to use mismatch modifiers with this notebook!\n",
+    "Feel free to add other models, be sure to update the index and give you credit ;)\n",
+    "\n",
+    "Table of contents:\n",
+    "1. [Setup](#setup)\n",
+    "1. [N. Martin & J. M. Ruiz Experimental Mismatch Modifier](#n-martin--j-m-ruiz-experimental-mismatch-modifier)\n",
+    "\n",
+    "Authors:\n",
+    "* Echedey Luis (@echedey-ls), 2023 Feb"
+   ]
+  },
+  {
+   "attachments": {},
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "## Setup\n",
+    "Let's prepare the environment:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 1,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Show matplotlib's figures in the notebook\n",
+    "%matplotlib inline\n",
+    "import matplotlib.pyplot as plt\n",
+    "\n",
+    "import pandas as pd\n",
+    "\n",
+    "# And pvlib\n",
+    "import pvlib"
+   ]
+  },
+  {
+   "attachments": {},
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "### N. Martin & J. M. Ruiz Experimental Mismatch Modifier\n",
+    "This modifier takes into account the responsivities to different spectrums,\n",
+    "characterized by the airmass and the clearness index, as two independent\n",
+    "variables. In fact, it is 3 different modifiers, each one for each component\n",
+    "(``poa_direct``, ``poa_sky_diffuse``, ``poa_ground_diffuse``)\n",
+    "\n",
+    "The formula for each component has three coefficients; we are lucky the authors\n",
+    "of this model computed fitting values for m-Si, p-Si and a-Si!\n",
+    "However, if you would like to compute and/or use your own values, keep reading."
+   ]
+  },
+  {
+   "attachments": {},
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "First step is get to the effective irradiance. For simplicity, we will copy the procedure explained in the tutorial ``tmy_to_power.ipynb`` [TODO: HOW CAN I LINK THIS?]. Please refer to it to get a more in depth explanation."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 2,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "site = pvlib.location.Location(40.4534, -3.7270, altitude=664,\n",
+    "                               name='IES-UPM, Madrid', tz='CET')\n",
+    "\n",
+    "surface_tilt = 40\n",
+    "surface_azimuth = 180  # Pointing South\n",
+    "\n",
+    "tmy_data, _, _, _ = pvlib.iotools.get_pvgis_tmy(site.latitude, site.longitude, map_variables=True,\n",
+    "                                                startyear=2005, endyear=2015)\n",
+    "tmy_data.index = [ts.replace(year=2022) for ts in tmy_data.index]\n",
+    "\n",
+    "solar_pos = site.get_solarposition(tmy_data.index)\n",
+    "\n",
+    "extra_rad = pvlib.irradiance.get_extra_radiation(tmy_data.index)\n",
+    "\n",
+    "poa_sky_diffuse = pvlib.irradiance.haydavies(surface_tilt, surface_azimuth, tmy_data['dhi'],\n",
+    "                            tmy_data['dni'], extra_rad,\n",
+    "                            solar_pos['apparent_zenith'], solar_pos['azimuth'])\n",
+    "\n",
+    "poa_ground_diffuse = pvlib.irradiance.get_ground_diffuse(surface_tilt, tmy_data['ghi'])\n",
+    "\n",
+    "aoi = pvlib.irradiance.aoi(surface_tilt, surface_azimuth, solar_pos['apparent_zenith'], solar_pos['azimuth'])\n",
+    "\n",
+    "# Let's consider this the irradiances without modifiers\n",
+    "poa_irrad = pvlib.irradiance.poa_components(aoi, tmy_data['dni'], poa_sky_diffuse, poa_ground_diffuse)\n",
+    "\n",
+    "# Following part will be needed later\n",
+    "thermal_params = pvlib.temperature.TEMPERATURE_MODEL_PARAMETERS['sapm']['open_rack_glass_polymer']\n",
+    "pvtemps = pvlib.temperature.sapm_cell(poa_irrad['poa_global'], tmy_data['temp_air'], tmy_data['wind_speed'], **thermal_params)\n",
+    "\n",
+    "sandia_modules = pvlib.pvsystem.retrieve_sam(name='SandiaMod')\n",
+    "sandia_module = sandia_modules['Canadian_Solar_CS5P_220M___2009_']"
+   ]
+  },
+  {
+   "attachments": {},
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Here comes the modifier. Let's calculate it and examine the introduced\n",
+    "difference.\n"
+   ]
+  },
+  {
+   "attachments": {},
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "That was a lot, yeah. But don't worry, now we can find the effective irradiance, the mismatch modifier (with the airmass and clearness index)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 9,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "Module type is: c-Si\n"
+     ]
+    }
+   ],
+   "source": [
+    "# First, let's find the airmass and the clearness index\n",
+    "# Little caution: default values for this model were fitted obtaining the airmass through the kasten1966 method, not used by default\n",
+    "airmass = site.get_airmass(solar_position=solar_pos, model='kasten1966')\n",
+    "clearness = pvlib.irradiance.clearness_index(ghi=tmy_data['ghi'],\n",
+    "                                             solar_zenith=solar_pos['zenith'],\n",
+    "                                             extra_radiation=extra_rad)\n",
+    "# Check module is m-Si (monocrystalline silicon)\n",
+    "print('Module type is: ' + sandia_module['Material'])  #-Reports 'c-Si'\n",
+    "\n",
+    "# Get the mismatch modifiers\n",
+    "modifiers = pvlib.spectrum.martin_ruiz_spectral_modifier(clearness,\n",
+    "                                                         airmass['airmass_absolute'],\n",
+    "                                                         cell_type='monosi')\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 16,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "poa_irrad_modified = pd.Series(dtype=pd.Float64Dtype)\n",
+    "poa_irrad_modified['poa_direct'] = poa_irrad['poa_direct'] * modifiers['direct']\n",
+    "poa_irrad_modified['poa_sky_diffuse'] = poa_irrad['poa_sky_diffuse'] * modifiers['sky_diffuse']\n",
+    "poa_irrad_modified['poa_ground_diffuse'] = poa_irrad['poa_ground_diffuse'] * modifiers['ground_diffuse']"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 15,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "22.054707712730973\n"
+     ]
+    }
+   ],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "venv",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.10.4"
+  },
+  "orig_nbformat": 4,
+  "vscode": {
+   "interpreter": {
+    "hash": "e7b76f25baca03aa641c501db0912de76daa352e2d97ceb6fcf3025f206f2928"
+   }
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 2
+}