CN103633658B - A kind of energy-storage system output calculation method and device based on adjusting time constant filter in real time - Google Patents

Abstract

The invention provides a kind of energy-storage system output calculation method and device based on adjusting time constant filter in real time, be applied to wind storing cogeneration technical field, the method comprises: the energy-storage system state-of-charge calculating current control period; Calculate wind farm grid-connected 1min fluctuation ratio; Calculate wind farm grid-connected 10min fluctuation ratio; Determine the charging and discharging state of current control period energy-storage system; Adopt the time constant filter regulating formulae discovery current control period; First-order low-pass ripple is carried out to wind power output power, obtains the energy-storage system power output standard of current control period; Amplitude limit is carried out to the energy-storage system power output standard of described current control period, the energy-storage system power output standard of the current control period after being optimized; Utilize the energy-storage system power output standard of the current control period after described optimization, instruct energy-storage system to exert oneself.The present invention, under the prerequisite ensureing wind power smooth effect, improves the energy-storage system life-span.

Description

A kind of energy-storage system output calculation method and device based on adjusting time constant filter in real time

Technical field

The present invention relates to wind storing cogeneration technical field, particularly, relating to a kind of energy-storage system output calculation method and device based on adjusting time constant filter in real time.

Background technology

The feature such as the intrinsic interval of wind energy, Stochastic sum are uncontrollable, can cause negative effect to the quality of power supply of electrical network and stability.At present, State Grid Corporation of China has put into effect wind power system access electric power network technique regulation, proposes strict demand to the power fluctuation rate of grid connected wind power.At wind energy turbine set exit configuration certain capacity and the energy-storage system of power, can smooth wind power power fluctuation effectively, improve stability of power system.The current extensive energy-storage battery being suitable for wind energy turbine set has lithium battery, flow battery and sodium-sulphur battery etc.But energy-storage system cost is higher, the energy-storage battery capacity causing wind energy turbine set to configure and power are all very limited, and reasonable employment and protection energy-storage system, extending battery life have important economy and engineering significance.

In the level and smooth control strategy of existing energy-storage system, be mainly the methods such as first-order low-pass ripple, wavelet filtering.In first-order low-pass ripple algorithm, for the smooth effect fluctuated to wind-powered electricity generation, common method is simpler: if need smooth effect better, increase time constant filter, if less demanding to smooth effect, then reduces time constant filter.

Also can consider to protect the energy-storage system life-span in existing research.Common methods is on the basis of first-order low-pass ripple algorithm, under different energy-storage system state-of-charge (SOC, State of Charge), reasonably interval to ensure the energy-storage system life-span by regulating time constant filter to make energy-storage system be operated in as far as possible.But the impact on wind-powered electricity generation smooth effect under not considering the prerequisite in guarantee energy-storage system life-span in existing method.

Summary of the invention

The main purpose of the embodiment of the present invention is to provide a kind of energy-storage system output calculation method and device based on adjusting time constant filter in real time, to solve a kind of technology that wind power smooth effect can take into account again the raising energy-storage system life-span that can ensure.

To achieve these goals, the embodiment of the present invention provides a kind of energy-storage system output calculation method based on adjusting time constant filter in real time, comprising:

Go out the energy-storage system state-of-charge of force data, control cycle, last control cycle according to the energy-storage system of current control period, calculate the energy-storage system state-of-charge of current control period;

According to the energy storage reference output power of control cycle each in 1min before wind energy turbine set rated output power, current control period, calculate wind farm grid-connected 1min fluctuation ratio;

According to the energy storage reference output power of control cycle each in 10min before wind energy turbine set rated output power, current control period, calculate wind farm grid-connected 10min fluctuation ratio;

By the wind power of current control period and the energy storage reference output power of last control cycle, determine the charging and discharging state of current control period energy-storage system;

Adopt the following time constant filter regulating formulae discovery current control period:

T(t)＝T(t-1)+a[k _p1×(f ₁(t)-f _ref1)+k _d1×[(f ₁(t)-f ₁(t-1))-(f ₁(t-1)-f ₁(t-2))]]

+b[k _p10×(f ₁₀(t)-f _ref10)+k _d10×[(f ₁₀(t)-f ₁₀(t-1))-(f ₁₀(t-1)-f ₁₀(t-2))]]

Utilize the time constant filter of described current control period, first-order low-pass ripple is carried out to wind power output power, obtain the energy-storage system power output standard of current control period;

Amplitude limit is carried out to the energy-storage system power output standard of described current control period, the energy-storage system power output standard of the current control period after being optimized;

Utilize the energy-storage system power output standard of the current control period after described optimization, instruct energy-storage system to exert oneself;

In described adjustment formula, t is control cycle sequence number; T (t) is the time constant filter of t control cycle; f ₁t () is wind farm grid-connected 1min fluctuation ratio; f ₁₀t () is wind farm grid-connected 10min fluctuation ratio; f _ref1for 1min fluctuation ratio maximum permissible value; f _ref10for 10min fluctuation ratio maximum permissible value; A is 1min fluctuation ratio control weight parameter; B is 10min fluctuation ratio control weight parameter; k _p1for 1min fluctuation ratio proportional control parameter; k _p10for 10min fluctuation ratio proportional control parameter; k _d1for 1min fluctuation ratio differential controling parameters; k _d10for 10min fluctuation ratio differential controling parameters;

Wherein, a, b carry out value according to following setting rule:

Work as f ₁(t) > f _ref1and f ₁₀(t) > f _ref10or f ₁(t)≤f _ref1and f ₁₀(t)≤f _ref10time, a=b=0.5;

Work as f ₁(t) > f _ref1and f ₁₀(t)≤f _ref10time, a=0.7, b=0.3;

Work as f ₁(t)≤f _ref1and f ₁₀(t) > f _ref10time, a=0.3, b=0.7;

K _p1value is carried out according to following setting rule:

Work as f ₁(t) > f _ref1, SOC (t) < S _low, when energy-storage system is charged state, k _p1=W ₁× K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) < S _low, when energy-storage system is discharge condition, k _p1=W ₂× K _{p1_yes};

Work as f ₁(t) > f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p1=K _{p1_yes};

Work as f ₁(t) > f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p1=K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) > S _high, when energy-storage system is charged state, k _p1=W ₂× K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) > S _high, when energy-storage system is discharge condition, k _p1=W ₁× K _{p1_yes};

Work as f ₁(t)≤f _ref1, SOC (t) < S _low, when energy-storage system is charged state, k _p1=W ₂× K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) < S _low, when energy-storage system is discharge condition, k _p1=W ₁× K _{p1_no};

Work as f ₁(t)≤f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p1=K _{p1_no};

Work as f ₁(t)≤f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p1=K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) > S _high, when energy-storage system is charged state, k _p1=W ₁× K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) > S _high, when energy-storage system is discharge condition, k _p1=W ₂× K _{p1_no};

K _p10value is carried out according to following setting rule:

Work as f ₁₀(t) > f _ref10, SOC (t) < S _low, when energy-storage system is charged state, k _p10=W ₁× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) < S _low, when energy-storage system is discharge condition, k _p10=W ₂× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p10=K _{p10_yes};

