Embodiment one
The embodiment of the present invention provides a kind of electrical network automatic electricity generation control system controller parameter optimization method, as it is shown in figure 1,
This electrical network automatic electricity generation control system controller parameter optimization method includes:
Step S1, set up electrical network automatic electricity generation control system phantom, electrical network automatic electricity generation control system phantom
Including controller.
Exemplarily, the electrical network automatic electricity generation control system phantom containing fuzzy controller can be set up, specifically set up
The method of the electrical network automatic electricity generation control system phantom containing fuzzy controller and existing foundation are containing fuzzy-adaptation PID control
The method of the electrical network automatic electricity generation control system phantom of device is identical, and those skilled in the art can refer to existing method and set up
Electrical network automatic electricity generation control system phantom, the embodiment of the present invention is not defined.
Step S2, according to electrical network automatic electricity generation control system phantom, set up the optimization of controller parameter and adjust model.
Exemplarily, the conventional integral square error index of control system, time can be used to take advantage of integral square error index
Or the time take advantage of Error Absolute Value integration index etc. to adjust as the optimization of controller noted above parameter the object function of model, this area
Technical staff can be arranged according to practical situation, and the embodiment of the present invention is not defined.
Step S3, model of adjusting according to optimization, by microhabitat bacterial foraging algorithm, it is thus achieved that controller parameter the most excellent
Change result.
Step S4, according to initial optimization result, pass through pattern search algorithm, it is thus achieved that the final optimization pass of controller parameter knot
Really.
In the technical scheme of the present embodiment, according to electrical network automatic electricity generation control system phantom, establish electrical network certainly
The optimization of dynamic power-generating control system simulation model controller parameter is adjusted after model, first passes through microhabitat bacterial foraging algorithm,
Obtain the initial optimization result of this controller parameter, after obtaining initial optimization result, continue to use pattern search algorithm,
The initial optimization result of this controller parameter is optimized further, thus obtains the controller parameter that more optimizes
Whole optimum results, by being combined with pattern search algorithm by microhabitat bacterial foraging algorithm, substantially increases controller parameter
The optimization precision of optimum results, thus improve the dynamic control performance of electrical network automatic electricity generation control system, it is ensured that electric power
The safe operation of system.
Exemplarily, electrical network automatic electricity generation control system phantom as shown in Figure 2, this electrical network automatic generation can be set up
Control System Imitation model includes all including in first area and second area, first area and second area fuzzy-adaptation PID control
Device 1, speed regulator 2, steam turbine 3 etc..
In figure, Δ f1For the frequency departure of first area, Δ f2Frequency departure for second area;ΔPt12It it is first area
Flow to the dominant eigenvalues deviation of second area;B1For the frequency bias coefficient of first area, B2Frequency departure for second area
Coefficient;ACE1For the district control deviation of first area, ACE2District control deviation for second area;u1For first area
The output of fuzzy controller, u2Output for the fuzzy controller of second area;R1Poor for the tune of unit in first area
Coefficient, R2For the difference coefficient of unit in second area;TG1For speed regulator time constant in first area, TG2For in second area
Speed regulator time constant;TT1For steam turbine time constant in first area, TT2For steam turbine time constant in second area;Δ
PG1For unit output variable quantity in first area, Δ PG2For unit output variable quantity in second area;ΔPT1For in first area
Steam turbine valve variable quantity, Δ PT2For steam turbine valve variable quantity in second area;ΔPL1For load variations in first area
Amount, Δ PL2For load variations amount in second area;KPS1Gain merit frequency transform coefficient for system in first area, KPS2It it is the secondth district
In territory, system is gained merit frequency transform coefficient;TPS1For first area system time constant, TPS2For second area system time constant;
T12Interconnection synchronization factor for first area Yu second area;α12Dominant eigenvalues for first area with second area turns
Change coefficient, when two zonal basis power are identical, α12=-1.
Specifically, the structure of above-mentioned fuzzy controller 1 is as it is shown on figure 3, this fuzzy controller includes K1、K2、K3、
K4Four scale factors (i.e. the controller parameter of this fuzzy controller), differentiator 11 and integrator 12, wherein, fuzzy
The district control deviation (ACE) that input is above-mentioned control system of controller;x1And x2The input being respectively fuzzy logic control becomes
Amount, y1And y2It is respectively the output variable of fuzzy logic control;The regulation that output u is above-mentioned control system of fuzzy controller
Power.
Exemplarily, the time can be used to take advantage of Error Absolute Value integration index to imitate as above-mentioned electrical network automatic electricity generation control system
The optimization of true mode middle controller parameter is adjusted the object function of model.Specifically, the optimization of the parameter of this controller is adjusted mould
Type includes object function and inequality constraints condition:
Object function is:
Wherein, T is the simulation time of electrical network automatic electricity generation control system phantom, Δ f1Frequency for first area is inclined
Difference, Δ f2For the frequency departure of second area, Δ Pt12The dominant eigenvalues deviation of second area is flowed to for first area.
Inequality constraints condition includes:
Wherein, K1、K2、K3And K4It is controller parameter.It should be noted that " min " in the embodiment of the present invention refers both to
Be minima, " max " refers to maximum.
Exemplarily, step S3 is adjusted model according to optimization, by microhabitat bacterial foraging algorithm, it is thus achieved that controller is joined
The concrete steps of the initial optimization result of number include:
Step S31, obtain the initial value of many group controllers parameter.
Step S32, adjust according to optimization model and the initial value of every group controller parameter, it is thus achieved that optimize mesh in model of adjusting
The initial value of scalar functions.
Exemplarily, above-mentioned electrical network automatic electricity generation control system emulation mould can be built in Matlab/Simulink software
Type, for each group of (K1,K2,K3,K4) initial value, be assigned in phantom corresponding controller module, passed through
Matlab/Simulink simulation calculation goes out corresponding Δ f1、Δf2、ΔPt12Value, and then calculate the initial of object function J
Value.
Step S33, by microhabitat bacterial foraging algorithm, the initial value of combined objective function, to often organizing scale factor
Initial value is iterated calculating, it is thus achieved that initial optimization result.
Exemplarily, the basic step of above-mentioned microhabitat bacterial foraging algorithm as shown in Figure 4, look for food calculation by this microhabitat antibacterial
The basic step of method includes:
Step 1: initialize the relevant parameter of microhabitat bacterial foraging algorithm.
Specifically, the relevant parameter of above-mentioned microhabitat bacterial foraging algorithm includes: dimension p to be optimized, and the present invention is real
Execute in example, have four scale factors to be optimized, now, p=4, antibacterial number S (be taken scale factor group
Number), chemotactic times Nc, maximum travelling times Ns, breed times Nre, migrate times Ned, transition probability Ped, attracting factor daWith
ωa, rejection factor hrAnd ωr, microhabitat radius δ0With spacing distance L etc..
You need to add is that, in the embodiment of the present invention, initial flora position can be randomly generated.
Step 2: judge whether bacterial migration number of times reaches maximum and migrate number of times, the most then export optimum results;If it is not,
Then enter next step.
Step 3: judge whether flora breeding number of times reaches maximum breeding number of times, the most then forward step 8 to, enter and migrate
Operation;If it is not, then enter next step.
Step 4: judge whether antibacterial maximum chemotactic number of times reaches maximum chemotactic number of times, the most then forward step 6 to, enters
Breeding operation;If it is not, then enter next step.
Step 5: each organisms in flora is carried out chemotactic operation, then goes to step 4.
Step 6: the flora obtained after chemotactic is passed through health degree shared mechanism and its health degree of limit competition Developing Tactics.
Step 7: use roulette method choice a new generation bacterial population, then go to step 3.
Step 8: antibacterial is performed migration operation, then goes to step 2.
By above-mentioned microhabitat bacterial foraging algorithm, can increase the multiformity of population, the overall situation that improve this algorithm is sought
Excellent ability, it is to avoid problem that local optimum occurs.
Additionally, after obtaining above-mentioned initial optimization result, according to initial optimization result in above-mentioned steps S4, pass through pattern
Searching algorithm, it is thus achieved that the concrete steps of the final optimization pass result of controller parameter include:
Step S41, according to initial optimization result, it is thus achieved that the initial optimization value of a group controller parameter.
Step S42, initial optimization value according to controller parameter, it is thus achieved that optimize the most excellent of object function in model of adjusting
Change value.
Exemplarily, the concrete grammar of the initial optimization value of above-mentioned acquisition object function and step S32 obtain target letter
The concrete grammar of the initial value of number is similar to, and those skilled in the art can refer to obtain in step S32 the initial value of object function
Method obtains the initial optimization value of object function, the most no longer repeats.
Step S43, by pattern search algorithm, the initial optimization value of combined objective function, to controller noted above parameter
Initial optimization value is iterated calculating, it is thus achieved that final optimization pass result.
Exemplarily, above-mentioned pattern search algorithm basic step as it is shown in figure 5, in figure A represent by above-mentioned microhabitat
The initial optimization result that bacterial foraging algorithm obtains, " PS " is writing a Chinese character in simplified form of " pattern search " pattern search, " poll " table
Show structure dragnet Mesh centered by starting point, find the process of the point making object function decline.Specifically, this pattern search
The basic step of algorithm includes:
Step 1:
Given initial point x0, and calculate the target function value f (x of initial point0), given initialization step-size in search Δ0, expand
Factor lambda and contraction factor β, and set maximum iteration time gmax, make g=0.
Step 2:
Current search starting point x is setk。
Step 3:
Structural model vector ViAnd site Mk,i, i=1,2 ..., n.
Step 4:
Calculate the target function value f (M of each sitek,i)。
Step 5:
If searching any one some Mk,i, have f (Mk,i) < f (xk), then poll success, makes xk+1=Mk,i, and press formula
Δk+1=λ ΔkUpdate Δk+1, wherein, λ > 1;
If searching any one some Mk,i, have f (Mk,i) > f (xk), then poll is unsuccessful, makes xk+1=xk, and press formula
Δk+1=β ΔkUpdate Δk+1, wherein, 0 < β < 1.
Step 6:
Iterations g=g+1, if g is < gmax, forward step 2 to, enter next iteration;Otherwise, algorithm terminates, and exports excellent
Change result.
The present invention using the initial optimization result that obtained by microhabitat bacterial foraging algorithm as at the beginning of pattern search algorithm
Initial point, by performing fine search at microhabitat bacterial foraging algorithm, it is possible to carry near the initial optimization result that obtains further
The precision of the optimum results of high gained.
For the ease of skilled artisan understands that and implementing, below the embodiment of the present invention be given one use above-mentioned control
The concrete application example of device parameter optimization method.
Specifically, in Matlab/Simulink software, build electrical network automatic electricity generation control system emulation as shown in Figure 2
Model, for simplifying the analysis, it is assumed that two district system parameters are identical, and the parameter of fuzzy controller 1 is the most identical.Set:
K1、K2、K3、K4Span be [0,2].
The value of the initial parameter of microhabitat bacterial foraging algorithm is as follows:
P=4, S=20, Nc=8, Ns=3, Nre=4, Ned=2, Ped=0.25, da=0.01, ωa=0.04, hr=
0.01, ωr=10, δ0=1, L=0.1.
The value of the initial parameter of pattern search algorithm is as follows:
λ=2, β=0.5, Δ0=1, gmax=50.
The value of the initial parameter of grid control system is as follows:
TG1=TG2=0.08s, Tt1=Tt2=0.3s, TPS1=TPS2=20s, KPS1=KPS2=120Hz/p.u., T12=
0.545p.u./Hz, B1=B2=0.425p.u./Hz, R1=R2=2.4Hz/p.u., α12=-1.
Assume that first area load increases 0.1p.u., be respectively adopted pattern search algorithm, microhabitat bacterial foraging algorithm and
Algorithm microhabitat bacterial foraging algorithm and pattern search algorithm combined provided in the embodiment of the present invention, to fuzzy
The scale factor of controller is optimized, and maximum iteration time is all set to 50 times.Concrete optimization process is as follows:
Step 1: randomly generate 20 initial population in the range of search volume [0,2], the individual generation of each in population
One group of (K of table1,K2,K3,K4) value.
Step 2: the value of each individuality is updated in built electrical network automatic electricity generation control system phantom,
Each individual corresponding target function value is gone out by Matlab/Simulink simulation calculation.Exemplarily, simulation time can be set
It is set to 15 seconds.
Step 3: be respectively adopted in pattern search algorithm, microhabitat bacterial foraging algorithm and the embodiment of the present invention provide general
Population at individual is iterated calculating by the algorithm that microhabitat bacterial foraging algorithm and pattern search algorithm combine.
Step 4: judge whether to reach the condition of convergence.The most then export optimal result;If it is not, then forward under step 2 carries out
An iteration optimizing.
Step 5: calculate the electrical network automatic electricity generation control system response under optimal result.
As shown in Figure 6, " PS " in Fig. 6 represents pattern search algorithm to each convergence of algorithm curve, and " NBFO " represents your pupil
Border bacterial foraging algorithm, " hNBFO-PS " represent in the embodiment of the present invention provide microhabitat bacterial foraging algorithm and pattern are searched
The algorithm that rope algorithm combines, " ITAE " is to be multiplied by Error Absolute Value integration the time, and the optimization that i.e. embodiment of the present invention provides is whole
The object function of cover half type.It will be appreciated from fig. 6 that use microhabitat bacterial foraging algorithm and pattern search algorithm to be iterated excellent respectively
During change, after iteration 30 times, the most substantially enter convergence state.Therefore, in use, microhabitat bacterial foraging algorithm and pattern are searched
When the algorithm that rope algorithm combines is iterated optimizing, it is arranged on use microhabitat bacterial foraging algorithm and is iterated 30 and takes second place
After, proceed to use pattern search algorithm to proceed iteration optimization.
Use pattern search algorithm, microhabitat bacterial foraging algorithm and the embodiment of the present invention provide by microhabitat antibacterial
The algorithm that foraging algorithm and pattern search algorithm combine is iterated calculating optimum results such as table 1 institute of gained to population at individual
Show.
Table 1 is by the best proportion factor of the fuzzy controller of algorithms of different gained
|
K1 |
K2 |
K3 |
K4 |
Pattern search algorithm |
0.7317 |
0.6477 |
0.4509 |
0.5470 |
Microhabitat bacterial foraging algorithm |
0.9393 |
0.4998 |
0.7191 |
0.8834 |
Microhabitat bacterial foraging algorithm and pattern search algorithm combination algorithm |
1.5393 |
0.4998 |
0.8691 |
1.0834 |
The best proportion factor of the fuzzy controller of algorithms of different gained is applied to above-mentioned electrical network Automatic Generation Control
After in system, dynamic time-domain response curve such as Fig. 7 institute of the frequency departure of first area in this electrical network automatic electricity generation control system
Show, the dynamic time-domain response curve of the frequency departure of second area as shown in Figure 8, the contact between first area and second area
As it is shown in figure 9, wherein, " PS " in Fig. 7, Fig. 8, Fig. 9 all represents pattern search to the dynamic time-domain response curve of linear heat generation rate deviation
Algorithm, " NBFO " all represents microhabitat bacterial foraging algorithm, " hNBFO-PS " all represent in the embodiment of the present invention provide by little
The algorithm that habitat bacterial foraging algorithm and pattern search algorithm combine.
As shown in Figure 7, Figure 8 and Figure 9, by what the optimization method provided in the use embodiment of the present invention obtained, there is ratio of greater inequality
After the controller of the example factor is applied to electrical network automatic electricity generation control system, it is possible to make this electrical network automatic electricity generation control system every
Dominant eigenvalues deviation between the frequency departure in individual region and two regions returns to zero soon, substantially increases electrical network
The dynamic time domain response performance of automatic electricity generation control system, it is ensured that the safe operation of power system.
The above, the only detailed description of the invention of the present invention, but protection scope of the present invention is not limited thereto, and any
Those familiar with the art, in the technical scope that the invention discloses, can readily occur in change or replace, should contain
Cover within protection scope of the present invention.Therefore, protection scope of the present invention should be as the criterion with described scope of the claims.