From 37a525cb4f78b3e6c27537bb866dd656f3fc720a Mon Sep 17 00:00:00 2001 From: Dusan Vojacek Date: Mon, 25 May 2026 03:09:33 +0200 Subject: [PATCH] predvybiti baterky --- backend/services/planning_engine.py | 66 +++++++++++++++++++- backend/tests/test_planning_dispatch_milp.py | 6 +- docs/04-modules/planning.md | 2 +- docs/planning-changelog.md | 8 +++ 4 files changed, 76 insertions(+), 6 deletions(-) diff --git a/backend/services/planning_engine.py b/backend/services/planning_engine.py index b881e04..a5488d3 100644 --- a/backend/services/planning_engine.py +++ b/backend/services/planning_engine.py @@ -66,7 +66,9 @@ EXTREME_BUY_DUMP_PREWINDOW_SLOTS = 12 NEG_SELL_BAT_DUMP_SHORTFALL_PENALTY_CZK_KWH = 80.0 NEG_BUY_CHARGE_SHORTFALL_PENALTY_CZK_KWH = 100.0 PRE_NEG_CHARGE_PENALTY_CZK_KWH = 400.0 -PLANNER_BUILD_TAG = "2026-05-28-buy-sell-split-v22b" +PRE_NEG_BATT_EXPORT_SHORTFALL_PENALTY_CZK_KWH = 80.0 +PRE_NEG_BATT_EXPORT_MIN_SELL_CZK_KWH = 1.0 +PLANNER_BUILD_TAG = "2026-05-28-pre-neg-batt-discharge-v23" CORRECTION_WINDOW_H = 1 # hodina zpět pro výpočet korekčního faktoru CORRECTION_MIN_CLAMP = 0.5 # spodní limit korekčního faktoru CORRECTION_MAX_CLAMP = 1.5 # horní limit korekčního faktoru @@ -877,6 +879,38 @@ def _morning_pre_neg_export_indices( return out +def _pre_neg_buy_discharge_indices( + slots: list[PlanningSlot], + first_neg_buy_idx: int | None, + *, + charge_acquisition_czk_kwh: float, + min_spread: float, + fixed_tariff: bool, +) -> set[int]: + """ + Sloty před prvním buy<0: výboj baterie do sítě při kladném sell (včetně noci). + Bez rozšíření discharge_export_slots (v19b — jinak w_arb → Infeasible). + """ + if first_neg_buy_idx is None or first_neg_buy_idx <= 0: + return set() + out: set[int] = set() + for i in range(first_neg_buy_idx): + s = slots[i] + if float(s.buy_price) < 0.0: + continue + if float(s.sell_price) < PRE_NEG_BATT_EXPORT_MIN_SELL_CZK_KWH: + continue + if not _slot_profitable_battery_export( + s, + charge_acquisition_czk_kwh=charge_acquisition_czk_kwh, + min_spread=min_spread, + fixed_tariff=fixed_tariff, + ): + continue + out.add(i) + return out + + def _evening_peak_export_indices( slots: list[PlanningSlot], *, @@ -1295,6 +1329,15 @@ def solve_dispatch( min_spread_pre = float(degradation_cost_effective) purchase_fixed_pre = _purchase_pricing_fixed(grid) fixed_tariff_like_pre = purchase_fixed_pre or _horizon_fixed_tariff_like(slots) + pre_neg_buy_discharge_ts: set[int] = set() + if om == "AUTO" and first_neg_buy_idx is not None and first_neg_buy_idx > 0: + pre_neg_buy_discharge_ts = _pre_neg_buy_discharge_indices( + slots, + first_neg_buy_idx, + charge_acquisition_czk_kwh=charge_acquisition_czk_kwh, + min_spread=min_spread_pre, + fixed_tariff=fixed_tariff_like_pre, + ) neg_sell_bat_dump_slots = _neg_sell_bat_dump_slots( slots, operating_mode=om, @@ -1392,6 +1435,7 @@ def solve_dispatch( neg_sell_bat_dump_shortfall: list[tuple[int, pulp.LpVariable, float]] = [] neg_sell_soc_underfill: list[tuple[int, pulp.LpVariable]] = [] neg_buy_charge_shortfall: list[tuple[int, pulp.LpVariable, float]] = [] + pre_neg_batt_export_shortfall: list[tuple[int, pulp.LpVariable, float]] = [] fixed_tariff_like = fixed_tariff_like_pre block_export_neg_sell = bool(getattr(grid, "block_export_on_negative_sell", False)) if om == "AUTO": @@ -1411,6 +1455,10 @@ def solve_dispatch( )) sf = pulp.LpVariable(f"export_shortfall_{t}", 0, cap_w) peak_export_shortfall.append((t, sf, cap_w)) + export_cap_w = _battery_export_cap_w(battery, grid) + for t_pnd in sorted(pre_neg_buy_discharge_ts): + sf_pnd = pulp.LpVariable(f"pre_neg_bat_export_sf_{t_pnd}", 0, export_cap_w) + pre_neg_batt_export_shortfall.append((t_pnd, sf_pnd, export_cap_w)) if not relaxed_neg_buy_charge: neg_buy_slot_indices = [ t for t, s in enumerate(slots) if float(s.buy_price) < 0.0 @@ -1597,6 +1645,10 @@ def solve_dispatch( sf * NEG_BUY_CHARGE_SHORTFALL_PENALTY_CZK_KWH * INTERVAL_H / 1000.0 for _t, sf, _cap in neg_buy_charge_shortfall ) + + pulp.lpSum( + sf * PRE_NEG_BATT_EXPORT_SHORTFALL_PENALTY_CZK_KWH * INTERVAL_H / 1000.0 + for _t, sf, _cap in pre_neg_batt_export_shortfall + ) + pulp.lpSum( (bc_pv[t] + bc_gi[t]) * PRE_NEG_CHARGE_PENALTY_CZK_KWH @@ -1627,6 +1679,8 @@ def solve_dispatch( prob += us >= float(battery.soc_max_wh) - soc[t_us] for t_sf, sf, cap_w in neg_buy_charge_shortfall: prob += sf >= cap_w - (bc_gi[t_sf] + bc_pv[t_sf]) + for t_sf, sf, cap_w in pre_neg_batt_export_shortfall: + prob += sf >= cap_w - ge_bat[t_sf] preneg_export_min_soc_wh = float(min_soc_wh) + max( float(battery.max_discharge_power_w) * float(battery.discharge_efficiency) @@ -1639,6 +1693,8 @@ def solve_dispatch( for t_peak in morning_pre_neg_export_ts: if t_peak in profitable_export_ts: prob += ge_bat[t_peak] >= export_push_w * z_export[t_peak] + for t_pnd in pre_neg_buy_discharge_ts: + prob += ge_bat[t_pnd] >= export_push_w * z_export[t_pnd] evening_push_ts = _evening_battery_export_push_indices( slots, profitable_export_ts=profitable_export_ts, @@ -1903,6 +1959,8 @@ def solve_dispatch( and floor_pct is not None ): export_soc_floor_t = float(planner_floor_effective_wh) + elif om == "AUTO" and t in pre_neg_buy_discharge_ts: + export_soc_floor_t = float(min_soc_wh) elif ( om == "AUTO" and t in morning_pre_neg_export_ts @@ -2004,7 +2062,11 @@ def solve_dispatch( prob += bc_pv[t] == 0 else: prob += bc_pv[t] <= float(pv_surplus_w) - if t not in discharge_export_slots and t not in neg_sell_bat_dump_slots: + if ( + t not in discharge_export_slots + and t not in neg_sell_bat_dump_slots + and t not in pre_neg_buy_discharge_ts + ): prob += ge_bat[t] == 0 prob += z_export[t] == 0 diff --git a/backend/tests/test_planning_dispatch_milp.py b/backend/tests/test_planning_dispatch_milp.py index 24b1000..0bd44de 100644 --- a/backend/tests/test_planning_dispatch_milp.py +++ b/backend/tests/test_planning_dispatch_milp.py @@ -1230,7 +1230,7 @@ class NegativeSellPvChargeTests(unittest.TestCase): 50.0, operating_mode="AUTO", ) - self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-buy-sell-split-v22b") + self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-pre-neg-batt-discharge-v23") self.assertGreater( results[0].battery_setpoint_w, 5_500, @@ -1380,7 +1380,7 @@ class NegativeSellPvChargeTests(unittest.TestCase): 50.0, operating_mode="AUTO", ) - self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-buy-sell-split-v22b") + self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-pre-neg-batt-discharge-v23") self.assertEqual(len(results), len(slots)) def test_gen_cutoff_full_soc_neg_sell_with_pv_b_feasible(self) -> None: @@ -1444,7 +1444,7 @@ class NegativeSellPvChargeTests(unittest.TestCase): 55.0, operating_mode="AUTO", ) - self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-buy-sell-split-v22b") + self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-pre-neg-batt-discharge-v23") self.assertEqual(len(results), len(slots)) def test_fixed_tariff_neg_sell_no_grid_export(self) -> None: diff --git a/docs/04-modules/planning.md b/docs/04-modules/planning.md index 291dc6a..010d8f6 100644 --- a/docs/04-modules/planning.md +++ b/docs/04-modules/planning.md @@ -45,7 +45,7 @@ 1. **Globální rozpočet Wh** (`discharge_slot_buffer × exportovatelná kapacita`): sloty podle `sell_price desc`. Před prvním `sell < 0` se z rozpočtu **vynechají** sloty, kde **později tentýž den** existuje `sell` vyšší o více než `degradation` (OTE, ne pevné hodiny 00–04). 2. **Večerní špičky per den:** `sell ≥ max(sell) − degradation` jen pro hodiny **≥ 17** (Prague), ne globální max horizontu (jinak by vyhrála půlnoc 3,7 Kč místo večera). 3. **Ranní pásmo před prvním `sell < 0`:** hodiny **5–11** téhož kalendářního dne — všechny sloty s `sell ≥ lokální_max_ráno − degradation`; ostatní sloty mezi ranním pásmem a prvním `sell < 0` s nižším sell mají export **zakázán** (žádný dump v 07:30 za 2 Kč). **`charge_acquisition`:** vážený `buy` před prvním exportem **téhož dne** jako záporné výkupní okno. - **Planner tag v22b:** oddělení **záporného nákupu** (`buy<0` → import/nabíjení, `ge=0`) a **záporného výkupu** (`sell<0`, `buy≥0` → export pole B / curtail A). Před `buy<0` měkká penalizace nabíjení; při Infeasible druhý retry bez `neg_buy` shortfall. Viz changelog v22b. + **Planner tag v23:** v22b + **výboj baterie do sítě** před `buy<0` (`_pre_neg_buy_discharge_indices`, sell≥1 Kč/kWh, push `ge_bat` z DB limitů). Viz changelog v23. V `solve_dispatch` (AUTO): **`charge_slots`** = `allow_charge` z DB + **`buy < 0`** + všechny sloty **`sell < 0`** s PV přebytkem > 500 W (i bez `block_export_on_negative_sell`, BA81). **`pv_charge_shortfall`** / **`NEG_SELL_CURTAIL_PENALTY`** platí v těchto slotech. Při **`sell < 0`**: safety deficit cílí **`soc_max_wh`** (plný planner strop). Po posledním **`sell < 0`** tentýž den: **`post_neg_pv_topup`** dobije z FVE na `soc_max` před exportem (kladný sell, ne high-sell peak). U **fixního tarifu** s polem B: **`ge_pv ≤ pv_b`** (ne pv_store **`ge_pv = 0`**). Při **`deye_gen_microinverter_cutoff_enabled`**: **`ge == 0` jen** pokud **`block_export_on_negative_sell`** (KV1), ne kvůli samotnému `z_gen_cutoff` (BA81 musí moci exportovat B při plné baterii). Vstupní **`soc_wh`** z telemetrie se před MILP omezí přes **`_planner_soc_for_solver`** (rezerva ~650 Wh pod `soc_max`, jinak Infeasible při 100 % SoC a dlouhém záporném výkupu). **`planner_build_tag`** v `solver_params`. Changelog: [`docs/planning-changelog.md`](../planning-changelog.md). - **Záporná nákupní cena:** - horní mez `grid_import` zahrnuje `load_baseline_w` + nabíjení/EV/TČ (bez nekonečného importu). diff --git a/docs/planning-changelog.md b/docs/planning-changelog.md index 5e64418..de1c7fb 100644 --- a/docs/planning-changelog.md +++ b/docs/planning-changelog.md @@ -5,6 +5,14 @@ Formát: **datum (ISO)** · stručný důvod · soubory · chování / ověřen --- +## 2026-05-28 — noční/ranní výboj baterie před buy<0 (v23) + +**Požadavek:** Před ranním oknem záporných cen **vybít baterii do sítě** (ne jen ~500 W do domu), aby zůstala kapacita na levný import v `buy<0`. + +**Oprava (tag `2026-05-28-pre-neg-batt-discharge-v23`):** `_pre_neg_buy_discharge_indices` — sloty `t < first_neg_buy_idx`, `sell ≥ 1` Kč/kWh, marže exportu z baterie; **`ge_bat`** + shortfall + push na DB export cap, **bez** přidání do `discharge_export_slots` (v19b). Výjimka z `ge_bat=0` v pre-selection; exportní SoC podlaha `min_soc`. + +--- + ## 2026-05-28 — rozlišení buy<0 vs sell<0 (v22 / v22b) **v22b — Infeasible:** Tvrdý `is_daytime_pv_surplus` + `ge_pv=0` z pv_store blokoval export před buy<0. Oprava: jen měkká `PRE_NEG_CHARGE_PENALTY`; u `buy<0` přeskočit sell<0 ventil; `neg_buy` shortfall jen na **posledním** buy<0 slotu; retry `relaxed_neg_buy_charge`. Tag `2026-05-28-buy-sell-split-v22b`.