Work as f ₁₀(t) > f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p10=K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) > S _high, when energy-storage system is charged state, k _p10=W ₂× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) > S _high, when energy-storage system is discharge condition, k _p10=W ₁× K _{p10_yes};

Work as f ₁₀(t)≤f _ref10, SOC (t) < S _low, when energy-storage system is charged state, k _p10=W ₂× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) < S _low, when energy-storage system is discharge condition, k _p10=W ₁× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p10=K _{p10_no};

Work as f ₁₀(t)≤f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p10=K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) > S _high, when energy-storage system is charged state, k _p10=W ₁× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) > S _high, when energy-storage system is discharge condition, k _p10=W ₂× K _{p10_no};

K _d1value is carried out according to following setting rule:

Work as f ₁(t)≤f _ref1time, k _d1=0;

Work as f ₁(t) > f _ref1time, k _d1get the first set point;

K _d10value is carried out according to following setting rule:

Work as f ₁₀(t)≤f _ref10time, k _d10=0;

Work as f ₁₀(t) > f _ref10time, k _d10get the second set point;

Wherein, SOC (t) the energy-storage system state-of-charge that is t control cycle; S _lowfor the Low threshold of energy-storage system state-of-charge; S _highfor the high threshold of energy-storage system state-of-charge; K _{p1_yes}for f ₁(t) > f _ref1time proportional control parameter base value; K _{p1_no}for f ₁(t)≤f _ref1time proportional control parameter base value; K _{p10_yes}for f ₁₀(t) > f _ref10time proportional control parameter base value; K _{p10_no}for f ₁₀(t)≤f _ref10time proportional control parameter base value; W ₁, W ₂for proportional control parameter regulating parameter, and W ₁> 1 > W ₂.

Accordingly, the present invention also provides a kind of energy-storage system output calculation device based on adjusting time constant filter in real time, comprising:

State-of-charge computing module, for going out the energy-storage system state-of-charge of force data, control cycle, last control cycle according to the energy-storage system of current control period, calculates the energy-storage system state-of-charge of current control period;

1min fluctuation ratio computing module, for the energy storage reference output power according to control cycle each in 1min before wind energy turbine set rated output power, current control period, calculates wind farm grid-connected 1min fluctuation ratio;

10min fluctuation ratio computing module, for the energy storage reference output power according to control cycle each in 10min before wind energy turbine set rated output power, current control period, calculates wind farm grid-connected 10min fluctuation ratio;

Discharge and recharge determination module, for by the wind power of current control period and the energy storage reference output power of last control cycle, determines the charging and discharging state of current control period energy-storage system;

Time constant filter computing module, for adopting the time constant filter of following adjustment formulae discovery current control period:

T(t)＝T(t-1)+a[k _p1×(f ₁(t)-f _ref1)+k _d1×[(f ₁(t)-f ₁(t-1))-(f ₁(t-1)-f ₁(t-2))]]

+b[k _p10×(f ₁₀(t)-f _ref10)+k _d10×[(f ₁₀(t)-f ₁₀(t-1))-(f ₁₀(t-1)-f ₁₀(t-2))]]

First-order low-pass mode block, for utilizing the time constant filter of described current control period, carries out first-order low-pass ripple to wind power output power, obtains the energy-storage system power output standard of current control period;

Optimize module, for carrying out amplitude limit to the energy-storage system power output standard of described current control period, the energy-storage system power output standard of the current control period after being optimized;

Instructing module, for utilizing the energy-storage system power output standard of the current control period after described optimization, instructing energy-storage system to exert oneself;

In described adjustment formula, t is control cycle sequence number; T (t) is the time constant filter of t control cycle; f ₁t () is wind farm grid-connected 1min fluctuation ratio; f ₁₀t () is wind farm grid-connected 10min fluctuation ratio; f _ref1for 1min fluctuation ratio maximum permissible value; f _ref10for 10min fluctuation ratio maximum permissible value; A is 1min fluctuation ratio control weight parameter; B is 10min fluctuation ratio control weight parameter; k _p1for 1min fluctuation ratio proportional control parameter; k _p10for 10min fluctuation ratio proportional control parameter; k _d1for 1min fluctuation ratio differential controling parameters; k _d10for 10min fluctuation ratio differential controling parameters;

Wherein, a, b carry out value according to following setting rule:

Work as f ₁(t) > f _ref1and f ₁₀(t) > f _ref10or f ₁(t)≤f _ref1and f ₁₀(t)≤f _ref10time, a=b=0.5;

Work as f ₁(t) > f _ref1and f ₁₀(t)≤f _ref10time, a=0.7, b=0.3;

Work as f ₁(t)≤f _ref1and f ₁₀(t) > f _ref10time, a=0.3, b=0.7;

K _p1value is carried out according to following setting rule:

Work as f ₁(t) > f _ref1, SOC (t) < S _low, when energy-storage system is charged state, k _p1=W ₁× K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) < S _low, when energy-storage system is discharge condition, k _p1=W ₂× K _{p1_yes};

Work as f ₁(t) > f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p1=K _{p1_yes};

Work as f ₁(t) > f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p1=K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) > S _high, when energy-storage system is charged state, k _p1=W ₂× K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) > S _high, when energy-storage system is discharge condition, k _p1=W ₁× K _{p1_yes};

Work as f ₁(t)≤f _ref1, SOC (t) < S _low, when energy-storage system is charged state, k _p1=W ₂× K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) < S _low, when energy-storage system is discharge condition, k _p1=W ₁× K _{p1_no};

Work as f ₁(t)≤f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p1=K _{p1_no};

Work as f ₁(t)≤f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p1=K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) > S _high, when energy-storage system is charged state, k _p1=W ₁× K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) > S _high, when energy-storage system is discharge condition, k _p1=W ₂× K _{p1_no};

K _p10value is carried out according to following setting rule:

Work as f ₁₀(t) > f _ref10, SOC (t) < S _low, when energy-storage system is charged state, k _p10=W ₁× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) < S _low, when energy-storage system is discharge condition, k _p10=W ₂× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p10=K _{p10_yes};

Work as f ₁₀(t) > f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p10=K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) > S _high, when energy-storage system is charged state, k _p10=W ₂× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) > S _high, when energy-storage system is discharge condition, k _p10=W ₁× K _{p10_yes};

Work as f ₁₀(t)≤f _ref10, SOC (t) < S _low, when energy-storage system is charged state, k _p10=W ₂× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) < S _low, when energy-storage system is discharge condition, k _p10=W ₁× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p10=K _{p10_no};

Work as f ₁₀(t)≤f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p10=K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) > S _high, when energy-storage system is charged state, k _p10=W ₁× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) > S _high, when energy-storage system is discharge condition, k _p10=W ₂× K _{p10_no};

K _d1value is carried out according to following setting rule:

Work as f ₁(t)≤f _ref1time, k _d1=0;

Work as f ₁(t) > f _ref1time, k _d1get the first set point;

K _d10value is carried out according to following setting rule:

Work as f ₁₀(t)≤f _ref10time, k _d10=0;

Work as f ₁₀(t) > f _ref10time, k _d10get the second set point;

