|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
含风电场的潮流计算,我用c++编程
2 [- Y4 U9 J) @- H+ v, g) Z现在程序编好了,关键是处理风电场对应雅克比矩阵元素修改; `4 H9 ?- h& C2 Y; H
( e/ t+ {. ~+ i, l! F2 ~8 K5 W$ w
我是这样处理风电场节点的
' a& o& g$ `3 a5 [% M& q) s9 G
4 T8 f5 X# k: j7 ^, b8 x/ h看做PQ节点,动态修改无功,修改雅克比矩阵对应元素
* X: `$ ~# H6 B. R6 x
( ]- y$ B2 M' q8 |, H' t+ E把普通异步感应发电机当成PQ节点' u* {7 J: ]( @3 z0 _- i7 v7 C
) J; S2 @1 ]) v5 L
这就需要修改雅克比矩阵 来完成含风电场的电力系统潮流计算。' |/ ~4 h$ g9 `/ C" F( T6 F. \
& q0 ~- r4 Q7 d5 M2 r首先其实有多种算法来计算潮流分布
1 g" R; `! ?$ O' `2 q+ M' v& x* Q1不修改雅克比矩阵- w% e" v* r$ Y* d) v4 t# {0 I9 `. I
就是把风机等值电路画出来 并简化 步骤是这样的
9 f7 Z) q- f1 `0 \$ |6 x' ^" t(1)给定风速,求有功p 再设定风电场节点电压u$ {. Y4 p- H3 Y6 ~- _
(2)根据等效简化电路 导出无功与有功和电压的关系 进而可以求出无功
( J. [3 D/ H3 V, r' I% X(3)再把风电场节点看做PQ节点 求解雅克比矩阵 完成潮流计算& {& w5 Z" s: P7 a6 t9 X. A
" _" S$ N) [9 y! ?9 `9 p2是修改雅克比矩阵
; R) v$ j& `; e" \' _$ ~由于无功是电压和有功的函数
: _6 m% C% c2 y9 Z0 V- P: Q, F对其求导可以得出dQ/dU + N5 W& R4 D$ h. Y! B& M `
(1)设各个节点电压初值4 I) ?4 ?" S5 s) N8 N
(2)给定风速,从而知道有功,进而也可以求出无功(和第一种方法一样,导出无功和电压有功的关系式)% o! \6 `" T) |* N. a" v+ B- Q
(3)修改雅克比矩阵,对风电场节点的雅克比矩阵元素按照dQ/dU来修正,意思好像是说,原来的L加上上面求导的dQ/dU
7 G% `% z( }' r% [. [: F(3)求解方程 修改节点电压和相角& N& `/ c" H# A
9 i! A; ^# I: J( K! Q4 }我现在不明白的是 同样都是处理成PQ节点,为什么一个修改了矩阵 另一个却不修改矩阵 : z) q+ S( p+ ?& R
为什么修改矩阵的这方法 修改元素要加上dQ/dU 而不是替换成dQ/dU 应该是有什么公式吧?: v' `; e8 L3 C6 {3 m/ X
修改对应元素的时候, T1 Y) o- o3 b2 y: r+ u% f
. S$ @% h, r" X是这样修改吗?YKB[i+countpq+1][k+count1+1]=LL[i][k]+dQ/dU *u[k]
6 r/ P# q% f0 E: ]2 Z- B5 `" [k=i 为风电场节点编号 countpq是总的pq节点个数7 M: F# B( ]- g; i# Q
只是修改这一个元素 对不对啊
& `; Z2 F7 A X2 C/ A7 u }5 f谢谢
0 V3 i$ U/ u- j7 T0 H7 V. u谢谢了 |
|