CN111799842A - Multi-stage power transmission network planning method and system considering flexibility of thermal power generating unit - Google Patents

Multi-stage power transmission network planning method and system considering flexibility of thermal power generating unit Download PDF

Info

Publication number
CN111799842A
CN111799842A CN202010758007.5A CN202010758007A CN111799842A CN 111799842 A CN111799842 A CN 111799842A CN 202010758007 A CN202010758007 A CN 202010758007A CN 111799842 A CN111799842 A CN 111799842A
Authority
CN
China
Prior art keywords
generating unit
power
year
planning
transmission line
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.)
Granted
Application number
CN202010758007.5A
Other languages
Chinese (zh)
Other versions
CN111799842B (en
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.)
State Grid Corp of China SGCC
Economic and Technological Research Institute of State Grid Shandong Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Economic and Technological Research Institute of State Grid Shandong Electric Power 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 State Grid Corp of China SGCC, Economic and Technological Research Institute of State Grid Shandong Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202010758007.5A priority Critical patent/CN111799842B/en
Publication of CN111799842A publication Critical patent/CN111799842A/en
Application granted granted Critical
Publication of CN111799842B publication Critical patent/CN111799842B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/004—Generation forecast, e.g. methods or systems for forecasting future energy generation
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/76—Power conversion electric or electronic aspects
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • Y—GENERAL 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

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

Abstract

The invention provides a multi-stage power transmission network planning method and system considering flexibility of a thermal power generating unit, and the method and system are based on time value of capital and time sequence characteristics of load, wind power and the like, a multi-stage power transmission network planning optimization model considering flexibility modification of the thermal power generating unit is constructed, mixed integer nonlinear constraint conditions in the optimization model are converted into mixed integer linear constraint conditions, and the converted mixed integer linear planning model is solved by adopting a mixed integer linear planning method, so that a multi-stage joint decision scheme of flexibility modification of the thermal power generating unit and power transmission network planning is obtained, and the economy of planning operation investment of a power system is improved.

Description

