CN115173488B

CN115173488B - 基于驻点启停策略的电力系统机组组合可行解计算方法

Info

Publication number: CN115173488B
Application number: CN202210906422.XA
Authority: CN
Inventors: 吴晨; 谈健; 陈琛; 薛贵元; 牛文娟; 吴垠
Original assignee: Economic and Technological Research Institute of State Grid Jiangsu Electric Power Co Ltd
Current assignee: Economic and Technological Research Institute of State Grid Jiangsu Electric Power Co Ltd
Priority date: 2022-07-29
Filing date: 2022-07-29
Publication date: 2023-11-24
Anticipated expiration: 2042-07-29
Also published as: CN115173488A

Abstract

本发明公开了一种基于驻点启停策略的电力系统机组组合可行解计算方法，包括：收集待计算电力系统基础数据；将机组分为A、B类；建立电力系统机组组合驻点简化模型；采用混合整数规划求解算法求解机组组合驻点简化模型中四个驻点处A类机组启停状态值，求解过程中B类机组运行状态值为运行日初始值；依据机组组合的运行特点，辨识出A类机组中必开必停机组，形成机组组合整数变量削减模型；采用混合整数规划求解算法求解机组组合整数变量削减模型，获取运行日任意时段机组组合结果；本发明依据驻点时段机组启停值，确定了部分离散变量取值，降低了问题规模，能够快速计算出满足系统安全稳定经济运行需求的近优解乃至最优解。

Description

基于驻点启停策略的电力系统机组组合可行解计算方法

技术领域

本发明属于电气工程领域，尤其涉及基于驻点启停策略的电力系统机组组合可行解快速计算方法。

背景技术

在机组组合问题中，每台火电机组都是一个独立的建模对象，其运行状态采用二进制整数变量描述。目前，我国省级电网的机组数目在数十至百台范围内，若干个省级电网组成的机组规模可超过三百台。若对电力系统的全部机组进行逐一建模，模型将包含庞大数量的二进制整数变量和约束，求解问题变得异常的庞大和复杂，很难求解出全局最优的解。

快速计算可行解与耗时获取最优解的权衡，实际应用中更偏向于快速计算可行解。为此，一些文献对模型中约束条件进行了辨识，目的是降低模型复杂度。工程上，某些情况会依据调度人员多年工程经验，提前筛选出起作用的约束，以减少计算过程中的变量，但是人工经验的执行手段简单且缺乏全面量化的依据。

发电机组的启停状态(0-1离散变量)是电力系统机组组合模型中重要约束条件，它既制约目标函数又影响模型求解速度。由于锅炉和汽机本身的技术条件要求，机组一旦启动，必须运行一段时间而不能马上关闭。同理，机组停机后，也应停运一段时间后才能运行。实际运行中每一天启停变化的机组数目很少，大多数机组处于必开或者必停状态。因此，根据机组运行特点，辨识出必开或者必停机组，大幅减少全时段机组组合模型中的整数变量数目是降低电力系统机组组合问题求解难度的有效途径。

发明内容

为了解决现有技术中存在的问题，本发明提供一种基于驻点启停策略的电力系统机组组合可行解快速计算方法，以解决现有庞大和复杂电力系统机组组合问题求解速度慢的问题，具体由以下技术方案实现：

所述基于驻点启停策略的电力系统机组组合可行解计算方法，该方法包括：

步骤1)收集待计算电力系统基础数据；

步骤2)依据运行日前一天机组启停变化将机组分为A、B类，所述A类机组是指运行日前一天没有发生启停变化的机组，B类机组是指运行日前一天发生过由开到停或者由停到开变化的机组；

步骤3)建立目标函数公式，增设包含：启动成本约束、系统功率平衡约束、机组出力约束、机组爬坡速率约束、旋转备用约束、机组启停约束的约束条件，形成电力系统机组组合驻点简化模型；

步骤4)采用混合整数规划求解算法求解所述电力系统机组组合驻点简化模型中四个驻点A类机组启停状态值，求解过程中B类机组运行状态值为运行日初始值，所述四个驻点分别对应电力系统净负荷最大时段、电力系统净负荷最小时段、电力系统净负荷向上爬坡最快时段、电力系统净负荷向下爬坡最快时段；

步骤5)依据机组组合的运行特点，辨识出A类机组中必开必停类机组；

步骤6)依据所述必开必停机组信息，在电力系统机组组合驻点简化模型基础上增设A类机组必开必停类机组启停状态约束条件，建立电力系统机组组合整数变量削减模型；

