CN106786784A - A kind of wind-powered electricity generation group of planes Poewr control method and system - Google Patents

Abstract

This application discloses a kind of wind-powered electricity generation group of planes Poewr control method and system, the method includes：The population adaptable with the scale of a wind-powered electricity generation group of planes is pre-created and initialized, initial population is obtained；Formula is determined using penalty function, it is determined that per the penalty function of typhoon group of motors；Using target function value computing formula, each individual target function value is calculated；Based on differential evolution algorithm, and each individual target function value is combined, Evolution of Population treatment is carried out to initial population, untill convergence, obtain global optimum individual；It is individual according to the global optimum, corresponding control is launched to the wind-powered electricity generation group of planes, to control the propeller pitch angle of every typhoon group of motors and the frequency of frequency-changing AC bus.The application can reduce the wind energy loss amount during wind turbine group control, and the power output of Wind turbines can be avoided to exceed power-handling capability.

Description

A kind of wind-powered electricity generation group of planes Poewr control method and system

Technical field

The present invention relates to technical field of wind power generation, more particularly to a kind of wind-powered electricity generation group of planes Poewr control method and system.

Background technology

Coastal waters wind-powered electricity generation has the advantages that wind speed is high, turbulence intensity is small, wind speed and direction is stable, is the main of wind-powered electricity generation industry development Trend.According to wind turbine group's centralized Control substitution discrete control of unit, fault rate is advantageously reduced, so as to reduce marine maintenance Work, increase effective generating dutation, while system synthesis sheet can be reduced.However, when wind speed profile is uneven in wind power plant, wind-powered electricity generation Group of planes centralized Control is easily caused the relatively strong winds loss of energy, and power output easily exceeds power-handling capability, how uneven in wind speed profile In the case of, reduce the wind energy loss amount during wind turbine group control and avoid the power output of Wind turbines from exceeding specified work( Rate value is the problem for needing further solution at present.

The content of the invention

In view of this, it is an object of the invention to provide a kind of wind-powered electricity generation group of planes Poewr control method and system, in wind speed point In the case of cloth inequality, the wind energy loss amount during wind turbine group control is reduced, and the defeated of Wind turbines can be avoided Go out power beyond power-handling capability.Its concrete scheme is as follows：

A kind of wind-powered electricity generation group of planes Poewr control method, is applied to a wind-powered electricity generation group of planes, wherein, per typhoon motor in the wind-powered electricity generation group of planes Rotating speed of the group in power control process is consistent, and methods described includes：

The population adaptable with the scale of the wind-powered electricity generation group of planes is pre-created and initialized, initial population is obtained；Wherein, institute The vector expression for stating i-th individuality in initial population is：

x_i,0=[x_1,i,0,x_2,i,0,…x_D,i,0,x_bus,i,0], i=1,2 ..., N；

In formula, D represents the quantity of the wind-powered electricity generation group of planes apoplexy group of motors, and N represents the quantity of individuality in the initial population, x_j,i,0Represent the propeller pitch angle of the jth typhoon group of motors in i-th individuality of the initial population, j=1,2 ..., D, x_bus,i,0 Represent the frequency of frequency-changing AC bus in i-th individuality of the initial population；

Formula is determined using penalty function, it is determined that per the penalty function of typhoon group of motors；Wherein, the penalty function determines formula For：

In formula, P_j,iRepresent the power output of the jth typhoon group of motors in i-th individuality, P_NRepresent the specified of Wind turbines Performance number, P_PENj,iRepresent the penalty function of the jth typhoon group of motors in i-th individuality；

Using target function value computing formula, each individual target function value is calculated；Wherein, the target function value meter Calculating formula is：

In formula, P_iI-th target function value of individuality is represented, C represents default penalty factor；

Based on differential evolution algorithm, and each individual target function value is combined, population is carried out to the initial population and is entered Change is processed, and untill convergence, obtains global optimum individual；Wherein, the individual vector expression of the global optimum is：

x_g=[x_1g,x_2g,…x_Dg,x_busg]；

In formula, x_jgRepresent the optimal propeller pitch angle of jth typhoon group of motors in the wind-powered electricity generation group of planes, x_busgRepresent that frequency-changing AC is female The optimal frequency of line；

It is individual according to the global optimum, corresponding control is launched to the wind-powered electricity generation group of planes, to control every typhoon group of motors Propeller pitch angle and frequency-changing AC bus frequency.

Optionally, the process for being pre-created and initializing the population adaptable with the scale of a wind-powered electricity generation group of planes, including：

Using initialization of population formula, the population adaptable with the scale of the wind-powered electricity generation group of planes is pre-created and initializes, Obtain the initial population；Wherein, the initialization of population formula is：

In formula, rand [0,1] represents the uniform random number for producing between zero and one, x_j ^URepresent that the propeller pitch angle of Wind turbines is adjusted The higher limit of section value, x_j ^LRepresent the lower limit of the propeller pitch angle regulated value of Wind turbines, x_bus ^URepresent the frequency of frequency-changing AC bus The higher limit of regulated value, x_bus ^LRepresent the lower limit of the frequency adjustment value of frequency-changing AC bus.

Optionally, the utilization penalty function determines formula, it is determined that the process of the penalty function per typhoon group of motors, including：

Formula is calculated using power output, the power output per typhoon group of motors is calculated；Wherein, the power output is calculated Formula is：

In formula, P_j,iThe power output of the jth typhoon group of motors in i-th individuality is represented, ρ represents atmospheric density, and A is represented Blade swept area, C_j(V_j,x_j,i,x_bus,i) represent jth typhoon group of motors power coefficient, V_jRepresent jth typhoon group of motors Real-time wind speed, x_j,iRepresent the propeller pitch angle of the jth typhoon group of motors in i-th individuality, x_bus,iFrequency conversion is handed in representing i-th individuality Flow the frequency of bus；

