二阶锥松弛化的无功优化代码 考虑储能系统ESS，有载调压分接头OLTC，离散电容器CB和静止无功补偿SVC设备约束 以购电和网损成本最小为目标 采用IEEE 33节点算例仿真 使用二阶锥松弛将模型线性化 使用yalmip-cplex/gurobi求解器编写

一、代码概述

本代码基于二阶锥松弛化技术，构建了考虑储能系统（ESS）、有载调压分接头（OLTC）、离散电容器（CB）及静止无功补偿器（SVC）等多元设备约束的配电网24小时无功优化模型。通过多目标优化（购电成本最小化、网损最小化、电压偏差最小化），实现配电网运行状态的精细化调控，同时满足潮流平衡、设备运行边界、电压质量等约束条件，适用于含新能源（风电、光伏）接入的现代配电网调度决策场景。

二、核心功能模块

（一）系统参数初始化模块

该模块负责完成配电网基础参数、设备参数及时间维度参数的标准化配置，为后续优化计算提供数据支撑，所有参数均基于标幺值体系进行统一转换，确保计算精度与工程实用性。

1. 基础参数配置
   - 时间维度：设定优化周期为24小时（NT=24），覆盖一天内负荷与新能源出力的典型波动周期。
   - 基准值体系：以10MW为基准功率（SB）、12.66kV为基准电压（VB），推导基准阻抗（ZB=VB²/SB）与基准电流（IB=SB/(√3×VB)），作为所有电气量标幺化的基础。
   - 配电网拓扑：定义32条线路（NL=32）、33个节点（NB=NL+1）的拓扑关系，包含线路首末节点编号、负荷数据（有功Pload、无功Qload）、阻抗参数（电阻R、电抗X）及线路潮流限额（Pmax、Pmin），并通过标幺化转换为计算用单位。
2. 负荷与新能源参数配置
   - 负荷特性：基于各节点负荷占系统总负荷的比例（Pdlevel、Qdlevel），结合24小时系统总负荷曲线（Pdday、Qdday），生成各节点24小时有功/无功负荷时序数据（Pd、Qd），并完成标幺化处理。
   - 风电参数：在节点3接入风电，输入24小时风电出力预测曲线（Wg1），设定风电无功出力范围（QgWTmin、QgWTmax），明确有功出力上下限（WTPgmin、WTPgmax）。
   - 光伏参数：在节点6接入光伏，输入24小时光伏出力预测曲线（PVg1），考虑光伏通常不提供无功出力的特性，仅设定有功出力上下限（PVPgmin、PVPgmax）。
3. 无功补偿与储能设备参数配置
   - OLTC：在16、23号线路配置OLTC，设定调节范围（0.95-1.05pu）、调节步长（0.01pu），共11个调节档位，形成时序化调节矩阵（T_OLTC）。
   - 离散电容器（CB）：在5、10、17、30号节点配置离散电容器，每台设备最大补偿容量0.2MVar、最小0容量、调节步长0.05MVar，采用3档位（0-1-2档）线性化处理策略。
   - SVC：在6、18号节点配置连续无功补偿器，设定最大补偿容量分别为2MVar、1MVar，最小补偿容量分别为-0.2MVar、-0.1MVar，支持无功功率的连续平滑调节。
   - 储能系统（ESS）：在7号节点接入储能，设定最大充电功率0.5MW、最大放电功率0.4MW，荷电状态（SoC）范围0-1MWh、初始SoC 0.1MWh，充放电效率分别为0.98、0.95，通过二进制变量实现充放电状态互斥控制。

（二）决策变量定义模块

围绕配电网运行状态与设备调控需求，定义多类型决策变量，涵盖功率流、电压、设备出力及状态变量，全面覆盖优化问题的控制维度：

1. 功率流变量：包含线路首端有功功率（P）、无功功率（Q），下游子线路首端功率（PkidLine、QkidLine），线路电流幅值平方（I2），覆盖潮流计算的核心电气量。
2. 节点电压变量：定义节点电压幅值平方（U2），作为电压质量约束的核心变量，避免平方根运算带来的非线性复杂度。
3. 新能源出力变量：风电有功/无功出力（PgWT、QgWT）、光伏有功出力（Pg_PV），匹配新能源发电的调控特性。
4. 无功补偿设备变量：SVC无功出力（QgSVC）、OLTC档位选择变量（eij、y_j）、离散电容器档位决策变量（xd、dlt），其中OLTC与离散电容器通过二进制变量实现离散状态的线性化表达。
5. 储能系统变量：充电功率（Pchar）、放电功率（Pdisc），充放电状态控制变量（Achar、Adisc），以及荷电状态（E_ES），确保储能运行的时序连续性与安全性。
6. 电网交互变量：上级电网有功/无功注入功率（Pgrid、Qgrid），作为购电成本计算与系统功率平衡的关键变量。

（三）约束条件构建模块

基于配电网运行物理规律与设备安全准则，构建多维度约束体系，确保优化结果的可行性与工程合规性，约束体系分为以下类别：

1. 储能系统约束（ConsES）
   - 状态互斥约束：通过Achar + Adisc ≤ 1确保同一时刻储能不发生充放电冲突。
   - 荷电状态时序约束：EES(t) = EES(t-1) + Pchar(t)×ηchar - Pdisc(t)/ηdisc，且初始与末时刻荷电状态相等（EES(1)=EES0、EES(24)=EES0），保证储能周期运行的一致性。
   - 边界约束：荷电状态在EESmin-EESmax范围内，充放电功率分别不超过pcharmax、pdisc_max，且非储能节点的充放电变量置零。