Wherein, SOC (t) the energy-storage system state-of-charge that is t control cycle; S _lowfor the Low threshold of energy-storage system state-of-charge; S _highfor the high threshold of energy-storage system state-of-charge; K _{p1_yes}for f ₁(t) > f _ref1time proportional control parameter base value; K _{p1_no}for f ₁(t)≤f _ref1time proportional control parameter base value; K _{p10_yes}for f ₁₀(t) > f _ref10time proportional control parameter base value; K _{p10_no}for f ₁₀(t)≤f _ref10time proportional control parameter base value; W ₁, W ₂for proportional control parameter regulating parameter, and W ₁> 1 > W ₂.

By means of technique scheme, the present invention makes wind-electricity integration fluctuation ratio decline by energy-storage system, meet the requirement that grid company fluctuates to wind-powered electricity generation, improve stability of power system, under the prerequisite ensureing wind power smooth effect, invention increases the energy-storage system life-span, energy storage can be utilized to greatest extent, compared to prior art, the invention solves the contradiction between smooth effect and energy-storage system life-span in wind power smoothing process.

Accompanying drawing explanation

In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, below the accompanying drawing used required in describing embodiment is briefly described, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.

Fig. 1 is the schematic flow sheet of the energy-storage system output calculation method based on real-time adjustment time constant filter provided by the invention;

Fig. 2 is the wind storage associating output smoothing effect that the embodiment of the present invention one provides;

Fig. 3 is that the 1min fluctuation ratio that the embodiment of the present invention one provides improves situation;

Fig. 4 is that the 10min fluctuation ratio that the embodiment of the present invention one provides improves situation;

Fig. 5 is the energy storage power curve that the embodiment of the present invention one provides;

Fig. 6 is the state-of-charge variation effect comparison diagram that the embodiment of the present invention one provides;

Fig. 7 is the structural representation of the energy-storage system output calculation device based on real-time adjustment time constant filter provided by the invention.

Embodiment

Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.

Time constant filter is first-order low-pass wave parameter, and the larger wind farm grid-connected power of time constant filter is more level and smooth, and corresponding energy-storage system is exerted oneself larger, otherwise grid-connected power fluctuation increases, and energy-storage system is exerted oneself reduction.

Energy-storage system exert oneself excessive or long-time running at state-of-charge the excessive or too small stage, the energy-storage system life-span can be affected.For making wind farm grid-connected power fluctuation rate meet certain requirement, wishing that the energy-storage system life-span increases simultaneously, needing to adjust time constant filter in real time.

Need to follow following principle to the adjustment of time constant filter:

1. do not meet grid-connected requirement when wind energy turbine set fluctuation ratio, i.e. wind farm grid-connected 1min fluctuation ratio, the 10min fluctuation ratio respective limit value that had to exceed or both exceeded limit value, now need to increase time constant filter and go out force value to improve energy storage;

2., when energy-storage system state-of-charge is too high, limit its charging, encourage electric discharge; When energy-storage system state-of-charge is too low, limits its electric discharge, encourage charging.

According to above principle, the invention provides the following energy-storage system output calculation method based on adjusting time constant filter in real time, by adjusting time constant filter in real time, under the prerequisite ensureing wind power smooth effect, improve the life-span of energy-storage system, to solve the contradiction in wind power smoothing process between smooth effect and energy-storage system life-span.As shown in Figure 1, the method comprises:

Step S11, goes out the energy-storage system state-of-charge of force data, control cycle, last control cycle according to the energy-storage system of current control period, calculate the energy-storage system state-of-charge of current control period.

In a kind of preferred embodiment, this step adopts following formula:

SOC(t)＝SOC(t-1)+P _B(t)×Δt；

Wherein, t is control cycle sequence number; P _bt () goes out force data for the energy-storage system of t control cycle; The energy-storage system state-of-charge that SOC (t) is t control cycle; Δ t is control cycle, preferably, and Δ t=1s.

Step S12, according to the energy storage reference output power of control cycle each in 1min before wind energy turbine set rated output power, current control period, calculates wind farm grid-connected 1min fluctuation ratio.

In a kind of preferred embodiment, this step adopts following formula:

$f_1\left(t\right)=\frac{\underset{i=\mathrm{0,1,2},...,59}{\max }{P}_{\mathrm{ref}}(t-i)-\underset{i=\mathrm{0,1,2},...,59}{\min }{P}_{\mathrm{ref}}(t-i)}{P_{w\_\mathrm{rate}}};$

Wherein, f ₁t () is wind farm grid-connected 1min fluctuation ratio; P _reft energy storage reference output power that () is t control cycle; P _{w_rate}for wind energy turbine set rated output power.

Step S13, according to the energy storage reference output power of control cycle each in 10min before wind energy turbine set rated output power, current control period, calculates wind farm grid-connected 10min fluctuation ratio.

In a kind of preferred embodiment, this step adopts following formula:

$f_{10}\left(t\right)=\frac{\underset{i=\mathrm{0,1,2},...,599}{\max }{P}_{\mathrm{ref}}(t-i)-\underset{i=\mathrm{0,1,2},...,599}{\min }{P}_{\mathrm{ref}}(t-i)}{P_{w\_\mathrm{rate}}};$

Wherein, f ₁₀t () is wind farm grid-connected 10min fluctuation ratio; P _reft energy storage reference output power that () is t control cycle; P _{w_rate}for wind energy turbine set rated output power.

Step S14, by the wind power of current control period and the energy storage reference output power of last control cycle, determines the charging and discharging state of current control period energy-storage system.

In a kind of preferred embodiment, this step is according to following rule:

If when the wind power of current control period is greater than the energy storage reference output power of last control cycle, determine that current control period energy-storage system is charged state;

If when the wind power of current control period is less than the energy storage reference output power of last control cycle, determine that current control period energy-storage system is discharge condition.

Step S15, adopts the following time constant filter regulating formulae discovery current control period:

T(t)＝T(t-1)+a[k _p1×(f ₁(t)-f _ref1)+k _d1×[(f ₁(t)-f ₁(t-1))-(f ₁(t-1)-f ₁(t-2))]]

+b[k _p10×(f ₁₀(t)-f _ref10)+k _d10×[(f ₁₀(t)-f ₁₀(t-1))-(f ₁₀(t-1)-f ₁₀(t-2))]]

In this adjustment formula, T (t) is the time constant filter of t control cycle; f _ref1for 1min fluctuation ratio maximum permissible value; f _ref10for 10min fluctuation ratio maximum permissible value; A is 1min fluctuation ratio control weight parameter; B is 10min fluctuation ratio control weight parameter; k _p1for 1min fluctuation ratio proportional control parameter; k _p10for 10min fluctuation ratio proportional control parameter; k _d1for 1min fluctuation ratio differential controling parameters; k _d10for 10min fluctuation ratio differential controling parameters.

Wherein, a, b carry out value according to following setting rule:

Work as f ₁(t) > f _ref1and f ₁₀(t) > f _ref10or f ₁(t)≤f _ref1and f ₁₀(t)≤f _ref10time, a=b=0.5;

Work as f ₁(t) > f _ref1and f ₁₀(t)≤f _ref10time, a=0.7, b=0.3;

Work as f ₁(t)≤f _ref1and f ₁₀(t) > f _ref10time, a=0.3, b=0.7.