Multi-stage power transmission network planning method and system considering flexibility of thermal power generating unit
Technical Field
The invention relates to the technical field of power grid planning, in particular to a multi-stage power transmission grid planning method and system considering flexibility of a thermal power generating unit.
Background
The rapid development of new energy is a necessary choice for solving the problems of shortage of fossil energy and environmental pollution at present. In the new energy power generation technology, the wind power generation technology is the most mature and has the most large-scale development value, and the wind power generation in the world is in a large-scale development situation. Under the influence of natural factors, the wind power generation has obvious anti-peak-shaving characteristics, namely, the wind power output is low in the daytime peak load period and is high in the nighttime valley load period. Large-scale wind power integration also brings great challenges to the planning and operation of power systems.
The output power of the thermal power generating unit can be adjusted within a certain range, and the thermal power generating unit plays an important role in balancing large-scale wind power change and promoting wind power consumption. The flexibility of the thermal power generating unit is embodied in that the system power unbalance amount caused by load change or wind power change can be tracked at a certain rate, and when the wind power output power is increased, the output power of the thermal power generating unit is reduced to allow clean energy to generate electricity so as to promote wind power consumption. In China, the thermal power generating unit generally has the problems of abundant total amount and insufficient flexibility, and the wind power grid-connected consumption is severely limited. In order to solve the problem, the flexibility modification of the thermal power generating units needs to be implemented, but how to make a scientific and reasonable flexibility modification time sequence of the thermal power generating units is worthy of deep research.
The goal of power grid planning is to meet power transmission requirements by making decisions on the construction of transmission lines, given load predictions and power supply layouts. In the traditional power grid planning, a target grid structure is used as a research object for optimization decision making, and the time sequence characteristics of the power grid structure are difficult to reflect. The power grid planning is a long-term and rolling decision process, the transition of schemes among multiple planning stages and the restriction of decision variables on the value of each stage need to be considered, and the time value of capital needs to be considered in the planning process to carry out multi-period coordinated dynamic planning. At present, the simulation of the operation situation of each planning stage in the multi-stage power grid planning decision usually selects an operation mode corresponding to the peak load, and the change scenes of the load and the renewable energy power in each planning stage are less involved.
Disclosure of Invention
The invention aims to provide a multi-stage power transmission network planning method and system considering the flexibility of a thermal power generating unit, aims to solve the problems of lack of load and wind power consideration in multi-stage power network planning decision-making in the prior art, realizes a multi-stage combined decision-making scheme, and improves the economy of planning operation investment of a power system.
In order to achieve the technical purpose, the invention provides a multi-stage power transmission network planning method considering the flexibility of a thermal power generating unit, which comprises the following operations:
constructing a multi-stage power transmission network planning optimization model considering flexibility transformation of the thermal power generating unit, wherein the optimization model takes the minimum sum of the investment cost and the present value of the operation cost in the whole planning period as a target and comprises a plurality of constraint conditions;
and converting the mixed integer nonlinear constraint condition in the optimization model into a mixed integer linear constraint condition, and solving the converted mixed integer linear programming model by using a mixed integer linear programming algorithm to obtain a final thermal power unit flexibility modification and transmission grid multi-stage combined programming scheme.
Preferably, the method further comprises the steps of inputting parameters of a traditional thermal power unit and parameters of a power transmission element of the power grid in the current situation before the model is built, inputting maximum load, maximum wind power, 24-hour load power and wind power change per month of each year in a planning period, setting a range of the thermal power unit which can participate in flexibility modification and a modification cost according to the operation condition of the thermal power unit, setting a range of a candidate new power transmission line and a construction cost according to the corridor condition of the power transmission line, and setting a wind power curtailment cost coefficient and a system allowable curtailment rate.
Preferably, the objective function expression of the optimization model is:
Figure BDA0002612203900000031
wherein nw is the total number of the wind power plants; ny is the total year of the planning cycle; ng is the number of conventional thermal power generating units; nl is the number of transmission line corridors planned and constructed; alpha is alphawRepresenting wind curtailment penalty cost;
Figure BDA0002612203900000032
the electric power curtailment of the wind power plant w is the typical day h period of the mth year, mth month in the planning cycle; beta is agRepresenting the flexibility modification cost of the thermal power generating unit g;
Figure BDA0002612203900000033
the binary variable represents whether the thermal power generating unit g is subjected to flexible modification in the y year in the planning period,
Figure BDA0002612203900000034
the thermal power generating unit g is flexibly modified in the y year in the planning period,
Figure BDA0002612203900000035
showing that the thermal power generating unit g does not perform flexible modification in the y year in the planning period; clInvestment cost for construction of a single-circuit power transmission line on a power transmission corridor l;
Figure BDA0002612203900000036
representing the number of newly added power transmission lines in the y year on the power transmission line corridor l; ρ is the discount rate.
Preferably, the constraint condition includes:
1) the upper and lower limits of active power of a conventional thermal power generating unit are constrained:
Figure BDA00026122039000000314
Figure BDA00026122039000000315
Figure BDA00026122039000000316
wherein,
Figure BDA0002612203900000037
and
Figure BDA0002612203900000038
the active power upper and lower limits of the thermal power generating unit g are respectively set;
Figure BDA0002612203900000039
the active power of the thermal power generating unit g is the typical day h period of the mth year, mth month in the planning cycle;
Figure BDA00026122039000000310
representing the increased deep peak regulation power of the thermal power generating unit through flexible modification until the y year; delta PgThe method comprises the steps that the depth peak regulation power which can be increased by flexible modification of a thermal power generating unit g is shown;
Figure BDA00026122039000000311
the binary variable represents whether the thermal power generating unit g is subjected to flexible modification in the y year in the planning period,
Figure BDA00026122039000000312
the thermal power generating unit g is flexibly modified in the y year in the planning period,
Figure BDA00026122039000000313
representing that the thermal power generating unit g is not subjected to flexibility modification in the y year in a planning period; zgIs a binary variable which indicates whether the live generator g is flexibly modified in a planning period, Zg1 denotes the flexible modification of the generator set g in the planning cycle, ZgWhen the number is 0, the live-wire generator g is not flexibly transformed in a planning period;
2) conventional thermal power generating unit climbing restraint:
Figure BDA0002612203900000041
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgThe method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure BDA0002612203900000042
the active power of the thermal power generating unit g is the active power of the thermal power generating unit g in a typical day h-1 of mth month in the planned cycle;
3) node power balance constraint:
Figure BDA0002612203900000043
wherein, Pl y,m,hThe transmission active power of a power transmission line corridor l in a typical day h period of a month m in the y year in a planning cycle;
Figure BDA0002612203900000044
the active power predicted value of the wind power plant w in a typical day h period of the mth year m month in the planning cycle;
Figure BDA0002612203900000045
the active power predicted value of the load d in a typical day h period of the mth year m month in the planning cycle; siAnd EiThe total number of the power transmission lines with the node i as a head node and the tail end node respectively; gi、WiAnd DiRespectively representing the total number of thermal power generating units, wind power plants and loads on the node i; nb is the number of the nodes of the power grid;
4) and (5) abandoned wind power constraint:
Figure BDA0002612203900000046
Figure BDA0002612203900000047
wherein, γyRepresenting a set wind power curtailment rate threshold allowed in the y year;
5) transmission capacity constraint of the transmission line:
Figure BDA0002612203900000051
Figure BDA0002612203900000052
wherein, BlThe susceptance of a single-circuit power transmission line on a power transmission line corridor l; thetai y,m,hA node i voltage phase angle is a typical day h period of a mth month in the y year in a planning cycle;
Figure BDA0002612203900000053
the number of the existing transmission lines on the transmission line corridor l is shown;
Figure BDA0002612203900000054
representing the number of newly added transmission lines in the t year on the transmission line corridor l; pl y,m,hThe transmission active power of a power transmission line corridor l in a typical day h period of a month m in the y year in a planning cycle; pl maxThe transmission capacity of a single-circuit transmission line on a transmission line corridor l;
6) the maximum extensible line number constraint of the power transmission line corridor is as follows:
Figure BDA0002612203900000055
wherein Z isl maxThe maximum number of the extensible power transmission lines is 1;
7) node voltage phase angle constraint:
Figure BDA0002612203900000056
wherein,
Figure BDA0002612203900000057
a node i voltage phase angle is a typical day h period of a mth month in the y year in a planning cycle;
8) n-1, the upper and lower limits of active power of the thermal power generating unit are constrained under the condition of an expected accident:
Figure BDA0002612203900000058
wherein,
Figure BDA0002612203900000059
active power of the thermal power generating unit g under an expected accident k in a typical day h period of m month in the y year in a planning cycle, wherein nk is the total number of the expected accidents of the power transmission line N-1;
9) n-1 conventional thermal power generating unit climbing restraint under the condition of anticipated accidents:
Figure BDA00026122039000000510
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgThe method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure BDA0002612203900000061
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h period of mth month in the y year in a planning cycle;
Figure BDA0002612203900000062
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h-1 of mth month in the y year in a planning cycle;
10) n-1 node power balance constraint under the condition of an expected accident:
Figure BDA0002612203900000063
wherein, Pl y,m,h,kThe transmission active power of the transmission line corridor l under the expected accident k is predicted for the typical day h period of the mth month in the y year in the planning cycle;
Figure BDA0002612203900000064
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h period of mth month in the y year in a planning cycle;
Figure BDA0002612203900000065
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth year, mth month in the planning cycle;
11) n-1 abandon wind power restraint under the condition of anticipated accidents:
Figure BDA0002612203900000066
wherein,
Figure BDA0002612203900000067
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth year, mth month in the planning cycle;
12) n-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure BDA0002612203900000068
Figure BDA0002612203900000069
wherein nl is the total number of the transmission line corridors; b islThe susceptance of a single-circuit power transmission line on a power transmission line corridor l; ny is the total year of the planning cycle;
Figure BDA00026122039000000610
the number of the existing transmission lines on the transmission line corridor l is shown;
Figure BDA0002612203900000071
representing the number of newly added transmission lines in the t year on the transmission line corridor l;
Figure BDA0002612203900000072
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure BDA0002612203900000073
to representIt is expected that the wire on the k power transmission line corridor is shut down under the accident,
Figure BDA0002612203900000074
representing that no line is stopped on the transmission line corridor under the expected accident k; pl y,m,h,kThe transmission active power of the transmission line corridor l under the expected accident k is predicted for the typical day h period of the mth month in the y year in the planning cycle;
Figure BDA0002612203900000075
the transmission capacity of a single-circuit transmission line on a transmission line corridor l;
13) n-1 node voltage phase angle constraint under the expected accident condition:
Figure BDA0002612203900000076
wherein,
Figure BDA0002612203900000077
and forecasting the voltage phase angle of the node i under the accident k for the typical day h period of the mth year m month in the planning cycle.
Preferably, the non-linear constraint expression in the transmission capacity constraint of the power transmission line is converted into an equivalent linear expression:
Figure BDA0002612203900000078
Figure BDA0002612203900000079
wherein M is a constant.
Preferably, the nonlinear constraint expression in the transmission capacity constraint of the transmission line under the N-1 expected accident condition is converted into an equivalent linear expression, that is:
Figure BDA00026122039000000710
Figure BDA00026122039000000711
the invention also provides a multi-stage power transmission network planning system considering the flexibility of the thermal power generating unit, which comprises the following components:
the planning model construction module is used for constructing a multi-stage power transmission network planning optimization model considering flexibility transformation of the thermal power generating unit, the optimization model takes the minimum sum of the investment cost and the present value of the operation cost in the whole planning period as a target and comprises a plurality of constraint conditions;
and the model solving module is used for converting the mixed integer nonlinear constraint condition in the optimization model into a mixed integer linear constraint condition, and solving the converted mixed integer linear programming model by using a mixed integer linear programming algorithm to obtain a final thermal power unit flexibility modification and power transmission network multi-stage combined programming scheme.
Preferably, the system further comprises a model parameter setting module, which is used for inputting parameters of a traditional thermal power unit and parameters of a power transmission element of the power grid in the current situation before the model is built, inputting maximum load, maximum wind power, 24-hour load power and wind power change per month in a planning period, setting a range of the thermal power unit which can participate in flexibility modification and a modification cost according to the operation condition of the thermal power unit, setting a range of a candidate newly-built power transmission line and a construction cost according to the corridor condition of the power transmission line, and setting a wind power curtailment cost coefficient and a power curtailment rate allowed by the system.
Preferably, the objective function expression of the optimization model is:
Figure BDA0002612203900000081
wherein nw is the total number of the wind power plants; ny is the total year of the planning cycle; ng is the number of conventional thermal power generating units; nl is the number of transmission line corridors planned and constructed; alpha is alphawRepresenting wind curtailment penalty cost;
Figure BDA0002612203900000082
the electric power curtailment of the wind power plant w is the typical day h period of the mth year, mth month in the planning cycle; beta is agRepresenting the flexibility modification cost of the thermal power generating unit g;
Figure BDA0002612203900000083
the binary variable represents whether the thermal power generating unit g is subjected to flexible modification in the y year in the planning period,
Figure BDA0002612203900000084
the thermal power generating unit g is flexibly modified in the y year in the planning period,
Figure BDA0002612203900000085
showing that the thermal power generating unit g does not perform flexible modification in the y year in the planning period; clInvestment cost for construction of a single-circuit power transmission line on a power transmission corridor l;
Figure BDA0002612203900000091
representing the number of newly added power transmission lines in the y year on the power transmission line corridor l; ρ is the discount rate.
Preferably, the constraint condition includes:
1) the upper and lower limits of active power of a conventional thermal power generating unit are constrained:
Figure BDA0002612203900000092
Figure BDA0002612203900000093
Figure BDA0002612203900000094
wherein,
Figure BDA0002612203900000095
and
Figure BDA0002612203900000096
respectively of thermal power generating units gAn upper and lower active power limit;
Figure BDA0002612203900000097
the active power of the thermal power generating unit g is the typical day h period of the mth year, mth month in the planning cycle;
Figure BDA0002612203900000098
representing the increased deep peak regulation power of the thermal power generating unit through flexible modification until the y year; delta PgThe method comprises the steps that the depth peak regulation power which can be increased by flexible modification of a thermal power generating unit g is shown;
Figure BDA0002612203900000099
the binary variable represents whether the thermal power generating unit g is subjected to flexible modification in the y year in the planning period,
Figure BDA00026122039000000910
the thermal power generating unit g is flexibly modified in the y year in the planning period,
Figure BDA00026122039000000911
representing that the thermal power generating unit g is not subjected to flexibility modification in the y year in a planning period; zgIs a binary variable which indicates whether the live generator g is flexibly modified in a planning period, Zg1 denotes the flexible modification of the generator set g in the planning cycle, ZgWhen the number is 0, the live-wire generator g is not flexibly transformed in a planning period;
2) conventional thermal power generating unit climbing restraint:
Figure BDA00026122039000000912
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgThe method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure BDA00026122039000000913
the active power of the thermal power generating unit g is the active power of the thermal power generating unit g in a typical day h-1 of mth month in the planned cycle;
3) node power balance constraint:
Figure BDA0002612203900000101
wherein, Pl y,m,hThe transmission active power of a power transmission line corridor l in a typical day h period of a month m in the y year in a planning cycle;
Figure BDA0002612203900000102
the active power predicted value of the wind power plant w in a typical day h period of the mth year m month in the planning cycle;
Figure BDA0002612203900000103
the active power predicted value of the load d in a typical day h period of the mth year m month in the planning cycle; siAnd EiThe total number of the power transmission lines with the node i as a head node and the tail end node respectively; gi、WiAnd DiRespectively representing the total number of thermal power generating units, wind power plants and loads on the node i; nb is the number of the nodes of the power grid;
4) and (5) abandoned wind power constraint:
Figure BDA0002612203900000104
Figure BDA0002612203900000105
wherein, γyRepresenting a set wind power curtailment rate threshold allowed in the y year;
5) transmission capacity constraint of the transmission line:
Figure BDA0002612203900000106
Figure BDA0002612203900000107
wherein, BlFor power transmission lineSusceptance of a single-circuit power transmission line on the corridor l;
Figure BDA0002612203900000108
a node i voltage phase angle is a typical day h period of a mth month in the y year in a planning cycle;
Figure BDA0002612203900000109
the number of the existing transmission lines on the transmission line corridor l is shown;
Figure BDA00026122039000001010
representing the number of newly added transmission lines in the t year on the transmission line corridor l; pl y,m,hThe transmission active power of a power transmission line corridor l in a typical day h period of a month m in the y year in a planning cycle; pl maxThe transmission capacity of a single-circuit transmission line on a transmission line corridor l;
6) the maximum extensible line number constraint of the power transmission line corridor is as follows:
Figure BDA0002612203900000111
wherein Z isl maxThe maximum number of the extensible power transmission lines is 1;
7) node voltage phase angle constraint:
Figure BDA0002612203900000112
wherein,
Figure BDA0002612203900000113
a node i voltage phase angle is a typical day h period of a mth month in the y year in a planning cycle;
8) n-1, the upper and lower limits of active power of the thermal power generating unit are constrained under the condition of an expected accident:
Figure BDA0002612203900000114
wherein,
Figure BDA0002612203900000115
active power of the thermal power generating unit g under an expected accident k in a typical day h period of m month in the y year in a planning cycle, wherein nk is the total number of the expected accidents of the power transmission line N-1;
9) n-1 conventional thermal power generating unit climbing restraint under the condition of anticipated accidents:
Figure BDA0002612203900000116
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgThe method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure BDA0002612203900000117
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h period of mth month in the y year in a planning cycle;
Figure BDA0002612203900000118
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h-1 of mth month in the y year in a planning cycle;
10) n-1 node power balance constraint under the condition of an expected accident:
Figure BDA0002612203900000119
wherein, Pl y,m,h,kThe transmission active power of the transmission line corridor l under the expected accident k is predicted for the typical day h period of the mth month in the y year in the planning cycle;
Figure BDA00026122039000001110
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h period of mth month in the y year in a planning cycle;
Figure BDA0002612203900000121
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth year, mth month in the planning cycle;
11) n-1 abandon wind power restraint under the condition of anticipated accidents:
Figure BDA0002612203900000122
wherein,
Figure BDA0002612203900000123
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth year, mth month in the planning cycle;
12) n-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure BDA0002612203900000124
Figure BDA0002612203900000125
wherein nl is the total number of the transmission line corridors; b islThe susceptance of a single-circuit power transmission line on a power transmission line corridor l; ny is the total year of the planning cycle;
Figure BDA0002612203900000126
the number of the existing transmission lines on the transmission line corridor l is shown;
Figure BDA0002612203900000127
representing the number of newly added transmission lines in the t year on the transmission line corridor l;
Figure BDA0002612203900000128
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure BDA0002612203900000129
indicating a wire outage in the k power line corridor under an expected accident,
Figure BDA00026122039000001210
representing that no line is stopped on the transmission line corridor under the expected accident k; pl y,m,h,kThe transmission active power of the transmission line corridor l under the expected accident k is predicted for the typical day h period of the mth month in the y year in the planning cycle; pl maxThe transmission capacity of a single-circuit transmission line on a transmission line corridor l;
13) n-1 node voltage phase angle constraint under the expected accident condition:
Figure BDA00026122039000001211
wherein,
Figure BDA00026122039000001212
and forecasting the voltage phase angle of the node i under the accident k for the typical day h period of the mth year m month in the planning cycle.
The effect provided in the summary of the invention is only the effect of the embodiment, not all the effects of the invention, and one of the above technical solutions has the following advantages or beneficial effects:
compared with the prior art, the method is based on the time value of capital, the time sequence characteristics of load, wind power and the like are considered, the multi-stage power transmission network planning optimization model considering thermal power unit flexibility modification is constructed, the mixed integer nonlinear constraint condition in the optimization model is converted into the mixed integer linear constraint condition, the converted mixed integer linear programming model is solved by adopting the mixed integer linear programming method, the multi-stage joint decision scheme of thermal power unit flexibility modification and power transmission network planning is obtained, and the economy of power system planning operation investment is improved.
Drawings
Fig. 1 is a flowchart of a multi-stage power transmission network planning method considering flexibility of a thermal power generating unit according to an embodiment of the present invention;
fig. 2 is a block diagram of a multi-stage power transmission network planning system considering flexibility of a thermal power generating unit according to an embodiment of the present invention.
Detailed Description
In order to clearly explain the technical features of the present invention, the following detailed description of the present invention is provided with reference to the accompanying drawings. The following disclosure provides many different embodiments, or examples, for implementing different features of the invention. To simplify the disclosure of the present invention, the components and arrangements of specific examples are described below. Furthermore, the present invention may repeat reference numerals and/or letters in the various examples. This repetition is for the purpose of simplicity and clarity and does not in itself dictate a relationship between the various embodiments and/or configurations discussed. It should be noted that the components illustrated in the figures are not necessarily drawn to scale. Descriptions of well-known components and processing techniques and procedures are omitted so as to not unnecessarily limit the invention.
The following describes in detail a multi-stage power transmission network planning method and system considering flexibility of a thermal power generating unit, which are provided by the embodiments of the present invention, with reference to the accompanying drawings.
As shown in fig. 1, the invention discloses a multi-stage power transmission network planning method considering flexibility of a thermal power generating unit, which comprises the following operations:
s1, inputting parameters of a traditional thermal power generating unit and parameters of a power transmission element of a current power grid, inputting maximum load, maximum wind power, 24-hour load power and wind power change per typical day of each year in a planning period, setting a thermal power generating unit range capable of participating in flexibility modification and modification cost according to the operation condition of the thermal power generating unit, setting a candidate newly-built power transmission line range and construction cost according to the corridor condition of the power transmission line, and setting a wind power curtailment cost coefficient and a system allowable curtailment rate.
And S2, constructing a multi-stage power transmission network planning optimization model considering flexibility transformation of the thermal power generating unit, wherein the optimization model aims at minimizing the sum of the investment cost and the present value of the operation cost in the whole planning period and comprises a plurality of constraint conditions.
The objective function expression of the optimization model is as follows:
Figure BDA0002612203900000141
wherein nw is the total number of the wind power plants; ny is the total year of the planning cycleCounting; ng is the number of conventional thermal power generating units; nl is the number of transmission line corridors planned and constructed; alpha is alphawRepresenting wind curtailment penalty cost;
Figure BDA0002612203900000142
the electric power curtailment of the wind power plant w is the typical day h period of the mth year, mth month in the planning cycle; beta is agRepresenting the flexibility modification cost of the thermal power generating unit g;
Figure BDA0002612203900000143
the binary variable represents whether the thermal power generating unit g is subjected to flexible modification in the y year in the planning period,
Figure BDA0002612203900000144
the thermal power generating unit g is flexibly modified in the y year in the planning period,
Figure BDA0002612203900000145
showing that the thermal power generating unit g does not perform flexible modification in the y year in the planning period; clInvestment cost for construction of a single-circuit power transmission line on a power transmission corridor l;
Figure BDA0002612203900000146
representing the number of newly added power transmission lines in the y year on the power transmission line corridor l; ρ is the discount rate.
The optimization model includes 13 constraints, which are as follows:
1) the upper and lower limits of active power of a conventional thermal power generating unit are constrained:
Figure BDA0002612203900000151
Figure BDA0002612203900000152
Figure BDA0002612203900000153
wherein,
Figure BDA0002612203900000154
and
Figure BDA0002612203900000155
the active power upper and lower limits of the thermal power generating unit g are respectively set;
Figure BDA0002612203900000156
the active power of the thermal power generating unit g is the typical day h period of the mth year, mth month in the planning cycle;
Figure BDA0002612203900000157
representing the increased deep peak regulation power of the thermal power generating unit through flexible modification until the y year; delta PgThe method comprises the steps that the depth peak regulation power which can be increased by flexible modification of a thermal power generating unit g is shown;
Figure BDA0002612203900000158
the binary variable represents whether the thermal power generating unit g is subjected to flexible modification in the y year in the planning period,
Figure BDA0002612203900000159
the thermal power generating unit g is flexibly modified in the y year in the planning period,
Figure BDA00026122039000001510
representing that the thermal power generating unit g is not subjected to flexibility modification in the y year in a planning period; zgIs a binary variable which indicates whether the live generator g is flexibly modified in a planning period, Zg1 denotes the flexible modification of the generator set g in the planning cycle, ZgAnd (5) indicating that the live generator set g is not flexibly transformed in a planning period as 0.
2) Conventional thermal power generating unit climbing restraint:
Figure BDA00026122039000001511
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgIs a thermal power generating unit gThe climbing rate variable quantity after the flexibility modification;
Figure BDA00026122039000001512
the active power of the thermal power generating unit g is the active power of the thermal power generating unit g in a typical day h-1 of the mth year m and month in the planning period.
3) Node power balance constraint:
Figure BDA0002612203900000161
wherein, Pl y,m,hThe transmission active power of a power transmission line corridor l in a typical day h period of a month m in the y year in a planning cycle;
Figure BDA0002612203900000162
the active power predicted value of the wind power plant w in a typical day h period of the mth year m month in the planning cycle;
Figure BDA0002612203900000163
the active power predicted value of the load d in a typical day h period of the mth year m month in the planning cycle; siAnd EiThe total number of the power transmission lines with the node i as a head node and the tail end node respectively; gi、WiAnd DiRespectively representing the total number of thermal power generating units, wind power plants and loads on the node i; and nb is the number of the nodes of the power grid.
4) And (5) abandoned wind power constraint:
Figure BDA0002612203900000164
Figure BDA0002612203900000165
wherein, γyAnd indicating the set allowable wind power curtailment rate threshold value in the y year.
5) Transmission capacity constraint of the transmission line:
Figure BDA0002612203900000166
Figure BDA0002612203900000167
wherein, BlThe susceptance of a single-circuit power transmission line on a power transmission line corridor l;
Figure BDA0002612203900000168
a node i voltage phase angle is a typical day h period of a mth month in the y year in a planning cycle;
Figure BDA0002612203900000169
the number of the existing transmission lines on the transmission line corridor l is shown;
Figure BDA00026122039000001610
representing the number of newly added transmission lines in the t year on the transmission line corridor l; pl y,m,hThe transmission active power of a power transmission line corridor l in a typical day h period of a month m in the y year in a planning cycle; pl maxThe transmission capacity of a single-circuit transmission line on a transmission line corridor l.
6) The maximum extensible line number constraint of the power transmission line corridor is as follows:
Figure BDA0002612203900000171
wherein Z isl maxAnd the maximum number of the extensible transmission lines is the transmission line corridor l.
7) Node voltage phase angle constraint:
Figure BDA0002612203900000172
wherein,
Figure BDA0002612203900000173
the voltage phase angle of a node i is a typical day h period of a month m in the y year in the planning cycle.
8) N-1, the upper and lower limits of active power of the thermal power generating unit are constrained under the condition of an expected accident:
Figure BDA0002612203900000174
wherein,
Figure BDA0002612203900000175
the active power of the thermal power generating unit g under the expected accident k in the typical day h period of the mth year m month in the planning period is shown, and nk is the total number of the expected accidents of the power transmission line N-1.
9) N-1 conventional thermal power generating unit climbing restraint under the condition of anticipated accidents:
Figure BDA0002612203900000176
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgThe method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure BDA0002612203900000177
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h period of mth month in the y year in a planning cycle;
Figure BDA0002612203900000178
active power of the thermal power generating unit g under the accident k is predicted for a typical day h-1 period of the mth year, mth month in the planning period.
10) N-1 node power balance constraint under the condition of an expected accident:
Figure BDA0002612203900000179
wherein, Pl y,m,h,kThe transmission active power of the transmission line corridor l under the expected accident k is predicted for the typical day h period of the mth month in the y year in the planning cycle;
Figure BDA00026122039000001710
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h period of mth month in the y year in a planning cycle;
Figure BDA00026122039000001711
electric power curtailment of the wind farm w under the accident k is envisioned for a typical day h period of mth year, mth month in the planning cycle.
11) N-1 abandon wind power restraint under the condition of anticipated accidents:
Figure BDA0002612203900000181
wherein,
Figure BDA0002612203900000182
electric power curtailment of the wind farm w under the accident k is envisioned for a typical day h period of mth year, mth month in the planning cycle.
12) N-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure BDA0002612203900000183
Figure BDA0002612203900000184
wherein nl is the total number of the transmission line corridors; b islThe susceptance of a single-circuit power transmission line on a power transmission line corridor l; ny is the total year of the planning cycle;
Figure BDA0002612203900000185
the number of the existing transmission lines on the transmission line corridor l is shown;
Figure BDA0002612203900000186
representing the number of newly added transmission lines in the t year on the transmission line corridor l;
Figure BDA0002612203900000187
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure BDA0002612203900000188
indicating a wire outage in the k power line corridor under an expected accident,
Figure BDA0002612203900000189
representing that no line is stopped on the transmission line corridor under the expected accident k; pl y,m,h,kThe transmission active power of the transmission line corridor l under the expected accident k is predicted for the typical day h period of the mth month in the y year in the planning cycle; pl maxAnd (3) the transmission capacity of a single-circuit transmission line on the transmission line corridor l.
13) N-1 node voltage phase angle constraint under the expected accident condition:
Figure BDA00026122039000001810
wherein,
Figure BDA00026122039000001811
and forecasting the voltage phase angle of the node i under the accident k for the typical day h period of the mth year m month in the planning cycle.
And S3, converting the mixed integer nonlinear constraint condition in the optimization model into a mixed integer linear constraint condition, and solving the converted mixed integer linear programming model by using a mixed integer linear programming algorithm to obtain a final thermal power unit flexibility modification and power transmission network multi-stage joint programming scheme.
The mixed integer nonlinear constraint expression is converted into an equivalent mixed integer linear expression, for example, the nonlinear constraint expression in transmission capacity constraint of the power transmission line is converted into an equivalent linear expression:
Figure BDA0002612203900000191
Figure BDA0002612203900000192
wherein M is a constant.
Similarly, the nonlinear constraint expression in the transmission capacity constraint of the transmission line under the condition of N-1 expected accident is converted into an equivalent linear expression, namely:
Figure BDA0002612203900000193
Figure BDA0002612203900000194
in order to improve the solving calculation efficiency after model conversion, the following auxiliary constraint conditions are introduced:
Figure BDA0002612203900000195
therefore, the optimization model is converted into a mixed integer linear programming model, and the mixed integer linear programming model can be directly used for solving to obtain a final thermal power generating unit flexibility modification and power transmission network multi-stage combined programming scheme.
The method comprises the steps of constructing a multi-stage power transmission network planning optimization model considering flexibility transformation of the thermal power generating unit based on time value of capital and time sequence characteristics of load, wind power and the like, converting mixed integer nonlinear constraint conditions in the optimization model into mixed integer linear constraint conditions, and solving the converted mixed integer linear programming model by adopting a mixed integer linear programming method, so that a multi-stage joint decision scheme of flexibility transformation of the thermal power generating unit and power transmission network planning is obtained, and the economy of planning operation investment of a power system is improved.
As shown in fig. 2, an embodiment of the present invention further discloses a multi-stage power transmission network planning system considering flexibility of a thermal power generating unit, where the system includes:
the planning model construction module is used for constructing a multi-stage power transmission network planning optimization model considering flexibility transformation of the thermal power generating unit, the optimization model takes the minimum sum of the investment cost and the present value of the operation cost in the whole planning period as a target and comprises a plurality of constraint conditions;
and the model solving module is used for converting the mixed integer nonlinear constraint condition in the optimization model into a mixed integer linear constraint condition, and solving the converted mixed integer linear programming model by using a mixed integer linear programming algorithm to obtain a final thermal power unit flexibility modification and power transmission network multi-stage combined programming scheme.
The system further comprises a model parameter setting module, wherein the model parameter setting module is used for inputting parameters of a traditional thermal power generating unit and parameters of a power transmission element of the power grid in the current situation before the model is built, inputting maximum load, maximum wind power, 24-hour load power and wind power change standard values of typical days of each year in a planning period, setting a range of the thermal power generating unit capable of participating in flexibility modification and a modification cost according to the operation condition of the thermal power generating unit, setting a range of a candidate new power transmission line according to the corridor condition of the power transmission line and a construction cost, and setting a wind power curtailment cost coefficient and a power curtailment rate allowed by the system.
The objective function expression of the optimization model is as follows:
Figure BDA0002612203900000201
wherein nw is the total number of the wind power plants; ny is the total year of the planning cycle; ng is the number of conventional thermal power generating units; nl is the number of transmission line corridors planned and constructed; alpha is alphawRepresenting wind curtailment penalty cost;
Figure BDA0002612203900000202
the electric power curtailment of the wind power plant w is the typical day h period of the mth year, mth month in the planning cycle; beta is agRepresenting the flexibility modification cost of the thermal power generating unit g;
Figure BDA0002612203900000211
the binary variable represents whether the thermal power generating unit g is subjected to flexible modification in the y year in the planning period,
Figure BDA0002612203900000212
the thermal power generating unit g is flexibly modified in the y year in the planning period,
Figure BDA0002612203900000213
showing that the thermal power generating unit g does not perform flexible modification in the y year in the planning period; clInvestment cost for construction of a single-circuit power transmission line on a power transmission corridor l;
Figure BDA0002612203900000214
representing the number of newly added power transmission lines in the y year on the power transmission line corridor l; ρ is the discount rate.
The optimization model includes 13 constraints, which are as follows:
1) the upper and lower limits of active power of a conventional thermal power generating unit are constrained:
Figure BDA0002612203900000215
Figure BDA0002612203900000216
Figure BDA0002612203900000217
wherein,
Figure BDA0002612203900000218
and
Figure BDA0002612203900000219
the active power upper and lower limits of the thermal power generating unit g are respectively set;
Figure BDA00026122039000002110
the active power of the thermal power generating unit g is the typical day h period of the mth year, mth month in the planning cycle;
Figure BDA00026122039000002111
representing the increased deep peak regulation power of the thermal power generating unit through flexible modification until the y year; delta PgThe method comprises the steps that the depth peak regulation power which can be increased by flexible modification of a thermal power generating unit g is shown;
Figure BDA00026122039000002112
as binary variablesIndicating whether the thermal power generating unit g carries out flexible modification in the y year in the planning period,
Figure BDA00026122039000002113
the thermal power generating unit g is flexibly modified in the y year in the planning period,
Figure BDA00026122039000002114
representing that the thermal power generating unit g is not subjected to flexibility modification in the y year in a planning period; zgIs a binary variable which indicates whether the live generator g is flexibly modified in a planning period, Zg1 denotes the flexible modification of the generator set g in the planning cycle, ZgAnd (5) indicating that the live generator set g is not flexibly transformed in a planning period as 0.
2) Conventional thermal power generating unit climbing restraint:
Figure BDA0002612203900000221
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgThe method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure BDA0002612203900000222
the active power of the thermal power generating unit g is the active power of the thermal power generating unit g in a typical day h-1 of the mth year m and month in the planning period.
3) Node power balance constraint:
Figure BDA0002612203900000223
wherein, Pl y,m,hThe transmission active power of a power transmission line corridor l in a typical day h period of a month m in the y year in a planning cycle;
Figure BDA0002612203900000224
the active power predicted value of the wind power plant w in a typical day h period of the mth year m month in the planning cycle;
Figure BDA0002612203900000225
the active power predicted value of the load d in a typical day h period of the mth year m month in the planning cycle; siAnd EiThe total number of the power transmission lines with the node i as a head node and the tail end node respectively; gi、WiAnd DiRespectively representing the total number of thermal power generating units, wind power plants and loads on the node i; and nb is the number of the nodes of the power grid.
4) And (5) abandoned wind power constraint:
Figure BDA0002612203900000226
Figure BDA0002612203900000227
wherein, γyAnd indicating the set allowable wind power curtailment rate threshold value in the y year.
5) Transmission capacity constraint of the transmission line:
Figure BDA0002612203900000228
Figure BDA0002612203900000229
wherein, BlThe susceptance of a single-circuit power transmission line on a power transmission line corridor l;
Figure BDA00026122039000002210
a node i voltage phase angle is a typical day h period of a mth month in the y year in a planning cycle;
Figure BDA0002612203900000231
the number of the existing transmission lines on the transmission line corridor l is shown;
Figure BDA0002612203900000232
representing the number of newly added transmission lines in the t year on the transmission line corridor l; pl y,m,hFor planning typical days of mth month in the y year in the periodThe transmission active power of the transmission line corridor l in the h period; pl maxThe transmission capacity of a single-circuit transmission line on a transmission line corridor l.
6) The maximum extensible line number constraint of the power transmission line corridor is as follows:
Figure BDA0002612203900000233
wherein Z isl maxAnd the maximum number of the extensible transmission lines is the transmission line corridor l.
7) Node voltage phase angle constraint:
Figure BDA0002612203900000234
wherein,
Figure BDA0002612203900000235
the voltage phase angle of a node i is a typical day h period of a month m in the y year in the planning cycle.
8) N-1, the upper and lower limits of active power of the thermal power generating unit are constrained under the condition of an expected accident:
Figure BDA0002612203900000236
wherein,
Figure BDA0002612203900000237
the active power of the thermal power generating unit g under the expected accident k in the typical day h period of the mth year m month in the planning period is shown, and nk is the total number of the expected accidents of the power transmission line N-1.
9) N-1 conventional thermal power generating unit climbing restraint under the condition of anticipated accidents:
Figure BDA0002612203900000238
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgThe method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure BDA0002612203900000239
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h period of mth month in the y year in a planning cycle;
Figure BDA00026122039000002310
active power of the thermal power generating unit g under the accident k is predicted for a typical day h-1 period of the mth year, mth month in the planning period.
10) N-1 node power balance constraint under the condition of an expected accident:
Figure BDA0002612203900000241
wherein, Pl y,m,h,kThe transmission active power of the transmission line corridor l under the expected accident k is predicted for the typical day h period of the mth month in the y year in the planning cycle;
Figure BDA0002612203900000242
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h period of mth month in the y year in a planning cycle;
Figure BDA0002612203900000243
electric power curtailment of the wind farm w under the accident k is envisioned for a typical day h period of mth year, mth month in the planning cycle.
11) N-1 abandon wind power restraint under the condition of anticipated accidents:
Figure BDA0002612203900000244
wherein,
Figure BDA0002612203900000245
electric power curtailment of the wind farm w under the accident k is envisioned for a typical day h period of mth year, mth month in the planning cycle.
12) N-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure BDA0002612203900000246
Figure BDA0002612203900000247
wherein,nl is the total number of the power transmission line corridors; b islThe susceptance of a single-circuit power transmission line on a power transmission line corridor l; ny is the total year of the planning cycle;
Figure BDA0002612203900000248
the number of the existing transmission lines on the transmission line corridor l is shown;
Figure BDA0002612203900000249
representing the number of newly added transmission lines in the t year on the transmission line corridor l;
Figure BDA00026122039000002410
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure BDA00026122039000002411
indicating a wire outage in the k power line corridor under an expected accident,
Figure BDA00026122039000002412
representing that no line is stopped on the transmission line corridor under the expected accident k; pl y,m,h,kThe transmission active power of the transmission line corridor l under the expected accident k is predicted for the typical day h period of the mth month in the y year in the planning cycle; pl maxAnd (3) the transmission capacity of a single-circuit transmission line on the transmission line corridor l.
13) N-1 node voltage phase angle constraint under the expected accident condition:
Figure BDA0002612203900000251
wherein,
Figure BDA0002612203900000252
for in the planning periodA voltage phase angle of a node i under an accident k is predicted in a period of m months, typical days and h in the y year.
And converting the mixed integer nonlinear constraint condition in the optimization model into a mixed integer linear constraint condition, and solving the converted mixed integer linear programming model by using a mixed integer linear programming algorithm to obtain a final thermal power unit flexibility modification and transmission grid multi-stage combined programming scheme.
The mixed integer nonlinear constraint expression is converted into an equivalent mixed integer linear expression, for example, the nonlinear constraint expression in transmission capacity constraint of the power transmission line is converted into an equivalent linear expression:
Figure BDA0002612203900000253
Figure BDA0002612203900000254
wherein M is a constant.
Similarly, the nonlinear constraint expression in the transmission capacity constraint of the transmission line under the condition of N-1 expected accident is converted into an equivalent linear expression, namely:
Figure BDA0002612203900000255
Figure BDA0002612203900000256
in order to improve the solving calculation efficiency after model conversion, the following auxiliary constraint conditions are introduced:
Figure BDA0002612203900000257
therefore, the optimization model is converted into a mixed integer linear programming model, and the mixed integer linear programming model can be directly used for solving to obtain a final thermal power generating unit flexibility modification and power transmission network multi-stage combined programming scheme.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting the same, and although the present invention is described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.

