本发明涉及一种海岛智能微电网多源协调控制方法。
背景技术：
社会经济发展对于能源的需求越来越大，越来越严重的环境污染使得传统能源在应用过程中面临瓶颈，如何解决能源匮乏和环境污染是当今社会面临的重要问题。十八大报告提出“提高海洋资源开发能力，发展海洋经济，保护海洋生态，坚决维护国家海洋权益，建设海洋强国”。海岛作为海洋的重要组成部分，建设开发海岛是中国实施海洋战略的重要组成部分，是壮大海洋经济、拓展发展空间的重要依托。我国海岛众多，能源保障是海岛开发的基本要求，以此，开发海岛及其周围丰富的可再生能源，构建基于可再生能源分布式发电系统的海岛智能微电网，可以解决海岛开发过程中的能源问题。海洋能源的开发对于海洋环境的保护，促进节能减排也具有重要意义。
智能微电网是为了解决分布式发电系统并网接入而提出了一种新的电网组织形式。微电网集分布式发电系统、储能系统和负荷于一体，通过自身的控制实现孤岛运行和并网运行，并可在两种运行方式间平滑切换。微电网是提高分布式发电系统利用效率的有效途径，避免了大量分布式发电系统接入对电力系统的冲击，并能够提高电网的供电能力。海洋蕴含大量的可开发的可再生能源，其中风能、太阳能、波浪能、潮流能、潮汐能等能源相对开发难度较低。目前基于上述能源的分布式发电系统不断涌现为海岛智能微电网的应用提供了条件。
随着海岛智能微电网的建设，基于分布式发电系统运行特性的微电网优化运行提上日程。当多种能源接入微电网时，需要对各种可再生能源根据其运行特性进行综合调度。国内外学者在这方面进行了大量的研究，从经济性、可靠性、环保性等不同角度对微电网多种能源进行协调控制，提高微电网运行的可靠性和经济性，从而提高可再生能源分布式发电系统的利用效率。
技术实现要素：
本发明为了解决上述问题，提出了一种海岛智能微电网多源协调控制方法，本发明综合考虑各种可再生能源分布式发电系统的运行特性，以光伏、潮流能、柴油发电机和储能进行协调控制，提供一种海岛智能微电网多源协调控制方法。该方法能够为微电网的综合运行提供指导，提高微电网运行的经济性和可靠性。
为了实现上述目的，本发明采用如下技术方案：
一种海岛智能微电网多源协调控制方法，考虑海岛可再生能源分布式发电系统的运行特性，确定微电网网络拓扑结构，进行输出功率的预测，构建经济性、可靠性为目标函数的优化运行模型，基于粒子群算法，在约束条件下求解优化运行模型的全局最优解和个体最优解，确定最优值以对微电网多源优化调度。
优选的，所述可再生能源分布式发电系统具体包括光伏发电系统和潮流能发电系统。
优选的，当可再生能源分布式发电系统的能源不足时，通过柴油发电系统以及储能系统进行调节优化控制。
优选的，进行输出功率的预测时，利用人工神经网络方法光伏、潮流能和负荷预测模型预测光伏系统、潮流能发电系统未来一段时间的输出功率和负荷功率。
优选的，微电网优化运行模型的目标函数为微电网运行成本、微电网运行环境折算成本和供电可靠性成本之和最小。
进一步的，微电网运行成本包括运行维护成本、能源消耗成本和各机组的启动成本。
进一步的，微电网运行环境折算成本为所有污染物的排放量与其对应的环境折算成本乘积之和。
进一步的，供电可靠性成本包含微电网优化运行过程中用户失去负荷的损失，由电价、用户费用补偿系数和相应时刻的负荷中断功率确定。
优选的，所述约束条件包括：微电网各分布式电源发出的功率之和与用户所需求的功率相互匹配，每个微电源的实际输出功率任意时刻都要在其容量的上下限之内，每个电源的开机、关机时间均大于等于各自的最小开机和关机时间，调度时段内的启停次数要小于等于最大启停次数，或/和储能蓄电池在单位时间充放电过程中充电、放电功率不超过最大容量与设定阈值的乘积。
优选的，所述基于粒子群算法的方法为，将微电网中需要优化的分布式发电系统输出功率为粒子群每个粒子的每一维，确定最大迭代次数、优化参数搜索范围和变邻域搜索区间，进行目标函数的适应度计算，获得粒子群算法个体最优解和全局最优解，并进行粒子速度和位置的更新，对全局最优解和个体最优解进行局部的混沌变邻域搜索，根据邻域搜索确定新的个体最优解和全局最优解位置，直到满足迭代次数。
进一步的，混沌变邻域搜索时，利用logistic函数构造混沌变量，将产生的混沌变量映射为变量变化值，确定邻域半径和个体最优解或者全局最优解位置偏移值。
进一步的，混沌变邻域搜索时，邻域半径根据粒子位置的最大和最小值变化而变化，邻域半径随着迭代次数的增加而非线性的减小。
进一步的，确定新的个体最优解位置和全局最优解位置时，分别对比邻域搜索前后个体最优解和全局最优解的适应度，将适应度高的粒子确定为最新的个体最优解和全局最优解。
本发明的有益效果为：
(1)本发明提供了一种优化协调控制方法，综合利用海岛各种可再生能源，同时又不破坏海岛生态，为解决海岛供电和提高海岛可再生能源的综合开发提供保障；
(2)本发明能够充分考虑海岛可再生能源分布式发电系统的运行特性，结合柴油发电系统和储能系统，实现微电网内部分布式发电系统的协调控制，达到微电网优化运行的目的，从而提高系统运行的可靠性和经济性；
(3)本发明计算过程简单，且综合了经济型、可靠性和环保性的角度对为电网多种能源进行调控，提高微电网运行的可靠性和经济性，从而提高可再生能源分布式发电系统的利用效率。
附图说明
图1为本发明的流程示意图。
具体实施方式：
下面结合附图与实施例对本发明作进一步说明。
一种海岛智能微电网多源协调控制方法，包括以下步骤：
(1)构建微电网网络拓扑结构，对微电网中的分布式发电系统进行功率预测，确定光伏系统、潮流能系统以及负荷的预测情况；
(2)构建光伏系统、潮流能系统的模型，建立以微电网运行的经济性、可靠性为目标函数的微电网优化运行模型；
(3)构建功率平衡、输出功率、最小关机和最小运行时间、启停次数和储能蓄电池状态的约束条件；
(4)基于粒子群算法，在约束条件下求解微电网优化运行模型的全局最优解和个体最优解，进行粒子更新，输出最优值；
(5)将最优值作为微电网运行各分布式发电系统的最优出力进行微电网运行优化调度，实现微电网的多源协调优化控制。
所述步骤(1)中，微电网网络的分布式发电系统，具体包括光伏发电系统、潮流能发电系统和柴油机发电系统，在优化调度时考虑光伏发电系统、潮流能发电系统，当电源不足时再通过柴油发电系统以及储能系统进行调节优化控制。
所述步骤(1)中，利用人工神经网络方法光伏、潮流能和负荷预测模型预测光伏系统、潮流能发电系统未来一段时间的输出功率和负荷功率。
所述步骤(2)中，微电网优化运行模型的目标函数为微电网运行成本、微电网运行环境折算成本和供电可靠性成本之和最小。
所述步骤(2)中，微电网运行成本包括运行维护成本、能源消耗成本和各机组的启动成本。
所述步骤(2)中，微电网运行环境折算成本为所有污染物的排放量与其对应的环境折算成本乘积之和。
所述步骤(2)中，可靠性成本包含微电网优化运行过程中用户失去负荷的损失，由电价、用户费用补偿系数和相应时刻的负荷中断功率确定。
所述步骤(3)中，约束条件具体包括：
(I)微电网各分布式电源发出的功率之和与用户所需求的功率相互匹配；
(II)每个微电源的实际输出功率任意时刻都要在其容量的上下限之内；
(III)每个电源的开机、关机时间均大于等于各自的最小开机和关机时间；
(IV)调度时段内的启停次数要小于等于最大启停次数；
(V)储能蓄电池在单位时间充放电过程中充电、放电功率不超过最大容量与设定阈值的乘积。
所述步骤(4)中，具体步骤包括：
(4-1)微电网中需要优化的分布式发电系统输出功率为粒子群每个粒子的每一维，确定最大迭代次数、优化参数搜索范围、变邻域搜索区间；
(4-2)根据初始粒子位置和飞行速度进行目标函数的适应度计算，获得粒子群算法个体最优解和全局最优解，并进行粒子速度和位置的更新；
(4-3)对全局最优解和个体最优解进行局部的混沌变邻域搜索，根据邻域搜索确定新的个体最优解和全局最优解位置；
(4-4)根据确定的个体最优解和全局最优解进行粒子更新，根据最大迭代次数判断算法是否终止，如果未达到最大迭代次数，返回步骤(4-2)，如果达到迭代次数，计算终止输出结果。
所述步骤(4-3)中，利用logistic函数构造混沌变量，将产生的混沌变量映射为变量变化值，确定邻域半径和个体最优解或者全局最优解位置偏移值。
所述步骤(4-3)中，邻域半径根据粒子位置的最大和最小值变化而变化，邻域半径随着迭代次数的增加而非线性的减小。
所述步骤(4-4)中，确定新的个体最优解位置和全局最优解位置，分别对比邻域搜索前后个体最优解和全局最优解的适应度，将适应度高的粒子确定为最新的个体最优解和全局最优解。
如图1所示，本发明提供了一种智能海岛微电网多源协调控制方法，所述方法包含以下步骤：
步骤1：构建基于分布式发电系统运行特性的微电网多源协调优化控制模型。包含的分布式发电系统有光伏发电系统、潮流能发电系统、柴油机发电系统三种分布式发电系统。其中，光伏发电系统、潮流能发电系统为可再生能源分布式发电系统，在优化调度时优先考虑，当电源不足时通过柴油发电系统以及储能系统进行调节优化控制，另外在紧急情况下可切除一般负荷来满足电网运行的需求。
步骤2：采用基于混沌变量改进粒子群优化算法对分布式发电系统进行优化协调控制，达到系统最优运行
在步骤1中，首先确定微电网网络拓扑结构，光伏系统、潮流能系统以及负荷的预测情况。在实际运行过程中需要对微电网中的分布式发电系统进行功率预测，利用人工神经网络光伏、潮流能和负荷预测模型预测光伏系统、潮流能发电系统未来一天的输出功率，未来一天的负荷功率。
建立微电网优化运行模型，构建以微电网运行的经济性、可靠性为目标函数的微电网优化运行模型。
光伏发电系统的模型为：
潮流能发电系统模型为：
其中，PPV光伏发电系统功率，PCT为潮流能发电系统输出功率，G(t)为光照强度，T(t)为温度，GSTC、PSTC、TSTC分别为标准测试环境下(1000W/m2,25℃)下的光照强度、最大输出功率和温度。Cp为水轮机获能系数，ρ为海水的密度，S为水轮机叶轮的扫截面积，v为潮流上游流速。
微电网优化协调优化模型目标函数为：min F(x)＝[FCost(x),FEnvi(x),FReli(x)]T，其中FCost(x)为微电网运行成本，FEnvi(x)为微电网运行环境折算成本，FReli(x)为供电可靠性成本。
微电网运行成本由COM为运行维护成本，CF为能源消耗成本，CSC为各机组的启动成本组成，可表示为：其中，运行维护成本COM正比于机组的装机容量柴油发电机燃料费用成本表示不间断电源机组的启停费用成本取决于上一次机组停运的时间，机组的启动成本CSC表示为：
其中，NG为各发电单元的数量，为第i种电源的装机容量，Komi为第i种电源的运行维护系数；为t时段内的第i个柴油机所需燃料，为发电机启停状态，Ci为燃料价格；NSTOP为电源的启动次数，σi为电源i的热启动费用，δi为电源i的冷启动费用，τi为电源i的冷却时间常数，Toff,i为i电源已停运时间；
考虑到微网的环境效益，将环境损失作为微网优化运行的目标，表示为
可靠性成本包含微电网优化运行过程中用户失去负荷的损失，
其中，CE为环保折算费用，Yj为第j种污染物的排放量，Dj为第j种污染物的环境折算成本；M为中断负荷数量，C0为电价，βi为用户费用补偿系数，为t时刻负荷中断功率。负荷中断功率由负荷容量SIL,i和负荷关断标记ITIL,i决定，并满足
3)功率平衡约束。微电网各分布式电源发出的功率之和与用户所需求的功率相互匹配，其中，为t时刻用户需求功率；分别是潮流、光伏、储能、柴油机在t时刻所发出的功率。
输出功率约束。每个微电源的实际输出功率任意时刻都要满足其容量的上下限约束：Pi,min≤Pi≤Pi,max,Pi,min为电源i的最小输出功率，Pi,max为电源i的最大输出功率。
最小关机和最小运行时间约束。Ti,off≥Ti,off,min，Ti,on≥Ti,on,min。Ti,off，Ti,on分别为电源i的关机、开机时间，Ti,off,min，Ti,on,min分别为电源i的最小关机和开机时间。
启停次数约束。Ni,max为调度时段内最大启停次数。
储能蓄电池状态约束。考虑到其使用寿命的影响，其单位时间充放电过程中，不能超过其最大容量20％，要满足其充放电约束及能量状态约束：
P+为单位小时内的充电功率，P-为单位小时内的放电功率，EBT为蓄电池的最大容量，Δt为单位时间，为时段内的剩余容量，ηch为蓄电池充电效率，tch为蓄电池充电时间，ηdis为蓄电池放电效率，tdis为蓄电池放电时间。
建立系统优化模型后，利用步骤2进行求解：
步骤2-1：首先，构造粒子群结构，每个粒子的每一维为微电网中需要优化的分布式发电系统输出功率X＝[x1,x2,...xn]。初始化粒子群优化算法参数包括最大迭代次数、优化参数搜索范围、变邻域搜索区间；
步骤2-2：根据初始粒子位置和飞行速度进行目标函数的适应度计算，获得粒子群算法个体最优解xi,pbest和全局最优解xgbest，并进行粒子速度和位置的更新，更新的方程为：
其中，是第i个粒子在t时刻的速度，是第i个粒子在t时刻的位置，xi,pbest第i个粒子在t时刻的个体最优极值，gbest为所有粒子的全局最优极值，r1，r2是介于(0，1)的随机数，c1，c2是学习因子，w是粒子群的动态权重值。
步骤2-3：利用logistic函数ui,j+1＝4ui,j(1-ui,j)产生一个混沌变量，其中ui,j初始混沌变量取值在(0,1)之间，并不等于0.25、0.5、0.75，ui,j+1位下一代混沌变量。
步骤2-4：将产生的混沌变量映射为变量变化值Δxi＝-β+2β·ui,j+1，其中β为邻域半径，Δxi为个体最优解或者全局最优解位置偏移值。其中，邻域半径的变化为b＝(xi,max-xi,min)e-t，xi,min、xi,max为粒子位置的最大和最小值，邻域半径随着迭代次数的增加而非线性的减小。
步骤2-5：根据邻域搜索确定新的个体最优解和全局最优解位置。利用式确定新的个体最优解位置和全局最优解位置，分别对比邻域搜索前后个体最优解和全局最优解的适应度，将适应度高的粒子确定为最新的个体最优解和全局最优解。
步骤2-6：根据确定的个体最优解和全局最优解进行粒子更新，根据最大迭代次数判断算法是否终止，如果未达到最大迭代次数，返回步骤2，如果达到迭代次数，计算终止输出结果。
步骤3：根据计算得到的未来一天微电网运行各分布式发电系统的最优出力进行微电网运行优化调度，实现微电网的多源协调优化控制。
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。