CN112258076B - 一种多时段高维投影机组组合模型的构建方法和装置 - Google Patents
一种多时段高维投影机组组合模型的构建方法和装置 Download PDFInfo
- Publication number
- CN112258076B CN112258076B CN202011208689.9A CN202011208689A CN112258076B CN 112258076 B CN112258076 B CN 112258076B CN 202011208689 A CN202011208689 A CN 202011208689A CN 112258076 B CN112258076 B CN 112258076B
- Authority
- CN
- China
- Prior art keywords
- unit
- constraint
- period
- state
- power
- 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
- 238000010276 construction Methods 0.000 title abstract description 15
- 230000006870 function Effects 0.000 claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 30
- 230000009194 climbing Effects 0.000 claims abstract description 18
- 238000005516 engineering process Methods 0.000 claims abstract description 11
- 230000007704 transition Effects 0.000 claims description 10
- 230000001174 ascending effect Effects 0.000 claims description 5
- 229910052799 carbon Inorganic materials 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 3
- 230000000750 progressive effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 description 18
- 230000006872 improvement Effects 0.000 description 8
- 238000010248 power generation Methods 0.000 description 8
- 230000008901 benefit Effects 0.000 description 6
- 238000004134 energy conservation Methods 0.000 description 6
- 239000000446 fuel Substances 0.000 description 6
- 230000009467 reduction Effects 0.000 description 6
- 230000009466 transformation Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 239000005431 greenhouse gas Substances 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
Classifications
-
- 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/067—Enterprise or organisation modelling
-
- 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
- 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (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
本发明适用于电力技术领域，提供了一种多时段高维投影机组组合模型的构建方法和装置，该方法通过采集电力系统电机机组的基础数据，利用变量投影技术将机组功率变量及与其相关的参数投影到区间[0,1]，并引入电机机组在多个时段内单个机组的所有运行状态，并基于机组的实际运行需求利用引入的状态变量系统地推导出多时段内的机组出力上界约束和爬坡约束的理想表达形式；利用引入的新状态变量与其他状态变量之间的关系消除新状态变量以提升模型的简洁性，加入其他约束和目标函数后最终得到多时段的机组组合模型。根据该建模方法，一个多时段高维投影机组组合模型被建立。该模型更紧凑，且计算效率更高。
Description
技术领域
本发明属于计算机技术领域，尤其涉及一种多时段高维投影机组组合模型的构建方法和装置。
背景技术
能源是现代社会的血液，是人类社会生存和发展的重要物质基础。然而，化石能源的大量消耗也导致能源短缺、环境污染等问题愈加严重。电能作为我国社会生产和生活的重要能源消费形式，在能源消费总量中占据的比重较大。然而，由于火力发电需要消耗大量的化石能源，并产生大量的二氧化碳等温室气体，使得电力行业成为了化石能源的主要“消费者”和温室气体的主要“生产者”。在此背景下，使用可再生能源替代化石能源进行发电和提高能源的利用效率是缓解能源紧张、改善环境质量的重要措施。然而，可再生能源受天气等不确定因素影响较大，提高其利用率仍有许多技术难点要解决，所以当下的最为有效的手段是提高化能源的利用效率。因此，优化火力发电机组的调度计划，对减少化石能源的使用和温室气体的排放具有十分重要的意义。
优化发电机组的调度计划，就是在保证电力系统安全运行的前提下，以总发电费用最少或者其它需求为目标优化发电机组的启停和出力计划。为了便于研究，人们根据数学规划理论将发电机组调度计划的优化问题转化为一类混合整数规划(mixed integerprogramming，MIP)问题，即机组组合(the unit commitment，UC)问题。然而，UC问题是一种具有非凸、高维和离散等特性的混合整数非线性规划(mixed integer nonlinearprogramming，MINLP)问题，求解十分困难。
因此，目前迫切需要一种计算效率高的UC模型，用于控制和调度电网发电资源，同时满足机组完成合同电量的实际需求，达到资源优化配置和节能减排的目标。
发明内容
本发明实施例的目的在于提供一种多时段高维投影机组组合模型的构建方法，旨在解决现有的UC模型求解困难，计算效率低，不能同时满足机组完成合同电量的实际需求，从而达到资源优化配置和节能减排的目标的问题。
本发明实施例是这样实现的，一种多时段高维投影机组组合模型的构建方法，所述方法包括：
采集电力系统发电机机组的基础数据，所述基础数据包括发电机机组的运行特性数据、负荷预测数据以及各时段的旋转备用数据；
利用变量投影技术将电机机组功率变量以及与其相关的参数进行投影，得到投影后的电机机组功率变量以及与其相关的参数，其中投影区间为[0,1]；
根据所述基础数据、投影后的电机机组功率变量以及与其相关的参数以及电机机组在运行、启动和关机三个时段内单个机组的所有运行状态变量，构建多时段高维投影机组组合模型。
本发明实施例的另一目的在于提供一种多时段高维投影机组组合模型的构建装置，所述装置包括：
数据采集单元，用于采集电力系统发电机机组的基础数据，所述基础数据包括发电机机组的运行特性数据、负荷预测数据以及各时段的旋转备用数据；
变量投影单元，用于利用变量投影技术将电机机组功率变量以及与其相关的参数进行投影，得到投影后的电机机组功率变量以及与其相关的参数，其中投影区间为[0,1]；以及模型构建单元，用于根据所述基础数据、投影后的电机机组功率变量以及与其相关的参数以及电机机组在运行、启动和关机三个时段内单个机组的所有运行状态变量，构建多时段高维投影机组组合模型。
本发明实施例提供的多时段高维投影机组组合模型的构建方法，通过将传统电力系统UC模型中的机组出力变量由半连续变量转化为连续变量，使UC模型中的出力约束式的数量至少减少了一半，降低了机组模型求解的复杂度；并且利用投影技术对目标函数进行变化，使目标函数值域范围更广；基于CPLEX12.7.1工具的模型求解表明，采用本发明的构建方法构建得到的三时段机组组合模型，在目标函数无分段线性化的情况下模型的计算效率至少提高了10％，在目标函数分段线性化的情况下，设置求解的精度为0.5％，模型的计算效率至少提高了20％。可见，本发明所构建得到的多时段高维投影机组组合模型的求解难度降低，且计算效率高，可用于控制和调度电网发电资源，能够同时满足机组完成合同电量的实际需求，以达到资源优化配置和节能减排的目标。
附图说明
图1为本发明实施例提供的一种多时段高维投影机组组合模型的构建方法的流程图；
图2为本发明实施例提供的一种多时段高维投影机组组合模型的构建装置的结构框图；
图3为本发明实施例提供的多时段高维投影机组组合模型的构建装置中的模型构建单元的结构框图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明实施例提供的多时段高维投影机组组合模型的构建方法，通过将传统电力系统UC模型中的机组出力变量由半连续变量转化为连续变量，使UC模型中的出力约束式的数量至少减少了一半，降低了机组模型求解的复杂度；并且利用投影技术对目标函数进行变化，使目标函数值域范围更广，且经验证，本发明构建的模型在目标函数无分段线性化的情况下模型的计算效率至少提高了10％，在目标函数分段线性化的情况下，设置求解的精度为0.5％，模型的计算效率至少提高了20％。
图1示出了本发明实施例提供的一种多时段高维投影机组组合模型的构建方法的实现流程。为了便于说明，仅示出了与本发明实施例相关的部分，详述如下：
如图1所示，本发明实施例提供了一种多时段高维投影机组组合模型的构建方法，包括步骤S102、步骤S104和步骤S106。
步骤S102，采集电力系统发电机机组的基础数据，所述基础数据包括发电机机组的运行特性数据、负荷预测数据以及各时段的旋转备用数据。
步骤S104，利用变量投影技术将电机机组功率变量以及与其相关的参数进行投影，得到投影后的电机机组功率变量以及与其相关的参数，其中投影区间为[0,1]。
在本发明实施例中，利用变量投影技术将机组功率变量投影到区间[0,1]。机组i在t时段的发电功率Pi，t是一个半连续变量，即其中P i表示机组i的出力下界，表示机组i的出力上界。现令对功率变量Pi，t进性投影得到且为了便于理解，可以将看作为机组i在t时段的出力在中所占的比例。同时，其他与机组功率相关的参数也需要投影，如：其中，Pup，i和Pdown，i分别表示机组i向上或向下爬坡的限制，Pstart，i和Pshut，i分别表示机组i开启或关闭的最小功率。
步骤S106，根据所述基础数据、投影后的电机机组功率变量以及与其相关的参数以及电机机组在运行、启动和关机三个时段内单个机组的所有运行状态变量，构建多时段高维投影机组组合模型。
在本发明实施例中，电机机组在运行、启动和关机三个时段内单个机组的所有运行状态变量之间的关系如下表1所示。
表1
其中，ui，t表示机组i在t时段的运行状态(0表示关机，1表示开机)，si，t表示机组i在t时段的启动状态(1表示机组出现由关机到开机的状态转变，0表示其他情况)，di，t表示t时段机组i关机状态(1表示机组出现由开机到关机的状态转变，0表示其他情况)。由于表1仅展示单个机组的状态变量，故变量中的机组号i被省略。这八类0-1变量表示机组i在t-1、t和t+1三个时段的八种不同的状态，例如：若ui，t-1＝ui，t＝ui，t+1＝0，则否则，若ui，t-1＝ui，t＝0，ui，t+1＝1，则否则，同时，需要注意的是，由于表1仅展示单个机组的状态变量，故变量中的机组号i被省略。
本发明实施例提供的多时段高维投影机组组合模型的构建方法，通过将传统电力系统UC模型中的机组出力变量由半连续变量转化为连续变量，使UC模型中的出力约束式的数量至少减少了一半，降低了机组模型求解的复杂度；并且利用投影技术对目标函数进行变化，使目标函数值域范围更广；基于CPLEX12.7.1工具的模型求解表明，采用本发明的构建方法构建得到的三时段机组组合模型，在目标函数无分段线性化的情况下模型的计算效率至少提高了10％，在目标函数分段线性化的情况下，设置求解的精度为0.5％，模型的计算效率至少提高了20％。可见，本发明所构建得到的多时段高维投影机组组合模型的求解难度降低，且计算效率高，可用于控制和调度电网发电资源，能够同时满足机组完成合同电量的实际需求，以达到资源优化配置和节能减排的目标。
在本发明实施例中，上述步骤S106包括：
根据所述基础数据、投影后的电机机组功率变量以及与其相关的参数以及电机机组在运行、启动和关机三个时段内单个机组的所有运行状态变量，构建所述多时段高维投影机组组合模型中的目标函数；确定所述多时段高维投影机组组合模型中的所述目标函数的约束条件。
在本发明实施例中，基于机组的实际运行需求利用新状态变量系统地推导出三时段内的机组出力上界约束和爬坡约束的理想表达形式。根据机组的实际运行要求，分别得到下表2、表3和表4的数据。
表2三时段单机组出力上界
表3三时段单机组上爬坡上限
表4三时段单机组下爬坡上限
根据表1和表2可知，在t时段机组的状态变量与功率上界之间存在如下线性关系：
即，
根据表1中各变量之间的关系将和消去，得到
由于以上式子是根据机组在t时段的实际运行要求推导而成，因此该式子是一个十分理想的机组出力上界约束表达式。根据以上思路，可以推导出机组在t-1和t+1时段的功率上界约束的理想表达式。同样，三时段的机组上爬坡限制约束的理想表达式也可以按照上述思路得出。
利用新状态变量与其他状态变量之间的关系消除新状态变量以提升模型的简洁性。例如，根据表1和表2，机组在t时段的上界约束表达式中的可以被消去，并得到如下表达式：
根据表1和表3可知，在机组的状态变量与上爬坡上限之间存在如下线性关系：
即
同样，由表1可知，和可以消去，可得
同样也可以推导出机组其他时段的下爬坡约束(和)的表达式。
根据表1和表4可知，在机组的状态变量与下爬坡上限之间存在如下线性关系：
即
根据以上思路和可以消去，
同样也可以推导出机组其他时段的上下爬坡约束( 和)的表达式。
利用新状态变量与其他状态变量之间的关系消除新状态变量以提升模型的简洁性。根据表1可知，
由此，机组出力上界约束：
其中，
机组出力下界约束：
向上爬坡约束：
向下爬坡约束：
其中，
在构建本发明的多时段高维投影机组组合模型中，还需要加入初始状态约束、机组开机费用约束、最小启/停约束、状态约束、系统旋转热备用约束和功率平衡约束；同时，需要加入二次目标函数。
其中，与本发明相关的约束条件如下：
机组状态约束：si，t-di，t＝ui，t-uy，t-1，其中，si，t表示机组i在t时刻开机状态(1表示机组出现由关机到开机的状态转变，0表示其他情况)，di，t表示t时段机组i关机状态(1表示机组出现由开机到关机的状态转变，0表示其他情况)。
最小启停时间约束包括开机时间约束和关机时间约束。
开机时间约束：
关机时间约束：
上述约束中，其中表示机组i在时刻开机。
机组初始启动状态，构建二类状态变量的机组初始状态约束：ui，t＝ui，0，t∈[1，...，Ui+Li]，上述约束中，Ui＝[min[T，ui，0(T on，i-Ti，0)]]+，Li＝[min[T，(1-ui，0)(T off，i+Ti，0)]]+。其中，T表示调度总时段数，T on，i表示机组i的最小开机时间，T off，i表示机组i的最小停机时间，ui，0表示机组i的初始状态，Ti，0表示机组i初始时已经运行或停机的时间，Ui表示机组i在初始时刻仍需运行时间，Li表示机组i在初始时刻仍需停机时间。
根据负荷预测数据和旋转备用数据，构建机组旋转备用约束如下：其中i表示机组下标，t表示时段下标，N表示机组总数，ui，t表示机组i在t时刻的运行状态，表示机组i的出力上界，PD，t表示t时段时电力系统负荷，Rt表示t时段电力系统所需旋转备用值。
根据机组组合基础数据和投影的机组出力连续变量，建立机组功率平衡约束：其中P i表示机组i的出力下界。
利用机组冷热启动费用差异特性，建立机组启动费用约束：
上述表达式中，若t-T off，i-Tcold，i-1≤0且[-Ti，0]+＜|t-T off，i-Tcold，i-1|+1，则f′init，i，t＝1，否则f′init，i，t＝0；其中表示机组i在t时刻热启动费用超出的部分，且Chot，i表示机组i的热启动费用，Ccold，i表示机组i的冷启动费用，Tcold，i表示计算机组i冷启动时间，f′init，i，t表示机组i在t时刻计及热启动费用超出的部分。
结合启动费用约束，建立二类状态变量机组模型的目标表达式如下：
表达式中 和其中FC是优化目标，表示机组i的出力成本，即机组燃料费用，αi，βi，γi表示机组i的二次燃料费用函数系数，表示机组i进行投影变换的二次燃料费用函数系数，表示机组i在t时刻热启动费用超出的部分，Chot，i表示机组i的热启动费用。
经过上述步骤之后，构建得到本发明实施例提供的多时段高维投影机组组合模型的目标函数：
其中，所述目标函数中的约束条件如下：
(1)机组出力约束：
其中，Pi，t表示机组i在t时刻的出力大小， Pup，i和Pdown，i分别表示机组i向上或向下爬坡的限制，Pstart，i和Pshut，i分别表示机组i开启或关闭的最小功率，
(2)功率平衡约束：
(3)旋转备用约束：
(4)爬坡速率约束：
其中，[·]+表示max(0，·)，Pi，t表示机组i在t时刻的出力大小， Pup，i和Pdown，i分别表示机组i向上或向下爬坡的限制，Pstart，i和Pshut，i分别表示机组i开启或关闭的最小功率，表示机组i的出力上界，P i表示机组i的出力下界，ui，0表示机组i的初始状态，Ti，0表示机组i初始时已经运行或停机的时间，ui，t表示机组i在t时刻的运行状态(1表示机组处于运行状态，0表示停机状态)，si，t表示机组i在t时刻开机状态(1表示机组出现由关机到开机的状态转变，0表示其他情况)，di，t表示t时段机组i关机状态(1表示机组出现由开机到关机的状态转变，0表示其他情况)。
(5)机组状态约束：
ui，t-ui，t-1≤si，t
(5)机组初始状态约束：ui，t＝ui，0，t∈[1，...，Ui+Li]，上述约束中，Ui＝[min[T，ui，0(T on，i-Ti，0)]]+，Li＝[min[T，(1-ui，0)(T off，i+Ti，0)]]+。
(7)机组最小启停时间约束：
开机时间约束：
关机时间约束：
(8)机组启动费用约束：
上述约束公式中，FC是优化目标，表示机组i的出力成本，即机组燃料费用，i表示机组下标，t表示时段下标，N表示机组总数，T表示调度总时段数，αi，βi，γi表示机组i的二次燃料费用函数系数，表示机组i进行投影变换的二次燃料费用函数系数，Chot，i表示机组i的热启动费用，Ccold，i表示机组i的冷启动费用，T on，i表示机组i的最小开机时间，T off，i表示机组i的最小停机时间，Tcold，i表示计算机组i冷启动时间，f′init，i，t表示机组i在t时刻计及热启动费用超出的部分，[·]+表示max(0，·)，表示机组i的出力上界，P i表示机组i的出力下界，表示机组模型的机组i在t时刻的出力大小，PD，t表示t时段时电力系统所需负荷，Rt表示t时段电力系统所需旋转备用值，Pup，i表示机组i的向上爬坡速率，表示机组i进行投影变换后的向上爬坡速率，Pdown，i表示机组i的向下爬坡速率，表示机组i进行投影变换后的向下爬坡速率，Pstart，i表示机组i开机时的最小出力值，表示机组i进行投影变换后的在开机时的最小出力值，Pshut，i表示机组i关机时的最大出力值，表示机组i进行投影变换后的在关机时的最大出力值，ui，0表示机组i的初始状态，Ti，0表示机组i初始时已经运行或停机的时间，ui，t表示机组i在t时刻的运行状态，si，t表示机组i在t时刻开机，表示机组i在t时刻热启动费用超出的部分，表示机组i在时刻开机。
综上所述，本发明实施例提供的多时段高维投影机组组合模型的构建方法，通过将传统电力系统UC模型中的机组出力变量由半连续变量转化为连续变量，使传统的UC模型中的出力约束式的数量至少减少了一半，降低了机组模型求解的复杂度；并利用投影技术对目标函数进行变化，使目标函数值域范围更广；新构建的爬坡速率约束减少了冗余约束，并提出了一种更紧凑的三时段机组组合出力约束和爬坡约束。另外，基于CPLEX12.7.1工具对模型进行求解表明，本发明的三时段机组组合模型，在目标函数无分段线性化的情况下模型的计算效率至少提高了10％，在目标函数分段线性化的情况下，设置求解的精度为(0.5％),模型的计算效率至少提高了20％。因此，本发明所构建得到的多时段高维投影机组组合模型的求解难度降低，且计算效率高，可用于控制和调度电网发电资源，能够同时满足机组完成合同电量的实际需求，以达到资源优化配置和节能减排的目标。
以下通过具体的实例来证明本发明构建得到的模型的可行性。
采集电力系统发电机机组的基础数据，包括发电机机组的运行特性数据、负荷预测数据如下表5、6、7。
表5发电机数据(一)
表6发电机数据(二)
表7各时段负荷
| 时段 | 有功负荷(MW) |
| 1 | 178.69 |
| 2 | 168.45 |
| 3 | 161.84 |
| 4 | 157.83 |
| 5 | 158.16 |
| 6 | 163.69 |
| 7 | 176.86 |
| 8 | 198.21 |
| 9 | 209.67 |
| 10 | 223.54 |
| 11 | 233.18 |
| 12 | 240.80 |
| 13 | 247.03 |
| 14 | 248.47 |
| 15 | 253.83 |
| 16 | 270.90 |
| 17 | 290.12 |
| 18 | 251.68 |
| 19 | 250.89 |
| 20 | 242.10 |
| 21 | 242.05 |
| 22 | 231.68 |
| 23 | 198.07 |
| 24 | 190.69 |
各时段的旋转备用数据为所有机组总功率的10％。
按照本发明的上述建模步骤，建立出多时段高维投影机组组合模型，并经基于CPLEX12.7.1工具的模型求解表明，采用本发明的构建方法构建得到的三时段机组组合模型，在目标函数无分段线性化的情况下模型的计算效率至少提高了10％，在目标函数分段线性化的情况下，设置求解的精度为0.5％，模型的计算效率至少提高了20％。
图2是本发明实施例提供的一种多时段高维投影机组组合模型的构建装置的结构框图，为了便于说明，图中仅示出与本发明实施例相关的部分，详述如下：
参见图2，本发明实施例提供的多时段高维投影机组组合模型的构建装置包括数据采集单元210、变量投影单元220和模型构建单元230。
数据采集单元210，用于采集电力系统发电机机组的基础数据，所述基础数据包括发电机机组的运行特性数据、负荷预测数据以及各时段的旋转备用数据。
变量投影单元220，用于利用变量投影技术将电机机组功率变量以及与其相关的参数进行投影，得到投影后的电机机组功率变量以及与其相关的参数，其中投影区间为[0,1]。
模型构建单元230，用于根据所述基础数据、投影后的电机机组功率变量以及与其相关的参数以及电机机组在运行、启动和关机三个时段内单个机组的所有运行状态变量，构建多时段高维投影机组组合模型。
本发明实施例提供的多时段高维投影机组组合模型的构建装置，通过将传统电力系统UC模型中的机组出力变量由半连续变量转化为连续变量，使UC模型中的出力约束式的数量至少减少了一半，降低了机组模型求解的复杂度；并且利用投影技术对目标函数进行变化，使目标函数值域范围更广；基于CPLEX12.7.1工具的模型求解表明，采用本发明的构建方法构建得到的三时段机组组合模型，在目标函数无分段线性化的情况下模型的计算效率至少提高了10％，在目标函数分段线性化的情况下，设置求解的精度为0.5％，模型的计算效率至少提高了20％。可见，本发明所构建得到的多时段高维投影机组组合模型的求解难度降低，且计算效率高，可用于控制和调度电网发电资源，能够同时满足机组完成合同电量的实际需求，以达到资源优化配置和节能减排的目标。
图3是本发明实施例提供的一种模型构建单元的结构框图，为了便于说明，图中仅示出与本发明实施例相关的部分，详述如下：
所述模型构建单元230包括目标函数构建模块231和约束条件确定模块232。
目标函数构建模块231，用于根据所述基础数据、投影后的电机机组功率变量以及与其相关的参数以及电机机组在运行、启动和关机三个时段内单个机组的所有运行状态变量，构建所述多时段高维投影机组组合模型中的目标函数。
约束条件确定模块232，用于确定所述多时段高维投影机组组合模型中的所述目标函数的约束条件。
关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。
应该理解的是，虽然本发明各实施例的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，各实施例中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一非易失性计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。
Claims (1)
1.多时段高维投影的机组组合模型的构建方法和装置，包括如下步骤：
(1)收集电力系统发电机机组基础数据；
所述发电机机组基础数据包括发电机机组的运行特性数据、负荷预测数据以及各时段的旋转备用数据；
(2)利用变量投影技术将机组功率变量投影到[0,1]区间，机组i在t时段的发电功率Pi,t是一个半连续变量，即其中Pi 表示机组i的出力下界，表示机组i的出力上界，现令对功率变量Pi,t进性投影得到且 可视为机组i在t时段的出力在中所占的比例，其他与机组功率相关的参数也需要投影， 其中， Pup,i和Pdown,i分别表示机组i向上或向下爬坡的限制，Pstart,i和Pshut,i分别表示机组i开启或关闭的最小功率；
(3)引入新的状态变量来表达多时段内单个机组的所有运行状态，表1中给出了三时段内单个机组的所有运行状态和新引入的状态变量
表1三时段单机组各个状态变量之间的关系
其中，ui,t表示机组i在t时段的运行状态，0表示关机，1表示开机，si,t表示机组i在t时段的启动状态，1表示机组出现由关机到开机的状态转变，0表示其他情况，di,t表示t时段机组i关机状态，1表示机组出现由开机到关机的状态转变，0表示其他情况，这八类0-1变量表示机组i在t-1、t和t+1三个时段的八种不同的状态；
(4)基于机组的实际运行需求利用新状态变量系统地推导出多时段内的机组出力上界约束和爬坡约束的理想表达形式，表2给出了三时段的机组出力上界，表3给出了三时段的机组爬坡上限；
表2三时段单机组出力上界
表3三时段单机组上爬坡限制
根据表1和表2可知，在t时段机组的状态变量与功率上界之间存在如下线性关系：
即，
根据表1中各变量之间的关系将和消去，得到
由以上推导过程可知，所得表达式为理想的机组出力上界约束表达式；
根据以上思路，可以推导出机组在t-1和t+1时段的功率上界约束的理想表达式，并且三时段的机组上爬坡限制约束的理想表达式也可以按照上述思路得出；
(5)利用新状态变量与其他状态变量之间的关系消除新状态变量以提升模型的简洁性，根据表1和表2，机组在t时段的上界约束表达式中的可以被消去，并得到如下机组出力约束和爬坡速率约束表达式：
出力约束表达式：
爬坡速率约束：
上述约束公式中,i表示机组下标，t表示时段下标,N表示机组总数，T表示调度总时段数，Chot,i表示机组i的热启动费用，Ccold,i表示机组i的冷启动费用，T on,i表示机组i的最小开机时间，T off,i表示机组i的最小停机时间，Tcold,i表示计算机组i冷启动时间，[·]+表示max(0,·)，Pi,t表示机组i在t时刻的出力大小，Pup,i和Pdown,i分别表示机组i向上或向下爬坡的限制，Pstart,i和Pshut,i分别表示机组i开启或关闭的最小功率，表示机组i的出力上界，Pi 表示机组i的出力下界，ui,0表示机组i的初始状态，Ti,0表示机组i初始时已经运行或停机的时间，ui,t表示机组i在t时刻的运行状态，1表示机组处于运行状态，0表示停机状态，si,t表示机组i在t时刻开机状态，1表示机组出现由关机到开机的状态转变，0表示其他情况，di,t表示t时段机组i关机状态，1表示机组出现由开机到关机的状态转变，0表示其他情况；
(6)加入其他约束和目标函数：
为了构建一个完整的机组组合模型，还需要加入初始状态约束、机组开机费用约束、最小启/停约束、状态约束、系统旋转热备用约束和功率平衡约束，同时，还要加入二次目标函数；
(7)最终得到一个多时段高维投影的机组组合模型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011208689.9A CN112258076B (zh) | 2020-11-03 | 2020-11-03 | 一种多时段高维投影机组组合模型的构建方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011208689.9A CN112258076B (zh) | 2020-11-03 | 2020-11-03 | 一种多时段高维投影机组组合模型的构建方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112258076A CN112258076A (zh) | 2021-01-22 |
| CN112258076B true CN112258076B (zh) | 2023-12-12 |
Family
ID=74268603
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011208689.9A Active CN112258076B (zh) | 2020-11-03 | 2020-11-03 | 一种多时段高维投影机组组合模型的构建方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112258076B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106655153A (zh) * | 2016-10-11 | 2017-05-10 | 广西大学 | 二类状态变量的电力系统机组组合模型及其建立方法 |
| CN106712116A (zh) * | 2017-03-23 | 2017-05-24 | 广西大学 | 完全分布式的电力系统机组投入配置方法及系统 |
| CN109149571A (zh) * | 2018-09-21 | 2019-01-04 | 国网福建省电力有限公司 | 一种考虑系统燃气和火电机组特性的储能优化配置方法 |
| CN109830991A (zh) * | 2019-04-02 | 2019-05-31 | 广西大学 | 两阶段分布式电力系统机组投入配置方法及系统 |
| CN111327052A (zh) * | 2020-03-12 | 2020-06-23 | 清华大学 | 一种电力系统机组组合加速优化方法及装置 |
-
2020
- 2020-11-03 CN CN202011208689.9A patent/CN112258076B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106655153A (zh) * | 2016-10-11 | 2017-05-10 | 广西大学 | 二类状态变量的电力系统机组组合模型及其建立方法 |
| CN106712116A (zh) * | 2017-03-23 | 2017-05-24 | 广西大学 | 完全分布式的电力系统机组投入配置方法及系统 |
| CN109149571A (zh) * | 2018-09-21 | 2019-01-04 | 国网福建省电力有限公司 | 一种考虑系统燃气和火电机组特性的储能优化配置方法 |
| CN109830991A (zh) * | 2019-04-02 | 2019-05-31 | 广西大学 | 两阶段分布式电力系统机组投入配置方法及系统 |
| CN111327052A (zh) * | 2020-03-12 | 2020-06-23 | 清华大学 | 一种电力系统机组组合加速优化方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112258076A (zh) | 2021-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | A reserve capacity model of AA-CAES for power system optimal joint energy and reserve scheduling | |
| CN105528466B (zh) | 考虑电力系统适应性和经济性的风电优化规划建模方法 | |
| CN107292766B (zh) | 面向风电消纳的电力系统调峰手段经济性评估方法与系统 | |
| CN109472401B (zh) | 分布式供能运行策略的确定方法、装置、设备和存储介质 | |
| CN111355265B (zh) | 一种微电网能量两阶段鲁棒优化方法及系统 | |
| CN107944733B (zh) | 一种含风电、需求响应和压缩空气储能的电网调度方法 | |
| CN103855718B (zh) | 抽水蓄能电站参与含风电电力系统的调度方法 | |
| CN107800153B (zh) | 一种电蓄热消纳风电的电热能量滚动鲁棒调度方法 | |
| Zheng et al. | Optimal design for a multi-level energy exploitation unit based on hydrogen storage combining methane reactor and carbon capture, utilization and storage | |
| CN115528674A (zh) | 基于减载运行的海上风电场虚拟惯量最优分配方法及系统 | |
| CN111126675A (zh) | 多能互补微网系统优化方法 | |
| CN112668751B (zh) | 一种机组优化调度模型的建立方法及装置 | |
| Zhang et al. | Modeling and optimization of integrated energy system for renewable power penetration considering carbon and pollutant reduction systems | |
| CN115081700A (zh) | 基于综合储能技术的数据中心多能协同优化方法及系统 | |
| CN114928117A (zh) | 一种基于场景聚焦的新能源承载能力评估方法及装置 | |
| Li et al. | Multi-mode optimal operation of advanced adiabatic compressed air energy storage: Explore its value with condenser operation | |
| BR102022003674A2 (pt) | Método e dispositivos para controlar um complexo de instalação de gás industrial | |
| CN112258076B (zh) | 一种多时段高维投影机组组合模型的构建方法和装置 | |
| CN117613957A (zh) | 一种电-氢-氨的最优动态调节控制方法及装置 | |
| Li et al. | Modeling of integrated energy system of offshore oil and gas platforms considering couplings between energy supply system and oil and gas production system | |
| CN113036817A (zh) | 一种基于槽式光热电站的电力系统优化调度方法及系统 | |
| Klymchuk et al. | Enhancing efficiency of using energy resources in heat supply systems of buildings with variable operation mode | |
| CN112307603A (zh) | 考虑大规模风电接入的混合储能容量优化配置方法及系统 | |
| CN116613829A (zh) | 含规模化多类型调节资源参与的电力系统运行优化方法及系统 | |
| CN114389262B (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 |