CN115173489B - 一种基于二分法的火电集群调度方法及系统 - Google Patents
一种基于二分法的火电集群调度方法及系统 Download PDFInfo
- Publication number
- CN115173489B CN115173489B CN202210907585.XA CN202210907585A CN115173489B CN 115173489 B CN115173489 B CN 115173489B CN 202210907585 A CN202210907585 A CN 202210907585A CN 115173489 B CN115173489 B CN 115173489B
- Authority
- CN
- China
- Prior art keywords
- thermal power
- sum
- target
- output
- coal consumption
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 56
- 239000003245 coal Substances 0.000 claims abstract description 89
- 230000008859 change Effects 0.000 claims abstract description 23
- 230000009194 climbing Effects 0.000 claims description 11
- 238000010586 diagram Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000012360 testing method Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000001965 increasing effect Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 230000003442 weekly effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
-
- 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/10—The dispersed energy generation being of fossil origin, e.g. diesel generators
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Power Engineering (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于二分法的火电集群调度方法及系统,获取煤耗微增率上限和下限,经二分处理得到煤耗微增率目标;初始化火电出力变量之和和负荷改变量,利用得到的煤耗微增率目标确定火电出力变量之和;将火电机组出力变化量之和与负荷改变量的大小进行判断,安排火电机组出力,对获得的火电机组出力情况进行传递并返回初始化,进行火电机组下一时刻的出力安排,实现火电集群调度。本发明结合煤耗微增率建立基于二分法的火电经济调度方法,能够快速计算出火电经济调度的出力安排,找出一种较为经济的电力系统运行方式。
Description
技术领域
本发明属于电力系统技术领域,具体涉及一种基于二分法的火电集群调度方法及系统。
背景技术
随着社会的快速发展,人们对能源的需求不断上涨。而电力作为主要的终端能源形式,承担着能源供给的重要工作。为了使能源被高效地利用,研究者们也投入了大量的精力在解决如何找到一种较为经济的电力系统运行方式。
实际上对于电力系统的经济调度,已经有了大量的研究。简单来说,电力系统的经济调度就是通过调整火电机组的出力,找到一种在满足电力电量平衡的情况下系统运行成本最低的出力大小,这就是我们常说的以社会福利最大化的经济调度。工业或学界通常是将电力系统运行调度建立成为一个线性优化问题,利用单纯型法或内点法进行求解。
但对于线性规划,随着问题规模的增大,其求解的难度也会变得更大,而这将产生没有必要的算力负担和时间成本。此外,当一种规划方案需要利用经济调度的办法进行评估的时候,规划者同样不希望花费大量的时间和算力在结果评估的工作上。从算力消耗上和计算时间的角度上来看,求解线性优化问题就会存在一定的不足。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足,提供一种基于二分法的火电集群调度方法及系统,能够快速计算出火电经济调度的出力安排,克服由于所求解的问题规模过大而造成不必要的算力负担和时间成本的提高。
本发明采用以下技术方案:
一种基于二分法的火电集群调度方法,获取煤耗微增率上限IRupper和下限IRlower,经二分处理得到煤耗微增率目标IRtarget;初始化火电出力变量之和ΔPt g,sum和负荷改变量ΔPd,sum,利用得到的煤耗微增率目标IRtarget得到火电出力变量之和ΔPt g,sum;将火电机组出力变化量之和ΔPt g,sum与负荷改变量ΔPt d,sum的大小进行判断,当|ΔPt g,sum-ΔPt d,sum|<limit时,limit为设定阈值,安排火电机组出力,对获得的火电机组出力情况进行传递并返回初始化,进行火电机组下一时刻的出力安排,实现火电集群调度。
具体的,煤耗微增率目标IRtarget为:
其中,IRtarget为煤耗微增率目标。
具体的,火电出力变量之和ΔPt g,sum为:
其中,Ωg为火电集合,为时刻t火电机组i按照煤耗微增率目标IRtarget进行调整的出力变化量。
进一步的,时刻t火电机组i按照煤耗微增率目标IRtarget进行调整的出力变化量为:
其中,为时刻t-1火电机组i的煤耗微增率;为火电机组i的出力上限;为时刻t-1火电机组i的出力;RUi为火电机组i的上爬坡上限;Pi g为火电机组i的出力下限;RDi为火电机组i的下爬坡上限。
具体的,负荷改变量ΔPd,sum为:
其中,Ωd为负荷集合;为时刻t负荷i的负荷需求;为时刻t-1负荷i的负荷需求。
具体的,当|ΔPt g,sum-ΔPt d,sum|≥limit时,对煤耗微增率目标IRtarget进行修正,根据修正后的煤耗微增率目标IRtarget重新计算火电出力变量之和ΔPt g,sum,返回重新进行调整。
进一步的,对煤耗微增率目标进行修正IRtarget具体为:
其中,IRupper为设定的煤耗微增率上限,IRlower为设定的煤耗微增率下限。
进一步的,修正后的煤耗微增率目标IRtarget为:
具体的,t时刻火电机组i的出力安排为:
其中,为时刻t-1火电机组i的出力;为时刻t-1火电机组i的煤耗微增率;为火电机组i的出力上限;RUi为火电机组i的上爬坡上限;Pi g为火电机组i的出力下限;RDi为火电机组i的下爬坡上限。
第二方面,本发明实施例提供了一种基于二分法的火电集群调度计算系统,包括:
二分模块,获取煤耗微增率上限IRupper和下限IRlower,经二分处理得到煤耗微增率目标IRtarget;
初始模块,初始化火电出力变量之和ΔPt g,sum和负荷改变量ΔPd,sum,利用二分模块得到的煤耗微增率目标IRtarget得到火电出力变量之和ΔPt g,sum;
调度模块,将初始模块得到的火电机组出力变化量之和ΔPt g,sum与负荷改变量ΔPt d,sum的大小进行判断,当|ΔPt g,sum-ΔPt d,sum|<limit时,安排火电机组出力,对获得的火电机组出力情况进行传递,返回初始模块,进行火电机组下一时刻的出力安排,实现火电集群调度。
与现有技术相比,本发明至少具有以下有益效果:
本发明一种基于二分法的火电集群调度方法,基于等煤耗微增率点为火电机组最优运行点的基本思路,设计一种搜寻最优煤耗微增率点的方法以实现对火电最优运行点的寻找;利用二分法,结合后校验的思路,通过迭代收敛找出最优的煤耗微增率点;根据得到的煤耗微增率反解出火电机组的具体出力,达到实现最优调度的目的。
进一步的,煤耗微增率目标IRtarget是诱导实现最优火电机组出力的手段;煤耗微增率目标可以反推出机组出力的调整量,从而确定出机组的出力;另外,所以机组的出力调整后的煤耗微增率最终是需要对煤耗微增率目标进行“趋同”,而等煤耗微增率是机组出力最优运行点,因此设置煤耗微增率目标IRtarget可以通过寻找该目标,从而寻找到最优运行点。
进一步的,火电出力变量之和ΔPt g,sum是煤耗微增率目标搜寻方向的盘踞,煤耗微增率目标过大造成系统出力大于负荷;反之过小造成系统出力小于负荷,无法实现功率的平衡。
进一步的,煤耗微增率是煤耗曲线的一阶导数;对于二次函数形式的煤耗曲线,通过计算出的煤耗微增率目标,可以知道火电机组在向煤耗微增率目标IRtarget调整需要的出力变化量根据能够计算出此时此刻火电的机组出力。
进一步的,负荷改变量ΔPd,sum是负荷的变化量,火电机组出力需要根据负荷的变化量对其机组的进行调整。
进一步的,判断此时的出力调整量和负荷该变量是否满足精度要求,不满足则需要重新进行煤耗微增率目标的调整,以实现机组出力调整量等于负荷改变量的结果。
进一步的,修正煤耗微增率目标,是迭代过程中的行为。能够逐渐修正找到满足出力调整量等于负荷变化量的最优火电机组运行的煤耗微增率目标
进一步的,修正后煤耗微增率目标,用于校验修正后煤耗微增率目标是否满足调整精度;满足条件,修正后的煤耗微增率目标作为最终的煤耗微增率目标知道火电机组出力调整,不满足条件则作为新的迭代边界
进一步的,时刻t火电机组i的出力是最终计算得到的火电机组出力的结果,是目标结果。
可以理解的是,上述第二方面的有益效果可以参见上述第一方面中的相关描述,在此不再赘述。
综上所述,本发明结合煤耗微增率建立基于二分法的火电经济调度方法,能够快速计算出火电经济调度的出力安排,找出一种较为经济的电力系统运行方式。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为二分法收敛速度图;
图2为周电力电量平衡图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要理解的是,术语“包括”和“包含”指示所描述特征、整体、步骤、操作、元素和/或组件的存在,但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。
还应当理解,在本发明说明书中所使用的术语仅仅是出于描述特定实施例的目的而并不意在限制本发明。如在本发明说明书和所附权利要求书中所使用的那样,除非上下文清楚地指明其它情况,否则单数形式的“一”、“一个”及“该”意在包括复数形式。
还应当进一步理解,在本发明说明书和所附权利要求书中使用的术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合,并且包括这些组合,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。另外,本文中字符“/”,一般表示前后关联对象是一种“或”的关系。
应当理解,尽管在本发明实施例中可能采用术语第一、第二、第三等来描述预设范围等,但这些预设范围不应限于这些术语。这些术语仅用来将预设范围彼此区分开。例如,在不脱离本发明实施例范围的情况下,第一预设范围也可以被称为第二预设范围,类似地,第二预设范围也可以被称为第一预设范围。
取决于语境,如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”或“响应于检测”。类似地,取决于语境,短语“如果确定”或“如果检测(陈述的条件或事件)”可以被解释成为“当确定时”或“响应于确定”或“当检测(陈述的条件或事件)时”或“响应于检测(陈述的条件或事件)”。
在附图中示出了根据本发明公开实施例的各种结构示意图。这些图并非是按比例绘制的,其中为了清楚表达的目的,放大了某些细节,并且可能省略了某些细节。图中所示出的各种区域、层的形状及它们之间的相对大小、位置关系仅是示例性的,实际中可能由于制造公差或技术限制而有所偏差,并且本领域技术人员根据实际所需可以另外设计具有不同形状、大小、相对位置的区域/层。
输入数据涵盖如下内容:
①净负荷曲线(特指:已经完成新能源带负荷曲线和库容式水电带负荷曲线的净负荷);
②火电机组基本参数(包括:火电机组最大(最小)技术出力、火电机组最大上(下)爬坡速率;火电机组煤耗曲线);
输出数据涵盖如下内容:
①火电机组出力;
②火电机组机组组合状态;
火电典型煤耗曲线:
其中,Fi g为火电机组i的煤耗曲线;为时刻t火电机组i的出力;为火电机组i煤耗曲线的二次项/一次项/常数项。
对火电典型煤耗曲线求微分即得对应火电机组的煤耗微增率:
此时,获得了针对本发明最重要的对象为火电典型煤耗曲线/火电典型煤耗微增率。
本发明一种基于二分法的火电集群调度方法,包括以下步骤:
S1、对煤耗微增率的上限和下限进行初始化,得到煤耗微增率上限IRupper和下限IRlower,对煤耗微增率上限IRupper和下限IRlower进行二分,得到煤耗微增率目标IRtarget;
煤耗微增率上限IRupper和下限IRlower具体为:
IRupper=1×107 (3)
IRlower=0×107 (4)
其中,IRupper为煤耗微增率上限;IRlower为煤耗微增率下限。
煤耗微增率目标IRtarget具体为:
其中,IRtarget为煤耗微增率目标。
S2、对火电机组出力变化量之和ΔPt g,sum和负荷改变量ΔPt d,sum初始化,得到火电出力对应变量之和ΔPt g,sum和负荷改变量ΔPd,sum,将步骤S1得到的煤耗微增率目标IRtarget进行计算,得到火电出力之和ΔPt g,sum;
火电出力对应变量之和ΔPt g,sum为:
ΔPt g,sum=0.0 (6)
负荷改变量ΔPd,sum为:
其中,ΔPt g,sum为时刻t火电出力之和;ΔPt d,sum为火电出力之和;Ωd为负荷集合;为时刻t负荷i的负荷需求;为时刻t-1负荷i的负荷需求。
火电出力之和ΔPt g,sum为:
按如下公式计算得到:
其中,Ωg为火电集合;为时刻t火电机组i按照煤耗微增率目标IRtarget进行调整的出力变化量;为时刻t-1火电机组i的煤耗微增率;为火电机组i的出力上限;为时刻t-1火电机组i的出力;RUi为火电机组i的上爬坡上限;Pi g为火电机组i的出力下限;RDi为火电机组i的下爬坡上限。
S3、对步骤S2得到的火电机组出力变化量之和ΔPt g,sum和负荷改变量ΔPt d,sum的大小进行判断,对煤耗微增率目标进行修正IRtarget;
对煤耗微增率目标进行修正IRtarget如下:
S4、利用步骤S3得到的火电出力变化量之和ΔPt g,sum对计算精度进行判断,当|ΔPt g,sum-ΔPt d,sum|<limit不成立,返回步骤S2变化量重新进行调整;若成立,则安排机组出力;
火电机组出力具体为:
其中,为时刻t火电机组i的出力。
S6、对步骤S5获得的火电机组新的出力情况进行传递并返回步骤S3,进行下一时刻的出力安排。
本发明再一个实施例中,提供一种基于二分法的火电集群调度计算系统,该系统能够用于实现上述基于二分法的火电集群调度方法,具体的,该基于二分法的火电集群调度计算系统包括二分模块、初始模块以及调度模块。
其中,二分模块,获取煤耗微增率上限IRupper和下限IRlower,经二分处理得到煤耗微增率目标IRtarget;
初始模块,初始化火电出力变量之和ΔPt g,sum和负荷改变量ΔPd,sum,利用二分模块得到的煤耗微增率目标IRtarget得到火电出力变量之和ΔPt g,sum;
调度模块,将初始模块得到的火电机组出力变化量之和ΔPt g,sum与负荷改变量ΔPt d,sum的大小进行判断,当|ΔPt g,sum-ΔPt d,sum|<limit时,安排火电机组出力,对获得的火电机组出力情况进行传递,返回初始模块,进行火电机组下一时刻的出力安排,实现火电集群调度。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中的描述和所示的本发明实施例的组件可以通过各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
首先理论上证明二分法能够保证收敛,收敛速度可以稳定保持为线性收敛效果,收敛效果如图1所示。
为验证计算方法的有效性,选取ROTS测试系统对所提方案进行测试。测试系统、计算平台和计算环境基本情况如下:
ROTS测试系统包含:
1)32处火电场站;
2)4处光伏场站;
3)5处风电场站;
4)3处库容式水电;
5)全年8736小时负荷曲线和资源曲线;
测试平台:
1)CPU:AMD Ryzen 7 5700G with Radeon Graphics 3.8GHz;
2)RAM:16.0GB(13.9GB可用);
计算环境:
1)操作系统:Windows 11专业版22H2;
2)集成开发环境:Visual studio 2019;
3)CPLEX 12.80(用于进行对比)。
对比传统求解线性规划和所提方法求解效率进行对比,对比结果如表1所示。
表1 算法效率对照
请参阅图2,将本发明方法应用到ROTS测试系统汇总(以8736计算时间尺度为准,截取部分运行结果),获得如下结果(注:如无特殊说明,本发明中涉及电力功率的单位统一为100MW)。
火电机组前24小时出力情况(电力功率单位:100MW)
上述表格是ROTS测试系统全部火电机组24h的出力结果。
综上所述,本发明一种基于二分法的火电集群调度方法及系统,在基本不影响计算精度的情况下,能够实现对火电集群调度结果的快速计算,迅速获得调度方案。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。
Claims (10)
1.一种基于二分法的火电集群调度方法,其特征在于,获取煤耗微增率上限IRupper和下限IRlower,经二分处理得到煤耗微增率目标IRtarget;初始化火电出力变量之和ΔPt g,sum和负荷改变量ΔPd,sum,利用得到的煤耗微增率目标IRtarget得到火电出力变量之和ΔPt g,sum;将火电机组出力变化量之和ΔPt g,sum与负荷改变量ΔPt d,sum的大小进行判断,当|ΔPt g,sum-ΔPt d ,sum|<limit时,limit为设定阈值,安排火电机组出力,对获得的火电机组出力情况进行传递并返回初始化,进行火电机组下一时刻的出力安排,实现火电集群调度。
2.根据权利要求1所述的基于二分法的火电集群调度方法,其特征在于,煤耗微增率目标IRtarget为:
其中,IRtarget为煤耗微增率目标。
3.根据权利要求1所述的基于二分法的火电集群调度方法,其特征在于,火电出力变量之和ΔPt g,sum为:
其中,Ωg为火电集合,为时刻t火电机组i按照煤耗微增率目标IRtarget进行调整的出力变化量。
4.根据权利要求3所述的基于二分法的火电集群调度方法,其特征在于,时刻t火电机组i按照煤耗微增率目标IRtarget进行调整的出力变化量为:
其中,为时刻t-1火电机组i的煤耗微增率;为火电机组i的出力上限;为时刻t-1火电机组i的出力;RUi为火电机组i的上爬坡上限;P i g为火电机组i的出力下限;RDi为火电机组i的下爬坡上限。
5.根据权利要求1所述的基于二分法的火电集群调度方法,其特征在于,负荷改变量ΔPd,sum为:
其中,Ωd为负荷集合;为时刻t负荷i的负荷需求;为时刻t-1负荷i的负荷需求。
6.根据权利要求1所述的基于二分法的火电集群调度方法,其特征在于,当|ΔPt g,sum-ΔPt d,sum|≥limit时,对煤耗微增率目标IRtarget进行修正,根据修正后的煤耗微增率目标IRtarget重新计算火电出力变量之和ΔPt g,sum,返回初始化火电出力变量之和ΔPt g,sum和负荷改变量ΔPd,sum重新进行调整。
7.根据权利要求6所述的基于二分法的火电集群调度方法,其特征在于,对煤耗微增率目标进行修正IRtarget具体为:
其中,IRupper为设定的煤耗微增率上限,IRlower为设定的煤耗微增率下限。
8.根据权利要求6所述的基于二分法的火电集群调度方法,其特征在于,修正后的煤耗微增率目标IRtarget为:
9.根据权利要求1所述的基于二分法的火电集群调度方法,其特征在于,t时刻火电机组i的出力安排为:
其中,为时刻t-1火电机组i的出力;为时刻t-1火电机组i的煤耗微增率;为火电机组i的出力上限;RUi为火电机组i的上爬坡上限;P i g为火电机组i的出力下限;RDi为火电机组i的下爬坡上限。
10.一种基于二分法的火电集群调度计算系统,其特征在于,包括:
二分模块,获取煤耗微增率上限IRupper和下限IRlower,经二分处理得到煤耗微增率目标IRtarget;
初始模块,初始化火电出力变量之和ΔPt g,sum和负荷改变量ΔPd,sum,利用二分模块得到的煤耗微增率目标IRtarget得到火电出力变量之和ΔPt g,sum;
调度模块,将初始模块得到的火电机组出力变化量之和ΔPt g,sum与负荷改变量ΔPt d,sum的大小进行判断,当|ΔPt g,sum-ΔPt d,sum|<limit时,安排火电机组出力,对获得的火电机组出力情况进行传递,返回初始模块,进行火电机组下一时刻的出力安排,实现火电集群调度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210907585.XA CN115173489B (zh) | 2022-07-29 | 2022-07-29 | 一种基于二分法的火电集群调度方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210907585.XA CN115173489B (zh) | 2022-07-29 | 2022-07-29 | 一种基于二分法的火电集群调度方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115173489A CN115173489A (zh) | 2022-10-11 |
CN115173489B true CN115173489B (zh) | 2023-10-31 |
Family
ID=83476743
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210907585.XA Active CN115173489B (zh) | 2022-07-29 | 2022-07-29 | 一种基于二分法的火电集群调度方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115173489B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116454890B (zh) * | 2023-04-20 | 2024-02-06 | 中国南方电网有限责任公司 | 基于scuc模型的机组组合控制方法、装置和设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103178552A (zh) * | 2013-02-07 | 2013-06-26 | 大连理工大学 | 一种动态机组组合与等微增率法相结合的火电机组节能负荷分配方法 |
CN103956779A (zh) * | 2014-05-19 | 2014-07-30 | 大唐七台河发电有限责任公司 | 一种基于能量平衡的火电机组全厂负荷经济调度技术 |
CN104993524A (zh) * | 2015-07-17 | 2015-10-21 | 三峡大学 | 一种基于改进离散粒子群算法的含风电电力系统动态调度方法 |
CN106998080A (zh) * | 2017-03-20 | 2017-08-01 | 国网浙江省电力公司 | 一种agc增量指令厂级节能优化分配方法 |
CN113937819A (zh) * | 2021-06-30 | 2022-01-14 | 国网甘肃省电力公司电力科学研究院 | 一种多能源短期优化调度方法 |
WO2022089304A1 (zh) * | 2020-10-27 | 2022-05-05 | 国家电网有限公司 | 自动发电控制方法和系统 |
-
2022
- 2022-07-29 CN CN202210907585.XA patent/CN115173489B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103178552A (zh) * | 2013-02-07 | 2013-06-26 | 大连理工大学 | 一种动态机组组合与等微增率法相结合的火电机组节能负荷分配方法 |
CN103956779A (zh) * | 2014-05-19 | 2014-07-30 | 大唐七台河发电有限责任公司 | 一种基于能量平衡的火电机组全厂负荷经济调度技术 |
CN104993524A (zh) * | 2015-07-17 | 2015-10-21 | 三峡大学 | 一种基于改进离散粒子群算法的含风电电力系统动态调度方法 |
CN106998080A (zh) * | 2017-03-20 | 2017-08-01 | 国网浙江省电力公司 | 一种agc增量指令厂级节能优化分配方法 |
WO2022089304A1 (zh) * | 2020-10-27 | 2022-05-05 | 国家电网有限公司 | 自动发电控制方法和系统 |
CN113937819A (zh) * | 2021-06-30 | 2022-01-14 | 国网甘肃省电力公司电力科学研究院 | 一种多能源短期优化调度方法 |
Non-Patent Citations (2)
Title |
---|
Study on Unit Optimal Scheduling Considering the joint constraint of "deep peak load regulation and coal consumption";Dan Song et al.;《2020 Chinese Automation Congress (CAC)》;第1909-1913页 * |
基于等综合煤耗微增率的火电机组节能发电调度算法;吕素等;《中国电机工程学报》;第32卷(第32期);第1-7、S1页 * |
Also Published As
Publication number | Publication date |
---|---|
CN115173489A (zh) | 2022-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106786807B (zh) | 一种基于模型预测控制的风电场有功功率控制方法 | |
CN106655280B (zh) | 一种基于电量控制的梯级水电短期调峰模型及求解方法 | |
CN109447393A (zh) | 一种电力系统经济负荷分配的改进的粒子群算法 | |
CN107317353B (zh) | 一种含分布式光伏发电配电网电压控制方法及系统 | |
CN105207253B (zh) | 考虑风电及频率不确定性的agc随机动态优化调度方法 | |
CN110400232B (zh) | 考虑电网断面约束的水电站群月度交易计划电量分解方法 | |
CN115173489B (zh) | 一种基于二分法的火电集群调度方法及系统 | |
CN113537562B (zh) | 一种新能源场站储能配置计算方法及系统 | |
CN109193814A (zh) | 一种新能源有功功率分配方法、装置和计算机设备 | |
CN112803494A (zh) | 一种含风光水火的多目标agc协调优化方法及系统 | |
CN108154266B (zh) | 一种含风电电力系统经济调度方法 | |
CN115618540A (zh) | 基于三层次动态变异率的风力发电机组优化布局方法 | |
CN117540767A (zh) | 一种水风光蓄短期多目标调度模型求解的pso-dp耦合嵌套算法 | |
CN109711605B (zh) | 一种多类型电源并网容量确定方法和装置 | |
CN104269886B (zh) | 用于风电并网系统容量效益裕度计算及分配的方法 | |
CN111798060A (zh) | 一种基于机组爬坡率估计的功率指令优化分配方法 | |
CN115663849B (zh) | 一种水电厂与工业温控负荷协同快速调控的方法及控制器 | |
CN117094746A (zh) | 一种虚拟发电厂优化方法、装置、设备及介质 | |
CN105511270A (zh) | 一种基于协同进化的pid控制器参数优化方法和系统 | |
CN115203861A (zh) | 多级机站通风系统风机变频优化方法、装置及设备 | |
CN113595153B (zh) | 一种新能源集群的输出功率优化方法及装置 | |
CN113036769B (zh) | 一种电力系统静态电压稳定分析方法及系统 | |
CN114944659A (zh) | 一种风光储联合电站自动控制方法及系统 | |
CN110717694B (zh) | 基于新能源消纳期望值的储能配置随机决策方法及装置 | |
CN108879764B (zh) | 一种抑制风电光伏不确定性鲁棒优化方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |