Reactive power optimization system and method of power grid based on the double-fish-swarm algorithm
10439398 ยท 2019-10-08
Assignee
Inventors
- Huaguang Zhang (Shenyang, CN)
- Jun Yang (Shenyang, CN)
- Qiuye Sun (Shenyang, CN)
- Feiye Wu (Shenyang, CN)
- Xinrui Liu (Shenyang, CN)
- Dongsheng Yang (Shenyang, CN)
- Zhiliang Wang (Shenyang, CN)
- Jian Feng (Shenyang, CN)
- Bonan Huang (Shenyang, CN)
- Yanhong Luo (Shenyang, CN)
- Guotao Hui (Shenyang, CN)
- Zhenwei Liu (Shenyang, CN)
Cpc classification
Y04S40/20
GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
Y02E40/30
GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
G05F1/14
PHYSICS
Y02E40/70
GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
H02J2203/20
ELECTRICITY
Y04S10/50
GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
H02J3/1885
ELECTRICITY
H02J3/18
ELECTRICITY
Y02E60/00
GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
International classification
H02J3/18
ELECTRICITY
H02J3/00
ELECTRICITY
G05F1/14
PHYSICS
Abstract
The invention provides a reactive power optimization system and a method of a power grid based on a double-fish-swarm algorithm. The system includes a power grid state data acquiring module, a reactive power regulating module and a reactive power executing module. The power grid state data acquiring module includes a power grid state data acquiring processor and a relay transmitter. The reactive power regulating module is a control terminal. The reactive power executing module includes generator terminal voltage regulators, transformer tap regulators and reactive power compensation regulators. The method is used for acquiring the initial data to be optimized in the current network; and optimizing the initial data to be optimized in the current network based on a double-fish-swarm algorithm so as to obtain optimal value of control variables in the power grid. According to the method, the distribution network to be optimized can realize reasonable reactive power flow distribution.
Claims
1. A reactive power optimization system of a power grid based on a double-fish-swarm algorithm, comprising a power grid state data acquiring module, a reactive power regulating module and a reactive power executing module, the power grid state data acquiring module connected to the reactive power regulating module, the reactive power regulating module connected to the reactive power executing module, wherein the power grid state data acquiring module comprises a power grid state data acquiring processor and a relay transmitter; the reactive power regulating module is a control terminal; the reactive power executing module comprises generator terminal voltage regulators, transformer tap regulators and reactive power compensation regulators; an input end of the power grid state data acquiring processor is connected with a power grid; an output end of the power grid state data acquiring processor is connected with the input end of the relay transmitter; the output end of the relay transmitter is connected with the input end of the control terminal; the output end of the control terminal is connected with the input end of each of the generator terminal voltage regulators, the input end of each of the transformer tap regulators and the input end of each of the reactive power compensation regulators; the output end of each of the generator terminal voltage regulators is connected with a corresponding generator in the power grid; the output end of each of the transformer tap regulators is connected with a corresponding transformer in the power grid; and the output end of each of the reactive power compensation regulators is connected with a corresponding reactive power compensation device in the power grid; the power grid state data acquiring processor is used for acquiring the current network information of the power grid and judging whether the current network information meets the optimal state required by the power grid or not; and if the current network information cannot meet the optimal state required by the power grid, the current network information is transmitted to the relay transmitter, wherein the current network information includes power grid node information, branch information, generator information, transformer information and reactive power compensation device information; the relay transmitter is used for transmitting the power grid node information, the branch information, the generator information, the transformer information and the reactive power compensation device information which are acquired by the power grid state data acquiring processor and required for reactive power optimization in the network information to the control terminal; the reactive power regulating module comprises a processor configured to function as a parameter acquiring unit, a double-fish-swarm algorithm-base reactive power optimization unit, and an optimization decision-making control unit, the parameter acquiring unit connected to the reactive power optimization unit, and the reactive power optimization unit connected to the optimization decision-making control unit; the parameter acquiring unit is used for acquiring the power grid node information, the branch information, the generator information, the transformer information and the reactive power compensation device information transmitted by the relay transmitter and used as initial data to be optimized in the current network; the double-fish-swarm algorithm-base reactive power optimization unit is used for establishing a mathematical model for reactive power optimization of a power system, and the acquired initial data to be optimized in the current network is optimized based on the double-fish-swarm algorithm, so that an optimal value of each of control variables in the power grid is obtained, wherein the control variables comprise generator terminal voltage amplitudes, transformer adjustable ratios, and the reactive power capacity of reactive power compensation devices; the optimization decision-making control unit is used for transmitting optimal values of the control variables to the reactive power executing module; the generator terminal voltage regulators are used for regulating generator terminal voltages according to the optimal values of the generator terminal voltage amplitudes, obtained by the reactive power regulating module; the transformer tap regulators are used for regulating transformer taps according to the optimal values of the transformer adjustable ratios, obtained by the reactive power regulating module; and the reactive power compensation regulators are used for regulating the compensation capacity of reactive power compensators according to the optimal values of the reactive power capacity of the reactive power compensation devices, obtained by the reactive power regulating module.
2. The reactive power optimization system of a power grid based on the double-fish-swarm algorithm according to claim 1, wherein the mathematical model for establishing the reactive power optimization of the power system is shown as follows:
3. The reactive power optimization system of a power grid based on the double-fish-swarm algorithm according to claim 2, wherein a specific process for optimizing the acquired initial data to be optimized in the current network based on the double-fish-swarm algorithm so as to obtain the optimal values of the control variables in the power grid comprises the steps of: generating a small fish-swarm and a ferocious fish-swarm in a random manner according to the upper limit value and the lower limit value of the control variables to be optimized in the power grid; initializing the optimal value of current food concentration, and displaying the initialized optimal value of current food concentration in a bulletin board: by using the reciprocal 1/F of the objective function of the mathematical model for reactive power optimization of a power grid as a value FC of food concentration, calculating the value of food concentration in the current network and of small fishes and ferocious fishes under an initial data condition, taking the maximum value of food concentration as the optimal value of current food concentration, displaying the optimal value of current food concentration in the bulletin board, and saving the state and the current value FC; enabling the small fish-swarm to act: enabling small fish individuals to perform a foraging clustering behavior, a rear-end behavior and a protective clustering behavior according to the distance between each of the small fish individuals and ferocious fish individuals as well as the distance between the small fish individuals in the small fish-swarm, and updating the value in the bulletin board Y1; enabling the ferocious fish-swarm to act: enabling each of the ferocious fish individuals to perform a predation behavior, a tracking behavior and a clustering behavior according to the distance between each of the ferocious fishes and the small fish individuals as well as the distance between the ferocious fish individuals in the ferocious fish-swarm, and updating the value in the bulletin board Y2; and when iteration times achieve the maximum iteration times, taking the greater values of food concentration displayed in the bulletin boards Y1 and Y2 as well as the state thereof as optimal results so as to obtain the optimal values of the control variables, namely the optimal values of the generator terminal voltage amplitudes, the transformer adjustable ratios, and the reactive power capacity of the reactive power compensation devices.
4. A method for reactive power optimization of a power grid by adopting the reactive power optimization system of a power grid based on the double-fish-swarm algorithm according to claim 3, comprising the steps of: step 1: acquiring the current network information of the power grid state data and judging whether the current network information meets the optimal state required by the power grid or not by the power grid state data acquisitor; and if the current network information cannot meet the optimal state required by the power grid, transmitting the current network information to the relay transmitter; step 2: transmitting the power grid node information, the branch information, the generator information, the transformer information and the reactive power compensation device information which are acquired by the power grid state data acquiring processor and required for reactive power optimization in the network information to the control terminal, by the relay transmitter; step 3: acquiring the power grid node information, the branch information, the generator information, the transformer information and the reactive power compensation device information transmitted by the relay transmitter, by the control terminal as initial data to be optimized in the current network; step 4: establishing the mathematical model for the reactive power optimization of the power system through the control terminal, optimizing the acquired initial data to be optimized in the current network based on the double-fish-swarm algorithm so as to obtain the optimal values of the control variables in the power grid, and transmitting the optimal values of the control variables to the reactive power executing module; and step 5: regulating generator terminal voltages by using the generator terminal voltage regulators according to the optimal values of the generator terminal voltage amplitudes, regulating the transformer taps by using the transformer tap regulators according to optimal values of the transformer adjustable ratios, and regulating the compensation capacity of the reactive power compensators by using the reactive power compensation regulators according to the optimal values of the reactive power capacity of the reactive power compensation devices.
5. The method for reactive power optimization of a power grid based on the double-fish-swarm algorithm according to claim 4, wherein the step 4 comprises the steps of: step 4.1: establishing a mathematical model for the reactive power optimization of the power system by using an alpha-method, with the minimum network loss and the optimal voltage level as the optimization objectives:
6. The method for reactive power optimization of a power grid based on the double-fish-swarm algorithm according to claim 5, wherein the step 4.3 comprises the steps of: step 4.3.1: initializing the number t.sub.1=0 of the generated small fishes; step 4.3.2: randomly generating the small fish individuals according to the value ranges of the control variables to be optimized in the power grid:
X=[U.sub.G11,U.sub.G12, . . . U.sub.G1g,K.sub.T11,K.sub.T12, . . . ,K.sub.T1t,Q.sub.C11,Q.sub.C12, . . . Q.sub.C1c]=[x.sub.1,x.sub.2, . . . ,x.sub.s, . . . ,x.sub.S]; wherein, X is the value sequence of the control variables of the small fish individuals, x.sub.s=x.sub.s min+rand(g)(x.sub.s maxx.sub.s min), rand(g) is a random number of an interval (0, 1), x.sub.s min and x.sub.s max are respectively the lower limit value and the upper limit value of the corresponding control variables, U.sub.G1g, K.sub.T1t and Q.sub.C1c are respectively the generator terminal voltage amplitude, the transformer adjustable ratio, and the reactive power capacity of the reactive power compensation devices in the randomly generated small and artificial fish individuals, x.sub.s is the control variables of the randomly generated small fish individuals, 1ss, S=N.sub.g+N+N.sub.c, N.sub.g is the number of the generators in the power grid, N.sub.i is the number of the transformers in the power grid, and N.sub.c is the number of the reactive power compensation devices in the power grid; step 4.3.3: performing power flow calculation on the control variables in the randomly generated small fish individuals in the step 4.3.2 by using a PQ decomposition method, if the power flow value of the current small fish individual is converged, retaining the small fish individual, generating a small fish number t.sub.1=t.sub.1+1 of the current small fish individual, and executing the step 4.3.4; or else, not retaining the small fish individual, and returning to the step 4.3.2; and step 4.3.4: if t.sub.1N1, obtaining a current small fish-swarm, and executing the step 4.4; or else, returning to the step 4.3.2.
7. The method for reactive power optimization of a power grid based on the double-fish-swarm algorithm according to claim 5, wherein the step 4.4 comprises the steps of: step 4.4.1: initializing the number t.sub.2=0 of the generated ferocious fishes; step 4.4.2: randomly generating the ferocious fish individuals according to the value ranges of the control variables to be optimized in the power grid:
W=[U.sub.G21,U.sub.G22, . . . U.sub.G2g,K.sub.T21,K.sub.T22, . . . ,K.sub.T2t,Q.sub.C21,Q.sub.C22, . . . ,Q.sub.C2c]=[w.sub.1,w.sub.2, . . . ,w.sub.s, . . . ,w.sub.S]; wherein, W is the value sequence of the control variables of the ferocious fish individuals, w.sub.s=x.sub.s max+rand(g)(x.sub.s maxx.sub.s min), U.sub.G2g, K.sub.T2t and Q.sub.C2c are respectively the generator terminal voltage amplitude, the transformer adjustable ratio, and the reactive power capacity of the reactive power compensation devices in the randomly generated ferocious fish individuals, and w.sub.s is the control variables of the randomly generated ferocious fish individuals; step 4.4.3: performing power flow calculation on the control variables in the randomly generated ferocious fish individuals in the step 4.4.2 by using a PQ decomposition method, if the power flow value of the current ferocious fish individual is converged, retaining the ferocious fish individual, generating a ferocious fish number t.sub.2=t.sub.2+1 of the current ferocious fish individual, and executing the step 4.4.4; or else, not retaining the ferocious fish individual, and returning to the step 4.4.2; and step 4.4.4: if t.sub.2N2, obtaining a current ferocious fish-swarm, and executing the step 4.5; or else, returning to the step 4.4.2.
8. The method for reactive power optimization of a power grid based on the double-fish-swarm algorithm according to claim 5, wherein the step 4.6 comprises the steps of: step 4.6.1: determining the distance D={D.sub.i1, D.sub.i2, . . . , D.sub.ij, . . . , D.sub.iN2} between the current i small fish individual X.sub.i and each of the ferocious fish individuals W.sub.j, wherein, D.sub.ij=X.sub.iW.sub.j, 1<j<N2, and N2 is the number of the ferocious fish individuals; step 4.6.2: if the distance D.sub.ijVisual1 between the current i th small fish individual X.sub.i and each of the ferocious fish individuals W.sub.j, executing the step 4.6.3; or else, executing the step 4.6.5; step 4.6.3: determining the total number n of the small fish individuals within the perception range of the i th small fish individual X.sub.i, to determine a protective clustering safe position X.sub.safe, and enabling the i th small fish individual X.sub.i to perform the protective clustering behavior so as to obtain an update value X.sub.inext1 of the i th small fish individual performing the protective clustering behavior, wherein the calculation formula for the protective clustering safe position X.sub.safe is as follows:
X.sub.safe=(X.sub.i1+X.sub.i2+ . . . +X.sub.in)/n, wherein, =(n+p)/n is an escape factor, p is the total number of the ferocious fish individuals within the perception range Visual1 of the small fish individual X.sub.i, and X.sub.in is the small fish individuals within the perception range of the small fish individual X.sub.i; the calculation formula for the update value X.sub.inext1 of the i th small fish individual performing the protective clustering behavior is as follows:
X.sub.inext1=X.sub.i+rand(g)Step1.sub.1(X.sub.safeX.sub.i)/X.sub.safeX.sub.i; step 4.6.4: judging whether the update value X.sub.inext1 of the current small fish individual performing the protective clustering behavior meets a constraint condition and the power flow is converged; if yes, recording the state of the update value X.sub.inext1, calculating the value of food concentration, updating the value in the bulletin board Y1, and executing the step 4.7; or else, returning to the step 4.6.3; step 4.6.5: determining the small fish individual X.sub.j with the maximum value FC.sub.j of food concentration in the small fish-swarm within the perception range of the i th small fish individual X.sub.i, and enabling the i th small fish individual X.sub.i to perform the rear-end behavior so as to obtain an update value X.sub.inext2 of the i th small fish individual X.sub.i performing the rear-end behavior, wherein the calculation formula for the update value X.sub.inext2 of the i th small fish individual X.sub.i performing the rear-end behavior is as follows:
X.sub.inext2=X.sub.i+rand(g)Step1.sub.1(X.sub.jX.sub.i)/X.sub.jX.sub.i; step 4.6.6: determining the total number n of the small fish individuals within the perception range of the i th small fish individual X.sub.i to determine the central position X.sub.c1 of the small fish individuals performing foraging clustering behavior, and enabling the i th small fish individual X.sub.i to perform the foraging clustering behavior so as to obtain an update value X.sub.inext3 of the i th small fish individual performing the foraging clustering behavior, wherein the calculation formula for the central position X.sub.c1 of the small fish individuals performing foraging clustering behavior is as follows:
X.sub.c1=(X.sub.i1+X.sub.i2+ . . . +X.sub.in)/n; the calculation formula for the update value X.sub.inext3 of the i th small fish individual performing the foraging clustering behavior is as follows:
X.sub.inext3=X.sub.i+rand(g)Step1.sub.1(X.sub.c1X.sub.i)/X.sub.c1X.sub.i; and step 4.6.7: judging whether the update value X.sub.inext2 of the current small fish individual performing the rear-end behavior and the update value X.sub.inext3 of the current small fish individual performing foraging clustering behavior meet constraint conditions and the power flow is converged; if yes, recording the state of the update value X.sub.inext2 and the state of the update value X.sub.inext3, calculating the value of food concentration, updating the value in the bulletin board Y1, and executing the step 4.7; or else, returning to the step 4.6.5.
9. The method for reactive power optimization of a power grid based on the double-fish-swarm algorithm according to claim 5, wherein the step 4.7 comprises the steps of: step 4.7.1: determining the distance L={L.sub.i1, L.sub.i2, . . . , L.sub.ij, . . . , L.sub.iN1} between the current ith ferocious fish individual W.sub.i and each of the small fish individuals X.sub.j, wherein L.sub.ij=W.sub.iX.sub.j, 1<j<N1, and N1 is the number of the small fish individuals; step 4.7.2: if the distance L.sub.ijVisual2 between the current ith ferocious and artificial fish individual W.sub.i and each of the small and artificial fish individuals X.sub.j, executing the step 4.7.3; or else, executing the step 4.7.6; step 4.7.3: determining the total number m of the small fish individuals within the perception range of the current ith ferocious fish individual W.sub.i, enabling the ith ferocious fish individual W.sub.i to perform the predation behavior with a central position X.sub.c2 of small fish-swarm within the perception range of the ith ferocious fish individual W.sub.i as the maximum value of current food concentration, and updating the ith ferocious fish individual W.sub.i as the central position X.sub.c2 of the small fish-swarm so as to obtain an update value X.sub.1next1=X.sub.c2 of the ith ferocious fish individual W.sub.i performing the predation behavior; the calculation formula for the central position X.sub.c2 of small fish-swarms within the perception range of the ith ferocious fish individual W.sub.i is as follows:
X.sub.c2=(X.sub.i1+X.sub.i2+ . . . +X.sub.im)/m; wherein X.sub.im is the small fishes within the perception range of the ferocious fish individual W.sub.i; step 4.7.4: enabling the ith ferocious fish individual W.sub.i to perform the tracking behavior with the central position X.sub.c2 of small fish-swarms within the perception range of the current ith ferocious fish individual W.sub.i as the maximum value of current food concentration so as to obtain an update value W.sub.inext2 of the ith ferocious fish individual W.sub.i performing the tracking behavior; the calculation formula for the update value W.sub.inext2 of the ith ferocious fish individual W.sub.i performing the tracking behavior is as follows:
W.sub.inext2=W.sub.j+rand(g)Step2.sub.2(X.sub.c2W.sub.i)/X.sub.c2W.sub.i; step 4.7.5: judging whether the update value W.sub.inext1 of the current ferocious fish individual W.sub.i performing the predation behavior and the update value W.sub.inext2 of the ferocious fish individual W.sub.i preforming the tracking behavior meet the constraint conditions and the power flow is converged; if yes, updating the value in the bulletin board Y2, and executing the step 4.7.6; or else, returning to the step 4.7.3; step 4.7.6: determining the total number r of the ferocious fish individuals within the perception range of the ith ferocious fish individual W.sub.i to determine the central position W.sub.c of the ferocious fish individuals preforming the clustering behavior, and enabling the i ferocious fish individual W.sub.i to perform the clustering behavior so as to obtain an update value W.sub.inext3 of the ith ferocious fish individual W.sub.i performing the clustering behavior, wherein the calculation formula for the central position W.sub.c of the fish individual performing the clustering behavior is as follows:
W.sub.c=(W.sub.i1+W.sub.i2+ . . . +W.sub.ir)/r; wherein, W.sub.ir is the ferocious fish individuals within the perception range of the ith ferocious fish individual W.sub.i; the calculation formula for the update value W.sub.inext3 of the ith ferocious fish individual W.sub.i performing the clustering behavior is as follows:
W.sub.inext1=W.sub.i+rand(g)Step2.sub.2(W.sub.cW.sub.i)/W.sub.cW.sub.i; and step 4.7.7: judging whether the update value W.sub.inext3 of the current ferocious fish individual W.sub.i performing the clustering behavior meets the constraint condition and the power flow is converged; if yes, recording the state of the update value W.sub.inext3, calculating the value of the food concentration, updating the value in the bulletin board Y2, and executing the step 4.8; or else, returning to the step 4.7.6.
10. A method for reactive power optimization of a power grid by adopting the reactive power optimization system of a power grid based on the double-fish-swarm algorithm according to claim 2, comprising the steps of: step 1: acquiring the current network information of the power grid state data and judging whether the current network information meets the optimal state required by the power grid or not by the power grid state data acquiring processor; and if the current network information cannot meet the optimal state required by the power grid, transmitting the current network information to the relay transmitter; step 2: transmitting the power grid node information, the branch information, the generator information, the transformer information and the reactive power compensation device information which are acquired by the power grid state data acquiring processor and required for reactive power optimization in the network information to the control terminal, by the relay transmitter; step 3: acquiring the power grid node information, the branch information, the generator information, the transformer information and the reactive power compensation device information transmitted by the relay transmitter, by the control terminal as initial data to be optimized in the current network; step 4: establishing the mathematical model for the reactive power optimization of the power system through the control terminal, optimizing the acquired initial data to be optimized in the current network based on the double-fish-swarm algorithm so as to obtain the optimal values of the control variables in the power grid, and transmitting the optimal values of the control variables to the reactive power executing module; and step 5: regulating generator terminal voltages by using the generator terminal voltage regulators according to the optimal values of the generator terminal voltage amplitudes, regulating the transformer taps by using the transformer tap regulators according to optimal values of the transformer adjustable ratios, and regulating the compensation capacity of the reactive power compensators by using the reactive power compensation regulators according to the optimal values of the reactive power capacity of the reactive power compensation devices.
11. The method for reactive power optimization of a power grid based on the double-fish-swarm algorithm according to claim 10, wherein the step 4 comprises the steps of: step 4.1: establishing a mathematical model for the reactive power optimization of the power system by using an alpha-method, with the minimum network loss and the optimal voltage level as the optimization objectives:
12. The method for reactive power optimization of a power grid based on the double-fish-swarm algorithm according to claim 11, wherein the step 4.3 comprises the steps of: step 4.3.1: initializing the number t.sub.1=0 of the generated small fishes; step 4.3.2: randomly generating the small fish individuals according to the value ranges of the control variables to be optimized in the power grid:
X=[U.sub.G11,U.sub.G12, . . . U.sub.G1g,K.sub.T11,K.sub.T12, . . . ,K.sub.T1t,Q.sub.C11,Q.sub.C12, . . . Q.sub.C1c]=[x.sub.1,x.sub.2, . . . ,x.sub.s, . . . ,x.sub.S]; wherein, X is the value sequence of the control variables of the small fish individuals, x.sub.s=x.sub.s min+rand(g)(x.sub.s maxx.sub.s min), rand(g) is a random number of an interval (0, 1), x.sub.s min and x.sub.s max are respectively the lower limit value and the upper limit value of the corresponding control variables, U.sub.G1g, K.sub.T1t and Q.sub.C1c are respectively the generator terminal voltage amplitude, the transformer adjustable ratio, and the reactive power capacity of the reactive power compensation devices in the randomly generated small and artificial fish individuals, x.sub.s is the control variables of the randomly generated small fish individuals, 1ss, S=N.sub.g+N+N.sub.c, N.sub.g is the number of the generators in the power grid, N.sub.i is the number of the transformers in the power grid, and N.sub.c is the number of the reactive power compensation devices in the power grid; step 4.3.3: performing power flow calculation on the control variables in the randomly generated small fish individuals in the step 4.3.2 by using a PQ decomposition method, if the power flow value of the current small fish individual is converged, retaining the small fish individual, generating a small fish number t.sub.1=t.sub.1+1 of the current small fish individual, and executing the step 4.3.4; or else, not retaining the small fish individual, and returning to the step 4.3.2; and step 4.3.4: if t.sub.1N1, obtaining a current small fish-swarm, and executing the step 4.4; or else, returning to the step 4.3.2.
13. The method for reactive power optimization of a power grid based on the double-fish-swarm algorithm according to claim 11, wherein the step 4.4 comprises the steps of: step 4.4.1: initializing the number t.sub.2=0 of the generated ferocious fishes; step 4.4.2: randomly generating the ferocious fish individuals according to the value ranges of the control variables to be optimized in the power grid:
W=[U.sub.G21,U.sub.G22, . . . U.sub.G2g,K.sub.T21,K.sub.T22, . . . ,K.sub.T2t,Q.sub.C21,Q.sub.C22, . . . ,Q.sub.C2c]=[w.sub.1,w.sub.2, . . . ,w.sub.s, . . . ,w.sub.S]; wherein, W is the value sequence of the control variables of the ferocious fish individuals, w.sub.s=x.sub.s max+rand(g)(x.sub.s maxx.sub.s min), U.sub.G2g, K.sub.T2t and Q.sub.C2c are respectively the generator terminal voltage amplitude, the transformer adjustable ratio, and the reactive power capacity of the reactive power compensation devices in the randomly generated ferocious fish individuals, and w.sub.s is the control variables of the randomly generated ferocious fish individuals; step 4.4.3: performing power flow calculation on the control variables in the randomly generated ferocious fish individuals in the step 4.4.2 by using a PQ decomposition method, if the power flow value of the current ferocious fish individual is converged, retaining the ferocious fish individual, generating a ferocious fish number t.sub.2=t.sub.2+1 of the current ferocious fish individual, and executing the step 4.4.4; or else, not retaining the ferocious fish individual, and returning to the step 4.4.2; and step 4.4.4: if t.sub.2N2, obtaining a current ferocious fish-swarm, and executing the step 4.5; or else, returning to the step 4.4.2.
14. The method for reactive power optimization of a power grid based on the double-fish-swarm algorithm according to claim 11, wherein the step 4.6 comprises the steps of: step 4.6.1: determining the distance D={D.sub.i1, D.sub.i2, . . . , D.sub.ij, . . . , D.sub.iN2} between the current i small fish individual X.sub.i and each of the ferocious fish individuals W.sub.j, wherein, D.sub.ij=X.sub.iW.sub.j, 1<j<N2, and N2 is the number of the ferocious fish individuals; step 4.6.2: if the distance D.sub.ijVisual1 between the current i th small fish individual X.sub.i and each of the ferocious fish individuals W.sub.j, executing the step 4.6.3; or else, executing the step 4.6.5; step 4.6.3: determining the total number n of the small fish individuals within the perception range of the i th small fish individual X.sub.i, to determine a protective clustering safe position X.sub.safe, and enabling the i th small fish individual X.sub.i to perform the protective clustering behavior so as to obtain an update value X.sub.inext1 of the i th small fish individual performing the protective clustering behavior, wherein the calculation formula for the protective clustering safe position X.sub.safe is as follows:
X.sub.safe=(X.sub.i1+X.sub.i2+ . . . +X.sub.in)/n, wherein, =(n+p)/n is an escape factor, p is the total number of the ferocious fish individuals within the perception range Visual1 of the small fish individual X.sub.i, and X.sub.in is the small fish individuals within the perception range of the small fish individual X.sub.i; the calculation formula for the update value X.sub.inext1 of the i th small fish individual performing the protective clustering behavior is as follows:
X.sub.inext1=X.sub.i+rand(g)Step1.sub.1(X.sub.safeX.sub.i)/X.sub.safeX.sub.i; step 4.6.4: judging whether the update value X.sub.inext1 of the current small fish individual performing the protective clustering behavior meets a constraint condition and the power flow is converged; if yes, recording the state of the update value X.sub.inext1, calculating the value of food concentration, updating the value in the bulletin board Y1, and executing the step 4.7; or else, returning to the step 4.6.3; step 4.6.5: determining the small fish individual X.sub.j with the maximum value FC.sub.j of food concentration in the small fish-swarm within the perception range of the i th small fish individual X.sub.i, and enabling the i th small fish individual X.sub.i to perform the rear-end behavior so as to obtain an update value X.sub.inext2 of the i th small fish individual X.sub.i performing the rear-end behavior, wherein the calculation formula for the update value X.sub.inext2 of the i th small fish individual X.sub.i performing the rear-end behavior is as follows:
X.sub.inext2=X.sub.i+rand(g)Step1.sub.1(X.sub.jX.sub.i)/X.sub.jX.sub.i; step 4.6.6: determining the total number n of the small fish individuals within the perception range of the i th small fish individual X.sub.i to determine the central position X.sub.c1 of the small fish individuals performing foraging clustering behavior, and enabling the i th small fish individual X.sub.i to perform the foraging clustering behavior so as to obtain an update value X.sub.inext3 of the i th small fish individual performing the foraging clustering behavior, wherein the calculation formula for the central position X.sub.c1 of the small fish individuals performing foraging clustering behavior is as follows:
X.sub.c1=(X.sub.i1+X.sub.i2+ . . . +X.sub.in)/n; the calculation formula for the update value X.sub.inext3 of the i th small fish individual performing the foraging clustering behavior is as follows:
X.sub.inext3=X.sub.i+rand(g)Step1.sub.1(X.sub.c1X.sub.i)/X.sub.c1X.sub.i; and step 4.6.7: judging whether the update value X.sub.inext2 of the current small fish individual performing the rear-end behavior and the update value X.sub.inext3 of the current small fish individual performing foraging clustering behavior meet constraint conditions and the power flow is converged; if yes, recording the state of the update value X.sub.inext2 and the state of the update value X.sub.inext3, calculating the value of food concentration, updating the value in the bulletin board Y1, and executing the step 4.7; or else, returning to the step 4.6.5.
15. The method for reactive power optimization of a power grid based on the double-fish-swarm algorithm according to claim 11, wherein the step 4.7 comprises the steps of: step 4.7.1: determining the distance L={L.sub.i1, L.sub.i2, . . . , L.sub.ij, . . . , L.sub.iN1} between the current ith ferocious fish individual W.sub.i and each of the small fish individuals X.sub.j, wherein L.sub.ij=W.sub.iX.sub.j, 1<j<N1, and N1 is the number of the small fish individuals; step 4.7.2: if the distance L.sub.ijVisual2 between the current ith ferocious and artificial fish individual W.sub.i and each of the small and artificial fish individuals X.sub.j, executing the step 4.7.3; or else, executing the step 4.7.6; step 4.7.3: determining the total number m of the small fish individuals within the perception range of the current ith ferocious fish individual W.sub.i, enabling the ith ferocious fish individual W.sub.i to perform the predation behavior with a central position X.sub.c2 of small fish-swarm within the perception range of the ith ferocious fish individual W.sub.i as the maximum value of current food concentration, and updating the ith ferocious fish individual W.sub.i as the central position X.sub.c2 of the small fish-swarm so as to obtain an update value X.sub.1next1=X.sub.c2 of the ith ferocious fish individual W.sub.i performing the predation behavior; the calculation formula for the central position X.sub.c2 of small fish-swarms within the perception range of the ith ferocious fish individual W.sub.i is as follows:
X.sub.c2=(X.sub.i1+X.sub.i2+ . . . +X.sub.im)/m; wherein X.sub.im is the small fishes within the perception range of the ferocious fish individual W.sub.i; step 4.7.4: enabling the ith ferocious fish individual W.sub.i to perform the tracking behavior with the central position X.sub.c2 of small fish-swarms within the perception range of the current ith ferocious fish individual W.sub.i as the maximum value of current food concentration so as to obtain an update value W.sub.inext2 of the ith ferocious fish individual W.sub.i performing the tracking behavior; the calculation formula for the update value W.sub.inext2 of the ith ferocious fish individual W.sub.i performing the tracking behavior is as follows:
W.sub.inext2=W.sub.j+rand(g)Step2.sub.2(X.sub.c2W.sub.i)/X.sub.c2W.sub.i; step 4.7.5: judging whether the update value W.sub.inext1 of the current ferocious fish individual W.sub.i performing the predation behavior and the update value W.sub.inext2 of the ferocious fish individual W.sub.i preforming the tracking behavior meet the constraint conditions and the power flow is converged; if yes, updating the value in the bulletin board Y2, and executing the step 4.7.6; or else, returning to the step 4.7.3; step 4.7.6: determining the total number r of the ferocious fish individuals within the perception range of the ith ferocious fish individual W.sub.i to determine the central position W.sub.c of the ferocious fish individuals preforming the clustering behavior, and enabling the i ferocious fish individual W.sub.i to perform the clustering behavior so as to obtain an update value W.sub.inext3 of the ith ferocious fish individual W.sub.i performing the clustering behavior, wherein the calculation formula for the central position W.sub.c of the fish individual performing the clustering behavior is as follows:
W.sub.c=(W.sub.i1+W.sub.i2+ . . . +W.sub.ir)/r; wherein, W.sub.ir is the ferocious fish individuals within the perception range of the ith ferocious fish individual W.sub.i; the calculation formula for the update value W.sub.inext3 of the ith ferocious fish individual W.sub.i performing the clustering behavior is as follows:
W.sub.inext1=W.sub.i+rand(g)Step2.sub.2(W.sub.cW.sub.i)/W.sub.cW.sub.i; and step 4.7.7: judging whether the update value W.sub.inext3 of the current ferocious fish individual W.sub.i performing the clustering behavior meets the constraint condition and the power flow is converged; if yes, recording the state of the update value W.sub.inext3, calculating the value of the food concentration, updating the value in the bulletin board Y2, and executing the step 4.8; or else, returning to the step 4.7.6.
16. The reactive power optimization system of a power grid based on the double-fish-swarm algorithm according to claim 1, wherein the specific process for optimizing the acquired initial data to be optimized in the current network based on the double-fish-swarm algorithm so as to obtain the optimal values of the control variables in the power grid comprises the steps of: generating a small fish-swarm and a ferocious fish-swarm in a random manner according to the upper limit value and the lower limit value of the control variables to be optimized in the power grid; initializing the optimal value of current food concentration, and displaying the initialized optimal value of current food concentration in a bulletin board: by using the reciprocal 1/F of the objective function of the mathematical model for reactive power optimization of a power grid as a value FC of food concentration, calculating the value of food concentration in the current network and of small fishes and ferocious fishes under an initial data condition, taking the maximum value of food concentration as the optimal value of current food concentration, displaying the optimal value of current food concentration in the bulletin board, and saving the state and the current value FC; enabling the small fish-swarm to act: enabling small fish individuals to perform a foraging clustering behavior, a rear-end behavior and a protective clustering behavior according to the distance between each of the small fish individuals and ferocious fish individuals as well as the distance between the small fish individuals in the small fish-swarm, and updating the value in the bulletin board Y1; enabling the ferocious fish-swarm to act: enabling each of the ferocious fish individuals to perform a predation behavior, a tracking behavior and a clustering behavior according to the distance between each of the ferocious fishes and the small fish individuals as well as the distance between the ferocious fish individuals in the ferocious fish-swarm, and updating the value in the bulletin board Y2; and when iteration times achieve the maximum iteration times, taking the greater values of food concentration displayed in the bulletin boards Y1 and Y2 as well as the state thereof as optimal results so as to obtain the optimal values of the control variables, namely the optimal values of the generator terminal voltage amplitudes, the transformer adjustable ratios, and the reactive power capacity of the reactive power compensation devices.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1)
(2)
(3)
(4)
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT
(5) The specific embodiments disclosed by the present invention are described in details below with reference to the accompanying drawings.
(6) In the embodiment, an IEEE10 node system as shown in
(7) The reactive power optimization system of a power grid based on a double-fish-swarm algorithm, as shown in
(8) the power grid state data acquiring module includes a power grid state data acquiring processor and a relay transmitter;
(9) the reactive power regulating module is a control terminal;
(10) the reactive power executing module includes generator terminal voltage regulators, transformer tap regulators and reactive power compensation regulators;
(11) the input end of the power grid state data acquiring processor is connected with a power grid; the output end of the power grid state data acquiring processor is connected with the input end of the relay transmitter; the output end of the relay transmitter is connected with the input end of the control terminal; the output end of the control terminal is connected with the input end of each of the generator terminal voltage regulators, the input end of each of the transformer tap regulators and the input end of each of the reactive power compensation regulators; the output end of each of the generator terminal voltage regulators is connected with a corresponding generator in the power grid; the output end of each of the transformer tap regulators is connected with a corresponding transformer in the power grid; and the output end of each of the reactive power compensation regulators is connected with a corresponding reactive power compensation device in the power grid.
(12) In the embodiment, the power grid state data acquiring processor of which the type is DCZL23-CL156N is used for acquiring the current network information of the power grid and judging whether the current network information meets the optimal state required by the power grid or not; and if the current network information cannot meet the optimal state required by the power grid, the current network information is transmitted to the relay transmitter, wherein the current network information includes power grid node information, branch information, generator information, transformer information and reactive power compensation device information.
(13) In the embodiment, the transformer information and the corresponding branch information obtained according to the IEEE10 node system as shown in
(14) TABLE-US-00001 TABLE 1 Two-winding transformer information and corresponding branch information Transformer Transformation Branch R X ratio 2-6 0 0.0625 1.0 3-8 0 0.0586 1.0 In the table, R is resistance and X is reactance.
(15) TABLE-US-00002 TABLE 2 Three-winding transformer information and corresponding branch information Transformer High-Voltage Medium-Voltage Low-Voltage Branch Impedance (R/X) Impedance (R/X) Impedance (R/X) 1-4-10 0.0011 0.085 0.0037 0.015 0.0015 0.095
(16) TABLE-US-00003 TABLE 3 Load parameters Node Number P.sub.load Q.sub.load 5 1.25 0.80 7 1.00 0.55 9 0.90 1.20 In the table, P.sub.load is the active power of the loads and Q.sub.load is the reactive power of the loads.
(17) TABLE-US-00004 TABLE 4 Branch parameters Branch First/Last Number Node Number R X B/2 1 4-5 0.0100 0.085 0.0440 2 5-6 0.0320 0.161 0.0765 3 6-7 0.0085 0.072 0.0373 4 7-8 0.0119 0.1008 0.0523 5 8-9 0.0390 0.170 0.0895 6 9-4 0.017 0.092 0.0395
(18) In the embodiment, the relay transmitter of which the type is RD 980 is used for transmitting the power grid node information, the branch information, the generator information, the transformer information and the reactive power compensation device information which are acquired by the power grid state data acquiring processor and required for reactive power optimization in the network information to the control terminal.
(19) The reactive power regulating module includes a processor that is configured to perform as a parameter acquiring unit, a double-fish-swarm algorithm-base reactive power optimization unit, and an optimization decision-making control unit, the parameter acquiring unit is connected to the reactive power optimization unit, and the reactive power optimization unit is connected to the optimization decision-making control unit;
(20) The parameter acquiring unit is used for acquiring the power grid node information, the branch information, the generator information, the transformer information and the reactive power compensation device information transmitted by the relay transmitter and used as initial data to be optimized in the current network;
(21) The double-fish-swarm algorithm-base reactive power optimization unit is used for establishing a mathematical model for reactive power optimization of a power system, and the acquired initial data to be optimized in the current network is optimized based on the double-fish-swarm algorithm, so that an optimal value of each of control variables in the power grid is obtained, wherein the control variables includes generator terminal voltage amplitudes, transformer adjustable ratios, and the reactive power capacity of reactive power compensation devices;
(22) In the embodiment, a mathematical model for the reactive power optimization of the power system is established in a penalty function manner by using an alpha-method, with the minimum network loss and the optimal voltage level as the optimization objectives, as shown in the formula (1):
(23)
(24) wherein, F is an objective function of the mathematical model for reactive power optimization of a power grid,
(25)
1<h<nl, nl is the total number of branches in a power grid system, G.sub.l is the electrical conductance of a connecting branch l, .sub. is the phase angle of a node , .sub.l is the phase angle of a node l, .sub.1 is the weight coefficient of an objective function for minimum network loss, .sub.2 is the weight coefficient of an objective function for optimum voltage, u.sub.i is the voltage of the node l, u.sub. is the voltage of the node , 1<<N.sub.d, N.sub.d is the total number of load nodes in the power grid system, 1<k<N.sub.g, N.sub.g is the number of generators in the power grid, V.sub. is the voltage of the node , V.sub..sup.sp is the set voltage value of the node , V.sub..sup.max is the maximum deviation value of voltage of the node , V.sub.,max is the upper limit value of the voltage of the node , V.sub.,min is the lower limit value of the voltage of the node , .sub.u is a penalty coefficient for voltage outside tolerances of the load nodes, .sub.q is a penalty coefficient for the reactive power output outside tolerances of the generators, Q.sub.k is the reactive power output of a generator node k, Q.sub.k,max is the upper limit value of the reactive power output of the generator node k, and Q.sub.k,min is the lower limit value of the reactive power output of the generator node k.
(26) In the embodiment, the special process for optimizing the initial data to be optimized in the current network based on the double-fish-swarm algorithm so as to obtain the optimal values of control variables in the power grid is as follows:
(27) generating a small fish-swarm and a ferocious fish-swarm in a random manner according to the upper limit value and the lower limit value of the control variables to be optimized in the power grid, wherein the control variables include generator terminal voltage amplitudes, transformer adjustable ratios, and reactive power capacity of reactive power compensation devices;
(28) generating a ferocious fish-swarm in a random manner according to the upper limit value and the lower limit value of the control variables to be optimized in the power grid;
(29) initializing the optimal value of current food concentration, and displaying the initialized optimal value of current food concentration in a bulletin board: by using the reciprocal of the objective function of the mathematical model for reactive power optimization of a power grid as a value FC of food concentration, by calculating the value of food concentration in the current network in the state of initial data to be optimized, the value of the food concentration in the state of the small and artificial fish individuals, and the value of the food concentration in the state of the ferocious and artificial fish individuals, taking the maximum value of food concentration as the optimal value of current food concentration, displaying the optimal value of current food concentration in the bulletin board, and saving the state and the optimal value FC of the current food concentration;
(30) in the embodiment, the calculation formula for the value FC of the the food concentration is as shown in the formula (2):
FC=1/F(2)
(31) enabling the small fish-swarm to act: enabling small fish individuals to perform a foraging clustering behavior, a rear-end behavior and a protective clustering behavior according to the distance between each of the small fish individuals and ferocious fish individuals as well as the distance between the small fish individuals in the small fish-swarm, and updating the value in the bulletin board Y1;
(32) enabling the ferocious fish-swarm to act: enabling each of the ferocious fish individuals to perform a predation behavior, a tracking behavior and a clustering behavior according to the distance between each of the ferocious fishes and the small fish individuals as well as the distance between the ferocious fish individuals in the ferocious fish-swarm, and updating Y2; and
(33) when the termination condition is achieved, taking the greater value of food concentration displayed in the bulletin boards Y1 and Y2 as well as the state thereof as optimal values so as to obtain the optimal values of the control variables, namely the optimal values of the generator terminal voltage amplitudes, the transformer adjustable ratios, and the reactive power capacity of the reactive power compensation devices.
(34) The optimization decision-making control unit is used for transmitting the optimal values of the control variables to the reactive power executing module.
(35) In the embodiment, the generator terminal voltage regulators of which the type is EQ1512 are used for regulating generator terminal voltages according to the optimal values of the generator terminal voltage amplitudes, obtained by the reactive power regulating module.
(36) In the embodiment, the transformer tap regulators of which the type is TSGC2J are used for regulating transformer taps according to the optimal values of the transformer adjustable ratios, obtained by the reactive power regulating module.
(37) In the embodiment, the reactive power compensation regulators of which the type is GLTSC are used for regulating the compensation capacity of reactive power compensators according to the optimal values of the reactive power capacity of the reactive power compensation devices, obtained by the reactive power regulating module.
(38) A method for reactive power optimization of a power grid by adopting the reactive power optimization system of a power grid based on the double-fish-swarm algorithm, as shown in
(39) step 1: acquiring the current network information of the power grid and judging whether the current network information meets the optimal state required by the power grid or not by the power grid state data acquiring processor; and if the current network information cannot meet the optimal state required by the power grid, transmitting the current network information to the relay transmitter;
(40) step 2: transmitting the power grid node information, the branch information, the generator information, the transformer information and the reactive power compensation device information which are acquired by the power grid state data acquiring processor and required for reactive power optimization in the network information to the control terminal, by the relay transmitter;
(41) step 3: acquiring the power grid node information, the branch information, the generator information, the transformer information and the reactive power compensation device information transmitted by the relay transmitter and used as initial data to be optimized in the current network, by the control terminal, wherein
(42) in the embodiment, the initial data to be optimized in the current network includes the node information, the branch information, the number of generators as well as the upper limit and the lower limit of each of the generator terminal voltage amplitudes U.sub.G, the number of transformers as well as the upper limit and the lower limit of each of the transformer adjustable ratios K.sub.t, the number of switching groups of the reactive power compensation devices, and the upper limit and the lower limit of each of the reactive power capacity Q.sub.c of the reactive power compensation devices;
(43) step 4: establishing the mathematical model for the reactive power optimization of the power system through the control terminal, optimizing the acquired initial data to be optimized in the current network based on the double-fish-swarm algorithm so as to obtain the optimal values of the control variables in the power grid, and transmitting the optimal values of the control variables to the reactive power executing module, as shown in
(44) step 4.1: establishing a mathematical model for the reactive power optimization of the power system by using an alpha-method in a weighted sum method, with the minimum network loss and the optimal voltage level as the optimization objectives;
(45) step 4.2: setting parameters for the double-fish-swarm algorithm: the size N1 of the small fish-swarm, the size N2 of the ferocious fish-swarm, the perception range Visual1 of the small fish-swarm, the perception range Visual2 of the ferocious fish-swarm, the shift step length Step1 of the small fish-swarm, the shift step length Step2 of the ferocious fish-swarm, the crowding factor .sub.1 of the small fish-swarm, the crowding factor .sub.2 of the ferocious fish-swarm, the attempt times try_number1 of the small fish-swarm, the attempt times try_number2 of the ferocious fish-swarm, and the maximum iteration times K>0,
(46) and in the embodiment, the maximum iteration times K is 100;
(47) step 4.3: generating a small fish-swarm in a random manner according to the upper limit value and the lower limit value of the control variables to be optimized in the power grid:
(48) step 4.3.1: initializing the number t.sub.1=0 of the generated small fishes;
(49) step 4.3.2: generating small and artificial fish individuals in a random manner according to the upper limit value and the lower limit value of the control variables to be optimized in the power grid, as shown in the formula (3):
X=[U.sub.G11,U.sub.G12, . . . U.sub.G1g,K.sub.T11,K.sub.T12, . . . ,K.sub.T1t,Q.sub.C11,Q.sub.C12, . . . ,Q.sub.C1c]=[x.sub.1,x.sub.2, . . . x.sub.s, . . . ,X.sub.S](3)
(50) wherein, X is the value sequence of the control variables of the small and artificial fish individuals, x.sub.s=x.sub.s min+rand(g)(x.sub.x maxx.sub.s min), rand(g) is a random number of an interval (0, 1), x.sub.s min the lower limit value of the corresponding control variables, x.sub.s max is the upper limit value of the corresponding control variable, U.sub.G1g is the generator terminal voltage amplitude in the randomly generated small and artificial fish individuals, K.sub.T1t is the transformer adjustable ratio in the randomly generated small and artificial fish individuals, Q.sub.C1c is the reactive power capacity of the reactive power compensation devices in the randomly generated small and artificial fish individuals, x.sub.s is the control variables of the randomly generated small fish individuals, 1sS, S=N.sub.g+N.sub.t+N.sub.c, N.sub.g is the number of the generators in the power grid, N.sub.t is the number of the transformers in the power grid, and N.sub.c is the number of the reactive power compensation devices in the power grid;
(51) step 4.3.3: performing power flow calculation on the control variables in the randomly generated small fish individuals in the step 4.3.2 by using a PQ decomposition method. If the power flow value of the current small fish individual is converged, retaining the small fish individual, enabling t.sub.1=t.sub.1+1, and executing the step 4.3.4; or else, not retaining the small fish individual, and returning to the step 4.3.2;
(52) step 4.3.4: if the number t.sub.1 of the current generated small fish individuals achieves the scale N1 of the small fish-swarm, obtaining a current small fish-swarm, and executing the step 4.4; or else, returning to the step 4.3.2;
(53) step 4.4: generating a ferocious fish-swarm in a random manner according to the upper limit value and the lower limit value of the control variables to be optimized in the power grid:
(54) step 4.4.1: initializing the number t.sub.2=0 of the generated ferocious fishes;
(55) step 4.4.2: generating ferocious and artificial fish individuals in a random manner according to the upper limit value and the lower limit value of the control variables to be optimized in the power grid, as shown in the formula (4):
W=[U.sub.G21,U.sub.G22, . . . ,U.sub.G2g,K.sub.T21,K.sub.T22, . . . ,K.sub.T2t,Q.sub.C21,Q.sub.C22, . . . ,Q.sub.C2c]=[w.sub.1,w.sub.2, . . . ,w.sub.s, . . . ,w.sub.S] (4)
(56) step 4.4.3: performing power flow calculation on the control variables in the randomly generated small ferocious individuals in the step 4.4.2 by using a PQ decomposition method. If the power flow value of the current ferocious fish individual is converged, retaining the ferocious fish individual, enabling t.sub.2=t.sub.2+1, and executing the step 4.4.4; or else, not retaining the ferocious fish individual, and returning to the step 4.4.2;
(57) step 4.4.4: if the number t.sub.2 of the current ferocious fish individuals achieves the scale N2 of the ferocious fish-swarm, obtaining a current small fish-swarm, and executing the step 4.5; or else, returning to the step 4.4.2;
(58) step 4.5: initializing the current iteration times k=0, initializing the optimal value of current food concentration, and displaying the initialized optimal value of current food concentration in a bulletin board: by using the reciprocal of the objective function of the mathematical model for reactive power optimization of a power grid as a value FC of food concentration, by calculating the value of food concentration in the current network in the state of initial data to be optimized, the value of the food concentration in the state of the small and artificial fish individuals, and the value of the food concentration in the state of the ferocious and artificial fish individuals, taking the maximum value of food concentration as the optimal value of current food concentration, displaying the optimal value of current food concentration in the bulletin board, and saving the state and the optimal value FC of the current food concentration;
(59) step 4.6: enabling the small fish-swarm to act: enabling small fish individuals to perform a foraging clustering behavior, a rear-end behavior and a protective clustering behavior, and updating the value in the bulletin board Y1;
(60) step 4.6.1: determining the distance between the current i th small and artificial fish individual X.sub.i and each of the ferocious and artificial fish individuals W.sub.j, as shown in the formula (5):
D={D.sub.i1,D.sub.i2, . . . ,D.sub.ij, . . . ,D.sub.iN2}(5)
(61) wherein, D.sub.ij=X.sub.iW.sub.j, 1<j<N2, N2 is the number of the ferocious and artificial fish individuals;
(62) step 4.6.2: if the distance D.sub.ijVisual1 between the current i th small and artificial fish individual X.sub.i and each of the ferocious and artificial fish individuals W.sub.j, executing the step 4.6.3; or else, executing the step 4.6.5;
(63) step 4.6.3: determining the total number n of the small fish individuals within the perception range of the i th small fish individual X.sub.i, to determine a protective clustering safe position X.sub.safe, and enabling the i th small fish individual X.sub.i to perform the protective clustering behavior so as to obtain an update value X.sub.inext1 of the i th small fish individual performing the protective clustering behavior, wherein
(64) in the embodiment, the small fish individuals X.sub.i within the perception range of the i th small fish individual, is determined, namely the small fish individuals X.sub.in meeting the formula (6):
d.sub.in=X.sub.iX.sub.inVisual1(6)
(65) in the embodiment, the calculation formula for the protective clustering safe position X.sub.safe is as shown in the formula (7):
X.sub.safe=(X.sub.i1+X.sub.i2+ . . . +X.sub.in)/n(7)
(66) wherein =(n+p)/n is escape factors, p is the total number of the ferocious fish individuals within the perception range Visual1 of the small fish individual X.sub.i, X.sub.in is the number of the small fish individuals within the perception range of the small fish individuals X.sub.i, the escape factors are introduced to enable the small fish-swarm to be collected to the safe position away from the ferocious fish-swarm as possible when performing the protective clustering behavior;
(67) in the embodiment, the update value X.sub.inext1 of the i th small fish individual performing the protective clustering behavior is shown as the formula (8):
X.sub.inext1=X.sub.i+rand(g)Step1.sub.1(X.sub.safeX.sub.i)/X.sub.safeX.sub.i(8)
(68) step 4.6.4: judging whether the update value X.sub.inext1 of the current small and artificial fish individual performing the protective clustering behavior meets a constraint condition and the power flow is converged; if yes, recording the state of the update value X.sub.inext1, calculating the value of food concentration, updating the value in the bulletin board Y1, and executing the step 4.7; or else, returning to the step 4.6.3;
(69) in the embodiment, the constraint conditions are the corresponding number of generators as well as the upper limit and the lower limit of each of the generator terminal voltage amplitudes U.sub.G, the number of transformers as well as the upper limit and the lower limit of each of the transformer adjustable ratios K.sub.i, the number of switching groups of the reactive power compensation devices, and the upper limit and the lower limit of each of the reactive power capacity Q.sub.c of the reactive power compensation devices;
(70) step 4.6.5: determining the small fish individual X.sub.j with the maximum value FC.sub.j of food concentration in the small fish-swarm within the perception range of the i th small fish individual X.sub.i, and enabling the i th small fish individual X.sub.i to perform the rear-end behavior so as to obtain an update value X.sub.inext2 of the i th small fish individual X.sub.i performing the rear-end behavior, wherein
(71) in the embodiment, the update value X.sub.inext2 of the i th small and artificial fish individual X.sub.i performing the rear-end behavior is shown as the formula (9):
X.sub.inext2=X.sub.i+rand(g)Step1.sub.1(X.sub.jX.sub.i)/X.sub.jX.sub.i(9)
(72) step 4.6.6: determining the total number n of the small fish individuals within the perception range of the i th small fish individual X.sub.i to determine the central position X.sub.c1 of the small fish individuals performing foraging clustering behavior, and enabling the i th small fish individual X.sub.i to perform the foraging clustering behavior so as to obtain an update value X.sub.inext3 of the i th small fish individual performing the foraging clustering behavior, wherein
(73) in the embodiment, the calculation formula for the central position X.sub.c1 of the fish individuals performing the foraging clustering behavior is as shown in the formula (10):
X.sub.c1=(X.sub.i1+X.sub.i2+ . . . +X.sub.in)/n(10)
(74) in the embodiment, the update value X.sub.inext3 of the i th small fish individual performing the foraging clustering behavior is shown as the formula (11):
X.sub.inext3=X.sub.i+rand(g)Step1.sub.1(X.sub.c1X.sub.i)/X.sub.c1X.sub.i(11)
(75) step 4.6.7: judging whether the update value X.sub.inext2 of the current small and artificial fish individual and the update value X.sub.inext3 performing the foraging clustering behavior of the small and artificial fish individual performing the foraging clustering behavior meet constraint conditions and the power flow is converged; if yes, recording the state of the update value X.sub.inext2 and the state of the update value X.sub.inext3, calculating the value of food concentration, updating the value in the bulletin board Y1, and executing the step 4.7; or else, returning to the step 4.6.5.
(76) step 4.7: enabling the ferocious fish-swarm to act: enabling each of the ferocious fish individuals to perform a predation behavior, a tracking behavior and a clustering behavior according to the distance between each of the ferocious fish individuals and the small fish individuals in the ferocious fish-swarm as well as the distance between of the ferocious fish individuals, and updating the value Y2 of the food concentration of the ferocious fish-swarm performing the corresponding behavior in the bulletin board;
(77) step 4.7.1: determining the distance between the current ith ferocious fish individual W.sub.i and each of the small fish individuals X.sub.j, as shown in the formula (12):
L={L.sub.i1,L.sub.i2, . . . ,L.sub.ij, . . . ,L.sub.iN1}(12)
(78) wherein, L.sub.ij=W.sub.iX.sub.j, 1<j<N1, N1 is the number of small fish individuals;
(79) step 4.7.2: if L.sub.ijVisual2, executing the step 4.7.3; or else, executing the step 4.7.6;
(80) step 4.7.3: determining the total number m of the small fish individuals within the perception range of the current ith ferocious fish individual W.sub.i, enabling the ith ferocious fish individual W.sub.i to perform the predation behavior with a central position X.sub.c2 of small fish-swarm within the perception range of the ith ferocious fish individual W.sub.i as the maximum value of current food concentration, and updating the ith ferocious fish individual W.sub.i as the central position X.sub.c2 of the small fish-swarm so as to obtain an update value W.sub.inext1=X.sub.c2 of the ith ferocious fish individual W.sub.i performing the predation behavior, wherein
(81) in the embodiment, the calculation formula for the central position X.sub.c2 of small fish-swarms within the perception range of the ith ferocious fish individual W.sub.i is shown as the formula (13):
X.sub.c2=(X.sub.i1+X.sub.i2+ . . . +X.sub.im)/m(13)
(82) wherein, X.sub.im is the small fish individuals within the perception range of the ferocious fish individual W.sub.i;
(83) step 4.7.4: enabling the ith ferocious fish individual W.sub.i to perform the tracking behavior with the central position X.sub.c2 of small fish-swarms within the perception range of the current ith ferocious fish individual W.sub.i as the maximum value of current food concentration so as to obtain an update value W.sub.inext2 of the ith ferocious fish individual W.sub.i performing the tracking behavior;
(84) in the embodiment, the update value W.sub.inext2 of the ith ferocious fish individual W.sub.i performing the tracking behavior is shown as the formula (14):
W.sub.inext2=W.sub.i+rand(g)Step2.sub.2(X.sub.c2W.sub.i)/X.sub.c2W.sub.i(14)
(85) step 4.7.5: judging whether the update value W.sub.inext1 of the current ferocious fish individual W.sub.i performing the foraging behavior and the update value W.sub.inext2 of the ferocious fish individual W.sub.i performing the tracking behavior meet constraint conditions and the power flow is converged; if yes, recording the state of the update value W.sub.inext1 and the state of the update value W.sub.inext2, calculating the value of the food concentration, updating the value in the bulletin board Y2, and executing the step 4.7.6; or else, returning to the step 4.7.3;
(86) step 4.7.6: determining the total number r of the ferocious fish individuals within the perception range of the ith ferocious fish individual W.sub.i to determine the central position W.sub.c of the ferocious fish individuals preforming the clustering behavior, and enabling the ith ferocious fish individual W.sub.i to perform the clustering behavior so as to obtain an update value W.sub.inext3 of the ith ferocious fish individual W.sub.i performing the clustering behavior, wherein
(87) in the embodiment, the calculation formula for the central position W.sub.c of the fish individuals performing the foraging clustering behavior is as shown in the formula (15):
W.sub.c=(W.sub.i1+W.sub.i2+ . . . +W.sub.ir)/r(15)
(88) wherein, W.sub.ir is the ferocious fish individuals within the perception range of the ith ferocious fish individual W.sub.i,
(89) in the embodiment, the update value W.sub.inext3 of the ith ferocious fish individual W.sub.i performing the clustering behavior is shown as the formula (16):
W.sub.inext1=W.sub.i+rand(g)Step2.sub.2(W.sub.cW.sub.i)/W.sub.cW.sub.i(16)
(90) step 4.7.7: judging whether the update value W.sub.inext3 of the current ferocious fish individual W.sub.i performing the clustering behavior meets the constraint condition and the power flow is converged; if yes, recording the state of the update value W.sub.inext3, calculating the value of the food concentration, updating the value in the bulletin board Y2, and executing the step 4.8; or else, returning to the step 4.7.6.
(91) step 4.8: judging whether the current iteration times k reach the maximum iteration times K, if yes, executing the step 4.9; or else, k=k+1, and returning to the step 4.6;
(92) step 4.9: taking the greater value of food concentration displayed in the bulletin boards Y1 and Y2 as well as the state thereof as optimal values so as to obtain the optimal values of the control variables, namely the optimal values of the generator terminal voltage amplitudes, the transformer adjustable ratios, and the reactive power capacity of the reactive power compensation devices; and
(93) step 4.10: transmitting the optimal values of the control variables to the reactive power executing module; and
(94) step 5: regulating generator terminal voltages by using the generator terminal voltage regulators according to the optimal values of the generator terminal voltage amplitudes, regulating the transformer taps by using the transformer tap regulators according to optimal values of the transformer adjustable ratios, and regulating the compensation capacity of the reactive power compensators by using the reactive power compensation regulators according to the optimal values of the reactive power capacity of the reactive power compensation devices.
(95) In the embodiment, a particle swarm algorithm, a genetic algorithm and the double-fish-swarm algorithm in the present invention are respectively used to perform reactive power optimization simulation on the system. 100 iterations and 10 simulation experiments are performed for each of the three algorithms. The comparison of simulation results are as shown in Table 5 and Table 6:
(96) TABLE-US-00005 TABLE 5 Comparison of node voltages before and after optimization by using particle swarm algorithm, the genetic algorithm and the double-fish-swarm algorithm in the present invention Voltage after Optimization Voltage after Voltage after by Particle Optimization Optimization by Node Voltage before Swarm by Genetic Double-Fish-Swarm Name Optimization Algorithm Algorithm Algorithm 1 1.030 1.045 1.056 1.100 2 1.025 1.036 1.040 1.100 3 1.025 1.036 1.040 1.100 4 0.980 1.011 1.015 1.031 5 1.000 1.022 1.055 1.127 6 1.017 1.045 1.056 1.142 7 0.991 1.023 1.061 1.133 8 1.007 1.015 1.026 1.047 9 0.910 1.012 1.035 1.074 10 1.003 1.015 1.020 1.037
(97) TABLE-US-00006 TABLE 6 Comparison of system network loss before and after optimization by using the particle swarm algorithm, the genetic algorithm and the double-fish-swarm algorithm in the present invention Network Network Network Losses after Losses after Losses after Network Optimization by Optimization Optimization by Loss before Particle Swarm by Genetic Double-Fish-Swarm Optimization Algorithm Algorithm Algorithm (MW) (MW) (MW) (MW) 6.482560 5.852618 5.224537 4.607137
(98) According to Table 5 and Table 6, through the comparative analysis after compensation, compared with the particle swarm algorithm and the genetic algorithm, the double-fish-swarm algorithm achieves a better effect in the process of reactive power optimization, the voltage quality is obviously improved, and the network loss is reduced to the maximum, so that the reactive power distribution is improved, the safety and stability of the system are guaranteed, and the results show that the double-fish-swarm algorithm has good feasibility and practicability.