最优潮流编程 节点导纳矩阵编程.pdf

《最优潮流编程 节点导纳矩阵编程.pdf》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、一、实验项目名称最优潮流计算实验二、实验目的与要求：目的：电力系统分析的潮流计算是电力系统分析的一个重要的部分。通过对电力系统潮流分布的分析和计算，可进一步对系统运行的安全性，经济性进行分析、评估，提出改进措施。电力系统潮流的计算和分析是电力系统运行和规划工作的基础。潮流计算是指对电力系统正常运行状况的分析和计算。通常需要已知系统参数和条件，给定一些初始条件，从而计算出系统运行的电压和功率等。本实验采用简化梯度法最优潮流进行电力系统分析的潮流计算程序的编制与调试，获得电力系统中各节点电压，为进一步进行电力系统分析作准备。要求：熟悉各种常用应用软件，熟悉硬件设备

2、的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。编制调试电力系统潮流计算的计算机程序。程序要求根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，完成该电力系统的潮流计算，要求计算出节点电压、功率等参数。三、主要仪器设备及耗材每组计算机1台、相关计算软件1套四、实验内容（一)最优潮流的概念最优潮流(OptimalPowerFlow，OPF)是指当系统的结构参数和负荷情况都已给定时，调节可利用的控制变量(如发电机输出功率、可调变压器抽头等)来找到能满足所有运行约束条件的，并使系统的某一性能指标(如发电成

3、本或网络损耗)达到最优值下的潮流分布。经典最优潮流常常在满足可行性约束和安全性约束的条件下追求最小运行费用，或者最小网损、最小负荷、最高电压水平等等。(二)最优潮流的变量：最优潮流的变量分为控制变量（u）及状态变量（x）。一般常用的控制变量有：（1）除平衡节点外，其它发电机的有功出力；（2）所有发电机节点及具有可调无功补偿设备节点的电压模值；（3）移相器抽头位置（4）带负荷调压变压器的变比。（5）并联电抗器/电容器容量状态变量常见的有：（1）除平衡节点外，其它所有节点的电压相角；（2）除发电机节点以及具有可调无功补偿设备节点之外，其它所有节点的电压模值。(三)

4、最优潮流的目标函数电力系统经济调度运行中的最优潮流计算一般以系统运行成本最小为目标，其数学模型如下：（1）全系统发电燃料总耗量（或总费用）f=∑KPi(Gi)iNG∈式中：NG为全系统发电机的集合，其中包括平衡节点s的发电机组。Ki(PGi)是发电机组Gi的耗量特性。由于平衡节点s的电源有功出力不是控制变量，其节点注入功率必须通过潮流计算才能决定，是节点电压模值U及相角θ的函数，于是P=PU(,)θ+PGssLs式中：PS(U,θ)为注入节点s而通过与节点相关的线路输出的有功功率；PLS为节点s的负荷功率。所以（1）式可写成：f=∑Ki(PGi)+Ks(PGs

5、)i∈NGi≠s（2）有功网损，即：f=∑(Pij+Pji)ijNL,∈式中：NL为所有支路的集合。可直接采用平衡节点的有功注入作有功网损最小化的目标函数minf=minP(U,θ)s由上可见，最优潮流的目标函数不仅与控制变量u有关，同时和状态变量x有关。因此可用简洁的形式表示f＝f（u,x）(四)最优潮流的约束条件及其数学模型（1）等式约束条件：最优潮流分布必须满足基本潮流方程，即n⎫PGi−PDi−Vi∑VGj(ijcosθij+Bijsinθij)=0⎪j=1⎪⎬i∈SBnQ−Q−V∑VG(sinθ−Bcosθ)=0⎪GiDiijijijijij⎪j=1

6、⎭该式可简化为：g(x,u)＝0（2）不等式约束条件：1）有功电源出力上下限约束；2）可调无功电源出力上下限约束；3）带负荷调压变压器变比K调整范围约束；4）节点电压模值上下限约束；5）输电线路或变压器元件中通过的最大电流或视在功率约束；6）线路通过的最大有功潮流或无功潮流约束；7）线路两端节点电压相角差约束，等等。部分用不等式表示如下P≤P≤P（i∈S）GiGiGiGQ≤Q≤Q（i∈S）RiRiRiRV≤V≤V（i∈S）iiiB2Pl=Pij=VVGij(ijcosθij+Bijsinθij)−VGiij≤Pl(l∈Sl)SSSSB、G、R、l分别是系统所有

7、节点集合、所有发电机集合、所有无功源PQPQ集合、所有支路集合。Gi、Gi为发电机i的有功、无功出力；Di、Di为节点i的有功、无功负荷；Vi、θi为节点i的电压幅值和相角，其中θi−θj=θij。G、ijB为节点导纳矩阵第i行第j列约束的实部和虚部；P为线路l的有功潮流、设ijl线路l两端节点为i、j。该模型采用的是节点电压极坐标表示形式，当然也可以采用节点电压直角坐标的表示形式。统一表示为h(u,x)<=0则电力系统最优潮流的数学模型可表示为minf(,)ux⎫u⎪⎪st..gux(,)=0⎬(hux,)≤0⎪⎪⎭五、实验数据及处理结果简化梯度法的迭代计算

8、步骤：1）令迭代计数k＝0；(0)2）