CN106097141A - 光伏发电系统置信容量的三元区间数评估方法 - Google Patents

Abstract

本申请公开一种光伏发电系统置信容量的三元区间数评估方法，首先建立常规发电机组、光伏发电系统及负荷的三元区间数状态及对应的三元区间数概率向量；并得到常规发电机组、光伏发电系统及负荷的三元区间数通用生成函数模型；以电力不足期望值为电力系统可靠性指标分别计算不含光伏的发电系统可靠性和加入光伏后发电系统可靠性；利用基于三元区间数的布谷鸟搜索算法计算光伏发电系统的置信容量。本发明克服了普通点值评估及二元区间评估的不足，考虑了光伏出力和负荷的不确定性、随机性对评估结果的影响，在给出评估结果范围的同时也指明最大可能取值，便于生产决策及规划设计。

Description

光伏发电系统置信容量的三元区间数评估方法

技术领域

本发明涉及光伏发电系统置信容量评估领域，更具体地说是涉及一种考虑光伏发电系统、负荷的不确定性及随机性的光伏发电系统置信容量三元区间数的评估方法。

背景技术

现有针对光伏发电系统置信容量的评估方法中，求得的光伏发电系统置信容量均为一个点值，如中文文献《光伏发电系统可靠性分析及其置信容量计算》(王秀丽,武泽辰,曲翀.中国电机工程学报,2014,34(1):15-21.)中，利用分层抽样的序贯蒙特卡洛方法进行光伏发电系统置信容量的评估，且研究了天气变化、设备故障等随机因素对光伏发电系统置信容量计算的影响；外文文献《Madaeni S H,Sioshansi R,Denholm P.Comparingcapacity value estimation techniques for photovoltaic solar power[J],IEEEJournal of Photovoltaics,2013,3(1):407-415.》(Madaeni S.H.，Sioshansi R.,和Denholm P.(“光伏发电系统的容量价值评估方法比较分析”，IEEE光伏期刊，2013年第3卷第1期，ISNN 21563381,第407-415页)总结了光伏发电系统置信容量计算的各种近似解析法，其中包括基于光伏发电容量因子的近似方法、Garver近似计算方法、光伏发电系统多状态机组等效方法及Z-method方法等，这些近似方法采用解析法去逼近仿真算法得到的精确值，其在计算光伏发电系统置信容量中，只能进行简单的解析运算。

中国发明专利《一种确定光伏发电容量可信度的方法》(CN103218757B)建立了分时段的光伏发电系统的置信容量评估模型，能够计算白天和全时段的光伏发电置信容量。上述各种方法虽然均能计算光伏发电系统的置信容量，但是计算等到的点值却并非能表示光伏置信容量的精确值。实际中，由于光伏发电系统出力的波动性、间歇性及负荷需求的不确定性，均会导致光伏发电系统置信容量在一定范围内波动，故上述计算得到的点值只是光伏发电系统置信容量可能的一个取值。

综上所述，现有技术中还存在以下不足：

1.基于分层抽样的序贯蒙特卡洛方法计算时间长、成本大；该方法虽然能考虑光伏发电系统出力及负荷的不确定性，但其最终求得的光伏发电系统置信容量只是考虑随机性后的一个期望值，并不能给出光伏发电系统置信容量的波动范围及可能的最大取值，所求的期望值对系统规划及后期调度运行的指导作用有限；

2.所述的近似解析法评估虽然计算量小，但其精确度低，且随着光伏安装容量、负荷的变化不同，其评估结果差异较大，使得该方法对在实际中将产生较大不确定性，且该方法鲁棒性较差。

因此，为了更合理的表示光伏发电系统置信容量的波动性及可能的取值范围，目前需要提出一种表征光伏发电系统置信容量波动范围的三元区间数评估方法，不仅能衡量光伏发电系统出力及负荷的不确定性对光伏发电系统置信容量的影响程度，同时也能指出光伏发电系统置信容量最可能取值。在此基础上提出基于三元区间数的通用生成函数用于计算光伏发电系统的可靠性指标具有建模简单、求解方便、计算量小的优点。

发明内容

针对现有技术的不足，本发明的目的在于，提供一种基于三元区间数的光伏发电系统置信容量评估方法，该方法能克服常规点值评估方法的不足。能描述光伏发电系统置信容量因光伏出力、负荷波动性和元件可靠性数据非精确性导致的波动范围，同时三元区间数不同于常规的二元区间数，其区间宽度能表示光伏发电系统置信容量波动的范围，三元区间数还能指定光伏发电系统置信容量最可能的取值。

本发明的目的是这样实现的，本发明提供了一种光伏发电系统置信容量的三元区间数评估方法，应用于包含光伏发电系统的电力系统，所述光伏发电系统的额定容量为C_PV，C_PV>0；所述三元区间数其中a^-、a⁺和a⁰分别为三元区间数的下边界、上边界和特元；所述方法包括以下步骤：

步骤1，通过采集历史数据获取包含光伏发电系统的电力系统中常规发电机组、光伏发电系统及负荷数据，包括：各常规发电机组的额定容量G_i、强迫停运率p_i，i＝1,2,…,N，N为常规发电机组数目；负荷一年8760小时的数据并按负荷值从大到小进行降序排列以得到负荷年小时需求序列L_s＝{L₁,L₂,…,L₈₇₆₀}；光伏发电系统一年8760小时的出力数据并按光伏出力从大到小进行降序排列以得到光伏发电系统年小时出力序列C_s＝{C₁,C₂,…,C₈₇₆₀}；

步骤2，将步骤1中获取的常规发电机组额定容量G_i、强迫停运率p_i、负荷年小时需求序列L_s及光伏发电系统年小时出力序列C_s分别从点值转化为三元区间数状态及对应的三元区间数概率的形式，具体按如下步骤进行：

步骤2.1，常规发电机组额定容量G_i及强迫停运率p_i从点值到三元区间数状态及对应的三元区间数概率的转化；

设第i台常规发电机组的额定容量G_i的三元区间数状态表示为G_i]，强迫停运率p_i的点值对应的三元区间数概率通过重复n次贝努利试验得到，随机生成“0”和“1”，若生成的数据为“0”，表示强迫停运率p_i为非精确值，若生成的数据为“1”，表示强迫停运率p_i为精确值，此时强迫停运率p_i的点值对应的三元区间数概率为：

$\[{\tilde{p}}_i\]=\[p_i-\frac{S}{\sqrt{n}}{t}_{\mathrm{\α}/2}(n-1),p_i,p_i+\frac{S}{\sqrt{n}}{t}_{\mathrm{\α}/2}(n-1)\]---\left(1\right)$

