CN114243677A - New energy power system coordinated operation simulation method - Google Patents

New energy power system coordinated operation simulation method Download PDF

Info

Publication number
CN114243677A
CN114243677A CN202111294116.7A CN202111294116A CN114243677A CN 114243677 A CN114243677 A CN 114243677A CN 202111294116 A CN202111294116 A CN 202111294116A CN 114243677 A CN114243677 A CN 114243677A
Authority
CN
China
Prior art keywords
new energy
output
unit
power
constraint
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111294116.7A
Other languages
Chinese (zh)
Inventor
邓娇娇
王鹏
郑晓明
李旭霞
王尧
李佳
刘红丽
胡迎迎
王凯凯
荆永明
陈洁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Economic and Technological Research Institute of State Grid Shanxi Electric Power Co Ltd
Original Assignee
Economic and Technological Research Institute of State Grid Shanxi 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 Shanxi Electric Power Co Ltd filed Critical Economic and Technological Research Institute of State Grid Shanxi Electric Power Co Ltd
Priority to CN202111294116.7A priority Critical patent/CN114243677A/en
Publication of CN114243677A publication Critical patent/CN114243677A/en
Pending legal-status Critical Current

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
    • 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
    • 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
    • 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/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • 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]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/50—Photovoltaic [PV] energy
    • Y02E10/56—Power conversion systems, e.g. maximum power point trackers
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/76—Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a coordination operation simulation method of a new energy power system, which comprises a power system deterministic operation simulation model based on a time sequence load curve, wherein the core of the model is a scheduling simulation model taking daily operation as the core. The invention has the beneficial effects that: the method is innovated on the basis of the traditional operation simulation model, is oriented to the access problem of large-scale new energy, and establishes the new energy power system coordinated operation simulation method which fully reflects the intermittence and the randomness of the new energy output. The method lays theoretical and technical foundation for improving the technical level of power system planning and promoting the optimal operation of the power system and the medium-long-term sustainable development of energy in the future.

Description