Power output and the penalty function using every typhoon group of motors determine formula, determine every typhoon group of motors Penalty function.

Optionally, the power coefficient C of jth typhoon group of motors_j(V_j,x_j,i,x_bus,i) computing formula be：

Wherein, index S is：

In formula, R represents the blade radius of Wind turbines, n_pRepresent motor number of pole-pairs, n_gRepresent gear-box no-load voltage ratio, K₁To K₆ To be in advance based on the coefficient of Wind turbines aerofoil profile determination.

Optionally, the wind-powered electricity generation group of planes be using the central controlled wind-powered electricity generation group of planes of high voltage direct current transmission converter, or use The frequency dividing central controlled wind-powered electricity generation group of planes of ac transmission current transformer, or use the central controlled wind-powered electricity generation group of planes of variable frequency transformer.

Optionally, it is described based on differential evolution algorithm, and each individual target function value is combined, to the initial population Evolution of Population treatment is carried out, untill convergence, the individual process of global optimum is obtained, including：

Step A1：The minimum individuality of target function value is filtered out from the initial population, and the individuality is defined as to work as Preceding optimum individual x_History；

Step A2：Using the formula that makes a variation, to current population Y_GIn each individuality carry out mutation operation, produce corresponding change Different individuality；Wherein, G represents the algebraically numbering of population, and the variation formula is：

v_i,G+1=x_r1,G+F*(x_r2,G-x_r3,G)；

In formula, v_i,G+1Represent and current population Y_GIn i-th individual corresponding variation it is individual, x_r1,G、x_r2,GAnd x_r3,G Represent in current population Y_GIn the r1 that randomly selects, r2 and the r3 individuality, and r1 ≠ r2 ≠ r3 ≠ i, F represent Mutation operator；

Step A3：Crossover operation is carried out using crossing formula, obtains testing vectorial u_i,G+1；

Wherein, u_i,G+1=(u_1,i,G+1,u_2,i,G+1,...u_D,i,G+1,u_bus,i,G+1), the crossing formula is：

In formula, rand (j) represent between zero and one it is random be uniformly distributed, CR represents crossover operator, rnb (i) represent [1, 2 ..., D+1] in a randomly selected sequence；

Step A4：Using individual replacement formula, the individuality to not meeting edge-restraint condition in the experiment vector is carried out Update and replace, obtain interim population；

Wherein, the edge-restraint condition includes：x_j ^L≤u_j,i,G+1≤x_j ^U, x_bus ^L≤u_bus,i,G+1≤x_bus ^U；The individuality Replacement formula is：

In formula, rand [0,1] represents the uniform random number for producing between zero and one, x_j ^URepresent that the propeller pitch angle of Wind turbines is adjusted The higher limit of section value, x_j ^LRepresent the lower limit of the propeller pitch angle regulated value of Wind turbines, x_bus ^URepresent the frequency of frequency-changing AC bus The higher limit of regulated value, x_bus ^LRepresent the lower limit of the frequency adjustment value of frequency-changing AC bus；

Step A5：Calculate each individual target function value in interim population, and by target function value in interim population compared with Big individuality accordingly replaces current population Y_GThe middle less individuality of target function value, the population Y after currently being updated_G+1；

Step A6：From current population Y_G+1In filter out the maximum individuality of target function value, and judge the individual target letter Whether numerical value is more than current optimum individual x_HistoryTarget function value, if it is, using the individuality to current optimum individual x_HistoryIt is replaced renewal；

Step A7：Reenter step A2, untill evolutionary generation is restrained, and by evolve terminate after obtain it is current most Excellent individual x_HistoryIt is defined as the global optimum individual.

The invention also discloses a kind of wind-powered electricity generation group of planes power control system, a wind-powered electricity generation group of planes is applied to, wherein, the wind turbine It is consistent per rotating speed of the typhoon group of motors in power control process in group, the system includes：

Initialization of population module, for being pre-created and initializes the population adaptable with the scale of the wind-powered electricity generation group of planes, Obtain initial population；Wherein, the vector expression of i-th individuality in the initial population is：

x_i,0=[x_1,i,0,x_2,i,0,…x_D,i,0,x_bus,i,0], i=1,2 ..., N；

In formula, D represents the quantity of the wind-powered electricity generation group of planes apoplexy group of motors, and N represents the quantity of individuality in the initial population, x_j,i,0Represent the propeller pitch angle of the jth typhoon group of motors in i-th individuality of the initial population, j=1,2 ..., D, x_bus,i,0 Represent the frequency of frequency-changing AC bus in i-th individuality of the initial population；

Penalty function determining module, for determining formula using penalty function, it is determined that per the penalty function of typhoon group of motors；Wherein, The penalty function determines that formula is：

In formula, P_j,iRepresent the power output of the jth typhoon group of motors in i-th individuality, P_NRepresent the specified of Wind turbines Performance number, P_PENj,iRepresent the penalty function of the jth typhoon group of motors in i-th individuality；

Target function value computing module, for utilizing target function value computing formula, calculates each individual object function Value；Wherein, the target function value computing formula is：

In formula, P_iI-th target function value of individuality is represented, C represents default penalty factor；

Evolution of Population processing module, for based on differential evolution algorithm, and combines each individual target function value, to institute Stating initial population carries out Evolution of Population treatment, untill convergence, obtains global optimum individual；Wherein, the global optimum The vector expression of body is：

x_g=[x_1g,x_2g,…x_Dg,x_busg]；

In formula, x_jgRepresent the optimal propeller pitch angle of jth typhoon group of motors in the wind-powered electricity generation group of planes, x_busgRepresent that frequency-changing AC is female The optimal frequency of line；

