CN114884134A - Thermal power generating unit flexibility adjusting and scheduling method based on interval optimization - Google Patents

Thermal power generating unit flexibility adjusting and scheduling method based on interval optimization Download PDF

Info

Publication number
CN114884134A
CN114884134A CN202210574276.5A CN202210574276A CN114884134A CN 114884134 A CN114884134 A CN 114884134A CN 202210574276 A CN202210574276 A CN 202210574276A CN 114884134 A CN114884134 A CN 114884134A
Authority
CN
China
Prior art keywords
generating unit
power generating
thermal power
cost
flexibility
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210574276.5A
Other languages
Chinese (zh)
Inventor
李明扬
张明瑞
陈丰
林忠伟
王林
刘晓航
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Xian Thermal Power Research Institute Co Ltd
Huaneng Group Technology Innovation Center Co Ltd
Huaneng Nanjing Jinling Power Generation Co Ltd
Original Assignee
North China Electric Power University
Xian Thermal Power Research Institute Co Ltd
Huaneng Group Technology Innovation Center Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by North China Electric Power University, Xian Thermal Power Research Institute Co Ltd, Huaneng Group Technology Innovation Center Co Ltd filed Critical North China Electric Power University
Priority to CN202210574276.5A priority Critical patent/CN114884134A/en
Publication of CN114884134A publication Critical patent/CN114884134A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381Dispersed generators
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • H02J2300/28The renewable source being wind energy
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/40Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation wherein a plurality of decentralised, dispersed or local energy generation technologies are operated simultaneously
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Power Engineering (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Wind Motors (AREA)

Abstract

The invention discloses a thermal power generating unit flexibility adjusting and scheduling method based on interval optimization, and belongs to the technical field of power system optimization. Step 1: establishing a flexibility regulation certainty model of the thermal power generating unit; step 2: establishing a flexibility adjusting interval optimization model of the thermal power generating unit; and step 3: solving the thermal power generating unit flexibility adjusting interval optimization model established in the step 2 by using a GSOMP algorithm; using a GSOMP algorithm to cooperatively optimize the mean value and the deviation value of the target interval and generating a candidate scheme with the characteristics of the comprehensive cost intervals in different years; and 4, step 4: and (3) evaluating the annual comprehensive cost intervals in the candidate schemes in the step (3) by using a multi-attribute decision method, and selecting a group of optimal schemes according to the comprehensive performance of the characteristic quantities of each interval. The method can better adapt to the uncertainty of the power flow under the high wind power access proportion, effectively relieves the blocking phenomenon of the power flow, and reduces the abandoned wind volume under each uncertain scene.

Description

Thermal power generating unit flexibility adjusting and scheduling method based on interval optimization
Technical Field
The invention relates to the technical field of power system optimization, in particular to a thermal power generating unit flexibility adjusting and scheduling method based on interval optimization.
Background
Under the promotion of rapid development of economy, the proportion of renewable energy sources is continuously increased in consideration of environmental protection and persistent development, and the utilization and development of the renewable energy sources become important problems in the world, so that the energy consumption of users is more economical and low-carbon. Renewable energy is accessed into a power grid in a large scale, higher requirements are provided for the peak regulation capacity of the power grid, the existing power supply structure in China is used, the installed capacity and the generated energy are mainly thermal power, and the main mode for realizing the renewable energy consumption of the power grid is also the regulation of a thermal power generating unit at present. The thermal power generating unit has more flexibility and better peak regulation and frequency modulation capability and has years of peak regulation operation experience of the coal-fired unit.
In order to solve the problem of renewable energy consumption, improve the flexibility of a power system, assist in achieving the carbon neutralization goal, and China actively explores a thermal power generating unit flexibility adjusting and optimizing scheduling method. The flexibility of the power system not only comprises definition and analysis of the flexibility, but also comprises various contents such as related evaluation indexes of the flexibility of the power system, related products of the flexibility and the like, and the existing related research mainly takes large-scale renewable intermittent power supply grid connection as a background, discusses how to better maintain the active balance of the power system in the aspects of concepts, evaluation systems and the like, and whether the flexibility regulation and optimization scheduling is economical and reasonable, and how to have difference of comprehensive benefits among different methods, and further research and discussion are needed. The economic evaluation system for the flexibility regulation and optimization scheduling technology of the thermal power generating unit is researched, a flexibility regulation comprehensive evaluation model of the thermal power generating unit is constructed, the reasonable adjustment of the power supply structure layout in China is facilitated, the consumption and storage capacity of renewable energy sources are improved, and the efficient utilization of the renewable energy sources is promoted. The research on the flexibility of the power system not only has academic value, but also has important significance on the grid-connected access and market consumption of renewable energy sources in China, and a scientific and reasonable flexibility evaluation system of the power system is constructed, so that the transformation of the energy sources in China can be better assisted.
When the thermal power generating unit participates in deep peak shaving, flexibility adjustment and optimization of the thermal power generating unit are required to be performed in advance so as to reduce the lowest output limit of the thermal power generating unit. In the past, only a flexible adjustment optimization model suitable for a thermal power generating unit in a large-scale wind power access scene is considered, the system economy of the thermal power generating unit under different flexible adjustment schemes is analyzed, and only main factors influencing the peak regulation capacity of the thermal power generating unit and influences of the thermal power generating unit on the system under different peak regulation depths are analyzed. Because the unit can constantly carry out large-scale regulation and reduction power and frequent start and stop, the running cost is increased, the service life is reduced, and the problems faced by the current peak regulation are difficult to solve only by depending on the current peak regulation capability of the system. Therefore, the influence of a network topological structure on flexibility adjustment and optimization scheduling needs to be deeply excavated, and peak shaving of the thermal power generating unit is assisted by combining the adaptability of the scheme under multiple scenes, so that the wind and light abandoning amount is reduced.
Aiming at the problem that the unit power generation cost is greatly fluctuated due to the fact that the power generation cost of a thermal power generating unit is greatly increased in the deep peak regulation stage, the operation cost optimization of multiple scenes needs to be achieved through a full-scene optimization method, and a scheme with better comprehensive economy is provided. In the aspect of the trend adaptability, considering that a power grid needing flexibility modification usually contains high-permeability wind power, a full-scene optimization method is adopted to process the trend uncertainty under the high-permeability wind power, so that a scheme with better trend adaptability is provided.
Disclosure of Invention
The invention aims to provide a thermal power generating unit flexibility adjusting and scheduling method based on interval optimization, which is characterized by comprising the following steps of:
step 1: comprehensively considering flexibility adjusting cost and operating cost of the thermal power generating unit, and establishing a flexibility adjusting certainty model of the thermal power generating unit by taking annual comprehensive cost as an optimization target;
step 2: the method comprises the steps of limiting the boundary of a scene set by limiting the fluctuation range of uncertain scenes of each node, and establishing a thermal power generating unit flexibility adjusting interval optimization model;
and step 3: solving the thermal power generating unit flexibility adjusting interval optimization model established in the step 2 by using a GSOMP algorithm; using a GSOMP algorithm to cooperatively optimize the mean value and the deviation value of the target interval and generating a candidate scheme with the characteristics of the comprehensive cost intervals in different years;
and 4, step 4: and (3) evaluating the annual comprehensive cost intervals in the candidate schemes in the step (3) by using a multi-attribute decision method, and selecting a group of optimal schemes according to the comprehensive performance of the characteristic quantities of each interval.
The thermal power unit flexibility adjustment certainty model in the step 1 is specifically as follows:
an objective function: the minimum annual comprehensive cost is taken as an optimization target
min(C B +C O ) In the formula (1), C B And C O Respectively adjusting the annual flexibility adjusting cost and the annual operation cost of the thermal power generating unit;
the annual transformation cost of the thermal power generating unit is as follows (2):
Figure BDA0003661442030000021
in the formula, P i G,min And P i G,min0 Respectively representing the lowest output of the thermal power generating unit i before and after modification; c y_rebuild Annual flexibility modification costs per unit capacity; m G For system medium fireThe number of the motor sets;
the annual flexibility transformation cost per unit capacity is converted into the following formula (3):
Figure BDA0003661442030000031
in the formula, C rebuild Cost for flexibility modification per unit capacity; y is the unit life; eta is annual interest rate;
the middle-aged operating cost of the system is as follows:
Figure BDA0003661442030000032
in the formula, D m Days of month m; t is the number of hours of operation in each typical day;
Figure BDA0003661442030000033
the output of the thermal power generating unit i in a typical day period t in a month m is obtained; lambda [ alpha ] w Punishment cost of unit air volume abandon; m W The number of wind generating sets in the system;
Figure BDA0003661442030000034
the wind power is the wind curtailment quantity of the wind turbine generator j in the typical day period t in the month m;
and power balance constraint:
Figure BDA0003661442030000035
in the formula (I), the compound is shown in the specification,
Figure BDA0003661442030000036
the output of the wind turbine generator j is the output of the wind turbine generator j in a typical day period t in a month m, and the wind turbine generator output refers to the predicted output of the wind turbine generator;
Figure BDA0003661442030000037
for load, node k is active power at typical day time t in month m; m D Is the number of load nodes in the system;
output restraint of the thermal power generating unit:
Figure BDA0003661442030000038
in the formula, P i G,max The maximum output of the thermal power generating unit i is obtained;
flexibility modification constraint of the thermal power generating unit:
Figure BDA0003661442030000039
Figure BDA00036614420300000311
in the formula, P i G,min1 And P i G,min2 Respectively carrying out non-oil-feeding peak regulation transformation and oil-feeding peak regulation transformation on the thermal power generating unit i to obtain the lowest output of the thermal power generating unit; x is the number of i,0 、x i,1 、x i,2 The variable of the transformation scheme decision of the unit i is 0-1 variable, x i,0 =1、x i,1 =1、x i,2 When the unit i is 1, the unit i is not transformed, is not transformed in an oil feeding mode, and is transformed in an oil feeding mode;
the climbing of thermal power generating unit restricts:
Figure BDA00036614420300000310
in the formula (I), the compound is shown in the specification,
Figure BDA0003661442030000041
and
Figure BDA0003661442030000042
the maximum upward climbing speed and the maximum downward climbing speed of the thermal power generating unit i are respectively set;
and (3) line power flow constraint:
Figure BDA0003661442030000043
in the formula, G n-i 、G n-j 、G n-k Power transfer factors of a thermal power generating unit i, a wind power generating unit j and a load node k to a line n are respectively set; m L Is the number of lines in the system;
Figure BDA0003661442030000044
is the maximum allowed transmission power of line n;
abandoned air volume constraint of wind turbine generator
Figure BDA0003661442030000045
The thermal power generating unit flexibility adjusting interval optimization model in the step 2 is specifically as follows:
the system operating cost is as follows:
Figure BDA0003661442030000046
in the formula: c D The system operating cost;
Figure BDA0003661442030000047
and
Figure BDA0003661442030000048
respectively representing the lower limit and the upper limit of the output of the thermal power generating unit i in a typical day period t in a month m;
Figure BDA0003661442030000049
and
Figure BDA00036614420300000410
respectively setting a lower limit and an upper limit of the wind curtailment quantity of the wind turbine generator j in a typical day period t in a month m;
lower limit and upper limit of annual running cost interval:
Figure BDA00036614420300000411
Figure BDA00036614420300000412
in the formula:X D and
Figure BDA00036614420300000413
respectively the lower limit and the upper limit of the system operation cost;
power balance equation:
Figure BDA00036614420300000414
output restraint of the thermal power generating unit:
Figure BDA00036614420300000415
and (3) climbing restraint:
Figure BDA00036614420300000416
and (3) line power flow constraint:
Figure BDA00036614420300000417
wind abandon restraint:
Figure BDA0003661442030000051
the invention has the beneficial effects that:
the invention can effectively improve the comprehensive economy and the trend adaptability of the scheme; the model adopts a full scene optimization method, so that the running cost under the scene of large wind power output is obviously improved, the low running cost under the conventional wind power output scene is maintained, and the comprehensive economy of the scheme is improved; the full-scene optimization mode can better adapt to the tidal current uncertainty under the high wind power access proportion, effectively relieves the tidal current blocking phenomenon, and reduces the air abandonment rate under each uncertain scene, thereby obtaining a group of thermal power generating unit flexibility adjusting schemes with optimal comprehensive economy.
Drawings
Fig. 1 is a flowchart of a thermal power generating unit flexibility adjustment scheduling method based on interval optimization.
Detailed Description
The invention provides a thermal power generating unit flexibility adjusting and scheduling method based on interval optimization, and the method is further explained by combining an attached drawing and a specific embodiment.
Fig. 1 is a flowchart of a thermal power generating unit flexibility adjustment scheduling method based on interval optimization. The thermal power generating unit flexibility adjusting model based on interval optimization provided by the invention takes the annual comprehensive cost interval of thermal power generating unit flexibility adjusting cost and multi-scene operation cost as an optimization target, extracts the mean value and the deviation value of the target interval as the characteristic quantity of the interval, and selects the characteristic quantity as 2 optimization sub-targets in the model. In order to comprehensively measure the overall characteristics of a target interval and screen out a scheme which can meet the user requirements to the maximum extent, 2 steps of multi-target optimization and multi-attribute decision are used in model solution. Firstly, a GSOMP algorithm is used for cooperatively optimizing the mean value and the deviation value of a target interval, so that a plurality of candidate schemes are generated, and each scheme has different annual comprehensive cost interval characteristics. And secondly, evaluating annual comprehensive cost intervals in each scheme by using a multi-attribute decision method, and selecting an optimal scheme according to comprehensive performance of characteristic quantities in each interval. The method comprises the following specific steps:
1. establishing a flexibility regulation certainty model of a traditional thermal power generating unit
The established thermal power unit flexibility regulation deterministic model comprehensively considers flexibility regulation cost and operation cost of the thermal power unit. In order to facilitate the unification of long-term planning and calculation periods in operation, the optimization period of the model is selected to be 1 year, and the annual comprehensive cost is used as the optimization target of the model. The model is specifically expressed as follows:
min(C B +C O ) (1)
Figure BDA0003661442030000052
Figure BDA0003661442030000061
Figure BDA0003661442030000062
Figure BDA0003661442030000063
Figure BDA0003661442030000064
Figure BDA0003661442030000065
Figure BDA0003661442030000066
Figure BDA0003661442030000067
Figure BDA0003661442030000068
Figure BDA0003661442030000069
wherein: c B And C O Respectively adjusting the annual flexibility adjusting cost and the annual operation cost of the thermal power generating unit; p i G,min And P i G ,min0 Respectively representing the lowest output of the thermal power generating unit i before and after modification; c y_rebuild Cost for annual flexibility modification per unit volume; c rebuild Transforming the cost for the flexibility of unit capacity, and converting the cost into annual cost; y is the unit life; eta is annual interest rate; d m Days of month m; t is the number of hours of operation in each typical day;
Figure BDA00036614420300000610
and
Figure BDA00036614420300000611
respectively outputting power of a thermal power generating unit i and a wind power generating unit j in a typical day time t in a month m, wherein the output of the wind power generating unit refers to the predicted output of the wind power generating unit;
Figure BDA00036614420300000612
for load, node k is active power at a typical day period t in month m;
Figure BDA00036614420300000613
the wind power is the wind curtailment quantity of the wind turbine generator j in the typical day period t in the month m; p i G,max The maximum output of the thermal power generating unit i is obtained; p i G,min1 And P i G,min2 The minimum output of the thermal power generating unit i after non-oil-feeding peak regulation transformation and oil-feeding peak regulation transformation are respectively carried out; x is the number of i,0 、x i,1 、x i,2 The variable of the transformation scheme decision of the unit i is 0-1 variable, x i,0 =1、x i,1 =1、x i,2 When the unit i is 1, the unit i is not transformed, is not transformed in an oil feeding mode, and is transformed in an oil feeding mode; m is a group of G 、M W 、M D 、M L Respectively as a system medium-temperature motorThe number of groups, wind turbines, load nodes and lines; lambda [ alpha ] w Punishment cost of unit air volume abandon;
Figure BDA00036614420300000614
and
Figure BDA00036614420300000615
the maximum upward climbing speed and the maximum downward climbing speed of the thermal power generating unit i are respectively set;
Figure BDA00036614420300000616
is the maximum allowed transmission power of line n; g n-i 、G n-j 、G n-k And power transfer factors of the thermal power generating unit i, the wind power generating unit j and the load node k to the line n are respectively.
The optimization goal of the model represented by formula (1) is to minimize the annual combined cost, including the annual modification cost and the annual operation cost. And the formula (2) gives a formula of annual transformation cost of the thermal power generating unit. Equation (3) gives a conversion equation for the annual reform cost per unit volume. Equation (4) gives the equation for the annual operating costs in the system. Equation (5) is the power balance constraint in the system. And the formula (6) is the output constraint of the thermal power generating unit. The formula (7) and the formula (8) are flexibility modification constraints of the thermal power generating unit. Formula (9) is the climbing restraint of thermal power unit. Equation (10) is the line flow constraint. And the formula (11) is the wind abandoning amount constraint of the wind turbine generator.
The model can be represented in the general form:
Figure BDA0003661442030000071
wherein U is a scheme decision variable
Figure BDA0003661442030000072
The method is used for expressing the optimal flexibility transformation scheme of each thermal power generating unit; c is an optimization target; x is a state variable and is a state variable,
Figure BDA0003661442030000073
Figure BDA0003661442030000074
the method comprises the steps of obtaining various operation variables in optimal power flow calculation, such as the output of each thermal power generating unit and the air curtailment quantity of each wind power plant; g (X, U) and h (X, U) are equality constraint functions and inequality constraint functions respectively. The optimal lower output limit P of each thermal power generating unit can be obtained according to the decision variable U i G,min0
2. Establishing thermal power generating unit flexibility adjusting interval optimization model
Compared with the traditional thermal power generating unit, the output of the wind power generating unit has stronger fluctuation, which causes higher uncertainty of the output of the wind power in actual calculation. When a plurality of wind power plants fluctuate within the output range interval of each time period, a plurality of uncertain scenes are generated, wherein the uncertain scene s can be expressed as:
Figure BDA0003661442030000075
in the formula:
Figure BDA0003661442030000076
the output of a wind generating set j in a month M under an uncertain scene s in a time period t, wherein j is 1,2, …, M W ,t=1,2,…,T;P s,m And (4) the wind turbine generator set output vector in the month m under the uncertain scene s.
From a given scene s expression P s,m It can be seen that a large number of uncertain scenes can be derived after considering the volatility of a plurality of wind power plants. Under the condition that a large number of uncertain scenes exist in the system, the interval method limits the boundary of the scene set by limiting the fluctuation range of the uncertain scenes of each node, and the value in all transition scenes can be covered in the boundary range. The proposed model aims to find a group of thermal power generating unit flexibility modification schemes, so that the feasible power flow of a power grid in a wind power fluctuation range in multiple operation days can be met, and the optimal economy is achieved.
The expressions (4) to (6) and (9) to (11) in the original model can be converted into interval forms, and the expression forms and the corresponding scene reduction methods are as follows.
1) System operating costs
Figure BDA0003661442030000077
In the formula: c D The system operating cost;
Figure BDA0003661442030000078
and
Figure BDA0003661442030000079
respectively representing the lower limit and the upper limit of the output of the thermal power generating unit i in a typical day period t in a month m;
Figure BDA00036614420300000710
and
Figure BDA00036614420300000711
respectively is the lower limit and the upper limit of the wind curtailment quantity of the wind turbine generator j in the typical day period t in the month m.
The upper and lower limits of the annual running cost interval are shown in the formula (15) and the formula (16), respectively.
Figure BDA0003661442030000081
Figure BDA0003661442030000082
In the formula:C D and
Figure BDA0003661442030000083
respectively, the lower limit and the upper limit of the system operating cost.
2) Power balance equation
Figure BDA0003661442030000084
The power balance equation in the system consists of 4 parts, namely thermal power unit output, wind power unit output, load and air abandon amount. Because the output and the load power of the wind turbine generator are both known quantities, the known quantities can be transferred to the equal sign right side of an equation, and the output range of unknown quantities is reflected.
Figure BDA0003661442030000085
By respectively combining the upper limit value and the lower limit value of the 2 intervals on the left side of the equal sign of the equation (18), 2 limit scenes can be obtained, and the interval equation can be converted into the following power balance equation under the 2 limit scenes.
Figure BDA0003661442030000086
Figure BDA0003661442030000087
3) Thermal power unit output constraint
Figure BDA0003661442030000088
Equation (21) can be converted into the following 2 inequalities by strict constraints on the upper and lower limits of the interval inequality.
Figure BDA0003661442030000089
Figure BDA00036614420300000810
4) Climbing restraint
Figure BDA00036614420300000811
By restricting the upper and lower limits of these 2 intervals, equation (24) can be converted into the following 2 inequalities.
Figure BDA00036614420300000812
Figure BDA00036614420300000813
5) Line flow constraint
Figure BDA0003661442030000091
Transferring the known quantity and the unknown quantity of the interval inequality to two sides of the inequality respectively, and converting the two-sided inequality into 2 single-sided inequalities to solve conveniently, thus obtaining:
Figure BDA0003661442030000092
Figure BDA0003661442030000093
for the known quantity on the right side of the equation, the interval boundary is shrunk according to a low boundary theory, the minimum quantity on the right side of the inequality is selected as a shrinking condition of the boundary, the condition is the tightest constraint added to the line power flow, and the power flow under all possible scenes can be guaranteed to be feasible. And for the unknown quantity on the left side of the equation, carrying out constraint according to the upper limit and the lower limit of the interval. Through the above process, the 2 inequalities can be converted as follows.
Figure BDA0003661442030000094
Figure BDA0003661442030000095
6) Wind curtailment
Figure BDA0003661442030000096
For the interval inequality, it can be ensured that all scenes can meet the requirements as long as the upper limit and the lower limit of the 2 intervals are respectively constrained. This inequality can be converted into the following 2 constraints.
Figure BDA0003661442030000097
Figure BDA0003661442030000098
In summary, the constraint conditions and the objective function in the original model may be expressed in an interval form, and full scene reduction may be performed through rigorous mathematical derivation, so as to ensure that each interval range covers all possible transition scenes. The protomodel can be expressed in the following general form:
Figure BDA0003661442030000101
in the formula: w, W,P w
Figure BDA0003661442030000102
Wind power output, a wind power output lower limit and a wind power output upper limit vector are respectively;Xand
Figure BDA0003661442030000103
the lower and upper limits of X, respectively. The optimization target of the proposed model comprises 2 parts in total, wherein one part is the year of the thermal power generating unitThe flexibility modification cost is a single-value target; the other part is an annual operation cost interval formed under a given wind power output prediction interval, and is an interval target. For ease of calculation, the cost of combining these 2 parts results in an annual combined cost interval. In order to comprehensively measure the overall characteristics of the target interval, the interval mean value and the deviation value are used as the characteristic quantity of the interval to carry out collaborative optimization. The master model can be converted into:
Figure BDA0003661442030000104
in the formula: a (C (X, U, W)) and J (C (X, U, W)) are the average value and the deviation value of the target section C (X, U, W), respectively;C(X, U, W) and
Figure BDA0003661442030000105
(X, U, W) are the lower and upper limits of the target interval C (X, U, W), respectively.
3. Solving interval optimization model by GSOMP algorithm
The GSOMP algorithm is a multi-objective optimization method based on group search, and the method simulates a swarm intelligent search technology in the optimization process, so that the GSOMP algorithm has good performance in the application of searching a non-dominated solution set. And (3) solving a group of non-dominated optimal solution sets with different sub-targets in focus by adopting a GSOMP algorithm, wherein the solution sets comprise a plurality of candidate schemes with different annual comprehensive cost interval characteristics.
4. Target interval evaluation method based on evidence reasoning for evaluation
An optimal scheme is selected from a plurality of candidate schemes by adopting an evidence reasoning method, the evidence reasoning method is a three-layer heuristic evaluation method based on attributes, and the basic attributes to be evaluated, namely the mean value and the deviation value of the annual comprehensive cost interval, are defined at the bottom layer of an evaluation frame. In the middle layer of the evaluation framework, 5 evaluation levels are set, wherein the evaluation levels are excellent, good, general, marginal and poor, the layer divides each basic attribute into 5 evaluation levels and generates corresponding level confidence, and comprehensive evaluation results of each scheme on the 5 evaluation levels can be obtained. And combining the evaluation results of all the grades and the corresponding grade demand degrees at the top layer of the evaluation framework to obtain the comprehensive evaluation result of the scheme.
The model in the embodiment can effectively improve the comprehensive economy and the trend adaptability of the scheme, the model adopts a full-scene optimization method, the running cost under the scene with large wind power output is obviously improved, the low running cost under the conventional wind power output scene is maintained, and the comprehensive economy of the scheme is improved. In the aspect of flow adaptability, the full-scene optimization mode can better adapt to the flow uncertainty under the high wind power access proportion, effectively relieves the flow blocking phenomenon and reduces the air abandonment rate under each uncertain scene, thereby obtaining a group of thermal power generating unit flexibility adjusting schemes with optimal comprehensive economy.
The present invention is not limited to the above embodiments, and any changes or substitutions that can be easily made by those skilled in the art within the technical scope of the present invention are also within the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (3)

