电力系统潮流计算课程设计报告

课程设计报告学生姓名:学号：学院:电气工程学院班级:题目:电力系统潮流计算指导教师：职称:副教授指导教师：李翠萍职称:副教授2014年01月10日初壮1潮流计算的目的与意义潮流计算的目的：已知电网的接线方式与参数及运行条件，计算电力系统稳态运行各母线电压、个支路电流与功率及网损。对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。潮流计算的意义：(1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。(2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。(3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。2潮流计算数学模型1.变压器的数学模型：变压器忽略对地支路等值电路：2.输电线的数学模型：π型等值电路：3数值方法与计算流程利用牛顿拉夫逊法进行求解，用MATLAB软件编程，可以求解系统潮流分布根据题目的不同要求对参数进行调整，通过调节变压器变比和发电厂的电压，求解出合理的潮流分布，最后用matpower进行潮流分析，将两者进行比较。牛顿—拉夫逊法1、牛顿—拉夫逊法概要首先对一般的牛顿—拉夫逊法作一简单的说明。已知一个变量X函数为：0)(Xf到此方程时，由适当的近似值)0(X出发，根据：,......)2,1()()()()()()1(nXfXfXXnnnn反复进行计算，当)(nX满足适当的收敛条件就是上面方程的根。这样的方法就是所谓的牛顿—拉夫逊法。这一方法还可以做下面的解释，设第n次迭代得到的解语真值之差，即)(nX的误差为时，则：0)()(nXf把)()(nXf在)(nX附近对用泰勒级数展开0......)(!2)()()()(2)()()(nnnnXfXfXfXf上式省略去2以后部分0)()()()(nnXfXf)(nX的误差可以近似由上式计算出来。)()()()(nnXfXf比较两式，可以看出牛顿—拉夫逊法的休整量和)(nX的误差的一次项相等。用同样的方法考虑，给出n个变量的n个方程：0),,,(0),,,(0),,,(21212211nnnnXXXfXXXfXXXf对其近似解1X得修正量1X可以通过解下边的方程来确定：nnnnnnnnnnnXXXxfxfxfxfxfxfxfxfxfXXXfXXXfXXXf2121222121211121212211),,,(),,,(),,,(式中等号右边的矩阵nnxf都是对于nXXX,,,21的值。这一矩阵称为雅可比（JACOBI）矩阵。按上述得到的修正向量nXXX,,,21后，得到如下关系nnnXXX这比nXXX,,,21更接近真实值。这一步在收敛到希望的值以前重复进行，一般要反复计算满足1112121111,,,maxnnnnnnnnXXXXXX为预先规定的小正数，1nnX是第n次迭代nX的近似值。2、牛顿法的框图及求解过程（1）用牛顿法计算潮流时，有以下的步骤：①给这各节点电压初始值)0()0(,fe；②将以上电压初始值代入公式，求修正方程的常数项向量)0(2)0()0()(,,VQP；③将电压初始值在带入上述公式，求出修正方程中系数矩阵的各元素。④解修正方程式)0()0(,fe；⑤修正各节点电压)0()0()1(eee，)0()0()1(fff；⑥将)1(e，)1(f在带入方程式，求出)1(2)1()1()(,,VQP；⑦检验是否收敛，即)()(,maxkikiQP如果收敛，迭代到此结束，进一步计算各线路潮流和平衡节点功率，并打印输出结果。如果不收敛，转回②进行下次迭代计算，直到收敛为止。（2）程序框图如下：启动输入原始数据形成节点导纳矩阵分解各节点初始电压的实部和虚部迭代次数K=0求PQ节点的)(kiP,)(kiQ，求PV节点的)(kiP，)(kiU置节点号i=0雅克比矩阵是否形成，in求得雅克比矩阵各元素)(kijH)(kijN)(kijJ)(kijL)(kijR)(kijS增大节点号i=i+1把雅克比矩阵单位化求解修正方程，得)(kie，)(kif求解最大修正量max)(||kie，max)(||kif是否收敛回带各电压新值，K=K+1计算输出电压大小及相角，节电功率及支路损耗停止否是否是4算例分析4.1第一问4.1.1节点设置及分类根据系统图可知此系统为两端供电网路，将母线1，2设为节点1，10，将变电所1、2、3、4的高低压侧分别设为节点2、3、4、5、6、7、8、9。并且，将节点1设为平衡节点，将节点10设为PV节点，其余节点设为PQ节点。4.1.2参数求取设定基准值MVASB100，KVUB220，所以4842BBBUSZ根据题目原始资料，计算发电厂、变压器及线路的参数。（1）运用下列公式计算变压器参数：SUPRNNkT221000SUUXNNKT1002%XRZTTTZZZBTT（2）计算线路参数LjxrjXRZ)(（3）变电所负荷分别为：变电所1LS=50+j30.987变电所2LS=40+j27.79变电所3LS=50+j30.987变电所4LS=60+j37.18将参数整理，见下表：首端号末端号阻抗有名值阻抗标幺值电纳有名值电纳标幺值（4）计算变压器分接头变比变压器有5个抽头，电压调节范围为NU2*2.5%，NU对应的分接头开始时设变压器高压侧接主接头,降压变压器5个分接头时的非标准变比*k以备调压时选用955.02201110%)51(2201*k932.02201110%)5.21(2202*k909.022011102203*k886.02201110%)5.21(2204*k864.02201110%)51(2205*k对变电所低压母线为35KV时，非标准变比与10KV时相同。所以得到的电力系统数学模型如下图所示：128.5+j20.10.0176+j0.0415j0.000556j0.26911413.6+j32.160.0218+j0.0664j0.0002224j0.10761613.6+j32.160.0218+j0.0664j0.0002224j0.1076231.49+j40.3330.0031+j0.083300451.78+j53.8850.037+j0.1113004610.2+j24.120.0211+j0.0498j0.0001668j0.0807671.49+j40.3330.0031+j0.083300686.8+j16.080.0140+j0.0332j0.0004448j0.2153891.78+j53.8850.0037+j0.1113008108.5+j20.10.0176+j0.0415j0.000556j0.26914.2第二问当变电所1,2,3,4负荷分别为：40MW，50MW，40MW，60MW时，假设初始的branch，bus，generator矩阵如下图：S=40+j24.8S=50+j31S=40+j24.8S=60+j37.2发电厂一发电厂二1:k11:k21:k31:k4110248635798.5+j20.113.6+j32.161.49+j40.3331.78+j53.8851.49+j40.3331.78+j53.88510.2+j24.1213.6+j32.166.8+j16.088.5+j20.1j0.0001112j0.0001112j0.0005564j0.0000834j0.0004448δ=0U=231VP=200MWU=231V得到的结果如下图：节点1电压节点2电压节点3电压节点4电压节点5电压节点6电压节点7电压节点8电压节点9电压节点10电压有功功率损耗无功功率损耗1.00.9871.0110.9710.9600.9740.9970.9740.9731.014.05543.75由上图分析诸多参数都不符合要求，且网络的无功有功损耗过大。因此对母线电压，发电厂发电功率，变压器变比逐个调节，得到如下关系：1，母线电压1的电压升高会导致以下参数的改变：2~9母线的电压会提高，网络总损耗变大。2，发电厂2的发电功率减少会导致以下参数的改变：4~9母线电压会降低，网络总损耗变小。3，变压器变比的减少会导致以下参数的改变：相应的母线电压提高，网络总损耗变小。综上述分析，经过调试参数得到较为满足条件且网络损耗较小的情况：得到branch，bus，generator矩阵如下图结果如下图：节点1电压节点2电压节点3电压节点4电压节点5电压节点6电压节点7电压节点8电压节点9电压节点10电压有功功率损耗无功功率损耗1.010.9971.0220.9730.9890.9670.9890.9730.9981.05.40623.54.3第三问由于要随着负荷变化而使电力系统满足给类要求，所以我们先初步给定各可控参数：（由于发电厂1的母线设定为参考节点，所以暂且将电厂1的发电量设为0）发电厂2发电量发电厂1母线电压发电厂2母线电压变压器1变比变压器2变比变压器3变比变压器4变比601.0110.9550.9320.9550.932结果如下表：节点1电压节点2电压节点3电压节点4电压节点5电压节点6电压节点7电压节点8电压节点9电压节点10电压有功功率损耗无功功率损耗1.010.9971.0220.9730.9890.9670.9890.9730.9981.05.40623.5由此表中所给的可控参数，我们得出的各电压参数不符合题目要求，且有无功损耗较高，因此我们还需要调整各可控参数来获得更小的有无功损耗。4.4第四问经过调试各可控参数我们得到以下结果：当可调参数的值为下表时：得到的结果如下表所示：节点1电压节点2电压节点3电压节点4电压节点5电压节点6电压节点7电压节点8电压节点9电压节点10电压有功功率损耗无功功率损耗四个负1.010.9971.0230.9861.0030.9911.0150.9911.0171.013.82219.09荷均增大2%四个负荷均减小2%1.010.9981.0240.9871.0070.9921.0180.9921.0211.013.49417.451,4下降2,3上升2%1.010.9981.0240.9861.0030.9911.0160.9911.021.013.70818.23满足题目要求，且有无功损耗较小4.5第五问设变电所1，2，3，4的负荷分别为40MW，50MW，40MW，60MW，此时若不断开任何一条回线且可控参数为第四问得到的较优解，则运行结果如下：节点1电压节点2电压节点3电压节点4电压节点5电压节点6电压节点7电压节点8电压节点9电压节点10电压有功功率损耗无功功率损耗未断开回线1.010.9971.0220.9871.0050.9911.0150.9911.0191