WIND POWER CONSUMPTION METHOD OF VIRTUAL POWER PLANT WITH CONSIDERATION OF COMPREHENSIVE DEMAND RESPONSES OF ELECTRICAL LOADS AND HEAT LOADS

Abstract

The present invention discloses a wind power consumption method of a virtual power plant with consideration of comprehensive demand responses of electrical loads and heat loads, which comprises: establishing a wind turbine output model, so as to obtain a wind power prediction curve; establishing heat load demand models before/after demand responses and heat supply equipment output models before/after the demand responses, so as to obtain the abandoned wind quantities per moment before/after the demand responses and the total abandoned wind quantities before/after the demand responses; judging that whether consumption is promoted or not according to the total abandoned wind quantities before/after the demand responses; and establishing a storage battery capacity model and judging the charging/discharging state and the charging/discharging capacity of a storage battery.

Claims

1. A wind power consumption method of a virtual power plant with consideration of comprehensive demand responses of electrical loads and heat loads, comprising the following steps: (1) establishing a wind turbine output model, so as to obtain a wind power prediction curve by the wind turbine output model; (2) establishing a heat load demand model before demand responses and a heat supply equipment output model before the demand responses, so as to obtain the demand of the heat loads before the demand responses according to the heat load demand model before the demand responses; taking the heat supply equipment output model before the demand responses as an electrical boiler output model before the demand responses in the virtual power plant and calculating the wind power quantity consumed by heat supply equipment before the demand responses according to the electrical boiler output model before the demand responses and the demand of the heat loads before the demand responses; and calculating the abandoned wind quantity per moment before the demand responses and the total abandoned wind quantity before the demand responses according to the wind power prediction curve, the demand of the electrical loads before the demand responses and the wind power quantity consumed by the heat supply equipment before the demand responses, wherein the virtual power plant comprises the electrical boiler, a storage battery and a wind turbine; (3) establishing a heat load demand model after the demand responses and a heat supply equipment output model after the demand responses, so as to obtain the demand of the heat loads after the demand responses according to the heat load demand model after the demand responses; and taking the heat supply equipment output model after the demand responses as an electrical boiler output model after the demand responses in the virtual power plant and calculating the wind power quantity consumed by the heat supply equipment after the demand responses according to the electrical boiler output model after the demand responses and the demand of the heat loads after the demand responses; (4) calculating the demand of the electrical loads after the demand responses; calculating the abandoned wind quantity per moment after the demand responses and the total abandoned wind quantity after the demand responses according to the wind power prediction curve, the demand of the electrical loads after the demand responses and the wind power quantity consumed by the heat supply equipment after the demand responses; then calculating a difference value between the total abandoned wind quantity before the demand responses and the total abandoned wind quantity after the demand responses; if the difference value is more than 0, judging that the wind power consumption is promoted; and if the difference value is less than 0, judging that the wind power consumption is not promoted; and (5) establishing a storage battery capacity model; and judging the charging/discharging state and the charging/discharging capacity of the storage battery according to the storage battery capacity model and the abandoned wind quantity per moment after the demand responses.

2. The wind power consumption method of the virtual power plant with consideration of comprehensive demand responses of electrical loads and heat loads according to claim 1, wherein in the step (1), the wind turbine output model is: $g_{\mathrm{WPP}}\left(t\right)=\{\begin{array}{cc}0& ,v_t\le v_{\mathrm{in}}\\ \frac{v_t-v_{\mathrm{in}}}{v_R-v_{\mathrm{in}}}g& ,v_{\mathrm{in}}\le v_i\ge v_R\\ g& ,v_R\le v_t\le v_{\mathrm{out}}\\ 0& ,v_t\ge v_{\mathrm{out}}\end{array}$ wherein in the formula, g represents the rated power of the wind turbine; v.sub.in represents the cut-in wind speed of the wind turbine; v.sub.R represents the rated wind speed of the wind turbine; v.sub.out represents the cut-out wind speed of the wind turbine; v.sub.t represents the real-time wind speed of the wind turbine at the moment t; t represents the moment t; and g.sub.WPP (t) represents the actual output of the wind turbine at the moment t; and the actual output of the wind turbine meets the following constraint condition:
g.sub.WPP.sup.min≤g.sub.WPP(t)≤g.sub.WPP.sup.max wherein in the formula, g.sub.WPP.sup.min represents the lower limit of power of the wind turbine; and g.sub.WPP.sup.max represents the upper limit of power of the wind turbine.

3. The wind power consumption method of the virtual power plant with consideration of comprehensive demand responses of electrical loads and heat loads according to claim 1, wherein in the step (2), the heat load demand model before the demand responses is: $Q_{\mathrm{heart}}^1\left(t\right)=\underset{t=1}{\overset{24}{.Math.}}\left(C_{\mathrm{air}}{\rho }_{\mathrm{air}}\mathrm{NSH}+\partial \mathrm{KA}\right)\left(T_{\mathrm{PMV}}^1-T_{\mathrm{out}}\left(t\right)\right)-Q_{\mathrm{inc}}-Q_{\mathrm{inh}}$ wherein in the formula, Q.sub.heart.sup.1(t) represents the demand of the heat loads before the demand responses; ∂ represents the modified temperature difference factor of an envelope enclosure; K represents the heat transfer coefficient of the envelope enclosure; A represents the area of the envelope enclosure; T.sub.out (t) represents the outdoor temperature at the moment t; C.sub.air represents the specific heat capacity of air; ρ.sub.air represents the air density; N represents the air change rate; S represents the area of a building; H represents the indoor height of the building; Q.sub.ine represents the calorific value of electrical equipment; Q.sub.inh represents the calorific value of a human body; and T.sub.PMV.sup.1 represents the set indoor temperature before the demand responses; the set indoor temperature T.sub.PMV.sup.1 before the demand responses is calculated by the following formula: $T_{\mathrm{PMV}}^1=T_S-\frac{M\left(2.43-{\lambda }_{\mathrm{PMV}}^1\right)\left(I_{\mathrm{cl}}+0.1\right)}{3.76};$ in the formula, T.sub.S represents the human skin temperature in a normal temperature state; M represents the human energy metabolism rate; λ.sub.PMV.sup.1 represents the initial PMV (Predicted Mean Vote) index; and I.sub.c1 represents the dress heat resistance; the electrical boiler output model before the demand responses is:
Q.sup.1.sub.EB(t)=g.sup.1.sub.EB(t).Math.η.sub.EB wherein in the formula, Q.sup.1.sub.EB(t) represents the power of heat supply of an electrical boiler before the demand responses at the moment t; g.sup.1.sub.EB(t) represents the wind power quantity consumed by the work of the electrical boiler before the demand responses at the moment t; and η.sub.EB represents the electricity to heat conversion efficiency; the actual output of the electrical boiler before the demand responses meets the following constraint condition:
Q.sup.1min.sub.EB≤Q.sup.1.sub.EB(t)≤Q.sup.1max.sub.EB wherein in the formula, Q.sup.1min.sub.EB represents the lower limit of power of the electrical boiler before the demand responses; Q.sup.1max.sub.EB represents the upper limit of power of the electrical boiler before the demand responses; and Q.sup.1.sub.EB(t) represents the actual output of the electrical boiler before the demand responses at the moment t.

4. The wind power consumption method of the virtual power plant with consideration of comprehensive demand responses of electrical loads and heat loads according to claim 3, wherein the heat supply equipment before the demand responses only comprises the electrical boiler before the demand responses, so the numerical value of Q.sup.1.sub.EB(t) is equal to that of Q.sub.heart.sup.1(t); and the output of the electrical boiler before the demand responses at the moment t is obtained according to the demand Q.sub.heart.sup.1(t) of the heat loads before the demand responses, so as to obtain the wind power quantity g.sup.1.sub.EB(t) consumed by the heat supply equipment before the demand responses according to the electrical boiler output model before the demand responses.

5. The wind power consumption method of the virtual power plant with consideration of comprehensive demand responses of electrical loads and heat loads according to claim 1, wherein in the step (2), the abandoned wind quantity per moment before the demand responses and the total abandoned wind quantity before the demand responses are respectively obtained by the following formulas: $g_W^1\left(t\right)=g_{\mathrm{WPP}}\left(t\right)-g_E^1\left(t\right)-g_{\mathrm{EB}}^1\left(t\right)g_W^1=\underset{t=1}{\overset{24}{.Math.}}{g}_W^1\left(t\right)$ wherein in the formulas, g.sub.W.sup.1 represents the total abandoned wind quantity before the demand responses; g.sub.W.sup.1(t) represents the abandoned wind quantity before the demand responses at the moment t; g.sub.E.sup.1(t) represents the demand of the electrical loads before the demand responses at the moment t; and g.sub.EB.sup.1(t) represents the wind power quantity consumed by the heat supply equipment before the demand responses.

6. The wind power consumption method of the virtual power plant with consideration of comprehensive demand responses of electrical loads and heat loads according to claim 1, wherein in the step (3), the heat load demand model after the demand responses is: $Q_{\mathrm{heart}}^2\left(t\right)=\underset{t=1}{\overset{24}{.Math.}}\left(C_{\mathrm{air}}{\rho }_{\mathrm{air}}\mathrm{NSH}+\partial \mathrm{KA}\right)\left(T_{\mathrm{PMV}}^2-T_{\mathrm{out}}\left(t\right)\right)-Q_{\mathrm{inc}}-Q_{\mathrm{inh}}$ wherein in the formula, Q.sub.heart.sup.2(t) represents the demand of the heat loads after the demand responses; T.sub.PMV.sup.2 represents the set indoor temperature after the demand responses; ∂ represents the modified temperature difference factor of the envelope enclosure; K represents the heat transfer coefficient of the envelope enclosure; A represents the area of the envelope enclosure; T.sub.out (t) represents the outdoor temperature at the moment t; C.sub.air represents the specific heat capacity of air; ρ.sub.air represents the air density; N represents the air change rate; S represents the area of a building; H represents the indoor height of the building; Q.sub.ine represents the calorific value of the electrical equipment; and Q.sub.inh represents the calorific value of the human body; the set indoor temperature T.sub.PMV.sup.2 after the demand responses is calculated by the following formula: $T_{\mathrm{PMV}}^2=T_S-\frac{M\left(2.43-{\lambda }_{\mathrm{PMV}}^2\right)\left(I_{\mathrm{cl}}+0.1\right)}{3.76},$ wherein in the formula, T.sub.S represents the human skin temperature in a normal temperature state; M represents the human energy metabolism rate; λ.sub.PMV.sup.2 represents the initial PMV index; I.sub.c1 represents the dress heat resistance; the electrical boiler output model after the demand responses is:
Q.sup.2.sub.EB≤Q.sup.2.sub.EB(t)≤η.sub.EB wherein in the formula, Q.sup.2.sub.EB(t) represents the power of heat supply of the electrical boiler after the demand responses at the moment t; g.sup.2.sub.EB(t) represents the wind power quantity consumed by the work of the electrical boiler after the demand responses at the moment t; and η.sub.EB a represents the electricity to heat conversion efficiency; the actual output of the electrical boiler after the demand responses meets the following constraint condition:
Q.sup.2min.sub.EB≤Q.sup.2.sub.EB(t)≤Q.sup.2max.sub.EB wherein in the formula, Q.sup.2min.sub.EB represents the lower limit of power of the electrical boiler after the demand responses; Q.sup.2max.sub.EB represents the upper limit of power of the electrical boiler after the demand responses; and Q.sup.2.sub.EB(t) represents the actual output of the electrical boiler after the demand responses at the moment t.

7. The wind power consumption method of the virtual power plant with consideration of comprehensive demand responses of electrical loads and heat loads according to claim 6, wherein the heat supply equipment after the demand responses only comprises the electrical boiler after the demand responses, so the numerical value of Q.sup.2.sub.EB(t) is equal to that of Q.sub.heart.sup.2(t); and the output of the electrical boiler after the demand responses at the moment t is obtained according to the demand Q.sub.heart.sup.2(t) of the heat loads after the demand responses, so as to obtain the wind power quantity g.sup.2.sub.EB(t) consumed by the heat supply equipment after the demand responses according to the electrical boiler output model after the demand responses.

8. The wind power consumption method of the virtual power plant with consideration of comprehensive demand responses of electrical loads and heat loads according to claim 1, wherein in the step (4), the demand of the electrical loads after the demand responses is obtained specifically by adopting the following manners: firstly, the variation of the electrical loads per moment after the demand responses is calculated by the following formulas: $\Delta g_{\mathrm{on}}\left(t\right)=e_{\mathrm{on}}\times g_{\mathrm{on}}\left(t\right)\times \frac{\Delta P_{\mathrm{on}}\left(t\right)}{P_{\mathrm{on}}\left(t\right)},t\in T_{\mathrm{on}}\Delta g_{\min d}\left(t\right)=e_{\min d}\times g_{\min d}\left(t\right)\times \frac{\Delta P_{\min d}\left(t\right)}{P_{\min d}\left(t\right)},t\in T_{\mathrm{mind}}\Delta g_{\mathrm{off}}\left(t\right)=e_{\mathrm{off}}\times g_{\mathrm{off}}\left(t\right)\times \frac{\Delta P_{\mathrm{off}}\left(t\right)}{P_{\mathrm{off}}\left(t\right)},t\in T_{\mathrm{off}}$ wherein in the formulas, Δg.sub.on(t), Δg.sub.mind(t) and Δg.sub.off(t) represent the variations of the electrical loads at the peak period of power consumption, the flat period of power consumption and the valley period of power consumption after the demand responses; g.sub.on(t), g.sub.mind(t) and g.sub.off(t) represent the electrical loads at the peak period of power consumption, the flat period of power consumption and the valley period of power consumption before the demand responses; P.sub.on(t), P.sub.mind(t) and P.sub.off(t) represent the energy consumption at the peak period of power consumption, the flat period of power consumption and the valley period of power consumption before the demand responses; ΔP.sub.on(t), ΔP.sub.mind(t) and ΔP.sub.off(t) represent the variations of energy consumption at the peak period of power consumption, the flat period of power consumption and the valley period of power consumption after the demand responses; e.sub.on, e.sub.mind and e.sub.off represent the elastic coefficients of energy consumption at the peak period of power consumption, the flat period of power consumption and the valley period of power consumption; and T.sub.on represents the peak period of energy consumption, T.sub.mind represents the flat period of energy consumption, and T.sub.off represents the valley period of energy consumption; and then, the demand g.sub.E.sup.2(t) of the electrical loads per moment after the demand responses is obtained by the variations of the electrical loads per moment after the demand responses plus the demand g.sub.E.sup.1(t) of the electrical loads per moment before the demand responses.

9. The wind power consumption method of the virtual power plant with consideration of comprehensive demand responses of electrical loads and heat loads according to claim 1, wherein in the step (4), the abandoned wind quantity per moment after the demand responses and the total abandoned wind quantity after the demand responses are respectively obtained by the following formulas: $g_W^2\left(t\right)=g_{\mathrm{WPP}}\left(t\right)-g_E^2\left(t\right)-g_{\mathrm{EB}}^2\left(t\right)g_W^2=\underset{t=1}{\overset{24}{.Math.}}{g}_W^2\left(t\right)$ wherein in the formulas, g.sub.W.sup.1 represents the total abandoned wind quantity after the demand responses; g.sub.W.sup.2(t) represents the abandoned wind quantity after the demand responses at the moment t; g.sub.E.sup.2(t) represents the demand of the electrical loads after the demand responses at the moment t; and g.sub.EB.sup.2(t) represents the quantity of electricity consumed by heat supply of the electrical boiler after the demand responses at the moment t.

10. The wind power consumption method of the virtual power plant with consideration of comprehensive demand responses of electrical loads and heat loads according to claim 1, wherein the step (5) specifically is: the storage battery capacity model is:
S.sub.soc(t)=S.sub.soc(t−1)+(S.sub.ch(t)−S.sub.dis(t))
S.sub.ch(t)=g.sub.ch(t)η.sub.ch
S.sub.dis(t)=g.sub.dis(t)η.sub.dis wherein in the formulas, S.sub.SOC.sup.min (t) represents the capacitance of the storage battery at the moment t; S.sub.soc (t−1) represents the capacitance of the storage battery at the moment (t−1); S.sub.ch (t) represents the charging capacity of the storage battery at the moment t; S.sub.dis (t) represents the discharging capacity of the storage battery at the moment t; g.sub.ch (t) represents the charging power of the storage battery at the moment t; η.sub.ch represents the charging efficiency of the storage battery; g.sub.dis (t) represents the discharging power of the storage battery at the moment t; and η.sub.dis represents the discharging efficiency of the storage battery; the capacity of the storage battery meets the following constraint condition:
S.sub.SOC.sup.min≤S.sub.SOC(t)≤S.sub.SOC.sup.max wherein in the formula, S.sub.SOC.sup.min represents the minimum charging capacity of the storage battery, and the S.sub.SOC.sup.max represents the maximum charging capacity of the storage battery; the output constraint of the storage battery meets the following constraint condition:
g.sub.ch.sup.min≤g.sub.ch(t)≤g.sub.ch.sup.max,
g.sub.dis.sup.min≤g.sub.dis(t)≤g.sub.dis.sup.max wherein in the formulas, g.sub.ch.sup.min represents the minimum charging power of the storage battery; g.sub.ch.sup.max represents the maximum charging power of the storage battery; g.sub.dis.sup.min represents the minimum discharging power of the storage battery; and g.sub.dis.sup.max represents the maximum discharging power of the storage battery; when the value of the actual output g.sub.WPP (t) of the wind turbine at the moment t is less than the sum of the value of the wind power quantity g.sup.2.sub.EB(t) consumed by the heat supply equipment after the demand responses at the moment t and the value of the demand g.sub.E.sup.2(t) of the electrical loads after the demand responses at the moment t, the value of the obtained abandoned wind quantity g.sub.W.sup.2(t) after the demand responses at the moment t is less than 0; and when the value of the actual output g.sub.WPP (t) of the wind turbine at the moment t is more than or equal to the sum of the value of the wind power quantity g.sup.2.sub.EB(t) consumed by the heat supply equipment after the demand responses at the moment t and the value of the demand g.sub.E.sup.2(t) of the electrical loads after the demand responses at the moment t, the value of the obtained abandoned wind quantity g.sub.W.sup.2(t) after the demand responses at the moment t is more than or equal to 0; and when the value of the abandoned wind quantity g.sub.W.sup.2(t) after the demand responses at the moment t is less than 0, the storage battery is configured to discharge to assist the wind turbine to supply power; when the value of the abandoned wind quantity g.sub.W.sup.2(t) after the demand responses at the moment t is more than or equal to 0, the storage battery is configured to be charged; the charging quantity S.sub.ch (t) is equal to the abandoned wind quantity g.sub.W.sup.2(t) after the demand responses at the moment t; and the discharging quantity S.sub.dis (t) is equal to the absolute value of the abandoned wind quantity g.sub.W.sup.2(t) after the demand responses at the moment t.

Description

DESCRIPTION OF THE DRAWINGS

[0055] FIG. 1 is a schematic diagram of wind power consumption of a virtual power plant of the present invention;

[0056] FIG. 2 is a diagram of a demand curve of electrical loads and heat loads in an embodiment of the present invention;

[0057] FIG. 3 is a diagram of the condition of wind power consumption in a scenario 1 of the embodiment of the present invention;

[0058] FIG. 4 is a diagram of a demand change curve after demand responses of the electrical loads in the embodiment of the present invention;

[0059] FIG. 5 is a diagram of a demand change curve after demand responses of the heat loads in the embodiment of the present invention; and

[0060] FIG. 6 is a diagram of the condition of wind power consumption in a scenario 2 of the embodiment of the present invention.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

[0061] The present invention is further described in details hereinafter through combination with the drawings and specific embodiments.

[0062] In the implementation of the present invention, the present invention is implemented according to the specific steps in the contents of the description and the contents of the claims, and the specific step process is not described here.

[0063] The embodiments of the present invention are described as follows:

[0064] In the embodiments, a virtual power plant which comprises a wind turbine, a storage battery and an electrical boiler is taken as an example. The output condition of the wind turbine is shown in FIG. 2. The initial capacity of the storage battery is 0 kWh. A demand curve of electrical loads and heat loads is shown in FIG. 3. The uniform energy consumption cost before demand responses is 0.5/kWh; after peak-valley energy consumption is adopted, the peak period of energy consumption is 11:00-15:00 and 19:00-21:00; the cost at the peak period of energy consumption is 0.8/kWh; and the corresponding elastic coefficient e.sub.on of energy consumption is −0.07. The flat period of energy consumption is 08:00-10:00, 16:00-18:00 and 22:00-23:00; the cost at the flat period of energy consumption is 0.4/kWh; and the corresponding elastic coefficient e.sub.mind of energy consumption is 0.05. The valley period of energy consumption is 24:00-07:00; the cost at the valley period of energy consumption is 0.2/kWh; and the corresponding elastic coefficient e.sub.off of energy consumption is 0.08. In the embodiments, two scenarios are considered. In the first scenario, energy storage equipment is not considered, and the demand responses are not implemented on the electrical loads and the heat loads. In the second scenario, the energy storage equipment is considered, and comprehensive demand responses are implemented on the electrical loads and the heat loads. In the embodiment, the abandoned wind quantities of the virtual power plant in the two scenarios are respectively calculated and are compared with each other.

[0065] The virtual power plant in the embodiment is mainly used for wind power consumption for the power utilization of users and the electricity to heat conversion of the electrical boiler. In the scenario 2, the demand of the electrical loads after the demand responses is obtained according to the elastic cost coefficients of the electrical loads at different periods. It can be seen from FIG. 4 that peak load shifting of the electrical loads is realized by the demand responses. At the demand peak period of the heat loads, incentive demand responses are implemented by lowering the set indoor temperature T.sub.PMV, and the compensation cost for the demand responses of the virtual power plant is considered to be minimum. In the embodiment, the incentive demand responses of the heat loads are implemented at the common peak period (19:00-23:00) of the electrical loads and the heat loads. It can be seen from FIG. 5 that the curtailment of the heat loads is realized by the demand responses; meanwhile, when the abandoned wind quantity is larger, the storage battery is configured to store energy, so as to promote wind power consumption. Through comparison on the conditions of wind power consumption in the two scenarios, as shown in Tab. 1, it is found that the abandoned wind quantity in the scenario 2 is reduced, and the consumption capacity is higher.

TABLE-US-00001 TABLE 1 Table of conditions of wind power consumption of a virtual power plant in two scenarios Scenario 1 Scenario 2 Wind power generation/kW 4886 4886 Electrical load/kW 2177 2184.3 Output of electrical boiler/kW 2656.66 2653.89 Output of storage battery/kW 0 12.76 Abandoned wind quantity/kW 52.34 35.05

[0066] Therefore, the comprehensive responses of the electrical loads and the heat loads are considered in the present invention, the capacity of wind power consumption of the virtual power plant is effectively improved; the reference is provided for a research that the comprehensive demand responses are implemented on various loads; and an effective way is provided for promoting wind power consumption of a power generation and distribution system, such as the virtual power plant and the like and relieving the problem of wind abandonment.

[0067] Finally, it should be noted that the above example is only used for describing the technical solution of the present invention, rather than the limit to technical solution of the present invention. Although the present invention is described with reference to the above example, those skilled in the art should understand that a specific implementation manner of the present invention can still be modified or equivalently replaced, and any modification or equivalent replacement made without departing from the spirit and the scope of the present invention shall be covered in the scope of the claims of the present invention.