CN111799793B - 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
CN111799793B
CN111799793B CN202010757955.7A CN202010757955A CN111799793B CN 111799793 B CN111799793 B CN 111799793B CN 202010757955 A CN202010757955 A CN 202010757955A CN 111799793 B CN111799793 B CN 111799793B
Authority
CN
China
Prior art keywords
power
period
load
typical day
generating unit
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010757955.7A
Other languages
Chinese (zh)
Other versions
CN111799793A (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 problem that the total installed quantity of a thermal power generating unit is rich but the flexibility is insufficient is solved, and by the end of 2019, the flexibility of a traditional thermal power generating unit is improved by a coal-fired thermal power generating unit 1040GW in China, so that the flexibility space of a source side can be improved to a great extent. 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 power demand side in real time. 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 optimized power grid planning method is to improve the traditional power grid planning method so that the optimized 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 the 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 alpha w Representing wind abandon 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 a g Representing the flexibility modification cost of the thermal power generating unit g; u. of g Is a binary variable and indicates whether the thermal power generating unit g is subjected to flexible transformation in a planned target year or not, u g 1 represents that the thermal power generating unit g is subjected to flexibility modification in a planning target year, and u represents g When the thermal power generating unit g is not transformed flexibly in the planning target year, 0 is shown; c l Investment cost for construction of a single-circuit power transmission line on a power transmission corridor l; z l Indicating new addition to line corridor lThe number of the transmission lines is counted;
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
the active power of the interruptible load r interrupted in a typical day h period of the mth month of the planning target year; eta r The marginal price of the load r can be interrupted for a 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,
Figure BDA0002612204280000037
and
Figure BDA0002612204280000038
respectively representing the upper limit and the lower limit of the active power of the thermal power generating unit g;
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 P g The method comprises the steps that the depth peak regulation power which can be increased by modifying the thermal power unit g through flexibility is shown;
2) conventional thermal power generating unit climbing restraint:
Figure BDA0002612204280000041
wherein R is g The ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ R g The 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,
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
planning an active power predicted value of a load d in a typical day h period of the mth month of a 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; s i And E i The total number of the power transmission lines with the node i as a head node and the tail end node respectively; g i 、W i And D i Respectively representing the total number of thermal power generating units, wind power plants and loads on the node i; c i And R i Respectively 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) abandon 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,
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, B l The susceptance of a single-circuit power transmission line on a power transmission line corridor l;
Figure BDA0002612204280000057
planning a voltage phase angle of a node i in 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:
Z l ≤Z l max ,l＝1,L,nl
wherein Z is l max The maximum number of the extensible power transmission lines is 1;
8) node voltage phase angle constraint:
Figure BDA00026122042800000511
wherein,
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) and (2) constraining upper and lower limits of active power of the thermal power generating unit under the condition of an N-1 anticipated accident:
Figure BDA00026122042800000513
wherein,
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 is g The ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ R g The 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
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 expected accident:
Figure BDA0002612204280000064
wherein,
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
active power of the thermal power generating unit g under an expected accident k for planning a typical day h period of the mth month of a 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
for planning the target yearActive power which can stimulate the load c to increase under an expected accident k in a typical day h period of m months;
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,
Figure BDA00026122042800000611
forecasting the abandoned electric power of the wind power plant w under the accident k for the planning of the mth month, the typical day h period;
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 expected accidents:
Figure BDA0002612204280000071
Figure BDA0002612204280000072
wherein,
Figure BDA0002612204280000073
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure BDA0002612204280000074
indicating that the wire in the k line corridor is shut down in anticipation of the 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,
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 nonlinear 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 source network load cooperative power transmission network planning optimization model, and 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 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 power transmission elements of the power grid in the current situation before the model is built, 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 flexible 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 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 gaugeDividing the number of constructed power transmission line corridors; nc and nr represent the number of energizable loads and interruptible loads respectively; alpha is alpha w Representing 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 a g Representing the flexibility modification cost of the thermal power generating unit g; u. of g Is a binary variable and indicates whether the thermal power generating unit g is subjected to flexible transformation in a planned target year or not, u g 1 represents the thermal power generating unit g is subjected to flexibility transformation in a planned target year, u g When the thermal power generating unit g is not subjected to flexibility transformation in a planned target year, the thermal power generating unit g is represented as 0; c l Investment cost for construction of a single-circuit power transmission line on a power transmission corridor l; z is a linear or branched member l Representing 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; eta r The 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,
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 P g The 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 is g The ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ R g The 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,
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; s i And E i The total number of the power transmission lines with the node i as a head node and the tail end node respectively; g i 、W i And D i Respectively representing the total number of thermal power generating units, wind power plants and loads on the node i; c i And R i Respectively 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) abandon 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,
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, B l The susceptance of a single-circuit power transmission line on a power transmission line corridor l is obtained;
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; z is a linear or branched member l Representing 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:
Z l ≤Z l max ,l＝1,L,nl
wherein Z is l max The maximum number of the extensible power transmission lines is 1;
8) node voltage phase angle constraint:
Figure BDA0002612204280000111
wherein,
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) and (2) constraining upper and lower limits of active power of the thermal power generating unit under the condition of an N-1 anticipated accident:
Figure BDA0002612204280000113
wherein,
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 is g The ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ R g The climbing rate variable quantity after the flexibility of the thermal power generating unit g is modified;
Figure BDA0002612204280000116
active power of the thermal power generating unit g under an expected accident k for planning a typical day h period of the mth month of a 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,
Figure BDA0002612204280000119
the active power transmission of the transmission line corridor l under the forecast accident k in the typical day h period of the mth month of the planned target year is planned;
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,
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,
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,
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 source-network-load-coordinated power transmission network planning method and system provided by the 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 takes the minimum sum of investment cost and operation cost in the whole planning period as a target and contains constraint conditions.
The objective function expression of the optimization model is:
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 alpha w Representing 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 a g Representing the flexibility modification cost of the thermal power generating unit g; u. of g Is a binary variable and indicates whether the thermal power generating unit g is subjected to flexible transformation in a planned target year or not, u g 1 represents the thermal power generating unit g is subjected to flexibility transformation in a planned target year, u g When the thermal power generating unit g does not operate in the planning target year as 0Modifying flexibility; c l Investment cost for construction of a single-circuit power transmission line on a power transmission corridor l; z is a linear or branched member l Representing 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
active power of interruptible load r interrupted for planning a typical day h period of the mth month of the target year; eta r The 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) the thermal power unit active power upper and lower limit constraint:
Figure BDA0002612204280000153
wherein,
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 P g The 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 is g The ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ R g The climbing rate variable quantity after flexibility transformation of the thermal power generating unit g is obtained.
3) Node power balance constraint:
Figure BDA0002612204280000161
wherein,
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 excitable 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 for planning a typical day h period of the mth month of a target year; s. the i And E i The total number of the power transmission lines taking the node i as a head node and the tail end node respectively; g i 、W i And D i Respectively representing the total number of thermal power generating units, wind power plants and loads on the node i; c i And R i Respectively 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,
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
load interruptible for planning typical day h period of mth month of target yearThe upper and lower limits of the active power of r.
6) Transmission capacity constraint of the transmission line:
Figure BDA00026122042800001614
Figure BDA0002612204280000171
wherein, B l The 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; z l Representing the number of newly added transmission lines on the transmission line corridor l;
Figure BDA0002612204280000174
the transmission active power of the power transmission line corridor l in the typical day h period of the mth month of the planned target year is planned;
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:
Z l ≤Z l max ,l＝1,L,nl
wherein Z is l max And the maximum number of the extensible transmission lines is the transmission line corridor l.
8) Node voltage phase angle constraint:
Figure BDA0002612204280000176
wherein,
Figure BDA0002612204280000177
the node i voltage phase angle is planned for 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,
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 is g The ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ R g The method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure BDA00026122042800001711
active power of the thermal power generating unit g under an expected accident k for planning a typical day h period of the mth month of a 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,
Figure BDA0002612204280000182
the active power transmission of the transmission line corridor l under the forecast accident k in the typical day h period of the mth month of the planned target year is planned;
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,
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,
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,
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 the 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 source network load cooperative power transmission network planning optimization model, and 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 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 the 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 power transmission elements of a 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 typical day of each year 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 running 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 new 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 alpha w Representing 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 a g Representing the flexibility modification cost of the thermal power generating unit g; u. of g Is a binary variable and indicates whether the thermal power generating unit g is subjected to flexible transformation in a planned target year or not, u g 1 represents that the thermal power generating unit g is subjected to flexibility modification in a planning target year, and u represents g When the thermal power generating unit g is not transformed flexibly in the planning target year, 0 is shown; c l Investment cost for construction of a single-circuit power transmission line on a power transmission corridor l; z l Representing 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; eta r The marginal price of the load r can be interrupted for a 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,
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 P g The method shows the depth peak regulation power which can be increased by the flexible modification of the thermal power generating unit g.
2) Climbing restraint of a conventional thermal power generating unit:
Figure BDA0002612204280000225
wherein R is g The ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ R g The 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,
Figure BDA0002612204280000227
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 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; s i And E i The total number of the power transmission lines with the node i as a head node and the tail end node respectively; g i 、W i And D i Respectively representing the total number of thermal power generating units, wind power plants and loads on the node i; c i And R i Respectively 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,
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, B l The 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; z l Representing the number of newly added transmission lines on the transmission line corridor l;
Figure BDA00026122042800002312
for planning the typical day of the mth month of the target yearThe transmission active power of the transmission line corridor l in the h period;
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:
Z l ≤Z l max ,l＝1,L,nl
wherein Z is l max And the maximum number of the extensible transmission lines is the transmission line corridor l.
8) Node voltage phase angle constraint:
Figure BDA0002612204280000241
wherein,
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,
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 is g The ramp rate is the ramp rate of the thermal power generating unit g without modification; deltaR g The 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,
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
active power of the thermal power generating unit g under an expected accident k for planning a typical day h period of the mth month of a 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 of an interruptible load r interrupted under an expected accident k for a planning target year mth month in a typical day h period.
12) N-1 abandon wind power restraint under the condition of anticipated accidents:
Figure BDA0002612204280000251
wherein,
Figure BDA0002612204280000252
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 BDA0002612204280000253
Figure BDA0002612204280000254
14) N-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure BDA0002612204280000255
Figure BDA0002612204280000256
wherein,
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
means 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,
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 (6)