New energy power system coordinated operation simulation method
Technical Field
The invention relates to a power system operation simulation method, in particular to a new energy power system coordinated operation simulation method, and belongs to the technical field of power system operation.
Background
With the continuous development of clean and green energy, the installed capacity of renewable energy sources such as wind power and photovoltaic is continuously improved.
The new energy may have explosive growth in the future, impact and challenge which cannot be ignored are caused to the reliable and stable operation of the power system, the complexity of the operation of the power system is increased, and the multi-link coordinated operation problem is increasingly prominent. It can be expected that the high-proportion situation will be gradually formed by the power generation of the wind power, the photovoltaic power and other energy sources with uncertain output in most areas in China, and in order to meet the challenge, the flexible resources in the new energy sources will be fully mobilized and utilized, so that the electric power system in most areas in the future will have the characteristics of diversified main elements, complicated operation mode and high-dimensional planning scheme, and research and development of the key technology for coordinating and planning the new energy electric power system are urgently needed.
In the aspect of research on new energy consumption problems, currently, research works such as consumption of new energy power generation in China are relatively limited, and research in this aspect at home and abroad mainly focuses on grid-connected management and power generation technologies, and documents [ Pei philosophy, dingjie, lie morning, zhongchung, liangxing, bright in, xu spring, zhang. The document [ cuhuan, discussion about a new energy power generation wind power generation technology, science and technology wind 2020,135-36(2020) ], researches and analyzes key technologies of a wind turbine generator and a wind power plant, and proposes a suggestion for enhancing innovation and application of the wind power generation technology;
in the aspect of peak shaving research of an electric power system, the current peak shaving modes in China mainly focus on hydroelectric and thermal power units, and documents [ Zhang Chang, Yang Jianhua, Shuaihang, Shukang and ai Xiao Meng ] market auxiliary service research based on cross-regional electric power trading of a power grid in China, Chinese power 50,139-45(2017) ] propose a compensation calculation method for peak shaving of the thermal power unit in order to improve peak shaving willingness of the thermal power unit, documents [ Tangpolitics, Shijun, yellow bright, Marangwen, Liu Yan ] propose a peak shaving load distribution research based on a moth flame optimization algorithm, and results of further reducing the total peak shaving cost of the system are obtained by introducing MFO as a model solution method. Thermal power generation is a main power supply form in China, but the adjusting capacity of the thermal power generation is limited by the operation of a unit, so that how to play the adjusting role of high-proportion thermal power in China still needs to be further researched;
in the aspect of power system production simulation research, the main production simulation methods at present comprise power and electric quantity balance, random production simulation, Monte Carlo simulation and the like. Although a lot of researches are carried out at home and abroad aiming at the problems of unit combination, economic dispatch and the like, documents [ Zhu Zeli, Zhou Jing Yang, Pan Yi, Yan Cui Hui, Cui Hui, Liwei rigid, safety constraint economic dispatch considering power and electric quantity balance, China Motor engineering reports 33,168-76(2013) ] and documents [ Zhang Ning, Chen Hui Kun, Luo Xiao Ming, Lijialong, Xiuqing, Kangqing, Guangdong power grid energy-saving power generation dispatching model and algorithm, power grid technologies 32,11-15(2008) ] provide various methods for optimal balance of power and electric quantity, but the methods cannot consider the influence of load characteristics and power supply composition. The calculation method of random production simulation has been widely used due to high calculation efficiency and flexible problem handling, and the literature [ Kangchongqing, Baili super, Xianhui, Xiangniande, implementation of random production simulation algorithm based on sequence operation theory, Chinese electro-mechanical engineering newspaper 2002,7-12(2002) ] proposes to calculate the reliability and economic index of the system in the supply and demand balance process by using the sequence operation theory, and obtains the simulation result with time sequence, but
The error of this method is difficult to control. The Monte Carlo method is relatively simple, and documents [ R.Billington, HuaChen, R.Ghajar.A. sampling simulation for frequency generation system including simulation of power system production of renewable power source are proposed in IEEE11, 728-34 (1996) ].
Disclosure of Invention
The invention aims to provide a method for simulating coordinated operation of a new energy power system to solve the problems.
The invention realizes the purpose through the following technical scheme: a coordination operation simulation method for a new energy power system comprises the following steps:
step one, according to an installation schedule, considering unit commissioning, decommissioning, technical improvement and the like, and determining a commissioning unit; next, removing the maintenance unit according to the maintenance plan, and finally determining the operable unit and parameters thereof;
step two, arranging all thermal power generating units capable of determining output, including external protocol power transmission, the thermal power generating units and units considering specified output, and correcting corresponding load curves according to the area where the power supply is located;
step three, according to the new energy simulation output of the new energy unit generated by the new energy operation simulation module in a random simulation mode, arranging the new energy output and correcting a corresponding load curve;
fourthly, on the basis of the corrected multi-region load curve, arranging peak clipping and valley filling for the pumped storage of the hydroelectric and pumped storage units and the conventional hydroelectric generating sets, setting the pumped storage mode to be a flat pumping or full pumping mode, meeting the constraints of capacity, electric quantity and the like of the units, and correcting the corresponding load curve again according to the region where the power supply is located;
and step five, performing optimization simulation operation on the rest units, wherein the following parameters are required to be prepared besides the calculation result: the system comprises manually specified unit states, positive and negative spare amount of each time interval, partition spare amount, cross-partition spare amount, time-sharing quotation or cost of the unit, start-stop cost of starting and stopping the unit, network constraint and the like.
As a still further scheme of the invention: the model adopts a power system deterministic operation simulation model based on a time sequence load curve, the core of the model is a scheduling simulation model taking daily operation as the core, and the objective function of the model is expressed as:
Figure BDA0003336001350000041
in the formula: t is the total time interval in the optimization cycle; c (P)t) For each type of unit, the output power is P in t periodtThe cost of operation; subscripts c, f, h, p and w respectively represent that thermal power cannot be started or stopped in a day, thermal power can be started or stopped in a day, water and electricity can be started or stopped in a day, the pumping storage is carried out, and new energy is obtained; cwIn order to cut off the cost of new energy; pwd tSwitching new energy power for t time period; dd tCutting load power for t time period; vdLoad loss is cut for each node; cf、CcThe unit is charged for starting and stopping; θ, η, and γ are weighting coefficients, and are usually 1, and may be adjusted as necessary.
As a still further scheme of the invention: in the second step, the output constraint for the output unit specifically includes:
the thermal power unit output upper and lower limits and the thermal power unit climbing restraint.
Figure BDA0003336001350000051
In the formula: pcmin,Pcmax,Pfmin,PfmaxRespectively the minimum output and the maximum output of the unit; ic is the state variable of the unit which can not be started or stopped in the day,Itfthe state variable is a state variable of a unit which can be started and stopped at a time t within a day;
Figure BDA0003336001350000052
Figure BDA0003336001350000053
the down-climbing speed and the up-climbing speed of the unit are respectively.
As a still further scheme of the invention: in the third step, the output constraint of the new energy unit specifically comprises:
new energy output prediction variable introduced into new energy modeling
Figure BDA0003336001350000054
And a new energy cutting mechanism is introduced into the daily operation simulation model, so that the model can cut off part of the new energy output under the conditions that the system cannot provide peak regulation capacity, the system spare capacity is insufficient or the new energy is blocked to be sent out.
Figure BDA0003336001350000055
In the formula:
Figure BDA0003336001350000061
new energy output for t time period;
Figure BDA0003336001350000062
switching new energy power for t time period;
Figure BDA0003336001350000063
and predicting a force value for the new energy in the t period.
Constraint of power abandonment rate of new energy
Figure BDA0003336001350000064
Wherein the content of the first and second substances,rwdthe upper limit proportion of the power abandoning rate of the new energy is determined,
Figure BDA0003336001350000065
in order to switch the power sum of new energy,
Figure BDA0003336001350000066
the output of new energy is integrated.
As a still further scheme of the invention: in the fourth step, the restriction of the water, electricity and pumped storage output specifically comprises:
the hydroelectric generating set optimizes and distributes the output of each time period in the model according to the output range and daily generated energy of the hydroelectric generating set given by the intermediate-term and long-term cross-basin cascade hydropower optimization scheduling result; the pumping and storage unit considers that the daily pumping amount is balanced with the generated energy.
Figure BDA0003336001350000067
In the formula: phmin、PhmaxThe first line represents the upper and lower limit constraints of the output of the hydroelectric generating set; qhydroThe daily generated energy of the hydroelectric generating set is represented by the second row, the daily generated energy constraint of the hydroelectric generating set is represented by the second row, if the equation of the optimization result is not satisfied, the hydropower station generates water abandon, and the water abandon is converted into electric quantity
Figure BDA0003336001350000068
Pp，pump、Pp，genThe third line represents the upper and lower limit constraints of the output of the pumping and storage unit for the maximum pumping amount and the generated energy of the pumping and storage unit in a unit time period; i ist ppump，、It pgen，In order to describe the state variable of pumping or generating of the pumping and storage unit in the period t, the fourth line represents the mutual exclusion constraint of the pumping and generating states of the pumping and storage unit; lambda [ alpha ]pFor the efficiency of the pumping and storage unit, the fifth element represents the balance constraint of the pumping amount and the generating capacity of the pumping and storage unit each day.
As a still further scheme of the invention: in the fifth step, the positive and negative reserve volume constraint in each time interval specifically includes:
Figure BDA0003336001350000071
in the formula: r isuFor the positive standby rate of the system required in time t, in the formula, the contribution of new energy to the system standby should be predicted to output P according to the predicted outputwf tCalculating that the cut-off part of the spare capacity is counted even if the spare capacity is cut off.
Figure BDA0003336001350000072
In the formula: r isd tThe minimum output of the new energy is considered to be 0 equivalently, and the minimum output of the new energy can be cut off at any time.
As a still further scheme of the invention: in the fifth step, the partition spare amount constraint specifically includes:
the different regions need to ensure that the startup in each region reserves sufficient spare capacity for the load in the region, and the partition spare constraint can be expressed as:
Figure BDA0003336001350000073
in the formula: z is the total number of the areas; dzt,Is the load of the z region at time t; r isu zt,、rd zt，Respectively the positive standby rate and the negative standby rate of the z area at the time period t; z+And Z-The contact line sets sent in and out for the z area respectively; f. ofl tIs the flow of the l lines. The first row of constraints represents the partition positive spare constraint and the second row represents the partition negative spare constraint.
As a still further scheme of the invention: in the fifth step, the network constraint specifically includes:
establishing a line and section transmission constraint based on the direct current power flow model, as shown in the following formula:
Figure BDA0003336001350000081
in the formula: fl t，Fs tRespectively a t-period circuit and a section flow matrix; w is a generator transfer distribution factor; a. thengc，Angf，Angh，Angp，AngwRespectively node-unit incidence matrixes of different types of units; a. theslIs a section-line correlation matrix.
As a still further scheme of the invention: in the fifth step, the cross-region power transmission constraint specifically includes:
in order to simulate a specific power transmission plan in actual power grid operation, cross-region power transmission constraint is introduced, namely, the transmission power flow of each time section of a section is limited.
Figure BDA0003336001350000082
Wherein Θ issThe line set contained in the cross-region power transmission section s is provided.
As a still further scheme of the invention: and step five, dynamic constraint is further included, wherein the thermal power starting and stopping is carried out at most once in the running state in one day, and the large-capacity thermal power generating unit is not allowed to be started or stopped.
The invention has the beneficial effects that: the method is innovated on the basis of the traditional operation simulation model, is oriented to the access problem of large-scale new energy, and establishes the new energy power system coordinated operation simulation method which fully reflects the intermittence and the randomness of the new energy output. The method lays theoretical and technical foundation for improving the technical level of electric power system planning and promoting the optimized operation of the power system and the medium-long-term sustainable development of energy.
Drawings
FIG. 1 is a schematic modeling flow diagram of the present invention;
FIG. 2 is a schematic diagram of installed capacity of various units according to the second embodiment of the present invention;
FIG. 3 is a diagram illustrating the operation of each type of unit in 8760 hours in a second power system according to an embodiment of the present invention;
fig. 4 is a schematic diagram of peak-reduction margin and wind-abandoning and light-abandoning in the second embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Example one
Referring to fig. 1, a method for simulating coordinated operation of a new energy power system includes the following steps:
step one, according to an installation schedule, considering unit commissioning, decommissioning, technical improvement and the like, and determining a commissioning unit; next, removing the maintenance unit according to the maintenance plan, and finally determining the operable unit and parameters thereof;
step two, arranging all thermal power generating units capable of determining output, including external protocol power transmission, the thermal power generating units and units considering specified output, and correcting corresponding load curves according to the area where the power supply is located;
step three, according to the new energy simulation output of the new energy unit generated by the new energy operation simulation module in a random simulation mode, arranging the new energy output and correcting a corresponding load curve;
fourthly, on the basis of the corrected multi-region load curve, arranging peak clipping and valley filling for the pumped storage of the hydroelectric and pumped storage units and the conventional hydroelectric generating sets, setting the pumped storage mode to be a flat pumping or full pumping mode, meeting the constraints of capacity, electric quantity and the like of the units, and correcting the corresponding load curve again according to the region where the power supply is located;
and step five, performing optimization simulation operation on the rest units, wherein the following parameters are required to be prepared besides the calculation result: the system comprises manually specified unit states, positive and negative spare amount of each time interval, partition spare amount, cross-partition spare amount, time-sharing quotation or cost of the unit, start-stop cost of starting and stopping the unit, network constraint and the like.
In the embodiment of the invention, the model adopts a power system deterministic operation simulation model based on a time sequence load curve, the core of the model is a scheduling simulation model taking daily operation as the core, and the objective function of the model is expressed as:
Figure BDA0003336001350000101
in the formula: t is the total time interval in the optimization cycle; c (P)t) For each type of unit, the output power is P in t periodtThe cost of operation; subscripts c, f, h, p and w respectively represent that thermal power cannot be started or stopped in a day, thermal power can be started or stopped in a day, water and electricity can be started or stopped in a day, the pumping storage is carried out, and new energy is obtained; cwIn order to cut off the cost of new energy; pwd tSwitching new energy power for t time period; dd tCutting load power for t time period; vdLoad loss is cut for each node; cf、CcThe unit is charged for starting and stopping; θ, η, and γ are weighting coefficients, and are usually 1, and may be adjusted as necessary.
In the embodiment of the present invention, in the second step, the output constraint on the output unit specifically includes:
the thermal power unit output upper and lower limits and the thermal power unit climbing restraint.
Figure BDA0003336001350000111
In the formula: pcmin，Pcmax，Pfmin，PfmaxRespectively the minimum output and the maximum output of the unit; ic is the state variable of the unit which can not be started or stopped in the day,Itfthe state variable is a state variable of a unit which can be started and stopped at a time t within a day;
Figure BDA0003336001350000112
Figure BDA0003336001350000113
the down-climbing speed and the up-climbing speed of the unit are respectively.
In the embodiment of the present invention, in the third step, the new energy unit output constraint specifically includes:
new energy output prediction variable introduced into new energy modeling
Figure BDA0003336001350000114
And a new energy cutting mechanism is introduced into the daily operation simulation model, so that the model can cut off part of the new energy output under the conditions that the system cannot provide peak regulation capacity, the system spare capacity is insufficient or the new energy is blocked to be sent out.
Figure BDA0003336001350000115
In the formula:
Figure BDA0003336001350000116
new energy output for t time period;
Figure BDA0003336001350000117
switching new energy power for t time period;
Figure BDA0003336001350000118
and predicting a force value for the new energy in the t period.
Constraint of power abandonment rate of new energy
Figure BDA0003336001350000119
Wherein r iswdThe upper limit proportion of the power abandoning rate of the new energy is determined,
Figure BDA00033360013500001110
in order to switch the power sum of new energy,
Figure BDA00033360013500001111
for new energy exportAnd (4) force integration.
In the fourth step of the present invention, the constraints of the hydropower and the pumped storage output specifically include:
the hydroelectric generating set optimizes and distributes the output of each time period in the model according to the output range and daily generated energy of the hydroelectric generating set given by the intermediate-term and long-term cross-basin cascade hydropower optimization scheduling result; the pumping and storage unit considers that the daily pumping amount is balanced with the generated energy.
Figure BDA0003336001350000121
In the formula: phmin、PhmaxThe first line represents the upper and lower limit constraints of the output of the hydroelectric generating set; qhydroThe daily generated energy of the hydroelectric generating set is represented by the second row, the daily generated energy constraint of the hydroelectric generating set is represented by the second row, if the equation of the optimization result is not satisfied, the hydropower station generates water abandon, and the water abandon is converted into electric quantity
Figure BDA0003336001350000122
Pp,pump、Pp,genThe third line represents the upper and lower limit constraints of the output of the pumping and storage unit for the maximum pumping amount and the generated energy of the pumping and storage unit in a unit time period; i ist ppump,、It pgen,In order to describe the state variable of pumping or generating of the pumping and storage unit in the period t, the fourth line represents the mutual exclusion constraint of the pumping and generating states of the pumping and storage unit; lambda [ alpha ]pFor the efficiency of the pumping and storage unit, the fifth element represents the balance constraint of the pumping amount and the generating capacity of the pumping and storage unit each day.
In the embodiment of the present invention, in the fifth step, the positive/negative spare amount constraint in each time interval specifically includes:
Figure BDA0003336001350000123
in the formula: r isuFor the positive standby rate of the system required in time t, in the formula, the contribution of new energy to the system standby should be predicted to output P according to the predicted outputwf tIs calculated, i.e.So that the cut-out portion thereof also takes into account the spare capacity.
Figure BDA0003336001350000124
In the formula: r isd tThe minimum output of the new energy is considered to be 0 equivalently, and the minimum output of the new energy can be cut off at any time.
In the fifth step, the partition spare amount constraint specifically includes:
the different regions need to ensure that the startup in each region reserves sufficient spare capacity for the load in the region, and the partition spare constraint can be expressed as:
Figure BDA0003336001350000131
in the formula: z is the total number of the areas; dzt,Is the load of the z region at time t; r isu zt,、rd zt,Respectively the positive standby rate and the negative standby rate of the z area at the time period t; z+And Z-The contact line sets sent in and out for the z area respectively; f. ofl tIs the flow of the l lines. The first row of constraints represents the partition positive spare constraint and the second row represents the partition negative spare constraint.
In the fifth step, the network constraint specifically includes:
establishing a line and section transmission constraint based on the direct current power flow model, as shown in the following formula:
Figure BDA0003336001350000134
in the formula: fl t，Fs tRespectively a t-period circuit and a section flow matrix; w is a generator transfer distribution factor; a. thengc，Angf，Angh，Angp，AngwRespectively node-unit incidence matrixes of different types of units; a. theslIs a section-line correlation matrix.
In the embodiment of the present invention, in the fifth step, the cross-region power transmission constraint specifically includes:
in order to simulate a specific power transmission plan in actual power grid operation, cross-region power transmission constraint is introduced, namely, the transmission power flow of each time section of a section is limited.
Figure BDA0003336001350000141
Wherein Θ issThe line set contained in the cross-region power transmission section s is provided.
In the embodiment of the invention, in the fifth step, dynamic constraint is further included, wherein the starting and stopping of the thermal power generating unit are performed at most once in the operating state in one day, and the large-capacity thermal power generating unit is not allowed to be started and stopped.
Example two
Referring to fig. 2 to 4, a new energy power system coordinated operation simulation method adopts data of 220kV and 500kV main racks of a certain power saving network, including 220kV substation load data, ac/dc outgoing data, industrial heavy load data of 8760 hours a year round in 2019, and output data of a wind power plant and a photovoltaic power plant mounted on the main racks.
The line and substation data for the main grid are as follows:
798 lines of 220 kV;
177 lines of 550kV lines;
main transformer 117 seat, total 107390MW of transformation capacity.
Results of examples
(1) The installed capacity of each type of unit is shown in fig. 2: 129 coal-fired units are provided, and the total capacity is 47615 MW; 2 gas units are provided, and the total capacity is 435 MW; 3 hydroelectric generating sets are provided, and the total capacity is 868 MW; the total number of the photovoltaic units is 10, and the total capacity is 3520 MW; the total number of the wind turbines is 50, and the total capacity is 6561.4 MW;
(2) the results of the annual output situation of each type of unit are shown in fig. 3:
the generated energy of the coal-fired unit is 195433911.00MWh, the utilization hours are 4104.46 hours, and the annual output accounts for 88.88%;
the generated energy of the gas turbine unit is 2095830.00MWh, the utilization hours are 4818 hours, and the annual output accounts for 0.95%;
the generating capacity of the hydroelectric generating set is 3896589.45MWh, the utilization hours are 4104.46 hours, and the annual output accounts for 1.77%;
the generated energy of the wind turbine generator is 13329774.16MWh, the utilization hours are 2031.54 hours, and the annual output accounts for 6.06%;
the generated energy of the photovoltaic unit is 5132932.98MWh, the utilization hours are 1458.22 hours, and the annual output accounts for 2.33 percent.
(3) The operating cost of the power system is as follows:
as shown in the table below, the annual operating cost of the power saving network system is 164.32 billion yuan.
Figure BDA0003336001350000151
(4) New energy consumption situation
According to an annual electric quantity balance meter of operation simulation output, the wind power curtailment electric quantity of the province is 26355.87MWh, the curtailment rate is 0.20%, and the wind power utilization hours are 2031.54 hours; the photovoltaic electricity rejection was 17946.82MWh, the electricity rejection was 0.35%, and the photovoltaic utilization hours was 1458.22 hours, with the results shown in fig. 4.
Figure BDA0003336001350000152
Figure BDA0003336001350000161
(5) Network congestion situations are shown in the following table:
Figure BDA0003336001350000162
the working principle is as follows: the method comprises the following steps of adopting a power system deterministic operation simulation model based on a time sequence load curve, wherein the core of the model is a scheduling simulation model taking daily operation as the core, and determining a unit to be put into operation according to an installation schedule by considering the operation, the retirement, the technical improvement and the like of the unit; next, the maintenance unit is removed according to the maintenance plan, and finally, the operable unit and parameters thereof are determined; then arranging all units capable of determining output, including external protocol power transmission, nuclear power units and units considering specified output, and correcting corresponding load curves according to the region where the power supply is located; according to the new energy simulated output generated by the new energy operation simulation module in a random simulation mode, arranging the new energy output and correcting a corresponding load curve; based on the corrected multi-region load curve, the pumped storage and conventional hydroelectric generating sets are arranged to carry out peak clipping and valley filling, the pumped storage can be set to be a flat pumping or full pumping mode, the constraints of capacity, electric quantity and the like of the generating sets are met, and the corresponding load curve is corrected again according to the region where the power supply is located. And finally, performing optimization simulation operation on the rest units, wherein the following parameters are required to be prepared besides the calculation result: the system comprises a manually-assigned unit state, positive and negative standby quantities in each time interval, time-sharing quotation or cost of the unit, start-stop cost of starting and stopping the unit, network constraint and the like.
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 attributes 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.
Furthermore, it should be understood that although the present description refers to embodiments, not every embodiment may contain only a single embodiment, and such description is for clarity only, and those skilled in the art should integrate the description, and the embodiments may be combined as appropriate to form other embodiments understood by those skilled in the art.

