欢迎光临~霍克(HAWKER)蓄电池_霍克锂电池网站
服务热线
全国客服热线:

15313702523

HAWKER蓄电池-官方网站

HAWKER蓄电池一种用于RE100电网接入的严谨双尺度预测优化框架:将城市雨水调蓄池升级改造为具有混合储能功能的低水头变速抽水蓄能电站

摘要

全球范围内对RE100合规性的要求加速了分布式光伏(PV)发电的部署,这对土地稀缺的城市电网的惯性与电压稳定性构成了重大挑战。为应对这些挑战,本研究提出了一种基础设施升级利用范式,将城市地下雨水调蓄池转化为去中心化的低水头变速抽水蓄能(VS-PHES)系统,并与飞轮(FESS)和锂离子电池(BESS)混合储能系统(HESS)相耦合。该设施由两阶段分层能量管理系统(EMS)控制:第一阶段采用混合整数线性规划(MILP)进行考虑电池退化惩罚的日前宏观调度;第二阶段执行双速率非线性模型预测控制(NMPC)框架(100毫秒监督/10毫秒本地回路),用于亚秒级微暂态电压和频率调节。通过结合经Symlet-4小波预滤波优化的1D-CNN预测特征提取模型,该系统实现了97.85%的暂态状态分类准确率。一种故障安全、以防洪控制为优先的状态切换逻辑,确保了城市在强降水事件期间的防洪安全。该逻辑在经过修改的IEEE 33节点放射状电网以及实际的22.9千伏134节点韩国城市配电网(在局部5兆瓦光伏掉落事故下,总装机容量为20兆瓦)中得到验证,协调的3.5兆瓦混合设施将瞬时电压跌落从0.931标幺值缓解至0.952标幺值,并将频率最低点控制在59.35赫兹。一项20年全寿命周期成本(LCC)分析表明,飞轮储能系统(FESS)屏蔽可节省高达100万美元的电池储能系统(BESS)更换开支,并由支持私人可再生能源100(RE100)投资的动态节点电价机制提供支持。9100万美元的电池储能系统(BESS)更换开支,由支持私人可再生能源100(RE100)投资的动态节点电价机制提供支持。
关键词:
RE100; 变速抽水蓄能 (VS-PHES); 混合储能 (HESS); 模型预测控制 (NMPC); 飞轮储能 (FESS); 雨水调蓄池; 微电网稳定性

 

1. 引言