K _p1and k _p10for proportional control parameter, represent the proportional component that fluctuation ratio regulates time constant filter, its setting principle is: 1., when fluctuation ratio exceedes maximum fluctuation permissible value, proportional control parameter is larger; When fluctuation ratio is no more than maximum fluctuation permissible value, proportional control parameter is less; 2. on the basis of principle 1, when fluctuation ratio exceedes maximum fluctuation permissible value, because scale parameter unquote internal symbol is just, so when energy-storage system state-of-charge is lower, encourage charging (increase of proportional control parameter), restriction electric discharge (reduction of proportional control parameter), when energy-storage system state-of-charge is higher, encourage electric discharge (increase of proportional control parameter), restriction charging (reduction of proportional control parameter); When fluctuation ratio is no more than maximum fluctuation permissible value, because scale parameter unquote internal symbol is negative, so when energy-storage system state-of-charge is lower, encourage charging (reduction of proportional control parameter), restriction electric discharge (increase of proportional control parameter), when energy-storage system state-of-charge is higher, encourage electric discharge (reduction of proportional control parameter), restriction charging (increase of proportional control parameter).

If S _lowfor the Low threshold of energy-storage system state-of-charge; S _highfor the high threshold of energy-storage system state-of-charge; K _{p1_yes}for f ₁(t) > f _ref1time proportional control parameter base value; K _{p1_no}for f ₁(t)≤f _ref1time proportional control parameter base value; K _{p10_yes}for f ₁₀(t) > f _ref10time proportional control parameter base value; K _{p10_no}for f ₁₀(t)≤f _ref10time proportional control parameter base value; W ₁, W ₂for proportional control parameter regulating parameter, and W ₁> 1 > W ₂, then:

K _p1value is carried out according to following setting rule:

Work as f ₁(t) > f _ref1, SOC (t) < S _low, when energy-storage system is charged state, k _p1=W ₁× K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) < S _low, when energy-storage system is discharge condition, k _p1=W ₂× K _{p1_yes};

Work as f ₁(t) > f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p1=K _{p1_yes};

Work as f ₁(t) > f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p1=K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) > S _high, when energy-storage system is charged state, k _p1=W ₂× K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) > S _high, when energy-storage system is discharge condition, k _p1=W ₁× K _{p1_yes};

Work as f ₁(t)≤f _ref1, SOC (t) < S _low, when energy-storage system is charged state, k _p1=W ₂× K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) < S _low, when energy-storage system is discharge condition, k _p1=W ₁× K _{p1_no};

Work as f ₁(t)≤f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p1=K _{p1_no};

Work as f ₁(t)≤f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p1=K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) > S _high, when energy-storage system is charged state, k _p1=W ₁× K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) > S _high, when energy-storage system is discharge condition, k _p1=W ₂× K _{p1_no};

K _p10value is carried out according to following setting rule:

Work as f ₁₀(t) > f _ref10, SOC (t) < S _low, when energy-storage system is charged state, k _p10=W ₁× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) < S _low, when energy-storage system is discharge condition, k _p10=W ₂× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p10=K _{p10_yes};

Work as f ₁₀(t) > f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p10=K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) > S _high, when energy-storage system is charged state, k _p10=W ₂× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) > S _high, when energy-storage system is discharge condition, k _p10=W ₁× K _{p10_yes};

Work as f ₁₀(t)≤f _ref10, SOC (t) < S _low, when energy-storage system is charged state, k _p10=W ₂× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) < S _low, when energy-storage system is discharge condition, k _p10=W ₁× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p10=K _{p10_no};

Work as f ₁₀(t)≤f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p10=K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) > S _high, when energy-storage system is charged state, k _p10=W ₁× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) > S _high, when energy-storage system is discharge condition, k _p10=W ₂× K _{p10_no};

K _d1and k _d10for differential controling parameters, represent that its Controlling principle is: when fluctuation ratio does not exceed maximum permission undulating value, do not carry out differential control, and namely differential controling parameters is 0 because when fluctuation ratio is greater than maximum permission undulating value, needs additionally increase the adjustment to fluctuation ratio tendency; When fluctuation ratio exceedes maximum permission undulating value, differential controls to be set as definite value.

Concrete, k _d1value is carried out according to following setting rule:

Work as f ₁(t)≤f _ref1time, k _d1=0;

Work as f ₁(t) > f _ref1time, k _d1get the first set point.

Concrete, k _d10value is carried out according to following setting rule:

Work as f ₁₀(t)≤f _ref10time, k _d10=0;

Work as f ₁₀(t) > f _ref10time, k _d10get the second set point.

Step S16, utilizes the time constant filter of described current control period, carries out first-order low-pass ripple to wind power output power, obtains the energy-storage system power output standard of current control period.

First-order low-pass ripple frequency-domain expression is: $P_{\mathrm{ref}}\left(s\right)=\frac{1}{1+\mathrm{sT}}{P}_w\left(s\right);$

Wherein, P _refs () is the wind storing cogeneration reference output power in frequency domain.

Frequency-domain expression is converted to discrete time-domain expression formula to calculate:

$P_{\mathrm{ref}}\left(t\right)=\frac{T\left(t\right)}{T\left(t\right)+\mathrm{\&Delta;t}}{P}_{\mathrm{ref}}(t-1)+\frac{\mathrm{\&Delta;t}}{T\left(t\right)+\mathrm{\&Delta;t}}{p}_W\left(t\right);$

And then obtain energy-storage system power output standard:

$P_B^{\&prime;}\left(t\right)=\frac{T\left(t\right)}{T\left(t\right)+\mathrm{\&Delta;t}}[P_{\mathrm{ref}}(t-1)-P_w\left(t\right)];$

Wherein, P ' _bt energy-storage system power output standard that () is t control cycle; P _wt () is the wind power of t control cycle.

Step S17, carries out amplitude limit to the energy-storage system power output standard of described current control period, the energy-storage system power output standard of the current control period after being optimized.

Because energy-storage system has maximum charge-discharge electric power to limit P _max, so energy-storage system power output instruction demand fulfillment scope [-P _max, P _max], concrete operations are as follows:

As P ' _b(t) > P _maxtime, P _{b_cmd}(t)=P _max;

As-P _max≤ P ' _b(t)≤P _maxtime, P _{b_cmd}(t)=P ' _b(t);

As P ' _b(t) <-P _maxtime, P _{b_cmd}(t)=-P _max;

Wherein, P _{b_cmd}(t) for optimize after t control cycle energy-storage system power output standard.

Step S18, utilizes the energy-storage system power output standard of the current control period after described optimization, instructs energy-storage system to exert oneself.

Embodiment one

In the present embodiment, wind energy turbine set capacity is 100MW, and energy-storage system is 20MW/100MWh, 1min fluctuation ratio maximum permissible value is f _ref1=3%, 10min fluctuation ratio maximum permissible value is f _ref10=10%, the Low threshold S of energy-storage system state-of-charge _low=30%, the high threshold S of energy-storage system state-of-charge _high=70%.

Fig. 2 is Power Output for Wind Power Field smooth effect, compared with exerting oneself separately with wind-powered electricity generation as seen, wind storage combine exert oneself more level and smooth; Fig. 3 is corresponding 1min fluctuation ratio situation of change, and visible 1min fluctuation ratio has clear improvement, and within reaching limit value 3%; Fig. 4 is corresponding 10min fluctuation ratio situation of change, and visible 10min fluctuation ratio also has clear improvement, and within reaching limit value 10%; Fig. 5 is corresponding energy storage power curve; The state-of-charge variation effect using this method to regulate time constant filter and existing first-order low-pass ripple to adopt fixing time constant filter to obtain in real time contrasts by Fig. 6, visible, this method makes energy-storage system state-of-charge more operate in 30%-70% interval, makes the operation that energy-storage system is more healthy.

Accordingly, the present invention also provides a kind of energy-storage system output calculation device based on adjusting time constant filter in real time, and as shown in Figure 7, this device comprises:

State-of-charge computing module 71, for going out the energy-storage system state-of-charge of force data, control cycle, last control cycle according to the energy-storage system of current control period, calculates the energy-storage system state-of-charge of current control period;

1min fluctuation ratio computing module 72, for the energy storage reference output power according to control cycle each in 1min before wind energy turbine set rated output power, current control period, calculates wind farm grid-connected 1min fluctuation ratio;

10min fluctuation ratio computing module 73, for the energy storage reference output power according to control cycle each in 10min before wind energy turbine set rated output power, current control period, calculates wind farm grid-connected 10min fluctuation ratio;

Discharge and recharge determination module 74, for by the wind power of current control period and the energy storage reference output power of last control cycle, determines the charging and discharging state of current control period energy-storage system;

Time constant filter computing module 75, for adopting the time constant filter of following adjustment formulae discovery current control period:

T(t)＝T(t-1)+a[k _p1×(f ₁(t)-f _ref1)+k _d1×[(f ₁(t)-f ₁(t-1))-(f ₁(t-1)-f ₁(t-2))]]

+b[k _p10×(f ₁₀(t)-f _ref10)+k _d10×[(f ₁₀(t)-f ₁₀(t-1))-(f ₁₀(t-1)-f ₁₀(t-2))]]

First-order low-pass mode block 76, for utilizing the time constant filter of described current control period, carries out first-order low-pass ripple to wind power output power, obtains the energy-storage system power output standard of current control period;

Optimize module 77, for carrying out amplitude limit to the energy-storage system power output standard of described current control period, the energy-storage system power output standard of the current control period after being optimized;

Instructing module 78, for utilizing the energy-storage system power output standard of the current control period after described optimization, instructing energy-storage system to exert oneself;

In described adjustment formula, t is control cycle sequence number; T (t) is the time constant filter of t control cycle; f ₁t () is wind farm grid-connected 1min fluctuation ratio; f ₁₀t () is wind farm grid-connected 10min fluctuation ratio; f _ref1for 1min fluctuation ratio maximum permissible value; f _ref10for 10min fluctuation ratio maximum permissible value; A is 1min fluctuation ratio control weight parameter; B is 10min fluctuation ratio control weight parameter; k _p1for 1min fluctuation ratio proportional control parameter; k _p10for 10min fluctuation ratio proportional control parameter; k _d1for 1min fluctuation ratio differential controling parameters; k _d10for 10min fluctuation ratio differential controling parameters;

Wherein, a, b carry out value according to following setting rule:

Work as f ₁(t) > f _ref1and f ₁₀(t) > f _ref10or f ₁(t)≤f _ref1and f ₁₀(t)≤f _ref10time, a=b=0.5;

Work as f ₁(t) > f _ref1and f ₁₀(t)≤f _ref10time, a=0.7, b=0.3;

Work as f ₁(t)≤f _ref1and f ₁₀(t) > f _ref10time, a=0.3, b=0.7;

K _p1value is carried out according to following setting rule:

Work as f ₁(t) > f _ref1, SOC (t) < S _low, when energy-storage system is charged state, k _p1=W ₁× K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) < S _low, when energy-storage system is discharge condition, k _p1=W ₂× K _{p1_yes};

Work as f ₁(t) > f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p1=K _{p1_yes};

Work as f ₁(t) > f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p1=K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) > S _high, when energy-storage system is charged state, k _p1=W ₂× K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) > S _high, when energy-storage system is discharge condition, k _p1=W ₁× K _{p1_yes};

Work as f ₁(t)≤f _ref1, SOC (t) < S _low, when energy-storage system is charged state, k _p1=W ₂× K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) < S _low, when energy-storage system is discharge condition, k _p1=W ₁× K _{p1_no};

Work as f ₁(t)≤f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p1=K _{p1_no};

Work as f ₁(t)≤f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p1=K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) > S _high, when energy-storage system is charged state, k _p1=W ₁× K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) > S _high, when energy-storage system is discharge condition, k _p1=W ₂× K _{p1_no};

K _p10value is carried out according to following setting rule:

Work as f ₁₀(t) > f _ref10, SOC (t) < S _low, when energy-storage system is charged state, k _p10=W ₁× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) < S _low, when energy-storage system is discharge condition, k _p10=W ₂× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p10=K _{p10_yes};

Work as f ₁₀(t) > f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p10=K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) > S _high, when energy-storage system is charged state, k _p10=W ₂× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) > S _high, when energy-storage system is discharge condition, k _p10=W ₁× K _{p10_yes};

Work as f ₁₀(t)≤f _ref10, SOC (t) < S _low, when energy-storage system is charged state, k _p10=W ₂× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) < S _low, when energy-storage system is discharge condition, k _p10=W ₁× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p10=K _{p10_no};

Work as f ₁₀(t)≤f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p10=K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) > S _high, when energy-storage system is charged state, k _p10=W ₁× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) > S _high, when energy-storage system is discharge condition, k _p10=W ₂× K _{p10_no};

K _d1value is carried out according to following setting rule:

Work as f ₁(t)≤f _ref1time, k _d1=0;

Work as f ₁(t) > f _ref1time, k _d1get the first set point;

K _d10value is carried out according to following setting rule:

Work as f ₁₀(t)≤f _ref10time, k _d10=0;

Work as f ₁₀(t) > f _ref10time, k _d10get the second set point;

Wherein, SOC (t) the energy-storage system state-of-charge that is t control cycle; S _lowfor the Low threshold of energy-storage system state-of-charge; S _highfor the high threshold of energy-storage system state-of-charge; K _{p1_yes}for f ₁(t) > f _ref1time proportional control parameter base value; K _{p1_no}for f ₁(t)≤f _ref1time proportional control parameter base value; K _{p10_yes}for f ₁₀(t) > f _ref10time proportional control parameter base value; K _{p10_no}for f ₁₀(t)≤f _ref10time proportional control parameter base value; W ₁, W ₂for proportional control parameter regulating parameter, and W ₁> 1 > W ₂.

In said apparatus, state-of-charge computing module 71 adopts the energy-storage system state-of-charge of following formulae discovery current control period:

SOC(t)＝SOC(t-1)+P _B(t)×Δt；

Wherein, P _bt () goes out force data for the energy-storage system of t control cycle; Δ t is control cycle.

In said apparatus, 1min fluctuation ratio computing module 72 adopts the wind farm grid-connected 1min fluctuation ratio of following formulae discovery:

$f_1\left(t\right)=\frac{\underset{i=\mathrm{0,1,2},...,59}{\max }{P}_{\mathrm{ref}}(t-i)-\underset{i=\mathrm{0,1,2},...,59}{\min }{P}_{\mathrm{ref}}(t-i)}{P_{w\_\mathrm{rate}}};$

