一种储能系统典型工况曲线的提取方法及其评价系统与流程

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

技术特征:

1.一种储能系统典型工况曲线的提取方法,其特征在于,所述方法包括如下步骤:

步骤1、读取新能源发电系统和储能系统的运行数据,所述运行数据包括新能源实际功率或储能系统充放电功率数据;

步骤2、确定或读取所述储能系统在一个采集时长的充放电功率数据,并组成充放电功率矩阵M;

步骤3、根据所述储能系统的功率值的区间分布特性,得到k时刻的特征充放电功率值向量;

步骤4、确定所述储能系统在k时刻的典型功率值:P1k

步骤5、剔除幅值较小的储能系统功率,确定所述储能系统在k时刻的典型功率值:P2k

步骤6、整合P1k和P2k,确定所述储能系统在k时刻的典型功率值:Pk

步骤7、分析并整合所有时刻的典型功率值,得到储能系统典型工况曲线。

步骤8、取得评价典型工况曲线。

2.根据权利要求1所述提取方法,其特征在于,所述步骤2包括如下步骤:

步骤2-1、当所述运行数据为新能源实际功率时,则采用基于斜率控制的平滑控制策略方法,确定所述储能系统的一个采集时长的充放电功率数据;

步骤2-2、当所述运行数据为储能系统充放电功率数据时,则直接读取所述充放电功率数据;

步骤2-3、将所述充放电功率数据组成充放电功率矩阵Mm×n

式中,m为所述一个采集时长的天数;n为采集时刻的总数;矩阵Mm×n中的元素pgk为第g天k时刻的所述储能系统的功率值,且k=1,2,…,n;g=1,2,…m。

3.根据权利要求2所述提取方法,其特征在于,所述步骤2-1包括如下步骤:

步骤2-1-1、根据所述新能源实际功率P(t),计算平滑前功率时间变化率rp(t)的:

<mrow> <msub> <mi>r</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>P</mi> <mi>h</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&Delta;</mi> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&Delta;</mi> <mi>t</mi> </mrow> </mfrac> </mrow>

式中,Ph(t)为平滑后功率;Δt为时间间隔;Ph(t-Δt)为t-Δt时刻平滑后功率;t为时刻;

步骤2-1-2、比较rp(t)与时间变化率临界值k时刻的大小,判断所述储能系统的充放电操作,并确定充放电功率。

4.根据权利要求3所述提取方法,其特征在于,所述步骤2-1-2包括如下步骤:

步骤2-1-2-1、当rp(t)>k时,所述储能系统充电,且所述储能系统在t时刻的充电功率为:

<mrow> <msubsup> <mi>P</mi> <mrow> <mi>b</mi> <mi>e</mi> <mi>s</mi> <mi>s</mi> </mrow> <mi>c</mi> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>&Delta;</mi> <mi>t</mi> <mo>&lsqb;</mo> <mi>k</mi> <mo>-</mo> <msub> <mi>r</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&rsqb;</mo> </mrow>

步骤2-1-2-2、当rp(t)<-k时,所述储能系统放电,且所述储能系统在t时刻的放电功率为:

<mrow> <msubsup> <mi>P</mi> <mrow> <mi>b</mi> <mi>e</mi> <mi>s</mi> <mi>s</mi> </mrow> <mrow> <mi>d</mi> <mi>i</mi> <mi>s</mi> <mi>c</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>&Delta;</mi> <mi>t</mi> <mo>&lsqb;</mo> <mo>-</mo> <mi>k</mi> <mo>-</mo> <msub> <mi>r</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&rsqb;</mo> </mrow>

步骤2-1-2-3、当-k<rp(t)<k时,所述储能系统不出力,并减少充放电次数。

5.根据权利要求1所述提取方法,其特征在于,所述步骤3包括如下步骤:

步骤3-1、提取m天内k时刻的各所述储能系统充放电功率向量p1k,p2k,…,pmk,组成向量Pk1=[p1k,p2k,…,pmk]T

步骤3-2、从向量Pk1中确定pgk分布最集中的区间,并将该区间作为最大概率分布区间Lopt1