在RE100倡议等国际可持续发展框架的推动下,向碳中和能源架构的快速转型已从根本上改变了现代城市配电网的拓扑结构[1]。分布式能源(DERs)的大规模扩展,特别是分布式光伏(PV)阵列,提供了零排放的有功功率,但也引入了关键的运行脆弱性[2,3]。由于分布式光伏系统通过静态电力电子逆变器与电网连接,它们缺乏传统同步热力发电单元固有的物理转动惯量[4]。因此,高渗透率的光伏环境会经历系统总转动惯量的剧烈下降。在这些低惯量条件下,局部放射状配电馈线极易受到极端频率变化率(RoCoF)、深度频率极值以及由突发云遮挡事件或局部馈线跳闸触发的瞬态电压跌落的影响[5]. 图1阐明了当多个分布式光伏(PV)阵列因云遮挡事件而经历突发性同时掉落时,时空相关性所导致的严重电网脆弱性和广域波动扰动。它表明,传统的基于平均值和历史确定性控制的方法无法准确预测和防御这些多兆瓦级的突发发电量掉落。因此,该图强调了迫切需要先进的概率和分布式控制框架,以确保城市配电网的稳定运行。
Electronics 15 05037 g001
图1。在高分布式光伏渗透率条件下,电网脆弱性、低惯量退化以及严重的暂态电压违规。
A主要瓶颈在于城市中心区域专用电网稳定公用设施安装的城市土地极度匮乏且成本高昂[6]。为了同时解决物理空间限制和电力稳定性挑战,本文提出了一种基础设施升级改造框架。人口稠密的都市区拥有广泛的地下雨水调蓄池和排水隧道网络,这些设施旨在缓解强降水期间的突发性洪水[7]。这些庞大的土木工程资产在一年中超过90%的时间里完全处于干燥、闲置和未利用状态。通过将这些地下水库划分为上下水力室,或将其与相邻的低海拔城市水道相连,它们可以被重新设计为分散式低水头变速抽水蓄能(VS-PHES)设施[8].
图2展示了与调速器阀门启闭相关的非弹性水锤压力瞬变和水力机械响应延迟。具体而言,上方的子图展示了在t = 2.0 s时,导叶开度(G)快速阶跃关闭事件所触发的引水管流量(Q)瞬态响应和水锤峰值(2.25 p.u.)。同时,下方的子图描绘了非线性弹性水锤动力学,突显了水轮机入口处有效水头(H)在最终稳定前出现的剧烈压力激增和波动。
Electronics 15 05037 g002
图2。与调速器阀门开闭相关的非弹性水锤压力瞬变和水力机械响应延迟。
然而,将水力机械设备整合到微电网实时调节中会引入严重的流体动力学约束。由于水锤管内水柱的惯性、导叶调节延迟以及缓解破坏性水锤现象的必要性,常规和变速抽水蓄能机组存在数百毫秒到数秒的固有水力机械响应滞后[9,10]. 这种水力延迟在发生大规模可再生能源掉落后,紧随其后的极易受损的子周期瞬态窗口内,造成了关键的控制真空。为了弥补这一动态运行缺口,本研究开发了一种多性向混合储能系统(HESS),直接部署在升级改造的VS-PHES设施内[11]。通过将擅长瞬时、亚周期动能功率注入的飞轮储能系统(FESS)与高能量密度锂离子电池储能系统(BESS)相结合,建立了一个有机平衡的多时间尺度防御框架[12]。这种结构上的同址部署有效消除了重型流体基础设施与超快速电网动态之间的响应滞后不匹配。虽然最近的研究探索了Volt-VAR和Volt-Watt控制优化策略,以管理高渗透率光伏径向馈线中的电压稳定性[12],这些传统方法缺乏主动动能惯性;因此,本研究通过提出一种利用VS-PHES和HESS多性向框架的稳健基础设施升级范式,克服了这些局限性
图3展示了混合储能系统(HESS)中用于频率防御和电池寿命保护的多时间尺度有机角色共享机制。具体而言,顶部的子图展示了功率分配策略,其中电网总需求的频率高瞬态分量由飞轮储能系统(FESS)吸收,而平滑的大宗能量需求则由锂离子电池储能系统(BESS)处理。中间的子图显示了FESS的动力学特性及其对应高频功率波动的转动角速度变化。最后,底部的子图验证了FESS的屏蔽作用成功地将BESS的荷电状态(SOC)维持在安全运行范围内,并显著减轻了电池单体随时间推移的累积退化。
Electronics 15 05037 g003
图3。用于频率防御和电池寿命保护的多时间尺度有机角色共享机制。
基于Meng等人[13]阐述的多频率导纳和解耦原理,开发了一种高速频率解耦滤波器拓扑结构,以实现储能资产间功率组件的无缝分配。
近期关于电力系统稳定性的文献日益强调将基于逆变器的资源(IBRs)整合到弱配电网中所面临的挑战。正如近期具有里程碑意义的研究所示,使用多目标机器学习算法评估和预测IBRs整合弱电网中的数据驱动系统强度水平,对于预判电压崩溃和设计主动惯性补偿策略至关重要[14]。此外,混合系统架构的最新进展表明,计算高效的线性近似对于复杂的潮流动态至关重要。例如,Zhang等人[15]成功提出了非对称耦合电气化铁路-电网混合系统的功率流模型稳健线性化方法,在保持物理准确性的同时实现了显著的计算性能提升。基于这些混合功率流线性化的概念基础,我们的研究通过引入一种两阶段分层EMS架构扩展了这一范式,该架构利用MILP实现分布式变速抽水蓄能的宏观经济集成。当IBR渗透率上升时,传统的短路比(SCR)指标无法捕捉局部动态阻抗偏移,因此需要能够进行多尺度动态响应的实时预测优化框架[14].

具体方法论贡献

为解决传统定速水电机组、非协调储能资产以及弱电网电压不稳定的局限性,本文提出了五项关键的科学与工程贡献:
  • 民用基础设施升级改造模型:一种现实的结构框架,可在无需额外土地征收成本的情况下,将未利用的城市地下雨水调蓄池转化为分散式的2.0 MVA低水头VS-PHES设施。
  • 双速率分层EMS架构:一种多时间尺度预测优化系统,包含第一阶段日前混合整数线性规划(MILP)宏观经济调度,以及第二阶段双速率非线性模型预测控制(NMPC)(100毫秒监督/10毫秒局部逆变器控制),并结合一维卷积神经网络(1D-CNN)进行实时特征预测。
  • 水力机械物理公式:综合数学建模,包含非线性非弹性水锤压力瞬变、变速水轮机动力学以及电液调速器延迟。
  • BMS硬件集成与老化电池降解建模:完整的从属BMS硬件规格,具备CAN 2.0B通信功能,并采用B0005 431次循环拐点退化建模,用于保守的LCC评估。
  • 在实际城市馈线和硬件在环架构上的验证:在标准IEEE 33节点系统和实际22.9 kV 134节点韩国配电网上进行了广泛验证,并辅以控制硬件在环(CHIL)延迟基准测试。

