|
|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
含风电场的潮流计算,我用c++编程
+ O8 L# g; i0 k4 ^5 M现在程序编好了,关键是处理风电场对应雅克比矩阵元素修改5 F6 f( K/ S# G/ [! s" a: H
f9 i# U! V6 u" E我是这样处理风电场节点的
: N# `/ M2 Z) S' _9 g5 x v# z5 P
看做PQ节点,动态修改无功,修改雅克比矩阵对应元素; A- Y6 ^3 a/ j2 a& v
: g, L: u1 u8 ]# w/ z
把普通异步感应发电机当成PQ节点' l3 p8 j6 s2 T$ ~$ b
) G. S9 Q- x# [$ Y$ b这就需要修改雅克比矩阵 来完成含风电场的电力系统潮流计算。" I- ]# X+ F( S" Q) k- f$ k6 N7 B
# B8 _% M. d3 j& _
首先其实有多种算法来计算潮流分布9 C8 h) W7 K) t9 {5 l
1不修改雅克比矩阵 R7 e+ o4 ~! \4 |6 k
就是把风机等值电路画出来 并简化 步骤是这样的% n. L- S/ y# |. Y% M: Y+ X) S6 n
(1)给定风速,求有功p 再设定风电场节点电压u
7 D" N- G, z( g% K(2)根据等效简化电路 导出无功与有功和电压的关系 进而可以求出无功
5 A- c# ^: n* j q$ o! [(3)再把风电场节点看做PQ节点 求解雅克比矩阵 完成潮流计算 ?8 Q- B9 g% |
: A. J+ m! J7 @
2是修改雅克比矩阵3 f) n) w3 R O, ?2 K6 R0 }* ]
由于无功是电压和有功的函数+ E% D+ @6 r' L( P- G1 z3 \
对其求导可以得出dQ/dU
; ^! R# q4 I6 |* F2 t8 d* i! Z(1)设各个节点电压初值5 L2 V7 s2 a9 A6 v4 W$ k1 D: ~2 Y
(2)给定风速,从而知道有功,进而也可以求出无功(和第一种方法一样,导出无功和电压有功的关系式)
: g/ \! g& N j- I" G(3)修改雅克比矩阵,对风电场节点的雅克比矩阵元素按照dQ/dU来修正,意思好像是说,原来的L加上上面求导的dQ/dU
! G6 ]8 i6 m4 h(3)求解方程 修改节点电压和相角0 h0 Z1 I" g% q3 I
. E5 L6 {. c% j T* |我现在不明白的是 同样都是处理成PQ节点,为什么一个修改了矩阵 另一个却不修改矩阵
$ o3 i8 B8 B: U8 w! F& w" T为什么修改矩阵的这方法 修改元素要加上dQ/dU 而不是替换成dQ/dU 应该是有什么公式吧?
- b5 {# H% z+ @8 _9 s修改对应元素的时候,
! J6 I6 X5 q5 G k# n: E- u% R5 B" S6 I8 i" n8 t
是这样修改吗?YKB[i+countpq+1][k+count1+1]=LL[i][k]+dQ/dU *u[k]( j/ s. ^9 H: K$ j z# U
k=i 为风电场节点编号 countpq是总的pq节点个数1 X j4 k+ i% _8 n0 r
只是修改这一个元素 对不对啊) a9 _+ V1 I+ m7 D+ O& W
谢谢
* J1 b7 A) A; H# c谢谢了 |
|