CN103490449A

CN103490449A - Optimization method for multi-energy combined power generation system operation simulation

Info

Publication number: CN103490449A
Application number: CN201310469621.XA
Authority: CN
Inventors: 赵冬梅; 李星宇; 牟澎涛; 俞勤政; 李龙龙; 王建锋; 郭威; 金小明; 张东辉; 魏国清; 胡剑琛; 祁永福; 吴锋
Original assignee: China South Power Grid International Co ltd; North China Electric Power University; Hainan Power Grid Co Ltd
Current assignee: China South Power Grid International Co ltd; North China Electric Power University; Hainan Power Grid Co Ltd
Priority date: 2013-10-10
Filing date: 2013-10-10
Publication date: 2014-01-01
Anticipated expiration: 2033-10-10
Also published as: CN103490449B

Abstract

The invention discloses the technical field of power system operation and scheduling, and in particular relates to an optimization method for operation simulation of a multi-energy combined power generation system. Select a few days or a week as the scheduling cycle, and establish the expected scenario of wind power and photovoltaic energy output according to the predicted deviation value of wind power, photovoltaic and other intermittent power sources; in the multi-energy joint power generation system, first determine the hydropower unit whose output can be flexibly adjusted The working position in the load curve, and then starting from the initial period of the dispatch cycle, arrange the combination of thermal power units and allocate the load reasonably to various types of units with the local goal of optimizing the operation cost, and then form a sufficient load regulation for the operating units in each period Finally, by continuously selecting and updating the load regulation adequacy interval of units in each period, the unit combination arrangement and load distribution among units in the entire dispatching cycle are finally realized. The invention ensures the power balance and stable operation of the system.

Description

An optimization method for operation simulation of multi-energy combined power generation system

技术领域technical field

本发明属于电力系统运行和调度技术领域，尤其涉及一种多能源联合发电系统运行模拟的优化方法。The invention belongs to the technical field of power system operation and scheduling, and in particular relates to an optimization method for operation simulation of a multi-energy combined power generation system.

背景技术Background technique

近年来，随着国内电力工业的长足发展和市场化改革的进行，我国电力系统迅猛发展，系统规模日益庞大，运行环境越发复杂。与此同时，风电场、光伏电站等新能源机组在电力系统中的比例逐步增大，这使得电源结构发生明显改变，电力系统正过渡为含有火电、水电、核电、风电、光伏、储能等多类型电源的联合发电系统。部分地区的风电的渗透率已不容忽视，由于风电、光伏出力的随机性和间歇性，在接入系统后，会使系统面临更多的不确定性因素。因此，传统的电力系统调度运行理论需要加以改进和完善，以适应系统电源结构的改变。在这种形势下，提出一种多能源联合发电系统运行模拟的优化方法，在数日或一周的调度周期内研究系统机组组合与调度问题，使系统内各类型电源的出力能够满足系统用电负荷的需求，保证系统功率平衡与稳定运行。In recent years, with the rapid development of the domestic electric power industry and the market-oriented reform, my country's electric power system has developed rapidly, the scale of the system has become larger and larger, and the operating environment has become more and more complex. At the same time, the proportion of new energy units such as wind farms and photovoltaic power stations in the power system is gradually increasing, which has significantly changed the power structure. The power system is transitioning to include thermal power, hydropower, nuclear power, wind power, photovoltaics, energy storage, etc. Combined power generation system of multiple types of power sources. The penetration rate of wind power in some areas cannot be ignored. Due to the randomness and intermittent nature of wind power and photovoltaic output, after connecting to the system, the system will face more uncertain factors. Therefore, the traditional power system dispatching operation theory needs to be improved and perfected to adapt to the change of system power structure. In this situation, an optimization method for the operation simulation of a multi-energy combined power generation system is proposed, and the system unit combination and scheduling problems are studied within a scheduling period of several days or a week, so that the output of various types of power sources in the system can meet the power consumption of the system. Load demand, to ensure system power balance and stable operation.

发明内容Contents of the invention

本发明的目的在于，提供一种多能源联合发电系统运行模拟的优化方法，用于解决风电机组和光伏发电机组接入系统后，由于二者的随机性和间歇性，导致的发电系统不稳定的问题，从而保证系统功率平衡与稳定运行。The purpose of the present invention is to provide an optimization method for the operation simulation of a multi-energy combined power generation system, which is used to solve the instability of the power generation system caused by the randomness and intermittent nature of the wind turbines and photovoltaic generators after they are connected to the system problems, so as to ensure the system power balance and stable operation.

为了实现上述目的，本发明提出的技术方案是，一种多能源联合发电系统运行模拟的优化方法，其特征是所述方法包括：In order to achieve the above object, the technical solution proposed by the present invention is an optimization method for operation simulation of a multi-energy combined power generation system, which is characterized in that the method includes:

步骤1：设定调度周期并将所述调度周期等分成T个时段；Step 1: Set a scheduling cycle and divide the scheduling cycle into T time periods;

步骤2：设定每个时段抽样场景数目N，并确定各个时段每个抽样场景的风电机组和光电机组混合出力值

其中，m=1，2，...，N，j=1，2，...，T；Step 2: Set the number N of sampling scenarios in each period, and determine the combined output value of wind turbines and photovoltaic units in each sampling scenario in each period

Among them, m=1, 2,..., N, j=1, 2,..., T;

步骤3：根据公式P_L1，j=P_L，j-P_f，j计算每个时段系统中可调节型发电机组应承担的负荷P_L1，j，并将P_L1，j记为第一等值负荷；Step 3: According to the formula P _{L1, j} = P _{L, j} - P _{f, j,} calculate the load P _L1, j that the adjustable generator set in the system should bear in each period, and record P _{L1, j} as the first class value load;

所述系统中可调节型发电机组包括可调节型水电机组和可调节型火电机组；The adjustable generating unit in the system includes an adjustable hydroelectric unit and an adjustable thermal power unit;

P_L，j为第j个时段系统应承担的总负荷；P _{L, j} is the total load that the system should bear in the jth time period;

P_f，j为第j个时段系统中除可调节型发电机组之外的发电机组提供的负荷；P _{f, j} is the load provided by the generating units other than the adjustable generating units in the jth period of time;

步骤4：根据公式

计算每个时段系统中可调节型火电机组应承担的负荷P_L2，j，并将P_L2，j记为第二等值负荷；同时，计算每个时段可调节型水电机组向上调节极限值和向下调节极限值

Step 4: According to the formula

Calculate the load P _L2,j that the adjustable thermal power unit in the system should bear in each period, and record P _L2,j as the second equivalent load; at the same time, calculate the upward adjustment limit value of the adjustable hydroelectric unit in each period and down adjustment limit

μ_Hi，j为第j个时段第Hi个可调节型水电机组的工作状态，当第j个时段第Hi个可调节型水电机组投入运行时，μ_Hi，j=1；当第j个时段第Hi个可调节型水电机组未投入运行时，μ_Hi，j=0；μHi _,j is the working state of the Hi adjustable hydroelectric unit in the jth period, when the Hi adjustable hydroelectric unit is put into operation in the jth period, μHi _,j =1; when the jth period When the Hi-th adjustable hydroelectric unit is not in operation, μHi _{, j} =0;

P_Hi，j为第j个时段第Hi个可调节型水电机组的输出功率且P_Hi，j=η_Hi×g×(H_up，Hi，j-H_{down，Hi，j})×Q_Hi，j；P _Hi,j is the output power of the Hi-th adjustable hydroelectric unit in the j-th period and P _Hi,j =η _Hi ×g×(H _up,Hi,j -H _down,Hi,j )×Q _{Hi, j} ;

η_Hi为第Hi个可调节型水电机组的发电效率系数； _ηHi is the power generation efficiency coefficient of the Hi-th adjustable hydroelectric unit;

g为重力常数；g is the gravitational constant;

H_up，Hi，j为第j个时段第Hi个可调节型水电机组的坝上水位；H _{up, Hi, j} is the water level above the dam of the Hi-th adjustable hydroelectric unit in the j-th period;

H_{down，Hi，j}为第j个时段第Hi个可调节型水电机组的下游水位；H _{down, Hi, j} is the downstream water level of the Hi-th adjustable hydroelectric unit in the j-th period;

Q_Hi，j为第j个时段第Hi个可调节型水电机组的发电平均流量；Q _Hi,j is the average flow rate of the Hi-th adjustable hydroelectric unit in the j-th time period;

N_H为可调节型水电机组的数量；N _H is the number of adjustable hydroelectric units;

μ_Hp，j为第j个时段抽水蓄能发电机组的抽水状态，当第j个时段抽水蓄能发电机组投入运行时，μ_Hp，j=1，当第j个时段抽水蓄能发电机组未投入运行时，μ_Hp，j=0；μ _Hp,j is the pumping state of the pumped storage generator set in the jth period. When the pumped storage generator set is put into operation in the jth period, μ _{Hp, j} = 1. When put into operation, μ _{Hp, j} =0;

P_Hp，j为第j个时段抽水蓄能发电机组的抽水功率；P _Hp,j is the pumping power of the pumped storage generator set in the jth period;

步骤5：确定可调节型火电机组的初始投入运行状态；具体是，分别计算调度周期前5个时段的第二等值负荷最大值和前5个时段的第二等值负荷最小值，按照第一设定顺序投入可调节型火电机组直至投入运行的可调节型火电机组的总装机容量大于调度周期前5个时段的第二等值负荷最大值，并且投入运行的可调节型火电机组的最小技术出力总和小于调度周期前5个时段的第二等值负荷最小值；Step 5: Determine the initial operation status of the adjustable thermal power unit; specifically, calculate the second equivalent load maximum value and the first 5 period second equivalent load minimum value in the first 5 periods of the dispatch cycle, according to the first A set order to put the adjustable thermal power units into operation until the total installed capacity of the adjustable thermal power units put into operation is greater than the second equivalent load maximum value in the first 5 periods of the scheduling cycle, and the minimum value of the adjustable thermal power units put into operation The sum of technical output is less than the minimum value of the second equivalent load in the first 5 periods of the dispatch cycle;

所述可调节型火电机组包括不可停机的可调节型火电机组、可停机的可调节型火电机组和燃气轮机机组；The adjustable thermal power unit includes an adjustable thermal power unit that cannot be shut down, an adjustable thermal power unit that can be shut down, and a gas turbine unit;

步骤6：设定初始时段t=5；Step 6: Set the initial time period t=5;

步骤7：根据公式

计算当前时段t可调节型火电机组在每个抽样场景下应承担的负荷

并将

记为第三等值负荷；Step 7: According to the formula

Calculate the load that the adjustable thermal power unit should bear in each sampling scenario in the current period t

and will

Recorded as the third equivalent load;

步骤8：为当前时段t的每个抽样场景下已投入运行的可调节型火电机组分配负荷

同时确定当前时段t每个抽样场景下已投入运行的可调节型火电机组的负荷调节充裕度区间；Step 8: Assign loads to the adjustable thermal power units that have been put into operation in each sampling scenario of the current period t

At the same time, determine the load regulation adequacy interval of the adjustable thermal power units that have been put into operation under each sampling scenario in the current period t;

步骤9：计算当前时段t的每个抽样场景下已投入运行的可调节型火电机组的运行成本，选择运行成本的最小值对应的抽样场景τ下已投入运行的可调节型火电机组及其分配的负荷

作为当前时段t可调节型火电机组的最佳投入运行方案；同时，选择运行成本的最小值对应的抽样场景τ下已投入运行的可调节型火电机组的负荷调节充裕度区间作为当前时段t可调节型火电机组的最佳负荷调节充裕度区间；Step 9: Calculate the operating cost of the adjustable thermal power units that have been put into operation under each sampling scenario of the current period t, and select the adjustable thermal power units that have been put into operation under the sampling scenario τ corresponding to the minimum value of the operating cost and their allocation load

As the optimal operation plan for the adjustable thermal power unit in the current period t; at the same time, the load regulation adequacy range of the adjustable thermal power unit that has been put into operation under the sampling scenario τ corresponding to the minimum value of the operating cost is selected as the current period t can be The optimal load regulation margin range of regulating thermal power units;

步骤10：令t=t+1，计算当前时段和当前时段之前连续4个时段的第二等值负荷最大值，以及当前时段和当前时段之前连续4个时段的第二等值负荷最小值；Step 10: Let t=t+1, calculate the second equivalent load maximum value of the current period and the 4 consecutive periods before the current period, and the second equivalent load minimum value of the current period and the 4 consecutive periods before the current period;

如果当前时段的前一个时段的可调节型火电机组投入运行的总装机容量小于当前时段和当前时段之前连续4个时段的第二等值负荷最大值，则按照第一设定顺序投入还未投入运行的可调节型火电机组，直至投入运行的可调节型火电机组的总装机容量大于当前时段和当前时段之前连续4个时段的第二等值负荷最大值；If the total installed capacity of the adjustable thermal power units put into operation in the previous period of the current period is less than the second equivalent load maximum value of the current period and the 4 consecutive periods before the current period, it will be put into use according to the first setting order. Adjustable thermal power units in operation until the total installed capacity of the adjustable thermal power units put into operation is greater than the second equivalent load maximum value of the current period and the 4 consecutive periods before the current period;

