CN111799793A - Source-grid-load cooperative power transmission network planning method and system - Google Patents

Source-grid-load cooperative power transmission network planning method and system Download PDF

Info

Publication number
CN111799793A
CN111799793A CN202010757955.7A CN202010757955A CN111799793A CN 111799793 A CN111799793 A CN 111799793A CN 202010757955 A CN202010757955 A CN 202010757955A CN 111799793 A CN111799793 A CN 111799793A
Authority
CN
China
Prior art keywords
power
period
load
typical day
target year
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
CN202010757955.7A
Other languages
Chinese (zh)
Other versions
CN111799793B (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 CN202010757955.7A priority Critical patent/CN111799793B/en
Publication of CN111799793A publication Critical patent/CN111799793A/en
Application granted granted Critical
Publication of CN111799793B publication Critical patent/CN111799793B/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/003—Load forecast, e.g. methods or systems for forecasting future load demand
    • 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • 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
    • 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)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Power Engineering (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • Quality & Reliability (AREA)
  • Development Economics (AREA)
  • Water Supply & Treatment (AREA)
  • Educational Administration (AREA)
  • Primary Health Care (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a power transmission network planning method and system with source-grid-load coordination, which are based on the time sequence characteristics of flexible modification of a source-side thermal power generating unit, a demand-side response mechanism, typical daily load, wind power and the like, by taking economy as a target and combining reasonable wind curtailment allowed by system operation practice, a power transmission network planning optimization model with source network load coordination is constructed, converting the mixed integer nonlinear constraint condition in the optimization model into a mixed integer linear constraint condition, solving the converted mixed integer linear programming model by adopting a mixed integer linear programming method, therefore, the final thermal power generating unit flexibility transformation, load demand response mechanism and power transmission network planning scheme is obtained, the combined decision of thermal power generating unit flexibility transformation, load demand response mechanism and power transmission network planning is carried out, and the economy of power system planning investment is improved while wind power consumption is promoted.

Description

Source-grid-load cooperative power transmission network planning method and system
Technical Field
The invention relates to the technical field of power grid planning, in particular to a power transmission network planning method and system with source-grid-load coordination.
Background
In recent years, with the increasing energy crisis and environmental pollution problems, the traditional energy consumption mode based on fossil energy utilization is difficult to continue, and under the background, clean energy represented by wind power is rapidly developed in a large scale. However, under the influence of natural factors, wind power generation has random fluctuation, higher requirements are provided for planning and operating of a power system, and in order to deal with the change of wind power output and reduce the condition of wind abandonment, the flexibility of the source network load side of the power system needs to be improved.
The mode of improving the flexibility of the source side comprises three main modes of building a pumped storage power station, building a flexible gas power plant and performing flexibility transformation on a traditional thermal power generating unit. The former two modes need to increase occupied land resources, and have high investment cost and long construction period. The method has the advantages that the thermal power generating unit is flexibly transformed without increasing land resources, the investment cost is low, the traditional coal-fired thermal power generating unit is only technically transformed, the period is short, and the method is a feasible mode for improving the flexibility of the source side. In China, the total installed amount of the thermal power generating unit is rich but the flexibility is insufficient, and by the end of 2019, the flexibility of the traditional thermal power generating unit is improved by 1040GW (coal-fired thermal power generating unit) in China, so that the flexibility space of the source side is greatly improved. However, when the thermal power generating units are modified, how to make a time sequence of a flexibility modification plan of the thermal power generating units is worthy of deep research.
The mode for improving the flexibility of the load side comprises the steps of building an energy storage power station and managing the real-time power demand side. The former needs to increase occupied land resources and has high investment cost. The latter only needs to guide the user to improve the electricity utilization behavior through an incentive mechanism, and the cost is relatively low. The demand side management can optimize a power grid load demand curve, improve the load rate of power grid equipment, reduce the capacity demand of power transmission and transformation equipment and save the construction investment of a power grid.
The method for improving the flexibility of the grid side comprises the steps of applying a novel flexible power transmission technology and optimizing a power grid planning method. The novel flexible power transmission technology is applied, so that the investment cost of power grid construction is greatly increased, and the reliability of the novel technology is relatively low. The optimization power grid planning method is to improve the traditional power grid planning method so that the optimization power grid planning method can adapt to diversified source load situations. Traditional power system planning meets the ever-increasing peak load demands by simply expanding the power and transmission capacity, making the utilization of power transmission and transformation equipment low. In the power grid planning method under the new situation, the actual operation condition of the source grid load needs to be simulated in the planning process, so that the power grid planning scheme can be matched with the actual operation condition of the power system.
Disclosure of Invention
The invention aims to provide a source-grid-load cooperative power transmission network planning method and system, and aims to solve the problem that power grid planning in the prior art cannot take account of low utilization rate caused by source-grid-load side, so that wind power consumption is promoted, and the economy of planning investment of a power system is improved.
In order to achieve the technical purpose, the invention provides a source-grid-load coordinated power transmission network planning method, which comprises the following operations:
constructing a power transmission network planning optimization model with source network load coordination, wherein the optimization model takes the minimum sum of investment cost and operation cost in the whole planning period as a target and comprises constraint conditions;
and converting the mixed integer nonlinear constraint condition in the optimization model into a mixed integer linear constraint condition, converting the optimization model into a mixed integer linear programming model, and directly solving by using a mixed integer linear programming algorithm to obtain a final thermal power unit flexibility modification, load demand response mechanism and power transmission network 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 capable of participating in flexibility modification and a modification cost according to the operation condition of the thermal power unit, setting a response quantity range and a response price of demand response load participating in demand response, setting a candidate range of a 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 system-allowed curtailment rate.
Preferably, the objective function expression of the optimization model is:
Figure BDA0002612204280000031
in the formula, nw is the total number of the wind power plants; ng is the number of conventional thermal power generating units; nl is the number of transmission line corridors planned and constructed; nc and nr represent the number of energizable loads and interruptible loads, respectively; alpha is alphawRepresenting wind curtailment penalty cost;
Figure BDA0002612204280000032
the electric power curtailment of the wind power plant w is planned for the typical day h period of the mth month of the target year; beta is agRepresenting the flexibility modification cost of the thermal power generating unit g; u. ofgIs a binary variable and indicates whether the thermal power generating unit g is subjected to flexible transformation in a planned target year or not, ug1 represents that the thermal power generating unit g is subjected to flexibility modification in a planning target year, and u representsgWhen the thermal power generating unit g is not transformed flexibly in the planning target year, 0 is shown; clInvestment cost for construction of a single-circuit power transmission line on a power transmission corridor l; zlRepresenting the number of newly added transmission lines on the transmission line corridor l;
Figure BDA0002612204280000033
active power which can stimulate the load c to increase for planning the typical day h period of the mth month of the target year;
Figure BDA0002612204280000034
marginal price for the energizable load c;
Figure BDA0002612204280000035
for planning the typical day of the mth month of the target yearThe active power of the load r interrupted in the h time period can be interrupted; etarThe marginal price of the interruptible load r for a time period t.
Preferably, the constraint condition includes:
1) and (3) constraining the upper limit and the lower limit of the active power of the thermal power generating unit:
Figure BDA0002612204280000036
wherein the content of the first and second substances,
Figure BDA0002612204280000037
and
Figure BDA0002612204280000038
the active power upper and lower limits of the thermal power generating unit g are respectively set;
Figure BDA0002612204280000039
the active power of the thermal power generating unit g in a typical day h period of the mth month of a planned target year is planned; 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;
2) conventional thermal power generating unit climbing restraint:
Figure BDA0002612204280000041
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;
3) node power balance constraint:
Figure BDA0002612204280000042
wherein the content of the first and second substances,
Figure BDA0002612204280000043
the transmission active power of the transmission line corridor l is planned in a typical day h period of the mth month of the target year;
Figure BDA0002612204280000044
the active power predicted value of the wind power plant w in a typical day h period of the mth month of the planning target year;
Figure BDA0002612204280000045
an active power predicted value of a load d in a typical day h period of the mth month of a planning target year;
Figure BDA0002612204280000046
an active power predicted value of an energizable load c for planning a typical day h period of the mth month of a target year;
Figure BDA0002612204280000047
an active power predicted value of an interruptible load r in a typical day h period of the mth month of a planned target year; 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; ciAnd RiRespectively representing the number of excitable loads and interruptible loads on a node i; nb is the number of the nodes of the power grid;
4) and (5) abandoned wind power constraint:
Figure BDA0002612204280000048
wherein gamma represents a set wind power curtailment rate threshold value;
5) demand response load response range constraints
Figure BDA0002612204280000049
Figure BDA00026122042800000410
Wherein the content of the first and second substances,
Figure BDA0002612204280000051
and
Figure BDA0002612204280000052
the upper and lower limits of active power of the energizable load c are planned in a typical day h period of the mth month of the target year;
Figure BDA0002612204280000053
and
Figure BDA0002612204280000054
the upper and lower limits of the active power of the interruptible load r in a typical day h period of the mth month of the planned target year are planned;
6) transmission capacity constraint of the transmission line:
Figure BDA0002612204280000055
Figure BDA0002612204280000056
wherein, BlThe susceptance of a single-circuit power transmission line on a power transmission line corridor l;
Figure BDA0002612204280000057
a node i voltage phase angle is planned for a typical day h period of the mth month of a target year;
Figure BDA0002612204280000058
the number of the existing transmission lines on the transmission line corridor l is shown; zl represents the number of newly added transmission lines on the transmission line corridor l;
Figure BDA0002612204280000059
the transmission active power of the transmission line corridor l is planned in a typical day h period of the mth month of the target year;
Figure BDA00026122042800000510
the transmission capacity of a single-circuit transmission line on a transmission line corridor l;
7) the maximum extensible line number constraint of the power transmission line corridor is as follows:
Zl≤Zl max,l＝1,L,nl
wherein Z isl maxThe maximum number of the extensible power transmission lines is 1;
8) node voltage phase angle constraint:
Figure BDA00026122042800000511
wherein the content of the first and second substances,
Figure BDA00026122042800000512
a node i voltage phase angle is planned for a typical day h period of the mth month of a target year;
9) 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 BDA00026122042800000513
wherein the content of the first and second substances,
Figure BDA00026122042800000514
in order to plan the active power of the thermal power generating unit g under the expected accident k in the typical day h of the mth month of the target year, nk is the total number of the expected accidents of the power transmission line N-1;
10) n-1 conventional thermal power generating unit climbing restraint under the condition of anticipated accidents:
Figure BDA0002612204280000061
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 BDA0002612204280000062
forecasting the active power of the thermal power generating unit g under an accident k for planning the typical day h period of the mth month of the target year;
Figure BDA0002612204280000063
for planning the mth month dictionary of the target yearActive power of the thermal power generating unit g under an expected accident k in a model day h-1 period;
11) n-1 node power balance constraint under the condition of an expected accident:
Figure BDA0002612204280000064
wherein the content of the first and second substances,
Figure BDA0002612204280000065
the transmission active power of a transmission line corridor l under an expected accident k is planned for the mth month typical day h period of a target year;
Figure BDA0002612204280000066
forecasting the active power of the thermal power generating unit g under an accident k for planning the typical day h period of the mth month of the target year;
Figure BDA0002612204280000067
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth month of the planned target year;
Figure BDA0002612204280000068
active power increased by an energizable load c under an expected accident k is planned for a typical day h period of the mth month of the planned target year;
Figure BDA0002612204280000069
active power of an interruptible load r interrupted under an expected accident k for planning a typical day h period of an mth month of a target year;
12) n-1 abandon wind power restraint under the condition of anticipated accidents:
Figure BDA00026122042800000610
wherein the content of the first and second substances,
Figure BDA00026122042800000611
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth month of the planned target year;
13) demand response load response range constraints under N-1 anticipated accident conditions
Figure BDA00026122042800000612
Figure BDA00026122042800000613
14) N-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure BDA0002612204280000071
Figure BDA0002612204280000072
wherein the content of the first and second substances,
Figure BDA0002612204280000073
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure BDA0002612204280000074
indicating a wire outage in the k power line corridor under an expected accident,
Figure BDA0002612204280000075
representing that no line is stopped on the transmission line corridor under the expected accident k;
Figure BDA0002612204280000076
the transmission active power of a transmission line corridor l under an expected accident k is planned for the mth month typical day h period of a target year;
15) n-1 node voltage phase angle constraint under the expected accident condition:
Figure BDA0002612204280000077
wherein the content of the first and second substances,
Figure BDA0002612204280000078
a voltage phase angle of a node i under an accident k is predicted for a typical day h period of an mth month of a planned target year.
Preferably, the non-linear constraint expression in the transmission capacity constraint of the power transmission line is converted into an equivalent linear expression:
Figure BDA0002612204280000079
Figure BDA00026122042800000710
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 BDA00026122042800000711
Figure BDA0002612204280000081
the invention also provides a source network and load coordinated power transmission network planning system, which comprises:
the planning model construction module is used for constructing a power transmission network planning optimization model with source network load coordination, and the optimization model takes the minimum sum of investment cost and operation cost in the whole planning period as a target and contains 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, converting the optimization model into a mixed integer linear programming model, and directly solving by using a mixed integer linear programming algorithm to obtain a final thermal power unit flexibility modification, load demand response mechanism and power transmission network 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 modification cost according to the operation condition of the thermal power unit, setting a response quantity range and response price of demand response load participating in demand response, setting a candidate range of a newly-built power transmission line and 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 BDA0002612204280000082
in the formula, nw is the total number of the wind power plants; ng is the number of conventional thermal power generating units; nl is the number of transmission line corridors planned and constructed; nc and nr represent the number of energizable loads and interruptible loads, respectively; alpha is alphawRepresenting wind curtailment penalty cost;
Figure BDA0002612204280000083
the electric power curtailment of the wind power plant w is planned for the typical day h period of the mth month of the target year; beta is agRepresenting the flexibility modification cost of the thermal power generating unit g; u. ofgIs a binary variable and indicates whether the thermal power generating unit g is subjected to flexible transformation in a planned target year or not, ug1 represents that the thermal power generating unit g is subjected to flexibility modification in a planning target year, and u representsgWhen the thermal power generating unit g is not transformed flexibly in the planning target year, 0 is shown; clInvestment cost for construction of a single-circuit power transmission line on a power transmission corridor l; zlRepresenting the number of newly added transmission lines on the transmission line corridor l;
Figure BDA0002612204280000091
active power which can stimulate the load c to increase for planning the typical day h period of the mth month of the target year;
Figure BDA0002612204280000092
marginal price for the energizable load c;
Figure BDA0002612204280000093
the active power of the interruptible load r interrupted in a typical day h period of the mth month of the planning target year; etarThe marginal price of the interruptible load r for a time period t.
Preferably, the constraint condition includes:
1) and (3) constraining the upper limit and the lower limit of the active power of the thermal power generating unit:
Figure BDA0002612204280000094
wherein the content of the first and second substances,
Figure BDA0002612204280000095
and
Figure BDA0002612204280000096
the active power upper and lower limits of the thermal power generating unit g are respectively set;
Figure BDA0002612204280000097
the active power of the thermal power generating unit g in a typical day h period of the mth month of a planned target year is planned; 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;
2) conventional thermal power generating unit climbing restraint:
Figure BDA0002612204280000098
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;
3) node power balance constraint:
Figure BDA0002612204280000099
wherein the content of the first and second substances,
Figure BDA00026122042800000910
the transmission active power of the transmission line corridor l is planned in a typical day h period of the mth month of the target year;
Figure BDA00026122042800000911
the active power predicted value of the wind power plant w in a typical day h period of the mth month of the planning target year;
Figure BDA00026122042800000912
an active power predicted value of a load d in a typical day h period of the mth month of a planning target year;
Figure BDA00026122042800000913
an active power predicted value of an energizable load c for planning a typical day h period of the mth month of a target year;
Figure BDA0002612204280000101
an active power predicted value of an interruptible load r in a typical day h period of the mth month of a planned target year; 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; ciAnd RiRespectively representing the number of excitable loads and interruptible loads on a node i; nb is the number of the nodes of the power grid;
4) and (5) abandoned wind power constraint:
Figure BDA0002612204280000102
wherein gamma represents a set wind power curtailment rate threshold value;
5) demand response load response range constraints
Figure BDA0002612204280000103
Figure BDA0002612204280000104
Wherein the content of the first and second substances,
Figure BDA0002612204280000105
and
Figure BDA0002612204280000106
the upper and lower limits of active power of the energizable load c are planned in a typical day h period of the mth month of the target year;
Figure BDA0002612204280000107
and
Figure BDA0002612204280000108
the upper and lower limits of the active power of the interruptible load r in a typical day h period of the mth month of the planned target year are planned;
6) transmission capacity constraint of the transmission line:
Figure BDA0002612204280000109
Figure BDA00026122042800001010
wherein, BlThe susceptance of a single-circuit power transmission line on a power transmission line corridor l;
Figure BDA00026122042800001011
a node i voltage phase angle is planned for a typical day h period of the mth month of a target year;
Figure BDA00026122042800001012
the number of the existing transmission lines on the transmission line corridor l is shown; zlRepresenting the number of newly added transmission lines on the transmission line corridor l;
Figure BDA00026122042800001013
the transmission active power of the transmission line corridor l is planned in a typical day h period of the mth month of the target year;
Figure BDA00026122042800001014
the transmission capacity of a single-circuit transmission line on a transmission line corridor l;
7) the maximum extensible line number constraint of the power transmission line corridor is as follows:
Zl≤Zl max,l＝1,L,nl
wherein Z isl maxThe maximum number of the extensible power transmission lines is 1;
8) node voltage phase angle constraint:
Figure BDA0002612204280000111
wherein the content of the first and second substances,
Figure BDA0002612204280000112
a node i voltage phase angle is planned for a typical day h period of the mth month of a target year;
9) 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 BDA0002612204280000113
wherein the content of the first and second substances,
Figure BDA0002612204280000114
in order to plan the active power of the thermal power generating unit g under the expected accident k in the typical day h of the mth month of the target year, nk is the total number of the expected accidents of the power transmission line N-1;
10) n-1 conventional thermal power generating unit climbing restraint under the condition of anticipated accidents:
Figure BDA0002612204280000115
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 BDA0002612204280000116
for planningActive power of the thermal power generating unit g under an expected accident k in a typical day h period of the mth month of the target year;
Figure BDA0002612204280000117
active power of the thermal power generating unit g under an expected accident k is planned in a typical day h-1 period of the mth month of a planned target year;
11) n-1 node power balance constraint under the condition of an expected accident:
Figure BDA0002612204280000118
wherein the content of the first and second substances,
Figure BDA0002612204280000119
the transmission active power of a transmission line corridor l under an expected accident k is planned for the mth month typical day h period of a target year;
Figure BDA00026122042800001110
forecasting the active power of the thermal power generating unit g under an accident k for planning the typical day h period of the mth month of the target year;
Figure BDA0002612204280000121
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth month of the planned target year;
Figure BDA0002612204280000122
active power increased by an energizable load c under an expected accident k is planned for a typical day h period of the mth month of the planned target year;
Figure BDA0002612204280000123
active power of an interruptible load r interrupted under an expected accident k for planning a typical day h period of an mth month of a target year;
12) n-1 abandon wind power restraint under the condition of anticipated accidents:
Figure BDA0002612204280000124
wherein the content of the first and second substances,
Figure BDA0002612204280000125
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth month of the planned target year;
13) demand response load response range constraints under N-1 anticipated accident conditions
Figure BDA0002612204280000126
Figure BDA0002612204280000127
14) N-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure BDA0002612204280000128
Figure BDA0002612204280000129
wherein the content of the first and second substances,
Figure BDA00026122042800001210
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure BDA00026122042800001211
indicating a wire outage in the k power line corridor under an expected accident,
Figure BDA00026122042800001212
representing that no line is stopped on the transmission line corridor under the expected accident k;
Figure BDA00026122042800001213
the transmission active power of a transmission line corridor l under an expected accident k is planned for the mth month typical day h period of a target year;
15) n-1 node voltage phase angle constraint under the expected accident condition:
Figure BDA00026122042800001214
wherein the content of the first and second substances,
Figure BDA00026122042800001215
a voltage phase angle of a node i under an accident k is predicted for a typical day h period of an mth month of a planned target year.
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 sequential characteristics of flexibility modification of the thermal power unit at the source side, a demand side response mechanism, a typical daily load, wind power and the like, the economy is taken as a target, reasonable wind abandon is allowed by combining with the actual operation of the system, a power grid load cooperative power transmission network planning optimization model is constructed, mixed integer nonlinear constraint conditions in the optimization model are converted into mixed integer linear constraint conditions, a mixed integer linear programming method is adopted to solve the converted mixed integer linear programming model, and therefore the final thermal power unit flexibility modification, load demand response mechanism and power transmission network planning scheme is obtained, the combined decision of the thermal power unit flexibility modification, the load demand response mechanism and the power transmission network planning is carried out, and the economy of planning investment of a power system is improved while wind power consumption is promoted.
Drawings
Fig. 1 is a flowchart of a source-grid-load coordinated power transmission network planning method provided in an embodiment of the present invention;
fig. 2 is a block diagram of a source-grid-load coordinated power transmission network planning system provided in 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 a power transmission network planning method and system with source-grid-load coordination according to an embodiment of the present invention in detail with reference to the accompanying drawings.
As shown in fig. 1, the invention discloses a source-grid-load coordinated power transmission network planning method, which includes 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 range of the thermal power generating unit capable of participating in flexibility transformation and transformation cost according to the operation condition of the thermal power generating unit, setting a response quantity range and response price of demand response load participating in demand response, setting a range of a candidate newly-built power transmission line and construction cost according to the corridor condition of the power transmission line, and setting a wind power curtailment cost coefficient and a curtailment rate allowed by the system.
And S2, constructing a source network load coordinated power transmission network planning optimization model, wherein the optimization model aims at minimizing the sum of investment cost and operation cost in the whole planning period and comprises constraint conditions.
The objective function expression of the optimization model is as follows:
Figure BDA0002612204280000141
in the formula, nw is the total number of the wind power plants; ng is the number of conventional thermal power generating units; nl is the number of transmission line corridors planned and constructed; nc and nr represent the number of energizable loads and interruptible loads, respectively; alpha is alphawRepresenting wind curtailment penalty cost;
Figure BDA0002612204280000142
the electric power curtailment of the wind power plant w is planned for the typical day h period of the mth month of the target year; beta is agRepresenting the flexibility modification cost of the thermal power generating unit g; u. ofgIs a binary variable and indicates whether the thermal power generating unit g is subjected to flexible transformation in a planned target year or not, ug1 represents that the thermal power generating unit g is subjected to flexibility modification in a planning target year, and u representsgWhen the thermal power generating unit g is not transformed flexibly in the planning target year, 0 is shown; clInvestment cost for construction of a single-circuit power transmission line on a power transmission corridor l; zlRepresenting the number of newly added transmission lines on the transmission line corridor l;
Figure BDA0002612204280000143
active power which can stimulate the load c to increase for planning the typical day h period of the mth month of the target year;
Figure BDA0002612204280000151
marginal price for the energizable load c;
Figure BDA0002612204280000152
the active power of the interruptible load r interrupted in a typical day h period of the mth month of the planning target year; etarThe marginal price of the interruptible load r for a time period t.
The constraint conditions of the optimization model comprise thermal power unit active power upper and lower limit constraints, conventional thermal power unit climbing constraints, node power balance constraints, abandoned wind power quantity constraints, demand response load response range constraints, transmission capacity constraints of the transmission line, maximum extensible line number constraints of a transmission line corridor, node voltage phase angle constraints, thermal power unit active power upper and lower limit constraints in the case of N-1 expected accidents, conventional thermal power unit climbing constraints in the case of N-1 expected accidents, node power balance constraints in the case of N-1 expected accidents, abandoned wind power quantity constraints in the case of N-1 expected accidents, demand response load response range constraints in the case of N-1 expected accidents, transmission capacity constraints in the case of N-1 expected accidents and node voltage phase angle constraints in the case of N-1 expected accidents, the method comprises the following specific steps:
1) and (3) constraining the upper limit and the lower limit of the active power of the thermal power generating unit:
Figure BDA0002612204280000153
wherein the content of the first and second substances,
Figure BDA0002612204280000154
and
Figure BDA0002612204280000155
the active power upper and lower limits of the thermal power generating unit g are respectively set;
Figure BDA0002612204280000156
the active power of the thermal power generating unit g in a typical day h period of the mth month of a planned target year is planned; delta PgThe method shows the depth peak regulation power which can be increased by the flexible modification of the thermal power generating unit g.
2) Conventional thermal power generating unit climbing restraint:
Figure BDA0002612204280000157
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgThe method is the climbing speed variable quantity after the flexibility of the thermal power generating unit g is improved.
3) Node power balance constraint:
Figure BDA0002612204280000161
wherein the content of the first and second substances,
Figure BDA0002612204280000162
the transmission active power of the transmission line corridor l is planned in a typical day h period of the mth month of the target year;
Figure BDA0002612204280000163
the active power predicted value of the wind power plant w in a typical day h period of the mth month of the planning target year;
Figure BDA0002612204280000164
an active power predicted value of a load d in a typical day h period of the mth month of a planning target year;
Figure BDA0002612204280000165
an active power predicted value of an energizable load c for planning a typical day h period of the mth month of a target year;
Figure BDA0002612204280000166
an active power predicted value of an interruptible load r in a typical day h period of the mth month of a planned target year; 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; ciAnd RiRespectively representing the number of excitable loads and interruptible loads on a node i; and nb is the number of the nodes of the power grid.
4) And (5) abandoned wind power constraint:
Figure BDA0002612204280000167
wherein gamma represents a set wind power curtailment rate threshold value.
5) Demand response load response range constraints
Figure BDA0002612204280000168
Figure BDA0002612204280000169
Wherein the content of the first and second substances,
Figure BDA00026122042800001610
and
Figure BDA00026122042800001611
the upper and lower limits of active power of the energizable load c are planned in a typical day h period of the mth month of the target year;
Figure BDA00026122042800001612
and
Figure BDA00026122042800001613
the upper and lower limits of the active power of the interruptible load r are planned for the typical day h period of the mth month of the target year.
6) Transmission capacity constraint of the transmission line:
Figure BDA00026122042800001614
Figure BDA0002612204280000171
wherein, BlThe susceptance of a single-circuit power transmission line on a power transmission line corridor l;
Figure BDA0002612204280000172
a node i voltage phase angle is planned for a typical day h period of the mth month of a target year;
Figure BDA0002612204280000173
the number of the existing transmission lines on the transmission line corridor l is shown; zlRepresenting the number of newly added transmission lines on the transmission line corridor l;
Figure BDA0002612204280000174
the transmission active power of the transmission line corridor l is planned in a typical day h period of the mth month of the target year;
Figure BDA0002612204280000175
the transmission capacity of a single-circuit transmission line on a transmission line corridor l.
7) The maximum extensible line number constraint of the power transmission line corridor is as follows:
Zl≤Zl max,l＝1,L,nl
wherein Z isl maxAnd the maximum number of the extensible transmission lines is the transmission line corridor l.
8) Node voltage phase angle constraint:
Figure BDA0002612204280000176
wherein the content of the first and second substances,
Figure BDA0002612204280000177
the voltage phase angle of the node i is planned in the typical day h period of the mth month of the target year.
9) 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 BDA0002612204280000178
wherein the content of the first and second substances,
Figure BDA0002612204280000179
and in order to plan the active power of the thermal power generating unit g under the expected accident k in the typical day h of the mth month of the target year, nk is the total number of the expected accidents of the power transmission line N-1.
10) N-1 conventional thermal power generating unit climbing restraint under the condition of anticipated accidents:
Figure BDA00026122042800001710
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 BDA00026122042800001711
forecasting the active power of the thermal power generating unit g under an accident k for planning the typical day h period of the mth month of the target year;
Figure BDA00026122042800001712
active power of the thermal power generating unit g under the accident k is expected for planning the typical day h-1 period of the mth month of the target year.
11) N-1 node power balance constraint under the condition of an expected accident:
Figure BDA0002612204280000181
wherein the content of the first and second substances,
Figure BDA0002612204280000182
the transmission active power of a transmission line corridor l under an expected accident k is planned for the mth month typical day h period of a target year;
Figure BDA0002612204280000183
forecasting the active power of the thermal power generating unit g under an accident k for planning the typical day h period of the mth month of the target year;
Figure BDA0002612204280000184
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth month of the planned target year;
Figure BDA0002612204280000185
active power increased by an energizable load c under an expected accident k is planned for a typical day h period of the mth month of the planned target year;
Figure BDA0002612204280000186
active power at which the interruptible load r is interrupted under the accident k is envisioned for planning the typical day h period of the mth month of the target year.
12) N-1 abandon wind power restraint under the condition of anticipated accidents:
Figure BDA0002612204280000187
wherein the content of the first and second substances,
Figure BDA0002612204280000188
electric power curtailment of the wind farm w under the accident k is envisioned for planning the typical day h period of the mth month of the target year.
13) Demand response load response range constraints under N-1 anticipated accident conditions
Figure BDA0002612204280000189
Figure BDA00026122042800001810
14) N-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure BDA00026122042800001811
Figure BDA00026122042800001812
wherein the content of the first and second substances,
Figure BDA00026122042800001813
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure BDA00026122042800001814
indicating a wire outage in the k power line corridor under an expected accident,
Figure BDA0002612204280000191
representing that no line is stopped on the transmission line corridor under the expected accident k;
Figure BDA0002612204280000192
and (4) forecasting the transmission active power of the transmission line corridor l under the accident k for planning the typical day h period of the mth month of the target year.
15) N-1 node voltage phase angle constraint under the expected accident condition:
Figure BDA0002612204280000193
wherein the content of the first and second substances,
Figure BDA0002612204280000194
a voltage phase angle of a node i under an accident k is predicted for a typical day h period of an mth month of a planned target year.
And S3, converting the mixed integer nonlinear constraint condition in the optimization model into a mixed integer linear constraint condition, converting the optimization model into a mixed integer linear programming model, and directly solving by using a mixed integer linear programming algorithm to obtain the final thermal power unit flexibility modification, load demand response mechanism and power transmission network 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 transmission line is converted into an equivalent linear expression:
Figure BDA0002612204280000195
Figure BDA0002612204280000196
wherein M is a constant.
Similarly, a nonlinear constraint expression in transmission capacity constraint of the transmission line under the condition of the N-1 expected accident is converted into an equivalent linear expression:
Figure BDA0002612204280000197
Figure BDA0002612204280000201
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 the final thermal power generating unit flexibility modification, load demand response mechanism and power transmission network planning scheme.
The method comprises the steps of constructing a power grid load-coordinated power transmission network planning optimization model based on sequential characteristics of flexibility modification, a demand side response mechanism, a typical daily load, wind power and the like of a power source side thermal power unit, taking economy as a target, combining reasonable wind abandon allowed by system operation reality, converting mixed integer nonlinear constraint conditions in the optimization model into mixed integer linear constraint conditions, and solving the converted mixed integer linear constraint conditions by adopting a mixed integer linear programming method, so that a final thermal power unit flexibility modification, load demand response mechanism and power transmission network planning scheme is obtained, joint decision of the thermal power unit flexibility modification, the load demand response mechanism and the power transmission network planning is carried out, and the economy of planning investment of a power system is improved while wind power consumption is promoted.
As shown in fig. 2, an embodiment of the present invention further discloses a source-grid-load coordinated power transmission network planning system, where the system includes:
the planning model construction module is used for constructing a power transmission network planning optimization model with source network load coordination, and the optimization model takes the minimum sum of investment cost and operation cost in the whole planning period as a target and contains 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, converting the optimization model into a mixed integer linear programming model, and directly solving by using a mixed integer linear programming algorithm to obtain a final thermal power unit flexibility modification, load demand response mechanism and power transmission network 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 a 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 generating unit capable of participating in flexibility modification and modification cost according to the operation condition of the thermal power generating unit, setting a response quantity range and response price of demand response load participating in demand response, setting a candidate range of a newly-built power transmission line and 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.
The objective function expression of the optimization model is as follows:
Figure BDA0002612204280000211
in the formula, nw is the total number of the wind power plants; ng is the number of conventional thermal power generating units; nl is the number of transmission line corridors planned and constructed; nc and nr represent the number of energizable loads and interruptible loads, respectively; alpha is alphawRepresenting wind curtailment penalty cost;
Figure BDA0002612204280000212
the electric power curtailment of the wind power plant w is planned for the typical day h period of the mth month of the target year; beta is agRepresenting the flexibility modification cost of the thermal power generating unit g; u. ofgIs a binary variable and indicates whether the thermal power generating unit g is subjected to flexible transformation in a planned target year or not, ug1 represents that the thermal power generating unit g is subjected to flexibility modification in a planning target year, and u representsgWhen the thermal power generating unit g is not transformed flexibly in the planning target year, 0 is shown; clInvestment cost for construction of a single-circuit power transmission line on a power transmission corridor l; zlRepresenting the number of newly added transmission lines on the transmission line corridor l;
Figure BDA0002612204280000213
active power which can stimulate the load c to increase for planning the typical day h period of the mth month of the target year;
Figure BDA0002612204280000214
marginal price for the energizable load c;
Figure BDA0002612204280000215
the active power of the interruptible load r interrupted in a typical day h period of the mth month of the planning target year; etarThe marginal price of the interruptible load r for a time period t.
The constraint conditions of the optimization model comprise thermal power unit active power upper and lower limit constraints, conventional thermal power unit climbing constraints, node power balance constraints, abandoned wind power quantity constraints, demand response load response range constraints, transmission capacity constraints of the transmission line, maximum extensible line number constraints of a transmission line corridor, node voltage phase angle constraints, thermal power unit active power upper and lower limit constraints in the case of N-1 expected accidents, conventional thermal power unit climbing constraints in the case of N-1 expected accidents, node power balance constraints in the case of N-1 expected accidents, abandoned wind power quantity constraints in the case of N-1 expected accidents, demand response load response range constraints in the case of N-1 expected accidents, transmission capacity constraints in the case of N-1 expected accidents and node voltage phase angle constraints in the case of N-1 expected accidents, the method comprises the following specific steps:
1) and (3) constraining the upper limit and the lower limit of the active power of the thermal power generating unit:
Figure BDA0002612204280000221
wherein the content of the first and second substances,
Figure BDA0002612204280000222
and
Figure BDA0002612204280000223
the active power upper and lower limits of the thermal power generating unit g are respectively set;
Figure BDA0002612204280000224
the active power of the thermal power generating unit g in a typical day h period of the mth month of a planned target year is planned; delta PgThe method shows the depth peak regulation power which can be increased by the flexible modification of the thermal power generating unit g.
2) Conventional thermal power generating unit climbing restraint:
Figure BDA0002612204280000225
wherein R isgThe ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ RgThe method is the climbing speed variable quantity after the flexibility of the thermal power generating unit g is improved.
3) Node power balance constraint:
Figure BDA0002612204280000226
wherein the content of the first and second substances,
Figure BDA0002612204280000227
for planning power transmission in typical day h period of mth month of target yearThe transmission active power of the line corridor l;
Figure BDA0002612204280000228
the active power predicted value of the wind power plant w in a typical day h period of the mth month of the planning target year;
Figure BDA0002612204280000229
an active power predicted value of a load d in a typical day h period of the mth month of a planning target year;
Figure BDA00026122042800002210
an active power predicted value of an energizable load c for planning a typical day h period of the mth month of a target year;
Figure BDA00026122042800002211
an active power predicted value of an interruptible load r in a typical day h period of the mth month of a planned target year; 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; ciAnd RiRespectively representing the number of excitable loads and interruptible loads on a node i; and nb is the number of the nodes of the power grid.
4) And (5) abandoned wind power constraint:
Figure BDA0002612204280000231
wherein gamma represents a set wind power curtailment rate threshold value.
5) Demand response load response range constraints
Figure BDA0002612204280000232
Figure BDA0002612204280000233
Wherein the content of the first and second substances,
Figure BDA0002612204280000234
and
Figure BDA0002612204280000235
the upper and lower limits of active power of the energizable load c are planned in a typical day h period of the mth month of the target year;
Figure BDA0002612204280000236
and
Figure BDA0002612204280000237
the upper and lower limits of the active power of the interruptible load r are planned for the typical day h period of the mth month of the target year.
6) Transmission capacity constraint of the transmission line:
Figure BDA0002612204280000238
Figure BDA0002612204280000239
wherein, BlThe susceptance of a single-circuit power transmission line on a power transmission line corridor l;
Figure BDA00026122042800002310
a node i voltage phase angle is planned for a typical day h period of the mth month of a target year;
Figure BDA00026122042800002311
the number of the existing transmission lines on the transmission line corridor l is shown; zlRepresenting the number of newly added transmission lines on the transmission line corridor l;
Figure BDA00026122042800002312
the transmission active power of the transmission line corridor l is planned in a typical day h period of the mth month of the target year;
Figure BDA00026122042800002313
the transmission capacity of a single-circuit transmission line on a transmission line corridor l.
7) The maximum extensible line number constraint of the power transmission line corridor is as follows:
Zl≤Zl max,l＝1,L,nl
wherein Z isl maxAnd the maximum number of the extensible transmission lines is the transmission line corridor l.
8) Node voltage phase angle constraint:
Figure BDA0002612204280000241
wherein the content of the first and second substances,
Figure BDA0002612204280000242
the voltage phase angle of the node i is planned in the typical day h period of the mth month of the target year.
9) 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 BDA0002612204280000243
wherein the content of the first and second substances,
Figure BDA0002612204280000244
and in order to plan the active power of the thermal power generating unit g under the expected accident k in the typical day h of the mth month of the target year, nk is the total number of the expected accidents of the power transmission line N-1.
10) N-1 conventional thermal power generating unit climbing restraint under the condition of anticipated accidents:
Figure BDA0002612204280000245
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 BDA0002612204280000246
forecasting the active power of the thermal power generating unit g under an accident k for planning the typical day h period of the mth month of the target year;
Figure BDA0002612204280000247
active power of the thermal power generating unit g under the accident k is expected for planning the typical day h-1 period of the mth month of the target year.
11) N-1 node power balance constraint under the condition of an expected accident:
Figure BDA0002612204280000248
wherein the content of the first and second substances,
Figure BDA0002612204280000249
the transmission active power of a transmission line corridor l under an expected accident k is planned for the mth month typical day h period of a target year;
Figure BDA00026122042800002410
forecasting the active power of the thermal power generating unit g under an accident k for planning the typical day h period of the mth month of the target year;
Figure BDA00026122042800002411
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth month of the planned target year;
Figure BDA00026122042800002412
active power increased by an energizable load c under an expected accident k is planned for a typical day h period of the mth month of the planned target year;
Figure BDA00026122042800002413
active power at which the interruptible load r is interrupted under the accident k is envisioned for planning the typical day h period of the mth month of the target year.
12) N-1 abandon wind power restraint under the condition of anticipated accidents:
Figure BDA0002612204280000251
wherein the content of the first and second substances,
Figure BDA0002612204280000252
for planning a target yearThe m-th month typical day h period envisions the abandoned electric power of the wind farm w under the accident k.
13) Demand response load response range constraints under N-1 anticipated accident conditions
Figure BDA0002612204280000253
Figure BDA0002612204280000254
14) N-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure BDA0002612204280000255
Figure BDA0002612204280000256
wherein the content of the first and second substances,
Figure BDA0002612204280000257
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure BDA0002612204280000258
indicating a wire outage in the k power line corridor under an expected accident,
Figure BDA0002612204280000259
representing that no line is stopped on the transmission line corridor under the expected accident k;
Figure BDA00026122042800002510
and (4) forecasting the transmission active power of the transmission line corridor l under the accident k for planning the typical day h period of the mth month of the target year.
15) N-1 node voltage phase angle constraint under the expected accident condition:
Figure BDA00026122042800002511
wherein the content of the first and second substances,
Figure BDA00026122042800002512
a voltage phase angle of a node i under an accident k is predicted for a typical day h period of an mth month of a planned target year.
And converting the mixed integer nonlinear constraint condition in the optimization model into a mixed integer linear constraint condition, converting the optimization model into a mixed integer linear programming model, and directly solving by using a mixed integer linear programming algorithm to obtain a final thermal power unit flexibility modification, load demand response mechanism and power transmission network 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 transmission line is converted into an equivalent linear expression:
Figure BDA0002612204280000261
Figure BDA0002612204280000262
wherein M is a constant.
Similarly, a nonlinear constraint expression in transmission capacity constraint of the transmission line under the condition of the N-1 expected accident is converted into an equivalent linear expression:
Figure BDA0002612204280000263
Figure BDA0002612204280000264
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 the final thermal power generating unit flexibility modification, load demand response mechanism and power transmission network planning 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 source-grid-load coordinated power transmission network planning method is characterized by comprising the following operations:
constructing a power transmission network planning optimization model with source network load coordination, wherein the optimization model takes the minimum sum of investment cost and operation cost in the whole planning period as a target and comprises constraint conditions;
and converting the mixed integer nonlinear constraint condition in the optimization model into a mixed integer linear constraint condition, converting the optimization model into a mixed integer linear programming model, and directly solving by using a mixed integer linear programming algorithm to obtain a final thermal power unit flexibility modification, load demand response mechanism and power transmission network programming scheme.
2. The method for power transmission network planning with source-grid-load coordination according to claim 1, further comprising inputting parameters of a traditional thermal power unit and parameters of a current power transmission element of a power network before constructing the model, inputting maximum load, maximum wind power, and 24-hour load power and wind power change per typical day of each year in a planning period, setting a range of the thermal power unit capable of participating in flexibility modification and a modification cost according to the operation condition of the thermal power unit, setting a response quantity range and a response price of demand response load participating in demand response, 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 curtailment cost coefficient and a system-allowable curtailment rate.
3. The source-grid-load coordinated power transmission network planning method according to claim 1, wherein an objective function expression of the optimization model is as follows:
Figure FDA0002612204270000011
in the formula, nw is the total number of the wind power plants; ng is the number of conventional thermal power generating units; nl is the number of transmission line corridors planned and constructed; nc and nr represent the number of energizable loads and interruptible loads, respectively; alpha is alphawRepresenting wind curtailment penalty cost;
Figure FDA0002612204270000012
the electric power curtailment of the wind power plant w is planned for the typical day h period of the mth month of the target year; beta is agRepresenting the flexibility modification cost of the thermal power generating unit g; u. ofgIs a binary variable and indicates whether the thermal power generating unit g is subjected to flexible transformation in a planned target year or not, ug1 represents that the thermal power generating unit g is subjected to flexibility modification in a planning target year, and u representsgWhen the thermal power generating unit g is not transformed flexibly in the planning target year, 0 is shown; clInvestment cost for construction of a single-circuit power transmission line on a power transmission corridor l; zlRepresenting the number of newly added transmission lines on the transmission line corridor l;
Figure FDA0002612204270000021
active power which can stimulate the load c to increase for planning the typical day h period of the mth month of the target year;
Figure FDA00026122042700000211
marginal price for the energizable load c;
Figure FDA0002612204270000022
the active power of the interruptible load r interrupted in a typical day h period of the mth month of the planning target year; etarThe marginal price of the interruptible load r for a time period t.
4. The source-grid-load coordinated power transmission network planning method according to claim 1, wherein the constraint condition comprises:
1) and (3) constraining the upper limit and the lower limit of the active power of the thermal power generating unit:
Figure FDA0002612204270000023
wherein the content of the first and second substances,
Figure FDA0002612204270000024
and
Figure FDA0002612204270000025
the active power upper and lower limits of the thermal power generating unit g are respectively set;
Figure FDA0002612204270000026
the active power of the thermal power generating unit g in a typical day h period of the mth month of a planned target year is planned; 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;
2) conventional thermal power generating unit climbing restraint:
Figure FDA0002612204270000027
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;
3) node power balance constraint:
Figure FDA0002612204270000028
wherein, Pl m,hThe transmission active power of the transmission line corridor l is planned in a typical day h period of the mth month of the target year;
Figure FDA0002612204270000029
the active power predicted value of the wind power plant w in a typical day h period of the mth month of the planning target year;
Figure FDA00026122042700000210
an active power predicted value of a load d in a typical day h period of the mth month of a planning target year;
Figure FDA0002612204270000031
an active power predicted value of an energizable load c for planning a typical day h period of the mth month of a target year;
Figure FDA0002612204270000032
an active power predicted value of an interruptible load r in a typical day h period of the mth month of a planned target year; 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; ciAnd RiRespectively representing the number of excitable loads and interruptible loads on a node i; nb is the number of the nodes of the power grid;
4) and (5) abandoned wind power constraint:
Figure FDA0002612204270000033
Figure FDA0002612204270000034
wherein gamma represents a set wind power curtailment rate threshold value;
5) demand response load response range constraints
Figure FDA0002612204270000035
Figure FDA0002612204270000036
Wherein the content of the first and second substances,
Figure FDA0002612204270000037
and
Figure FDA0002612204270000038
the upper and lower limits of active power of the energizable load c are planned in a typical day h period of the mth month of the target year;
Figure FDA0002612204270000039
and
Figure FDA00026122042700000310
the upper and lower limits of the active power of the interruptible load r in a typical day h period of the mth month of the planned target year are planned;
6) transmission capacity constraint of the transmission line:
Figure FDA00026122042700000311
Figure FDA00026122042700000312
wherein, BlThe susceptance of a single-circuit power transmission line on a power transmission line corridor l;
Figure FDA00026122042700000313
a node i voltage phase angle is planned for a typical day h period of the mth month of a target year;
Figure FDA00026122042700000314
the number of the existing transmission lines on the transmission line corridor l is shown; zlRepresenting the number of newly added transmission lines on the transmission line corridor l; pl m,hThe transmission active power of the transmission line corridor l is planned in a typical day h period of the mth month of the target year; pl maxThe transmission capacity of a single-circuit transmission line on a transmission line corridor l;
7) the maximum extensible line number constraint of the power transmission line corridor is as follows:
Zl≤Zl max,l＝1,L,nl
wherein Z isl maxThe maximum number of the extensible power transmission lines is 1;
8) node voltage phase angle constraint:
Figure FDA0002612204270000041
wherein the content of the first and second substances,
Figure FDA0002612204270000042
a node i voltage phase angle is planned for a typical day h period of the mth month of a target year;
9) 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 FDA0002612204270000043
wherein the content of the first and second substances,
Figure FDA0002612204270000044
in order to plan the active power of the thermal power generating unit g under the expected accident k in the typical day h of the mth month of the target year, nk is the total number of the expected accidents of the power transmission line N-1;
10) n-1 conventional thermal power generating unit climbing restraint under the condition of anticipated accidents:
Figure FDA0002612204270000045
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 FDA0002612204270000046
forecasting the active power of the thermal power generating unit g under an accident k for planning the typical day h period of the mth month of the target year;
Figure FDA0002612204270000047
active power of the thermal power generating unit g under an expected accident k is planned in a typical day h-1 period of the mth month of a planned target year;
11) n-1 node power balance constraint under the condition of an expected accident:
Figure FDA0002612204270000048
wherein, Pl m,h,kThe transmission active power of a transmission line corridor l under an expected accident k is planned for the mth month typical day h period of a target year;
Figure FDA0002612204270000051
forecasting the active power of the thermal power generating unit g under an accident k for planning the typical day h period of the mth month of the target year;
Figure FDA0002612204270000052
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth month of the planned target year;
Figure FDA0002612204270000053
active power increased by an energizable load c under an expected accident k is planned for a typical day h period of the mth month of the planned target year;
Figure FDA0002612204270000054
active power of an interruptible load r interrupted under an expected accident k for planning a typical day h period of an mth month of a target year;
12) n-1 abandon wind power restraint under the condition of anticipated accidents:
Figure FDA0002612204270000055
wherein the content of the first and second substances,
Figure FDA0002612204270000056
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth month of the planned target year;
13) demand response load response range constraints under N-1 anticipated accident conditions
Figure FDA0002612204270000057
Figure FDA0002612204270000058
14) N-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure FDA0002612204270000059
Figure FDA00026122042700000510
wherein the content of the first and second substances,
Figure FDA00026122042700000511
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure FDA00026122042700000512
indicating a wire outage in the k power line corridor under an expected accident,
Figure FDA00026122042700000513
representing that no line is stopped on the transmission line corridor under the expected accident k; pl m,h,kThe transmission active power of a transmission line corridor l under an expected accident k is planned for the mth month typical day h period of a target year;
15) n-1 node voltage phase angle constraint under the expected accident condition:
Figure FDA00026122042700000514
wherein the content of the first and second substances,
Figure FDA0002612204270000061
a voltage phase angle of a node i under an accident k is predicted for a typical day h period of an mth month of a planned target year.
5. The source-grid-load coordinated power transmission network planning method according to claim 4, wherein the non-linear constraint expression in the transmission capacity constraint of the transmission line is converted into an equivalent linear expression:
Figure FDA0002612204270000062
Figure FDA0002612204270000063
wherein M is a constant.
6. The source-grid-load coordinated power transmission network planning method 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:
Figure FDA0002612204270000064
Figure FDA0002612204270000065
7. a source-grid-load coordinated power transmission network planning system, the system comprising:
the planning model construction module is used for constructing a power transmission network planning optimization model with source network load coordination, and the optimization model takes the minimum sum of investment cost and operation cost in the whole planning period as a target and contains 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, converting the optimization model into a mixed integer linear programming model, and directly solving by using a mixed integer linear programming algorithm to obtain a final thermal power unit flexibility modification, load demand response mechanism and power transmission network programming scheme.
8. The power transmission network planning system with source-grid-load coordination according to claim 7, further comprising a model parameter setting module, configured to input parameters of a conventional thermal power unit and parameters of power transmission elements of a current power grid before building a model, input maximum load, maximum wind power, and 24-hour load power and wind power change per typical day of each year in a planning period, set a range of the thermal power unit that can participate in flexibility modification and set modification costs according to an operation condition of the thermal power unit, set a response quantity range and a response price of demand response load participating in demand response, set a range of a candidate new power transmission line and set construction costs according to a corridor condition of the power transmission line, and set a wind curtailment cost coefficient and a power curtailment rate allowed by the system.
9. The source-grid-load coordinated power transmission network planning system according to claim 7, wherein the objective function expression of the optimization model is as follows:
Figure FDA0002612204270000071
in the formula, nw is the total number of the wind power plants; ng is the number of conventional thermal power generating units; nl is the number of transmission line corridors planned and constructed; nc and nr represent the number of energizable loads and interruptible loads, respectively; alpha is alphawRepresenting wind curtailment penalty cost;
Figure FDA0002612204270000072
the electric power curtailment of the wind power plant w is planned for the typical day h period of the mth month of the target year; beta is agRepresenting the flexibility modification cost of the thermal power generating unit g; u. ofgIs a binary variable and represents the thermal power generating unit in a planned target yearg if flexibility modification is to be performed, ug1 represents that the thermal power generating unit g is subjected to flexibility modification in a planning target year, and u representsgWhen the thermal power generating unit g is not transformed flexibly in the planning target year, 0 is shown; clInvestment cost for construction of a single-circuit power transmission line on a power transmission corridor l; zlRepresenting the number of newly added transmission lines on the transmission line corridor l;
Figure FDA0002612204270000073
active power which can stimulate the load c to increase for planning the typical day h period of the mth month of the target year;
Figure FDA0002612204270000075
marginal price for the energizable load c;
Figure FDA0002612204270000074
the active power of the interruptible load r interrupted in a typical day h period of the mth month of the planning target year; etarThe marginal price of the interruptible load r for a time period t.
10. The source-grid-load coordinated power transmission network planning system according to claim 7, wherein the constraint condition comprises:
1) and (3) constraining the upper limit and the lower limit of the active power of the thermal power generating unit:
Figure FDA0002612204270000081
wherein the content of the first and second substances,
Figure FDA0002612204270000082
and
Figure FDA0002612204270000083
the active power upper and lower limits of the thermal power generating unit g are respectively set;
Figure FDA0002612204270000084
for planning typical day h of mth month of target yearActive power of the thermal power generating unit g is segmented; 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;
2) conventional thermal power generating unit climbing restraint:
Figure FDA0002612204270000085
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;
3) node power balance constraint:
Figure FDA0002612204270000086
wherein, Pl m,hThe transmission active power of the transmission line corridor l is planned in a typical day h period of the mth month of the target year;
Figure FDA0002612204270000087
the active power predicted value of the wind power plant w in a typical day h period of the mth month of the planning target year;
Figure FDA0002612204270000088
an active power predicted value of a load d in a typical day h period of the mth month of a planning target year;
Figure FDA0002612204270000089
an active power predicted value of an energizable load c for planning a typical day h period of the mth month of a target year;
Figure FDA00026122042700000810
an active power predicted value of an interruptible load r in a typical day h period of the mth month of a planned target year; 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 of thermal power generating units, wind power plants and loads on a node iCounting; ciAnd RiRespectively representing the number of excitable loads and interruptible loads on a node i; nb is the number of the nodes of the power grid;
4) and (5) abandoned wind power constraint:
Figure FDA0002612204270000091
Figure FDA0002612204270000092
wherein gamma represents a set wind power curtailment rate threshold value;
5) demand response load response range constraints
Figure FDA0002612204270000093
Figure FDA0002612204270000094
Wherein the content of the first and second substances,
Figure FDA0002612204270000095
and
Figure FDA0002612204270000096
the upper and lower limits of active power of the energizable load c are planned in a typical day h period of the mth month of the target year;
Figure FDA0002612204270000097
and
Figure FDA0002612204270000098
the upper and lower limits of the active power of the interruptible load r in a typical day h period of the mth month of the planned target year are planned;
6) transmission capacity constraint of the transmission line:
Figure FDA0002612204270000099
Figure FDA00026122042700000910
wherein, BlThe susceptance of a single-circuit power transmission line on a power transmission line corridor l;
Figure FDA00026122042700000911
a node i voltage phase angle is planned for a typical day h period of the mth month of a target year;
Figure FDA00026122042700000912
the number of the existing transmission lines on the transmission line corridor l is shown; zlRepresenting the number of newly added transmission lines on the transmission line corridor l; pl m,hThe transmission active power of the transmission line corridor l is planned in a typical day h period of the mth month of the target year; pl maxThe transmission capacity of a single-circuit transmission line on a transmission line corridor l;
7) the maximum extensible line number constraint of the power transmission line corridor is as follows:
Zl≤Zl max,l＝1,L,nl
wherein Z isl maxThe maximum number of the extensible power transmission lines is 1;
8) node voltage phase angle constraint:
Figure FDA0002612204270000101
wherein the content of the first and second substances,
Figure FDA0002612204270000102
a node i voltage phase angle is planned for a typical day h period of the mth month of a target year;
9) 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 FDA0002612204270000103
wherein the content of the first and second substances,
Figure FDA0002612204270000104
in order to plan the active power of the thermal power generating unit g under the expected accident k in the typical day h of the mth month of the target year, nk is the total number of the expected accidents of the power transmission line N-1;
10) n-1 conventional thermal power generating unit climbing restraint under the condition of anticipated accidents:
Figure FDA0002612204270000105
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 FDA0002612204270000106
forecasting the active power of the thermal power generating unit g under an accident k for planning the typical day h period of the mth month of the target year;
Figure FDA0002612204270000107
active power of the thermal power generating unit g under an expected accident k is planned in a typical day h-1 period of the mth month of a planned target year;
11) n-1 node power balance constraint under the condition of an expected accident:
Figure FDA0002612204270000108
wherein, Pl m,h,kThe transmission active power of a transmission line corridor l under an expected accident k is planned for the mth month typical day h period of a target year;
Figure FDA0002612204270000109
forecasting the active power of the thermal power generating unit g under an accident k for planning the typical day h period of the mth month of the target year;
Figure FDA00026122042700001010
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth month of the planned target year;
Figure FDA00026122042700001011
active power increased by an energizable load c under an expected accident k is planned for a typical day h period of the mth month of the planned target year;
Figure FDA00026122042700001012
active power of an interruptible load r interrupted under an expected accident k for planning a typical day h period of an mth month of a target year;
12) n-1 abandon wind power restraint under the condition of anticipated accidents:
Figure FDA0002612204270000111
wherein the content of the first and second substances,
Figure FDA0002612204270000112
forecasting the abandoned electric power of the wind power plant w under the accident k for the typical day h period of the mth month of the planned target year;
13) demand response load response range constraints under N-1 anticipated accident conditions
Figure FDA0002612204270000113
Figure FDA0002612204270000114
14) N-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure FDA0002612204270000115
Figure FDA0002612204270000116
wherein the content of the first and second substances,
Figure FDA0002612204270000117
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure FDA0002612204270000118
indicating a wire outage in the k power line corridor under an expected accident,
Figure FDA0002612204270000119
representing that no line is stopped on the transmission line corridor under the expected accident k; pl m,h,kThe transmission active power of a transmission line corridor l under an expected accident k is planned for the mth month typical day h period of a target year;
15) n-1 node voltage phase angle constraint under the expected accident condition:
Figure FDA00026122042700001110
wherein the content of the first and second substances,
Figure FDA00026122042700001111
a voltage phase angle of a node i under an accident k is predicted for a typical day h period of an mth month of a planned target year.
CN202010757955.7A 2020-07-31 2020-07-31 Source-grid-load cooperative power transmission network planning method and system Active CN111799793B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010757955.7A CN111799793B (en) 2020-07-31 2020-07-31 Source-grid-load cooperative power transmission network planning method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010757955.7A CN111799793B (en) 2020-07-31 2020-07-31 Source-grid-load cooperative power transmission network planning method and system

