fix solveru
This commit is contained in:
@@ -59,7 +59,7 @@ NEG_SELL_CURTAIL_PENALTY_CZK_KWH = 1.0
|
|||||||
NEG_SELL_PV_CHARGE_REWARD_CZK_KWH = 0.8
|
NEG_SELL_PV_CHARGE_REWARD_CZK_KWH = 0.8
|
||||||
# Měkký tlak: v okně sell<0 dobít na soc_max (ne zastavit na ~94 % kvůli curtail).
|
# Měkký tlak: v okně sell<0 dobít na soc_max (ne zastavit na ~94 % kvůli curtail).
|
||||||
NEG_SELL_SOC_UNDERFILL_PENALTY_CZK_PER_WH = 0.35
|
NEG_SELL_SOC_UNDERFILL_PENALTY_CZK_PER_WH = 0.35
|
||||||
PLANNER_BUILD_TAG = "2026-05-24-ba81-pv-b-export-v5"
|
PLANNER_BUILD_TAG = "2026-05-24-ba81-gen-cutoff-v6"
|
||||||
CORRECTION_WINDOW_H = 1 # hodina zpět pro výpočet korekčního faktoru
|
CORRECTION_WINDOW_H = 1 # hodina zpět pro výpočet korekčního faktoru
|
||||||
CORRECTION_MIN_CLAMP = 0.5 # spodní limit korekčního faktoru
|
CORRECTION_MIN_CLAMP = 0.5 # spodní limit korekčního faktoru
|
||||||
CORRECTION_MAX_CLAMP = 1.5 # horní limit korekčního faktoru
|
CORRECTION_MAX_CLAMP = 1.5 # horní limit korekčního faktoru
|
||||||
@@ -1607,11 +1607,12 @@ def solve_dispatch(
|
|||||||
)
|
)
|
||||||
if block_pv_export_neg_sell:
|
if block_pv_export_neg_sell:
|
||||||
prob += ge_pv[t] == 0
|
prob += ge_pv[t] == 0
|
||||||
# Tvrdý zákaz celého vývozu (GEN / fixní nákup bez pole B).
|
# Tvrdý zákaz vývozu jen při block_export_on_negative_sell (KV1).
|
||||||
|
# GEN cut-off (z_gen_cutoff) nesmí vynutit ge==0 — jinak nelze odvést pole B při plné baterii (BA81).
|
||||||
block_neg_sell_export = bool(
|
block_neg_sell_export = bool(
|
||||||
getattr(grid, "block_export_on_negative_sell", False)
|
getattr(grid, "block_export_on_negative_sell", False)
|
||||||
)
|
)
|
||||||
if z_gen_cutoff is not None or block_neg_sell_export:
|
if block_neg_sell_export:
|
||||||
prob += ge[t] == 0
|
prob += ge[t] == 0
|
||||||
prob += ge_pv[t] == 0
|
prob += ge_pv[t] == 0
|
||||||
prob += ge_bat[t] == 0
|
prob += ge_bat[t] == 0
|
||||||
@@ -1812,7 +1813,10 @@ def solve_dispatch(
|
|||||||
and not getattr(grid, "block_export_on_negative_sell", False)
|
and not getattr(grid, "block_export_on_negative_sell", False)
|
||||||
)
|
)
|
||||||
if fixed_pv_b_export_cap:
|
if fixed_pv_b_export_cap:
|
||||||
prob += ge_pv[t] <= max(0.0, float(s.pv_b_forecast_w))
|
if z_gen_cutoff is not None:
|
||||||
|
prob += ge_pv[t] <= float(s.pv_b_forecast_w) * (1 - z_gen_cutoff[t])
|
||||||
|
else:
|
||||||
|
prob += ge_pv[t] <= max(0.0, float(s.pv_b_forecast_w))
|
||||||
if (
|
if (
|
||||||
not allow_pre_neg_pv_export
|
not allow_pre_neg_pv_export
|
||||||
and not skip_pv_store_block
|
and not skip_pv_store_block
|
||||||
|
|||||||
@@ -1222,7 +1222,7 @@ class NegativeSellPvChargeTests(unittest.TestCase):
|
|||||||
50.0,
|
50.0,
|
||||||
operating_mode="AUTO",
|
operating_mode="AUTO",
|
||||||
)
|
)
|
||||||
self.assertEqual(snap.get("planner_build_tag"), "2026-05-24-ba81-pv-b-export-v5")
|
self.assertEqual(snap.get("planner_build_tag"), "2026-05-24-ba81-gen-cutoff-v6")
|
||||||
self.assertGreater(
|
self.assertGreater(
|
||||||
results[0].battery_setpoint_w,
|
results[0].battery_setpoint_w,
|
||||||
5_500,
|
5_500,
|
||||||
@@ -1366,7 +1366,71 @@ class NegativeSellPvChargeTests(unittest.TestCase):
|
|||||||
50.0,
|
50.0,
|
||||||
operating_mode="AUTO",
|
operating_mode="AUTO",
|
||||||
)
|
)
|
||||||
self.assertEqual(snap.get("planner_build_tag"), "2026-05-24-ba81-pv-b-export-v5")
|
self.assertEqual(snap.get("planner_build_tag"), "2026-05-24-ba81-gen-cutoff-v6")
|
||||||
|
self.assertEqual(len(results), len(slots))
|
||||||
|
|
||||||
|
def test_gen_cutoff_full_soc_neg_sell_with_pv_b_feasible(self) -> None:
|
||||||
|
"""BA81: 100 % SoC + sell<0 + GEN cut-off — dříve ge==0 → Infeasible."""
|
||||||
|
slots = [
|
||||||
|
PlanningSlot(
|
||||||
|
interval_start=datetime(2026, 5, 24, 9, 0, tzinfo=timezone.utc)
|
||||||
|
+ timedelta(minutes=15 * i),
|
||||||
|
buy_price=3.088,
|
||||||
|
sell_price=-1.5,
|
||||||
|
pv_a_forecast_w=8_000,
|
||||||
|
pv_b_forecast_w=2_800,
|
||||||
|
load_baseline_w=500,
|
||||||
|
ev1_connected=False,
|
||||||
|
ev2_connected=False,
|
||||||
|
allow_charge=True,
|
||||||
|
allow_discharge_export=False,
|
||||||
|
charge_acquisition_buy_czk_kwh=3.61,
|
||||||
|
is_daytime_pv_surplus_slot=True,
|
||||||
|
safety_soc_target_wh=6_250.0,
|
||||||
|
)
|
||||||
|
for i in range(8)
|
||||||
|
]
|
||||||
|
battery = _battery(uc_wh=12_500.0, terminal_soc_value_factor=0.0)
|
||||||
|
battery.soc_max_wh = 12_500.0
|
||||||
|
battery.max_charge_power_w = 6_250
|
||||||
|
battery.max_discharge_power_w = 6_250
|
||||||
|
battery.degradation_cost_czk_kwh = 0.3
|
||||||
|
battery.planner_daytime_charge_target_enabled = True
|
||||||
|
hp = SimpleNamespace(
|
||||||
|
rated_heating_power_w=0,
|
||||||
|
tuv_min_temp_c=0.0,
|
||||||
|
tuv_target_temp_c=55.0,
|
||||||
|
)
|
||||||
|
grid = SimpleNamespace(
|
||||||
|
max_import_power_w=17_000,
|
||||||
|
max_export_power_w=16_000,
|
||||||
|
block_export_on_negative_sell=False,
|
||||||
|
deye_gen_microinverter_cutoff_enabled=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,
|
||||||
|
),
|
||||||
|
]
|
||||||
|
results, _ms, snap = solve_dispatch(
|
||||||
|
slots,
|
||||||
|
battery,
|
||||||
|
hp,
|
||||||
|
grid,
|
||||||
|
[None, None],
|
||||||
|
vehicles,
|
||||||
|
12_500.0,
|
||||||
|
55.0,
|
||||||
|
operating_mode="AUTO",
|
||||||
|
)
|
||||||
|
self.assertEqual(snap.get("planner_build_tag"), "2026-05-24-ba81-gen-cutoff-v6")
|
||||||
self.assertEqual(len(results), len(slots))
|
self.assertEqual(len(results), len(slots))
|
||||||
|
|
||||||
|
|
||||||
|
|||||||
@@ -44,7 +44,7 @@
|
|||||||
1. **Globální rozpočet Wh** (`discharge_slot_buffer × exportovatelná kapacita`): sloty podle `sell_price desc`. Před prvním `sell < 0` se z rozpočtu **vynechají** sloty, kde **později tentýž den** existuje `sell` vyšší o více než `degradation` (OTE, ne pevné hodiny 00–04).
|
1. **Globální rozpočet Wh** (`discharge_slot_buffer × exportovatelná kapacita`): sloty podle `sell_price desc`. Před prvním `sell < 0` se z rozpočtu **vynechají** sloty, kde **později tentýž den** existuje `sell` vyšší o více než `degradation` (OTE, ne pevné hodiny 00–04).
|
||||||
2. **Večerní špičky per den:** `sell ≥ max(sell) − degradation` jen pro hodiny **≥ 17** (Prague), ne globální max horizontu (jinak by vyhrála půlnoc 3,7 Kč místo večera).
|
2. **Večerní špičky per den:** `sell ≥ max(sell) − degradation` jen pro hodiny **≥ 17** (Prague), ne globální max horizontu (jinak by vyhrála půlnoc 3,7 Kč místo večera).
|
||||||
3. **Ranní pásmo před prvním `sell < 0`:** hodiny **5–11** téhož kalendářního dne — všechny sloty s `sell ≥ lokální_max_ráno − degradation`; ostatní sloty mezi ranním pásmem a prvním `sell < 0` s nižším sell mají export **zakázán** (žádný dump v 07:30 za 2 Kč). **`charge_acquisition`:** vážený `buy` před prvním exportem **téhož dne** jako záporné výkupní okno.
|
3. **Ranní pásmo před prvním `sell < 0`:** hodiny **5–11** téhož kalendářního dne — všechny sloty s `sell ≥ lokální_max_ráno − degradation`; ostatní sloty mezi ranním pásmem a prvním `sell < 0` s nižším sell mají export **zakázán** (žádný dump v 07:30 za 2 Kč). **`charge_acquisition`:** vážený `buy` před prvním exportem **téhož dne** jako záporné výkupní okno.
|
||||||
V `solve_dispatch` (AUTO): **`charge_slots`** = `allow_charge` z DB + **`buy < 0`** + všechny sloty **`sell < 0`** s PV přebytkem > 500 W (i bez `block_export_on_negative_sell`, BA81). **`pv_charge_shortfall`** / **`NEG_SELL_CURTAIL_PENALTY`** platí v těchto slotech. Při **`sell < 0`**: safety deficit cílí **`soc_max_wh`** (plný planner strop). Po posledním **`sell < 0`** tentýž den: **`post_neg_pv_topup`** dobije z FVE na `soc_max` před exportem (kladný sell, ne high-sell peak). U **fixního tarifu** s polem B: **`ge_pv ≤ pv_b`** místo tvrdého **`ge_pv = 0`** z pv_store (jinak Infeasible při plné baterii). **`planner_build_tag`** v `solver_params`. Changelog: [`docs/planning-changelog.md`](../planning-changelog.md).
|
V `solve_dispatch` (AUTO): **`charge_slots`** = `allow_charge` z DB + **`buy < 0`** + všechny sloty **`sell < 0`** s PV přebytkem > 500 W (i bez `block_export_on_negative_sell`, BA81). **`pv_charge_shortfall`** / **`NEG_SELL_CURTAIL_PENALTY`** platí v těchto slotech. Při **`sell < 0`**: safety deficit cílí **`soc_max_wh`** (plný planner strop). Po posledním **`sell < 0`** tentýž den: **`post_neg_pv_topup`** dobije z FVE na `soc_max` před exportem (kladný sell, ne high-sell peak). U **fixního tarifu** s polem B: **`ge_pv ≤ pv_b`** (ne pv_store **`ge_pv = 0`**). Při **`deye_gen_microinverter_cutoff_enabled`**: **`ge == 0` jen** pokud **`block_export_on_negative_sell`** (KV1), ne kvůli samotnému `z_gen_cutoff` (BA81 musí moci exportovat B při plné baterii). **`planner_build_tag`** v `solver_params`. Changelog: [`docs/planning-changelog.md`](../planning-changelog.md).
|
||||||
- **Záporná nákupní cena:**
|
- **Záporná nákupní cena:**
|
||||||
- horní mez `grid_import` zahrnuje `load_baseline_w` + nabíjení/EV/TČ (bez nekonečného importu).
|
- horní mez `grid_import` zahrnuje `load_baseline_w` + nabíjení/EV/TČ (bez nekonečného importu).
|
||||||
- **Záporná prodejní cena → tvrdý zákaz vývozu (`ge = 0`)** (`planning_engine.solve_dispatch`): platí ve slotu kde `sell_price < 0`, pokud lokality zapne některou z opcí —
|
- **Záporná prodejní cena → tvrdý zákaz vývozu (`ge = 0`)** (`planning_engine.solve_dispatch`): platí ve slotu kde `sell_price < 0`, pokud lokality zapne některou z opcí —
|
||||||
|
|||||||
@@ -5,15 +5,18 @@ Formát: **datum (ISO)** · stručný důvod · soubory · chování / ověřen
|
|||||||
|
|
||||||
---
|
---
|
||||||
|
|
||||||
## 2026-05-24 (j) — BA81: Solver Infeasible (plná baterie + pole B)
|
## 2026-05-24 (j) — BA81: Solver Infeasible (plná baterie + pole B + GEN cut-off)
|
||||||
|
|
||||||
**Problém:** Po deployi `evening-export-v4` u BA81 plánování padá na **`Solver: Infeasible`** (KV1 OK), typicky při vysokém SoC.
|
**Problém:** Po deployi večerních oprav u BA81 plánování padá na **`Solver: Infeasible`** (KV1 OK), typicky při **SoC ≈ 100 %** během dlouhého okna `sell < 0` (dnešní OTE).
|
||||||
|
|
||||||
**Příčina:** U fixního tarifu s **polem B** (`pv_b > 0`) guard `ge_pv = 0` (pv_store) blokoval export B při `sell < future_sell − degrad`, zatímco přebytek B nelze curtailovat — při plné baterii neřešitelné.
|
**Příčiny (dvě vrstvy):**
|
||||||
|
|
||||||
**Oprava:** tag **`2026-05-24-ba81-pv-b-export-v5`** v `planning_engine.py` — `ge_pv ≤ pv_b` místo `ge_pv = 0`; večerní tvrdý push jen top 3 sloty/den; bez rozšíření `profitable_export_ts` na všechny discharge masky.
|
1. **v5:** `ge_pv = 0` z pv_store při `pv_b > 0` → oprava `ge_pv ≤ pv_b`.
|
||||||
|
2. **v6 (skutečný blocker u BA81):** `deye_gen_microinverter_cutoff_enabled` společně s `sell < 0` vynucovalo **`ge == 0`** (podmínka `z_gen_cutoff is not None`). Při plné baterii nelze nabít ani exportovat přebytek pole B → Infeasible. BA81 má v kontextu `soc_wh = soc_max_wh = 12 500`.
|
||||||
|
|
||||||
**Ověření:** `pytest …::NegativeSellPvChargeTests`; MCP: `planner_build_tag = 2026-05-24-ba81-pv-b-export-v5`.
|
**Oprava:** tag **`2026-05-24-ba81-gen-cutoff-v6`** — `ge == 0` jen při `block_export_on_negative_sell`; `ge_pv ≤ pv_b × (1 − z_gen_cutoff)`; v5 večerní push + pv_b cap zůstávají.
|
||||||
|
|
||||||
|
**Ověření:** `pytest backend/tests/test_planning_dispatch_milp.py::NegativeSellPvChargeTests`; MCP po deployi: `planner_build_tag = 2026-05-24-ba81-gen-cutoff-v6`.
|
||||||
|
|
||||||
---
|
---
|
||||||
|
|
||||||
|
|||||||
Reference in New Issue
Block a user