CN112260268B - Method for obtaining system flexibility domain - Google Patents

Method for obtaining system flexibility domain Download PDF

Info

Publication number
CN112260268B
CN112260268B CN202011003447.6A CN202011003447A CN112260268B CN 112260268 B CN112260268 B CN 112260268B CN 202011003447 A CN202011003447 A CN 202011003447A CN 112260268 B CN112260268 B CN 112260268B
Authority
CN
China
Prior art keywords
domain
wind field
period
max
flexibility
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011003447.6A
Other languages
Chinese (zh)
Other versions
CN112260268A (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.)
Economic and Technological Research Institute of State Grid Hebei Electric Power Co Ltd
Original Assignee
Economic and Technological Research Institute of State Grid Hebei 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 Economic and Technological Research Institute of State Grid Hebei Electric Power Co Ltd filed Critical Economic and Technological Research Institute of State Grid Hebei Electric Power Co Ltd
Priority to CN202011003447.6A priority Critical patent/CN112260268B/en
Publication of CN112260268A publication Critical patent/CN112260268A/en
Application granted granted Critical
Publication of CN112260268B publication Critical patent/CN112260268B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas 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/004—Generation forecast, e.g. methods or systems for forecasting future energy generation
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

The invention discloses a method for acquiring a system flexibility domain, which comprises the following steps of 1) acquiring a basic operation point; 2) A flexibility domain solving model; 3) Optimizing the flexible domain solving model to obtain a flexible domain simple solving model; 4) Flexibility domain acquisition. According to the invention, the workload and the working time are greatly reduced by improving the solving mode of the step (2), and the time for acquiring the system flexibility domain is greatly shortened.

Description

Method for obtaining system flexibility domain
Technical Field
The invention relates to the field of power system analysis, in particular to a method for acquiring a system flexibility domain.
Background
Under the consensus of the global main countries on the climate change, the rapid development of renewable energy power generation technology represented by wind and light is promoted. However, inherent fluctuation and uncertainty of wind and light power generation output cause the problem of large-scale wind and light power generation digestion to be increasingly remarkable, and the risk of safe and stable operation of the system is continuously increased. The concept of the flexibility domain of the system is given birth, which refers to the upper and lower boundaries of renewable energy sources which can be consumed by the system under the condition of the existing basic operation point, namely the maximum fluctuation interval of renewable energy source output which can be tolerated by the system.
The size of the flexibility domain of the system is affected by flexibility resources, such as minimum technical output and climbing speed of a traditional conventional unit. At present, the system flexibility domain is studied in a plurality of ways, but is mostly calculated by constructing a scheduling model, but the problem is a two-stage robust planning problem, which is generally obtained by solving through Benders decomposition or CC & G algorithm, and has high calculation complexity and long time.
Disclosure of Invention
The present invention is directed to a method for obtaining a system flexibility domain, so as to solve the problems set forth in the background art.
In order to achieve the above purpose, the present invention provides the following technical solutions: a method for obtaining a system flexibility domain, comprising the steps of:
1) Obtaining a basic operation point;
its objective function:
wherein P is g,t b The basic operation point of the thermal power unit g in the period t is set; p (P) k,t b The basic operation point of the wind field k in the period t is set; p (P) k,t,pre Predicted output for wind field k during period t; c (C) g (. Cndot.) is the thermal power unit g cost function; c (C) w Punishment coefficients for the wind curtailment; g is a thermal power unit set; k is a wind field set; t is the set of time periods.
Wherein D is t For system period tLoad demand.
In SF n,l A power transfer distribution factor of the node n to the line l; f (F) l,max Is the maximum capacity of line l; d (D) n,t Is the load demand of node n during period t.
P g,min ≤P g,t b ≤P g,max
Wherein P is g,min 、P g,max The upper and lower limits of the output of the thermal power unit g are respectively set.
0≤P k,t b ≤P k,t,pre
2) A flexibility domain solving model;
its objective function:
in sigma k The wind field k flexibility domain weight is given; p (P) k,t,UB For wind field k to absorb upper bound, P k,t,LB The lower bound is absorbed for the wind field k, and the lower bound and the wind field k are the flexibility domain; p (P) k,rate Is the installed capacity of the wind farm k.
Wherein P is g,t * The method comprises the steps of (1) setting an operating point of a thermal power unit g period t in a flexibility domain solution; p (P) k,t * And (3) operating points in the flexibility domain solving for the wind field k period t.
P g,t b -Δ g,dn ≤P g,t * ≤P g,t b +Δ g,up ，P g,min ≤P g,t * ≤P g,max
In the formula delta g,up 、Δ g,dn The up-and-down climbing capacity of the unit g.
0≤P k,LB ≤P k,t b ≤P k,UB ≤P k,rate
3) Optimizing the flexible domain solving model to obtain a flexible domain simple solving model;
its upper bound objective function:
P g,t b -Δ g,dn ≤P g,t * ≤P g,t b +Δ g,up ，P g,min ≤P g,t * ≤P g,max
P k,t b ≤P k,t * ≤P k,rate
its lower bound objective function:
P g,t b -Δ g,dn ≤P g,t * ≤P g,t b +Δ g,up ，P g,min ≤P g,t * ≤P g,max
P k,t b ≤P k,t * ≤P k,rate
4) Flexible domain acquisition
Wherein T is N Is the number of time periods.
Compared with the prior art, the invention has the beneficial effects that: the flexibility domain solving model in the step (2) is a two-stage robust planning model, the calculating complexity is high, the solving is complex, the flexibility domain solving model is generally obtained by solving through a Benders decomposition or CC & G algorithm, the calculating complexity is high, the cost is long, the working efficiency is low, and the workload and the working time are greatly reduced by improving the solving mode in the step (2).
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a system flexibility domain of an embodiment;
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Taking the IEEE-RTS79 as an example of a power generation reliability test system for illustration, adding wind fields to the node 1 and the node 2, wherein the wind abandoning penalty is 0.3 per kWh, and other parameters are referred to the IEEE-RTS79 test system.
1) Basic operating point acquisition
Through the foundationOperating Point model acquisition P g,t b 、P k,t b 。
2) Flexible domain upper bound acquisition
Inputting basic operation point and up-down climbing parameter delta of unit g,up 、Δ g,dn Obtaining through a flexible domain upper bound model
3) Flexibility domain lower bound acquisition
Inputting basic operation point and up-down climbing parameter delta of unit g,up 、Δ g,dn Obtaining through a flexible domain lower bound model
4) Obtain system flexibility domain [ P ] LB ；P UB ]The results are shown in FIG. 2.
Thus, the method provided by the invention is implemented.
It is worth mentioning that the invention can calculate not only the flexibility domain of renewable energy sources, but also the flexibility domain of other uncertainty resources, such as electric loads.
It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential characteristics thereof. The present embodiments are, therefore, to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference sign in a claim should not be construed as limiting the claim concerned.