如果当前时段和当前时段之前连续4个时段的第二等值负荷最大值小于当前时段的前一个时段的可调节型火电机组的最小技术出力总和，则按照第一设定顺序的倒序停止已经投入运行的可调节型火电机组，直至投入运行的可调节型火电机组的最小技术出力总和小于当前时段和当前时段之前连续4个时段的第二等值负荷最小值；If the maximum value of the second equivalent load in the current period and the 4 consecutive periods before the current period is less than the sum of the minimum technical output of the adjustable thermal power unit in the previous period of the current period, stop the input in the reverse order of the first setting order Adjustable thermal power units in operation until the sum of the minimum technical output of the adjustable thermal power units put into operation is less than the second equivalent load minimum value of the current period and the 4 consecutive periods before the current period;

步骤11：判断当前时刻投入运行的可调节型火电机组与当前时刻的前一时刻投入运行的可调节型火电机组是否相同，如果当前时刻投入运行的可调节型火电机组与当前时刻的前一时刻投入运行的可调节型火电机组相同，则执行步骤12；否则，执行步骤15；Step 11: Determine whether the adjustable thermal power unit put into operation at the current moment is the same as the adjustable thermal power unit put into operation at the moment before the current moment. If the adjustable thermal power units put into operation are the same, go to step 12; otherwise, go to step 15;

步骤12：计算当前时段t每个抽样场景下可调节型火电机组负荷变化量

如果ΔP>0，则执行步骤13；如果ΔP<0，则执行步骤14；如果ΔP=0，则执行步骤15；Step 12: Calculate the load variation of the adjustable thermal power unit under each sampling scenario in the current period t

If ΔP>0, execute step 13; if ΔP<0, execute step 14; if ΔP=0, execute step 15;

步骤13：按照第二设定顺序调整当前时段已投入运行且尚有可调节容量的可调节型发电机组，直至参与调整的可调节型发电机组满足 $ΔP = Σ_{Hi = 1}^{N_{H}} λ_{Hi} P_{Hi, t, u} + Σ_{Gk = 1}^{N_{G}} λ_{Gk} P_{Gk, t, u}^{m},$ 而后执行步骤15；Step 13: Adjust the adjustable generator sets that have been put into operation in the current period and still have adjustable capacity according to the second setting sequence until the adjustable generator sets participating in the adjustment meet the $ΔP = Σ_{hi = 1}^{N_{h}} λ_{hi} P_{hi, t, u} + Σ_{K = 1}^{N_{G}} λ_{K} P_{K, t, u}^{m},$ Then execute step 15;

其中，λ_Hi为当前时段第Hi台可调节型水电机组参与调整的状态，当前时段第Hi台可调节型水电机组参与调整时，λ_Hi=1；当前时段第Hi台可调节型水电机组未参与调整时，λ_Hi=0；Among them, _λHi is the state that the Hi adjustable hydroelectric unit participates in the adjustment in the current period. When the Hi adjustable hydroelectric unit participates in the adjustment in the current period, _λHi =1; the Hi adjustable hydroelectric unit in the current period does not When participating in the adjustment, λ _Hi =0;

P_Hi，t，u为当前时段第Hi台可调节型水电机组的出力增加变化量；P _{Hi, t, u} is the output increase variation of the Hi adjustable hydroelectric unit in the current period;

λ_Gk为当前时段投入运行的火电机组中第Gk台火电机组参与调整的状态，当前时段投入运行的火电机组中第Gk台火电机组参与调整时，λ_GK=1；当前时段投入运行的火电机组中第Gk台火电机组未参与调整时，λ_Gk=0；λ _Gk is the status of the Gkth thermal power unit participating in the adjustment among the thermal power units put into operation _in the current period. When the Gkth thermal power unit does not participate in the adjustment, λ _Gk =0;

为当前时段第m个抽样背景下投入运行的火电机组中第Gk台火电机组出力增加变化量；

Increase the variation in the output of the Gkth thermal power unit in the thermal power unit put into operation under the mth sampling background in the current period;

N_G为已投入运行的可调节型火电机组的数量；N _G is the number of adjustable thermal power units that have been put into operation;

步骤14：按照第二设定顺序调整当前时段已投入运行且尚有可调节容量的可调节型发电机组，直至参与调整的可调节型发电机组满足Step 14: Adjust the adjustable generator sets that have been put into operation in the current period and still have adjustable capacity according to the second setting order until the adjustable generator sets participating in the adjustment meet the

$| | ΔP ΔP | | = = {Σ Σ}_{Hi hi = = 11}^{{N N}_{H h}} {λ λ}_{Hi hi} {P P}_{Hi hi,, t t,, d d} + + {Σ Σ}_{Gk K = = 11}^{{N N}_{G G}} {λ λ}_{Gk K} {P P}_{Gk K,, t t,, d d}^{m m};;$

P_Hi，t，d为当前时段第Hi台可调节型水电机组的出力减少变化量；P _{Hi, t, d} is the output reduction variation of the Hi-th adjustable hydroelectric unit in the current period;

为当前时段第m个抽样背景下投入运行的火电机组中第Gk台火电机组出力减少变化量；

Reduce the variation in the output of the Gkth thermal power unit in the thermal power unit put into operation under the background of the mth sampling in the current period;

步骤15：为当前时段的每个抽样场景下已投入运行的可调节型火电机组分配负荷

同时计算当前时段每个抽样场景下已投入运行的可调节型火电机组的负荷调节充裕度区间；Step 15: Assign loads to the adjustable thermal power units that have been put into operation under each sampling scenario in the current period

At the same time, calculate the load regulation adequacy interval of the adjustable thermal power units that have been put into operation under each sampling scenario in the current period;

步骤16：计算当前时段t的每个抽样场景下已投入运行的可调节型火电机组的运行成本，选择运行成本的最小值对应的抽样场景τ下已投入运行的可调节型火电机组及其分配的负荷

作为当前时段t可调节型火电机组的最佳投入运行方案；同时，选择运行成本的最小值对应的抽样场景τ下已投入运行的可调节型火电机组的负荷调节充裕度区间作为当前时段t可调节型火电机组的最佳负荷调节充裕度区间；Step 16: Calculate the operating cost of the adjustable thermal power units that have been put into operation under each sampling scenario of the current period t, and select the adjustable thermal power units that have been put into operation under the sampling scenario τ corresponding to the minimum value of the operating cost and their allocation load

步骤17：判断t>N是否成立，如果t>N，则执行步骤18；否则，返回步骤10；Step 17: Judging whether t>N holds true, if t>N, execute step 18; otherwise, return to step 10;

步骤18：结束。Step 18: End.

所述确定各个时段每个抽样场景的风电机组和光电机组混合出力值包括如下子步骤：The determination of the mixed output value of wind turbines and photoelectric units for each sampling scene in each time period Including the following sub-steps:

子步骤101：分别根据风电机组和光电机组出力的历史数据，做出风电机组出力预测偏差的频率分布图和光电机组出力预测偏差的频率分布图；Sub-step 101: According to the historical data of the output of the wind turbine and the photoelectric unit, respectively, make a frequency distribution diagram of the output prediction deviation of the wind turbine and a frequency distribution diagram of the output prediction deviation of the photovoltaic unit;

子步骤102：通过拟合方法分别得到各个时段风电机组出力预测偏差概率密度函数和光电机组出力预测偏差概率密度函数；Sub-step 102: Obtain the probability density function of wind turbine output forecast deviation and the probability density function of photoelectric unit output forecast deviation probability density function in each period respectively by fitting method;

子步骤103：对各个时段风电机组出力预测偏差概率密度函数和光电机组出力预测偏差概率密度函数分别进行积分运算，得到各个时段的风电机组出力预测偏差概率密度分布函数和光电机组出力预测偏差概率密度分布函数；Sub-step 103: Integrate the wind turbine output forecast deviation probability density function and the photoelectric unit output forecast deviation probability density function in each period to obtain the wind turbine output forecast deviation probability density function and photovoltaic unit output forecast deviation probability density in each period Distribution function;

子步骤104：利用公式计算各个时段每个抽样场景下的基础采样值；其中，m=1，2，...，N，N为每个时段抽样场景数目，U为区间(0，1)上的随机数；Sub-step 104: using the formula Calculate the basic sampling value under each sampling scene in each time period; wherein, m=1, 2, ..., N, N is the number of sampling scenes in each time period, and U is a random number on the interval (0, 1);

子步骤105：利用公式

计算各个时段每个抽样场景下风电机组出力预测偏差，并利用公式

计算各个时段每个抽样场景下光电机组出力预测偏差；其中，为风电机组出力预测偏差概率密度分布函数的反函数，

为光电机组出力预测偏差概率密度分布函数的反函数；j=1，2，...，T，T为调度周期的时段数；Sub-step 105: Using the formula

Calculate the output forecast deviation of wind turbines under each sampling scenario in each period, and use the formula

Calculate the output prediction deviation of photoelectric units in each sampling scene in each period; among them, is the inverse function of the wind turbine output prediction deviation probability density distribution function,

is the inverse function of the output prediction deviation probability density distribution function of the photoelectric unit; j=1, 2,..., T, T is the number of time periods in the scheduling cycle;

子步骤106：根据公式的风电机组和光电机组混合出力值；Sub-step 106: according to the formula The combined output value of the wind turbine and photoelectric unit;

其中，P_W，j为各个时段风电机组出力预测值；P_P，j为各个时段光电机组出力预测值。Among them, P _W,j is the predicted value of wind turbine output in each period; P _P,j is the predicted output of photovoltaic unit in each period.

所述计算每个时段可调节型水电机组向上调节极限值

采用公式The calculation of the upward adjustment limit value of the adjustable hydroelectric unit in each period

use the formula

${P P}_{H h,, j j}^{op op} = = {Σ Σ}_{Hi hi = = 11}^{{N N}_{H h}} {μ μ}_{Hi hi,, j j} {P P}_{Hi hi,, j j}^{op op};;$

其中，

为第j个时段第Hi个可调节型水电机组的向上调节容量，且in,

is the upward adjustment capacity of the Hi-th adjustable hydroelectric unit in the j-th period, and

${P P}_{Hi hi,, j j}^{op op} = = {P P}_{Hi hi,, j j}^{max max} - - {P P}_{Hi hi,, j j};;$

为第j个时段第Hi个可调节型水电机组的最大输出功率；

is the maximum output power of the Hi-th adjustable hydroelectric unit in the j-th period;

P_Hi，j为第j个时段第Hi个可调节型水电机组的输出功率；P _{Hi, j} is the output power of the Hi-th adjustable hydroelectric unit in the j-th period;

μ_Hi，j为第j个时段第Hi个可调节型水电机组的工作状态，当第j个时段第Hi个可调节型水电机组投入运行时，μ_Hi，j=1；当第j个时段第Hi个可调节型水电机组未投入运行时，μ_Hi，j=0。μHi _,j is the working state of the Hi adjustable hydroelectric unit in the jth period, when the Hi adjustable hydroelectric unit is put into operation in the jth period, μHi _,j =1; when the jth period When the Hi-th adjustable hydroelectric unit is not in operation, μHi _,j =0.

所述计算每个时段可调节型水电机组向下调节极限值

采用公式The calculation of the downward adjustment limit value of the adjustable hydroelectric unit in each period

use the formula

${P P}_{H h,, j j}^{ne ne} = = {Σ Σ}_{Hi hi = = 11}^{{N N}_{H h}} {μ μ}_{Hi hi,, j j} {P P}_{Hi hi,, j j}^{ne ne};;$

其中，

为第j个时段第Hi个可调节型水电机组的向下调节容量，且in,

is the downward adjustment capacity of the Hi-th adjustable hydroelectric unit in the j-th period, and

${P P}_{Hi hi,, j j}^{ne ne} = = {P P}_{Hi hi,, j j} - - {P P}_{Hi hi,, j j}^{min min};;$

为第j个时段第Hi个可调节型水电机组的最小强迫输出功率；

is the minimum forced output power of the Hi-th adjustable hydroelectric unit in the j-th period;

所述按照第一设定顺序投入可调节型火电机组具体是，先按最小运行比耗量由低到高的顺序投入不可停机的可调节型火电机组，再按最小比耗量由低到高的顺序投入可停机的可调节型火电机组，最后按输出功率调节速率由高到低的顺序投入燃气轮机机组。The said putting into operation the adjustable thermal power units according to the first setting sequence is specifically to start the non-stop adjustable thermal power units according to the order of the minimum operating specific consumption from low to high, and then according to the minimum specific consumption from low to high The adjustable thermal power units that can be shut down are put into operation in the order of the output power adjustment rate, and the gas turbine units are put into use in the order of output power adjustment rate from high to low.

所述为当前时段t的每个抽样场景下已投入运行的可调节型火电机组分配负荷

具体包括如下子步骤：The load allocation for the adjustable thermal power units that have been put into operation in each sampling scenario of the current period t

Specifically include the following sub-steps:

子步骤201：设定已投入运行的可调节型火电机组出力越限标志nFlag=1；Sub-step 201: set the output limit flag nFlag=1 of the adjustable thermal power unit that has been put into operation;

子步骤202：根据公式

计算火电机组的耗量微增率；其中，

为第三等值负荷且m=1，2，...，N；Sub-step 202: according to the formula

Calculate the slight increase rate of consumption of thermal power units; where,

is the third equivalent load and m=1, 2,..., N;

N为每个时段抽样场景数目；N is the number of sampling scenes in each period;

N′_G为当前时段t已投入运行且出力未超限的可调节型火电机组的数量；N′ _G is the number of adjustable thermal power units that have been put into operation in the current period t and whose output has not exceeded the limit;

A_Gk′为当前时段t第Gk′台已投入运行且出力未超限的可调节型火电机组的耗量特性二次函数的二次项系数；A _Gk' is the quadratic term coefficient of the quadratic function of the consumption characteristic of the adjustable thermal power unit whose Gk'th unit has been put into operation in the current period t and the output has not exceeded the limit;