Publications (2)

Publication Number Publication Date
CN111799793A true CN111799793A (en) 2020-10-20
CN111799793B CN111799793B (en) 2022-08-30

Family

ID=72828051

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010757955.7A Active CN111799793B (en) 2020-07-31 2020-07-31 Source-grid-load cooperative power transmission network planning method and system

Country Status (1)

Country Link
CN (1) CN111799793B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113722940A (en) * 2021-11-04 2021-11-30 国网湖北省电力有限公司经济技术研究院 Power grid unit combination and technical improvement plan combined optimization method
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 (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN109858774A (en) * 2019-01-09 2019-06-07 燕山大学 Improve the source net lotus planing method of security of system and harmony
CN110852565A (en) * 2019-10-10 2020-02-28 国家电网有限公司 Power transmission network frame planning method considering different functional attributes
WO2020063144A1 (en) * 2018-09-30 2020-04-02 中国电力科学研究院有限公司 Method and system for evaluating energy delivery capacity in flexible dc electrical grid

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
WO2020063144A1 (en) * 2018-09-30 2020-04-02 中国电力科学研究院有限公司 Method and system for evaluating energy delivery capacity in flexible dc electrical grid
CN109858774A (en) * 2019-01-09 2019-06-07 燕山大学 Improve the source net lotus planing method of security of system and harmony
CN110852565A (en) * 2019-10-10 2020-02-28 国家电网有限公司 Power transmission network frame planning method considering different functional attributes

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
樊金柱 等: "考虑网源协同的输电网适应性扩展规划", 《电网技术》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113722940A (en) * 2021-11-04 2021-11-30 国网湖北省电力有限公司经济技术研究院 Power grid unit combination and technical improvement plan combined optimization method
CN113722940B (en) * 2021-11-04 2022-02-08 国网湖北省电力有限公司经济技术研究院 Power grid unit combination and technical improvement plan combined optimization method
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
CN111799793B (en) 2022-08-30

Similar Documents

Publication Publication Date Title
Xu et al. Carbon emission reduction and reliable power supply equilibrium based daily scheduling towards hydro-thermal-wind generation system: A perspective from China
Franco et al. A network flow model for short-term hydro-dominated hydrothermal scheduling problems
Chen et al. Key technologies for integration of multitype renewable energy sources—Research on multi-timeframe robust scheduling/dispatch
CN109767078B (en) Multi-type power supply maintenance arrangement method based on mixed integer programming
Ahmadi et al. Multi-objective economic emission dispatch considering combined heat and power by normal boundary intersection method
Mc Garrigle et al. Quantifying the value of improved wind energy forecasts in a pool-based electricity market
Xie et al. Optimal capacity and type planning of generating units in a bundled wind–thermal generation system
Delarue et al. Adaptive mixed-integer programming unit commitment strategy for determining the value of forecasting
CN104809545B (en) A kind of virtual plant runs modeling method
CN111799793B (en) Source-grid-load cooperative power transmission network planning method and system
CN105243600A (en) Grid power generation adjustment method
Kumar et al. Application of BARON solver for solution of cost based unit commitment problem
CN111799841B (en) Combined decision method and system for thermal power generating unit flexibility transformation and power transmission planning
Lin et al. An interval parameter optimization model for sustainable power systems planning under uncertainty
CN113363976A (en) Scene graph-based mid-term optimized scheduling method for wind, light and water complementary power generation system
Sun et al. Interval mixed-integer programming for daily unit commitment and dispatch incorporating wind power
Contaxis et al. Optimal power flow considering operation of wind parks and pump storage hydro units under large scale integration of renewable energy sources
CN112308411A (en) Comprehensive energy station random planning method and system based on dynamic carbon transaction model
CN111799842B (en) Multi-stage power transmission network planning method and system considering flexibility of thermal power generating unit
CN108039739B (en) Dynamic random economic dispatching method for active power distribution network
Fang et al. Parallel improved DPSA algorithm for medium-term optimal scheduling of large-scale cascade hydropower plants
Wang et al. Multi-agent interaction of source, load and storage to realize peak shaving and valley filling under the guidance of the market mechanism
Lin et al. Three-stage optimization model of medium and long-term contract energy decomposition considering equilibrium of carbon emissions
Liming et al. Research on the Hydropower Coupling-Based Hydropower Station Scheduling Optimization Model
Lin et al. Real-time dispatch with flexible ramping products in a power system with high penetration of renewable energy generation

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