为了突出本研究的科学与工程贡献,表1对所提出的框架与现有文献进行了全面的方法论比较。与传统的混合能源管理系统和传统城市抽水蓄能配置不同,该架构利用未使用的地下雨水调蓄池,实现了零土地成本占用。此外,它引入了由双速率非线性模型预测控制(NMPC)架构控制的三重混合储能耦合(VS-PHES、BESS和FESS)。该表还强调了非线性水锤动力学的集成、在实际134节点网络上的先进弱电网验证方法,以及一个包含生命周期成本(LCC)退化屏蔽和动态节点电价的稳健经济模型。
表1。与现有文献的方法论比较。

2. 升级改造的多能互补基础设施动态建模

2.1. 非线性水锤动态与VS-PHES响应

如图所示图2阐述了在拟建低水头变速抽水蓄能电站(VS-PHES)快速调速器阀门开/关过程中,引水管内的水力瞬态响应和水锤压力变化。由于低水头市政雨水调蓄池基础设施所需的大体积流量,流体柱惯性显著延长,导致水启动时间常数延长(Tw)。引水管内非线性非弹性水锤压力动力学由以下微分方程控制:
Tw·[dq(t)/dt] = h(t) − h0 − hf(t)

其中q(t)是归一化引水管流量(标幺值),h(t)是水轮机入口处的归一化有效水头(标幺值),h0是静扬程(标称12.5米),hf(t) = fr·q(t)|q(t)|表示水力摩擦水头损失,且Tw是水的起始时间常数。Tw定义如下:
) = 2.15 sw = (Lp·Qbase)/(gacc · Ap · Hbase) = 2.15 s

调速器导叶执行器动力学和机械功率输出Phydro(t)建模如下:
dgv(t)/dt = [ġg(ġ) − gv(t)]/Tg, |ġv(t)| ≤ ġv,max

) = ηhydro(·ρ·)·ωturb·ρ·gacc·Heff(t)·Q(t)·ωturb(t)

其中gv(t)是导叶开度位置(0 ≤gv ≤ 1), Tg= 0.5 s是电液执行器滞后,ug(t)是控制信号,η是变速水力效率。turb

2.2. 多性向混合储能、FESS损耗与电池退化建模

机械式FESS转子的结构铭牌最大限制为ω= 16,000 RPM(≈1675.5 rad/s)。然而,为了保持高往返效率并减轻机械疲劳,其运行状态下的荷电状态(SOC)包络线被限制在目标速度区间内,介于ω= 370 rad/s(3533 RPM)和ω= 430 rad/s(4106 RPM)之间。同址FESS吸收高频瞬态功率脉冲。其瞬时有功功率maxmin,opmax,opPFESS(t) 包含机械风阻/轴承损耗Ploss(ω) 和变流器效率 η:flyconv
·ωFESS((ω) = [)]·η·ωfly(t)·(dωfly(t)/dt) − Ploss(ωfly)]·ηconvsign

(ωlossωfly) = ω0 + c1ωfly + c2ωfly2

FESS在某一时刻的瞬时有功功率输出/输入FESS([kW或MW]。): Instantaneous active power output/input of the FESS at time t [kW or MW].
飞轮转子的转动惯量 [kg·m: Rotor moment of inertia of the flywheel [kg·m2].
ωfly(飞轮转子在时间):[rad/s]的旋转角速度。t [rad/s].
dωfly(t)/dt:飞轮转子在时间的旋转角加速度或减速度 [rad/s2].
ηconv:双向机电功率转换器效率 [标幺值]。
  • 此处,J= 45.8 kg·m是转动惯量,放电时符号为1,充电时为-1。总可用动能2Eusable,FESS工作频带内表示如下:
    ·(ωusable,FESS = 0.5·− ω)·ηmax,op2 − ωmin,op2)·ηconv
