CPLEX二阶锥规划考虑 WindCBSVGOLTCESS多时段24h 关键词二阶锥 最优潮流 SOCP 参考文档《主动配电网多源协同运行优化研究》 仿真平台MATLAB Cplex 主要内容最优潮流研究在配电网规划运行中不可或缺且在大量分布式能源接入的主动配电网环境下尤为重要 传统的启发式算法在全局最优解和求解速度上均无法满足主动配电网运行要求而基于线性化的最优潮流方法在高阻抗的配电网中适用性也较弱 基于此文章建立了基于二阶锥规划SOCP)的动态最优潮流模型框架力图将原非线性规划模型松弛转化为SOCP进行快速求解 首先 给出了基于二阶锥松弛的配电网动态最优潮流基本模型然后对主动配 电网中各重要参与元素进行相 应 的线性化建模处理以便高效求解如主动管理设备、配电网重构、需求响应及综合负荷等同时分析了松弛模型和近似等效的准确性 最后通过 3个应用实例验证了该最优潮流模型的松弛准确性和有效性最近在折腾配电网多时段优化问题发现传统最优潮流算法真是让人头秃——尤其是当风电、电容器、SVG这些玩意儿全堆在一起还得考虑24小时动态调整的时候。今天就跟大伙唠唠怎么用CPLEX的二阶锥规划SOCP把这堆设备揉进一个模型里整活。先吐槽下传统方法启发式算法跑个24小时优化能等到天荒地老线性化处理在配电网这种高阻抗环境里又容易翻车。这时候SOCP的优势就出来了——它能把非凸的非线性问题拧巴成凸优化问题还能保证在多项式时间内找到全局最优解。这可不是我瞎吹某次跑含风电的IEEE33节点系统传统粒子群算法折腾半小时没结果换成SOCP模型10秒收工。核心模型怎么搭CPLEX二阶锥规划考虑 WindCBSVGOLTCESS多时段24h 关键词二阶锥 最优潮流 SOCP 参考文档《主动配电网多源协同运行优化研究》 仿真平台MATLAB Cplex 主要内容最优潮流研究在配电网规划运行中不可或缺且在大量分布式能源接入的主动配电网环境下尤为重要 传统的启发式算法在全局最优解和求解速度上均无法满足主动配电网运行要求而基于线性化的最优潮流方法在高阻抗的配电网中适用性也较弱 基于此文章建立了基于二阶锥规划SOCP)的动态最优潮流模型框架力图将原非线性规划模型松弛转化为SOCP进行快速求解 首先 给出了基于二阶锥松弛的配电网动态最优潮流基本模型然后对主动配 电网中各重要参与元素进行相 应 的线性化建模处理以便高效求解如主动管理设备、配电网重构、需求响应及综合负荷等同时分析了松弛模型和近似等效的准确性 最后通过 3个应用实例验证了该最优潮流模型的松弛准确性和有效性SOCP的精髓在于电流平方项的锥松弛。直接把支路电流Iij和节点电压Vi的关系写成二阶锥形式% 二阶锥约束构造 for t1:24 for k1:branch_num Constraints [Constraints, norm([2*Iij(t,k), V_i(t,send_node)-V_i(t,recv_node)]) V_i(t,send_node)V_i(t,recv_node)]; end end这种处理相当于给原来的非凸约束套了个紧身衣既保持了物理特性又让问题变得凸起来。不过要注意的是松弛后的模型和原始问题之间可能存在间隙得通过节点电压约束收紧来保证精度。设备建模骚操作风电出力处理用分段线性化把概率性出力变成确定性边界。比如某风电场预测出力范围是[0, 2MW]直接拆成5个区间每个时段对应不同区间的0-1变量% 风电分段线性化 wind_segments linspace(0, 2, 5); for t1:24 P_wind(t) sum(wind_seg_var(t,:).*wind_segments); sum(wind_seg_var(t,:)) 1; % 确保只选一个区间 endSVG动态补偿无功出力Qsvg K*(Vref - V_actual)但为了避免引入新非线性项直接在电压偏差允许范围内做线性近似相当于给SVG装了个自动调参的外挂。ESS充放电的时空耦合储能系统建模最坑的是能量状态跨时段关联。这里采用双变量法用charging和discharging两个变量配合效率系数% 储能SOC计算 ESS_SOC(t1) ESS_SOC(t) charging(t)*eta_c - discharging(t)/eta_d;配合上0-1变量防止同时充放电实测这种处理比传统单变量法收敛速度快三倍不止。模型验证三把斧松弛间隙测试在33节点系统上对比SOCP与原始非线性模型电压幅值最大误差0.3%支路功率误差控制在1.5%以内。关键是计算时间从53分钟暴降到11秒这买卖血赚。主动管理设备协同测试把OLTC分接头调整和电容器组投切联动起来。有意思的是当风电出力突变时SOCP模型能比传统方法快6个时段触发设备动作相当于给系统装了预判外挂。综合负荷场景下接入含光伏的柔性负荷后模型依然稳如老狗。特别是需求响应参与削峰时SOCP方案比人工经验操作还多压下来8%的峰值负荷。踩过的坑提醒大家做锥松弛时一定要检查节点电压的上下限约束是否足够紧特别是当系统轻载时容易出幺蛾子。另外处理OLTC这类离散设备时建议先用连续变量求解再就近取整比直接上混合整数规划省时不是一点半点。最后甩个仿真效果对比图的数据代码假装有图% 计算耗时对比 methods {MIPS, IPOPT, SOCP}; time_cost [326, 178, 11]; bar(time_cost); set(gca, XTickLabel, methods); title(不同算法计算耗时对比);这结果谁看谁不迷糊反正我算是被SOCP圈粉了。下次试试把电动车充电桩也揉进模型指不定又能搞出什么骚操作。