电力系统分析习题集(第四章)

电力系统分析习题集（第四章）【例4-1】试求解图2-8的简单系统的最优潮流。【解】除由图2-8提供的系统母线负荷功率数据、线路参数和变压器支路参数数据、变压器变比数据（非标准变比在首端）之外，以下顺序给出线路传输功率边界、发电机有功、无功出力上下界和燃料耗费曲线参数(燃料耗费曲线所用有功功率变量为标幺值)。若不作说明所有数据都是以标幺值形式给出，功率基准值为100MVA，母线电压上下界分别为1.1和0.9。表4-1线路传输功率边界支路号首末端母线号线路传输功率边界11-2221-30.6532-3242-4653-55表4-2发电机数据发电机序号母线号出力上界出力下界燃料耗费曲线参数有功无功有功无功二次系数一次系数常数14831-350.4395200.43351200.648525851-2.1200.55500.7461857.201首先，我们先列出该算例的数学模型和有关计算公式。在该算例中，共有节点5个，相应的状态变量为：5544332211~VVVVVx系统中有2台发电机，没有其它无功源，因此控制变量为：2121~RRGGQQPPu应该指出，此处发电机和无功源的编号与节点编号无关，是独立编号的。这是因为系统中一个节点可能接有多台发电机的缘故。因此系统中总变量共14个：55443322112121VVVVVQQPPRRGGx。最优潮流的数学模型为：目标函数：.min)()(022122222011112121aPaPaaPaPaGGGG约束条件：每个节点有2个潮流方程，共有10个等式约束条件。对非发电机节点：0)sincos(51ijjijijijjiDiiBGVVPP0)cossin(51ijjijijijjiDiiBGVVQQ（3,2,1i）对发电机节点：0)sincos(51ijjijijijjiDiikGkiBGVVPPP0)cossin(51ijjijijijjiDiikGkiBGVVQQQ（5,4i）（ik表示第k台发电机接在节点i上，41k，52k）不等式约束条件共有14个，分别为：GiGiGiPPP（2,1i）RiRiRiQQQ（2,1i）iiiVVV（5,...,1i）ijijijPPP（对所有5条支路）其中：)sincos(2ijijijijjiijiijBGVVGVP根据以上模型可以形成式（4-30）的修正方程。构成该方程式包括形成等式左边的系数矩阵和等式右边的常数项两部分。（1）形成系数矩阵式（4-30）中修正方程的系数矩阵主要由4大部分组成：等式约束雅可比矩阵)(xhx、不等式约束雅可比矩阵)(xgx、对角矩阵ZL1，WU1和海森矩阵H。以下分别进行讨论。等式约束的雅可比矩阵：1014~)(xhQhPhxhRGx式中右端矩阵包含3个子矩阵：1022525212115151111......GGGGGGGGPQPPPQPPPQPPPQPPGPh其中：jijiPPPQGijGij0101022525212115151111......RRRRRRRRQPQPQPQPQQQPQQQPRQh其中：jijiQQQPRijRij010式中i为发电机的序号，j为节点号，ji表示第i台发电机是接在节点j上的，反之用ji表示。101055555151555551511515111115151111............~VQVPVQVPQPQPVQVPVQVPQPQPxh（潮流计算中的雅可比矩阵）不等式约束的雅可比矩阵：1414432143214321~~~~xgxgxgxgQgQgQgQgPgPgPgPggRRRRGGGGx式中1g、2g、3g和4g依次表示电源有功出力的上下界约束，无功电源无功出力的上下界约束，节点电压幅值的上下界约束和线路潮流约束。221IPgG，2210QgR，2101~0xg；2220PgG，222IQgR，2102~0xg；5230PgG，5230QgR，5240PgG，5240QgR510310...0000...0000...0100...00~xg（其中第i2行i列元素为1，其余元素均为0）5105455425415455425411451421411451421414............~VgVgVggggVgVgVggggxg矩阵中的元素为：)cossin(ijijijijjiiijBGVVp)cossin(ijijijijjijijBGVVp)sincos(2ijijijijjijiiijBGVGVVp)sincos(ijijijijijijBGVVp对角矩阵),...,(1414111lzlzdiagZL，),...,(1414111uwuwdiagWU。海森矩阵)(])[())(()()(222xgWUZLxgwzxgyxhxfH11Txxxxx。这是最复杂的部分，共包含4项。由以上推导已经可以得到其中第4项：)(])[(xgWUZLxg11Txx而其余3项是：目标函数的海森矩阵)(2xfx、等式约束海森矩阵与拉格朗日乘子y的乘积yxh)(2x和不等式约束海森矩阵与拉格朗日乘子wz的乘积))((2wzxgx，现分别讨论如下。①目标函数的海森矩阵：1414222)(00000000Axfx，其中2A是以机组燃料用二次系数ia2（GSi）为对角线的矩阵22R。②等式约束海森矩阵yxh)(2x与拉格朗日乘子y的乘积，yxh)(2x可表示为：niiiiiiiiiiiiiiiiiiRiiRiiiiiiiiiiiiiiGiiGiyQyPyQyPyQyPyQyPyQyPyQyPyQyPyQyPyQyP1222122222122221222212222212222212222122221222221222)~~()~~()~~()~~()()()~~()()(xxQxQxPxPxxQxQQQPQPQxPxPQPQPPPRRGGRRGRGRGGRGRG1414101021021010222221022222A00000000因此只需求其中)(25112QiiPiiiyyAAA，为此首先应求出piA和QiA：252552152152512512112212VPVPVVPVPVPPVPPiiiiiiiipiA根据iP的表达式（见模型）不难得到矩阵中的元素，如：)sincos(22ijijijijijjiiiBGVVP)sincos(2ijijijijjijiiBGVVP)cossin(2ijijijijijjiiiBGVVP)cossin(2ijijijijijiiBGVVP等等。同理，对于252552152152512512112212......VQVQVVQVQVQQVQQiiiiiiiiQiA也可根据iQ的表达式（见模型）不难得到矩阵中的元素，如：)cossin(22ijijijijijjiiiBGVVQ)cossin(2ijijijijjijiiBGVVQ)sincos(2ijijijijijjiiiBGVVQ)sincos(2ijijijijijiiBGVVQ等等。综合以上公式，即可得到A中各元素为：)]cossincos(sin)sincossin(cos[)(212212212212222125122jijjijiijiijijjijjijiijiijijijjijijjjijyyyyByyyyGVVyQyP)]sincossincos()cossincos(sin[)(2122122122122212512jijjijiijiijijjijjijiijiijijijjjiijjjiijyyyyByyyyGVyVQyVP)cossincossin()sincossincos([)(2122122122122212512jikjikiikiikikjikjikiikiikikkijkijjjkijyyyyByyyyGVVyQyP)sincossincos()cossincos(sin[)(2122122122122212512jikjikiikiikikjikjikiikiikikijkijjjkijyyyyByyyyGVyVQyVP)(2)(212222125122iiiiiijijjjijyByGyVQyVP)]sincossincos()cossincos(sin[)(2122122122122212512jikjikiikiikikjikijjikiikiikikjjiijjjiijyyyyByyyyGVyVQyVP)sincossin(cos)cossincossin([)(2122122122122212512jikjikiikiikikjikjikiikiikikjjkijjjkijyyyyByyyyGVyVQyVP)cossincos(sin)sincossin(cos[)(2122122122122212512jikjikiikiikikjikjikiikiikikjkijjjkijyyyyByyyyGyVVQyVVP③不等式约束海森矩阵与拉格朗日乘子wz的乘积))((2wzxg设cwz，我们有：