|
|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
含风电场的潮流计算,我用c++编程3 ]# }8 J" f2 q* g* d/ J. H
现在程序编好了,关键是处理风电场对应雅克比矩阵元素修改
$ Y( ]1 q3 V- j% r4 y, P% o. X/ ~" z, w+ M
我是这样处理风电场节点的
" P% n/ J! l0 |; u) J- {5 P2 c
) [4 e6 p3 J- M G I: l看做PQ节点,动态修改无功,修改雅克比矩阵对应元素& F7 L- b7 ?' k, u
8 G+ l, K0 g6 d% {: f$ M
把普通异步感应发电机当成PQ节点$ e4 Y" A" j: }; B8 m4 Q/ m
' G8 @) M4 J( m+ g# Z
这就需要修改雅克比矩阵 来完成含风电场的电力系统潮流计算。1 V& ?- E. n; c* w
& S/ u# _% H% k
首先其实有多种算法来计算潮流分布
+ u5 K) r* d9 m' n" P) _$ W1 s7 q1不修改雅克比矩阵
" J% G1 U W- N. L$ Z- ~5 K就是把风机等值电路画出来 并简化 步骤是这样的6 N1 X* o; O, e5 _- k
(1)给定风速,求有功p 再设定风电场节点电压u
: \' m0 {8 C% E" a3 n7 v8 ^(2)根据等效简化电路 导出无功与有功和电压的关系 进而可以求出无功
6 ?' `* f1 y3 ?' d(3)再把风电场节点看做PQ节点 求解雅克比矩阵 完成潮流计算# n' b; `; k: |$ B
( [1 t. M/ v" {/ r& b1 ?8 m* ^; _; H7 [2是修改雅克比矩阵$ |9 M; ?. z" v0 ^8 j [, [3 ?
由于无功是电压和有功的函数% q+ N1 b* ?4 e
对其求导可以得出dQ/dU % b: ?% X) j* p/ U
(1)设各个节点电压初值
9 i9 j& A% J `# r4 R: t5 ]: E3 n6 c(2)给定风速,从而知道有功,进而也可以求出无功(和第一种方法一样,导出无功和电压有功的关系式)( _1 T, Q8 D( }* e4 F6 S' X% K
(3)修改雅克比矩阵,对风电场节点的雅克比矩阵元素按照dQ/dU来修正,意思好像是说,原来的L加上上面求导的dQ/dU , Q/ S" V. ]& o% R$ u) \: h
(3)求解方程 修改节点电压和相角) R8 b, f) p6 [$ a3 t8 j
5 C, M1 g1 U& _0 {6 r. ~! q: E
我现在不明白的是 同样都是处理成PQ节点,为什么一个修改了矩阵 另一个却不修改矩阵 6 d9 |: y' a6 k1 L- D3 ?& r
为什么修改矩阵的这方法 修改元素要加上dQ/dU 而不是替换成dQ/dU 应该是有什么公式吧?
* z" y' r5 n6 E' v1 p# o8 Y修改对应元素的时候,
) A3 [: ~- p1 Y" A: k, Y$ ~4 m: }& o/ o. C/ Q$ V. n/ ]+ u
是这样修改吗?YKB[i+countpq+1][k+count1+1]=LL[i][k]+dQ/dU *u[k]5 d! F! l/ x6 p/ u* }, j" K* t
k=i 为风电场节点编号 countpq是总的pq节点个数! d7 J+ |8 F6 L# f
只是修改这一个元素 对不对啊$ D5 X: }) V4 U: O' T, v9 t
谢谢2 z1 ^- h; {1 J1 v
谢谢了 |
|