在公式(1)中，S为n次贝努利试验的标准差，由下式计算得到：

$S=\sqrt{\frac{m}{n}(1-\frac{m}{n})}---\left(2\right)$

在公式(2)中，m为n次贝努利试验中发生“1”的次数，α＝0.05，t_α/2(n-1)的数值由t分布统计表查询得到；

重复上述过程得到N台常规发电机组的三元区间数状态向量表示为：

$\{\[{\tilde{G}}_c\]\}=\{\[{\tilde{G}}_1\],\[{\tilde{G}}_2\],\mathrm{...},\[{\tilde{G}}_i\],\mathrm{...},\[{\tilde{G}}_N\]\}---\left(3\right)$

对应的三元区间数概率向量表示为：

$\{\[{\tilde{p}}_c\]\}=\{\[{\tilde{p}}_1\],\[{\tilde{p}}_2\],\mathrm{...},\[{\tilde{p}}_i\],\mathrm{...},\[{\tilde{p}}_N\]\}---\left(4\right)$

步骤2.2，负荷年小时需求序列L_s从点值到三元区间数状态及对应的三元区间数概率的转化；

利用等宽划分的方法将所述负荷年小时需求序列L_s等分为M_L个等宽区间，每一等宽区间均采用模糊C均值聚类方法计算该区间的聚类中心，则等宽区间的左端点、右端点和聚类中心分别对应该区间三元区间数的下边界、上边界和特元，得到负荷年小时需求序列L_s的三元区间数状态向量为：

$\{\[\tilde{L}\]\}=\{\[{\tilde{L}}_{e1}\],\[{\tilde{L}}_{e2}\],\mathrm{...},\[{\tilde{L}}_{ej}\],\mathrm{...},\[{\tilde{L}}_{{\mathrm{eM}}_L}\]\}---\left(5\right)$

对应的负荷年小时需求序列L_s的三元区间数概率向量为：

$\{\[{\tilde{p}}_L\]\}=\{\[{\tilde{p}}_{e1}\],\[{\tilde{p}}_{e2}\],\mathrm{...},\[{\tilde{p}}_{ej}\],\mathrm{...},\[{\tilde{p}}_{{\mathrm{eM}}_L}\]\},j=1,2,\mathrm{...},M_L---\left(6\right)$

在公式(5)、(6)中，表示第j个等宽区间的负荷三元区间数状态，表示对应的负荷三元区间数概率；

步骤2.3，光伏发电系统年小时出力序列C_s从点值到三元区间数状态及对应的三元区间数概率的转化；

利用等宽划分的方法将所述光伏发电系统年小时出力序列C_s等分为M_c个等宽区间，每一等宽区间均采用模糊C均值聚类方法计算该区间的聚类中心，则等宽区间的左端点、右端点和聚类中心分别对应该区间三元区间数的下边界、上边界和特元，得到光伏发电系统年小时出力序列C_s的三元区间数状态向量表示为：

$\{\[\tilde{V}\]\}=\{\[{\tilde{V}}_1\],\[{\tilde{V}}_2\],\mathrm{...},\[{\tilde{V}}_k\],\mathrm{...},\[{\tilde{V}}_{M_c}\]\}---\left(7\right)$

对应的光伏发电系统年小时出力序列的三元区间数概率向量表示为：

$\{\[{\tilde{p}}_V\]\}=\{\[{\tilde{p}}_{V_1}\],\[{\tilde{p}}_{V_2}\],\mathrm{...},\[{\tilde{p}}_{V_k}\],\mathrm{...},\[{\tilde{p}}_{V_{M_c}}\]\},k=1,2,\mathrm{...},M_c;---\left(8\right)$

在公式(7)、(8)中，表示第k个等宽区间的光伏发电系统的三元区间数状态，表示对应的光伏发电系统三元区间数概率；

步骤3，由步骤2得到的常规发电机组、负荷及光伏发电系统的三元区间数状态向量及对应的三元区间数概率向量建立其对应的三元区间数通用生成函数，具体步骤如下：

步骤3.1，由步骤2.1中式(3)和(4)得到第i台常规发电机组的三元区间数通用生成函数其表达式为：

${\tilde{U}}_i\left(z\right)=\[{\tilde{p}}_i\]\×z^{\[\tilde{0}\]}+(\[\tilde{1}\]-\[{\tilde{p}}_i\])\×z^{\[{\tilde{G}}_i\]},i=1,2,\mathrm{...},N---\left(9\right)$

则由N台常规发电机组构成的常规发电系统的三元区间数通用生成函数表示为：

$\begin{array}{c}{\tilde{U}}_G\left(z\right)=\⊕\left({\tilde{U}}_1\right(z),{\tilde{U}}_2(z),\mathrm{...},{\tilde{U}}_i(z),\mathrm{...},{\tilde{U}}_N(z\left)\right)\\ =\underset{m=1}{\overset{M_G}{\Σ}}\[{\tilde{p}}_m^G\]\×z^{\left({\tilde{G}}_m^G\right)}\end{array}---\left(10\right)$

在公式(9)、(10)中，算子表示两个及多个三元区间数通用生成函数的“和”运算，即将各个三元区间数状态相加，对应的三元区间数概率相乘；z为通用生成函数变换因子，M_G表示由N台常规发电机组构成的常规发电系统的总状态数，m表示所述常规发电系统的第m个状态，表示所述常规发电系统的第m个三元区间数状态，表示对应的三元区间数概率；

步骤3.2，根据步骤2.2中公式(5)和(6)得到负荷的三元区间数通用生成函数其表达式为：

${\tilde{U}}_L\left(z\right)=\underset{j=1}{\overset{M_L}{\Σ}}\[{\tilde{p}}_{ej}\]\×z^{\[{\tilde{L}}_{ej}\]}---\left(11\right)$

步骤3.3，根据步骤2.3中公式(7)和(8)得到光伏发电系统的三元区间数通用生成函数其表达式为：

${\tilde{U}}_{PV}\left(z\right)=\underset{k=1}{\overset{M_c}{\Σ}}\[{\tilde{p}}_{V_k}\]\×z^{\[{\tilde{V}}_k\]}---\left(12\right)$

步骤4，分别计算不含光伏发电系统及包含光伏发电系统的电力系统电力不足期望值可靠性指标即LOLE指标，具体计算步骤如下：

步骤4.1，按下式计算不含光伏发电系统的电力系统电力不足期望值可靠性指标

$\&lsqb;{\tilde{R}}_0\&rsqb;=8760\&times;\underset{j=1}{\overset{M_L}{\&Sigma;}}\underset{m=1}{\overset{M_G}{\&Sigma;}}p(\&lsqb;{\tilde{G}}_m^G\&rsqb;<\&lsqb;{\tilde{L}}_{ej}\&rsqb;)---\left(13\right)$

