CN104063808B - Trans-provincial power transmission cascade hydropower station group peak-shaving dispatching two-phase search method - Google Patents

Abstract

The invention belongs to the field of hydropower dispatching operation, and discloses a trans-provincial power transmission cascade hydropower station group peak-shaving dispatching two-phase search method. The characteristics of load complementation between receiving-end power grids can be fully utilized and different peak power demands of multiple receiving-end power grids can be taken into consideration to realize reasonable allocation of quality resources between the multiple power grids. The technical scheme is as follows: the trans-provincial power transmission cascade hydropower station group peak-shaving goal is realized in two phases according to the characteristics of the trans-provincial power transmission peak-shaving issue; in the optimization calculation process, load curves with great difference of a plurality of provincial power grids are reconstructed into an overall load by adoption of a load reconstruction strategy, and the hydropower station power output process is determined on the basis by adoption of a variable-scale optimization method; and then the allocation process of output power of the station between the receiving-end power grids is optimized by adoption of an exterior point search method, an initial solution is generated through successive load shedding, and the peak power sizes of the receiving-end power grids are corrected in an iterative mode by taking load level as heuristic information to meet the peak-shaving demands of the multiple power grids. Peak-shaving demands of multiple provincial power grids can be taken into consideration, and a peak-shaving result better than actual dispatching can be obtained. The method is feasible.

Description

One power transmission Hydropower Stations peak regulation scheduling transprovincially two benches searching method

Technical field

The present invention relates to water power management and running field, particularly one power transmission Hydropower Stations peak regulation scheduling transprovincially two benches searching method.

Technical background

Through the basin step rolling development of more than 20 years, China has formed especially big basin and master stream Hydro Power Systems with Cascaded Reservoirs general layout, the Wujiang River, large Hydropower Base master stream, Hongsuihe River two step are substantially all built up, Jinsha jiang River middle and lower reaches, Middle-lower Reaches of The Lancang River, Yalongjiang River, large Hydropower Base especially big basin, Dadu River four Hydropower Stations also entered and be completed and concentrated operation period comprehensively, and Chinese water power has marched toward extensive transregionally to dissolve scheduling and operation phase transprovincially.Simultaneously, that " gives direct current again " with south electric network " two cross direct current " and national grid puts into operation, water power transprovincially transregional transmission of electricity scale achieves swift and violent growth, has become the important electric power resource of the developed area such as Yangtze River Delta, Delta of the Pearl River receiving end electrical network.Different from on-site elimination mode, present stage single basin step need to transmit electric power for multiple provincial power network or regional power grid simultaneously, meet the electricity consumption particularly peak power demands in these receiving end areas.How to utilize load complementary characteristic between each receiving end electrical network, give full play to the adjusting function of these high-quality hydroelectric resources, alleviate the provincial power network ubiquitous peak regulation pressure such as the regional power grid such as East China, south and Shanghai, Jiangsu, Zhejiang, Guangdong better and become theory and practice problem urgently to be resolved hurrily in China's electrical network and Model for Cascade Hydroelectric Stations, its Study on Problems is the key point of China's flood TV university wide resource configuration.

Current this respect research there is not yet pertinent literature report both at home and abroad.Transprovincially power transmission peaking problem relates to multiple electrical network, need the more complicated interconnected constraint of electric power of process and optimization object function, face complex conditions and the difficult problem such as distributing electric power between multiple network load of very different and net, Optimization Solution difficulty is larger, needs effective method for solving.Achievement of the present invention with master stream, CHINA SOUTHERN POWER Hongsuihe River Hydropower Stations transprovincially power transmission problem for background, taken into account many peaking demand of power grids that Hydropower Stations faces in power transmission transprovincially, there is important popularizing value.

Summary of the invention

The technical problem to be solved in the present invention is to provide one power transmission Hydropower Stations peak regulation scheduling transprovincially two benches searching method, the load complementary characteristic between receiving end electrical network can be made full use of, take into account the differentiation peak power demand of multiple receiving end electrical network, realize the Reasonable Regulation And Control configuration of high-quality resource between multiple electrical network.

Technical scheme of the present invention is: present invention is disclosed one power transmission Hydropower Stations peak regulation scheduling transprovincially two benches searching method, comprise load reconstruct and output distribution two large divisions, complete power transmission peak regulation scheduling process transprovincially according to following step (1)-(8):

