oprava KV1 nabijeni rano misto prodeje

2026-05-22 15:36:56 +02:00
parent f960e08307
commit c5525c729f
3 changed files with 97 additions and 19 deletions
--- a/backend/services/planning_engine.py
+++ b/backend/services/planning_engine.py
@@ -619,21 +619,31 @@ def _slots_with_charge_acquisition(
    ]


-def _pv_store_value_czk_kwh(
-    slot: PlanningSlot,
-    charge_acquisition_czk_kwh: float,
-    min_spread: float,
-) -> float:
+def _pv_store_value_czk_kwh(slot: PlanningSlot, min_spread: float) -> float:
    """
-    Minimální efektivní sell [Kč/kWh], pod kterým je FVE→síť horší než uložení
-    (večerní peak / náklad zásoby z levného nákupu).
+    Minimální sell [Kč/kWh], pod kterým je FVE→síť horší než uložení na večerní peak.
+    Používá jen future_sell_opportunity (ne charge_acquisition — u fixního tarifu KV1
+    by jinak blokoval export i při kladném sell 2 Kč).
    """
    future = float(
        slot.future_sell_opportunity_czk_kwh
        if slot.future_sell_opportunity_czk_kwh is not None
        else slot.sell_price
    )
-    return max(future, float(charge_acquisition_czk_kwh)) - min_spread
+    return future - min_spread
+
+
+def _pre_negative_sell_export_window(
+    slots: list[PlanningSlot],
+) -> tuple[int | None, int | None]:
+    """Index prvního sell<0 a posledního slotu před ním (pro strategii „vyvézt dřív“)."""
+    first_neg = next(
+        (i for i, s in enumerate(slots) if float(s.sell_price) < 0),
+        None,
+    )
+    if first_neg is None or first_neg <= 0:
+        return first_neg, None
+    return first_neg, first_neg - 1


 def _pv_forced_vent_export_allowed(
@@ -917,7 +927,7 @@ 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 = next((i for i, s in enumerate(slots) if float(s.sell_price) < 0), None)
+    first_neg_sell_idx, pre_neg_export_last_t = _pre_negative_sell_export_window(slots)
    if first_neg_sell_idx is not None and first_neg_sell_idx > 0 and floor_pct is not None:
        t_anchor = first_neg_sell_idx - 1
        soc_anchor_slack = pulp.LpVariable("soc_anchor_slack_wh", 0, float(battery.usable_capacity_wh))
@@ -1311,16 +1321,25 @@ def solve_dispatch(
                0.0,
                float(s.pv_a_forecast_w) + float(s.pv_b_forecast_w) - load_t,
            )
-            # FVE export jen pokud sell ≥ hodnota uložení (večerní peak / acquisition − degradace).
-            pv_store_val = _pv_store_value_czk_kwh(
-                s, charge_acquisition_czk_kwh, min_spread
+            # FVE export: před prvním sell<0 smí jít přebytek do sítě (kladný sell), pak nabít
+            # v záporném okně z PV. Jinak držet energii na future_sell peak.
+            allow_pre_neg_pv_export = (
+                first_neg_sell_idx is not None
+                and pre_neg_export_last_t is not None
+                and t <= pre_neg_export_last_t
+                and sell_t >= 0
            )
-            if sell_t < pv_store_val and not _pv_forced_vent_export_allowed(
+            pv_store_val = _pv_store_value_czk_kwh(s, min_spread)
+            if (
+                not allow_pre_neg_pv_export
+                and sell_t < pv_store_val
+                and not _pv_forced_vent_export_allowed(
                    t,
                    current_soc_wh=current_soc_wh,
                    battery=battery,
                    soc_headroom_wh=soc_headroom_wh,
                    pv_surplus_w=pv_surplus_w,
+                )
            ):
                prob += ge_pv[t] == 0
            # Drahý nákup: dům + TČ z baterie (ne import ze sítě); síť jen EV (+ případně TČ).
--- a/backend/tests/test_planning_dispatch_milp.py
+++ b/backend/tests/test_planning_dispatch_milp.py
@@ -1818,6 +1818,65 @@ class LoadFirstDispatchTests(unittest.TestCase):
        self.assertEqual(r.export_mode, "PV_SURPLUS")


+class PreNegativeSellExportTests(unittest.TestCase):
+    """Před prvním sell<0: export přebytku (BA81/KV1 strategie), ne nabíjení + pozdní vývoz."""
+
+    def test_kv1_like_morning_exports_before_negative_sell_window(self) -> None:
+        base = datetime(2026, 5, 22, 6, 45, tzinfo=timezone.utc)
+        slots = [
+            PlanningSlot(
+                interval_start=base + timedelta(minutes=15 * i),
+                buy_price=6.35,
+                sell_price=2.2,
+                pv_a_forecast_w=5000,
+                pv_b_forecast_w=0,
+                load_baseline_w=400,
+                ev1_connected=False,
+                ev2_connected=False,
+                allow_charge=False,
+                allow_discharge_export=False,
+                charge_acquisition_buy_czk_kwh=6.35,
+                future_sell_opportunity_czk_kwh=5.5,
+            )
+            for i in range(8)
+        ] + [
+            PlanningSlot(
+                interval_start=base + timedelta(hours=2),
+                buy_price=6.35,
+                sell_price=-0.3,
+                pv_a_forecast_w=6000,
+                pv_b_forecast_w=0,
+                load_baseline_w=400,
+                ev1_connected=False,
+                ev2_connected=False,
+                allow_charge=False,
+                allow_discharge_export=False,
+                charge_acquisition_buy_czk_kwh=6.35,
+                future_sell_opportunity_czk_kwh=-0.3,
+            ),
+        ]
+        battery = _battery(uc_wh=12_500.0, min_pct=10.0, arb_pct=30.0)
+        battery.max_charge_power_w = 6250
+        hp = SimpleNamespace(rated_heating_power_w=0, tuv_min_temp_c=45.0, tuv_target_temp_c=55.0)
+        grid = SimpleNamespace(
+            max_import_power_w=17_000,
+            max_export_power_w=8000,
+            block_export_on_negative_sell=True,
+        )
+        vehicles = [
+            SimpleNamespace(max_charge_power_w=0, battery_capacity_kwh=1.0, default_target_soc_pct=80.0),
+            SimpleNamespace(max_charge_power_w=0, battery_capacity_kwh=1.0, default_target_soc_pct=80.0),
+        ]
+        soc0 = 0.85 * battery.soc_max_wh
+        results, _, _ = solve_dispatch(
+            slots, battery, hp, grid, [None, None], vehicles, soc0, 50.0, operating_mode="AUTO"
+        )
+        self.assertLess(results[0].grid_setpoint_w, -500, "ráno: přebytek FVE do sítě před sell<0")
+        neg = results[8]
+        self.assertGreater(neg.battery_setpoint_w, 500, "záporný sell: PV do baterie")
+        self.assertEqual(neg.export_mode, "NONE")
+
+
 class Home01PvStoreValueTests(unittest.TestCase):
    """FVE (zejména pole B) nesmí jít do sítě pod hodnotou uložení / večerní peak."""

--- a/docs/04-modules/planning.md
+++ b/docs/04-modules/planning.md
@@ -16,7 +16,7 @@
  - **LP (AUTO):** objective explicitně `−ge_pv×sell − ge_bat×sell + ge_bat×acquisition` v exportních slotech; **bez** cross-slot vynucení `ge_pv ≥ surplus`. Guard FVE: `ge_pv=0` jen pokud `sell < charge_acquisition − degrad` (ne `sell < buy` ve slotu). Viz [`planning-arbitrage-accounting.md`](planning-arbitrage-accounting.md).
  - **Load-first (Deye, AUTO):** proměnné `pv_ld` (PV → load+EV+TČ), `pv_sp` (přebytek), `bc_pv` / `bc_gi`. Plná bilance `pv_a + pv_b + gi + bd = load + ev + hp + bc + ge`; `bc_pv + ge_pv ≤ pv_sp`; `gi ≤ load + bc_gi`; mimo `allow_discharge_export`: `bd ≤ load − pv_ld` a **`pv_ld ≥ load − gi − bd`**. Snapshot: `load_first_enabled=true`. Test `LoadFirstDispatchTests`.
  - **Tvrdé výkonové limity site/baterie:** `gi ≤ site_grid_connection.max_import_power_w` (breaker); **`bc_pv + bc_gi ≤ asset_battery.max_charge_power_w`**; **`ge ≤ max_export_power_w`** (proměnná `ge`, platí `ge = ge_pv + ge_bat`); **`bd + ge_bat ≤ asset_battery.max_discharge_power_w`** (vybíjení do domu + export z baterie nesmí současně překročit BMS). Dříve LP dovoloval import+nabíjení a dvojnásobné nabíjení; u prodeje hrozilo současné `bd` a `ge_bat` až 2× max discharge — viz `SitePowerCapTests`.
-  - **Hodnota FVE (PV store value):** `ge_pv = 0`, pokud `sell < max(future_sell_opportunity, charge_acquisition) − degradation` — přebytek jde do baterie (`bc_pv`), ne do sítě za haléře. Výjimka jen **nucený vent** (kotovací slot `t=0`, SoC u `soc_max`, přebytek PV). Dříve výjimka „jakékoli pole B s přebytkem“ obcházela guard → export i při sell 0,05 nebo −0,2 Kč/kWh (`Home01PvStoreValueTests`). Při **`sell < 0`** je `ge_pv = 0` (a `ge_bat = 0`) stejně, pokud není vent.
+  - **Hodnota FVE (PV store value):** `ge_pv = 0`, pokud `sell < future_sell_opportunity − degradation` (ne `charge_acquisition` — u fixního KV1 by jinak blokoval export při sell 2 Kč). **Před prvním `sell < 0` v horizontu:** při `sell ≥ 0` smí `ge_pv` až do `pv_sp` (strategie BA81: vyvézt přes poledne, pak nabít z FVE v záporném okně). Výjimka **nucený vent** jen plná baterie. Testy `Home01PvStoreValueTests`, `PreNegativeSellExportTests`.
  - **Drahý nákup → vlastní spotřeba z baterie:** mimo `allow_charge` platí `bd + pv_ld ≥ load_baseline + TČ` a `gi ≤ EV (+ TČ)` — ne import celého domu ze sítě při `buy > min(buy horizontu)` **nebo** `buy > charge_acquisition + degradace` (fixní tarif KV1, kde je každý slot „stejně drahý“). Testy `AutoPassiveSelfConsumptionTests`, `test_fixed_tariff_expensive_slot_discharges_not_grid_load`.
  - **`ref_buy_min` (brána exportu):** `min(buy_price)` horizontu — jen „existuje levný nákup?“, **ne** průměrná cena nabití přes hodiny. Export sloty: `sell > ref_buy_min + degradation` (spot). Viz [`planning-arbitrage-accounting.md`](planning-arbitrage-accounting.md).
  - Pokud `energy_to_fill <= 0` nebo `charge_slot_buffer = 0`: všechny sloty povoleny.