Merge h2 csp

.github/workflows/ci.yml

@@ -18,7 +18,7 @@ jobs:
           python-version: ${{ matrix.python-version }}
       - name: Install dependencies
         run: |
-          sudo apt-get update && sudo apt-get install -y libglpk-dev glpk-utils
+          sudo apt-get update && sudo apt-get install -y libglpk-dev glpk-utils coinor-cbc
           python -m pip install --upgrade pip
           pip install shapely==1.7.1
           pip install -r requirements.txt

.gitignore

@@ -6,16 +6,40 @@ __pycache__/
 *$py.class
 
 log/
+examples/HOPP_examples/
 HOPP_examples/
 HOPP-demos/
-tests/analysis/results/
+# needed for h2 test inputs: tests/analysis/results/
+WB_case_study/
 examples/H2 Analysis/results/reopt_precomputes/
 
+tests/hybrid/REoptResultsNoExportAboveLoad.json
+
+# Results
+examples/analysis/results/*.csv
+resource_files/wind
+resource_files/solar
+
+resource_files/grid/tunisia_est_grid_prices.csv
+resource_files/solar/Hazoua
+resource_files/solar/Beni_Miha
+
+# clean up after merge
+examples/Powerflow Analysis
+examples/analysis/Engie
+examples/As_Built_Analysis_Blue Creek.csv
+examples/Blue_Creek_Analysis.py
+examples/dudgeon_wind_extraction_validation.py
+examples/powerflow_details.csv
+examples/reopt_load.py
+examples/data
+
+
 .idea/*
 .idea/workspace.xml
 .vscode/
 
-.vscode/
+.vscode/*
 
 # Distribution / packaging
 .Python
@@ -109,22 +133,11 @@ ENV/
 *.png
 *.xlsx
 
-# Results
-examples/analysis/results/*.csv
-resource_files/wind
-resource_files/solar
 *.pkl
-
 .DS_Store
+*.val
+*.db
+*.db-shm
+*.gz
+*.db-wal
 
-examples/HOPP_examples/
-
-examples/Powerflow Analysis
-examples/analysis/Engie
-examples/As_Built_Analysis_Blue Creek.csv
-examples/Blue_Creek_Analysis.py
-examples/dudgeon_wind_extraction_validation.py
-examples/powerflow_details.csv
-examples/reopt_load.py
-examples/data
-tests/hybrid/REoptResultsNoExportAboveLoad.json

2030_SETO_targets.json

@@ -0,0 +1,74 @@
+{
+  "cost_info": {
+    "wind_installed_cost_mw": 0,
+    "solar_installed_cost_mw": 585860,
+    "storage_installed_cost_mw": 200000,
+    "storage_installed_cost_mwh": 150000,
+    "wind_bos_cost_mw": 0,
+    "solar_bos_cost_mw": 0,
+    "storage_bos_cost_mw": 0,
+    "storage_bos_cost_mwh": 0
+  },
+  "csp_costs": {
+    "site_spec_cost": 10,
+    "cost_sf_fixed": 0,
+    "heliostat_spec_cost": 50,
+    "tower_fixed_cost": 1818300,
+    "tower_exp": 0.0113,
+    "rec_ref_cost": 62428300,
+    "rec_ref_area": 1571,
+    "rec_cost_exp": 0.7,
+    "tes_spec_cost": 10,
+    "plant_spec_cost": 700,
+    "bop_spec_cost": 0,
+    "fossil_spec_cost": 0,
+    "contingency_rate": 7,
+    "land_spec_cost": 10000,
+    "csp.pt.cost.plm.percent": 0,
+    "csp.pt.cost.plm.per_watt": 0,
+    "csp.pt.cost.plm.fixed": 0,
+    "csp.pt.cost.epc.per_acre": 0,
+    "csp.pt.cost.epc.percent": 13,
+    "csp.pt.cost.epc.per_watt": 0,
+    "csp.pt.cost.epc.fixed": 0,
+    "sales_tax_frac": 80,
+    "sales_tax_rate": 5,
+    "tower_rec_cost_per_kwt": 100
+    },
+   "SystemCosts": {
+    "om_capacity_escal": 3.0,
+    "om_fixed_escal": 0,
+    "om_production_escal": 0,
+    "pv": {
+      "om_capacity": [4.800000],
+      "om_fixed": [0.000000],
+      "om_production": [0.000000]
+    },
+    "tower": {
+      "om_capacity": [44.000000],
+      "om_fixed": [0.000000],
+      "om_production": [3.500000]
+    },
+    "trough": {
+      "om_capacity": [44.000000],
+      "om_fixed": [0.000000],
+      "om_production": [4.000000]
+    },
+    "battery": {
+      "en_batt": 1,
+      "om_capacity": [4.8],
+      "batt_replacement_option": 2,
+      "batt_bank_replacement": [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
+                                0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
+                                0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
+                                0, 0, 0, 0, 0
+      ],
+      "batt_replacement_schedule_percent": [0, 0, 0, 0, 0, 0, 0, 0, 0, 100,
+                                            0, 0, 0, 0, 0, 0, 0, 0, 0, 100,
+                                            0, 0, 0, 0, 0, 0, 0, 0, 0, 50,
+                                            0, 0, 0, 0, 0
+      ],
+      "om_batt_replacement_cost": [200]
+    }
+  }
+}

RELEASE.md

@@ -1,6 +1,29 @@
 # Release Notes
 
-## Version 0.1.0, ~
+## Version 0.1.0.dev3, Mar. 11, 2022
+* Include CBC package data for Windows
+
+## Version 0.1.0.dev2, Mar. 11, 2022
+* Add CBC open-source dispatch solver as default (faster than GLPK)
+* Add Xpress and Gurobi AMPL commercial dispatch solvers
+* Add clustering of simulation days with some tuning (use with caution)
+* Remove taxes from grid dispatch model
+* Add weighting of hybrid financials by production or cost ratio
+* Update to csp dispatch cost parameters for absolute grid prices
+* Add csv file output option to driver
+* Add PV charging only and grid charging options for dispatch
+* Add key parameter scaling when scaling CSP tower capacities
+* Add capacity credit payment calculations based on dispatch
+* Reformulate grid dispatch model
+* Add desired schedule option to dispatch model to follow schedule at least cost
+* Improve class documentation
+* Fix general bugs and clean-up code
+
+## Version 0.1.0.dev1, Oct. 11, 2021
+* Updated requirements
+* Added presolve to GLPK LP solver
+
+## Version 0.1.0.dev, Oct. 10, 2021
 * Battery dispatch and simulation with example in examples/simulate_hybrid_wbattery_dispatch.py 
 * Separate layout-related functions from optimization code into Layout classes
 * Refactor Optimizer classes to use HybridSimulation with examples in examples/optimization/hybrid_npv.py and examples/optimization/hybrid_sizing_problem.py
@@ -9,6 +32,9 @@
 * Rename "Solar" classes to "PV"
 * Add ElectricityPrices class with some example input files in resource_files/grid
 * Add storage costs to CostCalculator
+* Add concentrating solar power (CSP) tower and trough configurations through pySSC wrapper
+* Add dispatch optimization model for CSP models
+* Add design evaluation methods to iterate on design variables through space sampling, single objective derivative-free optimization, and multi-objective optimization
 
 ## Version 0.0.5, Apr 30, 2021
 * Update PySAM requirements

alt_dev/Price Stats.ipynb

@@ -0,0 +1,147 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": 1,
+   "id": "4dbcf4e7",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import ipywidgets as widgets\n",
+    "from ipywidgets import interact, Layout\n",
+    "import plotly.express as px\n",
+    "import pandas as pd\n",
+    "import numpy as np\n",
+    "import os\n",
+    "\n",
+    "# Plotting imports\n",
+    "import plotly.graph_objects as go\n",
+    "from plotly.subplots import make_subplots\n",
+    "\n",
+    "import plotly.io as pio\n",
+    "PAPER_LAYOUT=dict(font=dict(family=\"Computer Modern\", size=18),\n",
+    "                  margin=dict(t=40))\n",
+    "pio.templates[\"paper\"] = go.layout.Template(layout=PAPER_LAYOUT)\n",
+    "pio.templates.default = \"paper\"\n",
+    "\n",
+    "def chunks(lst, n):\n",
+    "    for i in range(0, len(lst), n):\n",
+    "        yield lst[i:i + n]\n",
+    "    "
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 2,
+   "id": "3998d12c",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "path = '../resource_files/grid'\n",
+    "\n",
+    "options = os.listdir(path)\n",
+    "data = [pd.read_csv(os.path.join(path, file), header=None)[0] for file in options]\n",
+    "\n",
+    "def dayPlot(x):\n",
+    "    idx = options.index(x)\n",
+    "    xvals = np.arange(24)\n",
+    "    y = data[idx]\n",
+    "    s = None\n",
+    "    \n",
+    "    fig = go.Figure()\n",
+    "    \n",
+    "    for chunk in chunks(y, 24):\n",
+    "        d = np.array(chunk, ndmin=2)\n",
+    "        \n",
+    "        if s is None:\n",
+    "            s = d            \n",
+    "        else:\n",
+    "            s = np.concatenate([s, d], axis=0)\n",
+    "            \n",
+    "        fig.add_trace(go.Scatter(x=xvals, y=chunk, line=dict(color='gray'), showlegend=False))\n",
+    "    \n",
+    "    fig.add_trace(go.Scatter(x=xvals, y=s.min(axis=0), line=dict(color='black', dash='dash'), showlegend=False))\n",
+    "    fig.add_trace(go.Scatter(x=xvals, y=s.max(axis=0), line=dict(color='black', dash='dash'), showlegend=False))\n",
+    "    fig.add_trace(go.Scatter(x=xvals, y=s.mean(axis=0), line=dict(color='black'), showlegend=False))\n",
+    "    \n",
+    "    fig.update_layout(title=x)\n",
+    "    fig.show()\n",
+    "    \n",
+    "    for i in range(2):\n",
+    "        print('hour', ''.join([f\"{i:7.0f}\" for i in xvals[i*12:12+i*12]]))\n",
+    "        print(' max', ''.join([f\"{i:7.0f}\" for i in s[:,i*12:12+i*12].max(axis=0)]))\n",
+    "        print('mean', ''.join([f\"{i:7.0f}\" for i in s[:,i*12:12+i*12].mean(axis=0)]))\n",
+    "#         print('mode', ''.join([f\"{i:7.0f}\" for i in s[:,i*12:12+i*12].mode(axis=0)]))\n",
+    "        print(' min', ''.join([f\"{i:7.0f}\" for i in s[:,i*12:12+i*12].min(axis=0)]))\n",
+    "        print()\n",
+    "        "
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 3,
+   "id": "34ade27e",
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "f3a2047c82ae4733aff8ad6a6cddad5a",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "interactive(children=(Dropdown(description='x', options=('pricing-data-2020-IronMtn-002.csv', 'test_prices.csv…"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "text/plain": [
+       "<function __main__.dayPlot(x)>"
+      ]
+     },
+     "execution_count": 3,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "interact(dayPlot, x=options) #('test_prices.csv')"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "ca3ac299",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "[np.isnan(d).sum() for d in data]"
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3",
+   "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.8.10"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}