基于风光储能和需求响应的微电网日前经济调度(Python代码实现)

2023-05-16

目录

0 引言

1 计及风光储能和需求响应的微电网日前经济调度模型

1.1风光储能需求响应都不参与的模型

1.2风光参与的模型

1.3风光和储能参与模型

1.4 风光和需求响应参与模型

1.5 风光储能和需求响应都参与模型 

2 需求侧响应评价

2.1  负载率

2.2 可再生能源消纳率

2.3 用户舒适度

2.4 日最高负荷与日最高-最低负荷比率

3 算例实现

3.1算例介绍

3.2风光参与的模型求解

3.3 风光和储能参与的模型求解

3.5 风光储能和需求响应都参与模型求解

3.6 结果分析对比 

4 Python代码及算例数据


摘要:以实现经济性最优为目标,在得到新能源出力、负荷出力的日前预测上,考虑电网侧的实时电价。对风光储能需求响应都不参与、风光参与、风光和储能参与、风光和需求响应参与、风光储能和需求响应都参与五种调度方案进行探讨。对于需求响应调度,首先基于价格伸缩系数法的价格型需求响应,再加入激励型需求响应调度。对比各种调度情况下负载系数、新能源消纳率及用户舒适度。得到各种方案的供电构成图和供电成本,实验结果表明:风光储能和需求响应都参与微电网经济调度降低微电网的运行成本,并起到很好的消峰填谷的作用。

关键词:实时电价、风机、光伏、储能、需求响应、微电网经济调度


0 引言

近年来,微电网、清洁能源等已成为全球关注的热点。清洁能源在我国可持续发展战略中具有日益重要的地位,政府、企业、学术界均在清洁能源技术及其评价上做了大量工作[1]。

 微电网(Micro-Grid)日前经济调度问题是指考虑电网的分时电价基础上,对常规负荷、光伏出力、风机出力进行日前(未来 24 小时)预测,然后充分利用微网中的储能等可调控手段,使微电网运行的经济性最优[2-5]。

需求响应是电力需求侧管理在电力市场中的最新发展[6]。

 众多学者对此做了大量研究,文献[7]在考虑电价与负荷响应量相关性的基础上,建立了以运行成本最低、可再生能源消纳比例最高和用户满意度最好为目标的日前调度模型。文献[8]考虑风电功率以及可再生能源停运的不确定性,提出一种考虑风险的微电网日前随机优化调度方法。文献[9]以运行成本最小为优化目标、同时考虑储能,提出了一种微电网日前调度调度模型。文献[10]针对微电网新能源出力不确定的问题,提出了考虑电价激励需求响应下微电网日前优化调度方法。文献[11]以提高供电可靠性和能源利润为目的,构建考虑激励型DR和电池储能的优化模型,算例分析了激励型DR和电池储能对微网可靠性与盈利能力的影响。

本文以微电网风光、储能、与主网交换成本、需求侧响应成本为目标函数,考虑功率平衡约束、设备约束以及系统规划约束条件对此微电网进行优化调度,并通过峰值消减指数、可再生能源消纳率以及用户舒适度三个需求侧响应指标来量化需求侧响应效益。在允许微电网与主网进行功率交互的前提下,研究了微电网在不同情况下的调度方法。最后,通过实际算例分析了各种情况下的微电网经济优化及峰谷差最小运行方案。


1 计及风光储能和需求响应的微电网日前经济调度模型

1.1风光储能需求响应都不参与的模型

微网中储能和需求响应不作用,微电网与电网交换功率无约束,且无可再生能源,所以微电网中所需的电量只能从主网购买。这时候只需要满足功率平衡约束即可。

目标函数见式子(1)

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 

1.2风光参与的模型

若较负荷所需功率不足,不足部分向主网购买;若剩余,则剩余部分售卖给主网。光伏和风机可以舍弃部分。同时微电网与主网交换功率有约束。则目标函数见式子(3):

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 

