From c43bd0a6c619cdea86b32577aa4e347cf15b4ac6 Mon Sep 17 00:00:00 2001
From: Dusan Vojacek <dusan.vojacek@direct.cz>
Date: Sun, 24 May 2026 12:11:37 +0200
Subject: [PATCH] dalsi fix

---
 backend/services/planning_engine.py          | 30 ++++++++++++++++++--
 backend/tests/test_planning_dispatch_milp.py |  6 ++--
 docs/04-modules/planning.md                  |  2 +-
 docs/planning-changelog.md                   | 12 ++++++++
 4 files changed, 44 insertions(+), 6 deletions(-)

diff --git a/backend/services/planning_engine.py b/backend/services/planning_engine.py
index f8258d6..6e47a16 100644
--- a/backend/services/planning_engine.py
+++ b/backend/services/planning_engine.py
@@ -59,7 +59,7 @@ NEG_SELL_CURTAIL_PENALTY_CZK_KWH = 1.0
 NEG_SELL_PV_CHARGE_REWARD_CZK_KWH = 0.8
 # Měkký tlak: v okně sell<0 dobít na soc_max (ne zastavit na ~94 % kvůli curtail).
 NEG_SELL_SOC_UNDERFILL_PENALTY_CZK_PER_WH = 0.35
-PLANNER_BUILD_TAG = "2026-05-24-ba81-gen-cutoff-v6"
+PLANNER_BUILD_TAG = "2026-05-24-ba81-soc-headroom-v7"
 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
