stopadesaty fix

2026-05-25 01:52:08 +02:00
parent 5fb4c10ff6
commit a2a35981a1
5 changed files with 24 additions and 14 deletions
--- 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-v19"
+PLANNER_BUILD_TAG = "2026-05-27-pre-neg-buy-strategy-v19b"
 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
@@ -1253,10 +1253,11 @@ def solve_dispatch(
                st = slots[t]
                if float(st.sell_price) < _disch_sell_thr:
                    continue
-                # Jen sloty, kde SQL už povolilo export, nebo je výrazný kladný sell
-                # (rozšíření discharge_export_slots na celé dopoledne dělá krátké horizonty infeasible).
-                if st.allow_discharge_export or float(st.sell_price) >= 1.0:
-                    discharge_export_slots.add(t)
+                # 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.
+                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),
@@ -1539,8 +1540,6 @@ def solve_dispatch(
            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:
-            if t not in discharge_export_slots:
-                continue
            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))
@@ -1677,7 +1676,6 @@ def solve_dispatch(
        + pulp.lpSum(
            -35.0 * z_export[t]
            for t in pre_neg_buy_discharge_ts
-            if t in discharge_export_slots
        )
    )

@@ -2090,7 +2088,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

--- 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-v19")
+        self.assertEqual(snap.get("planner_build_tag"), "2026-05-27-pre-neg-buy-strategy-v19b")
        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-v19")
+        self.assertEqual(snap.get("planner_build_tag"), "2026-05-27-pre-neg-buy-strategy-v19b")
        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-v19")
+        self.assertEqual(snap.get("planner_build_tag"), "2026-05-27-pre-neg-buy-strategy-v19b")
        self.assertEqual(len(results), len(slots))

    def test_fixed_tariff_neg_sell_no_grid_export(self) -> None:
--- 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 v19, `solve_dispatch`):** sloty se **`sell < 0` + PV** v `charge_slots` mají **`bc_pv = 0`** (kapacita pro import); noční/dopolední výboj přes **`pre_neg_buy_discharge_ts`** (SQL `allow_discharge_export` nebo **`sell ≥ 1` Kč/kWh**) + shortfall **`ge_bat`**; **`neg_sell_soc_underfill`** jen od **`first_neg_buy_idx`**; v **`buy < 0`** slotech **`neg_buy_soc_underfill`** + **`neg_buy_grid_shortfall`** (headroom z `soc_panel_min`). Viz [`docs/planning-changelog.md`](../planning-changelog.md) v19.
+    **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.
    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).
--- 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&lt;0 / sell&lt;0 (home-01):** v19 — noční výboj, bc_pv=0 před buy&lt;0, import v buy&lt;0. Ověřit po deploy + replan. Refactor: `planning_engine.py` ~3100 řádků → vyčlenit constraints modul. | `planning_engine.py`, `docs/planning-changelog.md` v19 | programátor |
+| **Strategie buy&lt;0 / sell&lt;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 |
 | **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 |

--- 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-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`).
+
+**Oprava (tag `2026-05-27-pre-neg-buy-strategy-v19b`):** `pre_neg_buy_discharge_ts` pouze povolí **`ge_bat`** (+ shortfall), **bez** rozšíření `discharge_export_slots`. Baterie dál může vybíjet k domu (`bd`) a paralelně exportovat (`ge_bat`).
+
+---
+
 ## 2026-05-27 (i) — strategie před buy&lt;0: noční výboj, bez PV→bat, import v záporném nákupu (v19)

 **Problém (home-01 run 16662, tag v18):** Večerní/ranní export OK. Zbývá: (1) noc jen ~500 W do domu, žádný `ge_bat` výboj před `buy<0`; (2) 08:45–11:30 nabíjení z PV A do ~98 % ještě před `buy<0` (13:00); (3) v `buy<0` baterie plná → žádný import; (4) `neg_sell_soc_underfill` tlačilo na soc_max už v ranním `sell<0` okně.