Wind turbine group control module, for individual according to the global optimum, corresponding control is launched to the wind-powered electricity generation group of planes System, to control the propeller pitch angle of every typhoon group of motors and the frequency of frequency-changing AC bus.

Optionally, the initialization of population module, specifically for utilizing initialization of population formula, is pre-created and initializes The population adaptable with the scale of the wind-powered electricity generation group of planes, obtains the initial population；Wherein, the initialization of population formula is：

In formula, rand [0,1] represents the uniform random number for producing between zero and one, x_j ^URepresent that the propeller pitch angle of Wind turbines is adjusted The higher limit of section value, x_j ^LRepresent the lower limit of the propeller pitch angle regulated value of Wind turbines, x_bus ^URepresent the frequency of frequency-changing AC bus The higher limit of regulated value, x_bus ^LRepresent the lower limit of the frequency adjustment value of frequency-changing AC bus.

Optionally, the penalty function determining module, including：

Power output calculating sub module, for calculating formula using power output, calculates the output work per typhoon group of motors Rate；Wherein, the power output calculating formula is：

In formula, P_j,iThe power output of the jth typhoon group of motors in i-th individuality is represented, ρ represents atmospheric density, and A is represented Blade swept area, C_j(V_j,x_j,i,x_bus,i) represent jth typhoon group of motors power coefficient, V_jRepresent jth typhoon group of motors Real-time wind speed, x_j,iRepresent the propeller pitch angle of the jth typhoon group of motors in i-th individuality, x_bus,iFrequency conversion is handed in representing i-th individuality Flow the frequency of bus；

Penalty function determination sub-module, it is public for being determined using the power output and the penalty function per typhoon group of motors Formula, determines the penalty function of every typhoon group of motors.

Optionally, the power output calculating sub module, including：

Power coefficient computing unit, for utilizing power coefficient computing formula, calculates the power coefficient per typhoon group of motors； Wherein, the power coefficient computing formula is：

Wherein, index S is：

In formula, R represents the blade radius of Wind turbines, n_pRepresent motor number of pole-pairs, n_gRepresent gear-box no-load voltage ratio, K₁To K₆ To be in advance based on the coefficient of Wind turbines aerofoil profile determination.

In the present invention, wind-powered electricity generation group of planes Poewr control method, including：It is pre-created and initializes the scale phase with a wind-powered electricity generation group of planes The population of adaptation, obtains initial population；Formula is determined using above-mentioned penalty function, it is determined that per the penalty function of typhoon group of motors；Utilize Above-mentioned target function value computing formula, calculates each individual target function value；Based on differential evolution algorithm, and combine per each and every one The target function value of body, Evolution of Population treatment is carried out to the initial population, untill convergence, obtains global optimum individual； It is individual according to the global optimum, corresponding control is launched to the wind-powered electricity generation group of planes.

It can be seen that, the present invention reduces a wind-powered electricity generation group of planes to a certain extent using the differential evolution algorithm for incorporating penalty function Wind energy is lost, and in the case where the power output of Wind turbines exceeds power-handling capability, it can be ensured that Wind turbines are maintained In invariable power generating state, the purpose of maximal wind-energy output control of the wind-powered electricity generation group of planes in full blast speed scope is achieved in, thus It can be seen that, the present invention can reduce the wind energy loss amount during wind turbine group control, and can avoid the output of Wind turbines Power exceeds power-handling capability.

Brief description of the drawings

In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this Inventive embodiment, for those of ordinary skill in the art, on the premise of not paying creative work, can also basis The accompanying drawing of offer obtains other accompanying drawings.

Fig. 1 is a kind of wind-powered electricity generation group of planes Poewr control method flow chart disclosed in the embodiment of the present invention；

Fig. 2 is a kind of wind-powered electricity generation group of planes power control system structural representation disclosed in the embodiment of the present invention.

Specific embodiment

Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.It is based on Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made Embodiment, belongs to the scope of protection of the invention.

The embodiment of the invention discloses a kind of wind-powered electricity generation group of planes Poewr control method, a wind-powered electricity generation group of planes is applied to, wherein, wind turbine It is consistent per rotating speed of the typhoon group of motors in power control process in group, shown in Figure 1, the method includes：

Step S11：The population adaptable with the scale of a wind-powered electricity generation group of planes is pre-created and initialized, initial population is obtained；Its In, the vector expression of i-th individuality in initial population is：

x_i,0=[x_1,i,0,x_2,i,0,…x_D,i,0,x_bus,i,0], i=1,2 ..., N；

In formula, D represents the quantity of wind-powered electricity generation group of planes apoplexy group of motors, and N represents the quantity of individuality in initial population, x_j,i,0Represent The propeller pitch angle of the jth typhoon group of motors in i-th individuality of initial population, j=1,2 ..., D, x_bus,i,0Represent initial population I-th individuality in frequency-changing AC bus frequency；

Step S12：Formula is determined using penalty function, it is determined that per the penalty function of typhoon group of motors；Wherein, penalty function determines public Formula is：

In formula, P_j,iRepresent the power output of the jth typhoon group of motors in i-th individuality, P_NRepresent the specified of Wind turbines Performance number, P_PENj,iRepresent the penalty function of the jth typhoon group of motors in i-th individuality；

Step S13：Using target function value computing formula, each individual target function value is calculated；Wherein, object function Value computing formula is：

In formula, P_iI-th target function value of individuality is represented, C represents default penalty factor；

Step S14：Based on differential evolution algorithm, and each individual target function value is combined, initial population is planted Group's evolution treatment, untill convergence, obtains global optimum individual；Wherein, the individual vector expression of global optimum is：

x_g=[x_1g,x_2g,…x_Dg,x_busg]；

In formula, x_jgRepresent the optimal propeller pitch angle of jth typhoon group of motors in a wind-powered electricity generation group of planes, x_busgRepresent frequency-changing AC bus Optimal frequency；

