交直流输电网潮流计算这两年问的人越来越多。以前做潮流分析清一色都是交流节点顶多挂个负荷、设个PV节点现在风电、光伏大规模接入直流汇集、直流外送成为常态常规直流和柔性直流一起往交流电网里怼整个系统变成交流、直流强耦合的混联网络。很多手里攥着成熟交流潮流程序的人突然发现直接拿这套东西去算交直流系统要么算不出来要么结果连功率平衡都对不上。原因不复杂直流线路、换流站、控制方式这些对象的方程传统交流潮流模型里根本没有位置。这就要回到交直流输电网潮流计算的两种经典思路——统一法和交替法。统一法把交流、直流、换流器方程全部塞进一个大雅可比矩阵里联立求解数学上最干净收敛性也最可靠交替法把交流潮流和直流子系统拆开各自求解、来回迭代编程上最省事能继续复用成熟交流潮流程序。我在实际项目里两种方法都写过也踩过不少坑这篇文章就把我对这两种方法的理解、实现细节、参数设置和排错经验全部倒出来给正在做交直流系统研究、规划或者程序开发的朋友一个能直接上手的参照。1. 交直流混联系统潮流计算的本质与两条技术路线1.1 为什么交直流混联系统不能拿纯交流潮流程序硬算要理解统一法和交替法先得搞清楚交直流混联系统给自己加了什么麻烦。交流潮流的核心是节点功率平衡方程对每个节点列写有功和无功两个方程未知量是节点电压幅值和相角。到了交直流系统里换流站接入交流母线后那根母线上的注入功率不再只由发电机和负荷决定还取决于换流站从交流侧吸取还是注入多少功率。更关键的是换流站的行为不是随便定的它受直流侧运行状态和控制策略约束。举个具体例子一个两端常规直流系统整流站和逆变站的直流电压、直流电流通过直流线路耦合在一起整流站可能定功率逆变站可能定熄弧角。要算清楚这两个换流站各自从交流电网吸收多少无功、注入多少有功必须联立求解直流网络方程和换流器方程而这些方程在纯交流潮流框架里根本不存在。有人可能会说那我干脆把换流站等效成一个PQ节点或者PV节点丢进去不就行了在小规模、单端接入、耦合不强的场景下这么近似一下还有一定精度但遇到多端直流、直流电网、柔性直流电压下垂控制这类情况换流站之间的功率分配、电压分布是全网耦合的你连它该等效成什么节点类型都说不清楚更别提算准了。1.2 统一法值得花力气的一次性大联立统一法的核心思想非常直接把交流节点方程、直流节点方程、换流器方程以及各种控制方式方程全部放进同一个非线性方程组用牛顿-拉夫逊法统一迭代求解。也就是说在每一次迭代里不仅要修正交流节点电压的幅值和相角还要修正直流节点电压、换流器直流电流、控制角等直流侧变量。所有方程同时收敛没有谁等谁的问题。这种做法的优点很突出。首先收敛性最好理论上具有牛顿法二阶收敛特性尤其是在交直流强耦合的场景下非常稳定。其次结果一致性高因为交流侧和直流侧是在同一个代数系统里求解不会出现交替迭代那种两边各自收敛但整体不匹配的尴尬情况。缺点也很明显编程复杂度高需要把所有方程和变量理清楚雅可比矩阵的分块、解析求导、稀疏存储都要自己处理对初值敏感初值给不好矩阵直接奇异或者迭代发散调试难度大出了问题你很难一眼看出是交流侧发散还是换流器方程写错了。但说实话如果目标是做科研、写高水平程序、算直流电网这种强耦合场景统一法值得下功夫。我现在的项目里很多核心计算模块就是统一法稳定性和精度都让人放心。1.3 交替法一种成熟的折中和妥协交替法的思路则是典型的工程思维——不要试图一步到位把大问题拆成两个能分别求解的小问题。交替法的做法是初始化时先给交流系统设定换流站的注入功率解一次交流潮流得到换流站接入母线的电压然后把交流电压当作已知量带到直流子系统求解直流网络和换流器方程得到新的换流站直流电压、电流和交流侧注入功率接着再用更新后的交流注入功率去解交流潮流如此反复循环直到交流侧和直流侧变量的变化量都小于给定阈值。这种方法的吸引力在于模块化。交流潮流可以继续用你熟悉的程序或开源工具直流子系统单独写一个模块两边通过换流站这个边界交换数据。对于只改造已有交流程序的工程任务交替法是最稳妥的选择。代价是收敛性没有保证尤其当交流系统比较弱、直流馈入容量占比高的时候交替迭代经常出现振荡甚至发散。此外交替法本质上把交流系统和直流系统当作彼此的黑箱两边各自收敛并不代表整体满足所有方程所以某些极端工况下即使收敛了结果也需要额外校验。1.4 两种方法的直观对比对比维度统一法交替法方程处理方式交流、直流、换流器方程组成统一非线性方程组交流潮流与直流子系统分别求解、边界交互迭代收敛速度迭代次数少通常5~8次左右收敛迭代次数多常见20~50次强耦合时可能不收敛编程复杂度高需要处理矩阵分块和解析求导较低可在成熟交流潮流程序基础上增量开发初值敏感度敏感初值不合理容易矩阵奇异或发散相对不敏感每个子问题都有独立收敛机制对系统规模的适应性适合直流电网、多端直流等强耦合场景适合单端或少量直流馈入的大交流电网运行速度单次迭代耗时长但总迭代次数少单次迭代快但总迭代次数多通信开销大看过这个对比就会明白不是哪个方法绝对好而是看你的场景。后面我会详细讲两种方法的实现细节把关键步骤和公式摆出来这样你可以直接对照搭建自己的程序。2. 统一法把所有方程打包成一个大雅可比矩阵2.1 待求解的方程与未知量梳理做统一法第一步不是写代码而是把未知量和方程数对齐。这一步做不好后面雅可比矩阵必然奇异我见过太多人栽在这里。交直流混联系统的未知量大致分三类第一类是交流系统的常规未知量包括所有交流节点的电压幅值V和相角θPQ节点算多出来的未知量平衡节点只有相角为已知约束。第二类是直流系统的未知量包括每条直流节点的电压V_dc、每条直流线路或换流器的电流I_dc。第三类是换流器的控制变量包括换流器直流电压、直流电流、控制角整流侧触发角α或逆变侧熄弧角γ如果是柔性直流还要包括调制比m和移相角δ。方程也要对应分三类。第一类是交流节点功率平衡方程和纯交流潮流完全相同只是PQ节点中如果接有换流站注入功率会多出换流器交流侧的交换功率项。第二类是直流网络方程根据直流线路电阻列写节点电压方程或支路潮流方程。第三类是换流器方程这部分是统一法的精华也是最容易写错的地方。以常规直流换流站为例需要列写的方程式包括换流器电压方程、直流电流与交流电流的关系方程以及控制方式决定的约束方程。换流器电压方程长这样U_dc (3√2/π) · T · U_ac · cosθ - (3/π) · X_c · I_dc其中T是换流变压器变比U_ac是换流变压器阀侧空载电压θ是控制角X_c是换流变压器漏抗折算到阀侧的等效电抗。这个公式可以在任何一本直流输电教材里找到但真正实现时会发现标幺值换算是一个大坑。如果交流侧和直流侧的基准值选取不一致公式里的系数要重新推导很多初学者的雅可比矩阵奇异就是栽在这里。另外不同控制方式对应不同方程。定功率控制意味着换流站吸收的有功功率等于给定值这个方程可以表达为P_spec - U_dc·I_dc 0定直流电压控制意味着V_dc - V_spec 0定熄弧角控制意味着γ - γ_spec 0。这些控制方程必须和前面的物理方程放在一起联立才能唯一确定换流器的运行状态。2.2 牛顿-拉夫逊统一迭代的完整步骤统一法的求解流程可以按下面这些步骤拆开步骤一初始化。交流节点电压幅值取1.0标幺相角取0度平衡节点除外直流节点电压初始化为额定值换流器控制角初始化为典型运行值比如整流侧取14度左右逆变侧取18度左右。如果知道一种合理运行方式用那种运行方式的值做初值更好。步骤二按当前变量值计算全部方程的不平衡量。这一步要分别计算交流节点有功不平衡量ΔP、无功不平衡量ΔQ、直流节点功率不平衡量ΔP_dc、换流器方程残差ΔR。这些不平衡量的含义是方程左右两边当前差值牛顿法的目标就是通过修正变量把差值压到零。步骤三组装雅可比矩阵。矩阵整体结构是分块的左上角是交流功率方程对交流电压、相角的偏导即传统交流潮流的雅可比矩阵右下角是直流方程和换流器方程对直流变量、控制变量的偏导右上角和左下角是交直流耦合项的偏导比如换流器交流侧功率对直流电压的导数、直流方程对交流电压的导数。步骤四求解线性方程组J·Δx -ΔF得到变量的修正量Δx。这里的Δx是所有待求变量的增量。步骤五更新变量并检查收敛条件通常取变量修正量绝对值的最大值与不平衡量的范数同时小于给定阈值比如1e-6标幺。伪代码大概长这样# 统一法伪代码框架 def unified_pf(ac_system, dc_system, converters, max_iter20, tol1e-6): x initialize_variables(ac_system, dc_system, converters) for k in range(max_iter): F compute_mismatch(x, ac_system, dc_system, converters) if norm(F, np.inf) tol: return x, True J compute_jacobian(x, ac_system, dc_system, converters) dx solve_linear_system(J, -F) x x dx return x, False这段伪代码看起来简单实际工程量全在compute_jacobian里。我的建议是雅可比矩阵中的每一项尽量推导解析表达式不要用数值差分。因为交直流耦合项的数值差分在变量量级差距大的时候误差会放大容易导致迭代假收敛。为了推导方便可以把所有方程先写成关于统一变量向量x的标量函数再用自动微分或者符合作业要求的简单中心差分去验证前几次迭代的Jacobian最后部署时替换成解析表达式。2.3 换流器方程的标幺值处理与雅可比修正换流器方程的标幺化处理我单独拎出来说因为这是统一法最容易翻车的地方。实际计算中交流侧使用的基准值通常是交流系统额定容量S_b和额定线电压U_b而直流侧如果也采用同样的功率基准值那么直流电压基准值需要单独选取通常是直流额定电压直流电流基准值则由功率基准和直流电压基准推导出来即I_dc_b S_b / U_dc_b。在这个基准体系下换流器电压方程的标幺值形式和有名值略有不同重点在于换流变压器变比T的取值。换流变压器分接头调整时T会在一定范围内变化如果程序里把T当作常数和实际情况有偏差如果当作变量需要多设一个未知量并增加分接头控制方程。从潮流计算的工程实践看对于规划研究一般把T设为固定分接头或者简化处理但对于运行分析必须考虑分接头调节对交直流耦合的影响。雅可比矩阵修正量中换流站交流侧注入功率P_ac、Q_ac对换流器变量U_dc、I_dc、θ的偏导是整个矩阵里数值最敏感的地方。功率表达式如下P_ac U_dc · I_dc Q_ac U_ac · I_ac · sinφ其中I_ac是换流器交流侧电流φ是功率因数角。要注意Q_ac并不是U_dc·I_dc的简单函数它还和换流器吸收的无功功率有关通常用Q_ac P_ac · tanφ来计算。φ由控制角和换流器等效电抗共同决定这就在雅可比矩阵里引入了更多耦合项。处理这一块时我习惯先把换流器的某种控制方程比如定电压、定功率代入消去一部分中间变量再对剩下的自变量求偏导矩阵会更紧凑数值性质也更好。2.4 变量编号顺序和初值选择对收敛性的影响统一法收敛调不动的时候首先怀疑变量编号和初值这两个因素远比想象中影响大。变量编号顺序方面我强烈建议把交流变量和直流变量分块排序比如前一部分放交流节点的相角和幅值中间放直流节点电压最后放换流器控制变量。这样做的好处是雅可比矩阵近似块状带状稀疏性更明显直接稀疏LU分解也好、预处理迭代求解也好性能都更好。如果混在一起乱序编号矩阵非零元分布完全不可预知条件数也会变差。程序实现上可以用一个索引字典记录每个变量在总向量中的位置这样组装矩阵时不容易搞错。初值选择方面我踩过一个印象很深的坑。做某次多端柔性直流接入弱交流系统的算例统一法怎么调都不收敛矩阵定性检查也看不出问题后来发现是整流站的直流电压初值给得太离谱导致换流器方程残差巨大雅可比矩阵中相应行的数值完全失衡。换成一个相对合理的初值后一次收敛。从那以后我的习惯是测试新算例前先用交替法跑一遍哪怕不收敛也能拿到一个差不多的近似解再用这个解做统一法的初值。这种方法几乎成了我的标准流程。3. 交替法拆成两个子问题来回迭代3.1 交替法的整体流程设计交替法的流程看起来比统一法简单实际设计时需要想清楚两个子系统之间的数据交换协议。我的实现步骤是这样的第一步初始化换流站交流侧注入功率P_ac、Q_ac。对于定功率站直接按给定功率初始化对于定直流电压站需要先估算一条直流电流再换算出交流侧有功。这一步的粗估算不需要太准确后面会迭代修正。第二步求解交流子系统。把换流站当作特殊的负荷或电源节点交流潮流程序跑一次得到所有交流节点电压特别注意换流站接入母线的电压U_ac。第三步求解直流子系统。把第二步得到的交流母线电压作为已知量代入所有换流器方程和直流网络方程联立求解得到直流电压、直流电流和换流器控制角。这一步的直流网络方程一般是线性方程组因为直流网络是纯电阻网络。第四步更新换流站交流侧注入功率。根据第三步得到的直流功率和换流器无功特性重新计算各换流站注入交流节点的P_ac、Q_ac。第五步判断收敛性。比较本轮和上一轮的直流功率、直流电压和交流母线电压如果变化量都小于阈值则停止否则返回第二步继续迭代。从这个流程可以看出交替法的核心问题是第三步和第四步之间的耦合处理。理论上只需要交换少量变量但实际计算中经常发现两个子系统的接口变量如果更新太激进迭代会震荡如果更新太保守收敛又会很慢。所以接口变量的平滑处理非常重要。3.2 直流子系统的矩阵方程列写与求解直流子系统在交替法中通常是一个纯电阻网络它的节点电压方程可以用直流导纳矩阵来列写。设直流网络有n个节点那么有G_dc · V_dc I_dc其中G_dc是直流网络的节点导纳矩阵V_dc是直流节点电压向量I_dc是节点注入电流向量。这个方程形式上和交流系统牛顿法中的线性方程非常相似但注意直流网络是线性的不需要迭代求解一次线性方程求解就够了。不过I_dc并不是简单给定的它由各换流站的运行状态决定。整流站在定功率控制下P_spec和直流母线电压V_dc之间的关系是I_dc P_spec / V_dc这是一个非线性关系。如果多个换流站都这样列方程直流子系统的求解就变成解一个非线性方程组了。处理办法是对直流子系统再做一个内部的小迭代或者直接把I_dc的表达式中与V_dc相关的部分搬到方程左边形成线性方程组。我自己在实际项目里更喜欢后一种做法比如定功率站的方程写成G_dc(i,:)·V_dc P_spec / V_dc_i然后移项变成(V_dc_i的平方)与其它电压的线性组合再用牛顿法解这个小规模方程。对于柔性直流的多端系统直流网络方程还需要考虑各端控制模式。定直流电压站提供直流电压的基准定功率站按给定功率运行如果有下垂控制站还要把下垂系数写进方程里。比如某端电压变化时它的功率会按P P_0 k·(V_ref - V_dc)调整这个方程要作为约束加入直流子系统。否则整个直流电压没有参考点方程组会出现零特征值表现为矩阵奇异。3.3 收敛判据、松弛因子与振荡抑制交替法的收敛判据通常取接口变量的相对变化量。我一般取三个指标同时满足换流站交流侧有功注入变化量小于0.0001标幺直流节点电压变化量小于0.0001标幺交流母线电压幅值变化量小于0.0001标幺。之所以要同时看交流量和直流量是因为只盯一侧变量容易误判某些工况下直流功率已经稳定但交流电压还在缓慢漂移。振荡抑制是交替法实操中过不去的坎。最常用的手段是松弛因子也就是不直接采用最新计算值更新接口变量而是用旧值和新值的加权平均x_new x_old λ · (x_calculated - x_old)λ松弛因子的取值一般在0.5到0.9之间。系统耦合较强时取0.5左右耦合较弱时取0.9左右。λ太大了容易震荡λ太小了收敛慢。需要说明的是这没有放之四海而皆准的值每换一个算例最好先扫一遍λ看收敛趋势。另外还可以对交流母线电压做低通滤波类似地让接口电压变化更平滑这在弱交流系统中效果明显。我曾经处理过一个弱系统场景交流系统短路容量比SCR只有2.2交替法怎么都收敛不了波形呈现出明显的持续等幅振荡。后来把松弛因子调到0.45同时对整流站的直流电流做了限幅平滑迭代40多步后终于稳定下来。这个经历让我意识到交替法的调试过程更像调PID参数要有点耐心也要有点手感。3.4 为什么交替法在实际工程中依然很有市场统一法理论上更好为什么很多工程场景还是选择交替法答案在于代码复用和软件生态。大多数电网计算软件和现有潮流程序对交流系统模型非常完善发电机模型、负荷模型、变压器调压、无功补偿设备这些都能精细模拟。如果改用统一法就需要把这些模型全部重新实现一遍工作量巨大。交替法则不同交流侧完全不用动只需要额外写一个直流子系统模块每次算完交流潮流后调用一次把结果导来导去就行。在PSASP、BPA、PSS/E这些软件里交直流系统的处理大多也是交替或解耦迭代的思路只是把细节封装在产品里了。所以理解交替法对新版软件里交直流混合仿真卡片参数的理解也有帮助。比如直流控制方式的填卡逻辑、换流站无功上下限设置本质上就是在设定交替迭代时的接口变量更新规则。4. 收敛性、计算效率与场景匹配的实际对照4.1 一个典型算例下的表现对比我用一个比较典型的测试算例说明两种方法的实际表现。测试系统采用修改后的IEEE 14节点交流系统接入一个双端常规直流整流侧定功率控制逆变侧定熄弧角控制直流功率约占系统总负荷的20%。两种方法我都用Python实现同样是1e-6收敛精度初值取平启动。统一法迭代7次收敛每次迭代解18阶线性方程组总耗时大约0.25秒。交替法收敛需要46次迭代每次迭代包括一次快速交流潮流计算和一次直流子系统求解总耗时大约0.6秒。这个对比符合预期统一法迭代次数少但单次迭代开销大交替法单次迭代便宜但次数多。有意思的是如果把这个直流功率占比从20%提高到40%交替法的迭代次数跳跃到80多步统一法依然是8步左右差距进一步拉大。需要说明的是这个数值只是我自己测试算例下的表现不代表所有情况但它反映了两种方法在收敛性上的普遍差异规律。在做方案对比时最好不要只依赖单一算例至少换三四个不同强弱系统的算例看趋势。4.2 大电网场景下怎么组合使用实际工程里单一方法往往不是最优解。我现在的做法是先用交替法做初值探索跑一个大致的潮流分布观察哪些换流站接入点电压偏低、哪些直流线路潮流较重等心里有数了再切换到统一法做最终精确计算。这个组合策略兼顾了交替法的易用性和统一法的稳定性。如果系统非常大比如数千个交流节点的区域电网全项目统一法会面临巨大的雅可比矩阵组装和求解压力。这时候合理的选择是交流系统用交替法、把直流子系统放到更大的联合迭代里或者在统一法的框架下考虑区域解耦把交直流耦合界面上的变量用外部网络等值的方法消去以减小联合矩阵规模。另外对于多端柔性直流和直流电网我倾向于直接上统一法。因为多端直流网络里任何一端的运行方式变化都会迅速传导到整个直流网络换流站之间的耦合非常强交替法几乎必然面临收敛性问题。一个折中的方案是直流网络内部的方程仍然保留联立但交流系统和直流系统的边界继续采用交替迭代这种半统一半交替的做法在某些工程场景下既保证了精度又控制住了计算规模。4.3 校验计算结果合理性的几种方法不管用哪种方法计算结果都必须经过合理性校验。我每次跑完潮流后的固定动作是检查五件事第一全系统功率平衡。交流发电机总出力、交流负荷、直流输送功率、换流站损耗和直流线损所有功率加总应该为0在给定精度内。这一步最容易暴露方程漏项或者正负号错误。第二直流网络KCL。每个直流节点的注入电流应该等于流出的线路电流之和用来检查直流节点方程是否正确。第三换流站运行参数范围。整流侧触发角一般在5到20度之间逆变侧熄弧角一般在15到20度之间超出这个范围说明运行点不合理。第四直流电压水平。直流节点电压应该在额定值附近不远如果某节点电压严重偏离可能是定电压站和定功率站的配置有问题。第五交流母线电压越限检查。换流站接入点的交流电压应在0.95到1.05标幺之间否则可能导致换流器无法正常工作。这些校验项写起来虽然琐碎但确实是防止算例出问题的最有效防线。我每次程序有改动都会重新过一遍省去了很多后面排查异常结果的麻烦。5. 常见问题与排查技巧实录5.1 统一法雅可比矩阵奇异怎么办矩阵奇异的直接表现是解线性方程组时得到NaN或者迭代一步就满天飞。原因是未知量和方程数量不匹配或者方程之间存在线性相关。排查思路是先检查方程数量和未知量数量是否相等用手工统计表把每个未知量对应的真实方程列出来逐一核对。比较隐蔽的问题是控制方程冲突。比如两个换流站同时定直流电压在两端直流系统里还有可能分配电压差在多端直流系统里如果直流网络连接紧密两个定电压站之间可能引起方程冗余。解决方法是把其中一个改成一个定功率或定电流站只保留一个定电压站。另外换流变压器变比和分接头处理的简化也可能导致方程线性相关这时考虑把T固定为常数。初值远离可行域也会造成数值上的“瞬时奇异”。例如整流侧触发角初值给20度实际可行域在5到20度之间虽然还没有越限但换流器方程对控制角的导数可能接近零导致矩阵条件数极差。这种情况可以观察雅可比矩阵对角线元素有没有突然变小快速定位到对应变量。5.2 交替法不收敛和振荡的处理路径交替法不收敛通常表现为几种模式发散型变量一轮比一轮大等幅振荡型两个子系统的结果在两个固定值之间来回跳慢收敛型迭代几十步还在缓慢漂移。前两种优先考虑加松弛因子和低通滤波。发散了把λ降到0.3左右同时限制单次迭代的变量变化幅度。等幅振荡则要检查换流站控制方式的设置是否和潮流方向一致特别是有没有出现整流和逆变角色判定错误。慢收敛型的问题往往出在交流系统与直流系统的耦合强度上。如果换流站接入点是弱交流母线电压对功率变化非常敏感接口电压的小波动会引发直流子系统的功率大变化形成正反馈。处理方式是给接口电压加一阶惯性环节或者先固定换流站交流母线电压一个中间值等直流子系统稳定后再放开。这个方法我试过多次对弱系统场景非常有效。5.3 控制模式切换的边界处理潮流计算过程中换流站控制模式可能会发生切换比如定功率运行时直流电流达到额定上限此时控制模式应该变为定电流定直流电压站电压偏差过大导致电流越限可能切换到定电流或者定功率。模式切换本质上改变了控制方程这在统一法和交替法里都是一个麻烦。我的处理方式是分层迭代策略。外层迭代固定一组控制模式内层迭代求解当前模式下的潮流方程。每一次外层收敛后检查所有换流站的越限情况如果某个换流站越限则更新它的控制模式重新开始内层迭代。这虽然增加了几层循环但逻辑清晰不会出现在一次迭代中边算边切导致的数值振荡。另外要特别注意切换边界在有些算例里存在多个可行解的情况两个模式可能都满足方程这时候需要结合物理合理性判断到底应该落在哪个模式上。5.4 一个容易被忽略的细节直流线路电阻的标幺值最后一个坑也是我经常在帮别人看程序时发现的问题——直流线路电阻标幺值算错。有些人在交流系统标幺值体系下直接用有名电阻除以交流阻抗基准这是不对的。直流线路的电阻标幺值应该用直流功率基准和直流电压基准来归算即R_dc_pu R_dc_ohm / (U_dc_b² / S_b)。如果误用了交流侧基准值直流网络的潮流结果差得不是一点半点而且这种错误特别隐蔽因为直流电压还是能算出一个看似合理的结果。这种问题最有效的排查办法是把直流子系统单独拿出来做一个小规模的纯直流手算测试。直流线路电阻、两端电压都给定算一次功率分布跟程序输出对比。对了再接入交直流联合计算错了先在最小例子上把基准值体系彻底校准。最后的一些体会做了这些年交直流混联系统计算我的切身体会是统一法和交替法没有高下之分只有场景匹配之分。如果你在写研究性代码、需要稳定的收敛性和精确的结果统一法值得花时间啃下来如果你是在现有交流程序基础上延展功能、快速得到工程近似解交替法会更顺手。我自己现在的工作习惯是两类方法都保留先用交替法摸清算例脾性再切换到统一法出最终结果。另外想提醒一点这两种方法看似是潮流计算问题本质上是换流器模型与控制方式的理解深度问题。把换流器方程吃透了无论用哪种方法代码写起来都会顺畅很多。