1.3风光和储能参与模型

蓄电池的荷电状态应满足上下限约束,同时蓄电池的单位时间蓄电池充放电功率也有上下限。微电网与主网功率交换有约束,风光可以部分舍弃。在此种情况下:目标函数见式子(5)

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 

1.4 风光和需求响应参与模型

1.4.1 基于价格型需求响应

在经济学领域,通常采用价格伸缩系数来描述商品的需求对于价格敏感程度[12]。

客户用电量与电价之间的关系采用电价伸缩系数\xi来表示。\xi可以定义为客户用电量波动率与价格波动率的比值,如式(7):

1.4.2 基于激励型的需求响应

负荷可分为可控负荷和不可控负荷,不可控负荷为一些基础性用电设备,如照明用电、通信用电等,在调度时段内优先满足其供电;可控负荷为一些用电时段具有一定灵活性的用电设备,如洗衣机、可定时电饭煲等,在调度时段内采用激励型需求响应直接控制方式进行负荷转移,并给予补偿。

本文所研究的需求侧响应模型主要是以可控制负荷为研究对象。

可转移负荷主要是根据供需双方事先签订好的协议,由调度中心向用户发出信号,将用电高峰时期的某些比较灵活的负荷转移到用电低谷期或者新能源出力的高峰期,经过用户的响应之后,主动转移部分负荷。[15]

目标函数为式子(13)

 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

1.5 风光储能和需求响应都参与模型 


2 需求侧响应评价

电力需求侧响应是通过改变用户用电方式、提高终端用电效率等手段对终端电能使用进行控制,在满足用户电能使用效用的同时降低电力消耗微电网优化调度的主要目标是消减高峰负荷同时提高新能源消纳率。[16]因此可以提出三个指标来量化需求侧响应结果。

2.1  负载率

需求侧管理如果合理规划,可以减少一些高峰时段负荷,从而减轻操作员的压力。负载率可表示为:

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

2.2 可再生能源消纳率

需求侧管理如果合理规划,可以在可再生能源出力峰值时刻多转移一些负荷来消纳风光等能源,所以可再在能源的消纳率应该会增加可再生能源消纳率指数M_{rei}可表示为:

2.3 用户舒适度

用户各时段的用电量同用户的舒适度密切相关,一般来讲,考虑到用户对电价的自适应性,可以认为,在相对稳定的电价水平下,用户将会自主选择其舒适度最大的用电方式需求侧响应的投人会改变系统中负荷结构,使用户舒适度降低.因此用户舒适度也是需求侧响应的一个评价指标,其可表示为:

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 

2.4 日最高负荷与日最高-最低负荷比率

对于微电网电网而言,如果在其允许的范围内峰值负荷过大,会导致其转带备用容量不足,负荷波动率较高,容易影响该配电网的稳定性 。同时, 通过柯西-施瓦茨不等式证明了网损与最高-最低负荷比率成正比,即日内负荷的峰谷差越小,系统线损率越低。可见,缩小电网的高低负荷比率能够有效降低网损,提高系统负荷率,从而提升设备利用率。

因为这两指标可以从直观上判断出系统的供电能力,故也可以定义这两个规模目标为系统的显性供电能力,其值越小,系统的显性供电能力越好。日最高负荷与日最高-最低负荷比率如式(2.4.1):

3 算例实现

3.1算例介绍

一个含有风机、光伏、蓄电池以及负荷的微电网系统见示意图3.1。

风机的装机容量360kW,单位运维成本0.52元/kWh。光伏的装机容量 260kW,单位运维成本0.75元/kWh。蓄电池额定容量为700kWh,电池 SOC 运行范围 为[0.4,0.9],初始 SOC 值为 0.4,由充电至放电成本为 0.1 元/kWh,1个小时充放电功率最大为储能的20%。微网与主网允许交换功率不超过 200kW。一天的售电和购电价格见表1。

