风电光伏与储能互补调度优化及Matlab实现

1. 风电、光伏与储能互补调度研究背景

可再生能源的快速发展正在重塑全球能源格局，其中风电和光伏发电作为两大主力清洁能源，装机容量呈现爆发式增长。然而，这些能源的间歇性和波动性给电网调度带来了巨大挑战。以中国西北某省为例，2022年风电和光伏的弃电率曾高达15%，主要原因就是缺乏有效的调节手段。

储能技术的引入为解决这一问题提供了关键突破口。特别是电池储能和废弃矿井抽水蓄能（UPSH）的组合，能够实现分钟级到周级的多时间尺度能量调节。这种"快慢结合"的储能方式，既能够平抑风光出力的快速波动，又可以实现能量的季节性转移。

2. 系统建模与优化框架

2.1 目标函数设计

在Matlab建模中，我们构建了一个双层优化目标函数。上层以系统运行经济性最大化为目标，包含以下关键项：

matlab复制f = sum(
    C.*P_w +          % 风电收益
    C.*P_pv +         % 光伏收益  
    C.*P_h +          % 抽蓄发电收益
    C.*P_d -          % 电池放电收益
    C_p.*P_p -        % 抽蓄抽水成本
    C_c.*P_c -        % 电池充电成本
    yibuxil_lack*P_lack -  % 缺电惩罚
    yibuxil_DL*P_DL -      % 弃电惩罚
    k_ps_h*P_h -      % 抽蓄运行成本
    k_ba_d*P_d        % 电池运行成本
);

这个目标函数充分考虑了各种收益和成本项，通过经济杠杆引导系统优化运行。

2.2 约束条件实现

系统约束通过YALMIP工具箱的SDPSET函数实现，主要包括：

matlab复制F = [];
for t = 1:n
    % 功率平衡约束
    F = [F, P_w(t) + P_pv(t) + P_h(t) + P_d(t) + P_lack(t) == P_load(t)];
    
    % 抽蓄运行约束
    F = [F, P_p(t)*P_h(t) == 0];  % 抽水和发电不同时
    F = [F, E(t+1) == E(t) + eta_p*P