B_Gk′为当前时段t第Gk′台已投入运行且出力未超限的可调节型火电机组的耗量特性二次函数的一次项系数；B _Gk' is the first-order coefficient of the consumption characteristic quadratic function of the adjustable thermal power unit whose Gk'th unit has been put into operation in the current period t and the output has not exceeded the limit;

子步骤203：令Gk=1；Sub-step 203: let Gk=1;

子步骤204：根据公式

计算当前时段t的每个抽样场景下第Gk台已投入运行的可调节型火电机组出力值；Sub-step 204: according to the formula

Calculate the output value of the Gkth adjustable thermal power unit that has been put into operation under each sampling scenario in the current period t;

A_Gk为当前时段t第Gk台已投入运行的可调节型火电机组的耗量特性二次函数的二次项系数；A _Gk is the quadratic term coefficient of the quadratic function of the consumption characteristic of the adjustable thermal power unit that has been put into operation in the current period t;

B_Gk为当前时段t第Gk台已投入运行的可调节型火电机组的耗量特性二次函数的一次项系数；B _Gk is the first term coefficient of the quadratic function of the consumption characteristic of the adjustable thermal power unit that has been put into operation in the current period t;

子步骤205：如果

超限，则令

为该台已投入运行的可调节型火电机组的不越限时的最大出力值，且令nFlag=0；Sub-step 205: if

overrun, then order

is the maximum output value of the adjustable thermal power unit that has been put into operation without exceeding the limit, and nFlag=0;

子步骤206：判断Gk>N_G是否成立，如果Gk>N_G，则执行子步骤207；否则，令Gk=Gk+1，返回子步骤204；Sub-step 206: judge whether Gk>N _G is established, if Gk> _NG , then execute sub-step 207; otherwise, let Gk=Gk+1, return to sub-step 204;

N_G为当前时段t已投入运行的可调节型火电机组的数量；N _G is the number of adjustable thermal power units that have been put into operation in the current period t;

子步骤207：判断nFlag=0是否成立，如果nFlag=0，返回子步骤201；否则，执行子步骤208；Sub-step 207: judge whether nFlag=0 is established, if nFlag=0, return to sub-step 201; otherwise, execute sub-step 208;

子步骤208：令当前时段t的每个抽样场景下第Gk台已投入运行的可调节型火电机组提供的负荷为

Sub-step 208: let the load provided by the Gkth adjustable thermal power unit that has been put into operation under each sampling scenario of the current period t be

所述计算当前时段t的每个抽样场景下已投入运行的可调节型火电机组的运行成本采用公式 $F_{\cos t} (P_{G, t}^{m}) = Σ_{Gk = 1}^{N_{T}} [μ_{Gk, t}^{m} \times f_{k} (P_{Gk, t}^{m}) + μ_{Gk, t}^{m} (1 - μ_{Gk, t - 1}^{m}) S_{Gk}];$ The calculation of the operating cost of the adjustable thermal power unit that has been put into operation under each sampling scenario of the current period t adopts the formula $f_{\cos t} (P_{G, t}^{m}) = Σ_{K = 1}^{N_{T}} [μ_{K, t}^{m} \times f_{k} (P_{K, t}^{m}) + μ_{K, t}^{m} (1 - μ_{K, t - 1}^{m}) S_{K}];$

其中，

为当前时段t的第m个抽样场景下已投入运行的可调节型火电机组的运行成本，m=1，2，...，N，N为每个时段中的抽样场景数目；in,

is the operating cost of the adjustable thermal power unit that has been put into operation in the mth sampling scenario of the current period t, m=1, 2,..., N, N is the number of sampling scenarios in each period;

N_T为可调节型火电机组数量；N _T is the number of adjustable thermal power units;

为当前时段t第Gk台可调节型火电机组投入运行状态，当前时段t第Gk台可调节型火电机组投入运行，则

当前时段t第Gk台可调节型火电机组未投入运行，则

The Gkth adjustable thermal power unit is put into operation in the current period t, and the Gkth adjustable thermal power unit is put into operation in the current period t, then

The Gkth adjustable thermal power unit in the current period t is not put into operation, then

为当前时段t第Gk台可调节型火电机组的耗量特性二次函数，且

is the quadratic function of the consumption characteristic of the Gkth adjustable thermal power unit in the current period t, and

${f f}_{k k} (({P P}_{Gk K,, t t}^{m m})) = = {A A}_{Gk K} \times \times {(({P P}_{Gk K,, t t}^{m m}))}^{22} + + {B B}_{Gk K} \times \times {P P}_{Gk K,, t t}^{m m} + + {C C}_{Gk K};;$

A_Gk为当前时段t第Gk台可调节型火电机组的耗量特性二次函数的二次项系数；A _Gk is the quadratic term coefficient of the quadratic function of the consumption characteristic of the Gkth adjustable thermal power unit in the current period;

B_Gk为当前时段t第Gk台可调节型火电机组的耗量特性二次函数的一次项系数；B _Gk is the first term coefficient of the quadratic function of the consumption characteristic of the Gkth adjustable thermal power unit in the current period;

C_Gk为当前时段t第Gk台可调节型火电机组的耗量特性二次函数的常数项系数；C _Gk is the coefficient of the constant term of the quadratic function of the consumption characteristic of the Gkth adjustable thermal power unit in the current period;

S_Gk为当前时段t第Gk台可调节型火电机组的启动费用。S _Gk is the start-up cost of the Gkth adjustable thermal power unit in the current period.

所述当前时段t每个抽样场景下已投入运行的可调节型火电机组的负荷调节充裕度区间为 $[P_{L 3, t}^{m} - P_{H, t}^{ne} - P_{G, t}^{ne, m}, P_{L 3, t}^{m} + P_{H, t}^{op} + P_{G, t}^{op, m}];$ The load adjustment adequacy range of the adjustable thermal power units that have been put into operation under each sampling scenario in the current period t is $[P_{L 3, t}^{m} - P_{h, t}^{ne} - P_{G, t}^{ne, m}, P_{L 3, t}^{m} + P_{h, t}^{op} + P_{G, t}^{op, m}];$

其中，

为当前时段t可调节型水电机组向下调节极限值；in,

Adjust the limit value downward for the adjustable hydroelectric unit in the current period t;

为当前时段t每个抽样场景下可调节型火电机组运行下界，且

is the lower bound of the operation of adjustable thermal power units in each sampling scenario in the current period t, and

${P P}_{G G,, t t}^{ne ne,, m m} = = {Σ Σ}_{Gk K = = 11}^{{N N}_{T T}} {μ μ}_{Gk K,, t t}^{m m} {P P}_{Gk K,, t t}^{ne ne,, m m};;$

为当前时段t第Gk台可调节型火电机组的运行下界，且

is the lower bound of the operation of the Gkth adjustable thermal power unit in the current period t, and

${P P}_{Gk K,, t t}^{ne ne,, m m} = = min min (({P P}_{Gk K,, t t}^{m m} - - {P P}_{Gk K,, min min},, - - {R R}_{Gk K,, d d} \times \times Δt Δt));;$

P_Gk，min为第Gk台可调节型火电机组出力下限；P _{Gk, min} is the output lower limit of the Gkth adjustable thermal power unit;

R_Gk，d为第Gk台可调节型火电机组向下爬坡速率；R _{Gk, d} is the downward climbing rate of the Gkth adjustable thermal power unit;

当前时段t第Gk台可调节型火电机组未投入运行，则

Δt为当前时段t的时间长度；Δt is the time length of the current period t;

N_T为可调节型火电机组的数量；N _T is the number of adjustable thermal power units;

为当前时段t可调节型水电机组向上调节极限值； Adjust the limit value upward for the adjustable hydroelectric unit in the current period t;

为当前时段t每个抽样场景下可调节型火电机组运行上界，且

is the upper bound of the operation of adjustable thermal power units in each sampling scenario in the current period t, and

${P P}_{G G,, t t}^{op op,, m m} = = {Σ Σ}_{Gk K = = 11}^{{N N}_{T T}} {μ μ}_{Gk K,, t t}^{m m} {P P}_{Gk K,, t t}^{op op,, m m};;$

为当前时段t第Gk台可调节型火电机组运行上界，且 is the upper bound of the operation of the Gkth adjustable thermal power unit in the current period t, and

${P P}_{Gk K,, t t}^{op op,, m m} = = min min (({P P}_{Gk K,, max max} - - {P P}_{Gk K,, t t}^{m m},, - - {R R}_{Gk K,, u u} \times \times Δt Δt));;$

P_Gk，max为第Gk台可调节型火电机组的出力上限；P _Gk,max is the output upper limit of the Gkth adjustable thermal power unit;

R_Gk，u为第Gk台可调节型火电机组向上爬坡速率；R _{Gk, u} is the upward climbing rate of the Gkth adjustable thermal power unit;

m=1，2，...，N，N为每个时段中的抽样场景数目。m=1, 2, . . . , N, where N is the number of sampled scenes in each period.

所述按照第二设定顺序调整当前时段已投入运行且尚有可调节容量的可调节型发电机组具体为，先调整当前时段已投入运行且尚有可调节容量的可调节型水电机组，当前时段已投入运行的所有可调节型水电机组没有可调节容量时，按输出功率调节速率由高到低的顺序调整已投入运行的燃气轮机机组，当前时段已投入运行的所有燃气轮机机组没有可调节容量时，按最小比耗量由低到高的顺序调整投入运行的可停机的可调节型火电机组，当前时段已投入运行的所有可停机的可调节型火电机组没有可调节容量时，按最小运行比耗量由低到高的顺序调整不可停机的可调节型火电机组。The adjustment according to the second setting sequence of the adjustable generating units that have been put into operation in the current period and still have adjustable capacity is specifically, firstly adjust the adjustable hydroelectric units that have been put into operation in the current period and still have adjustable capacity, and the current When all the adjustable hydropower units that have been put into operation during the period have no adjustable capacity, adjust the gas turbine units that have been put into operation according to the order of output power adjustment rate from high to low. When all the gas turbine units that have been put into operation in the current period have no adjustable capacity , adjust the adjustable thermal power units that are put into operation and can be stopped according to the order of minimum specific consumption from low to high. The adjustable thermal power units that cannot be shut down are adjusted in order of consumption from low to high.

本发明提供的方法，通过研究在数日或一周的调度周期内系统机组组合与调度问题，使系统内各类型电源的出力能够满足系统用电负荷的需求，保证系统功率平衡与稳定运行。The method provided by the invention enables the output of various types of power sources in the system to meet the demand of the system's power load by studying the combination and scheduling of system units within a scheduling period of several days or a week, ensuring system power balance and stable operation.

附图说明Description of drawings

图1是本发明提供的一种多能源联合发电系统运行模拟的优化方法流程图；Fig. 1 is a flow chart of an optimization method for operation simulation of a multi-energy combined power generation system provided by the present invention;

图2是各个时段每个抽样场景下风力发电机组出力示意图；Figure 2 is a schematic diagram of the output of wind turbines under each sampling scenario in each time period;

图3是可调节型水电机组工作位置示意图；Fig. 3 is a schematic diagram of the working position of the adjustable hydroelectric unit;

图4是各个时段每个抽样场景下已投入运行的可调节型火电机组分配负荷流程图。Figure 4 is a flow chart of the load distribution of adjustable thermal power units that have been put into operation under each sampling scenario in each period.

具体实施方式Detailed ways

下面结合附图，对优选实施例作详细说明。应该强调的是，下述说明仅仅是示例性的，而不是为了限制本发明的范围及其应用。The preferred embodiments will be described in detail below in conjunction with the accompanying drawings. It should be emphasized that the following description is only exemplary and not intended to limit the scope of the invention and its application.

图1是本发明提供的一种多能源联合发电系统运行模拟的优化方法流程图。如图1所示，本发明提供的多能源联合发电系统运行模拟的优化方法包括Fig. 1 is a flowchart of an optimization method for operation simulation of a multi-energy combined power generation system provided by the present invention. As shown in Figure 1, the optimization method of the operation simulation of the multi-energy combined power generation system provided by the present invention includes

如下步骤：Follow the steps below:

步骤1：设定调度周期并将所述调度周期等分成T个时段。Step 1: Set a scheduling cycle and divide the scheduling cycle into T time periods.

其中，m=1，2，...，N，j=1，2，...，T。Step 2: Set the number N of sampling scenarios in each period, and determine the combined output value of wind turbines and photovoltaic units in each sampling scenario in each period

Wherein, m=1, 2, . . . , N, j=1, 2, . . . , T.

风力发电和光伏发电都是新能源发电，通常根据风力发电机组和光伏发电机组预测出力值及其预测偏差值（即实际出力值与预测出力之差）的统计数据，在每个时段内生成若干出力场景，以反映风力、光伏等新能源发电机组的出力的间歇与波动性。因此，确定各个时段每个抽样场景的风电机组和光电机组混合出力值

包括如下子步骤：Both wind power generation and photovoltaic power generation are new energy power generation. Usually, according to the statistical data of the predicted output value of the wind power generator set and the photovoltaic generator set and its predicted deviation value (that is, the difference between the actual output value and the predicted output value), several Output scenarios to reflect the intermittent and fluctuating output of new energy generating units such as wind power and photovoltaics. Therefore, determine the mixed output value of wind turbines and photovoltaic units for each sampling scene in each period

Including the following sub-steps:

子步骤101：分别根据风电机组和光电机组出力的历史数据，做出风电机组出力预测偏差的频率分布图和光电机组出力预测偏差的频率分布图。Sub-step 101: According to the historical data of the output of the wind turbine and the photovoltaic unit, respectively, make a frequency distribution diagram of the output prediction deviation of the wind turbine and a frequency distribution diagram of the deviation of the photovoltaic unit output prediction.

子步骤102：通过拟合方法分别得到各个时段风电机组出力预测偏差概率密度函数和光电机组出力预测偏差概率密度函数。Sub-step 102: The probability density function of the output prediction deviation of the wind power unit and the probability density function of the output prediction deviation of the photovoltaic unit in each period are respectively obtained by a fitting method.

子步骤103：对各个时段风电机组出力预测偏差概率密度函数和光电机组出力预测偏差概率密度函数分别进行积分运算，得到各个时段的风电机组出力预测偏差概率密度分布函数和光电机组出力预测偏差概率密度分布函数。Sub-step 103: Integrate the wind turbine output forecast deviation probability density function and the photoelectric unit output forecast deviation probability density function in each period to obtain the wind turbine output forecast deviation probability density function and photovoltaic unit output forecast deviation probability density in each period Distribution function.

上述子步骤101-103是本领域技术人员常用的发电机组出力预测偏差计算方法。由于在各个时段每个抽样场景下只能有一个预测值，因此上述两个概率密度分布函数都是一一对应的可逆函数。The above sub-steps 101-103 are commonly used methods for calculating the output prediction deviation of generator sets by those skilled in the art. Since there can only be one predicted value for each sampling scenario in each time period, the above two probability density distribution functions are reversible functions corresponding to each other.

子步骤104：计算各个时段每个抽样场景下的基础采样值，其计算公式如下：Sub-step 104: Calculate the basic sampling value under each sampling scene in each time period, the calculation formula is as follows:

${U u}_{m m} = = \frac{U u}{N N} + + \frac{m m - - 11}{N N} - - - - - - ((11))$

公式（1）中，m=1，2，...，N，N为每个时段抽样场景数目，U为区间(0，1)上的随机数。In formula (1), m=1, 2, ..., N, N is the number of sampling scenes in each period, and U is a random number on the interval (0, 1).

子步骤105：计算各个时段每个抽样场景下风电机组出力预测偏差和光电机组出力预测偏差。其中，各个时段每个抽样场景下风电机组出力预测偏差的计算公式为：Sub-step 105: Calculate the output prediction deviation of the wind power unit and the output prediction deviation of the photovoltaic unit under each sampling scenario in each time period. Among them, the calculation formula of wind turbine output forecast deviation in each sampling scenario in each period is:

${X x}_{W W,, j j}^{m m} = = {F f}_{W W,, j j}^{- - 11} (({U u}_{m m})) - - - - - - ((22))$

各个时段每个抽样场景下光电机组出力预测偏差的计算公式为：The formula for calculating the output prediction deviation of photoelectric units under each sampling scenario in each period is:

${X x}_{P P,, j j}^{m m} = = {F f}_{P P,, j j}^{- - 11} (({U u}_{m m})) - - - - - - ((33))$

公式（2）中，为风电机组出力预测偏差概率密度分布函数的反函数。公式（3）中，

为光电机组出力预测偏差概率密度分布函数的反函数。j=1，2，...，T，T为调度周期的时段数。In formula (2), It is the inverse function of the probability density distribution function of wind turbine output prediction deviation. In formula (3),

It is the inverse function of the output prediction deviation probability density distribution function of the photoelectric unit. j=1, 2, ..., T, T is the number of time slots in the scheduling cycle.

子步骤106：计算各个时段每个抽样场景的风电机组和光电机组混合出力值，该混合出力值的计算公式为：Sub-step 106: Calculate the mixed output value of wind turbines and photovoltaic units for each sampling scene in each time period. The calculation formula of the mixed output value is:

${P P}_{New new,, j j}^{m m} = = {P P}_{W W,, j j} + + {X x}_{W W,, j j}^{m m} + + {P P}_{P P,, j j} + + {X x}_{P P,, j j}^{m m} - - - - - - ((44))$

公式（4）中，P_W，j为各个时段风电机组出力预测值，P_P，j为各个时段光电机组出力预测值。图2是各个时段每个抽样场景下风力发电机组出力示意图。各个时段每个抽样场景下风力发电机组出力值应为各个时段每个抽样场景下光电机组出力值应为

In formula (4), P _W,j is the predicted value of wind turbine output in each period, and P _P,j is the predicted output of photovoltaic unit in each period. Figure 2 is a schematic diagram of the output of wind turbines in each sampling scenario in each time period. The output value of wind turbines in each sampling scenario in each time period should be The output value of the photoelectric unit in each sampling scene in each period should be

步骤3：计算每个时段系统中可调节型发电机组应承担的负荷P_L1，j，并将P_L1，j记为第一等值负荷。Step 3: Calculate the load P _L1,j that the adjustable generator set in the system should bear in each time period, and record P _L1,j as the first equivalent load.

在多能源联合发电系统中，发电机组通常分为不可调节型发电机组和可调节型发电机组。不可调节型发电机组的出力在整个发电过程中都是无法变动的稳定出力。不可调节型发电机组的出力包括无调节型水电机组的出力、可调节型水电机组的强迫出力、热电机组的强迫出力和核电机组的出力，这些出力由于是系统运行的平稳出力，工作位置在负荷曲线基部，其出力不参与功率平衡的调节过程。而系统中的可调节型发电机组则包括可调节型水电机组和可调节型火电机组。In a multi-energy combined power generation system, generator sets are usually divided into non-adjustable generator sets and adjustable generator sets. The output of the non-adjustable generating set is a stable output that cannot be changed during the entire power generation process. The output of non-adjustable generating units includes the output of non-adjustable hydroelectric units, the forced output of adjustable hydroelectric units, the forced output of thermoelectric units and the output of nuclear power units. Since these outputs are the stable output of the system operation, the working position is under load. At the base of the curve, its output does not participate in the adjustment process of power balance. The adjustable generating units in the system include adjustable hydroelectric units and adjustable thermal power units.

因此，每个时段系统中可调节型发电机组应承担的负荷P_L1，j的计算公式为：Therefore, the calculation formula of the load P _L1,j that the adjustable generator set should bear in the system at each time period is:

P_L1，j=P_L，j-P_f，j （5）P _{L1, j} = P _{L, j} - P _{f, j} (5)

公式（5）中，P_L，j为第j个时段系统应承担的总负荷，P_f，j为第j个时段系统中除可调节型发电机组之外的发电机组提供的负荷。In the formula (5), _{PL, j} is the total load that the system should bear in the jth period, and P _{f, j} is the load provided by the generators in the system except the adjustable generators in the jth period.

步骤4：计算每个时段系统中可调节型火电机组应承担的负荷P_L2，j，并将P_L2，j记为第二等值负荷。Step 4: Calculate the load P _{L2,j that the adjustable thermal power unit in the system should bear in each time period,} and record P _L2,j as the second equivalent load.

由于P_L1，j是每个时段系统中可调节型发电机组应承担的负荷，因此每个时段系统中可调节型火电机组应承担的负荷，也就是每个时段系统中可调节型发电机组应承担的负荷减去每个时段系统中可调节型水电机组应承担的负荷。Since P _{L1, j} is the load that the adjustable generator set in the system should bear in each period, so the load that the adjustable thermal power unit in the system should bear in each period, that is, the adjustable generator set in the system in each period should be The load borne minus the load that the adjustable hydroelectric unit in the system should bear in each period.

