基于全纯嵌入法的非迭代电力系统最优潮流计算通过求解电力系统最优潮流获得各发电机的有功出力对于电力系统经济运行具有重要意义。全纯嵌入法是一种用于潮流计算的非迭代算法。该文将全纯嵌入法应用于最优潮流的求解，首先推导出约束条件和KKT条件的全纯嵌入形式；然后通过上述两组方程交叉得到求解通式，递归地计算出节点电压和发电机有功出力的表达式；最后根据解析表达式直接得到各发电机的有功出力，进而计算系统的最低发电成本。在IEEE 3节点、4节点、30节点、118节点和新英格兰39节点电力系统上的仿真结果验证了所提方法的可行性和准确性。
我国电力行业的快速发展，对电力系统资源的优化配置提出了更高的要求[1]。在确保系统安全稳定的运行条件下，优化电力系统的调度运行，降低系统的能量损耗，得到发电成本最低的潮流分布，正是最优潮流（Optimal Power Flow, OPF）所要解决的问题之一。因此，最优潮流自诞生之日起便是电力系统中研究的热点问题，对于电力系统的节能环保低碳和经济性运行有着重要的理论和实践意义[2]开云网址·(中国)官方网站。
最优潮流的概念在1962年由J. Carpentier首次提出[3]，为了适应电力系统的飞速发展，最优潮流的优化模型与求解方法也在不断更新。根据潮流的类型不同分为直流最优潮流[4]和交流最优潮流；根据网络的拓扑结构可以分为辐射型网络最优潮流[5]和环型网络最优潮流[6]。当前求解最优潮流的方法主要可以分为三类：①线]。第①类方法将交流潮流约束近似为直流潮流约束，求解直流最优潮流，同时也避免了OPF模型中的非凸问题，但将交流潮流近似为线性问题得到的结果可能不是交流潮流中的可行解。第②类方法主要采用传统的内点法[9]、梯度法[10]、牛顿法[11]等数学上常用的求解算法对非线性优化问题的进行求解，由于这些方法容易陷入局部最优，近年来针对约束条件进行松弛的凸松弛算法（如半正定规划松弛[12]，二阶锥规划松弛[13-15]等）一直是研究的热点，对非凸问题进行松弛，将优化问题松弛为凸问题，从而能够应用内点法等求解算法得到问题的最优解，在精确松弛的前提下，可以保证得到的解为全局最优解，但不同的优化问题对于松弛有不同的精确要求，在松弛不精确时所得到的解可能为系统的不可运行点[2]。第③类方法为近几年发展迅速的智能算法，如粒子群算法[16]、进化算法[17-18]、遗传算法[19]等，这些启发式算法不采用传统的内点法作为求解器，而是通过经验选择行之有效的方法。全纯嵌入法是由西班牙学者A. Trias于2012年提出的。该方法最初被应用在潮流计算中，它避免了许多传统方法（如牛顿-拉夫逊法[20]、快速分解法[21]
[22-26]，但在系统有功最优潮流的研究方面，目前相关的文献研究较少。由于全纯嵌入法存在非迭代、递归性、解析性、不依赖初始点的选取、所得到的解为可运行点等优点，因此该方法在有功最优潮流的计算中有较大的应用潜力。本文将全纯嵌入的方法用于有功最优潮流的求解，提出了基于全纯嵌入法的电力系统有功最优潮流的计算方法。首先构造出电力系统最优潮流约束条件和KKT条件的全纯嵌入形式；然后推导出节点电压和发电机有功出力解析表达式幂级数系数的递推方程组；最后通过上述两组方程组交叉递归求得解析表达式中的系数，得到优化后系统发电机的出力以及发电成本。本文通过IEEE 3节点、4节点、30节点、118节点和新英格兰39节点系统有功最优潮流的算例仿真，验证了该方法的可行性及准确性。1 最优潮流的模型
在传统的最优潮流问题中在针对经济调度问题时一般不考虑PV节点，但发电机还是有自身的电压控制功能，所以为方便推导及理解，本文将发电机节点考虑为PV节点。
同时为了便于理解且不失一般性，本节以3节点系统模型为例，包含1个平衡节点（节点1）、一个PV节点（节点2）和一个PQ节点（节点3），如图1所示。其中方程中所有的大写字母表示常数，小写字母表示变量。同时，附录给出了由
p2分别为节点1和节点2处的发电机有功出力，发电成本表示为发电机有功出力的二次函数；Cij为成本系数，视为常数项。平衡节点和PV节点的有功出力p1和p2为控制变量，而PQ节点的有功输出P3和无功输出Q3为已知量。1.2 约束条件以直角坐标系下的潮流方程组以及全网的有功平衡方程作为约束条件。
i的电压；ViSP为固定值，表示平衡节点或PV节点的电压值；q1和q2分别为节点1和节点2处的发电机无功出力；Gij为节点i到节点j之间的互电导；Gii为节点i的对地电导；Bij为节点i到节点j之间的互电纳；Bii为节点i的对地电纳；Yi,shunt为节点i对地的导纳；Yij,trans为节点i和节点j之间的互导纳；Si为节点i处的视在功率，如图1所示；上标*表示共轭算子。式（7）中左端第一项为网络损耗，分为节点对地的导纳产生的损耗和节点之间互导纳产生的损耗。同时目标函数里的p1和p2应是发电机的有功出力，而约束里的P是指对应节点的有功注入，即发电机和负荷的共同效果。图1中所示的3节点系统中节点1、2只有发电机，节点3只有负荷，故不需要考虑此类问题。目前为止本文只考虑了等式约束，没有涉及不等式约束，在数学上求取考虑不等式约束条件目标函数最小值的方法已有很多[27]，同时文献[28]解决了全纯嵌入法在潮流计算中考虑控制上下限的问题，上述解决不等式约束的方法均可以和本文所提方法相结合。本文侧重于首度探索全纯嵌入法在有功优化方面的可行性，因此考虑不等式约束的有功优化将可以作为后续研究内容。
虽然式（7）中左端第一项括号中的线路损耗没有表示出与发电机有功出力的相关性，但由于约束条件式（2）～式（5）的存在使得发电机有功出力的改变会改变节点处的电压，所以网络的有功损耗既是所有节点电压的函数，也是所有节点有功输入的函数
s0的邻域U内复可微，则称函数x(s)是复数域U上的全纯函数。全纯函数能够在该复数域内展开成幂级数的形式，有（10）
x)=0按照一定的规则将其变量x（如节点电压v和发电机的有功出力p）构建成全纯函数x(s)，可以得到复合函数g(x(s))=0，其中只要计算出式（10）中各个变量的幂级数系数x[n]，则能得到变量的解析表达式。本文所有幂级数系数均使用方括号表示，函数使用小括号表示。全纯嵌入函数的构造可按照如下规则：1）当s
x(0))=0，如在潮流计算中可以表示为发电机在没有有功输入、系统没有负荷时的状态。2）当s=1时表示所要计算的目标解，同时也能够满足g
g(x(s))=0左右两侧s同阶系数相同构造的递推关系求得[7]。将s=0代入式（10）中，由上述构造规则1）可知，得到每个变量的第0阶幂级数x
s=1代入式（10）中，根据构造规则2）可知，所有幂级数系数之和即为变量的目标解，当求得的幂级数阶数越高，目标解就越精确，越接近变量的实际值。2.2 约束条件的全纯嵌入形式根据约束条件式（2）～式（7）与全纯嵌入的基本构造规则，节点电压的初始值设为v
s=1时，上述方程则变为前文中的约束条件。所以本文提出的方法不依赖于系统实际初始点的选取，与系统实际的初始点无关。将式（10）代入约束方程的全纯嵌入方程中，采用附录中的方式，将电压的实部e和虚部f
s在不同阶数时的电压的实部与虚部表达式（12）。式中开云网址·(中国)官方网站，Re(·)表示取实部；Im(·)表示取虚部；
n]可以通过附录中PV节点的电压幅值恒定公式推导得到。（13）将表达式右端项对应PV节点的节点2的无功变量
n]移到等式左侧，将左侧节点2的电压实部e2[n]移到等式的右侧，可以得到节点电压幂级数系数的递推公式为由式（14）可以看出，等式右侧均为低阶幂级数系数或能够通过式（13）计算得到的系数，所以能够通过式（14）逐渐递推出节点电压高阶的幂级数系数。
[0]，fi[0]），后续的系数（ei[n]，fi[n]）可以通过式（13）、式（14）进行递推求解，直至达到所设定的计算精度，节点电压的解析表达式可以由式（15）表示。将
根据拉格朗日方程式（8）、式（9）与全纯嵌入的基本数学规定，设置系统初始运行条件下发电机没有有功出力，将复变量s
(0)=0、λ(0)=0；同时在s=1时，上述方程则变为KKT条件的原始方程。若发电机的有功出力初始值不同，只需要改变全纯嵌入的形式，再次体现了本文所提方法不依赖于系统实际初始点的选取。将式（10）代入式（16）中，根据附录中介绍的方法，可以得到发电机的输出功率pi
[0]，λ[0]），后续系数（pi[n]，λ[n]）根据式（17）递推得到，直至达到所设定的计算精度。发电机的有功出力解析表达式为（18）同理，递推的过程就是在逐步更新发电机有功出力值，使结果逐渐接近有功出力实际值的过程。
将潮流计算的修正方程式代入式（19）中，能够得到采用转置雅克比矩阵法求解网损微增率的计算公式为（20）
本文所提的基于全纯嵌入的最优潮流求解流程如图2所示，其中KKT条件与约束条件的递推计算交替进行。首先根据网络的数据构建约束条件递推方程（14）与KKT条件递推方程（17）左侧矩阵，为计算节点电压和发电机有功出力的幂级数展开系数做准备。给定系统的初始状态，根据前文的推导过程开云网址·(中国)官方网站，令所有节点电压实部e
fi[0]=0，发电机有功输入pi[0]=0，此时由于系统中没有潮流流过，所以初始状态下系统也没有有功损耗。图2 应用全纯嵌入法交叉求解有功最优潮流流程
的第一阶幂级数系数；再由得到的所有变量的第一阶系数应用式（17）计算发电机出力的第二阶系数，然后再应用式（13）、式（14）得到节点电压的第二阶幂级数系数，循环递推。递推的过程就是不断更新发电机出力和节点电压的过程，当计算的递推阶数越高，由式（15）、式（18）计算得到的变量实际值越接近目标解。当变量实际值变化达到所设置的计算精度要求时，停止递推，并计算实际的发电成本。
由于在进行递推求解的过程中，每一步都包含KKT条件和约束条件的计算，所以得到的解满足约束条件，即为系统的可运行点，同时由全纯嵌入的规则可知，该方法不依赖于初始点的选取。应用全纯嵌入法求解有功最优潮流伪代码见表1。
算法：基于全纯嵌入法的优化潮流算法输入：系统参数输出：发电成本cost 开始 初始化，准备式（14）、式（17）中左侧矩阵 给定初始条件ei[0], fi[0], pi[0] While(网损变化值＜设置精度)应用式（17）计算得到pi[n], λ[n]应用式（13）、式（14）计算得到ei[n], fi[n]应用式（15）、式（18）更新节点电压与有功出力的实际值更新雅克比矩阵J应用式（19）网损微增率应用式（7）计算网络损耗n = n+1end 应用式（1）计算发电成本 结束
3 仿线节中推导的全纯嵌入方法在有功最优潮流中的可行性与准确性，本文采用IEEE 3节点、4节点、30节点、新英格兰39节点电力系统和IEEE 118电力系统算例进行仿真计算，同时与内点法的计算结果做比较，仿真软件采用Matlab 2020a，在配置为Intel Core i5-9400 CPU @ 2.90 GHz、8 GB内存的计算机上运行，内点法的优化潮流结果应用Matpower4.1工具箱计算得到。
为了验证本文所提的基于全纯嵌入的潮流计算方法计算有功最优潮流的可行性，全纯嵌入法在IEEE 3节点和4节点电力系统中的发电机有功出力和节点电压的仿线分别表示全纯嵌入法在IEEE 30节点和新英格兰39节点电力系统的仿线表示全纯嵌入法在四个系统的仿真过程收敛情况。有功网络损耗变化的精度设置为1×10
图4 新英格兰39节点电力系统仿线 Simulation results of 39-bus power system in New England
x*)/x*计算，其中x为全纯嵌入法的计算结果，x*表示内点法的计算结果。由表4和表5可知，全纯嵌入法计算前两个系统的相对误差均为10-6数量级，计算结果满足精度要求，误差来源为递推的阶数，每次递推的过程都是使结果逐渐逼近精确解的过程。误差随着系统规模的增大而增大，在新英格兰39节点系统中误差较大，是因为在该系统中没有考虑不等式约束条件，使得有功出力较为理想化，从而使发电成本和有功网络损耗都偏小。所以接下来要对考虑不等式约束的全纯嵌入法进行深一步的研究。但本文通过仿真实例证实了该方法的可行性，为电力系统有功优化潮流计算提供了新思路。表4 发电成本Tab.4 Power generating costIEEE节点系统全纯嵌入法/（$/h）内点法/（$/h）相对误差 IEEE 37 861.407 861.391.27×10-6 IEEE 49 936.139 936.103.02×10-6 IEEE 303.1193.3025.54×10-2 IEEE 3939.12245.7291.44×10-1表5 有功网络损耗
由2.4节可知，全纯嵌入法计算有功最优潮流的时间与递推的阶数相关。图6展示了全纯嵌入法在对四个系统进行仿真时递推阶数与计算时长的关系，过多的递推次数在增加计算精度的同时必定会耗费更多的时间，因此在满足精度的情况下选择合适的递推次数即可，在合适的递推精度下全纯嵌入法有着较好的计算效率。在本文所设置精度下全纯嵌入法与内点法的计算时长见表7。
本文所提方法的优越性在于构造全纯嵌入函数时所选取的初始点与系统的当前运行状态无关，且能够在递推计算时考虑以潮流计算为约束条件，进而确保计算结果为可运行解，同时约束条件和KKT条件的全纯嵌入方程左侧矩阵固定，可以离线推导得到，所以能够减少在线节点系统仿真
上述结果表明，所提的基于全纯嵌入的最优潮流计算方法在大规模系统中具有可行性，展现了良好的准确性和计算效率。4 结论本文建立了以全纯嵌入法为基础的有功最优潮流的计算方法，从有功最优潮流的基本数学模型入手，在约束条件和KKT条件中嵌入复变量，构造约束条件和KKT条件的全纯嵌入方程组，再通过递推求解的方式得到节点电压与发电机有功出力的幂级数系数，从而得到二者的解析解，最终得到优化后的发电成本。在IEEE 3节点、4节点、30节点、118节点和新英格兰39节点系统中的仿真结果验证了所提方法的可行性与准确性。主要结论如下：1）应用全纯嵌入法能够得出节点电压和发电机有功出力的解析表达式，通过KKT条件和约束条件两部分交叉进行非迭代的递推计算，求得解析表达式幂级数的系数。
3）采用全纯嵌入法可以将有功最优潮流的求解变为递推求解系统节点电压和发电机有功出力的过程，逐步更新节点电压和有功出力，计算结果满足精度要求，同时证明了本文所提方法的可行性。
本文所提出的基于全纯嵌入的方法为大型电力系统最优潮流求解提供了新的思路，具有一定的参考价值。但本文在构建全纯嵌入函数时，发电机节点不含有负荷，如果考虑在发电机节点侧的恒定负荷，只需要根据全纯嵌入的规则将潮流约束的全纯嵌入方程组稍加修改。同时本文没有考虑不等式约束条件，在文献[25]中提供了全纯嵌入法考虑不等式约束的潮流计算应用，所以考虑不等式约束的基于全纯嵌入法的有功最优潮流计算值得进一步深入研究。
分别为PV节点和PQ节点的集合。1）平衡节点对于平衡节点，由于其电压固定，根据式（A1）中平衡节点公式可得其节点电压全纯嵌入的表达式为
将等式左右两侧所有变量写成式（10）表示的幂级数形式，再令等式两侧s相同阶数的系数相同，由初始条件p