Claims (10)

1. A new energy power system coordinated operation simulation method is characterized in that: the method comprises the following steps:
step one, according to an installation schedule, considering unit commissioning, decommissioning, technical improvement and the like, and determining a commissioning unit; next, removing the maintenance unit according to the maintenance plan, and finally determining the operable unit and parameters thereof;
arranging all thermal power generating units capable of determining output, including external protocol power transmission, nuclear power generating units and units considering specified output, and correcting corresponding load curves according to the region where the power supply is located;
step three, according to the new energy simulation output of the new energy unit generated by the new energy operation simulation module in a random simulation mode, arranging the new energy output and correcting a corresponding load curve;
fourthly, on the basis of the corrected multi-region load curve, arranging peak clipping and valley filling for the pumped storage of the hydroelectric and pumped storage units and the conventional hydroelectric generating sets, setting the pumped storage mode to be a flat pumping or full pumping mode, meeting the constraints of capacity, electric quantity and the like of the units, and correcting the corresponding load curve again according to the region where the power supply is located;
and step five, performing optimization simulation operation on the rest units, wherein the following parameters are required to be prepared besides the calculation result: the system comprises a manually-assigned unit state, positive and negative spare amount of each time interval, partition spare amount, cross-partition spare amount, time-sharing quotation or cost of the unit, start-stop cost of starting and stopping the unit, network constraint and the like.
2. The method for simulating the coordinated operation of the new energy power system according to claim 1, characterized in that: the model adopts a power system deterministic operation simulation model based on a time sequence load curve, the core of the model is a scheduling simulation model taking daily operation as the core, and the objective function of the model is expressed as:
Figure FDA0003336001340000011
in the formula: t is the total time interval in the optimization cycle; c (P)t) For each type of unit, the output power is P in t periodtThe cost of operation; subscriptc、f、h、p、wRespectively indicating that thermal power cannot be started or stopped in a day, thermal power can be started or stopped in a day, water and electricity can be started or stopped in a day, and new energy can be pumped and stored; cwIn order to cut off the cost of new energy; pwd tSwitching new energy power for t time period; dd tCutting load power for t time period; vdLoad loss is cut for each node; cf、CcThe cost for starting and stopping the unit; θ, η, and γ are weighting coefficients, and are usually 1, and may be adjusted as necessary.
3. The method for simulating the coordinated operation of the new energy power system according to claim 1, characterized in that: in the second step, the output constraint on the output unit specifically includes:
the thermal power unit output upper and lower limits and the thermal power unit climbing restraint.
Figure FDA0003336001340000021
In the formula: pcmin,Pcmax,Pfmin,PfmaxRespectively the minimum output and the maximum output of the unit; ic is the state variable of the unit which can not be started or stopped in the day,Itfthe state variable is a state variable of a unit which can be started and stopped at a time t within a day;
Figure FDA0003336001340000022
Figure FDA0003336001340000023
the down-climbing speed and the up-climbing speed of the unit are respectively.
4. The method for simulating the coordinated operation of the new energy power system according to claim 1, characterized in that: in the third step, the output constraint of the new energy unit specifically comprises:
new energy output prediction variable introduced into new energy modeling
Figure FDA0003336001340000024
And a new energy cutting mechanism is introduced into the daily operation simulation model, so that the model can cut off part of new energy output under the conditions that the system cannot provide peak regulation capacity, the system spare capacity is insufficient or the new energy is blocked to be sent out.
Figure FDA0003336001340000031
In the formula:
Figure FDA0003336001340000032
new energy output for t time period;
Figure FDA0003336001340000033
switching new energy power for t time period;
Figure FDA0003336001340000034
and predicting a force value for the new energy in the t period.
Constraint of power abandonment rate of new energy
Figure FDA0003336001340000035
Wherein r iswdThe upper limit proportion of the power abandoning rate of the new energy is determined,
Figure FDA0003336001340000036
in order to switch the power sum of new energy,
Figure FDA0003336001340000037
the output of new energy is integrated.
5. The method for simulating the coordinated operation of the new energy power system according to claim 1, characterized in that: in the fourth step, the restriction of the water, electricity and pumped storage output specifically comprises:
the hydroelectric generating set optimizes and distributes the output of each time period in the model according to the output range and daily generated energy of the hydroelectric generating set given by the intermediate-term and long-term cross-basin cascade hydropower optimization scheduling result; the pumping and storage unit considers the balance between the daily pumping amount and the generated energy.
Figure FDA0003336001340000038
In the formula: phmin、PhmaxThe first line represents the upper and lower limit constraints of the output of the hydroelectric generating set; qhydroThe daily generated energy of the hydroelectric generating set is represented by the second row, the daily generated energy constraint of the hydroelectric generating set is represented by the second row, if the equation of the optimization result is not satisfied, the hydropower station generates water abandon, and the water abandon is converted into electric quantity
Figure FDA0003336001340000041
Pp,pump、Pp,genThe third line represents the upper and lower limit constraints of the output of the pumping and storage unit for the maximum pumping amount and the generated energy of the pumping and storage unit in a unit time period;
Figure FDA0003336001340000044
in order to describe the state variable of pumping or generating of the pumping and storage unit in the period t, the fourth line represents the mutual exclusion constraint of the pumping and generating states of the pumping and storage unit; lambda [ alpha ]pFor the efficiency of the pumping and storage unit, the fifth element represents the balance constraint of the pumping amount and the generating capacity of the pumping and storage unit each day.
6. The method for simulating the coordinated operation of the new energy power system according to claim 1, characterized in that: in the fifth step, the positive and negative reserve volume constraint in each time interval specifically includes:
Figure FDA0003336001340000042
in the formula: r isuFor the positive standby rate of the system required in time t, in the formula, the contribution of new energy to the system standby should be predicted to output P according to the predicted outputwf tCalculating that the cut-off part of the spare capacity is counted even if the spare capacity is cut off.
Figure FDA0003336001340000043
In the formula: r isd tThe minimum output of the new energy is considered to be 0 equivalently, and the minimum output of the new energy can be cut off at any time.
7. The method for simulating the coordinated operation of the new energy power system according to claim 1, characterized in that: in the fifth step, the partition spare amount constraint specifically includes:
the different regions need to ensure that the startup in each region reserves sufficient spare capacity for the load in the region, and the partition spare constraint can be expressed as:
Figure FDA0003336001340000051
in the formula: z is the total number of the areas; dzt,Is the load of the z region at time t; r isu zt,、rd zt,Respectively the positive standby rate and the negative standby rate of the z area at the time period t; z+And Z-Connecting line sets which are respectively sent in and sent out for the z area; f. ofl tIs the flow of the l lines. The first row of constraints represents the partition positive spare constraint and the second row represents the partition negative spare constraint.
8. The method for simulating the coordinated operation of the new energy power system according to claim 1, characterized in that: in the fifth step, the network constraint specifically includes:
establishing a line and section transmission constraint based on the direct current power flow model, as shown in the following formula:
Figure FDA0003336001340000052
in the formula:
Figure FDA0003336001340000053
respectively a t-period circuit and a section flow matrix; w is a generator transfer distribution factor; a. thengc，Angf，Angh，Angp，AngwRespectively node-unit incidence matrixes of different types of units; a. theslIs a section-line correlation matrix.
9. The method for simulating the coordinated operation of the new energy power system according to claim 1, characterized in that: in the fifth step, the cross-region power transmission constraint specifically includes:
in order to simulate a specific power transmission plan in actual power grid operation, cross-region power transmission constraint is introduced, namely, the transmission power flow of each time section of a section is limited.
Figure FDA0003336001340000061
Wherein Θ issThe line set contained in the cross-region power transmission section s is provided.
10. The method for simulating the coordinated operation of the new energy power system according to claim 1, characterized in that: and step five, dynamic constraint is further included, wherein the thermal power starting and stopping is carried out at most once in the running state in one day, and the large-capacity thermal power generating unit is not allowed to be started or stopped.
CN202111294116.7A 2021-11-03 2021-11-03 New energy power system coordinated operation simulation method Pending CN114243677A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111294116.7A CN114243677A (en) 2021-11-03 2021-11-03 New energy power system coordinated operation simulation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111294116.7A CN114243677A (en) 2021-11-03 2021-11-03 New energy power system coordinated operation simulation method

