主权项 |
一种基于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>φ</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>·</mo><msqrt><msqrt><mrow><mo>(</mo><mfrac><mn>2</mn><mrow><mi>π</mi><mo>·</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>π</mi><mo>·</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函数的抽样数据,用于计算微电网的可靠性。 |