步骤7)采用混合整数规划求解算法快速求解电力系统机组组合整数变量削减模型，获取运行日任意时段机组组合可行解。

所述基于驻点启停策略的电力系统机组组合可行解计算方法的进一步设计在于，所述步骤2)中所述A类机组数目和B类机组数目满足式(1)：

N_A+N_B＝N_T (1)

式中，N_A为A类机组数目，N_B为B类机组数目,N_T为机组总数；

设定A类机组所组成的机组集合为J₁，B类机组所组成的机组集合J₂，A类机组t时段所对应的机组启停数据集为A_i∈J₁，B类机组t时段所对应的机组启停数据集为B_i∈J₂。

所述基于驻点启停策略的电力系统机组组合可行解计算方法的进一步设计在于，所述步骤3)中所述电力系统机组组合驻点简化模型中目标函数如式(2)：

式中，F_c为总发电成本；为机组i在第t时段的发电成本；/>为机组i在第t时段的启动成本；K为时段集合；J为机组集合。

所述基于驻点启停策略的电力系统机组组合可行解计算方法的进一步设计在于，所述步骤3)中采用分段线性函数来近似表示如式(3)所示的目标函数的发电成本，

f_i ^min＝c_i P_i ²+b_i P_i +a_i (3)

式(3)中，c_i，b_i，a_i分别为机组发电成本二次项系数、一次及常数项；为机组i在第t时段的状态(“1”表示机组处于运行状态；“0”表示机组处于停机状态)；P_i 为机组i最小技术出力；m_i为线性煤耗系数；

电力系统机组组合驻点简化模型目标函数中启动成本约束计算公式如式(4)：

式(4)中，分别为机组i的热起动和冷起动费用，t^off(i,t-1)为机组i在第t-1个时段连续停机时段数，/>为机组i最小连续停机时间，/>为机组i冷启动时间；

电力系统机组组合驻点简化模型中系统功率平衡约束计算公式如式(5)：

式(5)中，为新能源机组w第t时段的出力，/>为系统在t时段的总负荷需求，W为新能源机组集合，P_i ^t为机组i在t时段的出力；

电力系统机组组合驻点简化模型中机组出力约束计算公式如式(6)：

式(6)中，为机组i最大技术出力，

电力系统机组组合驻点简化模型中机组爬坡速率约束计算公式如式(7)：

式(7)中，P_i ^up、P_i ^down分别为机组i上、下爬坡功率变化量限值；

电力系统机组组合驻点简化模型中旋转备用约束计算如式(8)：

式(8)中，为第t时段的系统总备用功率；

电力系统机组组合驻点简化模型中机组启停约束计算公式如式(9)：

式(9)中，为机组i最小连续运行时间，/>为机组i最小连续停机时间，T^on(i,t)为机组i到第t时段连续运行的时段数。

所述基于驻点启停策略的电力系统机组组合可行解计算方法的进一步设计在于，所述步骤4)中设定电力系统净负荷功率计算式如式(10)，

式(10)中，为第t时刻系统净负荷功率；

相邻时间段净负荷变化量为：当/>则t_Δt1＝k+1；当/>则t_Δt2＝m+1；/>T为整个计算周期的时段数；

求解四个驻点A类机组启停状态值的具体步骤为：

步骤4-1)首先将4个驻点对应时段按照从小到大的顺序进行排列，并且依次用k₁、k₂、k₃、k₄表示，即t_Pmax表示为电力系统净负荷最大时段，t_Pmin表示为电力系统净负荷最小时段，t_Δt1表示为电力系统净负荷向上爬坡最快时段为运行日相邻时间段净负荷变化量最大所对应时段，t_Δt2表示为电力系统净负荷向下爬坡最快时段为运行日相邻时间段净负荷变化量最小所对应时段；

步骤4-2)根据式(7)所示的爬坡约束计算公式，A类机组四个驻点按照如下关系式确定：

步骤4-3)根据式(9)所示的启停约束计算公式，A类机组四个驻点T^on(A_i,t-1)、T^off(A_i，t-1)、按照如下关系式确定：

T^on(A_i,k₂-1)＝T^on(A_i,k₁)

T^on(A_i,k₃-1)＝T^on(A_i,k₂)

T^on(A_i,k₄-1)＝T^on(A_i,k₃)

T^off(A_i,k₂-1)＝T^off(A_i,k₁)