Claims (10)

1. A multi-stage power transmission network planning method considering flexibility of a thermal power generating unit is characterized by comprising the following operations:
constructing a multi-stage power transmission network planning optimization model considering flexibility transformation of the thermal power generating unit, wherein the optimization model takes the minimum sum of the investment cost and the present value of the operation cost in the whole planning period as a target and comprises a plurality of constraint conditions;
and converting the mixed integer nonlinear constraint condition in the optimization model into a mixed integer linear constraint condition, and solving the converted mixed integer linear programming model by using a mixed integer linear programming algorithm to obtain a final thermal power unit flexibility modification and transmission grid multi-stage combined programming scheme.
2. The method for planning a multi-stage power transmission network considering the flexibility of the thermal power generating unit according to claim 1, further comprising inputting parameters of a traditional thermal power generating unit and parameters of a power transmission element of a current power grid before constructing the model, inputting maximum load, maximum wind power, and 24-hour load power and wind power change standard values of typical days per month of each year in a planning period, setting a range of the thermal power generating unit capable of participating in the flexibility modification and a modification cost according to the operation condition of the thermal power generating unit, setting a range of a candidate new power transmission line and a construction cost according to the corridor condition of the power transmission line, and setting a wind power curtailment cost coefficient and a system-allowed curtailment rate.
3. The method for multi-stage power transmission network planning with consideration of thermal power generating unit flexibility according to claim 1, wherein an objective function expression of the optimization model is as follows:
Figure FDA0002612203890000011
wherein nw is the total number of the wind power plants; ny is the total year of the planning cycle; ng is the number of conventional thermal power generating units; nl is the number of transmission line corridors planned and constructed; alpha is alphawRepresenting wind curtailment penalty cost;
Figure FDA0002612203890000012
the electric power curtailment of the wind power plant w is the typical day h period of the mth year, mth month in the planning cycle; beta is agRepresenting the flexibility modification cost of the thermal power generating unit g;
Figure FDA0002612203890000013
the binary variable represents whether the thermal power generating unit g is subjected to flexible modification in the y year in the planning period,
Figure FDA0002612203890000021
the thermal power generating unit g is flexibly modified in the y year in the planning period,
Figure FDA0002612203890000022
showing that the thermal power generating unit g does not perform flexible modification in the y year in the planning period; clInvestment cost for construction of a single-circuit power transmission line on a power transmission corridor l;
Figure FDA0002612203890000023
representing the number of newly added power transmission lines in the y year on the power transmission line corridor l; ρ is the discount rate.
4. The method for multi-stage power transmission network planning with consideration of thermal power generating unit flexibility according to claim 1, wherein the constraint condition comprises:
1) the upper and lower limits of active power of a conventional thermal power generating unit are constrained:
Figure FDA0002612203890000024
Figure FDA0002612203890000025
Figure FDA0002612203890000026
wherein,
Figure FDA0002612203890000027
and
Figure FDA0002612203890000028
the active power upper and lower limits of the thermal power generating unit g are respectively set;
Figure FDA0002612203890000029
the active power of the thermal power generating unit g is the typical day h period of the mth year, mth month in the planning cycle;
Figure FDA00026122038900000210
representing the increased deep peak regulation power of the thermal power generating unit through flexible modification until the y year; delta PgThe method comprises the steps that the depth peak regulation power which can be increased by flexible modification of a thermal power generating unit g is shown;
Figure FDA00026122038900000211
the binary variable represents whether the thermal power generating unit g is subjected to flexible modification in the y year in the planning period,
Figure FDA00026122038900000212
the thermal power generating unit g is flexibly modified in the y year in the planning period,
Figure FDA00026122038900000213
representing that the thermal power generating unit g is not subjected to flexibility modification in the y year in a planning period; zgIs a binary variable which indicates whether the live generator g is flexibly modified in a planning period, Zg1 denotes the flexible modification of the generator set g in the planning cycle, ZgWhen the number is 0, the live-wire generator g is not flexibly transformed in a planning period;
2) conventional thermal power generating unit climbing restraint:
Figure FDA0002612203890000031
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgThe method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure FDA0002612203890000032
the active power of the thermal power generating unit g is the active power of the thermal power generating unit g in a typical day h-1 of mth month in the planned cycle;
3) node power balance constraint:
Figure FDA0002612203890000033
wherein, Pl y,m,hThe transmission active power of a power transmission line corridor l in a typical day h period of a month m in the y year in a planning cycle;
Figure FDA0002612203890000034
the active power predicted value of the wind power plant w in a typical day h period of the mth year m month in the planning cycle;
Figure FDA0002612203890000035
the active power predicted value of the load d in a typical day h period of the mth year m month in the planning cycle; siAnd EiThe total number of the power transmission lines with the node i as a head node and the tail end node respectively; gi、WiAnd DiRespectively representing thermal power generating units and wind on node iThe total number of electric fields and loads; nb is the number of the nodes of the power grid;
4) and (5) abandoned wind power constraint:
Figure FDA0002612203890000036
Figure FDA0002612203890000037
wherein, γyRepresenting a set wind power curtailment rate threshold allowed in the y year;
5) transmission capacity constraint of the transmission line:
Figure FDA0002612203890000038
Figure FDA0002612203890000039
wherein, BlThe susceptance of a single-circuit power transmission line on a power transmission line corridor l;
Figure FDA00026122038900000310
a node i voltage phase angle is a typical day h period of a mth month in the y year in a planning cycle;
Figure FDA0002612203890000041
the number of the existing transmission lines on the transmission line corridor l is shown;
Figure FDA0002612203890000042
representing the number of newly added transmission lines in the t year on the transmission line corridor l; pl y,m,hThe transmission active power of a power transmission line corridor l in a typical day h period of a month m in the y year in a planning cycle; pl maxThe transmission capacity of a single-circuit transmission line on a transmission line corridor l;
6) the maximum extensible line number constraint of the power transmission line corridor is as follows:
Figure FDA0002612203890000043
wherein Z isl maxThe maximum number of the extensible power transmission lines is 1;
7) node voltage phase angle constraint:
Figure FDA0002612203890000044
wherein,
Figure FDA00026122038900000410
a node i voltage phase angle is a typical day h period of a mth month in the y year in a planning cycle;
8) n-1, the upper and lower limits of active power of the thermal power generating unit are constrained under the condition of an expected accident:
Figure FDA0002612203890000045
wherein,
Figure FDA0002612203890000046
active power of the thermal power generating unit g under an expected accident k in a typical day h period of m month in the y year in a planning cycle, wherein nk is the total number of the expected accidents of the power transmission line N-1;
9) n-1 conventional thermal power generating unit climbing restraint under the condition of anticipated accidents:
Figure FDA0002612203890000047
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgThe method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure FDA0002612203890000048
forecast accident k for typical day h period of mth month m of the y year in planning cycleThe active power of the lower thermal power generating unit g;
Figure FDA0002612203890000049
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h-1 of mth month in the y year in a planning cycle;
10) n-1 node power balance constraint under the condition of an expected accident:
Figure FDA0002612203890000051
wherein, Pl y,m,h,kThe transmission active power of the transmission line corridor l under the expected accident k is predicted for the typical day h period of the mth month in the y year in the planning cycle;
Figure FDA0002612203890000052
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h period of mth month in the y year in a planning cycle;
Figure FDA0002612203890000053
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth year, mth month in the planning cycle;
11) n-1 abandon wind power restraint under the condition of anticipated accidents:
Figure FDA0002612203890000054
wherein,
Figure FDA0002612203890000055
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth year, mth month in the planning cycle;
12) n-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure FDA0002612203890000056
Figure FDA0002612203890000057
wherein nl is the total number of the transmission line corridors; b islThe susceptance of a single-circuit power transmission line on a power transmission line corridor l; ny is the total year of the planning cycle;
Figure FDA0002612203890000058
the number of the existing transmission lines on the transmission line corridor l is shown;
Figure FDA0002612203890000059
representing the number of newly added transmission lines in the t year on the transmission line corridor l;
Figure FDA00026122038900000510
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure FDA00026122038900000511
indicating a wire outage in the k power line corridor under an expected accident,
Figure FDA00026122038900000512
representing that no line is stopped on the transmission line corridor under the expected accident k; pl y,m,h,kThe transmission active power of the transmission line corridor l under the expected accident k is predicted for the typical day h period of the mth month in the y year in the planning cycle; pl maxThe transmission capacity of a single-circuit transmission line on a transmission line corridor l;
13) n-1 node voltage phase angle constraint under the expected accident condition:
Figure FDA0002612203890000061
wherein,
Figure FDA0002612203890000066
for planning typical days of mth month in the y year in the periodThe h period envisions the node i voltage phase angle at incident k.
5. The method for multi-stage power transmission network planning with consideration of thermal power generating unit flexibility according to claim 4, wherein a non-linear constraint expression in the transmission capacity constraint of the transmission line is converted into an equivalent linear expression:
Figure FDA0002612203890000062
Figure FDA0002612203890000063
wherein M is a constant.
6. The method for multi-stage power transmission network planning considering thermal power generating unit flexibility according to claim 4, wherein the nonlinear constraint expression in the transmission capacity constraint of the transmission line under the N-1 expected accident condition is converted into an equivalent linear expression, namely:
Figure FDA0002612203890000064
Figure FDA0002612203890000065
7. a multi-stage grid planning system that takes into account thermal power plant flexibility, the system comprising:
the planning model construction module is used for constructing a multi-stage power transmission network planning optimization model considering flexibility transformation of the thermal power generating unit, the optimization model takes the minimum sum of the investment cost and the present value of the operation cost in the whole planning period as a target and comprises a plurality of constraint conditions;
and the model solving module is used for converting the mixed integer nonlinear constraint condition in the optimization model into a mixed integer linear constraint condition, and solving the converted mixed integer linear programming model by using a mixed integer linear programming algorithm to obtain a final thermal power unit flexibility modification and power transmission network multi-stage combined programming scheme.
8. The multi-stage power transmission network planning system considering the flexibility of the thermal power generating unit according to claim 7, further comprising a model parameter setting module, configured to input parameters of a conventional thermal power generating unit and parameters of a power transmission element of a current power grid before a model is constructed, input maximum load, maximum wind power, and 24-hour load power and wind power change per unit value of a typical day per month of each year in a planning period, set a range of the thermal power generating unit that can participate in the flexibility modification according to an operation condition of the thermal power generating unit and set a modification cost, set a range of a candidate new power transmission line according to a corridor condition of the power transmission line and set a construction cost, and set a wind power curtailment cost coefficient and a power curtailment rate allowed by the system.
9. The multi-stage power transmission network planning system considering thermal power generating unit flexibility according to claim 7, wherein an objective function expression of the optimization model is as follows:
Figure FDA0002612203890000071
wherein nw is the total number of the wind power plants; ny is the total year of the planning cycle; ng is the number of conventional thermal power generating units; nl is the number of transmission line corridors planned and constructed; alpha is alphawRepresenting wind curtailment penalty cost;
Figure FDA0002612203890000072
the electric power curtailment of the wind power plant w is the typical day h period of the mth year, mth month in the planning cycle; beta is agRepresenting the flexibility modification cost of the thermal power generating unit g;
Figure FDA0002612203890000073
for binary variables, representingWhether the thermal power generating unit g is flexibly modified in the y year in the planning period,
Figure FDA0002612203890000074
the thermal power generating unit g is flexibly modified in the y year in the planning period,
Figure FDA0002612203890000075
showing that the thermal power generating unit g does not perform flexible modification in the y year in the planning period; clInvestment cost for construction of a single-circuit power transmission line on a power transmission corridor l;
Figure FDA0002612203890000076
representing the number of newly added power transmission lines in the y year on the power transmission line corridor l; ρ is the discount rate.
10. The multi-stage power transmission network planning system considering thermal power generating unit flexibility according to claim 7, wherein the constraint condition comprises:
1) the upper and lower limits of active power of a conventional thermal power generating unit are constrained:
Figure FDA0002612203890000081
Figure FDA0002612203890000082
Figure FDA0002612203890000083
wherein,
Figure FDA0002612203890000084
and
Figure FDA0002612203890000085
the active power upper and lower limits of the thermal power generating unit g are respectively set;
Figure FDA0002612203890000086
the active power of the thermal power generating unit g is the typical day h period of the mth year, mth month in the planning cycle;
Figure FDA0002612203890000087
representing the increased deep peak regulation power of the thermal power generating unit through flexible modification until the y year; delta PgThe method comprises the steps that the depth peak regulation power which can be increased by flexible modification of a thermal power generating unit g is shown;
Figure FDA0002612203890000088
the binary variable represents whether the thermal power generating unit g is subjected to flexible modification in the y year in the planning period,
Figure FDA0002612203890000089
the thermal power generating unit g is flexibly modified in the y year in the planning period,
Figure FDA00026122038900000810
representing that the thermal power generating unit g is not subjected to flexibility modification in the y year in a planning period; zgIs a binary variable which indicates whether the live generator g is flexibly modified in a planning period, Zg1 denotes the flexible modification of the generator set g in the planning cycle, ZgWhen the number is 0, the live-wire generator g is not flexibly transformed in a planning period;
2) conventional thermal power generating unit climbing restraint:
Figure FDA00026122038900000811
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgThe method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure FDA00026122038900000812
a typical day h-1 period of m months in the y year in a planning cycleActive power of the thermal power generating unit g;
3) node power balance constraint:
Figure FDA0002612203890000091
wherein, Pl y,m,hThe transmission active power of a power transmission line corridor l in a typical day h period of a month m in the y year in a planning cycle;
Figure FDA0002612203890000092
the active power predicted value of the wind power plant w in a typical day h period of the mth year m month in the planning cycle;
Figure FDA0002612203890000093
the active power predicted value of the load d in a typical day h period of the mth year m month in the planning cycle; siAnd EiThe total number of the power transmission lines with the node i as a head node and the tail end node respectively; gi、WiAnd DiRespectively representing the total number of thermal power generating units, wind power plants and loads on the node i; nb is the number of the nodes of the power grid;
4) and (5) abandoned wind power constraint:
Figure FDA0002612203890000094
Figure FDA0002612203890000095
wherein, γyRepresenting a set wind power curtailment rate threshold allowed in the y year;
5) transmission capacity constraint of the transmission line:
Figure FDA0002612203890000096
Figure FDA0002612203890000097
wherein, BlThe susceptance of a single-circuit power transmission line on a power transmission line corridor l;
Figure FDA0002612203890000098
a node i voltage phase angle is a typical day h period of a mth month in the y year in a planning cycle;
Figure FDA0002612203890000099
the number of the existing transmission lines on the transmission line corridor l is shown;
Figure FDA00026122038900000910
representing the number of newly added transmission lines in the t year on the transmission line corridor l; pl y,m,hThe transmission active power of a power transmission line corridor l in a typical day h period of a month m in the y year in a planning cycle; pl maxThe transmission capacity of a single-circuit transmission line on a transmission line corridor l;
6) the maximum extensible line number constraint of the power transmission line corridor is as follows:
Figure FDA0002612203890000101
wherein Z isl maxThe maximum number of the extensible power transmission lines is 1;
7) node voltage phase angle constraint:
Figure FDA0002612203890000102
wherein,
Figure FDA0002612203890000103
a node i voltage phase angle is a typical day h period of a mth month in the y year in a planning cycle;
8) n-1, the upper and lower limits of active power of the thermal power generating unit are constrained under the condition of an expected accident:
Figure FDA0002612203890000104
wherein,
Figure FDA0002612203890000105
active power of the thermal power generating unit g under an expected accident k in a typical day h period of m month in the y year in a planning cycle, wherein nk is the total number of the expected accidents of the power transmission line N-1;
9) n-1 conventional thermal power generating unit climbing restraint under the condition of anticipated accidents:
Figure FDA0002612203890000106
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgThe method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure FDA0002612203890000107
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h period of mth month in the y year in a planning cycle;
Figure FDA0002612203890000108
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h-1 of mth month in the y year in a planning cycle;
10) n-1 node power balance constraint under the condition of an expected accident:
Figure FDA0002612203890000109
wherein, Pl y,m,h,kThe transmission active power of the transmission line corridor l under the expected accident k is predicted for the typical day h period of the mth month in the y year in the planning cycle;
Figure FDA00026122038900001010
active power of the thermal power generating unit g under an expected accident k is predicted for a typical day h period of mth month in the y year in a planning cycle;
Figure FDA0002612203890000111
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth year, mth month in the planning cycle;
11) n-1 abandon wind power restraint under the condition of anticipated accidents:
Figure FDA0002612203890000112
wherein,
Figure FDA0002612203890000113
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth year, mth month in the planning cycle;
12) n-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure FDA0002612203890000114
Figure FDA0002612203890000115
wherein,nl is the total number of the power transmission line corridors; b islThe susceptance of a single-circuit power transmission line on a power transmission line corridor l; ny is the total year of the planning cycle;
Figure FDA0002612203890000116
the number of the existing transmission lines on the transmission line corridor l is shown;
Figure FDA0002612203890000117
representing the number of newly added transmission lines in the t year on the transmission line corridor l;
Figure FDA0002612203890000118
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure FDA0002612203890000119
indicating a wire outage in the k power line corridor under an expected accident,
Figure FDA00026122038900001110
representing that no line is stopped on the transmission line corridor under the expected accident k; pl y,m,h,kThe transmission active power of the transmission line corridor l under the expected accident k is predicted for the typical day h period of the mth month in the y year in the planning cycle; pl maxThe transmission capacity of a single-circuit transmission line on a transmission line corridor l;
13) n-1 node voltage phase angle constraint under the expected accident condition:
Figure FDA00026122038900001111
wherein,
Figure FDA00026122038900001112
and forecasting the voltage phase angle of the node i under the accident k for the typical day h period of the mth year m month in the planning cycle.
CN202010758007.5A 2020-07-31 2020-07-31 Multi-stage power transmission network planning method and system considering flexibility of thermal power generating unit Active CN111799842B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010758007.5A CN111799842B (en) 2020-07-31 2020-07-31 Multi-stage power transmission network planning method and system considering flexibility of thermal power generating unit

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010758007.5A CN111799842B (en) 2020-07-31 2020-07-31 Multi-stage power transmission network planning method and system considering flexibility of thermal power generating unit

