一种含风电场的电力系统调度优化方法与流程

文档序号:12602119阅读:来源:国知局

技术特征:

1.一种含风电场的电力系统调度优化方法,

将停电损失引入到目标函数,通过最优模型求解出最优的机组组合及最小的成本,具体模型如下:

<mrow> <msub> <mi>F</mi> <mn>1</mn> </msub> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>T</mi> </munderover> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mo>{</mo> <msub> <mi>a</mi> <mi>n</mi> </msub> <mo>+</mo> <msub> <mi>b</mi> <mi>n</mi> </msub> <msub> <mi>P</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>c</mi> <mi>n</mi> </msub> <msup> <msub> <mi>P</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>Q</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <msub> <mi>u</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mo>|</mo> <msub> <mi>e</mi> <mi>n</mi> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mo>&lsqb;</mo> <msub> <mi>f</mi> <mi>n</mi> </msub> <mrow> <mo>(</mo> <msubsup> <mi>P</mi> <mi>n</mi> <mi>min</mi> </msubsup> <mo>-</mo> <msub> <mi>P</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>&rsqb;</mo> <mo>|</mo> <mo>}</mo> <mo>+</mo> <mi>J</mi> <mn>1</mn> <mo>+</mo> <mi>J</mi> <mn>2</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>

<mrow> <mi>J</mi> <mn>1</mn> <mo>=</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msubsup> <mi>P</mi> <mi>t</mi> <mi>L</mi> </msubsup> <mi>L</mi> <mi>%</mi> <mo>+</mo> <msubsup> <mi>P</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>w</mi> </msubsup> <mi>w</mi> <mi>u</mi> <mi>%</mi> <mo>)</mo> </mrow> </mrow>

<mrow> <mi>J</mi> <mn>2</mn> <mo>=</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <msubsup> <mi>P</mi> <mi>t</mi> <mi>L</mi> </msubsup> <mi>L</mi> <mi>%</mi> <mo>+</mo> <mo>(</mo> <mrow> <msub> <mi>P</mi> <mi>w</mi> </msub> <mo>-</mo> <msubsup> <mi>P</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>w</mi> </msubsup> </mrow> <mo>)</mo> <mi>w</mi> <mi>d</mi> <mi>%</mi> <mo>)</mo> </mrow> </mrow>

上式中J1,J2分别表示正负备用容量的损失成本,q1,q2分别为正负备用容量损失成本系数,(7)式中F1为加入停运损失后的总发电费用,

式中,T为时段数,t为时段,N为发电机台数,an,bn,cn分别表示第n台发电机组燃料成本系数,Pn,t为第n台发电机组第t时段输出的有功功率,Qn,t为机组n的启停费用,un,t为机组n在t时段的启停状态,0表示停机,1表示开机,en,fn为阀点效应系数,表示t时段机组n发电的最小出力,

对于式(1)有如下约束条件:

a)有功功率平衡约束。

<mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>u</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <msub> <mi>P</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>+</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>w</mi> </msub> </munderover> <msubsup> <mi>P</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>w</mi> </msubsup> <mo>=</mo> <msubsup> <mi>P</mi> <mi>t</mi> <mi>L</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>

式中为风电场k在第t段时段的有功功率输出,Nw为风电场总数,为第t时段的负荷总需求,负荷需求由常规机组出力和加上新能源风电场的功率输出,

b)常规机组发电出力约束

<mrow> <msubsup> <mi>P</mi> <mi>n</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msubsup> <mo>&le;</mo> <msub> <mi>P</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>&le;</mo> <msubsup> <mi>P</mi> <mi>n</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>

式中,分别为第n发电机组最小出力和最大出力,

c)发电机爬坡率约束

URn≤Pn,t-Pn,t-1≤DRn (4)

式中URn、DRn分别为发电机下调和上调爬坡速率,

d)系统正负旋转备用容量约束

<mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>P</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>max</mi> </msubsup> <msub> <mi>u</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>&GreaterEqual;</mo> <msubsup> <mi>P</mi> <mi>t</mi> <mi>L</mi> </msubsup> <mo>+</mo> <msubsup> <mi>P</mi> <mi>t</mi> <mi>L</mi> </msubsup> <mi>L</mi> <mi>%</mi> <mo>+</mo> <msubsup> <mi>P</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>w</mi> </msubsup> <mi>w</mi> <mi>u</mi> <mi>%</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>

负备用容量约束条件:

<mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>P</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>min</mi> </msubsup> <msub> <mi>u</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>&le;</mo> <msubsup> <mi>P</mi> <mi>t</mi> <mi>L</mi> </msubsup> <mo>-</mo> <msubsup> <mi>P</mi> <mi>t</mi> <mi>L</mi> </msubsup> <mi>L</mi> <mi>%</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mi>w</mi> </msub> <mo>-</mo> <msubsup> <mi>P</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>w</mi> </msubsup> <mo>)</mo> </mrow> <mi>w</mi> <mi>d</mi> <mi>%</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>

式(5)、(6)中分别表示t时段机组n发电的最大和最小出力,L%为系统总负荷对正负备用需求系数,wu%为风电出力预测误差对正备用容量的需求系数,wd%为风电出力预测出力对负备用容量的需求系数,Pw为风电场的额定装机容量总和,

采用粒子群算法,假设在一个D维的目标搜索空间中,有m个粒子组成一个群体,其中第i个粒子(i=1,2,...,m)位置表示为即第i个粒子在D维搜索空间中的位置是Xi,也就是说,每一个粒子的位置就是一个潜在解,将Xi代入目标函数就可以计算出其适应值,根据适应值的大小衡量其优劣,

粒子个体经历过的最好位置记为整个群体所有粒子经历过的最好位置记为粒子i的速度记为粒子群算法采用下列公式对粒子所在的位置不断更新:

<mrow> <msubsup> <mi>v</mi> <mi>i</mi> <mrow> <mi>d</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <msubsup> <mi>wv</mi> <mi>i</mi> <mi>d</mi> </msubsup> <mo>+</mo> <msub> <mi>&delta;</mi> <mn>1</mn> </msub> <msub> <mi>c</mi> <mn>1</mn> </msub> <msub> <mi>r</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>p</mi> <mi>i</mi> <mi>d</mi> </msubsup> <mo>-</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mi>d</mi> </msubsup> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&delta;</mi> <mn>2</mn> </msub> <msub> <mi>c</mi> <mn>2</mn> </msub> <msub> <mi>r</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>p</mi> <mi>g</mi> <mi>d</mi> </msubsup> <mo>-</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mi>d</mi> </msubsup> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>3</mn> </msub> <msub> <mi>r</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mrow> <msup> <mi>h</mi> <mi>d</mi> </msup> <mo>-</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mi>d</mi> </msubsup> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>

<mrow> <msubsup> <mi>x</mi> <mi>i</mi> <mrow> <mi>d</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mi>d</mi> </msubsup> <mo>+</mo> <msubsup> <mi>&alpha;v</mi> <mi>i</mi> <mrow> <mi>d</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>

公式中hd为第d次迭代中的最优解,δ1、δ2为影响因子,r3是[0,1]区间的随机数,

同时加速系数根据如下公式更新:

<mrow> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>c</mi> <mrow> <mn>1</mn> <mi>min</mi> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <msub> <mi>c</mi> <mrow> <mn>1</mn> <mi>max</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>c</mi> <mrow> <mn>1</mn> <mi>min</mi> </mrow> </msub> </mrow> <msub> <mi>d</mi> <mi>max</mi> </msub> </mfrac> <mi>d</mi> <mo>,</mo> <msub> <mi>c</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>c</mi> <mrow> <mn>2</mn> <mi>min</mi> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <msub> <mi>c</mi> <mrow> <mn>2</mn> <mi>max</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>c</mi> <mrow> <mn>2</mn> <mi>min</mi> </mrow> </msub> </mrow> <msub> <mi>d</mi> <mi>max</mi> </msub> </mfrac> <mi>d</mi> <mo>,</mo> <msub> <mi>c</mi> <mn>3</mn> </msub> <mo>=</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mrow> <mo>-</mo> <msub> <mi>c</mi> <mn>2</mn> </msub> <mi>d</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> </mrow>

其中c1min=c2min=0.5,c1max=c2max=2,dmax为最大迭代次数,

在搜索一开始,选择较大的c1和较小的c2,随着迭代次数的增加,c1线性递减,c2线性递增,使得接下去c1变小,c2变大,从而提高了全局最优值的收敛能力,

其特征在于,含风电场的电力系统调度优化方法包括以下步骤:

S101,种群初始化;

S102,参数代入经济调度模型公式(7);

S103,计算每个粒子的适应度值更新每个粒子的最佳适应度值、全局最佳值和所有最佳值中的最优值;

S104,由公式(10)、公式(11)更新粒子位置坐标和速度;

S105,判断是否达到最大迭代次数,如果是,则转S104,如果否,则执行下一步骤;

S106,获得最佳模型。

当前第2页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1