基于杂交PSO算法的风光储微网日前优化调度(MATLAB实现)

2024-01-16 18:36

本文主要是介绍基于杂交PSO算法的风光储微网日前优化调度(MATLAB实现),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

        微网中包含:风电、光伏、储能、微型燃气轮机,以最小化电网购电成本、光伏风机的维护成本、蓄电池充放电维护成本、燃气轮机运行成本及污染气体治理成本为目标,综合考虑:功率平衡约束、燃气轮机爬坡约束、电网交换功率约束、储能装置约束、可控微电源出力约束等约束条件建立优化调度模型。

1.系统结构模型

1.1 光伏

         光伏电池的空载电压与光谱辐照度有关,与电池面积无关。在100MW/cm2的太阳光谱辐照度下,单晶硅光伏电池的空载电压为450——600mV,最高可达690mV。当入射光谱辐照度变化时,光伏电池的空载电压与入射光谱辐照度的对数成正比。环境温度升高时,光伏电池的空载电压将下降,一般温度每升高1℃,Uoc约下降2——3mV。
        由上述定义,可列出光伏电池等效电路中各变量的方程式如下:

式中,IO为光伏电池内部等效二极管PN结反向饱和电流,一般为常数;UD为等效二极管端电压;q为电子电荷;k为玻尔兹曼常量;T为热力学温度,A为PN结曲线常数。
        光伏电池的伏安特性与光照强度和电池温度有关,可以认为光伏电池的数学物理模型方程以光照强度和电池温度作为参变量。以下是环境因素对光伏电池数学物理模型的修正。
        当电池温度T、光照强度S不是参考光照强度和参考电池温度时,必须考虑它们的影响,并需要对光伏电池的数学物理模型进行修正。根据参考光照强度和参考电池温度下的Isc、Uoc、Im和Um,再将它们代入实用表达式,即可得到新的光照强度和新的电池温度下的伏安特性曲线。

1.2 风电

        风力发电机的输出功率可以用下式表示:

式中:Pwt为风力发电机的输出功率,\rho为空气密度,V为风速,R为风力发电机组的叶片半径,

Cp为风能利用系数。
        风力发电机的出力和风速的关系为:

式中:Prated为额定功率;Vrated为额定风速,Vcutin为切入风速,Vcutout为切出风速

1.3 微型燃气轮机

        选取燃氢微燃机作为氢电转换部分的能量耦合设备进行发电,该设备以氢⽓为燃料,燃烧产物只有水,实现零排放。
        其数学模型为:

式中:P_{MG}为微燃机输出功率;\eta _{M}^{}为微燃机热效率;H_{MT}为消耗的天然气量;C_{L}为低位热值常量。

1.4 储能装置模型

式中:S_{ES}为储能装置容量;\tau为储能损失系数;\eta _{ch}\eta _{dis}分别为充放电效率。

2. 目标函数

        目标函数为微网运行经济成本,主要包括:与主网交换功率成本,光伏、风机、蓄电池的维护成本,燃气轮机运行成本,污染气体治理成本,燃气轮机开机成本。

2.1 与主网交换功率成本

式中:C_{buy}C_{sell}分别对应与电网购电售电的分时电价;Nt为调度周期。

2.2 分布式电源维护成本

式中:P_{i}为第i个分布式电源的输出功率;K_{om,i}为第i个分布式电源的运行维护系数,单位为元/kWh;Ni为分布式电源总个数;

2.3 微型燃气轮机运行费用

式中:Cmt为燃气轮机燃料成本系数,单位为元/m3;fmt为燃气轮机燃料消耗系数,单位为m3/kWh;Cmts为燃气轮机开机成本,单位为元/次;S(t)为燃气轮机t时段的启停状态.。

2.4 污染气体治理费用

3.约束条件

3.1 系统功率平衡约束

3.2 分布式电源输出功率约束

3.3 燃气轮机的爬坡约束

式中:Rup、Rdown为燃气轮机输出功率最大上升速率与下降速率,单位为kW/h。

3.4 微网与主网交换功率约束

        微电网与主网的连接线有最大功率约束条件,故应满足该约束。

3.5 主网购售电状态互斥约束

        主网购售电的状态量,取值为0或1。

3.6 燃气轮机启停时间约束

式中:Ut-k+1为燃气轮机在t-k+1时段的启停状态;Tmup、Tmdown分别为燃气轮机的最小启动时间和最小关停时间。

3.7 储能装置约束

式中:\varepsilon为一个调度周期后蓄电池SOC变化的最大范围;\eta _{B}为自放电系数。

4. 求解算法

        采用PSO粒子群算法,用一种粒子来表示一个个体,每个粒子可视为N维搜索空间中的一个搜索个体,粒子的当前位置即为对应优化问题的一个候选解,粒子的飞行过程即为该个体的搜索过程.粒子的飞行速度可根据粒子历史最优位置和种群历史最优位置进行动态调整.粒子仅具有两个属性:速度和位置,速度代表移动的快慢,位置代表移动的方向。每个粒子单独搜寻的最优解叫做个体极值,粒子群中最优的个体极值作为当前全局最优解。不断迭代,更新速度和位置。最终得到满足终止条件的最优解。
        算法流程如下:
        ①初始化
        首先,我们设置最大迭代次数,目标函数的自变量个数,粒子的最大速度,位置信息为整个搜索空间,我们在速度区间和搜索空间上随机初始化速度和位置,设置粒子群规模,每个粒子随机初始化一个速度。
        ②个体极值与全局最优解
        定义适应度函数,个体极值为每个粒子找到的最优解,从这些最优解找到一个全局值,叫做本次全局最优解。与历史全局最优比较,进行更新。
        ③更新粒子的速度和位置

