最近在整理交直流混合配电网的一组算例顺手把潮流计算的 Matlab 实现重新撸了一遍。做这个方向的人应该都有体会系统里一旦出现换流器 VSC 和多端直流线路传统纯交流潮流的思路就不太够用了因为换流器交流侧的注入功率不再是提前给定的常数它跟直流电压、直流潮流、控制模式全都耦合在一起。我最后采用的是统一求解法也就是把交流节点电压、直流节点电压、换流器功率变量全部放进同一个非线性方程组里用牛顿法一次迭代全部更新。这篇就把背后的建模、方程、代码框架和调试经验完整拆一遍希望能给正在做交直流混合配电网、微电网或者分布式电源接入分析的朋友省点时间。1. 交直流混合配电网潮流计算到底难在哪1.1 纯交流潮流程序遇到的新问题传统配电网潮流计算无论是牛顿法还是前推回代法处理的都是纯交流模型节点分 PQ、PV、平衡三类支路就是线路和变压器。这套体系很成熟Matpower、OpenDSS 里都有现成实现直接调用就行。但交直流混合配电网出现了三样纯交流思路里没有的东西。第一直流网络本身。直流线路没有电抗只有电阻直流节点的功率方程和交流节点完全不是一个形式。交流节点的有功方程是P_i U_i Σ U_j (G cosθ B sinθ)直流节点的方程则简单粗暴得多P U_dc Σ G_dc (U_dc_i - U_dc_j)电压差驱动功率流动。写程序的时候不能复用同一套支路数据。第二VSC 换流器是强耦合接口。它落在交流母线和直流节点之间交流侧的注入有功不是固定值。比如说一个换流器采用定直流电压控制它的交流注入有功完全取决于直流网络到底需要抽出多少功率、换流器自身损耗多大。这类变量在迭代过程中会变而且变的逻辑跟交流系统状态、直流系统状态都相关。第三直流电压是全新的状态变量。它没有相位只有幅值而且标幺值基准和交流电压基准不一样。当你把交流和直流放在同一套方程组里求解首先要解决量纲统一、标幺折算这些细节问题。这些还不是最麻烦的。最麻烦的是换流器的控制模式。配电网里的 VSC 常见的控制方式有定功率、定直流电压、定交流电压、下垂控制等每种控制方式都会改变方程组的变量类型和方程个数。比如定直流电压的直流节点电压值已知不需要直流功率平衡方程而它的换流器交流注入功率反而变成了待求量。这种“角色互换”在纯交流程序里根本不会出现所以刚开始写的时候很容易懵。1.2 交替求解法和统一求解法怎么选解决交直流混合潮流有两条经典路线交替求解法和统一求解法。交替求解法的思路很直接先假定换流器交流侧注入功率已知把 VSC 当作某个交流节点的 PQ 或 PV 注入用常规交流潮流程序解出交流系统状态再根据交流侧结果计算换流器直流侧注入功率解直流网络得到直流电压和直流功率然后把新的直流侧功率折算回交流侧再次解交流潮流。如此往复直到边界功率收敛。这个方法最大的优点是工程上实现快因为可以直接复用成熟的交流潮流程序模块化程度高。我一开始也是走这条路线但很快就撞了墙交直流耦合比较强的时候比如多个 VSC 同时承担直流电压控制、或者带有下垂控制外环迭代经常出现振荡怎么调阻尼都不收敛。后来翻了一些文献才意识到交替法本质上是在两个子系统之间做不动点迭代它没有利用交直流耦合的梯度信息收敛性依赖于交流系统和直流系统的解耦程度。配电网本身 R/X 比大、电压支撑弱比输电网更容易出问题。统一求解法是另一条路把交流节点电压幅值/相角、直流节点电压、每台换流器的功率注入变量全部列为未知量把所有不平衡方程写成一个大的非线性方程组用牛顿法同时迭代修正。这样每一次迭代都是对全系统状态的同时更新交直流接口的耦合关系被完整地保留在雅可比矩阵里收敛性和鲁棒性比交替法高一个量级。两种方法的对比如下对比维度交替求解法统一求解法编程难度较低模块化清晰较高需要统一列写方程和雅可比收敛性弱耦合时尚可强耦合时易振荡牛顿法特性收敛性好对现有程序复用可以直接复用交流潮流程序需要重新建立统一框架处理下垂控制外环嵌套复杂自然融入方程组适合场景早期验证、解耦较强的场合多端直流、强耦合、控制复杂的配电网我的建议是如果只是临时算一个交直流耦合很弱的简单算例交替法可以快速出结果但如果要做配电网规划、多端直流、下垂控制或者反复调用潮流计算做优化直接上统一法别在交替法上浪费时间。1.3 配电网场景的特殊性统一求解法在输电网交直流潮流里用得很多但这套方法论搬到配电网还有两个额外的问题要注意。第一配电网的 R/X 比通常很高线路电阻不能忽略甚至有电阻大于电抗的情况。这种情况下交流潮流方程的病态程度比输电网严重牛顿法对初值更敏感。第二配电网里的分布式电源和负荷种类多样光伏、储能、直流充电桩、变频设备各自控制特性不同负荷往往还要考虑电压静特性。这些都会增加潮流方程的复杂度但不会改变统一求解法的基本框架只需要在节点功率方程里叠加相应的功率项即可。另外配电网的规模通常不大几十到几百个节点系统维度完全可以支撑完整的牛顿法。这跟输电网动辄上千节点还不一样配电网用统一求解法矩阵求逆、迭代计算这些压力并不大反而是建模灵活性更重要。所以我的结论很明确配电网交直流混合潮流统一求解法是最合适的选择。2. 换流器建模统一方程里的核心接口2.1 VSC 稳态模型的基本方程统一求解法的核心是 VSC 换流器的稳态建模。交流侧和直流侧之间的功率、电压关系如果不写对后面的雅可比矩阵全都是错的。在潮流计算这个尺度上我们不需要关心换流器内部的开关过程和谐波只需要保留基波相量层面的等效关系。VSC 可以等效为交流母线经过一个等效阻抗换流变压器漏抗加桥臂电抗连接到换流器桥臂的端口电压源。这个电压源的幅值和相位是可以由 PWM 控制的。设交流母线电压为U_s ∠θ_s换流器端口电压为U_c ∠θ_c等效阻抗为Z_c R_c jX_c定义g R_c / (R_c² X_c²)b -X_c / (R_c² X_c²)那么 VSC 注入交流母线的功率可以写成P_s U_s² g - U_s U_c [g cos(θ_s - θ_c) b sin(θ_s - θ_c)] Q_s -U_s² b U_s U_c [b cos(θ_s - θ_c) - g sin(θ_s - θ_c)]这个方程是后面整个建模的基础。换流器端口电压幅值U_c与直流电压和调制比有关在标幺值体系下可以简化表示为U_c m U_dc其中m是等效调制比。不同调制方式下的系数会有一点差别但潮流计算关心的是功率平衡系数差异经过基准值折算后不影响方程组结构所以代码实现时不需要过于纠结。换流器本身的损耗不能忽略。常见的有功损耗模型是二次多项式P_loss a b |I_conv| c |I_conv|²|I_conv|是换流器交流侧电流幅值可以由视在功率和电压幅值计算|I_conv| sqrt(P_s² Q_s²) / U_s。这部分损耗是潮流计算里经常被忽略但又很重要的细节因为交直流两侧的有功功率差就是靠它连接的。2.2 直流侧功率平衡VSC 交流侧注入有功P_s、损耗P_loss、直流侧注入直流网络的有功P_dc三者之间的关系是一个纯有功平衡方程P_dc P_s - P_loss注意符号方向以交流母线流向换流器为正换流器流向直流网络为正。整流器从交流吸收有功注入直流P_s 0P_dc 0逆变器从直流抽有功送入交流P_s 0P_dc 0。这个式子必须作为独立方程写进统一方程组里因为对于定直流电压控制的换流器P_s是未知的它要靠这条方程衔接上直流网络的功率平衡才能被确定下来。如果漏掉这条方程未知量个数就会多于方程个数雅可比矩阵直接奇异。2.3 常见控制方式的方程组封装换流器控制方式决定了哪些量是已知常数、哪些量是待求变量。配电网 VSC 常见的控制模式有以下几种控制模式已知量待求量直流节点等效类型附加方程定 P 定 QP_s, Q_sU_dc定功率节点无P_s 直接给定定 U_dc 定 QU_dc, Q_sP_s定电压节点直流节点电压已知P_s 由直流平衡决定定 P 定 U_dcP_s, U_dcQ_s定电压节点Q_s 由交流无功平衡决定或由交流电压控制决定下垂控制Q_s, 下垂系数P_s, U_dc下垂节点P_s - P_ref - K(U_dc - U_dc_ref) 0实际工程里一个直流网络必须保证功率平衡所以至少要有一台换流器采用定直流电压控制或者下垂控制否则直流电压会无约束漂移方程组无解或者解无意义。这和纯交流潮流里必须有一个平衡节点是类似的道理交流系统的频率和电压参考点需要有人维持直流系统的电压水平也需要有人维持。2.4 直流网络的节点类型要求这里有一个很隐蔽的坑如果直流网络是分区运行的比如中间通过不同电压等级的直流线路隔开那么每个电气连通的直流孤岛都至少需要一个定电压控制节点。我遇到过这样的情况——整个系统只有一个换流器采用定直流电压控制但直流网络实际上是两个不相连的岛结果雅可比矩阵奇异导致牛顿法直接崩掉。排查了很久才发现是直流拓扑和控制方式配置不匹配。另外直流节点的功率平衡方程形式也要区分节点类型。对定功率直流节点方程是P_dc_k - U_dc_k Σ_j G_dc,kj (U_dc_k - U_dc_j) 0这里G_dc是直流网络的节点电导矩阵P_dc_k是第 k 个 VSC 注入该直流节点的有功。对定电压直流节点这个方程被拆掉直流电压直接写成已知量节点的注入功率通过其它方程最终确定。3. 统一求解法的方程列写与雅可比矩阵3.1 变量与方程怎么“合并同类项”写统一求解法代码前最重要的一步是先明确未知量集合和方程集合做到一一对应。我用一个简单的例子来说明。假设系统有 3 个交流节点节点 1 是交流平衡节点节点 2、3 是 PQ 节点有 2 个直流节点101 和 102通过一条直流线路相连有 2 台换流器VSC1 接在交流节点 2 和直流节点 101采用定直流电压定无功控制VSC2 接在交流节点 3 和直流节点 102采用定有功定无功控制。交流未知量节点 2、3 的相角θ2、θ3和电压幅值U2、U3共 4 个。节点 1 是平衡节点电压幅值和相角都是已知量。直流未知量节点 101 是定电压节点U_dc101已知节点 102 是定功率节点U_dc102未知共 1 个。换流器未知量VSC1 是定直流电压控制交流侧有功注入P_s1是待求量VSC2 是定功率控制P_s2、Q_s2都是已知量不新增未知量。共 1 个。总未知量θ2, θ3, U2, U3, P_s1, U_dc102一共 6 个。对应的方程是交流节点 2 有功和无功平衡方程2 个交流节点 3 有功和无功平衡方程2 个VSC1 的换流器功率平衡方程P_s1 - P_dc1 - P_loss1 01 个直流节点 102 的功率平衡方程P_dc2 - P_net102 01 个正好 6 个方程对应 6 个未知量。这里最大的认知障碍在于VSC1 的交流注入有功P_s1是未知量它出现在交流节点 2 的有功平衡方程里而它的值又要通过换流器功率平衡方程和直流网络功率平衡方程确定。这就体现出“统一”的含义交流变量、直流变量、换流器变量在同一个方程组里互相制约谁也不能单独决定谁。3.2 雅可比矩阵的块结构统一求解法的雅可比矩阵比纯交流潮流大一圈但结构依然清晰。按变量分块可以分成三部分一是交流变量之间的偏导块。这部分和传统交流牛顿法是一样的包括 H、N、M、L 四个子矩阵块对应 P 对 θ 的偏导、P 对 U 的偏导、Q 对 θ 的偏导、Q 对 U 的偏导。唯一的不同是VSC 接入的交流节点其功率方程里可能包含P_s这个换流器变量所以对应位置会多一列非零元素。二是直流变量之间的偏导块。直流功率方程P_dc_k - U_dc_k Σ G_dc (U_dc_k - U_dc_j) 0对直流电压的偏导可以解析推导它依赖于直流电导矩阵元素计算量很小。这部分雅可比相对稀疏也容易验证。三是交直流耦合块。包括换流器功率平衡方程对交流和直流变量的偏导、直流节点功率方程对换流器功率变量的偏导。这些非零元素数量不多但正是它们把两个系统“扣”在一起。我曾经在一篇论文里看到个说法说统一法的雅可比矩阵本质上是一个带边界的块稀疏矩阵内部是交流子矩阵和直流子矩阵边界是换流器耦合块。实际写代码的时候也是这个感觉——先写纯交流部分再写直流部分最后把换流器的接口块补上就能拼出完整的矩阵。3.3 牛顿法迭代流程统一求解法的整体迭代流程和传统牛顿法一致初始化所有未知量计算不平衡向量F(x)检查收敛判据计算雅可比矩阵J(x)求解修正方程J Δx -F更新x x Δx回到第 2 步收敛判据我习惯用双条件电压修正量max|Δx|和功率不平衡量max|F|都小于设定容差。前者保证状态量稳定后者保证功率平衡方程真正满足。有些代码只检查功率不平衡量导致迭代次数很多但电压还在缓慢漂移工程上不够严谨。3.4 初值选择统一法成败的关键统一求解法理论上拥有牛顿法的二次收敛特性但前提是初值要足够好。交直流混合系统里初值问题比纯交流系统更突出因为不仅交流电压幅值、相角需要初值直流节点电压和换流器的待求功率也需要初值。我推荐的初值策略是分层初始化。交流部分采用 flat start平衡节点电压幅值用真实设定值相角为 0其它节点电压幅值取 1.0 p.u.相角取 0。直流部分把定电压节点设为设定值其它直流节点初值统一取 1.0 p.u.。换流器自由注入功率的初值可以先设成 0再通过迭代逐步收敛。如果系统负荷较重或者换流器传输功率较大冷启动可能会不收敛。这时候有个技巧先跑一次不考虑换流器功率注入的纯交流潮流把交流电压结果作为初值再用交流结果计算一个粗略的直流功率初值喂给统一迭代。这样相当于给牛顿法一个“半热启动”实际效果立竿见影。4. Matlab 代码实现与算例验证4.1 数据结构和输入格式写 Matlab 代码的第一步是定义数据结构。参考 Matpower 的风格我用结构体会话来组织数据同时扩展直流部分和换流器部分。% 交流节点矩阵: [编号, 类型, 有功负荷, 无功负荷, 电压幅值初值, 电压相角初值] % 类型: 1平衡节点, 2PQ节点, 3PV节点 bus [ 1, 1, 0.00, 0.00, 1.05, 0.0; 2, 2, 0.50, 0.30, 1.00, 0.0; 3, 2, 0.30, 0.10, 1.00, 0.0; ]; % 交流支路: [起点, 终点, R, X] 标幺值 branch [ 1, 2, 0.020, 0.060; 2, 3, 0.040, 0.080; ]; % 直流节点: [编号, 类型, 电压初值] % 类型: 1定电压节点, 2定功率节点 dcBus [ 101, 1, 1.00; 102, 2, 1.00; ]; % 直流支路: [起点, 终点, R] 标幺值 dcBranch [ 101, 102, 0.020; ]; % 换流器: [编号, 交流节点, 直流节点, 控制模式, 有功设定, 无功设定/直流电压设定] % 控制模式: 1定P定Q, 2定Udc定Q vsc [ 1, 2, 101, 2, 0.00, 1.00; % VSC1: 定Udc1.0, 定Q 2, 3, 102, 1, 0.20, 0.10; % VSC2: 定P0.2, 定Q0.1 ];这套数据结构的好处是扩展性好。以后想加第三台换流器、加直流负荷、加下垂控制参数只需要在对应矩阵里加行或者加列不用改主程序框架。实际上后续如果要跑 IEEE 33 节点改造算例完全可以把该算例的交流 bus、branch 数据塞进来再附加上直流部分和换流器部分。唯一的注意点交流节点编号和直流节点编号不要冲突我习惯把直流节点编号设成 101、102、103 这样的三位数交流节点用 1、2、3 这种个位数一眼就能区分。4.2 不平衡量函数与雅可比构造这套代码里最核心的是两个函数不平衡量函数mismatch和雅可比矩阵函数。不平衡量函数负责根据当前状态计算所有方程残差雅可比矩阵负责计算灵敏度。不平衡量函数的核心逻辑是先由交流节点电压计算交流功率不平衡量再根据 VSC 控制模式计算换流器功率平衡残差和直流节点功率残差最后拼接成一个列向量。function F mismatch(x, bus, branch, dcBus, dcBranch, vsc) % 状态向量 x 的组成: % x [theta_ac; U_ac; P_s_free; U_dc_free] % 这里为了可读性把状态先拆回各物理量 [theta, U, P_s_free, U_dc_free] unpack_x(x, bus, dcBus, vsc); % 1. 交流节点功率不平衡量 S_calc calc_ac_power(bus, branch, theta, U); S_spec -bus(:,3) - 1i*bus(:,4); % 负荷为负注入这里把负荷取负号 F_ac [real(S_calc) - bus(:,3); imag(S_calc) - bus(:,4)]; % 2. 换流器交流注入功率要填进交流功率方程 % 根据控制模式组装 P_s、Q_s然后叠加到对应交流节点上 % 具体代码见下文说明 % 3. 换流器功率平衡方程和直流节点功率方程 % 组装 P_dc、P_loss计算直流网络注入功率 ... end雅可比矩阵我不建议一上来就手推解析表达式因为交直流耦合项的偏导数很容易推导出错而且排查起来很浪费时间。先写一个有限差分版本可以快速验证整体逻辑是否正确function J fd_jacobian(x, bus, branch, dcBus, dcBranch, vsc) n length(x); J zeros(n, n); F0 mismatch(x, bus, branch, dcBus, dcBranch, vsc); h 1e-6; for j 1:n xp x; xp(j) xp(j) h; xm x; xm(j) xm(j) - h; Fp mismatch(xp, bus, branch, dcBus, dcBranch, vsc); Fm mismatch(xm, bus, branch, dcBus, dcBranch, vsc); J(:, j) (Fp - Fm) / (2*h); end end这个函数虽然用了两层循环看起来不高级但对配电网这种节点规模来说完全够用而且调试极其友好。先把有限差分版跑通确认算例结果正确再决定要不要优化成解析雅可比。我在实际项目中也是这么做的有限差分版验证正确之后再逐块替换成解析偏导替换一块验证一块最后整体性能大幅提升同时正确性也有保障。4.3 主循环与收敛控制主循环部分相对简单就是标准的牛顿法加上一个步长限制逻辑function [V, iter, history] unified_pf(bus, branch, dcBus, dcBranch, vsc, tol, maxIter) x init_state(bus, dcBus, vsc); history zeros(maxIter, 1); for iter 1:maxIter F mismatch(x, bus, branch, dcBus, dcBranch, vsc); history(iter) max(abs(F)); if history(iter) tol break; end J fd_jacobian(x, bus, branch, dcBus, dcBranch, vsc); dx -J \ F; % 限制单步最大修正量防止过大的直流电压修正导致负电压 stepLim 0.2; dx(abs(dx) stepLim) sign(dx(abs(dx) stepLim)) * stepLim; x x dx; end end这里限幅操作要特别注意如果直接把超限的元素剪切到 0.2可能会改变修正方向导致残差反而变大。更稳妥的做法是引入一个全局阻尼系数 α先试 α1如果更新后的残差没有下降则 α 减半重试直到残差确实下降。这是牛顿法里标准的阻尼策略。不过对于绝大多数配电网算例限幅在 0.2~0.3 以内就已经足够稳定了。4.4 测试算例与结果分析我用上面这套代码跑了一个小算例3 个交流节点、2 个直流节点、2 台换流器。交流部分节点 1 是平衡节点节点 2、3 带负荷。VSC2 是定功率整流器从交流系统吸收 0.2j0.1 的功率注入直流网络VSC1 是定直流电压逆变器维持直流节点 101 的电压为 1.0 p.u.向交流节点 2 注入有功。粗看功率流向是这样的交流系统通过 VSC2 向直流网络注入有功直流电流从电压较高的节点 102 流向节点 101VSC1 再把直流功率送回交流系统。直流线路电阻 0.02 p.u.在 0.2 p.u. 左右的传输功率下直流电压降落大约在 0.004 p.u. 左右即 U_dc102 略高于 1.0 p.u.。典型收敛结果如下节点电压幅值 p.u.相角 deg备注交流节点 11.0500.0平衡节点交流节点 2约 1.020约 -1.6VSC1 逆变接入交流节点 3约 0.992约 -3.1VSC2 整流接入直流节点 1011.000-定电压节点直流节点 102约 1.004-定功率节点换流器功率方面VSC2 的交流侧有功保持 0.20 p.u. 不变扣除换流器损耗之后注入直流节点 102 的功率约为 0.193 p.u.VSC1 从直流网络抽出约 0.191 p.u.扣除自身损耗后向交流节点 2 注入约 0.184 p.u.。这些数值的精确结果取决于损耗系数的大小但量级和物理直觉完全对得上。迭代次数方面从 flat start 出发通常在 6~8 次迭代内可以收敛到 1e-8 的容差符合牛顿法的预期。比交替法普遍需要的 20~50 次外环迭代要快得多。有意思的是如果把 VSC1 改成下垂控制收敛速度会略慢一点大约 9~10 次但过程的稳定性依然很好不会出现交替法那种外环振荡。这个特性在后续做多端直流和新能源接入优化时非常有用。4.5 扩展从 3 节点到 IEEE 33 节点改造小算例验证通过后往大算例扩展就顺理成章了。IEEE 33 节点配电网是配电网研究里最常用的基准算例把它改造成交直流混合系统通常的操作是选一条支路比如节点 8 到节点 9 这一段替换成直流线路在两端分别接入 VSC或者是在某个节点附近接入一个直流微网通过 VSC 和主网相连。改造时只需要注意三个点。第一交流侧的 bus 和 branch 数据按标准 IEEE 33 节点的 Matpower 格式填好线路阻抗要保持标幺值一致。第二直流部分的功率基准要和交流部分统一否则换流器接口的功率平衡方程会出问题。第三断开的交流支路要把数据从 branch 矩阵里删掉不然会出现交流网络拓扑错误。这套代码结构对这一类改造非常友好因为交流数据和直流数据分开存放换流器矩阵只负责两者的连接关系改拓扑只需要改数据不用动方程逻辑。我实际测试过 IEEE 33 改造为含两段直流馈线的系统总共 35 个交流节点、 4 个直流节点、 4 台换流器规模不大有限差分雅可比版跑一次潮流只需要零点几秒。5. 调试经验、常见坑与性能优化5.1 常见问题速查表统一求解法的代码写完很容易但跑通需要一段调试时间。我把实际踩过的坑整理成一张速查表方便排查。现象可能原因处理方法牛顿法前几步就发散初值太差直流功率方向不合理用交替法粗解做初值或者加阻尼限幅雅可比矩阵奇异直流孤岛缺少定电压节点检查直流拓扑每个孤岛至少配一个定电压或下垂节点直流电压迭代出负数功率方向反了或者初值设得太低检查 VSC 功率符号约定限制直流电压下限迭代次数很多但残差不降步长限制把修正方向改了改为阻尼牛顿法重新计算残差判断下降交流 PV 节点电压失稳VSC 定交流电压接入同一节点两个电压源冲突检查控制模式组合避免同一节点同时两种电压控制换流器损耗为负损耗多项式系数设置不合理检查 a、b、c 系数确保损耗始终为正其中直流孤岛缺定电压节点这个问题最隐蔽。我最初测试时只跑一个简单的单直流网络没发现问题后来把系统改成双直流馈线才发现雅可比奇异排查了很久才定位到是第二个直流孤岛没有下悬挂控制直接把方程组顶爆了。5.2 收敛性调优实操如果遇到收敛困难我的诊断顺序是先看初值再看阻尼最后才怀疑方程列写错误。初值问题是最常见的尤其当换流器传输功率较大时flat start 的P_s0初值可能落在收敛域之外。这时候用前面提到的“交替法粗解做初值”策略最有效代价很小但收益很大。阻尼策略要说得再细一点。简单限幅的问题在于它不保证下降我最终采用的是带线搜索的阻尼牛顿法dx_full -J \ F; alpha 1.0; while alpha 1e-4 x_trial x alpha * dx_full; F_trial mismatch(x_trial, bus, branch, dcBus, dcBranch, vsc); if max(abs(F_trial)) max(abs(F)) break; end alpha alpha * 0.5; end x x alpha * dx_full;这个策略相当于每次迭代都试探性地减小步长确保残差单调下降。虽然多算了几次mismatch但对中等规模的配电网来说完全可接受而且大大提高了算法稳定性。还有一个隐蔽但很实用的诊断技巧打印每次迭代的残差范数观察是否呈现规律的下降趋势。正常牛顿法在收敛域内应该是残差连续下降并在最后两次迭代出现平方级跳变。如果残差出现波动大概率是方程写法或者控制模式处理有误如果残差一直缓慢下降但降不到容差以下大概率是雅可比的某个导数方向算错了。5.3 代码性能优化配电网规模不大但如果在优化算法里反复调用潮流计算性能也需要优化。我的优化顺序是第一步把雅可比从满矩阵改成稀疏矩阵。利用sparse函数预分配非零元素位置矩阵填充完后再每次更新数值。仅这一项就能让较大算例的速度提升不少。代码里只需要把J zeros(n,n)换成J sparse(n,n)后续操作几乎不用改。第二步把有限差分雅可比换成解析雅可比。这一步提升最明显。有限差分版每轮要调用约2n次mismatch函数解析版只需要一次对 200 节点的系统速度差异能达到接近一个数量级。代价是需要逐块推导偏导公式建议边替换边用有限差分验证。我强烈建议保留一个调试开关解析版算完之后和有限差分结果对比一下确认没有笔误再关掉。第三步如果还有性能需求可以用向量化方法同时处理多个潮流的初值或灵敏度但这属于进阶优化了普通场景不需要。5.4 一些更贴近工程的建议这套统一求解法框架的可扩展性其实很好工程上常见的一些需求都能在这个框架里直接加。负荷模型方面配电网节点可能包含 ZIP 负荷也就是恒定阻抗、恒定电流、恒定功率的加权组合。处理方式是在节点功率方程里改成P_load P0 * (a bU cU²)其中abc1对应的雅可比偏导数也要增加一项。这个修改对代码的侵入很小但能显著提高潮流计算对实际负荷的适应性。三相不平衡方面配电网有些场景需要做三相潮流。统一求解法的思路可以平移到三相模型把每个节点扩展成 A、B、C 三相三个节点VSC 模型相应改成三相接口方程数量翻三倍核心逻辑不变。这也是将来做低压交直流混合系统时的主要扩展方向。下垂控制方面多端直流系统里一般会有两台以上换流器参与直流电压支撑这时候下垂控制比单点定电压控制更合理。实现很简单在换流器方程里加一条P_s - P_ref - K(U_dc - U_dc_ref) 0把下垂系数 K 作为参数再让该直流节点电压成为待求变量即可。K 的大小对雅可比矩阵的条件数有直接影响K 过大容易让矩阵病态实际计算时一般取 0.01~0.05 每单位。另外如果要做配电网规划优化潮流计算往往被嵌入到迭代优化里这时候统一求解法还有一个额外好处它输出的雅可比矩阵可以直接用来计算灵敏度比如某个节点注入功率变化对另一个节点电压的影响系数。规划问题里的布点优化、容量优化都依赖这类灵敏度信息等于一次计算两处收益。最后再分享一个小技巧无论你用这套代码做论文还是做工程报告我都会建议把迭代过程中的残差范数变化曲线保存下来。原因很简单这条曲线能直观证明你的算法收敛速度确实符合牛顿法理论审稿人和评审专家都很吃这一套。另外代码里建议把容差设置和迭代上限做成可配置的选项不要写死在脚本里后面做批量仿真的时候就知道省多少事了。我个人的体会是统一求解法从编程角度看起来比交替法要复杂但它的代码结构其实是更整齐的——所有相互作用都在一个雅可比矩阵里你不需要在两个子系统之间反复横跳。这个框架搭好之后加新的控制模式、新的设备模型都只是在原有方程组里加块不需要动整体结构。做交直流混合配电网研究值得把这一步功夫花下去。
交直流混合配电网潮流计算：VSC换流器建模与统一求解法Matlab实现
最近在整理交直流混合配电网的一组算例顺手把潮流计算的 Matlab 实现重新撸了一遍。做这个方向的人应该都有体会系统里一旦出现换流器 VSC 和多端直流线路传统纯交流潮流的思路就不太够用了因为换流器交流侧的注入功率不再是提前给定的常数它跟直流电压、直流潮流、控制模式全都耦合在一起。我最后采用的是统一求解法也就是把交流节点电压、直流节点电压、换流器功率变量全部放进同一个非线性方程组里用牛顿法一次迭代全部更新。这篇就把背后的建模、方程、代码框架和调试经验完整拆一遍希望能给正在做交直流混合配电网、微电网或者分布式电源接入分析的朋友省点时间。1. 交直流混合配电网潮流计算到底难在哪1.1 纯交流潮流程序遇到的新问题传统配电网潮流计算无论是牛顿法还是前推回代法处理的都是纯交流模型节点分 PQ、PV、平衡三类支路就是线路和变压器。这套体系很成熟Matpower、OpenDSS 里都有现成实现直接调用就行。但交直流混合配电网出现了三样纯交流思路里没有的东西。第一直流网络本身。直流线路没有电抗只有电阻直流节点的功率方程和交流节点完全不是一个形式。交流节点的有功方程是P_i U_i Σ U_j (G cosθ B sinθ)直流节点的方程则简单粗暴得多P U_dc Σ G_dc (U_dc_i - U_dc_j)电压差驱动功率流动。写程序的时候不能复用同一套支路数据。第二VSC 换流器是强耦合接口。它落在交流母线和直流节点之间交流侧的注入有功不是固定值。比如说一个换流器采用定直流电压控制它的交流注入有功完全取决于直流网络到底需要抽出多少功率、换流器自身损耗多大。这类变量在迭代过程中会变而且变的逻辑跟交流系统状态、直流系统状态都相关。第三直流电压是全新的状态变量。它没有相位只有幅值而且标幺值基准和交流电压基准不一样。当你把交流和直流放在同一套方程组里求解首先要解决量纲统一、标幺折算这些细节问题。这些还不是最麻烦的。最麻烦的是换流器的控制模式。配电网里的 VSC 常见的控制方式有定功率、定直流电压、定交流电压、下垂控制等每种控制方式都会改变方程组的变量类型和方程个数。比如定直流电压的直流节点电压值已知不需要直流功率平衡方程而它的换流器交流注入功率反而变成了待求量。这种“角色互换”在纯交流程序里根本不会出现所以刚开始写的时候很容易懵。1.2 交替求解法和统一求解法怎么选解决交直流混合潮流有两条经典路线交替求解法和统一求解法。交替求解法的思路很直接先假定换流器交流侧注入功率已知把 VSC 当作某个交流节点的 PQ 或 PV 注入用常规交流潮流程序解出交流系统状态再根据交流侧结果计算换流器直流侧注入功率解直流网络得到直流电压和直流功率然后把新的直流侧功率折算回交流侧再次解交流潮流。如此往复直到边界功率收敛。这个方法最大的优点是工程上实现快因为可以直接复用成熟的交流潮流程序模块化程度高。我一开始也是走这条路线但很快就撞了墙交直流耦合比较强的时候比如多个 VSC 同时承担直流电压控制、或者带有下垂控制外环迭代经常出现振荡怎么调阻尼都不收敛。后来翻了一些文献才意识到交替法本质上是在两个子系统之间做不动点迭代它没有利用交直流耦合的梯度信息收敛性依赖于交流系统和直流系统的解耦程度。配电网本身 R/X 比大、电压支撑弱比输电网更容易出问题。统一求解法是另一条路把交流节点电压幅值/相角、直流节点电压、每台换流器的功率注入变量全部列为未知量把所有不平衡方程写成一个大的非线性方程组用牛顿法同时迭代修正。这样每一次迭代都是对全系统状态的同时更新交直流接口的耦合关系被完整地保留在雅可比矩阵里收敛性和鲁棒性比交替法高一个量级。两种方法的对比如下对比维度交替求解法统一求解法编程难度较低模块化清晰较高需要统一列写方程和雅可比收敛性弱耦合时尚可强耦合时易振荡牛顿法特性收敛性好对现有程序复用可以直接复用交流潮流程序需要重新建立统一框架处理下垂控制外环嵌套复杂自然融入方程组适合场景早期验证、解耦较强的场合多端直流、强耦合、控制复杂的配电网我的建议是如果只是临时算一个交直流耦合很弱的简单算例交替法可以快速出结果但如果要做配电网规划、多端直流、下垂控制或者反复调用潮流计算做优化直接上统一法别在交替法上浪费时间。1.3 配电网场景的特殊性统一求解法在输电网交直流潮流里用得很多但这套方法论搬到配电网还有两个额外的问题要注意。第一配电网的 R/X 比通常很高线路电阻不能忽略甚至有电阻大于电抗的情况。这种情况下交流潮流方程的病态程度比输电网严重牛顿法对初值更敏感。第二配电网里的分布式电源和负荷种类多样光伏、储能、直流充电桩、变频设备各自控制特性不同负荷往往还要考虑电压静特性。这些都会增加潮流方程的复杂度但不会改变统一求解法的基本框架只需要在节点功率方程里叠加相应的功率项即可。另外配电网的规模通常不大几十到几百个节点系统维度完全可以支撑完整的牛顿法。这跟输电网动辄上千节点还不一样配电网用统一求解法矩阵求逆、迭代计算这些压力并不大反而是建模灵活性更重要。所以我的结论很明确配电网交直流混合潮流统一求解法是最合适的选择。2. 换流器建模统一方程里的核心接口2.1 VSC 稳态模型的基本方程统一求解法的核心是 VSC 换流器的稳态建模。交流侧和直流侧之间的功率、电压关系如果不写对后面的雅可比矩阵全都是错的。在潮流计算这个尺度上我们不需要关心换流器内部的开关过程和谐波只需要保留基波相量层面的等效关系。VSC 可以等效为交流母线经过一个等效阻抗换流变压器漏抗加桥臂电抗连接到换流器桥臂的端口电压源。这个电压源的幅值和相位是可以由 PWM 控制的。设交流母线电压为U_s ∠θ_s换流器端口电压为U_c ∠θ_c等效阻抗为Z_c R_c jX_c定义g R_c / (R_c² X_c²)b -X_c / (R_c² X_c²)那么 VSC 注入交流母线的功率可以写成P_s U_s² g - U_s U_c [g cos(θ_s - θ_c) b sin(θ_s - θ_c)] Q_s -U_s² b U_s U_c [b cos(θ_s - θ_c) - g sin(θ_s - θ_c)]这个方程是后面整个建模的基础。换流器端口电压幅值U_c与直流电压和调制比有关在标幺值体系下可以简化表示为U_c m U_dc其中m是等效调制比。不同调制方式下的系数会有一点差别但潮流计算关心的是功率平衡系数差异经过基准值折算后不影响方程组结构所以代码实现时不需要过于纠结。换流器本身的损耗不能忽略。常见的有功损耗模型是二次多项式P_loss a b |I_conv| c |I_conv|²|I_conv|是换流器交流侧电流幅值可以由视在功率和电压幅值计算|I_conv| sqrt(P_s² Q_s²) / U_s。这部分损耗是潮流计算里经常被忽略但又很重要的细节因为交直流两侧的有功功率差就是靠它连接的。2.2 直流侧功率平衡VSC 交流侧注入有功P_s、损耗P_loss、直流侧注入直流网络的有功P_dc三者之间的关系是一个纯有功平衡方程P_dc P_s - P_loss注意符号方向以交流母线流向换流器为正换流器流向直流网络为正。整流器从交流吸收有功注入直流P_s 0P_dc 0逆变器从直流抽有功送入交流P_s 0P_dc 0。这个式子必须作为独立方程写进统一方程组里因为对于定直流电压控制的换流器P_s是未知的它要靠这条方程衔接上直流网络的功率平衡才能被确定下来。如果漏掉这条方程未知量个数就会多于方程个数雅可比矩阵直接奇异。2.3 常见控制方式的方程组封装换流器控制方式决定了哪些量是已知常数、哪些量是待求变量。配电网 VSC 常见的控制模式有以下几种控制模式已知量待求量直流节点等效类型附加方程定 P 定 QP_s, Q_sU_dc定功率节点无P_s 直接给定定 U_dc 定 QU_dc, Q_sP_s定电压节点直流节点电压已知P_s 由直流平衡决定定 P 定 U_dcP_s, U_dcQ_s定电压节点Q_s 由交流无功平衡决定或由交流电压控制决定下垂控制Q_s, 下垂系数P_s, U_dc下垂节点P_s - P_ref - K(U_dc - U_dc_ref) 0实际工程里一个直流网络必须保证功率平衡所以至少要有一台换流器采用定直流电压控制或者下垂控制否则直流电压会无约束漂移方程组无解或者解无意义。这和纯交流潮流里必须有一个平衡节点是类似的道理交流系统的频率和电压参考点需要有人维持直流系统的电压水平也需要有人维持。2.4 直流网络的节点类型要求这里有一个很隐蔽的坑如果直流网络是分区运行的比如中间通过不同电压等级的直流线路隔开那么每个电气连通的直流孤岛都至少需要一个定电压控制节点。我遇到过这样的情况——整个系统只有一个换流器采用定直流电压控制但直流网络实际上是两个不相连的岛结果雅可比矩阵奇异导致牛顿法直接崩掉。排查了很久才发现是直流拓扑和控制方式配置不匹配。另外直流节点的功率平衡方程形式也要区分节点类型。对定功率直流节点方程是P_dc_k - U_dc_k Σ_j G_dc,kj (U_dc_k - U_dc_j) 0这里G_dc是直流网络的节点电导矩阵P_dc_k是第 k 个 VSC 注入该直流节点的有功。对定电压直流节点这个方程被拆掉直流电压直接写成已知量节点的注入功率通过其它方程最终确定。3. 统一求解法的方程列写与雅可比矩阵3.1 变量与方程怎么“合并同类项”写统一求解法代码前最重要的一步是先明确未知量集合和方程集合做到一一对应。我用一个简单的例子来说明。假设系统有 3 个交流节点节点 1 是交流平衡节点节点 2、3 是 PQ 节点有 2 个直流节点101 和 102通过一条直流线路相连有 2 台换流器VSC1 接在交流节点 2 和直流节点 101采用定直流电压定无功控制VSC2 接在交流节点 3 和直流节点 102采用定有功定无功控制。交流未知量节点 2、3 的相角θ2、θ3和电压幅值U2、U3共 4 个。节点 1 是平衡节点电压幅值和相角都是已知量。直流未知量节点 101 是定电压节点U_dc101已知节点 102 是定功率节点U_dc102未知共 1 个。换流器未知量VSC1 是定直流电压控制交流侧有功注入P_s1是待求量VSC2 是定功率控制P_s2、Q_s2都是已知量不新增未知量。共 1 个。总未知量θ2, θ3, U2, U3, P_s1, U_dc102一共 6 个。对应的方程是交流节点 2 有功和无功平衡方程2 个交流节点 3 有功和无功平衡方程2 个VSC1 的换流器功率平衡方程P_s1 - P_dc1 - P_loss1 01 个直流节点 102 的功率平衡方程P_dc2 - P_net102 01 个正好 6 个方程对应 6 个未知量。这里最大的认知障碍在于VSC1 的交流注入有功P_s1是未知量它出现在交流节点 2 的有功平衡方程里而它的值又要通过换流器功率平衡方程和直流网络功率平衡方程确定。这就体现出“统一”的含义交流变量、直流变量、换流器变量在同一个方程组里互相制约谁也不能单独决定谁。3.2 雅可比矩阵的块结构统一求解法的雅可比矩阵比纯交流潮流大一圈但结构依然清晰。按变量分块可以分成三部分一是交流变量之间的偏导块。这部分和传统交流牛顿法是一样的包括 H、N、M、L 四个子矩阵块对应 P 对 θ 的偏导、P 对 U 的偏导、Q 对 θ 的偏导、Q 对 U 的偏导。唯一的不同是VSC 接入的交流节点其功率方程里可能包含P_s这个换流器变量所以对应位置会多一列非零元素。二是直流变量之间的偏导块。直流功率方程P_dc_k - U_dc_k Σ G_dc (U_dc_k - U_dc_j) 0对直流电压的偏导可以解析推导它依赖于直流电导矩阵元素计算量很小。这部分雅可比相对稀疏也容易验证。三是交直流耦合块。包括换流器功率平衡方程对交流和直流变量的偏导、直流节点功率方程对换流器功率变量的偏导。这些非零元素数量不多但正是它们把两个系统“扣”在一起。我曾经在一篇论文里看到个说法说统一法的雅可比矩阵本质上是一个带边界的块稀疏矩阵内部是交流子矩阵和直流子矩阵边界是换流器耦合块。实际写代码的时候也是这个感觉——先写纯交流部分再写直流部分最后把换流器的接口块补上就能拼出完整的矩阵。3.3 牛顿法迭代流程统一求解法的整体迭代流程和传统牛顿法一致初始化所有未知量计算不平衡向量F(x)检查收敛判据计算雅可比矩阵J(x)求解修正方程J Δx -F更新x x Δx回到第 2 步收敛判据我习惯用双条件电压修正量max|Δx|和功率不平衡量max|F|都小于设定容差。前者保证状态量稳定后者保证功率平衡方程真正满足。有些代码只检查功率不平衡量导致迭代次数很多但电压还在缓慢漂移工程上不够严谨。3.4 初值选择统一法成败的关键统一求解法理论上拥有牛顿法的二次收敛特性但前提是初值要足够好。交直流混合系统里初值问题比纯交流系统更突出因为不仅交流电压幅值、相角需要初值直流节点电压和换流器的待求功率也需要初值。我推荐的初值策略是分层初始化。交流部分采用 flat start平衡节点电压幅值用真实设定值相角为 0其它节点电压幅值取 1.0 p.u.相角取 0。直流部分把定电压节点设为设定值其它直流节点初值统一取 1.0 p.u.。换流器自由注入功率的初值可以先设成 0再通过迭代逐步收敛。如果系统负荷较重或者换流器传输功率较大冷启动可能会不收敛。这时候有个技巧先跑一次不考虑换流器功率注入的纯交流潮流把交流电压结果作为初值再用交流结果计算一个粗略的直流功率初值喂给统一迭代。这样相当于给牛顿法一个“半热启动”实际效果立竿见影。4. Matlab 代码实现与算例验证4.1 数据结构和输入格式写 Matlab 代码的第一步是定义数据结构。参考 Matpower 的风格我用结构体会话来组织数据同时扩展直流部分和换流器部分。% 交流节点矩阵: [编号, 类型, 有功负荷, 无功负荷, 电压幅值初值, 电压相角初值] % 类型: 1平衡节点, 2PQ节点, 3PV节点 bus [ 1, 1, 0.00, 0.00, 1.05, 0.0; 2, 2, 0.50, 0.30, 1.00, 0.0; 3, 2, 0.30, 0.10, 1.00, 0.0; ]; % 交流支路: [起点, 终点, R, X] 标幺值 branch [ 1, 2, 0.020, 0.060; 2, 3, 0.040, 0.080; ]; % 直流节点: [编号, 类型, 电压初值] % 类型: 1定电压节点, 2定功率节点 dcBus [ 101, 1, 1.00; 102, 2, 1.00; ]; % 直流支路: [起点, 终点, R] 标幺值 dcBranch [ 101, 102, 0.020; ]; % 换流器: [编号, 交流节点, 直流节点, 控制模式, 有功设定, 无功设定/直流电压设定] % 控制模式: 1定P定Q, 2定Udc定Q vsc [ 1, 2, 101, 2, 0.00, 1.00; % VSC1: 定Udc1.0, 定Q 2, 3, 102, 1, 0.20, 0.10; % VSC2: 定P0.2, 定Q0.1 ];这套数据结构的好处是扩展性好。以后想加第三台换流器、加直流负荷、加下垂控制参数只需要在对应矩阵里加行或者加列不用改主程序框架。实际上后续如果要跑 IEEE 33 节点改造算例完全可以把该算例的交流 bus、branch 数据塞进来再附加上直流部分和换流器部分。唯一的注意点交流节点编号和直流节点编号不要冲突我习惯把直流节点编号设成 101、102、103 这样的三位数交流节点用 1、2、3 这种个位数一眼就能区分。4.2 不平衡量函数与雅可比构造这套代码里最核心的是两个函数不平衡量函数mismatch和雅可比矩阵函数。不平衡量函数负责根据当前状态计算所有方程残差雅可比矩阵负责计算灵敏度。不平衡量函数的核心逻辑是先由交流节点电压计算交流功率不平衡量再根据 VSC 控制模式计算换流器功率平衡残差和直流节点功率残差最后拼接成一个列向量。function F mismatch(x, bus, branch, dcBus, dcBranch, vsc) % 状态向量 x 的组成: % x [theta_ac; U_ac; P_s_free; U_dc_free] % 这里为了可读性把状态先拆回各物理量 [theta, U, P_s_free, U_dc_free] unpack_x(x, bus, dcBus, vsc); % 1. 交流节点功率不平衡量 S_calc calc_ac_power(bus, branch, theta, U); S_spec -bus(:,3) - 1i*bus(:,4); % 负荷为负注入这里把负荷取负号 F_ac [real(S_calc) - bus(:,3); imag(S_calc) - bus(:,4)]; % 2. 换流器交流注入功率要填进交流功率方程 % 根据控制模式组装 P_s、Q_s然后叠加到对应交流节点上 % 具体代码见下文说明 % 3. 换流器功率平衡方程和直流节点功率方程 % 组装 P_dc、P_loss计算直流网络注入功率 ... end雅可比矩阵我不建议一上来就手推解析表达式因为交直流耦合项的偏导数很容易推导出错而且排查起来很浪费时间。先写一个有限差分版本可以快速验证整体逻辑是否正确function J fd_jacobian(x, bus, branch, dcBus, dcBranch, vsc) n length(x); J zeros(n, n); F0 mismatch(x, bus, branch, dcBus, dcBranch, vsc); h 1e-6; for j 1:n xp x; xp(j) xp(j) h; xm x; xm(j) xm(j) - h; Fp mismatch(xp, bus, branch, dcBus, dcBranch, vsc); Fm mismatch(xm, bus, branch, dcBus, dcBranch, vsc); J(:, j) (Fp - Fm) / (2*h); end end这个函数虽然用了两层循环看起来不高级但对配电网这种节点规模来说完全够用而且调试极其友好。先把有限差分版跑通确认算例结果正确再决定要不要优化成解析雅可比。我在实际项目中也是这么做的有限差分版验证正确之后再逐块替换成解析偏导替换一块验证一块最后整体性能大幅提升同时正确性也有保障。4.3 主循环与收敛控制主循环部分相对简单就是标准的牛顿法加上一个步长限制逻辑function [V, iter, history] unified_pf(bus, branch, dcBus, dcBranch, vsc, tol, maxIter) x init_state(bus, dcBus, vsc); history zeros(maxIter, 1); for iter 1:maxIter F mismatch(x, bus, branch, dcBus, dcBranch, vsc); history(iter) max(abs(F)); if history(iter) tol break; end J fd_jacobian(x, bus, branch, dcBus, dcBranch, vsc); dx -J \ F; % 限制单步最大修正量防止过大的直流电压修正导致负电压 stepLim 0.2; dx(abs(dx) stepLim) sign(dx(abs(dx) stepLim)) * stepLim; x x dx; end end这里限幅操作要特别注意如果直接把超限的元素剪切到 0.2可能会改变修正方向导致残差反而变大。更稳妥的做法是引入一个全局阻尼系数 α先试 α1如果更新后的残差没有下降则 α 减半重试直到残差确实下降。这是牛顿法里标准的阻尼策略。不过对于绝大多数配电网算例限幅在 0.2~0.3 以内就已经足够稳定了。4.4 测试算例与结果分析我用上面这套代码跑了一个小算例3 个交流节点、2 个直流节点、2 台换流器。交流部分节点 1 是平衡节点节点 2、3 带负荷。VSC2 是定功率整流器从交流系统吸收 0.2j0.1 的功率注入直流网络VSC1 是定直流电压逆变器维持直流节点 101 的电压为 1.0 p.u.向交流节点 2 注入有功。粗看功率流向是这样的交流系统通过 VSC2 向直流网络注入有功直流电流从电压较高的节点 102 流向节点 101VSC1 再把直流功率送回交流系统。直流线路电阻 0.02 p.u.在 0.2 p.u. 左右的传输功率下直流电压降落大约在 0.004 p.u. 左右即 U_dc102 略高于 1.0 p.u.。典型收敛结果如下节点电压幅值 p.u.相角 deg备注交流节点 11.0500.0平衡节点交流节点 2约 1.020约 -1.6VSC1 逆变接入交流节点 3约 0.992约 -3.1VSC2 整流接入直流节点 1011.000-定电压节点直流节点 102约 1.004-定功率节点换流器功率方面VSC2 的交流侧有功保持 0.20 p.u. 不变扣除换流器损耗之后注入直流节点 102 的功率约为 0.193 p.u.VSC1 从直流网络抽出约 0.191 p.u.扣除自身损耗后向交流节点 2 注入约 0.184 p.u.。这些数值的精确结果取决于损耗系数的大小但量级和物理直觉完全对得上。迭代次数方面从 flat start 出发通常在 6~8 次迭代内可以收敛到 1e-8 的容差符合牛顿法的预期。比交替法普遍需要的 20~50 次外环迭代要快得多。有意思的是如果把 VSC1 改成下垂控制收敛速度会略慢一点大约 9~10 次但过程的稳定性依然很好不会出现交替法那种外环振荡。这个特性在后续做多端直流和新能源接入优化时非常有用。4.5 扩展从 3 节点到 IEEE 33 节点改造小算例验证通过后往大算例扩展就顺理成章了。IEEE 33 节点配电网是配电网研究里最常用的基准算例把它改造成交直流混合系统通常的操作是选一条支路比如节点 8 到节点 9 这一段替换成直流线路在两端分别接入 VSC或者是在某个节点附近接入一个直流微网通过 VSC 和主网相连。改造时只需要注意三个点。第一交流侧的 bus 和 branch 数据按标准 IEEE 33 节点的 Matpower 格式填好线路阻抗要保持标幺值一致。第二直流部分的功率基准要和交流部分统一否则换流器接口的功率平衡方程会出问题。第三断开的交流支路要把数据从 branch 矩阵里删掉不然会出现交流网络拓扑错误。这套代码结构对这一类改造非常友好因为交流数据和直流数据分开存放换流器矩阵只负责两者的连接关系改拓扑只需要改数据不用动方程逻辑。我实际测试过 IEEE 33 改造为含两段直流馈线的系统总共 35 个交流节点、 4 个直流节点、 4 台换流器规模不大有限差分雅可比版跑一次潮流只需要零点几秒。5. 调试经验、常见坑与性能优化5.1 常见问题速查表统一求解法的代码写完很容易但跑通需要一段调试时间。我把实际踩过的坑整理成一张速查表方便排查。现象可能原因处理方法牛顿法前几步就发散初值太差直流功率方向不合理用交替法粗解做初值或者加阻尼限幅雅可比矩阵奇异直流孤岛缺少定电压节点检查直流拓扑每个孤岛至少配一个定电压或下垂节点直流电压迭代出负数功率方向反了或者初值设得太低检查 VSC 功率符号约定限制直流电压下限迭代次数很多但残差不降步长限制把修正方向改了改为阻尼牛顿法重新计算残差判断下降交流 PV 节点电压失稳VSC 定交流电压接入同一节点两个电压源冲突检查控制模式组合避免同一节点同时两种电压控制换流器损耗为负损耗多项式系数设置不合理检查 a、b、c 系数确保损耗始终为正其中直流孤岛缺定电压节点这个问题最隐蔽。我最初测试时只跑一个简单的单直流网络没发现问题后来把系统改成双直流馈线才发现雅可比奇异排查了很久才定位到是第二个直流孤岛没有下悬挂控制直接把方程组顶爆了。5.2 收敛性调优实操如果遇到收敛困难我的诊断顺序是先看初值再看阻尼最后才怀疑方程列写错误。初值问题是最常见的尤其当换流器传输功率较大时flat start 的P_s0初值可能落在收敛域之外。这时候用前面提到的“交替法粗解做初值”策略最有效代价很小但收益很大。阻尼策略要说得再细一点。简单限幅的问题在于它不保证下降我最终采用的是带线搜索的阻尼牛顿法dx_full -J \ F; alpha 1.0; while alpha 1e-4 x_trial x alpha * dx_full; F_trial mismatch(x_trial, bus, branch, dcBus, dcBranch, vsc); if max(abs(F_trial)) max(abs(F)) break; end alpha alpha * 0.5; end x x alpha * dx_full;这个策略相当于每次迭代都试探性地减小步长确保残差单调下降。虽然多算了几次mismatch但对中等规模的配电网来说完全可接受而且大大提高了算法稳定性。还有一个隐蔽但很实用的诊断技巧打印每次迭代的残差范数观察是否呈现规律的下降趋势。正常牛顿法在收敛域内应该是残差连续下降并在最后两次迭代出现平方级跳变。如果残差出现波动大概率是方程写法或者控制模式处理有误如果残差一直缓慢下降但降不到容差以下大概率是雅可比的某个导数方向算错了。5.3 代码性能优化配电网规模不大但如果在优化算法里反复调用潮流计算性能也需要优化。我的优化顺序是第一步把雅可比从满矩阵改成稀疏矩阵。利用sparse函数预分配非零元素位置矩阵填充完后再每次更新数值。仅这一项就能让较大算例的速度提升不少。代码里只需要把J zeros(n,n)换成J sparse(n,n)后续操作几乎不用改。第二步把有限差分雅可比换成解析雅可比。这一步提升最明显。有限差分版每轮要调用约2n次mismatch函数解析版只需要一次对 200 节点的系统速度差异能达到接近一个数量级。代价是需要逐块推导偏导公式建议边替换边用有限差分验证。我强烈建议保留一个调试开关解析版算完之后和有限差分结果对比一下确认没有笔误再关掉。第三步如果还有性能需求可以用向量化方法同时处理多个潮流的初值或灵敏度但这属于进阶优化了普通场景不需要。5.4 一些更贴近工程的建议这套统一求解法框架的可扩展性其实很好工程上常见的一些需求都能在这个框架里直接加。负荷模型方面配电网节点可能包含 ZIP 负荷也就是恒定阻抗、恒定电流、恒定功率的加权组合。处理方式是在节点功率方程里改成P_load P0 * (a bU cU²)其中abc1对应的雅可比偏导数也要增加一项。这个修改对代码的侵入很小但能显著提高潮流计算对实际负荷的适应性。三相不平衡方面配电网有些场景需要做三相潮流。统一求解法的思路可以平移到三相模型把每个节点扩展成 A、B、C 三相三个节点VSC 模型相应改成三相接口方程数量翻三倍核心逻辑不变。这也是将来做低压交直流混合系统时的主要扩展方向。下垂控制方面多端直流系统里一般会有两台以上换流器参与直流电压支撑这时候下垂控制比单点定电压控制更合理。实现很简单在换流器方程里加一条P_s - P_ref - K(U_dc - U_dc_ref) 0把下垂系数 K 作为参数再让该直流节点电压成为待求变量即可。K 的大小对雅可比矩阵的条件数有直接影响K 过大容易让矩阵病态实际计算时一般取 0.01~0.05 每单位。另外如果要做配电网规划优化潮流计算往往被嵌入到迭代优化里这时候统一求解法还有一个额外好处它输出的雅可比矩阵可以直接用来计算灵敏度比如某个节点注入功率变化对另一个节点电压的影响系数。规划问题里的布点优化、容量优化都依赖这类灵敏度信息等于一次计算两处收益。最后再分享一个小技巧无论你用这套代码做论文还是做工程报告我都会建议把迭代过程中的残差范数变化曲线保存下来。原因很简单这条曲线能直观证明你的算法收敛速度确实符合牛顿法理论审稿人和评审专家都很吃这一套。另外代码里建议把容差设置和迭代上限做成可配置的选项不要写死在脚本里后面做批量仿真的时候就知道省多少事了。我个人的体会是统一求解法从编程角度看起来比交替法要复杂但它的代码结构其实是更整齐的——所有相互作用都在一个雅可比矩阵里你不需要在两个子系统之间反复横跳。这个框架搭好之后加新的控制模式、新的设备模型都只是在原有方程组里加块不需要动整体结构。做交直流混合配电网研究值得把这一步功夫花下去。
