fix max limitu
This commit is contained in:
@@ -789,10 +789,8 @@ def solve_dispatch(
|
|||||||
)
|
)
|
||||||
|
|
||||||
# --- Proměnné ---
|
# --- Proměnné ---
|
||||||
# gi[t] horní mez: site breaker (max_import_power_w) je fyzický strop.
|
# Import ze sítě: tvrdý strop = site breaker (max_import_power_w).
|
||||||
# Pro robustnost (forecast PV/load nemusí sedět) používáme měkký cap: dovolíme gi nominálně
|
gi_upper = float(grid.max_import_power_w)
|
||||||
# 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 = [pulp.LpVariable(f"gi_{t}", 0, gi_upper) for t in range(T)]
|
||||||
gi_over = [
|
gi_over = [
|
||||||
pulp.LpVariable(f"gi_over_{t}", 0, max(0.0, gi_upper - float(grid.max_import_power_w)))
|
pulp.LpVariable(f"gi_over_{t}", 0, max(0.0, gi_upper - float(grid.max_import_power_w)))
|
||||||
@@ -1018,6 +1016,11 @@ def solve_dispatch(
|
|||||||
)
|
)
|
||||||
pv_total_ub = float(s.pv_a_forecast_w) + float(s.pv_b_forecast_w)
|
pv_total_ub = float(s.pv_a_forecast_w) + float(s.pv_b_forecast_w)
|
||||||
|
|
||||||
|
# Součet nabíjení z FVE + ze sítě nesmí překročit max_charge_power_w baterie.
|
||||||
|
prob += bc_pv[t] + bc_gi[t] <= battery.max_charge_power_w
|
||||||
|
# Breaker: import ze site je tvrdě omezen (gi_over jen numerická pojistka).
|
||||||
|
prob += gi[t] <= gi_upper
|
||||||
|
|
||||||
if om == "AUTO":
|
if om == "AUTO":
|
||||||
load_site_expr = float(s.load_baseline_w) + ev_total_t + hp[t]
|
load_site_expr = float(s.load_baseline_w) + ev_total_t + hp[t]
|
||||||
prob += pv_ld[t] + pv_sp[t] == pv_a_net + pv_b_effective
|
prob += pv_ld[t] + pv_sp[t] == pv_a_net + pv_b_effective
|
||||||
@@ -1088,13 +1091,14 @@ def solve_dispatch(
|
|||||||
if z_gen_cutoff is not None and not allow_gen_cutoff:
|
if z_gen_cutoff is not None and not allow_gen_cutoff:
|
||||||
prob += z_gen_cutoff[t] == 0
|
prob += z_gen_cutoff[t] == 0
|
||||||
|
|
||||||
# Záporná nákupní cena → cap import (baseline domu + akumulace + řízené zátěže)
|
# Záporná nákupní cena → import jen na load + nabíjení + EV + TČ (stále ≤ breaker).
|
||||||
if s.buy_price < 0:
|
if s.buy_price < 0:
|
||||||
prob += gi[t] <= (
|
prob += gi[t] <= min(
|
||||||
s.load_baseline_w
|
gi_upper,
|
||||||
|
float(s.load_baseline_w)
|
||||||
+ battery.max_charge_power_w
|
+ battery.max_charge_power_w
|
||||||
+ sum(v.max_charge_power_w for v in vehicles)
|
+ sum(v.max_charge_power_w for v in vehicles)
|
||||||
+ heat_pump.rated_heating_power_w
|
+ heat_pump.rated_heating_power_w,
|
||||||
)
|
)
|
||||||
|
|
||||||
# Záporný prodej (sell < 0): baterii v tomhle okně nevybíjíme (dump má proběhnout předtím).
|
# Záporný prodej (sell < 0): baterii v tomhle okně nevybíjíme (dump má proběhnout předtím).
|
||||||
|
|||||||
@@ -1782,5 +1782,62 @@ class LoadFirstDispatchTests(unittest.TestCase):
|
|||||||
self.assertEqual(r.export_mode, "PV_SURPLUS")
|
self.assertEqual(r.export_mode, "PV_SURPLUS")
|
||||||
|
|
||||||
|
|
||||||
|
class SitePowerCapTests(unittest.TestCase):
|
||||||
|
"""Tvrdé limity site import a součtu nabíjení baterie."""
|
||||||
|
|
||||||
|
def test_grid_charge_respects_import_and_battery_caps(self) -> None:
|
||||||
|
"""home-01 typ: CHARGE slot nesmí překročit 17 kW import ani 18 kW do baterie."""
|
||||||
|
base = datetime(2026, 5, 22, 8, 45, tzinfo=timezone.utc)
|
||||||
|
slots = [
|
||||||
|
PlanningSlot(
|
||||||
|
interval_start=base,
|
||||||
|
buy_price=0.7,
|
||||||
|
sell_price=2.5,
|
||||||
|
pv_a_forecast_w=5000,
|
||||||
|
pv_b_forecast_w=0,
|
||||||
|
load_baseline_w=1961,
|
||||||
|
ev1_connected=False,
|
||||||
|
ev2_connected=False,
|
||||||
|
is_predicted_price=False,
|
||||||
|
allow_charge=True,
|
||||||
|
allow_discharge_export=False,
|
||||||
|
)
|
||||||
|
]
|
||||||
|
battery = _battery(uc_wh=64_000.0)
|
||||||
|
battery.max_charge_power_w = 18_000
|
||||||
|
battery.max_discharge_power_w = 18_000
|
||||||
|
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=13_500)
|
||||||
|
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.5 * battery.usable_capacity_wh
|
||||||
|
results, _, _ = solve_dispatch(
|
||||||
|
slots,
|
||||||
|
battery,
|
||||||
|
hp,
|
||||||
|
grid,
|
||||||
|
[None, None],
|
||||||
|
vehicles,
|
||||||
|
soc0,
|
||||||
|
50.0,
|
||||||
|
operating_mode="AUTO",
|
||||||
|
)
|
||||||
|
r = results[0]
|
||||||
|
self.assertLessEqual(
|
||||||
|
r.grid_setpoint_w,
|
||||||
|
17_000,
|
||||||
|
msg="import ze site ≤ max_import_power_w",
|
||||||
|
)
|
||||||
|
self.assertGreaterEqual(r.grid_setpoint_w, 0)
|
||||||
|
self.assertLessEqual(
|
||||||
|
r.battery_setpoint_w,
|
||||||
|
18_000,
|
||||||
|
msg="nabíjení baterie ≤ max_charge_power_w",
|
||||||
|
)
|
||||||
|
self.assertGreater(r.battery_setpoint_w, 0)
|
||||||
|
|
||||||
|
|
||||||
if __name__ == "__main__":
|
if __name__ == "__main__":
|
||||||
unittest.main()
|
unittest.main()
|
||||||
|
|||||||
@@ -14,7 +14,8 @@
|
|||||||
- **Grid ze sítě (vrstva B, před FVE):** spot, výchozí **AM/PM 50/50** z `grid_target × charge_slot_buffer` (do `soc_max`); **nevyčerpaný AM Wh přejde do PM** (`R__063`). Výběr: **nejlevnější `buy`** v pásmu (den plánu → před exportním oknem → `buy ASC`). Cap slotů: `ceil(budget/per_slot_wh) × charge_slot_buffer`. **`charge_acquisition`:** vážený `buy` u `allow_grid_charge` před 1. exportem; two-pass v `planning_engine.py`.
|
- **Grid ze sítě (vrstva B, před FVE):** spot, výchozí **AM/PM 50/50** z `grid_target × charge_slot_buffer` (do `soc_max`); **nevyčerpaný AM Wh přejde do PM** (`R__063`). Výběr: **nejlevnější `buy`** v pásmu (den plánu → před exportním oknem → `buy ASC`). Cap slotů: `ceil(budget/per_slot_wh) × charge_slot_buffer`. **`charge_acquisition`:** vážený `buy` u `allow_grid_charge` před 1. exportem; two-pass v `planning_engine.py`.
|
||||||
- **PV vrstva A:** jen pokud `sell ≥ future_sell_opportunity − degradation` (držet FVE na večerní peak, ne „nabíjet z FVE“ při nízkém sell).
|
- **PV vrstva A:** jen pokud `sell ≥ future_sell_opportunity − degradation` (držet FVE na večerní peak, ne „nabíjet z FVE“ při nízkém sell).
|
||||||
- **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).
|
- **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`** (jinak LP raději exportuje celou FVE a dům kryje baterie kvůli výnosu z `ge_pv`). Snapshot: `load_first_enabled=true`. Implementace: `planning_engine.py`; test `LoadFirstDispatchTests`.
|
- **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, ne měkká penalizace nad 17 kW); **`bc_pv + bc_gi ≤ asset_battery.max_charge_power_w`** (ne dva kanály po 18 kW). Dříve LP dovoloval `gi` až import+nabíjení a `battery_setpoint` až 2× max charge — viz `SitePowerCapTests`.
|
||||||
- **`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).
|
- **`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.
|
- Pokud `energy_to_fill <= 0` nebo `charge_slot_buffer = 0`: všechny sloty povoleny.
|
||||||
- **LP ekonomické guardy** (`solve_dispatch`, AUTO): `ge_pv=0` pokud `sell < charge_acquisition − degradation` (výjimka: plná baterie, přebytek **pv_b**). Pokud `buy > min(buy)+degradation` mimo charge masku → `gi` jen na load+EV+TČ. Viz `planning_engine.py` po slot pre-selection.
|
- **LP ekonomické guardy** (`solve_dispatch`, AUTO): `ge_pv=0` pokud `sell < charge_acquisition − degradation` (výjimka: plná baterie, přebytek **pv_b**). Pokud `buy > min(buy)+degradation` mimo charge masku → `gi` jen na load+EV+TČ. Viz `planning_engine.py` po slot pre-selection.
|
||||||
|
|||||||
Reference in New Issue
Block a user