In said apparatus, 10min fluctuation ratio computing module 73 adopts the wind farm grid-connected 10min fluctuation ratio of following formulae discovery:

$f_{10}\left(t\right)=\frac{\underset{i=\mathrm{0,1,2},...,599}{\max }{P}_{\mathrm{ref}}(t-i)-\underset{i=\mathrm{0,1,2},...,599}{\min }{P}_{\mathrm{ref}}(t-i)}{P_{w\_\mathrm{rate}}};$

Wherein, P _reft energy storage reference output power that () is t control cycle; P _{w_rate}for wind energy turbine set rated output power.

In said apparatus, discharge and recharge determination module 74 specifically for:

If when the wind power of current control period is greater than the energy storage reference output power of last control cycle, determine that current control period energy-storage system is charged state;

If when the wind power of current control period is less than the energy storage reference output power of last control cycle, determine that current control period energy-storage system is discharge condition.

In said apparatus, first-order low-pass mode block 76 adopts following formula to carry out first-order low-pass ripple to wind power output power, obtains the energy-storage system power output standard of current control period:

$P_B^{\&prime;}\left(t\right)=\frac{T\left(t\right)}{T\left(t\right)+\mathrm{\&Delta;t}}[P_{\mathrm{ref}}(t-1)-P_w\left(t\right)]$

Wherein, P ' _bt energy-storage system power output standard that () is t control cycle; P _wt () is the wind power of t control cycle.

In said apparatus, optimize module 77 concrete according to following rule, amplitude limit carried out to the energy-storage system power output standard of described current control period, the energy-storage system power output standard of the current control period after being optimized:

As P ' _b(t) > P _maxtime, P _{b_cmd}(t)=P _max;

As-P _max≤ P ' _b(t)≤P _maxtime, P _{b_cmd}(t)=P ' _b(t);

As P ' _b(t) <-P _maxtime, P _{b_cmd}(t)=-P _max;

Wherein, P _maxfor the maximum charge-discharge electric power limits value of energy-storage system; P _{b_cmd}(t) for optimize after t control cycle energy-storage system power output standard.

The present invention makes wind-electricity integration fluctuation ratio decline by energy-storage system, meet the requirement that grid company fluctuates to wind-powered electricity generation, improve stability of power system, under the prerequisite ensureing wind power smooth effect, invention increases the energy-storage system life-span, energy storage can be utilized to greatest extent, compared to prior art, the invention solves the contradiction between smooth effect and energy-storage system life-span in wind power smoothing process.

Above-described specific embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only specific embodiments of the invention; the protection range be not intended to limit the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (12)

1., based on the energy-storage system output calculation method adjusting time constant filter in real time, it is characterized in that, comprising:

Go out the energy-storage system state-of-charge of force data, control cycle, last control cycle according to the energy-storage system of current control period, calculate the energy-storage system state-of-charge of current control period;

According to the energy storage reference output power of control cycle each in 1min before wind energy turbine set rated output power, current control period, calculate wind farm grid-connected 1min fluctuation ratio;

According to the energy storage reference output power of control cycle each in 10min before wind energy turbine set rated output power, current control period, calculate wind farm grid-connected 10min fluctuation ratio;

By the wind power of current control period and the energy storage reference output power of last control cycle, determine the charging and discharging state of current control period energy-storage system;

Adopt the following time constant filter regulating formulae discovery current control period:

T(t)＝T(t-1)+a[k _p1×(f ₁(t)-f _ref1)+k _d1×[(f ₁(t)-f ₁(t-1))-(f ₁(t-1)-f ₁(t-2))]]

+b[k _p10×(f ₁₀(t)-f _ref10)+k _d10×[(f ₁₀(t)-f ₁₀(t-1))-(f ₁₀(t-1)-f ₁₀(t-2))]]

Utilize the time constant filter of described current control period, first-order low-pass ripple is carried out to wind power output power, obtain the energy-storage system power output standard of current control period;

Amplitude limit is carried out to the energy-storage system power output standard of described current control period, the energy-storage system power output standard of the current control period after being optimized;

Utilize the energy-storage system power output standard of the current control period after described optimization, instruct energy-storage system to exert oneself;

In described adjustment formula, t is control cycle sequence number; T (t) is the time constant filter of t control cycle; f ₁t () is wind farm grid-connected 1min fluctuation ratio; f ₁₀t () is wind farm grid-connected 10min fluctuation ratio; f _ref1for 1min fluctuation ratio maximum permissible value; f _ref10for 10min fluctuation ratio maximum permissible value; A is 1min fluctuation ratio control weight parameter; B is 10min fluctuation ratio control weight parameter; k _p1for 1min fluctuation ratio proportional control parameter; k _p10for 10min fluctuation ratio proportional control parameter; k _d1for 1min fluctuation ratio differential controling parameters; k _d10for 10min fluctuation ratio differential controling parameters;

Wherein, a, b carry out value according to following setting rule:

Work as f ₁(t) > f _ref1and f ₁₀(t) > f _ref10or f ₁(t)≤f _ref1and f ₁₀(t)≤f _ref10time, a=b=0.5;

Work as f ₁(t) > f _ref1and f ₁₀(t)≤f _ref10time, a=0.7, b=0.3;

Work as f ₁(t)≤f _ref1and f ₁₀(t) > f _ref10time, a=0.3, b=0.7;

K _p1value is carried out according to following setting rule:

Work as f ₁(t) > f _ref1, SOC (t) < S _low, when energy-storage system is charged state, k _p1=W ₁× K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) < S _low, when energy-storage system is discharge condition, k _p1=W ₂× K _{p1_yes};

Work as f ₁(t) > f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p1=K _{p1_yes};

Work as f ₁(t) > f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p1=K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) > S _high, when energy-storage system is charged state, k _p1=W ₂× K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) > S _high, when energy-storage system is discharge condition, k _p1=W ₁× K _{p1_yes};

Work as f ₁(t)≤f _ref1, SOC (t) < S _low, when energy-storage system is charged state, k _p1=W ₂× K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) < S _low, when energy-storage system is discharge condition, k _p1=W ₁× K _{p1_no};

Work as f ₁(t)≤f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p1=K _{p1_no};

Work as f ₁(t)≤f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p1=K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) > S _high, when energy-storage system is charged state, k _p1=W ₁× K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) > S _high, when energy-storage system is discharge condition, k _p1=W ₂× K _{p1_no};

K _p10value is carried out according to following setting rule:

Work as f ₁₀(t) > f _ref10, SOC (t) < S _low, when energy-storage system is charged state, k _p10=W ₁× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) < S _low, when energy-storage system is discharge condition, k _p10=W ₂× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p10=K _{p10_yes};

Work as f ₁₀(t) > f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p10=K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) > S _high, when energy-storage system is charged state, k _p10=W ₂× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) > S _high, when energy-storage system is discharge condition, k _p10=W ₁× K _{p10_yes};

Work as f ₁₀(t)≤f _ref10, SOC (t) < S _low, when energy-storage system is charged state, k _p10=W ₂× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) < S _low, when energy-storage system is discharge condition, k _p10=W ₁× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p10=K _{p10_no};

