简介交直流混联电力系统潮流计算与机组组合相关专业文献面向电力系统运行分析与调度优化领域的研究人员和工程技术人员旨在解决大规模交直流混联网络中传统交流潮流算法难以直接适用、机组组合需兼顾直流输电约束等实际问题。整个文档共1个PDF文件压缩包约2.13MB便于下载后直接阅读。文中系统梳理了两类主要潮流算法统一迭代法及其通过简化Jacobian矩阵、改进牛顿法来降低计算量的优化方向交替迭代法结合快速解耦或牛顿法改善收敛性的处理手段同时针对安全约束机组组合介绍了基于Benders分解的主子问题交替迭代以及分层抽样场景修正、惩罚项处理负荷不确定性等思路。阅读后可掌握交直流混联系统在建模变量、算法选择及收敛性改进上的关键要点为相关研究或工程实践提供参考。已有113人学习/下载。1. 交直流混联电力系统潮流计算与机组组合先搞清楚为什么这两件事必须放在一起算交直流混联电力系统潮流计算与机组组合这类方案真正要解决的是一个大电网调度里的脱节问题直流通道的功率是可控的却经常在潮流计算里被当成固定注入、在机组组合里被当成固定边界结果算出来的机组启停计划一遇到直流功率调整或闭锁就翻车。适合它的场景很明确——送端新能源大规模外送、受端多馈入直流、LCC和VSC混合接入的工程。这里从换流站稳态模型讲起给出可复现的交直流潮流计算与机组组合建模路径并把参数设置和踩坑经验一并说明。读完你就能用一套小规模交直流混联算例跑通从潮流计算、机组组合到结果验证的完整闭环。2. 交直流潮流计算换流站模型、控制方式优先级与等效注入法2.1 交直流潮流计算的关键换流站是交流侧的“节点”不是直流侧的“负荷”潮流计算的本质是求一组满足节点功率平衡的电压幅值和相角。交流系统里每个节点只有有功和无功两个方程但直流通道一旦接入换流站就引入了直流电压、直流电流、触发角或熄弧角这些新状态量。你不能把一条直流线路简单塞进交流节点里当负荷——虽然很多人一开始就是这样干的算出来的结果平时看着差不多换流站无功特性一暴露就偏差很大。LCC电网换相换流器的稳态特性通常用准稳态方程描述直流电压与交流侧电压、触发角、换流变漏抗和直流电流直接相关。落到潮流计算里你真正需要关心的是四个量直流电流 Id、直流电压 Ud、触发角 α逆变侧为熄弧角 γ和交流侧电压 E。定功率方式下换流站通过控制 Id 维持 PUd×Id同时交流侧会吸收无功吸收量与直流功率和熄弧角强相关。所以交直流潮流计算必须把换流站的无功消耗建模进去否则受端母线电压会被明显拉低。VSC电压源换流器则简单很多有功和无功可以独立控制潮流建模时直接当 PQ 或 PV 节点处理不需要解熄弧角和无功的耦合关系。这就是为什么柔性直流接入后潮流计算反而比常规直流更顺手。如果你的系统是 LCC 和 VSC 混合接入建议把两类换流站分开建模不要用统一等效模型糊弄过去。2.2 LCC与VSC的控制方式怎么选四类经典控制与定功率优先原则LCC 工程里最常见的控制方式有四类定功率、定电流、定熄弧角、定电压。两端必须配套一般送端定功率或定电流受端定熄弧角或定电压组成一组完整的控制对。控制方式不匹配是潮流发散的头号原因。比如送端定功率、受端也定功率直流电压就失去参考点迭代过程会出现冗余或矛盾送端定电流、受端定熄弧角则是经典组合收敛性最稳。VSC 的控制方式选择逻辑类似两端 VSC 必须有一端定直流电压另一端才能定有功功率如果两端都定有功直流电压飘到天上去也没人管。实际工程里送端 VSC 常设定有功功率受端 VSC 常设定直流电压和定交流电压。这个原则和纯交流系统里必须有平衡节点的逻辑是一回事。我的习惯是先把控制方式按“定功率优先”排好——只要系统里有功率计划就把定功率放在送端受端负责定电压或定熄弧角。如果算例里出现两端都是定功率的情况直接调整控制方式而不是去调初值这是最省时间的做法。2.3 交替迭代法实现等效注入法的最小可跑流程交直流潮流的求解思路分两类统一迭代法和交替迭代法。统一迭代法把直流方程并入交流雅可比矩阵一起求解收敛性好但程序实现复杂商业软件里普遍采用。交替迭代法把直流系统当作送端的等值负荷和受端的等值注入先跑交流潮流再用直流方程修正直流功率和无功消耗两边交替迭代直到收敛。做规划分析或快速评估时我一般用交替迭代法代码改动量小物理概念也清楚。下面是最小化的等效注入法演示。这里用 Pypower 的 IEEE 9 节点算例做底网在第 4 号母线和第 7 号母线之间模拟一条直流通道。# 交直流潮流计算的等效注入法第一步Pypower import pypower.api as pp # 1) 以IEEE 9节点算例为基础 ppc pp.case9() # 2) 指定换流母线编号母线4为送端母线7为受端 bus_sender 4 bus_receiver 7 P_dc 70.0 # 直流通道传输功率单位MW # 3) 送端增加等值负荷受端减小等值负荷 # bus矩阵第3列(Pd)加负荷第4列(Qd)加无功 ppc[bus][bus_sender - 1, 2] P_dc ppc[bus][bus_receiver - 1, 2] -P_dc # 4) 换流站无功消耗近似按直流功率的40%估算 Q_comp 0.4 * P_dc ppc[bus][bus_sender - 1, 3] Q_comp ppc[bus][bus_receiver - 1, 3] Q_comp # 5) 交流牛顿-拉夫逊潮流计算 res pp.runpf(ppc) print(success:, res[success]) print(iterations:, res[iterations])这段代码完成的是等效注入法的第一步把直流功率折算成交流侧的注入。逻辑不复杂但有两个关键点要说明。送端母线加 P_dc 的等值负荷含义是直流通道从送端“拿走”了这么多功率受端母线减 P_dc 的等值负荷相当于直流通道向受端“注入”了这么多功率。这样一来直流通道对交流系统的影响就被嵌入到节点功率平衡方程里。Q_comp 按直流功率的 40% 估算换流站无功消耗这是工程上的粗算经验值精确做法是用准稳态公式计算实际项目中应按换流变漏抗、触发角和直流电流逐点计算。更完整的交替迭代要求你把这段代码放进循环里每次跑完交流潮流后用直流方程反算 P_dc 和 Q_comp再回到第 3 步继续迭代直到直流功率修正量小于 0.1MW 或换流母线电压变化小于 0.1kV。如果迭代振荡就在修正量前乘一个 0.3~0.5 的阻尼因子。2.4 收敛判据与参数设置残差、阻尼因子和无功补偿的工程经验交替迭代法最大的问题是容易出现“交流收敛、直流不收敛”的假象。只盯着交流潮流残差不够换流站边界的直流功率、直流电压和控制角都必须在收敛判据里。常见的收敛阈值组合是交流节点有功无功不平衡量小于 1e-5 p.u.直流功率修正量小于 0.1MW换流母线电压变化量小于 0.1kV。三个条件同时满足再判收敛。参数作用工程经验值阻尼因子控制直流功率和无功消耗的修正步长0.3~0.5振荡时降到 0.2换流站无功消耗初值等效注入法中 Qd 的估算比例0.3~0.5 × P_dc熄弧角LCC 逆变站控制目标12°~15°收敛阈值直流功率修正量0.1 MW无功补偿这块值得多说一句。LCC 换流站本身消耗大量无功实际工程里两侧都会配交流滤波器组和并联电容器潮流计算时不能只建模无功消耗还要把补偿装置同时建进去。等效注入法里最简单的处理是把换流母线当成 PQ 节点在 Qd 中先计入换流站无功消耗再按补偿容量在母线导纳矩阵里加等值电抗。如果计算出来的受端母线电压低于 0.95 p.u.先查无功补偿容量够不够再查控制方式对不对顺序不要反。3. 机组组合建模直流功率从固定边界变成决策变量的关键改动3.1 目标函数怎么改直流功率进入成本项与调整惩罚机组组合原本是纯交流系统的优化问题给定负荷曲线决定每个时段开哪些机组、出多少力。直流通道一进来问题就变了——直流功率是可控的而且可控范围、调整速率、故障后的功率归零特性都与常规机组完全不同。把直流功率当作固定边界条件去算结果就是调度方案在边界条件稍微一变时就不可行。在送端直流功率相当于电能外送所以在目标函数里可以把它视为“送端负荷”的组成部分在受端直流功率等效为电源降低受端本地机组的出力需求。直流功率本身不直接进入目标函数的成本项它只通过平衡约束影响机组出力。若要让直流通道承担调节作用可以对直流功率偏离计划值的偏差加二次惩罚项让求解器尽量贴近交易计划而不是任由直流功率在相邻时段跳来跳去。需要提醒一点不要轻易把“输电价格”加进目标函数除非有明确的跨区交易合同。我见过很多次在目标函数里加了单位输电成本之后优化结果为了压低直流功率而多开受端高价机组调度员完全不认可这种方案。成本项只放机组发电成本和启停成本直流功率通过约束来管理是更稳妥的建模方式。3.2 约束条件怎么补送受端平衡、直流上下限与调整速率多时段机组组合的目标函数统一写成 min ∑t∑i(ci·P_i,t SU_i·y_i,t SD_i·z_i,t)约束条件则要按送端和受端分别列写。送端功率平衡为 ∑P_S,i,t L_S,t P_dc,t P_loss,t受端功率平衡为 ∑P_R,i,t L_R,t − P_dc,t。P_dc,t 是同一个变量把两个子系统的平衡方程串在一起这就是交直流混联机组组合区别于纯交流问题的核心。直流功率上下限要区分静态和动态静态上限是额定容量动态上限还要考虑交流母线电压水平和交流滤波器投退组数。实际运行中直流通道的可用功率常常低于额定值比如受端交流电压偏低时逆变站可能被迫降功率运行。建模时把动态上限做成 P_dc,t ≤ P_dc,max,t每个时段的 P_dc,max,t 由外部数据表给定比常数上限可靠得多。调整速率约束也很容易被忽略|P_dc,t − P_dc,t−1| ≤ ΔP_dc。虽然直流的功率调节速度比火电快得多但日前计划里一般按保守值设置工程上有取 ±100MW/h 的也有取 ±150MW/h 的取决于直流工程的实际运行规程。这个约束不写优化出的直流功率曲线可能出现锯齿状跳变实际调度根本执行不了。3.3 安全约束机组组合SCUC直流闭锁预想故障的建模方式纯交流系统的 SCUC 预想故障集主要是线路和变压器 N-1交直流混联系统则必须把直流闭锁、换相失败后的功率恢复、多回直流同时馈入场景下的同送同受闭锁加进去。直流闭锁对送端意味着窝电对受端意味着大功率缺额两端都需要重新安排出力。工程上最常见的落地方式是故障后再调度可行性检查固定正常态的机组启停状态对每个预想故障检查是否存在一套再调度出力或由备用直接覆盖使系统恢复功率平衡且不越限。这组约束不改变机组启停决策变量只对出力变量的可行域做限制求解器可以一次性处理。严格的两阶段做法是 Benders 分解或鲁棒优化把故障态再调度作为第二阶段决策嵌入机组组合。对大多数工程项目先做“固定启停下的再调度可行性检查”就够用了项目周期紧的时候不要一上来就上两阶段鲁棒。一个容易遗漏的场景是多回直流同时馈入受端电网某一回直流闭锁后功率会转移到其它直流通道和交流联络线上可能造成新的断面越限。预想故障集里不仅要放单回闭锁还要放“直流闭锁后功率转移导致另一回直流过载”的组合故障。这个细节直接决定 SCUC 结果能不能通过调度运行方式校核。3.4 最小可跑示例两节点送受端机组组合的Gurobi实现下面是一个两节点交直流机组组合的最小示例送端 S 和受端 R 之间只有一条直流通道送受端各有两台机组目标是同时优化机组启停、出力和直流功率。用 Gurobi 求解。# 两节点交直流机组组合最小示例Gurobi import gurobipy as gp from gurobipy import GRB # 机组数据送端S受端Rc为边际成本ub为出力上限min为最小技术出力 gen { S1: {bus: S, c: 80, ub: 500, min: 50}, S2: {bus: S, c: 120, ub: 300, min: 30}, R1: {bus: R, c: 100, ub: 450, min: 40}, R2: {bus: R, c: 150, ub: 250, min: 20}, } # 开机固定成本单时段示例中用于体现启停代价 startup {S1: 800, S2: 600, R1: 500, R2: 700} L_S 400.0 # 送端负荷 L_R 600.0 # 受端负荷 m gp.Model(acdc_uc_2bus) # 出力变量下界设为0最小出力由启停约束决定 p {} u {} for g, d in gen.items(): p[g] m.addVar(lb0, ubd[ub], namefp_{g}) u[g] m.addVar(vtypeGRB.BINARY, namefu_{g}) # 直流功率优化变量上限350MW P_dc m.addVar(lb0, ub350, nameP_dc) # 目标发电成本 开机成本 m.setObjective( gp.quicksum(gen[g][c] * p[g] for g in gen) gp.quicksum(startup[g] * u[g] for g in gen), GRB.MINIMIZE ) # 送端平衡本地发电 本地负荷 直流外送 m.addConstr( gp.quicksum(p[g] for g in gen if gen[g][bus] S) L_S P_dc, balance_S ) # 受端平衡本地发电 直流受入 本地负荷 m.addConstr( gp.quicksum(p[g] for g in gen if gen[g][bus] R) L_R - P_dc, balance_R ) # 出力上限与启停变量联动最小出力约束 for g, d in gen.items(): m.addConstr(p[g] d[ub] * u[g], fcap_ub_{g}) m.addConstr(p[g] d[min] * u[g], fcap_lb_{g}) m.optimize()这个示例里有个很多新手会踩的坑出力变量的下界直接设成了机组最小技术出力。如果 u0则 p ≤ 0 与 p ≥ min 直接冲突模型不可行。正确做法是变量下界设 0最小出力通过 p ≥ min×u 这个约束来强制。这样机组停机时出力被压到 0开机时才受最小技术出力限制。送端平衡约束把直流外送当作送端负荷受端平衡约束把直流受入当作受端电源同一个 P_dc 变量在两条约束中符号相反这就把两端的优化问题耦合在了一起。单时段模型的开机成本退化为固定成本扩展到 96 时段时启动成本要写成 SU_i×(1−u_i,t−1)×u_i,t同时补爬坡约束和最小启停时间约束。这个最小模型的价值在于把交直流机组组合的骨架搭清楚工业级模型无非是在这个骨架上不断加约束。4. 算例设计与结果验证能通过调度员质询才算完成4.1 基准算例怎么搭IEEE-RTS加一条直流通道没有完整电网数据时我一般用公开的 IEEE RTS-24 或 IEEE 39 节点系统做底网从中选两个母线作为送端和受端加入一条直流通道。RTS-24 的负荷分布清晰火电水电比例合理加一条直流后既有交流联络线又有直流通道能同时验证两种输电方式下的潮流和机组组合行为。项目数值/方式基准网架IEEE RTS-24发电机数按 RTS-24 原数据直流通道送端母线至受端母线额定 1000MW直流控制方式送端定功率受端定熄弧角无功补偿换流站两侧各配 200Mvar 交流滤波器备用率送端 8%受端 10%直流调整速率±100MW/h母线编号这里不写死因为你手上的 RTS-24 算例文件版本不同母线编号可能有差异。关键是送端要选在电源富余区域受端要选在负荷密集区域否则直流功率优化结果缺乏物理意义。如果送端选了纯负荷节点受端选了纯电源节点模型会为了降低成本让直流功率跑满完全背离实际工程的物理约束。4.2 结果校验的四条铁律机组组合算完不是看目标函数值就完事要过一遍结果校验而且标准要定得足够严格否则方案到调度员手里根本站不住。校验项通过标准不通过的排查路径全局有功平衡偏差 ≤ 0.1MW检查直流损耗公式、送受端负荷数据潮流残差母线有功/无功不平衡 ≤ 1e-5 p.u.检查无功补偿与换流站控制方式机组出力范围无机组低于最小出力运行检查启停状态与出力约束联动直流功率轨迹相邻时段调整量 ≤ ±100MW/h检查直流爬坡约束或目标函数缺少惩罚项第一条铁律是全局有功平衡。把送受端所有机组出力加总减去所有负荷、直流损耗和网损偏差必须在 0.1MW 以内。超过这个量基本可以断定模型中某个等式约束写错了方向最常见的是直流损耗重复计算或送受端负荷数据抄错。第二条铁律是机组出力范围。优化结果里任何一台机组都不能低于最小技术出力运行。出现这种情况通常是启停状态与出力约束的联动写漏了。第三条铁律是机组不能频繁启停同一台机组在相邻几个时段内反复启停说明目标函数里缺少启停惩罚或启停成本设置过低。第四条铁律是直流功率轨迹要平滑不能出现锯齿状跳变这也是判断直流调整速率约束是否生效的标志。4.3 极限方式测试满送、零送、半送三种边界各跑一遍常规校验通过之后还要做极限方式测试。把直流功率的上下限分别设为 0 和额定值加上半功率工况三种方式各跑一遍机组组合对比结果差异。零送时受端所有负荷由本地机组承担送端富余功率全部本地消纳此时能看出受端本地电源是否足够满送时受端本地机组可以压到最小开机方式能看出受端对直流的依赖程度半送则是最贴近日常运行的中间状态考验直流功率和机组出力的协同调节能力。把 24 小时按 96 点展开看曲线这个习惯我建议保留。看三样东西直流功率曲线是否平滑、机组启停状态是否稳定、送受端边际电价是否出现异常尖峰。如果直流功率曲线在某几个时段出现剧烈波动但负荷曲线没有对应的剧烈波动说明目标函数或约束条件里有地方写得不合理。调度员质询的时候通常最关心的就是“这台机组为什么必须开”“直流为什么不满送”“闭锁后受端靠什么顶住”这三类问题都能从 96 点曲线和校验结果里找到答案。5. 避坑交直流混联建模与计算的6条血泪经验5.1 潮流不收敛换流站控制方式冲突的排查顺序现象加入直流通道后交流潮流迭代发散残差停在 1e-2 级别不再下降怎么调初值都没用。原因八成是换流站控制方式冲突。最常见的是送端定功率、受端也定功率直流电压失去参考点或者受端交流电压偏低导致熄弧角越限换流站无法维持定功率。解决先查控制方式是否满足“一端定功率、另一端定电压或定熄弧角”的配套原则。再查无功补偿容量换流站无功不足会把受端母线电压拉低导致迭代不收敛。最后才考虑调阻尼因子和初值。排查顺序不要反控制方式的问题不解决调初值只是碰运气这就是为什么有时候潮流不收敛看起来像玄学其实每一步都有迹可循。5.2 直流功率结果在相邻时段跳变MIP的振荡解现象机组组合结果里直流功率在相邻两个时段从满功率跳到半功率再跳回来受端机组跟着启停运行方式完全不可执行。原因目标函数对直流功率不敏感。直流功率本身不进成本项只要满足平衡约束取哪个值对目标函数影响不大MIP 求解器在多个近似最优解之间随机跳动就产生了振荡。解决给目标函数加直流功率偏差惩罚项偏离计划值越多成本越高同时补直流调整速率约束从可行域上限制相邻时段的跳变量。两招同时用基本能消除振荡。如果仍然有轻微波动可以在求解后做一个“最小调整量”后处理把直流功率曲线进一步平滑。5.3 把直流额定容量当成可用容量动态上限表现象优化结果里直流功率长期顶在额定值但实际运行中直流通道根本无法长时间满功率运行。原因直流可用功率不是常数它受交流母线电压、换流变抽头位置、交流滤波器投退组数等多方面限制。额定容量只是理想条件下的上限直接当约束用等于给模型开了绿灯。解决建一张直流可用功率限值表按送受端电压水平和滤波器投退组数索引每个时段读取对应的动态上限。这个表可以从直流输电工程的运行规程里提取也可以根据历史运行数据统计。建模时把常数上限替换为动态上限求解结果立刻贴近实际。5.4 直流爬坡速率缺失新能源反调峰场景下的翻车现场现象送端光伏大发时直流满送傍晚光伏出力快速下降直流功率却还在高位受端出现电力缺口。原因模型里没有直流功率调整速率约束优化出的直流功率轨迹在现实中跟不上。直流功率调整本身很快但交流滤波器投切和换流站控制方式切换需要时间日前计划里如果不留调整空间实际运行必然翻车。解决在 SCUC 约束中补上 |P_dc,t − P_dc,t−1| ≤ ΔP_dcΔP_dc 按直流工程运行规程设置。校验环节要单独把直流功率轨迹画出来与送端新能源出力曲线叠加对比看是否存在“新能源骤降而直流功率还没降下来”的时段。这个检查十秒钟就能看出问题但很多人根本不做。5.5 预想故障集只有闭锁没有恢复受端两次冲击现象N-1 校核通过了但实际发生直流闭锁后受端电网在故障恢复过程中出现新的断面越限。原因直流闭锁后功率快速降为零但换相失败引起的功率波动和直流功率恢复过程会对受端造成两次乃至三次冲击。模型里只建模了闭锁终态没有建模功率恢复过程的中间态。解决预想故障集加入“直流闭锁 功率恢复”的时间序列至少做两个时间断面的校核闭锁瞬间和恢复中期。受端在闭锁瞬间需要本地备用快速顶上在恢复中期要防止交流断面因功率转移而过载。多回直流馈入的受端电网还要把“一回恢复期间另一回再闭锁”的组合故障一并建模。5.6 数据版本不同步潮流和机组组合算例对不上现象潮流计算用的直流参数和机组组合用的直流参数不一致比如电阻值差 0.5Ω直流损耗差异导致送受端功率不平衡优化结果被调度部门质疑。原因交流网架数据和直流参数通常由不同专业维护一个存在方式数据文件里一个存在直流台账里更新节奏不一致。版本不同步导致两个算例各说各话。解决建一个统一的算例基准文件交流网架、直流参数、负荷水平全部放在一个数据字典里潮流计算和机组组合从一个版本的基准文件读取数据。我习惯在生成算例时打印一份参数指纹包含直流额定功率、电阻、控制方式、送受端母线号两个算例的指纹一致再往下走。这个坑朴实无华但造成的返工量一点不比技术问题少。6. 进阶从单时段优化走向滚动协同与安全校核闭环6.1 多时间尺度滚动日前定直流计划日内修正机组组合从单时段扩展到 96 时段之后下一个值得投入的方向是滚动优化。日前阶段用小时级模型确定机组启停和直流功率计划日内阶段每 15 分钟滚动修正一次直流功率和机组出力。直流功率在日前只做小时级计划在日内做 15 分钟级修正这样既能给送端新能源留出消纳空间又能给受端留出调节余地。两个尺度之间的直流功率偏差需要设置惩罚项否则日内修正幅度过大会与日前计划脱节。6.2 模型验证习惯先在小网络上验证逻辑再上全系统不管做多少进阶功能我的习惯是先在小网络上验证逻辑再把模型搬到完整系统上。两个节点加一条直流的模型能跑通再扩展到 IEEE RTS-24RTS-24 的结果能通过调度员质询再上真实电网数据。规模一大就出现奇怪结果八成是数据版本问题不是模型问题。交直流混联系统的模型复杂度比纯交流系统高了一个量级调试时把问题暴露在小算例里省下来的时间远超建小算例本身的花费。希望这些从实际项目里磨出来的方法和教训能帮你在做潮流计算和机组组合时少走几步弯路。本文还有配套的精品资源点击获取
