CN102355006A - Joint wind-solar-storage joint grid-connected power generation coordination control method - Google Patents

Abstract

The invention provides a wind-solar-storage joint grid-connected power generation coordination control method. The control method is characterized by comprising the following steps: establishing a wind energy conversion unit mathematical model which comprises a wind energy conversion unit aerodynamic mathematical model, a two-block shafting mathematical model equation for a wind turbine and a generator and a voltage equation for a doubly-fed induction motor under a synchronous rotation coordinate system; establishing a solar energy PV power generation unit mathematical model; establishing an energy storage unit mathematical model; and carrying out joint grid-connected power generation coordination control including active power coordination control and reactive power coordination control. In the invention, on the basis of scientific and reasonable establishment and characteristics analysis of the wind energy conversion unit mathematical model, the solar energy PV power generation unit mathematical model and the energy storage unit mathematical model, space-time complementarity of wind and solar resource is fully considered so as to improve the capability of a power grid in new energy absorption. Therefore, the wind-solar-storage joint grid-connected power generation coordination control method provided by the invention has the advantages of strong adaptability, high actual application value, capability of efficiently utilizing renewable energy and the like.

Description

The control method for coordinating that generates electricity by way of merging two or more grid systems is united in a kind of scene storage

Technical field

The present invention relates to power system operation and control, is that the control method for coordinating that generates electricity by way of merging two or more grid systems is united in a kind of scene storage, is applied to honourable storing cogeneration modeling and simulating, the specificity analysis that is incorporated into the power networks, studies and the optimal economic capacity ratio with the electrical network interaction mechanism.

Background technology

In recent years, along with socioeconomic fast development, energy demand is growing.The carbon emission that a large amount of fossil fuels produce causes global warming, severe contamination human environment of depending on for existence.Adopt environmentally friendly new forms of energy and clean reproducible energy to replace traditional fossil fuel, people-oriented, walk continuable scientific development strategic thought to meet country's proposition.Enriching, distribute extensively, be easy to exploitation in China's wind energy and solar energy resources, is the ideal green alternative energy source.Therefore, scene storage being united to generate electricity by way of merging two or more grid systems and coordinates control and have important and practical meanings to the efficient development and use of wind energy and solar energy resources and with the research of electrical network interaction mechanism.

The control method for coordinating that generates electricity by way of merging two or more grid systems is united in traditional scene storage: be that honourable storage unit compiles the energy through direct current (DC) bus that public inverter connects on the one hand, its shortcoming is the system that is not suitable for scale Expansion Planning at a specified future date; Be to coordinate control through regulating inverter direct-flow side voltage and power on the other hand, its shortcoming is to be unfavorable for that extensive scene stores up the sporadic development utilization, has limited the efficient utilization of regenerative resource.

Summary of the invention

The objective of the invention is; Providing a kind of is analyzing on wind energy converting unit, solar energy power generating unit and the energy-storage units Mathematical Modeling characteristic basis; Consider that honourable resource has the space-time complementarity; Improve electrical network to the new forms of energy ability to arrange jobs, adaptability is strong, and the control method for coordinating that generates electricity by way of merging two or more grid systems is united in the scene storage with higher actual application value.

The objective of the invention is to realize by following technical scheme:

The control method for coordinating that generates electricity by way of merging two or more grid systems is united in a kind of scene storage, it is characterized in that it may further comprise the steps:

1) sets up wind energy converting unit Mathematical Modeling

(1.1) wind energy converting unit aerodynamics Mathematical Modeling is:

$P_M={\mathrm{\ρ}}_{\mathrm{air}}{C}_p(\mathrm{\λ},\mathrm{\β})\mathrm{\π}{R}^2{V}_w^3/2---\left(1\right)$

Wherein: P _MThe wind-powered electricity generation unit mechanical output that the energy of from wind, catching for the wind energy converting unit changes into,

ρ _AirBe atmospheric density,

R is the wind turbine impeller radius,

λ is a tip speed ratio,

β is a propeller pitch angle,

C _pFor the wind energy conversion efficiency coefficient of blade, be the function of λ and β,

V _wBe wind speed;

(1.2) two matter piece axle coefficients of wind turbine and generator model equation is:

