最优潮流问题的应用与计算.doc

学校代码： 10128学 号： 200621202089 本科毕业设计说明书（题 目：最优潮流问题的应用与计算学生姓名：史界铭学 院：电力学院系 别：电力工程系专 业：电气工程及其自动化班 级：电气063班指导教师：韩如月讲师/硕士 二 一年 六 月 摘要电力系统的运行必须在任何情况下，都应该尽可能的保证电力系统运行的可靠性，保证电能的良好质量，保证运行的最大经济性。因此，做好规划设计就十分重要。在规划领域，需要进行潮流分析验证规划方案的合理性。潮流是确定电力网络运行状态的基本因素，潮流问题是研究电力系统稳态问题的基础和前提。无论进行电力系统规划设计，还是对各种运行状态的研究分析，都须进行潮流计算。同时，为了求使电力系统的某一指标达到最优值的潮流分布，法国学者在Carpentier在20世纪60年代提出了最优潮流的概念。电力系统最优潮流，简称OPF（Optimal Power Flow）。最优潮流是一个复杂的非线性规划问题，要求在满足特定的电力系统运行和安全约束条件下，通过调整系统中可利用控制手段实现预定目标最优的系统稳定运行状态。发展到今天，最优潮流应用领域以十分广泛，针对不同的应用，最优潮流模型可以选择不同的控制变量、状态变量集合，不同的目标函数以及不同的约束条件。本次毕业设计将以媒耗和网损为目标函数，进行研究。关键词：电力系统潮流计算；牛顿拉夫逊法潮流计算； MATLAB；最优潮流；煤耗；网损abstractPower system operation must be under any circumstances, should as far as possible the operation reliability of power system, ensure the quality, ensure good electrical maximum efficiency. Therefore, the plan design is very important.In the planning area, the need for current analysis verify the rationality of the plan. Tidal power network operation is to determine the basic factors, current state of research power system is the basis and premise of steady-state problems. Whether for power system planning and design, or for the operational status of research and analysis, are calculated to tide. At the same time, in order to make the power system for an index to achieve optimal value of the trend, French scholars in Carpentier distribution in the 1960s put forward the concept of optimal fashion.Power system Optimal fashion, referred to as Optimal Power Flow OPF (.). The optimal fashion is a complex nonlinear programming problem, in particular requirements of power system operation and safety conditions, and by adjusting control system can realize the optimal objectives of stable operation system. Today, the development of optimal power flow in a very wide application field, according to the different application, the optimal power flow model can choose different control variables, the state variables set, different target function and constraints. The graduation design in coal consumption and net loss for target function, were studied.Keywords: the power flow calculation system, Thenewton-raphson method computation, MATLAB, OPF, Coal, Network loss目 录第一章 电力系统常规潮流计算61.1潮流计算简介61.2潮流计算的意义及其发展61.3潮流计算的数学模型潮流方程71.4电力系统节点分类81.5潮流计算的约束条件91.6牛顿-拉夫逊法101.7电力系统潮流计算算例：17第二章 电力系统最优潮流计算242.1最优潮流概述242.2算例的最优潮流计算：252.3最优潮流与常规潮流的比较：37第三章 最优潮流的应用403.1最优潮流的应用40总结43参考文献44总 结45参考文献46谢 辞47第一章 电力系统常规潮流计算1.1潮流计算简介电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和基础的计算。 电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。 利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。对潮流计算的要求可以归纳为下面几点：（1）计算方法的可靠性或收敛性；（2）对计算机内存量的要求；（3）计算速度；（4）计算的方便性和灵活性。电力系统潮流计算问题在数学上是一组多元非线性方程式求解问题，其解法都离不开迭代。因此，对潮流计算方法，首先要求它能可靠地收敛，并给出正确答案。由于电力系统结构及参数的一些特点，并且随着电力系统不断扩大，潮流计算的方程式阶数也越来越高，对这样的方程式并不是任何数学方法都能保证给出正确答案的。这种情况成为促使电力系统计算人员不断寻求新的更可靠方法的重要因素。1.2潮流计算的意义及其发展电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采用牛顿-拉夫逊法。在运行方式管理中，潮流是确定电网运行方式的基本出发点；在规划领域，需要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了多个在预想操作情况下电网的潮流分布以校验运行可靠性。在电力系统调度运行的多个领域都涉及到电网潮流计算。潮流是确定电力网络运行状态的基本因素，潮流问题是研究电力系统稳态问题的基础和前提。在用数字计算机解电力系统潮流问题的开始阶段，普遍采取以节点导纳矩阵为基础的逐次代入法。这个方法的原理比较简单，要求的数字计算机内存量比较下，适应50年代电子计算机制造水平和当时电力系统理论水平。但它的收敛性较差，当系统规模变大时，迭代次数急剧上升，在计算中往往出现迭代不收敛的情况。这就迫使电力系统计算人员转向以阻抗矩阵为基础的逐次代入法。阻抗法改善了系统潮流计算问题的收敛性，解决了导纳法无法求解的一些系统的潮流计算，在60年代获得了广泛的应用。阻抗法的主要缺点是占用计算机内存大，每次迭代的计算量大。当系统不断扩大时，这些缺点就更加突出。为了克服阻抗法在内存和速度方面的缺点，60年代中期发展了以阻抗矩阵为基础的分块阻抗法。这个方法把一个大系统分割为几个小的地区系统，在计算机内只需要存储各个地区系统的阻抗矩阵及它们之间联络线的阻抗，这样不仅大幅度地节省了内存容量，同时也提高了计算速度。克服阻抗法缺点的另一途径是采用牛顿-拉夫逊法。这是数学中解决非线性方程式的典型方法，有较好的收敛性。在解决电力系统潮流计算问题时，是以导纳矩阵为基础的，因此，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿法潮流程序的效率。自从60年代中期，在牛顿法中利用了最佳顺序消去法以后，牛顿法在收敛性。内存要求。速度方面都超过了阻抗法，成为60年代末期以后广泛采用的优秀方法。1.3潮流计算的数学模型潮流方程它具有如下特点：（1）它是一组代数方程，因而表征的是电力系统的稳定运行特性。（2）它是一组非线性方程，因而只能用迭代方法求其数值解。（3）由于方程中的电压和导纳既可以表为直角坐标，又可表为极坐标，因而潮流方程有多种表达形式-极坐标形式，直角坐标形式。取 ， ，得到潮流方程的直角坐标形式：（1-1）取， ，得到潮流方程的极坐标形式： （1-2）不同坐标形式的潮流方程适用于不同的迭代解法。1.4电力系统节点分类用一般的电路理论求解网络方程，目的是给出电压源(或电流源)研究网络内的电流(或电压)分布，作为基础的方程式，一般用线性代数方程式表示。然而在电力系统中，给出发电机或负荷连接母线上电压或电流(都是向量)的情况是很少的，一般是给出发电机母线上发电机的有功功率(P)和母线电压的幅值(U)，给出负荷母线上负荷消耗的有功功率(P)和无功功率(Q)。主要目的是由这些已知量去求电力系统内的各种电气量。所以，根据电力系统中各节点性质的不同，很自然地把节点分成三类： PQ节点对这一类点，事先给定的是节点功率(P，Q)，待求的未知量是节点电压向量(U，)，所以叫PQ节点。通常变电所母线都是PQ节点，当某些发电机的输出功率P。Q给定时，也作为PQ节点。PQ节点上的发电机称之为PQ机(或PQ给定型发电机)。在潮流计算中，系统大部分节点属于PQ节点。 PU节点这类节点给出的参数是该节点的有功功率P及电压幅值U，待求量为该节点的无功功率Q及电压向量的相角。这类节点在运行中往往要有一定可调节的无功电源。用以维持给定的电压值。通常选择有一定无功功率储备的发电机母线或者变电所有无功补偿设备的母线做PU节点处理。PU节点上的发电机称为PU机(或PU给定型发电机) 平衡节点在潮流计算中，这类节点一般只设一个。对该节点，给定其电压值，并在计算中取该节点电压向量的方向作为参考轴，相当于给定该点电压向量的角度为零。也就是说，对平衡节点给定的运行参数是U和，因此有城为U节点，而待求量是该节点的P。Q，整个系统的功率平衡由这一节点承担。关于平衡节点的选择，一般选择系统中担任调频调压的某一发电厂(或发电机)，有时也可能按其他原则选择，例如，为提高计算的收敛性。可以选择出线数多或者靠近电网中心的发电厂母线作平衡节点。以上三类节点4个运行参数P。Q。U。中，已知量都是两个，待求量也是两个，只是类型不同而已。1.5潮流计算的约束条件电力系统运行必须满足一定技术和经济上的要求。这些要求够成了潮流问题中某些变量的约束条件，常用的约束条件如下：1. 节点电压应满足 (1-3)从保证电能质量和供电安全的要求来看，电力系统的所有电气设备都必须运行在额定电压附近。PU节点电压幅值必须按上述条件给定。因此，这一约束条件对PQ节点而言。2. 节点的有功功率和无功功率应满足 (1-4) PQ节点的有功功率和无功功率，以及PU节点的有功功率，在给定是就必须满足上述条件，因此，对平衡节点的P和Q以及PU节点的Q应按上述条件进行检验。3. 节点之间电压的相位差应满足 (1-5) 为了保证系统运行的稳定性，要求某些输电线路两端的电压相位不超过一定的数值。这一约束的主要意义就在于此。 因此，潮流计算可以归结为求解一组非线性方程组，并使其解答满足一定的约束条件。常用的方法是迭代法和牛顿法，在计算过程中，或得出结果之后用约束条件进行检验。如果不能满足要求，则应修改某些变量的给定值，甚至修改系统的运行方式，重新进行计算。1.6牛顿-拉夫逊法电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采用牛顿-拉夫逊法。以下讨论的是用直角坐标形式的牛顿拉夫逊法潮流的求解过程。当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量由于平衡节点的电压向量是给定的，因此待求两共需要2(n-1)个方程式。事实上，除了平衡节点的功率方程式在迭代过程中没有约束作用以外，其余每个节点都可以列出两个方程式。对PQ节点来说，是给定的，因而可以写出 （1-6）对PV节点来说，给定量是，因此可以列出 （1-7）求解过程大致可以分为以下步骤：（1）形成节点导纳矩阵（2）将各节点电压设初值U，（3）将节点初值代入相关求式，求出修正方程式的常数项向量（4）将节点电压初值代入求式，求出雅可比矩阵元素（5）求解修正方程，求修正向量（6）求取节点电压的新值（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步（8）计算支路功率分布，PV节点无功功率和平衡节点柱入功率。以直角坐标系形式表示. 迭代推算式 采用直角坐标时,节点电压相量及复数导纳可表示为: (1-8)将以上二关系式代入上式中,展开并分开实部和虚部;假定系统中的第1,2,m号为PQ节点,第m+1,m+2,n-1为PV节点,根据节点性质的不同,得到如下迭代推算式: 对于PQ节点 (1-9)对于PV节点 (1-10)对于平衡节点 平衡节点只设一个,电压为已知,不参见迭代,其电压为: (1-11). 修正方程式(2-3-5)和(2-3-6)两组迭代式工包括2(n-1)个方程.选定电压初值及变量修正量符号之后代入式(2-3-5)和(2-3-6),并将其按泰勒级数展开,略去二次方程及以后各项,得到修正方程如下: (1-12) (1-13).雅可比矩阵各元素的算式式(3-2-8)中, 雅可比矩阵中的各元素可通过对式(3-2-4)和(3-2-5)进行偏导而求得.当时, 雅可比矩阵中非对角元素为 (1-14)当时,雅可比矩阵中对角元素为: (1-15)由式(3-2-9和(3-2-10)看出,雅可比矩阵的特点:矩阵中各元素是节点电压的函数,在迭代过程中,这些元素随着节点电压的变化而变化;导纳矩阵中的某些非对角元素为零时,雅可比矩阵中对应的元素也是为零.若,则必有;雅可比矩阵不是对称矩阵;雅可比矩阵各元素的表示如下: （1-16） 牛顿拉夫逊法的程序框图1.7电力系统潮流计算算例：系统如下图所示计算节点导纳矩阵：节点导纳矩阵： 1.37874 -0.62402 -0.75471 0 0 -j6.29166 j3.90015 j2.64150 0 0 -0.62402 1.45390 -0.82987 0 0 j3.90015 -j66.98082 j3.11203 j63.19206 0 -0.75471 -0.92987 1.58459 0 0 j2.64150 j3.11203 -j35.73786 -j66.66667 0 0 0 0 0 0 0 j63.49206 0 -j66.66667 0 0 0 0 0 0 0 0 j31.74603 0 -j33.33333MATLAB程序：clearG(1,1)=1.37874;B(1,1)=-6.29166;G(1,2)=-0.62402;B(1,2)=3.90015;G(1,3)=-0.754 71;B(1,3)=2.64150;G(1,4)=0;B(1,4)=0;G(1,5)=0;B(1,5)=0;G(2,1)=-0.62402;B(2,1)=3.90015;G(2,2)=1.45390;B(2,2)=-66.98082;G(2,3)=-0.82987;B(2,3)=3.11203;G(2,4)=0;B(2,4)=63.49206;G(2,5)=0;B(2,5)=0;G(3,1)=-0.75471;B(3,1)=2.64150;G(3,2)=-0.82987;B(3,2)=3.11203;G(3,3)=1.58459;B(3,3)=-35.73786;G(3,4)=0;B(3,4)=0;G(3,5)=0;B(3,5)=31.74603;G(4,1)=0;B(4,1)=0;G(4,2)=0;B(4,2)=63.49206;G(4,3)=0;B(4,3)=0;G(4,4)=0;B(4,4)=-66.66667;G(4,5)=0;B(4,5)=0;G(5,1)=0;B(5,1)=0;G(5,2)=0;B(5,2)=0;G(5,3)=0;B(5,3)=31.74603;G(5,4)=0;B(5,4)=0;G(5,5)=0;B(5,5)=-33.33333;Y=G+j*Bdelt(1)=0;u(1)=1.0;delt(2)=0;u(2)=1.0;delt(3)=0;u(3)=1.0;delt(4)=0;p(1)=-1.69;p(2)=-2.2;p(3)=-4.5;p(4)=5.5;q(1)=-0.8;q(2)=-1;q(3)=-1.3;k=0;precision=1;N1=3;%the N1 is the amount of the PQ busN2=1;while precision>0.00001u(4)=1.05; delt(5)=0; u(5)=1.05; for m=1:N1+N2 for n=1:N1+N2+1 pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n); end pp(m)=p(m)-sum(pt);endfor m=1:N1 for n=1:N1+N2+1 qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n); end qq(m)=q(m)-sum(qt);endfor m=1:N1+N2 for n=1:N1+N2+1 h0(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n); end H(m,m)=sum(h0)-u(m)2*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n); H(m,m)=sum(h0)-u(m)2*(G(m,m)*sin(delt(m)-delt(m)-B(m,m)*cos(delt(m)-delt(m); JJ(2*m-1,2*m-1)=H(m,m);end%for m=1:N1 for n=1:N1+N2+1 n0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n); j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n); l0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n); end N(m,m)=sum(n0)-2*u(m)2*G(m,m)+u(m)2*(G(m,m)*cos(delt(m)-delt(m)+B(m,m)*sin(delt(m)-delt(m); J(m,m)=sum(j0)+u(m)2*(G(m,m)*cos(delt(m)-delt(m)+B(m,m)*sin(delt(m)-delt(m); L(m,m)=sum(l0)+2*u(m)2*B(m,m)+u(m)2*(G(m,m)*sin(delt(m)-delt(m)-B(m,m)*cos(delt(m)-delt(m); JJ(2*m-1,2*m)=N(m,m); JJ(2*m,2*m-1)=J(m,m); JJ(2*m,2*m)=L(m,m);end for m=1:N1+N2 for n=1:N1+N2 if m=n else H(m,n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n); JJ(2*m-1,2*n-1)=H(m,n); end end end for m=1:N1+N2 for n=1:N1 if m=n else N(m,n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n); JJ(2*m-1,2*n)=N(m,n); end end endfor m=1:N1 for n=1:N1+N2 if m=n else J(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n); JJ(2*m,2*n-1)=J(m,n); end endend for m=1:N1 for n=1:N1 if m=n else L(m,n)=-H(m,n); JJ(2*m,2*n)=L(m,n); end endend%for m=1:N1+N2 PP(2*m-1)=pp(m); end for m=1:N1 PP(2*m)=qq(m); end uu=-inv(JJ)*PP' precision=max(abs(uu); for n=1:N1 delt(n)=delt(n)+uu(2*n-1); u(n)=u(n)+uu(2*n); end n=N1+N2; delt(n)=delt(n)+uu(n); delt(n)=delt(n)+uu(7); k=k+1;endk-1,delt',u'% the following program is used to calculate the S5 for n=1:5 U(n)=u(n)*(cos(delt(n)+j*sin(delt(n);endfor m=1:5 I(m)=conj(Y(5,m)*conj(U(m);endS5=U(5)*sum(I)for m=1:5 I(m)=conj(Y(4,m)*conj(U(m);endS4=U(4)*sum(I)% the following program is used to calculate the Smny0(1,2)=j*0.25;y0(2,1)=j*0.25;y0(1,3)=0;y0(3,1)=0;y0(1,4)=0;y0(1,5)=0;y0(2,3)=j*0.25;y0(3,2)=j*0.25;y0(2,4)=-j*0.3308;y0(4,2)=j*0.315;y0(2,2)=0;y0(3,5)=-j*0.6615;y0(5,3)=j*0.63;y0(2,4)=0;y0(2,5)=0;y0(3,3)=0;y0(3,4)=0;y0(4,1)=0;y0(4,3)=0;y0(4,4)=0;y0(4,5)=0;y0(5,1)=0;y0(5,2)=0;y0(5,4)=0;y0(5,5)=0;for m=1:5 for n=1:5 S(m,n)=U(m)*(conj(U(m)-conj(U(n)*conj(-Y(m,n)+u(m)2*conj(y0(m,n); endendSPG1=real(S4);PG2=real(S5);f=(50.4395*PG12+200.4335*PG1+1200.6485)+(200.55*PG22+500.746*PG2+1857.201)计算结果：发电机有功及无功出力发电机序号有功出力无功出力45.49992.071153.23612.5965各节点电压向量母线序号电压幅值电压相角10.8373-0.095621.07460.344831.0292-0.094541.05000.421751.05000支路号首末端母线号支路有功功率PijPji112-1.46621.5845213-0.13380.15693231.4156-1.2774424-55535-2.57942.5794支路有功功率第二章 电力系统最优潮流计算2.1最优潮流概述电力系统最优潮流，简称OPF（Optimal Power Flow）。最优潮流是一个复杂的非线性规划问题，要求在满足特定的电力系统运行和安全约束条件下，通过调整系统中可利用控制手段实现预定目标最优的系统稳定运行状态。发展到今天，最优潮流应用领域以十分广泛，针对不同的应用，最优潮流模型可以选择不同的控制变量、状态变量集合，不同的目标函数以及不同的约束条件。2.2最优潮流模型最优潮流模型是在以下前提条件下提出的：(1)各火电核电投入运行的机组己知(不解决机组开停问题)。(2)各水电机组的出力己定(由水库经济调度确定)。(3)电力网络结构确定(不受接线方式影响，不考虑网络重构问题)。 最优潮流在数学上是一个带约束条件的优化问题，其中主要构成包括变量集合，约束条件和目标函数。最优潮流模型中，变量主要分为两大类。一类是控制变量，是可以控制的自变量，通常包括各机组的有功出力，各发电机或同步补偿机无功出力；移相器抽头位置、可调变压器抽头位置、并联电抗器或电容器容量；在某些紧急状况下，水电机组快速启动、某些负荷的卸载也可以作为控制手段。另一类是状态变量，是控制变量的因变量，通常包括各节点电压和各支路功率等。最优潮流考虑的系统约束条件有：(1)各节点有功功率和无功功率平衡约束(2) 各发电机有功出力上下界约束。（3）各发电机/同步补偿机无功出力上下界约束。(4)并联电抗器/电容器容量约束。(5) 移相器抽头位置约束。（6）可调变压器抽头位置约束。(7)各节点电压幅值上下界约束。(8) 各支路传输功率约束。从数学观点来看，以上约束条件中，（1）为等式约束，其余为不等式约束；从约束的物理特性来看，（2）-（6）为控制变量约束，（7）（8）称为状态变量约束。最优潮流有各式各样的目标函数，最常用的形式有以下两种：（1） 系统运行成本最小。该目标函数一般表示为火电机组燃料费用最小，不考虑机组的启动、停机等费用。其中机组成本耗费曲线是模型的关键问题，它不仅影响解的最优性，还制约求解方法的选取。（2） 有功传输损耗最小。无功优化潮流通常以有功传输损耗最小为目标函数，它在减少系统有功损耗的同时，还能改善电压质量。2.2算例的最优潮流计算：系统如图所示示：线路传输功率边界支路号首末端母线号线路传输功率边界11222130.65323242465355发电机数据发电机序号母线号出力上界出力下界燃料耗费参数有功无功有功无功二次系数一次系数常数14831-350.4395200.43351200.648525851-2.1200.55500.7461857.201状态变量10个：X=(Ø1 V1 Ø2 V2 Ø3 V3 Ø4 V4 Ø5 V5 ) (2-1)系统中两台发电机为控制变量：U=（PG1 PG2 QG1 QG2） (2-2)系统总变量为14个以煤耗为目标函数的最优潮流：目标函数： (2-3)约束条件: 非发电机节点 (i=1,2,3) (2-4)发电机节点 （i=4,5） (2-5)式中：表示第k台发电机接在节点i上，不等式约束条件共有14个，分别为 （i=1,2） (2-6)