PyPSA: Python for Power System Analysis

1 About

PyPSA stands for "Python for Power System Analysis". It is pronounced "pipes-ah".

PyPSA is a free software toolbox for simulating and optimising modern power systems that include features such as conventional generators with unit commitment, variable wind and solar generation, storage units, sector coupling and mixed alternating and direct current networks. PyPSA is designed to scale well with large networks and long time series.

As of 2017 PyPSA is under heavy development and therefore it is recommended to use caution when using it in a production environment. Some APIs may change - the changes in each PyPSA version are listed in the release notes.

PyPSA was initially developed by the Renewable Energy Group at FIAS to carry out simulations for the CoNDyNet project, financed by the German Federal Ministry for Education and Research (BMBF) as part of the Stromnetze Research Initiative.

2 Download

See the installation instructions.

You can also see/download the code directly from github.

4 What PyPSA does and does not do (yet)

PyPSA can calculate:

  • static power flow (using both the full non-linear network equations and the linearised network equations)
  • linear optimal power flow (optimisation of power plant and storage dispatch within network constraints, using the linear network equations, over several snapshots)
  • security-constrained linear optimal power flow
  • total electricity system investment optimisation (using linear network equations, over several snapshots simultaneously for optimisation of generation and storage dispatch and investment in the capacities of generation, storage and transmission)

It has models for:

  • meshed multiply-connected AC and DC networks, with controllable converters between AC and DC networks
  • standard types for lines and transformers following the implementation in pandapower
  • conventional dispatchable generators with unit commitment
  • generators with time-varying power availability, such as wind and solar generators
  • storage units with efficiency losses
  • simple hydroelectricity with inflow and spillage
  • coupling with other energy carriers
  • basic components out of which more complicated assets can be built, such as Combined Heat and Power (CHP) units, heat pumps, resistive Power-to-Heat (P2H), Power-to-Gas (P2G), battery electric vehicles (BEVs), etc.; each of these is demonstrated in the examples

Functionality that will definitely be added soon:

  • Multi-year investment optimisation
  • Simple RMS simulations with the swing equation
  • Distributed active power slack
  • Non-linear power flow solution using analytic continuation in the complex plane following GridCal

Functionality that may be added in the future:

  • Short-circuit current calculations
  • Dynamic RMS simulations
  • Small signal stability analysis
  • Interactive web-based GUI with SVG
  • OPF with the full non-linear network equations
  • Dynamic EMT simulations
  • Unbalanced load flow
  • Port to Julia

5 Example scripts as Jupyter/iPython notebooks

There are extensive examples available as Jupyter/iPython notebooks.

6 Screenshots

line-loading.png

Figure 1: Line loading with high wind feed-in in North Germany.

lmp.png

Figure 2: Nodal prices with high wind feed-in in North Germany.

reactive-power.png

Figure 3: Reactive power in Germany.

stacked-gen.png

Figure 4: Aggregated feed-in over a day.

storage-scigrid.png

Figure 5: Storage operation.

scigrid-curtailment.png

Figure 6: Curtailment of wind.

meshed-ac-dc.png

Figure 7: Meshed AC-DC hybrid nework.

euro-pie-pre-7-branch_limit-1-256.png

Figure 8: Optimised capacities of generation and storage for a 95% reduction in CO2 emissions in Europe compare to 1990 levels.

legend-flat.png

7 What PyPSA uses under the hood

PyPSA is written and tested to be compatible with both Python 2.7 and Python 3.5.

It leans heavily on the following Python packages:

  • pandas for storing data about components and time series
  • numpy and scipy for calculations, such as linear algebra and sparse matrix calculations
  • pyomo for preparing optimisation problems (currently only linear)
  • networkx for some network calculations
  • py.test for unit testing
  • logging for managing messages

The optimisation uses pyomo so that it is independent of the preferred solver (you can use e.g. the free software GLPK or the commercial software Gurobi).

The time-expensive calculations, such as solving sparse linear equations, are carried out using the scipy.sparse libraries.

8 Mailing list

PyPSA has a Google Group forum / mailing list.

9 Citing PyPSA

If you use PyPSA for your research, we would appreciate it if you would cite the following preprint paper (which has not yet been through peer review):

If you want to cite a specific PyPSA version, each release of PyPSA is stored on Zenodo with a release-specific DOI. This can be found linked from the overall PyPSA Zenodo DOI: zenodo.786605.svg.

10 Licence

PyPSA is released as free software under the GPLv3.