1. A source-grid-load coordinated power transmission network planning method is characterized by comprising the following operations:
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, load power in 24 hours 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 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 range of a candidate newly-built power transmission line according to the condition of a power transmission line corridor and setting construction cost, and setting a wind power curtailment cost coefficient and a system-allowed curtailment rate;
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;
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 objective function expression of the optimization model is as follows:
Figure FDA0003763659450000011
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 alpha w Representing wind curtailment penalty cost;
Figure FDA0003763659450000012
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 a g Representing the flexibility modification cost of the thermal power generating unit g; ug is a binary variable and indicates whether the thermal power generating unit g is subjected to flexible transformation in a planned target year, and u is a variable g 1 represents that the thermal power generating unit g is subjected to flexibility modification in a planning target year, and u represents g When the thermal power generating unit g is not subjected to flexibility transformation in a planned target year, the thermal power generating unit g is represented as 0; c l Investment cost for construction of a single-circuit power transmission line on a power transmission corridor l; z l Representing the number of newly added transmission lines on the transmission line corridor l;
Figure FDA0003763659450000013
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 FDA0003763659450000021
marginal price for the energizable load c;
Figure FDA0003763659450000022
for planning the mth month of the target yearActive power with interrupted interruptible load r in typical day h period; eta r The marginal price of the interruptible load r for a time period t.
2. The source-grid-load coordinated power transmission network planning method according to claim 1, wherein the constraint condition comprises:
1) the thermal power unit active power upper and lower limit constraint:
Figure FDA0003763659450000023
wherein,
Figure FDA0003763659450000024
and
Figure FDA0003763659450000025
respectively representing the upper limit and the lower limit of the active power of the thermal power generating unit g;
Figure FDA0003763659450000026
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 P g The 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 FDA0003763659450000027
wherein R is g The ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ R g The 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 FDA0003763659450000028
wherein, P l m,h 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 FDA0003763659450000029
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 FDA00037636594500000210
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 FDA00037636594500000211
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 FDA00037636594500000212
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; s. the i And E i The total number of the power transmission lines with the node i as a head node and the tail end node respectively; g i 、W i And D i Respectively representing the total number of thermal power generating units, wind power plants and loads on the node i; c i And R i Respectively 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 FDA0003763659450000031
wherein gamma represents a set wind power curtailment rate threshold value;
5) demand response load response range constraints
Figure FDA0003763659450000032
Figure FDA0003763659450000033
Wherein,
Figure FDA0003763659450000034
and
Figure FDA0003763659450000035
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 FDA0003763659450000036
and
Figure FDA0003763659450000037
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 FDA0003763659450000038
Figure FDA0003763659450000039
wherein, B l The susceptance of a single-circuit power transmission line on a power transmission line corridor l;
Figure FDA00037636594500000310
a node i voltage phase angle is planned for a typical day h period of the mth month of a target year;
Figure FDA00037636594500000311
the number of the existing transmission lines on the transmission line corridor l is shown; z l Representing the number of newly added transmission lines on the transmission line corridor l; p l m,h 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; p is l max 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:
Z l ≤Z l max ,l＝1,…,nl
wherein Z is l max The maximum number of the extensible power transmission lines is set for the power transmission line corridor l;
8) node voltage phase angle constraint:
Figure FDA0003763659450000041
wherein,
Figure FDA0003763659450000042
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 FDA0003763659450000043
wherein,
Figure FDA0003763659450000044
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 FDA0003763659450000045
wherein R is g The ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ R g The method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure FDA0003763659450000046
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 FDA0003763659450000047
active power of the thermal power generating unit g under an expected accident k for planning a typical day h-1 period of the mth month of a target year;
11) n-1 node power balance constraint under the condition of expected accident:
Figure FDA0003763659450000048
wherein, P l m,h,k 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 FDA0003763659450000049
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 FDA00037636594500000410
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 FDA00037636594500000411
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 FDA00037636594500000412
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 FDA0003763659450000051
wherein,
Figure FDA0003763659450000052
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 constraint under N-1 forecast accident situation
Figure FDA0003763659450000053
Figure FDA0003763659450000054
14) N-1 transmission capacity constraint of the transmission line under the condition of anticipated accidents:
Figure FDA0003763659450000055
Figure FDA0003763659450000056
wherein,
Figure FDA0003763659450000057
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure FDA0003763659450000058
indicating a wire outage in the k power line corridor under an expected accident,
Figure FDA0003763659450000059
representing that no line is stopped on the transmission line corridor under the expected accident k;
Figure FDA00037636594500000510
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 FDA00037636594500000511
wherein,
Figure FDA00037636594500000512
and forecasting the voltage phase angle of the node i under the accident k for the planning of the mth month, the typical day h period.
3. The source-network-load-coordinated power transmission network planning method according to claim 2, wherein a nonlinear constraint expression in the transmission capacity constraint of the transmission line is converted into an equivalent linear expression:
Figure FDA00037636594500000513
Figure FDA0003763659450000061
wherein M is a constant.
4. The source-network-load coordinated power transmission network planning method according to claim 2, wherein a nonlinear constraint expression in transmission capacity constraint of the transmission line under the N-1 expected accident condition is converted into an equivalent linear expression:
Figure FDA0003763659450000062
5. a grid-source 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;
the model parameter setting module is used for inputting parameters of a traditional thermal power unit and parameters of a power transmission element of a current power grid before a model is built, 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 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 newly-built power transmission line range and construction cost according to the power transmission line corridor condition, and setting a wind power curtailment cost coefficient and a system-allowed curtailment rate;
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 objective function expression of the optimization model is as follows:
Figure FDA0003763659450000063
in the formula, nw is the total number of the wind power plants; ng is the number of the 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;α w representing wind abandon penalty cost;
Figure FDA0003763659450000071
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 a g Representing the flexibility modification cost of the thermal power generating unit g; u. of g Is a binary variable and indicates whether the thermal power generating unit g is subjected to flexible transformation in a planned target year or not, u g 1 represents that the thermal power generating unit g is subjected to flexibility modification in a planning target year, and u represents g When the thermal power generating unit g is not transformed flexibly in the planning target year, 0 is shown; c l Investment cost for construction of a single-circuit power transmission line on a power transmission corridor l; z l Representing the number of newly added transmission lines on the transmission line corridor l;
Figure FDA0003763659450000072
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 FDA0003763659450000073
marginal price for the energizable load c;
Figure FDA0003763659450000074
the active power of the interruptible load r interrupted in a typical day h period of the mth month of the planning target year; eta r The marginal price of the interruptible load r for a time period t.
6. The source-grid-load coordinated power transmission network planning system according to claim 5, 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 FDA0003763659450000075
wherein,
Figure FDA0003763659450000076
and
Figure FDA0003763659450000077
respectively representing the upper limit and the lower limit of the active power of the thermal power generating unit g;
Figure FDA0003763659450000078
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 P g The method comprises the steps that the depth peak regulation power which can be increased by modifying the thermal power unit g through flexibility is shown;
2) conventional thermal power generating unit climbing restraint:
Figure FDA0003763659450000079
wherein R is g The ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ R g The 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 FDA0003763659450000081
wherein, P l m,h 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 FDA0003763659450000082
planning an active power predicted value of a wind power plant w in a typical day h period of the mth month of a target year;
Figure FDA0003763659450000083
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 FDA0003763659450000084
for planning typical day h period of mth month of target yearAn active power predicted value of the excitation load c;
Figure FDA0003763659450000085
an active power predicted value of an interruptible load r for planning a typical day h period of the mth month of a target year; s i And E i The total number of the power transmission lines with the node i as a head node and the tail end node respectively; g i 、W i And D i Respectively representing the total number of thermal power generating units, wind power plants and loads on the node i; c i And R i Respectively 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) abandon wind power constraint:
Figure FDA0003763659450000086
Figure FDA0003763659450000087
wherein gamma represents a set wind power curtailment rate threshold value;
5) demand response load response range constraints
Figure FDA0003763659450000088
Figure FDA0003763659450000089
Wherein,
Figure FDA00037636594500000810
and
Figure FDA00037636594500000811
method for stimulating load c for planning typical day h period of mth month of target yearUpper and lower power limits;
Figure FDA00037636594500000812
and
Figure FDA00037636594500000813
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 FDA00037636594500000814
Figure FDA00037636594500000815
wherein, B l The susceptance of a single-circuit power transmission line on a power transmission line corridor l;
Figure FDA0003763659450000091
planning a voltage phase angle of a node i in a typical day h period of the mth month of a target year;
Figure FDA0003763659450000092
the number of the existing transmission lines on the transmission line corridor l is shown; z l Representing the number of newly added transmission lines on the transmission line corridor l; p is l m,h 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; p l max 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:
Z l ≤Z l max ,l＝1,…,nl
wherein Z is l max The maximum number of the extensible power transmission lines is 1;
8) node voltage phase angle constraint:
Figure FDA0003763659450000093
wherein,
Figure FDA0003763659450000094
a node i voltage phase angle is planned for a typical day h period of the mth month of a target year;
9) and (2) constraining upper and lower limits of active power of the thermal power generating unit under the condition of an N-1 anticipated accident:
Figure FDA0003763659450000095
wherein,
Figure FDA0003763659450000096
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 FDA0003763659450000097
wherein R is g The ramp rate is the ramp rate of the thermal power generating unit g without modification; Δ R g The method comprises the steps of changing the climbing rate of a thermal power generating unit g after flexibility modification;
Figure FDA0003763659450000098
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 FDA0003763659450000099
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 FDA0003763659450000101
wherein, P l m,h,k 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 FDA0003763659450000102
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 FDA0003763659450000103
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 FDA0003763659450000104
active power increased by an excitable load c under an forecast accident k for a typical day h period of the mth month of a planned target year;
Figure FDA0003763659450000105
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 FDA0003763659450000106
wherein,
Figure FDA0003763659450000107
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 FDA0003763659450000108
Figure FDA0003763659450000109
14) N-1 transmission capacity constraint of the transmission line under the condition of expected accidents:
Figure FDA00037636594500001010
Figure FDA00037636594500001011
wherein,
Figure FDA00037636594500001012
indicating whether the power transmission line corridor is shut down or not under the expected accident k,
Figure FDA00037636594500001013
indicating a wire outage in the k power line corridor under an expected accident,
Figure FDA00037636594500001014
representing that no line is stopped on the transmission line corridor under the expected accident k; p l m,h,k 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 FDA0003763659450000111
wherein,
Figure FDA0003763659450000112
and forecasting the voltage phase angle of the node i under the accident k for the planning of the mth month, the typical day h period.
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 CN111799793A (en) 2020-10-20
CN111799793B true 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)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113722940B (en) * 2021-11-04 2022-02-08 国网湖北省电力有限公司经济技术研究院 Power grid unit combination and technical improvement plan combined optimization method
CN115062817A (en) * 2022-05-09 2022-09-16 国网山东省电力公司电力科学研究院 Source network load cooperative scheduling method and system considering air conditioner load
CN116703084B (en) * 2023-06-05 2024-03-22 东北林业大学 Multi-technology flexibility transformation planning method for coal-fired unit based on high wind power permeability

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109102115B (en) * 2018-08-03 2021-09-14 国网山东省电力公司经济技术研究院 Reference power grid opportunity constraint planning method suitable for large-scale wind power 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
CN109713737B (en) * 2018-09-30 2021-12-03 中国电力科学研究院有限公司 New energy delivery capacity evaluation method and system for accessing flexible direct current power grid
CN109858774B (en) * 2019-01-09 2020-10-27 燕山大学 Source network load planning method for improving system safety and coordination
CN110852565B (en) * 2019-10-10 2023-05-12 国家电网有限公司 Grid frame planning method considering different functional attributes