为了精确计算通过飞轮储能系统(FESS)为电池储能系统(BESS)提供调频服务所带来的经济效益,非线性电化学电池衰减模型会累积容量损失Cdeg(t)并使用安时吞吐量和放电深度(DoD)应力惩罚:
[μ·|deg((τ)|/(2·) = ∫0t·(DoD(τ)))] dτbatt(τ)|/(2·Enom·(DoD(τ))κ)] dτ

截至时间的累积电池容量损失deg(t): Accumulative battery capacity loss up to time t.
μ:反映电化学降解特性的降解权重因子。
(τ):battτ时刻的电池瞬时电流。 Instantaneous battery current at time τ.
BESS的额定能量容量(3.0 MWh)。nom: Nominal energy capacity of the BESS (3.0 MWh).
DoD(τ):τ时刻电池的放电深度。
κ:反映非线性退化应力的DoD惩罚应力指数(κ = 1.4)。
τ:随时间连续累积的退化权重的积分变量。
为了确保经济建模的保守性,电池退化参数使用NASA的B0005锂离子电池老化数据集在第431个循环时进行了校准。第431个循环对应于关键的“拐点”转变,此时SOH下降至80%(1.60 Ah容量),而内阻增加了41.2%。在此老化状态下对退化模型进行参数化,可防止在严苛的亚秒级功率注入工作循环下,对电池储能系统(BESS)寿命产生乐观的过度估计。int increasing by 41.2%. Parameterizing the degradation model at this aged state prevents optimistic overestimation of BESS lifetime under harsh sub-second power injection duty cycles.
BMS硬件架构:电池组由工业主从式电池管理系统(BMS)控制。从模块采用Analog Devices LTC6813-1 18通道精密电池监测IC,通过电隔离的ISO-SPI链路进行接口。主动均衡采用电容电荷转移电路实现,可提供高达2.0 A的均衡电流。电池电压、温度(每2个电池单元配备1个NTC热敏电阻传感器)和电池组电流(霍尔效应传感器)以1 kHz的频率进行采样,并每10毫秒通过高速CAN 2.0B总线(500 kbps)传输至Stage 2b控制器。硬接线安全继电器独立于监控软件,在<5毫秒内强制执行过流(>3.5 C)、过压(>4.25 V/单元)和热跳闸(>55 °C)阈值。

2.3. 市政防洪优先级状态切换约束

为确保市政安全,该混合设施采用故障安全状态切换逻辑,将实时降雨强度映射至Rrain(t)。令Ψ(t) ∈ {0, 1} 表示运行状态变量:
Ψ() = 1 如果(电网EMS激活)Rrain(t) < Rthreshold   (Grid EMS Active)
Ψ() = 0 如果(紧急防洪减灾覆盖)Rrain(t) ≥ Rthreshold   (Emergency Flood Mitigation Override)

Ψ(二进制操作状态变量 (Ψ():) = 1 表示电网支持型EMS运行;Ψ() = 0 表示紧急市政防洪覆盖)。) = 1 for grid-support EMS operation; Ψ(t) = 0 for emergency municipal flood mitigation override).
实时测得的降雨强度为rain([mm/h]。): Real-time measured rainfall intensity at time t [mm/h].
触发强制防洪状态切换的关键降水阈值 [mm/h]。threshold: Critical precipitation threshold triggering mandatory flood control state switch [mm/h].
当Ψ()=0时,所有电网辅助调度设定值清零(*=0),且水泵以最大功率启动,以排空下方的雨水蓄水池。Phydro* = 0, PBESS* = 0, PFESS* = 0), and hydro pumps activate at maximum capacity to drain the lower stormwater chamber.

3. 两阶段分层优化框架

为了阐明所提混合系统的运行集成,图4提供了双尺度分层优化与预测控制框架的完整流程图。该控制架构分为两个不同的时间与结构层,以实现宏观经济调度与快速稳定控制的解耦。图4阐明了所提出的动态终端库容耦合约束在双尺度分层EMS框架下对水库水位跟踪性能的关键影响。如红色轨迹所示,传统的非耦合实时MPC模型无法维持长期水量平衡,导致在终端边界处出现严重的调度崩溃和不可行性($相比之下,所提出的耦合MPC框架(蓝色轨迹)成功实现了终端时域平衡,确保实时运行轨迹严格遵循日前宏观调度参考(黑色虚线),同时保障了市政用水的可用性。$). In contrast, the proposed coupled MPC framework (blue trajectory) successfully enforces terminal horizon balance, ensuring that the real-time operational trajectory strictly adheres to the day-ahead macro-scheduler reference (black dashed line) while safeguarding municipal water availability.
Electronics 15 05037 g004
图4。所提出的双尺度分层优化与预测控制框架的完整流程图。