Step S15：It is individual according to global optimum, corresponding control is launched to a wind-powered electricity generation group of planes.

Wherein, the detailed process of above-mentioned steps S15 can include：By x_busgThe frequency regulator of a wind-powered electricity generation group of planes is assigned to, with Make frequency regulator according to x_busgCarry out corresponding frequency regulation, and by x_1g,x_2g,…x_DgIt is assigned to the oar of corresponding Wind turbines Elongation controller, so that each award setting device carries out corresponding award setting according to the optimal propeller pitch angle for each getting Treatment.

It should be noted that the wind-powered electricity generation group of planes in the present embodiment is specifically as follows being concentrated using high voltage direct current transmission converter A wind-powered electricity generation group of planes for control, or using the frequency dividing central controlled wind-powered electricity generation group of planes of ac transmission current transformer, or use variable frequency transformation The central controlled wind-powered electricity generation group of planes of device.In the present embodiment, per rotating speed of the typhoon group of motors in power control process in a wind-powered electricity generation group of planes It is consistent, it is necessary to all Wind turbines in a wind-powered electricity generation group of planes are uniformly controlled in above-mentioned power control process, and Any Wind turbines will not individually be controlled.

It can be seen that, the embodiment of the present invention reduces wind-powered electricity generation to a certain extent using the differential evolution algorithm for incorporating penalty function The wind energy loss of a group of planes, and in the case where the power output of Wind turbines exceeds power-handling capability, it can be ensured that wind turbine Group maintains invariable power generating state, is achieved in the mesh of maximal wind-energy output control of the wind-powered electricity generation group of planes in full blast speed scope , it can be seen that, the embodiment of the present invention can reduce the wind energy loss amount during wind turbine group control, and can avoid wind The power output of group of motors exceeds power-handling capability.

The embodiment of the invention discloses a kind of specific wind-powered electricity generation group of planes Poewr control method, comprise the following steps：

Step S21：Using initialization of population formula, the kind adaptable with the scale of a wind-powered electricity generation group of planes is pre-created and initialized Group, obtains initial population；Wherein, the vector expression of i-th individuality in initial population is：

x_i,0=[x_1,i,0,x_2,i,0,…x_D,i,0,x_bus,i,0], i=1,2 ..., N；

In formula, D represents the quantity of wind-powered electricity generation group of planes apoplexy group of motors, and N represents the quantity of individuality in initial population, x_j,i,0Represent The propeller pitch angle of the jth typhoon group of motors in i-th individuality of initial population, j=1,2 ..., D, x_bus,i,0Represent initial population I-th individuality in frequency-changing AC bus frequency.

In the present embodiment, above-mentioned initialization of population formula is specially：

In formula, rand [0,1] represents the uniform random number for producing between zero and one, x_j ^URepresent that the propeller pitch angle of Wind turbines is adjusted The higher limit of section value, x_j ^LRepresent the lower limit of the propeller pitch angle regulated value of Wind turbines, x_bus ^URepresent the frequency of frequency-changing AC bus The higher limit of regulated value, x_bus ^LRepresent the lower limit of the frequency adjustment value of frequency-changing AC bus.

Step S22：Formula is calculated using power output, the power output per typhoon group of motors is calculated；Wherein, power output Computing formula is：

In formula, P_j,iThe power output of the jth typhoon group of motors in i-th individuality is represented, ρ represents atmospheric density, and A is represented Blade swept area, C_j(V_j,x_j,i,x_bus,i) represent jth typhoon group of motors power coefficient, V_jRepresent jth typhoon group of motors Real-time wind speed, x_j,iRepresent the propeller pitch angle of the jth typhoon group of motors in i-th individuality, x_bus,iFrequency conversion is handed in representing i-th individuality Flow the frequency of bus.

In the present embodiment, the power coefficient C of jth typhoon group of motors_j(V_j,x_j,i,x_bus,i) computing formula be：

Wherein, index S is：

In formula, R represents the blade radius of Wind turbines, n_pRepresent motor number of pole-pairs, n_gRepresent gear-box no-load voltage ratio, K₁To K₆ To be in advance based on the coefficient of Wind turbines aerofoil profile determination.

Step S23：Power output and penalty function using every typhoon group of motors determine formula, determine every typhoon motor The penalty function of group.

Wherein, penalty function determines that formula is：

In formula, P_j,iRepresent the power output of the jth typhoon group of motors in i-th individuality, P_NRepresent the specified of Wind turbines Performance number, P_PENj,iRepresent the penalty function of the jth typhoon group of motors in i-th individuality.

Step S24：Based on differential evolution algorithm, and each individual target function value is combined, initial population is planted Group's evolution treatment, untill convergence, obtains global optimum individual；Wherein, the individual vector expression of global optimum is：

x_g=[x_1g,x_2g,…x_Dg,x_busg]；

In formula, x_jgRepresent the optimal propeller pitch angle of jth typhoon group of motors in a wind-powered electricity generation group of planes, x_busgRepresent frequency-changing AC bus Optimal frequency.

In the present embodiment, the process of above-mentioned steps S24 can be seen and specifically include below step A1 to A7：

Step A1：The minimum individuality of target function value is filtered out from initial population, and the individuality is defined as currently most Excellent individual x_History；

Step A2：Using the formula that makes a variation, to current population Y_GIn each individuality carry out mutation operation, produce corresponding change Different individuality；Wherein, G represents the algebraically numbering of population, and variation formula is：

v_i,G+1=x_r1,G+F*(x_r2,G-x_r3,G)；

In formula, v_i,G+1Represent and current population Y_GIn i-th individual corresponding variation it is individual, x_r1,G、x_r2,GAnd x_r3,G Represent in current population Y_GIn the r1 that randomly selects, r2 and the r3 individuality, and r1 ≠ r2 ≠ r3 ≠ i, F represent Mutation operator；