1. The thermal power generating unit flexibility adjusting and scheduling method based on interval optimization is characterized by comprising the following steps of:
step 1: comprehensively considering flexibility adjusting cost and operating cost of the thermal power generating unit, and establishing a flexibility adjusting certainty model of the thermal power generating unit by taking annual comprehensive cost as an optimization target;
and 2, step: the method comprises the steps of limiting the boundary of a scene set by limiting the fluctuation range of uncertain scenes of each node, and establishing a thermal power generating unit flexibility adjusting interval optimization model;
and step 3: solving the thermal power generating unit flexibility adjusting interval optimization model established in the step 2 by using a GSOMP algorithm; using a GSOMP algorithm to cooperatively optimize the mean value and the deviation value of the target interval and generating a candidate scheme with the characteristics of the comprehensive cost intervals in different years;
and 4, step 4: and (3) evaluating the annual comprehensive cost intervals in the candidate schemes in the step (3) by using a multi-attribute decision method, and selecting a group of optimal schemes according to the comprehensive performance of the characteristic quantities of each interval.
2. The thermal power generating unit flexibility regulation scheduling method based on interval optimization according to claim 1, wherein the thermal power generating unit flexibility regulation certainty model in the step 1 is specifically as follows:
an objective function: the minimum annual comprehensive cost is taken as an optimization target
min(C B +C O ) (1)
In the formula, C B And C O Respectively adjusting the annual flexibility adjusting cost and the annual operation cost of the thermal power generating unit;
the annual transformation cost of the thermal power generating unit is as follows (2):
Figure FDA0003661442020000011
in the formula, P i G,min And P i G,min0 Respectively representing the lowest output of the thermal power generating unit i before and after modification; c y_rebuild Cost for annual flexibility modification per unit volume; m G The number of the thermoelectric generator sets in the system;
the annual flexibility transformation cost per unit capacity is converted into the following formula (3):
Figure FDA0003661442020000012
in the formula, C rebuild Cost for flexibility modification per unit capacity; y is the unit life; eta is annual interest rate;
the middle-aged operating cost of the system is as follows:
Figure FDA0003661442020000021
in the formula, D m Days of month m; t is the number of hours of operation in each typical day;
Figure FDA0003661442020000022
the output of the thermal power generating unit i in a typical day period t in a month m is obtained; lambda [ alpha ] w Punishment cost of unit air volume abandon; m is a group of W The number of wind generating sets in the system;
Figure FDA0003661442020000023
the wind power is the wind curtailment quantity of the wind turbine generator j in the typical day period t in the month m;
and power balance constraint:
Figure FDA0003661442020000024
in the formula (I), the compound is shown in the specification,
Figure FDA0003661442020000025
the output of the wind turbine generator j in a typical day time t in a month m refers to the predicted output of the wind turbine generator;
Figure FDA0003661442020000026
for load, node k is active power at a typical day period t in month m; m D Is the number of load nodes in the system;
output restraint of the thermal power generating unit:
Figure FDA0003661442020000027
in the formula, P i G,max The maximum output of the thermal power generating unit i is obtained;
flexibility modification constraint of the thermal power generating unit:
Figure FDA0003661442020000028
Figure FDA0003661442020000029
in the formula, P i G,min1 And P i G,min2 Respectively carrying out non-oil-feeding peak regulation transformation and oil-feeding peak regulation transformation on the thermal power generating unit i to obtain the lowest output of the thermal power generating unit; x is the number of i,0 、x i,1 、x i,2 The variable of the transformation scheme decision of the unit i is 0-1 variable, x i,0 =1、x i,1 =1、x i,2 When the unit i is 1, the unit i is not transformed, is not transformed in an oil feeding mode, and is transformed in an oil feeding mode;
and (3) climbing restraint of the thermal power generating unit:
Figure FDA00036614420200000210
in the formula (I), the compound is shown in the specification,
Figure FDA00036614420200000211
and
Figure FDA00036614420200000212
the maximum upward climbing speed and the maximum downward climbing speed of the thermal power generating unit i are respectively set;
and (3) line power flow constraint:
Figure FDA0003661442020000031
in the formula, G n-i 、G n-j 、G n-k Power transfer factors of a thermal power generating unit i, a wind power generating unit j and a load node k to a line n are respectively set; m L Is the number of lines in the system;
Figure FDA0003661442020000032
is the maximum allowed transmission power of line n;
abandoned air volume constraint of wind turbine generator
Figure FDA0003661442020000033
3. The thermal power generating unit flexibility adjusting and scheduling method based on interval optimization according to claim 1, wherein the thermal power generating unit flexibility adjusting interval optimization model in the step 2 is specifically as follows:
the system operating cost is as follows:
Figure FDA0003661442020000034
in the formula: c D The system operating cost;
Figure FDA0003661442020000035
and
Figure FDA0003661442020000036
respectively representing the lower limit and the upper limit of the output of the thermal power generating unit i in a typical day period t in a month m;
Figure FDA0003661442020000037
and
Figure FDA0003661442020000038
respectively setting a lower limit and an upper limit of the wind curtailment quantity of the wind turbine generator j in a typical day period t in a month m;
lower limit and upper limit of annual running cost interval:
Figure FDA0003661442020000039
Figure FDA00036614420200000310
in the formula:C D and
Figure FDA00036614420200000311
respectively the lower limit and the upper limit of the system operation cost;
power balance equation:
Figure FDA00036614420200000312
output restraint of the thermal power generating unit:
Figure FDA00036614420200000313
and (3) climbing restraint:
Figure FDA00036614420200000314
and (3) line power flow constraint:
Figure FDA0003661442020000041
wind abandon restraint:
Figure FDA0003661442020000042
CN202210574276.5A 2022-05-25 2022-05-25 Thermal power generating unit flexibility adjusting and scheduling method based on interval optimization Pending CN114884134A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210574276.5A CN114884134A (en) 2022-05-25 2022-05-25 Thermal power generating unit flexibility adjusting and scheduling method based on interval optimization

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210574276.5A CN114884134A (en) 2022-05-25 2022-05-25 Thermal power generating unit flexibility adjusting and scheduling method based on interval optimization

