From b8e47e262342e6a3cbfef97366c1c359f1fea371 Mon Sep 17 00:00:00 2001 From: Dusan Vojacek Date: Mon, 25 May 2026 03:00:51 +0200 Subject: [PATCH] n --- backend/services/planning_engine.py | 57 ++++++++++++++------ backend/tests/test_planning_dispatch_milp.py | 6 +-- docs/04-modules/planning.md | 2 +- docs/planning-changelog.md | 8 ++- scripts/diagnose_home01_infeasible.py | 4 +- 5 files changed, 55 insertions(+), 22 deletions(-) diff --git a/backend/services/planning_engine.py b/backend/services/planning_engine.py index d0d1299..b881e04 100644 --- a/backend/services/planning_engine.py +++ b/backend/services/planning_engine.py @@ -65,8 +65,8 @@ NEG_SELL_PV_B_VENT_PENALTY_CZK_KWH = 4.0 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 = 250.0 -PLANNER_BUILD_TAG = "2026-05-28-buy-sell-split-v22" +PRE_NEG_CHARGE_PENALTY_CZK_KWH = 400.0 +PLANNER_BUILD_TAG = "2026-05-28-buy-sell-split-v22b" 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 @@ -1061,11 +1061,13 @@ def solve_dispatch( charge_commitment_prev_w: Optional[list[Optional[float]]] = None, planner_version: str | None = None, relaxed_expensive_import: bool = False, + relaxed_neg_buy_charge: bool = False, ) -> tuple[list[DispatchResult], int, dict[str, Any]]: """ LP solver pro dispatch optimalizaci. Vrátí (výsledky, solver_duration_ms, solver_debug_snapshot). relaxed_expensive_import: nouzový režim po Infeasible — síť smí krmit baseload v drahých slotech. + relaxed_neg_buy_charge: druhý nouzový retry bez neg_buy charge shortfall. """ T = len(slots) if T < 1: @@ -1409,12 +1411,15 @@ def solve_dispatch( )) sf = pulp.LpVariable(f"export_shortfall_{t}", 0, cap_w) peak_export_shortfall.append((t, sf, cap_w)) - for t in range(T): - if float(slots[t].buy_price) >= 0.0: - continue - cap_w = float(battery.max_charge_power_w) - sf_nb = pulp.LpVariable(f"neg_buy_charge_sf_{t}", 0, cap_w) - neg_buy_charge_shortfall.append((t, sf_nb, 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 + ] + if neg_buy_slot_indices: + t_nb_last = max(neg_buy_slot_indices) + cap_w = float(battery.max_charge_power_w) + sf_nb = pulp.LpVariable(f"neg_buy_charge_sf_{t_nb_last}", 0, cap_w) + neg_buy_charge_shortfall.append((t_nb_last, sf_nb, cap_w)) for t in range(T): if float(slots[t].sell_price) >= 0: continue @@ -1771,6 +1776,9 @@ def solve_dispatch( if s.sell_price < 0: prob += w_arb[t] == 0 prob += bd[t] <= pulp.lpSum(ev_via_bat[e][t] for e in range(EV)) + # buy<0: export už zakázán výše; neaplikovat sell<0 ventil (bilance / infeasible). + if float(s.buy_price) < 0.0: + continue block_neg_sell_export_t = bool( getattr(grid, "block_export_on_negative_sell", False) ) @@ -1986,14 +1994,7 @@ def solve_dispatch( before_neg_buy = ( first_neg_buy_idx is not None and t < first_neg_buy_idx ) - if ( - before_neg_buy - and float(s.buy_price) >= 0.0 - and bool(getattr(s, "is_daytime_pv_surplus_slot", False)) - ): - prob += bc_pv[t] == 0 - prob += bc_gi[t] == 0 - elif before_neg_buy and sell_t_pre < 0.0 and pv_surplus_w > 0: + 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). prob += bc_pv[t] == 0 elif t not in charge_slots: @@ -2040,6 +2041,10 @@ def solve_dispatch( and pre_neg_export_last_t is not None and t <= pre_neg_export_last_t and sell_t >= 0 + and ( + first_neg_buy_idx is None + or t < first_neg_buy_idx + ) ) pv_store_val = _pv_store_value_czk_kwh(s, min_spread) skip_pv_store_block = ( @@ -2189,6 +2194,26 @@ def solve_dispatch( planner_version=planner_version, relaxed_expensive_import=True, ) + if not relaxed_neg_buy_charge: + logger.warning( + "solve_dispatch still Infeasible, retry without neg_buy_charge_shortfall" + ) + return solve_dispatch( + slots, + battery, + heat_pump, + grid, + ev_sessions, + vehicles, + current_soc_wh, + current_tuv_temp_c, + tuv_delta_stats=tuv_delta_stats, + operating_mode=operating_mode, + charge_commitment_prev_w=charge_commitment_prev_w, + planner_version=planner_version, + relaxed_expensive_import=True, + relaxed_neg_buy_charge=True, + ) raise RuntimeError(f"Solver: {pulp.LpStatus[status]}") # --- Post-processing --- diff --git a/backend/tests/test_planning_dispatch_milp.py b/backend/tests/test_planning_dispatch_milp.py index da2b05e..24b1000 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-v22") + self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-buy-sell-split-v22b") 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-v22") + self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-buy-sell-split-v22b") 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-v22") + self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-buy-sell-split-v22b") 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 6636281..291dc6a 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 v22:** 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` při `buy≥0` žádné nabíjení (`bc_pv=bc_gi=0`). Viz changelog v22. + **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. 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 f721c79..5e64418 100644 --- a/docs/planning-changelog.md +++ b/docs/planning-changelog.md @@ -5,7 +5,13 @@ Formát: **datum (ISO)** · stručný důvod · soubory · chování / ověřen --- -## 2026-05-28 — rozlišení buy<0 vs sell<0 (v22) +## 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`. + +--- + +## 2026-05-28 — rozlišení buy<0 vs sell<0 (v22, superseded by v22b) **Problém (MCP run 16706, v21b):** Znaménka v objective OK (`grid<0` = export, `bat>0` = nabíjení). Chování ale „opačně“: - **Před buy<0** (05:30–07:00, buy≥0): nabíjení z PV/sítě místo přípravy kapacity. diff --git a/scripts/diagnose_home01_infeasible.py b/scripts/diagnose_home01_infeasible.py index 7de9583..61b0194 100644 --- a/scripts/diagnose_home01_infeasible.py +++ b/scripts/diagnose_home01_infeasible.py @@ -13,7 +13,9 @@ sys.path.insert(0, str(Path(__file__).resolve().parents[1] / "backend")) from services.planning_engine import PlanningSlot, solve_dispatch, solve_dispatch_two_pass, PLANNER_BUILD_TAG # Export z MCP: planning_interval run_id=16674 + fn_planning_site_context(2) -SLOTS_JSON = Path(__file__).with_name("home01_run16674_slots.json") +SLOTS_JSON = Path(__file__).with_name("home01_run16706_slots.json") +if not SLOTS_JSON.exists(): + SLOTS_JSON = Path(__file__).with_name("home01_run16674_slots.json") SOC_WH = 37120.0