Skip to content

pythermogis

pythermogis is a Python package that provides API access to the ThermoGIS doublet simulations and economic calculations, ThermoGIS is developed and maintained by the Geological Survey of the Netherlands which is part of TNO.

Image 1
Image 1

pythermogis is currently used by the following research programs:

Image 1
  1. GoForward
  2. Geothermal Atlas for Africa (which is a part of the LEAP-RE project).

Usage

The intended usage of this pacakge is to import into your own projects so you can run your own doublet simulations.

Detailed examples can be found in the usage documentation page.

This package allows a user to simulate geothermal doublets providing the following parameters:

  • Top depth [m +ive downwards]
  • Thickness [m]
  • Porosity [0-1]
  • Permeability [0-1]
  • Temperature [°C] (Optional: if not provided, a temperature gradient will be used)

The code will simulate a Geothermal doublet, utilizing ThermoGIS with DoubletCalc1D as the engine to produce values of:

  • power [MW]
  • heat pump power [MW]
  • capex (Capital expenditure) [Million €]
  • opex (Operational expenditure) in first year [Million €]
  • utc (Unit Technical Cost [cts/kWH])
  • npv (Net-present-value) [Million €]
  • hprod (Discounted Heat Produced) [MWh]
  • cop (coefficient of performance of system) [-]
  • cophp (coefficient of performance of heat pump) [-]
  • pressure for driving the thermal loop (wells+reservoir) [bar]
  • flow rate [m³/hr]
  • reservoir depth [m]

The (geothermal) power and hprod, and utc only take into account geothermal heat, and do not include the heat pump contribution. Therefore the estimated power and hprod are lower than the actual power and hprod of the system, which includes heat from electricity consumption in the ESP and heatpump. The reason to do so is to avoid attributing produced energy which is not renewable, as renewable energy, such that utc for the goeothermal heat can be used as reference for feed-in tariffs and subsidies.

For details on how these parameters are calculated we refer users to the theory section as well as Thermogis calculation webpage.

For energy conversion to electricity, and chill the power concerns net electricty production and net chill respectively

The simulations are conducted in Java, and this package uses JPype to create a Python-Java binding.

It uses xarray to handle input and output parameters, a User should inform themselves and learn the usage of the Xarray package.

This module can be combined with the pygridsio package to read and process 2D raster data.

What can this package not do?

This package provides a narrow access point to the doublet simulations and UTC (Unit technical costs) economic calculations from ThermoGIS, it does not implement the following ThermoGIS processes/methodologies:

  • Stacking of Aquifers
  • Calculation of Potential maps
  • Calculation of Resources
  • Calculation of Overview maps
  • HTO-ATES simulations
  • Property Modelling
  • Temperature Modelling
  • Reading, and writing of grid files (Checkout pygridsio)