Publications (1)

Publication Number Publication Date
CN114884134A true CN114884134A (en) 2022-08-09

Family

ID=82677091

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210574276.5A Pending CN114884134A (en) 2022-05-25 2022-05-25 Thermal power generating unit flexibility adjusting and scheduling method based on interval optimization

Country Status (1)

Country Link
CN (1) CN114884134A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116454890A (en) * 2023-04-20 2023-07-18 中国南方电网有限责任公司 Combined control method, device and equipment for unit based on SCUC model
CN117744410A (en) * 2024-02-19 2024-03-22 中汽研汽车检验中心(天津)有限公司 Interval multi-objective optimization method, equipment and medium based on vehicle stability requirement

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116454890A (en) * 2023-04-20 2023-07-18 中国南方电网有限责任公司 Combined control method, device and equipment for unit based on SCUC model
CN116454890B (en) * 2023-04-20 2024-02-06 中国南方电网有限责任公司 Combined control method, device and equipment for unit based on SCUC model
CN117744410A (en) * 2024-02-19 2024-03-22 中汽研汽车检验中心(天津)有限公司 Interval multi-objective optimization method, equipment and medium based on vehicle stability requirement
CN117744410B (en) * 2024-02-19 2024-05-07 中汽研汽车检验中心(天津)有限公司 Interval multi-objective optimization method, equipment and medium based on vehicle stability requirement