T^off(A_i,k₃-1)＝T^off(A_i,k₂)

T^off(A_i,k₄-1)＝T^off(A_i,k₃)

式中，A_i∈J₁,

步骤4-4)采用混合整数规划求解算法求解驻点简化模型中四个驻点，求解结果表示为：

所述基于驻点启停策略的电力系统机组组合可行解计算方法的进一步设计在于，求解四个驻点A类机组启停状态值时B类机组运行状态值为运行日初始值即

所述基于驻点启停策略的电力系统机组组合可行解计算方法的进一步设计在于，根据式(10)建立A类机组中必开必停类机组启停状态约束计算公式，

所述基于驻点启停策略的电力系统机组组合可行解计算方法的进一步设计在于，所述步骤1)中电力系统基础数据包括：机组总数、各机组启动成本、各机组发电成本系数、运行日各时段系统总负荷、运行日各时段新能源总出力、各机组最小技术出力、各机组最大技术出力、各时段系统总备用功率、各机组最小连续运行时间、各机组最小连续停机时间、各机组上、下爬坡功率变化量限值、各机组初始出力、各机组运行日前一天启停变化情况以及各机组运行日初始启停状。

本发明还涉及一种电子设备，包括存储器、处理器以及计算机程序，其中所述计算机程序存储于所述存储器中，并被配置为由所述处理器执行以实现如所述基于驻点启停策略的电力系统机组组合可行解计算方法。

本发明还涉及一种计算机可读存储介质，所述存储介质存储有计算机程序，所述计算机程序被处理器执行以实现所述基于驻点启停策略的电力系统机组组合可行解计算方法。

本发明的有益效果：

本发明的基于驻点启停策略的电力系统机组组合可行解快速计算方法利用大型机组不能频繁启动，启停机都是24小时以上的运行特点，确定大部分离散变量取值，大幅减少模型中引入的整数变量个数，显著提升了模型求解效率；

通过计算4个驻点可行解时B类机组启停值为已知量，整数变量削减模型中B类机组启停值为未知量，充分体现计算过程考虑了每一个机组启停要求约束，获取的可行解绝大多数是最优的，结果贴近实际的日前发电计划，保障了计算结果的有效性。

通过挖掘机组运行特点，科学地辨识出必开必停类机组，缩小了模型求解的可行域，得到更加紧凑的求解空间，模型整数变量数减少方法计算过程清晰，具有较高的适用性和推广性。

附图说明

图1为本发明基于驻点启停策略的电力系统机组组合可行解计算方法的流程示意图。

图2为10机24节点算例拓扑图。

图3为运行日系统总负荷曲线示意图。

图4为运行日新能源出力曲线示意图。

图5为运行日系统净负荷曲线示意图。

具体实施方式

下面结合附图对本发明实施例中的技术方案进行清楚、完整地描述。

如图1，本发明的基于驻点启停策略的电力系统机组组合可行解快速计算方法包括如下步骤：

步骤1)收集待计算电力系统基础数据。电力系统基础数据包括：机组总数、各机组启动成本、各机组发电成本系数、运行日各时段系统总负荷、运行日各时段新能源总出力、各机组最小技术出力、各机组最大技术出力、各时段系统总备用功率、各机组最小连续运行时间、各机组最小连续停机时间、各机组上、下爬坡功率变化量限值、各机组初始出力、各机组运行日前一天启停变化情况、各机组运行日初始启停状态。

步骤2)依据运行日前一天机组启停变化将机组分为A、B类。运行日前一天机组启停变化可根据步骤1)中基础数据获取。其中，A类机组是指运行日前一天没有发生启停变化的机组，B类机组是指运行日前一天发生过由开到停或者由停到开变化的机组。

A类机组数目和B类机组数目满足如下式子：

N_A+N_B＝N_T

式中，N_A为A类机组数目，N_B为B类机组数目，N_T为机组总数。

A类机组所组成的机组集合为J₁，B类机组所组成的机组集合J₂。

A类机组i在t时段所对应的机组启停决策变量集为A_i∈J₁.

B类机组i在t时段所对应的机组启停决策变量集为B_i∈J₂。

步骤3)建立电力系统机组组合驻点简化模型。电力系统机组组合驻点简化模型中目标函数计算公式X₁为：

式中，F_c为总发电成本；C^P(i,t)为机组i在第t时段的发电成本；C^U(i,t)为机组i在第t时段的启停成本；K为时段集合；J为机组集合。

