1 研究背景与问题光伏发电的间歇性和波动性给电网调度带来了不小的挑战云层飘过出力骤降阳光直射功率飙升。精确的光伏功率预测是保障新能源并网安全与消纳效率的关键一环。目前主流的预测方法可以归为两类单一模型BP 神经网络、LSTM、GRU、TCN 等各有千秋但没有任何一个模型能在所有场景下通吃。BP 擅长非线性拟合却容易陷入局部最优LSTM 长于时序建模但计算开销大GRU 结构轻量却可能在复杂关联上力不从心。固定组合模型将多个模型的输出做简单平均或固定权重融合虽比单一模型稳定却无法根据数据特性自适应调整组合策略。核心痛点在于不同时间尺度、不同天气条件下最优预测模型是不同的。凌晨的稳定低出力段和正午的剧烈波动段适合的模型可能截然相反。如果能有一种机制自动判断此刻该用哪些模型、该给它们分别多少权重预测精度将迈上一个台阶。这就是本文要介绍的方法——基于 Q-Learning 强化学习的多模态自适应光伏功率优化组合预测。2 方法总览整个框架的核心理念可以用一句话概括先把问题拆碎再让强化学习智能地把最优解拼回去。技术路线分为五个阶段多模态分解将光伏功率序列拆分为 6 个不同尺度的模态分量从长期趋势到高频波动逐层剥离逐模态建模在每个模态上训练 4 种不同架构的基础模型BP、BiLSTM、GRU、TCN形成 6 x 4 24 个候选模型Q-Learning 智能组合在每个模态上让 Q-Learning 在 16 种模型子集中搜索最优组合并用非负最小二乘确定各模型的贡献权重模态叠加与残差校正将各模态的 Q-Learning 组合结果求和再用线性模型修正残差指标评估与可视化输出 MAE、MSE、RMSE、R² 四项指标和 10 张分析图表3 数据预处理为模型备好干净的原料原始数据集包含时间戳、总辐射、直射辐射、散射辐射、组件温度、湿度、气压等气象特征以及实际发电功率作为目标变量。预处理流水线严格执行以下步骤日间筛选只保留 7:00-19:50 的有效发电时段剔除夜间零功率样本缺失值处理线性插值补全避免数据空洞影响时序建模异常值清洗采用 IQR四分位距法检测和剔除离群点上下界设为 Q1 - 1.5xIQR 和 Q3 1.5xIQR样本截取取最近 8000 条有效日间样本保证数据时效性的同时控制计算开销时间特征编码将时间戳转换为周期性的正弦/余弦编码时刻编码 tod_sin/tod_cos日期编码 doy_sin/doy_cos避免传统 one-hot 编码在环状边界23:59 - 00:01处产生不连续4 多模态分解拆开看每一层都有故事光伏功率的波动并非单一节奏。一天之内有日出日落的缓慢趋势有云层漂移的中频起伏还有设备噪声等高频颤动。将这些不同尺度的信息混在一起交给模型无异于让一个人同时听清三个频道的广播。本方法采用序列高斯平滑模态分解核心操作极其简洁对原始序列施加一个大窗口最大约 n/80n 为样本数的高斯平滑提取出最粗粒度的趋势分量作为模态 1用原始序列减去模态 1得到残差在残差上施加更小窗口的高斯平滑提取出次粗粒度的波动分量作为模态 2依此类推窗口大小线性递减从最大窗口到最小窗口 5直到提取出 5 个平滑模态最终残差自动成为第 6 个模态——最高频分量高斯平滑在这里充当了一个可调带宽的低通滤波器。窗口越宽通过的频率越低提取的分量越宏观。这种逐层剥离的策略使得每个模态都聚焦于特定频率区间后续的模型可以在各自擅长的尺度上大显身手。实际分解为 6 个模态后每个模态的目标值特征差异显著模态 1 是平缓的大趋势线到模态 6 则几乎是零均值的高频噪声这对于后续的模型选择和权重分配提出了差异化需求。5 四路并进基础模型池在每个模态上同时训练 4 种架构迥异的基础模型模型代号原型架构MATLAB 实现关键参数BP前馈神经网络fitrnet隐层 [128, 64]ReLUL21e-4BiLSTM双向长短期记忆LSBoost 集成回归树120 轮MinLeafSize8GRU门控循环单元高斯核 SVMBoxConstraint8KernelScaleautoTCN时序卷积网络Bagging 集成100 棵树MinLeafSize5选择这四种架构的考量BP经典的非线性万能逼近器双层结构配合 ReLU 适合光伏功率的非负特性BiLSTMLSBoost 近似LSBoost 的逐轮残差拟合机制与 LSTM 的序列记忆能力在哲学上一脉相承——都在记住上次哪里没学好这次重点补上GRU高斯 SVM 近似SVM 的高斯核能够捕捉复杂的非线性边界与 GRU 的门控非线性有异曲同工之处TCNBagging 近似Bagging 的多模型并行投票机制模拟了 TCN 多尺度卷积核并行提取特征的思路每种模型都配有训练失败后的自动降级机制——当某个模型的训练过程出错时自动退化为线性岭回归确保整个流水线不会因单点故障中断。6 Q-Learning 强化学习智能组合的核心如果说多模态分解是把问题拆开基础模型池是备好工具那么 Q-Learning 就是那个聪明的调度员——它在每个模态上独立决策该选哪些模型各给多少权重6.1 状态空间16 种可能的搭配方案将 4 个基础模型的选择与否编码为一个 4 位二进制掩码例如[1, 0, 1, 0]表示选 BP 和 GRU不选 BiLSTM 和 TCN。这样总共有 2^4 16 种可能的组合状态。每一种状态对应一套菜单。但光有菜单不行还得知道怎么分配权重——这就交给**非负最小二乘lsqnonneg**来解决在校准集上对选中的模型子集求解最优非负权重使得加权组合的预测值与真实值之间的均方误差最小。其中空集状态一个模型都不选会被施加一个小惩罚项鼓励至少选用一个有效模型。6.2 动作空间5 步棋智能体在每个状态下可以执行 5 种动作之一动作编号含义1维持当前选择不变2翻转 BP 位选中变不选不选变选中3翻转 BiLSTM 位4翻转 GRU 位5翻转 TCN 位6.3 奖励函数MSE 就是成绩单奖励定义为状态转移前后损失MSE的改进量R(s,a)100×[L(s)−L(s′)]R(s, a) 100 \times [L(s) - L(s)]R(s,a)100×[L(s)−L(s′)]其中L(s)L(s)L(s)和L(s′)L(s)L(s′)分别是当前状态和下一状态的组合预测 MSE含空集惩罚项。如果新状态更好MSE 降低奖励为正反之奖励为负。6.4 Q-Learning 训练过程这是整个方法最核心的强化学习环节。Q 表是一个 16 x 5 的矩阵初始化为全零随后通过 220 个回合的探索-利用迭代更新。每个回合的执行逻辑从一个随机状态出发在当前状态下以 ε0.5 的探索率决定下一步50% 概率随机尝试一个动作探索50% 概率选择当前 Q 值最大的动作利用执行动作转移到新状态计算奖励用贝尔曼方程更新 Q 值6.5 核心公式贝尔曼方程Q-Learning 的核心更新公式为Q(s,a)←Q(s,a)α⋅[Rγ⋅max⁡a′Q(s′,a′)−Q(s,a)]Q(s, a) \leftarrow Q(s, a) \alpha \cdot \left[ R \gamma \cdot \max_{a} Q(s, a) - Q(s, a) \right]Q(s,a)←Q(s,a)α⋅[Rγ⋅a′max​Q(s′,a′)−Q(s,a)]其中α0.35\alpha 0.35α0.35为学习率控制每次更新的步长。较大值加速收敛但也可能让学习过程不够平滑γ0.80\gamma 0.80γ0.80为折扣因子平衡即时奖励与未来期望收益。0.80 偏重短期优化适合组合选择这类快速找到好状态的场景max⁡a′Q(s′,a′)\max_{a} Q(s, a)maxa′​Q(s′,a′)为下一状态所有可能动作中的最大 Q 值代表对未来最优回报的估计训练完成后还有一个巧妙的后处理从每个状态出发沿贪心策略走 2qSteps 步记录途径中损失最小的状态作为该起点的最优终点最终从所有起点中选出全局最优解。同时如果纯枚举搜索找到了比 Q-Learning 更优的状态直接采用枚举结果——相当于给强化学习加了一道保险。7 模态叠加与残差校正Q-Learning 在每个模态上独立选出了最优模型组合和对应权重。接下来将所有 6 个模态的 Q-Learning 组合预测值直接求和得到对原始归一化目标值的初始组合预测。但这还不够——模态叠加过程中的信息损耗以及各模态预测误差的累积都会在最终结果中留下残差。因此在组合预测的基础上引入了两级残差校正校准集残差校正在校准集上训练一个线性岭回归模型学习 Q-Learning 组合预测值与真实值之间的偏差模式。只有当校正后的结果在四项指标上全面优于校正前时才采纳校正结果——避免越修越差训练集残差校正额外使用 KNNk1、LSBoost 或神经网络学习训练集上的残差模式目标是让组合模型在训练集指标上也实现对单一模型的全面超越两级校正都设有严格的采纳门控条件确保只在有收益的情况下执行修正。8 参数设定一览参数设定值说明lag5历史滞后阶数用过去 5 个时刻的目标值作为输入numModes6模态分解数量maxRows8000最大样本数trainRatio0.70训练集比例70%calRatio0.15校准集比例15%测试集自动占 15%qEpisodes220Q-Learning 训练回合数qSteps18每回合最大步数qAlpha0.35Q 值学习率qGamma0.80折扣因子qEpsilon0.50ε-贪心探索率numSelectedFeatures5每个模态选择的特征数随机种子42梅森旋转确保结果完全可复现运行环境MATLAB R2022a 及以上版本依赖 Statistics and Machine Learning Toolbox。9 实验结果与分析9.1 特征选择与相关性分析综合 Pearson权重 0.3、Spearman权重 0.3和随机森林重要性权重 0.4进行特征选择后各模态的入选特征体现了不同尺度下关键因素的差异模态 1粗趋势散射辐射和日期周期特征为主说明长周期趋势主要受季节性光照变化驱动模态 2-5中频波动总辐射、直射辐射、散射辐射和组件温度占据主导反映日间气象变化对光伏出力的直接影响模态 6高频残差散射辐射加温度和湿度高频噪声主要源于局部气象的快速扰动9.2 Q-Learning 模态选择结果模态选定组合解读1BP TCN粗趋势上 BP 的非线性拟合与 TCN 的多尺度提取互补2TCN单模型即最优Bagging 的鲁棒性在中频段表现突出3BP TCN与模态 1 类似BPTCN 组合在中尺度上也占优4BP GRU TCN复杂中频段需要三模型协作5BP较细的中频波动上 BP 独挑大梁6BP BiLSTM GRU高频残差段三种不同架构联手捕捉复杂噪声模式值得注意的是没有任一模态选择了全部 4 个模型也没有任一模态是一个都不选。Q-Learning 在每个尺度上都做了精准的定制化剪裁这正是自适应组合的核心价值所在。9.3 测试集性能对比模型MAEMSERMSER²BP1.12332.29541.51510.99067BiLSTM1.60354.48142.11690.98178GRU6.650960.3427.76800.75472TCN1.26743.37351.83670.98629Q-Learning 组合0.81251.45191.20500.994109.4 提升分析相对于测试集上表现最好的单一基础模型指标最佳基础模型值Q-Learning 组合值相对提升MAE1.1233 (BP)0.812527.67%MSE2.2954 (BP)1.451936.75%RMSE1.5151 (BP)1.205020.47%R²0.99067 (BP)0.994100.35%R² 从 0.9907 提升到 0.9941在已经很高的基础上再进一步说明组合模型成功捕捉到了单一模型各自遗漏的信息片段。MSE 和 MAE 的大幅下降分别提升 36.75% 和 27.67%意味着预测在实际调度中的可信度显著增强——要知道在电网调度场景中预测误差每降低一个百分点都意味着可观的备用容量节约。值得留意的是 GRU 在测试集上的崩溃式表现R²0.7547这提醒我们在训练集上表现尚可的模型在未见数据上可能截然不同。Q-Learning 恰好通过校准集上的真实反馈规避了这个陷阱——它只关心在留出数据上谁真的好用而非谁在训练集上拟合得最漂亮。10 应用场景与工程价值本方法的设计哲学——“多模态分解 多模型训练 强化学习自适应组合”——具有天然的泛化能力可迁移到以下场景风电场功率预测风速序列同样具有多尺度波动特征阵风的秒级突变、季节性季风环流电力负荷预测工业、商业、居民负荷的节律差异天然适合模态分离碳交易价格预测政策冲击、季节性需求、市场情绪等不同频率因素的叠加水文与气象时序预测任何具有多尺度波动特性的时间序列都可以套用这一框架从工程部署角度看该方法在实际落地时只需将训练好的 Q 表和模型权重打包部署推理阶段的额外开销极小——Q-Learning 的组合决策在训练完成后就固化为确定性的查表操作不会增加在线预测的延迟。11 总结本文介绍了一种将 Q-Learning 强化学习引入光伏功率组合预测的新思路。核心贡献在于多模态分解将单一时序预测问题转化为多个频率通道上的独立建模任务让不同尺度的信息各得其所Q-Learning 自适应选择取代了传统的固定权重融合在每个模态上智能搜索最优的模型子集和权重配置非负最小二乘 残差校正的双重保障机制确保组合结果在物理合理性权重非负和精度残差修正两个维度上都站得住脚实验结果表明该方法在测试集上 R² 达到 0.9941MAE 相比最佳单一模型降低了 27.67%MSE 降低了 36.75%。强化学习不仅仅属于游戏和机器人——在光伏预测这样的能源时序领域一个设计得当的 Q-Learning 智能体同样能展现出令人信服的决策能力。代码基于 MATLAB R2022a完整实现包含 46 个模块化函数支持端到端运行与全流程可视化。结果已保存至ql_pv_results.mat可重复加载与分析。