CN109726894B - New energy active instruction calculation method for guaranteeing spot transaction and medium and long-term electric quantity - Google Patents

New energy active instruction calculation method for guaranteeing spot transaction and medium and long-term electric quantity Download PDF

Info

Publication number
CN109726894B
CN109726894B CN201811463917.XA CN201811463917A CN109726894B CN 109726894 B CN109726894 B CN 109726894B CN 201811463917 A CN201811463917 A CN 201811463917A CN 109726894 B CN109726894 B CN 109726894B
Authority
CN
China
Prior art keywords
new energy
power
active
ith
electric quantity
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
CN201811463917.XA
Other languages
Chinese (zh)
Other versions
CN109726894A (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
State Grid Shanxi Electric Power Co Ltd
Nari Technology Co Ltd
Nanjing NARI Group Corp
Original Assignee
State Grid Corp of China SGCC
State Grid Shanxi Electric Power Co Ltd
Nari Technology Co Ltd
Nanjing NARI Group Corp
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, State Grid Shanxi Electric Power Co Ltd, Nari Technology Co Ltd, Nanjing NARI Group Corp filed Critical State Grid Corp of China SGCC
Priority to CN201811463917.XA priority Critical patent/CN109726894B/en
Publication of CN109726894A publication Critical patent/CN109726894A/en
Application granted granted Critical
Publication of CN109726894B publication Critical patent/CN109726894B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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
    • Y04S50/00Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
    • Y04S50/10Energy trading, including energy flowing from end-user application to grid

Abstract

The invention discloses a new energy active instruction calculation method for guaranteeing spot transactions and medium and long term electric quantity, and belongs to the technical field of operation and control of electric power systems. Aiming at the grid-connected active real-time control of a new energy station comprising complex electric quantity targets such as temporary spot transaction, basic guarantee electric quantity, medium and long term electric power transaction and the like, the invention provides that an active output instruction of the new energy station is decomposed into two stages for solving, the first stage considers comprehensive factors such as the power grid acceptance capacity, station performance indexes, sensitivity, temporary spot application quantity and the like and distributes optimal temporary spot application batch quantity for each station, when the final instruction of the station is solved in the second stage, the temporary spot transaction electric power calculated in the first stage is used as the lower limit constraint of the instruction, the medium and long term electric quantity is introduced into an objective function to complete the progress index, the temporary spot transaction electric quantity is stripped from the generated energy, the medium and long term electric quantity part is not influenced, and the fairness of market transaction is effectively guaranteed.

Description

New energy active instruction calculation method for guaranteeing spot transaction and medium and long-term electric quantity
Technical Field
The invention relates to a new energy station active control instruction two-stage calculation method for guaranteeing temporary spot transactions and medium and long-term electric quantity, and belongs to the technical field of operation and control of electric power systems.
Background
At present, with the deepening of the structure and the system of the electric power market in recent years, temporary spot market transaction becomes an important way for promoting the cross-regional consumption of new energy in various places. Although the new energy participating in the medium and long term electricity trading already has an electricity trading mechanism and a regulation and control principle under a corresponding long time scale, a corresponding processing strategy is lacked, the coordination of various regulation and control requirements such as temporary spot trading, medium and long term trading, guaranteed purchasing and the like is not well realized, and the regulation and control difficulty under the electric power market environment is improved by the temporary spot trading. The two electric quantity components of the medium-term and long-term output plan and the temporary spot application are coordinated, so that the balance of the electric power market is facilitated, and the new energy consumption level is improved.
With the rapid development of new energy, the construction of an electric power market system is deepened continuously, when extra-high voltage transaction is considered, the regulation and control object of a conventional new energy active control method is medium and long term electric power transaction generally, for the application of the appeared temporary spot electric quantity, the patent application number (201810245441.6) of 'an active real-time control method for considering medium and long term transaction and temporary spot transaction constraint' proposes that the transaction electric quantity execution index of each power plant is counted in real time, the grid-connected active power of each power plant is controlled in real time according to the transaction electric quantity execution index, and then fair index distribution is carried out on each power plant to complete the transaction plan electric quantity to the maximum extent. However, the patent does not consider the problem of coordination control between the temporary spot transaction application and the medium-and-long-term plan, and the supply-demand balance control of the power market is not refined enough.
According to the method, distribution of temporary spot transaction amount is used as a priority principle, then output instruction distribution is carried out on medium and long-term transaction targets, factors such as the predicted performance index, the adjusting performance index and the economic and environmental protection index of a new energy station and safety and stability constraints, peak regulation constraints, power plant adjusting plans, power grid power transmission channel constraints and the like which need to be met by power grid operation are considered, and the dispatching plans are corrected to meet the requirements of a real-time power market.
Disclosure of Invention
The purpose is as follows: in order to overcome the defects in the prior art and respond to the requirement of temporary spot electric quantity transaction application and give consideration to the long-term electric quantity fairness of the output of the new energy station, the invention provides a new energy active instruction calculation method for guaranteeing spot transaction and medium and long-term electric quantity.
The technical scheme is as follows: in order to solve the technical problems, the technical scheme adopted by the invention is as follows:
the basic principle of the invention is as follows: when the temporary current electric quantity transaction application is not considered, the optimization model of the new energy station grid-connected active power control achieves the purpose of maximum consumption of new energy output under the conditions of considering the requirements of power grid safety and stability, peak regulation constraint, power plant plan, transmission channel constraint and the like, and also considers economic and environmental protection factors. If temporary spot goods electric quantity transaction application needs to be considered, an active power output instruction of the new energy station is decomposed into two parts according to the unit output condition, the conventional load prediction requirement, the temporary spot goods transaction application, the power grid constraint and the electric quantity on-line quotation, wherein the first part is the temporary spot goods application batch quantity distributed by the station; and the second part is the active power output planned by the station aiming at the medium and long-term electric quantity, the optimization model is subjected to two-stage linear programming calculation, and then the optimal solutions calculated in the two stages are added to obtain the final active power output instruction of the station.
Specifically, the method for calculating the active instruction of the new energy for guaranteeing spot transaction and medium and long-term electric quantity comprises the following steps:
a new energy active instruction calculation method for guaranteeing spot transaction and medium and long-term electric quantity comprises the following steps:
step 1: let the current time be t 0 If at present t 0 The moment is the moment of the active control of the power grid, and the step 2 is entered; otherwise, continuing to wait until the current time t 0 The moment of active control of the power grid;
step 2: will t 0 The set of new energy stations participating in the output plan adjustment at any moment is recorded as A, and the power grid is divided byThe set of the conventional units except the new energy station is marked as B; defining the electric quantity calculation period as T, and recording the jth electric quantity calculation period as [ T w,j ,t w,j+1 ]Wherein t is w,j+1 =t w,j +T,t w,j For the starting point of the j-th electricity calculation cycle, t w,j+1 Calculating an end point of a cycle for the jth electrical quantity; recording the target value of the middle-long term electric quantity of the ith new energy station in the A within the middle-long term electric quantity progress check period as W p,i The initial time t of the 1 st electric quantity calculation cycle w,1 To the beginning time t of the j electric quantity calculation period w,j The sent power is recorded as
Figure BDA0001887429380000031
The set 1 st electric quantity calculation period starting time t w,1 To the beginning time t of the j electric quantity calculation period w,j The amount of electricity for the completed temporary spot transaction is recorded
Figure BDA0001887429380000032
Recording the middle-long term electric quantity progress assessment period of the ith new energy station in the A as T p,i Calculating by using a formula (1) to obtain the starting time t of the ith new energy station from the set 1 st electric quantity calculation period w,1 To the beginning time t of the j electric quantity calculation period w,j The medium and long term electric quantity completion progress index
Figure BDA0001887429380000033
The predicted performance index beta of the new energy station is obtained through calculation of a calculation formula (2) i (ii) a The regulation performance index gamma of the new energy station is obtained through calculation of a calculation formula (3) i Then entering step 3;
and step 3: for t 0 The station which has the temporary spot power application but is not repeated at the moment, an optimization model represented by a formula (4) is solved by adopting a linear programming method, the temporary spot power application repetition value of each new energy station is calculated, if the optimal solution can be obtained, the optimal solution is used as the basis of the second-stage instruction calculation, and then the step 4 is carried out;
and 4, step 4: calculating a grid-connected active control instruction of each power plant in the current control moment A and B by solving an optimization function represented by a formula (5), and then entering the step 5;
and 5: issuing the new energy station grid-connected active control instruction obtained by calculation to a corresponding new energy station for execution, and then counting [ t [ [ t ] w,j ,t w,j+1 ]And updating the medium-and-long-term generated electricity quantity and the temporary spot transaction electricity quantity of the ith new energy station in the time period according to formulas (6) and (7), updating the medium-and-long-term generated electricity quantity and the temporary spot transaction completed electricity quantity of the ith new energy station at the moment of finishing the jth electricity quantity calculation cycle, and iterating formulas according to formulas (8) and (9).
Preferably, the formulas (1), (2) and (3) are as follows:
Figure BDA0001887429380000041
wherein, T o,i Calculating the starting time t of the period for the set 1 st electric quantity w,1 To the beginning time t of the j electric quantity calculation period w,j The time of (d); t is a unit of p,i A middle-long term electric quantity examination period;
Figure BDA0001887429380000042
Figure BDA0001887429380000043
wherein M is i Is a distance t 0 The total number of data sampling points P of a prediction curve of the grid-connected active power in the free power generation time period of the grid-connected active power in the latest set period pre,i,k 、P i,t,k And P i,c Respectively obtaining a grid-connected active power predicted value, a grid-connected active power actual value and a starting capacity of a kth data sampling point of the ith new energy station; n is a radical of hydrogen i Is a distance t 0 Data of prediction curve of grid-connected active power in limited time period of grid-connected active power in set period nearest to momentTotal number of sample points, P i,t,k 、P i,t,k And P i,c Respectively is the grid-connected active instruction value, the grid-connected active power actual value and the starting capacity of the ith new energy station at the kth data sampling point.
Preferably, the formula (4) is as follows:
Figure BDA0001887429380000051
wherein A is a set of new energy stations; b is the set of the conventional units; n is the number of tie lines; delta P i,1 Applying a repeat instruction value for the temporary spot goods of the ith new energy station in the A for a superior dispatching mechanism; beta is a i Is the predicted performance index, gamma, of the ith new energy station i Is the regulation performance index, lambda, of the ith new energy station i Is an economic and environmental index of the ith new energy station; delta P d,i Temporary spot volume, P' submitted to superordinate scheduling authority for ith new energy site pre,i Is a predicted value, P, of the active power output of the ith new energy station at the next moment acc The temporary spot total transaction power which is approved by the upper dispatching mechanism is provided; s 1,j,i The active sensitivity of the variable quantity of the temporary spot order to the section j is obtained; s 2,j,i Active sensitivity of the ith conventional unit output variation in B to section j, P g,i Is the current active output value, P' of the ith conventional unit g,i The output instruction value of the ith conventional unit at the next moment is obtained; s 3,j,i The active sensitivity of the power delivered to section j for the ith link, P n,i For the current power value, P', of the power sent out for the ith link n,i Sending out the power value of the power at the next moment for the ith connecting line; l is the set of load nodes in the power grid, S 4,j,i Active sensitivity of load change to section j, P l,i For the active actual value, P' of the i-th node load in the grid l,i The active predicted value of the ith node load in the power grid is obtained; p is tl,j And P tl,j,max The actual active power value and the safety and stability limit of the section are respectively.
Preferably, the formula (5) is as follows:
Figure BDA0001887429380000061
wherein, P i 'is the active power output total command, P' of the ith new energy station at the next moment A pre,i The active output predicted value of the ith new energy station at the next moment A is obtained; b is t 0 The grid loss rate of the power grid at any moment; wherein, P g,i,max And P g,i,min The maximum technical output and the minimum technical output of the ith conventional unit in the B are respectively, v is the maximum allowable adjusting speed of the conventional unit, and delta t is an active control period; s 2,j,i The active sensitivity of the output change of the ith conventional unit in the B to the section j, P g,i Is the current active output value, P' of the ith conventional unit g,i The output instruction value of the ith conventional unit at the next moment is obtained; s. the 3,j,i The active sensitivity of the power delivered to section j for the ith link, P n,i For the current power value, P', of the power sent out for the ith link n,i Sending out the power value of the power at the next moment for the ith connecting line; l is the set of load nodes in the power grid, S 4,j,i Active sensitivity of load change to section j, P l,i For the active actual value, P' of the i-th node load in the grid l,i The active predicted value of the ith node load in the power grid is obtained; s 5,j,i And the active sensitivity of the total output variation of the ith new energy station to the section j is obtained. P tl,j And P tl,j,max Respectively the actual value of active power of the section and the safety and stability limit.
Preferably, the formula (6) is as follows:
Figure BDA0001887429380000071
preferably, the formula (7) is as follows:
Figure BDA0001887429380000072
preferably, the formula (8) is as follows:
Figure BDA0001887429380000073
preferably, the formula (9) is as follows:
Figure BDA0001887429380000074
has the beneficial effects that: compared with the effect of the technical scheme of the power generation dispatching plan management at the current stage, the new energy active instruction calculation method for guaranteeing spot transaction and medium and long term electric quantity has the advantages that two-stage solving strategies are set in the real-time control optimization of the grid-connected active power of the new energy station, the temporary spot transaction part is calculated, then the total instruction is calculated, and the temporary spot transaction electric power is used as the lower limit of the output of the new energy station, so that the temporary spot transaction is fully guaranteed to be executed.
In addition, after the temporary spot-shipment transaction electric quantity in the sent electric quantity of the new energy power station is partially deducted, the long-term electric quantity target completion progress is checked, the fairness of the new energy station in the electric quantity transaction check period is guaranteed, and the enthusiasm of the new energy station for participating in the temporary spot-shipment market is improved. In addition, the prediction performance index, the regulation performance index and the economic and environmental protection index of the new energy station are considered in the temporary spot goods approval stage and the real-time control stage, so that the technical level and the regulation performance of the new energy station are promoted to be improved.
Drawings
FIG. 1 is a schematic flow diagram of the process of the present invention.
Detailed Description
The present invention will be further described with reference to the accompanying drawings.
As shown in fig. 1, a method for calculating a new energy active instruction for guaranteeing spot transactions and medium and long term electric quantities includes the following steps:
step 1: let the current time be t 0 If at the present t 0 The moment is the active control moment of the power grid, and the step 2 is entered; otherwise, continuing to wait until the current time t 0 The moment of active control of the power grid;
step 2: will t 0 And a set of new energy stations participating in the output plan adjustment at any moment is marked as A, and a set of conventional units except the new energy stations in the power grid is marked as B. Defining the electric quantity calculation period as T, and recording the jth electric quantity calculation period as [ T w,j ,t w,j+1 ]Wherein t is w,j+1 =t w,j +T,t w,j For the starting point of the j-th electricity calculation cycle, t w,j+1 Calculating a termination point of a cycle for the jth electrical quantity; recording the medium-long term electric quantity target value of the ith new energy station in the A in the medium-long term electric quantity progress check period as W p,i The set 1 st electric quantity calculation period starting time t w,1 To the beginning time t of the j electric quantity calculation period w,j The sent power is recorded as
Figure BDA0001887429380000081
The initial time t of the 1 st electric quantity calculation period w,1 To the beginning time t of the j electric quantity calculation period w,j The completed temporary spot transaction electric quantity is recorded as
Figure BDA0001887429380000082
Recording the middle-long term electric quantity progress check period of the ith new energy station in the A as T p,i Calculating by using a formula (1) to obtain the starting time t of the ith new energy station from the set 1 st electric quantity calculation period w,1 To the beginning time t of the j electric quantity calculation period w,j The medium and long term electric quantity completion progress index
Figure BDA0001887429380000083
The predicted performance index beta of the new energy station is obtained through calculation of a calculation formula (2) i (ii) a The regulation performance index gamma of the new energy station is obtained through calculation of a calculation formula (3) i Then, go to step 3;
Figure BDA0001887429380000084
wherein, T o,i Calculating the starting time t of the period for the set 1 st electric quantity w,1 To the beginning time t of the j electric quantity calculation period w,j The time of (d); t is a unit of p,i The electric quantity examination period is a middle-long term electric quantity examination period.
Figure BDA0001887429380000091
Figure BDA0001887429380000092
Wherein M is i Is a distance t 0 The total number of data sampling points P of a prediction curve of the grid-connected active power in the free power generation time period of the grid-connected active power in the latest set period pre,i,k 、P i,t,k And P i,c The predicted value, the actual value and the starting capacity of the grid-connected active power of the kth data sampling point of the ith new energy station are respectively; n is a radical of hydrogen i Is a distance t 0 Total number of data sampling points, P' ″ of prediction curve of grid-connected active power in grid-connected active power limited time period in set cycle nearest to time i,t,k 、P i,t,k And P i,c Respectively is the grid-connected active instruction value, the grid-connected active power actual value and the starting capacity of the ith new energy station at the kth data sampling point.
And step 3: for t 0 The method comprises the steps that stations which have temporary spot power application but are not repeated at any moment exist, an optimization model represented by a formula (4) is solved by a linear programming method, the temporary spot power application repetition value of each new energy station is calculated, if the optimal solution can be obtained, the optimal solution is used as the basis of second-stage instruction calculation, and then the step 4 is carried out;
Figure BDA0001887429380000093
wherein A is a set of new energy stations; b is the set of the conventional units; n is the number of tie lines; delta P i,1 And applying a repeat instruction value for the temporary spot goods of the ith new energy station in the A for the superior dispatching mechanism. Beta is a i Is a predicted performance index, gamma, of the ith new energy station i Is the regulation performance index, lambda, of the ith new energy station i Is an economic and environmental protection index of the ith new energy station. Delta P d,i Temporary spot volume, P' submitted to superordinate scheduling authority for ith new energy site pre,i The predicted value P of the active power output of the ith new energy station at the next moment acc The temporary spot total transaction electric power which is batched by the upper-level dispatching mechanism is provided; s. the 1,j,i The active sensitivity of the variable quantity of the temporary spot order to the section j is obtained; s 2,j,i The active sensitivity of the output change of the ith conventional unit in the B to the section j, P g,i Is the current active power value, P g,i The output instruction value of the ith conventional unit at the next moment is obtained; s 3,j,i Sending out the active sensitivity of power to the section j for the ith link, P n,i For the current power value of the transmitted power of the ith link, P ″ n,i Sending a power value of the power at the next moment for the ith connecting line; l is the set of load nodes in the power grid, S 4,j,i Active sensitivity of load change to section j, P l,i For the active actual value, P' of the i-th node load in the grid l,i The active predicted value of the ith node load in the power grid is obtained; p tl,j And P tl,j,max Respectively the actual value of active power of the section and the safety and stability limit.
And 4, step 4: calculating a grid-connected active control instruction of each power plant in the current round of control time A and B by solving an optimization function represented by a formula (5), and then entering a step 5;
Figure BDA0001887429380000101
wherein, P i ' is the next momentActive power output total instruction, P pre,i The active power output predicted value of the ith new energy station at the next moment A is obtained; b is t 0 The grid loss rate of the power grid at any moment; wherein, P g,i,max And P g,i,min Respectively the maximum technical output and the minimum technical output of the ith conventional unit in the B, v is the maximum allowable adjusting speed of the conventional unit, and delta t is an active control period;
S 2,j,i active sensitivity of the ith conventional unit output variation in B to section j, P g,i Is the current active output value, P 'of the ith conventional unit' g,i The output instruction value of the ith conventional unit at the next moment is obtained; s. the 3,j,i Sending out the active sensitivity of power to the section j for the ith link, P n,i For the current power value, P', of the power sent out for the ith link n,i Sending out the power value of the power at the next moment for the ith connecting line; l is the set of load nodes in the power grid, S 4,j,i Active sensitivity of load change to section j, P l,i Active actual value, P' of the i-th node load in the grid l,i The active predicted value of the ith node load in the power grid is obtained; s. the 5,j,i And the active sensitivity of the total output variation of the ith new energy station to the section j is obtained. P tl,j And P tl,j,max The actual active power value and the safety and stability limit of the section are respectively.
And 5: issuing the new energy station grid-connected active control instruction obtained by calculation to a corresponding new energy station for execution, and then counting [ t [ [ t ] w,j ,t w,j+1 ]And (3) updating the medium-and-long-term generated electricity quantity and the temporary spot transaction electricity quantity of the ith new energy station in the time period, such as formulas (6) and (7), updating the medium-and-long-term generated electricity quantity and the temporary spot transaction completed electricity quantity of the ith new energy station at the moment of finishing the electricity quantity calculation period of the jth round, and iterating formulas, such as formulas (8) and (9), thereby finishing the method.
Figure BDA0001887429380000111
Figure BDA0001887429380000112
Figure BDA0001887429380000113
Figure BDA0001887429380000121
The above description is only of the preferred embodiments of the present invention, and it should be noted that: it will be apparent to those skilled in the art that various modifications and adaptations can be made without departing from the principles of the invention, and such modifications and adaptations are intended to be within the scope of the invention.

Claims (6)

1. A new energy active instruction calculation method for guaranteeing spot transaction and medium and long-term electric quantity is characterized by comprising the following steps: the method comprises the following steps:
step 1: let the current time be t 0 If at the present t 0 The moment is the active control moment of the power grid, and the step 2 is entered; otherwise, continuing to wait until the current time t 0 The active control moment is the power grid;
step 2: will t 0 The set of the new energy stations participating in the output plan adjustment at any moment is marked as A, and the set of the conventional units except the new energy stations in the power grid is marked as B; defining the electric quantity calculation period as T, recording the jth electric quantity calculation period as [ T ] w,j ,t w,j+1 ]Wherein t is w,j+1 =t w,j +T,t w,j For the starting point of the j-th electricity calculation cycle, t w,j+1 Calculating a termination point of a cycle for the jth electrical quantity; recording the target value of the middle-long term electric quantity of the ith new energy station in the A within the middle-long term electric quantity progress check period as W p,i The set 1 st electric quantity calculation period starting time t w,1 To the beginning time t of the j electric quantity calculation period w,j The amount of electricity sent is recorded as
Figure FDA0003708015000000011
The set 1 st electric quantity calculation period starting time t w,1 To the beginning time t of the j electric quantity calculation period w,j The completed temporary spot transaction electric quantity is recorded as
Figure FDA0003708015000000012
Recording the middle-long term electric quantity progress check period of the ith new energy station in the A as T p,i And calculating to obtain the starting time t of the ith new energy station from the set 1 st electric quantity calculation period w,1 To the beginning time t of the j electric quantity calculation period w,j The medium and long term electric quantity completion progress index
Figure FDA0003708015000000013
Calculating to obtain the predicted performance index beta of the new energy station i (ii) a Calculating to obtain the regulation performance index gamma of the new energy station i Then entering step 3;
and step 3: for t 0 The station which has the temporary spot power application but is not repeated at the moment, an optimization model represented by a formula (4) is solved by adopting a linear programming method, the temporary spot power application repetition value of each new energy station is calculated, if the optimal solution can be obtained, the optimal solution is used as the basis of the second-stage instruction calculation, and then the step 4 is carried out;
the formula (4) is as follows:
Figure FDA0003708015000000021
wherein A is a set of new energy stations; b is the set of the conventional units; n is the number of the connecting lines; delta P i,1 Applying a repeat instruction value for the temporary spot goods of the ith new energy station in the A for the upper dispatching mechanism; beta is a i Is a predicted performance index, gamma, of the ith new energy station i Is the regulation performance index, lambda, of the ith new energy station i Is an economic and environmental index of the ith new energy station; delta P d,i For the ith new energy station to the upper levelProvisional spot requisition volume, P' submitted by a scheduling authority pre,i The predicted value P of the active power output of the ith new energy station at the next moment acc The temporary spot total transaction power which is approved by the upper dispatching mechanism is provided; s 1,j,i The active sensitivity of the variable quantity of the temporary spot order to the section j is obtained; s 2,j,i The active sensitivity of the output change of the ith conventional unit in the B to the section j, P g,i Is the current active output value, P' of the ith conventional unit g,i The output instruction value at the next moment of the ith conventional unit is obtained; s 3,j,i The active sensitivity of the power delivered to section j for the ith link, P n,i For the current power value of the transmitted power of the ith link, P ″ n,i Sending a power value of the power at the next moment for the ith connecting line; l is the set of load nodes in the power grid, S 4,j,i Active sensitivity of load change to section j, P l,i Is the active actual value, P' of the ith node load in the power grid l,i The active predicted value of the ith node load in the power grid is obtained; p tl,j And P tl,j,max Respectively the actual value of active power and the safety and stability limit of the section;
and 4, step 4: calculating a grid-connected active control instruction of each power plant in the current round of control time A and B by solving an optimization function represented by a formula (5), and then entering a step 5; the formula (5) is as follows:
Figure FDA0003708015000000031
wherein, P i Is the active power output total command of the ith new energy station at the next moment A, P pre,i The active output predicted value of the ith new energy station at the next moment A is obtained; b is t 0 The grid loss rate of the power grid at the moment; wherein, P g,i,max And P g,i,min The maximum technical output and the minimum technical output of the ith conventional unit in the B are respectively, v is the maximum allowable adjusting speed of the conventional unit, and delta t is an active control period; s. the 2,j,i The active sensitivity of the output change of the ith conventional unit in the B to the section j is obtained,P g,i is the current active output value, P' of the ith conventional unit g,i The output instruction value of the ith conventional unit at the next moment is obtained; s. the 3,j,i Sending out the active sensitivity of power to the section j for the ith link, P n,i For the current power value, P', of the power sent out for the ith link n,i Sending out the power value of the power at the next moment for the ith connecting line; l is the set of load nodes in the power grid, S 4,j,i Active sensitivity of load change to section j, P l,i Is the active actual value, P' of the ith node load in the power grid l,i The active predicted value of the ith node load in the power grid is obtained; s 5,j,i And the active sensitivity of the total output variation of the ith new energy station to the section j is obtained. P is tl,j And P tl,j,max Respectively the actual value of active power and the safety and stability limit of the section;
and 5: issuing the new energy station grid-connected active control instruction obtained by calculation to a corresponding new energy station for execution, and then counting [ t [ [ t ] w,j ,t w,j+1 ]Middle-long term power generation capacity of ith new energy station in time period
Figure FDA0003708015000000041
And temporary spot transaction power
Figure FDA0003708015000000042
And updating the power already sent by the ith new energy station for a medium and long term till the end of the power calculation period of the jth round
Figure FDA0003708015000000043
And the temporary spot transaction has completed the amount of electricity
Figure FDA0003708015000000044
2. The method for calculating the active instruction of new energy for guaranteeing spot transactions and medium and long term electric quantity according to claim 1, wherein the method comprises the following steps: the above-mentioned
Figure FDA0003708015000000045
β i 、γ i The calculation formula is as follows:
Figure FDA0003708015000000046
wherein, T o,i Calculating the starting time t of the period for the set 1 st electric quantity w,1 To the beginning time t of the j electric quantity calculation period w,j The time of (d); t is p,i The electric quantity examination period is a medium-long term electric quantity examination period;
Figure FDA0003708015000000047
Figure FDA0003708015000000048
wherein, M i Is a distance t 0 The total number of data sampling points P of a prediction curve of the grid-connected active power in the grid-connected active power free power generation time period in the nearest set period pre,i,k 、P i,t,k And P i,c Respectively obtaining a grid-connected active power predicted value, a grid-connected active power actual value and a starting capacity of a kth data sampling point of the ith new energy station; n is a radical of i Is a distance t 0 The total number of data sampling points, P' of the prediction curve of the grid-connected active power in the limited grid-connected active power time period in the set period closest to the moment i,t,k 、P i,t,k And P i,c Respectively is the grid-connected active instruction value, the grid-connected active power actual value and the starting capacity of the ith new energy station at the kth data sampling point.
3. The method for calculating the active command of new energy for guaranteeing spot transactions and medium and long term electric quantity according to claim 2, wherein the method comprises the following steps: the above-mentioned
Figure FDA0003708015000000049
The calculation formula is as follows:
Figure FDA0003708015000000051
4. the method for calculating the active instruction of new energy for guaranteeing spot transactions and medium and long term electric quantity according to claim 2, wherein the method comprises the following steps: the above-mentioned
Figure FDA0003708015000000052
The calculation formula is as follows:
Figure FDA0003708015000000053
5. the method for calculating the active instruction of new energy for guaranteeing spot transactions and medium and long term electric quantity according to claim 2, wherein the method comprises the following steps: the above-mentioned
Figure FDA0003708015000000054
The calculation formula is as follows:
Figure FDA0003708015000000055
6. the method for calculating the active instruction of new energy for guaranteeing spot transactions and medium and long term electric quantity according to claim 2, wherein the method comprises the following steps: the described
Figure FDA0003708015000000056
The calculation formula is as follows:
Figure FDA0003708015000000057
CN201811463917.XA 2018-11-30 2018-11-30 New energy active instruction calculation method for guaranteeing spot transaction and medium and long-term electric quantity Active CN109726894B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811463917.XA CN109726894B (en) 2018-11-30 2018-11-30 New energy active instruction calculation method for guaranteeing spot transaction and medium and long-term electric quantity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811463917.XA CN109726894B (en) 2018-11-30 2018-11-30 New energy active instruction calculation method for guaranteeing spot transaction and medium and long-term electric quantity

Publications (2)

Publication Number Publication Date
CN109726894A CN109726894A (en) 2019-05-07
CN109726894B true CN109726894B (en) 2022-08-19

Family

ID=66294823

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811463917.XA Active CN109726894B (en) 2018-11-30 2018-11-30 New energy active instruction calculation method for guaranteeing spot transaction and medium and long-term electric quantity

Country Status (1)

Country Link
CN (1) CN109726894B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111478373B (en) * 2020-04-28 2021-07-13 国电南瑞科技股份有限公司 Active control method and system for new energy in consideration of medium-long-term transaction and real-time replacement
CN111835039B (en) * 2020-07-08 2021-09-14 国网宁夏电力有限公司 Cross-region increment spot goods access-oriented new energy active closed-loop control method and device
CN113469410B (en) * 2021-05-26 2023-10-31 国电南瑞科技股份有限公司 Online calculation method and device for maximum power consumption capacity of power grid in consideration of uncertainty of new energy

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106655279B (en) * 2016-08-03 2019-11-05 国电南瑞科技股份有限公司 The grid-connected active power real-time control method in power plant of meter and electricity trading program
CN108173261B (en) * 2017-12-25 2021-03-12 国电南瑞科技股份有限公司 New energy grid-connected active high-reliability automatic control system
CN108471145B (en) * 2018-03-20 2021-04-02 国网吉林省电力有限公司 Wind power plant active power control method based on multiple transaction plan virtual load rates

Also Published As

Publication number Publication date
CN109726894A (en) 2019-05-07

Similar Documents

Publication Publication Date Title
CN109492861B (en) Method for decomposing medium-term electricity quantity trading plan of cascade hydropower station group
CN109726894B (en) New energy active instruction calculation method for guaranteeing spot transaction and medium and long-term electric quantity
CN103683326B (en) A kind of regional power grid wind-powered electricity generation multiple spot accesses the computational methods of best receiving ability
CN108173293B (en) New energy active output maximum value on-line calculation method considering power grid operation constraint
CN108520315B (en) Power grid active real-time control method considering medium and long-term transaction and spot transaction constraints
CN109886567B (en) Short-term load prediction method considering somatosensory temperature and radiation intensity
CN110543696B (en) Method for small unmodeled unit to participate in electric power market clearing and safety check
CN107959285B (en) Optimization method, system and device for day-ahead power generation and transmission plan among regions
CN110472825B (en) Multistage scheduling mechanism coordinated step hydropower station real-time scheduling water abandoning reduction method
CN109245184B (en) Multi-source cooperative active control method suitable for multi-type constraint and multi-control interval
CN108054790B (en) Wind-solar power generation cluster active real-time optimization control method based on predicted output successive approximation
CN108964121B (en) Wind, light and water real-time control method considering water and power planning and power target in day before water and power
CN112909933B (en) Intraday rolling optimization scheduling method containing pumped storage unit under spot market environment
Li et al. A coordinated peak shaving strategy using neural network for discretely adjustable energy-intensive load and battery energy storage
CN114820046A (en) Regional power grid hybrid energy storage auxiliary frequency modulation economic optimization and compensation pricing method
CN114493222A (en) Wind power plant energy storage power station multi-market participation strategy considering output prediction and price
CN115441437A (en) Regional power grid day-ahead-day joint scheduling method based on deep reinforcement learning
CN111769599A (en) New energy active control method and device in electric power spot market
CN109829624B (en) Wind power cooperative game climbing control method and device
CN111130149A (en) Power grid power generation active control method and system considering performance distribution characteristics
CN107947246B (en) Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power
CN111008463A (en) Capacity allocation optimization method, device and equipment considering energy storage at power generation side
CN113794279A (en) Intelligent power grid system
CN105631549A (en) Virtual power plant distributed model prediction control method under active power distribution network environment
CN108471145B (en) Wind power plant active power control method based on multiple transaction plan virtual load rates

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