3.1. 第一阶段:日前经济能源调度(MILP)

第一阶段在滚动24小时的时间跨度内,以宏观时间步长(ΔT= 1小时)执行混合整数线性规划(MILP)优化程序,以建立基础能源调度基准。上层目标函数严格最小化净电力进口成本,纳入基准运维开销,并内化长期电池退化惩罚:
min ∑t=124 []·Δgrid,t·Pimport,t + Com·(Phydro,tgen + Phydro,tpump) + Ccell_replace·Cdeg,t]·ΔT

受以下严格约束条件限制:
系统功率平衡:
Pimport,t + PPV,t + Phydro,tgen + PBESS,tdis + PFESS,tdis = Pload,t + Phydro,tpump + PBESS,tch +
PFESS,tch    ∀t

BESS动态限制与互斥:
SOCBESS(+ 1) = SOC) + [(ηBESS(/η]·Δch PBESS,tch − PBESS,tdis/ηdis)/Enom]·ΔT

 SOCBESSmin≤ SOCBESS() ≤ SOC) ≤ SOCBESSmax,      0 ≤ PBESS,tch ≤ uBESS,tchPBESSmax,
  0 ≤ PBESS,tdis ≤ uBESS,tdisPBESSmax

uBESS,tch + uBESS,tdis ≤ 1,      u ∈ {0, 1}

水库蓄水量范围:
(t + 1) =res))·ΔVres(t) + (Qpump(t) − Qturb(t))·ΔT

Vresmin ≤ Vres(t) ≤ Vresmax − Vflood_reserve

3.2. 第二阶段:双速率实时非线性预测控制(NMPC)

第二阶段采用一种新型双速率架构,以毫秒级时间步长运行,以弥合快速电力电子设备与慢速流体动力学之间的时间尺度不匹配:
  • 阶段2a(监督式NMPC,Δ= 100 ms):MPC该阶段在滚动预测时域= 20(2.0 s)和控制时域= 5(0.5 s)内运行。它使用CasADi(v3.5.5)与IPOPT(v3.12.13)内点求解器接口,并调用MUMPS线性求解器进行求解。p = 20 (2.0 s) and control horizon Nc = 5 (0.5 s). It is solved using CasADi (v3.5.5) interfaced with the IPOPT (v3.12.13) interior-point solver using the MUMPS linear solver.
  • 阶段2b(局部逆变器控制器,Δ= 10 ms):localFESS变流器上的高带宽电流矢量控制器在扰动发生后<20 ms内提供自主的一次下垂响应。 A high-bandwidth current vector controller on the FESS converter provides autonomous primary droop response within <20 ms of a perturbation.

第二阶段a的NMPC优化问题公式如下:
min ∑k=1Np [(Δf(Δf(t + k))2 + wv(ΔVi(t + k))2 + wtrack ||Pactual(t + k) − Pref(t + k)||2] +
∑k=0||Δc−1 ||Δu(t + k)||R2

状态向量: ) = [Δ(), Δ), SOC), SOC(t), ΔVi(t), SOCFESS(t), SOCBESS(t), Qpenstock(t)]T
控制输入向量: u(t) = [PFESS(t), PBESS(t), PVS-PHES(t)]T
惩罚权重: wf = 10, 4wv = 10, 3wtrack = 10, 2R= diag(10, 1.0, 0.1)
+ 1) − εres(t + Np) ≥ Vref(t + 1) − εslack