Step A3：Crossover operation is carried out using crossing formula, obtains testing vectorial u_i,G+1；

Wherein, u_i,G+1=(u_1,i,G+1,u_2,i,G+1,…u_D,i,G+1,u_bus,i,G+1), crossing formula is：

In formula, rand (j) represent between zero and one it is random be uniformly distributed, CR represents crossover operator, rnb (i) represent [1, 2 ..., D+1] in a randomly selected sequence；

Step A4：Using individual replacement formula, the individuality to not meeting edge-restraint condition in experiment vector is updated Replace, obtain interim population；

Wherein, edge-restraint condition includes：x_j ^L≤u_j,i,G+1≤x_j ^U, x_bus ^L≤u_bus,i,G+1≤x_bus ^U；Individual replacement formula For：

In formula, rand [0,1] represents the uniform random number for producing between zero and one, x_j ^URepresent that the propeller pitch angle of Wind turbines is adjusted The higher limit of section value, x_j ^LRepresent the lower limit of the propeller pitch angle regulated value of Wind turbines, x_bus ^URepresent the frequency of frequency-changing AC bus The higher limit of regulated value, x_bus ^LRepresent the lower limit of the frequency adjustment value of frequency-changing AC bus；

Step A5：Calculate each individual target function value in interim population, and by target function value in interim population compared with Big individuality accordingly replaces current population Y_GThe middle less individuality of target function value, the population Y after currently being updated_G+1；

Step A6：From current population Y_G+1In filter out the maximum individuality of target function value, and judge the individual target letter Whether numerical value is more than current optimum individual x_HistoryTarget function value, if it is, using the individuality to current optimum individual x_HistoryIt is replaced renewal；

Step A7：Reenter step A2, untill evolutionary generation is restrained, and by evolve terminate after obtain it is current most Excellent individual x_HistoryIt is defined as global optimum individual.

Step S25：It is individual according to global optimum, corresponding control is launched to a wind-powered electricity generation group of planes.

Accordingly, the embodiment of the invention also discloses a kind of wind-powered electricity generation group of planes power control system, a wind-powered electricity generation group of planes is applied to, its In, it is consistent per rotating speed of the typhoon group of motors in power control process in a wind-powered electricity generation group of planes, shown in Figure 2, the system Including：

Initialization of population module 11, for being pre-created and initializes the population adaptable with the scale of a wind-powered electricity generation group of planes, obtains To initial population；Wherein, the vector expression of i-th individuality in initial population is：

x_i,0=[x_1,i,0,x_2,i,0,…x_D,i,0,x_bus,i,0], i=1,2 ..., N；

In formula, D represents the quantity of wind-powered electricity generation group of planes apoplexy group of motors, and N represents the quantity of individuality in initial population, x_j,i,0Represent The propeller pitch angle of the jth typhoon group of motors in i-th individuality of initial population, j=1,2 ..., D, x_bus,i,0Represent initial population I-th individuality in frequency-changing AC bus frequency；

Penalty function determining module 12, for determining formula using penalty function, it is determined that per the penalty function of typhoon group of motors；Its In, penalty function determines that formula is：

In formula, P_j,iRepresent the power output of the jth typhoon group of motors in i-th individuality, P_NRepresent the specified of Wind turbines Performance number, P_PENj,iRepresent the penalty function of the jth typhoon group of motors in i-th individuality；

Target function value computing module 13, for utilizing target function value computing formula, calculates each individual target letter Numerical value；Wherein, target function value computing formula is：

In formula, P_iI-th target function value of individuality is represented, C represents default penalty factor；

Evolution of Population processing module 14, for based on differential evolution algorithm, and combines each individual target function value, right Initial population carries out Evolution of Population treatment, untill convergence, obtains global optimum individual；Wherein, global optimum it is individual to Measuring expression formula is：

x_g=[x_1g,x_2g,…x_Dg,x_busg]；

In formula, x_jgRepresent the optimal propeller pitch angle of jth typhoon group of motors in a wind-powered electricity generation group of planes, x_busgRepresent frequency-changing AC bus Optimal frequency；

Wind turbine group control module 15, for individual according to global optimum, corresponding control is launched to a wind-powered electricity generation group of planes.

In the present embodiment, above-mentioned initialization of population module specifically can be used for, using initialization of population formula, being pre-created And the population adaptable with the scale of a wind-powered electricity generation group of planes is initialized, obtain initial population；Wherein, initialization of population formula is：

In formula, rand [0,1] represents the uniform random number for producing between zero and one, x_j ^URepresent that the propeller pitch angle of Wind turbines is adjusted The higher limit of section value, x_j ^LRepresent the lower limit of the propeller pitch angle regulated value of Wind turbines, x_bus ^URepresent the frequency of frequency-changing AC bus The higher limit of regulated value, x_bus ^LRepresent the lower limit of the frequency adjustment value of frequency-changing AC bus.

Further, the penalty function determining module in the present embodiment, specifically can include power output calculating sub module with And penalty function determination sub-module；Wherein,

Power output calculating sub module, for calculating formula using power output, calculates the output work per typhoon group of motors Rate；Wherein, power output calculating formula is：

In formula, P_j,iThe power output of the jth typhoon group of motors in i-th individuality is represented, ρ represents atmospheric density, and A is represented Blade swept area, C_j(V_j,x_j,i,x_bus,i) represent jth typhoon group of motors power coefficient, V_jRepresent jth typhoon group of motors Real-time wind speed, x_j,iRepresent the propeller pitch angle of the jth typhoon group of motors in i-th individuality, x_bus,iFrequency conversion is handed in representing i-th individuality Flow the frequency of bus；

Penalty function determination sub-module, for determining formula using the power output and penalty function per typhoon group of motors, really Make the penalty function of every typhoon group of motors.

