From 161b463367b2d293719202ac62fb67b1bae369e5 Mon Sep 17 00:00:00 2001 From: Dusan Vojacek Date: Mon, 25 May 2026 01:59:33 +0200 Subject: [PATCH] a dalsi --- backend/services/planning_engine.py | 99 +++++++++++++------- backend/tests/test_planning_dispatch_milp.py | 6 +- docs/04-modules/planning.md | 2 +- docs/05-todo.md | 2 +- docs/planning-changelog.md | 11 +++ 5 files changed, 82 insertions(+), 38 deletions(-) diff --git a/backend/services/planning_engine.py b/backend/services/planning_engine.py index 42605a7..cce2f07 100644 --- a/backend/services/planning_engine.py +++ b/backend/services/planning_engine.py @@ -68,7 +68,7 @@ NEG_SELL_BAT_DUMP_SHORTFALL_PENALTY_CZK_KWH = 80.0 NEG_BUY_GRID_CHARGE_SHORTFALL_PENALTY_CZK_KWH = 120.0 # Měkký tlak: v buy<0 okně dobít na soc_max (ne zastavit na ~94 %). NEG_BUY_SOC_UNDERFILL_PENALTY_CZK_PER_WH = 0.45 -PLANNER_BUILD_TAG = "2026-05-27-pre-neg-buy-strategy-v19b" +PLANNER_BUILD_TAG = "2026-05-27-pre-neg-buy-strategy-v19c" 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 @@ -937,6 +937,16 @@ def _evening_battery_export_push_indices( return sorted(out) +def _negative_buy_slot_indices(slots: list[PlanningSlot]) -> list[int]: + return [t for t, s in enumerate(slots) if float(s.buy_price) < 0.0] + + +def _neg_buy_soc_pressure_slots(slots: list[PlanningSlot]) -> list[int]: + """Tlak na soc_max jen na konci buy<0 okna (ne každý 15min slot zvlášť).""" + neg = _negative_buy_slot_indices(slots) + return [neg[-1]] if neg else [] + + def _pre_neg_buy_discharge_push_indices( slots: list[PlanningSlot], pre_neg_buy_discharge_ts: set[int], @@ -1088,11 +1098,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_pressure: 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_pressure: vypne měkké v19 shortfall (pre-neg export, neg-buy soc/grid). """ T = len(slots) if T < 1: @@ -1253,12 +1265,9 @@ def solve_dispatch( st = slots[t] if float(st.sell_price) < _disch_sell_thr: continue - # Jen výrazný kladný sell; NE přidávat do discharge_export_slots (w_arb by - # zablokoval bd → load v noci → Infeasible). ge_bat povolíme přes pre_neg_buy_discharge_ts. + # Jen kladný sell (sell<0 + z_export by vynutilo ge_bat≥1 zároveň s ge_bat=0). if float(st.sell_price) >= 1.0: pre_neg_buy_discharge_ts.add(t) - elif st.allow_discharge_export: - pre_neg_buy_discharge_ts.add(t) # SELF_SUSTAIN dřív vynucoval ge[t] == 0, což umí udělat MILP infeasible v okamžiku, kdy: # - baterie je na max SoC (nelze nabíjet), # - PV pole B není curtailable, @@ -1515,34 +1524,35 @@ def solve_dispatch( float(battery.usable_capacity_wh), ) neg_sell_soc_underfill.append((t, us)) - for t in range(T): - if float(slots[t].buy_price) >= 0: - continue - us_buy = pulp.LpVariable( - f"neg_buy_soc_under_{t}", - 0, - float(battery.usable_capacity_wh), - ) - neg_buy_soc_underfill.append((t, us_buy)) - # Grid charge shortfall jen pokud je v slotu reálně headroom (soc panel min pod max). - headroom_wh = float(battery.soc_max_wh) - float(soc_panel_min[t]) - if headroom_wh < 500.0: - continue - cap_gi = float( - min( - battery.max_charge_power_w, - grid.max_import_power_w, - headroom_wh / max(INTERVAL_H * battery.charge_efficiency, 1e-6), + if not relaxed_neg_buy_pressure: + for t in _neg_buy_soc_pressure_slots(slots): + us_buy = pulp.LpVariable( + f"neg_buy_soc_under_{t}", + 0, + float(battery.usable_capacity_wh), ) - ) - if cap_gi < 500.0: - continue - sf_gi = pulp.LpVariable(f"neg_buy_gi_shortfall_{t}", 0, cap_gi) - neg_buy_grid_shortfall.append((t, sf_gi, cap_gi)) - for t in pre_neg_buy_discharge_ts: - cap_pre = float(_battery_export_cap_w(battery, grid)) - sf_pre = pulp.LpVariable(f"preneg_buy_disch_sf_{t}", 0, cap_pre) - pre_neg_buy_export_shortfall.append((t, sf_pre, cap_pre)) + neg_buy_soc_underfill.append((t, us_buy)) + for t in range(T): + if float(slots[t].buy_price) >= 0: + continue + headroom_wh = float(battery.soc_max_wh) - float(soc_panel_min[t]) + if headroom_wh < 500.0: + continue + cap_gi = float( + min( + battery.max_charge_power_w, + grid.max_import_power_w, + headroom_wh / max(INTERVAL_H * battery.charge_efficiency, 1e-6), + ) + ) + if cap_gi < 500.0: + continue + sf_gi = pulp.LpVariable(f"neg_buy_gi_shortfall_{t}", 0, cap_gi) + neg_buy_grid_shortfall.append((t, sf_gi, cap_gi)) + for t in pre_neg_buy_discharge_ts: + cap_pre = float(_battery_export_cap_w(battery, grid)) + sf_pre = pulp.LpVariable(f"preneg_buy_disch_sf_{t}", 0, cap_pre) + pre_neg_buy_export_shortfall.append((t, sf_pre, cap_pre)) for t in neg_sell_bat_dump_slots: dump_target_w = _battery_export_cap_w(battery, grid) sf_dump = pulp.LpVariable(f"neg_bat_dump_shortfall_{t}", 0, dump_target_w) @@ -1843,7 +1853,7 @@ def solve_dispatch( block_neg_sell_export_t = bool( getattr(grid, "block_export_on_negative_sell", False) ) - if t not in neg_sell_bat_dump_slots: + if t not in neg_sell_bat_dump_slots and t not in pre_neg_buy_discharge_ts: prob += ge_bat[t] == 0 ev_cap_neg = sum( float(vehicles[e].max_charge_power_w) @@ -2275,6 +2285,28 @@ def solve_dispatch( charge_commitment_prev_w=charge_commitment_prev_w, planner_version=planner_version, relaxed_expensive_import=True, + relaxed_neg_buy_pressure=relaxed_neg_buy_pressure, + ) + if not relaxed_neg_buy_pressure: + logger.warning( + "solve_dispatch Infeasible, retry with relaxed_neg_buy_pressure " + "(skip v19 soft shortfalls)" + ) + 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_pressure=True, ) raise RuntimeError(f"Solver: {pulp.LpStatus[status]}") @@ -2506,6 +2538,7 @@ def solve_dispatch( }, "load_first_enabled": om == "AUTO", "relaxed_expensive_import": relaxed_expensive_import, + "relaxed_neg_buy_pressure": relaxed_neg_buy_pressure, "charge_acquisition_buy_czk_kwh": charge_acquisition_czk_kwh, "charge_acquisition_cutoff_at": ( slots[0].charge_acquisition_cutoff_at.isoformat() diff --git a/backend/tests/test_planning_dispatch_milp.py b/backend/tests/test_planning_dispatch_milp.py index 0b00582..a423e65 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-27-pre-neg-buy-strategy-v19b") + self.assertEqual(snap.get("planner_build_tag"), "2026-05-27-pre-neg-buy-strategy-v19c") 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-27-pre-neg-buy-strategy-v19b") + self.assertEqual(snap.get("planner_build_tag"), "2026-05-27-pre-neg-buy-strategy-v19c") 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-27-pre-neg-buy-strategy-v19b") + self.assertEqual(snap.get("planner_build_tag"), "2026-05-27-pre-neg-buy-strategy-v19c") 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 948c719..661e821 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. - **Před oknem `buy < 0` (tag v19b, `solve_dispatch`):** sloty se **`sell < 0` + PV** v `charge_slots` mají **`bc_pv = 0`** (kapacita pro import); noční výboj přes **`pre_neg_buy_discharge_ts`** (`sell ≥ 1` nebo SQL `allow_discharge_export`) + shortfall **`ge_bat`** — **bez** přidání do **`discharge_export_slots`** (jinak Infeasible: `bd` k loadu). **`neg_sell_soc_underfill`** jen od **`first_neg_buy_idx`**; v **`buy < 0`** **`neg_buy_soc_underfill`** + **`neg_buy_grid_shortfall`**. Viz changelog v19/v19b. + **Před oknem `buy < 0` (tag v19c, `solve_dispatch`):** sloty se **`sell < 0` + PV** v `charge_slots` mají **`bc_pv = 0`** (kapacita pro import); noční výboj přes **`pre_neg_buy_discharge_ts`** (`sell ≥ 1` nebo SQL `allow_discharge_export`) + shortfall **`ge_bat`** — **bez** přidání do **`discharge_export_slots`** (jinak Infeasible: `bd` k loadu). **`neg_sell_soc_underfill`** jen od **`first_neg_buy_idx`**; v **`buy < 0`** **`neg_buy_soc_underfill`** + **`neg_buy_grid_shortfall`**. Viz changelog v19/v19b. 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/05-todo.md b/docs/05-todo.md index d59c065..b8d1768 100644 --- a/docs/05-todo.md +++ b/docs/05-todo.md @@ -29,7 +29,7 @@ Shrnutí otevřených bodů z `docs/06-open-questions.md`, checklistů v modulec | ~~**`charge_acquisition` po solve (two-pass):**~~ hotovo — `solve_dispatch_two_pass` v `planning_engine.py` (AUTO daily/rolling). | `planning_engine.py`, [`planning-arbitrage-accounting.md`](04-modules/planning-arbitrage-accounting.md) §6 | — | | ~~**Grid maska B (nejlevnější sloty):**~~ hotovo — `buy ASC` v AM/PM do Wh rozpočtu; cap z `ceil(budget/per_slot_wh)`. | `R__063` | — | | **Self-konzistentní filtr B + acquisition bez `buy<0`:** iterativní filtr v `R__063` (v12); vážená acquisition pro filtr i `charge_acquisition_buy_czk_kwh` jen z `allow_grid_charge` s `buy>=0` (záporný OTE buy zůstává `allow_charge`, ale neřítí exportní marži). Two-pass `_recompute_charge_acquisition_from_results` také přeskočí `buy<0`. Ověřit po deploy: `two_pass_converged=true` na home-01. | `R__063`, `planning_engine.py` | programátor | -| **Strategie buy<0 / sell<0 (home-01):** v19b — oprava Infeasible (pre-neg mimo discharge_export_slots). Deploy + replan, tag `…-v19b`. | `planning_engine.py`, `docs/planning-changelog.md` v19b | programátor | +| **Strategie buy<0 / sell<0 (home-01):** v19c — Infeasible fix (soc jen poslední buy<0 slot). Deploy + replan, tag `…-v19c`. | `planning_engine.py`, `docs/planning-changelog.md` v19c | programátor | | **KV1 replan timeout (~120 s):** ruční/rolling replan občas spadne na timeout; 5. pokus prošel. Profilovat `fn_load_planning_slots_full` (iterativní filtr) + MILP délku horizontu; případně zkrátit horizont pro test nebo zvýšit limit API. | backend replan endpoint, APScheduler | programátor | | **home-01 export při `sell<0` (26 slotů):** záměrně **ne** `block_export_on_negative_sell` (neriditelné PV B + zelený bonus). Plán stále může dávat `PV_SURPLUS` ~6–7 kW od ~10:30 když je SoC ~97 %+ — jiná osa než noční grid 4,8 Kč. Review ventilu `w_pv_b_vent_neg` / nabíjení před exportem, ne stejný fix jako KV1. | `planning_engine.py`, `planning-arbitrage-accounting.md` | programátor | diff --git a/docs/planning-changelog.md b/docs/planning-changelog.md index bdca038..7136c77 100644 --- a/docs/planning-changelog.md +++ b/docs/planning-changelog.md @@ -5,6 +5,17 @@ Formát: **datum (ISO)** · stručný důvod · soubory · chování / ověřen --- +## 2026-05-27 (k) — Infeasible: soc na každém buy<0 slotu + sell<0 v pre-neg (v19c) + +**Problém:** (1) **`neg_buy_soc_underfill`** na **každém** `buy<0` slotu vyžadoval `soc = soc_max` každých 15 min — při startu pod max fyzicky nemožné. (2) **`pre_neg_buy_discharge_ts`** mohlo zahrnout `sell<0` + `allow_discharge_export` → **`ge_bat=0`** (sell<0) vs **`z_export` → ge_bat≥1** → Infeasible. + +**Oprava (tag `2026-05-27-pre-neg-buy-strategy-v19c`):** +- `neg_buy_soc_underfill` jen na **posledním** `buy<0` slotu horizontu. +- `pre_neg_buy_discharge_ts` jen při **`sell ≥ 1`** (ne SQL discharge maska se záporným sell). +- Třetí retry: **`relaxed_neg_buy_pressure`** (vypne měkké shortfall, ponechá `bc_pv=0` před buy<0). + +--- + ## 2026-05-27 (j) — Infeasible: pre-neg export mimo discharge_export_slots (v19b) **Problém:** v19 přidávalo noční sloty (`sell ≥ 1`) do **`discharge_export_slots`** → režim **`w_arb`**: `bd` jen při exportu, ne k loadu → v noci nešlo pokrýt baseload → **Solver: Infeasible** (i po `relaxed_expensive_import`).