可调节型水电机组和抽水蓄能发电机组在负荷高峰时期发电。抽水蓄能发电机组安排在调度周期内每日的低谷负荷时段进行抽水，在发电工况下则视为可调节型水电机组。根据可调节型水电机组的水文条件以及每日的计划发电量，确定可调节型水电机组在负荷曲线上的工作位置，如图3所示。在图3中，P_H，up为相应时段可调节型水电机组在负荷曲线上的工作位置上限，P_H，down为相应时段可调节型水电机组在负荷曲线上的工作位置下限，则实际系统中可调节型水电机组的预期出力值为： $P_{H, j} = \{\begin{matrix} P_{H, up} - P_{H, down}, & P_{L 1, j} &GreaterEqual; P_{H, up} \\ P_{L 1, j} - P_{H, down}, & P_{H, down} < P_{L 1, j} < P_{H, up} \\ 0, & P_{L 1, j} \leq P_{H, down} \end{matrix},$ 也即第j个时段的可调节型水电机组的预期出力值。Adjustable hydroelectric units and pumped storage generating units generate electricity during peak load periods. The pumped storage power generation unit is arranged to pump water during the daily low load period in the dispatch cycle, and it is regarded as an adjustable hydroelectric unit under the power generation condition. According to the hydrological conditions of the adjustable hydroelectric unit and the daily planned power generation, determine the working position of the adjustable hydroelectric unit on the load curve, as shown in Figure 3. In Fig. 3, P _H,up is the upper limit of the working position of the adjustable hydroelectric unit on the load curve of the corresponding period, and P _H,down is the lower limit of the working position of the adjustable hydroelectric unit on the load curve of the corresponding period, then the actual system The expected output value of the adjustable hydroelectric unit is: $P_{h, j} = \{\begin{matrix} P_{h, up} - P_{h, down}, & P_{L 1, j} &Greater Equal; P_{h, up} \\ P_{L 1, j} - P_{h, down}, & P_{h, down} < P_{L 1, j} < P_{h, up} \\ 0, & P_{L 1, j} \leq P_{h, down} \end{matrix},$ That is, the expected output value of the adjustable hydroelectric unit in the jth period.

调度周期内可调节型水电机组的预期发电量，是根据电站水文条件、库容调节计划确定的定值。因此，可以利用从负荷曲线的日负荷峰值处定步长下移的方式，确定可调节型水电机组的工作上下限的位置，保证上下限功率的差值尽量接近可调节型水电机组的预想出力。The expected power generation of the adjustable hydroelectric unit in the dispatch period is a fixed value determined according to the hydrological conditions of the power station and the storage capacity adjustment plan. Therefore, the position of the upper and lower limits of the adjustable hydroelectric unit can be determined by moving down from the daily load peak of the load curve to ensure that the difference between the upper and lower limit power is as close as possible to the expected output of the adjustable hydroelectric unit .

那么，对于每个时段系统中可调节型火电机组应承担的负荷P_L2，j，其计算公式应为：Then, for the load P _L2,j that the adjustable thermal power unit in the system should bear in each time period, the calculation formula should be:

${P P}_{L L 22,, j j} = = {P P}_{L L 11,, j j} - - {Σ Σ}_{Hi hi = = 11}^{{N N}_{H h}} {μ μ}_{Hi hi,, j j} {P P}_{Hi hi,, j j} + + {μ μ}_{Hp HP,, j j} {P P}_{Hp HP,, j j} - - - - - - ((66))$

公式（6）中，μ_Hi，j为第j个时段第Hi个可调节型水电机组的工作状态，当第j个时段第Hi个可调节型水电机组投入运行时，μ_Hi，j=1；当第j个时段第Hi个可调节型水电机组未投入运行时，μ_Hi，j=0。N_H为可调节型水电机组的数量。μ_Hp，j为第j个时段抽水蓄能发电机组的抽水状态，当第j个时段抽水蓄能发电机组投入运行时，μ_Hp，j=1，当第j个时段抽水蓄能发电机组未投入运行时，μ_Hp，j=0。P_Hp，j为第j个时段抽水蓄能发电机组的抽水功率。P_Hi，j为第j个时段第Hi个可调节型水电机组的输出功率，其计算公式为：In formula (6), μHi _,j is the working state of the Hi adjustable hydroelectric unit in the jth period, when the Hi adjustable hydroelectric unit is put into operation in the jth period, μHi _,j =1 ; When the Hi-th adjustable hydroelectric unit is not in operation in the j-th period, μ _Hi,j =0. N _H is the number of adjustable hydroelectric units. μ _Hp,j is the pumping state of the pumped storage generator set in the jth period. When the pumped storage generator set is put into operation in the jth period, μ _{Hp, j} = 1. When put into operation, μ _{Hp, j} =0. P _Hp,j is the pumping power of the pumped storage generator set in the jth period. P _{Hi, j} is the output power of the Hi-th adjustable hydroelectric unit in the j-th period, and its calculation formula is:

P_Hi，j=η_Hi×g×(H_up，Hi，j-H_{down，Hi，j})×Q_Hi，j （7）P _{Hi, j} = η _Hi × g × (H _{up, Hi, j} - H _{down, Hi, j} ) × Q _{Hi, j} (7)

公式（7）中，η_Hi为第Hi个可调节型水电机组的发电效率系数，g为重力常数且取g=9.81，单位：牛顿/千克，H_up，Hi，j为第j个时段第Hi个可调节型水电机组的坝上水位，H_{down，Hi，j}为第j个时段第Hi个可调节型水电机组的下游水位，Q_Hi，j为第j个时段第Hi个可调节型水电机组的发电平均流量。In formula (7), _ηHi is the power generation efficiency coefficient of the Hi-th adjustable hydroelectric unit, g is the gravitational constant and takes g=9.81, unit: Newton/kg, H _{up, Hi, j} is the j-th time interval The water level above the dam of the Hi adjustable hydroelectric unit, H _down,Hi,j is the downstream water level of the Hi adjustable hydroelectric unit in the jth period, Q _Hi,j is the Hi adjustable hydroelectric unit in the jth period The average flow of electricity generated by the hydroelectric unit.

与此同时，分别计算每个时段可调节型水电机组向上调节极限值

和向下调节极限值

上述两个值在后续步骤的计算中会用到。At the same time, calculate the upward adjustment limit value of the adjustable hydroelectric unit in each period

and down adjustment limit

The above two values will be used in the calculations in the subsequent steps.

计算每个时段可调节型水电机组向上调节极限值

采用的公式为：Calculate the upward adjustment limit value of the adjustable hydroelectric unit in each period

The formula used is:

${P P}_{H h,, j j}^{op op} = = {Σ Σ}_{Hi hi = = 11}^{{N N}_{H h}} {μ μ}_{Hi hi,, j j} {P P}_{Hi hi,, j j}^{op op} - - - - - - ((88))$

公式（8）中，

为第j个时段第Hi个可调节型水电机组的向上调节容量且

为第j个时段第Hi个可调节型水电机组的最大输出功率，P_Hi，j为第j个时段第Hi个可调节型水电机组的输出功率。μ_Hi，j为第j个时段第Hi个可调节型水电机组的工作状态，当第j个时段第Hi个可调节型水电机组投入运行时，μ_Hi，j=1；当第j个时段第Hi个可调节型水电机组未投入运行时，μ_Hi，j=0。In formula (8),

is the upward adjustment capacity of the Hi-th adjustable hydroelectric unit in the j-th period and

P _{Hi, j} is the output power of the Hi adjustable hydroelectric unit in the jth period. μHi _,j is the working state of the Hi adjustable hydroelectric unit in the jth period, when the Hi adjustable hydroelectric unit is put into operation in the jth period, μHi _,j =1; when the jth period When the Hi-th adjustable hydroelectric unit is not in operation, μHi _,j =0.

计算每个时段可调节型水电机组向下调节极限值

采用的公式为：Calculate the downward adjustment limit value of the adjustable hydroelectric unit in each period

The formula used is:

${P P}_{H h,, j j}^{ne ne} = = {Σ Σ}_{Hi hi = = 11}^{{N N}_{H h}} {μ μ}_{Hi hi,, j j} {P P}_{Hi hi,, j j}^{ne ne} - - - - - - ((99))$

公式（9）中，

为第j个时段第Hi个可调节型水电机组的向下调节容量且为第j个时段第Hi个可调节型水电机组的最小强迫输出功率，P_Hi，j为第j个时段第Hi个可调节型水电机组的输出功率。μ_Hi，j为第j个时段第Hi个可调节型水电机组的工作状态，当第j个时段第Hi个可调节型水电机组投入运行时，μ_Hi，j=1；当第j个时段第Hi个可调节型水电机组未投入运行时，μ_Hi，j=0。In formula (9),

is the downward adjustment capacity of the Hi-th adjustable hydroelectric unit in the j-th period and P _{Hi, j} is the output power of the Hi adjustable hydroelectric unit in the jth period. μHi _,j is the working state of the Hi adjustable hydroelectric unit in the jth period, when the Hi adjustable hydroelectric unit is put into operation in the jth period, μHi _,j =1; when the jth period When the Hi-th adjustable hydroelectric unit is not in operation, μHi _,j =0.

步骤5：确定可调节型火电机组的初始投入运行状态。Step 5: Determine the initial operation status of the adjustable thermal power unit.

在混合发电系统中，可调节型火电机组一般包括不可停机的可调节型火电机组、可停机的可调节型火电机组和燃气轮机机组。不可停机的可调节型火电机组在调度周期内不进行启停调峰（不考虑故障检修特例）。可停机的可调节型火电机组可以受调度命令安排启停机。燃气轮机机组容量小、响应速度快，做调节电源使用。In a hybrid power generation system, adjustable thermal power units generally include adjustable thermal power units that cannot be shut down, adjustable thermal power units that can be shut down, and gas turbine units. Adjustable thermal power units that cannot be shut down do not perform peak shaving during the dispatch period (the special case of troubleshooting is not considered). The adjustable thermal power unit that can be stopped can be scheduled to start and stop according to the dispatch command. The gas turbine unit has small capacity and fast response speed, so it is used for regulating power supply.

考虑到参与启停调峰的小型火电和燃气机组，有最小启停时间的限制，因此考虑当前时段的机组组合方案时，需要联合之后的预测负荷情况，以便让机组有一定缓冲时间进行启动或关停，本方法选用包括当前时段在内的向后连续5个时段。所以，在确定可调节型火电机组的初始投入运行状态时，先分别计算调度周期前5个时段的第二等值负荷最大值和前5个时段的第二等值负荷最小值，然后按照第一设定顺序投入可调节型火电机组。Considering that the small thermal power and gas-fired units involved in start-stop peak regulation have a minimum start-stop time limit, so when considering the unit combination plan for the current period, it is necessary to combine the predicted load conditions so that the units have a certain buffer time for start-up or Shutdown, this method selects 5 consecutive backward periods including the current period. Therefore, when determining the initial operating state of the adjustable thermal power unit, first calculate the second equivalent load maximum value and the first five periods of the first five periods of the second equivalent load minimum value of the dispatch cycle, and then according to the first A set sequence is put into the adjustable thermal power unit.

此处，按照第一设定顺序投入可调节型火电机组具体为，先按最小运行比耗量由低到高的顺序投入不可停机的可调节型火电机组，再按最小比耗量由低到高的顺序投入可停机的可调节型火电机组，最后按输出功率调节速率由高到低的顺序投入燃气轮机机组。陆续投入可调节型火电机组，当投入运行的可调节型火电机组的总装机容量大于调度周期前5个时段的第二等值负荷最大值，并且投入运行的可调节型火电机组的最小技术出力总和小于调度周期前5个时段的第二等值负荷最小值时，停止投入可调节型火电机组，可调节型火电机组的初始投入运行状态确定完毕。Here, according to the first setting sequence, the adjustable thermal power units are put into operation specifically, firstly, the adjustable thermal power units that cannot be shut down are put into operation according to the order of the minimum operating specific consumption from low to high, and then according to the minimum specific consumption from low to high. The adjustable thermal power units that can be shut down are put into operation in the highest order, and the gas turbine units are put into operation according to the order of output power adjustment rate from high to low. When adjustable thermal power units are put into operation one after another, when the total installed capacity of the adjustable thermal power units put into operation is greater than the second equivalent load maximum value in the first 5 periods of the dispatch cycle, and the minimum technical output of the adjustable thermal power units put into operation When the sum is less than the minimum value of the second equivalent load in the first 5 periods of the scheduling cycle, the adjustable thermal power unit is stopped, and the initial operation status of the adjustable thermal power unit is determined.

上述，可调节型火电机组的最小技术出力是指，保证锅炉维持稳定燃烧（不至于停火或偏离稳态的温度与压强）的情况下，机组可以发出的最小功率值。每台火电机组的最小技术出力都由机组厂家制造完成后提供。As mentioned above, the minimum technical output of the adjustable thermal power unit refers to the minimum power value that the unit can generate under the condition of ensuring that the boiler maintains stable combustion (without stopping the fire or deviating from the steady-state temperature and pressure). The minimum technical output of each thermal power unit is provided by the unit manufacturer after it is manufactured.

步骤6：设定初始时段t=5。Step 6: Set the initial time period t=5.

步骤7：计算当前时段t可调节型火电机组在每个抽样场景下应承担的负荷

并将

记为第三等值负荷。第三等值负荷

的计算公式为：Step 7: Calculate the load that the current period t adjustable thermal power unit should bear in each sampling scenario

and will

Recorded as the third equivalent load. third equivalent load

The calculation formula is:

${P P}_{L L 33,, t t}^{m m} = = {P P}_{L L 22,, t t} - - {P P}_{New new,, t t}^{m m} - - - - - - ((1010))$

同时确定当前时段t每个抽样场景下已投入运行的可调节型火电机组的负荷调节充裕度区间。Step 8: Assign loads to the adjustable thermal power units that have been put into operation in each sampling scenario of the current period t

At the same time, determine the load regulation adequacy interval of the adjustable thermal power units that have been put into operation under each sampling scenario in the current period t.

第三等值负荷

是当前时段t可调节型火电机组在每个抽样场景下应承担的总负荷，具体到每台可调节型火电机组，在当前时段t每个抽样场景下应承担的多少负荷，应当通过下述方法计算。third equivalent load

is the total load that the adjustable thermal power unit should bear in each sampling scenario in the current period t. Specifically, for each adjustable thermal power unit, how much load it should bear in each sampling scenario in the current period t should be determined by the following method calculation.

首先，每台火电机组都有各自的耗量特性二次函数，该耗量特性二次函数有时火电机组的生产厂家会直接给出。通过对火电机组耗量特性曲线拟合也可得到。我们设定第Gk台火电机组的耗量特性二次函数为

其中，A_k为第k台火电机组的耗量特性二次函数的二次项系数，B_k为第k台火电机组的耗量特性二次函数的一次项系数，C_k为第k台火电机组的耗量特性二次函数的常数项系数，P_k为第k台火电机组的负荷。First of all, each thermal power unit has its own quadratic function of consumption characteristics, which is sometimes directly given by the manufacturer of the thermal power unit. It can also be obtained by fitting the consumption characteristic curve of thermal power units. We set the consumption characteristic quadratic function of the Gkth thermal power unit as

Among them, A _k is the coefficient of the quadratic term of the quadratic function of the consumption characteristic of the kth thermal power unit, B _k is the coefficient of the first term of the quadratic function of the consumption characteristic of the kth thermal power unit, and C _k is the coefficient of the quadratic term of the kth thermal power unit The coefficient of the constant term of the quadratic function of the consumption characteristic of the unit, P _k is the load of the kth thermal power unit.

当火电机组在出力变化时，若耗量微增率均等，则响应负荷变化的调节效果最优，用数学方程组描述，则为When the output of the thermal power unit changes, if the consumption micro-increase rate is equal, then the adjustment effect in response to the load change is optimal, which is described by a set of mathematical equations, which is

$\{\begin{matrix} {f f}_{op op} = = {22 A A}_{11} {P P}_{11} + + {B B}_{11} = = {22 A A}_{22} {P P}_{22} + + {B B}_{22} = = \cdot &Center Dot; \cdot &Center Dot; \cdot &Center Dot; {22 A A}_{M m} {P P}_{M m} + + {B B}_{M m} \\ {Σ Σ}_{k k = = 11}^{M m} {P P}_{k k} = = {P P}_{L L} \end{matrix} - - - - - - ((1111))$

公式（11）中，M为火电机组的数量，P_L为所有火电机组的总负荷，f_op为耗量微增率。依照公式（11），可以解得：In formula (11), M is the number of thermal power units, _PL is the total load of all thermal power units, and f _op is the consumption micro-increase rate. According to formula (11), it can be solved as follows:

${f f}_{op op} = = \frac{{P P}_{L L} + + {Σ Σ}_{k k = = 11}^{M m} (({0.5 0.5 B B}_{k k} / / {A A}_{k k}))}{{Σ Σ}_{k k = = 11}^{M m} {(({22 A A}_{k k}))}^{- - 11}} - - - - - - ((1212))$

每台可调节型火电机组在当前时段每个抽样场景下应承担的负荷，应为不越限的值。因此，本发明采用循环迭代的方式，通过已投入运行的可调节型火电机组耗量微增率，调节已投入运行的可调节型火电机组的出力值，从而确保每台可调节型火电机组的出力均不越限。因此，根据公式（12），每台可调节型火电机组在当前时段t每个抽样场景下应承担的负荷的具体计算过程如图4所示，包括：The load that each adjustable thermal power unit should bear in each sampling scenario during the current period should be a value that does not exceed the limit. Therefore, the present invention adopts a cyclic and iterative method to adjust the output value of the adjustable thermal power unit that has been put into operation through the slight increase rate of the consumption of the adjustable thermal power unit that has been put into operation, thereby ensuring the output of each adjustable thermal power unit. The output does not exceed the limit. Therefore, according to formula (12), the specific calculation process of the load that each adjustable thermal power unit should bear in each sampling scenario in the current period t is shown in Figure 4, including:

子步骤201：设定已投入运行的可调节型火电机组出力越限标志nFlag=1。Sub-step 201: Set the output limit flag nFlag=1 of the adjustable thermal power unit that has been put into operation.

子步骤202：计算已投入运行的可调节型火电机组耗量微增率。开始时，可能会有已投入运行的可调节型火电机组的出力值处于越限的状态。计算耗量微增率时，应当只考虑那些出力值不越限的已投入运行的可调节型火电机组。Sub-step 202: Calculate the consumption micro-increase rate of the adjustable thermal power unit that has been put into operation. At the beginning, the output value of the adjustable thermal power unit that has been put into operation may be in a state of exceeding the limit. When calculating the slight increase rate of consumption, only those adjustable thermal power units that have been put into operation and whose output value does not exceed the limit should be considered.

采用公式（13）计算已投入运行的可调节型火电机组耗量微增率：Use formula (13) to calculate the slight increase rate of consumption of adjustable thermal power units that have been put into operation:

${f f}_{op op} = = \frac{{P P}_{L L 33,, t t}^{m m} + + {Σ Σ}_{{Gk K}^{' '} = = 11}^{{N N}_{G G}^{' '}} ((0.5 0.5 {B B}_{Gk K}' ' / / {A A}_{G G {k k}^{' '}}))}{{Σ Σ}_{{Gk K}^{' '} = = 11}^{{N N}_{G G}^{' '}} {(({22 A A}_{{Gk K}^{' '}}))}^{- - 11}} - - - - - - ((1313))$

公式（13）中，

为第三等值负荷且m=1，2，...，N，N为每个时段抽样场景数目，N′_G为当前时段t已投入运行且出力未超限的可调节型火电机组的数量，A_Gk′为当前时段t第Gk′台已投入运行且出力未超限的可调节型火电机组的耗量特性二次函数的二次项系数，B_Gk′为当前时段t第Gk′台已投入运行且出力未超限的可调节型火电机组的耗量特性二次函数的一次项系数。In formula (13),

is the third equivalent load and m=1, 2,..., N, N is the number of sampling scenes in each period, N′ _G is the adjustable thermal power unit that has been put into operation in the current period t and the output has not exceeded the limit A _Gk′ is the quadratic term coefficient of the quadratic function of the consumption characteristic of the adjustable thermal power unit whose Gk′th unit has been put into operation in the current period t and its output is not exceeded, and B _Gk′ is the Gk′th unit in the current period t The coefficient of the first term of the quadratic function of the consumption characteristic of an adjustable thermal power unit that has been put into operation and the output has not exceeded the limit.

子步骤203：令Gk=1。Sub-step 203: Let Gk=1.

子步骤204：利用耗量微增率f_op，更新每一个已投入运行的可调节型火电机组出力值，其公式为：Sub-step 204: Utilize the consumption micro-increase rate f _op to update the output value of each adjustable thermal power unit that has been put into operation, and the formula is:

${P P}_{Gk K,, t t}^{m m} = = \frac{{f f}_{op op} - - {B B}_{Gk K}}{{22 A A}_{Gk K}} - - - - - - ((1414))$

公式（14）中，A_Gk为当前时段t第Gk台已投入运行的可调节型火电机组的耗量特性二次函数的二次项系数，B_Gk为当前时段t第Gk台已投入运行的可调节型火电机组的耗量特性二次函数的一次项系数。In formula (14), A _Gk is the quadratic term coefficient of the quadratic function of the consumption characteristic of the adjustable thermal power unit that has been put into operation at the Gkth unit in the current period t, and B _Gk is the coefficient of the quadratic term of the quadratic function of the consumption characteristic of the Gkth unit that has been put into operation in the current period t The coefficient of the first term of the quadratic function of the consumption characteristic of the adjustable thermal power unit.

子步骤205：如果

超限，则令

为该台已投入运行的可调节型火电机组的不越限时的最大出力值，且令nFlag=0。Sub-step 205: if

overrun, then order

is the maximum output value of the adjustable thermal power unit that has been put into operation without exceeding the limit, and nFlag=0.

子步骤206：判断Gk>N_G是否成立，如果Gk>N_G，则执行子步骤207；否则，令Gk=Gk+1，返回子步骤204。N_G为当前时段t已投入运行的可调节型火电机组的数量。Sub-step 206: Determine whether Gk>N _G holds, if Gk> _NG , execute sub-step 207; otherwise, set Gk=Gk+1, return to sub-step 204. N _G is the number of adjustable thermal power units that have been put into operation in the current period t.

在所有已投入运行的可调节型火电机组的出力值更新完成后，可能还会有超限的情况存在，此时将越限标志nFlag置为0，通过子步骤207进行进一步的迭代。After the output values of all the adjustable thermal power units that have been put into operation are updated, there may still be an over-limit situation. At this time, the over-limit flag nFlag is set to 0, and further iterations are performed through sub-step 207 .

子步骤207：判断nFlag=0是否成立，如果nFlag=0，返回子步骤201；否则，执行子步骤208。Sub-step 207: judge whether nFlag=0 is true, if nFlag=0, return to sub-step 201; otherwise, execute sub-step 208.

如果在某次更新过程中，越限标志nFlag为1，则说明所有已投入运行的可调节型火电机组的出力值都不越限，此时当前时段t的每个抽样场景下第Gk台已投入运行的可调节型火电机组提供的负荷为即为其出力值

If in a certain update process, the over-limit flag nFlag is 1, it means that the output values of all adjustable thermal power units that have been put into operation do not exceed the limit. The load provided by the adjustable thermal power unit put into operation is its output value

步骤9：计算当前时段t的每个抽样场景下已投入运行的可调节型火电机组的运行成本，其计算公式为：Step 9: Calculate the operating cost of the adjustable thermal power units that have been put into operation under each sampling scenario in the current period t, and the calculation formula is:

${F f}_{cos cos t t} (({P P}_{G G,, t t}^{m m})) = = {Σ Σ}_{Gk K = = 11}^{{N N}_{T T}} [[{μ μ}_{Gk K,, t t}^{m m} \times \times {f f}_{k k} (({P P}_{Gk K,, t t}^{m m})) + + {μ μ}_{Gk K,, t t}^{m m} ((11 - - {μ μ}_{Gk K,, t t - - 11}^{m m})) {S S}_{Gk K}]] - - - - - - ((1515))$

公式（15）中，为当前时段t的第m个抽样场景下已投入运行的可调节型火电机组的运行成本，m=1，2，...，N，N为每个时段中的抽样场景数目，N_T为可调节型火电机组数量。

当前时段t第Gk台可调节型火电机组未投入运行，则

为当前时段t第Gk台可调节型火电机组的耗量特性二次函数的二次项系数，B_Gk为当前时段t第Gk台可调节型火电机组的耗量特性二次函数的一次项系数，C_Gk为当前时段t第Gk台可调节型火电机组的耗量特性二次函数的常数项系数，S_Gk为当前时段t第Gk台可调节型火电机组的启动费用。In formula (15), is the operating cost of the adjustable thermal power unit that has been put into operation in the mth sampling scenario of the current period t, m=1, 2,..., N, N is the number of sampling scenarios in each period, _NT is The number of adjustable thermal power units.

B Gk is the quadratic coefficient of the quadratic function of the consumption characteristic of the Gkth adjustable thermal power unit in the current period, and B _Gk is the primary coefficient of the quadratic function of the consumption characteristic of the Gkth adjustable thermal power unit in the current period , C _Gk is the constant term coefficient of the quadratic function of the consumption characteristic of the Gkth adjustable thermal power unit in the current period t, and S _Gk is the start-up cost of the Gkth adjustable thermal power unit in the current period t.

同时，确定当前时段t每个抽样场景下已投入运行的可调节型火电机组的负荷调节充裕度区间。该区间实际为

其中，

为当前时段t可调节型水电机组向下调节极限值，

为当前时段t每个抽样场景下可调节型火电机组运行下界且

为当前时段t第Gk台可调节型火电机组的运行下界且

P_Gk，min为第Gk台可调节型火电机组出力下限，RGk，d为第Gk台可调节型火电机组向下爬坡速率。

为当前时段t可调节型水电机组向上调节极限值，

为当前时段t每个抽样场景下可调节型火电机组运行上界且

为当前时段t第Gk台可调节型火电机组运行上界且P_Gk，max为第Gk台可调节型火电机组的出力上限，R_Gk，u为第Gk台可调节型火电机组向上爬坡速率。为当前时段t第Gk台可调节型火电机组投入运行状态，当前时段t第Gk台可调节型火电机组投入运行，则

当前时段t第Gk台可调节型火电机组未投入运行，则

Δt为当前时段t的时间长度，N_T为可调节型火电机组的数量，m=1，2，...，N，N为每个时段中的抽样场景数目。At the same time, determine the load regulation adequacy interval of the adjustable thermal power units that have been put into operation under each sampling scenario in the current period t. The interval is actually

in,

Adjust the limit value downward for the adjustable hydroelectric unit in the current period t,

is the lower bound of the operation of adjustable thermal power units in each sampling scenario in the current period t and

is the lower bound of the operation of the Gkth adjustable thermal power unit in the current period and

P _Gk,min is the output lower limit of the Gkth adjustable thermal power unit, RGk,d is the downward climbing rate of the Gkth adjustable thermal power unit.

is the upward adjustment limit value of the adjustable hydroelectric unit in the current period t,

is the upper bound of the operation of the adjustable thermal power unit in each sampling scenario in the current period t and

is the upper bound of the operation of the Gkth adjustable thermal power unit in the current period and P _Gk,max is the output upper limit of the Gkth adjustable thermal power unit, R _Gk,u is the upward climbing rate of the Gkth adjustable thermal power unit. The Gkth adjustable thermal power unit is put into operation in the current period t, and the Gkth adjustable thermal power unit is put into operation in the current period t, then

Δt is the time length of the current period t, _NT is the number of adjustable thermal power units, m=1, 2,..., N, N is the number of sampling scenes in each period.

当前时段t的每个抽样场景下，都会计算得到一个值，选择所有

值（m=1，2，...，N）中最小的值，当

取最小值时的抽样场景值设为τ，则抽样场景τ下已投入运行的可调节型火电机组及其分配的负荷

就是当前时段t可调节型火电机组的最佳投入运行方案。同时，抽样场景τ下已投入运行的可调节型火电机组的负荷调节充裕度区间

即为当前时段t可调节型火电机组的最佳负荷调节充裕度区间。In each sampling scenario of the current period t, a value, select all

The smallest value among the values (m=1, 2, ..., N), when

The sampling scene value when taking the minimum value is set to τ, then the adjustable thermal power units that have been put into operation under the sampling scene τ and their distributed load

It is the best operation plan for the t-adjustable thermal power unit in the current period. At the same time, the load regulation adequacy range of the adjustable thermal power units that have been put into operation under the sampling scenario τ

That is, the optimal load regulation adequacy interval of the adjustable thermal power unit in the current period t.

步骤10：令t=t+1，计算当前时段和当前时段之前连续4个时段的第二等值负荷最大值，以及当前时段和当前时段之前连续4个时段的第二等值负荷最小值。Step 10: Let t=t+1, calculate the second equivalent load maximum value of the current period and the 4 consecutive periods before the current period, and the second equivalent load minimum value of the current period and the 4 consecutive periods before the current period.

如果当前时段的前一个时段的可调节型火电机组投入运行的总装机容量小于当前时段和当前时段之前连续4个时段的第二等值负荷最大值，则按照第一设定顺序投入还未投入运行的可调节型火电机组，直至投入运行的可调节型火电机组的总装机容量大于当前时段和当前时段之前连续4个时段的第二等值负荷最大值。If the total installed capacity of the adjustable thermal power units put into operation in the previous period of the current period is less than the second equivalent load maximum value of the current period and the 4 consecutive periods before the current period, it will be put into use according to the first setting order. The operating adjustable thermal power units until the total installed capacity of the adjustable thermal power units put into operation is greater than the second equivalent load maximum value of the current period and the 4 consecutive periods before the current period.

如果当前时段和当前时段之前连续4个时段的第二等值负荷最大值小于当前时段的前一个时段的可调节型火电机组的最小技术出力总和，则按照第一设定顺序的倒序停止已经投入运行的可调节型火电机组，直至投入运行的可调节型火电机组的最小技术出力总和小于当前时段和当前时段之前连续4个时段的第二等值负荷最小值。If the maximum value of the second equivalent load in the current period and the 4 consecutive periods before the current period is less than the sum of the minimum technical output of the adjustable thermal power unit in the previous period of the current period, stop the input in the reverse order of the first setting order The adjustable thermal power units in operation until the sum of the minimum technical output of the adjustable thermal power units put into operation is less than the second equivalent load minimum value of the current period and the 4 consecutive periods before the current period.

按照第一设定顺序的倒序停止已经投入运行的可调节型火电机组实际就是，先按输出功率调节速率由低到高的顺序停止燃气轮机机组，再按最小比耗量由高到低的顺序停止可停机的可调节型火电机组，最后按最小运行比耗量由高到低的顺序停止不可停机的可调节型火电机组。Stopping the adjustable thermal power units that have been put into operation in the reverse order of the first setting order is actually to stop the gas turbine units in the order of output power adjustment rate from low to high, and then stop in the order of minimum specific consumption from high to low The adjustable thermal power units that can be stopped, and finally stop the adjustable thermal power units that cannot be stopped according to the order of the minimum operating specific consumption from high to low.

步骤11：判断当前时刻投入运行的可调节型火电机组与当前时刻的前一时刻投入运行的可调节型火电机组是否相同。Step 11: Determine whether the adjustable thermal power unit put into operation at the current moment is the same as the adjustable thermal power unit put into operation at the moment before the current moment.

如果当前时刻投入运行的可调节型火电机组与当前时刻的前一时刻投入运行的可调节型火电机组相同，则说明没有新增加的发电机组，也没有新关闭的发电机组，此时需要对以投入运行的可调节型发电机组进行微调，直接跳至步骤12。If the adjustable thermal power unit put into operation at the current moment is the same as the adjustable thermal power unit put into operation at the previous moment at the current moment, it means that there is no newly added generator set, and there is no newly closed generator set. At this time, the following Fine-tuning the adjustable generator set put into operation, skip to step 12 directly.

如果当前时刻投入运行的可调节型火电机组与当前时刻的前一时刻投入运行的可调节型火电机组不相同，则说明有新增加的机组或者新关闭的机组，直接步骤15，选择当前时刻最佳投入运行方案和最佳负荷调节充裕度区间。If the adjustable thermal power unit put into operation at the current moment is different from the adjustable thermal power unit put into operation at the previous moment at the current moment, it means that there is a newly added unit or a newly closed unit, go to step 15 and select the most The best put into operation scheme and the best load adjustment margin range.

步骤12：计算当前时段t每个抽样场景下可调节型火电机组负荷变化量Step 12: Calculate the load variation of the adjustable thermal power unit under each sampling scenario in the current period t

$ΔP ΔP = = {P P}_{L L 33,, t t}^{m m} - - {P P}_{L L 33,, t t - - 11}^{τ τ} . .$

如果ΔP>0，则说明负荷增加了，需要对投入运行的可调节型发电机组进行微调，此时跳至步骤13，使得已经投入运行的可调节型发电机组承担增加的负荷变化量。If ΔP>0, it means that the load has increased, and it is necessary to fine-tune the adjustable generator set that has been put into operation. At this time, skip to step 13, so that the adjustable generator set that has been put into operation can bear the increased load change.

如果ΔP<0，则说明负荷减少了，需要对投入运行的可调节型发电机组进行微调，此时跳至步骤14，使得已经投入运行的可调节型发电机组承担减少的负荷变化量。If ΔP<0, it means that the load has decreased, and it is necessary to fine-tune the adjustable generator set that has been put into operation. At this time, skip to step 14, so that the adjustable generator set that has been put into operation bears the reduced load change.

如果ΔP=0，则无需对投入运行的可调节型发电机组进行微调，直接跳至步骤15，选择当前时刻最佳投入运行方案和最佳负荷调节充裕度区间。If ΔP=0, there is no need to fine-tune the adjustable generator set that is put into operation, and skip directly to step 15 to select the best put-in-run scheme and the best load adjustment margin range at the current moment.

步骤13：按照第二设定顺序调整当前时段已投入运行且尚有可调节容量的可调节型发电机组，直至参与调整的可调节型发电机组满足 $ΔP = Σ_{Hi = 1}^{N_{H}} λ_{Hi} P_{Hi, t, u} + Σ_{Gk = 1}^{N_{G}} λ_{Gk} P_{Gk, t, u}^{m},$ 而后执行步骤15。Step 13: Adjust the adjustable generator sets that have been put into operation in the current period and still have adjustable capacity according to the second setting sequence until the adjustable generator sets participating in the adjustment meet the $ΔP = Σ_{hi = 1}^{N_{h}} λ_{hi} P_{hi, t, u} + Σ_{K = 1}^{N_{G}} λ_{K} P_{K, t, u}^{m},$ Then execute step 15.

其中，λ_Hi为当前时段第Hi台可调节型水电机组参与调整的状态，当前时段第Hi台可调节型水电机组参与调整时，λ_Hi=1；当前时段第Hi台可调节型水电机组未参与调整时，λ_Hi=0。P_Hi，t，u为当前时段第Hi台可调节型水电机组的出力增加变化量。λ_Gk为当前时段投入运行的火电机组中第Gk台火电机组参与调整的状态，当前时段投入运行的火电机组中第Gk台火电机组参与调整时，λ_GK=1；当前时段投入运行的火电机组中第Gk台火电机组未参与调整时，λ_Gk=0。

为当前时段第m个抽样背景下投入运行的火电机组中第Gk台火电机组出力增加变化量，N_H为可调节型水电机组的数量，N_G为已投入运行的可调节型火电机组的数量。Among them, _λHi is the state that the Hi adjustable hydroelectric unit participates in the adjustment in the current period. When the Hi adjustable hydroelectric unit participates in the adjustment in the current period, _λHi =1; the Hi adjustable hydroelectric unit in the current period does not When participating in the adjustment, λ _Hi =0. P _{Hi, t, u} is the output increase variation of the Hi-th adjustable hydroelectric unit in the current period. λ _Gk is the status of the Gkth thermal power unit participating in the adjustment among the thermal power units put into operation _in the current period. When the Gkth thermal power unit does not participate in the adjustment, λ _Gk =0.

is the output increase variation of the Gkth thermal power unit put into operation under the background of the mth sampling in the current period, N _H is the number of adjustable hydropower units, and N _G is the number of adjustable thermal power units that have been put into operation .

按照第二设定顺序调整当前时段已投入运行且尚有可调节容量的可调节型发电机组具体为，先调整当前时段已投入运行且尚有可调节容量的可调节型水电机组，当前时段已投入运行的所有可调节型水电机组没有可调节容量时，按输出功率调节速率由高到低的顺序调整已投入运行的燃气轮机机组，当前时段已投入运行的所有燃气轮机机组没有可调节容量时，按最小比耗量由低到高的顺序调整投入运行的可停机的可调节型火电机组，当前时段已投入运行的所有可停机的可调节型火电机组没有可调节容量时，按最小运行比耗量由低到高的顺序调整不可停机的可调节型火电机组。According to the second setting sequence, adjust the adjustable hydropower units that have been put into operation and still have adjustable capacity in the current period. When all the adjustable hydropower units put into operation have no adjustable capacity, adjust the gas turbine units that have been put into operation according to the order of output power adjustment rate from high to low. When all the gas turbine units that have been put into operation in the current period have no adjustable capacity, press The order of the minimum specific consumption from low to high is to adjust the adjustable thermal power units that are put into operation and can be stopped. When all the adjustable thermal power units that have been put into operation in the current period have no adjustable capacity, the minimum operating specific consumption Adjust non-stop adjustable thermal power units in sequence from low to high.

$| | ΔP ΔP | | = = {Σ Σ}_{Hi hi = = 11}^{{N N}_{H h}} {λ λ}_{Hi hi} {P P}_{Hi hi,, t t,, d d} + + {Σ Σ}_{Gk K = = 11}^{{N N}_{G G}} {λ λ}_{Gk K} {P P}_{Gk K,, t t,, d d}^{m m} . .$

其中，λ_Hi为当前时段第Hi台可调节型水电机组参与调整的状态，当前时段第Hi台可调节型水电机组参与调整时，λ_Hi=1；当前时段第Hi台可调节型水电机组未参与调整时，λ_Hi=0。P_Hi，t，d为当前时段第Hi台可调节型水电机组的出力减少变化量。λ_Gk为当前时段投入运行的火电机组中第Gk台火电机组参与调整的状态，当前时段投入运行的火电机组中第Gk台火电机组参与调整时，λ_GK=1；当前时段投入运行的火电机组中第Gk台火电机组未参与调整时，λ_Gk=0。

为当前时段第m个抽样背景下投入运行的火电机组中第Gk台火电机组出力减少变化量，N_H为可调节型水电机组的数量，N_G为已投入运行的可调节型火电机组的数量。Among them, _λHi is the state that the Hi adjustable hydroelectric unit participates in the adjustment in the current period. When the Hi adjustable hydroelectric unit participates in the adjustment in the current period, _λHi =1; the Hi adjustable hydroelectric unit in the current period does not When participating in the adjustment, λ _Hi =0. P _{Hi, t, d} is the output reduction variation of the Hi-th adjustable hydroelectric unit in the current period. λ _Gk is the status of the Gkth thermal power unit participating in the adjustment among the thermal power units put into operation _in the current period. When the Gkth thermal power unit does not participate in the adjustment, λ _Gk =0.

is the output reduction variation of the Gkth thermal power unit put into operation under the background of the mth sampling in the current period, N _H is the number of adjustable hydropower units, and N _G is the number of adjustable thermal power units that have been put into operation .

同时计算当前时段每个抽样场景下已投入运行的可调节型火电机组的负荷调节充裕度区间。Step 15: Assign loads to the adjustable thermal power units that have been put into operation under each sampling scenario in the current period

At the same time, the load regulation adequacy interval of the adjustable thermal power units that have been put into operation under each sampling scenario in the current period is calculated.

该步骤的具体过程与步骤8相同，此处不再赘述。The specific process of this step is the same as step 8, and will not be repeated here.

作为当前时段t可调节型火电机组的最佳投入运行方案。同时，选择运行成本的最小值对应的抽样场景τ下已投入运行的可调节型火电机组的负荷调节充裕度区间作为当前时段t可调节型火电机组的最佳负荷调节充裕度区间。Step 16: Calculate the operating cost of the adjustable thermal power units that have been put into operation under each sampling scenario of the current period t, and select the adjustable thermal power units that have been put into operation under the sampling scenario τ corresponding to the minimum value of the operating cost and their allocation load

As the best operation plan for the current period t adjustable thermal power unit. At the same time, the load regulation adequacy interval of the adjustable thermal power unit that has been put into operation under the sampling scenario τ corresponding to the minimum value of the operating cost is selected as the optimal load regulation adequacy interval of the adjustable thermal power unit in the current period t.

该步骤的具体过程与步骤9相同，此处不再赘述。The specific process of this step is the same as that of step 9, and will not be repeated here.

步骤17：判断t>N是否成立，如果t>N，则说明整个调度周期结束，此时执行步骤18；否则，返回步骤10，选择下一个时刻最佳投入运行方案和最佳负荷调节充裕度区间。Step 17: Judging whether t>N holds true, if t>N, it means that the entire dispatching cycle is over, and then execute step 18; otherwise, return to step 10, and select the best put into operation plan and the best load adjustment adequacy at the next time interval.

步骤18：结束。Step 18: End.

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。The above is only a preferred embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any person skilled in the art can easily conceive of changes or modifications within the technical scope disclosed in the present invention. Replacement should be covered within the protection scope of the present invention. Therefore, the protection scope of the present invention should be determined by the protection scope of the claims.

Claims

1. a multi-energy resource combined electric generating system is moved the optimization method of simulating, and it is characterized in that described method comprises:

Step 1: the configuration scheduling cycle also will be divided into T period described dispatching cycle;

Step 2: set each period sampling number of scenes N, and determine that each period wind-powered electricity generation unit and photoelectricity unit of each sampling scene blends the power value

wherein, m=1,2 ..., N, j=1,2 ..., T;

Step 3: according to formula P _{l1, j}=P _{l, j}-P _{f, j}calculate the load P that in each period system, the adjustable type generating set be should bear _{l1, j}, and by P _{l1, j}be designated as the first duty value such as grade;

In described system, the adjustable type generating set comprises adjustable type Hydropower Unit and adjustable type fired power generating unit;

P _{l, j}be j the total load that the period system be should bear;

P _{f, j}for the load that in j period system, the generating set except the adjustable type generating set provides;

Step 4: according to formula

calculate the load P that in each period system, the adjustable type fired power generating unit be should bear _{l2, j}, and by P _{l2, j}be designated as the second duty value such as grade; Simultaneously, calculate each period adjustable type Hydropower Unit accommodation limit value that makes progress

with downward accommodation limit value

μ _{hi, j}be the operating state of j period Hi adjustable type Hydropower Unit, when j period, Hi adjustable type Hydropower Unit put into operation, μ _{hi, j}=1; When j period, Hi adjustable type Hydropower Unit do not put into operation, μ _{hi, j}=0;

P _{hi, j}be power output and the P of j period Hi adjustable type Hydropower Unit _{hi, j}=η _hi* g * (H _{up, Hi, j}-H _{down, Hi, j}) * Q _{hi, j};

η _hiit is the generating efficiency coefficient of Hi adjustable type Hydropower Unit;

G is gravity constant;

H _{up, Hi, j}it is water level on the dam of j period Hi adjustable type Hydropower Unit;

H _{down, Hi, j}it is the level of tail water of j period Hi adjustable type Hydropower Unit;

Q _{hi, j}it is the generating average discharge of j period Hi adjustable type Hydropower Unit;

N _hquantity for the adjustable type Hydropower Unit;

μ _{hp, j}be the state that draws water of j period pumped storage machine, when j the period, pumped storage machine put into operation, μ _{hp, j}=1, when j the period, pumped storage machine did not put into operation, μ _{hp, j}=0;

P _{hp, j}be the power that draws water of j period pumped storage machine;

Step 5: the state that initially puts into operation of determining the adjustable type fired power generating unit; Specifically, calculate respectively the duty value minimum value such as second of the second duty value maximums such as grade of front 5 periods of dispatching cycle and front 5 periods, sequentially drop into the adjustable type fired power generating unit until the total installation of generating capacity of the adjustable type fired power generating unit put into operation is greater than the second maximum such as duty value such as grade of front 5 periods of dispatching cycle according to the first setting, and the minimum technology of the adjustable type fired power generating unit put into operation the summation of exerting oneself is less than the duty value minimum values such as second of front 5 periods of dispatching cycle;

Described adjustable type fired power generating unit comprises the adjustable type fired power generating unit that can not shut down, adjustable type fired power generating unit and the gas turbine unit that can shut down;

Step 6: set initial period t=5;

Step 7: according to formula

calculate the load that current period t adjustable type fired power generating unit be should bear under each sampling scene

and will

be designated as the C grade duty value;

Step 8: the adjustable type fired power generating unit distribution load put into operation under each the sampling scene for current period t

determine the Load Regulation abundant intensity interval of the adjustable type fired power generating unit put into operation under each sampling scene of current period t simultaneously;

Step 9: calculate the operating cost of the adjustable type fired power generating unit that each sampling of current period t put into operation under scene, select the adjustable type fired power generating unit that put into operation under sampling scene τ corresponding to the minimum value of operating cost and the load of distribution thereof

as the best of current period t adjustable type fired power generating unit scheme that puts into operation; Simultaneously, select the Load Regulation abundant intensity interval of the adjustable type fired power generating unit that put into operation under sampling scene τ corresponding to the minimum value of operating cost to regulate the abundant intensity interval as the optimum load of current period t adjustable type fired power generating unit;

Step 10: make t=t+1, calculate the duty value maximum such as second of continuous 4 periods before current period and current period, and the duty value minimum value such as second of continuous 4 periods before current period and current period;

If the total installation of generating capacity that the adjustable type fired power generating unit of the previous period of current period puts into operation is less than the duty value maximum such as second of continuous 4 periods before current period and current period, according to the first setting, sequentially drop into the adjustable type fired power generating unit also do not put into operation, until the total installation of generating capacity of the adjustable type fired power generating unit put into operation is greater than the duty value maximum such as second of continuous 4 periods before current period and current period;

The summation if the minimum technology of adjustable type fired power generating unit that the duty value maximum such as second of continuous 4 periods is less than the previous period of current period before current period and current period is exerted oneself, the adjustable type fired power generating unit that stops having put into operation according to the inverted order of the first setting order, until the minimum technology of the adjustable type fired power generating unit put into operation is exerted oneself, summation is less than the duty value minimum value such as second of continuous 4 periods before current period and current period;

Step 11: judge that whether the adjustable type fired power generating unit that adjustable type fired power generating unit that current time puts into operation puts into operation with the previous moment of current time is identical, if the adjustable type fired power generating unit that the adjustable type fired power generating unit that current time puts into operation puts into operation with the previous moment of current time is identical, perform step 12; Otherwise, perform step 15;

Step 12: calculate adjustable type fired power generating unit load variations amount under each sampling scene of current period t

if Δ P>0, perform step 13; If Δ P<0, perform step 14; If Δ P=0, perform step 15;

Step 13: sequentially adjust according to the second setting the adjustable type generating set that the current period put into operation and adjustable capacity is still arranged, until the adjustable type generating set that participates in adjusting meets

ΔP = Σ_{Hi = 1}^{N_{H}} λ_{Hi} P_{Hi, t, u} + Σ_{Gk = 1}^{N_{G}} λ_{Gk} P_{Gk, t, u}^{m};

Then perform step 15;

Wherein, λ _hifor current period Hi platform adjustable type Hydropower Unit participates in the state of adjusting, when current period Hi platform adjustable type Hydropower Unit participates in adjusting, λ _hi=1; Current period Hi platform adjustable type Hydropower Unit has neither part nor lot in while adjusting, λ _hi=0;

P _{hi, t, u}the increase variable quantity of exerting oneself for current period Hi platform adjustable type Hydropower Unit;

λ _gkin the fired power generating unit put into operation for the current period, Gk platform fired power generating unit participates in the state of adjusting, when in the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit participates in adjusting, and λ _gK=1; In the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit has neither part nor lot in while adjusting, λ _gk=0;

for the increase variable quantity of exerting oneself of Gk platform fired power generating unit in the fired power generating unit put into operation under m of current period sampling background;

N _hquantity for the adjustable type Hydropower Unit;

N _gquantity for the adjustable type fired power generating unit that put into operation;

Step 14: sequentially adjust according to the second setting the adjustable type generating set that the current period put into operation and adjustable capacity is still arranged, until the adjustable type generating set that participates in adjusting meets

| ΔP | = Σ_{Hi = 1}^{N_{H}} λ_{Hi} P_{Hi, t, d} + Σ_{Gk = 1}^{N_{G}} λ_{Gk} P_{Gk, t, d}^{m};

P _{hi, t, d}the minimizing variable quantity of exerting oneself for current period Hi platform adjustable type Hydropower Unit;

for the minimizing variable quantity of exerting oneself of Gk platform fired power generating unit in the fired power generating unit put into operation under m of current period sampling background;

N _hquantity for the adjustable type Hydropower Unit;

Step 15: the adjustable type fired power generating unit distribution load put into operation under each the sampling scene for the current period

calculate the Load Regulation abundant intensity interval of the adjustable type fired power generating unit put into operation under each sampling scene of current period simultaneously;

Step 16: calculate the operating cost of the adjustable type fired power generating unit that each sampling of current period t put into operation under scene, select the adjustable type fired power generating unit that put into operation under sampling scene τ corresponding to the minimum value of operating cost and the load of distribution thereof

Step 17: judge t > whether N sets up, if t > N, perform step 18; Otherwise, return to step 10;

Step 18: finish.

2. method according to claim 1, is characterized in that wind-powered electricity generation unit and the photoelectricity unit of described each sampling scene of definite each period blends the power value comprise following sub-step:

Sub-step 101: respectively according to the historical data of wind-powered electricity generation unit and photoelectricity unit output, make the histogram of wind-powered electricity generation unit output prediction deviation and the histogram of photoelectricity unit output prediction deviation;

Sub-step 102: by approximating method, obtain respectively each period wind-powered electricity generation unit output prediction deviation probability density function and photoelectricity unit output prediction deviation probability density function;

Sub-step 103: each period wind-powered electricity generation unit output prediction deviation probability density function and photoelectricity unit output prediction deviation probability density function are carried out respectively to integral operation, obtain wind-powered electricity generation unit output prediction deviation probability density function and the photoelectricity unit output prediction deviation probability density function of each period;

Sub-step 104: utilize formula

calculate the basic sampled value under each sampling scene of each period; Wherein, m=1,2 ..., N, N is each period sampling number of scenes, U is the random number on interval (0,1);

Sub-step 105: utilize formula

calculate each sampling scene leeward group of motors of each period prediction deviation of exerting oneself, and utilize formula

calculate photoelectricity unit output prediction deviation under each sampling scene of each period; Wherein,

for the inverse function of wind-powered electricity generation unit output prediction deviation probability density function,

inverse function for photoelectricity unit output prediction deviation probability density function; J=1,2 ..., T, the time hop count that T is dispatching cycle;

Sub-step 106: according to formula

the wind-powered electricity generation unit and the photoelectricity unit that calculate each sampling scene of each period blend the power value;

Wherein, P _{w, j}for each period wind-powered electricity generation unit output predicted value; P _{p, j}for each period photoelectricity unit output predicted value.

3. method according to claim 1, is characterized in that each period adjustable type Hydropower Unit of described calculating accommodation limit value that makes progress

adopt formula

Wherein, be the upwards pondage of j period Hi adjustable type Hydropower Unit, and

P_{Hi, j}^{op} = P_{Hi, j}^{\max} - P_{Hi, j};

it is the peak power output of j period Hi adjustable type Hydropower Unit;

P _{hi, j}it is the power output of j period Hi adjustable type Hydropower Unit;

μ _{hi, j}be the operating state of j period Hi adjustable type Hydropower Unit, when j period, Hi adjustable type Hydropower Unit put into operation, μ _{hi, j}=1; When j period, Hi adjustable type Hydropower Unit do not put into operation, μ _{hi, j}=0.

4. method according to claim 1, is characterized in that the downward accommodation limit value of each period adjustable type Hydropower Unit of described calculating

adopt formula

Wherein, be the downward pondage of j period Hi adjustable type Hydropower Unit, and

P_{Hi, j}^{ne} = P_{Hi, j} - P_{Hi, j}^{\min};

the minimum that is j period Hi adjustable type Hydropower Unit is forced power output;

5. method according to claim 1; it is characterized in that describedly according to the first setting, sequentially dropping into the adjustable type fired power generating unit specifically; first by minimum operating ratio consumption, order from low to high drops into the adjustable type fired power generating unit that can not shut down; by minimum specific consumption, order from low to high drops into the adjustable type fired power generating unit that can shut down again, and finally by the power output regulations speed, order from high to low drops into gas turbine unit.

6. method according to claim 1, is characterized in that the adjustable type fired power generating unit distribution load put into operation under described each sampling scene for current period t

specifically comprise following sub-step:

Sub-step 201: set the adjustable type fired power generating unit put into operation the out-of-limit sign nFlag=1 that exerts oneself;

Sub-step 202: according to formula calculate the micro-gaining rate of consumption of fired power generating unit; Wherein,

for C grade duty value and m=1,2 ..., N;

N is each period sampling number of scenes;

N ' _gfor current period t put into operation and the exerted oneself quantity of the adjustable type fired power generating unit that do not transfinite;

A _{gk '}for current period t the Gk ' platform put into operation and the exerted oneself quadratic term coefficient of consumption characteristic quadratic function of the adjustable type fired power generating unit that do not transfinite;

B _{gk '}for current period t the Gk ' platform put into operation and the exerted oneself Monomial coefficient of consumption characteristic quadratic function of the adjustable type fired power generating unit that do not transfinite;

Sub-step 203: make Gk=1;

Sub-step 204: according to formula

calculate adjustable type fired power generating unit value of exerting oneself of having put into operation of Gk platform under each sampling scene of current period t;

A _gkthe quadratic term coefficient of the consumption characteristic quadratic function of the adjustable type fired power generating unit put into operation for current period t Gk platform;

B _gkthe Monomial coefficient of the consumption characteristic quadratic function of the adjustable type fired power generating unit put into operation for current period t Gk platform;

Sub-step 205: if

transfinite, order

maximum output value during for this adjustable type fired power generating unit put into operation not out-of-limit, and make nFlag=0;

Sub-step 206: judgement Gk>N _gwhether set up, if Gk N _g, carry out sub-step 207; Otherwise, make Gk=Gk+1, return to sub-step 204;

N _gthe quantity of the adjustable type fired power generating unit put into operation for current period t;

Sub-step 207: judge whether nFlag=0 sets up, if nFlag=0 returns to sub-step 201; Otherwise, carry out sub-step 208;

Sub-step 208: make the load that under each sampling scene of current period t, the adjustable type fired power generating unit that put into operation of Gk platform provides be

7. method according to claim 1, is characterized in that the operating cost of the adjustable type fired power generating unit that put into operation under each sampling scene of the current period t of described calculating adopts formula

F_{\cos t} (P_{G, t}^{m}) = Σ_{Gk = 1}^{N_{T}} [μ_{Gk, t}^{m} \times f_{k} (P_{Gk, t}^{m}) + μ_{Gk, t}^{m} (1 - μ_{Gk, t - 1}^{m}) S_{Gk}];

Wherein, the operating cost of the adjustable type fired power generating unit put into operation under m the sampling scene for current period t, m=1,2 ..., N, N is the sampling number of scenes in each period;

N _tfor adjustable type fired power generating unit quantity;

for the current period t Gk platform adjustable type fired power generating unit state that puts into operation, current period t Gk platform adjustable type fired power generating unit puts into operation,

current period t Gk platform adjustable type fired power generating unit does not put into operation,

for the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit, and

f_{k} (P_{Gk, t}^{m}) = A_{Gk} \times {(P_{Gk, t}^{m})}^{2} + B_{Gk} \times P_{Gk, t}^{m} + C_{Gk};

A _gkquadratic term coefficient for the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit;

B _gkmonomial coefficient for the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit;

C _gkconstant term coefficient for the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit;

S _gkfor the payment for initiation of current period t Gk platform adjustable type fired power generating unit is used.

8. method according to claim 1, is characterized in that the Load Regulation abundant intensity interval of the adjustable type fired power generating unit that put into operation under each sampling scene of described current period t is

[P_{L 3, t}^{m} - P_{H, t}^{ne} - P_{G, t}^{ne, m}, P_{L 3, t}^{m} + P_{H, t}^{op} + P_{G, t}^{op, m}];

Wherein,

for the downward accommodation limit value of current period t adjustable type Hydropower Unit;

for adjustable type thermal power unit operation lower bound under each sampling scene of current period t, and

P_{G, t}^{ne, m} = Σ_{Gk = 1}^{N_{T}} μ_{Gk, t}^{m} P_{Gk, t}^{ne, m};

for the operation lower bound of current period t Gk platform adjustable type fired power generating unit, and

P_{Gk, t}^{ne, m} = \min (P_{Gk, t}^{m} - P_{Gk, \min}, - R_{Gk, d} \times Δt);

P _{gk, min}it is the Gk platform adjustable type fired power generating unit lower limit of exerting oneself;

R _{gk, d}it is the downward creep speed of Gk platform adjustable type fired power generating unit;

for the current period t Gk platform adjustable type fired power generating unit state that puts into operation, current period t Gk platform adjustable type fired power generating unit puts into operation, current period t Gk platform adjustable type fired power generating unit does not put into operation,

The time span that Δ t is current period t;

N _tquantity for the adjustable type fired power generating unit;

for the current period t adjustable type Hydropower Unit accommodation limit value that makes progress;

for the adjustable type thermal power unit operation upper bound under each sampling scene of current period t, and

P_{G, t}^{op, m} = Σ_{Gk = 1}^{N_{T}} μ_{Gk, t}^{m} P_{Gk, t}^{op, m};

for the current period t Gk platform adjustable type thermal power unit operation upper bound, and

P_{Gk, t}^{op, m} = \min (P_{Gk, \max} - P_{Gk, t}^{m}, - R_{Gk, u} \times Δt);

P _{gk, max}it is the upper limit of exerting oneself of Gk platform adjustable type fired power generating unit;

R _{gk, u}it is the Gk platform adjustable type fired power generating unit creep speed that makes progress;

M=1,2 ..., N, N is the sampling number of scenes in each period.

9. method according to claim 1, it is characterized in that describedly sequentially adjusting the current period according to the second setting and having put into operation and still had the adjustable type generating set of adjustable capacity to be specially, first adjust the adjustable type Hydropower Unit that the current period put into operation and adjustable capacity is still arranged, when all adjustable type Hydropower Unit of having put into operation of current period do not have adjustable capacity, the gas turbine unit that order adjustment has from high to low put into operation by the power output regulations speed, when all gas turbine units of having put into operation of current period do not have adjustable capacity, the adjustable type fired power generating unit of shutting down that order adjustment from low to high puts into operation by minimum specific consumption, when all adjustable type fired power generating unit of shutting down of having put into operation of current period do not have adjustable capacity, the adjustable type fired power generating unit that order adjustment from low to high can not be shut down by minimum operating ratio consumption.