a zase dalsi
This commit is contained in:
@@ -64,7 +64,9 @@ NEG_SELL_PV_B_VENT_PENALTY_CZK_KWH = 4.0
|
|||||||
# Výboj baterie při sell<0 jen těsně před extrémně záporným buy (round-trip arbitráž).
|
# Výboj baterie při sell<0 jen těsně před extrémně záporným buy (round-trip arbitráž).
|
||||||
EXTREME_BUY_DUMP_PREWINDOW_SLOTS = 12
|
EXTREME_BUY_DUMP_PREWINDOW_SLOTS = 12
|
||||||
NEG_SELL_BAT_DUMP_SHORTFALL_PENALTY_CZK_KWH = 80.0
|
NEG_SELL_BAT_DUMP_SHORTFALL_PENALTY_CZK_KWH = 80.0
|
||||||
PLANNER_BUILD_TAG = "2026-05-28-morning-hold-soc-v21b"
|
NEG_BUY_CHARGE_SHORTFALL_PENALTY_CZK_KWH = 100.0
|
||||||
|
PRE_NEG_CHARGE_PENALTY_CZK_KWH = 250.0
|
||||||
|
PLANNER_BUILD_TAG = "2026-05-28-buy-sell-split-v22"
|
||||||
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
|
||||||
@@ -1387,6 +1389,7 @@ def solve_dispatch(
|
|||||||
pv_charge_shortfall: list[tuple[int, pulp.LpVariable, float]] = []
|
pv_charge_shortfall: list[tuple[int, pulp.LpVariable, float]] = []
|
||||||
neg_sell_bat_dump_shortfall: list[tuple[int, pulp.LpVariable, float]] = []
|
neg_sell_bat_dump_shortfall: list[tuple[int, pulp.LpVariable, float]] = []
|
||||||
neg_sell_soc_underfill: list[tuple[int, pulp.LpVariable]] = []
|
neg_sell_soc_underfill: list[tuple[int, pulp.LpVariable]] = []
|
||||||
|
neg_buy_charge_shortfall: list[tuple[int, pulp.LpVariable, float]] = []
|
||||||
fixed_tariff_like = fixed_tariff_like_pre
|
fixed_tariff_like = fixed_tariff_like_pre
|
||||||
block_export_neg_sell = bool(getattr(grid, "block_export_on_negative_sell", False))
|
block_export_neg_sell = bool(getattr(grid, "block_export_on_negative_sell", False))
|
||||||
if om == "AUTO":
|
if om == "AUTO":
|
||||||
@@ -1406,9 +1409,17 @@ def solve_dispatch(
|
|||||||
))
|
))
|
||||||
sf = pulp.LpVariable(f"export_shortfall_{t}", 0, cap_w)
|
sf = pulp.LpVariable(f"export_shortfall_{t}", 0, cap_w)
|
||||||
peak_export_shortfall.append((t, sf, cap_w))
|
peak_export_shortfall.append((t, sf, cap_w))
|
||||||
|
for t in range(T):
|
||||||
|
if float(slots[t].buy_price) >= 0.0:
|
||||||
|
continue
|
||||||
|
cap_w = float(battery.max_charge_power_w)
|
||||||
|
sf_nb = pulp.LpVariable(f"neg_buy_charge_sf_{t}", 0, cap_w)
|
||||||
|
neg_buy_charge_shortfall.append((t, sf_nb, cap_w))
|
||||||
for t in range(T):
|
for t in range(T):
|
||||||
if float(slots[t].sell_price) >= 0:
|
if float(slots[t].sell_price) >= 0:
|
||||||
continue
|
continue
|
||||||
|
if float(slots[t].buy_price) < 0.0:
|
||||||
|
continue
|
||||||
if t not in charge_slots:
|
if t not in charge_slots:
|
||||||
continue
|
continue
|
||||||
# Před buy<0: nepenalizovat / netlačit PV→bat (jinak 98 % v 09:15 a export v sell<0).
|
# Před buy<0: nepenalizovat / netlačit PV→bat (jinak 98 % v 09:15 a export v sell<0).
|
||||||
@@ -1444,6 +1455,8 @@ def solve_dispatch(
|
|||||||
for t in range(T):
|
for t in range(T):
|
||||||
if float(slots[t].sell_price) >= 0:
|
if float(slots[t].sell_price) >= 0:
|
||||||
continue
|
continue
|
||||||
|
if float(slots[t].buy_price) < 0.0:
|
||||||
|
continue
|
||||||
if t not in charge_slots:
|
if t not in charge_slots:
|
||||||
continue
|
continue
|
||||||
if first_neg_buy_idx is not None and t < first_neg_buy_idx:
|
if first_neg_buy_idx is not None and t < first_neg_buy_idx:
|
||||||
@@ -1575,6 +1588,22 @@ def solve_dispatch(
|
|||||||
sf * NEG_SELL_BAT_DUMP_SHORTFALL_PENALTY_CZK_KWH * INTERVAL_H / 1000.0
|
sf * NEG_SELL_BAT_DUMP_SHORTFALL_PENALTY_CZK_KWH * INTERVAL_H / 1000.0
|
||||||
for _t, sf, _cap in neg_sell_bat_dump_shortfall
|
for _t, sf, _cap in neg_sell_bat_dump_shortfall
|
||||||
)
|
)
|
||||||
|
+ pulp.lpSum(
|
||||||
|
sf * NEG_BUY_CHARGE_SHORTFALL_PENALTY_CZK_KWH * INTERVAL_H / 1000.0
|
||||||
|
for _t, sf, _cap in neg_buy_charge_shortfall
|
||||||
|
)
|
||||||
|
+ pulp.lpSum(
|
||||||
|
(bc_pv[t] + bc_gi[t])
|
||||||
|
* PRE_NEG_CHARGE_PENALTY_CZK_KWH
|
||||||
|
* INTERVAL_H
|
||||||
|
/ 1000.0
|
||||||
|
for t in range(T)
|
||||||
|
if (
|
||||||
|
first_neg_buy_idx is not None
|
||||||
|
and t < first_neg_buy_idx
|
||||||
|
and float(slots[t].buy_price) >= 0.0
|
||||||
|
)
|
||||||
|
)
|
||||||
+ pulp.lpSum(
|
+ pulp.lpSum(
|
||||||
-25.0 * z_export[t]
|
-25.0 * z_export[t]
|
||||||
for t in range(T)
|
for t in range(T)
|
||||||
@@ -1591,6 +1620,8 @@ def solve_dispatch(
|
|||||||
prob += sf >= cap_w - ge_bat[t_sf]
|
prob += sf >= cap_w - ge_bat[t_sf]
|
||||||
for t_us, us in neg_sell_soc_underfill:
|
for t_us, us in neg_sell_soc_underfill:
|
||||||
prob += us >= float(battery.soc_max_wh) - soc[t_us]
|
prob += us >= float(battery.soc_max_wh) - soc[t_us]
|
||||||
|
for t_sf, sf, cap_w in neg_buy_charge_shortfall:
|
||||||
|
prob += sf >= cap_w - (bc_gi[t_sf] + bc_pv[t_sf])
|
||||||
preneg_export_min_soc_wh = float(min_soc_wh) + max(
|
preneg_export_min_soc_wh = float(min_soc_wh) + max(
|
||||||
float(battery.max_discharge_power_w)
|
float(battery.max_discharge_power_w)
|
||||||
* float(battery.discharge_efficiency)
|
* float(battery.discharge_efficiency)
|
||||||
@@ -1731,6 +1762,9 @@ def solve_dispatch(
|
|||||||
+ 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,
|
||||||
)
|
)
|
||||||
|
prob += ge[t] == 0
|
||||||
|
prob += ge_pv[t] == 0
|
||||||
|
prob += ge_bat[t] == 0
|
||||||
|
|
||||||
# Záporný prodej (sell < 0): výboj baterie jen před extrémně záporným buy (v11).
|
# Záporný prodej (sell < 0): výboj baterie jen před extrémně záporným buy (v11).
|
||||||
# Export FVE při sell<0: spot = nabíjení/curtail A; ventil jen pole B při plné baterii.
|
# Export FVE při sell<0: spot = nabíjení/curtail A; ventil jen pole B při plné baterii.
|
||||||
@@ -1952,7 +1986,14 @@ def solve_dispatch(
|
|||||||
before_neg_buy = (
|
before_neg_buy = (
|
||||||
first_neg_buy_idx is not None and t < first_neg_buy_idx
|
first_neg_buy_idx is not None and t < first_neg_buy_idx
|
||||||
)
|
)
|
||||||
if before_neg_buy and sell_t_pre < 0.0 and pv_surplus_w > 0:
|
if (
|
||||||
|
before_neg_buy
|
||||||
|
and float(s.buy_price) >= 0.0
|
||||||
|
and bool(getattr(s, "is_daytime_pv_surplus_slot", False))
|
||||||
|
):
|
||||||
|
prob += bc_pv[t] == 0
|
||||||
|
prob += bc_gi[t] == 0
|
||||||
|
elif before_neg_buy and sell_t_pre < 0.0 and pv_surplus_w > 0:
|
||||||
# Ranní sell<0 před buy<0: PV do sítě/curtail, ne do baterie (kapacita na import).
|
# Ranní sell<0 před buy<0: PV do sítě/curtail, ne do baterie (kapacita na import).
|
||||||
prob += bc_pv[t] == 0
|
prob += bc_pv[t] == 0
|
||||||
elif t not in charge_slots:
|
elif t not in charge_slots:
|
||||||
@@ -2005,6 +2046,7 @@ def solve_dispatch(
|
|||||||
float(s.pv_b_forecast_w) > 0
|
float(s.pv_b_forecast_w) > 0
|
||||||
and not getattr(grid, "block_export_on_negative_sell", False)
|
and not getattr(grid, "block_export_on_negative_sell", False)
|
||||||
and sell_t < 0
|
and sell_t < 0
|
||||||
|
and buy_t >= 0.0
|
||||||
and not purchase_fixed_pre
|
and not purchase_fixed_pre
|
||||||
) or (
|
) or (
|
||||||
# KV1: plná baterie + kladný sell — neblokovat ge_pv==0 (jinak masivní curtail).
|
# KV1: plná baterie + kladný sell — neblokovat ge_pv==0 (jinak masivní curtail).
|
||||||
@@ -2046,6 +2088,7 @@ def solve_dispatch(
|
|||||||
# zablokovalo legitimní pre-neg-pv export pole A z testů).
|
# zablokovalo legitimní pre-neg-pv export pole A z testů).
|
||||||
if (
|
if (
|
||||||
sell_t < 0
|
sell_t < 0
|
||||||
|
and buy_t >= 0.0
|
||||||
and float(s.pv_b_forecast_w) > 0
|
and float(s.pv_b_forecast_w) > 0
|
||||||
and not getattr(grid, "block_export_on_negative_sell", False)
|
and not getattr(grid, "block_export_on_negative_sell", False)
|
||||||
):
|
):
|
||||||
|
|||||||
@@ -1230,7 +1230,7 @@ class NegativeSellPvChargeTests(unittest.TestCase):
|
|||||||
50.0,
|
50.0,
|
||||||
operating_mode="AUTO",
|
operating_mode="AUTO",
|
||||||
)
|
)
|
||||||
self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-morning-hold-soc-v21b")
|
self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-buy-sell-split-v22")
|
||||||
self.assertGreater(
|
self.assertGreater(
|
||||||
results[0].battery_setpoint_w,
|
results[0].battery_setpoint_w,
|
||||||
5_500,
|
5_500,
|
||||||
@@ -1380,7 +1380,7 @@ class NegativeSellPvChargeTests(unittest.TestCase):
|
|||||||
50.0,
|
50.0,
|
||||||
operating_mode="AUTO",
|
operating_mode="AUTO",
|
||||||
)
|
)
|
||||||
self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-morning-hold-soc-v21b")
|
self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-buy-sell-split-v22")
|
||||||
self.assertEqual(len(results), len(slots))
|
self.assertEqual(len(results), len(slots))
|
||||||
|
|
||||||
def test_gen_cutoff_full_soc_neg_sell_with_pv_b_feasible(self) -> None:
|
def test_gen_cutoff_full_soc_neg_sell_with_pv_b_feasible(self) -> None:
|
||||||
@@ -1444,7 +1444,7 @@ class NegativeSellPvChargeTests(unittest.TestCase):
|
|||||||
55.0,
|
55.0,
|
||||||
operating_mode="AUTO",
|
operating_mode="AUTO",
|
||||||
)
|
)
|
||||||
self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-morning-hold-soc-v21b")
|
self.assertEqual(snap.get("planner_build_tag"), "2026-05-28-buy-sell-split-v22")
|
||||||
self.assertEqual(len(results), len(slots))
|
self.assertEqual(len(results), len(slots))
|
||||||
|
|
||||||
def test_fixed_tariff_neg_sell_no_grid_export(self) -> None:
|
def test_fixed_tariff_neg_sell_no_grid_export(self) -> None:
|
||||||
|
|||||||
@@ -45,7 +45,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.
|
||||||
**Planner tag v21b:** drží SoC před prvním `buy<0` (`bc_pv=0`, shortfally jen od `buy<0`); před `buy<0` v sell<0 povolen export **pole B** (`ge_pv≤pv_b` bez ventilu plné baterie — jinak Infeasible). Stále **v17** + **v18**. Viz changelog v21/v21b.
|
**Planner tag v22:** oddělení **záporného nákupu** (`buy<0` → import/nabíjení, `ge=0`) a **záporného výkupu** (`sell<0`, `buy≥0` → export pole B / curtail A). Před `buy<0` při `buy≥0` žádné nabíjení (`bc_pv=bc_gi=0`). Viz changelog v22.
|
||||||
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). Vstupní **`soc_wh`** z telemetrie se před MILP omezí přes **`_planner_soc_for_solver`** (rezerva ~650 Wh pod `soc_max`, jinak Infeasible při 100 % SoC a dlouhém záporném výkupu). **`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). Vstupní **`soc_wh`** z telemetrie se před MILP omezí přes **`_planner_soc_for_solver`** (rezerva ~650 Wh pod `soc_max`, jinak Infeasible při 100 % SoC a dlouhém záporném výkupu). **`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).
|
||||||
|
|||||||
@@ -5,6 +5,22 @@ Formát: **datum (ISO)** · stručný důvod · soubory · chování / ověřen
|
|||||||
|
|
||||||
---
|
---
|
||||||
|
|
||||||
|
## 2026-05-28 — rozlišení buy<0 vs sell<0 (v22)
|
||||||
|
|
||||||
|
**Problém (MCP run 16706, v21b):** Znaménka v objective OK (`grid<0` = export, `bat>0` = nabíjení). Chování ale „opačně“:
|
||||||
|
- **Před buy<0** (05:30–07:00, buy≥0): nabíjení z PV/sítě místo přípravy kapacity.
|
||||||
|
- **Při buy<0** (12:15–12:45): export do sítě místo importu — ventil `w_pv_b_vent` u sell<0 platil i když buy<0.
|
||||||
|
|
||||||
|
**Oprava (tag `2026-05-28-buy-sell-split-v22`):**
|
||||||
|
- **Před `first_neg_buy_idx` a buy≥0:** tvrdé `bc_pv=bc_gi=0` jen v `is_daytime_pv_surplus_slot` (SQL); jinak měkká penalizace `PRE_NEG_CHARGE_PENALTY`.
|
||||||
|
- **sell<0 a buy≥0:** export pole B / curtail A (v21b), bez `neg_sell_soc` shortfallu v buy<0 slotech.
|
||||||
|
- **buy<0:** tvrdě `ge=ge_pv=ge_bat=0` + měkký **`neg_buy_charge_shortfall`** (tlak na `bc_gi+bc_pv`).
|
||||||
|
- **sell<0 + buy<0:** žádný větev ventilu plné baterie → jen nabíjení/curtail.
|
||||||
|
|
||||||
|
**Ověření:** replan home-01 → tag v22; 11:00–11:45 import+nabíjení, 12:15 bez exportu při buy<0.
|
||||||
|
|
||||||
|
---
|
||||||
|
|
||||||
## 2026-05-28 — ranní sell<0: držet SoC před buy<0 (v21 / v21b)
|
## 2026-05-28 — ranní sell<0: držet SoC před buy<0 (v21 / v21b)
|
||||||
|
|
||||||
**Problém (MCP run 16692, tag v20):** Od ~05:30 nabíjení z PV; v 09:15 už **98,3 %** SoC; od 09:15 masivní **export při sell<0** (−7 kW). V **11:00–12:45** `buy<0`, ale baterie plná → **žádný import**.
|
**Problém (MCP run 16692, tag v20):** Od ~05:30 nabíjení z PV; v 09:15 už **98,3 %** SoC; od 09:15 masivní **export při sell<0** (−7 kW). V **11:00–12:45** `buy<0`, ale baterie plná → **žádný import**.
|
||||||
|
|||||||
Reference in New Issue
Block a user