在公式(13)中，表示三元区间数小于三元区间数的概率，其按如下过程计算：

令由点(c^-,0),(c⁰,1)和(c⁺,0)构成的三角形面积为该三角形在平面直角坐标系中y轴左侧的面积记为S_left，则表示为：

$p(\&lsqb;{\tilde{G}}_m^G\&rsqb;<\&lsqb;{\tilde{L}}_{ej}\&rsqb;)=\frac{S_{left}}{S_{\&lsqb;\tilde{c}\&rsqb;}}---\left(14\right)$

步骤4.2，计算包含光伏发电系统后发电系统的电力不足期望值可靠性指标按如下步骤计算：

光伏发电系统及所述的常规发电系统构成的总发电系统的三元区间数通用生成函数如下：

${\tilde{U}}_E\left(z\right)={\tilde{U}}_G\left(z\right)\&CirclePlus;{\tilde{U}}_{PV}\left(z\right)=\underset{q=1}{\overset{M}{\&Sigma;}}\&lsqb;{\tilde{p}}_q^E\&rsqb;z^{\&lsqb;{\tilde{E}}_q\&rsqb;}---\left(15\right)$

在公式(15)中，M为所述的总发电系统的状态数，为所述的总发电系统的第q个三元区间数状态，为第q个三元区间数状态对应的三元区间数概率；令加入额定容量为C_PV的光伏发电系统后其可额外承担的负荷增加量为待求变量则包含光伏发电系统后发电系统的电力不足期望值可靠性指标表示如下：

$\&lsqb;{\tilde{R}}_{PV}\&rsqb;=8760\&times;\underset{j=1}{\overset{M_L}{\&Sigma;}}\underset{q=1}{\overset{M}{\&Sigma;}}p(\&lsqb;{\tilde{E}}_q\&rsqb;<\&lsqb;{\tilde{L}}_{ej}\&rsqb;+\&lsqb;\mathrm{\&Delta;}\tilde{L}\&rsqb;)---\left(16\right)$

步骤5，运用基于三元区间数的布谷鸟搜索算法求解式(16)中的待求变量使其满足下式:

$\&lsqb;{\tilde{R}}_0\&rsqb;=\&lsqb;{\tilde{R}}_{PV}\&rsqb;---\left(17\right)$

公式(17)得到的负荷增加量即为所述光伏发电系统置信容量的三元区间数评估结果，具体按以下步骤实现：

步骤5.1，在区间[0,C_PV]范围内随机生成w个三元区间数0≤S_i ^-≤S_i ⁰≤S_i ⁺≤C_PV，i＝1，2，...，w；

步骤5.2，令代入公式(16)中，计算在该负荷增加量下的电力系统LOLE指标计算与的相似度重复该过程计算w个三元区间数下的电力系统LOLE指标及其与的相似度，并记录相似度最大值，判断其是否大于预先给定的阈值，如果大于则停止运行，且此时所述光伏发电系统的置信容量即为相似度最大时对应的负荷增加量执行步骤5.4；否则，执行步骤5.3；

步骤5.3，对步骤5.1中所述的w个三元区间数均进行列维飞行生成新的三元区间数并令代入公式(16)中，计算在该负荷增加量下的电力系统的LOLE指标计算与的相似度比较F_i和F_i′的大小，如果F_i′>F_i，则用三元区间数替换否则，保留原三元区间数不变；针对w个新的三元区间数，记录相似度最大值，判断其是否大于预先给定的阈值，如果大于则停止运行，且此时所述光伏发电系统的置信容量即为相似度最大时对应的负荷增加量执行步骤5.4；否则，继续执行步骤5.3；

步骤5.4，记录步骤5.2或步骤5.3中得到的相似度最大时对应的负荷增加量即此时所述光伏发电系统的置信容量为

优选的，步骤5中所述的列维飞行实现过程如下：

设三元区间数初始位置为通过列维飞行后新的三元区间数更新为具体按如下进行：

$\&lsqb;{\tilde{S}}_i^{\&prime;}\&rsqb;=\&lsqb;{\tilde{S}}_i\&rsqb;+\&lsqb;\widetilde{\mathrm{\&alpha;}}\&rsqb;\&CircleTimes;\mathrm{\&gamma;}\left(\mathrm{\&epsiv;}\right)---\left(18\right)$

γ(ε)～|ε|^-λ,(1＜λ≤3) (19)

在公式(18)和(19)中，表示随机三元区间数步长，符号表示点乘计算，γ(ε)为列维随机飞行步长，ε为列维分布参数，λ为指数参数。

优选的，步骤5中所述的相似度按如下定义进行：

设两个三元区间数其相似度定义为：