Lyapunov稳定性证明:通过利用最优代价函数值)) 作为Lyapunov候选函数,引入松弛变量ε*(和正定二次状态惩罚项,确保Δ()) < 0,从而证明了所有运行状态轨迹的递归可行性和渐近稳定性。)) = J*(x(t)) as a Lyapunov candidate function, the inclusion of slack variable εslack and positive-definite quadratic state penalties guarantees ΔV*(x(t)) < 0, proving recursive feasibility and asymptotic stability across all operational state trajectories.
图5展示了所提出的动态终端库容耦合约束在双尺度分层能量管理系统(EMS)框架下,对水库水位跟踪性能的关键影响。如红色轨迹所示,传统的非耦合实时模型预测控制(MPC)模型无法维持长期水量平衡,导致在终端边界(T)处出现严重的调度失效和不可行性。相比之下,所提出的耦合MPC框架(蓝色轨迹)成功实现了终端时域平衡,确保实时运行轨迹严格遵循日前宏观调度参考值(黑色虚线),同时保障了市政用水的可用性。
Electronics 15 05037 g005
图5。双尺度分层EMS框架下动态终端容量耦合约束对油藏水位跟踪的影响。
图6展示了在第二阶段频率解耦MPC分配器框架内实现的快速频率解耦滤波器架构和功率分配路径。如图所示,原始功率失配信号(定义为ΔP grid (t) = P PV (t) − P load (t))被送入双速率滤波回路,以根据频率分离瞬态分量。高通滤波器路径(G HPF (s) = (τ f s)/(τ f s + 1),处理频率> 01 Hz的高频波动)通过亚周期动力学回路(Δt = 100 ms)将快速功率峰值路由至FESS快速执行器,从而产生瞬时有功功率输出P FESS (t)。
Electronics 15 05037 g006
图 6。高速频率解耦滤波器架构与功率分配路径。
相反,低通滤波器路径(G LPF (s) = 1/(τ f s + 1),处理 < 0.1 Hz 的低频分量)将缓慢变化的整体能量引导至 BESS SOC 管理器,以实现平滑的动态调节(产生 P BESS (t)),并引导至 VS-PHES 调速器,后者严格执行水锤约束(T w = 2.15 s)来调节整体水电出力 P hydro (t)。这种多时间尺度解耦策略有效防止了重型流体机械的高频机械疲劳,同时最大限度地延长了电化学储能资产的使用寿命。

4. 模拟设置与系统参数

所提出的框架被映射到一个经过修改的22.9 kV径向IEEE 33节点配电系统中。整个网络中安装的分布式光伏总容量为20 MW。使用MATLAB R2023b (MathWorks)/Simulink Simscape Electrical与CasADi集成,并采用Δ= 1 ms的电力系统微分方程求解器积分步长,进行了动态联合仿真。sim = 1 ms.
在本研究中,“系统频率”指在高IBR弱网条件下运行的局部22.9 kV配电馈线中心惯性(COI)频率。为了保守评估暂态稳定性极限,外部电网由一个具有故意设定的低惯性常数(H = 1.2 s)和一次调频特性(R = 5%)的等效同步发电机模型表示。如动态频率响应曲线所示,频率在5兆瓦光伏脱网后稳定在57.90 Hz的稳态偏差;这种未恢复的频率偏差是初级调频控制在没有二次调频(AGC)干预下作用的物理表现。
图7展示了在光伏(PV)掉落事故后,包含调速器下垂特性和低惯量转子动力学的等效电网频率响应模型的区块图。如闭环传递函数所示,光伏掉落引起的突发有功功率不平衡(ΔPL)直接扰动转子和负载动力学区块(1/(2Heqs + D)),导致频率偏差(Δω)。为了抵消这种不平衡,反馈回路通过一次下垂控制(1/R)、调速器执行器(1/(Tgs + 1))和汽轮机动力学(1/(Tts + 1))区块处理频率误差,以生成补偿性的机械功率注入(ΔPm),从而评估低惯量弱电网条件下系统的暂态频率稳定性。
Electronics 15 05037 g007
图7。包含调速器下垂特性和光伏掉落事故后低惯量转子动力学的等效电网频率响应模型框图。
为了阐明所提出的双速率控制架构的时间和功能层级,表2详细阐述了两级分层能量管理系统的运行规范。第一阶段负责日前能量调度,以1小时(ΔT = 1 h)的宏观时间步长运行,采用混合整数线性规划(MILP)优化宏观经济调度,降低电网购电成本,并缓解电池退化。相比之下,第二阶段的实时瞬态控制被划分为两个高速子回路。监督瞬态控制(第二阶段a)采用100毫秒分辨率的非线性模型预测控制(NMPC)算法,管理微瞬态动态设定点分配并执行终端约束。最后,本地执行器回路(第二阶段b)以10毫秒时间步长运行,利用一次下垂控制和矢量控制,提供亚周期(<20毫秒)动能功率注入,以实现超快速频率稳定。
表2。两阶段运行规范(修订后的双速率架构)。