@@ -836,6 +836,28 @@ def _evening_battery_export_push_indices(
     return sorted(out)
 
 
+def _planner_soc_for_solver(
+    current_soc_wh: float,
+    battery,
+) -> tuple[float, float | None]:
+    """
+    SoC pro MILP. Při telemetrii na soc_max a dlouhém sell<0 s vysokou FVE bez rezervy pod stropem
+    je model neřešitelný (nelze nabít / odvést přebytek). Necháme min. ~650 Wh pod soc_max.
+    """
+    soc_max = float(battery.soc_max_wh)
+    soc_min = float(battery.min_soc_wh)
+    soc = max(soc_min, min(float(current_soc_wh), soc_max))
+    charge_slot_wh = (
+        float(battery.max_charge_power_w)
+        * INTERVAL_H
+        / max(float(battery.charge_efficiency), 1e-6)
+    )
+    headroom = max(650.0, 0.382 * charge_slot_wh)
+    if soc > soc_max - headroom:
+        return max(soc_min, soc_max - headroom), headroom
+    return soc, None
+
+
 def _pv_forced_vent_export_allowed(
     t: int,
     *,
@@ -1002,7 +1024,10 @@ def solve_dispatch(
             else float(sm)
             for i, sm in enumerate(soc_min_series)
         ]
-    current_soc_wh = float(current_soc_wh)
+    soc_headroom_applied_wh: float | None = None
+    current_soc_wh, soc_headroom_applied_wh = _planner_soc_for_solver(
+        current_soc_wh, battery
+    )
     current_soc_wh = max(soc_min_series[0], min(current_soc_wh, float(battery.soc_max_wh)))
     arb_base_wh = max(
         float(getattr(battery, "arb_floor_wh", battery.reserve_soc_wh)),
@@ -2078,6 +2103,7 @@ def solve_dispatch(
         "planner_build_tag": PLANNER_BUILD_TAG,
         "inputs": {
             "current_soc_wh": float(current_soc_wh),
+            "soc_headroom_applied_wh": soc_headroom_applied_wh,
             "operating_mode": operating_mode,
             "planner_version": planner_version_resolved,
             "battery": {
diff --git a/backend/tests/test_planning_dispatch_milp.py b/backend/tests/test_planning_dispatch_milp.py
index 204e016..4853b58 100644
--- a/backend/tests/test_planning_dispatch_milp.py
+++ b/backend/tests/test_planning_dispatch_milp.py
@@ -1222,7 +1222,7 @@ class NegativeSellPvChargeTests(unittest.TestCase):
             50.0,
             operating_mode="AUTO",
         )
-        self.assertEqual(snap.get("planner_build_tag"), "2026-05-24-ba81-gen-cutoff-v6")
+        self.assertEqual(snap.get("planner_build_tag"), "2026-05-24-ba81-soc-headroom-v7")
         self.assertGreater(
             results[0].battery_setpoint_w,
             5_500,
@@ -1366,7 +1366,7 @@ class NegativeSellPvChargeTests(unittest.TestCase):
             50.0,
             operating_mode="AUTO",
         )
-        self.assertEqual(snap.get("planner_build_tag"), "2026-05-24-ba81-gen-cutoff-v6")
+        self.assertEqual(snap.get("planner_build_tag"), "2026-05-24-ba81-soc-headroom-v7")
         self.assertEqual(len(results), len(slots))
 
     def test_gen_cutoff_full_soc_neg_sell_with_pv_b_feasible(self) -> None:
@@ -1430,7 +1430,7 @@ class NegativeSellPvChargeTests(unittest.TestCase):
             55.0,
             operating_mode="AUTO",
         )
-        self.assertEqual(snap.get("planner_build_tag"), "2026-05-24-ba81-gen-cutoff-v6")
+        self.assertEqual(snap.get("planner_build_tag"), "2026-05-24-ba81-soc-headroom-v7")
         self.assertEqual(len(results), len(slots))
 
 
diff --git a/docs/04-modules/planning.md b/docs/04-modules/planning.md
index 67cb4c4..0539cee 100644
--- a/docs/04-modules/planning.md
+++ b/docs/04-modules/planning.md
@@ -44,7 +44,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.
-    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). **`planner_build_tag`** v `solver_params`. Changelog: [`docs/planning-changelog.md`](../planning-changelog.md).
+    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).
 - **Záporná prodejní cena → tvrdý zákaz vývozu (`ge = 0`)** (`planning_engine.solve_dispatch`): platí ve slotu kde `sell_price < 0`, pokud lokality zapne některou z opcí —
diff --git a/docs/planning-changelog.md b/docs/planning-changelog.md
index 4892eaf..d5a6afc 100644
--- a/docs/planning-changelog.md
+++ b/docs/planning-changelog.md
@@ -5,6 +5,18 @@ Formát: **datum (ISO)** · stručný důvod · soubory · chování / ověřen
 
 ---
 
+## 2026-05-24 (k) — BA81: Infeasible při SoC = 100 % (telemetrie = soc_max)
+
+**Problém:** Po v6 stále `Solver: Infeasible` při replanu, když `fn_planning_site_context` vrátí `soc_wh = soc_max_wh` (12 500).
+
+**Příčina:** Při dlouhém `sell < 0` a vysoké FVE MILP potřebuje alespoň ~**650 Wh** rezervy pod `soc_max` pro modelování PV→baterie / export B. Na přesně 100 % SoC je model neřešitelný (reprodukce na datech runu 16184).
+
+**Oprava:** tag **`2026-05-24-ba81-soc-headroom-v7`** — `_planner_soc_for_solver()` sníží vstupní SoC na `soc_max − max(650 Wh, 0,382×slot_nabíjení)`; v `solver_params.inputs.soc_headroom_applied_wh` je audit.
+
+**Ověření:** `pytest …::NegativeSellPvChargeTests`; replan BA81 s telemetrií 100 % → tag v7, bez Infeasible.
+
+---
+
 ## 2026-05-24 (j) — BA81: Solver Infeasible (plná baterie + pole B + GEN cut-off)
 
 **Problém:** Po deployi večerních oprav u BA81 plánování padá na **`Solver: Infeasible`** (KV1 OK), typicky při **SoC ≈ 100 %** během dlouhého okna `sell < 0` (dnešní OTE).