Also Published As

Publication number Publication date
CN111799793A (en) 2020-10-20

Similar Documents

Publication Publication Date Title
CN111799793B (en) Source-grid-load cooperative power transmission network planning method and system
Chen et al. Key technologies for integration of multitype renewable energy sources—Research on multi-timeframe robust scheduling/dispatch
Franco et al. A network flow model for short-term hydro-dominated hydrothermal scheduling problems
CN109767078B (en) Multi-type power supply maintenance arrangement method based on mixed integer programming
Kuzle et al. Virtual power plant dispatch optimization using linear programming
Mc Garrigle et al. Quantifying the value of improved wind energy forecasts in a pool-based electricity market
Soroudi et al. Stochastic real-time scheduling of wind-thermal generation units in an electric utility
Xie et al. Optimal capacity and type planning of generating units in a bundled wind–thermal generation system
Aliari et al. Planning for integration of wind power capacity in power generation using stochastic optimization
CN104809545B (en) A kind of virtual plant runs modeling method
CN113988714B (en) Multi-uncertainty-based dynamic planning method, equipment and medium for park comprehensive energy system
CN105243600A (en) Grid power generation adjustment method
CN108767855B (en) Random production simulation method for time sequence continuous hybrid power system
CN103633641B (en) A kind ofly consider the medium and long-term transaction operation plan acquisition methods that wind-powered electricity generation is received
Li et al. Modeling and optimal energy operation considering probabilistic and sustainable renewable energy models and demand side management
CN118114469A (en) Virtual power plant participation power distribution network collaborative optimization scheduling method and system
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
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
CN111799842B (en) Multi-stage power transmission network planning method and system considering flexibility of thermal power generating unit
Rendroyoko et al. Integration method of unit commitment using PL-GA binary dispatch algorithm for intermittent RES in isolated microgrids system
Mohammadhosseini et al. A combination of genetic algorithm and particle swarm optimization for power systems planning subject to energy storage
Yang et al. Microgrid's generation expansion planning considering lower carbon economy
Kamalanathan et al. A Novel Approach Based Optimal Power Scheduling of GENCOs to Improve the Profit in Electricity Market Considering Wind Power 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