From 0a0668000bc915c4d74601a6911c9264f350b699 Mon Sep 17 00:00:00 2001 From: Dusan Vojacek Date: Mon, 25 May 2026 02:41:36 +0200 Subject: [PATCH] x --- backend/services/planning_engine.py | 44 ++++++++++++-------- backend/tests/test_planning_dispatch_milp.py | 6 +-- docs/04-modules/planning.md | 2 +- docs/planning-changelog.md | 14 +++---- 4 files changed, 36 insertions(+), 30 deletions(-) diff --git a/backend/services/planning_engine.py b/backend/services/planning_engine.py index 58de671..5d0a9ee 100644 --- a/backend/services/planning_engine.py +++ b/backend/services/planning_engine.py @@ -64,7 +64,7 @@ NEG_SELL_PV_B_VENT_PENALTY_CZK_KWH = 4.0 # Výboj baterie při sell<0 jen těsně před extrémně záporným buy (round-trip arbitráž). EXTREME_BUY_DUMP_PREWINDOW_SLOTS = 12 NEG_SELL_BAT_DUMP_SHORTFALL_PENALTY_CZK_KWH = 80.0 -PLANNER_BUILD_TAG = "2026-05-28-morning-hold-soc-v21" +PLANNER_BUILD_TAG = "2026-05-28-morning-hold-soc-v21b" 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 @@ -1262,6 +1262,10 @@ def solve_dispatch( # Kotva: poslední slot před prvním sell<0 by měl končit u planner floor (pokud relaxace existuje). # Slack penalizujeme v objective; samotné omezení přidáme až po definici soc. first_neg_sell_idx, pre_neg_export_last_t = _pre_negative_sell_export_window(slots) + first_neg_buy_idx = next( + (t for t, s in enumerate(slots) if float(s.buy_price) < 0.0), + None, + ) last_neg_sell_by_prague_date: dict[object, int] = {} for t_ln, st_ln in enumerate(slots): if float(st_ln.sell_price) < 0: @@ -1402,16 +1406,13 @@ def solve_dispatch( )) sf = pulp.LpVariable(f"export_shortfall_{t}", 0, cap_w) peak_export_shortfall.append((t, sf, cap_w)) - first_neg_buy_idx_pen = next( - (t for t, s in enumerate(slots) if float(s.buy_price) < 0.0), None - ) for t in range(T): if float(slots[t].sell_price) >= 0: continue if t not in charge_slots: continue # Před buy<0: nepenalizovat / netlačit PV→bat (jinak 98 % v 09:15 a export v sell<0). - if first_neg_buy_idx_pen is not None and t < first_neg_buy_idx_pen: + if first_neg_buy_idx is not None and t < first_neg_buy_idx: continue pv_surplus_w = max( 0.0, @@ -1445,7 +1446,7 @@ def solve_dispatch( continue if t not in charge_slots: continue - if first_neg_buy_idx_pen is not None and t < first_neg_buy_idx_pen: + if first_neg_buy_idx is not None and t < first_neg_buy_idx: continue pv_surplus_w = max( 0.0, @@ -1697,6 +1698,7 @@ def solve_dispatch( om == "AUTO" and float(s.sell_price) < 0.0 and t in charge_slots + and (first_neg_buy_idx is None or t >= first_neg_buy_idx) ): # Záporný výkup: dobít na planner soc_max (typicky 95–100 %), ne jen SQL safety ~50 %. eff_tgt_s = max(eff_tgt_s, float(battery.soc_max_wh)) @@ -1781,15 +1783,24 @@ def solve_dispatch( prob += ge_pv[t] == 0 elif not purchase_fixed_pre: # Spot (home-01): ge_pv=0 dokud není plná baterie; pak jen ventil pole B (ne celý surplus). - soc_prev_neg = current_soc_wh if t == 0 else soc[t - 1] - w_pv_b_vent = pulp.LpVariable(f"w_pv_b_vent_neg_{t}", cat=pulp.LpBinary) - m_soc_neg = float(battery.soc_max_wh) - prob += soc_prev_neg >= ( - m_soc_neg - - soc_headroom_wh - - m_soc_neg * (1 - w_pv_b_vent) + # Před buy<0 + bc_pv=0: přebytek pole B musí jít do sítě (ge_pv≤pv_b), jinak Infeasible. + before_first_neg_buy = ( + first_neg_buy_idx is not None and t < first_neg_buy_idx ) - prob += ge_pv[t] <= float(s.pv_b_forecast_w) * w_pv_b_vent + if before_first_neg_buy and float(s.pv_b_forecast_w) > 0: + prob += ge_pv[t] <= float(s.pv_b_forecast_w) + else: + soc_prev_neg = current_soc_wh if t == 0 else soc[t - 1] + w_pv_b_vent = pulp.LpVariable( + f"w_pv_b_vent_neg_{t}", cat=pulp.LpBinary + ) + m_soc_neg = float(battery.soc_max_wh) + prob += soc_prev_neg >= ( + m_soc_neg + - soc_headroom_wh + - m_soc_neg * (1 - w_pv_b_vent) + ) + prob += ge_pv[t] <= float(s.pv_b_forecast_w) * w_pv_b_vent soc_prev_expr = current_soc_wh if t == 0 else soc[t - 1] arb_t = arb_floor_series[t] @@ -1920,9 +1931,6 @@ def solve_dispatch( prob += bd[t] == 0 # Slot pre-selection (z DB fn_load_planning_slots_full → allow_*) - first_neg_buy_idx_main = next( - (t for t, s in enumerate(slots) if float(s.buy_price) < 0.0), None - ) if om == "AUTO": for t in range(T): s = slots[t] @@ -1942,7 +1950,7 @@ def solve_dispatch( ): prob += bc_gi[t] == 0 before_neg_buy = ( - first_neg_buy_idx_main is not None and t < first_neg_buy_idx_main + first_neg_buy_idx is not None and t < first_neg_buy_idx ) if before_neg_buy and sell_t_pre < 0.0 and pv_surplus_w > 0: # Ranní sell<0 před buy<0: PV do sítě/curtail, ne do baterie (kapacita na import). diff --git a/backend/tests/test_planning_dispatch_milp.py b/backend/tests/test_planning_dispatch_milp.py index 827e784..84ca410 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-morning-hold-soc-v21") + self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-morning-hold-soc-v21b") 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-morning-hold-soc-v21") + self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-morning-hold-soc-v21b") 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-morning-hold-soc-v21") + self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-morning-hold-soc-v21b") 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 650f50e..82da683 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 v21:** po revertu v20 znovu drží SoC před prvním `buy<0` — bez `neg_sell_soc_underfill` / `pv_charge_shortfall` a bez `bc_pv` z FVE v ranním `sell<0` okně před `buy<0` (MCP run 16692: v20 nabíjelo od 05:30 na 98 %, pak export v sell<0). Stále **v17** + **v18**; noční výboj / import v `buy<0` → **`R__063`**. Viz changelog v21. + **Planner tag v21b:** drží SoC před prvním `buy<0` (`bc_pv=0`, shortfally jen od `buy<0`); před `buy<0` v sell<0 povolen export **pole B** (`ge_pv≤pv_b` bez ventilu plné baterie — jinak Infeasible). Stále **v17** + **v18**. Viz changelog v21/v21b. 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 21d8442..ec27928 100644 --- a/docs/planning-changelog.md +++ b/docs/planning-changelog.md @@ -5,19 +5,17 @@ Formát: **datum (ISO)** · stručný důvod · soubory · chování / ověřen --- -## 2026-05-28 — ranní sell<0: držet SoC před buy<0 (v21) +## 2026-05-28 — ranní sell<0: držet SoC před buy<0 (v21 / v21b) -**Problém (MCP run 16692, tag v20):** Od ~05:30 nabíjení z PV; v 09:15 už **98,3 %** SoC; od 09:15 masivní **export při sell<0** (−7 kW). V **11:00–12:45** `buy<0`, ale baterie plná → **žádný import**. Opak požadovaného chování (ráno prodat FVE při kladném sell, nabíjet až v `buy<0`). +**Problém (MCP run 16692, tag v20):** Od ~05:30 nabíjení z PV; v 09:15 už **98,3 %** SoC; od 09:15 masivní **export při sell<0** (−7 kW). V **11:00–12:45** `buy<0`, ale baterie plná → **žádný import**. -**Příčina v20:** Po revertu v19 zůstaly **`neg_sell_soc_underfill`** a **`pv_charge_shortfall`** pro všechny ranní `sell<0` charge sloty (06:45–10:45) **před** prvním `buy<0` (11:00). LP tlačil PV→bat na soc_max ještě před levným importem. +**v21:** `neg_sell_soc_underfill` / `pv_charge_shortfall` jen od `first_neg_buy_idx`; **`bc_pv=0`** před buy<0 v sell<0. -**Oprava (tag `2026-05-28-morning-hold-soc-v21`):** -- `neg_sell_soc_underfill` a `pv_charge_shortfall` (sell<0) jen pro sloty **`t ≥ first_neg_buy_idx`** -- Tvrdé **`bc_pv[t]=0`** pro `sell<0` + PV surplus **před** `first_neg_buy_idx` (bez rozšíření `discharge_export_slots` — neopakuje v19b Infeasible) +**v21b — Infeasible:** `bc_pv=0` + ventil **`w_pv_b_vent`** (export jen při plné baterii) → přebytek **pole B** (`pv_b` > load) nemá kam (bilance). **Oprava:** před `first_neg_buy_idx` povolit **`ge_pv ≤ pv_b`** bez ventilu; safety `soc_max` u sell<0 charge jen od `first_neg_buy_idx`. -**Bez:** pre-neg discharge setů, per-slot `neg_buy_soc` shortfall, retry vrstev v19. +**Tag:** `2026-05-28-morning-hold-soc-v21b` -**Ověření:** `scripts/diagnose_home01_infeasible.py` (realistic masks OK); po deploy replan home-01 → tag v `solver_params`; MCP: ráno nižší SoC, import v 11:00–12:45. +**Ověření:** `scripts/diagnose_home01_infeasible.py`; replan home-01 → tag v `solver_params`. ---