From a03b45d4a958fade3d70de619eadde0814b0f5cb Mon Sep 17 00:00:00 2001 From: Dusan Vojacek Date: Sat, 30 May 2026 23:23:17 +0200 Subject: [PATCH] oprava zbytecneho curtailu A --- backend/services/control/setpoints.py | 8 ++ backend/services/planning_engine.py | 78 ++++++++++++++----- backend/tests/test_control_exporter_reg340.py | 19 +++++ backend/tests/test_planning_dispatch_milp.py | 60 +++++++++++++- docs/04-modules/modbus-registers.md | 2 +- docs/04-modules/planning.md | 3 +- docs/planning-changelog.md | 30 +++++++ 7 files changed, 177 insertions(+), 23 deletions(-) diff --git a/backend/services/control/setpoints.py b/backend/services/control/setpoints.py index d86c63c..e0b1ce1 100644 --- a/backend/services/control/setpoints.py +++ b/backend/services/control/setpoints.py @@ -137,7 +137,15 @@ def _build_setpoints( buy_raw = pi.get("effective_buy_price") buy_f: float | None = float(buy_raw) if buy_raw is not None else None pv_b = int(pi.get("pv_b_forecast_solver_w") or 0) + # Slabý úsvit: neposílat reg 340 — forecast nepřesný, Deye řídí sám (108/109/142). + _low_pv_no_reg340_w = 1500 if ( + forecast < _low_pv_no_reg340_w + and curtail <= 0 + and pv_b > 0 + ): + pv_a_allowed = None + elif ( buy_f is not None and sell_f is not None and float(buy_f) < 0.0 diff --git a/backend/services/planning_engine.py b/backend/services/planning_engine.py index 2d9aef5..d4c7748 100644 --- a/backend/services/planning_engine.py +++ b/backend/services/planning_engine.py @@ -71,7 +71,7 @@ NEG_BUY_CHARGE_SHORTFALL_PENALTY_CZK_KWH = 100.0 PRE_NEG_CHARGE_PENALTY_CZK_KWH = 400.0 PRE_NEG_BATT_EXPORT_SHORTFALL_PENALTY_CZK_KWH = 80.0 PRE_NEG_BATT_EXPORT_MIN_SELL_CZK_KWH = 1.0 -PLANNER_BUILD_TAG = "2026-05-31-evening-push-budget-primary-night-v49" +PLANNER_BUILD_TAG = "2026-05-31-ba81-dawn-no-micro-curtail-v51" # Ranní slabá FVE: neaplikovat pv_store ge_pv=0 (jinak curtail při sell < večerní peak). DAWN_LOW_PV_NO_CURTAIL_W = 1500 # Mimo evening_push: preferovat bd pro dům místo gi, když buy >> acq (účinná cena importu). @@ -1564,8 +1564,9 @@ def _slot_pv_surplus_w(slot: PlanningSlot) -> float: def _battery_export_push_defer_to_pv(slot: PlanningSlot) -> bool: """ - Při kladném sell a PV přebytku nevnucovat plný ge_bat push (pre-neg / ranní větve). - Exportní cap má pokrýt ge_pv; baterii řeší večerní push a sell<0 okna. + Při kladném sell a PV přebytku nevnucovat vývoz z baterie (ge_bat / z_export). + Platí pro ranní pre-neg, večerní push i KV1 odpoledne (block_export + fixní tarif): + přetok řeší ge_pv / Deye PASSIVE, ne BATTERY_SELL. """ if float(slot.sell_price) < 0.0: return False @@ -2629,6 +2630,11 @@ def solve_dispatch( slots, evening_push_ts ) night_self_consume_discourage_ts |= post_evening_push_night_ts + battery_export_defer_pv_ts = { + t for t in range(T) if _battery_export_push_defer_to_pv(slots[t]) + } + else: + battery_export_defer_pv_ts = set() pre_neg_buy_soc_ceiling_wh = _pre_neg_buy_soc_ceiling_wh( slots, first_neg_buy_idx=first_neg_buy_idx, @@ -3249,6 +3255,8 @@ def solve_dispatch( for t_peak in sorted(evening_push_ts): if t_peak not in discharge_export_slots: continue + if t_peak in battery_export_defer_pv_ts: + continue push_floor_w = _evening_push_battery_export_w( slots[t_peak], battery, grid ) @@ -3256,6 +3264,17 @@ def solve_dispatch( prob += z_export[t_peak] == 1 prob += ge_bat[t_peak] >= push_floor_w prob += soc[t_peak] >= float(discharge_floor_wh) + for t_pv in sorted(battery_export_defer_pv_ts): + if t_pv in evening_push_ts: + continue + if t_pv in morning_pre_neg_export_ts: + continue + if t_pv in pre_neg_buy_discharge_ts: + continue + if t_pv in pre_neg_buy_empty_ts: + continue + prob += ge_bat[t_pv] == 0 + prob += z_export[t_pv] == 0 # Ostatní profitable sloty: měkká shortfall penalizace (ne večerní push). if ( last_pos_sell_pre_neg_buy is not None @@ -3738,35 +3757,56 @@ def solve_dispatch( fixed_pre_neg_pv_export = ( purchase_fixed_pre and sell_t >= 0.0 - and pv_surplus_w > 500.0 + and pv_surplus_w > NIGHT_EXPORT_PV_SUNRISE_SURPLUS_W and ( first_neg_sell_idx is None or t < first_neg_sell_idx ) ) - skip_pv_store_block = ( - float(s.pv_b_forecast_w) > 0 + fixed_block_pv_surplus_export = ( + purchase_fixed_pre + and bool(getattr(grid, "block_export_on_negative_sell", False)) + and sell_t >= 0.0 + and pv_surplus_w > NIGHT_EXPORT_PV_SUNRISE_SURPLUS_W + ) + # BA81: ge_pv≤pv_b jen při významném poli A — při úsvitu nechat Deye bez plného curtail A. + fixed_mi_low_pv_surplus_export = ( + purchase_fixed_pre + and float(s.pv_b_forecast_w) > 0 and not getattr(grid, "block_export_on_negative_sell", False) - and sell_t < 0 - and buy_t >= 0.0 - and not purchase_fixed_pre - and ( - first_neg_buy_idx is None - or t < first_neg_buy_idx + and sell_t >= 0.0 + and int(s.pv_a_forecast_w) < DAWN_LOW_PV_NO_CURTAIL_W + and pv_surplus_w > 0.0 + ) + skip_pv_store_block = ( + ( + float(s.pv_b_forecast_w) > 0 + and not getattr(grid, "block_export_on_negative_sell", False) + and sell_t < 0 + and buy_t >= 0.0 + and not purchase_fixed_pre + and ( + first_neg_buy_idx is None + or t < first_neg_buy_idx + ) ) - ) or ( - # Spot: při sell>=0 neblokovat ge_pv (export vs bc_pv; večerní peak = ge_bat). - not purchase_fixed_pre - and sell_t >= 0 - and pv_surplus_w > 500 - ) or fixed_pre_neg_pv_export - # BA81: export pole B jen při kladném sell mimo pre-neg okno (jinak jen ge_pv≤pv_b → curtail A). + or ( + # Spot: při sell>=0 neblokovat ge_pv (export vs bc_pv; večerní peak = ge_bat). + not purchase_fixed_pre + and sell_t >= 0 + and pv_surplus_w > NIGHT_EXPORT_PV_SUNRISE_SURPLUS_W + ) + or fixed_pre_neg_pv_export + or fixed_block_pv_surplus_export + or fixed_mi_low_pv_surplus_export + ) fixed_pv_b_export_cap = ( purchase_fixed_pre and float(s.pv_b_forecast_w) > 0 and not getattr(grid, "block_export_on_negative_sell", False) and sell_t >= 0 and not fixed_pre_neg_pv_export + and int(s.pv_a_forecast_w) >= DAWN_LOW_PV_NO_CURTAIL_W ) if fixed_pre_neg_pv_export: prob += ge_pv[t] <= max(0.0, pv_surplus_w) diff --git a/backend/tests/test_control_exporter_reg340.py b/backend/tests/test_control_exporter_reg340.py index 57d7d07..b61936f 100644 --- a/backend/tests/test_control_exporter_reg340.py +++ b/backend/tests/test_control_exporter_reg340.py @@ -112,6 +112,25 @@ class BuildSetpointsReg340Tests(unittest.TestCase): assert sp is not None self.assertEqual(sp.pv_a_allowed_w, 0) + def test_skipped_low_pv_forecast_with_mi_no_curtail(self) -> None: + """BA81 úsvit: slabý forecast, bez curtail — EMS neposílá reg 340.""" + sp = _build_setpoints( + _auto_mode(), + _pi_base( + pv_a_forecast_solver_w=405, + pv_b_forecast_solver_w=49, + pv_a_curtailed_w=0, + grid_setpoint_w=-100, + battery_setpoint_w=0, + export_mode="PV_SURPLUS", + export_limit_w=100, + ), + pv_a_cap_w=32_000, + reg340_pv_a_control_enabled=True, + ) + assert sp is not None + self.assertIsNone(sp.pv_a_allowed_w) + def test_skipped_when_no_export_no_charge_no_curtail(self) -> None: sp = _build_setpoints( _auto_mode(), diff --git a/backend/tests/test_planning_dispatch_milp.py b/backend/tests/test_planning_dispatch_milp.py index 33c82b9..7bce326 100644 --- a/backend/tests/test_planning_dispatch_milp.py +++ b/backend/tests/test_planning_dispatch_milp.py @@ -3664,6 +3664,54 @@ class PreNegativeSellExportTests(unittest.TestCase): self.assertLess(r0.pv_a_curtailed_w, 500, "pole A nesmí jít do curtail při sell>0 před neg") self.assertLess(r0.grid_setpoint_w, -4000, "export přebytku A+B do site") + def test_ba81_dawn_low_pv_no_full_curtail_for_mi_cap(self) -> None: + """Úsvit: malé pole A + MI — ne plný curtail A kvůli ge_pv≤pv_b (reg 340).""" + prague = ZoneInfo("Europe/Prague") + slot = PlanningSlot( + interval_start=datetime(2026, 5, 31, 5, 15, tzinfo=prague).astimezone(timezone.utc), + buy_price=3.088, + sell_price=2.95, + pv_a_forecast_w=405, + pv_b_forecast_w=49, + load_baseline_w=400, + ev1_connected=False, + ev2_connected=False, + allow_charge=True, + allow_discharge_export=True, + charge_acquisition_buy_czk_kwh=3.088, + future_sell_opportunity_czk_kwh=6.5, + ) + 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=False, + purchase_pricing_mode="fixed", + ) + 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), + ] + res, _, _ = solve_dispatch( + [slot], + battery, + hp, + grid, + [None, None], + vehicles, + 0.12 * battery.soc_max_wh, + 50.0, + operating_mode="AUTO", + ) + r = res[0] + self.assertLess( + r.pv_a_curtailed_w, + int(slot.pv_a_forecast_w) - 50, + "nesmí useknout celé A kvůli fixed_pv_b_export_cap", + ) + def test_rolling_horizon_drains_to_reserve_before_first_neg(self) -> None: """Rolling bez D−1 večera: výboj před 1. sell<0 na reserve (+ slack).""" prague = ZoneInfo("Europe/Prague") @@ -3722,7 +3770,7 @@ class PreNegativeSellExportTests(unittest.TestCase): interval_start=(base + timedelta(minutes=15 * i)).astimezone(timezone.utc), buy_price=6.35, sell_price=sell, - pv_a_forecast_w=800 if sell < 4 else 200, + pv_a_forecast_w=4500 if i < 4 else (800 if sell < 4 else 200), pv_b_forecast_w=0, load_baseline_w=400, ev1_connected=False, @@ -3759,7 +3807,15 @@ class PreNegativeSellExportTests(unittest.TestCase): operating_mode="AUTO", ) self.assertLess(res[10].grid_setpoint_w, -500, "20:15 sell 0` (řiditelné pole A + nenulový strop střídače z `deye_reg340_max_solar_w` / `max_dc_input_w`). Bez bonusu nebo cap **0** EMS reg 340 **nezapisuje**. - **Hodnota:** z `ControlSetpoints.pv_a_allowed_w` (AUTO): bez curtailmentu = plný cap; při `pv_a_curtailed_w > 0` viz tabulka výše. Režimy **SELF_SUSTAIN / PRESERVE / CHARGE_CHEAP** mají `pv_a_allowed_w = None` → žádný zápis 340 z EMS v daném ticku. -- **Bez zápisu 340 (2026-05):** pokud plán má **bez exportu** (`export_mode = NONE` nebo `grid_setpoint_w ≥ 0` a `export_limit_w = 0`), **bez nabíjení baterie** (`battery_setpoint_w ≤ 0`) a **bez curtailu A** (`pv_a_curtailed_w = 0`), EMS reg 340 **neposílá** — Deye řídí PV A přes **108/109/142** a při **plné baterii** typicky **solar sell off** (hardware). Funkce `plan_skips_deye_reg340_write` v `setpoints.py`. **Plánovač v32:** škrcení A v okně `sell < 0` jde přes `pv_a_curtailed_w` → reg 340; registry 108/109 se kvůli fázím nemění. +- **Bez zápisu 340 (2026-05):** pokud plán má **bez exportu** (`export_mode = NONE` nebo `grid_setpoint_w ≥ 0` a `export_limit_w = 0`), **bez nabíjení baterie** (`battery_setpoint_w ≤ 0`) a **bez curtailu A** (`pv_a_curtailed_w = 0`), EMS reg 340 **neposílá** — Deye řídí PV A přes **108/109/142** a při **plné baterii** typicky **solar sell off** (hardware). Funkce `plan_skips_deye_reg340_write` v `setpoints.py`. **v51:** navíc při **`pv_a_forecast_solver_w < 1500`** a **`pv_a_curtailed_w = 0`** (úsvit + MI) → **bez reg 340** i při malém exportu. **Plánovač v32:** škrcení A v okně `sell < 0` jde přes `pv_a_curtailed_w` → reg 340; registry 108/109 se kvůli fázím nemění. - **Výjimka:** explicitní curtail v plánu nebo záporné buy+sell s PV B → `pv_a_allowed_w` se dopočítá / vynuluje jako dřív. - **Živé čtení:** `read_deye_registers_live` vrací **`reg340_max_solar_power_w`** (integer) jen pokud je přepínač zapnutý; jinak **`null`** (bez extra FC3 čtení reg 340). diff --git a/docs/04-modules/planning.md b/docs/04-modules/planning.md index 05b5f3b..950e032 100644 --- a/docs/04-modules/planning.md +++ b/docs/04-modules/planning.md @@ -18,8 +18,9 @@ - **Load-first (Deye, AUTO, tvrdý od v34):** 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 ≤ bc_gi + max(0, max_load − pv_forecast)`** (při vysoké FVE žádný fiktivní import = load); při `pv ≥ load + 500 W` **`pv_ld ≥ load`**; mimo `allow_discharge_export`: `bd ≤ load − pv_ld`, `pv_ld ≥ load − bd`. Tag `2026-05-28-load-first-hard-v34`. 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):** tvrdé `ge_pv = 0` jen pokud `sell < future_sell_opportunity − degradation` **a** `sell < 0` (spot), nebo u fixního tarifu dle `fixed_pv_b_export_cap`. Při **`sell ≥ 0` (spot home-01, KV1):** `ge_pv` **neblokuje** pv_store — solver volí export vs. `bc_pv` podle `−ge_pv×sell` a degradace; **baterii** na večerní peak drží `ge_bat` (`evening_early` / push), ne curtail FVE. **v31:** při `sell ≥ 0` + PV přebytek **není** plný `ge_bat` push z `pre_neg_buy_discharge` / ranních shortfallů (export cap pro FVE). **Před prvním `sell < 0`:** `allow_pre_neg_pv_export`. Tag `2026-05-28-morning-pv-export-priority-v31`. Testy `Home01PvStoreValueTests`, `PreNegativeSellExportTests`. + - **BA81 úsvit + MI (v51):** `fixed_pv_b_export_cap` (`ge_pv ≤ pv_b`) jen pokud **`pv_a_forecast ≥ 1500 W`** (`DAWN_LOW_PV_NO_CURTAIL_W`); při slabším A + přebytku → `fixed_mi_low_pv_surplus_export` (bez pv_store bloku). Exporter: při `forecast < 1500` a bez curtail A → **bez reg 340** (`setpoints.py`). Tag `2026-05-31-ba81-dawn-no-micro-curtail-v51`. Test `test_ba81_dawn_low_pv_no_full_curtail_for_mi_cap`. - **Drahý nákup → vlastní spotřeba z baterie:** mimo `allow_charge` platí `bd + pv_ld ≥ load_baseline + hp[t]` a `gi ≤ EV + hp[t]` (ne `hp_rated`). **Spot:** drahý slot = `buy > min(buy≥0) + degradace`. **Fixní nákup (DB `purchase_pricing_mode=fixed` nebo heuristika rozptylu buy < 0,25):** navíc `buy > charge_acquisition + degradace`. Na spotu **nesmí** `charge_acquisition` (~0,9 Kč) označit všechny sloty jako drahé → Infeasible (home-01). Při **Infeasible** solver jednou opakuje s `relaxed_expensive_import` (síť smí krmit baseload v drahých slotech; v `solver_params.inputs.relaxed_expensive_import=true`). Testy `AutoPassiveSelfConsumptionTests`, `test_spot_low_acquisition_does_not_mark_all_slots_expensive`, `test_negative_buy_in_horizon_does_not_block_all_grid_import`. - - **Záporný výkup (`sell < 0`) bez exportu:** `block_export_on_negative_sell` (KV1) **nebo** `purchase_pricing_mode=fixed` (BA81). **Spot (home-01):** `ge_pv=0` dokud není plná baterie; při plné jen ventil pole B (`ge_pv ≤ pv_b`, `w_pv_b_vent_neg`); výboj baterie při `sell<0` jen **12 slotů** před `buy ≤ planner_extreme_buy_threshold` (default −2), pokud spread do budoucna dává smysl — tag `2026-05-26-neg-sell-bat-dump-extreme-buy-v11`. Večerní discharge maska u spotu: denní peak ≥17:00 (ne `sell > ref_buy` v slotu). + - **Záporný výkup (`sell < 0`) bez exportu:** `block_export_on_negative_sell` (KV1) **nebo** `purchase_pricing_mode=fixed` (BA81). **Spot (home-01):** `ge_pv=0` dokud není plná baterie; při plné jen ventil pole B (`ge_pv ≤ pv_b`, `w_pv_b_vent_neg`); výboj baterie při `sell<0` jen **12 slotů** před `buy ≤ planner_extreme_buy_threshold` (default −2), pokud spread do budoucna dává smysl — tag `2026-05-26-neg-sell-bat-dump-extreme-buy-v11`. Večerní discharge maska u spotu: denní peak ≥17:00 (ne `sell > ref_buy` v slotu). **v50:** u **KV1** při `sell≥0` a PV přebytku >500 W i **po** 1. `sell<0` → `ge_pv` (PV_SURPLUS), ne tvrdý `ge_bat` z večerního peak/push. - **Pole B při sell<0 (home-01):** pokud `block_export_on_negative_sell = false`, LP nesmí vynutit `ge_pv = 0` (přebytek neriťitelného PV B). KV1 s `block_export = true` jen curtail A / nabíjení. - **`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. diff --git a/docs/planning-changelog.md b/docs/planning-changelog.md index b58b71c..6adfb58 100644 --- a/docs/planning-changelog.md +++ b/docs/planning-changelog.md @@ -5,6 +5,36 @@ Formát: **datum (ISO)** · stručný důvod · soubory · chování / ověřen --- +## 2026-05-31 — BA81 úsvit: žádný plný curtail A / zápis reg 340 (v51) + +**Problém:** Při malém ranním PV (např. **405 W** A, **49 W** B) LP kvůli `fixed_pv_b_export_cap` (`ge_pv ≤ pv_b`) **usekl celé pole A** (`curt_a = pv_a`) a exporter posílal **reg 340** z nepřesného forecastu — zbytečný HW zápis, baterie prázdná. + +**Změna (v51):** +- `fixed_pv_b_export_cap` jen když **`pv_a_forecast ≥ 1500 W`** (`DAWN_LOW_PV_NO_CURTAIL_W`). +- **`fixed_mi_low_pv_surplus_export`:** úsvit + MI + přebytek → neblokovat `ge_pv` přes pv_store. +- **`setpoints.py`:** při `forecast < 1500 W` a `curt_a = 0` → **`pv_a_allowed_w = None`** (bez reg 340). + +Tag **`2026-05-31-ba81-dawn-no-micro-curtail-v51`**. + +**Ověření:** `pytest … -k ba81_dawn_low_pv`; MCP BA81 05:15: `curt_a ≪ pv_a`. + +--- + +## 2026-05-31 — KV1/BA81: při PV přebytku FVE→síť, ne bat→síť (v50) + +**Problém:** **KV1** (`block_export`, fixní buy) odpoledne s FVE (~4 kW, sell>0) plán **BATTERY_SELL** místo **PV_SURPLUS** (home-01 OK). Příčiny: `skip_pv_store_block` jen před 1. `sell<0`; večerní **push** bez `defer_to_pv`; **z_export/ge_bat** u profitable peak. + +**Změna (v50):** +- **`fixed_block_pv_surplus_export`:** KV1 + `sell≥0` + PV přebytek → neblokovat `ge_pv` (pv_store). +- **`battery_export_defer_pv_ts`:** `ge_bat=0`, `z_export=0` (výjimky: morning pre-neg / pre-neg buy větve). +- **`evening_push_ts`:** přeskočit push, když platí defer. + +Tag **`2026-05-31-kv1-pv-surplus-over-bat-export-v50`**. + +**Ověření:** `pytest … -k kv1_evening_battery_push`; MCP KV1 17:00–18:00: `export_mode=PV_SURPLUS`, `curt_a` malé. + +--- + ## 2026-05-31 — Večerní push: celý Wh rozpočet jen pro dnešní noc (v49) **Problém (v43):** `push_budget / počet_kalendářních_večerů` dělil **aktuální SoC** mezi dnešní a **zítřejší** večer v horizontu — přes den FVE / neg nabíjení. Dnes večer dostal ~polovinu rozpočtu → chyběly sloty (např. 23:15); zítra večer push z dnešní SoC nedává smysl.