为确保动态模拟的物理有效性和计算可重复性,表3详细列出了所提系统的全面技术参数配置矩阵。该表规定了升级后的变速抽水蓄能(VS-PHES)的水力机械约束,包括2.0 MVA的额定容量、12.5米的有效水头、2.15秒的水柱启动时间以及0.50秒的调速器执行器时间常数。此外,它还概述了混合储能资产的规格,为锂离子电池储能系统(BESS)设定了3.0 MWh的额定能量容量和92.0%的双向效率。机械飞轮储能系统(FESS)的参数设定为45.8 kg-m2转子转动惯量,最大运行速度为16,000 RPM,有效运行速度范围为370至430 rad/s。最后,数值模拟环境受到1.0 ms微分方程求解器步长的严格约束,以准确捕捉亚周期电学和流体动力学。
表3。用于系统模拟与可复现性的综合技术参数配置矩阵。

5. 案例研究与对比性能分析

5.1. 机器学习模型性能与执行延迟评估

为了增强NMPC回路内的实时扰动分类,研究人员评估了经Symlet-4离散小波变换(DWT)预滤波的1D卷积神经网络(1D-CNN)特征提取器与LSTM和SVM模型的性能。如表4所示,虽然在噪声环境下(SNR = 20 dB),1D-CNN比LSTM的准确率提升了1.45%,但其核心优势在于嵌入式TI TMS320F28379D DSP硬件上仅1.45毫秒的极低推理延迟,使其能够完美满足低于10毫秒的实时控制时限要求。
表4。机器学习模型性能、准确率与DSP硬件延迟矩阵。

5.2. 韩国22.9 kV 134节点城市网络实际案例验证

为证明其超越IEEE 33节点基准的可扩展性,该框架被部署于一个装机光伏容量为45 MW的韩国22.9 kV 134节点实际城市配电网络中。在商业用电高峰期发生12 MW光伏复合损失的故障情况下,所提出的双速率NMPC无缝协调了3.5 MW的升级混合设施。馈线电压跌落被控制在0.951 p.u.(相比未补偿时的0.912 p.u.),频率最低点维持在59.38 Hz,证实了其稳健的公用事业级性能。

5.3. 控制硬件在环(CHIL)与实验室原型路线图

控制硬件在环(CHIL)测试通过将运行第二阶段代码的双核TI C2000 DSP与RTDS实时电网模拟器进行接口连接来执行。CHIL结果证实,CAN总线通信延迟(2毫秒)和ADC量化噪声引起的性能下降可忽略不计(电压偏差<0.4%)。作为未来的研究路线图,计划在松实大学建立一个10 kW/30 kWh的物理实验原型,以验证硬件约束下的电液调速器响应。
为了严格基准测试系统性能,评估了四个对比基准框架:
  • 基准A(无补偿电网):馈线在零辅助储能支持下运行。
  • 基准B(传统仅BESS):电网仅由独立的BESS支持,不具备FESS频率解耦或VS-PHES能力。
  • 基准C(基于规则/PI拆分的混合硬件):共址的VS-PHES、BESS和FESS,由传统的PI频率分频滤波器控制,不具备NMPC预测能力。
  • 基准D(提议框架):采用提议的两阶段双速率NMPC(带动态终端容量耦合)来协调VS-PHES、BESS和FESS。

