VSC交直流混合系统潮流计算：快速灵活全纯嵌入法原理与工程实践
1. 项目概述：当潮流计算遇上全纯嵌入法
在电力系统这个庞大而精密的“神经网络”中，潮流计算扮演着“体检医生”的角色。它的核心任务，就是根据给定的发电、负荷和网络拓扑，计算出电网中每个节点的电压幅值、相角，以及每条线路上的功率流动。这听起来像是解一道复杂的多元非线性方程组——没错，它确实是。几十年来，以牛顿-拉夫逊法为代表的迭代算法一直是解决这道题的标准答案，工程师们对它又爱又恨：爱其成熟高效，恨其对初始值“脾气”太大，一旦给的初值离真实解太远，就可能直接“罢工”不收敛。
随着以风电、光伏为代表的新能源大规模并网，以及基于电压源换流器的柔性直流输电技术快速发展，现代电网正演变为复杂的交直流混合系统。这种系统里，交流侧和直流侧通过VSC紧密耦合，VSC本身又具备多种灵活的控制模式。这给传统的牛顿-拉夫逊法带来了新挑战：每次VSC控制模式切换，都意味着潮流方程的结构发生变化，需要重新构建雅可比矩阵并求解，计算复杂且容易因矩阵奇异而失败。更重要的是，在规划或运行方式分析中，我们常常需要从一个“冷启动”状态（比如全网电压设为1.0∠0°）开始计算，牛顿法在这种情况下的收敛性就像开盲盒，缺乏理论保证。
正是在这样的背景下，一种基于复分析理论的“非迭代”方法——全纯嵌入法进入了我们的视野。它不像牛顿法那样“猜”一个初值然后反复迭代逼近，而是将实数域的潮流方程“嵌入”到复平面，构造一个关于复变量的全纯函数。通过求解这个函数在特定点（对应目标运行状态）的值，直接获得潮流解。其最大的理论魅力在于，只要系统存在高电压可行解，HE方法就能保证收敛到这个解，彻底摆脱了对初值的依赖。然而，早期的HE方法计算速度较慢，且每次计算都必须从一个固定的“零注入”状态开始，在需要频繁重算的交直流混合系统场景中，效率瓶颈明显。
我们今天要深入探讨的，正是针对这一痛点提出的“快速灵活全纯嵌入法”（FFHE）。它继承了HE法收敛可靠的理论优点，同时通过引入更灵活的嵌入形式，允许我们使用更接近真实解的中间结果作为新的计算起点。这就好比在解一个复杂迷宫时，传统HE法每次都必须从唯一的入口重新开始走；而FFHE法则允许你在任何一个已经到达的、靠近出口的位置，直接开始新的探索，大大减少了“绕路”的时间。本文将结合VSC交直流混合系统的具体场景，拆解FFHE方法的模型构建、求解框架和实现技巧，分享如何将这一前沿理论转化为稳定高效的工程计算工具。
2. 核心思路：如何为交直流混合系统“定制”FFHE模型
要将FFHE方法成功应用于VSC交直流混合系统，不能简单地将交流系统和直流系统的模型拼凑在一起。我们需要一套量身定制的建模策略，既要准确刻画VSC站这个“桥梁”的物理特性与控制逻辑，又要保持FFHE方法本身“构造性”和“确定性”的数学美。这里的核心挑战在于，直流子系统中的节点类型（如定直流电压母线、定功率母线、下垂控制母线）与交流系统截然不同，且受VSC控制模式直接影响。
2.1 系统稳态模型的统一描述
首先，我们必须清晰地定义系统的三个组成部分：纯交流子系统、直流子系统和VSC换流站。VSC站通常被建模为一个包含变压器、相电抗器和滤波器的等效π型电路。关键在于明确几个连接点：与交流系统连接的公共连接点母线（PCC bus），与换流器连接的交流连接点母线（PAC bus），以及与直流系统连接的直流连接点母线（PDC bus）。功率就在这几个点之间流动，并伴随着换流器损耗。
VSC的控制策略决定了PCC母线和PDC母线的“身份”。例如，当VSC采用定有功功率和定无功功率控制时，其PCC母线在潮流计算中就视为PQ节点；若采用定有功功率和定交流电压控制，则PCC母线视为PV节点。直流侧亦然，定直流电压的PDC母线类似于交流系统的平衡节点，而定功率或下垂控制的PDC母线则需要建立相应的功率-电压方程。这种多样性是传统统一潮流算法雅可比矩阵构造复杂的重要原因。
2.2 交流子系统FFHE模型的构建
对于交流子系统（包含纯交流母线和作为PQ/PV节点的PCC母线），我们可以直接借鉴纯交流系统的FFHE模型。其核心思想是引入一个复变量α，将原始的潮流方程（8）和（9）嵌入为一个关于α的全纯函数方程（13）和（14）。当α=0时，方程有一个我们预先设定的、已知的“胚解”；当α=1时，方程则恢复为原始的潮流方程，此时的解V͂_i(1)就是我们要求的节点电压。
这个模型的灵活性就体现在“胚解”的选择上。传统HE方法强制要求胚解必须是“平启动”（即所有电压初值为1∠0°）。而FFHE方法允许我们指定任意非零复数vi作为V͂_i(0)的初值。这意味着，我们可以把上一次计算得到的、接近真实解的电压值直接作为下一次计算的起点，这是提升效率的关键。
2.3 直流子系统FFHE模型的创新设计
直流子系统的建模是本文的难点和创新点。直流潮流方程本质上是实数域的非线性方程，但为了应用基于复分析的解析延拓技术，嵌入参数α必须保持为复数。我们需要设计一个新的直流FFHE模型，使其满足三个条件：1）在α=0时有已知且可调的胚解；2）在α=1时能恢复原始直流潮流方程；3）模型结构能统一容纳定功率、定电压和下垂控制等多种VSC控制模式，且切换控制参数时无需改变模型结构。
最终给出的模型（16）巧妙地实现了这些目标。对于采用定功率或下垂控制的PDC母线，其模型如（16a）所示。通过设置下垂系数Ki=0，该模型即退化为定功率母线模型；Ki≠0时，则对应下垂控制母线。对于定直流电压母线，模型则简化为（16b）的形式。无论Ki、Vref等参数如何变化，模型的基本数学结构保持不变，这保证了后续求解过程的通用性。同样，我们为直流电压Vdc_i(α)设定一个非零实数初值vdc_i作为胚解。
关键理解：FFHE模型的“构造性”体现在，通过将非线性方程嵌入复平面，我们将其转化为一个关于α的幂级数求解问题。系数矩阵Aac和Adc在给定胚解后是固定不变的，只需要在递归求解开始时构造一次。这与牛顿法每次迭代都需要重新计算并分解雅可比矩阵形成了鲜明对比，后者在控制模式切换时会产生大量额外计算。
3. 求解引擎：递归过程与帕德逼近的加速魔法
建立了AC/DC FFHE模型，下一步就是如何高效地求解出那个幂级数，并最终得到α=1时的数值解。这个过程分为两个核心步骤：递归求解幂级数系数，以及使用帕德逼近加速级数求和。
3.1 幂级数系数的递归求解流程
我们将节点电压的全纯函数V͂_i(α)和Vdc_i(α)在α=0处展开成幂级数，如式（15）和（17）所示。我们的目标就是求出这些系数uiq和udc_iq。
求解过程是递归进行的。首先，第0阶系数就是我们的胚解：ui0 = vi, udc_i0 = vdc_i。然后，将幂级数表达式代入我们构建的FFHE模型方程（13）、（14）和（16）。通过比较方程两边α的同次幂系数，我们可以得到关于第q阶系数uiq和udc_iq的线性方程组（21）和（36）。
这里有一个精妙之处：在求解第q阶系数时，所有低于q阶的系数（即ui0, ui1, ..., ui(q-1)）都已经是已知量，它们构成了线性方程组的右端项Bac和Bdc。而系数矩阵Aac和Adc仅由网络导纳参数和最初设定的胚解vi、vdc_i决定，在整个递归求解过程中保持不变！这意味着，我们只需要在第一步对系数矩阵进行一次LU分解或求逆，后续每一步求解新的高阶系数时，都只是进行一次前代回代或矩阵乘向量的运算，计算量极小。
递归求解的流程可以清晰地表示为：
- 初始化：设定胚解vi, vdc_i，令q=1。
- 构造方程：根据当前阶数q，利用已知的低阶系数，形成线性方程组A * X[q] = B。
- 求解系数：解该线性方程组，得到第q阶系数向量X[q]（包含uiq或udc_iq）。
- 收敛判断：利用已求得的系数，通过帕德逼近估算α=1时的电压值，并计算功率偏差。若偏差小于预设精度ε，则停止；否则，令q=q+1，返回步骤2。
这个过程是确定性的，只要胚解选择得当（位于高电压解支路上），理论上总可以通过增加级数项数来达到任意精度，不存在牛顿法中的发散问题。
3.2 帕德逼近：让级数快速“收敛”的加速器
理论上，我们需要计算无穷多项幂级数的和才能得到精确解，这显然不现实。实际上，我们计算有限项（例如前L项），然后用一个有理分式函数——帕德逼近来近似这个无穷级数在α=1处的值，如式（43）所示。
帕德逼近的精髓在于，它通常比原始幂级数本身具有更好的收敛性和更大的收敛半径。对于许多函数，即使其泰勒级数收敛很慢甚至发散，帕德逼近也能给出很好的近似值。在FFHE中，我们采用了加速的帕德逼近计算技巧，无需显式计算出有理分式的分子分母多项式系数，而是通过计算一个行列式的比值来直接得到[L/M]_i(1)的近似值。这一步是FFHE相比早期HE方法计算效率大幅提升的关键技术之一。
实操心得：在编程实现时，递归求解系数和帕德逼近是两个相对独立的模块。建议将系数矩阵A的分解（或求逆）计算一次并存储，这会极大提升递归步骤的速度。对于帕德逼近，通常取L和M在10-20之间（即计算10-20阶系数）就能达到极高的精度（如10^-12）。过高的阶数对精度提升有限，但会增加计算量。
4. 整体框架：顺序求解与灵活初始化的协同
有了针对AC和DC子系统的FFHE求解器，我们还需要一个顶层的调度框架，将它们和VSC站的计算模块有机结合起来，处理AC/DC之间的数据交换，并应对VSC控制模式切换的挑战。本文采用的是经典的顺序求解框架，但其与FFHE的结合产生了“1+1>2”的效果。
4.1 顺序计算流程与数据交换机制
整个AC/DC FFHE潮流计算的流程如图2所示，是一个典型的“交流-换流站-直流”顺序迭代过程：
- 初始化：读取网络数据，为AC和DC FFHE模型设定胚解vi和vdc_i（首次计算通常使用平启动）。
- 执行交流FFHE计算：根据当前VSC控制模式，将PCC母线作为PQ或PV节点纳入交流子系统，调用AC FFHE求解器，得到所有交流母线（包括PCC母线）的电压V͂_s。
- 换流站计算：根据求得的PCC母线电压V͂_s和已知功率，利用公式（1）-（3）计算PAC母线的电压V͂_c和功率S͂_c，进而通过公式（5）、（6）计算换流器损耗Ploss。
- 数据交换（AC→DC）：根据公式（4），由PAC母线功率和损耗计算注入PDC母线的有功功率Pdc。对于定功率母线，此Pdc即作为已知量Pref。
- 执行直流FFHE计算：根据已知的Pref、Vref、K等参数，调用DC FFHE求解器，求解直流网络，得到所有PDC母线的电压Vdc。
- 数据交换（DC→AC）：对于定电压母线和下垂母线，由上一步求得的Vdc反推其注入功率Pdc，再通过公式（4）更新PAC母线的有功功率P_c(new)。
- 收敛判断：利用更新后的P_c(new)，重新计算PCC母线的注入功率P_s(new)。检查所有PCC母线前后两次迭代的功率偏差最大值是否小于收敛精度ε_ac/dc。若满足，则整体计算结束；否则，将当前轮次计算得到的AC和DC电压解赋值给下一轮的胚解vi和vdc_i，然后返回步骤2，开始新一轮计算。
4.2 VSC控制模式切换的智能处理
VSC的运行受到其容量和母线电压限值的约束。在潮流计算过程中，可能需要根据计算结果切换其控制模式。例如，一个原本控制交流电压的VSC，若计算发现其无功功率输出越限，则需要切换到定无功功率控制模式。
传统的HE方法在控制模式切换后，必须“重置”到平启动状态重新开始整个AC/DC潮流计算，这造成了巨大的计算浪费。而FFHE方法的优势在此刻凸显。在我们的框架中，当检测到需要切换控制模式时：
- 保存当前轮次计算得到的AC和DC电压解。
- 根据越限情况，更新VSC的控制模式（如将PV节点改为PQ节点）。
- 将上一步保存的电压解直接设置为新一轮AC/DC FFHE计算的胚解vi和vdc_i。
- 在新的控制模式下，重新开始从步骤2起的顺序计算。
由于模式切换前后的系统状态通常变化不大，保存的电压解非常接近新模式下的真实解，以此为胚解可以极大地减少FFHE达到收敛所需的级数项数，从而快速得到新模式下的潮流解。
避坑指南：实现控制模式切换逻辑时，关键在于正确更新相关母线的类型和已知量。例如，从Vs-control切换到Qs-control，对应的PCC母线就从PV节点变为PQ节点，其已知量从（Ps, Vs）变为（Ps, Qs），其中Qs应设置为越限的边界值（如上限或下限）。同时，在构建新一轮的AC FFHE模型时，系数矩阵Aac会根据新的母线类型和已知量自动调整，但得益于FFHE模型的结构化特性，这个调整是局部的、高效的。
5. 性能实测：精度、效率与鲁棒性分析
理论再优美，也需要实际的测试数据来验证。我们通过在两个测试系统（一个简单的5交流/3直流系统和一个基于IEEE RTS-96的复杂系统）上进行大量仿真，将所提出的AC/DC FFHE方法与传统的牛顿-拉夫逊法以及基础HE方法进行了全面对比。
5.1 计算精度验证
首先在测试系统A上，对比了所提方法与N-R方法的结果。设定统一的收敛精度（如10^-6），并以N-R法的结果作为基准。如表IV和表V所示，两种方法计算得到的交流母线电压幅值、相角以及直流母线电压完全一致，最大误差出现在VSC站的功率计算上，但相对误差低于0.1%。这充分证明了AC/DC FFHE方法在数学上是严谨的，能够收敛到与N-R法相同的物理可行解。
5.2 计算效率优势
效率对比是重点。我们比较了完成第一轮AC/DC潮流计算时，HE法和FFHE法所需的级数项数和计算时间。如图5和图6所示，得益于加速帕德逼近技术的应用，FFHE法达到相同精度所需的级数项数远少于HE法，因此单轮计算时间显著缩短。
更重要的效率提升体现在多轮计算和控制模式切换场景中。如图7所示，当需要进行多轮迭代以满足整体收敛，或在控制模式切换后需要重算时，FFHE法可以利用上一轮的解作为“默认启动”胚解，而HE法只能被迫回到“平启动”。测试结果表明，在测试系统B的复杂场景下，采用默认启动的FFHE法总计算时间比采用平启动的HE法减少了38.1%，其效率甚至接近或优于N-R法。这是因为N-R法在每次模式切换后也需要重新构建和求解雅可比矩阵，而FFHE法仅需以更接近的解为起点进行新一轮快速递归。
5.3 初始化策略的影响分析
FFHE法虽然对初值选择比HE法灵活，但并非“全局收敛”。它的收敛域是围绕高电压解的一个局部区域。我们通过改变测试系统A中两个PQ母线（Bus 4, Bus 5）的电压初值vi，观察了AC FFHE达到收敛所需级数项数的变化，结果以三维曲面形式展示在图8和图9中。
可以发现：
- 存在最优启动点：当选择的初值恰好等于或非常接近最终的真实解时，FFHE仅需很少的级数项（如4项）即可收敛。
- 对初值仍敏感：当初值偏离真实解较远时，所需的级数项数会急剧增加，甚至可能不收敛。例如，Bus 4的收敛域比Bus 5更“狭窄”，说明系统不同节点对初值的敏感度不同。
- 工程意义重大：在顺序计算框架中，虽然第一轮计算可能从平启动开始（需要较多项数，如12项），但一旦完成一轮计算，得到的中间解已经非常接近最终解。在后续轮次或模式切换后的计算中，以此中间解为“默认启动”，可以稳定地将所需级数项数降至最低水平（如4项），如表VII所示。这正是FFHE在工程应用中效率飞跃的核心原因。
经验总结：在实际编程和应用中，一个稳健的策略是：首次计算采用平启动（vi=1+j0, vdc_i=1）以保证可行性；之后的每一次重算（无论是顺序迭代未收敛，还是控制模式切换），都毫无例外地将前一次计算得到的全网电压解作为新一轮FFHE的胚解。这一简单的策略能最大化发挥FFHE的灵活性和效率优势。
6. 实现细节与工程化考量
将论文中的数学模型转化为可运行的代码，还需要处理许多工程细节。这里分享一些在复现和实践过程中积累的经验。
6.1 系数矩阵的构建与存储
AC FFHE的系数矩阵Aac（式22a）和DC FFHE的系数矩阵Adc（式37a）是方法的核心。它们的元素由网络导纳矩阵和选定的胚解vi、vdc_i计算而来。注意，对于交流系统，矩阵元素是复数，但在构建方程（21）时，我们将其拆分为实部和虚部，形成了实数的增广矩阵。在编程时，建议预先计算好导纳矩阵Y，并根据当前的母线类型（平衡节点、PQ、PV）和胚解，一次性构建好Aac和Adc。由于它们在递归过程中不变，可以对其进行LU分解并存储分解结果，这样在求解每一阶的线性方程组时，计算复杂度将从O(n^3)降为O(n^2)。
6.2 帕德逼近的稳健实现
帕德逼近的数值稳定性至关重要。建议使用成熟的数值线性代数库（如LAPACK）中的函数来处理行列式计算或直接求解线性方程组来获得帕德逼近值。对于接近奇异的情况，可以增加级数项数（L, M）或采用更高阶的逼近。在实际测试中，对于规模在数百节点的系统，取L=M=15通常足以达到双精度下的机器精度。
6.3 收敛判据的设置
在递归求解过程中，我们需要判断何时停止计算更高阶的系数。如式（44）所示，判据是基于帕德逼近得到的电压解计算出的功率不平衡量。通常设置ε_ac和ε_dc在10^-6到10^-9之间。需要注意的是，由于我们是在AC和DC子系统分别求解，因此还需要一个全局的收敛判据ε_ac/dc来检查AC/DC接口功率（即PCC母线功率）的迭代变化，如式（45）。只有当三者都满足时，整个计算才算收敛。
6.4 处理雅可比矩阵奇异或病态的场景
这是传统牛顿法的传统难题，但在FFHE框架下被自然规避了。FFHE不涉及雅可比矩阵的求逆，其递归求解过程是求解一个固定的、良态的系数矩阵。即使系统运行在接近电压稳定极限的“病态”点，牛顿法可能因雅可比矩阵奇异而失败，FFHE理论上仍能通过增加级数项数来获得解。这使其在电压稳定分析、寻找最大输电能力等场景中具有独特优势。
7. 总结与展望
基于快速灵活全纯嵌入法的VSC交直流混合系统潮流计算，代表了一种从根本原理上革新传统迭代算法的思路。它通过将实数域的潮流问题“嵌入”复平面，利用解析函数的优良性质，将非线性方程组的求解转化为一系列线性方程组的递归求解和有理函数逼近问题。这种方法不仅从理论上保证了向高电压可行解的收敛，更通过灵活的胚解选择机制，在需要频繁重算的复杂场景（如顺序迭代、控制模式切换）中，实现了计算效率的质的提升。
从我个人的实现经验来看，将FFHE方法工程化的最大收获在于思维模式的转变。我们不再需要小心翼翼地“猜测”初值、调试迭代参数、处理迭代发散，而是专注于构建正确的数学模型和高效的数据结构。一旦模型建立，求解过程几乎是确定性的、自动化的。当然，这也对编程的严谨性提出了更高要求，因为递归和帕德逼近中的任何小错误都可能导致结果偏差。
展望未来，这一方向仍有丰富的探索空间。一方面，可以进一步研究代码的深度优化和并行化技术，例如将不同阶系数的求解任务并行处理，以应对未来超大规模交直流混合电网的在线分析需求。另一方面，FFHE方法在理论上能够追踪解支路的特性，使其非常适合于静态电压稳定极限的精确求取。如何利用FFHE直接、高效地定位交直流混合系统的功率传输极限点，将是一个极具理论和实用价值的研究课题。对于从事电力系统分析工具开发的工程师和研究者而言，掌握并发展这类基于现代数学理论的非迭代算法，无疑是保持技术领先性的关键。