feat(planner): EV účtování v2 — headroom fix, deadline boundary, min. výkon WB, via-bat reporting
Hloubková diagnóza EV potvrdila: oportunitní ekonomika via-baterie je v LP
správně, ale okraje lhaly nebo byly nevykonatelné:
- V099 + R__039: ems.ev_session.opportunistic_value_czk_kwh (NULL = zdědit
z asset_vehicle, 0 = vypnout pro session); headroom_wh z max(target_soc,
soc_at_connect) — „nenabíjet" (nízký target) už paradoxně NEzvětšuje
oportunistickou vrstvu; vehicles JSON nese min_power_w wallboxu.
- R__015: patch klíč opportunistic_value_czk_kwh (validace >= 0).
- solver_v2: (a) deadline suma range(t_dl) — slot začínající v deadline už
nepatří „do deadline"; (b) Σ ev_direct <= gi + PV (fyzikální split);
(c) binárka ev_on → setpoint ∈ {0} ∪ [min_power_w, max] (konec 400–900 W
nevykonatelných setpointů); (d) bez session EV == 0 (stop-session i golden
fixtures — žádné pumpování při buy<0); dekompozice total == needed − unmet
+ opp i pro needed = 0; (e) battery_arbitrage_czk = via_bat kWh × oportunitní
cena (min sell exportního slotu téhož pražského dne, jinak terminal value)
místo konstantní 0. Oportunismus PO deadline zůstává POVOLENÝ (rozhodnutí:
auto často doma, odjezd řeší rolling replan).
- R__033: fn_plan_current_bundle.intervals + ev1/ev2_via_bat_w (UI nemá cenit
EV kWh z baterie slotovým buy).
Golden gate beze změny snapshotů (v1 nedotčen, fixtures bez EV sessions);
solver_v2_eval před/po identický (CELKEM −1283.5 Kč, Δ −221.9 vs v1);
tests/test_solver_v2.py +7 testů; plná sada 310 passed / 4 xfailed.
Co-Authored-By: Claude Opus 4.8 (1M context) <noreply@anthropic.com>
This commit is contained in:
Dusan Vojacek
@@ -28,9 +28,17 @@
|
||||
# vede k "nabít plným výkonem hned, pak řezat A" — emergentně, bez rampy.
|
||||
# - oportunistické EV („měkký cíl"): nad tvrdý target smí auto vzít až
|
||||
# headroom_wh (do 100 %), oceněno opportunistic_value_czk_kwh (= budoucí
|
||||
# ušetřené nabíjení, DB) — kupuje jen velmi levnou/zápornou energii.
|
||||
# Dekompozice Σ(EV energie) == needed − unmet + opp zároveň stropuje
|
||||
# celkovou energii do auta (dřív při buy<0 bez stropu).
|
||||
# ušetřené nabíjení, session override → vozidlo, DB) — kupuje jen velmi
|
||||
# levnou/zápornou energii. Dekompozice Σ(EV energie) == needed − unmet + opp
|
||||
# zároveň stropuje celkovou energii do auta (dřív při buy<0 bez stropu);
|
||||
# opp vrstva NENÍ vázaná deadline (auto bývá doma dál, odjezd řeší rolling
|
||||
# replan); bez session je EV == 0 (stop-session). Deadline suma jde po
|
||||
# slot PŘED deadline (slot začínající v deadline už nepatří „do deadline").
|
||||
# - min. výkon wallboxu (asset_ev_charger.min_power_w, 6 A ≈ 1380 W):
|
||||
# binárka ev_on → setpoint ∈ {0} ∪ [min_power_w, max]; ev_direct ≤ gi + PV
|
||||
# (fyzikální split direct/via_bat). Reporting: kWh přes ev_via_bat plní
|
||||
# battery_arbitrage_czk oportunitní cenou (min sell exportního slotu dne,
|
||||
# jinak terminal value) — slotový buy pro ně neplatí.
|
||||
# - denní SoC rampa: deficit pod slot.safety_soc_target_wh (R__063: reserve →
|
||||
# reserve+noc, 6–19 h) platí za slot nájem buy×faktor (DB
|
||||
# planner_safety_soc_risk_factor) — ráno se nejdřív dotáhne rezerva
|
||||
@@ -58,13 +66,14 @@ from services.planning.constants import (
|
||||
from services.planning.types import (
|
||||
DispatchResult,
|
||||
PlanningSlot,
|
||||
_prague_calendar_date,
|
||||
_prague_dow_hour,
|
||||
)
|
||||
from services.planning.heuristics import _dispatch_grid_setpoint_w
|
||||
|
||||
logger = logging.getLogger(__name__)
|
||||
|
||||
V2_BUILD_TAG = "v2-clean-2026-06-11"
|
||||
V2_BUILD_TAG = "v2-ev-accounting-2026-06-12"
|
||||
|
||||
# Cena za vypnutí GEN portu (mikroinvertory pole B): reálné riziko/opotřebení
|
||||
# cyklování stykače — drobná, ale nenulová, aby cutoff platil jen při sell < 0.
|
||||
@@ -166,6 +175,10 @@ def solve_dispatch_v2(
|
||||
]
|
||||
ev_unmet: list = [] # slack Wh per session (cena V2_EV_UNMET_CZK_KWH)
|
||||
ev_opp: list = [] # (var, value_czk_kwh) — energie nad target (měkký cíl)
|
||||
# min. výkon wallboxu (IEC 61851: 6 A ≈ 1380 W) — setpoint ∈ {0} ∪ [min, max]
|
||||
ev_min_w = [
|
||||
max(0.0, float(getattr(vehicles[e], "min_power_w", 0) or 0)) for e in range(EV)
|
||||
]
|
||||
nb_buffer_wh = [max(0.0, float(s.night_baseload_buffer_wh or 0.0)) for s in slots]
|
||||
safety_risk = float(getattr(battery, "planner_safety_soc_risk_factor", 0.0) or 0.0)
|
||||
safety_tgt_wh = [
|
||||
@@ -222,6 +235,11 @@ def solve_dispatch_v2(
|
||||
prob += bc_gi[t] <= gi[t], f"bcgi_src_{t}"
|
||||
# vybíjení kryje dům + EV-via-bat + export z baterie
|
||||
prob += ge_bat[t] + pulp.lpSum(ev_via_bat[e][t] for e in range(EV)) <= bd[t], f"bd_split_{t}"
|
||||
# ev_direct fyzicky jen ze sítě + PV (ne z baterie) — split direct/via_bat
|
||||
# není arbitrární, ekonomiku nemění (bilance platí stejně)
|
||||
prob += (
|
||||
pulp.lpSum(ev_direct[e][t] for e in range(EV)) <= gi[t] + pv_a_net + pv_b_eff
|
||||
), f"evd_src_{t}"
|
||||
|
||||
# zákaz současného importu a exportu
|
||||
prob += gi[t] <= max_imp * y_imp[t], f"imp_excl_{t}"
|
||||
@@ -245,15 +263,20 @@ def solve_dispatch_v2(
|
||||
if float(s.sell_price) < 0.0 and block_neg_sell:
|
||||
prob += ge_pv[t] + ge_bat[t] == 0, f"neg_sell_block_{t}"
|
||||
|
||||
# EV dostupnost
|
||||
# EV dostupnost + min. výkon wallboxu (binárka jen kde je min > 0)
|
||||
for e in range(EV):
|
||||
if not _connected(e, t):
|
||||
prob += ev_direct[e][t] == 0
|
||||
prob += ev_via_bat[e][t] == 0
|
||||
else:
|
||||
prob += ev_direct[e][t] + ev_via_bat[e][t] <= float(
|
||||
vehicles[e].max_charge_power_w
|
||||
)
|
||||
ev_max_w = float(vehicles[e].max_charge_power_w)
|
||||
ev_total = ev_direct[e][t] + ev_via_bat[e][t]
|
||||
if 0 < ev_min_w[e] <= ev_max_w:
|
||||
on = pulp.LpVariable(f"evon_{e}_{t}", cat=pulp.LpBinary)
|
||||
prob += ev_total >= ev_min_w[e] * on, f"ev_min_{e}_{t}"
|
||||
prob += ev_total <= ev_max_w * on, f"ev_max_{e}_{t}"
|
||||
else:
|
||||
prob += ev_total <= ev_max_w
|
||||
|
||||
# provozní režimy (tvrdé constraints dle operating-modes.md)
|
||||
if om == "SELF_SUSTAIN":
|
||||
@@ -266,28 +289,40 @@ def solve_dispatch_v2(
|
||||
prob += ge_pv[t] + ge_bat[t] == 0
|
||||
prob += bd[t] == 0
|
||||
|
||||
# EV deadline (s placeným slackem místo infeasibility)
|
||||
# EV deadline (s placeným slackem místo infeasibility) + měkký cíl.
|
||||
# Bez session není mandát nabíjet: připojené auto bez session (stop-session,
|
||||
# golden fixtures s vynulovanými sessions) nesmí při buy<0 „pumpovat" energii.
|
||||
for e in range(EV):
|
||||
sess = ev_sessions[e] if e < len(ev_sessions) else None
|
||||
if sess is None or not getattr(sess, "energy_needed_wh", 0):
|
||||
if sess is None:
|
||||
for t in range(T):
|
||||
if _connected(e, t):
|
||||
prob += ev_direct[e][t] == 0, f"ev_nosess_d_{e}_{t}"
|
||||
prob += ev_via_bat[e][t] == 0, f"ev_nosess_b_{e}_{t}"
|
||||
continue
|
||||
t_dl = next(
|
||||
(t for t in range(T) if slots[t].interval_start >= sess.target_deadline),
|
||||
T - 1,
|
||||
)
|
||||
unmet = pulp.LpVariable(f"ev_unmet_{e}", 0, float(sess.energy_needed_wh))
|
||||
needed = max(0.0, float(getattr(sess, "energy_needed_wh", 0.0) or 0.0))
|
||||
unmet = pulp.LpVariable(f"ev_unmet_{e}", 0, needed)
|
||||
ev_unmet.append(unmet)
|
||||
prob += (
|
||||
pulp.lpSum(
|
||||
(ev_direct[e][t] + ev_via_bat[e][t]) * INTERVAL_H
|
||||
for t in range(t_dl + 1)
|
||||
if _connected(e, t)
|
||||
if needed > 0:
|
||||
# první slot s interval_start >= deadline už do deadline NEPATŘÍ
|
||||
# (slot [deadline, deadline+15min) dodává energii až po odjezdu)
|
||||
t_dl = next(
|
||||
(t for t in range(T) if slots[t].interval_start >= sess.target_deadline),
|
||||
T,
|
||||
)
|
||||
+ unmet
|
||||
>= float(sess.energy_needed_wh)
|
||||
), f"ev_deadline_{e}"
|
||||
prob += (
|
||||
pulp.lpSum(
|
||||
(ev_direct[e][t] + ev_via_bat[e][t]) * INTERVAL_H
|
||||
for t in range(t_dl)
|
||||
if _connected(e, t)
|
||||
)
|
||||
+ unmet
|
||||
>= needed
|
||||
), f"ev_deadline_{e}"
|
||||
|
||||
# měkký cíl: dekompozice celkové energie == needed − unmet + opp
|
||||
# měkký cíl: dekompozice celkové energie == needed − unmet + opp.
|
||||
# Oportunistická vrstva NENÍ omezená deadline — auto bývá doma dál,
|
||||
# odjezd řeší rolling replan (rozhodnutí 2026-06-12).
|
||||
headroom = max(0.0, float(getattr(sess, "headroom_wh", 0.0) or 0.0))
|
||||
opp_val = float(getattr(sess, "opportunistic_value_czk_kwh", 0.0) or 0.0)
|
||||
opp = pulp.LpVariable(f"ev_opp_{e}", 0, headroom if opp_val > 0 else 0.0)
|
||||
@@ -298,7 +333,7 @@ def solve_dispatch_v2(
|
||||
for t in range(T)
|
||||
if _connected(e, t)
|
||||
)
|
||||
== float(sess.energy_needed_wh) - unmet + opp
|
||||
== needed - unmet + opp
|
||||
), f"ev_total_{e}"
|
||||
|
||||
# TUV look-ahead (převzato z v1 — komfortní constraint, ne heuristika)
|
||||
@@ -384,9 +419,32 @@ def solve_dispatch_v2(
|
||||
v = pulp.value(var)
|
||||
return float(v) if v is not None else 0.0
|
||||
|
||||
# Reporting EV-via-bat: kWh do auta z baterie neplatí slotový buy (jdou
|
||||
# z baterie), ale ušlou příležitost. Aproximace oportunitní ceny: nejnižší
|
||||
# sell slotu, kde plán exportuje, v témže pražském dni; bez exportu ten den
|
||||
# terminal value (Kč/kWh). Plní battery_arbitrage_czk (dřív konstantní 0).
|
||||
day_min_export_sell: dict[Any, float] = {}
|
||||
for t in range(T):
|
||||
if _val(ge_pv[t]) + _val(ge_bat[t]) >= 1.0:
|
||||
d_key = _prague_calendar_date(slots[t])
|
||||
sp = float(slots[t].sell_price)
|
||||
if d_key not in day_min_export_sell or sp < day_min_export_sell[d_key]:
|
||||
day_min_export_sell[d_key] = sp
|
||||
|
||||
results: list[DispatchResult] = []
|
||||
for t in range(T):
|
||||
s = slots[t]
|
||||
via1_w = _val(ev_via_bat[0][t]) if EV > 0 else 0.0
|
||||
via2_w = _val(ev_via_bat[1][t]) if EV > 1 else 0.0
|
||||
via_kwh = (via1_w + via2_w) * wh
|
||||
if via_kwh > 1e-9:
|
||||
opp_price = max(
|
||||
0.0,
|
||||
day_min_export_sell.get(_prague_calendar_date(s), terminal * 1000.0),
|
||||
)
|
||||
arb_czk = via_kwh * opp_price
|
||||
else:
|
||||
arb_czk = 0.0
|
||||
bc_tot = _val(bc_pv[t]) + _val(bc_gi[t])
|
||||
bd_v = _val(bd[t])
|
||||
batt_w = round(bc_tot - bd_v)
|
||||
@@ -434,8 +492,8 @@ def solve_dispatch_v2(
|
||||
if EV > 1 and s.ev2_connected
|
||||
else None
|
||||
),
|
||||
ev1_via_bat_w=round(_val(ev_via_bat[0][t])) if EV > 0 else 0,
|
||||
ev2_via_bat_w=round(_val(ev_via_bat[1][t])) if EV > 1 else 0,
|
||||
ev1_via_bat_w=round(via1_w),
|
||||
ev2_via_bat_w=round(via2_w),
|
||||
heat_pump_enabled=hp_on,
|
||||
heat_pump_setpoint_w=int(rated_hp) if hp_on else 0,
|
||||
pv_a_curtailed_w=round(_val(ca[t])),
|
||||
@@ -444,7 +502,7 @@ def solve_dispatch_v2(
|
||||
effective_sell_price=float(s.sell_price),
|
||||
is_predicted_price=bool(s.is_predicted_price),
|
||||
cashflow_czk=round(cash_t, 4),
|
||||
battery_arbitrage_czk=0.0,
|
||||
battery_arbitrage_czk=round(arb_czk, 4),
|
||||
penalty_czk=round(pen_t, 4),
|
||||
green_bonus_czk=float(getattr(s, "green_bonus_czk_per_slot", 0.0) or 0.0),
|
||||
)
|
||||
@@ -462,6 +520,7 @@ def solve_dispatch_v2(
|
||||
"gen_cutoff_available": gen_cutoff_avail,
|
||||
"slot_count": T,
|
||||
"ev_sessions": sum(1 for x in ev_sessions if x is not None),
|
||||
"ev_min_power_w": ev_min_w,
|
||||
"masks_ignored": True,
|
||||
"night_buffer_slots": sum(1 for b in nb_buffer_wh if b > 0),
|
||||
"pv_risk_frontload_czk_kwh": frontload if neg_idx else 0.0,
|
||||
|
||||
Reference in New Issue
Block a user