Publications (2)

Publication Number Publication Date
CN111799842A true CN111799842A (en) 2020-10-20
CN111799842B CN111799842B (en) 2023-11-10

Family

ID=72827921

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010758007.5A Active CN111799842B (en) 2020-07-31 2020-07-31 Multi-stage power transmission network planning method and system considering flexibility of thermal power generating unit

Country Status (1)

Country Link
CN (1) CN111799842B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116703084A (en) * 2023-06-05 2023-09-05 东北林业大学 Multi-technology flexibility transformation planning method for coal-fired unit based on high wind power permeability

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105528466A (en) * 2014-09-28 2016-04-27 国家电网公司 Wind power optimal planning modeling method considering adaptability and economy of power system
CN106953354A (en) * 2017-03-10 2017-07-14 国网山东省电力公司经济技术研究院 Consider the method for Unit Commitment containing wind-powered electricity generation of voltage support
CN108847667A (en) * 2018-08-03 2018-11-20 国网山东省电力公司经济技术研究院 A kind of method for expansion planning of power transmission network considering electric network composition optimization
CN109102115A (en) * 2018-08-03 2018-12-28 国网山东省电力公司经济技术研究院 A kind of reference power grid chance constrained programming method adapting to wind-powered electricity generation large-scale grid connection
CN109146706A (en) * 2018-08-14 2019-01-04 国网四川省电力公司经济技术研究院 A kind of Transmission Expansion Planning in Electric method considering the flexibility equilibrium of supply and demand
CN109165773A (en) * 2018-08-03 2019-01-08 国网山东省电力公司经济技术研究院 A kind of Transmission Expansion Planning in Electric evolutionary structural optimization
CN110429591A (en) * 2019-08-02 2019-11-08 西安交通大学 A kind of power transmission network utilization rate appraisal procedure based on electric system timing coupling
CN110808597A (en) * 2019-11-06 2020-02-18 山东电力工程咨询院有限公司 Distributed power supply planning method considering three-phase imbalance in active power distribution network
CN110852565A (en) * 2019-10-10 2020-02-28 国家电网有限公司 Power transmission network frame planning method considering different functional attributes
CN110880791A (en) * 2019-12-24 2020-03-13 国网节能服务有限公司 Coordination optimization method for hybrid alternating current-direct current power distribution network

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105528466A (en) * 2014-09-28 2016-04-27 国家电网公司 Wind power optimal planning modeling method considering adaptability and economy of power system
CN106953354A (en) * 2017-03-10 2017-07-14 国网山东省电力公司经济技术研究院 Consider the method for Unit Commitment containing wind-powered electricity generation of voltage support
CN108847667A (en) * 2018-08-03 2018-11-20 国网山东省电力公司经济技术研究院 A kind of method for expansion planning of power transmission network considering electric network composition optimization
CN109102115A (en) * 2018-08-03 2018-12-28 国网山东省电力公司经济技术研究院 A kind of reference power grid chance constrained programming method adapting to wind-powered electricity generation large-scale grid connection
CN109165773A (en) * 2018-08-03 2019-01-08 国网山东省电力公司经济技术研究院 A kind of Transmission Expansion Planning in Electric evolutionary structural optimization
CN109146706A (en) * 2018-08-14 2019-01-04 国网四川省电力公司经济技术研究院 A kind of Transmission Expansion Planning in Electric method considering the flexibility equilibrium of supply and demand
CN110429591A (en) * 2019-08-02 2019-11-08 西安交通大学 A kind of power transmission network utilization rate appraisal procedure based on electric system timing coupling
CN110852565A (en) * 2019-10-10 2020-02-28 国家电网有限公司 Power transmission network frame planning method considering different functional attributes
CN110808597A (en) * 2019-11-06 2020-02-18 山东电力工程咨询院有限公司 Distributed power supply planning method considering three-phase imbalance in active power distribution network
CN110880791A (en) * 2019-12-24 2020-03-13 国网节能服务有限公司 Coordination optimization method for hybrid alternating current-direct current power distribution network

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
孙东磊;李雪亮;韩学山;李山;杨思;杨金洪;: "考虑运行方式优化和拓扑校正控制的参考电网优化方法", 电力***保护与控制, no. 13, pages 97 - 102 *
杨波;李政;: "火电机组热力***主导因素变工况建模方法研究", 中国电机工程学报, no. 24, pages 96 - 102 *
王静;廖胜利;程春田;蔡华祥;: "中期火电开机优化的混合搜索算法", 中国电机工程学报, no. 01, pages 94 - 100 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116703084A (en) * 2023-06-05 2023-09-05 东北林业大学 Multi-technology flexibility transformation planning method for coal-fired unit based on high wind power permeability
CN116703084B (en) * 2023-06-05 2024-03-22 东北林业大学 Multi-technology flexibility transformation planning method for coal-fired unit based on high wind power permeability