$\left\{\begin{array}{c}{2H}_{T}d{\mathrm{\ω}}_T/\mathrm{dt}=T_M-K_S{\mathrm{\θ}}_S-D_T{\mathrm{\ω}}_T\\ {2H}_{G}d{\mathrm{\ω}}_G/\mathrm{dt}=K_S{\mathrm{\θ}}_S-T_E-D_G{\mathrm{\ω}}_G\\ {\mathrm{d\θ}}_S/\mathrm{dt}={\mathrm{\ω}}_0({\mathrm{\ω}}_T-{\mathrm{\ω}}_G)\end{array}\right.---\left(2\right)$

Wherein: H _T, H _GBe respectively wind turbine and generator inertia constant,

K _SBe the stiffness coefficient of axle,

D _T, D _GBe respectively the damping coefficient of wind turbine rotor and generator amature,

θ _SBe relative angular displacement between the two matter pieces,

T _M, T _EBe respectively wind turbine machine torque and generator electromagnetic torque,

ω _T, ω _GBe respectively wind turbine and generator amature rotating speed,

ω ₀Be synchronous speed;

(1.3) voltage equation of doubly fed induction generator does under the synchronous rotating frame

$\left\{\begin{array}{c}{u}_{\mathrm{sd}}=d{\mathrm{\ψ}}_{\mathrm{sd}}/\mathrm{dt}-{\mathrm{\ω}}_s{\mathrm{\ψ}}_{\mathrm{sq}}+R_s{i}_{\mathrm{sd}}\\ u_{\mathrm{sq}}=d{\mathrm{\ψ}}_{\mathrm{sq}}/\mathrm{dt}-{\mathrm{\ω}}_s{\mathrm{\ψ}}_{\mathrm{sd}}+R_s{i}_{\mathrm{sq}}\\ u_{\mathrm{rd}}=d{\mathrm{\ψ}}_{\mathrm{rd}}/\mathrm{dt}-{\mathrm{\ω}}_s{\mathrm{\ψ}}_{\mathrm{rq}}+R_r{i}_{\mathrm{rd}}\\ u_{\mathrm{rq}}=d{\mathrm{\ψ}}_{\mathrm{rq}}/\mathrm{dt}-{\mathrm{\ω}}_s{\mathrm{\ψ}}_{\mathrm{rd}}+R_r{i}_{\mathrm{rq}}\end{array}\right.---\left(3\right)$

The magnetic linkage equation does

$\left\{\begin{array}{c}{\mathrm{\ψ}}_{\mathrm{sd}}=L_s{i}_{\mathrm{sd}}+L_m{i}_{\mathrm{rd}}\\ {\mathrm{\ψ}}_{\mathrm{sq}}=L_s{i}_{\mathrm{sq}}+L_m{i}_{\mathrm{rq}}\\ {\mathrm{\ψ}}_{\mathrm{rd}}=L_r{i}_{\mathrm{rd}}+L_m{i}_{\mathrm{sd}}\\ {\mathrm{\ψ}}_{\mathrm{rq}}=L_r{i}_{\mathrm{rq}}+L_m{i}_{\mathrm{sq}}\end{array}\right.---\left(4\right)$

Wherein: ω _sBe coordinate system rotation angular speed,

U, i, ψ are respectively voltage, electric current and the magnetic linkage of winding,

R is the resistance of winding,

L _s, L _rBe respectively the self-induction of stator winding and rotor winding,

L _mBe the mutual inductance between the stator and rotor winding,

Subscript s, r represent the stator and the rotor of motor respectively,

Subscript d, q represent d, the q axle winding of motor respectively,

S is the revolutional slip of motor;

2) set up solar energy power generating unit Mathematical Modeling

Under the reference conditions, consider that the photovoltaic array V-I equation of intensity of solar radiation and variation of ambient temperature does

$I=N_p^{\mathrm{PV}}{I}_{\mathrm{sc}}\·\{1-C_1[\exp \left(\frac{V-\mathrm{dV}}{C_2{N}_s^{\mathrm{PV}}{V}_{\mathrm{oc}}}\right)-1\left]\right\}+\mathrm{dI}---\left(5\right)$

In the formula: $\left\{\begin{array}{c}{\mathrm{<figure-callout\; id="1"\; label="C"\; filenames="HDA0000095634810000011.png"\; state="\{\{state\}\}">C</figure-callout>}}_1=(1-I_m/I_{\mathrm{sc}})\&CenterDot;\exp ({-V}_m/C_2{V}_{\mathrm{oc}})\\ {\mathrm{<figure-callout\; id="2"\; label="C"\; filenames="HDA0000095634810000011.png"\; state="\{\{state\}\}">C</figure-callout>}}_2=(V_m/V_{\mathrm{oc}}-1)/\left(\ln (1-I_m/I_{\mathrm{sc}})\right)\\ \mathrm{dI}=-\mathrm{\&alpha;}\&CenterDot;G/G_{\mathrm{ref}}\&CenterDot;(T_c-T_{\mathrm{ref}})+(G/G_{\mathrm{ref}}-1)\&CenterDot;N_p^{\mathrm{PV}}{I}_{\mathrm{sc}}\\ \mathrm{dV}=\mathrm{\&beta;}\&CenterDot;\mathrm{dT}-R_s\&CenterDot;\mathrm{dI}\end{array}\right.$

Wherein: V, I are respectively photovoltaic array terminal voltage and end electric current,

V _Oc, I _ScBe respectively open circuit voltage and short circuit current that photovoltaic module producer provides,

V _m, I _mBe respectively maximum power point voltage and electric current,

α, β are respectively under the reference radiation intensity, current temperature variation coefficient and voltage temperature variation coefficient,

R _sBe the series resistance of assembly,

is respectively the parallelly connected number and the serial number of assembly in the photovoltaic array;

3) set up the energy-storage units Mathematical Modeling

The batteries Mathematical Modeling does

$V_B=N_{B-s}{E}_0-\frac{\mathrm{KSOC}}{\mathrm{SOC}-N_{B-s}{Q}_n{\&Integral;}_0^t{i}_B\left(\mathrm{\&tau;}\right)\mathrm{d\&tau;}}+\mathrm{Aexp}(-B{\&Integral;}_0^t{i}_B\left(\mathrm{\&tau;}\right)\mathrm{d\&tau;}+C_P)-\frac{{\mathrm{Ri}}_B}{N_{B-s}}---\left(6\right)$

In the formula: $\left\{\begin{array}{c}{C}_P=C_t(T_b-25)\\ C_R=1-0.025(T_b-25)\\ R=N_{B-s}{V}_n(1-\mathrm{\&eta;})C_R/\left(0.2N_{B-p}{Q}_n\right)\\ \mathrm{SOC}=\frac{N_{B-s}{N}_{B-p}{Q}_n-{\&Integral;}_0^t{i}_B\left(\mathrm{\&tau;}\right)\mathrm{d\&tau;}}{N_{B-s}{N}_{B-p}{Q}_n}\&times;100\%\end{array}\right.$

Wherein: SOC is the state-of-charge of batteries,

A, B are respectively change in voltage coefficient and volume change coefficient,

K is the polarizing voltage constant,

C _pBe polarization reaction temperature compensation factor,

C _RBe the resistance temperature compensation factor,

T _bBe battery temp,

C _tBe the polarization reaction constant,

Q _nBe battery rating,

V _nBe the battery rated voltage,

E ₀Be the storage battery initial potential,

V is an accumulator voltage,

η is a battery efficiency,

N _B-sWith N _B-pBe respectively storage battery series connection and parallelly connected number;

The fuel cell Mathematical Modeling does

$V_{\mathrm{FC}}=N_{\mathrm{cell}}[E-\ln \left(\frac{I_{\mathrm{FC}}+A_{\mathrm{cell}}{i}_n}{A_{\mathrm{cell}}{i}_o}\right)-r\left(\frac{I_{\mathrm{FC}}}{A_{\mathrm{cell}}}\right)-\mathrm{mexp}\left(n\frac{I_{\mathrm{FC}}}{A_{\mathrm{cell}}}\right)]---\left(7\right)$

Wherein: V _FCBe the fuel cell terminal voltage,

I _FCBe fuel cell end electric current,

N _CellBe fuel cell series quantity,

A _CellBe the fuel cell electrode area,

R is a sheet resistance,

i _nBe the internal operation current density,

i _oBe inner nominal current density,

M and n are the junction loss constants;

4) unite the coordination control of generating electricity by way of merging two or more grid systems

(4.1) active power is coordinated control

Meritorious expression formula between wind light generation unit and the network load demand does

P _net＝P _wind+P _PV-P _load-P _sc (8)

Wherein: P _NetBe the active power amount of unbalance,

P _WindExert oneself for the wind energy converting unit is meritorious,

P _PVExert oneself for solar power generation unit is meritorious,

P _LoadBe the burden with power demand of dispatching of power netwoks appointment,

P _ScFor scene storage association system from electricity consumption, less because of accounting for the generating ratio, can ignore;

The superfluous P of wind light generation unit active power _Net＞0 o'clock, the power balance equation formula did

P _wind+P _PV＝P _load+P _ES-i _n (9)

Wherein: P _ES-inFor energy-storage units absorbs active power, promptly storage battery absorbs meritorious power P _Battery-in

Wind light generation cell power vacancy P _Net＜0 o'clock, the power balance equation formula did

P _wind+P _PV+P _ES-out＝P _load (10)

In the formula: P _ES-out=P _{Battery_out}+ P _FC

Wherein: P _ES-outThe active power of sending for energy-storage units,

P _{Battery_out}The active power of sending for storage battery,

P _FCThe active power of sending for fuel cell, when storage battery state-of-charge η＜25%, the starting fluid battery unit;

(4.2) reactive power is coordinated control

Wind light generation cell power factor range of operation is 0.85 to lag behind 0.85 in advance, is 0.98 to lag behind when normally moving, and the meritorious expression formula between wind light generation unit and the network load demand does

Q _net＝Q _Wind-out+Q _PV-out-Q _load (11)

Wherein: Q _NetBe the reactive power amount of unbalance,

Q _{Wind_out}Exert oneself for the wind energy converting unit is idle,

Q _{PV_out}Exert oneself for solar power generation unit is idle,

Q _LoadLoad or burden without work demand for the dispatching of power netwoks appointment;

The superfluous Q of wind light generation unit reactive power _Net＞0 and the point voltage that is incorporated into the power networks be higher than its rated voltage V＞V _RefThe time, the power balance equation formula does

Q _Wind-out+Q _PV-out＝Q _load+Q _ES-in+Q _STATCOM-in (12)

Wherein: Q _ES-inBe the reactive power of energy-storage units absorption,

Q _STATCOM-inReactive power for STATCOM STATCOM absorption;

Wind light generation cell power vacancy Q _Net＜0 o'clock, the power balance equation formula did

Q _load＝Q _Wind-out+Q _PV-out+Q _ES-out+Q _STATCOM-out V＜V _ref (13)

Or

Q _load+Q _STATCOM-in＝Q _Wind-out+Q _PV-out+Q _ES-out V＞V _ref (14)

Wherein: Q _ES-outThe reactive power of sending for energy-storage units,

Q _STATCOM-outThe reactive power of sending for STATCOM STATCOM;

When honourable unit runs on power factor 0.98 when leading, i.e. Q _Net＜0, the reactive power equilibrium equation

Q _Wind-in+Q _PV-in+Q _load＝Q _ES-out+Q _STATCOM-out V＜V _ref (15)

Or

Q _Wind-in+Q _PV-in+Q _load+Q _STATCOM-in＝Q _ES-out V＞V _ref (16)

The control method for coordinating that generates electricity by way of merging two or more grid systems is united in a kind of scene storage of the present invention; Analyze in scientific and reasonable foundation on the basis of wind energy converting unit, solar energy power generating unit and energy-storage units Mathematical Modeling characteristic; It is complementary to have taken into full account the space-time that honourable resource had, and then has improved electrical network to the new forms of energy ability to arrange jobs, and it is strong to have adaptability; Actual application value is high, can make the renewable energy resources advantage such as efficiently utilized.

Description of drawings

Fig. 1 is based on the sample calculation analysis single line schematic diagram of IEEE 3 machine 1-9 node test systems.

Fig. 2 is somewhere typical case's day actual daily load demand curve sketch map.

Fig. 3 is the wind speed sketch map of actual monitoring.

Fig. 4 is the intensity of solar radiation sketch map of actual acquisition.

Fig. 5 is a wind energy converting unit power power curve sketch map.

Fig. 6 is a solar energy power generating cell power power curve sketch map.

Fig. 7 is meritorious, the reactive power sketch map that batteries absorbs and sends.

Fig. 8 is that cell of fuel cell is meritorious, the idle sketch map of exerting oneself.

Embodiment

The control method for coordinating that generates electricity by way of merging two or more grid systems is united in a kind of scene storage of the present invention, may further comprise the steps:

1) sets up wind energy converting unit Mathematical Modeling

(1.1) wind energy converting unit aerodynamics Mathematical Modeling is:

$P_M={\mathrm{\&rho;}}_{\mathrm{air}}{C}_p(\mathrm{\&lambda;},\mathrm{\&beta;})\mathrm{\&pi;}{R}^2{\mathrm{<figure-callout\; id="3"\; label="V\; w"\; filenames="HDA0000095634810000011.png"\; state="\{\{state\}\}">V</figure-callout>}}_w^{}/2---\left(1\right)$

Wherein: P _MThe wind-powered electricity generation unit mechanical output that the energy of from wind, catching for the wind energy converting unit changes into, ρ _AirBe atmospheric density,

R is the wind turbine impeller radius,

λ is a tip speed ratio,

β is a propeller pitch angle,

C _pFor the wind energy conversion efficiency coefficient of blade, be the function of λ and β,

V _wBe wind speed;

(1.2) two matter piece axle coefficients of wind turbine and generator model equation is:

$\left\{\begin{array}{c}{2H}_{T}d{\mathrm{\&omega;}}_T/\mathrm{dt}=T_M-K_S{\mathrm{\&theta;}}_S-D_T{\mathrm{\&omega;}}_T\\ {2H}_{G}d{\mathrm{\&omega;}}_G/\mathrm{dt}=K_S{\mathrm{\&theta;}}_S-T_E-D_G{\mathrm{\&omega;}}_G\\ {\mathrm{d\&theta;}}_S/\mathrm{dt}={\mathrm{\&omega;}}_0({\mathrm{\&omega;}}_T-{\mathrm{\&omega;}}_G)\end{array}\right.---\left(2\right)$

Wherein: H _T, H _GBe respectively wind turbine and generator inertia constant,

K _SBe the stiffness coefficient of axle,

D _T, D _GBe respectively the damping coefficient of wind turbine rotor and generator amature,

θ _SBe relative angular displacement between the two matter pieces,

T _M, T _EBe respectively wind turbine machine torque and generator electromagnetic torque,

ω _T, ω _GBe respectively wind turbine and generator amature rotating speed,

ω ₀Be synchronous speed;

(1.3) voltage equation of doubly fed induction generator does under the synchronous rotating frame

$\left\{\begin{array}{c}{u}_{\mathrm{sd}}=d{\mathrm{\&psi;}}_{\mathrm{sd}}/\mathrm{dt}-{\mathrm{\&omega;}}_s{\mathrm{\&psi;}}_{\mathrm{sq}}+R_s{i}_{\mathrm{sd}}\\ u_{\mathrm{sq}}=d{\mathrm{\&psi;}}_{\mathrm{sq}}/\mathrm{dt}-{\mathrm{\&omega;}}_s{\mathrm{\&psi;}}_{\mathrm{sd}}+R_s{i}_{\mathrm{sq}}\\ u_{\mathrm{rd}}=d{\mathrm{\&psi;}}_{\mathrm{rd}}/\mathrm{dt}-{\mathrm{\&omega;}}_s{\mathrm{\&psi;}}_{\mathrm{rq}}+R_r{i}_{\mathrm{rd}}\\ u_{\mathrm{rq}}=d{\mathrm{\&psi;}}_{\mathrm{rq}}/\mathrm{dt}-{\mathrm{\&omega;}}_s{\mathrm{\&psi;}}_{\mathrm{rd}}+R_r{i}_{\mathrm{rq}}\end{array}\right.---\left(3\right)$

The magnetic linkage equation does

$\left\{\begin{array}{c}{\mathrm{\&psi;}}_{\mathrm{sd}}=L_s{i}_{\mathrm{sd}}+L_m{i}_{\mathrm{rd}}\\ {\mathrm{\&psi;}}_{\mathrm{sq}}=L_s{i}_{\mathrm{sq}}+L_m{i}_{\mathrm{rq}}\\ {\mathrm{\&psi;}}_{\mathrm{rd}}=L_r{i}_{\mathrm{rd}}+L_m{i}_{\mathrm{sd}}\\ {\mathrm{\&psi;}}_{\mathrm{rq}}=L_r{i}_{\mathrm{rq}}+L_m{i}_{\mathrm{sq}}\end{array}\right.---\left(4\right)$

Wherein: ω _sBe coordinate system rotation angular speed,

U, i, ψ are respectively voltage, electric current and the magnetic linkage of winding,

R is the resistance of winding,

L _s, L _rBe respectively the self-induction of stator winding and rotor winding,

L _mBe the mutual inductance between the stator and rotor winding,

Subscript s, r represent the stator and the rotor of motor respectively,

Subscript d, q represent d, the q axle winding of motor respectively,

S is the revolutional slip of motor;

2) set up solar energy power generating unit Mathematical Modeling

Under the reference conditions, consider that the photovoltaic array V-I equation of intensity of solar radiation and variation of ambient temperature does

$I=N_p^{\mathrm{PV}}{I}_{\mathrm{sc}}\&CenterDot;\{1-C_1[\exp \left(\frac{V-\mathrm{dV}}{C_2{N}_s^{\mathrm{PV}}{V}_{\mathrm{oc}}}\right)-1\left]\right\}+\mathrm{dI}---\left(5\right)$

In the formula: $\left\{\begin{array}{c}{\mathrm{<figure-callout\; id="1"\; label="C"\; filenames="HDA0000095634810000011.png"\; state="\{\{state\}\}">C</figure-callout>}}_1=(1-I_m/I_{\mathrm{sc}})\&CenterDot;\exp ({-V}_m/C_2{V}_{\mathrm{oc}})\\ {\mathrm{<figure-callout\; id="2"\; label="C"\; filenames="HDA0000095634810000011.png"\; state="\{\{state\}\}">C</figure-callout>}}_2=(V_m/V_{\mathrm{oc}}-1)/\left(\ln (1-I_m/I_{\mathrm{sc}})\right)\\ \mathrm{dI}=-\mathrm{\&alpha;}\&CenterDot;G/G_{\mathrm{ref}}\&CenterDot;(T_c-T_{\mathrm{ref}})+(G/G_{\mathrm{ref}}-1)\&CenterDot;N_p^{\mathrm{PV}}{I}_{\mathrm{sc}}\\ \mathrm{dV}=\mathrm{\&beta;}\&CenterDot;\mathrm{dT}-R_s\&CenterDot;\mathrm{dI}\end{array}\right.$

Wherein: V, I are respectively photovoltaic array terminal voltage and end electric current,

V _Oc, I _ScBe respectively open circuit voltage and short circuit current that photovoltaic module producer provides,

V _m, I _mBe respectively maximum power point voltage and electric current,

α, β are respectively under the reference radiation intensity, current temperature variation coefficient and voltage temperature variation coefficient,

R _sBe the series resistance of assembly,

is respectively the parallelly connected number and the serial number of assembly in the photovoltaic array;

3) set up the energy-storage units Mathematical Modeling

The batteries Mathematical Modeling does

$V_B=N_{B-s}{E}_0-\frac{\mathrm{KSOC}}{\mathrm{SOC}-N_{B-s}{Q}_n{\&Integral;}_0^t{i}_B\left(\mathrm{\&tau;}\right)\mathrm{d\&tau;}}+\mathrm{Aexp}(-B{\&Integral;}_0^t{i}_B\left(\mathrm{\&tau;}\right)\mathrm{d\&tau;}+C_P)-\frac{{\mathrm{Ri}}_B}{N_{B-s}}---\left(6\right)$

In the formula: $\left\{\begin{array}{c}{C}_P=C_t(T_b-25)\\ C_R=1-0.025(T_b-25)\\ R=N_{B-s}{V}_n(1-\mathrm{\&eta;})C_R/\left(0.2N_{B-p}{Q}_n\right)\\ \mathrm{SOC}=\frac{N_{B-s}{N}_{B-p}{Q}_n-{\&Integral;}_0^t{i}_B\left(\mathrm{\&tau;}\right)\mathrm{d\&tau;}}{N_{B-s}{N}_{B-p}{Q}_n}\&times;100\%\end{array}\right.$

Wherein: SOC is the state-of-charge of batteries,

B, B are respectively change in voltage coefficient and volume change coefficient,

K is the polarizing voltage constant,

C _pBe polarization reaction temperature compensation factor,

C _RBe the resistance temperature compensation factor,

T _bBe battery temp,

C _tBe the polarization reaction constant,

Q _nBe battery rating,

V _nBe the battery rated voltage,

E ₀Be the storage battery initial potential,

V is an accumulator voltage,

η is a battery efficiency,

N _B-sWith N _B-pBe respectively storage battery series connection and parallelly connected number;

The fuel cell Mathematical Modeling does

$V_{\mathrm{FC}}=N_{\mathrm{cell}}[E-\ln \left(\frac{I_{\mathrm{FC}}+A_{\mathrm{cell}}{i}_n}{A_{\mathrm{cell}}{i}_o}\right)-r\left(\frac{I_{\mathrm{FC}}}{A_{\mathrm{cell}}}\right)-\mathrm{mexp}\left(n\frac{I_{\mathrm{FC}}}{A_{\mathrm{cell}}}\right)]---\left(7\right)$

Wherein: V _FCBe the fuel cell terminal voltage,

I _FCBe fuel cell end electric current,

N _CellBe fuel cell series quantity,

A _CellBe the fuel cell electrode area,

R is a sheet resistance,

i _nBe the internal operation current density,

i _oBe inner nominal current density,

M and n are the junction loss constants;

4) unite the coordination control of generating electricity by way of merging two or more grid systems

(4.1) active power is coordinated control

Meritorious expression formula between wind light generation unit and the network load demand does

P _net＝P _wind+P _PV-P _load-P _sc (8)

Wherein: P _NetBe the active power amount of unbalance,

P _WindExert oneself for the wind energy converting unit is meritorious,

P _PVExert oneself for solar power generation unit is meritorious,

P _LoadBe the burden with power demand of dispatching of power netwoks appointment,

P _ScFor scene storage association system from electricity consumption, less because of accounting for the generating ratio, can ignore;

The superfluous P of wind light generation unit active power _Net＞0 o'clock, the power balance equation formula did

P _wind+P _PV＝P _load+P _ES-in (9)

Wherein: P _ES-inFor energy-storage units absorbs active power, promptly storage battery absorbs meritorious power P _Battery-in

Wind light generation cell power vacancy P _Net＜0 o'clock, the power balance equation formula did

P _wind+P _PV+P _ES-out＝P _load (10)

In the formula: P _ES-out=P _{Battery_out}+ P _FC

Wherein: P _ES-outThe active power of sending for energy-storage units,

P _{Battery_out}The active power of sending for storage battery,

P _FCThe active power of sending for fuel cell, when storage battery state-of-charge η＜25%, the starting fluid battery unit;

(4.2) reactive power is coordinated control

Wind light generation cell power factor range of operation is 0.85 to lag behind 0.85 in advance, is 0.98 to lag behind when normally moving, and the meritorious expression formula between wind light generation unit and the network load demand does

Q _net＝Q _Wind-out+Q _PV-out-Q _load (11)

Wherein: Q _NetBe the reactive power amount of unbalance,

Q _{Wind_out}Exert oneself for the wind energy converting unit is idle,

Q _{PV_out}Exert oneself for solar power generation unit is idle,

Q _LoadLoad or burden without work demand for the dispatching of power netwoks appointment;

The superfluous Q of wind light generation unit reactive power _Net＞0 and the point voltage that is incorporated into the power networks be higher than its rated voltage V＞V _RefThe time, the power balance equation formula does

Q _Wind-out+Q _PV-out＝Q _load+Q _ES-in+Q _STATCOM-in (12)

Wherein: Q _ES-inBe the reactive power of energy-storage units absorption,

Q _STATCOM-inReactive power for STATCOM STATCOM absorption;

Wind light generation cell power vacancy Q _Net＜0 o'clock, the power balance equation formula did

Q _load＝Q _Wind-out+Q _PV-out+Q _ES-out+Q _STATCOM-out V＜V _ref (13)

Or

Q _load+Q _STATCOM-in＝Q _Wind-out+Q _PV-out+Q _ES-out V＞V _ref (14)

Wherein: Q _ES-outThe reactive power of sending for energy-storage units,

Q _STATCOM-outThe reactive power of sending for STATCOM STATCOM;

When honourable unit runs on power factor 0.98 when leading, i.e. Q _Net＜0, the reactive power equilibrium equation

Q _Wind-in+Q _PV-in+Q _load＝Q _ES-out+Q _STATCOM-out V＜V _ref (15)

Or

Q _Wind-in+Q _PV-in+Q _load+Q _STATCOM-in＝Q _ES-out V＞V _ref (16)

Instantiation:

Be the basis with somewhere actual measurement load curve and Practical Meteorological Requirements condition, honourable storing cogeneration control method for coordinating is analyzed.Fig. 1 is the sample calculation analysis single line schematic diagram based on IEEE 3 machine 1-9 node test systems; Wind energy converting unit, solar energy power generating unit, energy-storage units exchange collection bus with STATCOM (STATCOM) through 20kV electric energy are pooled together, and incorporate high-voltage fence through step-up transformer at 6 nodes.Fig. 2 is typical case's day actual load demand curve in somewhere, northeast summer, summer typical case's day wind speed data and intensity of solar radiation data such as Fig. 3 and shown in Figure 4.On this basis honourable storing cogeneration control method for coordinating is carried out simulation analysis, Fig. 5 is under wind friction velocity shown in Figure 3, meritorious, the idle situation of exerting oneself (power factor is 0.98, lags behind) of wind energy converting unit.Can know that by figure the wind energy converting unit is exerted oneself consistent with the wind speed variation tendency, and under the conciliation of propeller pitch angle controller, realize its maximum power tracing.Fig. 6 is under solar radiation condition shown in Figure 4, meritorious, the idle change curve (power factor is 0.98, and is leading) of solar energy power generating unit.Analysis shows that intensity of solar radiation is the major influence factors of solar energy power generating cell power output, and the photovoltaic cell ambient temperature is less to its influence, and the photovoltaic power curve is almost identical with the solar radiation variation trends.Fig. 7 and Fig. 8 are respectively storage battery and the fuel battery power change curve on honourable storing cogeneration control method basis.Can know that by Fig. 7 when control method for coordinating required balance power during greater than the energy-storage units capacity, energy-storage units is exerted oneself by its maximum capacity plan that allows, and when the batteries state-of-charge was lower than certain value 25%, the starting fluid battery unit substituted storage battery.Fig. 8 has shown that the batteries state-of-charge is lower than at 25% o'clock, and cell of fuel cell is meritorious, the idle situation of exerting oneself.Show that through simulating, verifying it is efficient and practical that the control method for coordinating that generates electricity by way of merging two or more grid systems is united in the scene storage.

Claims (1)

1. the control method for coordinating that generates electricity by way of merging two or more grid systems is united in a scene storage, it is characterized in that it may further comprise the steps:

1) sets up wind energy converting unit Mathematical Modeling

(1.1) wind energy converting unit aerodynamics Mathematical Modeling is:

Wherein: P _MThe wind-powered electricity generation unit mechanical output that the energy of from wind, catching for the wind energy converting unit changes into, ρ _AirBe atmospheric density,

R is the wind turbine impeller radius,

λ is a tip speed ratio,

β is a propeller pitch angle,

C _pFor the wind energy conversion efficiency coefficient of blade, be the function of λ and β,

V _wBe wind speed;

(1.2) two matter piece axle coefficients of wind turbine and generator model equation is:

Wherein: H _T, H _GBe respectively wind turbine and generator inertia constant,

K _SBe the stiffness coefficient of axle,

D _T, D _GBe respectively the damping coefficient of wind turbine rotor and generator amature,

θ _SBe relative angular displacement between the two matter pieces,

T _M, T _EBe respectively wind turbine machine torque and generator electromagnetic torque,

ω _T, ω _GBe respectively wind turbine and generator amature rotating speed,

ω ₀Be synchronous speed;

(1.3) voltage equation of doubly fed induction generator does under the synchronous rotating frame

The magnetic linkage equation does

Wherein: ω _sBe coordinate system rotation angular speed,

U, i, ψ are respectively voltage, electric current and the magnetic linkage of winding,

R is the resistance of winding,

L _s, L _rBe respectively the self-induction of stator winding and rotor winding,

L _mBe the mutual inductance between the stator and rotor winding,

Subscript s, r represent the stator and the rotor of motor respectively,

Subscript d, q represent d, the q axle winding of motor respectively,

S is the revolutional slip of motor;

2) set up solar energy power generating unit Mathematical Modeling

Under the reference conditions, consider that the photovoltaic array V-I equation of intensity of solar radiation and variation of ambient temperature does

In the formula:

Wherein: V, I are respectively photovoltaic array terminal voltage and end electric current,

V _Oc, I _ScBe respectively open circuit voltage and short circuit current that photovoltaic module producer provides,

V _m, I _mBe respectively maximum power point voltage and electric current,

α, β are respectively under the reference radiation intensity, current temperature variation coefficient and voltage temperature variation coefficient,

R _sBe the series resistance of assembly,

is respectively the parallelly connected number and the serial number of assembly in the photovoltaic array;

3) set up the energy-storage units Mathematical Modeling

