protazeni exportniho okna az do rana pred vyrobu FVE (prichazeli sme o prilezitosti mezi pulnoci a ranem)
Some checks failed
CI and deploy / migration-check (push) Failing after 23s
CI and deploy / deploy (push) Has been skipped

This commit is contained in:
Dusan Vojacek
2026-05-26 00:00:06 +02:00
parent b03f08d3a0
commit 25c864db61
4 changed files with 228 additions and 44 deletions

View File

@@ -68,11 +68,15 @@ 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-28-pv-positive-sell-solver-v29"
PLANNER_BUILD_TAG = "2026-05-28-night-export-window-midnight-v30"
POS_SELL_PRE_NEG_SOC_SHORTFALL_PENALTY_CZK_PER_WH = 0.30
PRE_NEG_BUY_SOC_CEILING_SLACK_PENALTY_CZK_PER_WH = 0.25
PRE_NEG_BUY_EMPTY_EXPORT_SHORTFALL_PENALTY_CZK_KWH = 80.0
EVENING_PEAK_SELL_EPS_CZK_KWH = 0.05
# Noční výprodej baterie: večer (≥17h) + ráno do východu FVE (05h Prague), jedna špička přes půlnoc.
NIGHT_EXPORT_EVENING_START_HOUR = 17
NIGHT_EXPORT_MORNING_END_HOUR = 5
NIGHT_EXPORT_PV_SUNRISE_SURPLUS_W = 500.0
# Převáží terminal SoC shadow price při krátkém večerním horizontu (home-01).
EVENING_PUSH_Z_EXPORT_BONUS_CZK = 2500.0
# buy<0: preferovat import před PV A→bat (měkké; tvrdé bc_pv=0 láme bilanci s polem B).
@@ -961,28 +965,71 @@ def _pre_neg_buy_discharge_indices(
return out
def _slot_pv_surplus_w(slot: PlanningSlot) -> float:
load_w = float(slot.load_baseline_w)
pv_w = float(slot.pv_a_forecast_w) + float(slot.pv_b_forecast_w)
return max(0.0, pv_w - load_w)
def _in_night_battery_export_window(slot: PlanningSlot) -> bool:
"""
Noční okno pro večerní push / peak sell: >=17h Prague, nebo 05h (přes půlnoc).
Končí prvním slotem s významným PV přebytkem (východ FVE), ne kalendářním dnem.
"""
if _slot_pv_surplus_w(slot) > NIGHT_EXPORT_PV_SUNRISE_SURPLUS_W:
return False
h = _prague_hour(slot)
if h >= NIGHT_EXPORT_EVENING_START_HOUR:
return True
return h <= NIGHT_EXPORT_MORNING_END_HOUR
def _night_export_window_segments(slots: list[PlanningSlot]) -> list[list[int]]:
"""Souvislé úseky nočního okna v horizontu (oddělené denní pauzou / východem FVE)."""
segments: list[list[int]] = []
current: list[int] = []
for t, s in enumerate(slots):
if _in_night_battery_export_window(s):
current.append(t)
else:
if current:
segments.append(current)
current = []
if current:
segments.append(current)
return segments
def _night_peak_sell_czk_kwh(slots: list[PlanningSlot], slot_index: int) -> float:
"""Max sell v nočním úseku, do kterého slot patří (pro evening_early)."""
for seg in _night_export_window_segments(slots):
if slot_index in seg:
return max(float(slots[t].sell_price) for t in seg)
return 0.0
def _evening_peak_export_indices(
slots: list[PlanningSlot],
*,
degrad_czk_kwh: float,
evening_start_hour: int = 17,
) -> list[int]:
"""Večerní špičky per den (shodně s R__063, hour >= 17 Prague)."""
peak_by_day: dict = {}
for s in slots:
if _prague_hour(s) < evening_start_hour:
continue
d = _prague_calendar_date(s)
peak_by_day[d] = max(peak_by_day.get(d, 0.0), float(s.sell_price))
"""
Noční špičky sell: jeden peak na souvislý úsek (17h → půlnoc → ráno do východu FVE),
ne per kalendářní den (oprava 23:30 vs 00:00).
"""
_ = evening_start_hour # kompatibilita volání; okno řídí NIGHT_EXPORT_* konstanty
out: list[int] = []
for t, s in enumerate(slots):
if _prague_hour(s) < evening_start_hour:
for seg in _night_export_window_segments(slots):
if not seg:
continue
d = _prague_calendar_date(s)
peak = peak_by_day.get(d, 0.0)
if peak > 0 and float(s.sell_price) >= peak - degrad_czk_kwh:
out.append(t)
return out
peak = max(float(slots[t].sell_price) for t in seg)
if peak <= 0.0:
continue
for t in seg:
if float(slots[t].sell_price) >= peak - degrad_czk_kwh:
out.append(t)
return sorted(out)
def _evening_push_discharge_budget_wh(
@@ -1017,11 +1064,10 @@ def _evening_battery_export_push_indices(
evening_start_hour: int = 17,
) -> list[int]:
"""
Večerní push: plný ge_bat na top sell sloty (≥17h Prague).
Noční push: plný ge_bat na top sell sloty v nočním okně (≥17h + 05h do východu FVE).
Ne jeden slot — kolik slotů unese Wh rozpočet (v24), seřazených sell desc.
Kandidáti jen u denního večerního max EVENING_PEAK_SELL_EPS (úzké pásmo),
ne celé široké peakdegrad. Ráno / odpoledne řeší jiné větve solveru.
Ne jeden slot — kolik slotů unese Wh rozpočet, seřazených sell desc.
Peak sell je max v celém nočním úseku (přes půlnoc), ne per kalendářní den.
"""
if per_slot_discharge_wh <= 0.0:
return []
@@ -1569,21 +1615,12 @@ def solve_dispatch(
discharge_slot_buffer=discharge_buf_pre,
)
)
max_evening_sell_by_day: dict[object, float] = {}
for t_ev, s_ev in enumerate(slots):
if _prague_hour(s_ev) < 17:
continue
d_ev = _prague_calendar_date(s_ev)
max_evening_sell_by_day[d_ev] = max(
max_evening_sell_by_day.get(d_ev, 0.0),
float(s_ev.sell_price),
)
# Zákaz ge_bat jen *před* prvním push slotem (ne po něm — jinak terminal SoC + load
# drží energii pro 1921 h bez prodeje, home-01 téměř neexportuje).
first_evening_push_t = min(evening_push_ts) if evening_push_ts else None
if first_evening_push_t is not None:
for t_ev, s_ev in enumerate(slots):
if _prague_hour(s_ev) < 17:
if not _in_night_battery_export_window(s_ev):
continue
if t_ev >= first_evening_push_t:
continue
@@ -1591,8 +1628,7 @@ def solve_dispatch(
continue
if t_ev in evening_push_ts:
continue
d_ev = _prague_calendar_date(s_ev)
peak_sell = max_evening_sell_by_day.get(d_ev, 0.0)
peak_sell = _night_peak_sell_czk_kwh(slots, t_ev)
if float(s_ev.sell_price) < peak_sell - EVENING_PEAK_SELL_EPS_CZK_KWH:
evening_early_export_penalty_ts.add(t_ev)
last_pos_sell_pre_neg_buy = _last_non_negative_sell_before_neg_buy(