目标函数中发电成本可采用分段线性函数来近似表示，计算公式X₂为：

f_i ^min＝c_i P_i ²+b_i P_i +a_i

式中，c_i，b_i，a_i分别为机组i发电成本二次项系数、一次及常数项；为机组i在第t时段的状态(“1”表示机组处于运行状态；“0”表示机组处于停机状态)；P_i 为机组i最小技术出力；m_i为线性煤耗系数。

电力系统机组组合驻点简化模型目标函数中启动成本计算公式X₃为：

式中，分别为火电机组i的热起动和冷起动费用，t^off(i,t-1)为机组i在第t-1个时段连续停机的时段数，T_i ^off 为机组i最小连续停机时间，/>为机组i冷启动时间。

电力系统机组组合驻点简化模型中系统功率平衡约束计算公式X₄为：

式中，为新能源机组w第t时段的出力，/>为系统在t时段的总负荷需求，W为新能源机组集合，P_i ^t为机组i在t时段的出力。

电力系统机组组合驻点简化模型中机组出力约束计算公式X₅为：

式中，为机组i最大技术出力。

电力系统机组组合驻点简化模型中爬坡速率约束计算公式X₆为：

式中，P_i ^up、P_i ^down分别为机组i上、下爬坡功率变化量限值。

电力系统机组组合驻点简化模型中旋转备用约束计算公式X₇为：

式中，为第t时段的系统总备用功率。

电力系统机组组合驻点简化模型中机组启停约束计算公式X₈为:

式中，为机组i最小连续运行时间，/>为机组i最小连续停机时间，T^on(i,t-1)为机组i到第t-1时段连续运行时段数。

步骤4)采用混合整数规划求解算法求解机组组合驻点简化模型中4个驻点A类机组启停状态值，求解过程中B类机组运行状态值为运行日初始值。

电力系统净负荷功率计算式为：

式中，为第t时刻系统净负荷功率。

四个驻点为电力系统净负荷最大时段、电力系统净负荷最小时段、电力系统净负荷向上爬坡最快时段、电力系统净负荷向下爬坡最快时段。电力系统净负荷最大时段用t_Pmax表示，电力系统净负荷最小时段用t_Pmin表示，电力系统净负荷向上爬坡最快时段为运行日相邻时间段净负荷向上变化量最大的时段用t_Δt1表示，电力系统净负荷向下爬坡最快时段为运行日相邻时间段净负荷向下变化量最大的时段用t_Δt2表示。

则相邻时间段净负荷变化量为:当/>则t_Δt1＝k+1；当/>则t_Δt2＝m+1；/>式中，T为整个计算周期的时段数。

求解四个驻点A类机组启停状态值时B类机组运行状态值为运行日初始值即

接着，采用混合整数规划求解算法求解机组组合驻点简化模型中A类机组四个驻点可行解时，求解过程如下：

步骤4-1)首先将4个驻点对应时段按照从小到大的顺序进行排列，并且依次用k₁、k₂、k₃、k₄表示即步骤4-2)爬坡约束计算公式中A类机组四个驻点前一时刻出力值按照如下方法确定。

步骤4-3)启停约束计算公式中A类机组四个驻点T^on(A_i,t-1)、T^off(A_i,t-1)、按照如下方法确定：

T^on(A_i,k₂-1)＝T^on(A_i,k₁)

T^on(A_i,k₃-1)＝T^on(A_i,k₂)

T^on(A_i,k₄-1)＝T^on(A_i,k₃)

T^off(A_i,k₂-1)＝T^off(A_i,k₁)

T^off(A_i,k₃-1)＝T^off(A_i,k₂)

T^off(A_i,k₄-1)＝T^off(A_i,k₃)

式中，A_i∈J₁。

步骤4-4)采用混合整数规划求解算法求解驻点简化模型中4个驻点。求解结果表示为：

步骤5)依据机组组合的运行特点，辨识出A类机组中必开必停类机组。

必开必停类机组为电网运行中必开发电机组或者必停发电机组，不需要在机组模型中进行机组状态的优化决策。由于锅炉和汽机本身的技术条件要求，大型机组不能频繁启动，启停机都是24小时以上，所以不会在一天内出现两次及以上的启停，最多出现一次启停。运用机组该特点，将部分机组启停值辨识出来，降低模型中变量数。

运行日A类机组初始时刻状态可根据电力系统基础数据获取，为数据集如果/>则A类机组中第A_i台为运行日内必开机组即为如果/>则A类机组中第A_i台机组为运行日内必停机组，即为/>

