Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a power distribution network cooperative operation method for multi-investment subject and multi-element interaction, which comprises the following steps:
the method comprises the following steps: judging whether the power distribution network has unbalanced electric quantity, and entering a second step if the power distribution network has unbalanced electric quantity;
step two: the power distribution network operator sends interaction coordination information to each distributed energy investment operation main body;
step three: each distributed energy investment operation main body determines interactive electric quantity and interactive effect according to self operation characteristic constraint;
step four: the power distribution network operator collects the interactive response information returned by each distributed energy investment operation main body and determines the response sequence according to the benefit priority;
step five: the power distribution network operator collects the interactive response information returned by each distributed energy investment operation main body and determines the response sequence according to the benefit priority;
step six: the power distribution network operator determines a cooperative scheduling plan according to the response sequence and calculates the residual unbalanced electric quantity;
step seven: and when the residual unbalanced electric quantity is zero, the operation main bodies of all the distributed energy investment main bodies execute scheduling arrangement according to the collaborative interaction plan.
Further, the power distribution network operator in the second step sends interaction coordination information to each distributed energy investment operation subject, and sends incentive price to each distributed energy investment operation subject before sending the coordination information, wherein the incentive price is obtained by the following formula:
λt=(1+μt)Υbuy
wherein, Pt PVGThe total output of the distributed photovoltaic power generation units in the power distribution network in the t time period; pt WTGThe total output of the distributed wind power generation units in the power distribution network in the t period.
Furthermore, the cooperative operation and the interactive response of the power distribution network are realized by three types of distributed energy sources, namely energy storage, flexible load and interruptible load.
Further, the state of charge of the stored energy in the t period is represented as:
therein, SOCtThe state of charge of the stored energy in a time period t; pt EESThe absolute value of the power of the stored energy in the t period; capN EESCapacity for stored energy; bt CAnd bt DRespectively representing the 0-1 variable of the charge-discharge state of the stored energy in the t period, 1 represents that the state is effective, and 0 tableIndicating that the state is invalid ηEESCharge-discharge efficiency for energy storage; n is a radical ofmaxC&DThe maximum charge and discharge times allowed in one day for energy storage; t ist CAnd Tt DA set of time periods during which the stored energy is in a charged state and a discharged state within 24 hours; gamma raybuyTime-of-use electricity prices executed for distribution networks; lambda [ alpha ]t maxCAnd λt minDAre respectively at Tt CAnd Tt DThe maximum and minimum distribution network time of day prices within, which need to be updated at each time interval.
Further, the flexible load Pt FLThe actual power usage was determined using the following equation:
Pt FL=Pt planFL(λt/Υbuy)-ε,Pt minFL≤Pt FL≤Pt maxFL
wherein, Pt planFLPlanned power consumption, epsilon, for flexible loads in time period t>0 is a load elasticity factor, and the larger the value of the load elasticity factor is, the greater the flexibility of the flexible load is; pt minFLAnd Pt maxFLThe minimum and maximum electricity consumption of the flexible load in the period t are distinguished.
Further, said interruptible load Pt SLThe following formula is adopted:
πt SL=1/nSL,t∈Tt SL
wherein, PtotalSLIs the total power demand of the interruptible load during the day; pit SLIs the ratio of interruptible load electric quantity in total required electric quantity in the t period; t ist SLIs a set of interruptible load selectable time periods; n isSLThe number of time periods may be selected for the interruptible load.
Further, the interactive benefits of the two resource investment main bodies of the flexible load resource and the interrupt load resource are respectively as follows:
prt FL=Pt planFL(1-(λt/Υbuy)-ε)(λt-Υbuy)
prt SL=ΔPt SL|λt-Υbuy(tshifted)|
wherein, prt FLIs the interactive benefit of flexible load resources; pr (total reflection)t SLIs the interactive benefit of interruptible load resources; t is tshiftedIs a target period during which load shifting can be interrupted.
The invention has the beneficial effects that: the invention not only provides a brief interactive benefit quantification method for distributed energy investors, but also realizes the maximum consumption of renewable energy and the smoothness of a system net load curve through the multi-element interactive response among all distributed energy, thereby providing a method for the benefit distribution problem among multiple investment subjects in market reformation and promotion and the safe and stable operation problem of the power distribution network under the condition that the proportion of distributed energy is continuously increased.
Detailed Description
The technical solutions of the present invention are further described in detail below with reference to the accompanying drawings, but the scope of the present invention is not limited to the following.
As shown in figure 1, the invention provides a distributed energy multi-investment subject and multi-element interaction oriented power distribution network cooperative operation method, which mainly comprises a distributed energy investment subject interactive benefit quantification method and a power distribution network cooperative scheduling method based on price incentive.
Because the output upper limit of distributed renewable power generation such as wind power, photovoltaic and the like is completely controlled by natural factors such as real-time wind speed and illumination intensity of the area where the distributed renewable power generation is located, the distributed renewable power generation cannot be scheduled and arranged without operations such as wind abandoning and light abandoning, and the like, the cooperative operation and interactive response of the power distribution network need to be realized by means of three types of distributed energy such as energy storage, flexible load and interruptible load.
In order to maximize the benefit under the constraint of satisfying the maximum charging and discharging times per se, the energy storage needs to perform discharging and charging operations in the peak-valley period of the time-of-use electricity price, so as to earn a certain price difference, and the state of charge (SOC) of the energy storage in the t period is represented as:
therein, SOCtThe state of charge of the stored energy in a time period t; pt EESThe absolute value of the power of the stored energy in the t period; capN EESCapacity for stored energy; bt CAnd bt Dη are respectively a variable 0-1 representing the charging and discharging state of the stored energy in the period t, wherein 1 represents the state effective, and 0 represents the state ineffectiveEESCharge-discharge efficiency for energy storage; n is a radical ofmaxC&DThe maximum charge and discharge times allowed in one day for energy storage; t ist CAnd Tt DA set of time periods during which the stored energy is in a charging state and a discharging state within one day (24 h); gamma raybuyTime-of-use electricity prices executed for distribution networks; lambda [ alpha ]t maxCAnd λt minDAre respectively at Tt CAnd Tt DMaximum and minimum distribution network time-of-use electricity prices in the grid, which need to be updated at each time interval;
since the charging and discharging plan made by the investor of the stored energy for pursuing the benefit may slow down or aggravate the peak-valley difference of the power utilization curve, the power distribution network operator will provide a price signal to motivate the investor to participate in the blade interaction coordination operation in order to mobilize the stored energy to participate in the unified scheduling arrangement of the system.
In this case, the energy storage needs to first determine whether the energy storage can participate in the interaction and the corresponding interaction benefit according to the self-constraint conditions listed in equations (1) - (7) and in combination with the price signal provided by the power distribution network operator.
Specifically, assuming that the distribution network operator raises the incentive price for reducing the power usage during a certain period of time,
if the stored energy is in a discharge state within a certain period of time, the stored energy investment operator needs to consider whether the operator has residual capacity to participate in interaction, and if so, the operator has lambdat>λt minDThen the remaining portion of the discharge schedule would be executed ahead of time for this period to gain further benefit;
if the stored energy is in a charging state, the charging cost is higher than lambda due to the periodt maxCSo the energy storage investment operator will plan the charging for this period of time towardsAnd delaying and updating the subsequent charging and discharging plan. Thus, the interactive benefit of the energy storage investment entity can be calculated by the following formula:
wherein λ istProviding an interactive excitation price to the energy storage unit for the power distribution network; delta Pt CAnd Δ Pt DRespectively the response charging power and the response discharging power of the stored energy in the t period.
For load power demand, the load power demand can be decomposed into a basic fixed load P according to the controlled flexibility degreeFBLDaily random load Pt FDLFlexible load Pt FLAnd interruptible load Pt SLAnd fourthly, the method comprises the following steps.
The first part is the amount of power supply which must be provided to meet the basic life needs of the user, and can be approximately considered as a constant, such as lighting power consumption, monitoring equipment power consumption and the like;
the second part is a major source of uncertainty in the load, which appears to be constantly changing with the user's daily behavior, and is therefore a random quantity, such as television, elevator, computer power usage, etc.;
the third part refers to a type of controllable load which can be changed within a certain range, the total amount of the controllable load in each time period can be different, but the controllable load can be flexibly adjusted according to instructions, such as the electricity consumption of a thermal control air conditioner;
the last part is a type of load which has a fixed total power consumption in one day but can be flexibly arranged in the power consumption period, such as an intelligent washing machine, an electric automobile and the like. Thus, the demand for electricity during the period t may be expressed as:
Pt ED=PFBL+Pt FDL+Pt FL+Pt SL(9)
Pt FL=Pt planFL(λt/Υbuy)-ε,Pt minFL≤Pt FL≤Pt maxFL(10)
wherein, Pt planFLThe planned power consumption of the flexible load in the period t; epsilon>0 is a load elasticity factor, and the larger the value of the load elasticity factor is, the greater the flexibility of the flexible load is; pt minFLAnd Pt maxFLThe minimum and maximum electricity consumption of the flexible load in the time period t are distinguished; ptotalSLIs the total power demand of the interruptible load during the day; pit SLIs the ratio of interruptible load electric quantity in total required electric quantity in the t period; t ist SLIs a set of interruptible load selectable time periods; n isSLThe number of time periods may be selected for the interruptible load.
After obtaining the incentive price of the time interval from the power distribution network operator, the flexible load resource and interruptible load investment operator determines whether the flexible load resource and interruptible load investment operator can participate in the interactive response according to the formulas (9) to (11) or not.
The flexible load resource determines the actual power consumption according to the response characteristic formula (10) of the flexible load resource, the interactive behavior of the interruptible load resource is similar to the energy storage, the power utilization plan in the period is delayed to be implemented in the period with the lowest power price in the subsequent period after the power distribution network operator increases the incentive power price, and the power utilization plan in the period with the highest power price in the subsequent period is advanced to be implemented in the short time after the power distribution network operator decreases the incentive power price. Therefore, the interactive benefits of the two resource investment subjects are respectively as follows:
prt FL=Pt planFL(1-(λt/Υbuy)-ε)(λt-Υbuy) (13)
prt SL=ΔPt SL|λt-Υbuy(tshifted)| (14)
wherein, prt FLIs the interactive benefit of flexible load resources; pr (total reflection)t SLIs the interactive benefit of interruptible load resources; t is tshiftedIs a target period during which load shifting can be interrupted.
In order to improve the operation characteristics of the power distribution network, reduce reverse trend, smooth a net load curve and simultaneously realize optimized scheduling on the basis of guaranteeing the benefits of distributed energy investors and respecting the subjective intentions of the distributed energy investors, a power distribution network operator sends incentive price information to each distributed energy investment operation main body before each interactive coordination, and the power distribution network operator autonomously judges whether to participate in the interactive coordination and corresponding interactive benefits according to the provided information and the self operation characteristic constraint. Considering the inverse proportional relationship between incentive price and supply-demand ratio (SDR), the incentive price used can be calculated by the following formula:
λt=(1+μt)Υbuy(15)
wherein, Pt PVGThe total output of the distributed photovoltaic power generation units in the power distribution network in the t time period; pt WTGThe total output of the distributed wind power generation units in the power distribution network in the t period.
Therefore, the power distribution network cooperative operation problem can be constructed into the following double-layer planning model:
Lower:min|Pt EES+Pt ED-Pt PVG-Pt WTG|for t∈[24d-23,24d]
subjectto:(1)-(20) (21)
in order to meet the actual operation speed requirement of the engineering, the approximate optimal solution of the model is quickly solved by the benefit priority ranking criterion under the cost of solving precision within an acceptable range. In each time interval, after receiving interactive response information and incentive price from a power distribution network operator, each distributed energy investment operation main body determines whether to participate in the interaction according to the operation characteristics of the distributed energy investment operation main body, if so, the corresponding interactive electric quantity and interactive benefit are calculated and returned to the power distribution network operator, but after the power distribution network operator acquires feedback information of all the distributed energy investment operation main bodies, the distributed energy investment operation main bodies which wish to participate in the collaborative interaction are subjected to priority ranking according to the feedback interactive benefits from large to small, then the distributed energy investment operation main bodies respond to the interactive requests in sequence until all the distributed energy investment operation main bodies complete traversal or all unbalanced electric quantities obtain all satisfied positions, and the rest unbalanced electric quantities are still balanced to the medium-high voltage power distribution network.
The foregoing is illustrative of the preferred embodiments of this invention, and it is to be understood that the invention is not limited to the precise form disclosed herein and that various other combinations, modifications, and environments may be resorted to, falling within the scope of the concept as disclosed herein, either as described above or as apparent to those skilled in the relevant art. And that modifications and variations may be effected by those skilled in the art without departing from the spirit and scope of the invention as defined by the appended claims.