Publications (1)

Publication Number Publication Date
CN114243677A true CN114243677A (en) 2022-03-25

Family

ID=80743705

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111294116.7A Pending CN114243677A (en) 2021-11-03 2021-11-03 New energy power system coordinated operation simulation method

Country Status (1)

Country Link
CN (1) CN114243677A (en)

Similar Documents

Publication Publication Date Title
CN103580063B (en) A kind of method of large-scale grid connection wind-powered electricity generation of dissolving based on demanding party's response
CN107153885B (en) Real-time power generation plan optimization method considering deep peak shaving of thermal power generating unit
CN108092324B (en) AGC control system and control method for wind power participating in peak shaving frequency modulation
CN104377726B (en) A kind of peak regulating method of extensive new energy grid-connected power
CN107276122B (en) Peak-shaving resource calling decision method suitable for large-scale renewable energy grid connection
CN112395748B (en) Power system rotation reserve capacity optimization method considering supply and demand double-side flexible resources
CN115764927A (en) Power grid peak regulation method and system based on wind, light, water and fire multi-energy complementary characteristics
Wei et al. Energy-efficient distribution in smart grid
CN111030101A (en) Clean energy consumption linkage regulation and control method and system based on diversified big data
CN114243677A (en) New energy power system coordinated operation simulation method
CN115049431A (en) Pricing method of water and electricity in electric power spot market
CN114139362A (en) Intelligent soft switch optimal configuration method considering permeability of renewable energy
CN105244879A (en) Nuclear power water pumping energy-storage power station combined operation optimization method considering nuclear power unit maintenance
CN112801816A (en) Resource optimization scheduling method for total benefits of wind, light and water complementary system
CN111985844A (en) Day-ahead economic dispatching method for wind power and light energy comprehensive energy system
Xiao et al. Complementary coordinated dispatch of multi-energy systems considering power flexibility margin
CN112332456A (en) Optimal control method and device for heat accumulating type electric heating load to participate in power grid operation
Yun et al. Optimization of Operation Strategy of Virtual Power Plants Involved in Peak Shaving
CN112884191A (en) Thermoelectric day-ahead scheduling model based on network source coordination and calculation method
Weichun et al. Research on Joint Planning of the Source-Grid-Load with New Energy System Energy Balance and Economy
Zhang et al. Study on the microgrid pluripotent complementary and comprehensive economic optimization planning based on virtual energy theory
Feng et al. Research on Multi-Energy Complementary Clean DC Delivery Considering Sending-End Resources and the Matching Mode of Receiving-End Peak Load Regulation
Wang et al. Research on Source-load Synergy Multi-time Scale Optimal Scheduling Considering Large-scale Renewable Energy Access
Ming et al. Study on multi-source coordinated peak shaving strategy in Jiajiu area with high proportion of renewable energy
Liang et al. Study on deep-peak regulation of coal-gas-hydro power system in Portugal when considering high wind penetration

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