CN105719059B - 光伏发电系统置信容量评估方法 - Google Patents
光伏发电系统置信容量评估方法 Download PDFInfo
- Publication number
- CN105719059B CN105719059B CN201610025199.2A CN201610025199A CN105719059B CN 105719059 B CN105719059 B CN 105719059B CN 201610025199 A CN201610025199 A CN 201610025199A CN 105719059 B CN105719059 B CN 105719059B
- Authority
- CN
- China
- Prior art keywords
- photovoltaic
- load
- generating system
- capacity
- year
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000004364 calculation method Methods 0.000 claims abstract description 32
- 238000005070 sampling Methods 0.000 claims abstract description 30
- 238000013528 artificial neural network Methods 0.000 claims abstract description 19
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 18
- 230000035699 permeability Effects 0.000 claims abstract description 13
- 230000005611 electricity Effects 0.000 claims description 12
- 238000010248 power generation Methods 0.000 claims description 12
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 5
- 240000002853 Nelumbo nucifera Species 0.000 claims description 5
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 5
- 230000006870 function Effects 0.000 claims description 5
- 238000003646 Spearman's rank correlation coefficient Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000012614 Monte-Carlo sampling Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 210000002569 neuron Anatomy 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims 1
- 238000013480 data collection Methods 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 9
- 238000004146 energy storage Methods 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000010977 unit operation Methods 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Development Economics (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Photovoltaic Devices (AREA)
Abstract
本发明公开一种光伏发电系统置信容量评估方法,包括提出影响光伏发电系统置信容量的四个因素:光伏渗透率、数据采样时间间隔、光伏‑负荷出力波动性及光伏出力与负荷波动的时序相关性;利用序贯蒙特卡洛仿真计算电力系统可靠性,通过弦截法求解置信容量;建立四个因素与置信容量之间的三层误差反向传播神经网络;利用不同辐照模式下得到的输入输出数据集训练神经网络;利用训练完成的神经网络的泛化能力获得给定条件下的光伏发电系统置信容量;本发明提出的光伏发电系统置信容量评估方法可用于含光伏的发输电系统的规划设计,且无需再单独对每个特定的光伏发电系统进行序贯蒙特卡洛仿真,解决了现有技术中通用性不强的问题。
Description
技术领域
本发明涉及一种光伏发电系统置信容量评估方法,属于光伏发电领域。
背景技术
随着光伏发电成本的降低、技术的进步及其对环境的友好性,大型地面光伏发电系统得到快速发展,安装容量日益增加,光伏发电在电力系统中的渗透率逐渐增加。现有研究表明,光伏发电不仅具有能量价值,而且具有容量价值,其容量价值体现在光伏发电系统有效可信容量,即为光伏发电系统的置信容量。为了避免重复建设及资源的浪费,光伏发电系统置信容量的评估成为目前热门研究课题。
现有的关于光伏发电系统置信容量的研究主要分为两大类:基于近似解析法的光伏发电系统置信容量评估和基于蒙特卡洛仿真的光伏发电系统置信容量评估。目前,光伏发电系统置信容量的研究主要集中于对特定的光伏发电系统进行置信容量的评估,如王秀丽,武泽辰,曲翀三人于2014年在中国电机工程学报上发表的文献《光伏发电系统可靠性分析及其置信容量计算》中,利用分层抽样的序贯蒙特卡洛方法进行光伏发电系统置信容量的评估,且研究了天气变化、设备故障等随机因素对光伏发电系统置信容量计算的影响。基于蒙特卡洛的仿真法虽然精确度高,但是,该计算方法存在以下不足:
1)针对特定的光伏发电系统均需要重新进行新的序贯蒙特卡洛仿真,对同样的光伏发电系统而言,当外界环境如辐照强度、负荷变化模式不同时,均需重新进行计算,该算法的适应性和鲁棒性不强;
2)该算法时间消耗大,当常规发电系统规模大,光伏发电系统数目众多时,计算开销将明显增加,甚至在一些极端情况下导致算法不收敛的现象发生。
为了降低蒙特卡洛仿真带来的高时间消耗,Madaeni S.H.,Sioshansi R.,和Denholm P.三人于2013年在期刊IEEE Journal of Photovoltaics上发表的文献“Comparing capacity value estimation techniques for photovoltaic solar power,2013,3(1):407-415.”(Madaeni S.H.,Sioshansi R.,和Denholm P.(“光伏发电系统的容量价值评估方法比较分析”,IEEE光伏期刊,2013年第3卷第1期,ISNN 21563381,第407-415页)一文总结了光伏发电系统置信容量计算的各种近似解析法,其中包括基于光伏发电容量因子的近似方法、Garver近似计算方法、光伏发电系统多状态机组等效方法及Z-method方法等,这些近似方法采用解析法去逼近仿真算法得到的精确值,其在计算光伏发电系统置信容量中,只需进行简单的解析运算。但是,这种近似解析法存在如下不足:
1)近似解析法评估虽然计算量小,但其精确度低,不同的解析法其计算准确度各不相同,与蒙特卡洛仿真法得到的精确解误差各不相同,在不同的场景下选用何种近似解析法则需要具体计算后方能确定,这在实际中不利于工程应用;
2)对于同一种近似解析法,随着光伏渗透率、年辐照强度、负荷的变化不同,其评估得到的结果的误差也将变化,导致该方法在使用中将产生不确定性。
中国发明专利申请公开说明书CN104077664A于2014年10月01日公开的《一种风电储能发电系统的置信容量评估方法》介绍了与光伏发电类似的风力发电/储能系统的置信容量评估方法,该发明专利考虑了风电和储能捆绑下的风电场置信容量评估,解决了以风电置信容量最大化为目标下的储能系统的容量配置问题。该方法能计算光伏发电系统的置信容量,但是针对每一个特定的系统都得重新搭建系统,计算光伏发电系统的置信容量,通用性不强,其模型和方法的通用性和鲁棒性较差。同时,未深层次分析影响置信容量的因素及其对置信容量的影响程度。
综上所述,为了更好的评估光伏发电系统的置信容量,不仅需要考虑影响光伏发电置信容量的主要因素,如光伏渗透率、数据采样时间间隔、光伏-负荷出力波动性、光伏出力与负荷波动的时序相关性,而且还需要提出一种具有通用性的光伏发电系统置信容量评估方法。
发明内容
本发明要解决的技术问题为克服现有技术中存在的问题,提出一种综合考虑各种因素,并利用训练完成的BP神经网络实现计算任意给定条件下的光伏发电系统置信容量评估方法。
本发明的目的是这样实现的,本发明提供了一种光伏发电系统置信容量评估方法,其主要步骤如下:
步骤1,通过采样获取发电系统常规机组额定容量及可靠性数据,包括每台常规机组额定容量Gi、平均正常工作时间tMTTF_i、平均故障时间tMTTR_i(i=1,2,…,N),其中N为常规机组台数,所有常规机组额定容量和为光伏发电系统额定容量CPV、1分钟级的年辐照强度数据、1分钟级的年负荷数据;
步骤2,在仿真时段T中,以Δt为数据采样时间间隔,Δt≥1min;通过等宽采样法得到数据采样时间间隔Δt下的光伏年出力序列和负荷年需求序列VLoad={LΔt,L2Δt,...,LiΔt,...,LnΔt}(i=1,2,…,n),其中n为仿真时段T中的时间间隔数目;采用序贯蒙特卡洛抽样法得到常规机组的年出力状态序列G;先将所述年出力状态序列G以1分钟为时间间隔采样得到1分钟级的出力序列,再以数据采样时间间隔Δt为宽度通过等宽采样法,获取以Δt为间隔的常规机组运行及停运两状态下的年出力序列VGen;
步骤3,计算光伏渗透率r,数据采样时间间隔Δt下的光伏‐负荷出力波动性υ,光伏出力与负荷波动的时序相关性Cr,将其记为输入向量(r,Δt,υ,Cr),具体按下式进行:
(1)光伏渗透率r:
(2)光伏-负荷出力波动性υ:
在(2)式中,n为仿真时段T中的时间间隔数目,和分别表示光伏发电系统在(i+1)Δt和iΔt时刻的出力,L(i+1)Δt和LiΔt分别表示(i+1)Δt和iΔt时刻的负荷需求,且
(3)光伏出力与负荷波动的时序相关性Cr:
在(3)式中,ρs为斯皮尔曼秩相关系数,按下式计算得到:
在(4)式中,di=xi-yi,xi和yi分别表示变量和LiΔt在相应的光伏年出力序列VPV和负荷年需求序列VLoad中的秩次;
δa为弗雷歇平均距离,按下式计算得到:
在(5)式中,Li′=LiΔt/max(VPV,VLoad),max(VPV,VLoad)为序列VPV和VLoad的最大值;
步骤4,利用序贯蒙特卡洛仿真计算不含光伏发电系统时的发电系统电力不足期望值可靠性指标R0,计算公式如下:
假定新增光伏发电系统额定容量CPV可额外承担的有效载荷容量为ΔL,此时加入安装容量为CPV的光伏发电系统后的发电系统可靠性指标发电系统可靠性RΔL计算公式如下:
在(6)、(7)式中,Gt表示常规机组年出力序列VGen在t时刻的出力,Lt为负荷年需求序列VLoad在t时刻的负荷需求,表示光伏年出力序列VPV在t时刻的出力,表示在仿真时段T内,t时刻的常规机组出力Gt小于负荷需求Lt的概率之和;
步骤5,通过弦截法求解R0=RΔL时的新增光伏发电系统额定容量CPV可额外承担的有效载荷容量ΔL,并利用下式计算光伏发电系统的置信容量kc和光伏发电系统的容量置信度kp,并将其记录为输出向量(kc,kp):
kc=ΔL (8)
步骤6,在保持步骤1中常规机组容量及可靠性数据不变的前提下改变光伏年出力序列及负荷年需求序列,重复步骤2至步骤5,计算得到在多个不同数据采样时间间隔Δt下的光伏渗透率r、光伏‐负荷出力波动性υ、光伏出力与负荷波动的时序相关性Cr、光伏发电系统的置信容量kc和光伏发电系统的容量置信度kp的数值,记录多个以(r,Δt,υ,Cr)为输入向量,以(kc,kp)为输出向量的数据集Sd;建立含输入层、隐含层、输出层的三层误差反向传播BP神经网络,输入层以(r,Δt,υ,Cr)为四个输入,隐含层有4个神经元,输出层以(kc,kp)为两个输出,利用得到的数据集Sd完成BP神经网络的训练;
步骤7,利用训练完成的BP神经网络的泛化能力得到给定条件(r,Δt,υ,Cr)下的光伏发电系统的置信容量kc和容量置信度kp。
优选的,步骤5中所述的弦截法按如下步骤进行:
通过序贯蒙特卡洛仿真计算得到在负荷峰值为Lpk0、电源侧加入新增光伏发电系统额定容量CPV时的发电系统可靠性指标RA,并记为点A(Lpk0,RA);通过序贯蒙特卡洛仿真计算得到在负荷侧加入与新增光伏发电系统额定容量CPV等容量的负荷后的发电系统可靠性指标RB,并记为点B(Lpk0+CPV,RB);
连接A、B点记为线段lAB,求解lAB与y=R0的直线交点为P1(L1,R0),将L1作为新的负荷加入原发电系统并通过序贯蒙特卡洛仿真计算此时的发电系统可靠性指标RC,记为C(L1,RC)点;
连接B、C点记为线段lBC,求解lBC与y=R0的直线交点为P2(L2,R0),将L2作为新的负荷加入原发电系统通过序贯蒙特卡洛仿真计算此时的发电系统可靠性指标RD,记为D(L2,RC)点;
继续迭代上述过程,直到相邻两次计算得到的发电系统可靠性指标小于收敛阈值ε=0.01,停止迭代,记录此时对应的峰值负荷LE,则新增光伏发电系统额定容量CPV可额外承担的有效载荷容量ΔL=LE-Lpk0。优选的,步骤6中所述的BP神经网络的训练采用列文伯格‐马夸尔特优化算法,训练目标最小误差取0.05,训练次数取600次,学习率取0.05,隐含层和输出层的传递函数均采用S型正切函数。
与已有技术相比,本发明的有益效果体现在:
1、本发明的评估方法分析了影响光伏发电系统置信容量的四个因素,便于研究各因素对光伏发电系统置信容量评估的影响,与已有技术中只是单纯计算光伏发电系统置信容量的结果相比,本发明提出的方法能研究影响光伏发电系统置信容量评估的主要因素,从而为挖掘提高光伏发电系统置信容量的有效措施提供参考价值。
2、本发明的评估方法提出的光伏-负荷出力波动性指标υ反映了不同数据采样时间间隔Δt下的光伏和负荷的出力波动情况,提出的光伏出力与负荷波动的时序相关性Cr指标结合斯皮尔曼秩相关系数和弗雷歇平均距离,不仅表征了光伏出力与负荷波动趋势的关系,而且反映了光伏出力序列和负荷需求序列的平均距离,能综合反映光伏出力和负荷曲线之间的匹配关系。
3、本发明中采用三层BP神经网络训练多种辐照模式下的光伏发电系统置信容量,利用训练完成的神经网络的泛化能力计算其他任意给定条件下的光伏发电系统的置信容量,较其他对特定系统必须重新计算的方法具有较强的通用性和较短的计算时间,为光伏发电系统置信容量的评估提供了更便捷的方法。
附图说明
图1为本发明光伏发电系统置信容量评估方法的流程图;
图2为本发明弦截法迭代求解过程示意图;
图3为本发明光伏发电系统置信容量评估BP神经网络结构图。
具体实施方式
下面结合附图对本发明的优选方式作进一步详细的描述。
如图1所示,本发明的主要步骤如下:
步骤1,通过采样获取发电系统常规机组额定容量及可靠性数据,包括每台常规机组额定容量Gi、平均正常工作时间tMTTF_i、平均故障时间tMTTR_i,(i=1,2,…,N),其中N为常规机组台数,所有常规机组额定容量和为光伏发电系统额定容量CPV、1分钟级的年辐照强度数据、1分钟级的年负荷数据。
实际中获取的负荷数据通常以小时为间隔,为了得到1分钟级的负荷数据可以通过两数据点之间的线性插值得到;
步骤2,在仿真时段T中,以Δt为数据采样时间间隔,Δt≥1min;通过等宽采样法得到数据采样时间间隔Δt下的光伏年出力序列和负荷年需求序列VLoad={LΔt,L2Δt,...,LiΔt,...,LnΔt}(i=1,2,…,n),其中n为仿真时段T中的时间间隔数目;采用序贯蒙特卡洛抽样法得到常规机组的年出力状态序列G;先将所述年出力状态序列G以1分钟为时间间隔采样得到1分钟级的出力序列,再以数据采样时间间隔Δt为宽度通过等宽采样法,获取以Δt为间隔的常规机组运行及停运两状态下的年出力序列VGen。
步骤3,计算光伏渗透率r,数据采样时间间隔Δt下的光伏‐负荷出力波动性υ,光伏出力与负荷波动的时序相关性Cr,将其记为输入向量(r,Δt,υ,Cr),具体按下式进行:
(1)光伏渗透率r:
(2)光伏-负荷出力波动性υ:
在(2)式中,n为仿真时段T中的时间间隔数目,和分别表示光伏发电系统在(i+1)Δt和iΔt时刻的出力,L(i+1)Δt和LiΔt分别表示(i+1)Δt和iΔt时刻的负荷需求,且
(3)光伏出力与负荷波动的时序相关性Cr:
在(3)式中,ρs为斯皮尔曼秩相关系数,按下式计算得到:
在(4)式中,di=xi-yi,xi和yi分别表示变量和LiΔt在相应的光伏年出力序列VPV和负荷年需求序列VLoad中的秩次;
δa为弗雷歇平均距离,按下式计算得到:
在(5)式中,Li′=LiΔt/max(VPV,VLoad),max(VPV,VLoad)为序列VPV和VLoad中的最大值。
光伏出力与负荷波动的时序相关性Cr指标不仅表征了光伏年出力序列和负荷年需求序列的波动趋势,同时也反应了两者之间的平均距离大小,当光伏年出力序列和负荷年需求序列趋势一致其平均距离很小时,Cr指标取值较大;反之,当Cr指标取值较小时,说明光伏年出力序列和负荷年需求序列平均距离大且变化趋势不一致。
步骤4,利用序贯蒙特卡洛仿真计算不含光伏发电系统时的发电系统电力不足期望值可靠性指标R0,计算公式如下:
假定新增光伏发电系统额定容量CPV可额外承担的有效载荷容量为ΔL,此时加入安装容量为CPV的光伏发电系统后的发电系统可靠性指标发电系统可靠性RΔL计算公式如下:
在(6)、(7)式中,Gt表示常规机组年出力序列VGen在t时刻的出力,Lt为负荷年需求序列VLoad在t时刻的负荷需求,表示光伏年出力序列VPV在t时刻的出力,表示在仿真时段T内,t时刻的常规机组出力Gt小于负荷需求Lt的概率之和。
因负荷年需求序列及光伏年出力序列均可从实际中获取,故本发明只针对常规机组进行序贯蒙特卡洛采样得到常规机组年出力序列。
所述的序贯蒙特卡洛仿真实现过程如下:
由程序生成位于[0,1]之间的随机数γ1和γ2,此时第i台常规机组的运行时间和故障停运时间分别为:
其中tMTTF_i和tMTTR_i分别为第i台常规机组的平均正常工作时间和平均故障时间,此时每次抽样得到一对运行和故障停运状态;
继续上述过程,直到总的状态累计时间和为8760小时,停止抽样,得到第i台常规机组一年的运行和停运两状态年序列;继续对N台常规机组均按照该步骤进行序贯蒙特卡洛仿真,将所有机组一年的运行和停运两状态年序列累加即得到所有常规机组的运行及停运两状态下的年出力序列VGen。
步骤5,通过弦截法求解R0=RΔL时的新增光伏CPV可额外承担的有效载荷容量ΔL,即在保证原先发电系统可靠性指标不变的前提下,新增光伏发电系统能够额外承担的有效载荷容量ΔL大小即为该光伏发电系统的置信容量;并利用下式计算光伏发电系统的置信容量kc和光伏发电系统的容量置信度kp,并将其记录为输出向量(kc,kp):
kc=ΔL (8)
所述的弦截法具体按如下步骤进行:
通过序贯蒙特卡洛仿真计算得到在负荷年需求序列VLoad={LΔt,L2Δt,...,LiΔt,...,LnΔt}(i=1,2,…,n)的负荷峰值为Lpk0、电源侧加入新增光伏发电系统额定容量CPV时的发电系统可靠性指标RA,并记为点A(Lpk0,RA),其计算公式如下:
通过序贯蒙特卡洛仿真计算得到在负荷侧加入与新增光伏发电系统额定容量CPV等容量的负荷后的发电系统可靠性指标RB,并记为点B(Lpk0+CPV,RB),其计算公式如下:
A点和B点的位置如图2所示。
连接A、B点记为线段lAB,求解lAB与y=R0的直线交点为P1(L1,R0),将L1作为新的负荷加入原发电系统通过序贯蒙特卡洛仿真计算此时的发电系统可靠性指标RC,图2中记为C(L1,RC)点。
连接B、C点记为线段lBC,求解lBC与y=R0的直线交点为P2(L2,R0),将L2作为新的负荷加入原发电系统通过序贯蒙特卡洛仿真计算此时的发电系统可靠性指标RD,图2中记为D(L2,RC)点。
继续迭代上述过程,直到相邻两次计算得到的发电系统可靠性指标小于收敛阈值ε=0.01,停止迭代,记录此时对应的峰值负荷LE,则新增光伏CPV可额外承担的有效载荷容量ΔL=LE-Lpk0。
步骤6,在保持步骤1中常规机组容量及可靠性数据不变的前提下改变光伏年出力序列及负荷年需求序列,重复步骤2至步骤5,计算得到在多个不同数据采样时间间隔Δt下的光伏渗透率r、光伏‐负荷出力波动性υ、光伏出力与负荷波动的时序相关性Cr、光伏发电系统的置信容量kc和光伏发电系统的容量置信度kp的数值,记录多个以(r,Δt,υ,Cr)为输入向量,以(kc,kp)为输出向量的数据集Sd;建立含输入层、隐含层、输出层的三层误差反向传播BP神经网络,输入层以(r,Δt,υ,Cr)为四个输入,隐含层有4个神经元,输出层以(kc,kp)为两个输出,利用得到的数据集Sd完成BP神经网络的训练;
所述的BP神经网络的训练采用列文伯格‐马夸尔特优化算法,训练目标最小误差取0.05,训练次数取600次,学习率取0.05,隐含层和输出层的传递函数均采用S型正切函数。BP神经网络结构图如图3所示。
步骤7,利用训练完成的BP神经网络的泛化能力得到给定条件(r,Δt,υ,Cr)下的光伏发电系统的置信容量kc和容量置信度kp。
Claims (3)
1.一种光伏发电系统置信容量评估方法,其特征在于,本评估方法的主要步骤如下:
步骤1,通过采样获取发电系统常规机组额定容量及可靠性数据,包括每台常规机组额定容量Gi、平均正常工作时间tMTTF_i、平均故障时间tMTTR_i,(i=1,2,…,N),其中N为常规机组台数,所有常规机组额定容量和为光伏发电系统额定容量CPV、1分钟级的年辐照强度数据、1分钟级的年负荷数据;
步骤2,在仿真时段T中,以Δt为数据采样时间间隔,Δt≥1min;通过等宽采样法得到数据采样时间间隔Δt下的光伏年出力序列和负荷年需求序列VLoad={LΔt,L2Δt,...,LiΔt,...,LnΔt}(i=1,2,…,n),其中n为仿真时段T中的时间间隔数目;采用序贯蒙特卡洛抽样法得到常规机组的年出力状态序列G;先将所述年出力状态序列G以1分钟为时间间隔采样得到1分钟级的出力序列,再以数据采样时间间隔Δt为宽度通过等宽采样法,获取以Δt为间隔的常规机组运行及停运两状态下的年出力序列VGen;
步骤3,计算光伏渗透率r,数据采样时间间隔Δt下的光伏-负荷出力波动性υ,光伏出力与负荷波动的时序相关性Cr,将其记为输入向量(r,Δt,υ,Cr),具体按下式进行:
(1)光伏渗透率r:
(2)光伏-负荷出力波动性υ:
在(2)式中,n为仿真时段T中的时间间隔数目,和分别表示光伏发电系统在(i+1)Δt和iΔt时刻的出力,L(i+1)Δt和LiΔt分别表示(i+1)Δt和iΔt时刻的负荷需求,且
(3)光伏出力与负荷波动的时序相关性Cr:
在(3)式中,ρs为斯皮尔曼秩相关系数,按下式计算得到:
在(4)式中,di=xi-yi,xi和yi分别表示变量和LiΔt在相应的光伏年出力序列VPV和负荷年需求序列VLoad中的秩次;
δa为弗雷歇平均距离,按下式计算得到:
在(5)式中,Li′=LiΔt/max(VPV,VLoad),max(VPV,VLoad)为序列VPV和VLoad中的最大值;
步骤4,利用序贯蒙特卡洛仿真计算不含光伏发电系统时的发电系统电力不足期望值可靠性指标R0,计算公式如下:
假定新增光伏发电系统额定容量CPV可额外承担的有效载荷容量为ΔL,此时加入安装容量为CPV的光伏发电系统后的发电系统可靠性指标发电系统可靠性RΔL计算公式如下:
在(6)、(7)式中,Gt表示常规机组年出力序列VGen在t时刻的出力,Lt为负荷年需求序列VLoad在t时刻的负荷需求,表示光伏年出力序列VPV在t时刻的出力,表示在仿真时段T内,t时刻的常规机组出力Gt小于负荷需求Lt的概率之和;
步骤5,通过弦截法求解R0=RΔL时的新增光伏发电系统额定容量CPV可额外承担的有效载荷容量ΔL,并利用下式计算光伏发电系统的置信容量kc和光伏发电系统的容量置信度kp,并将其记录为输出向量(kc,kp):
kc=ΔL (8)
步骤6,在保持步骤1中常规机组容量及可靠性数据不变的前提下改变光伏年出力序列及负荷年需求序列,重复步骤2至步骤5,计算得到在多个不同数据采样时间间隔Δt下的光伏渗透率r、光伏-负荷出力波动性υ、光伏出力与负荷波动的时序相关性Cr、光伏发电系统的置信容量kc和光伏发电系统的容量置信度kp的数值,记录多个以(r,Δt,υ,Cr)为输入向量,以(kc,kp)为输出向量的数据集Sd;建立含输入层、隐含层、输出层的三层误差反向传播BP神经网络,输入层以(r,Δt,υ,Cr)为四个输入,隐含层有4个神经元,输出层以(kc,kp)为两个输出,利用得到的数据集Sd完成BP神经网络的训练;
步骤7,利用训练完成的BP神经网络的泛化能力得到给定条件(r,Δt,υ,Cr)下的光伏发电系统的置信容量kc和容量置信度kp。
2.根据权利要求1所述的一种光伏发电系统置信容量评估方法,其特征在于,步骤5中所述的弦截法按如下步骤进行:
步骤5.1,通过序贯蒙特卡洛仿真计算得到在负荷峰值为Lpk0、电源侧加入新增光伏发电系统额定容量CPV时的发电系统可靠性指标RA,并记为点A(Lpk0,RA);通过序贯蒙特卡洛仿真计算得到在负荷侧加入与新增光伏发电系统额定容量CPV等容量的负荷后的发电系统可靠性指标RB,并记为点B(Lpk0+CPV,RB);
步骤5.2,连接A、B点记为线段lAB,求解lAB与y=R0的直线交点为P1(L1,R0),将L1作为新的负荷加入原发电系统并通过序贯蒙特卡洛仿真计算此时的发电系统可靠性指标RC,记为C(L1,RC)点;
步骤5.3,连接B、C点记为线段lBC,求解lBC与y=R0的直线交点为P2(L2,R0),将L2作为新的负荷加入原发电系统通过序贯蒙特卡洛仿真计算此时的发电系统可靠性指标RD,记为D(L2,RC)点;
步骤5.4,继续迭代步骤5.2-步骤5.3中的上述过程,直到相邻两次计算得到的发电系统可靠性指标小于收敛阈值ε=0.01,停止迭代,记录此时对应的峰值负荷LE,则新增光伏发电系统额定容量CPV可额外承担的有效载荷容量ΔL=LE-Lpk0。
3.根据权利要求1所述的一种光伏发电系统置信容量评估方法,其特征在于,步骤6中所述的BP神经网络的训练采用列文伯格-马夸尔特优化算法,训练目标最小误差取0.05,训练次数取600次,学习率取0.05,隐含层和输出层的传递函数均采用S型正切函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610025199.2A CN105719059B (zh) | 2016-01-13 | 2016-01-13 | 光伏发电系统置信容量评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610025199.2A CN105719059B (zh) | 2016-01-13 | 2016-01-13 | 光伏发电系统置信容量评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105719059A CN105719059A (zh) | 2016-06-29 |
CN105719059B true CN105719059B (zh) | 2019-04-09 |
Family
ID=56147189
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610025199.2A Expired - Fee Related CN105719059B (zh) | 2016-01-13 | 2016-01-13 | 光伏发电系统置信容量评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105719059B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106777487B (zh) * | 2016-11-18 | 2019-11-05 | 清华大学 | 一种含有储能系统的光伏电站可信容量计算方法及系统 |
CN107316115A (zh) * | 2017-07-28 | 2017-11-03 | 国网山东省电力公司经济技术研究院 | 分布式电源不同渗透率下网架规划方法 |
CN112909969B (zh) * | 2019-11-19 | 2023-05-16 | 电力规划总院有限公司 | 应用于储能系统的容量置信度确定方法与装置 |
CN111461479B (zh) * | 2020-01-19 | 2023-07-14 | 湖北省电力勘测设计院有限公司 | 一种远洋海岛换电船舶的置信容量的评估方法 |
CN111293681B (zh) * | 2020-01-22 | 2023-04-25 | 重庆大学 | 一种基于RankBoost的光伏场站出力波动性量化评价方法 |
CN111666664B (zh) * | 2020-05-22 | 2023-10-27 | 湖北省电力勘测设计院有限公司 | 电-气综合能源系统新增供能元件置信容量的评估方法 |
CN113609686B (zh) * | 2021-08-09 | 2023-08-29 | 国网浙江省电力有限公司经济技术研究院 | 一种新能源置信容量分析方法及系统 |
CN113644669B (zh) * | 2021-08-17 | 2023-10-13 | 山东电力工程咨询院有限公司 | 一种基于储能容量利用率的储能容量配置方法及系统 |
CN113705103A (zh) * | 2021-08-31 | 2021-11-26 | 武汉大学 | 一种基于改进bp神经网络的分布式光伏功率分离方法 |
CN115018282A (zh) * | 2022-05-24 | 2022-09-06 | 广东电力交易中心有限责任公司 | 一种间断性可再生能源的置信容量计算方法 |
CN115065052B (zh) * | 2022-07-15 | 2024-07-19 | 华北电力大学 | 基于互联网数据中心需求的电力系统置信容量确定方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103218757A (zh) * | 2013-04-28 | 2013-07-24 | 重庆大学 | 一种确定光伏发电容量可信度的方法 |
CN104077664A (zh) * | 2014-07-09 | 2014-10-01 | 东南大学 | 一种风电储能发电系统的置信容量评估方法 |
-
2016
- 2016-01-13 CN CN201610025199.2A patent/CN105719059B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103218757A (zh) * | 2013-04-28 | 2013-07-24 | 重庆大学 | 一种确定光伏发电容量可信度的方法 |
CN104077664A (zh) * | 2014-07-09 | 2014-10-01 | 东南大学 | 一种风电储能发电系统的置信容量评估方法 |
Non-Patent Citations (2)
Title |
---|
光伏发电系统可靠性分析及其置信容量计算;王秀丽 等;《中国机电工程学报》;20140105;第34卷(第1期);1-14页 |
大规模光伏发电对电力系统影响综述;丁明 等;《中国电机工程学报》;20140105;第34卷(第1期);15-21页 |
Also Published As
Publication number | Publication date |
---|---|
CN105719059A (zh) | 2016-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105719059B (zh) | 光伏发电系统置信容量评估方法 | |
Li et al. | Long-term complementary operation of a large-scale hydro-photovoltaic hybrid power plant using explicit stochastic optimization | |
Jamil et al. | Uncertainty analysis of energy production for a 3× 50 MW AC photovoltaic project based on solar resources | |
Henning et al. | A comprehensive model for the German electricity and heat sector in a future energy system with a dominant contribution from renewable energy technologies—Part I: Methodology | |
CN108734331B (zh) | 基于lstm的短期光伏发电功率预测方法及系统 | |
Yona et al. | Determination method of insolation prediction with fuzzy and applying neural network for long-term ahead PV power output correction | |
CN103218757B (zh) | 一种确定光伏发电容量可信度的方法 | |
CN102385654B (zh) | 一种基于vsc并网光伏、储能电站的机电暂态建模方法 | |
Han et al. | Tabu search algorithm optimized ANN model for wind power prediction with NWP | |
CN103823504B (zh) | 一种基于最小二乘支持向量机的最大功率跟踪控制方法 | |
CN110543720B (zh) | 基于sdae-elm伪量测模型的状态估计方法 | |
Yuan et al. | Research on dynamic load modelling based on power quality monitoring system | |
Shen et al. | Generation scheduling of a hydrodominated provincial system considering forecast errors of wind and solar power | |
CN102479347A (zh) | 基于数据驱动的风电场短期风速预测方法和系统 | |
CN116581755B (zh) | 功率预测方法、装置、设备及存储介质 | |
Varzaneh et al. | Optimal energy management for PV‐integrated residential systems including energy storage system | |
Wang et al. | Analysis of wind farm output characteristics based on descriptive statistical analysis and envelope domain | |
Kim et al. | A probabilistic approach to potential estimation of renewable energy resources based on augmented spatial interpolation | |
Chen et al. | N‐1 security‐constrained coordinated scheduling of integrated electricity and natural gas system considering gas dynamics and wind power uncertainty | |
Liu et al. | Effect of the quality of streamflow forecasts on the operation of cascade hydropower stations using stochastic optimization models | |
CN103543327A (zh) | 一种测定光伏组件的发电功率的方法及设备 | |
CN110752622A (zh) | 一种配电网仿射状态估计方法 | |
CN103091581B (zh) | 一种基于电压崩溃指数的极限传输容量计算方法 | |
Karki et al. | Reliability modeling for evaluating the contribution of photovoltaics in electric power systems | |
CN103107533A (zh) | 适应大规模间隙式能源接入的输电网短期扩展规划方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190409 |
|
CF01 | Termination of patent right due to non-payment of annual fee |