步骤3-3、提取向量Pk1中落在所述最大概率分布区间Lopt1内的功率值p’1k,p’2k,…,p’tk,组成向量Ik1=[p’1k,p’2k,…,p’tk]T,并计算分布概率Fk1

<mrow> <msub> <mi>F</mi> <mrow> <mi>k</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mi>t</mi> <mi>m</mi> </mfrac> </mrow>

式中,t为向量Ik1中功率值的点数;

步骤3-4、采用遗传算法,以选取特征功率Ek1到向量Ik1中其他元素的欧氏距离最短为原则,根据Ek1的取值区间[p’min,p’max],得到用于获取所述选取特征功率Ek1的寻优目标函数min T:

<mrow> <mi>min</mi> <mi> </mi> <mi>T</mi> <mo>=</mo> <mi>m</mi> <mi>i</mi> <mi>n</mi> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>t</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mrow> <mi>k</mi> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <msup> <mi>p</mi> <mo>&prime;</mo> </msup> <mrow> <mi>j</mi> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow>

式中,p’max为向量Ik1的最大值;p’min为向量Ik1的最小值;p'jk为向量Ik1中的各功率值,j为向量Ik1中的某功率值,且j=1、2,…,tp

得到所述选取特征功率Ek1

步骤3-5、从向量Pk1中减去向量Ik1中的功率值,得到新的向量Pk2,得到最大概率分布区间Lopt2和向量Ik2,并计算分布概率Fk2

步骤3-6、判断Fk2与阈值ε大小;

若Fk2小于阈值ε,则k=k+1,即进入寻找下一时刻的典型功率值;

若Fk2大于阈值ε,则返回步骤3-2确定特征功率值Ek2及所有特征功率值Ek,直到Fki小于阈值ε,得到k时刻的特征充放电功率值向量[Ek1,Ek2,…,Ekh]。

6.根据权利要求5所述方法,其特征在于,所述步骤3-2包括如下步骤:

步骤3-2-1、确定分布区间L:

L=[x,x+d]

式中,x为所述分布区间的下限,d为所述区间长度;

步骤3-2-2、采用遗传算法,根据x的取值区间[pmin,pmax],得到用于求得优化所述分布区间的下限x取值后的最优值xopt1的优化目标函数maxJ:

<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>max</mi> <mi> </mi> <mi>J</mi> <mo>=</mo> <mi>max</mi> <mi>&Sigma;</mi> <mi>f</mi> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mrow> <mi>g</mi> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>f</mi> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mrow> <mi>g</mi> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mn>1</mn> <mo>,</mo> <mi>x</mi> <mo>&le;</mo> <msub> <mi>p</mi> <mrow> <mi>g</mi> <mi>k</mi> </mrow> </msub> <mo>&le;</mo> <mi>x</mi> <mo>+</mo> <mi>d</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>0</mn> <mo>,</mo> <mi>e</mi> <mi>l</mi> <mi>s</mi> <mi>e</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> </mtable> </mfenced>

式中,pmax为向量Pk1最大值,pmin为向量Pk1最小值;

得到优化所述分布区间的下限x取值后的最优值xopt1

步骤3-2-3、根据所述最优值xopt1,确定所述最大概率分布区间Lopt1

Lopt1=[xopt1,xopt1+d]。

7.根据权利要求1所述方法,其特征在于,所述步骤4中,确定所述储能系统在k时刻的典型功率值P1k,公式如下:

<mrow> <msub> <mi>P</mi> <mrow> <mn>1</mn> <mi>k</mi> </mrow> </msub> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>h</mi> </munderover> <mfrac> <msub> <mi>F</mi> <mrow> <mi>k</mi> <mi>i</mi> </mrow> </msub> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>h</mi> </munderover> <msub> <mi>F</mi> <mrow> <mi>k</mi> <mi>i</mi> </mrow> </msub> </mrow> </mfrac> <mo>&CenterDot;</mo> <msub> <mi>E</mi> <mrow> <mi>k</mi> <mi>i</mi> </mrow> </msub> </mrow>