Work as f ₁₀(t)≤f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p10=K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) > S _high, when energy-storage system is charged state, k _p10=W ₁× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) > S _high, when energy-storage system is discharge condition, k _p10=W ₂× K _{p10_no};

K _d1value is carried out according to following setting rule:

Work as f ₁(t)≤f _ref1time, k _d1=0;

Work as f ₁(t) > f _ref1time, k _d1get the first set point;

K _d10value is carried out according to following setting rule:

Work as f ₁₀(t)≤f _ref10time, k _d10=0;

Work as f ₁₀(t) > f _ref10time, k _d10get the second set point;

Wherein, SOC (t) the energy-storage system state-of-charge that is t control cycle; S _lowfor the Low threshold of energy-storage system state-of-charge; S _highfor the high threshold of energy-storage system state-of-charge; K _{p1_yes}for f ₁(t) > f _ref1time proportional control parameter base value; K _{p1_no}for f ₁(t)≤f _ref1time proportional control parameter base value; K _{p10_yes}for f ₁₀(t) > f _ref10time proportional control parameter base value; K _{p10_no}for f ₁₀(t)≤f _ref10time proportional control parameter base value; W ₁, W ₂for proportional control parameter regulating parameter, and W ₁> 1 > W ₂.

2. method according to claim 1, it is characterized in that, the described energy-storage system state-of-charge going out force data, control cycle, last control cycle according to the energy-storage system of current control period, calculates the energy-storage system state-of-charge of current control period, adopts following formula:

SOC(t)＝SOC(t-1)+P _B(t)×Δt；

Wherein, P _bt () goes out force data for the energy-storage system of t control cycle; Δ t is control cycle.

3. method according to claim 2, is characterized in that,

The described energy storage reference output power according to control cycle each in 1min before wind energy turbine set rated output power, current control period, calculates wind farm grid-connected 1min fluctuation ratio, adopts following formula:

$f_1\left(t\right)=\frac{\underset{i=\mathrm{0,1,2},\&CenterDot;\&CenterDot;\&CenterDot;,59}{\max }{P}_{\mathrm{ref}}(t-i)-\underset{i=\mathrm{0,1,2},\&CenterDot;\&CenterDot;\&CenterDot;,59}{\min }{P}_{\mathrm{ref}}(t-i)}{P_{w\_\mathrm{rate}}};$

The described energy storage reference output power according to control cycle each in 10min before wind energy turbine set rated output power, current control period, calculates wind farm grid-connected 10min fluctuation ratio, adopts following formula:

$f_{10}\left(t\right)=\frac{\underset{i=\mathrm{0,1,2},\&CenterDot;\&CenterDot;\&CenterDot;,599}{\max }{P}_{\mathrm{ref}}(t-i)-\underset{i=\mathrm{0,1,2},\&CenterDot;\&CenterDot;\&CenterDot;,599}{\min }{P}_{\mathrm{ref}}(t-i)}{P_{w\_\mathrm{rate}}};$

Wherein, P _reft energy storage reference output power that () is t control cycle; P _{w_rate}for wind energy turbine set rated output power.

4. method according to claim 1, is characterized in that, described passes through the wind power of current control period and the energy storage reference output power of last control cycle, determines the charging and discharging state of current control period energy-storage system, is specially:

If when the wind power of current control period is greater than the energy storage reference output power of last control cycle, determine that current control period energy-storage system is charged state;

If when the wind power of current control period is less than the energy storage reference output power of last control cycle, determine that current control period energy-storage system is discharge condition.

5. method according to claim 4, it is characterized in that, the described time constant filter utilizing described current control period, carries out first-order low-pass ripple to wind power output power, obtain the energy-storage system power output standard of current control period, adopt following formula:

$P_B^{\&prime;}\left(t\right)=\frac{T\left(t\right)}{T\left(t\right)+\mathrm{\&Delta;t}}[P_{\mathrm{ref}}(t-1)-P_w\left(t\right)];$

Wherein, P ' _bt energy-storage system power output standard that () is t control cycle; P _ref(t-1) be the energy storage reference output power of t-1 control cycle; P _wt () is the wind power of t control cycle.

6. method according to claim 5, is characterized in that, described carries out amplitude limit to the energy-storage system power output standard of described current control period, and the energy-storage system power output standard of the current control period after being optimized, is specially:

As P ' _b(t) > P _maxtime, P _{b_cmd}(t)=P _max;

As-P _max≤ P ' _b(t)≤P _maxtime, P _{b_cmd}(t)=P ' _b(t);

As P ' _b(t) <-P _maxtime, P _{b_cmd}(t)=-P _max;

Wherein, P _maxfor the maximum charge-discharge electric power limits value of energy-storage system; P _{b_cmd}(t) for optimize after t control cycle energy-storage system power output standard.

7., based on the energy-storage system output calculation device adjusting time constant filter in real time, it is characterized in that, comprising:

State-of-charge computing module, for going out the energy-storage system state-of-charge of force data, control cycle, last control cycle according to the energy-storage system of current control period, calculates the energy-storage system state-of-charge of current control period;

1min fluctuation ratio computing module, for the energy storage reference output power according to control cycle each in 1min before wind energy turbine set rated output power, current control period, calculates wind farm grid-connected 1min fluctuation ratio;

10min fluctuation ratio computing module, for the energy storage reference output power according to control cycle each in 10min before wind energy turbine set rated output power, current control period, calculates wind farm grid-connected 10min fluctuation ratio;

Discharge and recharge determination module, for by the wind power of current control period and the energy storage reference output power of last control cycle, determines the charging and discharging state of current control period energy-storage system;

Time constant filter computing module, for adopting the time constant filter of following adjustment formulae discovery current control period:

T(t)＝T(t-1)+a[k _p1×(f ₁(t)-f _ref1)+k _d1×[(f ₁(t)-f ₁(t-1))-(f ₁(t-1)-f ₁(t-2))]]

+b[k _p10×(f ₁₀(t)-f _ref10)+k _d10×[(f ₁₀(t)-f ₁₀(t-1))-(f ₁₀(t-1)-f ₁₀(t-2))]]

First-order low-pass mode block, for utilizing the time constant filter of described current control period, carries out first-order low-pass ripple to wind power output power, obtains the energy-storage system power output standard of current control period;

Optimize module, for carrying out amplitude limit to the energy-storage system power output standard of described current control period, the energy-storage system power output standard of the current control period after being optimized;

Instructing module, for utilizing the energy-storage system power output standard of the current control period after described optimization, instructing energy-storage system to exert oneself;

In described adjustment formula, t is control cycle sequence number; T (t) is the time constant filter of t control cycle; f ₁t () is wind farm grid-connected 1min fluctuation ratio; f ₁₀t () is wind farm grid-connected 10min fluctuation ratio; f _ref1for 1min fluctuation ratio maximum permissible value; f _ref10for 10min fluctuation ratio maximum permissible value; A is 1min fluctuation ratio control weight parameter; B is 10min fluctuation ratio control weight parameter; k _p1for 1min fluctuation ratio proportional control parameter; k _p10for 10min fluctuation ratio proportional control parameter; k _d1for 1min fluctuation ratio differential controling parameters; k _d10for 10min fluctuation ratio differential controling parameters;

