发明名称 基于pair-copula函数的微电网可靠性评估数据抽样方法
摘要 本发明公开了一种基于pair-copula函数方法的微电网可靠性评估数据抽样方法,利用Spearsman秩参数方法对微电网内采集到的风速、光照强度和负荷数据进行统计分析处理,得到微电网内风速、光照强度和负荷统计数据的随机相关性参数,将随机相关性参数带入基于pair-copula函数的微电网可靠性评估算法风速、光照强度和负荷数据的抽样过程中,由函数F(x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>)抽样风速、光照强度和负荷的小时数据,用于微电网可靠性评估算法中,由风速与风机出力关系、光照强度与光伏出力关系计算微电网中风机和光伏设备的实时出力,得到基于pair-copula函数的抽样数据。本发明解决了传统概率模型中无法考虑风速、光照强度和负荷数据随机相关性的问题;提高了微网可靠性评估算法中抽样数据的准确性。
申请公布号 CN104504618A 申请公布日期 2015.04.08
申请号 CN201410848222.9 申请日期 2014.12.29
申请人 天津大学 发明人 王守相;张兴友;葛磊蛟
分类号 G06Q50/06(2012.01)I 主分类号 G06Q50/06(2012.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 李素兰
主权项 一种基于pair‑copula函数方法的微电网可靠性评估数据抽样方法,其特征在于,该方法包括以下步骤:步骤(1)、输入微电网内每隔1小时采集到的风速、光照强度、负荷的年度数据;步骤(2)、运用Spearman秩参数方法计算微电网内风速、光照强度和负荷的随机秩相关性参数,利用风速与光照强度之间、风速与负荷之间、光照强度与负荷统计数据之间的Spearman秩参数ρ<sub>1</sub>、ρ<sub>2</sub>、ρ<sub>3</sub>,得到随机相关性参数ρ为ρ=(ρ<sub>1</sub>,ρ<sub>2</sub>,ρ<sub>3</sub>);步骤(3)、根据风速、光照强度和负荷的统计特性,计算风速、光照强度和负荷统计数据的边缘分布;步骤(4)、将Spearsman秩参数方法得出的随机相关性参数ρ代入pair‑copula函数,得到微电网内风速、光照强度和负荷的随机相关性表达式:F(x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>)=C(F<sub>1</sub>(x<sub>1</sub>),F<sub>2</sub>(x<sub>2</sub>),F<sub>3</sub>(x<sub>3</sub>),ρ)            (1)式中:x<sub>1</sub>、x<sub>2</sub>、x<sub>3</sub>分别代表风速、光照强度和负荷变量;F<sub>1</sub>(x<sub>1</sub>)、F<sub>2</sub>(x<sub>2</sub>)、F<sub>3</sub>(x<sub>3</sub>)分别为变量x<sub>1</sub>、x<sub>2</sub>、x<sub>3</sub>的边缘分布;ρ为步骤(2)中计算得到的Spearsman秩相关参数;C为pair‑copula函数;步骤(5)、由函数F(x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>)抽样风速、光照强度和负荷的小时数据,用于微电网可靠性评估算法中,由风速与风机出力关系、光照强度与光伏出力关系计算微电网中风机和光伏设备的实时出力;抽样过程包括以下步骤:在区间[0,1]产生均匀分布随机数w<sub>1</sub>、w<sub>2</sub>、w<sub>3</sub>;令x<sub>1</sub>=w<sub>1</sub>;将w<sub>2</sub>带入公式x<sub>2</sub>=h<sup>‑1</sup>(w<sub>2</sub>,x<sub>1</sub>,ρ<sub>1</sub>),得到x<sub>2</sub>;将w<sub>2</sub>、w<sub>3</sub>带入公式x<sub>3</sub>=h<sup>‑1</sup>[h<sup>‑1</sup>(w<sub>3</sub>,x<sub>2</sub>,ρ<sub>3</sub>),x<sub>1</sub>,p<sub>2</sub>],得到x<sub>3</sub>;(x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>)即为满足要求的随机数,代入风速、光照强度、负荷的边缘分布函数F<sub>1</sub>(x<sub>1</sub>)、F<sub>2</sub>(x<sub>2</sub>)、F<sub>3</sub>(x<sub>3</sub>)即得到风速、光照强度和负荷的抽样值;函数F(x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>)由其反函数h<sub>‑</sub><sup>1</sup>(x,y,ρ)计算得到:<img file="FDA0000650710310000021.GIF" wi="908" he="101" />其中:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mi>&phi;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mn>2</mn></msqrt><msup><mi>erf</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mn>2</mn><mi>x</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000650710310000022.GIF" wi="526" he="85" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>erf</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>sgn</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msqrt><msqrt><mrow><mo>(</mo><mfrac><mn>2</mn><mrow><mi>&pi;</mi><mo>&CenterDot;</mo><mi>a</mi></mrow></mfrac><mo>+</mo><mfrac><mrow><mi>ln</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>x</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><msup><mo>)</mo><mn>2</mn></msup><mo>-</mo><mfrac><mrow><mi>ln</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>x</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow><mi>a</mi></mfrac></mrow></msqrt><mo>-</mo><mrow><mo>(</mo><mfrac><mn>2</mn><mrow><mi>&pi;</mi><mo>&CenterDot;</mo><mi>a</mi></mrow></mfrac><mo>+</mo><mfrac><mrow><mi>ln</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>x</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>)</mo></mrow></msqrt><mo>;</mo></mrow>]]></math><img file="FDA0000650710310000023.GIF" wi="1499" he="172" /></maths>a=0.147;至此得到基于pair‑copula函数的抽样数据,用于计算微电网的可靠性。
地址 300072 天津市南开区卫津路92号