diff --git a/backend/services/planning_engine.py b/backend/services/planning_engine.py index cb11111..a9e4dd5 100644 --- 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( - t, - current_soc_wh=current_soc_wh, - battery=battery, - soc_headroom_wh=soc_headroom_wh, - pv_surplus_w=pv_surplus_w, + 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Č). diff --git a/backend/tests/test_planning_dispatch_milp.py b/backend/tests/test_planning_dispatch_milp.py index f345105..9502a22 100644 --- 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.""" diff --git a/docs/04-modules/planning.md b/docs/04-modules/planning.md index b4b34dd..ab5fdbf 100644 --- 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.