292 lines
11 KiB
Ruby
292 lines
11 KiB
Ruby
require 'open-uri'
|
|
EASY_ENERGY_TARIFFS = {}
|
|
|
|
# See https://www.belastingdienst.nl/wps/wcm/connect/bldcontentnl/belastingdienst/zakelijk/overige_belastingen/belastingen_op_milieugrondslag/tarieven_milieubelastingen/tabellen_tarieven_milieubelastingen
|
|
# Without VAT
|
|
ENERGY_TAX_KWH = { 2020 => 0.09770, 2021 => 0.09428, 2022 => 0.03679}
|
|
ODE_KWH = { 2020 => 0.0273, 2021 => 0.0300, 2022 => 0.0305}
|
|
# merge by adding values
|
|
TAX_KWH = ENERGY_TAX_KWH.merge(ODE_KWH){|key, energy_tax, ode| energy_tax + ode}
|
|
|
|
|
|
class Cost
|
|
|
|
attr_accessor :max_charge_kwh
|
|
attr_reader :battery, :entsoe, :zone
|
|
|
|
def initialize(zone="Amsterdam", battery_capacity=10.0, max_charge=5.0, storage_cost=0.05)
|
|
@zone = zone
|
|
@entsoe = Entsoe.new(zone, storage_cost)
|
|
@max_charge_kwh = max_charge
|
|
@battery = Battery.new(battery_capacity)
|
|
end
|
|
|
|
# Government reduced VAT to 9% from 1 July 2022 until 31 Dec 2022
|
|
def vat_at(date)
|
|
jul22 = Date.parse("2022-09-01")
|
|
dec22 = Date.parse("2022-12-31")
|
|
if (date >= jul22 && date <= dec22)
|
|
0.09
|
|
else
|
|
0.21
|
|
end
|
|
end
|
|
|
|
def format_cost(cost)
|
|
cost ? "EUR %0.03f" % cost : "EUR ?"
|
|
end
|
|
|
|
# Assume: usage_kwh_cost, return_kwh_cost already includes vat
|
|
def add_tax(formatted_hour,usage_kwh,usage_kwh_cost,return_kwh, return_kwh_cost)
|
|
return nil if (usage_kwh.nil? || usage_kwh_cost.nil? || return_kwh.nil? || return_kwh_cost.nil?)
|
|
|
|
vat = 1 + vat_at(Date.parse(formatted_hour))
|
|
year = Date.parse(formatted_hour).year
|
|
# calculate tax per kwh - which the sum of energy_tax and ode
|
|
tax = TAX_KWH[year]*vat
|
|
|
|
usage_kwh * (usage_kwh_cost + tax) - return_kwh * (return_kwh_cost + tax)
|
|
|
|
end
|
|
|
|
######################################################
|
|
# Easy Energy API (proprietary) - better to use entsoe
|
|
######################################################
|
|
|
|
# returns a hash with keys formatted "yyyy-mm-dd-hr" and values [usage, return]
|
|
# e.g. { 2021-12-16-06"=>[0.36058, 0.298] }
|
|
def easy_energy_tariffs(date)
|
|
p "Fetching EasyEnergy tariffs for %s" % date.strftime("%F-%H")
|
|
|
|
# calculate offset (date is in UTC, and we want to have tariffs in Amsterdam zone)
|
|
zone = 'Amsterdam'
|
|
offset = DateTime.now.in_time_zone(zone).utc_offset
|
|
date = date.beginning_of_day.advance(seconds: offset)
|
|
|
|
url = "https://mijn.easyenergy.com/nl/api/tariff/getapxtariffs?startTimestamp=%s&endTimestamp=%s&grouping=" % [date.strftime("%F %T"),date.advance(:hours => 24).strftime("%F %T")]
|
|
#p urlvandaag
|
|
json = JSON.load(URI::open(url))
|
|
# advancing with 1 hrs (to offset for something?)
|
|
json.map{|t| [DateTime.parse(t["Timestamp"]).strftime("%F-%H"), [t["TariffUsage"], t["TariffReturn"]]]}.to_h
|
|
end
|
|
|
|
def easy_energy_rate(formatted_hour)
|
|
unless EASY_ENERGY_TARIFFS.key?(formatted_hour)
|
|
EASY_ENERGY_TARIFFS.merge!(easy_energy_tariffs(Date.parse(formatted_hour)))
|
|
end
|
|
EASY_ENERGY_TARIFFS[formatted_hour]
|
|
end
|
|
|
|
def easy_energy_cost(formatted_hour, usage_kwh, return_kwh)
|
|
return nil if (usage_kwh.nil? || return_kwh.nil?)
|
|
|
|
usage_kwh_cost, return_kwh_cost = easy_energy_rate(formatted_hour)
|
|
add_tax(formatted_hour, usage_kwh, usage_kwh_cost, return_kwh, return_kwh_cost)
|
|
end
|
|
|
|
def easy_energy_high_low(date)
|
|
hour_start = date.beginning_of_day
|
|
day_end = hour_start.advance(days: 1)
|
|
max_rate = 0.0
|
|
min_rate = 100.0
|
|
min_hour = nil
|
|
max_hour = nil
|
|
while(hour_start < day_end) do
|
|
easy_usage_rate, easy_return_rate = easy_energy_rate(hour_start.strftime("%F-%H"))
|
|
unless easy_usage_rate.nil?
|
|
# determine max
|
|
if easy_usage_rate > max_rate
|
|
max_rate = easy_usage_rate
|
|
max_hour = hour_start.strftime("%F-%H")
|
|
end
|
|
# determine min
|
|
if easy_usage_rate < min_rate
|
|
min_rate = easy_usage_rate
|
|
min_hour = hour_start.strftime("%F-%H")
|
|
end
|
|
end
|
|
hour_start = hour_start.advance(:hours => 1)
|
|
end
|
|
return min_hour,max_hour
|
|
end
|
|
|
|
######################################################
|
|
# Oxxio rates and cost
|
|
######################################################
|
|
|
|
def oxxio_rate(formatted_hour, high_tariff)
|
|
year = Date.parse(formatted_hour).year
|
|
case year
|
|
when 2020
|
|
high_tariff ? 0.07865 : 0.06215
|
|
when 2021
|
|
high_tariff ? 0.06782 : 0.05259
|
|
when 2022
|
|
high_tariff ? 0.23665 : 0.19408
|
|
end
|
|
end
|
|
|
|
def oxxio_energy_cost(formatted_hour, normaal_kwh, dal_kwh, year_shift=0)
|
|
return nil if (normaal_kwh.nil? || dal_kwh.nil?)
|
|
#year = Date.parse(formatted_hour).year+year_shift
|
|
date = Date.parse(formatted_hour).advance(years: year_shift)
|
|
case date.to_time.to_i
|
|
# Date.parse("2019-12-08").to_time.to_i
|
|
# From 8 Dec 2019 until 7 Dec 2020
|
|
when 1575759600..1607295600
|
|
#normaal_kwh * (0.07865 + 0.11822 + 0.03303) + dal_kwh * (0.06215 + 0.11822 + 0.03303)
|
|
normaal_kwh_cost = 0.07865
|
|
dal_kwh_cost = 0.06215
|
|
# From 8 Dec 2020 until 7 Dec 2021
|
|
when 1607382000..1638831600
|
|
#normaal_kwh * (0.06782 + 0.11408 + 0.03630) + dal_kwh * (0.05259 + 0.11408 + 0.03630)
|
|
normaal_kwh_cost = 0.06782
|
|
dal_kwh_cost = 0.05259
|
|
# From 8 Dec 2021 until 7 Sept 2022
|
|
when 1638918000..1662501600
|
|
# normaal_kwh * (0.23665 + 0.04452 + 0.03691) + dal_kwh * (0.19408 + 0.04452 + 0.03691)
|
|
normaal_kwh_cost = 0.23665
|
|
dal_kwh_cost = 0.19408
|
|
# From 8 Sept 2022 until 7 Dec 2022
|
|
when 1662588000..1670367600
|
|
normaal_kwh_cost = 0.60824
|
|
dal_kwh_cost = 0.43701
|
|
end
|
|
normaal_cost = add_tax(formatted_hour, normaal_kwh,normaal_kwh_cost,0,0) # return_kwh already accounted for
|
|
dal_cost = add_tax(formatted_hour, dal_kwh, dal_kwh_cost,0,0)
|
|
# result
|
|
normaal_cost + dal_cost
|
|
end
|
|
|
|
#
|
|
# Entsoe
|
|
#
|
|
|
|
def entsoe_energy_cost(formatted_hour, usage_kwh, return_kwh)
|
|
return nil if (usage_kwh.nil? || return_kwh.nil?)
|
|
|
|
usage_kwh_cost = return_kwh_cost = entsoe.price_at(formatted_hour)*(1+vat_at(Date.parse(formatted_hour)))
|
|
add_tax(formatted_hour, usage_kwh, usage_kwh_cost, return_kwh, return_kwh_cost)
|
|
end
|
|
|
|
######################################################
|
|
# Calculate the per_hour usage and costs
|
|
######################################################
|
|
|
|
def hours(date, year_shift=0)
|
|
hour_start = date.in_time_zone("Amsterdam").beginning_of_day
|
|
day_end = hour_start.advance(days: 1)
|
|
result = []
|
|
lowest_hour, highest_hour,high_hours = entsoe.high_low_hours(date)
|
|
|
|
while(hour_start < day_end) do
|
|
hour_end = hour_start.end_of_hour
|
|
formatted_hour = hour_start.strftime("%F %H")
|
|
|
|
#p "Fetching meter readings between %s and %s" % [hour_start,hour_end]
|
|
hour_readings = Reading.where("created_at > :begin AND created_at < :end", {:begin => hour_start, :end => hour_end})
|
|
hour_diff = hour_readings.last ? hour_readings.last.diff(hour_readings.first) : UNKNOWN_READING
|
|
# calculate cost of this hour
|
|
usage_kwh = hour_diff[:total_kwh_consumed_high] + hour_diff[:total_kwh_consumed_low] rescue nil
|
|
return_kwh = hour_diff[:total_kwh_produced_high] + hour_diff[:total_kwh_produced_low] rescue nil
|
|
entsoe_cost = entsoe_energy_cost(formatted_hour, usage_kwh, return_kwh) # without battery use
|
|
|
|
#
|
|
# Make battery work
|
|
#
|
|
if !usage_kwh.nil?
|
|
# charge battery with return_kwh
|
|
return_kwh -= battery.charge(return_kwh)
|
|
|
|
if (lowest_hour.eql?(formatted_hour) || entsoe.price_at(formatted_hour) < 0) # lowest_hour = "" if small difference between high/low
|
|
# charge battery during lowest hour, or when prices are negative
|
|
usage_kwh += battery.charge(max_charge_kwh)
|
|
else
|
|
# if during expensive hours || more than <max_charge_kwh> kwh then discharge_battery
|
|
# if (battery.battery_kwh > max_charge_kwh) || high_hours.include?(formatted_hour)
|
|
if high_hours.include?(formatted_hour)
|
|
usage_kwh -= battery.discharge(usage_kwh)
|
|
end
|
|
end
|
|
end
|
|
entsoe_cost_with_battery = entsoe_energy_cost(formatted_hour, usage_kwh, return_kwh) # with battery use
|
|
|
|
#
|
|
# end of battery work
|
|
#
|
|
|
|
normal_usage_kwh = hour_diff[:total_kwh_consumed_low] - hour_diff[:total_kwh_produced_low] rescue nil
|
|
low_usage_kwh = hour_diff[:total_kwh_consumed_high] - hour_diff[:total_kwh_produced_high] rescue nil
|
|
oxxio_cost = oxxio_energy_cost(formatted_hour,normal_usage_kwh, low_usage_kwh, year_shift)
|
|
oxxio_rate = oxxio_rate(formatted_hour, (hour_readings.first.high_tarif rescue true))
|
|
|
|
one_hour = [formatted_hour,
|
|
hour_diff[:total_kwh_consumed_high],
|
|
hour_diff[:total_kwh_consumed_low],
|
|
hour_diff[:total_kwh_produced_high],
|
|
hour_diff[:total_kwh_produced_low],
|
|
usage_kwh,
|
|
return_kwh,
|
|
battery.battery_kwh,
|
|
entsoe.price_at(formatted_hour), # entsoe rate
|
|
entsoe_cost_with_battery,
|
|
entsoe_cost,
|
|
oxxio_rate,
|
|
oxxio_cost]
|
|
|
|
p "%s,%s,%s,%s,%s,%s,%s,%0.1f,%s,%s,%s,%s,%s" % (one_hour[0..7] + one_hour[8..12].map{|c| format_cost(c)})
|
|
result << one_hour
|
|
|
|
hour_start = hour_start.advance(:hours => 1)
|
|
end
|
|
result
|
|
end
|
|
|
|
def summarize_hours(hours)
|
|
usage_kwh = hours.map{|e| (e[1] && e[2]) ? (e[1] + e[2]) : nil}.compact.sum
|
|
return_kwh = hours.map{|e| (e[3] && e[4]) ? (e[3] + e[4]) : nil}.compact.sum
|
|
usage_kwh_with_battery = hours.map{|e| e[5]}.compact.sum
|
|
return_kwh_with_battery = hours.map{|e| e[6]}.compact.sum
|
|
entsoe_cost_with_battery = hours.map{|e| e[9]}.compact.sum
|
|
entsoe_cost = hours.map{|e| e[10]}.compact.sum
|
|
oxxio_cost = hours.map{|e| e[12]}.compact.sum
|
|
|
|
return usage_kwh, return_kwh, usage_kwh_with_battery, return_kwh_with_battery, format_cost(entsoe_cost_with_battery), format_cost(entsoe_cost), format_cost(oxxio_cost)
|
|
end
|
|
|
|
def hours_in_month(date, year_shift=0)
|
|
day = date.beginning_of_month
|
|
last_day = date.end_of_month
|
|
all_hours = []
|
|
while (day <= last_day)
|
|
all_hours += hours(day, year_shift)
|
|
day = day.advance(:days => 1)
|
|
end
|
|
all_hours
|
|
end
|
|
|
|
# heel2021 = hours_in(Date.parse("2021-01-01"), Date.parse("2021-12-31"), 1)
|
|
|
|
def hours_in(from,to, year_shift=0)
|
|
day = from.to_date
|
|
last_day = to.to_date
|
|
all_hours = []
|
|
while (day <= last_day)
|
|
all_hours += hours(day, year_shift)
|
|
day = day.advance(:days => 1)
|
|
end
|
|
all_hours
|
|
end
|
|
|
|
def write_csv(filename,hours)
|
|
CSV.open(filename, "wb") do |csv|
|
|
csv << ["hour", "kwh_consumed_low", "kwh_consumed_high", "kwh_produced_low", "kwh_produced_high", "easy_usage_rate", "easy_return_rate", "easy_cost", "oxxio_rate", "oxxio_cost"]
|
|
hours.each do |row|
|
|
csv << row
|
|
end
|
|
end
|
|
end
|
|
|
|
|
|
end
|