图8展示了使用多元高斯分布对分布式光伏(PV)下调掉落进行经验表征,捕捉了城市电网馈线集群的空间-时间相关性和联合概率密度分布。具体而言,左上角的子图展示了不同光伏发电场址(例如相隔2公里)之间强烈的空间相关性,而右上角的子图则描绘了用于建模输出不确定性的校对边缘概率密度函数(α = 6.0, β = 1.5)。底部的子图展示了合成的空间相关光伏事故场景(S raw),突出了高峰时段同步发生的严重掉落事件,该事件会触发局部电压和频率暂态。
Electronics 15 05037 g008
图8。使用多元高斯分布对分布式光伏降功率掉落事件进行经验表征,展示了城市电网馈线集群的时空相关性及联合概率密度分布。
馈线事故物理故障(5.0 MW光伏掉落事件):
5.0 MW 总有功功率损失(本地馈线光伏掉落事件)
Electronics 15 05037 i001
图9展示了在发生5.0兆瓦光伏掉落后,动态频率和局部电压分布的曲线= 2.0秒。基准A(未补偿)遭受严重的频率掉落,低至57.90 Hz,电压跌至0.808 p.u.。基准C(PI频率解耦)部分缓解了扰动,但由于水力机械闸门延迟,表现出暂态超调和振荡。相比之下,所提出的双速率NMPC(基准D)利用快速FESS动能注入(<20毫秒)来遏制初始RoCoF骤降,在1.5秒内将频率恢复至59.35 Hz,终端电压恢复至0.952 p.u.。 = 2.0 s. Baseline A (uncompensated) suffers severe frequency drops down to 57.90 Hz and voltage sags to 0.808 p.u. Baseline C (PI frequency-decoupling) mitigates the disturbance partially but exhibits transient overshoot and hunting due to hydro-mechanical gate latency. In contrast, the proposed dual-rate NMPC (Baseline D) leverages fast FESS kinetic injection (<20 ms) to arrest the initial RoCoF slump, restoring frequency to 59.35 Hz and terminal voltage to 0.952 p.u. within 1.5 s.
Electronics 15 05037 g009
图9。不同控制基准下的瞬态频率响应曲线和局部电压分布。
如图图10所示,在关键馈线末端(Bus 19),基准B(仅BESS)突破了法定电压限制(0.95 p.u.),由于变流器电流限制,电压跌至0.931 p.u.。基准C缓慢恢复,但经历了长时间的低电压压力。基准D(提出的双尺度框架)成功将电压维持在0.952 p.u.的法定边界,防止了低电压负荷切除。
Electronics 15 05037 g010
图10。在不同控制策略下,局部5.0 MW分布式光伏掉落故障时关键终端节点(节点19)的暂态电压动态响应。

6. 生命周期成本(LCC)与辅助服务市场定价评估

为了量化长期经济回报,制定了20年生命周期成本(LCC)矩阵。为消除监管合规障碍,我们设计了市场驱动的动态节点电价。Pr,total(i, t在提议的动态节点电价下,归属于RE100私人投资者的补偿定价制定如下:
) = SMP(r,total() + α·Δ, ) + β·[))] + γ·t) + α·ΔPLoss(i, t) + β·[Zi-slack/(Vi(t)·Rpen(t))] + γ·CInertia(i, t)

为了验证所提出的多性向储能架构的长期经济和财务可行性,针对三种不同的部署场景(终端、均匀和负荷中心)进行了为期20年的生命周期成本(LCC)分析和辅助服务市场结算评估,综合结果展示在表5.
表5。不同布局下的综合生命周期成本(LCC)矩阵和辅助服务市场结算评估。

根据分析,在系统运行波动极大的情景1(终端节点集中)下,20年累计电池单元退化成本增加至42万美元。然而,飞轮储能系统(FESS)提供的动态动能屏蔽带来了191万美元的BESS更换节省费用。此外,它还产生了128万美元的辅助服务市场结算收益,清楚地证实了所提出的多能互补架构的长期经济优势。
在情景1(终端节点集中)下,高运行波动使20年累计单元退化成本增加至42万美元。然而,FESS提供的动态动能屏蔽带来了191万美元的BESS更换节省费用,并产生了128万美元的辅助服务市场结算收益,证实了多能互补架构的长期经济优势。

7. 结论

本研究提出了一种数学严谨且土木工程可行的结构升级框架,将利用率不足的城市防洪基础设施转化为去中心化的2.0 MVA低水头变速抽水蓄能(VS-PHES)枢纽。通过在单一市政流域设施内共址部署高速机械飞轮储能系统(FESSs)和高能量密度锂离子电池储能系统(BESSs),所提出的架构建立了一种高度韧性的多能流储能配置。整个框架采用一种新型两阶段分层能源管理系统(HEMS)进行动态管理,该系统由高速频率解耦滤波器拓扑结构增强。
该系统由采用双速率非线性模型预测控制(NMPC)控制器(100毫秒监督/10毫秒局部回路)且具有动态终端电压耦合的两阶段分层能源管理系统(EMS)控制,在维持第一阶段宏观经济能源计划的同时,提供低于20毫秒的一次调频和电压调节。在发生局部5兆瓦光伏掉落事故的改进型IEEE 33节点配电网上的验证表明,3.5兆瓦混合资产将终端电压从0.931标幺值恢复至0.952标幺值,并将频率最低点控制在59.35赫兹,性能优于传统基于规则的控制器42%。通过标准IEEE 33节点和实际134节点韩国城市网络(辅以控制硬件在环测试)的验证,证实了其高度的可扩展性。此外,飞轮储能系统(FESS)屏蔽功能在20年内可节省191万美元的电池储能系统(BESS)更换费用。集成自动化市政防洪优先级状态切换逻辑