(1) initial calculation condition, comprises hydroelectric station operation condition and constraint, and each network load process and by electrical demand;

(2) power station numbering m=1 is made;

(3) what build m power station faces workload demand { C ₁, C ₂..., C _t, wherein C ₁, C ₂..., C _tbe respectively the period 1,2 ..., the load value of T;

(4) with { C ₁, C ₂..., C _tbe workload demand, under given objective function and constraint condition, adopt DFP-method to determine day 96 processes of exerting oneself in m power station;

(5) under given controlled condition, outer point searching method is adopted to determine the output distribution process of m power station at G electrical network, i.e. { p _{m, 1,1}, p _{m, 1,2}..., p _{m, 1, T}, { p _{m, 2,1}, p _{m, 2,2}..., p _{m, 2, T}..., { p _{m, G, 1}, p _{m, G, 2}..., p _{m, G, T}, wherein, { p _{m, 1,1}, p _{m, 1,2}..., p _{m, 1, T}, { p _{m, 2,1}, p _{m, 2,2}..., p _{m, 2, T}..., { p _{m, G, 1}, p _{m, G, 2}..., p _{m, G, T}represent that m power station is the 1st, 2 respectively ..., the 1st, 2 of G electrical network ..., the load value of T period;

(6) each electrical network residue load is recalculated;

(7) make m=m+1, if m is less than or equal to M, then skip to step (3); Otherwise go to step (8);

(8) calculate the desired value of current solution, judge whether target is improved, if so, using the initial solution of current solution as next round iteration, and skip to step (2); Otherwise, judge whether to reach to maximum iteration time, if so, algorithm search stops, and exports optimal result, otherwise, skip to step (2).

The present invention contrasts prior art following beneficial effect: one of the present invention power transmission Hydropower Stations peak regulation scheduling transprovincially two benches searching method, realizes power transmission Hydropower Stations peak regulation target transprovincially in two stages; In optimization computation process, first adopt load reconstruction strategy that the load curve of multiple provincial power network very different is reconstructed into a total load, exerting oneself process in power station to adopt DFP-method to determine on this basis; Then outer point searching method is adopted to optimize the assigning process of output of power station between receiving end electrical network, in the assignment procedure, cutting load is adopted successively to generate initial solution, and using load height as heuristic information, the peak power size of each receiving end electrical network of iterated revision, to meet multiple peak load regulation network requirement.Contrast prior art, the present invention can take into full account load complementary characteristic between each receiving end electrical network, takes into account the electricity needs of many receiving ends electrical network in transprovincially power transmission peak regulation, realizes the Reasonable Regulation And Control configuration of high-quality resource between many electrical networks.

Accompanying drawing explanation

Fig. 1 is that the inventive method totally solves block diagram.

Fig. 2 is single power station power transmission two provincial power network schematic diagram.

Fig. 3 is distributing electric power problem initial solution schematic diagram between the net that obtains of outer point searching method.

Fig. 4 (a) be in the Hydropower Stations optimum results of master stream, Hongsuihe River Guangxi Power Grid by electric process.

Fig. 4 (b) be in the Hydropower Stations optimum results of master stream, Hongsuihe River Guangdong Power Grid by electric process.

Embodiment

Below in conjunction with drawings and Examples, the invention will be further described.

The extensive transregional conveying transprovincially of water power is to regulate and control high-quality resource between the multiple electrical network of configuration to alleviate day by day serious peak regulation pressure, this problem has obviously different from the water power scheduling problem of single electrical network in the past, need the more complicated interconnected constraint of electric power of process and optimization object function, face complex conditions and the difficult problem such as distributing electric power between multiple network load of very different and net, Optimization Solution difficulty is larger.For how making full use of load complementary characteristic between each receiving end electrical network, taking into account the electricity needs of many receiving ends electrical network in transprovincially power transmission peak regulation, rarely having bibliographical information both at home and abroad.The present invention discloses one power transmission Hydropower Stations peak regulation scheduling transprovincially two benches searching method, will transprovincially be decomposed into two crucial subproblems (one is consider multiple provincial power network workload demand, optimizes the process of exerting oneself determining power station by power transmission peaking problem; Two is that stationary electric power plant exerts oneself process and each electrical network by electricity, determine the distributing electric power process of power station between receiving end electrical network), and for these two subproblems, adopt stage by stage based on load reconstruct DFP-method and outer point searching method solve respectively, with coordinate Hydropower Stations exert oneself and each electrical network by electric process, meet multiple peak load regulation network requirement.

The present invention is reconstructed by load and output distribution dual stage process forms.First stage adopts load reconstruction strategy that the load curve of multiple provincial power network very different is reconstructed into a total load, and exerting oneself process in power station to adopt DFP-method to determine on this basis; Subordinate phase proposes outer point searching method and optimizes the assigning process of output of power station between receiving end electrical network, cutting load is adopted successively to generate initial solution, and using load height as heuristic information, the peak power size of each receiving end electrical network of iterated revision, to meet multiple peak load regulation network requirement.

The concrete operation method in each stage is achieved according to following process (a)-(g):

A DFP-method that () reconstructs based on load

DFP-method based on load reconstruct mainly solves first subproblem, and the Hydropower Stations generating optimization problem namely under multiple network load demand, solution procedure is mainly divided into following two steps.

The first step: reconstruct power station faces workload demand

In multiple provincial power network load curve is between the order of magnitude, peak and low ebb number, peak are current etc., usually there is very big-difference, directly determines that the peaking operation plan in power station is very difficult.So, adopt load reconstruction strategy that many load curves are reconstructed into a total load herein, to reduce output of power station search difficulty; For making total load retain each network load Variation Features as far as possible, especially by electricity proportion compared with the load peak of bulk power grid and low ebb information, in literary composition using power transmission amount ratio as weight, adopt method for normalizing to carry out load reconstruct, see following formula.

$C_t=\underset{g=1}{\overset{G}{\mathrm{\Σ}}}{R}_{m,g}\×C_{g,t}\×\frac{C_{g,t}-C_{g,\min }}{C_{g,\max }-C_{g,\min }},$ Wherein, 1≤t≤T

In formula: C _tfor the residue load in the t period of electrical network, MW; C _g,tfor the load value of g electrical network t period; C _{g, max}for the load maximal value of g electrical network; C _{g, min}for the load minimum value of g electrical network; G is electrical network numbering, and 1≤g≤G, G is electrical network sum; T is period sequence number, and 1≤t≤T, T is period sum.

Second step: determine output of power station process

Based on the total load that the first step obtains, this subproblem can be converted into GROUP OF HYDROPOWER STATIONS generating optimization problem in the single electrical network that is easy to solve, then DFP-method is adopted to solve, see document: the GROUP OF HYDROPOWER STATIONS short-term mutative scale Optimization Scheduling under complicated period coupling constraints, Shen Jianjian etc., Proceedings of the CSEE .2014,34 (1): 87-95.

The outer point searching method of distributing electric power between (b) net

Second subproblem is the fixing output of power station process obtained on last stage, utilizes the load difference between provincial power network, optimizes its assigning process at each receiving end electrical network, meets the peak regulation demand of these electrical networks.Rarely have the research about this respect to report both at home and abroad at present, institute, in the process of the present invention according to problematic features, proposes a kind of outer point searching method and carries out problem solving.The method take peak load regulation network as target, and using load height as heuristic information, from an initial infeasible solutions, the allocated size of iterated revision day part output of power station between receiving end electrical network, until converge to the feasible solution of problem.Because method search procedure progressively converges to feasible solution from an infeasible solutions, therefore be referred to as outer point searching method, its solution procedure comprises following two major parts:

(1) initial solution is determined.Stationary electric power plant exert oneself process condition under, adopt following formula can obtain each electrical network by electricity, so with by electricity for control objectives, face by power station based on load, adopt successively cutting load method to determine the powered process of exerting oneself of the ideal of each electrical network.

E _m,g＝E _m×R _m,g

In formula: E _m,gfor m power station is to the day power transmission amount of g electrical network; E _mrepresent the total electricity of power station m in schedule periods; R _m,gfor the electricity ratio of g electrical network is sent in m power station,

Without loss of generality, Fig. 2 gives the examples of problems of single power station power transmission two provincial power networks simultaneously adopting said method to obtain, due in search procedure to reach given target electricity for the condition of convergence, retrain so the initial solution in Fig. 3 meets electrical network by electricity demand for control, but violate period power balance constraint, namely the powered sum of exerting oneself of all electrical networks corresponding period is not equal to this period output of power station, as can be seen from Figure 3.

(2) based on the initial solution obtained by the first step, reduce the low remaining lotus period and exert oneself according to preferential and preferentially increase the rule that the high remaining lotus period exerts oneself, redistribute the day part electrical network not meeting period power balance constraint poweredly to exert oneself, until all periods all reach power balance.

In short term scheduling runs; the peak regulation task of system is born in power station mostly; object utilizes the advantages such as Hydropower Unit start and stop are rapid, Ramp Rate is fast; carry out tracking to electrical network peak value and load process to regulate; to reduce the remaining load peak-valley difference of electrical network, and level and smooth as far as possible remaining load process, and then the frequent startup-shutdowns of power supply such as the poor thermoelectricity of adjusting function can be reduced; save generating to consume, realize power grid security, energy-conservation, economical operation.So this method adopts electrical network remaining load mean square deviation minimum expression Hydropower Stations peak regulation optimization aim, sees following formula:

$\min F=\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{((C_t-\underset{m=1}{\overset{M}{\mathrm{\Σ}}}{p}_{m,t})-\underset{t=1}{\overset{T}{\mathrm{\Σ}}}(C_t-\underset{m=1}{\overset{M}{\mathrm{\Σ}}}{p}_{m,t})/T)}^2$

In formula: F is the remaining load mean square deviation of electrical network; P _m,tfor m power station exerting oneself in the t period, MW; C _tfor the residue load in the t period of electrical network, MW; 1≤m≤M, M is power station sum; 1≤t≤T, T is period sum.

For transprovincially power transmission Hydropower Stations, single power station needs the peak regulation requirement simultaneously meeting multiple provincial power network, this kind of GROUP OF HYDROPOWER STATIONS scheduling problem is typical multi-objective optimization question, and need the remaining load mean square deviation simultaneously minimizing multiple electrical network, its objective function can be expressed as:

$\min F_g=\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{((C_{g,t}-\underset{m=1}{\overset{M}{\mathrm{\Σ}}}{p}_{m,g,t})-\underset{i=1}{\overset{T}{\mathrm{\Σ}}}(C_{g,t}-\underset{m=1}{\overset{M}{\mathrm{\Σ}}}{p}_{m,g,t})/T)}^2$

g＝1，2，...，G

In formula: F _gfor the remaining load mean square deviation of g electrical network; P _{m, g, t}for exerting oneself of g electrical network is delivered in the t period in m power station, MW; C _g,tfor the load value of g electrical network t period; 1≤g≤G, G is electrical network sum.

For reducing problem solving difficulty, adopting conventional target weight method that multi-objective optimization question is converted into single-objective problem, seeing following formula:

$\min F=\underset{g=1}{\overset{G}{\mathrm{\Σ}}}{w}_g\×\frac{F_g}{C_{g,\max }^2}$

In formula: w _gfor the target weight of g electrical network, consider the equal importance of each peaking demand of power grid, in literary composition, get w _g=1/G; C _{g, max}for the load maximal value of g electrical network, the object introducing this variable is normalized the residue load in objective function, causes optimization solution deviation to avoid each network load level difference.

The constraint condition expression formula solving power transmission peaking problem demand fulfillment is transprovincially as follows:

(1) water balance constraint

V _m,t+1＝V _m,t+3600×(Q _m,t-q _m,t-Qd _m,t)Δ _t

In formula: V _{m, t+1}for power station m is at the storage capacity of period t+1; V _m,tfor power station m is at the storage capacity of period t; Q _m,tfor power station m is at the reservoir inflow of period t, k is total number in power station, m upstream, power station; after during for considering stagnant, power station, upstream k flows into the flow of power station m at period t; Qn _m,tfor power station m is at the interval flow of period t; Qd _m,tfor power station m abandons discharge at period t; q _m,tfor power station m is at the generating flow of period t; Δ _tfor t period hourage.

(2) power station day controls electrical demand

E _m＝E' _m

In formula: E _m, E' _mrepresent the total electricity of power station m in schedule periods and given demand electricity respectively.

(3) electrical network is by electricity demand for control

$\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{p}_{m,g,t}\×{\mathrm{\Δ}}_t=E_m\×R_{m,g}$

In formula: R _m,gfor the electricity ratio of g electrical network is sent in m power station,

(4) period power balance constraint

$\underset{g=1}{\overset{G}{\mathrm{\Σ}}}{p}_{m,g,t}=p_{m,t}$

(5) generating flow constraint

${\underset{\‾}{q}}_{m,t}\≤q_{m,t}\≤{\stackrel{\‾}{q}}_{m,t}$

In formula: q _m,t, q _m,trepresent the generating flow at period t and the bound thereof of power station m respectively.

(6) storage outflow constraint

${\underset{\‾}{S}}_{m,t}\≤S_{m,t}\≤{\stackrel{\‾}{S}}_{m,t}$

In formula: S _m,t, s _m,trepresent that power station m is in the storage outflow of period t and bound thereof respectively.

(7) output of power station constraint

${\underset{\‾}{p}}_{m,t}\≤p_{m,t}\≤{\stackrel{\‾}{p}}_{m,t}$

In formula: p _m,trepresent the average output bound of power station m at period t respectively.

(8) reservoir level constraint

${\underset{\‾}{Z}}_{m,t}\≤Z_{m,t}\≤{\stackrel{\‾}{Z}}_{m,t}$

In formula: Z _m,t, z _m,trepresent that power station m is in the upper pond level of period t and bound thereof respectively.

(9) output of power station climbing restriction

$\left\{\begin{array}{ccc}(1+{\mathrm{\μ}}_m){\stackrel{\‾}{p}}_{m,t-1}-p_{m,t}\≥0& \mathrm{if}& p_{m,t}>p_{m,t-1}\\ (1-{\mathrm{\μ}}_m){\stackrel{\‾}{p}}_{m,t-1}-p_{m,t}\≤0& \mathrm{else}& \end{array}\right.$

In formula: for power station m is in the average output upper limit of period t-1; μ _mfor the climbing capacity of power station m.

(10) power station startup-shutdown sustained periods of time requirement

$\left\{\begin{array}{ccc}{p}_{m,t}>0& \mathrm{if}& p_{m,t-t{g}_m}=0\mathrm{and}{p}_{m,t-1}>0\\ p_{m,t}=0& \mathrm{if}& p_{m,t-t{s}_m}>0\mathrm{and}{p}_{m,t-1}=0\\ p_{m,t}\≥0& \mathrm{else}& \end{array}\right.$

In formula: for power station m is at period t-tg _maverage output; for power station m is at period t-ts _maverage output; t _gm, t _smrepresent the minimum start of power station m respectively and shut down sustained periods of time.

(11) vibrating area constraint

$(p_{m,t}-{\stackrel{\‾}{\mathrm{ps}}}_{m,t,k})(p_{m,t}-{\underset{\‾}{\mathrm{ps}}}_{m,t,k})>0$

In formula: ps _{m, t, k}represent the bound of power station m at a kth vibrating area of period t respectively.

(12) minimum start units limits

(p _m,t-pmin _m)p _m,t≥0

In formula: pmin _mfor the minimum service capacity of power station m.

According to above-mentioned thought, once complete Optimized Operation process, is achieved according to following step (1)-(8):

(1) initial calculation condition, comprises hydroelectric station operation condition and constraint, and each network load process and by electrical demand;

(2) power station numbering m=1 is made;

(3) what build m power station faces workload demand { C ₁, C ₂..., C _t, wherein C ₁, C ₂..., C _tbe respectively the period 1,2 ..., the load value of T;

(4) with { C ₁, C ₂..., C _tbe workload demand, under given objective function and constraint condition, adopt DFP-method to determine day 96 processes of exerting oneself in m power station;

(5) under given controlled condition, outer point searching method is adopted to determine the output distribution process of m power station at G electrical network, i.e. { p _{m, 1,1}, p _{m, 1,2}..., p _{m, 1, T}, { p _{m, 2,1}, p _{m, 2,2}..., p _{m, 2, T}..., { p _{m, G, 1}, p _{m, G, 2}..., p _{m, G, T}, wherein, { p _{m, 1,1}, p _{m, 1,2}..., p _{m, 1, T}, { p _{m, 2,1}, p _{m, 2,2}..., p _{m, 2, T}..., { p _{m, G, 1}, p _{m, G, 2}..., p _{m, G, T}represent that m power station is the 1st, 2 respectively ..., the 1st, 2 of G electrical network ..., the load value of T period;

(6) each electrical network residue load is recalculated;

(7) make m=m+1, if m is less than or equal to M, then skip to step (3); Otherwise go to step (8);

(8) calculate the desired value of current solution, judge whether target is improved, if so, using the initial solution of current solution as next round iteration, and skip to step (2); Otherwise, judge whether to reach to maximum iteration time, if so, algorithm search stops, and exports optimal result, otherwise, skip to step (2).

Now using Southwestern China area Korean pine bark rot Tianshengqiao-I (being called for short sky one), natural bridge secondary (being called for short sky two), Ping Ban and all the people present power station, imperial beach as research object, the inventive method is adopted to make its operational plan a few days ago.It one, sky two, power station, three, imperial beach be the controlling power station of whole basin step, directly adjusts pipe by CHINA SOUTHERN POWER, to Guangdong, Guangxi two province provides electric power, bears the peak regulation task of two provincial power networks.Power station basic data is as shown in table 1, and result of calculation as shown in Figure 4.From map analysis, when meeting given power transmission amount, the inventive method, by distributing electric power between coordinating net, makes Guangdong and Guangxi Power Grid is powered all concentrates on respective peak period, effectively meet the differentiation peak-load requirement of two electrical networks, achieve good peak regulation effect.Due to about 6.3 times that Guangdong Power Grid payload is Guangxi Power Grid, thus when identical by electricity, Guangxi Power Grid peak regulation range is obviously greater than Guangdong Power Grid, and the more former load of the former remaining load peak-valley difference reduces 44%, and the latter reduces 6%., contrast with the actual running results, the peak regulation range of two electrical networks on average improves about 5%, and wherein Guangxi Power Grid peak regulation range is up to 12% meanwhile, and electrical network residue load is more level and smooth, and peak regulation successful improves.

Table 1

Claims (1)

1. the two benches of a power transmission Hydropower Stations peak regulation scheduling transprovincially searching method, its feature comprises the steps:

(1) initial calculation condition, comprises hydroelectric station operation condition and constraint, and each network load process and by electrical demand;

(2) power station numbering m=1 is made;

(3) what build m power station faces workload demand { C ₁, C ₂..., C _t, wherein C ₁, C ₂..., C _tbe respectively the period 1,2, the load value of T;

(4) with { C ₁, C ₂..., C _tbe workload demand, under given objective function and constraint condition, adopt DFP-method to determine day 96 processes of exerting oneself in m power station;

(5) under given controlled condition, outer point searching method is adopted to determine the output distribution process of m power station at G electrical network, i.e. { p _{m, 1,1}, p _{m, 1,2}..., p _{m, 1, T}, { p _{m, 2,1}, p _{m, 2,2}..., p _{m, 2, T}..., { p _{m, G, 1}, p _{m, G, 2}..., p _{m, G, T}, wherein, { p _{m, 1,1}, p _{m, 1,2}..., p _{m, 1, T}, { p _{m, 2,1}, p _{m, 2,2}..., p _{m, 2, T}..., { p _{m, G, 1}, p _{m, G, 2}..., p _{m, G, T}represent that m power station is the 1st, 2 respectively, the 1st, 2 of G electrical network, the load value of T period;

Wherein, outer point searching method solution procedure comprises following two major parts:

A () determines initial solution; Stationary electric power plant exert oneself process condition under, adopt following formula to obtain each electrical network by electricity, so with by electricity for control objectives, face by power station based on load, adopt successively cutting load method to determine the powered process of exerting oneself of the ideal of each electrical network;

E _m，g＝E _m×R _m，g

In formula: E _m,gfor m power station is to the day power transmission amount of g electrical network; E _mrepresent the total electricity of power station m in schedule periods; R _m,gfor the electricity ratio of g electrical network is sent in m power station,

Based on b initial solution that () obtains by step (a), reduce the low remaining lotus period and exert oneself according to preferential and preferentially increase the rule that the high remaining lotus period exerts oneself, redistribute the day part electrical network not meeting period power balance constraint poweredly to exert oneself, until all periods all reach power balance;

(6) each electrical network residue load is recalculated;

(7) make m=m+1, if m is less than or equal to M, then skip to step (3); Otherwise go to step (8); Wherein, M is power station sum;

(8) calculate the desired value of current solution, judge whether target is improved, if so, using the initial solution of current solution as next round iteration, and skip to step (2); Otherwise, judge whether to reach to maximum iteration time, if so, algorithm search stops, and exports optimal result, otherwise, skip to step (2).