2. 潮流平衡约束（ConsP、ConsQ）
   - 有功平衡：线路首端有功功率扣除线路电阻损耗（R×I2），叠加新能源有功出力、储能充放电功率后，等于下游子线路功率与节点负荷之和，确保系统有功功率守恒。
   - 无功平衡：分两种场景构建约束——安装离散电容器的节点，需叠加电容器补偿容量；未安装节点则叠加SVC与风电无功出力，最终与下游子线路无功功率、节点无功负荷平衡，覆盖不同无功补偿配置的场景。
3. 电压约束（ConsU）
   - 根节点电压固定：设定节点1（上级电网接入端）电压为1.05²pu，作为电压参考基准。
   - 线路电压降落约束：不含OLTC的线路，基于阻抗参数与功率流计算电压降落（U2(Lj)=U2(Li)-2×(R×P+X×Q)+(R²+X²)×I2）；含OLTC的线路，通过档位变量（eij、y_j）线性化处理变比影响，确保电压降落计算与OLTC调节特性匹配。
   - 电压质量约束：所有节点电压幅值平方在0.95²-1.05²pu范围内，满足配电网电压质量标准。
4. 二阶锥约束（Cons_SOC）：通过I2×U2(Li) ≥ P² + Q²将潮流方程的非线性约束转化为二阶锥标准形式，避免传统潮流计算的非线性迭代问题，提升求解效率与收敛性。
5. 设备出力边界约束
   - 线路潮流：线路有功功率在Pmin-Pmax范围内，避免线路过载。
   - 新能源：风电、光伏有功出力不超过预测上限，风电无功出力在设定范围内。
   - 无功补偿：SVC无功出力在SVCmin-SVCmax范围内，离散电容器档位通过xd变量约束在0-(Cjmax-Cjmin)/Sj范围内。
   - 电网交互：上级电网注入功率（Pgrid、Qgrid）非负，符合配电网从上级电网受电的运行模式。

（四）目标函数构建模块

以配电网经济运行与运行质量为核心，构建多目标优化函数，通过加权整合实现综合效益最优：

1. 购电成本最小化（Objgrid）：基于24小时分时电价曲线（Eleprice），计算上级电网有功注入功率（Pgrid）的总购电成本，电价单位转换为元/(pu·h)，确保成本计算与标幺值体系一致。
2. 网损最小化（Obj_Ploss）：通过线路电流幅值平方（I2）与电阻（R）的乘积求和，计算24小时线路总有功损耗，实现配电网能效优化。
3. 电压偏差最小化（Obj_Vol）：以1.05²pu为理想电压基准，计算所有节点24小时电压幅值平方（U2）与基准值的偏差总和，提升系统电压质量的稳定性。
4. 目标函数整合：默认采用Obj = Objgrid + ObjPloss，优先保障经济运行与能效；可通过加权系数调整，兼顾电压偏差目标，适应不同调度需求。

（五）求解与结果可视化模块

1. 求解配置：选择Gurobi作为求解器，设定求解精度（MIPGapAbs=1e-3、MIPGap=1e-3），平衡求解速度与结果精度；支持切换至CPLEX、MOSEK等求解器，适应不同计算环境。
2. 求解结果判断：通过sol.problem == 0判断求解是否成功，成功则输出总运行成本（C_Opr），失败则提示错误信息并输出求解日志。
3. 结果可视化：生成多维度图表直观展示优化结果：
   - 有功出力源曲线：对比风电、光伏、上级电网24小时有功出力变化，清晰呈现新能源消纳与电网供电的协同关系。
   - 无功出力源曲线：展示风电、两台SVC及上级电网的无功出力时序特性，反映无功补偿设备的调控效果。
   - 节点电压三维图：通过mesh图呈现33个节点24小时电压幅值分布，直观评估电压质量空间与时序变化。
   - OLTC档位图：以阶梯图展示两台OLTC24小时档位调节轨迹，反映分接头的动态调压过程。
   - 离散电容器动作图：展示4台离散电容器的无功补偿容量时序变化，体现离散设备的档位切换逻辑。
   - 储能运行图：通过双Y轴图同步展示储能充放电功率（柱状图）与荷电状态（折线图），清晰呈现储能的充放电策略与状态演变。

三、代码核心价值与应用场景

1. 技术价值：采用二阶锥松弛化处理潮流非线性问题，结合离散设备（OLTC、CB）的线性化建模，实现混合整数二阶锥规划问题的高效求解，突破传统配电网优化中“非线性+离散变量”的求解瓶颈。
2. 应用场景：适用于含新能源接入、多元无功补偿设备配置的配电网日前调度，可支撑调度人员制定购电计划、设备调控策略，同时为配电网规划（如无功补偿设备选址定容）提供量化分析工具。
3. 扩展性：支持新增设备模型（如分布式电源、电动汽车充电桩）、优化目标（如碳排放成本），通过模块化设计可快速适配不同配电网拓扑与运行需求，具备较强的工程推广潜力。

二阶锥松弛化的无功优化代码 考虑储能系统ESS，有载调压分接头OLTC，离散电容器CB和静止无功补偿SVC设备约束 以购电和网损成本最小为目标 采用IEEE 33节点算例仿真 使用二阶锥松弛将模型线性化 使用yalmip-cplex/gurobi求解器编写