四川大学电力系统计辅期末考试13

四川大学期末考试试题（开卷）（2015——2016学年第1学期）A卷课程号：303033040课序号：课程名称：电力系统分析的计算机算法任课教师：成绩：适用专业年级：电气工程及其自动化专业13级学生人数：印题份数：学号：姓名：考生承诺我已认真阅读并知晓《四川大学考场规则》和《四川大学本科学生考试违纪作弊处分规定（修订）》，郑重承诺：1、已按要求将考试禁止携带的文具用品或与考试有关的物品放置在指定地点；2、不带手机进入考场；3、考试期间遵守以上两项规定，若有违规行为，同意按照有关条款接受处理。考生签名：1、如下图1所示，一个有三台发电机和五个节点的系统的等值电路，图中标明了各元件的参数（均已归算至统一的基准值）。试列写节点导纳矩阵。（20分）2G10.2j1.05:13G1G23450.1j0.2j0.2j0.4j1.50.9j3.51.2j3.51.8j0.2j0.2j0.1j图1系统图和等值电路解：逐个计算导纳矩阵元素如下：111112131221121331132222212324233223244224333331(0.23.51.8)(0.2)(0.4)3.52.20.2,0.4(1.50.9)(0.2)(0.2)0.11.51.20.2,0.1YyyyjjjjjYYyjYYyjYyyyyjjjjjYYyjYYyjYyyy32352353553354444425555530.1(3.51.2)(0.4)(0.2)3.51.711.050.10.0951.050.20.10.30.20.10.3jyjjjjyjYYjkYyyjjjYyyjjj3.52.20.20.4000.21.51.20.20.100.40.23.51.7100.09500.100.30000.09500.3NjjjjjjjYjjjjjjjj评分标准：对角元素计算每个2分共5*2=10分；非对角元素计算每个1分，共5*1=5；列写出节点导纳矩阵的每一列，各1分，共5分。2、牛顿-拉夫逊法的基本原理是什么？其潮流计算的修正方程式是什么？用直角坐标形式表示的与用极坐标形式表示的不平衡方程式的个数有何不同？（20分）a.牛顿法其实就是切线法，通过逐步线性化的方法，反复迭代使解接近准确值。在每一步迭代过程当中，建立修正方程式并求解修正方程式，反复迭代直至满足条件为止(2分)对于多变量非线性方程组，()0FX，在第k次迭代时，建立修正方程式()()()()kkkFXJX和(1)()()kkkXXX。（2分，写成单变量方程扣1分）式中：X和ΔX分别由n个变量和修正量组成的n维列向量；F(X)是由n个多元函数组成的维列向量；J是n*n阶方阵，称为雅可比矩阵，它第i行、第j列的元素iijjfJx第i个函数12(,,,)infxxx对第j个变量jx的偏导数；上角标（k）表示J阵的每一个元素都在点()()()12(,,,)kkknxxx处取值。（1分）迭代过程一直进行到满足收敛判据()()()121max(,,,)kkkinfxxx或()2kix。（1分）b.潮流计算的修正方程式有两种，分别为直角坐标系下和极坐标系下。直角坐标系下，节点电压可表示为iiiijijijVeifYGjB。修正方程式表示为WJV。（3分）其中22111111[,,,,,,,,,]TmmmmnnWPQPQPVPV，111111[,,,,,,,,,]TmmmmnnVefefefPf，雅克比矩阵J中的个元素通过iP，iQ和2iV分别对je和jf求偏导数获得。（2分）极坐标下修正方程式为12DPHNVVQKL（3分）其中121nPPPP，12mQQQQ，121n，12mVVVV，122DmVVVV；H是(n-1)*(n-1)阶方阵，其元素为iijjPH；N是(n-1)*m阶矩阵，其元素为iijjjPNVV；K是m*(n-1)阶矩阵，其元素为iijjQK；L是m*m阶方阵，其元素为iijjjQLVV。（2分）c.直角坐标系下修正方程式个数为2n-2个，极坐标系下为n-1+m个（m为系统PQ节点个数）。极坐标系下修正方程式个数比直角坐标系下少n-1-m个。（4分）3、试说明PQ分解法是如何从牛顿法简化而来的？其修正方程与牛顿法的相比有何特点？（20分）a、PQ分解法根据电力系统实际运行状态的物理特点，对极坐标形式的牛顿法修正方程式进行了合理简化。极坐标形式的牛顿法修正方程式表示为12DPHNVVQKL。由于输电线路中电阻远小于电抗，系统中母线有功功率的变化主要受电压相位影响，无功功率的变化主要受电压幅值变化的影响。近似认为PV和Q为0，即将极坐标形式的牛顿法修正方程中的子块N、K忽略不计，认为他们的元素都等于零。（3分）这样，n-1+m阶的修正方程式便分解为一个n-1阶和一个m阶方程12DPHQLVV。这样，有功不平衡量只修正电压相位，无功不平衡量只修正电压幅值。这就是所谓的有功-无功功率分解法，也成为快速解耦法。（3分）但是矩阵H和L的元素都是节点电压幅值和相位差的函数，其数值在迭代过程中是不断变化的算法最关键的一步简化就在于，把系数矩阵H、L化简为常数对称矩阵。（2分）其根据如下：1）一般情况下，线路两端电压的相位差不大，因此可以认为cos1,sinijijijijGB（1分）2）与系统个节点无功功率相对应的导纳2/iiQV通常远小于该节点自导纳的虚部iiB，即2iiiiQBV或2iiiiQVB（1分）计及以上两点的关系，矩阵H和L个元素的表达式可简化写成(,1,2,,1)(,1,2,,)ijijijijijijHVVBijnLVVBijm（2分）从而得到简化了的修正方程式11112'''DDDVPBVVQBV（2分）b、1）用一个n-1阶和一个m阶的线性方程组代替了一个n-1+m阶线性方程组（2分）2）系数矩阵B’和B’’为常数矩阵。（2分）3）系数矩阵是对称矩阵（2分）4、试问电力系统暂态稳定计算的数学模型主要包括哪些微分方程和代数方程？（15分）发电机的数学模型包括四组方程式：1）机械暂态过程方程式：即发电机转子运动方程，描述转子角度与转速岁原动机与发电机间不平衡功率e()TPP的动态变化过程。（2分）2）电磁暂态过程方程式：描述发电机内电动势（磁通）和电流的暂态变化过程。（2分）3）原动机调速器方程式：描述原动机输出的机械功率TP岁发电机转速的动态变化过程（2分）4）励磁调节系统方程式：描述发电机励磁调节系统的输出电压fV随发电机端电压tV的动态变化过程。（2分）负荷可以根据具体情况选用恒定阻抗模型、考虑感应电动机机械暂态过程或考虑感应电动机机电暂态过程的综合负荷动态模型。（2分）电力网络一般采用节点导纳矩阵形式的代数方程描述。发电机方程和负荷方程通过他们的端电压GV、LV和注入电流GI、LI（或注入功率）与网络方程联系起来。（3分）于是，暂态稳定计算全系统的数学模型包含描述电力系统有关元件动态特性的微分方程和电力网络的代数方程。其一般形式可以写成(,)(,)0xfxygxy。式中：x，y分别为状态变量和代数变量；代数变量y表示电力网络的运行参数，由于在暂态稳定计算中一般不考虑电力网络的暂态过程，因此，y在计算过程中是可以突变的量。（2分）5、两机电力系统等值电路如下图所示。设发电机均装有比例式励磁调节器，采用dX后电势E恒定的模型，试计算发电机G1的功率极限及稳定储备系数。系统参数如图所示。给定运行条件为：0.10V，20201.50.7PjQj。（25分）2.51.0LDSj10.8Xj20.3Xj0.10VG1G21010PjQ2020PjQ解：（1）给定方式的潮流计算101020202.51.01.50.710.3LDPjQSPjQjjj（1分）22'''10110110110110000221100.30.810.80.301.476PRQXPXQREVVV（2分）22''22'2022022000010.70.31.50.31.291QXPXEVVV（2分）（1分）（2）计算负荷阻抗和输入阻抗、转移阻抗20222212.510.3450.1380.37221.82.51LDLDLDLDLDLDLDVZRjXPjQjjPQ（2分）''11112//0.80.3450.138//0.30.10.9730.97884.13LDZRjXZjXjjjj（2分）（1分）''''121212//0.61.341.468114.12LDjXjXZjXjXjZ（2分）1224.12（1分）（3）各发电机的功率特性（4分）（4）计算发电机G1的功率极限及稳定储备系数112121.297cos(24.12)0GdPd（2分）129024.1265.88m（1分）10.2281.2971.525GmP（1分）22101010101201GPQPPRV115.8721010201111212111212sinsin0.2281.29()7sin(24.|Z||Z12)|GEEEP110()101.52481.0100%100%1.052.5%GmGsmPGPPKP（3分）第1页，共2页试卷编号：第1页，共2页