课设论文-电力系统潮流计算课程设计

电力系统潮流计算课程设计任务书（1）设计名称：电力系统潮流计算课程设计设计性质：理论计算，计算机仿真与验证计划学时：两周一、 设计目的 1. 培养学生独立分析问题、解决问题的能力；2. 培养学生的工程意识，灵活运用所学知识分析工程问题的能力3. 编制程序或利用电力系统分析计算软件进行电力系统潮流分析。二、 原始资料1、系统图：IEEE30节点。2、原始资料：见IEEE30节点标准数据库三、 课程设计基本内容：1. 采用PSAT仿真工具中的潮流计算软件计算系统潮流；1) 熟悉PSAT仿真工具的功能；2) 掌握IEEE标准数据格式内容；3) 将IEEE标准数据转化为PSAT计算数据；2. 分别采用NR法和PQ分解法计算潮流，观察NR法计算潮流中雅可比矩阵的变化情况，分析两种方法计算潮流的优缺点；3. 分析系统潮流情况，包括电压幅值、相角，线路过载情况以及全网有功损耗情况。4. 选择以下内容之一进行分析：1) 找出系统中有功损耗最大的一条线路，给出减小该线路损耗的措施，比较各种措施的特点，并仿真验证；2) 找出系统中电压最低的节点，给出调压措施，比较各种措施的特点，并仿真验证；3) 找出系统中流过有功功率最大的一条线路，给出减小该线路有功功率的措施，比较各种措施的特点，并仿真验证；5. 任选以下内容之一作为深入研究：（不做要求）1) 找出系统中有功功率损耗最大的一条线路，改变发电机有功出力，分析对该线路有功功率损耗灵敏度最大的发电机有功功率，并进行有效调整，减小该线路的损耗；2) 找出系统中有功功率损耗最大的一条线路，进行无功功率补偿，分析对该线路有功功率损耗灵敏度最大的负荷无功功率，并进行有效调整，减小该线路的损耗；3) 找出系统中电压最低的节点，分析对该节点电压幅值灵敏度最大的发电机端电压，并有效调整发电机端电压，提高该节点电压水平；四、 课程设计成品基本要求：1. 绘制系统潮流图，潮流图应包括：1) 系统网络参数2) 节点电压幅值及相角3) 线路和变压器的首末端有功功率和无功功率2. 撰写设计报告，报告内容应包括以下几点：1) 本次设计的目的和设计的任务；2) 电力系统潮流计算的计算机方法原理，分析NR法和PQ分解法计算潮流的特点；3) 对潮流计算结果进行分析，评价该潮流断面的运行方式安全性和经济性；4) 找出系统中运行的薄弱环节，如电压较低点或负载较大线路，给出调整措施；5) 分析各种调整措施的特点并比较它们之间的差异；6) 结论部分以及设计心得；五、 考核形式1. 平时表现：3人一组，考查本组学生平时纪律、态度等；2. 报告质量：设计成品的完成质量、撰写水平等；3. 答辩考核：参照设计成品，对计算机方法进行电力系统潮流计算的相关问题等进行答辩；4. 考核最终成绩由平时表现、报告成绩、答辩成绩组成。电力系统潮流计算摘要：电力系统潮流计算是电力系统分析的基础内容。通过对电力网络的潮流计算我们可以对系统的安全性，经济型，稳定性进行分析判断。并且能够准确的，快速的对电力系统中出现的各种问题进行调整。潮流计算是电力系统最基本，最常用的计算。Newton-Raphson法是数学上解非线性方程式的有效方法，有较好的收敛性。本次课程设计的主要内容就是基于PSAT的潮流计算和分析。针对问题一：通过题目中给的ieee格式的数据进行分析和理解，将其转换为PSAT标准格式的数据。然后使用PSAT利用N-R法对题目中的6机30节点的电力网络进行潮流计算。针对问题二：通过对PSAT中的N-R法的中断控制可以得到N-R法每次迭代的雅克比矩阵。通过对这些雅克比矩阵的分析得到雅克比矩阵的变化规律，即每次迭代出来的雅克比矩阵的个元素都在减小并且变化率也越来越小。利用P-Q法再次对题目中的数据进行潮流计算，得到的潮流的结果与N-R法求得的结果一样，计算的时间短于N-R法。然而，利用N-R法计算潮流其收敛性要好于P-Q分解法。针对问题三：利用PSAT中的simulink功能画出此6机30节点的系统。将目中的数据以PSAT标准格式填入此网络中。并对此网络进行仿真运行得到voltage mangnitudes图，voltage angles图，line flows图。通过对这三个图的分析得到此网络的电压水平和潮流情况。针对问题四：通过对voltage mangntudes图的观察可以看出30节点的电压较低，需要提高30节点的电压。因此我采用了4种方法来提高30节点的电压。(1)提高8节点发电机的机端电压，提高30节点的电压。(2)提高27和28节点之间的变压器变比，提高30节点的电压。(3)在27节点和30节点之间串联电容，提高30节点的电压。(4)在母线10或24处并联电容器，提高30节点的电压。关键词：PSAT 潮流计算 牛顿-拉夫逊法 P-Q分解法二课程设计准备工作2.1牛顿-拉夫逊法原理及应用牛顿-拉夫逊法按照电压的表示方法不同，又分为直角坐标形式和极坐标形式，牛顿-拉夫逊法潮流计算具有二阶收敛特性，计算中收敛速度较快，但是当导纳矩阵阶数较高时，初值敏感性问题突出，如果初值选取的合适，则计算会非常快速准确。若初值选取不合适则计算有可能会不收敛。牛顿-拉夫逊法设计原理：设有非线性方程组 ； ； ；其近似解。近似解与精确解分别相差则有如下：；对上述每一项式按泰勒级数展开忽略的高次方，可得如下修正方程组：；简写为： ；其中：J称为函数的雅可比矩阵。求解功率方程步骤如下：；代入及得：；将上式的实部与虚部分开列写如下：；及分别为i节点的电压幅值及相角，及分别为i节点注入的有功功率、无功功率。修正方程式如下：有功功率不平衡量、无功功率不平衡量、分别由下式计算：1、 当时雅可比矩阵各个元素分别为：；2、 当时雅可比矩阵各个元素为：；牛顿-拉夫逊法潮流计算步骤有以下几步：(1)形成节点导纳阵。(2)设各节点的电压的初值。(3)利用各节点的电压初值计算修正方程的不平衡量。(4)利用各节点电压初值求出雅克比矩阵的各个元素(5)解修正方程，求出各节点电压的变化量，即修正量。如果符合要求则跳出，如果不符合要求则继续下一步。(6)计算各节点电压的新值，即修正后的值。(7)御用各节点电压的新值自第三步开始下一次迭代。(8)计算平衡节点的功率和线路功率。牛顿-拉夫逊法计算框图：输入原始数据形成节点导纳阵启动设节点电压初值，置迭代次数对PQ节点：计算、对PV节点：计算、置节点号i=1雅克比矩是否已全部形成计算雅克比矩阵元素增大节点号i=i+1否解修正方程式，由和J求各节点电压变量、；i=1,2,n,is是求出、迭代是否收敛，、？计算平衡节点功率和线路功率停止是计算各节点电压新值i=1,2,n;is增大迭代次数，k=k+1否图一：牛顿-拉夫逊法潮流计算流程图2.2 P-Q分解法原理P-Q分解法是极坐标牛顿-拉夫逊法的一种简化算法快速分解法，有两个主要特点: (1)降阶在潮流计算的修正方程中利用了有功功率主要与节点电压相位有关,无功功率主要与节点电压幅值有关的特点,实现P-Q分解,使系数矩阵由原来的2N2N阶降为NN阶,N为系统的节点数(不包括缓冲节点)。(2)因子表固定化利用了线路两端电压相位差不大的假定,使修正方程系数矩阵元素变为常数,并且就是节点导纳的虚部。由于以上两个特点,使快速分解法每一次迭代的计算量比牛顿法大大减少。P-Q分解法只具有一次收敛性,因此要求的迭代次数比牛顿法多,但总体上快速分解法的计算速度仍比牛顿法快。快速分解法只适用于高压网的潮流计算,对中、低压网,因线路电阻与电抗的比值大,线路两端电压相位差不大的假定已不成立,用快速分解法计算,会出现不收敛问题。P-Q分解法的计算步骤与牛顿-拉夫逊法的计算步骤略有不同，在开始迭代之前就形成了系数矩阵和并得到了他们的逆阵。P-Q分解法计算步骤如下：(1)形成节点导纳矩阵。(2)形成系数矩阵,。(3)给定电压初值，。(4)求功率不平衡量。(5)接修正方程式得。(6)求功率不平衡量。(7)解修正方程式得。(8)修正电压初值。(9)判断若满足条件则计算功率,若不满足条件则返回第3步重新求解。P-Q分解法的计算框图如下：启动输入原始数据形成节点导纳阵形成阵，并求其逆阵形成阵，并求其逆阵设节点电压：和（i=1,2,3,n）设置迭代次k=0置=0，=0计算，i=1,2,n is迭代是否收敛，？是置=0=0？是计算平衡节点功率和线路功率计算，(i=1,2,n,is)计算新的值，增大迭代次数kk+1置置增大迭代次数kk+1计算,(i=1,2,m,is)迭代是否收敛，？计算，(i=1,2,n,is)计算新的值，置?结束是否是否否图二：P-Q分解法流程图三、课程设计问题分析与解答3.1采用PSAT仿真工具中的潮流计算软件计算系统潮流3.1.1对原始数据进行分析将我的数据30IEEE.DAT格式的数据导入到psat软件中生成一个m-file。这个m-file是psat标准数据格式。具体数据见表一至表六。表一：母线参数母线号电压基准值电压标幺值电压相角11001.060.0000021001.03-0.0869231001.02-0.2218341001.018-0.1801251001.02-0.1532761001.05-0.2481971001.062-0.2331881001.05-0.2331891001.056-0.26075101001.051-0.26354111001.057-0.25813121001.055-0.26302131001.05-0.26442141001.035-0.27978表二：平衡节点参数节点号功率基准值(MVA)电压基准值(KV)电压幅值电压相角最大无功功率最小无功功率最大电压最小电压有功功率11001001.0609999-999.91.10.90.50386表三：pv节点参数节点号功率基准值电压基准值有功功率电压幅值最大无功最小无功最大电压最小电压21001000.41.030.5-0.41.10.931001000.51.020.401.10.961001000.51.050.24-0.061.10.9810010011.050.24-0.061.10.9表四：pq节点的参数节点号功率基准值电压基准值有功功率无功功率最大电压最小电压21001000.2170.1271.10.931001000.5420.091.10.941001000.1780.0391.10.951001000.5760.1161.10.961001000.1120.0751.10.991001000.2950.1661.10.9101001000.090.0581.10.9111001000.5350.3181.10.9121001000.1610.1161.10.9131001000.1350.0881.10.9141001000.1490.071.10.9表五：对地电纳参数节点号功率基准值电压基准值频率电导电感91001006000.19表六：线路参数线路始端线路末端功率基准值电压基准值频率线路长度单位电阻单位电抗单位电纳线路变比5610010060000.2520200.9324710010060000.2091200.9784910010060000.5561800.969121001006000.019380.059170.05280231001006000.046990.197970.04380241001006000.058110.176320.03740151001006000.054030.223040.04920251001006000.056950.173880.0340341001006000.067010.171030.03460451001006000.013350.042110.012807810010060000.17615007910010060000.11001009101001006000.031810.0845006111001006000.094980.1989006121001006000.122910.25581006131001006000.066150.13027009141001006000.127110.270380010111001006000.082050.192070012131001006000.220920.199880013141001006000.170930.34802003.1.2对数据进行求解通过psat生成这些数据之后，我利用牛顿-拉夫逊法进行求解得到的结果如下：GLOBAL SUMMARY REPORT总发电机功率视在功率p.u. 3.0452 有功功率p.u. 1.1919 总负荷功率视在功率p.u. 2.99 有功功率p.u. 1.0675 总损耗视在功率p.u. 0.05524 有功功率p.u. 0.1244 表7：POWER FLOW RESULTSBusVPhaseP genQ genP loadQ loadBUS-11.0600.645240.5501300BUS-100.99364-0.09121000.090.058BUS-110.96366-0.13733000.5350.318BUS-121.0074-0.12373000.1610.116BUS-131.0175-0.11778000.1350.088BUS-140.99376-0.10597000.1490.07BUS-21.03-0.011920.4-0.226180.2170.127BUS-31.02-0.030750.50.056480.5420.09BUS-41.0104-0.03924000.1780.039BUS-51.0075-0.0506000.5760.116BUS-61.05-0.106790.50.569380.1120.075BUS-71.02320.013210000BUS-81.050.1779110.2421100BUS-91.0143-0.066000.295-0.02948表八：LINE FLOWSFROM BUSTO BUSLINEP FLOWQ FLOWP LOSSQ LOSSBUS-5BUS-610.219740.121800.01567BUS-4BUS-72-0.253460.053200.01374BUS-4BUS-930.047780.0505500.00264BUS-1BUS-240.358030.391820.00528-0.04157BUS-2BUS-350.106480.004460.00054-0.04376BUS-2BUS-460.180170.037750.00196-0.03298BUS-1BUS-570.287210.158310.00563-0.02937BUS-2BUS-580.24910.037990.0035-0.02461BUS-3BUS-490.063950.01470.00033-0.03481BUS-4BUS-5100.26951-0.02250.00095-0.01002BUS-7BUS-811-1-0.0729700.16914BUS-7BUS-9120.746540.1124300.05989BUS-9BUS-10130.346470.121680.004170.01108BUS-6BUS-11140.304640.31270.016420.03438BUS-6BUS-12150.125370.115180.003230.00672BUS-6BUS-13160.177730.172630.003680.00725BUS-9BUS-14170.152850.008260.002890.00616BUS-10BUS-11180.25230.05260.005520.01292BUS-12BUS-1319-0.03886-0.007550.000340.00031BUS-13BUS-1420-0.000150.069520.00080.00163表九：LINE FLOWFROM BUSTO BUSLINEP FLOWQ FLOWP LOSSQ LOSSBUS-6BUS-51-0.21974-0.1061300.01567BUS-7BUS-420.25346-0.0394600.01374BUS-9BUS-43-0.04778-0.0479100.00264BUS-2BUS-14-0.35276-0.433380.00528-0.04157BUS-3BUS-25-0.10595-0.048220.00054-0.04376BUS-4BUS-26-0.17821-0.070740.00196-0.03298BUS-5BUS-17-0.28158-0.187690.00563-0.02937BUS-5BUS-28-0.2456-0.06260.0035-0.02461BUS-4BUS-39-0.06362-0.049510.00033-0.03481BUS-5BUS-410-0.268560.012480.00095-0.01002BUS-8BUS-71110.2421100.16914BUS-9BUS-712-0.74654-0.0525400.05989BUS-10BUS-913-0.3423-0.11060.004170.01108BUS-11BUS-614-0.28822-0.278320.016420.03438BUS-12BUS-615-0.12214-0.108450.003230.00672BUS-13BUS-616-0.17405-0.165380.003680.00725BUS-14BUS-917-0.14995-0.00210.002890.00616BUS-11BUS-1018-0.24678-0.039680.005520.01292BUS-13BUS-12190.03920.007850.000340.00031BUS-14BUS-13200.00095-0.06790.00080.001633.2观察NR法计算潮流中雅可比矩阵的变化情况，分别采用NR法和PQ分解法计算潮流，分析两种方法计算潮流的优缺点3.2.1观察牛顿-拉夫逊法雅克比矩阵找到psat里fm_spf.m这段子程序，通过对“if Settings.distrsw, DAE.kg = DAE.kg + inc(end); end”这段命令的中断控制得到每次迭代后的雅克比矩阵。雅克比矩阵见下：第一次迭代雅克比矩阵（部分，下同）：10000000032.33471-5.13856-5.50541-5.563230000-4.817799.990482-5.1726900000-5.17635-5.3446739.64437-22.20750-5.0489834100-5.326240-22.588936.21387-3.94286000000-3.9428618.8337100000-5.048980021.56978277-6.3304第二次迭代雅克比矩阵：10000000032.56461-5.04192-5.41028-5.488520000-5.004510.29105-5.2865500000-5.31507-5.2401939.67242-22.35420-4.938339200-5.355050-22.20435.9915-3.93372000000-3.9337218.1701100000-4.938340020.73693015-6.11274第三次迭代的雅克比矩阵：10000000032.46863-5.04534-5.37035-5.454660000-5.0005410.24205-5.2415100000-5.27457-5.2066839.14478-22.04790-4.830063100-5.318290-21.8935.57703-3.90649000000-3.9064918.0510500000-4.830060020.2584549-6.01845第四次迭代雅克比矩阵：10000000032.46738-5.04538-5.3698-5.454210000-5.0004910.24139-5.240900000-5.27396-5.2061639.13689-22.04360-4.8282732300-5.31770-21.885435.57086-3.90606000000-3.9060618.0490500000-4.828270020.24990738-6.01676已知雅克比矩阵的基本形式为：其中各个元素为， ，。这些元素都是关于电压的函数，所以雅克比矩阵的变化时随着电压的变化而变化的。因为每次迭代电压都会更接近真实值，各个元素的数值也是在不断减小知道满足精度要求。通过对比分析雅克比矩阵的各个元素，找到了以下规律：(1)雅克比矩阵为稀疏阵，但不是对称阵。这是因为。(2)雅克比矩阵的各个元素的数值在不断的减小，直至迭代结束满足精度要求。(3)各个元素的变化率在不断减小，这是因为每一次迭代都会更加接近真实值，在运行刚刚开始的时候数值变化较大，随着迭代的次数不断增加变化越来越小。3.2.2 运用P-Q分解法进行潮流计算并与牛顿-拉夫逊法进行对比P-Q分解法是源自于牛顿-拉夫逊法的极坐标表示形式，它是基于对修正方程的两个简化假设得到的。对修正方程的第一个简化是：计及电力网络中各元件的电抗远大于电阻，以至于节点电压相位的改变主要影响各元件的有功功率潮流，从而影响为节点的注入有功功率；各节点电压大小的影响主要影响注入的无功功率。因此将雅克比矩阵中的N、J略去。对修正方程的第二个假简化是：基于对状态变量的的约束条件，从而可以得出,。与牛顿-拉夫逊法相比P-Q分解法有以下特点(1)以一个(n-1)阶和一个(m-1)阶的系数矩阵和代替原有的(n+m-2)阶的雅克比矩阵，提高了计算的速度。(2)在迭代过程中保持不变的系数矩阵和代替变化的雅克比矩阵，显著的提高了计算的速度。(3)以对称的系数矩阵和代替不对称的雅克比矩阵，是求逆等运算量和所需的储存容量都大为减少。首先在psat中导入原始数据，将潮流计算方法改为XB fast decoupled。然后进行潮流计算。通过与牛顿-拉夫逊法得到的结果相对比，发现两种计算结果是一样的。运用牛顿-拉夫逊法迭代次数与迭代时间为：PF solver: Newton-Raphson methodSingle slack bus modelIteration = 1 Maximum Convergency Error = 0.40573Iteration = 2 Maximum Convergency Error = 0.015461Iteration = 3 Maximum Convergency Error = 0.00040656Iteration = 4 Maximum Convergency Error = 2.0894e-007Power Flow completed in 0.047 s用P-Q分解法计算迭代次数与迭代时间为：PF solver: XB fast decoupled methodSingle slack bus modelIteration = 1 Maximum Convergency Error = 0.6258Iteration = 2 Maximum Convergency Error = 0.074761Iteration = 3 Maximum Convergency Error = 0.010438Iteration = 4 Maximum Convergency Error = 0.0005484Iteration = 5 Maximum Convergency Error = 0.00013099Iteration = 6 Maximum Convergency Error = 2.1114e-005Iteration = 7 Maximum Convergency Error = 1.1852e-006Power Flow completed in 0.046 s将以上数据导入到EXCEL中，利用其绘图功能，绘制出两种方法的对数收敛特性曲线进行直观的分析。图三：牛顿-拉夫逊法和P-Q分解法收敛特性图通过对以上信息的分析可以得到以下结论：(1)P-Q分解法与牛顿-拉夫逊法所算出潮流结果是一样的。(2)P-Q分解法迭代次数要多于N-R法，但是总需要时间与N-R法相差无几。故P-Q分解法每次迭代所要的时间要小于N-R法。(3)通过对P-Q和N-R收敛特性图的仔细比较，会发现P-Q分解法收敛速度要好于N-R法。这是因为N-R法事具有二阶收敛特性，而P-Q分解法事具有一阶收敛特性的。(4) P-Q分解法是利用了电网的高电压特性所以只适用于高压网的潮流计算,对中、低压网,因线路电阻与电抗的比值大,线路两端电压相位差不大的假定已不成立,用快速分解法计算,会出现不收敛问题。3.3 分析系统潮流情况，包括电压幅值、相角，线路过载情况以及全网有功损耗情况3.3.1利用psat绘制出5机14节点网络图在psat中利用psat simulink library功能绘制出网络图，并将原始数据按照psat的标准数据格式填入到此网络中。simulink library仿真图可以很直观的表示出整个网络的运行状况，有助于分析和解决问题。网络见下图：图四：simulink library仿真图及潮流运算结果3.3.2分析电压幅值、相角、和潮流分布情况通过对图五的仿真运行，得到各节点电压幅值并且可以绘制出电压幅值network bisualisation图。 表十：电压幅值表节点号电压幅值Bus011.06Bus021.03Bus031.02Bus041.0104Bus051.0075Bus061.05Bus071.0232Bus081.05Bus091.0143Bus100.99364Bus110.96366Bus121.0074Bus131.0175Bus140.99376图五：全网电压幅值图结论：通过对全网电压图和电压幅值表的分析，很容易看出11节点的电压最低。1节点的电压高同样的方法绘制出电压相角图： 表十一：电压相角表节点号电压相角Bus010Bus02-0.01192Bus03-0.03075Bus04-0.03924Bus05-0.0506Bus06-0.10679Bus070.01321Bus080.17791Bus09-0.066Bus10-0.09121Bus11-0.13733Bus12-0.12373Bus13-0.11778Bus14-0.10597图六：电压相角图结论：通过对电压相角图和电压相角表的分析，可以看出6节点电压相角最小，11节点电压相角最大。做出全网功率潮流图：图七：全网潮流图表十二：各条线路的有功和无功线路始端线路末端有功功率无功功率线路始端线路末端有功功率无功功率Bus02Bus01-0.35276-0.43338Bus01Bus020.358030.39182Bus02Bus030.106480.00446Bus03Bus02-0.10595-0.04822Bus06Bus110.304640.3127Bus11Bus06-0.28822-0.27832Bus06Bus120.125370.11518Bus12Bus06-0.12214-0.10845Bus06Bus130.177730.17263Bus13Bus06-0.17405-0.16538Bus09Bus140.152850.00826Bus14Bus09-0.14995-0.0021Bus11Bus10-0.24678-0.03968Bus10Bus110.25230.0526Bus12Bus13-0.03886-0.00755Bus13Bus120.03920.00785Bus13Bus14-0.000150.06952Bus14Bus130.00095-0.0679Bus02Bus040.180170.03775Bus04Bus02-0.17821-0.07074Bus01Bus050.287210.15831Bus05Bus01-0.28158-0.18769Bus02Bus050.24910.03799Bus05Bus02-0.2456-0.0626Bus03Bus040.063950.0147Bus04Bus03-0.06362-0.04951Bus05Bus04-0.268560.01248Bus04Bus050.26951-0.0225Bus08Bus0710.24211Bus07Bus08-1-0.07297Bus07Bus090.746540.11243Bus09Bus07-0.74654-0.05254Bus10Bus09-0.3423-0.1106Bus09Bus100.346470.12168Bus05Bus060.219740.1218Bus06Bus05-0.21974-0.10613Bus04Bus07-0.253460.0532Bus07Bus040.25346-0.03946Bus04Bus090.047780.05055Bus09Bus04-0.04778-0.04791结论：可以看出BUS07到BUS08流过的有功功率最大，BUS12到BUS13流过的功率最小。BUS07到BUS08之间的线路将会过载。3.4找出系统中电压最低的节点，给出调压措施，比较各种措施的特点，并仿真验证通过上面的分析得出11节点的电压最低，我们知道调压的措施总共有四种：(1)借改变发电机机端电压调压(2)借改变变压器变比调压(3)串联电容调压(4)并联补偿机调压通过观察整个网络，发现全网共有发电机5台其中一台为平衡节点;有3台变压器;11节点连接着6节点和10节点。因此我的调压方案就是分别将4台pv发电机的机端电压提升到原来的1.05倍，但不超过其发电机的运行极限；分别将调整3台变压器的变比；分别将11节点所连接的两个之路的电抗减少0.03；在11节点并联一个0.3的电容。经过调压之后要对以上四类措施进行一下横向比较和纵向比较，选出一种更加经济、安全、稳定的调压措施。3.4.1借改变发电机极端电压调压在simulink 仿真网络中对每台发电机分别进行调压，调压幅度为增大5%并且不使其超过发电机运行极限。然后对其仿真运行，记录数据，观察调压的效果。数据见下表：表十三：改变发电机端电压调压节点号原始电压幅值6节点发电机调至1.12节点发电机调制1.0858节点发电机调制1.13节点发电机调至1.071Bus011.061.061.061.061.06Bus021.031.031.08151.031.03Bus031.021.021.021.021.071Bus041.01041.01841.03061.01581.0228Bus051.00751.01731.02831.01111.0154Bus061.051.11.051.051.05Bus071.02321.03581.03241.04791.0288Bus081.051.051.051.11.05Bus091.01431.03661.02311.03171.0197Bus100.993641.02181.0011.00830.99816Bus110.963661.00490.967670.971620.96612Bus121.00741.05721.00811.00881.0079Bus131.01751.06461.01881.02011.0183Bus140.993761.02760.999421.0050.99724并且运用excel画出电压幅值对比图以便于观察分析电压的变化情况：图八：4台发电机调压效果散点图通过对数据表和数据图的分析，可以看出只有调节6节点的发电机会明显使11节点的电压升高，并且调节完之后会使全网的电压有明显的改善。使原来其它不符合要求的节点的电压也能有一定的提高。而调节其它发电机并没有使11节点的电压有明显改善，并且对全网电压的也没有更好的提高，仅使它周围的电压有所改善。因此调发电机机端电压方法应选择改变6节点发电机的机端电压。3.4.2借改变变压器变比调压在simulink仿真图中分别将每台变压器的变比分别提高1.05倍，仿真运行之后，记录数据并分析结果。调压前后数据见下表：表十四：改变变比调压节点号原始电压幅值将5、6节点变压器变比提高至1.05将4、7节点变压器变比提高至1.05将4、9节点变压器变比提高至1.05将4、7节点变压器变比降至0.85将4、9节点变压器变比降至0.85Bus011.061.061.061.061.061.06Bus021.031.031.031.031.031.03Bus031.021.021.021.021.021.02Bus041.01041.02331.01821.01040.994121.0045Bus051.00751.02871.01231.00750.997551.004Bus061.051.051.051.051.051.05Bus071.02321.0290.999881.02321.06051.0303Bus081.051.051.051.051.051.05Bus091.01431.01970.999411.01431.03771.0322Bus100.993640.99820.981060.993641.01331.0087Bus110.963660.966130.956790.963660.974380.97187Bus121.00741.00791.00621.00741.00921.0088Bus131.01751.01831