风机出力、光伏出力、常规负荷日前(未来 24 小时)预测见图2,新能源预测出力和净负荷,净负荷表示的是风机、光伏满足负荷后的剩余功率,如图3所示。

         

3.2风光参与的模型求解

在无风光情况下:微电网所需功率全部来自主网,功率等于负荷,在已知负荷和分时电价情况下,可以直接相乘求解。但是由图7我们可以看出,超过了微网与主网交换功率限额。

在有风光参与的情况下,由于风机功率、光伏功率、微电网与主网交换功率都是可变的,故直接相乘得不到结果,通过粒子群算法求解。

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

3.3 风光和储能参与的模型求解

要尽量满足用户负荷需求要,我们考虑未满足负荷用电量和,并与系统运行成本取权重作为新的目标函数。

通过制定策略:风力发电成本最低,所以对于可再生能源发电,优先发风机。然后我们采用允许弃风、弃光的方法,得到可再生能源的最优出力后,若较负荷所需功率不足,则由主网供电。还对风、光、储能和主网出力越限做了惩罚项处理。

3.4 风光和需求响应参与模型求解

对于需求响应参与的模型,我们先对负荷进行K-Means聚类算法,把负荷分为高峰、平段、低谷三个时段,如表2所示。通过1.4.1节的基于价格型需求策略得到需求响应后的负荷P',然后在调度时段内再采用激励型需求响应直接控制方式进行负荷转移,并给予补偿。

这里的出力策略和风光和储能参与的模型求解的策略方法一样,把储能出力换成可转入转出负荷。

 

3.5 风光储能和需求响应都参与模型求解

通过把3.3节和3.4节结合起来,得到如下风光储能和需求响应都参与模型

3.6 结果分析对比 

             

                    

                 

                 

由图7-12可以看出,当只有风光参与供电的模型负荷缺额量很大,当加入储能后负荷缺额量降低了很多。在只考虑风光和需求响应时【电价型】可以很明显的看出,很好的起到了消峰填谷的作用。在风光储能和需求响应都参与的模型中,各方面的效果都有很大改善,如表3-7所示。

可视化结果:

4 Python代码及算例数据

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