步骤6)依据辨识得到的必开必停机组信息，建立电力系统机组组合整数变量削减模型。

电力系统机组组合整数变量削减模型为电力系统机组组合驻点简化模型基础上增加A类机组必开必停类机组启停状态约束条件。

A类机组中必开必停类机组启停状态约束计算公式X₉为：

电力系统机组组合整数变量削减模型包含计算公式X₁-X₉。

图2为所选取的10机24节点测试系统的系统图，母线6、14处各有1个风电场，以此为例进行实例说明。

步骤1)收集待计算电力系统基础数据。

如图3所示，机组集中机组总数N_T＝10，运行日24时段系统总负荷曲线如图3所示、运行日24时段新能源出力曲线如图4所示，各机组运行参数如表1所示。

表1各机组运行参数

步骤2)依据运行日前一天机组启停变化将机组分为A、B类

根据表1中数据运行日前一天10机组中1、3、5机组启停发生过变化，归为B类机组，其他机组归为A类机组。

因此，A类机组数目为N_A＝7，B类机组数目为N_B＝3

N_A+N_B＝10

A类机组集合为

J₁＝{机组2、机组4、机组6、机组7、机组8、机组9、机组10}

B类机组集合为J₂＝{机组1、机组3、机组5}

步骤3)建立电力系统机组组合驻点简化模型

电力系统机组组合驻点简化模型中目标函数计算公式X₁为：

目标函数中发电成本可采用分段线性函数来近似表示,计算公式X₂为：

f_i ^min＝c_i P_i ²+b_i P_i +a_i

电力系统机组组合驻点简化模型中目标函数中启动成本计算公式X₃为：

步骤4)采用混合整数规划求解算法求解初始模型中4个驻点A类机组启停状态值，求解过程中B类机组运行状态值为运行日初始值。

根据图3、图4的数据及系统净负荷变化量公式得到运行日24时段系统净负荷大小，如图5所示。

系统净负荷最大时段系统净负荷最小时段/>20时段系统净负荷变化量最大，23时段系统净负荷变化量最小，所以t_Δt1＝20，t_Δt2＝23

求解四个驻点过程中1、3、5机组为B类机组，其运行状态值为运行日初始值，可根据表1中数据获取，如表2所示1、3、5机组运行状态值。

表2 B类机组运行状态值

采用混合整数规划求解算法求解初始模型中A类机组4个驻点可行解时，求解过程如下：

步骤4-1)4个离散驻点t_Pmin＝1、t_Pmax＝10、t_Δt1＝20、t_Δt2＝23按照从小到大的顺序进行排列为t_Pmin＜t_Pmax＜t_Δt1＜t_Δt2，依次用k₁、k₂、k₃、k₄表示即k₁＝1、k₂＝10、k₃＝20、k₄＝23

步骤4-2)爬坡约束计算公式X6中A类机组四个驻点按照如下关系式确定,

步骤4-3)启停约束计算公式X8中A类机组四个驻点T^on(A_i,t-1)、T^off(A_i，t-1)、按照如下方法确定

T^on(A_i,10-1)＝T^on(A_i,1)

T^on(A_i,20-1)＝T^on(A_i,10)

T^on(A_i,23-1)＝T^on(A_i,20)

T^off(A_i,10-1)＝T^off(A_i,1)

T^off(A_i,20-1)＝T^off(A_i,10)

T^off(A_i,23-1)＝T^off(A_i,20)

步骤4-4)采用混合整数规划求解算法求解驻点简化模型中4个驻点求解结果可表示为：

步骤5)依据机组组合的运行特点，辨识出A类机组中必开必停类机组

运行日A类机组初始时段状态启停值可根据表1中数据获取如果/>则A类机组中第A_i台为运行日必开机组即为/>如果/>则A类机组中第A_i台机组为运行日必停组即为/>

步骤6)依据辨识得到的必开必停机组信息，建立电力系统机组组合整数变量削减的模型。

电力系统机组组合整数变量削减模型为电力系统机组组合驻点简化模型基础上增加A类机组必开必停类机组启停状态约束条件

A类机组中必开必停类机组启停状态约束计算公式X₉为：

电力系统机组组合整数变量削减模型包含计算公式X₁-X₉。

本发明方案所公开的技术手段不仅限于上述实施方式所公开的技术手段，还包括由以上技术特征任意组合所组成的技术方案。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。

Claims