The batteries Mathematical Modeling does

In the formula:

Wherein: SOC is the state-of-charge of batteries,

A, B are respectively change in voltage coefficient and volume change coefficient,

K is the polarizing voltage constant,

C _pBe polarization reaction temperature compensation factor,

C _RBe the resistance temperature compensation factor,

T _bBe battery temp,

C _tBe the polarization reaction constant,

Q _nBe battery rating,

V _nBe the battery rated voltage,

E ₀Be the storage battery initial potential,

V is an accumulator voltage,

η is a battery efficiency,

N _B-sWith N _B-pBe respectively storage battery series connection and parallelly connected number;

The fuel cell Mathematical Modeling does

Wherein: V _FCBe the fuel cell terminal voltage,

I _FCBe fuel cell end electric current,

N _CellBe fuel cell series quantity,

A _CellBe the fuel cell electrode area,

R is a sheet resistance,

i _nBe the internal operation current density,

i _oBe inner nominal current density,

M and n are the junction loss constants;

4) unite the coordination control of generating electricity by way of merging two or more grid systems

(4.1) active power is coordinated control

Meritorious expression formula between wind light generation unit and the network load demand does

P _net＝P _wind+P _PV-P _load-P _sc (8)

Wherein: P _NetBe the active power amount of unbalance,

