fix max grid kw

backend/services/planning_engine.py

@@ -318,6 +318,11 @@ def solve_dispatch(
 
     prob = pulp.LpProblem("ems_dispatch", pulp.LpMinimize)
 
+    # Penalizace překročení breakeru (Kč/kWh importu nad max_import_power_w).
+    # Záměr: breaker je fyzický strop, ale kvůli chybám forecastu a krátkým „extrémním“ oknům
+    # (např. záporná nákupní cena) umožníme solveru nominálně jít nad breaker, ovšem pouze za cenu.
+    IMPORT_OVER_BREAKER_PENALTY_CZK_KWH = 10.0
+
     min_soc_wh = float(getattr(battery, "min_soc_wh", battery.reserve_soc_wh))
     arb_base_wh = max(
         float(getattr(battery, "arb_floor_wh", battery.reserve_soc_wh)),
@@ -331,14 +336,15 @@ def solve_dispatch(
         )
 
     # --- Proměnné ---
-    # gi[t] horní mez: site breaker (max_import_power_w) je fyzický strop, ale o jeho dodržení
-    # se v reálném čase stará **Deye reg 128** (grid charge current) + firmware throttling —
-    # dynamicky sníží nabíjení baterie, když aktuální `load + bc` přesáhne breaker. Proto LP
-    # povolí nominálně import až **breaker + BMS max charge**, aby mohl plánovat `bc = BMS max`
-    # i v slotech s vyšší baseline zátěží (jinak tvrdý strop zbytečně osekává arbitráž v cenově
-    # nejlepších 15min oknech). Reálný hardware nikdy víc než breaker nenatáhne.
+    # gi[t] horní mez: site breaker (max_import_power_w) je fyzický strop.
+    # Pro robustnost (forecast PV/load nemusí sedět) používáme měkký cap: dovolíme gi nominálně
+    # až ~breaker + BMS max charge, ale překročení breakeru je penalizované (viz gi_over).
     gi_upper = float(grid.max_import_power_w) + float(battery.max_charge_power_w)
     gi  = [pulp.LpVariable(f"gi_{t}",  0, gi_upper)                        for t in range(T)]
+    gi_over = [
+        pulp.LpVariable(f"gi_over_{t}", 0, max(0.0, gi_upper - float(grid.max_import_power_w)))
+        for t in range(T)
+    ]
     ge  = [pulp.LpVariable(f"ge_{t}",  0, grid.max_export_power_w)         for t in range(T)]
     bc  = [pulp.LpVariable(f"bc_{t}",  0, battery.max_charge_power_w)      for t in range(T)]
     bd  = [pulp.LpVariable(f"bd_{t}",  0, battery.max_discharge_power_w)   for t in range(T)]
@@ -381,6 +387,7 @@ def solve_dispatch(
         pulp.lpSum(
             gi[t]  * slots[t].buy_price  * INTERVAL_H / 1000
             - ge[t] * slots[t].sell_price * INTERVAL_H / 1000
+            + gi_over[t] * IMPORT_OVER_BREAKER_PENALTY_CZK_KWH * INTERVAL_H / 1000
             + 0.5 * (bc[t] + bd[t]) * degradation_cost_effective * INTERVAL_H / 1000
             + pulp.lpSum(
                 ev_direct[e][t]  * slots[t].buy_price * INTERVAL_H / 1000
@@ -412,6 +419,9 @@ def solve_dispatch(
             == s.load_baseline_w + ev_total_t + hp[t] + bc[t] + ge[t]
         )
 
+        # Měkký breaker cap: gi_over[t] >= max(0, gi[t] - breaker).
+        prob += gi_over[t] >= gi[t] - float(grid.max_import_power_w)
+
         # SoC kontinuita
         soc_prev = current_soc_wh if t == 0 else soc[t - 1]
         prob += soc[t] == (