METHOD AND DEVICE FOR ANALYZING VIRTUAL POWER PLANT OPERATION RISK

Abstract

Disclosed are a virtual power plant operation risk analysis method and a device. The method comprises: establishing a multi-state model of wind turbine output, analyzing influence of wind speed on wind turbine failure rate based on the multi-state model of wind turbine output, and establishing a wind turbine failure model considering wind turbine time-varying failure rate; establishing a multi-state model of wind turbine output considering the wind speed and the wind turbine time-varying failure rate by an improved general generating function method based on the multi-state model of wind turbine output and the wind turbine failure model considering the wind turbine time-varying failure rate; establishing a multi-state output model of virtual power plant based on the multi-state model of wind turbine output considering wind speed and wind turbine time-varying failure rate; and calculating operation risk indicators of virtual power plant through the multi-state output model of virtual power plant.

Claims

1. A method for analyzing virtual power plant operation risks, comprising: establishing a multi-state model of wind turbine output; analyzing influence of a wind speed on a wind turbine time-varying failure rate based on the multi-state model of the wind turbine output, and establishing a wind turbine failure model considering a wind turbine time-varying failure rate; establishing a multi-state model of the wind turbine output considering the wind speed and the wind turbine time-varying failure rates by a general generating function method based on the multi-state model of the wind turbine output and the wind turbine failure model considering the wind turbine time-varying failure rate; establishing a multi-state output model of a virtual power plant based on the multi-state model of the wind turbine output considering the wind speed and the wind turbine time-varying failure rate; and calculating operation risk indicators of the virtual power plant through the multi-state output model of the virtual power plant, and completing analysis of virtual power plant operation risks; wherein the step of establishing the multi-state model of the wind turbine output comprises: analyzing a relationship between the wind speed and the wind turbine output without considering wind turbine failure, building the multi-state model of the wind turbine output based on the relationship between the wind speed and the wind turbine output, dividing a wind speed s(t) into K.sub.s states, modelling the wind turbine output with Markov process, and dividing wp.sub.i.sup.1(t) into K.sub.s states, obtaining a time-varying probability value q.sub.i,k(t) of the wind turbine output wp.sub.i,k.sup.1 of a k.sub.s-th state, and building the multi-state model of wind turbine output by using the improved general generating function method, wherein wp.sub.i.sup.1(t) represents an output of a wind turbine i at the wind speed S(t) at time t.

2. (canceled)

3. The method of claim 1, wherein the relationship between the wind speeds and wind turbine output is: $w{p}_i^1\left(t\right)=\{\begin{array}{c}0,0s\left(t\right)s_i^{\mathrm{ci}}\mathrm{or}s\left(t\right)>s_i^{\mathrm{co}}\\ a_i{s\left(t\right)}^3+b_i,s_i^{\mathrm{ci}}s\left(t\right)s_i^c\\ {\mathrm{wp}}_i^r,s_i^{c}s\left(t\right)s_i^{\mathrm{co}}\end{array},$ wherein t represents time, i represents a serial number of a wind turbine, wp.sub.i.sup.1(t) represents an output of a wind turbine i at the wind speed S(t) at time t, s.sub.i.sup.ci, s.sub.i.sup.c, s.sub.i.sup.co represent a cut-in wind speed, a rated wind speed and a cut-out wind speed of the wind turbine i, respectively, wp.sub.i.sup.r represents a rated power of the wind turbine i; a.sub.i and b.sub.i are correlation coefficients between the output of the wind turbine and the wind speed, respectively, $a_i=\frac{{\mathrm{wp}}_i^r}{{\left(s_i^c\right)}^3-{\left(s_i^{ci}\right)}^3},$ $b_i=\frac{{\left({\mathrm{wp}}_i^r{s}_i^{\mathrm{ci}}\right)}^3}{{\left(s_i^{\mathrm{ci}}\right)}^3-{\left(s_i^c\right)}^3}.$ wherein the time-varying probability value q.sub.i,k(t) of the wind turbine output wp.sub.i,k.sup.1 of the k.sub.s-th state is: $\{\begin{array}{c}\frac{{\mathrm{dq}}_{i,k}\left(t\right)}{\mathrm{dt}}=\left[\underset{k=1,kl}{\overset{K_s}{.Math.}}{q}_{i,l}\left(t\right){}_{k,l}^s\right]-q_{i,k}\left(t\right)\underset{l=1,lk}{\overset{K_s}{.Math.}}{}_{k,l}^{s}k=1,.Math.,K_s\\ q_{i,k}\left(t_0\right)=1,q_{i,l}\left(t_0\right)=0,kl\end{array}.$ wherein .sub.k,l.sup.s represents a state transition rate of the wind turbine output from the k.sub.s-th state to a l.sub.s-th state, q.sub.i,k(t.sub.0) is a time-varying probability value of the wind turbine output in the k.sub.s-th state of the wind turbine i at time t.sub.0, q.sub.i,l(t.sub.0) represents a time-varying probability value of the wind turbine output in the l.sub.s-th state of the wind turbine i at time t.sub.0, and q.sub.i,l(t) represents a time-varying probability value of the wind turbine output in a l state of the wind turbine i at time t; wherein the multi-state model of the wind turbine output is: $u_l^1\left(z,t\right)=\underset{k_s=1}{\overset{K_s}{.Math.}}{q}_{i,k}\left(t\right).Math.z^{{\mathrm{wp}}_{i,k}^1},$ wherein u.sub.i.sup.1(z,t) represents an improved general generating function representation method of output of the wind turbine i regardless of the wind turbine failure, Z represents a state value of a random variable, and z.sup.wp.sup.i,k.sup.1 represents an output value of the wind turbine i is wp.sub.i,k.sup.1.

4. The method of claim 1, wherein the wind turbine failure model considering the wind turbine time-varying failure rate comprises: analyzing the influence of the wind speed on the wind turbine failure rate, and establishing the wind turbine failure model;
.sub.I(t)=.sub.i,0+.sub.i,s(t), wherein .sub.I (t) represents the time-varying failure rates of the wind turbine i at time t, .sub.i,0 represents a basic failure rate of the wind turbine i, and .sub.i,s(t) represents a variable failure rate of the wind turbine i caused by the wind speed at time t; wherein a relationship model between the variable failure rate of the wind turbine i caused by the wind speed at time t and the wind speed s(t) is as follows:
.sub.i,s(t)=(.sub.i,maxs(t).sup.2.sub.i,mins(t).sup.2)/(s.sub.i.sup.co2s.sub.i.sup.ci2)+c.sub.s, wherein .sub.i,max represents the wind turbine failure rate corresponding to the cut-out wind speed s.sub.i.sup.co of the wind turbine i, .sub.i,min represents the wind turbine failure rate corresponding to the cut-in wind speed S.sub.i.sup.ci of the wind turbine i, and C.sub.s represents a constant related to the cut-in wind speed and the cut-out wind speed.

5. The method of claim 4, wherein the basic failure rate of the wind turbine at different wind speed is described by the multi-state model, and the variable failure rate of the wind turbine at different wind speed is considered, and a failure probability of the wind turbine is obtained as follows:
q.sub.i,k.sup.wt(t)=1e.sup.-.sup.i,k.sup.t, wherein t is time, q.sub.i,k.sup.wt(t) represents failure probability of the wind turbine iat a k.sub.s-th wind speed at time t , and .sub.i,k failure rate of wind turbine i at the k.sub.s-th wind speed at time t; based on the failure probability of the wind turbine, the failure model of the wind turbine i in the k.sub.s-th wind speed state is established by using the improved general generating function method:
u.sub.i.sup.2(z,t)=(1q.sub.i,k.sup.wt(t).Math.z.sup.1+q.sub.i,k.sup.wt(t).Math.z.sup.0=e.sup.-.sup.i,k.sup.t.Math.z.sup.1+(1e.sup.-.sup.i,k.sup.t).Math.z.sup.0, wherein u.sub.i.sup.2(z,t) represents the improved general generating function representation method of the failure model of the wind turbine i considering a influence of wind speed on wind turbine failure probability, z.sup.1 indicates that wind turbine i is in normal operation and z.sup.0 represents the wind turbine i is in failure.

6. The method of claim 1, wherein the multi-state model of wind turbine output considering the wind speed and the wind turbine time-varying failure rate is: $\begin{array}{c}{u}_i^w\left(z,t\right)={}_{\mathrm{ser}}\left\{u_i^1\left(z,t\right),u_i^2\left(z,t\right)\right\}\\ ={}_{\mathrm{ser}}\left\{\underset{k_s=1}{\overset{K_s}{.Math.}}{q}_{i,k}\left(t\right).Math.z^{{\mathrm{wp}}_{i,k}^1},\left(1-q_{i,k}^{\mathrm{wt}}\left(t\right)\right).Math.z^1+q_{i,k}^{\mathrm{wt}}\left(t\right).Math.z^0\right\}\\ ={}_{\mathrm{ser}}\left\{\underset{k_s=1}{\overset{K_s}{.Math.}}{q}_{i,k}\left(t\right).Math.z^{{\mathrm{wp}}_{i,k}^1},e^{-{}_{i,k}t}.Math.z^1+\left(1-e^{-{}_{i,k}t}\right).Math.z^0\right\}\\ =\underset{k_s=1}{\overset{K_s}{.Math.}}{q}_{i,k}\left(t\right).Math.\left(1-q_{i,k}^{\mathrm{wt}}\left(t\right)\right).Math.z^{{\mathrm{wp}}_{i,k}^1}+\underset{k_s=1}{\overset{K_s}{.Math.}}{q}_{i,k}\left(t\right).Math.q_{i,k}^{\mathrm{wt}}\left(t\right).Math.z^0\\ =\underset{j=1}{\overset{n_w}{.Math.}}{q}_{i,j}^w\left(t\right).Math.z^{{\mathrm{wp}}_{i,j}^1}\end{array},$ wherein u.sub.i.sup.w(z,t) represents the general generating function representation method of the output model of the wind turbine i considering the wind speed and the wind turbine time-varying failure rates, .sub.ser represents a series operator, q.sub.i,j.sup.w(t) represents a probability of the wind turbine i in state j, and Z.sup.wp.sup.i,j.sup.w represents an output value of the wind turbine i in state j.

7. The method of claim 1, wherein establishing the multi-state model of the virtual power plant composed of a plurality of distributed wind powers by the multi-state model of wind turbine output considering the wind speeds and the wind turbine time-varying failure rate: $\begin{array}{c}{U}^{\mathrm{VPP}}\left(z,t\right)={}_{\mathrm{par}}\left\{u_1^w\left(z,t\right),.Math.,u_i^w\left(z,t\right),.Math.,u_{N_w}^w\left(z,t\right)\right\}\\ =\underset{i=1}{\overset{N_w}{.Math.}}\underset{{j=1}^{}}{\overset{n_w}{.Math.}}{q}_{i,j}^w\left(t\right).Math.z^{{\mathrm{wp}}_j^w}\\ =\underset{m=1}{\overset{N_{\mathrm{vpp}}}{.Math.}}{q}_m^{\mathrm{vpp}}\left(t\right).Math.z^{{\mathrm{VPP}}_m}\end{array},$ wherein u.sup.VPP(z,t) represents the general generating function representation method for the output model of the virtual power plant by aggregating N.sub.w independent wind turbines, and .sub.par represents a parallel operator, q.sub.m.sup.vpp(t) represents a probability of the virtual power plant in state m, and z.sup.VPP.sup.m in represents that an output value of the virtual power plant in state m is VPP.sub.m.

8. The method of claim 1, wherein calculating the operation risk indicators of the virtual power plant: $D\left(t\right)=\underset{m}{.Math.}{q}_m^{\mathrm{vpp}}\left(t\right),{\mathrm{VPP}}_{m}L,$ $E\left(t\right)=\underset{w}{.Math.}{q}_m^{\mathrm{vpp}}\left(t\right).Math.{\mathrm{VPP}}_m,$ $A\left(T\right)=\underset{t=0}{\overset{T}{.Math.}}\left(L-{\mathrm{VPP}}_m\right)q_m^{\mathrm{vpp}}\left(t\right)\mathrm{cdf}\left(\right),{\mathrm{VPP}}_m<L,$ wherein D(t) is a power supply shortage probability, E(t) is an expected power supply shortage and A(t) is a power supply shortage loss for industrial users, L represents load values of industrial users powered by the virtual power plant, T represents total power supply duration of the virtual power plant, t represents the time and t [0,T], represents duration of power outage, and cdf () represents loss function of power supply shortage of the industrial users, and is related to the duration of power outage.

9. A device for analyzing virtual power plant operation risks, used to realize the method of claim 1, comprising: a wind speed and wind turbine output module, used to construct a relationship model between the wind turbine output and the wind speed, divide the wind speed into multiple states, establish a multi-state wind speed model; calculate wind turbine output values and corresponding probability values without considering wind turbine failure according to the relationship model between the wind turbine output and the wind speed and the multi-state output model of the wind turbine; a wind turbine time-varying failure rate acquisition module, used for acquiring variable failure rate of the wind turbine, and obtaining the wind turbine time-varying failure rate by adding the basic failure rate of the wind turbine and variable failure rate of the wind turbine caused by the wind speed; a wind turbine output module considering wind speed and wind turbine time-varying failure rate, used for constructing a wind turbine failure model considering wind turbine time-varying failure rate based on the wind turbine time-varying failure rate acquisition module; obtaining the wind turbine output value and a corresponding probability value considering the wind speed and the wind turbine time-varying failure rate based on the wind turbine output obtained by the wind speed and wind turbine output module, and a virtual power plant operation risk assessment module, used for constructing a virtual power plant output model comprising a plurality of distributed wind powers; establishing a virtual power plant operation risk indicator system comprising a power supply shortage probability, an expected power supply shortage and a power supply shortage loss of industrial users in the virtual power plant, and calculating the power supply shortage probability, the expected power supply shortage and the power supply shortage loss of industrial users in the virtual power plant.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0043] In order to more clearly explain the embodiments of the present application or the technical solutions in the prior art, the following will briefly introduce the drawings that need to be used in the embodiments. Obviously, the drawings in the following description are only some embodiments of the present application. For those of ordinary skill in the art, other drawings may be obtained according to these drawings without any creative labor.

[0044] FIG. 1 is a flow chart of a method for analyzing virtual power plant operation considering time-varying failure of distributed wind power provided by an embodiment of the present application.

[0045] FIG. 2 is a schematic diagram of a virtual power plant structure model in the virtual power plant operation analysis method considering the time-varying failure rate of distributed wind power provided by an embodiment of the present application.

[0046] FIG. 3 is a probability result diagram of power supply shortage of virtual power plant in an embodiment of the present application, whether time-varying failure rate of distributed wind power is considered or not.

[0047] FIG. 4 is a structural schematic diagram of a virtual power plant operation risk analysis device considering time-varying failure rate of distributed wind power provided by an embodiment of the present application.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0048] The technical solutions in embodiments of the present application will be clearly and completely described below with reference to the drawings in embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, but not all of them. Based on the embodiment of the present application, all other embodiments obtained by ordinary technicians in the field without creative labor are within the scope of the present application.

[0049] In order to make the above objects, features and advantages of the present application more obvious and understandable, the present application will be explained in further detail below with reference to the drawings and detailed description.

[0050] A flow chart of a virtual power plant operation risk analysis method shown in FIG. 1 specifically includes:

[0051] S1, A multi-state model of wind turbine output without considering wind turbine failure is established.

[0052] The S1 specifically includes: [0053] (1) relationship between wind speed and wind turbine output is analyzed without considering wind turbine failure.

[0054] When the wind turbine is running well, relationship between output of wind turbine i and the wind speed is expressed by the following formula:

[00008]$w{p}_i^1\left(t\right)=\{\begin{array}{c}0,0s\left(t\right)s_i^{\mathrm{ci}}\mathrm{or}s\left(t\right)>s_i^{\mathrm{co}}\\ a_i{s\left(t\right)}^3+b_i,s_i^{\mathrm{ci}}s\left(t\right)s_i^c\\ {\mathrm{wp}}_i^r,s_i^{c}s\left(t\right)s_i^{\mathrm{co}}\end{array},$ [0055] where t represents time, i represents a serial number of the wind turbines, wp.sub.i.sup.1(t) indicates output of the wind turbine i at the wind speed S(t) at time t, s.sub.i.sup.ci, s.sub.i.sup.c, s.sub.i.sup.co represent cut-in wind speed, rated wind speed and cut-out wind speed of the wind turbine i, respectively, wp.sub.i.sup.r represents rated power of the wind turbine i; a.sub.i and b.sub.i are correlation coefficients between the output of the wind turbine and the wind speed, respectively,

[00009]$a_i=\frac{{\mathrm{wp}}_i^r}{{\left(s_i^c\right)}^3-{\left(s_i^{ci}\right)}^3},$$b_i=\frac{{\left({\mathrm{wp}}_i^r{s}_i^{\mathrm{ci}}\right)}^3}{{\left(s_i^{\mathrm{ci}}\right)}^3-{\left(s_i^c\right)}^3}.$ [0056] (2) The multi-state model of wind turbine output is established without considering wind turbine failure.

[0057] The wind speed S(t) is divided into K.sub.s states, and a multi-state wind speed model is established. The wind speed of k.sub.s(k.sub.s=1, . . . ,K.sub.s)-th state is s.sub.k, and s.sub.k<s.sub.k+1. Transition rate between different wind speed states k.sub.s, l.sub.s is .sub.k,l.sup.s where k.sub.s and l.sub.s represent serial numbers of different wind speed states, respectively.

[0058] The multi-state model of wind turbine output is established according to the relationship between wind turbine output and wind speed in normal operation of wind turbine. Markov process is used to model the wind turbine output, wp.sub.i.sup.1(t) is divided into K.sub.s states, and the wind turbine output in the k.sub.s(k.sub.s=1, . . . , K.sub.s)-th state is wp.sub.i,k.sup.1, and wp.sub.i,k.sup.1<wp.sub.i,k+1.sup.1. Time-varying probability value q.sub.i,k(t) of the wind turbine output wp.sub.i,k.sup.1 in the k.sub.s-th state is obtained according to differential equations of multi-state Markov process as follows:

[00010]$\{\begin{array}{c}\frac{{\mathrm{dq}}_{i,k}\left(t\right)}{\mathrm{dt}}=\left[\underset{k=1,kl}{\overset{K_s}{.Math.}}{q}_{i,l}\left(t\right){}_{k,l}^s\right]-q_{i,k}\left(t\right)\underset{l=1,lk}{\overset{K_s}{.Math.}}k=1,.Math.,K_s\\ q_{i,k}\left(t_0\right)=1,q_{i,l}\left(t_0\right)=0,kl\end{array}.$

[0059] A multi-state model of the output of the wind turbine i is established by using an improved general generating function method and is expressed by following formula.

[00011]$u_l^1\left(z,t\right)=\underset{k_s=1}{\overset{K_s}{.Math.}}{q}_{i,k}\left(t\right).Math.z^{{\mathrm{wp}}_{i,k}^1},$ [0060] where u.sub.i.sup.1(z,t) represents improved general generating function representation method of the output of the wind turbine i regardless of wind turbine failure, Z represents a state value of random variable, and z.sup.wp.sup.i,k.sup.1 represents an output value of wind turbine i is wp.sub.i,k.sup.1.

[0061] S2, Influence of wind speed on wind turbine failure rate is analyzed, and a wind turbine failure model is established considering the wind turbine time-varying failure rate.

[0062] The S2 specifically includes: [0063] (1) the influence of wind speed on wind turbine failure rate is analyzed, and a wind turbine time-varying failure model is established.

[0064] As the wind turbine failure in the virtual power plant is closely related to the wind speed, the application considers the influence of the wind speed change on the wind turbine failure rate, and the wind turbine time-varying failure rate consists of basic failure rate of the wind turbine and variable failure rate of the wind turbine caused by the wind speed, as shown in the following formula:

.sub.I(t)=.sub.i,0+.sub.i,s(t), [0065] where .sub.I (t) represents the wind turbine time-varying failure rate i at time t, .sub.i,0 represents basic failure rate of wind turbine i, and .sub.i,s(t) represents variable failure rate of wind turbine i caused by wind speed at time t.

[0066] A relationship model between the variable failure rate of wind turbine i caused by wind speed at time t and wind speed s(t) is as follows:

.sub.i,s(t)=(.sub.i,maxs(t).sup.2.sub.i,mins(t).sup.2)/(s.sub.i.sup.co2s.sub.i.sup.ci2)+c.sub.s, [0067] where .sub.i,max represents a wind turbine failure rate corresponding to the cut-out wind speed s.sub.i.sup.co of the wind turbine i, .sub.i,min represents a wind turbine failure rate corresponding to the cut-in wind speed S.sub.i.sup.ci of the wind turbine i, and c.sub.s represents constants related to the cut-in wind speed and the cut-out wind speed. It can be seen from the above formula that the greater the wind speed, the greater the wind turbine failure rate. [0068] (2) a wind turbine failure model considering the wind turbine time-varying failure rate is established.

[0069] According to the application, the influence of the change of wind speed on the failure rate of the wind turbine is considered, and a multi-state model is adopted to describe the failure rate of the wind turbine at different wind speeds. Accordingly, the failure rate .sub.i (t) of the wind turbine i is divided into K.sub.s states according to divided s(t) states of wind speed s(t), and the failure rate of the k.sub.s(k.sub.s=1, . . . , K.sub.s)-th state is .sub.i,k and .sub.i,k<.sub.i,k+1. The failure model of the wind turbine usually adopts two-state model. The two-state refers to normal operation state and complete failure state. Considering the wind turbine time-varying failure rate .sub.i,k under different wind speeds, the failure probability of the wind turbine under the k.sub.s(k.sub.s=1, . . . , K.sub.s)-th wind speed state is obtained as follows:

q.sub.i,k.sup.wt(t)=1e.sup.-.sup.i,k.sup.t, [0070] where q.sub.i,k.sup.wt(t) represents the failure probability of wind turbine i at k.sub.s-th wind speed at time t.

[0071] Based on the failure probability of the wind turbine, the failure model of the wind turbine i in the k.sub.s-th wind speed state is established by using the improved general generating function method:

u.sub.i.sup.2(z,t)=(1q.sub.i,k.sup.wt(t).Math.z.sup.1+q.sub.i,k.sup.wt(t).Math.z.sup.0=e.sup.-.sup.i,k.sup.t.Math.z.sup.1+(1e.sup.-.sup.i,k.sup.t).Math.z.sup.0, [0072] where u.sub.i.sup.2(z,t) represents the improved general generating function representation method of failure model of wind turbine i considering the influence of wind speed on wind turbine failure probability, z.sup.1 indicates that wind turbine i is in normal operation and z.sup.0 indicates wind turbine i is in failure.

[0073] S3, On the basis of S1 and S2, a multi-state output model of distributed wind power considering wind speed and wind turbine time-varying failure rate is established.

[0074] The S3 specifically includes: [0075] based on the multi-state model of wind turbine output without considering wind turbine failure in Si and the wind turbine failure model considering wind turbine time-varying failure rate in S2, a multi-state model of output of the wind turbine i considering wind speed and wind turbine time-varying failure rate is established by using the improved general generating function method and is expressed by the following formula:

[00012]$\begin{array}{c}{u}_i^w\left(z,t\right)={}_{\mathrm{ser}}\left\{u_i^1\left(z,t\right),u_i^2\left(z,t\right)\right\}\\ ={}_{\mathrm{ser}}\left\{\underset{k_s=1}{\overset{K_s}{.Math.}}{q}_{i,k}\left(t\right).Math.z^{{\mathrm{wp}}_{i,k}^1},\left(1-q_{i,k}^{\mathrm{wt}}\left(t\right)\right).Math.z^1+q_{i,k}^{\mathrm{wt}}\left(t\right).Math.z^0\right\}\\ ={}_{\mathrm{ser}}\left\{\underset{k_s=1}{\overset{K_s}{.Math.}}{q}_{i,k}\left(t\right).Math.z^{{\mathrm{wp}}_{i,k}^1},e^{-{}_{i,k}t}.Math.z^1+\left(1-e^{-{}_{i,k}t}\right).Math.z^0\right\}\\ =\underset{k_s=1}{\overset{K_s}{.Math.}}{q}_{i,k}\left(t\right).Math.\left(1-q_{i,k}^{\mathrm{wt}}\left(t\right)\right).Math.z^{{\mathrm{wp}}_{i,k}^1}+\underset{k_s=1}{\overset{K_s}{.Math.}}{q}_{i,k}\left(t\right).Math.q_{i,k}^{\mathrm{wt}}\left(t\right).Math.z^0\\ =\underset{j=1}{\overset{n_w}{.Math.}}{q}_{i,j}^w\left(t\right).Math.z^{{\mathrm{wp}}_{i,j}^1}\end{array},$ [0076] where u.sub.i.sup.w(z,t) represents the improved general generating function representation method of wind turbine i output model considering wind speed and wind turbine time-varying failure rate, .sub.ser represents series operator, q.sub.i,j.sup.w(t) represents a probability of wind turbine i in state j is represented, and z.sup.wp.sup.i,j.sup.w represents the output value of wind turbine i in state j is wp.sub.i,j.sup.w.

[0077] S4, on the basis of S3, a multi-state output model of a virtual power plant composed of a plurality of distributed wind power is established.

[0078] The S4 specifically include:

[0079] For N.sub.w independent wind turbines in the virtual power plant, the output of N.sub.w wind turbines is shown in the following formula:

[00013]$\begin{array}{c}{U}^{\mathrm{VPP}}\left(z,t\right)={}_{\mathrm{par}}\left\{u_1^w\left(z,t\right),.Math.,u_i^w\left(z,t\right),.Math.,u_{N_w}^w\left(z,t\right)\right\}\\ =\underset{i=1}{\overset{N_w}{.Math.}}\underset{{j=1}^{}}{\overset{n_w}{.Math.}}{q}_{i,j}^w\left(t\right).Math.z^{{\mathrm{wp}}_j^w}\\ =\underset{m=1}{\overset{N_{\mathrm{vpp}}}{.Math.}}{q}_m^{\mathrm{vpp}}\left(t\right).Math.z^{{\mathrm{VPP}}_m}\end{array},$ [0080] where u.sup.VPP(z,t) represents an improved general generating function representation method for the output model of virtual power plant by aggregating N.sub.w independent wind turbines, .sub.par represents parallel operator, q.sub.m.sup.vpp(t) represents a probability of virtual power plant in state m, and z.sup.VPP.sup.m represents that the output value of virtual power plant in state m is VPP.sub.m.

[0081] S5, operation risk indicators of the virtual power plant are calculated according to the multi-state output model of the virtual power plant established in S4.

[0082] The S5 specifically includes: [0083] according to the multi-state output model of the virtual power plant obtained in S4, the operation risk indicators of the virtual power plant are calculated. The operation risk indicators include the power supply shortage probability D(t) , the expected power supply shortage E(t) and the power supply shortage loss A(t) of industrial users. The specific calculation formula is as follows:

[00014]$D\left(t\right)=\underset{m}{.Math.}{q}_m^{\mathrm{vpp}}\left(t\right),{\mathrm{VPP}}_{m}L,$$E\left(t\right)=\underset{w}{.Math.}{q}_m^{\mathrm{vpp}}\left(t\right).Math.{\mathrm{VPP}}_m,$$A\left(T\right)=\underset{t=0}{\overset{T}{.Math.}}\left(L-{\mathrm{VPP}}_m\right)q_m^{\mathrm{vpp}}\left(t\right)\mathrm{cdf}\left(\right),{\mathrm{VPP}}_m<L,$ [0084] where D(t) is power supply shortage probability of the virtual power plant, E(t) is expected power supply shortage of the virtual power plant and A(t) is power supply shortage loss for industrial users in the virtual power plant; L represents load value of industrial users powered by the virtual power plant, T represents total power supply time of the virtual power plant, and t represents the time and t [0,T]; indicates duration of power outage, cdf() indicates loss function of power supply shortage of industrial users, and is related to the duration of power outage.

[0085] Embodiments of the present application are as follows.

[0086] The virtual power plant in the embodiment is composed of 10 distributed wind turbines of 2 Million Watt (MW) and 2 industrial users with electricity demand of 5 MW. The schematic diagram of the virtual power plant structure model shown in FIG. 2 is established. The multi-state model of wind turbine output in Si without considering wind turbine failure is built, and the wind turbine output value and the time-varying probability value at different output values are gotten. Secondly, the wind turbine failure model considering the wind turbine time-varying failure rate in S2 is established, and the failure probability of the wind turbine under the influence of different wind speeds is obtained. Then the multi-state model of single independent wind turbine output considering wind speed and wind turbine time-varying failure rate in S3 is established by using the improved general generating function method and the output model of the virtual power plant composed of 8 distributed wind turbines is established by using the method in S4. Finally, when the virtual power plant runs for 100 hours, the operation risk indicators of the virtual power plant are calculated. The operation risk indicators of the virtual power plant include the power supply shortage probability, the expected power shortage, and the power supply shortage loss of industrial users, in which the loss function of power supply shortage of industrial users is expressed by a piecewise function. The unit losses of industrial users in virtual power plant with different power outage durations are shown in Table 1.

TABLE-US-00001 TABLE 1 Unit loss Duration of power outage (hours) (RMB/MWh) 0.01 0.5 1 4 8 Industrial 16.8 3.5 2.2 1.2 0.9 user 1 Industrial 96.5 22.6 15.3 13.0 10.6 user 2

[0087] From the above steps, it can be calculated that the power supply shortage probability in a virtual power plant operation risk indicator system whether to consider the time-varying failure rate of distributed wind power is shown in FIG. 3. After 100-hours operation, the expected power supply shortage of virtual power plant without considering the time-varying failure rate of distributed wind and with considering the time-varying failure rate of distributed wind is 422.93 MWh and 435.98 MWh, respectively. The power supply shortage losses of industrial users of virtual power plant without considering the time-varying failure rate of distributed wind power and with considering the time-varying failure rate of distributed wind power are 3700.66 Chinese yuan and 3814.85 Chinese yuan, respectively. Through comparison, it can be found that the power supply shortage probability, expected power shortage and power shortage loss of virtual power plant considering the time-varying failure rate of distributed wind power are all higher than the operation risk indicators of virtual power plant without considering the time-varying failure rate of distributed wind power. When t=100 hours, the power supply shortage probability of virtual power plant considering the time-varying failure rate of distributed wind power is 6.28% higher than that of virtual power plant without considering the time-varying failure rate of distributed wind power. The expected power shortage and power shortage loss of virtual power plant considering the time-varying failure rate of distributed wind power are 3.09% higher than those of virtual power plant without considering the time-varying failure rate of distributed wind power. To sum up, considering the wind turbine time-varying failure rate has certain influence on the operation risk of the virtual power plant, and quantitatively analyzing the influence of the wind turbine time-varying failure rate on the operation risk of the virtual power plant provides reference for the operation of the virtual power plant.

[0088] In this embodiment of the application, a virtual power plant operation risk analysis device considering time-varying failure rate of distributed wind power is also constructed, as shown in FIG. 4 and mainly includes a wind speed and wind turbine output module 10, a wind turbine time-varying failure rate acquisition module 20, a wind turbine failure probability acquisition module 30, a wind turbine output module 40 considering wind speed and wind turbine time-varying failure rate, and a virtual power plant operation risk assessment module 50.

[0089] The wind speed and wind turbine output module is used to construct the relationship model between wind turbine output and wind speed when the wind turbine is running well, to divide the wind speed into a plurality of states, and establish a multi-state wind speed model, and to calculate the wind turbine output value and the corresponding probability value without considering the wind turbine failure according to the relationship model between wind turbine output and wind speed and the multi-state output model of the wind turbine.

[0090] The wind turbine time-varying failure rate acquisition module is used for acquiring the variable failure rate of the wind turbine according to the relationship model between the variable failure rate of the wind turbine caused by wind speed and wind speed, and obtaining the wind turbine time-varying failure rate by adding the basic failure rate of the wind turbine and the variable failure rate of the wind turbine caused by the wind speed.

[0091] The wind turbine output module considering wind speed and wind turbine time-varying failure rate is used to build a wind turbine failure model considering wind turbine time-varying failure rate based on the wind turbine time-varying failure rate acquisition module, and to obtain the wind turbine output value and the corresponding probability value considering the wind speed and the wind turbine time-varying failure rate based on the wind turbine output obtained by the wind speed and wind turbine output module.

[0092] The virtual power plant operation risk assessment module is used to construct a output model of virtual power plant output model including a plurality of distributed wind power, to establish a virtual power plant operation risk indicator system including power supply shortage probability, expected power supply shortage and power supply shortage loss of industrial users in virtual power plant, and to calculate the power supply shortage probability, expected power supply shortage and power supply shortage loss of industrial users in virtual power plant.

[0093] It should be understood by those skilled in the art that the embodiments of the present application can provide methods, systems, or computer program products. Therefore, this application can take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product implemented on one or more computer usable storage media (including but not limited to disk storage, Compact Disc Read-Only Memory (CD-ROM), optical storage, etc.). The computer usable storage media contain computer usable program codes.

[0094] The application is described with reference to flowcharts and/or block diagrams of methods, devices (systems), and computer program products according to embodiments of the application. It should be understood that each flow and/or block in flowchart and/or block diagram, and combinations of flows and/or blocks in flowchart and/or block diagram can be realized by computer program instructions. These computer program instructions may be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing equipment to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing equipment produce a device for implementing the functions specified in one or more flow charts and/or one or more blocks of the block diagram.

[0095] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing equipment to work in a specific way, so that the instructions stored in the computer-readable memory produce an article of manufacture including instruction devices that implement the functions specified in one or more flow charts and/or one or more blocks of the block diagrams.

[0096] These computer program instructions may also be loaded on a computer or other programmable data processing equipment, so that a series of operation steps are executed on the computer or other programmable equipment to produce a computer-implemented process, so that the instructions executed on the computer or other programmable equipment provide steps for realizing the functions specified in one or more flows of the flowchart and/or one or more blocks of the block diagram.

[0097] The above-mentioned embodiments only describe the preferred mode of the application, but do not limit the scope of the application. On the premise of not departing from the design spirit of the application, all kinds of modifications and improvements made by ordinary technicians in the field to the technical scheme of the application shall fall within the scope of protection determined by the claims of the application.