In addition, above-mentioned power output calculating sub module, can specifically include：

Power coefficient computing unit, for utilizing power coefficient computing formula, calculates the power coefficient per typhoon group of motors； Wherein, power coefficient computing formula is：

Wherein, index S is：

In formula, R represents the blade radius of Wind turbines, n_pRepresent motor number of pole-pairs, n_gRepresent gear-box no-load voltage ratio, K₁To K₆ To be in advance based on the coefficient of Wind turbines aerofoil profile determination.

The corresponding contents disclosed in previous embodiment are may be referred on the more detailed course of work of above-mentioned modules, No longer repeated herein.

It can be seen that, the embodiment of the present invention reduces wind-powered electricity generation to a certain extent using the differential evolution algorithm for incorporating penalty function The wind energy loss of a group of planes, and in the case where the power output of Wind turbines exceeds power-handling capability, it can be ensured that wind turbine Group maintains invariable power generating state, is achieved in the mesh of maximal wind-energy output control of the wind-powered electricity generation group of planes in full blast speed scope , it can be seen that, the embodiment of the present invention can reduce the wind energy loss amount during wind turbine group control, and can avoid wind The power output of group of motors exceeds power-handling capability.

Finally, in addition it is also necessary to explanation, herein, such as first and second or the like relational terms be used merely to by One entity or operation make a distinction with another entity or operation, and not necessarily require or imply these entities or operation Between there is any this actual relation or order.And, term " including ", "comprising" or its any other variant meaning Covering including for nonexcludability, so that process, method, article or equipment including a series of key elements not only include that A little key elements, but also other key elements including being not expressly set out, or also include for this process, method, article or The intrinsic key element of equipment.In the absence of more restrictions, the key element limited by sentence "including a ...", does not arrange Except also there is other identical element in the process including the key element, method, article or equipment.

A kind of wind-powered electricity generation group of planes Poewr control method provided by the present invention and system are described in detail above, herein In apply specific case principle of the invention and implementation method be set forth, the explanation of above example is only intended to side Assistant solves the method for the present invention and its core concept；Simultaneously for those of ordinary skill in the art, according to think of of the invention Think, will change in specific embodiments and applications, in sum, it is right that this specification content should not be construed as Limitation of the invention.

Claims (10)

1. a kind of wind-powered electricity generation group of planes Poewr control method, it is characterised in that be applied to a wind-powered electricity generation group of planes, wherein, in the wind-powered electricity generation group of planes It is consistent per rotating speed of the typhoon group of motors in power control process, methods described includes：

The population adaptable with the scale of the wind-powered electricity generation group of planes is pre-created and initialized, initial population is obtained；Wherein, it is described first The vector expression of i-th individuality in beginning population is：

x_i,0=[x_1,i,0,x_2,i,0,…x_D,i,0,x_bus,i,0], i=1,2 ..., N；

In formula, D represents the quantity of the wind-powered electricity generation group of planes apoplexy group of motors, and N represents the quantity of individuality in the initial population, x_j,i,0 Represent the propeller pitch angle of the jth typhoon group of motors in i-th individuality of the initial population, j=1,2 ..., D, x_bus,i,0Represent The frequency of frequency-changing AC bus in i-th individuality of the initial population；

Formula is determined using penalty function, it is determined that per the penalty function of typhoon group of motors；Wherein, the penalty function determines that formula is：

In formula, P_j,iRepresent the power output of the jth typhoon group of motors in i-th individuality, P_NRepresent the rated power of Wind turbines Value, P_PENj,iRepresent the penalty function of the jth typhoon group of motors in i-th individuality；

Using target function value computing formula, each individual target function value is calculated；Wherein, the target function value calculates public Formula is：

$P_i=\underset{j=1}{\overset{D}{\Σ}}{P}_{j,i}-C\underset{j=1}{\overset{D}{\Σ}}{P}_{PENj,i};$

In formula, P_iI-th target function value of individuality is represented, C represents default penalty factor；

Based on differential evolution algorithm, and each individual target function value is combined, the initial population is carried out at Evolution of Population Reason, untill convergence, obtains global optimum individual；Wherein, the individual vector expression of the global optimum is：

x_g=[x_1g,x_2g,…x_Dg,x_busg]；

In formula, x_jgRepresent the optimal propeller pitch angle of jth typhoon group of motors in the wind-powered electricity generation group of planes, x_busgRepresent frequency-changing AC bus Optimal frequency；

It is individual according to the global optimum, corresponding control is launched to the wind-powered electricity generation group of planes, to control the oar of every typhoon group of motors The frequency of elongation and frequency-changing AC bus.

2. wind-powered electricity generation group of planes Poewr control method according to claim 1, it is characterised in that described to be pre-created and initialize The process of the population adaptable with the scale of a wind-powered electricity generation group of planes, including：

Using initialization of population formula, the population adaptable with the scale of the wind-powered electricity generation group of planes is pre-created and initialized, obtain The initial population；Wherein, the initialization of population formula is：