1.一种基于驻点启停策略的电力系统机组组合可行解计算方法，其特征在于，该方法包括：步骤1)收集待计算电力系统基础数据；

2.根据权利要求1所述的基于驻点启停策略的电力系统机组组合可行解计算方法，其特征在于，所述步骤2)中所述A类机组数目和B类机组数目满足式(1)：

N_A+N_B＝N_T (1)

式中，N_A为A类机组数目，N_B为B类机组数目,N_T为机组总数；

3.根据权利要求1所述的基于驻点启停策略的电力系统机组组合可行解计算方法，其特征在于，所述步骤3)中所述电力系统机组组合驻点简化模型中目标函数如式(2)：

4.根据权利要求1所述的基于驻点启停策略的电力系统机组组合可行解计算方法，其特征在于，所述步骤3)中采用分段线性函数来近似表示如式(3)所示的目标函数的发电成本，

式(4)中，分别为机组i的热起动和冷起动费用，t^off(i,t-1)为机组i在第t-1个时段连续停机时段数，T_i ^o ^ff 为机组i最小连续停机时间，/>为机组i冷启动时间；

式(6)中，为机组i最大技术出力，

式(8)中，为第t时段的系统总备用功率；

5.根据权利要求1所述的基于驻点启停策略的电力系统机组组合可行解计算方法，其特征在于，所述步骤4)中设定电力系统净负荷功率计算式如式(10)，

式(10)中，为第t时刻系统净负荷功率；

相邻时间段净负荷变化量为：当/>则t_Δt1＝k+1；当/>则t_Δt2＝m+1；/>k、m∈0、2...T-1，T为整个计算周期的时段数；

求解四个驻点A类机组启停状态值的具体步骤为：

步骤4-1)首先将4个驻点对应时段按照从小到大的顺序进行排列，并且依次用k₁、k₂、k₃、k₄表示，即k₁、k₂、k₃、k₄∈{t_Pmax、t_Pmin、t_Δt1、t_Δt2}，t_Pmax表示为电力系统净负荷最大时段，t_Pmin表示为电力系统净负荷最小时段，t_Δt1表示为电力系统净负荷向上爬坡最快时段为运行日相邻时间段净负荷变化量最大所对应时段，t_Δt2表示为电力系统净负荷向下爬坡最快时段为运行日相邻时间段净负荷变化量最小所对应时段；

步骤4-3)根据式(9)所示的启停约束计算公式，A类机组四个驻点T^on(A_i,t-1)、T^off(A_i,t-1)、按照如下关系式确定：

T^on(A_i,k₂-1)＝T^on(A_i,k₁)

T^on(A_i,k₃-1)＝T^on(A_i,k₂)

T^on(A_i,k₄-1)＝T^on(A_i,k₃)

T^off(A_i,k₂-1)＝T^off(A_i,k₁)

T^off(A_i,k₃-1)＝T^off(A_i,k₂)

T^off(A_i,k₄-1)＝T^off(A_i,k₃)

式中，A_i∈J₁,

6.根据权利要求5所述的基于驻点启停策略的电力系统机组组合可行解计算方法，其特征在于，求解四个驻点A类机组启停状态值时B类机组运行状态值为运行日初始值即

7.根据权利要求5所述的基于驻点启停策略的电力系统机组组合可行解计算方法，其特征在于，根据式(10)建立A类机组中必开必停类机组启停状态约束计算公式，

8.根据权利要求1所述的基于驻点启停策略的电力系统机组组合可行解计算方法，其特征在于，所述步骤1)中电力系统基础数据包括：机组总数、各机组启动成本、各机组发电成本系数、运行日各时段系统总负荷、运行日各时段新能源总出力、各机组最小技术出力、各机组最大技术出力、各时段系统总备用功率、各机组最小连续运行时间、各机组最小连续停机时间、各机组上、下爬坡功率变化量限值、各机组初始出力、各机组运行日前一天启停变化情况以及各机组运行日初始启停状态值。

9.一种电子设备，其特征在于，包括存储器、处理器以及计算机程序，其中所述计算机程序存储于所述存储器中，并被配置为由所述处理器执行以实现如权利要求1-8任一项所述的基于驻点启停策略的电力系统机组组合可行解计算方法。

10.一种计算机可读存储介质，其特征在于所述存储介质存储有计算机程序，所述计算机程序被处理器执行以实现权利要求1-8任一项所述的基于驻点启停策略的电力系统机组组合可行解计算方法。