Claims (1)

1. A method for obtaining a system flexibility domain, comprising the steps of:
1) Obtaining a basic operation point;
its objective function:
wherein P is g,t b The basic operation point of the thermal power unit g in the period t is set; p (P) k,t b The basic operation point of the wind field k in the period t is set; p (P) k,t,pre Predicted output for wind field k during period t; c (C) g (. Cndot.) is the thermal power unit g cost function; c (C) w Punishment coefficients for the wind curtailment; g is a thermal power unit set; k is a wind field set; t is a time period set;
wherein D is t Load demand for system period t;
in SF n,l A power transfer distribution factor of the node n to the line l; f (F) l,max Is the maximum capacity of line l; d (D) n,t The load requirement of the node n in the period t is met;
P g,min ≤P g,t b ≤P g,max
wherein P is g,min 、P g,max The upper and lower output limits of the thermal power unit g are respectively set;
0≤P k,t b ≤P k,t,pre
2) A flexibility domain solving model;
its objective function:
in sigma k The wind field k flexibility domain weight is given; p (P) k,t,UB For wind field k to absorb upper bound, P k,t,LB The lower bound is absorbed for the wind field k, and the lower bound and the wind field k are the flexibility domain; p (P) k,rate Is the installed capacity of the wind field k;
wherein P is g,t * The method comprises the steps of (1) setting an operating point of a thermal power unit g period t in a flexibility domain solution; p (P) k,t * The method comprises the steps of (1) setting an operating point of a wind field k period t in a flexibility domain solution;
P g,t b -Δ g,dn ≤P g,t * ≤P g,t b +Δ g,up ，P g,min ≤P g,t * ≤P g,max
in the formula delta g,up 、Δ g,dn The up-down climbing capacity of the unit g;
0≤P k,LB ≤P k,t b ≤P k,UB ≤P k,rate
3) Optimizing the flexible domain solving model to obtain a flexible domain simple solving model;
its upper bound objective function:
P g,t b -Δ g,dn ≤P g,t * ≤P g,t b +Δ g,up ，P g,min ≤P g,t * ≤P g,max
P k,t b ≤P k,t * ≤P k,rate
its lower bound objective function:
P g,t b -Δ g,dn ≤P g,t * ≤P g,t b +Δ g,up ，P g,min ≤P g,t * ≤P g,max
P k,t b ≤P k,t * ≤P k,rate
4) Flexible domain acquisition
[P LB ；P UB ]，Wherein T is N Is the number of time periods.
CN202011003447.6A 2020-09-22 2020-09-22 Method for obtaining system flexibility domain Active CN112260268B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011003447.6A CN112260268B (en) 2020-09-22 2020-09-22 Method for obtaining system flexibility domain

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011003447.6A CN112260268B (en) 2020-09-22 2020-09-22 Method for obtaining system flexibility domain