基于风光储能和需求响应的微电网日前经济调度(Python代码实现) 的相关文章

  • vscode代码提交到gittee码云 第一次提交方法

    学习3 xff1a 今天是第一次将vscode代码提交到gittee xff0c 废话不多说 xff0c 直接上方法 xff1a 查看git仓库 gt git status 将当前项目文件初始化为仓库 如果当前文件夹不是git仓库 xff0
  • 明火烟雾目标检测项目部署(YoloV5+Flask)

    明火烟雾目标检测项目部署 文章目录 明火烟雾目标检测项目部署1 拉取Docker PyToch镜像2 配置系统环境2 1 更换软件源2 2 下载vim2 3 解决vim中文乱码问题 3 运行项目3 1 拷贝项目到容器中3 2 安装项目所需的
  • 操作系统实践课作业(南航)

    操作系统实践课作业 xff08 南航 xff09 文章目录 操作系统实践课作业 xff08 南航 xff09 1 job21 1 main c1 2 math c1 3 Makefile 2 job32 1 myecho c2 2 myca
  • 在Linux系统下安装Neo4j图数据库

    在Linux系统下安装Neo4j图数据库 文章目录 在Linux系统下安装Neo4j图数据库1 Java JDK1 1 安装1 2 查看安装路径 2 Neo4j2 1 下载2 2 拷贝到容器中2 3 修改neo4j conf配置文件2 4
  • 大数定律 与 中心极限定理 的理解

    目录 1 大数定律 2 中心极限定理 1 大数定律 当样本的数量足够大时 xff0c 样本的统计特性就可以近似代表总体的统计特性 大数 是指样本的数量足够大或者试验的次数足够多 2 中心极限定理 设总体为 为总体的 N 个样本集 xff0c
  • 操作系统实践05—文件描述符和系统调用

    操作系统实践05 文件描述符和系统调用 文章目录 操作系统实践05 文件描述符和系统调用1 概念1 1 文件描述符1 2 系统调用1 3 例子 2 内核实现2 1 file结构体2 2 文件描述符表2 3 进程控制块2 4 私有的文件描述符
  • 医疗问答机器人项目部署

    医疗问答机器人项目部署 文章目录 医疗问答机器人项目部署1 拉取TensorFlow镜像2 配置系统环境2 1 更换软件源2 2 下载vim2 3 解决vim中文乱码问题2 4 安装Neo4J图数据库2 5 安装网络工具包 3 运行项目3
  • SimpleITK学习

    SimpleITK学习 文章目录 SimpleITK学习1 SimpleITK ReadImage path 2 SimpleITK GetArrayFromImage itk img 3 itk img GetOrigin 4 itk i
  • 【Docker】服务器部署项目

    服务器部署项目 文章目录 服务器部署项目1 远程连接服务器2 在Linux系统上安装Docker2 1 卸载旧版本2 2 使用 APT 安装2 3 安装Docker2 4 使用脚本自动安装2 5 启动Docker2 6 测试 Docker
  • 计算机网络04—网络层

    网络层 学习参考资料 xff1a 湖南科技大学 计算机网络谢希仁 计算机网络 xff08 第7版 xff09 文章目录 网络层1 概述1 1 IP协议及配套协议 2 两种服务2 1 面向连接的虚电路服务2 2 无连接的数据报服务2 3 对比
  • torch.nn学习

    torch nn学习 文章目录 torch nn学习1 卷积层1 1 Conv2d 2 池化层2 1 MaxPool2d2 2 MaxUnpool2d2 3 AvgPool2d 3 代码实践3 1 Inception Module3 2 R
  • 深度学习基础知识点【更新中】

    深度学习基础知识点 文章目录 深度学习基础知识点1 数据归一化2 数据集划分3 混淆矩阵4 模型文件5 权重矩阵初始化6 激活函数7 模型拟合8 卷积操作9 池化操作10 深度可分离卷积11 转置卷积 1 数据归一化 过大的输入数据未归一化
  • VS Code配置C/C++环境

    VS Code配置C C 43 43 环境 文章目录 VS Code配置C C 43 43 环境1 下载Visual Studio Code2 下载MinGW3 VS Code设置3 1 下载插件3 2 新建工作区3 3 C 43 43 环
  • 计算机网络05—运输层

    运输层 学习参考资料 xff1a 湖南科技大学 计算机网络谢希仁 计算机网络 xff08 第7版 xff09 文章目录 运输层1 概述1 1 两个主要协议1 2 端口 2 用户数据报协议UDP3 传输控制协议TCP3 1 概述3 2 可靠运
  • 【2022春招研发】字节笔试记录(测试方向)

    20220410字节笔试 测试方向 文章目录 20220410字节笔试 测试方向一 编程题2道 xff08 50分 xff09 二 单选题10道 xff08 20分 xff09 三 多选题10道 xff08 30分 xff09 一 编程题2
  • 浏览器主页被劫持篡改了怎么办

    就想下载个驱动 xff0c 结果一通操作把我的 Edge 浏览器主页篡改成了 桔梗网 xff0c 就下面这个网站 算了不喷它了 xff0c 来说说怎么改回去吧 其他浏览器的修改方式相同 找到 Microsoft Edge 浏览器的桌面快捷方
  • 【2022春实习】百度笔试记录(机器学习/数据挖掘/自然语言)

    20220412百度笔试 机器学习 数据挖掘 自然语言 文章目录 20220412百度笔试 机器学习 数据挖掘 自然语言一 选择题30道 xff08 60分 xff09 二 问答题1道 xff08 20分 xff09 三 系统设计题1道 x
  • 【算法工程师】华为技术面面试记录

    20220419华为技术面 面试岗位是算法工程师 文章目录 20220419华为技术面1 自我介绍 2 算法题3 专业知识3 1 数据结构3 2 计算机网络3 3 操作系统3 4 设计模式3 5 机器学习3 6 其他 4 提问环节 1 自我
  • 操作系统实践06—线程

    操作系统实践06 线程 文章目录 操作系统实践06 线程1 创建线程1 1 原型1 2 线程参数1 3 参数类型1 4 例子一1 5 例子二 2 等待线程2 1 原型2 2 线程返回值2 3 例子一2 4 例子二 3 线程互斥3 1 初始化
  • VS Code指定扩展安装位置

    VS Code指定扩展安装位置 默认情况下 xff0c Windows vscode的安装路径为C Users 用户名 vscode extensions 如果想要自定义扩展的安装路径 xff0c 无法直接在vscode中修改 但是 xff