P _WindExert oneself for the wind energy converting unit is meritorious,

P _PVExert oneself for solar power generation unit is meritorious,

P _LoadBe the burden with power demand of dispatching of power netwoks appointment,

P _ScFor scene storage association system from electricity consumption, less because of accounting for the generating ratio, can ignore;

The superfluous P of wind light generation unit active power _Net＞0 o'clock, the power balance equation formula did

P _wind+P _PV＝P _load+P _ES-in (9)

Wherein: P _ES-inFor energy-storage units absorbs active power, promptly storage battery absorbs meritorious power P _Battery-in

Wind light generation cell power vacancy P _Net＜0 o'clock, the power balance equation formula did

P _wind+P _PV+P _ES-out＝P _load (10)

In the formula: P _ES-out=P _{Battery_out}+ P _FC

Wherein: P _ES-outThe active power of sending for energy-storage units,

P _{Battery_out}The active power of sending for storage battery,

P _FCThe active power of sending for fuel cell, when storage battery state-of-charge η＜25%, the starting fluid battery unit;

(4.2) reactive power is coordinated control

Wind light generation cell power factor range of operation is 0.85 to lag behind 0.85 in advance, is 0.98 to lag behind when normally moving, and the meritorious expression formula between wind light generation unit and the network load demand does

