dalsi pokusy
Some checks failed
CI and deploy / migration-check (push) Failing after 11s
CI and deploy / deploy (push) Has been skipped

This commit is contained in:
Dusan Vojacek
2026-05-23 00:06:30 +02:00
parent 8845350c0b
commit a52be1b792
5 changed files with 119 additions and 12 deletions

View File

@@ -661,17 +661,26 @@ def _pre_negative_sell_export_window(
return first_neg, first_neg - 1 return first_neg, first_neg - 1
def _prague_calendar_date(slot: PlanningSlot):
dt = slot.interval_start
if dt.tzinfo is None:
dt = dt.replace(tzinfo=timezone.utc)
return dt.astimezone(ZoneInfo("Europe/Prague")).date()
def _pre_neg_peak_sell_idx( def _pre_neg_peak_sell_idx(
slots: list[PlanningSlot], slots: list[PlanningSlot],
first_neg_sell_idx: int | None, first_neg_sell_idx: int | None,
) -> int | None: ) -> int | None:
"""Index slotu s nejvyšším kladným sell před prvním sell<0.""" """Nejvyšší kladný sell před prvním sell<0 ve stejném kalendářním dni (Prague)."""
if first_neg_sell_idx is None or first_neg_sell_idx <= 0: if first_neg_sell_idx is None or first_neg_sell_idx <= 0:
return None return None
neg_day = _prague_calendar_date(slots[first_neg_sell_idx])
positive = [ positive = [
(i, float(slots[i].sell_price)) (i, float(slots[i].sell_price))
for i in range(first_neg_sell_idx) for i in range(first_neg_sell_idx)
if float(slots[i].sell_price) >= 0.0 if float(slots[i].sell_price) >= 0.0
and _prague_calendar_date(slots[i]) == neg_day
] ]
if not positive: if not positive:
return None return None

View File

@@ -211,6 +211,10 @@ def _select_charge_slots(
grid_pm += 1 grid_pm += 1
grid_filled_wh += cum grid_filled_wh += cum
for t, s in enumerate(slots):
if float(s.buy_price) < 0:
selected.add(t)
pv_layer_cap = max(charge_target_wh - grid_filled_wh, 0.0) pv_layer_cap = max(charge_target_wh - grid_filled_wh, 0.0)
pv_candidates: list[tuple[int, float, float]] = [] pv_candidates: list[tuple[int, float, float]] = []
for t, s in enumerate(slots): for t, s in enumerate(slots):
@@ -279,6 +283,33 @@ def _select_discharge_export_slots(
break break
selected.add(t) selected.add(t)
cum += per_slot_wh cum += per_slot_wh
max_sell = max((float(s.sell_price) for s in slots), default=0.0)
if max_sell > 0:
for t, s in enumerate(slots):
if float(s.sell_price) >= max_sell - degrad and float(s.sell_price) > sell_min:
selected.add(t)
first_neg = next(
(i for i, s in enumerate(slots) if float(s.sell_price) < 0),
None,
)
preneg_min_soc = min_soc_wh + max(per_slot_wh, 1000.0)
if (
first_neg is not None
and first_neg > 0
and current_soc_wh >= preneg_min_soc
):
neg_day = _prague_date(slots[first_neg])
positive = [
i
for i in range(first_neg)
if float(slots[i].sell_price) >= 0 and _prague_date(slots[i]) == neg_day
]
if positive:
peak_t = max(positive, key=lambda i: (float(slots[i].sell_price), i))
selected.add(peak_t)
return selected return selected

View File

@@ -2198,6 +2198,44 @@ class SitePowerCapTests(unittest.TestCase):
class PlannerArbitrageImprovementsTests(unittest.TestCase): class PlannerArbitrageImprovementsTests(unittest.TestCase):
"""Regrese: záporný buy, peak sell před sell<0, večerní export cap.""" """Regrese: záporný buy, peak sell před sell<0, večerní export cap."""
def test_pre_neg_peak_ignores_earlier_day_in_horizon(self) -> None:
"""Horizont přes půlnoc: peak je na dni záporného sell, ne včerejší večer."""
base = datetime(2026, 5, 22, 18, 0, tzinfo=timezone.utc)
slots = [
PlanningSlot(
interval_start=base,
buy_price=4.0,
sell_price=4.6,
pv_a_forecast_w=0,
pv_b_forecast_w=0,
load_baseline_w=1000,
ev1_connected=False,
ev2_connected=False,
),
PlanningSlot(
interval_start=base + timedelta(hours=12),
buy_price=4.0,
sell_price=3.06,
pv_a_forecast_w=3000,
pv_b_forecast_w=0,
load_baseline_w=1000,
ev1_connected=False,
ev2_connected=False,
),
PlanningSlot(
interval_start=base + timedelta(hours=12, minutes=45),
buy_price=0.5,
sell_price=-0.1,
pv_a_forecast_w=5000,
pv_b_forecast_w=0,
load_baseline_w=1000,
ev1_connected=False,
ev2_connected=False,
),
]
first_neg = 2
self.assertEqual(_pre_neg_peak_sell_idx(slots, first_neg), 1)
def test_pre_neg_peak_idx_is_highest_positive_sell(self) -> None: def test_pre_neg_peak_idx_is_highest_positive_sell(self) -> None:
base = datetime(2026, 5, 23, 4, 0, tzinfo=timezone.utc) base = datetime(2026, 5, 23, 4, 0, tzinfo=timezone.utc)
slots = [ slots = [

View File

@@ -78,7 +78,9 @@ declare
v_grid_slots_pm int := 0; v_grid_slots_pm int := 0;
v_acquisition_cutoff timestamptz; v_acquisition_cutoff timestamptz;
v_first_neg_sell_ord int; v_first_neg_sell_ord int;
v_first_neg_prague_date date;
v_pre_neg_peak_sell_ord int; v_pre_neg_peak_sell_ord int;
v_preneg_export_min_soc_wh numeric;
v_max_sell_czk_kwh numeric; v_max_sell_czk_kwh numeric;
v_charge_acquisition numeric; v_charge_acquisition numeric;
v_est_grid_wh numeric; v_est_grid_wh numeric;
@@ -518,6 +520,23 @@ begin
end loop; end loop;
end if; end if;
-- První záporný výkup v horizontu (od p_from = „budoucí“ sloty od replanu).
-- Peak před ním a acquisition cutoff se vážou na STEJNÝ kalendářní den (Prague),
-- ne na včerejší večer v tomže horizontu.
select min(wk.slot_ord)
into v_first_neg_sell_ord
from _ems_plan_slot_wk wk
where wk.sell_price < 0;
if v_first_neg_sell_ord is not null then
select (wk.interval_start at time zone 'Europe/Prague')::date
into v_first_neg_prague_date
from _ems_plan_slot_wk wk
where wk.slot_ord = v_first_neg_sell_ord;
end if;
v_preneg_export_min_soc_wh := v_min_soc_wh + greatest(v_per_slot_discharge_wh, 1000::numeric);
-- discharge-export mask -- discharge-export mask
if v_discharge_buf <= 0 then if v_discharge_buf <= 0 then
update _ems_plan_slot_wk wk set allow_discharge_export = true; update _ems_plan_slot_wk wk set allow_discharge_export = true;
@@ -565,18 +584,17 @@ begin
); );
end if; end if;
-- Před prvním sell<0: export baterie v lokálním maximu kladného sell (ne jen v posledních slotech). -- Před prvním sell<0 téhož dne: export v lokálním max kladného sell (např. 07:00, ne včerejší 20:45).
select min(wk.slot_ord) if v_first_neg_sell_ord is not null
into v_first_neg_sell_ord and v_first_neg_prague_date is not null
from _ems_plan_slot_wk wk and p_current_soc_wh >= v_preneg_export_min_soc_wh
where wk.sell_price < 0; then
if v_first_neg_sell_ord is not null then
select wk.slot_ord select wk.slot_ord
into v_pre_neg_peak_sell_ord into v_pre_neg_peak_sell_ord
from _ems_plan_slot_wk wk from _ems_plan_slot_wk wk
where wk.slot_ord < v_first_neg_sell_ord where wk.slot_ord < v_first_neg_sell_ord
and wk.sell_price >= 0 and wk.sell_price >= 0
and (wk.interval_start at time zone 'Europe/Prague')::date = v_first_neg_prague_date
order by wk.sell_price desc, wk.slot_ord order by wk.sell_price desc, wk.slot_ord
limit 1; limit 1;
@@ -587,12 +605,23 @@ begin
end if; end if;
end if; end if;
-- Vážená acquisition cena zásoby (grid + FVE opportunity) jen pro sloty PŘED prvním -- Acquisition: grid nabíjení před prvním exportem ve STEJNÝ den jako záporné výkupní okno
-- plánovaným exportem z baterie — nepočítá nákup po večerním/nočním vybití do sítě. -- (ne dřívější večerní export v horizontu rolling replanu).
select min(wk.interval_start) select min(wk.interval_start)
into v_acquisition_cutoff into v_acquisition_cutoff
from _ems_plan_slot_wk wk from _ems_plan_slot_wk wk
where wk.allow_discharge_export; where wk.allow_discharge_export
and (
v_first_neg_prague_date is null
or (wk.interval_start at time zone 'Europe/Prague')::date = v_first_neg_prague_date
);
if v_acquisition_cutoff is null then
select min(wk.interval_start)
into v_acquisition_cutoff
from _ems_plan_slot_wk wk
where wk.allow_discharge_export;
end if;
-- Acquisition: jen grid vrstva B (ne odpolední FVE z vrstvy A) před 1. exportem. -- Acquisition: jen grid vrstva B (ne odpolední FVE z vrstvy A) před 1. exportem.
select select

View File

@@ -39,7 +39,7 @@
- měkký cíl na konci 24h přes `_soc_security_profile` + tvrdé dvouúrovňové pravidlo výše. - měkký cíl na konci 24h přes `_soc_security_profile` + tvrdé dvouúrovňové pravidlo výše.
- **Dynamická ekonomická podlaha (fáze 2):** - **Dynamická ekonomická podlaha (fáze 2):**
- `_dynamic_arb_floor_wh_series`: podle součtu FVE výkonu v dalších ~8 h (`ARB_LOOKAHEAD_SLOTS`) se `arb_floor_wh[t]` posouvá mezi `min_soc_wh` a rezervou z DB silné očekávané slunce ji sníží (ráno / po obloze); vynutit konstantní chování lze `battery.disable_dynamic_arb_floor=True` jen pro testy / ladění. - `_dynamic_arb_floor_wh_series`: podle součtu FVE výkonu v dalších ~8 h (`ARB_LOOKAHEAD_SLOTS`) se `arb_floor_wh[t]` posouvá mezi `min_soc_wh` a rezervou z DB silné očekávané slunce ji sníží (ráno / po obloze); vynutit konstantní chování lze `battery.disable_dynamic_arb_floor=True` jen pro testy / ladění.
- **Výběr exportních slotů (`allow_discharge_export`):** `ems.fn_load_planning_slots_full` označí jen sloty, kde smí solver **úmyslně** vybíjet baterii do sítě (SELL). Výběr je **globálně** podle `sell_price desc` (ne AM/PM 50/50), doplněno o: (1) **všechny sloty** s `sell ≥ max(sell) degradation` (večerní špičky vždy exportovatelná), (2) **lokální maximum kladného `sell` před prvním `sell < 0`** (ráno prodat v nejlepším slotu, ne až v posledním před záporným oknem). **Spot nákup:** `sell_price > ref_buy + degradation_cost_czk_kwh` (`ref_buy` = **min `buy` horizontu**). **Fixní nákup** (`purchase_pricing_mode = fixed`): `sell_price > degradation_cost_czk_kwh`. V `solve_dispatch` (AUTO): **`ge_pv`** / **`ge_bat`**; v **high-sell** exportních slotech měkká penalizace **`export_shortfall`** (`ge` pod `max_export_power_w`, koef. **`PEAK_EXPORT_SHORTFALL_PENALTY_CZK_KWH`** v `planning_engine.py`). **Kotva před `sell < 0`:** SoC ≤ planner floor v **posledním slotu před prvním `sell < 0`**; export baterie smí být v **lokálním maximu kladného sell** (`allow_discharge_export` + `ge_bat` push v `solve_dispatch`). Mimo exportní sloty: **`ge_bat = 0`**; **`bc_gi = 0`** mimo masku, **výjimka `buy < 0`** (grid nabíjení i bez `allow_charge`). **`deye_physical_mode`** = PASSIVE kromě CHARGE/SELL. - **Výběr exportních slotů (`allow_discharge_export`):** `ems.fn_load_planning_slots_full` označí jen sloty, kde smí solver **úmyslně** vybíjet baterii do sítě (SELL). Výběr je **globálně** podle `sell_price desc` (ne AM/PM 50/50), doplněno o: (1) **všechny sloty** s `sell ≥ max(sell) degradation` (večerní špičky vždy exportovatelná), (2) **lokální maximum kladného `sell` před prvním `sell < 0` ve stejném kalendářním dni (Europe/Prague)** — horizont od `p_from` může zahrnovat víc dní (rolling večer + ráno); peak **není** včerejší večerní špička. Povoleno jen pokud `p_current_soc_wh ≥ min_soc + 1 slot discharge` (SoC). **`charge_acquisition`:** vážený `buy` u `allow_grid_charge` před **prvním exportem téhož dne** jako záporné výkupní okno (ne před včerejším exportem v horizontu). **Spot nákup:** `sell_price > ref_buy + degradation_cost_czk_kwh`. V `solve_dispatch` (AUTO): **`ge_pv`** / **`ge_bat`**; v **high-sell** exportních slotech měkká penalizace **`export_shortfall`**. **Kotva před `sell < 0`:** SoC ≤ planner floor v posledním slotu před prvním `sell < 0`; **`ge_bat` push** v peak slotu (Python, shodný den). Mimo exportní sloty: **`ge_bat = 0`**; **`bc_gi = 0`** mimo masku, **výjimka `buy < 0`**. **`deye_physical_mode`** = PASSIVE kromě CHARGE/SELL.
- **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í —