$F_{\&lsqb;\tilde{A}\&rsqb;,\&lsqb;\hat{B}\&rsqb;}=\left\{\begin{array}{ccc}\frac{S_{\&lsqb;\tilde{A}\&rsqb;\&cap;\&lsqb;\tilde{B}\&rsqb;}}{\max (S_{\&lsqb;\tilde{A}\&rsqb;},S_{\&lsqb;\tilde{B}\&rsqb;})},& (S_{\&lsqb;\tilde{A}\&rsqb;\&cap;\&lsqb;\tilde{B}\&rsqb;}\&GreaterEqual;0)& \left(a\right)\\ -\frac{|b^{-}-a^{+}|}{|b^{+}-a^{-}|},& (b^{-}>a^{+})& \left(b\right)\\ -\frac{|a^{-}-b^{+}|}{|a^{+}-b^{-}|},& (a^{-}>b^{+})& \left(c\right)\end{array}\right.---\left(20\right)$

在公式(20)中，表示由点(a^-,0),(a⁰,1)和(a⁺,0)构成的三角形面积，表示由点(b^-,0),(b⁰,1)和(b⁺,0)构成的三角形面积，表示和的重叠部分面积。

优选的，步骤5中所述的预先给定的阈值的取值范围为(0.95，1.00)。

与已有技术相比，本发明的有益效果体现在：

1、本发明方法能考虑光伏发电系统出力及负荷不确定性，与已有技术相比较，其在光伏发电系统置信容量评估中具有更好的鲁棒性；

2、本发明采用了三元区间数对光伏发电系统的置信容量进行了评估，三元区间数不仅能给出评估结果的波动范围，特别的三元区间数的特元还能指明光伏发电系统置信容量最可能的取值，为含光伏发电系统的电力系统的规划提供价值参考；

3、本发明将常规的基于点值或二元区间数的通用生成函数扩展至三元区间数的通用生成函数，为考虑不确定性系统的可靠性评估提供了新的手段和方法；

4、本发明中采用基于三元区间数的布谷鸟搜索算法，实现了三元区间数的智能搜索算法，为其他三元区间数的问题提供了新的求解方法。

附图说明

图1为本发明光伏发电系统置信容量的三元区间数评估方法的流程图；

图2为本发明中基于三元区间数的布谷鸟搜索算法的执行过程图。

具体实施方式

下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。

参考图1，示出了本发明光伏发电系统置信容量的三元区间数评估方法的流程图，本发明所述方法应用于包含光伏发电系统的电力系统，所述光伏发电系统的额定容量为C_PV(C_PV>0)；所述的三元区间数指的是形如 的数，其中a^-,a⁺和a⁰分别为三元区间数的下边界、上边界和特元；所述方法包括以下步骤：

步骤1，通过采集历史数据获取包含光伏发电系统的电力系统中常规发电机组、光伏发电系统及负荷数据，包括：各常规发电机组的额定容量G_i、强迫停运率p_i(i＝1,2,…,N)，N为常规发电机组数目；负荷一年8760小时的数据并按负荷值从大到小进行降序排列以得到负荷年小时需求序列L_s＝{L₁,L₂,…,L₈₇₆₀}；光伏发电系统一年8760小时的出力数据并按光伏出力从大到小进行降序排列以得到光伏发电系统年小时出力序列C_s＝{C₁,C₂,…,C₈₇₆₀}。

对于部分常规发电机组其强迫停运率未给出，但是可由其平均正常工作时间和平均故障时间计算得到该常规发电机组的强迫停运率。如第i台常规发电机组的平均正常工作时间为t_{MTTF_i}，平均故障时间为t_{MTTR_i}，则第i台常规发电机组的故障率p_i按下式计算得到：

$p_i=\frac{t_{MTTR\_i}}{t_{MTTR\_i}+t_{MTTF\_i}}---\left(1\right)$

步骤2，将步骤1中获取的常规发电机组额定容量G_i、强迫停运率p_i、负荷年小时需求序列L_s及光伏发电系统年小时出力序列C_s分别从点值转化为三元区间数状态及对应的三元区间数概率的形式，具体按如下步骤进行：

步骤2.1，常规发电机组额定容量G_i及强迫停运率p_i从点值到三元区间数状态及对应的三元区间数概率的转化。

本申请中常规发电机组采用运行和停运两状态可靠性模型，即常规发电机组以额定功率运行或者处于完全故障状态，且对该两种状态其出力状态均为确定的数值，无波动范围，故设第i台常规发电机组的额定容量G_i的三元区间数状态表示为强迫停运率p_i的点值对应的三元区间数概率通过重复n次贝努利试验得到，随机生成“0”和“1”，若生成的数据为“0”表示强迫停运率p_i为非精确值，若生成的数据为“1”表示强迫停运率p_i为精确值，此时强迫停运率p_i的点值对应的三元区间数概率为：

$\&lsqb;{\tilde{p}}_i\&rsqb;=\&lsqb;p_i-\frac{S}{\sqrt{n}}{t}_{\mathrm{\&alpha;}/2}(n-1),p_i,p_i+\frac{S}{\sqrt{n}}{t}_{\mathrm{\&alpha;}/2}(n-1)\&rsqb;---\left(2\right)$

上式中，S为n次贝努利试验的标准差，由下式计算得到：

$S=\sqrt{\frac{m}{n}(1-\frac{m}{n})}---\left(3\right)$

上式中，m为n次贝努利试验中发生“1”的次数，α＝0.05，t_α/2(n-1)的数值由t分布统计表查询得到。

重复上述过程得到N台常规发电机组的三元区间数状态向量表示为：

$\{\&lsqb;{\tilde{G}}_c\&rsqb;\}=\{\&lsqb;{\tilde{G}}_1\&rsqb;,\&lsqb;{\tilde{G}}_2\&rsqb;,\mathrm{...},\&lsqb;{\tilde{G}}_i\&rsqb;,\mathrm{...},\&lsqb;{\tilde{G}}_N\&rsqb;\}---\left(4\right)$

对应的三元区间数概率向量表示为：

$\{\&lsqb;{\tilde{p}}_c\&rsqb;\}=\{\&lsqb;{\tilde{p}}_1\&rsqb;,\&lsqb;{\tilde{p}}_2\&rsqb;,\mathrm{...},\&lsqb;{\tilde{p}}_i\&rsqb;,\mathrm{...},\&lsqb;{\tilde{p}}_N\&rsqb;\};---\left(5\right)$

步骤2.2，负荷年小时需求序列L_s从点值到三元区间数状态及对应的三元区间数概率的转化。

利用等宽划分的方法将所述负荷年小时需求序列L_s等分为M_L个等宽区间，每一等宽区间均采用模糊C均值聚类方法计算该区间的聚类中心，则等宽区间的左端点、右端点和聚类中心分别对应该区间三元区间数的下边界、上边界和特元，得到负荷年小时需求序列L_s的三元区间数状态向量为：

$\{\&lsqb;\tilde{L}\&rsqb;\}=\{\&lsqb;{\tilde{L}}_{e1}\&rsqb;,\&lsqb;{\tilde{L}}_{e2}\&rsqb;,\mathrm{...},\&lsqb;{\tilde{L}}_{ej}\&rsqb;,\mathrm{...},\&lsqb;{\tilde{L}}_{{\mathrm{eM}}_L}\&rsqb;\}---\left(6\right)$

对应的负荷年小时需求序列的三元区间数概率向量为：

$\{\&lsqb;{\tilde{p}}_L\&rsqb;\}=\{\&lsqb;{\tilde{p}}_{e1}\&rsqb;,\&lsqb;{\tilde{p}}_{e2}\&rsqb;,\mathrm{...},\&lsqb;{\tilde{p}}_{ej}\&rsqb;,\mathrm{...},\&lsqb;{\tilde{p}}_{{\mathrm{eM}}_L}\&rsqb;\},(j=1,2,\mathrm{...},M_L)---\left(7\right)$

在公式(6)、(7)中，表示第j个等宽区间的负荷三元区间数状态，表示对应的负荷三元区间数概率。

步骤2.3，光伏发电系统年小时出力序列C_s从点值到三元区间数状态及对应的三元区间数概率的转化。

利用等宽划分的方法将所述光伏发电系统年小时出力序列C_s等分为M_c个等宽区间，每一等宽区间均采用模糊C均值聚类方法计算该区间的聚类中心，则等宽区间的左端点、右端点和聚类中心分别对应该区间三元区间数的下边界、上边界和特元，得到光伏发电系统年小时出力序列C_s的三元区间数状态向量表示为：

$\{\&lsqb;\tilde{V}\&rsqb;\}=\{\&lsqb;{\tilde{V}}_1\&rsqb;,\&lsqb;{\tilde{V}}_2\&rsqb;,\mathrm{...},\&lsqb;{\tilde{V}}_k\&rsqb;,\mathrm{...},\&lsqb;{\tilde{V}}_{M_c}\&rsqb;\}---\left(8\right)$

对应的光伏发电系统年小时出力序列的三元区间数概率向量表示为：

$\{\&lsqb;{\tilde{p}}_V\&rsqb;\}=\{\&lsqb;{\tilde{p}}_{V_1}\&rsqb;,\&lsqb;{\tilde{p}}_{V_2}\&rsqb;,\mathrm{...},\&lsqb;{\tilde{p}}_{V_k}\&rsqb;,\mathrm{....},\&lsqb;{\tilde{p}}_{V_{M_c}}\&rsqb;\},(k=1,2,\mathrm{...},M_c);---\left(9\right)$

在公式(8)、(9)中，表示第k个等宽区间的光伏发电系统的三元区间数状态，表示对应的光伏发电系统三元区间数概率。

步骤3，由步骤2得到的常规发电机组、负荷及光伏发电系统的三元区间数状态向量及对应的三元区间数概率向量建立其对应的三元区间数通用生成函数，具体步骤如下：

步骤3.1，由步骤2.1中式(4)和(5)得到第i台常规发电机组的三元区间数通用生成函数其表达式为：

${\tilde{U}}_i\left(z\right)=\&lsqb;{\tilde{p}}_i\&rsqb;\&times;z^{\&lsqb;\tilde{0}\&rsqb;}+(\&lsqb;\tilde{1}\&rsqb;-\&lsqb;{\tilde{p}}_i\&rsqb;)\&times;z^{\&lsqb;{\tilde{G}}_i\&rsqb;}---\left(10\right)$

则由N台常规发电机组构成的常规发电系统的三元区间数通用生成函数表示为：

$\begin{array}{c}{\tilde{U}}_G\left(z\right)=\&CirclePlus;\left({\tilde{U}}_1\right(z),{\tilde{U}}_2(z),\mathrm{...},{\tilde{U}}_i(z),\mathrm{...},{\tilde{U}}_N(z\left)\right)\\ =\underset{m=1}{\overset{M_G}{\&Sigma;}}\&lsqb;{\tilde{p}}_m^G\&rsqb;\&times;z^{\left({\tilde{G}}_m^G\right)}\end{array}---\left(11\right)$

在公式(10)、(11)中，算子表示两个及多个三元区间数通用生成函数的“和”运算，即将各个三元区间数状态相加，对应的三元区间数概率相乘；z为通用生成函数变换因子，M_G表示由N台常规发电机组构成的常规发电系统的总状态数，m表示所述常规发电系统的第m个状态，表示所述常规发电系统的第m个三元区间数状态，表示对应的三元区间数概率。

步骤3.2，根据步骤2.2中公式(6)和(7)得到负荷的三元区间数通用生成函数其表达式为：

${\tilde{U}}_L\left(z\right)=\underset{j=1}{\overset{M_L}{\&Sigma;}}\&lsqb;{\tilde{p}}_{ej}\&rsqb;\&times;z^{\&lsqb;{\tilde{L}}_{ej}\&rsqb;}---\left(12\right)$

步骤3.3，根据步骤2.3中公式(8)和(9)得到光伏发电系统的三元区间数通用生成函数其表达式为：

${\tilde{U}}_{PV}\left(z\right)=\underset{k=1}{\overset{M_c}{\&Sigma;}}\&lsqb;{\tilde{p}}_{V_k}\&rsqb;\&times;z^{\&lsqb;{\tilde{V}}_k\&rsqb;}---\left(12\right)$

步骤4，分别计算不含光伏发电系统及包含光伏发电系统的电力系统电力不足期望值可靠性指标(LOLE)，具体计算步骤如下：

步骤4.1，按下式计算计算不含光伏发电系统的电力系统电力不足期望值可靠性指标

$\&lsqb;{\tilde{R}}_0\&rsqb;=8760\&times;\underset{j=1}{\overset{M_L}{\&Sigma;}}\underset{m=1}{\overset{M_G}{\&Sigma;}}p(\&lsqb;{\tilde{G}}_m^G\&rsqb;<\&lsqb;{\tilde{L}}_{ej}\&rsqb;)---\left(14\right)$

在公式(14)中，表示三元区间数小于三元区间数的概率，其按如下过程计算：

令由点(c^-,0),(c⁰,1)和(c⁺,0)构成的三角形面积为该三角形在平面直角坐标系中y轴左侧的面积记为S_left，则表示为：

$p(\&lsqb;{\tilde{G}}_m^G\&rsqb;<\&lsqb;{\tilde{L}}_{ej}\&rsqb;)=\frac{S_{left}}{S_{\&lsqb;\tilde{c}\&rsqb;}}---\left(15\right)$

步骤4.2，计算包含光伏发电系统后发电系统的电力不足期望值可靠性指标按如下步骤计算：

光伏发电系统及所述的常规发电系统构成的总发电系统的三元区间数通用生成函数如下：

${\tilde{U}}_E\left(z\right)={\tilde{U}}_G\left(z\right)\&CirclePlus;{\tilde{U}}_{PV}\left(z\right)=\underset{q=1}{\overset{M}{\&Sigma;}}\&lsqb;{\tilde{p}}_q^E\&rsqb;z^{\&lsqb;{\tilde{E}}_q\&rsqb;}---\left(16\right)$

在公式(16)中，M为所述的总发电系统的状态数，为所述的总发电系统的第q个三元区间数状态，为第q个三元区间数状态对应的三元区间数概率；令加入额定容量为C_PV的光伏发电系统后其可额外承担的负荷增加量为待求变量则包含光伏发电系统后发电系统的电力不足期望值可靠性指标表示如下：

$\&lsqb;{\tilde{R}}_{PV}\&rsqb;=8760\&times;\underset{j=1}{\overset{M_L}{\&Sigma;}}\underset{q=1}{\overset{M}{\&Sigma;}}p(\&lsqb;{\tilde{E}}_q\&rsqb;<\&lsqb;{\tilde{L}}_{ej}\&rsqb;+\&lsqb;\mathrm{\&Delta;}\tilde{L}\&rsqb;)---\left(17\right)$

步骤5，运用基于三元区间数的布谷鸟搜索算法求解式(17)中的待求变量使其满足下式:

$\&lsqb;{\tilde{R}}_0\&rsqb;=\&lsqb;{\tilde{R}}_{PV}\&rsqb;---\left(18\right)$

公式(18)得到的负荷增加量即为所述光伏发电系统置信容量的三元区间数评估结果，具体按以下步骤实现：

步骤5.1，在区间[0,C_PV]范围内随机生成w个三元区间数0≤S_i ^-≤S_i ⁰≤S_i ⁺≤C_PV，i＝1，2，...，w。

步骤5.2，令代入公式(17)中，计算在该负荷增加量下的电力系统LOLE指标计算与的相似度重复该过程计算w个三元区间数下的电力系统LOLE指标及其与的相似度，并记录相似度最大值，判断其是否大于预先给定的阈值，如果大于则停止运行，且此时所述光伏发电系统的置信容量即为相似度最大时对应的负荷增加量执行步骤5.4；否则，执行步骤5.3。

所述的相似度按如下定义进行：

设两个三元区间数其相似度定义为：

$F_{\&lsqb;\tilde{A}\&rsqb;,\&lsqb;\hat{B}\&rsqb;}=\left\{\begin{array}{ccc}\frac{S_{\&lsqb;\tilde{A}\&rsqb;\&cap;\&lsqb;\tilde{B}\&rsqb;}}{\max (S_{\&lsqb;\tilde{A}\&rsqb;},S_{\&lsqb;\tilde{B}\&rsqb;})},& (S_{\&lsqb;\tilde{A}\&rsqb;\&cap;\&lsqb;\tilde{B}\&rsqb;}\&GreaterEqual;0)& \left(a\right)\\ -\frac{|b^{-}-a^{+}|}{|b^{+}-a^{-}|},& (b^{-}>a^{+})& \left(b\right)\\ -\frac{|a^{-}-b^{+}|}{|a^{+}-b^{-}|},& (a^{-}>b^{+})& \left(c\right)\end{array}\right.---\left(19\right)$

在公式(19)中，表示由点(a^-,0),(a⁰,1)和(a⁺,0)构成的三角形面积，表示由点(b^-,0),(b⁰,1)和(b⁺,0)构成的三角形面积，表示和的重叠部分面积。

步骤5.3，参考图2，基于三元区间数的布谷鸟搜索算法可按如下过程执行：对步骤5.1中所述的w个三元区间数均进行列维飞行生成新的三元区间数并令代入公式(17)中，计算在该负荷增加量下的电力系统的LOLE指标计算与的相似度比较F_i和F_i′的大小，如果F_i′>F_i，则用三元区间数替换否则，保留原三元区间数不变；针对w个新的三元区间数，记录相似度最大值，判断其是否大于预先给定的阈值。如果大于该阈值则停止运行，且此时所述光伏发电系统的置信容量即为相似度最大时对应的负荷增加量执行步骤5.4；否则，继续执行步骤5.3。

在本发明中，在保证求解结果准确性的前提下并加快三元区间数的布谷鸟搜索算法收敛速度，预先给定的阈值取值范围为(0.95，1.00)，本实施例中，预先给定的阈值为0.95。

所述的列维飞行实现过程如下：

设三元区间数初始位置为通过列维飞行后新的三元区间数更新为具体按如下进行：

$\&lsqb;{\tilde{S}}_i^{\&prime;}\&rsqb;=\&lsqb;{\tilde{S}}_i\&rsqb;+\&lsqb;\widetilde{\mathrm{\&alpha;}}\&rsqb;\&CircleTimes;\mathrm{\&gamma;}\left(\mathrm{\&epsiv;}\right)---\left(20\right)$

γ(ε)～|ε|^-λ,(1＜λ≤3) (21)

在公式(20)和(21)中，表示随机三元区间数步长，符号表示点乘计算，γ(ε)为列维随机飞行步长，ε为列维分布参数，λ为指数参数。

步骤5.4，记录步骤5.2或步骤5.3中得到的相似度最大时对应的负荷增加量即此时所述光伏发电系统的置信容量为

Claims (4)

1.一种光伏发电系统置信容量的三元区间数评估方法，其特征在于，应用于包含光伏发电系统的电力系统，所述光伏发电系统的额定容量为C_PV，C_PV>0；所述三元区间数a^-≤a⁰≤a⁺，其中a^-、a⁺和a⁰分别为三元区间数的下边界、上边界和特元；所述方法包括以下步骤：

步骤1，通过采集历史数据获取包含光伏发电系统的电力系统中常规发电机组、光伏发电系统及负荷数据，包括：各常规发电机组的额定容量G_i、强迫停运率p_i，i＝1,2,…,N，N为常规发电机组数目；负荷一年8760小时的数据并按负荷值从大到小进行降序排列以得到负荷年小时需求序列L_s＝{L₁,L₂,…,L₈₇₆₀}；光伏发电系统一年8760小时的出力数据并按光伏出力从大到小进行降序排列以得到光伏发电系统年小时出力序列C_s＝{C₁,C₂,…,C₈₇₆₀}；

步骤2，将步骤1中获取的常规发电机组额定容量G_i、强迫停运率p_i、负荷年小时需求序列L_s及光伏发电系统年小时出力序列C_s分别从点值转化为三元区间数状态及对应的三元区间数概率的形式，具体按如下步骤进行：

步骤2.1，常规发电机组额定容量G_i及强迫停运率p_i从点值到三元区间数状态及对应的三元区间数概率的转化；

设第i台常规发电机组的额定容量G_i的三元区间数状态表示为 强迫停运率p_i的点值对应的三元区间数概率通过重复n次贝努利试验得到，随机生成“0”和“1”，若生成的数据为“0”，表示强迫停运率p_i为非精确值，若生成的数据为“1”，表示强迫停运率p_i为精确值，此时强迫停运率p_i的点值对应的三元区间数概率为：

$\&lsqb;{\tilde{p}}_i\&rsqb;=\&lsqb;p_i-\frac{S}{\sqrt{n}}{t}_{\mathrm{\&alpha;}/2}(n-1),p_i,p_i+\frac{S}{\sqrt{n}}{t}_{\mathrm{\&alpha;}/2}(n-1)\&rsqb;---\left(1\right)$

在公式(1)中，S为n次贝努利试验的标准差，由下式计算得到：

$S=\sqrt{\frac{m}{n}(1-\frac{m}{n})}---\left(2\right)$

在公式(2)中，m为n次贝努利试验中发生“1”的次数，α＝0.05，t_α/2(n-1)的数值由t分布统计表查询得到；

重复上述过程得到N台常规发电机组的三元区间数状态向量表示为：

$\{\&lsqb;{\tilde{G}}_c\&rsqb;\}=\{\&lsqb;{\tilde{G}}_1\&rsqb;,\&lsqb;{\tilde{G}}_2\&rsqb;,...,\&lsqb;{\tilde{G}}_i\&rsqb;,...,\&lsqb;{\tilde{G}}_N\&rsqb;\}---\left(3\right)$

对应的三元区间数概率向量表示为：

$\{\&lsqb;{\tilde{p}}_c\&rsqb;\}=\{\&lsqb;{\tilde{p}}_1\&rsqb;,\&lsqb;{\tilde{p}}_2\&rsqb;,...,\&lsqb;{\tilde{p}}_i\&rsqb;,...,\&lsqb;{\tilde{p}}_N\&rsqb;\}---\left(4\right)$

步骤2.2，负荷年小时需求序列L_s从点值到三元区间数状态及对应的三元区间数概率的转化；

利用等宽划分的方法将所述负荷年小时需求序列L_s等分为M_L个等宽区间，每一等宽区间均采用模糊C均值聚类方法计算该区间的聚类中心，则等宽区间的左端点、右端点和聚类中心分别对应该区间三元区间数的下边界、上边界和特元，得到负荷年小时需求序列L_s的三元区间数状态向量为：

$\{\&lsqb;\tilde{L}\&rsqb;\}=\{\&lsqb;{\tilde{L}}_{e1}\&rsqb;,\&lsqb;{\tilde{L}}_{e2}\&rsqb;,...,\&lsqb;{\tilde{L}}_{ej}\&rsqb;,...,\&lsqb;{\tilde{L}}_{{\mathrm{eM}}_L}\&rsqb;\}---\left(5\right)$

对应的负荷年小时需求序列L_s的三元区间数概率向量为：

$\{\&lsqb;{\tilde{p}}_L\&rsqb;\}=\{\&lsqb;{\tilde{p}}_{e1}\&rsqb;,\&lsqb;{\tilde{p}}_{e2}\&rsqb;,...,\&lsqb;{\tilde{p}}_{ej}\&rsqb;,...,\&lsqb;{\tilde{p}}_{{\mathrm{eM}}_L}\&rsqb;\},j=1,2,...,M_L---\left(6\right)$

在公式(5)、(6)中，表示第j个等宽区间的负荷三元区间数状态，表示对应的负荷三元区间数概率；

步骤2.3，光伏发电系统年小时出力序列C_s从点值到三元区间数状态及对应的三元区间数概率的转化；

利用等宽划分的方法将所述光伏发电系统年小时出力序列C_s等分为M_c个等宽区间，每一等宽区间均采用模糊C均值聚类方法计算该区间的聚类中心，则等宽区间的左端点、右端点和聚类中心分别对应该区间三元区间数的下边界、上边界和特元，得到光伏发电系统年小时出力序列C_s的三元区间数状态向量表示为：

$\{\&lsqb;\tilde{V}\&rsqb;\}=\{\&lsqb;{\tilde{V}}_1\&rsqb;,\&lsqb;{\tilde{V}}_2\&rsqb;,...,\&lsqb;{\tilde{V}}_k\&rsqb;,...,\&lsqb;{\tilde{V}}_{M_c}\&rsqb;\}---\left(7\right)$

对应的光伏发电系统年小时出力序列的三元区间数概率向量表示为：

$\{\&lsqb;{\tilde{p}}_V\&rsqb;\}=\{\&lsqb;{\tilde{p}}_{V_1}\&rsqb;,\&lsqb;{\tilde{p}}_{V_2}\&rsqb;,...,\&lsqb;{\tilde{p}}_{V_k}\&rsqb;,...,\&lsqb;{\tilde{p}}_{V_{M_c}}\&rsqb;\},k=1,2,...,M_c;---\left(8\right)$

在公式(7)、(8)中，表示第k个等宽区间的光伏发电系统的三元区间数状态，表示对应的光伏发电系统三元区间数概率；

步骤3，由步骤2得到的常规发电机组、负荷及光伏发电系统的三元区间数状态向量及对应的三元区间数概率向量建立其对应的三元区间数通用生成函数，具体步骤如下：

步骤3.1，由步骤2.1中式(3)和(4)得到第i台常规发电机组的三元区间数通用生成函数其表达式为：

${\tilde{U}}_i\left(z\right)=\&lsqb;{\tilde{p}}_i\&rsqb;\&times;z^{\&lsqb;\tilde{0}\&rsqb;}+(\&lsqb;\tilde{1}\&rsqb;-\&lsqb;{\tilde{p}}_i\&rsqb;)\&times;z^{\&lsqb;{\tilde{G}}_i\&rsqb;},i=1,2,...,N---\left(9\right)$

则由N台常规发电机组构成的常规发电系统的三元区间数通用生成函数表示为：

$\begin{array}{c}{\tilde{U}}_G\left(z\right)=\&CirclePlus;\left({\tilde{U}}_1\right(z),{\tilde{U}}_2(z),...,{\tilde{U}}_i(z),...,{\tilde{U}}_N(z\left)\right)\\ =\underset{m=1}{\overset{M_G}{\&Sigma;}}\&lsqb;{\tilde{p}}_m^G\&rsqb;\&times;z^{\&lsqb;{\tilde{G}}_m^G\&rsqb;}\end{array}---\left(10\right)$

在公式(9)、(10)中，算子表示两个及多个三元区间数通用生成函数的“和”运算，即将各个三元区间数状态相加，对应的三元区间数概率相乘；z为通用生成函数变换因子，M_G表示由N台常规发电机组构成的常规发电系统的总状态数，m表示所述常规发电系统的第m个状态，表示所述常规发电系统的第m个三元区间数状态，表示对应的三元区间数概率；

步骤3.2，根据步骤2.2中公式(5)和(6)得到负荷的三元区间数通用生成函数其表达式为：

${\tilde{U}}_L\left(z\right)=\underset{j=1}{\overset{M_L}{\&Sigma;}}\&lsqb;{\tilde{p}}_{ej}\&rsqb;\&times;z^{\&lsqb;{\tilde{L}}_{ej}\&rsqb;}---\left(11\right)$

步骤3.3，根据步骤2.3中公式(7)和(8)得到光伏发电系统的三元区间数通用生成函数其表达式为：

${\tilde{U}}_{PV}\left(z\right)=\underset{k=1}{\overset{M_c}{\&Sigma;}}\&lsqb;{\tilde{p}}_{V_k}\&rsqb;\&times;z^{\&lsqb;{\tilde{V}}_k\&rsqb;}---\left(12\right)$

步骤4，分别计算不含光伏发电系统及包含光伏发电系统的电力系统电力不足期望值可靠性指标即LOLE指标，具体计算步骤如下：

步骤4.1，按下式计算不含光伏发电系统的电力系统电力不足期望值可靠性指标

$\&lsqb;{\tilde{R}}_0\&rsqb;=8760\&times;\underset{j=1}{\overset{M_L}{\&Sigma;}}\underset{m=1}{\overset{M_G}{\&Sigma;}}p(\&lsqb;{\tilde{G}}_m^G\&rsqb;<\&lsqb;{\tilde{L}}_{ej}\&rsqb;)---\left(13\right)$

在公式(13)中，表示三元区间数小于三元区间数的概率，其按如下过程计算：

令由点(c^-,0),(c⁰,1)和(c⁺,0)构成的三角形面积为该三角形在平面直角坐标系中y轴左侧的面积记为S_left，则表示为：

$p(\&lsqb;{\tilde{G}}_m^G\&rsqb;<\&lsqb;{\tilde{L}}_{ej}\&rsqb;)=\frac{S_{left}}{S_{\&lsqb;\tilde{c}\&rsqb;}}---\left(14\right)$

步骤4.2，计算包含光伏发电系统后发电系统的电力不足期望值可靠性指标按如下步骤计算：

光伏发电系统及所述的常规发电系统构成的总发电系统的三元区间数通用生成函数如下：

${\tilde{U}}_E\left(z\right)={\tilde{U}}_G\left(z\right)\&CirclePlus;{\tilde{U}}_{PV}\left(z\right)=\underset{q=1}{\overset{M}{\&Sigma;}}\&lsqb;{\tilde{p}}_q^E\&rsqb;z^{\&lsqb;{\tilde{E}}_q\&rsqb;}---\left(15\right)$

在公式(15)中，M为所述的总发电系统的状态数，为所述的总发电系统的第q个三元区间数状态，为第q个三元区间数状态对应的三元区间数概率；令加入额定容量为C_PV的光伏发电系统后其可额外承担的负荷增加量为待求变量则包含光伏发电系统后发电系统的电力不足期望值可靠性指标表示如下：

$\&lsqb;{\tilde{R}}_{PV}\&rsqb;=8760\&times;\underset{j=1}{\overset{M_L}{\&Sigma;}}\underset{q=1}{\overset{M}{\&Sigma;}}p(\&lsqb;{\tilde{E}}_q\&rsqb;<\&lsqb;{\tilde{L}}_{ej}\&rsqb;+\&lsqb;\mathrm{\&Delta;}\tilde{L}\&rsqb;)---\left(16\right)$

步骤5，运用基于三元区间数的布谷鸟搜索算法求解式(16)中的待求变量使其满足下式:

$\&lsqb;{\tilde{R}}_0\&rsqb;=\&lsqb;{\tilde{R}}_{PV}\&rsqb;---\left(17\right)$

公式(17)得到的负荷增加量即为所述光伏发电系统置信容量的三元区间数评估结果，具体按以下步骤实现：

步骤5.1，在区间[0,C_PV]范围内随机生成w个三元区间数

步骤5.2，令代入公式(16)中，计算在该负荷增加量下的电力系统LOLE指标计算与的相似度重复该过程计算w个三元区间数下的电力系统LOLE指标及其与的相似度，并记录相似度最大值，判断其是否大于预先给定的阈值，如果大于则停止运行，且此时所述光伏发电系统的置信容量即为相似度最大时对应的负荷增加量执行步骤5.4；否则，执行步骤5.3；

步骤5.3，对步骤5.1中所述的w个三元区间数均进行列维飞行生成新的三元区间数并令代入公式(16)中，计算在该负荷增加量下的电力系统的LOLE指标计算与的相似度比较F_i和F′_i的大小，如果F′i>F_i，则用三元区间数替换否则，保留原三元区间数不变；针对w个新的三元区间数，记录相似度最大值，判断其是否大于预先给定的阈值，如果大于则停止运行，且此时所述光伏发电系统的置信容量即为相似度最大时对应的负荷增加量执行步骤5.4；否则，继续执行步骤5.3；

步骤5.4，记录步骤5.2或步骤5.3中得到的相似度最大时对应的负荷增加量即此时所述光伏发电系统的置信容量为

2.根据权利要求1所述的光伏发电系统置信容量的三元区间数评估方法，其特征在于，步骤5中所述的列维飞行实现过程如下：

设三元区间数初始位置为通过列维飞行后新的三元区间数更新为具体按如下进行：

$\&lsqb;{\tilde{S}}_i^{\&prime;}\&rsqb;=\&lsqb;{\tilde{S}}_i\&rsqb;+\&lsqb;\widetilde{\mathrm{\&alpha;}}\&rsqb;\&CircleTimes;\mathrm{\&gamma;}\left(\mathrm{\&epsiv;}\right)---\left(18\right)$

γ(ε)～|ε|^-λ,(1＜λ≤3) (19)

在公式(18)和(19)中，表示随机三元区间数步长，符号表示点乘计算，γ(ε)为列维随机飞行步长，ε为列维分布参数，λ为指数参数。

3.根据权利要求1所述的光伏发电系统置信容量的三元区间数评估方法，其特征在于，步骤5中所述的相似度按如下定义进行：

设两个三元区间数其相似度定义为：

$F_{\&lsqb;\tilde{A}\&rsqb;,\&lsqb;\hat{B}\&rsqb;}=\left\{\begin{array}{ccc}\frac{S_{\&lsqb;\tilde{A}\&rsqb;\&cap;\&lsqb;\tilde{B}\&rsqb;}}{max(S_{\&lsqb;\tilde{A}\&rsqb;},S_{\&lsqb;\tilde{B}\&rsqb;})},& (S_{\&lsqb;\tilde{A}\&rsqb;\&cap;\&lsqb;\tilde{B}\&rsqb;}\&GreaterEqual;0)& \left(a\right)\\ -\frac{|b^{-}-a^{+}|}{|b^{+}-a^{-}|},& (b^{-}>a^{+})& \left(b\right)\\ -\frac{|a^{-}-b^{+}|}{|a^{+}-b^{-}|},& (a^{-}>b^{+})& \left(c\right)\end{array}\right.---\left(20\right)$

在公式(20)中，表示由点(a^-,0),(a⁰,1)和(a⁺,0)构成的三角形面积，表示由点(b^-,0),(b⁰,1)和(b⁺,0)构成的三角形面积，表示和的重叠部分面积。

4.根据权利要求1所述的光伏发电系统置信容量的三元区间数评估方法，其特征在于，步骤5中所述的预先给定的阈值的取值范围为(0.95，1.00)。