式中:ω为惯性因子,C1和C2为加速度常数,前者为每个粒子的个体学习因子,后者为每个粒子的社会学习因子。
        ④终止条件
(1)达到设定迭代次数;(2)代数之间的差值满足最小界限。
        本文在PSO粒子群算法的基础上进行改进,首先将传统固定惯性权重系数更改为根据全局最优点自适应调整惯性权重系数,即权重根据粒子的位置不同而动态变化。采用的非线性动态惯性权重系数公式为:

其中,f表示粒子实时的目标函数值;favgfmin分别表示当前所有粒子的平均值和最小目标值。从上面公式可以看出,惯性权重随粒子目标函数值的变化而变化,当粒子目标值分散时,减小惯性权重;当粒子目标值集中时,增大惯性权重。
        借鉴遗传算法中的杂交的概念,在每次迭代中,根据杂交率选取指定数量的粒子放入杂交池内,池内的粒子随机地两两杂交,产生同样数目的子代粒子,并用子代粒子代替父代粒子,具体计算公式和计算步骤,见《MATLAB智能算法》P132。

5.仿真结果

6.代码下载

        本文相关全部代码下载地址:基于杂交PSO算法的风光储微网日前优化调度(MATLAB实现)

这篇关于基于杂交PSO算法的风光储微网日前优化调度(MATLAB实现)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/613525

相关文章

SpringBoot3实现Gzip压缩优化的技术指南

《SpringBoot3实现Gzip压缩优化的技术指南》随着Web应用的用户量和数据量增加,网络带宽和页面加载速度逐渐成为瓶颈,为了减少数据传输量,提高用户体验,我们可以使用Gzip压缩HTTP响应,... 目录1、简述2、配置2.1 添加依赖2.2 配置 Gzip 压缩3、服务端应用4、前端应用4.1 N

SpringBoot实现数据库读写分离的3种方法小结

《SpringBoot实现数据库读写分离的3种方法小结》为了提高系统的读写性能和可用性,读写分离是一种经典的数据库架构模式,在SpringBoot应用中,有多种方式可以实现数据库读写分离,本文将介绍三... 目录一、数据库读写分离概述二、方案一:基于AbstractRoutingDataSource实现动态

Python FastAPI+Celery+RabbitMQ实现分布式图片水印处理系统

《PythonFastAPI+Celery+RabbitMQ实现分布式图片水印处理系统》这篇文章主要为大家详细介绍了PythonFastAPI如何结合Celery以及RabbitMQ实现简单的分布式... 实现思路FastAPI 服务器Celery 任务队列RabbitMQ 作为消息代理定时任务处理完整

Java枚举类实现Key-Value映射的多种实现方式

《Java枚举类实现Key-Value映射的多种实现方式》在Java开发中,枚举(Enum)是一种特殊的类,本文将详细介绍Java枚举类实现key-value映射的多种方式,有需要的小伙伴可以根据需要... 目录前言一、基础实现方式1.1 为枚举添加属性和构造方法二、http://www.cppcns.co

使用Python实现快速搭建本地HTTP服务器

《使用Python实现快速搭建本地HTTP服务器》:本文主要介绍如何使用Python快速搭建本地HTTP服务器,轻松实现一键HTTP文件共享,同时结合二维码技术,让访问更简单,感兴趣的小伙伴可以了... 目录1. 概述2. 快速搭建 HTTP 文件共享服务2.1 核心思路2.2 代码实现2.3 代码解读3.

MySQL双主搭建+keepalived高可用的实现

《MySQL双主搭建+keepalived高可用的实现》本文主要介绍了MySQL双主搭建+keepalived高可用的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,... 目录一、测试环境准备二、主从搭建1.创建复制用户2.创建复制关系3.开启复制,确认复制是否成功4.同

Java实现文件图片的预览和下载功能

《Java实现文件图片的预览和下载功能》这篇文章主要为大家详细介绍了如何使用Java实现文件图片的预览和下载功能,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... Java实现文件(图片)的预览和下载 @ApiOperation("访问文件") @GetMapping("

Spring Boot + MyBatis Plus 高效开发实战从入门到进阶优化(推荐)

《SpringBoot+MyBatisPlus高效开发实战从入门到进阶优化(推荐)》本文将详细介绍SpringBoot+MyBatisPlus的完整开发流程,并深入剖析分页查询、批量操作、动... 目录Spring Boot + MyBATis Plus 高效开发实战:从入门到进阶优化1. MyBatis

MyBatis 动态 SQL 优化之标签的实战与技巧(常见用法)

《MyBatis动态SQL优化之标签的实战与技巧(常见用法)》本文通过详细的示例和实际应用场景,介绍了如何有效利用这些标签来优化MyBatis配置,提升开发效率,确保SQL的高效执行和安全性,感... 目录动态SQL详解一、动态SQL的核心概念1.1 什么是动态SQL?1.2 动态SQL的优点1.3 动态S

使用Sentinel自定义返回和实现区分来源方式

《使用Sentinel自定义返回和实现区分来源方式》:本文主要介绍使用Sentinel自定义返回和实现区分来源方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录Sentinel自定义返回和实现区分来源1. 自定义错误返回2. 实现区分来源总结Sentinel自定