式中,i为特征充放电功率值的某点数,h为特征充放电功率值的点数的总数;Eki为第i个特征功率值;Fki为第i个分布概率。

8.根据权利要求1所述提取方法,其特征在于,所述步骤5包括如下步骤:

步骤5-1、计算向量Pk1中元素的概率分布;

步骤5-2、根据概率分布,得到储能功率为0的概率α;同时获得分布区间[-x,x],选取分布区间的原则为向量Pk1中元素落入该区间内的概率大于阈值β,且区间长度最短;

步骤5-3、从向量Pk1中剔除落入区间[-x,x]内的元素,组成新的向量Pk1

步骤5-4、根据所述步骤3到4的方法,确定所述储能系统在k时刻的典型功率值P2k

9.根据权利要求1所述方法,其特征在于,所述步骤6中,确定所述储能系统在k时刻的典型功率值Pk,通过产生0~1的随机数r,若r≤α,则Pk=0,若α<r≤β,则Pk=P1k,否则Pk=P2k

10.根据权利要求1所述提取方法,其特征在于,所述步骤7中,分析并依次计算各个时刻的所述储能系统的典型充放功率值P1,P2,…,Pn;得到所述储能系统典型工况曲线Ptypical,公式如下:

Ptypical=[P1,P2,…,Pn]。

11.根据权利要求1所述方法,其特征在于,所述步骤8包括如下步骤:

步骤8-1、为定量评价储能系统典型功率曲线与实际功率曲线的偏差,定义偏差评价指标IR

<mrow> <msub> <mi>I</mi> <mi>R</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>P</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> </mfrac> <msqrt> <mfrac> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>n</mi> <mn>1</mn> </msub> </munderover> <msup> <mrow> <mo>&lsqb;</mo> <msub> <mi>P</mi> <mrow> <mi>t</mi> <mi>y</mi> <mi>p</mi> <mi>i</mi> <mi>c</mi> <mi>a</mi> <mi>l</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>P</mi> <mo>(</mo> <mi>i</mi> <mo>)</mo> <mo>&rsqb;</mo> </mrow> <mn>2</mn> </msup> </mrow> <msub> <mi>n</mi> <mn>1</mn> </msub> </mfrac> </msqrt> </mrow>

式中,IR为储能系统典型功率曲线与实际功率曲线的均方差,Ptypical为储能系统典型功率曲线,P为储能系统实际功率曲线,Pmax为典型功率曲线最大功率值,n1为功率曲线点数;

步骤8-2、为定量评价典型功率曲线与实际功率曲线概率分布的偏差,定义偏差评价指标IF

<mrow> <msub> <mi>I</mi> <mi>F</mi> </msub> <mo>=</mo> <msup> <mi>d</mi> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </msup> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>n</mi> <mn>2</mn> </msub> </munderover> <mo>|</mo> <msub> <mi>f</mi> <mrow> <mi>t</mi> <mi>y</mi> <mi>p</mi> <mi>i</mi> <mi>c</mi> <mi>a</mi> <mi>l</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>|</mo> </mrow>

式中,IF为储能系统典型功率曲线功率分布的拟合曲线与实际功率曲线功率分布拟合曲线的偏差评价指标,ftypical是储能系统典型功率曲线功率分布的拟合曲线,f是储能系统实际功率曲线功率分布拟合曲线,n2是拟合点数,d为储能系统功率分布直方图组距。

12.一种储能系统典型工况曲线的评价系统,其特征在于,所述评价系统包括:所述系统包括相互通信的通讯模块、数据存储与管理模块、控制策略执行模块、典型工况曲线挖掘模块和典型工况曲线分析与评估模块;

所述通讯模块用数据通讯网络与监控平台及所述数据存储与管理模块通信;

所述数据存储与管理模块向所述控制策略执行模块发送信号及数据;

所述控制策略执行模块向所述数据存储与管理模块及所述典型工况曲线挖掘模块发送数据信号;

所述典型工况曲线挖掘模块生成储能系统典型工况曲线,并向所述数据存储与管理模块和典型工况曲线分析与评估模块发送数据信号;

所述典型工况曲线分析与评估模块用于评估典型工况曲线的偏差。

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