@@ -50,10 +50,10 @@ DEFAULT_PLANNER_DISCHARGE_RELAX_PREWINDOW_SLOTS = 8
|
||||
# Penalizace je v Kč/Wh (např. 0.20 = 200 Kč/kWh). Musí být dost velká, aby přebila
|
||||
# bezpečnostní SoC buffer + terminal shadow cenu a solver skutečně „dovylil“ před sell<0.
|
||||
PRENEG_SELL_SOC_ANCHOR_SLACK_PENALTY_CZK_PER_WH = 0.20
|
||||
# Penalita za překročení SoC capu před prvním buy<0 slotem. Musí být VÝRAZNĚ VYŠŠÍ
|
||||
# než marginal arbitrage (peak_sell − avg_neg_buy ~ 5 Kč/kWh), aby LP cap dodržoval
|
||||
# i přes alternativní cesty.
|
||||
PRE_NEG_BUY_SOC_SLACK_PENALTY_CZK_PER_WH = 0.05
|
||||
# Penalita za překročení SoC capu před prvním buy<0 slotem. Měla by být VYŠŠÍ než
|
||||
# alternativní marginal arbitrage (acquisition - avg_neg_buy ~ 1 Kč/kWh) aby LP
|
||||
# preferoval buy<0 nabíjení před ranním PV.
|
||||
PRE_NEG_BUY_SOC_SLACK_PENALTY_CZK_PER_WH = 0.005
|
||||
PEAK_EXPORT_SHORTFALL_PENALTY_CZK_KWH = 80.0
|
||||
# Měkký tlak: v okně sell<0 + block_export využít PV přebytek do baterie (ne curtail).
|
||||
PV_CHARGE_SHORTFALL_PENALTY_CZK_KWH = 120.0
|
||||
@@ -68,7 +68,7 @@ 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áž).
|
||||
EXTREME_BUY_DUMP_PREWINDOW_SLOTS = 12
|
||||
NEG_SELL_BAT_DUMP_SHORTFALL_PENALTY_CZK_KWH = 80.0
|
||||
PLANNER_BUILD_TAG = "2026-05-27-hard-pv-mask-v15"
|
||||
PLANNER_BUILD_TAG = "2026-05-27-pre-neg-buy-soc-cap-v14"
|
||||
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_MAX_CLAMP = 1.5 # horní limit korekčního faktoru
|
||||
@@ -1204,33 +1204,6 @@ def solve_dispatch(
|
||||
discharge_export_slots = {
|
||||
t for t, s in enumerate(slots) if s.allow_discharge_export
|
||||
}
|
||||
# Pokud je v horizontu buy<0 okno a startovní SoC je nad cap (rezervační target),
|
||||
# potřebujeme cestu jak baterii do okna vybít — rozšířit discharge_export okno o
|
||||
# vysoké-sell sloty před prvním buy<0. R__063 v noci typicky nepovoluje export
|
||||
# (allow_discharge_export=false), ale tady chceme zaplatit za vyklízení baterky,
|
||||
# abychom potom v buy<0 okně nasáli max z OTE záporných cen.
|
||||
_neg_buy_indices_pre = [
|
||||
t for t, s in enumerate(slots) if float(s.buy_price) < 0
|
||||
]
|
||||
_first_neg_buy_idx_pre = _neg_buy_indices_pre[0] if _neg_buy_indices_pre else None
|
||||
_usable_cap_pre = float(battery.soc_max_wh) - float(min_soc_wh)
|
||||
_cap_target_raw_wh_pre = float(min_soc_wh) + 0.15 * _usable_cap_pre
|
||||
if (
|
||||
_first_neg_buy_idx_pre is not None
|
||||
and _first_neg_buy_idx_pre > 0
|
||||
and float(current_soc_wh) > _cap_target_raw_wh_pre
|
||||
):
|
||||
# Práh sell: alespoň 80 % maxima sell v okně před buy<0 (peak hours).
|
||||
_pre_sells = [
|
||||
float(slots[t].sell_price) for t in range(_first_neg_buy_idx_pre)
|
||||
if float(slots[t].sell_price) > 0
|
||||
]
|
||||
if _pre_sells:
|
||||
_sell_peak_pre = max(_pre_sells)
|
||||
_sell_thr_pre = 0.8 * _sell_peak_pre
|
||||
for t in range(_first_neg_buy_idx_pre):
|
||||
if float(slots[t].sell_price) >= _sell_thr_pre:
|
||||
discharge_export_slots.add(t)
|
||||
# SELF_SUSTAIN dřív vynucoval ge[t] == 0, což umí udělat MILP infeasible v okamžiku, kdy:
|
||||
# - baterie je na max SoC (nelze nabíjet),
|
||||
# - PV pole B není curtailable,
|
||||
@@ -1296,17 +1269,24 @@ def solve_dispatch(
|
||||
first_neg_buy_idx = neg_buy_indices[0] if neg_buy_indices else None
|
||||
pre_neg_buy_soc_cap_wh: float | None = None
|
||||
if first_neg_buy_idx is not None and first_neg_buy_idx > 0:
|
||||
# Cap = min_soc + 15 % usable kapacity (~25 % SoC celkem pro home-01).
|
||||
# Strukturální hodnota: do buy<0 okna LP musí dorazit s prakticky prázdnou
|
||||
# baterií, aby tam mohl maximálně nasát z OTE záporných cen + PV (+ následné
|
||||
# sell<0 sloty taky doplní PV → odpolední přebytek se nevyhodí v curtail).
|
||||
# Spočítat max nabíjecí kapacitu v buy<0 oknu (per slot ≤ max_charge_power_w,
|
||||
# PV přebytek + grid import). Konzervativně: jen sloty se buy<0, nikoli okolní.
|
||||
neg_window_capacity_wh = 0.0
|
||||
for t in neg_buy_indices:
|
||||
s_neg = slots[t]
|
||||
pv_sur = max(
|
||||
0.0,
|
||||
float(s_neg.pv_a_forecast_w) + float(s_neg.pv_b_forecast_w) - float(s_neg.load_baseline_w),
|
||||
)
|
||||
slot_charge_pot_w = min(
|
||||
float(battery.max_charge_power_w),
|
||||
pv_sur + float(grid.max_import_power_w),
|
||||
)
|
||||
neg_window_capacity_wh += slot_charge_pot_w * INTERVAL_H * float(battery.charge_efficiency)
|
||||
# Konzervativní headroom 10 % kapacity (PV forecast error, load fluctuation)
|
||||
usable_capacity = float(battery.soc_max_wh) - float(min_soc_wh)
|
||||
pre_neg_buy_soc_cap_wh = float(min_soc_wh) + 0.15 * usable_capacity
|
||||
# Pokud start SoC je už nad capem (zdědili jsme „plnou" baterii od minulého
|
||||
# plánu / replanu), zvedneme cap na current SoC. Constraint je pak „nepřibývat"
|
||||
# místo „prudce klesnout" — LP není infeasible / megapenalty když nemá kam vybít.
|
||||
if float(current_soc_wh) > pre_neg_buy_soc_cap_wh:
|
||||
pre_neg_buy_soc_cap_wh = float(current_soc_wh)
|
||||
reserve_target_wh = min(neg_window_capacity_wh, usable_capacity * 0.9)
|
||||
pre_neg_buy_soc_cap_wh = max(float(min_soc_wh), float(battery.soc_max_wh) - reserve_target_wh)
|
||||
last_neg_sell_by_prague_date: dict[object, int] = {}
|
||||
for t_ln, st_ln in enumerate(slots):
|
||||
if float(st_ln.sell_price) < 0:
|
||||
|
||||
Reference in New Issue
Block a user