随机推荐

  • Ubuntu 网络配置顺序:(Ubuntu 16.4)

    网络配置顺序 xff1a xff08 Ubuntu 16 4 xff09 1 xff0c 网卡硬件 xff08 硬件 vm DHCP用NAT直接到物理网 xff0c 静态用桥接通过本地网络链接转发 xff09 xff0c 2 xff0c 系
  • C语言实现汉诺塔问题

    目录 一 程序 1 实现代码 2 程序执行结果 二 背景 1 汉诺塔问题描述 2 直观理解 3 思考过程 三 函数执行过程 举例 n 61 2 四 总结 递归问题 一 程序 1 实现代码 include lt stdio h gt 函数名
  • ubuntu 系统下运行DQN Flappy bird

    github链接 xff1a xff08 两个项目均可运行成功 xff0c 推荐第二个 xff0c 只有200多行 xff0c 好入门 xff09 https github com songrotek DRL FlappyBird http
  • eclipse中如何查看程序源码

    在eclipse中如何查看程序源码 1 查看方法 xff0c 按下Ctrl键同时 xff0c 用鼠标指向自己想要查看源码的关键字 xff0c 该字会出现下划线 xff0c 单击左键 xff0c 就可以进入源码查看了 2 在此前 xff0c
  • 基于51单片机的智能窗帘proteus仿真数码管显示

    硬件设计 该硬件设计是基于51单片机为MCU xff0c ADC采用ADC0804 xff0c 电机驱动芯片是L298 xff0c 显示部分采用的是4位数码管 ADC0804芯片的简介 xff1a 工作电压 xff1a 43 5V xff0
  • Mac VSCode 配置Gitee |使用Sourcetree关联Gitee仓库

    目录 一 使用Git的基本操作 xff08 如果只是想看如何上传到gitee的直接跳过这个 xff09 xff1a 1 安装插件GitLens 2 初始化 3 git文件 4 创建一个示例文件 5 操作日志 6 操作日志查看 7 版本回退
  • 马尔可夫链 以及 隐马尔可夫模型(HMM)

    背景 xff1a 马尔可夫过程 xff08 Markov process xff09 是一类随机过程 它的原始模型马尔科夫链 xff0c 由俄国数学家A A Markov于1907年提出 马尔可夫过程是研究离散时间动态系统状态空间的重要方法
  • AIOps探索:基于VAE模型的周期性KPI异常检测方法——VAE异常检测

    AIOps探索 xff1a 基于VAE模型的周期性KPI异常检测方法 VAE异常检测 参考文章 xff1a xff08 1 xff09 AIOps探索 xff1a 基于VAE模型的周期性KPI异常检测方法 VAE异常检测 xff08 2 x
  • 51单片机系列--蜂鸣器

    工作原理 蜂鸣器发声原理是电流通过电磁线圈 xff0c 使电磁线圈产生磁场来驱动振动膜发声的 xff0c 因此需要一定的电流才能驱动它 51单片机IO口输出的TTl电流无法驱动蜂鸣器 xff0c 故而蜂鸣器内部需要一个三极管来进行电流放大
  • debian9.8安装网卡驱动

    一 挂载本地镜像源 xff08 参考下面博客 xff09 debian9 8添加iso为本地源 weixin 46027366的博客 CSDN博客 二 配置安装所需环境 1 xff09 安装gcc和cmake xff0c 命令行输入如下命令
  • 浮点型数据的输入和输出(C语言)

    目录 1 浮点型数据的输入 1 1 单精度输入 1 2 双精度和长双精度 2 浮点型数据的输出 2 1 浮点数的默认输出 2 2 指定输出格式 m n f 2 3 输出示例 3 案例 3 1 案例 1 浮点型数据的输入 1 1 单精度输入
  • 开学送给她的礼物(Python实现)

    目录 1 卿为朝朝暮暮 2 情感起伏 3 礼物赠送 4 Python之实现turtle 1 卿为朝朝暮暮 先手抄一遍 xff0c 然后再键盘敲出来 xff1a 飞鸟集中的一句话 xff0c 改编的一首诗是过样的 浮世万千 xff0c 吾爱有
  • 最详细matlab 2018a安装教程步骤.

    链接 xff1a https pan baidu com s 1XjfAKeFY otNy7HfGhYQCw 提取码 xff1a cmzv 来自百度网盘超级会员V3的分享 1 鼠标右击 Matlab R2018a Win64 压缩包 xff
  • 牵着她——表白不成功算我输(Python实现)

    目录 1 牵着她的手一直走下去 2 一首小情诗送给甜甜的她 3 历史总结的哲学想法 4 表白不成功算我输 xff08 Python代码 xff09 1 牵着她的手一直走下去 今天牵着她的手 xff0c 她很贴心 一起并肩赏樱花 x1f338
  • Python|十五个超级炫酷代码

    x1f96c x1f96c x1f96c 欢迎来到本博客 x1f60a x1f60a x1f60a 本次博客内容将分享几个超级炫酷的Python代码和前端 x1f4dd 目前更新 xff1a x1f31f x1f31f x1f31f 炫酷炫
  • 【无人机】四轴无人机的轨迹进行可视化和动画处理(Matlab代码实现)

    x1f4cb x1f4cb x1f4cb 本文目录如下 xff1a 目录 1 概述 2 运行结果 3 参考文献 4 Matlab代码实现 1 概述 随着传感器检测技术 智能控制技术和材料技术的快速发展 四轴无人机及其配套系统的发展越来越成熟
  • 【无人机】四旋翼飞行器控制、路径规划和轨迹优化(Matlab代码实现)

    x1f4a5 x1f4a5 x1f4a5 x1f49e x1f49e x1f49e 欢迎来到本博客 x1f4a5 x1f4a5 x1f4a5 x1f3c6 博主优势 xff1a x1f31e x1f31e x1f31e 博客内容尽量做到思维
  • Viterbi算法(维特比算法)

    维特比算法背景 xff1a 安德鲁 维特比 xff08 Andrew J Viterbi xff09 xff0c CDMA之父 xff0c IEEE Fellow xff0c 高通公司创始人之一 xff0c 高通首席科学家 他开发了卷积码编
  • 【信号处理】卡尔曼(Kalman)滤波(Matlab代码实现)

    x1f468 x1f393 个人主页 xff1a 研学社的博客 x1f4a5 x1f4a5 x1f49e x1f49e 欢迎来到本博客 x1f4a5 x1f4a5 x1f3c6 博主优势 xff1a x1f31e x1f31e x1f31e
  • 基于风光储能和需求响应的微电网日前经济调度(Python代码实现)

    目录 0 引言 1 计及风光储能和需求响应的微电网日前经济调度模型 1 1风光储能需求响应都不参与的模型 1 2风光参与的模型 1 3风光和储能参与模型 1 4 风光和需求响应参与模型 1 5 风光储能和需求响应都参与模型 2 需求侧响应评