Wherein, a, b carry out value according to following setting rule:

Work as f ₁(t) > f _ref1and f ₁₀(t) > f _ref10or f ₁(t)≤f _ref1and f ₁₀(t)≤f _ref10time, a=b=0.5;

Work as f ₁(t) > f _ref1and f ₁₀(t)≤f _ref10time, a=0.7, b=0.3;

Work as f ₁(t)≤f _ref1and f ₁₀(t) > f _ref10time, a=0.3, b=0.7;

K _p1value is carried out according to following setting rule:

Work as f ₁(t) > f _ref1, SOC (t) < S _low, when energy-storage system is charged state, k _p1=W ₁× K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) < S _low, when energy-storage system is discharge condition, k _p1=W ₂× K _{p1_yes};

Work as f ₁(t) > f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p1=K _{p1_yes};

Work as f ₁(t) > f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p1=K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) > S _high, when energy-storage system is charged state, k _p1=W ₂× K _{p1_yes};

Work as f ₁(t) > f _ref1, SOC (t) > S _high, when energy-storage system is discharge condition, k _p1=W ₁× K _{p1_yes};

Work as f ₁(t)≤f _ref1, SOC (t) < S _low, when energy-storage system is charged state, k _p1=W ₂× K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) < S _low, when energy-storage system is discharge condition, k _p1=W ₁× K _{p1_no};

Work as f ₁(t)≤f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p1=K _{p1_no};

Work as f ₁(t)≤f _ref1, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p1=K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) > S _high, when energy-storage system is charged state, k _p1=W ₁× K _{p1_no};

Work as f ₁(t)≤f _ref1, SOC (t) > S _high, when energy-storage system is discharge condition, k _p1=W ₂× K _{p1_no};

K _p10value is carried out according to following setting rule:

Work as f ₁₀(t) > f _ref10, SOC (t) < S _low, when energy-storage system is charged state, k _p10=W ₁× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) < S _low, when energy-storage system is discharge condition, k _p10=W ₂× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p10=K _{p10_yes};

Work as f ₁₀(t) > f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p10=K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) > S _high, when energy-storage system is charged state, k _p10=W ₂× K _{p10_yes};

Work as f ₁₀(t) > f _ref10, SOC (t) > S _high, when energy-storage system is discharge condition, k _p10=W ₁× K _{p10_yes};

Work as f ₁₀(t)≤f _ref10, SOC (t) < S _low, when energy-storage system is charged state, k _p10=W ₂× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) < S _low, when energy-storage system is discharge condition, k _p10=W ₁× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is charged state, k _p10=K _{p10_no};

Work as f ₁₀(t)≤f _ref10, S _low≤ SOC (t)≤S _high, when energy-storage system is discharge condition, k _p10=K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) > S _high, when energy-storage system is charged state, k _p10=W ₁× K _{p10_no};

Work as f ₁₀(t)≤f _ref10, SOC (t) > S _high, when energy-storage system is discharge condition, k _p10=W ₂× K _{p10_no};

K _d1value is carried out according to following setting rule:

Work as f ₁(t)≤f _ref1time, k _d1=0;

Work as f ₁(t) > f _ref1time, k _d1get the first set point;

K _d10value is carried out according to following setting rule:

Work as f ₁₀(t)≤f _ref10time, k _d10=0;

Work as f ₁₀(t) > f _ref10time, k _d10get the second set point;

Wherein, SOC (t) the energy-storage system state-of-charge that is t control cycle; S _lowfor the Low threshold of energy-storage system state-of-charge; S _highfor the high threshold of energy-storage system state-of-charge; K _{p1_yes}for f ₁(t) > f _ref1time proportional control parameter base value; K _{p1_no}for f ₁(t)≤f _ref1time proportional control parameter base value; K _{p10_yes}for f ₁₀(t) > f _ref10time proportional control parameter base value; K _{p10_no}for f ₁₀(t)≤f _ref10time proportional control parameter base value; W ₁, W ₂for proportional control parameter regulating parameter, and W ₁> 1 > W ₂.

8. device according to claim 7, is characterized in that, described state-of-charge computing module adopts the energy-storage system state-of-charge of following formulae discovery current control period:

SOC(t)＝SOC(t-1)+P _B(t)×Δt；

Wherein, P _bt () goes out force data for the energy-storage system of t control cycle; Δ t is control cycle.

9. device according to claim 8, is characterized in that,

Described 1min fluctuation ratio computing module adopts the wind farm grid-connected 1min fluctuation ratio of following formulae discovery:

$f_1\left(t\right)=\frac{\underset{i=\mathrm{0,1,2},\&CenterDot;\&CenterDot;\&CenterDot;,59}{\max }{P}_{\mathrm{ref}}(t-i)-\underset{i=\mathrm{0,1,2},\&CenterDot;\&CenterDot;\&CenterDot;,59}{\min }{P}_{\mathrm{ref}}(t-i)}{P_{w\_\mathrm{rate}}};$

Described 10min fluctuation ratio computing module adopts the wind farm grid-connected 10min fluctuation ratio of following formulae discovery:

$f_{10}\left(t\right)=\frac{\underset{i=\mathrm{0,1,2},\&CenterDot;\&CenterDot;\&CenterDot;,599}{\max }{P}_{\mathrm{ref}}(t-i)-\underset{i=\mathrm{0,1,2},\&CenterDot;\&CenterDot;\&CenterDot;,599}{\min }{P}_{\mathrm{ref}}(t-i)}{P_{w\_\mathrm{rate}}};$

Wherein, P _reft energy storage reference output power that () is t control cycle; P _{w_rate}for wind energy turbine set rated output power.

10. device according to claim 7, is characterized in that, described discharge and recharge determination module specifically for:

If when the wind power of current control period is greater than the energy storage reference output power of last control cycle, determine that current control period energy-storage system is charged state;

If when the wind power of current control period is less than the energy storage reference output power of last control cycle, determine that current control period energy-storage system is discharge condition.

11. devices according to claim 10, is characterized in that, described first-order low-pass mode block adopts following formula to carry out first-order low-pass ripple to wind power output power, obtain the energy-storage system power output standard of current control period:

$P_B^{\&prime;}\left(t\right)=\frac{T\left(t\right)}{T\left(t\right)+\mathrm{\&Delta;t}}[P_{\mathrm{ref}}(t-1)-P_w\left(t\right)]$

Wherein, P ' _bt energy-storage system power output standard that () is t control cycle; P _ref(t-1) be the energy storage reference output power of t-1 control cycle; P _wt () is the wind power of t control cycle.

12. devices according to claim 11, it is characterized in that, described optimization module is concrete according to following rule, carries out amplitude limit, the energy-storage system power output standard of the current control period after being optimized to the energy-storage system power output standard of described current control period:

As P ' _b(t) > P _maxtime, P _{b_cmd}(t)=P _max;

As-P _max≤ P ' _b(t)≤P _maxtime, P _{b_cmd}(t)=P ' _b(t);

As P ' _b(t) <-P _maxtime, P _{b_cmd}(t)=-P _max;

Wherein, P _maxfor the maximum charge-discharge electric power limits value of energy-storage system; P _{b_cmd}(t) for optimize after t control cycle energy-storage system power output standard.