Q _net＝Q _Wind-out+Q _PV-out-Q _load (11)

Wherein: Q _NetBe the reactive power amount of unbalance,

Q _{Wind_out}Exert oneself for the wind energy converting unit is idle,

Q _{PV_out}Exert oneself for solar power generation unit is idle,

Q _LoadLoad or burden without work demand for the dispatching of power netwoks appointment;

The superfluous Q of wind light generation unit reactive power _Net＞0 and the point voltage that is incorporated into the power networks be higher than its rated voltage V＞V _RefThe time, the power balance equation formula does

Q _Wind-out+Q _PV-out＝Q _load+Q _ES-in+Q _STATCOM-in (12)

Wherein: Q _ES-inBe the reactive power of energy-storage units absorption,

Q _STATCOM-inReactive power for STATCOM STATCOM absorption;

Wind light generation cell power vacancy Q _Net＜0 o'clock, the power balance equation formula did

Q _load＝Q _Wind-out+Q _PV-out+Q _ES-out+Q _STATCOM-out V＜V _ref (13)

Or

Q _load+Q _STATCOM-in＝Q _Wind-out+Q _PV-out+Q _ES-out V＞V _ref (14)

Wherein: Q _ES-outThe reactive power of sending for energy-storage units,

Q _STATCOM-outThe reactive power of sending for STATCOM STATCOM;

When honourable unit runs on power factor 0.98 when leading, i.e. Q _Net＜0, the reactive power equilibrium equation

Q _Wind-in+Q _PV-in+Q _load＝Q _ES-out+Q _STATCOM-out V＜V _ref (15)

Or

Q _Wind-in+Q _PV-in+Q _load+Q _STATCOM-in＝Q _ES-out V＞V _ref (16)

Images