$\left\{\begin{array}{c}{x}_{j,i,0}=rand\[0,1\]({x_j}^U-{x_j}^L)+{x_j}^L\\ x_{bus,i,0}=rand\[0,1\]({x_{bus}}^U-{x_{bus}}^L)+{x_{bus}}^L\end{array}\right.;$

In formula, rand [0,1] represents the uniform random number for producing between zero and one, x_j ^URepresent the propeller pitch angle regulated value of Wind turbines Higher limit, x_j ^LRepresent the lower limit of the propeller pitch angle regulated value of Wind turbines, x_bus ^URepresent the frequency regulation of frequency-changing AC bus The higher limit of value, x_bus ^LRepresent the lower limit of the frequency adjustment value of frequency-changing AC bus.

3. wind-powered electricity generation group of planes Poewr control method according to claim 1, it is characterised in that the utilization penalty function determines public Formula, it is determined that the process of the penalty function per typhoon group of motors, including：

Formula is calculated using power output, the power output per typhoon group of motors is calculated；Wherein, the power output calculates formula For：

$P_{j,i}=\frac{1}{2}{\mathrm{\ρAC}}_j(V_j,x_{j,i},x_{bus,i}){V_j}^3;$

In formula, P_j,iThe power output of the jth typhoon group of motors in i-th individuality is represented, ρ represents atmospheric density, and A represents blade Swept area, C_j(V_j,x_j,i,x_bus,i) represent jth typhoon group of motors power coefficient, V_jRepresent the real-time of jth typhoon group of motors Wind speed, x_j,iRepresent the propeller pitch angle of the jth typhoon group of motors in i-th individuality, x_bus,iFrequency-changing AC is female in representing i-th individuality The frequency of line；

Power output and the penalty function using every typhoon group of motors determine formula, and that determines every typhoon group of motors penalizes letter Number.

4. wind-powered electricity generation group of planes Poewr control method according to claim 3, it is characterised in that the power of jth typhoon group of motors Coefficient C_j(V_j,x_j,i,x_bus,i) computing formula be：

$C_j\left(V_j,x_{j,i},x_{bus,i}\right)=K_1\left(\frac{K_2}{\frac{2{\mathrm{\πx}}_{bus,i}R}{n_p{n}_g{V}_j}+0.08x_{j,i}}-\frac{0.035K_2}{{x_{j,i}}^3+1}-K_3{x}_{j,i}-K_4\right)e^S+K_5\frac{2{\mathrm{\πx}}_{bus,i}R}{n_p{n}_g{V}_j};$

Wherein, index S is：

$S=\frac{-K_6}{\frac{2{\mathrm{\πx}}_{bus,i}R}{n_p{n}_g{V}_j}+0.008x_{j,i}}+\frac{0.035K_6}{{x_{j,i}}^3+1};$

In formula, R represents the blade radius of Wind turbines, n_pRepresent motor number of pole-pairs, n_gRepresent gear-box no-load voltage ratio, K₁To K₆For advance Based on the coefficient that Wind turbines aerofoil profile determines.

5. wind-powered electricity generation group of planes Poewr control method according to claim 1, it is characterised in that the wind-powered electricity generation group of planes is using height The central controlled wind-powered electricity generation group of planes of straightening stream current transmission converter, or using the frequency dividing central controlled wind turbine of ac transmission current transformer Group, or use the central controlled wind-powered electricity generation group of planes of variable frequency transformer.

6. the wind-powered electricity generation group of planes Poewr control method according to any one of claim 1 to 5, it is characterised in that described based on difference Divide evolution algorithm, and combine each individual target function value, Evolution of Population treatment is carried out to the initial population, until convergence Untill, the individual process of global optimum is obtained, including：

Step A1：The minimum individuality of target function value is filtered out from the initial population, and the individuality is defined as currently most Excellent individual x_History；

Step A2：Using the formula that makes a variation, to current population Y_GIn each individuality carry out mutation operation, produce corresponding variation Body；Wherein, G represents the algebraically numbering of population, and the variation formula is：

v_i,G+1=x_r1,G+F*(x_r2,G-x_r3,G)；

In formula, v_i,G+1Represent and current population Y_GIn i-th individual corresponding variation it is individual, x_r1,G、x_r2,GAnd x_r3,GRepresent In current population Y_GIn the r1 that randomly selects, r2 and the r3 individuality, and r1 ≠ r2 ≠ r3 ≠ i, F represent variation Operator；

Step A3：Crossover operation is carried out using crossing formula, obtains testing vectorial u_i,G+1；

Wherein, u_i,G+1=(u_1,i,G+1,u_2,i,G+1,...u_D,i,G+1,u_bus,i,G+1), the crossing formula is：

In formula, rand (j) represent between zero and one it is random be uniformly distributed, CR represents crossover operator, rnb (i) represent [1, 2 ..., D+1] in a randomly selected sequence；

Step A4：Using individual replacement formula, the individuality to not meeting edge-restraint condition in the experiment vector is updated Replace, obtain interim population；

Wherein, the edge-restraint condition includes：x_j ^L≤u_j,i,G+1≤x_j ^U, x_bus ^L≤u_bus,i,G+1≤x_bus ^U；The individual replacement Formula is：

$\left\{\begin{array}{c}{x}_{j,i,G+1}=rand\[0,1\]({x_j}^U-{x_j}^L)+{x_j}^L\\ x_{bus,i,G+1}=rand\[0,1\]({x_{bus}}^U-{x_{bus}}^L)+{x_{bus}}^L\end{array}\right.;$

In formula, rand [0,1] represents the uniform random number for producing between zero and one, x_j ^URepresent the propeller pitch angle regulated value of Wind turbines Higher limit, x_j ^LRepresent the lower limit of the propeller pitch angle regulated value of Wind turbines, x_bus ^URepresent the frequency regulation of frequency-changing AC bus The higher limit of value, x_bus ^LRepresent the lower limit of the frequency adjustment value of frequency-changing AC bus；

Step A5：Each individual target function value in interim population is calculated, and target function value in interim population is larger Individuality accordingly replaces current population Y_GThe middle less individuality of target function value, the population Y after currently being updated_G+1；

Step A6：From current population Y_G+1In filter out the maximum individuality of target function value, and judge the individual target function value Whether current optimum individual x is more than_HistoryTarget function value, if it is, using the individuality to current optimum individual x_History It is replaced renewal；

Step A7：Reenter step A2, untill evolutionary generation is restrained, and by evolve terminate after obtain it is current optimal Body x_HistoryIt is defined as the global optimum individual.

7. a kind of wind-powered electricity generation group of planes power control system, it is characterised in that be applied to a wind-powered electricity generation group of planes, wherein, in the wind-powered electricity generation group of planes It is consistent per rotating speed of the typhoon group of motors in power control process, the system includes：

Initialization of population module, for being pre-created and initializes the population adaptable with the scale of the wind-powered electricity generation group of planes, obtains Initial population；Wherein, the vector expression of i-th individuality in the initial population is：

x_i,0=[x_1,i,0,x_2,i,0,…x_D,i,0,x_bus,i,0], i=1,2 ..., N；

In formula, D represents the quantity of the wind-powered electricity generation group of planes apoplexy group of motors, and N represents the quantity of individuality in the initial population, x_j,i,0 Represent the propeller pitch angle of the jth typhoon group of motors in i-th individuality of the initial population, j=1,2 ..., D, x_bus,i,0Represent The frequency of frequency-changing AC bus in i-th individuality of the initial population；

Penalty function determining module, for determining formula using penalty function, it is determined that per the penalty function of typhoon group of motors；Wherein, it is described Penalty function determines that formula is：

In formula, P_j,iRepresent the power output of the jth typhoon group of motors in i-th individuality, P_NRepresent the rated power of Wind turbines Value, P_PENj,iRepresent the penalty function of the jth typhoon group of motors in i-th individuality；

Target function value computing module, for utilizing target function value computing formula, calculates each individual target function value；Its In, the target function value computing formula is：

$P_i=\underset{j=1}{\overset{D}{\Σ}}{P}_{j,i}-C\underset{j=1}{\overset{D}{\Σ}}{P}_{PENj,i};$

In formula, P_iI-th target function value of individuality is represented, C represents default penalty factor；

Evolution of Population processing module, for based on differential evolution algorithm, and combines each individual target function value, to it is described just Beginning population carries out Evolution of Population treatment, untill convergence, obtains global optimum individual；Wherein, the global optimum is individual Vector expression is：

x_g=[x_1g,x_2g,…x_Dg,x_busg]；

In formula, x_jgRepresent the optimal propeller pitch angle of jth typhoon group of motors in the wind-powered electricity generation group of planes, x_busgRepresent frequency-changing AC bus Optimal frequency；

Wind turbine group control module, for individual according to the global optimum, corresponding control is launched to the wind-powered electricity generation group of planes, with The propeller pitch angle of the every typhoon group of motors of control and the frequency of frequency-changing AC bus.

8. wind-powered electricity generation group of planes power control system according to claim 7, it is characterised in that

The initialization of population module, specifically for utilizing initialization of population formula, is pre-created and initializes and the wind-powered electricity generation The adaptable population of the scale of a group of planes, obtains the initial population；Wherein, the initialization of population formula is：

$\left\{\begin{array}{c}{x}_{j,i,0}=rand\[0,1\]({x_j}^U-{x_j}^L)+{x_j}^L\\ x_{bus,i,0}=rand\[0,1\]({x_{bus}}^U-{x_{bus}}^L)+{x_{bus}}^L\end{array}\right.;$

In formula, rand [0,1] represents the uniform random number for producing between zero and one, x_j ^URepresent the propeller pitch angle regulated value of Wind turbines Higher limit, x_j ^LRepresent the lower limit of the propeller pitch angle regulated value of Wind turbines, x_bus ^URepresent the frequency regulation of frequency-changing AC bus The higher limit of value, x_bus ^LRepresent the lower limit of the frequency adjustment value of frequency-changing AC bus.

9. wind-powered electricity generation group of planes power control system according to claim 7, it is characterised in that the penalty function determining module, Including：

Power output calculating sub module, for calculating formula using power output, calculates the power output per typhoon group of motors；Its In, the power output calculates formula and is：

$P_{j,i}=\frac{1}{2}{\mathrm{\ρAC}}_j(V_j,x_{j,i},x_{bus,i}){V_j}^3;$

In formula, P_j,iThe power output of the jth typhoon group of motors in i-th individuality is represented, ρ represents atmospheric density, and A represents blade Swept area, C_j(V_j,x_j,i,x_bus,i) represent jth typhoon group of motors power coefficient, V_jRepresent the real-time of jth typhoon group of motors Wind speed, x_j,iRepresent the propeller pitch angle of the jth typhoon group of motors in i-th individuality, x_bus,iFrequency-changing AC is female in representing i-th individuality The frequency of line；

Penalty function determination sub-module, for determining formula using the power output and the penalty function per typhoon group of motors, really Make the penalty function of every typhoon group of motors.

10. wind-powered electricity generation group of planes power control system according to claim 9, it is characterised in that the output power meter operator Module, including：

Power coefficient computing unit, for utilizing power coefficient computing formula, calculates the power coefficient per typhoon group of motors；Its In, the power coefficient computing formula is：

$C_j\left(V_j,x_{j,i},x_{bus,i}\right)=K_1\left(\frac{K_2}{\frac{2{\mathrm{\πx}}_{bus,i}R}{n_p{n}_g{V}_j}+0.08x_{j,i}}-\frac{0.035K_2}{{x_{j,i}}^3+1}-K_3{x}_{j,i}-K_4\right)e^S+K_5\frac{2{\mathrm{\πx}}_{bus,i}R}{n_p{n}_g{V}_j};$

Wherein, index S is：

$S=\frac{-K_6}{\frac{2{\mathrm{\πx}}_{bus,i}R}{n_p{n}_g{V}_j}+0.08x_{j,i}}+\frac{0.035K_6}{{x_{j,i}}^3+1};$

In formula, R represents the blade radius of Wind turbines, n_pRepresent motor number of pole-pairs, n_gRepresent gear-box no-load voltage ratio, K₁To K₆For advance Based on the coefficient that Wind turbines aerofoil profile determines.