Similar Documents

Publication Publication Date Title
JP7261507B2 (en) Electric heat pump - regulation method and system for optimizing cogeneration systems
CN109256799B (en) New energy power system optimal scheduling method based on sample entropy
CN114884134A (en) Thermal power generating unit flexibility adjusting and scheduling method based on interval optimization
Yang et al. Hybrid policy-based reinforcement learning of adaptive energy management for the Energy transmission-constrained island group
CN110620397B (en) Peak regulation balance evaluation method for high-proportion renewable energy power system
CN112701687B (en) Robust optimization operation method of gas-electricity distribution network system considering price type combined demand response
CN112564109B (en) Frequency modulation optimization operation method based on participation of energy storage system in large-scale offshore wind power
CN109711706A (en) Consider the active distribution network substation planning method of distributed generation resource and demand response
CN115081693A (en) Two-stage capacity configuration method of water-wind-light complementary system with uncertain energy
CN109886446A (en) Based on the Electrical Power System Dynamic economic load dispatching method for improving Chaos particle swarm optimization algorithm
CN110796373A (en) Wind power consumption-oriented multi-stage scene generation electric heating system optimal scheduling method
CN115759610A (en) Multi-target planning method for source-grid and storage cooperation of power system and application thereof
CN116523240A (en) Micro-grid evaluation feedback scheduling method based on MCMC method prediction
CN115051388A (en) Distribution robustness-based 'source-network-load-storage' two-stage scheduling optimization method
CN110991927B (en) Power supply planning method for improving intermittent power supply complementation effect of different areas of regional power grid
CN110992206B (en) Optimal scheduling method and system for multi-source electric field
CN116341881B (en) Robust advanced scheduling method and system for electric-thermal system considering flexibility of heat supply network
CN116468168A (en) Distributed power supply multi-target hierarchical planning method based on improved beluga optimization algorithm
CN115758763A (en) Multi-energy flow system optimal configuration method and system considering source load uncertainty
CN115860169A (en) Multi-objective optimization planning method and system for deep peak regulation transformation of thermal power generating unit
CN114398777A (en) Power system flexibility resource allocation method based on Bashi game theory
CN114400652A (en) Multi-energy power generation optimization scheduling method considering active participation of nuclear power in peak shaving
CN113807566A (en) Community energy scheduling method based on renewable energy and uncertain internal load
CN113642165B (en) Distributed optimization method of electric-gas coupling distribution network considering uncertainty of demand side
CN107633330A (en) System and method are coordinated in the multipotency interconnection that optimization is searched based on tree

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20231208

Address after: 102206 No. 2 Nong Road, Changping District, Beijing

Applicant after: NORTH CHINA ELECTRIC POWER University

Applicant after: HUANENG NANJING JINLING POWER GENERATION Co.,Ltd.

Applicant after: HUANENG GROUP TECHNOLOGY INNOVATION CENTER Co.,Ltd.

Applicant after: Xi'an Thermal Power Research Institute Co.,Ltd.

Address before: 102206 No. 2 Nong Road, Changping District, Beijing

Applicant before: NORTH CHINA ELECTRIC POWER University

Applicant before: HUANENG GROUP TECHNOLOGY INNOVATION CENTER Co.,Ltd.

Applicant before: Xi'an Thermal Power Research Institute Co.,Ltd.