Also Published As

Publication number Publication date
CN111799842B (en) 2023-11-10

Similar Documents

Publication Publication Date Title
CN110739725B (en) Optimal scheduling method for power distribution network
CN106327091B (en) Multi-region asynchronous coordination dynamic economic dispatching method based on robust tie line plan
CN107968439B (en) Active power distribution network joint optimization algorithm based on mixed integer linear programming
CN109508499A (en) Multi-period more optimal on-positions of scene distribution formula power supply and capacity research method
CN102496968A (en) Generation plan optimizing method in intermittent energy and conventional energy coordinated dispatching mode
CN112072710B (en) Source network load integrated economic dispatching method and system considering demand response
CN111144655A (en) Combined optimization method for site selection, volume fixing and power distribution network frame of distributed power supply
CN101989743B (en) Energy-saving power generation dispatching optimization method based on direct current power flow
CN115422758A (en) Active power distribution network intelligent soft switch planning method and device based on incremental linearization
CN110867907B (en) Power system scheduling method based on multi-type power generation resource homogenization
CN111799793B (en) Source-grid-load cooperative power transmission network planning method and system
CN111799841B (en) Combined decision method and system for thermal power generating unit flexibility transformation and power transmission planning
CN116402210A (en) Multi-objective optimization method, system, equipment and medium for comprehensive energy system
CN116526469A (en) Long-term random dynamic scheduling method for water-wind-solar complementary system
Hu et al. Adaptive time division power dispatch based on numerical characteristics of net loads
CN111799842B (en) Multi-stage power transmission network planning method and system considering flexibility of thermal power generating unit
CN107392350B (en) Comprehensive optimization method for power distribution network extension planning containing distributed energy and charging stations
CN113363976A (en) Scene graph-based mid-term optimized scheduling method for wind, light and water complementary power generation system
CN112039126A (en) Multi-time scale coordinated scheduling method and system for power distribution network containing distributed power supply
CN116540545A (en) Photovoltaic power generation hydrogen production cluster random optimization scheduling method based on ember process
CN108718084B (en) Power supply and power grid coordination planning method suitable for electric power market reformation
CN109213104B (en) Scheduling method and scheduling system of energy storage system based on heuristic dynamic programming
CN115425650A (en) Power supply station microgrid configuration method, device, equipment and medium
CN114884101A (en) Pumped storage dispatching method based on self-adaptive model control prediction
CN114358430A (en) Novel rural power distribution network planning strategy based on robust optimization

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant