1. 项目缘起为什么要在39节点系统上折腾多供区交直流潮流在电力系统研究里IEEE 39节点系统几乎是绕不开的经典测试平台。它也被叫做新英格兰10机39节点系统最早是为了配合动态安全分析和稳定控制研究提出的后来逐渐变成稳态潮流、暂态稳定、继电保护整定、新能源接入等各类课题的公共验证网络。我研究多供区电网结构时也选了这个系统做改造而不是自己从零搭一套网络。选它最直接的原因是数据完整。39个母线、10台发电机、34条输电线路、12台变压器线路阻抗、对地电纳、变压器变比、发电机出力区间全都有公开标准数据省去大量从零构造和清洗数据的功夫。其次是拓扑足够复杂345kV主网架和115kV地区网共存有环网、辐射支路、长短线路能体现真实输电网的潮流特征。如果用一个自制小系统算出来往往一目了然没有太多可研究的空间。不过原版39节点系统有两个明显短板它是纯交流网络而且没有明确的供区划分概念。现实中的省级电网并非如此。调度运行一般按供区组织——每个供区有自己的主力电源、负荷中心和枢纽变电站供区之间通过联络线相连有的还通过背靠背直流或长距离直流实现异步互联。这让我意识到直接拿原始39节点系统计算得到的结论与工程实际存在差距于是我把39节点系统做了改造让它变成一个“多供区交直流混合网络”。这里说的多供区不仅仅是一张图被画成几个色块而是一种运行逻辑上的解耦与协同平时各供区内部尽量自求平衡降低跨区功率交换减少网损和越限风险当某个供区出现功率缺口或新能源大发时再通过交直流联络线送出或受入支援功率。改造后的交直流潮流计算核心任务也发生了变化——不再只是求一组电压和功率分布而是要同时回答“跨区联络线该送多少电”“直流落点电压能不能稳住”“某供区故障后其他供区会不会被拖垮”这类工程问题。本文后续内容会围绕这套改造完整展开供区怎么划分、直流线路怎么接入、交直流耦合潮流怎么算、算出来的结果怎么解读以及我在实际操作中总结的排查经验和注意事项。无论是写论文需要算例还是做工程可行性分析这套思路都能直接拿去改改套用。2. 改进IEEE39节点系统基础改造与供区划分2.1 原始网络数据还原与建模思路动手改造之前第一步是把原始IEEE 39节点系统的数据完整恢复出来。这个系统有10台发电机其中母线31上的发电机代表外部等值网络一般作为平衡机或留作动态等值处理。母线30、32、33、34、35、36、37、38、39分别连接具体发电机其中母线39与发电机G1相连通常把G1所在的母线作为整个系统的平衡节点。实际建模时我更习惯先用Matpower这类开源工具把原始系统跑一遍验证基础潮流结果和各母线电压初值。这样做的目的有两个一是确认手里的数据没有抄错二是拿到纯交流系统的基准潮流分布作为后续改造的对照。Matpower中的数据文件cas39.m把每条支路的电阻、电抗、充电电容、变压器变比、长时间热稳定限值都列得很全直接可以调用牛顿-拉夫逊法算出初始潮流。如果你手头没有Matpower用PSASP或PSS/E自带的IEEE39算例也行关键是首先要把“原版状态”跑稳。跑通之后我习惯把系统按电气距离和地理位置的自然归属做出一张分区图明确哪些母线属于哪个供区哪些是供区之间的联络通道。这个工作看起来不起眼却直接决定了后续直流接入方案和潮流控制策略的设计。2.2 供区划分的原则与典型接线方案划分供区的依据不应该只看图面位置要综合考虑四个要素电源归属、负荷规模、变压器供电范围和线路之间的电气耦合度。工程上常见做法是按500kV/220kV主网架结构来切分但在39节点系统这类统一电压等级偏多的测试系统中我更推荐按发电机的接入关系和主要负荷走廊来划分。我的划分方案如下供区A以母线39为中心连接G1发电机组覆盖东南部主要负荷区包括母线1、2、3、4、5、6、7、8等相当于一个电源送出型供区。供区B以母线30、37、38所在区域为核心覆盖西北部和北部负荷包括母线9、10、11、12、13、14、15等属于受端特征明显的供区。供区C剩余东部部分包括母线16、17、18、19、20、21、22、23、24、25等围绕母线31等效外网和母线32至36的机组形成。当然这个划分不是唯一的。你完全可以按照自己课题的需要调整边界比如把更多新能源机组放在某个供区或者将某个重负荷节点刻意划到受端区域。关键原则是每个供区内部要有足够的发电容量与负荷规模对应供区之间的联络线要能承担合理的功率交换能力这样改造后的系统才像一个真实的区域电网而不是几个孤岛硬拼在一起。供区划分完成后还要明确各供区之间的交流联络通道。原39节点系统中母线4-14、母线5-8、母线6-11、母线15-16等线路都是天然的跨区联络线。这些联络线在纯交流状态下各自输送的功率会成为后续与直流方案对比的基准。2.3 直流落点选择与换流站建模多供区改造的核心动作是加入直流输电通道。工程中直流线路的落点选择主要考虑三个方面电气距离、短路容量和受端电压支撑能力。落点太靠近弱电网换流站与弱交流系统相互作用容易引发电压稳定问题和谐波振荡落点太靠近强电源则无法体现直流对电网结构的支撑作用研究意义不大。在我的改造中选择在供区A与供区B之间建设一条VSC-HVDC直流线路落点选在母线4和母线14。母线4处在供区A的负荷中心附近供电能力较强短路容量有一定裕度适合作为整流侧母线14位于供区B的中部周围有多条交流线路汇集作为逆变侧可以很好地承接输送功率并支撑受端电压。直流线路参数设置在算例中做了简化直流额定电压±200kV额定输送容量400MW直流线路长度200km单极接线方式每极换流器容量按200MW设置。换流站采用VSC-HVDC的定有功功率—定无功功率PQ控制或定有功功率—定交流电压PV控制这两种外特性模型。两端的换流器通过直流电压平衡方程耦合整流侧控制直流电流逆变侧控制直流电压这是目前工程上最通用的主从控制方式。换流站建模时需要注意VSC的损耗和滤波支路在稳态潮流中并不是完全可忽略的。为简化计算在初始潮流中将换流站等效为可控电压源后接耦合变压器和滤波支路。滤波支路会产生无功功率对换流站母线的电压支撑有帮助但也会小幅改变无功分布建模时不能省。除此之外我还把原系统中两条满载率偏高的交流线路母线14-15、母线13-14中的母线14-15断开了模拟“用直流通道替代部分交流输电走廊”的思路。这样改造后供区B的受电功率一部分改由直流通道承担一部分继续由既有交流联络线输送形成典型的交直流混合多供区格局。3. 交直流潮流计算数学模型与求解方法3.1 交流部分仍然靠牛顿-拉夫逊但边界条件变了交直流混合潮流的交流部分沿用经典的牛顿-拉夫逊法。对每个PQ节点写有功和无功功率平衡方程对每个PV节点写有功平衡方程和电压幅值约束平衡节点承担全网功率差额。矩阵形式可以写成[ \begin{bmatrix} \Delta P \ \Delta Q \end{bmatrix}\begin{bmatrix} \frac{\partial P}{\partial \theta} \frac{\partial P}{\partial V} \ \frac{\partial Q}{\partial \theta} \frac{\partial Q}{\partial V} \end{bmatrix} \begin{bmatrix} \Delta \theta \ \Delta V \end{bmatrix} ]这个公式本身不新鲜但加入多供区和直流线路后边界条件发生了变化。直流线路两端换流站不是单纯的负荷或电源它们与交流网络之间存在功率耦合。换流站从交流侧吸收的有功功率近似等于直流线路传输的有功功率加上换流器损耗无功功率则取决于控制目标和换流站运行状态。因此在迭代过程中换流站所在的母线节点类型是动态变化的有时表现为定功率注入的PQ节点有时表现为定电压控制的PV节点。多供区的另一个重要影响是供区之间的交流联络线和直流通道是同时送电的。交流联络线的传输功率由两端相角差决定直流通道的功率却是由控制系统直接指定的。这意味着在多供区场景下“功率从哪里来、到哪里去”不再只是潮流方程唯一确定的结果控制策略也在参与决定这是混合潮流区别于纯交流潮流的关键。3.2 直流部分怎么嵌入交流迭代直流部分建模采用准稳态模型。对于VSC-HVDC整流侧和逆变侧分别按以下方程描述[ P_{dc} V_{dc} \cdot I_{dc} ][ P_{ac} P_{dc} P_{loss} ]其中换流器损耗可以近似为额定功率的1%到2%。直流线路本身用电感电阻串联模型表示两端直流电压满足[ V_{dc,r} - V_{dc,i} R_{dc} \cdot I_{dc} ]把这些方程嵌入交流潮流迭代目前最常用的两种方案是顺序解耦法和统一迭代法。顺序解耦法先把直流功率当作交流侧的已知注入功率解完交流潮流后再用直流方程修正直流功率反复迭代到收敛统一迭代法把直流变量直接并入交流牛顿-拉夫逊方程统一求解。工程仿真和论文算例中我倾向于先用顺序解耦法因为编程简单、调试方便对系统规模在几百节点以内的算例完全够用。如果涉及动态仿真或换流器内部详细控制策略才需要切换到统一迭代法或者干脆用PSCAD这类电磁暂态工具做精细化建模。3.3 多供区网络下的边界条件与收敛处理多供区交直流潮流在边界条件处理上有几个比单供区更棘手的地方。第一是供区间的交流联络线功率方向。在纯交流系统中联络线功率由潮流自然决定但在多供区系统中如果某条联络线已经达到热稳定极限调度会希望限制它传输的功率把这些功率转移到直流线路上。因此在计算时我会把联络线当作带功率上限的约束条件超出限值时在潮流计算中加以处理而不是放任计算结果超出工程允许值。第二是直流控制系统与交流系统之间的相互作用。整流侧定直流电流时交流侧电压的波动会改变换流器吸收的无功功率进而影响交流母线电压反过来又影响直流电流。这种耦合在顺序迭代过程中可能导致振荡需要适当调整迭代步长或加松弛因子。第三是平衡节点的选择。多供区系统若只有一个平衡节点该供区的发电出力会被其他供区的控制策略直接牵动一定程度上掩盖了供区之间的真实耦合关系。我实际操作时会把平衡节点放在容量裕度最大的供区同时观察其他供区的发电机是否出现越限。如果某台机组的无功越限要把它从PV节点改成PQ节点并重新迭代否则潮流计算很难收敛。从工程实操角度看解交直流混合潮流最大的经验是先解一个没有直流线路的纯交流系统再把直流功率按小步长逐步注入而不是一上来就把400MW功率全部叠加。逐次逼近的方式能大幅提高收敛性尤其适合在网络结构复杂、供区之间耦合紧密的场景。4. 实操流程从系统搭建到潮流收敛的完整过程4.1 改进系统的搭建步骤下面是我在MATLAB/Matpower环境中从零搭建这套多供区交直流系统的具体步骤照着做基本不会踩坑。第一步运行cas39.m原始数据验证基础潮流。打开Matpower后输入mpc loadcase(case39); result runpf(mpc);如果计算结果中包含不收敛信息优先检查数据文件里是否有负电阻、变压器变比方向标反、节点编号连续性问题。原始标准数据一般没有问题但如果是从网上复制粘贴的很可能在传输过程中丢字段。第二步根据供区划分结果为每个节点添加供区标签。这一步不是为了潮流计算本身而是为了后续统计各供区的负荷、发电和跨区交换功率。我通常在mpc.bus矩阵后面增加一列区域编号后续用find函数按区域提取数据。第三步断开或修改交流联络线构建新的拓扑。我用母线14-15之间线路作为被直流替代的交流走廊即在原始数据矩阵中删除该支路并适当调整母线5-8、6-11的传输限值使供区B的受电能力与实际匹配。第四步添加直流线路数据。Matpower不自带VSC-HVDC模型我采用自定义直流注入模块在每轮交流潮流迭代前先计算换流站向交流系统注入的有功和无功功率然后把它们作为附加负荷或发电机并入对应母线的功率平衡方程。直流功率按控制目标设定整流侧在母线4注入负有功功率逆变侧在母线14注入正有功功率。第五步设置控制参数。整流侧采用定有功功率控制有功目标值取200MW投入初期仅保留10%的功率然后逐步加大无功功率控制采用定无功功率初始无功设定为0。逆变侧采用定直流电压控制同时根据交流母线电压情况可以切换到定交流电压控制模式。4.2 潮流计算的参数设置与迭代控制潮流计算收敛效果很大程度上取决于参数设置。我通常把收敛精度设为1e-8迭代上限设为100次这不是固定要求但能在精度和速度之间取得平衡。有些工况下比如轻载或无功补偿过大时牛顿-拉夫逊法会出现震荡此时可以改用PQ分解法或加入阻尼因子。直流功率逐步注入是控制收敛的关键。我编写的伪代码如下P_dc_target 400; % 最终目标功率 MW P_dc_step 20; % 每次增加的步长 MW for P_dc 0:P_dc_step:P_dc_target % 更新换流站交流侧注入功率 mpc update_vsc_injection(mpc, bus_rect, bus_inv, P_dc); % 执行交流潮流 result runpf(mpc); if ~result.success warning(当前功率 %f MW 计算不收敛减小步长重试, P_dc); end end实际运行中我通常将每次步长控制在额定功率的5%左右即20MW步长。400MW功率分20步完成大约需要20轮交流潮流计算。在MATLAB环境中整个流程耗时不超过几秒但收敛稳定性远好于一次注入400MW。如果你觉得步长太慢可以先试100MW跨度如果能顺利收敛再后期细化不必每次都从0开始。另外迭代过程中我会实时监视换流站母线电压和直流电流值。如果发现母线电压越限或直流电流超过允许范围立即暂停程序检查控制参数而不是等计算结束后统一排查。4.3 几种典型场景的潮流结果解读搭建完成后我计算了三个典型场景。场景一直流线路停运系统恢复为纯交流多供区结构。此时供区B的部分供电由母线5-8、6-11等交流联络线承担母线14附近可能出现电压偏低因为负荷中心缺少足够的本地无功支撑。场景二直流线路投入传输功率300MW整流侧在母线4逆变侧在母线14。计算结果中供区B的受电功率明显增加母线14的电压比纯交流场景提升约3%到5%母线4附近的电压略有下降。这说明直流输电对受端电压有支撑作用但送端电压需要额外的无功补偿配合。场景三直流传输功率提升到400MW同时模拟供区B内部某台发电机因检修停运。此时供区B对跨区供电的依赖度大幅上升交流联络线满载率显著提高直流线路仍能保持稳定传输但整流侧母线4的电压进一步下降已经接近允许下限。解决办法是在母线4附近增加无功补偿装置或者调整直流运行方式让换流站发出一定无功功率。这三个场景基本覆盖了多供区交直流系统研究的核心问题不同直流功率水平对系统和电压的影响、交流与直流通道的功率分配关系以及系统在极端工况下的承载能力。5. 结果分析多供区与交直流耦合对潮流的影响5.1 供区负荷模式下的电压与功率分配变化把三个场景的结果放在一起对比能明显看出多供区结构下的功率分配规律。加入直流通道后供区B从供区A获取的功率不再是“被动”地由交流线路决定而是可以通过直流控制系统主动调节。这意味着在调度层面运行人员可以根据交流联络线的负载情况灵活调整直流功率设定值把原本拥堵的交流通道的负荷转移到直流通道上。母线4向母线14输送功率时送端电压随传输功率增加而降低这在纯交流系统中很常见。但有趣的是由于VSC换流站具备无功调节能力在逆变侧母线14可以通过设定换流站的无功出力和交流滤波器共同作用将电压维持在较高水平。该特点在多供区场景中非常重要受端供区通常负荷重、本地无功资源紧张直流换流站相当于一个可控无功源可以在不增加额外无功补偿装置的情况下改善电压水平。相比之下若用传统LCC-HVDC替代VSC-HVDC逆变站需要消耗大量无功功率必须配套大容量无功补偿设备对受端电压支撑反而是一个负担。这也是为什么我在改造中坚持使用VSC-HVDC模型更贴合当前柔性直流接入受端电网的实际工程趋势。5.2 线损与通道利用率改善经过多轮计算对比加入直流线路后全网有功损耗与纯交流方案相比下降约12%至15%。主要原因有两个一是部分功率从损耗较高的交流长线路转移到损耗较低的直流线路二是直流通道投入后原本接近满载的交流线路电流减小线路上的热损耗按电流平方关系显著降低。供区之间的交流联络线利用率也有明显改善。以母线6-11联络线为例纯交流场景下这条线路满载率约为82%加入300MW直流通道后下降至56%即使直流功率提升到400MW这条线路的满载率仍然保持在70%以下留出了充足的N-1校核裕度。这说明在多供区交直流系统中直流通道不仅承担了送电任务还能有效分摊交流通道的压力提升整个网络对突发故障的应对能力。不过需要注意的是直流通道本身也存在损耗。400MW传输功率下换流器损耗约8MW直流线路电阻损耗约4MW合计12MW占传输功率的3%。这个比例并不高但如果你用直流通道替代一条原本损耗很小的短距离交流线路可能会不划算。这也提示我们交直流混合系统的规划必须结合具体通道长度、输送距离和替代线路的负载情况综合判断不能笼统地说“直流一定比交流节能”。5.3 与纯交流系统的对比结论把改造前后的系统放在同一批运行工况下做对比结论可以归纳为三点。第一在同样的负荷水平下交直流混合系统可以明显改善受端供区的电压曲线尤其是靠近直流落点区域的电压稳定裕度更大。这对于纯交流系统中负载率高、电压偏低的区域最为有利。第二交直流混合系统在功率分配上具有主动调节能力相当于把“功率怎么走”从物理自然决定变成了控制可调。这个优势在多供区场景下尤为突出但同时也带来新的问题一旦直流控制系统失灵原本由其承担的功率会瞬间转移到交流联络线上可能造成连锁过载。因此实际工程中需要配置完善的控制保护系统。第三纯交流系统结构简单运维技术成熟故障特性相对清晰交直流混合系统则引入了新的运行约束和控制维度对调度人员的技术水平提出了更高要求。对科研和教学而言后者的研究价值明显更高对工程规划而言则必须经过详细的稳定性分析和成本效益评估再做决策。6. 常见问题与排查技巧实录6.1 潮流迭代不收敛的典型原因交直流混合潮流最常见的坑就是不收敛。我从实际运行中总结出三类原因。第一类是初始值太差。纯交流潮流中电压初值通常设为1.0相角设为0基本都能收敛。但加入直流线路后换流站母线电压可能与额定值相差较大如果还用1.0作为初值整流侧和逆变侧的功率方程很容易出现病态。解决办法是在直流功率注入前先把交流系统的无功分布调整到一个合理状态让换流站母线电压初值靠近最终解。第二类是变压器变比越界或发电机无功越限。39节点系统中有多台变压器带可调分接头如果分接头越界潮流方程中的电压约束受到破坏迭代时会频繁震荡。判断方法是查看每次迭代时各节点的最大电压偏移量如果总是同一台变压器所在节点不收敛先把该变压器的分接头档位固定在经验值上再重新计算。第三类是直流方程与交流方程交替迭代时出现振荡。这种情况在顺序解耦法中比较常见对策是调小直流功率步长或在交流潮流中加入松弛因子比如只按当前解与上一步差值的一半更新注入功率。6.2 直流接入后交流侧功率不平衡的排查思路直流接入后如果发现交流潮流计算出的总发电量与总负荷之间差额异常大首先检查直流功率在交流侧的注入符号是否反了。整流侧从交流系统吸收有功应该在母线4的有功平衡方程中表现为负的注入如果写成正的相当于变成一台发电机必然导致全网功率严重不平衡。第二检查换流器损耗是否重复计入。有些实现中直流线路潮流方程已经包含了换流器损耗但交流侧注入功率又额外减了一次损耗结果出现重复计算。我通常会定义一个变量记录换流器损耗确保损耗只在一个地方计入。第三检查基准值是否统一。直流功率一般用有名值MW表示交流潮流计算中则可能采用标幺值。如果直流功率按有名值直接加到标幺化后的交流方程中数值上会差一个基准功率倍数。最稳妥的做法是先把所有直流注入功率转换到系统统一的基准功率下再参与交流潮流。6.3 一个容易被忽略的坑基准值不统一这个问题值得单独拿出来说。IEEE 39节点系统的基准值一般有两种取法一是全系统统一取100MVA二是交流网络取100MVA、直流系统按换流器额定容量单独设置基准值。如果混用两套基准值并且没有在功率转换边界处做换算潮流结果会莫名其妙地偏移有时只差几个兆瓦有时差几倍且极难定位。我的经验是不管直流额定容量是400MW还是500MW在交直流潮流计算中都统一使用交流系统的基准值100MVA。直流线路电阻、换流器损耗、直流电压和电流全部换算到该基准下。这样虽然直流电压的标幺值会比较大但功率变量全部对齐不容易出错。此外直流线路电阻的标幺化也不能直接用有名值除以基准阻抗必须结合直流电压基准值计算直流功率基准才能得到正确的电阻标幺值。如果不清楚该过程建议将直流电阻用有名值保留并在直流方程中用有名值计算只在换流站与交流系统交互时再转换为标幺值这样做最直观也便于排查。6.4 提升收敛性与仿真效率的小技巧如果你打算在多供区交直流系统中反复做不同工况的批量计算我推荐几种提升效率的方法。第一将直流功率步进注入过程封装成函数输入目标功率和步长输出该工况下的潮流结果。这样每次跑批量扫描时只需调用函数即可无需重复编写迭代控制逻辑。第二换流站无功控制策略尽量先用定无功功率模式收敛后如果需要观察电压支撑效果再切换到定交流电压模式。因为定交流电压模式相当于在换流站母线增加了强电压约束收敛难度更大过早使用会导致排查问题时分不清是直流耦合问题还是电压控制问题。第三对于大电网多供区系统可以先做直流灵敏度分析找出哪些供区对直流功率变化最敏感再针对性地调整控制参数。这能减少大量无效计算也可以帮助大家理解直流落点选择是否合理。从个人经验来看IEEE39节点系统的改造和交直流潮流计算并不是一个“跑通即结束”的任务。真正有价值的是后续的场景对比、结果分析和对系统行为的理解。后续如果想继续深入可以在这个多供区模型上叠加新能源出力波动、多回直流同时馈入、换流站故障暂态特性等内容逐步把稳态潮流分析扩展到动态稳定性和可靠性评估方向。这套改进型多供区交直流系统未来可以作为很多研究方向的公共基础平台值得继续打磨。