Publications (2)

Publication Number Publication Date
CN112260268A CN112260268A (en) 2021-01-22
CN112260268B true CN112260268B (en) 2023-08-29

Family

ID=74232597

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011003447.6A Active CN112260268B (en) 2020-09-22 2020-09-22 Method for obtaining system flexibility domain

Country Status (1)

Country Link
CN (1) CN112260268B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109687531A (en) * 2019-01-18 2019-04-26 华北电力大学 A kind of fired power generating unit flexibility remodeling method under large-scale wind power access
CN110247426A (en) * 2019-06-12 2019-09-17 国网山西省电力公司电力科学研究院 A kind of robust Unit Combination method based on the uncertain set of multiband
CN110797919A (en) * 2019-12-05 2020-02-14 国网四川省电力公司经济技术研究院 Clean energy power supply planning method based on Wasserstein distance and distribution robust optimization
CN110991857A (en) * 2019-11-28 2020-04-10 华北电力大学 Method for evaluating wind power consumption capability of electric heating integrated energy system
CN111341424A (en) * 2020-02-26 2020-06-26 杭州电子科技大学 Operation scheduling optimization method and system based on two-stage robust optimization model

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109687531A (en) * 2019-01-18 2019-04-26 华北电力大学 A kind of fired power generating unit flexibility remodeling method under large-scale wind power access
CN110247426A (en) * 2019-06-12 2019-09-17 国网山西省电力公司电力科学研究院 A kind of robust Unit Combination method based on the uncertain set of multiband
CN110991857A (en) * 2019-11-28 2020-04-10 华北电力大学 Method for evaluating wind power consumption capability of electric heating integrated energy system
CN110797919A (en) * 2019-12-05 2020-02-14 国网四川省电力公司经济技术研究院 Clean energy power supply planning method based on Wasserstein distance and distribution robust optimization
CN111341424A (en) * 2020-02-26 2020-06-26 杭州电子科技大学 Operation scheduling optimization method and system based on two-stage robust optimization model

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
考虑灵活弃风的电力系统优化调度;王卓群;刘蕾;马平;;青岛农业大学学报(自然科学版)(第01期);全文 *

Also Published As

Publication number Publication date
CN112260268A (en) 2021-01-22

Similar Documents

Publication Publication Date Title
CN105811407B (en) A kind of micro-capacitance sensor primary frequency modulation control method based on distributed Newton method
CN103269091B (en) Wind-solar energy storage capacity configuring method based on wind-solar energy average output curve
CN108695857B (en) Automatic voltage control method, device and system for wind power plant
CN107276122B (en) Peak-shaving resource calling decision method suitable for large-scale renewable energy grid connection
CN112994115B (en) New energy capacity configuration method based on WGAN scene simulation and time sequence production simulation
Yao et al. A control configuration of wind farm for load-following and frequency support by considering the inertia issue
CN108960642B (en) New energy power plant dynamic polymerization method and system
Li et al. Coordinated control scheme of a hybrid renewable power system based on hydrogen energy storage
CN116402210A (en) Multi-objective optimization method, system, equipment and medium for comprehensive energy system
CN110601184A (en) Wind power plant extended photovoltaic multi-objective optimization method considering boost main transformer capacity
CN112583051B (en) Optimized scheduling model construction method of variable-speed pumped storage unit in regional power grid
CN112260268B (en) Method for obtaining system flexibility domain
CN106786766B (en) A method of the raising wind-powered electricity generation maximum grid connection capacity based on P2G technology
CN103972923A (en) Multi-combination device of large solar photovoltaic grid-connected power generation system
CN103473618A (en) System and method for predicting short-term wind power of wind power plant
Kazemi et al. A novel method for estimating wind turbines power output based on least square approximation
CN114400660A (en) Method and system for determining proportion of wind power photovoltaic hydrogen production system
Ma et al. Multi-Point Layout Planning of Multi-Energy Power Supplies Based on Time-series Production Simulation
CN106026197A (en) Wind power and thermal power coupling delivery optimization control method
CN108919878A (en) A kind of cooperation control system of power in power plant
CN115563815B (en) Method and device for simulating and calculating space-time evolution of carbon emission flow of double-high power system
CN113644669B (en) Energy storage capacity configuration method and system based on energy storage capacity utilization rate
CN110601263B (en) Wind power plant access point voltage risk assessment method based on node type transformation method
CN109245133B (en) Energy storage regulation and control method for microgrid, electronic equipment and computer readable storage medium
CN117674112A (en) Multi-target electro-hydrogen fusion energy hub method and device for participating in power grid peak shaving

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